Diffusion maps for high-dimensional single-cell analysis of differentiation data

A single-cell and tissue-scale analysis suite resolves Mixl1's role in heart development

Perturbation studies using gene knockouts have become a key tool for understanding the roles of regulatory genes in development. However, large-scale studies dissecting the molecular role of development master regulators in every cell type throughout the embryo are technically challenging and scarce. Here, we systematically characterize the knockout effects of the key developmental regulators T/Brachyury and Mixl1 in gastrulation and early organogenesis using single-cell profiling of chimeric mouse embryos. For the analysis of these experimental data, we present COSICC, an effective suite of statistical tools to characterize perturbation effects in complex developing cell populations. We gain insights into T's role in lateral plate mesoderm, limb development, and posterior intermediate mesoderm specification. Furthermore, we generate Mixl1 -/- embryonic chimeras and reveal the role of this key transcription factor in discrete mesoderm lineages, in particular concerning developmental dysregulation of the recently identified juxta-cardiac field.

Induced B cell receptor diversity predicts PD-1 blockade immunotherapy response

Immune checkpoint inhibitors such as anti-Programmed Death-1 antibodies (aPD-1) can be effective in treating advanced cancers. However, many patients do not respond, and the mechanisms underlying these differences remain incompletely understood. In this study, we profile a cohort of patients with locally advanced or metastatic basal cell carcinoma undergoing aPD-1 therapy using single-cell RNA sequencing, high-definition spatial transcriptomics in tumors and draining lymph nodes, and spatial immunoreceptor profiling, with long-term clinical follow-up. We find that successful responses to PD-1 inhibition are characterized by an induction of B cell receptor (BCR) clonal diversity after treatment initiation. These induced BCR clones spatially colocalize with T cell clones, facilitate their activation, and traffic alongside them between tumor and draining lymph nodes to enhance tumor clearance. Furthermore, we validated aPD-1-induced BCR diversity as a predictor of clinical response in a larger cohort of glioblastoma, melanoma, and head and neck squamous cell carcinoma patients, suggesting that this is a generalizable predictor of treatment response across many types of cancers. We find that pretreatment tumors harbor a characteristic gene expression signature that portends a higher probability of inducing BCR clonal diversity after aPD-1 therapy, and we develop a machine learning model that predicts PD-1-induced BCR clonal diversity from baseline tumor RNA sequencing. These findings underscore a dynamic role of B cell diversity during immunotherapy, highlighting its importance as a prognostic marker and a potential target for intervention in non-responders.

PLXDC1+ Tumor-Associated Pancreatic Stellate Cells Promote Desmoplastic and Immunosuppressive Niche in Pancreatic Ductal Adenocarcinoma

Pancreatic stellate cells (PSCs) contribute to pancreatic ductal adenocarcinoma (PDAC) progression and therapeutic resistance, yet their detailed functions remain unclear. This study combined RNA sequencing and assay for transposase-accessible chromatin using sequencing (ATAC-seq) on sorted PSCs from adjacent normal and PDAC tissues to investigate their transcriptional and epigenetic activation. PSCs heterogeneity and functions are characterized through bulk, single-cell, and spatial transcriptomes, as well as in situ sequencing. The clinical relevance of PSCs in immunotherapy is assessed using an in-house immune-checkpoint blockade (ICB) treatment cohort. Findings showed that stress and hypoxia signaling activated PSCs in PDAC. Three common PSCs (CPSCs) and four tumor-associated PSCs (TPSCs) are identified, each with distinct functions. CPSCs differentiated into CCL19+ TPSCs in immune-enriched regions, MYH11+ TPSCs in the stromal region, and PLXDC1+ TPSCs, which exhibited cancer-associated myofibroblasts (myCAFs) phenotype linked to poor prognosis. Notably, PLXDC1+ TPSCs, located near aggressive LRRC15+ myCAFs and SPP1+ macrophages, formed a desmoplastic and immunosuppressive niche around the tumor boundary, promoting CD8 T cell exhaustion. Single-cell transcriptomics of PDAC patients treated with ICB revealed that PLXDC1+ TPSCs correlated with poor immunotherapy efficacy. Overall, this study provides key insights into PSCs in PDAC and potential therapeutic targets.

Single-cell analysis of bidirectional reprogramming between early embryonic states identify mechanisms of differential lineage plasticities in mice

Two distinct lineages, pluripotent epiblast (EPI) and primitive (extra-embryonic) endoderm (PrE), arise from common inner cell mass (ICM) progenitors in mammalian embryos. To study how these sister identities are forged, we leveraged mouse embryonic stem (ES) cells and extra-embryonic endoderm (XEN) stem cells-in vitro counterparts of the EPI and PrE. Bidirectional reprogramming between ES and XEN coupled with single-cell RNA and ATAC-seq analyses showed distinct rates, efficiencies, and trajectories of state conversions, identifying drivers and roadblocks of reciprocal conversions. While GATA4-mediated ES-to-iXEN conversion was rapid and nearly deterministic, OCT4-, KLF4-, and SOX2-induced XEN-to-induced pluripotent stem (iPS) reprogramming progressed with diminished efficiency and kinetics. A dominant PrE transcriptional program, safeguarded by GATA4, alongside elevated chromatin accessibility and reduced DNA methylation of the EPI underscored the differential plasticities of the two states. Mapping in vitro to embryo trajectories tracked reprogramming cells in either direction along EPI and PrE in vivo states, without transitioning through the ICM.

Reconstructing signaling histories of single cells via perturbation screens and transfer learning

Manipulating the signaling environment is an effective approach to alter cellular states for broad-ranging applications, from engineering tissues to treating diseases. Such manipulation requires knowing the signaling states and histories of the cells in situ , for which high-throughput discovery methods are lacking. Here, we present an integrated experimental-computational framework that learns signaling response signatures from a high-throughput in vitro perturbation atlas and infers combinatorial signaling activities in in vivo cell types with high accuracy and temporal resolution. Specifically, we generated signaling perturbation atlas across diverse cell types/states through multiplexed sequential combinatorial screens on human pluripotent stem cells. Using the atlas to train IRIS, a neural network-based model, and predicting on mouse embryo scRNAseq atlas, we discovered global features of combinatorial signaling code usage over time, identified biologically meaningful heterogeneity of signaling states within each cell type, and reconstructed signaling histories along diverse cell lineages. We further demonstrated that IRIS greatly accelerates the optimization of stem cell differentiation protocols by drastically reducing the combinatorial space that needs to be tested. This framework leads to the revelation that different cell types share robust signal response signatures, and provides a scalable solution for mapping complex signaling interactions in vivo to guide targeted interventions.

Discrepancies between human and murine model cerebral aneurysms at single-cell resolution

Background

The murine model of cerebral aneurysm (CA) serves as a prevalent tool for investigating the molecular underpinnings of CA. However, the extent to which the CA murine model aligns with that of human remains elusive.

Methods

The present study employed a comprehensive integration and exploration of the single-cell RNA-seq (scRNA-seq) datasets, along with multiple trajectory and gene regulatory network analyses, to investigate the cellular and molecular discrepancies between human and murine model CAs.

Results

The uniform manifold approximation and projection (umap) embedding exhibits that the primary discrepancies between human and murine model CAs reside in the cells of modifiable phenotype, encompassing vascular smooth muscle cell (vSMC), monocyte/macrophage, and neutrophil. The vSMCs from human CA tissue exhibit a fibroblast-like phenotype in comparison to that of murine model. Distinct patterns of neutrophil recruitment are observed in human and murine models, with the former characterized by neutrophil-derived CXCL8 and the latter by monocyte/macrophage-derived CCLs. In addition, macrophages originated from human unruptured CA express higher levels of M2 gene markers. Moreover, the inflammatory status of the CA tissue differs between humans and mouse models, with the former exhibiting a more acute and intense inflammation.

Conclusion

These findings demonstrate subtle but important disparities between human and murine model CAs, and may shed light upon an optimization of murine CA model.

DMGAT: predicting ncRNA-drug resistance associations based on diffusion map and heterogeneous graph attention network

Non-coding RNAs (ncRNAs) play crucial roles in drug resistance and sensitivity, making them important biomarkers and therapeutic targets. However, predicting ncRNA-drug associations is challenging due to issues such as dataset imbalance and sparsity, limiting the identification of robust biomarkers. Existing models often fall short in capturing local and global sequence information, limiting the reliability of predictions. This study introduces DMGAT (diffusion map and heterogeneous graph attention network), a novel deep learning model designed to predict ncRNA-drug associations. DMGAT integrates diffusion maps for sequence embedding, graph convolutional networks for feature extraction, and GAT for heterogeneous information fusion. To address dataset imbalance, the model incorporates sensitivity associations and employs a random forest classifier to select reliable negative samples. DMGAT embeds ncRNA sequences and drug SMILES using the word2vec technique, capturing local and global sequence information. The model constructs a heterogeneous network by combining sequence similarity and Gaussian Interaction Profile kernel similarity, providing a comprehensive representation of ncRNA-drug interactions. Evaluated through five-fold cross-validation on a curated dataset from NoncoRNA and ncDR, DMGAT outperforms seven state-of-the-art methods, achieving the highest area under the receiver operating characteristic curve (0.8964), area under the precision-recall curve (0.8984), recall (0.9576), and F1-score (0.8285). The raw data are released to Zenodo with identifier 13929676. The source code of DMGAT is available at https://github.com/liutingyu0616/DMGAT/tree/main.

Advanced Carriers for Precise Delivery and Therapeutic Mechanisms of Traditional Chinese Medicines: Integrating Spatial Multi-Omics and Delivery Visualization

The complex composition of traditional Chinese medicines (TCMs) has posed challenges for in-depth study and global application, despite their abundance of bioactive compounds that make them valuable resources for disease treatment. To overcome these obstacles, it is essential to modernize TCMs by focusing on precise disease treatment. This involves elucidating the structure-activity relationships within their complex compositions, ensuring accurate in vivo delivery, and monitoring the delivery process. This review discusses the research progress of TCMs in precision disease treatment from three perspectives: spatial multi-omics technology for precision therapeutic activity, carrier systems for precise in vivo delivery, and medical imaging technology for visualizing the delivery process. The aim is to establish a novel research paradigm that advances the precision therapy of TCMs.

Mosquito Cell Atlas: A single-nucleus transcriptomic atlas of the adult Aedes aegypti mosquito

The female mosquito's remarkable ability to hunt humans and transmit pathogens relies on her unique biology. Here, we present the Mosquito Cell Atlas (MCA), a comprehensive single-nucleus RNA sequencing dataset of more than 367,000 nuclei from 19 dissected tissues of adult female and male Aedes aegypti, providing cellular-level resolution of mosquito biology. We identify novel cell types and expand our understanding of sensory neuron organization of chemoreceptors to all sensory tissues. Our analysis uncovers male-specific cells and sexually dimorphic gene expression in the antenna and brain. In female mosquitoes, we find that glial cells in the brain, rather than neurons, undergo the most extensive transcriptional changes following blood feeding. Our findings provide insights into the cellular basis of mosquito behavior and sexual dimorphism. The MCA aims to serve as a resource for the vector biology community, enabling systematic investigation of cell-type specific expression across all mosquito tissues.

Fokker-Planck diffusion maps of microglial transcriptomes reveal radial differentiation into substates associated with Alzheimer's pathology

The identification of microglia subtypes is important for understanding the role of innate immunity in neurodegenerative diseases. Current methods of unsupervised cell type identification assume a small noise-to-signal ratio of transcriptome measurements to produce well-separated cell clusters. However, identification of subtypes can be obscured by gene expression noise, which diminishes the distances in transcriptome space between distinct cell types, blurs boundaries, and reduces reproducibility. Here we use Fokker-Planck (FP) diffusion maps to model cellular differentiation as a stochastic process whereby cells settle into local minima that correspond to cell subtypes, in a potential landscape constructed from transcriptome data using a nearest neighbor graph approach. By applying critical transition fields, we identify individual cells on the verge of transitioning between subtypes, revealing microglial cells in an inactivated, homeostatic state before radially transitioning into various specialized subtypes. Specifically, we show that cells from Alzheimer's disease patients are enriched in a microglia subtype associated to antigen presentation and T-cell recruitment, and are depleted in an anti-inflammatory subtype.

Inferring cell trajectories of spatial transcriptomics via optimal transport analysis

The integration of cell transcriptomics and spatial position to organize differentiation trajectories remains a challenge. Here, we introduce SpaTrack, which leverages optimal transport to reconcile both gene expression and spatial position from spatial transcriptomics into the transition costs, thereby reconstructing cell differentiation. SpaTrack can construct detailed spatial trajectories that reflect the differentiation topology and trace cell dynamics across multiple samples over temporal intervals. To capture the dynamic drivers of differentiation, SpaTrack models cell fate as a function of expression profiles influenced by transcription factors over time. By applying SpaTrack, we successfully disentangle spatiotemporal trajectories of axolotl telencephalon regeneration and mouse midbrain development. Diverse malignant lineages expanding within a primary tumor are uncovered. One lineage, characterized by upregulated epithelial mesenchymal transition, implants at the metastatic site and subsequently colonizes to form a secondary tumor. Overall, SpaTrack efficiently advances trajectory inference from spatial transcriptomics, providing valuable insights into differentiation processes.

A model-based factorization method for scRNA data unveils bifurcating transcriptional modules underlying cell fate determination

Manifold-learning is particularly useful to resolve the complex cellular state space from single-cell RNA sequences. While current manifold-learning methods provide insights into cell fate by inferring graph-based trajectory at cell level, challenges remain to retrieve interpretable biology underlying the diverse cellular states. Here, we described MGPfactXMBD, a model-based manifold-learning framework and capable to factorize complex development trajectories into independent bifurcation processes of gene sets, and thus enables trajectory inference based on relevant features. MGPfactXMBD offers a more nuanced understanding of the biological processes underlying cellular trajectories with potential determinants. When bench-tested across 239 datasets, MGPfactXMBD showed advantages in major quantity-control metrics, such as branch division accuracy and trajectory topology, outperforming most established methods. In real datasets, MGPfactXMBD recovered the critical pathways and cell types in microglia development with experimentally valid regulons and markers. Furthermore, MGPfactXMBD discovered evolutionary trajectories of tumor-associated CD8+ T cells and yielded new subtypes of CD8+ T cells with gene expression signatures significantly predictive of the responses to immune checkpoint inhibitor in independent cohorts. In summary, MGPfactXMBD offers a manifold-learning framework in scRNA-seq data which enables feature selection for specific biological processes and contributing to advance our understanding of biological determination of cell fate.

Multiplexing cortical brain organoids for the longitudinal dissection of developmental traits at single-cell resolution

Dissecting human neurobiology at high resolution and with mechanistic precision requires a major leap in scalability, given the need for experimental designs that include multiple individuals and, prospectively, population cohorts. To lay the foundation for this, we have developed and benchmarked complementary strategies to multiplex brain organoids by pooling cells from different pluripotent stem cell (PSC) lines either during organoid generation (mosaic models) or before single-cell RNA sequencing (scRNA-seq) library preparation (downstream multiplexing). We have also developed a new computational method, SCanSNP, and a consensus call to deconvolve cell identities, overcoming current criticalities in doublets and low-quality cell identification. We validated both multiplexing methods for charting neurodevelopmental trajectories at high resolution, thus linking specific individuals' trajectories to genetic variation. Finally, we modeled their scalability across different multiplexing combinations and showed that mosaic organoids represent an enabling method for high-throughput settings. Together, this multiplexing suite of experimental and computational methods provides a highly scalable resource for brain disease and neurodiversity modeling.

Computing hematopoiesis plasticity in response to genetic mutations and environmental stimulations

Cell plasticity (CP), describing a dynamic cell state, plays a crucial role in maintaining homeostasis during organ morphogenesis, regeneration, and trauma-to-repair biological process. Single-cell-omics datasets provide an unprecedented resource to empower CP analysis. Hematopoiesis offers fertile opportunities to develop quantitative methods for understanding CP. In this study, we generated high-quality lineage-negative single-cell RNA-sequencing datasets under various conditions and introduced a working pipeline named scPlasticity to interrogate naïve and disturbed plasticity of hematopoietic stem and progenitor cells with mutational or environmental challenges. Using embedding methods UMAP or FA, a continuum of hematopoietic development is visually observed in wild type where the pipeline confirms a low proportion of hybrid cells ( P hc , with bias range: 0.4∼0.6) on a transition trajectory. Upon Tet2 mutation, a driver of leukemia, or treatment of DSS, an inducer of colitis, P hc is increased and plasticity of hematopoietic stem and progenitor cells was enhanced. We prioritized several transcription factors and signaling pathways, which are responsible for P hc alterations. In silico perturbation suggests knocking out EGR regulons or pathways of IL-1R1 and β-adrenoreceptor partially reverses P hc promoted by Tet2 mutation and inflammation.

Starfysh integrates spatial transcriptomic and histologic data to reveal heterogeneous tumor-immune hubs

Spatially resolved gene expression profiling provides insight into tissue organization and cell-cell crosstalk; however, sequencing-based spatial transcriptomics (ST) lacks single-cell resolution. Current ST analysis methods require single-cell RNA sequencing data as a reference for rigorous interpretation of cell states, mostly do not use associated histology images and are not capable of inferring shared neighborhoods across multiple tissues. Here we present Starfysh, a computational toolbox using a deep generative model that incorporates archetypal analysis and any known cell type markers to characterize known or new tissue-specific cell states without a single-cell reference. Starfysh improves the characterization of spatial dynamics in complex tissues using histology images and enables the comparison of niches as spatial hubs across tissues. Integrative analysis of primary estrogen receptor (ER)-positive breast cancer, triple-negative breast cancer (TNBC) and metaplastic breast cancer (MBC) tissues led to the identification of spatial hubs with patient- and disease-specific cell type compositions and revealed metabolic reprogramming shaping immunosuppressive hubs in aggressive MBC.

Diffusive topology preserving manifold distances for single-cell data analysis

Manifold learning techniques have emerged as crucial tools for uncovering latent patterns in high-dimensional single-cell data. However, most existing dimensionality reduction methods primarily rely on 2D visualization, which can distort true data relationships and fail to extract reliable biological information. Here, we present DTNE (diffusive topology neighbor embedding), a dimensionality reduction framework that faithfully approximates manifold distance to enhance cellular relationships and dynamics. DTNE constructs a manifold distance matrix using a modified personalized PageRank algorithm, thereby preserving topological structure while enabling diverse single-cell analyses. This approach facilitates distribution-based cellular relationship analysis, pseudotime inference, and clustering within a unified framework. Extensive benchmarking against mainstream algorithms on diverse datasets demonstrates DTNE's superior performance in maintaining geodesic distances and revealing significant biological patterns. Our results establish DTNE as a powerful tool for high-dimensional data analysis in uncovering meaningful biological insights.

Progressive plasticity during colorectal cancer metastasis

As cancers progress, they become increasingly aggressive-metastatic tumours are less responsive to first-line therapies than primary tumours, they acquire resistance to successive therapies and eventually cause death1,2. Mutations are largely conserved between primary and metastatic tumours from the same patients, suggesting that non-genetic phenotypic plasticity has a major role in cancer progression and therapy resistance3-5. However, we lack an understanding of metastatic cell states and the mechanisms by which they transition. Here, in a cohort of biospecimen trios from same-patient normal colon, primary and metastatic colorectal cancer, we show that, although primary tumours largely adopt LGR5+ intestinal stem-like states, metastases display progressive plasticity. Cancer cells lose intestinal cell identities and reprogram into a highly conserved fetal progenitor state before undergoing non-canonical differentiation into divergent squamous and neuroendocrine-like states, a process that is exacerbated in metastasis and by chemotherapy and is associated with poor patient survival. Using matched patient-derived organoids, we demonstrate that metastatic cells exhibit greater cell-autonomous multilineage differentiation potential in response to microenvironment cues compared with their intestinal lineage-restricted primary tumour counterparts. We identify PROX1 as a repressor of non-intestinal lineage in the fetal progenitor state, and show that downregulation of PROX1 licenses non-canonical reprogramming.

Expansion of peripheral cytotoxic CD4+ T cells in Alzheimer's disease: New insights from multi-omics evidence

The significance of the adaptive immune response in Alzheimer's disease (AD) is increasingly recognized. We analyzed scRNA-Seq data from AD patients, revealing a notable rise in CD4 cytotoxic T cells (CD4-CTLs) in peripheral blood mononuclear cells (PBMCs), validated in vivo and in vitro. This rise correlates with cognitive decline in AD patients. We also identified transcription factors TBX21 and MYBL1 as key drivers of CD4-CTL expansion. Further analyses indicate these cells are terminally differentiated, showing clonal expansion, metabolic changes, and unique communication patterns. Mendelian randomization identified risk genes SRGN and ITGB1, suggesting their genetic regulation in CD4-CTLs may contribute to AD. To summarize, our findings characterize the expansion of CD4-CTLs in the PBMCs of AD patients, providing valuable understanding into the possible mechanisms involved in the expansion of CD4-CTLs in AD.

Oncogene activated human breast luminal progenitors contribute basally located myoepithelial cells

BACKGROUND

Basal-like breast cancer originates in luminal progenitors, frequently with an altered PI3K pathway, and focally in close association with genetically altered myoepithelial cells at the site of tumor initiation. The exact trajectory behind this bi-lineage phenomenon remains poorly understood.

METHODS AND RESULTS

Here we used a breast cancer relevant transduction protocol including hTERT, shp16, shp53, and PIK3CAH1047R to immortalize FACS isolated luminal cells, and we identified a candidate multipotent progenitor. Specifically, we identified a keratin 23 (K23)+/ALDH1A3+/CALML5- ductal-like progenitor with the potential to differentiate into CALML5+ lobular-like cells. We found that the apparent luminal phenotype of these oncogene transduced progenitors was metastable giving rise to basal-like cells dependent on culture conditions. In 3D organoid culture and upon transplantation to mice the bipotent progenitor cell line organized into a bi-layered acinus-like structure reminiscent of that of the normal breast gland.

CONCLUSIONS

These findings provide proof of principle that progenitors within the human breast luminal epithelial compartment may serve as a source of correctly positioned myoepithelial cells. This may prove useful in assessing the role of myoepithelial cells in breast tumor progression.

Induced B-Cell Receptor Diversity Predicts PD-1 Blockade Immunotherapy Response

Immune checkpoint inhibitors such as anti-PD-1 antibodies (aPD1) can be effective in treating advanced cancers. However, many patients do not respond and the mechanisms underlying these differences remain incompletely understood. In this study, we profile a cohort of patients with locally-advanced or metastatic basal cell carcinoma undergoing aPD-1 therapy using single-cell RNA sequencing, high-definition spatial transcriptomics in tumors and draining lymph nodes, and spatial immunoreceptor profiling, with long-term clinical follow-up. We find that successful responses to PD-1 inhibition are characterized by an induction of B-cell receptor (BCR) clonal diversity after treatment initiation. These induced BCR clones spatially co-localize with T-cell clones, facilitate their activation, and traffic alongside them between tumor and draining lymph nodes to enhance tumor clearance. Furthermore, we validated aPD1-induced BCR diversity as a predictor of clinical response in a larger cohort of glioblastoma, melanoma, and head and neck squamous cell carcinoma patients, suggesting that this is a generalizable predictor of treatment response across many types of cancers. We discover that pre-treatment tumors harbor a characteristic gene expression signature that portends a higher probability of inducing BCR clonal diversity after aPD-1 therapy, and we develop a machine learning model that predicts PD-1-induced BCR clonal diversity from baseline tumor RNA sequencing. These findings underscore a dynamic role of B cell diversity during immunotherapy, highlighting its importance as a prognostic marker and a potential target for intervention in non-responders.

Complex hierarchical structures analysis in single-cell data with Poincaré deep manifold transformation

Single-cell RNA sequencing (scRNA-seq) offers remarkable insights into cellular development and differentiation by capturing the gene expression profiles of individual cells. The role of dimensionality reduction and visualization in the interpretation of scRNA-seq data has gained widely acceptance. However, current methods face several challenges, including incomplete structure-preserving strategies and high distortion in embeddings, which fail to effectively model complex cell trajectories with multiple branches. To address these issues, we propose the Poincaré deep manifold transformation (PoincaréDMT) method, which maps high-dimensional scRNA-seq data to a hyperbolic Poincaré disk. This approach preserves global structure from a graph Laplacian matrix while achieving local structure correction through a structure module combined with data augmentation. Additionally, PoincaréDMT alleviates batch effects by integrating a batch graph that accounts for batch labels into the low-dimensional embeddings during network training. Furthermore, PoincaréDMT introduces the Shapley additive explanations method based on trained model to identify the important marker genes in specific clusters and cell differentiation process. Therefore, PoincaréDMT provides a unified framework for multiple key tasks essential for scRNA-seq analysis, including trajectory inference, pseudotime inference, batch correction, and marker gene selection. We validate PoincaréDMT through extensive evaluations on both simulated and real scRNA-seq datasets, demonstrating its superior performance in preserving global and local data structures compared to existing methods.

Systematic analysis identifies a connection between spatial and genomic variations of chromatin states

Chromatin states play important roles in the maintenance of cell identities, yet their spatial patterns remain poorly characterized at the organism scale. We developed a systematic approach to analyzing spatial epigenomic data and then applied it to a recently published spatial-CUT&Tag dataset that was obtained from a mouse embryo. We identified a set of spatial genes whose H3K4me3 patterns delineate tissue boundaries. These genes are enriched with tissue-specific transcription factors, and their corresponding genomic loci are marked by broad H3K4me3 domains. Integrative analysis with H3K27me3 profiles showed coordinated spatial transitions across tissue boundaries, which is marked by the continuous shortening of H3K4me3 domains and expansion of H3K27me3 domains. Motif-based analysis identified transcription factors whose activities change significantly during such transitions. Taken together, our systematic analyses reveal a strong connection between the genomic and spatial variations of chromatin states. A record of this paper's transparent peer review process is included in the supplemental information.

scPair: Boosting single cell multimodal analysis by leveraging implicit feature selection and single cell atlases

Multimodal single-cell assays profile multiple sets of features in the same cells and are widely used for identifying and mapping cell states between chromatin and mRNA and linking regulatory elements to target genes. However, the high dimensionality of input features and shallow sequencing depth compared to unimodal assays pose challenges in data analysis. Here we present scPair, a multimodal single-cell data framework that overcomes these challenges by employing an implicit feature selection approach. scPair uses dual encoder-decoder structures trained on paired data to align cell states across modalities and predict features from one modality to another. We demonstrate that scPair outperforms existing methods in accuracy and execution time, and facilitates downstream tasks such as trajectory inference. We further show scPair can augment smaller multimodal datasets with larger unimodal atlases to increase statistical power to identify groups of transcription factors active during different stages of neural differentiation.

Learning glass transition temperatures via dimensionality reduction with data from computer simulations: Polymers as the pilot case

Machine learning methods provide an advanced means for understanding inherent patterns within large and complex datasets. Here, we employ the principal component analysis (PCA) and the diffusion map (DM) techniques to evaluate the glass transition temperature (Tg) from low-dimensional representations of all-atom molecular dynamic simulations of polylactide (PLA) and poly(3-hydroxybutyrate) (PHB). Four molecular descriptors were considered: radial distribution functions (RDFs), mean square displacements (MSDs), relative square displacements (RSDs), and dihedral angles (DAs). By applying Gaussian Mixture Models (GMMs) to analyze the PCA and DM projections and by quantifying their log-likelihoods as a density-based metric, a distinct separation into two populations corresponding to melt and glass states was revealed. This separation enabled the Tg evaluation from a cooling-induced sharp increase in the overlap between log-likelihood distributions at different temperatures. Tg values derived from the RDF and MSD descriptors using DM closely matched the standard computer simulation-based dilatometric and dynamic Tg values for both PLA and PHB models. This was not the case for PCA. The DM-transformed DA and RSD data resulted in Tg values in agreement with experimental ones. Overall, the fusion of atomistic simulations and DMs complemented with the GMMs presents a promising framework for computing Tg and studying the glass transition in a unified way across various molecular descriptors for glass-forming materials.

An integrated transcriptomic cell atlas of human neural organoids

Human neural organoids, generated from pluripotent stem cells in vitro, are useful tools to study human brain development, evolution and disease. However, it is unclear which parts of the human brain are covered by existing protocols, and it has been difficult to quantitatively assess organoid variation and fidelity. Here we integrate 36 single-cell transcriptomic datasets spanning 26 protocols into one integrated human neural organoid cell atlas totalling more than 1.7 million cells1-26. Mapping to developing human brain references27-30 shows primary cell types and states that have been generated in vitro, and estimates transcriptomic similarity between primary and organoid counterparts across protocols. We provide a programmatic interface to browse the atlas and query new datasets, and showcase the power of the atlas to annotate organoid cell types and evaluate new organoid protocols. Finally, we show that the atlas can be used as a diverse control cohort to annotate and compare organoid models of neural disease, identifying genes and pathways that may underlie pathological mechanisms with the neural models. The human neural organoid cell atlas will be useful to assess organoid fidelity, characterize perturbed and diseased states and facilitate protocol development.

Root-specific theanine metabolism and regulation at the single-cell level in tea plants (Camellia sinensis)

Root-synthesized secondary metabolites are critical quality-conferring compounds of foods, plant-derived medicines, and beverages. However, information at a single-cell level on root-specific secondary metabolism remains largely unexplored. L-Theanine, an important quality component of tea, is primarily synthesized in roots, from which it is then transported to new shoots of tea plant. In this study, we present a single-cell RNA sequencing (scRNA-seq)-derived map for the tea plant root, which enabled cell-type-specific analysis of glutamate and ethylamine (two precursors of theanine biosynthesis) metabolism, and theanine biosynthesis, storage, and transport. Our findings support a model in which the theanine biosynthesis pathway occurs via multicellular compartmentation and does not require high co-expression levels of transcription factors and their target genes within the same cell cluster. This study provides novel insights into theanine metabolism and regulation, at the single-cell level, and offers an example for studying root-specific secondary metabolism in other plant systems.

Cross-tissue human fibroblast atlas reveals myofibroblast subtypes with distinct roles in immune modulation

Fibroblasts, known for their functional diversity, play crucial roles in inflammation and cancer. In this study, we conduct comprehensive single-cell RNA sequencing analyses on fibroblast cells from 517 human samples, spanning 11 tissue types and diverse pathological states. We identify distinct fibroblast subpopulations with universal and tissue-specific characteristics. Pathological conditions lead to significant shifts in fibroblast compositions, including the expansion of immune-modulating fibroblasts during inflammation and tissue-remodeling myofibroblasts in cancer. Within the myofibroblast category, we identify four transcriptionally distinct subpopulations originating from different developmental origins, with LRRC15+ myofibroblasts displaying terminally differentiated features. Both LRRC15+ and MMP1+ myofibroblasts demonstrate pro-tumor potential that contribute to the immune-excluded and immune-suppressive tumor microenvironments (TMEs), whereas PI16+ fibroblasts show potential anti-tumor functions in adjacent non-cancerous regions. Fibroblast-subtype compositions define patient subtypes with distinct clinical outcomes. This study advances our understanding of fibroblast biology and suggests potential therapeutic strategies for targeting specific fibroblast subsets in cancer treatment.

Distinct Immune Homeostasis Remodeling Patterns after HLA-Matched and Haploidentical Transplantation

Allogeneic hematopoietic stem cell transplantation (allo-HSCT) is a widely used treatment for a variety of hematopoietic disorders, and also provides a valuable platform for investigating the development of donor-derived immune cells in recipients post-HSCT. The immune system remodels from the donor to the recipient during allo-HSCT. However, little is known about the cell profile alterations as donor homeostasis rebalances to recipient homeostasis following HSCT. Here, multi-omics technology is applied at both the single cell and bulk sample levels, as well as spectrum flow cytometry and fluorescent transgenic mouse models, to dissect the dynamics of the rebalanced homeostatic immune system in recipients after allo-HSCT. The data reveal that all immune subpopulations observed in donors are successfully restored in recipients, though with varying levels of abundance. The remodeling of immune homeostasis exhibits different patterns in HLA-matched and haploidentical HSCT, highlighting distinct biases in T cell reconstitution from the central and peripheral pathways. Furthermore, ZNF683 is critical for maintaining the persistence and quiescence of CD8 T-cell in haploidentical HSCT. The research can serve as a foundation for developing novel strategies to induce immune tolerance.

An integrated approach identifies the molecular underpinnings of murine anterior visceral endoderm migration

The anterior visceral endoderm (AVE) differs from the surrounding visceral endoderm (VE) in its migratory behavior and ability to restrict primitive streak formation to the opposite side of the mouse embryo. To characterize the molecular bases for the unique properties of the AVE, we combined single-cell RNA sequencing of the VE prior to and during AVE migration with phosphoproteomics, high-resolution live-imaging, and short-term lineage labeling and intervention. This identified the transient nature of the AVE with attenuation of "anteriorizing" gene expression as cells migrate and the emergence of heterogeneities in transcriptional states relative to the AVE's position. Using cell communication analysis, we identified the requirement of semaphorin signaling for normal AVE migration. Lattice light-sheet microscopy showed that Sema6D mutants have abnormalities in basal projections and migration speed. These findings point to a tight coupling between transcriptional state and position of the AVE and identify molecular controllers of AVE migration.

Single-cell and spatial transcriptome analyses reveal tertiary lymphoid structures linked to tumour progression and immunotherapy response in nasopharyngeal carcinoma

Tertiary lymphoid structures are immune cell aggregates linked with cancer outcomes, but their interactions with tumour cell aggregates are unclear. Using nasopharyngeal carcinoma as a model, here we analyse single-cell transcriptomes of 343,829 cells from 77 biopsy and blood samples and spatially-resolved transcriptomes of 31,316 spots from 15 tumours to decipher their components and interactions with tumour cell aggregates. We identify essential cell populations in tertiary lymphoid structure, including CXCL13+ cancer-associated fibroblasts, stem-like CXCL13+CD8+ T cells, and B and T follicular helper cells. Our study shows that germinal centre reaction matures plasma cells. These plasma cells intersperse with tumour cell aggregates, promoting apoptosis of EBV-related malignant cells and enhancing immunotherapy response. CXCL13+ cancer-associated fibroblasts promote B cell adhesion and antibody production, activating CXCL13+CD8+ T cells that become exhausted in tumour cell aggregates. Tertiary lymphoid structure-related cell signatures correlate with prognosis and PD-1 blockade response, offering insights for therapeutic strategies in cancers.

StableMate: a statistical method to select stable predictors in omics data

Identifying statistical associations between biological variables is crucial to understanding molecular mechanisms. Most association studies are based on correlation or linear regression analyses, but the identified associations often lack reproducibility and interpretability due to the complexity and variability of omics datasets, making it difficult to translate associations into meaningful biological hypotheses. We developed StableMate, a regression framework, to address these challenges through a process of variable selection across heterogeneous datasets. Given datasets from different environments, such as experimental batches, StableMate selects environment-agnostic (stable) and environment-specific predictors in predicting the response of interest. Stable predictors represent robust functional dependencies with the response, and can be used to build regression models that make generalizable predictions in unseen environments. We applied StableMate to (i) RNA sequencing data of breast cancer to discover genes that consistently predict estrogen receptor expression across disease status; (ii) metagenomics data to identify microbial signatures that show persistent association with colon cancer across study cohorts; and (iii) single-cell RNA sequencing data of glioblastoma to discern signature genes associated with the development of pro-tumour microglia regardless of cell location. Our case studies demonstrate that StableMate is adaptable to regression and classification analyses and achieves comprehensive characterization of biological systems for different omics data types.

Pan-cancer single-cell dissection reveals phenotypically distinct B cell subtypes

Characterizing the compositional and phenotypic characteristics of tumor-infiltrating B cells (TIBs) is important for advancing our understanding of their role in cancer development. Here, we establish a comprehensive resource of human B cells by integrating single-cell RNA sequencing data of B cells from 649 patients across 19 major cancer types. We demonstrate substantial heterogeneity in their total abundance and subtype composition and observe immunoglobulin G (IgG)-skewness of antibody-secreting cell isotypes. Moreover, we identify stress-response memory B cells and tumor-associated atypical B cells (TAABs), two tumor-enriched subpopulations with prognostic potential, shared in a pan-cancer manner. In particular, TAABs, characterized by a high clonal expansion level and proliferative capacity as well as by close interactions with activated CD4 T cells in tumors, are predictive of immunotherapy response. Our integrative resource depicts distinct clinically relevant TIB subsets, laying a foundation for further exploration of functional commonality and diversity of B cells in cancer.

Tumour vasculature at single-cell resolution

Tumours can obtain nutrients and oxygen required to progress and metastasize through the blood supply1. Inducing angiogenesis involves the sprouting of established vessel beds and their maturation into an organized network2,3. Here we generate a comprehensive atlas of tumour vasculature at single-cell resolution, encompassing approximately 200,000 cells from 372 donors representing 31 cancer types. Trajectory inference suggested that tumour angiogenesis was initiated from venous endothelial cells and extended towards arterial endothelial cells. As neovascularization elongates (through angiogenic stages SI, SII and SIII), APLN+ tip cells at the SI stage (APLN+ TipSI) advanced to TipSIII cells with increased Notch signalling. Meanwhile, stalk cells, following tip cells, transitioned from high chemokine expression to elevated TEK (also known as Tie2) expression. Moreover, APLN+ TipSI cells not only were associated with disease progression and poor prognosis but also hold promise for predicting response to anti-VEGF therapy. Lymphatic endothelial cells demonstrated two distinct differentiation lineages: one responsible for lymphangiogenesis and the other involved in antigen presentation. In pericytes, endoplasmic reticulum stress was associated with the proangiogenic BASP1+ matrix-producing pericytes. Furthermore, intercellular communication analysis showed that neovascular endothelial cells could shape an immunosuppressive microenvironment conducive to angiogenesis. This study depicts the complexity of tumour vasculature and has potential clinical significance for anti-angiogenic therapy.

Polyploidisation pleiotropically buffers ageing in hepatocytes

BACKGROUND & AIMS

Polyploidy in hepatocytes has been proposed as a genetic mechanism to buffer against transcriptional dysregulation. Here, we aim to demonstrate the role of polyploidy in modulating gene regulatory networks in hepatocytes during ageing.

METHODS

We performed single-nucleus RNA sequencing in hepatocyte nuclei of different ploidy levels isolated from young and old wild-type mice. Changes in the gene expression and regulatory network were compared to three independent strains that were haploinsufficient for HNF4A, CEBPA or CTCF, representing non-deleterious perturbations. Phenotypic characteristics of the liver section were additionally evaluated histologically, whereas the genomic allele composition of hepatocytes was analysed by BaseScope.

RESULTS

We observed that ageing in wild-type mice results in nuclei polyploidy and a marked increase in steatosis. Haploinsufficiency of liver-specific master regulators (HFN4A or CEBPA) results in the enrichment of hepatocytes with tetraploid nuclei at a young age, affecting the genomic regulatory network, and dramatically suppressing ageing-related steatosis tissue wide. Notably, these phenotypes are not the result of subtle disruption to liver-specific transcriptional networks, since haploinsufficiency in the CTCF insulator protein resulted in the same phenotype. Further quantification of genotypes of tetraploid hepatocytes in young and old HFN4A-haploinsufficient mice revealed that during ageing, tetraploid hepatocytes lead to the selection of wild-type alleles, restoring non-deleterious genetic perturbations.

CONCLUSIONS

Our results suggest a model whereby polyploidisation leads to fundamentally different cell states. Polyploid conversion enables pleiotropic buffering against age-related decline via non-random allelic segregation to restore a wild-type genome.

IMPACT AND IMPLICATIONS

The functional role of hepatocyte polyploidisation during ageing is poorly understood. Using single-nucleus RNA sequencing and BaseScope approaches, we have studied ploidy dynamics during ageing in murine livers with non-deleterious genetic perturbations. We have identified that hepatocytes present different cellular states and the ability to buffer ageing-associated dysfunctions. Tetraploid nuclei exhibit robust transcriptional networks and are better adapted to genomically overcome perturbations. Novel therapeutic interventions aimed at attenuating age-related changes in tissue function could be exploited by manipulation of ploidy dynamics during chronic liver conditions.

Pan-cancer profiling of tumor-infiltrating natural killer cells through transcriptional reference mapping

The functional diversity of natural killer (NK) cell repertoires stems from differentiation, homeostatic, receptor-ligand interactions and adaptive-like responses to viral infections. In the present study, we generated a single-cell transcriptional reference map of healthy human blood- and tissue-derived NK cells, with temporal resolution and fate-specific expression of gene-regulatory networks defining NK cell differentiation. Transfer learning facilitated incorporation of tumor-infiltrating NK cell transcriptomes (39 datasets, 7 solid tumors, 427 patients) into the reference map to analyze tumor microenvironment (TME)-induced perturbations. Of the six functionally distinct NK cell states identified, a dysfunctional stressed CD56bright state susceptible to TME-induced immunosuppression and a cytotoxic TME-resistant effector CD56dim state were commonly enriched across tumor types, the ratio of which was predictive of patient outcome in malignant melanoma and osteosarcoma. This resource may inform the design of new NK cell therapies and can be extended through transfer learning to interrogate new datasets from experimental perturbations or disease conditions.

Diffusion on PCA-UMAP Manifold: The Impact of Data Structure Preservation to Denoise High-Dimensional Single-Cell RNA Sequencing Data

Single-cell transcriptomics (scRNA-seq) is revolutionizing biological research, yet it faces challenges such as inefficient transcript capture and noise. To address these challenges, methods like neighbor averaging or graph diffusion are used. These methods often rely on k-nearest neighbor graphs from low-dimensional manifolds. However, scRNA-seq data suffer from the 'curse of dimensionality', leading to the over-smoothing of data when using imputation methods. To overcome this, sc-PHENIX employs a PCA-UMAP diffusion method, which enhances the preservation of data structures and allows for a refined use of PCA dimensions and diffusion parameters (e.g., k-nearest neighbors, exponentiation of the Markov matrix) to minimize noise introduction. This approach enables a more accurate construction of the exponentiated Markov matrix (cell neighborhood graph), surpassing methods like MAGIC. sc-PHENIX significantly mitigates over-smoothing, as validated through various scRNA-seq datasets, demonstrating improved cell phenotype representation. Applied to a multicellular tumor spheroid dataset, sc-PHENIX identified known extreme phenotype states, showcasing its effectiveness. sc-PHENIX is open-source and available for use and modification.

Supervised discovery of interpretable gene programs from single-cell data

Factor analysis decomposes single-cell gene expression data into a minimal set of gene programs that correspond to processes executed by cells in a sample. However, matrix factorization methods are prone to technical artifacts and poor factor interpretability. We address these concerns with Spectra, an algorithm that combines user-provided gene programs with the detection of novel programs that together best explain expression covariation. Spectra incorporates existing gene sets and cell-type labels as prior biological information, explicitly models cell type and represents input gene sets as a gene-gene knowledge graph using a penalty function to guide factorization toward the input graph. We show that Spectra outperforms existing approaches in challenging tumor immune contexts, as it finds factors that change under immune checkpoint therapy, disentangles the highly correlated features of CD8+ T cell tumor reactivity and exhaustion, finds a program that explains continuous macrophage state changes under therapy and identifies cell-type-specific immune metabolic programs.

Quantifying cell-state densities in single-cell phenotypic landscapes using Mellon

Cell-state density characterizes the distribution of cells along phenotypic landscapes and is crucial for unraveling the mechanisms that drive diverse biological processes. Here, we present Mellon, an algorithm for estimation of cell-state densities from high-dimensional representations of single-cell data. We demonstrate Mellon's efficacy by dissecting the density landscape of differentiating systems, revealing a consistent pattern of high-density regions corresponding to major cell types intertwined with low-density, rare transitory states. We present evidence implicating enhancer priming and the activation of master regulators in emergence of these transitory states. Mellon offers the flexibility to perform temporal interpolation of time-series data, providing a detailed view of cell-state dynamics during developmental processes. Mellon facilitates density estimation across various single-cell data modalities, scaling linearly with the number of cells. Our work underscores the importance of cell-state density in understanding the differentiation processes, and the potential of Mellon to provide insights into mechanisms guiding biological trajectories.

Building and analyzing metacells in single-cell genomics data

The advent of high-throughput single-cell genomics technologies has fundamentally transformed biological sciences. Currently, millions of cells from complex biological tissues can be phenotypically profiled across multiple modalities. The scaling of computational methods to analyze and visualize such data is a constant challenge, and tools need to be regularly updated, if not redesigned, to cope with ever-growing numbers of cells. Over the last few years, metacells have been introduced to reduce the size and complexity of single-cell genomics data while preserving biologically relevant information and improving interpretability. Here, we review recent studies that capitalize on the concept of metacells-and the many variants in nomenclature that have been used. We further outline how and when metacells should (or should not) be used to analyze single-cell genomics data and what should be considered when analyzing such data at the metacell level. To facilitate the exploration of metacells, we provide a comprehensive tutorial on the construction and analysis of metacells from single-cell RNA-seq data ( https://github.com/GfellerLab/MetacellAnalysisTutorial ) as well as a fully integrated pipeline to rapidly build, visualize and evaluate metacells with different methods ( https://github.com/GfellerLab/MetacellAnalysisToolkit ).

CellRank 2: unified fate mapping in multiview single-cell data

Single-cell RNA sequencing allows us to model cellular state dynamics and fate decisions using expression similarity or RNA velocity to reconstruct state-change trajectories; however, trajectory inference does not incorporate valuable time point information or utilize additional modalities, whereas methods that address these different data views cannot be combined or do not scale. Here we present CellRank 2, a versatile and scalable framework to study cellular fate using multiview single-cell data of up to millions of cells in a unified fashion. CellRank 2 consistently recovers terminal states and fate probabilities across data modalities in human hematopoiesis and endodermal development. Our framework also allows combining transitions within and across experimental time points, a feature we use to recover genes promoting medullary thymic epithelial cell formation during pharyngeal endoderm development. Moreover, we enable estimating cell-specific transcription and degradation rates from metabolic-labeling data, which we apply to an intestinal organoid system to delineate differentiation trajectories and pinpoint regulatory strategies.

Single-cell-resolved interspecies comparison shows a shared inflammatory axis and a dominant neutrophil-endothelial program in severe COVID-19

A key issue for research on COVID-19 pathogenesis is the lack of biopsies from patients and of samples at the onset of infection. To overcome these hurdles, hamsters were shown to be useful models for studying this disease. Here, we further leverage the model to molecularly survey the disease progression from time-resolved single-cell RNA sequencing data collected from healthy and severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)-infected Syrian and Roborovski hamster lungs. We compare our data to human COVID-19 studies, including bronchoalveolar lavage, nasal swab, and postmortem lung tissue, and identify a shared axis of inflammation dominated by macrophages, neutrophils, and endothelial cells, which we show to be transient in Syrian and terminal in Roborovski hamsters. Our data suggest that, following SARS-CoV-2 infection, commitment to a type 1- or type 3-biased immunity determines moderate versus severe COVID-19 outcomes, respectively.

Recent Developments in Machine Learning for Mass Spectrometry

Statistical analysis and modeling of mass spectrometry (MS) data have a long and rich history with several modern MS-based applications using statistical and chemometric methods. Recently, machine learning (ML) has experienced a renaissance due to advents in computational hardware and the development of new algorithms for artificial neural networks (ANN) and deep learning architectures. Moreover, recent successes of new ANN and deep learning architectures in several areas of science, engineering, and society have further strengthened the ML field. Importantly, modern ML methods and architectures have enabled new approaches for tasks related to MS that are now widely adopted in several popular MS-based subdisciplines, such as mass spectrometry imaging and proteomics. Herein, we aim to provide an introductory summary of the practical aspects of ML methodology relevant to MS. Additionally, we seek to provide an up-to-date review of the most recent developments in ML integration with MS-based techniques while also providing critical insights into the future direction of the field.

Identification of a myofibroblast differentiation program during neonatal lung development

Alveologenesis is the final stage of lung development in which the internal surface area of the lung is increased to facilitate efficient gas exchange in the mature organism. The first phase of alveologenesis involves the formation of septal ridges (secondary septae) and the second phase involves thinning of the alveolar septa. Within secondary septa, mesenchymal cells include a transient population of alveolar myofibroblasts (MyoFBs) and a stable but poorly described population of lipid-rich cells that have been referred to as lipofibroblasts or matrix fibroblasts (MatFBs). Using a unique Fgf18CreER lineage trace mouse line, cell sorting, single-cell RNA sequencing and primary cell culture, we have identified multiple subtypes of mesenchymal cells in the neonatal lung, including an immature progenitor cell that gives rise to mature MyoFB. We also show that the endogenous and targeted ROSA26 locus serves as a sensitive reporter for MyoFB maturation. These studies identify a MyoFB differentiation program that is distinct from other mesenchymal cell types and increases the known repertoire of mesenchymal cell types in the neonatal lung.

Prediction of protein-RNA interactions from single-cell transcriptomic data

Proteins are crucial in regulating every aspect of RNA life, yet understanding their interactions with coding and noncoding RNAs remains limited. Experimental studies are typically restricted to a small number of cell lines and a limited set of RNA-binding proteins (RBPs). Although computational methods based on physico-chemical principles can predict protein-RNA interactions accurately, they often lack the ability to consider cell-type-specific gene expression and the broader context of gene regulatory networks (GRNs). Here, we assess the performance of several GRN inference algorithms in predicting protein-RNA interactions from single-cell transcriptomic data, and propose a pipeline, called scRAPID (single-cell transcriptomic-based RnA Protein Interaction Detection), that integrates these methods with the catRAPID algorithm, which can identify direct physical interactions between RBPs and RNA molecules. Our approach demonstrates that RBP-RNA interactions can be predicted from single-cell transcriptomic data, with performances comparable or superior to those achieved for the well-established task of inferring transcription factor-target interactions. The incorporation of catRAPID significantly enhances the accuracy of identifying interactions, particularly with long noncoding RNAs, and enables the identification of hub RBPs and RNAs. Additionally, we show that interactions between RBPs can be detected based on their inferred RNA targets. The software is freely available at https://github.com/tartaglialabIIT/scRAPID.

In vivo single-cell high-dimensional mass cytometry analysis to track the interactions between Klebsiella pneumoniae and myeloid cells

In vivo single-cell approaches have transformed our understanding of the immune populations in tissues. Mass cytometry (CyTOF), that combines the resolution of mass spectrometry with the ability to conduct multiplexed measurements of cell molecules at the single cell resolution, has enabled to resolve the diversity of immune cell subsets, and their heterogeneous functionality. Here we assess the feasibility of taking CyTOF one step further to immuno profile cells while tracking their interactions with bacteria, a method we term Bac-CyTOF. We focus on the pathogen Klebsiella pneumoniae interrogating the pneumonia mouse model. Using Bac-CyTOF, we unveil the atlas of immune cells of mice infected with a K. pneumoniae hypervirulent strain. The atlas is characterized by a decrease in the populations of alveolar and monocyte-derived macrophages. Conversely, neutrophils, and inflammatory monocytes are characterized by an increase in the subpopulations expressing markers of less active cells such as the immune checkpoint PD-L1. These are the cells infected. We show that the type VI secretion system (T6SS) contributes to shape the lung immune landscape. The T6SS governs the interaction with monocytes/macrophages by shifting Klebsiella from alveolar macrophages to interstitial macrophages and limiting the infection of inflammatory monocytes. The lack of T6SS results in an increase of cells expressing markers of active cells, and a decrease in the subpopulations expressing PD-L1. By probing Klebsiella, and Acinetobacter baumannii strains with limited ability to survive in vivo, we uncover that a heightened recruitment of neutrophils, and relative high levels of alveolar macrophages and eosinophils and the recruitment of a characteristic subpopulation of neutrophils are features of mice clearing infections. We leverage Bac-CyTOF-generated knowledge platform to investigate the role of the DNA sensor STING in Klebsiella infections. sting-/- infected mice present features consistent with clearing the infection including the reduced levels of PD-L1. STING absence facilitates Klebsiella clearance.

Data-driven selection of analysis decisions in single-cell RNA-seq trajectory inference

Single-cell RNA sequencing (scRNA-seq) experiments have become instrumental in developmental and differentiation studies, enabling the profiling of cells at a single or multiple time-points to uncover subtle variations in expression profiles reflecting underlying biological processes. Benchmarking studies have compared many of the computational methods used to reconstruct cellular dynamics; however, researchers still encounter challenges in their analysis due to uncertainty with respect to selecting the most appropriate methods and parameters. Even among universal data processing steps used by trajectory inference methods such as feature selection and dimension reduction, trajectory methods' performances are highly dataset-specific. To address these challenges, we developed Escort, a novel framework for evaluating a dataset's suitability for trajectory inference and quantifying trajectory properties influenced by analysis decisions. Escort evaluates the suitability of trajectory analysis and the combined effects of processing choices using trajectory-specific metrics. Escort navigates single-cell trajectory analysis through these data-driven assessments, reducing uncertainty and much of the decision burden inherent to trajectory inference analyses. Escort is implemented in an accessible R package and R/Shiny application, providing researchers with the necessary tools to make informed decisions during trajectory analysis and enabling new insights into dynamic biological processes at single-cell resolution.

Single-cell multi-omics analysis identifies context-specific gene regulatory gates and mechanisms

There is a growing interest in inferring context specific gene regulatory networks from single-cell RNA sequencing (scRNA-seq) data. This involves identifying the regulatory relationships between transcription factors (TFs) and genes in individual cells, and then characterizing these relationships at the level of specific cell types or cell states. In this study, we introduce scGATE (single-cell gene regulatory gate) as a novel computational tool for inferring TF-gene interaction networks and reconstructing Boolean logic gates involving regulatory TFs using scRNA-seq data. In contrast to current Boolean models, scGATE eliminates the need for individual formulations and likelihood calculations for each Boolean rule (e.g. AND, OR, XOR). By employing a Bayesian framework, scGATE infers the Boolean rule after fitting the model to the data, resulting in significant reductions in time-complexities for logic-based studies. We have applied assay for transposase-accessible chromatin with sequencing (scATAC-seq) data and TF DNA binding motifs to filter out non-relevant TFs in gene regulations. By integrating single-cell clustering with these external cues, scGATE is able to infer context specific networks. The performance of scGATE is evaluated using synthetic and real single-cell multi-omics data from mouse tissues and human blood, demonstrating its superiority over existing tools for reconstructing TF-gene networks. Additionally, scGATE provides a flexible framework for understanding the complex combinatorial and cooperative relationships among TFs regulating target genes by inferring Boolean logic gates among them.

Plant biotechnology research with single-cell transcriptome: recent advancements and prospects

KEY MESSAGE

Single-cell transcriptomic techniques have emerged as powerful tools in plant biology, offering high-resolution insights into gene expression at the individual cell level. This review highlights the rapid expansion of single-cell technologies in plants, their potential in understanding plant development, and their role in advancing plant biotechnology research. Single-cell techniques have emerged as powerful tools to enhance our understanding of biological systems, providing high-resolution transcriptomic analysis at the single-cell level. In plant biology, the adoption of single-cell transcriptomics has seen rapid expansion of available technologies and applications. This review article focuses on the latest advancements in the field of single-cell transcriptomic in plants and discusses the potential role of these approaches in plant development and expediting plant biotechnology research in the near future. Furthermore, inherent challenges and limitations of single-cell technology are critically examined to overcome them and enhance our knowledge and understanding.

A fast, scalable and versatile tool for analysis of single-cell omics data

Single-cell omics technologies have revolutionized the study of gene regulation in complex tissues. A major computational challenge in analyzing these datasets is to project the large-scale and high-dimensional data into low-dimensional space while retaining the relative relationships between cells. This low dimension embedding is necessary to decompose cellular heterogeneity and reconstruct cell-type-specific gene regulatory programs. Traditional dimensionality reduction techniques, however, face challenges in computational efficiency and in comprehensively addressing cellular diversity across varied molecular modalities. Here we introduce a nonlinear dimensionality reduction algorithm, embodied in the Python package SnapATAC2, which not only achieves a more precise capture of single-cell omics data heterogeneities but also ensures efficient runtime and memory usage, scaling linearly with the number of cells. Our algorithm demonstrates exceptional performance, scalability and versatility across diverse single-cell omics datasets, including single-cell assay for transposase-accessible chromatin using sequencing, single-cell RNA sequencing, single-cell Hi-C and single-cell multi-omics datasets, underscoring its utility in advancing single-cell analysis.

Dictionary learning for integrative, multimodal and scalable single-cell analysis

Mapping single-cell sequencing profiles to comprehensive reference datasets provides a powerful alternative to unsupervised analysis. However, most reference datasets are constructed from single-cell RNA-sequencing data and cannot be used to annotate datasets that do not measure gene expression. Here we introduce 'bridge integration', a method to integrate single-cell datasets across modalities using a multiomic dataset as a molecular bridge. Each cell in the multiomic dataset constitutes an element in a 'dictionary', which is used to reconstruct unimodal datasets and transform them into a shared space. Our procedure accurately integrates transcriptomic data with independent single-cell measurements of chromatin accessibility, histone modifications, DNA methylation and protein levels. Moreover, we demonstrate how dictionary learning can be combined with sketching techniques to improve computational scalability and harmonize 8.6 million human immune cell profiles from sequencing and mass cytometry experiments. Our approach, implemented in version 5 of our Seurat toolkit ( http://www.satijalab.org/seurat ), broadens the utility of single-cell reference datasets and facilitates comparisons across diverse molecular modalities.