葡糖酸杆菌属(Gluconobacter)中16S rRNA相似性与基于基因组的物种界限的分歧(Discordance Between 16S rRNA Similarity and Genome-Based Species Boundaries in Gluconobacter)

《MicrobiologyOpen》：Discordance Between 16S rRNA Similarity and Genome-Based Species Boundaries in Gluconobacter

　　

摘要：16S rRNA基因被广泛用于细菌鉴定，但在若干类群中其分辨能力有限。在葡糖酸杆菌属(Gluconobacter)中，尽管依赖基于16S的分类，物种界限仍不明确。为考察这些局限性，研究人员分析了所有公开可用的完整基因组(n=12)及经整理的16S rRNA序列，并将16S rRNA序列相似度与基于基因组的亲缘关系进行比较。在所取样的完整基因组内，16S rRNA序列高度保守，种内与种间相似度存在大量重叠。相反，基于平均核苷酸一致性(average nucleotide identity, ANI)的离散度显著：ANI值在79.7%至100%之间，且氧化葡糖酸杆菌(G. oxydans)包含多于一个基因组谱系。值得注意的是，两株菌株(H24和1.637)与泰国葡糖酸杆菌(G. thailandicus)的ANI值(97.2%–97.4%)高于与其他G. oxydans菌株的ANI值。即便是超可变区V4、V6和V8也缺乏分支特异性特征，表明该组内物种水平分辨力有限。此外，G. oxydans、弗拉特葡糖酸杆菌(G. frateurii)和日本葡糖酸杆菌(G. japonicus)的模式菌株16S rRNA序列相似度≥99%，并在SILVA NR99数据库中归入同一非冗余代表序列，同样表明物种水平分辨力不足。这些结果显示葡糖酸杆菌属中16S rRNA序列相似度与基因组结构之间存在明显不匹配。准确的分类学及基于微生物组的分析离不开以基因组标准为依据、辅以精细化参考数据库和验证基因标记的方法。

论文解读：葡糖酸杆菌属(Gluconobacter)中16S rRNA相似性与基于基因组的物种界限的分歧

该研究发表于《MicrobiologyOpen》。

研究背景与立项依据

16S rRNA基因因具普遍保守性及含可变与保守区的组合，是细菌分类学与微生物组研究中最常用的分子标记，通常可有效实现属(level)水平鉴定。然而其分辨能力在各分类单元间差异显著，近缘种常共享几乎完全相同的序列，阻碍物种(species)水平区分。这种局限在临床与环境研究中已有详载，仅凭16S易致过度归类或误鉴定。此外许多细菌基因组含多个16S rRNA操纵子(operon)， intragenomic变异虽细微却可致ASV(amplified sequence variant)坍缩或夸大感知多样性。葡糖酸杆菌属(Gluconobacter, Acetobacteraceae)是好氧醋酸菌，广泛用于食品发酵与氧化生物工艺，具有重要工业价值，但该属分类数十年来受困于不同物种间16S rRNA序列极度保守。微生物组数据集中reads常被归为G. oxydans或Gluconobacter sp.，而G. frateurii、G. cerinus、G. japonicus等公认种无法仅凭16S可靠分离。SILVA、RDP等参考数据库中多模式菌株序列共享>99%相似度，常被归入同一非冗余参考簇，掩盖物种界限。虽有研究用recA、rpoB、adhA或多位点序列分析(multilocus sequence analysis, MLSA)获改善之分辨力，且16S–23S ITS区分析提示某些种含需重新考量的多重谱系，但关联intragenomic异质性与基因组水平离散度的系统性分析仍缺。因此16S在多大程度上反映物种界限、当前参考数据库可靠性如何等根本问题尚无答案。本研究定量评估现有Gluconobacter完整基因组内16S rRNA基因保守程度与多样性，整合16S rRNA系统发育分析与全基因组平均核苷酸一致性(average nucleotide identity, ANI)，评价基于16S分类的局限及其对分类学与微生物组研究的影响。

主要关键技术方法概述

研究人员选取截至2025年10月NCBI中所有Gluconobacter属完整基因组(n=12，代表9个种)，用DFAST重注释并提取16S rRNA特征；补充LPSN所列18个有效发表种模式菌株16S rRNA序列及SILVA SSU Ref NR99 release 138.2中长度≥1350 bp的Gluconobacter标注序列。用FastANI(v1.33)计算12个完整基因组两两ANI值。用MAFFT(--auto)比对所有16S rRNA序列，snp-dists算SNP，全局序列相似度按多序列比对中排除gap后相同核苷酸占比计算。以Escherichia coli K-12 MG1655 16S rRNA基因为可变区(V1–V9)位置参照。取75条独特16S rRNA序列（含Acetobacter aceti NCIB 8621^T作外类群）用IQ-TREE v2.4.0（GTR+G模型，ModelFinder选模，100次bootstrap）建最大似然(maximum-likelihood)系统发育树，iTOL可视化。

研究结果

3.1 Genomic Overview and 16S rRNA Dataset Characteristics（基因组概览与16S rRNA数据集特征）

分析涵盖12个完整基因组（9个种）、18个模式种16S rRNA序列及SILVA中111条序列。完整基因组来源之16S rRNA基因为每基因组4–5拷贝（共55拷贝）。全株内各拷贝至多33个SNP，单个基因组内16S rRNA拷贝间仅0–2个SNP，intragonomic异质性极低；G. oxydans种内16S rRNA序列离散度最高（达27–29个SNP）。FastANI显示12株间ANI值79.7%–100%，以95% ANI为物种界限 cutoff，G. oxydans至少分两个不同基因组簇：DSM 3504、ATCC 9937、621H相互ANI 96.0%–100%为一簇；H24与1.637相互ANI 99.9%但与上述簇约~80%，反而与G. thailandicus HD924有更高ANI值(97.2%–97.4%)，超过常见物种阈值，提示H24和1.637可能更接近G. thailandicus而非严格意义上的G. oxydans。被分析完整基因组分两大组：第一组含G. sphaericus、G. oxydans、G. roseus、G. albidus（组内ANI 86.2%–100%）；第二组含G. thailandicus、G. japonicus、G. cerinus、G. frateurii（组内ANI 83.3%–100%），对应既往报道"Gluconobacter clade I"。Clade I与G. albidus具5拷贝16S rRNA/基因组，其余种为4拷贝。

3.2 Sequence Identity Analysis of the 16S rRNA Gene（16S rRNA基因序列相似度分析）

模式菌株种间16S rRNA相似度为97.36%–100%，大量比较超过常用98.7%物种阈值——G. oxydans与G. roseus甚至100%相同，多数种间比较>99%。扩大至含参考序列的数据集，同种内相似度96.3%–100%，不同种间大量配对≥98.7%（如G. oxydans与G. roseus、G. albidus、G. sphaericus等多达≥99%），而G. cerinus、G. frateurii、G. japonicus部分种内配对反而低于98.7%。G. oxydans与除G. kondonii和G. morbifer外几乎所有种交叉超阈值，呼应ANI发现之G. thailandicus相关株。整体呈"高度保守却分类不一致"悖论：16S rRNA适合属水平指派，但无足够分辨力做Gluconobacter内物种水平分类。

3.3 Phylogenetic Relationships and Variable-Region Characteristics（系统发育关系与可变区特征）

基于74条非冗余序列的最大似然树重建出3个主分支(clades 1–3)，但种名与系统分支明显不符——多模式菌株未形成单系种特异性分支，G. oxydans、G. frateurii、G. japonicus序列分散于数枝；部分分支跨越多种边界（部分G. oxydans近G. thailandicus或G. frateurii，个别G. thailandicus嵌于G. oxydans分支），种内变异与种间分化大量重叠。可变区分析显示大多数碱基替换集中于V4(585–633 bp)、V6(984–1035 bp)、V8(1251–1292 bp)区，其余区高度保守；即便此三区也无清晰分支特异性特征，组间差仅数个核苷酸，证实16S rRNA基因无Gluconobacter物种水平系统发育分辨力。

讨论与结论总结

本研究目的非提出新分类修订，而是在核糖体操纵子完全解析条件下定量评估16S rRNA相似度与全基因组离散度的对应关系。与既往观察一致，Gluconobacter 16S rRNA序列分化极有限，种内种间相似度大幅重叠，反映强功能约束与慢进化速率，故适于属级鉴定却无法划定近缘种。ANI分析揭示G. oxydans含多个相异基因组谱系，部分株更近G. thailandicus；既往已提示此种并系性(polyphyly)，本研究直接将其与16S rRNA近乎一致相关联。即便超可变区(V4、V6、V8)差异微小且无谱系特征，强调基于16S rRNA物种分类的内在局限——凭16S做物种水平指派有相当大误分类风险，可能延续过时命名、掩盖真实基因组关系。此局限不独见于Gluconobacter而广布多类细菌。明确鉴定应采基因组指标如ANI或digital DNA–DNA hybridization(dDDH)，16S rRNA基因序列更适合初步宽泛筛查。对微生物组研究而言，16S rRNA调查所报"G. oxydans"可能是无法物种水平解析的多物种复合体混合物，故属水平报告（如Gluconobacter sp.）更可靠保守。当前Gluconobacter分类未充分反映基因组结构；G. oxydans或需分类复评，与邻种界限待进一步评估。未来工作应纳入基因组物种界限、开发验证单基因标记（如rpoB、groEL），完善SILVA与GTDB等数据库注释，增扩高质量基因组组装以延伸观察至全属多样性。总之，虽16S rRNA基因在细菌系统学具历史核心地位，在Gluconobacter中因其极端保守无法可靠推断物种界限——即便基因组水平分化显著。本研究通过对现有完整基因组中16S rRNA多样性与基因组离散度的定量比较，阐明为何16S rRNA序列不支持该群组物种水平分类。精确的Gluconobacter分类学需以基因组框架为基础，辅以精选参考数据库与验证单基因标记，并对物种界限重新审视。