Histone Methyltransferase SsDim5 Regulates Fungal Virulence through H3K9 Trimethylation in Sclerotinia sclerotiorum

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 H₂O₂, 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.

3. Results

3.1. Identification of S. sclerotiorum Histone H3K9 Methyltransferase

In present, H3K9me3 modification is considered as an epigenetic mark that manipulates gene expression by regulating chromosomal accessibility [30]. To identify the H3K9 methyltransferase in the S. sclerotiorum and investigate its biological functions, we utilized the B. cinerea Dim5 as a blast template and identified the homologous gene SS1G_01550 in the S. sclerotiorum (hereafter referred to as SsDim5). Gene structure analysis revealed that SsDim5 consists of three exons and two introns, encoding a peptide of 308 amino acid residues. Protein BLAST analysis indicated that Dim5 orthologs were widely distributed among species of the Sclerotinia, Monilinia, Fusarium, Colletotrichum, Alternaria, and Rutstroemia P. To analyze the phylogenetic relationships among Dim5 orthologs, we constructed a phylogenetic tree using Dim5 orthologs from the aforementioned species. The results indicated a close evolutionary relationship between SsDim5 and its ortholog in the B. cinerea (Figure 1A). Furthermore, protein structure analysis revealed that Dim5 orthologs from these species all contained the typical domains of H3K9 methyltransferase, including the PreSET domain, SET domain, and PostSET domain (Figure 1B). This suggests that these orthologs may possess histone modification activity, and the protein encoded by SS1G_01550 in S. sclerotiorum is an ortholog of the B. cinerea Dim5.

3.2. Expression Patterns of SsDim5 during Development and Infection Stages

In order to elucidate the regulatory functions of SsDim5 during the developmental and invasive stages of the S. sclerotiorum, quantitative real-time PCR was employed to examine the accumulation of its transcripts at different stages. The expression profile revealed distinct expression patterns of SsDim5 at various stages. Overall, during the developmental stage of the S. sclerotiorum, the expression of SsDim5 showed an increasing trend, reaching its highest level during the formation and maturation stages of the sclerotia (Figure 1C). Conversely, during different stages of S. sclerotiorum infection on B. napus leaves, the expression of SsDim5 was significantly suppressed initially at 9-24 h post-infection, followed by a significant upregulation at 48 h post-infection (Figure 1D). This indicates that SsDim5 might play a crucial regulatory role in the formation of sclerotia and the infection process on the host.

3.3. Generation of SsDim5 Knockout Mutants and Genetic Complementation Strains

To investigate the function of SsDim5 within the S. sclerotiorum, a homologous recombination strategy was employed to generate two independent knockout mutants, Ssdim5-1 and Ssdim5-2 (Figure S1A). PCR analysis confirmed the absence of the SsDim5 gene segment in both mutants, replaced by the HYG gene segment (Figure S1B). The results of qRT-PCR also indicated a lack of SsDim5 transcript in the knockout mutants (Figure S1C). Subsequently, the accuracy of gene complementation in the Ssdim5-1 mutant was validated by transforming the wild type SsDim5 gene segment, including a 1500 bp native-promoter, into the mutant strain. PCR and qRT-PCR further confirmed the accuracy of the gene complementation strain, designated as SsDim5-C (Figure S1B,C).

WT, two SsDim5 knockout mutants (Ssdim5-1/Ssdim5-2), and the complementation strain SsDim5-C exhibited similar nutritional growth and developmental phenotypes on PDA medium. Further statistical analysis indicated no significant differences in growth rates among the strains (Figure 2A, B). Likewise, the number of sclerotium produced and the average weight of each sclerotium by Ssdim5-1/Ssdim5-2 showed no significant differences compared to the WT and complementation strain on each culture plate (Figure 2C, D). These results imply that SsDim5 does not impact the mycelial growth and sclerotium formation of the S. sclerotiorum.

3.4. Deletion of SsDim5 Impairs the Virulence of S. sclerotiorum

To investigate the role of SsDim5 in pathogenicity, isolated leaves of B. napus were inoculated with the WT strain, SsDim5 knockout strains (Ssdim5-1/Ssdim5-2), and the complemented strain (SsDim5-C). After 24 h of infection, the lesion areas caused by the knockout strains Ssdim5-1/Ssdim5-2 were significantly smaller than that induced by the WT, whereas the lesion area caused by the complementation strain SsDim5-C resembled that of the WT (Figure 3A,B). Subsequently, each strain was employed to infect detached N. benthamiana leaves, and the infection outcomes mirrored those observed in B. napus leaves (Figure S2A,B). Notably, the virulence of the knockout strains Ssdim5-1/Ssdim5-2 was significantly diminished both on B. napus and N. benthamiana leaves, suggesting that the reduced virulence is not specific to a particular host species.

Appressorium formation and oxalic acid (OA) production are key strategies for regulating S. sclerotiorum virulence. Since the deletion of SsDim5 exhibits a reduced virulence phenotype, the OA production and appressoria of the knockout strain Ssdim5-1/Ssdim5-2 was tested. The results showed that the morphology and number of appressoria formed by the knockout strain Ssdim5-1/Ssdim5-2 were similar to those of the WT and complemented strain (SsDim5-C), with no significant difference, whether on a glass slide (Figure 3C) or onion epidermal cells (Figure 3D). In addition, by detecting OA production on PDA supplemented with bromophenol blue, as shown in Figure 3E, the knockout strain Ssdim5-1/Ssdim5-2 had the same phenotype as the WT and complemented strain SsDim5-C. These experiments indicate that SsDim5 plays a critical role in virulence, but not by affecting appressorium formation and OA production.

3.5. H3K9 Trimethylation Levels Are Significantly Reduced in SsDim5 Knockout Mutants

To investigate the role of SsDim5 in histone modification, Western blot analysis was employed to assess the levels of histone methylation in the WT strain, Ssdim5-1/Ssdim5-2, and SsDim5-C. When probed with an antibody against histone H3, all strains exhibited bands of similar size (Figure 4). However, noteworthy differences emerged when examining trimethylation at H3K9. Abundant specific bands corresponding to trimethylated H3K9 were detected in the WT strain and SsDim5-C, while the knockout strains Ssdim5-1/Ssdim5-2 showed a complete absence of signals for trimethylated H3K9 (Figure 4), indicative of a crucial role of SsDim5 in H3K9 trimethylation.

3.6. SsDim5 is Related to the Synthesis of Mycotoxins

To further investigate the biological functions of SsDim5 in S. sclerotiorum, whole-genome expression profiling of the WT and SsDim5 knockout strain (Ssdim5-1) mycelia was conducted through RNA sequencing (RNA-Seq). The results revealed 544 differentially expressed genes (DEGs), with 205 upregulated and 339 downregulated (p-value < 0.01 and fold change ≥2) (Figure 5A). Subsequently, an enrichment analysis based on Gene Ontology (GO) was performed on the DEGs, resulting in the assignment of 22 GO terms across three categories: biological processes (10 terms), cellular components (3 terms), and molecular functions (8 terms) (Figure S3). Notably, the biological process of “mycotoxins biosynthesis process” was significantly enriched in Ssdim5-1 (Figure 5B). According to the transcriptional profile, among the 6 genes involved in “mycotoxins biosynthesis process”, SS1G_13355, SS1G_07229, SS1G_02251, SS1G_13850, and SS1G_13358 exhibited significantly lower expression levels in Ssdim5-1 compared to the WT, while SS1G_13636 showed significantly increased expression in Ssdim5-1 (Figure 5C). qRT-PCR validation of these 6 genes in the WT, Ssdim5-1/Ssdim5-2 and SsDim5-C further confirmed the consistent modulation pattern. Specifically, SS1G_13355, SS1G_07229, SS1G_02251, SS1G_13850, and SS1G_13358 were suppressed in Ssdim5-1 and Ssdim5-2, while SS1G_13636 was enhanced in both Ssdim5-1 and Ssdim5-2 (Figure 5D). These findings suggest that SsDim5 regulates the synthesis of mycotoxins in S. sclerotiorum.

3.7. Regulation of SsDim5 in Response to Environmental Stresses

The synthesis of mycotoxins is predominantly induced by factors such as oxidative stress, nutritional stress and various external environmental stimuli. To assess the role of SsDim5 in the resistance of S. sclerotiorum to external environmental stimuli, the growth performance of the WT, Ssdim5-1/Ssdim5-2 knockout strains, and the complemented strain SsDim5-C was tested under different stress conditions. The growth inhibition rates of the Ssdim5-1/Ssdim5-2 knockout strains were significantly higher than those of the WT strain and the complemented strain SsDim5-C on PDA plates containing 1M NaCl, 1M KCl, and 1M Sorbitol (Figure 6A, B). This indicates that SsDim5 is associated with the high osmotic tolerance of S. sclerotiorum. Similarly, the growth inhibition rates of the Ssdim5-1/Ssdim5-2 strains were higher than those of the WT and the complemented strain SsDim5-C on PDA plates containing 10mM H₂O₂ and 0.02% SDS, suggesting that SsDim5 also plays a crucial regulatory role in oxidative stress and cell integrity in S. sclerotiorum (Figure 6A, B). Together, these results indicate that SsDim5 plays a crucial regulatory role in the response of S. sclerotiorum to environmental stressors.

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.

5. Conclusions

In summary, our results indicate that SsDim5 possesses histone H3 lysine methyltransferase activity, and plays a crucial regulatory role in the pathogenicity of S. sclerotiorum through the regulation of fungal mycotoxins biosynthesis, and the response to external stressors, thereby provides theoretical guidance for the development of new target points for the prevention and control of Sclerotinia stem rot.