《Frontiers in Plant Science》：Selective regime shifts and intracellular gene transfer in the organellar genomes of the mycoheterotrophic orchid Galeola nudifolia Lour. (Orchidaceae)
编辑推荐：
由光合自养向菌异养（mycoheterotrophy）营养方式的转变在兰科（Orchidaceae）中多次发生，并预期会对质体基因组（plastome）与线粒体基因组（mitogenome）产生相反的进化信号。然而，这两个细胞器区室之间的定量比较仍然稀少，尤其
由光合自养向菌异养（mycoheterotrophy）营养方式的转变在兰科（Orchidaceae）中多次发生，并预期会对质体基因组（plastome）与线粒体基因组（mitogenome）产生相反的进化信号。然而，这两个细胞器区室之间的定量比较仍然稀少，尤其在香荚兰亚科（Vanilloideae）内更是如此。研究人员报道了香荚兰族（Vanilleae）内完全菌异养兰类无叶珊瑚兰（Galeola nudifolia Lour.）的首套成对质体基因组与线粒体基因组，并将这些数据与已发表的菌异养及光合兰科细胞器基因组整合分析。该线粒体基因组为852,490 bp，由24条线性重叠群（contig）组成，编码38个蛋白编码基因、8个转运RNA（tRNA）基因和3个核糖体RNA（rRNA）基因；rpl2的丢失经靶向重新组装确认。101,856 bp的质体基因组保留典型四部分结构（quadripartite structure），但光合作用相关基因广泛丢失。研究人员在 mitochondrial genome 中检测到42个质体来源区域，对应20个质体蛋白编码基因；这些区域均无完整开放阅读框（open reading frame，ORF），表明物理转移广泛但功能转移缺失。分支模型（codeml）与RELAX分析显示两基因组选择机制相反：质体基因组中纯化选择（purifying selection）的放松集中于光合作用与ATP合成相关位点（atpI经多重检验校正后显著，padj = 0.033，更多ATP合酶基因获RELAX支持），而线粒体基因组大体维持纯化选择。预测的RNA编辑密度在营养型组间无实质差异，尽管该分析为探索性。在质体系统发育中，Cyrtosia相对于G. nudifolia呈并系，但取样有限且支持弱，需更密取样与核数据验证。
研究背景方面，陆生植物中向非光合生活方式的转变反复独立发生，兰科中菌异养起源超过30次。质体基因组编码光合作用相关基因，丧失光合能力后这些基因受放松的纯化选择并逐步丢失；线粒体基因组编码呼吸与ATP合成基因，不论碳获取方式如何都须维持功能。然而，菌异养类群中两个细胞器基因组不对称进化的定量检验稀缺，且已有细胞器成对数据在香荚兰亚科等类群中极不充分。香荚兰亚科同时含光合与完全菌异养属，但此前仅有质体基因组和一种未充分注释的Cyrtosia lindleyana线粒体草图，无法严谨检验菌异养下选择机制与系统发育一致性。无叶珊瑚兰（Galeola nudifolia）是香荚兰族广布完全菌异养种，研究人员因此测序组装其成对细胞器基因组，并与已发表兰科数据整合，相关论文发表于《Frontiers in Plant Science》。
主要关键技术方法方面，样本采自马来西亚霹雳州玲珑村一带无叶珊瑚兰全株，为避免内生真菌污染优先用硅干花组织与苞片提基因组DNA，经Illumina NovaSeq 6000双端测序。细胞器组装采用GetOrganelle与NOVOPlasty种子组装、BLASTN筛选、Geneious合并去冗余、SPAdes混合组装，得到24条线性线粒体contig与环状质体序列；rpl2丢失用靶向NOVOPlasty重装验证。系统发育用MAFFT比对、trimAl裁剪、ModelFinder选模、IQ-TREE2最大似然及MrBayes贝叶斯推断。选择压力用PAML codeml分支模型比较光合与菌异养分支非同义/同义替代比（ω = dN/dS），并用HyPhy RELAX检验选择强度参数K。质体到线粒体转移用BLASTN同源比对、BEDTools合并区域、ORF完整性判断。RNA编辑用Deepred-Mt预测，简单序列重复（SSR）与散在重复分别用Phobos和Geneious重复查找。
3.1 细胞器基因组特征
研究人员通过混合组装与注释发现，无叶珊瑚兰线粒体基因组总长达852,490 bp、GC含量46.7%、平均覆盖651.98倍，含38个蛋白编码基因、8个tRNA基因、3个rRNA基因；rpl2在29个其他兰类线粒体中也于11种缺失，靶向重装证明确为真实丢失。质体基因组101,856 bp、GC 33.7%、覆盖7,627.66倍，保留大单拷贝（LSC）、小单拷贝（SSC）与两段倒置重复（IR）四部分结构，注释64个唯一基因，光合作用相关ndh、psa、psb、pet、rbcL等大量丢失或假基因化。线粒体预测RNA编辑位541个且全在nad基因，菌异养与光合组编辑密度无实质差异。散在重复36处中34处仅在线粒体、2处仅在质体，无跨细胞器共享。质体来源区域42处对应20个质体蛋白编码基因，但均截断或含内部终止密码子，无完整读框，说明物理转移广泛而功能转移缺失。
3.2 系统发育关系
研究人员用38个线粒体蛋白编码基因加3个rRNA构建最大似然与贝叶斯树，二者拓扑一致：无叶珊瑚兰与Cyrtosia lindleyana成支，该香荚兰亚科枝姊妹于Chamaegastrodia shikokiana（树兰亚科），菌异养属具长枝但主干支持率高。质体树取样更广，Cyrtosia nana与无叶珊瑚兰成高支持枝，C. lindleyana以低支持位于其姊枝，使Cyrtosia相对Galeola呈并系；研究人员指出该安排因仅3个Cyrtosia和1个Galeola且无核数据而属初步。
3.3 细胞器内替代率变异与选择压力
分支模型下，线粒体基因前台ω略升但类别间Δω近零且无多重检验显著者；质体基因明显放松纯化选择，atpI前台ω 0.41、后台0.10且padj = 0.033，atpB、atpE、atpH、psbE等ATP合酶与光系统基因Δω为正。RELAX中11个基因显著：质体atpA、atpB、atpE、atpH、atpI、psbE、matK放松（K < 1），ycf2与rps7强化（K > 1）；线粒体仅nad6放松、atp4强化。结论为菌异养下质体光合与ATP合成基因纯化选择放松，线粒体呼吸基因大体维持纯化选择。
讨论部分总结中，研究人员指出兰科线粒体结构高度异质，contig数与大小不随营养型变化，无叶珊瑚兰与C. lindleyana共有rpl2丢失且该丢失跨光合与菌异养类群出现。质体大小在Vanilleae菌异养梯度中由Galeola、Cyrtosia到Lecanorchis递减，支持营养型与质体退化相关；但IR收缩/扩张在光合Vanilla、Pogonia及树兰亚科Chamaegastrodia中同样发生，故结构重排是谱系特异而非营养型驱动。质体到线粒体转移（MTPT）中，无叶珊瑚兰42区对应20个质体基因，部分光合作用基因在其质体已失，却仍存于线粒体为残段；trnN-GUU与trnR-ACG邻接质体同源区，且两tRNA分化程度不同，支持tRNA介导的多次独立细胞内基因转移（intracellular gene transfer，IGT）。所有质体同源序列均无完整ORF，物理转移广泛但无功能 relocation。
研究结论部分译为：香荚兰亚科完全菌异养兰无叶珊瑚兰的首套成对质体与前所报道不同的线粒体数据显示，质体基因组因光合能力丧失而发生基因丢失与纯化选择放松，尤其atpI及ATP合酶相关位点；线粒体基因组在菌异养下仍保留呼吸功能相关基因的纯化选择，仅少数基因如nad6放松、atp4强化。质体来源序列大量转入线粒体，但均无完整开放阅读框，表明细胞内基因转移以物理残留为主而无功能性细胞器间基因置换。质体系统发育受长枝吸引偏向，线粒体系统发育提供更独立类群关系检验；Cyrtosia对Galeola并系为取样限制下的初步结论。该研究填补了Vanilloideae成对细胞器基因组空白，证明菌异养转型中两细胞器基因组呈结构、功能与选择机制上的解耦进化，为兰科非光合转型下的细胞器基因组演化提供了量化证据。
