Massively Parallel Reporter Assay Confirms Regulatory Potential of hQTLs and Reveals Important Variants in Lupus and Other Autoimmune Diseases
Ontology highlight
ABSTRACT:
ORGANISM(S): synthetic construct Homo sapiens
PROVIDER: GSE254502 | GEO | 2024/03/10
REPOSITORIES: GEO
Similar Datasets
Project description:Genome-wide association studies (GWAS) have successfully identified 145 genomic regions that contribute to schizophrenia risk, but linkage disequilibrium (LD) makes it challenging to discern causal variants. Computational finemapping prioritized thousands of credible variants, ~98% of which lie within poorly characterized non-coding regions. To functionally validate their regulatory effects, we performed a massively parallel reporter assay (MPRA) on 5,173 finemapped schizophrenia GWAS variants in primary human neural progenitors (HNPs). We identified 440 variants with allelic regulatory effects (MPRA-positive variants), with 72% of GWAS loci containing at least one MPRA-positive variant. Transcription factor binding had modest predictive power for predicting the allelic activity of MPRA-positive variants, while GWAS association, finemap posterior probability, enhancer overlap, and evolutionary conservation failed to predict MPRA-positive variants. Furthermore, 64% of MPRA-positive variants did not exhibit eQTL signature, suggesting that MPRA could identify yet unexplored variants with regulatory potentials. MPRA-positive variants differed from eQTLs, as they were more frequently located in distal neuronal enhancers. Therefore, we leveraged neuronal 3D chromatin architecture to identify 273 genes that physically interact with MPRA-positive variants. These genes annotated by chromatin interactome displayed higher mutational constraints and regulatory complexity than genes annotated by eQTLs, recapitulating a recent finding that eQTL- and GWAS-detected variants map to genes with different properties. Finally, we propose a model in which allelic activity of multiple variants within a GWAS locus can be aggregated to predict gene expression by taking chromatin contact frequency and accessibility into account. In conclusion, we demonstrate that MPRA can effectively identify functional regulatory variants and delineate previously unknown regulatory principles of schizophrenia.
2023-08-07 | GSE211045 | GEO
Project description:Genome-wide association studies of Systemic Lupus Erythematosus (SLE) nominate 3,073 genetic variants at 91 risk loci. To systematically screen these variants for allelic transcriptional enhancer activity, we constructed a massively parallel reporter assays (MPRA) library comprising 12,396 DNA oligonucleotides containing the genomic context around every allele of each SLE variant. Transfection into Epstein-Barr virus-transformed B cell line GM12878 revealed 482 variants with enhancer activity, with 51 variants showing genotype-dependent (allelic) enhancer activity at 27 risk loci. Combined with allelic enhancer activity analyses in Jurkat cell line, we identified shared and unique allelic transcriptional regulatory mechanisms at SLE risk loci. In-depth analysis of allelic transcription factor (TF) binding at and around 51 allelic variants identified one class of TFs whose DNA-binding motif tends to be directly altered by the risk variant and a second, larger class of TFs that also bind allelically but do not have their motifs directly altered by the variant. Collectively, our approach provides a blueprint for the discovery of allelic gene regulation at risk loci for any disease and offers insight into the transcriptional regulatory mechanisms underlying SLE.
2021-01-31 | GSE143792 | GEO
Project description:The most recent genome-wide association study (GWAS) of cutaneous melanoma identified 54 risk-associated loci, but functional variants and their target genes for most have not been established. Here, we performed massively parallel reporter assays (MPRA) using malignant melanoma and normal melanocyte cells and further integrated multi-layer annotation to systematically prioritize functional variants and susceptibility genes from these GWAS loci. Of 1,992 risk-associated variants tested in MPRA, we identified 285 from 42 loci (78% of the known loci) displaying significant allelic transcriptional activities in either cell type (FDR < 1%). We further characterized MPRA-significant variants by motif prediction, epigenomic annotation, and statistical/functional fine-mapping to create integrative variant scores, which prioritized one to six plausible candidate variants per locus for the 42 loci and nominated a single variant for 43% of these loci. Overlaying the MPRA-significant variants with genome-wide significant expression or methylation quantitative trait loci (QTLs) from melanocytes or melanomas identified candidate susceptibility genes for 60% of variants (172 of 285 variants). CRISPRi of top-scoring variants validated their cis-regulatory effect on the eQTL target genes, MAFF (22q13.1) and GPRC5A (12p13.1). Finally, we identified 36 melanoma-specific and 45 melanocyte-specific MPRA-significant variants, a subset of which are linked to cell-type-specific target genes. Analyses of transcription factor availability in MPRA datasets and variant-transcription factor interaction in eQTL datasets highlighted the roles of transcription factors in cell-type-specific variant functionality. In conclusion, MPRA along with variant scoring effectively prioritized plausible candidates for most melanoma GWAS loci and highlighted cellular contexts where the susceptibility variants are functional.
2022-11-28 | GSE210354 | GEO
Project description:The most recent genome-wide association study (GWAS) of cutaneous melanoma identified 54 risk-associated loci, but functional variants and their target genes for most have not been established. Here, we performed massively parallel reporter assays (MPRA) using malignant melanoma and normal melanocyte cells and further integrated multi-layer annotation to systematically prioritize functional variants and susceptibility genes from these GWAS loci. Of 1,992 risk-associated variants tested in MPRA, we identified 285 from 42 loci (78% of the known loci) displaying significant allelic transcriptional activities in either cell type (FDR < 1%). We further characterized MPRA-significant variants by motif prediction, epigenomic annotation, and statistical/functional fine-mapping to create integrative variant scores, which prioritized one to six plausible candidate variants per locus for the 42 loci and nominated a single variant for 43% of these loci. Overlaying the MPRA-significant variants with genome-wide significant expression or methylation quantitative trait loci (QTLs) from melanocytes or melanomas identified candidate susceptibility genes for 60% of variants (172 of 285 variants). CRISPRi of top-scoring variants validated their cis-regulatory effect on the eQTL target genes, MAFF (22q13.1) and GPRC5A (12p13.1). Finally, we identified 36 melanoma-specific and 45 melanocyte-specific MPRA-significant variants, a subset of which are linked to cell-type-specific target genes. Analyses of transcription factor availability in MPRA datasets and variant-transcription factor interaction in eQTL datasets highlighted the roles of transcription factors in cell-type-specific variant functionality. In conclusion, MPRA along with variant scoring effectively prioritized plausible candidates for most melanoma GWAS loci and highlighted cellular contexts where the susceptibility variants are functional.
2022-11-28 | GSE210355 | GEO
Project description:We applied a massively parallel reporter assay (MPRA) in lymphoblastoid cells to functionally evaluate 49,256 allelic pairs, representing 30,893 genetic variants in high, local linkage disequilibrium for 744 independent cis-expression quantitative trait loci (eQTLs) assessed for colocalization across 114 traits.
2021-05-24 | GSE174534 | GEO
Project description:Genome-wide association studies (GWAS) have identified multiple lung cancer risk loci including those that are distinct in the major histological types. However, most of these loci have not been functionally characterized. Here we employed massively parallel reporter assays (MPRA) to assess allelic transcriptional activity of risk-associated variants en masse. We tested 2,245 variants in 42 loci from 3 recent GWASs of East Asian and European populations in the context of lung adenocarcinoma and squamous cell carcinoma cells while incorporating exposure to a tobacco-smoke carcinogen, benzo[a]pyrene. At FDR < 0.01, we identified 844 MPRA-significant variants with allelic transcriptional activity across 88% of lung cancer loci. Further variant scoring using lung-specific epigenomic annotation demonstrated that 63% of the loci harbored multiple equally functional variants in linkage disequilibrium. Cell-type-specific variants were observed in 72% of the loci, a subset of which aligned with histology-specific association in the GWAS. Distinct subsets of transcription factors were predicted to bind to cell-type-specific variants and those affecting multiple GWAS loci in trans. Linking MPRA-significant variants to target genes based on four different approaches identified candidate susceptibility genes including essential genes for lung cancer cell growth. CRISPR-interference of a high-scoring MPRA-significant variant validated multiple variant-gene connections from different datasets, including RTEL1, SOX18, and ARFRP1. Our data provide a comprehensive catalog of functional characterization of lung cancer GWAS loci and the molecular basis of heterogeneity and polygenicity of lung cancer susceptibility.
2024-05-16 | GSE245699 | GEO
Project description:Genome-wide association studies (GWAS) have boosted our knowledge of genetic risk variants in autoimmune diseases (AIDs). Most of the risk variants are located within or near genes with immunological functions, and the majority is found to be non-coding, pointing towards a regulatory role. We have performed a cis expression quantitative trait locus (eQTL) screen to investigate whether single nucleotide polymorphisms (SNPs) associated with AIDs influence gene expression in thymus. Genotyping was performed using the Immunochip and 353 AID associated SNPs were tested against expression of surrounding genes (+/- 1 Mb) from human thymic tissue (N=42). We identified eight genes where the expression was associated with AID risk SNPs at a study-wide level of significance (P < 2.57x10-5). Five genes (FCRL3, RNASET2, C2orf74, SIRPG and SYS1) displayed cis eQTL signals also in other tissues, while for two loci (NPIPB8 and LOC388814), the eQTL signal appear to be thymus-specific. Since many AID risk variants from GWAS have been subsequently fine-mapped in recent Immunochip projects, we explored the overlap between these novel AID risk variants and the thymic eQTL regions. Moreover, we examined the functional annotation of the seven expression altering SNPs (eSNPs). Our study reveals autoimmune risk variants that act as eQTLs in thymus. We have highlighted functional variants within these genetic regions that potentially can represent causal autoimmune risk variants. Total RNA from 42 human thymic samples were obtained from children undergoing cardiac surgery.
2016-01-21 | E-GEOD-77052 | biostudies-arrayexpress
Project description:Systemic lupus erythematosus (SLE) is a heterogeneous autoimmune disease. Knowledge of circulating immune cell types and cell states associated with SLE remains incomplete. We profiled over 1.2 million PBMCs (162 cases, 99 controls) with multiplexed single-cell RNA-seq (mux-seq). Cases exhibited prominent expression of type-1 interferon-stimulated genes (ISG) in monocytes, reduction of naïve CD4+ T cells that correlated with monocyte ISG expression, and expansion of repertoire-restricted cytotoxic GZMH+ CD8+ T cells. Cell-type-specific expression features accurately predicted case-control status and stratified patients into two molecular subtypes. We integrated dense genotyping data, mapping cell-type-specific cis-eQTLs and linked known and novel SLE-associated variants to cell-type-specific gene expression. These results demonstrate mux-seq as a systematic approach to characterize cellular composition, identify transcriptional signatures, and annotate genetic variants associated with SLE.
2021-12-01 | GSE174188 | GEO
Project description:Although genetic studies have identified many hundreds of loci associated with human traits and diseases, pinpointing the causal alleles remains difficult, particularly for non-coding variants. To address this challenge, we have enhanced the sensitivity and reproducibility of the massively parallel reporter assay (MPRA), adapting it to identify variants that directly modulate gene expression. We then applied it to over 29,000 single nucleotide and insertion/deletion polymorphisms from 3,965 cis-expression quantitative trait loci (eQTL). We demonstrate strong correlation between our MPRA approach and existing measures of regulatory function, and determine an approximate sensitivity of ~20% with a positive predictive value of 60-65% to detect an eQTL causal allele. We identify 842 variants showing differential expression between alleles, including 53 well-annotated variants associated with diseases and traits. Thus, we have created a resource of concrete leads for understanding the genetic basis of specific phenotypes and illustrate the promise of this kind of approach for comprehensively interrogating how non-coding polymorphism shapes human biology. The study consists of two separate MPRA experiments, a 78,958 oligo (79k study) and a 7,500 oligo library (7.5k). For each library we processed independent transfections into two lymphoblastoid cell lines; in total we completed 5 replicates of NA12878 and 3 replicates of NA19239. For the 79k library we also performed 5 replicate transfections into the hepatocarcinoma cell line HepG2. Raw data is provided as Illumina reads of the 20 bp barcode from the RNA extracted 24 hours post transfection as well as from the plasmid library used for transfection. We also provide oligo/barcode combinations from the âgfp vector in the form of a tab delimited file containing the raw sequence reads of the barcode (column 1) and genomic sequence (column 2). Processed count files are unnormalized counts for each oligo acquired by summing all barcode matches together for each replicate.
2016-06-01 | E-GEOD-75661 | biostudies-arrayexpress