Dataset Information

Massively parallel reporter assays with multi-layer annotation identified cell-type-specific functional variants and genes associated with melanoma (RNA-Seq)

ABSTRACT: 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.

ORGANISM(S): Homo sapiens

PROVIDER: GSE210355 | GEO | 2022/11/28

REPOSITORIES: GEO

ACCESS DATA

Similar Datasets

Project description:Background: The identification of causal variants responsible for disease associations from genome-wide association studies (GWAS) facilitates functional understanding of the disease mechanisms implicated by GWAS. One of the earliest GWAS associations to COPD spans an intragenic region within FAM13A, but the causal variants at this loci have not yet been identified. Massively parallel reporter assays (MPRA) can be used to prioritize functional regulatory variants in a high-throughput manner. Methods: We used an integrated approach using fine-mapping in over 10,000 subjects from COPD GWAS studies, two MPRA experiments, traditional reporter assays, chromatin conformation capture, and CRISPR-based gene editing to characterize COPD-associated regulatory variants in FAM13A in human bronchial epithelial cell lines. Results: Conditional genetic association and fine mapping analyses identified two independent COPD association signals in FAM13A. MPRA identified 45 common functional regulatory variants, and six COPD-associated putative functional variants were prioritized for further functional investigation. Three variants demonstrated significant activity in traditional reporter assays, and one variant, rs2013701, was selected for further testing in the endogenous genomic context based on a direction of effect consistent with postulated mechanisms of FAM13A-mediated COPD susceptibility. CRISPR-based genome editing for this variant confirmed allele-specific effects on FAM13A expression and altered rates of cellular proliferation, providing multiple levels of functional characterization for this COPD-associated variant. Conclusions: Comprehensive screening for regulatory variants near FAM13A identified the presence of extensive functional regulatory variation within a 250kb window of FAM13A in HBECs. Focused functional evaluation of the COPD-associated functional variants in LD with the two independent association signals in this region prioritized the common variant rs2013701, for which multiple parallel lines of functional evidence confirm allelic effects on FAM13A regulation.

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.

Project description:Recent genome-wide association studies (GWAS) have identified a number of novel genetic associations with complex human diseases. In spite of these successes, results from GWAS generally explain only a small proportion of disease heritability, an observation termed the M-bM-^@M-^\missing heritability problem.M-bM-^@M-^] Several sources for the missing heritability have been proposed, including the contribution of many common variants with small individual effect sizes, which cannot be reliably found using the standard GWAS approach. The goal of our study was to explore a complementary approach, which combines GWAS results with functional data in order to identify novel genetic associations with small effect sizes. To do so, we conducted a GWAS for lymphocyte count, a physiologic quantitative trait associated with asthma, in 462 Hutterites. In parallel, we performed a genome-wide gene expression study in lymphoblastoid cell lines (LCLs) from 96 Hutterites. We found significant support for genetic associations using the GWAS data when we considered variants near the 193 genes whose expression levels across individuals were most correlated with lymphocyte counts. Interestingly, these variants are also enriched with signatures of an association with asthma susceptibility, an observation we were able to replicate. The associated loci include genes previously implicated in asthma susceptibility, as well as novel candidate genes enriched for functions related to T cell receptor signaling and ATP synthesis. Our results, therefore, establish a new set of asthma susceptibility candidate genes. More generally, our observations support the notion that many loci of small effects influence variation in lymphocyte count and asthma susceptibility. 96 RNA samples were collected (1 subsequently excluded) from lymphoblastoid cell lines derived from Hutterite subjects chosen to represent the extremes of absolute lymphocyte count. The high and low absolute lymphocyte count groups were balanced with respect to age, gender, and relatedness.

Dataset Information

Massively parallel reporter assays with multi-layer annotation identified cell-type-specific functional variants and genes associated with melanoma (RNA-Seq)

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure