IF50.5｜六倍体燕麦泛基因组与泛转录组揭示结构变异对育种的影响

33个燕麦基因组+23个品系6组织表达图谱，揭示染色体倒位与易位如何驱动早花和高产。

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

研究背景

燕麦是世界第七大谷物，富含膳食纤维，近年来因植物基奶制品需求激增而备受关注。燕麦是异源六倍体（A、C、D三个亚基因组），基因组大且重复序列丰富，进化过程中经历了亚基因组间的交换。与小麦、大麦等作物相比，燕麦基因组学研究起步较晚，首个染色体级别参考基因组近两年才发表。此前对燕麦结构基因组多样性和基因表达变异的系统性调查十分缺乏，限制了基因组学辅助育种的应用。

科学问题

卡在哪里： 已有单个参考基因组无法代表燕麦物种的遗传多样性，尤其是大型结构变异（如倒位、易位）在栽培品种中的分布、频率及其对农艺性状的影响尚不清楚。同时，六倍体中基因丢失后同源基因的补偿表达机制及其受亚基因组分化程度的约束也缺乏系统研究。

本文要回答： 构建栽培燕麦及其近缘种的泛基因组和泛转录组，系统解析基因存在/缺失变异、亚基因组表达动态及大型结构变异对育种的影响。

技术路线

作者主要用了：PacBio HiFi组装、Hi-C挂载、RNA-seq、Iso-Seq、kmerGWAS、C显带核型分析。

作者选取33个代表性燕麦品系（包括栽培品种、野生种、合成六倍体等），利用PacBio HiFi长读长测序和Hi-C染色体构象捕获技术组装染色体级别基因组，并整合已发表的短读长组装。对其中23个品系的6个组织（胚、叶、根、茎、穗、颖果）进行RNA-seq和Iso-Seq测序，构建基因表达图谱。通过OrthoFinder构建直系同源框架，划分核心、壳层和云基因，分析基因存在/缺失变异。同时利用全基因组比对和C显带核型分析鉴定大型结构变异，并结合重测序数据和田间试验数据解析其遗传效应。

图 1　Fig. 1。图1展示PanOat品系在9,111份燕麦多样性空间中的位置（MDS图），覆盖了大部分遗传多样性。b图通过单拷贝k-mer曲线估计泛基因组复杂度，显示随着样本增加非冗余序列持续增长。c图展示染色体结构变异概览，包括1A/1C易位、3C倒位等已知事件以及新发现的4C倒位和7D倒位。

核心亮点

亮点 1｜核心与可变基因的功能分化

在30个燕麦品系中，核心基因组包含12,671个直系同源群（占基因总数32.7%），壳层基因占66.2%，云基因占1.1%。核心基因富集于花形成、营养吸收等必需生理过程，且表达水平显著高于壳层和云基因。壳层基因富集于防御反应和种子储存过程，包括MYB、WRKY、NAC等转录因子家族。云基因则富集于磷酸肌醇信号和P型ATP酶。C亚基因组基因的平均表达显著低于A和D亚基因组（Fisher精确检验P=5.46×10⁻⁴⁵），证实了先前在单一品种中的发现。

图 2　Fig. 2。图2a展示6个取样组织。b图显示20个A. sativa品系中稳定和可变基因的比例及表达类别。c图三元图显示颖果组织中稳定（左）和可变（右）三联体的亚基因组表达比例分布，稳定三联体多呈平衡表达，可变三联体分布更分散。d图显示不同组织间表达保守性差异。f图展示基因丢失后剩余同源基因的补偿表达频率，A/D间补偿高于涉及C的补偿。

亮点 2｜同源基因丢失的补偿表达受亚基因组分化约束

在20个A. sativa品系中，分析5,965个1:1:1同源三联体，发现49.4%为稳定三联体，其中94%呈平衡表达。对于缺失一个亚基因组拷贝的944个不完全三联体，A或D拷贝丢失后，剩余D或A拷贝在13.3%和23.2%的病例中显著上调补偿；而C拷贝丢失后，A或D拷贝补偿上调仅发生在4.3%和4%的病例中。A和D亚基因组分化时间约430万年前，远短于与C亚基因组的800万年前，表明进化距离越近，补偿表达越频繁。

图 3　Fig. 3。图3a比较Bilby（祖先型）和GS7（倒位型）的7D染色体，显示约450 Mb的大倒位和50 Mb的远端片段。b图显示307个品系基于SNP的单倍型聚类，绿色为祖先型，黄色为倒位型。c图kmerGWAS曼哈顿图显示7A和7D上的显著关联峰。d图显示不同单倍型的抽穗期差异。e图显示FT1-7D在倒位品系中表达升高。

亮点 3｜7D染色体大倒位与早花关联

在染色体7D上发现一个约420 Mb的着丝粒周围倒位，在PanOat品系和295个北美CORE品种中呈现两种单倍型。携带祖先单倍型（haplotype 1）的品系平均比携带衍生单倍型（haplotype 2）的品系早开花3.7天。kmerGWAS定位到7D上468.7–469.8 Mb区间内的FT1/VRN3同源基因，其在倒位品系中表达显著升高。倒位品系中FT1-7D存在18 bp缺失，FT1-7A存在12 bp缺失导致提前终止密码子。该倒位可能通过抑制重组或改变开花调控基因表达影响抽穗期。

图 4　Fig. 4。图4a比较HiFi（非易位）和Bannister（易位）的2A和2C染色体核型，显示短臂交换。b图展示2A/2C易位的起源和系谱传播，绿色标注携带易位且以OT207为共同亲本的品种。c图显示澳大利亚NVT试验中携带易位品种的产量优势。d图展示ND060432中6D被6A替换的核型证据。e图显示读深分布支持6Al替换6Ds。f图显示6号染色体组中D抑制表达模式显著富集。

亮点 4｜2A/2C易位源自诱变育种并在澳大利亚品种中富集

在澳大利亚品种Bannister中发现2A和2C短臂的相互易位。通过诊断性PCR标记筛查全球564份种质，发现17份（3%）携带该易位，主要来自澳大利亚，包括Dolphin、Echidna、Bannister等高产品种。系谱分析表明这些品种共享亲本OT207，该品系是1970年快中子诱变获得的半矮秆突变体，携带Dw6等位基因。易位与Dw6不连锁，推测为诱变的非预期副作用。在2017-2022年澳大利亚国家品种试验中，携带易位的品种平均产量3,991 kg/ha，显著高于非易位品种的3,667 kg/ha（t检验P=0.03）。

亮点 5｜合成六倍体衍生系中6D被6A替换导致表达偏移

在Amagalon（合成六倍体）衍生的北美育种系中，发现4个品系在6D染色体上读深下降，6A上读深增加。核型分析证实ND060432中6Ds被6Al替换，后者源自A. longiglumis。在叶片组织中，ND060432的6D染色体平均表达2.22 TPM，而6A为14.81 TPM（Wilcoxon检验FDR校正P=2.8×10⁻⁶）。亚基因组表达偏倚分析显示D抑制模式在6号染色体组中显著富集（卡方检验P<0.001，标准化残差32.8）。这些品系因农艺性状优良被选入CORE群体，表明同源染色体替换可被育种利用。

亮点 6｜结构变异影响重组和分离

对13个双亲群体的重组和分离模式分析显示，在1A、1C、3C、4C和7D等染色体区域存在分离畸变和重组抑制，与已鉴定的倒位和易位一致。一个半同胞杂交组合表现出1A和1C的假连锁（与1A/1C相互易位杂合子相关）以及7D上重组抑制（与倒位杂合子相关）。另一个组合在3C上表现出重组抑制。这些结果表明结构变异在育种群体中广泛影响遗传重组，可能导致连锁累赘或有利等位基因的固定。

生信可带走

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

• 方法栈： PacBio HiFi组装、Hi-C挂载、RNA-seq、Iso-Seq、kmerGWAS、C显带核型分析
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 从事作物基因组学、多倍体进化、分子育种或结构变异分析的研究人员。

带走一句

做作物泛基因组时，别只盯着SNP和PAV，大型结构变异（倒位、易位）对重组、表达和农艺性状的影响可能被严重低估，应纳入常规分析流程。