整合基因组结构方程模型解析胰岛素抵抗相关代谢性状的共享遗传结构及候选基因

《PLOS One》：Integrative genomic SEM analysis of shared genetic architecture and candidate genes for insulin resistance–related metabolic traits

　　

胰岛素抵抗（insulin resistance, IR）是一种复杂代谢状态，涉及胰岛素敏感性受损、代偿性胰岛素分泌、脂质代谢异常及脂肪分布改变。然而，这些相关表型背后的共享遗传结构仍未完全明确。本研究应用基因组结构方程模型（Genomic Structural Equation Modeling, Genomic SEM）对空腹胰岛素（fasting insulin, FI）、胰岛素敏感性指数（insulin sensitivity index, ISI）、甘油三酯（triglycerides, TG）和腰臀比（waist-to-hip ratio, WHR）的全基因组关联研究（genome-wide association study, GWAS）汇总统计量进行分析，构建IR相关潜在遗传因子；随后对该潜在因子开展单核苷酸多态性（single-nucleotide polymorphism, SNP）水平关联分析，并进行风险位点注释、精细定位、转录组范围关联分析（transcriptome-wide association analysis, TWAS）、基因水平优先排序、富集分析及对所选错义替换的探索性结构生物学评估。单因子模型显示近似拟合可接受但非最优，并呈现生物学一致的遗传模式：FI、TG和WHR载荷为正，ISI载荷为负。该潜在因子具有可检测的SNP遗传力，并揭示3,067个全基因组显著候选SNP，同时将Q_SNP保留为补充性SNP水平异质性诊断指标。功能注释鉴定出100个基因组风险位点及MAGMA中1,530个错误发现率（false discovery rate, FDR）显著基因。精细定位优先排序多个候选区域，包括PPARG、RSPO3、SLC38A11、SLC39A8、BCL2、GDF5和SPATA5附近位点。TWAS与FOCUS进一步突出参与脂肪细胞生物学、糖脂代谢、胰岛素分泌和表达调控的候选基因，包括PPARG、MLXIPL、TCF7L2、FAM13A、CCDC92和SLC22A3。富集分析提示，与IR相关潜在因子有关的遗传信号主要与脂肪分布、心脏代谢性状、脂肪组织、间充质细胞谱系及活跃调控基因组区域相关。在一项单独的探索性分析中，AlphaFold3建模、ThermoMPNN预测和分子动力学模拟提示，提名候选蛋白中的三个超罕见错义替换存在结构扰动：BCL2-W188G、GDF5-C400Y和RSPO3-C129Y。这些替换依据功能影响注释和建模可行性选择，并非GWAS lead变异或精细定位因果变异，且不能确立GWAS关联背后的机制。这些发现为理解IR相关代谢异常的共享机制提供了整合遗传与结构的框架，同时强调仍需进一步实验验证。

研究背景：胰岛素抵抗（insulin resistance, IR）指靶组织对胰岛素生物学作用反应性降低，主要表现为骨骼肌葡萄糖摄取受损、肝脏葡萄糖生成增加、脂肪组织代谢调控失调及代偿性高胰岛素血症；其发生发展与脂肪分布异常、脂质代谢紊乱、慢性低度炎症、氧化应激和胰岛素信号受损相关。IR是2型糖尿病、肥胖、代谢综合征、非酒精性脂肪性肝病及心血管疾病风险增加的重要共同病理基础，常涉及多器官系统。尽管临床评估、机制研究和治疗管理已有进展，IR的遗传基础及相关表型共享的生物学机制仍不完全清楚；单一临床指标不能充分捕获IR，因此研究人员开展本研究，应用基因组结构方程模型（Genomic Structural Equation Modeling, Genomic SEM）整合IR相关表型的GWAS汇总统计量，构建IR相关潜在遗传因子，并连接遗传统计发现与结构生物学解释。研究人员基于空腹胰岛素（fasting insulin, FI）、胰岛素敏感性指数（insulin sensitivity index, ISI）、甘油三酯（triglycerides, TG）和腰臀比（waist-to-hip ratio, WHR）构建潜在因子，开展SNP水平关联、风险位点注释、精细定位、转录组整合、富集、细胞类型与功能注释、多基因风险评分及候选错义替换的探索性结构分析。结论提示，共享IR相关遗传信号与脂肪分布、脂质代谢、脂肪组织功能及调控基因组区域密切相关；精细定位和TWAS/FOCUS优先排序出参与脂肪细胞生物学、代谢调控和表达控制的候选位点与基因；单独探索性结构分析提示所选超罕见错义替换可能改变蛋白稳定性与构象动力学，但这些替换并非GWAS lead或精细定位因果变异，不能解释或验证常见变异GWAS关联。该研究的意义在于提供从遗传统计发现到结构生物学解释的整合分析框架，为后续功能验证和优先候选基因的机制研究提供理论基础。

关键技术方法（≤250字）：研究人员基于四项欧洲血统公开GWAS汇总统计量（FI、ISI、TG、WHR；原始研究已获伦理批准和知情同意），采用多变量连锁不平衡得分回归（linkage disequilibrium score regression, LDSC）估计遗传协方差，并以Genomic SEM构建单因子潜在遗传模型及SNP水平关联；随后用FUMA/MAGMA进行风险位点与基因水平注释，SuSiE和FINEMAP进行精细定位，FUSION结合GTEx v8 sCCA模型开展TWAS并用FOCUS进行基因水平贝叶斯精细定位；进一步采用CELLECT、分区LDSC和PRS-CS分析细胞类型、功能注释与多基因风险；最后对候选错义替换进行AlphaFold3建模、ThermoMPNN ΔΔG预测及GROMACS 100 ns分子动力学（molecular dynamics, MD）模拟。

研究结果：
3.1 基因组结构方程模型潜在因子构建与模型拟合：研究人员通过多变量LDSC和单因子Genomic SEM分析发现，四个输入GWAS均具有可检测SNP遗传力，表型间存在遗传相关；模型拟合指标为χ²=23.486、df=2、P=7.94×10^?6、AIC=39.486、CFI=0.9399、SRMR=0.0494，属可接受但非最优近似拟合；标准化载荷为FI 0.924、TG 0.467、WHR 0.557、ISI ?0.822，方向与IR生物学模式一致。
3.2 SNP水平基因组结构方程模型GWAS与基因组控制：通过SNP水平Genomic SEM GWAS发现，潜在因子平均χ²为1.6344、λ_GC为1.4034、LDSC截距为0.9755，未见明显残余膨胀；观察尺度SNP遗传力h²=0.0599，以P<5×10^?8为阈值鉴定3,067个候选SNP，P<5×10^?12和P<5×10^?16分别鉴定501个和9个信号；Q_SNP作为补充性SNP水平诊断单独考察。
3.3 Lead SNP、风险位点与基因水平结果：通过FUMA和MAGMA分析鉴定100个基因组风险位点、1,530个FDR<0.05候选基因及152个lead SNP，后者多位于非编码区；与四个单表型GWAS显著位点进行排除式比较，进一步识别142个可能主要由共享遗传结构驱动的候选SNP。
3.4 精细定位分析：通过SuSiE和FINEMAP对主要风险区域进行精细定位，发现多个区域存在后验包含概率（posterior inclusion probability, PIP）较高且两方法支持一致或可比的重要候选变异，代表位点涉及AC006458.3、DNAH9、PPARG、RNA5SP111、RSPO3、SLC38A11和SLC39A8附近区域；BCL2、BCL2L11、MECOM、GDF5、MAFB、SPATA5及HOXC6/MIR196A2附近区域亦获不同程度后验概率支持。
3.5 TWAS/FOCUS候选基因优先排序：通过FUSION TWAS和FOCUS贝叶斯基因水平精细定位发现402个Bonferroni显著基因-组织信号，以TWAS显著性和FOCUS可信集联合筛选得到221个重叠基因-面板信号、147个唯一候选基因；CCDC92、SNX10、MLXIPL、JAZF1和PPARG预测表达与潜在IR因子正相关，BHLHE41、FAM13A和TCF7L2呈负相关，SLC22A3和VEGFB在不同sCCA面板方向不一致。
3.6 通路、细胞类型与孟德尔疾病基因富集：基于MAGMA Bonferroni显著基因开展基因集、通路和表型本体富集，发现候选基因主要富集于腰臀比、BMI校正腰臀比、腰围、臀围和体型指数等脂肪分布性状，以及甘油三酯、HDL胆固醇、空腹胰岛素、2型糖尿病、代谢综合征、血压和冠状动脉疾病等心脏代谢表型；GO和GO-slim分析涉及碳水化合物代谢、蛋白糖基化、受体配体活性、囊泡介导运输、内体和高尔基运输、内质网功能及染色质相关过程；细胞和组织富集指向脂肪组织及间充质细胞谱系，单细胞分析中肢体肌肉间充质干细胞、气管间充质细胞和脂肪来源间充质干细胞通过FDR校正，皮下和内脏脂肪组织及脂肪核中H3K4me1、H3K9ac、H3K27ac等活跃调控标记显著富集；MendelVar未发现经FDR校正的疾病本体或人类表型本体条目。
3.7 基因组功能区域遗传力贡献：通过分区LDSC分析发现，潜在IR因子的SNP遗传力主要富集于保守元件、DHS/DGF开放染色质、增强子、超级增强子、转录因子结合位点及转录起始位点附近区域，H3K27ac、H3K4me1、H3K4me3和H3K9ac等活跃调控标记亦显著富集；编码区、扩展编码区、内含子、扩展启动子区和UTR区呈不同程度富集，而抑制性染色质注释富集值低于1。
3.8 染色体分层遗传信号分布：通过染色体分层汇总PRS分析发现，IR相关遗传信号在染色体间分布不均，19、22、16、12和4号染色体呈主要正贡献，其中19号染色体贡献最大；5、18、13、20、9和6号染色体呈负贡献；该结果反映聚合遗传趋势，不能直接解释为特定因果位点或功能基因组区域。
3.9 代表性错义替换的结构分析选择：结合精细定位、位点注释和整合遗传证据提名候选基因，并依据功能预测注释、AlphaFold置信度、氨基酸替换类型、蛋白大小、结构可解释性及建模可行性，从BCL2、GDF5和RSPO3中各选择1个代表性替换进入探索性结构分析，即BCL2-W188G、GDF5-C400Y和RSPO3-C129Y；其Ensembl注释CADD分别为32、31、32，REVEL分别为0.756、0.986、0.816；在gnomAD v4.1.1中BCL2-W188G和GDF5-C400Y未检出，RSPO3-C129Y为单例，合并等位基因频率6.203×10^?7；这些超罕见替换并非GWAS lead变异或精细定位因果变异。
3.10 结构建模分析：通过AlphaFold3预测发现，BCL2、GDF5和RSPO3的突变位点及邻近区域处于可接受至较高pLDDT置信范围，提示预测结构适合后续比较性结构评估和探索性MD模拟。
3.11 分子动力学分析：通过100 ns MD模拟发现，BCL2-W188G在平衡阶段RMSD高于野生型，Rg和SASA升高，远端残基RMSF增加，提示整体构象偏差增大、紧凑性降低、溶剂暴露增加及远端灵活性增强；GDF5-C400Y的RMSD未高于野生型，但Rg持续升高，模拟中后期SASA增加，局部及远端RMSF呈区域重分布；RSPO3-C129Y的RMSD和SASA明显升高，突变位点及多个远端区域RMSF增强，Rg后期未升高，提示构象重排或局部重堆积而非单纯整体松散。
3.12 蛋白ΔΔG预测：通过ThermoMPNN预测发现，BCL2-W188G、GDF5-C400Y和RSPO3-C129Y在突变敏感性热图中呈不同程度红色信号，提示这些位点替换可能增加蛋白自由能差；GDF5-C400Y和RSPO3-C129Y红色信号更明显，提示对蛋白稳定性影响可能更强，BCL2-W188G为轻度至中度结构失稳倾向。

讨论部分总结：研究人员将潜在因子解释为IR相关共享代谢遗传表型，而非IR临床诊断变量；其载荷方向符合IR预期，支持IR相关代谢异常具有多基因基础及连接脂肪分布、血脂异常、代偿性胰岛素分泌和胰岛素敏感性降低的共享遗传结构。阈值分层信号涉及脂肪分布、脂质代谢和血管代谢调节，其中6p21.1区域应解释为与脂肪分布、脂质相关表型和代谢风险遗传相关的区域，不能据此确立VEGFA或POLR1C通路；其他信号功能注释有限，主要提示Genomic SEM可识别单表型GWAS不易突出的跨表型共享遗传结构候选区域。精细定位和转录组整合提示PPARG、RSPO3、SLC39A8等位点涉及脂肪组织功能、体脂分布和代谢应激调节，PPARG、MLXIPL、TCF7L2、FAM13A和CCDC92具有较明确代谢可解释性，VEGFB和SLC22A3因方向不一致应谨慎解释；TWAS/FOCUS仅反映遗传预测表达与潜在因子的关联，不能单独确立表达变化的因果作用。富集分析支持脂肪组织及调控基因组区域是共享IR遗传信号的重要背景，分区LDSC提示遗传力集中于保守、开放、增强子、超级增强子、转录因子结合位点和转录起始位点附近区域；染色体分层结果支持共享IR相关代谢遗传易感性的多基因特征。探索性结构分析提示三个超罕见错义替换可能造成不同程度构象动力学改变，BCL2-W188G可能削弱局部疏水堆积和空间支撑，GDF5-C400Y可能影响构象紧凑性和区域动态分布，RSPO3-C129Y可能影响局部结构稳定性及远端动态性质并诱导构象重排；ThermoMPNN ΔΔG预测与MD结果总体方向互补，但均属计算预测，需蛋白相互作用实验、细胞功能实验和体内模型验证。局限包括：输入GWAS主要来自欧洲血统人群，外推需验证；潜在因子不能替代IR临床诊断或用于个体风险评估；新位点识别和精细定位受LD参考面板、区域定义和统计模型影响；TWAS/FOCUS不能单独确立表达变化因果作用；AlphaFold3、ThermoMPNN和100 ns MD仅提供探索性结构证据，候选基因和错义变异需进一步实验验证。

研究结论（翻译）：本研究使用Genomic SEM构建了整合空腹胰岛素、胰岛素敏感性指数、甘油三酯和腰臀比的胰岛素抵抗相关潜在遗传因子。结果提示，共享IR相关遗传信号与脂肪分布、脂质代谢、脂肪组织功能和调控基因组区域密切相关。精细定位和TWAS/FOCUS分析优先排序出若干参与脂肪细胞生物学、代谢调节和表达控制的候选位点与基因。单独的探索性结构分析提示，所选超罕见错义替换可能改变蛋白稳定性和构象动力学；这些替换是依据预测功能影响和建模适用性而非变异水平关联或精细定位支持选择的。遗传发现提名了有待进一步研究的候选对象，结构结果产生了替换特异性假设，并不能解释或验证常见变异GWAS关联。