IF50.5｜茄属泛基因组揭示旁系同源基因是作物工程中的背景依赖因素

基因复制与旁系同源基因分化是跨物种基因型-表型可预测性的主要障碍。

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

研究背景

泛基因组学和基因编辑技术正在革新全球作物的育种。一个变革性机遇在于将主要作物（全球栽培）与本土作物（局部栽培）之间的基因型-表型知识进行交换，以增强粮食系统。然而，物种特异的遗传变异及其与期望的自然或工程突变之间的相互作用，即使在近缘作物之间，也会阻碍实现可预测的表型效应。茄属（Solanum）包含数十种作物和许多可食用的野生种，是研究旁系同源基因进化的理想系统。此前对番茄泛基因组的研究表明，基因复制是背景修饰因子的重要来源，但旁系同源基因在较短时间尺度上的分化动态及其对性状工程的影响仍不清楚。

科学问题

卡在哪里： 以往对旁系同源基因及其相互作用的研究主要集中在物种内或远缘谱系之间的比较，中间尺度的序列、表达和功能变化在很大程度上未被探索。此外，本土作物如非洲茄（Solanum aethiopicum）的基因组和遗传资源有限，背景修饰因子对性状工程可预测性的影响尚不清楚。

本文要回答： 建立茄属泛基因组，整合功能基因组学和泛遗传学，解析旁系同源基因的进化动态，并评估其对跨物种基因型-表型可预测性的影响。

技术路线

作者主要用了：PacBio HiFi测序、光学图谱、Hi-C、RNA-seq、CRISPR-Cas9、QTL-seq。

作者选取22个茄属物种，涵盖野生、本土食用、观赏和驯化作物，利用PacBio HiFi测序和长程支架技术组装染色体水平参考基因组，构建茄属泛基因组。通过整合多组织RNA-seq数据，注释基因模型并分析旁系同源基因的表达分化模式。进一步以非洲茄为对象，构建其泛基因组并进行QTL-seq定位果实心室数相关位点，结合CRISPR-Cas9基因编辑验证候选基因功能，从而揭示旁系同源基因分化对性状可预测性的影响。

图 1　Fig. 1。图1展示了茄属泛基因组的物种选择、表型多样性和基因组特征。a图显示22个物种的原产地分布，按用途分为野生、本土食用、观赏和驯化。b图展示部分物种的枝条和果实表型多样性。c图的系统发育树将物种分为grade I和clade II两大分支。d图显示基因组大小与重复序列比例的关系，Lasiocarpa亚分支拥有最大的基因组。e图显示基因宏共线性在染色体1、2、6和9上最保守，而clade II内存在大规模染色体重排。

核心亮点

亮点 1｜茄属泛基因组揭示广泛基因复制与动态旁系同源基因

作者组装了22个茄属物种的染色体水平参考基因组，平均QV>53，contig N50平均66.7 Mb，BUSCO完整性96.9%。共注释825,493个基因模型，其中约60%为核心基因。基因复制广泛存在，共鉴定575,464个重复基因，分为全基因组复制（WGD）和单基因复制（串联、邻近、转座、分散）。串联和邻近复制基因的顺式调控序列保守性高，而WGD、分散和转座复制基因在编码序列分歧增大时，顺式调控序列保守性反而增强，提示表达模式可能保持保守。

图 2　Fig. 2。图2揭示了旁系同源基因的表达剂量约束和分化模式。a图示意剂量约束与非约束旁系同源基因对的概念。b图的PCA分析显示5,146个单拷贝基因的表达数据按组织类型聚类。c图显示剂量约束的旁系同源基因对组织特异性较低。d-g图展示四种表达分化模式的定义及其与序列相似性、基因家族大小、表达广度和转录水平的关系。h-i图显示表达保守的旁系同源基因对具有更高的顺式调控序列保守性，且更可能起源于全基因组复制。

亮点 2｜旁系同源基因表达剂量约束与分化模式

通过分析240个RNA-seq样本，作者将旁系同源基因对分为剂量约束和剂量非约束两类。剂量约束的旁系同源基因对表现出较低的非同义突变率和较低的组织特异性。进一步将表达关系分为四组：剂量平衡、旁系同源显性、特化和分化。剂量平衡组（I）具有更高的序列相似性、更小的基因家族、更广的表达范围和更高的转录水平，且主要起源于WGD；而特化和分化组（III和IV）则富集于小规模复制（SSD）。约67%的表达旁系同源基因对属于可能影响作物工程的可变冗余关系。

图 3　Fig. 3。图3展示了茄属中与驯化性状相关的基因变异，特别是CLV3位点的进化动态。a图显示17个与开花时间、花序结构和果实大小相关的直系同源群存在广泛的拷贝数变异和假基因化。b图详细展示了东半球有刺茄分支中CLV3旁系同源基因的单倍型多样性，包括完整和假基因化拷贝的分布。c-e图显示CRISPR-Cas9敲除不同物种CLV3后的表型，单拷贝物种出现严重带化，而S. prinophyllum中单拷贝突变体三心室果实比例增加。

亮点 3｜CLV3旁系同源基因的复杂进化历史

在茄属中，经典果实大小调控基因CLV3的旁系同源基因CLE9在除番茄和Solanum americanum外的所有物种中均假基因化或丢失。然而，在东半球有刺茄分支中，部分物种存在CLV3的局部串联复制，形成完整和假基因化拷贝的复杂单倍型。例如，S. prinophyllum有两个完整拷贝，而S. aethiopicum及其祖先S. anguivi各有一个完整和一个假基因化拷贝。CRISPR-Cas9敲除实验表明，在S. prinophyllum中敲除任一CLV3拷贝均导致三心室果实比例从野生型的5%增加到30%，表明两个拷贝不能完全补偿。

图 4　Fig. 4。图4聚焦于非洲茄的泛基因组和CLV3位点的结构变异。a图显示非洲茄在乌干达和美国的田间种植情况。b图的系统发育树显示Gilo和Aculeatum栽培品种形成两个主要分支，而Shum品种不形成单系群。c图显示超过68%的结构变异位于基因上下游5 kb内，7,234个SV与外显子重叠。d-g图展示一个Aculeatum品种中CLV3位点的约300 kb缺失，导致SaetCLV3a和SaetCLV3b融合为单一功能性拷贝SaetCLV3DEL。

亮点 4｜非洲茄CLV3位点的结构变异与果实心室数

在非洲茄（S. aethiopicum）中，作者发现一个Aculeatum栽培品种（804750136）仅含有一个完整的CLV3拷贝，且存在一个约300 kb的大缺失，将SaetCLV3a和SaetCLV3b融合为一个功能性拷贝SaetCLV3DEL。QTL-seq分析在Gilo×Shum群体中定位到一个位于2号染色体3.9 Mb区域的主效QTL，候选基因为丝氨酸羧肽酶SaetSCPL25-like，其5 bp外显子移码缺失与约两个额外心室相关。CRISPR-Cas9敲除番茄和S. prinophyllum中的同源基因均导致心室数增加，验证了该基因的功能。

图 5　Fig. 5。图5展示了非洲茄心室数的遗传解析。a图显示用于QTL-seq的F2群体构建。b图的QTL-seq分析在2号染色体上鉴定出一个主效QTL（3.9 Mb），在5号和10号染色体上鉴定出微效QTL。c图显示三个位点的等位基因在不同栽培品种中的分布，SaetCLV3DEL仅存在于一个Aculeatum品种中。d图显示CRISPR-Cas9敲除番茄和S. prinophyllum中的SaetSCPL25-like同源基因均增加心室数。e图总结了番茄和非洲茄心室数调控的异同。

亮点 5｜非洲茄心室数调控的独特遗传结构

非洲茄心室数由至少三个位点控制：SaetSCPL25-like的移码突变（主效）、5号染色体上的微效QTL和10号染色体上的SaetCLV3DEL结构变异（微效）。与番茄中CLV3启动子变异导致心室数大幅增加不同，非洲茄中SaetCLV3DEL减少心室数，可能作为上位性修饰因子减弱其他位点的增加效应。这些结果表明，尽管CLV信号通路在番茄和非洲茄中均参与心室数调控，但具体的基因、等位基因和上位性相互作用截然不同。

生信可带走

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

• 方法栈： PacBio HiFi测序、光学图谱、Hi-C、RNA-seq、CRISPR-Cas9、QTL-seq
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 从事比较基因组学、作物遗传改良和基因编辑研究的科研人员。

带走一句

做跨物种基因功能预测时，必须考虑旁系同源基因的剂量约束和分化状态，否则可能误判编辑效果。