Paper · 2026 · The Plant Cell

Cruising with an epigenetic brake: role of DNA methylation in maintaining immune homeostasis in Soybean

Meenu Singla-Rastogi; Pei Qin Ng

20
claims
39
verbatim quotes
0
citations
3
research areas
✓
Open access · cc-by

Abstract

Plants constantly walk a fine line between growing and defending themselves. Strong immune responses are essential for survival, but they come at a cost. When defense pathways are switched on too early or for too long, it results in severe penalties to plant growth and overall fitness (Deleris et al. 2016). Epigenetic regulation mechanisms, such as DNA methylation, have long been thought to help manage this growth–defense tradeoff, yet clear links between DNA methylation and immune gene control have been surprisingly hard to pin down. This is especially true in Arabidopsis thaliana, as changes in DNA methylation often prompt limited effects on gene expression owing to the relatively low transposable element (TE) density in its genome (Dowen et al. 2012). In plants, RNA-directed DNA methylation (RdDM) is a key epigenetic regulatory pathway that typically relies on RNA Polymerase (Pol) IV-derived 24-nt small RNA (sRNA) to direct DNA methylation toward regions such as TEs, while RdDM is also achievable using Pol V-derived siRNAs via transcription of noncoding RNAs (Gallego-Bartolomé 2020), likely indirectly. In Glycine max (soybean), the function of these 2 plant-specific polymerases in driving RdDM has yet to be explored together in an RNA-silencing context. Here we present 2 independent studies recently published in The Plant Cell investigating Pol IV- and Pol V-derived sRNA-mediated RdDM in soybean, a crop with a large, transposon-rich genome composed of more than 50% TEs (Schmutz et al. 2010; Wang et al. 2023). In their new work, Zhang and colleagues (2026) used CRISPR/Cas9 genome editing to generate knockout mutants of GmNRPD1 (Pol IV largest subunit) and GmNRPE1 (Pol V largest subunit) to uncover a central role for RdDM in balancing soybean immune responses. In a separate study, Wang and colleagues (2026) generated RNAi lines targeting GmNRPD1, GmNRPE1, and the shared second-largest subunit GmNRP(D/E)2. Together, these complementary genetic approaches reveal a central role for RdDM in maintaining immune homeostasis in soybean. In both studies, significant developmental phenotypes were noted in respective lines, including pronounced developmental phenotypes upon disruption of Pol IV or Pol V. The CRISPR null mutants of GmNRPD1 (gmnrpd1-1, gmnrpd1-2) and GmNRPE1 (gmnrpe1-1, gmnrpe1-2) displayed severe dwarfism and growth retardation. Similarly, RNAi lines targeting GmNRPE1 or GmNRP(D/E)2 exhibited delayed growth under long-day conditions, and GmNRPE1-RNAi plants showed reduced seed set and smaller seeds, while notably GmNRPD1 RNAi lines did not show developmental phenotypes. These phenotypes underscore the importance of Pol IV and Pol V in soybean growth and development. Next, to assess how Pol IV and Pol V affect sRNA production, Zhang et al. (2026) performed sRNA-seq using the null mutants of GmNRPD1 and GmNRPE1. As expected, loss of GmNRPD1 caused a dramatic depletion of 24-nt siRNAs, while GmNRPE1 mutants showed a milder decrease, consistent with their respective roles in the RdDM pathway (Matzke and Mosher 2014). Interestingly, both mutants also accumulated 21- and 22-nt sRNAs, suggesting compensatory or noncanonical sRNA responses at RdDM target loci (Fig. 1a). Effects of pol IV and V knockout and knockdown on sRNA accumulation in soybean mutants. a) Identification of hyper- and hypo-DSRs for each size class (21 to 24 nt) in CRISPR mutants of GmNRPD1 and GmNRPE1 relative to wild type (WT). The genome was divided into non-overlapping 100-base windows, and rRNA-normalized sRNA abundance within each window was compared between mutants and WT. A fold change (FC) ≥ 1.5, and a false discovery rate (FDR) < 0.05 was used as thresholds. Panel from Fig. 2C of Zhang et al. 2026. b) Identification of DSRs for each size class between transgenic RNAi (T) and nontransgenic (NT) lines for GmNRPD1, GmNRP(D/E)2, and GmNRPE1. The genome was tiled into 100-bp non-overlapping windows, and rRNA-normalized sRNA abundance was compared between the T and NT samples for each window. A fold change (FC) ≥ 2 and adjusted P-value <0.01 was used as thresholds. Panel from Fig. 2A of Wang et al., 2026. Consistently, Wang et al. (2026) reported similar observations with sRNA-seq analysis of their RNAi lines. While 24-nt siRNA levels were globally reduced in GmNRPD1 and GmNRP(D/E)2 RNAi lines, the GmNRPE1 RNAi lines were not affected (Fig. 1b). After further categorization of hypo-differential small RNA regions (DSRs) into Pol IV-dependent, Pol V-dependent, and Pol IV-V shared subsets, the authors gained clearer insights into the sRNA accumulation at different DSRs. The authors also reported enriched 21- and 22-nt sRNAs in GmNRP(D/E)2 and GmNRPE1 but not in GmNRPD1 RNAi lines. This was not observed in nrpe1 mutants of other plants such as Arabidopsis, Brasica rapa, and Solanum lycopersicum, suggesting this phenomenon is unique to soybean Pol V (Wang et al. 2026). These 21- to 22-nt GmNRPE1 dependent hyper-DSRs corresponded to TE (213 DSRs) and coding regions (66 DSRs). Genes associated with the 21- to 22-nt GmNRPE1 dependent hyper-DSRs were over-represented in Gene Ontologies (GOs) for functions in “defense response” and “Leucine-rich-repeats.” In Zhang et al. (2026), the impact of RdDM disruption on the methylome was examined using whole-genome bisulfite sequencing. This analysis revealed widespread loss of non-CG DNA methylation, particularly CHH methylation, at thousands of RdDM-controlled regions in null mutants of GmNRPD1 and GmNRPE1. A substantial subset of the differentially methylated regions was enriched in TEs—especially LTR retrotransposons—located close to genes. Similarly, Wang et al. (2026) also reported a significant reduction of non-CG DNA methylation in GmNRPD1 RNAi lines, with significant reduction at CG, CHG, and CHH sites in GmNRP(D/E)2 and GmNRPE1 RNAi lines at 24-siRNA hypo-DSRs. Taken together, both studies suggest that soybean Pol V is key in DNA methylation of CG, CHG and CHH. In both studies, the physiological relevance of this epigenetic misregulation in the mutants was tested using pathogen infection assays. Both GmNRPD1 and GmNRPE1 null mutants exhibited severe dwarfism and growth retardation, reminiscent of autoimmune phenotypes. Unsurprisingly, both mutants displayed enhanced resistance to the oomycete pathogen Phytophthora sojae. Thus, RdDM-mediated repression of immune genes is not merely a passive silencing mechanism but also a crucial safeguard that allows plants to grow normally while retaining the capacity for rapid defense activation. Importantly, this study revealed that RdDM is dynamically modulated during pathogen infection. In wild-type soybean, P. sojae challenge led to transcriptional downregulation of core RdDM components, accompanied by loss of CHH methylation at RdDM-targeted defense genes to promote their robust induction. This suggests that soybean primarily employs a passive release of RdDM-mediated repression, rather than active DNA demethylation, to reprogram immune gene expression during infection (Zhang et al. 2026). Wang et al. (2026) similarly observed significant reductions in P. sojae colonization across all RNAi lines, with GmNRP(D/E)2-RNAi showing the strongest effect. Overall, both studies establish RdDM as a key regulator of immune homeostasis in soybean and illustrate how genome architecture amplifies the biological impact of epigenetic regulation. They also collectively demonstrate that RdDM acts as an “epigenetic brake,” which prevents immune signaling pathways from accelerating when they are not needed. These findings open the door to future strategies to manipulate RdDM to enhance disease resistance in crops while minimizing the growth penalties typically associated with constitutive immune activation. Vaucheret and Voinnet (2024) reviewed the current understanding of the biogenesis of diverse sRNAs in plants and their relative functions in different stages of plant growth, including 24-nt and 21- to 22-nt sRNAs featured in this In Brief. Dodds et al. 2024 reviewed the key principles and players of plant immunity, including LRRs, and highlighted the research that leads to these insights. Zhang et al. 2025 investigated the molecular basis of CG maintenance methylation by plant METHYLTRANSFERASE 1 (MET1), illuminating the molecular basis of MET1 autoinhibition and its preference for hemimethylated DNA substrates. None declared. No new data were generated or analyzed in support of this research. The following phenotypic, genotypic, and functional terms are of significance to the work described in this paper: MET1 Gramene: AT5G49160 MET1 Araport: AT5G49160 NRPE1 Gramene: AT2G40030 NRPE1 Araport: AT2G40030

Abstract from the paper’s OpenAlex record; open copy licensed cc-by. Read at the source ↗

The text the model read, with the source record’s inline markup (italics, sub- and superscripts) rendered and its character entities decoded; quote offsets refer to the raw text.

Claims 20

Mechanismauthor interpretationreview summarybasis n/a

Claim · model-written summary

Complementary genetic studies indicate that RdDM maintains soybean immune homeostasis.

Taxon
soybeanGlycine maxe0
Outcome
immune homeostasise0
Qualifier
RdDMe0

Organisms pictured · Representative pictures of the matched species, not figures from the paper.

Glycine max
Glycine max · TaxonPublic domain · Wikimedia Commons ↗

Verbatim

e0

“Together, these complementary genetic approaches reveal a central role for RdDM in maintaining immune homeostasis in soybean.”

✓ located verbatim in the source · abstract, characters 1940–2065 · show in text ↑
Observationobservedreview summarymeasured

Claim · model-written summary

GmNRPD1 and GmNRPE1 null mutants had severe dwarfism and growth retardation.

Taxon
SoybeanGlycine maxt0
Gene
GmNRPD1e0
Outcome
severe dwarfism and growth retardatione0

Organisms pictured · Representative pictures of the matched species, not figures from the paper.

Glycine max
Glycine max · TaxonPublic domain · Wikimedia Commons ↗

Verbatim

e0

“The CRISPR null mutants of GmNRPD1 (gmnrpd1-1, gmnrpd1-2) and GmNRPE1 (gmnrpe1-1, gmnrpe1-2) displayed severe dwarfism and growth retardation.”

✓ located verbatim in the source · abstract, characters 2234–2376 · show in text ↑
t0

“Soybean”

✓ located verbatim in the source · title, characters 96–103 · show in text ↑
Observationobservedreview summarymeasured

Claim · model-written summary

GmNRPE1 RNAi lines showed delayed growth under long-day conditions.

Taxon
SoybeanGlycine maxt0
Gene
GmNRPE1e0
Outcome
delayed growthe0
Qualifier
long-day conditionse0

Organisms pictured · Representative pictures of the matched species, not figures from the paper.

Glycine max
Glycine max · TaxonPublic domain · Wikimedia Commons ↗

Verbatim

e0

“Similarly, RNAi lines targeting GmNRPE1 or GmNRP(D/E)2 exhibited delayed growth under long-day conditions”

✓ located verbatim in the source · abstract, characters 2377–2482 · show in text ↑
t0

“Soybean”

✓ located verbatim in the source · title, characters 96–103 · show in text ↑
Observationobservedreview summarymeasured

Claim · model-written summary

GmNRP(D/E)2 RNAi lines showed delayed growth under long-day conditions.

Taxon
SoybeanGlycine maxt0
Gene
GmNRP(D/E)2e0
Outcome
delayed growthe0
Qualifier
long-day conditionse0

Organisms pictured · Representative pictures of the matched species, not figures from the paper.

Glycine max
Glycine max · TaxonPublic domain · Wikimedia Commons ↗

Verbatim

e0

“Similarly, RNAi lines targeting GmNRPE1 or GmNRP(D/E)2 exhibited delayed growth under long-day conditions”

✓ located verbatim in the source · abstract, characters 2377–2482 · show in text ↑
t0

“Soybean”

✓ located verbatim in the source · title, characters 96–103 · show in text ↑
Observationobservedreview summarynegativemeasured

Claim · model-written summary

GmNRPE1 RNAi plants had reduced seed set and smaller seeds.

Taxon
SoybeanGlycine maxt0
Gene
GmNRPE1e0
Outcome
seed set and smaller seedse0
Direction
reducede0

Organisms pictured · Representative pictures of the matched species, not figures from the paper.

Glycine max
Glycine max · TaxonPublic domain · Wikimedia Commons ↗

Verbatim

e0

“GmNRPE1-RNAi plants showed reduced seed set and smaller seeds”

✓ located verbatim in the source · abstract, characters 2488–2549 · show in text ↑
t0

“Soybean”

✓ located verbatim in the source · title, characters 96–103 · show in text ↑
Observationobservedreview summaryno effectmeasured

Claim · model-written summary

GmNRPD1 RNAi lines lacked developmental phenotypes.

Taxon
SoybeanGlycine maxt0
Gene
GmNRPD1e0
Outcome
did not show developmental phenotypese0

Organisms pictured · Representative pictures of the matched species, not figures from the paper.

Glycine max
Glycine max · TaxonPublic domain · Wikimedia Commons ↗

Verbatim

e0

“GmNRPD1 RNAi lines did not show developmental phenotypes”

✓ located verbatim in the source · abstract, characters 2565–2621 · show in text ↑
t0

“Soybean”

✓ located verbatim in the source · title, characters 96–103 · show in text ↑
Observationobservedreview summarynegativemeasured

Claim · model-written summary

Loss of GmNRPD1 dramatically depleted 24-nt siRNAs.

Taxon
SoybeanGlycine maxt0
Gene
GmNRPD1e0
Outcome
24-nt siRNAse0
Direction
depletione0

Organisms pictured · Representative pictures of the matched species, not figures from the paper.

Glycine max
Glycine max · TaxonPublic domain · Wikimedia Commons ↗

Verbatim

e0

“loss of GmNRPD1 caused a dramatic depletion of 24-nt siRNAs, while GmNRPE1 mutants showed a milder decrease”

✓ located verbatim in the source · abstract, characters 2881–2988 · show in text ↑
t0

“Soybean”

✓ located verbatim in the source · title, characters 96–103 · show in text ↑
Observationobservedreview summarynegativemeasured

Claim · model-written summary

GmNRPE1 mutants showed a milder decrease in 24-nt siRNAs.

Taxon
SoybeanGlycine maxt0
Gene
GmNRPE1e0
Outcome
24-nt siRNAse0
Direction
decreasee0

Organisms pictured · Representative pictures of the matched species, not figures from the paper.

Glycine max
Glycine max · TaxonPublic domain · Wikimedia Commons ↗

Verbatim

e0

“loss of GmNRPD1 caused a dramatic depletion of 24-nt siRNAs, while GmNRPE1 mutants showed a milder decrease”

✓ located verbatim in the source · abstract, characters 2881–2988 · show in text ↑
t0

“Soybean”

✓ located verbatim in the source · title, characters 96–103 · show in text ↑
Observationobservedreview summarymeasured

Claim · model-written summary

Both polymerase mutants accumulated 21- and 22-nt sRNAs.

Taxon
SoybeanGlycine maxt0
Outcome
21- and 22-nt sRNAse0
Direction
accumulatede0

Organisms pictured · Representative pictures of the matched species, not figures from the paper.

Glycine max
Glycine max · TaxonPublic domain · Wikimedia Commons ↗

Verbatim

e0

“both mutants also accumulated 21- and 22-nt sRNAs”

✓ located verbatim in the source · abstract, characters 3090–3139 · show in text ↑
t0

“Soybean”

✓ located verbatim in the source · title, characters 96–103 · show in text ↑
Observationobservedreview summarymeasured

Claim · model-written summary

GmNRPD1 RNAi lines had globally reduced 24-nt siRNA levels.

Taxon
SoybeanGlycine maxt0
Gene
GmNRPD1e0
Outcome
24-nt siRNA levelse0
Direction
reducede0

Organisms pictured · Representative pictures of the matched species, not figures from the paper.

Glycine max
Glycine max · TaxonPublic domain · Wikimedia Commons ↗

Verbatim

e0

“24-nt siRNA levels were globally reduced in GmNRPD1 and GmNRP(D/E)2 RNAi lines, the GmNRPE1 RNAi lines were not affected”

✓ located verbatim in the source · abstract, characters 4266–4386 · show in text ↑
t0

“Soybean”

✓ located verbatim in the source · title, characters 96–103 · show in text ↑
Observationobservedreview summarymeasured

Claim · model-written summary

GmNRP(D/E)2 RNAi lines had globally reduced 24-nt siRNA levels.

Taxon
SoybeanGlycine maxt0
Gene
GmNRP(D/E)2e0
Outcome
24-nt siRNA levelse0
Direction
reducede0

Organisms pictured · Representative pictures of the matched species, not figures from the paper.

Glycine max
Glycine max · TaxonPublic domain · Wikimedia Commons ↗

Verbatim

e0

“24-nt siRNA levels were globally reduced in GmNRPD1 and GmNRP(D/E)2 RNAi lines, the GmNRPE1 RNAi lines were not affected”

✓ located verbatim in the source · abstract, characters 4266–4386 · show in text ↑
t0

“Soybean”

✓ located verbatim in the source · title, characters 96–103 · show in text ↑
Observationobservedreview summaryno effectmeasured

Claim · model-written summary

GmNRPE1 RNAi lines had unchanged 24-nt siRNA levels.

Taxon
SoybeanGlycine maxt0
Gene
GmNRPE1e0
Outcome
24-nt siRNA levelse0
Direction
not affectede0

Organisms pictured · Representative pictures of the matched species, not figures from the paper.

Glycine max
Glycine max · TaxonPublic domain · Wikimedia Commons ↗

Verbatim

e0

“24-nt siRNA levels were globally reduced in GmNRPD1 and GmNRP(D/E)2 RNAi lines, the GmNRPE1 RNAi lines were not affected”

✓ located verbatim in the source · abstract, characters 4266–4386 · show in text ↑
t0

“Soybean”

✓ located verbatim in the source · title, characters 96–103 · show in text ↑
Observationobservedreview summarymeasured

Claim · model-written summary

GmNRP(D/E)2 RNAi lines were enriched in 21- and 22-nt sRNAs.

Taxon
SoybeanGlycine maxt0
Gene
GmNRP(D/E)2e0
Outcome
21- and 22-nt sRNAse0
Direction
enrichede0

Organisms pictured · Representative pictures of the matched species, not figures from the paper.

Glycine max
Glycine max · TaxonPublic domain · Wikimedia Commons ↗

Verbatim

e0

“The authors also reported enriched 21- and 22-nt sRNAs in GmNRP(D/E)2 and GmNRPE1 but not in GmNRPD1 RNAi lines.”

✓ located verbatim in the source · abstract, characters 4624–4736 · show in text ↑
t0

“Soybean”

✓ located verbatim in the source · title, characters 96–103 · show in text ↑
Observationobservedreview summarymeasured

Claim · model-written summary

GmNRPE1 RNAi lines were enriched in 21- and 22-nt sRNAs.

Taxon
SoybeanGlycine maxt0
Gene
GmNRPE1e0
Outcome
21- and 22-nt sRNAse0
Direction
enrichede0

Organisms pictured · Representative pictures of the matched species, not figures from the paper.

Glycine max
Glycine max · TaxonPublic domain · Wikimedia Commons ↗

Verbatim

e0

“The authors also reported enriched 21- and 22-nt sRNAs in GmNRP(D/E)2 and GmNRPE1 but not in GmNRPD1 RNAi lines.”

✓ located verbatim in the source · abstract, characters 4624–4736 · show in text ↑
t0

“Soybean”

✓ located verbatim in the source · title, characters 96–103 · show in text ↑
Observationobservedreview summaryno effectmeasured

Claim · model-written summary

GmNRPD1 RNAi lines lacked enrichment in 21- and 22-nt sRNAs.

Taxon
SoybeanGlycine maxt0
Gene
GmNRPD1e0
Outcome
21- and 22-nt sRNAse0
Direction
note0

Organisms pictured · Representative pictures of the matched species, not figures from the paper.

Glycine max
Glycine max · TaxonPublic domain · Wikimedia Commons ↗

Verbatim

e0

“The authors also reported enriched 21- and 22-nt sRNAs in GmNRP(D/E)2 and GmNRPE1 but not in GmNRPD1 RNAi lines.”

✓ located verbatim in the source · abstract, characters 4624–4736 · show in text ↑
t0

“Soybean”

✓ located verbatim in the source · title, characters 96–103 · show in text ↑
Observationobservedreview summarymeasured

Claim · model-written summary

GmNRPD1 and GmNRPE1 null mutants lost non-CG methylation at thousands of RdDM-controlled regions.

Taxon
SoybeanGlycine maxt0
Gene
GmNRPD1e0
Outcome
non-CG DNA methylatione0
Direction
losse0
Value
thousands of RdDM-controlled regionse0

Organisms pictured · Representative pictures of the matched species, not figures from the paper.

Glycine max
Glycine max · TaxonPublic domain · Wikimedia Commons ↗

Verbatim

e0

“This analysis revealed widespread loss of non-CG DNA methylation, particularly CHH methylation, at thousands of RdDM-controlled regions in null mutants of GmNRPD1 and GmNRPE1.”

✓ located verbatim in the source · abstract, characters 5339–5514 · show in text ↑
t0

“Soybean”

✓ located verbatim in the source · title, characters 96–103 · show in text ↑
Observationobservedreview summarypositivemeasured

Claim · model-written summary

Both polymerase null mutants had enhanced resistance to Phytophthora sojae.

Taxon
SoybeanGlycine maxt0
Interactor
Phytophthora sojaePhytophthora sojaee0pathogen
Outcome
resistancee0
Direction
enhancede0

Organisms pictured · Representative pictures of the matched species, not figures from the paper.

Glycine max
Glycine max · TaxonPublic domain · Wikimedia Commons ↗

Verbatim

e0

“both mutants displayed enhanced resistance to the oomycete pathogen Phytophthora sojae.”

✓ located verbatim in the source · abstract, characters 6265–6352 · show in text ↑
t0

“Soybean”

✓ located verbatim in the source · title, characters 96–103 · show in text ↑
Observationobservedreview summarymeasured

Claim · model-written summary

P. sojae challenge downregulated core RdDM components in wild-type soybean.

Taxon
SoybeanGlycine maxt0
Interactor
P. sojaePhytophthora sojaee0pathogen
Outcome
core RdDM componentse0
Direction
downregulatione0

Organisms pictured · Representative pictures of the matched species, not figures from the paper.

Glycine max
Glycine max · TaxonPublic domain · Wikimedia Commons ↗

Verbatim

e0

“P. sojae challenge led to transcriptional downregulation of core RdDM components, accompanied by loss of CHH methylation at RdDM-targeted defense genes to promote their robust induction.”

✓ located verbatim in the source · abstract, characters 6685–6871 · show in text ↑
t0

“Soybean”

✓ located verbatim in the source · title, characters 96–103 · show in text ↑
Observationobservedreview summarymeasured

Claim · model-written summary

P. sojae challenge was accompanied by loss of CHH methylation at RdDM-targeted defense genes.

Taxon
SoybeanGlycine maxt0
Interactor
P. sojaePhytophthora sojaee0pathogen
Outcome
CHH methylatione0
Direction
losse0

Organisms pictured · Representative pictures of the matched species, not figures from the paper.

Glycine max
Glycine max · TaxonPublic domain · Wikimedia Commons ↗

Verbatim

e0

“P. sojae challenge led to transcriptional downregulation of core RdDM components, accompanied by loss of CHH methylation at RdDM-targeted defense genes to promote their robust induction.”

✓ located verbatim in the source · abstract, characters 6685–6871 · show in text ↑
t0

“Soybean”

✓ located verbatim in the source · title, characters 96–103 · show in text ↑
Observationobservedreview summarypositivesignificant (as stated)measured

Claim · model-written summary

All RNAi lines reduced P. sojae colonization, with the strongest effect in GmNRP(D/E)2-RNAi.

Taxon
SoybeanGlycine maxt0
Interactor
P. sojaePhytophthora sojaee0pathogen
Gene
GmNRP(D/E)2e0
Outcome
P. sojae colonizatione0
Direction
reductionse0

Organisms pictured · Representative pictures of the matched species, not figures from the paper.

Glycine max
Glycine max · TaxonPublic domain · Wikimedia Commons ↗

Verbatim

e0

“Wang et al. (2026) similarly observed significant reductions in P. sojae colonization across all RNAi lines, with GmNRP(D/E)2-RNAi showing the strongest effect.”

✓ located verbatim in the source · abstract, characters 7076–7236 · show in text ↑
t0

“Soybean”

✓ located verbatim in the source · title, characters 96–103 · show in text ↑