《Journal of Pain Research》：Dissecting the Genetic Basis of Low Back Pain Independent of BMI Through Genomic Structural Equation Modeling
编辑推荐：
这篇综述运用基因组结构方程模型（GWAS-by-subtraction）创新性地剥离了体质指数（BMI）的影响，系统揭示了与BMI无关的腰痛（Low Back Pain, LBP）遗传架构。研究鉴定出三个全基因组显著位点，并通过多组学（multi-omics）证据（如转录组关联研究/TWAS、蛋白质组关联研究/PWAS）共同指向CHST3等关键候选基因，强调了神经免疫调节、细胞外基质稳态等独立于肥胖的病理机制，为理解腰痛的异质性提供了新视角。
方法
本研究基于公开的全基因组关联研究（GWAS）汇总统计数据，采用基因组结构方程模型（Genomic Structural Equation Modeling, Genomic SEM）进行二次分析。BMI的汇总统计数据来自遗传调查人体测量特征（Genetic Investigation of ANthropometric Traits, GIANT）联盟，涵盖了约70万名欧洲血统个体。腰痛（Low Back Pain, LBP）数据则来源于芬兰基因组（FinnGen）项目第12版（R12），包含60,099名病例和440,249名对照，均通过国际疾病分类第10版（ICD-10）编码M54（背痛）进行定义。通过GWAS-by-subtraction方法，研究者旨在识别每个单核苷酸多态性（SNP）与腰痛之间独立于BMI的关联。该模型（结构方程模型示意如 所示）将腰痛和BMI的汇总统计量回归到两个潜在变量上：BMI变量和非BMI变量。该分析采用了Genomic-SEM R软件包，并以P<5×10?8为全基因组显著性阈值。
结果
GWAS-by-subtraction分析识别出三个与BMI无关的腰痛相关的独立全基因组显著位点：rs6916321位于B7NZA1位点，rs2596501靠近HLA-B位点，以及位于SEC24C的一个新位点rs4148946。在曼哈顿图 中可以清楚看到这些显著位点。
随后进行的精细定位（fine-mapping）使用了Sum of Single Effects（SuSiE）统计框架。精细定位曼哈顿图 展示了各SNP的后验包含概率（posterior inclusion probability, PIP）。其中，rs2596501的PIP达到0.9999，表明其极有可能是该位点的因果变异。
转录组关联研究（Transcriptome-Wide Association Study, TWAS）整合了基因型-组织表达（Genotype-Tissue Expression, GTEx）v8参考面板，发现了与BMI无关的腰痛存在显著关联的基因-组织对。曼哈顿式TWAS图 突出了包括RNF212、FOXP2、NSUN5和CHST3在内的显著基因。
预测性蛋白质组关联研究（proteome-wide association study, PWAS）则识别出其预测蛋白表达水平与BMI无关的腰痛相关的基因。PWAS的曼哈顿图 标注了包括CHST3在内的关键基因，其P值为4.49×10?5，FDR为0.0350。
共定位分析（colocalization analysis）通过fastENLOC进行。热图 0) across 49 GTEx tissues. Each cell represents a unique gene-tissue pair, with the color intensity corresponding to the GLCP value. A GLCP value greater than 0.5 (represented by a darker shade) indicates strong evidence of colocalization between a GWAS locus and an eQTL signal for the given gene-tissue pair. Genes with multiple tissues showing strong colocalization are highlighted."> 展示了在49个GTEx组织中15个基因的基因水平共定位概率（Gene-level colocalization posterior probability, GLCP）。其中CHST3在包括多个大脑区域和淋巴细胞的八个组织中显示出强共定位证据。
一项跨组学收敛性分析（cross-omics convergence analysis）比较了TWAS、PWAS和共定位分析鉴定的位点。如韦恩图 所示，大多数显著发现在单一分析框架内。重要的是，CHST3是唯一在三种分析方法中均被一致识别的基因。
基因和通路水平的富集分析（Gene- and Pathway-Level Enrichment）使用了多标记基因组注释分析（Multi-marker Analysis of Genomic Annotation, MAGMA）工具。分析结果以条形图展示 ，揭示了显著关联的基因（如APOE、TOMM40、DCC、NFU1等）以及富集的通路，包括血清素能信号传导、突触后传递、染色质组织和细胞应激反应。
细胞类型富集分析（Cell-Type Specific Enrichment）旨在确定遗传效应可能介导的特定细胞类型和基因。UMAP图 展示了在神经外胚层（如感觉神经元、雪旺细胞）等群体中的富集。条形图 则清晰地表明，在众多细胞类型中，“分化神经元”的富集显著性最高（P=0.0063）。为了识别分化神经元细胞类型中具有影响力的基因，研究者进一步进行了分析。散点图 可视化展示了这些基因，其中GNG3、NHLH1和RND2等因其高基因风险和细胞类型特异性而被确定为具有高度影响力的基因。
讨论
本研究通过应用GWAS-by-subtraction框架，剥离了BMI相关效应，揭示了与BMI无关的腰痛的遗传学基础，鉴定出包括SEC24C新位点在内的三个全基因组显著位点，并汇聚了多组学证据，将CHST3确立为一个关键的候选因果基因。CHST3基因编码软骨素6-磺基转移酶，调节软骨细胞外基质中糖胺聚糖的硫酸化。其突变可导致严重的脊柱发育不良，而常见变异则与椎间盘退变相关。该基因在多种分析中的一致性信号，突显了细胞外基质稳态在独立于肥胖的腰痛机制中的重要性。
此外，富集分析指向了主要组织相容性复合体区域（如BTN2A1、NFKBIL1、PSORS1C2）的免疫相关基因，支持了免疫失调在疼痛产生中的作用。同时，通路富集和细胞类型分析强调了神经过程的作用，如血清素能信号传导、突触后传递，以及分化神经元的富集，并识别出GNG3、NHLH1、RND2等具有影响力的基因，这表明神经元发育和兴奋性在独立于BMI的腰痛易感性中扮演着核心角色。TWAS和PWAS中发现的MIPOL1、NFU1等基因则指向了发育和线粒体代谢过程。
结论
总而言之，这项研究首次大规模地剖析了BMI无关的腰痛遗传风险因素，识别了新的风险位点，强调了CHST3的汇聚作用，并揭示了独立于BMI发挥作用的免疫、神经元和代谢机制。这些发现细化了腰痛的遗传结构，并为未来的机制研究指明了新方向。
