A functional variant rs912304 for late-onset T1D risk contributes to islet dysfunction by regulating proinflammatory cytokine-responsive gene STXBP6 expression

Abstract

BackgroundOur previous genome‑wide association studies (GWAS) have suggested rs912304 in 14q12 as a suggestive risk variant for type 1 diabetes (T1D). However, the association between this risk region and T1D subgroups and related clinical risk features, the underlying causal functional variant(s), putative candidate gene(s), and related mechanisms are yet to be elucidated.MethodsWe assessed the association between variant rs912304 and T1D, as well as islet autoimmunity and islet function, stratified by the diagnosed age of 12. We used epigenome bioinformatics analyses, dual luciferase reporter assays, and expression quantitative trait loci (eQTL) analyses to prioritize the most likely functional variant and potential causal gene. We also performed functional experiments to evaluate the role of the causal gene on islet function and its related mechanisms.ResultsWe identified rs912304 as a risk variant for T1D subgroups with diagnosed age ≥ 12 but not < 12. This variant is associated with residual islet function but not islet-specific autoantibody positivity in T1D individuals. Bioinformatics analysis indicated that rs912304 is a functional variant exhibiting spatial overlaps with enhancer active histone marks (H3K27ac and H3K4me1) and open chromatin status (ATAC-seq) in the human pancreas and islet tissues. Luciferase reporter gene assays and eQTL analyses demonstrated that the biallelic sites of rs912304 had differential allele-specific enhancer activity in beta cell lines and regulated STXBP6 expression, which was defined as the most putative causal gene based on Open Targets Genetics, GTEx v8 and Tiger database. Moreover, Stxbp6 was upregulated by T1D-related proinflammatory cytokines but not high glucose/fat. Notably, Stxbp6 over-expressed INS-1E cells exhibited decreasing insulin secretion and increasing cell apoptosis through Glut1 and Gadd45β, respectively.ConclusionsThis study expanded the genomic landscape regarding late-onset T1D risk and supported islet function mechanistically connected to T1D pathogenesis.