Optimal-Transport Analysis of Single-Cell Gene Expression Identifies Developmental Trajectories in Reprogramming

Cytokines drive the formation of memory-like NK cell subsets via epigenetic rewiring and transcriptional regulation

Activation of natural killer (NK) cells with the cytokines interleukin-12 (IL-12), IL-15, and IL-18 induces their differentiation into memory-like (ML) NK cells; however, the underlying epigenetic and transcriptional mechanisms are unclear. By combining ATAC-seq, CITE-seq, and functional analyses, we discovered that IL-12/15/18 activation results in two main human NK fates: reprogramming into enriched memory-like (eML) NK cells or priming into effector conventional NK (effcNK) cells. eML NK cells had distinct transcriptional and epigenetic profiles and enhanced function, whereas effcNK cells resembled cytokine-primed cNK cells. Two transcriptionally discrete subsets of eML NK cells were also identified, eML-1 and eML-2, primarily arising from CD56bright or CD56dim mature NK cell subsets, respectively. Furthermore, these eML subsets were evident weeks after transfer of IL-12/15/18-activated NK cells into patients with cancer. Our findings demonstrate that NK cell activation with IL-12/15/18 results in previously unappreciated diverse cellular fates and identifies new strategies to enhance NK therapies.

Computational single-cell methods for predicting cancer risk

Despite recent biotechnological breakthroughs, cancer risk prediction remains a formidable computational and experimental challenge. Addressing it is critical in order to improve prevention, early detection and survival rates. Here, I briefly summarize some key emerging theoretical and computational challenges as well as recent computational advances that promise to help realize the goals of cancer-risk prediction. The focus is on computational strategies based on single-cell data, in particular on bottom-up network modeling approaches that aim to estimate cancer stemness and dedifferentiation at single-cell resolution from a systems-biological perspective. I will describe two promising methods, a tissue and cell-lineage independent one based on the concept of diffusion network entropy, and a tissue and cell-lineage specific one that uses transcription factor regulons. Application of these tools to single-cell and single-nucleus RNA-seq data from stages prior to invasive cancer reveal that they can successfully delineate the heterogeneous inter-cellular cancer-risk landscape, identifying those cells that are more likely to turn cancerous. Bottom-up systems biological modeling of single-cell omic data is a novel computational analysis paradigm that promises to facilitate the development of preventive, early detection and cancer-risk prediction strategies.

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.

Optimal transport for automatic alignment of untargeted metabolomic data

Untargeted metabolomic profiling through liquid chromatography-mass spectrometry (LC-MS) measures a vast array of metabolites within biospecimens, advancing drug development, disease diagnosis, and risk prediction. However, the low throughput of LC-MS poses a major challenge for biomarker discovery, annotation, and experimental comparison, necessitating the merging of multiple datasets. Current data pooling methods encounter practical limitations due to their vulnerability to data variations and hyperparameter dependence. Here, we introduce GromovMatcher, a flexible and user-friendly algorithm that automatically combines LC-MS datasets using optimal transport. By capitalizing on feature intensity correlation structures, GromovMatcher delivers superior alignment accuracy and robustness compared to existing approaches. This algorithm scales to thousands of features requiring minimal hyperparameter tuning. Manually curated datasets for validating alignment algorithms are limited in the field of untargeted metabolomics, and hence we develop a dataset split procedure to generate pairs of validation datasets to test the alignments produced by GromovMatcher and other methods. Applying our method to experimental patient studies of liver and pancreatic cancer, we discover shared metabolic features related to patient alcohol intake, demonstrating how GromovMatcher facilitates the search for biomarkers associated with lifestyle risk factors linked to several cancer types.

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.

Decoding changes in tumor-infiltrating leukocytes through dynamic experimental models and single-cell technologies

The ability to characterize immune cells and explore the molecular interactions that govern their functions has never been greater, fueled in recent years by the revolutionary advance of single-cell analysis platforms. However, precisely how immune cells respond to different stimuli and where differentiation processes and effector functions operate remain incompletely understood. Inferring cellular fate within single-cell transcriptomic analyses is now omnipresent, despite the assumptions typically required in such analyses. Recently developed experimental models support dynamic analyses of the immune response, providing insights into the temporal changes that occur within cells and the tissues in which such transitions occur. Here we will review these approaches and discuss how these can be combined with single-cell technologies to develop a deeper understanding of the immune responses that should support the development of better therapeutic options for patients.

Patterning and folding of intestinal villi by active mesenchymal dewetting

Tissue folds are structural motifs critical to organ function. In the intestine, bending of a flat epithelium into a periodic pattern of folds gives rise to villi, finger-like protrusions that enable nutrient absorption. However, the molecular and mechanical processes driving villus morphogenesis remain unclear. Here, we identify an active mechanical mechanism that simultaneously patterns and folds the intestinal epithelium to initiate villus formation. At the cellular level, we find that PDGFRA+ subepithelial mesenchymal cells generate myosin II-dependent forces sufficient to produce patterned curvature in neighboring tissue interfaces. This symmetry-breaking process requires altered cell and extracellular matrix interactions that are enabled by matrix metalloproteinase-mediated tissue fluidization. Computational models, together with in vitro and in vivo experiments, revealed that these cellular features manifest at the tissue level as differences in interfacial tensions that promote mesenchymal aggregation and interface bending through a process analogous to the active dewetting of a thin liquid film.

Wasserstein Wormhole: Scalable Optimal Transport Distance with Transformers

Optimal transport (OT) and the related Wasserstein metric ( W ) are powerful and ubiquitous tools for comparing distributions. However, computing pairwise Wasserstein distances rapidly becomes intractable as cohort size grows. An attractive alternative would be to find an embedding space in which pairwise Euclidean distances map to OT distances, akin to standard multidimensional scaling (MDS). We present Wasserstein Wormhole, a transformer-based autoencoder that embeds empirical distributions into a latent space wherein Euclidean distances approximate OT distances. Extending MDS theory, we show that our objective function implies a bound on the error incurred when embedding non-Euclidean distances. Empirically, distances between Wormhole embeddings closely match Wasserstein distances, enabling linear time computation of OT distances. Along with an encoder that maps distributions to embeddings, Wasserstein Wormhole includes a decoder that maps embeddings back to distributions, allowing for operations in the embedding space to generalize to OT spaces, such as Wasserstein barycenter estimation and OT interpolation. By lending scalability and interpretability to OT approaches, Wasserstein Wormhole unlocks new avenues for data analysis in the fields of computational geometry and single-cell biology. Software is available at http://wassersteinwormhole.readthedocs.io/en/latest/.

Single-Cell and Spatial Transcriptome Analysis of Dermal Fibroblast Development in Perinatal Mouse Skin: Dynamic Lineage Differentiation and Key Driver Genes

Several single-cell RNA studies of developing mouse skin have elucidated the molecular and cellular processes involved in skin development. However, they have primarily focused on either the fetal or early postnatal period, leaving a gap in our understanding of skin development. In this study, we conducted a comprehensive time-series analysis by combining single-cell RNA-sequencing datasets collected at different stages of development (embryonic days 13.5, 14.5, and 16.5 and postnatal days 0, 2, and 4) and validated our findings through multipanel in situ spatial transcriptomics. Our analysis indicated that embryonic fibroblasts exhibit heterogeneity from a very early stage and that the rapid determination of each lineage occurs within days after birth. The expression of putative key driver genes, including Hey1, Ebf1, Runx3, and Sox11 for the dermal papilla trajectory; Lrrc15 for the dermal sheath trajectory; Zfp536 and Nrn1 for the papillary fibroblast trajectory; and Lrrn4cl and Mfap5 for the fascia fibroblast trajectory, was detected in the corresponding, spatially identified cell types. Finally, cell-to-cell interaction analysis indicated that the dermal papilla lineage is the primary source of the noncanonical Wnt pathway during skin development. Together, our study provides a transcriptomic reference for future research in the field of skin development and regeneration.

Single-cell transcriptomics reveals evolutionary reconfiguration of embryonic cell fate specification in the sea urchin Heliocidaris erythrogramma

Altered regulatory interactions during development likely underlie a large fraction of phenotypic diversity within and between species, yet identifying specific evolutionary changes remains challenging. Analysis of single-cell developmental transcriptomes from multiple species provides a powerful framework for unbiased identification of evolutionary changes in developmental mechanisms. Here, we leverage a "natural experiment" in developmental evolution in sea urchins, where a major life history switch recently evolved in the lineage leading to Heliocidaris erythrogramma, precipitating extensive changes in early development. Comparative analyses of scRNA-seq developmental time courses from H. erythrogramma and Lytechinus variegatus (representing the derived and ancestral states respectively) reveals numerous evolutionary changes in embryonic patterning. The earliest cell fate specification events, and the primary signaling center are co-localized in the ancestral dGRN but remarkably, in H. erythrogramma they are spatially and temporally separate. Fate specification and differentiation are delayed in most embryonic cell lineages, although in some cases, these processes are conserved or even accelerated. Comparative analysis of regulator-target gene co-expression is consistent with many specific interactions being preserved but delayed in H. erythrogramma, while some otherwise widely conserved interactions have likely been lost. Finally, specific patterning events are directly correlated with evolutionary changes in larval morphology, suggesting that they are directly tied to the life history shift. Together, these findings demonstrate that comparative scRNA-seq developmental time courses can reveal a diverse set of evolutionary changes in embryonic patterning and provide an efficient way to identify likely candidate regulatory interactions for subsequent experimental validation.

A logic-incorporated gene regulatory network deciphers principles in cell fate decisions

Organisms utilize gene regulatory networks (GRN) to make fate decisions, but the regulatory mechanisms of transcription factors (TF) in GRNs are exceedingly intricate. A longstanding question in this field is how these tangled interactions synergistically contribute to decision-making procedures. To comprehensively understand the role of regulatory logic in cell fate decisions, we constructed a logic-incorporated GRN model and examined its behavior under two distinct driving forces (noise-driven and signal-driven). Under the noise-driven mode, we distilled the relationship among fate bias, regulatory logic, and noise profile. Under the signal-driven mode, we bridged regulatory logic and progression-accuracy trade-off, and uncovered distinctive trajectories of reprogramming influenced by logic motifs. In differentiation, we characterized a special logic-dependent priming stage by the solution landscape. Finally, we applied our findings to decipher three biological instances: hematopoiesis, embryogenesis, and trans-differentiation. Orthogonal to the classical analysis of expression profile, we harnessed noise patterns to construct the GRN corresponding to fate transition. Our work presents a generalizable framework for top-down fate-decision studies and a practical approach to the taxonomy of cell fate decisions.

Unlocking trophectoderm mysteries: In vivo and in vitro perspectives on human and mouse trophoblast fate induction

In recent years, the pursuit of inducing the trophoblast stem cell (TSC) state has gained prominence as a compelling research objective, illuminating the establishment of the trophoblast lineage and unlocking insights into early embryogenesis. In this review, we examine how advancements in diverse technologies, including in vivo time course transcriptomics, cellular reprogramming to TSC state, chemical induction of totipotent stem-cell-like state, and stem-cell-based embryo-like structures, have enriched our insights into the intricate molecular mechanisms and signaling pathways that define the mouse and human trophectoderm/TSC states. We delve into disparities between mouse and human trophectoderm/TSC fate establishment, with a special emphasis on the intriguing role of pluripotency in this context. Additionally, we re-evaluate recent findings concerning the potential of totipotent-stem-like cells and embryo-like structures to fully manifest the trophectoderm/trophoblast lineage's capabilities. Lastly, we briefly discuss the potential applications of induced TSCs in pregnancy-related disease modeling.

Gene trajectory inference for single-cell data by optimal transport metrics

Single-cell RNA sequencing has been widely used to investigate cell state transitions and gene dynamics of biological processes. Current strategies to infer the sequential dynamics of genes in a process typically rely on constructing cell pseudotime through cell trajectory inference. However, the presence of concurrent gene processes in the same group of cells and technical noise can obscure the true progression of the processes studied. To address this challenge, we present GeneTrajectory, an approach that identifies trajectories of genes rather than trajectories of cells. Specifically, optimal transport distances are calculated between gene distributions across the cell-cell graph to extract gene programs and define their gene pseudotemporal order. Here we demonstrate that GeneTrajectory accurately extracts progressive gene dynamics in myeloid lineage maturation. Moreover, we show that GeneTrajectory deconvolves key gene programs underlying mouse skin hair follicle dermal condensate differentiation that could not be resolved by cell trajectory approaches. GeneTrajectory facilitates the discovery of gene programs that control the changes and activities of biological processes.

q-Diffusion leverages the full dimensionality of gene coexpression in single-cell transcriptomics

Unlocking the full dimensionality of single-cell RNA sequencing data (scRNAseq) is the next frontier to a richer, fuller understanding of cell biology. We introduce q-diffusion, a framework for capturing the coexpression structure of an entire library of genes, improving on state-of-the-art analysis tools. The method is demonstrated via three case studies. In the first, q-diffusion helps gain statistical significance for differential effects on patient outcomes when analyzing the CALGB/SWOG 80405 randomized phase III clinical trial, suggesting precision guidance for the treatment of metastatic colorectal cancer. Secondly, q-diffusion is benchmarked against existing scRNAseq classification methods using an in vitro PBMC dataset, in which the proposed method discriminates IFN-γ stimulation more accurately. The same case study demonstrates improvements in unsupervised cell clustering with the recent Tabula Sapiens human atlas. Finally, a local distributional segmentation approach for spatial scRNAseq, driven by q-diffusion, yields interpretable structures of human cortical tissue.

Bioinformatics in urology - molecular characterization of pathophysiology and response to treatment

The application of bioinformatics has revolutionized the practice of medicine in the past 20 years. From early studies that uncovered subtypes of cancer to broad efforts spearheaded by the Cancer Genome Atlas initiative, the use of bioinformatics strategies to analyse high-dimensional data has provided unprecedented insights into the molecular basis of disease. In addition to the identification of disease subtypes - which enables risk stratification - informatics analysis has facilitated the identification of novel risk factors and drivers of disease, biomarkers of progression and treatment response, as well as possibilities for drug repurposing or repositioning; moreover, bioinformatics has guided research towards precision and personalized medicine. Implementation of specific computational approaches such as artificial intelligence, machine learning and molecular subtyping has yet to become widespread in urology clinical practice for reasons of cost, disruption of clinical workflow and need for prospective validation of informatics approaches in independent patient cohorts. Solving these challenges might accelerate routine integration of bioinformatics into clinical settings.

Data- and theory-driven approaches for understanding paths of epithelial-mesenchymal transition

Reversible transitions between epithelial and mesenchymal cell states are a crucial form of epithelial plasticity for development and disease progression. Recent experimental data and mechanistic models showed multiple intermediate epithelial-mesenchymal transition (EMT) states as well as trajectories of EMT underpinned by complex gene regulatory networks. In this review, we summarize recent progress in quantifying EMT and characterizing EMT paths with computational methods and quantitative experiments including omics-level measurements. We provide perspectives on how these studies can help relating fundamental cell biology to physiological and pathological outcomes of EMT.

Two distinct waves of transcriptome and translatome changes drive Drosophila germline stem cell differentiation

The tight control of fate transitions during stem cell differentiation is essential for proper tissue development and maintenance. However, the challenges in studying sparsely distributed adult stem cells in a systematic manner have hindered efforts to identify how the multilayered regulation of gene expression programs orchestrates stem cell differentiation in vivo. Here, we synchronised Drosophila female germline stem cell (GSC) differentiation in vivo to perform in-depth transcriptome and translatome analyses at high temporal resolution. This characterisation revealed widespread and dynamic changes in mRNA level, promoter usage, exon inclusion, and translation efficiency. Transient expression of the master regulator, Bam, drives a first wave of expression changes, primarily modifying the cell cycle program. Surprisingly, as Bam levels recede, differentiating cells return to a remarkably stem cell-like transcription and translation program, with a few crucial changes feeding into a second phase driving terminal differentiation to form the oocyte. Altogether, these findings reveal that rather than a unidirectional accumulation of changes, the in vivo differentiation of stem cells relies on distinctly regulated and developmentally sequential waves.

Partial label learning for automated classification of single-cell transcriptomic profiles

Single-cell RNA sequencing (scRNASeq) data plays a major role in advancing our understanding of developmental biology. An important current question is how to classify transcriptomic profiles obtained from scRNASeq experiments into the various cell types and identify the lineage relationship for individual cells. Because of the fast accumulation of datasets and the high dimensionality of the data, it has become challenging to explore and annotate single-cell transcriptomic profiles by hand. To overcome this challenge, automated classification methods are needed. Classical approaches rely on supervised training datasets. However, due to the difficulty of obtaining data annotated at single-cell resolution, we propose instead to take advantage of partial annotations. The partial label learning framework assumes that we can obtain a set of candidate labels containing the correct one for each data point, a simpler setting than requiring a fully supervised training dataset. We study and extend when needed state-of-the-art multi-class classification methods, such as SVM, kNN, prototype-based, logistic regression and ensemble methods, to the partial label learning framework. Moreover, we study the effect of incorporating the structure of the label set into the methods. We focus particularly on the hierarchical structure of the labels, as commonly observed in developmental processes. We show, on simulated and real datasets, that these extensions enable to learn from partially labeled data, and perform predictions with high accuracy, particularly with a nonlinear prototype-based method. We demonstrate that the performances of our methods trained with partially annotated data reach the same performance as fully supervised data. Finally, we study the level of uncertainty present in the partially annotated data, and derive some prescriptive results on the effect of this uncertainty on the accuracy of the partial label learning methods. Overall our findings show how hierarchical and non-hierarchical partial label learning strategies can help solve the problem of automated classification of single-cell transcriptomic profiles, interestingly these methods rely on a much less stringent type of annotated datasets compared to fully supervised learning methods.

Modeling single cell trajectory using forward-backward stochastic differential equations

Recent advances in single-cell sequencing technology have provided opportunities for mathematical modeling of dynamic developmental processes at the single-cell level, such as inferring developmental trajectories. Optimal transport has emerged as a promising theoretical framework for this task by computing pairings between cells from different time points. However, optimal transport methods have limitations in capturing nonlinear trajectories, as they are static and can only infer linear paths between endpoints. In contrast, stochastic differential equations (SDEs) offer a dynamic and flexible approach that can model non-linear trajectories, including the shape of the path. Nevertheless, existing SDE methods often rely on numerical approximations that can lead to inaccurate inferences, deviating from true trajectories. To address this challenge, we propose a novel approach combining forward-backward stochastic differential equations (FBSDE) with a refined approximation procedure. Our FBSDE model integrates the forward and backward movements of two SDEs in time, aiming to capture the underlying dynamics of single-cell developmental trajectories. Through comprehensive benchmarking on multiple scRNA-seq datasets, we demonstrate the superior performance of FBSDE compared to other methods, highlighting its efficacy in accurately inferring developmental trajectories.

DELVE: feature selection for preserving biological trajectories in single-cell data

Single-cell technologies can measure the expression of thousands of molecular features in individual cells undergoing dynamic biological processes. While examining cells along a computationally-ordered pseudotime trajectory can reveal how changes in gene or protein expression impact cell fate, identifying such dynamic features is challenging due to the inherent noise in single-cell data. Here, we present DELVE, an unsupervised feature selection method for identifying a representative subset of molecular features which robustly recapitulate cellular trajectories. In contrast to previous work, DELVE uses a bottom-up approach to mitigate the effects of confounding sources of variation, and instead models cell states from dynamic gene or protein modules based on core regulatory complexes. Using simulations, single-cell RNA sequencing, and iterative immunofluorescence imaging data in the context of cell cycle and cellular differentiation, we demonstrate how DELVE selects features that better define cell-types and cell-type transitions. DELVE is available as an open-source python package: https://github.com/jranek/delve .

Error modelled gene expression analysis (EMOGEA) provides a superior overview of time course RNA-seq measurements and low count gene expression

Temporal RNA-sequencing (RNA-seq) studies of bulk samples provide an opportunity for improved understanding of gene regulation during dynamic phenomena such as development, tumor progression or response to an incremental dose of a pharmacotherapeutic. Moreover, single-cell RNA-seq (scRNA-seq) data implicitly exhibit temporal characteristics because gene expression values recapitulate dynamic processes such as cellular transitions. Unfortunately, temporal RNA-seq data continue to be analyzed by methods that ignore this ordinal structure and yield results that are often difficult to interpret. Here, we present Error Modelled Gene Expression Analysis (EMOGEA), a framework for analyzing RNA-seq data that incorporates measurement uncertainty, while introducing a special formulation for those acquired to monitor dynamic phenomena. This method is specifically suited for RNA-seq studies in which low-count transcripts with small-fold changes lead to significant biological effects. Such transcripts include genes involved in signaling and non-coding RNAs that inherently exhibit low levels of expression. Using simulation studies, we show that this framework down-weights samples that exhibit extreme responses such as batch effects allowing them to be modeled with the rest of the samples and maintain the degrees of freedom originally envisioned for a study. Using temporal experimental data, we demonstrate the framework by extracting a cascade of gene expression waves from a well-designed RNA-seq study of zebrafish embryogenesis and an scRNA-seq study of mouse pre-implantation and provide unique biological insights into the regulation of genes in each wave. For non-ordinal measurements, we show that EMOGEA has a much higher rate of true positive calls and a vanishingly small rate of false negative discoveries compared to common approaches. Finally, we provide two packages in Python and R that are self-contained and easy to use, including test data.

Label-aware distance mitigates temporal and spatial variability for clustering and visualization of single-cell gene expression data

Clustering and visualization are essential parts of single-cell gene expression data analysis. The Euclidean distance used in most distance-based methods is not optimal. The batch effect, i.e., the variability among samples gathered from different times, tissues, and patients, introduces large between-group distance and obscures the true identities of cells. To solve this problem, we introduce Label-Aware Distance (LAD), a metric using temporal/spatial locality of the batch effect to control for such factors. We validate LAD on simulated data as well as apply it to a mouse retina development dataset and a lung dataset. We also found the utility of our approach in understanding the progression of the Coronavirus Disease 2019 (COVID-19). LAD provides better cell embedding than state-of-the-art batch correction methods on longitudinal datasets. It can be used in distance-based clustering and visualization methods to combine the power of multiple samples to help make biological findings.

DeST-OT: Alignment of Spatiotemporal Transcriptomics Data

Spatially resolved transcriptomics (SRT) measures mRNA transcripts at thousands of locations within a tissue slice, revealing spatial variations in gene expression and distribution of cell types. In recent studies, SRT has been applied to tissue slices from multiple timepoints during the development of an organism. Alignment of this spatiotemporal transcriptomics data can provide insights into the gene expression programs governing the growth and differentiation of cells over space and time. We introduce DeST-OT (Developmental SpatioTemporal Optimal Transport), a method to align SRT slices from pairs of developmental timepoints using the framework of optimal transport (OT). DeST-OT uses semi-relaxed optimal transport to precisely model cellular growth, death, and differentiation processes that are not well-modeled by existing alignment methods. We demonstrate the advantage of DeST-OT on simulated slices. We further introduce two metrics to quantify the plausibility of a spatiotemporal alignment: a growth distortion metric which quantifies the discrepancy between the inferred and the true cell type growth rates, and a migration metric which quantifies the distance traveled between ancestor and descendant cells. DeST-OT outperforms existing methods on these metrics in the alignment of spatiotemporal transcriptomics data from the development of axolotl brain.

Multi-condition and multi-modal temporal profile inference during mouse embryonic development

The emergence of single-cell time-series datasets enables modeling of changes in various types of cellular profiles over time. However, due to the disruptive nature of single-cell measurements, it is impossible to capture the full temporal trajectory of a particular cell. Furthermore, single-cell profiles can be collected at mismatched time points across different conditions (e.g., sex, batch, disease) and data modalities (e.g., scRNA-seq, scATAC-seq), which makes modeling challenging. Here we propose a joint modeling framework, Sunbear, for integrating multi-condition and multi-modal single-cell profiles across time. Sunbear can be used to impute single-cell temporal profile changes, align multi-dataset and multi-modal profiles across time, and extrapolate single-cell profiles in a missing modality. We applied Sunbear to reveal sex-biased transcription during mouse embryonic development and predict dynamic relationships between epigenetic priming and transcription for cells in which multi-modal profiles are unavailable. Sunbear thus enables the projection of single-cell time-series snapshots to multi-modal and multi-condition views of cellular trajectories.

scPerturb: harmonized single-cell perturbation data

Analysis across a growing number of single-cell perturbation datasets is hampered by poor data interoperability. To facilitate development and benchmarking of computational methods, we collect a set of 44 publicly available single-cell perturbation-response datasets with molecular readouts, including transcriptomics, proteomics and epigenomics. We apply uniform quality control pipelines and harmonize feature annotations. The resulting information resource, scPerturb, enables development and testing of computational methods, and facilitates comparison and integration across datasets. We describe energy statistics (E-statistics) for quantification of perturbation effects and significance testing, and demonstrate E-distance as a general distance measure between sets of single-cell expression profiles. We illustrate the application of E-statistics for quantifying similarity and efficacy of perturbations. The perturbation-response datasets and E-statistics computation software are publicly available at scperturb.org. This work provides an information resource for researchers working with single-cell perturbation data and recommendations for experimental design, including optimal cell counts and read depth.

Oscillatory differentiation dynamics fundamentally restricts the resolution of pseudotime reconstruction algorithms

The challenge to understand differentiation and cell lineages in development has resulted in many bioinformatics software tools, notably those working with gene expression data obtained via single-cell RNA sequencing obtained at snapshots in time. Reconstruction methods for trajectories often proceed by dimension reduction, data clustering and then computation of a tree graph in which edges indicate closely related clusters. Cell lineages can then be deduced by following paths through the tree. In the case of multi-potent cells undergoing differentiation, this trajectory reconstruction involves the reconstruction of multiple distinct lineages corresponding to commitment to each of a set of distinct fates. Recent work suggests that there may be cases in which the cell differentiation process involves trajectories that explore, in a dynamic and oscillatory fashion, propensity to differentiate into a number of possible cell fates before commitment finally occurs. Here, we show theoretically that the presence of such oscillations provides intrinsic constraints on the quality and resolution of the trajectory reconstruction process, even for idealized noise-free data. These constraints point to inherent common limitations of current methodologies and serve both to provide additional challenge in the development of software tools and also may help to understand features observed in recent experiments.

Cellular plasticity in reprogramming, rejuvenation and tumorigenesis: a pioneer TF perspective

The multistep process of in vivo reprogramming, mediated by the transcription factors (TFs) Oct4, Sox2, Klf4, and c-Myc (OSKM), holds great promise for the development of rejuvenating and regenerative strategies. However, most of the approaches developed so far are accompanied by a persistent risk of tumorigenicity. Here, we review the groundbreaking effects of in vivo reprogramming with a particular focus on rejuvenation and regeneration. We discuss how the activity of pioneer TFs generates cellular plasticity that may be critical for inducing not only reprogramming and regeneration, but also cancer initiation. Finally, we highlight how a better understanding of the uncoupled control of cellular identity, plasticity, and aging during reprogramming might pave the way to the development of rejuvenating/regenerating strategies in a nontumorigenic manner.

Spatially organized cellular communities form the developing human heart

The heart, which is the first organ to develop, is highly dependent on its form to function1,2. However, how diverse cardiac cell types spatially coordinate to create the complex morphological structures that are crucial for heart function remains unclear. Here we integrated single-cell RNA-sequencing with high-resolution multiplexed error-robust fluorescence in situ hybridization to resolve the identity of the cardiac cell types that develop the human heart. This approach also provided a spatial mapping of individual cells that enables illumination of their organization into cellular communities that form distinct cardiac structures. We discovered that many of these cardiac cell types further specified into subpopulations exclusive to specific communities, which support their specialization according to the cellular ecosystem and anatomical region. In particular, ventricular cardiomyocyte subpopulations displayed an unexpected complex laminar organization across the ventricular wall and formed, with other cell subpopulations, several cellular communities. Interrogating cell-cell interactions within these communities using in vivo conditional genetic mouse models and in vitro human pluripotent stem cell systems revealed multicellular signalling pathways that orchestrate the spatial organization of cardiac cell subpopulations during ventricular wall morphogenesis. These detailed findings into the cellular social interactions and specialization of cardiac cell types constructing and remodelling the human heart offer new insights into structural heart diseases and the engineering of complex multicellular tissues for human heart repair.

Forskolin-driven conversion of human somatic cells into induced neurons through regulation of the cAMP-CREB1-JNK signaling

Human somatic cells can be reprogrammed into neuron cell fate through regulation of a single transcription factor or application of small molecule cocktails. Methods: Here, we report that forskolin efficiently induces the conversion of human somatic cells into induced neurons (FiNs). Results: A large population of neuron-like phenotype cells was observed as early as 24-36 h post-induction. There were >90% TUJ1-, >80% MAP2-, and >80% NEUN-positive neurons at 5 days post-induction. Multiple subtypes of neurons were present among TUJ1-positive cells, including >60% cholinergic, >20% glutamatergic, >10% GABAergic, and >5% dopaminergic neurons. FiNs exhibited typical neural electrophysiological activity in vitro and the ability to survive in vitro and in vivo more than 2 months. Mechanistically, forskolin functions in FiN reprogramming by regulating the cAMP-CREB1-JNK signals, which upregulates cAMP-CREB1 expression and downregulates JNK expression. Conclusion: Overall, our studies identify a safer and efficient single-small-molecule-driven reprogramming approach for induced neuron generation and reveal a novel regulatory mechanism of neuronal cell fate acquisition.

Longitudinal single-cell transcriptional dynamics throughout neurodegeneration in SCA1

Neurodegeneration is a protracted process involving progressive changes in myriad cell types that ultimately results in the death of vulnerable neuronal populations. To dissect how individual cell types within a heterogeneous tissue contribute to the pathogenesis and progression of a neurodegenerative disorder, we performed longitudinal single-nucleus RNA sequencing of mouse and human spinocerebellar ataxia type 1 (SCA1) cerebellar tissue, establishing continuous dynamic trajectories of each cell population. Importantly, we defined the precise transcriptional changes that precede loss of Purkinje cells and, for the first time, identified robust early transcriptional dysregulation in unipolar brush cells and oligodendroglia. Finally, we applied a deep learning method to predict disease state accurately and identified specific features that enable accurate distinction of wild-type and SCA1 cells. Together, this work reveals new roles for diverse cerebellar cell types in SCA1 and provides a generalizable analysis framework for studying neurodegeneration.

Data integration and inference of gene regulation using single-cell temporal multimodal data with scTIE

Single-cell technologies offer unprecedented opportunities to dissect gene regulatory mechanisms in context-specific ways. Although there are computational methods for extracting gene regulatory relationships from scRNA-seq and scATAC-seq data, the data integration problem, essential for accurate cell type identification, has been mostly treated as a standalone challenge. Here we present scTIE, a unified method that integrates temporal multimodal data and infers regulatory relationships predictive of cellular state changes. scTIE uses an autoencoder to embed cells from all time points into a common space by using iterative optimal transport, followed by extracting interpretable information to predict cell trajectories. Using a variety of synthetic and real temporal multimodal data sets, we show scTIE achieves effective data integration while preserving more biological signals than existing methods, particularly in the presence of batch effects and noise. Furthermore, on the exemplar multiome data set we generated from differentiating mouse embryonic stem cells over time, we show scTIE captures regulatory elements highly predictive of cell transition probabilities, providing new potentials to understand the regulatory landscape driving developmental processes.

QOT: Efficient Computation of Sample Level Distance Matrix from Single-Cell Omics Data through Quantized Optimal Transport

Single-cell technologies have emerged as a transformative technology enabling high-dimensional characterization of cell populations at an unprecedented scale. The data's innate complexity and voluminous nature pose significant computational and analytical challenges, especially in comparative studies delineating cellular architectures across various biological conditions (i.e., generation of sample level distance matrices). Optimal Transport (OT) is a mathematical tool that captures the intrinsic structure of data geometrically and has been applied to many bioinformatics tasks. In this paper, we propose QOT (Quantized Optimal Transport), a new method enables efficient computation of sample level distance matrix from large-scale single-cell omics data through a quantization step. We apply our algorithm to real-world single-cell genomics and pathomics datasets, aiming to extrapolate cell-level insights to inform sample level categorizations. Our empirical study shows that QOT outperforms OT-based algorithms in terms of accuracy and robustness when obtaining a distance matrix at the sample level from high throughput single-cell measures. Moreover, the sample level distance matrix could be used in downstream analysis (i.e. uncover the trajectory of disease progression), highlighting its usage in biomedical informatics and data science.

Direct reprogramming of non-limb fibroblasts to cells with properties of limb progenitors

The early limb bud consists of mesenchymal limb progenitors derived from the lateral plate mesoderm (LPM). The LPM also gives rise to the mesodermal components of the flank and neck. However, the cells at these other levels cannot produce the variety of cell types found in the limb. Taking advantage of a direct reprogramming approach, we find a set of factors (Prdm16, Zbtb16, and Lin28a) normally expressed in the early limb bud and capable of imparting limb progenitor-like properties to mouse non-limb fibroblasts. The reprogrammed cells show similar gene expression profiles and can differentiate into similar cell types as endogenous limb progenitors. The further addition of Lin41 potentiates the proliferation of the reprogrammed cells. These results suggest that these same four factors may play pivotal roles in the specification of endogenous limb progenitors.

Molecular mechanisms of cellular dysfunction in testes from men with non-obstructive azoospermia

Male factor infertility affects 50% of infertile couples worldwide; the most severe form, non-obstructive azoospermia (NOA), affects 10-15% of infertile males. Treatment for individuals with NOA is limited to microsurgical sperm extraction paired with in vitro fertilization intracytoplasmic sperm injection. Unfortunately, spermatozoa are only retrieved in ~50% of patients, resulting in live birth rates of 21-46%. Regenerative therapies could provide a solution; however, understanding the cell-type-specific mechanisms of cellular dysfunction is a fundamental necessity to develop precision medicine strategies that could overcome these abnormalities and promote regeneration of spermatogenesis. A number of mechanisms of cellular dysfunction have been elucidated in NOA testicular cells. These mechanisms include abnormalities in both somatic cells and germ cells in NOA testes, such as somatic cell immaturity, aberrant growth factor signalling, increased inflammation, increased apoptosis and abnormal extracellular matrix regulation. Future cell-type-specific investigations in identifying modulators of cellular transcription and translation will be key to understanding upstream dysregulation, and these studies will require development of in vitro models to functionally interrogate spermatogenic niche dysfunction in both somatic and germ cells.

Tracking early mammalian organogenesis - prediction and validation of differentiation trajectories at whole organism scale

Early organogenesis represents a key step in animal development, during which pluripotent cells diversify to initiate organ formation. Here, we sampled 300,000 single-cell transcriptomes from mouse embryos between E8.5 and E9.5 in 6-h intervals and combined this new dataset with our previous atlas (E6.5-E8.5) to produce a densely sampled timecourse of >400,000 cells from early gastrulation to organogenesis. Computational lineage reconstruction identified complex waves of blood and endothelial development, including a new programme for somite-derived endothelium. We also dissected the E7.5 primitive streak into four adjacent regions, performed scRNA-seq and predicted cell fates computationally. Finally, we defined developmental state/fate relationships by combining orthotopic grafting, microscopic analysis and scRNA-seq to transcriptionally determine cell fates of grafted primitive streak regions after 24 h of in vitro embryo culture. Experimentally determined fate outcomes were in good agreement with computationally predicted fates, demonstrating how classical grafting experiments can be revisited to establish high-resolution cell state/fate relationships. Such interdisciplinary approaches will benefit future studies in developmental biology and guide the in vitro production of cells for organ regeneration and repair.

Single-cell multiomics decodes regulatory programs for mouse secondary palate development

Perturbations in gene regulation during palatogenesis can lead to cleft palate, which is among the most common congenital birth defects. Here, we perform single-cell multiome sequencing and profile chromatin accessibility and gene expression simultaneously within the same cells (n = 36,154) isolated from mouse secondary palate across embryonic days (E) 12.5, E13.5, E14.0, and E14.5. We construct five trajectories representing continuous differentiation of cranial neural crest-derived multipotent cells into distinct lineages. By linking open chromatin signals to gene expression changes, we characterize the underlying lineage-determining transcription factors. In silico perturbation analysis identifies transcription factors SHOX2 and MEOX2 as important regulators of the development of the anterior and posterior palate, respectively. In conclusion, our study charts epigenetic and transcriptional dynamics in palatogenesis, serving as a valuable resource for further cleft palate research.

Metabolic control of induced pluripotency

Pluripotent stem cells of the mammalian epiblast and their cultured counterparts-embryonic stem cells (ESCs) and epiblast stem cells (EpiSCs)-have the capacity to differentiate in all cell types of adult organisms. An artificial process of reactivation of the pluripotency program in terminally differentiated cells was established in 2006, which allowed for the generation of induced pluripotent stem cells (iPSCs). This iPSC technology has become an invaluable tool in investigating the molecular mechanisms of human diseases and therapeutic drug development, and it also holds tremendous promise for iPSC applications in regenerative medicine. Since the process of induced reprogramming of differentiated cells to a pluripotent state was discovered, many questions about the molecular mechanisms involved in this process have been clarified. Studies conducted over the past 2 decades have established that metabolic pathways and retrograde mitochondrial signals are involved in the regulation of various aspects of stem cell biology, including differentiation, pluripotency acquisition, and maintenance. During the reprogramming process, cells undergo major transformations, progressing through three distinct stages that are regulated by different signaling pathways, transcription factor networks, and inputs from metabolic pathways. Among the main metabolic features of this process, representing a switch from the dominance of oxidative phosphorylation to aerobic glycolysis and anabolic processes, are many critical stage-specific metabolic signals that control the path of differentiated cells toward a pluripotent state. In this review, we discuss the achievements in the current understanding of the molecular mechanisms of processes controlled by metabolic pathways, and vice versa, during the reprogramming process.

Prediction of single-cell RNA expression profiles in live cells by Raman microscopy with Raman2RNA

Single-cell RNA sequencing and other profiling assays have helped interrogate cells at unprecedented resolution and scale, but are inherently destructive. Raman microscopy reports on the vibrational energy levels of proteins and metabolites in a label-free and nondestructive manner at subcellular spatial resolution, but it lacks genetic and molecular interpretability. Here we present Raman2RNA (R2R), a method to infer single-cell expression profiles in live cells through label-free hyperspectral Raman microscopy images and domain translation. We predict single-cell RNA sequencing profiles nondestructively from Raman images using either anchor-based integration with single molecule fluorescence in situ hybridization, or anchor-free generation with adversarial autoencoders. R2R outperformed inference from brightfield images (cosine similarities: R2R >0.85 and brightfield <0.15). In reprogramming of mouse fibroblasts into induced pluripotent stem cells, R2R inferred the expression profiles of various cell states. With live-cell tracking of mouse embryonic stem cell differentiation, R2R traced the early emergence of lineage divergence and differentiation trajectories, overcoming discontinuities in expression space. R2R lays a foundation for future exploration of live genomic dynamics.

Transcription Factors in Brain Regeneration: A Potential Novel Therapeutic Target

Transcription factors play a crucial role in providing identity to each cell population. To maintain cell identity, it is essential to balance the expression of activator and inhibitor transcription factors. Cell plasticity and reprogramming offer great potential for future therapeutic applications, as they can regenerate damaged tissue. Specific niche factors can modify gene expression and differentiate or transdifferentiate the target cell to the required fate. Ongoing research is being carried out on the possibilities of transcription factors in regenerating neurons, with neural stem cells (NSCs) being considered the preferred cells for generating new neurons due to their epigenomic and transcriptome memory. NEUROD1/ASCL1, BRN2, MYTL1, and other transcription factors can induce direct reprogramming of somatic cells, such as fibroblasts, into neurons. However, the molecular biology of transcription factors in reprogramming and differentiation still needs to be fully understood.

Vision sculpts a continuum of L2/3 cell types in the visual cortex during the critical period

We previously reported that vision specifies Layer 2/3 (L2/3) glutamatergic cell-type identity in the primary visual cortex (V1). Using unsupervised clustering of single-nucleus RNA-sequencing data, we identified molecularly distinct L2/3 cell types in normal-reared (NR) and dark-reared (DR) mice, but the two sets exhibited poor correspondence. Here, we show that classification of cell types was confounded in DR by vision-dependent gene programs that are orthogonal to gene programs underlying cell-type identity. A focused clustering analysis successfully matches cell types between DR and NR, suggesting that cell identity-defining gene programs persist under vision deprivation but are overshadowed by vision-dependent transcriptomic variation. Using multi-tasking theory we show that L2/3 cell types form a continuum between three cell-archetypes. Visual deprivation markedly shifts this distribution along the continuum. Thus, dark-rearing markedly influences cell states thereby masking cell-type-identities and changes the distribution of L2/3 types along a transcriptomic continuum.

Murine trophoblast organoids as a model for trophoblast development and CRISPR-Cas9 screening

The placenta becomes one of the most diversified organs during placental mammal radiation. The main in vitro model for studying mouse trophoblast development is the 2D differentiation model of trophoblast stem cells, which is highly skewed to certain lineages and thus hampers systematic screens. Here, we established culture conditions for the establishment, maintenance, and differentiation of murine trophoblast organoids. Murine trophoblast organoids under the maintenance condition contain stem cell-like populations, whereas differentiated organoids possess various trophoblasts resembling placental ones in vivo. Ablation of Nubpl or Gcm1 in trophoblast organoids recapitulated their deficiency phenotypes in vivo, suggesting that those organoids are valid in vitro models for trophoblast development. Importantly, we performed an efficient CRISPR-Cas9 screening in mouse trophoblast organoids using a focused sgRNA (single guide RNA) library targeting G protein-coupled receptors. Together, our results establish an organoid model to investigate mouse trophoblast development and a practicable approach to performing forward screening in trophoblast lineages.

Deciphering the decisive factors driving fate bifurcations in somatic cell reprogramming

Single-cell studies have demonstrated that somatic cell reprogramming is a continuous process of cell fates transition. Only partial reprogramming intermediates can overcome the molecular bottlenecks to acquire pluripotency. To decipher the underlying decisive factors driving cell fate, we identified induced pluripotent stem cells or stromal-like cells (iPSCs/SLCs) and iPSCs or trophoblast-like cells (iPSCs/TLCs) fate bifurcations by reconstructing cellular trajectory. The mesenchymal-epithelial transition and the activation of pluripotency networks are the main molecular series in successful reprogramming. Correspondingly, intermediates diverge into SLCs accompanied by the inhibition of cell cycle genes and the activation of extracellular matrix genes, whereas the TLCs fate is characterized by the up-regulation of placenta development genes. Combining putative gene regulatory networks, seven (Taf7, Ezh2, Klf2, etc.) and three key factors (Cdc5l, Klf4, and Nanog) were individually identified as drivers of the successful reprogramming by triggering downstream pluripotent networks during iPSCs/SLCs and iPSCs/TLCs fate bifurcation. Conversely, 11 factors (Cebpb, Sox4, Junb, etc.) and four factors (Gata2, Jund, Ctnnb1, etc.) drive SLCs fate and TLCs fate, respectively. Our study sheds new light on the understanding of decisive factors driving cell fate, which is helpful for improving reprogramming efficiency through manipulating cell fates to avoid alternative fates.

Synthetic genetic circuits to uncover the OCT4 trajectories of successful reprogramming of human fibroblasts

Reprogramming human fibroblasts to induced pluripotent stem cells (iPSCs) is inefficient, with heterogeneity among transcription factor (TF) trajectories driving divergent cell states. Nevertheless, the impact of TF dynamics on reprogramming efficiency remains uncharted. We develop a system that accurately reports OCT4 protein levels in live cells and use it to reveal the trajectories of OCT4 in successful reprogramming. Our system comprises a synthetic genetic circuit that leverages noise to generate a wide range of OCT4 trajectories and a microRNA targeting endogenous OCT4 to set total cellular OCT4 protein levels. By fusing OCT4 to a fluorescent protein, we are able to track OCT4 trajectories with clonal resolution via live-cell imaging. We discover that a supraphysiological, stable OCT4 level is required, but not sufficient, for efficient iPSC colony formation. Our synthetic genetic circuit design and high-throughput live-imaging pipeline are generalizable for investigating TF dynamics for other cell fate programming applications.

Elucidating dynamic cell lineages and gene networks in time-course single cell differentiation

Single cell RNA sequencing (scRNA-seq) technologies provide researchers with an unprecedented opportunity to exploit cell heterogeneity. For example, the sequenced cells belong to various cell lineages, which may have different cell fates in stem and progenitor cells. Those cells may differentiate into various mature cell types in a cell differentiation process. To trace the behavior of cell differentiation, researchers reconstruct cell lineages and predict cell fates by ordering cells chronologically into a trajectory with a pseudo-time. However, in scRNA-seq experiments, there are no cell-to-cell correspondences along with the time to reconstruct the cell lineages, which creates a significant challenge for cell lineage tracing and cell fate prediction. Therefore, methods that can accurately reconstruct the dynamic cell lineages and predict cell fates are highly desirable. In this article, we develop an innovative machine-learning framework called Cell Smoothing Transformation (CellST) to elucidate the dynamic cell fate paths and construct gene networks in cell differentiation processes. Unlike the existing methods that construct one single bulk cell trajectory, CellST builds cell trajectories and tracks behaviors for each individual cell. Additionally, CellST can predict cell fates even for less frequent cell types. Based on the individual cell fate trajectories, CellST can further construct dynamic gene networks to model gene-gene relationships along the cell differentiation process and discover critical genes that potentially regulate cells into various mature cell types.

Reconstructing growth and dynamic trajectories from single-cell transcriptomics data

Time-series single-cell RNA sequencing (scRNA-seq) datasets provide unprecedented opportunities to learn dynamic processes of cellular systems. Due to the destructive nature of sequencing, it remains challenging to link the scRNA-seq snapshots sampled at different time points. Here we present TIGON, a dynamic, unbalanced optimal transport algorithm that reconstructs dynamic trajectories and population growth simultaneously as well as the underlying gene regulatory network from multiple snapshots. To tackle the high-dimensional optimal transport problem, we introduce a deep learning method using a dimensionless formulation based on the Wasserstein-Fisher-Rao (WFR) distance. TIGON is evaluated on simulated data and compared with existing methods for its robustness and accuracy in predicting cell state transition and cell population growth. Using three scRNA-seq datasets, we show the importance of growth in the temporal inference, TIGON's capability in reconstructing gene expression at unmeasured time points and its applications to temporal gene regulatory networks and cell-cell communication inference.

Proliferation history and transcription factor levels drive direct conversion

The sparse and stochastic nature of reprogramming has obscured our understanding of how transcription factors drive cells to new identities. To overcome this limit, we developed a compact, portable reprogramming system that increases direct conversion of fibroblasts to motor neurons by two orders of magnitude. We show that subpopulations with different reprogramming potentials are distinguishable by proliferation history. By controlling for proliferation history and titrating each transcription factor, we find that conversion correlates with levels of the pioneer transcription factor Ngn2, whereas conversion shows a biphasic response to Lhx3. Increasing the proliferation rate of adult human fibroblasts generates morphologically mature, induced motor neurons at high rates. Using compact, optimized, polycistronic cassettes, we generate motor neurons that graft with the murine central nervous system, demonstrating the potential for in vivo therapies.

Reconstructing human brown fat developmental trajectory in vitro

Brown adipocytes (BAs) represent a specialized cell type that is able to uncouple nutrient catabolism from ATP generation to dissipate energy as heat. In humans, the brown fat tissue is composed of discrete depots found throughout the neck and trunk region. BAs originate from a precursor common to skeletal muscle, but their developmental trajectory remains poorly understood. Here, we used single-cell RNA sequencing to characterize the development of interscapular brown fat in mice. Our analysis identified a transient stage of BA differentiation characterized by the expression of the transcription factor GATA6. We show that recapitulating the sequence of signaling cues identified in mice can lead to efficient differentiation of BAs in vitro from human pluripotent stem cells. These precursors can in turn be efficiently converted into functional BAs that can respond to signals mimicking adrenergic stimuli by increasing their metabolism, resulting in heat production.

scTOP: physics-inspired order parameters for cellular identification and visualization

Advances in single-cell RNA sequencing provide an unprecedented window into cellular identity. The abundance of data requires new theoretical and computational frameworks to analyze the dynamics of differentiation and integrate knowledge from cell atlases. We present 'single-cell Type Order Parameters' (scTOP): a statistical, physics-inspired approach for quantifying cell identity given a reference basis of cell types. scTOP can accurately classify cells, visualize developmental trajectories and assess the fidelity of engineered cells. Importantly, scTOP does this without feature selection, statistical fitting or dimensional reduction (e.g. uniform manifold approximation and projection, principle components analysis, etc.). We illustrate the power of scTOP using human and mouse datasets. By reanalyzing mouse lung data, we characterize a transient hybrid alveolar type 1/alveolar type 2 cell population. Visualizations of lineage tracing hematopoiesis data using scTOP confirm that a single clone can give rise to multiple mature cell types. We assess the transcriptional similarity between endogenous and donor-derived cells in the context of murine pulmonary cell transplantation. Our results suggest that physics-inspired order parameters can be an important tool for understanding differentiation and characterizing engineered cells. scTOP is available as an easy-to-use Python package.

IL-1β+ macrophages fuel pathogenic inflammation in pancreatic cancer

Pancreatic ductal adenocarcinoma (PDAC) is a lethal disease with high resistance to therapies1. Inflammatory and immunomodulatory signals co-exist in the pancreatic tumour microenvironment, leading to dysregulated repair and cytotoxic responses. Tumour-associated macrophages (TAMs) have key roles in PDAC2, but their diversity has prevented therapeutic exploitation. Here we combined single-cell and spatial genomics with functional experiments to unravel macrophage functions in pancreatic cancer. We uncovered an inflammatory loop between tumour cells and interleukin-1β (IL-1β)-expressing TAMs, a subset of macrophages elicited by a local synergy between prostaglandin E2 (PGE2) and tumour necrosis factor (TNF). Physical proximity with IL-1β+ TAMs was associated with inflammatory reprogramming and acquisition of pathogenic properties by a subset of PDAC cells. This occurrence was an early event in pancreatic tumorigenesis and led to persistent transcriptional changes associated with disease progression and poor outcomes for patients. Blocking PGE2 or IL-1β activity elicited TAM reprogramming and antagonized tumour cell-intrinsic and -extrinsic inflammation, leading to PDAC control in vivo. Targeting the PGE2-IL-1β axis may enable preventive or therapeutic strategies for reprogramming of immune dynamics in pancreatic cancer.