《Ecology and Evolution》：Exploring the Mitochondrial Genomes of Phoebe Species (Lauraceae): Structural Dynamics and Functional Conservation
编辑推荐：
植物线粒体基因组（mitogenome）在基因组大小和结构架构上存在显著差异，尽管其序列进化速率总体较慢。楠属（Phoebe）是樟科（Lauraceae）中一个具有重要生态和经济价值的属，但其线粒体基因组的多样性至今仍缺乏系统表征。本研究对来自中国的三种国家二
植物线粒体基因组（mitogenome）在基因组大小和结构架构上存在显著差异，尽管其序列进化速率总体较慢。楠属（Phoebe）是樟科（Lauraceae）中一个具有重要生态和经济价值的属，但其线粒体基因组的多样性至今仍缺乏系统表征。本研究对来自中国的三种国家二级保护野生植物——刨花楠（P. bournei）、浙江楠（P. chekiangensis）和桢楠（P. zhennan）的线粒体基因组进行了重新测序、组装和注释，并比较了它们的线粒体基因组特征。三个组装结果被解析为代表性的环状结构，大小范围为808至864 kb，具有相似的GC含量和保守的蛋白质编码能力。每个线粒体基因组包含41个不同的蛋白质编码基因（protein-coding genes, PCGs）、27–28个转运RNA（transfer RNAs, tRNAs）和3个核糖体RNA（ribosomal RNAs, rRNAs）。共线性分析（synteny analysis）显示，尽管三个物种之间存在显著的序列同源性，但其同源片段的顺序和方向发生了广泛变化。大量重复序列主要分布于非编码区，而叶绿体来源片段（plastid-derived fragments）则记录了历史上的细胞内基因转移事件。三个物种在密码子使用模式和预测的RNA编辑模式上表现相似，而较低的同义替换率限制了对配对比值的推断。基于线粒体蛋白质编码基因的系统发育分析支持楠属为一个高支持度的单系谱系。这些结果揭示了结构上的显著分化伴随核苷酸组成和编码能力的保守性，为进一步理解楠属及樟科植物线粒体基因组的进化变异提供了宝贵数据。
**论文解读：樟科楠属线粒体基因组的比较基因组学研究**
**1. 研究背景与科学问题**
线粒体作为真核细胞中负责细胞呼吸和能量代谢的关键细胞器，其自身基因组在进化过程中经历了广泛的缩减和与核基因组的整合。与真菌和动物线粒体基因组相比，植物线粒体基因组在基因组大小和结构架构上表现出极大的差异，其大小变异幅度超过200倍，从槲寄生属物种的66 kb至银杉的18.99 Mb不等。尽管植物线粒体编码序列相对保守，但其基因组常以动态的环状、线性、分支、多部分以及重组衍生等多种形式存在。这种结构上的高度动态性主要归因于非编码区的扩张，包括重复元件、细胞内基因转移、外源插入及未知起源序列的积累。与此形成鲜明对比的是，植物线粒体编码序列进化缓慢，其同义替换速率低于质体基因组和核基因组。因此，探究这种结构多变性与序列保守性并存的模式，对理解植物线粒体基因组的稳定性、结构进化及其系统发育应用具有重要意义。然而，由于结构复杂性和组装挑战，被子植物线粒体基因组的研究远不如质体基因组深入。
楠属隶属于樟科，包含约100种分布于热带和亚热带亚洲的物种，在中国约有30种，主要分布于长江流域及其南部地区。楠属物种多为高大乔木，其木材耐久、抗腐、芳香且纹理美观，长期以来被用于建筑、家具和雕刻，具有重要的经济价值。同时，楠属物种在亚热带常绿阔叶林中扮演着重要的生态角色，具有可观的碳汇能力和涵养水源、防止水土流失的生态功能。在中国，包括刨花楠、浙江楠、桢楠和闽楠在内的四种楠属植物被列为国家二级重点保护野生植物。尽管楠属物种的质体基因组已被广泛研究，且部分物种已获得染色体水平的核基因组，但目前经过严格验证的楠属线粒体基因组资源和相关比较分析仍然十分缺乏。尤其是已有公开的所谓刨花楠线粒体基因组序列，其原始样本并非真正的刨花楠，而是以刨花润楠的名义提交的，这给后续研究带来了分类学上的混乱。因此，开展基于严正确认样本的多种楠属植物线粒体基因组测序与精细比较，是填补这一空白、厘清物种间线粒体基因组演化规律的关键步骤。
**2. 主要技术方法**
研究人员采集了来自贵州省林业科学院（刨花楠和桢楠）以及南京中山植物园（浙江楠）的三种楠属植物新鲜叶片样本，样本物种身份经分类学专家确认。基因组DNA采用改良CTAB法提取。刨花楠和浙江楠采用PacBio Revio平台进行HiFi测序，桢楠则采用Oxford Nanopore Technologies平台的PromethION进行测序。随后运用PMAT工具进行线粒体基因组组装，并结合Bandage进行组装图谱的检查与手动调整。基因注释使用PMGA和GeSeq工具，并以香樟和浙江桂为参考序列。分析流程包括：使用CPStools、Tandem Repeats Finder和REPuter进行重复序列分析；CodonW分析密码子使用偏好；Deepred-Mt预测RNA编辑位点；PAML估算非同义/同义替换率；利用BLAST+鉴定线粒体中叶绿体来源的DNA片段（MTPTs），并通过PyGenomeViz进行共线性分析；最后基于线粒体蛋白质编码基因使用IQ-TREE构建最大似然系统发育树，包括13个樟目物种作为内群和2个木兰目物种作为外群。
**3. 研究结果**
**3.1 线粒体基因组特征**
三个新组装的线粒体基因组均被解析为环状结构，大小分别为刨花楠847,788 bp、浙江楠863,663 bp、桢楠807,955 bp。三者的总GC含量高度保守，介于47.04%至47.11%之间。每个线粒体基因组均包含41个不同的蛋白质编码基因（PCGs），其中包括24个经典核心基因。浙江楠的rps12基因存在额外拷贝，而其余基因均为单拷贝。三个基因组还编码27–28个转运RNA（tRNAs）和3个核糖体RNA（rRNAs）。tRNA含量的差异主要表现为基因存在与否及拷贝数变化。
**3.2 重复序列分析**
在三个线粒体基因组中共鉴定出850个简单重复序列（SSRs），以四核苷酸重复最为丰富，主要分布于基因间隔区。此外还鉴定出232个串联重复和4545个散在重复，其中正向和回文重复占主导。大多数散在重复位于基因间隔区，表明非编码区是重复序列的主要载体。
**3.3 密码子使用模式**
密码子使用分析显示，三个物种的GC3（39.02%–39.48%）均低于GC1和GC2，偏好使用以A或U结尾的密码子。ENC值在53.74至54.13之间，表明整体密码子使用偏倚较弱。RSCU分析显示三个物种的偏好密码子谱高度相似，其中丙氨酸的GCU和组氨酸的CAU均表现出持续较高的RSCU值。
**3.4 预测RNA编辑位点**
分别在刨花楠、浙江楠和桢楠的41个线粒体PCGs中预测到764、773和797个C-to-U RNA编辑位点。nad4基因包含最多的编辑位点（66–68个），其次为nad5（48个）。编辑位点主要集中于密码子第二位（占60.8%–61.7%），其次是第一位，第三位最少。
**3.5 核苷酸替换模式**
配对Ka/Ks比值估算中，由于许多基因对没有同义替换，仅分别有10、6和9个基因在三组两两比较中可计算Ka/Ks比值。部分基因如nad2（1.5522）、mttB（1.3395）和atp6（1.1346）在特定比较中呈现出大于1的比值，但这些仅可作为潜在适应性进化的候选基因，而非正选择的确定性证据。
**3.6 线粒体基因组共线性与结构重排**
成对共线性分析显示，三个楠属物种及其近缘种之间存在广泛的同源区域，但在片段顺序、方向和连续性上差异显著。刨花楠与浙江楠的比较中仅含11个共线性片段，累计长度覆盖了刨花楠线粒体基因组的99.93%，而其他比较中片段数从24到172不等，表明同源区域的排列呈现显著的结构变异。
**3.7 叶绿体来源DNA片段**
三个物种中分别鉴定出67、67和65个线粒体叶绿体来源DNA片段（MTPTs），累计长度占线粒体基因组全长的5.55%至5.76%。每个物种的MTPTs包含5个完整的叶绿体PCGs和15个完整的tRNA基因，且未检测到与叶绿体小单拷贝区对应的片段。
**3.8 系统发育关系**
基于36个共享线粒体PCGs的最大似然系统发育分析高度支持樟科为单系（BS=100）。三个新测序的楠属物种构成一个单系谱系（BS=85），并与润楠属互为姐妹类群。在楠属内部，刨花楠与浙江楠构成姐妹类群。值得注意的是，NCBI中标注为刨花楠的两个序列（OR168698.1和OR168699.1，原始鉴定为刨花润楠）并未与新组装的刨花楠聚在一起，而是高支持度地与润楠属物种聚为姐妹分支，表明这两个序列的物种归属确实存在分类学标注错误。
**4. 讨论与结论**
讨论部分指出，三个楠属线粒体基因组的组合特征体现了被子植物线粒体基因组结构多变性与功能基因含量保守性并存的普遍规律。组装图谱中浙江楠呈现重复介导的8字形拓扑，反映了植物线粒体基因组通过重复序列介导重组形成多种可能结构的复杂性质。tRNA含量的差异可能与细胞内序列转移相关，但仅有序列完整性并不能判断其功能性，转录组证据是必要的。重复序列为同源重组提供了底物，可能与观察到的结构变异有关，但重复丰度本身并不直接证明活跃重组的发生。RNA编辑模式的相似性、密码子使用的偏好性以及部分基因如atp6、nad2和mttB等呈现的较高Ka/Ks比值，共同揭示了线粒体基因在进化过程中的保守特征与潜在的谱系特异性进化压力。系统发育分析中存在的NCBI序列归属错误问题，通过严格的形态学鉴定和系统发育位置分析，研究者认为最简洁的解释是原始样本的物种鉴定错误，并建议重新核查相关凭证标本。总体而言，共线性分析表明快速的结构进化与相对保守的编码序列之间存在部分解耦，重复序列可能是结构变异的基础。
**研究结论**：本研究对刨花楠、浙江楠和桢楠的线粒体基因组进行了系统的比较特征分析。三个线粒体基因组具有相似的核苷酸组成和保守的蛋白质编码能力，但在组装图谱拓扑结构和同源片段排列上存在显著差异，揭示了近缘物种间大量的结构变异。重复序列和MTPTs是形成非编码序列多样性的贡献因素，但其在重组和线粒体功能中的具体作用仍需进一步研究。三个物种间大致相似的基因组成、密码子使用谱和预测的RNA编辑模式表明其线粒体基本特征具有保守性。线粒体系统发育分析支持三个新组装的楠属线粒体基因组之间以及楠属与润楠属之间的密切亲缘关系。这些发现扩大了这些国家保护楠属物种的线粒体基因组资源，为樟科线粒体基因组进化和分类学关系的进一步研究提供了有价值的数据。
