《Frontiers in Plant Science》：Large-scale organellar genomics reveals maternal evolutionary history and dynamism in the wolfberries
编辑推荐：
摘要专业翻译
引言： 枸杞属（Lycium）复合群的细胞器基因组多样性与母系演化历史仍未得到充分解析。
方法： 研究人员对201份枸杞属材料的细胞器基因组进行了比较分析，其中对宁夏枸杞（L. barbarum）和黑果枸杞（L. ruthenicum）进行了密集
摘要专业翻译
引言： 枸杞属（Lycium）复合群的细胞器基因组多样性与母系演化历史仍未得到充分解析。
方法： 研究人员对201份枸杞属材料的细胞器基因组进行了比较分析，其中对宁夏枸杞（L. barbarum）和黑果枸杞（L. ruthenicum）进行了密集采样，并纳入了具有代表性的近缘物种。研究人员以de novo方式组装了一个黑果枸杞线粒体基因组（mitogenome），并将其作为186份黑果枸杞与宁夏枸杞复合群材料比较分析的参考序列。
结果： 质体基因组（plastome）分析解析出七个母系遗传聚类（maternal genetic clusters），并鉴定出与红果枸杞复合群复杂母系祖先相符的模式。然而，仅凭质体数据无法区分杂交、基因渐渗（introgression）与不完全谱系分选（incomplete lineage sorting）。全部201个质体基因组均保留了保守的四分体结构（quadripartite structure）和79个独特的蛋白质编码基因。Illumina测序深度图谱鉴定出一个约5.5 kb的候选结构差异区间，该区间在宁夏枸杞代表材料中表现出显著降低的覆盖度。对应的参考区间包含orf27、orf92和orf93三个开放阅读框。该深度支持的区间界定了一个候选线粒体区域，但尚未确定其精确的分子构型或生物学功能。
讨论： 该细胞器基因组资源为检验演化假说和支持枸杞种质资源研究提供了框架。
论文解读文章
一、研究背景与科学问题
枸杞属（Lycium）隶属于茄科（Solanaceae），全球约有70至80个物种，广泛分布于干旱和半干旱地区。其中，宁夏枸杞（L. barbarum，红果枸杞）和黑果枸杞（L. ruthenicum，黑果枸杞）因富含类胡萝卜素、多糖和花青素等生物活性化合物而具有重要的经济和药用价值。尽管这两个物种在分类学上亲缘关系密切，但它们在果实颜色和生态适应性方面表现出显著的表型差异。因此，阐明这种分化背后的演化模式，需要对其遗传架构进行深入研究，其中细胞器基因组尤为重要——它们为母系遗传和谱系分化提供了宝贵的见解。
植物质体基因组通常具有保守的四分体结构，即两个反向重复序列（inverted repeats, IRs）将大单拷贝区（large single-copy region, LSC）和小单拷贝区（small single-copy region, SSC）分隔开。在此保守架构内，密码子使用模式和序列变异可以提供遗传多样性信息。大规模质体基因组比较数据集还能揭示与网状演化（reticulate evolution）一致的母系不一致模式。然而，仅凭质体数据无法区分杂交或基因渐渗与不完全谱系分选等过程。
与相对稳定的质体基因组相比，植物线粒体基因组以结构多样性和频繁重排著称。线粒体基因组可形成由重复序列介导的分支状或多染色体构型，重排可能产生新的开放阅读框（open reading frames, ORFs）。然而，ORF预测本身并不能确立生物学功能。短读长测序深度比较可以鉴定候选区间尺度的差异，而其精确的分子构型则需要靶向或长读长验证。在枸杞属中，质体基因组的研究远比群体水平线粒体结构变异的研究更为深入。
二、研究内容与主要结论
本研究分析了201份枸杞属材料的细胞器基因组数据。研究人员利用大规模质体基因组比较数据集来表征基因组架构和母系遗传结构，并鉴定与网状母系历史兼容的模式。同时，研究人员以de novo方式组装了一个黑果枸杞线粒体基因组，并将其作为LR（黑果枸杞）和LB（宁夏枸杞）复合群群体水平比较的参考序列。线粒体组分聚焦于从测序深度图谱中鉴定出的区间尺度结构变异。
研究的主要结论包括：质体基因组分析解析出七个母系遗传聚类，鉴定出与红果枸杞复合群复杂母系祖先相符的模式；全部201个质体基因组均保留了保守的四分体结构和79个独特的蛋白质编码基因；Illumina深度图谱鉴定出一个约5.5 kb的候选结构差异区间，该区间在LB代表材料中覆盖度显著降低，对应的参考区间包含orf27、orf92和orf93三个候选ORF。该深度支持的区间界定了一个候选线粒体区域，但其精确分子构型或生物学功能尚待确定。
三、主要关键技术方法
研究人员从201份枸杞属个体中提取基因组DNA，采用Illumina HiSeq X Ten平台进行双端测序，每份材料约产生10 Gb原始数据。质体基因组使用GetOrganelle v1.7.7进行de novo组装，使用Plastome Genome Annotator v2.0进行注释。变异检测使用Snippy v4.6.0。核苷酸多样性（π）使用DnaSP v6计算，非同义（dN）和同义（dS）替换速率使用PAML v4.9中的codeml估算。系统发育分析使用IQ-TREE基于79个共有质体CDS的串联比对进行最大似然（ML）推断。群体结构分析使用ADMIXTURE 1.3，单倍型网络使用PopART v1.7构建。重复序列分析使用MISA-web和REPuter。线粒体基因组使用PMAT v1.5基于单个黑果枸杞材料（SRR28640618）的PacBio高保真长读长进行de novo组装，使用PMGA进行注释，候选ORF使用AlphaFold 3建模和TMHMM进行跨膜螺旋预测。群体水平线粒体结构变异分析使用BWA-MEM将186份材料（LR组n=53，LB复合群n=133）的全基因组测序数据比对至黑果枸杞参考线粒体基因组，使用PanDepth v2.26生成测序深度图谱。
四、研究结果
3.1 大规模枸杞质体基因组数据集的一般特征与核苷酸变异
201份新组装的枸杞质体基因组长度为154,856至155,811 bp，GC含量稳定在37.8%–37.9%。每个组装均呈现典型的四分体结构。每个质体基因组包含113个独特基因，包括79个CDS、30个tRNA和4个rRNA；计入反向重复序列的重复拷贝后，每个质体基因组包含131个注释基因拷贝。共检测到49,250个变异，包括33,296个SNV、7,494个缺失、6,027个插入和2,433个复杂变异，SNV约占67.6%。LSC和SSC区域的变异密度高于IR区域。在77个保留的CDS中，ycf1变异最丰富，其次为rpoC2和ndhD。滑动窗口分析显示核苷酸多样性（π）范围为0至0.00728，均值为0.0012，IR区域多样性最低，LSC区域最高，突变热点包括atpI、rps4-trnT和cemA。
3.2 枸杞属内的母系关系与遗传结构
基于79个共有质体CDS串联比对构建的ML树解析出七个母系遗传聚类：CP（n=3）、LC（n=4）、LR（n=53）、LBC（n=53）、LCYN（n=8）、LB1（n=73）和LB2（n=7）。LB1、LB2和LBC彼此亲缘关系更近，其中LB1和LB2主要包含宁夏枸杞材料，LBC以宁夏枸杞var. chinense为主。LC和LCYN聚类在ML树和PCA中同样关系密切，共同代表与L. chinense相关的材料。LR聚类主要由黑果枸杞组成，与红果枸杞谱系距离较远。CP聚类包含三个北美物种，形成最分歧的类群。群体结构分析显示K=6时交叉验证误差最低。成对分化以CP和LBC之间最高（Fst=0.99）。单倍型网络解析出44个单倍型，支持LB1-LBC-LB2复合群的紧密分组。
3.3 重复序列与SSR多态性表征
MISA-web分析鉴定出三类SSR，包括单核苷酸、二核苷酸和三核苷酸重复。单核苷酸重复最为丰富，占SSR总数的97.00%，其中A/T和AT/AT基序占所有SSR的98.89%，表明枸杞属及更广泛的茄科质体基因组中SSR具有强烈的A/T富集偏好。散在重复序列分析共鉴定出121个重复元件，其中回文重复（55个，45.45%）最丰富，其次为正向重复（49个，40.50%）、反向重复（11个，9.09%）和互补序列（6个，4.96%）。散在重复序列长度从30至100 bp不等，以40–49 bp类别最为常见。
3.4 质体基因密码子使用模式与选择压力
RSCU分析显示枸杞属蛋白质编码基因偏向A/T结尾密码子。在64个密码子中，30个同义密码子（大多数以A或U结尾）的RSCU值大于1。亮氨酸对UUA的偏好最强（1.92），其次为丙氨酸对GCU（1.77）和精氨酸对AGA（1.73）。ENC-GC3s图中所有基因均位于标准曲线下方，表明突变偏差单独无法解释观察到的密码子使用模式。ENC值范围为22.97（petN）至51.74（ycf2）。基于77个保留CDS的替换速率估计显示，大多数基因的dS超过dN且dN/dS低于1，与纯化选择占主导一致。ycf1、rpl32、rpl36、rps19和cemA表现出相对较高的替换速率。
3.5 黑果枸杞与宁夏枸杞线粒体基因组比较分析
单个黑果枸杞参考线粒体基因组以具有替代路径的分支状组装图表示，与保守的质体四分体架构形成对比。该402.4 kb参考序列具有稳定的GC含量谱和广泛的非编码基因间区。LR和LB组表现出显著的表型差异。通过PanDepth v2.26生成的Illumina测序深度图谱，鉴定出两组之间约5.5 kb的候选线粒体结构差异，表现为LB代表材料中连续的覆盖度降低。对应的参考区段包含orf27、orf92和orf93。AlphaFold 3模型和TMHMM预测表明候选蛋白缺乏预测的跨膜螺旋。
五、讨论与结论
讨论部分指出，201个质体基因组的基因组大小、GC含量、基因顺序和四分体架构高度保守。ycf1的高变异性与其在被子植物中的快速进化一致。母系模式显示LB1-LBC-LB2复合群内部关系密切，LR分化更深，但质体拓扑结构本身无法确定分化时间或建立唯一的演化轨迹。CP聚类仅包含三个北美物种各一个材料，其分歧位置最好解释为地理和分类上独特的母系谱系，而非欧亚聚类。重复序列分析显示强结构保守性，鉴定的SSR和散在重复序列可为区分母系遗传组提供标记。密码子使用模式可能反映突变偏差和选择的共同贡献，大多数基因的低dN/dS比值与广泛纯化选择一致。线粒体参考基因组的分支状组装图与植物线粒体中重复相关重组兼容，但当前数据无法量化重组频率或确定每条路径是否代表稳定的in vivo分子。深度图谱鉴定出的5.5 kb候选区间支持区间尺度结构差异，但不确立其精确分子构型或功能效应。候选ORF的结构预测不确立其表达、翻译或生理活性。
研究结论部分指出，本研究为201份枸杞属材料提供了比较细胞器基因组资源。质体基因组数据解析出七个母系遗传聚类，并鉴定了保守的基因组架构以及具有谱系信息的序列变异。基于单个黑果枸杞参考线粒体基因组的深度比较鉴定出LR和LB之间约5.5 kb的候选结构差异。该深度支持的区间是枸杞属线粒体差异的候选区域。未来研究应调查候选ORF的生物学功能并检验其在枸杞属中的演化历史。
