《Plant Biology》：DNA methylome responses to biotic and abiotic stress in Arabidopsis thaliana: A multi-study analysis
编辑推荐：
研究背景：植物会经历多种生物和非生物胁迫，这些胁迫可诱导DNA甲基化发生变化。然而，现有研究之间的比较因分析方法和实验设计的差异而变得复杂。研究人员旨在通过多项研究的整合分析，识别共享的和胁迫特异性的DNA甲基化响应。
方法：研究人员采用统一的生物信息学流程
研究背景：植物会经历多种生物和非生物胁迫，这些胁迫可诱导DNA甲基化发生变化。然而，现有研究之间的比较因分析方法和实验设计的差异而变得复杂。研究人员旨在通过多项研究的整合分析，识别共享的和胁迫特异性的DNA甲基化响应。
方法：研究人员采用统一的生物信息学流程，重新分析了来自13项拟南芥（Arabidopsis thaliana）研究的16个全基因组亚硫酸盐测序（whole-genome bisulphite sequencing，WGBS）数据集。在CG、CHG和CHH序列背景下评估差异甲基化区域（differentially methylated regions，DMRs），并检查其与基因、转座子（transposable elements，TEs）、基因本体（Gene Ontology，GO）术语、基因邻近TE超家族及易发生表观突变位点的关系。
结果：全基因组甲基化水平总体稳定，而DMR的数量和基因组分布随胁迫类型和甲基化背景而变化。CG-DMR主要出现在基因体中，而CHG-和CHH-DMR富集于转座子中。功能分析确定了不同条件下共享的胁迫相关生物学过程。基因邻近的LINE/L1、RathE1_cons和DNA/HAT元件在胁迫响应甲基化变化中显著富集。一部分胁迫相关的CG-DMR与已知在世代间积累稳定表观突变的位点重叠。
结论：标准化的跨研究分析揭示了共享的和胁迫特异性的DNA甲基化模式。胁迫响应DMR与基因邻近转座子及易发生表观突变位点的关联，提示环境响应、基因组调控与长期表观遗传变异之间存在潜在联系。
植物固着生长，需应对干旱、盐碱、极端温度、病原侵染等生物和非生物胁迫。DNA甲基化是植物胁迫适应的重要表观遗传机制，在CG、CHG和CHH序列背景下分别由不同甲基转移酶催化，并参与转座子沉默、基因表达调控和基因组稳定维持。目前已有大量全基因组亚硫酸盐测序研究，但由于分析流程、差异甲基化区域判定阈值、实验设计（胁迫类型、强度、时长、组织、发育时期、基因型等）不同，各研究结果难以直接比较。为此，研究人员从NCBI-SRA中筛选出13项已发表拟南芥研究的16个WGBS数据集，用统一流程重新分析，以鉴定共享及胁迫特异性的甲基化响应，并重点关注胁迫响应的基因、转座子及可能具有长期稳定性的表观变异位点。
主要技术方法包括：使用MethylStar流程统一处理WGBS数据（Trimmomatic修剪、Bismark比对、Methimpute胞嘧啶甲基化调用）；采用jDMR在CG/CHG（500-bp滑动窗口）与CHH（100-bp滑动窗口）中鉴定DMR；基于TAIR10注释将DMR关联至基因、TE、启动子、基因间区和非编码RNA；用g:Profiler进行GO富集分析；计算TE与最近基因距离并按四分位划分，对1.2 kb内基因邻近TE超家族进行置换检验富集分析；按热胁迫、磷饥饿、生物胁迫和叶片组织等元数据特征进行分层亚组分析；利用AlphaBeta分析了野生型Col-0突变积累（mutation accumulation，MA）系17代WGBS数据（GEO GSE223810），估计CG甲基化获得和丢失速率，以评估胁迫相关CG-DMR与自发稳定表观突变位点的重叠。
研究结果如下：
全球DNA甲基化水平在胁迫下保持稳定。比较处理与对照的全基因组CG、CHG、CHH甲基化水平，除干旱在CHG背景下显著降低（P=0.003）外，其余研究均无显著差异，说明全球甲基化稳态受到维持，但局部DMR可能被全基因组平均所掩盖。
胁迫条件相关的DMRs。各研究DMR总数差异较大，从52个到约13.5千个不等。DMR高甲基化比例因胁迫类型和序列背景而异：CG背景下CO-CHL、HE-1、HE-1Rec、PST-S、PSY、TRV、UV偏向高甲基化，GR、DR、HE-2、PST-R偏向低甲基化；CHG背景下PSY高甲基化比例接近90%，DR则明显低甲基化；CHH背景下GR、HE-1、HE-1Rec、HE-2、OS、PST-R、PSY、UV显著富集高甲基化，而CO-CHL、CO-LL、DR、HE-3、SA显著富集低甲基化。DMR密度分析显示，CHG和CHH-DMR主要集中在着丝粒周围区域，CG-DMR在染色体臂上分布更均匀。
基于注释的DMR分析。CG-DMR主要与基因相关，平均70.5%的CG-DMR落在基因区域；CHG-和CHH-DMR则分别有33.9%和42.2%与TE相关，其次为基因间区。CHH-DMR的注释谱在研究间变异最大，但总体上序列背景决定基因组分布的规律较为保守。
基于基因的DMR分析。共鉴定到10,639个非冗余CG-DMR相关基因，其中59.4%至少在两个研究中出现，AT1G67500（REV3）在10个研究中重叠。CHG和CHH背景下也有大量共享基因；CHH背景下AT2G07240在所有16个研究中均出现。Bootstrap检验表明DMR相关基因重叠显著非随机。
功能富集分析。CG背景下共得到474个非冗余GO terms，其中258个至少在两个研究中出现，7个GO terms在15个研究中富集，主要涉及阴离子结合和ATP结合；CHG和CHH背景下富集terms较少，最常涉及细胞壁结构和胞外区。非基因体甲基化（non-gbM）基因比基因体甲基化（gbM）基因在CHG和CHH背景下富集更广泛，主要与细胞壁结构、胞外区和水解酶活性相关。网络分析显示UV、TRV、PST-R、PST-S和HE-1共享较多结合与酶活性相关terms，而SA和DR富集有限。
基于转座子的DMR分析。CG、CHG和CHH-DMR相关TE分别有7,051、10,179和19,504个；CHH-TE共享程度最高，其中8个CHH-TE（多为RC/Helitron）出现于全部16个研究。TE重叠的Bootstrap检验同样显著非随机。
基因-TE邻近DMR分析及邻近TE基因的GO分析。25%的TE位于基因上下游约1.2 kb以内，其中RC/Helitron最为丰富。在基因邻近TE中，DNA/HAT、DNA/MuDR和RC/Helitron在CG、CHG和CHH三种背景下均显著富集于DMR；LTR/Copia在CG背景下多个胁迫中富集，LINE/L1仅在特定处理（如HE-2的CHH）中显著。邻近TE基因的GO富集涉及结合、酶活性、运输、繁殖/发育、防御、细胞壁生物合成，以及CG背景下的染色质重塑和CHH背景下的ncRNA调控。
潜在跨代遗传能力分析。AlphaBeta分析显示，胁迫相关CG-DMR中的CG位点可在MA系世代中积累自发表观突变，提示这些区域可能具有长期表观遗传稳定性；但研究人员并未直接证明所分析胁迫数据中的甲基化变化发生了跨代遗传。
元数据分层亚组分析。在按热胁迫、磷饥饿、生物胁迫和叶片组织分组的亚组中，大多数DMR相关基因和TE仍为研究特异，但CHH背景和叶片亚组中可检测到重复出现的特征；亚组共识的基因邻近TE超家族在CG中为DNA/Pogo和DNA/HAT，在CHG中为DNA/HAT、RathE1_cons、DNA/Pogo和RC/Helitron，在CHH中为DNA/HAT、LTR/Gypsy、DNA/MuDR和RathE3_cons。
讨论指出，全球甲基化稳态在胁迫下未被广泛破坏，但局部DMR呈明显的背景依赖性分布；CG-DMR主要位于基因和染色体臂，非CG-DMR富集于异染色质和TE区域。多个高重现行基因（如REV3、LOX2和IAN相关基因）与DNA损伤耐受、茉莉酸信号和免疫响应相关。非gbM基因比gbM基因更倾向于富集于胁迫相关功能类别，细胞壁重塑和糖苷水解酶活性在热胁迫相关数据中突出。TE甲基化，尤其是LTR/Gypsy元件的反复出现，支持TE作为胁迫相关表观调控的重要靶点。基因邻近TE超家族如RC/Helitron、DNA/MuDR和DNA/HAT的鉴定，扩展了已知的可响应胁迫的TE类型。研究还通过AT5TE35120/CNI1位点的例子说明，跨研究重新分析可以恢复已知的胁迫相关甲基化事件，并提示TE邻近甲基化可能与邻基因调控相关。讨论强调，保守性主要体现在基因组注释、TE超家族和功能类别水平，而非单个位点。局限性在于纳入的WGBS数据集在胁迫强度、暴露时间、组织类型和实验设计上存在异质性，无法进行完全匹配的meta分析；未来需结合甲基化组与转录组匹配数据、定向甲基化编辑和多世代胁迫实验，以确定哪些甲基化变化具有直接调控功能并可能跨代持续。
研究结论：标准化跨研究分析揭示了拟南芥胁迫响应中既共享又具胁迫特异性的DNA甲基化模式。环境胁迫不会导致全基因组甲基化紊乱，而是诱导局部、背景依赖性的基因和TE甲基化变化；这些变化与胁迫相关通路、TE相关调控元件以及一部分具有潜在长期稳定性的CG甲基化位点相关联，提示DNA甲基化可能在环境响应和长期表观遗传变异之间发挥桥梁作用。
