突变压力与基因特异性选择塑造古老水生维管植物水韭属（Isoetes）叶绿体基因组密码子使用偏倚

《Ecology and Evolution》：Mutation Pressure and Gene-Specific Selection Shape Codon Usage Bias in the Chloroplast Genomes of Isoetes, an Ancient Aquatic Vascular Plant Lineage

　　

水韭属（Isoetes）是最古老的现存维管植物属之一，代表了石松类植物中的一个独特进化谱系。其叶绿体基因组保存了与早期陆地植物进化相关的分子特征。然而，塑造水韭属叶绿体基因组中密码子使用偏倚（Codon Usage Bias, CUB）的进化力量仍鲜有研究。研究人员分析了60个水韭属物种的叶绿体基因组，以刻画CUB模式，并厘清突变压力和自然选择对这一基部维管植物谱系叶绿体基因组进化的相对贡献。叶绿体基因组结构在物种间高度保守，基因组大小为142,880 bp至145,535 bp，GC含量为37.56%至38.18%。奇偶规则2（Parity Rule 2, PR2）图分析显示第三密码子位点存在显著的T/G偏倚，而中性图分析显示GC₃与GC₁₂之间无显著相关性。所有偏好密码子（RSCU > 1）均以A或U结尾，与AT结尾偏好一致。简单序列重复（Simple Sequence Repeat, SSR）分析显示，SSR主要分布于基因间区，其中单核苷酸重复占主导，表明其可作为系统发育研究的有价值分子标记。研究结果表明，水韭属叶绿体基因组表现出微弱但结构化的密码子使用偏倚。突变压力似乎是建立富含AT组成背景的主导力量，而自然选择以基因特异性的方式发挥作用，主要通过纯化选择和对特定叶绿体基因的功能约束。这些结果为基部维管植物的叶绿体基因组进化提供了新见解，并强调了突变偏倚与选择约束对密码子使用进化的共同作用。与以往大多数集中于被子植物的叶绿体CUB研究不同，这项对水韭属的属级全谱系分析为基部石松类谱系中的密码子使用进化提供了见解，并提示CUB模式、质体基因组保守性、功能约束和基因水平进化动态之间可能存在关联。

本论文针对古老水生维管植物水韭属（Isoetes）叶绿体基因组的密码子使用偏倚（Codon Usage Bias, CUB）进行了属级尺度的系统研究。水韭属是现存最古老的维管植物属之一，属于石松类谱系，化石记录可追溯至泥盆纪，被视为“活化石”。其叶绿体基因组保存了早期陆地植物进化的分子特征，但以往研究多集中于质体结构、系统发育、RNA编辑或少数代表物种，缺乏对CUB进化力量的系统评估。突变压力、自然选择与功能约束在该属叶绿体基因组进化中的相对贡献尚不清楚。为此，研究人员利用60个水韭属物种的叶绿体基因组，检验了三个假设：AT偏倚的突变压力塑造密码子使用背景；选择或翻译约束以基因特异性方式作用于光系统和核糖体蛋白基因；高度保守的质体结构与局部序列变异共存。研究结果表明，突变压力主导了富含AT的组成背景，而自然选择通过纯化选择对特定功能基因施加约束，二者共同塑造了微弱但结构化的CUB。该研究为基部维管植物叶绿体基因组进化提供了新见解，并凸显了突变偏倚与选择约束的协同作用。

本研究的关键技术方法包括：从美国国家生物技术信息中心（NCBI）数据库获取60个水韭属物种的完整叶绿体基因组序列，经严格质控后使用GeSeq重新注释并人工校正；计算密码子使用指标（ENC、GC₃s等）和相对同义密码子使用度（RSCU）；利用MISA-web鉴定简单序列重复（SSR）；通过IRscope比较反向重复（IR）边界；进行ENC-GC₃s图、奇偶规则2（PR2）图和中性图分析；使用R语言进行对应分析（COA）；基于31个共享叶绿体蛋白编码基因，以70个石松类分类群（60个水韭属物种和10个外类群）构建最大似然系统发育树（IQ-TREE）；通过GenoPlotR进行共线性分析；利用KaKs_Calculator计算Ka/Ks；采用mVISTA进行基因组相似性分析；并计算核苷酸多样性（π）、Watterson估计值（θ）和Tajima's D。

研究结果如下：

3.1 水韭叶绿体基因组特征：基因组大小为142,880 bp至145,535 bp，GC含量为37.56%至38.18%，结构高度保守，AT丰富；平均ENC为47.67，接近理论最大值61，表明整体CUB较弱。

3.2 RSCU值聚类分析：所有高RSCU密码子均以A或U结尾，亮氨酸密码子UUA强偏好，丙氨酸密码子GCU也表现出偏好，提示密码子使用模式在属内进化稳定。

3.3 重复序列分析：共鉴定5084个SSR，单核苷酸重复最丰富（占44.51%），主要分布于基因间区；保守型与物种特异型SSR共存，提示SSR是局部微变异的重要来源。

3.4 IR边界区比较：水韭属IR边界结构高度保守，ycf1和ycf2假基因在SSC/IRa边界稳定保留；与近缘属卷柏属（Selaginella）的明显边界扩张形成鲜明对比，表明水韭属质体结构演化策略保守。

3.5 ENC-GC₃s图分析：大多数基因位于高ENC、低GC₃s区域；22.91%的基因记录显著偏离Wright预期曲线，且主要集中于光系统II、光系统I和核糖体蛋白基因，提示这些功能模块受到更强的基因特异性约束。

3.6 PR2图分析：第三位点核苷酸呈现T/G偏倚，即T高于A、G高于C；A/T偏倚在基因组中相对稳定，而G/C偏倚在不同基因间变异更大，表明突变或链特异性过程作用于同义位点。

3.7 中性图分析：GC₃与GC₁₂回归斜率极低（-0.058至0.1668），相关性极弱，说明第三位点变异受突变压力驱动，而第一、二位点受氨基酸组成和蛋白质功能约束，二者去耦。

3.8 对应分析（COA）：前两个主轴仅解释总变异的9.71%-11.32%和7.45%-9.64%，基因在二维空间中连续分布，无强分化，提示主要密码子使用趋势整体一致，但变异分散于多个维度。

3.9 系统发育和基因组共线性分析：最大似然树支持水韭目（Isoetales）与卷柏目（Selaginellales）为姐妹群，水韭属为单系，但种间分支短、支持率有限；33个代表物种的质体基因组共线性高度保守，基因顺序和方向一致。

3.10 叶绿体基因选择压力模式：大多数基因的Ka/Ks远小于1，纯化选择主导；ndhC等NDH复合体基因出现相对升高的Ka/Ks（均值1.81），可能反映松弛选择或谱系特异速率变化，而非确定性正选择。

3.11 比较基因组相似性分析：以日本水韭（I. japonica）为参考，差异主要集中在非编码区；ndhB、rps2、atpF编码区分歧相对较高；I. nuttallii、云贵水韭（I. yunguiensis）和I. duriei与参考差异最大，而I. baodongii和I. orientalis极为保守。

3.12 核苷酸多样性分析：81个共享蛋白编码基因的平均π为0.0133，SSC区多样性最高，IR区最低；atpF、ycf1、ccsA、matK、ndhF和rpl32为高变区；76.5%的基因Tajima's D为负值，提示纯化选择或群体扩张迹象。

讨论部分总结了CUB的进化力量：突变压力建立AT富集背景，选择约束集中于光系统和核糖体蛋白基因，SSR分布与功能约束相关，NDH基因的异常Ka/Ks可能反映松弛选择而非正选择。结论部分指出，水韭属叶绿体基因组在基因组大小、GC含量、IR结构和基因组织上高度保守，CUB整体较弱但结构化，所有偏好密码子以A或U结尾。多项证据表明，AT偏倚的突变压力是密码子使用背景的主导因素，而基因特异性约束进一步塑造特定功能类别。Ka/Ks分析显示纯化选择普遍存在，部分NDH基因升高的Ka/Ks可能反映松弛选择。本研究鉴定的SSR、高变区和密码子使用特征为水韭属的系统发育和保护研究提供了基因组资源，但生态学解释仍需结合环境数据和功能验证。