《Diabetes, Metabolic Syndrome and Obesity》:Gestational Diabetes Mellitus and Birth Weight: A Bidirectional Mendelian Randomization and Cohort Study in China
编辑推荐:
背景:妊娠期糖尿病(gestational diabetes mellitus, GDM)被视为胎儿巨大儿的驱动因素,然而关于其对出生体重(birth weight, BW)影响的观察性估计易受混杂因素干扰。GDM是否因果性地塑造胎儿生长,以及通过何种共享生物
背景:妊娠期糖尿病(gestational diabetes mellitus, GDM)被视为胎儿巨大儿的驱动因素,然而关于其对出生体重(birth weight, BW)影响的观察性估计易受混杂因素干扰。GDM是否因果性地塑造胎儿生长,以及通过何种共享生物学机制起作用,仍悬而未决。方法:研究人员利用FinnGen R11 GDM数据(16,802例病例)和EGG联盟出生体重全基因组关联研究(genome-wide association study, GWAS)数据(n=286,868)进行双向两样本孟德尔随机化(Mendelian randomization, MR),采用逆方差加权(inverse-variance weighted, IVW,主要方法)、MR-Egger、加权中位数法、加权众数法、RAPS估计量、MR-PRESSO及贝叶斯共定位分析。多组学证据(GSE92772、ST001948、INTERVAL顺式pQTL)作为假说生成证据被整合。在24,448例单胎妊娠队列中,研究人员采用线性和logistic回归、自举中介分析及双重/去偏机器学习(double/debiased machine learning, DML)进行验证。结果:在极端异质性下,GDM至BW的主要IVW估计值为零(β=-0.004, P=0.889),而稳健估计量显示名义上的阳性信号。21个工具变量中有7个被标记为异常值,其中6个与共定位位点重叠。因此,无法排除可靠的因果解释。反向MR(BW至GDM)为探索性分析(5个工具变量,均为异常值)。共定位分析识别出6个位点(CCND2、MTNR1B、GPSM1、ADCY5、CDKAL1、TCF7L2)具有共享因果变异证据(PP.H4≥0.5)。所有已确立的胰岛β细胞/胰岛素分泌基因均提示共享上游生物学机制。在队列中,GDM与巨大儿发生率升高(OR 1.21)及低出生体重(low birth weight, LBW)发生率升高(OR 1.45)相关。调整母体BMI和血脂后,巨大儿关联大幅减弱(OR 1.01, P=0.913)。通过双向中介分析,研究人员识别出5个候选中介因子(BMI、TG、HDL-C、LDL-C和TG/HDL-C比值)在双向中均为显著中介。此外,在DML分析中,高血压疾病与LBW的关联最强(PE β=-0.204;GH β=-0.067)。结论:关于GDM对出生体重具有一致性因果效应的证据仍不确凿。巨大儿的观察性关联主要由母体代谢状态解释。共享的胰岛素分泌生物学更可能反映遗传多效性,而非直接因果通路。高血压疾病与出生体重降低的关联强于GDM,支持对血糖、BMI、血脂和血压进行综合管理。
妊娠期糖尿病与出生体重因果关联的三角验证:一项双向孟德尔随机化、多组学整合与中国大型队列研究
1. 研究背景与问题
妊娠期糖尿病(gestational diabetes mellitus, GDM)是指妊娠期间发生或首次发现的糖耐量异常,全球约15%–17%的妊娠受影响,亚洲人群负担更重。传统临床范式将GDM视为“巨大儿的驱动因素”,多项观察性荟萃分析报告GDM女性发生巨大儿的比值比(odds ratio, OR)为1.5~2.5。然而,GDM对出生体重(birth weight, BW)的影响呈现令人困惑的双向效应:根据Pedersen假说,母体高血糖穿过胎盘,刺激胎儿胰岛β细胞分泌胰岛素,促进胎儿脂肪合成,导致巨大儿;反之,严重GDM可能通过胎盘血管病变和功能障碍引起胎儿生长受限,增加低出生体重(low birth weight, LBW)风险。这两种看似矛盾的结果在GDM人群中并存,其潜在机制尚不明确。观察性研究难以解决因果推断问题,主要受三方面制约:一是严重混杂(母体BMI、年龄、产次、妊娠期增重、家族因素等);二是反向因果(胎儿生长状态可能反向影响母体糖代谢);三是中介效应(GDM可能通过血脂、炎症、胰岛素抵抗等中间路径间接影响BW)。孟德尔随机化(Mendelian randomization, MR)利用遗传变异作为工具变量,可很大程度上摆脱传统混杂和反向因果。尽管已有将两样本MR应用于母体高血糖与新生儿结局的研究,但针对GDM-BW轴的系统性双向两样本MR,并结合共定位、多组学三角验证和大队列验证的研究仍有限。为此,研究人员构建了三层三角验证框架,旨在估计母体GDM倾向对子代出生体重的因果效应,并识别共享生物学机制。
2. 研究设计与方法概述
研究采用“三层三角验证”框架:第一层,利用FinnGen R11 GDM GWAS(16,802例病例)和EGG联盟出生体重GWAS(n=286,868)进行双向两样本MR,运用逆方差加权(IVW-RE,主要方法)、MR-Egger、加权中位数、加权众数、RAPS等五种方法,并通过Cochran's Q、I2、MR-Egger截距检验和MR-PRESSO进行敏感性分析;随后对21个GDM先导信号进行贝叶斯共定位分析。第二层,整合多组学数据:转录组(GSE92772母体外周血RNA-Seq)、代谢组(ST001948孕早期血浆代谢组,20例GDM vs 20例对照)、蛋白质组(INTERVAL顺式pQTL,n=3,301),进行差异表达、WGCNA、代谢通路富集、顺式pQTL重叠和超几何富集分析。第三层,基于2019—2024年常州市妇幼保健院24,448例单胎妊娠队列,采用线性/逻辑回归(M1、M2、M3模型)、双向自举中介分析(2,000次迭代,评估BMI、TG、HDL-C、LDL-C、TC、TG/HDL-C六个候选中介)以及双重/去偏机器学习(double/debiased machine learning, DML,15协变量LassoCV和46协变量ridge规范)进行验证。该研究经常州市妇幼保健院伦理委员会批准(编号2026-K040-01)。
3. 研究结果
双向孟德尔随机化:正向(GDM→BW)经LD clumping后保留21个独立工具变量,主要IVW-RE估计为β=?0.004(P=0.889),无统计学意义;但加权中位数(β=0.051, P=2.005×10?11)、加权众数(β=0.051, P=3.031×10??)和RAPS(β=0.042, P=1.930×10?12)呈名义阳性。Cochran's Q=595(P=3×10?113)提示极端异质性,MR-PRESSO标记7个异常工具变量(其中6个与共定位位点重叠)。由于主要估计量无效、极端异质性及异常值众多,正向MR结果判定为“不确定”,而非“无效”或“阳性”。反向(BW→GDM)仅保留5个工具变量,IVW-RE功效不足(β=?0.511, P=0.872),但加权中位数(β=?2.386)和加权众数(β=?2.383)显示名义保护性关联(OR≈0.09/SD)。留一分析显示剔除MTNR1B先导变异rs10830963后IVW转为显著(β=?2.56, OR=0.078, P=4.4×10?2?),但所有5个工具变量均被MR-PRESSO标记为异常值,且先导变异未通过Steiger方向性检验,故反向结果仅视为探索性。
共定位分析:在21个GDM先导信号中,6个位点与BW GWAS显示共享因果变异证据(PP.H4≥0.5):CCND2(PP.H4=1.00)、MTNR1B(PP.H4=1.00)、GPSM1(PP.H4=0.99)、ADCY5(PP.H4=0.98)、CDKAL1(PP.H4=0.98)、TCF7L2(PP.H4=0.76)。这6个基因均为已确立的胰岛β细胞/胰岛素分泌基因。敏感性分析显示其中5个位点对先验假设稳健,TCF7L2对p12敏感,需谨慎解读。
多组学整合:转录组层,GSE92772中1,864个基因达名义显著(P<0.05),755个基因达联合显著(P<0.05, |log2FC|>0.5),WGCNA识别出hub基因PLA2G7、SGPP1、GNPDA2、KLF9,涉及脂质代谢、β细胞功能和炎症。代谢组层,GDM组中鞘磷脂类下调,黄嘌呤、次黄嘌呤、2-羟基丁酸和CAR 16:0上调;MSEA富集到6条通路(FDR<0.05);自举中介显示4种鞘磷脂为GDM→BW负向中介,4种嘌呤/酮酸代谢物为正向中介。蛋白质组层,185个GDM先导SNP中有36个与cis-pQTL重叠,形成5个位点簇,其中GCKR为最大簇(15/36,占41.7%,β=+0.697, P=6.46×10??)。跨层富集分析显示GO_INSULIN_SECRETION、GO_BETA_CELL和GO_KEGG_T2D三条通路均标记TCF7L2和SLC30A8。综上,多组学证据指向胰岛素分泌能力变异作为GDM易感性与胎儿生长的共享上游驱动因素。
队列验证:24,448例队列中,GDM组与非GDM组新生儿出生体重无显著差异(3.34±0.49 vs 3.34±0.45 kg, P=0.7249),但GDM组巨大儿和LBW发生率均显著升高。线性回归中,GDM系数从M1(β=?0.0028, P=0.725)到M2(β=?0.0041, P=0.609)再到M3(调整BMI和血脂后β=?0.0309, P=5.0×10??)转为显著,符合统计抑制效应。巨大儿OR从M1的1.21(P=0.0038)降至M3的1.01(P=0.9127),表明巨大儿关联主要由母体BMI和血脂解释;LBW的OR从1.45增至1.55(P<0.001)。双向中介分析显示12条路径中10条(83.3%)存在显著双向中介效应,主要中介为TG/HDL-C、TG、BMI、HDL-C和LDL-C。DML分析显示,高血压疾病对出生体重的负向效应显著大于GDM:子痫前期(preeclampsia, PE)β=?0.204(P=3.3×10?3?,为GDM的7.8倍),妊娠期高血压(gestational hypertension, GH)β=?0.067(P=8.5×10??,2.6倍),妊娠期肝内胆汁淤积症(intrahepatic cholestasis of pregnancy, ICP)β=?0.059(P=0.001,2.3倍),而GDM本身效应接近零。敏感性分析中,在可获得空腹生化数据的14,967例女性中,每SD升高的甘油三酯-葡萄糖-体质指数(triglyceride-glucose body mass index, TyG-BMI)与出生体重+0.108 kg相关(95% CI 0.100~0.115),大于单独BMI的效应(+0.096),且与巨大儿的关联更强(OR 1.93 vs 1.78)。
4. 讨论与结论总结
讨论部分指出,传统将非显著IVW估计解读为“无因果效应”并不恰当,因为极端异质性和异常工具变量使点估计不可靠,正向GDM→BW的因果效应属于“不确定”。这并不与随机对照试验中治疗GDM可减少巨大儿的证据矛盾,因为试验干预的是整个因果包(血糖、BMI、血脂、护理强度等),而非孤立的遗传效应。反向MR虽名义支持“母胎代谢共生”假说,但因Steiger检验失败及样本重叠风险,仅具探索价值。MR-PRESSO标记的正向异常值位点与共定位位点高度重叠(6/7),提示GDM与BW共享上游胰岛素分泌生物学,而非直接因果通路。队列中巨大儿关联在调整母体代谢状态后大幅衰减,而LBW关联增强,调和了GDM中巨大儿与LBW共存的悖论。高血压疾病对出生体重降低的效应远大于GDM,提示临床管理中血压监测应与血糖控制同等重要。
研究结论翻译:在这项三角验证研究中,研究人员未发现GDM对出生体重具有一致性因果效应。巨大儿的观察性关联在调整母体BMI和血脂后大幅减弱。反向MR为探索性分析,不支持遗传预测的出生体重对母体GDM风险具有保护效应。多组学证据与共定位分析共同指向共享的胰岛素分泌生物学,作为GDM-BW关联的假说生成解释,有待独立数据集验证。高血压疾病与出生体重降低的关联强于GDM。这些发现支持涵盖BMI、血脂、血压和血糖的综合围产期管理。