肠道微生物群与心血管疾病的因果关系评价:双向孟德尔随机化分析
1.
2.
Causal relationship between gut microbiota and cardiovascular diseases: a bidirectional Mendelian randomization analysis
1.
2.
通讯作者: 闫奎坡,主任医师,博士;电子信箱:ykp19821122@163.com。
编委: 瞿麟平
收稿日期: 2025-04-17 接受日期: 2025-06-17 网络出版日期: 2025-12-28
| 基金资助: |
|
Corresponding authors: YAN Kuipo, E-mail:ykp19821122@163.com.
Received: 2025-04-17 Accepted: 2025-06-17 Online: 2025-12-28
| Fund supported: |
|
目的·通过孟德尔随机化(Mendelian randomization,MR)方法探讨肠道微生物群与心血管疾病(cardiovascular disease,CVD)之间的因果关系。方法·使用MiBioGen联盟提供的肠道微生物群数据(n=18 340)和IEU Open GWAS数据库提供的4种CVD(心房颤动1 030 836例、冠状动脉疾病547 261例、高血压20 526例、心力衰竭977 323例)相关的遗传位点作为工具变量。研究采用逆方差加权法(inverse variance weighted,IVW)作为主要研究方法。同时,使用Cochran's Q检验评估遗传工具变量的异质性,MR-Egger截距检验评估水平多效性,留一法评估作为工具变量的单核苷酸多态性(single-nucleotide polymorphism,SNP)对暴露和结局因果关系影响的敏感性。采用MR Steiger检验验证肠道微生物群与CVD之间的因果方向。结果·IVW法的研究结果表明:Victivallales(OR=0.939)、霍氏菌属(OR=0.939)、厌氧链球菌属(OR=0.922)、双歧杆菌科(OR=0.916)、黏胶球形菌纲(OR=0.936)、臭气杆菌属(OR=0.909)、Intestinibacter(OR=0.933)、黏胶球形菌门(OR=0.926)和双歧杆菌目(OR=0.916)对心房颤动表现为保护因素,而链状杆菌属(OR=1.057)、毛螺菌科UCG008(OR=1.051)、链球菌属(OR=1.089)和Victivallis(OR=1.038)则为危险因素;乳杆菌目(OR=0.919)和副拟杆菌属(OR=0.866)是冠状动脉疾病的保护因素,而韦荣球菌科(OR=1.065)、Lachnoclostridium(OR=1.093)、毛螺菌科(OR=1.094)、草酸杆菌属(OR=1.062)、臭气杆菌属(OR=1.160)是危险因素;Mollicutes RF9(OR=0.851),红椿菌纲(OR=0.803)、目(OR=0.803)、科(OR=0.803),以及Intestinibacter(OR=0.819)是高血压的保护因素,而Christensenellaceae R7 group(OR=1.218)、脱硫弧菌属(OR=1.167)和消化球菌科(OR=1.230)是危险因素;芽孢杆菌目(OR=0.955)和厌氧链球菌属(OR=0.899)是心力衰竭的保护因素,而瘤胃球菌UCG009(OR=1.107)、Eubacterium oxidoreducens group(OR=1.117)、月形单胞菌目(OR=1.106)、阴性杆菌目(OR=1.107)、Eubacterium eligens group(OR=1.139)和解黄酮菌属(OR=1.144)是危险因素。Cochran's Q检验显示,与CVD存在因果关系的肠道微生物群的SNP之间不存在异质性(均P>0.05);基因多效性检验未发现多效性(均P>0.05);留一法敏感性分析证实研究结果的稳健性。MR Steiger方向性检验结果支持肠道微生物群作为暴露、CVD作为结局的因果方向。结论·部分肠道微生物群对CVD存在显著的因果效应;改变其丰度可能影响CVD风险,这为基于微生物群的干预策略提供了潜在靶点。
关键词:
Objective ·To investigate the causal relationship between gut microbiota and cardiovascular diseases (CVDs) using Mendelian randomization (MR). Methods ·Instrumental variables included genetic loci from gut microbiota data provided by the MiBioGen consortium (n=18 340) and CVD data from the IEU Open GWAS database, covering four CVD types: atrial fibrillation (n=1 030 836), coronary artery disease (n=547 261), hypertension (n=20 526), and heart failure (n=977 323). The inverse variance weighted (IVW) method was employed as the primary analytical approach. Additionally, Cochran's Q test was used to assess heterogeneity of genetic instruments, the MR-Egger intercept test to evaluate horizontal pleiotropy, and leave-one-out analysis to examine the sensitivity of single-nucleotide polymorphisms (SNPs) on the exposure-outcome causal relationship. The MR Steiger test was applied to validate the causal direction between gut microbiota and CVDs. Results ·The IVW analysis indicated that Victivallales (OR=0.939), Howardella (OR=0.939), Anaerostipes (OR=0.922), Bifidobacteriaceae (OR=0.916), Lentisphaeria (OR=0.936), Odoribacter (OR=0.909), Intestinibacter (OR=0.933), Lentisphaerae (OR=0.926), and Bifidobacteriales (OR=0.916) were protective factors against atrial fibrillation, while Catenibacterium (OR=1.057), Lachnospiraceae UCG008 (OR=1.051), Streptococcus (OR=1.089), and Victivallis (OR=1.038) were risk factors. For coronary artery disease, Lactobacillales (OR=0.919) and Parabacteroides (OR=0.866) were protective factors, while Veillonellaceae (OR=1.065), Lachnoclostridium (OR=1.093), Lachnospiraceae (OR=1.094), Oxalobacter (OR=1.062), and Odoribacter (OR=1.160) were risk factors. For hypertension, Mollicutes RF9 (OR=0.851), Coriobacteriia (OR=0.803), Coriobacteriales (OR=0.803), Coriobacteriaceae (OR=0.803), and Intestinibacter (OR=0.819) were protective factors, while Christensenellaceae R7 group (OR=1.218), Desulfovibrio (OR=1.167) ,and Peptococcaceae (OR=1.230) were risk factors. For heart failure, Bacillales (OR=0.955) and Anaerostipes (OR=0.899) were protective factors, while Ruminococcus UCG009 (OR=1.107), Eubacterium oxidoreducens group (OR=1.117), Selenomonadales (OR=1.106), Negativicutes (OR=1.107), Eubacterium eligens group (OR=1.139), and Flavonifractor (OR=1.144) were risk factors. Cochran's Q test showed no heterogeneity among SNPs of gut microbiota causally associated with CVDs (all P>0.05). The pleiotropy test found no evidence of horizontal pleiotropy (all P>0.05). Leave-one-out sensitivity analysis confirmed the robustness of the results. The MR Steiger directionality test supported the causal direction from gut microbiota (exposure) to CVDs (outcome). Conclusion ·Some gut microbiota have significant causal effects on CVDs; altering their abundance may influence CVD risk, providing potential targets for microbiota-based intervention strategies.
Keywords:
本文引用格式
马会华, 闫奎坡, 刘刚, 徐亚洲, 张磊, 孙彦琴.
MA Huihua, YAN Kuipo, LIU Gang, XU Yazhou, ZHANG Lei, SUN Yanqin.
心血管疾病(cardiovascular disease,CVD)是全球发病和死亡的首要原因。一项研究[1]结果发现,1990至2019年的全球流行病学数据显示,CVD的疾病负担呈显著上升趋势:发病人数从1990年的3 131万例增至2019年的5 545万例,增幅77.1%;死亡人数从1 207万例升到1 856万例,增幅53.8%。在我国,随着社会经济的快速发展和生活方式的显著变化,特别是人口老龄化和城镇化进程的加速,CVD危险因素对居民健康的影响日益突出。流行病学调查[2]显示,我国现有CVD患病人数约3.3亿。越来越多关于生活方式的研究[3-4]表明,不同饮食习惯可能导致肠道微生物群水平及其代谢物浓度的变化。TELLE-HANSEN等[5]的研究表明,在膳食中以不饱和脂肪酸代替饱和脂肪酸能够调节人体肠道微生物群中相关菌群的丰度,并进一步降低血清胆固醇水平,最终有助于降低患CVD的风险。其他几项研究[6-7]表明,由肠道细菌产生的三甲胺N-氧化物(trimethylamine-N-oxide,TAMO)的变化与心脏的不同病理状况密切相关,包括高血压及动脉粥样硬化。因此,TAMO不仅被认为是CVD的一种新的生物标志物,而且被认为是慢性疾病的治疗靶点[8]。
肠道微生物群(gut microbiota)主要由细菌构成,还包括真菌和病毒。它在维持宿主健康方面发挥多重生物学功能,参与肠道黏膜屏障的形成与调控,调节营养物质的摄取与代谢,促进免疫系统的发育与成熟,以及抑制病原微生物的定植与繁殖[9-12]。大量研究[13-14]报道了肠道微生物群与疾病之间的相互作用,如胃肠道肿瘤、脑缺血再灌注损伤、肝纤维化和CVD。事实上,适当的肠道微生物群结构和代谢物功能对于维持体内平衡至关重要。如果肠道微生物群失衡,就会导致动脉粥样硬化、高血压、心力衰竭(heart failure,HF)和心律失常[15-16]。尽管许多研究报道了肠道微生物群的影响,但尚未确定与CVD相关的所有肠道微生物群。由于研究中的仪器、设备和实验对象各不相同,因此从CVD研究中得出的结论没有达成共识。因此,肠道微生物群与CVD之间的因果关系需要进一步研究。
本研究旨在通过全基因组关联研究(genome-wide association study,GWAS)[17]和孟德尔随机化(Mendelian randomization,MR)方法[18],探讨肠道微生物群与CVD之间的因果关系。MR分析通过利用常见遗传变异与可改变环境暴露因素之间的内在关联特性,已成为探索环境暴露与疾病潜在因果关系的常用方法[19-21]。两样本MR分析可以整合来自独立GWAS的单核苷酸多态性(single nucleotide polymorphism,SNP)暴露与SNP结果的关联,将其合并为单一因果效应估计值。随着肠道微生物群和疾病的GWAS数量迅速增加,大规模汇总统计变得更为广泛,支持两样本MR分析,并显著提高统计功效[22-23]。
1 对象与方法
1.1 研究设计
本研究采用两样本MR分析方法,评估肠道微生物群与CVD之间的因果关系。基于双向MR分析框架,本研究以遗传变异作为工具变量(instrumental variable),系统探讨肠道微生物与CVD的潜在因果关系。为确保研究结果的可靠性,MR分析需满足以下3个核心假设。①相关性假设:最终纳入的工具变量必须与暴露因素(肠道微生物群)的关联具有统计学意义,以确保工具变量能够有效地反映暴露因素的变化。②独立性假设:工具变量与影响CVD的混杂因素无关。③排他性假设:工具变量仅能通过暴露因素(肠道微生物群)影响结局变量(CVD)。满足上述假设是确保MR分析结果有效性的关键前提。
1.2 数据来源和研究人群
1.2.1 暴露数据
本研究的暴露因素和结局变量均来自公共数据库。暴露因素为人类肠道微生物,其GWAS统计数据来源于MiBioGen联盟数据库(
1.2.2 结局数据
本研究中,CVD定义为心房颤动(atrial fibrillation,AF)、冠状动脉疾病(coronary artery disease,CAD)、高血压(hypertension)和HF 4种疾病,疾病数据均来源于IEU Open GWAS项目(
表1 CVD的GWAS数据集
Tab 1
| GWAS data ID | Disease | Sample size/n | Case/n | Control/n | Population |
|---|---|---|---|---|---|
| ebi-a-GCST006414 | AF | 1 030 836 | 60 620 | 970 216 | European |
| ebi-a-GCST005195 | CAD | 547 261 | 122 733 | 424 528 | European |
| ebi-a-GCST008036 | Hypertension | 20 526 | 11 863 | 8 663 | European |
| ebi-a-GCST009541 | HF | 977 323 | 47 309 | 930 014 | European |
1.3 工具变量选择
将细菌分类群分为5个层次(门、纲、目、科、属)进行分析,每个分类群作为一个特征。为了确保肠道微生物组与CVD风险因果关系结论的真实性和准确性,本研究采用以下质量控制步骤来选择最佳的工具变量。首先,筛选出GWAS汇总数据中与全基因组显著性(P<1×10-5)遗传相关的SNP作为工具变量。其次,基于千人基因组样本数据进行连锁不平衡(linkage disequilibrium,LD)分析,运行参数设置为r2<0.001,窗口大小10 000 kb。其次,通过计算F统计量确保遗传工具的暴露强度。F>10表示不存在弱工具变量偏差,剔除F≤10的工具变量。F值计算公式如下:F=β2/se2[24]。β为等位基因效应值,表示工具变量与暴露因素的回归系数,反映其关联强度;se为标准误(standard error),可衡量估计的精确度。最后,本研究将所有SNP的等位基因与GRCh37参考基因组进行对比,并剔除存在链方向歧义的回文SNP以及等位基因信息模糊或重复的位点,确保遗传工具变量的编码准确性。
1.4 MR分析
统计分析使用R语言软件(R studio 4.3.3)中的TwoSampleMR(0.6.2)和MR-PRESSO(1.0)软件包进行。P<0.05被认为是潜在因果关系的显著性证据。本研究采用了逆方差加权(inverse variance weighting,IVW)法、MR-Egger法、加权中位法、加权模型、简单模式以及多种敏感性分析。据报道,在某些条件下,IVW方法比其他方法更为稳健。因此,多个工具变量的结果主要基于IVW法。IVW法依据各遗传工具变量估计精度的差异进行加权汇总,对更精确的估计赋予更高的权重,以此优化整体因果效应的估计结果[25-26];其余4种方法作为补充来进行全面评估。如果只有一种SNP作为工具变量,则暴露与结局之间的因果关系将通过Wald比值来估计[24,27]。
1.5 敏感性分析
为了严谨地评估肠道微生物群与CVD之间因果效应的稳健性,本研究进行了敏感性分析。运用Cochran's Q检验和MR-Egger截距测试来评估异质性和水平多效性。此外,为进一步检验单个SNP对整体因果效应估计结果的影响程度,开展了留一法(leave-one-out)敏感性分析。通过依次排除每一个SNP,重新计算剩余SNP组合下肠道微生物与CVD因果关系的估计值,并与包含所有SNP时的估计结果进行对比,以此判断单个SNP是否对整体估计结果产生不成比例的影响。
1.6 MR Steiger方向性检验
1.7 生物信息学分析
为系统阐释与CVD存在潜在因果关系的肠道微生物群的生物学功能,本研究在MR分析的基础上,对筛选得到的菌群特征进行基因本体论(Gene Ontology,GO)和京都基因组百科全书(Kyoto Encyclopedia of Genes and Genomes,KEGG)通路富集分析,为肠道微生物群与CVD因果关联提供潜在的生物学机制解释。
2 结果
2.1 工具变量筛选
根据显著性水平P<1×10-5,本研究在门、纲、目、科和属的层次上分别发现1、3、7、5和20个(共36个)肠道微生物群的394个SNP作为工具变量,并与肠道微生物群和CVD之间的因果关系相关。此外,所有工具变量的F统计值均大于10,表明所选SNP均为强工具变量,本研究不存在弱工具变量偏倚。通过MR-PRESSO离群值检验和MR-Egger回归鉴定去除多效性SNP后,工具变量未显示水平多效性偏倚(表2)。
表2 基于IVW法分析肠道微生物与CVD的因果关系
Tab 2
| Outcome | Exposure | MR method | SNP/n | β | se | OR (95%CI) | P value | Correct causal direction | P value (Steiger test) |
|---|---|---|---|---|---|---|---|---|---|
| AF | |||||||||
| Catenibacterium | IVW | 5 | 0.05 | 0.025 | 1.057 (1.005‒1.112) | 0.030 | True | 4.85×10-6 | |
| Victivallales | IVW | 8 | -0.06 | 0.023 | 0.939 (0.896‒0.984) | 0.008 | True | 1.52×10-6 | |
| Howardella | IVW | 10 | -0.06 | 0.019 | 0.939 (0.904‒0.977) | 0.001 | True | 3.04×10-5 | |
| Lachnospiraceae UCG008 | IVW | 11 | 0.04 | 0.025 | 1.051 (1.001‒1.104) | 0.047 | True | 2.52×10-7 | |
| Anaerostipes | IVW | 13 | -0.08 | 0.037 | 0.922 (0.857‒0.993) | 0.030 | True | 5.20×10-6 | |
| Bifidobacteriaceae | IVW | 11 | -0.08 | 0.034 | 0.916 (0.855‒0.980) | 0.011 | True | 1.42×10-5 | |
| Lentisphaeria | IVW | 8 | -0.06 | 0.024 | 0.936 (0.887‒0.983) | 0.003 | True | 3.96×10-6 | |
| Streptococcus | IVW | 15 | 0.08 | 0.039 | 1.089 (1.008‒1.177) | 0.029 | True | 9.37×10-6 | |
| Victivallis | IVW | 10 | 0.03 | 0.018 | 1.038 (1.001‒1.077) | 0.043 | True | 4.62×10-6 | |
| Odoribacter | IVW | 7 | -0.09 | 0.046 | 0.909 (0.831‒0.996) | 0.040 | True | 1.81×10-5 | |
| Intestinibacter | IVW | 15 | -0.06 | 0.031 | 0.933 (0.877‒0.993) | 0.028 | True | 2.28×10-6 | |
| Lentisphaerae | IVW | 9 | -0.07 | 0.022 | 0.926 (0.886‒0.968) | <0.001 | True | 5.06×10-6 | |
| Bifidobacteriales | IVW | 11 | -0.08 | 0.034 | 0.916 (0.855‒0.980) | 0.011 | True | 1.35×10-5 | |
| CAD | |||||||||
| Lactobacillales | IVW | 15 | -0.08 | 0.032 | 0.919 (0.862‒0.980) | 0.010 | True | 9.28×10-7 | |
| Veillonellaceae | IVW | 19 | 0.06 | 0.026 | 1.065 (1.011‒1.122) | 0.018 | True | 1.21×10-5 | |
| Parabacteroides | IVW | 6 | -0.14 | 0.050 | 0.866 (0.784‒0.956) | 0.004 | True | 6.31×10-5 | |
| Lachnospiraceae | IVW | 17 | 0.09 | 0.038 | 1.094 (1.012‒1.183) | 0.023 | True | 7.23×10-6 | |
| Lachnoclostridium | IVW | 13 | 0.08 | 0.039 | 1.093 (1.011‒1.181) | 0.025 | True | 9.14×10-5 | |
| Oxalobacter | IVW | 11 | 0.05 | 0.020 | 1.062 (1.019‒1.106) | 0.004 | True | 6.05×10-6 | |
| Odoribacter | IVW | 7 | 0.14 | 0.047 | 1.160 (1.056‒1.275) | 0.001 | True | 6.42×10-7 | |
| Hypertension | |||||||||
| Mollicutes RF9 | IVW | 13 | -0.16 | 0.068 | 0.851 (0.743‒0.974) | 0.019 | True | 7.23×10-5 | |
| Peptococcaceae | IVW | 9 | 0.20 | 0.086 | 1.230 (1.038‒1.458) | 0.016 | True | 2.07×10-5 | |
| Christensenellaceae R7 group | IVW | 10 | 0.19 | 0.100 | 1.218 (1.001‒1.483) | 0.049 | True | 9.97×10-5 | |
| Coriobacteriales | IVW | 14 | -0.21 | 0.085 | 0.803 (0.679‒0.950) | 0.010 | True | 0.001 | |
| Coriobacteriia | IVW | 14 | -0.21 | 0.085 | 0.803 (0.679‒0.950) | 0.010 | True | 0.002 | |
| Desulfovibrio | IVW | 10 | 0.15 | 0.062 | 1.167 (1.033‒1.318) | 0.012 | True | 0.007 | |
| Coriobacteriaceae | IVW | 14 | -0.21 | 0.085 | 0.803 (0.679‒0.950) | 0.010 | true | 0.007 | |
| Intestinibacter | IVW | 15 | -0.19 | 0.083 | 0.819 (0.696‒0.965) | 0.010 | True | 0.015 | |
| HF | |||||||||
| Ruminococcaceae UCG009 | IVW | 12 | 0.06 | 0.030 | 1.107 (1.009‒1.137) | 0.022 | True | 6.12×10-6 | |
| Eubacterium oxidoreducens group | IVW | 4 | 0.11 | 0.043 | 1.117 (1.026‒1.215) | 0.010 | True | 3.26×10-5 | |
| Bacillales | IVW | 9 | -0.04 | 0.022 | 0.955 (0.913‒0.998) | 0.010 | True | 5.84×10-6 | |
| Selenomonadales | IVW | 12 | 0.10 | 0.044 | 1.106 (1.013‒1.208) | 0.023 | True | 1.53×10-5 | |
| Anaerostipes | IVW | 13 | -0.10 | 0.043 | 0.899 (0.825‒0.974) | 0.013 | True | 3.77×10-6 | |
| Negativicutes | IVW | 12 | 0.10 | 0.044 | 1.107 (1.014‒1.208) | 0.023 | True | 4.48×10-5 | |
| Eubacterium eligens group | IVW | 7 | 0.13 | 0.057 | 1.139 (1.019‒1.274) | 0.022 | True | 6.43×10-6 | |
| Flavonifractor | IVW | 5 | 0.13 | 0.053 | 1.144 (1.031‒1.270) | 0.011 | True | 1.06×10-5 |
2.2 肠道微生物群对4种CVD发展的影响
2.2.1 AF
IVW方法研究结果(表2)表明,13种肠道微生物群与AF有关。Victivallales(OR=0.939)、霍氏菌属(Howardella;OR=0.939)、厌氧链球菌属(Anaerostipes;OR=0.922)、双歧杆菌科(Bifidobacteriaceae;OR=0.916)、黏胶球形菌纲(Lentisphaeria;OR=0.936)、臭气杆菌属(Odoribacter;OR=0.909)、Intestinibacter(OR=0.933)、黏胶球形菌门(Lentisphaerae;OR=0.926)和双歧杆菌目(Bifidobacteriales;OR=0.916)对AF具有保护作用。另一方面,链状杆菌属(Catenibacterium;OR=1.057)、毛螺菌科UCG008(Lachnospiraceae UCG008;OR=1.051)、链球菌属(Streptococcus;OR=1.089)和Victivallis(OR=1.038)等菌群与AF存在反向因果关系,表明这些微生物菌群比例升高会增加AF的风险。此外,还采用MR-Egger、加权中位数、简单模式与加权模式4种方法加以验证。结果(图1)显示,4种辅助方法得出的效应方向与IVW法大体上保持一致。
图1
图1
基于4种MR法分析肠道微生物与AF因果关系的森林图
Fig 1
Forest plot of causal relationship between gut microbiota and AF based on 4 MR methods
2.2.2 CAD
图2
图2
基于4种MR法分析肠道微生物与CAD因果关系的森林图
Fig 2
Forest plot of causal relationship between gut microbiota and CAD based on 4 MR methods
2.2.3 高血压
从高血压与肠道微生物群的MR分析中可知Mollicutes RF9(OR=0.851),红椿菌纲(Coriobacteriia;OR=0.803)、目(Coriobacteriales;OR=0.803)、科(Coriobacteriaceae;OR=0.803),以及Intestinibacter(OR=0.819)是高血压的保护因素;而Christensenellaceae R7 group(OR=1.218)、脱硫弧菌属(Desulfovibrio;OR=1.167)和消化球菌科(Peptococcaceae;OR=1.230)菌群的增加是高血压的危险因素(表2)。值得注意的是,MR-Egger方法的结果显示,红椿菌纲、目、科和Intestinibacter的效应方向与IVW法相反,而加权中位数、简单模式与加权模式得到的效应方向与IVW方法一致。其余菌群采用4种方法得到的效应方向均与IVW法一致(图3)。
图3
图3
基于4种MR法分析肠道微生物与高血压因果关系的森林图
Fig 3
Forest plot of causal relationship between gut microbiota and hypertension based on 4 MR methods
2.2.4 HF
MR分析表明,芽孢杆菌目(Bacillales;OR=0.955)和厌氧链球菌属(OR=0.899)对HF具有保护作用(表2)。本研究还发现一些肠道菌群对HF而言是危险因素(表2),包括瘤胃球菌UCG009(Ruminococcaceae UCG009;OR=1.107)、Eubacterium oxidoreducens group(OR=1.117)、月形单胞菌目(Selenomonadales;OR=1.106)、阴性杆菌目(Negativicutes;OR=1.107)、Eubacterium eligens group(OR=1.139)和解黄酮菌属(Flavonifractor;OR=1.144)(表2)。另外4种方法得出的效应方向均与IVW法一致(图4)。
图4
图4
基于4种MR法分析肠道微生物与HF因果关系的森林图
Fig 4
Forest plot of causal relationship between gut microbiota and HF based on 4 MR methods
2.3 敏感性分析
表3 肠道微生物对CVD效应的异质性及水平多效性分析
Tab 3
| Outcome | Exposure | Cochran's Q | MR-Egger intercept | |||||
|---|---|---|---|---|---|---|---|---|
| MR-Egger | IVW | Intercept | P value | |||||
| Q value | P value | Q value | P value | |||||
| AF | ||||||||
| Catenibacterium | 0.370 | 0.946 | 0.488 | 0.974 | 0.010 | 0.753 | ||
| Victivallales | 3.948 | 0.683 | 4.155 | 0.761 | 0.005 | 0.665 | ||
| Howardella | 7.364 | 0.497 | 7.622 | 0.572 | -0.006 | 0.624 | ||
| Lachnospiraceae UCG008 | 8.051 | 0.528 | 8.162 | 0.612 | 0.004 | 0.747 | ||
| Anaerostipes | 12.334 | 0.339 | 12.561 | 0.402 | -0.003 | 0.661 | ||
| Bifidobacteriaceae | 8.821 | 0.453 | 9.111 | 0.521 | 0.004 | 0.603 | ||
| Lentisphaeria | 3.948 | 0.683 | 4.155 | 0.761 | 0.005 | 0.665 | ||
| Streptococcus | 20.126 | 0.092 | 20.223 | 0.123 | 0.002 | 0.807 | ||
| Victivallis | 7.007 | 0.535 | 8.153 | 0.518 | 0.019 | 0.315 | ||
| Odoribacter | 3.621 | 0.605 | 4.634 | 0.591 | 0.010 | 0.360 | ||
| Intestinibacter | 14.962 | 0.309 | 17.041 | 0.253 | 0.010 | 0.203 | ||
| Lentisphaerae | 7.455 | 0.383 | 7.841 | 0.449 | 0.007 | 0.566 | ||
| Bifidobacteriales | 8.821 | 0.454 | 9.111 | 0.521 | 0.004 | 0.603 | ||
| CAD | ||||||||
| Lactobacillales | 13.551 | 0.406 | 14.344 | 0.424 | -0.005 | 0.398 | ||
| Veillonellaceae | 14.934 | 0.528 | 1.741 | 0.495 | 0.005 | 0.241 | ||
| Parabacteroides | 4.608 | 0.329 | 4.742 | 0.448 | 0.005 | 0.750 | ||
| Lachnospiraceae | 21.818 | 0.112 | 24.485 | 0.079 | -0.008 | 0.195 | ||
| Lachnoclostridium | 5.378 | 0.911 | 6.793 | 0.870 | 0.011 | 0.259 | ||
| Oxalobacter | 3.473 | 0.942 | 3.692 | 0.960 | -0.006 | 0.650 | ||
| Odoribacter | 4.080 | 0.537 | 5.288 | 0.507 | 0.012 | 0.321 | ||
| Hypertension | ||||||||
| Mollicutes RF9 | 3.784 | 0.975 | 0.711 | 0.877 | <0.001 | 0.957 | ||
| Peptococcaceae | 5.131 | 0.643 | 0.331 | 0.178 | 0.055 | 0.358 | ||
| Christensenellaceae R7 group | 6.448 | 0.597 | 7.960 | 0.538 | 0.008 | 0.575 | ||
| Coriobacteriales | 8.827 | 0.717 | 9.577 | 0.728 | 0.012 | 0.575 | ||
| Coriobacteriia | 8.827 | 0.717 | 9.577 | 0.728 | 0.013 | 0.145 | ||
| Desulfovibrio | 3.323 | 0.912 | 7.960 | 0.538 | 0.222 | 0.253 | ||
| Coriobacteriaceae | 8.827 | 0.717 | 9.577 | 0.728 | -0.020 | 0.403 | ||
| Intestinibacter | 19.114 | 0.119 | 21.572 | 0.087 | -0.020 | 0.403 | ||
| HF | ||||||||
| Ruminococcaceae UCG009 | 8.939 | 0.537 | 10.956 | 0.446 | 0.017 | 0.185 | ||
| Eubacterium oxidoreducens group | 2.451 | 0.293 | 2.483 | 0.478 | -0.002 | 0.887 | ||
| Bacillales | 2.705 | 0.910 | 2.870 | 0.942 | 0.006 | 0.697 | ||
| Selenomonadales | 4.459 | 0.924 | 5.180 | 0.922 | -0.007 | 0.415 | ||
| Anaerostipes | 9.973 | 0.532 | 10.057 | 0.610 | -0.002 | 0.777 | ||
| Negativicutes | 4.459 | 0.924 | 5.180 | 0.922 | -0.007 | 0.415 | ||
| Eubacterium eligens group | 1.547 | 0.907 | 6.428 | 0.376 | -0.038 | 0.078 | ||
| Flavonifractor | 2.919 | 0.404 | 2.919 | 0.404 | 0.002 | 0.903 | ||
2.4 MR Steiger方向性检验
本研究通过MR Steiger方向性检验结果表明,所有纳入分析的遗传工具变量对肠道微生物的表型解释量均显著大于对CVD的解释量Steiger test P值<0.05,且全部支持正确的方向;该结果充分肯定了本研究中肠道微生物为暴露、CVD为结局这一因果关系假设(表2)。
2.5 KEGG和GO富集途径分析
图5
图5
影响CVD的肠道微生物群所携带基因的GO功能和KEGG通路富集分析
Note: A. Biological process (GO). B. Cell component (GO). C. Molecular function (GO). D. KEGG analysis.
Fig 5
GO function and KEGG pathway enrichment analysis of genes carried by gut microbiota affecting CVD
3 讨论
本研究进行了双向MR分析,以确定肠道微生物群和CVD之间的因果关系,为今后基于肠道微生物群预防或治疗CVD的研究提供参考。已有研究[29-31]证实,肠道中的微生物群落是饮食中的关键“代谢过滤器”,可以将常见营养素转化为代谢物,并且特定的微生物相关代谢物,如氧化三甲胺(trimethylamine-N-oxide,TMAO)、短链脂肪酸(short-chain fatty acids,SCFA)和次级胆汁酸,已被证明会影响CVD进展。GUO等[32]的研究表明,间歇性禁食可以改变肠道微生物群落的组成及其代谢途径,这与减少心脏代谢风险因素有关;这提示饮食习惯对维持心脏健康至关重要,不同的饮食方式能够改变肠道微生物的丰度及其代谢物的浓度。KOREN等[33]的早期测序研究表明,人类动脉粥样硬化斑块含有细菌DNA,推测微生物群可能与动脉粥样硬化有关,但尚未确定斑块中细菌DNA是否来自动脉壁内的活细菌。
AF是导致心脑血管疾病的重要因素之一。本研究的结果表明,9种肠道微生物群对AF具有保护作用,而4种肠道微生物群则是危险因素。因此,探索肠道微生物群在AF预防和治疗中的潜力具有重要意义。双歧杆菌(目、科)的增加可能降低AF风险。双歧杆菌能够增强胃肠道屏障功能并产生SCFA等有益代谢物[34]。此外,双歧杆菌是人体主要的益生菌之一,已被证明可以通过多种机制改善肠道屏障功能并预防高渗透性,包括Toll样受体2介导的免疫调节和抗炎作用[35]。值得注意的是,益生菌可能通过抑制与菌群失调相关的AF风险因素(如CAD、高血压和HF)提供保护作用[36-38],但是益生菌是否对患者AF具有直接的保护作用需要进一步的研究。近年来的几项研究[39-41]表明,肠道微生物群失调引起的肠道屏障损伤可能会诱发炎症反应,导致心房电活动失调和结构重塑,从而影响AF的发生和发展。而Victivallales、霍氏菌属、厌氧链球菌属、黏胶球形菌纲、黏胶球形菌门、臭气杆菌属和Intestinibacter对AF的影响作用目前没有直接的研究证据,但是其代谢产物(如SCFA、丁酸、丙酸)具有显著抗炎症作用,可能对AF产生影响。SCFA通过激活核因子κB(nuclear factor κB,NF-κB)信号通路,减少炎症因子的释放;同时还可促进调节性T细胞的分化,增强免疫耐受性,抑制过度的自身免疫反应[35-36]。丁酸盐作为一种具有抗炎特性的代谢产物,可通过调节巨噬细胞M1和M2表型的平衡,促进神经重塑,从而改善心脏功能并增强电稳定性,进而减少心律失常的发生[37]。在老年人群中,AF是一种常见病,不仅会导致HF,还会增加卒中的风险;如果可以通过检测肠道微生物群分布预测AF的发病,有利于及时干预,减少疾病和并发症的风险。
对于CAD,本研究发现乳杆菌和Parabacteroides具有保护作用。血脂异常是CAD的主要危险因素;已有动物研究[42]表明,乳杆菌能够改善高胆固醇和高血脂,在一定程度上减缓动脉粥样硬化的发展。在人体研究[43-45]中,补充2种益生菌菌株——弯曲乳杆菌HY 7601和植物乳杆菌KY 1032,能够降低甘油三酯水平,并提升载脂蛋白A-V水平和增大低密度脂蛋白胆固醇颗粒的尺寸。此外,本研究还发现,韦荣球菌科、Lachnoclostridium、毛螺菌科、草酸杆菌属和臭气杆菌属会增加CAD的风险。其中,韦荣球菌科、Lachnoclostridium及毛螺菌科均属于厚壁菌门。该门部分菌种是TMAO的关键生产者之一,而TMAO是已知的、可促进动脉粥样硬化形成的重要代谢物[46]。草酸杆菌作为革兰阴性菌,其携带的脂多糖(lipopolysaccharide,LPS)可透过受损或完整的肠道上皮屏障进入血液循环,通过驱动全身性低度炎症来促进动脉粥样硬化[47]。本研究通过MR分析发现,臭气杆菌属对AF的发生表现出保护性因果效应,却是CAD的危险因素。臭气杆菌属对AF发挥保护作用可能归因于SCFA丁酸盐[48]的作用,丁酸盐在动脉粥样硬化中同样发挥着积极作用[49];那么为何臭气杆菌属会促进CAD的发生发展,这有待进一步的实验性研究来证实。
高血压是全球范围内最常见的心血管疾病之一,然而肠道微生物群失调是否参与其发生发展目前尚不明确。在本研究中,我们鉴定出5种对心血管健康具有潜在保护作用的细菌。尽管目前尚无直接证据表明红椿菌的增加对高血压具有保护作用,但该菌在调节宿主葡萄糖代谢方面发挥重要作用,并对糖尿病患者表现出保护效应。在肠道中,红椿菌参与多种宿主生物学功能,包括胆汁酸与脂质代谢,促进脂肪的消化与吸收。此外,红椿菌还可作为信号分子调节肠道微生物群的组成与功能,协助维持宿主代谢平衡、肠道环境稳定及脂质代谢稳态,从而可能对血压调控产生积极影响[50]。此外,红椿菌所产生的SCFA可通过激活肠-血管轴中的SCFA-G蛋白偶联受体43/41信号通路,抑制血管平滑肌细胞内钙离子浓度,并调节肾素-血管紧张素系统。该过程能够诱导血管舒张,抑制血管紧张素转化酶活性,减少血管紧张素Ⅱ的生成,改善血管内皮功能,最终发挥降压作用[51]。
HF是多种CVD的终末期阶段,常伴随高致残率和死亡率。在本研究中,我们鉴定出8种与HF存在因果关系的细菌群。该结果与之前一项MR研究[52]的发现存在差异,可能是因为本研究采用的测量方法及样本来源不同。基于HF的肠道假说,心输出量降低和体循环改变可引起肠道灌注不足和黏膜缺血。肠道屏障受损会进一步增加肠道通透性,促进微生物移位,并使微生物代谢产物进入血液循环,最终在HF患者中诱发轻度慢性炎症[53]。本研究发现,6种肠道微生物可能通过产生TMAO对HF发挥促进作用。临床前研究[54]表明,TMAO可能通过多种机制直接损害心脏功能,包括诱导心肌肥大与纤维化、引发血管内皮损伤,并进一步触发全身性血管炎症反应,以及导致心脏线粒体功能障碍,从而加剧HF的进展。此外,有报道[55]指出TMAO可增加氧化应激并加速血管衰老,其特征是人脐静脉内皮细胞增殖和迁移能力受损。TMAO还能促进血栓形成、诱导细胞凋亡、加重心肌纤维化并降低心肌收缩力,从多途径对心血管功能产生不利影响[56]。这些发现为开发针对肠道微生物群的HF干预策略提供了重要契机,例如个性化饮食调节、益生菌和益生元的应用。同时,靶向TMAO代谢的植物化学物质也展现出潜在的治疗价值。
根据研究结果,本研究的MR分析提供了一种不依赖肠道代谢物作为媒介、直接建立肠道微生物群与疾病之间关联的方法,这在方法学上相较于以往研究具有创新性。然而,肠道微生物群组成复杂,其系统性的机制解析仍需高水平、深入的实验研究。
本研究仍存在一些局限性。第一,尽管采用MR方法有助于控制混杂因素、提高结果的准确性与可靠性,但微生物组在多个分类水平上的结构差异是一个重要的协变量,这种差异很大程度上源于人群异质性与技术方案的不同。第二,所使用的4项GWAS均采用meta分析,整合了多个队列数据;该方法虽有助于捕捉种族特异性遗传特征,但也可能因样本重叠而引入偏倚。第三,MR分析通常反映的是长期暴露效应,因此其对暴露因素的影响仍需通过随机对照试验进一步验证。第四,受原始数据所限,本研究未能进行生存分析。第五,这4项GWAS的参与者主要为非亚洲人群,因此本研究的结果主要适用于欧美人群,在推广至其他种族时应保持谨慎。
本研究揭示了肠道微生物群与CVD之间存在保护性与致病性的关联,并由此鉴定出若干肠道微生物可作为潜在的治疗靶点。总体而言,本研究结果表明,特定肠道细菌是CVD的保护因子,而另一些则作为风险因子参与CVD的进展。
作者贡献
闫奎坡负责文章研究设计;马会华负责数据收集及分析、论文撰写;刘刚、徐亚洲、张磊、孙彦琴参与论文的写作与修改。所有作者均阅读并同意最终稿件的提交。
Authors' Contributions
YAN Kuipo was responsible for the research design of the article. MA Huihua was responsible for data collection and analysis, and manuscript writing. LIU Gang, XU Yazhou, ZHANG Lei, and SUN Yanqin contributed to the writing and revision of the manuscript. All authors have read the final version of paper and consented to its submission.
利益冲突声明
所有作者声明不存在利益冲突。
Conflict of Interests
All authors declare no relevant conflict of interests.
参考文献
/
| 〈 |
|
〉 |

