《Applied and Environmental Microbiology》：Genomic divergence and phenotypic heterogeneity of Legionella longbeachae reveal a putative novel proinflammatory serogroup with genome reduction and mobilome expansion
编辑推荐：
长滩军团菌（*Legionella longbeachae*）是一种主要存在于土壤环境中的新兴呼吸道病原体，但其基因组结构和进化策略仍然缺乏系统表征。本研究通过对242株长滩军团菌分离株进行泛基因组学和功能分析，研究人员揭示了一个由广泛水平基因转移驱动的开放
长滩军团菌（*Legionella longbeachae*）是一种主要存在于土壤环境中的新兴呼吸道病原体，但其基因组结构和进化策略仍然缺乏系统表征。本研究通过对242株长滩军团菌分离株进行泛基因组学和功能分析，研究人员揭示了一个由广泛水平基因转移驱动的开放泛基因组，以及核心基因组和辅助基因组之间显著的功能分化。研究人员鉴定出一个系统发育上独特的谱系，命名为推定血清群3（putative serogroup 3, sg3），该谱系来源于中国环境样本。该谱系通过基因组缩减（平均大小为4.01 Mb）、质粒缺失以及缺少关键N-乙酰转移酶编码基因（*orf9*）的截短O-抗原生物合成簇，与经典血清群1和2发生分化。在多种人类细胞系的体外感染模型中，与sg1和sg2相比，sg3诱发更强的宿主细胞死亡和上皮细胞中促炎细胞因子转录水平的升高。经系统发育校正的关联分析显示，该表型与单个辅助毒力因子之间无显著关联。相反，秩转换的系统发育广义最小二乘（phylogenetically generalized least squares, PGLS）回归表明，该表型与移动基因组（COG X）扩张呈显著正相关，而细胞内转运（COG U）和防御（COG V）基因库则呈收缩状态。这些发现揭示了sg3与经典血清群不同的独特致病特征，其特点为炎症激活增强而非免疫逃逸。本研究深化了对长滩军团菌种群结构的理解，对传统血清分型的还原论应用提出了挑战，并强调了对表现出增强细胞毒性及促炎反应（可能与其移动基因组扩张相关）的谱系进行修订诊断考量和持续监测的必要性。
## 论文解读
### 一、研究背景与科学问题
嗜肺军团病（Legionnaires' disease, LD）是由军团菌属细菌引起的一种严重的非典型肺炎，长期以来构成全球公共卫生挑战。尽管嗜肺军团菌（*Legionella pneumophila*）是绝大多数已记录暴发事件的病原体，长滩军团菌（*Legionella longbeachae*）却作为一种临床和生态学上均独特的病原体逐渐浮现，尤其在澳大利亚和新西兰呈地方性流行，东南亚和欧洲的报道也日益增多。与嗜肺军团菌的水生生态位适应不同，长滩军团菌主要在土壤环境（尤其是堆肥、盆栽土和园艺覆盖物）中生存，这为园艺工人、景观设计师和免疫功能低下人群创造了独特的暴露途径。在临床表现上，长滩军团菌感染呈现典型的军团病症状，包括高热、肺部实变、神经系统受累和胃肠道不适。然而，积累的临床和免疫学证据表明，长滩军团菌可能诱发与嗜肺军团菌不同的宿主免疫应答、表现出差异化的组织嗜性和可变疾病严重程度，但这些差异的分子和基因组基础仍不清楚。
历史上，长滩军团菌的分类学和流行病学追踪依赖于基于巨噬细胞感染增强子（macrophage infectivity potentiator, *mip*）基因序列变异的分子血清分型。目前该物种被划分为两个公认的血清群（sg1和sg2），其中sg1占全球人类感染的绝大多数。然而，传统血清分型方法可能无法充分反映环境军团菌群体固有的广泛基因组可塑性和功能异质性。近年来比较基因组学研究已经确定军团菌属物种拥有高度开放的泛基因组，其由频繁的水平基因转移、移动遗传元件增殖和同源重组塑造。但公认的血清群框架是否能真正捕捉长滩军团菌多样性的全部范围仍然不确定，特别是当新兴环境分离株日益超出既有sg1/sg2边界时。谱系特异性分化在宿主细胞入侵、胞内复制和免疫调节能力方面的驱动因素仍不清楚。这种分化可能源于离散毒力效应因子的获得或丢失，也可能源于更高阶的基因组结构变化，包括移动基因组扩张、O-抗原簇重塑和适应性基因组缩减。解决这一问题需要整合比较基因组学、功能注释和宿主-病原体分析，并置于系统发育严格的框架中。
### 二、研究内容与主要结论
研究人员对242株长滩军团菌基因组开展了泛基因组学和功能学系统研究，将39株新测序的中国环境分离株与203个高质量公共基因组组装整合分析。研究 delineated 种群结构、O-抗原簇多样性和核心与辅助基因组的功能组成，鉴定出一个源自中国环境的独特谱系，该谱系表现出显著的基因组缩减和移动基因组富集。通过整合O-抗原簇结构分析、*mip*等位基因谱分析和遗传相关性分析，研究人员将该谱系与经典sg1和sg2分离株区分开来，将其命名为推定新型血清群（putative serogroup 3, sg3）。通过多种人类细胞类型的系统性体外感染实验，研究人员证明了在内化效率、宿主细胞毒性和促炎细胞因子诱导方面的显著表型异质性。通过将系统发育广义最小二乘（phylogenetically generalized least squares, PGLS）回归与功能分析相结合，研究发现表型分化与单个辅助毒力基因的差异分布无强烈关联，而是与协调的基因组重组和移动基因组扩张显著相关。
### 三、关键技术与方法
本研究采用了多项关键技术方法。泛基因组构建使用Bakta v1.11.4进行基因组注释，Panaroo v1.5.1（≥98%同一性和≥95%比对覆盖率）进行泛基因组推断，采用Heaps定律拟合评估泛基因组开放度。平均核苷酸一致性（average nucleotide identity, ANI）通过fastANI v1.34计算，物种界定通过in silico DNA-DNA杂交（iDDH）验证。质粒检测使用geNomad 1.8.0，毒力因子谱分析通过ABRicate比对毒力因子数据库（VFDB）。系统发育推断采用IQ-TREE v3.0.1最大似然法，共线性分析使用Mauve 2.3.1。功能表型评估采用多细胞系体外感染模型——包括PMA分化的THP-1巨噬细胞样细胞、A549肺泡上皮细胞和HEK293T胚胎肾上皮细胞，检测内化效率（流式荧光淬灭法）、胞内复制动力学（菌落计数）、细胞毒性（annexin V/PI双染流式细胞术）和细胞因子转录（qRT-PCR）。样本队列来源为39株中国环境分离株（来自上海和滁州等地土壤、盆栽混合基质和水源）。
### 四、研究结果要点
**1. 泛基因组分析揭示长滩军团菌是一个开放、功能异质性的物种**
泛基因组被划分为核心、软核心、壳层和云层基因集。云层基因存在于不足15%的分离株中，占总遗传库的40%以上，表明显著的基因组多态性。泛基因组动力学分析显示核心基因组快速饱和，而泛基因组随着基因组数量的增加持续扩张，计算开放度参数γ=0.11，符合开放泛基因组模型。基因频率分布呈偏斜长尾分布，由稀有低拷贝基因主导。功能比较显示核心基因组主要富集代谢和翻译等持家功能，而辅助基因组显著富集复制/重组/修复（COG L）、移动基因组（COG X）和防御机制（COG V）相关基因。
**2. 基因组异质性揭示长滩军团菌内的一个独特谱系**
以参考菌株NCTC11477为基准的平均核苷酸一致性分析解析出三个遗传上不同的簇。242株分离株的成对ANI分析确认，231株具有指定血清群状态的分离株严格划分为sg1和sg2，而11株中国环境分离株形成一个离散的系统发育分歧簇。进一步对*mip*基因的检测发现这11株分离株携带与sg1或sg2均不一致的*mip*多态性。
**3. O-抗原簇分歧支持推定新型血清群的认定**
O-抗原生物合成位点的比较分析在长滩军团菌中鉴定出27个基因，包括*rfbF*、*rfbG*和*wecA*以及24个开放阅读框（orf1-orf24）。sg1位点跨23.8kb含20个ORF，sg2位点跨17.1kb含15个ORF，新簇跨16.5kb含15个ORF但独特地缺少编码推定N-乙酰转移酶的*orf9*。结合独特的*mip*谱、遗传相关性和O-抗原位点结构，11株分离株被指定为推定新型血清群（sg3）。所有11株分离株与模式株共享>96% ANI，iDDH值范围70.6%–70.9%，均远高于70%物种边界，明确支持这些分离株属于长滩军团菌物种而非独立物种。
**4. 推定sg3谱系的系统发育独特性和基因组缩减特征**
核心基因组SNP和泛基因组系统发育一致地将sg1、sg2和sg3解析为三个离散的充分支持的分支，sg3形成单系簇。共线性比较显示广泛共线性但sg3基因组明显收缩。定量分析表明sg3平均基因组大小（4.01 Mb）显著小于sg1（4.12 Mb）和sg2（4.20 Mb）（单因素方差分析，P<0.0001）。质粒筛查显示sg3分离株缺乏可检测质粒，而绝大多数sg1（208/209，99.52%）和大多数sg2（16/22，72.73%）分离株携带一个或多个质粒序列。基于ANI的聚类显示sg2质粒与部分sg1分离株所携带质粒聚类，与地理来源或分离来源无关，表明这两个血清群之间活跃的质粒交换；sg3则与之形成鲜明对比，不携带质粒，提示其可能在质粒介导的基因交换方面与sg1/sg2存在遗传隔离。
**5. 推定血清群3与sg1和sg2相比表现出增强的细胞毒性和促炎能力**
基于COG类别丰度的主坐标分析（PCoA）显示sg3占据与sg1和sg2不同的功能空间。感染表型分析显示上皮细胞内化率在三个血清群间相当，但THP-1细胞中sg2的内化效率显著高于sg1和sg3。值得注意的是，sg3感染在三种细胞系中均诱导显著的宿主细胞死亡，该细胞病变效应严格依赖活菌——热灭活sg3（70°C，30 min）不能降低细胞活力。热灭活细菌制剂仅诱导低水平、血清群间无差异的细胞因子应答，表明增强的炎症应答需要活跃的细菌过程（如Dot/Icm分泌系统）。在293T细胞中，sg3在大部分时间点一致地诱导更高水平的IL1α、IL1β、IL6、IL8和TNFα转录；在6h和12h时五种细胞因子均显著上调；在24h时IL1α、IL1β、IL6和IL8维持显著升高。在A549细胞中，sg3的显著细胞因子诱导局限于早期（6h）。IL10在各时间点维持在基线水平。这些表型差异并非由胞内增殖差异解释，细菌增长率在各血清群间相当，表明增强的细胞毒性和炎症应答具有细胞类型和血清群特异性。
**6. 血清群特异性宿主互作表型的功能基因组相关性**
推定sg3表现出移动基因组（COG X）显著扩张伴随细胞内转运/分泌/囊泡运输（COG U）和防御机制（COG V）基因库收缩的独特功能特征。通过对176个精选毒力因子（包括已知Dot/Icm效应子和表面结构）的筛查鉴定出148个核心和28个谱系特异性位点，但PGLS分析显示经FDR校正后任何单个毒力基因的存在/缺失与宿主互作表型无显著关联。应用秩转换PGLS到COG类别丰度后，发现COG X丰度与293T细胞中两个促炎细胞因子水平呈显著正相关：6h IL1β（系数=4.51，P=0.00017，q=0.013）和24h IL1α（系数=2.02，P=0.00028，q=0.013）。无监督层次聚类显示sg3分离株聚为一类，较高的COG X丰度与促炎谱相聚，较高的COG U/COG V与宿主细胞存活相聚。这表明sg3增强的炎症表型可能归因于移动基因组扩张这一高阶基因组重组特征，而非单个毒力因子的获得或丢失。
### 五、讨论与结论
开放泛基因组（γ=0.11）将长滩军团菌归入具有高度基因组流动性的细菌物种之列。云层基因的高比例和辅助基因组中移动基因组、复制/重组和防御功能的富集强调了水平基因转移在塑造多样性中的作用，这可能反映了对复杂土壤环境的适应。sg3的发现为理解长滩军团菌多样性增添了新维度：其独特的*mip*多态性、O-抗原簇结构、系统发育分离和表型分歧支持其作为推定新型血清群的临时分类。sg3缺少*orf9*（编码推定N-乙酰转移酶）这一事实，考虑到N-乙酰化修饰深刻影响脂多糖（lipopolysaccharide, LPS）免疫原性和血清学特异性，暗示传统sg1/sg2特异性检测可能对推定sg3灵敏度降低，这可能在临床诊断环境中产生潜在影响。sg3的基因组缩减（4.01 Mb）和质粒完全缺失定义了其独特的基因组结构和进化轨迹，但缩减的基因组并未削弱其宿主互作能力，反而与特定细胞类型中增强的细胞毒性和促炎能力共存。移动基因组扩张与促炎表型的显著关联提示高阶基因组重组而非离散效应子获得可能与该谱系的独特致病策略相关。本研究对基因组大小与毒力关系的传统观点提出质疑，强调了对这类表现增强细胞毒性和促炎应答的谱系进行持续基因组监测的必要性。未来研究应优先推进sg3 LPS的结构表征、sg3特异性效应子的系统筛选、原代细胞和动物模型研究、sg3特异性诊断试剂的开发以及定义其流行病学意义的临床研究。
