Identification of mycoparasitism-related genes against the phytopathogen Sclerotinia sclerotiorum through transcriptome and expression profile analysis in Trichoderma harzianum
© Steindorff et al.; licensee BioMed Central Ltd. 2014
Received: 28 November 2013
Accepted: 6 March 2014
Published: 18 March 2014
The species of T. harzianum are well known for their biocontrol activity against plant pathogens. However, few studies have been conducted to further our understanding of its role as a biological control agent against S. sclerotiorum, a pathogen involved in several crop diseases around the world. In this study, we have used RNA-seq and quantitative real-time PCR (RT-qPCR) techniques in order to explore changes in T. harzianum gene expression during growth on cell wall of S. sclerotiorum (SSCW) or glucose. RT-qPCR was also used to examine genes potentially involved in biocontrol, during confrontation between T. harzianum and S. sclerotiorum.
Data obtained from six RNA-seq libraries were aligned onto the T. harzianum CBS 226.95 reference genome and compared after annotation using the Blast2GO suite. A total of 297 differentially expressed genes were found in mycelia grown for 12, 24 and 36 h under the two different conditions: supplemented with glucose or SSCW. Functional annotation of these genes identified diverse biological processes and molecular functions required during T. harzianum growth on SSCW or glucose. We identified various genes of biotechnological value encoding proteins with functions such as transporters, hydrolytic activity, adherence, appressorium development and pathogenesis. To validate the expression profile, RT-qPCR was performed using 20 randomly chosen genes. RT-qPCR expression profiles were in complete agreement with the RNA-Seq data for 17 of the genes evaluated. The other three showed differences at one or two growth times. During the confrontation assay, some genes were up-regulated during and after contact, as shown in the presence of SSCW which is commonly used as a model to mimic this interaction.
The present study is the first initiative to use RNA-seq for identification of differentially expressed genes in T. harzianum strain TR274, in response to the phytopathogenic fungus S. sclerotiorum. It provides insights into the mechanisms of gene expression involved in mycoparasitism of T. harzianum against S.sclerotiorum. The RNA-seq data presented will facilitate improvement of the annotation of gene models in the draft T. harzianum genome and provide important information regarding the transcriptome during this interaction.
Sclerotinia sclerotiorum (Lib.) de Bary is one of the most devastating and cosmopolitan plant pathogens. This fungus infects over 400 species of plants worldwide including important crops and numerous weeds . S. sclerotiorum poses a threat to crops such as sunflower, soybean, oilseed rape, edible dry bean, chickpea, peanut, dry pea, lentils, onion and tulip . It is capable of infecting flowers, leaves, fruits or stems  and its life cycle initiates by the germination of sclerotia and formation of infectious propagules. During the crop growing season, which dependson a set of environmental factors, fungal sclerotia germinate to form mycelia, which can directly infect host plants, or produce ascospores. Further, ascospores develop forming apothecia . Ascospores are the primary infective propagules of S. sclerotiorum on many crops, and also can act in disease scattering since they may be transported to neighboring fields as well as over longer distances .
Chemical treatment is today the main strategy employed worldwide for fungal disease control. In order to develop alternative and sustainable methods for control of white mold, which do not cause negative environmental or economic impacts during crop production, as observed with the routine use of fungicides, our research group has isolated strains of a number of Trichoderma species from diverse agro-ecosystems in Brazil and assessed their potential for biocontrol of S. sclerotiorum. This analysis has included evaluation of antagonistic capacity, production of cell wall-degrading enzymes and production of volatile antibiotics . Our results have identified T. harzianum Rifai (anamorph) strain TR274 as a promising biocontrol agent against S.sclerotiorum in vitro and under field conditions [5, 6].
Biological control is a complex process which requires the host to be recognised by the biocontrol agent, followed by hydrolytic enzyme and antibiotic production which is triggered by the direct attachment of the mycoparasite to the host fungus. This contact is mediated by lectins and proteins harboring cellulose binding modules from hyphae of the host and mycoparasitic fungus, respectively, thereby eliciting a signaling cascade comprising G-proteins and MAPKs that can modulate Trichoderma’s protein expression pattern [7, 8]. However, the detailed molecular mechanisms involved in this process remain unknown. This complex mechanism is influenced by the pathogen and Trichoderma isolates evaluated . In this context, studies conducted on different strains are necessary for increased understanding of the biocontrol mechanism.
The sequencing of expressed-sequence-tag (EST) libraries for different Trichoderma strains cultivated in the presence of host fungi has contributed significantly to the large-scale identification ofmycoparasitism-related genes [10–12]. Our research group has described gene mapping using EST and suppression subtractive hybridization (SSH) approaches [12, 13] during the interaction of T. harzianum with Fusarium solani and proteomic approaches for T. harzianum grown in liquid containing Rhizoctonia solani, Macrophomina phaseolina and Fusarium sp cell walls . Whilst DNA microarrays have been used to study the interaction among Trichoderma and host plants , only two studies have employed high-throughput transcriptomic approaches to investigate mycoparasitism mechanisms of Trichoderma[16, 17]. RNA sequencing (RNA-Seq), a high-throughput technology used to sequence complementary DNA, has been widely thought of as a revolutionary tool for transcriptomics [18–21]. The publically available whole genome sequence for T. harzianum CBS 226.95 , which was recently released by the Joint Genome Institute (JGI) (http://genome.jgi.doe.gov/Triha1/Triha1.home.html), now allows for use of RNA-Seq approaches and mapping of data to the reference sequence, which will likely contribute to identification of mycoparasite-related genes, as well as the molecular mechanisms by which this fungus is able to inhibit phytopathogen fungal growth.
The present study is the first initiative to use RNA-seq for identification of differentially expressed genes in T. harzianum strain TR274, in response to the phytopathogenic fungus S. sclerotiorum. T. harzianum TR274 was cultivated on liquid medium containing S. sclerotiorum cell wall (SSCW) to mimick fungal host presence. Quantitative real-time PCR (RT-qPCR) supported in silico-based evidence for differential gene expression of candidate genes involved in mycoparasitism.
Results and discussion
Illumina sequencing and mapping onto the T. harzianumreference genome
In this present study, an RNA-seq approach was used to map genes differentially expressed during T. harzianum growth on SSCW. The experimental design enabled comparison of gene expression in the presence of host cell wall or glucose as sole carbon source. Samples of mRNA from T. harzianum following three growth periods in the presence of SSCW (12, 24 and 36 h) were used to construct six Illumina libraries. A total of 171,442,148 sequence reads were obtained after quality trimming, varying from 25 to 100 bp in length. Each library was represented by at least 16,845,349 reads, representing a coverage of 66X when compared with the full transcriptome, a density regarded as adequate to perform gene expression analysis . Complete mapping information can be accessed in Additional file 1: Table S8.
Sequence reads were aligned onto the T. harzianum CBS 226.95 reference genome . This strain was isolated from garden soil in the UK, while strain TR274 was isolated from Brazilian cerrado soil. Despite the genetic differences among T. harzianum isolates described in the literature [24, 4] and the fact that the reference genome published is only a first draft, 82.6% ± 10.08% of the obtained reads were mapped onto the reference genome using the default settings of the Bowtie2 aligner. The high mapping percentage suggests a high similarity between these isolates, at least at the transcriptome level. Only 0.8% ± 0.15% of reads for each library was mapped in more than one region on the reference genome and these reads were not used in the expression analysis.
Gene expression analysis using “in silico” approach
The most notable increase in the gene expression was observed between 12 and 24 hours with a noteworthly change in expression pattern and protein functions detected between 12 and 36 hours. The main genes up or downregulated after 12 and 36 h of growth encoded CAZymes and transporters (Additional file 2: Table S1 and Additional file 3: Table S4). This suggests that some proteins required during early growth phases which are involved in cell wall degradation and sugar transport are no longer necessary after 36 hours of growth.
log 2 Fold Change (FC) of the top 10 differentially expressed genes after 12, 24 and 36 h
12 h up
Acid proteinase protein
Alkaline serine protease
Isotrichodermin c-15 hydroxylase
Hypothetical protein TRIVIDRAFT_62551
Hypothetical protein TRIVIDRAFT_130513
Uncharacterized serine-rich protein
24 h up
Transcriptional activator srcap-like protein
Conidiation-specific protein (con-13) protein
Hypothetical protein FOC4_g10000877
wsc domain-containing protein
Peptidase family m28
Related to spore coat protein sp96 precursor
36 h up
Fasciclin domain containing protein
Peptidase family m28
Glycoside hydrolase family 76
Proteinase inhibitor kazal
Endonuclease/exonuclease/phosphatase family protein
MFS carboxylic acid transport protein
Alkaline proteinase/serin endopeptidase
The down-regulated transcripts for all the stages of growth were categorized into oxidoreductase activity, oxidation-reduction process and some “binding” child categories, such as small molecule binding and nucleotide binding proteins. A hypothesis for repression could be the nature of basal metabolism of the genes belonging to these categories. This fact is consistent with the extensive metabolic activity expected for a filamentous fungus growing on a rich medium (glucose 2% medium) with an easily assimilable substrate . Under this culture condition up-regulation of a specific subset of oxidoreductases and nucleotide binding proteins endoding genes related to primary metabolism is commonly observed for Trichoderma species, which are down-regulated in the presence of complex substrates . Vieira et al. in 2013 showed the same pattern of repressed categories when T. harzianum was grown on Fusarium solani cell wall, suggesting that this result is not pathogen-dependent.
The functional category distribution frequency for each cluster was then calculated to identify differences in the distribution of genes among the three Trichoderma growth periods (Figure 3D). High percentages of hydrolases, peptidases and transporter activities were observed in clusters 1, 2 and 3, which included most of the genes up-regulated for all growth times, mainly after 12 h. Clusters 4 and 5 are represented by the lowest number of genes; however they presented a high diversity of functional categories and the smallest percentage of the hydrolase activity category. Cluster 4 is mainly represented by specific transporters, oxidoreductase and peroxidase activities, which are absent in the other clusters. Cluster 5 was less diverse, but contains a high percentage of peptidases and the appearance of the lyase activity category. Clusters 4 and 5 contain genes induced after 24 and 36 hours. This pattern of categorization suggests a late increasing in the expression levels of genes encoding substrate transporters and other CAZYmes such as lyases, instead of the classical cell wall hydrolases (chitinases, β-1,3-glucanases and proteases). Recently, we have described the production of polysaccharide lyases by T. harzianum during growth in the presence of different phytopathogen cell walls using a proteomic approach, which is in agreement with the presented data in this work (unpublished data). These results suggest a potential role of these enzymes, especially the lyase families 7 and 8, in mycoparasitism by Trichoderma independent of the host pathogen.
Validation of RNA-seq gene expression
As a set of peptidases were up-regulated in the presence of SSCW, two serine peptidases genes (tripsin-like peptidase and serine endopeptidase 33 kDa), one aspartyl protease gene, one aminopeptidase gene (Peptidase M28) and one carboxypeptidase gene (Peptidase M14) were chosen to perform qPCR experiments. The expression data over the five growth times revealed an upward trend along the time course with the highest expression values at 18 hours followed by a decrease after this time or 24 hours. Their diversity and uniform time course production indicated that these enzymes may form a synergistic proteolytic system related to mycoparasistim in T. harzianum. However, their exact role in mycoparasitism has not been clearly elucidated yet. The main accepted hypothesis presents these enzymes as contributing to the breakdown of the fungal host cell wall, constituted by chitin and glucan polymers embedded in, and covalently linked to, a protein matrix , and also as acting as proteolytic inactivators of pathogen enzymes which are typically involved in the plant infection process .
The Trichoderma species, T. reesei, T. atroviride, and T. virens, may have one of the largest sets of proteases among fungi, of which the total number of predicted proteases are 383 (4.2% of all predicted protein coding genes), 445 (3.75%), and 479 (3.85%), respectively . These authors described that the dominant groups are classified as aspartyl proteases, serine proteases, subtilisin-like proteases, and dipeptidyl and tripeptidyl peptidases. These enzymes have been described as performing a central role in the mycoparasitic activity of Trichoderma species, as they have been consistently identified during interaction against different phytopathogenic fungi using transcriptomic and proteomic approaches [12–14]. Their diversity and uniform time course production provide evidence that these enzymes may form a synergistic proteolytic system related to mycoparasitism in this Trichoderma species. Indeed, in our study, a number of genes encoding a serine peptidase, anaspartyl protease, a subtilisin-like, a trypsin-like beyond metallopeptidases (M28, M14) were also differentially expressed in the presence of SSCW. These results strongly suggest a role of these enzymes in T. harzianum mycoparasitism against S. sclerotiorum, and support a putative common action mechanism of mycoparasite fungi within the Trichoderma genus.
About 29% of the differentially expressed genes encode CAZy enzymes (Additional file 2: Table S1), suggesting a central role for them in Trichoderma mycoparasitism, probably during host fungal cell wall degradation. The kinetics expression for CAZy category enzymes presented as expression average showed a decrease trend along the growth timecourse, with a maximal expression at 24 h. An exception were the enzymes categorized as auxiliary activities and polysaccharide lyases (Additional file 5: Table S3). Glycosyl hydrolase family 18 and other enzymes which act mainly as fungal cell wall degrading enzymes (CWDE) have also been described as presenting a central role in mycoparasitism in Trichoderma atroviride. Genes encoding this family were also up-regulated in the present work based on “in silico” RNA-seq data analysis. As a consequence, eight genes encoding three GH18 (Chitinase 37 kDa, Chitinase 42 kDa and a not well characterized Chitinase), a α-1,3-glucanase, a β-glicosidase, a lyzosyme-like, a polyssacharide lyase Family 7 and a carbohydrate binding Family 13 protein were selected for expression validation by RT-qPCR. The expression levels of all CWDE genes followed the same trend, except for the β-glucosidase gene that was repressed from 12 hours onwards. The common expression profile was an increase in transcripts after 12 hours until 24 hours and a similar level or decrease of transcripts at 36 hours. This kinetic suggests use of small sugars during the first growth time and an expression of CWDE after 12 hours, indicating a role in degradation of the cell wall as carbon source to allow continued growth. KEGG pathway analysis (Additional file 6: Table S7) shows the mapped genes in starch and sucrose metabolism. All genes mapped are in the final stages of the pathway, mainly in the formation of small sugars such as D-glucose and D-xylose.
Alginate lyases are enzymes classified as belonging to the polyssacharide lyase Family 7 and are usually involved in the deconstruction of complex polyssacharides, such as polyguluronate and polymannuronate . Their expression observed in our study suggests a possible role in mycoparasitism, complementing the classical GH activity which is known to be involved. The carbohydrate binding module Family 13 gene encodes a protein with a domain related to lectins, which, in Rhizoctonia solani, is implicated in fungal insecticidal activity . Our expression data suggest that this gene is related to the lectins and may play a role mediating the physical contact with the host and elicitation of the signaling cascade comprising G-proteins and MAPKs.
Small secreted cysteine-rich proteins (SSCPs) have been described as up-regulated in Trichoderma species during mycoparasitism against phytopathogenic fungi . The present work is the first to report their probable role in mycoparasitism of T. harzianum against S. sclerotiorum.
In this work we also indentified genes enconding two predicted cistein-rich proteins and qid74, all up-regulated in the presence of SSCW. These genes were highly expressed in the later induction time periods (24 and 36 h) as showed by “in silico” analysis of RNA-seq and are in agreement with the results of validation by RT-qPCR. Hydrophobins I and II were also identified by the “in silico” analysis of RNA-seq data, as shown in Additional file 4: Table S2 with a high expression at 24 and 36 hours.
Small secreted cysteine-rich proteins (SSCPs) are one of the largest groups of proteins secreted by Trichoderma. Hydrophobins, probably the most widely known SSCPs, are characterized by the presence of eight positionally conserved cysteine residues of which four occur in doubles. They are found on the outer surfaces of cell walls of hyphae and conidia, where they mediate interactions between the fungus and the environment and also between the fungus and host plant roots . Class II hydrophobins represent the predominant class described for Trichoderma species . T. atroviride and T. virens have also class I-like hydrophobins, however they present differences in several aspects when compared to other fungi and form a separate clade in phylogenetic analysis within the Ascomycetes . As well as other cysteine-rich proteins, T. harzianum the qid74 gene encodes a cell wall protein which has an important role in adherence to hydrophobic surfaces and cellular protection . The gene RNaseT2, which has been described as a stress related protein and involved in permeability and stability of the plasmatic membrane in Saccharomyces cerevisiae, was also up-regulated after 12 hours growth of T.harzianum on SSCW, decreasing over time.
Among the 297 differentially expressed genes, 30 encode transporter proteins (Additional file 3: Table S4). MFS (Major Facilitator Superfamily) permeases are the most abundant proteins among transporter proteins over the three growth times. These proteins enable the transport of essential nutrients and ions, in addition to the excretion of end products of metabolism and cell-environment communication . Their expression levels vary according to the time of growth and culture growth condition (presence or absence of SSCW and glucose).
In summary, our results demonstrated a time course dependent expression of T. harzianum genes during growth on media with S. sclerotiorum cell wall as sole carbon and nitrogen source. The majority of the genes described in this work have also been reported in the literature during growth of Trichoderma spp. on cell-wall of phytopathogenic fungi [10, 12, 13], as well as during confrontation against R. solani.
RT-qPCR for dual cultures
The two small cystein-rich proteins were up-regulated during the interaction between T. harzianum and S. sclerotiorum although they showed differing expression patterns over the time period (Figure 6B). The gene cystein rich (511478) was highly expressed after contact, whilst the gene small cystein rich (518220) before contact. These data suggests that this group of proteins (SSCP) may present a synergistic time course dependent activity during the interaction. The RNAse T2 gene was also up-regulated during the contact stage, confirming its role in this interaction, possibly through conferring membrane stability during contact with the phytopathogen hyphae (Figure 6B).
Among the eight genes encoding putative glycoside hydrolases, three were up-regulated in the interaction: polyssacharidelyase Family 7, β-glicosidase and the CB module Family 13. This data suggests a role of these proteins in the mycoparasitism, complementing the classical GH activity in the interaction. The carbohydrate binding module Family 13 gene was induced during and after contact with S. sclerotiorum. The other genes identified as up-regulated by RNA-seq “in silico” analysis and validated by qRT-PCR using SCCW were not up-regulated in this condition. This “not-complete” agreement between cell wall induction with SSCW and direct confrontation with S. sclerotiorum is expected due to differences in growth conditions for T. harzianum. Investigation of mycoparasitisn using inactivated cell walls is a useful approach to evaluate the broad arsenal of induced genes in Trichoderma spp. and to identify candidate mycoparasitism related genes that can be further evaluated for expression during mycoparasitism through dual culture assays. The genes which were observed to be up-regulated in both interaction bioassays are certainly promising candidates for future biotechnological applications as well as further detailed investigations to unveil their precise function in T. harzianum mycoparasistism.
The RNA-seq data presented will not only facilitate improvement of the annotation of gene models in the draft T. harzianum genome, but also provide important information regarding the transcriptome during growth on SSCW and during “in vivo” interactions with S. sclerotiorum. Our data represent an important step towards understanding the mycoparasitic process of T. harzianum during its interaction with S. sclerotiorum. Further studies for functional characterization of candidate genes reported here are necessary in order to better define the exact pathways involved in mycoparasitism in T. harzianum.
Organism culture conditions
T. harzianum TR274 was isolated from soil samples from the Brazilian Cerrado region and identified to species level based upon ribosomal RNA ITS 1 and 2 sequence identities (Genbank number KC993076). This strain is deposited in the ICB Enzymology Group culture collection at the Universidade Federal de Goias. A strain of the phytopathogenic fungus S. sclerotiorum belonging to the EMBRAPA-CNPAF culture collection was originally isolated from Phaseolus vulgaris. These fungi were grown on MYG medium (w/v: 0.5% malt extract, 0.25% yeast extract, 1% glucose and 2% agar) for 2 days at 28°C. T. harzianum spores were collected from cultures by washing with sterile water and centrifugation at 2000 rpm. Following two rounds of washes, spore suspensions at a concentration of 107 spores mL−1 were used to inoculate flasks containing 50 mL of TLE medium . Cultures were incubated at 28°C with constant shaking at 120 rpm. After 24 hours growth, mycelia were collected and transferred to flasks containing 50 mL of minimal medium (TM) (w/v: 0.2% KH2PO4, 1.4% (NH4)2SO4, 0.03% MgSO4.7H2O) supplemented with 2% (w/v) of glucose or 0.5% (w/v) of inactivated SSCW (previously autoclaved at 120°C for 20 min). Cultures were incubated with constant shaking at 120 rpm at 28°C. After 6, 12, 18, 24 and 36 hours of growth, mycelium was harvested and immediately flash frozen in liquid nitrogen and stored at -80°C until RNA isolation. This experiment was carried out in triplicate for each growth period, with mycelia subsequently pooled to constitute composed samples.
RNA isolation, cDNA library preparation and sequencing
RNA was extracted from macerated frozen mycelia using TRI-Reagent (Sigma-Aldrich), according to manufacturer’s instructions. Integrity and quantity of isolated RNA were assessed using a RNA Pico chip on an Agilent Bioanalyzer 2100 (Agilent Technologies) (Additional file 7: Table S5). The RNA samples were stored at -80°C until library construction for sequencing and RT-qPCR. Normalized starting quantities of total RNA extracted from mycelia grown for 12, 24 and 36 h under the two different conditions (supplemented with glucose or SSCW), were then used to prepare six separate barcoded RNA-seq libraries using the TruSeq™ RNA sample preparation kit (Illumina, CA, USA). Three biological replicates were pooled in preparation of each of the six final samples. All library preparation and sequencing was carried out by Eurofins MWG Operon (Al, USA). Barcoded libraries were prepared according to the manufacturer’s instructions, then sequenced on a single lane on an Illumina HiSeq2000 sequencer, to generate 100 bp paired-end reads.
Mapping of sequenced reads and assessment of gene expression
FastQCfiles were used to visualize the libraries quality before and after trimming. For quality trimming and sequence filtering, the software Trimmomatic (version 0.30)  was employed to remove sequencing adapters, k-mers, and bases with a Phred quality score lower than 20 from the read ends. All reads shorter than 25 nucleotides were then discarded. Filtered reads were mapped onto the Trichoderma harzianum v1.0 genome sequence (http://genome.jgi.doe.gov/Triha1/Triha1.home.html) using TopHat 2.0.8 release with default settings . Gene expression values were determined using Cufflinks 2.1.1 release , and the FPKM (Fragments per kilobase mapped) values were calculated for each transcript. The gene expression levels in the T. harzianum genome were obtained using the Cuffdiff software within Cufflinks. Fungal transcript levels were calculated using uniquely mapped reads onto the genome. The expression profiles of differentially expressed genes were determined by SOTA (Self Organizing Tree Algorithm), cluster analysis was carried out using the MeV 4.9 software, and gene ontology assignment conducted according to Gene Ontology (GO) guidelines  using the Blast2GO platform .
Each T. harzianum sample grown in SSCW had a corresponding control sample of culture growth on glucose, enabling normalization using the relation FPKM SSCW/FPKM glucose. Positive values were considered as up-regulated transcripts in the presence of cell wall and negative values were considered as down-regulated transcripts. Expression level differences were judged to be significant and the expression level estimate status acceptable when a gene was identified as differentially expressed with an FDR of the Benjamini-Hochberg multiple tests of 5% (P < 0.05).
Gene regulatory network of T. harzianum
A regulatory network was generated using Cytoscape 3.0.2 software and a table of data containing the differentially expressed genes detected using the cuffdiff program, presenting Fold Change > 2, the interaction type (up- or down-regulated) and the target gene (i.e., the Protein ID of each gene affected). This analysis was carried out in order to reconstruct a T. harzianum time course network representation for all 297 identified genes (Figure 1B) .
Twenty genes were randomly selected between the differentially expressed genes in silico for the expression analysis by RT-qPCR. Real-time qPCR (Additional file 8: Table S6) primers were designed using the PerlPrimer v1.1.20 software. Total RNA from the above described preparations was digested with DNase I (Invitrogen) and a total of 5 μg from each pooled sample was reverse transcribed into cDNA using the Maxima™ First Strand cDNA synthesis kit for RT-qPCR in a final volume of 20 μL (Fermentas). The synthesized cDNA was diluted with 80 μL of water and used as a template for real-time PCR reactions using the instrument iQ5 real-time PCR system (Bio-Rad). Each reaction (20 μL) contained 10 μL of MAXIMA® SYBR-green PCR Master Mix (Fermentas), forward and reverse primers (500 nM each), cDNA template, and nuclease free water. PCR cycling conditions were: 10 min at 95°C (1 cycle), 15 s at 95°C followed by 1 min at 60°C (40 cycles), and a melting curve ramping from 60°C to 95°C with an increasing temperature of 0.2°C for 10 s (1 cycle) and continuous data collection to test for primer dimers and nonspecific amplification. Determination of the PCR efficiency was performed using triplicate reactions from a dilution series of cDNA (1, 0.1, 10−2, and 10−3). Amplification efficiency was then calculated from the given slopes in the IQ5 Optical System Software v2.0 (Additional file 8: Table S6). The α-tubulin and β-actin transcript were used as internal references to normalize the amount of total RNA present in each reaction . Gene expression levelswere calculated from the threshold cycle according to the 2-ΔΔCT method . All samples were analyzed in at least two independent experiments with three technical replicates in each run.
Analysis of expression of biocontrol-related genes
RT-qPCR was used to evaluate T. harzianum gene expression during confrontation against the fungal pathogen S. sclerotiorum. Confrontation bioassays were conducted as described in Steindorff et al. . T. harzianum TR274 and S. sclerotiorum circular plaques of 5 mm diameter were cut from 7-day-old cultures on MYG Plates. T. harzianum TR274 plaques were inoculated onto plates containing solid minimal medium supplemented with 0.2% of glucose and overlaid with cellophane at a distance of 7 cm from S. sclerotiorum plaque mycelia. A control confrontation assay was conducted following the same setup described above, except that T. harzianum was challenged against itself. Confrontation plates were incubated in the dark at 28°C and the mycelia were harvested at different growth stages: prior to fungal contact, at contact, and after contact (overlapping mycelia). The confrontation assays were conducted in triplicate with RNA extracted for all treatments and replicates. The RNA samples were used for RT-qPCR reactions as described above, with results compared by ANOVA coupled with the Dunnet’s test (α = 5%) using GraphPad Prism 5 software, to allow analysis of differences between confrontation assay gene expression patterns.
Availability of supporting data
Sequences have been deposited at the Sequence Read Archive (SRA) of the National Center for Biotechnology under BioProject number PRJNA216008. Raw sequence reads can be found in http://www.ncbi.nlm.nih.gov/sra/?term=PRJNA216008.
These sequence data (reference genome) were produced by the US Department of Energy Joint Genome Institute http://www.jgi.doe.gov/ in collaboration with the user community. The current study was funded by the National Council for Scientific and Technological Development (CNPq) (Process 559680/2009-0). EFN and CJU were supported by a biotechnology research grant (FAPEGO and CNPq). ASS and MHSR were supported by CNPq PhD scholarships.
- Boland GJ, Hall R: Index of plant hosts of Sclerotinia sclerotiorum. Can J Pl Pathol. 1994, 16: 93-108. 10.1080/07060669409500766.View ArticleGoogle Scholar
- Attanayake RN, Carter PA, Jiang D, Del Río-Mendoza L, Chen W: Sclerotinia sclerotiorum populations infecting canola from China and the United States are genetically and phenotypically distinct. Phytopathology. 2013, 103 (7): 750-761. 10.1094/PHYTO-07-12-0159-R.PubMedView ArticleGoogle Scholar
- Sun P, Yang XB: Light, temperature, and moisture effects on Apothecium production of Sclerotinia sclerotiorum. Plant Dis. 2000, 84: 1287-1293. 10.1094/PDIS.2000.84.12.1287.View ArticleGoogle Scholar
- Lopes FA, Steindorff AS, Geraldine AM, Brandão RS, Monteiro VN, Lobo M, Coelho AS, Ulhoa CJ, Silva RN: Biochemical and metabolic profiles of Trichoderma strains isolated from common bean crops in the Brazilian Cerrado, and potential antagonism against Sclerotinia sclerotiorum. Fungal Biol. 2012, 116 (7): 815-824. 10.1016/j.funbio.2012.04.015.PubMedView ArticleGoogle Scholar
- Qualhato TF, Lopes FA, Steindorff AS, Brandão RS, Jesuino RS, Ulhoa CJ: Mycoparasitism studies of Trichoderma species against three phytopathogenic fungi: evaluation of antagonism and hydrolytic enzyme production. Biotechnol Lett. 2013, 35 (9): 1461-1468. 10.1007/s10529-013-1225-3.PubMedView ArticleGoogle Scholar
- Geraldine AM, Lopes FAC, Carvalho DDC, Barbosa ET, Rodrigues AR, Brandão RS, Ulhoa CJ, Lobo-Junior M: Cell wall-degrading enzymes and parasitism of sclerotia are key factors on field biocontrol of white mold by Trichoderma spp. Biol Control. 2013, 67 (3): 308-316. 10.1016/j.biocontrol.2013.09.013.View ArticleGoogle Scholar
- Lorito M, Woo SL, Harman GE, Monte E: Translational research on Trichoderma: from Omics to the field. Ann Rev Phytopathol. 2010, 48: 395-417. 10.1146/annurev-phyto-073009-114314.View ArticleGoogle Scholar
- Hermosa R, Viterbo A, Chet I, Monte E: Plant-beneficial effects of Trichoderma and of its genes. Microbiology. 2012, 158: 17-25. 10.1099/mic.0.052274-0.PubMedView ArticleGoogle Scholar
- Druzhinina IS, Seidl-Seiboth V, Herrera-Estrella A, Horwitz BA, Kenerley CM, Monte E, Mukherjee PK, Zeilinger S, Grigoriev IV, Kubicek CP: Trichoderma: the genomics of opportunistic success. Nat Rev Microbiol. 2011, 16;9 (10): 749-759.View ArticleGoogle Scholar
- Vizcaíno JA, Redondo J, Suarez MB, Cardoza RE, Hermosa R, Gonzales FJ, Rey M, Monte E: Generation, annotation and analysis of ESTs from four different Trichoderma strains grown under conditions related to biocontrol. Appl Microbiol Biotechnol. 2007, 75: 853-862. 10.1007/s00253-007-0885-0.PubMedView ArticleGoogle Scholar
- Seidl V, Song L, Lindquist E, Gruber S, Koptchinskiy A, Zeilinger S, Schmoll M, Martínez M, Sun J, Grigoriev I, Herrera-Estrella A, Baker SE, Kubicek CP: Transcriptomic response of the mycoparasitic fungus Trichodermaatroviride to the presence of a fungal prey. BMC Genomics. 2009, 10: 567-10.1186/1471-2164-10-567.PubMed CentralPubMedView ArticleGoogle Scholar
- Steindorff AS, Silva RN, Coelho ASG, Noronha EF, Ulhoa CJ: Trichoderma harzianum expressed sequence tags for identification of genes with putative roles in mycoparasitism against F. solani. Biol Control. 2012, 61 (2): 134-140. 10.1016/j.biocontrol.2012.01.014.View ArticleGoogle Scholar
- Vieira PM, Coelho AS, Steindorff AS, De Siqueira SJ, Silva Rdo N, Ulhoa CJ: Identification of differentially expressed genes from Trichoderma harzianum during growth on cell wall of Fusarium solani as a tool for biotechnological application. BMC Genomics. 2013, 14: 177-10.1186/1471-2164-14-177.PubMed CentralPubMedView ArticleGoogle Scholar
- Monteiro VN, Silva RN, Steindorff AS, Costa FT, Noronha EF, Ricart CAO, Sousa MV, Vainstein MH, Ulhoa CJ: New insights in Trichoderma harzianum antagonism of fungal plant pathogens by secreted protein analysis. Curr Microb. 2010, 61: 298-305. 10.1007/s00284-010-9611-8.View ArticleGoogle Scholar
- Rubio MB, Domínguez S, Monte E, Hermosa R: Comparative study of Trichoderma gene expression in interactions with tomato plants using high-density oligonucleotide microarrays. Microbiology. 2012, 158 (Pt 1): 119-128.PubMedView ArticleGoogle Scholar
- Reithner B, Ibarra-Laclette E, Mach RL, Herrera-Estrella A: Identification of mycoparasitism-related genes in Trichoderma atroviride. Appl Environ Microbiol. 2011, 77 (13): 4361-4370. 10.1128/AEM.00129-11.PubMed CentralPubMedView ArticleGoogle Scholar
- Samolski I, de Luis A, Vizcaíno JA, Monte E, Suárez MB: Gene expression analysis of the biocontrol fungus Trichodermaharzianum in the presence of tomato plants, chitin, or glucose using a high-density oligonucleotide microarray. BMC Microbiol. 2009, 9: 217-10.1186/1471-2180-9-217.PubMed CentralPubMedView ArticleGoogle Scholar
- Wang WC, Lin FM, Chang WC, Lin KY, Huang HD, Lin NS: miRExpress: analyzing high-throughput sequencing data for profiling microRNA expression. BMC Bioinform. 2009, 10: 328-10.1186/1471-2105-10-328.View ArticleGoogle Scholar
- Marioni JC, Mason CE, Mane SM, Stephens M, Gilad Y: RNA-seq: an assessment of technical reproducibility and comparison with gene expression arrays. Genome Res. 2008, 18 (9): 1509-1517. 10.1101/gr.079558.108.PubMed CentralPubMedView ArticleGoogle Scholar
- Ozsolak F, Milos PM: RNA sequencing: advances, challenges and opportunities. Nat Rev Genet. 2011, 12 (2): 87-98. 10.1038/nrg2934.PubMed CentralPubMedView ArticleGoogle Scholar
- Martin JA, Wang Z: Next-generation transcriptome assembly. Nat Rev Genet. 2011, 12 (10): 671-682. 10.1038/nrg3068.PubMedView ArticleGoogle Scholar
- Grigoriev IV, Nordberg H, Shabalov I, Aerts A, Cantor M, Goodstein D, Kuo A, Minovitsky S, Nikitin R, Ohm RA, Otillar R, Poliakov A, Ratnere I, Riley R, Smirnova T, Rokhsar D, Dubchak I: The genome portal of the Department of energy joint genome Institute. Nucleic Acids Res. 2012, 40: 26-32.View ArticleGoogle Scholar
- Morin RD, O’Connor MD, Griffith M, Kuchenbauer F, Delaney A, Prabhu AL, Zhao Y, McDonald H, Zeng T, Hirst M, Eaves CJ, Marra MA: Application of massively parallel sequencing to microRNA profiling and discovery in human embryonic stem cells. Genome Res. 2008, 18: 610-621. 10.1101/gr.7179508.PubMed CentralPubMedView ArticleGoogle Scholar
- Druzhinina IS, Kubicek CP, Komoń-Zelazowska M, Mulaw TB, Bissett J: The Trichoderma harzianum demon: complex speciation history resulting in coexistence of hypothetical biological species, recent agamospecies and numerous relict lineages. BMC Evol Biol. 2010, 10: 94-10.1186/1471-2148-10-94.PubMed CentralPubMedView ArticleGoogle Scholar
- Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G: Gene ontology: tool for the unification of biology: the gene ontology consortium. Nat Genet. 2000, 25 (1): 25-29. 10.1038/75556.PubMed CentralPubMedView ArticleGoogle Scholar
- Conesa A, Götz S, Garcia-Gomez JM, Terol J, Talon M, Robles M: Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005, 21: 3674-3676. 10.1093/bioinformatics/bti610.PubMedView ArticleGoogle Scholar
- Free SJ: Fungal cell wall organization and biosynthesis. Adv Gene. 2013, 81: 33-82.View ArticleGoogle Scholar
- Kapteyn JC, Montijn RC, Vink E, de la Cruz J, Llobell A, Douwes JE, Shimoi H, Lipke PN, Klis FM: Retention of Saccharomyces cerevisiae cell wall proteins through a phosphodiester-linked beta-1,3-/beta-1,6-glucan heteropolymer. Glycobiology. 1996, 6 (3): 337-345. 10.1093/glycob/6.3.337.PubMedView ArticleGoogle Scholar
- Elad Y, Kapat A: The role of Trichoderma harzianum protease in the biocontrol of Botrytis cinerea. Eur J Plant Pathol. 1999, 105: 177-189. 10.1023/A:1008753629207.View ArticleGoogle Scholar
- Druzhinina IS, Shelest E, Kubicek CP: Novel traits of Trichoderma predicted through the analysis of its secretome. FEMS Microbiol Lett. 2012, 337 (1): 1-9. 10.1111/j.1574-6968.2012.02665.x.PubMed CentralPubMedView ArticleGoogle Scholar
- Atanasova L, Crom SL, Gruber S, Coulpier F, Seidl-Seiboth V, Kubicek CP, Druzhinina IS: Comparative transcriptomics reveals different strategies of Trichoderma mycoparasitism. BMC Genomics. 2013, 14: 121-10.1186/1471-2164-14-121.PubMed CentralPubMedView ArticleGoogle Scholar
- Wang Y, Guo EW, Yu WG, Han F: Purification and characterization of a new alginate lyase from a marine bacterium Vibrio sp. Biotech Let. 2013, 35 (5): 703-708. 10.1007/s10529-012-1134-x.View ArticleGoogle Scholar
- Skamnaki VT, Peumans WJ, Kantsadi AL, Cubeta MA, Plas K, Pakala S, Zographos SE, Smagghe G, Nierman WC, Van Damme EJ, Leonidas DD: Structural analysis of the Rhizoctonia solani agglutinin reveals a domain-swapping dimeric assembly. FEBS J. 2013, 280 (8): 1750-1763. 10.1111/febs.12190.PubMedView ArticleGoogle Scholar
- Samolski I, Rinco AM, Pinzo LM, Viterbo A, Monte E: The qid74 gene from Trichoderma harzianum has a role in root architecture and plant biofertilization. Microbiology. 2012, 158: 129-138. 10.1099/mic.0.053140-0.PubMedView ArticleGoogle Scholar
- Kubicek CP, Baker S, Gamauf C, Kenerley CM, Druzhinina IS: Purifying selection and birth-and-death evolution in the class II hydrophobin gene families of the ascomycete Trichoderma/Hypocrea. BMC Evol Biol. 2008, 8: 4-10.1186/1471-2148-8-4.PubMed CentralPubMedView ArticleGoogle Scholar
- Seidl-Seiboth V, Gruber S, Sezerman U, Schwecke T, Albayrak A, Neuhof T, Von Öhren H, Baker SE, Kubicek CP: Novel hydrophobins from Trichoderma define a new hydrophobin subclass: protein properties, evolution, regulation, and processing. J Mol Evol. 2011, 72: 339-351. 10.1007/s00239-011-9438-3.PubMedView ArticleGoogle Scholar
- MacIntosh GC, Bariola PA, Newbigin E, Green PJ: Characterization of Rny1, the saccharomyces cerevisiae member of the T2 RNase family of RNases: unexpected functions for ancient enzymes?. Proc Natl Acad Sci USA. 2001, 98 (3): 1018-1023. 10.1073/pnas.98.3.1018.PubMed CentralPubMedView ArticleGoogle Scholar
- Pao SS, Paulsen IT, Saier MH: Major facilitator superfamily. Microbiol Mol Biol R. 1998, 62: 1-34.Google Scholar
- Lohse M, Bolger AM, Nagel A, Fernie AR, Lunn JE, Stitt M, Usadel B: RobiNA: a user-friendly, integrated software solution for RNA-Seq-based transcriptomics. Nucleic Acids Res. 2012, 40: W622-W627. 10.1093/nar/gks540.PubMed CentralPubMedView ArticleGoogle Scholar
- Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL: TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013, 14: R36-10.1186/gb-2013-14-4-r36.PubMed CentralPubMedView ArticleGoogle Scholar
- Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T: Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13: 2498-2504. 10.1101/gr.1239303.PubMed CentralPubMedView ArticleGoogle Scholar
- Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, Salzberg SL, Wold BJ, Pachter L: Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010, 28: 511-515. 10.1038/nbt.1621.PubMed CentralPubMedView ArticleGoogle Scholar
- Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCt method. Methods. 2001, 25 (4): 402-408. 10.1006/meth.2001.1262.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.