《Brain Sciences》：NR3C1 Exon 1F Promoter Methylation and mRNA Expression in the Dorsolateral Prefrontal Cortex of Individuals Who Died by Suicide: A Post-Mortem Case–Control Study from Northern Mexico
编辑推荐：
背景：糖皮质激素受体基因（NR3C1）调控下丘脑–垂体–肾上腺（HPA）轴的负反馈调节。其在自杀中的表观遗传及转录改变的证据仍然有限且不一致，尤其是在代表性不足的人群中。方法：研究人员测量了来自墨西哥北部13例自杀死亡者和13例非精神科对照者的背外侧前额叶皮层
背景：糖皮质激素受体基因（NR3C1）调控下丘脑–垂体–肾上腺（HPA）轴的负反馈调节。其在自杀中的表观遗传及转录改变的证据仍然有限且不一致，尤其是在代表性不足的人群中。方法：研究人员测量了来自墨西哥北部13例自杀死亡者和13例非精神科对照者的背外侧前额叶皮层（BA9）中NR3C1外显子1F启动子区域内的区域甲基化水平及总NR3C1 mRNA表达。各组在年龄、性别和死后间隔（PMI）方面进行了频数匹配。甲基化通过甲基化敏感性限制酶（MSRE）qPCR评估，表达通过RT-qPCR评估。组间差异采用双侧精确Mann–Whitney U检验；效应量为优势概率（PS）和Cliff's δ。相关性和仅男性分析为探索性分析。所有p值均为名义p值。结果：自杀死亡者表现出更高的总NR3C1 mRNA表达（U = 45，p = 0.04412；PS = 0.734；δ = 0.467）和外显子1F区域甲基化（U = 29，p = 0.00349；PS = 0.828；δ = 0.657）。两项指标呈非显著正相关（ρ = 0.27，p = 0.179）。在男性中，甲基化在自杀死亡者中仍较高（p = 0.0147），而表达虽方向性升高但未达显著（p = 0.165）。结论：自杀死亡者在BA9中表现出更高的外显子1F区域甲基化和NR3C1 mRNA表达。两项指标无显著相关，应将其解释为独立的组间差异，而非甲基化–表达直接偶联的证据。表达结果属于初步发现。鉴于样本量小、临床信息有限以及与无精神科病史的对照比较，这些差异不能特异性地归因于自杀，需进一步重复验证。
论文解读
一、研究背景与目的
自杀是全球重大公共卫生问题，每年导致超过72万人死亡，是全球15–29岁人群的第三大死因。在墨西哥，自杀率在过去十年持续上升，对青少年和年轻成人影响尤为突出，给精神卫生保健系统带来日益沉重的负担。尽管情绪障碍的诊断和药物治疗已取得显著进展，但驱动自杀行为的生物学机制仍未完全阐明，这一认知空白持续阻碍基于生物学的有效预防策略的开发。
下丘脑–垂体–肾上腺（HPA）轴失调已被确定为应激相关精神疾病（包括重性抑郁障碍和自杀行为）的关键病理生理特征。HPA轴调控的核心是糖皮质激素受体（GR），由NR3C1基因编码，通过响应循环皮质醇调节负反馈。表观遗传改变（尤其是NR3C1启动子区域内的DNA甲基化）可影响基因转录和受体对糖皮质激素的敏感性。大量研究已将NR3C1启动子甲基化与早年逆境、创伤和自杀倾向相关联，凸显其作为应激脆弱性分子中介的潜在作用。
NR3C1在大脑中的调控具有高度区域和细胞类型特异性。例如，在有儿童虐待史的自杀死亡者海马组织中，NR3C1外显子1F启动子甲基化升高与GR-mRNA表达降低相关；另一项针对自杀死亡者前额叶皮层的研究则报告了NR3C1某一转录本表达降低。背外侧前额叶皮层（dlPFC），包括Brodmann 9区（BA9），在执行控制、决策和应激的认知评估中发挥核心作用——这些领域均与自杀行为密切相关。Rizavi等人曾检测青少年自杀死亡者死后BA9中NR3C1外显子1甲基化、羟甲基化及转录本特异性表达，为该区域NR3C1表观遗传与转录调控的整合研究确立了重要先例。然而，来自独立人群和替代分析方法的证据仍然有限。
本研究检测了墨西哥北部杜兰戈省自杀死亡者和非精神科对照者BA9样本中NR3C1外显子1F启动子区域甲基化和总NR3C1 mRNA表达。研究人员假设自杀死亡者会表现出更高的区域甲基化和更低的NR3C1表达，并探索了这些指标之间的关系。本研究通过提供来自独立墨西哥北部队列的数据，扩展了先前研究。
二、研究方法概述
本研究为回顾性病例对照设计，纳入13例自杀死亡者和13例非精神科对照者，两组在年龄、性别和死后间隔（PMI）方面进行组级频数匹配。所有样本均来自墨西哥杜兰戈法医服务中心（SEMEFO），经杜兰戈总医院研究伦理委员会批准，样本采集前获得每位死者近亲属的书面知情同意。尸检时由法医病理学家从左侧背外侧前额叶皮层解剖鉴定BA9区域，样本包含皮层灰质及邻近白质，以整体组织方式处理。
甲基化分析采用甲基化敏感性限制酶（MSRE）qPCR方法，使用OneStep? qMethyl PCR试剂盒。研究人员通过MethPrimer设计了一对引物，扩增覆盖NR3C1外显子1F及其上游启动子区域的264 bp片段（GRCh38/hg38，chr5:143,404,220–143,404,483，反向链），扩增子内包含29个CpG位点和16个MSRE识别基序。表达分析采用TaqMan?基因表达检测试剂盒，以B2M作为内参基因，通过2?ΔΔCt方法计算相对表达量。统计分析使用SPSS 25.0和R 4.5.0，采用双侧精确置换Mann–Whitney U检验比较组间差异，效应量以优势概率（PS）和Cliff's δ表示，相关性分析采用Spearman秩相关。
三、研究结果
3.1 队列特征与人口学及死因
研究共纳入13例自杀死亡者和13例非精神科对照者。每组各含10名男性和3名女性。两组在年龄、性别分布和PMI方面无统计学显著差异。自杀死亡者中，窒息为主要方式（n = 10，76.9%），中毒、枪伤和药物过量各1例（各占7.7%）。对照组死因包括交通事故相关（n = 7，53.8%）、多发创伤（n = 4，30.8%）、溺水（n = 1，7.7%）和枪伤（n = 1，7.7%）。
3.2 前额叶皮层甲基化–表达谱
自杀死亡者前额叶皮层样本表现出显著高于非精神科对照者的NR3C1 mRNA表达（U = 45，p = 0.04412）。同时，自杀死亡者NR3C1启动子区域（外显子1F）的DNA甲基化也显著升高（U = 29，p = 0.00349）。探索性Spearman分析显示NR3C1启动子甲基化与mRNA表达之间存在非显著正相关（ρ = 0.27，p = 0.1792）。
3.3 效应量与敏感性分析
自杀组个体NR3C1表达值高于随机对照个体的估计概率为0.734，对应Cliff's δ为0.467。外显子1F启动子区域甲基化的PS为0.828，对应Cliff's δ为0.657。敏感性分析表明，每组13人、双侧α = 0.05的条件下，80%效能可检测的最小效应对应PS = 0.817，即Cliff's δ = 0.634。因此，现有样本主要对较大的病例–对照差异敏感。
3.4 探索性仅男性分析
由于女性数量有限（每组n = 3），未进行女性特异性比较和组别×性别交互分析。在男性中（每组n = 10），NR3C1表达在自杀死亡者中仍较高但未达统计学显著（p = 0.165），而启动子甲基化在自杀死亡者中仍显著升高（p = 0.0147）。男性中区域外显子1F甲基化与总NR3C1 mRNA表达呈非显著正相关（Spearman's ρ = 0.215，p = 0.363）。
3.5 按自杀方式的探索性分析
NR3C1表达和启动子甲基化按自杀方式进行了描述性检查。窒息占10例，中毒、枪伤和药物过量各1例。鉴于亚组规模小且高度不平衡，未进行正式的方法间比较。
四、讨论与结论
在BA9区域，自杀死亡者表现出更高的NR3C1外显子1F启动子区域甲基化和更高的总NR3C1 mRNA表达。NR3C1 mRNA表达的升高应视为初步发现。甲基化与表达无显著相关，应解释为独立的组间差异。由于窒息是主要自杀方式，不能排除缺氧或其他方式相关的濒死效应的影响，但本研究设计无法确定其作用。
在糖皮质激素假说框架下，HPA轴反馈受损与GR/NR3C1信号减弱相关。McGowan等人在有儿童虐待史的自杀死亡者海马组织中报告了NR3C1外显子1F启动子甲基化升高和总及外显子1F特异性表达降低。Rizavi等人检测了青少年自杀死亡者的BA9和海马，在BA9中观察到GR-1F转录本表达降低，而海马中未观察到。其甲基化分析聚焦于BA9，发现与外显子1F启动子重叠的区域5mC升高、5hmC降低。与研究人员先验假设的甲基化部分一致，自杀死亡者的外显子1F区域甲基化较高。然而，表达结果与假设方向相反：自杀死亡者的总NR3C1 mRNA表达更高而非更低。这一意外发现的机制尚不清楚，因为未测量外显子1F特异性转录本、GR蛋白水平和受体功能。
可能导致这一意外组合的因素包括NR3C1复杂的5′调控结构，其包含多个由不同启动子控制的可变非翻译第一外显子。在成人人脑dlPFC中，含外显子1B、1C和1H的转录本高表达，而外显子1F转录本相对较少。由于本研究检测的是总NR3C1表达而非启动子特异性转录本，局限于外显子1F区域的甲基化可能不与总mRNA丰度平行。不同研究间的差异可能同样反映了所检测的特定启动子和转录本的不同。
整体组织样本的使用限制了结果的细胞学解释。甲基化和表达测量均代表神经元和多种胶质细胞群的聚合信号，这些细胞可能具有不同的分子特征。因此无法确定观察到的甲基化和表达差异是否源于相同细胞群，或细胞组成变异是否促成了这些指标间的表观分离。
NR3C1调控可能涉及靶向甲基化检测未捕获的分子过程。基因转录受三维染色质结构和远端增强子–启动子相互作用影响，但本研究未检测这些机制。此外，NR3C1 mRNA丰度不一定反映GR蛋白丰度或功能活性。FKBP5编码的共伴侣蛋白FKBP51可降低GR依赖性转录活性并延迟受体核转位，而不改变GR蛋白水平。
既往BA9研究曾报告自杀死亡者GR表达降低和多个NR3C1外显子1调控区域甲基化升高。相比之下，本研究观察到总NR3C1 mRNA表达升高与外显子1F区域甲基化升高并存。这些差异应谨慎解读，因为各研究考察了不同人群、调控区域和转录本靶标，并使用了不同的甲基化评估方法。转录组学综述也报告了HPA轴相关变化在不同皮层区域的方向和幅度差异。
在仅男性分析中，启动子甲基化在自杀死亡者中仍较高，而表达差异方向不变但不再具有统计学显著性。甲基化差异的持续存在与既往仅男性死后研究报告的NR3C1甲基化改变大体一致，但涉及的脑区、临床背景和分子变化方向在不同研究间存在差异。然而，该发现并不确立性别特异性效应。
这些发现提示，局限于单一NR3C1区域的甲基化测量可能无法捕捉GR调控的复杂性。未来整合启动子分辨甲基化、转录本特异性表达、GR蛋白水平、FKBP5、功能性HPA轴指标以及可及外周组织中相应评估的研究，可能有助于阐明这些改变的生物学和转化意义。
研究结论
自杀死亡者在dlPFC中表现出更高的NR3C1外显子1F启动子区域甲基化和更高的总NR3C1 mRNA表达。由于甲基化与表达无显著相关，这些发现应解释为平行的组间差异，而非检测区域内甲基化直接促进转录的证据。表达发现应视为初步结果，不确立GR蛋白丰度或受体信号增强。这些发现通过提供来自独立且代表性不足的墨西哥北部死后队列的分子数据，扩展了先前研究。然而，不能排除方式相关濒死效应和未测量的精神科或临床因素的贡献。需要在更大规模、临床特征明确的队列中，使用细胞类型和启动子分辨的甲基化、转录本特异性表达及GR蛋白测量进行重复验证，以阐明该模式的生物学意义。
