Genome-wide association studies have identified 2p13.1 as a prominent susceptibility locus for systemic lupus erythematosus (SLE)—a complex, multisystem autoimmune disease. However, the identity of underlying causal variant (s) and molecular mechanisms for increasing disease susceptibility are poorly understood. Using meta-analysis (cases = 10,252, controls = 21,604) followed by conditional analysis, bioinformatic annotation, and eQTL and 3D-chromatin interaction analyses, we computationally prioritized potential functional variants and subsequently experimentally validated their effects. Ethnicity-specific meta-analysis revealed striking allele frequency differences between Asian and European ancestries, but with similar odds ratios. We identified 20 genome-wide significant (p < 5 × 10⁻⁸) variants, and conditional analysis pinpointed two potential functional variants, rs6705628 and rs2272165, likely to explain the association. The two SNPs are near DGUOK, mitochondrial deoxyguanosine kinase, and its associated antisense RNA DGUOK-AS1. Using luciferase reporter gene assays, we found significant cell type- and allele-specific promoter activity at rs6705628 and enhancer activity at rs2272165. This is supported by ChIP-qPCR showing allele-specific binding with three histone marks (H3K27ac, H3K4me3, and H3K4me1), RNA polymerase II (Pol II), transcriptional coactivator p300, CCCTC-binding factor (CTCF), and transcription factor ARID3A. Transcriptome data across 28 immune cell types from Asians showed both SNPs are cell-type-specific but only in B-cells. Splicing QTLs showed strong regulation of DGUOK-AS1. Genotype-specific DGOUK protein levels are supported by Western blots. Promoter capture Hi-C data revealed long-range chromatin interactions between rs2272165 and several nearby promoters, including DGUOK. Taken together, we provide mechanistic insights into how two noncoding variants underlie SLE risk at the 2p13.1 locus.

Systemic lupus erythematosus (SLE) is a complex autoimmune disorder characterized by widespread inflammation and autoantibody production. Its development and progression involve genetic, epigenetic, and environmental factors. Although genome-wide association studies (GWAS) have repeatedly identified a susceptibility signal at 16p13, its fine-scale source and its functional and mechanistic role in SLE remain unclear. We used bioinformatics to prioritize likely functional variants and validated the top candidate through various experimental techniques, including clustered regularly interspaced short palindromic repeats (CRISPR)-based genome editing in B cells. To assess the functional impact of the proposed causal variant in C-type lectin domain family 16, member A (CLEC16A), we compared autophagy levels between wild-type (WT) and knock-out (KO) cells. Systematic bioinformatics analysis identified the highly conserved non-coding intronic variant rs17673553, with the risk allele apparently affecting enhancer function and regulating several target genes, including CLEC16A itself. Luciferase reporter assays followed by chromatin immunoprecipitation-quantitative polymerase chain reaction (ChIP-qPCR) validated this enhancer activity, demonstrating that the risk allele increases the binding of enhancer histone marks (H3K27ac and H3K4me1), the CTCF-binding factor, and key immune transcription factors (GATA3 and STAT3). Knock-down of GATA3 and STAT3 via siRNA led to a significant decrease in CLEC16A expression. These regulatory effects on the target gene were further confirmed using CRISPR-based genome editing and CRISPR-dCas9-based epigenetic activation/silencing. Functionally, WT cells exhibited higher levels of starvation-induced autophagy compared to KO cells, highlighting the role of CLEC16A and the rs17673553 locus in autophagy regulation. These findings suggest that the rs17673553 locus—particularly the risk allele—drives significant allele-specific chromatin modifications and binding of multiple transcription factors, thereby mechanistically regulating the expression of target autophagy-associated genes, including CLEC16A itself. This mechanism could potentially explain the association between rs17673553 and SLE, and could underlie the signal at 16p13.

Systemic lupus erythematosus (SLE) is a complex autoimmune disorder characterized by widespread inflammation and autoantibody production. Its development and progression involve genetic, epigenetic, and environmental factors. Although genome-wide association studies (GWAS) have repeatedly identified a susceptibility signal at 16p13, its fine-scale source and its functional and mechanistic role in SLE remain unclear. We used bioinformatics to prioritize likely functional variants and validated the top candidate through various experimental techniques, including CRISPR-based genome editing in B cells. To assess the functional impact of the proposed causal variant in CLEC16A , we compared autophagy levels between wild-type (WT) and knock-out (KO) cells. Systematic bioinformatics analysis identified the highly conserved non-coding intronic variant rs17673553, with the risk allele apparently affecting enhancer function and regulating several target genes, including CLEC16A itself. Luciferase reporter assays followed by ChIP-qPCR validated this enhancer activity, demonstrating that the risk allele increases the binding of enhancer histone marks (H3K27ac and H3K4me1), CTCF-binding factor, and key immune transcription factors (GATA3 and STAT3). Knock-down of GATA3 and STAT3 via siRNA led to a significant decrease in CLEC16A expression. These regulatory effects on the target gene were further confirmed using CRISPR-based genome editing and CRISPR-dCas9-based epigenetic activation/silencing. Functionally, WT cells exhibited higher levels of starvation-induced autophagy compared to KO cells, highlighting the role of CLEC16A and the rs17673553 locus in autophagy regulation. These findings suggest that the rs17673553 locus – particularly the risk allele – drives significant allele-specific chromatin modifications and binding of multiple transcription factors, thereby mechanistically regulating the expression of target autophagy-associated genes, including CLEC16A itself. This mechanism could potentially explain the association between rs17673553 and SLE, and underlie the signal at 16p13.