色素沉着绒毛结节性滑膜炎(PVNS)是一种罕见的、侵袭性的滑膜增生性疾病,主要特征为关节腱鞘和滑囊内的含铁血黄素沉积以及绒毛结节样增生
[1],在发达国家的发病率为1.8/100万
[2]。PVNS常可发生在手、肘、髋、膝和踝等关节部位,并且明显偏向于膝关节,由于其较强的侵袭性,可累及关节内外,导致关节肿胀、积液以及进行性骨质破坏,严重时需要进行关节置换手术,却仍有一定的复发率
[3]。目前对于PVNS多种治疗方法的疗效存在争议,其已经成为临床治疗的难点,备受关注;且对于PVNS的病因及发病机制尚不完全清楚,但研究推测其与脂质代谢紊乱、创伤出血、炎症过程和肿瘤相关。
肠道菌群(GM)作为人体内重要的微生物系统,其组成以及多样性可影响多种疾病(如炎症性疾病、肥胖等)进程,与人类健康息息相关
[4]。目前随着对GM研究的不断深入,国外有学者研究GM与骨骼健康之间的复杂机制,揭示了GM通过微生物衍生代谢产物在骨代谢与骨形成中充当着重要角色
[5]。例如,肠道微生物群与滑膜、软骨和骨关节的炎症密切联系,可以增强软骨愈合、降低滑膜炎症,是骨关节炎的潜在治疗与预防方式之一
[6, 7]。尚没有研究证实GM与PVNS发病机制之间的联系。
孟德尔随机化(MR)是一种利用遗传变异作为工具变量的有效分析方法,通过对微生物组群数据和全基因组关联研究数据(GWAS)的综合分析,从而筛选与GM组成紧密相关的遗传变异,克服传统研究的局限性,降低混杂因素及反向因果所造成的偏倚,为因果推论研究提供更可靠的证据
[8]。本研究旨在揭示GM与PVNS之间的因果关系,探究其潜在机制,为PVNS的预防和治疗方案提供新思路。
1 资料和方法
1.1 研究设计
本研究设计采用两样本双向MR分析数据,为GM和PVNS之间的因果关系提供强有力的支持。利用来自公共数据的GWAS荟萃分析,选择GM为暴露因素,筛选与GM显著相关的单核苷酸多态性(SNP)作为工具变量(IVs);PVNS作为结局变量,而后进行MR分析,并通过Cochran's Q与MR Egger法进行异质性和水平多效性检验,最后进行敏感性分析以验证结果的可靠性。为把控文章质量,获取更可靠的数据,严格遵循满足MR的3个基本假设:选择与GM显著关联的SNPs作为工具变量;SNPs应是相对独立的,与已知的混杂因素(影响GM和PVNS)无关;SNPs只通过暴露 (GM或PVNS)影响结局(PVNS或GM)(
图1)。
1.2 数据来源
本研究与GM相关的遗传信息数据来源于国际MiBioGen联盟的一项研究
[9],该研究包括24个队列的18 340个个体,主要分析了参与者粪便菌群16S rRNA 测序图谱来查明影响肠道微生物群相对丰度的遗传位点,最终检查评估了211个门、纲、目、科和属。与PVNS相关的遗传信息数据来源于芬兰数据库R9发布的数据集(
https://storage.googleapis.com/finngen-public-data-r9/summary_stats/finngen_R9_M13_VILLONODULAR.gz),其中包括212例PVNS病例和240 862例作为对照。为证明本研究结果的可靠性和可重复性,我们选取了来源于芬兰数据库R10发布的最新PVNS数据集按相同条件进行验证分析。
1.3 工具变量选择
严格控制数据质量,选择合适的工具变量,以保障有关GM组成与PVNS之间因果关系推测的完整性和准确性。 (1)本研究为满足相关性假设,在GWAS数据集中筛选与GM显著相关的SNPs(
P<1×10
-5)
[10,11],以确保SNPs与工具变量显著相关。(2)设置参数k=10 000,
r2=0.001,以去除连锁不平衡的影响,避免遗传变异残留所造成的误差
[12]。(3)为使数据协调,从分析中删除所有的回文SNP。(4)通过Pheno Scanner
[13](
http://www.phenoscanner.medschl.cam.ac.uk/),我们已经对所有SNP在二级表型进行检索,以尽可能排除混杂因素的影响。(5)使用
F统计值评估IV的强度,设置
F>10作为强IV的统计阈值,否则视为IV与暴露的相关性较差,应予以排除
[14]。
1.4 统计学分析
本研究使用RStudio软件(4.3.2版)中TwoSampleMR、 MR-PRESSO和MRcML软件包进行数据分析。
1.4.1 孟德尔随机化分析
本研究使用逆方差加权分析方法(IVW)为分析GM与PVNS之间因果关系的主要方法
[15],该方法是根据方差的倒数计算加权平均值,当所有IV有效且不存在异质性及多效性时,选定IVW的结果是最可靠的。同时采用MR-Egger回归和加权中位数(WME)作为补充分析方法。MR-Egger回归也可以用来检验因果关系,但本研究结合MR-Egger回归排除一些结果
[16]。WME只有当存在至少50%有效的IV时,才能够获得较为准确的估计结果
[17]。为了更好的分析GM与PVNS之间的因果关系,本研究还采用同样的数据源和参数设置,将正向MR中发现的与结局显著相关菌群和PVNS进行反向MR分析。
1.4.2 敏感性分析
本研究进行敏感性分析以评估结果的稳健性,包括逐一剔除检验、异质性检验和水平多效性检验。采用留一法逐个剔除每一个SNP,然后重新进行IVW分析,以评估结果是否存在偏倚或单个SNP是否对结果产生特异性影响
[18]。使用Cochran's Q统计量和MR-Egger回归评估异质性,若
P<0.05则认为检验具有统计学意义,存在异质性。MR-Egger截距用来评估多效性,以确定IV不会通过混杂因素影响结局变量,若MR-Egger截距分析中
P<0.05,则认为检验具有统计学意义,具有水平多效性。同时本研究还采用了MR- PRESSO法,以分析和校准IVs中显著异常值,以降低水平多效性对因果推论的影响
[10];采用一种基于约束最大似然和模型平均的MR方法cML-MA,用于控制相关和不相关的多效性效应,有助于提供更可靠的结果
[19]。
2 结果
2.1 工具变量筛选结果
根据IVW结果显示,9种不同的肠道菌群(1目、1科和7属)可能与PVNS之间可能存在因果关联,同时使用MR-Egger回归和加WME进行补充分析。其中放线菌门(Actinobacteria)和Family XIII AD3011 group由于可能存在潜在的离群值,MR-Egger回归法与IVW法分析结果方向相悖,为确保数据的准确性故排除。最终筛选出7种(1科和6属)可能与PVNS之间可能存在因果关联。最后SNPs的筛选结果和结果如
表1所示,最终筛选出80个SNPs,且统计量
F均>10。
2.2 肠道菌群与色素沉着绒毛结节性滑膜炎的因果效应
进行孟德尔随机化分析,3种分析方法结果(
图2、
表2)。
本研究发现了个是PVNS风险因素的微生物群,包括巴恩斯氏菌属(Barnesiella)(OR=3.12,95% CI:1.15~8.41,P=0.025)和瘤胃球菌科UCG010(RuminococcaceaeUCG010)(OR=4.03,95% CI:1.19~13.68,P=0.025)。此外。我们还发现5个微生物群是PVNS的保护因素,包括毛螺菌科(Lachnospiraceae)(OR=0.33,95% CI:0.12~0.91,P=0.032)、另枝菌属(Alistipes)(OR=0.16,95% CI:0.05~0.53,P=0.003)、经黏液真杆菌属(Blautia)(OR=0.20,95% CI:0.06~0.61,P=0.005)、毛螺菌科FCS020(LachnospiraceaeFCS020group)(OR=0.38,95% CI:0.15~0.94,P=0.036)和瘤胃球菌科UCG014(RuminococcaceaeUCG014)(OR=0.36,95% CI:0.14~0.94,P=0.037)。GM与PVNS的MR分析的3种主要分析方法结果的散点图(图3)。
2.3 敏感性分析
结果显示所有的Q_pval均>0.05,表明各微生物群分类组的全部SNP均不存在异质性。用MR-Egger回归评估SNPs与结果之间的多效性,MR-Egger回归中截距项的P值始终超过0.05,表明不存在水平多效性。同时,在MR-PRESSO分析中,如果相关微生物群群落的全局检验P值大于0.05,则表明不存在水平多效性。多效性检验是MR分析最重要的环节,其确保了结果的准确性和稳健性,MR-Egger截距回归分析较为合理的解决了不相关水平多效性的偏倚,但仍无法控制相关水平多效性的影响。故本研究最后采用RStudio软件的MRcML包对GM与PVNS进行分析,结果显示GM与PVNS结果不存在相关水平多效性(表3)。留一法检验结果证实了本研究MR结果的稳健性,未见单个SNP对研究结果产生较大偏倚(图4)。
2.4 反向MR分析的结果
采用IVW方法进行反向MR分析。结果显示在相同的筛选条件下(P<1×10-5,k=10 000,r2=0.001),未鉴定出与PVNS密切相关的工具变量。
2.5 可重复性验证
选择芬兰数据库最新的PVNS数据集再次分析,结果显示,毛螺菌科(OR=0.32,95%
CI:0.12~0.85,
P=0.022)、另枝菌属(OR=0.25,95%
CI:0.08~0.76,
P=0.014)、经黏液真杆菌属(OR=0.19,95%
CI:0.07~0.54,
P=0.002)是PVNS的保护因素;瘤胃球菌科UCG010(OR=6.57,95%
CI:2.03~21.67,
P=0.002)是风险因素(
表4)。
3 讨论
本研究旨在探究GM 与PVNS之间的因果关系,通过MR研究全面评估211种肠道微生物,鉴别出9种GM与PVNS相关,对研究结果进行严格的质量控制,避免混杂因素及反向因果所造成的影响,最终显示与PVNS患病风险存在因果关系的GM有7种。其中,增加PVNS患病风险的菌群有2种,包括巴恩斯氏菌属和瘤胃球菌科UCG010;有5种菌群对PVNS的发病存在保护作用,包括毛螺菌科、另枝菌属、经黏液真杆菌属、毛螺菌科FCS020和瘤胃球菌科UCG014。
目前尚未有研究能完全阐述GM与PVNS之间的作用机制,但仍有现存的临床研究表明两者存在一定的联系,这可以为本研究的结果提供一定的证据支撑。巴恩斯氏菌属归属于拟杆菌门,是一种近些年新发现的菌属,大多研究认为其作为一类肠道益生菌,可以帮助维持肠道菌群的平衡,促进消化和吸收营养物质,它还可以抑制有害菌的生长,提高免疫系统的功能
[20-22]。有研究发现其可以作为一种“抗肿瘤益生菌”存在,通过促进特异性T细胞的浸润,改善抗肿瘤调节剂的疗效
[23]。然而,另外有研究表明巴恩斯氏菌属可能与癌症的发生有关。国内有研究表明,对正常小鼠和原发性肝癌小鼠肠道菌群进行多样性研究,二者间巴恩斯氏菌属存在显著性差异,推测其与肝癌发生存在可能的相关性
[24]。而恰恰有“肿瘤源”学说认为PVNS可能属于一种肿瘤样病变,可广泛侵袭滑膜组织,并且恶变的PVNS与肿瘤疾病的一般特性更加相似
[25, 26],这表明巴恩斯氏菌属在PVNS中可能也具有尚未发现的影响,增加其潜在的发生风险。另枝菌属也属于拟杆菌门,是一类常见的肠道细菌,有研究证实其被视为潜在的短链脂肪酸(SCFAs)生产者而具有潜在的抗炎作用,可以在癌症免疫治疗中通过调节肿瘤微环境发挥有益作用
[27]。Iida等
[28]进行的一项动物研究认为,另枝菌属与Toll样受体4(TLR4)启动和肿瘤坏死因子(TNF)产生的作用之间存在正相关联系,可与TLR4相结合从而诱导肿瘤相关髓样细胞启动TNF表达,从而改变微环境,使肿瘤发生细胞凋亡。同时结合我们的研究结果,我们推测另枝菌属可通过调节免疫应答从而抑制滑膜肿瘤样病变,成为PVNS的保护因素。但是目前现存研究并未有明确的结论,要明确两者之间的因果关系还需要更多的研究。在我们的研究中显示巴恩斯氏菌属对罹患PVNS起促进作用,另枝菌属则降低PVNS的发生风险,这可能为我们以后阐述PVNS的发病机制及更加有效的治疗方案提供新思路。
瘤胃球菌科UCG010和瘤胃球菌科UCG014同属于瘤胃球菌科,作为一种革兰阳性厌氧菌,可以通过酵解宿主消化系统中的纤维素、葡萄糖以及木糖醇等获取能量,在人体新陈代谢中发挥着至关重要的作用。本研究结果显示两个菌属对于PVNS的发生产生不同的影响,这可能是因为瘤胃球菌科既包括有益菌又包括有害菌。目前尚未有研究明确瘤胃球菌科与PVNS之间的因果关系。但有研究显示,瘤胃球菌科UCG010可与血脂异常呈正相关,并且在肥胖易感的小鼠中,发现富集的瘤胃球菌科可以促进脂肪的合成
[29, 30]。如今有学者认为脂质代谢紊乱可能是PVNS的病因之一
[31],所以我们推测肠道中异常的瘤胃球菌科UCG010积累,可导致人体脂肪代谢的紊乱,从而增加发生PVNS的风险。另外有研究显示,瘤胃球菌科能够产生SCFAs,尤其是丁酸盐,是肠上皮细胞的主要能量来源,同时抑制促炎细胞因子的信号通路,具有抗炎作用,在调节肠道炎症方面具有重要作用
[32];同时,也有研究发现瘤胃球菌科与关节炎类疾病,如膝骨关节炎等有一定的相关性,在一项观察艾灸治疗膝骨关节炎大鼠肠道菌群及炎症影响的研究中,发现艾灸可以明显调节体内炎症因子水平,降低促炎因子白细胞介素-1β(IL-1β)和肿瘤坏死因子-α(TNF-α),同时相应地发现瘤胃球菌科UCG014丰度增加,证实瘤胃球菌科UCG014与炎症因子水平的相关性
[33]。有国内学者探究PVNS的发病机制,发现PVNS患者滑液中的炎症因子(IL-1β和TNF-α)显著上调和钙粘蛋白-11的高表达,肯定了炎症因子在PVNS发病中的重要影响
[34]。所以我们认为瘤胃球菌科UCG014可能通过降低炎症因子表达来减小PVNS的易感风险。
经黏液真杆菌属和毛螺菌科FCS020同属于毛螺菌科,是人体最丰富的肠道菌群之一,存在于大多数健康人群的肠道中,同样可以参与人体碳水化合物的代谢,水解淀粉以及其他多糖产生SCFAs和丁酸盐等
[35,36]。根据现有研究显示
[37],丁酸盐能够显著降低软骨炎症因子的水平,包括 IL-1β、基质金属蛋白酶-13 (MMP-13)和 TNF-α 等,还能够降低受SCFAs调节的iNOS水平,抑制氧化应激反应,并恢复紧密连接蛋白、改善肠道组织形态,抑制肠道损伤,通过减少炎症细胞死亡和改善自噬过程保护软骨细胞。这在维持骨与软骨的完整性和防御机制方面起着至关重要的作用。若肠道菌群紊乱,肠道炎症发生,部分肠道免疫炎症因子如IL-1β、MMP13 和 TNF-α 等可以通过体液循环直接从肠道转移到骨关节,从而影响骨系统的免疫活性。此外,肠道微生物组对免疫细胞的成长和分化有一定影响,肠道细菌及其代谢产物可以与免疫细胞相互作用。一项动物研究显示
[38],丁酸可以调节调节性T细胞(Tregs)的数量和功能,并减少辅助性T细胞17(Th17)的数量且不促进其极化。合理的肠道菌群组成有助于提升对骨关节疾病产生有效免疫反应的能力并维持全身免疫稳态。根据本研究结果,我们可以深入探究毛螺菌科降低PVNS潜在发生风险的保护机制,探索新的干预或治疗措施如增加益生菌或粪便微生物群移植等,这可能成为一种新的有效预防或治疗PVNS的临床策略。
本研究仍存在一定的局限性。首先,MR基于作为工具变量的遗传变异仅通过研究的暴露因素影响结果这个假设。虽然尽可能地剔除混淆因素影响,并进行异质性和敏感性分析来进一步确保研究结果的稳定性,但仍然可能存在未纳入的混淆因素。其次,本研究数据集主要来自于欧洲人种,对亚洲人种可能不具有广适性。对于亚洲人群而言,可能存在未知的肠道菌群对色素沉着绒毛结节性滑膜炎存在影响。此外,本研究为了确保结果的准确性和外推性,拟采用多队列GWAS数据对本研究结果进行验证,但现在鲜针对有PVNS的全基因组关联测序研究,故本研究验证队列数据来自于芬兰R10数据库,但两版本的结局变量数据集的遗传相关性研究在统计学准许范围内,样本重叠对本研究主要结果不存在较大偏倚影响。
综上所述,本研究通过MR探讨GM与PVNS之间的因果关系,发现7种特定肠道微生物的丰度变化与PVNS的发病风险增加或降低相关,可能有助于阐明GM对炎症因子和免疫调节作用的影响,并激发对PVNS预防和治疗领域医学的发展。