科学问题
作物驯化常伴随剧烈的性状重塑，但传统研究多聚焦于线性DNA序列变异，忽略了更高层级的空间拓扑结构。植物缺乏动物核心架构蛋白CTCF，其三维基因组在群体水平如何演变、受何种规律支配仍是未解之谜。本研究旨在阐明异源四倍体棉花驯化过程中染色质三维折叠如何分化，并探明其对纤维品质性状的调控机制。
背景与领域脉络
真核生物基因组折叠形成区室、拓扑关联结构域与染色质环，精细控制基因转录。在哺乳动物中，CTCF与黏连蛋白介导环挤压并界定结构域边界。植物缺乏CTCF同源物，主要依赖特定转录因子与表观修饰维持构象。前人对植物三维基因组的探索多停留在单一参考株系或特定发育时期，无法捕捉群体层面的构象多态性。由于缺乏群体三维基因组学工具，科学界难以系统评估空间结构变异对作物驯化表型的贡献，这一技术瓶颈长期阻碍了隐性调控靶点的发掘。
方法与技术路线
团队采集34份半野生和267份栽培陆地棉在开花后12天的胚珠纤维，构建了337个原位Hi-C高深度文库。测序平均有效互作率达93.7%，绝大多数达到4 kb分辨率，部分达2 kb。研究结合ATAC-seq与CUT&Tag技术绘制染色质开放性及组蛋白修饰图谱，划分出四类染色质状态。分析中整合重测序与泛基因组结构变异数据，建立染色质互作组全基因组关联分析新框架，并利用序列建模和定点突变对关键基序展开功能验证。
主要结果
绘制高分辨率群体空间图谱
团队获得了301份材料共337个Hi-C数据集，平均有效互作率达93.7%。研究者在98.3%的文库中实现4 kb分辨率，并在每个样本中平均鉴定出9,268个结构域与13,849个染色质环。此外，全群体共汇总出86,468个高置信度染色质环，精准刻画了三维构象的多样性。
结构变异驱动拓扑结构重组
群体中检测到13,100个结构域分裂位点，其中86.60%的结构域保持稳定构象。在30个具有参考基因组的材料中，半野生系结构变异重叠边界的比例为20.18%，显著高于栽培种的6.57%。统计表明，结构变异介导了31.36%的边界破坏事件，有力驱动了构象分化。
空间构象关联纤维优良品质
团队建立的互作组关联分析检测到105个结构域重构与58个染色质环重排，证实它们直接参与调控纤维发育。引入三维基因组特征后，纤维强度的遗传力估算值提升了16%，纤维长度遗传力估算值提高20%。该方法显著弥补了传统全基因组关联分析无法解释的遗传缺失。
驯化选择与保守转录因子
驯化选择扫荡区固定了57%的三维构象特征，表明表观构象选择可与底层序列选择发生解耦。这一过程还改变了39对同源基因的亚基因组表达偏好。序列特征建模与突变实验进一步证实，具有C2H2锌指结构的YY1转录因子在维持植物染色质空间折叠中发挥着保守的调控作用。
创新点
本研究从单基因组静态观察跨越到群体规模动态解析，建立了包含301份材料的棉花全三维基因组资源。研究原创性开发了染色质互作组全基因组关联分析方法，打破了传统关联分析依赖二维线性序列的局限。此外，研究揭示了空间构象选择与DNA序列选择的解耦现象，并突破性发现植物中YY1转录因子具备类似于动物架构蛋白的拓扑组织功能。
意义与影响
该工作证明三维染色质重塑是作物表型演化的独立调控维度，极大拓宽了数量遗传学与演化基因组学的理论视野。所建立的互作组分析框架可广泛迁移至玉米、水稻等大宗作物，显著提高复杂农艺性状的遗传解释度。发现的构象变异为作物精准设计育种提供了大量以往无法识别的非编码区隐性靶点，推动作物分子育种迈向三维空间时代。
局限与待验证
当前研究聚焦于开花后12天单一时间节点的纤维细胞，未能覆盖纤维发育的全周期与其他重要组织。群体中半野生棉样本仅34份，相比267份栽培品种数量较少，可能遗漏部分野生特异性构象。另外，YY1介导染色质折叠的生化互作网络仍需依赖体内蛋白交联与高分辨率显微镜成像进一步确证。
对我们的相关性
对于基因组学研究者，本文展示了群体Hi-C解析复杂非编码变异的范式，所提CIWAS算法可直接集成到现有QTL与GWAS分析流程中。单细胞与空间组学研究者可借鉴其在单一细胞类型中消除组织异质性干扰的策略，探索空间转录组与三维构象的协同演变。AI领域学者则可参考其序列特征建模方法，训练预测非编码突变对染色质折叠及基因表达影响的深度学习模型。
术语速查
全文 PDF 已归档 Google Drive
04-学术文献/Paper/；原文版权归 Cell 及作者所有