The Dis1/Stu2/XMAP215 Family Gene FgStu2 Is Involved in Vegetative Growth, Morphology, Sexual and Asexual Reproduction, Pathogenicity and DON Production of Fusarium graminearum

The conserved Dis1/Stu2/XMAP215 microtubule association proteins (MAPs) family plays an important role in microtubule dynamics, nucleation, and kinetochore-microtubule attachments. However, function of Dis1/Stu2/XMAP215 homolog in plant pathogenic fungi has not been determined. Here, we identified and investigated the Dis1/Stu2/XMAP215 homolog (FGSG_10528) in Fusarium graminearum (FgStu2p). Co-localization experiment and co-immunoprecipitation (Co-IP) assay demonstrated that FgStu2p is a microtubule associated protein. Besides, FgStu2 could also interact with Fgγ-tubulin and presumed FgNdc80, which suggested that the FgStu2 gene might associate with microtubule nucleation and kinetochore-microtubule attachments like Dis1/Stu2/XMAP215 homologs in other species. Moreover, the FgStu2 promoter replacement mutants (FgStu2-Si mutants) produced twisted hyphae and decreased growth rate. Microscope examination further showed that the microtubule polymerization was reduced in FgStu2-Si mutants, which could account for the aberrant morphology. Although the microtubule polymerization was affected in FgStu2-Si mutants, the FgStu2-Si mutants didn’t show highly increased sensitivity to anti-microtubule fungicide carbendazim (methyl benzimidazol-2-ylcarbamate [MBC]). In addition, the FgStu2-Si mutants exhibited curved conidia, decreased number of conidial production, blocked ability of perithecia production, decreased pathogenicity and deoxynivalenol (DON) production. Taken together, these results indicate that the FgStu2 gene plays a crucial role in vegetative growth, morphology, sexual reproduction, asexual reproduction, virulence and deoxynivalenol (DON) production of F. graminearum, which brings new insights into the functions of Dis1/Stu2/XMAP215 homolog in plant pathogenic fungi.

Introduction

Microtubules are hollow cylindrical polymers assembled from α/β-tubulin heterodimers that have a specific behavior: microtubules switch between growing or shrinking states, a property known as dynamic instability, which is vital to many intracellular activities such as mitosis, material transportation and cell morphology (Nogales, 2001; Tischfield et al., 2010; Aher and Akhmanova, 2018). Individual microtubule has two different tips with distinct conformations, the minus end and the plus end. The minus end is the nucleation site of microtubule. While the plus end is the growing end which could recruit or disassociate α/β-tubulin heterodimers (Schuyler and Pellman, 2001).

Microtubule association proteins (MAPs) are a great kind of proteins that could attach to microtubule or tubulin and regulate microtubule dynamics (Walczak, 2000). One of the most important MAPs families is the Dis1/Stu2/XMAP215 family which was first identified in Xenopus laevis (Gard and Kirschner, 1987). In the present, the Dis1/Stu2/XMAP215 family has been widely studied in X. laevis (XMAP215), Human (ch-TOG, human XMAP215 homolog), budding yeast (Stu2, budding yeast XMAP215 homolog), and fission yeast (Dis1, fission yeast XMAP215 homolog) (Al-Bassam et al., 2006; Al-Bassam and Chang, 2011). The Dis1/Stu2/XMAP215 family was famous as they have 2-5 conserved tumor-overexpressed gene (TOG) domains at amino-terminus and a microtubule binding domain at carboxy-terminus (Thawani et al., 2018). TOG domain could bind to α/β-tubulin heterodimers and is essential for microtubule polymerization activity of Dis/Stu2/XMAP215 family (Ayaz et al., 2012; Geyer et al., 2018). Previous studies showed that knockdown of these proteins often leads to small or abnormally organized spindles and short astral microtubules (Al-Bassam and Chang, 2011). Moreover, the Dis1/Stu2/XMAP215 family could also act together with other MAPs such as end binding 1 (EB1) to regulate microtubule dynamics (Zanic et al., 2013). At kinetochores, the Dis1/stu2/XMAP215 family associates with the Ndc80 kinetochore complex and regulates the kinetochore-microtubule attachment which is vital for cytokinesis (Hsu and Toda, 2011; Miller et al., 2016). Recent studies showed that the Dis1/stu2/XMAP215 family is a microtubule nucleation factor that functions synergistically with the γ-tubulin ring complex (Thawani et al., 2018). In addition, in budding yeast, Stu2p promote oligomerization of the γ-tubulin ring complex (γ-TuRC) and cytoplasmic microtubule nucleation via interaction with Spc72 (Chen et al., 1998). Taken together, Dis/Stu2/XMAP215 family is essential for microtubule polymerization and function of kinetochore and disrupt these gene is fatal for the cell. Although Dis/Stu2/XMAP215 family has been well studied in human, yeast and X. laevis, its role in plant pathogenic fungi has not been determined. Therefore, in this study, we investigated the functions of Dis/Stu2/XMAP215 family homolog in Fusarium graminearum (FgStu2).

Fusarium graminearum (teleomorph: Gibberella zeae), a universal plant pathogenic fungi which causes the Fusarium head blight (FHB), a destructive disease in cereal crops (Bai and Shaner, 2004). In addition to give rise to serious yield loss, FHB can also contaminate grain by and mycotoxin such as deoxynivalenol (DON) and zearalenone (ZEA), which posed severe threaten to human and livestock health (Forsyth et al., 1977; Chen et al., 2019). DON is the most widely detected mycotoxin in cereal corps with a 56% average incidence rate worldwide compared with 14% of zearalenone and other mycotoxin (Lee and Ryu, 2017). Due to the highest occurrence and toxicity, DON has gained more attention and has been widely studied (Lee and Ryu, 2017; Chen et al., 2019). Although many genes of F. graminearum have been demonstrated to regulate the DON production, the role of microtubule-associated proteins such as the Dis/Stu2/XMAP215 family in DON production has not been widely characterized (Chen et al., 2019). Therefore, we are interested in the functional role of FgStu2 gene in DON production.

In the present, control of FHB is mainly depended on chemical fungicides due to lack of fully resistant cultivars (Blandino et al., 2012). One of the most widely used fungicides for control FHB is carbendazim (methyl benzimidazol-2-ylcarbamate [MBC]), which binds to tubulin heterodimers and disrupt microtubule dynamics (Davidse, 1986; Chen et al., 2009). Previous studies have shown that F. graminearum contains two β-tubulin genes: β₁- and β₂-tubulin, and mutations such as F167Y, E198K, and F200Y in β₂-tubulin gene confer resistance to MBC (Chen et al., 2009; Qiu et al., 2012). In this study, we also investigated the relationship between FgStu2 and Fgβ1- and Fgβ2-tubulin and the role of FgStu2 in regulation of MBC sensitivity. By the means of promoter substitution, Co-immunoprecipitation (Co-IP), EGFP (Enhanced green fluorescent protein)-tagging, and RFP (Red fluorescent protein)-tagging, we revealed the functional importance of FgStu2 gene in F. graminearum. Our results will provide information about the functional role of FgStu2 gene and its relationship with tubulin.

Materials and Methods

Sequence Alignment and Analysis

We used the amino acid sequences of Saccharomyces cerevisiae Stu2 (Gene bank ID: 856736) and S. pombe Dis1 (Gene bank ID: 2539505) to find the Dis1/Stu2/XMAP215 homologs in F. graminearum (FgStu2 gene bank ID: 23557431). The protein blast was carried out by online blast services of the National Center for Biotechnology information (NCBI¹) and the Ensembl Fungi². The presumed FgNdc80 gene (FGSG_09262, gene bank ID: 23556224) was also identified by online blast services of the NCBI. The phylogenetic analysis (phylogenetic tree) of Dis1/Stu2/XMAP215 homologs of Aspergillus nidulans (Protein ID: XP_663125.1), F. graminearum (Protein ID: XP_011319522.1), Neurospora crassa (Protein ID: XP_956946.3), S. cerevisiae (Protein ID: NP_013146.1) and S. pombe (Protein ID: NP_587785.1) was finished by the Molecular Evolutionary Genetics Analysis software (MEGA, version: 7.0) (Supplementary Figure 1). The domain prediction was conducted by the online service of Simple Modular Architecture Research Tool (SMART³) (Letunic et al., 2014). Multiple sequences alignments of tumor-overexpressed gene (TOG) domains and terminal areas in Dis/Stu2/XMAP215 homologs of A. nidulans (Protein ID: XP_663125.1), F. graminearum (Protein ID: XP_011319522.1), N. crassa (Protein ID: XP_956946.3), S. cerevisiae (Protein ID: NP_013146.1), and S. pombe (Protein ID: NP_587785.1) were performed using the online service⁴ (Figures 1B,C and Supplementary Figure 2). The amino acids sequences alignment of presumed FgNdc80 (protein ID: XP_011328506.1) and budding yeast Ndc80 (protein ID: NP_012122.3) was performed using the same online service (Supplementary Figure 3). The alignment images were generated by the online service of ESPript 3.0⁵ (Robert and Gouet, 2014).

FIGURE 1

Figure 1. Domain analysis of Stu2 homologs and sequences alignment. (A) Domain analysis of Dis1/Stu2/XMAP215 homologs in Aaspergillus nidulans (Protein ID: XP_663125.1), Fusarium graminearum (Protein ID: XP_011319522.1), Neurospora crassa (Protein ID: XP_956946.3), Saccharomyces cerevisiae (Protein ID: NP_013146.1) and Schizosaccharomyces pombe (Protein ID: NP_587785.1). The Domain prediction was conducted by the online service of the Simple Modular Architecture Research Tool (SMART, https://smart.embl.de). (B) Sequences alignment of TOG1 domain. Multiple sequences alignment was performed using the online service (www.genome.jp/tools-bin/clustalw). The alignment images were generated by the online service of ESPript 3.0 (espript.ibcp.fr/ESPript/cgi-bin/ ESPript.cgi). (C) Sequences alignment of TOG2 domain. The red background indicates the same amino acid residues.

Strains, Cultural Condition, and Fungicides

All strains used in this study are listed in Table 1. The PH-1 strain is the standard strain of F. graminearum which was completely sequenced in Cuomo et al. (2007). In the recent years, PH-1 strain has become a model strain in plant pathology and fungal genetics (Liu et al., 2017; Chen et al., 2019; Jiang et al., 2020). More importantly, the genetic transformation and fluorescence labeling technologies of PH-1 are now perfect and convenient. Taken together, we would like to use PH-1 strain to study the functional role of Dis1/Stu2/XMAP215 homolog in plant pathogenic fungi and the results would be more representative than the other wild-type strain. Conidia and mycelia were produced in mung bean liquid (MBL) medium and yeast extract peptone dextrose (YEPD) medium at 25°C respectively. For colonial morphology, strains were grown on complete medium (CM) and minimal medium (MM) plants at 25°C (Gu et al., 2015). For sensitivity assays, all strains were grown on potato dextrose agar (PDA) containing carbendazim (MBC) at 25°C. DON mycotoxin was produced in trichothecene biosynthesis induction (TBI) medium (Gardiner et al., 2009). MBC were dissolved in 0.1 M HCl at 10 mg/ml as stock solutions (Chen et al., 2009).

TABLE 1

Table 1. Strains used in this study.

Construction of EGFP-, mRFP-, 3 × FLAG-Tagging Mutants

The Fgβ₁-tubulin-EGFP, Fgβ₂-tubulin-EGFP single labeling strains (Table 1) were stored in our lab. The FgStu2-3 × FLAG strain was constructed using homologous recombination strategy (Supplementary Figure 4), first a 1.0 kb fragment of FgStu2 3′-terminus coding sequence was amplified from genomic DNA of PH-1 by primer containing 1 × FLAG-tag sequence. This fragment was used as a template for another two rounds of PCR to add 2 × FLAG-tag. Second, a 1.2 kb geneticin (G418) fragment was amplified from pNEO. Third, a 1.0 kb FgStu2 downstream fragment was amplified from genomic DNA of PH-1. Finally, the three fragments were fused by DJ-PCR and used as a template for amplify the FgStu2-3 × FLAG vector. Transformates were identified by PCR, sequencing and western blot (Figure 2 and Supplementary Figure 4). For constructing the FgStu2-Si-β1EGFP and FgStu2-Si-β2EGFP strains, we cloned the Fgβ1-tubulin and Fgβ2-tubulin coding sequences (CDS) into PYF11 plasmid and then transformed these plasmids into protoplasts of FgStu2-Si mutant to generate the FgStu2-Si-β1EGFP and FgStu2-Si-β2EGFP strains (Bruno et al., 2004). For constructing the Fgβ2EGFP-Stu2mRFP strain, we cloned RFP-Stu2 DNA fragment into PYF11 plasmid and then transformed this plasmid into protoplasts of Fgβ₂-tubulin-EGFP strain to generate the Fgβ2EGFP-Stu2mRFP strain. The Fgγ-tubulin-EGFP and FgNdc80-EGFP strains were constructed using homologous recombination strategy. First, a 1.0 kb 3′-terminal CDS of Fgγ-tubulin (Gene bank ID: 23556916) and presumed FgNdc80 gene (Gene bank ID: 23556224) was amplify from the genomic DNA of wild-type strain PH-1; a 0.7 kb EGFP CDS and 1.7 kb hygromycin resistance gene (HPH) fragment were amplified from PYF11 plasmid and pNEO plasmid respectively, and then fused the two fragments by double-joint PCR to generate a 2.4 kb GFP-HPH fragment; a 1.0 kb downstream fragment of Fgγ-tubulin and presumed FgNdc80 gene was amplified from the genomic DNA of PH-1. Finally, fused the three fragments by double-joint PCR to generate the Fgγ-tubulin-EGFP and FgNdc80-EGFP vectors and then transferred the vectors into protoplasts of PH-1. The transformates were identified by PCR, sequencing, and western blot (Figure 3 and Supplementary Figure 5). Protoplasts transformations were performed as previously described (Zheng et al., 2014).

FIGURE 2

Figure 2. Microtubule-associated properties of FgStu2p (A) Co-localization of RFP-FgStu2 and Fgβ-tubulin-EGFP. Fresh spores (10⁵) of all double-labeling mutants were first incubated into 20 ml of yeast extract-peptone-dextrose medium (YEPD) on an orbital shaker (175 rpm) at 25°C to germinate for 10 h and then spotted 5 μl of culture onto the slides for micro examination. Bar = 25 μm. (B) Co-immunoprecipitation (Co-IP) analysis of the interaction of FgStu2p and Fgβ1/β2-tubulin. The INPUT (total protein) indicates that the double labeling strains (Fgβ1EGFP-Stu2-3 × FLAG and Fgβ2EGFP-Stu2-3 × FLAG) and single labeling strains (FgStu2-3 × FLAG, Fgβ1-tubulin-EGFP, Fgβ2-tubulin-EGFP) were successfully labeled with FLAG and EGFP tags respectively; the protein samples were also detected with anti-actin antibody as a reference; the IP (Immunoprecipitation, polyclonal anti-EGFP rabbit antibody) indicates the Fgβ1-tubulin-EGFP and Fgβ2-tubulin-EGFP were precipitated by Protein-G magnetic beads which had been incubated with polyclonal anti-EGFP rabbit antibody; the IB (Immune blotting, monoclonal anti-FLAG mouse antibody) indicates that the FgStu2-3 × FLAG was pulled down by the Fgβ1-tubulin-EGFP or Fgβ2-tubulin-EGFP.

FIGURE 3

Figure 3. Co-immunoprecipitation (Co-IP) analysis of interactions of FgStu2p-Fgγ-tubulin and FgStu2p-FgNdc80p. The INPUT (total protein) indicates that the double labeling strains (FgγEGFP-Stu2-3 × FLAG and FgNdc80-EGFP-Stu2-3 × FLAG) and single labeling strains (FgStu2-3 × FLAG, Fgγ-tubulin-EGFP, FgNdc80-EGFP) were successfully labeled with FLAG and EGFP tags respectively; the IP (Immunoprecipitation, polyclonal anti-EGFP rabbit antibody) indicates the Fgγ-tubulin-EGFP, FgNdc80-EGFP were precipitated by Protein-G magnetic beads which had been incubated with polyclonal anti-EGFP antibody; the IB (Immune blotting, monoclonal anti-FLAG mouse antibody) indicates that the FgStu2-3 × FLAG was pulled down by the Fgγ-tubulin-EGFP, FgNdc80-EGFP.

Co-immunoprecipitation (Co-IP) Assays

Strains containing EGFP or 3 × FLAG tag were inocubated into MBL medium for 3 days. After that, fresh spores were inoculated into 100 ml YEPD medium for 24 h at 25°C with agitation (175 rpm) and then fresh mycelia were collected and ready for use. For total protein extraction, fresh dried pressed mycelia were quickly frozen in liquid nitrogen and ground into powder. For each strain, 0.15 g frozen mycelial powder was suspended in 1.5 mL of extraction buffer (20 mM Tris, 200 mM NaCl, 0.1% Triton X-100, 1 × protease inhibitor cocktail [20124ES10, YiSheng Biotechnology Inc., Shanghai, China], PH = 7.5). After incubation on ice for 30 min, cell debris was spun down at 12000g for 20 min at 4°C, and the clear supernatant was transferred to a new tube. The Co-IP assay was conducted follow the standard immunoprecipitation protocol of protein G magnetic beads (SureBeads^TM, Bio-Rad Laboratories, lnc).

Western Blotting Assays

Protein samples was added 5 × Loading buffer and boiled for 5 min before loaded into Sodium dodecyl sulfate-polyacrylamide gel. Protein separated by Sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) was transferred to polyvinylidene fluoride membrane (PVDF) with a Bio-Rad electroblotting device. After that, the PVDF membrane was incubated with 5% defatted milk (5 g non-fat dry milk diluted into 100 ml TBST buffer [10 mM Tris, 100 mM NaCl, 0.1% Tween 20, PH = 7.4]) and blocking for 1 h at 25°C. After Blocking, antibody (polyclonal Anti-EGFP Rabbit antibody, monoclonal Anti-FLAG mouse antibody and monoclonal Anti-Actin antibody were purchased form Zhengneng Biotechnology, Inc., Chengdu, China) was added into the defatted milk. After incubation at 25°C for 2 h, the PVDF membrane was washed with TBST buffer for 5 min with 3 buffer changes and then incubated for 1 h using horse radish peroxidase (HRP)-conjugated secondary antibody diluted into 5% defatted milk. Finally, the PVDF membrane was washed three times by TBST and ready for blotting.

Promoter Replacement Strategy

To study the functional role of FgStu2 gene, we substituted original FgStu2 promoter with an inducible promoter Pzear (Promoter of the FGSG_04581.3 gene) to generate FgStu2 silent mutant (FgStu-Si) as previously described (Lee et al., 2010). The schematic of promoter replacement was shown in Figure 4. First, a 1.2 kb geneticin resistance gene (G418) and a 0.8 kb Pzear fragment were amplified from plasmid pNEO and genomic DNA of PH-1 respectively and fused the two fragments to obtain the G418-Pzear fragment by double joint polymerase chain reaction (DJ-PCR) as previously described (Yu et al., 2004; Ren et al., 2014). Second, a 1.0 kb upstream region of FgStu2 original promoter (Pstu2) and a 1.0 kb FgStu2 coding sequence were amplified from genomic DNA of PH-1. Finally, the three fragments (upstream region of Pstu2, G418-Pzear and 1.0 kb FgStu2 coding sequence) were fused by DJ-PCR to obtain the Pstu2 substitution vector which was used for protoplast transformation. To induce Pzear replacement, 30 μM zearalenone (ZEA, Sigma Aldrich, St. Louis, MO, United States) was added to the medium during the regeneration (Lee et al., 2010).

FIGURE 4

Figure 4. Promoter replacement strategy and identification of promoter replacement mutant of FgStu2 gene. (A) Schematic of promoter replacement of FgStu2 gene. (B) Polymerase chain reaction (PCR) identification of FgStu2 promoter replacement mutants. (a) a 1.9 kb fragment was amplified by primers P9/P10 indicates G418-Pzear fragment integrated at the 3′-terminal of up-stream non-coding sequence of FgStu2 gene. (b) a 0.9 kb fragment could not amplified by primers P11/P12 indicates original promoter of FgStu2 has been substituted with G418-Pzear. (C) a 1.8 kb fragment was amplified by primers P13/P14 indicates that G418-Pzear integrated at the 5′-terminal of FgStu2 coding sequence (CDS). (C) Southern blot identification of FgStu2 promoter replacement mutants. (D) Real-time quantitative PCR identification of FgStu2 promoter replacement mutants.

RNA Extraction and Quantitative Real-Time PCR

To test the expression level of FgStu2 gene, the mycelium was harvested after cultivation for 36 h in YEPD medium. This time is sufficient for gaining enough mycelium and distinguishing the expression level of FgStu2 gene in wild-type strain and FgStu2-Si mutants. To test the expression level of FgTri5 and FgTri6 genes, the mycelium was harvested after three-days cultivation in TBI medium at 28°C in the dark (Gardiner et al., 2009). The reason why we chose three days is that the level of DON increases obviously between 48 and 72 h and then increase slowly after 72 h and reaches the stationary phase (Tang et al., 2018). Total RNA extraction was follow the protocol of Total RNA Extraction Kit (Tiangen, China, DP419). First cDNA was synthesized by HiScript^® II Reverse Transcriptase (Vazyme, Nanjing, China, R223-01). The ChamQTM SYBR^® qPCR Master Mix (Vazyme, Nanjing, China, Q311-02/03) was used for Quantitative real-time PCR (qRT-PCR), which was conducted at CFX Connect Real-Time System (Bio-Rad) with the following procedure: 95°C 30 s; following 40 cycles of 95°C 10 s and 60°C 30 s; and a melt curve step of 95°C for 15 s and 60°C for 60 s, followed by 71 cycles of gradual increase in temperature from 60 to 95°C for 13 s/cycle. The primers used for qRT-PCR analysis are listed in Supplementary Table 1. The endogenous housekeeping gene Fgactin was used for normalization.

Morphology Observation

For photograph the colonial morphology, 5-mm mycelial plugs from the margin of a 3-day-old colony (on PDA plates) of wild-type strain PH-1 and FgStu2 silent mutant were transferred on PDA, MM, and CM plates with or without 30 mM ZEA and keep it at 25°C for three days. After that, the colonial morphology was photographed. For hyphal morphology, wild-type strain PH-1 and FgStu2 silent mutants were inoculated on PDA plates with or without 30 μM ZEA for 24 h. After that, the margin of the small colony was photographed using an Olympus IX-71 inverted microscope.

Fluorescence Observation of Fgβ1-Tubulin, Fgβ2-Tubulin, and FgStu2-mRFP

To observe the localizations of FgStu2 and Fgβ1-tubulin, Fgβ2-tubulin, the fresh spores of Fgβ2EGFP-Stu2mRFP were inoculated into 20 ml YEPD medium with or without 30 μM ZEA to germinate at 25°C for 10-12 h and then pipette 10 μl medium onto the slide for micro examination. This experiment was conducted by a LEICA TCS SP8 confocal laser-scanning microscope (LEICA).

Growth Rate and MBC Sensitivity Tests

The growth rate was tested as previously described (Zhu et al., 2018). In brief, 5 mm mycelial plugs taken from the margin of a 3-day-old colony of wild-type strain PH1 and FgStu2 silent mutants were transferred into PDA plates with or without 30 μM ZEA and incubated at 25°C. After three days cultivation, the colony diameters of wild-type strain and FgStu2-Si mutants were significantly different. Then the diameters of each strain were measured to calculate the mean growth rates. This experiment was conducted three independent times.

Radial growth was used to test the sensitivity of the wild-type strain PH-1 and FgStu2 mutants to MBC as previously described (Zhu et al., 2018). Colony diameter was measured after 3 days for PH1; after 7 days for FgStu2-Si mutants. The MBC was added to PDA plates at 0.2, 0.4, 0.6, 0.8, 1.0, and 1.2 μg/ml for PH-1 and FgStu2-Si mutant; at 5, 10, 15, 20, and 25 μg/ml for Fgβ2F167Y and F167Y-FgStu2-Si mutant.

Assessment of Asexual Reproduction

In order to assess the ability of asexual reproduction of FgStu2-Si mutants, we observed the conidial morphology and tested the number of conidial production of FgStu2-Si mutant as previously described (Zhu et al., 2018). For observe the septa and nuclei, fresh spores were stained with calcofluor white (CFW) and 4′,6-diamidino-2-phenylindole (DAPI) and examined with a LEICA TCS SP8 confocal laser-scanning microscope (LEICA). The numbers of conidia produced by the wild-type strain PH-1 and FgStu2-Si mutants were calculated using hemocytometer as previously described (Zheng et al., 2014). This experiment was conducted three independent times.

Assessment of Sexual Reproduction

Sexual reproduction test was performed as previous described (Bui et al., 2016). Briefly, inoculate fresh mycelia of each stain to carrot medium plates with or without 30 μM ZEA and keep it at 25°C incubator until the plates are full of mycelia. Then clean out the mycelia and add 2.5% Tween 60 solution to the plates. Keep the plates into 25°C incubator with near-UV light for about one or two weeks until the black perithecia are appeared and then photograph the plates by stereomicroscope. This experiment was conducted three independent times.

Pathogenicity Assay

In order to compare the pathogenicity of the wild-type strain PH-1 and FgStu2-Si mutants we preformed wheat coleoptile infection assay as previous described (Zhang et al., 2012). In brief, Huaimai 33 wheat seeds were placed in culture dishes and placed it into a 25°C chamber with light. After 3 days, remove the 2-3 mm tip of all coleoptiles and the wound was inoculated with 2.5 μl (1 × 10⁶ spores/ml) fresh conidia, which was collected from MBL medium with or without 30 μM ZEA. Coleoptiles inoculated with MBL medium with or without 30 μM ZEA were used as control. After 7 to 10 days cultivation, the infected coleoptiles were photographed and the length of lesion was measured. At least 20 wheat seedlings were measured for each strain. This experiment was conducted three independent times.

Assessment of DON Reproduction

Deoxynivalenol production of FgStu2-Si mutants were determined by enzyme linked immunosorbent assay (ELISA) followed the instruction of DON detection kit (Weisai Technology, Jiangsu, China). Briefly, two fresh 5 mm mycelial plugs taken from the margin of a 3-day-old colony of wild-type strain PH-1 and FgStu2-Si mutants were transferred into 30 ml TBI medium with or without 30 μM ZEA and cultured at 28°C, 175 rpm in the dark (3 repeats were set for each strain). After 7 days cultivation, the medium and mycelium was collected. The medium was used for DON determination, the mycelium was dried in the dryer for 48h. After that the dry weight of the mycelium was used for calculating the DON production (μg DON/g mycelium). The DON production of wild-type strain PH-1 was used as reference to generate relative DON production of FgStu2-Si mutants. The experiment was repeated three independent times.

For visualize the localization FgTri1p (also known as “toxisomes”) in FgStu2-Si mutant and wild-type strain PH-1, we cloned FgTri1 gene into PDL2 plasmid as previous described (Liu et al., 2017). The recombinant plasmid was identified by PCR, and sequencing. Then we transformed the recombinant plasmid into protoplasm of FgStu2-Si mutant and wild-type strain PH-1 to generate FgStu2-Si-Tri1-EGFP and PH-1-Tri1-EGFP strains. The transformants were identified by PCR, sequencing and fluorescence microscope. The identified FgStu2-Si-Tri1-EGFP and PH-1-Tri1-EGFP strains were inoculated into TBI medium and incubated for 2 days at 28°C, 175 rpm in the dark, the fresh mycelium were used for fluorescence observation. This experiment was conducted by a LEICA TCS SP8 confocal laser-scanning microscope (LEICA). This experiment was repeated three independent times.

Results

Identification and Amino Acid Sequence Analysis of FgStu2 Gene

The protein blast results of NCBI and Ensembl Fungi website suggested that the FGSG_10528 (gene bank ID: 23557431) is the Dis1/Stu2/XMAP215 family homolog gene in F. graminearum. The FgStu2 gene is predicted to encode an 888 amino acids protein (FgStu2p) which is similar with the lengths of S. cerevisiae Stu2 and S. pombe Dis1 (Figure 1A). Sequence analysis showed that FgStu2p contains two TOG domains at the amino terminal, TOG1 domain (Position 1 to 233) and TOG2 domain (position 273 to 505) (Figure 1A). The FgStu2-TOG1 domain shares only 32% identity with that of S. cerevisiae Stu2 (ScStu2) and S. pombe Dis1 (SpDis1). The FgStu2-TOG2 domain shares only 26 and 33% identity with that of S. cerevisiae Stu2 and S. pombe Dis1. However, the TOG1 domain of FgStu2p share 63 and 71% identify with that of A. nidulans alpA and Dis1/Stu2/XMAP215 homolog of N. crassa respectively (Figure 1B). The TOG2 domain of FgStu2 share 68 and 74% identity with that of A. nidulans alpA and Dis1/Stu2/Xmap215 homolog of N. crassa respectively, which are two folds higher than that of S. cerevisiae Stu2 and S. pombe Dis1 (Figure 1C). Moreover, The carboxyl terminal part of FgStu2p (position 506 to 877) shares 46 and 52% identity with that of A. nidulans alpA (position 510 to 891) and N. crassa stu2 homolog (position 509 to 941) which are nearly three times higher than the 18% identity with that of S. cerevisiae Stu2 (position 567 to 888) and S. pombe Dis1 (position 539 to 882) (Supplementary Figure 2).

FgStu2p Is a Microtubule- and Kinetochore-Associated Protein

The Fgβ-tubulin-GFP and RFP-FgStu2 double labeling experiment indicated that most of the FgStu2p was localized at cytoplasm and some of them were co-localized with free tubulin heterodimers which were not polymerized into microtubules. Another part of FgStu2p was co-localized with microtubules which contain Fgβ1-tubulin or Fgβ2-tubulin (Figure 2A). However, we could not identify the specific location of FgStu2p in the microtubule such as the microtubule plus end. The result of co-immunoprecipitation (Co-IP) assay demonstrated that the FgStu2p binds to tubulin heterodimers which contain Fgβ1-tubulin and Fgβ2-tubulin (Figure 2B). In addition, our Co-IP assay also suggested that FgStu2p binds to Fgγ-tubulin and FgNdc80p, component of kinetochore, which was previously reported in fission yeast (Figure 3; Hsu and Toda, 2011).

Construction of FgStu2 Promoter Replacement Mutant in F. graminearum

After attempted for several times, we failed to obtain FgStu2 gene deletion mutant by homologous recombination method (data not shown). We therefore used promoter replacement strategy to study the genetic function of FgStu2 as previously described (Figure 4A). The transformants were identified by PCR and southern blot. The results indicate that the promoter of FgStu2 gene has been replaced by G418 resistance gene and zear promoter (Figures 4B,C). The expression level of FgStu2 promoter replacement mutants (FgStu2-Si mutants) was significantly decreased compared with those of wild-type strain PH-1 and Fgβ2F167Y strain. Addition of 30 μM ZEA partially restored the FgStu2 gene expression level (Figure 4D). These results indicate that we have successfully obtained the FgStu2-Si mutants.

FgStu2 Gene Is Involved in Hyphal Morphology, Vegetative Growth, and Microtubule Polymerization of F. graminearum

The FgStu2-Si mutants produced twisted hyphae compared with that of the wild-type strain PH-1. Addition of 30 μM ZEA restored the hyphal morphology (Figure 5A). Moreover, the FgStu2-Si mutants grew significantly slower than the wild-type strain PH-1 on PDA, MM, and CM plates. Addition of 30 μM ZEA partially restored the growth rate (Figure 5B and Table 2). Besides, microtubule network was sparse in FgStu2-Si mutants, which was recovered by addition of 30 μM ZEA (Figure 6).

TABLE 2

Table 2. Hyphal growth rate and sensitivity of wild-type strain PH1, FgStu2 promoter replacement mutants to MBC on PDA with or without 30 μM zearalenone (ZEA).

FIGURE 5

Figure 5. Roles of FgStu2 gene in hyphal and colonial morphology. (A) Hyphal morphology was photographed after strains were incubated on potato dextrose agar (PDA) plants with or without 30 μM zearalenone (ZEA) for 24 h at 25°C. Bar = 50 μm. (B) Colony morphology was photographed after strains were inoculated on potato dextrose agar (PDA), complete medium (CM), and minimal medium (MM) plants with or without 30 μM ZEA for 3 days at 25°C.

FIGURE 6

Figure 6. Roles of FgStu2 gene in microtubule polymerization. Fresh spores (10⁵) of FgStu2-Si-β1EGFP and FgStu2-Si-β2EGFP were first incubated into 20 ml of yeast extract-peptone-dextrose medium (YEPD) with or without 30 μM zearalenone (ZEA) and then put it on an orbital shaker (175 rpm) at 25°C to germinate for 10 h and then spotted 5 μl of culture onto the slides for confocal micro examination. Bar = 25 μm.

FgStu2-Si Mutants Had Slightly Decreased MBC Sensitivity

We tested MBC sensitivity of FgStu2 promoter replacement mutants. The effect concentration against 50% mycelia growth (EC₅₀ value) of FgStu2-Si mutants was slightly decreased compared with those of wild-type strain PH-1 (Table 2). However, the FgStu2-Si mutants cannot grow on PDA containing 0.8 μg/ml MBC. When the PDA plates was added 30 μM ZEA, the EC₅₀ values FgStu2-Si mutants were restored and even a bit higher than that of wild-type strain PH-1 (Supplementary Figure 6A). The same phenomenon was observed in FgStu2 promoter replacement mutant generated from Fgβ2F167Y (Fgβ2F167Y-Stu2-Si), an MBC resistant strain reported previously. The EC₅₀ value of Fgβ2F167Y-Stu2-Si (4.28 ± 0.32 μg/ml) was decreased about 30% compared with that of parental strain Fgβ2F167Y (6.48 ± 0.27 μg/ml) (Table 3). Moreover, the Fgβ2F167Y-FgStu2-Si cannot grow on PDA plates containing 20 μg/ml MBC (Supplementary Figure 6B).

TABLE 3

Table 3. Hyphal growth rate and sensitivity of Fgβ2F167Y and Fgβ2F167Y-Stu2-Simutant to MBC^∗ on PDA with or without 30 μM zearalenone (ZEA).

The FgStu2 Gene Is Involved in Sexual and Asexual Reproduction of F. graminearum

The FgStu2 gene promoter replacement mutant showed curved conidia and decreased number of conidial production compared with those of the wild-type strain PH-1. While the conidial morphology and number of conidial production were restored after addition of 30 μM ZEA (Figures 7A,C). Moreover, the FgStu2-Si mutants lost ability to produce perithecia. Addition of 30 μM ZEA restored the ability of perithecia production (Figure 7B).

FIGURE 7

Figure 7. Roles of FgStu2 gene in sexual and asexual reproduction. (A) Conidial morphology of wild-type strain PH-1, FgStu2-Si mutants. The conidia of all strains were harvested after 3-days cultivation in mung bean liquid (MBL) medium with or without 30 μM zearalenone (ZEA) at 25°C. Fresh conidia were first stained with calcofluor white (CFW) and 4’,6-diamidino-2-phenylindole (DAPI) and then photographed using a confocal laser-scanning microscope. Bar = 50 μm. (B) Ability of perithecia production of wild-type strain PH-1, FgStu2-Si-1, FgStu2-Si-2. Inoculate fresh mycelia of each stain to carrot medium plates with or without 30 M m ZEA and keep it at 25°C incubator until the plates are full of mycelia. Then clean out the mycelia and add 2.5% Tween 60 solution to the plates. Keep the plates into 25°C incubator with near-UV light for about one or two weeks until the black perithecia were appeared and then photograph the plates by stereomicroscope. Bar = 1000 μm. (C) Numbers of conidial production of wild-type strain PH-1, FgStu2-Si-1, FgStu2-Si-2. Conidial production was measured by hemacytometer after three days incubation in mung bean liquid (MBL) medium with or without 30 μM ZEA. Means of the column with the same letter are not significantly different (P < 0.05).

The FgStu2 Gene Is Involved in Virulence of F. graminearum

We performed the coleoptile infection assay to assess the rule of FgStu2 gene in pathogenicity of F. graminearum. The lesion length of coleoptiles inoculated with conidia of FgStu2-Si mutants was decreased about 35% compared with those of the wild-type strain PH-1 (Figure 8A). However, the conidia of FgStu2-Si mutants collected from MBL medium containing 30 μM ZEA did not cause the same size of lesion (Figures 8A,B).

FIGURE 8

Figure 8. Roles of FgStu2 gene in pathogenicity. (A) Infected wheat coleoptiles were photographed at 7 days after inoculation. The small solid arrows indicate the inoculation sites and the hollow arrow indicates the lesion development direction. (B) Lengths of brown lesions of wheat coleoptiles. Values are means ± standard deviation of three repeated experiments. Means of the column with the same letter are not significantly different (P < 0.05).

The FgStu2 Gene Is Involved in DON Production of F. graminearum

The DON production was significantly decreased in FgStu2-Si mutants compared with that of wild-type strain PH-1 (Figure 9A). However, addition of ZEA didn’t restore the DON production of FgStu2-Si mutants and decreased the DON production of PH-1. Moreover, expression levels of trichothecene biosynthesis pathway genes (Tri gene) including FgTri5 and FgTri6 were significantly decreased in FgStu2-Si mutants (Figure 9B). The localization of Tri1-EGFP, also known as “toxisome,” in FgStu2-Si mutant was similar with that of wild-type strain PH-1 (Figure 9C).

FIGURE 9

Figure 9. Roles of FgStu2 gene in DON production. (A) Relative DON production of FgStu2-Si mutants. Wild-type strain PH-1 and FgStu2-Si mutants were inoculated into TBI medium with or without 30 μM zearalenone (ZEA) and cultivate for 7 days at 28 °C in the dark. DON production was tested using enzyme linked immunosorbent assay (ELISA). Means come from three independent assays. (B) Relative expression levels of FgTri5 and FgTri6 genes in FgStu2-Si mutants. Wild-type strain PH-1 was used as control. Wild-type strain PH-1 and FgStu2-Si mutants were inoculated into TBI medium with or without 30 μM ZEA and cultivate for 3 days at 28 °C in the dark. The fresh mycelium was used for total RNA extraction. Means come from three independent assays. (C) FgTri1-EGFP localization in wild-type strain PH-1 and FgStu2-Si mutants. PH-1-Tri1-EGFP and FgStu2-Si-Tri1EGFP strains were inoculated into TBI medium with or without 30 μM ZEA and cultivate for 2 days at 28 °C in the dark. The fresh mycelium was observed using a Leica SP8 laser confocal microscope. Bar = 25 μm. This experiment was repeated three times.

Discussion

In this study, we identified the FgStu2 gene (FGSG_10528) is the Dis1/Stu2/XMAP215 microtubule associated protein family homolog of F. graminearum. The amino acid sequences of two FgStu2-TOG domains share 63 and 74% identity with that of A. nidulans and N. crassa, which is much higher than the 26 and 33% identity with ScStu2 and SpDis1 (Figure 1B). The similar result was also reflected in the sequence alignment of carboxyl terminal area. These results indicate that the FgStu2 is more homologous to Dis1/Stu2/XMAP215 homologs in filamentous fungi than yeast (Supplementary Figure 1). Moreover, we tried to substitute FgStu2 in wild-type strain PH-1 with ScStu2. However, after several attempts, we didn’t acquire the substitution mutant which means ScStu2 might not exert the function of FgStu2. This phenomenon may due to the different tubulin structures and sequences of F. graminearum and yeast.

Previous studies demonstrated that Dis1/Stu2/XMAP215 homologs binds to tubulin and microtubule (Al-Bassam and Chang, 2011; Widlund et al., 2011; Thawani et al., 2018). We therefore examined whether FgStu2p interacts with Fgβ-tubulin physically. The co-localization experiment and co-immunoprecipitation demonstrated that FgStu2p binds to Fgβ1-tubulin and Fgβ2-tubulin and partly localizes to microtubule in hyphae of F. graminearum (Figures 2A,B). However, due to the microtubule network in F. graminearum is more complicated than that of yeast, we couldn’t confirm the microtubule plus end tracking property of FgStu2 under our experiment condition. At least, we demonstrated the microtubule associated property of FgStu2p in F. graminearum.

In order to investigate the functional role of FgStu2 gene in F. graminearum, we firstly tried to disrupt FgStu2 gene by gene replacement strategy. However, we cannot acquire correct transformants after several attempts (data not shown), which means FgStu2 gene might be essential for F. graminearum. This phenomenon is consistent with the result of Stu2 disruption experiment in S. cerevisiae (Wang and Huffaker, 1997; Al-Bassam et al., 2006). However, deletion of XMAP215 homologs in fission yeast and A. nidulans is viable (Enke et al., 2007; Hsu and Toda, 2011). These results indicate that functional importance of XMAP215 homologs in different fungal species is divergent. Another possibility is that some species may contain more than one XMAP215 homolog. For example, fission yeast contains two XMAP215 homologs: Dis1 and Alp14 genes, delete any one of them would not be lethal (Hsu and Toda, 2011; Flor-Parra et al., 2018). Therefore, we preformed promoter replacement strategy to study the FgStu2 gene. The promoter replacement mutant FgStu2-Si exhibited drastically decreased growth rate compared with that of wild-type strain PH-1. According to a previous review, loss functions of the XMAP215 homologs leads to decreased growth rates in many species (Al-Bassam and Chang, 2011). Our results are identical with these studies which indicate that FgStu2 gene is essential for the viability of F. graminearum.

Besides exhibited decreased growth rate, the FgStu2-Si also produced twisted hyphae and curved conidia (Figures 5A, 7A). According to previous studies, benomyl (a kind of benzimidazole fungicide, which can destabilize microtubules) treatment can destabilize microtubule, affect hyphal tip growth of A. nidulans, and develop twisted hyphae (Horio and Oakley, 2005); disrupt the FgEB1 (Ending binding one protein homolog in F. graminearum, another kind of conserved microtubule-associated protein) gene resulted in less organized microtubule and therefore partly caused the twisted hyphae and curved conidia (Liu et al., 2017); the Fgα1-tubulin and Fgβ2-tubulin gene deletion mutants displayed twisted hyphae and abnormal conidia (Hu et al., 2015; Zhu et al., 2018). These results suggest that microtubule plays a vital role in morphology of filamentous fungi. Moreover, the alpA gene (Dis1/Stu2/XMAP215 homolog in A. nidulans) deletion mutant exhibited reduced microtubule polymerization and twisted hyphae (Enke et al., 2007). We therefore wonder whether the aberrant hyphae and conidia of FgStu2-Si mutant are caused by abnormal microtubule polymerization. As our expected, The FgStu2-Si mutant exhibited less organized microtubule (Figure 6). This result is consistent with other studies about XMAP215 homologs in other species (Al-Bassam and Chang, 2011; Thawani et al., 2018). We therefore conclude that FgStu2 regulates the morphology of F. graminearum via microtubule associated function. Although the microtubule polymerization was decreased in FgStu2-Si mutants, the FgStu2-Si mutant didn’t show highly increased sensitivity to microtubule destabilizing fungicide MBC (Tables 2, 3 and Supplementary Figure 6). This result indicates that the existed microtubule in FgStu2-Si did not exhibit highly increased sensitivity to microtubule-destabilizing effect of MBC. Previous study showed that AlpA deletion mutant of A. nidulans exhibited increased sensitivity to benomyl (Enke et al., 2007). However, in this study, we didn’t acquire the deletion mutant of FgStu2. The FgStu2 expression was not completely eliminated in FgStu2-Si mutants (Figure 4D), which may still contribute to stabilize microtubule of FgStu2-Si mutant. Another possibility is that the microtubule polymerization promoting activity of FgStu2p is too limited to resist the microtubule-depolymerizing activity of MBC. Although FgStu2-Si mutants are not hypersensitive to MBC, they was unable to grow on PDA plates containing 0.8 μg/ml MBC and Fgβ2F167Y-Stu2-Si also showed increased sensitivity to MBC (Supplementary Figure 6). Therefore, FgStu2 plays a role, but not a major role, in regulation of MBC sensitivity of F. graminearum.

Besides the abnormal conidial morphology, the numbers of conidia produced by FgStu2-Si mutants were decreased compared with that of wild-type strain PH-1. Moreover, FgStu2-Si mutants did not produce perithecia. These results indicate that repression of FgStu2 gene affects both the sexual and asexual reproduction. Therefore, the cytokinesis of F. graminearum must be disturbed in FgStu2-Si mutants. According to previous studies, microtubule nucleation is vital to formation of spindle microtubules and the XMAP215 involves microtubule nucleation by interaction with γ-tubulin ring complex (Thawani et al., 2018); interaction between spindle microtubule and kinetochore is important for chromosome segregation during cytokinesis and the Stu2 regulates this interaction through interacting with Ndc80 complex of kinetochore (Miller et al., 2016); the Dis1 gene is required for proper meiotic chromosome disjunction in fission yeast (Rockmill and Fogel, 1988); In Drosophila, the msps gene (XMAP215 homolog in Drosophila) mutant displayed defect in chromosome segregation (Cullen and Ohkura, 2001); In human HELA cell line, deletion of XMAP215 homolog TOGp affects the integrity of centrosomes and spindle poles and delays mitosis (Cassimeris and Morabito, 2004). These studies indicate that loss functions of Dis1/Stu2/XMAP215 homologs in different species will disturb cytokinesis which is similar with our results. Moreover, we also demonstrated that the FgStu2 interacts with Fgγ-tubulin and presumed FgNdc80 which suggests that FgStu2 may involve the cytokinesis associated functions of FgNdc80 and Fgγ-tubulin just like XMAP215 homologs in other species.

Previous studies showed that deletion mutants of many genes in F. graminearum, such as Fgα1-tubulin, Fgβ2-tubulin, and FgTBCA (Tubulin folding cofactor A homolog in F. graminearum) etc., have dramatically decreased growth rates and are almost incapable of infecting host (Hu et al., 2015; Zhang et al., 2015; Zhu et al., 2018). Similar with these studies, the mean growth rates of FgStu2-Si mutants (0.5 cm/day) were decreased about 80% compared with that of wild-type strain PH-1 (2.6 cm/day) (Figure 5B and Table 2). However, the mean lengths of lesion caused by the FgStu2-Si mutants (1 cm) were decreased only about 35% compared with that of wild-type strain PH-1 (1.5 cm), which indicates that the FgStu2-Si mutants are still capable of infecting host (Figure 8). We considered the growth of FgStu2-Si mutants might be recovered during invasive growth in coleoptile. A previous study showed that the expression level of many genes of F. graminearum were upregulated during invasive growth in planta. Among them, the expression level of FGSG_04581 gene of PH-1 growing in planta has upregulated 6.79 times compared with that of PH-1 growing in vitro (Zhang et al., 2012). According to another study, the expression level of the FGSG_04581 gene is regulated by the zearalenone (Lee et al., 2010). Therefore, the zearalenone concentration in the hyphae maybe increased during infection. Interestingly, the expression of the FgStu2 gene in FgStu2-Si mutants is control by the same promoter of FGSG_04581 gene, which indicates that the expression of FgStu2 gene was likely to be restored by the increased zearalenone concentration and therefore restored the growth and pathogenicity. Nonetheless, we still observed the decreased pathogenicity of FgStu2-Si mutant. The dramatically decreased growth rate would be detrimental to the virulence of F. graminearum. Besides the growth rates, microtubules also play a crucial role in pathogenicity. Previous studies showed that, microtubule is important for transportation of secondary metabolites of plant pathogenic fungi (Conesa et al., 2001; Chanda et al., 2009). In addition, FgEB1 deletion mutants showed less organized microtubule and decreased pathogenicity (Liu et al., 2017). In this study, microtubule polymerization was decreased in FgStu2-Si mutant, which may affect the transportation of virulence-associated factors such as effective protein and DON and therefore contributed to the decreased pathogenicity. However, spores of FgStu2-Si mutant collected from MBL medium containing zearalenone did not caused the same size of lesion compared with that of PH-1 which is similar with the FgHSP90 knockdown mutant (Bui et al., 2016). This is likely due to the FgStu2p was not enough during invasion and therefore showed decreased pathogenicity. Future study should find another way to replace the method of addition of zearalenone.

The DON was wildly considered as an important virulence factor of F. graminearum (Chen et al., 2019). In this study, the DON production of FgStu2-Si mutants were significantly decreased which means the FgStu2 gene is involved the DON production (Figure 9A). Previous studies showed that the expression level of Tri genes were associated with the DON production (Chen et al., 2019). In this study, the expression level of FgTri5 and FgTri6 were decreased in FgStu2-Si mutants, which is consistent with previous studied. However, the FgTri1 localization, which is known as “toxisomes,” has not been affected in FgStu2-Si mutant (Menke et al., 2013), which is different from the FgTri1 localization of Fgβ2-tubulin and FgEB1 deletion mutants (Liu et al., 2017; Wang et al., 2019). This result indicates that different microtubule-associated proteins may play divergent roles in DON synthesis pathway, which needs further investigation. Nonetheless, the decreased microtubule polymerization FgStu2-Si mutant would disturb the transportation of DON, which is similar with FgEB1 deletion mutants. It is strange that addition of zearalenone decreased the DON production of wild-type strain PH-1 and didn’t restore the DON production of FgStu2-Si mutant and PH-1. This phenomenon indicates that addition of zearalenone might result in a negative regulation of DON synthesis of F. graminearum.

In summary, this study revealed the key roles of FgStu2 gene in phenotypes, such as vegetative growth, morphology, sexual reproduction, asexual reproduction, pathogenicity and DON production of F. graminearum. We also demonstrated the microtubule-associated property of FgStu2, which is vital to the vegetative growth and morphology. Our results supplemented the functional roles of Dis1/Stu2/XMAP215 homologs in plant pathogenic fungi, which enriched the theory of Dis1/Stu2/XMAP215 family. Future studies should go deep into the role of FgStu2 and other microtubule-associated proteins in DON and ZEA production, which could be useful for developing new kinds of fungicide to control the mycotoxin contamination.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

Author Contributions

YYZ, YSZ, and MZ conceived and designed the experiments. NL, WR, YH, and XS guided the methods. YYZ and YSZ performed the experiments, analyzed the data, and prepared the figures and tables. YYZ and MZ wrote the manuscript. YSZ, WR, YH, YD, and MZ revised the manuscript. All authors have read and approved the final manuscript.

Funding

This research was supported by the Key Program of National Natural Science Foundation of China (31730072).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2020.545015/full#supplementary-material

Supplementary Figure 1 | The phylogenetic analysis of Dis1/Stu2/XMAP215 homologs. The phylogenetic tree was generated by MEGA 7.0 software. The whole sequences of Dis1/Stu2/XMAP215 homologs (Aaspergillus nidulans protein ID: XP_663125.1; Fusarium graminearum protein ID: XP_011319522.1; Neurospora crassa protein ID: XP_956946.3; Saccharomyces cerevisiae protein ID: NP_013146.1; and Schizosaccharomyces pombe protein ID: NP_587785.1) were acquired from the National Center for Biotechnology information website (NCBI, https://www.ncbi.nlm.nih.gov/).

Supplementary Figure 2 | Multiple Sequences alignment of C-terminal areas of Stu2 homologs. Multiple sequences (Aaspergillus nidulans protein ID: XP_663125.1; Fusarium graminearum protein ID: XP_011319522.1; Neurospora crassa protein ID: XP_956946.3; Saccharomyces cerevisiae protein ID: NP_013146.1; and Schizosaccharomyces pombe protein ID: NP_587785.1) alignment was performed using the online service (www.genome.jp/tools-bin/ clustalw). The alignment images were generated by the online service of ESPript 3.0 (espript.ibcp.fr/ESPript/cgi-bin/ESPript.cgi).

Supplementary Figure 3 | Sequences alignment of budding yeast Ndc80 and presumed FgNdc80. Amino acid sequences (budding yeast Ndc80 protein ID: NP_012122.3; FgNdc80 protein ID: XP_011328506.1) alignment was performed using the online service (www.genome.jp/tools-bin/clustalw). The alignment images were generated by the online service of ESPript 3.0 (espript.ibcp.fr/ESPript/cgi-bin/ESPript.cgi).

Supplementary Figure 4 | Construction of the FgStu2-3 × FLAG strain. (A) Schematic of vector construction. (B) Polymerase chain reaction (PCR) identification of FgStu2-3 × FLAG strains. a 1.8 kb fragment was amplified by primers P29/P30 indicates that 3 × FLAG-G418 fragment was successfully inserted at 3′-terminal of FgStu2 coding sequence.

Supplementary Figure 5 | Construction of the Fgγ-tubulin-EGFP and FgNdc80-EGFP strains. (A) Schematic of vector construction. (B) Polymerase chain reaction (PCR) identification of Fgγ-tubulin-EGFP and FgNdc80-EGFP strains. a 2.7 kb fragment was amplified by primers P35/P36 indicates that the EGFP-HPH fragment was successfully inserted at 3′-terminal of Fgγ-tubulin coding sequence. a 2.7 kb fragment was amplified by primers P41/P42 indicates that the EGFP-HPH fragment was successfully inserted at 3′-terminal of presumed FgNdc80 coding sequence.

Supplementary Figure 6 | MBC Sensitivity of wild-type strain PH-1, Fgβ 2F167Y, and FgStu2-Si mutants. (A) MBC Sensitivity of wild-type strain PH-1 and FgStu2-Si-1 on PDA plates with or without 30 μM zearalenone (ZEA). (B) MBC Sensitivity of Fgβ2F167Y and F167Y-Stu2-Si on PDA plates with or without 30 μM ZEA. Plugs were transferred from the margins of 3-day old PDA colonies to fresh PDA plants containing different concentrations of MBC. The timing for photography was based on the growth rates of each strain.

Supplementary Table 1 | Primers used in this study.

References

Aher, A., and Akhmanova, A. (2018). Tipping microtubule dynamics, one protofilament at a time. Curr. Opin. Cell Biol. 50, 86–93. doi: 10.1016/j.ceb.2018.02.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Al-Bassam, J., and Chang, F. (2011). Regulation of microtubule dynamics by TOG-domain proteins XMAP215/Dis1 and CLASP. Trends Cell Biol. 21, 604–614. doi: 10.1016/j.tcb.2011.06.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Al-Bassam, J., van Breugel, M., Harrison, S. C., and Hyman, A. (2006). Stu2p binds tubulin and undergoes an open-to-closed conformational change. J. Cell Biol. 172, 1009–1022. doi: 10.1083/jcb.200511010

PubMed Abstract | CrossRef Full Text | Google Scholar

Ayaz, P., Ye, X. C., Huddleston, P., Brautigam, A. C., and Rice, M. L. (2012). A TOG: αβ-tubulin complex structure reveals conformation-based mechanisms for a microtubule polymerase. Science 337, 857–860. doi: 10.1126/science.1221698

PubMed Abstract | CrossRef Full Text | Google Scholar

Bai, G. H., and Shaner, G. (2004). Management and resistance in wheat and barley to Fusarium head blight. Annu. Rev. Phytopathol. 42, 135–161. doi: 10.1146/annurev.phyto.42.040803.140340

PubMed Abstract | CrossRef Full Text | Google Scholar

Blandino, M., Haidukowski, M., Pascale, M., Plizzari, L., Scudellari, D., and Reyneri, A. (2012). Integrated strategies for the control of Fusarium head blight and deoxynivalenol contamination in winter wheat. Field Crop Res. 133, 139–149. doi: 10.1016/j.fcr.2012.04.004

CrossRef Full Text | Google Scholar

Bruno, K. S., Tenjo, F., Li, L., Hamer, J. E., and Xu, J. R. (2004). Cellular localization and role of kinase activity of PMK1 in Magnaporthe grisea. Eukaryot. Cell. 3, 1525–1532. doi: 10.1128/EC.3.6.1525-1532.2004

PubMed Abstract | CrossRef Full Text | Google Scholar

Bui, D., Lee, Y. J., Lim, J. Y., Fu, M. M., Kim, J., Choi, G. J., et al. (2016). Heat shock protein 90 is required for sexual and asexual development, virulence, and heat shock response in Fusarium graminearum. Sci. Rep. 6:28154. doi: 10.1038/srep28154

PubMed Abstract | CrossRef Full Text | Google Scholar

Cassimeris, L., and Morabito, J. (2004). TOGp, the human homolog of XMAP215/Dis1, is required for centrosome integrity, spindle pole organization, and bipolar spindle assembly. Mol. Biol. Cell. 15, 1580–1590. doi: 10.1091/mbc.e03-07-0544

PubMed Abstract | CrossRef Full Text | Google Scholar

Chanda, A., Roze, L. V., Kang, S., Artymovich, K. A., Hicks, G. R., Raikhel, N. V., et al. (2009). A key role for vesicles in fungal secondary metabolism. Proc. Natl. Acad. Sci. USA. 106, 19533–19538. doi: 10.1073/pnas.0907416106

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, C. J., Yu, J. J., Bi, C. W., Zhang, Y. N., Xu, J. Q., Wang, J. X., et al. (2009). Mutations in a β-tubulin confer resistance of Gibberella zeae to benzimidazole fungicides. Phytopathol. 99, 1403–1411. doi: 10.1094/PHYTO-99-12-1403

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, P. X. Y., Yin, H. W., and Huffaker, T. C. (1998). The yeast spindle pole body component Spc72p interacts with Stu2p and is required for proper microtubule assembly. J. Cell. Biol. 141, 1169–1179. doi: 10.1083/jcb.141.5.1169

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Y., Kistler, H. C., and Ma, Z. H. (2019). Fusarium graminearum trichothecene mycotoxins: biosynthesis, regulation, and management. Annu. Rev. Phytopathol. 57, 1.1–1.25. doi: 10.1146/annurev-phyto-082718-100318

PubMed Abstract | CrossRef Full Text | Google Scholar

Conesa, A., Punt, P. J., Van, L. N., and Ca, V. D. H. (2001). The secretion pathway in filamentous fungi: a biotechnological view. Fungal Genet. Biol. 33, 155–171. doi: 10.1006/fgbi.2001.1276

PubMed Abstract | CrossRef Full Text | Google Scholar

Cullen, C. F., and Ohkura, H. (2001). Msps protein is localized to acentrosomal poles to ensure bipolarity of Drosophila meiotic spindles. Nat. Cell Biol. 3, 637–642. doi: 10.1038/35083025

PubMed Abstract | CrossRef Full Text | Google Scholar

Cuomo, C. A., Güldener, U., Xu, J. R., Trail, F., Turgeon, B. G., Pietro, A. D., et al. (2007). The Fusarium graminearum genome reveals a link between localized polymorphism and pathogen specialization. Science 317, 1400–1402. doi: 10.1126/science.1143708

PubMed Abstract | CrossRef Full Text | Google Scholar

Davidse, L. C. (1986). Benzimidazole fungicides: Mechanism of action and biological impact. Annu. Rev. Phytopathol. 24, 43–65. doi: 10.1146/annurev.py.24.090186.000355

CrossRef Full Text | Google Scholar

Enke, C., Zekert, N., Veith, D., Schaaf, C., Konzack, S., and Fisher, R. (2007). Aspergillus nidulans Dis1/XMAP215 protein alpA localizes to spindle pole bodies and microtubule plus ends and contributes to growth directionality. Eukaryot. Cell. 6, 555–562. doi: 10.1128/EC.00266-06

PubMed Abstract | CrossRef Full Text | Google Scholar

Flor-Parra, P., Iglesias-Romero, A. B., and Chang, F. (2018). The XMAP215 ortholog Alp14 promotes microtubule nucleation in fission yeast. Curr. Biol. 28, 1681–1691. doi: 10.1016/j.cub.2018.04.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Forsyth, D., Yoshizawa, T., Morooka, N., and Tuite, J. (1977). Emetic and refusal activity of 632 deoxynivalenol to swine. Appl. Environ. Microbiol. 34, 547–552. doi: 10.1007/BF02313762

PubMed Abstract | CrossRef Full Text | Google Scholar

Gard, D., and Kirschner, M. (1987). A microtubule-associated protein from Xenopus eggs that specifically promotes assembly at the plus-end. J. Cell Biol. 105, 2203–2215. doi: 10.1083/jcb.105.5.2203

PubMed Abstract | CrossRef Full Text | Google Scholar

Gardiner, D., Kazan, K., and Manners, J. (2009). Nutrient profiling reveals potent inducers of trichothecene biosynthesis in Fusarium graminearum. Fungal Genet. Biol. 46, 604–613. doi: 10.1016/j.fgb.2009.04.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Geyer, A. E., Miller, P. M., Brautigam, A. C., Biggins, S., and Rice, M. L. (2018). Design principles of a microtubule polymerase. eLife 7:e34574. doi: 10.7554/eLife.34574

PubMed Abstract | CrossRef Full Text | Google Scholar

Gu, Q., Zhang, C. Q., Yu, F. W., Yin, Y. N., Shim, W., and Ma, Z. H. (2015). Protein kinase FgSch9 serves as a mediator of the target of rapamycin and high osmolarity glycerol pathways and regulates multiple stress responses and secondary metabolism in Fusarium graminearum. Environ. Microbiol. 17, 2661–2676. doi: 10.1111/1462-2920.12522

PubMed Abstract | CrossRef Full Text | Google Scholar

Horio, T., and Oakley, B. R. (2005). The role of microtubules in rapid hyphal tip growth of Aspergillus nidulans. Mol. Biol. Cell. 16, 918–926. doi: 10.1091/mbc.e04-09-0798

PubMed Abstract | CrossRef Full Text | Google Scholar

Hsu, K., and Toda, T. (2011). Ndc80 internal loop interacts with Dis1/TOG to ensure proper kinetochore-spindle attachment in fission yeast. Curr. Biol. 21, 214–220. doi: 10.1016/j.cub.2010.12.048

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, W. Q., Zhang, X. P., Chen, X., Zheng, J. W., Yin, Y. N., and Ma, Z. H. (2015). The α1-tubulin FaTuA1 plays crucial roles in vegetative growth and conidiation in Fusarium asiaticum. Res. Microbiol. 166, 132–142. doi: 10.1016/j.resmic.2015.01.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, C., Hei, R. N., Yang, Y., Zhang, S. J., Wang, Q. H., Wang, W., et al. (2020). An orphan protein of Fusarium graminearum modulates host immunity by mediating proteasomal degradation of TaSnRK1α. Nat. Commun. 11, 4382. doi: 10.1038/s41467-020-18240-y

CrossRef Full Text | Google Scholar

Lee, H. J., and Ryu, D. (2017). Worldwide occurrence of mycotoxins in cereals and cereal-derived food products: public health perspectives of their co-occurrence. J. Agric. Food Chem. 65, 7034–7051. doi: 10.1021/acs.jafc.6b04847

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, J., Son, H., Lee, S., Park, A. R., and Lee, Y. (2010). Development of a conditional gene expression system using a zearalenone-inducible promoter for the Ascomycete Fungus Gibberella zeae. Appl. Environ. Microbiol. 76, 3089–3096. doi: 10.1128/AEM.02999-09

PubMed Abstract | CrossRef Full Text | Google Scholar

Letunic, I., Doerks, T., and Bork, P. (2014). SMART: recent updates, new developments and status in 2015. Nucleic Acids Res. 43, 257–260. doi: 10.1093/nar/gku949

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Z. Y., Wu, S. S., Chen, Y., Han, X. Y., Gu, Q., Yin, Y. N., et al. (2017). The microtubule end-binding protein FgEB1 regulates polar growth and fungicide sensitivity via different interactors in Fusarium graminearum. Environ. Microbiol. 19, 1791–1807. doi: 10.1111/1462-2920.13651

PubMed Abstract | CrossRef Full Text | Google Scholar

Menke, J., Weber, J., Broz, K., and Kistler, H. C. (2013). Cellular development associated with induced mycotoxin synthesis in the filamentous fungus fusarium graminearum. PloS One. 8:e63077. doi: 10.1371/journal.pone.006307

CrossRef Full Text | Google Scholar

Miller, P. M., Asbury, L. C., and Biggins, S. (2016). A TOG protein confers tension sensitivity to kinetochore-microtubule attachments. Cell 165, 1428–1439. doi: 10.1016/j.cell.2016.04.030

PubMed Abstract | CrossRef Full Text | Google Scholar

Nogales, E. (2001). Structural insights into microtubule function. Annu. Rev. Biochem. 69, 397–420. doi: 10.1146/annurev.biophys.30.1.397

PubMed Abstract | CrossRef Full Text | Google Scholar

Qiu, J. B., Huang, T. T., Xu, J. Q., Bi, C. W., Chen, C. J., and Zhou, M. G. (2012). β-Tubulins in Gibberella zeae: Their characterization and contribution to carbendazim resistance. Pest Manage. Sci. 68, 1191–1198. doi: 10.1002/ps.3283

PubMed Abstract | CrossRef Full Text | Google Scholar

Ren, W. C., Zhang, Z. H., Shao, W. Y., Yang, Y. L., Zhou, M. G., and Chen, C. J. (2014). The autophagy-related gene BcATG1 is involved in fungal development and pathogenesis in Botrytis cinerea. Mol. Plant. Pathol. 18, 238–248. doi: 10.1111/mpp.12396

PubMed Abstract | CrossRef Full Text | Google Scholar

Robert, X., and Gouet, P. (2014). Deciphering key features in protein structures with the new ENDscript server. Nucleic Acids Res. 42, 320–324. doi: 10.1093/nar/gku316

PubMed Abstract | CrossRef Full Text | Google Scholar

Rockmill, B., and Fogel, S. (1988). DIS1: A Yeast Gene Required for Proper Meiotic Chromosome Disjunction. Genetics 119, 261–272.

Google Scholar

Schuyler, S., and Pellman, D. (2001). Microtubule “Plus-End-Tracking Proteins”: The end is just the beginning. Cell 205, 421–424. doi: 10.1016/S0092-8674(01)00364-6

CrossRef Full Text | Google Scholar

Tang, G. F., Chen, Y., Xu, J. R., Kistler, H. C., and Ma, Z. H. (2018). The fungal myosin I is essential for Fusarium toxisome formation. PLoS Pathog. 14:e1006827. doi: 10.1371/journal.ppat.1006827

PubMed Abstract | CrossRef Full Text | Google Scholar

Thawani, A., Kadzik, S. R., and Petry, S. (2018). XMAP215 is a microtubule nucleation factor that functions synergistically with the γ-tubulin ring complex. Nat. Cell. Biol. 20, 575–585. doi: 10.1038/s41556-018-0091-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Tischfield, M. A., Baris, H. N., Wu, C., Rudolph, G., Van Maldergem, L., He, W., et al. (2010). Human TUBB3 mutations perturb microtubule dynamics, kinesin interactions, and axon guidance. Cell 140, 74–87. doi: 10.1016/j.cell.2009.12.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Walczak, C. (2000). Microtubule dynamics and tubulin interacting proteins. Curr. Opin. Cell Biol. 12, 52–56. doi: 10.1016/S0955-0674(99)00056-3

CrossRef Full Text | Google Scholar

Wang, H., Chen, D. P., Li, C. L., Tian, N., Zhang, J., Xu, J. R., et al. (2019). Stage-specific functional relationships between Tub1 and Tub2 beta-tubulins in the wheat scab fungus Fusarium graminearum. Fungal Genet. Biol. 132:103251. doi: 10.1016/j.fgb.2019.103251

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, P. J., and Huffaker, T. C. (1997). Stu2p: a microtubule-binding protein that is an essential component of the yeast spindle pole body. J. Cell. Biol. 139, 1271–1280. doi: 10.1083/jcb.139.5.1271

PubMed Abstract | CrossRef Full Text | Google Scholar

Widlund, P. O., Stear, J. H., Pozniakovsky, A., Zanic, M., Reber, S., Brouhard, G. J., et al. (2011). XMAP215 polymerase activity is built by combining multiple tubulin-binding TOG domains and a basic lattice-binding region. PNAS. 108, 2741–2746. doi: 10.1073/pnas.1016498108

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, J. H., Hamari, Z., Han, K. H., Seo, J. A., Reyes-Dominguez, Y., and Scazzocchio, C. (2004). Double-joint PCR: a PCR-based molecular tool for gene manipulations in filamentous fungi. Fungal Genet. Biol. 41, 973–981. doi: 10.1016/j.fgb.2004.08.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Zanic, M., Widlund, O. P., Hyman1, A. A., and Howard, J. (2013). Synergy between XMAP215 and EB1 increases microtubule growth rates to physiological levels. Nat. Cell Biol. 15, 688–693. doi: 10.1038/ncb2744

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, X. P., Chen, X., Jiang, J. H., Yu, M. H., Yin, Y. N., and Ma, Z. H. (2015). The tubulin cofactor A is involved in hyphal growth, conidiation and cold sensitivity in Fusarium asiaticum. BMC Microbiol. 15:35. doi: 10.1186/s12866-015-0374-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, X. W., Jia, L. J., Zhang, Y., Jiang, G., Li, X., Zhang, D., et al. (2012). In planta stage-specific fungal gene profiling elucidates the molecular strategies of Fusarium graminearum growing inside wheat coleoptiles. Plant Cell 24, 5159–5176. doi: 10.1105/tpc.112.105957

PubMed Abstract | CrossRef Full Text | Google Scholar

Zheng, Z. T., Gao, T., Zhang, Y., Hou, Y. P., Wang, J. X., and Zhou, M. G. (2014). FgFim, a key protein regulating resistance to the fungicide JS399-19, asexual and sexual development, stress responses and virulence in Fusarium graminearum. Mol. Plant Pathol 15, 488–499. doi: 10.1111/mpp.12108

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, Y. Y., Liang, X. Y., Li, Y. J., Duan, Y. B., Zheng, Z. T., Wang, J. X., et al. (2018). F240 of β2-tubulin explains why Fusarium graminearum is less sensitive to carbendazim than Botrytis cinerea. Phytopathology 108, 352–361. doi: 10.1094/PHYTO-09-17-0295-R

PubMed Abstract | CrossRef Full Text | Google Scholar