scSLAM-seq reveals core features of transcription dynamics in single cells

Current single-cell RNA sequencing approaches gives a snapshot of a cellular phenotype but convey no information on the temporal dynamics of transcription. Moreover, the stochastic nature of transcription at molecular level is not recovered. Here, we present single-cell SLAM-seq (scSLAM-seq), which integrates metabolic RNA labeling, biochemical nucleoside conversion and single-cell RNA-seq to directly measure total transcript levels and transcriptional activity by differentiating newly synthesized from pre-existing RNA for thousands of genes per single cell. scSLAM-seq recovers the earliest virus-induced changes in cytomegalovirus infection and reveals a so far hidden phase of viral gene expression comprising promiscuous transcription of all kinetic classes. It depicts the stochastic nature of transcription and demonstrates extensive gene-specific differences. These range from stable transcription rates to on-off dynamics which coincide with gene-/promoter-intrinsic features (Tbp-TATA-box interactions and DNA methylation). Gene but not cell-specific features thus explain the heterogeneity in transcriptomes between individual cells and the transcriptional response to perturbations.

[1]  Donald Macleod,et al.  A fraction of the mouse genome that is derived from islands of nonmethylated, CpG-rich DNA , 1985, Cell.

[2]  T. Shenk,et al.  A subset of viral transcripts packaged within human cytomegalovirus particles. , 2000, Science.

[3]  Vitaly Epshtein,et al.  Cooperation Between RNA Polymerase Molecules in Transcription Elongation , 2003, Science.

[4]  P. Khavari,et al.  Divergent gene regulation and growth effects by NF-κB in epithelial and mesenchymal cells of human skin , 2003, Oncogene.

[5]  T. Shenk,et al.  RNAs Are Packaged into Human Cytomegalovirus Virions in Proportion to Their Intracellular Concentration , 2004, Journal of Virology.

[6]  D. Brutlag,et al.  A genome-wide analysis of CpG dinucleotides in the human genome distinguishes two distinct classes of promoters , 2006, Proceedings of the National Academy of Sciences of the United States of America.

[7]  A. Oudenaarden,et al.  Nature, Nurture, or Chance: Stochastic Gene Expression and Its Consequences , 2008, Cell.

[8]  R. Zimmer,et al.  High-resolution gene expression profiling for simultaneous kinetic parameter analysis of RNA synthesis and decay. , 2008, RNA.

[9]  R. Kalejta Functions of human cytomegalovirus tegument proteins prior to immediate early gene expression. , 2008, Current topics in microbiology and immunology.

[10]  Israel Steinfeld,et al.  BMC Bioinformatics BioMed Central , 2008 .

[11]  Carlos Bustamante,et al.  Nucleosomal Fluctuations Govern the Transcription Dynamics of RNA Polymerase II , 2009, Science.

[12]  Achim Tresch,et al.  Dynamic transcriptome analysis measures rates of mRNA synthesis and decay in yeast , 2011, Molecular systems biology.

[13]  B. Roizman,et al.  HSV-1 gene expression from reactivated ganglia is disordered and concurrent with suppression of latency-associated transcript and miRNAs , 2011, Proceedings of the National Academy of Sciences.

[14]  Thomas Bonfert,et al.  Real-time Transcriptional Profiling of Cellular and Viral Gene Expression during Lytic Cytomegalovirus Infection , 2012, PLoS pathogens.

[15]  A. van Oudenaarden,et al.  Using Gene Expression Noise to Understand Gene Regulation , 2012, Science.

[16]  I. Mohr,et al.  Transient Reversal of Episome Silencing Precedes VP16-Dependent Transcription during Reactivation of Latent HSV-1 in Neurons , 2012, PLoS pathogens.

[17]  Thomas Höfer,et al.  Multi-layered stochasticity and paracrine signal propagation shape the type-I interferon response , 2012, Molecular systems biology.

[18]  Peter A. Jones Functions of DNA methylation: islands, start sites, gene bodies and beyond , 2012, Nature Reviews Genetics.

[19]  S. Kaufmann,et al.  Ultrashort and progressive 4sU-tagging reveals key characteristics of RNA processing at nucleotide resolution , 2012, Genome research.

[20]  Aleksandra A. Kolodziejczyk,et al.  Accounting for technical noise in single-cell RNA-seq experiments , 2013, Nature Methods.

[21]  Mikhail Pachkov,et al.  SwissRegulon, a database of genome-wide annotations of regulatory sites: recent updates , 2012, Nucleic Acids Res..

[22]  P. Kharchenko,et al.  Bayesian approach to single-cell differential expression analysis , 2014, Nature Methods.

[23]  Cole Trapnell,et al.  The dynamics and regulators of cell fate decisions are revealed by pseudotemporal ordering of single cells , 2014, Nature Biotechnology.

[24]  L. Dölken,et al.  Murine Cytomegalovirus Virion-Associated Protein M45 Mediates Rapid NF-κB Activation after Infection , 2014, Journal of Virology.

[25]  Björn Usadel,et al.  Trimmomatic: a flexible trimmer for Illumina sequence data , 2014, Bioinform..

[26]  Jason D. Buenrostro,et al.  Transient Acquisition of Pluripotency During Somatic Cell Transdifferentiation with iPSC Reprogramming Factors , 2015, Nature Biotechnology.

[27]  Allon M. Klein,et al.  Droplet Barcoding for Single-Cell Transcriptomics Applied to Embryonic Stem Cells , 2015, Cell.

[28]  Evan Z. Macosko,et al.  Highly Parallel Genome-wide Expression Profiling of Individual Cells Using Nanoliter Droplets , 2015, Cell.

[29]  Felix Naef,et al.  Structure of silent transcription intervals and noise characteristics of mammalian genes , 2015, Molecular systems biology.

[30]  R. Shamir,et al.  ElemeNT: a computational tool for detecting core promoter elements , 2015, Transcription.

[31]  Thomas M. Norman,et al.  Perturb-Seq: Dissecting Molecular Circuits with Scalable Single-Cell RNA Profiling of Pooled Genetic Screens , 2016, Cell.

[32]  Thomas M. Norman,et al.  A Multiplexed Single-Cell CRISPR Screening Platform Enables Systematic Dissection of the Unfolded Protein Response , 2016, Cell.

[33]  A. Regev,et al.  Revealing the vectors of cellular identity with single-cell genomics , 2016, Nature Biotechnology.

[34]  Joseph L. Herman,et al.  Characterizing transcriptional heterogeneity through pathway and gene set overdispersion analysis , 2015, Nature Methods.

[35]  C. Benedict,et al.  cGAS-STING Signaling Regulates Initial Innate Control of Cytomegalovirus Infection , 2016, Journal of Virology.

[36]  Christophe Zimmer,et al.  A single-molecule view of transcription reveals convoys of RNA polymerases and multi-scale bursting , 2016, Nature Communications.

[37]  Aafke A. van den Berg,et al.  Crowding-induced transcriptional bursts dictate polymerase and nucleosome density profiles along genes , 2016, bioRxiv.

[38]  Johannes Zuber,et al.  Thiol-linked alkylation of RNA to assess expression dynamics , 2017, Nature Methods.

[39]  Jesse J. Lipp,et al.  SLAM-seq defines direct gene-regulatory functions of the BRD4-MYC axis , 2018, Science.

[40]  Yi Zhang,et al.  Two-phase differential expression analysis for single cell RNA-seq , 2018, Bioinform..

[41]  Paul Hoffman,et al.  Integrating single-cell transcriptomic data across different conditions, technologies, and species , 2018, Nature Biotechnology.

[42]  Meaghan C. Sullivan,et al.  TimeLapse-seq: Adding a temporal dimension to RNA sequencing through nucleoside recoding , 2018, Nature Methods.

[43]  R. Satija,et al.  Single-cell RNA sequencing to explore immune cell heterogeneity , 2017, Nature Reviews Immunology.

[44]  Christopher J. Cronin,et al.  Dynamics and Spatial Genomics of the Nascent Transcriptome by Intron seqFISH , 2018, Cell.

[45]  Tim De Meyer,et al.  Analysis of DNA methylation in cancer: location revisited , 2018, Nature Reviews Clinical Oncology.

[46]  Florian Erhard,et al.  Dissecting newly transcribed and old RNA using GRAND-SLAM , 2018, Bioinform..

[47]  Andrew E. Teschendorff,et al.  Statistical and integrative system-level analysis of DNA methylation data , 2017, Nature Reviews Genetics.

[48]  Extended data : Figures , 2020 .