多倍体在植物进化和农业中至关重要，其源于全基因组复制（WGD），能催生新表型并增强适应性。但新形成的多倍体（如人工诱导的四倍体）常出现随机的表观遗传和基因表达变化，这些变化能否在长期进化中保留仍是未解之谜。沙生拟南芥（Arabidopsis arenosa，以下简称”拟南芥”）因同时存在天然二倍体、进化四倍体及可人工诱导的新四倍体，可作为理想的研究模型。

瑞士苏黎世联邦理工学院Kirsten Bomblies 团队利用 ATAC-seq 和 RNA-seq 技术，首次对比分析这三类种群的染色质可及性与基因表达，发现WGD 初期引发的染色质差异区域中，70% 与多倍化直接相关，但多数基因表达变化在进化四倍体中“回归” 二倍体水平，揭示多倍体进化中表观遗传与转录组的动态关联以及独立演化特性。

文章题目：Chromatin Accessibility and Gene Expression Vary Between a New and Evolved Autopolyploid of Arabidopsis arenosa

影响因子：11.0

涉及技术：ATAC-seq、RNA-seq

发表期刊：Molecular Biology and Evolution

发表时间：2024.10

文章概览

主要结论

1. 三种拟南芥间染色质可及性的组织特异性聚类与动态变化趋势

通过 ATAC-seq 在叶片和花瓣中分别鉴定出5,046 个和5,674 个差异染色质可及区域（dACRs），约70% 的 dACRs 由 WGD 直接引发，且具有明显的组织特异性。通过 k-means 聚类将 dACRs 分为 4 个组（k-groups），其中68.7%（叶）和 73.3%（花瓣）的 dACRs 在 WGD 初期（新四倍体 vs 二倍体）即发生显著变化，约38%（叶）至 70%（花瓣）的 dACRs 在进化四倍体中出现可及性趋势逆转。

图1 莲座叶和花瓣样本中发现的 dACRs

2. 在进化的四倍体中，WGD 带来的基因表达变化发生逆转

全基因组复制（WGD）后，新四倍体中差异表达的基因（DEGs）中，约60% 在进化四倍体中恢复至二倍体表达水平，表明多数 WGD 引发的基因表达变化在进化四倍体中逆转至接近二倍体水平；但仍然存在约37.1%-42.9% 的 DEGs在进化四倍体中持续或强化这种表达变化，作者认为这些基因可能参与多倍体特有的表型适应过程。接着，该团队探究了表达变化是否存在组织特异性，发现花瓣 DEGs 总数（3,665 个）远高于叶片（1,967 个），且花瓣中73.3% 的 dACRs 在 WGD 初期变化，其中约70% 在进化中逆转，高于叶片的38% 逆转比例，作者认为进化修正效应在花瓣中更显著。

图2 不同细胞型(天然二倍体、进化四倍体及新四体)的转录组差异

3. WGD过程中可及性变化与表达变化间的关系错综复杂

染色质可及性区域通常是TF的结合热点，参与调控细胞类型或组织特异性基因网络，因此，作者通过整合 ATAC-seq 和 RNA-seq 数据集，研究了染色质可及性变化与多倍体背景下基因表达变化的相关程度。结果发现，仅有少数DEGs与dACRs共定位，且超 40% 的 dACR-DEG 对呈现“不一致关联”，表明二者在多数情况下独立变化。接下来作者研究了 dACRs 是否富集了TF 结合基序，试图证明是否存在由特定先驱 TF 介导的 "倍性响应"，结果显示dACRs 富集 24 个已知 TF 结合基序，但这些 TF 的靶基因中仅少数为 DEGs，提示 TF 介导的协同调控有限。即染色质可及性变化对基因表达的直接调控有限，二者可能通过不同机制影响多倍体表型。

图3 dACR and DEG关联

4. 顺式转座元件与三种拟南芥间染色质可及性动态的定量关联

基于已有的发现，该团队又进一步探究了转座元件（TEs）是否通过顺式作用影响 WGD引发的染色质可及性动态变化。结果显示，超 60% 的 dACRs 与 TEs（含重复序列）共定位，其中 Mutator、hAT、PIF Harbinger 等 DNA 转座子（TIR 型）显著富集，提示 TEs 是驱动可及性动态变化的关键靶点。同时，还发现邻近 TEs 越多的 dACRs，其可及性在细胞型间（天然二倍体、进化四倍体及新四倍体）的变异越强，其中长末端重复转座子（LTR 型 Gypsy）和 DNA 转座子（CACTA）相关区域的可及性变化最显著。最终，该团队推断：TEs 富集区域的可及性变化与邻近基因表达无直接关联，TEs可能主要通过改变染色质结构而非基因转录发挥作用。

图4 dACR 的顺式特征与不同细胞型的可及性变化有关

研究总结

本研究利用 ATAC-seq 和 RNA-seq 技术，分析了拟南芥二倍体、新四倍体及进化四倍体的染色质可及性与基因表达差异，发现全基因组复制初期会引发显著的组织特异性染色质可及性变化，且多数变化在长期进化中逐渐逆转至接近二倍体水平。染色质可及性与基因表达的关联性较弱，而转座元件富集区域的可及性变异更为显著，表明转座元件可能通过重塑染色质结构间接影响多倍体进化，而非直接调控基因表达。研究揭示了多倍化过程中表观遗传与转录组的独立演化动态，为理解基因组复制的适应性机制提供了新视角。

文献延伸——给科研人的实用建议

Ø 借鉴本文 “二倍体→新多倍体→进化多倍体” 的对比模型， 开展“跨代际多组学比较” 

例如：结合 ATAC-seq、RNA-seq 等技术，在其他多倍体物种中解析全基因组复制（WGD）的短期扰动与长期适应机制。

Ø 聚焦转座元件（TEs）与表观遗传的互作网络

例如：结合甲基化测序（如ONT 甲基化）和三维基因组技术（如Pore-C），探究 TEs 如何通过 DNA 甲基化、染色质高级结构重塑影响多倍体的表观遗传稳定性。

Ø 挖掘多倍体适应性的核心调控模块

例如：结合 GWAS、QTL 定位等方法，筛选与多倍体优势性状（如抗逆性）直接相关的稳定可及染色质区域，建立 “表观标记 - 基因表达 - 表型” 的因果关系。

参考文献：

Srikant T, Gonzalo A, Bomblies K. Chromatin Accessibility and Gene Expression Vary Between a New and Evolved Autopolyploid of Arabidopsis arenosa. Mol Biol Evol. 2024.