系统发育与生态系统共同塑造原核微生物组中结构化非编码RNA介导的代谢物感知

《Genome》：Phylogeny and ecosystem jointly shape metabolite sensing by structured noncoding RNAs

非编码RNA（ncRNAs）通过多样的结构和序列依赖性机制调控基因表达，且常无需蛋白质中间体参与。核糖开关（riboswitches）是直接感知小分子代谢物以控制基因表达的顺式调控非编码RNA（cis-regulatory ncRNAs）。许多核糖开关通过动力学控制（kinetic control）发挥作用，即配体（ligand）与适配体结构域（aptamer domain）结合的速度决定基因表达被开启（ON）还是关闭（OFF）。在此，研究人员对八类微生物组中的核糖开关谱进行了比较调查，包括人类肠道、口腔、皮肤、阴道、小鼠肠道、猪肠道、羊瘤胃和海洋环境，覆盖13 946个原核基因组。利用协方差模型（covariance models），研究人员在12 916个基因组（92.6%）中鉴定出86 485个核糖开关，代表23个代谢物感知类别，其中6个优势类别约占全部核糖开关命中的70%。核糖开关组成在生态系统间差异显著，但方差分解（variance partitioning）显示该结构主要由系统发育（phylogeny）驱动：门（phylum）解释了25.5%的组成方差，而在控制分类学（taxonomy）后，独属于生态系统的方差比例仅为2.2%。总体而言，核糖开关分布反映了系统发育和生态学的双重影响，突出其在塑造微生物生命代谢物响应性调控策略中的作用。

研究背景：非编码RNA（ncRNAs，不翻译为蛋白质的RNA）可通过结构性和序列依赖性机制调控基因表达。核糖开关（riboswitches）是位于mRNA 5′-UTR中的顺式调控RNA，由配体结合适配体结构域（aptamer domain，配体感应区）和相邻表达平台（expression platform，将配体结合转化为表达调控的结构区）组成；配体结合后可通过提前终止转录、阻断翻译或促进mRNA降解等方式调控基因表达。已知核糖开关可感知氨基酸、辅因子、离子、核苷酸衍生物和信号分子等，已报道近55类，并存在大量配体未明的孤儿核糖开关。尽管单数据集研究已揭示其分类分布，但跨多种微生物组的核糖开关谱组成究竟主要由系统发育继承还是生态系统选择塑造，以及在控制分类学后生态签名是否仍保留，尚缺乏系统比较。该研究因此比较八类微生物组，以解析系统发育与生态对代谢物感知策略的相对贡献，并为后续实验验证提供可检验假设。

研究人员整合8个微生物组、13 946个物种级代表基因组，开展跨生态系统的核糖开关普查，共鉴定86 485个核糖开关，覆盖12 916个基因组（92.6%）和23个代谢物感知类别。结论是核糖开关组成在生态系统间显著不同，但该结构主要由系统发育驱动；控制门后生态系统仍有小但显著贡献。总体分布反映系统发育和生态双重影响，突出结构化非编码RNA在微生物代谢响应性调控策略中的作用。

技术方法（≤250字）：整合MGnify参考目录：人肠道4744（UHGG v2.0.2）、口腔452、皮肤579、阴道280，小鼠肠道2847、猪肠道1376、羊瘤胃2172、海洋1496，共13 946个物种级代表基因组；分类采用GTDB r220并设质量阈值。用Infernal v1.1.5 cmscan和Rfam 15.0的42个协方差模型（CMs，同时建模RNA序列与二级结构）按gathering threshold扫描。以23代谢物类别和41家族统计，采用二值流行率控制多家族偏差；用PERMANOVA（置换多元方差分析）、方差分解、Fisher精确检验、系统发育信号、NRSI（Niche-Riboswitch Specificity Index，生态位-核糖开关特异性指数）、共现网络、PCoA（主坐标分析）/Shannon多样性和bootstrap敏感性分析。

研究结果：（1）核糖开关在八类微生物组中近乎普遍：对13 946个基因组扫描检测到86 485个核糖开关，12 916个基因组（92.6%）至少含一个；检测率81.9%—96.2%，校正基因组完整性和系统发育后生态仍显著；六类占69.5%，钴胺素（cobalamin）命中最多，硫胺素焦磷酸（TPP）基因组流行率最高；人肠和皮肤负担高、海洋最低。（2）核糖开关组成跨生态系统显著不同：PERMANOVA显示组成差异显著（R²=0.060；二值Jaccard R²=0.055），Shannon多样性人肠最高、海洋最低，PCoA分离；海洋钴胺素比例高，口腔和阴道TPP比例高，肠更均衡。（3）门水平富集揭示不同代谢策略：对6门×23类做Fisher检验，124/138显著；Bacillota_A富集四氢叶酸（THF）、镍/钴（Ni/Co）、c-di-GMP等而贫锰/SAH（S-腺苷同型半胱氨酸），Bacteroidota多数贫化但富钴胺素，Pseudomonadota富SAH等，其他门各有签名。（4）区分系统发育与生态系统：方差分解显示门独有25.5%，控制门后生态仅2.2%；但10个门内生态效应均显著，中位R²=0.118；Pagel λ中位0.913、Fritz–Purvis D中位0.385，41家族去分组后生态调整R²=0.020，结论为系统发育主导、生态小幅显著。（5）样本量敏感性：多bin重抽样bootstrap显示PERMANOVA R²围绕0.060、Kruskal ε²约0.10，Fisher结果保持显著；仅人/鼠肠在2000 bin因组成相似出现例外，说明数据量不平衡未改变主要结论。（6）代谢物感知核糖开关显示生态位特异分布：NRSI显示人肠道7类、皮肤6类、阴道和海洋各3类等最高；TPP对皮肤、赖氨酸对人肠、黄素单核苷酸（FMN）对阴道NRSI最高，锰对皮肤、谷氨酰胺和SAH对海洋、孤儿对人肠；但NRSI需结合系统发育解释。（7）共现网络揭示模块化代谢感知架构：二值矩阵Jaccard>0.1构建55节点135边网络，Leiden识别7社区；S-腺苷甲硫氨酸（SAM）为高连接枢纽，镁具最高中介中心性，FMN-赖氨酸等共现强，c-di-AMP（环二腺苷酸）和氨基糖苷位于低度外围。（8）人身体部位核糖开关谱不同：5708个人类相关基因组有命中，肠/皮肤负担高于口腔/阴道；TPP各部位最流行，钴胺素在肠高阴道低，c-di-GMP（环二鸟苷酸）肠高阴道低，皮肤富甘氨酸、锰和c-di-AMP。

讨论与结论：从生态视角看，一个基因组的核糖开关集合如同代谢签名，揭示生物对环境的适应；营养丰富的宿主相关环境，尤其哺乳动物肠道，倾向于选择具有更广泛代谢物感知调控网络的生物。辅因子感知核糖开关占主导（钴胺素、TPP、SAM、FMN合计56%命中），并在八个生态系统中保持，提示为细菌调控的普遍特征。TPP核糖开关广泛分布支持其可能为最早演化的核糖开关功能之一；海洋基因组核糖开关多样性较低可能反映贫营养下基因组简化。系统发育与生态通过不同机制塑造核糖开关谱，门水平富集反映谱系特异调控策略，方差分解表明系统发育主导，而NRSI捕获栖息地偏好且控制系统发育后仍有小但显著的生态效应；核糖开关类别缺失常反映替代调控机制而非通路缺失。共现网络支持核心感知模块和低连接外围节点；人体部位负担差异可能反映代谢条件差异，皮肤甘氨酸和锰核糖开关富集提示生态位特异代谢需求。研究优势包括物种级代表基因组去冗余、Rfam family-specific gathering thresholds、二值流行率控制多家族偏差并同时报告效应量；局限包括每种仅一个代表遗漏株水平变异、序列检测不等同功能验证、单一CM可能涵盖不同配体亚类、数据量不等仍可能影响部分统计、横断面设计无法完全分离系统发育约束与环境选择。结论部分译述为：核糖开关集合可作为代谢适应签名；营养丰富的宿主相关环境选择更广泛的代谢物感知网络；辅因子感知核糖开关在八类生态系统中占主导；系统发育是组成的主要驱动因素，而生态贡献较小但显著；这些模式为优先实验验证特定核糖开关-环境关联及界定生态位特异非编码RNA调控提供基础。