Anaeroglobus作为儿童龋病候选唾液生物标志物的多队列外部验证

《Frontiers in Cellular and Infection Microbiology》：Multicohort external validation of Anaeroglobus as a candidate salivary biomarker for childhood caries

　　

目的：评估Anaeroglobus作为儿童龋病候选唾液微生物组信号（salivary microbiome signal）的跨队列可重复性，并将该信号置于牙菌斑群落、功能及宿主反应背景下。方法：研究人员使用公共儿童唾液16S rRNA（16S ribosomal RNA）队列，采用发现–初级外部验证–第二验证框架。经中心对数比变换（centered log-ratio transformation, CLR）和协变量调整后评估属级特征，并以牙菌斑微生物组、牙菌斑转录组和深龋牙髓单细胞数据集（single-cell datasets）作为支持证据层。结果：在发现队列中，乳杆菌属（Lactobacillus）与龋病状态相关，厌氧球属（Anaeroglobus）与总龋负担相关。在初级外部队列中，Anaeroglobus和Lactobacillus仍与龋病呈正相关[比值比（odds ratio, OR）=3.452（95%置信区间（confidence interval, CI）：1.851–6.449）和OR=2.118（95% CI：1.049–4.282）]。在第二验证队列中，Anaeroglobus方向一致[OR=3.009（95% CI：1.241–7.312）]。初级外部微生物组模型显示出中等判别能力[曲线下面积（area under the curve, AUC）=0.782；近似95% CI：0.703–0.861]。探索性合并效应分析（pooled-effect synthesis）支持正向方向一致性，而支持性牙菌斑及多组学层表明更广泛的群落重构、碳水化合物/氧化应激适应及宿主免疫-基质重塑。结论：Anaeroglobus代表一个方向可重复的候选唾液风险相关信号，而非直接的机制性或即时诊断性生物标志物。36样本的第二验证队列应解读为支持性年龄外推证据，合并效应分析应视为探索性结果；需进一步开展前瞻性、配对唾液-牙菌斑及多模态验证。

**论文解读：Anaeroglobus作为儿童龋病候选唾液生物标志物的多队列外部验证**

**研究背景与问题**

儿童龋病，尤其是早期儿童龋病（early childhood caries, ECC），是全球范围内最普遍的口腔疾病之一，未经治疗的龋病负担依然沉重。世界卫生组织（World Health Organization, WHO）的实施指南强调以预防为导向的早期干预。病因学层面，对龋病的解释已从传统单一病原体观点转向以牙菌斑生物膜（dental plaque biofilms）、群落失衡（community imbalance）和生态压力（ecological pressure）为核心的框架。生态学观点进一步被后续概念性工作所完善，维持口腔微生物恢复力（oral microbial resilience）被认为对疾病预防至关重要。随着高通量测序和口腔多组学（oral multi-omics）的快速发展，越来越多的研究致力于识别候选微生物属、构建风险预测模型，并评估微生物组相关信号是否能在队列间复制。然而，儿童口腔微生物组本身具有明显的年龄依赖性发育特征和显著的队列间异质性，生态演替持续受出生后因素和发育阶段影响。因此，当前ECC研究日益强调多物种生物膜（multispecies biofilms）和基于生态位（niche）的解释，而非单一属的绝对决定论。基于此，本研究遵循三层唾液分析路径（发现、初级外部验证、第二独立外部验证），评估Anaeroglobus及相关候选属的方向一致性和外部可重复性，同时引入同一生态位牙菌斑证据、牙菌斑转录组和深龋牙髓单细胞发现作为支持层，以提供群落背景、功能背景和宿主反应背景。

**研究内容与结论**

研究人员开展了多队列外部验证研究，利用公共儿童唾液16S rRNA队列，采用发现–初级外部验证–第二验证框架。在发现队列中，Lactobacillus与龋病状态独立相关，Anaeroglobus与总龋负担呈正相关，但双胞胎对内分析中，经错误发现率（false discovery rate, FDR）校正后无属级特征显著，提示严格家族约束下单个属效应有限。在初级外部验证队列（130样本）中，Anaeroglobus和Lactobacillus均保持方向一致且统计显著，Anaeroglobus的调整后OR为3.452（95% CI：1.851–6.449），Lactobacillus的OR为2.118（95% CI：1.049–4.282）。在第二独立外部验证队列（36样本，5岁儿童）中，Anaeroglobus仍方向一致且显著[OR=3.009（95% CI：1.241–7.312）]，而Lactobacillus方向一致但未达到统计显著性。初级外部微生物组模型的AUC为0.782（近似95% CI：0.703–0.861），表明中等判别能力。区域亚组分析和基于ComBat的探索性批次效应敏感性分析显示，Anaeroglobus和Lactobacillus的效应方向未受实质性改变，ComBat调整后AUC为0.746（95% CI：0.662–0.831）。同一生态位牙菌斑公共数据库支持显示，跨队列牙菌斑研究发现群落级差异和模型判别能力较单属更稳定，更稳定的属为Alloprevotella和Megasphaera，而非Anaeroglobus。牙菌斑转录组功能支持表明，碳水化合物利用转录本约占总转录本的10%，单糖/二糖代谢相关转录本约占15%，氧化应激相关转录本占应激反应转录本的50%–75%，揭示了龋病相关生物膜的碳水化合物代谢、酸适应和氧化应激适应背景。深龋牙髓单细胞数据支持显示，深龋（而非牙釉质龋）改变了牙髓免疫细胞组成，并上调了促炎、抗炎和矿化相关基因，提示宿主免疫-基质重塑。探索性合并效应分析显示，Anaeroglobus的合并随机效应OR为3.335（95% CI：2.087–5.329，I2=0.0%），Lactobacillus的合并OR为2.082（95% CI：1.310–3.308，I2=0.0%），留一法敏感性分析支持信号稳定性。该研究发表在《Frontiers in Cellular and Infection Microbiology》上，其重要意义在于识别了Anaeroglobus作为可重复的候选风险相关信号，强调了方向复制而非模型可移植性，并指出需在生态失调（dysbiosis）框架下理解该信号，而非单一病原体主导。

**主要关键技术方法**

本研究主要依赖公共数据库中的三个唾液16S rRNA队列：发现队列来自Kasimoglu等（2020年，198样本，包括99对双胞胎），初级外部验证队列来自Li等（2021年，SCC0006，130样本，6–8岁混合牙列期），第二验证队列来自Yang等（2023年，36样本，5岁儿童）。属级特征经中心对数比变换（CLR）和协变量调整（年龄、性别、双胞胎类型、牙列期）后，采用逻辑回归模型分析，双胞胎对内分析以家庭相关性为聚类单元。支持证据来自同一生态位牙菌斑跨队列整合研究（Khan等，2024年，5个ECC牙菌斑研究，223样本），牙菌斑转录组数据集GSE31439（Peterson等，2014年，38样本），以及深龋牙髓单细胞数据集GSE185222（Opasawatchai等，2022年，6810个牙髓细胞）。统计方法包括基于DerSimonian–Laird随机效应模型的探索性合并效应分析、ComBat批次效应校正和Benjamini–Hochberg多重检验校正。

**研究结果**

**队列组成与基线特征**：主要分析共纳入327份儿童唾液样本，包括发现队列197份（151例龋病，46例无龋），初级外部验证队列130份（65例龋病，65例对照），第二验证队列36份。发现队列中，龋病组和无龋组在年龄、牙菌斑指数和牙龈指数上存在差异，但α多样性（Shannon指数和观测属数）无显著差异。

**发现队列中的关键龋病相关属**：经调整年龄、性别、双胞胎类型和牙列期后，Lactobacillus与龋病状态独立相关（β=0.566，q=0.017），Candidatus_Saccharibacteria_bacterium_UB2523也呈正相关（β=0.398，q=0.013），Anaeroglobus与总龋负担正相关（β=0.029，q=0.021）。在18对不一致双胞胎中，FDR校正后无属级特征显著。

**初级和第二独立外部验证结果**：在初级外部验证队列中，Anaeroglobus每增加一个标准差，调整后OR为3.452（95% CI：1.851–6.449，P<0.001）；Lactobacillus的OR为2.118（95% CI：1.049–4.282，P=0.035）。在第二独立验证队列中，Anaeroglobus的OR为3.009（95% CI：1.241–7.312，P=0.015），Lactobacillus的OR为1.929（95% CI：0.866–4.262，P=0.106）。初级外部模型AUC为0.782（近似95% CI：0.703–0.861）。

**区域亚组和探索性批次效应敏感性分析**：按区域分层后，Anaeroglobus在青岛亚组（n=96，OR=3.381）和广州亚组（n=34，OR=3.658）均方向一致；Lactobacillus在青岛亚组显著（OR=2.079），广州亚组未达显著（OR=2.310）。ComBat调整后AUC从0.782降至0.746，效应方向未变。

**同一生态位牙菌斑公共数据库：群落级和特征级支持**：跨队列牙菌斑研究（5个研究，223样本）显示，批次校正前ANOSIM R值在无龋组为0.64，在ECC组为0.49；校正后分别降至0.075和0.091，表明批次效应强烈影响群落解释。随机森林AUROC在物种级约0.85，属级约0.83，更稳定的属为Alloprevotella和Megasphaera。

**牙菌斑转录组功能支持**：牙菌斑转录组数据集GSE31439（38样本）显示，碳水化合物利用转录本约占总转录本10%，单糖/二糖代谢相关转录本约占15%，氧化应激相关转录本占应激反应转录本的50%–75%，超氧化物歧化酶相关转录本约占氧化应激转录本的22%。

**深龋牙髓单细胞数据的第三层宿主反应支持**：深龋牙髓单细胞数据集GSE185222（6810个细胞）显示，深龋改变牙髓免疫细胞组成，上调促炎、抗炎和矿化相关基因，并增强与巨噬细胞的细胞间通讯，提示宿主免疫-基质重塑。

**探索性合并效应分析**：基于青岛亚组、广州亚组和5岁队列的合并效应分析显示，Anaeroglobus的合并OR为3.335（95% CI：2.087–5.329，I2=0.0%），Lactobacillus的合并OR为2.082（95% CI：1.310–3.308，I2=0.0%），留一法分析支持稳定性。

**讨论与结论**

讨论部分强调，本研究的主要贡献并非提供一个可稳定跨场景应用的单一属模型，而是识别了Anaeroglobus在三个唾液层（发现阶段、初级外部验证队列和第二独立外部验证队列）中的相对一致方向复制信号。该发现表明该属在不同儿童队列中有一定概率重复出现，但“方向复制”不等同于“模型可迁移性”或立即临床可用性。基于属级唾液16S特征的模型仍仅表现出中等外部AUC，这与最新系统综述中儿童龋病风险预测模型普遍缺乏充分外部验证、仅中等判别能力和有限泛化能力的结果一致。因此，本研究更合适的定位是提供一个跨多个队列反复观察到的候选风险信号，而非已建立直接可推广的微生物组生物标志物模型。从病因学角度，该发现更适合在龋病相关生态失衡框架内解释，而非对单一优势致龋菌的重新确认。Anaeroglobus的意义更接近龋病生态转变的可重复信号，而非独立且充分的病因学证据。其生物学合理性应理解为辅助风险标志物，而非决定性病原体。研究局限性包括：未对牙菌斑数据集进行统一原始数据重处理；主要分析基于唾液，与局部病灶生态位存在差异；发现与验证队列在年龄、牙列期和数据预处理流程上存在差异；第二验证队列仅36样本，统计效力不足，应视为年龄外推支持；支持性证据层描述了一般性龋病生物膜和宿主反应背景，不能推断Anaeroglobus特异功能或宿主机制；模型评估主要基于AUC，校准和临床效用指标仍需未来个体水平和前瞻性验证。

**研究结论翻译**：综合公共儿童口腔微生物组数据集的分析表明，Anaeroglobus在发现阶段、130样本初级外部验证队列和36样本第二独立外部验证队列中方向可重复，探索性合并效应分析描述性提示合并随机效应OR为3.335（95% CI：2.087–5.329），Lactobacillus也显示稳定合并效应（OR=2.082，95% CI：1.310–3.308），而仅基于唾液16S属级特征的模型外部泛化能力仍有限。第二验证队列因样本量小应视为支持性年龄外推证据，合并效应估计应视为完成验证子集的探索性综合，而非完全独立外部可重复性的确定性证据。需明确区分，本研究中“多队列外部验证”仅指唾液发现层、初级外部层和第二独立外部层；添加的牙菌斑、牙菌斑转录组和深龋牙髓单细胞层属于支持性证据层，仅用于表明所观察到的信号更可能位于牙菌斑群落重构、增强的碳水化合物代谢和氧化应激适应以及宿主免疫-基质重塑的更广泛背景中。因此，Anaeroglobus更适合解释为需要进一步验证的候选风险相关辅助信号，而非即时可翻译的单一属诊断工具。