Modeling collective cell behavior in cancer: Perspectives from an interdisciplinary conversation.

Collective cell behavior contributes to all stages of cancer progression. Understanding how collective behavior emerges through cell-cell interactions and decision-making will advance our understanding of cancer biology and provide new therapeutic approaches. Here, we summarize an interdisciplinary discussion on multicellular behavior in cancer, draw lessons from other scientific disciplines, and identify future directions.

Systems epidemiology and cancer: A review of the National Institutes of Health extramural grant portfolio 2013-2018.

OBJECTIVES: Systems epidemiology approaches may lead to a better understanding of the complex and dynamic multi-level constellation of contributors to cancer risk and outcomes and help target interventions. This grant portfolio analysis aimed to describe the National Institutes of Health (NIH) and the National Cancer Institute (NCI) investments in systems epidemiology and to identify gaps in the cancer systems epidemiology portfolio. METHODS: The analysis examined grants funded (2013-2018) through seven NIH systems science Funding Opportunity Announcements (FOAs) as well as cancer-specific systems epidemiology grants funded by NCI during that same time. Study characteristics were extracted from the grant abstracts and specific aims and coded. RESULTS: Of the 137 grants awarded under the NIH FOAs, 52 (38%) included systems epidemiology. Only five (4%) were focused on cancer systems epidemiology. The NCI-wide search (N = 453 grants) identified 35 grants (8%) that included cancer systems epidemiology in their specific aims. Most of these grants examined epidemiology and surveillance-based questions (60%); fewer addressed clinical care or clinical trials (37%). Fifty-four percent looked at multiple scales within the individual (e.g., cell, tissue, organ), 49% looked beyond the individual (e.g., individual, community, population), and few (9%) included both. Across all grants examined, the systems epidemiology grants primarily focused on discovery or prediction, rather than on impacts of intervention or policy. CONCLUSIONS: The most notable finding was that grants focused on cancer versus other diseases reflected a small percentage of the portfolio, highlighting the need to encourage more cancer systems epidemiology research. Opportunities include encouraging more multiscale research and continuing the support for broad examination of domains in these studies. Finally, the nascent discipline of systems epidemiology could benefit from the creation of standard terminology and definitions to guide future progress.

A lineage-resolved molecular atlas of C. elegans embryogenesis at single-cell resolution.

Caenorhabditis elegans is an animal with few cells but a wide diversity of cell types. In this study, we characterize the molecular basis for their specification by profiling the transcriptomes of 86,024 single embryonic cells. We identify 502 terminal and preterminal cell types, mapping most single-cell transcriptomes to their exact position in C. elegans' invariant lineage. Using these annotations, we find that (i) the correlation between a cell's lineage and its transcriptome increases from middle to late gastrulation, then falls substantially as cells in the nervous system and pharynx adopt their terminal fates; (ii) multilineage priming contributes to the differentiation of sister cells at dozens of lineage branches; and (iii) most distinct lineages that produce the same anatomical cell type converge to a homogenous transcriptomic state.

Gene expression distribution deconvolution in single-cell RNA sequencing.

Single-cell RNA sequencing (scRNA-seq) enables the quantification of each gene's expression distribution across cells, thus allowing the assessment of the dispersion, nonzero fraction, and other aspects of its distribution beyond the mean. These statistical characterizations of the gene expression distribution are critical for understanding expression variation and for selecting marker genes for population heterogeneity. However, scRNA-seq data are noisy, with each cell typically sequenced at low coverage, thus making it difficult to infer properties of the gene expression distribution from raw counts. Based on a reexamination of nine public datasets, we propose a simple technical noise model for scRNA-seq data with unique molecular identifiers (UMI). We develop deconvolution of single-cell expression distribution (DESCEND), a method that deconvolves the true cross-cell gene expression distribution from observed scRNA-seq counts, leading to improved estimates of properties of the distribution such as dispersion and nonzero fraction. DESCEND can adjust for cell-level covariates such as cell size, cell cycle, and batch effects. DESCEND's noise model and estimation accuracy are further evaluated through comparisons to RNA FISH data, through data splitting and simulations and through its effectiveness in removing known batch effects. We demonstrate how DESCEND can clarify and improve downstream analyses such as finding differentially expressed genes, identifying cell types, and selecting differentiation markers.

SAVER: gene expression recovery for single-cell RNA sequencing.

In single-cell RNA sequencing (scRNA-seq) studies, only a small fraction of the transcripts present in each cell are sequenced. This leads to unreliable quantification of genes with low or moderate expression, which hinders downstream analysis. To address this challenge, we developed SAVER (single-cell analysis via expression recovery), an expression recovery method for unique molecule index (UMI)-based scRNA-seq data that borrows information across genes and cells to provide accurate expression estimates for all genes.

Rare Cell Detection by Single-Cell RNA Sequencing as Guided by Single-Molecule RNA FISH.

Although single-cell RNA sequencing can reliably detect large-scale transcriptional programs, it is unclear whether it accurately captures the behavior of individual genes, especially those that express only in rare cells. Here, we use single-molecule RNA fluorescence in situ hybridization as a gold standard to assess trade-offs in single-cell RNA-sequencing data for detecting rare cell expression variability. We quantified the gene expression distribution for 26 genes that range from ubiquitous to rarely expressed and found that the correspondence between estimates across platforms improved with both transcriptome coverage and increased number of cells analyzed. Further, by characterizing the trade-off between transcriptome coverage and number of cells analyzed, we show that when the number of genes required to answer a given biological question is small, then greater transcriptome coverage is more important than analyzing large numbers of cells. More generally, our report provides guidelines for selecting quality thresholds for single-cell RNA-sequencing experiments aimed at rare cell analyses.

PIVOT: platform for interactive analysis and visualization of transcriptomics data.

BACKGROUND: Many R packages have been developed for transcriptome analysis but their use often requires familiarity with R and integrating results of different packages requires scripts to wrangle the datatypes. Furthermore, exploratory data analyses often generate multiple derived datasets such as data subsets or data transformations, which can be difficult to track. RESULTS: Here we present PIVOT, an R-based platform that wraps open source transcriptome analysis packages with a uniform user interface and graphical data management that allows non-programmers to interactively explore transcriptomics data. PIVOT supports more than 40 popular open source packages for transcriptome analysis and provides an extensive set of tools for statistical data manipulations. A graph-based visual interface is used to represent the links between derived datasets, allowing easy tracking of data versions. PIVOT further supports automatic report generation, publication-quality plots, and program/data state saving, such that all analysis can be saved, shared and reproduced. CONCLUSIONS: PIVOT will allow researchers with broad background to easily access sophisticated transcriptome analysis tools and interactively explore transcriptome datasets.

Primary Cell Culture of Live Neurosurgically Resected Aged Adult Human Brain Cells and Single Cell Transcriptomics.

Investigation of human CNS disease and drug effects has been hampered by the lack of a system that enables single-cell analysis of live adult patient brain cells. We developed a culturing system, based on a papain-aided procedure, for resected adult human brain tissue removed during neurosurgery. We performed single-cell transcriptomics on over 300 cells, permitting identification of oligodendrocytes, microglia, neurons, endothelial cells, and astrocytes after 3 weeks in culture. Using deep sequencing, we detected over 12,000 expressed genes, including hundreds of cell-type-enriched mRNAs, lncRNAs and pri-miRNAs. We describe cell-type- and patient-specific transcriptional hierarchies. Single-cell transcriptomics on cultured live adult patient derived cells is a prime example of the promise of personalized precision medicine. Because these cells derive from subjects ranging in age into their sixties, this system permits human aging studies previously possible only in rodent systems.

Assessing characteristics of RNA amplification methods for single cell RNA sequencing.

BACKGROUND: Recently, measurement of RNA at single cell resolution has yielded surprising insights. Methods for single-cell RNA sequencing (scRNA-seq) have received considerable attention, but the broad reliability of single cell methods and the factors governing their performance are still poorly known. RESULTS: Here, we conducted a large-scale control experiment to assess the transfer function of three scRNA-seq methods and factors modulating the function. All three methods detected greater than 70% of the expected number of genes and had a 50% probability of detecting genes with abundance greater than 2 to 4 molecules. Despite the small number of molecules, sequencing depth significantly affected gene detection. While biases in detection and quantification were qualitatively similar across methods, the degree of bias differed, consistent with differences in molecular protocol. Measurement reliability increased with expression level for all methods and we conservatively estimate measurements to be quantitative at an expression level greater than ~5-10 molecules. CONCLUSIONS: Based on these extensive control studies, we propose that RNA-seq of single cells has come of age, yielding quantitative biological information.

Single-cell transcriptomics and functional target validation of brown adipocytes show their complex roles in metabolic homeostasis.

Brown adipocytes (BAs) are specialized for adaptive thermogenesis and, upon sympathetic stimulation, activate mitochondrial uncoupling protein (UCP)-1 and oxidize fatty acids to generate heat. The capacity for brown adipose tissue (BAT) to protect against obesity and metabolic disease is recognized, yet information about which signals activate BA, besides ß3-adrenergic receptor stimulation, is limited. Using single-cell transcriptomics, we confirmed the presence of mRNAs encoding traditional BAT markers (i.e., UCP1, expressed in 100% of BAs Adrb3, expressed in <50% of BAs) in mouse and have shown single-cell variability (>1000-fold) in their expression at both the mRNA and protein levels. We further identified mRNAs encoding novel markers, orphan GPCRs, and many receptors that bind the classic neurotransmitters, neuropeptides, chemokines, cytokines, and hormones. The transcriptome variability between BAs suggests a much larger range of responsiveness of BAT than previously recognized and that not all BAs function identically. We examined the in vivo functional expression of 12 selected receptors by microinjecting agonists into live mouse BAT and analyzing the metabolic response. In this manner, we expanded the number of known receptors on BAs at least 25-fold, while showing that the expression of classic BA markers is more complex and variable than previously thought.

Variation is function: Are single cell differences functionally important?: Testing the hypothesis that single cell variation is required for aggregate function.

There is a growing appreciation of the extent of transcriptome variation across individual cells of the same cell type. While expression variation may be a byproduct of, for example, dynamic or homeostatic processes, here we consider whether single-cell molecular variation per se might be crucial for population-level function. Under this hypothesis, molecular variation indicates a diversity of hidden functional capacities within an ensemble of identical cells, and this functional diversity facilitates collective behavior that would be inaccessible to a homogenous population. In reviewing this topic, we explore possible functions that might be carried by a heterogeneous ensemble of cells; however, this question has proven difficult to test, both because methods to manipulate molecular variation are limited and because it is complicated to define, and measure, population-level function. We consider several possible methods to further pursue the hypothesis that variation is function through the use of comparative analysis and novel experimental techniques.

Deep sequencing reveals cell-type-specific patterns of single-cell transcriptome variation.

BACKGROUND: Differentiation of metazoan cells requires execution of different gene expression programs but recent single-cell transcriptome profiling has revealed considerable variation within cells of seeming identical phenotype. This brings into question the relationship between transcriptome states and cell phenotypes. Additionally, single-cell transcriptomics presents unique analysis challenges that need to be addressed to answer this question. RESULTS: We present high quality deep read-depth single-cell RNA sequencing for 91 cells from five mouse tissues and 18 cells from two rat tissues, along with 30 control samples of bulk RNA diluted to single-cell levels. We find that transcriptomes differ globally across tissues with regard to the number of genes expressed, the average expression patterns, and within-cell-type variation patterns. We develop methods to filter genes for reliable quantification and to calibrate biological variation. All cell types include genes with high variability in expression, in a tissue-specific manner. We also find evidence that single-cell variability of neuronal genes in mice is correlated with that in rats consistent with the hypothesis that levels of variation may be conserved. CONCLUSIONS: Single-cell RNA-sequencing data provide a unique view of transcriptome function; however, careful analysis is required in order to use single-cell RNA-sequencing measurements for this purpose. Technical variation must be considered in single-cell RNA-sequencing studies of expression variation. For a subset of genes, biological variability within each cell type appears to be regulated in order to perform dynamic functions, rather than solely molecular noise.

IVT-seq reveals extreme bias in RNA sequencing.

BACKGROUND: RNA-seq is a powerful technique for identifying and quantifying transcription and splicing events, both known and novel. However, given its recent development and the proliferation of library construction methods, understanding the bias it introduces is incomplete but critical to realizing its value. RESULTS: We present a method, in vitro transcription sequencing (IVT-seq), for identifying and assessing the technical biases in RNA-seq library generation and sequencing at scale. We created a pool of over 1,000 in vitro transcribed RNAs from a full-length human cDNA library and sequenced them with polyA and total RNA-seq, the most common protocols. Because each cDNA is full length, and we show in vitro transcription is incredibly processive, each base in each transcript should be equivalently represented. However, with common RNA-seq applications and platforms, we find 50% of transcripts have more than two-fold and 10% have more than 10-fold differences in within-transcript sequence coverage. We also find greater than 6% of transcripts have regions of dramatically unpredictable sequencing coverage between samples, confounding accurate determination of their expression. We use a combination of experimental and computational approaches to show rRNA depletion is responsible for the most significant variability in coverage, and several sequence determinants also strongly influence representation. CONCLUSIONS: These results show the utility of IVT-seq for promoting better understanding of bias introduced by RNA-seq. We find rRNA depletion is responsible for substantial, unappreciated biases in coverage introduced during library preparation. These biases suggest exon-level expression analysis may be inadvisable, and we recommend caution when interpreting RNA-seq results.

Transcriptome in vivo analysis (TIVA) of spatially defined single cells in live tissue.

Transcriptome profiling of single cells resident in their natural microenvironment depends upon RNA capture methods that are both noninvasive and spatially precise. We engineered a transcriptome in vivo analysis (TIVA) tag, which upon photoactivation enables mRNA capture from single cells in live tissue. Using the TIVA tag in combination with RNA sequencing (RNA-seq), we analyzed transcriptome variance among single neurons in culture and in mouse and human tissue in vivo. Our data showed that the tissue microenvironment shapes the transcriptomic landscape of individual cells. The TIVA methodology is, to our knowledge, the first noninvasive approach for capturing mRNA from live single cells in their natural microenvironment.

Serotonergic neuron regulation informed by in vivo single-cell transcriptomics.

Despite the recognized importance of the dorsal raphe (DR) serotonergic (5-HT) nuclei in the pathophysiology of depression and anxiety, the molecular components/putative drug targets expressed by these neurons are poorly characterized. Utilizing the promoter of an ETS domain transcription factor that is a stable marker of 5-HT neurons (Pet-1) to drive 5-HT neuronal expression of YFP, we identified 5-HT neurons in live acute slices. We isolated RNA from single 5-HT neurons in the ventromedial and lateral wings of the DR and performed single-cell RNA-Seq analysis identifying >500 G-protein coupled receptors (GPCRs) including receptors for classical transmitters, lipid signals, and peptides as well as dozens of orphan-GPCRs. Using these data to inform our selection of receptors to assess, we found that oxytocin and lysophosphatidic acid 1 receptors are translated and active in costimulating, with the α1-adrenergic receptor, the firing of DR 5-HT neurons, while the effects of histamine are inhibitory and exerted at H3 histamine receptors. The inhibitory histamine response provides evidence for tonic in vivo histamine inhibition of 5-HT neurons. This study illustrates that unbiased single-cell transcriptomics coupled with functional analyses provides novel insights into how neurons and neuronal systems are regulated.

Quantitative biology of single neurons.

The building blocks of complex biological systems are single cells. Fundamental insights gained from single-cell analysis promise to provide the framework for understanding normal biological systems development as well as the limits on systems/cellular ability to respond to disease. The interplay of cells to create functional systems is not well understood. Until recently, the study of single cells has concentrated primarily on morphological and physiological characterization. With the application of new highly sensitive molecular and genomic technologies, the quantitative biochemistry of single cells is now accessible.