1. Introduction
Sclerotinia sclerotiorum is a widespread fungal pathogen that parasitizes various hosts, causing severe diseases in over 600 plant species globally [
1,
2], which have significantly impact on the yield and quality of important economic crops such as canola and soybeans [
3]. To successfully infect and parasitize its hosts,
S. sclerotiorum has evolved complex and sophisticated infection strategies. In the early stages of infection, the fungus overcomes host defense responses through various pathways, establishing a short biotrophic phase, including the secretion of effectors and the formation of infection cushions [
2]. Notably,
SsCmu1 is significantly upregulated during infection of
Brassica napus [
4], which encodes a secreted cutinase, capable of inhibiting host salicylic acid synthesis, thereby promoting infection [
5]. Other putative effectors, such as
SsITL [
6],
SsCM1 [
7],
SsCVNH [
8],
Ssv263 [
9], and fungal effector proteins containing LysM domains, also exhibit high upregulation in the early infection stages, impacting the virulence of
S. sclerotiorum [
10]. After a brief biotrophic phase, the fungus induces host cell necrosis and cell wall degradation by producing numerous toxins, necrosis-inducing secreted proteins, and cell wall-degrading enzymes [
2].
For pathogens, in order to adapt to changing host environments and defenses, infection processes require rapid and subtle adjustments to their gene expression programs. [
11,
12,
13]. The regulation of gene expression is achieved at the transcriptional level through various mechanisms, include histone modification, as one of the primary epigenetic regulatory mechanisms, playing a particularly important role in shaping fungal pathogenicity [
14]. Histones can undergo various covalent modifications, including methylation, acetylation, phosphorylation, and ubiquitination [
15]. Histone methylation primarily occurs on the side chains of lysine and arginine. Lysine can undergo monomethylation, dimethylation, or trimethylation, while arginine may undergo monomethylation, symmetric dimethylation, or asymmetric demethylation [
15,
16]. Histone lysine methylation, a process orchestrated by histone lysine methyltransferases (HKMTs) with the SET domain, involves the transfer of a methyl group from S-adenosyl-L-methionine (SAM) to lysine residues at the N-terminus of H3 or H4 histones [
17]. Typically, lysine 9 of histone H3 (H3K9) is linked to transcriptionally inactive heterochromatin, serving as a highly conserved epigenetic marker for heterochromatin formation and transcriptional silencing [
18]. As a result, histone modifications play a crucial role in the regulation of gene expression, especially during the infection process. Through these modification mechanisms, fungi can flexibly respond to dynamic changes in the host environment.
H3K9 methylation is facilitated by specific proteins, including Clr4 (Cryptic loci regulator 4) in
Schizosaccharomyces pombe and Dim5 (Defective in methylation 5) in
Neurospora crassa [
19]. In
Botrytis cinerea, a plant pathogenic fungus, the absence of the
Dim5 gene results in a significant reduction of H3K9me3, causing the downregulation of pathogenic genes related to host signal perception, host tissue colonization, stress response, toxin synthesis, and host immune response [
20]. In the maize pathogen
Fusarium verticillioides, disruption of
Dim5 significantly reduces H3K9me3 levels, leading to a pronounced decrease in fungal virulence, accompanied by an unexpected increase in osmotic stress tolerance and expression of melanin synthesis genes [
21]. In
Fusarium mangiferae, the absence of the
KMT1 gene significantly impedes the biosynthesis of fumonisin and deoxynivalenol toxins in mango [
22]. These indicate that the maintenance of H3K9me3 is crucial for the virulence of pathogenic fungi. However, currently there is no report on the biological function of H3K9me3 in
S. sclerotiorum.
In this study, we characterized the histone H3 lysine methyltransferase SsDim5 by exploring its roles in H3 lysine trimethylation, regulation of mycotoxins synthesis, and pathogenicity as well as response to external stress, aimed to gain deeper insights into the physiological function of epigenetic regulation in the fungal pathogen S. sclerotiorum.
2. Materials and Methods
2.1. Fungal Strains, Plants, and Culture Conditions
The wild-type (WT) strain 1980 of
S. sclerotiorum [
23] was cultured on Potato Dextrose Agar (PDA), while the knockout mutants
Ssdim5-1/Ssdim5-2 were grown on PDA supplemented with 150 μg/mL hygromycin (Roche). The genetic complementation strain
SsDim5-C was cultured on PDA containing 100 μg/mL Geneticin (G418). All these strains were incubated at 20 °C.
Nicotiana benthamiana and B. napus (ZS11) plants used for pathogenicity tests were obtained from our laboratory stock and cultivated at 22 °C with a 16-h light/8-h dark photoperiod.
2.2. Bioinformatics Analysis of SsDim5
First, the genomic sequence of
BcDim5 from
B. cinerea was obtained by accessing the NCBI database (
http://www.ncbi.nlm.nih.gov/). BlastP analysis was then conducted to identify its orthologs in
S. sclerotiorum and other species. Subsequently, multiple sequence alignment was performed using DNAMAN 6.0 (Lynnon BioSoft, Quebec, Canada). Finally, a phylogenetic tree was constructed using the Maximum Likelihood method with MEGA 6.0 software [
24].
2.3. Gene Knockout and Genetic Complementation of SsDim5
Refer to the previous method [
25], sequences of the
SsDim5 gene were amplified using genomic DNA from WT strain 1980 as a template. Subsequently, these sequences were fused with the left and right portions of the hygromycin expression cassette, generating the gene knockout fragment. The knockout fragment was then introduced into WT protoplasts using PEG-mediated protoplast transformation [
26]. Selection was performed by screening on PDA medium containing 150 mg/L hygromycin, resulting in the putative knockout transformants of
SsDim5. Finally, through successive sub-culturing of mycelial tips, pure knockout strains were obtained and confirmed by PCR and qRT-PCR.
Simultaneously, using genomic DNA from the WT as a template, the sequence containing the full-length SsDim5 gene, including its native promoter, was amplified. After digestion with KpnⅠ and EcoRⅠ enzymes, this fragment was ligated into the linearized pCH-NEO1 vector. The vector containing the SsDim5 gene fragment was introduced into the knockout mutant Ssdim5-1 through Agrobacterium tumefaciens-mediated transformation, thereby obtaining a genetically complemented strain.
2.4. DNA Extraction, RNA Extraction and cDNA Synthesis
Fresh mycelia were inoculated onto PDA plates covered with glass paper. Following a 2-day incubation, the mycelia were harvested, rapidly frozen in liquid nitrogen, and subsequently pulverized into powder. Genomic DNA extraction was performed using the cetyltrimethylammonium bromide (CTAB) method, as described by Allen et al. [
27].
A commercial RNA extraction kit (AG21019, Accurate Biotechnology (Hunan) Co., Ltd., Changsha, China) was utilized for RNA extraction. Following the manufacturer’s guidelines, first-strand cDNA synthesis was conducted using the Evo M-MLV reverse transcription kit (AG11705, Accurate Biotechnology (Hunan) Co., Ltd., Changsha, China).
2.5. Quantitative Real-Time PCR (qRT-PCR) Analysis
The qRT-PCR experiments were conducted using the StepOne™ Real-Time PCR System and the SYBR
® Green Premix Pro Taq HS qPCR Kit II (AG11702, Accurate Biotechnology (Hunan) Co., Ltd., Changsha, China). The PCR program consisted of 40 cycles, including an initial denaturation at 94 °C for 2 min, denaturation at 94 °C for 15 s, and annealing at 58 °C for 1 min.
SsTubulin1 was employed as the reference gene. The analysis of relative gene expression levels utilized the 2^
(-ΔΔCT) method [
28].
2.6. Inoculation and Virulence Determination
Inoculation was performed following previously established protocols [
29]. Mycelial plugs obtained from the actively growing colony edge were used for inoculating
N. benthamiana and
B. napus leaves (diameter 6 mm). Following inoculation, the leaves were cultured at 22 °C with a relative humidity of 95–100%. Photographic documentation was carried out 24 h post-infection.
2.7. Appressorium Observation and Oxalic Acid Analysis
Mycelial plugs (diameter 6 mm) were obtained from the actively growing edges. The mycelia were embedded in agar blocks and placed on glass slides. After 16 h of incubation, the morphology and quantity of adherent cells were observed. Additionally, Mycelial plugs (diameter 1 mm) were inoculated onto the epidermis of onion epidermis. After 16 h of invasion, the onion epidermis was stained in a 0.5% trypan blue solution for 30 min. Subsequently, bleaching solution (ethanol: acetic acid: glycerol = 3:1:1) was used for decolorization. Samples were observed and photographed using an optical microscope (Axio Imager 2, ZEISS, Oberkochen, Germany).
Furthermore, mycelial plugs (diameter 6 mm) were inoculated on PDA medium containing 100 μg/mL bromophenol blue to assess oxalic acid secretion.
2.8. Western Blot Analysis of H3K9 Trimethylates
The mycelia were collected from PDA plates covered with glass paper membrane, approximately 500 mg in weight. After rapid freezing in liquid nitrogen and grinding, the mycelial powder was added to 100 μL of protein extraction buffer (200 mM pH = 6.8 Tris-HCl, 40% glycerol, 20% β-mercaptoethanol, 8% SDS, 0.4% bromophenol blue) for total protein extraction. The obtained proteins were then separated on a 10% denaturing polyacrylamide gel using sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) and subsequently transferred to a nitrocellulose membrane using a Bio-Rad blotting apparatus.
After incubation with the primary antibody, incubation with the secondary antibody followed, and protein detection was carried out using the SuperSignal West Pico PLUS (product number: 34095, Thermo Scientific). The antibodies used, along with their sources and dilutions, were as follows: monoclonal antibody against histone H3 (Immunoway, 1:2000), monoclonal antibody against histone H3K9 (Abcam, 1:2000), and HRP-conjugated Rabbit anti-mouse IgG (1:5000).
2.9. Abiotic Stress Response Assay
To assess the response of the Ssdim5-1/Ssdim-2 mutants to various stresses, we cultured them on PDA media supplemented with 1M NaCl, 1M KCl, 1M sorbitol, 0.005% SDS and 10 mM H2O2, respectively. After incubation for 48 h, we measured the mycelial diameter and calculated the growth inhibition rate using the formula: Growth Inhibition Rate (%) = 100 × (Colony diameter on pure PDA—Colony diameter under different stress conditions) / (Colony diameter on pure PDA).
2.10. RNA Sequencing and Data Analysis
Cultures of WT and Ssdim5-1 mutant mycelia, grown in PDB medium for 2 days, were harvested for transcriptome sequencing at Biomarker Technologies Corporation. Sequencing was performed on the Illumina NovaSeq platform according to the manufacturer’s instructions, generating 150 bp paired-end sequences. During data cleaning, sequences containing adapters, poly-N sequences, and low-quality sequences were excluded to obtain high-quality clean data. Parameters such as Q20, Q30, GC content, and sequence duplication levels were computed. Following data validation, alignment to the reference genome sequence was executed. Gene functions were annotated through sequence alignment using various databases: Nr (NCBI non-redundant protein sequences), Pfam (protein families), KOG/COG (clusters of orthologous groups for eukaryotic complete genomes), Swiss-Prot (manually annotated and reviewed protein sequence database), KO (KEGG Ortholog database), and GO (gene ontology). Original counts were normalized using fragments per kilobase of transcript per million mapped reads (FPKM), and differential expression analysis was carried out with DESeq2. Genes meeting the criteria of a corrected p-value < 0.01 and a fold change ≥ 2, as determined by DESeq2 analysis, were considered differentially expressed. GO enrichment analysis of differentially expressed genes utilized the clusterProfiler package, employing the Wallenius non-central hypergeometric distribution.
4. Discussion
In eukaryotes, lysine methylation modification of histones is mediated by histone lysine methyltransferases (KMT) [
31], and H3K9 methylation relies on the activity of the KMT1 family [
32]. Here, we identified a putative H3K9 methyltransferase, named
SsDim5, by performing blast analysis on the protein sequence of the KMT1 family member
BcDim5 in
B. cinerea.
SsDim5 is predicted to have a SET (Su(var)3–9, Enhancer-of-zeste, and Trithorax) domain, commonly associated with histone N-terminal lysine methyltransferase activity [
33]. Previous studies have shown that human SUV39H1 and mouse Suv39h1 possess typical SET domains and catalyze directed H3K9 methylation modifications [
34]. The Dim5 in
N. crassa can specifically catalyze trimethylation of H3K9 [
19]. In
S. pombe, the histone methyltransferase Clr4 is responsible for monomethylation, dimethylation, and trimethylation of H3K9 [
35]. Moreover, deletion of
BcDim5 in
B. cinerea resulted in the loss of H3K9me3 [
20]. Similar to previous results, western blot analysis indicates that H3K9me3 is nearly completely abolished in the
SsDim5 knockout mutant, and restored in the complemented strain, suggesting that
SsDim5 indeed functions as an H3K9 methyltransferase in vivo and affects H3K9me3 in
S. sclerotiorum. However, some studies have suggested that disruption of
Dim5 in
Beauveria bassiana leads to the loss of H3K9me3 and a significant reduction in H3K4me1/me2, H3K9me1/me2, and H3K36me2 [
36]. Therefore, it cannot be ruled out that
SsDim5 in
S. sclerotiorum may simultaneously possess other histone methylation modifications.
At different stages of growth, development, and infection, pathogens undergo extensive transcriptional reprogramming. H3K9 methylation is typically closely associated with heterochromatin formation, influencing various gene expressions by modulating chromosomal accessibility [
30]. Through a gene knockout study of
SsDim5 in
S. sclerotiorum, we discovered an indispensable role for
SsDim5 in the virulence of the fungus, without affecting its hyphal development and sclerotium formation. In the case of
N. crassa, the absence of
TrDim5 resulted in impaired vegetative growth and conidiation [
37]. The deletion mutant of
FvDim5 in
Fusarium graminearum (
ΔFvDim5) exhibited significant defects in conidiation, perithecium formation, and fungal virulence [
21]. Similarly, the disruption of
BcDim5 in
B. cinerea led to a significant reduction in hyphal growth, conidiophore production, and sclerotium yield, accompanied by a decrease in virulence [
20]. Undoubtedly, our research results once again confirm the conservation of Dim5 in facilitating the physiological function of pathogen infection. However, there are also differences existed in Dim5 function among different species. In
S. sclerotiorum, Dim5 is expressed at various stages of growth and development, yet no significant growth inhibition was observed through the observation of knockout mutants. Therefore, we speculate that
SsDim5 may have redundant functions in the regulation of growth and development processes in
S. sclerotiorum.
Genome sequencing has revealed that gene clusters involved in fungal secondary metabolism are often located near telomeres or heterochromatin regions [
38]. H3K9 methylation, closely associated with the establishment of heterochromatin [
39], has been demonstrated as an effective regulatory mechanism for disrupting or maintaining heterochromatin, impacting the generation of fungal secondary metabolites [
40,
41]. Here, RNA-seq analysis unveiled dysregulation in the expression of metabolism-related genes in the
SsDim5 knockout mutant, with significant enrichment observed in the pathway of mycotoxins biosynthesis through Gene Ontology (GO) analysis. Secondary metabolites (SM) play a crucial role in the virulence, development, and overall lifestyle of fungal pathogens. In the interaction between fungi and plant hosts, these metabolites may function as effectors while infecting the host and when recognized by the plant host during the infection process [
42]. Mycotoxins biosynthesis is a subset of secondary metabolites often overproduced in response to external stressors. The main factors that enhance mycotoxin production include oxidative stress, nutritional stress, light stress, as well as environmental factors such as pH, temperature, water activity, fungicides, and plant secondary metabolites [
43]. Studies have demonstrated that fungi regulate oxidative bursts through mycotoxins, enhancing ecological benefits [
44]. qRT-PCR results confirmed a significant decrease in the expression of genes related to mycotoxins biosynthesis in
SsDim5 knockout mutant. Additionally, the
SsDim5 knockout strain exhibited increased sensitivity to osmotic stress, oxidative stress, and substances damaging cell integrity. Intriguingly, disturbances in the redox state of
S. sclerotiorum affected the accumulation of OA [
45], impacting the levels of osmotic stress, high salt, and cell wall stress-related functional genes that influence virulence [
46,
47,
48,
49,
50]. Therefore, it can be inferred that in
S. sclerotiorum, the normal function of Dim5 is crucial for the synthesis of fungal secondary metabolites, especially mycotoxins, directly or indirectly influencing the virulence of
S. sclerotiorum to host plants.
Figure 1.
Phylogenetic tree, conserved protein domain and expression analysis of SsDim5. (A) Phylogenetic analysis of SsDim5 was conducted using the neighbour-joining method, aligning protein sequences with ClustalW based on JTT in MEGAX. Statistical confidence in the phylogenetic relationships was evaluated through 1000 bootstrap replicates. SsDim5 was marked with a black triangle, and UniProt database entry numbers were provided in brackets. (B) Conserved protein domains of SsDim5 and its orthologues were predicted using TBtool. (C) Expression levels of SsDim5 at different growth stages of S. sclerotiorum. Utilizing SsTubulin1 as the reference gene, average values and standard deviations were computed based on data from three independent biological replicates. Differences were evaluated using the one-way ANOVA test. * denotes P<0.05, ** denotes p < 0.01. (D) Expression levels of SsDim5 at different infection stages in B. napus leaves. Utilizing SsTubulin1 as the reference gene, average values and standard deviations were computed based on data from three independent biological replicates. Differences were evaluated using the one-way ANOVA test. * denotes p < 0.05, ** denotes p < 0.01.
Figure 1.
Phylogenetic tree, conserved protein domain and expression analysis of SsDim5. (A) Phylogenetic analysis of SsDim5 was conducted using the neighbour-joining method, aligning protein sequences with ClustalW based on JTT in MEGAX. Statistical confidence in the phylogenetic relationships was evaluated through 1000 bootstrap replicates. SsDim5 was marked with a black triangle, and UniProt database entry numbers were provided in brackets. (B) Conserved protein domains of SsDim5 and its orthologues were predicted using TBtool. (C) Expression levels of SsDim5 at different growth stages of S. sclerotiorum. Utilizing SsTubulin1 as the reference gene, average values and standard deviations were computed based on data from three independent biological replicates. Differences were evaluated using the one-way ANOVA test. * denotes P<0.05, ** denotes p < 0.01. (D) Expression levels of SsDim5 at different infection stages in B. napus leaves. Utilizing SsTubulin1 as the reference gene, average values and standard deviations were computed based on data from three independent biological replicates. Differences were evaluated using the one-way ANOVA test. * denotes p < 0.05, ** denotes p < 0.01.
Figure 2.
Phenotypes of SsDim5 knockout and complemented strains. (A) Phenotypes of the WT, SsDim5-1, SsDim5-2, and SsDim5-C. One representative biological replicate was shown. (B) Radial growth length assessment on PDA of WT, Ssdim5-1, Ssdim5-2, and SsDim5-C strains. Average values and standard deviations were computed based on data from three independent biological replicates. (C) Number of sclerotia per plate. Average values and standard deviations were computed based on data from ten independent biological replicates. Differences were evaluated using the one-way ANOVA test. (D) Sclerotial weight per plate. Average values and standard deviations were computed based on data from ten independent biological replicates. Differences were evaluated using the one-way ANOVA test.
Figure 2.
Phenotypes of SsDim5 knockout and complemented strains. (A) Phenotypes of the WT, SsDim5-1, SsDim5-2, and SsDim5-C. One representative biological replicate was shown. (B) Radial growth length assessment on PDA of WT, Ssdim5-1, Ssdim5-2, and SsDim5-C strains. Average values and standard deviations were computed based on data from three independent biological replicates. (C) Number of sclerotia per plate. Average values and standard deviations were computed based on data from ten independent biological replicates. Differences were evaluated using the one-way ANOVA test. (D) Sclerotial weight per plate. Average values and standard deviations were computed based on data from ten independent biological replicates. Differences were evaluated using the one-way ANOVA test.
Figure 3.
Pathogenicity assays of individual strains. (A) Disease phenotype of WT, Ssdim5-1, Ssdim5-2 and SsDim5-C leaves of B. napus. Photographs were taken at 24 hpi. One representative biological replicate was shown. Bar = 1 cm. (B) Statistical analysis of the lesion area in panels. WT, Ssdim5-1, Ssdim5-2 and SsDim5-C. Average values and standard deviations were computed based on data from three independent biological replicates. Differences were evaluated using the one-way ANOVA test. * denotes P<0.05, ** denotes p < 0.01. Bar = 500 μm. (C) The morphology of the appressoria of each strain after 16 h of inoculation on a glass slide under a stereomicroscope. (D) The morphology of the appressoria of each strain after 16 h of inoculation on the onions epidermis under a stereomicroscope. Bar = 200 μm. (E) Mycelium of WT, SsDim5-1, Ssdim5-2 and Ssdim5-C strains cultivated on PDA medium supplemented with bromophenol blue.
Figure 3.
Pathogenicity assays of individual strains. (A) Disease phenotype of WT, Ssdim5-1, Ssdim5-2 and SsDim5-C leaves of B. napus. Photographs were taken at 24 hpi. One representative biological replicate was shown. Bar = 1 cm. (B) Statistical analysis of the lesion area in panels. WT, Ssdim5-1, Ssdim5-2 and SsDim5-C. Average values and standard deviations were computed based on data from three independent biological replicates. Differences were evaluated using the one-way ANOVA test. * denotes P<0.05, ** denotes p < 0.01. Bar = 500 μm. (C) The morphology of the appressoria of each strain after 16 h of inoculation on a glass slide under a stereomicroscope. (D) The morphology of the appressoria of each strain after 16 h of inoculation on the onions epidermis under a stereomicroscope. Bar = 200 μm. (E) Mycelium of WT, SsDim5-1, Ssdim5-2 and Ssdim5-C strains cultivated on PDA medium supplemented with bromophenol blue.
Figure 4.
Western-blot analysis of S. sclerotiorum using antibodies histone H3 or trimethylated H3K9.
Figure 4.
Western-blot analysis of S. sclerotiorum using antibodies histone H3 or trimethylated H3K9.
Figure 5.
RNA-Seq analysis of WT and SsDim5 knockout strains. (A) The volcano plot of differentially expressed genes (DEGs) (SsDim5 knockout strains vs. WT). (B) Significantly enriched Gene Ontology (GO) terms for all DEGs. (C) The expression profiles of mycotoxins biosynthesis genes. (D) The expression levels of mycotoxins biosynthesis genes by qRT-PCR. Utilizing SsTubulin1 as the reference gene, average values and standard deviations were computed based on data from three independent biological replicates. Differences were evaluated using the one-way ANOVA test. * denotes p < 0.05, ** denotes p < 0.01.
Figure 5.
RNA-Seq analysis of WT and SsDim5 knockout strains. (A) The volcano plot of differentially expressed genes (DEGs) (SsDim5 knockout strains vs. WT). (B) Significantly enriched Gene Ontology (GO) terms for all DEGs. (C) The expression profiles of mycotoxins biosynthesis genes. (D) The expression levels of mycotoxins biosynthesis genes by qRT-PCR. Utilizing SsTubulin1 as the reference gene, average values and standard deviations were computed based on data from three independent biological replicates. Differences were evaluated using the one-way ANOVA test. * denotes p < 0.05, ** denotes p < 0.01.
Figure 6.
The growth of individual strains in the presence of various stressors. (A) Phenotypes of the WT, Ssdim5-1, Ssdim5-2 and SsDim5-C supplemented with 1 M NaCl, 1 M KCl, 1 M Sorbitol, 0.005% SDS, 10 mM H2O2, respectively. (B) Statistics of inhibition rates under various stresses. Average values and standard deviations were computed based on data from three independent biological replicates. Differences were evaluated using the one-way ANOVA test. * denotes p < 0.05, ** denotes p < 0.01.
Figure 6.
The growth of individual strains in the presence of various stressors. (A) Phenotypes of the WT, Ssdim5-1, Ssdim5-2 and SsDim5-C supplemented with 1 M NaCl, 1 M KCl, 1 M Sorbitol, 0.005% SDS, 10 mM H2O2, respectively. (B) Statistics of inhibition rates under various stresses. Average values and standard deviations were computed based on data from three independent biological replicates. Differences were evaluated using the one-way ANOVA test. * denotes p < 0.05, ** denotes p < 0.01.