# Transcriptional Profiling of ESTs from the Biocontrol FungusChaetomium cupreum

Authors: Haiyan Zhang; Min Li
Journal: The Scientific World Journal
 (2012)

---

## Abstract

Comparative analysis was applied to two cDNA/ESTs libraries (C1 and C2) fromChaetomium 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.


---

## Body

## 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. cupreumAmes 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 ofT. 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 fromC. 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

TheC. cupreumisolate was kindly provided by King Mongkut’s Institute of Technology of Ladkrabang, Thailand.Rhizoctonia solaniwas stored in our laboratory.C1 was constructed from mycelia ofC. 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 withR. 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 ofC. 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 (http://www.phrap.org/) 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.


## 2.1. Fungal Strains and Culture Conditions

TheC. cupreumisolate was kindly provided by King Mongkut’s Institute of Technology of Ladkrabang, Thailand.Rhizoctonia solaniwas stored in our laboratory.C1 was constructed from mycelia ofC. 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 withR. 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 ofC. 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 (http://www.phrap.org/) 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 anE-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 Table1. 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.Table 1
The difference of high-redundancy genes expression between cDNA libraries fromC. cupreum.
Contig no.AnnotationC1C2ESTs%ESTs%Contig525Glyceraldehyde-3-phosphate dehydrogenase1097.41141.69Contig531Coproporphyrinogen oxidase724.8950.6Contig30Predicted protein563.8100Contig433C-4 sterol methyl oxidase543.6700Contig313Predicted protein362.4500Contig2Pyruvate decarboxylase271.8400Contig39Predicted protein231.5600Contig507Xylulose-5-phosphate phosphoketolase211.4330.36Contig 152Predicted protein201.36161.93Contig523EF1-alpha translation elongation factor191.29121.45Contig42Thiazole biosynthetic enzyme191.2920.24Contig520Actin171.1640.48Contig495Glutamine synthetase171.1610.12Contig517Aspartic protease161.0930.36Contig317Ammonium transporter161.0900Contig419ATP synthase protein161.0900Contig4Alcohol dehydrogenase151.0240.48Contig474β-1,3-exoglucanase140.9520.24Contig494ATP citrate lyase140.9510.12Contig412Histone H2B130.8800Contig298Pyruvate kinase130.8800Contig310Aspartate aminotransferase110.7510.12Contig408Heat shock protein 30110.7500Contig489ADP-ATP translocase110.7580.97Contig165MFS monosaccharide transporter00202.42Contig162Phosphoenolpyruvate carboxykinase00121.45Relative 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 (Table2).Table 2
The metabolic pathways difference between cDNA libraries fromC. cupreum.
Map no.PathwaysESTs in C1ESTs in C200010Glycolysis/gluconeogenesis1883100020Citrate cycle47900030Pentose phosphate cycle35600040Pentose and glucuronate interconversions6400051Fructose and mannose metabolism41000052Galactose metabolism7700053Ascorbate and aldarate metabolism5500061Fatty acid biosynthesis (path 1)3100062Fatty acid biosynthesis (path 2)0100071Fatty acid metabolism241200100Sterol, vitamin K, vitamin E, carotenoids biosynthesis1000120Bile acid biosynthesis191200130Ubiquinone biosynthesis14300150Androgen and estrogen metabolism41200190Electron transport and oxidative phosphorylation34600220Urea cycle and metabolism of amino groups10000230Purine metabolism41100240Pyrimidine metabolism6100251Glutamate metabolism24400252Alanine and aspartate metabolism14100260Glycine, serine, and threonine metabolism2800271Methionine metabolism12500272Cysteine metabolism6100280Valine, leucine, and isoleucine degradation6600290Valine, leucine, and isoleucine biosynthesis12100300Lysine biosynthesis4100310Lysine degradation6500330Arginine and proline metabolism15600340Lysine degradation4400350Tyrosine metabolism25900360Phenylalanine metabolism9600361Gamma-hexachlorocyclohexane degradation5600380Tryptophan metabolism221400400Phenylalanine, tyrosine, and tryptophan biosynthesis131300410B-alanine metabolism5500440Aminophosphonate metabolism4100450Selenoamino acid metabolism14500472D-arginine and ornithine metabolism21500480Glutathione metabolism8200500Starch and sucrose metabolism161100520Nucleotide sugars metabolism7400530Aminosugars metabolism7500550Peptide protein biosynthesis17000561Glycerolipid metabolism5900562Inositol phosphate metabolism4200600Sphingoglycolipid metabolism4300620Pyruvate metabolism351000626Nitrobenzene degradation1100630Glyoxylate and dicarboxylate metabolism9400640Propanoate metabolism9600643Styrene degradation0100650Butanoate metabolism13600670One-carbon pool by folate6300680Methane metabolism35300740Riboflavin metabolism5500750Vitamine B6 metabolism0100760Nicotinate and nicotinamide metabolism4300770Pantothenate and CoA biosynthesis9100790Folate biosynthesis4200860Porphyrin and chlorophyll metabolism180600920Sulfur metabolism: reduction and fixation2000950Alkaloid biosynthesis I7000970Aminoacyl-tRNA biosynthesis11203020RNA polymerase3103030DNA polymerase4003130Cytochrome C oxidase1003140Cytochrome C reductase10403150Succinate dehydrogenase31Glycolysis/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 (Table3), the latter appearing to be induced by the mycoparasitic process directly. Two sequences similar to the 48-KDa endochitinase (GenBank accession nos. DV546055, DV544732,andDV544989) of Aspergillus nidulansand 46-KDa endochitinase (DV547883andDV547485) of Hypocrea virenswere identified in libraries C1 and C2, respectively. Four ESTs from library C1 (DV546459, DV546294, DV544423, andDV546484) 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).Table 3
The different expression of genes related to biocontrol between cDNA libraries fromC. cupreum.
AnnotationC1C2ESTsESTs48 KDa endochitinase3046 KDa endochitinase02β-N-acetylglucosaminidase01β-1,3-exoglucanase222β-1,3-endoglucanase10β-glucosidase D11β-glucosidase 560β-glucosidase 610Serine protease41Aspartic proteinase205Immune antigen 101C-8 sterol isomerase22Sterol C-22 desaturase11Sterol-C5-desaturase22β-endo-1,4-xylanase03Peptide transporter 201MAP kinase A01ABC transporters01Figure 1
DV546294_shared 90% identity with serine proteases ofAcremonium sp.Figure 2
DV548513 shared 98% identity with MAP kinase ofGibberella intermedia.


## 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 anE-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 Table1. 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.Table 1
The difference of high-redundancy genes expression between cDNA libraries fromC. cupreum.
Contig no.AnnotationC1C2ESTs%ESTs%Contig525Glyceraldehyde-3-phosphate dehydrogenase1097.41141.69Contig531Coproporphyrinogen oxidase724.8950.6Contig30Predicted protein563.8100Contig433C-4 sterol methyl oxidase543.6700Contig313Predicted protein362.4500Contig2Pyruvate decarboxylase271.8400Contig39Predicted protein231.5600Contig507Xylulose-5-phosphate phosphoketolase211.4330.36Contig 152Predicted protein201.36161.93Contig523EF1-alpha translation elongation factor191.29121.45Contig42Thiazole biosynthetic enzyme191.2920.24Contig520Actin171.1640.48Contig495Glutamine synthetase171.1610.12Contig517Aspartic protease161.0930.36Contig317Ammonium transporter161.0900Contig419ATP synthase protein161.0900Contig4Alcohol dehydrogenase151.0240.48Contig474β-1,3-exoglucanase140.9520.24Contig494ATP citrate lyase140.9510.12Contig412Histone H2B130.8800Contig298Pyruvate kinase130.8800Contig310Aspartate aminotransferase110.7510.12Contig408Heat shock protein 30110.7500Contig489ADP-ATP translocase110.7580.97Contig165MFS monosaccharide transporter00202.42Contig162Phosphoenolpyruvate carboxykinase00121.45Relative 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 (Table2).Table 2
The metabolic pathways difference between cDNA libraries fromC. cupreum.
Map no.PathwaysESTs in C1ESTs in C200010Glycolysis/gluconeogenesis1883100020Citrate cycle47900030Pentose phosphate cycle35600040Pentose and glucuronate interconversions6400051Fructose and mannose metabolism41000052Galactose metabolism7700053Ascorbate and aldarate metabolism5500061Fatty acid biosynthesis (path 1)3100062Fatty acid biosynthesis (path 2)0100071Fatty acid metabolism241200100Sterol, vitamin K, vitamin E, carotenoids biosynthesis1000120Bile acid biosynthesis191200130Ubiquinone biosynthesis14300150Androgen and estrogen metabolism41200190Electron transport and oxidative phosphorylation34600220Urea cycle and metabolism of amino groups10000230Purine metabolism41100240Pyrimidine metabolism6100251Glutamate metabolism24400252Alanine and aspartate metabolism14100260Glycine, serine, and threonine metabolism2800271Methionine metabolism12500272Cysteine metabolism6100280Valine, leucine, and isoleucine degradation6600290Valine, leucine, and isoleucine biosynthesis12100300Lysine biosynthesis4100310Lysine degradation6500330Arginine and proline metabolism15600340Lysine degradation4400350Tyrosine metabolism25900360Phenylalanine metabolism9600361Gamma-hexachlorocyclohexane degradation5600380Tryptophan metabolism221400400Phenylalanine, tyrosine, and tryptophan biosynthesis131300410B-alanine metabolism5500440Aminophosphonate metabolism4100450Selenoamino acid metabolism14500472D-arginine and ornithine metabolism21500480Glutathione metabolism8200500Starch and sucrose metabolism161100520Nucleotide sugars metabolism7400530Aminosugars metabolism7500550Peptide protein biosynthesis17000561Glycerolipid metabolism5900562Inositol phosphate metabolism4200600Sphingoglycolipid metabolism4300620Pyruvate metabolism351000626Nitrobenzene degradation1100630Glyoxylate and dicarboxylate metabolism9400640Propanoate metabolism9600643Styrene degradation0100650Butanoate metabolism13600670One-carbon pool by folate6300680Methane metabolism35300740Riboflavin metabolism5500750Vitamine B6 metabolism0100760Nicotinate and nicotinamide metabolism4300770Pantothenate and CoA biosynthesis9100790Folate biosynthesis4200860Porphyrin and chlorophyll metabolism180600920Sulfur metabolism: reduction and fixation2000950Alkaloid biosynthesis I7000970Aminoacyl-tRNA biosynthesis11203020RNA polymerase3103030DNA polymerase4003130Cytochrome C oxidase1003140Cytochrome C reductase10403150Succinate dehydrogenase31Glycolysis/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 (Table3), the latter appearing to be induced by the mycoparasitic process directly. Two sequences similar to the 48-KDa endochitinase (GenBank accession nos. DV546055, DV544732,andDV544989) of Aspergillus nidulansand 46-KDa endochitinase (DV547883andDV547485) of Hypocrea virenswere identified in libraries C1 and C2, respectively. Four ESTs from library C1 (DV546459, DV546294, DV544423, andDV546484) 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).Table 3
The different expression of genes related to biocontrol between cDNA libraries fromC. cupreum.
AnnotationC1C2ESTsESTs48 KDa endochitinase3046 KDa endochitinase02β-N-acetylglucosaminidase01β-1,3-exoglucanase222β-1,3-endoglucanase10β-glucosidase D11β-glucosidase 560β-glucosidase 610Serine protease41Aspartic proteinase205Immune antigen 101C-8 sterol isomerase22Sterol C-22 desaturase11Sterol-C5-desaturase22β-endo-1,4-xylanase03Peptide transporter 201MAP kinase A01ABC transporters01Figure 1
DV546294_shared 90% identity with serine proteases ofAcremonium sp.Figure 2
DV548513 shared 98% identity with MAP kinase ofGibberella intermedia.


## 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 ofC. 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 atin 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 andSclerotinia 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 Trichodermastrains revealed the involvement of MAP-kinases in the mycoparasitic interaction, including production of hydrolytic enzymes such as chitinases and secretion of antibiotic substances [17].Figure 3
DV548480 shared 73% identity with ABC transporter ofMyceliophthora thermophila ATCC 42464.In this study, we sequenced and analyzed two independent cDNA libraries, providing the first comparative analysis of the transcriptome ofC. 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.



---

*Source: 340565-2012-02-01.xml*