全基因组关联联合分析揭示99个疼痛敏感性风险位点及其与精神、代谢和免疫性状的多效性关系

Genome wide association joint analysis reveals 99 risk loci for pain susceptibility and pleiotropic relationships with psychiatric, metabolic, and immunological traits

作者信息Evelina Mocci, Kathryn Ward, James A Perry, Angela Starkweather, Laura S Stone, Siobhan M Schabrun, Cynthia Renn, Susan G Dorsey, Seth A Ament
PMID37844115
发布时间2023-10-16
DOI10.1371/journal.pgen.1010977

实验完整度

研究为大规模GWAS联合分析及后续生物信息学分析,包含多层级验证(MTAG、独立队列复制、遗传相关、MR),但缺乏直接的功能实验验证。

主要模型

UK Biobank队列(194,174例女性,167,020例男性) 人类疼痛转录组数据集(64份血液样本,包括急性、慢性腰痛患者和健康对照)

重点核对

17个疼痛相关表型的病例/对照数量及定义 PLEIO联合分析中使用的MAF≥0.01和HWE p>0.00001的过滤标准 用于基因定位和功能注释的FUMA参数(r²阈值、LD blocks) 用于孟德尔随机化的工具变量选择标准(F统计量>10)

摘要

慢性疼痛在美国呈流行趋势,对公共健康系统造成重大负担,并与日益严重的阿片类药物危机相伴。虽然已有大量针对特定疼痛相关性状的全基因组关联研究(GWAS)报道,但其中许多研究统计功效不足,且这些性状之间的遗传关系仍知之甚少。在此,我们对英国生物银行(UK Biobank)中17个疼痛敏感性性状的GWAS汇总统计进行了联合分析。该分析揭示了99个全基因组显著风险位点,其中65个与早期研究中发现的位点重叠,其余34个为新位点。我们采用留一性状法荟萃分析评估每个性状对联合分析的影响,结果表明这些位点可分为四类:与几乎所有疼痛相关性状相关的位点;主要与单一性状相关的位点;与多种骨骼肌疼痛相关的位点;以及与头痛相关疼痛相关的位点。总体而言,通过基因组邻近性、eQTL和染色质互作,将664个基因定位到99个位点上,其中约15%的基因在急性或慢性疼痛个体中与健康对照相比表现出差异表达。风险位点富集于涉及神经和炎症通路的基因。遗传相关性和双样本孟德尔随机化表明,精神、代谢和免疫性状介导了部分效应。

实验结论

提炼研究问题、关键发现与证据,快速把握文章的核心贡献。

研究问题
疼痛敏感性是否存在共享的遗传结构?哪些风险位点与多种疼痛性状相关?
核心机制
疼痛风险位点富集于神经(如神经发生)和免疫(如趋化因子分泌、T细胞分化)相关基因,提示疼痛敏感性涉及神经和免疫双重机制。
主要证据
基于UK Biobank 17个性状的PLEIO联合分析鉴定99个风险位点,通过MTAG、留一分析及独立队列(偏头痛、芬兰疼痛)复制验证。
研究意义
这些风险位点为开发新的非成瘾性疼痛药物提供了潜在治疗靶点,并支持疼痛易感性由神经和免疫过程共同介导。

研究路径

按研究推进顺序梳理实验设计、验证步骤与关键观察。

1

GWAS表型选择与质量控制

选择具有显著遗传力的疼痛相关性状,并进行数据质量控制。

从UK Biobank中筛选17个疼痛相关性状,基于样本量>10,000和LDSC遗传力z>4,并剔除遗传相关rg<0.35的性状;对常染色体SNP进行MAF和HWE过滤。

2

PLEIO联合分析

利用多性状联合分析提高统计功效,发现多效性风险位点。

使用PLEIO对17个性状进行联合分析,并分别进行男性和女性分层分析,以及13性状敏感性分析。

3

风险位点鉴定与基因定位

定义LD独立的基因组风险位点,并定位候选基因。

使用FUMA SNP2GENE定义风险位点,并进行基因定位(位置、eQTL、染色质互作)。

4

基因集富集分析

识别风险位点富集的基因功能类别和组织表达特异性。

使用MAGMA进行基因集富集分析和组织特异性表达分析。

5

遗传相关与孟德尔随机化

评估疼痛易感性与非疼痛性状的遗传相关及潜在因果关系。

使用LDHub进行遗传相关分析,对8个性状进行双样本孟德尔随机化分析(包括反向MR和MVMR)。

6

独立队列复制

在独立队列中验证风险位点的关联。

分别在偏头痛GWAS(29,209例病例)和芬兰疼痛GWAS(151,484例病例)中进行复制分析。

研究方法

按研究目的归类文中使用的方法,便于定位所需技术。

产品清单

实验环节名称品牌货号
FUMA----
MAGMA----
TwoSampleMR----
MR-PRESSO----
MVMR----
RNA测序----

关键环节

汇总复现实验时建议重点确认的条件及原文阅读提示。

环节核对要点
GWAS表型选择
17种疼痛相关表型的病例/对照数量、定义(如疼痛部位、持续时间、药物使用)及排除标准(如重叠病例)。
阅读提示:Methods: Sample and phenotype selection; Table 1
PLEIO联合分析
分析软件参数(如MAF、HWE阈值、效应量转换),性染色体是否排除,以及样本重叠处理方法。
阅读提示:Methods: Joint-analysis of 17 pain-related UKBB GWASs
风险位点鉴定
LD参考面板(1000 Genomes Phase 3)、独立显著SNP和lead SNP的r²阈值、LD块合并距离。
阅读提示:Methods: Risk loci and gene mapping
基因定位与功能注释
基因映射方法(位置、eQTL、染色质互作)的具体定义和数据库版本(如PsychENCODE、GTEx)。
阅读提示:Methods: Risk loci and gene mapping
孟德尔随机化
工具变量选择标准(如P值、LD剪除)、来源GWAS(独立于UKBB)、F统计量阈值、MR方法(IVW、MR-Egger等)。
阅读提示:Methods: Genetic correlation and Mendelian Randomization