IF31.7｜玉米转录因子结合位点遗传变异解释大部分表型遗传力

利用MOA-seq在25个玉米F1中构建叶片泛顺反组，发现bQTL捕获了约72%表型的多数加性遗传方差。

这篇按生信解读号的读法拆：先看背景和问题，再看作者怎么做、做出了哪几条硬发现，最后把数据和代码收走。

研究背景

过去二十年，全基因组关联分析（GWAS）极大推动了玉米等作物复杂性状遗传基础的研究。已有研究估计，非编码变异贡献了植物表型多样性中约50%的加性遗传方差。然而，由于GWAS分辨率有限，难以直接定位因果SNP或INDEL，且无法提供变异功能的分子证据。转录因子（TF）结合位点（顺式元件）的变异被认为是连接基因型与表型的关键环节，但植物中缺乏全基因组、高分辨率的顺式元件功能变异图谱。MOA-seq作为一种MNase定义的顺反组占据分析方法，能够以较高分辨率（<100 bp）在全基因组范围内鉴定假定的TF结合位点，且相比ATAC-seq能捕获更多未被注释的调控区域。

科学问题

卡在哪里： 尽管MOA-seq等技术能够鉴定TF结合位点，但如何在群体水平上系统地将顺式元件上的遗传或表观遗传变异与TF结合占据的差异关联起来，并评估这些变异对表型遗传力的贡献，仍是未解决的问题。此前缺乏在植物中大规模、高分辨率地构建顺反组功能变异图谱的方法。

本文要回答： 本文旨在利用25个玉米F1杂交种，通过单倍型特异性的MOA-seq，构建玉米叶片在正常和干旱条件下的泛顺反组，鉴定与顺式元件占据相关的遗传和表观遗传变异（bQTL），并评估这些变异对复杂性状遗传力的解释程度。

技术路线

作者主要用了：MOA-seq、RNA-seq、全基因组测序、DNA甲基化分析、双亲本参考基因组、局部关联分析。

作者以B73为共同母本，与25个多样性玉米自交系杂交获得F1，利用F1中同一细胞核内两套单倍型共享反式环境的特点，通过MOA-seq定量每个单倍型的TF足迹占据。测序数据比对到拼接的双亲本参考基因组，消除参考基因组偏好。对MOA峰内的SNP（MP）进行等位基因特异性分析，鉴定出等位基因特异性占据位点（AMP）。进一步结合亲本基因型、DNA甲基化数据，对每个MP进行局部关联分析，鉴定与占据变异相关的bQTL。同时进行RNA-seq获取等位基因特异性表达，并整合GWAS、eQTL等公共数据验证bQTL功能。

图 1　Fig. 1。图1展示了单倍型特异性MOA-seq流程图及验证结果。a部分示意从F1细胞核纯化、MOA-seq产生蛋白-DNA足迹，到利用SNP鉴定等位基因特异性占据，再到等位基因特异性mRNA-seq的完整流程。b-d部分显示B73×Mo17 F1中AMP与亲本等位基因偏倚高度一致（约90%方向相同），且与ZmBZR1 ChIP-seq数据一致性超过70%，验证了方法的可靠性。

核心亮点

亮点 1｜构建玉米叶片泛顺反组

对25个F1杂交种进行单倍型特异性MOA-seq，平均每个F1鉴定到约237,000个MOA峰，覆盖约2%的杂交基因组。在所有MOA峰内的SNP（MP）中，平均19.9%表现出显著的等位基因偏倚（AMP），与ChIP-seq检测单个TF等位基因特异性结合的比例（18.3%）接近。AMP与已发表的ZmBZR1等位基因特异性ChIP-seq数据高度一致（>70%方向相同），验证了MOA-seq的准确性。通过累积分析表明，25个F1已接近饱和，捕获了B73叶片顺反组的大部分。

图 2　Fig. 2。图2比较了三种比对策略对等位基因特异性占据分析的影响。仅用B73参考基因组或SNP替换为N的伪基因组均导致明显的参考基因组偏倚（高密度偏向B73等位基因），而双亲本拼接参考基因组策略消除了这种偏倚，呈现对称分布。该图还展示了bQTL在vgt1位点的分布，包括新发现的vgt1-MOA增强子，以及Hi-C数据支持的远程染色质互作。

亮点 2｜bQTL与已知因果调控位点重合

通过局部关联分析，共鉴定出176,613个bQTL（FDR<0.05），其中93,682个仅由基因型变异解释，51,192个仅由DNA甲基化变异解释，31,739个由两者共同解释。bQTL与多个已知因果调控位点重合，包括控制开花时间的vgt1、vgt1-DMR，以及新发现的第三个增强子vgt1-MOA（位于ZmRAP2.7上游>100 kb），Hi-C数据证实其与vgt1及启动子存在物理互作。此外，bQTL还定位到ZmTRE1启动子中的DOF结合基序、ZmSUB11启动子等，且等位基因特异性占据与基因表达水平显著相关。

图 3　Fig. 3。图3展示了bQTL与基因表达的相关性。a部分显示ZmTRE1和ZmSUB11启动子区bQTL的等位基因分布，b部分显示对应F1中单倍型特异性mRNA计数，表明bQTL等位基因与表达水平显著相关。c部分显示ASE基因启动子区AMP富集（比非ASE基因高74%）。d部分显示bQTL与cis-eQTL的连锁不平衡高于背景SNP。e-f部分以ZmPGM1为例，展示bQTL占据与转录水平的关系，以及转座子插入和甲基化对表达的影响。

亮点 3｜DNA甲基化差异预测MOA占据

在F1中，差异甲基化（CG和/或CHG）区域与14.8%的MP重叠，但在AMP中重叠率升至38.1%，在强偏倚AMP（≥85%偏向一个等位基因）中达到51.5%。98.2%的AMP表现出低甲基化等位基因具有更高的MOA覆盖。在共享差异甲基化等位基因的F1中，71.2%也共享单倍型特异性MOA足迹，而共享同等甲基化等位基因的F1中该比例仅为42.9%。这些结果表明DNA甲基化差异在决定玉米TF结合中起重要作用。

图 4　Fig. 4。图4揭示了DNA甲基化与MOA占据的关系。a部分显示差异甲基化区域与MP、AMP及强偏倚AMP的重叠比例逐步升高。b部分显示98.2%的AMP中低甲基化等位基因具有更高的MOA覆盖。c部分显示共享差异甲基化等位基因的F1中，71.2%也共享单倍型特异性MOA足迹，而共享同等甲基化等位基因的F1中该比例仅为42.9%，表明甲基化差异是TF结合占据的重要预测因子。

亮点 4｜bQTL解释大部分表型遗传力

将bQTL与GWAS命中位点（lead SNP±100 bp）进行共定位分析，发现bQTL富集约2倍（1.75倍和2.17倍，取决于GWAS数据集）。利用NAM群体143个性状的加性遗传方差分区分析，发现仅由基因型关联的bQTL（78,398个）在约72%的性状中解释了大部分加性遗传方差，且优于匹配的背景SNP（81.1%性状中bQTL表现更优）。植物高度、叶片大小/形状和抗病性等性状中bQTL解释的遗传方差比例最高，而维生素E相关性状则主要由背景SNP或基因组其他SNP解释。

图 5　Fig. 5。图5展示了bQTL对表型遗传力的贡献。a部分显示bQTL在GWAS命中位点附近富集约2倍，且该富集不依赖于与基因的距离。b部分显示在143个性状中，bQTL在约72%的性状中解释了大部分加性遗传方差，且优于匹配的背景SNP。c部分按性状类别展示bQTL解释遗传方差的比例，植物高度、叶片性状和抗病性等比例最高，而维生素E相关性状最低。

亮点 5｜干旱响应顺反组及ZmTINY验证

比较正常和干旱条件下的MOA-seq，鉴定出124,504个干旱特异bQTL（DS-bQTL）。整合干旱响应GWAS和eQTL，获得高置信度候选基因，包括已知耐旱基因ZmNAC111（其84 bp MITE插入位点附近存在DS-bQTL）和ZmSO。对ZmTINY的启动子进行双荧光素酶实验，发现CML333等位基因启动子在ABA处理下活性显著高于B73等位基因，与MOA占据和mRNA水平一致。对干旱响应AMP的聚类分析显示，83%的AMP在干旱下改变占据水平而非等位基因偏倚方向，表明干旱诱导的TF占据动态主要源于TF丰度变化而非等位基因偏倚改变。

生信可带走

这一块是给做分析的人用的，不是科普点缀。

• 方法栈： MOA-seq、RNA-seq、全基因组测序、DNA甲基化分析、双亲本参考基因组、局部关联分析
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 从事植物调控基因组学、数量遗传学或作物遗传改良的研究人员，尤其是关注非编码变异功能和表型预测的学者。

带走一句

利用F1杂交种和双亲本参考基因组进行单倍型特异性MOA-seq，可高效构建顺反组功能变异图谱，bQTL能解释大部分表型遗传力，为复杂性状遗传解析和精准编辑提供靶点。