Next Article in Journal
Combining the Peptide RWQWRWQWR and an Ethanolic Extract of Bidens pilosa Enhances the Activity against Sensitive and Resistant Candida albicans and C. auris Strains
Previous Article in Journal
In Vivo Microevolutionary Analysis of a Fatal Case of Rhinofacial and Disseminated Mycosis Due to Azole-Drug-Resistant Candida Species
Previous Article in Special Issue
Multiple Clonostachys rosea UDP-Glycosyltransferases Contribute to the Production of 15-Acetyl-Deoxynivalenol-3-O-Glycoside When Confronted with Fusarium graminearum
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Epistatic Relationship between MGV1 and TRI6 in the Regulation of Biosynthetic Gene Clusters in Fusarium graminearum

by
Kristina Shostak
1,†,
Dianevys González-Peña Fundora
2,3,†,
Christopher Blackman
1,4,
Tom Witte
1,5,
Amanda Sproule
1,
David Overy
1,
Anas Eranthodi
2,
Nehal Thakor
3,
Nora A. Foroud
2,* and
Rajagopal Subramaniam
1,4,*
1
Ottawa Research and Development Centre, Agriculture and Agri-Food Canada, Ottawa, ON K1A 0C6, Canada
2
Lethbridge Research and Development Centre, Agriculture and Agri-Food Canada, Lethbridge, AB T1J 4B1, Canada
3
Department of Chemistry and Biochemistry, University of Lethbridge, Lethbridge, AB T1K 4M4, Canada
4
Department of Cell and System Biology, University of Toronto, Toronto, ON M5S 3B2, Canada
5
Department of Biology, University of Ottawa, Ottawa, ON K1N 6N5, Canada
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
J. Fungi 2023, 9(8), 816; https://doi.org/10.3390/jof9080816
Submission received: 1 July 2023 / Revised: 28 July 2023 / Accepted: 31 July 2023 / Published: 2 August 2023
(This article belongs to the Special Issue Genome Editing Tools in Fungi)

Abstract

:
Genetic studies have shown that the MAP kinase MGV1 and the transcriptional regulator TRI6 regulate many of the same biosynthetic gene clusters (BGCs) in Fusarium graminearum. This study sought to investigate the relationship between MGV1 and TRI6 in the regulatory hierarchy. Transgenic F. graminearum strains constitutively expressing MGV1 and TRI6 were generated to address both independent and epistatic regulation of BGCs by MGV1 and TRI6. We performed a comparative transcriptome analysis between axenic cultures grown in nutrient-rich and secondary metabolite-inducing conditions. The results indicated that BGCs regulated independently by Mgv1 included genes of BGC52, whereas genes uniquely regulated by TRI6 included the gene cluster (BGC49) that produces gramillin. To understand the epistatic relationship between MGV1 and TRI6, CRISPR/Cas9 was used to insert a constitutive promoter to drive TRI6 expression in the Δmgv1 strain. The results indicate that BGCs that produce deoxynivalenol and fusaoctaxin are co-regulated, with TRI6 being partially regulated by MGV1. Overall, the findings from this study indicate that MGV1 provides an articulation point to differentially regulate various BGCs. Moreover, TRI6, embedded in one of the BGCs provides specificity to regulate the expression of the genes in the BGC.

1. Introduction

The regulatory hierarchy governing the expression of BGCs has been limited to a few genetic pathways in F. graminearum [1,2,3,4,5]. Studies have predominantly focused on the regulation of the BGC that produces the trichothecene mycotoxin deoxynivalenol (DON). The BGC23 harbors the majority of the genes involved in the production of DON, including the two regulatory genes TRI6 and TRI10 [6]. Targeted disruption of TRI6 and TRI10 significantly reduced the expression of genes within the cluster and abolished DON production [7]. Chromatin immunoprecipitation followed by Illumina sequencing (ChIP-seq) revealed that TRI6 is a global transcription factor that regulates more than 200 additional targets involved in carbohydrate metabolism and signal transduction [8]. Furthermore, TRI6 autoregulates its expression under nutrient-rich conditions and binds to promoter elements of some of the genes in the cluster under nutrient-limiting conditions. TRI10 encodes a protein with no known functional domains, and unlike in F. sporotrichioides, TRI6 expression levels were not significantly reduced in the ΔTRI10 mutant in F. graminearum PH-1, although deletion of TRI10 significantly reduced DON production [5].
Environmental and physiological cues, such as nutrient sources, light, and pH of the culture media, play an important role in the regulation of BGCs, and several common signal transducers and transcriptional regulators of many fungi have been characterized [9]. As an example, the target of the rapamycin (TOR) complex, conserved in all eukaryotes, integrates many of the nutrient cues and in two species of Fusarium (F. fujikuroi and F. graminearum); members of the TOR complex define the specificity to the expression of BGC genes [3,10]. The signals emanating from the TOR complexes are amplified by the mitogen-activated protein kinase modules (MAPK) [11]. The three MAPK pathways present in Fusarium graminearum have been ascribed distinct roles, with the HOG1 MAPK responsible for sensing and responding to the osmotic stress, while the MAP kinase module, MGV1, is linked to both monitoring and maintenance of the cell wall integrity (CWI) pathway [12,13]. An orthologue of MGV1, GPMK1 represents the third MAPK pathway that is responsible for filamentous growth and sexual reproduction [14,15]. All three of these pathways are important for vegetative growth, reproduction, and pathogenicity in F. graminearum; however, evidence suggests that only HOG1 and MGV1 are major contributors to the regulation of BGC in F. graminearum [12]. Although targeted deletion of HOG1 and MGV1 shows reduced virulence, a phenotype related to reduced fitness, the impact on the production of secondary metabolites is distinct [12]. Both exhibit contrasting phenotypes in the production of DON; genetic studies show that MGV1 positively impacts DON production, while HOG1 has the opposite effect. Similar contrasting phenotypes have been observed with aurofusarin and zearalenone biosynthesis [12]. It is not currently known how the distinct regulations are achieved, except to suggest that many of the targets of the various MAPKs are regulated post-translationally and, thus, complex regulatory mechanisms are in play.
The signals from the MAPKs and other signaling modules are parsed through many transcriptional regulators that act either independently, synergistically, or epistatically to the various signaling modules [16,17]. One of the well-characterized families of transcription regulators, AreA, plays a pivotal role in integrating nitrogen metabolism [18]. AreA is responsive to nitrogen sources; regulated by the cyclic adenosine monophosphate (cAMP)—protein kinase A (PKA) pathway, it has been shown to interact with TRI10 and regulate the DON pathway in F. graminearum [5]. The conspicuous absence of regulation in zearalenone production by AreA underscores the distinctiveness of the regulation of BGCs by the various signaling modules. Additional support for this distinctiveness comes from studies of other transcriptional regulators, such as the velvet complex proteins that perceive light signals and PacC, which perceive changes in pH in the culture medium to regulate an array of BGCs, including the biosynthesis of DON [19,20].
To distinguish the regulation of BGCs by the MGV1 signaling pathway and the transcriptional regulator TRI6 in F. graminearum, we undertook comparative transcriptomic studies with deletion and constitutive expression of MGV1 and TRI6. Our results show that MGV1 broadly regulates many BGCs, while TRI6 distinguishes the BGC’s regulation towards the production of various metabolites.

2. Materials and Methods

2.1. F. graminearum Strains and Culture Conditions

F. graminearum (NRRL29169) was the wild-type (WT) and parental strain used for the construction of the transgenic strains. The construction of ∆mgv1 and MGV1Oex6 strains have been previously described [21,22]. Briefly, overexpression was achieved by replacing the native promoter of MGV1 with the constitutive promoter (pGPD) of the glyceraldehyde 3-phosphate dehydrogenase gene from A. nidulans [22]. The construction of other MGV1Oex lines was constructed similar to MGV1Oex6 [22,23]. The macroconidia were produced in carboxymethylcellulose media and used as the inoculum for all cultures at a concentration of 5 × 103 spores mL−1. Mycelia were grown in either preferred nutrient medium (PN) (56 mM NH4Cl, 8.1 mM MgSO4 7H2O, 0.23 mM FeSO4·7H2O, 14.7 mM KH2PO4, 2 g L−1 peptone, 2 g L−1 yeast extract, 2 g L−1 malt extract, and 111 mM glucose) or in nonpreferred nutrient (NPN) medium (6.2 mM putrescine dihydrochloride, 22 mM KH2PO4, 0.8 mM MgSO4·7H2O, 85.6 mM NaCl, 116.8 mM sucrose, 108.6 mM glycerol, pH 4.0). Unless otherwise specified, all strains in the study were cultured in the dark, at 28 °C, with shaking at 160 rpm for 24 h. For nonpreferred nutrient (NPN) conditions, strains were first cultured in PN for 24 h, as described, followed by washing and growth in a DON-producing medium, as previously described [24].

2.2. Construction of the Δmgv1/TRI6OX Strain by CRISPR-Cas9

Homologous-based recombination (HDR) using CRISPR/Cas9 in F. graminearum was based on previous studies [25,26]. Briefly, microhomology primers flanking MGV1 were used for the construction of a hygromycin repair template to replace the MGV1 gene (Figure 1A and Figure S1). Another set of microhomology primers flanking a 225 bp portion of the TRI6 promoter was used to construct the repair template containing 1016 bp of GAPDH promoter to functionally replace the native TRI6 promoter (Figure 1B and Figure S2). The upstream and downstream gRNAs associated with each genomic modification and 4 µg of each of the repair templates were incubated with 100 µL of protoplasts. Protoplast preparation, regeneration, and transformant selection were performed as before [24,25]. Deletion of the mgv1 gene and the TRI6 promoter replacement was confirmed using whole genome sequencing (Figure S3). The transformants were sequenced with Illumina Novaseq Shotgun DNA sequencing (https://www.genomequebec.com, accessed on 12 September 2021). Sequencing reads were aligned to the reference genome of NRRL29169 in CLC-Genomics Workbench (v20) using the default settings and examined for large InDels and the absence of the gene of interest [25]. Gene expression was confirmed using reverse transcriptase quantitative PCR (RT-qPCR) under PN conditions (Figure S4).

2.3. Gene Expression Analysis by RNA Sequencing

F. graminearum wild-type, ∆mgv1, MGV1Oex1, MGV1Oex6, and ∆mgv1::TRI6OX strains were grown in 4 mL of liquid culture under NPN conditions as described [24]. Total RNA from three biological replicates per strain was extracted using Trizol reagent, as described previously [8,25]. Total RNA was sequenced using Illumina NovaSeq 6000 platform (https://www.genomequebec.com accessed on 12 September 2021). The raw reads were trimmed and aligned to the F. graminearum gene coding sequences predicted in NRRL29169 assembly (GWAS Accession# SPRZ00000000) using the RNA-Seq Analysis feature in CLC Genomics Workbench, version 12, with the following parameters: mismatch cost 2, indel cost 3, length fraction 0.9, similarity fraction 0.8, the maximum number of hits for a read 10. Unique gene read counts were manually curated for genes with a minimum of 50 unique gene reads. Filtered counts were then imported into the R environment and normalized using default parameters in DESeq2, which uses negative binomial models to represent the number of reads assigned to a gene in a sample. Differential expression analysis was also carried out using DESeq2, with the threshold fold change ≥2 and the FDR-adjusted p-value ≤ 0.05. The normalized reads and differential expression genes (DEGs) are included in the Supplementary Materials (Tables S1 and S4).

2.4. Quantitative PCR (RT-qPCR) Analyses

The gene expression from the RNA-sequencing (RNA-seq) data was confirmed using RT-qPCR. One microgram of total RNA was reverse-transcribed using a high-capacity cDNA reverse transcription kit with random hexamer primers (Applied Biosystems, Waltham, MA, USA). The qPCR was performed using the Applied Biosystems PowerUp SYBR Green reaction mix and QuantStudio 3 Real-Time PCR system (Thermo Fisher, Nepean, ON, Canada). The relative expression was calculated using the Pfaffl method with EF1α (FGSG_08811) and β-tubulin (FGSG_09530) as reference genes [27]. Significance was determined using the Student’s t-test at p ≤ 0.05. The primers used for the amplification are listed in Table S5.

2.5. Mycelial Growth and Reproductive Structures Assays

Vegetative growth, spore germination, and perithecia formation were measured for each of the MGV1Oex transformants and the WT. For vegetative growth, the fungus was grown on potato dextrose agar (PDA) plates with a mycelial plug as the inoculum. The growth diameter was measured every 24 h until the mycelia reached the edge of the plate. Eight plates were evaluated for each strain, and the experiment was repeated twice (n = 16/experiment).
Macroconidia spore germination was measured on synthetic nutrient agar (Spezieller Nährstoffarmer Agar; SNA) at 27 °C. Each plate was marked with 0.7 × 0.7 cm2, and 8 squares (replicates) were inoculated with 3 μL of 104 macroconidia mL−1. Germination was observed in each square using a light microscope (Leica DM500, Concord, ON, Canada) at 100× magnification. At 3, 6, and 8 h, the number of germinated spores (when the germ tube exceeded the size of the conidia) was counted under a light microscope (Leica DM500) at 100× magnification. The experiment was performed in triplicate (n = 24). Data were reported as the percent macroconidia germination and analyzed with one-way ANOVA (analysis of variance), and the means were compared using the Tukey test, p ≤ 0.05.
Perithecia formation was determined on carrot agar inoculated with mycelial plugs, and after 4 d of incubation at 28 °C, the mycelia were removed. Perithecia formation was induced with the application of 1.0 mL (v/v) of 2.5% Tween 60 to the agar surface, and the plates were incubated at 25 °C with continuous light; this step was repeated if after 48 h the mycelia reappeared. The presence of perithecia was determined after 3–5 d by disrupting the perithecium with a cover slide to release the ascospores, which were then observed under a light microscope (Leica DM500) at 100× magnification. The experiment was conducted in triplicate.

2.6. Quantification of 15-Acetylated DON (15-ADON) Accumulation in Axenic Cultures

The 15-ADON accumulation analyses in Δmgv1, MGV1Oex strains, and WT were performed using a two-stage liquid media protocol [24]. A suspension of 10 μL containing 2 × 104 spores was inoculated into a Falcon Multiwell 6-well culture tray containing 4 mL of first-stage growth media at pH 7.0. Each well contained an autoclaved nylon net filter (100 μm NY1H type, Millipore, Burlington, MA, USA). The trays were sealed with parafilm and affixed to an orbital shaker, for 24 h in the dark, at 170 rpm. After 24 h, the liquid was removed using a sterile transfer pipette without disturbing the filter with the mycelial growth, and it was resuspended in sterile water to remove traces of first-stage media. This step was repeated, and the mycelia/filter was then resuspended in 4 mL of second-stage media, pH 4.0 [8]. The mycelium was grown in the second stage under the same conditions described above, and the supernatant was collected after 48 h.
Mycelial solids were collected and dried under a vacuum for 24 h before weighing. A nylon filter was used as a blank measurement. The sample supernatant was first filtered (0.2 μm) and combined 450 μL sample + 150 μL MeOH in the HPLC vial. Trichothecenes were analyzed on a Shimadzu prominence LC-20AD (Mandel) with 100 μL injection on a Shimadzu SIL-20A HT prominence autosampler. The samples were run on a Restek Pinnacle DB C18 Column (5 um, 150 × 4.6 mm) using a 22.5% isocratic MeOH: H2O flow over 20 min at a rate of 1 mL min−1. Trichothecenes were monitored by UV 220 nm [8,23].

2.7. Immunodetection of MAPKs

Protein was extracted from the mycelia of Δmgv1, MGV1Oex, and WT strains collected from four-day cultures on PDA of or from liquid culture under NPN conditions as described above. For each strain and condition, samples were collected in triplicate for one protein extraction, and the experiment was repeated three times. Mycelia were ground to a powder in liquid nitrogen using a mortar and pestle and transferred to a tube with 10 mL MAPK extraction buffer per gram of mycelia. MAPK extraction buffer was composed of 1 M NaCl, 50 mM NaF, 1 mM phenylmethylsulfonyl fluoride (PMSF), 1 mM Na3VO4, protease inhibitor cocktail (Thermofisher, Halt™ Protease Inhibitor Cocktail, 100 µL per 10 mL buffer), 0.2% β-mercaptoethanol (BME), and 50 mM sodium phosphate buffer, pH 8. The samples were incubated at 4 °C with gentle agitation for 1 h. The cell debris was removed by centrifugation in 1.5 mL tubes at 13,751× g for 5 min at 4 °C. The supernatant was filtered through a 0.45 µm syringe filter. The protein concentration was measured using the Bradford reagent (BioRad, USA).
Ten micrograms of total protein were denatured with 1 × NuPAGETM LDS Sample Buffer (ThermoFisher) and 0.5% (v/v) β-mercaptoethanol at 70 °C for 10 min. The samples were resolved on 12-well BoltTM 4–12% Bis-Tris-Plus gels (Invitrogen) in BoltTM MOPS SDS Running Buffer (Invitrogen) at 150 V for 1.5 h. The gels were run in duplicate, where one was stained with OmniPur Coomassie Blue R-250 (EMD Millipore) to visualize the resolved proteins, and the other was used to transfer proteins onto Immun-Blot® PVDF Membranes (BioRad, Saint-Laurent, QC, Canada). Protein transfer was carried out by electrophoresis in 25 mM Tris, 192 mM glycine, pH-8.3, and 20% methanol at 100 V for 1.5 h at room temperature. The uniform transfer was verified with the PierceTM Reversible Protein Stain Kit for PVDF Membranes, according to the manufacturer’s instructions. After stain removal, membranes were blocked with 3% bovine serum albumin in Tris-buffered saline (TBS) for 1 h and then probed overnight at 4 °C with 10 mL of 1:1000 dilutions of primary antibody in 5% BSA in TBST (1x TBS buffer with 0.05% Tween20). Four different primary polyclonal antibodies were used to detect different MAPKs: (1) anti-phospho-p44/42 MAPK (ERK1/2) (Cell Signaling Technology, Danvers, MA. USA) detects phosphorylated forms of the Fusarium MAPKs MGV1 (p-MGV1) and GPMK1 (p-GPMK1); (2) anti-MAP kinase (ERK1/2) (Sigma) detects total MGV1 and GPMK1; (3) anti-phospho-P38 MAPK (Cell Signaling Technology, USA) detects phosphorylated forms of the Fusarium p38 MAPK, p-HOG1; (4) anti-P38 MAPK (Cell Signal, USA) detects unphosphorylated forms of HOG1. Anti-α tubulin YOL 1/34 (Santa Cruz Biotechnology) from rats and anti-actin (Abcam, Waltham, MA, USA) from mice were used to detect the Fusarium housekeeping protein tubulin and actin, respectively. Following three washes in TBST, the membrane was incubated for 1 h at room temperature with 12 mL of Goat anti-rabbit-HRP, Goat anti-rat-HRP, or Goat anti-mouse-HRP conjugated secondary antibodies (BioRad), and 1:10,000 Precision ProteinTM StrepTactin-HRP conjugate (BioRad) in 3% BSA-TBST. The blots were developed with SuperSignalTM West Pico PLUS Chemiluminescent Substrate.

2.8. FHB Disease Assays

FHB assays were carried out in Canadian wheat cultivars, ‘Roblin’ (highly susceptible) and ‘Penhold’ (moderately resistant). Plants were grown at 22 °C/18 °C (day/night) in a greenhouse with a 16 h photoperiod. At anthesis, one spike per plant was point inoculated by pipetting 10 µL of inoculum (1 × 105 macroconidia mL−1 in 0.2% Tween20) into a single floret. After inoculation, plants were incubated in a 25 °C (day/night) misting growth cabinet (95% relative humidity) with a 16 h photoperiod for 3 days and then returned to their original growing environments. Disease rating was performed at 7-, 12-, and 18-days post inoculation (dpi) by counting the number of diseased spikelets below and including the inoculated spikelet. The experiment was repeated three times over four weeks. The mean values of five plants per inoculum were calculated for each repetition using one-way ANOVA followed by Tukey’s HSD test (α = 0.05).

3. Results

The NPN (nonpreferred nutrient) conditions have previously been documented to activate a plethora of BGCs, including BGC23, which is responsible for DON production (24). An outline of the scheme is presented to navigate the various RNA-seq datasets that were used to distinguish the expression of BGCs influenced by TRI6 and MGV1 under NPN conditions (Figure 2). The NPN acts as a stimulus to trigger the MAPK, MGV1, which can activate BGCs (Pathways A and B, Figure 2) or indirectly through the action of TRI6 (Pathway B, Figure 2). The stimulus can also activate BGCs that are independent of MGV1, but dependent on TRI6 (Pathway C, Figure 2). Finally, the BGCs can be co-regulated by both MGV1 and TRI6 (Pathway D, Figure 2)
To obtain the genes specifically regulated by MGV1, we constructed four MGV1Oex strains (Oex1, Oex2; Oex6, and Oex 8), and under noninducing (or PN) conditions, all four strains exhibited constitutive expression and phosphorylation of MGV1 compared to the WT strain (Figure 3A,B). In response to stimulus (NPN condition), MGV1 is activated in both the WT and MGV1Oex6, and these conditions do not alter the expression of the other two F. graminearum MAPKs, Gpmk1 and Hog1 (Figure 3C). The four strains with active MGV1 exhibited similar growth properties as the WT (mycelia and germination; Figure S5A,B). In contrast, all four strains produced higher levels of 15-ADON in the culture medium (Figure S6). Interestingly, the increased levels of 15-ADON in the four strains did not translate into an increase in disease spread in either a susceptible wheat variety (cv. Roblin) or an intermediate resistance wheat variety (cv. Penhold) (Figure S7A,B).
We performed RNA-seq analysis with WT, Δmgv1, and the two MGV1Oex strains (Oex1 and Oex6) grown in NPN conditions. Normalized reads and DEGs are shown in Tables S1 and S2. To obtain MGV1-specific regulated genes, we developed a common dataset of differentially regulated genes (DEGs) of the two MGV1Oex strains and the Δmgv1 strain (Figure 4A). The overlap dataset with 1120 DEGs represented MGV1-regulated genes of the two Δmgv1-MGV1Oex strains (1207 in MGV1Oex1 and 1185 in MGV1Oex6) (Figure 4B, Table S2). A FunCat analysis of 1120 DEGs encompassed genes from both primary and secondary metabolisms, cellular transport, and genes involved in virulence (Table S2). The RNA-seq was validated by RT-qPCR experiments with six genes that were either positively or negatively regulated by MGV1 (Figure 5).
The 1120 DEGs comprised 107 BGC genes (Figure 6A; Table S2). These included positive regulation of BGCs 23, 49, and 64 by MGV1, involved in the production of 15-ADON, butanolide, and fusaoctaxin, respectively, while aurofusarin (BGC13) and fusarin (BGC42) and BGCs 06, 09, 12, 37, 51, 52, and 62 with no known identified products were negatively regulated by MGV1 (Table S2). It is noteworthy to indicate that there was a 5-fold downregulation of TRI6 and concomitant reduction in TRI5 expression in the Δmgv1 mutant strain, suggesting that MGV1 partially regulates the expression of TRI6, which will impact the expression of BGC23 and other BGCs regulated by TRI6 (Figure 5) [24].
Similarly, to obtain the genes specifically regulated by TRI6, a dataset of DEGs from the RNA-seq analysis between Δtri6 and WT grown in NPN conditions was compared to the dataset of DEGs from RNA-seq analysis between TRI6OX and the ΔTri6 strain grown under the same conditions. This dataset has been published and the analysis of the 687 DEGs will only be briefly summarized to suit this study (Figure 6A, Table S3) [24]. We observed 110 BGCs that are differentially expressed (Figure 6A; Table S3). Similar to MGV1, TRI6 positively regulates BGCs that produce 15-ADON, butanolide, and fusaoctaxin and negatively influences the production of gramillin, aurofusarin, and fusarin (Table S3).
To obtain genes co-regulated by MGV1 and TRI6, we compared the 1120 DEGs regulated by MGV1 and 687 DEGs regulated by TRI6 (Figure 6, Table S3). The overlap of two datasets represented 228 co-regulated genes, representing 51 BGCs that included the positively regulated trichothecene gene cluster (BGC23), butanolide gene cluster (BGC49), and fusaoctaxin (BGC64), the negatively regulated aurofusarin gene cluster (BGC13). In addition to genes that are co-regulated, we also identified 892 DEGs with 56 BGCs that are MGV1 dependent but TRI6 independent (Table S3). This category contained only ‘orphan’ genes of a BGC; we defined ‘orphan’ genes as those that are part of the BGC but have different co-expression patterns. For example, only three (FGSG_03529, FGSG_03531, FGSG_03533) of the ~15 genes of the trichothecene gene cluster (BGC23) are represented in this category and their biochemical function in the cluster have not been properly elucidated. These ‘orphan’ genes are likely regulated independently of the MGV1 pathway. The Venn diagram also showed 459 DEGs comprising 59 BGCs that are regulated independently of MGV1 but are TRI6 dependent (Table S3). The majority of genes included in this category are represented in BGC02, (responsible for gramillin production), BGC16, and BGC66 (chrysogine production).
Since epistasis between MGV1 and TRI6 was alluded to, we were interested to differentiate the regulation of BGCs by the two genes. Therefore, we constitutively overexpressed TRI6 in locus in the Δmgv1 mutant strain background (Δmgv1::TRI6Oex). The strain was constructed by simultaneous insertion of the hygromycin resistance marker in the MGV1 locus and a ~1 kb promoter of glyceraldehyde 3-phosphate dehydrogenase (pGPD), a constitutive promoter, upstream of the TRI6 start site. CRISPR/Cas9 was used to make the double insertion construct, which was verified by both PCR and whole genome sequencing (Figure S4); we used RT-qPCR to confirm that the expression TRI6 was constitutive (Figure S5). We characterized this strain with respect to growth properties in three different media and, as observed, the constitutive expression of TRI6 did not significantly alter the phenotype observed with that from the Δmgv1 strain (Figure S8). This demonstrates that the fitness loss linked to MGV1 disruption cannot be rescued by overexpression of TRI6 and, therefore, is not linked to TRI6 [13].
An RNA-seq experiment was performed with the Δmgv1::TRI6Oex strain in the NPN condition and was compared to the WT strain grown in the same conditions. A total of 969 DEGs represented genes regulated by TRI6 that were independent of MGV1 (Table S4). Among these DEGs, 108 represented genes of the BGCs involved in the production of gramillin (BGC02), aurofusarin (BGC13), orcinol (BGC18), fusarin (BGC42), fusaristatin (BGC47), butanolide (BGC49), and chrysogine (BGC66) (Table S3).
To establish epistasis, the dataset was compared to the 1120 DEGs that represented MGV1-specific regulated genes (Figure 4). The Venn diagram (Figure 7) differentiated those DEGs that are regulated by MGV1 (Pathways A and B) and those that are regulated by TRI6 (Pathways C and D) (Figure 2). The 24 BGCs regulated by MGV1 alone represent the ‘orphan’ genes of a BGC. In contrast, the 84 BGCs regulated by TRI6 (independently or through MGV1) represent a full complement of genes of a given BGC, including those involved in the production of gramillin, aurofusarin, fusarin, fusaristatin, butanolide, etc. (Table S4).

4. Discussion

MAPKs are essential for the transduction and amplification of environmental signals to the transcriptional regulatory network. Disruption mutants of the MAPK, MGV1, indicate that this gene is involved in the regulation of DON biosynthesis in F. graminearum [21]. The transcription factor, TRI6, regulates the expression of genes of the trichothecene biosynthesis genes (TRI genes) [9]. A connection between MGV1 and TRI6 in this regulation has not previously been investigated. Notably, both genes have also been implicated in other metabolic pathways and BGCs [21,24]. Here, we explored the interaction between MGV1 and TRI6 in the regulation of BGCs in F. graminearum. Comparative transcriptomics was used with various genetic mutants of MGV1 and TRI6 to reveal BGCs that they are independently regulated by the two genes and those BGCs that are dependent on both regulators. Our previous study [24] that combined transcriptome and metabolome suggested that TRI6-regulated BGCs were involved in the production of 15-ADON, fusaoctaxin A (BGC64), and gramillin (BGC02). The regulatory hierarchy governing the biosynthesis of secondary metabolites constructed with the RNA-seq datasets was able to differentiate BGCs controlled by MGV1 and TRI6. The present study was able to discern the regulation of the aforementioned BGCs with BGCs 23 and 64 being co-regulated, while BGC02 was regulated by TRI6, independent of MGV1 (Figure 8, Table S5).
As related in the hierarchy scheme (Figure 2), a combination of A and B, comprising of the MGV1-activated genes, include BGC18 involved in the production of Orcinol and other uncharacterized BGCs. The role of the products synthesized by these uncharacterized BGCs in plant pathogenesis or against other microbes is yet to be determined, but CRISPR technology is facilitating our endeavor in an accelerated fashion. The TRI6 independent pathway (pathway C) includes BGCs that produce gramillin (BGC02), fusaristatin (BGC47), and an uncharacterized BGC16 (Figure 8, Table S5). The BGCs that are epistatic include clusters that produce aurofusarin (BGC13), triacetylfusarinine (BGC21), trichothecene (BGC23), butenolide (BGC49), chrysogine (BGC66), an unknown product synthesized by an NPS (Non-ribosomal peptide synthetase) belonging to BGC72 (Table S5; Pathway D, Figure 2). Since the expression of TRI6 is downregulated by 5–6 fold in the Δmgv1 mutant strain (Figure 5), it is more than likely that the epistasis is through the regulation of TRI6 by MGV1. The study also underpins the importance of TRI6 as a global regulator by influencing other BGCs [8]. The next step would be to understand the mechanism by which TRI6 regulates these BGCs. Our previous study using chromatin immunoprecipitation demonstrated that under PN and NPN conditions, TRI6 is only bound to the promoters of the BGC23 and not to other co-regulated clusters such as BGC47 or BGC02 [24]. This suggested that TRI6 imposes its regulation of other BGCs indirectly through its interactions with other regulators such as GRA2 [24].
It is unclear how the regulation of TRI6 is achieved; post-translational modification through phosphorylation of TRI6, is unlikely as phosphoproteome studies in F. graminearum did not reveal any such modifications and no predictive MAPK phosphorylation domains are present in the translated sequence [28]. TRI6 has been proposed to be regulated partially by TRI10, another regulator in the BGC23, involved in the production of 15-ADON, however, a recent study demonstrated that this may be due to the structural changes of BGC23 and not necessarily through the action of TRI10 [29].
Overall, the analysis of the BGCs reveals that many BGCs are regulated through MGV1, and a subset is finely tuned by the regulator TRI6 embedded within a BGC. It is more than likely that other BGCs in the A/B pathways are similarly fine-tuned by the other transcriptional regulators embedded with that respective BGC. As noted, the three MAPKs in F. graminearum have overlapping functions with respect to the regulation of 15-ADON production (BGC23) [14]. The observations reinforce the idea that MAPKs including MGV1 serve as an articulation point to various pathway-specific regulators. The molecular mechanism, either through differential phosphorylation of transcriptional regulators or through protein interactions, will be the next step to understanding how MGV1 differentiates various signaling pathways.
We also observed that both MGV1 and TRI6 act as positive and negative regulators. For example, TRI6 positively regulated genes in the BGC16, while negatively regulating genes of the BGC that produce gramillin (BGC02), and fusaristatin (BGC47) (Table S5). This may be reflective of the reallocation of resources within the fungus but may also represent antagonistic interactions that may be required to acquire full pathogenicity potential. Efforts are underway to delete the representative core gene in each BGC and characterize their virulence function individually and in combination.
In conclusion, our efforts to study the epistatic relationship between MGV1 and TRI6 were achieved through the manipulation of genes by the CRISPR/Cas9 approach. A traditional approach through sequential deletion/insertion would have been problematic as the Δmgv1 mutant strain with its fitness defect was recalcitrant to protoplasting (weak cell wall) and subsequent regeneration steps. Thus, the CRISPR technique allowed us to achieve the dual manipulation of MGV1 and TRI6 expression using WT protoplasts, sidestepping the fitness defect of the Δmgv1 mutant strain.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/jof9080816/s1, Figure S1: The region of insertion of hygromycin in the coding region of MGV1 in F. graminearum; Figure S2: The region of insertion of GAPDH promoter in the TRI6 promoter in F. graminearum; Figure S3: Confirmation of ∆mgv1/TRI6OX F. graminearum strain by whole genome sequencing; Figure S4: RT-qPCR confirmation of MGV1 and TRI6 expression in the ∆mgv1::TRI6 strains; Figure S5: Growth properties of WT F. graminearum and MGV1Oex strains; Figure S6: 15-ADON accumulation in Δmgv1 and MGV1Oex strains compared with WT; Figure S7: FHB disease assays in wheat cultivars Roblin and Penhold; Figure S8: Constitutive expression of TRI6 affects does not phenocopy MGV1-dependent fitness of F. graminearum; Table S1: Normalized Expression Levels of MGV1 regulated genes; Table S2: DEGs and BGCs regulated by MGV1; Table S3: DEGs and BGCs regulated by TRI6; Table S4: expression levels, DEGs, and BGCs of epistasis; Table S5: Summary of pathway and BGCs; Table S6: primers used in this study.

Author Contributions

Conceptualization, K.S., D.G.-P.F., R.S. and N.A.F.; methodology, K.S., D.G.-P.F., R.S. and N.A.F.; software, K.S., T.W. and C.B.; validation, K.S., D.G.-P.F. and A.E.; formal analysis, K.S., D.G.-P.F., R.S. and N.A.F.; investigation, D.G.-P.F., R.S. and N.A.F.; resources, R.S. and N.A.F.; data curation, K.S., D.G.-P.F., T.W. and C.B.; writing—original draft preparation, K.S. and D.G.-P.F.; writing, K.S., D.G.-P.F., R.S., A.S., D.O. and N.A.F.; review and editing, R.S., N.A.F. and N.T.; visualization, K.S., D.G.-P.F., R.S. and N.A.F.; supervision, R.S. and N.A.F.; project administration, R.S. and N.A.F.; funding acquisition, R.S. and N.A.F. All authors have read and agreed to the published version of the manuscript.

Funding

This project was funded by an internal grant within Agriculture and Agri-Food Canada; grant #514243. A. Eranthodi and D. González-Peña Fundora held Alberta Innovates–Technology Futures graduate student scholarships.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is provided in the Supplementary Tables.

Acknowledgments

We are grateful to Ravinder Goyal, and Daria Ryabova for technical guidance.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ochiai, N.; Tokai, T.; Nishiuchi, T. Involvement of the osmosensor histidine kinase and osmotic stress-activated protein kinases in the regulation of secondary metabolism in Fusarium graminearum. Biochem. Biophys. Res. Commun. 2007, 363, 639–644. [Google Scholar] [CrossRef] [PubMed]
  2. Van Nguyen, T.; Kröger, C.; Bönnighausen, J.; Schäfer, W.; Bormann, J. The ATF/CREB transcription factor Atf1 is essential for full virulence, deoxynivalenol production, and stress tolerance in the cereal pathogen Fusarium graminearum. Mol. Plant-Microbe Interact. 2013, 26, 1379–1394. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Yu, F.; Gu, Q.; Yun, Y.; Yin, Y.; Xu, J.-R.; Shim, W.-B.; Ma, Z. The TOR signaling pathway regulates vegetative development and virulence in Fusarium graminearum. New Phytol. 2014, 203, 219–232. [Google Scholar] [CrossRef] [PubMed]
  4. Yun, Y.; Liu, Z.; Yin, Y.; Jiang, J.; Chen, Y.; Xu, J.-R.; Ma, S. Functional analysis of the Fusarium graminearum phosphatome. New Phytol. 2015, 207, 119–134. [Google Scholar] [CrossRef] [PubMed]
  5. Jiang, C.; Zhang, C.; Wu, C.; Sun, P.; Hou, R.; Liu, H.; Wang, C.; Xu, J.-R. TRI6 and TRI10 play different roles in the regulation of deoxynivalenol (DON) production by cAMP signalling in Fusarium graminearum. Environ. Microbiol. 2016, 18, 3689–3701. [Google Scholar] [CrossRef]
  6. Sieber, C.M.K.; Lee, W.; Wong, P.; Münsterkötter, M.; Mewes, H.-W.; Schmeitzl, C.; Varga, E.; Berthiller, F.; Adam, G.; Güldener, U. The Fusarium graminearum Genome Reveals More Secondary Metabolite Gene Clusters and Hints of Horizontal Gene Transfer. PLoS ONE 2014, 9, e110311. [Google Scholar] [CrossRef] [Green Version]
  7. Seong, K.Y.; Pasquali, M.; Zhou, X.; Song, J.; Hilburn, K.; McCormick, S. Global gene regulation by Fusarium transcription factors Tri6 and Tri10 reveals adaptations for toxin biosynthesis. Mol. Microbiol. 2009, 72, 354–367. [Google Scholar] [CrossRef]
  8. Nasmith, C.G.; Walkowiak, S.; Wang, L.; Leung, W.W.; Gong, Y.; Johnston, A.; Harris, L.J.; Guttman, D.; Subramaniam, R. Tri6 is a global transcription regulator in the phytopathogen Fusarium graminearum. PLoS Pathog. 2011, 7, e1002266. [Google Scholar] [CrossRef]
  9. Brauer, E.K.; Subramaniam, R.; Harris, L.J. Regulation and Dynamics of Gene Expression During the Life Cycle of Fusarium graminearum. Phytopathology 2020, 110, 1368–1374. [Google Scholar] [CrossRef]
  10. Mogg, C.; Bonner, C.; Wang, L.; Schernthaner, J.; Smith, M.; Desveaux, D.; Subramaniam, R. Genomic Identification of the TOR Signaling Pathway as a Target of the Plant Alkaloid Antofine in the Phytopathogen Fusarium graminearum. Mbio 2019, 10, e00792-19. [Google Scholar] [CrossRef] [Green Version]
  11. Turrà, D.; Segorbe, D.; Di Pietro, A. Protein kinases in plant-pathogenic fungi: Conserved regulators of infection. Annu. Rev. Phytopathol. 2014, 52, 267–288. [Google Scholar] [CrossRef]
  12. Zheng, D.; Zheng, D.; Zhang, S.; Zhou, X.; Wang, C.; Xiang, P.; Zheng, Q.; Xu, J.-R. The FgHOG1 pathway regulates hyphal growth, stress responses, and plant infection in Fusarium graminearum. PLoS ONE 2012, 7, e49495. [Google Scholar] [CrossRef]
  13. Hou, Z.; Hou, Z.; Xue, C.; Peng, Y.; Katan, T.; Kistler, H.C.; Xu, J.-R. A mitogen-activated protein kinase gene (MGV1) in Fusarium graminearum is required for female fertility, heterokaryon formation, and plant infection. Mol. Plant-Microbe Interact. 2002, 15, 1119–1127. [Google Scholar] [CrossRef] [Green Version]
  14. Ren, J.; Zhang, Y.; Wang, Y.; Li, C.; Bian, Z.; Zhang, X.; Liu, H.; Xu, J.-R.; Jiang, C. Deletion of all three MAP kinase genes results in severe defects in stress responses and pathogenesis in Fusarium graminearum. Stress Biol. 2022, 2, 2–13. [Google Scholar] [CrossRef]
  15. Yun, Y.; Liu, Z.; Zhang, J.; Shim, W.-B.; Chen, Y.; Ma, Z. The MAPKK FgMkk1 of Fusarium graminearum regulates vegetative differentiation, multiple stress response, and virulence via the cell wall integrity and high-osmolarity glycerol signaling pathways. Environ. Microbiol. 2014, 16, 2023–2037. [Google Scholar] [CrossRef]
  16. Van Nguyen, T.; Schöfer, W.; Bormann, J. The stress-activated protein kinase FgOS-2 is a key regulator in the life cycle of the cereal pathogen Fusarium graminearum. Mol. Plant-Microbe Interact. 2012, 25, 1142–1156. [Google Scholar] [CrossRef] [Green Version]
  17. Zhao, X.; Mehrabi, R.; Xu, J.-R. Mitogen-activated protein kinase pathways and fungal pathogenesis. Fungal Genet. Biol. 2007, 152, 137–152. [Google Scholar] [CrossRef] [Green Version]
  18. Tudzynski, B. Nitrogen regulation of fungal secondary metabolism in fungi. Front. Microbiol. 2014, 5, 656. [Google Scholar] [CrossRef] [Green Version]
  19. Li, B.; Chen, Y.; Tian, S. Function of pH-dependent transcription factor PacC in regulating development, pathogenicity, and mycotoxin biosynthesis of phytopathogenic fungi. FEBS J. 2021, 289, 1723–1730. [Google Scholar] [CrossRef]
  20. Merhej, J.; Richard-Forget, F.; Barreau, C. The pH regulatory factor Pac1 regulates Tri gene expression and trichothecene production in Fusarium graminearum. Fungal Genet. Biol. 2011, 48, 275–284. [Google Scholar] [CrossRef]
  21. Rampitsch, C.; Leung, W.; Blackwell, B.; Subramaniam, R. MAP kinase MGV1: A potential shared control point of butenolide and deoxynivalenol biosynthesis in Fusarium graminearum. Plant Breed. Seed Sci. 2011, 64, 81–88. [Google Scholar] [CrossRef]
  22. Sridhar, P.S.; Trofimova, D.; Subramaniam, R.; González-Peña Fundora, D.; Foroud, N.A.; Allingham, J.A.; Loewen, M.C. Ste2 receptor-mediated chemotropism of Fusarium graminearum contributes to its pathogenicity against wheat. Sci. Rep. 2020, 10, 10770. [Google Scholar] [CrossRef] [PubMed]
  23. Frandsen, R.J.N. A guide to binary vectors and strategies for targeted genome modification in fungi using Agrobacterium tumefaciens-mediated transformation. J. Microbiol. Methods 2011, 87, 247–262. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Shostak, K.; Bonner, C.; Sproule, A.; Thapa, I.; Shield, S.W.J.; Blackwell, B.; Vierula, J.; Overy, D.; Subramaniam, R. Activation of biosynthetic gene clusters by the global transcriptional regulator TRI6 in Fusarium graminearum. Mol. Microbiol. 2020, 114, 664–680. [Google Scholar] [CrossRef]
  25. Miltenburg, M.G.; Bonner, C.; Hepworth, S.; Huang, M.; Rampitsch, C.; Subramaniam, R. Proximity-dependent biotinylation identifies a suite of candidate effector proteins from Fusarium graminearum. Plant J. 2022, 112, 369–382. [Google Scholar] [CrossRef]
  26. Abdallah, Q.; Wenbo, G.; Fortwendel, J. A simple and universal system for gene manipulation in Aspergillus fumigatus: In-vitro- Assembled Cas9- Guide RNA Ribonucleoproteins Coupled with Microhomology Repair Templates. mSphere 2017, 2, e00446-17. [Google Scholar] [CrossRef] [Green Version]
  27. Pfaffl, M.W. A new mathematical model for relative quantification in real-time RT–PCR. Nucleic Acid Res. 2001, 29, e45. [Google Scholar] [CrossRef]
  28. Rampitsch, C.; Subramaniam, R.; Djuric-Ciganovic, S.; Bykova, N.V. The phosphoproteome of Fusarium graminearum at the onset of nitrogen starvation. Proteomics 2009, 10, 124–140. [Google Scholar] [CrossRef]
  29. Liew, M.X.X.; Nakajima, Y.; Maeda, K.; Kitamura, N.; Kimura, M. Regulatory mechanism of trichothecene biosynthesis in Fusarium graminearum. Front. Microbiol. 2023, 14, 1148771. [Google Scholar] [CrossRef]
Figure 1. A scheme to construct F. graminearum Δmgv1/TRI6OX strain. Two repair templates were constructed: (A) first repair template consisted of the selection marker hygromycin that will replace the coding region of MGV1; (B) second repair template consisted of a ~1 kb promoter of GAPDH (pGPD) that will replace the TRI6 promoter; (C) Fusarium protoplasts were incubated with the two repair templates with the ribonucleoprotein (RNP) complex to promote a homologous recombination event.
Figure 1. A scheme to construct F. graminearum Δmgv1/TRI6OX strain. Two repair templates were constructed: (A) first repair template consisted of the selection marker hygromycin that will replace the coding region of MGV1; (B) second repair template consisted of a ~1 kb promoter of GAPDH (pGPD) that will replace the TRI6 promoter; (C) Fusarium protoplasts were incubated with the two repair templates with the ribonucleoprotein (RNP) complex to promote a homologous recombination event.
Jof 09 00816 g001
Figure 2. Schematic representation of hierarchical regulation of BGCs by MGV1 and TRI6. The NPN acts as a stimulus to trigger the MAPK, MGV1, which can activate BGCs (Pathways A and B) or indirectly through the action of TRI6 (Pathway B, Figure 1). Similarly, the stimulus can also regulate TRI6 function through a pathway that is independent of MGV1 (Pathway C, Figure 1) or dependent on MGV (Pathway D, Figure 1).
Figure 2. Schematic representation of hierarchical regulation of BGCs by MGV1 and TRI6. The NPN acts as a stimulus to trigger the MAPK, MGV1, which can activate BGCs (Pathways A and B) or indirectly through the action of TRI6 (Pathway B, Figure 1). Similarly, the stimulus can also regulate TRI6 function through a pathway that is independent of MGV1 (Pathway C, Figure 1) or dependent on MGV (Pathway D, Figure 1).
Jof 09 00816 g002
Figure 3. Differential expression and activation of MAPKs in WT, Δmgv1, and MGV1Oex strains. (A) RT-qPCR confirmation of increased Mgv1 expression in MGV1Oex strains. Relative expression was calculated using the Pfaffl method with EF1α (FGSG_08811) and β-tubulin (FGSG_09530) as reference genes. Bars represent a mean of three biological replicates with error bars representing standard error. Significance was determined by Student’s t-test at p ≤ 0.05; *** p < 0.001 (B) Activation of the Mgv1 enzyme in the MGV1Oex strains. Immunodetection of phosphorylated Mgv1 (pMgv1) and total Mgv1 (pMgv1 and unphosphorylated forms) and the housekeeping protein tubulin from mycelium-grown PDB. One of three replicates was presented. (C) Activation of MAPKs in different strains cultured in NPN. Immunodetection of phosphorylated Mgv1 (pMgv1) and total Mgv1 (pMgv1 and unphosphorylated forms), phosphorylated Gpmk1 (pGpmk1) and total Gpmk1 (pGpmnk1 and unphosphorylated forms), phosphorylated Hog1 (p-Hog1), unphosphorylated Hog1, and the housekeeping protein tubulin from mycelium grown in PNP. Triplicate results are presented in one blot. Note that the Mgv1 protein in Mgv1Oex strains carries a His-tag, which increased the molecular weight of Mgv1 from 46.9 kDa to 47.9 kDa in the overexpression strains.
Figure 3. Differential expression and activation of MAPKs in WT, Δmgv1, and MGV1Oex strains. (A) RT-qPCR confirmation of increased Mgv1 expression in MGV1Oex strains. Relative expression was calculated using the Pfaffl method with EF1α (FGSG_08811) and β-tubulin (FGSG_09530) as reference genes. Bars represent a mean of three biological replicates with error bars representing standard error. Significance was determined by Student’s t-test at p ≤ 0.05; *** p < 0.001 (B) Activation of the Mgv1 enzyme in the MGV1Oex strains. Immunodetection of phosphorylated Mgv1 (pMgv1) and total Mgv1 (pMgv1 and unphosphorylated forms) and the housekeeping protein tubulin from mycelium-grown PDB. One of three replicates was presented. (C) Activation of MAPKs in different strains cultured in NPN. Immunodetection of phosphorylated Mgv1 (pMgv1) and total Mgv1 (pMgv1 and unphosphorylated forms), phosphorylated Gpmk1 (pGpmk1) and total Gpmk1 (pGpmnk1 and unphosphorylated forms), phosphorylated Hog1 (p-Hog1), unphosphorylated Hog1, and the housekeeping protein tubulin from mycelium grown in PNP. Triplicate results are presented in one blot. Note that the Mgv1 protein in Mgv1Oex strains carries a His-tag, which increased the molecular weight of Mgv1 from 46.9 kDa to 47.9 kDa in the overexpression strains.
Jof 09 00816 g003
Figure 4. Venn diagram showing differentially expressed genes (DEGs) to identify MGV1-specific regulated genes. Venn diagram showing differentially expressed genes (DEGs) to identify MGV1-specific regulated genes. (A) DEGs from wild-type and Δmgv1 strains grown in NPN conditions (1495) were compared to DEGs of the Δmgv1 and the two in-locus MGV1 overexpression strains (Oex1; 1842; and Oex6; 2101) strains grown in NPN conditions. (B) An overlap (1120 DEGs) that are common to both overexpression strains, representing MGV1 regulated genes.
Figure 4. Venn diagram showing differentially expressed genes (DEGs) to identify MGV1-specific regulated genes. Venn diagram showing differentially expressed genes (DEGs) to identify MGV1-specific regulated genes. (A) DEGs from wild-type and Δmgv1 strains grown in NPN conditions (1495) were compared to DEGs of the Δmgv1 and the two in-locus MGV1 overexpression strains (Oex1; 1842; and Oex6; 2101) strains grown in NPN conditions. (B) An overlap (1120 DEGs) that are common to both overexpression strains, representing MGV1 regulated genes.
Jof 09 00816 g004
Figure 5. MGV1 affects the expression of BGC genes. The wild-type (WT), ∆mgv1, MGV1Oex1, and MGV1Oex6 strains were grown in NPN conditions for 24 h, total RNA was isolated, and RT-qPCR was performed with six genes that are differentially expressed in the RNA-seq analysis. TRI5, TRI6, NPS5, and NPS9 are positively regulated by MGV1 (AD), PKS12 and GIP1 are negatively regulated (E,F) Relative expression was calculated using the Pfaffl method with EF1α (FGSG_08811) and β-tubulin (FGSG_09530) as reference genes. Bars represent a mean of three biological replicates with error bars representing standard deviation. Significance (*) was determined with Student’s t-test at p ≤ 0.05.
Figure 5. MGV1 affects the expression of BGC genes. The wild-type (WT), ∆mgv1, MGV1Oex1, and MGV1Oex6 strains were grown in NPN conditions for 24 h, total RNA was isolated, and RT-qPCR was performed with six genes that are differentially expressed in the RNA-seq analysis. TRI5, TRI6, NPS5, and NPS9 are positively regulated by MGV1 (AD), PKS12 and GIP1 are negatively regulated (E,F) Relative expression was calculated using the Pfaffl method with EF1α (FGSG_08811) and β-tubulin (FGSG_09530) as reference genes. Bars represent a mean of three biological replicates with error bars representing standard deviation. Significance (*) was determined with Student’s t-test at p ≤ 0.05.
Jof 09 00816 g005
Figure 6. A Venn diagram to distinguish DEGs and BGCs regulated by TRI6 and MGV1. (A) datasets from MGV-specific regulated genes (Δmgv1/WT; 1120 DEGs) and TRI6-specific regulated genes (Δtri6/WT; 687 DEGs) were compared and identified both DEGs and BGCs regulated either independently or dependent of MGV and TRI6; (B) for example, 228 DEGs are co-regulated by both MGV1 and TRI6 also identifies the pathways and corresponds to pathway D. The details on the DEGs and BGCs are presented in Table S2.
Figure 6. A Venn diagram to distinguish DEGs and BGCs regulated by TRI6 and MGV1. (A) datasets from MGV-specific regulated genes (Δmgv1/WT; 1120 DEGs) and TRI6-specific regulated genes (Δtri6/WT; 687 DEGs) were compared and identified both DEGs and BGCs regulated either independently or dependent of MGV and TRI6; (B) for example, 228 DEGs are co-regulated by both MGV1 and TRI6 also identifies the pathways and corresponds to pathway D. The details on the DEGs and BGCs are presented in Table S2.
Jof 09 00816 g006
Figure 7. Venn diagram to identify epistatic regulation of BGCs. Epistatic regulation of BGCs between MGV1 and TRI6 was established by constitutive overexpression of TRI6 in the Δmgv1 mutant strain (Δmgv1::TRI6OX/WT). DEGs and BGCs were compared between the Δmgv1::TRI6OX/WT and the Δmgv1/WT strains. A, B, C, and D refers to the signaling pathways regulated by MGV1 and TRI6 as outlined in Figure 1.
Figure 7. Venn diagram to identify epistatic regulation of BGCs. Epistatic regulation of BGCs between MGV1 and TRI6 was established by constitutive overexpression of TRI6 in the Δmgv1 mutant strain (Δmgv1::TRI6OX/WT). DEGs and BGCs were compared between the Δmgv1::TRI6OX/WT and the Δmgv1/WT strains. A, B, C, and D refers to the signaling pathways regulated by MGV1 and TRI6 as outlined in Figure 1.
Jof 09 00816 g007
Figure 8. A summary of BGCs regulated by the regulatory hierarchy of MGV1 and TRI6 in F. graminearum. The BGCs activated through MGV1 signaling (A/B) by TRI6 (B) that is either independent of MGV1 (C) or dependent on MGV1 (D) are indicated.
Figure 8. A summary of BGCs regulated by the regulatory hierarchy of MGV1 and TRI6 in F. graminearum. The BGCs activated through MGV1 signaling (A/B) by TRI6 (B) that is either independent of MGV1 (C) or dependent on MGV1 (D) are indicated.
Jof 09 00816 g008
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Shostak, K.; González-Peña Fundora, D.; Blackman, C.; Witte, T.; Sproule, A.; Overy, D.; Eranthodi, A.; Thakor, N.; Foroud, N.A.; Subramaniam, R. Epistatic Relationship between MGV1 and TRI6 in the Regulation of Biosynthetic Gene Clusters in Fusarium graminearum. J. Fungi 2023, 9, 816. https://doi.org/10.3390/jof9080816

AMA Style

Shostak K, González-Peña Fundora D, Blackman C, Witte T, Sproule A, Overy D, Eranthodi A, Thakor N, Foroud NA, Subramaniam R. Epistatic Relationship between MGV1 and TRI6 in the Regulation of Biosynthetic Gene Clusters in Fusarium graminearum. Journal of Fungi. 2023; 9(8):816. https://doi.org/10.3390/jof9080816

Chicago/Turabian Style

Shostak, Kristina, Dianevys González-Peña Fundora, Christopher Blackman, Tom Witte, Amanda Sproule, David Overy, Anas Eranthodi, Nehal Thakor, Nora A. Foroud, and Rajagopal Subramaniam. 2023. "Epistatic Relationship between MGV1 and TRI6 in the Regulation of Biosynthetic Gene Clusters in Fusarium graminearum" Journal of Fungi 9, no. 8: 816. https://doi.org/10.3390/jof9080816

APA Style

Shostak, K., González-Peña Fundora, D., Blackman, C., Witte, T., Sproule, A., Overy, D., Eranthodi, A., Thakor, N., Foroud, N. A., & Subramaniam, R. (2023). Epistatic Relationship between MGV1 and TRI6 in the Regulation of Biosynthetic Gene Clusters in Fusarium graminearum. Journal of Fungi, 9(8), 816. https://doi.org/10.3390/jof9080816

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