- Research article
- Open Access
The use of global transcriptional analysis to reveal the biological and cellular events involved in distinct development phases of Trichophyton rubrum conidial germination
BMC Genomicsvolume 8, Article number: 100 (2007)
Conidia are considered to be the primary cause of infections by Trichophyton rubrum.
We have developed a cDNA microarray containing 10250 ESTs to monitor the transcriptional strategy of conidial germination. A total of 1561 genes that had their expression levels specially altered in the process were obtained and hierarchically clustered with respect to their expression profiles. By functional analysis, we provided a global view of an important biological system related to conidial germination, including characterization of the pattern of gene expression at sequential developmental phases, and changes of gene expression profiles corresponding to morphological transitions. We matched the EST sequences to GO terms in the Saccharomyces Genome Database (SGD). A number of homologues of Saccharomyces cerevisiae genes related to signalling pathways and some important cellular processes were found to be involved in T. rubrum germination. These genes and signalling pathways may play roles in distinct steps, such as activating conidial germination, maintenance of isotropic growth, establishment of cell polarity and morphological transitions.
Our results may provide insights into molecular mechanisms of conidial germination at the cell level, and may enhance our understanding of regulation of gene expression related to the morphological construction of T. rubrum.
Trichophyton rubrum is a worldwide pathogen causing various superficial infections, accounting for at least 60% of dermatophytosis, such as tinea capitis, tinea corporis, tinea inguinalis, tinea manus, tinea unguium and tinea pedis [1–4]. This pathogen, which normally causes well-characterized superficial infections, also produces skin infections in unusual parts of the body in immunodepressed patients [5–7]. T. rubrum has a relatively simple and non-sexual stage life-cycle. Similar to numerous other filamentous fungi, it forms a mycelial colony via hyphal tip growth, branching and hyphal fusion. T. rubrum produces numerous asexual multinucleate spores, called macroconidia or arthroconidia, which are considered to be the primary cause of skin and nail infections in the host. The infection is triggered by the adherence of T. rubrum arthrospores upon contact with the status corneum (horny layer) of the skin. The spores then germinate and the mycelium forms. The longitudinal growth of hyphae can penetrate into the deep part of the horny layer, while breadthwise growth can aggravate skin damage [8–10].
Eukaryotic genomes contain various non-coding regions, including regulatory components, introns and repeat sequences, making these genomes much larger and more complicated than those of bacteria. In recent years, more and more eukaryotic genomes have been completed or assembled. T. rubrum had been shortlisted for sequencing by the Fungal Genome Initiative (FGD white paper, 2004), yet expressed sequence tag (EST)-based mRNA transcription profiling and microarray technology are used extensively to characterize gene functions and discover functionally related genes needed for developmental and behavioural processes [11, 12]. T. rubrum consists of five chromosomes which range in size from 3.0 M bp to 5.8 M bp, and correspond to a total genome size at least 22.05 M bp . We commenced the ESTs project of T. rubrum in 2002 and part of ESTs of the genes were obtained by our group in 2003 . The project was completed recently. Extensive transcriptome data from eight life-cycle stages containing 11085 ESTs, which represent 80% of predicted genes, were obtained and annotated .
The prevalence of infections caused by T. rubrum and its human specific nature make it a good model for the study of human pathogenic filamentous fungi. The epidemiology, clinical case reports, strain relatedness, and drug susceptibilities of the organism are well documented . However, little is known about its genetic and biological characteristics for T. rubrum that cause the most common fungal infection in humans . In this study, we prepare a cDNA microarray consisting of 10250 ESTs to monitor the transcriptional pattern in a time serial process during T. rubrum germination. By functional analysis, we provide a global view of important biological system related to conidial germination. In addition, we hypothesize that the expression levels of genes related to some fundamental biological events, such as germination process, protein synthesis and other metabolic pathways activation, may be involved in the distinct steps of T. rubrum germination and sequential morphological transitions by impacting on distinguished cellular processes. Since signal transduction pathways and cell processes have been elucidated in the budding S. cerevisae, establishing a morphogenetic network that links cell-cycle events with cellular morphogenesis, we matched the ESTs sequences to the Saccharomyces Genome Database (SGD) and a number of homologues of S. cerevisae genes related to signalling transduction pathways and different cellular events, such as response to environmental conditions, maintenance of cell wall integrity and establishment of cell polarity were found to be involved in T. rubrum conidial germination. Our results at the cell level may provide stepping stones on the way to further study of the molecular mechanisms of T. rubrum conidial germination, and may enhance our understanding of some similarities between the morphogenetic networks of yeast and the filamentous fungi. Since it has been proposed that the identification and targeting of conidial germination-specific processes provides an excellent strategy for drug and fungicide development in these pathogenic fungi , this study may provide some clues for further approaches to finding new drug targets.
Morphological transitions during conidial germination
In T. rubrum, there are two visible morphological transitions in the course of germination. The first morphological transition occurred at 3–4 hr after incubation. In this stage, the conidia appear swollen and bright under microscopy. The second transition occurred at 9–10 hr after incubation. Most conidia budded and the germ tube emerged at the end of the conidia in this stage. The distinct morphological states in germination are shown in Figure 1.
Prominent transcriptional changes during T. rubrum conidial germination
In order to monitor the transcriptional changes that occur during T. rubrum conidial germination, we assembled an expression microarray that consists of 10250 ESTs. Using the microarray, the abundance of mRNA in samples collected at 1 hr intervals throughout the 15 hr time-course of conidial germination was evaluated. In the test procedure, we use genomic DNA as the reference control to normalize the relative expression levels under the mRNA hybridization signals. The genomic normalization procedure provided reproducible hybridization signals for 96% of the 10250 ESTs.
One-thousand five hundred seventy six putative genes which expression levels strongly changed during germination were selected to identify the correlation between expression profile and the events that occurred during germination. The 1576 genes were subjected to hierarchical clustering analysis using TIGR MultiExperiment Viewer (MeV) software . Five representative clusters were chosen visually, each of which showed distinctive profiles (Figure 2). The genes enrichment within cluster and the statistical significance of gene distribution across the clusters were estimated by post-hoc tests (detailed results see Additional files 1). The majority of genes (1561 out of the total of 1576 genes) belonged to one of the five chosen clusters. In order to verify the microarray result, the relative expression levels of 8 genes at several time points (0, 4, 10, and 15 hr) were estimated by Quantitative real-time RT-PCR. The results showed a strong positive correlation between the two techniques (detailed results see Additional file 1).
Cluster I contains 560 genes that showed maximum expression levels in dormant conidia, most of them with expression levels dropping quickly during conidial germination. A sub-cluster containing 80 genes shared a lowest expression at 4 hr but was induced at subsequent time-points. These data indicated that cluster I mRNAs are present in dormant conidia before the initiation of germination and must be stored in conidia before encapsulation. The 71 genes contained in cluster II were obviously induced during 0–2 hr. Their average expression levels were lower than those of genes in cluster IV and fluctuated several times throughout the process of germination. Cluster III contains 153 genes that had low and constant expression levels until 4 hr, but had increased expression levels subsequently. Cluster IV contains 404 genes and cluster V contains 373 genes. Transcripts of the genes in these two clusters were induced to increase from the beginning of germination. The genes in cluster V with a peak expression at 3–4 hr showed decreased levels of expression at the succeeding time-points, whereas the expression levels of genes within cluster IV were increased significantly during the first 3 hr but showed a higher or relatively constant trend thereafter.
Analysis of expression profiling associated with putative biological and physiological events during T. rubrum conidial germination
To identify the correlation between gene expression profile and physiological transitions during germination, ESTs of the selected 1576 genes were matched to the sequences in the Gene Ontology (GO) database by BLASTX searches[18, 19]. The GO identifiers for each target were determined on the basis of the BLASTX score and the GO evidence code (Table 1, for detailed results, see Additional file 2).
As mRNA was stored in conidia, the expression levels of genes in cluster I declined from the onset of germination. Function annotation showed that about 70 genes out of 1561 were predicted with function related to "response to external or abiotic stimulus" and 30 of those fall into cluster I. These results were consistent with the expectation that genes in cluster I may be involved in perceiving the changes of environmental conditions and play important roles in triggering conidial germination. Cluster I had 15 of 29 genes assigned the function "transcription" and the largest group assigned the function "regulation of biological process" and "cell communication" containing about half of the two categories of genes. The results implied that the transition from the dormant state to vegetative growth was a very complicated process, and needed various regulation pathways to respond to the changes of environmental conditions and initiate transcription to induce germination.
There are 15 genes devoted to "development" and 12 of them are in cluster I. Similar results were reported Qikai Xu for spore germination in Dictyostelium discoidium, who suggested that a major event in germination is the erasure of the developmental state . We suppose that those genes may be involved in maintaining the dormant state. Expression levels of a sub-cluster of genes in cluster I increased after a decrease from 0 to 4 hr. These genes may be stored in the conidia and required for later growth during germination.
When germination is induced, the conidia initiate isotropic growth. Morphologic investigation showed that the isotropic growth phase continued from 0 hr to 9 hr. In this phase, conidia became brighter and the diameter was larger than that in the dormant phase. Genes in cluster II, IV and V were induced and increased their expression levels immediately following the onset of germination. The difference between genes in these three clusters was that the expression level of cluster V genes was highest at 3–4 hr and decreased thereafter, whereas those in clusters II and IV increased continuously. This result implies that genes in all three clusters are involved in isotropic growth but only those in clusters II and IV were needed for germ tube emergence and vegetative growth.
Function annotation indicated that 401 out of all 1561 genes function in "metabolism". The largest group of 137 of those genes belonged to cluster V, 17 belonged to cluster II and 95 belonged to cluster IV. Apart from 117 of the genes in cluster I and 35 genes in cluster III, the expression levels of most metabolism-related genes increased immediately at the onset of germination and approached the highest expression levels.
About 149 of all 1561 genes functioned in "protein synthesis" and 88 of those are in cluster V, 5 are in cluster II and 17 are in cluster IV. This is consistent with the results that had been obtained for A. nidulans: cloning of genes that encode ribosomal proteins (rp genes) was shown to be induced rapidly upon activation of spore germination, and their mRNA levels were lower than those in the subsequent polarized growth . In N. crassa, previous work had showed that a group of genes that function in protein synthesis and as ribosome components were over-represented and had maximum expression levels in the first 4 hr after induction of germination . A total of 173 of all 1561 genes predicted in "biosynthesis" have the expression profiling similar to that in "protein synthesis".
Genes in cluster III were induced at 6 hr and may not be required in the early stage of germination but are involved in polarity growth and germ tube formation. Most genes involved in biological process in cluster III are devoted to "metabolism". Gene expression profiling analysis showed that genes in cluster II and IV also have high expression levels and may be involved in later phase of germination. In polarity growth phase, a large group of genes function in "protein synthesis" belongs to cluster V, which expression levels were fallen from 4 hr during germination. However more than half of genes related to "catabolism" belonged to cluster III and IV, their expression levels were increased as time extended. These results suggest that catabolism may become more active at polarity growth stage.
Comparison of signalling transduction pathways and cellular events involved in conidial germination in T. rubrum and yeast
To better understand the cellular processes and signalling pathways involved in conidial germination, the EST sequences of all 1561 genes were matched to GO terms at the Saccharomyces Genome Database (for detailed results, see Additional file 2). A total of 37 homologues of genes involved in several signalling pathways or regulation modules, such as the cAMP/PKA signalling pathway, the Ras-GTPase-related signalling pathway, the MAPK signalling pathways and two-component signal transduction systems were found and their expression levels were induced during the process (Table 2, more detailed results are also provided in Additional file 3). These results suggest that these signalling modulus may be involved in regulation of germination of T. rubrum.
cAMP/PKA signalling pathways and RAS-like GTPase signalling have been characterized, though are not always equivalent in function, playing important roles to perceive the external signal and initiate the germination in S. cerevisae, A. nidulans and other fungi [23, 24]. To T. rubrum, a sch9 homologue, DW685958, was found, and it belongs to cluster I. sch9 is sch9A in A. nidulans and seems to be involved in cAMP-independent activation of a PKA target in response to nitrogen or amino acids when a fermentable carbon source is abundant [23–27], and is important in spore germination in a.nidulans and in S. cerevisiae. Three Ras2 homologues, DW680378, DW683586, and DW679797, were identified during conidial germination of T. rubrum. DW683586/ras2 is in cluster I, and DW680378/ras2 and DW679797/ras2 are fallen in cluster V. These results suggest that Ras-related GTPase and cAMP/PKA signalling pathways may also play roles in of conidia germination in T. rubrum.
Three homologues of S. cerevisae genes sln1, ypd1 and skn7 related to two-component signal transduction systems were found, and DW695308 and DW705117 belong to cluster I. Sln1p regulates the HOG1 mitogen-activated protein kinase (MAPK) in response to changes in extracellular osmolarity in S. cerevisae . Sln1p also mediates a phosphorelay to a second response regulator, Skn7p, in the cell . Since two-component signal transduction systems are very conservative, they may have similar effects in conidial germination of T. rubrum.
In order to characterize the regulation systems involved in different phases during germination, the homologue genes of S. cerevisae found in all 1561 genes were also compared with the genes in "cell wall" and "polarity" related GO terms at Saccharomyces Genome Database (Table 3 and 4, more detailed results are also provided in Additional files 4 and 56 respectively).
There are 20 homologues of S. cerevisae genes related to "cell wall" induced during conidial germination and 15 of them fall in clusters II, IV and V (Table 3). These genes had a high level of expression in the isotropic growth phase. These results are consistent with the expectation that "cell wall"-related cellular processes were induced to avoid the increased volume resulting in cell lysis/bursting in isotropic growth . A total of 10 out of those 20 homologue genes were characterized as being related to the Rho-pkc-MAPK pathway, the MAPK signalling pathways or the septation machinery (related to the P-loop GTPases), or their interacting proteins, which indicates that these signalling pathways may be involved in the regulation of isotropic growth.
The isotropic growth phase finally leads to polarity growth and results in the formation of a germ cell in both yeast and filamentous fungi . There are 15 homologues of S. cerevisae genes devoted to the development of cell polarity found in the germination process and 14 of those belonging to cluster II, III and IV, respectively, having maximal expression levels in the polarity growth stage (Table 4). Out of those genes, 9 of their homologue genes in S. cerevisae are related to the Rho-GTPase, Ras-GTPase, MAPK and cAMP/PKA signalling pathways, and the septation machinery or their interacting proteins; 8 of them are in cluster IV or cluster III and have increased expression levels in polarity growth. Most of these genes are related to Rho-type GTPase modules, septation machinery and the MAPK signalling pathways. These results suggest that the Rho-type GTPase modules, the septation machinery and MAPK signalling pathways may play roles in regulation of cell polarity development during conidial germination in T. rubrum.
Discussion and Conclusion
As in other filamentous fungi, the conidial germination process in T. rubrum can be divided into three steps: (1) an activation step triggered by environmental factors; (2) an isotropic growth phase representing the first morphological event, referred to as swelling; and (3) a polarized growth phase. In T. rubrum, conidial germination is accompanied by two morphological transitions, swelling and germ tube emergence. The swelling appears at 3–4 hr and continues until germ tube emergence at 9–10 hr after induction of germination. In order to monitor the correlation between the pattern of gene expression with the biological and physiological events that occur during germination, a cDNA microarray containing 10250 ESTs was developed and the gene expression profiling was evaluated in a 15 hr time-series experiment throughout the process of germination.
Traditional protocols for microarray data normalization use a 'control' RNA sample from a particular tissue or time point (RNA normalization), a pool of 'grouped' RNA samples from different tissues or different time points. There are several problems with these approaches. For example, at some growth conditions or time point, the transcription levels of some genes will be undetectable (or very low). Furthermore, for microbial systems, the 'grouped RNA normalization' procedure may require pooling RNA from 20 or 30 experimental conditions at different growth phases. In response to these problems, genomic DNA standards for gene expression profiling had been tested and proposed by Talaat et al . In this report, we used genomic DNA as the reference control to normalize the relative expression levels under the mRNA hybridization signals. The genomic normalization procedure provided reproducible hybridization signals for 96% of the 10250 ESTs. Our results suggest that genomic DNA standards can also expediently be used to evaluate the gene expression profiling in T. rubrum.
A total of 1561 genes with strong statistical estimates on expression levels were obtained and clustered. The results show that the alteration of gene expression levels does not correspond to morphological transitions at certain time-points. There are three turning points, at 0 hr, 2–3 hr and 6–7 hr, on the curves of gene expression profiling during germination. Gene clustering and function analysis suggested that isotropic growth is induced immediately after the initiation of germination, and polarity growth may start at 6 hr or earlier during germination. These results imply that changes of genes expression profiling precede the morphological transitions.
The ability of fungal spores to store pre-packaged mRNA has been revealed in S. cerevisae, A. nidulans and N. crassa [32, 33]. These stored mRNA are primed for rapid activation and translation in the presence of nutrients. The decay of spores containing mRNA is induced soon after germination . In this study, mRNA of genes in cluster I existed in spores before induction of germination and the expression levels of those genes were down-regulated immediately at the onset of germination. Functional annotation showed that more than half of the genes related to "transcription", "cell communication" and "regulation of biological process" are in cluster I, and there are many genes in this cluster that function in "response to external or abiotic stimulus". These results imply that stored mRNA is important for the activation of germination and the transition of spores from the dormant state to vegetative growth is complicated, in that it is induced and regulated by both exterior factors and intra-cell signalling.
Spores are not simply quiescent cells, a basal level of RNA and protein synthesis is required for spore survival. So, we supposed that cluster I may include two types of genes. One type may be stored in conidia and induce decay after translation during germination, whereas the second type may have specific expression in conidia that decayed immediately after induction of germination. Most genes that function in "development" are in cluster I and are down-regulated at the onset of germination. Similar results have been reported by Xu et al . We suppose that genes may be involved in maintaining the dormant state. Since we have little knowledge about the events that occur in dormant conidia and the activation of germination in T. rubrum, the identity of the spore-contained mRNA may provide some clues for further research on the molecular mechanism of conidial survival and germination in T. rubrum.
The expression levels of genes in clusters II, IV and V that increased rapidly after induction of germination may be involved in isotropic growth. Gene functions in "metabolism" also are immediately induced and have high expression levels in isotropic growth indicate that dormant conidia breaking dormant state is a very quick process.
Protein synthesis has been shown to be required for germination [35–37]. In our data, most genes function in "protein biosynthesis" belongs to those clusters and achieve the maximum expression levels in isotropic growth. These results are consistent with previous reports that conidia contain fully active ribosomes and it is therefore likely that the activation of the rp genes during conidial germination, and hence de novo synthesis of ribosomes, is a prerequisite to achieve a growth rate appropriate for germ tube emergence that could not be obtained by the sole use of the translation machinery stored in the conidia . With the passage of time, transcription of the genes related to "catabolism" was induced and became more active, and the conidia gradually enter the vegetative growth state.
The budding yeast S. cerevisiae, with the powerful tools applicable in this organism, offers the opportunity to rapidly characterize signal transduction pathways of a eukaryotic cell in great detail. Therefore, yeast may be used as a reference library-much as its genome sequence has been-for the analysis of conserved signaling pathways, e.g., in the more complex, multicellular filamentous fungi . To better understand the mechanism and cellular events be involved in T. rubrum conidial germination, we matched The ESTs of the selected 1561 genes to GO terms at the Saccharomyces Genome Database. Some homologues to S. cerevisae genes related to signalling modules and important cellular processes were found and distributed in various clusters. Our results indicated that though the molecular mechanism and biological processes are differ in morphological construction, there are several conservative signalling pathways and similar cellular events between budding yeast and T. rubrum germination.
The primary requirement for initiation of germination and completion of the subsequent steps is the sensing of external signals. In S. cerevisae and A. nidulans, glucose or other fermentable carbon sources as external signals are necessary and sufficient to trigger spore germination. In S. cerevisiae, glucose sensing is mediated by the G-protein-coupled receptor (GPCR) Gpr1p that in turn activates the heterotrimeric G-protein α-subunit encoded by the GPA2 gene [39, 40]. The Gpr1p-Gpa2p system mediates glucose-dependent activation of the cAMP-dependent protein kinase (PKA) pathway that is associated with mobilization of trehalose, decreased stress resistance, and expression of ribosomal protein (rp) genes . The Ras pathway is the rate-limiting step of spore germination, since an elevated level of activated Ras protein increased the rate of germination [37, 41]. In A. nidulans, a very recent report has revealed GanB regulates conidial germination within the heterotrimeric G-protein GanB(α)-SfaD(β)-GpgA(γ) through activation of the cAMP/PKA pathway in response to glucose. A previous study revealed that RasA from A. nidulans regulates conidial germination via an undefined signalling pathway in parallel to the cAMP/PKA pathway. cAMP/PKA signalling controls early events of conidial germination in response to sensing a carbon source, and plays a critical, but not essential, role in the germination process. Indeed, inactivation of adenylate cyclase results in a severe delay in, but not a complete arrest of, germ tube emergence. Ras signalling may be a response to carbon source signalling initiating germination, but plays more important roles in inducing conidial germination . In our data, genes in cluster I are present as mRNA in the conidia and may respond to external signals triggering the germination. Several homologues of genes related to the cAMP/PKA and RAS GTPase signalling pathways were found related to conidial germination in T. rubrum. However, we also found that glucose only could not trigger the conidia germination and conidia could germinate and grow in nitrogen abundant but fermentable C source absent media such as soy protein liquid medium (SP) and keratin liquid medium (KSP) [42–44]. So, we suppose that the signal transduction modules are different molecular mechanism involved in sensing the external signal and initiation of germination in T. rubrum.
Isotropic growth results in a drastically increased volume. In our data, the "cell wall"-related process genes were significantly induced and most of them had the maximum expression levels in the isotropic growth phase. These results indicate that cell wall biosynthesis and integrity maintenance are also very important for isotropic growth of T. rubrum. Cell wall biosynthesis and monitoring of cell wall integrity is important to avoid cell lysis/bursting and may account for the increase of gene expression levels related to "cell wall". Two genes related to the MAPK and the Ras/Rho-type GTPases signalling pathways were found and shown to have similar expression profiling. Since Rho-GTPases and the MAPK signalling pathways had been characterized as being involved in the regulation of cell-wall biosynthesis in S. cerevisae and some other filamentous fungi, these signal transduction modules may play roles in cell-wall biosynthesis and be involved in the regulation of isotropic growth.
Polarity growth is also essential for germination and hyphal growth. In Ashbya gossypii, Agbem2 mutant loss of cell polarity eventually caused completely isotropic growth, resulting in large, balloon-shaped tip cells . Comparison with the GO terms at the Saccharomyces Genome Database indicated that there are 23 homologues of S. cerevisae genes devoted to "cell polarity development" found be involved in the germination process and half of those share high expression levels at polarity growth phase. Our data also suppose that Rho-type GTPase modules, septation machinery and MAPK signalling pathways may take part in regulation of cell polarity development and play roles in conidial germination of T. rubrum.
Two-component signal transduction pathways were found that may be involved in the conidial germination of T. rubrum. Two-component signal transduction pathways are used extensively to mediate prokaryotic signalling events. In recent years, these signalling systems have been found in eukaryotes, including plants, yeasts, filamentous fungi and slime moulds [45, 46]. In C. albicans, these signal transduction pathways regulate cell-wall biosynthesis (and, therefore, adherence to host cells), osmotic and oxidant adaptation, white-opaque switching, morphogenesis, and virulence of the organism . Since the presence of two-component signal transduction pathways has not been demonstrated in mammals, further study on these signalling pathways may lead to the discovery of new drug targets for T. rubrum control.
Strain and culture conditions
T. rubrum strain BMU01672 was isolated from nail scrapes of a patient suffering from tinea unguium. The strain was confirmed as T. rubrum by morphologic identification, as well as by PCR amplification and sequencing of the 18 S ribosomal DNA and internal transcribed spacer (ITS) regions. Strain reference samples are stored at the Research Centre for Medical Mycology, Beijing, China.
T. rubrum was grown on potato glucose agar (Difco) at 28°C for 10 days to produce conidia . The conidia were washed from the medium at 4°C with distilled water and passed through a 70 μm pore size nylon filter twice to remove hyphal fragments. The conidia sample was centrifuged at 1800 g for 10 min, and then resuspended to give a final concentration of 108 conidia/mL.
For germination analysis, 20 mL of conidia suspension was introduced into each of a series of 250-mL flasks containing 80 mL of Sabouraud liquid medium (containing 49 g of glucose, and 10 g of Difco Bacto peptone in 1 L of distilled water), and immediately incubated at 28°C with constant shaking (200 rpm, Innova 4230 Refrigerated Incubator Shaker, New Brunswick Scientific, Edison NJ). Two samples at time-point 0 hr were taken directly from the stored conidia. Subsequent samples were harvested every hour from the first to the fifteenth hour during the germination process. Two samples were collected from the flasks as independent biological duplicates, so a total of 32 samples were obtained during the whole process of germination from 0 hr to 15 hr
Morphological transitions during germination were examined using an OLYMPUS CKX41SF microscope (OLYMPUS Optical CO., LTD), and viewed under a 40 × /0.55 PH2 objective lens. Micrographs were taken with an OLYMPUS digital camera and captured images were processed using Adobe Photoshop software.
cDNA microarray construction
PCR fragments used for printing the microarray chip were amplified from the EST library in 96-well plates using vector-PCR amplification with T7 and SP6 universal primers on the GeneAmp PCR system 9600 (Perkin Elmer, Foster City, CA). A single EST cDNA clone was selected to represent each assembled sequence (i.e. putatively unique transcript). For an assembled contig represented by multiple ESTs, the rule followed was to select the clone with the most extensive DNA sequencing read length. PCRs were initiated with DNA denaturation (5 min at 96°C) followed by annealing for 30 s at 55°C and extension (1 min and 30 s at 72°C) steps carried out for 35 cycles. PCR products were puried using the QIAquick PCR Puri®cation kit (Qiagen) and eluted in 100 μl sterile water. Purified PCR products were dried and then resuspended with a concentration of 100 ng/μL.
A set of microarrays containing a total of 11,232 spots, including 10, 250 clones in the form of PCR products and 982 controls in each block, including blank controls, negative controls, and positive controls, were spotted onto CMT-GAPSII-coated slides (Corning Glass) using a 16-pin configured MicroGrid II array printer (BioRobotics) controlled by the MicroGrid TAS Application Suite (version 126.96.36.199).
The spotted cDNA was cross-linked to the surface of the slides (at 65 mJ) by using a StrataLinker instrument and washed with a 1% SDS solution to minimize the background. Slides were subsequently placed in a blocking solution containing 200 mM succinic anhydride and 50 mM sodium borate prepared in 1-methyl-2-pyrrolidinone for 20 min, washed for 2 min in 95 °C distilled water, and rinsed five times in 95% ethanol. Slides were spin dried at 500 rpm for 5 min and stored for future hybridizations.
RNA extraction and cDNA labelling
Samples were centrifuged at 1800 g for 10 min at 4°C and immediately grounded into mortars with liquid nitrogen. Total RNA were isolated and purified following the protocols of the RNeasy® Plant Mini Kit (Qiagen, Valencia, CA). The RNA concentration and purity were determined spectrophotometrically by measuring absorbance at 230, 260, 280, and 320 nm. The purity and integrity of the RNA samples were confirmed by agarose gel electrophoresis. Poly (A)+ mRNA was isolated with Oligotex mRNA Mini Kit (Qiagen).
Test mRNA enriched from 50 μg of total RNA was reverse transcribed into cDNA and labeled with a fluorescent dye (Cy3 or Cy5). Also, 5 μg of reference genome DNA isolated in the hyphal stage was labelled. Test samples and reference samples were dye-swapped to avoid variation arising from the dyes.
cDNA microarray hybridization
Labelled cDNA and reference genomic DNA were purified using the QIAquick PCR Puri®cation kit (Qiagen), then mixed and resuspended in 10 μL of distilled water, mixed with 1.5 μL of 50 × Denhardts, 2.25 μLof 20 × SSC, 1.125 μL of 500 μg/mL yeast tRNA, 0.375 μL of 1 M HEPES (pH 7.0), 0.375 μL of 10% SDS, 2 μL of poly (A) (5 mg/μL). The mixture was heated 100°C for 3 min. Hybridization was carried out as described by Hayes et al .
Image acquisition and data analysis
The processed slides were scanned with a GenePix 4100B scanner (Axon). Fluorescent spots and the local background intensities were quantified with GenePix Pro 6.0 software (Axon). The local background value was subtracted from the intensity of each spot. The mean of the signal intensities of the control spots hybridized with labeled reference genomic DNA in each experiment was calculated. Ten different A. thaliana genes and human β- actin gene from SpotReport™ cDNA Array Validation System (Stratagene) were used for the controls. The spots that showed intensity with labeled reference genomic DNA that was lower than the mean value of the control spots were excluded from further analysis. The mean log2 (sample/reference) ratios of signal intensity were calculated for analysis .
The raw data consisted of 64 hybridizations (32 dye-swap pairs) on duplicate printed slides from the same batch. At every time-point, two biological replicates containing 4 hybridized data were used for further analysis. These multiple data sets were flagged, fitting all the features: spot Dia. > = 80 μm, %B (532 or 635) + 2SD > 55, SNR635 (or 532) > = 3 and normalized (the ratio of medians of all features equal to 1) by the GenePix Pro software (version 6.0). The data sets were further normalized in two steps: total intensity and Lowess normalization using Tiger MIDAS V2.19 .
After spots flag and normalization process, 9470 spots (about 93.7% of all) in a microarray were used for the gene expression profiling analysis. The expression variations during the time-course were analyzed by ANOVA (p < 0.01) for each gene, using Tiger TMEV 3.1 software . 2772 genes were obtained that their expression levels were significantly altered during the germination. In order to identify genes whose expression levels were altered dramatically during germination, 1576 out of the 2772 genes fitting expression levels altered exceeding twofold were used for further analysis. The estimated expression levels of these genes were clustered by Tiger TMEV 3.1 software using a hierarchical clustering method, in which similarity in expression patterns between genes is measured as Pearson's correlation coefficient, and the closest two genes or clusters are successively joined. Distances between clusters represent the average distances between genes in the clusters. Five representative clusters were chosen visually, each of which showed distinctive profiles (Figure 2). The genes enrichment within cluster and the statistical significance of gene distribution across the clusters were estimated by post-hoc tests(p < 0.05) using SPSS 10.0 For Windows (detailed results see Additional file 1) .
To annotate the putative functions of selected genes and predicate the biological process involved in germination, we used a controlled vocabulary for describing gene function [18, 19]. The ESTs of 1576 selected genes were matched to the sequences in the GO database. GO identifiers were determined on the basis of the BLASTX score and the GO evidence code of the homologous genes (see Additional file 2 for details).
Quantitative real-time RT-PCR
In order to verify the microarray result, the relative expression levels of 8 genes at several time points (0, 4, 10, and 15) were estimated by Quantitative real-time RT-PCR. First-strand cDNAs were synthesized from 2 μg of total RNA in a 100-μl reaction volume using the SuperScript First-Strand Synthesis System for RT-PCR (Invitrogen, Carlsbad, CA) in accordance with the manufacturer's instructions. Quantitative real-time PCR experiments were performed in triplicate using the 7000 Sequence Detection System (Applied Biosystems, Foster City, CA). Independent PCRs were performed using the same cDNA for both the selected genes and the 18S rRNA, using the SYBR Green PCR Master Mix (Applied Biosystems). Gene-specific primers were designed for the genes and 18S rRNA using Primer Express software (Applied Biosystems) and are shown in Additional file 1. The PCR cycle consisted of AmpliTaq Gold activation at 95°C for 10 min, followed by 40 cycles of denaturation at 95°C for 15 s and annealing/extension at 58°C for 1 min. A dissociation curve was generated at the end of each PCR cycle to verify that a single product was amplified using software provided with the 7000 Sequence Detection System. The changes in fluorescence of SYBR Green I dye in each cycle were monitored by the system software, and the calculated threshold cycle (Ct) for each gene amplification was normalized to Ct of the 18S rRNA gene amplified from the corresponding sample before calculating the fold change from a selected time point to 0 time point using the following formula:
fold change = 2ΔΔCt
where ΔΔCt for gene j = (Ct,j- Ct,18S rRNA)a time point - (Ct,j- Ct,18S rRNA)time point 0.
EST sequences used for cDNA microarray preparation were deposited in GenBank under accession numbers: [DW405580–DW407270 and DW678211–DW711189]. The microarray related data were submitted to Gene Expression Omnibus (GEO) under accession number: [GSE5083]
Gpaser Y Kuijpers AFA, Prersber W, de HOOG GS: Molecular Taxonomy of the Trichophyton rubrum Complex. J Clin Microbiol. 2000, 38: 3329–3336-
Costa M Passos, XS, Souza LKH, Miranda ATB, Lemos JÁ, Oliveira Júnior JG, Silva MRR: Epidemiology and etiology of dermatophytosis in Goiania, GO, Brazil. Rev Soc Bras Med Trop. 2002, 35: 19-22. 10.1590/S0037-86822002000100004.
Jennings MB Weinberg JM, Koestenblatt EK, Lesczczynski C: Study of clinically suspected onychomycosis in a podiatric population. J Am Podiatr Med Assoc. 2002, 92: 327-330.
Monod M Jaccoud S, Zaugg C, Lechenn B, Baudraz F, Panizzon R: Survey of dermatophyte infections in the Lauseanne area (Switzerland). Dermatology. 2002, 205: 201-203. 10.1159/000063913.
Sentamilselvi C Janaki C, Kamalam A, Thambiah AS: Deep dermatophytosis caused by Trichophyton rubrum: a case report. Mycopathologia. 1998, 142: 9-11. 10.1023/A:1006984116046.
Smith KJ Welsh M, Skelton H: Trichophyton rubrum showing deep dermal invasion directly from the epidermis in immunosuppressed patients. Br J Dermatol. 2001, 145: 344-348. 10.1046/j.1365-2133.2001.04331.x.
Squeo RF, Beer R, Silvers D, Weitzman I, Grossman M: Invasive Trichophyton rubrum resembling blastomycosis infection in the immunocompromised host. J Am Acad Dermatol. 1998, 39: 379-380. 10.1016/S0190-9622(98)70396-2.
Gupta AK Beer R, Silvers D, Weitzman I, Grossman M: Arthroconidial formation in Trichophyton raubitschekii. Mycoses. 2003, 46: 332-338.
Kane J Summerbell RC, Sigler L, Krajden S, Land G: Laboratory Handbook of Dermatophytes: A Clinical Guide and Laboratory Manual of Dermatophytes and Other Filamentous Fungi from Skin, Hair, and Nails. 1997, Belmont, Calif: Star Publishing Company.
Zurita J Hay RJ: Adherence of dermatophyte microconidia and arthroconidia to human keratinocytes in vitro. J Invest Dermatol. 1987, 89: 529–534.-10.1111/1523-1747.ep12461067.
Giaever G Chu AM, Ni L, Connelly C, Riles L, Veronneau S, Dow S, Lucau-Danila A, Anderson K, Andre B, et al.: Functional profiling of the Saccharomyces cerevisiae genome. Nature. 2002, 418: 387-391. 10.1038/nature00935.
Hughes TR Marton MJ, Jones AR, Roberts CJ, Stoughton R, Armour CD, Bennett HA, Coffey E, Dai H, He YD, et al.: Functional discovery via a compendium of expression profiles. Cell. 2000, 102: 119–126.-10.1016/S0092-8674(00)00015-5.
Cervelatti EP Ferreira-Nozawa MS, Aquino-Ferreira R, Fachin AL, Martinez-Rossi NM: Electrophoretic molecular karyotype of the dermatophyte Trichophyton rubrum. Genet Mol Biol. 2004, 271: 99-102.
Wang LL Ma L, Leng WC, Yang J, Zhu JP, Dong J, Xue Y, Zhe W, Li RY, Jin Q: Analysis of part of the Trichophyton rubrum ESTs. Sci in Chi Ser C. 2004, 47: 389—395-
Wang LL Ma L, Leng WC, Liu T, Yu L,Yang J, Yang L, Zhang WL, Zhang Q, Dong J, Xue Y, Zhu YF, Xu XY, Wan Z, Ding GH, Yu FD, Tu K, Li YX, Li RY, Shen Y, Jin Q: Analysis of the dermatophyte Trichophyton rubrum expressed sequence tags. BMC Genomics. 2006, 7: 225-10.1186/1471-2164-7-225.
Denning DW Anderson MJ, Turner G, Latge JP, Bennet JW: Sequencing the Aspergillus fumigatus genome. Lancet Infect Dis. 2002, 2: 251-253. 10.1016/S1473-3099(02)00243-8.
Saeed AI Sharov V, White J, Li J, Liang W, Bhagabat N, Braisted J, Klapa M, Currier T, Thiagarajan M, et al: TM4: a free, open-source system for microarray data management and analysis. Biotechniques. 2003, 342: 374-378.
Ashburner M Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al.: Gene Ontology: tool for the unification of biology. Nat Genet. 2000, 25: 25–29.-
Kellam P: Microarray gene expression database: progress towards an international repository of gene expression data. Genome Biology. 2001, 2: reports4011.1-reports4011.3. 10.1186/gb-2001-2-5-reports4011.
Xu Q Ibarra M, Mahadeo D, Shaw C, Huang E, Kuspa A, Cotter D, Shaulsky G: Transcriptional Transitions during Dictyostelium Spore Germination. Eukaryot Cell. 2004, 3: 1101-1110. 10.1128/EC.3.5.1101-1110.2004.
D’Enfert C: Fungal Spore Germination: Insights from the Molecular Genetics of Aspergillus nidulans and Neurospora crassa. Fungal Genet & Biol. 1997, 21: 163–172-
Kasuga T Townsend JP, Tian CG, Gilbert LB, Mannhaupt G, Taylor J W, Glass NL: Long-oligomer microarray profiling in Neurospora crassa reveals the transcriptional program underlying biochemical and physiological events of conidial germination. Nucleic Acids Res. 2005, 33: 6469–6485.-10.1093/nar/gki953.
Fillinger S Chaveroche MK., Shimizu K, Keller N, d’Enfert C: cAMP and ras signalling independently control spore germination in the filamentous fungus Aspergillus nidulans. Mol Microbiol. 2002, 44: 1001-1016. 10.1046/j.1365-2958.2002.02933.x.
Yu J Corey D, James RB: Efficient transition to growth on fermentable carbon sources in Saccharomyces cerevisiae requires signaling through the Ras pathway. EMBO J. 1998, 17: 6942-6951. 10.1093/emboj/17.23.6942.
Crauwels M Donaton MC, Pernambuco MB, Winderickx J, de Winde JH, Thevelein JM: The Sch9 protein kinase in the yeast Saccharomyces cerevisiae controls cAPK activity and is required for nitrogen activation of the fermentable-growth-medium-induced (FGM) pathway. Microbiolgy. 1997, 143: 2627–2637.-
Lorenz MC Pan X, Harashima T, Cardenas ME, Xue Y, Hirsch JP, Heitman J: The G proteincoupled receptor Gpr1 is a nutrient sensor that regulates pseudohyphal differentiation in Saccharomyces cerevisiae. Genetics. 2000, 154: 609–622.-
Fabrizio P Pozza F, Pletcher SD, Gendron CM, Longo VD: Regulation of longevity and stress resistance by Sch9 in Yeast. Science. 2001, 292: 288–291.-10.1126/science.1059497.
Posas F Wurgler-Murphy SM, Maeda T, Witten EA, Thai TC, Saito H: Yeast HOG1 MAP kinase cascade is regulated by a multistep phosphorelay mechanism in the SLN1-YPD1-SSK1 "two-component" osmosensor. Cell. 1996, 86: 865-875. 10.1016/S0092-8674(00)80162-2.
Li S AA Malone CL, Raitt D, Dean S, Johnston LH, Deschenes RJ, Fassler JS: The yeast histidine protein kinase, Sln1p, mediates phosphotransfer to two response regulators, Ssk1p and Skn7p. EMBO J. 1998, 17: 6952-6962. 10.1093/emboj/17.23.6952.
Wendland J Philippsen P: Cell polarity and hyphal morphogenesis are controlled by multiple Rho-protein modules in the filamentous ascomycete Ashbya gossypii. Genetics. 2001, 157: 601–610.-
Talaat AM Howard ST, Hale W, Lyons R, Garner H, Johnstonb SA: Genomic DNA standards for gene expression profiling in Mycobacterium tuberculosis. Nucleic Acids Res. 2002, 30: e104. -10.1093/nar/gnf103.
Hollomon DW: RNA synthesis during fungal spore germination. J Gen Microbiol. 1970, 62: 75-87.
Osherov N May SG: The molecular mechanisms of conidial germination. FEMS Microbiol Lett. 2001, 199: 153-160. 10.1111/j.1574-6968.2001.tb10667.x.
Brengues M Pintard L, Lapeyre B: mRNA decay is rapidly induced after spore germination of Saccharomyces cerevisiae. J Biol Chem. 2002, 277: 40505-40512. 10.1074/jbc.M206700200.
Bonnen A Brambl R: Germination physiology of Neurospora crassa conidia. Exp Mycol. 1983, 7: 197–207.-10.1016/0147-5975(83)90040-3.
Brambl R Plesofsky-Vig N, Hammett JR, Russell PL: Preserved polyadenylated ribonucleic acid in dormant conidia of Neurospora crassa and new RNA synthesis during spore germination. Exp Mycol. 1987, 11: 317–330.-10.1016/0147-5975(87)90020-X.
Herman PK Rine J: Yeast spore germination: A requirement for Ras protein activity during re-entry into the cell cycle. EMBO J. 1997, 16: 6171–6181.-10.1093/emboj/16.20.6171.
Wendland J: Comparison of morphogenetic networks of filamentous fungi and yeast. Fungal Genet & Biol. 2001, 34: 63 -682. 10.1006/fgbi.2001.1290.
Kraakman L Lemaire K, Ma PS, Teunissen A, Donaton MCV, Van Dijck P, Winderickx J, De Winde J, Thevelein J: A Saccharomyces cerevisiae G-protein coupled receptor, Gpr1, is specifically required for glucose activation of the cAMP pathway during the transition to growth on glucose. Mol Microbiol. 1999, 32: 1002–1012.-10.1046/j.1365-2958.1999.01413.x.
Xue Y Batlle M, Hirsch J P: GPR1 encodes a putative G protein-coupled receptor that associates with the Gpa2p Galpha subunit and functions in a Ras-independent pathway. EMBO J. 1998, 17: 1996–2007.-10.1093/emboj/17.7.1996.
Jiang Y Davis C, Broach JR: Efficient transition to growth on fermentable carbon sources in Saccharomyces cerevisiae requires signaling through the Ras pathway. EMBO J. 1998, 17: 6942-6951. 10.1093/emboj/17.23.6942.
Brasch J MBS Christophers E: Enzyme release by Trichophyton rubrum depends on nutritional conditions. Mycoses. 1991, 34: 365-368.
Jousson O Lé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.
Michel M Barbara L, Olivier J. Daniela G, Christophe Z, Reto S, Eric G: Aminopeptidases and dipeptidyl-peptidases secreted by the dermatophyte Trichophyton rubrum. Microbiology. 2005, 151: 145–155-
Stock AM Robinson VL, Goudreau PN: Two-component signal transduction. Annu Rev Biochem. 2000, 69: 183–215.-10.1146/annurev.biochem.69.1.183.
Wolanin PM Thomason PA, Stock JB: Histidine protein kinases: key signal transducers outside the animal kingdom. Genome Biol. 2002, 3: 3013.1–3013.8.-10.1186/gb-2002-3-10-reviews3013.
Michael K Richard,C.: Two-component signal transduction in human fungal pathogens. FEMS Yeast Res. 2006, 6: 149–159.-
Jessup CJ Waner J, Isham N, Hasan I, Ghannoum MA: Antifungal susceptibility testing of dermatophytes: establishing a medium for inducing conidial growth and evaluation of susceptibility of clinical isolates. J Clin Microbiol. 2000, 38: 341–344-
Hayes A Zhang N, Wu J, Butler PR, Hauser NC, Hoheisel JD, Lim FL, Sharrocks AD, Oliver SG: Hybridization array technology coupled with chemostat culture: Tools to interrogate gene expression in Saccharomyces cerevisiae. Methods. 2002, 26: 281-290. 10.1016/S1046-2023(02)00032-4.
Peng JP Zhang XB, Yang J, Wang J, Yang J, Bin W, Wei CD, Sun MS, Jin Q: The use of comparative genomic hybridization to characterize genome dynamics and diversity among the serotypes of Shigella. BMC Genomics. 2006, 7: 218-10.1186/1471-2164-7-218.
Pan KH Lih CJ, Cohen SN: Analysis of DNA microarrays using algorithms that employ rule-based expert knowledge. Proc Natl Acad Sci USA. 2002, 99: 2118-2123. 10.1073/pnas.251687398.
The work was supported by The National High Technology Research and Development Program of China (Accession number: 2006AA020504) and National Key Technologies R&D Programme (Accession Number: 2002BA711A14).
TL carried out the T. rubrum functional genomics studies, EST sequence analysis, cDNA microarray preparation, the gene expression analysis and participated in drafting the manuscript. QZ contributed to cDNA microarray preparation and the gene expression analysis. LW, LY, WL, JP, JY, LM contributed to EST sequencing, data analysis and cDNA microarray preparation. LC, JD, XX, YX, WZ, YZ, LY, LQ, DY, GD, KT participated in EST sequence analysis. LS, WL, ZW contributed to T. rubrum material preparation. RL, YxL, YS participated in design of the studs. QJ proposed the research goal, supervised the whole studies and provided a critical review of the manuscript. All authors read and approved the final manuscript.
Tao Liu, Qian Zhang, Lingling Wang contributed equally to this work.