Comparative and functional genomics provide insights into the pathogenicity of dermatophytic fungi
- Anke Burmester†1, 2,
- Ekaterina Shelest†3,
- Gernot Glöckner†4,
- Christoph Heddergott†1, 2,
- Susann Schindler5, 6,
- Peter Staib7,
- Andrew Heidel4,
- Marius Felder4, 8,
- Andreas Petzold4,
- Karol Szafranski4,
- Marc Feuermann9,
- Ivo Pedruzzi9,
- Steffen Priebe3,
- Marco Groth4,
- Robert Winkler6, 10,
- Wenjun Li11,
- Olaf Kniemeyer1,
- Volker Schroeckh1,
- Christian Hertweck6, 10,
- Bernhard Hube6, 12,
- Theodore C White13,
- Matthias Platzer4,
- Reinhard Guthke3,
- Joseph Heitman11,
- Johannes Wöstemeyer2,
- Peter F Zipfel5, 6,
- Michel Monod14 and
- Axel A Brakhage1, 2Email author
© Burmester et al.; licensee BioMed Central Ltd. 2011
Received: 20 July 2010
Accepted: 19 January 2011
Published: 19 January 2011
Millions of humans and animals suffer from superficial infections caused by a group of highly specialized filamentous fungi, the dermatophytes, which exclusively infect keratinized host structures. To provide broad insights into the molecular basis of the pathogenicity-associated traits, we report the first genome sequences of two closely phylogenetically related dermatophytes, Arthroderma benhamiae and Trichophyton verrucosum, both of which induce highly inflammatory infections in humans.
97% of the 22.5 megabase genome sequences of A. benhamiae and T. verrucosum are unambiguously alignable and collinear. To unravel dermatophyte-specific virulence-associated traits, we compared sets of potentially pathogenicity-associated proteins, such as secreted proteases and enzymes involved in secondary metabolite production, with those of closely related onygenales (Coccidioides species) and the mould Aspergillus fumigatus. The comparisons revealed expansion of several gene families in dermatophytes and disclosed the peculiarities of the dermatophyte secondary metabolite gene sets. Secretion of proteases and other hydrolytic enzymes by A. benhamiae was proven experimentally by a global secretome analysis during keratin degradation. Molecular insights into the interaction of A. benhamiae with human keratinocytes were obtained for the first time by global transcriptome profiling. Given that A. benhamiae is able to undergo mating, a detailed comparison of the genomes further unraveled the genetic basis of sexual reproduction in this species.
Our results enlighten the genetic basis of fundamental and putatively virulence-related traits of dermatophytes, advancing future research on these medically important pathogens.
The two dermatophyte species Arthroderma benhamiae and Trichophyton verrucosum are both zoophilic, yet the natural reservoir of T. verrucosum is almost exclusively cattle, whereas A. benhamiae is usually found on rodents, in particular guinea pigs [12, 13]. The two species also differ in their ability to grow under laboratory conditions, with T. verrucosum being very difficult to cultivate at all . Conversely, A. benhamiae is comparatively fast growing and produces abundant microconidia. As a teleomorphic species, the fungus is even able to undergo sexual development, including the formation of sexual fructifications (cleistothecia) [15, 16]. These characteristics, together with the recent establishment of a guinea pig infection model and a genetic system for targeted gene deletion (P Staib and colleagues, manuscript submitted) for this species, suggest A. benhamiae is a useful model organism to investigate the fundamental biology and pathogenicity of dermatophytes . Despite the above mentioned phenotypic differences, A. benhamiae and T. verrucosum are phylogenetically very closely related, and both induce highly inflammatory cutaneous infections in humans, such as tinea corporis [15, 17]. Therefore, a genome comparison of the two species should reveal common basic pathogenicity-associated traits.
In the present study, we report and compare the genome sequences of A. benhamiae and T. verrucosum and refer to potential dermatophyte-specific pathogenicity-associated factors, as revealed by comparisons with groups of proteins important for pathogenicity in other species of the Onygenales (Coccidioides posadasii and Coccidioides immitis) and in the mould A. fumigatus. Applying our insights thereof, we used secretome analysis to reveal secreted factors of A. benhamiae that mediate extracellular in vitro keratin degradation. The interaction between A. benhamiae and the human host was monitored by global transcriptome profiling of the fungal cells in contact with human keratinocytes. Investigating the molecular basis of sexual reproduction, we inspected in detail the A. benhamiae mating type locus.
Results and discussion
Comparative genomics of A. benhamiae and T. verrucosum
Genome data of A. benhamiae and T. verrucosum
Mean CDS length
Genes with introns
However, we also identified notable dissimilarities between the genomes of A. benhamiae and T. verrucosum. After having detected the orthologous pairs with best bidirectional hits, we came up with lists of proteins that presumably were unique for either species. Since the best bidirectional hits were identified using protein Blast, we next applied BlastN to correct for possible gene prediction errors. We used a filter threshold for significant hits of 80% identity between sequences over less than 50% of the query length. There were 238 A. benhamiae sequences that gave no hits or non-significant hits in T. verrucosum, and 219 T. verrucosum genes were not found in A. benhamiae. Of these, 83 and 78 genes (A. benhamiae and T. verrucosum, respectively) have assigned names and/or functional domains. A list of the predictions is provided in Additional file 4. Given the overall strong genome sequence similarity, a future functional investigation of these distinctions appears to be of interest, in particular with respect to the tremendous differences between the two species in terms of in vitro growth ability and animal host preference (see also the 'Other interesting genes' section).
We analyzed the A. benhamiae fast-evolving genes in comparison to T. verrucosum. Using the dN/dS ratio as a measure for selective pressure, we obtained a list of positively selected genes (dN/dS >1) (Additional file 5). In total we found 132 positively selected genes with assigned functions, enabling assumptions about their roles in the cell and, hence, the reasons for their accelerated evolution. Of particular interest are the two most abundant groups of these genes, those encoding transcription factors (18 genes) and MFS transporters (5 genes). The latter are known to be usual constituents of secondary metabolite (SM) gene clusters.
Both dermatophyte genomes encode the basic metabolic machinery for glycolysis, tricarboxylic acid cycle, glyoxylate cycle, pentose phosphate shunt, and synthesis of all 20 standard amino acids and the five nucleic acid bases. Moreover, dermatophytes appear to be capable of producing a wide range of SMs, which is reflected by the presence of polyketide synthase (PKS)- and non-ribosomal peptide synthetase (NRPS)-encoding genes (see the 'Genetic basis for secondary metabolism gene clusters' section). The outstanding ability of dermatophytes to specifically infect superficial host structures may be supported by the possession of a broad repertoire of genes encoding hydrolytic enzymes, the expression of many of which was also proven experimentally (see the next paragraph and the 'Identification of secreted fungal proteins during keratin degradation by secretome analysis' section). In addition, the ability of dermatophytes to assimilate lipids, major constituents of the skin, is putatively reflected by the presence of 16 lipase genes in either genome. A putative link between the possession of lipases and fungus-induced skin disease has previously been revealed for basidiomycetes of the genus Malassezia .
Of particular note is the apparent relative paucity of tRNA genes in both dermatophytes in comparison with other closely related ascomycetes. The genomes of A. benhamiae and T. verrucosum contain 80 and 77 tRNA genes, respectively, whereas the number of tRNA genes varies between approximately 100 to 130 in Coccidioides species and 150 to 370 in aspergilli. However, some strains of H. capsulatum, representing a comparatively closely related pathogen, also possess only 83 to 89 tRNA genes, suggesting that the low number of tRNA genes is not specific to dermatophytes.
Identification of a broad repertoire of protease genes in dermatophyte genomes
Dermatophytes are keratinophilic fungi, sharing the ability to utilize compact hard keratin as a sole source of carbon and nitrogen. In line with this knowledge, the two sequenced genomes reflect a remarkable metabolic capability for protein degradation. They contain 235 predicted protease-encoding genes, 87 of the deduced proteins possessing a secretion signal (Table S3 in Additional file 6). We did not detect any protease in A. benhamiae or T. verrucosum unique to either species, a finding that may reflect similar life styles and/or host adaptation mechanisms, especially with respect to their common keratinophilic growth. In general, deviations in the number of proteases per genome are rather large in the fungal kingdom, ranging from approximately 90 in Ustilago maydis to approximately 350 in Gibberella zeae (according to the MEROPS database ). Dermatophytes belong to the most protease-rich species.
The protein sequence of each protease is highly conserved across dermatophyte species . Collections of predicted secreted proteases of A. benhamiae and T. verrucosum as well as Coccidioides spp. (Onygenales) were compared to those of A. fumigatus as a member of the Eurotiales, for which many secreted proteases have previously been characterized. Most A. fumigatus proteases in A1 (pepsins), M28 (leucine aminopeptidases), S9 (dipeptidylpeptidases), S10 (carboxypeptidases) and S53 (tripeptidylpeptidases) families have an orthologue in dermatophytes and Coccidioides spp. (Table S4 in Additional file 7). The major striking differences found between the secreted protease batteries of A. fumigatus and Onygenales are the following: subtilisin (S8), deuterolysin (M35), and fungalysin (M36), which belong to endoprotease gene families, have expanded in Onygenales (Table S4 in Additional file 7); the same is true for exopeptidases of the M14 family (metallocarboxypeptidases) and the M28 family (aminopeptidases) - a major carboxypeptidase (McpA) homologous to the human pancreatic carboxypeptidase A was previously characterized in dermatophytes , and of particular note, Aspergillus spp. have no McpA orthologue; and genes encoding acidic glutamic proteases (G1 family) were not detected in either dermatophytes or Coccidioides spp.
Major differences between dermatophytes and Coccidioides spp. proteases were found in M35, M36 and S8 proteases families (see the phylogenetic trees in Additional file 8). Proteases of these three families of dermatophytes and Coccidioides spp. form distinct clades in phylogenetic trees (Additional file 8). Members of the S8 and M36 families have undergone additional amplifications in the dermatophyte lineage, and expansion of the M35 family appears to be different in Coccidioides spp. and dermatophytes. In the latter, a clade was apparently lost. In addition, three genes encoding proteases of the S41 family were found in the dermatophyte genomes while only one gene encoding a protease of this family was identified in Coccidioides spp.
Recent comparative genomic analyses of Coccidioides species with other members of the Onygenales showed gene family sizes are associated with a host/substrate shift from plants to animals in these microorganisms . Experimentally, the expression of genes encoding fungalysins and subtilisins was recently monitored in A. benhamiae by cDNA microarray analysis during growth on keratin, and also during cutaneous infection of guinea pigs . Interestingly, the prominent keratin induced A. benhamiae subtilisin-encoding genes, such as SUB3 and SUB4, were not observed in this former analysis to be strongly activated in vivo, in contrast to others that conversely were not found to be induced during in vitro growth on keratin. A role for Sub3 was recently observed in adhesion of the dermatophyte Microsporum canis to feline epidermis, but not for the invasion thereof . These findings suggest additional functions of secreted proteases during host adaptation other than keratin degradation. Since the formerly used cDNA microarray does not comprise the full genome of A. benhamiae, the future identification of in vivo specific dermatophyte proteases on the basis of the presented genome appears to be of major interest.
Identification of secreted fungal proteins during keratin degradation by secretome analysis
Endo- and exoproteases secreted by microorganisms cooperate very efficiently in protein digestion to produce oligopeptides and free amino acids that can be incoporated via transporters. During the process of protein digestion the main function of endoproteases is to produce a large number of free end peptides on which exoproteases may act. At neutral and alkaline pH, synergistic action of Lap and DppIV was shown in Aspergillus spp. [13, 24]. Laps degrade peptides from the amino terminus until reaching an X-Pro sequence, which acts as a stop. In a complementary manner, the X-Pro sequences can be removed by DppIV, thus allowing Laps access to the next residue. Dermatophyte and Aspergillus spp. Lap1, Lap2, DppIV and DppV have shown comparable substrate specificity . Therefore, our proteomics approach allows us to hypothesize common basic mechanisms in dermatophytes during extracellular protein digestion. However, the presence of large protease gene families in dermatophytes reflects selection during evolution and the ability of these fungi to adapt to different environmental conditions during infection and saprophytic growth.
Differential gene expression in A. benhamiaeduring infection of keratinocytes
Growth of A. benhamiae on keratin might mimic selected in vivo growth substrates, yet may not reflect the entire process of infection. In order to gain more insights into basic host adaptation mechanisms, we studied the global transcriptional response of A. benhamiae during infection of human keratinocytes. After 12 h of co-cultivation, germinating A. benhamiae microconidia were observed to be localized and concentrated on the host cells, suggesting that the fungus actively adheres to the keratinocytes (Figure 1c,d). To perform 454 RNA sequencing, the fungal cells were harvested after incubation for 96 h with and without keratinocytes. About 50 A. benhamiae genes showed differential expression with a fold change >5 (P-value < 0.05; Table S7a in Additional file 6); 45 genes encoding putatively secreted proteins (Table S5 in Additional file 6) and 13 genes coding for proteins involved in the biosynthesis of SMs are expressed either only with or without keratinocytes, or under both conditions. Of the 235 predicted protease-encoding genes, 158 are expressed under both conditions. Sixteen potentially secreted proteins, including three proteases, are differentially expressed (Table S7b in Additional file 6). In particular, the expression profile of the genes encoding carboxypeptidase S1 and dipeptidyl-peptidase DppV implies their potential involvement in the infection process. The transcript levels of two NRPS genes were reduced during co-cultivation with keratinocytes, a finding that is noticeable but cannot be explained at this stage.
To confirm the RNAseq results, we selected several genes that were predicted to be differentially expressed and tested them by Northern blotting. We used two housekeeping genes, actin (ARB_04092) and glyceraldehyde 3-phosphate dehydrogenase (GAPDH, ARB_00831), as controls as they are not expected to be differentially regulated between the control and co-incubation conditions. All tested genes were regulated as expected from the RNAseq data (Figure S4 in Additional file 9). The expression level alterations of metabolic enzymes (ARB_07891, ARB_04156, ARB_01650 and ARB_04856) and membrane transporters (ARB_01027) reflect the adaptation of the fungus to the different nutrition provided by keratinocytes and their remnants, whereas the strong up-regulation of the hydrophobin ARB_06975 indicates altered binding properties and adhesivity during growth on epithelial cells and during infection. In conclusion, this independent experimental method shows that the accuracy of the RNAseq data was exemplary.
Genetic basis for secondary metabolism gene clusters
Putative PKS and NRPS genes of A. benhamiae, T. verrucosum, and C. immitis
LocusLink Arthroderma benhamiae
LocusLink Trichophyton verrucosum
LocusLink Coccidioides immitis
All nine PKS genes detected in A. benhamiae have unequivocal counterparts in the T. verrucosum genome (Table 2). An interesting feature of the dermatophyte PKS set is the unusual proportions of reducing and non-reducing PKSs. Whereas in all other closely related ascomycetes (such as aspergilli) most of the PKSs are non-reducing, in dermatophytes most are reducing PKSs. A comparison with the closest sequenced relative, C. immitis (Table 2; see more details below), also revealed substantial differences in the composition of the PKS set: the ratio of reducing to non-reducing in dermatophytes is 2:1, whereas in C. immitis it is 2:3. This observation suggests dermatophytes have an uncommon SM profile, which deserves future investigation. Particular attention should be paid to the fact that these fungi are characterized by intense pigmentation, a phenotype that may be related to their pathogenicity. For the related species T. rubrum, the polyketide-derived mycotoxin xanthomegnin has been suggested to be responsible for the characteristic red colony reverse pigment. Most interestingly, xanthomegnin production has even been detected in epidermal material infected by T. rubrum, in contrast to non-infected controls . A putative link between SM production and host adaptation of A. benhamiae might also be reflected by our observation that several genes associated with the synthesis of such molecules were found to be differentially regulated during infection of human keratinocytes (see the 'Differential gene expression in A. benhamiae during infection of keratinocytes' section).
To get an impression of possible expansions of families and evolutionary relationships, we compared the sets of SM producers in dermatophytes with that of C. immitis (Table 2; Figure S5.1 and S5.2 in Additional file 10). As mentioned above, the total number of SM gene clusters is higher in dermatophytes, mainly due to the more abundant NRPSs. However, we observe differences also in the PKS set as well as in the number of PKS/NRPS hybrids: C. immitis possesses only one hybrid, whereas each dermatophyte has three. The higher number of non-reducing PKSs in C. immitis is mainly due to the expansion of one clade; most likely we are seeing the results of duplication of some ancestor genes with a domain architecture of a beta-ketoacyl synthase domain, an acetyltransferase domain, an acyl carrier protein domain, and a methyltransferase domain (KS-AT-ACP-ME). Four of six C. immitis non-reducing PKSs belong to this clade. Of the other two, one has a clear ortholog in dermatophytes, and the other has an unusual structure (AT-KS-ACP-thioesterase domain (TE)) without an orthologous dermatophyte gene. In comparison to C. immitis, dermatophytes possess two additional non-reducing clades, which means that, in spite of the lower number of non-reducing PKSs, they have more various potential capacities. The reducing C. immitis PKSs also cannot boast great variety: two of four C. immitis genes are most likely the result of a duplication (they form a separate clade and do not have dermatophyte orthologs), one PKS has orthologs in dermatophytes, and one is only a probable homolog (see below). On the other hand, in dermatophytes we see an expansion of the group with a fumonisin synthase-like structure (KS-AT-ME-enoyl reductase domain (ER)-ketoacyl reductase domain (KR)-ACP): three orthologous pairs formed by out-paralogs in each species have only one close homolog in C. immitis. Since the C. immitis gene lacks one of the domains (methyltransferase), we cannot consider it as a fumonisin-like ortholog. Besides the 6-methyl-salicylic acid synthase, completely lacking in C. immitis, another not completely reducing PKS (KS-AT-ME-KR-ACP), as well as two PKS/NRPS hybrids, do not have homologs in C. immitis. Taken together, these data agree with our hypothesis that highly adapted parasites such as Coccidioides do not require a large arsenal of SMs.
Sexuality in dermatophytes
We did not identify a striking defect in the T. verrucosum mating type locus, which appears to be intact. Several strains of T. verrucosum were found to be of the same mating type as the sequenced strains, suggesting a strong disequilibrium towards mating type +.
In Aspergillus (Eurotiales), Coccidioides and Histoplasma (Onygenales) the mating type (MAT) loci are flanked by APN2 and the SLA2 genes encoding a DNA lyase and a cytoskeleton protein, respectively . The MAT idiomorphs and flanking regions described here for A. benhamiae and T. verrucosum are essentially identical to those of other closely related dermatophytes .
Other interesting genes
Of particular interest are the genes of A. benhamiae that have no obvious counterpart in T. verrucosum (Additional file 4) and whose predicted functions suggest their potential involvement in basic biological phenotypes and/or pathogenicity. Two such genes, ARB_04713 and ARB_02149, encoding a phosphopantetheine-binding domain and an NRPS, respectively, were found in the transcriptome analysis, although not expressed differentially. The expression pattern of the A. benhamiae-specific NRPS ARB_02149 further suggests that its as yet unidentified product is produced during infection by the fungal cells.
Another gene of particular interest encodes hydrophobin. In A. fumigatus, surface hydrophobin was shown to prevent immune recognition . The A. benhamiae hydrophobin gene (ARB_06975) shows 99% similarity with the respective T. verrucosum gene (TRV_00350) and displays moderate overexpression (1.6×) under co-cultivation conditions (Table S7b in Additional file 6). The analysis of a potential role of dermatophyte hydrophobins in immune response functions and/or adhesion to host surfaces will be part of future research.
Numerous examples in microbial pathogenicity research still need to be explained at the genomic level, thus requiring genome sequences to be made available. Here, we present the first genomes of dermatophyte species, filamentous fungi that cause most superficial infections in humans and animals. The presence of putative pathogenicity-related factors, such as numerous secreted proteases, was revealed at the genome level and also experimentally confirmed during keratin degradation by A. benhamiae. Although keratin utilization is traditionally supposed to be of major relevance for the pathogenicity of these microorganisms, the entire process of host adaptation during infection seems to be more complex. Transcriptome analysis showed that only some of the typically keratin-induced proteases were found to be strongly expressed during fungus-keratinocyte interaction. Instead, genome and transcriptome analyses draw attention to so far hardly noticed dermatophyte factors - for example, putative SMs - the role of which should be addressed in the future. Our research on dermatophytes was strongly facilitated by the selection of A. benhamiae as a model species, which provides practical advantages such as comparatively fast growth and the production of abundant microconidia. Moreover, future basic studies on the regulation of mating, dermatophyte evolution and host preference will profit from the ability of A. benhamiae to undergo sexual reproduction. In conclusion, by presenting dermatophyte genomes and global insights into major processes of host adaptation, we intend to advance molecular studies on these medically important microorganisms.
Materials and methods
A. benhamiae and T. verrucosumstrains and growth conditions
A clinical isolate of A. benhamiae strain 2354 was used (isolate LAU2354) . T. verrucosum strain 44  was kindly provided by Yvonne Gräser (Charité, Berlin, Germany). Strains were cultivated at 28°C on Sabouraud 2% (w/v) glucose agar (SG, Merck, Darmstadt, Germany) for 12 days; liquid cultures were shaken at 180 rpm at 30°C for 5 to 7 days. Hyphae and conidia were separated by filtration using a 40 μm cell strainer (BD Bioscience, Heidelberg, Germany). Conidia were counted with a cell counter (Beckman, Coulter, Krefeld, Germany) or manually using a Thoma chamber. For crossing experiments of A. benhamiae LAU2354 with the opposite mating type CBS 809.72, MAT medium  (1/10 SG, 0.1% (w/v) MgSO4 and 0.1% (w/v) KH2PO4) was used.
DNA and RNA preparation for DNA sequencing and cDNA library
For DNA preparation, mycelia were separated from the medium by filtration through Miracloth (Calbiochem, Darmstadt, Germany) and ground in a mortar under liquid nitrogen. After evaporation, the powder was suspended in a solution containing 150 mM EDTA, 50 mM Tris-HCl, pH 8.0, 1% (w/v) SDS, 20 mM NaCl and 100 μg/ml proteinase K (Merck). After incubation for 1 h at 55°C, the solution was gently mixed with 1/4 volume of 4 M NaCl and kept on ice for 30 minutes. After centrifugation for 10 minutes at 6,000 rpm and 4°C, polyethylene glycol 6000 (Serva, Heidelberg, Germany) was added to the supernatant to a final concentration of 10% (w/v). The DNA was precipitated for 1 h on ice and centrifuged for 10 minutes at 10,000 rpm at 4°C. The pellet was dissolved in a solution containing 25 mM Tris-HCl, pH 8.0, 5 mM EDTA, 10 mM NaCl and 1% (v/v) Triton X100. For density centrifugation, 1g CsCl and 12 μl bisbenzimide (10 mg/ml) for each milliliter of solution were added [41, 42]. Ultracentrifugation was performed in a vertical rotor at 44,000 rpm for 24 h at 25°C. DNA was separated into two bands of different density according to the AT-content of the DNA. The upper band contained a DNA fraction highly enriched for mitochondrial DNA. For T. verrucosum, two rounds of density gradient centrifugation were necessary. In the first round, ethidium bromide was used instead of bisbenzimide. For RNA preparation, SG medium was inoculated with conidia to a final concentration of 3 × 104 conidia/ml and shaken at 180 rpm for three days at 30°C. Total RNA was isolated using a commercial kit as described by the manufacturer (Qiagen, Hilden, Germany). After RNA extraction, a cDNA library was constructed from this material according to the manufacturer's protocols (MINT cDNA synthesis kit, Evrogen, Moscow, Russia).
Plasmid/fosmid libraries and sequencing
Nuclear DNA of A. benhamiae and T. verrucosum was sheared, size fractionated (3 to 4 kb), end-repaired, and cloned into the SmaI site of pUC18. For both fungal species, two fosmid libraries each were prepared in pCC1FOS (Epicentre Biotechnologies, Madison, WI, USA) as described by the manufacturer, one for the high-GC chromosomal DNA fraction and one for the AT-rich mitochondrial DNA fraction. For T. verrucosum, 40,000 fosmids from GC-rich and 80,000 fosmids from AT-rich DNA were obtained. For A. benhamiae, the corresponding yields were around 50,000 (GC-rich) and 20,000 (AT-rich), respectively. End sequences of plasmid and fosmid clones were obtained using dye terminator chemistry and a 3730×l sequencer (Applied Biosystems, Foster City, CA, USA). Moreover, a fosmid library was generated with a GC-rich DNA fraction of the A. benhamiae strain CBS 809.72 encoding the opposite mating type locus. We tested 1,000 fosmids by colony filter hybridization and in PCR experiments. Fosmids were identified by hybridization experiments with a digoxygenin-labeled part of the apn2 gene (ARB_07318) and in PCR experiments using apn2 amplifying primers (5'-CTTCTAGTGACTCGCCACAGG-3' forward and 5'-GAGTTGGAGGTTGAGATGCTGAC-3' reverse). Three clones were positively identified by both methods. To test whether the fosmids contained the full length MAT region, the clones were tested in PCR experiments amplifying parts of other flanking genes, such as the sla2 gene (ARB_07317) and the rps4 gene (ARB_07322). For sla2, PCR primer pair 5'-CTTGTTCAGGAGAGCTATGG-3' and 5'-CAGCTTCTCGAGCTCCTCCC-3' was used; for rps4, PCR primer pair 5'-CAGCGCCTGGTCAAGGTCGACG-3' and 5'-GGTCACGCTCCTCAGCAATGG-3' was used. DNA of a positive fosmid was shotgun sequenced using dye terminator chemistry (ABI).
In addition, genomic 454 libraries were generated according to the manufacturer's protocol and sequenced using a GS FLX (Roche, Mannheim, Germany). The nucleotide sequences were assembled species-specific using the newbler software. Clone gaps were filled using a primer walking strategy with custom primers. Isolation, quantification and quality control of total RNA was performed as described . A cDNA-library was constructed according to the manufacturer's protocols (Evrogen) and 1,411 ABI dye terminator sequences were obtained mainly from the 5' end. The sequences were matched to the assembled genomic sequences to determine exon/intron structures and to obtain an intron signature for the species.
Next generation sequencing and assembly
The same DNA as for the preparation of the plasmid/fosmid libraries was used for the preparation of genomic libraries for the 454/FLX system (Roche) according to the manufacturer's protocols. Three runs each were performed on the 454/FLX sequencing machine. All 454/FLX sequence data were assembled species-specific using the newbler software. The Sanger based sequencing reads were assembled onto this 'backbone'. Clone gaps were filled using a primer walking strategy with custom primers.
Both genomes are deposited in NCBI with accession codes [DDBJ/EMBL/GenBank:ABSU00000000] for A. benhamiae and [DDBJ/EMBL/GenBank:ACYE00000000] for T. verrucosum.
Gene models of both fungi were generated by using in silico predictions and sequence data from an EST library constructed from cultured A. benhamiae cells. We matched 1,411 ABI dye terminator sequences obtained from the cDNA library sequencing to the assembled genomic sequences to determine exon/intron structures and to obtain an intron signature for the species. The alignments of the cDNA sequences to the genomic backbone yielded evidence for 861 introns and at least 653 protein-coding open reading frames (coding sequences), which were validated by BLAST. These data were used to train the gene prediction program geneid . To validate the accuracy of the gene prediction, 47 gene structures in one genomic region were annotated manually and compared to the automated predictions, indicating a specificity of 82% at a sensitivity of 97%. For the annotation and comparative analyses of the genomes a web based genome browser was set up using the GenColors database/software system .
Best bidirectional hits and BlastN analysis
Blast analysis of all coding sequences of one genome against the other yielded best bidirectional hits. We used a filter threshold for significant hits of 30% identity between amino acid sequences over at least 50% of the protein.
A BlastN analysis of the genomic sequences was performed for all protein coding genes of T. verrucosum against all A. benhamiae contigs. A filter threshold for significant hits was 80% identity between sequences over at least 60% of the query length; 239 T. verrucosum sequences gave no hits or non-significant hits.
The human keratinocyte line HaCaT was obtained from Prof. Fusenig (Deutsches Krebsforschungszentrum, Heidelberg, Germany). The cells were cultivated in DMEM supplemented with 10% (v/v) fetal calf serum, gentamycin (28 μg/ml) and 1% (w/v) ultraglutamine at 37°C in a humidified atmosphere and 5% (v/v) CO2 for 2 days. Medium and supplements were purchased from Lonza (Basel, Belgium). Human keratinocytes were infected by A. benhamiae conidia with a multiplicity of infection (MOI) of 6. Infected human cells were cultivated in fetal calf serum-free DMEM supplemented with both gentamycin and ultraglutamine for 96 h at 28°C. As a control, A. benhamiae conidia were grown in the absence of keratinocytes under the same conditions. After infection, the human keratinocytes were lysed by addition of 0.03% (v/v) Triton X for 2 minutes and A. benhamiae was harvested. Fungal cells were collected by centrifugation for 3 minutes at 3,500 g. A. benhamiae cells were washed twice in Dulbecco's phosphate buffered saline (Lonza) and stored in aliquots at -80°C. For RNA sequencing, total RNA was isolated using RiboPure™-Yeast Kit (Ambion Europe, Huntingdon, UK) according to the manufacturer's instructions from keratinocytes co-incubated with conidia and conidia only grown in cell culture medium for 96 h.
RNA was reverse transcribed using a SMART technique (Evrogen). The single-stranded DNA was then amplified using SMART primers for 20 cycles to produce double-stranded DNA in sufficient quantity for GS-FLX sequencing (Roche). We generated A. benhamiae transcriptome data by sequencing parts of individual cDNAs after fragmentation by nebulization using 454/FLX sequencing technology. For postprocessing of these sequences, SMART adapters were identified and clipped using a combination of perl scripts plus cross_match. After further cleaning with seqclean (removal of polyA tails and low complexity regions), 682,580 ESTs (98.8 Mb) remained for mapping.
Mapping of the ESTs to the repeat-masked A. benhamiae genome as a backbone was done in two major steps. First, we used Blat  to assign each EST to its most probable position in the genome allowing a maximum intron length of 10 kb. A valid hit required a minimum length of 30 bp and a minimum identity of 90% to the backbone sequence. In the second step, each EST was realigned to its most probable position utilizing a slightly modified version of Exalin  that implements the Smith-Waterman algorithm and information theory for better alignments and intron prediction. Using this approach, we were able to align 571,963 ESTs to the genome. Finally, EST positions were translated to positions of known gene models if possible. In this way, we determined for each gene a set of ESTs and thereby its raw expression level. The data were normalized to the total number of mapping ESTs. Table S9 in Additional file 6 shows the total numbers of generated reads, the reads mapped to a genome, and the reads in gene models for each technical replicate of infection and control samples.
The raw counts for the transcripts were analyzed using the R Statistical Computing Environment and the Bioconductor packages DESeq  and edgeR . Both packages provide statistical routines for determining differential expression in digital gene expression data using a model based on the negative binomial distribution. The resulting P-values were adjusted using the Benjamini and Hochberg's approach for controlling the false discovery rate . Genes with an adjusted P-value <0.05 found by both packages were assigned as differentially expressed.
The RNAseq data are submitted to the Sequence read archive of NCBI and are available with the accession numbers [NCBI:SRR070551] and [NCBI:SRR070552] (sample runs) and [NCBI:SRR070553] and [NCBI:SRR070554] (control runs).
Total RNA from mycelial samples was isolated using RiboPure™-Yeast Kit (Ambion) according to the manufacturer's instructions. Total RNA was denatured (15 minutes, 60°C; 5% (v/v) formaldehyde, 50% (v/v) formamide, 40 mM MOPS, pH 7) and separated by agarose gel electrophoresis (1.2% agarose, 40 mM MOPS, 10 mM sodium acetate, 2 mM EDTA, 2% (v/v) formaldehyde, pH 7). Blotting, hybridization and chemoluminescent signal detection were performed according to the manufacturer's instructions (DIG Application Manual for Filter Hybridization, Roche). Gel load and blot signal strength were quantified and normalized using Bio-Rad (Munich, Germany) Quantity One (v4.6.7) software.
For cultivation of A. benhamiae, medium was prepared as follows: 10 g/l keratin (MP Biomedicals Europe, Illkirch, France) was autoclaved in water and subsequently 20 mM potassium phosphate pH 5.5, 0.4 mM MgSO4, 77 mM NaCl, 5 mM glucose and 0.5% (v/v) SL-8 trace elements  were added. Microconidia obtained from A. benhamiae cultivated for 7 days on MAT agar at 30°C were used to inoculate shaking flasks at a final spore concentration of 106 per milliliter. After cultivation for 2 days at 200 rpm and 30°C, cultures were filtered through Miracloth (Calbiochem, Darmstadt, Germany) and the supernatant was centrifuged at 4,000 g for 20 minutes at 4°C. Secreted proteins were precipitated with 10% (w/v) trichloroacetic acid/6.5 mM DTT overnight at 4°C. The precipitate was pelleted at 4,000 g for 20 minutes at 4°C and resuspended twice in ice-cold acetone/water (9:1)/6.5 mM DTT followed by subsequent centrifugation steps. The air-dried pellet was dissolved in lysis buffer 3, as described . Immobiline DryStrips of 11 cm covering a pH range from 3 to 10 (GE Healthcare Life Sciences) were rehydrated overnight according to the manufacturer's instructions. Isoelectric focusing was carried out in an Ettan IPGphor II using a 0 to 1 kV gradient for 11 h, 1 to 8 kV for 3 h and finally 8 kV for 24 kVh. Afterwards, strips were incubated for 15 minutes in equilibration buffer (6 M urea, 2% (w/v) SDS, 75 mM Tris.Cl pH 8.8, 30% (v/v) glycerol) with 65 mM DTT, followed by an alkylation step of the proteins with 135 mM iodoacetamide in equilibration buffer under the same conditions. Separation of proteins by the second dimension was carried out using pre-cast Criterion gels (12.5% (w/v), Tris-HCl; Bio-Rad) according to the manufacturer's instructions. Proteins were visualized by Colloidal Coomassie Brilliant Blue G-250 staining .
Protein spots were excised from the gels and digested with sequencing-grade Trypsin (Promega, Mannheim, Germany) as described elsewhere . Eluted peptides were mixed with an equal amount of a saturated alpha-cyano-hydroxycinnamic acid (Bruker Daltonics, Bremen, Germany) solution in aqueous 30% (v/v) acetonitrile and spotted on an MTP anchor-chip 800/384 (Bruker Daltonics). Mass spectrometry spectra were acquired with an Ultraflex I TOF/TOF (Bruker Daltonics) mass spectrometer using Peptide Mass Standard II (Bruker Daltonics) as calibrant. The five most intense mass spectrometry signals were selected for tandem mass spectrometry analysis. MASCOT (version 2.1.02; Matrix Science, London, UK) searching against protein predictions from the A. benhamiae genome and the NCBI database (taxa fungi) was used for protein identification with the following the parameters: fixed modification of cysteine to S-carbamidomethyl derivatives, variable methionine oxidation, no missed cleavage and a peptide mass tolerance of 200 ppm.
PKS and NRPS domain architecture prediction
Generation of the phylogenetic tree
For genome-based phylogeny, 23 proteins from 28 fully sequenced fungal genomes were used for the reconstruction of the phylogenetic relationships of A. benhamiae and T. verrucosum (Additional file 3). The 23 ortholog groups were selected based on the KOG (clusters of orthologous groups for eukaryotes) assignments, as described by Xu et al. . Only KOGs without paralogs, that is, proteins represented by a single protein in a species, were taken into consideration. Five proteins from the publication of Xu et al.  were not confirmed as fulfilling this requirement. Thus, they were not included. The genome set selected for the survey was non-redundant, that is, we did not consider four closely related Candida species as well as six Saccharomyces species, but only representatives of each clade, that is, C. albicans and S. cerevisiae, respectively. By contrast, we included all available Pezizomycetes, since A. benhamiae and T. verrucosum presumably belong to this phylum. A representative of Zygomycota (Rhizopus oryzae) was used as an outgroup. The considered genomes were as follws. Eurotiomycetes: Arthroderma benhamiae, Trichophyton verrucosum, Aspergillus clavatus, Aspergillus flavus, Aspergillus fumigatus, Aspergillus nidulans, Aspergillus oryzae, Aspergillus terreus, Botrytis cinerea, Coccidioides immitis, Histoplasma capsulatum, Paracoccidioides brasiliensis, Sclerotinia sclerotiorum, Stagonospora nodorum, Uncinocarpus reesii. Sordariomycetes: Chaetomium globosum, Fusarium graminearum, Magnaporthe grisea, Neurospora crassa. Saccharomycotina: Candida albicans, Lodderomyces elongisporus, Saccharomyces cerevisiae. Taphrinomycotina: Schizosaccharomyces japonicus. Basidiomycota: Coprinus cinereus, Cryptococcus neoformans, Puccinia graminis, Ustilago maydis. Zygomycota: Rhizopus oryzae.
The protein sets for each KOG protein shared among the 28 genomes were collected. Each set was then aligned by ClustalX, and the conserved blocks were extracted using the Gblocks tool  with allowance of smaller final blocks (five amino acids) and gap positions within the final blocks using otherwise default parameters. The extracted blocks were concatenated for each species. The phylogenetic analysis was performed using PHYML  for the construction of the maximal likelihood tree, and PHYLIP for the construction the neighbor joining tree, with the Jones-Taylor-Thornton model of the amino acid substitution in both cases. The neighbor joining and maximal likelihood trees had identical architecture.
The phylogenetic trees for proteases and enzymes involved in SM production were obtained using PHYLIP for the construction of the neighbor joining tree, with the Jones-Taylor-Thornton model of the amino acid substitution.
acyl carrier protein domain
Dulbecco's Modified Eagle's medium
expressed sequence tag
ketoacyl reductase domain
beta-ketoacyl synthase domain
- MAT locus:
mating type locus
non-ribosomal peptide synthetase
- SG medium:
Sabouraud 2% glucose medium
We thank Nancy Hannwacker, Silke Förster, Christin Heinrich, Sophia Keller, Ingrid Richter, Maria Pötsch and Silke Steinbach (HKI) for their technical assistance and advice. We are very grateful to the electron microscopy center at the University Hospital Jena for electron microscopic photographs, Yvonne Gräser (Berlin) for providing strains and Christina Cuomo (Broad Institute) for helpful discussions. This research was supported by the 'Pakt für Forschung und Innovation' of the Free State of Thuringia and the Federal Ministry of Science and Technology (BMBF, Germany), the HKI, the DFG funded excellence graduate school Jena School for Microbial Communication (JSMC) and the International Leibniz Research School for Microbial and Biomolecular Interactions Jena (ILRS). The Swiss-Prot group is part of the Swiss Institute of Bioinformatics (SIB) and of the UniProt Consortium. Swiss-Prot group activities are supported by the Swiss Federal Government through the Federal Office of Education and Science and by the National Institutes of Health (NIH) grant 2 U01 HG02712-04. Additional support comes from the European Commission contract SLING (226073).
- Weitzman I, Summerbell RC: The dermatophytes. Clin Microbiol Rev. 1995, 8: 240-259.PubMedPubMed CentralGoogle Scholar
- Monod M: Secreted proteases from dermatophytes. Mycopathologia. 2008, 166: 285-294. 10.1007/s11046-008-9105-4.PubMedView ArticleGoogle Scholar
- White TC, Oliver BG, Gräser Y, Henn MR: Generating and testing molecular hypotheses in the dermatophytes. Eukaryot Cell. 2008, 7: 1238-1245. 10.1128/EC.00100-08.PubMedPubMed CentralView ArticleGoogle Scholar
- Giddey K, Monod M, Barblan J, Potts A, Waridel P, Zaugg C, Quadroni M: Comprehensive analysis of proteins secreted by Trichophyton rubrum and Trichophyton violaceum under in vitro conditions. J Proteome Res. 2007, 6: 3081-3092. 10.1021/pr070153m.PubMedView ArticleGoogle Scholar
- Liu T, Zhang Q, Wang L, Yu L, Leng W, Yang J, Chen L, Peng J, Ma L, Dong J, Xu X, Xue Y, Zhu Y, Zhang W, Yang L, Li W, Sun L, Wan Z, Ding G, Yu F, Tu K, Qian Z, Li R, Shen Y, Li Y, Jin Q: The use of global transcriptional analysis to reveal the biological and cellular events involved in distinct development phases of Trichophyton rubrum conidial germination. BMC Genomics. 2007, 8: 100-10.1186/1471-2164-8-100.PubMedPubMed CentralView ArticleGoogle Scholar
- Zaugg C, Monod M, Weber J, Harshman K, Pradervand S, Thomas J, Bueno M, Giddey K, Staib P: Gene expression profiling in the human pathogenic dermatophyte Trichophyton rubrum during growth on proteins. Eukaryot Cell. 2009, 8: 241-250. 10.1128/EC.00208-08.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhang W, Yu L, Yang J, Wang L, Peng J, Jin Q: Transcriptional profiles of response to terbinafine in Trichophyton rubrum. Appl Microbiol Biotechnol. 2009, 82: 1123-1130. 10.1007/s00253-009-1908-9.PubMedView ArticleGoogle Scholar
- Staib P, Zaugg C, Mignon B, Weber J, Grumbt M, Pradervand S, Harshman K, Monod M: Differential gene expression in the pathogenic dermatophyte Arthroderma benhamiae in vitro versus during infection. Microbiology. 2010, 156: 884-895. 10.1099/mic.0.033464-0.PubMedView ArticleGoogle Scholar
- Yamada T, Makimura K, Abe S: Isolation, characterization, and disruption of dnr1, the areA/nit-2-like nitrogen regulatory gene of the zoophilic dermatophyte, Microsporum canis. Med Mycol. 2006, 44: 243-252. 10.1080/13693780500410909.PubMedView ArticleGoogle Scholar
- Fachin AL, Ferreira-Nozawa MS, Maccheroni W, Martinez-Rossi NM: Role of the ABC transporter TruMDR2 in terbinafine, 4-nitroquinoline N-oxide and ethidium bromide susceptibility in Trichophyton rubrum. J Med Microbiol. 2006, 55: 1093-1099. 10.1099/jmm.0.46522-0.PubMedView ArticleGoogle Scholar
- Ferreira-Nozawa MS, Silveira HC, Ono CJ, Fachin AL, Rossi A, Martinez-Rossi NM: The pH signaling transcription factor PacC mediates the growth of Trichophyton rubrum on human nail in vitro. Med Mycol. 2006, 44: 641-645. 10.1080/13693780600876553.PubMedView ArticleGoogle Scholar
- Chermette R, Ferreiro L, Guillot J: Dermatophytoses in animals. Mycopathologia. 2008, 166: 385-405. 10.1007/s11046-008-9102-7.PubMedView ArticleGoogle Scholar
- Drouot S, Mignon B, Fratti M, Roosje P, Monod M: Pets as the main source of two zoonotic species of the Trichophyton mentagrophytes complex in Switzerland, Arthroderma vanbreuseghemii and Arthroderma benhamiae.. Vet Dermatol. 2009, 20: 13-18. 10.1111/j.1365-3164.2008.00691.x.PubMedView ArticleGoogle Scholar
- Kane J, Smitka C: Early detection and identification of Trichophyton verrucosum. J Clin Microbiol. 1978, 8: 740-747.PubMedPubMed CentralGoogle Scholar
- Fumeaux J, Mock M, Ninet B, Jan I, Bontems O, Léchenne B, Lew D, Panizzon RG, Jousson O, Monod M: First report of Arthroderma benhamiae in Switzerland. Dermatology. 2004, 208: 244-250. 10.1159/000077311.PubMedView ArticleGoogle Scholar
- Takahashi Y, Sano A, Takizawa K, Fukushima K, Miyaji M, Nishimura K: The epidemiology and mating behaviour of Arthroderma benhamiae var. erinacei in household four-toed hedgehogs (Atelerix albiventris) in Japan. Nippon Ishinkin Gakkai Zasshi. 2003, 44: 31-38. 10.3314/jjmm.44.31.View ArticleGoogle Scholar
- Grunewald S, Paasch U, Gräser Y, Glander HJ, Simon JC, Nenoff P: Scarring tinea profunda in the pubic area due to Trichophyton verrucosum.. Hautarzt. 2006, 57: 811-813. 10.1007/s00105-005-1066-7.PubMedView ArticleGoogle Scholar
- Broad Institute database: Dermatophyte Comparative Database. [http://www.broadinstitute.org/annotation/genome/dermatophyte_comparative/MultiHome.html]
- Makimura K, Tamura Y, Mochizuki T, Hasegawa A, Tajiri Y, Hanazawa R, Uchida K, Saito H, Yamaguchi H: Phylogenetic classification and species identification of dermatophyte strains based on DNA sequences of nuclear ribosomal internal transcribed spacer 1 regions. J Clin Microbiol. 1999, 37: 920-924.PubMedPubMed CentralGoogle Scholar
- Summerbell RC, Haugland RA, Li A, Gupta AK: rRNA gene internal transcribed spacer 1 and 2 sequences of asexual, anthropophilic dermatophytes related to Trichophyton rubrum. J Clin Microbiol. 1999, 37: 4005-4011.PubMedPubMed CentralGoogle Scholar
- Gräser Y, El Fari M, Vilgalys R, Kuijpers AF, De Hoog GS, Presber W, Tietz H: Phylogeny and taxonomy of the family Arthrodermataceae (dermatophytes) using sequence analysis of the ribosomal ITS region. Med Mycol. 1999, 37: 105-114.PubMedView ArticleGoogle Scholar
- Kurtz S, Phillippy A, Delcher AL, Smoot M, Shumway M, Antonescu C, Salzberg SL: Versatile and open software for comparing large genomes. Genome Biol. 2004, 5: R12-10.1186/gb-2004-5-2-r12.PubMedPubMed CentralView ArticleGoogle Scholar
- Xu J, Saunders CW, Hu P, Grant RA, Boekhout T, Kuramae EE, Kronstad JW, Deangelis YM, Reeder NL, Johnstone KR, Leland M, Fieno AM, Begley WM, Sun Y, Lacey MP, Chaudhary T, Keough T, Chu L, Sears R, Yuan B, Dawson TL: Dandruff-associated Malassezia genomes reveal convergent and divergent virulence traits shared with plant and human fungal pathogens. Proc Natl Acad Sci USA. 2007, 104: 18730-18735. 10.1073/pnas.0706756104.PubMedPubMed CentralView ArticleGoogle Scholar
- Rawlings ND, Barrett AJ, Bateman A: MEROPS: the peptidase database. Nucleic Acids Res. 2010, 38: D227-D233. 10.1093/nar/gkp971.PubMedPubMed CentralView ArticleGoogle Scholar
- Giddey K, Favre B, Quadroni M, Monod M: Closely related dermatophyte species produce different patterns of secreted proteins. FEMS Microbiol Lett. 2007, 267: 95-101. 10.1111/j.1574-6968.2006.00541.x.PubMedView ArticleGoogle Scholar
- Zaugg C, Jousson O, Léchenne B, Staib P, Monod M: Trichophyton rubrum secreted and membrane-associated carboxypeptidases. Int J Med Microbiol. 2008, 298: 669-682. 10.1016/j.ijmm.2007.11.005.PubMedView ArticleGoogle Scholar
- Sharpton TJ, Stajich JE, Rounsley SD, Gardner MJ, Wortman JR, Jordar VS, Maiti R, Kodira CD, Neafsey DE, Zeng Q, Hung CY, McMahan C, Muszewska A, Grynberg M, Mandel MA, Kellner EM, Barker BM, Galgiani JN, Orbach MJ, Kirkland TN, Cole GT, Henn MR, Birren BW, Taylor JW: Comparative genomic analyses of the human fungal pathogens Coccidioides and their relatives. Genome Res. 2009, 19: 1722-1731. 10.1101/gr.087551.108.PubMedPubMed CentralView ArticleGoogle Scholar
- Baldo A, Mathy A, Tabart J, Camponova P, Vermout S, Massart L, Maréchal F, Galleni M, Mignon B: Secreted subtilisin Sub3 from Microsporum canis is required for adherence to but not for invasion of the epidermis. Br J Dermatol. 2010, 162: 990-997. 10.1111/j.1365-2133.2009.09608.x.PubMedView ArticleGoogle Scholar
- Segers R, Butt TM, Carder JH, Keen JN, Kerry BR, Peberdy JF: The subtilisins of fungal pathogens of insects, nematodes and plants: Distribution and variation. Mycol Res. 1999, 103: 395-402. 10.1017/S0953756298007345.View ArticleGoogle Scholar
- Henderson IR, Nataro JP: Virulence functions of autotransporter proteins. Infect Immun. 2001, 69: 1231-1243. 10.1128/IAI.69.3.1231-1243.2001.PubMedPubMed CentralView ArticleGoogle Scholar
- da Silva BA, dos Santos ALS, Barreto-Bergter E, Pinto MR: Extracellular peptidase in the fungal pathogen Pseudallescheria boydii. Curr Microbiol. 2006, 53: 18-22. 10.1007/s00284-005-0156-1.PubMedView ArticleGoogle Scholar
- Woodfolk JA, Slunt JB, Deuell B, Hayden ML, Platts-Mills TA: Definition of a Trichophyton protein associated with delayed hypersensitivity in humans. Evidence for immediate (IgE and IgG4) and delayed hypersensitivity to a single protein. J Immunol. 1996, 156: 1695-701.PubMedGoogle Scholar
- Monod M, Léchenne B, Jousson O, Grand D, Zaugg C, Stöcklin R, Grouzmann E: Aminopeptidases and dipeptidyl-peptidases secreted by the dermatophyte Trichophyton rubrum. Microbiology. 2005, 151: 145-155. 10.1099/mic.0.27484-0.PubMedView ArticleGoogle Scholar
- Gupta AK, Ahmad I, Borst I, Summerbell RC: Detection of xanthomegnin in epidermal materials infected with Trichophyton rubrum. J Invest Dermatol. 2000, 115: 901-905. 10.1046/j.1523-1747.2000.00150.x.PubMedView ArticleGoogle Scholar
- Ajello L, Cheng SL: The perfect state of Trichophyton mentagrophytes. Sabouraudia. 1967, 5: 230-234.PubMedView ArticleGoogle Scholar
- Symoens F, Jousson O, Planard C, Fratti M, Staib P, Mignon B, Monod M: Molecular analysis and mating behaviour of the Trichophyton mentagrophytesspecies complex. Int J Med Microbiol. 2010, Google Scholar
- Fraser JA, Stajich JE, Tarcha EJ, Cole GT, Inglis DO, Sil A, Heitman J: Evolution of the mating type locus: Insights gained from the dimorphic primary fungal pathogens Histoplasma capsulatum, Coccidioides immitis and Coccidioides posadasii.. Eukaryot Cell. 2007, 6: 622-629. 10.1128/EC.00018-07.PubMedPubMed CentralView ArticleGoogle Scholar
- Li W, Metin B, White TC, Heitman J: Organization and evolutionary trajectory of the mating type (MAT) locus in the dermatophyte and dimorphic fungal pathogens. Eukaryotic Cell. 2010, 9: 46-58. 10.1128/EC.00259-09.PubMedPubMed CentralView ArticleGoogle Scholar
- Aimanianda V, Bayry J, Bozza S, Kniemeyer O, Perruccio K, Elluru SR, Clavaud C, Paris S, Brakhage AA, Kaveri SV, Romani L, Latgé JP: Surface hydrophobin prevents immune recognition of airborne fungal spores. Nature. 2009, 460: 1117-1121. 10.1038/nature08264.PubMedView ArticleGoogle Scholar
- Takashio M: Sexual reproduction of some Arthroderma and Nannizia on diluted sabouraud agar with and without salts. Mykosen. 1972, 15: 11-17.PubMedView ArticleGoogle Scholar
- Burmester A, Wöstemeyer J: Cloned mitochondrial DNA from the zygomycete Absidia glauca promotes autonomous replication in Saccharomyces cerevisiae. Curr Genet. 1986, 10: 435-441. 10.1007/BF00419870.PubMedView ArticleGoogle Scholar
- Erhard M, Hipler UC, Burmester A, Brakhage AA, Wöstemeyer J: Identification of dermatophyte species causing onychomycosis and tinea pedis by MALDI-TOF mass spectrometry. Exp Dermatol. 2008, 17: 356-361. 10.1111/j.1600-0625.2007.00649.x.PubMedView ArticleGoogle Scholar
- Schroeckh V, Scherlach K, Nützmann HW, Shelest E, Schmidt-Heck W, Schuemann J, Martin K, Hertweck C, Brakhage AA: Intimate bacterial-fungal interaction triggers biosynthesis of archetypical polyketides in Aspergillus nidulans. Proc Natl Acad Sci USA. 2009, 106: 14558-14563. 10.1073/pnas.0901870106.PubMedPubMed CentralView ArticleGoogle Scholar
- Parra G, Blanco E, Guigo R: GeneID in Drosophila. Genome Res. 2000, 10: 511-515. 10.1101/gr.10.4.511.PubMedPubMed CentralView ArticleGoogle Scholar
- Romualdi A, Felder M, Rose D, Gausmann U, Schilhabel M, Glöckner G, Platzer M, Sühnel J, GenColors: Annotation and comparative genomics of prokaryotes made easy. Methods Mol Biol. 2007, 395: 75-96.PubMedView ArticleGoogle Scholar
- Kent WJ: BLAT-the BLAST-like alignment tool. Genome Res. 2002, 12: 656-664.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhang M, Gish W: Improved spliced alignment from an information theoretic approach. Bioinformatics. 2006, 22: 13-20. 10.1093/bioinformatics/bti748.PubMedView ArticleGoogle Scholar
- Anders S, Huber W: Differential expression analysis for sequence count data. Genome Biol. 2010, 11: R106-10.1186/gb-2010-11-10-r106.PubMedPubMed CentralView ArticleGoogle Scholar
- Robinson MD, McCarthy DJ, Smyth GK: edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010, 26: 139-140. 10.1093/bioinformatics/btp616.PubMedPubMed CentralView ArticleGoogle Scholar
- Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Series B. 1995, 57: 289-300.Google Scholar
- Biebl H, Pfennig N: Growth yields of green sulfur bacteria in mixed cultures with sulfur and sulfate reducing bacteria. Arch Microbiol. 1978, 117: 9-16. 10.1007/BF00689344.View ArticleGoogle Scholar
- Kniemeyer O, Lessing F, Scheibner O, Hertweck C, Brakhage AA: Optimisation of a 2-D gel electrophoresis protocol for the human-pathogenic fungus Aspergillus fumigatus. Curr Genet. 2006, 49: 178-189. 10.1007/s00294-005-0047-9.PubMedView ArticleGoogle Scholar
- Neuhoff V, Arold N, Taube D, Ehrhardt W: Improved staining of proteins in polyacrylamide gels including isoelectric focusing gels with clear background at nanogram sensitivity using Coomassie Brilliant Blue G-250 and R-250. Electrophoresis. 1988, 9: 255-262. 10.1002/elps.1150090603.PubMedView ArticleGoogle Scholar
- Vödisch M, Albrecht D, Lessing F, Schmidt AD, Winkler R, Guthke R, Brakhage AA, Kniemeyer O: Two-dimensional proteome reference maps for the human pathogenic filamentous fungus Aspergillus fumigatus. Proteomics. 2009, 9: 1407-1415.PubMedView ArticleGoogle Scholar
- InterProScan. [http://www.ebi.ac.uk/Tools/InterProScan]
- NRPS-PKS: [http://www.nii.res.in/nrps-pks.html]
- Talavera G, Castresana J: Improvement of phylogenies after removing divergent and ambiguously aligned blocks from protein sequence alignments. Syst Biol. 2007, 56: 564-577. 10.1080/10635150701472164.PubMedView ArticleGoogle Scholar
- Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52: 696-704. 10.1080/10635150390235520.PubMedView ArticleGoogle Scholar
- PHYLIP: [http://evolution.genetics.washington.edu/phylip.html]
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.