A relay velocity model infers cell-dependent RNA velocity

Adult Hippocampal Neurogenesis in the Human Brain: Updates, Challenges, and Perspectives

The existence of neurogenesis in the adult human hippocampus has been under considerable debate within the past three decades due to the diverging conclusions originating mostly from immunohistochemistry studies. While some of these reports conclude that hippocampal neurogenesis in humans occurs throughout physiologic aging, others indicate that this phenomenon ends by early childhood. More recently, some groups have adopted next-generation sequencing technologies to characterize with more acuity the extent of this phenomenon in humans. Here, we review the current state of research on adult hippocampal neurogenesis in the human brain with an emphasis on the challenges and limitations of using immunohistochemistry and next-generation sequencing technologies for its study.

DelaySSA: stochastic simulation of biochemical systems and gene regulatory networks with or without time delays

Stochastic Simulation Algorithm (SSA) is crucial for modeling biochemical reactions and gene regulatory networks. Traditional SSA is characterized by Markovian property and cannot naturally model systems with time delays. Several algorithms have already been designed to handle delayed reactions, yet few easy-to-use implementations exist. To address these challenges, we have developed DelaySSA, an R package that implements currently available algorithms for SSA with or without delays. Meanwhile, we also provided Matlab and Python versions to support wider applications. We demonstrated its accuracy and validity by simulating two classical models: the Bursty model and Refractory model. We then tested its capability to simulate the RNA Velocity model, where it successfully reproduced both the up- and down-regulation stages in the phase portrait. Finally, we extended its application to simulate a gene regulatory network of lung cancer adeno-to-squamous transition (AST) and qualitatively analyzed its bistability behavior by approximating the Waddington's landscape. Modeling the therapeutic intervention of a SOX2 degrader as a delayed degradation reaction, AST is effectively blocked and reprogrammed back to the adenocarcinoma state, providing a useful clue for targeting drug-resistant AST in the future. Taken together, DelaySSA is a powerful and easy-to-use software suite, facilitating accurate modeling of various kinds of biological systems and broadening the scope of stochastic simulations in systems biology.

Cell2fate infers RNA velocity modules to improve cell fate prediction

RNA velocity exploits the temporal information contained in spliced and unspliced RNA counts to infer transcriptional dynamics. Existing velocity models often rely on coarse biophysical simplifications or numerical approximations to solve the underlying ordinary differential equations (ODEs), which can compromise accuracy in challenging settings, such as complex or weak transcription rate changes across cellular trajectories. Here we present cell2fate, a formulation of RNA velocity based on a linearization of the velocity ODE, which allows solving a biophysically more accurate model in a fully Bayesian fashion. As a result, cell2fate decomposes the RNA velocity solutions into modules, providing a biophysical connection between RNA velocity and statistical dimensionality reduction. We comprehensively benchmark cell2fate in real-world settings, demonstrating enhanced interpretability and power to reconstruct complex dynamics and weak dynamical signals in rare and mature cell types. Finally, we apply cell2fate to the developing human brain, where we spatially map RNA velocity modules onto the tissue architecture, connecting the spatial organization of tissues with temporal dynamics of transcription.

Intestinal secretory differentiation reflects niche-driven phenotypic and epigenetic plasticity of a common signal-responsive terminal cell

Enterocytes and four classic secretory cell types derive from intestinal epithelial stem cells. Based on morphology, location, and canonical markers, goblet and Paneth cells are considered distinct secretory types. Here, we report high overlap in their transcripts and sites of accessible chromatin, in marked contrast to those of their enteroendocrine or tuft cell siblings. Mouse and human goblet and Paneth cells express extraordinary fractions of few antimicrobial genes, which reflect specific responses to local niches. Wnt signaling retains some ATOH1+ secretory cells in crypt bottoms, where the absence of BMP signaling potently induces Paneth features. Cells that migrate away from crypt bottoms encounter BMPs and thereby acquire goblet properties. These phenotypes and underlying accessible cis-elements interconvert in post-mitotic cells. Thus, goblet and Paneth properties represent alternative phenotypic manifestations of a common signal-responsive terminal cell type. These findings reveal exquisite niche-dependent cell plasticity and cis-regulatory dynamics in likely response to antimicrobial needs.

How to view the female reproductive tract through single-cell looking glasses

Single-cell technologies have emerged as an unprecedented tool for biologists and clinicians, allowing them to assess organs and tissues at the level of individual cells. In the field of women's reproductive biology, single-cell studies have provided insights into the cellular and molecular processes that regulate reproductive and obstetrical functions in health and disease. The knowledge that these studies generate is helping clinicians to improve the understanding and diagnosis of infertility related issues or pregnancy complications and to find new avenues for their treatment. However, navigating the expansive landscape of this type of transcriptomic data analysis represents a pivotal challenge in current research. Single cell RNA sequencing involves isolating cells into droplets, reverse transcribing RNA to generate complementary DNA, with each droplet content uniquely labeled by a barcode. Upon sequencing the complementary DNAs, the barcodes enable the reassignment of sequencing reads to individual droplets, facilitating the reconstruction of the cellular landscape of the sample obtained from a tissue or organ and beyond. Researchers, equipped with the metaphorical 'single-cell glasses,' must adequately choose from a plethora of strategies to dissect and interpret cellular information. Sophisticated algorithms and the decision-making process are often underestimated, resulting in artefactual or cumbersome interpreted results. Computational biologists apply and innovate computational tools designed to process, model, and interpret expansive datasets. The ramifications of their work extend far beyond the realm of data processing; they give shape to the outcome of analyses, playing a pivotal role in drawing meaningful conclusions from the wealth of information garnered. In this review, we describe the wide variety of approaches and analytical steps available with enough detail to gain a concise picture of what a complete examination of a single-cell dataset would be. We commence with a discussion on key points in experimental design, highlighting crucial questions one should consider. Following this, we delve into the various preprocessing and quality control steps essential for any single-cell dataset. The subsequent section offers a detailed guide on constructing a single-cell atlas, exploring nuances such as differential characteristics in visualization and clustering techniques, as well as strategies for assigning identity to cell populations through gene marker annotations. Moving beyond the creation of an atlas, we explore methods for investigating pathological conditions. This involves conducting cell population comparison tests between conditions and analyzing specific cell-to-cell communications and cellular differentiation trajectories in both health and disease scenarios. This work aims to furnish a newcomer researcher and/or clinician with essential guidelines to embark on a single-cell adventure without succumbing to common pitfalls. By bridging the gap between theory and practice, it facilitates the translation of single-cell technologies into clinically relevant applications. Throughout the manuscript, practical examples of its usage in women's reproductive health studies are provided. Various sections delve into specific clinical scenarios, demonstrating how these guidelines can be instrumental in unraveling the molecular landscapes of diseases and physiological processes related to women's reproduction.

In-depth functional analysis of BRD9 in fetal hematopoiesis reveals context-dependent roles

The hierarchical organization of hematopoietic stem cells (HSCs) governing adult hematopoiesis has been extensively investigated. However, the dynamic epigenomic transition from fetal to adult hematopoiesis remains incompletely understood, particularly regarding the involvement of epigenetic factors. In this study, we investigate the roles of BRD9, an essential component of the non-canonical BAF (ncBAF) complex known to govern the fate of adult HSCs, in fetal hematopoiesis. Consistent with observations in adult hematopoiesis, BRD9 loss impairs fetal HSC stemness and disturbs erythroid maturation. Intriguingly, the impact on myeloid lineage was discrepant: BRD9 loss inhibited and promoted myeloid differentiation in fetal and adult models, respectively. Through comprehensive transcriptomic and epigenomic analysis, we elucidate the differential roles of BRD9 in a context- and lineage-dependent manner. Our data uncover how BRD9/ncBAF complex modulates transcription in a stage-specific manner, providing deeper insights into the epigenetic regulation underlying the transition from fetal to adult hematopoiesis.

The landscape of cell lineage tracing

Cell fate changes play a crucial role in the processes of natural development, disease progression, and the efficacy of therapeutic interventions. The definition of the various types of cell fate changes, including cell expansion, differentiation, transdifferentiation, dedifferentiation, reprogramming, and state transitions, represents a complex and evolving field of research known as cell lineage tracing. This review will systematically introduce the research history and progress in this field, which can be broadly divided into two parts: prospective tracing and retrospective tracing. The initial section encompasses an array of methodologies pertaining to isotope labeling, transient fluorescent tracers, non-fluorescent transient tracers, non-fluorescent genetic markers, fluorescent protein, genetic marker delivery, genetic recombination, exogenous DNA barcodes, CRISPR-Cas9 mediated DNA barcodes, and base editor-mediated DNA barcodes. The second part of the review covers genetic mosaicism, genomic DNA alteration, TCR/BCR, DNA methylation, and mitochondrial DNA mutation. In the final section, we will address the principal challenges and prospective avenues of enquiry in the field of cell lineage tracing, with a particular focus on the sequencing techniques and mathematical models pertinent to single-cell genetic lineage tracing, and the value of pursuing a more comprehensive investigation at both the spatial and temporal levels in the study of cell lineage tracing.

Time-resolved single-cell secretion analysis via microfluidics

Revealing how individual cells alter their secretions over time is crucial for understanding their responses to environmental changes. Key questions include: When do cells modify their functions and states? What transitions occur? Insights into the kinetic secretion trajectories of various cell types are essential for unraveling complex biological systems. This review highlights seven microfluidic technologies for time-resolved single-cell secretion analysis: 1. Microwell real-time electrical detection: uses microelectrodes for precise, cell-specific, real-time measurement of secreted molecules. 2. Microwell real-time optical detection: employs advanced optical systems for real-time, multiplexed monitoring of cellular secretions. 3. Microvalve real-time optical detection: dynamically analyzes secretions under controlled in situ stimuli, enabling detailed kinetic studies at the single-cell level. 4. Droplet real-time optical detection: provides superior throughput by generating droplets containing single cells and sensors for high-throughput screening. 5. Microwell time-barcoded optical detection: utilizes sequential barcoding techniques to facilitate scalable assays for tracking multiple secretions over time. 6. Microvalve time-barcoded optical detection: incorporates automated time-barcoding via micro-valves for robust and scalable analysis. 7. Microwell time-barcoded sequencing: captures and labels secretions for sequencing, enabling multidimensional analysis, though currently limited to a few time points and extended intervals. This review specifically addresses the challenges of achieving high-resolution timing measurements with short intervals while maintaining scalability for single-cell screening. Future advancements in microfluidic devices, integrating innovative barcoding technologies, advanced imaging technologies, artificial intelligence-powered decoding and analysis, and automations are anticipated to enable highly sensitive, scalable, high-throughput single-cell dynamic analysis. These developments hold great promise for deepening our understanding of biosystems by exploring single-cell timing responses on a larger scale.

Identification of Meibomian gland stem cell populations and mechanisms of aging

Meibomian glands secrete lipid-rich meibum, which prevents tear evaporation. Aging-related Meibomian gland shrinkage may result in part from stem cell exhaustion and is associated with evaporative dry eye disease, a common condition lacking effective treatment. The identities and niche of Meibomian gland stem cells and the signals controlling their activity are poorly defined. Using snRNA-seq, in vivo lineage tracing, ex vivo live imaging, and genetic studies in mice, we identify markers for stem cell populations that maintain distinct regions of the gland and uncover Hedgehog (Hh) signaling as a key regulator of stem cell proliferation. Consistent with this, we show that human Meibomian gland carcinoma exhibits increased Hh signaling. Aged glands display decreased Hh and EGF signaling, deficient innervation, and loss of collagen I in niche fibroblasts, indicating that alterations in both glandular epithelial cells and their surrounding microenvironment contribute to age-related degeneration. These findings suggest new approaches to treat aging-associated Meibomian gland loss.

Towards understanding cancer dormancy over strategic hitching up mechanisms to technologies

Delving into cancer dormancy has been an inherent task that may drive the lethal recurrence of cancer after primary tumor relief. Cells in quiescence can survive for a short or long term in silence, may undergo genetic or epigenetic changes, and can initiate relapse through certain contextual cues. The state of dormancy can be induced by multiple conditions including cancer drug treatment, in turn, undergoes a life cycle that generally occurs through dissemination, invasion, intravasation, circulation, immune evasion, extravasation, and colonization. Throughout this cascade, a cellular machinery governs the fate of individual cells, largely affected by gene regulation. Despite its significance, a precise view of cancer dormancy is yet hampered. Revolutionizing advanced single cell and long read sequencing through analysis methodologies and artificial intelligence, the most recent stage in the research tool progress, is expected to provide a holistic view of the diverse aspects of cancer dormancy.

Multi-omics-based mapping of decidualization resistance in patients with a history of severe preeclampsia

Endometrial decidualization resistance (DR) is implicated in various gynecological and obstetric conditions. Here, using a multi-omic strategy, we unraveled the cellular and molecular characteristics of DR in patients who have suffered severe preeclampsia (sPE). Morphological analysis unveiled significant glandular anatomical abnormalities, confirmed histologically and quantified by the digitization of hematoxylin and eosin-stained tissue sections. Single-cell RNA sequencing (scRNA-seq) of endometrial samples from patients with sPE (n = 11) and controls (n = 12) revealed sPE-associated shifts in cell composition, manifesting as a stromal mosaic state characterized by proliferative stromal cells (MMP11 and SFRP4) alongside IGFBP1+ decidualized cells, with concurrent epithelial mosaicism and a dearth of epithelial-stromal transition associated with decidualization. Cell-cell communication network mapping underscored aberrant crosstalk among specific cell types, implicating crucial pathways such as endoglin, WNT and SPP1. Spatial transcriptomics in a replication cohort validated DR-associated features. Laser capture microdissection/mass spectrometry in a second replication cohort corroborated several scRNA-seq findings, notably the absence of stromal to epithelial transition at a pathway level, indicating a disrupted response to steroid hormones, particularly estrogens. These insights shed light on potential molecular mechanisms underpinning DR pathogenesis in the context of sPE.

GraphVelo allows for accurate inference of multimodal omics velocities and molecular mechanisms for single cells

RNA velocities and generalizations emerge as powerful approaches for extracting time-resolved information from high-throughput snapshot single-cell data. Yet, several inherent limitations restrict applying the approaches to genes not suitable for RNA velocity inference due to complex transcriptional dynamics, low expression, or lacking splicing dynamics, or data of non-transcriptomic modality. Here, we present GraphVelo, a graph-based machine learning procedure that uses as input the RNA velocities inferred from existing methods and infers velocity vectors lying in the tangent space of the low-dimensional manifold formed by the single cell data. GraphVelo preserves vector magnitude and direction information during transformations across different data representations. Tests on multiple synthetic and experimental scRNA-seq data including viral-host interactome and multi-omics datasets demonstrate that GraphVelo, together with downstream generalized dynamo analyses, extends RNA velocities to multi-modal data and reveals quantitative nonlinear regulation relations between genes, virus and host cells, and different layers of gene regulation.

GraphVelo allows for accurate inference of multimodal velocities and molecular mechanisms for single cells

RNA velocities and generalizations emerge as powerful approaches for extracting time-resolved information from high-throughput snapshot single-cell data. Yet, several inherent limitations restrict applying the approaches to genes not suitable for RNA velocity inference due to complex transcriptional dynamics, low expression, or lacking splicing dynamics, or data of non-transcriptomic modality. Here, we present GraphVelo, a graph-based machine learning procedure that uses as input the RNA velocities inferred from existing methods and infers velocity vectors lying in the tangent space of the low-dimensional manifold formed by the single cell data. GraphVelo preserves vector magnitude and direction information during transformations across different data representations. Tests on multiple synthetic and experimental scRNA-seq data including viral-host interactome and multi-omics datasets demonstrate that GraphVelo, together with downstream generalized dynamo analyses, extends RNA velocities to multi-modal data and reveals quantitative nonlinear regulation relations between genes, virus and host cells, and different layers of gene regulation.

Control of Cellular Differentiation Trajectories for Cancer Reversion

Cellular differentiation is controlled by intricate layers of gene regulation, involving the modulation of gene expression by various transcriptional regulators. Due to the complexity of gene regulation, identifying master regulators across the differentiation trajectory has been a longstanding challenge. To tackle this problem, a computational framework, single-cell Boolean network inference and control (BENEIN), is presented. Applying BENEIN to human large intestinal single-cell transcriptome data, MYB, HDAC2, and FOXA2 are identified as the master regulators whose inhibition induces enterocyte differentiation. It is found that simultaneous knockdown of these master regulators can revert colorectal cancer cells into normal-like enterocytes by synergistically inducing differentiation and suppressing malignancy, which is validated by in vitro and in vivo experiments.

Linking transcriptome and morphology in bone cells at cellular resolution with generative AI

Recent advancements in deep learning (DL) have revolutionized the capability of artificial intelligence (AI) by enabling the analysis of large-scale, complex datasets that are difficult for humans to interpret. However, large amounts of high-quality data are required to train such generative AI models successfully. With the rapid commercialization of single-cell sequencing and spatial transcriptomics platforms, the field is increasingly producing large-scale datasets such as histological images, single-cell molecular data, and spatial transcriptomic data. These molecular and morphological datasets parallel the multimodal text and image data used to train highly successful generative AI models for natural language processing and computer vision. Thus, these emerging data types offer great potential to train generative AI models that uncover intricate biological processes of bone cells at a cellular level. In this Perspective, we summarize the progress and prospects of generative AI applied to these datasets and their potential applications to bone research. In particular, we highlight three AI applications: predicting cell differentiation dynamics, linking molecular and morphological features, and predicting cellular responses to perturbations. To make generative AI models beneficial for bone research, important issues, such as technical biases in bone single-cell datasets, lack of profiling of important bone cell types, and lack of spatial information, needs to be addressed. Realizing the potential of generative AI for bone biology will also likely require generating large-scale, high-quality cellular-resolution spatial transcriptomics datasets, improving the sensitivity of current spatial transcriptomics datasets, and thorough experimental validation of model predictions.

Multivariate stochastic modeling for transcriptional dynamics with cell-specific latent time using SDEvelo

Recently, RNA velocity has driven a paradigmatic change in single-cell RNA sequencing (scRNA-seq) studies, allowing the reconstruction and prediction of directed trajectories in cell differentiation and state transitions. Most existing methods of dynamic modeling use ordinary differential equations (ODE) for individual genes without applying multivariate approaches. However, this modeling strategy inadequately captures the intrinsically stochastic nature of transcriptional dynamics governed by a cell-specific latent time across multiple genes, potentially leading to erroneous results. Here, we present SDEvelo, a generative approach to inferring RNA velocity by modeling the dynamics of unspliced and spliced RNAs via multivariate stochastic differential equations (SDE). Uniquely, SDEvelo explicitly models inherent uncertainty in transcriptional dynamics while estimating a cell-specific latent time across genes. Using both simulated and four scRNA-seq and spatial transcriptomics datasets, we show that SDEvelo can model the random dynamic patterns of mature-state cells while accurately detecting carcinogenesis. Additionally, the estimated gene-shared latent time can facilitate many downstream analyses for biological discovery. We demonstrate that SDEvelo is computationally scalable and applicable to both scRNA-seq and sequencing-based spatial transcriptomics data.

A statistical approach for systematic identification of transition cells from scRNA-seq data

Decoding cellular state transitions is crucial for understanding complex biological processes in development and disease. While recent advancements in single-cell RNA sequencing (scRNA-seq) offer insights into cellular trajectories, existing tools primarily study expressional rather than regulatory state shifts. We present CellTran, a statistical approach utilizing paired-gene expression correlations to detect transition cells from scRNA-seq data without explicitly resolving gene regulatory networks. Applying our approach to various contexts, including tissue regeneration, embryonic development, preinvasive lesions, and humoral responses post-vaccination, reveals transition cells and their distinct gene expression profiles. Our study sheds light on the underlying molecular mechanisms driving cellular state transitions, enhancing our ability to identify therapeutic targets for disease interventions.

Unlocking biological insights from differentially expressed genes: Concepts, methods, and future perspectives

BACKGROUND

Identifying differentially expressed genes (DEGs) is a core task of transcriptome analysis, as DEGs can reveal the molecular mechanisms underlying biological processes. However, interpreting the biological significance of large DEG lists is challenging. Currently, gene ontology, pathway enrichment and protein-protein interaction analysis are common strategies employed by biologists. Additionally, emerging analytical strategies/approaches (such as network module analysis, knowledge graph, drug repurposing, cell marker discovery, trajectory analysis, and cell communication analysis) have been proposed. Despite these advances, comprehensive guidelines for systematically and thoroughly mining the biological information within DEGs remain lacking.

AIM OF REVIEW

This review aims to provide an overview of essential concepts and methodologies for the biological interpretation of DEGs, enhancing the contextual understanding. It also addresses the current limitations and future perspectives of these approaches, highlighting their broad applications in deciphering the molecular mechanism of complex diseases and phenotypes. To assist users in extracting insights from extensive datasets, especially various DEG lists, we developed DEGMiner (https://www.ciblab.net/DEGMiner/), which integrates over 300 easily accessible databases and tools.

KEY SCIENTIFIC CONCEPTS OF REVIEW

This review offers strong support and guidance for exploring DEGs, and also will accelerate the discovery of hidden biological insights within genomes.

Statistical inference with a manifold-constrained RNA velocity model uncovers cell cycle speed modulations

Across biological systems, cells undergo coordinated changes in gene expression, resulting in transcriptome dynamics that unfold within a low-dimensional manifold. While low-dimensional dynamics can be extracted using RNA velocity, these algorithms can be fragile and rely on heuristics lacking statistical control. Moreover, the estimated vector field is not dynamically consistent with the traversed gene expression manifold. To address these challenges, we introduce a Bayesian model of RNA velocity that couples velocity field and manifold estimation in a reformulated, unified framework, identifying the parameters of an explicit dynamical system. Focusing on the cell cycle, we implement VeloCycle to study gene regulation dynamics on one-dimensional periodic manifolds and validate its ability to infer cell cycle periods using live imaging. We also apply VeloCycle to reveal speed differences in regionally defined progenitors and Perturb-seq gene knockdowns. Overall, VeloCycle expands the single-cell RNA sequencing analysis toolkit with a modular and statistically consistent RNA velocity inference framework.

scGRN-Entropy: Inferring cell differentiation trajectories using single-cell data and gene regulation network-based transfer entropy

Research on cell differentiation facilitates a deeper understanding of the fundamental processes of life, elucidates the intrinsic mechanisms underlying diseases such as cancer, and advances the development of therapeutics and precision medicine. Existing methods for inferring cell differentiation trajectories from single-cell RNA sequencing (scRNA-seq) data primarily rely on static gene expression data to measure distances between cells and subsequently infer pseudotime trajectories. In this work, we introduce a novel method, scGRN-Entropy, for inferring cell differentiation trajectories and pseudotime from scRNA-seq data. Unlike existing approaches, scGRN-Entropy improves inference accuracy by incorporating dynamic changes in gene regulatory networks (GRN). In scGRN-Entropy, an undirected graph representing state transitions between cells is constructed by integrating both static relationships in gene expression space and dynamic relationships in the GRN space. The edges of the undirected graph are then refined using pseudotime inferred based on cell entropy in the GRN space. Finally, the Minimum Spanning Tree (MST) algorithm is applied to derive the cell differentiation trajectory. We validate the accuracy of scGRN-Entropy on eight different real scRNA-seq datasets, demonstrating its superior performance in inferring cell differentiation trajectories through comparative analysis with existing state-of-the-art methods.

DeepKINET: a deep generative model for estimating single-cell RNA splicing and degradation rates

Messenger RNA splicing and degradation are critical for gene expression regulation, the abnormality of which leads to diseases. Previous methods for estimating kinetic rates have limitations, assuming uniform rates across cells. DeepKINET is a deep generative model that estimates splicing and degradation rates at single-cell resolution from scRNA-seq data. DeepKINET outperforms existing methods on simulated and metabolic labeling datasets. Applied to forebrain and breast cancer data, it identifies RNA-binding proteins responsible for kinetic rate diversity. DeepKINET also analyzes the effects of splicing factor mutations on target genes in erythroid lineage cells. DeepKINET effectively reveals cellular heterogeneity in post-transcriptional regulation.

Identification of Meibomian gland stem cell populations and mechanisms of aging

Meibomian glands secrete lipid-rich meibum, which prevents tear evaporation. Aging-related Meibomian gland shrinkage may result in part from stem cell exhaustion and is associated with evaporative dry eye disease, a common condition lacking effective treatment. The identities and niche of Meibomian gland stem cells and the signals controlling their activity are poorly defined. Using snRNA-seq, in vivo lineage tracing, ex vivo live imaging, and genetic studies in mice, we identified markers for stem cell populations that maintain distinct regions of the gland and uncovered Hh signaling as a key regulator of stem cell proliferation. Consistent with this, human Meibomian gland carcinoma exhibited increased Hh signaling. Aged glands displayed decreased Hh and EGF signaling, deficient innervation, and loss of collagen I in niche fibroblasts, indicating that alterations in both glandular epithelial cells and their surrounding microenvironment contribute to age-related degeneration. These findings suggest new approaches to treat aging-associated Meibomian gland loss.

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.

Forseti: a mechanistic and predictive model of the splicing status of scRNA-seq reads

MOTIVATION

Short-read single-cell RNA-sequencing (scRNA-seq) has been used to study cellular heterogeneity, cellular fate, and transcriptional dynamics. Modeling splicing dynamics in scRNA-seq data is challenging, with inherent difficulty in even the seemingly straightforward task of elucidating the splicing status of the molecules from which sequenced fragments are drawn. This difficulty arises, in part, from the limited read length and positional biases, which substantially reduce the specificity of the sequenced fragments. As a result, the splicing status of many reads in scRNA-seq is ambiguous because of a lack of definitive evidence. We are therefore in need of methods that can recover the splicing status of ambiguous reads which, in turn, can lead to more accuracy and confidence in downstream analyses.

RESULTS

We develop Forseti, a predictive model to probabilistically assign a splicing status to scRNA-seq reads. Our model has two key components. First, we train a binding affinity model to assign a probability that a given transcriptomic site is used in fragment generation. Second, we fit a robust fragment length distribution model that generalizes well across datasets deriving from different species and tissue types. Forseti combines these two trained models to predict the splicing status of the molecule of origin of reads by scoring putative fragments that associate each alignment of sequenced reads with proximate potential priming sites. Using both simulated and experimental data, we show that our model can precisely predict the splicing status of many reads and identify the true gene origin of multi-gene mapped reads.

AVAILABILITY AND IMPLEMENTATION

Forseti and the code used for producing the results are available at https://github.com/COMBINE-lab/forseti under a BSD 3-clause license.

Spatial transition tensor of single cells

Spatial transcriptomics and messenger RNA splicing encode extensive spatiotemporal information for cell states and transitions. The current lineage-inference methods either lack spatial dynamics for state transition or cannot capture different dynamics associated with multiple cell states and transition paths. Here we present spatial transition tensor (STT), a method that uses messenger RNA splicing and spatial transcriptomes through a multiscale dynamical model to characterize multistability in space. By learning a four-dimensional transition tensor and spatial-constrained random walk, STT reconstructs cell-state-specific dynamics and spatial state transitions via both short-time local tensor streamlines between cells and long-time transition paths among attractors. Benchmarking and applications of STT on several transcriptome datasets via multiple technologies on epithelial-mesenchymal transitions, blood development, spatially resolved mouse brain and chicken heart development, indicate STT's capability in recovering cell-state-specific dynamics and their associated genes not seen using existing methods. Overall, STT provides a consistent multiscale description of single-cell transcriptome data across multiple spatiotemporal scales.

Reconstructing developmental trajectories using latent dynamical systems and time-resolved transcriptomics

The snapshot nature of single-cell transcriptomics presents a challenge for studying the dynamics of cell fate decisions. Metabolic labeling and splicing can provide temporal information at single-cell level, but current methods have limitations. Here, we present a framework that overcomes these limitations: experimentally, we developed sci-FATE2, an optimized method for metabolic labeling with increased data quality, which we used to profile 45,000 embryonic stem (ES) cells differentiating into neural tube identities. Computationally, we developed a two-stage framework for dynamical modeling: VelvetVAE, a variational autoencoder (VAE) for velocity inference that outperforms all other tools tested, and VelvetSDE, a neural stochastic differential equation (nSDE) framework for simulating trajectory distributions. These recapitulate underlying dataset distributions and capture features such as decision boundaries between alternative fates and fate-specific gene expression. These methods recast single-cell analyses from descriptions of observed data to models of the dynamics that generated them, providing a framework for investigating developmental fate decisions.

PhyloVelo enhances transcriptomic velocity field mapping using monotonically expressed genes

Single-cell RNA sequencing (scRNA-seq) is a powerful approach for studying cellular differentiation, but accurately tracking cell fate transitions can be challenging, especially in disease conditions. Here we introduce PhyloVelo, a computational framework that estimates the velocity of transcriptomic dynamics by using monotonically expressed genes (MEGs) or genes with expression patterns that either increase or decrease, but do not cycle, through phylogenetic time. Through integration of scRNA-seq data with lineage information, PhyloVelo identifies MEGs and reconstructs a transcriptomic velocity field. We validate PhyloVelo using simulated data and Caenorhabditis elegans ground truth data, successfully recovering linear, bifurcated and convergent differentiations. Applying PhyloVelo to seven lineage-traced scRNA-seq datasets, generated using CRISPR-Cas9 editing, lentiviral barcoding or immune repertoire profiling, demonstrates its high accuracy and robustness in inferring complex lineage trajectories while outperforming RNA velocity. Additionally, we discovered that MEGs across tissues and organisms share similar functions in translation and ribosome biogenesis.

SOXC are critical regulators of adult bone mass

Pivotal in many ways for human health, the control of adult bone mass is governed by complex, incompletely understood crosstalk namely between mesenchymal stem cells, osteoblasts and osteoclasts. The SOX4, SOX11 and SOX12 (SOXC) transcription factors were previously shown to control many developmental processes, including skeletogenesis, and SOX4 was linked to osteoporosis, but how SOXC control adult bone mass remains unknown. Using SOXC loss- and gain-of-function mouse models, we show here that SOXC redundantly promote prepubertal cortical bone mass strengthening whereas only SOX4 mitigates adult trabecular bone mass accrual in early adulthood and subsequent maintenance. SOX4 favors bone resorption over formation by lowering osteoblastogenesis and increasing osteoclastogenesis. Single-cell transcriptomics reveals its prevalent expression in Lepr+ mesenchymal cells and ability to upregulate genes for prominent anti-osteoblastogenic and pro-osteoclastogenic factors, including interferon signaling-related chemokines, contributing to these adult stem cells' secretome. SOXC, with SOX4 predominantly, are thus key regulators of adult bone mass.

Forseti: A mechanistic and predictive model of the splicing status of scRNA-seq reads

Motivation

Short-read single-cell RNA-sequencing (scRNA-seq) has been used to study cellular heterogeneity, cellular fate, and transcriptional dynamics. Modeling splicing dynamics in scRNA-seq data is challenging, with inherent difficulty in even the seemingly straightforward task of elucidating the splicing status of the molecules from which sequenced fragments are drawn. This difficulty arises, in part, from the limited read length and positional biases, which substantially reduce the specificity of the sequenced fragments. As a result, the splicing status of many reads in scRNA-seq is ambiguous because of a lack of definitive evidence. We are therefore in need of methods that can recover the splicing status of ambiguous reads which, in turn, can lead to more accuracy and confidence in downstream analyses.

Results

We develop Forseti, a predictive model to probabilistically assign a splicing status to scRNA-seq reads. Our model has two key components. First, we train a binding affinity model to assign a probability that a given transcriptomic site is used in fragment generation. Second, we fit a robust fragment length distribution model that generalizes well across datasets deriving from different species and tissue types. Forseti combines these two trained models to predict the splicing status of the molecule of origin of reads by scoring putative fragments that associate each alignment of sequenced reads with proximate potential priming sites. Using both simulated and experimental data, we show that our model can precisely predict the splicing status of reads and identify the true gene origin of multi-gene mapped reads.

Availability

Forseti and the code used for producing the results are available at https://github.com/COMBINE-lab/forseti under a BSD 3-clause license.

Resolving cellular dynamics using single-cell temporal transcriptomics

Cellular dynamics, the transition of a cell from one state to another, is central to understanding developmental processes and disease progression. Single-cell transcriptomics has been pushing the frontiers of cellular dynamics studies into a genome-wide and single-cell level. While most single-cell RNA sequencing approaches are disruptive and only provide a snapshot of cell states, the dynamics of a cell could be reconstructed by either exploiting temporal information hiding in the transcriptomics data or integrating additional information. In this review, we describe these approaches, highlighting their underlying principles, key assumptions, and the rationality to interpret the results as models. We also discuss the recently emerging nondisruptive live-cell transcriptomics methods, which are highly complementary to the computational models for their assumption-free nature.

Studying temporal dynamics of single cells: expression, lineage and regulatory networks

Learning how multicellular organs are developed from single cells to different cell types is a fundamental problem in biology. With the high-throughput scRNA-seq technology, computational methods have been developed to reveal the temporal dynamics of single cells from transcriptomic data, from phenomena on cell trajectories to the underlying mechanism that formed the trajectory. There are several distinct families of computational methods including Trajectory Inference (TI), Lineage Tracing (LT), and Gene Regulatory Network (GRN) Inference which are involved in such studies. This review summarizes these computational approaches which use scRNA-seq data to study cell differentiation and cell fate specification as well as the advantages and limitations of different methods. We further discuss how GRNs can potentially affect cell fate decisions and trajectory structures.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12551-023-01090-5.

DeepVelo: deep learning extends RNA velocity to multi-lineage systems with cell-specific kinetics

Existing RNA velocity estimation methods strongly rely on predefined dynamics and cell-agnostic constant transcriptional kinetic rates, assumptions often violated in complex and heterogeneous single-cell RNA sequencing (scRNA-seq) data. Using a graph convolution network, DeepVelo overcomes these limitations by generalizing RNA velocity to cell populations containing time-dependent kinetics and multiple lineages. DeepVelo infers time-varying cellular rates of transcription, splicing, and degradation, recovers each cell's stage in the differentiation process, and detects functionally relevant driver genes regulating these processes. Application to various developmental and pathogenic processes demonstrates DeepVelo's capacity to study complex differentiation and lineage decision events in heterogeneous scRNA-seq data.

Improving the RNA velocity approach with single-cell RNA lifecycle (nascent, mature and degrading RNAs) sequencing technologies

We presented an experimental method called FLOUR-seq, which combines BD Rhapsody and nanopore sequencing to detect the RNA lifecycle (including nascent, mature, and degrading RNAs) in cells. Additionally, we updated our HIT-scISOseq V2 to discover a more accurate RNA lifecycle using 10x Chromium and Pacbio sequencing. Most importantly, to explore how single-cell full-length RNA sequencing technologies could help improve the RNA velocity approach, we introduced a new algorithm called 'Region Velocity' to more accurately configure cellular RNA velocity. We applied this algorithm to study spermiogenesis and compared the performance of FLOUR-seq with Pacbio-based HIT-scISOseq V2. Our findings demonstrated that 'Region Velocity' is more suitable for analyzing single-cell full-length RNA data than traditional RNA velocity approaches. These novel methods could be useful for researchers looking to discover full-length RNAs in single cells and comprehensively monitor RNA lifecycle in cells.

Transflammation in tissue regeneration and response to injury: How cell-autonomous inflammatory signaling mediates cell plasticity

Inflammation is a first responder against injury and infection and is also critical for the regeneration and repair of tissue after injury. The role of professional immune cells in tissue healing is well characterized. Professional immune cells respond to pathogens with humoral and cytotoxic responses; remove cellular debris through efferocytosis; secrete angiogenic cytokines and growth factors to repair the microvasculature and parenchyma. However, non-immune cells are also capable of responding to damage or pathogens. Non-immune somatic cells express pattern recognition receptors (PRRs) to detect pathogen-associated molecular patterns (PAMPs) and damage-associated molecular patterns (DAMPs). The PRRs activation leads to the release of inflammatory cytokines required for tissue defense and repair. Notably, the activation of PRRs also triggers epigenetic changes that promote DNA accessibility and cellular plasticity. Thus, non-immune cells directly respond to the local inflammatory cues and can undergo phenotypic modifications or even cell lineage transitions to facilitate tissue regeneration. This review will focus on the novel role of cell-autonomous inflammatory signaling in mediating cell plasticity, a process which is termed transflammation. We will discuss the regulation of this process by changes in the functions and expression levels of epigenetic modifiers, as well as metabolic and ROS/RNS-mediated epigenetic modulation of DNA accessibility during cell fate transition. We will highlight the recent technological developments in detecting cell plasticity and potential therapeutic applications of transflammation in tissue regeneration.

Integrative Analysis of Machine Learning and Molecule Docking Simulations for Ischemic Stroke Diagnosis and Therapy

Due to the narrow therapeutic window and high mortality of ischemic stroke, it is of great significance to investigate its diagnosis and therapy. We employed weighted gene coexpression network analysis (WGCNA) to ascertain gene modules related to stroke and used the maSigPro R package to seek the time-dependent genes in the progression of stroke. Three machine learning algorithms were further employed to identify the feature genes of stroke. A nomogram model was built and applied to evaluate the stroke patients. We analyzed single-cell RNA sequencing (scRNA-seq) data to discern microglia subclusters in ischemic stroke. The RNA velocity, pseudo time, and gene set enrichment analysis (GSEA) were performed to investigate the relationship of microglia subclusters. Connectivity map (CMap) analysis and molecule docking were used to screen a therapeutic agent for stroke. A nomogram model based on the feature genes showed a clinical net benefit and enabled an accurate evaluation of stroke patients. The RNA velocity and pseudo time analysis showed that microglia subcluster 0 would develop toward subcluster 2 within 24 h from stroke onset. The GSEA showed that the function of microglia subcluster 0 was opposite to that of subcluster 2. AZ_628, which screened from CMap analysis, was found to have lower binding energy with Mmp12, Lgals3, Fam20c, Capg, Pkm2, Sdc4, and Itga5 in microglia subcluster 2 and maybe a therapeutic agent for the poor development of microglia subcluster 2 after stroke. Our study presents a nomogram model for stroke diagnosis and provides a potential molecule agent for stroke therapy.

Epigenetic Induction of Smooth Muscle Cell Phenotypic Alterations in Aortic Aneurysms and Dissections

BACKGROUND

Smooth muscle cell (SMC) phenotypic switching has been increasingly detected in aortic aneurysm and dissection (AAD) tissues. However, the diverse SMC phenotypes in AAD tissues and the mechanisms driving SMC phenotypic alterations remain to be identified.

METHODS

We examined the transcriptomic and epigenomic dynamics of aortic SMC phenotypic changes in mice with angiotensin II-induced AAD by using single-cell RNA sequencing and single-cell sequencing assay for transposase-accessible chromatin. SMC phenotypic alteration in aortas from patients with ascending thoracic AAD was examined by using single-cell RNA sequencing analysis.

RESULTS

Single-cell RNA sequencing analysis revealed that aortic stress induced the transition of SMCs from a primary contractile phenotype to proliferative, extracellular matrix-producing, and inflammatory phenotypes. Lineage tracing showed the complete transformation of SMCs to fibroblasts and macrophages. Single-cell sequencing assay for transposase-accessible chromatin analysis indicated that these phenotypic alterations were controlled by chromatin remodeling marked by the reduced chromatin accessibility of contractile genes and the induced chromatin accessibility of genes involved in proliferation, extracellular matrix, and inflammation. IRF3 (interferon regulatory factor 3), a proinflammatory transcription factor activated by cytosolic DNA, was identified as a key driver of the transition of aortic SMCs from a contractile phenotype to an inflammatory phenotype. In cultured SMCs, cytosolic DNA signaled through its sensor STING (stimulator of interferon genes)-TBK1 (tank-binding kinase 1) to activate IRF3, which bound and recruited EZH2 (enhancer of zeste homolog 2) to contractile genes to induce repressive H3K27me3 modification and gene suppression. In contrast, double-stranded DNA-STING-IRF3 signaling induced inflammatory gene expression in SMCs. In Sting-/- mice, the aortic stress-induced transition of SMCs into an inflammatory phenotype was prevented, and SMC populations were preserved. Finally, profound SMC phenotypic alterations toward diverse directions were detected in human ascending thoracic AAD tissues.

CONCLUSIONS

Our study reveals the dynamic epigenetic induction of SMC phenotypic alterations in AAD. DNA damage and cytosolic leakage drive SMCs from a contractile phenotype to an inflammatory phenotype.

Neural cell diversity in the light of single-cell transcriptomics

The development of highly parallel and affordable high-throughput single-cell transcriptomics technologies has revolutionized our understanding of brain complexity. These methods have been used to build cellular maps of the brain, its different regions, and catalog the diversity of cells in each of them during development, aging and even in disease. Now we know that cellular diversity is way beyond what was previously thought. Single-cell transcriptomics analyses have revealed that cell types previously considered homogeneous based on imaging techniques differ depending on several factors including sex, age and location within the brain. The expression profiles of these cells have also been exploited to understand which are the regulatory programs behind cellular diversity and decipher the transcriptional pathways driving them. In this review, we summarize how single-cell transcriptomics have changed our view on the cellular diversity in the human brain, and how it could impact the way we study neurodegenerative diseases. Moreover, we describe the new computational approaches that can be used to study cellular differentiation and gain insight into the functions of individual cell populations under different conditions and their alterations in disease.