广东省农业科学院等团队在 Nat. Commun. (Q1/15.7) 发文!揭示茄子果色和青枯病抗性关键结构变异
茄子（Solanum melongena L.）是全球第三大茄科作物，在亚洲、地中海地区广泛种植，具有重要的经济和营养价值。然而，其驯化历史尚存争议，且与重要农艺性状（如果实颜色、青枯病抗性）相关的结构变异（Structural Variants, SVs）一直缺乏系统性解析。长读长测序和泛基因组学的兴起，为解决这一问题提供了强大工具。广东省农业科学院蔬菜研究所联合华南农业大学等单位在顶级期刊 《Nature Communications》 上发表了题为“Sub-pangenome analysis reveals structural variants associated with fruit color and bacterial wilt resistance in eggplant”的研究论文。该研究通过对226份茄子种质的群体基因组分析，结合11个新组装和6个已发表的基因组，构建了茄子亚泛基因组图谱，成功鉴定出与果色形成和青枯病抗性相关的关键结构变异，为茄子分子育种提供了重要基因资源和理论依据。
通讯作者及单位Vegetable Research Institute, Guangdong Academy of Agricultural Sciences, Guangzhou, Chinakuaileyouqianmeng@163.com
Baojuan Sun
Vegetable Research Institute, Guangdong Academy of Agricultural Sciences, Guangzhou, Chinasunbaojuan@hotmail.com
Vegetable Research Institute, Guangdong Academy of Agricultural Sciences, Guangzhou, Chinatianxing84@163.com
College of Horticulture, South China Agricultural University, Guangzhou, Chinayiliao@scau.edu.cn
核心内容核心结果1：茄子种质资源的群体结构与遗传多样性。研究团队对226份茄子种质（主要来自中国和东南亚）进行了全基因组重测序，平均覆盖深度达20.87×。群体结构分析（PCA和Admixture）将这些种质清晰地划分为四个类群（C1-C4），其分布与地理来源高度相关。其中，来自东南亚的C1群和来自南美/欧洲的C2群表现出较高的遗传多样性（π值分别为0.309和0.224），而主要来自中国的C3和C4群多样性相对较低，这可能反映了驯化和传播过程中的选择瓶颈效应。核心结果2：茄子亚泛基因组构建与基因家族分类。为了更全面地捕捉茄子的基因组多样性，研究团队利用PacBio HiFi长读长测序技术，从头组装了11个代表性种质的染色体级别基因组，并与6个已发表的基因组整合，构建了基于17个基因组的茄子亚泛基因组。分析显示，泛基因家族总数随着基因组增加而趋于饱和。其中，核心基因家族（在所有17个基因组中都存在）占41.15%，而可变基因家族（存在于部分基因组）占47.79%。值得注意的是，可变基因家族显著富集在“防御反应”等GO条目中，暗示它们可能参与了茄子对环境适应和抗病性的多样化调控。核心结果3：结构变异（SV）的鉴定、特征及功能影响。通过比较16个茄子基因组与参考基因组‘S076’，研究共鉴定到156,540个SVs（≥50 bp），包括插入、缺失、易位和倒位。其中，近60%的SVs为某个基因型所特有，且大多数SVs在群体中以低频存在。利用外群（S. insanum）推断祖先状态后，发现SVs的衍生等位基因频率显著低于SNPs，表明SVs平均对茄子适应性更具危害性。进一步分析表明，SVs会显著影响附近基因的表达波动。例如，一个418 bp的缺失和10,296 bp的缺失分别与两个基因（编码WAT1相关蛋白和BTB/POZ-MATH蛋白）的表达量显著降低相关。核心结果4：一个12.4 Mb的大片段倒位与茄子果色紧密关联。针对茄子果色（紫色 vs. 非紫色）的全基因组关联分析（GWAS）在10号染色体上检测到一个超显著的信号区间。该区间恰好与一个长度为12.4 Mb的大片段倒位重合。该倒位在群体中的频率高达50.44%，且在111个（96.52%）紫色果实的种质中存在。深入研究发现，该倒位本身并不直接影响已知的花青素合成基因SmMYB1的表达，但它通过其所在区域的共分离效应与一个关键的SmMYB1基因的6-bp缺失（导致功能丧失）形成“搭车效应”（hitchhiking），共同决定了果实的紫色表型。此外，该倒位区域表现出显著降低的核苷酸多样性、升高的Fst值和延长的连锁不平衡（LD）区块，表明它可能在育种过程中受到了选择。核心结果5：青枯病抗性相关基因及拷贝数变异（CNV）的发现。青枯病是影响茄子生产的主要病害之一。通过对197份种质进行青枯菌接种和GWAS分析，研究团队在4号和5号染色体上鉴定到多个显著关联信号。在4号染色体上，SmCYP82D47基因中的一个提前终止密码子（Premature stop codon）与青枯病抗性显著相关。在5号染色体上，GWAS信号富集在一个包含SmEPS1（参与水杨酸合成）和SmRog1（免疫受体）同源基因的区域。通过拷贝数变异（CNV）分析发现，这两个基因在抗病种质中往往具有更高的拷贝数。病毒诱导的基因沉默（VIGS）实验进一步证实，沉默SmCYP82D47、SmEPS1或SmRog1会显著降低茄子的青枯病抗性，揭示它们是抗病通路中的关键组分。
图解图1：226份茄子种质的地理分布与群体结构。a. 种质的地理来源及类群（C1-C4）分布图。b. 茄子果实、花和叶片的多样化表型展示。c-e. 主成分分析（PCA）和Admixture分析揭示了四个与地理来源一致的遗传类群。f. 基于SNPs构建的最大似然进化树。g. 各群体内的核苷酸多样性（π）和群体间的遗传分化指数（Fst）。图2：基于17个基因组的茄子亚泛基因组分析。a. 包含新组装及已发表基因组的230份种质进化树。b. 新组装基因组的连续性评估（BUSCO > 98.7%）。c. 泛基因家族（pan-genes）随基因组数量增加的增长曲线。d-e. 核心（Core）、软核心（Softcore）、可有可无（Dispensable）和私有（Private）基因家族的比例及在各基因组中的分布。f-g. 不同类别基因家族在表达水平及进化选择压力（Ka/Ks）上的差异。图3：茄子基因组中结构变异（SV）的特征。a. 各基因组相对于参考‘S076’的SV数量。b. 随样本增加，SV共享与特有数量的变化。c. 不同类型SV的占比（以TE插入为主）。d. SVs的衍生等位基因频率（dAFS）显著低于SNPs。e-f. SVs在基因区显著耗竭，且不同类型的SV距离基因的远近不同。g. 存在SVs的基因集表现出更高的表达变异性。h-i. 特定SVs与基因表达量变化的关联实例。图4：一个12.4 Mb的大片段倒位与茄子果色显著关联。a-b. GWAS曼哈顿图显示10号染色体的强烈信号与该倒位区域重合。c. 在14个基因组组装中，倒位（Inv）与紫色果实高度相关。d. 对227份种质的倒位基因型分析，显示携带倒位的个体绝大多数为紫色。e. 倒位区域内差异表达基因（DEGs）分析，仅有SmMYB1位于断点区且差异显著。f-g. SmMYB1基因中的一个6-bp缺失在非紫色种质中几乎固定。h. 表达量分析表明SmMYB1表达差异与果色相关，而非倒位本身。图5：12.4 Mb倒位区域伴随的遗传多样性降低与长单倍型区块。a. 该倒位在中国种质（C3, C4）中频率最高。b-c. 携带倒位的群体在倒位区域表现出更低的π值和更高的Fst值。d. Tajima’s D值在倒位区域出现显著分化。e. 倒位区域延伸出一个长单倍型区块和约6 Mb的强连锁不平衡（LD）区块。f. LD衰减分析显示倒位对不同区域LD模式的影响存在显著差异。图6：与青枯病抗性相关的基因组位点。a-b. 197份种质在青枯菌接种后的表型变异范围和动态变化。c. GWAS鉴定出包括Chr04和Chr05在内的多个显著关联峰。d. SmCYP82D47基因中的提前终止密码子（TGA）与感病表型显著相关。e. SmEPS1和SmRog1基因的拷贝数变异（CNV）分析。f. CNV与抗性表型的关联，显示高拷贝数倾向于抗病。g. VIGS实验证实沉默这些基因会降低抗病性。
展望这项研究不仅仅是一个基因组资源的汇编，更是对茄子重要性状遗传基础的一次深入挖掘。首先，所构建的茄子亚泛基因组图谱和超过18万个非冗余SVs数据集，为茄子功能基因组学和分子育种提供了宝贵的公共资源，有效减少了对单一参考基因组的偏倚。其次，12.4 Mb大片段倒位对果色的贡献，以及SmCYP82D47、SmEPS1和SmRog1在抗病性中的协同作用，明确了多个具有育种应用潜力的关键基因与变异。特别是对于青枯病抗性，研究揭示了不同抗性位点间的上位效应（SmCYP82D47是关键开关），并指出通过聚合有益等位基因（gene pyramiding）来培育广谱、持久抗性品种的可行路径。未来，可以将这些关键的SV标记（如倒位断点、CNV）开发成简便、经济的分子标记，直接应用于茄子标记辅助选择（MAS）育种。此外，对SmEPS1和SmRog1等基因家族的CNV如何精确调控抗性信号通路的机制研究，以及这些SVs在其他茄科作物（如番茄、辣椒）中的保守性，将是极具吸引力的研究方向。这项研究为利用结构变异进行作物改良提供了杰出的范例，预示着精准育种新时代的到来。
10个月宝宝每天需要喝多少奶粉？
