Deciphering molecular circuits from genetic variation underlying transcriptional responsiveness to stimuli

Individual genetic variation affects gene expression in response to stimuli, often by influencing complex molecular circuits. Here we combine genomic and intermediate-scale transcriptional profiling with computational methods to identify variants that affect the responsiveness of genes to stimuli (responsiveness QTLs; reQTLs) and to position these variants in molecular circuit diagrams. We apply this approach to study variation in transcriptional responsiveness to pathogen components in dendritic cells from recombinant inbred mouse strains. We identify reQTLs that correlate with particular stimuli and position them in known pathways. For example, in response to a virus-like stimulus, a trans-acting variant acts as an activator of the antiviral response; using RNAi, we identify Rgs16 as the likely causal gene. Our approach charts an experimental and analytic path to decipher the mechanisms underlying genetic variation in circuits that control responses to stimuli

High-Resolution Sequencing and Modeling Identifies Distinct Dynamic RNA Regulatory Strategies

Cells control dynamic transitions in transcript levels by regulating transcription, processing, and/or degradation through an integrated regulatory strategy. Here, we combine RNA metabolic labeling, rRNA-depleted RNA-seq, and DRiLL, a novel computational framework, to quantify the level; editing sites; and transcription, processing, and degradation rates of each transcript at a splice junction resolution during the LPS response of mouse dendritic cells. Four key regulatory strategies, dominated by RNA transcription changes, generate most temporal gene expression patterns. Noncanonical strategies that also employ dynamic posttranscriptional regulation control only a minority of genes, but provide unique signal processing features. We validate Tristetraprolin (TTP) as a major regulator of RNA degradation in one noncanonical strategy. Applying DRiLL to the regulation of noncoding RNAs and to zebrafish embryogenesis demonstrates its broad utility. Our study provides a new quantitative approach to discover transcriptional and posttranscriptional events that control dynamic changes in transcript levels using RNA sequencing data.National Human Genome Research Institute (U.S.) (Centers for Excellence in Genomics Science 1P50HG006193-01)Howard Hughes Medical InstituteNational Institutes of Health (U.S.) (Pioneer Award)Massachusetts Institute of Technology. William Asbjornsen Albert Memorial FellowshipXerox Fellowship Progra

Perturb-Seq: Dissecting Molecular Circuits with Scalable Single-Cell RNA Profiling of Pooled Genetic Screens

Genetic screens help infer gene function in mammalian cells, but it has remained difficult to assay complex phenotypes—such as transcriptional profiles—at scale. Here, we develop Perturb-seq, combining single-cell RNA sequencing (RNA-seq) and clustered regularly interspaced short palindromic repeats (CRISPR)-based perturbations to perform many such assays in a pool. We demonstrate Perturb-seq by analyzing 200,000 cells in immune cells and cell lines, focusing on transcription factors regulating the response of dendritic cells to lipopolysaccharide (LPS). Perturb-seq accurately identifies individual gene targets, gene signatures, and cell states affected by individual perturbations and their genetic interactions. We posit new functions for regulators of differentiation, the anti-viral response, and mitochondrial function during immune activation. By decomposing many high content measurements into the effects of perturbations, their interactions, and diverse cell metadata, Perturb-seq dramatically increases the scope of pooled genomic assays. Keywords: single-cell RNA-seq; pooled screen; CRISPR; epistasis; genetic interaction

Semiconductor-based DNA sequencing of histone modification states

The recent development of a semiconductor-based, non-optical DNA sequencing technology promises scalable, low-cost and rapid sequence data production. The technology has previously been applied mainly to genomic sequencing and targeted re-sequencing. Here we demonstrate the utility of Ion Torrent semiconductor-based sequencing for sensitive, efficient and rapid chromatin immunoprecipitation followed by sequencing (ChIP-seq) through the application of sample preparation methods that are optimized for ChIP-seq on the Ion Torrent platform. We leverage this method for epigenetic profiling of tumour tissues

Single-cell transcriptomics reveals bimodality in expression and splicing in immune cells

Recent molecular studies have shown that, even when derived from a seemingly homogenous population, individual cells can exhibit substantial differences in gene expression, protein levels and phenotypic output1, 2, 3, 4, 5, with important functional consequences4, 5. Existing studies of cellular heterogeneity, however, have typically measured only a few pre-selected RNAs1, 2 or proteins5, 6 simultaneously, because genomic profiling methods3 could not be applied to single cells until very recently7, 8, 9, 10. Here we use single-cell RNA sequencing to investigate heterogeneity in the response of mouse bone-marrow-derived dendritic cells (BMDCs) to lipopolysaccharide. We find extensive, and previously unobserved, bimodal variation in messenger RNA abundance and splicing patterns, which we validate by RNA-fluorescence in situ hybridization for select transcripts. In particular, hundreds of key immune genes are bimodally expressed across cells, surprisingly even for genes that are very highly expressed at the population average. Moreover, splicing patterns demonstrate previously unobserved levels of heterogeneity between cells. Some of the observed bimodality can be attributed to closely related, yet distinct, known maturity states of BMDCs; other portions reflect differences in the usage of key regulatory circuits. For example, we identify a module of 137 highly variable, yet co-regulated, antiviral response genes. Using cells from knockout mice, we show that variability in this module may be propagated through an interferon feedback circuit, involving the transcriptional regulators Stat2 and Irf7. Our study demonstrates the power and promise of single-cell genomics in uncovering functional diversity between cells and in deciphering cell states and circuits.National Institutes of Health (U.S.) (NIH Postdoctoral Fellowship (1F32HD075541-01))Charles H. Hood Foundation (Postdoctoral Fellowship)National Institutes of Health (U.S.) (NIH grant U54 AI057159)National Institutes of Health (U.S.) (NIH New Innovator Award (DP2 OD002230))National Institutes of Health (U.S.) (NIH CEGS Award (1P50HG006193-01))National Institutes of Health (U.S.) (NIH Pioneer Awards (5DP1OD003893-03))National Institutes of Health (U.S.) (NIH Pioneer Awards (DP1OD003958-01))Broad Institute of MIT and HarvardBroad Institute of MIT and Harvard (Klarman Cell Observatory

Single-cell RNA-seq reveals dynamic paracrine control of cellular variation

High-throughput single-cell transcriptomics offers an unbiased approach for understanding the extent, basis and function of gene expression variation between seemingly identical cells. Here we sequence single-cell RNA-seq libraries prepared from over 1,700 primary mouse bone-marrow-derived dendritic cells spanning several experimental conditions. We find substantial variation between identically stimulated dendritic cells, in both the fraction of cells detectably expressing a given messenger RNA and the transcript’s level within expressing cells. Distinct gene modules are characterized by different temporal heterogeneity profiles. In particular, a ‘core’ module of antiviral genes is expressed very early by a few ‘precocious’ cells in response to uniform stimulation with a pathogenic component, but is later activated in all cells. By stimulating cells individually in sealed microfluidic chambers, analysing dendritic cells from knockout mice, and modulating secretion and extracellular signalling, we show that this response is coordinated by interferon-mediated paracrine signalling from these precocious cells. Notably, preventing cell-to-cell communication also substantially reduces variability between cells in the expression of an early-induced ‘peaked’ inflammatory module, suggesting that paracrine signalling additionally represses part of the inflammatory program. Our study highlights the importance of cell-to-cell communication in controlling cellular heterogeneity and reveals general strategies that multicellular populations can use to establish complex dynamic responses.National Human Genome Research Institute (U.S.). Centers of Excellence in Genomic Science (1P50HG006193-01)National Institutes of Health (U.S.). Pioneer Award (DP1OD003958-01)Howard Hughes Medical InstituteBroad Institute of MIT and Harvard. Klarman Cell Observator

Single cell RNA Seq reveals dynamic paracrine control of cellular variation

High-throughput single-cell transcriptomics offers an unbiased approach for understanding the extent, basis, and function of gene expression variation between seemingly identical cells. Here, we sequence single-cell RNA-Seq libraries prepared from over 1,700 primary mouse bone marrow derived dendritic cells (DCs) spanning several experimental conditions. We find substantial variation between identically stimulated DCs, in both the fraction of cells detectably expressing a given mRNA and the transcript’s level within expressing cells. Distinct gene modules are characterized by different temporal heterogeneity profiles. In particular, a “core” module of antiviral genes is expressed very early by a few “precocious” cells, but is later activated in all cells. By stimulating cells individually in sealed microfluidic chambers, analyzing DCs from knockout mice, and modulating secretion and extracellular signaling, we show that this response is coordinated via interferon-mediated paracrine signaling. Surprisingly, preventing cell-to-cell communication also substantially reduces variability in the expression of an early-induced “peaked” inflammatory module, suggesting that paracrine signaling additionally represses part of the inflammatory program. Our study highlights the importance of cell-to-cell communication in controlling cellular heterogeneity and reveals general strategies that multicellular populations use to establish complex dynamic responses

Metabolic labeling of RNA uncovers principles of RNA production and degradation dynamics in mammalian cells

available in PMC 2011 November 01.Cellular RNA levels are determined by the interplay of RNA production, processing and degradation. However, because most studies of RNA regulation do not distinguish the separate contributions of these processes, little is known about how they are temporally integrated. Here we combine metabolic labeling of RNA at high temporal resolution with advanced RNA quantification and computational modeling to estimate RNA transcription and degradation rates during the response of mouse dendritic cells to lipopolysaccharide. We find that changes in transcription rates determine the majority of temporal changes in RNA levels, but that changes in degradation rates are important for shaping sharp 'peaked' responses. We used sequencing of the newly transcribed RNA population to estimate temporally constant RNA processing and degradation rates genome wide. Degradation rates vary significantly between genes and contribute to the observed differences in the dynamic response. Certain transcripts, including those encoding cytokines and transcription factors, mature faster. Our study provides a quantitative approach to study the integrative process of RNA regulation.Human Frontier Science Program (Strasbourg, France)Howard Hughes Medical InstituteBurroughs Wellcome Fund (Career Award at the Scientific Interface

A single-cell survey of the small intestinal epithelium

Intestinal epithelial cells (IECs) absorb nutrients, respond to microbes, provide barrier function and help coordinate immune responses. We profiled 53,193 individual epithelial cells from mouse small intestine and organoids, and characterized novel subtypes and their gene signatures. We showed unexpected diversity of hormone-secreting enteroendocrine cells and constructed their novel taxonomy. We distinguished between two tuft cell subtypes, one of which expresses the epithelial cytokine TSLP and CD45 (Ptprc), the pan-immune marker not previously associated with non-hematopoietic cells. We also characterized how cell-intrinsic states and cell proportions respond to bacterial and helminth infections. Salmonella infection caused an increase in Paneth cells and enterocytes abundance, and broad activation of an antimicrobial program. In contrast, Heligmosomoides polygyrus caused an expansion of goblet and tuft cell populations. Our survey highlights new markers and programs, associates sensory molecules to cell types, and uncovers principles of gut homeostasis and response to pathogens