Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Identification of a novel fungus, Trichoderma asperellum GDFS1009, and comprehensive evaluation of its biocontrol efficacy

  • Qiong Wu ,

    ‡ These authors are co-first authors on this work.

    Affiliations Department of Environment and Resource, School of Agriculture and Biology, Shanghai Jiao Tong University, Shanghai, China, State Key Laboratory of Microbial Metabolism, Shanghai Jiao Tong University, Shanghai, China, The Key laboratory of Urban (South) Agriculture, Ministry of Agriculture, Shanghai, China

  • Ruiyan Sun ,

    ‡ These authors are co-first authors on this work.

    Affiliations Department of Environment and Resource, School of Agriculture and Biology, Shanghai Jiao Tong University, Shanghai, China, State Key Laboratory of Microbial Metabolism, Shanghai Jiao Tong University, Shanghai, China, The Key laboratory of Urban (South) Agriculture, Ministry of Agriculture, Shanghai, China

  • Mi Ni ,

    ‡ These authors are co-first authors on this work.

    Affiliation State Key Laboratory of Plant Genomics and National Center for Plant Gene Research (Beijing), Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, Beijing, China

  • Jia Yu,

    Affiliations Department of Environment and Resource, School of Agriculture and Biology, Shanghai Jiao Tong University, Shanghai, China, State Key Laboratory of Microbial Metabolism, Shanghai Jiao Tong University, Shanghai, China, The Key laboratory of Urban (South) Agriculture, Ministry of Agriculture, Shanghai, China

  • Yaqian Li,

    Affiliations Department of Environment and Resource, School of Agriculture and Biology, Shanghai Jiao Tong University, Shanghai, China, State Key Laboratory of Microbial Metabolism, Shanghai Jiao Tong University, Shanghai, China, The Key laboratory of Urban (South) Agriculture, Ministry of Agriculture, Shanghai, China

  • Chuanjin Yu,

    Affiliations Department of Environment and Resource, School of Agriculture and Biology, Shanghai Jiao Tong University, Shanghai, China, State Key Laboratory of Microbial Metabolism, Shanghai Jiao Tong University, Shanghai, China, The Key laboratory of Urban (South) Agriculture, Ministry of Agriculture, Shanghai, China

  • Kai Dou,

    Affiliations Department of Environment and Resource, School of Agriculture and Biology, Shanghai Jiao Tong University, Shanghai, China, State Key Laboratory of Microbial Metabolism, Shanghai Jiao Tong University, Shanghai, China, The Key laboratory of Urban (South) Agriculture, Ministry of Agriculture, Shanghai, China

  • Jianhong Ren,

    Affiliation Suzhou BioNovoGene Metabolomics Platform, Suzhou, China

  • Jie Chen

    jiechen59@sjtu.edu.cn

    Affiliations Department of Environment and Resource, School of Agriculture and Biology, Shanghai Jiao Tong University, Shanghai, China, State Key Laboratory of Microbial Metabolism, Shanghai Jiao Tong University, Shanghai, China, The Key laboratory of Urban (South) Agriculture, Ministry of Agriculture, Shanghai, China

Abstract

Due to its efficient broad-spectrum antimicrobial activity, Trichoderma has been established as an internationally recognized biocontrol fungus. In this study, we found and identified a novel strain of Trichoderma asperellum, named GDFS1009. The mycelium of T. asperellum GDFS1009 exhibits a high growth rate, high sporulation capacity, and strong inhibitory effects against pathogens that cause cucumber fusarium wilt and corn stalk rot. T. asperellum GDFS1009 secretes chitinase, glucanase, and protease, which can degrade the cell walls of fungi and contribute to mycoparasitism. The secreted xylanases are good candidates for inducing plant resistance and enhancing plant immunity against pathogens. RNA sequencing (RNA-seq) and gas chromatography-mass spectrometry (GC-MS) showed that T. asperellum GDFS1009 produces primary metabolites that are precursors of antimicrobial compounds; it also produces a variety of antimicrobial secondary metabolites, including polyketides and alkanes. In addition, this study speculated the presence of six antimicrobial peptides via ultra-performance liquid chromatography quadrupole time of flight mass spectrometry (UPLC-QTOF-MS/MS). Future studies should focus on these antimicrobial metabolites for facilitating widespread application in the field of agricultural bio-control.

Introduction

Trichoderma spp. exhibits antagonistic effects against at least 18 genera and 29 species of pathogenic fungi, as well as a variety of pathogenic bacteria. The biocontrol mechanisms of Trichoderma spp. primarily include competition and mycoparasitism, followed by the stimulation of plant resistance and immunity [1].

Trichoderma spp. is highly adaptive to the environment, and their growth rates are generally much faster than those of plant pathogens. Trichoderma spp. could compete with nearby pathogens for limited space and nutrients, thereby inhibiting their growth. Competition mechanism plays an important role in the inhibition of soil-borne diseases and leaf pathogens [2, 3].

The most significant antagonistic mechanism of Trichoderma spp. against pathogens is mycoparasitism. After recognizing the pathogen, the Trichoderma mycelium grows alongside the pathogen mycelium in a spiral fashion, resulting in dissolution and death of the pathogen. This process is accompanied by the secretion of cell wall-degrading enzymes (CWDEs), such as chitinases, glucanase, and proteases, which penetrate the pathogen mycelium, absorbing its nutrients and ultimately dissolving the pathogen [46]. As is well known, the chitinase protein CHIT33 from Trichoderma harzianum CECT2413 exerts a significant inhibitory effect on Fusarium [7]. Overexpression of the endochitinase gene (chit36) from T. harzianum Rifai TM results in complete inhibition of Botrytis cinerea spore germination [8]. The alkaline protease gene prb1, isolated from T. harzianum, is associated with mycoparasitism [9]. Exo and endo-β-1, 3-glucanase are also involved in the mycoparasitism of Trichoderma [10, 11].

When interacting with plants, Trichoderma spp. produce elicitors, which induce systemic resistance in plants and enhance their immunity against pathogens. The xylanase protein XYN2 secreted by T. reesei QM6a causes tissue necrosis and induces systemic resistance in tobacco by stimulating the ethylene signaling pathway and regulating the overexpression of downstream resistance-associated proteins such as PR2b and ERF2 [12]. SM1, a small 14-kDa protein secreted by T. viride Gv 29–8, induces immune resistance in plants by stimulating the ethylene-jasmonic acid signaling pathway [13, 14]. In addition, several other elicitors have been discovered and found to induce plant systemic resistance [15].

Trichoderma spp. induces pathogen inhibition by secreting secondary metabolites. Different Trichoderma species secrete different substances, including isonitrile, diketopiperazines, sesquiterpenes, stemids, polyketides, alkylpyrones, and peptaibols [16]. It has been shown that 6-pentyl-α-pyrone, produced by T. harzianum, efficiently inhibits the growth of Rhizoctonia solani on rice [17]. At present, more than 300 compounds have been identified as antimicrobial peptaibols, including 11 from T. longibrachiatum and one peptaibol polymer with 20 peptaibol residues from T. brevicompactum [18, 19].

Previously, research of Trichoderma spp. mainly focused on only one or a few biocontrol factors in one strain. Our study reveals the resistance-related molecular mechanisms of T. asperellum GDFS1009, including almost all of the biocontrol factors related to mycoparasitism, induced resistance, and antibiosis, using genomics, transcriptomics, and metabolomics approaches, as well as a series of physiological and biochemical analyses. These data provide a good foundation for continued researches into T. asperellum GDFS1009 for facilitating widespread application in the field of agricultural bio-control.

Materials and methods

Isolation and morphological and molecular identification of a novel strain

Irrigated soil sample was collected in Foshan, Guangdong, China via the diagonal method. No specific permissions were required for the location and the related activities about the soil sample, because the sample was collected on private land and we confirm that the owner of the land gave permission to conduct the study on this site. Furthermore, we confirm that the field studies did not involve endangered or protected species.

Via dilution plating, rose bengal-agar culture plates containing 50 μg/mL streptomycin and 50 μg/mL chloramphenicol were used to isolate Trichoderma spp., among which a novel strain T. asperellum, named GDFS1009, was discovered by morphological observation. Then the classification position was confirmed by ITS sequencing.

T. asperellum GDFS1009 was cultured in potato-dextrose (PD) broth in a rotary shaker for 2 d at 28°C and 180 rpm until its mycelium was well developed. The supernatant was removed by vacuum filtration, the mycelium was collected as the pellet, and the genome was extracted using the hexadecyl trimethyl ammonium bromide (CTAB) method [20]. Using universal ITS primers (ITS1-F and ITS4) [21], ITS sequences were amplified by PCR using the T. asperellum GDFS1009 genome as the template. The PCR product was ligated into the pMD19-T vector (Takara, Shiga, Japan) and was then sequenced by the Shanghai Biosune Biotechnology Co. (Shanghai, China). The sequencing results were aligned with the ITS sequences of known strains using NCBI BLAST (http://www.ncbi.nlm.nih.gov/BLAST). MEGA 4.0 software was used to build a phylogenetic tree for ITS [22].

Detection of biocontrol efficacy

Discs were obtained from 4-day-old cultures of pathogens (Fusarium oxysporum f. sp. cucumerinum Owen and Fusarium graminearum) and T. asperellum GDFS1009 using sterile 7 mm hole-punchers. Each disc was inverted and placed on one side of a 9-cm Petri dish with PDA medium, with 4 cm between pathogen and T. asperellum GDFS1009 discs. The control was a dish containing only the pathogen disc. Dishes were incubated at 28°C for 5 d [23].

Then spore suspension (106 cfu/mL) of T. asperellum GDFS1009 was cultured in PD broth in a rotary shaker at 28°C and 180 rpm for 6 d. A total of 50 mL of filtered fermentation liquor from T. asperellum GDFS1009 or ddH2O (control) were mixed into 100mL PDA medium, and the plate was made in 9-cm Petri dish. 7 mm discs of 4-day-old cultures of pathogens (Fusarium oxysporum f. sp. cucumerinum Owen and Fusarium graminearum) were inverted and placed on the center of each plate, then were incubated at 28°C for 4 d. Additionally, related micro well dilution assay of pathogens inhibition was also done with filtered fermentation liquor from T. asperellum GDFS1009.

Cucumber (Jinyan no. 4) plants were grown to the 3–5-leaf stage in soil under a 10-h light/14-h dark cycle in a 28°C phytotron. A 20-mL spore suspension of F. oxysporum f. sp. cucumerinum Owen (106 cfu/mL) was irrigated onto the roots of cucumber seedlings, followed by a 20-mL spore suspension of T. asperellum GDFS1009 (107 cfu/mL) in the treatment group. The control group was irrigated with sterile water instead of T. asperellum GDFS1009. The inoculated plants were placed in the 28°C phytotron for at least 7 d, and biocontrol effects were observed and analyzed [24].

Assay of enzymes related to mycoparasitism and induced resistance

T. asperellum GDFS1009 was cultured in PD broth in a rotary shaker at 28°C and 180 rpm for 2 d until the mycelium was well developed. The supernatant was removed by vacuum filtration, and 1 g of mycelium pellet was transferred to xylanase induction medium includes 2% corncob powder, 0.4% NH4NO3, 0.05% K2HPO4, 0.05% MgSO4·7H2O, 0.05% KCl, and 0.001% FeSO4·7H2O. The activities of xylanase were measured using 3, 5-dinitrosalicylic acid (DNS). At 50°C and pH 5.0, 1 U of xylanase activity was defined as the amount of enzyme required to release 1 μmol xylose per minute. Xylose was used as standard for xylanase analysis [25].

Then spore suspension (104 cfu/mL) of T. asperellum GDFS1009 was cultured in protease induction medium (1% casein, 0.05% K2HPO4, 0.05% MgSO4·7H2O, 0.05% KCl, and 0.001% FeSO4·7H2O) in a rotary shaker at 28°C and 180 rpm for 5 d. The activities of protease were measured every day using foline-phenol reagent. At 50°C and pH 7.5, 1 U of protease activity was defined as the amount of enzyme required to release 1 μg L-tyrosine per minute. L-tyrosine was used as standard for protease analysis. Enzyme activity was statistically analyzed via SPSS software. Each measurement was performed in triplicate. Furthermore, discs were obtained from 2-day-old cultures of T. asperellum GDFS1009 using sterile 3 mm hole-punchers. Each disc was inverted and placed on the center of a 9-cm Petri dish with 1.5% skim milk and 2.5% agar. After incubation at 28°C for 36h, proteolytic circle was observed [2627].

Tobacco (Nicotiana benthamiana) plants and cucumber plants were grown to the 6–8-leaf stage in soil under a 10-h light/14-h dark cycle in a 24°C phytotron. A total of 20 μL of filtered xylanase preparation from T. asperellum GDFS1009 or blank broth (control) were injected into tobacco or cucumber leaves using 1 mL syringes. Injected leaves in each group were of the same size and position in the same tobacco or cucumber plant. The inoculated plant was placed in darkness in a 24°C phytotron, and tissue necrosis was observed every 12 h. For oxidative burst assay, injected leaves were stained in DAB solution (pH 3.8) for 8 h in darkness, boiled in ethanol, glycerin, and acetic acid (1:1:1, v/v) for 10–15 min, and washed twice with 60% ethanol [28].

Global analysis of gene expression based on RNA-seq

The spore suspension of T. asperellum GDFS1009 was transferred into PD medium with a final spore concentration of 106 cfu/mL. After culturing at 28°C and 180 rpm for 24 h, the mycelium was collected by vacuum filtration. Single-end (1×50) Solexa sequencing was performed, producing 20 million reads per sample. Additionally, 48-h samples were similarly sequenced, producing 10 million reads.

The reference genome for transcriptome analysis was acquired from http://genome.jgi.doe.gov/Trias1/Trias1.download.html. Conventional analysis was performed on the transcriptome data, including data pre-processing, genomic mapping, gene expression analysis, transcript expression analysis, alternative splicing analysis, analysis of differentially expressed genes, and GO/KEGG enrichment analysis of differentially expressed genes and non-differentially expressed genes [2931]. On the basis of these results, genes related to mycoparasitism, induced resistance, and antibiosis in other Trichoderma strains were downloaded, and then local BLAST was used to compare them to the public T. asperellum genomic data mentioned above. Then all homologous genes were assembled and matched to the RNA-seq data of T. asperellum GDFS1009 to determine their expression levels. When the FPKM value exceeded 100, the expression level was defined as “high”, while an FPKM value of less than 10 was defined as “low”, and a value less than 1 was considered “ultra-low” expression. A heat map was constructed in the R language [32].

Assay of primary metabolites based on GC-MS

The mycelia of T. asperellum GDFS1009 were inoculated into a PDA plate at 28°C for 4 d. Then, spores were collected in sterile water and diluted to 108 cfu/mL. Spores (1 mL) were added to 100 mL PD broth and incubated on a 180-rpm rotary shaker at 28°C for 1 d. Next, 50 mL of this sample was dropped directly into 200 mL 60% (v/v) methanol solution (taken from a -40°C deep freezer). The mixture was allowed to cool for 10 min to allow the temperature to return to -40°C. The mixture was then centrifuged at 5000 rpm for 5 min at -10°C, and the mycelium pellet was collected. Following this procedure, the mycelium pellet was freeze-dried at -40°C.

Subsequently, 500 μL 100% methanol at -80°C was added to 100 mg of sample, vortexed for 30 s, and placed in liquid nitrogen for 15 min. After thawing at 4°C for 15 min and vortexing for 1 min, the mixture was centrifuged at 13000 rpm for 5 min at 4°C, and the resulting supernatant was transferred to tube A. Milli-Q water (250 mL) was used to resuspend the pellet after centrifugation. The resuspension was then vortexed for 30 s and placed in liquid nitrogen for 15 min. After thawing at 4°C for 15 min, glass beads (G8772-100G, Sigma, St. Louis, MO, USA) were added and vortexed for 10 min, and then the mixture was centrifuged at 13000 rpm for 5 min at 4°C. The supernatant was also transferred to tube A. After drying with nitrogen, 60 μL methoxyl amine hydrochloride (15 mg/mL) was added to tube A and vortexed for 1 min, and then the mixture was allowed to react at 25°C for 12 h. After adding 60 μL MSTFA with 1% TMCS and allowing the reaction occur at room temperature for 1 h, the mixture was centrifuged at 13000 rpm for 5 min at 4°C.

The supernatant was then used for gas chromatography-mass spectrometry (GC-MS; Agilent 7890A/5975C, Agilent, Santa Clara, CA, USA). The derivatized samples (1 μL) were injected at a 1:20 split ratio onto a HP-5MS column (5% phenyl methyl silox; 30 m × 250 μm i.d., 0.25-μm; Agilent J&W Scientific, Folsom, CA, USA). During GC, helium was supplied at a constant rate of 1 mL/min. The temperature of the injection port was 280°C, the ion source temperature was 250°C, and the interface temperature was 150°C. The temperature program started at 40°C for 5 min and increased at a rate of 10°C/min until it reached 300°C, where it was held constant for 5 min. MS was determined using the full-scan method with a range from 35 to 780 m/z. Identification of metabolites was performed in six biological replicates. Bioinformatics analyses included data pre-processing (XCMS, www.bioconductor.org) and compound identification (NIST 2008, Wiley 9) [3335].

Assay of secondary metabolites based on GC-MS

The mycelia of T. asperellum GDFS1009 were inoculated into PDA plate at 28°C for 6 d. After conidia fully covered the plate, they were carefully scraped from the plate surface with a sterile spoon. Spores (6 g) were weighed and extracted with 50 volumes of dichloromethane at 4°C for 2 d. The extract was exposed to 5% activated carbon under oscillation for 2 h and filtered with four layers of sterilized gauze. The supernatant was washed twice with an equal volume of 3% sodium carbonate solution. The resulting dichloromethane fraction was dehydrated with anhydrous sodium sulfate and then filtered with fast filter paper. Finally, the dehydrated liquid was vacuum-evaporated at 40°C to obtain 1 mL of viscous crude extract [36]. GC-MS analysis was carried out using an Auto System XL GC/Turbo Mass MSUSA (Perkin Elmer, USA, DB-5MS column: 30 m×0.25 mm, 0.25 μm). The GC conditions were as follows: after injection of 1 μL, the initial column temperature was held at 50°C for 5 min and then raised to 300°C at 5°C/min, and the vaporizing chamber temperature was then maintained at 300°C. The GC was operated with a split ratio of 10:1 using helium as the carrier gas. The electron-ionization (EI) source was held at 230°C and the ionization voltage was 70 eV. The quadrupole mass filter was run at an effective scanning range of 290–500 amu. Metabolites were qualitatively analyzed by GC-MS. The peak area of each metabolite, acquired in SIM mode (selective ion monitoring; signal/noise > 5.0), was normalized to the peak area of 2-heptone prior to further data processing.

Then spore suspension (106 cfu/mL) of T. asperellum GDFS1009 was cultured in PD broth in a rotary shaker at 28°C and 180 rpm for 6 d. A total of 20 mL of filtered fermentation liquor from T. asperellum GDFS1009 (PD broth as control) was enriched to 5 mL, extracted three times with 5mL dichloromethane, and was washed twice with an equal volume of 3% sodium carbonate solution. The resulting dichloromethane fraction was vacuum-evaporated at 40°C to obtain 1 mL of viscous crude extract, and then was dehydrated with anhydrous sodium sulfate. Gas chromatography was performed on a HP-5MS capillary column (5% phenyl/95% methylpolysiloxane (30 m × 250 μm i.d., 0.25 μm film thickness, Agilent J & W Scientific, Folsom, CA, USA) to separate the derivatives at a constant flow of 1 mL/min helium. 1 μL of sample was injected in split mode in a 20:1 split ratio by the auto-sampler. Injection temperature was 280°C, the interface set to 150°C and the ion source adjusted to 230°C. The programs of temperature-rise was followed by initial temperature of 80°C for 5 min, 20°C/min rate up to 300°C and staying at 300°C for 6 min. Mass spectrometry was determined by full-scan method with range from 35 to 500 (m/z). Identification of metabolites using the Automatic Mass Spectral Deconvolution and Identification System (AMIDS) was searched against commercial available databases such as National Institute of Standards and Technology (NIST) and Wiley libraries [3335].

Assay of antimicrobial peptides based on UPLC-QTOF-MS/MS

Mycelia of T. asperellum GDFS1009 were inoculated onto a PDA plate at 28°C for 4 d. Then, spores were collected with sterile water to 108 cfu/mL. Spores (1 mL) were added to 100 mL antimicrobial peptide media [0.5% glucose, 0.08% KH2PO4, 0.07% KNO3, 0.02% Ca(H2PO4)2, 0.05% MgSO4.7H2O, 0.001% MnSO4.5H2O, 0.0005% CuSO4.5H2O, and 0.0001% FeSO4.7H2O], and the mixture was incubated on a 180-rpm rotary shaker at 28°C for 20 d. The broth was filtered using two layers of qualitative filter paper, extracted with a filtrate:butanol ratio of 3:1, and the lower layer liquid was transferred to a vacuum distillation flask. After distillation, the dry powder was completely dissolved in 80 mL of methanol:dichloromethane at a ratio of 1:1, and then the dissolved liquid was filtered with a 0.45-μm PTFE membrane. After evaporating to dryness in a vacuum distillation flask, the solution was completely dissolved in 5 mL of dichloromethane:methanol at a ratio of 85:15 in a column that had been previously washed with 85:15 dichloromethane:methanol to remove acetone. This was followed by the addition of 100 mL of 85:15 dichloromethane:methanol. The collected filtrate was again placed in a vacuum distillation flask and evaporated to dryness. The dry powder was completely dissolved in 5 mL of methanol:H2O at a ratio of 85:15, transferred to a 1.5-mL centrifuge tube, and centrifuged at 13000 rpm for 10 min. After centrifugation, the supernatant was stored in a 2.0-mL glass vial at 4°C [37].

UPLC-QTOF-MS/MS was carried out using an ACQUITY UPLC & SCIEX SelexION Triple QuadTM 5500 System (Waters, USA) equipped with an Acquity BEH C18 column (2.1 × 100 mm, 1.7 μm). The UPLC conditions were as follows: the mobile phase consisted of water (eluent A) and 0.1% acetonitrile (eluent B), the column temperature was 45°C, the injection volume was 3.00 μL, and the running time was 26.50 min. MS parameters were as follows: the spectrometer was operated in the positive ion mode using an electrospray-ionization source, the capillary voltage was 3.0 kV, the sampling cone voltage was 4.0 V, the extraction cone voltage was 3.0 V, the ion source temperature was 105°C, the atomization air temperature was 350°C, the atomizing gas flow was 700.0 L/h, the collision voltage was 3.0 eV, the scan time was 0.300 s, the scan range was 200–2200 m/z, and LockSpray online correction was performed using 200 ng/mL leucine enkephalin (556.2771 m/z). MassLynx software was used for the alignment, fitting, identification, and analysis of different antimicrobial peptides.

Results

Colony morphology, ITS identification, and phylogenetic tree

Morphological observations showed that after five days of culture, the front of the novel isolate’s colony was dark green (spore color) and the back was white (mycelium color). The morphology of the mycelium was coarse, and dark green spores began to form in the center of the colony at high rates (Fig 1A). Based on preliminary morphological observations, the strain was confirmed to be Trichoderma.

thumbnail
Fig 1. Morphological, ITS and phylogeny analysis of T. asperellum GDFS1009.

(A) Photograph of colonial morphology on PDA medium; (B) ITS amplification; M, Marker DL2000 (Takara, Japan); 1, Band of ITS; (C) Phylogenetic tree based on nrDNA ITS sequence of T. asperellum GDFS1009 (bold) and its closely related species and outgroup retrieved from the literature.

https://doi.org/10.1371/journal.pone.0179957.g001

The result of the PCR amplification of the ITS (KF367522) is showed in Fig 1B, with a product size of 642 bp. NCBI BLAST followed by phylogenetic analysis of ITS sequences revealed that this novel strain shares the highest homology with T. asperellum G, followed by T. harzianum, T. longibrachiatum, T. hamatum, T. reesei, and T. virens, with Protocrea pallida used as an outgroup (Fig 1C). Together with the morphological identification, we determined the strain to be T. asperellum and named it GDFS1009.

Biocontrol effects under in vitro and greenhouse conditions

In vitro dual-culture analysis showed that the inhibition rate of T. asperellum GDFS1009 against F. oxysporum f. sp. cucumerinum Owen was 80.82% (Fig 2A). The inhibition rate of filtered fermentation liquor from T. asperellum GDFS1009 against F. oxysporum f. sp. cucumerinum Owen was 67.59% (Fig 2B). Micro well dilution assay on 2 days showed consistent results, up to 76.28% (Fig 2C). Under greenhouse conditions, T. asperellum GDFS1009 exhibited 86.34% growth inhibition (Fig 2D). Based on the in vitro antagonism assay and greenhouse study, T. asperellum GDFS1009 exhibits significant biocontrol effects on F. oxysporum f. sp. cucumerinum Owen. Additionally, except in the case of a low incidence of fusarium wilt, after treatment with Trichoderma spore suspension, cucumber seedling growth was more robust, with bright green leaves and strong roots (Fig 2D).

thumbnail
Fig 2. Control efficacy of T. asperellum GDFS1009 on F. oxysporum f. sp. cucumerinum Owen on cucumber.

(A) Dual-culture assay, a1, Synergistic effect, a2, CK; (B) Resistant-dish assay, b1, Synergistic effect, b2, CK; (C) Micro well dilution assay; (D) Assay in green house, d1, Synergistic effect, d2, CK.

https://doi.org/10.1371/journal.pone.0179957.g002

Comparing with F. oxysporum f. sp. cucumerinum Owen, in vitro dual-culture analysis showed that the inhibition rate of T. asperellum GDFS1009 against F. graminearum was not so high. However, the inhibition rate of filtered fermentation liquor from T. asperellum GDFS1009 against F. graminearum was 100%. Micro well dilution assay on 2 days showed consistent results, nearly 100% (S1 Fig).

Mycoparasitism-related enzyme gene expression and activity

The relative expression levels of mycoparasitism genes were analyzed using RNA-seq analysis and illustrated with a heat map. Local BLAST identified 16 chitinase genes in the T. asperellum GDFS1009 genome and eight protease genes, with two being highly expressed after 24 h and one after 48 h. A total of 11 glucanase genes were discovered, with only one of these being highly expressed after 24 h and 48 h. Most of the non-highly expressed genes mentioned above were expressed at low levels, and a few did not appear to be expressed under the experimental conditions (Fig 3A and S1 Table).

thumbnail
Fig 3. Protease analysis related to mycoparasitism in T. asperellum GDFS1009.

(A) Heatmap of CWDEs expression based on RNA-seq; (B) Proteolytic circle assay; (C) Protease activity.

https://doi.org/10.1371/journal.pone.0179957.g003

After incubation at 28°C for 36h, proteolytic circle in skim milk plate was 1.656 ± 0.030 cm (Fig 3B). Activity of mycoparasitism-related protease was assessed in corresponding induction broth. And protease activity was highest after fermentation for 2 d, reaching 1.418 ± 0.023 U/mL at 50°C and pH 7.5 (Fig 3C). In addition, chitinase and glucanase were also discovered with high level expression in induction broth (Data not shown).

Induced resistance-related gene expression and induction of tissue necrosis and resistance in plants

In the T. asperellum GDFS1009 genome, 12 elicitor genes with different functions were found, including two endopolygalacturonase genes, two Epl protein genes, two hydrophobin genes, one polygalacturonase gene, one swollenin gene, and four xylanase genes. Most of the elicitor genes were expressed at very low levels or not at all, with only three being expressed highly in PD broth after 24 h and 48 h (Fig 4A and S2 Table). The xylanase activity of T. asperellum GDFS1009 was assessed in induction broth, and exhibiting a maximum activity of 11.477 ± 0.613 U/mL at 50°C and pH 5.0 after fermentation for 3 d (Fig 4B).

thumbnail
Fig 4. Induced resistance gene analysis of T. asperellum GDFS1009.

(A) Heatmap of induced resistance gene expression based on RNA-seq; (B) Xylanase acvitity; (C) Photographs of tissue necrosis and DAB staining in tobacco leaves after inoculation for 24 h, c1 and c3, Tobacco leaf injected with xylanase preparation of T. asperellum GDFS1009, c2 and c4, Tobacco leaf injected with supernatant of blank broth (control); (D) Photographs of tissue necrosis and DAB staining in cucumber leaves after inoculation for 24 h, d1 and d3, Cucumber leaf injected with xylanase preparation of T. asperellum GDFS1009, d2 and d4, Cucumber leaf injected with supernatant of blank broth (control).

https://doi.org/10.1371/journal.pone.0179957.g004

At 12 h after inoculation, necrotic lesions were seen on tobacco and cucumber leaves injected with a xylanase solution from T. asperellum GDFS1009, but no necrotic lesions were seen on leaves treated with control broth. This result was most apparent at 24 h. In addition, sites both near and far from the inoculation sites of tobacco and cucumber leaves injected with xylanase solution were stained by DAB, while those in tobacco and cucumber leaves injected with control broth exhibited nearly no staining (Fig 4C and 4D).

Biocontrol-related primary metabolites determined by GC-MS analysis

Using GC-MS analysis, nine were related to biocontrol as fungicides or a precursor/ intermediate of a fungicide, insecticide, or herbicide. Acetamide and diethylamine are fungicide precursors. Ethylene Glycol, ethanolamine, and o-toluic acid are fungicide intermediates. Citric acid is a fungicide. Ethylamine is a precursor of triazine herbicides, and glycine is a precursor of fungicides, insecticides, and herbicides (Table 1).

thumbnail
Table 1. Bio-control related primary metabolites in T. asperellum GDFS1009 based on GC-MS analysis.

https://doi.org/10.1371/journal.pone.0179957.t001

Secondary metabolites determined by GC-MS analysis and gene cluster prediction

In spore of T. asperellum GDFS1009, 68 compounds were obtained via GC-MS analysis, including 2 polyketides, 8 olefins, 24 alkanes, 2 acids, 25 esters, 1 aldehydes, 2 benzenes, 4 alcohols (Fig 5A and S3 Table). In fermentation liquor of T. asperellum GDFS1009, 28 compounds were obtained via GC-MS analysis, including 1 polyketides, 1 olefins, 23 alkanes, 2 acids, 1 nitriles (Fig 5B and S4 Table).

thumbnail
Fig 5. Antibiosis secondary metabolites in T. asperellum GDFS1009 based on GC-MS analysis.

(A) Secondary metabolites in spore of T. asperellum GDFS1009; (B) Secondary metabolites in fermentation liquor of T. asperellum GDFS1009; (C) Some structures of polyketides, alkanes and acids related to pathogen resistance and induced resistance.

https://doi.org/10.1371/journal.pone.0179957.g005

1-(4-bromobutyl)-2-piperidinone, 2,2,6,6-tetramethyl-4-piperidinone, 6-pentyl-2H-pyran-2-one, 2,6,10-trimethyl-tetradecane, 2,6,10-trimethyl-pentadecane, and oxalic acid, butyl 6-ethyloct-3-ylester in T. asperellum GDFS1009 were found, which may be related with pathogen inhibition and induced resistance (Fig 5C).

Local BLAST revealed 16 gene clusters of polyketide synthases (PKSs) in the T. asperellum GDFS1009 genome, among which 12 clusters were fully reducing PKSs and four were non-reducing PKSs. Based on the analysis of RNA-seq data, after 24 h and 48 h in PD broth, almost all of the gene clusters exhibited low levels of expression or were not expressed (Table 2).

thumbnail
Table 2. Gene cluster analysis of polyketide and nonribosomal peptides synthesis in T. asperellum GDFS1009.

https://doi.org/10.1371/journal.pone.0179957.t002

Antimicrobial peptides determined by UPLC-QTOF-MS/MS analysis and gene cluster prediction

Through UPLC-QTOF-MS/MS analysis, a total of six types of antimicrobial peptides were speculated from T. asperellum GDFS1009 with molecular weights ranging from 1092.6450 to 1734.0618. These were TBV (11 aa), trichotoxin_A-50_F (18 aa), hypomurocin_B_IIIa (18 aa), trichotoxin_A-50_I (18 aa), hypomurocin_B_I (18 aa), and hypomurocin_B_II (18 aa; Table 3). These antimicrobial peptides are formed through the catalysis of non-ribosomal peptide synthases (NRPS). Local BLAST revealed 16 NRPS gene clusters in the T. asperellum GDFS1009 genome, with two PKS/NRPS hetero-gene clusters. Based on the analysis of RNA-seq data, after 24 h and 48 h in PD broth, almost all of the gene clusters were lowly expressed or not expressed (Table 2).

thumbnail
Table 3. Peptaibols antimicrobial peptides in T. asperellum GDFS1009 based on UPLC-QTOF-MS/MS analysis.

https://doi.org/10.1371/journal.pone.0179957.t003

Discussion

In recent years, the United States, Canada, Australia, Poland, and many other countries have participated in collecting Trichoderma resources. At present, 141 Trichoderma species have been discovered in the world [3841]. Early studies on Trichoderma mostly focused on aspects of soil microbiology and general mycology. As our understanding of Trichoderma has deepened, research into the mechanism of biocontrol and application has gradually increased.

This study began with the identification and biocontrol assessment of a novel Trichoderma, named T. asperellum GDFS1009, which exhibits strong control effects on F. graminearum and F. oxysporum (Figs 1 and 2; S1 Fig). As is typical of soil-borne diseases, F. graminearum and F. oxysporum seriously harm vegetables and food crops. The biocontrol data from in vitro and greenhouse experiments provide an experimental basis for further investigations into spore granule production for use in agricultural settings.

Previously, research into the resistance mechanisms of Trichoderma spp. have been fairly broad, with most concentrated on a single or a few biocontrol factors in a single strain. Our study reveals the global resistance mechanisms of T. asperellum GDFS1009 at the molecular level, including almost all biocontrol factors related to mycoparasitism, induced resistance, and antibiosis using genomics, transcriptomics, and metabolomics analyses as well as a series of physiological and biochemical analyses.

Chitinases, involved in mycoparasitism, are diverse, inducible secretory proteins. Characteristics of chitinases differ in different Trichoderma species. Trichoderma koningii inhibits the growth of Sclerotium cepivorum in onion root via the production of exo- and endo-chitinases, which digest the cell walls of pathogens [42]. β-1, 3-glucanase, also associated with mycoparasitism, induces strong hydrolysis of pathogen cell walls, thereby inhibiting germination of pathogen spores and digestion of the mycelium [43]. In addition, proteases play an important role in mycoparasitism, as they are necessary for the complete hydrolysis of pathogen cell walls [44]. In our study, 16 chitinase genes, eight protease genes, and 11 glucanase genes were discovered in T. asperellum GDFS1009, only a few of which were highly expressed in PD broth without induction medium (Fig 3 and S1 Table). In future experiments, the specific functions of these enzymes, such their tolerances to heat, cold, or metal ions, could be predicted with bioinformatics analysis or via expression in E. coli or Saccharomyces cerevisiae in order to more fully understand their enzymatic properties. In the induction broths, the three types of enzymes were all detected and exhibited high activities, providing an experimental basis for subsequent development of enzyme preparations. Induction broths should be optimized for greater enzyme activity, and single or mixed enzyme preparations may be applied in agricultural production.

In recent years, it was discovered that Trichoderma and its formulations could elicit pathogen-related gene expression in host plants and induce systematic acquired resistance (SAR). In 1999, Yedidia found that when the Trichoderma mycelium crosses the epidermis layer of cucumber root, it induces a cascade of pathogen-related protein expression, similar to that caused by the chemical mutagen 2, 6-dichlorphenoxyacetic acid [45, 46]. Similarly, after treatment of the cucumber root with T. harzianum T-203, the JA/SA signal pathway was induced, and Pseudomonas lachrymans was inhibited [47, 48]. Numerous studies have shown that Trichoderma produces several types of inducible factors or elicitors, such as sm1 and xylanase. In this study, xylanase activity was also detected in induction broth. This may induce plant resistance due to the induction of tissue necrosis. Additionally, in the T. asperellum GDFS1009 genome, 12 elicitor genes with different functions were found, including two endopolygalacturonase genes, two Epl protein genes, two hydrophobin genes, one polygalacturonase gene, one swollenin gene, and four xylanase genes. The expression levels of most of these elicitor genes were very low or not detected, with only a few being highly expressed in PD broth after 24 h and 48 h (Fig 4 and S2 Table). In future experiments, elicitors which mechanisms are unknown could be expressed in E. coli or Saccharomyces and be purified for further investigation into the interactions with plants and pathogens.

Nine primary metabolites believed to be fungicides or precursors/intermediates of fungicides, insecticides, or herbicides were identified in T. asperellum GDFS1009 (Table 1). Acetamide is widely used as a pesticide. Dibromocyanoacetamide (DBNPA) is a novel potent sanitizer and excellent water treatment agent [49]. Ethylamine can be used as a raw material to produce triazine herbicides [50], while ethylene glycol can be used as an intermediate for pesticide production [51]. Glycine can be used to synthesize fungicides, and research has also shown that copper glycine inhibits B. cinerea and F. oxysporum [52]. Ethanolamine is an intermediate in the synthesis of pesticides. The combined application of Schiff bases of salicylaldehyde ethanolamine with copper acetate or zinc sulfate exhibited significant inhibition of Aspergillus niger and F. oxysporum [53]. As an important organic synthesis intermediate, o-toluic acid is also widely used in the production of pesticides [54]. L-malic acid can be used for pharmaceutical formulations and also for the synthesis of insect repellents [55]. Citric acid strongly inhibits Aeromonas hydrophila, Pseudomonas fluorescens, and other aquatic pathogens [56].

Trichoderma also produces a variety of antimicrobial secondary metabolites, including diketopiperazines, sesquiterpenes, polyketides, olefins, alkanes and so on [5758]. The biological activities of several antimicrobial substances have been identified [5962]. In 1987, Claydon purified alkyl pyrone from T. harzianum [63]. Horace reported that a main mechanism of the inhibition of Rhizoctonia solani by T. harzianum involves the production of 6-pentyl pyrone [17]. It was found that 2,2,6,6-tetramethylpiperidine also has antifungal activity [64]. Furthermore, oxalic acid-induced resistance to Rhizoctonia solani in rice is associated with induction of phenolics, peroxidase and pathogenesis-related proteins [65]. By comparing the secondary metabolites in T. asperellum GDFS1009 to standards and databases, in spore of T. asperellum GDFS1009, 68 compounds were obtained via GC-MS analysis, including 2 polyketides, 8 olefins, 24 alkanes, 2 acids, 25 esters, 1 aldehydes, 2 benzenes, 4 alcohols (Fig 5A and S3 Table). In fermentation liquor of T. asperellum GDFS1009, 28 compounds were obtained via GC-MS analysis, including 1 polyketides, 1 olefins, 23 alkanes, 2 acids, 1 nitriles (Fig 5B and S4 Table). Comparative analysis of spore and fermentation liquor showed that: (1) there is no great difference about numbers of polyketides, alkanes and acids; (2) the number of olefins in spore is more than in fermentation liquor; (3) esters, aldehydes, benzenes and alcohols are unique to spore, and nitriles is unique to fermentation liquor. 1-(4-bromobutyl)-2-piperidinone, 2,2,6,6-tetramethyl-4-piperidinone, 6-pentyl-2H-pyran-2-one, 2,6,10-trimethyl-tetradecane, 2,6,10-trimethyl-pentadecane, and oxalic acid, butyl 6-ethyloct-3-ylester in T. asperellum GDFS1009 were found, which are related with pathogen inhibition and induced resistance (Fig 5C). Additionally, via bioinformatics analysis, 16 gene clusters of polyketide synthases were discovered in T. asperellum GDFS1009, almost all of which are lowly expressed or non-detectable after 24 h and 48 h in PD broth (Table 2).

Antimicrobial peptides exhibit broad antimicrobial abilities and strong non-specific resistance to pathogens. In 2005, Nielsen isolated an antimicrobial peptide polymer with 20 residues from T. brevicompactum [18]. In 2006, Kada isolated a new, short antimicrobial peptide from Trichoderma [66]. In 2007, Ruiz identified 11 residual fragments from antimicrobial peptides in T. longibrachiatum [19]. In this study, six types of antimicrobial peptides were detected and speculated from T. asperellum GDFS1009: TBV, trichotoxin_A-50_F, hypomurocin_B_IIIa, trichotoxin_A-50_I, hypomurocin_B_I, and hypomurocin_B_II (Table 3). As antimicrobial peptides are formed through the catalysis of NRPS, related gene clusters were analyzed and 16 NRPSs were found in T. asperellum GDFS1009, with two PKS/NRPS hetero-gene clusters. After 24 h and 48 h in PD broth, almost all of these gene clusters were lowly expressed or not detected (Table 2).

Based on these results, we infer that there are many mycoparasitism-related enzymes, elicitors, primary and secondary metabolites in Trichoderma that have not yet been discovered, including chitinase, protease, polyketides, antimicrobial peptides and so on. Future experiments should focus on optimizing the ingredients of the induction broth and the incubation conditions during the fermentation process. In addition, the silent genes and gene clusters should be expressed in Saccharomycetes. In conclusion, Trichoderma spp. offer great promise for ongoing agricultural development, expressing many promising enzymes, elicitors, and secondary metabolites.

Conclusion

Due to its efficient broad-spectrum antimicrobial activity, Trichoderma has been established as an internationally recognized biocontrol fungus. Though researches about Trichoderma spp. were very broad, most just concentrated in one or several biocontrol factors in one strain. Our study reveals resistant molecular mechanisms of T. asperellum in global aspects, including almost all biocontrol factors related to mycoparasitism, induced resistance and antibiosis, based on genomics, transcriptomics, metabolomics and a series of physiological and biochemical analyses. Firstly, we found and identified a novel strain of T. asperellum, named GDFS1009. The mycelium of T. asperellum GDFS1009 has high growth rate, high sporulation capacity, and very strong inhibitory effects on pathogens that cause cucumber fusarium wilt and corn stalk rot. T. asperellum GDFS1009 secretes chitinase, glucanase and protease, which can degrade the cell walls of fungi and participate in the mycoparasitism. The secreted xylanase are good candidates in inducing plant resistance and enhancing plant immunity. Then, RNA-seq and GC-MS analysis showed that T. asperellum GDFS1009 can produce nine primary metabolites which are precursors of antimicrobial compounds, and a variety of antimicrobial secondary metabolites, including polyketides, alkanes and so on. In addition, this study speculated six antimicrobial peptides via UPLC-QTOF-MS/MS. However, through synthetic analysis we can infer that there are some enzymes and elicitors with specific function and lots of secondary metabolites in Trichoderma have not been discovered. Subsequently, these related genes or gene clusters could be predicted and expressed via E. coli and Saccharomyces systems for wider application in the field of agricultural biocontrol.

Supporting information

S1 Fig. Control efficacy of T. asperellum GDFS1009 on F. gramimearum.

(A) Dual-culture assay, a1, Synergistic effect, a2, CK; (B) Resistant-dish assay, b1, Synergistic effect, b2, CK; (C) Micro well dilution assay.

https://doi.org/10.1371/journal.pone.0179957.s001

(DOCX)

S1 Table. Amino acid sequence homologies between mycoparasitism related enzymes in T. asperellum GDFS1009 and other strains.

https://doi.org/10.1371/journal.pone.0179957.s002

(DOCX)

S2 Table. Amino acid sequence homologies between induced resistant related enzymes in T. asperellum GDFS1009 and other strains.

https://doi.org/10.1371/journal.pone.0179957.s003

(DOCX)

S3 Table. Antibiosis secondary metabolites analysis in spore of T. asperellum GDFS1009.

https://doi.org/10.1371/journal.pone.0179957.s004

(DOCX)

S4 Table. Antibiosis secondary metabolites analysis in fermentation liquor of T. asperellum GDFS1009.

https://doi.org/10.1371/journal.pone.0179957.s005

(DOCX)

Acknowledgments

We sincerely thank Suzhou BioNovoGene Metabolomics Platform for the assistance in GC-MS analysis.

Author Contributions

  1. Conceptualization: JC QW RS.
  2. Data curation: MN JY LY CY KD JR.
  3. Formal analysis: QW RS MN.
  4. Funding acquisition: JC.
  5. Investigation: QW RS MN.
  6. Methodology: JC QW RS.
  7. Project administration: QW RS MN JC.
  8. Resources: MN JY YL CY KD JR.
  9. Software: QW RS MN.
  10. Supervision: QW RS MN JC.
  11. Validation: JC QW RS.
  12. Visualization: QW RS MN JC.
  13. Writing – original draft: QW JC.
  14. Writing – review & editing: QW JC MN.

References

  1. 1. Harman GE. Overview of mechanisms and uses of Trichoderma spp. Phytopathology. 2006;96: 190–194. pmid:18943924
  2. 2. Howell CR. Mechanisms employed by Trichoderma species in the biological control of plant diseases: the history and evolution of current concepts. Plant Dis. 2007;87: 4–10.
  3. 3. Harman GE, Howell CR, Viterbo A, Chet I, Lorito M. Trichoderma species- opportunistic, avirulent plant symbionts. Nat Rev Microbiol. 2004;2: 43–56. pmid:15035008
  4. 4. Sharon E, Chet I, Viterbo A, Bareyal M, Nagan H, Samuels GJ, et al. Parasitism of Trichoderma on Meloidogyne javanica and role of the gelatinous matrix. Eur J Plant Pathol. 2007;118: 247–258.
  5. 5. Bailey BA, Bae H, Strem MD, Crozier J, Thomas SE, Samuels GJ, et al. Antibiosis, mycoparasitism, and colonization success for endophytic Trichoderma isolates with biological control potential in Theobroma cacao. Biol Control. 2008;46: 24–35.
  6. 6. Weindling R. Studies on lethal principle effective in the parasitic action of Trichoderma hamatum on Rhisoctonia solani and other soil fungi. Phytopathology. 1932;22: 837–845.
  7. 7. Limon MC, Pintor-Toro JA, Benitez T. Increased antifungal activity of Trichoderma harzianum transformants that overexpress a 33-KDa chitinase. Phytopathology. 1999;89: 254–261. pmid:18944767
  8. 8. Viterbo A, Haran S, Friesem D, Ramot O, Chet I. Antifungal activity of a novel endochitinase gene (chit36) from Trichoderma harzianum Rifai TM. FEMS Microbiol Lett. 2001;200: 169–174. pmid:11425470
  9. 9. Geremia RA, Goldman GH, Jacobs D, Ardiles W, Vila SB, Van Montagu M, et al. Molecular characterization of the proteinase-encoding gene, prb1, related to mycoparasitism by Trichoderma harzianum. Mol Microbiol. 1993;8: 603–613. pmid:8326868
  10. 10. Vazquea-Garcduenase S, Leal CA, Herrera EA. Analysis of the β-1, 3-glucanolytic system of the biocontrol agent Trichoderma harzianum. Appl Environ Microbiol. 1998;64: 1442–1446. pmid:16349547
  11. 11. de la Cruz J, Pintor-Toro JA, Benítez T, Llobell A, Romero LC. A novel endo-beta-1, 3-glucanase, BGN13.1, involved in the mycoparasitism of Trichoderma harzianum. J Bacteriol. 1995;177: 6937–6945. pmid:7592488
  12. 12. Dean JFD, Anderson JD. Ethylene biosynthesis-inducing xylanase II. Purification and physical characterization of the enzyme produced by Trichoderma viride. Plant Physiol. 1991;95: 316–323. pmid:16667971
  13. 13. Djonovic S, Pozo MJ, Dangott LJ, Howell CR, Kenerley CM. Sm1, a proteinaceous elicitor secreted by the biocontrol fungus Trichoderma virens induces plant defense responses and systemic resistance. Mol Plant Microbe Interact. 2006;19: 838–853. pmid:16903350
  14. 14. Djonovic S, Vargas WA, Kolomiets MV, Horndeski M, Wiest A, Kenerley CM. A proteinaceous elicitor Sm1 from the beneficial fungus Trichoderma virens is required for induced systemic resistance in maize. Plant Physiol. 2007;145: 875–889. pmid:17885089
  15. 15. Nawrocka J, Małolepsza U. Diversity in plant systemic resistance induced by Trichoderma. Biol Control. 2013;67: 149–156.
  16. 16. Degenkolb T, Dohren H, Nielsen KF, Samuels GJ, Brückner H. Recent advances and future prospects in peptaibiotics, hydrophobin, and mycotoxin research, and their importance for taxonomy of Trichoderma and Hypocrea. Chem Biodivers. 2008;5: 671–680. pmid:18493954
  17. 17. Horace GC, Rjchard HC, Fanist GC, Patsy DC. 6-pentyl-α-pyrone from Trichoderma harzianum: Its plant growth inhibitory and antimicrobial propenies. Agric Biol Chem. 1986;50: 2943–2945.
  18. 18. Nielsen KF, Gräfenhan T, Zafari D, Thrane U. Trichothecene production by Trichoderma brevicompactum. J Agric Food Chem. 2005;53: 8190–8196. pmid:16218663
  19. 19. Ruiz N, Wielgosz-Collin G, Poirier L, Grovel O, Petit KE, Mohamed-Benkada M, et al. New Trichobrachins, 11-residue peptaibols from a marine strain of Trichoderma longibrachiatum. Peptides. 2007;28: 1351–1358. pmid:17629355
  20. 20. Stewart CN, Via LE. A rapid CTAB DNA isolation technique useful for RAPD finger printing and other PCR applications. Biotechniques. 1993;14: 748–750. pmid:8512694
  21. 21. Manter DK, Vivanco JM. Use of the ITS primers, ITS1F and ITS4, to characterize fungal abundance and diversity in mixed-template samples by qPCR and length heterogeneity analysis. J Microbiol Meth. 2007;71: 7–14.
  22. 22. Tamura K, Dudley J, Nei M, Kumar S. MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol. 2007;24: 1596–1599. pmid:17488738
  23. 23. Selvaraj K, Panneerselvam A. Invitro assessment of antagonistic activity of Trichoderma sp. against Sarocladium Oryzae causing sheath rot disease in paddy. Int J App Biol Pharm Technol. 2011;2: 179–183.
  24. 24. Kerroum F, Noureddine K, Eddine HJ, Mebrouk K. Biological control of Fusarium crown and root rot disease of tomato by Trichoderma harzianum in the west of algeria. Int J Sci Nat. 2015;6: 141–146.
  25. 25. Wu Q, Li Y, Li Y, Gao S, Wang M, Zhang T, et al. Identification of a novel fungus, Leptosphaerulina chartarum SJTU59 and characterization of its xylanolytic enzymes. PLoS ONE. 2013;8: e73729. pmid:24040044
  26. 26. Chen LL, Liu LJ, Shi M, Song XY, Zheng CY, Chen XL,et al. Characterization and gene cloning of a novel serine protease with nematicidal activity from Trichoderma pseudokoningii SMF2. FEMS Microbiol Lett. 2009;299: 135–142. pmid:19702879
  27. 27. Wolf GA, Wirth SJ. Soluble, dye-labelled substrates for a micro-plate assay of proteinase activity. J Microbiol Meth. 1996;25: 337–342.
  28. 28. Thordal-Christensen H, Zhang ZG, Wei YD, Collinge DB. Subcellular localization of H2O2 in plants. H2O2 accumulation in papillae and hypersensitive response during the barley-powdery mildew interaction. Plant J. 1997;11: 1187–1194.
  29. 29. Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009;25: 1105–1111. pmid:19289445
  30. 30. Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, et al. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010;28: 511–515. pmid:20436464
  31. 31. Kanehisa M, Goto S, Sato Y, Furumichi M, Tanabe M. KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Res. 2011;40: D109–D114. pmid:22080510
  32. 32. R Development Core Team. R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing. Vienna, Austria 2011. Available: http://www.R-project.org/.
  33. 33. Smith CA, Want EJ, Omaille G, Abagyan R, Siuzdak G. XCMS: processing mass spectrometry data for metabolite profiling using nonlinear peak alignment, matching, and identification. Anal Chem. 2006;78: 779–787. pmid:16448051
  34. 34. Vanholme R, Storme V, Vanholme B, Sundin L, Christensen JH, Goeminne G, et al. A systems biology view of responses to lignin biosynthesis perturbations in Arabidopsis. Plant Cell. 2012;24: 3506–3529. pmid:23012438
  35. 35. Kopka J, Schauer N, Krueger S, Birkemeyer C, Usadel B, Bergmüller E, et al. GMD@CSB.DB: the Golm metabolome database. Bioinformatics. 2005;21: 1635–1638. pmid:15613389
  36. 36. Kadu S, Hande D. Characterization of different secondary metabolites from Trimatostromma scutellare by GC-MS. World J Pharm Life Sci. 2016;2: 222–229.
  37. 37. Krause C, Kirschbaum J, Brückner H. Peptaibiomics: an advanced, rapid and selective analysis of peptaibiotics/peptaibols by SPE/LC-ES-MS. Amino Acids. 2006;30: 435–443. pmid:16622603
  38. 38. Kubicek CP, Bissett J, Druzhinina I, Kullnig-Gradinger C, Szakacs G. Genetic and metabolic diversity of Trichoderma: a case study on South-East Asian isolates. Fungal Genet Biol. 2003;38: 310–319. pmid:12684020
  39. 39. Hoyos-Carvajal L,Orduz S, Bissett J. Genetic and metabolic biodiversity of Trichoderma from Colombia and adjacent neotropic regions. Fungal Genet Biol. 2009;46: 615–631. pmid:19439189
  40. 40. Sadfi-Zouaoui N, Hannachi I, Rouaissi M, Hajlaoui MR, Rubio MB, Monte E, et al. Biodiversity of Trichoderma strains in Tunisia. Can J Microbiol. 2009;55: 154–162. pmid:19295648
  41. 41. Błaszczyk L, Popiel D, Chełkowski J, Koczyk G, Samuels GJ, Sobieralski K, et al. Species diversity of Trichoderma in Poland. J Appl Genetics. 2011;52: 233–243.
  42. 42. Metcalf DA, Dennis JC, Wilson CR. Effect of inoculum density of Sclerotium cepivorum on the ability of Trichoderma koningii to suppress white rot of onion. Plant Dis. 2004;88: 287–291.
  43. 43. Sun H, Yang J, Lin C, Huang X, Xing R, Zhang K. Purification and properties of a β-1,3-glucanase from Chaetomium sp. that is involved in mycoparasitism. Biotechnology Letters. 2006;28: 131–135. pmid:16369697
  44. 44. Liu Y, Yang Q. Cloning and heterologous expression of serine protease SL41 related to biocontrol in Trichoderma harzianum. J Molecul Microbiol Biotechnol. 2013;23: 431–439.
  45. 45. Yedidia I, Benhamou N, Chet I. Induction of defense responses in cucumber plants (Cucumiso sativus L.) by the biocontrol agent Trichoderma harzianum. Appl Environ Microbiol. 1999;65: 1061–1070. pmid:10049864
  46. 46. Yedidia I, Benhamou N, Kapulnik Y, Chet I. Induction and accumulation of PR proteins activity during early stages of root colonization by the mycoparasite Trichoderma harzianum T-203. Plant Physiol Biochem. 2000;38: 863–873.
  47. 47. Yedidia I, Shoresh M, Kerem Z, Benhamou N, Kapulnik Y, Chet I. Concomitant induction of systemic resistance to Pseudomonas syringue pv. Lachrymans in cucumber by Trichoderma aspereiium (T-203) and accumulation of phytoalexins. Appl Microbiol Microbiol. 2003;69: 7343–7353.
  48. 48. Chet I, Shoresh M, Yadidia I, Viterbo A. Mechanisms involved in biocontrol ang plant induced resistance by Trichoderma aspereiium (T. harzianum T-203). J Zhejiang Univ (Agric & Life Sci). 2004;30: 385.
  49. 49. Lin J. Synthesis of 2, 2-Dibromo-3-Nitrilopionamide and its application to water treatment system. Environ Sci Technol. 2006;29: 82–83.
  50. 50. Hollink Emily, Simanek EE Bergbreiter DE. Strategies for protecting and manipulating triazine derivatives. Tetrahedron Lett. 2005;46: 2005–2008.
  51. 51. Yue H, Zhao Y, Ma Xinbin, Gong J. Ethylene glycol: properties, synthesis, and applications. Chem Soc Rev. 2012;41: 4218–4244. pmid:22488259
  52. 52. Indira DG, Smitha P. Synthesis and antifungal studies of glycine and glycine-metal complexes on Phytophthora capsici. Inter Res J Bio Sci. 2013;2: 16–21.
  53. 53. Kulkarni A1, Patil SA, Badami PS. Synthesis, characterization, DNA cleavage and in vitro antimicrobial studies of La(III), Th(IV) and VO(IV) complexes with Schiff bases of coumarin derivatives. Eur J Med Chem. 2009;44: 2904–2912. pmid:19155104
  54. 54. Yang H, Ding M. Determination of o-toluic acid and its micro amounts of impurities in industrial products by HPLC. J Liq Chromatogr R T. 2002;25: 2709–2715.
  55. 55. Norman BH, Morris ML. A stereospecific synthesis of (−)-Bestatin from L-malic acid. Tetrahedron Lett. 1992;33: 6803–6806.
  56. 56. Graham AF, Lund BM. The effect of citric acid on growth of proteolytic strains of Clostridium botulinum. J Appl Bacteriol. 1986;61: 39–49. pmid:3531139
  57. 57. Yuan J, Raza W, Shen Q, Huang Q. Antifungal activity of Bacillus amyloliquefaciens NJN-6 volatile compounds against Fusarium oxysporum f. sp. cubense. Appl Environ Microbiol. 2012; 78: 5942–5944. pmid:22685147
  58. 58. Sumayo M, Hahm MS, Ghim SY. Determinants of plant growth-promoting Ochrobactrum lupini kudc1013 involved in induction of systemic resistance against Pectobacterium carotovorum subsp. carotovorum in tobacco leaves. Plant Pathol J. 2013;29: 174–181. pmid:25288944
  59. 59. Reena P, Nagar K. Antimicrobial activity of Trichoderma harzianum against bacteria and fungi. Int J Curr Microbiol Appl Sci. 2014;3: 96–103.
  60. 60. Vinale F, Sivasithamparam K, Ghisalberti EL, Marra R, Barbetti MJ, Li H, et al. A novel role for Trichoderma secondary metabolites in the interactions with plants. Physiol Mol Plant Pathol. 2008;72: 80–86.
  61. 61. Vinale F, Sivasithamparam K, Ghisalberti EL, Marra R, Woo SL, Lorito M. Trichoderma-plant-pathogen interactions in soil agro-ecosystems. Soil Biol Biochem. 2008;40: 1–10.
  62. 62. Viterbo A, Wiest A, Brotman Y, Chet I, Kenerley C. The 18mer peptaibols from Trichoderma virens elicit plant defence responses. Mol Plant Pathol. 2007;8: 737–746. pmid:20507534
  63. 63. Claydon N, Allan M, Hanson JR, Avent AG. Antifungal alkyl pyrones of Trichoderma harzianum.Trans Brit Mycol Soc. 1987;88: 503–513.
  64. 64. Zakrzewski J, Krawczyk M. Reactions of nitroxides. Part X: antifungal activity of selected sulfur and selenium derivatives of 2,2,6,6-tetramethylpiperidine. Bioorg Med Chem Lett. 2011, 21: 514–516. pmid:21074995
  65. 65. Jayaraj J, Bhuvaneswari R, Rabindran R, Muthukrishnan S, Velazhahan R. Oxalic acid-induced resistance to Rhizoctonia solani in rice is associated with induction of phenolics, peroxidase and pathogenesis-related proteins. J Plant Interact. 2010;5: 147–157.
  66. 66. Mohamed-Benkada M, Montagu M, Biard JF, Mondeguer F, Verite P, Dalgalarrondo M. New short peptaibols from a marine Trichoderma strain. Rapid Commun in Mass Spectrom. 2006;20: 1176–1180.