尽管镇静催眠药使用障碍（SHUD）因其高度滥用潜力已成为一个主要的公共卫生优先事项，但个体易感性的潜在神经生物学因果机制在很大程度上仍属未知。研究人员利用来自英国生物银行（UK Biobank）（N = 33，000名参与者）的3935个脑影像学衍生表型（IDPs）的全基因组关联研究（GWAS）汇总统计数据，以及来自FinnGen R10版本（N = 2258例）的SHUD数据，开展了一项大规模的双向双样本孟德尔随机化（MR）研究。主要采用逆方差加权法（IVW）估算因果效应，并辅以一系列敏感性分析，包括MR-Egger回归、加权中位数法和MR-PRESSO。通过Benjamini-Hochberg错误发现率（FDR）校正控制多重检验。正向MR分析发现了34个与SHUD风险存在显著因果关联的脑IDPs。具体而言，易感性增高与颞叶和边缘区域的结构指标相关，包括左侧颞横沟（OR = 1.35， 95% CI 1.09–1.67， FDR = 0.04）和右侧颞中回（OR = 1.30， 95% CI 1.07–1.58， FDR = 0.04）皮质厚度增加，以及右侧脑岛（OR = 1.31， 95% CI 1.13–1.51， FDR = 0.02）和左侧海马旁回（OR = 1.29， 95% CI 1.08–1.54， FDR = 0.04）白质-灰质信号强度对比度（WG contrast）升高。相反，脑干体积（OR = 0.74， 95% CI 0.61–0.91， FDR = 0.04）和左侧扣带束（OR = 0.73， 95% CI 0.57–0.92， FDR = 0.04）的白质完整性则表现出保护效应。功能连接分析显示，突显网络（salience network）中较高的节点活动（OR = 0.70， 95% CI 0.53–0.92， FDR = 0.04）以及视觉网络与躯体运动网络之间增强的连接（OR范围：0.54–0.74）预示着风险降低。相比之下，默认模式网络（default mode network）和视觉网络内的特定功能连接边则与SHUD易感性呈正相关（OR分别为1.56和1.48）。反向MR分析提供的证据有限，不支持遗传预测的SHUD对右侧颞中回皮质厚度存在因果效应。虽然观察到名义上的关联（OR = 1.03， 95% CI 1.00–1.05， p = 0.02），但该发现在多重检验校正后并未存活（FDR = 0.56），并且在不同的敏感性模型中缺乏稳健性。本研究确立了SHUD易感性与特定皮质形态、皮质下体积和功能网络拓扑结构之间的因果联系，强调了时间-边缘系统（temporal-limbic systems）和突显系统的关键作用。这些发现为理解SHUD的神经生物学基础提供了新的机制见解，并为风险分层和治疗靶点的开发识别了潜在的影像学生物标志物。
**论文解读**
镇静催眠药物，尤其是苯二氮䓬类药物（BZDs）和Z类药物，是已知滥用潜力与依赖风险最高的药物类别之一。其作用机制主要源于对大脑抑制性γ-氨基丁酸A型（GABA
A）受体的调节。临床证据表明，尽管此类药物在短期治疗失眠及共病精神障碍中有效，但长期或不当使用极易导致依赖和成瘾。随着处方率和用量的持续上升，相关的急诊就诊和过量致死事件，特别是涉及BZD的事件，已成为严峻的公共卫生挑战。药物成瘾被认为是一种慢性、复杂的大脑疾病，涉及一个被称为“成瘾环路”的关键脑区网络。尽管传统神经影像学研究已在多种物质使用障碍中发现了执行控制网络、默认模式网络和突显网络的连接异常，但针对镇静催眠药使用障碍（SHUD）的系统性研究仍十分匮乏。传统观察性研究受限于样本量、混杂因素和因果方向性不明等问题，难以阐明影像学衍生表型（IDPs）与SHUD神经功能异常之间的因果关系。孟德尔随机化（MR）作为一种利用遗传变异作为工具变量（IVs）的强大流行病学框架，能够有效控制混杂和反向因果，为评估暴露与疾病结局间的因果关系提供了有力工具。此前，MR研究已成功揭示了物质滥用与大脑功能改变间的因果联系。然而，SHUD与大脑结构、功能特征间的因果关系尚未得到充分阐明。因此，本研究旨在利用双向双样本MR方法，系统分析3935个脑IDPs与SHUD之间的因果关联，以厘清大脑形态学特征与SHUD之间的潜在双向影响，为该疾病的预防和临床管理提供新的机制见解。该论文发表在《Addiction Biology》杂志上。
为开展此研究，研究人员主要运用了以下关键技术方法：首先，利用来自英国生物银行（UK Biobank）的GWAS汇总数据作为脑IDPs的暴露数据集（N ≈ 33,000），以及来自芬兰FinnGen研究（R10版本）的GWAS汇总数据作为SHUD的结局数据集（N = 2,258例）。其次，采用严格的流程筛选MR分析的遗传工具变量（IVs），包括设定显著性阈值（p < 1×10
-5）、进行连锁不平衡（LD）剔除（阈值 < 0.05），并计算F统计量以排除弱工具偏倚。核心分析采用双向双样本MR框架，以逆方差加权法（IVW）作为主要因果效应估计方法，并辅以MR-Egger、加权中位数法、加权众数法等多种敏感性分析方法进行稳健性检验。最后，通过MR-PRESSO全局检验、Cochran's Q统计量等评估异质性和水平多效性，并通过留一法分析评估单个SNP对结果的驱动性。
**3.1 基线特征**
研究队列包含2258名镇静催眠药使用者，其中男性1106人（49%），女性1152人（51.0%）。精神行为障碍的平均发病年龄为40.56岁（男性41.24岁，女性39.91岁）。发病年龄分布以20-30岁（26.3%）、30-40岁（22.7%）和40-50岁（19.4%）年龄段为主。在随访期间，与未暴露相比，使用镇静催眠药与多种神经系统和精神障碍的风险显著增高相关。其中，风险增加最显著的疾病包括阿尔茨海默病（HR = 5.29）、癫痫（HR = 5.00）、精神分裂症（HR = 18.29）和双相情感障碍（HR = 16.45）。这些结果提示SHUD与广泛的神经精神共病风险增高存在关联。
**3.2 工具变量（IV）选择**
对于正向MR分析，研究人员为每个脑IDP提取了有效的遗传工具变量（IVs）。筛选标准包括：单核苷酸多态性（SNP）与暴露（脑IDPs）显著关联（p < 10
-5）、LD剔除（r
2 < 0.01，窗口大小1 Mb）以及排除弱工具（F统计量 < 10）。最终为每个脑IDP获得了8至103个有效IVs。对于反向MR分析，应用类似的筛选标准（p < 10
-5，r
2 < 0.05，窗口大小1 Mb，F统计量 > 10），获得了31个与SHUD显著关联的IVs。
**3.3 正向MR分析：脑IDPs对SHUD的潜在因果效应**
在正向MR分析中，将脑IDPs作为暴露，SHUD作为结局。在评估的3935个IDPs中，有177个成功匹配了足够数量的有效IVs进行MR分析。经FDR校正（q < 0.05）后，34个脑IDPs显示出与SHUD的显著因果关联。其中，14个为结构IDPs，映射到特定的解剖区域，包括颞叶、顶叶和枕叶皮质、脑干、脑岛、边缘系统及多种白质束。其余20个显著IDPs为功能连接指标，分布于多个脑网络，包括突显网络、视觉网络、默认模式网络、听觉网络和躯体运动网络。
具体而言，颞叶皮质的改变与SHUD风险显著相关。两个皮质厚度指标显示出正向关联：左侧颞横沟（lh S-temporal-transverse）厚度每增加一个标准差（SD），与SHUD风险增加35%相关（OR = 1.35， 95% CI 1.09–1.67， FDR = 0.04）；右侧颞中回（rh G-temporal-middle）厚度每增加1 SD，与风险增加30%相关（OR = 1.30， 95% CI 1.07–1.58， FDR = 0.04）。相反，右侧颞平面（right planum temporale）的体积每增加1 SD，则与风险降低17%相关（OR = 0.83， 95% CI 0.73–0.95， FDR = 0.04）。
在边缘和旁边缘区域，白质-灰质信号强度对比度（WG contrast）与SHUD风险呈正向关联。左侧海马旁回（wg lh intensity-contrast parahippocampal）的WG contrast每增加1 SD，与SHUD的患病几率增加29%相关（OR = 1.29， 95% CI 1.08–1.54， FDR = 0.04）。此外，右侧脑岛（wg rh intensity-contrast insula）较高的WG contrast也与风险增高31%相关（OR = 1.31， 95% CI 1.13–1.51， FDR = 0.02），凸显了边缘系统结构完整性在SHUD易感性中的潜在作用。
在顶叶和枕颞区域，多个解剖图谱下的关联均保持稳健。在顶叶皮质，右侧楔前叶（right precuneus）较大的表面积与SHUD风险增加相关（OR = 1.29， 95% CI 1.08–1.53， FDR = 0.04）。舌回（lingual gyrus）的发现尤为一致；该区域体积增大在Desikan（左侧：OR = 1.29， FDR = 0.04）、DKT（左侧：OR = 1.30， FDR = 0.04）和a2009s（右侧：OR = 1.46， FDR = 0.03；左侧内侧枕颞/舌回：OR = 1.29， FDR = 0.04）分区中均被确定为显著的风险因素。此外，右侧中枕回和月状沟（rh S-oc-middle + Lunatus）皮质厚度的增加也与SHUD风险呈正向关联（OR = 1.36， 95% CI 1.09–1.69， FDR = 0.04）。
相比之下，在皮质下结构和白质束中观察到了保护效应。对于皮质下区域，脑干体积每增加1 SD，与SHUD风险降低26%相关（OR = 0.74， 95% CI 0.61–0.91， FDR = 0.04）。关于白质微观结构，左侧扣带束扣带回（cgc L）的各向异性模式（MO）和右侧听辐射（ar R）的第三特征值（L3）均与SHUD风险呈负相关（OR分别为0.73和0.71， FDR = 0.04）。
在功能层面，脑网络架构（包括网络内活动和网络间连接）在FDR校正后与SHUD风险显著相关。在多个主要系统中发现了网络内关联。就节点活动而言，突显网络（SN）中较高的rfMRI振幅（ICA25节点11）与SHUD风险降低相关（OR = 0.70， 95% CI 0.53–0.92， FDR = 0.04）。在听觉网络（AN）内，两个功能连接边（ICA100边329和531）表现出强烈的保护效应（均p < 0.001， FDR < 0.05）。对于视觉网络（VN）和默认模式网络（DMN），发现是双向的：特定边显示出风险降低（VN边9：OR = 0.69； DMN边45：OR = 0.79），而其他边则与风险增高相关（VN边395：OR = 1.48； DMN边255：OR = 1.56）。此外，在躯体运动网络（SMN）（边519：OR = 0.64）和额顶网络（FPN）（边263：OR = 0.61）内的连接也观察到了保护性关联。
网络间功能连接分析也揭示了与SHUD风险的广泛关联，特别是涉及视觉网络（VN）。VN与SMN之间的众多连接一致显示出保护效应，包括边834、1118、1122、1386、1453和1455（OR范围从0.54到0.74，所有FDR < 0.05）。VN与DMN之间的交互也显著，边862与风险增高相关（OR = 1.45），但边991和1466显示出保护性关联。此外，SHUD风险与VN-FPN连接（边776：OR = 1.34）呈正相关，而与VN-小脑网络（CBN）连接（边807：OR = 0.65）呈负相关。
**3.4 反向MR分析：SHUD对脑IDPs的因果效应**
为研究潜在的反向因果性，研究人员采用反向MR分析，以SHUD作为暴露，先前识别的34个IDPs作为结局。在五种MR方法中，仅IVW方法识别出显著关联。具体而言，遗传预测的SHUD与右侧颞中回（rMTG）皮质厚度增加呈名义上显著关联（OR = 1.03， 95% CI 1.00–1.05， p = 0.02）。值得注意的是，该发现在多重检验校正后并未存活（FDR = 0.56），并且在MR-Egger和加权中位数法等替代敏感性模型中缺乏一致性。这些结果表明，没有充分证据支持SHUD对已识别的皮质形态学特征存在稳健的反向因果效应。
**3.5 敏感性分析**
为评估双向MR发现的稳健性，进行了一系列敏感性分析。留一法敏感性分析显示，汇总效应估计保持一致，并未受到任何单个SNP的过度影响，所有结果均持续位于零效应线（β = 0）的同一侧。MR-Egger截距检验未显示水平多效性的证据，截距均接近零，所有p > 0.05。此外，通过Cochran's Q和Rücker's Q统计量（针对IVW和MR-Egger方法）未检测到显著异质性（所有p > 0.05）。这些结果共同强化了本研究基于MR的因果推断的稳定性和可靠性。
**讨论**
尽管脑IDPs已在多种物质使用障碍中被广泛涉及，但其与SHUD的因果关系仍不明确。利用双向双样本MR框架，研究人员系统评估了3935个IDPs，识别出34个与SHUD风险存在因果关联的特定大脑结构和功能特征。与正向MR发现的稳健性形成鲜明对比的是，反向MR分析提供的证据极少，不支持SHUD引起大脑改变。这种显著的不对称性表明，大脑形态学特征与SHUD之间的关系主要是单向的，即先天存在的大脑结构和功能架构主要作为疾病易感性的因果驱动因素，而非疾病的结果。
本研究的核心发现是SHUD风险与时间-边缘区域结构改变之间的因果关联。镇静催眠药（主要是BZDs和Z类药物）主要通过作用于广泛分布于大脑皮质和海马的GABA
A受体发挥效应。研究结果显示，颞横回和颞中回皮质厚度的增加与更高的SHUD风险相关。虽然皮质变薄通常见于神经退行性疾病的晚期阶段，但本研究发现皮质增厚与SHUD风险增高相关，指向了一种独特的神经生物学易感性。这一现象可能可以用异常的神经发育轨迹来解释，特别是在突触修剪不足方面。在青春期晚期和成年早期——这也是本队列中发病的高峰期——大脑经历着冗余突触的广泛修剪。这一过程的失败或延迟可能导致皮质增厚但计算效率低下，这与冲动性增加和对成瘾行为的易感性增加有关。这与近期生物物理模型提出的观点一致，即正常皮质成熟以兴奋-抑制比的下降为标志；这一成熟过程的失败，可能表现为异常增厚的皮质，可能是神经认知脆弱性的标志。此外，鉴于镇静催眠药的靶点是GABA
A受体，较厚的颞叶皮质可能反映了更高的GABA能中间神经元或受体复合物的基础密度。这种增加的生物学基础可能会放大镇静药物的药理作用，强化觅药行为，并增加发展为使用障碍的风险。
具体而言，海马旁回和脑岛中白质-灰质（WG）对比度与SHUD风险的正向关联值得注意。WG对比度是皮质髓鞘化和结构完整性的一个敏感标志。脑岛是突显网络的关键枢纽，其功能是区分内部和外部环境中最重要的刺激以指导行为。该区域升高的WG对比度和非典型髓鞘化可能会干扰显著性映射过程。这种神经结构特征可能会使个体倾向于高估与药物相关的线索，同时对内部稳态耗竭的敏感性降低，这是成瘾的核心机制。近期的神经影像学研究强调，主要由前脑岛和前扣带回组成的突显网络在检测显著刺激方面发挥着关键作用；其结构和功能的改变反复与物质使用障碍个体的异常显著性归因和认知控制受损相关联。
相反，研究人员在皮质下结构和白质束中发现了保护性因素。较大的脑干体积与SHUD风险降低相关（OR = 0.74）。考虑到脑干在调节上行网状激活系统（ARAS）中的作用及其高密度的GABA能神经元，强健的脑干体积可能作为对抗镇静药物诱导的觉醒系统脱敏的生物学缓冲。此外，扣带束（OR = 0.73）和额顶网络（边263， OR = 0.61）的保护作用，强化了自上而下抑制控制的重要性。研究结果表明，这些通路中更优的结构连接可能增强抑制控制和情绪韧性，从而减轻从初始镇静药使用到慢性使用障碍的进程。这与相关证据一致，即对成瘾的抵抗力以涉及习惯调节和自上而下抑制控制的环路超连接为特征。
在功能层面，MR分析凸显了SHUD易感性以功能脑架构失调为特征。所识别的涉及突显网络、默认模式网络和视觉网络的因果关系支持了成瘾的“三网络模型”。突显网络中较高的振幅（ICA25节点11）具有保护性（OR = 0.70），表明稳定的固有突显网络活动能够实现从默认模式网络的内部渴求到额顶网络执行控制的高效转换。这与近期将脑岛鉴定为介导物质使用障碍中网络活动间不协调转换的因果枢纽的证据相符。值得注意的是，视觉网络成为一个核心连接枢纽。尽管在成瘾研究中常被忽视，但视觉网络与躯体运动网络之间持续的保护效应（例如边834和1118）表明，高效的传感器整合是一种韧性因素。这种协调可能促进更高阶的感觉过滤，减少视觉药物线索转化为冲动运动行为的可能性。
反向MR分析揭示了SHUD易感性与右侧颞中回皮质厚度增加之间的名义上关联（OR = 1.03）。尽管右侧颞中回对认知处理至关重要，但这一微小效应需要谨慎对待。该关联在敏感性分析中未表现出一致性，表明其稳健性有限，并存在假阳性发现的风险。因此，这些有限的发现仅提供了微弱的反向因果证据，不支持双向反馈机制。相反，稳健反向信号的缺失证实了正向MR的发现，强化了所识别的IDPs代表先天的神经生物学脆弱性，而非慢性物质暴露的继发性后果。证据表明，尽管慢性物质使用可以改变大脑结构，但许多形态学偏差，特别是与皮质厚度相关的偏差，可能先于物质使用障碍的发生而存在。
本研究存在一些局限性。首先，反向MR中用于SHUD的IVs强度低于用于脑IDPs的IVs强度，这可能限制了检测细微反向因果效应的能力。其次，数据主要源自欧洲血统人群，需要进一步验证跨人群的普遍性。第三，MR无法完全排除基因-环境相互作用，而镇静催眠药的使用在很大程度上受环境和社会心理因素的影响。第四，本研究关注的是静态影像表型；整合纵向数据、转录组信息和多模态成像将为理解镇静催眠药使用障碍的神经生物学机制提供更全面的见解。
**结论**
总而言之，本研究提供了遗传学证据，表明大脑架构对SHUD风险存在主要的单向因果影响。通过识别涉及时间-边缘系统和功能网络的34个潜在因果IDPs，研究结果表明，先天存在的神经生物学基质可能主要驱动疾病易感性，而非仅仅是药物使用的结果。尽管这些基于MR的推断存在固有的局限性，但它们支持了SHUD的先天脆弱性模型。总体而言，这项工作提供了一个初步的因果框架，并识别了候选的神经影像学生物标志物，这些标志物需要进一步的纵向验证，以用于风险分层和预防策略。
