IF31.7｜花生T2T基因组与群体重测序揭示亚基因组不对称演化与关键农艺性状基因

六个端粒到端粒花生基因组揭示转座子与着丝粒不对称演化，群体分析锁定油分与种子大小关键基因。

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

研究背景

花生是全球重要的油料作物，其栽培种为异源四倍体（AABB），由两个二倍体祖先A. duranensis（AA）和A. ipaensis（BB）杂交后加倍形成。此前已发表多个花生基因组，但多数存在缺口和未定位支架，阻碍了结构变异和着丝粒演化研究。花生品种在种子大小、含油量、荚果形状和颜色等方面表现出丰富多样性，但受限于标记数量，许多关键农艺性状的因果基因尚未明确。近年来，端粒到端粒（T2T）基因组组装技术为解析复杂基因组提供了新机遇。

科学问题

卡在哪里： 已发表的花生基因组中，除Yuanza9102外，普遍存在缺口和未定位支架，无法精确鉴定结构变异和着丝粒演化。此外，var. hirsuta和var. fastigiata等品种类型缺乏基因组序列，限制了花生驯化和品种分化的深入研究。

本文要回答： 构建两个二倍体祖先和四个代表性异源四倍体花生的T2T基因组，结合521份种质重测序，解析基因组结构、演化动态和表型分化的遗传基础。

技术路线

作者主要用了：PacBio HiFi、ONT超长、Hi-C、RNA-seq、代谢组、GWAS、WGCNA。

作者选取两个二倍体祖先和四个代表性四倍体品种，利用PacBio HiFi、ONT超长和Hi-C数据组装T2T基因组，并注释转座子和着丝粒。随后通过比较基因组学鉴定结构变异，结合521份种质的重测序数据进行群体结构、选择信号和GWAS分析，锁定关键候选基因。最后整合转录组和代谢组数据，构建种子发育过程中脂质和花青素代谢的基因调控网络。

图 1　Fig. 1。图1展示六个花生品种的植株、荚果、种子形态及Circos图。Circos图从外到内显示GC含量、基因密度、TE密度和Gypsy/Copia/DNA TE分布。可见不同品种间染色体大小和TE分布存在差异，四倍体基因组明显大于二倍体，且Bt亚基因组TE密度更高。

核心亮点

亮点 1｜At与Bt亚基因组不对称演化

比较两个二倍体和四个四倍体基因组发现，Bt亚基因组TE含量（平均77.77%）始终高于At亚基因组（平均75.08%）。Gypsy元件在B(t)基因组中占58.91%，而在A(t)中仅占52.40%。LTR插入时间分析显示，A(t)基因组在约0.20百万年前有一次扩张，而B(t)基因组在0.69和0.27百万年前出现两次扩张峰。着丝粒长度在At亚基因组中平均增加至46.11 Mb，而Bt亚基因组显著缩短至平均25.70 Mb，主要由重复单元数量减少所致。这些结果表明四倍体形成后两个亚基因组经历了不同的演化轨迹。

图 2　Fig. 2。图2展示LTR逆转座子插入时间、Athila系统发育、CentO重复单元系统发育、CRM系统发育和Retand元件在着丝粒区的分布。A(t)和B(t)基因组在Gypsy扩张时间上明显不同，CentO单体按A/B基因组聚类，表明着丝粒重复单元独立演化。

亮点 2｜13.23 Mb插入与var. hirsuta表型相关

在chr. 14上发现一个13.23 Mb的大插入，存在于K30076、S83、HN51和HN873中，但在S245中缺失。群体分析显示该插入在所有161份var. hirsuta材料中存在，但在var. hypogaea中仅33/161份存在。插入区域包含109个基因，KEGG富集于脂质代谢、木质素生物合成和环境胁迫响应。其中chr.14.1804与水稻OsLAC同源，chr.14.1829与WOX3A/B同源，二者在花生侧枝、腋芽和下胚轴中表达，提示该插入可能参与var. hirsuta独特株型的形成。

图 3　Fig. 3。图3展示chr.14上13.23 Mb插入的验证。Illumina和HiFi reads覆盖度在插入区域出现明显差异，Hi-C热图显示S245与S83/K30076在该区域存在染色质互作差异。GO富集显示插入区基因涉及脂质代谢和胁迫响应，水稻同源基因OsLAC和WOX3A/B在花生特定组织中表达。

亮点 3｜AhWRI1启动子变异调控花生含油量

GWAS在chr. 8上定位到与含油量显著关联的位点，候选基因AhWRI1（chr.8.2620）距峰顶约54.9 kb。该基因启动子区存在A到C的碱基变异，第一外显子存在非同义SNP（T到G，精氨酸变为甲硫氨酸）。Hap2型（C变异）材料平均含油量54.10%，显著高于Hap1型（48.41%）。酵母单杂交和双荧光素酶实验表明，C变异破坏了AhFUS3与启动子中RY/E-box基序的结合，导致AhWRI1转录活性升高，进而上调下游脂肪酸合成基因AhACP1和AhKAS1的表达。

图 4　Fig. 4。图4展示521份花生种质的群体结构。系统发育树和ADMIXTURE分析将材料分为五个类群（G1-G5），PCA与之一致。FST和π分析显示G2与其他类群分化最大，G4与G5差异最小。LD衰减速率G1最低、G2最高。fd统计显示G2-G3和G4-G5间存在显著基因流，主要发生在At亚基因组。

亮点 4｜AhGSA1启动子插入缺失影响种子大小

GWAS在chr. 16上鉴定到与种子大小和重量相关的候选基因AhGSA1（chr.16.3093），其启动子区存在ATT/AT插入缺失。Hap1（ATT型）平均千粒重490.82 g，而Hap2（AT型）平均千粒重846.34 g（P=1.00×10^-20）。RT-qPCR显示Hap2型材料中AhGSA1表达量显著高于Hap1型。双荧光素酶和GUS染色实验证实AT启动子活性显著强于ATT启动子。Hap1主要分布在var. fastigiata群体中，表明该变异在品种分化中受到选择。

图 5　Fig. 5。图5展示含油量GWAS结果和AhWRI1功能验证。曼哈顿图显示chr.8上显著关联信号，LD block定位候选区间。AhWRI1基因结构显示启动子和外显子变异。Hap1和Hap2材料含油量存在显著差异。亚细胞定位显示AhWRI1定位于细胞核，过表达拟南芥和烟草叶片中脂滴增多。酵母单杂交和双荧光素酶实验证实AhFUS3结合AhWRI1启动子，C变异削弱结合并增强转录活性。

亮点 5｜转录组与代谢组揭示种子发育代谢网络

对S83和S245五个发育阶段的种子进行转录组和代谢组分析，共鉴定497种脂质代谢物。S83在各阶段脂质含量均高于S245，尤其在后期。WGCNA鉴定出与脂质积累高度相关的基因模块，包含脂肪酸去饱和酶和酮脂酰-ACP合酶等关键酶。花青素分析发现78种代谢物，S83富含矮牵牛素、锦葵素等，而S245富含原花青素、芍药素等，与种皮颜色差异一致。WGCNA鉴定出MYB和bHLH转录因子作为花青素生物合成的关键调控因子。

生信可带走

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

• 方法栈： PacBio HiFi、ONT超长、Hi-C、RNA-seq、代谢组、GWAS、WGCNA
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 从事作物基因组学、群体遗传学或油料作物分子育种的科研人员。

带走一句

T2T基因组结合群体重测序和组学分析，能精准锁定复杂基因组中的关键农艺性状基因，为花生分子育种提供范例。