This dataset contains the summary statistics from regional association plots from GWAS, and conditional eQTL and pQTL analysis in GTEx data, and functional data as described in Venema et al. 2023 "A cis-regulatory element regulates ERAP2 expression through autoimmune disease risk SNPs see preprint on BioRxiv
This dataverseNL dataset contains additional raw, processed, and metadata (see readme file and reproducible R notebooks (R script and Image) used for the analysis in the manuscript:
R scripts (markdown + R image) with step-by-step analyses
Figure_2.rmd (see "Figure_2.html")
Figure_3.rmd (see "Figure_3.html")
Figure_4.rmd (see "Figure_4.html")
Cell line authentification (including COSMIC reference genotypes) (see folder Appendix in files)
Experimental data
The qPCR (Ct value) data for each experiment are outlined the Supplemental Tables (see here)
The 4C seq data (FASTA and wig files) for the 3 birdshot LCLs can be found in subdirectory of Figure 4
Copy-number SNP-array data for CRISPR-edited cell lines (see folder "SNP_array_Data" and an example of THP-1 can be found here)
Other data used in the manuscript and R scripts
MPRA of ERAP2 eQTLs as reported by Abell et al. Science 2022
H3K27ac Chip combined with HiC in primary immune cell subsets generated by Chandra et al., Nat Genet. 2021
Supplemental data originally generated and provided upon request by Kerner et al., Cell Genomics 2023
Summary statistics of ERAP2 eQTL analysis in whole blood generated by GTEx Consortium
Summary statistics of ERAP2 pQTL analysis in human plasma generated in the INTERVAL study by Sun et al Nature 2018
Summary statistics of 5q15 of GWAS in Birdshot chorioretinopathy (full GWAS summary stats [>1GB] are available in folder "Figure 2/GWAS_summary_stats") generated by Kuiper et al Hum Mol Gen 2014
Summary statistics of GWAS in Crohn's Disease generated by Franke et al Nat Genet 2010
Summary statistics of 5q15 of GWAS in JIA generated by Hinks et al Nat Genet 2013
Single nucleotide polymorphisms (SNP) near the ERAP2 gene are associated with autoimmune conditions such as Crohn’s disease, and birdshot chorioretinopathy, as well as protection against lethal infections, including the Black Death. Due to high linkage disequilibrium (LD), a great number of trait-associated SNPs are correlated with ERAP2 expression, however their functional mechanisms remain unidentified. We used genome editing and functional genomics to identify causal variants that remain obscured by LD. We demonstrate by reciprocal allelic replacement that ERAP2 expression is directly controlled by the genotype of splice region SNP rs2248374. However, we demonstrate that autoimmune disease-risk SNPs located near the downstream LNPEP gene promoter are independently associated with ERAP2 expression. Allele-specific conformation capture assays revealed long-range chromatin contacts between the LNPEP promoter region and the ERAP2 promoter and showed that interactions were stronger in patients carrying the alleles that increase susceptibility to autoimmune diseases. Replacing the disease-associated SNPs in the LNPEP promoter by reference sequences lowered ERAP2 expression. These findings show that clustered GWAS signals associated with diverse autoimmune conditions and lethal infections act in concert to control ERAP2 expression and that disease-associated variants can convert a gene promoter region into a potent enhancer of a distal gene.
