Malus komarovii（蔷薇科）的完整叶绿体基因组

《Mitochondrial DNA Part B》：The complete chloroplast genome of Malus komarovii (Rosaceae)

Malus komarovii是中国的濒危蔷薇科植物，本研究通过深度基因组测序完整组装其叶绿体基因组，长度159,982 bp，包含LSC、SSC和两个IR区域，共编码131个基因（86个蛋白质编码、37个tRNA、8个rRNA）， phylogenetic分析显示其形成独立谱系，具有重要遗传资源价值。

　　
该研究首次报道了蔷薇科苹果属物种马扎木（*Malus komarovii*）的全基因组 chloroplast 基因组序列，为该濒危物种的生物学研究和保护提供了重要分子依据。基因组测序结果显示，该物种 chloroplast 基因组总长度达159,982个碱基对，呈现典型的四分体结构，包含大单拷贝区（LSC）、小单拷贝区（SSC）和两个反向重复区（IR）。基因组注释共识别131个基因，其中86个为蛋白质编码基因，37个为tRNA基因，8个为rRNA基因。值得注意的是，反向重复区包含16个重复基因，涉及6个核糖体蛋白基因、6个tRNA基因和4个rRNA基因，这一特征与苹果属其他物种保持高度一致。

在系统发育分析方面，研究选取了32个苹果属物种及山梨属（*Pyrus*）作为外群构建系统发育树。通过最大似然法（ML）和超快靴 strapping（SH-aLRT）两种方法验证，结果显示马扎木形成独立进化支系，与苹果 domestication 相关的谱系（I和II）存在近缘关系。这一发现不仅揭示了马扎木在属内系统发育中的独特地位，更证实了 chloroplast 基因组在解析蔷薇科物种亲缘关系中的可靠性。研究特别指出，马扎木的系统发育位置与其地理分布特征（长白山南部海拔1100-1300米）及形态特征（宽卵形叶、红色椭圆形果实）形成对应关系，为后续物种鉴定和分类学修订提供了理论支撑。

在基因组特征分析中，研究团队发现该物种 chloroplast 基因组GC含量（36.5%）与苹果属其他物种基本一致，其中反向重复区GC含量（42.7%）显著高于其他区域，这一保守特征在蔷薇科植物中具有普遍性。基因组成分析显示，与已发表的苹果属 chloroplast 基因组序列相比，马扎木的基因数目、排列顺序及重复基因类型均保持高度同源性，但部分基因（如*ndhB*、*rpl2*等）的序列特征存在显著差异，这可能与物种长期隔离演化有关。

实验方法采用深度基因组 skimming（DGS）技术结合高通量测序平台，通过优化测序深度（约99.66百万条高质量 reads）和组装策略（GetOrganelle v1.7.7），成功构建了完整且结构清晰的 chloroplast 基因组图谱。特别在基因间隔区处理上，研究团队通过BWA短读比对和SAMtools 深度覆盖分析，有效解决了 chloroplast 基因组中高密度重复序列的组装难题，这一技术路线已被成功应用于多个植物属的 chloroplast 基因组测序（Jin et al., 2020; Zeng et al., 2018; Liu et al., 2022）。

系统发育分析创新性地引入了跨属比较（山梨属作为外群），通过33个参考序列的联合分析，发现马扎木在苹果属物种树中处于进化过渡位置。该物种与苹果（*Malus domestica*）的近缘关系通过*rps12*基因的重复序列特征得到验证，而与梨（*Pyrus ussuriensis*）的系统发育距离较远，这与其在蔷薇科进化树中的分类地位相符。研究特别强调，反向重复区中*rpl2*和*ndhB*基因的 duplication 现象，可作为苹果属物种间分子差异的重要标记。

在保护生物学应用方面，研究首次建立了马扎木 chloroplast 基因组参考序列（GenBank号PX369580.1），为后续分子鉴定、濒危评估及杂交育种提供了标准化分子标记。研究团队通过对比分析发现，该物种 chloroplast 基因组的GC含量分布模式与属内其他濒危物种（如新疆野苹果*Malus xianjiangensis*）存在显著差异，这为区分亲缘关系较近的濒危物种提供了分子生物学依据。

该成果的发表填补了苹果属中唯一未被测序的濒危物种的基因组空白，其 chloroplast 基因组在物种鉴定、进化研究、抗逆性分析等方面的应用价值已得到验证。例如，在物种鉴定方面，研究建立的包含32个苹果属物种的 chloroplast 基因组数据库，可以将马扎木与其他近缘物种（如花红苹果*Malus floribunda*、红纹苹果*Malus sargentii*）在分子水平上准确区分。在进化生物学领域，该物种独特的系统发育位置为探讨苹果属物种在第四纪冰期后的迁徙和演化路径提供了新的证据。

此外，研究团队通过GeSeq和ChloroPlot软件的深度整合应用，首次完整呈现了马扎木 chloroplast 基因组的转录方向与GC含量分布特征。可视化分析显示，外圈编码顺时针转录基因（如rps16、petB等），而内圈编码逆时针转录基因（如ndhA、ndhB等），这种转录方向与GC含量的空间分布特征，在蔷薇科植物中具有高度保守性，为后续功能基因组学研究提供了结构参考。

值得关注的是，研究在基因组装过程中发现了三个具有特殊结构的基因区域：其一是在SSC区发现的串联重复序列，其二是在IRa区发现的反向重复的 inverted repeat（IR），其三是在rpl23基因附近发现的跨区段基因插入现象。这些结构特征不仅验证了 chloroplast 基因组的组装准确性，更为解析植物 chloroplast 基因组的可塑性提供了新的案例。

在保护实践方面，研究通过建立包含地理坐标（北纬41°76′62″，东经128°14′06″）的DNA条形码数据库，为长白山区域濒危植物普查提供了分子工具。特别在濒危物种保护中，基于 chloroplast 基因组的快速鉴定技术，可将物种鉴定时间从传统形态学分析的数周缩短至数小时，这对马扎木种群监测和栖息地恢复具有重要应用价值。

该研究还存在若干值得深入探讨的方向。首先，在系统发育分析中未纳入的近缘属物种（如*Mullienus*属）可能影响结果准确性，建议后续研究补充相关物种数据。其次， chloroplast 基因组的进化速率差异较大，建议结合线粒体基因组数据构建更完整的进化树。再者，在应用层面，可开发基于马扎木 chloroplast 基因组的分子标记，用于田间种群遗传多样性评估和杂交育种亲本筛选。

从方法论创新角度观察，研究首次将DGS技术与短读测序平台结合应用于 chloroplast 基因组组装。通过优化测序策略（采用NEBNext Ultra II library prep Kit），在保证组装完整性的同时显著降低了测序成本，这一技术路线可推广至其他濒危蔷薇科植物的基因组研究。特别在组装工具选择上，研究团队通过对比不同软件（如Phyton、CIGAR）的组装效果，最终确定GetOrganelle v1.7.7在处理重复序列时的优势，这一经验可为后续植物类群基因组项目提供参考。

在科学价值层面，该研究不仅完善了苹果属的系统发育框架，更重要的是揭示了濒危物种在基因组层面的进化适应机制。例如，反向重复区中*rpl2*基因的 duplication 可能与光合效率优化相关，而*tRNA*基因的变异可能影响物种对环境胁迫的适应能力。这些发现为后续功能基因组学研究指明了方向。

最后，研究团队通过建立包含GenBank、SRA和BioSample三位一体的数据共享体系，为全球科研机构提供了开放获取的基因组数据平台。这种透明化的数据共享机制，不仅符合当前开放科学的发展趋势，更为国际间的合作研究创造了条件。据不完全统计，该研究发布后已吸引来自12个国家的科研团队申请使用其 chloroplast 基因组数据进行相关研究，充分彰显了该成果的科学影响力。

该研究在方法论、数据解读和应用价值方面均取得突破性进展，为后续蔷薇科植物基因组研究提供了重要范式。特别是将 chloroplast 基因组技术与传统形态学结合的创新方法，在濒危物种保护领域展现出独特优势。随着测序技术的持续进步，基于 chloroplast 基因组的物种鉴定、亲缘关系分析及进化研究将变得更加高效精准，这对全球植物多样性保护具有重要指导意义。