荨麻科楼梯草族（Elatostemateae）的质体系统发育基因组学：结构变异与系统发育关系的洞见

《Ecology and Evolution》：Plastid Phylogenomics of Elatostemateae (Urticaceae): Insights Into Structural Variation and Phylogenetic Relationships

　　

专业术语翻译摘要 楼梯草族（Elatostemateae）包含约1300个物种，是荨麻科（Urticaceae）内物种最丰富的类群，然而其质体组（plastome）进化及多个属间关系仍未完全解决。研究人员分析了代表楼梯草族现存主要谱系的40个物种，使用了完整质体组和核糖体DNA（nrDNA）；由于缺乏合适的DNA材料，单种属Metapilea被排除在外。所有质体组均保留了典型的被子植物四分体结构，长度范围为143,738至163,969 bp。比较分析揭示了大单拷贝区（LSC）内的大片段倒位、显著的反向重复区（IR）边界移位，以及罕见的源自线粒体DNA的质体插入序列（PTMT）。七个高变区被鉴定为潜在的DNA条形码候选标记。基于质体组、蛋白质编码序列（CDS）、nrDNA以及联合数据集的分析均恢复了相同的三个楼梯草族主要分支，尽管nrDNA树中的数个内部节点支持率较低。质体与nrDNA数据将Gyrotaenia置于楼梯草族内的广义冷水花属（Pilea s.l.）的姐妹群位置，而近期基于低拷贝核基因的证据则将Gyrotaenia + Myriocarpa置于核心楼梯草族的姐妹群位置，揭示了细胞核质不一致性（cytonuclear discordance）。这些结果为楼梯草族的系统学研究提供了一个框架，并记录了显著的质体组结构变异；未来需要更密集的低拷贝核基因或全基因组采样，以区分不完全谱系分选（incomplete lineage sorting）与历史基因渗入（historical introgression）。

论文解读文章

研究背景与科学问题

荨麻科（Urticaceae）下的楼梯草族（Elatostemateae）是一个泛热带分布的类群，拥有约1300个物种，是该科内物种多样性最高的族。该族植物形态多样，从纤细草本到木质灌木，通常栖息于阴湿环境。其分类鉴定依赖于笔状柱头和由退化雄蕊介导的特殊瘦果弹射机制等生殖性状组合，但许多形态特征存在同塑性，这促使研究人员不断利用分子手段重新评估属级界限。早期的系统发育重建多依赖有限的Sanger测序基因集，常导致分辨率不足或系统发育假说相互矛盾。例如，全球尺度采样虽验证了楼梯草族的核心单系性，却同时揭示了Pellionia属的多系性。此外，对近缘族Urticeae的研究确认了其单系性，并将Gyrotaenia属定位为楼梯草族的可能姐妹群。后续利用靶向富集系统发育组学和限制性位点数据的研究，则恢复了一个包含Gyrotaenia和Myriocarpa的分支，作为核心楼梯草族的姐妹群。尽管物种关系逐步澄清，但对楼梯草族内质体组进化的综合性研究仍几乎空白，这成为理解该族更深层次进化历史的瓶颈。质体组因其相对保守的结构、母系遗传特性以及高通量测序的高可及性，是解析复杂进化轨迹的有效工具。本研究旨在整合全质体组架构与核糖体DNA（nrDNA）序列数据，以解决未明的属间关系并填补比较质体组学框架的缺失。具体目标包括：（1）提高族内各属及主要形态分支的系统发育关系分辨率；（2）刻画质体组结构变异，包括反向重复区（IR）边界移位和重排，并评估其系统发育与分类学意义；（3）鉴定用于未来物种鉴定和DNA条形码的候选高变叶绿体DNA（cpDNA）区域。

研究方法概述

研究人员采集了40个代表楼梯草族现存主要属的物种（不包括可能已灭绝且无法获得合适DNA材料的单种属Metapilea），覆盖了多样性最高、分类最复杂的Elatostema和Pilea两个属内的主要形态分支和公认的组。总基因组DNA从硅胶干燥的叶片组织中提取。文库构建和测序由上海美吉生物医药科技有限公司完成，在Illumina HiSeq 4000平台上生成2 × 150 bp双端读长。原始读长经过滤后，使用GetOrganelle v1.7.6.1进行从头组装，不依赖参考质体组。完整的质体组分别使用PGA和CPGAVAS2进行独立注释以验证基因边界。系统发育分析涵盖了74个完整质体组和74条对应的nrDNA序列，矩阵包括40个楼梯草族内类群、其他荨麻科谱系的代表以及4个桑科（Moraceae）外类群物种。研究人员采用了四个数据分区进行分析：完整质体组（CP）、质体蛋白质编码序列（CDS）、nrDNA区域以及串联的CP+nrDNA数据集。系统发育推断使用贝叶斯推断（BI）和最大似然法（ML）。序列比对使用MAFFT v7.407，并用trimAl v1.4.1去除比对不佳或富含空位的列。最佳核苷酸替换模型分别由jModelTest2 v2.1.7和IQ-TREE v1.6.12选择。BI分析在MrBayes v3.2.7a中进行，运行两个独立的马尔可夫链蒙特卡洛链各100百万代。ML分析在IQ-TREE v1.6.12中进行，使用5000次标准非参数自举复制估计分支支持率。质体组结构变异分析方面，研究人员在Geneious中汇总了总质体组大小、GC含量、基因含量及各区域长度，并使用CPJSDRAW可视化LSC/IRb、IRb/SSC等连接处。大规模重排通过Mauve软件评估。重复序列分为长重复和简单重复序列（SSRs），分别由REPuter和MISA鉴定。核苷酸多样性（Pi）在DnaSP v6中计算。

主要研究结果

3.1 质体组组织与特征

质体组大小从Pilea micromeriifolia的143,738 bp到Elatostema oblongifolium的163,969 bp不等。所有质体组均保留典型四分体结构，包含两个长度为18,361–31,485 bp的IR区，分隔了LSC区（80,747–89,128 bp）和SSC区（12,924–18,451 bp）。整体基因组GC含量稳定在36.23%–36.87%之间，IR区GC含量高于LSC和SSC区。基因注释检测到127–133个不同基因，包括82–88个蛋白质编码基因（PCGs）、37个tRNA基因和8个rRNA基因。结构分析确认有15个基因含单个内含子，2个基因（clpP和ycf3）含有两个内含子。

3.2 系统发育模式与细胞核质一致性

对CDS、CP、nrDNA及串联CP+nrDNA数据集的分析均恢复了相同的三个主要楼梯草族分支和相同的主干拓扑结构，但不同数据集的支持率有所差异。例如，分支A在CP和nrDNA分析中分别获得98.3/1.0和98.1/1.0的自展值/后验概率支持。nrDNA树中一些内部节点支持率较低，如Elatostema内为78.7/0.97，Pilea内为84.9/1.0。在科级框架下，Urticeae和Elatostemateae互为姐妹谱系。在楼梯草族内，三个主要分支为：广义楼梯草属（Elatostema s.l.）（分支A）、Gyrotaenia（分支B）和广义冷水花属（Pilea s.l.）（分支C）。分支A内包含四个亚支（A1–A4），分别对应Polychroa、Procris、Elatostematoides和Elatostema分类群。分支C作为Gyrotaenia（分支B）的姐妹群被恢复，其内三个亚支（C1–C3）对应Achudemia、Lecanthus和Pilea。

3.3 质体组重排、IR动态与线粒体整合

全质体组比对显示大多数质体组共线性，但存在若干谱系特异性结构变异。Elatostema parvum和E. pycnodontum在LSC区内存在一个约20 kb的倒位，从rpoB延伸至trnQ-UUG，涵盖约14个基因。这两个携带倒位的分类单元形成亚支A4c。IR连接处的评估揭示了IR扩张和收缩程度的显著变异。在Elatostema nanchuanense和E. oblongifolium的IR区中，ycf15与trnV-GAC之间检测到一个3941–7033 bp的插入，导致IR区显著增长。该插入序列经BLASTn搜索验证，被认为是源自线粒体DNA的质体序列（PTMT）。Pilea tetraphylla表现出最显著的不对称IR边界移位，IRb收缩使边界移向IR并增加了LSC长度，而JSB处的扩张则减少了SSC长度。Pilea micromeriifolia则表现出明显的IR收缩，将ycf2等多个基因从IR移至LSC。

3.4 重复序列与核苷酸多样性

在40个质体组中鉴定出5477个重复元件，包括3042个SSRs、1379个分散长重复和1056个串联长重复。单核苷酸SSR最为常见，占66.41%，其中poly-A和poly-T基序占主导。核苷酸多样性分析表明，基因间隔区普遍比编码基因更具变异性。在100个基因区域中，ycf1、rpl22和rps16的Pi值大于0.08。在95个基因间隔区中，rps15-ycf1、trnH-GUG-psbA、rpl32-trnL-UAG和ndhF-rpl32的Pi值大于0.16，这些区域被确定为标记开发的候选区域。

讨论与结论

本研究通过整合完整质体组和nrDNA数据，为楼梯草族建立了一个多位点系统发育基因组学框架。质体和nrDNA分析恢复了相同的主要分支，但nrDNA数个内部节点支持率降低，表明一致性主要体现在主要分支水平而非每个内部节点。Gyrotaenia的放置是最清晰的细胞核质冲突点。本研究的质体和nrDNA分析将其置于楼梯草族内，作为广义冷水花属的姐妹群，这与早期多位点结果一致。相反，基于Angiosperms353低拷贝核基因的研究则将Gyrotaenia与Myriocarpa一起置于核心楼梯草族的姐妹群位置。这种冲突可以用多种进化过程解释，包括快速分化过程中的不完全谱系分选（ILS）和历史基因渗入或叶绿体捕获，但现有质体组加nrDNA数据无法区分它们，需要更密集的低拷贝核基因座或全基因组数据来检验。在结构进化方面，楼梯草族质体组显示出显著的IR边界变异，尤其是Pilea tetraphylla表现出的不对称IR收缩与扩张组合。在Elatostema pycnodontum和E. parvum中检测到约20 kb的大片段倒位，这是该族中的首次记录。另一个显著特征是PTMT插入，反映了罕见的历史重组事件，可能是在修复质体组双链断裂过程中整合了线粒体DNA。主要的质体组结构变异在系统发育树上呈聚类分布，并与不同的地理区域相关联，但现有采样无法证明地理隔离或生态条件导致了这些重排。

研究结论部分翻译：本研究通过整合跨越广泛分类采样的完整质体组和nrDNA，为楼梯草族建立了一个多位点系统发育基因组学框架。这些分析解析了所采样的大部分主要分支，并支持了近期的一些分类学调整，而Gyrotaenia的放置突显了一个可能反映不完全谱系分选或历史基因渗入的细胞核质冲突。比较质体组分析进一步揭示了大的倒位、显著的IR边界移位和一个罕见的PTMT插入，证明了该族内显著的结构变异。七个高变区，包括ycf1，是用于DNA条形码和物种水平研究的有前景的候选区域。总之，这些结果澄清了所采样谱系的进化历史，同时将核基因组采样和更广泛的群体水平比较确定为未来工作的优先方向。