- Research article
- Open Access
Analysis of the dermatophyte Trichophyton rubrum expressed sequence tags
BMC Genomics volume 7, Article number: 255 (2006)
Dermatophytes are the primary causative agent of dermatophytoses, a disease that affects billions of individuals worldwide. Trichophyton rubrum is the most common of the superficial fungi. Although T. rubrum is a recognized pathogen for humans, little is known about how its transcriptional pattern is related to development of the fungus and establishment of disease. It is therefore necessary to identify genes whose expression is relevant to growth, metabolism and virulence of T. rubrum.
We generated 10 cDNA libraries covering nearly the entire growth phase and used them to isolate 11,085 unique expressed sequence tags (ESTs), including 3,816 contigs and 7,269 singletons. Comparisons with the GenBank non-redundant (NR) protein database revealed putative functions or matched homologs from other organisms for 7,764 (70%) of the ESTs. The remaining 3,321 (30%) of ESTs were only weakly similar or not similar to known sequences, suggesting that these ESTs represent novel genes.
The present data provide a comprehensive view of fungal physiological processes including metabolism, sexual and asexual growth cycles, signal transduction and pathogenic mechanisms.
Dermatophytes, consisting of organisms in the Trichophyton, Epidermophyton, and Microsporum genera, are the primary etiologic pathogens of various dermatophytoses, such as tinea capitis, tinea corporis, tinea inguinalis, tinea manus, tinea unguium and tinea pedis. These infections are widespread and increasing in prevalence on a global scale. Indeed, in some geographic regions, dermatophyte infection is now considered a major public health concern.
Unlike other fungi, dermatophytes can cause infections in healthy, immune-competent individuals. Estimates suggest that 30 to 70% of adults are asymptomatic carriers of these fungi . T. rubrum is the most common superficial fungus, accounting for at least 60% of all superficial fungal infections in humans. This organism may remain viable in the environment for over six months, thus accounting for widespread infections. Transmission occurs most often from person to person , for example, by shedding of infected skin cells and hair and by direct body contact .
Diagnosis of dermatophyte infections relies on clinical presentation, requiring successful isolation and culture and microscopy. It will take 2 to 4 weeks to culture and pleomorphic growth can lead to misidentification. Sometimes the morphological phenotype is not very stable. Rapid diagnostic tests using current molecular methodologies have been slow to develop for the dermatophytes .
Furthermore, despite the availability of new systemic antifungal therapies, nail infections are particularly difficult to eradicate, presenting a 25 to 40% recurrence rate . Another increasing problem in clinical treatments is growing resistance to antifungal drugs. In the past decade, more and more cases of azole- and terbinafine-resistant fungal infections have been reported [6–8]. Ryder et al [9, 10] characterized at the molecular level the cause of the resistant phenotype of some clinical isolates, considering the resistance in some T. rubrum clinical isolates to terbinafine appears to be due to a single amino acid substitution in the gene squalene epoxidase.
T. rubrum is anthropophilic, specialized to infect humans, but rare animal infections have been reported. It is urease negative and vitamins independent. On the Bromcresol purple-(BCP)-milk solids glucose media, it shows restricted growth, no indicator color change for 10 to 14 days, then profuse growth with indicator change to purple (alkaline). Its teleomorph has not yet been found. Arthroconidia and mycelial fragments are the infectious agents. Excreted extracellular proteinase, alkaline serine proteinase, plays a role in dermatophytes growth and multiplication as well as in the inflammatory reaction [2, 11–13]. Thus, they are considered to be critical virulence factors for dermatophytes. There are many reports of the isolation and characterization of one or two proteases from an individual species of dermatophyte [14–18]. In T. rubrum, Some keratinases have been isolated and a subtilisin gene family were identified at the genetic level [19–22]. Jousson et al isolated a five-member secreted metalloproteases (MEP) family from genomic libraries of T. rubrum, T. mentagrophytes and M. canis. Further phylogenetic analysis revealed that the metalloproteases secreted by the three species in vitro are encoded by orthologous genes, strongly suggesting that the multiplication of an ancestral metalloprotease gene occurred prior to the dermatophyte species divergence . However, only a small number of proteases from dermatophytes have been characterized at the genetic level.
Taken together, the lack of effective diagnostic and treatment strategies, the large number of individuals that experience dermatophyte infections, and the economic consequences highlight deficiencies in the research efforts aimed at understanding dermatophyte. T. rubrum is the most common superficial fungus, inducing dermatophytoses in various parts of the human skin, and can also cause deeper infections such as kerions, abscesses and granulomas. Research examining the pathogenesis of T. rubrum in numerous skin infections is needed in order to develop novel therapeutic agents and identify potential drug targets.
In the present study, we constructed various T. rubrum cDNA libraries and obtained 11,085 unique expressed sequence tags (ESTs). By analyzing these ESTs, we were able to further elucidate the physiological activities in metabolism, signal transduction, sexual life cycle, pathogenesis and cell wall composition of T. rubrum. These sequences provide an unparalleled resource for the future understanding of this remarkable fungus. Furthermore T. rubrum represents an attractive model species with which to study dermatophytes and other pathogenic filamentous fungi.
Results and discussion
Expressed sequence tags sequence determination
EST sequences were produced from 10 different cDNA libraries. After trimming for the vector and poor quality segments, a total of 34,670 sequences were identified and 11,085 unique ESTs were isolated. Isolated ESTs included 3,816 contigs and 7,269 singletons. All the constructed libraries, culture conditions and sequencing statistics employed in this study are summarized in table 1. The average sequence length of these unique ESTs was 586 nucleotides (nt). Approximately 8,639 (78%) of the unique ESTs were longer than 400 nt.
In the synthesis of the cDNA first strand, incomplete cDNA synthesis and/or 5' truncation of mRNA transcripts could lead to an overestimation of the gene number identified from our analysis. Accuracy of the estimation absolutely depends upon the quality of the library. Among the 7,764 matched ESTs, 1,439 could be matched to reference genes if the subjected starting site was set to be the first coding amino acid; meanwhile 2,325 matching genes could be found if the subjected starting site was allowed to be any of the first 10 amino acids. These data confirm that the quality of the libraries is reliable.
Comparison to Non-redundant and Clusters of Orthologous Groups databases
Identified ESTs were compared with the GenBank non-redundant (NR) protein database. This comparison enabled us to assign putative functions to or find homologs from other organisms for 7,764 (70%) of the ESTs. The remaining 3,321 (30%) ESTs were only weakly similar or not similar to known sequences (E ≥ 1E-05), suggesting that they represent novel genes.
The comparative results for the ESTs with respect to the Clusters of Orthologous Groups (COGs) database, and its eukaryotic counterpart termed the KOGs database, are summarized in table 2. The ESTs identified herein contained a broad range of genes, predominantly encoding putative proteins involved in primary metabolism, gene expression, post-translation processes and cell structure. A significant proportion of the identified ESTs were matched to genes involved in transcription and signal transduction, suggesting that T. rubrum has an elaborate regulation system.
The COG collection currently consists of 138,458 proteins, which form 4873 COGs and comprise 75% of the 185,505 (predicted) proteins encoded in 66 genomes of unicellular organisms, including 3 Eukaryota, Saccharomyces cerevisiae, Schizosaccharomyces pombe and Encephalitozoon cuniculi. The eukaryotic orthologous groups (KOGs) include proteins from 7 eukaryotic genomes: three animals (Caenorhabditis elegans, Drosophila melanogaster and Homo sapiens), one plant, Arabidopsis thaliana, two fungi (Saccharomyces cerevisiae and Schizosaccharomyces pombe), and the intracellular microsporidian parasite Encephalitozoon cuniculi. The current KOG set consists of 4852 clusters of orthologs, which include 59,838 proteins, or ~54% of the analyzed eukaryotic 110,655 gene products . Because the RNA processing in eukaryotic cells is more elaborated, and eukaryotic cells have some structures that prokaryotic cells lack, for example, nuclear and mitochondria, more genes involved in RNA processing and modification [A], chromatin structure and dynamics [B], cell cycle control, cell division, chromosome partitioning [D], nuclear structure [Y], cytoskeleton [Z] and intracellular trafficking, secretion, and vesicular transport [U] were identified in KOGs than COGs. The difference between gene number of class energy production and conversion [C] was mainly due to the mitochondrial related genes. On the contrary, more genes involved in cell motility [N], mainly encoding secretory pathway and cell motivity related proteins, were identified in COGs than KOGs. Meanwhile, the COG/KOG result was also affected by the subtle difference in functional classification to function similar proteins, especially transporters.
Comparison of ESTs to those already available in GenBank
To date 3,749 T. rubrum ESTs have been released to GenBank by other research groups, most of which (3,711) were obtained in condition mimicking virulence by Zaugg and Monod et al. When we compared these 3,749 previously reported ESTs to the NR database and to our EST pool, we found that 3,112 of the previously described ESTs in GenBank had matching homologs in other organisms in the NR database (E <1E-05) and that 637 ESTs were only weakly similar or did not match known sequences. Among the 3,749 ESTs identified previously, 2,738 ESTs matched 1,325 unique ESTs in our data (E <1E-05, identity ranges from 72% to 100%), suggesting these data be somewhat redundant. The remaining 1,011 did not have matching homologs in our EST database, including 263 ESTs which had no homology whatsoever in the NR database. Meanwhile, 748 of the 1,011 ESTs obtained only in their research had some homology with sequences in the NR database (E <1E-05), while 620 ESTs appeared to be hypothetical proteins, predicted proteins or unnamed protein. Some secreted proteinases that were identified by Monod et al. were not identified in our study (i.e. subtilisin-like protease SUB4, dipeptidyl-peptidase IV, leucine aminopeptidase 1, leucine aminopeptidase 2 and a putative secreted metalloprotease 4), suggesting that the condition in their research may be a more suitable environment to induce expression of these proteins. The proteinases subtilisin-like protease SUB1 and subtilisin-like protease SUB5 were expressed not only in the condition they mimicked but also in YPD media (see Additional file 1). By comparison with previous released ESTs, our data supply the information about gene expression in different conditions, advancing the current knowledge of T. rubrum transcriptome.
Metabolism and secondary metabolism
A large percentage of T. rubrum genes were annotated in glycolysis and oxidative phosphorylation systems, as expected from what is known from other aerobic filamentous fungi. Genes corresponding to the citrate cycle enzymes and to components of complexes I (NADH-CoQ reductase), II (Succinate-CoQ reductase), III (CoQ-Cytochrome C reductase), and IV (Cytochrome C oxidase) were identified. The presence of these genes reflects the fungus' ability to perform complete aerobic pyruvate degradation and oxidative phosphorylation.
The pathways involved in metabolizing mono- and disaccharides, such as glucose, fructose, mannose, and sucrose, as well as polysaccharide starch, were best represented in the identified ESTs. There were also a large number of membrane transporters for saccharides such as xylose and fucose. Among the identified ESTs, all aminoyl-tRNA synthases have been previously described with the exception of cyseine-tRNA synthase, of which only one contig was found to exhibit weak homology (C3495-Contig1, E value 0.00008). Comparison of our ESTs data to KEGG revealed that these proteins were probably involved in many amino acid metabolism pathways, including glutamate metabolism, alanine and aspartate metabolism, glycine, serine and threonine metabolism, methionine metabolism, cysteine metabolism, valine, leucine and isoleucine biosynthesis, lysine biosynthesis, arginine and proline metabolism, histidine metabolism, phenylalanine, tyrosine and tryptophan biosynthesis. Analysis of these pathways indicated that T. rubrum could itself synthesize lysine. It could synthesize some amino acids via amino acid synthetase including glutamate, glutamine, asparagine, cysteine, tryptophan and threonine, and it can also synthesize some other amino acids via amino transfer reactions. The primary lipids in T. rubrum are sterol and phospholipid. Our analysis of the metabolic pathways in T. rubrum revealed that it can also synthesize co-enzymes such as riboflavin, nicotinate, nicotinamide, coenzyme A, ubiquinone, and folate. The genes probably involved in the synthesis of thiamine included: NMT1 protein, THI4 protein, ThiC protein and a putative ThiG protein. Because T. rubrum's vitamin independence serves to separate it from species such as Trichophyton violaceum that have vitamin requirements, identification of the ESTs involved in vitamin synthesis is of great significance for molecular diagnosis. Some ESTs probably involved in metabolism are included in additional files (see Additional file 2).
Secondary metabolites are a remarkably diverse class of cellular products that often exhibit taxonomic specificity. Secondary metabolites are generally considered "nonessential" for organismal growth in culture. In addition to the metabolic pathways mentioned above, there are also several secondary metabolic pathways present in T. rubrum. It has long been known that dermatophytes can survive in toxic environments, presumably through the production of biological tools to resist these toxins. Indeed dermatophytes have even been isolated from sewer water and polluted soil containing organic compounds such as aromatic compounds. The pathways implicated by the EST analysis are consistent with T. rubrum's capability to enzymatically degrade various toxins such as phenanthrene, dibenzofuran, ethylbenzene, styrene, fluorene, and 1,1,1-Trichloro-2,2-Bis-(4'-Chlorophenyl)Ethane(DDT) (see Additional file 3).
Similar to other pathogenic filamentous fungi, such as Magnaorthe grisea  and Aspergillus fumigatus , several important secondary metabolic products, were also found in T. rubrum; these included non-ribosomal peptide synthases, polyketide synthases, two putative dimethylallyl tryptophan synthases, a putative arsenate reductase and a hydrophobin. These gene products may be related to T. rubrum growth and pathogenicity. Polyketides (derived from polyketones) are a class of secondary metabolites produced by most organisms, but they have been most extensively examined in bacteria and fungi. In fungi, numerous functions have been proposed for polyketides, including the production of toxins [27, 28] and spore pigments [29, 30]. Although it is well known that T. rubrum can produce pigments, the relationship between secondary metabolic pathways and pigments production remains unresolved. The various identified ESTs involved in secondary metabolism are listed in additional files (see Additional file 3).
Interestingly, T. rubrum harbors a putative sterigmatocystin biosynthesis monooxygenase StcW (C1113-Contig1), and a probable sterigmatocystin biosynthesis P450 monooxygenase STCL (Cytochrome P450 60B, EST000637, [GenBank:DW406216]), suggesting that it very likely may produce sterigmatocystin. However, confirmation of this possibility will require further investigation. Sterigmatocystin is of particular significance in evaluating the toxicity of T. rubrum products because prior evidence indicates that it is probably toxic to the human liver.
Although it is suspected that Trichophyton mentagrophytes can produce a penicillin-like substance , we did not find any evidence suggesting that T. rubrum generates antibiotic-like by-products.
We also identified a variety of signal transduction systems in T. rubrum, such as MAPK, cAMP-dependent pathways, G-protein pathways, Ras pathways, and a large number of serine/threonine protein kinase/phosphatases and signal histidine kinases. Although the first four systems are generally conserved among fungal and mammalian species, the numbers and functions of signal transduction of histidine kinases vary between fungi. For example, in S. cerevisiae there exists only one, but in Neurospora crassa  there are 11 different histidine kinases. Indeed, there are many differences in the numbers and types of histidine kinases between filamentous fungi and yeast. Some experimental results [33–36] obtained with Candida albicans and Aspergillus fumigatus suggest that histidine kinases are related to fungal pathogenesis. In T. rubrum many ESTs are homologous to histidine kinases, some of which may be related to osmotic and nutrient responses. Several identified histidine kinases contain PAS/PAC domains, suggesting that they are involved in oxygen and light responses.
In addition to histidine kinase-related ESTs, we also identified Calcium-Calmodulin homologues and related proteins, suggesting that Ca-CaM pathways may be utilized by T. rubrum, similar to the model organism Neurospora crassa. Commonly, Ca2+ release from internal stores is mediated by the second messengers inositol-1, 4, 5-trisphosphate (InsP3) and cADP ribose, or by Ca2+-induced Ca2+ release . Although InsP3 is present within Neurospora hyphae  and T. rubrum, Neurospora lack recognizable InsP3 receptors, ADP ribosyl cyclase and ryanodine receptor proteins, which are principal components of Ca2+ release mechanisms in plant and animal cells. These observations raise the question of whether there may be other second messenger systems that are responsible for Ca2+ release from internal stores that remain to be discovered in filamentous fungi. InsP3 receptors, ADP ribosyl cyclase or ryanodine receptor proteins were also not found in our EST data (the existence of these genes can only be determined after the whole genome is sequenced). Thus the mechanism of T. rubrum Ca-CaM pathways remains to be determined. The various identified ESTs involved in signal transduction are listed in additional files (see Additional file 4).
Sexual and asexual development
Although the sexual life cycle has been described for Trichophyton mentagrophytes, Trichophyton ajelloi and Microsporum canis [2, 39], which are all heterothallic species, the sexual cycle of T. rubrum remains to be elucidated. We compared sexual-cycle related genes with those in other fungi, and our findings suggest that T. rubrum may be capable of sexual reproduction.
Sexual reproduction in ascomycete filamentous fungi is governed, in part, by two different mating-type genes that establish sexual compatibility: one gene encodes a protein with a high mobility group (HMG) domain, and the other encodes a protein with an alpha box domain (MAT). Such MAT containing genes are termed alpha mating-type genes. Homothallic fungi typically possess both mating-type genes and are self-fertilizing. Heterothallic fungi possess only one mating-type gene and require a partner with a different complementary mating-type gene. Although the MAT locus has not yet been found in T. rubrum, HMG is present (4 ESTs: EST000702, [GenBank:DW406281]; EST001404, [GenBank:DW406983]; EST004330, [GenBank:DW680849]; EST004146, [GenBank:DW680665]), and a large number of meiosis-related genes were identified in the ESTs, many of which were also found in Aspergillus niduluns and/or A. fumigatus. An analysis of 215 genes implicated in the fungal mating process, pheromone response, meiosis and fruiting body development revealed that many genes present in A. nidulans [26, 40, 41] are also present in T. rubrum, including several genes for which the only known function is related to sexual reproduction (see Additional file 5). These results suggest that T. rubrum may also possess sexual cycles. However more in-depth exploration will be required in order to determine whether they indeed possess sexual cycles.
The current research on the fungal asexual development cycle primarily involves studies of Neurospora crassa and Aspergillus niduluns, and the literature contains two distinct models. Our comparison of the identified T. rubrum ESTs with the NR protein database, revealed a key enzyme (FlbD, EST023258, [GenBank:DW699777], E value 5E-52) that is present in A. niduluns , but no key enzymes present in Neurospora crassa. Because in A. nidulans four proteins (FlbC, BrlA, AbaA and WetA) besides FlbD are also required in macroconidiation pathway, it is still unclear whether T. rubrum is more similar to A. niduluns or N. crassa with regards to asexual development.
The most evident dermatophyte feature is the ability to digest keratin. Dermatophytes can degrade human and other animal keratin protein and utilize it. This represents the pathogenic feature that differentiates dermatophytes from other fungi. The secreted proteinases of dermatophytes play an important role in the process of infection and are thus considered the primary virulent factors. However, only a few proteinase sequences have been identified thus far in T. rubrum [19–22, 42–45]. In our data, secreted proteins were firstly analyzed by the NR comparison result. To the putative secreted protein identified in NR database, all of the corresponding ESTs were further subjected to SignalP prediction analysis. As to the ESTs listed in the Additional file 6, 48 of them were predicted to contain a signal peptide. Among the T. rubrum ESTs, identified secretory systems include type I, II, III and V, with type II, III being the most heavily represented. We also identified a large number of putative secretion related pathways and putative secreted proteinase, including some serine proteinases, aspartic protease, alkaline proteinase, peptidases and metalloprotease. (see Additional file 6). Among the serine proteinases, 10 unique ESTs were homologous to the known dermatophyte subtilisin-like serine protease family members SUB1, SUB5, and SUB6 ; meanwhile no identified ESTs were homologous to SUB2, SUB3, SUB4, or SUB7. The identity of the 6 serine protease ESTs with highest homology to T. rubrum serine proteinases varied from 80% to 93%, suggesting that the genes may be selectively spliced, and/or the family may have other members, and/or that these genes have many copies in the genome.
Secreted metalloproteases are thought to be associated with lesion extension . Recent reports  have identified a nitrogen regulating factor response region located upstream of the Microsporum canis MEP1 gene, consistent with the finding that, at the beginning of infection, there are only a few urea molecules, amino acids and glucose molecules available as a nutrient resource in sweat. Aspartic protease is an acid proteinase whereas secreted metal proteinase and subfamily members are neutral or basic proteinases. The presence of these proteinases supports the finding that initially T. rubrum can grow in either acidic or basic pH environments.
These proteinase and peptidase findings suggest that they may also be related to amino acid transportation. During T. rubrum infection, extracellular proteins are hydrolyzed into peptides by secreted proteinases then are further degraded into amino acids or dipeptides by peptidases and finally transported into cells. The identification of these proteinases and peptidases increases our understanding of the pathogenic mechanism underlying T. rubrum infection.
These secreted proteinases can not only degrade proteins such as keratin, elastin and collagen to supply nutrients to the fungi, but can also induce delayed-type hypersensitivity (such as with the SUB family and Tri r4) [1, 22]. However, although the T. rubrum allergen Tri r4 EST shares 98% homology with that in Trichophyton mentagrophytes, it is clear that the inflammation mediated by T. rubrum is not as severe as that induced by T. mentagrophytes . This is probably due, at least in part, to the presence of multiple inflammatory factors in T. mentagrophytes. When comparing the putative virulence factors present in T. rubrum with those in A. fumigatus  we also observed that many known virulence factors in A. fumigatus were not present in the identified T. rubrum ESTs. Likewise, many putative T. rubrum virulence factors were not found or were divergent in A. fumigatus. For example, the mep 1 gene exists in both A. fumigatus (Afu8g07080) and dermatophytes, but the two sequences were highly divergent (E value 5e-26).
The cell wall is a structure that humans lack but that fungi have. Therefore, the cell wall represents an ideal target for novel anti-fungal drugs. Many cell wall-related proteins were found among the presently identified ESTs, including chitin synthesis, chitinase, β1,3-glucan synthase, β1,6-glucan synthase, and 1,4-alpha-glucan branching enzyme. By analyzing the pathways involved in cell wall synthesis, it can be inferred that T. rubrum is probably capable of synthesizing peptidoglycan, the primary component of gram-positive bacteria. The principle enzymes involved in the peptidoglycan synthesis pathway UDP-N-acetylmuramoylalanyl-D-glutamyl-2, 6-diaminopimelate-D-alanyl-D-alanine ligase (EST020893, [GenBank:DW697412], GO:0008766) and phospho-N-acetylmuramoyl-pentapeptide-transferase (EST009048, [GenBank:DW685567], EC:22.214.171.124) were identified in our EST data. But it is still to be determined by experimental confirmation. We also identified some putative proteins that are likely to be involved in sterol synthesis, which is a primary target for clinically-available drugs. The components of the identified ESTs that are involved in cell wall synthesis are listed in table 3.
The estimated T. rubrum genome size is at least 22.05 Mb . And we estimate that a significant proportion of the gene content of T. rubrum is represented in this collection of sequences. Although ESTs can only represent genes that are actively expressed in particular phases, they complement the value of genomic sequencing through the functional identification of novel genes and provide information about gene structure and expression patterns [48–53]. Thirty percent of the 11,085 contigs and singletons identified in our analysis represent unique genes. Our Blast searches revealed that 33% of the 11,085 unique sequences possess matches in the yeast genome (<e-5), suggesting that there are clear distinctions (~67%) between yeast and filamentous fungi. These findings highlight the need for additional research on filamentous fungi.
In summary, although superficial fungi are the primary human pathogenic fungi, our knowledge about these organisms remains limited. T. rubrum represents an ideal model for the study of superficial fungi; our work in identifying ESTs in T. rubrum cDNA libraries will facilitate a greater understanding of the molecular mechanisms underlying its growth, metabolism, pathogenesis and drug resistance. In addition, our work may aid in the identification of novel effective drug targets and anti-fungal agents.
Preparation of strains and materials
T. rubrum (strain BMU01672) and medium have been described previously .
Preparation of cell cultures
Preparation of hyphae cultures
A few hyphae of T. rubrum were inoculated on potato glucose agar and incubated at 27.5°C for 2~3 weeks. In order to get cDNAs covering the entire growth phases, the hyphae were inoculated in liquid YPG medium (10 g/L yeast extract, 20 g/L peptone, 10 g/L D-glucose) and incubated in a 27.5°C bath shaker for 7, 10, 14, 15, 16, 20, 22, 26, 28, 34 and 36 d, respectively. Combine the mycelium from (7,10, 15, 20, 28) (7, 10, 16, 22, 26, 28), (16, 20, 22), (26, 28), and (34, 36) days, thus get 8 samples (or mixtures) from 7, 10, 14, (7, 10, 15, 20, 28), (7, 10, 16, 22, 26, 28), (16, 20, 22), (26, 28),(34, 36) day. All the samples or mixtures were used to construct unnormalized cDNA libraries except the mixture (7, 10, 15, 20, 28), which was used to prepare the tester cDNA of a subtracted cDNA library. A large number of sequenced plasmids in 7, 10, and 14 days libraries were in vitro transcripted (Riboprobe® in vitro Transcription Systems, Promega) and the RNA products were further reverse transcripted into the driver cDNA (Clontech PCR-Selected cDNA Subtraction Kit). The fungal samples (or mixtures) were centrifuged and the supernatant were discarded. The pelleted mycelium were washed twice with PBS.
Preparation of spore cultures
A few hyphae of T. rubrum were inoculated on potato glucose agar and incubated at 28°C for 2~3 weeks. The spores were then washed by liquid YPG medium and filtered by cell filter to get rid of the hyphae. Part of the spores obtained were incubated in YPG medium at 28°C for 12 hours. Finally we got two samples, spores and incubated spores, respectively.
Isolation and purification of total RNA and mRNA
Total RNA and mRNA were isolated and purified as described previously .
Construction of the cDNA library
Construction of the unnormalized cDNA library
The 9 cDNA libraries (7 representing mycelium and 2 representing spores) were constructed following the protocols of the SUPERSCRIPT™ Plasmid System with GATEWAY™ Technology for cDNA Synthesis and Cloning (Invitrogen).
Construction of the subtracted cDNA library
The subtracted library was constructed following the protocols of the Clontech PCR-Selected cDNA Subtraction Kit. According to the instruction of the Kit, driver cDNA will be subtracted from the tester cDNA therefore a normalized library is constructed.
Isolation and purification of plasmids, and sequencing of clones
The cDNA plasmids were isolated as described previously . Sequencing was performed with a generic T7 primer located 5' upstream of the inserted segments, following the protocol of the PRISM Big Dye Terminator Kit on an ABI3700 or MEGABACE automated sequencer.
Sequence data processing, EST clustering and re-assembly
The processing method has been described previously .
EST analysis and construction of the metabolic pathway
The clustered EST consensus sequences were assigned with potential functions through homologous comparisons by BLASTX searches of the GenBank non-redundant (NR) protein database (07/2005), Gene Ontology (GO, 11/2005) and Yeast Genome Database (YGD, 11/2005). ESTs were further classified according to the NCBI Clusters of Orthologous Groups of Proteins database (COGs, 06/2005) and the Eukaryotic Orthologous Groups (KOGs, 06/2005). The metabolic pathways of T. rubrum were partially reconstructed by searching for known pathway homologs found in the Kyoto Encyclopedia of Genes and Genomes database (KEGG, 11/2004). Secreted proteins were analyzed by comparison to NR. All the threshold cutoff were E<1E-05. As to the putative secreted protein identified in NR database, all of the corresponding ESTs were further subjected to SignalP prediction analysis.
The sequences have been submitted to GenBank (accession numbers from [GenBank:DW405580] to [GenBank:DW407270] and from [GenBank:DW678211] to [GenBank:DW711189]). All the information about GenBank numbers, contigs assembly, and ESTs annotation could be obtained form our Trichophyton rubrum database .
Woodfolk JA, Wheatley LM, Piyasena RV, Benjamin DC, Platts-Mills TE: Trichophyton antigens associated with IgE antibodies and delayed type hypersensitivity. Sequence homology to two families of serine proteinases. J Biol Chem. 1998, 273: 29489-29496. 10.1074/jbc.273.45.29489.
Kane J, Summerbell R, Sigler L, Krajden S, Land G: Laboratory Handbook of Dermatophytes. 1997, Star Publishing Company
Kwon-Chung KJ, Bennett JE: Medical Mycology. 1992, Lea & Febiger, Philadelphia
Kac G: Molecular approaches to the study of dermatophytes. Med Mycol. 2000, 38: 329-36.
Hay RJ: The future of onychomycosis therapy may involve a combination of approaches. Br J Dermatol. 2001, 145 (Suppl 60): 3-8.
Mukherjee PK, Leidich SD, Isham N, Leitner I, Ryder NS, Ghannoum MA: Clinical Trichophyton rubrum Strain Exhibiting Primary Resistance to Terbinafine. Antimicrob Agents Ch. 2003, 47: 82-86. 10.1128/AAC.47.1.82-86.2003.
Fernández-Torres B, Carrillo AJ, Martín E, Del Palacio A, Moore MK, Valverde A, Serrano M, Guarro J: In Vitro Activities of 10 Antifungal Drugs against 508 Dermatophyte Strains. Antimicrob Agents Ch. 2001, 45: 2524-2528. 10.1128/AAC.45.9.2524-2528.2001.
Pujol I, Capilla J, Ferna'ndez-Torres B, Ortoneda M, Guarro J: Use of the Sensititre Colorimetric Microdilution Panel for Antifungal Susceptibility Testing of Dermatophytes. J Clin Microbiol. 2002, 40: 2618-2621. 10.1128/JCM.40.7.2618-2621.2002.
Osborne CS, Leitner I, Favre B, Ryder NS: Amino Acid Substitution in Trichophyton rubrum Squalene Epoxidase Associated with Resistance to Terbinafine. Antimicrob Agents Ch. 2005, 49: 2840-2844. 10.1128/AAC.49.7.2840-2844.2005.
Osborne CS, Leitner I, Hofbauer B, Fielding CA, Favre B, Ryder NS: Biological, Biochemical, and Molecular Characterization of a New Clinical Trichophyton rubrum Isolate Resistant to Terbinafine. Antimicrob Agents Ch. 2006, 50: 2234-2236. 10.1128/AAC.01600-05.
Sanyal AK, Das SK, Banerjee AB: Purification and partial characterization of an exocellular proteinase from Trichophyton rubrum. Sabouraudia. 1985, 23: 165-178.
Collins JP, Grappel SF, Blank F: Role of keratinises in dermatophytosis. II. Fluorescent antibody studies with keratinase II of Trichophyton mentagrophytes. Dermatologica. 1973, 146: 95-100.
Grappel SF, Blank F: Role of keratinases in ermatophytosis. I. Immune responses of guinea pigs infected ith Trichophyton mentagrophytes and guinea pigs immunized with eratinases. Dermatologica. 1972, 145: 245-255.
Brouta F, Descamps F, Fett T, Losson B, Gerday C, Mignon B: Purification and characterization of a 43.5 kDa keratinolytic metalloprotease from Microsporum canis. Med Mycol. 2001, 39: 269-275.
Brouta F, Descamps F, Monod M, Vermout S, Losson B, Mignon B: Secreted metalloprotease gene family of Microsporum canis. Infect Immun. 2002, 70: 5676-5683. 10.1128/IAI.70.10.5676-5683.2002.
Tsuboi R, Ko IJ, Takamori K, Ogawa H: Isolation of a keratinolytic proteinase from Trichophyton mentagrophytes with enzymatic activity at acidic pH. Infect Immun. 1989, 57: 3479-3483.
Descamps F, Brouta F, Monod M, Zaugg C, Baar D, Losson B, Mignon B: Isolation of a Microsporum canis gene family encoding three subtilisin-like proteases expressed in vivo. J Invest Dermatol. 2002, 119: 830-835. 10.1046/j.1523-1747.2002.01784.x.
Pierard G, Gerday C, Losson B: Purification and characterization of a 31.5 kDa keratinolytic subtilisin-like serine protease from Microsporum canis and evidence of its secretion in aturally infected cats. Med Mycol. 1998, 36: 395-404. 10.1046/j.1365-280X.1998.00179.x.
Apodaca G, McKerrow JH: Purification and characterization of a 27,000-Mr extracellular proteinase from Trichophyton rubrum. Infect Immun. 1989, 57: 3072-3080.
Asahi M, Lindquist R, Fukuyama K, Apodaca G, Epstein WL, McKerrow JH: Purification and characterization of major extracellular proteinases from Trichophyton rubrum. Biochem J. 1985, 232: 139-144.
Lambkin I, Hamilton AJ, Hay RJ: Purification and characterisation of a novel 34,000-Mr cell-associated proteinase from the dermatophyte Trichophyton rubrum. FEMS Immunol Med Microbiol. 1996, 13: 131-140.
Joussona O, Le'chennea B, Bontemsa O, Mignonb B, Reichardc U, Barbland J, Quadronid M, Monod M: Secreted subtilisin gene family in Trichophyton rubrum. Gene. 2004, 339: 79-88. 10.1016/j.gene.2004.06.024.
Jousson O, Le' chenne B, Bontems O, Capoccia S, Mignon B, Barblan J, Quadroni M, Monod M: Multiplication of an ancestral gene encoding secreted fungalysin preceded species differentiation in the dermatophytes Trichophyton and Microsporum. Microbiology. 2004, 150: 301-310. 10.1099/mic.0.26690-0.
Tatusov RL, Fedorova ND, Jackson JD, Jacobs AR, Kiryutin B, Koonin EV, Krylov DM, Mazumder R, Mekhedov SL, Nikolskaya AN, Rao BS, Smirnov S, Sverdlov AV, Vasudevan S, Wolf YI, Yin JJ, Natale DA: The COG database: an updated version includes eukaryotes. BMC Bioinformatics. 2003, 4: 41-10.1186/1471-2105-4-41.
Dean RA, Talbot NJ, Ebbole DJ, Farman ML, Mitchell TK, Orbach MJ, Thon M, Kulkarni R, Xu JR, Pan H, Read ND, Lee YH, Carbone I, Brown D, Oh YY, Donofrio N, Jeong JS, Soanes DM, Djonovic S, Kolomiets E, Rehmeyer C, Li W, Harding M, Kim S, Lebrun MH, Bohnert H, Coughlan S, Butler J, Calvo S, Ma LJ, Nicol R, Purcell S, Nusbaum C, Galagan JE, Birren BW: The genome sequence of the rice blast fungus Magnaporthe grisea. Nature. 2005, 434: 980-986. 10.1038/nature03449.
Nierman WC, Pain A, Anderson MJ, Wortman JR, Kim HS, Arroyo J, Berriman M, Abe K, Archer DB, Bermejo C, Bennett J, Bowyer P, Chen D, Collins M, Coulsen R, Davies R, Dyer PS, Farman M, Fedorova N, Fedorova N, Feldblyum TV, Fischer R, Fosker N, Fraser A, Garcia JL, Garcia MJ, Goble A, Goldman GH, Gomi K, Griffith-Jones S, Gwilliam R, Haas B, Haas H, Harris D, Horiuchi H, Huang J, Humphray S, Jimenez J, Keller N, Khouri H, Kitamoto K, Kobayashi T, Konzack S, Kulkarni R, Kumagai T, Lafon A, Latge JP, Li W, Lord A, Lu C, Majoros WH, May GS, Miller BL, Mohamoud Y, Molina M, Monod M, Mouyna I, Mulligan S, Murphy L, O'Neil S, Paulsen I, Penalva MA, Pertea M, Price C, Pritchard BL, Quail MA, Rabbinowitsch E, Rawlins N, Rajandream MA, Reichard U, Renauld H, Robson GD, Rodriguez de Cordoba S, Rodriguez-Pena JM, Ronning CM, Rutter S, Salzberg SL, Sanchez M, Sanchez-Ferrero JC, Saunders D, Seeger K, Squares R, Squares S, Takeuchi M, Tekaia F, Turner G, Vazquez de Aldana CR, Weidman J, White O, Woodward J, Yu JH, Fraser C, Galagan JE, Asai K, Machida M, Hall N, Barrell B, Denning DW: Genomic sequence of the pathogenic and allergenic filamentous fungus Aspergillus fumigatus. Nature. 2005, 438: 1151-1156. 10.1038/nature04332.
Yu JH, Leonard TJ: Sterigmatocystin biosynthesis in Aspergillus nidulans requires a novel type I polyketide synthase. J Bacteriol. 1995, 177: 4792-4800.
Feng GH, Leonard TJ: Characterization of the polyketide synthase gene (pksL1) required for aflatoxin biosynthesis in Aspergillus parasiticus. J Bacteriol. 1995, 177: 6246-6254.
Linnemannstons P, Schulte J, del Mar Prado M, Proctor RH, Avalos J, Tudzynski B: The polyketide synthase gene pks4 from Gibberella fujikuroi encodes a key enzyme in the biosynthesis of the red pigment bikaverin. Fungal Genet Biol. 2002, 37: 134-148. 10.1016/S1087-1845(02)00501-7.
Mayorga ME, Timberlake WE: The developmentally regulated Aspergillus nidulans wA gene encodes a polypeptide homologous to polyketide and fatty acid synthases. Mol Gen Genet. 1992, 235: 205-212. 10.1007/BF00279362.
Fischer JB, Kane J: The detection of contamination in Trichophyton rubrum and T. mentagrophytes. Mycopathol Mycol Appl. 1971, 43: 169-180. 10.1007/BF02051718.
Galagan JE, Calvo SE, Borkovich KA, Selker EU, Read ND, Jaffe D, FitzHugh W, Ma LJ, Smirnov S, Purcell S, Rehman B, Elkins T, Engels R, Wang S, Nielsen CB, Butler J, Endrizzi M, Qui D, Ianakiev P, Bell-Pedersen D, Nelson MA, Werner-Washburne M, Selitrennikoff CP, Kinsey JA, Braun EL, Zelter A, Schulte U, Kothe GO, Jedd G, Mewes W, Staben C, Marcotte E, Greenberg D, Roy A, Foley K, Naylor J, Stange-Thomann N, Barrett R, Gnerre S, Kamal M, Kamvysselis M, Mauceli E, Bielke C, Rudd S, Frishman D, Krystofova S, Rasmussen C, Metzenberg RL, Perkins DD, Kroken S, Cogoni C, Macino G, Catcheside D, Li W, Pratt RJ, Osmani SA, DeSouza CP, Glass L, Orbach MJ, Berglund JA, Voelker R, Yarden O, Plamann M, Seiler S, Dunlap J, Radford A, Aramayo R, Natvig DO, Alex LA, Mannhaupt G, Ebbole DJ, Freitag M, Paulsen I, Sachs MS, Lander ES, Nusbaum C, Birren B: The genome sequence of the filamentous fungus Neurospora crassa. Nature. 2003, 422: 859-868. 10.1038/nature01554.
Clemons KV, Miller TK, Selitrennikoff CP, Stevens DA: A putative histidine kinase as a virulence factor for systemic aspergillosis. Med Mycol. 2002, 40: 259-262.
Pott GB, Miller TK, Bartlett JA, Palas JS, Selitrennikoff CP: The isolation of FOS-1, a gene encoding a putative two-component histidine kinase from Aspergillus fumigatus. Fungal Genet Biol. 2000, 31: 55-67. 10.1006/fgbi.2000.1225.
Selitrennikoff CP, Alex L, Miller TK, Clemons KV, Simon MI, Stevens DA: COS1, a putative two-component histidine kinase of Candida albicans, is an vivo virulence factor. Med Mycol. 2001, 39: 69-74.
Torosantucci A, Chiani P, De Bernardis F, Cassone A, Calera JA, Calderone R: Deletion of the two-component histidine kinase gene (CHK1) of Candida albicans contributes to enhanced growth inhibition and killing by human neutrophils in vitro. Infect Immun. 2002, 70: 985-987. 10.1128/IAI.70.2.985-987.2002.
Bootman MD, Collins TJ, Peppiatt CM, Prothero LS, MacKenzie L, De Smet P, Travers M, Tovey SC, Seo JT, Berridge MJ, Ciccolini F, Lipp P: Calcium signalling–an overview. Semin Cell Dev Biol. 2001, 12: 3-10. 10.1006/scdb.2000.0211.
Lakin-Thomas PL: Effects of inositol starvation on the levels of inositolphosphatesandinositol lipids in Neurospora crassa. Biochem J. 1993, 292: 805-811.
Weitzman I, Summerbell R: The Dermatophytes. Clin Microbiol Rev. 1995, 8 (2): 240-259.
Machida M, Asai K, Sano M, Tanaka T, Kumagai T, Terai G, Kusumoto K, Arima T, Akita O, Kashiwagi Y, Abe K, Gomi K, Horiuchi H, Kitamoto K, Kobayashi T, Takeuchi M, Denning DW, Galagan JE, Nierman WC, Yu J, Archer DB, Bennett JW, Bhatnagar D, Cleveland TE, Fedorova ND, Gotoh O, Horikawa H, Hosoyama A, Ichinomiya M, Igarashi R, Iwashita K, Juvvadi PR, Kato M, Kato Y, Kin T, Kokubun A, Maeda H, Maeyama N, Maruyama J, Nagasaki H, Nakajima T, Oda K, Okada K, Paulsen I, Sakamoto K, Sawano T, Takahashi M, Takase K, Terabayashi Y, Wortman JR, Yamada O, Yamagata Y, Anazawa H, Hata Y, Koide Y, Komori T, Koyama Y, Minetoki T, Suharnan S, Tanaka A, Isono K, Kuhara S, Ogasawara N, Kikuchi H: Genome sequenceing and analysis of Aspergillus oryzae. Nature. 2005, 438: 1157-1161. 10.1038/nature04300.
Galagan JE, Calvo SE, Cuomo C, Ma LJ, Wortman JR, Batzoglou S, Lee SI, Basturkmen M, Spevak CC, Clutterbuck J, Kapitonov V, Jurka J, Scazzocchio C, Farman M, Butler J, Purcell S, Harris S, Braus GH, Draht O, Busch S, D'Enfert C, Bouchier C, Goldman GH, Bell-Pedersen D, Griffiths-Jones S, Doonan JH, Yu J, Vienken K, Pain A, Freitag M, Selker EU, Archer DB, Penalva MA, Oakley BR, Momany M, Tanaka T, Kumagai T, Asai K, Machida M, Nierman WC, Denning DW, Caddick M, Hynes M, Paoletti M, Fischer R, Miller B, Dyer P, Sachs MS, Osmani SA, Birren BW: Sequencing of Aspergillus nidulans and comparative analysis with A. fumigatus and A. oryzae. Nature. 2005, 438: 1105-1115. 10.1038/nature04341.
Samdani AJ, Dykes PJ, Marks R: The proteolytic activity of strains of T. mentagrophytes and T. rubrum isolated from tinea pedis and tinea unguium infections. J Med Vet Mycol. 1995, 33: 167-170.
Apodaca G, McKerrow JH: Expression of proteolytic activity by cultures of Trichrophyton rubrum. J Med Vet Mycol. 1990, 28: 159-171.
Apodaca G, McKerrow JH: Regulation of Trichophyton rubrum proteolytic activity. Infect Immun. 1989, 57: 3081-3090.
Meevootisom V, Niederpruem DJ: Control of exocellular proteases in dermatophytes and especially Trichophyton rubrum. Sabouraudia. 1979, 17: 91-106.
Yamada T, Makimura K, Hirai A, Kano R, Hasegawa A, Uchida K, Yamaguchi H: Isolation of a region of a secreted metalloprotease gene from Microsporum canis. Jpn J Infect Dis. 2004, 57: 25-28.
Cervelatti EP, Ferreira-Nozawa MS, Aquino-Ferreira R, Fachin AL, Martinez-Rossi NM: Electrophoretic molecular karyotype of the dermatophytes Trichophyton rubrum. Genet Mol Biol. 2004, 27: 99-102. 10.1590/S1415-47572004000100016.
Goffeau A, Barrell BG, Bussey H, Davis RW, Dujon B, Feldmann H, Galibert F, Hoheisel JD, Jacq C, Johnston M, Louis EJ, Mewes HW, Murakami Y, Philippsen P, Tettelin H, Oliver SG: Life with 6000 genes. Science. 1996, 274: 546, 563-567. 10.1126/science.274.5287.546.
Posada-Buitrago ML, Frederick RD: Expressed sequence tag analysis of soybean rust pathogen Phakopsora pachyrhizi. Fungal Genet Biol. 2005, 42: 949-962.
Brown DW, Cheung F, Proctor RH, Butchko RA, Zheng L, Lee Y, Utterback T, Smith S, Feldblyum T, Glenn AE, Plattner RD, Kendra DF, Town CD, Whitelaw CA: Comparative analysis of 87,000 expressed sequence tags from the fumonisin-producing fungus Fusarium verticillioides. Fungal Genet Biol. 2005, 42: 848-861. 10.1016/j.fgb.2005.06.001.
Sacadura NT, Saville BJ: Gene expression and EST analysis of Ustilago maydis germination teliospores. Fungal Genet Biol. 2003, 40: 47-64. 10.1016/S1087-1845(03)00078-1.
Jones T, Federspiel NA, Chibana H, Dungan J, Kalman S, Magee BB, Newport G, Thorstenson YR, Agabian N, Magee PT, Davis RW, Scherer S: The diploid genome sequence of Candida albicans. Proc Natl Acad Sci. 2004, 101: 7329-7334. 10.1073/pnas.0401648101.
Braun BR, van Het Hoog M, d'Enfert C, Martchenko M, Dungan J, Kuo A, Inglis DO, Uhl MA, Hogues H, Berriman M, Lorenz M, Levitin A, Oberholzer U, Bachewich C, Harcus D, Marcil A, Dignard D, Iouk T, Zito R, Frangeul L, Tekaia F, Rutherford K, Wang E, Munro CA, Bates S, Gow NA, Hoyer LL, Kohler G, Morschhauser J, Newport G, Znaidi S, Raymond M, Turcotte B, Sherlock G, Costanzo M, Ihmels J, Berman J, Sanglard D, Agabian N, Mitchell AP, Johnson AD, Whiteway M, Nantel A: A human-curated annotation of the Candida albicans genome. PLoS Genet. 2005, 1: 36-57. 10.1371/journal.pgen.0010001.
Wang L, Ma L, Leng W, Yang J, Zhu J, Dong J, Xue Y, Zhe W, Li R, Jin Q: Analysis of part of the Trichophyton rubrum ESTs. Sci China C Life Sci. 2004, 47: 389-395. 10.1360/03yc0107.
Trichophyton rubrum database. [http://www.mgc.ac.cn/Trubrum/]
The work was supported by The National High Technology Research and Development Program of China (Accession number: 2001AA223021) and National Key Technologies R&D Programme (Accession Number: 2002BA711A14).
LW performed construction of the cDNA gene libraries, clone isolation, plasmids sequencing, data analysis and drafted the manuscript. LM, WL performed construction of the cDNA gene libraries, clone isolation, and plasmids sequencing. TL, LY participated in the construction of the cDNA libraries and the isolation of clones. JY carried out the bioinformatics analysis and constructed the T. rubrum database. LY, WZ, QZ, JD, YX, YZ and XX participated in the library construction, clone isolation and plasmids sequencing. ZW and RL were responsible for the strain identification, culture and growth conditions design. GD, FY, KT and YL participated in bioinformatics analysis. QJ designed the project, supervised the research and revised the manuscript. SY cooperated with QJ and implemented the project and supervise the research in the center of Chinese National Human Genome Center. All authors read and approved the final manuscript.
Lingling Wang, Li Ma, Wenchuan Leng contributed equally to this work.
Electronic supplementary material
About this article
Cite this article
Wang, L., Ma, L., Leng, W. et al. Analysis of the dermatophyte Trichophyton rubrum expressed sequence tags. BMC Genomics 7, 255 (2006) doi:10.1186/1471-2164-7-255
- Histidine Kinase
- Tinea Capitis
- Trichophyton Rubrum
- Trichophyton Mentagrophytes