NASC-seq monitors RNA synthesis in single cells

Single-cell new RNA sequencing reveals principles of transcription at the resolution of individual bursts

Analyses of transcriptional bursting from single-cell RNA-sequencing data have revealed patterns of variation and regulation in the kinetic parameters that could be inferred. Here we profiled newly transcribed (4-thiouridine-labelled) RNA across 10,000 individual primary mouse fibroblasts to more broadly infer bursting kinetics and coordination. We demonstrate that inference from new RNA profiles could separate the kinetic parameters that together specify the burst size, and that the synthesis rate (and not the transcriptional off rate) controls the burst size. Importantly, transcriptome-wide inference of transcriptional on and off rates provided conclusive evidence that RNA polymerase II transcribes genes in bursts. Recent reports identified examples of transcriptional co-bursting, yet no global analyses have been performed. The deep new RNA profiles we generated with allelic resolution demonstrated that co-bursting rarely appears more frequently than expected by chance, except for certain gene pairs, notably paralogues located in close genomic proximity. Altogether, new RNA single-cell profiling critically improves the inference of transcriptional bursting and provides strong evidence for independent transcriptional bursting of mammalian genes.

Global transcription regulation revealed from dynamical correlations in time-resolved single-cell RNA sequencing

Single-cell transcriptomics reveals significant variations in transcriptional activity across cells. Yet, it remains challenging to identify mechanisms of transcription dynamics from static snapshots. It is thus still unknown what drives global transcription dynamics in single cells. We present a stochastic model of gene expression with cell size- and cell cycle-dependent rates in growing and dividing cells that harnesses temporal dimensions of single-cell RNA sequencing through metabolic labeling protocols and cel lcycle reporters. We develop a parallel and highly scalable approximate Bayesian computation method that corrects for technical variation and accurately quantifies absolute burst frequency, burst size, and degradation rate along the cell cycle at a transcriptome-wide scale. Using Bayesian model selection, we reveal scaling between transcription rates and cell size and unveil waves of gene regulation across the cell cycle-dependent transcriptome. Our study shows that stochastic modeling of dynamical correlations identifies global mechanisms of transcription regulation. A record of this paper's transparent peer review process is included in the supplemental information.

Profiling cell identity and tissue architecture with single-cell and spatial transcriptomics

Single-cell transcriptomics has broadened our understanding of cellular diversity and gene expression dynamics in healthy and diseased tissues. Recently, spatial transcriptomics has emerged as a tool to contextualize single cells in multicellular neighbourhoods and to identify spatially recurrent phenotypes, or ecotypes. These technologies have generated vast datasets with targeted-transcriptome and whole-transcriptome profiles of hundreds to millions of cells. Such data have provided new insights into developmental hierarchies, cellular plasticity and diverse tissue microenvironments, and spurred a burst of innovation in computational methods for single-cell analysis. In this Review, we discuss recent advancements, ongoing challenges and prospects in identifying and characterizing cell states and multicellular neighbourhoods. We discuss recent progress in sample processing, data integration, identification of subtle cell states, trajectory modelling, deconvolution and spatial analysis. Furthermore, we discuss the increasing application of deep learning, including foundation models, in analysing single-cell and spatial transcriptomics data. Finally, we discuss recent applications of these tools in the fields of stem cell biology, immunology, and tumour biology, and the future of single-cell and spatial transcriptomics in biological research and its translation to the clinic.

Dissecting key regulators of transcriptome kinetics through scalable single-cell RNA profiling of pooled CRISPR screens

We present a combinatorial indexing method, PerturbSci-Kinetics, for capturing whole transcriptomes, nascent transcriptomes and single guide RNA (sgRNA) identities across hundreds of genetic perturbations at the single-cell level. Profiling a pooled CRISPR screen targeting various biological processes, we show the gene expression regulation during RNA synthesis, processing and degradation, miRNA biogenesis and mitochondrial mRNA processing, systematically decoding the genome-wide regulatory network that underlies RNA temporal dynamics at scale.

Stromal Rigidity Stress Accelerates Pancreatic Intraepithelial Neoplasia Progression and Chromosomal Instability via Nuclear Protein Tyrosine Kinase 2 Localization

Because the mechanotransduction by stromal stiffness stimulates the rupture and repair of the nuclear envelope in pancreatic progenitor cells, accumulated genomic aberrations are under selection in the tumor microenvironment. Analysis of cell growth, micronuclei, and phosphorylated Ser-139 residue of the histone variant H2AX (γH2AX) foci linked to mechanotransduction pressure in vivo during serial orthotopic passages of mouse KrasLSL-G12D/+;Trp53flox/flox;Pdx1-Cre (KPC) cancer cells in the tumor and in migrating through the size-restricted 3-μm micropores. To search for pancreatic cancer cell-of-origin, analysis of single-cell data sets revealed that the extracellular matrix shaped an alternate route of acinar-ductal transdifferentiation of acinar cells into topoisomerase II α (TOP2A)-overexpressing cancer cells and derived subclusters with copy number amplifications in MYC-PTK2 (protein tyrosine kinase 2) locus and PIK3CA. High-PTK2 expression is associated with 171 differentially methylated CpG loci, 319 differentially expressed genes, and poor overall survival in The Cancer Genome Atlas-Pancreatic Adenocarcinoma cohort. Abolished RGD-integrin signaling by disintegrin KG blocked the PTK2 phosphorylation, increased cancer apoptosis, decreased vav guanine nucleotide exchange factor 1 (VAV1) expression, and prolonged overall survival in the KPC mice. Reduction of α-smooth muscle actin deposition in the CD248 knockout KPC mice remodeled the tissue stroma and down-regulated TOP2A expression in the epithelium. In summary, stromal stiffness induced the onset of cancer cells-of-origin by ectopic TOP2A expression, and the genomic amplification of MYC-PTK2 locus via alternative transdifferentiation of pancreatic progenitor cells is the vulnerability useful for disintegrin KG treatment.

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.

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.

A dynamical perspective: moving towards mechanism in single-cell transcriptomics

As the field of single-cell transcriptomics matures, research is shifting focus from phenomenological descriptions of cellular phenotypes to a mechanistic understanding of the gene regulation underneath. This perspective considers the value of capturing dynamical information at single-cell resolution for gaining mechanistic insight; reviews the available technologies for recording and inferring temporal information in single cells; and explores whether better dynamical resolution is sufficient to adequately capture the causal relationships driving complex biological systems. This article is part of a discussion meeting issue 'Causes and consequences of stochastic processes in development and disease'.

Imprinted X chromosome inactivation at the gamete-to-embryo transition

In mammals, dosage compensation involves two parallel processes: (1) X inactivation, which equalizes X chromosome dosage between males and females, and (2) X hyperactivation, which upregulates the active X for X-autosome balance. The field currently favors models whereby dosage compensation initiates "de novo" during mouse development. Here, we develop "So-Smart-seq" to revisit the question and interrogate a comprehensive transcriptome including noncoding genes and repeats in mice. Intriguingly, de novo silencing pertains only to a subset of Xp genes. Evolutionarily older genes and repetitive elements demonstrate constitutive Xp silencing, adopt distinct signatures, and do not require Xist to initiate silencing. We trace Xp silencing backward in developmental time to meiotic sex chromosome inactivation in the male germ line and observe that Xm hyperactivation is timed to Xp silencing on a gene-by-gene basis. Thus, during the gamete-to-embryo transition, older Xp genes are transmitted in a "pre-inactivated" state. These findings have implications for the evolution of imprinting.

Cell-type-specific mRNA transcription and degradation kinetics in zebrafish embryogenesis from metabolically labeled single-cell RNA-seq

During embryonic development, pluripotent cells assume specialized identities by adopting particular gene expression profiles. However, systematically dissecting the relative contributions of mRNA transcription and degradation to shaping those profiles remains challenging, especially within embryos with diverse cellular identities. Here, we combine single-cell RNA-Seq and metabolic labeling to capture temporal cellular transcriptomes of zebrafish embryos where newly-transcribed (zygotic) and pre-existing (maternal) mRNA can be distinguished. We introduce kinetic models to quantify mRNA transcription and degradation rates within individual cell types during their specification. These models reveal highly varied regulatory rates across thousands of genes, coordinated transcription and destruction rates for many transcripts, and link differences in degradation to specific sequence elements. They also identify cell-type-specific differences in degradation, namely selective retention of maternal transcripts within primordial germ cells and enveloping layer cells, two of the earliest specified cell types. Our study provides a quantitative approach to study mRNA regulation during a dynamic spatio-temporal response.

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.

SLAM-Drop-seq reveals mRNA kinetic rates throughout the cell cycle

RNA abundance is tightly regulated in eukaryotic cells by modulating the kinetic rates of RNA production, processing, and degradation. To date, little is known about time‐dependent kinetic rates during dynamic processes. Here, we present SLAM‐Drop‐seq, a method that combines RNA metabolic labeling and alkylation of modified nucleotides in methanol‐fixed cells with droplet‐based sequencing to detect newly synthesized and preexisting mRNAs in single cells. As a first application, we sequenced 7280 HEK293 cells and calculated gene‐specific kinetic rates during the cell cycle using the novel package Eskrate. Of the 377 robust‐cycling genes that we identified, only a minor fraction is regulated solely by either dynamic transcription or degradation (6 and 4%, respectively). By contrast, the vast majority (89%) exhibit dynamically regulated transcription and degradation rates during the cell cycle. Our study thus shows that temporally regulated mRNA degradation is fundamental for the correct expression of a majority of cycling genes. SLAM‐Drop‐seq, combined with Eskrate, is a powerful approach to understanding the underlying mRNA kinetics of single‐cell gene expression dynamics in continuous biological processes.

Decoding Human Biology and Disease Using Single-cell Omics Technologies

Over the past decade, advances in single-cell omics (SCO) technologies have enabled the investigation of cellular heterogeneity at an unprecedented resolution and scale, opening a new avenue for understanding human biology and disease. In this review, we summarize the developments of sequencing-based SCO technologies and computational methods, and focus on considerable insights acquired from SCO sequencing studies to understand normal and diseased properties, with a particular emphasis on cancer research. We also discuss the technological improvements of SCO and its possible contribution to fundamental research of the human, as well as its great potential in clinical diagnoses and personalized therapies of human disease.

Transcriptomic forecasting with neural ordinary differential equations

Single-cell transcriptomics technologies can uncover changes in the molecular states that underlie cellular phenotypes. However, understanding the dynamic cellular processes requires extending from inferring trajectories from snapshots of cellular states to estimating temporal changes in cellular gene expression. To address this challenge, we have developed a neural ordinary differential-equation-based method, RNAForecaster, for predicting gene expression states in single cells for multiple future time steps in an embedding-independent manner. We demonstrate that RNAForecaster can accurately predict future expression states in simulated single-cell transcriptomic data with cellular tracking over time. We then show that by using metabolic labeling single-cell RNA sequencing (scRNA-seq) data from constitutively dividing cells, RNAForecaster accurately recapitulates many of the expected changes in gene expression during progression through the cell cycle over a 3-day period. Thus, RNAForecaster enables short-term estimation of future expression states in biological systems from high-throughput datasets with temporal information.

Understanding tumour endothelial cell heterogeneity and function from single-cell omics

Anti-angiogenic therapies (AATs) are used to treat different types of cancers. However, their success is limited owing to insufficient efficacy and resistance. Recently, single-cell omics studies of tumour endothelial cells (TECs) have provided new mechanistic insight. Here, we overview the heterogeneity of human TECs of all tumour types studied to date, at the single-cell level. Notably, most human tumour types contain varying numbers but only a small population of angiogenic TECs, the presumed targets of AATs, possibly contributing to the limited efficacy of and resistance to AATs. In general, TECs are heterogeneous within and across all tumour types, but comparing TEC phenotypes across tumours is currently challenging, owing to the lack of a uniform nomenclature for endothelial cells and consistent single-cell analysis protocols, urgently raising the need for a more consistent approach. Nonetheless, across most tumour types, universal TEC markers (ACKR1, PLVAP and IGFBP3) can be identified. Besides angiogenesis, biological processes such as immunomodulation and extracellular matrix organization are among the most commonly predicted enriched signatures of TECs across different tumour types. Although angiogenesis and extracellular matrix targets have been considered for AAT (without the hoped success), the immunomodulatory properties of TECs have not been fully considered as a novel anticancer therapeutic approach. Therefore, we also discuss progress, limitations, solutions and novel targets for AAT development.

Synergising single-cell resolution and 4sU labelling boosts inference of transcriptional bursting

Despite the recent rise of RNA-seq datasets combining single-cell (sc) resolution with 4-thiouridine (4sU) labelling, analytical methods exploiting their power to dissect transcriptional bursting are lacking. Here, we present a mathematical model and Bayesian inference implementation to facilitate genome-wide joint parameter estimation and confidence quantification (R package: burstMCMC). We demonstrate that, unlike conventional scRNA-seq, 4sU scRNA-seq resolves temporal parameters and furthermore boosts inference of dimensionless parameters via a synergy between single-cell resolution and 4sU labelling. We apply our method to published 4sU scRNA-seq data and linked with ChIP-seq data, we uncover previously obscured associations between different parameters and histone modifications.

grandR: a comprehensive package for nucleotide conversion RNA-seq data analysis

Metabolic labeling of RNA is a powerful technique for studying the temporal dynamics of gene expression. Nucleotide conversion approaches greatly facilitate the generation of data but introduce challenges for their analysis. Here we present grandR, a comprehensive package for quality control, differential gene expression analysis, kinetic modeling, and visualization of such data. We compare several existing methods for inference of RNA synthesis rates and half-lives using progressive labeling time courses. We demonstrate the need for recalibration of effective labeling times and introduce a Bayesian approach to study the temporal dynamics of RNA using snapshot experiments.

PerturbSci-Kinetics: Dissecting key regulators of transcriptome kinetics through scalable single-cell RNA profiling of pooled CRISPR screens

Here we described PerturbSci-Kinetics, a novel combinatorial indexing method for capturing three-layer single-cell readout (i.e., whole transcriptomes, nascent transcriptomes, sgRNA identities) across hundreds of genetic perturbations. Through PerturbSci-Kinetics profiling of pooled CRISPR screens targeting a variety of biological processes, we were able to decipher the complexity of RNA regulations at multiple levels (e.g., synthesis, processing, degradation), and revealed key regulators involved in miRNA and mitochondrial RNA processing pathways. Our technique opens the possibility of systematically decoding the genome-wide regulatory network underlying RNA temporal dynamics at scale and cost-effectively.

Spatiotemporally resolved transcriptomics reveals the subcellular RNA kinetic landscape

Spatiotemporal regulation of the cellular transcriptome is crucial for proper protein expression and cellular function. However, the intricate subcellular dynamics of RNA remain obscured due to the limitations of existing transcriptomics methods. Here, we report TEMPOmap-a method that uncovers subcellular RNA profiles across time and space at the single-cell level. TEMPOmap integrates pulse-chase metabolic labeling with highly multiplexed three-dimensional in situ sequencing to simultaneously profile the age and location of individual RNA molecules. Using TEMPOmap, we constructed the subcellular RNA kinetic landscape in various human cells from transcription and translocation to degradation. Clustering analysis of RNA kinetic parameters across single cells revealed 'kinetic gene clusters' whose expression patterns were shaped by multistep kinetic sculpting. Importantly, these kinetic gene clusters are functionally segregated, suggesting that subcellular RNA kinetics are differentially regulated in a cell-state- and cell-type-dependent manner. Spatiotemporally resolved transcriptomics provides a gateway to uncovering new spatiotemporal gene regulation principles.

Single-cell temporal dynamics reveals the relative contributions of transcription and degradation to cell-type specific gene expression in zebrafish embryos

During embryonic development, pluripotent cells assume specialized identities by adopting particular gene expression profiles. However, systematically dissecting the underlying regulation of mRNA transcription and degradation remains a challenge, especially within whole embryos with diverse cellular identities. Here, we collect temporal cellular transcriptomes of zebrafish embryos, and decompose them into their newly-transcribed (zygotic) and pre-existing (maternal) mRNA components by combining single-cell RNA-Seq and metabolic labeling. We introduce kinetic models capable of quantifying regulatory rates of mRNA transcription and degradation within individual cell types during their specification. These reveal different regulatory rates between thousands of genes, and sometimes between cell types, that shape spatio-temporal expression patterns. Transcription drives most cell-type restricted gene expression. However, selective retention of maternal transcripts helps to define the gene expression profiles of germ cells and enveloping layer cells, two of the earliest specified cell-types. Coordination between transcription and degradation restricts expression of maternal-zygotic genes to specific cell types or times, and allows the emergence of spatio-temporal patterns when overall mRNA levels are held relatively constant. Sequence-based analysis links differences in degradation to specific sequence motifs. Our study reveals mRNA transcription and degradation events that control embryonic gene expression, and provides a quantitative approach to study mRNA regulation during a dynamic spatio-temporal response.

SLAM-seq reveals early transcriptomic response mechanisms upon glutamine deprivation in Chinese hamster ovary cells

Mammalian cells frequently encounter subtle perturbations during recombinant protein production. Identifying the genetic factors that govern the cellular stress response can facilitate targeted genetic engineering to obtain production cell lines that demonstrate a higher stress tolerance. To simulate nutrient stress, Chinese hamster ovary (CHO) cells were transferred into a glutamine(Q)-free medium and transcriptional dynamics using thiol(SH)-linked alkylation for the metabolic sequencing of RNA (SLAM-seq) along with standard RNA-seq of stressed and unstressed cells were investigated. The SLAM-seq method allows differentiation between actively transcribed, nascent mRNA, and total (previously present) mRNA in the sample, adding an additional, time-resolved layer to classic RNA-sequencing. The cells tackle amino acid (AA) limitation by inducing the integrated stress response (ISR) signaling pathway, reflected in Atf4 overexpression in the early hours post Q deprivation, leading to subsequent activation of its targets, Asns, Atf3, Ddit3, Eif4ebp1, Gpt2, Herpud1, Slc7a1, Slc7a11, Slc38a2, Trib3, and Vegfa. The GCN2-eIF2α-ATF4 pathway is confirmed by a significant halt in transcription of translation-related genes at 24 h post Q deprivation. The downregulation of lipid synthesis indicates the inhibition of the mTOR pathway, further confirmed by overexpression of Sesn2. Furthermore, SLAM-seq detects short-lived transcription factors, such as Egr1, that would have been missed in standard experimental designs with RNA-seq. Our results describe the successful establishment of SLAM-seq in CHO cells and therefore facilitate its future use in other scenarios where dynamic transcriptome profiling in CHO cells is essential.

Well-TEMP-seq as a microwell-based strategy for massively parallel profiling of single-cell temporal RNA dynamics

Single-cell RNA sequencing (scRNA-seq) reveals the transcriptional heterogeneity of cells, but the static snapshots fail to reveal the time-resolved dynamics of transcription. Herein, we develop Well-TEMP-seq, a high-throughput, cost-effective, accurate, and efficient method for massively parallel profiling the temporal dynamics of single-cell gene expression. Well-TEMP-seq combines metabolic RNA labeling with scRNA-seq method Well-paired-seq to distinguish newly transcribed RNAs marked by T-to-C substitutions from pre-existing RNAs in each of thousands of single cells. The Well-paired-seq chip ensures a high single cell/barcoded bead pairing rate (~80%) and the improved alkylation chemistry on beads greatly alleviates chemical conversion-induced cell loss (~67.5% recovery). We further apply Well-TEMP-seq to profile the transcriptional dynamics of colorectal cancer cells exposed to 5-AZA-CdR, a DNA-demethylating drug. Well-TEMP-seq unbiasedly captures the RNA dynamics and outperforms the splicing-based RNA velocity method. We anticipate that Well-TEMP-seq will be broadly applicable to unveil the dynamics of single-cell gene expression in diverse biological processes.

A single-cell RNA labeling strategy for measuring stress response upon tissue dissociation

Tissue dissociation, a crucial step in single-cell sample preparation, can alter the transcriptional state of a sample through the intrinsic cellular stress response. Here we demonstrate a general approach for measuring transcriptional response during sample preparation. In our method, transcripts made during dissociation are labeled for later identification upon sequencing. We found general as well as cell-type-specific dissociation response programs in zebrafish larvae, and we observed sample-to-sample variation in the dissociation response of mouse cardiomyocytes despite well-controlled experimental conditions. Finally, we showed that dissociation of the mouse hippocampus can lead to the artificial activation of microglia. In summary, our approach facilitates experimental optimization of dissociation procedures as well as computational removal of transcriptional perturbation response.

Determining RNA Natural Modifications and Nucleoside Analog-Labeled Sites by a Chemical/Enzyme-Induced Base Mutation Principle

The natural chemical modifications of messenger RNA (mRNA) in living organisms have shown essential roles in both physiology and pathology. The mapping of mRNA modifications is critical for interpreting their biological functions. In another dimension, the synthesized nucleoside analogs can enable chemical labeling of cellular mRNA through a metabolic pathway, which facilitates the study of RNA dynamics in a pulse-chase manner. In this regard, the sequencing tools for mapping both natural modifications and nucleoside tags on mRNA at single base resolution are highly necessary. In this work, we review the progress of chemical sequencing technology for determining both a variety of naturally occurring base modifications mainly on mRNA and a few on transfer RNA and metabolically incorporated artificial base analogs on mRNA, and further discuss the problems and prospects in the field.

Application of nucleoside or nucleotide analogues in RNA dynamics and RNA-binding protein analysis

Cellular RNAs undergo dynamic changes during RNA biological processes, which are tightly orchestrated by RNA-binding proteins (RBPs). Yet, the investigation of RNA dynamics is hurdled by highly abundant steady-state RNAs, which make the signals of dynamic RNAs less detectable. Notably, the exert of nucleoside or nucleotide analogue-based RNA technologies has provided a remarkable platform for RNA dynamics research, revealing diverse unnoticed features in RNA metabolism. In this review, we focus on the application of two types of analogue-based RNA sequencing, antigen-/antibody- and click chemistry-based methodologies, and summarize the RNA dynamics features revealed. Moreover, we discuss emerging single-cell newly transcribed RNA sequencing methodologies based on nucleoside analogue labeling, which provides novel insights into RNA dynamics regulation at single-cell resolution. On the other hand, we also emphasize the identification of RBPs that interact with polyA, non-polyA RNAs, or newly transcribed RNAs and also their associated RNA-binding domains at genomewide level through ultraviolet crosslinking and mass spectrometry in different contexts. We anticipated that further modification and development of these analogue-based RNA and RBP capture technologies will aid in obtaining an unprecedented understanding of RNA biology. This article is categorized under: RNA Interactions with Proteins and Other Molecules > Protein-RNA Recognition RNA Structure and Dynamics > RNA Structure, Dynamics and Chemistry RNA Methods > RNA Analyses in Cells.

Reconstructing data-driven governing equations for cell phenotypic transitions: integration of data science and systems biology

Cells with the same genome can exist in different phenotypes and can change between distinct phenotypes when subject to specific stimuli and microenvironments. Some examples include cell differentiation during development, reprogramming for induced pluripotent stem cells and transdifferentiation, cancer metastasis and fibrosis progression. The regulation and dynamics of cell phenotypic conversion is a fundamental problem in biology, and has a long history of being studied within the formalism of dynamical systems. A main challenge for mechanism-driven modeling studies is acquiring sufficient amount of quantitative information for constraining model parameters. Advances in quantitative experimental approaches, especially high throughput single-cell techniques, have accelerated the emergence of a new direction for reconstructing the governing dynamical equations of a cellular system from quantitative single-cell data, beyond the dominant statistical approaches. Here I review a selected number of recent studies using live- and fixed-cell data and provide my perspective on future development.

Reconstructing data-driven governing equations for cell phenotypic transitions: integration of data science and systems biology

Cells with the same genome can exist in different phenotypes and can change between distinct phenotypes when subject to specific stimuli and microenvironments. Some examples include cell differentiation during development, reprogramming for induced pluripotent stem cells and transdifferentiation, cancer metastasis and fibrosis progression. The regulation and dynamics of cell phenotypic conversion is a fundamental problem in biology, and has a long history of being studied within the formalism of dynamical systems. A main challenge for mechanism-driven modeling studies is acquiring sufficient amount of quantitative information for constraining model parameters. Advances in quantitative experimental approaches, especially high throughput single-cell techniques, have accelerated the emergence of a new direction for reconstructing the governing dynamical equations of a cellular system from quantitative single-cell data, beyond the dominant statistical approaches. Here I review a selected number of recent studies using live- and fixed-cell data and provide my perspective on future development.

Inferring gene regulation from stochastic transcriptional variation across single cells at steady state

Regulatory relationships between transcription factors (TFs) and their target genes lie at the heart of cellular identity and function; however, uncovering these relationships is often labor-intensive and requires perturbations. Here, we propose a principled framework to systematically infer gene regulation for all TFs simultaneously in cells at steady state by leveraging the intrinsic variation in the transcriptional abundance across single cells. Through modeling and simulations, we characterize how transcriptional bursts of a TF gene are propagated to its target genes, including the expected ranges of time delay and magnitude of maximum covariation. We distinguish these temporal trends from the time-invariant covariation arising from cell states, and we delineate the experimental and technical requirements for leveraging these small but meaningful cofluctuations in the presence of measurement noise. While current technology does not yet allow adequate power for definitively detecting regulatory relationships for all TFs simultaneously in cells at steady state, we investigate a small-scale dataset to inform future experimental design. This study supports the potential value of mapping regulatory connections through stochastic variation, and it motivates further technological development to achieve its full potential.

Live-seq enables temporal transcriptomic recording of single cells

Single-cell transcriptomics (scRNA-seq) has greatly advanced our ability to characterize cellular heterogeneity1. However, scRNA-seq requires lysing cells, which impedes further molecular or functional analyses on the same cells. Here, we established Live-seq, a single-cell transcriptome profiling approach that preserves cell viability during RNA extraction using fluidic force microscopy2,3, thus allowing to couple a cell's ground-state transcriptome to its downstream molecular or phenotypic behaviour. To benchmark Live-seq, we used cell growth, functional responses and whole-cell transcriptome read-outs to demonstrate that Live-seq can accurately stratify diverse cell types and states without inducing major cellular perturbations. As a proof of concept, we show that Live-seq can be used to directly map a cell's trajectory by sequentially profiling the transcriptomes of individual macrophages before and after lipopolysaccharide (LPS) stimulation, and of adipose stromal cells pre- and post-differentiation. In addition, we demonstrate that Live-seq can function as a transcriptomic recorder by preregistering the transcriptomes of individual macrophages that were subsequently monitored by time-lapse imaging after LPS exposure. This enabled the unsupervised, genome-wide ranking of genes on the basis of their ability to affect macrophage LPS response heterogeneity, revealing basal Nfkbia expression level and cell cycle state as important phenotypic determinants, which we experimentally validated. Thus, Live-seq can address a broad range of biological questions by transforming scRNA-seq from an end-point to a temporal analysis approach.

Fast and highly sensitive full-length single-cell RNA sequencing using FLASH-seq

We present FLASH-seq (FS), a full-length single-cell RNA sequencing (scRNA-seq) method with increased sensitivity and reduced hands-on time compared to Smart-seq3. The entire FS protocol can be performed in ~4.5 hours, is simple to automate and can be easily miniaturized to decrease resource consumption. The FS protocol can also use unique molecular identifiers (UMIs) for molecule counting while displaying reduced strand-invasion artifacts. FS will be especially useful for characterizing gene expression at high resolution across multiple samples.

FLASH-seq speeds up high-sensitivity scRNA-seq while sparing resources.

Single-cell immunology: Past, present, and future

The immune system is a complex, dynamic, and plastic ecosystem composed of multiple cell types that constantly sense and interact with their local microenvironment to protect from infection and maintain homeostasis. For over a century, great efforts and ingenuity have been applied to the characterization of immune cells and their microenvironments, but traditional marker-based and bulk technologies left key questions unanswered. In the past decade, the advent of single-cell genomic approaches has revolutionized our knowledge of the cellular and molecular makeup of the immune system. In this perspective, we outline the past, present, and future applications of single-cell genomics in immunology and discuss how the integration of multiomics at the single-cell level will pave the way for future advances in immunology research and clinical translation.

Mapping transcriptomic vector fields of single cells

Single-cell (sc)-RNA-seq, together with RNA-velocity and metabolic labeling, reveals cellular states and transitions at unprecedented resolution. Fully exploiting these data, however, requires kinetic models capable of unveiling governing regulatory functions. Here, we introduce an analytical framework dynamo, that infers absolute RNA velocity, reconstructs continuous vector-field functions that predict cell fates, employs differential geometry to extract underlying regulations, and ultimately predicts optimal reprogramming paths and perturbation outcomes. We highlight dynamo’s power to overcome fundamental limitations of conventional splicing-based RNA velocity analyses to enable accurate velocity estimations on a metabolically-labeled human hematopoiesis scRNA-seq dataset. Furthermore, differential geometry analyses reveal mechanisms driving early megakaryocyte appearance and elucidate asymmetrical regulation within the PU.1–GATA1 circuit. Leveraging the Least-Action-Path method, dynamo accurately predicts drivers of numerous hematopoietic transitions. Finally, in silico perturbations predict cell-fate diversions induced by gene perturbations. Dynamo thus represents an important step in advancing quantitative and predictive theories of cell-state transitions.

Temporal modelling using single-cell transcriptomics

Methods for profiling genes at the single-cell level have revolutionized our ability to study several biological processes and systems including development, differentiation, response programmes and disease progression. In many of these studies, cells are profiled over time in order to infer dynamic changes in cell states and types, sets of expressed genes, active pathways and key regulators. However, time-series single-cell RNA sequencing (scRNA-seq) also raises several new analysis and modelling issues. These issues range from determining when and how deep to profile cells, linking cells within and between time points, learning continuous trajectories, and integrating bulk and single-cell data for reconstructing models of dynamic networks. In this Review, we discuss several approaches for the analysis and modelling of time-series scRNA-seq, highlighting their steps, key assumptions, and the types of data and biological questions they are most appropriate for.

In Situ Visualizing Nascent RNA by Exploring DNA-Templated Oxidative Amination of 4-Thiouridine

Tracking and mapping the nascent RNA molecules in cells is essential for deciphering embryonic development and neuronal differentiation. Here, we utilized 4-thiouridine (s4U) as a metabolic tag to label nascent RNA and developed a fluorescence imaging method based on the DNA-templated oxidative amination (DTOA) reaction of s4U. The DTOA reaction occurred between amine-modified DNA and s4U-containing RNA with high sequence specificity and chemical selectivity. Target nascent mRNAs in HeLa cells, including those encoding green fluorescent proteins (GFPs) and endogenous BAG-1, were thus lit up selectively by DTOA-based fluorescence in situ hybridization (DTOA FISH). We believe the DTOA conjugation chemistry shown in this study could be generally applied to investigate the spatial distribution of nascent transcription dynamics in cellular processes.

CellRank for directed single-cell fate mapping

Computational trajectory inference enables the reconstruction of cell state dynamics from single-cell RNA sequencing experiments. However, trajectory inference requires that the direction of a biological process is known, largely limiting its application to differentiating systems in normal development. Here, we present CellRank (https://cellrank.org) for single-cell fate mapping in diverse scenarios, including regeneration, reprogramming and disease, for which direction is unknown. Our approach combines the robustness of trajectory inference with directional information from RNA velocity, taking into account the gradual and stochastic nature of cellular fate decisions, as well as uncertainty in velocity vectors. On pancreas development data, CellRank automatically detects initial, intermediate and terminal populations, predicts fate potentials and visualizes continuous gene expression trends along individual lineages. Applied to lineage-traced cellular reprogramming data, predicted fate probabilities correctly recover reprogramming outcomes. CellRank also predicts a new dedifferentiation trajectory during postinjury lung regeneration, including previously unknown intermediate cell states, which we confirm experimentally.

CellRank infers directed cell state transitions and cell fates incorporating RNA velocity information into a graph based Markov process.

Non-radioactive In Vivo Labeling of RNA with 4-Thiouracil

RNA molecules and their expression dynamics play essential roles in the establishment of complex cellular phenotypes and/or in the rapid cellular adaption to environmental changes. Accordingly, analyzing RNA expression remains an important step to understand the molecular basis controlling the formation of cellular phenotypes, cellular homeostasis or disease progression. Steady-state RNA levels in the cells are controlled by the sum of highly dynamic molecular processes contributing to RNA expression and can be classified in transcription, maturation and degradation. The main goal of analyzing RNA dynamics is to disentangle the individual contribution of these molecular processes to the life cycle of a given RNA under different physiological conditions. In the recent years, the use of nonradioactive nucleotide/nucleoside analogs and improved chemistry, in combination with time-dependent and high-throughput analysis, have greatly expanded our understanding of RNA metabolism across various cell types, organisms, and growth conditions.In this chapter, we describe a step-by-step protocol allowing pulse labeling of RNA with the nonradioactive nucleotide analog, 4-thiouracil , in the eukaryotic model organism Saccharomyces cerevisiae and the model archaeon Haloferax volcanii .

Synthesis of 4-thiouridines with prodrug functionalization for RNA metabolic labeling

Metabolic labeling has emerged as a powerful tool to endow RNA with reactive handles allowing for subsequent chemical derivatization and processing. Recently, thiolated nucleosides, such as 4-thiouridine (4sU), have attracted great interest in metabolic labeling-based RNA sequencing approaches (TUC-seq, SLAM-seq, TimeLapse-seq) to study cellular RNA expression and decay dynamics. For these and other applications (e.g. PAR-CLIP), thus far only the naked nucleoside 4sU has been applied. Here we examined the concept of derivatizing 4sU into a 5′-monophosphate prodrug that would allow for cell permeation and potentially improve labeling efficiency by bypassing the rate-limiting first step of 5′ phosphorylation of the nucleoside into the ultimately bioactive 4sU triphosphate (4sUTP). To this end, we developed robust synthetic routes towards diverse 4sU monophosphate prodrugs. Using metabolic labeling assays, we found that most of the newly introduced 4sU prodrugs were well tolerated by the cells. One derivative, the bis(4-acetyloxybenzyl) 5′-monophosphate of 4sU, was also efficiently incorporated into nascent RNA.

Synthetic access to 4-thiouridine (4sU) derivatives with monophosphate prodrug patterns creates additional possibilities for metabolic labeling of RNA for different applications.

Inference of high-resolution trajectories in single-cell RNA-seq data by using RNA velocity

Trajectory inference (TI) methods infer cell developmental trajectory from single-cell RNA sequencing data. Current TI methods can be categorized into those using RNA velocity information and those using only single-cell gene expression data. The latter type of methods are restricted to certain trajectory structures, and cannot determine cell developmental direction. Recently proposed TI methods using RNA velocity information have limited accuracy. We present CellPath, a method that infers cell trajectories by integrating single-cell gene expression and RNA velocity information. CellPath overcomes the restrictions of TI methods that do not use RNA velocity information: it can find multiple high-resolution trajectories without constraints on the trajectory structure, and can automatically detect the direction of each trajectory path. We evaluate CellPath on both real and simulated datasets and show that CellPath finds more accurate and detailed trajectories than the state-of-the-art TI methods using or not using RNA velocity information.

Integrating single-cell and spatial transcriptomics to elucidate intercellular tissue dynamics

Single-cell RNA sequencing (scRNA-seq) identifies cell subpopulations within tissue but does not capture their spatial distribution nor reveal local networks of intercellular communication acting in situ. A suite of recently developed techniques that localize RNA within tissue, including multiplexed in situ hybridization and in situ sequencing (here defined as high-plex RNA imaging) and spatial barcoding, can help address this issue. However, no method currently provides as complete a scope of the transcriptome as does scRNA-seq, underscoring the need for approaches to integrate single-cell and spatial data. Here, we review efforts to integrate scRNA-seq with spatial transcriptomics, including emerging integrative computational methods, and propose ways to effectively combine current methodologies.

BRIE2: computational identification of splicing phenotypes from single-cell transcriptomic experiments

RNA splicing is an important driver of heterogeneity in single cells through the expression of alternative transcripts and as a determinant of transcriptional kinetics. However, the intrinsic coverage limitations of scRNA-seq technologies make it challenging to associate specific splicing events to cell-level phenotypes. BRIE2 is a scalable computational method that resolves these issues by regressing single-cell transcriptomic data against cell-level features. We show that BRIE2 effectively identifies differential disease-associated alternative splicing events and allows a principled selection of genes that capture heterogeneity in transcriptional kinetics and improve RNA velocity analyses, enabling the identification of splicing phenotypes associated with biological changes.

Dynamics of transcriptional and post-transcriptional regulation

Despite gene expression programs being notoriously complex, RNA abundance is usually assumed as a proxy for transcriptional activity. Recently developed approaches, able to disentangle transcriptional and post-transcriptional regulatory processes, have revealed a more complex scenario. It is now possible to work out how synthesis, processing and degradation kinetic rates collectively determine the abundance of each gene's RNA. It has become clear that the same transcriptional output can correspond to different combinations of the kinetic rates. This underscores the fact that markedly different modes of gene expression regulation exist, each with profound effects on a gene's ability to modulate its own expression. This review describes the development of the experimental and computational approaches, including RNA metabolic labeling and mathematical modeling, that have been disclosing the mechanisms underlying complex transcriptional programs. Current limitations and future perspectives in the field are also discussed.

A comparison of metabolic labeling and statistical methods to infer genome-wide dynamics of RNA turnover

Metabolic labeling of newly transcribed RNAs coupled with RNA-seq is being increasingly used for genome-wide analysis of RNA dynamics. Methods including standard biochemical enrichment and recent nucleotide conversion protocols each require special experimental and computational treatment. Despite their immediate relevance, these technologies have not yet been assessed and benchmarked, and no data are currently available to advance reproducible research and the development of better inference tools. Here, we present a systematic evaluation and comparison of four RNA labeling protocols: 4sU-tagging biochemical enrichment, including spike-in RNA controls, SLAM-seq, TimeLapse-seq and TUC-seq. All protocols are evaluated based on practical considerations, conversion efficiency and wet lab requirements to handle hazardous substances. We also compute decay rate estimates and confidence intervals for each protocol using two alternative statistical frameworks, pulseR and GRAND-SLAM, for over 11 600 human genes and evaluate the underlying computational workflows for their robustness and ease of use. Overall, we demonstrate a high inter-method reliability across eight use case scenarios. Our results and data will facilitate reproducible research and serve as a resource contributing to a fuller understanding of RNA biology.

Spatio-temporal mRNA tracking in the early zebrafish embryo

Early stages of embryogenesis depend on subcellular localization and transport of maternal mRNA. However, systematic analysis of these processes is hindered by a lack of spatio-temporal information in single-cell RNA sequencing. Here, we combine spatially-resolved transcriptomics and single-cell RNA labeling to perform a spatio-temporal analysis of the transcriptome during early zebrafish development. We measure spatial localization of mRNA molecules within the one-cell stage embryo, which allows us to identify a class of mRNAs that are specifically localized at an extraembryonic position, the vegetal pole. Furthermore, we establish a method for high-throughput single-cell RNA labeling in early zebrafish embryos, which enables us to follow the fate of individual maternal transcripts until gastrulation. This approach reveals that many localized transcripts are specifically transported to the primordial germ cells. Finally, we acquire spatial transcriptomes of two xenopus species and compare evolutionary conservation of localized genes as well as enriched sequence motifs.

Single-Cell Sequencing Methodologies: From Transcriptome to Multi-Dimensional Measurement

Cells are the basic building blocks of biological systems, with inherent unique molecular features and development trajectories. The study of single cells facilitates in-depth understanding of cellular diversity, disease processes, and organization of multicellular organisms. Single-cell RNA sequencing (scRNA-seq) technologies have become essential tools for the interrogation of gene expression patterns and the dynamics of single cells, allowing cellular heterogeneity to be dissected at unprecedented resolution. Nevertheless, measuring at only transcriptome level or 1D is incomplete; the cellular heterogeneity reflects in multiple dimensions, including the genome, epigenome, transcriptome, spatial, and even temporal dimensions. Hence, integrative single cell analysis is highly desired. In addition, the way to interpret sequencing data by virtue of bioinformatic tools also exerts critical roles in revealing differential gene expression. Here, a comprehensive review that summarizes the cutting-edge single-cell transcriptome sequencing methodologies, including scRNA-seq, spatial and temporal transcriptome profiling, multi-omics sequencing and computational methods developed for scRNA-seq data analysis is provided. Finally, the challenges and perspectives of this field are discussed.

Single-cell technologies and analyses in hematopoiesis and hematological malignancies

In recent years, single-cell technologies have emerged as breakthrough techniques that enable the characterization of hematopoietic cell populations of normal and malignant tissue samples and will be combined in the near future with bulk technologies, currently used in clinical practice, to improve diagnosis, prognosis, and the search for novel molecular targets. These single-cell methods have the advantage of not masking cell-to-cell variation features and involve the study of genetic, epigenetic, transcriptional, and proteomic landscapes from a single-cell perspective. Latest advances in this field have enabled the development of novel strategies that significantly increase both sensitivity and high throughput. In this review, we emphasize emerging techniques aimed at assessing individual or multiomic parameters at single-cell resolution and analyze how these technologies have helped us understand hematopoietic variability and identify unknown and/or rare subpopulations. We also summarize the impact of these single-cell profiling strategies on the characterization of cell diversity within the tumor and the clonal evolution of multiple hematological malignancies in samples from untreated and treated patients, which provide valuable information for diagnosis, prognosis, and future treatments and explain why current therapies may fail. However, despite these improvements, new challenges lie ahead.

Fast-acting chemical tools to delineate causality in transcriptional control

Multi-dimensional omics profiling continues to illuminate the complexity of cellular processes. Because of difficult mechanistic interpretation of phenotypes induced by slow perturbation, fast experimental setups are increasingly used to dissect causal interactions directly in cells. Here we review a growing body of studies that leverage rapid pharmacological perturbation to delineate causality in gene control. When coupled with kinetically matched readouts, fast chemical genetic tools allow recording of primary phenotypes before confounding secondary effects manifest. The toolbox encompasses directly acting probes, such as active-site inhibitors and proteolysis-targeting chimeras, as well as strategies using genetic engineering to render target proteins chemically tractable, such as analog-sensitive and degron systems. We anticipate that extrapolation of these concepts to single-cell setups will further transform our mechanistic understanding of transcriptional control in the future. Importantly, the concept of leveraging speed to derive causality should be broadly applicable to many aspects of biological regulation.

Dive into Single, Seek Out Multiple: Probing Cancer Metastases via Single-Cell Sequencing and Imaging Techniques

Metastasis is the cause of most cancer deaths and continues to be the biggest challenge in clinical practice and laboratory investigation. The challenge is largely due to the intrinsic heterogeneity of primary and metastatic tumor populations and the complex interactions among cancer cells and cells in the tumor microenvironment. Therefore, it is important to determine the genotype and phenotype of individual cells so that the metastasis-driving events can be precisely identified, understood, and targeted in future therapies. Single-cell sequencing techniques have allowed the direct comparison of the genomic and transcriptomic changes among different stages of metastatic samples. Single-cell imaging approaches have enabled the live visualization of the heterogeneous behaviors of malignant and non-malignant cells in the tumor microenvironment. By applying these technologies, we are achieving a spatiotemporal precision understanding of cancer metastases and clinical therapeutic translations.

Iterative epigenomic analyses in the same single cell

Gene expression in individual cells is epigenetically regulated by DNA modifications, histone modifications, transcription factors, and other DNA-binding proteins. It has been shown that multiple histone modifications can predict gene expression and reflect future responses of bulk cells to extracellular cues. However, the predictive ability of epigenomic analysis is still limited for mechanistic research at a single cell level. To overcome this limitation, it would be useful to acquire reliable signals from multiple epigenetic marks in the same single cell. Here, we propose a new approach and a new method for analysis of several components of the epigenome in the same single cell. The new method allows reanalysis of the same single cell. We found that reanalysis of the same single cell is feasible, provides confirmation of the epigenetic signals, and allows application of statistical analysis to identify reproduced reads using data sets generated only from the single cell. Reanalysis of the same single cell is also useful to acquire multiple epigenetic marks from the same single cells. The method can acquire at least five epigenetic marks: H3K27ac, H3K27me3, mediator complex subunit 1, a DNA modification, and a DNA-interacting protein. We can predict active signaling pathways in K562 single cells using the epigenetic data and confirm that the predicted results strongly correlate with actual active signaling pathways identified by RNA-seq results. These results suggest that the new method provides mechanistic insights for cellular phenotypes through multilayered epigenome analysis in the same single cells.

Decoding the function of an oncogenic transcription factor: finding the first responders

Transcription factors (TFs) are frequently altered in human diseases. Identifying the direct and immediate target genes of TFs is critical to understanding their role in pathophysiology. Stengel et al. (2020) applied chemogenetic and nascent transcriptome mapping technologies to define the core gene set regulated by AML1-ETO-an oncogenic TF fusion protein frequently found in acute myeloid leukemia (AML).

Fast and furious: insights of back splicing regulation during nascent RNA synthesis

Alternative splicing of eukaryotic precursor (messenger) RNAs in the nucleus not only increases transcriptomic complexity, but also expands proteomic and functional diversity. In addition to basic types of alternative splicing, recent transcriptome-wide analyses have also suggested other new types of non-canonical splicing, such as back splicing and recursive splicing, and their widespread expression across species Increasing lines of evidence have suggested mechanisms for back splicing, including insights from analyses of nascent RNA sequencing. In this review, we discuss our current understanding of back splicing regulation, and highlight its distinct characteristics in processing during nascent RNA synthesis by taking advantage of metabolic tagging nascent RNA sequencing. Features of recursive splicing are also discussed in the perspective of nascent RNA sequencing.