Direct identification of hundreds of expression-modulating variants using a multiplexed reporter assay
Ontology highlight
ABSTRACT:
ORGANISM(S): synthetic construct Homo sapiens
PROVIDER: GSE75661 | GEO | 2016/06/01
SECONDARY ACCESSION(S): PRJNA305000
REPOSITORIES: GEO
Similar Datasets
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
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 (GWAS-s) have linked thousands of genetic variants with complex diseases and traits. However, most genetic variants identified reside in non-coding regions of the genome making it hard to predict and understand the molecular mechanisms underlying causal alleles. We have previously performed molecular Quantitative Trait Locus (molQTL) analysis for transcription factor binding, chromatin accessibility, and H3K27 acetylation across primary Human Aortic Endothelial Cell (HAEC) samples. Here we expand this analysis to study more than 30,000 genetic variants using the massively parallel reporter assay (MPRA) STARR-Seq in immortalized HAEC cells (teloHAEC). We demonstrate that more than 5,000 variants exhibit differential expression between alleles and identify ETS and AP-1 motifs as the strongest sequence attributes—and cell type-specific chromatin accessibility as the strongest epigenetic attribute—associated with allele-specific regulatory activity. Using interleukin 1-beta (IL1b) as a modulator of cell state and environment, we observe robust evidence for context-specific SNP effects, thereby underscoring the prevalence of GxE effects on noncoding variant function. We integrate results from molQTL mapping and STARR-seq with eQTL data from HAECs and GTEx tissues to fine-map functional non-coding SNPs at GWAS loci for vascular diseases.
2021-12-21 | GSE180846 | 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: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 observation that variants regulating gene expression (expression quantitative trait loci, eQTL) are at a high frequency among SNPs associated with complex traits has made the genome-wide characterization of gene expression an important tool in genetic mapping studies of such traits. As part of a study to identify genetic loci contributing to bipolar disorder and a wide range of BP-related quantitative traits in members of 26 pedigrees from Costa Rica and Colombia, we measured gene expression in lymphoblastoid cell lines derived from 786 pedigree members. The study design enabled us to comprehensively reconstruct the genetic regulatory network in these families, provide estimates of heritability, identify eQTL, evaluate missing heritability for the eQTL, and quantify the number of different alleles contributing to any given locus.
2016-07-01 | GSE82042 | GEO
Project description:Genome-wide association studies (GWAS) have successfully identified thousands of associations between common genetic variants and human disease phenotypes, but the majority of these variants are non-coding, often requiring genetic fine-mapping, epigenomic profiling, and individual reporter assays to delineate potential causal variants. We employ a massively parallel reporter assay (MPRA) to simultaneously screen 2,756 variants in strong linkage disequilibrium with 75 sentinel variants associated with red blood cell traits. We show that this assay identifies elements with endogenous erythroid regulatory activity. Across 23 sentinel variants, we conservatively identified 32 MPRA functional variants (MFVs). We used targeted genome editing to demonstrate endogenous enhancer activity across 3 MFVs that predominantly affect the transcription of SMIM1, RBM38, and CD164. Functional follow-up of RBM38 delineates a key role for this gene in the alternative splicing program occurring during terminal erythropoiesis. Finally, we provide evidence for how common GWAS-nominated variants can disrupt cell-type-specific transcriptional regulatory pathways.
2016-10-07 | GSE87711 | GEO
Project description:Genome-wide association studies (GWAS) have successfully identified thousands of associations between common genetic variants and human disease phenotypes, but the majority of these variants are non-coding, often requiring genetic fine-mapping, epigenomic profiling, and individual reporter assays to delineate potential causal variants. We employ a massively parallel reporter assay (MPRA) to simultaneously screen 2,756 variants in strong linkage disequilibrium with 75 sentinel variants associated with red blood cell traits. We show that this assay identifies elements with endogenous erythroid regulatory activity. Across 23 sentinel variants, we conservatively identified 32 MPRA functional variants (MFVs). We used targeted genome editing to demonstrate endogenous enhancer activity across 3 MFVs that predominantly affect the transcription of SMIM1, RBM38, and CD164. Functional follow-up of RBM38 delineates a key role for this gene in the alternative splicing program occurring during terminal erythropoiesis. Finally, we provide evidence for how common GWAS-nominated variants can disrupt cell-type-specific transcriptional regulatory pathways.
2016-10-07 | GSE87710 | GEO
Project description:DNA variants that alter gene expression contribute to variation in many phenotypic traits. In particular, trans-acting variants, which are often located on different chromosomes from the genes they affect, are an important source of heritable gene expression variation. However, our knowledge about the identity and mechanism of causal trans-acting variants remains limited. Here, we developed a fine-mapping strategy called CRISPR-Swap and dissected three expression quantitative trait locus (eQTL) hotspots known to alter the expression of multiple genes in trans in the yeast Saccharomyces cerevisiae. Causal variants were identified by engineering recombinant alleles and quantifying the effects of these alleles on the expression of a green fluorescent protein-tagged gene affected by the given locus in trans. We validated the effect of each variant on the expression of multiple genes by RNA-sequencing. The three variants were strikingly different in their molecular mechanism, the type of genes they reside in, and their distribution in natural populations. While a missense leucine-to-serine variant at position 63 in the transcription factor Oaf1 (L63S) was almost exclusively present in the reference laboratory strain, the two other variants were frequent among S. cerevisiae isolates. A causal missense variant in the glucose receptor Rgt2 (V539I) occurred at a poorly conserved amino acid residue and its effect was strongly dependent on the concentration of glucose in the culture medium. A noncoding variant in the conserved fatty acid regulated (FAR) element of the OLE1 promoter influenced the expression of the fatty acid desaturase Ole1 in cis and, by modulating the level of this essential enzyme, other genes in trans. The OAF1 and OLE1 variants showed a non-additive genetic interaction, and affected cellular lipid metabolism. These results revealed remarkable diversity in the molecular basis of trans-regulatory variation, highlighting the challenges in predicting which natural genetic variants affect gene expression.
2019-08-12 | GSE134169 | GEO
Project description:Eosinophilic esophagitis (EoE) is a rare atopic disorder associated with esophageal dysfunction, including difficulty swallowing, food impaction, and inflammation. EoE develops in a small subset of people with food allergies under the influence of environmental and genetic risk factors. Genome wide association studies (GWAS) have identified 31 independent risk loci for the disease, and linkage disequilibrium (LD) expansion of these loci nominates a set of 531 variants that are potentially causal. These risk variants are non-coding, suggesting a likely role in altered gene regulatory mechanisms. To systematically interrogate the gene regulatory activity of these variants, we designed a massively parallel reporter assay (MPRA) containing the alleles of each variant within their 170 bp genomic sequence context cloned into a GFP reporter library. Transfection of the MPRA library into TE-7 esophageal epithelial cells, HaCaT skin keratinocytes, and Jurkat T cells revealed cell-type-specific gene regulation. We identify 32 allelic enhancer variants (allelic enVars) that regulate reporter gene expression in a genotype-dependent manner in at least one cellular context. By annotating these variants with expression quantitative trait loci (eQTL) and chromatin looping data in related tissues and cell types, we identify putative target genes affected by genetic variation in EoE patients, including TSLP and multiple genes at the HLA locus. Transcription factor enrichment analyses reveal possible roles for cell-type specific regulators, including GATA-3, a key regulator of type 2 inflammation. Collectively, our approach reduces the large set of EoE-associated variants to a set of 32 with allelic regulatory activity, providing new functional insights into the effects of genetic variation in this disease.
2023-11-26 | GSE232337 | GEO