《Nature Medicine》：The genetic architecture of fibromyalgia across 2.5 million individuals
编辑推荐：
纤维肌痛是一种常见且使人衰弱的慢性疼痛综合征，其病因尚不清楚。在此，研究人员开展了一项多祖先全基因组关联研究（GWAS）荟萃分析，涵盖来自11个队列的2,563,755名个体（54,629例病例和2,509,126例对照），确定了26个纤维肌痛风险位点。最强的
纤维肌痛是一种常见且使人衰弱的慢性疼痛综合征，其病因尚不清楚。在此，研究人员开展了一项多祖先全基因组关联研究（GWAS）荟萃分析，涵盖来自11个队列的2,563,755名个体（54,629例病例和2,509,126例对照），确定了26个纤维肌痛风险位点。最强的关联来自亨廷顿病（Huntington's disease, HD）致病基因HTT中的一个编码变异。基因优先化分析提示了HTT调节因子GPR52，以及多种具有神经功能的相关基因，包括DCC、DRD2/NCAM1、MDGA2和CELF4。纤维肌痛的遗传力（heritability）仅富集于脑组织和神经细胞类型。纤维肌痛与多种慢性疼痛、精神疾病和躯体疾病表现出强正遗传相关性，包括与下背痛、创伤后应激障碍（post-traumatic stress disorder, PTSD）和肠易激综合征（irritable bowel syndrome, IBS）的遗传相关性高于0.7。尽管纤维肌痛患病率存在显著的性别差异，但男性和女性之间的纤维肌痛遗传结构几乎相同。该研究提供了强有力的遗传证据，将纤维肌痛定义为一类中枢神经系统疾病，从而为其复杂的病理生理学和广泛的临床共病建立了生物学框架。
**研究背景与目的**：纤维肌痛是一种常见且致残的慢性疼痛综合征，以广泛性肌肉骨骼疼痛、疲劳、睡眠障碍、认知障碍和躯体症状为特征，常与多种疼痛状态、肠易激综合征、肌痛性脑脊髓炎/慢性疲劳综合征、自身免疫性和神经精神性疾病以及代谢综合征共病。尽管其被认为具有部分遗传基础，但此前寻找特异性遗传位点的尝试一直未能取得明确结果。该研究旨在通过大规模遗传分析，鉴定纤维肌痛的风险位点、优先候选致病基因、定位遗传力富集的组织与细胞类型、量化与共病的遗传重叠，并评估性别间遗传一致性，从而为这一长期争论的疾病建立坚实的生物学基础。该论文发表在《Nature Medicine》。
**主要关键技术与方法**：研究人员从11个队列（包括All of Us、UK Biobank、FinnGen、Estonian Biobank、Genes & Health、Mass General Brigham Biobank、Michigan Genomics Initiative、deCODE Genetics/Amgen Iceland、Copenhagen Hospital Biobank、Danish Blood Donor Study、Intermountain Health和Nashville Biosciences）收集数据，进行多祖先全基因组关联研究（GWAS）荟萃分析。采用固定效应逆方差加权荟萃分析，通过连锁不平衡（LD）评分回归截距校正膨胀，并利用LD评分回归（LDSC）进行遗传力估计、组织与细胞类型富集分析（LDSC-SEG）以及遗传相关性分析。使用FLAMES机器学习方法和deCODE遗传注释管道（包括表达数量性状位点（eQTL）和蛋白质数量性状位点（pQTL）分析）进行候选致病基因优先化。基于PRS-CS构建多基因风险评分（PRS）。
**研究结果**：
- **与26个变异的关联确立了纤维肌痛的遗传基础**：通过多祖先GWAS荟萃分析（54,629例病例，2,509,126例对照），研究人员发现26个独立全基因组显著（P<5×10
-8）风险位点，效应大小在队列间基本一致。性别分层分析显示，尽管女性病例占绝大多数（87.7%），但男女间遗传相关性为1.03，几乎完全一致，表明纤维肌痛的遗传结构在性别间高度相似。
- **纤维肌痛最强关联为亨廷顿基因中的编码变异**：最显著关联为HTT基因中的一个框内谷氨酸缺失变异（rs149109767-A，OR=1.09，P=2.2×10
-12）。该变异与亨廷顿病（HD）的关联在条件分析后消失，提示其与HD致病重复扩增不同。另一个可能关联涉及HD药物靶点GPR52。
- **其他关联涉及疼痛处理通路**：包括CELF4（调节突触形成和伤害感受器兴奋性）、STK31/NPY（神经肽Y参与疼痛调节）、ARHGAP15/KYNU（犬尿氨酸通路参与疼痛处理）等基因。
- **其他显著关联指向脑相关机制**：包括DRD2/NCAM1（多巴胺D2受体和神经细胞黏附分子）、CAMKV/MST1R（钙/钙调蛋白依赖性蛋白激酶和酪氨酸激酶受体）、MAML3（Notch信号通路共激活因子）、MDGA2（突触组织和轴突迁移）、DCC（轴突导向受体）、SRD5A2/SPAST（雄激素代谢和微管调节）、BCL11A（神经发育）、BPTF（染色质重塑）、PPP2R2B/GPR151（磷酸酶和孤儿G蛋白偶联受体）以及NPC1（胆固醇运输和自噬）等基因。
- **先导变异的表型组关联揭示了显著的多效性**：通过GWAS Catalog和MVP-Finngen-UKBB表型组关联研究（PheWAS），发现26个先导变异与教育程度、体重指数、2型糖尿病、疼痛强度、失眠、抑郁症状、精神分裂症等广泛性状相关，且与下背痛、胃食管反流病、镇痛药和抗炎药处方等关联最频繁。
- **纤维肌痛遗传力仅富集于神经组织和细胞类型**：使用LDSC-SEG分析，发现纤维肌痛遗传力在53个GTEx组织中的5个大脑区域（皮质、尾状核、额叶皮质、壳核和前扣带皮层）及119个小鼠单细胞类型中的13个细胞类型（主要为神经元，包括齿状回神经元、肠神经元和多种中间神经元）显著富集，且唯一显著富集的细胞谱系为神经谱系。
- **与共病障碍的普遍遗传相关性**：通过LDSC计算纤维肌痛与FinnGen中855种疾病的遗传相关性，发现337个Bonferroni显著相关。最强相关性见于疼痛和肌肉骨骼疾病（如颈臂综合征rg=0.87、肌痛rg=0.82）、精神行为障碍（如创伤后应激障碍rg=0.78、转换障碍rg=0.63、抑郁症rg=0.63）、消化系统疾病（如肠易激综合征rg=0.70、功能性消化不良rg=0.67）以及泌尿生殖系统疾病（如多囊卵巢综合征rg=0.59）。与自身免疫性疾病的遗传相关性相对较弱，且与哮喘的T2低亚型（rg=0.44）和类风湿关节炎的血清阴性亚型（rg=0.41）相关性更强，进一步支持纤维肌痛与典型外周驱动炎症性疾病的遗传结构差异。
- **遗传风险预测在患者分层方面效果有限**：基于多祖先留出UK Biobank的荟萃分析构建的PRS，在UK Biobank欧洲裔参与者中表现最佳（AUC=0.59），南亚裔（AUC=0.55）和非洲裔（AUC=0.55）准确性较低。尽管AUC中等，但PRS有效分层风险：欧洲裔最高PRS五分位数的纤维肌痛比值比为1.5（95%CI 1.4-1.7），最低五分位数为0.63（95%CI 0.57-0.68），患病率从最低五分位数的约1.0%增至最高五分位数的2.4%。
**总结与讨论部分**：该研究绘制了纤维肌痛的遗传结构，鉴定出26个风险位点，并为纤维肌痛主要是一种中枢神经系统疾病的概念提供了强有力的遗传验证。研究结果支持纤维肌痛的中枢敏化模型，并指出HTT、GPR52和CELF4等位点作为潜在治疗靶点。尽管存在样本主要为欧洲裔、依赖ICD编码等局限性，但该研究为纤维肌痛的病理生理学及其与共病共享病因的假设驱动研究提供了具体的分子起点。研究结论部分原文翻译如下：“我们的研究绘制了纤维肌痛的遗传结构，鉴定出26个风险位点，并为纤维肌痛主要是一种中枢神经系统疾病的概念提供了强有力的遗传验证。鉴定特异性风险位点为该领域提供了具体的分子起点，使得能够开展假设驱动的病理生理学及与共病共享病因的研究。”
