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

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 \u3e16,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

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\u27s 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

Reticuloendotheliosis Virus Strain T Induces miR-155, Which Targets JARID2 and Promotes Cell Survival ▿

The oncogenic microRNA miR-155 is upregulated by several oncogenic viruses. The precursor of miR-155, termed bic, was first observed to cooperate with myc in chicken B-cell lymphomas induced by avian leukosis proviral integrations. We identified another oncogenic retrovirus, reticuloendotheliosis virus strain T (REV-T), that upregulates miR-155 in chicken embryo fibroblasts. We also observed very high levels of miR-155 in REV-T-induced B-cell lymphomas. To study the role of miR-155 in these tumors, we identified JARID2/Jumonji, a cell cycle regulator and part of a histone methyltransferase complex, as a target of miR-155. The overexpression of miR-155 decreased levels of endogenous JARID2 mRNA. We confirmed that miR-155 directly targets both human and chicken JARID2 by assaying the repression of reporters containing the JARID2 3′-untranslated regions. Further, the overexpression of a sponge complementary to miR-155 in a tumor cell line increased endogenous JARID2 mRNA levels. The overexpression of JARID2 in chicken fibroblasts led to decreased cell numbers and an increase in apoptotic cells. The overexpression of miR-155 rescued cells undergoing cytopathic effect caused by infection with subgroup B avian retroviruses. Therefore, we propose that miR-155 has a prosurvival function that is mediated through the downregulation of targets including JARID2

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. Genome Res 2017 Feb; 27(2):208-222

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

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