Dataset Information

Single-cell genomics improves the discovery of risk variants and genes of Atrial Fibrillation [scATAC-seq]

ABSTRACT: Genome-wide association studies (GWAS) have linked hundreds of loci to cardiac diseases. However, in most loci the causal variants and their target genes remain unknown. We developed a combined experimental and analytical approach that integrates single cell epigenomics with GWAS to prioritize risk variants and genes. We profiled accessible chromatin in single cells obtained from human hearts and leveraged the data to study genetics of Atrial Fibrillation (AF), the most common cardiac arrhythmia. Enrichment analysis of AF risk variants using cell-type-resolved open chromatin regions (OCRs) implicated cardiomyocytes as the main mediator of AF risk. We then performed statistical fine-mapping, leveraging the information in OCRs, and identified putative causal variants in 122 AF-associated loci. Taking advantage of the fine-mapping results, our novel statistical procedure for gene discovery prioritized 45 high-confidence risk genes, highlighting transcription factors and signal transduction pathways important for heart development. We further leveraged our single-cell data to study genetics of gene expression. An unexpected finding from earlier studies is that expression QTLs (eQTLs) are often shared across tissues even though most regulatory elements are cell-type specific. We found that this sharing is largely driven by the limited power of eQTL studies using bulk tissues to detect cell-type-specific regulatory variants. This finding points to an important limitation of using eQTLs to interpret GWAS of complex traits. In summary, our analysis provides a comprehensive map of AF risk variants and genes, and a general framework to integrate single-cell genomics with genetic studies of complex traits.

ORGANISM(S): Homo sapiens

PROVIDER: GSE224996 | GEO | 2023/02/12

REPOSITORIES: GEO

ACCESS DATA

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.

Project description:Most loci identified in genome wide association studies (GWAS) of complex traits reside in non-coding DNA and may contribute to phenotype via changes in gene regulation. The discovery of expression quantitative trait loci (?eQTLs?) can thus be used to more precisely identify modest but real disease associations and provide insights into their underlying molecular mechanisms. This is particularly true for analyses of expression in non-transformed cells from tissues relevant to the complex traits of interest. We have conducted two independent studies to identify genetic, including both SNPs and copy-number variants, and environmental determinants of human liver gene expression variation. We analyzed two sets of primary livers (primary dataset: n=220; replication dataset: n=60) using Agilent and Illumina expression arrays and Illumina SNP genotyping (550K). At least 30% of genetic and non-genetic factors that meet genome-wide significance (p <1 x10-9) in one study fail to replicate in the second study, suggesting that artifacts, like unknown SNPs that affect RNA-probe hybridization or hidden confounding variables, often result in statistically significant but biologically irrelevant correlations. These data confirm the value of independent replications to enrich for truly predictive eQTLs, and given our study design we are able to identify hundreds of reproducible correlations. We show that such information can be used to provide insights into disease-relevant phenotypes, with specific examples including eQTLs related to lipid levels (e.g. LDL cholesterol), immune system function (e.g. HLA), and drug response (e.g. warfarin). Furthermore, in the interest of both fine-mapping and mechanistic annotation, we hypothesized that promoters and 3?UTRs are enriched for causal eQTL variants. Therefore, we re-sequenced the promoter and 3?UTR regions of 25 genes with eQTLs, cloned each discovered haplotype, and quantified their impact on transcription using a luciferase-based assay. These data reveal multiple examples of robust, haplotype-specific in vitro functional differences that correlate directly with in vivo expression levels. This suggests that many eQTLs can be rapidly fine-mapped to one or a few single-nucleotide variants and mechanistically characterized using such assays. Integration of functional assays with eQTL discovery, and eQTLs with complex trait associations, is a powerful means to exploit GWAS data and improve their biological interpretability. RNA expression levels were quantified on Illumina gene expression microarrays for 60 normal human livers. Expression quantitative trait loci were identified by genome wide association mapping.

Project description:Most loci identified in genome wide association studies (GWAS) of complex traits reside in non-coding DNA and may contribute to phenotype via changes in gene regulation. The discovery of expression quantitative trait loci (‘eQTLs’) can thus be used to more precisely identify modest but real disease associations and provide insights into their underlying molecular mechanisms. This is particularly true for analyses of expression in non-transformed cells from tissues relevant to the complex traits of interest. We have conducted two independent studies to identify genetic, including both SNPs and copy-number variants, and environmental determinants of human liver gene expression variation. We analyzed two sets of primary livers (primary dataset: n=220; replication dataset: n=60) using Agilent and Illumina expression arrays and Illumina SNP genotyping (550K). At least 30% of genetic and non-genetic factors that meet genome-wide significance (p <1 x10-9) in one study fail to replicate in the second study, suggesting that artifacts, like unknown SNPs that affect RNA-probe hybridization or hidden confounding variables, often result in statistically significant but biologically irrelevant correlations. These data confirm the value of independent replications to enrich for truly predictive eQTLs, and given our study design we are able to identify hundreds of reproducible correlations. We show that such information can be used to provide insights into disease-relevant phenotypes, with specific examples including eQTLs related to lipid levels (e.g. LDL cholesterol), immune system function (e.g. HLA), and drug response (e.g. warfarin). Furthermore, in the interest of both fine-mapping and mechanistic annotation, we hypothesized that promoters and 3’UTRs are enriched for causal eQTL variants. Therefore, we re-sequenced the promoter and 3’UTR regions of 25 genes with eQTLs, cloned each discovered haplotype, and quantified their impact on transcription using a luciferase-based assay. These data reveal multiple examples of robust, haplotype-specific in vitro functional differences that correlate directly with in vivo expression levels. This suggests that many eQTLs can be rapidly fine-mapped to one or a few single-nucleotide variants and mechanistically characterized using such assays. Integration of functional assays with eQTL discovery, and eQTLs with complex trait associations, is a powerful means to exploit GWAS data and improve their biological interpretability. RNA expression levels were quantified on Agilent gene expression microarrays for 206 normal human livers (464 unique arrays corresponding to those samples). Genotypes were also derived from 224 of these samples using Illumina SNP chips. Expression quantitative trait loci were identified by genome wide association mapping. This submission represents the primary transcriptome component of study.

Dataset Information

Single-cell genomics improves the discovery of risk variants and genes of Atrial Fibrillation [scATAC-seq]

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure