A Multilayered Post-Genome-Wide Association Study Analysis Pipeline Defines Functional Variants and Target Genes for Systemic Lupus Erythematosus.

OBJECTIVE: Systemic lupus erythematosus (SLE), an autoimmune disease with incompletely understood etiology, has a strong genetic component. Although genome-wide association studies (GWASs) have revealed multiple SLE susceptibility loci and associated single-nucleotide polymorphisms (SNPs), the precise causal variants, target genes, cell types, tissues, and mechanisms of action remain largely unknown. METHODS: Here, we report a comprehensive post-GWAS analysis using extensive bioinformatics, molecular modeling, and integrative functional genomic and epigenomic analyses to optimize fine-mapping. We compile and cross-reference immune cell-specific expression quantitative trait loci (cis- and trans-expression quantitative trait loci) with promoter capture high-throughput capture chromatin conformation (PCHi-C), allele-specific chromatin accessibility, and massively parallel reporter assay data to define predisposing variants and target genes. We experimentally validate a predicted locus using CRISPR/Cas9 genome editing, quantitative polymerase chain reaction, and Western blot. RESULTS: Anchoring on 452 index SNPs, we selected 9,931 high linkage disequilibrium (r2 > 0.8) SNPs and defined 182 independent non-human leukocyte antigen (HLA) SLE loci. The 3,746 SNPs from 143 loci were identified as regulating 564 unique genes. Target genes are enriched in lupus-related tissues and associated with other autoimmune diseases. Of these, 329 SNPs (106 loci) showed significant allele-specific chromatin accessibility and/or enhancer activity, indicating regulatory potential. Using CRISPR/Cas9, we validated reference SNP identifier 57668933 (rs57668933) as a functional variant regulating multiple targets, including SLE-risk gene ELF1 in B cells. CONCLUSION: We demonstrate and validate post-GWAS strategies for using multidimensional data to prioritize likely causal variants with cognate gene targets underlying SLE pathogenesis. Our results provide a catalog of significantly SLE-associated SNPs and loci, target genes, and likely biochemical mechanisms to guide experimental characterization.

A Non-Coding Variant in SLC15A4 Modulates Enhancer Activity and Lysosomal Deacidification Linked to Lupus Susceptibility.

Systemic lupus erythematosus (SLE) is a complex autoimmune disease with a strong genetic basis. Despite the identification of several single nucleotide polymorphisms (SNPs) near the SLC15A4 gene that are significantly associated with SLE across multiple populations, specific causal SNP(s) and molecular mechanisms responsible for disease susceptibility are unknown. To address this gap, we employed bioinformatics, expression quantitative trait loci (eQTLs), and 3D chromatin interaction analysis to nominate a likely functional variant, rs35907548, in an active intronic enhancer of SLC15A4 . Through luciferase reporter assays followed by chromatin immunoprecipitation (ChIP)-qPCR, we observed significant allele-specific enhancer effects of rs35907548 in diverse cell lines. The rs35907548 risk allele T is associated with increased regulatory activity and target gene expression, as shown by eQTLs and chromosome conformation capture (3C)-qPCR. The latter revealed long-range chromatin interactions between the rs35907548 enhancer and the promoters of SLC15A4, GLTLD1 , and an uncharacterized lncRNA. The enhancer-promoter interactions and expression effects were validated by CRISPR/Cas9 knock-out (KO) of the locus in HL60 promyeloblast cells. KO cells also displayed dramatically dysregulated endolysosomal pH regulation. Together, our data show that the rs35907548 risk allele affects multiple aspects of cellular physiology and may directly contribute to SLE.

A multilayered post-GWAS analysis pipeline defines functional variants and target genes for systemic lupus erythematosus (SLE).

Objectives: Systemic lupus erythematosus (SLE), an autoimmune disease with incompletely understood etiology, has a strong genetic component. Although genome-wide association studies (GWAS) have revealed multiple SLE susceptibility loci and associated single nucleotide polymorphisms (SNPs), the precise causal variants, target genes, cell types, tissues, and mechanisms of action remain largely unknown. Methods: Here, we report a comprehensive post-GWAS analysis using extensive bioinformatics, molecular modeling, and integrative functional genomic and epigenomic analyses to optimize fine-mapping. We compile and cross-reference immune cell-specific expression quantitative trait loci ( cis - and trans -eQTLs) with promoter-capture Hi-C, allele-specific chromatin accessibility, and massively parallel reporter assay data to define predisposing variants and target genes. We experimentally validate a predicted locus using CRISPR/Cas9 genome editing, qPCR, and Western blot. Results: Anchoring on 452 index SNPs, we selected 9,931 high-linkage disequilibrium (r 2 >0.8) SNPs and defined 182 independent non-HLA SLE loci. 3,746 SNPs from 143 loci were identified as regulating 564 unique genes. Target genes are enriched in lupus-related tissues and associated with other autoimmune diseases. Of these, 329 SNPs (106 loci) showed significant allele-specific chromatin accessibility and/or enhancer activity, indicating regulatory potential. Using CRISPR/Cas9, we validated rs57668933 as a functional variant regulating multiple targets, including SLE risk gene ELF1 , in B-cells. Conclusion: We demonstrate and validate post-GWAS strategies for utilizing multi-dimensional data to prioritize likely causal variants with cognate gene targets underlying SLE pathogenesis. Our results provide a catalog of significantly SLE-associated SNPs and loci, target genes, and likely biochemical mechanisms, to guide experimental characterization.

A Non-Coding Variant in SLC15A4 Modulates Enhancer Activity and Lysosomal Deacidification Linked to Lupus Susceptibility.

Systemic lupus erythematosus (SLE) is a complex autoimmune disease with a strong genetic basis. Despite the identification of several single nucleotide polymorphisms (SNPs) near the SLC15A4 gene that are significantly associated with SLE across multiple populations, specific causal SNP(s) and molecular mechanisms responsible for disease susceptibility are unknown. To address this gap, we employed bioinformatics, expression quantitative trait loci (eQTLs), and 3D chromatin interaction analysis to nominate a likely functional variant, rs35907548, in an active intronic enhancer of SLC15A4. Through luciferase reporter assays followed by chromatin immunoprecipitation (ChIP)-qPCR, we observed significant allele-specific enhancer effects of rs35907548 in diverse cell lines. The rs35907548 risk allele T is associated with increased regulatory activity and target gene expression, as shown by eQTLs and chromosome conformation capture (3C)-qPCR. The latter revealed long-range chromatin interactions between the rs35907548 enhancer and the promoters of SLC15A4, GLTLD1, and an uncharacterized lncRNA. The enhancer-promoter interactions and expression effects were validated by CRISPR/Cas9 knock-out (KO) of the locus in HL60 promyeloblast cells. KO cells also displayed dramatically dysregulated endolysosomal pH regulation. Together, our data show that the rs35907548 risk allele affects multiple aspects of cellular physiology and may directly contribute to SLE.

Connecting omics signatures and revealing biological mechanisms with iLINCS.

There are only a few platforms that integrate multiple omics data types, bioinformatics tools, and interfaces for integrative analyses and visualization that do not require programming skills. Here we present iLINCS ( http://ilincs.org ), an integrative web-based platform for analysis of omics data and signatures of cellular perturbations. The platform facilitates mining and re-analysis of the large collection of omics datasets (>34,000), pre-computed signatures (>200,000), and their connections, as well as the analysis of user-submitted omics signatures of diseases and cellular perturbations. iLINCS analysis workflows integrate vast omics data resources and a range of analytics and interactive visualization tools into a comprehensive platform for analysis of omics signatures. iLINCS user-friendly interfaces enable execution of sophisticated analyses of omics signatures, mechanism of action analysis, and signature-driven drug repositioning. We illustrate the utility of iLINCS with three use cases involving analysis of cancer proteogenomic signatures, COVID 19 transcriptomic signatures and mTOR signaling.

Discovery and Functional Characterization of Two Regulatory Variants Underlying Lupus Susceptibility at 2p13.1.

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−8) 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.

Lupus Susceptibility Region Containing CDKN1B rs34330 Mechanistically Influences Expression and Function of Multiple Target Genes, Also Linked to Proliferation and Apoptosis.

OBJECTIVE: In a recent genome-wide association study, a significant genetic association between rs34330 of CDKN1B and risk of systemic lupus erythematosus (SLE) in Han Chinese was identified. This study was undertaken to validate the reported association and elucidate the biochemical mechanisms underlying the effect of the variant. METHODS: We performed an allelic association analysis in patients with SLE, followed by a meta-analysis assessing genome-wide association data across 11 independent cohorts (n = 28,872). In silico bioinformatics analysis and experimental validation in SLE-relevant cell lines were applied to determine the functional consequences of rs34330. RESULTS: We replicated a genetic association between SLE and rs34330 (meta-analysis P = 5.29 × 10-22 , odds ratio 0.84 [95% confidence interval 0.81-0.87]). Follow-up bioinformatics and expression quantitative trait locus analysis suggested that rs34330 is located in active chromatin and potentially regulates several target genes. Using luciferase and chromatin immunoprecipitation-real-time quantitative polymerase chain reaction, we demonstrated substantial allele-specific promoter and enhancer activity, and allele-specific binding of 3 histone marks (H3K27ac, H3K4me3, and H3K4me1), RNA polymerase II (Pol II), CCCTC-binding factor, and a critical immune transcription factor (interferon regulatory factor 1 [IRF-1]). Chromosome conformation capture revealed long-range chromatin interactions between rs34330 and the promoters of neighboring genes APOLD1 and DDX47, and effects on CDKN1B and the other target genes were directly validated by clustered regularly interspaced short palindromic repeat (CRISPR)-based genome editing. Finally, CRISPR/dead CRISPR-associated protein 9-based epigenetic activation/silencing confirmed these results. Gene-edited cell lines also showed higher levels of proliferation and apoptosis. CONCLUSION: Collectively, these findings suggest a mechanism whereby the rs34330 risk allele (C) influences the presence of histone marks, RNA Pol II, and IRF-1 transcription factor to regulate expression of several target genes linked to proliferation and apoptosis. This process could potentially underlie the association of rs34330 with SLE.

LINCS Data Portal 2.0: next generation access point for perturbation-response signatures.

The Library of Integrated Network-Based Cellular Signatures (LINCS) is an NIH Common Fund program with the goal of generating a large-scale and comprehensive catalogue of perturbation-response signatures by utilizing a diverse collection of perturbations across many model systems and assay types. The LINCS Data Portal (LDP) has been the primary access point for the compendium of LINCS data and has been widely utilized. Here, we report the first major update of LDP (http://lincsportal.ccs.miami.edu/signatures) with substantial changes in the data architecture and APIs, a completely redesigned user interface, and enhanced curated metadata annotations to support more advanced, intuitive and deeper querying, exploration and analysis capabilities. The cornerstone of this update has been the decision to reprocess all high-level LINCS datasets and make them accessible at the data point level enabling users to directly access and download any subset of signatures across the entire library independent from the originating source, project or assay. Access to the individual signatures also enables the newly implemented signature search functionality, which utilizes the iLINCS platform to identify conditions that mimic or reverse gene set queries. A newly designed query interface enables global metadata search with autosuggest across all annotations associated with perturbations, model systems, and signatures.