- 1Key Lab of Integrated Crop Pest Management of Shandong, College of Plant Health and Medicine, Qingdao Agricultural University, Qingdao, China
- 2Shandong Province Key Laboratory of Applied Mycology, Qingdao Agricultural University, Qingdao, China
- 3State Key Laboratory of Microbial Metabolism, Shanghai Jiao Tong University, Shanghai, China
BcHpt is a core element of the high-osmolarity glycerol (HOG) transduction pathway in Botrytis cinerea. In contrast to other elements of the pathway, which have been characterized and proven to play important roles in vegetative differentiation, fungicide resistance, the multistress response, and virulence in B. cinerea, BcHpt (Histidine-containing phosphotransfer) is essential but uncharacterized in B. cinerea. Our previous study reported the first lysine acetylation site (Lys161) in BcHpt. In this study, the functions of this lysine acetylation site in BcHpt were characterized using site-directed mutagenesis. To mimic Lys161 acetylation, we generated the mutant strain ΔBcHPt + BcHptK161Q-GFP, which exhibited a slower growth rate; lower pathogenicity; higher sensitivity to multiple stresses, including osmotic and oxidative stresses, dicarboximides, and demethylation inhibitors (DMIs); and lower BcSak1 phosphorylation levels than wild-type B. cinerea. Constitutive acetylation of BcHpt Ly161 apparently inhibits hyphal growth, the multistress response, and sensitivity to fungicides in B. cinerea. Moreover, the lysine acetylation site affected phosphorylation of the MAPK BcSak1.
Introduction
Mitogen-activated protein kinase (MAPK) signaling pathways play important roles in the response of fungal pathogens to various extracellular stresses, including osmotic and oxidative stress (Kruppa and Calderone, 2006). In Saccharomyces cerevisiae, the high-osmolarity glycerol (HOG) pathway is one of the most well-characterized MAP kinase modules (Maeda et al., 1995; Posas and Saito, 1997; Barrett and Hoch, 1998; Koretke et al., 2000; Santos and Shiozaki, 2001). In S. cerevisiae, the HOG pathway contains two branches (Sln1 and Sho1), which converge on the MAPKK Pbs2 (Hohmann, 2002). The Sln1 branch consists of Ypd1, Skn7, Ssk1, and Ssk2/Ssk22. Sln1 can sense extracellular osmotic conditions and transfer different signals to downstream MAPK by sequential phosphorylation (Posas et al., 1996; Krantz et al., 2006). In many filamentous fungi phosphorylation of Ssk1 triggers activation of the downstream component. The Sho1 branch consists of Sho1, Cdc42, Ste20, Ste50, Ste11, and Pbs2 (Posas and Saito, 1997; O’Rourke and Herskowitz, 1998).
In filamentous fungi, several elements of the HOG pathway have been identified: an osmosensor histidine kinase (BcOs1, or Bcbos1 in Botrytis cinerea); a histidine phosphotransfer (Hpt) protein (Ypd1 in S. cerevisiae, BcHpt in B. cinerea); two response regulators (Ssk1 in S. cerevisiae, BRrg1 in B. cinerea); and downstream MAPK cascades containing (Ssk2/Ssk22 in S. cerevisiae, BcOs4 in B. cinerea; Pbs2 in S. cerevisiae, BcOs5 in B. cinerea and Hog1 in S. cerevisiae, BcSak1/BcOs2 in B. cinerea) (Viaud et al., 2006; Noguchi et al., 2007; Segmüller et al., 2007; Liu et al., 2008; Yan et al., 2010, 2011; Fillinger et al., 2012; Yang et al., 2012). Jiang et al. (2018) reviewed MAPK signaling in plant pathogenic fungi. The HOG pathway plays a species-specific and tissue-specific role in fungal virulence (Dixon et al., 1999; Mehrabi et al., 2006; Moriwaki et al., 2006; Jiang et al., 2018). In addition, the HOG pathway is important for mediating osmotic and oxidative stress responses and fungicide (phenylpyrrole and dicarboximide) sensitivity. Many researchers have suggested that the HOG pathway is the target of fludioxonil and fenpiclonil, which are phenylpyrrole fungicides (Jespers et al., 1993; Zhang et al., 2002; Fujimura et al., 2003; Motoyama et al., 2005; Noriyuki et al., 2007; Motoyama et al., 2008; Jespers and Waard, 2010). However, the role of this pathway in responses to other (e.g., cell wall) stresses varies among fungi (Jiang et al., 2018). In addition, the HOG pathway is also involved in the species-specific regulation of hyphal growth and development (Zhang et al., 2002; Jones et al., 2007; Vitalini et al., 2007).
The role of HPt proteins has been characterized in several fungi (Aoyama et al., 2000; Vargas-Pérez et al., 2007; Lee et al., 2011; Mavrianos et al., 2014; Jacob et al., 2015). Numerous fungi contain only one HPt protein, and it seems to integrate all signals from different sensor histidine kinases. Thus, Hpt is essential in many fungi, including S. cerevisiae, N. crassa and Aspergillus nidulans (Maeda et al., 1994; Banno et al., 2007; Furukawa et al., 2010). Some researchers have considered this protein a suitable target for novel antifungal drugs (Fassler and West, 2013).
Botrytis cinerea causes gray mold on over 400 plant species, leading to extreme financial losses worldwide (Williamson et al., 2007; Dean et al., 2012; Fillinger and Elad, 2016). Although most core elements are well characterized in B. cinerea, the function of BcHpt remains unclear. Our previous proteomics study in B. cinerea identified one acetylation site, Lys161, in BcHpt (Lv et al., 2016); this acetylation site was the first to be reported in BcHpt. To determine the role of lysine acetylation in BcHpt, we characterized Lys161 of BcHpt in B. cinerea using site-directed mutagenesis.
Materials and Methods
Strains and Culture Conditions
The standard reference strain B05.10 of B. cinerea Pers. Fr. [Botrytis fuckeliana (de Bary) Whetzel] was isolated from Vitis vinifera (Quidde et al., 1999). All B. cinerea strains used in this study were grown on potato dextrose agar (PDA: 200 g of potato, 20 g of dextrose, 20 g of agar, and 1 L of water).
Conidium and sclerotium formation was assessed after ten days or 4 weeks of incubation on PDA medium. Growth assays were conducted under different stress conditions, and the percentage of mycelial radial growth inhibition (RGI) was measured after 3 days of incubation on PDA as previously described (Yang et al., 2018).
Generation of the BcHpt Mutant Strains by Site-Directed Mutagenesis
The primers used in this study are listed in the Supplementary Table S1. Since the deletion of Hpt1 in B. cinerea is lethal (data not shown), generation of the BcHpt mutant strains was carried out by site-directed mutagenesis using the following protocol: First, the primers which contain the mutated site were designed and listed in Supplementary Table S1 (BcHpt-GFP-F + BcHpt-Q-R and BcHpt-Q-F + BcHpt-GFP-R for the BcHpt-Q-up and BcHpt-Q-down fragments, respectively; BcHpt-GFP-F + BcHpt-R-R and BcHpt-R-F + BcHpt-GFP-R for the BcHpt-R-up and BcHpt-R-down fragments, respectively) and used to amplify the BcHpt gene. Fusion PCR (BcHpt-Q-up and BcHpt-Q-down fragments; BcHpt-R-up and BcHpt-R-down fragments) was employed using BcHpt-GFP-F + BcHpt-GFP-R (Supplementary Table S1) to amplify the BcHptK161Q and BcHptK161R sequences (Yu et al., 2004). The resulting sequences were cotransformed with XhoI-digested pYF11 plasmid into the yeast strain XK1-25 to generate BcHptK161Q/R/K-GFP fusion vectors (Bruno et al., 2004). The resulting vectors: BcHptK161Q-GFP-pYF11, BcHptK161R-GFP-pYF11, and BcHptK161K-GFP-pYF11, were transformed into the B05.10 strain using protoplast formation and transformation of B. cinerea (Gronover et al., 2001; Jiang et al., 2011), and the resulting transformants (named B05.10 + BcHptK161Q-GFP, B05.10 + BcHptK161R-GFP, and B05.10 + BcHptK161K-GFP) were confirmed by PCR (GFP-F and GFP-R for detection of GFP gene), sequencing (BcHpt-SE for detection of site mutation) and Western blotting (using an anti-GFP antibody to confirm the expression of BcHptK161Q/R/-GFP).
Subsequently, the native BcHpt locus in the resulting transformants was deleted by a homologous recombination strategy to generate the mutant ΔBcHPt + BcHPtK161Q, ΔBcHPt + BcHptK161R, and ΔBcHPt + BcHptK161K strains (Supplementary Figure S1). The gene deletion vector was constructed by inserting two flanking sequences (BcHpt-up-F and BcHpt-up-R for BcHpt-up fragment; BcHpt-down-up and BcHpt-down-R for BcHpt-down fragment) of the BcHPT gene into two sides of the HPH (hygromycin resistance) gene in the pBS-HPH1 vector. The resulting vector, pBS-BcHPT-Del, was transformed into B05.10 + BcHptK161Q-GFP, B05.10 + BcHptK161R-GFP, and B05.10 + BcHptK161K-GFP strains using protoplast formation and transformation of B. cinerea. The gene deletion mutants were identified by PCR assays (BcHpt-out-F + BcHpt-out-R), and further confirmed by a Southern blot assay. The probe used in the Southern blot analysis was the BcHPT downstream fragment, which was amplified using BcHpt-down-F and BcHpt-down-R. DNA extracted from the mutants was digested with BamHI. In total of six ΔBcHPt + BcHptK161Q-GFP, four ΔBcHPt + BcHptK161R-GFP, and three ΔBcHPt + BcHptK161K-GFP were obtained, and showed similar phenotypes. Single spore mutants were isolated and transferred to selection medium.
Nucleic Acid Manipulations and qRT-PCR
Genomic DNA was extracted from B. cinerea as previously described (McDonald and Martinez, 1990). Plasmid miniprep purification kits (BioDev Co., Beijing, China) were used to purify plasmid DNA.
The expression levels of BcHpt were tested by qRT-PCR using the 2–ΔΔCt method (Livak and Schmittgen, 2001). Mycelia of the mutants were harvested after 2 days incubation in potato dextrose broth (PDB) in a shaker. RNA was extracted using a protocol described previously (Yang et al., 2018). Reverse transcription was carried out using Revert Aid H Minus First Strand cDNA Synthesis kits (Fermentas Life Sciences, Burlington, Canada). qRT-PCR was conducted using TAKARA SYBR Premix Ex Taq (TAKARA Bio Inc., Dalian, China) with the listed primers (Supplementary Table S1). β-tubulin gene was amplified as a reference. Three biological replicates were used for each sample.
Pathogenicity and Infection-Related Morphogenesis Assays
Pathogenicity testing was performed as previously described (Yang et al., 2013, 2018). An infection-related morphogenesis assay was performed on onion epidermis as previously described (Doehlemann et al., 2006; Viaud et al., 2006).
Western Blot Assay
To confirm the expression of the fusion proteins BcHptK161Q-GFP, BcHptK161R-GFP and BcHptK161K-GFP, Western blotting was carried out using an anti-GFP antibody (Beyotime, Shanghai, China) and anti-acetyl lysine mouse mAb (clone Kac-01, PTM-101) (PTM Biolabs, Hangzhou, China). Conidia were harvested from 10-day old cultures. The mycelia of the mutants were grown in yeast extract peptone dextrose (YEPD) at 25°C for 2 days in a shaker. After the cultures were treated with 20 mM H2O2, 0.2 mg/ml congo red, 1 M NaCl and 1 μg/ml iprodione for 2 h, mycelia were harvested. Protein extraction was carried out as previously described (Gu et al., 2015). The expression of BcSak1 and phosphorylated BcSak1 was examined by using an anti-Hog1 antibody (Santa Cruz Biotechnology, CA, United States) and an antibody against dually phosphorylated p38 (Thr180/Tyr182) (Cell Signaling Technology, MA, United States).
Results
Generation of the BcHpt Mutants by Site-Directed Mutagenesis
In our previous study, we obtained more than 800 ectopic mutants without deletion mutant, thus we believe BcHpt is essential (data not shown). Our previous proteomics study identified one acetylation site, Lys161, in BcHpt (Lv et al., 2016). The Lys161 acetylation site in BcHpt is conserved in 4 of the 11 fungal species analyzed (B. cinerea, Sclerotinia sclerotiorum, Fusarium Graminearum, and Penicillium marneffei) (Figure 1). In B. cinerea, Lys 161 is not inside the HPt domain of BcHpt (66–146) (Figure 1). Deletion of Hpt in B. cinerea and several other fungi was reported to be lethal (Maeda et al., 1994; Banno et al., 2007; Furukawa et al., 2010). Thus, to verify the function of the lysine acetylation site in BcHpt, BcHpt containing either the K161Q or K161R mutation was introduced into the wild-type strain B05.10. Glutamine (Q) and arginine (R) mimic acetylated and unacetylated lysine, respectively (Schwer et al., 2006; Li et al., 2007). In addition, BcHpt-GFP without any point mutations was used as a control. The resulting strains (B05.10 + BcHptK161Q-GFP, B05.10 + BcHptK161R-GFP, and B05.10 + BcHptK161K-GFP) were confirmed by PCR, sequencing and Western blotting using an anti-GFP antibody. Then, the native BcHpt locus was deleted in the transformants using a homologous recombination strategy to generate the mutant ΔBcHPt + BcHptK161Q-GFP, ΔBcHPt + BcHptK161R-GFP, and ΔBcHPt + BcHptK161K-GFP strains. The acetylation levels of the purified BcHptK161Q/R/K-GFP proteins were determined using a pan anti-acetyl lysine antibody. Western blotting showed that ΔBcHPt + BcHptK161Q-GFP and ΔBcHPt + BcHptK161R-GFP were not acetylated but that ΔBcHPt + BcHptK161K-GFP was, indicating that Lys161 is the only acetylation site in BcHpt (Figure 2). In addition, as shown in Supplementary Figure S2, the expression levels of BcHPT in ΔBcHPt + BcHptK161Q-GFP, ΔBcHPt + BcHptK161R-GFP, and ΔBcHPt + BcHptK161K-GFP were similar to that in B05.10. Thus the phenotypic changes were not due to differences in the expression levels.
Figure 1. Multiple sequence alignment of the HPt domain of several fungal Ypd1 orthologs. The genome database accession numbers are as follows: Saccharomyces cerevisiae NP_010046.1, Botrytis cinerea XP_024552997.1, N. crassa XP_011393451.1, F. graminearum XP_011321244.1, S. sclerotiorum SS1G_10394, Cryptococcus neoformans XP_568444.1, Magnaporthe oryzae XP_003715375.1, A. nidulans XP_659609.1, Eremothecium gossypii AAS50589.1, P. marneffei XP_002151180.1, and Trichophyton equinum EGE05000.1. BoxShade was used to highlight identical (black shading) or similar (gray shading) amino acids. The Hpt domain is underlined, and the His64 and Lys161 positions are indicated with arrows.
Figure 2. Acetylation levels of BcHPt in ΔBcHpt + BcHptK161K-GFP, ΔBcHpt + BcHptK161Q-GFP, and ΔBcHpt + BcHptK161R-GFP. BcHptK161Q-GFP, BcHptK161R-GFP, and BcHptK161K-GFP were transformed into B. cinerea, and BcHpt-GFP was pulled down in protein extracts from the mutants using anti-GFP antibody agarose beads. The acetylation and total protein levels were determined using an anti-acetyl lysine antibody and anti-GFP antibody, respectively.
Acetylation of BcHpt1-Lys161 Has Negative Effects on Hyphal Growth but Not on Conidium and Sclerotium Formation in B. cinerea
The mycelial growth rate of ΔBcHPt + BcHptK161K-GFP and ΔBcHPt + BcHptK161R-GFP was similar to that of the wild-type parent B05.10, while ΔBcHPt + BcHptK161Q-GFP showed minor growth retardation compared to B05.10 (Figure 3 and Supplementary Figure S3). After 3 days of incubation, the colony diameters of ΔBcHPt + BcHptK161Q were reduced by 15%. However, the formation of conidia and sclerotia in the three mutants was similar to that in B05.10 after one or 4 weeks of incubation on PDA (Figure 3). These results indicate that acetylation of BcHpt1 Lys161 down-regulates hyphal growth but does not affect conidium and sclerotium formation in B. cinerea.
Figure 3. BcHpt Lys161 is involved in hyphal growth but not in conidium and sclerotium formation in B. cinerea. Morphology of B05.10, ΔBcHpt + BcHptK161Q-GFP, ΔBcHpt + BcHptK161R-GFP, and ΔBcHpt + BcHptK161K-GFP on PDA medium after 3 days, 1 week and 4 weeks of incubation.
BcHpt Lys161 Is Involved in the Fungicide and Stress Sensitivity of B. cinerea
BcHpt is a core element of the HOG pathway, which can be activated by various stresses, including oxidative and osmotic stresses, fungicides and hormones; thus, the sensitivity of the mutants to these stresses was investigated. As shown in Figure 4 and Supplementary Figure S4, compared to the B05.10, ΔBcHPt + BcHptK161R-GFP, and ΔBcHPt + BcHptK161K-GFP strains, ΔBcHPt + BcHptK161Q-GFP exhibited increased sensitivity to iprodione and triadimefon, which belong to dicarboximides and ergosterol biosynthesis inhibitors (EBIs), respectively. Furthermore, ΔBcHPt + BcHptK161Q-GFP was also more sensitive to osmotic (NaCl and KCl), oxidative (H2O2), and cell wall stress (Congo red and caffeine) than the wild-type strain. These results indicate that constitutive acetylation of BcHpt Lys161 increases sensitivity of B. cinerea to fungicides (dicarboximides and DMIs) exposure, and multiple stresses, including oxidative, osmotic, and cell wall stresses.
Figure 4. Sensitivity of B05.10, ΔBcHpt + BcHptK161Q-GFP, ΔBcHpt + BcHptK161R-GFP, and ΔBcHpt + BcHptK161K-GFP to multiple stresses. Inhibition of mycelial growth in the treated strains compared with the growth of the untreated control strains. The bars denote the standard errors from three experiments, and statistical tests were carried out using Tukey’s test for multiple comparisons. Values on the bars followed by different letters are significantly different at P = 0.05.
Regulation of BcSak1 Phosphorylation by BcHpt Lys161
Phosphorylation of BcSak1 indicates the activation of the HOG pathway in B. cinerea; thus, we analyzed the effects of BcHpt Lys161 acetylation on the phosphorylation of BcSak1. Two bands were detected for the phosphorylated Sak1, consistent with the findings of previous studies (Liu et al., 2008; Yang et al., 2013). Western blot analysis showed that BcSak1 phosphorylation was strongly increased in response to 20 mM H2O2 and 0.2 mg/ml congo red in B05.10 and ΔBcHPt + BcHptK161R-GFP but not in ΔBcHPt + BcHptK161Q-GFP (Figure 5). In addition, in response to 1 M NaCl and 1 μg/ml iprodione, BcSak1 phosphorylation was strongly increased in B05.10 but not in ΔBcHPt + BcHptK161Q–GFP (Supplementary Figure S5). These results indicate that constitutive acetylation of BcHpt Lys161 affects the phosphorylation level of BcSak1 especially under activating conditions.
Figure 5. BcHpt Lys161 affects the phosphorylation levels of BcSak1. Comparison of BcSak1 phosphorylation in B05.10, ΔBcHpt + BcHptK161Q-GFP and ΔBcHpt + BcHptK161R-GFP. Phosphorylated and total BcSak1 proteins were detected using anti-phosphorylated p38 (Thr180/Tyr182) and anti-Hog1 antibodies, respectively.
Acetylation of BcHpt in Response to Various Stresses
Figure 5 indicates that BcHpt Lys161 affects the phosphorylation level of BcSak1 in response to various stresses; thus, it would be interesting to know whether the acetylation status of Lys161 can change in response to various stresses. As shown in Figure 6, BcHpt acetylation could be detected in the mycelium, but not in conidia. Furthermore, in response to osmotic (NaCl), oxidative (H2O2), and cell wall (Congo red) stresses, the acetylation levels of BcHpt were significantly decreased.
Figure 6. Acetylation of BcHpt in response to various stresses. Comparison of BcHpt acetylation in ΔBcHpt + BcHptK161K-GFP. Acetylated and total BcHpt proteins were detected using anti-Ac-K and anti-GFP antibodies, respectively.
Constitutive Acetylation of BcHpt Lys161 Affects the Virulence of B. cinerea
Since several components of the HOG pathway are involved in the virulence of B. cinerea, we also performed an infection test on tomato leaves to test the influence of BcHpt Lys161 on B. cinerea pathogenicity. As shown in Figure 7, 3 days after inoculation, ΔBcHPt + BcHptK161Q-GFP caused smaller (approximately 32.3%) lesions on tomato leaves than B05.10, ΔBcHPt + BcHptK161R-GFP and ΔBcHPt + BcHptK161K-GFP (Figure 7). The pathogenicity defects of ΔBcHPt + BcHptK161Q-GFP were analyzed in detail using onion epidermis penetration assays. However, the penetration efficiency of ΔBcHPt + BcHptK161Q-GFP was similar to that of B05.10 (Figure 8). The decreased virulence of ΔBcHPt + BcHptK161Q-GFP was likely due to its minor growth reduction, and increased sensitivity to osmotic and oxidative stresses.
Figure 7. Constitutive acetylation of BcHpt Lys161 affects the virulence of B. cinerea. Virulence tests were performed on tomato leaves following inoculation with B05.10, ΔBcHpt + BcHptK161Q-GFP, ΔBcHpt + BcHptK161R-GFP, and ΔBcHpt + BcHptK161K-GFP or agar plugs as negative controls (CK). (A) Disease symptoms caused by B05.10, ΔBcHpt + BcHptK161Q-GFP, ΔBcHpt + BcHptK161R-GFP, and ΔBcHpt + BcHptK161K-GFP on wounded tomato leaves 3 days after inoculation. (B) Diameter of disease lesions on tomato leaves after 3 days of incubation. The bars denote the standard errors of four replicates. Statistical tests were carried out using Tukey’s test for multiple comparisons, and values on the bars followed by different letters are significantly different at P = 0.05.
Figure 8. Onion epidermis penetration assay with B05.10 and ΔBcHpt + BcHptK161Q-GFP. Conidia of B05.10 and ΔBcHpt + BcHptK161Q-GFP were incubated with onion epidermis for 15 h before the cultures were photographed.
Discussion
Fungal HOG pathways play various roles in cells, including in osmotic, oxidative and other (e.g., cell wall) stress responses; fungicide (phenylpyrrole and dicarboximide) sensitivity; hyphal growth and development; virulence; sexual and asexual development; secondary metabolite production; and dimorphic phase transitions (Li et al., 2010; Jiang et al., 2018). HPt is essential in S. cerevisiae, N. crassa, A. nidulans, and B. cinerea (Maeda et al., 1994; Banno et al., 2007; Furukawa et al., 2010), but the function of BcHPt in B. cinerea remains unclear. However, we identified one acetylation site in BcHpt in our previous proteomics studies (Lv et al., 2016). To determine the role of lysine acetylation in BcHpt, we characterized Lys161 of BcHpt in B. cinerea by site directed mutagenesis. Since deletion of BcHPt is lethal, but ΔBcHPt + BcHptK161Q-GFP could be produced, acetylation at residue 161 is relevant for some but not all activities of the HPt protein. The regulatory mechanisms of this residue in the HPt protein remain to be further discovered.
To further characterize the function of BcHpt Lys161, we analyzed the phosphorylation levels of BcSak1 in the mutants. Under stresses, the phosphorylation levels of BcSak1 in the ΔBcHPt + BcHptK161Q-GFP mutant were significantly lower than that in the wild-type strain (Figure 5), indicating that acetylation of BcHpt Lys161 affects the phosphorylation level of BcSak1. The HPt domain of S. cerevisiae Ypd1 contains a four-helix bundle (αB-αC-αD-αG) as a core structure. In addition, His64 is predicted to play important roles in the phosphoryl transfer activity of Ypd1. The involved residues are conserved in many organisms and have been functionally characterized (Janiak and West, 2000; Fassler and West, 2013). In B. cinerea, this residue is His120, and the acetylation site is Lys161, which is not inside the HPt domain of BcHpt (66–146) (Figure 1).
ΔBcHPt + BcHptK161Q-GFP exhibited increased sensitivity to multiple stresses, including fungicides (dicarboximides and DMIs), oxidative, osmotic and cell wall stresses (Figure 3). However, its phenotype under oxidative stress is minor and very different from that observed for a BcSak1 deletion strain (Segmüller et al., 2007). It is reasonable since deletion of BcSak1 led to a complete loss of phosphorylated BcSak1; however, ΔBcHPt + BcHptK161Q-GFP, still has a relatively small amount of phosphorylated BcSak1. In addition, BcSak1 deletion mutants, but not ΔBcHPt + BcHptK161Q-GFP, showed conidium and sclerotium formation changes compared to wild type strain (Segmüller et al., 2007). This is likely because, unlike BcSak1 deletion mutants, ΔBcHPt + BcHptK161Q-GFP still has a low level of phosphorylated BcSak1.
ΔBcHPt + BcHptK161Q-GFP also showed increased sensitivity to triadimefon, which is an ergosterol biosynthesis inhibitor. This phenotype was consistent with those of BRrg1 and BcSkn7 (Yan et al., 2011; Yang et al., 2015). Considering phosphorylation levels of BcSak1 in three mutants (ΔBcHPt + BcHptK161Q-GFP, ΔBRrg1, and ΔBcSkn7) were lower than that in the wild-type strain (Yan et al., 2011; Yang et al., 2015), the triadimefon sensitivity was most likely a consequence of lower phosphorylation levels of BcSak1. However, the mechanisms of HOG pathway involvement in the response to an impaired ergosterol biosynthesis in B. cinerea have not been addressed so far.
Acetylation of BcHpt cannot be detected in response to osmotic (NaCl), oxidative (H2O2), and cell wall (Congo red) stresses, while the phosphorylation levels of BcSak1 were upregulated. These results indicate that acetylation could affect the phosphoryl transfer activity of BcHpt. Since acetylation and phosphorylation are both common post-translational protein modifications, further exploring the regulatory mechanism mediating the acetylation and phosphorylation of BcHpt would be of interest. Moreover, as a cell wall stressor, Congo red could also activate the HOG pathway (Figure 5) and cause deacetylation of BcHpt (Figure 6), which provides further evidence of the cross talk between the HOG and the CWI pathway.
Saccharomyces cerevisiae does not contain such a Lys residue, and only 4 out of 11 fungal species analyzed (B. cinerea, Sclerotinia sclerotiorum, F. graminearum, and P. marneffei) possess the Lys residue. These four fungal species belonged to Pezizomycotina, however, other Pezizomycotina fungi (N. crassa and Magnaporthe oryzae) do not have this residue. Acetylation of this site may impact the phosphorylation of BcSak1 in response to the above cited stresses, which means that the signal transduction is interrupted. Thus, our hypothesis is that acetylation/deacetylation of BcHpt plays important roles in its signal transduction function in the four cited species. As this residue is not conserved among fungal species, further studies are required to investigate signal transduction in the remaining species.
Conclusion
In conclusion, the acetylation of BcHpt Lys161 plays a significant role in the hyphal growth, osmotic and oxidative stress responses, and dicarboximide and DMI sensitivities of B. cinerea. BcHpt is deacetylation in response to osmotic, oxidative, and cell wall stresses. Constitutive acetylation of Lys161 interrupted signal transduction of the HOG pathway, which led to lower phosphorylation levels of BcSak1. The regulatory mechanism mediating these two common post-translational protein modifications requires further study.
Data Availability Statement
The raw data supporting the conclusions of this artilcle will be made available by the authors, without undue reservation, to any qualified researcher.
Author Contributions
QY and WL generated the hypothesis and planned the experiments. QY, LS, ZM, and MS performed the experiments. QY, WL, and YH wrote the manuscript. All other authors provided comments on the manuscript.
Funding
This research was supported by the National Science Foundation (31601624, 31701746, and 31722044), the Shandong Provincial Natural Science Foundation (ZR2016CQ03), the Taishan Scholar Construction Foundation of Shandong Province (tshw20130963), the Key Research and Development Program of Shandong Province (2017NC210004, 2018GNC111003, and 2018GNC113017), the Qingdao People’s Livelihood Technology Project (17-3-3-55-nsh), a project of the Shandong Province Higher Educational Science and Technology Program (J18KA118), Open Funding Project of State Key Laboratory of Microbial Metabolism (MMLKF17-05), First Class Grassland Science Discipline Program in Shandong Province, China, and the Scientific Research Fund for High-level Talents in Qingdao Agricultural University (No. 1116025).
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.2019.02965/full#supplementary-material
FIGURE S1 | Generation strategy of the BcHpt mutant strains of B. cinerea. The hygromycin-resistant cassette (HPH) is denoted by the large gray box, GFP is denoted by the large white box and the neomycin-resistant cassette (NEO) is denoted by the large shadowed area. The vectors: BcHptK161K-GFP- pYF11, BcHptK161R-GFP- pYF11, and BcHptK161Q-GFP- pYF11 (BcHptK161K/R/Q-GFP-pYF11), which contain the site-directed mutagenesis were transformed into the B05.10 strain. The resulting transformants were named B05.10 + BcHptK161Q-GFP, B05.10 + BcHptK161R-GFP, and B05.10 + BcHptK161K-GFP (B05.10 + BcHptK161K/R/Q-GFP), respectively. Subsequently, the native BcHpt locus in the resulting transformants was deleted by a homologous recombination strategy to generate the mutant ΔBcHpt + BcHptK161Q, ΔBcHpt + BcHptK161R, and ΔBcHpt + BcHptK161K (ΔBcHpt + BcHptK161K/R/Q) strains. Primer (codes 1–17) binding sites are indicated by arrows.
FIGURE S2 | Relative expression levels of BcHpt in B05.10, ΔBcHpt + BcHptK161Q-GFP, ΔBcHpt + BcHptK161R-GFP, and ΔBcHpt + BcHptK161K-GFP. Mycelia were incubated in PDB for 2 days. The bars denote the standard errors from three experiments, and statistical tests were carried out using Tukey’s test for multiple comparisons. Values on the bars followed by different letters are significantly different at P = 0.05.
FIGURE S3 | Morphology of B05.10 and two additional K161Q mutants: ΔBcHpt + BcHptK161Q-GFP-2, and ΔBcHpt + BcHptK161Q-GFP-3 on PDA medium after 3 days of incubation.
FIGURE S4 | Sensitivity of B05.10, ΔBcHpt + BcHptK161Q-GFP, ΔBcHpt + BcHptK161R-GFP, and ΔBcHpt + BcHptK161K-GFP grown on PDA containing different compounds at the concentrations noted in the figure. The cultures were photographed after 48 h of incubation at 25°C.
FIGURE S5 | Comparison of BcSak1 phosphorylation in B05.10 and ΔBcHpt + BcHptK161Q-GFP in response to 1 M NaCl and 1 μg/ml iprodione.
TABLE S1 | Primers used in the study.
References
Aoyama, K., Mitsubayashi, Y. H., and Mizuno, T. (2000). Spy1, a histidine-containing phosphotransfer signaling protein, regulates the fission yeast cell cycle through the Mcs4 response regulator. J. Bacteriol. 182, 4868–4874. doi: 10.1128/jb.182.17.4868-4874.2000
Banno, S., Noguchi, R., Yamashita, K., Fukumori, F., Kimura, M., Yamaguchi, I., et al. (2007). Roles of putative His-to-Asp signaling modules HPT-1 and RRG-2, on viability and sensitivity to osmotic and oxidative stresses in Neurospora crassa. Curr. Genet. 51, 197–208. doi: 10.1007/s00294-006-0116-8
Barrett, J. F., and Hoch, J. A. (1998). Two-component signal transduction as a target for microbial anti-infective therapy. Antimicrob. Agents Chemother. 42, 1529–1536. doi: 10.1128/aac.42.7.1529
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
Dean, R., Van Kan, J. A., Pretorius, Z. A., Hammond-Kosack, K. E., Di Pietro, A., Spanu, P. D., et al. (2012). The top 10 fungal pathogens in molecular plant pathology. Mol. Plant Pathol. 13, 414–430. doi: 10.1111/j.1364-3703.2011.00783.x
Dixon, K. P., Xu, J. R., Smirnoff, N., and Talbot, N. J. (1999). Independent signaling pathways regulate cellular turgor during hyperosmotic stress and appressorium-mediated plant infection by Magnaporthe grisea. Plant Cell 11, 2045–2058. doi: 10.1105/tpc.11.10.2045
Doehlemann, G., Berndt, P., and Hahn, M. (2006). Different signalling pathways involving a Gα protein, cAMP and a MAP kinase control germination of Botrytis cinerea conidia. Mol. Microbiol. 59, 821–835. doi: 10.1111/j.1365-2958.2005.04991.x
Fassler, J. S., and West, A. H. (2013). Histidine phosphotransfer proteins in fungal two-component signal transduction pathways. Eukaryot. Cell 12, 1052–1060. doi: 10.1128/EC.00083-13
Fillinger, S., Ajouz, S., Nicot, P., and Pierre, L. (2012). Functional and structural comparison of pyrrolnitrin- and iprodione-induced modifications in the class III histidine-kinase Bos1 of Botrytis cinerea. PLoS One 7:e42520. doi: 10.1371/journal.pone.0042520
Fillinger, S., and Elad, Y. (2016). Botrytis – the Fungus, the Pathogen and its Management in Agricultural Systems. New York, NY: Springer International Publishing, 413–486.
Fujimura, M., Ochiai, N. M., Motoyama, T., Ichiishi, A., and Usami, R. (2003). Putative homologs of SSK22 MAPKK kinase and PBS2 MAPK kinase of Saccharomyces cerevisiae encoded by os-4 and os-5 genes for osmotic sensitivity and fungicide resistance in Neurospora crassa. J. Agric. Chem. Soc. Japan 67, 186–191.
Furukawa, K., Hoshi, Y., Maeda, T., Nakajima, T., and Abe, K. (2010). Aspergillus nidulans HOG pathway is activated only by two-component signalling pathway in response to osmotic stress. Mol. Microbiol. 56, 1246–1261. doi: 10.1111/j.1365-2958.2005.04605.x
Gronover, C. S., Kasulke, D., Tudzynski, P., and Tudzynski, B. (2001). The role of G protein alpha subunits in the infection process of the gray mold fungus Botrytis cinerea. Mol. Plant Microbe. Interact 14, 1293–1302. doi: 10.1094/mpmi.2001.14.11.1293
Gu, Q., Zhang, C., Liu, X., and Ma, Z. (2015). A transcription factor FgSte12 is required for pathogenicity in Fusarium graminearum. Mol. Plant Pathol. 16, 1–13. doi: 10.1111/mpp.12155
Hohmann, S. (2002). Osmotic stress signaling and osmoadaptation in yeasts. Microbiol. Mol. Biol. R. 66, 300–372. doi: 10.1128/mmbr.66.2.300-372.2002
Jacob, S., Foster, A. J., Yemelin, A., and Thines, E. (2015). High osmolarity glycerol (HOG) signalling in Magnaporthe oryzae: identification of MoYPD1 and its role in osmoregulation, fungicide action, and pathogenicity. Fungal Biol. 119, 580–594. doi: 10.1016/j.funbio.2015.03.003
Janiak, S. F., and West, A. (2000). Functional roles of conserved amino acid residues surrounding the phosphorylatable histidine of the yeast phosphorelay protein YPD1. Mol. Microbiol. 37, 136–144. doi: 10.1046/j.1365-2958.2000.01973.x
Jespers, A. B. K., Davidse, L. C., and Dewaard, M. A. (1993). Biochemical effects of the phenylpyrrole fungicide fenpiclonil in Fusarium sulphureum (Schlecht). Pestic. Biochem. Phys. 45, 116–129. doi: 10.1006/pest.1993.1014
Jespers, A. B. K., and Waard, M. A. D. (2010). Effect of fenpiclonil on phosphorylation of glucose in Fusarium sulphureum. Pest Manag. Sci. 44, 167–175. doi: 10.1002/ps.2780440210
Jiang, C., Zhang, X., Liu, H., and Xu, J. R. (2018). Mitogen-activated protein kinase signaling in plant pathogenic fungi. PLoS Phathog. 14:e1006875. doi: 10.1371/journal.ppat.1006875
Jiang, J., Liu, X., Yin, Y., and Ma, Z. (2011). Involvement of a velvet protein FgVeA in the regulation of asexual development, lipid and secondary metabolisms and virulence in Fusarium graminearum. PLoS One 6:e28291. doi: 10.1371/journal.pone.0028291
Jones, C. A., Greer-Phillips, S. E., and Borkovich, K. A. (2007). The response regulator RRG-1 functions upstream of a mitogen-activated protein kinase pathway impacting asexual development, female fertility, osmotic stress, and fungicide resistance in Neurospora crassa. Mol. Biol. Cell 18, 2123–2136. doi: 10.1091/mbc.e06-03-0226
Koretke, K. K., Lupas, A. N., Warren, P. V., Rosenberg, M., and Brown, J. R. (2000). Evolution of two-component signal transduction. Mol. Biol. Evol. 17, 1956–1970.
Krantz, M., Becit, E., and Hohmann, S. (2006). Comparative analysis of HOG pathway proteins to generate hypotheses for functional analysis. Curr. Genet. 49, 152–165. doi: 10.1007/s00294-005-0039-9
Kruppa, M., and Calderone, R. (2006). Two-component signal transduction in human fungal pathogens. FEMS Yeast Res. 6, 149–159. doi: 10.1111/j.1567-1364.2006.00024.x
Lee, J. W., Ko, Y. J., Kim, S. Y., and Bahn, Y. S. (2011). Multiple roles of Ypd1 phosphotransfer protein in viability, stress response, and virulence factor regulation in Cryptococcus neoformans. Eukaryot. Cell 10, 998–1002. doi: 10.1128/EC.05124-11
Li, D., Agrellos, O. A., and Calderone, R. (2010). Histidine kinases keep fungi safe and vigorous. Curr. Opin. Microbiol. 13, 424–430. doi: 10.1016/j.mib.2010.04.007
Li, X., Zhang, S., Blander, G., Tse, J. G., Krieger, M., and Guarente, L. (2007). SIRT1 deacetylates and positively regulates the nuclear receptor LXR. Mol. Cell 28, 91–106. doi: 10.1016/j.molcel.2007.07.032
Liu, W., Leroux, P., and Fillinger, S. (2008). Th The HOG1-like MAP kinase Sak1 of Botrytis cinerea is negatively regulated by the upstream histidine kinase Bos1 and is not involved in dicarboximide- and phenylpyrrole-resistance. Fungal Genet. Biol. 45, 1062–1074. doi: 10.1016/j.fgb.2008.04.003
Livak, K. J., and Schmittgen, T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCT method. Methods 25, 402–408. doi: 10.1006/meth.2001.1262
Lv, B., Yang, Q., Li, D., Liang, W., and Song, L. (2016). Proteome-wide analysis of lysine acetylation in the plant pathogen Botrytis cinerea. Sci. Rep. 6, 29313. doi: 10.1038/srep29313
Maeda, T., Takekawa, M., and Saito, H. (1995). Activation of yeast PBS2 MAPKK by MAPKKKs or by binding of an SH3-containing osmosensor. Science 269, 554–558. doi: 10.1126/science.7624781
Maeda, T., Wurglermurphy, S. M., and Saito, H. (1994). A two-component system that regulates an osmosensing MAP kinase cascade in yeast. Nature 369, 242–245. doi: 10.1038/369242a0
Mavrianos, J., Desai, C., and Chauhan, N. (2014). Two-component histidine phosphotransfer protein Ypd1 is not essential for viability in Candida albicans. Eukaryot. Cell 13, 452–460. doi: 10.1128/ec.00243-13
McDonald, B., and Martinez, J. (1990). Restriction fragment length polymorphisms in Septoria tritici occur at a high frequency. Curr. Genet. 17, 133–138. doi: 10.1007/bf00312858
Mehrabi, R., Zwiers, L. H., De Waard, M. A., and Kema, G. H. (2006). MgHog1 regulates dimorphism and pathogenicity in the fungal wheat pathogen Mycosphaerella graminicola. Mol. Plant Microbe. Interact 19, 1262–1269. doi: 10.1094/mpmi-19-1262
Moriwaki, A., Kubo, E., Arase, S., and Kihara, J. (2006). Disruption of SRM1, a mitogen-activated protein kinase gene, affects sensitivity to osmotic and ultraviolet stressors in the phytopathogenic fungus Bipolaris oryzae. FEMS Microbiol. Lett. 257, 253–261. doi: 10.1111/j.1574-6968.2006.00178.x
Motoyama, T., Kadokura, K., Ohira, T., Ichiishi, A., Fujimura, M., Yamaguchi, I., et al. (2005). A two-component histidine kinase of the rice blast fungus is involved in osmotic stress response and fungicide action. Fungal Genet. Biol. 42, 200–212. doi: 10.1016/j.fgb.2004.11.002
Motoyama, T., Ochiai, N., Morita, M., Iida, Y., Usami, R., and Kudo, T. (2008). Involvement of putative response regulator genes of the rice blast fungus Magnaporthe oryzae in osmotic stress response, fungicide action, and pathogenicity. Curr. Genet. 54, 185–195. doi: 10.1007/s00294-008-0211-0
Noguchi, R., Banno, S., Ichikawa, R., Fukumori, F., Ichiishi, A., Kimura, M., et al. (2007). Identification of OS-2 MAP kinase-dependent genes induced in response to osmotic stress, antifungal agent fludioxonil, and heat shock in Neurospora crassa. Fungal Genet. Biol. 44, 208–218. doi: 10.1016/j.fgb.2006.08.003
Noriyuki, O., Takeshi, T., Takumi, N., Naoko, T. A., Makoto, F., and Makoto, K. (2007). Involvement of the osmosensor histidine kinase and osmotic stress-activated protein kinases in the regulation of secondary metabolism in Fusarium graminearum. Biochem. Bioph. Res. Commun. 363, 639–644. doi: 10.1016/j.bbrc.2007.09.027
O’Rourke, S. M., and Herskowitz, I. (1998). The Hog1 MAPK prevents cross talk between the HOG and pheromone response MAPK pathways in Saccharomyces cerevisiae. Genes Dev. 12, 2874–2886. doi: 10.1101/gad.12.18.2874
Posas, F., and Saito, H. (1997). Osmotic activation of the HOG MAPK pathway via Ste11p MAPKKK: scaffold role of Pbs2p MAPKK. Science 276, 1702–1705. doi: 10.1126/science.276.5319.1702
Posas, F., Wurgler-Murphy, S. M., Maeda, T., Witten, E. A., Thai, T. C., and Saito, H. (1996). Yeast HOG1 MAP kinase cascade is regulated by a multistep phosphorelay mechanism in the SLN1-YPD1-SSK1” two-component” osmosensor. Cell 86, 865–875. doi: 10.1016/s0092-8674(00)80162-2
Quidde, T., Büttner, P., and Tudzynski, P. (1999). Evidence for three different specific saponin-detoxifying activities in Botrytis cinerea and cloning and functional analysis of a gene coding for a putative avenacinase. Eur. J. Plant Pathol. 105, 273–283.
Santos, J. L., and Shiozaki, K. (2001). Fungal histidine kinases. Sci. Signal. 2001, re1. doi: 10.1126/scisignal.982001re1
Schwer, B., Bunkenborg, J., Verdin, R. O., Andersen, J. S., and Verdin, E. (2006). Reversible lysine acetylation controls the activity of the mitochondrial enzyme Acetyl-CoA synthetase 2. Proc. Natl. Acad. Sci. U.S.A. 103, 10224–10229. doi: 10.1073/pnas.0603968103
Segmüller, N., Ellendorf, U., Tudzynski, B., and Tudzynski, P. (2007). BcSAK1, a stress-activated mitogen-activated protein kinase, is involved in vegetative differentiation and pathogenicity in Botrytis cinerea. Eukaryot. Cell 6, 211–221. doi: 10.1128/ec.00153-06
Vargas-Pérez, I., Sánchez, O., Kawasaki, L., Georgellis, D., and Aguirre, J. (2007). Response regulators SrrA and SskA are central components of a phosphorelay system involved in stress signal transduction and asexual sporulation in Aspergillus nidulans. Eukaryot. Cell 6:1570. doi: 10.1128/ec.00085-07
Viaud, M., Fillinger, S., Liu, W., Polepalli, J. S., Le, P. P., Kunduru, A. R., et al. (2006). A class III histidine kinase acts as a novel virulence factor in Botrytis cinerea. Mol. Plant Microbe. Interact 19, 1042–1050. doi: 10.1094/mpmi-19-1042
Vitalini, M. W., Paula, R. M. D., Goldsmith, C. S., Jones, C. A., Borkovich, K. A., and Bell-Pedersen, D. (2007). Circadian rhythmicity mediated by temporal regulation of the activity of p38 MAPK. Proc. Natl. Acad. Sci. U.S.A. 104, 18223–18228. doi: 10.1073/pnas.0704900104
Williamson, B., Tudzynski, B., Tudzynski, P., and Van Kan, J. A. (2007). Botrytis cinerea: the cause of grey mould disease. Mol. Plant Pathol. 8, 561–580. doi: 10.1111/j.1364-3703.2007.00417.x
Yan, L., Yang, Q., Jiang, J., Michailides, T. J., and Ma, Z. (2011). Involvement of a putative response regulator Brrg-1 in the regulation of sporulation, sensitivity to fungicides, and osmotic stress in Botrytis cinerea. App. Microbiol. Biot. 90, 215–226. doi: 10.1007/s00253-010-3027-z
Yan, L., Yang, Q., Sundin, G. W., Li, H., and Ma, Z. (2010). The mitogen-activated protein kinase kinase BOS5 is involved in regulating vegetative differentiation and virulence in Botrytis cinerea. Fungal Genet. Biol. 47, 753–760. doi: 10.1016/j.fgb.2010.06.002
Yang, Q., Yan, L., Qin, G., and Ma, Z. (2012). The mitogen-activated protein kinase kinase kinase BcOs4 is required for vegetative differentiation and pathogenicity in Botrytis cinerea. App. Microbiol. Biot. 96, 481–492. doi: 10.1007/s00253-012-4029-9
Yang, Q., Yin, D., Yin, Y., Cao, Y., and Ma, Z. (2015). The response regulator BcSkn7 is required for vegetative differentiation and adaptation to oxidative and osmotic stresses in Botrytis cinerea. Mol. Plant Pathol. 16, 276–287. doi: 10.1111/mpp.12181
Yang, Q., Yu, F., Yin, Y., and Ma, Z. (2013). Involvement of protein tyrosine phosphatases BcPtpA and BcPtpB in regulation of vegetative development, virulence and multi-stress tolerance in Botrytis cinerea. PLoS One 8:e61307. doi: 10.1371/journal.pone.0061307
Yang, Q., Zhang, J., Hu, J., Wang, X., Lv, B., and Liang, W. (2018). Involvement of BcYak1 in the regulation of vegetative differentiation and adaptation to oxidative stress of Botrytis cinerea. Front. Microbiol. 9:281. doi: 10.3389/fmicb.2018.00281
Yu, J. H., Hamari, Z., Han, K. H., Seo, J. A., Reyes-Domínguez, 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
Keywords: Botrytis cinerea, lysine acetylation, histidine phosphotransfer protein, dicarboximides and demethylation inhibitor sensitivity, osmotic and oxidative stress, virulence
Citation: Yang Q, Song L, Miao Z, Su M, Liang W and He Y (2020) Acetylation of BcHpt Lysine 161 Regulates Botrytis cinerea Sensitivity to Fungicides, Multistress Adaptation and Virulence. Front. Microbiol. 10:2965. doi: 10.3389/fmicb.2019.02965
Received: 02 April 2019; Accepted: 09 December 2019;
Published: 08 January 2020
Edited by:
Daisuke Hagiwara, University of Tsukuba, JapanReviewed by:
Frank Ebel, Ludwig Maximilian University of Munich, GermanySabine Fillinger, Institut National de la Recherche Agronomique (INRA), France
Copyright © 2020 Yang, Song, Miao, Su, Liang and He. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Wenxing Liang, d2xpYW5nMUBxYXUuZWR1LmNu
†These authors have contributed equally to this work and share first authorship