Identifying disease-critical cell types and cellular processes by integrating single-cell RNA-sequencing and human genetics

Progress in single-cell sequencing of retinal vein occlusion or ischemic hypoxic retinopathy

Retinal vein occlusion (RVO) and ischemic hypoxic retinopathy (IHR) are leading cause of irreversible vision loss worldwide, compelled by complex microvascular dysfunction, neuroinflammation, and tissue hypoxia. Despite advances in imaging and treatment, a comprehensive understanding of cellular and molecular heterogeneity underlying these pathologies remains limited. Recently, single-cell RNA sequencing (scRNA-seq) has emerged as a transformative technology, enabling unprecedented resolution of cellular dynamics, transcriptomic landscapes, and intracellular communication within the retina. Single-cell technologies continue to evolve, they are poised to revolutionize our understanding of retinal vascular diseases, ultimately paving the way for precision diagnostics and targeted interventions. This technique has revolutionized our understanding regarding complex biological systems and enables proper analysis of cellular heterogeneity. This review highlights the recent progress for the application SCS to dissect the pathophysiology of RVO and IHR. Moreover, current study summarizes findings on altered gene expression endothelial cells, Muller glia, micro glia and photoreceptors under ischemic and hypoxic stress, shedding light on potential therapeutic targets and biomarkers. Furthermore, this study explores the integration of snRNA-seq, spatial transcriptomics, and multi-omics approaches to enhance the spatial and temporal mapping of retinal responses. Additionally, discuss the current challenges, including sample preservation, retinal cell-type annotation, and cross-species translation, while offering insights into future directions such as personalized medicine and regenerative strategies. This paper aims to provide clinicians and researchers with a comprehensive update on the rapidly expanding frontier of single-cell analysis in retinal ischemic diseases.

Integrating microbial GWAS and single-cell transcriptomics reveals associations between host cell populations and the gut microbiome

Microbial genome-wide association studies (GWAS) have uncovered numerous host genetic variants associated with gut microbiota. However, links between host genetics, the gut microbiome and specific cellular contexts remain unclear. Here we use a computational framework, scBPS (single-cell Bacteria Polygenic Score), to integrate existing microbial GWAS and single-cell RNA-sequencing profiles of 24 human organs, including the liver, pancreas, lung and intestine, to identify host tissues and cell types relevant to gut microbes. Analysing 207 microbial taxa and 254 host cell types, scBPS-inferred cellular enrichments confirmed known biology such as dominant communications between gut microbes and the digestive tissue module and liver epithelial cell compartment. scBPS also identified a robust association between Collinsella and the central-veinal hepatocyte subpopulation. We experimentally validated the causal effects of Collinsella on cholesterol metabolism in mice through single-nuclei RNA sequencing on liver tissue to identify relevant cell subpopulations. Mechanistically, oral gavage of Collinsella modulated cholesterol pathway gene expression in central-veinal hepatocytes. We further validated our approach using independent microbial GWAS data, alongside single-cell and bulk transcriptomic analyses, demonstrating its robustness and reproducibility. Together, scBPS enables a systematic mapping of the host-microbe crosstalk by linking cell populations to their interacting gut microbes.

Advances and applications in single-cell and spatial genomics

The applications of single-cell and spatial technologies in recent times have revolutionized the present understanding of cellular states and the cellular heterogeneity inherent in complex biological systems. These advancements offer unprecedented resolution in the examination of the functional genomics of individual cells and their spatial context within tissues. In this review, we have comprehensively discussed the historical development and recent progress in the field of single-cell and spatial genomics. We have reviewed the breakthroughs in single-cell multi-omics technologies, spatial genomics methods, and the computational strategies employed toward the analyses of single-cell atlas data. Furthermore, we have highlighted the advances made in constructing cellular atlases and their clinical applications, particularly in the context of disease. Finally, we have discussed the emerging trends, challenges, and opportunities in this rapidly evolving field.

Spatially resolved mapping of cells associated with human complex traits

Depicting spatial distributions of disease-relevant cells is crucial for understanding disease pathology1,2. Here we present genetically informed spatial mapping of cells for complex traits (gsMap), a method that integrates spatial transcriptomics data with summary statistics from genome-wide association studies to map cells to human complex traits, including diseases, in a spatially resolved manner. Using embryonic spatial transcriptomics datasets covering 25 organs, we benchmarked gsMap through simulation and by corroborating known trait-associated cells or regions in various organs. Applying gsMap to brain spatial transcriptomics data, we reveal that the spatial distribution of glutamatergic neurons associated with schizophrenia more closely resembles that for cognitive traits than that for mood traits such as depression. The schizophrenia-associated glutamatergic neurons were distributed near the dorsal hippocampus, with upregulated expression of calcium signalling and regulation genes, whereas depression-associated glutamatergic neurons were distributed near the deep medial prefrontal cortex, with upregulated expression of neuroplasticity and psychiatric drug target genes. Our study provides a method for spatially resolved mapping of trait-associated cells and demonstrates the gain of biological insights (such as the spatial distribution of trait-relevant cells and related signature genes) through these maps.

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.

Discovering genetically-supported drug targets for multisite chronic pain through multi-omics Mendelian randomization and single-cell RNA-sequencing

Multisite chronic pain (MCP) is a highly prevalent disorder with substantial unmet therapeutic needs.We conducted multi-omics Mendelian randomization and Bayesian colocalization to identify potential therapeutic targets for MCP. Summary-level data of gene expressions and protein abundance levels were obtained from corresponding quantitative trait loci studies, respectively. Summary-level data for MCP was leveraged from the UK Biobank. The transcriptome-wide association study (TWAS), Mendelian randomization, and Bayesian colocalization approaches were applied to investigate the potential causal effects of gene expressions and protein levels on MCP in both blood and brain tissues. Phenome-wide Mendelian randomization analysis (MR-PheWAS), single-cell sequencing data, protein-protein interaction (PPI), and reaction pathway analysis were further conducted to digging the underlying mechanisms. Our analysis identified and validated two plasma targets for MCP, namely KLC1 and LANCL1, at both gene expression levels and protein levels across multi-methodologies. Moreover, MR-PheWAS observed additional benefits associated with these two targets. Through analyses based on single-cell sequencing data, we identified critical cell types for KLC1, primarily megakaryocytes, and neurons, notably linked to the axon guidance pathway, while LANCL1 showed associations with B lymphocytes, neurons, and the electron transport pathway. In dorsal root ganglions, we identified enrichments of both LANCL1 and KLC1 in putative silent nociceptors. The effects are possibly mediated through axonal transport and the activation of NMDARs, supported by PPI and reaction pathway analysis. Our multi-dimensional analysis suggests that genetically determined KLC1 and LANCL1 are causally linked to MCP risk, holding promise as appealing drug targets for MCP.

Heterogeneity-preserving discriminative feature selection for disease-specific subtype discovery

Disease-specific subtype identification can deepen our understanding of disease progression and pave the way for personalized therapies, given the complexity of disease heterogeneity. Large-scale transcriptomic, proteomic, and imaging datasets create opportunities for discovering subtypes but also pose challenges due to their high dimensionality. To mitigate this, many feature selection methods focus on selecting features that distinguish known diseases or cell states, yet often miss features that preserve heterogeneity and reveal new subtypes. To overcome this gap, we develop Preserving Heterogeneity (PHet), a statistical methodology that employs iterative subsampling and differential analysis of interquartile range, in conjunction with Fisher's method, to identify a small set of features that enhance subtype clustering quality. Here, we show that this method can maintain sample heterogeneity while distinguishing known disease/cell states, with a tendency to outperform previous differential expression and outlier-based methods, indicating its potential to advance our understanding of disease mechanisms and cell differentiation.

Efficient count-based models improve power and robustness for large-scale single-cell eQTL mapping

Population-scale single-cell transcriptomic technologies (scRNA-seq) enable characterizing variant effects on gene regulation at the cellular level (e.g., single-cell eQTLs; sc-eQTLs). However, existing sc-eQTL mapping approaches are either not designed for analyzing sparse counts in scRNA-seq data or can become intractable in extremely large datasets. Here, we propose jaxQTL, a flexible and efficient sc-eQTL mapping framework using highly efficient count-based models given pseudobulk data. Using extensive simulations, we demonstrated that jaxQTL with a negative binomial model outperformed other models in identifying sc-eQTLs, while maintaining a calibrated type I error. We applied jaxQTL across 14 cell types of OneK1K scRNA-seq data (N=982), and identified 11-16% more eGenes compared with existing approaches, primarily driven by jaxQTL ability to identify lowly expressed eGenes. We observed that fine-mapped sc-eQTLs were further from transcription starting site (TSS) than fine-mapped eQTLs identified in all cells (bulk-eQTLs; P=1x10-4) and more enriched in cell-type-specific enhancers (P=3x10-10), suggesting that sc-eQTLs improve our ability to identify distal eQTLs that are missed in bulk tissues. Overall, the genetic effect of fine-mapped sc-eQTLs were largely shared across cell types, with cell-type-specificity increasing with distance to TSS. Lastly, we observed that sc-eQTLs explain more SNP-heritability (h2 ) than bulk-eQTLs (9.90 ± 0.88% vs. 6.10 ± 0.76% when meta-analyzed across 16 blood and immune-related traits), improving but not closing the missing link between GWAS and eQTLs. As an example, we highlight that sc-eQTLs in T cells (unlike bulk-eQTLs) can successfully nominate IL6ST as a candidate gene for rheumatoid arthritis. Overall, jaxQTL provides an efficient and powerful approach using count-based models to identify missing disease-associated eQTLs.

An Ocular Gene-Set Expression Library for Heritability Partition and Cell Line Enrichment Analyses

Purpose

Use of genome-wide association studies (GWASs) in combination with transcriptomic arrays of different tissues or cell lines can reveal relevant cellular profiles and significantly improve understanding of the mechanisms of diseases. However, due to difficulty of access, few ocular transcriptomics datasets are available. This work aimed to create and make available an expression library platform that can be used with popular and versatile tools such as the linkage disequilibrium score (LDSC) regression techniques to identify specific cell lines enriched in ocular diseases.

Methods

We used transcriptome information from six publicly available single-cell and single-nucleus RNA sequence datasets to obtain matrices of normalized gene expression and estimated enrichment of the 10% most expressed transcripts in each cell line. We tested for type 1 error using simulated GWAS datasets and then used LDSC regression analyses to study the enrichment in different eye diseases.

Results

Gene expression databases for over 197 ocular cell lines were generated. Simulations of 1000 random GWASs of varying sample sizes showed no genomic inflation. Cell line-specific analyses of GWAS results found that genes near single nucleotide polymorphisms (SNPs) associated with refractive error were significantly enriched in cones (P = 0.008), rods (P = 0.002) and peripheral retina Müller cells (P = 0.002), juxtacanalicular cribriform cells (P = 0.0017), stromal keratocytes (P = 0.0018), and one beam-cell subpopulation (P = 0.0028) in primary open-angle glaucoma, emphasizing the importance of intraocular pressure in its etiology.

Conclusions

We have built a structured ocular transcriptomics library to estimate cell line enrichment among association signals from genome-wide association studies that may be extended by incorporating other datasets. We identified cells that appear important in the genetics of common eye diseases.

Large-scale genome-wide association analyses identify novel genetic loci and mechanisms in hypertrophic cardiomyopathy

Hypertrophic cardiomyopathy (HCM) is an important cause of morbidity and mortality with both monogenic and polygenic components. Here, we report results from a large genome-wide association study and multitrait analysis including 5,900 HCM cases, 68,359 controls and 36,083 UK Biobank participants with cardiac magnetic resonance imaging. We identified 70 loci (50 novel) associated with HCM and 62 loci (20 novel) associated with relevant left ventricular traits. Among the prioritized genes in the HCM loci, we identify a novel HCM disease gene, SVIL, which encodes the actin-binding protein supervillin, showing that rare truncating SVIL variants confer a roughly tenfold increased risk of HCM. Mendelian randomization analyses support a causal role of increased left ventricular contractility in both obstructive and nonobstructive forms of HCM, suggesting common disease mechanisms and anticipating shared response to therapy. Taken together, these findings increase our understanding of the genetic basis of HCM, with potential implications for disease management.

GWAS meta-analysis of psoriasis identifies new susceptibility alleles impacting disease mechanisms and therapeutic targets

Psoriasis is a common, debilitating immune-mediated skin disease. Genetic studies have identified biological mechanisms of psoriasis risk, including those targeted by effective therapies. However, the genetic liability to psoriasis is not fully explained by variation at robustly identified risk loci. To refine the genetic map of psoriasis susceptibility we meta-analysed 18 GWAS comprising 36,466 cases and 458,078 controls and identified 109 distinct psoriasis susceptibility loci, including 46 that have not been previously reported. These include susceptibility variants at loci in which the therapeutic targets IL17RA and AHR are encoded, and deleterious coding variants supporting potential new drug targets (including in STAP2, CPVL and POU2F3). We conducted a transcriptome-wide association study to identify regulatory effects of psoriasis susceptibility variants and cross-referenced these against single cell expression profiles in psoriasis-affected skin, highlighting roles for the transcriptional regulation of haematopoietic cell development and epigenetic modulation of interferon signalling in psoriasis pathobiology.

Elucidating shared genetic association between female body mass index and preeclampsia

The prevalence of obesity is steadily rising and poses a significant challenge to women's health. Preeclampsia (PE), a leading cause of maternal and fetal mortality, is significantly linked to a high body mass index (BMI). However, the shared genetic architecture underlying these conditions remains poorly understood. In this study, we used summary-level data from large-scale genome-wide association studies of BMI (N = 434,794) and PE (Ncases = 8185; Ncontrols = 234,147) to assess the shared genetic architecture between them. Our findings revealed a significant genetic correlation between BMI and PE, with an estimated sample overlap of approximately 0.8%. We identified roughly 1100 shared genetic variants, with the most notable region of local genetic correlation located in 16q12.2. Enrichment analyses highlighted endothelial dysfunction as a key biological mechanism linking BMI and PE. Additionally, RABEP2 was identified as a novel shared risk gene. Mendelian randomization analysis demonstrated a bidirectional causal relationship between BMI and PE, with blood pressure identified as a key mediator. We identified the shared genetic foundation between BMI and PE, providing valuable insights into the comorbidity of these conditions and offering a new framework for future research into comorbidity.

Conditional similarity triplets enable covariate-informed representations of single-cell data

BACKGROUND

Single-cell technologies enable comprehensive profiling of diverse immune cell-types through the measurement of multiple genes or proteins per individual cell. In order to translate immune signatures assayed from blood or tissue into powerful diagnostics, machine learning approaches are often employed to compute immunological summaries or per-sample featurizations, which can be used as inputs to models for outcomes of interest. Current supervised learning approaches for computing per-sample representations are trained only to accurately predict a single outcome and do not take into account relevant additional clinical features or covariates that are likely to also be measured for each sample.

RESULTS

Here, we introduce a novel approach for incorporating measured covariates in optimizing model parameters to ultimately specify per-sample encodings that accurately affect both immune signatures and additional clinical information. Our introduced method CytoCoSet is a set-based encoding method for learning per-sample featurizations, which formulates a loss function with an additional triplet term penalizing samples with similar covariates from having disparate embedding results in per-sample representations.

CONCLUSIONS

Overall, incorporating clinical covariates enables the learning of encodings for each individual sample that ultimately improve prediction of clinical outcome. This integration of information disparate more robust predictions of clinical phenotypes and holds significant potential for enhancing diagnostic and treatment strategies.

Genetic and Plasma Proteomic Approaches to Identify Therapeutic Targets for Graves' Disease and Graves' Ophthalmopathy

Background

The blood proteome is a major source of biomarkers and therapeutic targets. We aimed to identify the causal proteins and potential targets for Graves' disease (GD) and Graves' ophthalmopathy (GO) via systematic genetic analyses.

Methods

Genome-wide association studies (GWASs) on the UK Biobank- Pharma Proteomics Project (UKB-PPP) collected 2923 Olink proteins from 54,219 participants. We conducted a proteome-wide Mendelian randomization (MR) study with cis-pQTLs to identify candidate proteins for GD and GO risk. Colocalization analysis and the Heidi test were used to examine whether the identified proteins and diseases shared the same variant. More proteins with potential causal associations were identified in Summary-data-based MR (SMR) analyses using trans-pQTLs. Then, downstream analyses were performed to detect protein interactions, gene function, cell type-specific expression and druggable information.

Results

This study genetically predicted levels of 62 plasma proteins were associated with GD risk. Four proteins (CD40, TINAGL1, GMPR and CXCL10) were prioritized with the evidence of sharing the same variants with GD. Specifically, some proteins had potential associations with GD with trans-pQTLs mapping in CD40. The four prioritized protein-coding genes were mainly enriched in the regulation of apoptotic and death processes. In addition, GMPR was associated with both GO and GD in a consistent direction. BTN1A1 and FCRL1 were prioritized as the causal proteins for GO onset and were not associated with GD.

Conclusion

By synthesizing proteomic and genetic data, we identified several protein biomarkers for GD, with one linked to both GD and GO and two other protein biomarkers specific to GO onset, which provides valuable insights into the etiology and potential therapeutic targets for the two diseases.

A general principle of neuronal evolution reveals a human-accelerated neuron type potentially underlying the high prevalence of autism in humans

The remarkable ability of a single genome sequence to encode a diverse collection of distinct cell types, including the thousands of cell types found in the mammalian brain, is a key characteristic of multicellular life. While it has been observed that some cell types are far more evolutionarily conserved than others, the factors driving these differences in evolutionary rate remain unknown. Here, we hypothesized that highly abundant neuronal cell types may be under greater selective constraint than rarer neuronal types, leading to variation in their rates of evolution. To test this, we leveraged recently published cross-species single-nucleus RNA-sequencing datasets from three distinct regions of the mammalian neocortex. We found a strikingly consistent relationship where more abundant neuronal subtypes show greater gene expression conservation between species, which replicated across three independent datasets covering >106 neurons from six species. Based on this principle, we discovered that the most abundant type of neocortical neurons-layer 2/3 intratelencephalic excitatory neurons-has evolved exceptionally quickly in the human lineage compared to other apes. Surprisingly, this accelerated evolution was accompanied by the dramatic down-regulation of autism-associated genes, which was likely driven by polygenic positive selection specific to the human lineage. In sum, we introduce a general principle governing neuronal evolution and suggest that the exceptionally high prevalence of autism in humans may be a direct result of natural selection for lower expression of a suite of genes that conferred a fitness benefit to our ancestors while also rendering an abundant class of neurons more sensitive to perturbation.

Single-cell omics in inflammatory bowel disease: recent insights and future clinical applications

Inflammatory bowel diseases (IBDs), which include ulcerative colitis (UC) and Crohn's disease (CD), are chronic conditions characterised by inflammation of the intestinal tract. Alterations in virtually all intestinal cell types, including immune, epithelial and stromal cells, have been described in these diseases. The study of IBD has historically relied on bulk transcriptomics, but this method averages signals across diverse cell types, limiting insights. Single-cell omic technologies overcome the intrinsic limitations of bulk analysis and reveal the complexity of multicellular tissues at a cell-by-cell resolution. Within healthy and inflamed intestinal tissues, single-cell omics, particularly single-cell RNA sequencing, have contributed to uncovering novel cell types and cell functions linked to disease activity or the development of complications. Collectively, these results help identify therapeutic targets in difficult-to-treat complications such as fibrostenosis, creeping fat accumulation, perianal fistulae or inflammation of the pouch. More recently, single-cell omics have gradually been adopted in studies to understand therapeutic responses, identify mechanisms of drug failure and potentially develop predictors with clinical utility. Although these are early days, such studies lay the groundwork for the implementation in clinical practice of new technologies in diagnostics, monitoring and prediction of disease prognosis. With this review, we aim to provide a comprehensive survey of the studies that have applied single-cell omics to the study of UC or CD, and offer our perspective on the main findings these studies contribute. Finally, we discuss the limitations and potential benefits that the integration of single-cell omics into clinical practice and drug development could offer.

Inferring gene regulatory networks from single-cell multiome data using atlas-scale external data

Existing methods for gene regulatory network (GRN) inference rely on gene expression data alone or on lower resolution bulk data. Despite the recent integration of chromatin accessibility and RNA sequencing data, learning complex mechanisms from limited independent data points still presents a daunting challenge. Here we present LINGER (Lifelong neural network for gene regulation), a machine-learning method to infer GRNs from single-cell paired gene expression and chromatin accessibility data. LINGER incorporates atlas-scale external bulk data across diverse cellular contexts and prior knowledge of transcription factor motifs as a manifold regularization. LINGER achieves a fourfold to sevenfold relative increase in accuracy over existing methods and reveals a complex regulatory landscape of genome-wide association studies, enabling enhanced interpretation of disease-associated variants and genes. Following the GRN inference from reference single-cell multiome data, LINGER enables the estimation of transcription factor activity solely from bulk or single-cell gene expression data, leveraging the abundance of available gene expression data to identify driver regulators from case-control studies.

HCNetlas: A reference database of human cell type-specific gene networks to aid disease genetic analyses

Cell type-specific actions of disease genes add a significant layer of complexity to the genetic architecture underlying diseases, obscuring our understanding of disease mechanisms. Single-cell omics have revealed the functional roles of genes at the cellular level, identifying cell types critical for disease progression. Often, a gene impact on disease through its altered network within specific cell types, rather than mere changes in expression levels. To explore the cell type-specific roles of disease genes, we developed HCNetlas (human cell network atlas), a resource cataloging cell type-specific gene networks (CGNs) for various healthy tissue cells. We also devised 3 network analysis methods to investigate cell type-specific functions of disease genes. These methods involve comparing HCNetlas CGNs with those derived from disease-affected tissue samples. These methods find that systemic lupus erythematosus genes predominantly function in myeloid cells, and Alzheimer's disease genes mainly play roles in inhibitory and excitatory neurons. Additionally, they suggest that many lung cancer-related genes may exert their roles in immune cells. These findings suggest that HCNetlas has the potential to link disease-associated genes to cell types of action, facilitating development of cell type-resolved diagnostics and therapeutic strategies for complex human diseases.

Causal modeling of gene effects from regulators to programs to traits: integration of genetic associations and Perturb-seq

Genetic association studies provide a unique tool for identifying causal links from genes to human traits and diseases. However, it is challenging to determine the biological mechanisms underlying most associations, and we lack genome-scale approaches for inferring causal mechanistic pathways from genes to cellular functions to traits. Here we propose new approaches to bridge this gap by combining quantitative estimates of gene-trait relationships from loss-of-function burden tests with gene-regulatory connections inferred from Perturb-seq experiments in relevant cell types. By combining these two forms of data, we aim to build causal graphs in which the directional associations of genes with a trait can be explained by their regulatory effects on biological programs or direct effects on the trait. As a proof-of-concept, we constructed a causal graph of the gene regulatory hierarchy that jointly controls three partially co-regulated blood traits. We propose that perturbation studies in trait-relevant cell types, coupled with gene-level effect sizes for traits, can bridge the gap between genetics and biology.

Utilizing sc-linker to integrate single-cell RNA sequencing and human genetics to identify cell types and driver genes associated with non-small cell lung cancer

BACKGROUND

Genome-wide association studies (GWAS) provide a powerful method for identifying the loci and genes that contribute to disease. However, in many cases, the specific cell types and states that confer disease risk through these genes remain unknown. Determining this relationship is crucial for identifying pathogenic processes and developing therapeutic strategies.

METHODS

In this study, we utilized the sc-linker framework developed by Jagadeesh, which is an integrated framework that combines single-cell RNA sequencing (scRNA-seq), epigenomic single nucleotide polymorphism (SNP)-to-gene mapping, and GWAS summary statistics to infer potential cell types and diseases affected by genetic variations.

RESULTS

Using normal cell type programs in the sc-linker, we identified type 2 alveolar cells in normal lung tissues that are closely associated with non-small cell lung cancer (NSCLC). Additionally, we identified cancer-associated fibroblasts (CAFs) associated with lung cancer using disease-dependent programs. By integrating extensive single-cell data from NSCLC, we discerned heterogeneity among CAFs subgroups. Finally, using MAGMA, we identified RAB31 as a driver gene in disease-related fibroblasts. Proteins from the RAB family are involved in the dynamic regulation of cell membrane compartments and are dysregulated in various tumor types, potentially altering biological properties such as the proliferation, migration, and invasion of cancer cells. We found that the Ras-related protein Rab-31 (RAB31) was significantly overexpressed in tumor-associated fibroblasts compared to that in normal fibroblasts and was closely associated with poor prognosis in patients with NSCLC.

CONCLUSIONS

By integrating scRNA-seq, epigenomic, and GWAS datasets, we found that ACT2 and CAFs have specific disease heritability in lung cancer and identified the driver gene RAB31 as a potential therapeutic target.

Machine learning-driven identification of critical gene programs and key transcription factors in migraine

BACKGROUND

Migraine is a complex neurological disorder characterized by recurrent episodes of severe headaches. Although genetic factors have been implicated, the precise molecular mechanisms, particularly gene expression patterns in migraine-associated brain regions, remain unclear. This study applies machine learning techniques to explore region-specific gene expression profiles and identify critical gene programs and transcription factors linked to migraine pathogenesis.

METHODS

We utilized single-nucleus RNA sequencing (snRNA-seq) data from 43 brain regions, along with genome-wide association study (GWAS) data, to investigate susceptibility to migraine. The cell-type-specific expression (CELLEX) algorithm was employed to calculate specific expression profiles for each region, while non-negative matrix factorization (NMF) was applied to decompose gene programs within the single-cell data from these regions. Following the annotation of brain region expression profiles and gene programs to the genome, we employed stratified linkage disequilibrium score regression (S-LDSC) to assess the associations between brain regions, gene programs, and migraine-related SNPs. Key transcription factors regulating critical gene programs were identified using a random forest model based on regulatory networks derived from the GTEx consortium.

RESULTS

Our analysis revealed significant enrichment of migraine-associated single nucleotide polymorphisms (SNPs) in the posterior nuclear complex-medial geniculate nuclei (PoN_MG) of the thalamus, highlighting this region's crucial role in migraine pathogenesis. Gene program 1, identified through NMF, was enriched in the calcium signaling pathway, a known contributor to migraine pathophysiology. Random forest analysis predicted ARID3A as the top transcription factor regulating gene program 1, suggesting its potential role in modulating calcium-related genes involved in migraine.

CONCLUSION

This study provides new insights into the molecular mechanisms underlying migraine, emphasizing the importance of the PoN_MG thalamic region, calcium signaling pathways, and key transcription factors like ARID3A. These findings offer potential avenues for developing targeted therapeutic strategies for migraine treatment.

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.

Connecting genomic results for psychiatric disorders to human brain cell types and regions reveals convergence with functional connectivity

Identifying cell types and brain regions critical for psychiatric disorders and brain traits is essential for targeted neurobiological research. By integrating genomic insights from genome-wide association studies with a comprehensive single-cell transcriptomic atlas of the adult human brain, we prioritized specific neuronal clusters significantly enriched for the SNP-heritabilities for schizophrenia, bipolar disorder, and major depressive disorder along with intelligence, education, and neuroticism. Extrapolation of cell-type results to brain regions reveals the whole-brain impact of schizophrenia genetic risk, with subregions in the hippocampus and amygdala exhibiting the most significant enrichment of SNP-heritability. Using functional MRI connectivity, we further confirmed the significance of the central and lateral amygdala, hippocampal body, and prefrontal cortex in distinguishing schizophrenia cases from controls. Our findings underscore the value of single-cell transcriptomics in understanding the polygenicity of psychiatric disorders and suggest a promising alignment of genomic, transcriptomic, and brain imaging modalities for identifying common biological targets.

Fine-mapping causal tissues and genes at disease-associated loci

Complex diseases often have distinct mechanisms spanning multiple tissues. We propose tissue-gene fine-mapping (TGFM), which infers the posterior inclusion probability (PIP) for each gene-tissue pair to mediate a disease locus by analyzing summary statistics and expression quantitative trait loci (eQTL) data; TGFM also assigns PIPs to non-mediated variants. TGFM accounts for co-regulation across genes and tissues and models uncertainty in cis-predicted expression models, enabling correct calibration. We applied TGFM to 45 UK Biobank diseases or traits using eQTL data from 38 Genotype-Tissue Expression (GTEx) tissues. TGFM identified an average of 147 PIP > 0.5 causal genetic elements per disease or trait, of which 11% were gene-tissue pairs. Causal gene-tissue pairs identified by TGFM reflected both known biology (for example, TPO-thyroid for hypothyroidism) and biologically plausible findings (for example, SLC20A2-artery aorta for diastolic blood pressure). Application of TGFM to single-cell eQTL data from nine cell types in peripheral blood mononuclear cells (PBMCs), analyzed jointly with GTEx tissues, identified 30 additional causal gene-PBMC cell type pairs.

The Human Cell Atlas from a cell census to a unified foundation model

With the convergence of notable advances in molecular and spatial profiling methods and new computational approaches taking advantage of artificial intelligence and machine learning, the construction of cell atlases is progressing from data collection to atlas integration and beyond. Here, we explore five ways in which cell atlases, including the Human Cell Atlas, are already revealing valuable biological insights, and how they are poised to provide even greater benefits in the coming years. In particular, we discuss cell atlases as censuses of cells; as 3D maps of cells in the body, across modalities and scales; as maps connecting genotype causes to phenotype effects; as 4D maps of development; and, ultimately, as foundation models of biology unifying all of these aspects and helping to transform medicine.

Multi-omics in exploring the pathophysiology of diabetic retinopathy

Diabetic retinopathy (DR) is a leading global cause of vision impairment, with its prevalence increasing alongside the rising rates of diabetes mellitus (DM). Despite the retina's complex structure, the underlying pathology of DR remains incompletely understood. Single-cell RNA sequencing (scRNA-seq) and recent advancements in multi-omics analyses have revolutionized molecular profiling, enabling high-throughput analysis and comprehensive characterization of complex biological systems. This review highlights the significant contributions of scRNA-seq, in conjunction with other multi-omics technologies, to DR research. Integrated scRNA-seq and transcriptomic analyses have revealed novel insights into DR pathogenesis, including alternative transcription start site events, fluctuations in cell populations, altered gene expression profiles, and critical signaling pathways within retinal cells. Furthermore, by integrating scRNA-seq with genetic association studies and multi-omics analyses, researchers have identified novel biomarkers, susceptibility genes, and potential therapeutic targets for DR, emphasizing the importance of specific retinal cell types in disease progression. The integration of scRNA-seq with metabolomics has also been instrumental in identifying specific metabolites and dysregulated pathways associated with DR. It is highly conceivable that the continued synergy between scRNA-seq and other multi-omics approaches will accelerate the discovery of underlying mechanisms and the development of novel therapeutic interventions for DR.

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.

Understanding genetic variants in context

Over the last three decades, human genetics has gone from dissecting high-penetrance Mendelian diseases to discovering the vast and complex genetic etiology of common human diseases. In tackling this complexity, scientists have discovered the importance of numerous genetic processes - most notably functional regulatory elements - in the development and progression of these diseases. Simultaneously, scientists have increasingly used multiplex assays of variant effect to systematically phenotype the cellular consequences of millions of genetic variants. In this article, we argue that the context of genetic variants - at all scales, from other genetic variants and gene regulation to cell biology to organismal environment - are critical components of how we can employ genomics to interpret these variants, and ultimately treat these diseases. We describe approaches to extend existing experimental assays and computational approaches to examine and quantify the importance of this context, including through causal analytic approaches. Having a unified understanding of the molecular, physiological, and environmental processes governing the interpretation of genetic variants is sorely needed for the field, and this perspective argues for feasible approaches by which the combined interpretation of cellular, animal, and epidemiological data can yield that knowledge.

Genome-wide association analysis provides insights into the molecular etiology of dilated cardiomyopathy

Dilated cardiomyopathy (DCM) is a leading cause of heart failure and cardiac transplantation. We report a genome-wide association study and multi-trait analysis of DCM (14,256 cases) and three left ventricular traits (36,203 UK Biobank participants). We identified 80 genomic risk loci and prioritized 62 putative effector genes, including several with rare variant DCM associations (MAP3K7, NEDD4L and SSPN). Using single-nucleus transcriptomics, we identify cellular states, biological pathways, and intracellular communications that drive pathogenesis. We demonstrate that polygenic scores predict DCM in the general population and modify penetrance in carriers of rare DCM variants. Our findings may inform the design of genetic testing strategies that incorporate polygenic background. They also provide insights into the molecular etiology of DCM that may facilitate the development of targeted therapeutics.

Mining single-cell data for cell type-disease associations

A robust understanding of the cellular mechanisms underlying diseases sets the foundation for the effective design of drugs and other interventions. The wealth of existing single-cell atlases offers the opportunity to uncover high-resolution information on expression patterns across various cell types and time points. To better understand the associations between cell types and diseases, we leveraged previously developed tools to construct a standardized analysis pipeline and systematically explored associations across four single-cell datasets, spanning a range of tissue types, cell types and developmental time periods. We utilized a set of existing tools to identify co-expression modules and temporal patterns per cell type and then investigated these modules for known disease and phenotype enrichments. Our pipeline reveals known and novel putative cell type-disease associations across all investigated datasets. In addition, we found that automatically discovered gene co-expression modules and temporal clusters are enriched for drug targets, suggesting that our analysis could be used to identify novel therapeutic targets.

Spatial multiomic landscape of the human placenta at molecular resolution

Successful pregnancy relies directly on the placenta's complex, dynamic, gene-regulatory networks. Disruption of this vast collection of intercellular and intracellular programs leads to pregnancy complications and developmental defects. In the present study, we generated a comprehensive, spatially resolved, multimodal cell census elucidating the molecular architecture of the first trimester human placenta. We utilized paired single-nucleus (sn)ATAC (assay for transposase accessible chromatin) sequencing and RNA sequencing (RNA-seq), spatial snATAC-seq and RNA-seq, and in situ sequencing and hybridization mapping of transcriptomes at molecular resolution to spatially reconstruct the joint epigenomic and transcriptomic regulatory landscape. Paired analyses unraveled intricate tumor-like gene expression and transcription factor motif programs potentially sustaining the placenta in a hostile uterine environment; further investigation of gene-linked cis-regulatory elements revealed heightened regulatory complexity that may govern trophoblast differentiation and placental disease risk. Complementary spatial mapping techniques decoded these programs within the placental villous core and extravillous trophoblast cell column architecture while simultaneously revealing niche-establishing transcriptional elements and cell-cell communication. Finally, we computationally imputed genome-wide, multiomic single-cell profiles and spatially characterized the placental chromatin accessibility landscape. This spatially resolved, single-cell multiomic framework of the first trimester human placenta serves as a blueprint for future studies on early placental development and pregnancy.

Uncovering cell-type-specific immunomodulatory variants and molecular phenotypes in COVID-19 using structurally resolved protein networks

Immunomodulatory variants that lead to the loss or gain of specific protein interactions often manifest only as organismal phenotypes in infectious disease. Here, we propose a network-based approach to integrate genetic variation with a structurally resolved human protein interactome network to prioritize immunomodulatory variants in COVID-19. We find that, in addition to variants that pass genome-wide significance thresholds, variants at the interface of specific protein-protein interactions, even though they do not meet genome-wide thresholds, are equally immunomodulatory. The integration of these variants with single-cell epigenomic and transcriptomic data prioritizes myeloid and T cell subsets as the most affected by these variants across both the peripheral blood and the lung compartments. Of particular interest is a common coding variant that disrupts the OAS1-PRMT6 interaction and affects downstream interferon signaling. Critically, our framework is generalizable across infectious disease contexts and can be used to implicate immunomodulatory variants that do not meet genome-wide significance thresholds.

M2 Macrophages are Major Mediators of Germline Risk of Endometriosis and Explain Pleiotropy with Comorbid Traits

Endometriosis is a common gynecologic condition that causes chronic life-altering symptoms including pain, infertility, and elevated cancer risk. There is an urgent need for new non-hormonal targeted therapeutics to treat endometriosis, but until very recently, the cellular and molecular signatures of endometriotic lesions were undefined, severely hindering the development of clinical advances. Integrating inherited risk data from analyses of >450,000 individuals with ∼350,000 single cell transcriptomes from 21 patients, we uncover M2-macrophages as candidate drivers of disease susceptibility, and nominate IL1 signaling as a central hub impacted by germline genetic variation associated with endometriosis. Extensive functional follow-up confirmed these associations and revealed a pleiotropic role for this pathway in endometriosis. Population-scale expression quantitative trail locus analysis demonstrated that genetic variation controlling IL1A expression is also associated with endometriosis risk variants. Manipulation of IL1 signaling in state-of-the-art in vitro decidualized assembloids impacted epithelial differentiation, and in an in vivo endometriosis model, treatment with anakinra (an interleukin-1 receptor antagonist) resulted in a significant, dose-dependent reduction in both spontaneous pain and evoked pain. Together these studies highlight non-diagnostic cell types as central to endometriosis susceptibility and support IL1 signaling as an important actionable pathway for this disease.

scPAS: single-cell phenotype-associated subpopulation identifier

Despite significant advancements in single-cell sequencing analysis for characterizing tissue sample heterogeneity, identifying the associations between cell subpopulations and disease phenotypes remains a challenging task. Here, we introduce scPAS, a new bioinformatics tool designed to integrate bulk data to identify phenotype-associated cell subpopulations within single-cell data. scPAS employs a network-regularized sparse regression model to quantify the association between each cell in single-cell data and a phenotype. Additionally, it estimates the significance of these associations through a permutation test, thereby identifying phenotype-associated cell subpopulations. Utilizing simulated data and various single-cell datasets from breast carcinoma, ovarian cancer, and atherosclerosis, as well as spatial transcriptomics data from multiple cancers, we demonstrated the accuracy, flexibility, and broad applicability of scPAS. Evaluations on large datasets revealed that scPAS exhibits superior operational efficiency compared to other methods. The open-source scPAS R package is available at GitHub website: https://github.com/aiminXie/scPAS.

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.

scPerb: Predict single-cell perturbation via style transfer-based variational autoencoder

INTRODUCTION

Traditional methods for obtaining cellular responses after perturbation are usually labor-intensive and costly, especially when working with multiple different experimental conditions. Therefore, accurate prediction of cellular responses to perturbations is of great importance in computational biology. Existing methodologies, such as graph-based approaches, vector arithmetic, and neural networks, either mix perturbation-related variances with cell-type-specific patterns or implicitly distinguish them within black-box models.

OBJECTIVES

This study aims to introduce and demonstrate a novel framework, scPerb, which explicitly extracts perturbation-related variances and transfers them from unperturbed to perturbed cells to accurately predict the effect of perturbation in single-cell level.

METHODS

scPerb utilizes a style transfer strategy by incorporating a style encoder into the architecture of a variational autoencoder. The style encoder captures the differences in latent representations between unperturbed and perturbed cells, enabling accurate prediction of post-perturbation gene expression data.

RESULTS

Comprehensive comparisons with existing methods demonstrate that scPerb delivers improved performance and higher accuracy in predicting cellular responses to perturbations. Notably, scPerb outperforms other methods across multiple datasets, achieving superior R2 values of 0.98, 0.98, and 0.96 on three benchmarking datasets.

CONCLUSION

scPerb offers a significant advancement in predicting cellular responses by effectively separating and transferring perturbation-related variances. This framework not only enhances prediction accuracy but also provides a robust tool for computational biology, addressing the limitations of current methodologies.

Interpretable GWAS by linking clinical phenotypes to quantifiable immune repertoire components

Bridging the gap between genotype and phenotype in GWAS studies is challenging. A multitude of genetic variants have been associated with immune-related diseases, including cancer, yet the interpretability of most variants remains low. Here, we investigate the quantitative components in the T cell receptor (TCR) repertoire, the frequency of clusters of TCR sequences predicted to have common antigen specificity, to interpret the genetic associations of diverse human diseases. We first developed a statistical model to predict the TCR components using variants in the TRB and HLA loci. Applying this model to over 300,000 individuals in the UK Biobank data, we identified 2309 associations between TCR abundances and various immune diseases. TCR clusters predicted to be pathogenic for autoimmune diseases were significantly enriched for predicted autoantigen-specificity. Moreover, four TCR clusters were associated with better outcomes in distinct cancers, where conventional GWAS cannot identify any significant locus. Collectively, our results highlight the integral role of adaptive immune responses in explaining the associations between genotype and phenotype.

Genes with differential expression across ancestries are enriched in ancestry-specific disease effects likely due to gene-by-environment interactions

Multi-ancestry genome-wide association studies (GWASs) have highlighted the existence of variants with ancestry-specific effect sizes. Understanding where and why these ancestry-specific effects occur is fundamental to understanding the genetic basis of human diseases and complex traits. Here, we characterized genes differentially expressed across ancestries (ancDE genes) at the cell-type level by leveraging single-cell RNA-sequencing data in peripheral blood mononuclear cells for 21 individuals with East Asian (EAS) ancestry and 23 individuals with European (EUR) ancestry (172,385 cells); then, we tested whether variants surrounding those genes were enriched in disease variants with ancestry-specific effect sizes by leveraging ancestry-matched GWASs of 31 diseases and complex traits (average n ∼ 90,000 and ∼ 267,000 in EAS and EUR, respectively). We observed that ancDE genes tended to be cell-type specific and enriched in genes interacting with the environment and in variants with ancestry-specific disease effect sizes, which suggests cell-type-specific, gene-by-environment interactions shared between regulatory and disease architectures. Finally, we illustrated how different environments might have led to ancestry-specific myeloid cell leukemia 1 (MCL1) expression in B cells and ancestry-specific allele effect sizes in lymphocyte count GWASs for variants surrounding MCL1. Our results imply that large single-cell and GWAS datasets from diverse ancestries are required to improve our understanding of human diseases.

Fine mapping of candidate effector genes for heart rate

An elevated resting heart rate (RHR) is associated with increased cardiovascular mortality. Genome-wide association studies (GWAS) have identified > 350 loci. Uniquely, in this study we applied genetic fine-mapping leveraging tissue specific chromatin segmentation and colocalization analyses to identify causal variants and candidate effector genes for RHR. We used RHR GWAS summary statistics from 388,237 individuals of European ancestry from UK Biobank and performed fine mapping using publicly available genomic annotation datasets. High-confidence causal variants (accounting for > 75% posterior probability) were identified, and we collated candidate effector genes using a multi-omics approach that combined evidence from colocalisation with molecular quantitative trait loci (QTLs), and long-range chromatin interaction analyses. Finally, we performed druggability analyses to investigate drug repurposing opportunities. The fine mapping pipeline indicated 442 distinct RHR signals. For 90 signals, a single variant was identified as a high-confidence causal variant, of which 22 were annotated as missense. In trait-relevant tissues, 39 signals colocalised with cis-expression QTLs (eQTLs), 3 with cis-protein QTLs (pQTLs), and 75 had promoter interactions via Hi-C. In total, 262 candidate genes were highlighted (79% had promoter interactions, 15% had a colocalised eQTL, 8% had a missense variant and 1% had a colocalised pQTL), and, for the first time, enrichment in nervous system pathways. Druggability analyses highlighted ACHE, CALCRL, MYT1 and TDP1 as potential targets. Our genetic fine-mapping pipeline prioritised 262 candidate genes for RHR that warrant further investigation in functional studies, and we provide potential therapeutic targets to reduce RHR and cardiovascular mortality.

Oral delivery of sodium alginate/chitosan bilayer microgels loaded with Lactobacillus rhamnosus GG for targeted therapy of ulcerative colitis

Probiotics regulate intestinal flora balance and enhance the intestinal barrier, which is useful in preventing and treating colitis. However, they have strict storage requirements. In addition, they degrade in a strongly acidic environment, resulting in a significant decrease in their activity when used as microbial agents. Lactobacillus rhamnosus GG (LGG) was loaded into acid-resistant and colon-targeting double-layer microgels. The inner layer consists of guar gum (GG) and low methoxyl pectin (LMP), which can achieve retention and degradation in the colon. To achieve colon localization, the outer layer was composed of chitosan (CS) and sodium alginate (SA). The formulation demonstrated favorable bio-responses across various pH conditions in vitro and sustained release of LGG in the colon lesions. Bare LGG survival decreased by 52.2 % in simulated gastric juice (pH 1.2) for 2 h, whereas that of encapsulated LGG decreased by 18.5 %. In the DSS-induced inflammatory model, LGG-loaded microgel significantly alleviated UC symptoms in mice and reduced inflammatory factor levels in the colon. Encapsulation of LGG improved its stability in acidic conditions, thus increasing its content at the colon lesions and reducing pathogenic bacteria. These findings provide an experimental basis and a technical reference for developing and applying probiotic microgel preparations.

An efficient, not-only-linear correlation coefficient based on clustering

Identifying meaningful patterns in data is crucial for understanding complex biological processes, particularly in transcriptomics, where genes with correlated expression often share functions or contribute to disease mechanisms. Traditional correlation coefficients, which primarily capture linear relationships, may overlook important nonlinear patterns. We introduce the clustermatch correlation coefficient (CCC), a not-only-linear coefficient that utilizes clustering to efficiently detect both linear and nonlinear associations. CCC outperforms standard methods by revealing biologically meaningful patterns that linear-only coefficients miss and is faster than state-of-the-art coefficients such as the maximal information coefficient. When applied to human gene expression data from genotype-tissue expression (GTEx), CCC identified robust linear relationships and nonlinear patterns, such as sex-specific differences, that are undetectable by standard methods. Highly ranked gene pairs were enriched for interactions in integrated networks built from protein-protein interactions, transcription factor regulation, and chemical and genetic perturbations, suggesting that CCC can detect functional relationships missed by linear-only approaches. CCC is a highly efficient, next-generation, not-only-linear correlation coefficient for genome-scale data. A record of this paper's transparent peer review process is included in the supplemental information.

Identifying genetic variants associated with chromatin looping and genome function

Here we present a comprehensive HiChIP dataset on naïve CD4 T cells (nCD4) from 30 donors and identify QTLs that associate with genotype-dependent and/or allele-specific variation of HiChIP contacts defining loops between active regulatory regions (iQTLs). We observe a substantial overlap between iQTLs and previously defined eQTLs and histone QTLs, and an enrichment for fine-mapped QTLs and GWAS variants. Furthermore, we describe a distinct subset of nCD4 iQTLs, for which the significant variation of chromatin contacts in nCD4 are translated into significant eQTL trends in CD4 T cell memory subsets. Finally, we define connectivity-QTLs as iQTLs that are significantly associated with concordant genotype-dependent changes in chromatin contacts over a broad genomic region (e.g., GWAS SNP in the RNASET2 locus). Our results demonstrate the importance of chromatin contacts as a complementary modality for QTL mapping and their power in identifying previously uncharacterized QTLs linked to cell-specific gene expression and connectivity.

Schizophrenia genomics: genetic complexity and functional insights

Determining the causes of schizophrenia has been a notoriously intractable problem, resistant to a multitude of investigative approaches over centuries. In recent decades, genomic studies have delivered hundreds of robust findings that implicate nearly 300 common genetic variants (via genome-wide association studies) and more than 20 rare variants (via whole-exome sequencing and copy number variant studies) as risk factors for schizophrenia. In parallel, functional genomic and neurobiological studies have provided exceptionally detailed information about the cellular composition of the brain and its interconnections in neurotypical individuals and, increasingly, in those with schizophrenia. Taken together, these results suggest unexpected complexity in the mechanisms that drive schizophrenia, pointing to the involvement of ensembles of genes (polygenicity) rather than single-gene causation. In this Review, we describe what we now know about the genetics of schizophrenia and consider the neurobiological implications of this information.

Recent developments and new directions in the use of natural products for the treatment of inflammatory bowel disease

BACKGROUND

Inflammatory bowel disease (IBD) represents a significant global health challenge, and there is an urgent need to explore novel therapeutic interventions. Natural products have demonstrated highly promising effectiveness in the treatment of IBD.

PURPOSE

This study systematically reviews the latest research advancements in leveraging natural products for IBD treatment.

METHODS

This manuscript strictly adheres to the PRISMA guidelines. Relevant literature on the effects of natural products on IBD was retrieved from the PubMed, Web of Science and Cochrane Library databases using the search terms "natural product," "inflammatory bowel disease," "colitis," "metagenomics", "target identification", "drug delivery systems", "polyphenols," "alkaloids," "terpenoids," and so on. The retrieved data were then systematically summarized and reviewed.

RESULTS

This review assessed the different effects of various natural products, such as polyphenols, alkaloids, terpenoids, quinones, and others, in the treatment of IBD. While these natural products offer promising avenues for IBD management, they also face challenges in terms of clinical translation and drug discovery. The advent of metagenomics, single-cell sequencing, target identification techniques, drug delivery systems, and other cutting-edge technologies heralds a new era in overcoming these challenges.

CONCLUSION

This paper provides an overview of current research progress in utilizing natural products for the treatment of IBD, exploring how contemporary technological innovations can aid in discovering and harnessing bioactive natural products for the treatment of IBD.

scParser: sparse representation learning for scalable single-cell RNA sequencing data analysis

The rapid rise in the availability and scale of scRNA-seq data needs scalable methods for integrative analysis. Though many methods for data integration have been developed, few focus on understanding the heterogeneous effects of biological conditions across different cell populations in integrative analysis. Our proposed scalable approach, scParser, models the heterogeneous effects from biological conditions, which unveils the key mechanisms by which gene expression contributes to phenotypes. Notably, the extended scParser pinpoints biological processes in cell subpopulations that contribute to disease pathogenesis. scParser achieves favorable performance in cell clustering compared to state-of-the-art methods and has a broad and diverse applicability.

Enhanced In Situ Spatial Proteomics by Effective Combination of MALDI Imaging and LC-MS/MS

Highly specialized cells are fundamental for the proper functioning of complex organs. Variations in cell-type-specific gene expression and protein composition have been linked to a variety of diseases. Investigation of the distinctive molecular makeup of these cells within tissues is therefore critical in biomedical research. Although several technologies have emerged as valuable tools to address this cellular heterogeneity, most workflows lack sufficient in situ resolution and are associated with high costs and extremely long analysis times. Here, we present a combination of experimental and computational approaches that allows a more comprehensive investigation of molecular heterogeneity within tissues than by either shotgun LC-MS/MS or MALDI imaging alone. We applied our pipeline to the mouse brain, which contains a wide variety of cell types that not only perform unique functions but also exhibit varying sensitivities to insults. We explored the distinct neuronal populations within the hippocampus, a brain region crucial for learning and memory that is involved in various neurological disorders. As an example, we identified the groups of proteins distinguishing the neuronal populations of the dentate gyrus (DG) and the cornu ammonis (CA) in the same brain section. Most of the annotated proteins matched the regional enrichment of their transcripts, thereby validating the method. As the method is highly reproducible, the identification of individual masses through the combination of MALDI-IMS and LC-MS/MS methods can be used for the much faster and more precise interpretation of MALDI-IMS measurements only. This greatly speeds up spatial proteomic analyses and allows the detection of local protein variations within the same population of cells. The method's general applicability has the potential to be used to investigate different biological conditions and tissues and a much higher throughput than other techniques making it a promising approach for clinical routine applications.

Scalable genetic screening for regulatory circuits using compressed Perturb-seq

Pooled CRISPR screens with single-cell RNA sequencing readout (Perturb-seq) have emerged as a key technique in functional genomics, but they are limited in scale by cost and combinatorial complexity. In this study, we modified the design of Perturb-seq by incorporating algorithms applied to random, low-dimensional observations. Compressed Perturb-seq measures multiple random perturbations per cell or multiple cells per droplet and computationally decompresses these measurements by leveraging the sparse structure of regulatory circuits. Applied to 598 genes in the immune response to bacterial lipopolysaccharide, compressed Perturb-seq achieves the same accuracy as conventional Perturb-seq with an order of magnitude cost reduction and greater power to learn genetic interactions. We identified known and novel regulators of immune responses and uncovered evolutionarily constrained genes with downstream targets enriched for immune disease heritability, including many missed by existing genome-wide association studies. Our framework enables new scales of interrogation for a foundational method in functional genomics.

Genetic variation across and within individuals

Germline variation and somatic mutation are intricately connected and together shape human traits and disease risks. Germline variants are present from conception, but they vary between individuals and accumulate over generations. By contrast, somatic mutations accumulate throughout life in a mosaic manner within an individual due to intrinsic and extrinsic sources of mutations and selection pressures acting on cells. Recent advancements, such as improved detection methods and increased resources for association studies, have drastically expanded our ability to investigate germline and somatic genetic variation and compare underlying mutational processes. A better understanding of the similarities and differences in the types, rates and patterns of germline and somatic variants, as well as their interplay, will help elucidate the mechanisms underlying their distinct yet interlinked roles in human health and biology.

RNA Sequencing in Disease Diagnosis

RNA sequencing (RNA-seq) enables the accurate measurement of multiple transcriptomic phenotypes for modeling the impacts of disease variants. Advances in technologies, experimental protocols, and analysis strategies are rapidly expanding the application of RNA-seq to identify disease biomarkers, tissue- and cell-type-specific impacts, and the spatial localization of disease-associated mechanisms. Ongoing international efforts to construct biobank-scale transcriptomic repositories with matched genomic data across diverse population groups are further increasing the utility of RNA-seq approaches by providing large-scale normative reference resources. The availability of these resources, combined with improved computational analysis pipelines, has enabled the detection of aberrant transcriptomic phenotypes underlying rare diseases. Further expansion of these resources, across both somatic and developmental tissues, is expected to soon provide unprecedented insights to resolve disease origin, mechanism of action, and causal gene contributions, suggesting the continued high utility of RNA-seq in disease diagnosis.

Non-coding variants impact cis-regulatory coordination in a cell type-specific manner

BACKGROUND

Interactions among cis-regulatory elements (CREs) play a crucial role in gene regulation. Various approaches have been developed to map these interactions genome-wide, including those relying on interindividual epigenomic variation to identify groups of covariable regulatory elements, referred to as chromatin modules (CMs). While CM mapping allows to investigate the relationship between chromatin modularity and gene expression, the computational principles used for CM identification vary in their application and outcomes.

RESULTS

We comprehensively evaluate and streamline existing CM mapping tools and present guidelines for optimal utilization of epigenome data from a diverse population of individuals to assess regulatory coordination across the human genome. We showcase the effectiveness of our recommended practices by analyzing distinct cell types and demonstrate cell type specificity of CRE interactions in CMs and their relevance for gene expression. Integration of genotype information revealed that many non-coding disease-associated variants affect the activity of CMs in a cell type-specific manner by affecting the binding of cell type-specific transcription factors. We provide example cases that illustrate in detail how CMs can be used to deconstruct GWAS loci, assess variable expression of cell surface receptors in immune cells, and reveal how genetic variation can impact the expression of prognostic markers in chronic lymphocytic leukemia.

CONCLUSIONS

Our study presents an optimal strategy for CM mapping and reveals how CMs capture the coordination of CREs and its impact on gene expression. Non-coding genetic variants can disrupt this coordination, and we highlight how this may lead to disease predisposition in a cell type-specific manner.