前驱期和确诊帕金森病肠道病毒组的鸟枪法宏基因组分析

Shotgun Metagenomic Profiling of the Gut Virome in Prodromal and Confirmed Parkinson's Disease

作者信息Deepika Dinesh, Xochitl C Morgan, Jordan Jensen, Kjetil Bjornevik, Michael A Schwarzschild, Alberto Ascherio, Curtis Huttenhower, Natalia Palacios
PMID42204574
发布时间2026-08
DOI10.1002/ana.78243

实验完整度

研究为前瞻性巢式病例对照设计,包含PD、pPD和健康对照的多组宏基因组分析,使用了新方法BAQLaVa进行病毒组特征分析,并通过多变量模型和随机森林验证,但未包含体外或动物等分子机制实验。

主要模型

PD患者粪便样本队列 pPD患者粪便样本队列

重点核对

样本来源:NHS和HPFS队列的粪便样本 病例定义:PD需经神经科医生确认,pPD依据便秘、pRBD和嗅觉评分等标准 样本采集:每位参与者提供2套粪便样本,平均间隔24.6个月 统计分析:MaAsLin2模型调整了批次、年龄、性别、饮食、吸烟、BMI、Bristol大便量表、采样日期等协变量 随机森林模型:采用交叉验证,将同一受试者的重复样本放在同一训练或测试集中

摘要

我们在护士健康研究和卫生专业人员随访研究中开展了一项巢式病例对照研究,以探讨肠道病毒组在帕金森病中的作用。我们采用了一种新的宏基因组病毒组分析流程BAQLaVa,对来自62名帕金森病患者、123名健康对照和90名前驱期帕金森病患者的粪便宏基因组数据进行了前瞻性分析。多变量线性模型鉴定出3个在帕金森病中丰度升高的病毒基因组bin:MVG081219(β=0.86,q=0.013)、MVG041501(β=0.95,q=0.048)、MVG081211(β=0.66,q=0.048),以及1个在帕金森病患者中丰度降低的病毒基因组bin MVG098915(β=-1.42,q=0.047)。这4个VGBs在前驱期帕金森病中也有类似的相关性。这项工作表明肠道病毒组有潜力作为未来帕金森病的生物标志物。

实验结论

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

研究问题
肠道病毒组的组成在帕金森病(PD)和前驱期帕金森病(pPD)中是否发生变化,以及这些变化是否具有作为PD生物标志物的潜力。
核心机制
文中未明确说明具体机制,但提示有尾噬菌体纲(Caudoviricetes)噬菌体可能通过与细菌宿主的裂解或溶原性相互作用,影响肠道微生物群落结构,从而参与肠道菌群失调及肠-脑轴信号调节。
主要证据
基于NHS和HPFS队列的鸟枪法宏基因组数据,使用BAQLaVa方法鉴定病毒组,通过MaAsLin2和随机森林模型,发现3个VGBs在PD中丰度升高,1个丰度降低,且在pPD中有类似关联。
研究意义
研究提示肠道病毒组可作为PD和pPD病因、风险和进展的潜在生物标志物,并强调了进一步研究Caudoviricetes与细菌相互作用及功能表征的必要性。

研究路径

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

1

样本收集与队列设计

建立具有明确病例定义的前瞻性PD、pPD及健康对照队列,收集粪便样本以分析肠道病毒组。

在NHS和HPFS队列中,纳入62名PD患者、123名健康对照和90名pPD参与者,每位参与者提供2套粪便样本,平均间隔24.6个月。PD由神经科医生确认,pPD根据便秘、pRBD和嗅觉评分等定义。

2

宏基因组测序与病毒组分析

从宏基因组测序数据中定量和分类肠道病毒基因组bin,获得病毒组丰度谱。

对粪便样本进行鸟枪法宏基因组测序,使用KneadData进行质控和宿主reads去除,随后使用BAQLaVa V0.1进行病毒组分析,基于核苷酸和蛋白质数据库进行双阶段映射,鉴定VGBs。

3

差异丰度分析

鉴定在PD或pPD与健康对照之间丰度有显著差异的VGBs。

使用MaAsLin2进行特征级关联分析,调整协变量,并以FDR q<0.05为显著性阈值。

4

总体病毒群落结构分析

评估PD、pPD和健康对照之间病毒组群落结构的差异。

使用Bray-Curtis距离和PERMANOVA检测β多样性,使用Shannon指数和线性混合模型检测α多样性。

5

随机森林分类器评估

评估病毒组作为生物标志物区分PD或pPD与健康对照的潜力。

使用scikit-learn的随机森林分类器,基于所有VGBs丰度构建分类模型,采用100次随机80/20交叉验证,并计算ROC AUC。

研究方法

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

产品清单

实验环节名称品牌货号
Megablast序列比对工具----
iPHoP噬菌体宿主预测工具----
肠道病毒组数据库----

关键环节

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

环节核对要点
样本采集与队列设计
样本来源:NHS和HPFS队列;PD诊断标准:需神经科医生确认;pPD定义标准:便秘、pRBD和嗅觉评分阈值;样本采集时间范围:2018年12月至2022年1月;两次样本采集间隔:平均24.6个月。
阅读提示:Methods中'PD Phenotypes'部分和Table 1
宏基因组测序与病毒组分析
测序平台和文库制备方法:文中未明确说明;数据质控流程:KneadData版本0.3;BAQLaVa版本0.1;数据库版本和存在性阈值:核苷酸和蛋白质特征覆盖≥50%。
阅读提示:Methods中'Bioinformatic Application for Quantification and Labeling of Viral Taxonomy'部分
差异丰度分析
协变量调整:批次、年龄、性别、地中海饮食评分、吸烟包年、BMI、Bristol大便量表、采样日期;随机效应:受试者ID;统计显著性阈值:q<0.05。
阅读提示:Methods中'Feature-Wise Analyses'部分
总体病毒群落结构分析
β多样性度量:Bray-Curtis距离;PERMANOVA的置换次数:文中未明确说明;α多样性指数:Shannon指数,并采用线性混合模型,随机效应为受试者。
阅读提示:Methods中'Overall Community Patterns of Microbial Variation'部分
随机森林分类器
实现包:scikit-learn默认参数;交叉验证:100次随机80/20分割,重复样本在同一折;性能指标:ROC AUC。
阅读提示:Methods中'Random Forest Classifier'部分