牛链球菌（Streptococcus uberis，S. uberis）是牛乳腺炎的主要病原体，然而，促进该病原体在不同宿主相关环境或选择压力下适应的基因组机制仍知之甚少。本研究分析了来自泰国三个奶牛场的全基因组序列数据，以探究重组和可移动遗传元件（MGEs）对S. uberis进化和适应的贡献。在138个S. uberis基因组中，鉴定出42个核心基因组序列类型（cgSTs），并频繁检测到MGEs，如质粒相关基因（81.1%的分离株）、前噬菌体（67.4%的分离株）和插入序列（26.1%的分离株）。来自农场A的分离株表现出最长的重组片段大小，但重组频率和重组与突变之比极低。相比之下，来自农场B的分离株具有最高的抗菌药物耐药性（AMR）基因流行率，表现出较高的重组与突变之比（R/θ = 4.42），并且更频繁地携带与AMR基因相关的MGEs。最后，来自农场C的分离株共享单一核心基因组和AMR谱，但携带有多种前噬菌体。多个前噬菌体与感染其他细菌属的噬菌体具有高度序列相似性（>99%），提示细菌物种间的生态重叠可能促进跨属基因交换，强调了微生物生态学对S. uberis进化的影响。总体而言，研究人员的结果说明了在奶牛养殖环境中促进S. uberis适应性进化的多种机制和遗传元件。
论文解读文章
**研究背景**
乳腺炎是奶牛养殖业中导致重大经济损失的疾病，表现为产奶量和质量下降、治疗成本增加以及患病动物淘汰。牛链球菌（Streptococcus uberis，S. uberis）是引起临床和亚临床牛乳腺炎的主要病原体之一，被归类为环境和接触性病原体，能在胃肠道、粪便、垫料、牧场、牛皮肤和乳腺等多种环境中存活。其显著的遗传多样性使其能够适应不同生态位，但这种适应性的基因组机制，特别是在不同宿主相关环境或选择压力下，仍不明确。全基因组测序已揭示S. uberis通过突变、重组和可移动遗传元件（MGEs）获取外源DNA来快速响应选择性压力，如宿主免疫、抗菌药物使用和环境条件。然而，在不同养殖环境或高抗菌药物选择压力下，这些过程如何表现尚不清楚，尤其是热带地区（如泰国）的奶牛场。此外，S. uberis中噬菌体的研究也十分有限。因此，本研究旨在探讨重组、MGEs及其与抗菌药物耐药性（AMR）基因的关联在S. uberis适应性进化中的作用，通过全基因组测序和比较基因组分析，解析泰国S. uberis分离株中重组事件和与AMR基因水平转移相关的MGEs，并探究MGEs（特别是前噬菌体）如何促进S. uberis的进化多样性和适应。该论文发表在《Veterinary Research》。
**主要关键技术方法**
研究人员从泰国三个奶牛场（农场A：2013年，29株；农场B：2019年，53株；农场C：2023年，56株）收集了138个S. uberis全基因组序列。关键方法包括：①核心基因组多位点序列分型（cgMLST）和多位点序列分型（MLST），通过PubMLST数据库进行等位基因谱分配和系统发育树构建；②重组检测，使用ClonalFrameML基于核心单核苷酸多态性（SNP）比对估算重组参数（重组与突变比R/θ、重组片段平均长度等），并用SplitsTree构建邻接网络图；③可移动遗传元件鉴定，采用Mobile Element Finder筛选MGEs（质粒相关基因、插入序列等），PHASTEST识别前噬菌体区域，并利用Phagonaute和Virus-Host DB进行宿主关联分析；④抗菌药物耐药基因检测，结合Mobile Element Finder和先前发表的数据集。所有分析均基于短读长Illumina组装基因组。
**研究结果**
**基于序列分型的核心基因组分析（Core genome analysis based on sequence typing）**
通过cgMLST和MLST分析，在138个基因组中鉴定出6个序列类型（STs）和42个核心基因组序列类型（cgSTs）。其中，ST2048是最常见的ST（27株），主要分布于农场A和B。农场C的分离株因等位基因缺失无法分配ST，但cgSTs显示其多样性较低。基于ST的Simpson's discrimination index（SID）在农场B最高（0.759），农场A中等（0.431），农场C为0（因无法分型）。基于AMR基因模式的SID在农场B最高（0.830），而基于AMR-MGE联合模式的SID在所有农场均较高（0.833–0.956），其中农场C高达0.956，表明MGEs分布驱动了显著变异，尽管核心基因组异质性低。系统发育分析显示，泰国分离株分散在全球多个谱系中，且按农场来源聚类，同时观察到克隆内替换和进化（如同一奶牛不同时间点分离株的附属基因组差异）。
**重组分析（Recombination analysis）**
ClonalFrameML分析揭示不同农场重组动态差异显著：农场B的重组与突变比（R/θ = 4.42）远高于农场A（0.00007）和农场C（0.00294），表明农场B中重组对遗传多样性的贡献远大于点突变。农场B的导入DNA分歧度（D = 0.238）也较高，与更大的基因组多样性和更广的MGEs及AMR基因谱一致。农场A的重组片段最长（约1000 bp），但频率极低；农场C重组频率低，但MGE多样性高，提示MGEs在遗传变异中起更大作用。邻接网络分析进一步支持这些发现，农场B分离株呈现网状连接，表明存在重组或水平基因转移信号，而农场A和C的分离株形成紧密簇。
**可移动遗传元件（Mobile genetic elements）**
共检测到三种插入序列（IS），其中ISSsu5最常见（32株）。35株分离株的AMR基因与MGEs位于同一重叠群，多数为IS。质粒相关基因（如repUS55）在农场C中普遍存在（41株中39株）。前噬菌体分析显示，农场A、B和C分别有16株（55%）、27株（51%）和51株（91%）携带前噬菌体。农场A仅有两种前噬菌体类型，而农场B和C各有11种和10种。尽管农场C的分离株共享单一核心基因组和AMR谱，但前噬菌体多样性高，表明噬菌体可驱动核心基因组稳定的种群内的微进化。多个前噬菌体与感染其他细菌属（如化脓链球菌、嗜热链球菌、乳酸乳球菌、粪肠球菌）的噬菌体具有高度序列相似性（>99%），提示跨属基因交换。此外，一些前噬菌体携带宿主样酶（如trans-prenyltransferase, mannosyl-glycoprotein endo-β-N-acetylglucosaminidase, alkylhydroperoxide reductase），可能参与宿主适应。
**讨论与结论总结**
讨论部分指出，本研究揭示了MGEs在S. uberis适应和进化中的核心作用，尤其是AMR基因与MGEs的频繁共定位，表明抗菌药物使用（AMU）与AMR和MGEs之间存在关联。插入序列（如ISSsu5和ISLgar5）邻近AMR基因（如ANT(6)-Ia和tet(M)），与先前研究一致。泰国S. uberis未形成单一独特谱系，而是分散在全球多个谱系中，支持该物种缺乏地理结构。在核心基因组稳定的种群中，MGEs是微进化和短期适应主要驱动因素。此外，IS和前噬菌体在细菌物种和属间的遗传交换证据，特别是与猪链球菌（S. suis）来源的ISSsu5和乳酸乳球菌/粪肠球菌来源的前噬菌体同源性，表明微生物生态重叠促进跨属基因转移，可能通过中间宿主（如牛胃肠道）实现。研究局限性包括仅三个农场样本、不同采样年份、数据库局限性以及未检测染色体突变。这些发现强调了MGEs介导的基因转移在S. uberis进化中的重要性。
**研究结论（翻译自原文Conclusions）**
本研究揭示了从不同农场环境分离的S. uberis中与AMR基因相关的重组事件和MGEs。来自AMR基因流行率较高农场的分离株还表现出频繁检测到与这些基因共定位的MGEs，其中IS与AMR决定因子密切相关。尽管在农场A和C的S. uberis种群中观察到低重组率和保守的核心基因组，但本研究发现了高水平的MGE多样性。数据表明，细菌属内和属间存在持续的局部基因交换。特别是，S. uberis的IS和前噬菌体与奶牛场环境中常见其他细菌属的同源元件之间存在同源性，提示微生物生态学可能通过水平基因转移影响S. uberis的进化。总体而言，不同细菌种群中噬菌体多样性和MGE介导的基因传播变异可能在S. uberis的进化和短期适应中发挥重要作用。
