A research on cerebral limbic-cortical circuit in patients with chronic obstructive pulmonary disease complicated with depression
-
摘要:目的
应用基于体素形态学分析(VBM)探讨慢性阻塞性肺疾病(COPD)合并抑郁患者的皮质边缘环路脑灰质体积改变。
方法前瞻性选取22例COPD合并抑郁患者为COPD合并抑郁组，另选取性别、年龄相匹配的20例健康受试者为对照组。实施T2-FLAIR、高分辨3D-T1WI检查，采用慢性阻塞性肺疾病评估测试(CAT)及SPM12软件分析2组脑灰质体积变化。
结果VBM结果显示, 与对照组相比, COPD合并抑郁患者左侧额叶、双侧中央旁小叶、双侧前扣带回、左中部扣带回、左后扣带回、左背侧丘脑、下丘脑、双侧杏仁体、双侧小脑灰质体积减少，差异有统计学意义(P < 0.05)。
结论COPD合并抑郁患者存在皮质边缘环路的多个脑区灰质体积减少情况，可为临床预防、治疗及干预提供参考。
Abstract:ObjectiveTo explore the change of gray matter volume in limbic-cortical circuit in patients with chronic obstructive pulmonary disease (COPD)complicated with depression by using voxel-based morphometry (VBM).
MethodsA total of 22 patients with COPD complicated with depression were prospectively selected as COPD complicated with depression group, and 20 healthy subjects matched by gender and age were selected as control group. T2-FLAIR and high-resolution 3d-T1WI were performed, and the change of cerebral gray matter volume in both groups was analyzed by COPDAssessment Test (CAT) and SPM12 software.
ResultsThe results of VBM showed that the volumes of gray matter in left frontal lobe, bilateral paracentral lobules, bilateral anterior cingulate gyrus, left central and left posterior cingulate gyrus, left dorsal thalamus, hypothalamus, bilateral amygdala and bilateral cerebellum significantly decreased in patients with COPD complicated with depression compared with the control group (P < 0.05).
ConclusionPatients with COPD complicated with depression present gray matter volume reduction in multiple regions of limbic-cortical circuitry, which can provide reference for clinical prevention, treatment and intervention.
-
社会压力是指社会环境或外部事件对个体精神层面造成的影响，是内在心理反应与外部环境相互作用的结果[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数据库，可能存在样本重叠的问题，导致结果出现偏倚; 其次，“暴露-结局”数据仅涉及欧洲人群，可能受地域和文化差异的影响，限制研究结果向其他人群推广; 最后，本研究未进一步分析反向因果关系。
综上所述，孟德尔随机化分析结果显示，社会压力与耳鸣的发生存在因果关系，且两者之间的遗传预测呈正相关。明确社会压力与耳鸣之间的关系，或可为未来耳鸣的临床防治提供新思路，即临床医师可通过调节社会压力水平对耳鸣患者进行干预。
-
![]()
图 1 COPD合并抑郁组与对照组脑灰质差异脑区
A: 下丘脑、双侧杏仁体、双侧小脑皮质、左背侧丘脑、左后扣带回灰质减少; B: 左额叶皮质、左侧中部扣带回、双侧前扣带回及双侧中央旁小叶灰质减少。蓝色表示T值有差异的灰质减低脑区。
表 1 2组临床资料比较(x±s)[n(%)][M(P25, P75)]
临床资料 COPD合并抑郁组(n=22) 对照组(n=20) 年龄/岁 66.14±6.78 63.05±6.32 抑郁评分/分 10.00(8.75, 16.50)* 5.00(3.00, 5.75) 高血压 12(54.55) 9(45.00) 与对照组比较, * P < 0.05。 
下载: 导出CSV
表 2 2组各脑区灰质体积变化比较
脑区 体积变化百分比/% T MNI坐标 X Y Z 下丘脑 3.73 7.62 -3 3 -18 右小脑 20.80 5.73 33 -63 -18 左小脑 20.33 4.41 -39 -57 -21 左背侧丘脑 4.49 5.32 -18 -30 -3 左杏仁体 3.32 5.40 -21 3 -21 右杏仁体 2.92 4.83 24 3 -21 左额上回 6.45 5.49 -12 24 63 左额中下回 14.97 6.20 -48 -18 39 左前中扣带回 4.43 3.95 -9 15 36 右前扣带回 3.75 3.92 6 24 36 左后扣带回 6.59 4.32 -12 -63 12 左中央旁小叶 5.32 3.92 -6 -21 66 右中央旁小叶 2.89 3.87 3 -18 66 MNI: 蒙特利尔神经病学研究所。
下载: 导出CSV
-
[1] 李思其, 高兴林. 《慢性阻塞性肺疾病全球倡议》2021年版更新解读[J]. 临床药物治疗杂志, 2021, 19(5): 36-42. doi: 10.3969/j.issn.1672-3384.2021.05.008 [2] 黄菲菲, 秦继山. 中老年慢性阻塞性肺疾病临床特征与病情进展的相关性[J]. 实用临床医药杂志, 2021, 25(1): 45-49, 53. https://www.cnki.com.cn/Article/CJFDTOTAL-XYZL202101012.htm [3] DENGX Y, FUJ, SONGY, et al. Glucocorticoid receptor dysfunction orchestrates inflammasome effects on chronic obstructive pulmonary disease-induced depression: a potential mechanism underlying the cross talk between lung and brain[J]. Brain Behav Immun, 2019, 79: 195-206. doi: 10.1016/j.bbi.2019.02.001
[4] 陈琰, 钱频, 袁琳, 等. 慢性阻塞性肺疾病肺外合并症的治疗现状和研究进展[J]. 中华肺部疾病杂志: 电子版, 2020, 13(1): 100-103. doi: 10.3877/cma.j.issn.1674-6902.2020.01.025 [5] 苏建花, 程克文, 陈小东. 脑钠肽在慢性阻塞性肺疾病中的应用进展[J]. 医学综述, 2017, 23(23): 4677-4682. doi: 10.3969/j.issn.1006-2084.2017.23.019 [6] 王芳, 沈文超, 胡笑楠, 等. 皮质-基底神经节-丘脑-皮质环路与抑郁症状调节的机制探讨[J]. 临床精神医学杂志, 2021, 31(4): 326-329. doi: 10.3969/j.issn.1005-3220.2021.04.024 [7] 俞元临, 肖泽彬, 张桦, 等. 阿尔茨海默病及轻度认知障碍所致脑灰质体积及脑功能改变[J]. 中国医学影像技术, 2021, 37(2): 200-205. https://www.cnki.com.cn/Article/CJFDTOTAL-ZYXX202102013.htm [8] 陈苑, 姜雨, 陈艺, 等. 首发早发性抑郁症患者脑灰质体积及结构协变网络的异常改变[J]. 中华放射学杂志, 2021, 55(9): 941-947. doi: 10.3760/cma.j.cn112149-20201126-01260 [9] VERBERKT C A, VANDEN BEUKEN-VAN EVERDINGEN M H J, DIRKSEN C D, et al. Cost-effectiveness of sustained-release morphine for refractory breathlessness in COPD: a randomized clinical trial[J]. Respir Med, 2021, 179: 106330. doi: 10.1016/j.rmed.2021.106330
[10] 彭兆辉, 杨林峰, 彭兆伟, 等. 基于VBM研究中度至重度COPD患者脑灰质体积变化[J]. 医学影像学杂志, 2019, 29(11): 1832-1836. https://www.cnki.com.cn/Article/CJFDTOTAL-XYXZ201911005.htm [11] 刘春梅, 王春美, 吴立山, 等. 基于体素的形态学方法对抑郁症患者的脑结构磁共振研究[J]. 中风与神经疾病杂志, 2021, 38(7): 606-608. https://www.cnki.com.cn/Article/CJFDTOTAL-ZFSJ202107006.htm [12] CHIN FATT C R, JHA M K, MINHAJUDDIN A, et al. Dysfunction of default mode network is associated with active suicidal ideation in youths and young adults with depression: findings from the T-RAD study[J]. J Psychiatr Res, 2021, 142: 258-262. doi: 10.1016/j.jpsychires.2021.07.047
[13] 鲁毅. 首发未治疗抑郁症患者脑结构的多模态影像学研究[D]. 昆明: 昆明医科大学, 2017. [14] BERBOTH S, MORAWETZ C. Amygdala-prefrontal connectivity during emotion regulation: a meta-analysis of psychophysiological interactions[J]. Neuropsychologia, 2021, 153: 107767. doi: 10.1016/j.neuropsychologia.2021.107767
[15] STOODLEY C J, SCHMAHMANN J D. Functional topography in the human cerebellum: A meta-analysis of neuroimaging studies[J]. Neuroimage, 2009, 44(2): 489-501. doi: 10.1016/j.neuroimage.2008.08.039
[16] 闫晓培, 石宝玉, 王慧, 等. 慢性阻塞性肺疾病伴高血压病的临床特征和危险因素分析[J]. 实用临床医药杂志, 2021, 25(3): 63-68. https://www.cnki.com.cn/Article/CJFDTOTAL-XYZL202103016.htm [17] 邹君惠, 高圆圆, 秦若梦, 等. 高血压病患者脑白质高信号体积与认知功能的相关性分析[J]. 中华神经科杂志, 2019, 52(3): 190-196. [18] YU L C, DE MAZANCOURT M, HESS A, et al. Functional connectivity and information flow of the respiratory neural network in chronic obstructive pulmonary disease[J]. Hum Brain Mapp, 2016, 37(8): 2736-2754. doi: 10.1002/hbm.23205
计量
- 文章访问数: 192
- HTML全文浏览量: 75
- PDF下载量: 15
苏公网安备 32100302010246号