The Scientific World Journal

The Scientific World Journal / 2012 / Article

Research Article | Open Access

Volume 2012 |Article ID 340565 | 7 pages |

Transcriptional Profiling of ESTs from the Biocontrol Fungus Chaetomium cupreum

Academic Editor: S. Liuni
Received24 Oct 2011
Accepted07 Dec 2011
Published01 Feb 2012


Comparative analysis was applied to two cDNA/ESTs libraries (C1 and C2) from Chaetomium cupreum. A total of 5538 ESTs were sequenced and assembled into 2162 unigenes including 585 contigs and 1577 singletons. BlastX analysis enabled the identification of 1211 unigenes with similarities to sequences in the public databases. MFS monosaccharide transporter was found as the gene expressed at the highest level in library C2, but no expression in C1. The majority of unigenes were library specific. Comparative analysis of the ESTs further revealed the difference of C. cupreum in gene expression and metabolic pathways between libraries. Two different sequences similar to the 48-KDa endochitinase and 46-KDa endochitinase were identified in libraries C1 and C2, respectively.

1. Introduction

One of the reasons for environmental disorder is that modern agriculture is an ecologically unbalanced system which has been destroyed by chemical fungicides. Biocontrol is highly interesting alternative method of chemical plant disease control. C. cupreum Ames is an ascomycete fungus with considerable biocontrol potential to plant fungal pathogens, especially several notorious examples belonging to the genera of Pythium, Rhizoctonia, and Pyricularia [1]. In Thailand and China, its biological products have been applied in agricultural disease management [1, 2]; nevertheless, the genetic basis of the defense mechanisms of C. cupreum is not well understood thus inhibiting its application.

ESTs analysis has been proven to be an efficient and valuable tool in obtaining coding gene information, understanding the pathways involved in a given physiological or environmental stimulus [3, 4]. To date, several ESTs studies have been carried out on fungal biocontrol agents, especially to species of Trichoderma. For example, analysis of 8,710 ESTs of T. harzianum CECT 2413 from eight cDNA libraries including those simulating mycoparasitism [5] and of ESTs from four different Trichoderma strains grown under conditions related to biocontrol [6].

In the present study, we performed comparative analysis of two cDNA/ESTs libraries from C. cupreum. An obvious difference in gene expression and metabolic pathways were detected between libraries. This research contributes to elucidating further the mycoparasitisic molecular mechanisms involved in C. cupreum and will help to develop novel strategies in fungal disease management.

2. Materials and Methods

2.1. Fungal Strains and Culture Conditions

The C. cupreum isolate was kindly provided by King Mongkut’s Institute of Technology of Ladkrabang, Thailand. Rhizoctonia solani was stored in our laboratory.

C1 was constructed from mycelia of C. cupreum grown on potato dextrose agar (PDA) culture. Library C2 was constructed with mycelia of Rhizoctonia solani cell wall (RsCW) as the carbon source, this was designed to resemble the phases of interaction with R. solani. For C1 library construction, mycelia of C. cupreum were cultured on potato dextrose (PD) medium for 60 h at 27°C and 150 rpm. The biomass was harvested and stored at −80°C until use. For C2 library construction, mycelia were initially grown in PD medium with shaking at 27°C and 150 rpm for 24 h then transferred to SM medium [2.8 g (NH4)2SO4 l−1; 0.6 Urea g l−1; 4 g KH2PO4 l−1; 0.6 g CaCl2·2H2O l−1; 0.2 g MgSO4 l−1; 0.01 g FeSO4·7H2O l−1; 0.0028 g ZnSO4·H2O l−1; 0.0032 g CoCl2·6H2O l−1; 5 g RsCW l−1] and incubated for a further 36 h at 27°C and 150 rpm.

2.2. Construction of the cDNA Libraries and DNA Sequencing

Total RNA was extracted using Trizol reagent from mycelia of C. cupreum. Polyadenylated RNA was purified using an Oligotex mRNA Kit (Qiagen). The course of cDNA library construction followed the procedures of Zhang as described before [7]. Unidirectional cDNA libraries were constructed using the pBluescript II plasmid system. Fragments of cDNA clones were sequenced using a T3 primer from the 5′ end with MegaBase1000 DNA sequencer.

2.3. Data Processing and Bioinformatics Analysis

Vector sequences, sequences shorter than 100 bp and containing more than 5% ambiguous bases, were discarded using the Crossmatch program. High-quality sequences were assembled using Phrap ( and accuracy of contigs was confirmed with Consed [8]. All unigenes were compared against public nonredundant (nr) protein databases using a BlastX search. According to KEGG (Kyoto Encyclopedia of Genes and Genomes) [9], unigenes were assigned to different metabolic pathways with the same criterion as described by Zhang [7]. All high-quality ESTs were submitted to the GenBank database under accession numbers DV544375-DV548659.

3. Results

3.1. ESTs Clustering and Function Assignment

A total of 5,538 cDNA clones with an insert size of more than 700 bp were selected for sequencing, resulting in 4285 ESTs (3066 from C1 and 1219 from C2) after removing sequences representing ribosomal, vector, and low-quality sequences. Minimum, average, and maximum lengths of ESTs were 102, 518, and 795 bp, respectively, with a large fraction falling between 500 and 700 bp (2110 from C1 and 795 from C2) in both libraries.

Using the Phrap and Consed programs, ESTs from both libraries were arranged into 585 contigs and 1577 singletons, giving a total of 2162 unigenes. Each unigene was subjected to analysis against homologous sequences in public protein databases using the BlastX algorithm. Approximately 1211 (56%) of the unigenes were assigned a function with an E-value of 10−5 or lower. The remaining 951 clones had no high homology to genes of known function. A total of 1138 (52.6%) and 691 (32%) unigenes were unique and only expressed in C1 or C2, respectively.

3.2. Exploration of Highly Expressed ESTs

Contigs containing 4 or more ESTs from each library are listed in Table 1. Of the 26 clusters, more than one third were (10/26) expressed only in the C1 library, 2 only in the C2 library, and half (14/29) in both but at a different level.

Contig no.AnnotationC1C2

Contig525Glyceraldehyde-3-phosphate dehydrogenase1097.41141.69
Contig531Coproporphyrinogen oxidase724.8950.6
Contig30Predicted protein563.8100
Contig433C-4 sterol methyl oxidase543.6700
Contig313Predicted protein362.4500
Contig2Pyruvate decarboxylase271.8400
Contig39Predicted protein231.5600
Contig507Xylulose-5-phosphate phosphoketolase211.4330.36
Contig 152Predicted protein201.36161.93
Contig523EF1-alpha translation elongation factor191.29121.45
Contig42Thiazole biosynthetic enzyme191.2920.24
Contig495Glutamine synthetase171.1610.12
Contig517Aspartic protease161.0930.36
Contig317Ammonium transporter161.0900
Contig419ATP synthase protein161.0900
Contig4Alcohol dehydrogenase151.0240.48
Contig494ATP citrate lyase140.9510.12
Contig412Histone H2B130.8800
Contig298Pyruvate kinase130.8800
Contig310Aspartate aminotransferase110.7510.12
Contig408Heat shock protein 30110.7500
Contig489ADP-ATP translocase110.7580.97
Contig165MFS monosaccharide transporter00202.42
Contig162Phosphoenolpyruvate carboxykinase00121.45

Relative values (%): ESTs numbers/all ESTs numbers in C1 or C2 library.

Glyceraldehyde-3-phosphate dehydrogenase was the most highly expressed transcript (109 ESTs) in the C1 library, occurring four times more than in the C2 library. The most highly represented transcripts in the C2 library coded for a putatively major facilitator superfamily (MFS) monosaccharide transporter; no such expression was observed in library C1. The expression of coproporphyrinogen oxidase, thiazole biosynthetic enzyme, glutamine synthetase, ATP citrate lyase, and aspartate aminotransferase were higher in the C1 library. It should be noted that many hits similar to predicted or unknown function proteins were detected in libraries. They are ideal candidates for future study.

3.3. Metabolic Pathways Analysis

Pathways analysis of KEGG was carried out using genes homologous to known functional sequences. A total of 65 and 61 different metabolic pathways were found in the C1 and C2 libraries, respectively. These results show evident difference in metabolic pathways between the libraries (Table 2).

Map no.PathwaysESTs in C1ESTs in C2

00020Citrate cycle479
00030Pentose phosphate cycle356
00040Pentose and glucuronate interconversions64
00051Fructose and mannose metabolism410
00052Galactose metabolism77
00053Ascorbate and aldarate metabolism55
00061Fatty acid biosynthesis (path 1)31
00062Fatty acid biosynthesis (path 2)01
00071Fatty acid metabolism2412
00100Sterol, vitamin K, vitamin E, carotenoids biosynthesis10
00120Bile acid biosynthesis1912
00130Ubiquinone biosynthesis143
00150Androgen and estrogen metabolism412
00190Electron transport and oxidative phosphorylation346
00220Urea cycle and metabolism of amino groups100
00230Purine metabolism411
00240Pyrimidine metabolism61
00251Glutamate metabolism244
00252Alanine and aspartate metabolism141
00260Glycine, serine, and threonine metabolism28
00271Methionine metabolism125
00272Cysteine metabolism61
00280Valine, leucine, and isoleucine degradation66
00290Valine, leucine, and isoleucine biosynthesis121
00300Lysine biosynthesis41
00310Lysine degradation65
00330Arginine and proline metabolism156
00340Lysine degradation44
00350Tyrosine metabolism259
00360Phenylalanine metabolism96
00361Gamma-hexachlorocyclohexane degradation56
00380Tryptophan metabolism2214
00400Phenylalanine, tyrosine, and tryptophan biosynthesis1313
00410B-alanine metabolism55
00440Aminophosphonate metabolism41
00450Selenoamino acid metabolism145
00472D-arginine and ornithine metabolism215
00480Glutathione metabolism82
00500Starch and sucrose metabolism1611
00520Nucleotide sugars metabolism74
00530Aminosugars metabolism75
00550Peptide protein biosynthesis170
00561Glycerolipid metabolism59
00562Inositol phosphate metabolism42
00600Sphingoglycolipid metabolism43
00620Pyruvate metabolism3510
00626Nitrobenzene degradation11
00630Glyoxylate and dicarboxylate metabolism94
00640Propanoate metabolism96
00643Styrene degradation01
00650Butanoate metabolism136
00670One-carbon pool by folate63
00680Methane metabolism353
00740Riboflavin metabolism55
00750Vitamine B6 metabolism01
00760Nicotinate and nicotinamide metabolism43
00770Pantothenate and CoA biosynthesis91
00790Folate biosynthesis42
00860Porphyrin and chlorophyll metabolism1806
00920Sulfur metabolism: reduction and fixation20
00950Alkaloid biosynthesis I70
00970Aminoacyl-tRNA biosynthesis112
03020RNA polymerase31
03030DNA polymerase40
03130Cytochrome C oxidase10
03140Cytochrome C reductase104
03150Succinate dehydrogenase31

Glycolysis/gluconeogenesis was the most represented pathway within each library. The second and third most enriched functional pathways in the C1 library were porphyrin and chlorophyll metabolism, which involved 180 ESTs (17.1%), and the citrate cycle, which involved 47 ESTs (4.5%). In the C2 library, the respective pathways were peptideprotein biosynthesis and d-arginine and ornithine metabolism.

It should be noted that the types of genes involved in the same metabolic pathways were greatly different between libraries. For instance, in the glycolysis/gluconeogenesis pathway, the enzymes in the C1 library were enolase, glyceraldehyde 3-phosphate dehydrogenase, fructose 1,6-biphosphate aldolase, and pyruvate decarboxylase, which were assigned to glycolysis; however, in the C2 library, they were fructose-1,6-bisphosphatase, pyruvate carboxylase, and phosphoenolpyruvate carboxykinase, which were assigned to gluconeogenesis. Because glucose is a very important source of nutrition, we speculate that the upregulated genes related to gluconeogenesis observed in the C2 library may be necessary for mycoparasitism, that is, maintenance of fast cell growth rate in response to the competition with the plant fungal pathogen.

3.4. Genes Induced by the Mycoparasitic Process

Differences were observed in gene groups associated with degradation of the cell wall, proteolysis, and toxins production. Seven contigs were presented in both libraries, four were specific to library C1, and eight to library C2 (Table 3), the latter appearing to be induced by the mycoparasitic process directly. Two sequences similar to the 48-KDa endochitinase (GenBank accession nos. DV546055, DV544732, and DV544989) of Aspergillus nidulans and 46-KDa endochitinase (DV547883 and DV547485) of Hypocrea virens were identified in libraries C1 and C2, respectively. Four ESTs from library C1 (DV546459, DV546294, DV544423, and DV546484) and 1 (DV548260) from library C2 shared similarity with serine proteases (Figure 1). One-gene homologue of MAP kinase A (DV548513) was identified in library C2 only (Figure 2).


48 KDa endochitinase30
46 KDa endochitinase02
β-glucosidase D11
β-glucosidase 560
β-glucosidase 610
Serine protease41
Aspartic proteinase205
Immune antigen 101
C-8 sterol isomerase22
Sterol C-22 desaturase11
Peptide transporter 201
MAP kinase A01
ABC transporters01

4. Discussion

It has been demonstrated that the cell wall of the fungal pathogen can simulate some aspects of the mycoparasitic interactions between biocontrol fungi and its targets [10].

Only a limited amount of overlap (333 unigenes) was observed in both libraries. A total of 1138 and 691 unigenes were unique and only expressed in C1 or C2, respectively. The lack of significant overlap between the individual libraries also suggests a high level of flexibility at the level of gene expression under the examined conditions, some of which may reflect particular requirements for phases of mycoparasitism.

The analysis of the frequency of specific ESTs that form individual contigs can give information about the expression levels of particular genes under different experimental conditions [11]. The most abundant transcripts in library C2 but no expression in C1 were MFS monosaccharide transporters. MFS transporters transport uni-, sym-, and antiporters of sugars, peptides, drugs, and organic and inorganic ions with 12 or 14 transmembrane spanners [12]. In the present study, the high proportion of ESTs expressing a homology to MFS monosaccharide transporters implies that they may be responsible for transport of monosaccharides derived from the degradation of RsCW. This was not consistent with the results of a study of T. harzianum CECT 2413 [13], in which abundant expression of peptide transporter 2 (PTR2) was found in a cDNA library of T. harzianum CECT 2413 when interacted directly with Botrytis cinerea. However, only one EST similar to PTR2 (DV547977) was detected in library C2. We speculate that this may have been caused by the different cultivation times of the two fungi.

Comparison analysis illustrated variations in the proportions of different pathways. Metabolic pathways of ubiquinone biosynthesis; electron transport and oxidative phosphorylation; purine metabolism; pyrimidine metabolism; alanine and aspartate metabolism; valine, leucine, and isoleucine biosynthesis; porphyrin and chlorophyll metabolism were proportionately more represented in library C1. In contrast, pentose and glucuronate interconversions; fructose and mannose metabolism; galactose metabolism; androgen and estrogen metabolism; glycine, serine, and threonine metabolism; valine, leucine, and isoleucine degradation; arginine and proline metabolism; histidine metabolism; tryptophan metabolism; d-arginine and ornithine metabolism; glycerolipid metabolism were overrepresented in library C2. Metabolic pathways of sterol, vitamin K, vitamin E, carotenoids biosynthesis; sulfur metabolism: reduction and fixation; DNA polymerase; cytochrome C oxidase were only observed in library C1, while those of fatty acid biosynthesis (path 2), styrene degradation, and Vitamine B6 metabolism were only observed in library C2.

The results showed that genes related to mycoparasitism were differentially expressed. Two different sequences similar to the 48-KDa endochitinase and 46-KDa endochitinase were identified in libraries C1 and C2, respectively. Since library C1 was obtained from cultivation on PDA medium, the 48 KDa endochitinase homolog might play a role in the dissolution and formation of the cell wall of C. cupreum. Similarly, because library C2 was constructed under conditions associated with mycoparasitism, the 46 KDa endochitinase homolog is expected to be involved in cell wall degradation of the fungal pathogen during the mycoparasitic process.

The conditions used for construction of library C2 were aimed at in vitro simulation of the mycoparasitic process, which is triggered by the recognition of the structural character of the pathogenic fungal cell wall. As a result, the genes involved in signal transduction pathways of mycoparasitism were acquired. Examples include homologue of gene encoding an ABC transporter (ATP-binding cassette transporter, DV548480, Figure 3) and MAP-kinase A (Tmk1 of T. atroviride). Four ESTs have sequence homology to an ABC transporter, it was also observed previously in other fungal pathogens (Gibberella pulicaris and Sclerotinia sclerotiorum) Mehrabi et al. [14] as potential pathogenicity factors responsible for tolerance to phytoalexins or a pathogenicity factor for the host Fleissner et al. [15] and Li et al. [16]. Studies on signal transduction pathways from Trichoderma strains revealed the involvement of MAP-kinases in the mycoparasitic interaction, including production of hydrolytic enzymes such as chitinases and secretion of antibiotic substances [17].

In this study, we sequenced and analyzed two independent cDNA libraries, providing the first comparative analysis of the transcriptome of C. cupreum under different conditions. The findings provide an entry point for understanding further the molecular mechanisms of this fungus and will also help to advance our efforts in developing novel strategies for biocontrol of fungal diseases.


  1. K. Soytong, Application of a New Broad Spectrum Biological Fungicide for Environmental Plant Protection, Biological Control and Biotechnology, HeiLongJiang Science and Technology Press, HeiLongJiang, China, 2003.
  2. Q. Yang, Evaluation of Chaetomium for Biological Control of Fusarium Wilt of Tomato in P. R. China, Biological Control and Biotechnology, HeiLongJiang Science and Technology Press, HeiLongJiang, China, 2003.
  3. M. D. Adams, J. M. Kelley, J. D. Gocayne et al., “Complementary DNA sequencing: expressed sequence tags and human genome project,” Science, vol. 252, no. 5013, pp. 1651–1656, 1991. View at: Google Scholar
  4. A. Bouck and T. Vision, “The molecular ecologist's guide to expressed sequence tags,” Molecular Ecology, vol. 16, no. 5, pp. 907–924, 2007. View at: Publisher Site | Google Scholar
  5. J. A. Vizcaíno, F. J. González, M. B. Suárez et al., “Generation, annotation and analysis of ESTs from Trichoderma harzianum CECT 2413,” BMC Genomics, vol. 7, p. 193, 2006. View at: Publisher Site | Google Scholar
  6. J. A. Vizcaíno, J. Redondo, M. B. Suárez et al., “Generation, annotation, and analysis of ESTs from four different Trichoderma strains grown under conditions related to biocontrol,” Applied Microbiology and Biotechnology, vol. 75, no. 4, pp. 853–862, 2007. View at: Publisher Site | Google Scholar
  7. H. Y. Zhang and Q. Yang, “Expressed sequence tags-based identification of genes in the biocontrol agent Chaetomium cupreum,” Applied Microbiology and Biotechnology, vol. 74, no. 3, pp. 650–658, 2007. View at: Publisher Site | Google Scholar
  8. D. Gordon, C. Abajian, and P. Green, “Consed: a graphical tool for sequence finishing,” Genome Research, vol. 8, no. 3, pp. 195–202, 1998. View at: Google Scholar
  9. M. Kanehisa and S. Goto, “KEGG: kyoto encyclopedia of genes and genomes,” Nucleic Acids Research, vol. 28, no. 1, pp. 27–30, 2000. View at: Google Scholar
  10. L. Sanz, M. Montero, J. Redondo, A. Llobell, and E. Monte, “Expression of an α-1,3-glucanase during mycoparasitic interaction of Trichoderma asperellum,” FEBS Journal, vol. 272, no. 2, pp. 493–499, 2005. View at: Publisher Site | Google Scholar
  11. D. J. Ebbole, Y. Jin, M. Thon et al., “Gene discovery and gene expression in the rice blast fungus, Magnaporthe grisea: analysis of expressed sequence tags,” Molecular Plant-Microbe Interactions, vol. 17, no. 12, pp. 1337–1347, 2004. View at: Google Scholar
  12. S. S. Pao, I. T. Paulsen, and M. H. Saier, “Major facilitator superfamily,” Microbiology and Molecular Biology Reviews, vol. 62, no. 1, pp. 1–34, 1998. View at: Google Scholar
  13. J. A. Vizcaíno, R. E. Cardoza, M. Hauser et al., “ThPTR2, a di/tri-peptide transporter gene from Trichoderma harzianum,” Fungal Genetics and Biology, vol. 43, no. 4, pp. 234–246, 2006. View at: Publisher Site | Google Scholar
  14. R. Mehrabi, T. van der Lee, C. Waalwijk, and G. H. J. Kema, “MgSlt2, a cellular integrity MAP kinase gene of the fungal wheat pathogen Mycosphaerella graminicola, is dispensable for penetration but essential for invasive growth,” Molecular Plant-Microbe Interactions, vol. 19, no. 4, pp. 389–398, 2006. View at: Publisher Site | Google Scholar
  15. A. Fleissner, C. Sopalla, and K. M. Weltring, “An ATP-binding cassette multidrug-resistance transporter is necessary for tolerance of Gibberella pulicaris to phytoalexins and virulence on potato tubers,” Molecular Plant-Microbe Interactions, vol. 19, no. 2, pp. 102–108, 2002. View at: Google Scholar
  16. R. Li, R. Roger, B. Lone et al., “Interaction of Sclerotinia sclerotiorum with a resistant Brassica napus cultivar: expressed sequence tag analysis identifies genes associated with fungal pathogenesis,” Fungal Genetics and Biology, vol. 41, no. 8, pp. 735–753, 2004. View at: Publisher Site | Google Scholar
  17. B. Reithner, R. Schuhmacher, N. Stoppacher, M. Pucher, K. Brunner, and S. Zeilinger, “Signaling via the Trichoderma atroviride mitogen-activated protein kinase Tmk1 differentially affects mycoparasitism and plant protection,” Fungal Genetics and Biology, vol. 44, no. 11, pp. 1123–1133, 2007. View at: Publisher Site | Google Scholar

Copyright © 2012 Haiyan Zhang and Min Li. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

927 Views | 687 Downloads | 1 Citation
 PDF  Download Citation  Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly and safely as possible. Any author submitting a COVID-19 paper should notify us at to ensure their research is fast-tracked and made available on a preprint server as soon as possible. We will be providing unlimited waivers of publication charges for accepted articles related to COVID-19. Sign up here as a reviewer to help fast-track new submissions.