Gefei Wang

Cell-cell communication is fundamental to tissue organization, homeostasis, and disease progression. Recent advances in spatial transcriptomics provide unprecedented opportunities to systematically characterize ligand-receptor interactions directly within intact tissues. However, robust inference of spatial ligand-receptor interactions remains challenging because intrinsic features of spatial transcriptomics data, including spatial autocorrelation, variation in total molecular counts, and measurement errors, can induce spurious spatial co-expression and lead to inflated false-positive results. Most existing methods do not adequately account for these confounding factors, limiting the reliability of inferred cellular communication. Here, we present CONCISE, a statistical method for spatially constrained co-expression and ligand-receptor interaction inference that jointly models spatial autocorrelation, variation in total molecular counts, measurement errors, and spatial proximity constraints. CONCISE combines efficient moment-based parameter estimation with analytical hypothesis testing, enabling fast and statistically rigorous inference without restrictive distributional assumptions. Through extensive simulations, real-data permutation experiments, and biologically motivated negative-control analyses across different spatial transcriptomics platforms, we show that most existing methods presented inflated false-positive rates, whereas CONCISE achieved well-calibrated inference, robust false-positive control, and improved detection power. Application of CONCISE to high-resolution MERFISH and CosMx datasets from intestinal inflammation and non-small cell lung cancer further highlights its biological utility in disease contexts. CONCISE uncovered inflammation-associated fibroblast-specific interactions during intestinal inflammation and delineated complex tumor-immune and tumor-stromal signaling networks within the tumor microenvironment.

Recent advances in spatial transcriptomics (ST) have generated an expanding collection of heterogeneous datasets, offering unprecedented opportunities to investigate tissue organizations and functions. However, effective interpretation and integration of data originating from diverse sources and conditions remain a major challenge. We present INSPIRE, a deep-learning method for interpretable, integrative analysis of multiple ST datasets. INSPIRE adopts an adversarial learning strategy with graph neural networks to achieve spatially informed and adaptive data integration. By incorporating non-negative matrix factorization, INSPIRE identifies interpretable spatial factors and associated gene programs that characterize tissue architecture, cell-type organization and biological processes. Across a broad range of applications, INSPIRE demonstrates superior performance in resolving fine-grained biological signals, integrating complementary strengths across technologies, capturing condition-specific variation, uncovering tumor microenvironment heterogeneity, elucidating developmental dynamics and facilitating three-dimensional tissue reconstruction. INSPIRE also scales to extremely large datasets, as demonstrated by applications to Xenium-profiled human breast cancer and Stereo-seq mouse organogenesis datasets.

Gradients provide low-dimensional representations of macroscale brain organization, yet how structural-functional gradient coupling develops and relates to behavioral and molecular features remains unclear. Here, we studied structural-functional gradient coupling across multiple metrics and spatial scales using high-resolution structural and functional connectivity from 5343 children in the Adolescent Brain Cognitive Development study and 875 adults from the Human Connectome Project. We find that gradient coupling shows developmental refinement from childhood to adulthood and distinct sex-specific patterns. Gradient coupling metrics are significantly associated with cognitive and mental health measures and enable robust out-of-sample prediction. Heritability analyses reveal that gradient coupling is strongly influenced by genetic factors. Transcriptomic analyses further demonstrate that highly heritable coupling patterns are enriched for genes expressed in deep-layer excitatory neurons. Together, our findings establish structural-functional gradient coupling as a biologically meaningful feature of brain organization that bridges macroscale connectivity, cognition, behavior, and molecular architecture.

Gradients are increasingly used to characterize the brain’s macroscale organization, offering low-dimensional representations of structural and functional connectivity. However, how structural-functional gradient coupling evolves during development and relates to behavioral and molecular features remains unclear. Here, we studied structural-functional gradient coupling across multiple metrics and spatial scales using high-resolution structural and functional connectivity from 7,025 children in the Adolescent Brain Cognitive Development study and 913 adults from the Human Connectome Project. We found that gradient coupling exhibits clear developmental refinement from childhood to adulthood and shows distinct sex-specific patterns. Gradient coupling metrics were significantly associated with a broad range of cognitive and mental health measures and enabled robust out-of-sample prediction under learning methods. Heritability analyses revealed that gradient coupling is strongly influenced by genetic factors. Transcriptomic analyses further demonstrated that highly heritable coupling patterns are enriched for genes expressed in deep-layer excitatory neurons, suggesting that gradient coupling reflects underlying cell-type-specific transcriptional architecture. Together, our findings establish structural-functional gradient coupling as a biologically meaningful feature of brain organization that bridges macroscale connectivity, cognition, behavior, and molecular architecture.

Large-scale biobanks provide comprehensive electronic health records (EHRs) that capture detailed clinical phenotypes, potentially enhancing disease risk prediction. However, traditional polygenic risk score (PRS) methods rely on simplified phenotype definitions or predefined trait sets, limiting their ability to represent the intricate phenotypic structures embedded within EHRs. To address this gap, we introduce a general framework, EEPRS, that leverages phenotype embeddings derived from EHRs to improve genetic risk prediction using only genome-wide association study (GWAS) summary statistics, enabling accurate, robust and interpretable risk prediction for a wide range of diseases. Employing embedding methods such as Word2Vec and GPT, we conducted EHR embedding-based GWAS and identified a distinct cardiovascular cluster via hierarchical clustering of genetic correlations. Across 41 clinical traits in the UK Biobank, our EEPRS framework consistently outperformed traditional single-trait PRS, particularly within this identified cluster. Validation using PRS-based phenome-wide association studies (PRS-PheWAS) further confirmed robust associations between EHR embedding-based PRS and circulatory system diseases. Furthermore, our data-adaptive method, EEPRS_optimal, employing cross-validation to select the best embedding method, leading to additional improvements in prediction. We further developed MTAG_EEPRS for multi-trait PRS, resulting in averaging 92.48% improvement in R ² for continuous traits and 24.06% in AUC for binary traits compared to single-trait PRS. Overall, EEPRS represents a robust and interpretable framework, enhancing genetic prediction accuracy through integrating EHR embeddings with single-trait and multi-trait PRS.

Recent advancements in single-cell technologies have enabled comprehensive characterization of cellular states through transcriptomic, epigenomic, and proteomic profiling at single-cell resolution. These technologies have significantly deepened our understanding of cell functions and disease mechanisms from various omics perspectives. As these technologies evolve rapidly and data resources expand, there is a growing need for computational methods that can integrate information from different modalities to facilitate joint analysis of single-cell multi-omics data. However, integrating single-cell omics datasets presents unique challenges due to varied feature correlations and technology-specific limitations. To address these challenges, we introduce scMODAL, a deep learning framework tailored for single-cell multi-omics data alignment using feature links. scMODAL integrates datasets with limited known positively correlated features, leveraging neural networks and generative adversarial networks to align cell embeddings and preserve feature topology. Our experiments demonstrate scMODAL’s effectiveness in removing unwanted variation, preserving biological information, and accurately identifying cell subpopulations across diverse datasets. scMODAL not only advances integration tasks but also supports downstream analyses such as feature imputation and feature relationship inference, offering a robust solution for advancing single-cell multi-omics research.

Understanding cellular responses to genetic perturbations is essential for understanding gene regulation and phenotype formation. While high-throughput single-cell RNA-sequencing has facilitated detailed profiling of heterogeneous transcriptional responses to perturbations at the single-cell level, there remains a pressing need for computational models that can decode the mechanisms driving these responses and accurately predict outcomes to prioritize target genes for experimental design. Here, we present scLAMBDA, a deep generative learning framework designed to model and predict single-cell transcriptional responses to genetic perturbations, including single-gene and combinatorial multi-gene perturbations. By leveraging gene embeddings derived from large language models, scLAMBDA effectively integrates prior biological knowledge and disentangles basal cell states from perturbation-specific salient representations. Through comprehensive evaluations on multiple single-cell CRISPR Perturb-seq datasets, scLAMBDA consistently outperformed state-of-the-art methods in predicting perturbation outcomes, achieving higher prediction accuracy. Notably, scLAMBDA demonstrated robust generalization to unseen target genes and perturbations, and its predictions captured both average expression changes and the heterogeneity of single-cell responses. Furthermore, its predictions enable diverse downstream analyses, including the identification of differentially expressed genes and the exploration of genetic interactions, demonstrating its utility and versatility.

Recent advancements in single-cell technologies have enabled comprehensive characterization of cellular states through transcriptomic, epigenomic, and proteomic profiling at single-cell resolution. These technologies have significantly deepened our understanding of cell functions and disease mechanisms from various omics perspectives. As these technologies evolve rapidly and data resources expand, there is a growing need for computational methods that can integrate information from different modalities to facilitate joint analysis of single-cell multi-omics data. However, integrating single-cell omics datasets presents unique challenges due to varied feature correlations and technology-specific limitations. To address these challenges, we introduce scMODAL, a deep learning framework tailored for single-cell multi-omics data alignment using feature links. scMODAL integrates datasets with limited known positively correlated features, leveraging neural networks and generative adversarial networks to align cell embeddings and preserve feature topology. Our experiments demonstrate scMODAL’s effectiveness in removing unwanted variation, preserving biological information, and accurately identifying cell subpopulations across diverse datasets. scMODAL not only advances integration tasks but also supports downstream analyses such as feature imputation and feature relationship inference, offering a robust solution for advancing single-cell multi-omics research.