克洛彭施泰特棒杆菌复合群的比较泛基因组学揭示物种结构与附属基因组多样性

《Frontiers in Microbiology》：Comparative pangenomics of the Corynebacterium kroppenstedtii complex reveals species structure and accessory genome diversity

　　

背景：克洛彭施泰特棒杆菌复合群（CKC）在乳腺炎中检出率上升，但其物种边界与基因组多样性尚不清楚。方法：经质量筛选后，研究人员分析了56个公开基因组，涵盖C. kroppenstedtii（Ck）、C. parakroppenstedtii（Cpak）与C. pseudokroppenstedtii（Cpse）。分类学指派由平均核苷酸一致性（ANI）、数字DNA–DNA杂交（dDDH）、基于平均氨基酸一致性（AAI）的蛋白相似性及核心基因组系统发育支持。随后考察泛基因组结构、移动元件、抗性基因、功能富集与保守基因邻域。结果：ANI与核心基因系统发育将Ck、Cpak、Cpse解析为三个物种。4706个聚类的泛基因组为开放性（拟合指数=0.232），含1751个核心聚类与2955个附属聚类。附属基因谱进一步解析出Cpak-P1至Cpak-P3三个谱系，与核心基因组系统发育一致。插入序列计数与附属聚类计数相关（Spearman ρ=0.640，p=1.12×10?7），附属聚类比核心聚类更频繁重叠预测移动遗传元件（MGEs）（402/2955 vs 98/1751；优势比=2.66，p=3.31×10?19），支持MGE–附属关联但不暗示因果。检出7类抗菌素抗性基因（ARGs），最常见为ermX、cmx、APH(6)-Id与tet(W)。AC410样序列见于35/56基因组，且与ARG携带强相关（34/35 vs 4/21；优势比=120.5，p=9.93×10?10）。8个完整基因组中有5个AC410样序列位于染色体而非独立质粒contig。功能富集与保守邻域鉴定出Cpak-P1与Cpak-P2共享的候选脂质–铁位点。结论：结果为CKC提供谱系解析的基因组框架，显示具实质物种与谱系相关附属多样性的开放泛基因组；进一步识别AC410样抗性相关基因组背景与候选脂质–铁位点供未来功能研究。

研究背景与立项依据

克洛彭施泰特棒杆菌复合群（Corynebacterium kroppenstedtii complex，CKC）自1998年首次描述以来，主要从非泌乳性乳腺炎、肉芽肿性炎症及复发性乳腺脓肿标本中恢复，但因嗜脂生长需求与常规鉴定方法对物种分辨率不足，临床与基因组研究中长期被归并为“C. kroppenstedtii”旧称。后续基因组与表型分析确立了C. parakroppenstedtii（Cpak）与C. pseudokroppenstedtii（Cpse）两个新种，但历史数据库标签混乱，导致既往关于耐药决定子、致病机制的结论难以归属到正确物种。已有研究提示CKC成员存在脂肪酸合成基因丢失、铁螯合糖脂分泌等谱系特异性性状，但缺乏在正确物种边界下对附属基因组、移动遗传元件（MGEs）、抗性组与代谢位点的系统性比较。该研究即在此背景下，对公开基因组做质量再筛选与物种重指派，构建谱系解析的泛基因组框架，论文发表于《Frontiers in Microbiology》。

主要技术方法概览

研究人员从GenBank检索标记为CKC的组装基因组及元数据，以Ck DSM 44385、Cpak MC-26、Cpse P15_C1为参考，用fastANI（≥95% ANI）与QUAST/MultiQC质控筛得56个高质量基因组（2 Ck、47 Cpak、7 Cpse，含8个完整基因组，来源以中国样本为主）。注释采用Bakta，物种边界辅以dDDH（GGDC 3.0）与AAI（FastAAI）；泛基因组用PGAP2严格模式（MMseqs2聚类）划分严格核心/软核心/壳/云四类；核心基因组系统发育用MAFFT比对+IQL-TREE最大似然树，并借Gubbins滤重组后重算cgSNP；附属基因组用Jaccard相异度PCoA与PERMANOVA/PERMDISP、PAM分群；MGEs用geNomad、ISEScan、MOB-suite预测；ARG用CARD（80/80阈值）定位；AC410样元件综合MOB-suite分类与完整基因组BLASTn共证；功能注释整合eggNOG-mapper、SignalP、TMHMM，富集用Fisher精确检验+BH校正。

研究结果

基因组重指派解析复合群内三物种

经ANI重指派，56个基因组归为Ck、Cpak、Cpse三类；与非CKC棒杆菌最高ANI仅80.87%，而群内≥98.16%；dDDH种内≥70%、种间45.2%–51.9%；AAI种内Jaccard 0.9869–0.9996、种间0.9168–0.9246；66分类单元核心基因组系统发育单系支持三物种框架，Cpak与Cpse为姐妹群。

CKC具以附属聚类为主的开放泛基因组

PGAP2得4706聚类（1496严格核心、255软核心、598壳、2357云），附属占62.79%；幂律拟合y=1788.064×x^0.232证实开放泛基因组。按核心系统发育排布时，壳/云存在缺失谱沿物种分支变化。QC限制子集（38基因组）仍开放（指数0.149），但云比例由50.08%降至34.14%，说明云聚类丰度受组装质量影响。

附属基因内容区分物种并解析Cpak三谱系

1596非单例附属聚类PCoA中物种解释54.2% Jaccard变异（PERMANOVA R²=0.542，p=0.001），PERMDISP不显著；Cpak内PAM最优三群Cpak-P1（36）、Cpak-P2（7）、Cpak-P3（4），与重组过滤后核心树一致，cgSNP中位数P1–P2间579、P3距P1/P2超8000。谱系富集分别获37/107/168个关联聚类（其中1/79/137为谱系限定），100次分层交叉验证宏F1=1.000（P3因样本少需谨慎）。

MGEs呈谱系结构并与附属变异关联

Cpak-P1插入序列（IS）负担中位最高（0.53%序列覆盖），组间IS负担/geNomad质粒样/前病毒Kruskal–Wallis显著。全56基因组IS计数与附属聚类计数Spearman ρ=0.640（p=1.12×10^?7），校正基因组大小后偏ρ=0.397（p=0.0027），QC子集与Cpak内校正后仍显著。MGE重叠CDS实例中附属CDS 8.31% vs 核心0.55%，聚类水平附属vs核心重叠优势比2.66（p=3.31×10^?19），CDS实例级逻辑回归OR=15.88，仅支持关联非因果。

AC410样区域与ARG携带强相关

AC410样证据（MOB-suite+完整基因组BLASTn）阳性35/56，集中于Cpak-P1（27/36）与Cpse（5/7），Cpak-P2 3/7，P3与Ck全阴。阳性组34/35携至少1个CARD定义ARG（ermX 35、cmx/APH(6)-Id 33、tet(W) 31等），阴性组仅4/21，OR=120.5（p=9.93×10^?10）。8完整基因组中5个AC410样定位于染色体；以Cpak SAMN57172177的35.8 kb染色体区段与C. resistens质粒pJA144188及他种同源染色体区做BLASTn共线性，确认sul1/cmx/aph(6)-Id/aph(3″)-Ib/aph(3′)-Ia等ARG与AC410参考块共定位，但无复制子/oriT/交配对形成系统，不能证其自主转移。

功能分化指向候选脂质–铁位点

Cpak-P2富集转运（OR=3.77，FDR=2.52×10^?5）、辅酶转运代谢（COG H）、无机离子转运（COG P）；Cpak-P3富集跨膜蛋白（OR=2.17）。预设15个脂质/铁相关聚类中，酰基转移酶–MMPL转运蛋白–短链脱氢酶/还原酶片段在Cpak-P1（36/36）与Cpak-P2（7/7）全保留，Cpak-P3/Ck/Cpse缺如；邻近严格核心锚定区间在装配连续代表中保守，命名为候选脂质–铁位点，生化产物未明。

讨论与结论翻译

讨论指出，旧“Ck”标签混淆了生物学各异的物种与谱系，未来临床研究须先做基因组分辨鉴定；开放泛基因组指数不指认机制，核心–附属一致性表明种群史贡献大；IS–附属关联与CDS级MGE重叠支持MGE相关成分但不证因果；AC410样是抗性背景而非已证可转移质粒；ARG检出不能替代药敏与表达数据；候选脂质–铁位点联系Cpak糖脂致肉芽肿乳腺炎假说但未证同分子。局限含地理采样偏倚（中国为主）、Ck/P3样本少、组装碎片影响云基因与邻域、横断面数据不能溯获取史。

结论（译）：所分析CKC基因组形成三个物种级群组，Cpak含三个具 distinct 附属谱的工作谱系；泛基因组开放，MGE相关测度与附属变异关联但未确立MGE为开放主驱；AC410样证据呈谱系偏倚且与ARG携带强相关；候选脂质–铁位点区分Cpak-P1/P2与P3/Ck/Cpse。这些模式为未来物种意识监测与机制研究提供框架，但临床、移动性与致病性主张需表型匹配与实验验证。