《Applied and Environmental Microbiology》：Fine-scale structural diversification among closely related Legionella pneumophila isolates from a bath-associated environmental reservoir
编辑推荐：
环境储库中的嗜肺军团菌（Legionella pneumophila）群体可表现基因组多样性，使暴发调查中的来源归因复杂化。基于核心基因组单核苷酸变异（core genome single-nucleotide variant, cgSNV）的方法被广泛用于分
环境储库中的嗜肺军团菌（Legionella pneumophila）群体可表现基因组多样性，使暴发调查中的来源归因复杂化。基于核心基因组单核苷酸变异（core genome single-nucleotide variant, cgSNV）的方法被广泛用于分子流行病学调查中鉴定紧密相关的谱系；然而，这些方法主要关注核苷酸水平的变异，可能忽视紧密相关分离株之间的结构变异。在本研究中，研究人员调查了一次军团菌病病例的流行病学调查期间获得的临床和环境嗜肺军团菌分离株的基因组多样性。通过分析来自单个样本的多个菌落，研究人员将cgSNV分析与比较基因组结构分析相结合，以表征高度相关群体内插入序列（insertion sequence, IS）介导的结构变异。共分析了属于序列型138（sequence type 138）的138株血清型1组分离株。cgSNV分析显示，临床分离株形成一个克隆谱系，而来自同一浴池系统的环境分离株尽管差异不超过20个cgSNV，却表现出显著多样性。对48株代表性分离株的完整基因组分析表明，基因组同线性几乎一致，大规模重排有限。相比之下，IS插入模式表现出相当大的变异性，特别是IS256，其在各分离株中的插入位点显示出广泛的多样性。即使在共享相同cgSNV基因型的分离株之间也观察到IS插入模式的差异，表明IS介导的结构变异可以独立于核心基因组SNV变异发生。相对于参考基因组，超过一半的编码序列差异与IS插入位点相关，且主要为IS256相关事件。这些发现表明，环境储库中紧密相关的嗜肺军团菌群体经历与动态IS活性相关的精细尺度多样化。
嗜肺军团菌（Legionella pneumophila）是一种广泛分布于天然及人工水系统中的机会性病原菌，通过吸入污染气溶胶导致军团菌病，公共浴池已被确认为日本军团菌感染的重要来源。全基因组测序技术的应用揭示了环境储库中嗜肺军团菌的遗传多样性，且暴发源内往往共存着多个紧密相关的基因型。在分子流行病学调查中，核心基因组单核苷酸变异（cgSNV）分析是鉴定紧密相关谱系的常用工具，但该方法侧重于核苷酸水平变异，难以识别由插入序列（IS）等移动遗传元件引发的结构变异。同时，以往研究中从单个样本分离的菌落数量有限，可能低估了源内基因组多样性。因此，有必要基于完整基因组对紧密相关分离株的结构变异进行系统解析。
本研究基于一起日本军团菌病病例的流行病学调查，研究人员从患者痰液和其使用过的浴池环境样本中分离了大量菌落，共获得138株血清型1组（SG1）且均为序列型138（ST138）的嗜肺军团菌分离株。通过整合cgSNV分析与完整基因组比较，以及蜡螟（Galleria mellonella）感染模型，研究人员刻画了该紧密相关群体内的结构变异和微进化特征。该成果发表于《Applied and Environmental Microbiology》。
主要关键技术方法包括：对所有分离株提取DNA后进行短读长测序（NovaSeq平台）和长读长测序（MinION平台），通过混合组装获得完整基因组；使用SNPcaster流程（整合BactSNP、Snippy和Gubbins）鉴定核心基因组SNV并构建网络；采用ISEScan预测插入序列，并基于侧翼基因定义锚位点以比较IS插入位点；通过BLAST比较编码序列（CDS）差异；代表性临床和环境菌株用于蜡螟幼虫感染试验，评估健康指数。样本队列来自2024年日本某浴池相关军团菌病病例调查，菌株分离自患者痰液、浴池水和浴缸表面拭子。
研究结果按照原文小标题总结如下：
**隔离概况**：研究人员从患者痰液、浴池水和浴缸表面拭子中分别获得60、57、22和1株分离株，共140株。其中138株为SG1 ST138，另2株为SG5 ST1631，后者在后续分析中被排除。该结果明确了菌株收集的来源和构成。
**cgSNV分析**：以临床分离株KLH309为参考基因组，138株SG1菌株的核心基因组大小为3,188,069 bp，共鉴定出136个cgSNV，最大两两差异为20个SNV。网络分析显示58株临床分离株与1株浴池水分离株形成同一节点，而大多数环境分离株分散在52个不同节点，表明环境菌株存在显著的精细遗传多样性。主要临床谱系包含62株，仅存在两个共享的非同义SNV。
**完整基因组测序**：对48株代表性菌株进行杂交组装，获得完整环化基因组。基因组大小介于3,401,904至3,423,096 bp之间，预测CDS数为3,015–3,035，rRNA和tRNA基因数量在所有菌株中一致。
**基因组结构变异**：与参考基因组KLH309相比，所有菌株保持了几乎一致的基因同线性，未发现大片段倒位、插入或缺失等大规模重排。BLAST比较显示CDS差异最多为19个/株，其中96.8%为片段化CDS，即序列相似性高但覆盖率降低。CDS差异数随cgSNV距离增加而增加（Spearman ρ=0.633，P<0.001）。将CDS差异并入cgSNV网络后总体拓扑结构不变，但部分共享相同cgSNV节点的菌株被区分开。
**IS插入模式多样性**：共鉴定出12个IS家族。IS256存在于所有菌株，插入位点广泛分布于137个参考锚位点，每个菌株携带11–20个拷贝，平均每个位点共享菌株数仅为4.8，显示出高度多样性。其他IS家族（如IS98、IS4、IS110等）分布相对固定，平均共享菌株数为15.1–48.0。相同cgSNV基因型的菌株间也观察到IS插入模式差异，表明IS介导的结构变异独立于核心基因组SNV积累。
**IS与CDS片段化的关联**：在342个CDS差异位点中，198个（57.9%）与IS插入位点重叠或位于10 bp范围内，其中IS256相关事件占96.5%。受影响的基因功能多样，无明显功能富集，提示IS介导的破坏在基因组中广泛发生。
**蜡螟感染试验**：选取6株代表性菌株（2株临床、2株与临床相关的环境、2株环境不同基因型）进行G. mellonella幼虫感染实验。所有菌株均引起幼虫健康指数随时间下降，但各时间点健康指数无显著差异（Kruskal–Wallis检验，P>0.05），对照组无变化。
讨论部分指出，临床分离株的高度克隆性提示单一谱系感染，而环境群体中多个基因型共存。尽管IS介导的结构变异广泛存在，但主要临床谱系中IS分布相对保守，且未观察到毒力差异，支持感染源自环境多样性群体的机会性暴露。本研究强调从单个样本中分析多个菌落对于准确表征源内多样性的重要性，并指出IS谱系分析为解析紧密相关菌株间关系提供了新的分辨率。局限性包括单次时间点采样、仅基于基因组推断等。
研究结论翻译为：总体而言，本研究结果强调了整合群体水平采样和结构基因组分析在嗜肺军团菌感染相关环境储库中解析基因组多样性的重要性。
