《Frontiers in Genetics》：Bridging bone and muscle: ontology-driven prioritization of candidate genes at the osteoporosis–sarcopenia interface by gene ontology reverse lookup
编辑推荐：
摘要专业翻译
引言： 骨肌少症（Osteosarcopenia）反映了骨质疏松症（Osteoporosis）与肌少症（Sarcopenia）的共存状态，但骨与肌肉稳态受损之间共享的分子过程仍未得到充分阐明。本研究应用一种透明、基于假设生成的注释框架，对骨质疏松
摘要专业翻译
引言： 骨肌少症（Osteosarcopenia）反映了骨质疏松症（Osteoporosis）与肌少症（Sarcopenia）的共存状态，但骨与肌肉稳态受损之间共享的分子过程仍未得到充分阐明。本研究应用一种透明、基于假设生成的注释框架，对骨质疏松—肌少症界面的候选基因进行优先级排序。
方法： 研究人员将GOReverseLookup应用于人工筛选的基因本体（Gene Ontology, GO）术语集，这些术语代表与骨质疏松症和肌少症相关的生物学机制。一组经筛选的367个GO术语作为输入，代表关键疾病机制。利用Fisher精确检验（Fisher's exact test）结合Benjamini–Hochberg校正，并采用预定义的直系同源（ortholog）传播设置，纳入智人（H. sapiens）、小鼠（M. musculus）、大鼠（R. norvegicus）、斑马鱼（D. rerio）和热带爪蟾（X. tropicalis）的注释，该分析流程对同时满足两个筛选状态显著性标准的基因进行优先级排序。优先候选基因随后通过靶向文献评估、DisGeNET疾病关联评分、Expression Atlas组织表达数据和KEGG通路分析进行情境化验证。
结果： 主要分析确定了37个候选基因，同时满足骨质疏松症和肌少症相关筛选状态的预定显著性标准。分析重点突出了部分研究较少的候选基因，包括RBCK1、LGALS9、ZBTB7A和RPS3，这些基因在选定的骨和肌肉相关组织数据集中具有可检测的表达。对已知骨和肌肉生物学介导因子（包括TNFSF11和STAT3）的恢复支持了该优先级排序框架的生物学合理性。通路特征分析主要提示炎症和固有免疫信号、受调控细胞死亡以及破骨细胞相关过程。
结论： GOReverseLookup生成了一个透明、基于注释的候选基因列表，用于研究骨质疏松症和肌少症共享机制。这些发现应被解释为假设生成的优先级排序结果，需经独立实验和临床验证。
论文解读文章
一、研究背景与问题
骨肌少症（Osteosarcopenia）是骨质疏松症与肌少症的共病综合征，在老年人群中日益受到关注。骨质疏松症以骨量减少和微结构退化为特征，导致骨骼脆性增加和骨折风险升高；肌少症则表现为骨骼肌质量和力量的进行性丧失，与身体功能障碍、生活质量下降和死亡率增加相关。尽管两者传统上被视为独立疾病，但越来越多的证据表明它们在发病机制和危险因素方面存在显著重叠，包括衰老、激素水平下降、营养缺乏和久坐行为等。骨与肌肉组织均起源于间充质干细胞，通过旁分泌、内分泌和机械刺激相互联系。肌因子（如肌生成抑制蛋白、胰岛素样生长因子1等）参与骨重塑，而骨因子（如骨钙素、转化生长因子-β等）影响肌肉代谢和收缩功能。慢性低度炎症同时加速破骨细胞生成和蛋白质分解，氧化应激和激素失调进一步放大这些病理过程。
然而，识别骨与肌肉共享生物学过程中的候选基因面临方法论挑战。全基因组关联研究（Genome-Wide Association Studies, GWAS）和转录组全关联研究（Transcriptome-Wide Association Studies, TWAS）虽能将遗传变异与复杂性状联系起来，但小效应信号和因果基因的解读仍具困难。表型组关联研究（Phenome-Wide Association Studies, PheWAS）易受表型错误分类和多重检验负担影响。近期基于表示学习的方法（如超图嵌入、蛋白质相互作用图神经网络和本体向量化）虽能预测疾病基因，但其预测缺乏生物学可审计性，且性能依赖于现有疾病—基因注释，对于骨肌少症这类新兴复合表型而言注释稀疏且有偏倚。因此，本研究采用GOReverseLookup这一反向富集分析工具，通过人工筛选基因本体（Gene Ontology, GO）术语定义表型，评估基因在用户定义状态中的统计过表达或欠表达情况，实现透明、可审计的候选基因优先级排序，适用于表型复杂、注释不完善且缺乏大型多组织数据集的早期假设生成阶段。
二、研究设计与方法概述
研究人员首先通过文献综述界定骨质疏松症和肌少症的关键病理生理机制，将每个机制组转化为代表性关键词，在QuickGO数据库中进行靶向人工检索，筛选出367个GO术语作为输入。随后使用GOReverseLookup（v1.0.56）进行反向富集分析，主要分析靶向智人，纳入斑马鱼、大鼠、小鼠和热带爪蟾的直系同源数据，使用间接注释（最多三级父项）并排除电子推断注释（IEA）。统计评估采用Fisher精确检验和Benjamini–Hochberg错误发现率校正，主要阈值为q < 5 × 10?8。候选基因进一步通过PubMed和Google Scholar靶向文献检索、DisGeNET数据库（版本25.1.1）基因—疾病关联评分、Expression Atlas组织表达数据（骨组织为颈椎椎体，肌肉组织包括臂肌、腿肌、骨骼肌等）以及ShinyGO v0.82的KEGG通路富集分析进行系统表征。
三、研究结果
3.1 候选基因的鉴定与分类
主要全跨物种分析中，群体集合包含30,768个GO术语，25,504个基因产物记录进入统计评分，其中2,461个具有足够合格注释。分析确定了37个候选基因，同时满足骨质疏松症和肌少症相关筛选状态的显著富集标准（q < 5 × 10?8）。其中14个基因此前未报道与任一疾病相关，9个基因此前文献报道与骨质疏松症和肌少症均相关且具有非零DisGeNET评分，14个基因仅与其中一种疾病相关或在一个或两个组织中表达低于阈值。在14个未报道基因中，RBCK1、TERF2IP、LGALS9、ZBTB7A和RPS3五个基因在骨和肌肉中均显示中等表达且无疾病关联评分（GDA≈0）。在9个已知基因中，TNFSF11和STAT3已被直接涉及骨肌少症。
3.2 显著性阈值和直系同源注释的敏感性分析
全跨物种模型在q < 5 × 10?8下确定37个候选基因，q < 0.01下为208个，q < 0.05下为200个。q < 0.05候选数低于q < 0.01反映了双向优先级排序规则，即基因若与负向状态显著相关可能被排除。仅人类注释模型在任何阈值下均无候选基因。在主要阈值下，加入小鼠和大鼠直系同源注释后获得31个候选基因，其中28个保留自主要37基因集合。10个研究较少的候选基因在哺乳动物扩展模型中保留，BCL2L11、BTK、IL18R1和RBCK1四个未被保留。
3.3 KEGG通路分析
37个优先候选基因的KEGG富集分析确定了15个显著过表达通路（q < 5 × 10?5，折叠富集>25）。非感染相关通路包括NF-κB信号通路（FE=57.9）、胞质DNA传感通路（FE=47.8）、NOD样受体信号通路（FE=39.5）、炎性肠病（FE=37.1）、催乳素信号通路（FE=34.4）、C型凝集素受体信号通路（FE=28.9）、破骨细胞分化（FE=28.7）、TNF信号通路（FE=26.9）、坏死性凋亡（FE=26.5）以及脂质和动脉粥样硬化（FE=25.3）。排除IL1B、IL18、MYD88、PYCARD、RIPK1、TRAF2、NOD2和BTK八个免疫相关基因后，无通路满足预定组合阈值。
四、讨论与结论总结
讨论部分指出，本研究是反向基因本体分析在复杂肌肉骨骼表型中的首次综合应用。研究恢复了已知骨和肌肉生物学相关基因，支持了优先级排序框架的生物学合理性和表面效度。研究较少的候选基因具有特殊价值：RBCK1编码HOIL-1，是线性泛素链组装复合物（LUBAC）的催化组分，调节NF-κB通路、凋亡和坏死性凋亡，其突变可导致多聚葡聚糖体肌病1型（PGBM1）；LGALS9编码半乳糖凝集素-9，通过TIM-3受体抑制破骨细胞生成并促进成骨细胞分化，还通过结合NLRP3诱导其自噬降解来抑制炎症小体；TERF2IP编码Rap1，是端粒庇护蛋白复合物组分，同时作为IκB激酶复合物中的衔接蛋白放大NF-κB信号，端粒功能障碍可触发DNA损伤反应和衰老相关分泌表型（SASP）；ZBTB7A编码锌指转录因子LRF，在破骨细胞谱系中发挥阶段特异性效应，在骨骼肌中受miR-182-5p抑制可促进肌源性分化；RPS3编码核糖体蛋白S3，具有多种核糖体外功能，包括作为NF-κB复合物的非经典亚基调控靶基因转录。
KEGG通路富集揭示了三个相互关联的机制主题：炎症信号与破骨细胞调控的整合、固有免疫传感与炎症小体相关细胞因子信号、以及炎症应激相关的受调控细胞死亡信号。坏死性凋亡通路的富集提示受调控细胞死亡机制可能构成骨和肌肉同时退化的更广泛炎症组织应激反应的一部分。
研究的方法学优势在于GOReverseLookup允许研究人员通过筛选GO术语明确定义感兴趣状态，产生透明、可审计和生物学可解释的基因列表，且对复合或多因素表型表现良好。局限性包括：文献来源的GO术语可能引入循环性；直系同源注释扩展可能增加假阳性风险；GO注释完整性限制可能导致假阴性；骨组织表达数据集有限，颈椎椎体样本为异质细胞混合物；KEGG富集受广泛注释的固有免疫和炎症基因驱动；缺乏实验验证和独立患者队列或外部遗传资源验证；未进行性别分层分析；无法区分因果驱动因子与关联生物标志物或炎症和退化的下游后果。研究结论为：GOReverseLookup生成了一个透明、基于注释的候选基因列表，用于研究骨质疏松症和肌少症共享机制，这些发现应被解释为假设生成的优先级排序结果，需经独立实验和临床验证。
