ScDALI: modeling allelic heterogeneity in single cells reveals context-specific genetic regulation.
Ontology highlight
ABSTRACT:
SUBMITTER: Heinen T
PROVIDER: S-EPMC8734213 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Skeletal muscle, the largest human organ by weight, is relevant to several polygenic metabolic traits and diseases including type 2 diabetes (T2D). Identifying genetic mechanisms underlying these traits requires pinpointing the relevant cell types, regulatory elements, target genes, and causal variants. Here, we used genetic multiplexing to generate population-scale single nucleus (sn) chromatin accessibility (snATAC-seq) and transcriptome (snRNA-seq) maps across 287 frozen human skeletal muscle biopsies representing 456,880 nuclei. We identified 13 cell types that collectively represented 983,155 ATAC summits. We integrated genetic variation to discover 6,866 expression quantitative trait loci (eQTL) and 100,928 chromatin accessibility QTL (caQTL) (5% FDR) across the five most abundant cell types, cataloging caQTL peaks that atlas-level snATAC maps often miss. We identified 1,973 eGenes colocalized with caQTL and used mediation analyses to construct causal directional maps for chromatin accessibility and gene expression. 3,378 genome-wide association study (GWAS) signals across 43 relevant traits colocalized with sn-e/caQTL, 52% in a cell-specific manner. 77% of GWAS signals colocalized with caQTL and not eQTL, highlighting the critical importance of population-scale chromatin profiling for GWAS functional studies. GWAS-caQTL colocalization showed distinct cell-specific regulatory paradigms. For example, a C2CD4A/B T2D GWAS signal colocalized with caQTL in muscle fibers and multiple chromatin loop models nominated VPS13C, a glucose uptake gene. Sequence of the caQTL peak overlapping caSNP rs7163757 showed allelic regulatory activity differences in a human myocyte cell line massively parallel reporter assay. These results illuminate the genetic regulatory architecture of human skeletal muscle at high-resolution epigenomic, transcriptomic, and cell state scales and serve as a template for population-scale multi-omic mapping in complex tissues and traits.
| S-EPMC10760134 | biostudies-literature
Project description:The host's gene expression and gene regulatory response to pathogen exposure can be influenced by a combination of the host's genetic background, the type of and exposure time to pathogens. Here we provide a detailed dissection of this using single-cell RNA-sequencing of 1.3M peripheral blood mononuclear cells from 120 individuals, longitudinally exposed to three different pathogens. These analyses indicate that cell-type-specificity is a more prominent factor than pathogen-specificity regarding contexts that affect how genetics influences gene expression (i.e., eQTL) and co-expression (i.e., co-expression QTL). In monocytes, the strongest responder to pathogen stimulations, 71.4% of the genetic variants whose effect on gene expression is influenced by pathogen exposure (i.e., response QTL) also affect the co-expression between genes. This indicates widespread, context-specific changes in gene expression level and its regulation that are driven by genetics. Pathway analysis on the CLEC12A gene that exemplifies cell-type-, exposure-time- and genetic-background-dependent co-expression interactions, shows enrichment of the interferon (IFN) pathway specifically at 3-h post-exposure in monocytes. Similar genetic background-dependent association between IFN activity and CLEC12A co-expression patterns is confirmed in systemic lupus erythematosus by in silico analysis, which implies that CLEC12A might be an IFN-regulated gene. Altogether, this study highlights the importance of context for gaining a better understanding of the mechanisms of gene regulation in health and disease.
| S-EPMC9174272 | biostudies-literature
Project description:Modern genetic mapping is plagued by the "missing heritability" problem, which refers to the discordance between the estimated heritabilities of quantitative traits and the variance accounted for by mapped causative variants. One major potential explanation for the missing heritability is allelic heterogeneity, in which there are multiple causative variants at each causative gene with only a fraction having been identified. The majority of genome-wide association studies (GWAS) implicitly assume that a single SNP can explain all the variance for a causative locus. However, if allelic heterogeneity is prevalent, a substantial amount of genetic variance will remain unexplained. In this paper, we take a haplotype-based mapping approach and quantify the number of alleles segregating at each locus using a large set of 7922 eQTL contributing to regulatory variation in the Drosophila melanogaster female head. Not only does this study provide a comprehensive eQTL map for a major community genetic resource, the Drosophila Synthetic Population Resource, but it also provides a direct test of the allelic heterogeneity hypothesis. We find that 95% of cis-eQTLs and 78% of trans-eQTLs are due to multiple alleles, demonstrating that allelic heterogeneity is widespread in Drosophila eQTL. Allelic heterogeneity likely contributes significantly to the missing heritability problem common in GWAS studies.
| S-EPMC4014434 | biostudies-literature
Project description:Many genes are known to regulate retinal regeneration after widespread tissue damage. Conversely, genes controlling regeneration after limited cell loss, as per degenerative diseases, are undefined. As stem/progenitor cell responses scale to injury levels, understanding how the extent and specificity of cell loss impact regenerative processes is important. Here, transgenic zebrafish enabling selective retinal ganglion cell (RGC) ablation were used to identify genes that regulate RGC regeneration. A single cell multiomics-informed screen of 100 genes identified seven knockouts that inhibited and 11 that promoted RGC regeneration. Surprisingly, 35 out of 36 genes known and/or implicated as being required for regeneration after widespread retinal damage were not required for RGC regeneration. The loss of seven even enhanced regeneration kinetics, including the proneural factors neurog1, olig2 and ascl1a. Mechanistic analyses revealed that ascl1a disruption increased the propensity of progenitor cells to produce RGCs, i.e. increased 'fate bias'. These data demonstrate plasticity in the mechanism through which Müller glia convert to a stem-like state and context specificity in how genes function during regeneration. Increased understanding of how the regeneration of disease-relevant cell types is specifically controlled will support the development of disease-tailored regenerative therapeutics.
| S-EPMC11361637 | biostudies-literature
Project description:Measuring multiple omics profiles from the same single cell opens up the opportunity to decode molecular regulation that underlies intercellular heterogeneity in development and disease. Here, we present co-sequencing of microRNAs and mRNAs in the same single cell using a half-cell genomics approach. This method demonstrates good robustness (~95% success rate) and reproducibility (R2 = 0.93 for both microRNAs and mRNAs), yielding paired half-cell microRNA and mRNA profiles, which we can independently validate. By linking the level of microRNAs to the expression of predicted target mRNAs across 19 single cells that are phenotypically identical, we observe that the predicted targets are significantly anti-correlated with the variation of abundantly expressed microRNAs. This suggests that microRNA expression variability alone may lead to non-genetic cell-to-cell heterogeneity. Genome-scale analysis of paired microRNA-mRNA co-profiles further allows us to derive and validate regulatory relationships of cellular pathways controlling microRNA expression and intercellular variability.
| S-EPMC6327095 | biostudies-literature
Project description:To determine the dynamics of allelic-specific expression during mouse development, we analyzed RNA-seq data from 23 F1 tissues from different developmental stages, including 19 female tissues allowing X chromosome inactivation (XCI) escapers to also be detected. We demonstrate that allelic expression arising from genetic or epigenetic differences is highly tissue-specific. We find that tissue-specific strain-biased gene expression may be regulated by tissue-specific enhancers or by post-transcriptional differences in stability between the alleles. We also find that escape from X-inactivation is tissue-specific, with leg muscle showing an unexpectedly high rate of XCI escapers. By surveying a range of tissues during development, and performing extensive validation, we are able to provide a high confidence list of mouse imprinted genes including 18 novel genes. This shows that cluster size varies dynamically during development and can be substantially larger than previously thought, with the Igf2r cluster extending over 10 Mb in placenta.
| S-EPMC5555720 | biostudies-literature
Project description:A functional allele of the mouse catechol-O-methyltransferase (Comt) gene is defined by the insertion of a B2 short interspersed repeat element in its 3'-untranslated region (UTR). This allele has been associated with a number of phenotypes, such as pain and anxiety. In comparison with mice carrying the ancestral allele (Comt+), Comt B2i mice show higher Comt mRNA and enzymatic activity levels. Here, we investigated the molecular genetic mechanisms underlying this allelic specific regulation of Comt expression. Insertion of the B2 element introduces an early polyadenylation signal generating a shorter Comt transcript, in addition to the longer ancestral mRNA. Comparative analysis and in silico prediction of Comt mRNA potential targets within the transcript 3' to the B2 element was performed and allowed choosing microRNA (miRNA) candidates for experimental screening: mmu-miR-3470a, mmu-miR-3470b, and mmu-miR-667. Cell transfection with each miRNA downregulated the expression of the ancestral transcript and COMT enzymatic activity. Our in vivo experiments showed that mmu-miR-667-3p is strongly correlated with decreasing amounts of Comt mRNA in the brain, and lentiviral injections of mmu-miR-3470a, mmu-miR-3470b, and mmu-miR-667 increase hypersensitivity in the mouse formalin model, consistent with reduced COMT activity. In summary, our data demonstrate that the Comt+ transcript contains regulatory miRNA signals in its 3'-untranslated region leading to mRNA degradation; these signals, however, are absent in the shorter transcript, resulting in higher mRNA expression and activity levels.
| S-EPMC4579042 | biostudies-literature
Project description:A majority of the variants identified in genome-wide association studies fall in non-coding regions of the genome, indicating their mechanism of impact is mediated via gene expression. Leveraging this hypothesis, transcriptome-wide association studies (TWAS) have assisted in both the interpretation and discovery of additional genes associated with complex traits. However, existing methods for conducting TWAS do not take full advantage of the intra-individual correlation inherently present in multi-context expression studies and do not properly adjust for multiple testing across contexts. We introduce CONTENT—a computationally efficient method with proper cross-context false discovery correction that leverages correlation structure across contexts to improve power and generate context-specific and context-shared components of expression. We apply CONTENT to bulk multi-tissue and single-cell RNA-seq data sets and show that CONTENT leads to a 42% (bulk) and 110% (single cell) increase in the number of genetically predicted genes relative to previous approaches. We find the context-specific component of expression comprises 30% of heritability in tissue-level bulk data and 75% in single-cell data, consistent with cell-type heterogeneity in bulk tissue. In the context of TWAS, CONTENT increases the number of locus-phenotype associations discovered by over 51% relative to previous methods across 22 complex traits. Transcriptome-wide association studies (TWAS) often ignore the specificity and sharing of effects across contexts (e.g., tissues). Here, the authors describe a method to split genetic effects into context-shared and context-specific terms. They apply their method to tissue and single-cell RNA-seq and show improved power in TWAS.
| S-EPMC9519579 | biostudies-literature
Project description:Uncovering the functional impact of genetic variation on gene expression is important in understanding tissue biology and the pathogenesis of complex traits. Despite large efforts to map expression quantitative trait loci (eQTLs) across many human tissues, our ability to translate those findings to understanding human disease has been incomplete, and the majority of disease loci are not explained by association with expression of a target gene. Cell-type specificity and the presence of multiple independent causal variants for many eQTLs are potential confounders contributing to the apparent discrepancy with disease loci. In this study, we investigate the tissue specificity of genetic effects on gene expression and the overlap with disease loci while considering the presence of multiple causal variants within and across tissues. We find evidence of pervasive tissue specificity of eQTLs, often masked by linkage disequilibrium that misleads traditional meta-analytic approaches. We propose CAFEH (colocalization and fine-mapping in the presence of allelic heterogeneity), a Bayesian method that integrates genetic association data across multiple traits, incorporating linkage disequilibrium to identify causal variants. CAFEH outperforms previous approaches in colocalization and fine-mapping. Using CAFEH, we show that genes with highly tissue-specific genetic effects are under greater selection, enriched in differentiation and developmental processes, and more likely to be involved in human disease. Last, we demonstrate that CAFEH can efficiently leverage the widespread allelic heterogeneity in genetic regulation of gene expression to prioritize the target tissue in genome-wide association complex trait loci, thereby improving our ability to interpret complex trait genetics.
| S-EPMC8874223 | biostudies-literature
Project description:Genetic regulation of gene expression is a complex process, with genetic effects known to vary across cellular contexts such as cell types and environmental conditions. We developed SURGE, a method for unsupervised discovery of context-specific expression quantitative trait loci (eQTLs) from single-cell transcriptomic data. This allows discovery of the contexts or cell types modulating genetic regulation without prior knowledge. Applied to peripheral blood single-cell eQTL data, SURGE contexts capture continuous representations of distinct cell types and groupings of biologically related cell types. We demonstrate the disease-relevance of SURGE context-specific eQTLs using colocalization analysis and stratified LD-score regression.
| S-EPMC10801966 | biostudies-literature