《BMC Genomics》：Characterization of the organelle genomes of mycoheterotrophic Thismia angustimitra and evidence suggesting horizontal gene transfer
编辑推荐：
背景：研究人员对植物线粒体基因组进化的了解远落后于对质体基因组进化的认识。在异养植物中，已知质体基因组退化与光合作用减弱或丧失密切相关，且菌异养植物与吸器寄生植物在质体基因组进化上几乎没有差异。线粒体基因组受异养生活方式影响较小，但也存在例外。与质体基因组进化
背景：研究人员对植物线粒体基因组进化的了解远落后于对质体基因组进化的认识。在异养植物中，已知质体基因组退化与光合作用减弱或丧失密切相关，且菌异养植物与吸器寄生植物在质体基因组进化上几乎没有差异。线粒体基因组受异养生活方式影响较小，但也存在例外。与质体基因组进化不同，线粒体进化可能在菌异养植物与吸器寄生植物之间存在差异，特别是在水平基因转移（horizontal gene transfer, HGT）的发生方面。然而，来自菌异养植物的数据十分匮乏，在十个分支中仅有三个分支拥有完整的线粒体基因组。
结果：研究人员在此展示了来自菌异养植物第四个分支——狭口杯菌（Thismia angustimitra）的完整质体基因组和线粒体基因组数据。其18,215 bp的质体基因组大幅缩减，仅包含少数基因。线粒体基因组呈多环状结构，由22条染色体组成，总长度为422,400 bp。基因含量和RNA编辑模式均正常。线粒体基因组中重复序列极少，重组活性有限。对外源DNA的搜索未发现真菌序列，质体DNA含量较低，但大部分基因间区DNA与远缘被子植物具有相似性。仅有少数质体插入片段似乎源于细胞内转移，而其他片段则极有可能源自水平基因转移。基因间序列的共线性模式表明，水平转移可能双向发生。
结论：T. angustimitra高度退化的质体基因组完全符合异养植物质体基因组退化的普遍规律。尽管其线粒体基因组在许多方面看似不受生活方式影响，但其重组活性异常低下。令人惊讶的是，研究人员发现了暗示水平基因转移的证据，表明这种现象在菌异养植物中可能比之前认为的更为常见。
一、研究背景与科学问题
植物细胞器基因组的研究长期侧重于质体基因组，而对线粒体基因组（mitogenome）的认知相对滞后。在异养植物（包括菌异养植物mycoheterotrophs和吸器寄生植物haustorial parasites）中，质体基因组（plastome）因光合作用功能丧失而发生显著退化，这一现象已被广泛证实。相比之下，线粒体基因组受异养生活方式的影响较小，但其进化模式在不同类群间可能存在差异，尤其是在水平基因转移（horizontal gene transfer, HGT）的发生频率和方向上。然而，目前菌异养植物中完整的线粒体基因组数据极为稀缺，在已知的十个进化分支中，仅有三个分支（杜鹃花科Ericaceae的Monotropa属、远志科Polygalaceae的Epirixanthes属以及兰科Orchidaceae的Gastrodia、Epipogium、Danxiaorchis属）完成了全基因组测序。为了填补这一空白，研究人员选取了隶属于薯蓣目（Dioscoreales）、水玉杯科（Thismiaceae）的狭口杯菌（Thismia angustimitra）作为研究对象。该物种是一种完全失去光合作用的菌异养植物，代表第四个尚未有线粒体基因组数据的进化分支。本研究旨在通过对其质体基因组和线粒体基因组的全面测序与分析，揭示其细胞器基因组的结构特征、基因含量、进化动态，并重点探究外源DNA的来源，特别是是否存在水平基因转移事件，从而加深对菌异养植物线粒体基因组进化规律的理解。
二、主要技术方法概述
研究人员采用Illumina HiSeq 2500平台对采自泰国汶干府普瓦野生动物保护区的一份硅胶干燥保存的T. angustimitra样本进行100 bp双端测序。质体基因组以T. tentaculata的质体基因组为参考进行组装和注释。线粒体基因组则利用GetOrganelle软件进行从头组装，并通过Bandage软件可视化评估组装结果。基因注释主要参考近缘种Trichopus zeylanicus的完整线粒体基因组以及T. rodwayi的部分基因序列，tRNA基因通过tRNAscan-SE 2.0检测。RNA编辑位点预测采用Deepred-Mt在线工具。重复序列分析使用Geneious Prime内置的重复查找器。为检测细胞内基因转移（intracellular gene transfer, IGT）和水平基因转移（HGT），研究人员将T. angustimitra的质体基因组与各线粒体染色体进行比对，并将所有线粒体基因间区序列与NCBI GenBank中的完整叶绿体基因组及全核苷酸数据库进行BLASTn相似性搜索。对于候选的外源插入片段，计算外星指数（Alien Index, AI）以区分其内源或外源属性，并对最长的psaA基因片段进行最大似然法系统发育分析以追溯其起源。此外，还专门针对真菌线粒体基因组进行了靶向BLAST搜索，以排除真菌来源的水平基因转移。
三、研究结果
质体基因组（plastome）与线粒体基因组（mitogenome）的特征
质体基因组被组装成一个18,215 bp的环状分子，平均测序深度为1,248×，具有标准的四分体结构，包含一个5,876 bp的反向重复区（inverted repeat, IR）。基因含量极低，仅含8个蛋白质编码基因、3个核糖体RNA（rRNA）基因和2个转运RNA（tRNA）基因。其中，rps12基因的三个外显子均存在，但其第二个内含子异常长，并完整包含了rrn5和rps2基因。
线粒体基因组被组装成22个环状染色体，总长度为422,400 bp，染色体大小从3,871 bp到94,341 bp不等。其中有5条染色体不含任何可识别的基因，但基于其与其它植物非编码线粒体序列的相似性、相似的测序深度（约235–489×）以及GC含量（41.2–43.1%），它们被确认为线粒体来源。基因补体包括39个蛋白质编码基因、3个rRNA基因和6个tRNA基因，其中rpl5为假基因。nad1、nad2和* nad5*三个基因为反式剪接（trans-spliced），其外显子分布在两至三条不同的染色体上。线粒体基因组中重复序列极少，仅占全基因组大小的0.8%，且重复介导的重组活性非常有限。唯一检测到参与重组的是一段99 bp的重复序列（Repeat 8），其替代构型仅获得8.6%的读段支持。
RNA编辑与外源DNA
RNA编辑预测结果显示，在T. angustimitra的38个潜在功能性蛋白质编码基因中共有499个编辑位点，低于自养近缘种Dioscorea oppositifolia的612个。在三者共享的27个基因中，Thismia的编辑位点总数（326个）也显著低于Dioscorea（476个）和Trichopus（385个），且未观察到可能指示反转录加工（retroprocessing）的C-to-T替换簇。
未检测到与真菌线粒体基因组相似的序列。线粒体基因组中质体来源的DNA含量很低，共鉴定出16个插入片段（MTPTs），总长4,068 bp，仅占线粒体基因组的0.96%。通过外星指数（AI）分析，其中5个片段为负值，提示可能源于细胞内基因转移（IGT）；另有9个片段为正值，提示可能源于水平基因转移（HGT），但仅有两个与漆树科（Anacardiaceae）最相似的片段AI值超过45。对其中最长的psaA片段（797 bp，AI=46）进行系统发育分析，结果强烈支持其嵌套在无患子目（Sapindales）的漆树科内部，与野漆（Toxicodendron vernicifluum）关系最近，支持率为100%（超快速自展值UB）。
进一步对全部线粒体基因间区序列进行全局BLAST搜索，发现155个长度超过200 bp的区域（总计87,450 bp，约占全基因组的21%）与其它植物线粒体基因间区序列高度相似。其中，绝大多数（119个）的最佳匹配序列来自广泛的真双子叶植物（eudicots），涉及16个目的31个科，尤以锦葵科（Malvaceae，26个）和列当科（Orobanchaceae，19个）最为突出。仅有9个片段的最佳匹配来自同属薯蓣目的物种。
系统发育分析
基于线粒体atp1基因、核18S rDNA和ITS区的联合系统发育分析，本研究的T. angustimitra标本以100%的支持率与另一份同种标本聚为一支，并与T. mirabilis、T. mucronata、T. okhaensis、T. puberula和T. nigricans共同构成一个支持率为90%的单系群。
四、讨论与结论
讨论总结：
研究人员指出，T. angustimitra高度退化的质体基因组完全符合菌异养植物质体基因组退化的普遍模式，其基因含量处于该属已知变异范围之内，且与同地采集的另一份标本近乎一致，仅存在一个碱基差异。关于线粒体基因组，其大小（422 kb）小于同目自养物种（如Trichopus 709 kb, Dioscorea 493 kb），不支持此前在Gastrodia中提出的“异养生活方式导致线粒体基因组扩张”的假说。相反，Thismia线粒体基因组中重复序列极少、重组活性极低的现象，与某些全寄生植物（如蛇菰科Balanophoraceae的Thonningia和Ombrophytum）类似。这可能归因于异养生活方式下，核编码的细胞器DNA复制、修复和重组机制功能失调，导致重复序列突变积累、同源性下降，进而抑制了重组。此外，多环状染色体结构的维持可能有助于减少突变负荷，甚至通过删除无基因的染色体来实现基因组精简。
在外源DNA方面，未发现真菌来源的序列，这与兰科菌异养植物中常见的古老真菌HGT事件形成对比。然而，研究人员发现了强有力的证据表明存在来自其它被子植物（尤其是真双子叶植物）的质体DNA通过HGT进入Thismia线粒体基因组。例如，一段psaA片段明确来源于漆树科，而大量基因间区序列与锦葵科、列当科等远缘类群的线粒体序列高度相似。这些发现挑战了以往认为菌异养植物中HGT罕见（主要局限于兰科）的观点，提示HGT在菌异养植物中可能比预期更为普遍。同时，基因间区序列的共线性模式暗示HGT可能双向发生。
研究结论翻译：
T. angustimitra高度退化的质体基因组完全符合异养植物质体基因组退化的普遍规律。尽管其线粒体基因组在许多方面看似不受生活方式影响，但其重组活性异常低下。令人惊讶的是，研究人员发现了暗示水平基因转移的证据，表明这种现象在菌异养植物中可能比之前认为的更为常见。
