A Mendelian randomization study on the causal association between social stress and tinnitus onset
-
摘要:目的
基于两样本孟德尔随机化方法探讨社会压力与耳鸣发病之间的因果关系。
方法从全基因组关联研究(GWAS)数据库中分别获取社会压力、耳鸣相关遗传数据, 筛选出相互独立且与社会压力高度相关的单核苷酸多态性(SNP)作为工具变量，以耳鸣作为结局变量。分别采用逆方差加权法(IVW)、加权中位数法和MR-Egger回归法进行孟德尔随机化分析。采用MR-Egger回归法的截距项评估水平多效性，采用Cochran's Q统计量评估异质性，并采用留一法进行敏感性分析。
结果社会压力数据包含459 742例样本，耳鸣数据包含117 882例样本，共筛选出10个与社会压力密切相关的SNP作为工具变量。随机效应IVW分析结果显示，社会压力与耳鸣发病存在因果效应(OR=1.251, P=0.036)。随机效应IVW、加权中位数法、MR-Egger回归法的分析结果相似(OR=1.251、1.274、1.438), 表明社会压力为耳鸣的危险因素，两者之间存在正向因果效应，且结果不存在异质性及多效性。
结论社会压力与耳鸣发病存在因果关系，或可为未来耳鸣的临床防治提供新思路。
Abstract:ObjectiveTo explore the causal relationship between social stress and tinnitus onset using a two-sample Mendelian randomization approach.
MethodsGenetic data pertaining to social stress and tinnitus were extracted from genome-wide association studies (GWAS) databases. Single nucleotide polymorphisms (SNP) that were independent and strongly correlated with social stress were selected as instrumental variables, with tinnitus serving as the outcome variable. Mendelian randomization analyses were performed using inverse-variance weighted (IVW) method, weighted median method, and Mendelian randomization-Egger regression. The intercept term from Mendelian randomization-Egger regression was utilized to assess horizontal pleiotropy, Cochran's Q statistic was employed to evaluate heterogeneity, and leave-one-out analysis was conducted for sensitivity assessment.
ResultsThe social stress dataset encompassed 459, 742 samples, while the tinnitus dataset comprised 117, 882 samples. A total of 10 SNPs tightly associated with social stress were identified as instrumental variables. The analysis results of random-effects inverse variance weighting (IVW), weighted median method, and Mendelian randomization-Egger regression were similar (OR=1.251, 1.274, 1.438), indicating that social stress was a risk factor for tinnitus, and there was a positive causal effect between them, with no heterogeneity or pleiotropy in the results.
ConclusionThere is a causal relationship between social stress and tinnitus onset, which may offer novel perspectives for the clinical prevention and treatment of tinnitus in the future.
-
社会压力是指社会环境或外部事件对个体精神层面造成的影响，是内在心理反应与外部环境相互作用的结果[1]。过高的社会压力难以缓解，可引起内分泌、免疫系统及心理健康等方面失衡，进而会对神经系统造成损害，导致失眠、焦虑、抑郁、耳鸣等神经系统疾病高发，因此社会压力与神经系统疾病间确切因果关系的研究日益受到重视[2-3]。耳鸣为神经系统常见疾病之一，是指在无外界声源刺激时感知到耳内有声音，具有高发病率和高复发率的特点，严重影响患者的生活质量[4]。耳鸣患者常伴有抑郁、敏感、紧张等情绪障碍，这些负面情绪可能会进一步加重耳鸣症状，但两者间的因果关系难以明确。既往研究[5]表明，社会压力是耳鸣发病的危险因素之一，且与耳鸣的严重程度呈正相关。然而，由于传统的观察性研究通常样本量有限，且易受到吸烟、熬夜等混杂因素的干扰，社会压力与耳鸣之间的确切关联性及因果关系目前仍存在争议。孟德尔随机化作为一种基于孟德尔独立分配定律的流行病学研究设计和数据分析方法，已被广泛用于评估暴露因素与疾病结局之间的因果效应[6-7]。孟德尔随机化利用与社会压力强相关的遗传变异，即单核苷酸多态性(SNP)作为工具变量，可以减少环境因素、疾病进程及其他混杂因素对遗传变异的影响，进而增强对因果关系的推断能力[8-9]。本研究采用孟德尔随机化方法评估社会压力与耳鸣发病之间的关联性和因果关系，并通过敏感性分析检验结果的稳健性，现报告如下。
1. 资料与方法
1.1 研究设计
本研究是一项基于全基因组关联研究(GWAS)数据库的公开数据进行的两样本孟德尔随机化研究，其中暴露因素为社会压力，结局变量为耳鸣。本研究中，孟德尔随机化分析要求工具变量需满足以下3个核心假设: ①相关性假设，即工具变量(SNP)与暴露因素(社会压力)具有强相关性。在全基因组水平上，选择与社会压力显著相关的SNP(相关性的阈值为P＜5×10-8), 并设置参数r2＜0.001, 遗传距离>10 000 kb以排除连锁不平衡的干扰。若符合条件的工具变量过少，则将相关性的阈值放宽至P＜5×10-6。工具变量的相关性强度通过计算统计量F进行衡量[10], 统计量F>10表示工具变量与暴露因素具有强相关性。计算公式为
$$ \begin{aligned} & F=\frac{R^2 \times(N-k-1)}{k \times\left(1-R^2\right)}, \\ & R^2=\frac{2 \times \mathrm{EAF} \times(1-\mathrm{EAF}) \times \beta^2}{2 \times \mathrm{EAF} \times(1-\mathrm{EAF}) \times \beta^2+2 \times \mathrm{EAF} \times(1-\mathrm{EAF}) \times N \times \mathrm{SE}^2}, \end{aligned} $$ 其中N为暴露因素的样本量, k为纳入的工具变量数目, R2为工具变量解释的暴露因素变异比例, EAF为效应等位基因频率, β为效应值, SE为β值的标准误差。②排他性假设，即工具变量(SNP)与已知的混杂因素不相关，仅通过暴露因素(社会压力)影响结局(耳鸣)。若工具变量通过其他途径对结局变量产生影响，则表示相关遗传变异存在多效性，违反排他性假设。本研究通过多种敏感性分析方法检测多效性，以确保结果的稳健性, P＞0.05表示结果不存在多效性。③独立性假设，即工具变量(SNP)与任何潜在影响暴露及结局的混杂因素无关，仅通过社会压力影响耳鸣，不与耳鸣直接相关。本研究通过Phenoscanner数据库查询并剔除与混杂因素相关的SNP[11]。剔除的混杂因素包括内分泌系统疾病(甲状腺功能亢进、甲状腺功能减退等)、消化系统疾病(胃溃疡、慢性胃炎等)、生活习惯(吸烟、饮酒等)以及与暴露相关的其他神经系统疾病(焦虑、抑郁、失眠、头痛等)。具体流程见图 1。
1.2 研究数据
从GWAS数据库(https://gwas.mrcieu.ac.uk)中分别获取社会压力、耳鸣相关的遗传数据，所有遗传数据均为合法公开的研究资料。暴露因素(社会压力)的数据(ukb-b-9667)来源于英国生物银行(UKB)数据库，该数据集共纳入459 742名欧洲人的信息，包含9 851 867个SNP; 结局变量(耳鸣)的数据(ukb-d-4803_14)亦来源于UKB数据库，该数据集共纳入117 882名欧洲人的信息，包含13 567 874个SNP。为了减少潜在的种族偏倚，本研究仅纳入欧洲人群的数据，且性别不限，数据的基本信息见表 1。
表 1 暴露因素与结局数据的基本信息暴露因素/结局 数据来源 样本量/例 单核苷酸多态性数量/个 样本种族来源 性别 数据公布年份 社会压力 ukb-b-9667 459 742 9 851 867 欧洲 男/女 2018 耳鸣 ukb-d-4803_14 117 882 13 567 874 欧洲 男/女 2018 1.3 工具变量筛选
以社会压力为暴露因素，应用R studio软件对社会压力的GWAS数据进行提取并筛选，共筛选出10个与社会压力相关的工具变量(SNP), 见表 2。筛选过程中，确保所选SNP均与表型密切相关(P＜5×10-8), 且排除连锁不平衡(r2＜0.001, 遗传距离>10 000 kb)。经计算，各SNP的统计量F均大于10, 有效避免了弱工具变量偏倚的发生。将获得的与社会压力显著相关的SNP作为本研究的工具变量，见表 3。
表 2 3种孟德尔随机化分析方法的结果方法 单核苷酸多态性数量/个 β SE P OR MR-Egger回归法 10 0.363 0.435 0.428 1.438 加权中位数法 10 0.242 0.137 0.078 1.274 逆方差加权法 10 0.224 0.107 0.036 1.251 表 3 与社会压力相关工具变量的特征单核苷酸多态性 基因位置 染色体 等位基因 SE β P EAF F EA OA rs72668117 3722721 1 A G 0.002 32 0.013 51 6.00E-09 0.021 06 33.832 87 rs359234 60471409 2 G C 0.000 69 0.004 43 1.40E-10 0.620 15 41.145 87 rs10935037 132650822 3 C A 0.000 67 -0.003 81 1.30E-08 0.510 09 32.383 25 rs73128428 39308367 7 T C 0.000 93 -0.005 29 1.30E-08 0.152 33 32.334 06 rs34883191 47261533 14 T G 0.000 67 -0.003 66 4.70E-08 0.521 63 29.827 52 rs4906365 104229230 14 C G 0.000 68 0.003 95 6.00E-09 0.419 09 33.821 93 rs6496749 91552003 15 T C 0.000 86 0.004 77 2.80E-08 0.186 64 30.813 01 rs2884737 31105554 16 C A 0.000 76 -0.004 25 2.70E-08 0.255 00 30.939 21 rs12923462 10328558 16 G A 0.000 70 0.003 97 1.60E-08 0.641 75 31.975 96 rs784220 53349931 18 C G 0.000 95 0.005 51 7.10E-09 0.855 89 33.499 32 1.4 统计学分析
本研究使用R 4.3.1软件及TwoSampleMR软件包(0.5.7)进行孟德尔随机化分析，采用逆方差加权法(IVW)、MR-Egger回归法、加权中位数法(WM)等方法检验社会压力与耳鸣之间是否存在因果关系，并以IVW的结果为主(P＜0.05表示具有统计学意义)。该分析拟合时使用的权重为结局方差的倒数，回归时不考虑截距项。基因的多效性通过MR-Egger回归法的截距项进行评估, P＜0.05表示存在水平多效性。通过Cochran′s Q统计量评估异质性, P＜0.05表示存在异质性。采用留一法(即依次去除每个SNP, 计算剩余SNP的合并效应)进行敏感性分析，评价单个SNP是否对暴露因素与结局的关系产生影响[12-14]。
2. 结果
2.1 工具变量的筛选及弱工具变量偏倚的判断
本研究从GWAS数据中最终筛选出10个与社会压力显著相关的SNP作为工具变量，所选SNP均具有显著相关性且相互独立(P＜5×10-8, r2＜0.001)。此外，计算结果显示所选工具变量与社会压力存在强相关关系(F>10), 表明弱工具变量不会导致结果发生偏倚，结果具有可靠性。
2.2 两样本孟德尔随机化分析结果
两样本孟德尔随机化分析结果显示，随机效应IVW得到的结果为β=0.224, P＜0.05, 见图 2, 表明社会压力是耳鸣的危险因素(OR=1.251), 社会压力的增加可升高耳鸣的发病风险。3种孟德尔随机化分析方法得到的结果相似，见表 2。各种孟德尔随机化分析方法得到的因果效应方向一致，见图 3。
2.3 敏感性分析
2.3.1 水平多效性检验
采用MR-Egger回归法的截距项对水平多效性进行评估，结果显示P=0.750(P＞0.05), 表明所选SNP间不存在水平多效性，即所选工具变量具有可靠性，见表 4。
表 4 社会压力与耳鸣的敏感性分析结果方法 异质性检验 水平多效性检验 Cochran′s Q P 截距 P 逆方差加权法 6.666 0.672 — — MR-Egger回归法 6.557 0.585 -0.001 0.750 2.3.2 异质性检验
Cochran′s Q检验结果显示, P=0.585(P＞0.05), 未发现异质性较高的SNP, 见表 4。漏斗图分析结果显示，所选工具变量的因果效应分布基本对称，未发现偏倚存在，见图 4。
2.3.3 留一法敏感性分析
采用留一法敏感性分析方法，逐一剔除所选SNP, 并分析其是否影响总体因果效应，结果显示，并未出现显著影响总体因果效应的SNP, 见图 5。
![]()
图 5 留一法敏感性分析结果黑点表示剔除此SNP后估计的社会压力对耳鸣的因果效应，黑实线表示相应的95%CI; 红点表示使用筛选过的SNP估计社会压力对耳鸣发病风险的整体因果效应，红线表示合并估计相应的95%CI。3. 讨论
耳鸣是神经系统常见疾病之一，其发病已被证实与社会压力有密切联系。潘芳芳等[5]的横断面研究中纳入183例耳鸣患者，压力评估结果显示耳鸣的严重性与压力呈正相关。BETZ L T等[15]对41例患者的研究结果表明，在压力增大的时期，耳鸣的严重程度随之增加。瑞典统计局开展的一项纳入12 166人的大规模流行病学研究[16]发现，社会心理压力大的个体发生耳鸣的概率与暴露于职业噪声环境中的个体相似，且社会心理压力可导致耳鸣症状恶化。以上研究提示，耳鸣与社会压力之间可能存在因果关联。然而，传统的观察性研究易受混杂因素的影响，因此社会压力与耳鸣之间的因果关系证据强度较弱。本研究采用两样本孟德尔随机化方法，从遗传学角度评估社会压力与耳鸣之间的因果关系，结果显示，在欧洲人群中，社会压力与耳鸣的发生发展风险呈正相关，且社会压力是耳鸣的危险因素之一，与既往研究结果相符。
目前，社会压力与耳鸣之间的机制研究主要集中在中枢神经可塑性方面。社会压力是指个体在社会环境中受到外界应激源刺激后，经大脑感知并进行认知评价后出现的生理及心理紧张应对状态[17-18]。现代社会中，人们普遍面临着较大的社会压力，若过多的压力无法缓解，可对神经系统造成损害，从而增加神经系统疾病发生风险。ZHANG G W等[19]予雄鼠社会压力源后发现，其内侧视前区谷氨酸能神经元钙活性显著增加。LIANG Y R等[20]在小鼠童年期给予GluN2B抑制剂注射，并在成年期给予社会压力刺激，通过蛋白质印迹法检测发现，前额叶皮层GluN2B蛋白表达量不减反增。由此可见，社会压力与谷氨酸信号代谢通路密切相关。谷氨酸作为中枢神经系统最常见的兴奋性神经递质，当社会压力增大时，会刺激神经元突触前膜释放谷氨酸，进而持续激活N-甲基-D-天冬氨酸受体(NMDAR), 导致大量钙离子内流，引发细胞内钙离子超载，破坏钙离子稳态，激活一系列酶和途径，引起神经突触强度变化，诱导突触可塑性[21-22]。随着研究的深入，越来越多的学者[20]开始关注耳鸣的中枢机制，并认为耳鸣不仅是耳鼻喉科疾病，也是神经系统的常见疾病。耳鸣的发生被认为与谷氨酸在听觉通路各层面的变化密切相关，既往研究[23-24]利用高剂量水杨酸盐诱导耳鸣时，通常伴随着NMDAR水平的升高。因此，社会压力与耳鸣在神经中枢分子水平上存在关联，且谷氨酸及NMDAR可能是两者关联的关键因素。然而，社会压力与耳鸣还可能受其他系统的影响，需要进一步深入研究[25]。
孟德尔随机化具有因果相关性显著、潜在混杂因素影响小、不受反向因果关系影响以及研究设计类似随机分配的特点，使得本研究的结果能够避免混杂因素或其他偏倚的影响，表明社会压力与耳鸣之间存在关联性。此外，本研究采用孟德尔随机化分析的3种统计方法，预测结果显示社会压力与耳鸣发生风险可能呈正相关，因此研究结果较为可靠。然而，本研究也存在一些局限性: 首先，两样本均来自UKB数据库，可能存在样本重叠的问题，导致结果出现偏倚; 其次，“暴露-结局”数据仅涉及欧洲人群，可能受地域和文化差异的影响，限制研究结果向其他人群推广; 最后，本研究未进一步分析反向因果关系。
综上所述，孟德尔随机化分析结果显示，社会压力与耳鸣的发生存在因果关系，且两者之间的遗传预测呈正相关。明确社会压力与耳鸣之间的关系，或可为未来耳鸣的临床防治提供新思路，即临床医师可通过调节社会压力水平对耳鸣患者进行干预。
-
![]()
图 5 留一法敏感性分析结果
黑点表示剔除此SNP后估计的社会压力对耳鸣的因果效应，黑实线表示相应的95%CI; 红点表示使用筛选过的SNP估计社会压力对耳鸣发病风险的整体因果效应，红线表示合并估计相应的95%CI。
表 1 暴露因素与结局数据的基本信息
暴露因素/结局 数据来源 样本量/例 单核苷酸多态性数量/个 样本种族来源 性别 数据公布年份 社会压力 ukb-b-9667 459 742 9 851 867 欧洲 男/女 2018 耳鸣 ukb-d-4803_14 117 882 13 567 874 欧洲 男/女 2018 
下载: 导出CSV
表 2 3种孟德尔随机化分析方法的结果
方法 单核苷酸多态性数量/个 β SE P OR MR-Egger回归法 10 0.363 0.435 0.428 1.438 加权中位数法 10 0.242 0.137 0.078 1.274 逆方差加权法 10 0.224 0.107 0.036 1.251
下载: 导出CSV
表 3 与社会压力相关工具变量的特征
单核苷酸多态性 基因位置 染色体 等位基因 SE β P EAF F EA OA rs72668117 3722721 1 A G 0.002 32 0.013 51 6.00E-09 0.021 06 33.832 87 rs359234 60471409 2 G C 0.000 69 0.004 43 1.40E-10 0.620 15 41.145 87 rs10935037 132650822 3 C A 0.000 67 -0.003 81 1.30E-08 0.510 09 32.383 25 rs73128428 39308367 7 T C 0.000 93 -0.005 29 1.30E-08 0.152 33 32.334 06 rs34883191 47261533 14 T G 0.000 67 -0.003 66 4.70E-08 0.521 63 29.827 52 rs4906365 104229230 14 C G 0.000 68 0.003 95 6.00E-09 0.419 09 33.821 93 rs6496749 91552003 15 T C 0.000 86 0.004 77 2.80E-08 0.186 64 30.813 01 rs2884737 31105554 16 C A 0.000 76 -0.004 25 2.70E-08 0.255 00 30.939 21 rs12923462 10328558 16 G A 0.000 70 0.003 97 1.60E-08 0.641 75 31.975 96 rs784220 53349931 18 C G 0.000 95 0.005 51 7.10E-09 0.855 89 33.499 32
下载: 导出CSV
表 4 社会压力与耳鸣的敏感性分析结果
方法 异质性检验 水平多效性检验 Cochran′s Q P 截距 P 逆方差加权法 6.666 0.672 — — MR-Egger回归法 6.557 0.585 -0.001 0.750
下载: 导出CSV
-
[1] 李菁林, 潘孝富, 毕文芬. 家庭收入何以影响个体社会压力——社会规范感知的中介作用[J]. 四川师范大学学报: 社会科学版, 2023, 50(2): 124-132. [2] 许骁, 江洪. 压力应激与自主神经的联系及在心血管疾病中的作用[J]. 中国循证心血管医学杂志, 2022, 14(11): 1406-1408. [3] 李高翔, 陈昱竹, 王贵平, 等. 长期精神压力与髓系来源抑制性细胞MDSC的肿瘤免疫抑制功能[J]. 西藏科技, 2020(5): 62-66. [4] 罗伊雯, 肖永涛, 高敏倩, 等. 慢性主观性耳鸣的声治疗研究进展[J]. 山东大学耳鼻喉眼学报, 2023, 37(4): 13-23. [5] 潘芳芳, 张恒欣, 陈君, 等. 压力诱发复杂性耳鸣的感知[J]. 温州医科大学学报, 2016, 46(7): 513-516. [6] 杨玉洁, 莫念, 郑旭威, 等. 精神心理障碍与银屑病发病风险的孟德尔随机化研究[J]. 现代预防医学, 2024, 51(3): 401-406, 411. [7] HEMANI G, BOWDEN J, DAVEY SMITH G. Evaluating the potential role of pleiotropy in Mendelian randomization studies[J]. Hum Mol Genet, 2018, 27(R2): R195-R208.
[8] GEORGAKIS M K, GILL D. Mendelian randomization studies in stroke: exploration of risk factors and drug targets with human genetic data[J]. Stroke, 2021, 52(9): 2992-3003.
[9] 刘波, 游俊杰, 郑思琳, 等. 基于孟德尔随机化分析胃食管反流病与慢性阻塞性肺疾病的因果关系研究[J]. 实用临床医药杂志, 2024, 28(1): 113-117. doi: 10.7619/jcmp.20232575 [10] CARTER A R, SANDERSON E, HAMMERTON G, et al. Mendelian randomisation for mediation analysis: current methods and challenges for implementation[J]. Eur J Epidemiol, 2021, 36(5): 465-478.
[11] KAMAT M A, BLACKSHAW J A, YOUNG R, et al. PhenoScanner V2: an expanded tool for searching human genotype-phenotype associations[J]. Bioinformatics, 2019, 35(22): 4851-4853.
[12] 张彬, 滕鑫, 吴彤, 等. 肠道微生物和前庭功能障碍之间的因果关系: 孟德尔随机化研究[J]. 中风与神经疾病杂志, 2023, 40(11): 983-988. [13] SUN H L, GAO X, QUE X M, et al. The causal relationships of device-measured physical activity with bipolar disorder and schizophrenia in adults: a 2-Sample Mendelian randomization study[J]. J Affect Disord, 2020, 263: 598-604.
[14] 罗斌玉, 徐仕晗, 李逸雯, 等. 阻塞性睡眠呼吸暂停与心血管疾病的因果关系: 一项两样本孟德尔随机化研究[J]. 中国循证医学杂志, 2023, 23(10): 1137-1142. [15] BETZ L T, MÜHLBERGER A, LANGGUTH B, et al. Stress reactivity in chronic tinnitus[J]. Sci Rep, 2017, 7: 41521.
[16] BAIGI A, ODEN A, ALMLID-LARSEN V, et al. Tinnitus in the general population with a focus on noise and stress: a public health study[J]. Ear Hear, 2011, 32(6): 787-789.
[17] SCHULKIN J, STERLING P. Allostasis: a brain-centered, predictive mode of physiological regulation[J]. Trends Neurosci, 2019, 42(10): 740-752.
[18] GUPTA S, BHARATHA A, COHALL D, et al. Aerobic exercise and endocannabinoids: a narrative review of stress regulation and brain reward systems[J]. Cureus, 2024, 16(3): e55468.
[19] ZHANG G W, SHEN L, TAO C, et al. Medial preoptic area antagonistically mediates stress-induced anxiety and parental behavior[J]. Nat Neurosci, 2021, 24(4): 516-528.
[20] LIANG Y R, ZHANG X H. Inhibition of GluN2B-containing NMDA receptors in early life combined with social stress in adulthood leads to alterations in prefrontal PNNs in mice[J]. Sheng Li Xue Bao, 2024, 76(1): 1-11.
[21] MARWAHA S, PALMER E, SUPPES T, et al. Novel and emerging treatments for major depression[J]. Lancet, 2023, 401(10371): 141-153.
[22] 毕翻, 姜俊杰, 潘恒恒, 等. N-甲基-D-天冬氨酸受体参与癫痫发病机制的研究进展[J]. 癫痫与神经电生理学杂志, 2023, 32(1): 37-41. [23] 易彬, 马青坪, 汪雪玲, 等. 水杨酸盐致大鼠耳鸣相关神经通路递质变化的初步研究[J]. 听力学及言语疾病杂志, 2023, 31(4): 1-7. [24] WALLACE M N, PALMER A R. Neural plasticity in tinnitus mechanisms[J]. Brain Sci, 2023, 13(12): 1615.
[25] SINGH A, SMITH P F, ZHENG Y W. Targeting the limbic system: insights into its involvement in tinnitus[J]. Int J Mol Sci, 2023, 24(12): 9889.
计量
- 文章访问数: 114
- HTML全文浏览量: 25
- PDF下载量: 11
苏公网安备 32100302010246号