- Review
- Open access
- Published:
Epigenetic weapons of plants against fungal pathogens
BMC Plant Biology volume 24, Article number: 175 (2024)
Abstract
In the natural environment, plants face constant exposure to biotic stress caused by fungal attacks. The plant’s response to various biotic stresses relies heavily on its ability to rapidly adjust the transcriptome. External signals are transmitted to the nucleus, leading to activation of transcription factors that subsequently enhance the expression of specific defense-related genes. Epigenetic mechanisms, including histone modifications and DNA methylation, which are closely linked to chromatin states, regulate gene expression associated with defense against biotic stress. Additionally, chromatin remodelers and non-coding RNA play a significant role in plant defense against stressors. These molecular modifications enable plants to exhibit enhanced resistance and productivity under diverse environmental conditions. Epigenetic mechanisms also contribute to stress-induced environmental epigenetic memory and priming in plants, enabling them to recall past molecular experiences and utilize this stored information for adaptation to new conditions. In the arms race between fungi and plants, a significant aspect is the cross-kingdom RNAi mechanism, whereby sRNAs can traverse organismal boundaries. Fungi utilize sRNA as an effector molecule to silence plant resistance genes, while plants transport sRNA, primarily through extracellular vesicles, to pathogens in order to suppress virulence-related genes. In this review, we summarize contemporary knowledge on epigenetic mechanisms of plant defense against attack by pathogenic fungi. The role of epigenetic mechanisms during plant-fungus symbiotic interactions is also considered.
Background
Due to their sedentary lifestyle, plants have to cope with constantly changing environmental conditions and also have to defend themselves against many biotic aggressors that threaten their development and reproduction. The response to various biotic stresses largely depends on the ability of the plant to rapidly modulate the transcriptome. External signals are translocated to the nucleus to activate transcription factors, resulting in increased expression of specific sets of defense-related genes. Among the mechanisms of transcription regulation in the plant response to stress, chromatin rearrangement (remodeling) emerges as a key process, based on modifications of histones by appropriate enzymes (deacetylases, acetylases, demethylases, methylases, enzymes carrying out ubiquitination) and ATP-dependent chromatin remodeling complexes. In addition, DNA methylation and non-coding RNAs also play an important role in modulating defense against stress factors.
Plant defense against fungal pathogens
The first line of plant defense against attack by pathogens is structural resistance, conditioned by the anatomical structure of the plant. The cell wall, made of polysaccharide (cellulose, hemicellulose and pectin) and non-polysaccharide (lignin) polymers and structural proteins and enzymes [1], constitutes the external barrier of plants against pathogen infections. The second line of defense is based on the production of constitutive secondary metabolites that are harmful to pathogens and activation of the hypersensitive response (HR), which is a mechanism employed by plants to prevent the spread of infection. The HR is characterized by the rapid death of cells in the local area surrounding the infection. The third line of defense is systemic acquired resistance (SAR) and induced systemic resistance (ISR). Plants often acquire systemic resistance to further infection after local infection with a pathogen [2, 3]. This requires accumulation of the plant hormone salicylic acid in tissue distal to the site of infection, which is part of the SAR, while the ISR is associated with JA/ET-dependent signaling [2]. ISR does not induce pathogenesis-related protein synthesis but can be induced by PAMP. Mobile immune molecules can be transferred from infected areas where a defense reaction has already started to distant healthy areas, triggering SAR there to protect the whole plant from infection [4]. These molecules affect the induction of pathogenesis-related genes (PR genes) and reprogram the cell to defend itself. Proper transcription reprogramming is believed to be critical for priming and plant defense against multiple pathogens [5,6,7].
The main goal of pathogens when infecting plants is to obtain nutrients necessary for their growth. For this purpose, in the first phase of infection, pathogens secrete plant cell wall-degrading enzymes (CWDE): polygalacturonases (PG), pectin lyases (PL) and pectin methylesterases (PME), as well as cellulases and hemicellulases [2]. In turn, the degradation of the plant cell wall, on the one hand, enables penetration and colonization of the host by pathogens, and on the other hand, as a result of the action of polygalacturonases secreted by pathogens, oligogalacturonans (pectin fragments) are released, which act as elicitors and activate plant defense mechanisms [8,9,10,11].
Pathogen-triggered molecular defense strategies of plants that lead to activation of the immune response and the signaling cascade have been divided into two main groups. The first of them is based on the perception of pathogen- and/or microbial-associated molecular patterns or damage-associated molecular patterns (PAMPs, MAMPs, DAMPs) and uses the recognition of a pathogen through plant cell surface anchored pattern recognition receptors (PRRs). Such recognition triggers a range of downstream defense mechanisms that lead to activation of pattern-triggered immunity (PTI). Examples of fungal MAMPs that interact with PRRs are chitin, chitosan, β-glucans, elicitins and ergosterol [12,13,14,15,16].
The second branch utilizes the recognition of microbial effectors or their activity, through resistance proteins (R proteins) to initiate effector-triggered immunity (ETI). Plant R genes typically encode intracellular receptors with nucleotide-binding leucine-rich repeat (NB-LRR) domains, organized with central NB and C-terminal LRRs. These NB-LRR-containing R proteins (NLRs) can be further classified into Toll/Interleukin1 receptor-like (TIR) or coiled-coil (CC) types according to their N-terminal sequences. NLRs can interact directly or indirectly with pathogen effectors to induce defense responses [12, 17,18,19,20]. It is generally accepted that, in contrast to PTI, ETI induces stronger and longer-lasting responses, and leads to a hypersensitive response (HR) resulting in synthesis of pathogenesis-related proteins (PR proteins). If this response fails, programmed cell death is induced [21, 22]. The other R gene-mediated resistance response is an oxidative burst that rapidly produces reactive oxygen species (ROS), which may have a direct antifungal effect or be a signal to activate other defense responses.
PTI and ETI affect the inhibition of pathogen multiplication by inducing various immune responses, including calcium ion signaling, nitric oxide and ROS production, alteration in membrane trafficking, activation of defense genes, cell death program limitation of nutrient transfer from the cytoplasmic matrix to the apoplast, biosynthesis of antimicrobial metabolites and defense phytohormones (salicylic acid (SA), jasmonic acid (JA) and ethylene (ET)), hydrolytic enzyme production, callose deposition at the infection site, stomatal closure and activation of the mitogen-activated protein kinase (MAPK) cascade, which leads to activation of WRKY type transcription factors, which activated by phosphorylation bind to the W-cassette of promoters of pathogenesis-related genes (PR genes), and changes in defense gene expression [23,24,25,26]. Epigenetic mechanisms such as DNA methylation, histone modifications and chromatin remodeling as well as epigenetic elements such as non-coding RNA contribute to modulation of gene transcription (Fig. 1) [27].
Histone modification
Gene expression is regulated by changes in chromatin structure. Less condensed regions of chromatin are transcriptionally active. Histones constituting the core of chromatin are subject to many modifications, e.g. acetylation, methylation, or ubiquitination (Fig. 2).
Histone acetylation disrupts interactions between nucleosomes, which leads to a looser state of chromatin and allows proteins involved in transcription to bind. In plants, lysine acetylation at the amino terminus of histone H3 and H4 tails has been linked to gene activation [28]. This modification reduces the ionic interaction between the positively charged lysine side chains and the negatively charged DNA backbone. Furthermore, lysine acetylation provides docking sites for transcriptional coactivator proteins containing bromodomains [29]. Histone acetyltransferases (HATs) are responsible for attaching the acetyl group of acetyl-CoA to the lysine amino group of histones, while histone deacetylases (HDACs) remove the acetyl group [27].
Both histone acetylation and deacetylation are involved in the plant response to pathogenic infections. Pathogen effectors can inhibit the activation of plant defense genes by e.g. interfering with the HAT function. During research on mechanisms underlying tomato responses to Botrytis cinerea, there was observed an increase in the H3K9 acetylation mark along the early induced genes SlyDES, SlyDOX1, and SlyLoxD encoding oxylipin-pathway enzymes, and SlyWRKY75 coding for a transcriptional regulator of hormonal signaling [30]. Song et al. identified the wheat histone acetyltransferase TaHAG1 as a positive regulator of resistance to powdery mildew caused by Blumeria graminis f. sp. tritici (Bgt). TaHAG1 acts as an epigenetic modulator and physically interacts with TaPLATZ5, a plant-specific zinc-binding protein. TaPLATZ5 directly binds to the promoter of TaPAD4 and together with TaHAG1 to potentiate the expression of TaPAD4 by increasing the levels of histone H3 acetylation. Acetylation can regulate expression of the key transducer gene TaPAD4 and promote accumulation of SA and reactive oxygen species for resistance to Bgt infection. A core transducer including the PAD4-EDS1 node is proposed as a convergence point between PTI and ETI [31].
One deacetylase that plays a role in immunity is histone deacetylase 19 (HDA19). HDA19 responds to injury and JA action. HDA19 overexpression in Arabidopsis thaliana (35S:HDA19) plants was associated with decreased histone acetylation levels, while HDA19 repressed (HDA19-RNAi) plants had increased histone acetylation levels. Compared to wild-type plants, the 35 S:HDA19 transgenic plants were more resistant to the pathogen Alternaria brassicicola. In these plants, increased expression of genes related to pathogenesis (glucanases and chitinases) was noted, supporting a role for HDA19 in the ethylene and JA-mediated defense response. The expression of these genes was downregulated in HDA19-RNAi plants [32]. HDA19 can bind directly to PR promoters and then deacetylate histones, ensuring low expression of PR genes under normal conditions and also avoiding overexpression of these genes during a defense response [33]. HDA6 can also inhibit the expression of pathogen-responsive genes, is induced by treatment with JA and an ET precursor, and interacts with the F-box protein, coronatine insensitive 1 (COI1), which mediates JA signaling [32, 34].
The histone deacetylase TaHDT701 was identified as a negative regulator of wheat defense responses. TaHDT701 associates with the RPD3 type histone deacetylase TaHDA6 and the WD40-repeat protein TaHOS15 to constitute a histone deacetylase complex. Knockdown of TaHDT701, TaHDA6, and TaHOS15 resulted in enhanced wheat powdery mildew resistance, suggesting that this histone deacetylase complex negatively regulates wheat defense responses. This was accompanied by increased histone acetylation and methylation [35]. Similar studies were performed on rice, where the transcription of H4-HDT701 of rice, a member of the HD2 deacetylase family, was increased in the compatible reaction and decreased in the incompatible reaction after infection by the fungal pathogen Magnaporthe oryzae. Overexpression of this deacetylase in transgenic rice decreased H4 acetylation levels and increased plant susceptibility to M. oryzae. HDT701 silencing led to increased levels of H4 acetylation and overexpression of the pattern recognition receptor (PRR) and other defense genes. Silenced transgenic lines showed increased M. oryzae as well as Xanthomonas oryzae pv oryzae (Xoo) resistance. HDT701 can also bind to defense-related genes to regulate their expression. The authors conclude that HDT701 negatively regulates innate immunity by modulating the levels of histone H4 acetylation of PRR and defense-related genes in rice [36].
Research on this deacetylase in wheat was carried out by the team of Zhi et al., who identified it as a negative regulator of wheat defense responses to Blumeria graminis f. sp. tritici (Bgt). They observed that TaHDT701 associated with the RPD3 type histone deacetylase TaHDA6 and the WD40-repeat protein TaHOS15 to constitute a histone deacetylase complex. Knockdown of this complex resulted in enhanced wheat powdery mildew resistance. In the absence of the TaHDT701-TaHDA6-TaHOS15 histone deacetylase complex, chromatin at the defense-related genes resides in a primed gene expression state marked by increased H4K16Ac, H3K9Ac, and H3K4me3, as well as reduced nucleosome occupancy, thereby stimulating defense-related transcription and defense responses to Bgt [35].
The HDAC proteins in the Sir2 family, which are NAD+-dependent HDACs, play various roles in plant physiological processes. In rice, OsSRT1 RNA interference induced an increase of histone H3K9 (lysine-9 of H3) acetylation and a decrease of H3K9 dimethylation, leading to H2O2 production, DNA fragmentation, cell death, and lesions mimicking plant hypersensitive responses during interactions with pathogens. Chromatin immunoprecipitation assays showed that OsSRT1 down-regulation induced histone H3K9 acetylation on the transposable elements and some of the hypersensitive response-related genes [37].
Research by Mengel et al. on A. thaliana indicated that NO also affects histone acetylation by targeting and inhibiting histone deacetylase complexes, resulting in hyperacetylation of specific genes, e.g. involved in plant defense against pathogens. Furthermore, salicylic acid, which is the major plant defense hormone against biotrophic pathogens, inhibited HDAC activity and increased histone acetylation by inducing endogenous NO production [38].
An important chromatin modifier that plays a role in histone acetylation and plant immunity is the Elongator complex. The Elongator complex is a large multi-subunit complex involved in one of the transcription stages – elongation – and many physiological processes in eukaryotes. The results of Defraia et al. indicate that the entire Elongator complex, in A. thaliana, is involved in basal immunity and ETI, but not in SAR, and suggest that it may play a role in facilitating transcriptional induction of defense genes through alterations to their chromatin, such as histone acetylation. The authors described Elongator Subunit 2 (AtELP2) as an accelerator of immune responses in A. thaliana since the Atelp2 mutant presents a delayed and reduced defense response. AtELP2 regulates cytosine methylation and histone acetylation levels on several defense genes, including several involved in responses to necrotrophic pathogens, such as B. cinerea and A. brassicicola [39, 40]. Research by Wang et al. showed that Elongator is required for full induction of the JA/ET defense pathway marker gene and for resistance to the necrotrophic fungal pathogens B. cinerea and A. brassicicola. The mutant with a loss-of-function mutation in the A. thaliana Elongator subunit 2 presents compromised resistance to the studied fungi, together with decreased histone acetylation and expression of JA/ET-defense genes WRKY33, ORA59 and PDF1.2 [41].
In the case of histone methylation, the situation is more complicated, because in addition to lysine residues, arginine residues can also be methylated in several places. Moreover, the matter is further complicated by the fact that specific methylation patterns are associated with both gene activation and repression [42]. In addition, there can occur mono-, di- and tri-methylation of histones, which present different physical properties. The strongest correlation between histone methylation and gene activity is found for trimethylation of Lys 4 on histone H3 (H3K4me3) on promoters and coding sequences of active genes [43]. H3K9me2, H3K27me3 and H2Aub1 are associated with transcriptional silencing/repression [44].
Many plant methylases and demethylases are involved in defense against pathogens. H3K36 methylation is involved in the regulation of PR gene transcription. In A. thaliana, set-domain group 8 (SDG8) is the main histone lysine methyltransferase catalyzing H3K36me2/3 [45]. SDG8 plays an important role in plant defense against fungal pathogens by modulating several genes in the JA and/or ET signaling pathways. Loss of SDG8 function causes mutant plants to show reduced resistance to the fungal pathogens A. brassicicola and B. cinerea [46]. In tomato, SDG8 orthologs were required for pathogen- and stress-induced enrichment of H3K36me3 and H3K4me3 at target genes, but their loss increased plant resistance to B. cinerea [47]. In addition to SDG8, SDG5 may also play role in plant immunity. Studies by Lee and co-authors indicate that in A. thaliana, SDG8 and SDG25 contribute to plant resistance to B. cinerea and A. brassicicola, either directly through histone lysine methylation or indirectly through H2B ubiquitination and by regulating plant resistance gene expression, lipid accumulation, carotenoid biosynthesis and maintaining integrity epidermis. However, SDG8 and SDG25 seem to have different molecular functions. SDG8 appears to perform the deposition of H3K36me2 and H3K36me3, while SDG25 may deposit H3K4me1 [48].
Histones can also be ubiquitinated on lysine residues, usually associated with transcriptional activation [44, 49]. A. thaliana RING E3 Histone Ligase Monoubiquitination 1 and 2 (HUB1 and HUB2) mediate the monoubiquitination of histone H2B. A. thaliana mutants deficient in HUB1 showed higher sensitivity to B. cinerea and A. brassicicola, while overexpressing HUB1 showed resistance to B. cinerea. The thickness of epidermal cell walls (acting as a physical barrier against invading fungal pathogens) is reduced in the hub1 mutant [50]. A thaliana H2Bub is also involved in regulating the dynamics of microtubules during the defense response to toxins produced by the necrotrophic pathogen Verticillium dahlia, likely through the protein tyrosine phosphatase-mediated signaling pathway [51]. HUB1/HUB2 homologs in tomatoes, SlHUB1/SlHUB2, have similar H2B monoubiquitination E3 ligase activity, and their silencing leads to increased plant susceptibility to B. cinerea [52].
In rice, histone lysine 2-hydroxyisobutyrylation (Khib) is involved in regulation of PR gene expression. The ascomycete Ustilaginoidea virens may affect regulation of Khib in rice, where most of the Khib sites in histone H3 were downregulated during infection. RPD3-like histone deacetylase-HDA705 is involved in the removal of Khib in rice, negatively regulating plant resistance to Ustilaginoidea virens, while HDA705 knockout increases resistance to these pathogens [53].
Chromatin remodelers
Chromatin remodeling not involving histone modifications is carried out by ATP-dependent chromosome remodeling complexes (CRC). Chromatin-remodeling factors can be divided into four subfamilies based on the catalytic Snf2 domain and other accessory domains: Switch/Sucrose on-fermentable (SWI/SNF), Inositol auxotrophy 80 (INO80), Imitation switch (ISWI), and Chromodomain helicase DNA-binding (CHD). In plants, members of these protein subfamilies have been associated with defense function [54] (Fig. 3).
Brahma (BRM) and Splayed (SYD) are the two most well-studied plant Snf2 proteins (belonging to the SWI/SNF subfamily). These proteins play roles in regulating plant developmental processes. SYD is required for expression of selected genes downstream of the jasmonate and ethylene signaling pathways. SYD is also directly recruited to the promoters of several of these genes. SYD contributes to defense against the necrotrophic pathogen B. cinerea in A. thaliana [55]. Research by Bezhani et al. indicated that A. thaliana SYD and BRM exhibit remarkable regulatory specificity by controlling the expression of a small number of targets involved in resistance [56].
Li et al. reported that the rice SWI/SNF2 ATPase gene BRHIS1 (BIT-responsive Histone-interacting SNF2 ATPase 1) was downregulated in response to a rice fungal pathogen (M. oryzae) or to the defense-priming-inducing compound BIT (1,2-benzisothiazol-3(2 h)-one,1, 1-dioxide). Their results showed that rice disease defense genes are initially organized in an expression-ready state by specific monoubiquitination of H2A and H2B variants deposited on their promoter regions, but are kept suppressed by the BRHIS1 complex, facilitating the prompt initiation of innate immune responses in response to infection through the stringent regulation of BRHIS1. Inhibition of the expression of BRHIS1 exhibited increased resistance to blast pathogens in rice, suggesting the negative regulatory role of BRHIS1 in plant immunity. After the plant receives a signal about the attacking pathogen, suppressed BRHIS1 expression provokes the inaccessible chromatin to become accessible for defense gene expression through enhanced monoubiquitination of those targeted histone variants. The expression profile of defense genes in plants with suppressed BRHIS1 expression indicates that there is a defense regulator that represses plant immunity in an SA-independent manner [57].
Another protein from the Snf2 family (INO80 subfamily) that may be involved in defense mechanisms is PIE1 (Photoperiod-Independent Early Flowering 1). PIE1 is part of the SWR1-like complex. One of the mechanisms involved in chromatin remodeling is the so-called ‘histone replacement’. An example of such a mechanism is the replacement of the canonical histone H2A with the variant histone H2A.Z, and in A. thaliana PIE1 is involved in this process. The H2A.Z variant interacts with PIE1. Loss of PIE1 function as well as depletion of H2A led to reduced basal resistance. There was a decrease in the level of expression of immunity genes; most of the mis-regulated genes were related to salicylic acid-dependent immunity. Functional knockout of these genes leads to constitutive activation of SA-mediated defense responses, and these mutants exhibit abnormal SA biosynthesis, overexpression of many SA-responsive genes, and spontaneous cell death under normal conditions, thus altering plant resistance to both biotrophic and necrotrophic pathogens. The SWR1-like complex is required for maintaining the repression of SA-dependent defense genes in unstressed plants [58, 59]. Cai et al. observed that the chromatin remodeling complex SWR1 enhanced resistance to the white mold fungus Sclerotinia sclerotiorum in A. thaliana via a process mediated by receptor-like kinase ERECTA (ER) signaling. The authors identified a series of WRKY33 target YODA DOWNSTREAM (YDD) genes and discovered that SWR1 and ER signaling was required to enrich H2A.Z histone variant and H3K4me3 histone modification at YDDs and the binding of WRKY33 to YDD promoters upon S. sclerotiorum infection. They also observed that the binding of WRKY33 to YDD promoters in turn promoted the enrichment of H2A.Z and H3K4me3 at YDD genes, thereby forming a positive regulatory loop to activate YDD expression. These results show the important role of chromatin structure in plant immune responses [60].
In A. thaliana, another ATPase of the INO80 subfamily, CHROMATIN REMODELER 19 (CHR19), has a conserved ATP-dependent nucleosome sliding activity. A variety of inducible genes, including several important genes in the salicylic acid and jasmonic acid pathways, were transcriptionally upregulated in the chr19 mutant under normal growth conditions, indicative of a role of CHR19 in transcriptional repression. In addition, the chr19 mutation triggered higher susceptibility to the JA pathway-defended necrotrophic fungal pathogen B. cinerea. CHR19 is involved in coordinating plant growth balance between development and stress response, and contributing to the improvement of plant resistance to fungal pathogens [61].
Decreased DNA Methylation 1 (DDM1) has a conserved SNF2 ATPase domain, which is required to maintain DNA methylation even though it has no methyltransferase activity itself [62] DDM1 controls R genes. DDM1 regulates expression of the A. thaliana plant resistance gene SNC1 (suppressor of npr1-1, constitutive1; one of the R genes) antagonistically to MOS1 (modifiers of snc1), by regulating the methylation levels in the upstream region of the gene and influencing the level of resistance to Hyaloperonospora arabidopsidis [63]. The repression of SNC1 expression in mos1 snc1 mutant plants of A. thaliana can be released by knocking out DDM1. It is thus likely that the repression of snc1 expression in mos1 snc1 is caused by altered chromatin structure rather than changes in DNA methylation. CHR5 (Chromatin-Remodeling Factor 5) acts as an antagonist to DDM1 in the regulation of SNC1. CHR5 also functions independently with the histone mono-ubiquitinase HUB1 in SNC1 regulation [64].
DNA methylation
DNA methylation and demethylation affect the resistance of plants to fungal diseases, probably due to the regulation of many genes of the immune response. DNA cytosine methylation is one of the major epigenetic mechanisms in higher eukaryotes and plays a key role in maintaining genome stability and regulating gene expression. DNA methylation targets tandem and interspersed repeats. At these loci, DNA methylation occurs in three different sequence contexts: symmetrical CG dinucleotides; symmetrical CHG, where H corresponds to A, T, or C; and asymmetrical. DNA methylation can also be found in gene bodies, exclusively at CG residues [65]. Cytosine methylation is established in all sequence contexts by de novo methyltransferases (DRM1/2) via the small interfering RNAs (RdDM) DNA methylation pathway. Here, DICER-dependent 21- to 24-nt siRNAs direct Argonaute proteins (AGO4/AGO6) to complementary sequences in the genome, possibly through a base-pairing mechanism and siRNA:nascent RNA base pairing mechanism to guide cytosine methylation. CG and CHG methylation is maintained through DNA replication by MET1, a homolog of the mammalian DNA methyltransferase DNMT1, and the plant-specific CMT3 methyltransferase. Active demethylation of methylcytosines is catalyzed by the DEMETER (DME) family of DNA glycosylases, while passive demethylation of methylcytosines is a consequence of DNA replication [65, 66] (Fig. 4).
RNA-directed DNA methylation (RdDM) is an epigenetic control mechanism driven by small interfering RNAs (siRNAs) that influence gene function. NRPD2 encodes the second largest subunit of the plant-specific RNA Polymerases IV and V (Pol IV and Pol V), which are crucial for the RdDM pathway. NRPD2 is the Overexpressor of Cationic Peroxidase 1. Lopez et al. evaluated the resistance of effective RdDM-deficient mutants of A. thaliana. Mutants of nrpe1, nrpd2, rdr2, drd1, ago4-2, and drm1drm2, except for nrpd1, exhibited increased susceptibility to B. cinerea. The results indicated that the RdDM pathway positively regulates resistance against necrotrophic fungi, presumably by promoting JA signaling during infection. The authors also suggested that while the Pol V complex is required for plant resistance, Pol IV appears to be redundant [67].
A. thaliana encodes four DNA demethylases: DEMETER (DME), Repressor Of Silencing 1 (ROS1), DEMETER-Like 2 (DML2), and DML3. Le et al. reported that ROS1, DML2 and DML3 play a role in fungal disease resistance in A. thaliana. A triple DNA demethylase mutant, rdd (ros1 dml2 dml3) shows increased susceptibility to the fungal pathogen Fusarium oxysporum. The expression of a significant number of genes involved in plant resistance was downregulated in mutants, suggesting that DNA demethylases maintain or positively regulate the expression of stress response genes required for F. oxysporum resistance. Stress response genes with downregulated rdd are enriched with short sequences of transposable elements in their promoters. Many of these transposable elements and their surrounding sequences show localized changes in DNA methylation in rdd and an overall reduction in CHH methylation, suggesting that the RNA-directed DNA methylation responsible for CHH methylation may be involved in DNA demethylase-mediated regulation of stress response genes. The results obtained by the authors also suggest that DNA demethylases target short promoter TEs (transposable elements) sequences to regulate these stress response genes. It is possible that the products of these genes, required for the stress response, are detrimental to plant development if accumulated at high levels and therefore must be suppressed under normal growth conditions via promoter methylation. Under stress conditions, these genes are temporarily activated by the action of DNA demethylases [68]. Many types of TEs have been identified in R-gene clusters. TE activity is usually suppressed by methylation or expressed at low levels as a result of methylation along their entire length [69]. TE methylation sometimes spreads into the genic regions of flanking genes, and suppresses their transcription [70]. Co-expression of genes and TEs is possible, in particular for long terminal repeat (LTR) elements. For example, the rice LTR retrotransposon Renovator is present in the promoter region of the Pit resistance gene in specific rice cultivars. The functional allele of PitK59 contains four amino acid substitutions and has the LTR retrotransposon Renovator inserted upstream, which enables expression of the rice blast resistance gene and confers resistance to Magnaporthe grisea. The 5′ side of Renovator is heavily methylated, whereas the 3′ side is only partially methylated. LTR methylation status correlates with transcriptional activity, suggesting that hypomethylation of the 3′ side of the Renovator sequence may be needed for higher Pit expression [71]. The above-mentioned and other examples confirming the role of DNA methylation in the response of plants to infections with pathogenic fungi are presented in Table 1.
Non-coding RNA
In addition to the role of protein-coding genes as key regulators of plant immunity, increasing evidence indicates the importance of non-coding RNAs (ncRNAs) in plant immune responses. Depending on their mode of biogenesis and their functions, ncRNAs have evolved into various forms that include microRNAs (miRNAs), small interfering RNAs (siRNAs), lncRNAs, circular RNAs (circRNAs), and derived ncRNAs (Fig. 5). NcRNAs may act as epigenetic modulators through chromatin remodeling or regulate gene expression at the transcriptional or posttranscriptional level [79, 80].
One of gene regulatory mechanism is RNA silencing. This process involves sRNA. sRNAs are divided into two major classes: miRNAs and siRNAs. miRNAs are usually 21–24 nucleotides (nt) long and are derived from RNAs with imperfectly base-paired hairpin structures. siRNAs are generated from perfectly complementary long double-stranded RNAs (dsRNAs) and may require RNA-dependent RNA polymerases [81].
One of the first reported miRNAs involved in plant resistance was miR393 from A. thaliana. It was induced by a PAMP, flg22, and silenced auxin receptors to turn down the auxin signaling pathway, resulting in activation of PTI [82]. Cassava anthracnose disease is caused by the fungus Colletotrichum gloeosporioides, and can cause tip die-backs and stem cankers. Quantitative RT-PCR analysis revealed differential expression of the two miRNAs and their target genes in the two cassava cultivars that had been subjected to fungal infection. The more resistant variant revealed upregulation of miR160 and miR393, which consequently led to low levels of ARF10 and TIR1 transcripts, respectively. The more susceptible variety showed the opposite pattern. Cis-regulatory elements relevant to defense and stress responses, fungal elicitor responses, and hormonal responses were most prevalent in the promoter regions of miRNA genes [83]. In A. thaliana miR396 levels gradually decreased during fungal infection, thus enabling its growth-regulating factor (GRF) transcription factor target genes to trigger host reprogramming, enabling defense against Plectosphaerella cucumerina, a destructive necrotrophic fungal pathogen [84].
Comparison using deep sequencing of small RNA libraries from the resistant and sensitive variety O. sativa allowed Li et al. to identify a group of miRNAs that were differentially expressed after infection with M. oryzae. Moreover, transgenic rice plants overexpressing miR160a and miR398b displayed enhanced resistance to M. oryzae, as demonstrated by decreased fungal growth, increased hydrogen peroxide accumulation at the infection site, and up-regulated expression of defense-related genes. miR398b overexpression reduced the transcript levels of genes encoding superoxide dismutases (CSD1, CSD2, SODX, and CCSD), leading to elevated ROS production and enhanced plant resistance against M. oryzae [85, 86]. Another example of non-coding RNAs involved in the regulation of ROS production in response to fungal infection is miR400 in A. thaliana. This miRNA directs the cleavage of two pentatricopeptide repeat genes, resulting in increased ROS accumulation and impaired resistance to B. cinerea [87]. In A. thaliana interference with miR773 activity by target mimics (in MIM773 plants) and concomitant upregulation of the miR773 target gene methyltransferase 2 (MET2) increased resistance to infection by necrotrophic (Plectosphaerella cucumerina) and hemibiotrophic (F. oxysporum, Colletotrichum higginsianum) fungal pathogens. Upon pathogen challenge, these plants accumulated higher levels of callose and reactive oxygen species than wild-type plants [88].
Silencing of individual cDNA clones in wheat challenged with Puccinia striiformis f. sp. tritici indicated that the chemocyanin-like protein gene TaCLP1 (the target gene of wheat miRNA R408) positively regulated resistance to stripe rust [89]. After infection of barley with the fungus Blumeria graminis f. sp. hordei, there was significant reprogramming of gene expression mediated by Mla- Mildew resistance locus a. Xu et al. utilized a proteomics-based approach, combined with barley Mla, required for Mla12 resistance1 (rar1), and restoration of Mla resistance1 (rom1) mutants, to identify components of Mla-directed signaling. Loss-of-function mutations in Mla and Rar1 both resulted in the reduced accumulation of chloroplast copper/zinc superoxide dismutase 1 (HvSOD1), whereas loss of function in Rom1 re-established HvSOD1 levels. Mla and Rom repress miR398-mediated SOD1 expression to change the hypersensitive reaction response to fungus [90]. Campo et al. investigated miRNAs that are regulated by elicitors from the blast fungus M. oryzae in rice (O. sativa). Elicitor treatment was accompanied by dynamic alterations in the expression of a significant number of miRNAs. They reported a rice miRNA, osa-miR7695, which negatively regulates an alternatively spliced transcript of OsNramp6 (Natural resistance-associated macrophage protein 6). Overexpression of this miRNA in rice confers pathogen resistance [91]. Li et al. found that miR169 acts as a negative regulator in rice immunity against the blast fungus M. oryzae by repressing the expression of nuclear factor Y-A (NF-YA) genes [92]. Another miRNA – osa-miR164a – is also involved in the defense of rice against this fungus. Wang et al. found that the expression of osa-miR164a decreased after infection with M. oryzae in both early and late stages, which was associated with induced expression of its target gene, OsNAC60, which encodes a transcription factor whose overexpression enhances defensive responses such as increased programmed cell death, greater ion leakage, greater accumulation of reactive oxygen species and callose deposition, and upregulation of defense-related genes. In transgenic rice where the miR164a/OsNAC60 regulatory module was dysfunctional, plants developed significant susceptibility to infection by M. oryzae [93].
Rice plants with activated, conserved, polycistronic, miR166 expression showed enhanced resistance to infection by the fungal pathogens M. oryzae and Fusarium fujikuroi. Disease resistance was associated with stronger expression of defense responses during pathogen infection. miR166k positively regulates rice immunity by controlling expression of the EIN2 (ethylene-insensitive 2) gene, which encodes a critical regulator of ethylene signaling [93].
Several miRNA families target genes encoding plant innate immunity receptors with a leucine-rich nucleotide binding site (NBS-LRR). Cotton plants can induce expression of NBS-LRR defense genes by suppressing the miRNA-mediated gene silencing pathway (miR482) after attack by the fungal pathogen V. dahlia [94]. Ouyang et al. investigated the role of miRNAs in tomato defense against F. oxysporum using comparative miRNA profiling of susceptible (Moneymaker) and resistant (Motelle) tomato cultivars. The authors focused on two miRNAs – slmiR482f and slmiR5300 – which are exceptionally downregulated in the resistant variety during fungal infection. Two predicted mRNA targets each of slmiR482f and slmiR5300 exhibited increased expression in Motelle. All predicted targets of these miRNAs encode proteins with nucleotide-binding (NB) domains and a motif associated with plant pathogen resistance. None of these targets correspond to I-2, the only known resistance (R) gene for F. oxysporum in tomatoes, which supports roles for additional R genes in the immune response. Fungus infection down-regulates the accumulation of miRNA to increase the expression of NB domain genes [95].
Analysis by Ma et al. showed that one of the miRNAs, Md-miRLn11 (Malus domestica microRNA Ln11), targeted an apple nucleotide-binding site (NBS)-leucine-rich repeat (LRR) class protein coding gene (Md-NBS). Analysis of the expression of Md-miRLn11 and Md-NBS during the optimum invasion period in 40 apple varieties showed that expression of the Md-NBS gene in resistant varieties is higher than in susceptible varieties, with an inverse pattern for Md-miRLn11. Infection with the fungus-causing apple leaf blotch on a seedling of a resistant apple cultivar showed a marked decrease in apple resistance. On the other hand, a susceptible apple cultivar with infiltration of the Md-NBS gene showed a significant increase in disease resistance [96]. In Brassica napus, 41 lncRNAs were identified as precursors for microRNAs, including miR156, miR169 and miR394, with significant roles in mediating plant responses to the fungal phytopathogen Sclerotinia sclerotiorum [97]. The abovementioned and other examples confirming the participation of miRNA in the plant defense response against infection with pathogenic fungi are presented in Table 2.
Gene expression mediated by siRNA plays an important role in regulation of plant growth, development and immunity. For example, Qiao et al. characterized the Tourist-miniature inverted-repeat transposable element (MITE)-derived siR109944 as having a conserved function that enhanced susceptibility to Rhizoctonia solani infection by affecting auxin homeostasis in rice and A. thaliana. One potential target of siR109944 is the F-Box domain and LRR-containing protein 55, which encode the transport inhibitor response 1 (TIR1)-like protein. The authors found that rice had significantly enhanced R. solani susceptibility when siR109944 was overexpressed. Overexpression of siR109944 in A. thaliana also increased susceptibility to B. cinerea, Sclerotinia sclerotium, and V. dahliae infection [120, 121]. In rice, the majority of characterized sRNAs are within the range of 21 to 24 nt long. Niu et al., using deep sequencing on rice plants, identified a group of small interfering RNAs between 25 and 40 nt in length, which constitute another class of sRNA. The results indicate that some rice this group of siRNA are differentially regulated after infection with the fungal pathogen R. solani, which causes rice sheath blight disease, and may also target genes related to defense [122].
Plants contain a large number of intracellular nucleotide-binding NLR immune receptors, encoded by resistance (R) genes, which recognize specific pathogen effectors and trigger resistance responses. Disease-related genes, particularly the NLR class of R plant genes, can be triggered by miRNAs to generate phasiRNAs (phased small interfering RNA), which could reduce the transcript levels of their targets. phasiRNAs are phased, secondary, small interfering RNAs, mostly triggered by the cleavage of the target mRNAs by 22-nucleotide (nt) miRNAs. phasiRNAs are widely present in plants, and they can regulate their target genes by cleaving mRNAs at the post-transcriptional level or directing DNA methylation at the transcriptional level. phasiRNAs can suppress the expression of NLR genes when there is no pathogen. Upon pathogen infection, the phasiRNAs are generally downregulated, which leads to upregulation of NLRs and activation of immune responses [123, 124]. For example, Liu et al. observed in barley that 22-nt miR9863 triggered 21-nt phasiRNA biogenesis and together repressed the expression of group I Mla alleles. Barley Mla alleles encode coiled-coil (CC), nucleotide binding, and NLR receptors that trigger isolate-specific immune responses against the powdery mildew fungus, Blumeria graminis f. sp. hordei [125]. Hunt et al. also identified barley phasiRNAs in response to Bgh infection, which overlap transcripts that encode receptor-like kinases (RLKs) and nucleotide-binding, leucine-rich domain proteins (NLRs). These include phasiRNA loci that overlap with a significant proportion of receptor-like kinases, suggesting that an additional sRNA control mechanism may be active in barley leaves as opposed to predominant R-gene phasiRNA overlap in many eudicots [126].
The expression of some NLRs is also regulated by transcriptional silencing by heterochromatic small interfering RNAs (hc-siRNAs) produced in the RDR2-DCL3-AGO4 pathway by plant-specific RNA polymerase IV and V. hc-siRNAs are usually 24–30 nt in length, derived mainly from transposons and repeats, and direct de novo DNA methylation and histone modifications at their target sites [127]. An example is the Rice Pygm locus, which confers lasting resistance to the fungus M. oryzae. PigmS expression and thus PigmR-mediated resistance are subjected to tight epigenetic regulation. The PigmS promoter contains two tandem miniature transposons, which associate with hc-siRNAs [128].
Not only short non-coding RNAs play an important role in the regulation of gene expression and response to biotic stress, but also long non-coding RNAs (lncRNAs) may perform such a function. lncRNAs are transcribed from a variety of genomic locations (introns, intergenic spaces, and coding regions) from the sense or antisense strand. lncRNAs function in cis or in trans and affect gene regulation transcriptionally or posttranscriptionally by diverse mechanisms, including recruiting factors that activate transcription or modify chromatin, serving as precursors of small RNAs, and even potentially affecting nuclear architecture. Plant lncRNAs function in RNA-directed DNA methylation (via production of small interfering RNAs) [129]. Bhatia et al., using a computational approach based on RNA-seq data, identified 71 powdery mildew-responsive Vitis vinifera lncRNAs. Further analysis revealed lncRNA association with Ca2+-binding proteins such as calmodulin/calmodulin-like proteins, enzymes involved in reactive oxygen species (ROS) metabolism, cell-wall modification/reinforcement, secondary metabolic pathways, phytoalexin (like resveratrol) production, pathogenesis-related proteins such as PR-1, PR-4 and PR-10, and phytohormone-based signal transduction [130]. Comparison of melon varieties sensitive and resistant to this pathogen also confirmed the role of lncRNA in plant defense against fungi [131].
With a strand-specific RNA-sequencing approach, Zhu et al. detected several lncRNAs induced in A. thaliana in response to infection with the fungus F oxysporum. Several noncoding natural antisense transcripts responsive and 20 responsive long noncoding transcriptionally active regions (lncTARs) to F. oxysporum infection were found in genes implicated in disease defense. Ten F. oxysporum-induced lncTARs were functionally characterized using T-DNA insertion or RNA-interference knockdown lines, and five were demonstrated to be related to disease development. Knockout mutants for these lncTARs presented higher susceptibility to this fungus, indicating their role in the regulation of defense. Regulation of the At2g30770 gene and its long noncoding natural antisense transcripts (lncNATs) was coordinated in response to F. oxysporum infection. The At2g30770 gene contains a TCA-element responsive to salicylic acid, and a TC-rich repeat responsive to stresses in its core promoter. Additionally, promoter analyses suggested that some of the F. oxysporum-induced lncTARs are direct targets of transcription factor(s) (TF) responsive to pathogen attack. Sense transcripts involved in the defense response and their lncNATs could be expressed in the same direction, which could be dependent or independent owing to the presence of similar TF binding sites upstream and downstream of protein-coding genes or similar pathogen-responsive elements in their promoter [132]. Sense and antisense strands were transcribed in the opposite direction and a negative correlation was observed between sense and antisense strand expression in both susceptible and resistant wheat cultivars, in response to powdery mildew infection. The comparison of lncRNA expression profiles between wheat species sensitive and resistant to Blumeria graminis f. sp. tritici exposed to this pathogen revealed that numerous lncRNAs were differentially expressed. Among them, some were precursors of small RNAs such as microRNAs and siRNAs, two long non-protein-coding RNAs (npcRNAs) were identified as signal recognition particle (SRP) 7 S RNA variants, and three were characterized as U3 small nucleolar RNAs (snoRNAs). In addition, lncRNAs were expressed in a tissue-specific manner [133]. In mulberry, the lncRNA MuLnc1 was found to be cleaved by mulmiR3954. It was observed that one of the siRNAs produced, si161579, could silence the expression of the calmodulin-like protein gene CML27 of mulberry (MuCML27). When MuCML27 was heterologously expressed in A. thaliana, the transgenic plants exhibited enhanced resistance to B. cinerea [134]. The most commonly cultivated grapevine species – V. vinifera – is susceptible to many pathogens, one of which is B. cinerea. The results obtained by Bhatia et al. indicate that lncRNAs, along with other regulators such as miRNAs, help modulate the basic plant defense response to necrotrophic fungi. The lncRNAs identified by the authors may coordinate in a grapevine defense response to B. cinerea that includes fungal chitin degradation, stilbenoid accumulation, ROS detoxification, and cell wall reinforcement [135]. Other research conducted on vines infected with Lasiodiplodia theobromae also confirmed the role of lncRNA in defense against fungal pathogens by affecting the target genes involved in cell wall organization and chitin signaling [136]. Zhau et al. identified 539 lncRNAs from powdery mildew-resistant (MR-1) and susceptible melon (Top Mark), of which 254 were significantly altered after fungal infection. These lncRNAs might be involved in the hydrolysis of chitin, callose deposition, plant-pathogen interaction pathway and the plant hormone signal transduction pathway [137]. The above-mentioned and other examples confirming the involvement of lncRNA in the plant defense response to fungal infections are presented in Table 3.
“Cross-kingdom RNAi”
Another important aspect is the fact that most sRNAs work endogenously, but some can travel across organismal boundaries between hosts and pathogens and silence genes in trans in interacting organisms, a mechanism called ‘’cross-kingdom RNAi’’. During the arms race between fungi and plants, some fungi send sRNA as an effector molecule to plant cells to silence plant resistance genes, while plants also transport sRNA mainly using extracellular vesicles to pathogens to suppress virulence-related genes. Cai et al. reported that A. thaliana host cells secreted exosome-like extracellular vesicles to deliver sRNA to the fungal pathogen B. cinerea. These sRNA-containing vesicles accumulate at sites of infection and are taken up by fungal cells. Transferred host sRNAs induced silencing of fungal genes critical for pathogenicity [148]. Zhang et al. noted that in response to V. dahliae infection (a vascular fungal pathogen responsible for wilt diseases in many crops), cotton plants increased the production of microRNAs 166 (miR166) and miR159 and exported both microRNAs to fungal hyphae for specific silencing of virulence genes (Ca2+-dependent cysteine protease calpain clp-1 (Clp-1) and isotrichodermin C-15 hydroxylase (HiC-15) genes). Clp-1 and HiC-15 transcript levels were indeed downregulated in recycled cotton shreds infected with V. dahliae. Clp-1 and HiC-15 transcript levels were downregulated in V. dahliae hyphae recovered from infected cotton [149]. Another example of a plant microRNA affecting a fungal pathogen is tomato miR1001. The study by Meng et al. demonstrated the regulatory role of this miRNA in the growth and development of B. cinerea. The results showed that miR1001 inhibited the virulence of B. cinerea in infected plants. In addition, miR1001 had a significant inhibitory effect on the germination of B. cinerea conidiospores. The results indicated that miR1001 can directly target the genes Bcin03g02170.1 and Bcin10g01400.1, which encode ATP-dependent metallopeptidase and cysteine-type endopeptidase, respectively, in B. cinerea [150].
Aggressive fungal pathogens such as B. cinerea and Verticillium spp. cause serious crop losses worldwide. B. cinerea has been found to deliver small RNA (Bc-sRNA) into plant cells to silence host resistance genes. Such sRNA effectors are mostly produced by B. cinerea Dicer-like protein 1 (Bc-DCL1) and Bc-DCL2. Wang et al. found that sRNA expression targeting Bc-DCL1 and Bc-DCL2 in A. thaliana and tomatoes silenced Bc-DCL genes and impaired fungal pathogenicity and growth, exemplifying bidirectional cross-kingdom RNAi and sRNA trafficking between plants and fungi [151].
Increasing evidence has shown that the regulatory function of ncRNAs plays a role in the association of plants with nonpathogenic symbiotic microorganisms. In research on the beneficial fungal endophyte Fusarium solani strain K (FsK), transformed FsK with a hairpin RNA (hpRNA) construct designed to target a reporter gene in its host Nicotiana benthamiana, after inoculation with the host plant, showed systemic RNA silencing and DNA methylation of the host reporter gene. The hpRNA was processed by FsK RNAi machinery presumably into 21-24-nt small RNAs that triggered RNA silencing but not DNA methylation in the fungal hyphae. These data suggest that RNAi signals can be translocated by endophytes to hosts and can modulate gene expression during mutualism [152].
Symbiotic fungi and epigenetic modification in plants
In addition to pathogenic fungi, a large proportion of fungal endophytes do not negatively affect the host plant (commensalism) or have a positive effect on plants (mutualism). These types of interactions are called symbiosis. Furthermore, the symbiotic relationship between plant roots and fungi is known as mycorrhiza. Symbiotic endophytes can counteract the development of host plant pathogens, for example by inducing host defense mechanisms or producing compounds that inhibit the growth of other microorganisms. To create and maintain symbiosis, constant communication between the mycobiome and the host plant is required. Furthermore, the plant will have to distinguish whether the microorganism is a friend or a foe. Environmental signals can induce epigenetic regulations that can modulate the host plant’s interaction with microorganisms or trigger defense [153,154,155,156].
Endosymbiotic relationships can be controlled by DNA methylation. Varga and Saulsbury carried out an experiment using Geranium robertianum in symbiosis with the arbuscular mycorrhizal fungus Funneliformis mosseae. This study demonstrated that colonization by an arbuscular mycorrhizal fungus can influence DNA methylation levels in hosts [157]. Moreover, the appropriate methylation status of the host’s DNA is important to create and regulate symbiotic interactions. Vigneaud et al. showed that host DNA hypomethylation limits the formation of ectomycorrhizas on the example of Populus spp. and the fungus Laccaria bicolor [158]. Research by Hubbard et al. showed that in the endosymbiotic seed-fungus relationship, colonization of the endophytic fungus SMCD 2206 was associated with changes in DNA methylation in wheat subjected to drought stress. Such epigenetic changes can be associated with increased plant resistance to drought. DNA methylation patterns observed in drought-stressed wheat seedlings co-cultured with SMCD 2206 resembled those of the unstressed control [159]. Research by Forte et al. on the mutualistic interaction between Epichloë sp. LpTG-3 strain AR37 and L. perenne showed that the presence of the endophyte led to a decrease in DNA methylation across genomic features in the host, with differentially methylated regions primarily located in intergenic regions and CHH contexts [160]. Etemadi et al. investigated the role of miR393, targeting several auxin receptors, during the colonization of Solanum lycopersicum, M. truncatula and O. sativa roots by arbuscular mycorrhizas. The results showed that miR393 is a negative regulator of arbuscule formation by hindering auxin perception in cells containing arbuscules [161]. Endosymbiosis of legumes and arbuscular mycorrhizal fungi is regulated by NSP2, which is a target of microRNA171h (miR171h). The spatiotemporal expression of miR171h and NSP2 is closely related to the nutritional status of the plant [162]. Results of research on M. truncatula conducted by Lauressergues et al. suggest that there is a regulatory mechanism involving miR171h-mediated negative regulation of NSP2, triggered by lipochito-oligosaccharides. This mechanism prevents excessive colonization of roots by arbuscular mycorrhizal fungi [163]. Furthermore, a microRNA targeting NSP2 (miR171h) is also rapidly induced by cytokinins and then shows an expression pattern anticorrelated with NSP2 [164]. Another M. truncatula miRNA, miR169a, was found to target MtHAP2 expression to control nodule cell differentiation and nodule development [165]. Analysis of the Phaseolus vulgaris genome allowed the identification of a set of small RNAs potentially important for regulating symbiosis. Among them, six, including miR-RH82, were differentially expressed in response to treatment with nodulation agents and may play a critical role during the early stages of symbiosis [166]. Studies in maize have shown that lncRNAs are involved in beneficial interactions between plants and the microfungus Rhizophagus irregularis. 63 lncRNAs were differentially expressed. The target genes of differentially expressed lncRNAs were mainly related to transmembrane transport of phosphate ions, the cellular response to potassium ion deficiency, and lipid catabolic processes [167]. It appears that suppression by epigenetic modifications of the plant defense system during symbiosis, as well as epigenetic changes in the symbiont genome, is crucial for promoting symbiosis and maintaining beneficial interactions. However, this still needs to be explored through further research.
Epigenetic memory
The immune memory is associated with the process of preparing cells (priming), consisting in strengthening the defensive response of plants previously in contact with the stress factor, activated only in response stress. Priming works at the phenotypic level do not change the DNA sequence, and can be reversible. The plant, under the influence of stress, corresponds to changes in gene expression: faster/earlier, more strongly, or expression of genes not expressed in a healthy plant. Epigenetic changes provide plants with immune memory [19, 44, 168, 169]. Among forms of epigenetic immune memory of plants, we can specify somatic memory (within the life cycle of a plant), intergenerational memory (stable for the first generation), and transgenerational memory (stable for at least two generations) [170].
Epigenetic memory leading to increased resistance to fungal pathogens can be induced by priming agents, including chemical compounds: β-aminobutyric acid (BABA), xenobiotic chemicals: benzothiadiazole (BTH) or (R)-β-homoserine (RBH), hormones (SA, JA, MeJA), microbes: Bacillus, Pseudomonas, plant-growth-promoting rhizobacteria (PGPR): Trichoderma spp., nonpathogenic strains of Fusarium spp., Piriformospora indica, and arbuscular mycorrhizal fungi (AMF) from the genus Glomeromycota and abiotic stress (salt, cold, heat) [7, 171, 172]. Seeds and plants could be primed. Seed priming increases seedlings’ tolerance of environmental stresses. Primed plants respond faster and/or more strongly to recurring defense stimuli. Some priming states are relatively short-term and disappear within a few days, while others are long lasting and can even be transmitted between plant generations (Fig. 6).
Exposure of plants to sub-lethal levels of salt, cold, or heat can enhance their resistance to infection by pathogens. In such cross-priming, PTI-responsive genes are enriched by histone acetylation associated with their transcriptional activation [23].
The existing literature provides a wealth of information on plant defense priming and immune memory, primarily focusing on the resistance of plants to bacterial or oomycete pathogens. However, there is a noticeable lack of information regarding the specific epigenetic mechanisms involved in these processes when it comes to fungal pathogens. In this review our purpose was to exclusively present examples demonstrating that seed or plant priming operates through epigenetic mechanisms, thereby enhancing plant resistance to fungal pathogen infections.
Epigenetic mechanisms underlying priming are not well understood. Catoni et al. investigated epigenetic changes in tomato plants primed for pathogen resistance by treatment with β-aminobutyric acid. Genomes of plants treated with BABA showed a significant reduction in global cytosine methylation, especially in CHH sequence contexts. Differentially methylated region (DMR) analysis showed that the CHHs in the DMR were almost exclusively hypomethylated. A large number of DMRs have been reported in gene promoters that are differentially expressed in response to infection with B. cinerea. However, most genes that showed priming did not contain DMR, and the overall distribution of methylated cytosines in priming genes was not altered by BABA treatment. It is possible that the BABA treatment of tomato seedlings causes characteristic changes in genome-wide DNA methylation, and that CHH hypomethylation only affects a minority of genes showing primed responses to pathogen infection. Probably methylation may affect priming via in-trans regulation, acting at a distance from defense genes, and by targeting a smaller group of regulatory genes controlling stress responses [173].
A lot of research has focused on seed priming [174], but the process can be successfully performed on whole plants. A good example is the study of Wojtasik et al., where flax seedlings were treated with the non-pathogenic strain F oxysporum Fo47 to increase resistance to the pathogenic strain of F oxysporum L. f. sp. lini. Two contrasting effects on the levels of methylation in flax were detected for both types of Fusarium strain infection: genome-wide hypermethylation and hypomethylation of resistance-related genes (β-1,3-glucanase and chitinase). Despite the differences in methylation profile, the expression of these genes increased. The peak of demethylation correlated with the alteration in gene expression induced by the non-pathogenic strain. In the case of pathogen infection, the expression peak lagged behind the gene demethylation. Plants pretreated with the non-pathogenic strain memorized the hypomethylation pattern and then reacted more efficiently upon pathogen infection [175].
Not only DNA methylation can affect priming but also processes related to non-coding RNAs play an important role. Wang et al. studied the small RNAs and transcriptome changes in maize leaves that were systemically induced by seed treatment with Trichoderma harzianum (strain T28) against Cochliobolus heterostrophus infection in leaves. Analysis of the sequencing data showed that genes involved in the plant hormone signal transduction pathway and the oxidation-reduction process were significantly enriched. In addition, 15 pairs of miRNA-mRNA interactions that played a role in T. harzianum-induced maize resistance to C. heterostrophus were identified; miR390, miR169j, miR408b, miR395a/p and the novel miRNA miRn5231 were most involved in the induction of maize resistance [176].
Conclusion
In this review, we have presented examples of how plants use epigenetic mechanisms to fight fungi, demonstrating how powerful and necessary this weapon is. All these mechanisms – DNA methylation, covalent histone modifications, ATP-dependent chromatin remodeling and ncRNAs – can operate independently or be coupled. Epigenetic mechanisms are activated in plants not only during pathogenic infection but also during the symbiotic relationship of plants with endophytic fungi, which aims to maintain beneficial interactions and promote symbiosis. Furthermore, epigenetic mechanisms are integral to both plant priming and the formation of stress-induced environmental epigenetic memory.
Data availability
All data generated or analyzed during this study are included in this published article.
Abbreviations
- A:
-
Adenine
- AGO:
-
Argonaute protein
- AMF:
-
Arbuscular mycorrhizal fungi
- AtELP:
-
Elongator Subunit
- BABA:
-
β-aminobutyric acid
- Bgt:
-
Blumeria graminis f. sp. tritici
- BIT:
-
1,2-benzisothiazol-3(2 h)-one,1, 1-dioxide
- BRM:
-
Brahma protein
- BTH:
-
Benzothiadiazole
- CC:
-
Coiled-coil
- C:
-
Cytosine
- CHD:
-
Chromodomain helicase DNA-binding
- CHR:
-
Chromatin remodeler
- COI1:
-
Coronatine insensitive 1
- CRC:
-
Chromosome remodeling complex
- CWDE:
-
Cell wall degrading enzyme
- DAMP:
-
Damage-associated molecular pattern
- DDM1:
-
Decreased DNA Methylation 1
- DME:
-
DEMETER (5-methylcytosine DNA glycosylase/demethylase)
- DML2:
-
DEMETER-Like protein 2
- DMR:
-
Differentially methylated region
- dsRNAs:
-
Double-stranded RNAs
- EIN2:
-
Ethylene-insensitive 2 gene
- ER:
-
ERECTA kinase
- ET:
-
Ethylene
- ETI:
-
Effector-triggered immunity
- G:
-
Guanine
- HAT:
-
Histone acetyltransferase
- hc-siRNAs:
-
Heterochromatic small interfering RNAs
- HDAC:
-
Histone deacetylase
- HR:
-
Hypersensitive response
- INO80:
-
Inositol auxotrophy 80
- ISR:
-
Induced systemic resistance
- ISWI:
-
Imitation switch
- JA:
-
Jasmonic acid
- lncNAT:
-
Long noncoding natural antisense transcript
- lncRNAs:
-
Long non-coding RNAs
- lncTARs:
-
Long noncoding transcriptionally active regions
- LTR:
-
Long terminal repeat
- MAMP:
-
Microbial associated molecular pattern
- MAPK:
-
Mitogen-activated protein kinases
- MET:
-
Methyltransferase
- miRNA:
-
MicroRNA
- MITE:
-
Tourist-miniature inverted-repeat transposable element
- Mla:
-
Mildew resistance locus
- MOS1:
-
Modifiers of snc1
- NB-LRR:
-
Nucleotide-binding leucine-rich repeat
- NB:
-
Nucleotide-binding
- NF:
-
Nuclear factor
- NLR:
-
Leucine-rich domain protein
- NLRs:
-
NB-LRR-containing R proteins
- NPR1:
-
Non-expressor of PR genes
- nt:
-
Nucleotide
- OGs:
-
Oligogalacturonides
- PAMP:
-
Pathogen-associated molecular pattern
- PG:
-
Polygalacturonase
- phasiRNA:
-
Phased small interfering RNA
- PL:
-
Pectin lyase
- PME:
-
Pectin methylesterase
- Pol:
-
Polymerase
- PR genes:
-
Pathogenesis-related genes
- PRR:
-
Pattern recognition receptor
- PTI:
-
Pattern-triggered immunity
- rar1:
-
For Mla12 resistance1
- RBH:
-
(R)-β-homoserine
- RdDM:
-
RNA-directed DNA methylation
- RLKs:
-
Receptor-like kinases
- rom1:
-
Restoration of Mla resistance1
- ROS:
-
Reactive oxygen species
- RT-PCR:
-
Reverse-transcription polymerase chain reaction
- SA:
-
Salicylic acid
- SAR:
-
Systemic acquired resistance
- siRNAs:
-
Small interfering RNAs
- SNC1:
-
Suppressor of npr1-1, constitutive1; one of the R genes
- SOD:
-
Superoxide dismutase
- sRNAs:
-
Small RNAs
- SWI/SNF:
-
Switch/Sucrose on-fermentable
- SYD:
-
Splayed protein
- TE:
-
Transposable element
- TF:
-
Transcription factor
- TIR:
-
Toll/Interleukin1 receptor-like
- T:
-
Thymine
- YDD:
-
YODA DOWNSTREAM
References
Keegstra K. Plant cell walls. Plant Physiol. 2010;154(2):483–6.
Choudhary DK, Prakash A, Johri BN. Induced systemic resistance (ISR) in plants: mechanism of action. Indian J Microbiol. 2007;47(4):289–97.
Durrant WE, Dong X. Systemic acquired resistance. Annu Rev Phytopathol. 2004;42:185–209.
Balakireva AV, Zamyatnin AA. Indispensable role of proteases in Plant Innate Immunity. Int J Mol Sci. 2018;19(2):629.
Tsuda K, Somssich IE. Transcriptional networks in plant immunity. New Phytol. 2015;206(3):932–47.
Birkenbihl RP, Liu S, Somssich IE. Transcriptional events defining plant immune responses. Curr Opin Plant Biol. 2017;38:1–9.
Mauch-Mani B, Baccelli I, Luna E, Flors V. Defense Priming: an adaptive part of Induced Resistance. Annu Rev Plant Biol. 2017;68:485–512.
Lionetti V, Cervone F, Bellincampi D. Methyl esterification of pectin plays a role during plant-pathogen interactions and affects plant resistance to diseases. J Plant Physiol. 2012;169(16):1623–30.
Lionetti V, Francocci F, Ferrari S, Volpi C, Bellincampi D, Galletti R, et al. Engineering the cell wall by reducing de-methyl-esterified homogalacturonan improves saccharification of plant tissues for bioconversion. Proc Natl Acad Sci U S A. 2010;107(2):616–21.
Berlin A, Maximenko V, Gilkes N, Saddler J. Optimization of enzyme complexes for lignocellulose hydrolysis. Biotechnol Bioeng. 2007;97(2):287–96.
Ferrari S, Savatin DV, Sicilia F, Gramegna G, Cervone F, Lorenzo GD. Oligogalacturonides: plant damage-associated molecular patterns and regulators of growth and development. Front Plant Sci. 2013;4:49.
Nishad R, Ahmed T, Rahman VJ, Kareem A. Modulation of Plant Defense System in response to Microbial interactions. Front Microbiol. 2020;11:1298.
Boutrot F, Zipfel C. Function, Discovery, and Exploitation of Plant Pattern Recognition receptors for Broad-Spectrum Disease Resistance. Annu Rev Phytopathol. 2017;55:257–86.
Jones JDG, Dangl JL. The plant immune system. Nature. 2006;444(7117):323–9.
Saijo Y, Loo EPI, Yasuda S. Pattern recognition receptors and signaling in plant-microbe interactions. Plant J. 2018;93(4):592–613.
Klemptner RL, Sherwood JS, Tugizimana F, Dubery IA, Piater LA. Ergosterol, an orphan fungal microbe-associated molecular pattern (MAMP). Mol Plant Pathol. 2014;15(7):747–61.
Windram O, Denby KJ. Modelling signaling networks underlying plant defence. Curr Opin Plant Biol. 2015;27:165–71.
Dangl JL, Jones JD. Plant pathogens and integrated defence responses to infection. Nature. 2001;411(6839):826–33.
Spoel SH, Dong X. How do plants achieve immunity? Defence without specialized immune cells. Nat Rev Immunol. 2012;12(2):89–100.
Nguyen QM, Iswanto ABB, Son GH, Kim SH. Recent advances in Effector-Triggered immunity in plants: New Pieces in the puzzle create a different paradigm. Int J Mol Sci. 2021;22(9):4709.
Peng Y, van Wersch R, Zhang Y. Convergent and divergent signaling in PAMP-Triggered immunity and effector-triggered immunity. Mol Plant Microbe Interact. 2018;31(4):403–9.
Zhang J, Zhou JM. Plant Immunity triggered by Microbial Molecular signatures. Mol Plant. 2010;3(5):783–93.
Kang H, Fan T, Wu J, Zhu Y, Shen WH. Histone modification and chromatin remodeling in plant response to pathogens. Front Plant Sci. 2022;13:986940.
Grant M, Brown I, Adams S, Knight M, Ainslie A, Mansfield J. The RPM1 plant disease resistance gene facilitates a rapid and sustained increase in cytosolic calcium that is necessary for the oxidative burst and hypersensitive cell death. Plant J. 2000;23(4):441–50.
Gao X, Chen X, Lin W, Chen S, Lu D, Niu Y, et al. Bifurcation of Arabidopsis NLR immune signaling via Ca2+-dependent protein kinases. PLoS Pathog. 2013;9(1):e1003127.
Caplan JL, Kumar AS, Park E, Padmanabhan MS, Hoban K, Modla S, et al. Chloroplast Stromules function during innate immunity. Dev Cell. 2015;34(1):45–57.
Ramirez-Prado JS, Abulfaraj AA, Rayapuram N, Benhamed M, Hirt H. Plant Immunity: from signaling to Epigenetic Control of Defense. Trends Plant Sci. 2018;23(9):833–44.
Eberharter A, Becker PB. Histone acetylation: a switch between repressive and permissive chromatin. Second in review series on chromatin dynamics. EMBO Rep. 2002;3(3):224–9.
Kanno T, Kanno Y, Siegel RM, Jang MK, Lenardo MJ, Ozato K. Selective recognition of acetylated histones by bromodomain proteins visualized in living cells. Mol Cell. 2004;13(1):33–43.
Crespo-Salvador Ó, Sánchez-Giménez L, López-Galiano MJ, Fernández-Crespo E, Schalschi L, García-Robles I, et al. The Histone Marks Signature in Exonic and Intronic regions is relevant in early response of Tomato genes to Botrytis Cinerea and in miRNA regulation. Plants (Basel). 2020;9(3):300.
Song N, Lin J, Liu X, Liu Z, Liu D, Chu W, et al. Histone acetyltransferase TaHAG1 interacts with TaPLATZ5 to activate TaPAD4 expression and positively contributes to powdery mildew resistance in wheat. New Phytol. 2022;236(2):590–607.
Zhou C, Zhang L, Duan J, Miki B, Wu K. HISTONE DEACETYLASE19 is involved in jasmonic acid and ethylene signaling of pathogen response in Arabidopsis. Plant Cell. 2005;17(4):1196–204.
Choi SM, Song HR, Han SK, Han M, Kim CY, Park J, et al. HDA19 is required for the repression of salicylic acid biosynthesis and salicylic acid-mediated defense responses in Arabidopsis. Plant J. 2012;71(1):135–46.
Wang Y, Hu Q, Wu Z, Wang H, Han S, Jin Y, et al. HISTONE DEACETYLASE 6 represses pathogen defence responses in Arabidopsis thaliana. Plant Cell Environ. 2017;40(12):2972–86.
Zhi P, Kong L, Liu J, Zhang X, Wang X, Li H, et al. Histone deacetylase TaHDT701 functions in TaHDA6-TaHOS15 complex to regulate wheat defense responses to Blumeria Graminis f.sp. Tritici. Int J Mol Sci. 2020;21(7):2640.
Ding B, Bellizzi M, del R, Ning Y, Meyers BC, Wang GL. HDT701, a histone H4 deacetylase, negatively regulates plant innate immunity by modulating histone H4 acetylation of defense-related genes in rice. Plant Cell. 2012;24(9):3783–94.
Huang L, Sun Q, Qin F, Li C, Zhao Y, Zhou DX. Down-regulation of a SILENT INFORMATION REGULATOR2-related histone deacetylase gene, OsSRT1, induces DNA fragmentation and cell death in rice. Plant Physiol. 2007;144(3):1508–19.
Mengel A, Ageeva A, Georgii E, Bernhardt J, Wu K, Durner J, et al. Nitric oxide modulates histone acetylation at stress genes by inhibition of Histone Deacetylases. Plant Physiol. 2017;173(2):1434–52.
DeFraia CT, Zhang X, Mou Z. Elongator subunit 2 is an accelerator of immune responses in Arabidopsis thaliana. Plant J. 2010;64(3):511–23.
Defraia CT, Wang Y, Yao J, Mou Z. Elongator subunit 3 positively regulates plant immunity through its histone acetyltransferase and radical S-adenosylmethionine domains. BMC Plant Biol. 2013;13:102.
Wang C, Ding Y, Yao J, Zhang Y, Sun Y, Colee J, et al. Arabidopsis Elongator subunit 2 positively contributes to resistance to the necrotrophic fungal pathogens Botrytis Cinerea and Alternaria Brassicicola. Plant J. 2015;83(6):1019–33.
Jaskiewicz M, Conrath U, Peterhänsel C. Chromatin modification acts as a memory for systemic acquired resistance in the plant stress response. EMBO Rep. 2011;12(1):50–5.
Ruthenburg AJ, Allis CD, Wysocka J. Methylation of lysine 4 on histone H3: intricacy of writing and reading a single epigenetic mark. Mol Cell. 2007;25(1):15–30.
Alvarez ME, Nota F, Cambiagno DA. Epigenetic control of plant immunity. Mol Plant Pathol. 2010;11(4):563–76.
Palma K, Thorgrimsen S, Malinovsky FG, Fiil BK, Nielsen HB, Brodersen P, et al. Autoimmunity in Arabidopsis acd11 is mediated by epigenetic regulation of an immune receptor. PLoS Pathog. 2010;6(10):e1001137.
Berr A, McCallum EJ, Alioua A, Heintz D, Heitz T, Shen WH. Arabidopsis histone methyltransferase SET DOMAIN GROUP8 mediates induction of the jasmonate/ethylene pathway genes in plant defense response to necrotrophic fungi. Plant Physiol. 2010;154(3):1403–14.
Bvindi C, Lee S, Tang L, Mickelbart MV, Li Y, Mengiste T. Improved pathogen and stress tolerance in tomato mutants of SET domain histone 3 lysine methyltransferases. New Phytol. 2022;235(5):1957–76.
Lee S, Fu F, Xu S, Lee SY, Yun DJ, Mengiste T. Global regulation of plant immunity by Histone Lysine Methyl Transferases. Plant Cell. 2016;28(7):1640–61.
Weake VM, Workman JL. Histone ubiquitination: triggering gene activity. Mol Cell. 2008;29(6):653–63.
Dhawan R, Luo H, Foerster AM, Abuqamar S, Du HN, Briggs SD, et al. HISTONE MONOUBIQUITINATION1 interacts with a subunit of the mediator complex and regulates defense against necrotrophic fungal pathogens in Arabidopsis. Plant Cell. 2009;21(3):1000–19.
Hu M, Pei BL, Zhang LF, Li YZ. Histone H2B monoubiquitination is involved in regulating the dynamics of microtubules during the defense response to Verticillium Dahliae toxins in Arabidopsis. Plant Physiol. 2014;164(4):1857–65.
Zhang Y, Li D, Zhang H, Hong Y, Huang L, Liu S, et al. Tomato histone H2B monoubiquitination enzymes SlHUB1 and SlHUB2 contribute to disease resistance against Botrytis Cinerea through modulating the balance between SA- and JA/ET-mediated signaling pathways. BMC Plant Biol. 2015;15:252.
Chen X, Xu Q, Duan Y, Liu H, Chen X, Huang J, et al. Ustilaginoidea virens modulates lysine 2-hydroxyisobutyrylation in rice flowers during infection. J Integr Plant Biol. 2021;63(10):1801–14.
Clapier CR, Iwasa J, Cairns BR, Peterson CL. Mechanisms of action and regulation of ATP-dependent chromatin-remodelling complexes. Nat Rev Mol Cell Biol. 2017;18(7):407–22.
Walley JW, Rowe HC, Xiao Y, Chehab EW, Kliebenstein DJ, Wagner D, et al. The chromatin remodeler SPLAYED regulates specific stress signaling pathways. PLoS Pathog. 2008;4(12):e1000237.
Bezhani S, Winter C, Hershman S, Wagner JD, Kennedy JF, Kwon CS, et al. Unique, Shared, and redundant roles for the Arabidopsis SWI/SNF chromatin remodeling ATPases BRAHMA and SPLAYED. Plant Cell. 2007;19(2):403–16.
Li X, Jiang Y, Ji Z, Liu Y, Zhang Q. BRHIS1 suppresses rice innate immunity through binding to monoubiquitinated H2A and H2B variants. EMBO Rep. 2015;16(9):1192–202.
March-Díaz R, García-Domínguez M, Lozano-Juste J, León J, Florencio FJ, Reyes JC. Histone H2A.Z and homologues of components of the SWR1 complex are required to control immunity in Arabidopsis. Plant J. 2008;53(3):475–87.
Berriri S, Gangappa SN, Kumar SV. SWR1 chromatin-remodeling complex subunits and H2A.Z have non-overlapping functions in immunity and gene regulation in Arabidopsis. Mol Plant. 2016;9(7):1051–65.
Cai H, Huang Y, Chen F, Liu L, Chai M, Zhang M, et al. ERECTA signaling regulates plant immune responses via chromatin-mediated promotion of WRKY33 binding to target genes. New Phytol. 2021;230(2):737–56.
Kang H, Liu Y, Fan T, Ma J, Wu D, Heitz T, et al. Arabidopsis CHROMATIN REMODELING 19 acts as a transcriptional repressor and contributes to plant pathogen resistance. Plant Cell. 2022;34(3):1100–16.
Ramirez-Prado JS, Piquerez SJM, Bendahmane A, Hirt H, Raynaud C, Benhamed M. Modify the histone to Win the battle: chromatin dynamics in Plant-Pathogen interactions. Front Plant Sci. 2018;9:355.
Li Y, Tessaro MJ, Li X, Zhang Y. Regulation of the expression of plant resistance gene SNC1 by a protein with a conserved BAT2 domain. Plant Physiol. 2010;153(3):1425–34.
Li Y, Dong OX, Johnson K, Li X, Zhang Y. MOS1 epigenetically regulates the expression of plant resistance gene SNC1. Plant Signal Behav. 2011;6(3):434–6.
Deleris A, Halter T, Navarro L. DNA methylation and demethylation in Plant Immunity. Annu Rev Phytopathol. 2016;54:579–603.
Dowen RH, Pelizzola M, Schmitz RJ, Lister R, Dowen JM, Nery JR, et al. Widespread dynamic DNA methylation in response to biotic stress. Proc Natl Acad Sci U S A. 2012;109(32):E2183–2191.
López A, Ramírez V, García-Andrade J, Flors V, Vera P. The RNA silencing enzyme RNA polymerase v is required for plant immunity. PLoS Genet. 2011;7(12):e1002434.
Le TN, Schumann U, Smith NA, Tiwari S, Au PCK, Zhu QH, et al. DNA demethylases target promoter transposable elements to positively regulate stress responsive genes in Arabidopsis. Genome Biol. 2014;15(9):458.
Hirochika H, Okamoto H, Kakutani T. Silencing of retrotransposons in arabidopsis and reactivation by the ddm1 mutation. Plant Cell. 2000;12(3):357–69.
Saze H, Kakutani T. Heritable epigenetic mutation of a transposon-flanked Arabidopsis gene due to lack of the chromatin-remodeling factor DDM1. EMBO J. 2007;26(15):3641–52.
Hayashi K, Yoshida H. Refunctionalization of the ancient rice blast disease resistance gene pit by the recruitment of a retrotransposon as a promoter. Plant J. 2009;57(3):413–25.
Xin Y, Ma B, Zeng Q, He W, Qin M, He N. Dynamic changes in transposable element and gene methylation in mulberry (Morus notabilis) in response to Botrytis Cinerea. Hortic Res. 2021;8(1):154.
Tirnaz S, Merce C, Bayer PE, Severn-Ellis AA, Edwards D, Batley J. Effect of Leptosphaeria maculans infection on promoter DNA methylation of Defence genes in Brassica napus. Agronomy. 2020;10(8):1072.
Zeng W, Huang H, Lin X, Zhu C, Kosami KI, Huang C, et al. Roles of DEMETER in regulating DNA methylation in vegetative tissues and pathogen resistance. J Integr Plant Biol. 2021;63(4):691–706.
Cui N, Chen X, Shi Y, Chi M, Hu J, Lai K, et al. Changes in the epigenome and transcriptome of rice in response to Magnaporthe oryzae infection. Crop J. 2021;9(4):843–53.
Wang S, Yan W, Yang X, Zhang J, Shi Q. Comparative methylome reveals regulatory roles of DNA methylation in melon resistance to Podosphaera Xanthii. Plant Sci. 2021;309:110954.
Sicilia A, Catara V, Scialò E, Lo Piero AR. Fungal infection induces anthocyanin biosynthesis and changes in DNA methylation configuration of blood Orange [Citrus sinensis L. (Osbeck)]. Plants (Basel). 2021;10(2):244.
Geng S, Kong X, Song G, Jia M, Guan J, Wang F, et al. DNA methylation dynamics during the interaction of wheat progenitor Aegilops tauschii with the obligate biotrophic fungus blumeria graminis f. sp. tritici. New Phytol. 2019;221(2):1023–35.
Song L, Fang Y, Chen L, Wang J, Chen X. Role of non-coding RNAs in plant immunity. Plant Commun. 2021;2(3):100180.
Bhogireddy S, Mangrauthia SK, Kumar R, Pandey AK, Singh S, Jain A, et al. Regulatory non-coding RNAs: a new frontier in regulation of plant biology. Funct Integr Genomics. 2021;21(3–4):313–30.
Huang J, Yang M, Zhang X. The function of small RNAs in plant biotic stress response. J Integr Plant Biol. 2016;58(4):312–27.
Navarro L, Dunoyer P, Jay F, Arnold B, Dharmasiri N, Estelle M, et al. A plant miRNA contributes to antibacterial resistance by repressing auxin signaling. Science. 2006;312(5772):436–9.
Pinweha N, Asvarak T, Viboonjun U, Narangajavana J. Involvement of miR160/miR393 and their targets in cassava responses to anthracnose disease. J Plant Physiol. 2015;174:26–35.
Soto-Suárez M, Baldrich P, Weigel D, Rubio-Somoza I, San Segundo B. The Arabidopsis miR396 mediates pathogen-associated molecular pattern-triggered immune responses against fungal pathogens. Sci Rep. 2017;7:44898.
Li Y, Lu YG, Shi Y, Wu L, Xu YJ, Huang F, et al. Multiple rice microRNAs are involved in immunity against the blast fungus Magnaporthe oryzae. Plant Physiol. 2014;164(2):1077–92.
Li Y, Cao XL, Zhu Y, Yang XM, Zhang KN, Xiao ZY, et al. Osa-miR398b boosts H2 O2 production and rice blast disease-resistance via multiple superoxide dismutases. New Phytol. 2019;222(3):1507–22.
Park YJ, Lee HJ, Kwak KJ, Lee K, Hong SW, Kang H. MicroRNA400-guided cleavage of Pentatricopeptide repeat protein mRNAs renders Arabidopsis thaliana more susceptible to pathogenic bacteria and fungi. Plant Cell Physiol. 2014;55(9):1660–8.
Salvador-Guirao R, Baldrich P, Weigel D, Rubio-Somoza I, San Segundo B. The MicroRNA miR773 is involved in the Arabidopsis Immune response to fungal pathogens. MPMI. 2018;31(2):249–59.
Feng H, Zhang Q, Wang Q, Wang X, Liu J, Li M, et al. Target of tae-miR408, a chemocyanin-like protein gene (TaCLP1), plays positive roles in wheat response to high-salinity, heavy cupric stress and stripe rust. Plant Mol Biol. 2013;83(4–5):433–43.
Xu W, Meng Y, Wise RP. Mla- and Rom1-mediated control of microRNA398 and chloroplast copper/zinc superoxide dismutase regulates cell death in response to the barley powdery mildew fungus. New Phytol. 2014;201(4):1396–412.
Campo S, Peris-Peris C, Siré C, Moreno AB, Donaire L, Zytnicki M, et al. Identification of a novel microRNA (miRNA) from rice that targets an alternatively spliced transcript of the Nramp6 (natural resistance-associated macrophage protein 6) gene involved in pathogen resistance. New Phytol. 2013;199(1):212–27.
Li Y, Zhao SL, Li JL, Hu XH, Wang H, Cao XL, et al. Osa-miR169 negatively regulates Rice Immunity against the Blast Fungus Magnaporthe oryzae. Front Plant Sci. 2017;8:2.
Wang Z, Xia Y, Lin S, Wang Y, Guo B, Song X et al. Osa-miR164a targets OsNAC60 and negatively regulates rice immunity against the blast fungus Magnaporthe oryzae. Plant J. 2018.
Zhu QH, Fan L, Liu Y, Xu H, Llewellyn D, Wilson I. miR482 regulation of NBS-LRR defense genes during fungal pathogen infection in cotton. PLoS ONE. 2013;8(12):e84390.
Ouyang S, Park G, Atamian HS, Han CS, Stajich JE, Kaloshian I, et al. MicroRNAs suppress NB domain genes in tomato that confer resistance to Fusarium oxysporum. PLoS Pathog. 2014;10(10):e1004464.
Ma C, Lu Y, Bai S, Zhang W, Duan X, Meng D, et al. Cloning and characterization of miRNAs and their targets, including a novel miRNA-targeted NBS-LRR protein class gene in apple (Golden Delicious). Mol Plant. 2014;7(1):218–30.
Joshi RK, Megha S, Basu U, Rahman MH, Kav NNV. Genome wide identification and functional prediction of long non-coding RNAs responsive to Sclerotinia sclerotiorum infection in Brassica napus. PLoS ONE. 2016;11(7):e0158784.
Dong W, Ren W, Wang X, Mao Y, He Y. MicroRNA319a regulates plant resistance to Sclerotinia stem rot. J Exp Bot. 2021;72(10):3540–53.
Regmi R, Newman TE, Kamphuis LG, Derbyshire MC. Identification of Brassica napus small RNAs responsive to infection by a necrotrophic pathogen. BMC Plant Biol. 2021;21:366.
Dey S, Sarkar A, Chowdhury S, Singh R, Mukherjee A, Ghosh Z, et al. Heightened miR6024-NLR interactions facilitate necrotrophic pathogenesis in tomato. Plant Mol Biol. 2022;109(6):717–39.
Wu F, Qi J, Meng X, Jin W. miR319c acts as a positive regulator of tomato against Botrytis cinerea infection by targeting TCP29. Plant Sci. 2020;300:110610.
Wu F, Xu J, Gao T, Huang D, Jin W. Molecular mechanism of modulating miR482b level in tomato with botrytis cinerea infection. BMC Plant Biol. 2021;21:496.
Mei J, Wu Y, Niu Q, Miao M, Zhang D, Zhao Y, et al. Integrative analysis of expression profiles of mRNA and MicroRNA provides insights of cotton response to Verticillium Dahliae. Int J Mol Sci. 2022;23(9):4702.
Miao Y, Chen K, Deng J, Zhang L, Wang W, Kong J, et al. miR398b negatively regulates cotton immune responses to Verticillium Dahliae via multiple targets. Crop J. 2022;10(4):1026–36.
Hu G, Wang B, Jia P, Wu P, Lu C, Xu Y, et al. The cotton miR530-SAP6 module activated by systemic acquired resistance mediates plant defense against Verticillium Dahliae. Plant Sci. 2023;330:111647.
Hu G, Hao M, Wang L, Liu J, Zhang Z, Tang Y, et al. The cotton miR477-CBP60A Module participates in Plant Defense against Verticillium dahlia. Mol Plant Microbe Interact. 2020;33(4):624–36.
Hu G, Lei Y, Liu J, Hao M, Zhang Z, Tang Y, et al. The ghr-miR164 and GhNAC100 modulate cotton plant resistance against Verticillium dahlia. Plant Sci. 2020;293:110438.
Wang X, Zhan S, Li M, Chen A, Lu Y, Wang P, et al. Identification and functional analysis of miRNAs in N. benthamiana in response to A. longipes infection. Physiol Mol Plant Pathol. 2022;117:101762.
Xu Y, Wang R, Ma P, Cao J, Cao Y, Zhou Z, et al. A novel maize microRNA negatively regulates resistance to Fusarium verticillioides. Mol Plant Pathol. 2022;23(10):1446–60.
Zhou Z, Cao Y, Li T, Wang X, Chen J, He H, et al. MicroRNAs are involved in Maize Immunity against Fusarium verticillioides ear rot. Genomics Proteom Bioinf. 2020;18(3):241–55.
Talesh Sasani S, Soltani M, Mehrabi B, Fereidoun Padasht-Dehkaei R. Expression alteration of candidate Rice MiRNAs in response to Sheath Blight Disease. Iran J Biotechnol. 2020;18(4):e2451.
Liu L, Chen H, Zhu J, Tao L, Wei C. miR319a targeting of CsTCP10 plays an important role in defense against gray blight disease in tea plant (Camellia sinensis). Tree Physiol. 2022;42(7):1450–62.
Wang S, Liu S, Liu L, Li R, Guo R, Xia X, et al. miR477 targets the phenylalanine ammonia-lyase gene and enhances the susceptibility of the tea plant (Camellia sinensis) to disease during pseudopestalotiopsis species infection. Planta. 2020;251(3):59.
Sun T, Zhou Q, Zhou Z, Song Y, Li Y, Wang HB, et al. SQUINT positively regulates resistance to the Pathogen Botrytis Cinerea via miR156-SPL9 Module in Arabidopsis. Plant Cell Physiol. 2022;63(10):1414–32.
Datta R, Mandal K, Boro P, Sultana A, Chattopadhyay S. Glutathione imparts stress tolerance against Alternaria brassicicola infection via miRNA mediated gene regulation. Plant Signal Behav. 2022;17(1):2047352.
Nandakumar M, Malathi P, Sundar AR, Rajadurai CP, Philip M, Viswanathan R. Role of miRNAs in the host-pathogen interaction between sugarcane and Colletotrichum Falcatum, the red rot pathogen. Plant Cell Rep. 2021;40(5):851–70.
Feng T, Zhang ZY, Gao P, Feng ZM, Zuo SM, Ouyang SQ. Suppression of Rice Osa-miR444.2 improves the resistance to Sheath Blight in Rice mediating through the Phytohormone Pathway. Int J Mol Sci. 2023;24(4):3653.
Yu X, Hou Y, Cao L, Zhou T, Wang S, Hu K, et al. MicroRNA candidate miRcand137 in apple is induced by Botryosphaeria dothidea for impairing host defense. Plant Physiol. 2022;189(3):1814–32.
Yu X, Gong H, Cao L, Hou Y, Qu S. MicroRNA397b negatively regulates resistance of Malus hupehensis to Botryosphaeria dothidea by modulating MhLAC7 involved in lignin biosynthesis. Plant Sci. 2020;292:110390.
Qiao L, Zhao H, Jin H, Niu D. Expression of rice siR109944 in Arabidopsis affects plant immunity to multiple fungal pathogens. Plant Signal Behav. 2020;15(4):1744347.
Qiao L, Zheng L, Sheng C, Zhao H, Jin H, Niu D. Rice siR109944 suppresses plant immunity to sheath blight and impacts multiple agronomic traits by affecting auxin homeostasis. Plant J. 2020;102(5):948–64.
Niu D, Zhang X, Song X, Wang Z, Li Y, Qiao L, et al. Deep sequencing uncovers Rice Long siRNAs and its involvement in immunity against Rhizoctonia solani. Phytopathology. 2018;108(1):60–9.
Dukowic-Schulze S, Sundararajan A, Ramaraj T, Kianian S, Pawlowski WP, Mudge J, et al. Novel meiotic miRNAs and indications for a role of PhasiRNAs in meiosis. Front Plant Sci. 2016;7:762.
Fei Q, Xia R, Meyers BC, Phased. Secondary, small interfering RNAs in Posttranscriptional Regulatory Networks. Plant Cell. 2013;25(7):2400–15.
Liu J, Cheng X, Liu D, Xu W, Wise R, Shen QH. The miR9863 family regulates distinct Mla alleles in barley to attenuate NLR receptor-triggered disease resistance and cell-death signaling. PLoS Genet. 2014;10(12):e1004755.
Hunt M, Banerjee S, Surana P, Liu M, Fuerst G, Mathioni S, et al. Small RNA discovery in the interaction between barley and the powdery mildew pathogen. BMC Genomics. 2019;20(1):610.
Weiberg A, Wang M, Bellinger M, Jin H. Small RNAs: a new paradigm in plant-microbe interactions. Annu Rev Phytopathol. 2014;52:495–516.
Deng Y, Zhai K, Xie Z, Yang D, Zhu X, Liu J, et al. Epigenetic regulation of antagonistic receptors confers rice blast resistance with yield balance. Science. 2017;355(6328):962–5.
Mach J. The long-noncoding RNA ELENA1 functions in Plant Immunity[OPEN]. Plant Cell. 2017;29(5):916.
Bhatia G, Upadhyay SK, Upadhyay A, Singh K. Investigation of long non-coding RNAs as regulatory players of grapevine response to powdery and downy mildew infection. BMC Plant Biol. 2021;21(1):265.
Gao C, Sun J, Dong Y, Wang C, Xiao S, Mo L, et al. Comparative transcriptome analysis uncovers regulatory roles of long non-coding RNAs involved in resistance to powdery mildew in melon. BMC Genomics. 2020;21(1):125.
Zhu QH, Stephen S, Taylor J, Helliwell CA, Wang MB. Long noncoding RNAs responsive to Fusarium oxysporum infection in Arabidopsis thaliana. New Phytol. 2014;201(2):574–84.
Xin M, Wang Y, Yao Y, Song N, Hu Z, Qin D, et al. Identification and characterization of wheat long non-protein coding RNAs responsive to powdery mildew infection and heat stress by using microarray analysis and SBS sequencing. BMC Plant Biol. 2011;11:61.
Gai YP, Yuan SS, Zhao YN, Zhao HN, Zhang HL, Ji XL. A novel LncRNA, MuLnc1, Associated with Environmental stress in Mulberry (Morus multicaulis). Front Plant Sci. 2018;9:669.
Bhatia G, Upadhyay SK, Singh K. Vitis vinifera (grapevine) lncRNAs are potential regulators of response to necrotrophic fungus, Botrytis cinerea infection. Physiol Mol Plant Pathol. 2020;112:101553.
Xing Q, Zhang W, Liu M, Li L, Li X, Yan J. Genome-wide identification of long non-coding RNAs responsive to Lasiodiplodia theobromae infection in Grapevine. Evol Bioinform Online. 2019;15:1176934319841362.
Zhou X, Cui J, Cui H, Jiang N, Hou X, Liu S, et al. Identification of lncRNAs and their regulatory relationships with target genes and corresponding miRNAs in melon response to powdery mildew fungi. Gene. 2020;735:144403.
Wang G, Wang X, Zhang Y, Yang J, Li Z, Wu L, et al. Dynamic characteristics and functional analysis provide new insights into long non-coding RNA responsive to Verticillium Dahliae infection in Gossypium hirsutum. BMC Plant Biol. 2021;21(1):68.
Zhang L, Liu J, Cheng J, Sun Q, Zhang Y, Liu J, et al. lncRNA7 and lncRNA2 modulate cell wall defense genes to regulate cotton resistance to Verticillium wilt. Plant Physiol. 2022;189(1):264–84.
Wang LL, Jin JJ, Li LH, Qu SH. Long non-coding RNAs responsive to Blast Fungus infection in Rice. Rice (N Y). 2020;13(1):77.
Zamora-Ballesteros C, Martín-García J, Suárez-Vega A, Diez JJ. Genome-wide identification and characterization of Fusarium circinatum-responsive lncRNAs in Pinus radiata. BMC Genomics. 2022;23(1):194.
Yin H, Zhang X, Zhang B, Luo H, He C. Revealing the dominant long noncoding RNAs responding to the infection with Colletotrichum gloeosporioides in Hevea brasiliensis. Biol Direct. 2019;14(1):7.
Feng S, Fang H, Liu X, Dong Y, Wang Q, Yang KQ. Genome-wide identification and characterization of long non-coding RNAs conferring resistance to Colletotrichum gloeosporioides in walnut (Juglans regia). BMC Genomics. 2021;22(1):15.
Zhou R, Dong Y, Wang C, Liu J, Liang Q, Meng X et al. LncRNA109897-JrCCR4-JrTLP1b forms a positive feedback loop to regulate walnut resistance against anthracnose caused by Colletotrichum gloeosporioides. Hortic Res. 2023;10.
Duan X, Song X, Wang J, Zhou M. Genome-wide identification and characterization of Fusarium graminearum-responsive lncRNAs in Triticum aestivum. Genes (Basel). 2020;11(10):1135.
Nie J, Wang H, Zhang W, Teng X, Yu C, Cai R, et al. Characterization of lncRNAs and mRNAs involved in Powdery Mildew Resistance in Cucumber. Phytopathology. 2021;111(9):1613–24.
Fu X, Li J, Min D, Zhao X, Liu J, Li F, et al. LncRNA4504 involved in methyl jasmonate-induced resistance to Botrytis cinerea in postharvest tomato fruit. Sci Hort. 2022;305:111381.
Cai Q, Qiao L, Wang M, He B, Lin FM, Palmquist J, et al. Plants send small RNAs in extracellular vesicles to fungal pathogen to silence virulence genes. Science. 2018;360(6393):1126–9.
Zhang T, Zhao YL, Zhao JH, Wang S, Jin Y, Chen ZQ, et al. Cotton plants export microRNAs to inhibit virulence gene expression in a fungal pathogen. Nat Plants. 2016;2(10):16153.
Meng X, Jin W, Wu F. Novel tomato miRNA miR1001 initiates cross-species regulation to suppress the conidiospore germination and infection virulence of Botrytis Cinerea in vitro. Gene. 2020;759:145002.
Wang M, Weiberg A, Lin FM, Thomma BPHJ, Huang HD, Jin H. Bidirectional cross-kingdom RNAi and fungal uptake of external RNAs confer plant protection. Nat Plants. 2016;2:16151.
Dalakouras A, Katsaouni A, Avramidou M, Dadami E, Tsiouri O, Vasileiadis S, et al. A beneficial fungal root endophyte triggers systemic RNA silencing and DNA methylation of a host reporter gene. RNA Biol. 2023;20(1):20–30.
Zeilinger S, Gupta VK, Dahms TES, Silva RN, Singh HB, Upadhyay RS, et al. Friends or foes? Emerging insights from fungal interactions with plants. FEMS Microbiol Rev. 2016;40(2):182–207.
Ågren GI, Hyvönen R, Baskaran P, Ectomycorrhiza. Friend or Foe? Ecosyst. 2019;22(7):1561–72.
Hardoim PR, van Overbeek LS, Berg G, Pirttilä AM, Compant S, Campisano A, et al. The Hidden World within plants: ecological and evolutionary considerations for defining functioning of Microbial endophytes. Microbiol Mol Biol Rev. 2015;79(3):293–320.
Lu H, Wei T, Lou H, Shu X, Chen Q. A critical review on communication mechanism within plant-endophytic Fungi interactions to cope with biotic and Abiotic Stresses. J Fungi (Basel). 2021;7(9):719.
Varga S, Soulsbury CD. Arbuscular mycorrhizal fungi change host plant DNA methylation systemically. Plant Biol (Stuttg). 2019;21(2):278–83.
Vigneaud J, Kohler A, Sow MD, Delaunay A, Fauchery L, Guinet F, et al. DNA hypomethylation of the host tree impairs interaction with mutualistic ectomycorrhizal fungus. New Phytol. 2023;238(6):2561–77.
Hubbard M, Germida J, Vujanovic V. Fungal endophyte colonization coincides with altered DNA methylation in drought-stressed wheat seedlings. Can J Plant Sci. 2014;94:223–34.
Forte FP, Malinowska M, Nagy I, Schmid J, Dijkwel P, Hume DE, et al. Methylome changes in Lolium perenne associated with long-term colonisation by the endophytic fungus Epichloë sp. LpTG-3 strain AR37. Front Plant Sci. 2023;14:1258100.
Etemadi M, Gutjahr C, Couzigou JM, Zouine M, Lauressergues D, Timmers A, et al. Auxin perception is required for arbuscule development in arbuscular mycorrhizal symbiosis. Plant Physiol. 2014;166(1):281–92.
Hofferek V, Mendrinna A, Gaude N, Krajinski F, Devers EA. MiR171h restricts root symbioses and shows like its target NSP2 a complex transcriptional regulation in Medicago truncatula. BMC Plant Biol. 2014;14:199.
Lauressergues D, Delaux PM, Formey D, Lelandais-Brière C, Fort S, Cottaz S, et al. The microRNA miR171h modulates arbuscular mycorrhizal colonization of Medicago truncatula by targeting NSP2. Plant J. 2012;72(3):512–22.
Ariel F, Brault-Hernandez M, Laffont C, Huault E, Brault M, Plet J, et al. Two direct targets of cytokinin signaling regulate symbiotic nodulation in Medicago truncatula. Plant Cell. 2012;24(9):3838–52.
Combier JP, Frugier F, de Billy F, Boualem A, El-Yahyaoui F, Moreau S, et al. MtHAP2-1 is a key transcriptional regulator of symbiotic nodule development regulated by microRNA169 in Medicago truncatula. Genes Dev. 2006;20(22):3084–8.
Formey D, Martín-Rodríguez JÁ, Leija A, Santana O, Quinto C, Cárdenas L, et al. Regulation of small RNAs and corresponding targets in nod factor-Induced Phaseolus vulgaris Root Hair cells. Int J Mol Sci. 2016;17(6):887.
Han G, Cheng C, Zheng Y, Wang X, Xu Y, Wang W, et al. Identification of long non-coding RNAs and the Regulatory Network Responsive to Arbuscular Mycorrhizal Fungi Colonization in Maize roots. Int J Mol Sci. 2019;20(18):4491.
Shah J, Zeier J. Long-distance communication and signal amplification in systemic acquired resistance. Front Plant Sci. 2013;4:30.
Lämke J, Bäurle I. Epigenetic and chromatin-based mechanisms in environmental stress adaptation and stress memory in plants. Genome Biol. 2017;18(1):124.
Alonso C, Ramos-Cruz D, Becker C. The role of plant epigenetics in biotic interactions. New Phytol. 2019;221(2):731–7.
Pieterse CMJ, Zamioudis C, Berendsen RL, Weller DM, Van Wees SCM, Bakker PAHM. Induced systemic resistance by beneficial microbes. Annu Rev Phytopathol. 2014;52:347–75.
Buswell W, Schwarzenbacher RE, Luna E, Sellwood M, Chen B, Flors V, et al. Chemical priming of immunity without costs to plant growth. New Phytol. 2018;218(3):1205–16.
Catoni M, Alvarez-Venegas R, Worrall D, Holroyd G, Barraza A, Luna E, et al. Long-lasting Defence Priming by β-Aminobutyric acid in Tomato is marked by genome-wide changes in DNA methylation. Front Plant Sci. 2022;13:836326.
Yang Z, Zhi P, Chang C. Priming seeds for the future: plant immune memory and application in crop protection. Front Plant Sci. 2022;13:961840.
Wojtasik W, Boba A, Preisner M, Kostyn K, Szopa J, Kulma A. DNA methylation Profile of β-1,3-Glucanase and chitinase genes in Flax shows specificity towards Fusarium Oxysporum strains differing in pathogenicity. Microorganisms. 2019;7(12):589.
Wang S, Wang X, Chen J. Identification of miRNAs involved in Maize-Induced systemic resistance primed by Trichoderma Harzianum T28 against Cochliobolus Heterostrophus. J Fungi (Basel). 2023;9(2):278.
Acknowledgements
Special thanks to Dr Kamil Kostyn for preparations of figures included in this manuscript.
Funding
This research was funded by the Polish National Science Centre, grant number 2019/33/B/NZ9/00408.
Author information
Authors and Affiliations
Contributions
JM and WW wrote and reviewed the main manuscript text.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Mierziak, J., Wojtasik, W. Epigenetic weapons of plants against fungal pathogens. BMC Plant Biol 24, 175 (2024). https://doi.org/10.1186/s12870-024-04829-8
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12870-024-04829-8