scGWAS: landscape of trait-cell type associations by integrating single-cell transcriptomics-wide and genome-wide association studies

Harnessing omics data for drug discovery and development in ovarian aging

BACKGROUND

Ovarian aging occurs earlier than the aging of many other organs and has a lasting impact on women's overall health and well-being. However, effective interventions to slow ovarian aging remain limited, primarily due to an incomplete understanding of the underlying molecular mechanisms and drug targets. Recent advances in omics data resources, combined with innovative computational tools, are offering deeper insight into the molecular complexities of ovarian aging, paving the way for new opportunities in drug discovery and development.

OBJECTIVE AND RATIONALE

This review aims to synthesize the expanding multi-omics data, spanning genome, transcriptome, proteome, metabolome, and microbiome, related to ovarian aging, from both tissue-level and single-cell perspectives. We will specially explore how the analysis of these emerging omics datasets can be leveraged to identify novel drug targets and guide therapeutic strategies for slowing and reversing ovarian aging.

SEARCH METHODS

We conducted a comprehensive literature search in the PubMed database using a range of relevant keywords: ovarian aging, age at natural menopause, premature ovarian insufficiency (POI), diminished ovarian reserve (DOR), genomics, transcriptomics, epigenomics, DNA methylation, RNA modification, histone modification, proteomics, metabolomics, lipidomics, microbiome, single-cell, genome-wide association studies (GWAS), whole-exome sequencing, phenome-wide association studies (PheWAS), Mendelian randomization (MR), epigenetic target, drug target, machine learning, artificial intelligence (AI), deep learning, and multi-omics. The search was restricted to English-language articles published up to September 2024.

OUTCOMES

Multi-omics studies have uncovered key mechanisms driving ovarian aging, including DNA damage and repair deficiencies, inflammatory and immune responses, mitochondrial dysfunction, and cell death. By integrating multi-omics data, researchers can identify critical regulatory factors and mechanisms across various biological levels, leading to the discovery of potential drug targets. Notable examples include genetic targets such as BRCA2 and TERT, epigenetic targets like Tet and FTO, metabolic targets such as sirtuins and CD38+, protein targets like BIN2 and PDGF-BB, and transcription factors such as FOXP1.

WIDER IMPLICATIONS

The advent of cutting-edge omics technologies, especially single-cell technologies and spatial transcriptomics, has provided valuable insights for guiding treatment decisions and has become a powerful tool in drug discovery aimed at mitigating or reversing ovarian aging. As technology advances, the integration of single-cell multi-omics data with AI models holds the potential to more accurately predict candidate drug targets. This convergence offers promising new avenues for personalized medicine and precision therapies, paving the way for tailored interventions in ovarian aging.

REGISTRATION NUMBER

Not applicable.

Digital twins as global learning health and disease models for preventive and personalized medicine

Ineffective medication is a major healthcare problem causing significant patient suffering and economic costs. This issue stems from the complex nature of diseases, which involve altered interactions among thousands of genes across multiple cell types and organs. Disease progression can vary between patients and over time, influenced by genetic and environmental factors. To address this challenge, digital twins have emerged as a promising approach, which have led to international initiatives aiming at clinical implementations. Digital twins are virtual representations of health and disease processes that can integrate real-time data and simulations to predict, prevent, and personalize treatments. Early clinical applications of DTs have shown potential in areas like artificial organs, cancer, cardiology, and hospital workflow optimization. However, widespread implementation faces several challenges: (1) characterizing dynamic molecular changes across multiple biological scales; (2) developing computational methods to integrate data into DTs; (3) prioritizing disease mechanisms and therapeutic targets; (4) creating interoperable DT systems that can learn from each other; (5) designing user-friendly interfaces for patients and clinicians; (6) scaling DT technology globally for equitable healthcare access; (7) addressing ethical, regulatory, and financial considerations. Overcoming these hurdles could pave the way for more predictive, preventive, and personalized medicine, potentially transforming healthcare delivery and improving patient outcomes.

MAAT: a new nonparametric Bayesian framework for incorporating multiple functional annotations in transcriptome-wide association studies

Transcriptome-wide association study (TWAS) has emerged as a powerful tool for translating the myriad variations identified by genome-wide association studies (GWAS) into regulated genes in the post-GWAS era. While integrating annotation information has been shown to enhance power, current annotation-assisted TWAS tools predominantly focus on epigenomic annotations. When including more annotations, the assumption of a positive correlation between annotation scores and SNPs' effect sizes, as adopted by current methods, often falls short. Here, we propose MAAT expanding the horizons of existing TWAS studies, generating a new model incorporating multiple annotations into TWAS and a new metric indicating the most important annotation.

sc2GWAS: a comprehensive platform linking single cell and GWAS traits of human

Identifying cell populations associated with risk variants is essential for uncovering cell-specific mechanisms that drive disease development and progression. Integrating genome-wide association studies (GWAS) with single-cell RNA sequencing (scRNA-seq) has become an effective strategy for detecting trait-cell relationships. The accumulation of trait-related single cell data has led to an urgent need for its comprehensively processing. To address this, we developed sc2GWAS (https://bio.liclab.net/sc2GWAS/), which aims to document large-scale GWAS trait-cell regulatory pairs at single-cell resolution and provide comprehensive annotations and enrichment analyses for these related pairs. The current version of sc2GWAS curates a total of 15 078 310 candidate trait-cell pairs from > 6 300 000 individual cells, offering a valuable resource for exploring complex regulatory relationships between traits and cells. We applied strict quality control measures on both scRNA-seq data and GWAS data, ensuring the reliability and accuracy of the datasets for the identification of trait-relevant cells and genes. In addition, sc2GWAS provides ranked lists of trait-relevant genes and extensive (epi) genetic annotations, making it a valuable resource for downstream analyses. We demonstrate the utility of the platform by investigating Alzheimer's disease, where we identified significant associations between the disease and microglial cells, with the APOE gene emerging as particularly significant. This platform facilitates detailed research into complex trait-cell and trait-gene interactions, we anticipate that sc2GWAS will become a comprehensive and valuable platform for exploring GWAS trait-cell regulatory mechanisms.

Evaluating methods for integrating single-cell data and genetics to understand inflammatory disease complexity

Background

Understanding genetic underpinnings of immune-mediated inflammatory diseases is crucial to improve treatments. Single-cell RNA sequencing (scRNA-seq) identifies cell states expanded in disease, but often overlooks genetic causality due to cost and small genotyping cohorts. Conversely, large genome-wide association studies (GWAS) are commonly accessible.

Methods

We present a 3-step robust benchmarking analysis of integrating GWAS and scRNA-seq to identify genetically relevant cell states and genes in inflammatory diseases. First, we applied and compared the results of three recent algorithms, based on pathways (scGWAS), single-cell disease scores (scDRS), or both (scPagwas), according to accuracy/sensitivity and interpretability. While previous studies focused on coarse cell types, we used disease-specific, fine-grained single-cell atlases (183,742 and 228,211 cells) and GWAS data (Ns of 97,173 and 45,975) for rheumatoid arthritis (RA) and ulcerative colitis (UC). Second, given the lack of scRNA-seq for many diseases with GWAS, we further tested the tools' resolution limits by differentiating between similar diseases with only one fine-grained scRNA-seq atlas. Lastly, we provide a novel evaluation of noncoding SNP incorporation methods by testing which enabled the highest sensitivity/accuracy of known cell-state calls.

Results

We first found that single-cell based tools scDRS and scPagwas called superior numbers of supported cell states that were overlooked by scGWAS. While scGWAS and scPagwas were advantageous for gene exploration, scDRS effectively accounted for batch effect and captured cellular heterogeneity of disease-relevance without single-cell genotyping. For noncoding SNP integration, we found a key trade-off between statistical power and confidence with positional (e.g. MAGMA) and non-positional approaches (e.g. chromatin-interaction, eQTL). Even when directly incorporating noncoding SNPs through 5' scRNA-seq measures of regulatory elements, non disease-specific atlases gave misleading results by not containing disease-tissue specific transcriptomic patterns. Despite this criticality of tissue-specific scRNA-seq, we showed that scDRS enabled deconvolution of two similar diseases with a single fine-grained scRNA-seq atlas and separate GWAS. Indeed, we identified supported and novel genetic-phenotype linkages separating RA and ankylosing spondylitis, and UC and crohn's disease. Overall, while noting evolving single-cell technologies, our study provides key findings for integrating expanding fine-grained scRNA-seq, GWAS, and noncoding SNP resources to unravel the complexities of inflammatory diseases.

Differential expression and co-expression reveal cell types relevant to genetic disorder phenotypes

MOTIVATION

Knowledge of the specific cell types affected by genetic alterations in rare diseases is crucial for advancing diagnostics and treatments. Despite significant progress, the cell types involved in the majority of rare disease manifestations remain largely unknown. In this study, we integrated scRNA-seq data from non-diseased samples with known genetic disorder genes and phenotypic information to predict the specific cell types disrupted by pathogenic mutations for 482 disease phenotypes.

RESULTS

We found significant phenotype-cell type associations focusing on differential expression and co-expression mechanisms. Our analysis revealed that 13% of the associations documented in the literature were captured through differential expression, while 42% were elucidated through co-expression analysis, also uncovering potential new associations. These findings underscore the critical role of cellular context in disease manifestation and highlight the potential of single-cell data for the development of cell-aware diagnostics and targeted therapies for rare diseases.

AVAILABILITY AND IMPLEMENTATION

All code generated in this work is available at https://github.com/SergioAlias/sc-coex.

A genome-wide meta-analysis of palmoplantar pustulosis implicates TH2 responses and cigarette smoking in disease pathogenesis

BACKGROUND

Palmoplantar pustulosis (PPP) is an inflammatory skin disorder that mostly affects smokers and manifests with painful pustular eruptions on the palms and soles. Although the disease can present with concurrent plaque psoriasis, TNF and IL-17/IL-23 inhibitors show limited efficacy. There is therefore a pressing need to uncover PPP disease drivers and therapeutic targets.

OBJECTIVES

We sought to identify genetic determinants of PPP and investigate whether cigarette smoking contributes to disease pathogenesis.

METHODS

We performed a genome-wide association meta-analysis of 3 North-European cohorts (n = 1,456 PPP cases and 402,050 controls). We then used the scGWAS program to investigate the cell-type specificity of the association signals. We also undertook genetic correlation analyses to examine the similarities between PPP and other immune-mediated diseases. Finally, we applied Mendelian randomization to analyze the causal relationship between cigarette smoking and PPP.

RESULTS

We found that PPP is not associated with the main genetic determinants of plaque psoriasis. Conversely, we identified genome-wide significant associations with the FCGR3A/FCGR3B and CCHCR1 loci. We also observed 13 suggestive (P < 5 × 10-6) susceptibility regions, including the IL4/IL13 interval. Accordingly, we demonstrated a significant genetic correlation between PPP and TH2-mediated diseases such as atopic dermatitis and ulcerative colitis. We also found that genes mapping to PPP-associated intervals were preferentially expressed in dendritic cells and often implicated in T-cell activation pathways. Finally, we undertook a Mendelian randomization analysis, which supported a causal role of cigarette smoking in PPP.

CONCLUSIONS

The first genome-wide association study of PPP points to a pathogenic role for deregulated TH2 responses and cigarette smoking.

Gene-environment interactions within a precision environmental health framework

Understanding the complex interplay of genetic and environmental factors in disease etiology and the role of gene-environment interactions (GEIs) across human development stages is important. We review the state of GEI research, including challenges in measuring environmental factors and advantages of GEI analysis in understanding disease mechanisms. We discuss the evolution of GEI studies from candidate gene-environment studies to genome-wide interaction studies (GWISs) and the role of multi-omics in mediating GEI effects. We review advancements in GEI analysis methods and the importance of large-scale datasets. We also address the translation of GEI findings into precision environmental health (PEH), showcasing real-world applications in healthcare and disease prevention. Additionally, we highlight societal considerations in GEI research, including environmental justice, the return of results to participants, and data privacy. Overall, we underscore the significance of GEI for disease prediction and prevention and advocate for integrating the exposome into PEH omics studies.

Prioritization of Kidney Cell Types Highlights Myofibroblast Cells in Regulating Human Blood Pressure

Introduction

Blood pressure (BP) is a highly heritable trait with over 2000 underlying genomic loci identified to date. Although the kidney plays a key role, little is known about specific cell types involved in the genetic regulation of BP.

Methods

Here, we applied stratified linkage disequilibrium score (LDSC) regression to connect BP genome-wide association studies (GWAS) results to specific cell types of the mature human kidney. We used the largest single-stage BP genome-wide analysis to date, including up to 1,028,980 adults of European ancestry, and single-cell transcriptomic data from 14 mature human kidneys, with mean age of 41 years.

Results

Our analyses prioritized myofibroblasts and endothelial cells, among the total of 33 annotated cell type, as specifically involved in BP regulation (P < 0.05/33, i.e., 0.001515). Enrichment of heritability for systolic BP (SBP) was observed in myofibroblast cells in mature human kidney cortex, and enrichment of heritability for diastolic BP (DBP) was observed in descending vasa recta and peritubular capillary endothelial cells as well as stromal myofibroblast cells. The new finding of myofibroblast, the significant cell type for both BP traits, was consistent in 8 replication efforts using 7 sets of independent data, including in human fetal kidney, in East-Asian (EAS) ancestry, using mouse single-cell RNA sequencing (scRNA-seq) data, and when using another prioritization method.

Conclusion

Our findings provide a solid basis for follow-up studies to further identify genes and mechanisms in myofibroblast cells that underlie the regulation of BP.

Disentangling associations between complex traits and cell types with seismic

Integrating single-cell RNA sequencing (scRNA-seq) with Genome-Wide Association Studies (GWAS) can help reveal GWAS-associated cell types, furthering our understanding of the cell-type-specific biological processes underlying complex traits and disease. However, current methods have technical limitations that hinder them from making systematic, scalable, interpretable disease-cell-type associations. In order to rapidly and accurately pinpoint associations, we develop a novel framework, seismic, which characterizes cell types using a new specificity score. We compare seismic with alternative methods across over 1,000 cell type characterizations at different granularities and 28 traits, demonstrating that seismic both corroborates findings and identifies trait-relevant cell groups which are not apparent through other methodologies. Furthermore, as part of the seismic framework, the specific genes driving cell type-trait associations can easily be accessed and analyzed, enabling further biological insights. The advantages of seismic are particularly salient in neurodegenerative diseases such as Parkinson's and Alzheimer's, where disease pathology has not only cell-specific manifestations, but also brain region-specific differences. Interestingly, a case study of Alzheimer's disease reveals the importance of considering GWAS endpoints, as studies relying on clinical diagnoses consistently identify microglial associations, while GWAS with a tau biomarker endpoint reveals neuronal associations. In general, seismic is a computationally efficient, powerful, and interpretable approach for identifying associations between complex traits and cell type-specific expression.

Integrating leiomyoma genetics, epigenomics, and single-cell transcriptomics reveals causal genetic variants, genes, and cell types

Uterine fibroids (UF), that can disrupt normal uterine function and cause significant physical and psychological health problems, are observed in nearly 70% of women of reproductive age. Although heritable genetics is a significant risk factor, specific genetic variations and gene targets causally associated with UF are poorly understood. Here, we performed a meta-analysis on existing fibroid genome-wide association studies (GWAS) and integrated the identified risk loci and potentially causal single nucleotide polymorphisms (SNPs) with epigenomics, transcriptomics, 3D chromatin organization from diverse cell types as well as primary UF patient's samples. This integrative analysis identifies 24 UF-associated risk loci that potentially target 394 genes, of which 168 are differentially expressed in UF tumors. Critically, integrating this data with single-cell gene expression data from UF patients reveales the causal cell types with aberrant expression of these target genes. Lastly, CRISPR-based epigenetic repression (dCas9-KRAB) or activation (dCas9-p300) in a UF disease-relevant cell type further refines and narrows down the potential gene targets. Our findings and the methodological approach indicate the effectiveness of integrating multi-omics data with locus-specific epigenetic editing approaches for identifying gene- and celt type-targets of disease-relevant risk loci.

CellSTAR: a comprehensive resource for single-cell transcriptomic annotation

Large-scale studies of single-cell sequencing and biological experiments have successfully revealed expression patterns that distinguish different cell types in tissues, emphasizing the importance of studying cellular heterogeneity and accurately annotating cell types. Analysis of gene expression profiles in these experiments provides two essential types of data for cell type annotation: annotated references and canonical markers. In this study, the first comprehensive database of single-cell transcriptomic annotation resource (CellSTAR) was thus developed. It is unique in (a) offering the comprehensive expertly annotated reference data for annotating hundreds of cell types for the first time and (b) enabling the collective consideration of reference data and marker genes by incorporating tens of thousands of markers. Given its unique features, CellSTAR is expected to attract broad research interests from the technological innovations in single-cell transcriptomics, the studies of cellular heterogeneity & dynamics, and so on. It is now publicly accessible without any login requirement at: https://idrblab.org/cellstar.

Dynamic regulatory elements in single-cell multimodal data implicate key immune cell states enriched for autoimmune disease heritability

In autoimmune diseases such as rheumatoid arthritis, the immune system attacks the body's own cells. Developing a precise understanding of the cell states where noncoding autoimmune risk variants impart causal mechanisms is critical to developing curative therapies. Here, to identify noncoding regions with accessible chromatin that associate with cell-state-defining gene expression patterns, we leveraged multimodal single-nucleus RNA and assay for transposase-accessible chromatin (ATAC) sequencing data across 28,674 cells from the inflamed synovial tissue of 12 donors. Specifically, we used a multivariate Poisson model to predict peak accessibility from single-nucleus RNA sequencing principal components. For 14 autoimmune diseases, we discovered that cell-state-dependent ('dynamic') chromatin accessibility peaks in immune cell types were enriched for heritability, compared with cell-state-invariant ('cs-invariant') peaks. These dynamic peaks marked regulatory elements associated with T peripheral helper, regulatory T, dendritic and STAT1+CXCL10+ myeloid cell states. We argue that dynamic regulatory elements can help identify precise cell states enriched for disease-critical genetic variation.

Craniofacial developmental biology in the single-cell era

The evolution of a unique craniofacial complex in vertebrates made possible new ways of breathing, eating, communicating and sensing the environment. The head and face develop through interactions of all three germ layers, the endoderm, ectoderm and mesoderm, as well as the so-called fourth germ layer, the cranial neural crest. Over a century of experimental embryology and genetics have revealed an incredible diversity of cell types derived from each germ layer, signaling pathways and genes that coordinate craniofacial development, and how changes to these underlie human disease and vertebrate evolution. Yet for many diseases and congenital anomalies, we have an incomplete picture of the causative genomic changes, in particular how alterations to the non-coding genome might affect craniofacial gene expression. Emerging genomics and single-cell technologies provide an opportunity to obtain a more holistic view of the genes and gene regulatory elements orchestrating craniofacial development across vertebrates. These single-cell studies generate novel hypotheses that can be experimentally validated in vivo. In this Review, we highlight recent advances in single-cell studies of diverse craniofacial structures, as well as potential pitfalls and the need for extensive in vivo validation. We discuss how these studies inform the developmental sources and regulation of head structures, bringing new insights into the etiology of structural birth anomalies that affect the vertebrate head.

Single-cell genomics meets human genetics

Single-cell genomic technologies are revealing the cellular composition, identities and states in tissues at unprecedented resolution. They have now scaled to the point that it is possible to query samples at the population level, across thousands of individuals. Combining single-cell information with genotype data at this scale provides opportunities to link genetic variation to the cellular processes underpinning key aspects of human biology and disease. This strategy has potential implications for disease diagnosis, risk prediction and development of therapeutic solutions. But, effectively integrating large-scale single-cell genomic data, genetic variation and additional phenotypic data will require advances in data generation and analysis methods. As single-cell genetics begins to emerge as a field in its own right, we review its current state and the challenges and opportunities ahead.

Traumatic brain injury induces an adaptive immune response in the meningeal transcriptome that is amplified by aging

Traumatic Brain Injury (TBI) is a major cause of disability and mortality, particularly among the elderly, yet our mechanistic understanding of how age renders the post-traumatic brain vulnerable to poor clinical outcomes and susceptible to neurological disease remains poorly understood. It is well established that dysregulated and sustained immune responses contribute to negative outcomes after TBI, however our understanding of the interactions between central and peripheral immune reservoirs is still unclear. The meninges serve as the interface between the brain and the immune system, facilitating important bi-directional roles in healthy and disease settings. It has been previously shown that disruption of this system exacerbates inflammation in age related neurodegenerative disorders such as Alzheimer's disease, however we have an incomplete understanding of how the meningeal compartment influences immune responses after TBI. Here, we examine the meningeal tissue and its response to brain injury in young (3-months) and aged (18-months) mice. Utilizing a bioinformatic approach, high-throughput RNA sequencing demonstrates alterations in the meningeal transcriptome at sub-acute (7-days) and chronic (1 month) timepoints after injury. We find that age alone chronically exacerbates immunoglobulin production and B cell responses. After TBI, adaptive immune response genes are up-regulated in a temporal manner, with genes involved in T cell responses elevated sub-acutely, followed by increases in B cell related genes at chronic time points after injury. Pro-inflammatory cytokines are also implicated as contributing to the immune response in the meninges, with ingenuity pathway analysis identifying interferons as master regulators in aged mice compared to young mice following TBI. Collectively these data demonstrate the temporal series of meningeal specific signatures, providing insights into how age leads to worse neuroinflammatory outcomes in TBI.