根管微生物组与根尖病变大小的关联:一项观察性鸟枪法宏基因组研究

Association Between the Root Canal Microbiome and Apical Lesion Size: An Observational Shotgun Metagenomic Study

作者信息Virginia Pérez-Carrasco, David Uroz-Torres, Ana Soriano-Lerma, Miguel Soriano, Jose Antonio García-Salcedo, Maria Teresa Arias-Moliz
PMID42227352
发布时间2026-10
DOI10.1111/iej.70190

实验完整度

观察性研究,含qPCR、鸟枪法宏基因组测序和生物信息分析,但无功能验证或机制实验。

主要模型

原发性根管感染患者的根管样本 小病变组(<3 mm) 大病变组(>7 mm)

重点核对

病变大小分组:小病变<3 mm,大病变>7 mm,由CBCT测量 样本量:最终纳入49例,小病变组27例,大病变组22例 细菌载量定量:16S rRNA基因qPCR,结果以Log10 DNA拷贝数表示 测序与注释:鸟枪法宏基因组测序,Bracken/Kraken2和HUMAnN3注释 差异分析:ALDEx2用于物种,edgeR用于KEGG通路,FDR校正,p<0.05

摘要

目的:旨在利用鸟枪法宏基因组测序表征原发性牙髓感染中微生物组的分类学和功能组成,并评估其与根尖周病变大小的关联。方法学:对诊断为根尖周炎的原发性根管感染样本进行鸟枪法测序分析。样本根据病变大小分为小病变(<3 mm)和大病变(>7 mm)。通过qPCR定量各组中的细菌DNA拷贝数。使用Bracken/Kraken2和HUMAnN3软件进行 taxonomic和功能注释。使用物种丰富度、Shannon、Simpson和Pielou指数测量α多样性。通过基于Bray-Curtis距离的主坐标分析评估研究组间细菌群落的相似性。使用ALDEx2包推断物种间差异,使用edgeR包分析KEGG通路。所有统计分析中,p<0.05被认为具有显著性。结果:共分析49个样本,27个为小病变,22个为大病变。物种丰富度和Shannon指数显示两组间存在差异,而Simpson和Pielou指数未见差异。两组间观察到不同的群落组成(PERMANOVA,p=0.0019)。三个物种在大病变样本中显著富集,即Filifactor alocis、Lachnospiraceae bacterium oral taxon 500和Olsenella uli,而另外三个物种在小病变样本中富集,即Acinetobacter baumannii、Acinetobacter pittii和Cutibacterium acnes。在功能上,苯甲酸盐、类黄酮和类固醇降解、鞘脂信号通路和蛋白酶体功能在大病变样本中富集。单萜生物合成、磷脂酶D信号、硫 relay系统和staurosporine生物合成在小病变中富集。结论:与大病变相比,大根尖周病变的牙齿携带更高的细菌载量并表现出更多样的微生物群落。两组间观察到物种水平分类组成的差异。在功能上,大病变富集与免疫逃逸和促炎活性相关的通路,而小病变的特征是与凋亡、代谢适应和抗炎过程相关的通路。这些发现表明,病变严重程度也受到微生物组调节宿主炎症的功能潜力的影响。

实验结论

提炼研究问题、关键发现与证据,快速把握文章的核心贡献。

研究问题
表征原发性牙髓感染中微生物组的分类学和功能组成,并评估其与根尖周病变大小的关联。
核心机制
大根尖周病变富集与免疫逃逸和促炎活性相关的细菌通路,而小病变富集与凋亡、代谢适应和抗炎过程相关的通路,提示微生物组功能潜力影响病变严重程度。
主要证据
49例样本的鸟枪法宏基因组测序和qPCR显示,大病变组细菌载量更高、α多样性更高,且物种和KEGG通路组成存在显著差异。
研究意义
表明根尖周疾病的进展和严重程度不仅由微生物组成决定,还受驻留微生物群放大或调节宿主炎症反应的功能能力影响。

研究路径

按研究推进顺序梳理实验设计、验证步骤与关键观察。

1

样本采集与分组

根据根尖周病变大小对原发性根管感染样本进行分组,以比较不同病变阶段的微生物组差异。

从原发性根管感染患者中采集样本,使用CBCT测量病变大小,将样本分为小病变(<3 mm)和大病变(>7 mm)两组。

2

细菌载量定量

比较小病变和大病变根管中的总细菌载量。

使用16S rRNA基因靶向qPCR定量各样本中的细菌DNA拷贝数。

3

鸟枪法宏基因组测序与生物信息分析

表征微生物组的分类学组成和功能通路,并比较两组间的差异。

对样本DNA进行鸟枪法测序,使用KneadData、Trimmomatic和Bowtie2进行质控和去宿主,使用Bracken/Kraken2进行分类注释,HUMAnN3进行功能注释,DIAMOND和UNIPROT进行翻译搜索。

4

多样性及群落组成分析

评估两组样本的α多样性和β多样性差异。

计算物种丰富度、Shannon、Simpson和Pielou指数,基于Bray-Curtis距离进行PCoA和PERMANOVA分析。

5

差异物种和功能通路分析

识别两组间显著差异的细菌物种和KEGG功能通路。

使用ALDEx2分析物种水平差异,edgeR分析KEGG通路差异,并进行Spearman相关性分析。

研究方法

按研究目的归类文中使用的方法,便于定位所需技术。

产品清单

实验环节名称品牌货号
Carestream CS Select CBCT设备Carestream Dental LLCCS 81003D
无菌车针Komet--
硫代硫酸钠Sigma-Aldrich--
无菌旋转锉Endogal--
无菌纸尖Endogal--
Tris-EDTA缓冲液Panreac Quimica--
QIAmp DNA迷你试剂盒Qiagen--
Nanodrop 2000紫外-可见分光光度计Thermo Fisher Scientific--
QuantStudio 6高通量实时PCR系统Thermo Fisher Scientific--
Power SYBR green PCR预混液Bio-Rad Laboratories--
MicroAmp Fast 96孔反应板Thermo Fisher Scientific--
DNA文库制备试剂盒Illumina1000000025416
Qubit 2.0荧光定量仪Thermo Fisher Scientific--
NovaSeq 6000测序仪Illumina--

关键环节

汇总复现实验时建议重点确认的条件及原文阅读提示。

环节核对要点
样本采集与分组
病变大小分类阈值:小病变<3 mm,大病变>7 mm;排除3–7 mm;最终样本量49例(小病变27例,大病变22例)
阅读提示:Methods 2.1 Subjects and Specimens;Results 3 及 Figure 1
细菌载量定量
qPCR引物序列、终浓度、循环条件、标准曲线及双重复设置
阅读提示:Methods 2.3 qPCR Assay
鸟枪法宏基因组测序
输入DNA量25–500 ng、测序平台NovaSeq 6000、总reads数1106 M、平均每样本20 M reads
阅读提示:Methods 2.4 Shotgun Sequencing and Bioinformatic Analysis
生物信息分析
分类注释软件Bracken/Kraken2,功能注释HUMAnN3,数据库UNIPROT,聚类阈值90%同源
阅读提示:Methods 2.4 Shotgun Sequencing and Bioinformatic Analysis
统计分析
α多样性指数、PCoA、PERMANOVA(9999次置换)、ALDEx2和edgeR差异分析、FDR校正、p<0.05
阅读提示:Methods 2.5 Statistical Analysis