《International Journal of Eating Disorders》：Integrative Genomic and Functional Investigation of the Multi-Layered Genetic Architecture Between Anorexia Nervosa and Bone Loss
编辑推荐：
**目的** 骨丢失是神经性厌食症（anorexia nervosa, AN）的一种严重且常不可逆的并发症，然而这种共病背后的遗传机制仍未得到充分探索。本研究旨在构建AN与估算跟骨骨矿物质密度（estimated calcaneal bone mineral
**目的** 骨丢失是神经性厌食症（anorexia nervosa, AN）的一种严重且常不可逆的并发症，然而这种共病背后的遗传机制仍未得到充分探索。本研究旨在构建AN与估算跟骨骨矿物质密度（estimated calcaneal bone mineral density, eBMD）之间的综合性遗传架构。
**方法** 研究人员应用了一个整合性框架，在单变异、多变异及基因表达水平上纳入遗传相关性、多效性关联及因果推断分析。通过体外（in vitro）实验进行功能验证，以探究关键候选基因的生物学作用。
**结果** 尽管缺乏显著的全局遗传相关性，局部遗传相关性分析仍鉴定了8p21.2和10q26.3处的显著信号。孟德尔随机化（Mendelian randomization, MR）分析提示，遗传易感的AN对eBMD具有提示性的负向因果效应。研究检测到广泛的多效性信号，尤其是在3p21.31和10q26.3位点，这些位点富含与两种性状均相关的基因。值得注意的是，研究人员在11q23.2的NCAM1附近识别出一个新的多效性信号，该信号得到了多层遗传证据的支持，并通过体外功能实验得以证实。NCAM1是一个公认的神经相关基因，当其在破骨前体细胞中过表达时，可促进破骨细胞分化和骨吸收，表明NCAM1在神经和骨骼组织中具有各自独特的功能作用。
**讨论** 本研究构建了AN与eBMD之间潜在的综合性遗传架构，并强调NCAM1是一个关键的多效性基因。
**神经性厌食症与骨丢失共享遗传架构的整合性研究解读**
**一、研究背景与问题提出**
骨密度（bone mineral density, BMD）的峰值主要在青春期和成年早期获得，是终身骨骼健康的关键决定因素。然而，神经性厌食症（anorexia nervosa, AN）在此关键发育时期的起病，常导致骨矿物质密度的显著下降，进而增加个体终生的骨质疏松和骨折风险。一项纳入43项研究的荟萃分析显示，AN患者在全身、脊柱、股骨和髋部等多个部位的BMD均显著低于健康对照。更为严峻的是，即使经过临床康复和体重恢复正常，BMD往往仍无法恢复至年龄匹配的正常水平，使其成为该疾病一种主要的不可逆并发症。纵向研究亦证实，有AN病史的个体在体重恢复后仍存在持续的骨量和骨质量缺陷。
为理解这一持续性骨骼损害的机制，既往研究多聚焦于下游的生理改变，如身体组分变化和内分泌系统失调，特别是下丘脑-垂体轴或激素平衡的紊乱。尽管这些研究提供了关于生理后果和潜在治疗靶点的重要见解，但它们主要处理的是继发性效应，而非上游的内在病因学机制。鉴于AN的高度多基因性和多器官受累特征，研究人员假设共享的遗传变异可能通过多效性效应同时影响神经和骨骼稳态。为验证此假设，需要足够的统计功效来检测微弱的共享遗传效应，因此选择目前骨性状中可用数据集最大、且与DXA来源BMD具有显著遗传一致性并对骨折预测具有临床相关性的估算跟骨骨矿物质密度（eBMD）作为主要骨骼表型。
大规模全基因组关联研究（genome-wide association studies, GWAS）已发现多个与各性状相关的遗传变异，为探索潜在共享遗传决定因素奠定了基础。然而，由于严格的基因组显著性阈值、显著的 phenotype 异质性以及技术限制，许多此类位点仍未被检测到。此外，多数GWAS鉴定的变异位于非编码或基因间区域，使其生物学解释变得复杂。为此，新兴的后GWAS分析方法，如跨性状荟萃分析、多变异检验和转录组wide关联研究（transcriptome-wide association study, TWAS），通过将遗传信号与基因水平功能和表达联系起来，为生物学机制提供了补充见解。
**二、研究内容与主要结论**
本研究旨在整合多种分析方法，阐明AN与eBMD之间的遗传架构和潜在多效性机制。研究有两个主要目标：(1) 通过在单变异、多变异和基因表达水平上整合分析，构建AN与eBMD之间全面的遗传架构，实现对遗传相关性、多效性关联和因果推断的系统评估；(2) 通过结合多层统计分析与体外验证的整合框架，优先筛选促成AN和骨丢失共享易感性的基因组位点和功能基因。
研究结果显示，AN和eBMD在基因组wide水平上无显著遗传相关性，但局部遗传相关性分析鉴定了8p21.2和10q26.3两个显著信号区域。孟德尔随机化分析提示AN对eBMD存在提示性的负向因果效应。跨性状分析检测到广泛的多效性信号，尤其集中于3p21.31和10q26.3区域。最值得注意的是，研究在11q23.2的NCAM1附近识别出一个新的多效性信号，并经体外功能实验证实：在破骨前体细胞中过表达NCAM1可显著促进破骨细胞分化和骨吸收功能。
**三、关键技术与方法**
本研究应用了多项关键分析技术。遗传相关性分析采用交叉性状连锁不平衡得分回归（LDSC）评估全局遗传相关，使用SUPERGNOVA进行局部遗传相关分析。因果推断采用双样本孟德尔随机化（MR）方法，以逆方差加权法为主分析，辅以加权中位数法和CAUSE进行敏感性分析。多效性位点识别运用了三种互补的跨性状meta分析：MTAG、CPASSOC和PLEIO。基因水平分析采用MAGMA进行基因关联分析，并用FUSION进行TWAS和SMR分析，整合eQTL数据。功能验证使用RAW264.7小鼠巨噬细胞系，通过质粒转染构建NCAM1过表达模型诱导破骨细胞分化。
主要数据来源包括：AN GWAS数据来自精神病学基因组学联盟进食障碍工作组（PGC-ED），包含16,992例病例和55,525例欧洲裔对照；eBMD GWAS数据来自英国生物银行（UKB）队列，涉及426,824名参与者。
**四、研究结果**
**4.1 全局与局部遗传相关性分析**：LDSC分析显示AN与eBMD之间无显著全局遗传相关性（rg=-0.01，p=0.56）。但SUPERGNOVA局部遗传相关分析鉴定出8p21.2和10q26.3两个显著信号区域。其中8p21.2位点已在既往eBMD GWAS中被报道，10q26.3区域包含与AN和eBMD均相关的MGMT基因。
**4.2 孟德尔随机化分析**：所有工具变量的F统计量均超过10，提示无弱工具变量偏倚。IVW方法未发现AN对eBMD的显著因果效应，但加权中位数法和排除回文SNP的分析均提示了微弱的负向因果效应。其他骨骼部位如股骨颈、腰椎和前臂BMD未观察到显著因果效应。反向MR分析显示eBMD对AN风险无因果影响。
**4.3 跨性状meta分析**：CPASSOC检测到7个全基因组显著多效性SNP，PLEIO识别出6个，两种方法有3个SNP重叠。排名前三位点均位于3p21.31和10q26.3。最显著的多效性SNP为位于3p21.31的rs62259939，10q26.3的rs12254582靠近MGMT。关键发现是位于11q23.2的rs11214474，这是唯一显著的“新”共享SNP，映射到NCAM1基因附近。
**4.4 基因水平验证**：MAGMA分析显示NCAM1在AN和eBMD中均达到Bonferroni校正显著性。TWAS分析发现NCAM1表达与eBMD在肾上腺和胫骨神经中显著相关。SMR分析基于BrainMeta v2 eQTL数据进一步确认了NCAM1表达与eBMD的关联。
**4.5 体外功能验证**：通过质粒转染构建NCAM1过表达的RAW264.7细胞，并诱导破骨细胞分化。TRAP染色显示NCAM1过表达增加了TRAP阳性破骨细胞数量和平均面积。F-肌动蛋白染色显示NCAM1过表达组的肌动蛋白环更大且核数量更多。骨吸收陷窝实验表明NCAM1过表达显著增加了骨吸收陷窝数。RT-qPCR结果显示NCAM1过表达上调了NFATc1、TRAP、Cathepsin K和MMP-9等破骨细胞分化相关基因的mRNA表达，同时下调了β-Catenin的表达。
**五、讨论与结论**
讨论部分指出，虽然AN和eBMD未显示显著的全局遗传相关性，但这一现象可能归因于两种疾病的高度多基因性和异质性，不同基因组区域的共享遗传效应在聚合成单一全基因组估计时可能相互抵消。此外，AN GWAS有限的样本量可能限制了检测微弱全基因组相关的统计功效。有关因果推断，本研究结果未提供AN对BMD因果影响的决定性证据，因此关联应被解读为提示性而非结论性。值得注意的是，10q26.3包含与AN和eBMD均相关的MGMT基因，3p21.31位点含有PRKAR2A基因，这两种基因此前均被发现与AN和骨发育相关。NCAM1与既往已确立的AN风险基因CADM1具有功能和基因组相似性，二者均编码免疫球蛋白超家族的细胞粘附分子，在突触形成和可塑性中起关键作用。在骨骼系统中，NCAM1在成骨细胞和破骨细胞分化中均发挥作用，本研究发现其过表达通过下调β-Catenin促进破骨细胞生成。研究还从公共卫生角度探讨了AN通常于青春期起病导致不可逆BMD降低，使骨质疏松和骨折风险显著增加，而目前的干预措施效果和安全性仍不确定。研究局限包括未探讨性别特异性效应、数据主要来自欧洲人群、未考虑年龄依赖性遗传效应及进一步机制实验的必要性。
研究结论指出，通过整合多层基因组分析与功能验证，研究人员鉴定了AN和eBMD之间的共享遗传架构，并突出NCAM1为关键多效性效应基因。这些发现为未来的机制研究提供了基础，并可能为探索解决AN个体骨骼并发症的策略提供后续研究参考。
