Publications

2025

1. 

   A computational framework for mapping isoform landscape and regulatory mechanisms from spatial transcriptomics data

   Jiayu Su^†, Yiming Qu, Megan Schertzer, and 10 more authors

   Nature Biotechnology (accepted in principle), 2025

   Abs DOI Code Website

   Transcript diversity including splicing and alternative 3’end usage is crucial for cellular identity and adaptation, yet its spatial coordination remains poorly understood. Here, we present SPLISOSM (SpatiaL ISOform Statistical Modeling), a computational framework for detecting isoform-resolution patterns from spatial transcriptomics data. SPLISOSM leverages multivariate testing to account for spot- and isoform-level dependencies, demonstrating robust and theoretically grounded performance on sparse data. In the mouse brain, we identify over 1,000 spatially variable transcript diversity events, primarily in synaptic signaling pathways linked to neuropsychiatric disorders, and uncover both known and novel regulatory relationships with region-specific RNA binding proteins. We further show that these patterns are evolutionarily conserved between mouse and human prefrontal cortex. Analysis of human glioblastoma highlights pervasive transcript diversity in antigen presentation and adhesion genes associated with specific microenvironmental conditions. Together, we present a comprehensive spatial splicing analysis in the brain under normal and neoplastic conditions.

2. A foundation model of transcription across human cell types

   Xi Fu, Shentong Mo, Alejandro Buendia, and 13 more authors

   Nature, 2025

   DOI
3. A niche driven mechanism determines response and a mutation-independent therapeutic approach for myeloid malignancies

   Ioanna Mosialou, Abdullah M Ali, Rossella Labella, and 31 more authors

   Cancer Cell, 2025

   DOI

2024

1. 

   Disentangling Interpretable Factors with Supervised Independent Subspace Principal Component Analysis

   Jiayu Su^†, David A. Knowles^†, and Raul Rabadan^†

   In Advances in Neural Information Processing Systems, 2024

   Abs arXiv Code Website

   The success of machine learning models relies heavily on effectively representing high-dimensional data. However, ensuring data representations capture human-understandable concepts remains difficult, often requiring the incorporation of prior knowledge and decomposition of data into multiple subspaces. Traditional linear methods fall short in modeling more than one space, while more expressive deep learning approaches lack interpretability. Here, we introduce Supervised Independent Subspace Principal Component Analysis (sisPCA), a PCA extension designed for multi-subspace learning. Leveraging the Hilbert-Schmidt Independence Criterion (HSIC), sisPCA incorporates supervision and simultaneously ensures subspace disentanglement. We demonstrate sisPCA’s connections with autoencoders and regularized linear regression and showcase its ability to identify and separate hidden data structures through extensive applications, including breast cancer diagnosis from image features, learning aging-associated DNA methylation changes, and single-cell analysis of malaria infection. Our results reveal distinct functional pathways associated with malaria colonization, underscoring the essentiality of explainable representation in high-dimensional data analysis.

2023

1. 

   Smoother: a unified and modular framework for incorporating structural dependency in spatial omics data

   Jiayu Su^†, Jean-Baptiste Reynier, Xi Fu, and 8 more authors

   Genome Biology, 2023

   Abs HTML Code Website

   Spatial omics technologies can help identify spatially organized biological processes, but existing computational approaches often overlook structural dependencies in the data. Here, we introduce Smoother, a unified framework that integrates positional information into non-spatial models via modular priors and losses. In simulated and real datasets, Smoother enables accurate data imputation, cell-type deconvolution, and dimensionality reduction with remarkable efficiency. In colorectal cancer, Smoother-guided deconvolution reveals plasma cell and fibroblast subtype localizations linked to tumor microenvironment restructuring. Additionally, joint modeling of spatial and single-cell human prostate data with Smoother allows for spatial mapping of reference populations with significantly reduced ambiguity.

2. 

   A transcriptome-based single-cell biological age model and resource for tissue-specific aging measures

   Shulin Mao^*, Jiayu Su^*, Longteng Wang, and 3 more authors

   Genome Research, 2023

   Abs HTML Code Website

   Accurately measuring biological age is crucial for improving healthcare for the elderly population. However, the complexity of aging biology poses challenges in how to robustly estimate aging and interpret the biological significance of the traits used for estimation. Here we present SCALE, a statistical pipeline that quantifies biological aging in different tissues using explainable features learned from literature and single-cell transcriptomic data. Applying SCALE to the “Mouse Aging Cell Atlas” (Tabula Muris Senis) data, we identified tissue-level transcriptomic aging programs for more than 20 murine tissues and created a multitissue resource of mouse quantitative aging-associated genes. We observe that SCALE correlates well with other age indicators, such as the accumulation of somatic mutations, and can distinguish subtle differences in aging even in cells of the same chronological age. We further compared SCALE with other transcriptomic and methylation “clocks” in data from aging muscle stem cells, Alzheimer’s disease, and heterochronic parabiosis. Our results confirm that SCALE is more generalizable and reliable in assessing biological aging in aging-related diseases and rejuvenating interventions. Overall, SCALE represents a valuable advancement in our ability to measure aging accurately, robustly, and interpretably in single cells.

3. Single-cell multi-omics defines the cell-type-specific impact of splicing aberrations in human hematopoietic clonal outgrowths

   Mariela Cortés-López, Paulina Chamely, Allegra G. Hawkins, and 31 more authors

   Cell Stem Cell, 2023

   Abs DOI HTML

   RNA splicing factors are recurrently mutated in clonal blood disorders, but the impact of dysregulated splicing in hematopoiesis remains unclear. To overcome technical limitations, we integrated genotyping of transcriptomes (GoT) with long-read single-cell transcriptomics and proteogenomics for single-cell profiling of transcriptomes, surface proteins, somatic mutations, and RNA splicing (GoT-Splice). We applied GoT-Splice to hematopoietic progenitors from myelodysplastic syndrome (MDS) patients with mutations in the core splicing factor SF3B1. SF3B1mut cells were enriched in the megakaryocytic-erythroid lineage, with expansion of SF3B1mut erythroid progenitor cells. We uncovered distinct cryptic 3? splice site usage in different progenitor populations and stage-specific aberrant splicing during erythroid differentiation. Profiling SF3B1-mutated clonal hematopoiesis samples revealed that erythroid bias and cell-type-specific cryptic 3? splice site usage in SF3B1mut cells precede overt MDS. Collectively, GoT-Splice defines the cell-type-specific impact of somatic mutations on RNA splicing, from early clonal outgrowths to overt neoplasia, directly in human samples.

2021

1. Immunotherapy for breast cancer using EpCAM aptamer tumor-targeted gene knockdown

   Ying Zhang, Xuemei Xie, Pourya Naderi Yeganeh, and 8 more authors

   Proceedings of the National Academy of Sciences, 2021

   Abs DOI HTML

   Immunotherapy benefits some aggressive breast cancers, but many breast tumors do not respond to checkpoint blockade. Novel strategies to increase breast cancer immunogenicity are needed to improve immunotherapy. Here, we used epithelial cell adhesion molecule (EpCAM) aptamer-linked small-interfering RNA chimeras (AsiC) to selectively knock down genes in mouse breast cancers to induce tumor neoantigens or overcome immune evasion. Individual gene knockdown markedly delayed tumor growth and enhanced antitumor immunity. Cd47 and Parp1 AsiCs outperformed anti-CD47 antibody and the PARP1 inhibitor Olaparib, respectively. Combining EpCAM-AsiCs targeting multiple pathways worked better than single agents and enhanced tumor inhibition by a checkpoint inhibitor. EpCAM-AsiCs have the potential to boost immunity to tumors that are poorly responsive to checkpoint blockade. New strategies for cancer immunotherapy are needed since most solid tumors do not respond to current approaches. Here we used epithelial cell adhesion molecule EpCAM (a tumor-associated antigen highly expressed on common epithelial cancers and their tumor-initiating cells) aptamer-linked small-interfering RNA chimeras (AsiCs) to knock down genes selectively in EpCAM+ tumors with the goal of making cancers more visible to the immune system. Knockdown of genes that function in multiple steps of cancer immunity was evaluated in aggressive triple-negative and HER2+ orthotopic, metastatic, and genetically engineered mouse breast cancer models. Gene targets were chosen whose knockdown was predicted to promote tumor neoantigen expression (Upf2, Parp1, Apex1), phagocytosis, and antigen presentation (Cd47), reduce checkpoint inhibition (Cd274), or cause tumor cell death (Mcl1). Four of the six AsiC (Upf2, Parp1, Cd47, and Mcl1) potently inhibited tumor growth and boosted tumor-infiltrating immune cell functions. AsiC mixtures were more effective than individual AsiC and could synergize with anti–PD-1 checkpoint inhibition.

2020

1. Single-cell transcriptome profiling reveals neutrophil heterogeneity in homeostasis and infection

   Xuemei Xie^*, Qiang Shi^*, Peng Wu, and 14 more authors

   Nature Immunology, 2020

   Abs DOI HTML Code

   The full neutrophil heterogeneity and differentiation landscape remains incompletely characterized. Here, we profiled >25,000 differentiating and mature mouse neutrophils using single-cell RNA sequencing to provide a comprehensive transcriptional landscape of neutrophil maturation, function and fate decision in their steady state and during bacterial infection. Eight neutrophil populations were defined by distinct molecular signatures. The three mature peripheral blood neutrophil subsets arise from distinct maturing bone marrow neutrophil subsets. Driven by both known and uncharacterized transcription factors, neutrophils gradually acquire microbicidal capability as they traverse the transcriptional landscape, representing an evolved mechanism for fine-tuned regulation of an effective but balanced neutrophil response. Bacterial infection reprograms the genetic architecture of neutrophil populations, alters dynamic transitions between subpopulations and primes neutrophils for augmented functionality without affecting overall heterogeneity. In summary, these data establish a reference model and general framework for studying neutrophil-related disease mechanisms, biomarkers and therapeutic targets at single-cell resolution.