阿片使用障碍易感性遗传变异影响的解析

《Frontiers in Cellular Neuroscience》：Deciphering the impact of genetic variants on vulnerability to Opioid Use Disorder

　　

引言：阿片使用障碍（Opioid Use Disorder，OUD）是一种以强迫性阿片摄入为特征的慢性疾病，造成广泛健康、社会和经济负担。方法：为阐明成瘾易感性的分子贡献因素，研究人员对OUD个体与匹配对照的死后核壳区（nucleus accumbens，NAc）组织进行了全面RNA测序（RNA-sequencing，RNA-seq）分析。结果：队列一致性过滤在OUD样本中鉴定出16个基因上的17个候选错义变异和1个候选无义获得变异；这些RNA-seq检出的表达外显子变异（RNA-detected expressed exonic variants，RDEVs）需DNA水平验证。错义变异可破坏关键蛋白结构域，其中FUT9、FMR1、MFN1、RYR3和DAG1五个变异既往与物质使用障碍（Substance Use Disorders，SUDs）性状相关。唯一无义获得变异ZNF117预计通过折叠受损产生截短蛋白。转录组分析后经Ingenuity通路分析（Ingenuity Pathway Analysis，IPA）预测神经发育程序激活，145个上调基因和29个下调基因共同涉及突触可塑性、神经元分化和轴突导向通路。重要的是，鉴定出LINC01554和LINC00996等长效非编码RNA（long non-coding RNAs，lncRNAs）差异表达，并与NPAS4和GADD45B等关键转录因子存在推定调控关联。讨论：综上，这些发现提供了OUD中候选遗传与转录组改变的整合视图，并将以lncRNA为中心的调控网络及候选携带变异基因作为未来功能和DNA水平验证研究的假设生成线索。

研究背景方面，OUD是一种慢性物质使用疾病，其从物质使用到成瘾的进程涉及伏隔核（nucleus accumbens，NAc）等奖赏回路区域。尽管双生子研究与大规模基因组分析显示OUD具有强遗传力，且全基因组关联研究（Genome-Wide Association Studies，GWAS）已报道OPRM1、FURIN等位点，但DNA水平遗传变异如何驱动NAc等奖赏脑区下游转录组改变的机制仍不清楚。此前仅有一项RNA测序（RNA-sequencing，RNA-seq）研究探索阿片暴露者腹侧中脑变异，NAc区域尚缺同类研究。因此，研究人员利用公开存档的死后NAc的RNA-seq数据，开展变异发现与转录组整合分析，以生成候选机制假说。

关键技术方法方面，研究人员使用Seney等2021年发布的40例死后NAc样本队列，其中OUD组20例、非使用对照20例。原始FastQ数据来自序列读段档案（Sequence Read Archive）项目PRJNA739548。表达定量采用nf-core/rnaseq流程，经FastQC、Trim Galore!、SortMeRNA、STAR比对、UMI-tools与Picard去重、Salmon定量后，用DESeq2进行差异表达分析，设计公式纳入年龄、性别、种族、死后间隔（Postmortem Interval，PMI）与RNA完整值（RNA Integrity Number，RIN）。变异发现采用队列一致性过滤，保留全部OUD样本存在而对照不存在的表达外显子变异（RNA-detected expressed exonic variants，RDEVs）。结构建模使用ColabFold版AlphaFold2对ZNF117野生型与截短型预测。结构域注释用NCBI保守结构域数据库（Conserved Domain Database，CDD）的CD-Search。lncRNA与mRNA互作用lncRRIsearch与RNAhybrid预测。通路分析用Ingenuity Pathway Analysis（IPA）。

3.1 Identification of cohort-concordance filtered RDEVs

研究人员对20例OUD与20例对照NAc样本做个体变异调用与合并，保留所有OUD个体共有、对照全无的RDEVs，得到16个基因上18个单核苷酸RDEVs。部分候选变异人群等位基因频率（Allele Frequency，AF）高于0.8，提示RNA-seq检测偏差可能，因此高AF候选仅作透明报告，机制解释侧重低AF候选；这些RDEVs是需DNA验证的假设生成发现。

3.2 Characterization of RDEV impact on protein domains

研究人员将18个RDEVs映射到蛋白结构域，发现ANKRD36变异位于无序区，DAG1影响 dystroglycan-1信号肽，EXOC1影响 pleckstrin同源域，FLNB影响INPPL1结合界面，FMR1影响RANBP9介导核输出无序区，FUT9影响糖基转移酶供体结合口袋，MFN1影响动力蛋白样GTP酶域，RYR3影响重复单元，ZNF117破坏C2H2锌指基序，USP37影响泛素羧基末端水解酶域等。结论为这些RDEVs可能扰乱蛋白稳定性、膜运输与蛋白互作。

3.3 Association of RDEVs with substance use disorders

研究人员从17个错义RDEVs中聚焦FUT9、FMR1、MFN1、RYR3、DAG1五个有文献支持的基因，指出FUT9关联焦虑与神经连接，FMR1关联奖赏冲动与突触可塑性，MFN1关联线粒体融合与钙依赖突触活动，RYR3关联神经元兴奋性与学习，DAG1关联抑制性突触形成与可塑性，说明这些基因可作为SUDs表型候选。

3.4 Functional analysis of a stop-gain variant in ZNF117

研究人员分析ZNF117染色体7:64978289的c.1282C>T、p.R428^*无义变异，AlphaFold2显示截短体保留约88.44%序列但丢失两个C端C2H2锌指，影响DNA结合与转录调控；ZNF117受FOSL2、JUND、NR1F1、OTUD2等即早转录因子调控，截断可能扰乱神经发育、奖赏处理与记忆形成相关网络。

3.5 Differentially expressed genes and pathway analysis in OUD vs. non-users

研究人员以|log₂折叠变化|>1且p<0.05为标准，鉴定145个上调与29个下调差异表达基因（Differentially Expressed Genes，DEGs），包括IL1B、IL6、SOCS3等SUDs相关基因。IPA显示上调基因富集免疫与炎症、细胞因子、S100家族、糖皮质激素受体、雌激素受体信号；IL6与IL6ST存在已知互作。结论为NAc免疫失调与炎症相关转录重编程可能是OUD神经生物学特征，但bulk RNA-seq无法区分细胞内在转录上调与免疫细胞比例变化。

3.6 lncRNA-mRNA regulatory networks underlying SUDs

研究人员预测lncRNA-mRNA互作并聚类，发现LINC00397、LINC01219、LINC01554、LINC02605上调时靶基因SERPINA3、SIGLEC9、ZFP36、DLX4等同上调；LINC00996、LINC01657、LINC02982下调时多个炎症与应激基因反向上调，NGFR为例外。LINC01554与NPAS4、LINC00996与GADD45B存在互作，二者均既往关联SUDs。结论为lncRNA介导网络是SUDs转录失调的重要调控层。

讨论部分总结方面，研究人员强调18个RDEVs是队列一致性过滤下的候选集，不是已验证胚系易感位点；FLNB、EXOC1、RYR3分别指向细胞骨架、囊泡栓系与钙信号模块。FUT9、FMR1、MFN1、RYR3、DAG1兼具神经发育与SUDs表型证据。ZNF117无义变异虽AF约0.88，若DNA验证成立，其截短会丧失锌指介导转录抑制并可能影响Notch相关神经可塑性。DEG分析显示免疫炎症上调伴调控机制抑制，IL6/IL6ST轴连接神经免疫与成瘾表型。lncRNA网络进一步把LINC01554、LINC00996、NPAS4、GADD45B带入SUDs调控模型。研究限制包括RNA-seq只能检测表达外显子区、不能区分胚系与体细胞、受RNA编辑与覆盖度影响、种族匹配不完全、bulk RNA-seq无细胞类型分辨力。

研究结论部分翻译为：本研究中，研究人员在OUD个体NAc中鉴定出一致存在而匹配非使用者中完全缺失的18个变异，注释了ANKRD36、DAG1、EXOC1、FLNB、FMR1、FUT9、HEATR5A、IL6ST、MFN1、NEMF、PTBP2、RYR3、TMCC1、TTN、ZNF117和USP37共16个基因；除一个无义获得变异外均为错义变异，预计中度影响蛋白功能与成瘾相关神经信号通路。研究人员继而检查受影响结构域、阐述变异与SUDs表型含义、表征无义变异对ZNF117蛋白形成的影响、分析DEGs及其与变异基因和SUDs基因互作，并检测DEGs中lncRNA及其调控互作。这些候选RDEVs将细胞骨架支架、囊泡栓系机器与钙信号模块列为阿片致病机制中值得功能随访的方向，但都需DNA水平确认。该论文发表于《Frontiers in Cellular Neuroscience》。