Skip to main content

Omics for understanding the tolerant mechanism of Trichoderma asperellum TJ01 to organophosphorus pesticide dichlorvos

Abstract

Backgroud

Though it is toxic to humans, dichlorvos is a widely used chemical pesticide and plays an important role in the control of plant pests. The application of a combination of the biocontrol agent Trichoderma with dichlorvos may reduce the need for chemical pesticides. Therefore, revealing the specific molecular mechanism of Trichoderma tolerance to dichlorvos has become particularly important.

Results

In this study, using transcriptome and metabolome analyses, changes in primary and secondary metabolisms in Trichoderma asperellum TJ01 were comprehensively studied in the presence of dichlorvos. A novel C2H2 zinc finger protein gene, zinc finger chimera 1 (zfc1), was discovered to be upregulated, along with a large number of oxidoreductase genes and ABC transporter genes under dichlorvos stress. In addition, gas chromatography-mass spectrometry (GC-TOF-MS), and liquid chromatography-mass spectrometry (LC-QQQ-MS) data revealed the global primary and secondary metabolic changes that occur in T. asperellum TJ01 under dichlorvos stress.

Conclusions

The tolerance mechanism of T. asperellum TJ01 to dichlorvos was proposed. In addition, the absorption and residue of dichlorvos were analyzed, laying the foundation for elucidation of the mechanism by which T. asperellum TJ01 degrades pesticide residues.

Background

Dichlorvos is an important organic phosphorus pesticide with fast-acting and broad-spectrum activity that has been in use in agricultural production for more than 50 years. However, it causes severe pollution of soil and water [1, 2]. The strong toxicity in non-target organisms poses a considerable threat to the environment and to humans [3]. Therefore, controlling pests with biological agents has garnered much attention in the field of plant protection.

Trichoderma is the most widely used and studied fungus in the biological control field. It offers the advantages of high activity and environmental safety. However, due to the limitations of its slow start-up and short shelf life, among others, Trichoderma cannot be used to completely replace chemical pesticides. At present, the ideal control strategy may be a combination of Trichoderma and dichlorvos, which would allow for a reduction in the amount of pesticides used while achieving the desired control effect [4]. Additionally, Trichoderma may contribute to the degradation of dichlorvos residues in the environment. As the tolerance of Trichoderma to dichlorvos directly affects the control afforded by the combination of the two pesticides, it is particularly important to investigate the mechanisms of Trichoderma tolerance to dichlorvos.

Until now, the tolerance mechanisms of Trichoderma to chemical pesticides have only been studied at the physiological and biochemical levels. These results have led to two hypotheses regarding the mechanisms of Trichoderma tolerance. The first is that Trichoderma is not recognized by chemical pesticides due to changes at target sites. For example, the main action of the fungicide benzimidazole is to bind the β-tubulin of a plant pathogen, thereby inhibiting its mitosis and morphogenesis. Trichoderma viride exhibits reduced affinity to pesticides due to mutations in the tubulin protein, resulting in its ability to tolerate benzimidazoles [5]. The second hypothesis is that Trichoderma effectively metabolizes or degrades some chemical pesticides. It was discovered that Trichoderma harzianum could degrade organochlorine insecticides via its oxidation system [6].

Transcription factors are important players in fungal tolerance to chemical pesticides and pharmaceutical drugs. In Saccharomyces cerevisiae, the zinc finger protein transcription factor pleiotropic drug resistance 1 (Pdr1) and its homologous protein Pdr3 are directly involved in the pleiotropic drug response and activate and suppress the expression of anti-drug-related genes, respectively. These regulators of drug-resistance-related genes include multiple ATP binding cassette (ABC) transporter genes, such as pdr5, pdr10, and pdr15 [7,8,9]. In Aspergillus fumigatus, the Zn2-Cys6 ABC-transporter regulating transcription factor (AtrR) provides resistance against azoles by regulating the expression of the drug target gene, Cytochrome P450 Family 51 Subfamily A (cyp51A) and the ABC transporter gene cdr1B [10]. In summary, transcription factors may regulate ABC transporter proteins for the removal of harmful substances from drug-resistant fungi and may also control the expression of the target protein recognized by the drug.

In this study, T. asperellum TJ01 was selected as it rapid growth, and high inhibition effect and tolerance to dichlorvos. The transcriptome and metabolome profiles were featured to study the tolerance mechanisms of T. asperellum TJ01 to dichlorvos. Meanwhile, the absorption and degradation of dichlorvos by T. asperellum TJ01 were also analyzed. The tolerance mechanisms discussed here should lay the foundation for use of T. asperellum TJ01 in the degradation of pesticide residues.

Results

Impact of dichlorvos on growth of T. asperellum TJ01

In this study, we determined the effect of dichlorvos on the growth of T. asperellum by solid medium cultrue. After culture on Potato Dextrose Aga (PDA) plates for 3 d, comparing with control, the fungi could grow well when exposured to dichlorvos (100 μg/mL), however, changes in the colony size and mycelium morphology of T. asperellum TJ01 were more and more evident as the concentration of dichlorvos increased (Fig. 1a). A nitrogen-poor medium, Burk medium, was also tested for this interaction. Similar tolerance phenotype was observed in a dichlorvos-dependent manner (Fig. 1b).

Fig. 1
figure 1

Morphology of T. asperellum TJ01 under dichlorvos stress after 3-d cultivation. a PDA plate; b Burk plate

Global impact of dichlorvos on T. asperellum TJ01

Based on RNA-seq analysis, in T. asperellum TJ01 treated with 100 μg/mL, 500 μg/mL, and 1000 μg/mL dichlorvos, 204, 490, and 872 genes were significantly upregulated, respectively, while 37, 177, and 383 genes were significantly downregulated, respectively (Fig. 2a).

Fig. 2
figure 2

Global impact of dichlorvos on T. asperellum TJ01 based on multiomics analysis. a Number of significantly differentially expressed genes based on RNA-seq analysis of T. asperellum TJ01 after treatment with dichlorvos for 24 h. D100: 100 μg/mL dichlorvos treatment/CK; D500: 500 μg/mL dichlorvos treatment/CK; D1000: 1000 μg/mL dichlorvos treatment/CK; b Compounds that were significantly altered in T. asperellum TJ01 after treatment with 500 μg/mL dichlorvos for 24 h. GC: GC-TOF-MS, LC(+): LC-QQQ-MS(+), LC(−): LC-QQQ-MS(−). Red: significantly upregulated, blue: significantly downregulated; c, The relationship of significantly upregulated genes in T. asperellum TJ01 treated with 100 μg/mL, 500 μg/mL, and 1000 μg/mL dichlorvos; d, The relationship of significantly downregulated genes in T. asperellum TJ01 treated with 100 μg/mL, 500 μg/mL, and 1000 μg/mL dichlorvos

Following GC-TOF-MS analysis of T. asperellum TJ01 for 24 h, 232 compounds were obtained. Compared with the control sample, the sample treated with dichlorvos exhibited 57 compounds that were significantly upregulated and 36 compounds that were significantly downregulated. Following LC-QQQ-MS (+) analysis of T. asperellum TJ01 for 24 h, 323 compounds were obtained. Compared with the control sample, the sample treated with dichlorvos treatment exhibited 69 compounds that were significantly upregulated and 71 compounds that were significantly downregulated. Finally, LC-QQQ-MS (−) analysis of T. asperellum TJ01 for 24 h resulted in the attainment of 173 compounds. Compared with the control sample, dichlorvos treatment led to significant upregulation of 48 compounds and significant downregulation of 22 compounds (Fig. 2b).

Then, the relationship of significantly changed genes in T. asperellum TJ01 treated with 100 μg/mL, 500 μg/mL, and 1000 μg/mL dichlorvos were illustrated by the venn diagram (Fig. 2c and d). While, the compounds detected via GC-TOF-MS, LC-QQQ-MS (+), and LC-QQQ-MS (−) have no coincident.

Impact of dichlorvos on primary metabolism of T. asperellum TJ01

Transcriptome profile analysis of T. asperellum TJ01 treated with 500 μg/mL dichlorvos showed that a citrate synthase, isocitrate dehydrogenase, and succinate dehydrogenase were significantly upregulated after treatment with dichlorvos (see Additional file 1). LC-QQQ-MS analysis of T. asperellum TJ01 treated with 500 μg/mL dichlorvos showed that citric acid and isocitrate were significantly downregulated, while succinic acid, fumaric acid, and malic acid were not significantly altered. Thus, comprehensive analysis of transcriptome and metabolome data suggested that dichlorvos has obvious effect on the tricarboxylic acid (TCA) cycle of T. asperellum TJ01 (Fig. 3a).

Fig. 3
figure 3

Impact of dichlorvos on primary metabolism of T. asperellum TJ01 based on multiomics analysis. a Impact of dichlorvos on TCA cycle based on RNA-seq, GC-TOF-MS, and LC-QQQ-MS analysis; b Impact of dichlorvos on metabolism of amino acids and derivatives based on GC-TOF-MS and LC-QQQ-MS analysis; c Impact of dichlorvos on lipid metabolism based on GC-TOF-MS and LC-QQQ-MS analysis; d Impact of dichlorvos on carbohydrate metabolism based on GC-TOF-MS and LC-QQQ-MS analysis. Red: significantly upregulated, blue: significantly downregulated, grey: not significantly altered, purple: rate-limiting step

KEGG analysis based on RNA-seq of T. asperellum TJ01 with 500 μg/mL dichlorvos treatment showed that 189 genes related to amino acid metabolism, 93 genes related to lipid metabolism, and 60 genes related to carbohydrate metabolism were significantly altered by dichlorvos. LC-QQQ-MS and GC-TOF-MS analysis of T. asperellum TJ01 treated with 500 μg/mL dichlorvos showed that, compared with the control, 16 (38%) amino acids and derivatives were significantly upregulated, while 26 (62%) were significantly downregulated (Fig. 3b; Additional file 2). Moreover, seven (39%) carbohydrates were significantly upregulated and 11 (61%) were significantly downregulated (Fig. 3c; Additional file 3), while 39 (89%) lipids were significantly upregulated and five (11%) were significantly downregulated (Fig. 3d; Additional file 4).

Impact of dichlorvos on secondary metabolism of T. asperellum TJ01

In T. asperellum TJ01 treated with 100 μg/mL, 500 μg/mL, and 1000 μg/mL dichlorvos, 32, 81, and 67 genes about secondary metabolism, respectively, were significantly altered based on RNA-seq analysis (Additional file 5). LC-QQQ-MS and GC-TOF-MS analysis of T. asperellum TJ01 treated with 500 μg/mL dichlorvos showed that, compared with the control, four flavonoids (5-methoxyflavanone, dihydromyricetin, afzelechin, and 4′-hydroxy-5,7-dimethoxyflavanone), one flavonol [isorhamnetin (3′-methoxyquercetin)], and three alkaloids (piperidine, pipecolate, and galactinol) were significantly upregulated, while one alkaloid (trigonelline) was significantly downregulated (Fig. 4; Additional file 6). In addition to these, many other secondary metabolites were significantly altered, such as griseofulvin, carnitine, and hycanthone.

Fig. 4
figure 4

Impact of dichlorvos on secondary metabolism of T. asperellum TJ01 based on GC-TOF-MS and LC-QQQ-MS analysis. a, Dihydromyricetin; b, 4'-hydroxy-5,7-dimethoxyflavanone; c, 5-methoxyflavanone; d, Trigonelline; e, Piperidine; f, Pipecolate. Red: significantly upregulated, blue: significantly downregulated

Impacts of dichlorvos on T. asperellum TJ01 transcription factors and stress resistance genes

In order to reveal the mechanisms of T. asperellum TJ01 tolerance to dichlorvos, detoxification- and tolerance-related genes were analyzed. Comparing with control, after 24 h of treatment with dichlorvos, five transcriptional regulators were significantly upregulated (Table 1). Meanwhile, a large number of redox proteins were significantly upregulated, including 49 peroxidases, 13 P450 monooxygenases, 23 glutathione metabolism proteins, two superoxide dismutases, five catalases, and 41 thioredoxins (Fig. 5; Additional file 7).

Table 1 Impact of dichlorvos with different concentrations on transcription regulator of T. asperellum TJ01 after 24 h treatment based on RNA-seq
Fig. 5
figure 5

Speculated mechanism of T. asperellum TJ01 resistance to dichlorvos. a Expression analysis of thioredoxin in T. asperellum TJ01 under the stress of different concentrations of dichlorvos. D100: log2 100 μg/mL dichlorvos treatment/CK, D500: log2 500 μg/mL dichlorvos treatment/CK, D1000: log2 1000 μg/mL dichlorvos treatment/CK; b Diagram of T. asperellum TJ01 resistance to dichlorvos. Red words represent significantly upregulated genes and their number

Additionally, after 24 h of treatment with dichlorvos, 13 ABC transporter genes were found to be significantly upregulated (Fig. 5b; Additional file 7).

Impact of dichlorvos on target genes of T. asperellum TJ01

Acetylcholinesterase (AChE) is a key enzyme that terminates nerve impulses by catalyzing the hydrolysis of the neurotransmitter acetylcholine in the nervous system. RNA-seq results showed that there were a total of 16 acetylcholinesterase genes in T. asperellum TJ01. Comparing with control, none were significantly altered following dichlorvos treatment, except for one under 1000 μg/mL dichlorvos (Table 2).

Table 2 Impact of dichlorvos with different concentrations on acetylcholinesterase of T. asperellum TJ01 based after 24 h treatment on RNA-seq

Absorption and residue degradation of dichlorvos by T. asperellum TJ01

UHPLC-MRM-MS/MS results showed that, after co-fermentation of T. asperellum TJ01 and 500 μg/mL dichlorvos in Burk broth for 3 d, the content of dichlorvos in T. asperellum TJ01 was 0.015 μg/mL (Fig. 6a) and the content of dichlorvos in the broth was 322.663 μg/mL (Fig. 6b). In contrast, after the pure fermentation of 500 μg/mL dichlorvos in Burk broth for 3 d, the content of dichlorvos was 415.908 μg/mL (Fig. 6c), while the self-degradation of dichlorvos was 84.092 μg/mL. This shows that T. asperellum TJ01 has a limited capacity to absorb dichlorvos, as most of the dichlorvos was still present in the fermentation broth.

Fig. 6
figure 6

Absorption and residue degradation of dichlorvos after co-fermentation with T. asperellum TJ01 for 3 d according to UHPLC-MS/MS detection. a Dichlorvos absorbed into T. asperellum TJ01; b Residual dichlorvos in broth after co-fermentation with T. asperellum TJ01; c Dichlorvos in broth following pure fermentation

Discussion

Trichoderma is one of the most widely used fungi in the field of biological control, as it exerts high activity and is environmentally friendly [11, 12]. However, due to the fact that it is slow-acting and has a short shelf life, in actual production, Trichoderma still cannot completely replace traditional chemical pesticides. At present, the most ideal control strategy may be a combination of Trichoderma and chemical pesticides, which would reduce the amount of pesticides used and achieve sustained control. Therefore, the tolerance of Trichoderma to dichlorvos is particularly critical. In previous work, it was found that T. atroviride had a good tolerance to dichlorvos [13], but there are no reports of the tolerance of T. asperellum. Recently, we isolated T. asperellum TJ01, which not only has a fast growth rate (strong competition), but also has high tolerance to dichlorvos. So this strain was chose as an in-depth study, hoping to lay a certain theoretical foundation for the production of complex microbial pesticides in agriculture.

There have been no thorough and comprehensive studies on the related mechanisms of Trichoderma tolerance to dichlorvos. Therefore, revealing the specific molecular mechanisms has become particularly important. In this study, transcriptome and metabolome data were used to analyze the mechanism underlying the response of Trichoderma asperellum TJ01 to dichlorvos. This represents the first detailed report of the global effects of dichlorvos on Trichoderma.

Primary metabolism is the process by which microorganisms absorb various nutrients from the outside world and produce energy and substances, including synthetic sugars, amino acids, common fatty acids, and nucleic acids, to sustain life activities through decomposition and anabolism. The TCA cycle is the basis of biological metabolism and is the hub of the carbohydrate, fatty acid, and amino acid metabolism pathways [14]. In this study, comprehensive analysis of transcriptome and metabolome data suggested that dichlorvos has obvious effect on primary metabolism of T. asperellum TJ01 (Fig. 3) [15]. Among the significantly upregulated lipid compounds, most were lysophosphatidylethanolamines and lysophosphatidylcholines. According to reports, most lysophosphatidylethanolamines and lysophosphatidylcholines are related to stress resistance, including LysoPC (lysophosphatidylcholine, 14:0), LysoPC (15:0), LysoPC (18:1), LysoPC (20:2), and LysoPC (20:1) (Additional file 4) [16, 17].

Secondary metabolism refers to the biosynthesis of non-essential substances and the storage of secondary metabolites. Trichoderma spp. secrete many secondary metabolites [18]. Some of these secondary metabolites are related to pathogen inhibition [19], while some are involved in stress resistance, such as flavonoids and alkaloids. Flavonoids act as a secondary antioxidant defense system in plant tissues exposed to different abiotic and biotic stresses [20]. It has also been reported that plant alkaloids may be impacted by drought stress and soil pH [21, 22]. However, there have been no similar reports in fungi and bacteria. In this study, four flavonoids and one alkaloid were significantly upregulated under the stress of dichlorvos (Fig. 4 and Additional file 6).

Transcription factors are important players in signal transduction pathways and determine the function of cells in many ways. There are 12 superfamily and 37 Protein family (PFAM) families containing transcription factor DNA binding domains in most fungal species [23, 24]. C2H2 zinc finger proteins are a large class of transcription factors in fungi, with 70–80 C2H2 zinc finger protein genes predicted in Trichoderma. It has been experimentally confirmed that carbon catabolite repression 1 (cre1) is involved in carbon metabolism [25], cseb1 is involved in response to osmotic stress, and Blue-light 7 (blu7) is involved in tolerance to light responsiveness [26]; however, a large number of genes remain to be characterized. In T. asperellum TJ01 treated with different concentrations of dichlorvos, a novel C2H2 zinc lipoprotein transcriptional regulator that may be involved in stress resistance (Table 1). Additionally, in order to reveal the mechanisms of T. asperellum TJ01 tolerance to dichlorvos, detoxification- and tolerance-related genes were analyzed [27, 28]. After dichlorvos treatment, a large number of redox proteins were significantly upregulated. This confirms that the ability of Trichoderma to withstand relatively high concentrations of a variety of synthetic and natural toxic compounds, including its own antibiotics, depends on efficient cell detoxification mechanisms, supported by a complex system of membrane pumps. It is well known that the genome of Trichoderma includes ABC transporters, which are members of a protein superfamily involved in the efflux of drugs from the cells of target organisms. These ABC transporters may explain the natural tolerance of Trichoderma to fungicides and their ability to successfully survive in extreme environments [28,29,30]. After dichlorvos treatment, lots of ABC transporter genes were found to be significantly upregulated. We speculate that all of the factors mentioned above contribute to the tolerance mechanism of T. asperellum TJ01 to dichlorvos (Fig. 5).

AChE is a key enzyme that terminates nerve impulses by catalyzing the hydrolysis of the neurotransmitter acetylcholine in the nervous system. Organophosphorus insecticides, such as dichlorvos, target AChE and irreversibly inhibit the enzyme by phosphorylating a serine hydroxyl group within the enzyme active site [31]. In this study, nearly none predicted AChE genes were significantly altered following dichlorvos treatment (Table 2). This may also be another reason for the fact that dichlorvos and T. asperellum TJ01 are able to co-exist at certain concentrations.

It has been reported that Trichoderma has a high tolerance and biosorption capacity of heavy metals like Zn2+, Pb2+, Ni3+, and Cu2+ [32]. However, the adsorption rate of dichlorvos by Trichoderma has not been reported. UHPLC-MRM-MS/MS results showed that T. asperellum TJ01 has a limited capacity to absorb dichlorvos. Interestingly, compared with the control, 93.23 μg/mL dichlorvos was missing from both the mycelium and the fermentation solution (Fig. 6). We speculated that this is being degraded by T. asperellum TJ01. It has already been reported that T. atroviride has the ability to degrade dichlorvos [33, 34]. Details of the degradation product of dichlorvos following degradation by T. asperellum TJ01 should be investigated in future studies.

This study provided a strong theoretical basis for the combination of dichlorvos and T. asperellum TJ01, which could be implemented in several ways, for example, dichlorvos with biocontrol factors or spores granules of T. asperellum TJ01. Therefore, future studies should evaluate the interaction of T. asperellum and dichlorvos in the presence of pathogens and plants.

Conclusions

In this study, association analyses of RNA-seq, GC-TOF-MS, and LC-QQQ-MS data revealed the global primary and secondary metabolic changes of T. asperellum TJ01 under dichlorvos stress. In addition, a novel C2H2 zinc finger protein gene, zfc1, was discovered and found to be upregulated in the presence of dichlorvos, along with a large number of oxidoreductase and ABC transporter genes, and then the tolerance mechanism of T. asperellum TJ01 to dichlorvos was proposed. Finally, the absorption and residue of dichlorvos were analyzed by LC-QQQ-MS, laying the foundation for the use of T. asperellum TJ01 in the degradation of pesticide residues.

Methods

Strain and reagent

The biocontrol fungus T. asperellum TJ01 was isolated from irrigated soil sample was collected in Fuyang, Anhui, China and stored at the School of Biotechnology and Food Engineering, Fuyang Normal College, Fuyang, Anhui, China. An analytical-grade standard of dichlorvos was purchased from the Shanghai Pesticide Research Institute Company.

Growth of T. asperellum TJ01 under dichlorvos stress

T. asperellum TJ01 was cultured in PDA dishes for 3 d, after which holes were punched in the T. asperellum TJ01 culture dish with a 7-mm punch. The dish was then turned upside-down over a fresh PDA or Burk (more suitable for studying mechanism as the clean background) plate [35] containing 0 μg/mL (control), 100 μg/mL, 500 μg/mL, or 1000 μg/mL dichlorvos. Plates were incubated at 28 °C for 3 d. The overall growth and sporulation of T. asperellum TJ01 were observed and photographed each day.

Effect of dichlorvos on transcription of T. asperellum TJ01 as determined by RNA-seq

The spore suspension of T. asperellum TJ01 was transferred to PD medium. After culturing at 28 °C and 180 rpm for 24 h, 2 g of mycelia was collected by vacuum filtration and transferred to Burk broth containing 0 μg/mL (control), 100 μg/mL, 500 μg/mL, or 1000 μg/mL dichlorvos. After culturing at 28 °C and 180 rpm for 24 h, the mycelia were collected by vacuum filtration and frozen with liquid nitrogen. RNA extraction, cDNA library construction, and Solexa sequencing were performed at the Beijing Genomics Institute (BGI, Beijing, China). Each experiment was repeated three times.

According to the reference genome (downloaded from http://genome.jgi.doe.gov/Trias1/Trias1.download.html) [36], conventional transcriptome analyses were performed, including GO (Gene Ontology) /KEGG (Kyoto Encyclopedia of Genes and Genomes) enrichment analysis of significantly differentially expressed genes [37, 38]. The NOISeq method was used to screen differently expressed miRNAs between two groups [39]. A log2(fold change) > 1 and probability > 0.8 indicated genes that were significantly upregulated, while a log2(fold change) < − 1 and probability > 0.8 indicated genes that were significantly downregulated.

Specific analyses of primary and secondary metabolic pathways were performed using the JGI website (https://genome.jgi.doe.gov/portal/).

Impact of dichlorvos on metabolome of T. asperellum TJ01 as determined by GC-TOF-MS analysis

The spore suspension of T. asperellum TJ01 was transferred to PD medium. After culturing at 28 °C and 180 rpm for 24 h, 2 g of mycelia was collected by vacuum filtration and transferred to Burk broth containing 0 μg/mL (control) or 500 μg/mL dichlorvos. After culturing at 28 °C and 180 rpm for 24 h, the mycelia were collected by vacuum filtration and frozen with liquid nitrogen.

GC-TOF-MS sample processing was performed as previously described [40]. GC-TOF-MS analysis was performed using an Agilent 7890 gas chromatography system coupled with a Pegasus HT time-of-flight mass spectrometer.

The system was equipped as previously described [40]. A 1-μL aliquot of the analyte was injected in splitless mode. Helium was used as the carrier gas, and the gas flow rate through the column was 1 mL min− 1. The initial temperature was kept at 50 °C for 1 min, then raised to 310 °C at a rate of 10 °C min− 1, then maintained at 310 °C for 8 min. The injection, transfer line, and ion source temperatures were 280, 270, and 220 °C, respectively. The energy was − 70 eV in electron impact mode. The mass spectrometry data were acquired in full-scan mode with an m/z range of 50–500 at a rate of 20 spectra per second after a solvent delay of 6.1 min [41, 42].

Chroma TOF 4.3X software (LECO Corporation) and the LECO-Fiehn Rtx5 database were used for raw peak exaction, baseline filtering and calibration, peak alignment, deconvolution analysis, peak identification, and integration of the peak area. Metabolic features detected in < 50% of QC samples were removed [42].

Global impact of dichlorvos on T. asperellum TJ01 metabolome as determined by LC-QQQ-MS analysis

The spore suspension of T. asperellum TJ01 was transferred to PD medium. After culturing at 28 °C and 180 rpm for 24 h, 2 g of mycelia was collected by vacuum filtration and transferred to Burk broth containing 0 μg/mL (control) or 500 μg/mL dichlorvos. After culturing at 28 °C and 180 rpm for 24 h, the mycelia were collected by vacuum filtration and frozen with liquid nitrogen.

The sample processing was performed as previously describeds [43]. Then the sample extracts were analyzed using an LC-ESI-MS/MS system (HPLC, Shim-pack UFLC SHIMADZU CBM30A system, www.shimadzu.com.cn/; MS, Applied Biosystems 4500 Q TRAP, https://sciex.com/products/mass-spectrometers/qtrap-systems/qtrap-4500-system). The analytical conditions for HPLC were as previously described [44], except that the column was that of Waters ACQUITY UPLC HSS T3 C18 (1.8 μm, 2.1 mm × 100 mm) and the gradient program was modified, 100:0 V/V at 0 min, 5:95 V/V at 11.0 min, 5:95 V/V at 12.0 min, 95:5 V/V at 12.1 min, 95:5 V/V at 15.0 min; flow rate, 0.40 mL/min. QQQ scans were acquired as MRM experiments and a specific set of MRM transitions were monitored for each period according to the metabolites eluted within this period [45].

Based on the self-built Metware database (MWDB) and the public database of metabolite information, qualitative analysis of the first-order and second-order spectra detected by mass spectrometry were performed. Some of substances were qualitatively analyzed by removing isotopic signals; repetitive signals containing K+, Na+, and NH4+ ions; and repeated signals of fragmented ions that themselves are of larger molecular weight. Metabolite structural elucidation was performed in reference to MassBank (https://massbank.eu/MassBank/), KNAPSAcK (http://kanaya.naist.jp/KNApSAcK/), HMDB (http://www.hmdb.ca/), and existing mass spectrometry public databases such as MoTo DB (http://www.ab.wur.nl/moto/) and METLIN (http://metlin.scripps.edu/index.php) [46].

Degradation and absorption of dichlorvos by T. asperellum TJ01 as determined by UHPLC-MS/MS

T. asperellum TJ01 spore suspension was transferred to PD medium to a final spore concentration of 106 cfu/mL. After culturing at 28 °C and 180 rpm for 2 d, 2 g of seed mycelium was transferred to Burk broth containing 500 μg/mL dichlorvos. After culturing at 28 °C and 180 rpm for 3 d, the culture was filtered through a 0.22-μm membrane filter. The mycelia of T. asperellum TJ01 were collected by vacuum filtration and washed twice with 4 °C sterile deionized water. Filtered Burk broth containing 500 μg/mL dichlorvos was used as a control.

For the liquid sample, a 10-μL aliquot of the sample was transferred to an EP tube with 990 μL of extraction solvent (acetonitrile:methanol:water, 2:2:1, precooled to − 20 °C). The mixture was vortexed for 30 s and centrifuged for 15 min at 12000 rpm and 4 °C. Then, the supernatant was transferred to an auto-sampler vial for UHPLC-MS/MS analysis (dilution factor = 100).

For the mycelim sample, 161.8 mg was precisely weighed and transferred to an Eppendorf tube with 1000 μL of extraction solvent (acetonitrile:methanol:water, 2:2:1, precooled to − 20 °C). The mixture was vortexed for 30 s, homogenized for 4 min at 45 Hz, and sonicated for 5 min in an ice-water bath. Homogenization and sonication were repeated three times, followed by incubation at − 20 °C for 1 h and centrifugation for 15 min at 12000 rpm and 4 °C. An 80-μL aliquot of the supernatant was transferred to an auto-sampler vial for LC-MS/MS analysis.

Abbreviations

ABC transporter:

ATP binding cassette transporter

AChE:

Acetylcholinesterase

ATP:

Adenosine triphosphate

AtrR:

ABC-transporter regulating transcription factor

BGI:

Beijing Genomics Institute

blu7:

Blue-light 7

cre1:

Carbon catabolite repression

cyp51A:

Cytochrome P450 Family 51 Subfamily A

GC-TOF-MS:

Gas chromatography-mass spectrometry

GO:

Gene Ontology

KEGG:

Kyoto Encyclopedia of Genes and Genomes

LC-QQQ-MS:

Liquid chromatography-mass spectrometry

LysoPC:

Lysophosphatidylcholine

PDA:

Potato Dextrose Aga

Pdr1:

Pleiotropic drug resistance 1

PFAM:

Protein family

TCA:

Tricarboxylic acid

UHPLC:

Ultra-high performance liquid chromatography

zfc1:

zinc finger chimera 1

References

  1. Lari SZ, Khan NA, Gandhi K, Meshram TS, Thacker NP. Comparison of pesticide residues in surface water and ground water of agriculture intensive areas. Iran J Environ Heal T. 2014;12(1):11.

    Article  CAS  Google Scholar 

  2. Anwar T, Ahmad I, Tahir S. Determination of pesticide residues in soil of Nawabshah district, Sindh, Pakistan. Pak J Zool. 2012;44(1):87–93.

    CAS  Google Scholar 

  3. Koutros S, Mahajan R, Zheng T, Hoppin JA, Ma X, Lynch CF, Blair A, Alavanja MCR. Dichlorvos exposure and human cancer risk: results from the agricultural health study. Cancer Cause Control. 2008;19(1):59–65.

    Article  Google Scholar 

  4. Wu Q, Zhang L, Xia H, Yu C, Dou K, Li Y, Chen J. Omics for understanding synergistic action of validamycin a and Trichoderma asperellum GDFS1009 against maize sheath blight pathogen. Sci Rep. 2017;7:40140.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  5. Goldman GH, Temmerman W, Jacobs D, Contreras R, Van Montagu M, Herrera-Estrella A. A nucleotide substitution in one of the β- tubulin genes of Trichoderma viride confers resistance to the antimitotic drug methyl benzimidazole-2-yl-carbamate. Mol Gen Genet. 1993;240(1):73–80.

    Article  PubMed  CAS  Google Scholar 

  6. Katayama A, Matsumura F. Degradation of organochlorine pesticides, particularly endosulfan, by Trichoderma harzianum. Environ Toxicol Chem. 1993;12(6):1059–65.

    Article  CAS  Google Scholar 

  7. Hagiwara D, Miura D, Shimizu K, Paul S, Ohba A, Gonoi T, Watanabe A, Kamei K, Shintani T, Moye-Rowley WS, et al. A novel Zn2-Cys6 transcription factor AtrR plays a key role in an azole resistance mechanism of Aspergillus fumigatus by co-regulating cyp51A and cdr1B expressions. PLoS Pathog. 2017;13(1):e1006096.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  8. Moye-Rowley WS. Transcriptional control of multidrug resistance in the yeast Saccharomyces. Prog Nucleic Acid Res Mol Biol. 2003;73:251–79.

    Article  PubMed  CAS  Google Scholar 

  9. Jungwirth H, Kuchler K. Yeast ABC transporters - a tale of sex, stress, drugs and aging. FEBS Lett. 2006;580(4):1131–8.

    Article  PubMed  CAS  Google Scholar 

  10. Delaveau T, Delahodde A, Carvajal E, Subik J, Jacq C. PDR3, a new yeast regulatory gene, is homologous to PDR1 and controls the multidrug resistance phenomenon. Mol Gen Genet. 1994;244(5):501–11.

    Article  PubMed  CAS  Google Scholar 

  11. Tewari L, Bhanu C. Evaluation of agro-industrial wastes for conidia based inoculum production of bio-control agent: Trichoderma harzianum. J Sci Ind Res India. 2004;63(10):807–12.

    CAS  Google Scholar 

  12. Brunner K, Zeilinger S, Ciliento R, Woo SL, Lorito M, Kubicek CP, Mach RL. Improvement of the fungal biocontrol agent Trichoderma atroviride to enhance both antagonism and induction of plant systemic disease resistance. Appl Environ Microb. 2005;71(7):3959–65.

    Article  CAS  Google Scholar 

  13. Tang J, Liu L, Huang X, Li Y, Chen Y, Chen J. Proteomic analysis of Trichoderma atroviride mycelia stressed by organophosphate pesticide dichlorvos. Can J Microbiol. 2010;56(2):121–7.

    Article  PubMed  CAS  Google Scholar 

  14. Mailloux RJ, Beriault R, Lemire J, Singh R, Chenier D, Hamel R, Appanna VD. The tricarboxylic acid cycle, an ancient metabolic network with a novel twist. PLoS One. 2007;2(8):e690.

  15. Arvas M, Pakula T, Lanthaler K, Saloheimo M, Valkonen M, Suortti T, Robson GD, Penttila M. Common features and interesting differences in transcriptional responses to secretion stress in the fungi Trichoderma reesei and Saccharomyces cerevisiae. BMC Genomics. 2006;7(1):32.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  16. Gao W, Li H, Xiao S, Chye M. Acyl-CoA-binding protein 2 binds lysophospholipase 2 and lysoPC to promote tolerance to cadmium-induced oxidative stress in transgenic Arabidopsis. Plant J. 2010;62(6):989–1003.

    PubMed  CAS  Google Scholar 

  17. Chen Z, Yoshida Y, Saito Y, Noguchi N, Niki E. Adaptive response induced by lipid peroxidation products in cell cultures. FEBS Lett. 2006;580(2):479–83.

    Article  PubMed  CAS  Google Scholar 

  18. Degenkolb T, Von Dohren H, Nielsen KF, Samuels GJ, Bruckner H. Recent advances and future prospects in peptaibiotics, hydrophobin, and mycotoxin research, and their importance for chemotaxonomy of Trichoderma and Hypocrea. Chem Biodivers. 2008;5(5):671–80.

    Article  PubMed  CAS  Google Scholar 

  19. Elhasan A, Walker F, Buchenauer H. Trichoderma harzianum and its metabolite 6-pentyl-α-pyrone suppress fusaric acid produced by Fusarium moniliforme. J Phytopathol. 2008;156(2):79–87.

    Article  CAS  Google Scholar 

  20. Agati G, Azzarello E, Pollastri S, Tattini M. Flavonoids as antioxidants in plants: location and functional significance. Plant Sci. 2012;196:67–76.

    Article  PubMed  CAS  Google Scholar 

  21. Jansen G, Jurgens H, Schliephake E, Ordon F. Effect of the soil pH on the alkaloid content of Lupinus angustifolius. Inter J Agron. 2012;2012:1–5.

    Article  CAS  Google Scholar 

  22. Phimchan P, Chanthai S, Bosland PW. Impact of drought stress on the accumulation of capsaicinoids in capsicum cultivars with different initial capsaicinoid levels. Hortscience. 2012;47(9):1204–9.

    Google Scholar 

  23. Shelest E. Transcription factors in fungi. FEMS Microbiol Lett. 2008;286(2):145–51.

    Article  PubMed  CAS  Google Scholar 

  24. Wingender E, Dietze P, Karas H, Knuppel R. TRANSFAC: a database on transcription factors and their DNA binding sites. Nucleic Acids Res. 1996;24(1):238–41.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  25. Portnoy T, Margeot A, Linke R, Atanasova L, Fekete E, Sandor E, Hartl L, Karaffa L, Druzhinina IS, Seiboth B, et al. The CRE1 carbon catabolite repressor of the fungus Trichoderma reesei: a master regulator of carbon assimilation. BMC Genomics. 2011;12:269.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  26. Cetz-Chel JE, Balcazar-Lopez E, Esquivel-Naranjo EU, Herrera-Estrella A. The Trichoderma atroviride putative transcription factor Blu7 controls light responsiveness and tolerance. BMC Genomics. 2016;17:327.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  27. Siddiquee S, Aishah SN, Azad SA, Shafawati SN, Naher L. Tolerance and biosorption capacity of Zn2+ , Pb2+ , Ni3+ and Cu2+ by filamentous fungi ( Trichoderma harzianum , T. aureoviride and T. virens ). Adv Biosci Biotechnol. 2013;04(04):570–83.

    Article  CAS  Google Scholar 

  28. Bae W, Chen X. Proteomic study for the cellular responses to Cd2+ in Schizosaccharomyces pombe through amino acid-coded mass tagging and liquid chromatography tandem mass spectrometry. Mol Cell Proteomics. 2004;3(6):596–607.

    Article  PubMed  CAS  Google Scholar 

  29. Ruocco M, Lanzuise S, Vinale F, Marra R, Turra D, Woo SL, Lorito M. Identification of a new biocontrol gene in Trichoderma atroviride: the role of an ABC transporter membrane pump in the interaction with different plant-pathogenic fungi. Mol Plant-Microbe Interact. 2009;22(3):291–301.

    Article  PubMed  CAS  Google Scholar 

  30. Waard DMA. Significance of ABC transporters in fungicide sensitivity and resistance. Pestic Sci. 1997;51(3):271–5.

    Article  Google Scholar 

  31. Assis CRD, Amaral IPG, Castro PF, Carvalho LB, Bezerra RDS. Effect of dichlorvos on the acetylcholinesterase from tambaqui (Colossoma macropomum) brain. Environ Toxicol Chem. 2007;26(7):1451–3.

    Article  PubMed  CAS  Google Scholar 

  32. Massaccesi G, Romero MC, Cazau MC, Bucsinszky AM. Cadmium removal capacities of filamentous soil fungi isolated from industrially polluted sediments, in La Plata (Argentina). World J Microb Biot. 2002;18(9):817–20.

    Article  CAS  Google Scholar 

  33. Sun W, Chen Y, Liu L, Tang J, Chen J, Liu P. Conidia immobilization of T-DNA inserted Trichoderma atroviride mutant AMT-28 with dichlorvos degradation ability and exploration of biodegradation mechanism. Bioresour Technol. 2010;101(23):9197–203.

    Article  PubMed  CAS  Google Scholar 

  34. Tang J, Liu L, Hu S, Chen Y, Chen J. Improved degradation of organophosphate dichlorvos by Trichoderma atroviride transformants generated by restriction enzyme-mediated integration (REMI). Bioresour Technol. 2009;100(1):480–3.

    Article  PubMed  CAS  Google Scholar 

  35. Wilson PW, Knight SC. Experiments in bacterial physiology. Mineapolis: Burguess; 1952.

    Google Scholar 

  36. Druzhinina IS, Chenthamara K, Zhang J, Atanasova L, Yang D, Miao Y, Rahimi MJ, Grujic M, Cai F, Pourmehdi S, et al. Massive lateral transfer of genes encoding plant cell wall-degrading enzymes to the mycoparasitic fungus Trichoderma from its plant-associated hosts. PLoS Genet. 2018;14(4):e1007322.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  37. Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009;25(9):1105–11.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  38. 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(5):511–5.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  39. Tarazona S, Garcia-Alcalde F, Dopazo J, Ferrer A, Conesa A. Differential expression in RNA-seq: a matter of depth. Genome Res. 2011;21(12):2213–23.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  40. Shang X, Zhong X, Tian X. Metabolomics of papillary thyroid carcinoma tissues: potential biomarkers for diagnosis and promising targets for therapy. Tumor Biol. 2016;37(8):11163–75.

    Article  CAS  Google Scholar 

  41. Kanehisa M, Sato Y, Kawashima M, Furumichi M, Tanabe M. KEGG as a reference resource for gene and protein annotation. Nucleic Acids Res. 2016;44:457–62.

    Article  CAS  Google Scholar 

  42. Wiklund S, Johansson E, Sjostrom L, Mellerowicz EJ, Edlund U, Shockcor JP, Gottfries J, Moritz T, Trygg J. Visualization of GC/TOF-MS-based metabolomics data for identification of biochemically interesting compounds using OPLS class models. Anal Chem. 2008;80(1):115–22.

    Article  PubMed  CAS  Google Scholar 

  43. Dong X, Gao Y, Chen W, Wang W, Gong L, Liu X, Luo J. Spatiotemporal distribution of phenolamides and the genetics of natural variation of hydroxycinnamoyl spermidine in rice. Mol Plant. 2015;8(1):111–21.

    Article  PubMed  CAS  Google Scholar 

  44. Chen W, Gong L, Guo Z, Wang W, Zhang H, Liu X, Yu S, Xiong L, Luo J. A novel integrated method for large-scale detection, identification, and quantification of widely targeted metabolites: application in the study of rice metabolomics. Mol Plant. 2013;6(6):1769–80.

    Article  PubMed  CAS  Google Scholar 

  45. Dunn WB, Broadhurst D, Begley P, Zelena E, Francismcintyre S, Anderson N, Brown M, Knowles JD, Halsall A, Haselden JN. Procedures for large-scale metabolic profiling of serum and plasma using gas chromatography and liquid chromatography coupled to mass spectrometry. Nat Protoc. 2011;6(7):1060–83.

    Article  PubMed  CAS  Google Scholar 

  46. Saccenti E, Hoefsloot HCJ, Smilde AK, Westerhuis JA, Hendriks MMWB. Reflections on univariate and multivariate analysis of metabolomics data. Metabolomics. 2014;10(3):361–74.

    Article  CAS  Google Scholar 

Download references

Funding

This work was supported by grants from the National Natural Science Foundation of China (31440025), Key Projects of Natural Science Research in Anhui Colleges and Universities (KJ2018A0339), Fuyang Normal University Research Project (2018HXXM21) and National Training Programs of Innovation and Entrepreneurship for Undergraduates (201610371025).

Availability of data and materials

The RNA-seq data can be found in the Gene Expression Omnibus (GEO) database at the National Center for Biotechnology Information (NCBI), under GEO accession number GSE114217.

Author information

Authors and Affiliations

Authors

Contributions

Conceived the study: JT. Generated the data: QW, MN, JT. Performed the analyses: QW, MN. Made substantial contributions to acquisition of data: GW, QL, MY. Wrote the manuscript: QW, MN, JT. Helped with writing and revision of the manuscript: GW, QL, MY. All authors read, performed critical revisions and approved the manuscript.

Corresponding author

Correspondence to Jun Tang.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Impact of dichlorvos with different concentrations on TCA cycle of T. asperellum TJ01 after 24 h treatment based on RNA-seq. Red represents up-regulated significantly, while blue represents down-regulated significantly. (XLS 7 kb)

Additional file 2:

Impact of 500 μg/mL dichlorvos on amino acid and derivatives metabolism of T. asperellum TJ01 after 24 h treatment based on LC-QQQ-MS and GC-TOF-MS analysis. Red represents up-regulated significantly, while blue represents down-regulated significantly. (XLS 11 kb)

Additional file 3:

Impact of 500 μg/mL dichlorvos on carbohydrate metabolism of T. asperellum TJ01 after 24 h treatment based on LC-QQQ-MS and GC-TOF-MS analysis. Red represents up-regulated significantly, while blue represents down-regulated significantly. (XLS 9 kb)

Additional file 4:

Impact of 500 μg/mL dichlorvos on lipid metabolism of T. asperellum TJ01 after 24 h treatment based on LC-QQQ-MS and GC-TOF-MS analysis. Red represents up-regulated significantly, while blue represents down-regulated significantly. (XLS 11 kb)

Additional file 5:

Number of significantly differentially expressed genes of secondary metabolism based on RNA-seq analysis of T. asperellum TJ01 after treatment with dichlorvos for 24 h. D100: 100 μg/mL dichlorvos treatment/CK; D500: 500 μg/mL dichlorvos treatment/CK; D1000: 1000 μg/mL dichlorvos treatment/CK. (TIF 504 kb)

Additional file 6:

Impact of 500 μg/mL dichlorvos on flavonoids and alkaloids of T. asperellum TJ01 after 24 h treatments based on LC-QQQ-MS and GC-TOF-MS analysis. Red represents up-regulated significantly, while blue represents down-regulated significantly. (XLS 9 kb)

Additional file 7:

Impact of dichlorvos with different concentrations on detoxification/tolerance mechanism related genes of T. asperellum TJ01 after 24 h treatment based on RNA-seq. Red represents up-regulated significantly, while blue represents down-regulated significantly. (XLS 24 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wu, Q., Ni, M., Wang, G. et al. Omics for understanding the tolerant mechanism of Trichoderma asperellum TJ01 to organophosphorus pesticide dichlorvos. BMC Genomics 19, 596 (2018). https://doi.org/10.1186/s12864-018-4960-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-018-4960-y

Keywords