Cancer-associated fibroblasts are the main contributors to epithelial-to-mesenchymal signatures in the tumor microenvironment.

Epithelial-to-mesenchymal transition (EMT) is associated with tumor initiation, metastasis, and drug resistance. However, the mechanisms underlying these associations are largely unknown. We studied several tumor types to identify the source of EMT gene expression signals and a potential mechanism of resistance to immuno-oncology treatment. Across tumor types, EMT-related gene expression was strongly associated with expression of stroma-related genes. Based on RNA sequencing of multiple patient-derived xenograft models, EMT-related gene expression was enriched in the stroma versus parenchyma. EMT-related markers were predominantly expressed by cancer-associated fibroblasts (CAFs), cells of mesenchymal origin which produce a variety of matrix proteins and growth factors. Scores derived from a 3-gene CAF transcriptional signature (COL1A1, COL1A2, COL3A1) were sufficient to reproduce association between EMT-related markers and disease prognosis. Our results suggest that CAFs are the primary source of EMT signaling and have potential roles as biomarkers and targets for immuno-oncology therapies.

Single-Cell RNA Sequencing Reveals Molecular Features of Heterogeneity in the Murine Retinal Pigment Epithelium.

Transcriptomic analysis of the mammalian retinal pigment epithelium (RPE) aims to identify cellular networks that influence ocular development, maintenance, function, and disease. However, available evidence points to RPE cell heterogeneity within native tissue, which adds complexity to global transcriptomic analysis. Here, to assess cell heterogeneity, we performed single-cell RNA sequencing of RPE cells from two young adult male C57BL/6J mice. Following quality control to ensure robust transcript identification limited to cell singlets, we detected 13,858 transcripts among 2667 and 2846 RPE cells. Dimensional reduction by principal component analysis and uniform manifold approximation and projection revealed six distinct cell populations. All clusters expressed transcripts typical of RPE cells; the smallest (C1, containing 1-2% of total cells) exhibited the hallmarks of stem and/or progenitor (SP) cells. Placing C1-6 along a pseudotime axis suggested a relative decrease in melanogenesis and SP gene expression and a corresponding increase in visual cycle gene expression upon RPE maturation. K-means clustering of all detected transcripts identified additional expression patterns that may advance the understanding of RPE SP cell maintenance and the evolution of cellular metabolic networks during development. This work provides new insights into the transcriptome of the mouse RPE and a baseline for identifying experimentally induced transcriptional changes in future studies of this tissue.

Cellular taxonomy and spatial organization of the murine ventral posterior hypothalamus.

The ventral posterior hypothalamus (VPH) is an anatomically complex brain region implicated in arousal, reproduction, energy balance, and memory processing. However, neuronal cell type diversity within the VPH is poorly understood, an impediment to deconstructing the roles of distinct VPH circuits in physiology and behavior. To address this question, we employed a droplet-based single-cell RNA sequencing (scRNA-seq) approach to systematically classify molecularly distinct cell populations in the mouse VPH. Analysis of >16,000 single cells revealed 20 neuronal and 18 non-neuronal cell populations, defined by suites of discriminatory markers. We validated differentially expressed genes in selected neuronal populations through fluorescence in situ hybridization (FISH). Focusing on the mammillary bodies (MB), we discovered transcriptionally-distinct clusters that exhibit neuroanatomical parcellation within MB subdivisions and topographic projections to the thalamus. This single-cell transcriptomic atlas of VPH cell types provides a resource for interrogating the circuit-level mechanisms underlying the diverse functions of VPH circuits.

Mapping systemic lupus erythematosus heterogeneity at the single-cell level.

Patients with systemic lupus erythematosus (SLE) display a complex blood transcriptome whose cellular origin is poorly resolved. Using single-cell RNA sequencing, we profiled ~276,000 peripheral blood mononuclear cells from 33 children with SLE with different degrees of disease activity and 11 matched controls. Increased expression of interferon-stimulated genes (ISGs) distinguished cells from children with SLE from healthy control cells. The high ISG expression signature (ISGhi) derived from a small number of transcriptionally defined subpopulations within major cell types, including monocytes, CD4+ and CD8+ T cells, natural killer cells, conventional and plasmacytoid dendritic cells, B cells and especially plasma cells. Expansion of unique subpopulations enriched in ISGs and/or in monogenic lupus-associated genes classified patients with the highest disease activity. Profiling of ~82,000 single peripheral blood mononuclear cells from adults with SLE confirmed the expansion of similar subpopulations in patients with the highest disease activity. This study lays the groundwork for resolving the origin of the SLE transcriptional signatures and the disease heterogeneity towards precision medicine applications.

Single-cell analyses reveal increased intratumoral heterogeneity after the onset of therapy resistance in small-cell lung cancer.

The natural history of small cell lung cancer (SCLC) includes rapid evolution from chemosensitivity to chemoresistance, although mechanisms underlying this evolution remain obscure due to scarcity of post-relapse tissue samples. We generated circulating tumor cell (CTC)-derived xenografts (CDXs) from SCLC patients to study intratumoral heterogeneity (ITH) via single-cell RNAseq of chemo-sensitive and -resistant CDXs and patient CTCs. We found globally increased ITH including heterogeneous expression of therapeutic targets and potential resistance pathways, such as EMT, between cellular subpopulations following treatment-resistance. Similarly, serial profiling of patient CTCs directly from blood confirmed increased ITH post-relapse. These data suggest that treatment-resistance in SCLC is characterized by coexisting subpopulations of cells with heterogeneous gene expression leading to multiple, concurrent resistance mechanisms. These findings emphasize the need for clinical efforts to focus on rational combination therapies for treatment-naïve SCLC tumors to maximize initial responses and counteract the emergence of ITH and diverse resistance mechanisms.

Cross-Species Single-Cell Analysis of Pancreatic Ductal Adenocarcinoma Reveals Antigen-Presenting Cancer-Associated Fibroblasts.

Cancer-associated fibroblasts (CAF) are major players in the progression and drug resistance of pancreatic ductal adenocarcinoma (PDAC). CAFs constitute a diverse cell population consisting of several recently described subtypes, although the extent of CAF heterogeneity has remained undefined. Here we use single-cell RNA sequencing to thoroughly characterize the neoplastic and tumor microenvironment content of human and mouse PDAC tumors. We corroborate the presence of myofibroblastic CAFs and inflammatory CAFs and define their unique gene signatures in vivo. Moreover, we describe a new population of CAFs that express MHC class II and CD74, but do not express classic costimulatory molecules. We term this cell population "antigen-presenting CAFs" and find that they activate CD4+ T cells in an antigen-specific fashion in a model system, confirming their putative immune-modulatory capacity. Our cross-species analysis paves the way for investigating distinct functions of CAF subtypes in PDAC immunity and progression. SIGNIFICANCE: Appreciating the full spectrum of fibroblast heterogeneity in pancreatic ductal adenocarcinoma is crucial to developing therapies that specifically target tumor-promoting CAFs. This work identifies MHC class II-expressing CAFs with a capacity to present antigens to CD4+ T cells, and potentially to modulate the immune response in pancreatic tumors.See related commentary by Belle and DeNardo, p. 1001.This article is highlighted in the In This Issue feature, p. 983.

The Genome of C57BL/6J "Eve", the Mother of the Laboratory Mouse Genome Reference Strain.

Isogenic laboratory mouse strains enhance reproducibility because individual animals are genetically identical. For the most widely used isogenic strain, C57BL/6, there exists a wealth of genetic, phenotypic, and genomic data, including a high-quality reference genome (GRCm38.p6). Now 20 years after the first release of the mouse reference genome, C57BL/6J mice are at least 26 inbreeding generations removed from GRCm38 and the strain is now maintained with periodic reintroduction of cryorecovered mice derived from a single breeder pair, aptly named Adam and Eve. To provide an update to the mouse reference genome that more accurately represents the genome of today's C57BL/6J mice, we took advantage of long read, short read, and optical mapping technologies to generate a de novo assembly of the C57BL/6J Eve genome (B6Eve). Using these data, we have addressed recurring variants observed in previous mouse genomic studies. We have also identified structural variations, closed gaps in the mouse reference assembly, and revealed previously unannotated coding sequences. This B6Eve assembly explains discrepant observations that have been associated with GRCm38-based analyses, and will inform a reference genome that is more representative of the C57BL/6J mice that are in use today.

Single-cell transcriptomic analysis of the lateral hypothalamic area reveals molecularly distinct populations of inhibitory and excitatory neurons.

The lateral hypothalamic area (LHA) coordinates an array of fundamental behaviors, including sleeping, waking, feeding, stress and motivated behavior. The wide spectrum of functions ascribed to the LHA may be explained by a heterogeneous population of neurons, the full diversity of which is poorly understood. We employed a droplet-based single-cell RNA-sequencing approach to develop a comprehensive census of molecularly distinct cell types in the mouse LHA. Neuronal populations were classified based on fast neurotransmitter phenotype and expression of neuropeptides, transcription factors and synaptic proteins, among other gene categories. We define 15 distinct populations of glutamatergic neurons and 15 of GABAergic neurons, including known and novel cell types. We further characterize a novel population of somatostatin-expressing neurons through anatomical and behavioral approaches, identifying a role for these neurons in specific forms of innate locomotor behavior. This study lays the groundwork for better understanding the circuit-level underpinnings of LHA function.

Detection of correlated hidden factors from single cell transcriptomes using Iteratively Adjusted-SVA (IA-SVA).

Single cell RNA-sequencing (scRNA-seq) precisely characterizes gene expression levels and dissects variation in expression associated with the state (technical or biological) and the type of the cell, which is averaged out in bulk measurements. Multiple and correlated sources contribute to gene expression variation in single cells, which makes their estimation difficult with the existing methods developed for batch correction (e.g., surrogate variable analysis (SVA)) that estimate orthogonal transformations of these sources. We developed iteratively adjusted surrogate variable analysis (IA-SVA) that can estimate hidden factors even when they are correlated with other sources of variation by identifying a set of genes associated with each hidden factor in an iterative manner. Analysis of scRNA-seq data from human cells showed that IA-SVA could accurately capture hidden variation arising from technical (e.g., stacked doublet cells) or biological sources (e.g., cell type or cell-cycle stage). Furthermore, IA-SVA delivers a set of genes associated with the detected hidden source to be used in downstream data analyses. As a proof of concept, IA-SVA recapitulated known marker genes for islet cell subsets (e.g., alpha, beta), which improved the grouping of subsets into distinct clusters. Taken together, IA-SVA is an effective and novel method to dissect multiple and correlated sources of variation in scRNA-seq data.

Author Correction: Single cell transcriptome profiling of retinal ganglion cells identifies cellular subtypes.

The original version of the Supplementary Information file associated with this Article contained an error in Supplementary Fig. 2. In panel c, the graph was inadvertently replaced with a duplicate of the graph in panel a. The error has now been fixed and the corrected version Supplementary Information PDF is available to download from the HTML version of the Article.

Single cell transcriptome profiling of retinal ganglion cells identifies cellular subtypes.

Retinal ganglion cells (RGCs) convey the major output of information collected from the eye to the brain. Thirty subtypes of RGCs have been identified to date. Here, we analyze 6225 RGCs (average of 5000 genes per cell) from right and left eyes by single-cell RNA-seq and classify them into 40 subtypes using clustering algorithms. We identify additional subtypes and markers, as well as transcription factors predicted to cooperate in specifying RGC subtypes. Zic1, a marker of the right eye-enriched subtype, is validated by immunostaining in situ. Runx1 and Fst, the markers of other subtypes, are validated in purified RGCs by fluorescent in situ hybridization (FISH) and immunostaining. We show the extent of gene expression variability needed for subtype segregation, and we show a hierarchy in diversification from a cell-type population to subtypes. Finally, we present a website for comparing the gene expression of RGC subtypes.

A Common Type 2 Diabetes Risk Variant Potentiates Activity of an Evolutionarily Conserved Islet Stretch Enhancer and Increases C2CD4A and C2CD4B Expression.

Genome-wide association studies (GWASs) and functional genomics approaches implicate enhancer disruption in islet dysfunction and type 2 diabetes (T2D) risk. We applied genetic fine-mapping and functional (epi)genomic approaches to a T2D- and proinsulin-associated 15q22.2 locus to identify a most likely causal variant, determine its direction of effect, and elucidate plausible target genes. Fine-mapping and conditional analyses of proinsulin levels of 8,635 non-diabetic individuals from the METSIM study support a single association signal represented by a cluster of 16 strongly associated (p < 10-17) variants in high linkage disequilibrium (r2 > 0.8) with the GWAS index SNP rs7172432. These variants reside in an evolutionarily and functionally conserved islet and ß cell stretch or super enhancer; the most strongly associated variant (rs7163757, p = 3 × 10-19) overlaps a conserved islet open chromatin site. DNA sequence containing the rs7163757 risk allele displayed 2-fold higher enhancer activity than the non-risk allele in reporter assays (p < 0.01) and was differentially bound by ß cell nuclear extract proteins. Transcription factor NFAT specifically potentiated risk-allele enhancer activity and altered patterns of nuclear protein binding to the risk allele in vitro, suggesting that it could be a factor mediating risk-allele effects. Finally, the rs7163757 proinsulin-raising and T2D risk allele (C) was associated with increased expression of C2CD4B, and possibly C2CD4A, both of which were induced by inflammatory cytokines, in human islets. Together, these data suggest that rs7163757 contributes to genetic risk of islet dysfunction and T2D by increasing NFAT-mediated islet enhancer activity and modulating C2CD4B, and possibly C2CD4A, expression in (patho)physiologic states.

Single-Cell Transcriptional Profiling Reveals Cellular Diversity and Intercommunication in the Mouse Heart.

Characterization of the cardiac cellulome, the network of cells that form the heart, is essential for understanding cardiac development and normal organ function and for formulating precise therapeutic strategies to combat heart disease. Recent studies have reshaped our understanding of cardiac cellular composition and highlighted important functional roles for non-myocyte cell types. In this study, we characterized single-cell transcriptional profiles of the murine non-myocyte cardiac cellular landscape using single-cell RNA sequencing (scRNA-seq). Detailed molecular analyses revealed the diversity of the cardiac cellulome and facilitated the development of techniques to isolate understudied cardiac cell populations, such as mural cells and glia. Our analyses also revealed extensive networks of intercellular communication and suggested prevalent sexual dimorphism in gene expression in the heart. This study offers insights into the structure and function of the mammalian cardiac cellulome and provides an important resource that will stimulate studies in cardiac cell biology.

Single-cell transcriptomes identify human islet cell signatures and reveal cell-type-specific expression changes in type 2 diabetes.

Blood glucose levels are tightly controlled by the coordinated action of at least four cell types constituting pancreatic islets. Changes in the proportion and/or function of these cells are associated with genetic and molecular pathophysiology of monogenic, type 1, and type 2 (T2D) diabetes. Cellular heterogeneity impedes precise understanding of the molecular components of each islet cell type that govern islet (dys)function, particularly the less abundant delta and gamma/pancreatic polypeptide (PP) cells. Here, we report single-cell transcriptomes for 638 cells from nondiabetic (ND) and T2D human islet samples. Analyses of ND single-cell transcriptomes identified distinct alpha, beta, delta, and PP/gamma cell-type signatures. Genes linked to rare and common forms of islet dysfunction and diabetes were expressed in the delta and PP/gamma cell types. Moreover, this study revealed that delta cells specifically express receptors that receive and coordinate systemic cues from the leptin, ghrelin, and dopamine signaling pathways implicating them as integrators of central and peripheral metabolic signals into the pancreatic islet. Finally, single-cell transcriptome profiling revealed genes differentially regulated between T2D and ND alpha, beta, and delta cells that were undetectable in paired whole islet analyses. This study thus identifies fundamental cell-type-specific features of pancreatic islet (dys)function and provides a critical resource for comprehensive understanding of islet biology and diabetes pathogenesis.

Determining exon connectivity in complex mRNAs by nanopore sequencing.

Short-read high-throughput RNA sequencing, though powerful, is limited in its ability to directly measure exon connectivity in mRNAs that contain multiple alternative exons located farther apart than the maximum read length. Here, we use the Oxford Nanopore MinION sequencer to identify 7,899 'full-length' isoforms expressed from four Drosophila genes, Dscam1, MRP, Mhc, and Rdl. These results demonstrate that nanopore sequencing can be used to deconvolute individual isoforms and that it has the potential to be a powerful method for comprehensive transcriptome characterization.

Regulation of alternative splicing in Drosophila by 56 RNA binding proteins.

Alternative splicing is regulated by RNA binding proteins (RBPs) that recognize pre-mRNA sequence elements and activate or repress adjacent exons. Here, we used RNA interference and RNA-seq to identify splicing events regulated by 56 Drosophila proteins, some previously unknown to regulate splicing. Nearly all proteins affected alternative first exons, suggesting that RBPs play important roles in first exon choice. Half of the splicing events were regulated by multiple proteins, demonstrating extensive combinatorial regulation. We observed that SR and hnRNP proteins tend to act coordinately with each other, not antagonistically. We also identified a cross-regulatory network where splicing regulators affected the splicing of pre-mRNAs encoding other splicing regulators. This large-scale study substantially enhances our understanding of recent models of splicing regulation and provides a resource of thousands of exons that are regulated by 56 diverse RBPs.

Genome-wide identification of zero nucleotide recursive splicing in Drosophila.

Recursive splicing is a process in which large introns are removed in multiple steps by re-splicing at ratchet points--5' splice sites recreated after splicing. Recursive splicing was first identified in the Drosophila Ultrabithorax (Ubx) gene and only three additional Drosophila genes have since been experimentally shown to undergo recursive splicing. Here we identify 197 zero nucleotide exon ratchet points in 130 introns of 115 Drosophila genes from total RNA sequencing data generated from developmental time points, dissected tissues and cultured cells. The sequential nature of recursive splicing was confirmed by identification of lariat introns generated by splicing to and from the ratchet points. We also show that recursive splicing is a constitutive process, that depletion of U2AF inhibits recursive splicing, and that the sequence and function of ratchet points are evolutionarily conserved in Drosophila. Finally, we identify four recursively spliced human genes, one of which is also recursively spliced in Drosophila. Together, these results indicate that recursive splicing is commonly used in Drosophila, occurs in humans, and provides insight into the mechanisms by which some large introns are removed.

Circuitous route to transcription regulation.

In this issue of Molecular Cell, Zhang and colleagues (2013) identify a new class of intron-derived circular RNAs (ciRNAs) and show that they have the potential to enhance transcription of their host gene.

Unexpected diversity and expression of avian endogenous retroviruses.

Endogenous retroviruses (ERVs) were identified and characterized in three avian genomes to gain insight into early retroviral evolution. Using the computer program RetroTector to detect relatively intact ERVs, we identified 500 ERVs in the chicken genome, 150 in the turkey genome, and 1,200 in the zebra finch genome. Previous studies suggested that endogenous alpharetroviruses were present in chicken genomes. In this analysis, a small number of alpharetroviruses were seen in the chicken and turkey genomes; however, these were greatly outnumbered by beta-like, gamma-like, and alphabeta proviruses. While the avian ERVs belonged to the same major groups as mammalian ERVs, they were more heterogeneous. In particular, the beta-like viruses revealed an evolutionary continuum with the gradual acquisition and loss of betaretroviral markers and a transition from beta to alphabeta and then to alpharetroviruses. Thus, it appears that birds may resemble a melting pot for early ERV evolution. Many of the ERVs were integrated in clusters on chromosomes, often near centromeres. About 25% of the chicken ERVs were in or near cellular transcription units; this is nearly random. The majority of these integrations were in the sense orientation in introns. A higher-than-random number of integrations were >100 kb from the nearest gene. Deep-sequencing studies of chicken embryo fibroblasts revealed that about 20% of the 500 ERVs were transcribed and translated. A subset of these were also transcribed in vivo in chickens, showing tissue-specific patterns of expression. IMPORTANCE Studies of avian endogenous retroviruses (ERVs) have given us a glimpse of an earlier retroviral world. Three different classes of ERVs were observed with many features of mammalian retroviruses, as well as some important differences. Many avian ERVs were transcribed and translated.