RNA velocity of single cells

Gioele La Manno, Zehua Liu
2018
Journal paper

Abstract

RNA abundance is a powerful indicator of the state of individual cells. Single-cell RNA sequencing can reveal RNA abundance with high quantitative accuracy, sensitivity and throughput 1 . However, this approach captures only a static snapshot at a point in time, posing a challenge for the analysis of time-resolved phenomena such as embryogenesis or tissue regeneration. Here we show that RNA velocity—the time derivative of the gene expression state—can be directly estimated by distinguishing between unspliced and spliced mRNAs in common single-cell RNA sequencing protocols. RNA velocity is a high-dimensional vector that predicts the future state of individual cells on a timescale of hours. We validate its accuracy in the neural crest lineage, demonstrate its use on multiple published datasets and technical platforms, reveal the branching lineage tree of the developing mouse hippocampus, and examine the kinetics of transcription in human embryonic brain. We expect RNA velocity to greatly aid the analysis of developmental lineages and cellular dynamics, particularly in humans.

Official source

https://infoscience.epfl.ch/record/265419?ln=en

About this result

This page is automatically generated and may contain information that is not correct, complete, up-to-date, or relevant to your search query. The same applies to every other page on this website. Please make sure to verify the information with EPFL's official sources.

Related concepts

Related publications

Cells live in ever-changing environments, thereby facing a variety of dynamic environmental signals. Environmental stimuli elicit intracellular responses through signaling pathways, which converge on transcriptional activation or repression of target genes. Despite intensive research, dissecting the complex interactions between pathway components that modulate mRNA and protein production still remains a difficult task. To this end, single-cell approaches provide unique insights into intracellular processes and cell responses to environmental stimuli, otherwise inaccessible with traditional bulk studies. Single-cell measurements have revealed that isogenic cells sharing the same environment display substantial heterogeneity in both transcript and protein level. As the initial step in gene expression, transcription is an important source of such variability in eukaryotic cells due to low number of molecules involved, transcription factor dynamics and the discrete nature of biochemical reactions. Notably, production of mRNA during transcription occurs in short periods of activity, termed as transcriptional bursts, followed by longer periods of inactivity. Despite widespread observations of transcriptional dynamics in mammals, upstream molecular mechanisms shaping the eukaryotic transcription remained elusive. In this study, we quantitatively linked upstream factors and transcriptional kinetics by developing an experimental system to temporally monitor, simultaneously, inputs (transcription factors) and outputs (target gene expression) in mammalian cells, in response to stimulus. We established two Tet-On inducible stable cell lines each expressing a fusion protein of a luminescence (Nanoluciferase) reporter to either SMAD4 or SMAD2, the two main transcrip- tion factors downstream of the TGF-Î² pathway. These cell lines that also contain a short-lived firefly luciferase reporter for the expression of the target endogenous connective tissue growth factor (ctgf ) gene allowed us to quantitatively link nuclear accumulation of SMADs upon TGF-Î² stimulation to the target gene expression in real-time single cell measurements. Time- lapse luminescence microscopy and image analysis with the custom developed CAST (Cell Automated Segmentation and Tracking platform) platform provided quantitative single-cell data to link the upstream transcription factor profile to its target gene activity. The data revealed weak single-cell correlations between translocation level and the target gene expression response which suggests a mechanism in which the translocation of SMADs initiates the transcriptional response but affects the response amplitude minimally. However, SMAD4 expression levels influenced the target gene response dynamics such that cells with high SMAD4 abundance favored to respond in a more sustained and even oscillatory manner. This suggests a mechanism consisting of a dynamic interplay between the transcriptional activators and feedback mechanisms in TGF-Î² signaling. We explored further the effect of different factors such as ligand concentration, ligand type on both signaling and target gene response. Taken together, this study proposes an experimental and quantitative framework that dissect mechanisms underlying transcriptional response to stimulus.

EPFL2018

Biological oscillators are pervasive in biology, covering all aspects of life from enzyme kinetics reactions to population dynamics. Although their behaviour has been intensively studied in the last decades, the recent advances of high-throughput experimental technologies in the fields of omics and microscopy has called for the development of new analysis methods. Among the many types of models and quantitative analyses, parametric approaches are promising as they enable for a mechanistic or physical explanation of the phenomena under study. In particular, dynamical systems theory seems particularly adapted as the vast majority of oscillators can be modelled through differential equations. Dynamical systems parameters can also be easily opti-mized via maximum likelihood approaches. The validity of the inferred model can then be assessed from the quality of its predictions. We here present three different scientific questions regarding noisy biological oscillators, which are answered using maximum-likelihood inference approaches applied to parametric models. We first take interest in the characterisation of the influence of the cell-cycle over the circadi-an clock in individual mammalian cells. To this end, we develop a method combining a Hidden Markov Model with an Expectation-Maximization algorithm to infer their coupling from single-cell microscopy traces. We show that this coupling predicts multiple phase-locked states exhib-iting different degrees of robustness against molecular fluctuations inherent to cellular scale bio-logical oscillators. We then try to understand how the mammalian transcriptome behaves in the liver. Thence, we use single-cell RNA sequencing (scRNA-seq) along with mixed-models to investigate the interplay between gene regulation in space and time. Categorising mRNA expression profiles using mixed-effect models and smFISH validations, we find that many genes in the liver are both zonated and rhythmic, most of them showing multiplicative space-time effects. Finally, we look more closely at the cell-cycle, as it is one of the main drivers of gene expres-sion cell-to-cell heterogeneity in otherwise homogeneous cell populations. Here, we would like to understand if and how cell-cycle velocity changes depending on the phase of the cycling cells. To that end, we formulate the problem in terms of an autonomous dynamical system and use this to infer consistent dynamics for the cell-cycle from scRNA-seq data. Phase inference being paramount in all of these three studies, a short technical review on the topic is also provided at the end of this thesis, along with Julia scripts for the main inference methods presented. Various computational tools assisting the understanding of the scientific questions at stake are also presented, including a Python Dash web-app, a D3 widget and many Matplotlib animations and widgets.

EPFL2020

The chronology of events in time-space is naturally available to the senses, and the spatial and temporal dimensions of events entangle in episodic memory when navigating the real world. The mapping of time-space during navigation in both animals and humans implicates the hippocampal formation. Yet, one arguably unique human trait is the capacity to imagine mental chronologies that have not been experienced but may involve real events-the foundation of causal reasoning. Herein, we asked whether the hippocampal formation is involved in mental navigation in time (and space), which requires internal manipulations of events in time and space from an egocentric perspective. To address this question, we reanalyzed a magnetoencephalography data set collected while participants self-projected in time or in space and ordered historical events as occurring before/after or west/east of the mental self [Gauthier, B., Pestke, K., & van Wassenhove, V. Building the arrow of time horizontal ellipsis Over time: A sequence of brain activity mapping imagined events in time and space.Cerebral Cortex,29, 4398-4414, 2019]. Because of the limitations of source reconstruction algorithms in the previous study, the implication of hippocampus proper could not be explored. Here, we used a source reconstruction method accounting explicitly for the hippocampal volume to characterize the involvement of deep structures belonging to the hippocampal formation (bilateral hippocampi [hippocampi proper], entorhinal cortices, and parahippocampal cortex). We found selective involvement of the medial temporal lobes (MTLs) with a notable lateralization of the main effects: Whereas temporal ordinality engaged mostly the left MTL, spatial ordinality engaged mostly the right MTL. We discuss the possibility of a top-down control of activity in the human hippocampal formation during mental time (and space) travels.

MIT PRESS2020

EPFL2018

EPFL2020