按HLA-DR3与HLA-DR4状态分层的遗传关联分析揭示1型糖尿病进展通路的异质性

《Diabetologia》：Genetic association stratified by HLA-DR3 and HLA-DR4 status reveals heterogeneity in pathways of progression to type 1 diabetes

　　

目的/假设：1型糖尿病（T1D）的临床表型因首现自身抗体不同存在显著差异，但机制不明。主要组织相容性复合体（MHC）上的DR3-DQ2（DR3）与DR4-DQ8（DR4）单倍型分别关联谷氨酸脱羧酶抗体（GADA）与胰岛素自身抗体（IAA）作为首现自身抗体，因此可作为大型队列中首现自身抗体的替代指标。方法：研究人员利用来自多队列的9091例T1D病例与14157例对照，开展首项按DR3和DR4状态分层的T1D全基因组关联研究（GWAS）。估计DR3-T1D与DR4-T1D及与其他免疫和血糖表型的遗传力与遗传相关，评估已知T1D位点在两组间的效应异质性，判定遗传力在细胞类型特异性顺式调控元件（cRE）与生物学通路中的富集，并注释cRE中风险变异。结果：DR4-T1D与DR3-T1D仅呈中度遗传相关（rg=0.6），低于按发病年龄和性别分层的相关，且与其他自身免疫病的相关模式各异。已知位点中，IL2位点在DR4中对T1D效应显著更大，TAGAP、KLRG1呈名义异质性。DR4-T1D变异更强富集于T细胞cRE与T细胞相关通路；DR3-T1D变异特异性富集于肥大细胞cRE与分泌相关通路。研究人员进一步在肥大细胞中优先定位对DR3个体T1D效应更强的位点。结论：首项按DR3/DR4分层的T1D GWAS揭示依赖高危HLA背景的遗传风险与生物学机制异质性。

研究背景与立项依据：1型糖尿病（T1D）以自身免疫性β细胞破坏为特征，进展时间从幼年起病至成人起病差异巨大。按首现自身抗体可划分内型（endotype）：谷氨酸脱羧酶抗体（GADA）首发者起病较晚，胰岛素自身抗体（IAA）首发者起病较早。主要组织相容性复合体（MHC）区域的高危单倍型HLA-DR3-DQ2（DR3）与HLA-DR4-DQ8（DR4）分别偏好关联GADA首发与IAA首发，且DR4携带者对靶向T细胞的干预（如teplizumab）反应更佳，提示不同HLA背景下进展机制不同。但由于队列中首现自身抗体数据稀疏且具瞬时性，直接按自身抗体分层不可行，因此研究人员以DR3/DR4状态作为首现自身抗体的遗传代理变量，开展首项按此分层的T1D全基因组关联研究（GWAS），以揭示异质性通路。该研究发表于《Diabetologia》。

主要关键技术方法：研究人员整合dbGaP、EGA及英国Exeter等11个欧洲祖源队列共17413例T1D病例与37570例对照，基于TOPMed与Michigan HLA参比面板进行基因型填补，按四位HLA等位基因定义DR3（DRB1 * 03:01–DQB1 * 02:01）与DR4（DRB1 * 04:01/02/04/05/08–DQB1 * 03:02/04）分组并排除DR3/DR4杂合子。病例-对照关联在DR3与DR4组内分别用EPACTS（Firth校正逻辑回归）与REGENIE分析，经固定效应逆方差荟萃合并。遗传力与遗传相关用LD评分回归（LDSC）估计；位点效应异质性用PLINK Firth交互检验与GWAMA随机效应荟萃；细胞类型特异性顺式调控元件（cRE）富集用LDSC分区遗传力（CATlas注释）；通路富集用MAGMA（GO/KEGG/REACTOME）；肥大细胞cRE内可信集变异效应比较基于已发表136个T1D信号99%可信集。

研究结果：

基因组-wide关联研究按DR3/DR4分层：在DR3-T1D（4057例/8043对照）与DR4-T1D（5034例/6114对照）中，分别有11个独立位点达全基因组显著（p<5×10^?8），均为既往T1D位点。其中PTPN22、IL2RA、INS、IKZF4、SH2B3、PTPN2与MHC两组均显著；CTLA4与HORMAD2仅DR3组显著，CLEC16A与IL2仅DR4组显著。

遗传力与遗传相关估计：DR3-T1D总遗传力h²=0.45，DR4-T1D h²=0.40（含HLA区），去除HLA区后分别为0.17与0.14，组间无统计差异。两组遗传相关r_g=0.68（去HLA区），显著低于按发病年龄（r_g=0.96）与性别（r_g=0.88）分层，水平近似克罗恩病-溃疡性结肠炎或精神分裂症-双相障碍间相关。对外周自身免疫病相关显示：两者均与类风湿关节炎、自身免疫性甲状腺病正相关；系统性红斑狼疮相关在DR4更强；原发性胆汁性胆管炎仅与DR3-T1D相关，乳糜泻与白癜风仅与DR4-T1D相关。

DR3与DR4个体中变异对T1D效应的异质性：在10个非MHC显著位点交互检验中，IL2位点在DR4效应更强（β_DR4=0.230 vs β_DR3=0.0878，FDR<0.10显著）；PTPN22在DR3略强但未过校正。扩至87个已知T1D位点后，TAGAP（DR4更强）与KLRG1（DR3更强）呈名义异质性（p<0.05），NOTCH2等呈Q检验名义异质，但均未过多重校正。

DR3与DR4中T1D关联的細胞类型特异性富集：LDSC分区显示两组均在CD4⁺与CD8⁺ T细胞cRE富集；DR3-T1D额外在肥大细胞cRE显著富集且强于DR4（p=0.043）。MAGMA通路分析示DR4-T1D富集T细胞活化、凋亡、免疫应答通路；DR3-T1D富集分泌囊泡、应激颗粒等分泌相关通路。在肥大细胞cRE重叠的7个T1D可信集位点中5个在DR3效应更强，例如FLI1（组胺合成调控）与10q24的SUFU（GLI转录因子调控）在DR3组有亚基因组显著而在DR4无关联。

讨论总结：研究人员指出DR3-T1D与DR4-T1D的遗传相关仅中度，低于传统年龄/性别分层，支持T1D在高危HLA背景下属遗传异质子型。IL2位点效应偏向DR4印证T细胞中枢作用，与DR4者对teplizumab反应佳一致；DR3-T1D特异富集肥大细胞cRE与分泌通路为既往GWAS未发现，因DR4信号掩盖。作者承认DR3/DR4仅为首现自身抗体代理、AAB瞬时性与队列量限制需前瞻队列（TEDDY、TrialNet）校核；DR3/DR4杂合子风险最高但对照不足未分入；非DR3/DR4 T1D约占10%且机制不同。结论部分翻译如下：研究人员开展了首项按DR3和DR4状态分层的T1D GWAS，揭示了依赖高危HLA背景的遗传风险与生物学机制异质性。