《Plant Diversity》:Chromosome-level genome assembly of
Nekemias grossedentata provides new insights into genome evolution and flavonoid biosynthesis
编辑推荐:
Nekemias grossedentata是一种葡萄科(Vitaceae)野生攀缘物种(染色体数2n = 40),其叶片和茎是兼具药用与食用价值的重要资源,素有“植物黄酮之王”之称。该物种富含约30种黄酮类化合物,其中二氢杨梅素(dihydromyricet
Nekemias grossedentata是一种葡萄科(Vitaceae)野生攀缘物种(染色体数2n = 40),其叶片和茎是兼具药用与食用价值的重要资源,素有“植物黄酮之王”之称。该物种富含约30种黄酮类化合物,其中二氢杨梅素(dihydromyricetin)约占总黄酮含量的70%–80%,是主要活性成分。然而,其黄酮含量异常高的分子机制此前尚未阐明。为填补该物种基因组资源匮乏的研究空白,并推动葡萄科系统发育与核型进化(karyotype evolution)的综合解析,研究人员构建了N. grossedentata的染色体水平参考基因组。研究整合DNBseq、Nanopore和Hi-C测序数据,最终获得大小为608.9 Mb、contig N50为5.17 Mb的基因组组装,并通过细胞学观察确认2n = 40的染色体数;BUSCO评估显示98.6%的完整度,共注释到29,173个蛋白质编码基因,重复序列占55.00%(其中LTR占47.74%)。比较基因组分析表明,N. grossedentata共享核心真双子叶植物祖先的全基因组三倍化事件(WGT),且其核型在所有已分析葡萄科物种中最接近祖先真双子叶核型(AEK)。基因家族进化分析鉴定出1265个扩张和2528个收缩的基因家族,扩张家族显著富集于类黄酮生物合成等代谢通路。进一步整合转录组、代谢组及全基因组变异图谱,共鉴定到63个类黄酮生物合成相关基因,并通过WGCNA筛选出四个与高含量二氢杨梅素显著相关的候选基因。研究为该物种的基因组进化与次生代谢调控提供了重要资源与新见解。
论文解读:Nekemias grossedentata染色体水平基因组组装与类黄酮生物合成研究
研究背景
Nekemias grossedentata(葡萄科)是一种野生攀缘植物,2n = 40,其叶和茎长期作为地方特色草药茶和药食同源资源利用,具有抗高血压、降脂、保肝、抗氧化及抗炎等多重药理活性。该植物被誉为“植物黄酮之王”,已分离出约30种黄酮类化合物,其中二氢杨梅素占总黄酮含量的70%–80%,是核心活性物质。然而,尽管其经济与药用价值突出,学界始终未能厘清其高黄酮含量背后的分子调控机制。在谱系演化方面,Nekemias属是葡萄科内最古老的支系之一,但其基因组信息长期缺失,严重制约了家庭尺度的比较基因组学和进化历史重建。以往虽有基于Hi-C组装的两个N. grossedentata基因组报道,但缺少独立的细胞学验证,且遗传背景代表性有限。进入泛基因组时代后,单一或少数地方居群的基因组已远不足以反映属内的基因组多样性。因此,构建并解析一个经核型验证的高质量染色体水平参考基因组,对于揭示葡萄科核型进化、物种分化及次生代谢产物积累的遗传基础具有重要的研究价值。
研究目的与主要科学问题
本研究旨在通过整合序列系统发育学和核型进化分析,解析葡萄科内部的系统发育关系;同时,在基因组水平上探索N. grossedentata高类黄酮含量的遗传基础,并筛选关键候选基因。具体科学问题包括:N. grossedentata的基因组结构特征与进化地位如何?其核型相较于近缘物种是否更保守?类黄酮生物合成相关基因的遗传变异如何贡献于组织特异性表达及代谢物积累差异?
研究内容与主要结论
研究人员利用DNBseq、Nanopore和Hi-C测序技术,组装获得了N. grossedentata染色体水平基因组,大小为608.9 Mb,contig N50为5.17 Mb,BUSCO完整度为98.6%,并通过显微观察确认其染色体数2n = 40。基因组注释共鉴定29,173个蛋白质编码基因,重复序列占55.00%,其中LTR反转录转座子占47.74%。在进化分析中,研究人员发现N. grossedentata仅经历核心真双子叶植物共有的古老全基因组三倍化(WGT),未发生谱系特异性WGD,并与Vitis vinifera及Muscadinia rotundifolia共享约56.9百万年前的共同祖先加倍事件。核型比较显示,N. grossedentata核型与祖先真双子叶核型(AEK)最为接近,比M. rotundifolia和V. vinifera更保守,其20对染色体起源于两条祖先染色体的融合事件,而M. rotundifolia额外发生了涉及4号和20号染色体的相互交叉重组。基因家族分析共将26,672个基因聚类为14,603个基因家族,其中扩张家族显著富集于类黄酮、苯丙烷及黄酮和黄酮醇生物合成等代谢通路。全基因组变异图谱显示,与V. vinifera和M. rotundifolia相比,N. grossedentata的SV相关基因富集于葡聚糖合成和纤维素代谢等生物学过程。针对黄酮类代谢,研究人员共鉴定63个类黄酮生物合成相关基因,其中24个在N. grossedentata与V. vinifera间存在遗传变异,16个与M. rotundifolia间存在变异;转录组分析显示N. grossedentata中19个基因在三个组织间差异表达,而V. vinifera中仅4个,且无遗传变异,提示两物种在类黄酮代谢调控上可能采取了不同的适应性进化策略。通过整合WGCNA(加权基因共表达网络分析)、遗传变异与差异表达分析,最终筛选出4个同时满足高黄酮积累相关模块、种间遗传变异、结构变异和全组织差异表达的候选基因,它们与N. grossedentata的高二氢杨梅素含量显著相关,是功能验证的优先靶点。
主要技术方法与样本来源(限250字概括)
本研究以Nekemias grossedentata叶和茎为材料,采用Kmerfreq估算基因组大小约594 Mb(基于122.35 Gb Illumina数据)。结合DNBseq、Nanopore(114.98 Gb)和Hi-C数据(约129.7 Gb,218×覆盖度)完成de novo组装与染色体挂载,并通过BUSCO评估组装与注释完整度,经细胞学观察验证染色体数目。基因功能注释整合从头预测、转录组辅助和同源搜索方式,并利用InterPro、KEGG、SwissProt等七个数据库进行功能注释。重复序列采用从头预测与同源比对相结合的方法。比较基因组分析包括基因家族聚类(OrthoFinder等)、系统发育树构建、分化时间估算、核型重排推断(基于AEK)以及全基因组变异(SNP、Indel、SV)谱构建。转录组测序涉及幼叶、老叶和茎三组织,代谢物测量包括总黄酮与二氢杨梅素含量(紫外分光光度法/液相色谱法),WGCNA用于模块-性状关联分析。
研究结果分述
1. 基因组组装与注释:最终组装608.9 Mb,锚定99.99%序列至20条染色体,染色体长度22.25–50.53 Mb。注释到29,173个基因,平均基因长度6398.64 bp,重复序列占55.00%,非编码RNA包括546个tRNA、150个miRNA、395个rRNA和621个snRNA。功能注释覆盖86.39%(25,202)的蛋白编码基因。
2. 比较共线分析与进化框架:N. grossedentata与V. vinifera和M. rotundifolia均展示高共线性,支持共享古老WGT事件,且无谱系特异性WGD。染色体4和20上的两个大的染色体内重排可能是种间分化的关键事件。
3. 基因家族进化与系统发育:N. grossedentata拥有1050个特有基因家族,扩张基因家族富集于氧化应激与渗透应激反应以及黄酮、苯丙烷、黄酮和黄酮醇、玉米素等生物合成通路。系统发育树位置与已有报道一致。
4. 核型进化重构:N. grossedentata核型源于21条AEK染色体融合为20条,其保留了最接近AEK的配置,优于M. rotundifolia和V. vinifera。AEK 2/5易位形成3号和12号染色体,AEK 9/17易位形成1号和10号染色体,AEK 3/21易位形成14号和中间染色体,后者再与AEK 18融合形成20号染色体。M. rotundifolia存在4/20号染色体相互易位,且7号和20号经端-端融合形成V. vinifera的7号染色体。
5. 全基因组变异图谱:N. grossedentata与V. vinifera间共检测997,449个SNP、133,532个Indel和1385个SV;与M. rotundifolia间为914,476个SNP、124,120个Indel和1231个SV。变异在染色体末端更集中,提示染色体臂进化更快。SV相关基因富集于葡聚糖与纤维素代谢,KEGG显示富集于生物合成和遗传信息处理通路。
6. 类黄酮生物合成相关基因分析:鉴定到63个类黄酮合成相关基因,其中24个与V. vinifera存在变异,16个与M. rotundifolia存在变异。转录组分析发现N. grossedentata中19个基因在幼叶、老叶、茎三组织间差异表达比例远高于V. vinifera,且结合遗传变异分析揭示N. grossedentata更依赖基因结构变异影响组织差异表达。代谢组与WGCNA结果显示MEturquoise模块与总黄酮(0.83)及二氢杨梅素(0.87)显著正相关,包含6225个基因,其中15个为已鉴定的黄酮合成基因。最终筛选出4个同时满足模块归属、种间遗传变异、结构变异及全组织差异表达的关键候选基因。
讨论与结论
讨论部分强调,将核型进化分析与序列系统发育学整合,可更精细地解析葡萄科内的分化关系。N. grossedentata核型的保守性表明,该谱系在演化过程中保持了更接近祖先的状态,为理解葡萄科染色体数目变异提供了新视角。同时,与V. vinifera相比,N. grossedentata中黄酮相关基因的遗传变异对表达差异贡献更强,提示二者在次生代谢调控上可能经历了不同的选择压力。结合转录组和代谢组数据筛选得到4个候选基因,为后续功能验证和分子育种提供了明确靶点。
研究结论:本研究成功构建了N. grossedentata染色体水平参考基因组,并证明其核型较V. vinifera与M. rotundifolia更为保守。鉴定出63个类黄酮生物合成相关基因,最终筛选出4个与二氢杨梅素高含量显著相关的候选基因。研究不仅加深了对葡萄科特殊物种的基因组理解,也为植物次生代谢进化及葡萄科遗传潜力挖掘提供了宝贵资源。