Massively parallel reporter assays identify functional enhancer variants at QT interval GWAS loci
Ontology highlight
ABSTRACT:
ORGANISM(S): Mus musculus Escherichia coli
PROVIDER: GSE286052 | GEO | 2026/07/07
REPOSITORIES: GEO
Dataset's files
| Action | DRS | |||
|---|---|---|---|---|
| Other | ||||
| GSE286052_RAW.tar | Raw | |||
| GSE286052_mpra_qtigc_pgl4v1_pool1.txt.gz | Txt | |||
| GSE286052_mpra_qtigc_pgl4v1_pool2.txt.gz | Txt | |||
| GSE286052_mpra_qtigc_pgl4v2_pool1.txt.gz | Txt |
Items per page: 5 1 - 5 of 7 |
Similar Datasets
Project description:Genome-wide association studies (GWAS) have suggested that sequence variation in cis-regulatory elements (CREs) modulate common disease risk and trait variation. However, identification of the majority of these causal variants and their mechanism of action have remained significant challenges. We used reporter assays for all common variants at the QT interval associated SCN5A GWAS locus with the goal of identifying causal variants for the association. A target region of ~500kb containing SCN5A was defined based on recombination hotspots (rate>10cM/Mb; estimated from HapMap) flanking the 5 independent QT interval GWAS hits. Within this region, all common variants (minor allele frequency >5%), from the 1000 Genomes European ancestry populations, in moderate linkage disequilibrium (LD; r2>0.3) with any of the 5 sentinel hits, were selected for analyses. Of a total 121 variants selected, 112 variants in 104 amplicons passed primer design (amplicon size 256-617bp; median 397bp), with both alleles of 106 variants successfully cloned along with flanking sequences into 98 amplicons upstream of a minimal promoter-driven firefly luciferase gene in pGL4.23. Cardiomyocyte cells, AC16 (human) and HL1 (mouse), were transfected with test constructs and Renilla luciferase vector (for transfection normalization) in triplicate; luciferase assays were performed 24h later. All cloning and reporter assays were performed in 96- and 24-well plates. In reporter assays across the two cell lines, compared to empty vector, 37 amplicons showed enhancer activity (z-score>99 %tile), with a concordance rate of ~70%. Of these 37 enhancer CREs, 9 overlapped open chromatin regions (DNase-seq peaks) observed in adult human heart tissue, largest among all human tissues evaluated by the RoadMap Epigenomics project. Of these 37 enhancer CREs, 12 showed allelic difference in reporter activity (P<0.05), thus identifying at least one enhancer CRE variant in high-to-moderate LD with each of the 5 sentinel hits. Using GTEx heart left ventricle (n=190) gene expression data, we showed correlation between SCN5A expression and the number of QT interval prolonging alleles across the 5 index SNPs. We are currently assessing the binding of cardiac nuclear factors at the 12 enhancer CRE variants by gel-shift assays and the effect of CRISPR-Cas9 mediated CRE deletion on SCN5A expression in AC16 and HL1 cells, an analyses helped by evaluation of AC16 and HL1 cells by RNA-seq, ATAC-seq and karyotyping. Thus, independent of the publicly available epigenomic data, which are of limited cell-type relevance, an unbiased in vitro reporter screen for CREs overlapping all common variants associated with QT interval at the SCN5A GWAS locus identified 12 common cis-regulatory variants that map to cardiac open chromatin regions and correlate with SCN5A cardiac expression.
2019-04-22 | GSE109716 | GEO
Project description:Our goal is to identify high-confidence candidate genes associated with sub-threshold QT interval loci. We predicted enhancer-promoter targets using two different methods.
2014-10-10 | GSE53567 | GEO
Project description:Our goal is to identify high-confidence candidate genes associated with sub-threshold QT interval loci. We predicted enhancer-promoter targets using two different methods. 1) We used the correlation in activity between enhancers and transcribed genes across 59 human cell lines and tissues to identify significantly correlated enhancer-promoter pairs that represent potential regulatory interactions. 2) We used chromosome conformation capture followed by high-throughput sequencing (4C-seq) to probe physical enhancer-promoter interactions.
2014-10-10 | E-GEOD-53567 | 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) 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: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:MULTIPLE SCN5A ENHANCERS MODULATE CARDIAC GENE EXPRESSION AND QT INTERVAL VARIATION
| PRJNA431786 | ENA
Project description:Genome-wide association studies (GWAS) have identified twenty melanoma susceptibility loci. To identify susceptibility genes and variants simultaneously from multiple GWAS loci, we integrated massively-parallel reporter assays (MPRA) with melanocyte-specific expression quantitative trait locus (eQTL) profiling. We identified thirty-nine candidate functional variants displaying allelic transcriptional activity, of which nine from four loci were also correlated with local gene expression in melanocytes (CTSS, CASP8, MX2, and MAFF). Among these, we further characterized the locus in MX2 gene on chromosome band Chr21q22.3 and validated a functional variant, rs398206. The functional variant mediates allelic transcriptional activity via binding of the transcription factor, YY1. This allelic transcriptional regulation largely accounts for a significant cis-eQTL of the HIV-1 restriction gene, MX2, in primary melanocytes, where the melanoma risk-associated A allele is correlated with higher MX2 levels. Melanocyte-specific transgenic expression of human MX2 in a zebrafish model demonstrated an accelerated melanoma formation in a BRAFV600E background. Thus, using an efficient scalable approach to streamline GWAS follow-up functional studies, we uncovered a pleiotropic function of MX2 in melanoma susceptibility.
2020-04-20 | GSE129242 | GEO
Project description:Genome-wide association studies (GWAS) have identified twenty melanoma susceptibility loci. To identify susceptibility genes and variants simultaneously from multiple GWAS loci, we integrated massively-parallel reporter assays (MPRA) with melanocyte-specific expression quantitative trait locus (eQTL) profiling. We identified thirty-nine candidate functional variants displaying allelic transcriptional activity, of which nine from four loci were also correlated with local gene expression in melanocytes (CTSS, CASP8, MX2, and MAFF). Among these, we further characterized the locus in MX2 gene on chromosome band Chr21q22.3 and validated a functional variant, rs398206. The functional variant mediates allelic transcriptional activity via binding of the transcription factor, YY1. This allelic transcriptional regulation largely accounts for a significant cis-eQTL of the HIV-1 restriction gene, MX2, in primary melanocytes, where the melanoma risk-associated A allele is correlated with higher MX2 levels. Melanocyte-specific transgenic expression of human MX2 in a zebrafish model demonstrated an accelerated melanoma formation in a BRAFV600E background. Thus, using an efficient scalable approach to streamline GWAS follow-up functional studies, we uncovered a pleiotropic function of MX2 in melanoma susceptibility.
2020-04-20 | GSE129237 | GEO