Next Article in Journal
Deepening the Whole Transcriptomics of Bovine Liver Cells Exposed to AFB1: A Spotlight on Toll-like Receptor 2
Next Article in Special Issue
Preparation of Monoclonal Antibody against Deoxynivalenol and Development of Immunoassays
Previous Article in Journal
Could the Presence of Thrips AFFECT the Yield Potential of Genetically Modified and Conventional Maize?
Previous Article in Special Issue
Detection of αB-Conotoxin VxXXIVA (αB-CTX) by ic-ELISA Based on an Epitope-Specific Monoclonal Antibody
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Toxins
Volume 14
Issue 7
10.3390/toxins14070503
toxins-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Deoxynivalenol Biosynthesis in Fusarium pseudograminearum Significantly Repressed by a Megabirnavirus

by
Ke Li
1,†,
Dongmei Liu
2,†,
Xin Pan
1,
Shuwei Yan
1,
Jiaqing Song
1,
Dongwei Liu
1,
Zhifang Wang
1,
Yuan Xie
1,
Junli Dai
1,
Jihong Liu
2,
Honglian Li
1,
Xiaoting Zhang
1,* and
Fei Gao
1,*
1
Department of Plant Protection, Henan Agricultural University, Zhengzhou 450002, China
2
Institute of Agricultural Quality Standards and Testing Technology, Henan Academy of Agricultural Sciences, Zhengzhou 450002, China
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Toxins 2022, 14(7), 503; https://doi.org/10.3390/toxins14070503
Submission received: 16 June 2022 / Revised: 8 July 2022 / Accepted: 15 July 2022 / Published: 19 July 2022
(This article belongs to the Special Issue Research on Pathogenic Fungi and Mycotoxins in China (Volume II))
Browse Figures
Versions Notes

Abstract

:
Deoxynivalenol (DON) is a mycotoxin widely detected in cereal products contaminated by Fusarium. Fusarium pseudograminearum megabirnavirus 1 (FpgMBV1) is a double-stranded RNA virus infecting Fusarium pseudograminearum. In this study, it was revealed that the amount of DON in F. pseudograminearum was significantly suppressed by FpgMBV1 through a high-performance liquid chromatography–tandem mass spectrometry (HPLC-MS/MS) assay. A total of 2564 differentially expressed genes were identified by comparative transcriptomic analysis between the FpgMBV1-containing F. pseudograminearum strain FC136-2A and the virus-free strain FC136-2A-V-. Among them, 1585 genes were up-regulated and 979 genes were down-regulated. Particularly, the expression of 12 genes (FpTRI1, FpTRI3, FpTRI4, FpTRI5, FpTRI6, FpTRI8, FpTRI10, FpTRI11, FpTRI12, FpTRI14, FpTRI15, and FpTRI101) in the trichothecene biosynthetic (TRI) gene cluster was significantly down-regulated. Specific metabolic and transport processes and pathways including amino acid and lipid metabolism, ergosterol metabolic and biosynthetic processes, carbohydrate metabolism, and biosynthesis were regulated. These results suggest an unrevealing mechanism underlying the repression of DON and TRI gene expression by the mycovirus FpgMBV1, which would provide new methods in the detoxification of DON and reducing the yield loss in wheat.
Keywords:
deoxynivalenol; Fusarium pseudograminearum megabirnavirus 1; mycoviruses; transcriptome
Key Contribution: The secondary metabolite and predominant mycotoxin deoxynivalenol (DON) of Fusarium pseudograminearum was significantly reduced by a megabirnavirus FpgMBV1. The expression of genes in the trichothecene biosynthetic (TRI) gene cluster was suppressed. These results provide an idea to detoxify of DON by directly using mycovirus FpgMBV1 or its regulation on related genes.

1. Introduction

Fusarium is a genus of filamentous fungi ubiquitously existing in agricultural and natural ecosystems [1,2]. Some species in Fusarium cause diseases including wilts, blights, rots, and cankers on crops and some horticultural, ornamental, and forest plants [3,4]. Wheat crown rot (FCR), one of the most destructive wheat diseases worldwide, is caused by species in Fusarium [5]. The dominant causal agent of FCR is F. pseudograminearum [6,7,8]. The symptoms of FCR include dry seedlings at the seedling stage, browning and rot at the base of the stem at the adult stage, withered white ears at the filling stage, and shriveled grains at the harvest stage [9]. Besides the significant decrease in wheat yield, F. pseudograminearum produces a diverse array of toxic secondary metabolites including mycotoxins [10,11]. Cereals contaminated with mycotoxins are unsuitable for food or feed [12,13]. The mycotoxins produced by Fusarium include trichothecenes, fumonisins, etc. The most common Fusarium mycotoxin is deoxynivalenol (DON), which causes the disease of emesis, oral lesions, dermatitis, and hemorrhaging in human and livestock [14,15]. The biosynthesis of DON starts from the cyclization of farnesyl pyrophosphate mediated by trichodiene synthases (encoded by TRI5) [16]. Then, multiple steps are followed by at least six additional enzymes encoded by TRI4, TRI101, TRI11, TRI3, TRI1, and TRI8. Moreover, two transcription factors, TRI6 and TRI10; a transmembrane transporter, TRI12; and genes with unknown functions, TRI9 and TRI14, may also be involved in DON biosynthesis and the virulence of Fusarium [17,18,19,20,21].
Mycoviruses are viruses hosted in fungi. They have been described in many fungal species. Some of the viruses are hypovirulence-related with an effect on reducing the virulence of the pathogenic fungi. Besides the hypovirulence effect, some mycoviruses have an influence on the cultivating features of host fungi including the morphology and development of the colonies and spores. For example, Cryphonectria hypovirus 1 (CHV1) decreased the pigmentation and spore production of Cryphonectria parasitica [22]. Sclerotinia sclerotiorum hypovirulence-associated DNA virus 1(SsHDV1) caused abnormal colony morphology and the reduction in colony growth and small sclerotia production to Sclerotinia sclerotiorum [23]. Rosellinia necatrix megabirnavirus 1 (RnMBV1) suppressed the growth and melanin biosynthesis in Rosellinia necatrix [24].
To date, 29 mycoviruses species have been reported in Fusarium spp. They belong to the families Hypoviridae, Chrysoviridae, Totiviridae, Partitiviridae, Tymoviridae, Alternaviridae, and Megabirnaviridae [25]. Among them, Fusarium graminearum virus 1 (FgV1) [26], Fusarium graminearum mycovirus-China 9 (FgV-ch9) [27], Fusarium graminearum hypovirus 2 (FgHV2) [28], and Fusarium pseudograminearum megabirnavirus 1 (FpgMBV1) were reported as hypovirulence-related viruses [29]. Hypovirulence-related viruses hosted in F. graminearum, specifically FgV1, FgV-ch9, and FgHV2, caused a significant reduction in the fungal vegetative growth [26,27,28]. FpgMBV1 hosted in F. pseudograminearum [29]. There were mild changes in the colony morphology and spore production, but significant reduction in the virulence of F. pseudograminearum to wheat caused by FpgMBV1 [30]. Moreover, there were two viruses in F. graminearum, FgV1 and Fusarium graminearum mycotymovirus 1 (FgMTV1/SX64), which were reported to have a reducing effect on the DON production. FgV1 had a hypovirulence effect on F. graminearum, while FgMTV1/SX64 had a mild effect on the virulence [26,31].
To reveal the mechanism underlying the regulation of mycoviruses on fungi, high-throughput mRNA sequencing (RNAseq) was applied. For instance, the differential expression of genes related to fungal metabolism, transcription, translation, and ribosomal RNA processing was shown to be related to FgV1 infection [32]. Host cell transport-related genes were down-regulated by Fusarium graminearum virus 3 (FgV3) [33]. Genes involved in RNA processing and ribosome transport assembly were down regulated during the infection of Fusarium graminearum virus 4 (FgV4) [33]. Genes related to glutamate metabolism, homoserine metabolism, cellular aldehyde metabolism, and lactate metabolism in F. graminearum were down-regulated by FgHV1 [34]. RNA silencing and virulence-related genes were considerably down-regulated in S. sclerotiorum strains infected by SsHADV-1 [35].
In this study, the DON production of the FpgMBV1-containing F. pseudograminearum strain was determined and compared with that of the virus-free strain by HPLC-MS/MS. The genome-wide transcriptional reprogramming in F. pseudograminearum under the infection of FpgMBV1 was outlined by RNA sequencing. Genes involved in the biosynthesis of DON and ergosterol and in RNA-silencing pathways were analyzed to unravel the interaction between FpgMBV1 and F. pseudograminearum. Data obtained in this study would provide a clue to the specific regulation mechanisms of mycovirus to DON biosynthesis and virulence of fungal pathogens including Fusarium.

2. Results

2.1. DON Synthesis Was Inhibited by Mycovirus FpgMBV1

To investigate the effect of FpgMBV1 on the content of DON in F. pseudograminearum, the isogenic strains FC136-2A and FC136-2A-V- were cultivated in wheat grains for 30 days at 25 °C. Then, 20mL acetonitrile/water/acetic acid (70:29:1, v/v/v) was used to extract DON from 5 g of these mixtures for each sample for HPLC-MS/MS analysis. Results showed that the content of DON in the FpgMBV1-containing F. pseudograminearum strain FC136-2A was 1.3 ± 0.1 μg/kg. The content of DON in the virus-free strain FC136-2A-V- was 10.6 ± 1.4 μg/kg. Compared to the virus-free strain FC136-2A-V-, the content of DON in the FpgMBV1-containing strain FC136-2A was reduced by about 87% (p < 0.05) (Figure 1). No significant difference was found for the content of 15-ADON in the FpgMBV1-containing strain and the virus-free strain (Figure S2).

2.2. Overview of RNA-Seq Data for FC136-2A and FC136-2A-V-

To reveal how the biosynthesis genes of DON were regulated by FpgMBV1, the FpgMBV1-containing F. pseudograminearum strain FC136-2A and the isogenic virus-free strain FC136-2A-V- were used in the transcriptome analysis with three biological replicates for each strain using the DNBSEQ platform (BGI, Shenzhen, China). Gene expression profiles of the strains FC136-2A and FC136-2A-V- were compared. An average of 7.24 Gb of data was yielded for each sample. The genome of F. pseudograminearum CS3096 was used as the reference genome. An average of 94.57% and 94.27% of total reads (56–57 million) were aligned to the genome, respectively. A total of 11,500 genes were detected as expressed, of which 11417 were known and 83 were predicted as unknown genes (Table S1). The absolute fold change ≥ 2 and adjusted p-value ≤ 0.001 were used to define differently expressed genes (DEGs) (Figure 2A). Compared to the gene expression of the strain FC136-2A-V-, a total of 2564 statistically significant DEGs were found in the strain FC136-2A. Among them, 1585 genes were up-regulated and 979 genes were down-regulated (Figure 2B).
The PHI database annotation was screened using the DIAMOND software. A total of 131 genes were associated with loss of pathogenicity, 479 genes with reduced virulence, and 32 genes with lethal consequences (Figure 2C). In this study, the DEGs with PHI identity > 80% were used for further screening and the heat map construction. There were 18 genes categorized as reduced virulence genes including 4 up-regulated genes and 14 down-regulated ones (Figure S1). There were four genes categorized as essential genes (Figure S1). All these essential genes were up-regulated. No genes associated with loss of pathogenicity were identified in this study. Moreover, 12 of these genes were associated with DON and ergosterol biosynthesis (Table 1).
Gene ontology (GO) terms significantly enriched in the three major functional ontologies: 57.5% for biological process (BP), 10.1% for cellular component (CC) and 32.4% for molecular function (MF) (Figure 2D). For the 1585 up-regulated genes, the most enriched three for the BP were carbohydrate metabolic process, transmembrane transport, catechol-containing compound metabolic process, and oxidoreductase activity, catalytic activity, and cofactor binding for the MF (Figure 3A, Table S2). For the 979 down-regulated genes, the most enriched three of the BP were sterol metabolic and biosynthetic process, steroid metabolic and biosynthetic process, ergosterol metabolic and biosynthetic process; and oxidoreductase activity, iron ion binding, and cofactor binding for the MF (Figure 3B, Table S3). The most enrichedthree of the CC were integral component of membrane, an intrinsic component of membrane, membrane part, and membrane in all DEGs. There were 20 significantly enriched KEGG pathways for the up-regulated genes (Figure 3C, Table S4) and down-regulated genes, respectively (Figure 3D, Table S5). Sixteen pathways were related to metabolism.

2.3. The Metabolic Balance in F. pseudograminearum Was Disturbed by FpgMBV1

For the up-regulated DEGs, ten of the enriched KEGG pathways were related to metabolism, including tryptophan metabolism (ko00380), phenylalanine metabolism (ko00360), fatty acid metabolism (ko01212), biotin metabolism (ko00780), galactose metabolism (ko00052), arginine and proline metabolism (ko00330), starch and sucrose metabolism (ko00500), sphingolipid metabolism (ko00600), propanoate metabolism (ko00640), and linoleic acid metabolism (ko00591) (Figure 3C, Table S4). In addition, the meiosis–yeast pathway (ko04113) related to hyphal growth and development was also enriched.
Among the down-regulated genes, six pathways were related to metabolism (Figure 3D, Table S5). They were mainly classified into sulfur metabolism (ko00920), glycine, serine and threonine metabolism (ko00260), ether lipid metabolism (ko00565), fatty acid metabolism (ko01212), tyrosine metabolism (ko00350), and nitrogen metabolism (ko00910). Moreover, the biosynthesis of antibiotics, autophagy–yeast, steroid biosynthesis, and fatty acid degradation were also enriched.
The KEGG enrichment heat map shows that the fungal metabolic balance of amino acids and lipids was disturbed by FpgMBV1 (Figure 4). The transcription and translation of cell membrane-associated genes was also regulated by FpgMBV1.

2.4. TRI Genes Regulating DON Biosynthesis Were Down-Regulated by FpgMBV1

The genes in the TRI gene cluster encoding key enzymes for DON biosynthesis were FgTRI1-FgTRI16 and FgTRI101 in F. graminearum [36]. Comparative transcriptomic results showed that 12 genes in the TRI gene cluster were down-regulated by FpgMBV1. They were FpTRI1, FpTRI3, FpTRI4, FpTRI5, FpTRI6, FpTRI8, FpTRI10, FpTRI11, FpTRI12, FpTRI14, FpTRI15, and FpTRI101 (Figure 5). These TRI genes regulated by FpgMBV1 accounted for about two-thirds of the total gene number in the TRI gene cluster. Additionally, they were all significantly repressed.

2.5. Ergosterol Biosynthesis and Metabolism Were Inhibited by FpgMBV1

Fungal ergosterol content is significantly and positively correlated with DON content. In this study, ergosterol synthesis and metabolism were significantly enriched according to the results of GO enrichment of down-regulated genes (Figure 6A). In the PHI database annotation, four DEGs related to ergosterol synthesis were significantly down-regulated by FpgMBV1. They were FgERG4 (PHI:2728), FgERG5B (PHI:3038), FgERG5A (PHI:3037), and FgERG3B (PHI:3036) (Figure 6B). These results showed inhibition in ergosterol biosynthesis and metabolism in the FpgMBV1-containing strain FC136-2A.

2.6. FpDicer1 and FpAGO1 Were Down-Regulated by FpgMBV1

Based the RNA-seq result of FC136-2A and FC136-2A-V-, two Dicer (DCL) genes, two argonaute (AGO) genes, and one RNA-dependent RNA polymerase (RdR) gene were found differentially expressed. They were key genes in RNA silencing, which is an adaptive defense mechanism against foreign nucleic acids, especially viruses in animals, fungi, and plants [37]. Among them, FpDicer2(FPSE_06330), FpAGO2(FPSE_07737), and FpRdR1(FPSE_07737) were significantly up-regulated and FpDicer1(FPSE_07072) and FpAGO1(FPSE_00006) were down-regulated in the FpgMBV1-containing strain FC136-2A compared to the virus-free strain FC136-2A-V- (Figure 7). These results suggest that these RNA-silencing-related genes in F. pseudograminearum were involved in the defensive immunity against FpgMBV1.

2.7. Gene Expression Level by Quantitive Real-Time RT-PCR

Through quantitative real-time RT-PCR, the expression levels of 15 DEGs involved in metabolism, RNA silencing, virulence, and DON biosynthesis were confirmed. Primers are listed in Table S6. The expression levels of these representative genes were consistent with those in the transcriptomic data (Figure 8).

3. Discussion

In this study, the production of key mycotoxins DON and 15-ADON was found to be inhibited by FpgMBV1 significantly. The expression of 12 out of 15 TRI genes was repressed in accordance with the HPLC-MS/MS data. Previously, mycovirus FgV1 was reported to cause decreased DON production (60-fold) in F. graminearum [26]. However, only the TRI12 gene in the TRI gene cluster was differently expressed in the FgV1-containing strains [32]. For FgMTV1/SX64, another mycovirus reducing the DON concentration significantly in F. graminearum, whether TRI gene expression was regulated is unknown [31]. The global regulation of most TRI genes by a mycovirus was first reported for FpgMBV1. Two transcription regulators, TRI6 and TRI10, were key regulators in the TRI gene cluster [17]. TRI6 has been identified as a global transcription regulator, not only enhancing the expression of the genes in the DON biosynthetic pathway but also involved in the upstream isoprenoid pathway for trichothecene accumulation [38,39]. TRI10 regulated the expression of TRI1, TRI3, TRI8, TRI11, TRI12, TRI14, and TRI15 in Fusarium spp. [40]. In this study, the significant down-regulation of TRI6 and TRI10 might be the reason for the down-regulation of other TRI genes. The significant disruption of TRI gene expression resulted in the low contents of DON mycotoxin. This suggests an unrevealing mechanism underlying the repression of DON and TRI gene expression by mycovirus FpgMBV1, which would provide new methods in the detoxification of DON and reducing the yield loss in wheat.
Other genes of biosynthetic pathways were differentially expressed under the infection of FpgMBV1. For example, four genes encoding key enzymes in ergosterol synthesis were down-regulated. They were FgERG3B, FgERG4, FgERG5A, and FgERG5B [41,42]. It has been reported that ergosterol synthesis has a strong, positive correlation with the content of mycotoxin DON in the infected grains [43,44]. Moreover, ergosterol is an important constituent of fungal membranes [45]. Membranes are essential for many cellular processes, including the defensing response against viruses in fungi [46,47]. Another significantly down-regulated gene, HMG-CoA-reductase (HMR1), is the key enzyme in the mevalonate pathway. HMR1 is involved in the biosynthesis of many primary and secondary metabolites [48]. Some mevalonate derivatives function in fungi–plant interaction, such as isoprenoids and gibberellins [49]. Another down-regulated gene, FgVELB, plays a critical role in regulating various cellular processes and acts as a negative regulator for lipid biosynthesis. The deletion mutant of the FgVELB gene in F. graminearum produces a very low level of DON [50,51]. These results demonstrated the specific genes involved in the regulation of DON production Further studies on these genes would clarify the virulence-related pathways in fungi, especially Fusarium.
In general, the gene expression profile of F. pseudograminearum was reprogrammed by FpgMBV1 with 2564 genes differentially expressed. The disrupted metabolism caused by FpgMBV1 included amino acid and lipid metabolism, ergosterol metabolic and biosynthetic processes, carbohydrate metabolism and biosynthesis. Some of these pathways were critical in the transcriptome of other fungi under the infection of mycoviruses. In C. parasitica, “biosynthesis of other secondary metabolites”, “amino acid metabolism”, “carbohydrate metabolism”, and “translation” were enriched among the DEGs after CHV1 infection, demonstrating that virus infection resulted in massive but specific changes in primary and secondary metabolism. Some of the highly induced metabolites played key roles in the growth, development and pathogenicity of fungi [52]. The complex interaction between leucine metabolism and the global regulator of mycotoxin biosynthesis, TRI6, and virulence in F. graminearum has been explored [53]. In this study, the KEGG enrichment results showed that tryptophan, phenylalanine, tyrosine, glycine, serine, and threonine-related genes were differentially expressed under the infection of FpgMBV1. Further studies on the regulation mechanism of FpgMBV1 would help in revealing the crosstalk between some primary metabolic pathways and mycotoxin biosynthesis and virulence in F. pseudograminearum.
Mycoviruses were triggers and targets of RNA silencing. Three genes (FpDCL2, FpAGO2, and FpRdR1) were up-regulated and two genes (FpDCL1 and FpAGO1) were down-regulated in F. pseudograminearum under the infection of FpgMBV1. They are genes for key components in the RNA-silencing machinery [54]. Gene expression of DCL2 and AGO2 were induced by CHV1 in C. parasitica [55]. Genes of SsAgl2, SsDcl1, and SsDCl2 were essential in defending against viruses in S. sclerotiorum [56]. The expression of FgDICER2 and FgAGO2 were suppressed by FgV1 in F. graminearum [57]. Most RNA-silencing genes in S. sclerotiorum were repressed by SsHADV-1 [35]. Considering the high diversity of these virus–fungi systems, the different regulation on the RNA-silencing pathway seems reasonable. The regulation network of FpgMBV1 is partially revealed in this study, while the underlying mechanism is to be explored.

4. Conclusions

Megabirnavirus FpgMBV1 significantly repressed the production of mycotoxin DON and the expression of 12 TRI genes in F. pseudograminearum. The expression of ergosterol biosynthesis and RNA-silencing-related genes and genes involved in metabolism were regulated by FpgMBV1(Figure 9). Prospectively, FpgMBV1 is valuable in the detoxification of DON and the management of diseases caused by F. pseudograminearum.

5. Materials and Methods

5.1. Fungal Material and Growth Conditions

FpgMBV1-containing F. pseudograminearum strain FC136-2A and its isogenic FpgMBV1-free strain FC136-2A-V- were maintained in the laboratory, department of plant protection, Henan Agricultural University, Henan province, China [27]. For RNA sequencing, these two strains were grown at 25 °C in the dark on potato dextrose agar (PDA) medium (Becton, Dickinson, and Company, Sparks, MD, USA). For HPLC-MS/MS analysis, these strains were cultured in wheat grain media. Typically, 180 g of wheat seeds was soaked in distilled water for 12 h and boiled for 30 min then air-dried and autoclaved at 120 °C for 30 min in triangle flasks to complete the preparation of the wheat grain medium. Three-day-old PDA plugs of FC136-2A or FC136-2A-V- mycelia were cultured in wheat grain media for 30 days in the dark at 25 °C.

5.2. HPLC-MS/MS Analysis of Type B Trichothecene

Finely ground wheat grain (5.00 g) was weighted and extracted with 20 mL acetonitrile/water/acetic acid (70:29:1, v/v/v). The suspension was blended for 5 min at 2500 rpm using multi-tube vortexer UMV-2. The homogenized solution was centrifuged for 10 min after cooling (T = 4 ℃, 8000 rpm). Then, 750 µL of the upper layer and 2 mL water were mixed in a 15 mL tube at 11,000 rpm for 10 min. The upper layer was cleaned via nylon filters and used for HPLC-MS/MS analysis.
HPLC analysis was performed using LC-MS-8060NX (Shimadzu, Japan) in gradient conditions. The separation of the toxin was performed using 100 nm × 2.1 mm, 1.9 µm (Thermo Fisher Scientific, Shanghai, China). The column temperature was set at 40 ℃ and the injection volume was 1 µL. Solvent A was 1 mmoL·L-1 ammonium formate in 0.1% formic acid and solvent B was methanol and 1 mmoL·L-1 ammonium formate in 0.1% formic acid. A gradient at a flow rate of 0.2 mL·min-1 was performed within 26 min (Table S7).
MS/MS was performed on an 8060NX triple quadrupole mass spectrometer equipped with a Turbo Ion-Spray electrospray ionization (ESI) source (Shimadzu, Japan) heated at 400 °C in the negative (DON, 3-ADON, and 15-ADON) and positive (ZEN, D3G, and NIV) ionization mode. Quantitation was performed using multiple-reaction monitoring (MRM) with a dwell time of 100 ms. The following transition reactions of DON, 3ADON,15ADON, ZEN, D3G, NIV, and the IS with the respective declustering potential (DP), collision energy (CE), and cell exit potential (CXP) in parentheses were recorded using the first mass transition for quantitation. DON: m/z 297.1 (DP −21.0 V, CE −11.0 V, CXP −25.0V), 3ADON: m/z 339.2 (DP −30V, CE −14V, CXP −25V), 15ADON: m/z 356.2 (DP −18V, CE −9.0V, CXP −23V), D3G: m/z 503.0 (DP +26V, CE +20V, CXP +29V), ZEN: m/z 317.2 (DP +16V, CE +24V, CXP +29V), and NIV: m/z 357.2 (DP +13V, CE +13V, CXP +28V).
Data are reported as mean values ± SD of three biological replications. Type B trichothecene (NIV, DON, FX, 3-ADON, and 15-ADON) yields are expressed as μg·kg -1 of dry fungal biomass. Values were compared at the 1% significance level using DPS software (control vs. treated).

5.3. Total RNA Extraction

The 0.5 mg mycelia of strains FC136-2A and FC136-2A-V- were harvested from PDA plates and ground in liquid nitrogen. Total RNA isolation was conducted using TRIzol (Invitrogen, Carlsbad, CA, USA). The extracted RNA was treated with RNase-free DNase I (Promega, Madison, WI, United Kingdom). The concentration of the extracted RNAs was determined with Nanodrop system (NanoDrop, Madison, WI, USA). The RNA integrity was examined by the RNA integrity number (RIN) using an Agilent 2100 bioanalyzer (Agilent, Santa Clara, CA, USA).

5.4. cDNA Library Preparation and Sequencing

MGI Easy RNA Library Preparation Kit (BGI-Shenzhen, China) and DNA Clean Beads MGI Easy DNA Cleanbeads Kit (BGI-Shenzhen, China) were used to purify mRNA and synthesize cDNA and PCR amplification; then, purification of the PCR product was carried out following the kit instructions. The product was validated on the Agilent Technologies 2100 bioanalyzer for quality control. The double-stranded PCR products from the previous step were heated, denatured, and circularized by splint oligo sequencing to obtain the final library. The single-stranded circular DNA was formatted as the last library and amplified with phi29 polymerase to make a DNA nanoball, with more than 300 copies per molecule. DNBs were loaded into the patterned nanoarray, and 100 paired base reads were generated on the BGI seq500 platform (BGI-Shenzhen, China).

5.5. RNA-Seq Data Analysis

The raw data were filtered with BGI’s filtering software SOAPnuke (version: v1.4.0) to remove reads containing splices (splice contamination), reads with unknown base N content greater than 10%, and reads in which bases with a quality value of lower than 15 accounted for more than 20% of the total number of bases. Using HISAT2 software (version: v2.1.0; http://www.ccb.jhu.edu/software/hisat), the clean reads were mapped to the F. pseudograminearum CS3096 genome (Taxonomy ID: 1028729). Bowtie2 (version: v2.2.5; http://bowtie-bio.sourceforge.net/Bowtie2/index.shtml) was used to compare the reference gene sequences and RSEM (version: v1.2.8) was applied to calculate the expression levels of genes and transcripts. Pearson correlation coefficients between every two samples were obtained using the core function in the R software, and principal component analysis (PCA) was performed using the princomp method. To identify significant virulence genes in fungus, DIAMOND software was used to annotate genes with query coverage of 50% and identity of 40% in the PHI database (version: V0.8.31). To improve the accuracy of DEGs, we defined genes with more than two-fold difference and Q-values ≤ 0.001 to be screened as significantly differentially expressed genes. Differential genes were functionally classified according to GO and KEGG annotation results and official classification. At the same time, enrichment analysis was performed using the phyper function in R software, with FDR correction for p-value, and functions with Q value ≤ 0.05 were considered significantly enriched.

5.6. Quantitative Real-Time RT-PCR (qRT-PCR) Analysis

Total RNA samples from FC136-2A and FC136-2A-V- were used for cDNA synthesis by Goldenstar® RT6 cDNA Synthesis kit ver.2 (Tsingke Biotechnology Co., Ltd., Beijing, China). The resulting cDNAs were used as templates for qRT-PCR.
The qRT-PCR was carried out in a QuantStudio 3 Real-Time PCR System (Thermo Fisher Scientific, Waltham, MA, USA) with 2×T5 Fast qPCR Mix (Tsingke Biotechnology Co., Ltd., Beijing, China). PCR amplification was performed under the following conditions: 95 °C for 1 min, followed by 40 cycles of 95 °C for 10 s, and 60 °C for 30 s. Meltcurve plots were analyzed for each gene tested after each PCR reaction. The ubiquitin gene of TEF1α (FPSE_11980) served as an internal reference gene. Primers for 15 of the DEGs were designed using Primer Premier 5 and are listed in Table S6.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/toxins14070503/s1, Figure S1: The virulence genes according to the PHI-base. (A) Heat map of reduced-virulence genes in DEGs. (B) Heat map of essential genes (according to PHI database) in DEGs, Figure S2: Accumulation of mycotoxin 15-ADON in wheat cultures of FC136-2A and FC136-2A-V-. Table S1: Statistics of the transcriptome-sequencing results, Table S2: The GO enrichment analysis of the up-regulated genes, Table S3: The GO enrichment analysis of the down-regulated genes, Table S4: The KEGG enrichment analysis of the up-regulated genes, Table S5: The KEGG enrichment analysis of the down-regulated genes, Table S6: qRT-PCR primers used in this study, Table S7: Gradient elution conditions of mobile phase.

Author Contributions

Conceptualization, F.G. and X.Z.; data curation, X.P., Y.X. and X.Z.; formal analysis, K.L. and J.S.; funding acquisition H.L., J.L. and X.Z.; investigation X.P., S.Y. and D.L. (Dongmei Liu); methodology, F.G., D.L. (Dongmei Liu) and J.D.; resources, Y.X.; supervision, J.L.; validation, D.L. (Dongwei Liu) and Z.W.; writing—original draft, K.L.; writing—review and editing, X.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This work was financially supported by the Major International (Regional) Joint Research NSFC-CGIAR Project (31961143018), Henan Province Agricultural Science and Technology Research Project (212102110138), Henan Agricultural University Science and Technology Innovation Project (KJCX2020A14), New discipline development project of Henan Academy of Agricultural Sciences (2020XK06, 2022XK08).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available in Supplementary Materials Tables S1–S7 and Figures S1 and S2.

Acknowledgments

We thank the editor for editing assistance, and the anonymous reviewers for constructive and helpful comments. We also acknowledge the BGI database development team for technical support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Burgess, L.; Bryden, W. Fusarium: A ubiquitous fungus of global significance. Microbiology 2012, 33, 22. [Google Scholar] [CrossRef]
  2. Ma, L.J.; Geiser, D.M.; Proctor, R.H.; Rooney, A.P.; O’Donnell, K.; Trail, F.; Gardiner, D.M.; Manners, J.M.; Kazan, K. Fusarium pathogenomics. Annu. Rev. Microbiol. 2013, 67, 399–416. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Dean, R.; Van Kan, J.A.L.; Pretorius, Z.A.; Hammond-Kosack, K.E.; Di Pietro, A.; Spanu, P.D.; Rudd, J.J.; Dickman, M.; Kahmann, R.; Ellis, J.; et al. The Top 10 fungal pathogens in molecular plant pathology. Mol. Plant Pathol. 2012, 13, 314–430. [Google Scholar] [CrossRef] [Green Version]
  4. Rampersad, S. Pathogenomics and management of Fusarium diseases in plants. Pathogens 2020, 9, 340. [Google Scholar] [CrossRef] [PubMed]
  5. Petronaitis, T.; Simpfendorfer, S.; Hüberli, D. Importance of Fusarium spp. in wheat to food security: A global perspective. In Plant Pathology in the 21st Century; Springer: Berlin/Hamburg, Germany, 2021; Volume 10. [Google Scholar] [CrossRef]
  6. Li, H.L.; Yuan, H.X.; Fu, B.; Xing, X.P.; Sun, B.J.; Tang, W.H. First Report of Fusarium pseudograminearum Causing Crown Rot of Wheat in Henan, China. Plant Dis. 2012, 103, 18–19. [Google Scholar] [CrossRef]
  7. Alahmad, S.; Simpfendorfer, S.; Bentley, A.R.; Hickey, L.T. Crown rot of wheat in Australia: Fusarium pseudograminearum taxonomy, population biology and disease management. Australas. Plant Pathol. 2018, 47, 285–299. [Google Scholar] [CrossRef]
  8. Kazan, K.; Gardiner, D.M. Fusarium crown rot caused by Fusarium pseudograminearum in cereal crops: Recent progress and future prospects. Mol. Plant Pathol. 2018, 19, 1547–1562. [Google Scholar] [CrossRef] [Green Version]
  9. Zhou, H.; He, X.; Wang, S.; Ma, Q.; Sun, B.; Ding, S.; Chen, L.; Zhang, M.; Li, H. Diversity of the Fusarium pathogens associated with crown rot in the Huanghuai wheat-growing region of China. Environ. Microbiol. 2019, 21, 2740–2754. [Google Scholar] [CrossRef]
  10. Bragard, C.; Baptista, P.; Chatzivassiliou, E.; DiSerio, F.; Gonthier, P.; JaquesMiret, J.A.; Justesen, A.F.; MacLeod, A.; Magnusson, C.S.; Milonas, P.; et al. Pest categorisation of Fusarium pseudograminearum. EFSA J. 2022, 29, e07399. [Google Scholar] [CrossRef]
  11. Perincherry, L.; Justyna, L.; Lukasz, S. Fusarium-Produced Mycotoxins in Plant-Pathogen Interactions. Toxins 2019, 11, 664. [Google Scholar] [CrossRef] [Green Version]
  12. Borbély, M.; Sipos, P.; Pelles, F.; Győri, Z. Mycotoxin contamination in cereals. J. Agroaliment. Process. Technol. 2010, 16, 96–98. Available online: https://www.researchgate.net/publication/265178683 (accessed on 10 June 2010).
  13. Lozowicka, B.; Iwaniuk, P.; Konecki, R.; Kaczynski, P.; Kuldybayev, N.; Dutbayev, Y. Impact of Diversified Chemical and Biostimulator Protection on Yield, Health Status, Mycotoxin Level, and Economic Profitability in Spring Wheat (Triticum aestivum L.) Cultivation. Agronomy 2022, 12, 258. [Google Scholar] [CrossRef]
  14. Skrzydlewski, P.; Twarużek, M.; Grajewski, J. Cytotoxicity of Mycotoxins and Their Combinations on Different Cell Lines: A Review. Toxins 2022, 14, 244. [Google Scholar] [CrossRef]
  15. Pierzgalski, A.; Bryła, M.; Kanabus, J.; Modrzewska, M.; Podolska, G. Updated Review of the Toxicity of Selected Fusarium Toxins and Their Modified Forms. Toxins 2021, 13, 768. [Google Scholar] [CrossRef]
  16. Flynn, C.M.; Broz, K.; Jonkers, W.; Schmidt-Dannert, C.; Kistler, H.C. Expression of the Fusarium graminearum terpenome and involvement of the endoplasmic reticulum-derived toxisome. Fungal Genet. Biol. 2019, 124, 78–87. [Google Scholar] [CrossRef] [PubMed]
  17. Seong, K.Y.; Pasquali, M.; Zhou, X.; Song, J.; Hilburn, K.; McCormick, S.; Dong, Y.; Xu, J.R.; Kistler, H.C. Global gene regulation by Fusarium transcription factors Tri6 and Tri10 reveals adaptations for toxin biosynthesis. Mol. Microbiol. 2009, 72, 354–367. [Google Scholar] [CrossRef]
  18. Menke, J.; Dong, Y.N.; Kistler, H.C. Fusarium graminearum Tri12p influences virulence to wheat and trichothecene accumulation. Mol. Plant-Microbe Interact. 2012, 25, 1408–1418. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Alexander, N.J.; McCormick, S.P.; Hohn, T.M. TRI12, a trichothecene efflux pump from Fusarium sporotrichioides: Gene isolation and expression in yeast. Mol. Gen. Genet. 1999, 261, 977–984. [Google Scholar] [CrossRef]
  20. Kimura, M.; Tokai, T.; O’Donnell, K.; Ward, T.J.; Fujimura, M.; Hamamoto, H.; Shibata, T.; Yamaguchi, I. The trichothecene biosynthesis gene cluster of Fusarium graminearum F15 contains a limited number of essential pathway genes and expressed non-essential genes. FEBS Lett. 2003, 539, 105–110. [Google Scholar] [CrossRef] [Green Version]
  21. Lee, T.; Lee, S.H.; Shin, J.Y.; Kim, H.K.; Yun, S.H.; Kim, H.Y.; Lee, S.; Ryu, J.G. Comparison of trichothecene biosynthetic gene expression between Fusarium graminearum and Fusarium asiaticum. Plant Pathol. J. 2014, 30, 33–42. [Google Scholar] [CrossRef] [Green Version]
  22. Lin, H.; Lan, X.; Liao, H.; Parsley, T.B.; Nuss, D.L.; Chen, B. Genome Sequence, Full-Length Infectious cDNA Clone, and Mapping of Viral Double-Stranded RNA Accumulation Determinant of Hypovirus CHV1-EP721. J. Virol. 2007, 81, 1813–1820. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Yu, X.; Li, B.; Fu, Y.; Jiang, D.; Ghabrial, S.A.; Li, G.; Peng, Y.; Xie, J.; Cheng, J.; Huang, J.; et al. A geminivirus-related DNA mycovirus that confers hypovirulence to a plant pathogenic fungus. Proc. Natl. Acad. Sci. USA 2010, 107, 8387–8392. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Chiba, S.; Salaipeth, L.; Lin, Y.-H.; Sasaki, A.; Kanematsu, S.; Suzuki, N. A Novel Bipartite Double-Stranded RNA Mycovirus from the White Root Rot Fungus Rosellinia necatrix: Molecular and Biological Characterization, Taxonomic Considerations, and Potential for Biological Control. J. Virol. 2009, 83, 12801–12812. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Li, P.; Bhattacharjee, P.; Wang, S.; Zhang, L.; Ahmed, I.; Guo, L. Mycoviruses in Fusarium species: An update. Front. Cell. Infect. Microbiol. 2019, 9, 257. [Google Scholar] [CrossRef] [Green Version]
  26. Chu, Y.M.; Jeon, J.J.; Yea, S.J.; Kim, Y.H.; Yun, S.H.; Lee, Y.W.; Kim, K.H. Double-stranded RNA mycovirus from Fusarium graminearum. Appl. Environ. Microbiol. 2002, 68, 2529–2534. [Google Scholar] [CrossRef] [Green Version]
  27. Darissa, O.; Adam, G.; Schäfer, W. A dsRNA mycovirus causes hypovirulence of Fusarium graminearum to wheat and maize. Eur. J. Plant Pathol. 2012, 134, 181–189. [Google Scholar] [CrossRef]
  28. Li, P.; Zhang, H.; Chen, X.; Qiu, D.; Guo, L. Molecular characterization of a novel hypovirus from the plant pathogenic fungus Fusarium graminearum. Virology 2015, 481, 151–160. [Google Scholar] [CrossRef] [Green Version]
  29. Zhang, X.; Gao, F.; Zhang, F.; Xie, Y.; Zhou, L.; Yuan, H.; Zhang, S.; Li, H. The complete genomic sequence of a novel megabirnavirus from Fusarium pseudograminearum, the causal agent of wheat crown rot. Arch. Virol. 2018, 163, 3173–3175. [Google Scholar] [CrossRef]
  30. Xie, Y.; Wang, Z.; Li, K.; Liu, D.; Jia, Y.; Gao, F.; Dai, J.; Zhang, S.; Zhang, X.; Li, H. A megabirnavirus alleviates the pathogenicity of Fusarium pseudograminearum to wheat. Phytopathology 2021, 112, 1175–1184. [Google Scholar] [CrossRef]
  31. Li, P.; Lin, Y.; Zhang, H.; Wang, S.; Qiu, D.; Guo, L. Molecular characterization of a novel mycovirus of the family Tymoviridae isolated from the plant pathogenic fungus Fusarium graminearum. Virology 2016, 489, 86–94. [Google Scholar] [CrossRef] [Green Version]
  32. Cho, W.K.; Yu, J.; Lee, K.M.; Son, M.; Min, K.; Lee, Y.W.; Kim, K.H. Genome-wide expression profiling shows transcriptional reprogramming in Fusarium graminearum by Fusarium graminearum virus 1-DK21 infection. BMC Genom. 2012, 13, 173. [Google Scholar] [CrossRef] [Green Version]
  33. Lee, K.M.; Cho, W.K.; Yu, J.; Son, M.; Choi, H.; Min, K.; Lee, Y.W.; Kim, K.H. A comparison of transcriptional patterns and mycological phenotypes following infection of Fusarium graminearum by four mycoviruses. PLoS ONE 2014, 9, e0100989. [Google Scholar] [CrossRef] [PubMed]
  34. Wang, S.; Zhang, J.; Li, P.; Qiu, D.; Guo, L. Transcriptome-based discovery of Fusarium graminearum stress responses to FgHV1 infection. Int. J. Mol. Sci. 2016, 17, 1922. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Qu, Z.; Fu, Y.; Lin, Y.; Zhao, Z.; Zhang, X.; Cheng, J.; Xie, J.; Chen, T.; Li, B.; Jiang, D. Transcriptional Responses of Sclerotinia sclerotiorum to the Infection by SsHADV-1. J. Fungi 2021, 7, 493. [Google Scholar] [CrossRef] [PubMed]
  36. Chen, Y.; Kistler, H.C.; Ma, Z. Fusarium graminearum Trichothecene Mycotoxins: Biosynthesis, Regulation, and Management. Annu. Rev. Phytopathol. 2019, 57, 15–39. [Google Scholar] [CrossRef] [Green Version]
  37. Ding, S.W.; Li, H.; Lu, R.; Li, F.; Li, W.X. RNA silencing: A conserved antiviral immunity of plants and animals. Virus Res. 2004, 102, 109–115. [Google Scholar] [CrossRef]
  38. Nasmith, C.G.; Walkowiak, S.; Wang, L.; Leung, W.W.Y.; Gong, Y.; Johnston, A.; Harris, L.J.; Guttman, D.S.; Subramaniam, R. Tri6 is a global transcription regulator in the phytopathogen Fusarium graminearum. PLoS Pathog. 2011, 7, e1002266. [Google Scholar] [CrossRef] [PubMed]
  39. Proctor, R.H.; Hohn, T.M.; McCormick, S.P.; Desjardins, A.E. Tri6 encodes an unusual zinc finger protein involved in regulation of trichothecene biosynthesis in Fusarium sporotrichioides. Appl. Environ. Microbiol. 1995, 61, 1923–1930. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Peplow, A.W.; Tag, A.G.; Garifullina, G.F.; Beremand, M.N. Identification of new genes positively regulated by Tri10 and a regulatory network for trichothecene mycotoxin production. Appl. Environ. Microbiol. 2003, 69, 2731–2736. [Google Scholar] [CrossRef] [Green Version]
  41. Yun, Y.; Yin, D.; Dawood, D.H.; Liu, X.; Chen, Y.; Ma, Z. Functional characterization of FgERG3 and FgERG5 associated with ergosterol biosynthesis, vegetative differentiation and virulence of Fusarium graminearum. Fungal Genet. Biol. 2014, 68, 60–70. [Google Scholar] [CrossRef]
  42. Liu, X.; Jiang, J.; Yin, Y.; Ma, Z. Involvement of FgERG4 in ergosterol biosynthesis, vegetative differentiation and virulence in Fusarium graminearum. Mol. Plant Pathol. 2013, 14, 71–83. [Google Scholar] [CrossRef] [PubMed]
  43. Gilbert, J.; Abramson, D.; Mccallum, B.; Clear, R. Comparison of Canadian Fusarium graminearum isolates for aggressiveness, vegetative compatibility, and production of ergosterol and mycotoxins. Mycopathologia 2002, 153, 209–215. [Google Scholar] [CrossRef] [PubMed]
  44. Lamper, C.; Teren, J.; Bartok, T.; Komoroczy, R.; Mesterhazy, A.; Sagi, F. Predicting DON contamination in Fusarium-infected wheat grains via determination of the ergosterol content. Cereal Res. Commun. 2000, 28, 337–344. [Google Scholar] [CrossRef]
  45. Rodrigues, M.L. The multifunctional fungal ergosterol. mBio 2018, 9, e01755-18. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Jacob-Wilk, D.; Turina, M.; Van Alfen, N.K. Mycovirus Cryphonectria hypovirus 1 elements cofractionate with trans-Golgi network membranes of the fungal host Cryphonectria parasitica. J. Virol. 2006, 80, 6588–6596. [Google Scholar] [CrossRef] [Green Version]
  47. Van De Sande, W.W.J.; Lo-Ten-Foe, J.R.; Van Belkum, A.; Netea, M.G.; Kullberg, B.J.; Vonk, A.G. Mycoviruses: Future therapeutic agents of invasive fungal infections in humans? Eur. J. Clin. Microbiol. Infect. Dis. 2010, 29, 755–763. [Google Scholar] [CrossRef] [Green Version]
  48. Seong, K.; Li, L.; Hou, Z.; Tracy, M.; Kistler, H.C.; Xu, J.R. Cryptic promoter activity in the coding region of the HMG-CoA reductase gene in Fusarium graminearum. Fungal Genet. Biol. 2006, 43, 34–41. [Google Scholar] [CrossRef]
  49. Stermer, B.A.; Bianchini, G.M.; Korth, K.L. Regulation of HMG-CoA reductase activity in plants. J. Lipid Res. 1994, 35, 1133–1140. [Google Scholar] [CrossRef]
  50. Lee, J.; Myong, K.; Kim, J.E.; Kim, H.K.; Yun, S.H.; Lee, Y.W. FgVelB globally regulates sexual reproduction, mycotoxin production and pathogenicity in the cereal pathogen Fusarium graminearum. Microbiology 2012, 158, 1723–1733. [Google Scholar] [CrossRef] [Green Version]
  51. Jiang, J.; Yun, Y.; Liu, Y.; Ma, Z. FgVeLB is associated with vegetative differentiation, secondary metabolism and virulence in Fusarium graminearum. Fungal Genet. Biol. 2012, 49, 653–662. [Google Scholar] [CrossRef]
  52. Chun, J.; Ko, Y.H.; Kim, D.H. Transcriptome analysis of Cryphonectria parasitica infected with Cryphonectria hypovirus 1 (CHV1) reveals distinct genes related to fungal metabolites, virulence, antiviral RNA-silencing, and their regulation. Front. Microbiol. 2020, 11, 1711. [Google Scholar] [CrossRef] [PubMed]
  53. Subramaniam, R.; Narayanan, S.; Walkowiak, S.; Wang, L.; Joshi, M.; Rocheleau, H.; Ouellet, T.; Harris, L.J. Leucine metabolism regulates TRI6 expression and affects deoxynivalenol production and virulence in Fusarium graminearum. Mol. Microbiol. 2015, 98, 760–769. [Google Scholar] [CrossRef] [PubMed]
  54. Leonetti, P.; Stuttmann, J.; Pantaleo, V. Regulation of plant antiviral defense genes via host RNA-silencing mechanisms. Virol. J. 2021, 18, 194. [Google Scholar] [CrossRef] [PubMed]
  55. Nuss, D.L. Mycoviruses, RNA Silencing, and Viral RNA Recombination. Adv. Virus Res. 2011, 80, 25–48. [Google Scholar] [CrossRef] [Green Version]
  56. Mochama, P.; Jadhav, P.; Neupane, A.; Marzano, S.Y.L. Mycoviruses as triggers and targets of RNA silencing in white mold fungus Sclerotinia sclerotiorum. Viruses 2018, 10, 214. [Google Scholar] [CrossRef] [Green Version]
  57. Yu, J.; Park, J.Y.; Heo, J.I.; Kim, K.H. The ORF2 protein of Fusarium graminearum virus 1 suppresses the transcription of FgDICER2 and FgAGO1 to limit host antiviral defences. Mol. Plant Pathol. 2020, 21, 230–243. [Google Scholar] [CrossRef] [Green Version]
Toxins 14 00503 g001 550
Figure 1. HPLC-MS/MS chromatograms showing DON production in wheat grain culture extracts of F. pseudograminearum strain FC136-2A harboring FpgMBV1 and the isogenic virus-free strain FC136-2A-V-.
Figure 1. HPLC-MS/MS chromatograms showing DON production in wheat grain culture extracts of F. pseudograminearum strain FC136-2A harboring FpgMBV1 and the isogenic virus-free strain FC136-2A-V-.
Toxins 14 00503 g001
Toxins 14 00503 g002 550
Figure 2. Effects of FpgMBV1 on the transcriptome profile of F. pseudograminearum as revealed by RNA-Seq. (A) The volcano plot showing gene signals detected in the strain FC136-2A comparing to FC136-2A-V-. Up-regulated (red), no change (gray) and down-regulated (blue). (B) The histogram of the number of differentially expressed genes (DEGs), including numbers of up-regulated (red) and down-regulated (blue) genes. (C) The number of DEGs annotated in the PHI database. (D) The rate of numbers of DEGs annotated by Gene ontology terms in molecular function (MF) (green), cellular component (CC) (red) and biological process (BP) (blue).
Figure 2. Effects of FpgMBV1 on the transcriptome profile of F. pseudograminearum as revealed by RNA-Seq. (A) The volcano plot showing gene signals detected in the strain FC136-2A comparing to FC136-2A-V-. Up-regulated (red), no change (gray) and down-regulated (blue). (B) The histogram of the number of differentially expressed genes (DEGs), including numbers of up-regulated (red) and down-regulated (blue) genes. (C) The number of DEGs annotated in the PHI database. (D) The rate of numbers of DEGs annotated by Gene ontology terms in molecular function (MF) (green), cellular component (CC) (red) and biological process (BP) (blue).
Toxins 14 00503 g002
Toxins 14 00503 g003 550
Figure 3. The Gene Ontology enrichment analysis of differently expressed genes, including up-regulated genes (A) and down-regulated genes (B) and the KEGG pathway enrichment analysis of the up-regulated genes (C) and down-regulated genes (D).
Figure 3. The Gene Ontology enrichment analysis of differently expressed genes, including up-regulated genes (A) and down-regulated genes (B) and the KEGG pathway enrichment analysis of the up-regulated genes (C) and down-regulated genes (D).
Toxins 14 00503 g003
Toxins 14 00503 g004 550
Figure 4. Heat map of metabolism-related genes differentially expressed in F. pseudograminearum strain FC136-2A compared to FC136-2A-V-.
Figure 4. Heat map of metabolism-related genes differentially expressed in F. pseudograminearum strain FC136-2A compared to FC136-2A-V-.
Toxins 14 00503 g004
Toxins 14 00503 g005 550
Figure 5. Heat map of TRI genes differentially expressed in F. pseudograminearum strain FC136-2A compared to FC136-2A-V-.
Figure 5. Heat map of TRI genes differentially expressed in F. pseudograminearum strain FC136-2A compared to FC136-2A-V-.
Toxins 14 00503 g005
Toxins 14 00503 g006 550
Figure 6. Heat map of genes related to ergosterol biosynthesis and metabolism differentially expressed in F. pseudograminearum strain FC136-2A compared to FC136-2A-V-, including genes annotated as ergosterol biosynthesis and metabolism based on NCBI-NR database (A), and based on PHI database (B).
Figure 6. Heat map of genes related to ergosterol biosynthesis and metabolism differentially expressed in F. pseudograminearum strain FC136-2A compared to FC136-2A-V-, including genes annotated as ergosterol biosynthesis and metabolism based on NCBI-NR database (A), and based on PHI database (B).
Toxins 14 00503 g006
Toxins 14 00503 g007 550
Figure 7. Heat map of RNA-silencing-related genes differentially expressed between F. pseudograminearum strain FC136-2A harboring FpgMBV1 and the isogenic virus-free strain FC136-2A-V-.
Figure 7. Heat map of RNA-silencing-related genes differentially expressed between F. pseudograminearum strain FC136-2A harboring FpgMBV1 and the isogenic virus-free strain FC136-2A-V-.
Toxins 14 00503 g007
Toxins 14 00503 g008 550
Figure 8. Gene expression comparison of some genes between F. pseudograminearum strain FC136-2A and FC136-2A-V- by quantitative real-time RT-PCR (orange) and RNA-seq (blue).
Figure 8. Gene expression comparison of some genes between F. pseudograminearum strain FC136-2A and FC136-2A-V- by quantitative real-time RT-PCR (orange) and RNA-seq (blue).
Toxins 14 00503 g008
Toxins 14 00503 g009 550
Figure 9. General view of genes and pathways regulated by FpgMBV1 in F. pseudograminearum revealed in this study.
Figure 9. General view of genes and pathways regulated by FpgMBV1 in F. pseudograminearum revealed in this study.
Toxins 14 00503 g009
Table
Table 1. DON and ergosterol biosynthesis and metabolism-related genes according to the PHI database.
Table 1. DON and ergosterol biosynthesis and metabolism-related genes according to the PHI database.
Gene NameGene IDSwissProt_Descriptionlog2(FC136_2A/FC136_2A_V-)Identity (%)E-ValuePHI
Accession
ERG3BFPSE_12291Probable Delta(7)-sterol 5(6)-desaturase−1.65692.37.3 × 10−179PHI:3036
ERG4FPSE_03901Delta(24(24(1)))-sterol reductase−1.33696.60PHI:2728
ERG5AFPSE_09181Cytochrome P450 61−1.58499.30PHI:3037
ERG5BFPSE_01847Sterol 22-desaturase−7.32094.31.6 × 10−292PHI:3038
TRI5FPSE_12160Trichodiene synthase−3.0301001.1 × 10−228PHI:6846
TRI6FPSE_12161Trichothecene biosynthesis transcription regulator 6−2.68695.91.2 × 10−126PHI:1362
TRI12FPSE_12157Trichothecene efflux pump TRI12−2.30591.74.5 × 10−307PHI:2704
TRI10FPSE_12159Trichothecene biosynthesis transcription regulator 10−2.70293.31.5 × 10−232PHI:2328
TRI15FPSE_02457Cys(2)-His(2) zinc finger protein−1.80395.78.3 × 10−183PHI:1363
HMR1FPSE_03466Hydroxymethylglutaryl CoA reductase gene−1.78194.70PHI:1006
VELBFPSE_11531Velvet complex subunit B−1.76083.31.1 × 10−159PHI:2427
GLXFPSE_04483WSC domain-containing protein ARB_07867−1.238870PHI:5393
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, K.; Liu, D.; Pan, X.; Yan, S.; Song, J.; Liu, D.; Wang, Z.; Xie, Y.; Dai, J.; Liu, J.; et al. Deoxynivalenol Biosynthesis in Fusarium pseudograminearum Significantly Repressed by a Megabirnavirus. Toxins 2022, 14, 503. https://doi.org/10.3390/toxins14070503

AMA Style

Li K, Liu D, Pan X, Yan S, Song J, Liu D, Wang Z, Xie Y, Dai J, Liu J, et al. Deoxynivalenol Biosynthesis in Fusarium pseudograminearum Significantly Repressed by a Megabirnavirus. Toxins. 2022; 14(7):503. https://doi.org/10.3390/toxins14070503

Chicago/Turabian Style

Li, Ke, Dongmei Liu, Xin Pan, Shuwei Yan, Jiaqing Song, Dongwei Liu, Zhifang Wang, Yuan Xie, Junli Dai, Jihong Liu, and et al. 2022. "Deoxynivalenol Biosynthesis in Fusarium pseudograminearum Significantly Repressed by a Megabirnavirus" Toxins 14, no. 7: 503. https://doi.org/10.3390/toxins14070503

APA Style

Li, K., Liu, D., Pan, X., Yan, S., Song, J., Liu, D., Wang, Z., Xie, Y., Dai, J., Liu, J., Li, H., Zhang, X., & Gao, F. (2022). Deoxynivalenol Biosynthesis in Fusarium pseudograminearum Significantly Repressed by a Megabirnavirus. Toxins, 14(7), 503. https://doi.org/10.3390/toxins14070503

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop