海洋硫氧化内共生菌转座子含量与基因组大小及宿主分类的关联性研究

《mSystems》：Transposable element content is associated with total genome size and host taxonomy in sulfur-oxidizing bacterial endosymbionts

　　

摘要 转座子（Transposable elements, TEs）是可移动的遗传序列，能够破坏、重定位或创造基因及调控元件，从而促进基因组进化。其在塑造适应性（包括共生关系）中的潜在作用日益受到关注。迄今为止，大多数关于TEs的研究集中于陆地环境，而对海洋生境中TEs的功能和分布知之甚少。本研究调查了海洋化能合成共生体系中TEs的存在与分布，特别是与生活在多种化学还原环境中的动物相关的硫氧化γ-变形菌纲（Gammaproteobacteria）基因组。研究数据表明，系统发育关系和生境类型与TE丰度无显著关联，而基因组大小、宿主科和生活方式等因素则与之显著相关。几乎所有囊蛤科（Vesicomyidae）蛤类共生菌缺乏TEs，这可能归因于还原性基因组进化；而贻贝科（Mytilidae）共生菌具有最高的TE丰度，可能与其同源重组率增加或与适应性基因组精简相关的TE增殖有关。转座酶家族多样性在各共生菌类群间差异显著，但与基因组大小无关，表明其他因素影响转座酶组成。值得注意的是，大部分转座酶基因侧翼为参与DNA复制与修复、膜转运及原核防御系统的基因，暗示其在环境或宿主适应中的潜在作用。本研究揭示了共生细菌中TE含量的动态特性，并为TEs在多样化学还原环境中细菌进化中的功能假说提供了依据。

论文解读：海洋硫氧化内共生菌转座子含量与基因组大小及宿主分类的关联性研究

一、研究背景与科学问题

转座子（Transposable elements, TEs）是一类可在基因组内移动的DNA序列，其通过转座酶实现自我切除与重新整合。传统上被视为"垃圾DNA"的TEs，如今被广泛认为是细菌基因组进化与环境适应的重要驱动因素，能够通过破坏、重定位或创造基因及调控元件来促进基因组进化。尽管陆地环境中TEs与细菌适应性之间的因果关系已有充分证据，但海洋微生物中TEs的基因组丰度与生态功能研究仍相对匮乏。特别是在深海热液喷口、烃类冷泉及富含有机质的硫化物泥滩等化学还原环境中，化能合成共生体系内TEs的分布规律与生态意义几乎未被系统探讨。已有研究表明，垂直传播的共生菌因缺乏自由生活阶段及重组机会减少，通常携带较少的TEs；而水平传播的共生菌则因环境暴露而可能保有较高的TE含量。然而，针对深海化能合成内共生菌的跨类群比较基因组研究此前尚属空白。本研究旨在系统评估海洋硫氧化γ-变形菌纲内共生菌基因组中TEs的丰度、多样性及其影响因素，填补海洋共生体系中TEs研究的空白。

二、研究内容与主要结论

研究人员收集了38个海洋硫氧化γ-变形菌纲内共生菌基因组及36个自由生活近缘菌基因组，涵盖囊蛤科（Vesicomyidae）、贻贝科（Mytilidae）、西伯加虫科（Siboglinidae）、索伦蛤科（Solemyidae）和露辛蛤科（Lucinidae）五个宿主科。通过转座酶基因检测、TE全基因组BLAST比对、转座酶家族组成分析、系统发育广义最小二乘（PGLS）模型及侧翼基因功能注释等手段，研究人员发现：宿主科与TE丰度存在显著关联，贻贝科共生菌TE含量最高（平均1.529%），囊蛤科共生菌几乎完全缺乏TEs（平均0.046%）；基因组大小与TE丰度在共生菌和自由生活菌中均呈显著正相关，但转座酶多样性与基因组大小的关联仅存在于自由生活菌中；生活方式显著影响TE含量，自由生活菌的TE百分比显著高于共生菌；生境类型（热液喷口vs冷泉）对TE丰度无显著影响，但在排除垂直传播的囊蛤科共生菌后，深水共生菌的TE含量显著高于浅水种类。侧翼基因功能分析显示，转座酶侧翼基因主要富集于DNA复制与修复、膜转运及原核防御系统等KEGG通路，暗示TEs可能通过携带这些功能基因参与宿主或环境适应。

三、关键技术方法

研究人员从NCBI微生物基因组数据库获取38个共生菌基因组及对应注释文件，并从RefSeq数据库下载36个自由生活硫氧化菌基因组作为对照。转座酶基因检测采用grep搜索注释文件结合bedtools合并计算基因组百分比，同时利用BLAST比对TnCentral、ISFinder和Integrall数据库进行TE全基因组鉴定（序列一致性≥80%，e值≤1e–10）。转座酶家族组成通过解析GFF注释文件的product字段并结合InterProScan对未分类序列进行功能预测。系统发育树使用GToTree v1.8.12基于74个保守单拷贝基因构建，PGLS模型采用Martins和Hansen模型控制系统发育信号。侧翼基因功能注释采用KofamScan工具比对KOfam数据库进行KEGG通路分配。所有统计分析均在R v4.4.3中完成。

四、研究结果

1. 化能合成共生菌的TE含量因宿主科而异

Kruskal-Wallis检验显示宿主科与TE丰度显著相关（P=0.000102）。两两比较表明，贻贝科共生菌TE百分比显著高于其他所有宿主科（中位数1.3%），而囊蛤科共生菌显著低于其他类群（中位数0%）。转座酶假基因含量的模式与之类似，但贻贝科与露辛蛤科之间无显著差异。转座酶多样性方面，贻贝科（Shannon指数1.56±0.342）与其他水平传播共生菌相近，但含垂直传播成分的索伦蛤科多样性显著较低（0.879±0.259）。IS3家族转座酶在所有宿主科中占主导（11%–52%），贻贝科富集IS5（23%）和IS4（14%）家族，索伦蛤科则以IS30家族为主（31%）。

2. 基因组大小影响TE含量与多样性

TE百分比与基因组大小在共生菌（P=0.0274）、自由生活菌（P=0.0000403）及合并数据集（P=0.000165）中均显著正相关，自由生活菌的关联强度（ρ=0.640）高于共生菌（ρ=0.358）。Shannon多样性指数与基因组大小仅在自由生活菌中显著相关（ρ=0.802，P=4.207e–09），共生菌中无此关联（ρ=0.0452，P=0.816）。

3. 细菌生活方式影响TE丰度

TE百分比与水深在全部共生菌中无显著相关（P=0.716），但排除囊蛤科后呈显著正相关（P=0.0165，ρ=0.442），仅含水平传播共生菌时结果类似（P=0.0348）。热液喷口与冷泉生境间TE丰度无显著差异（P=0.922）。自由生活菌TE百分比（平均1.08%）显著高于共生菌（平均0.61%，P=0.000699），排除囊蛤科后差异仍显著（P=0.0396）。但转座酶假基因百分比在两类群间无显著差异。

4. 系统发育广义最小二乘分析

PGLS模型显示，囊蛤科共生菌的TE百分比显著低于参考组露辛蛤科（P=0.0459）。生活方式对TE百分比的影响显著（P=0.0405），但加入基因组大小后，基因组大小解释了大部分变异（P=0.00160）。TE百分比的系统发育信号较弱（Blomberg's K=0.0130–0.00555），Pagel's lambda为中等水平（0.451–0.467）。

5. TE侧翼基因功能

侧翼基因功能注释显示，"复制与修复"通路占比最高（共生菌28.6%，自由生活菌24.3%），其次为"转运蛋白"（8.8% vs 8.4%）、"原核防御系统"（3.6% vs 7.8%）和"双组分系统"（4.8% vs 5.4%）。

五、讨论与结论

本研究表明，海洋硫氧化内共生菌的TE含量与基因组大小和宿主分类显著相关，而与系统发育关系和生境类型无明确关联。自由生活菌TE丰度高于共生菌，反映了其在更大基因组和更多重组机会下对TEs的较低选择压力。囊蛤科共生菌几乎完全缺乏TEs，符合垂直传播共生菌还原性基因组进化的特征。贻贝科共生菌的高TE含量可能与其高重组率及适应性基因组精简过程相关。深水共生菌（排除囊蛤科后）TE含量较高，可能反映了TEs在应对深海环境胁迫中的适应性作用。转座酶多样性与丰度在共生菌中缺乏关联，提示不同宿主生态位对转座酶家族存在差异性选择。侧翼基因中复制修复、膜转运和防御系统基因的富集，暗示TEs可能通过携带这些功能基因参与宿主或环境适应。转座酶基因在复制起点附近的聚集模式，提示其可能参与复制相关的基因组结构重排。总体而言，TEs的多样性、密度和分布为理解宿主-共生菌互作、微生物自由生活阶段的适应机制及微生物谱系进化提供了重要视角。