普陀杜鹃（Rhododendron simsii var. putuoense）完整线粒体基因组特征解析及其比较进化分析

《Frontiers in Plant Science》：Characterization of the complete mitochondrial genome of Rhododendron simsii var. putuoense and its comparative evolutionary analysis

　　

摘要专业翻译 普陀杜鹃（Rhododendron simsii var. putuoense）是一种具有重要生态和遗传资源价值的沿海特有种。尽管其叶绿体基因组（cpDNA）近期已被解析，但其线粒体基因组（mtDNA）仍属空白。本研究利用PacBio HiFi长读长测序技术，首次完成了普陀杜鹃mtDNA的全长组装与注释。该基因组序列全长为658,161 bp，由两个亚基因组线性分子（一个主要分子633,105 bp和一个次要分子25,056 bp）构成，共编码33个蛋白质编码基因（PCGs）、26个tRNA、3个rRNA及323个开放阅读框（ORFs）。重复序列分析显示散在重复占主导地位（72%），另有128个简单序列重复（SSRs）和20个串联重复。该mtDNA表现出强烈的A/U结尾密码子使用偏好，预测到463个RNA编辑位点。6个同源片段（占mtDNA的0.91%）为细胞器间基因转移（IGT）提供了证据。值得注意的是，成对Ka/Ks分析显示matR基因的比值升高，提示其存在谱系特异性进化动态，有待进一步开展位点特异性密码子模型检验。系统发育和共线性分析进一步确认了其与杜鹃（R. simsii）的密切亲缘关系，同时揭示了与生态特化相一致的结构分化。总体而言，这些发现加深了对杜鹃属（Rhododendron）线粒体基因组多样性的理解，并为这一具有重要生态意义的变种谱系特异性资源的保护建立了遗传学框架。

论文解读

一、研究背景与意义

线粒体基因组（mtDNA）作为植物能量代谢的核心细胞器基因组，其研究对于理解植物遗传、进化和适应具有重要意义。与核基因组的高密度编码和叶绿体基因组（cpDNA）的相对保守性不同，植物线粒体基因组最显著的特征在于其高度动态的结构和多样化的内容来源。随着高通量测序技术的发展，分析非模式植物的mtDNA已成为探索植物适应性进化的重要切入点。

杜鹃属（Rhododendron）作为杜鹃花科（Ericaceae）中进化最活跃的类群之一，其复杂的分类系统长期是植物系统学研究的难点。目前该属的细胞器基因组研究主要集中在叶绿体基因组层面，为系统发育重建提供了重要证据。然而，该属mtDNA的结构特征、进化压力及其与生态适应的关系尚未得到系统解决。这一空白在一定程度上限制了对杜鹃属进化机制的全面理解。

普陀杜鹃（Rhododendron simsii var. putuoense）是杜鹃（R. simsii）的一个变种，特生于中国东部沿海岛屿，地理分布范围狭窄。其独特的沿海栖息地以高盐度和强风等环境胁迫为特征，在基因组层面塑造了独特的进化性状。尽管其叶绿体基因组近期已被解析，但其mtDNA仍完全未被探索，对该分类群的细胞器基因组综合分析仍然缺乏。因此，开展普陀杜鹃mtDNA的组装与注释研究，对于理解杜鹃属mtDNA的结构可塑性和进化动态，以及为这一生态重要变种的保护与利用奠定遗传学基础，具有重要的科学意义。

二、研究内容与主要结论

研究人员利用PacBio HiFi长读长测序技术，首次完成了普陀杜鹃mtDNA的全长组装与注释，并通过比较基因组学、共线性分析、系统发育重建和选择压力分析等手段，系统揭示了其基因组结构、进化压力及与近缘种的系统发育关系。

普陀杜鹃mtDNA全长为658,161 bp，GC含量为45.84%，由两个线性亚基因组分子构成，其中主要分子为633,105 bp，次要分子为25,056 bp。基因组共注释到64个基因，包括33个蛋白质编码基因（PCGs）、26个tRNA基因、3个rRNA基因和2个假基因，另鉴定到323个开放阅读框（ORFs）。33个PCGs按功能可分为11类，其中7类属于呼吸链核心基因，包括ATP合酶（5个基因）、细胞色素c生物发生（4个基因）、泛醇-细胞色素c还原酶（1个基因）、细胞色素c氧化酶（3个基因）、成熟酶（1个基因）、转运膜蛋白（1个基因）和NADH脱氢酶（9个基因）；其余4类为非核心功能基因，包括大核糖体蛋白（LSU）（3个基因）、小核糖体蛋白（SSU）（6个基因）、琥珀酸脱氢酶（2个基因）和核糖体RNA（3个基因）。

RNA编辑分析共预测到463个编辑位点分布于33个mt PCGs中，其中mtB基因拥有最多的RNA编辑位点（39个），其次为ccmB和nad2，各含38个位点。编辑事件中49.67%导致氨基酸从亲水性转变为疏水性，0.22%导致提前终止密码子。

密码子使用偏好分析共检测到9,987个密码子。亮氨酸（Leu）是最丰富的氨基酸，占11.43%。31个密码子的相对同义密码子使用度（RSCU）值大于1，其中90.32%以A或U结尾，表明该mtDNA中高频使用密码子明显偏好A/U结尾。

重复序列分析共鉴定到527个重复序列。散在重复最为丰富，占72%，包括170对正向重复（F）和209对回文重复（P）；简单序列重复（SSRs）128个，占24%；串联重复20个，仅占4%。

叶绿体与线粒体基因组间同源片段分析共鉴定到6个同源片段，总长5,967 bp，占mtDNA的0.91%。其中5个片段完全位于同源序列内且均为tRNA基因（trnA-TGC、trnW-CCA、trnD-GTC、trnN-GTT和trnM-CAT），1个片段部分位于同源序列内，属于核糖体RNA基因trn18。

K_a/K_s分析表明，普陀杜鹃与杜鹃（R. simsii）比较时总体K_a/K_s值为0.65，与R. × pulchrum比较时为0.61，与R. decorum比较时为0.47，与R. vialii比较时为0.32，与R. nivale subsp. boreale比较时为0.31。几乎所有PCGs的K_a/K_s比值均小于1，表明大多数PCGs在进化过程中经历了纯化选择。值得注意的是，matR基因的K_a/K_s比值大于1，提示其可能存在升高的进化速率或正选择信号。

系统发育分析显示，普陀杜鹃与杜鹃（R. simsii）和R. × pulchrum聚为一个独立分支，自展支持率为100%，表明三者间具有极为密切的系统发育关系。共线性分析进一步支持了这一结论，普陀杜鹃与杜鹃（R. simsii）的共线性比例最高（82.97%和81.78%），其次为R. × pulchrum（75.85%和71.82%），而与R. nivale subsp. boreale的共线性比例最低（47.64%和62.45%），反映了广泛的基因组重排和更大的进化分歧。

三、讨论与结论

讨论部分围绕四个主题展开。第一，关于转录后调控特征与mt基因组进化稳定性，463个RNA编辑位点的数量与其他杜鹃属物种相当，表明该变种具有完整的RNA编辑系统；编辑后亲水性到疏水性转变占主导，有助于维持线粒体呼吸链跨膜蛋白的疏水性；A/U结尾密码子偏好与大多数被子植物mtDNA的典型特征一致；散在重复（72%）的主导地位为多分支构象结构和种间共线性变异提供了机制性解释。第二，关于进化压力与栖息地适应，matR基因K_a/K_s比值大于1可能暗示谱系特异性分化，但研究人员谨慎指出该结果需要位点特异性密码子模型、群体水平取样和功能实验加以验证。第三，关于细胞器间基因转移（IGT）及其进化意义，0.91%的转移比例表明基因组边界相对稳定，5个完整tRNA基因的富集反映了mt基因组的功能补偿机制。第四，关于系统发育关系与分类地位，线粒体基因组证据为普陀杜鹃作为杜鹃变种的分类处理提供了分子层面的支持，但mtDNA单亲遗传的特性决定了该证据主要反映母系谱系，未来需结合双亲核标记进行更全面的分析。

研究结论指出，本研究首次完成了普陀杜鹃mt基因组的完整组装与注释，揭示了由两个线性亚基因组分子构成的多分支构象结构，总长658,161 bp。该基因组在基因组成、RNA编辑、密码子使用偏好和重复序列分布方面表现出被子植物mt基因组的典型特征。更重要的是，研究揭示了包括IGT和matR基因序列分化在内的进化动态，为该沿海特有种的线粒体基因组变异性提供了初步见解。系统发育和共线性分析进一步确认了其与杜鹃（R. simsii）的密切关系，而独特的选择压力揭示了生态特化的轨迹。这些发现不仅阐明了该mt基因组的结构和进化复杂性，还突显了谱系特异性基因组多样性的保护价值，为保护该变种的适应潜力以及指导杜鹃属内细胞质-核共进化研究提供了关键的遗传学基础。