Single-cell profiling of the developing mouse brain and spinal cord with split-pool barcoding

Interfacial Nanoinjection-Based Nanoliter Single-Cell Analysis

Single-cell analysis offers unprecedented resolution for the investigation of cellular heterogeneity and the capture of rare cells from large populations. Here, described is a simple method named interfacial nanoinjection (INJ), which can miniaturize various single-cell assays to be performed in nanoliter water-in-oil droplets on standard microwell plates. The INJ droplet handler can adjust droplet volumes for multistep reactions on demand with high precision and excellent monodispersity, and consequently enables a wide range of single-cell assays. Importantly, INJ can be coupled with fluorescence-activated cell sorting (FACS), which is currently the most effective and accurate single-cell sorting and isolation method. FACS-INJ pipelines for high-throughput plate well-based single-cell analyses, including single-cell proliferation, drug-resistance testing, polymerase chain reaction (PCR), reverse-transcription PCR, and whole-genome sequencing are introduced. This FACS-INJ pipeline is compatible with a wide range of samples and can be extended to various single-cell analysis applications in microbiology, cell biology, and biomedical diagnostics.

A deep adversarial variational autoencoder model for dimensionality reduction in single-cell RNA sequencing analysis

BACKGROUND

Single-cell RNA sequencing (scRNA-seq) is an emerging technology that can assess the function of an individual cell and cell-to-cell variability at the single cell level in an unbiased manner. Dimensionality reduction is an essential first step in downstream analysis of the scRNA-seq data. However, the scRNA-seq data are challenging for traditional methods due to their high dimensional measurements as well as an abundance of dropout events (that is, zero expression measurements).

RESULTS

To overcome these difficulties, we propose DR-A (Dimensionality Reduction with Adversarial variational autoencoder), a data-driven approach to fulfill the task of dimensionality reduction. DR-A leverages a novel adversarial variational autoencoder-based framework, a variant of generative adversarial networks. DR-A is well-suited for unsupervised learning tasks for the scRNA-seq data, where labels for cell types are costly and often impossible to acquire. Compared with existing methods, DR-A is able to provide a more accurate low dimensional representation of the scRNA-seq data. We illustrate this by utilizing DR-A for clustering of scRNA-seq data.

CONCLUSIONS

Our results indicate that DR-A significantly enhances clustering performance over state-of-the-art methods.

Single-cell technologies - studying rheumatic diseases one cell at a time

Cells, the basic units of life, have striking differences at transcriptomic, proteomic and epigenomic levels across tissues, organs, organ systems and organisms. The coordination of individual immune cells is essential for the generation of effective immune responses to pathogens while immune tolerance is maintained to protect the host. In rheumatic diseases, when immune responses are dysregulated, pathologically important cells might represent only a small fraction of the immune system. Interrogation of the contributions of individual immune cells to pathogenesis and disease progression should therefore reveal important insights into the complicated aetiology of rheumatic diseases. Technological advances are enabling the high-dimensional dissection of single cells at multiple omics levels, which could facilitate the identification of dysregulated molecular mechanisms in patients with rheumatic diseases and the discovery of new therapeutic targets and biomarkers. The single-cell technologies that have been developed over the past decade and the experimental platforms that enable multi-omics integrative analyses have already made inroads into immunology-related fields of study and have potential for use in rheumatology. Layers of omics data derived from single cells are likely to fundamentally change our understanding of the molecular pathways that underpin the pathogenesis of rheumatic diseases.

CellTagging: combinatorial indexing to simultaneously map lineage and identity at single-cell resolution

Single-cell technologies are offering unparalleled insight into complex biology, revealing the behavior of rare cell populations that are masked in bulk population analyses. One current limitation of single-cell approaches is that lineage relationships are typically lost as a result of cell processing. We recently established a method, CellTagging, permitting the parallel capture of lineage information and cell identity via a combinatorial cell indexing approach. CellTagging integrates with high-throughput single-cell RNA sequencing, where sequential rounds of cell labeling enable the construction of multi-level lineage trees. Here, we provide a detailed protocol to (i) generate complex plasmid and lentivirus CellTag libraries for labeling of cells; (ii) sequentially CellTag cells over the course of a biological process; (iii) profile single-cell transcriptomes via high-throughput droplet-based platforms; and (iv) generate a CellTag expression matrix, followed by clone calling and lineage reconstruction. This lentiviral-labeling approach can be deployed in any organism or in vitro culture system that is amenable to viral transduction to simultaneously profile lineage and identity at single-cell resolution.

Determining sequencing depth in a single-cell RNA-seq experiment

An underlying question for virtually all single-cell RNA sequencing experiments is how to allocate the limited sequencing budget: deep sequencing of a few cells or shallow sequencing of many cells? Here we present a mathematical framework which reveals that, for estimating many important gene properties, the optimal allocation is to sequence at a depth of around one read per cell per gene. Interestingly, the corresponding optimal estimator is not the widely-used plug-in estimator, but one developed via empirical Bayes.

For single-cell RNA-seq experiments the sequencing budget is limited, and how it should be optimally allocated to maximize information is not clear. Here the authors develop a mathematical framework to show that, for estimating many gene properties, the optimal allocation is to sequence at the depth of one read per cell per gene.

Eleven grand challenges in single-cell data science

The recent boom in microfluidics and combinatorial indexing strategies, combined with low sequencing costs, has empowered single-cell sequencing technology. Thousands-or even millions-of cells analyzed in a single experiment amount to a data revolution in single-cell biology and pose unique data science problems. Here, we outline eleven challenges that will be central to bringing this emerging field of single-cell data science forward. For each challenge, we highlight motivating research questions, review prior work, and formulate open problems. This compendium is for established researchers, newcomers, and students alike, highlighting interesting and rewarding problems for the coming years.

Mammalian gene expression variability is explained by underlying cell state

Gene expression variability in mammalian systems plays an important role in physiological and pathophysiological conditions. This variability can come from differential regulation related to cell state (extrinsic) and allele-specific transcriptional bursting (intrinsic). Yet, the relative contribution of these two distinct sources is unknown. Here, we exploit the qualitative difference in the patterns of covariance between these two sources to quantify their relative contributions to expression variance in mammalian cells. Using multiplexed error robust RNA fluorescent in situ hybridization (MERFISH), we measured the multivariate gene expression distribution of 150 genes related to Ca2+ signaling coupled with the dynamic Ca2+ response of live cells to ATP. We show that after controlling for cellular phenotypic states such as size, cell cycle stage, and Ca2+ response to ATP, the remaining variability is effectively at the Poisson limit for most genes. These findings demonstrate that the majority of expression variability results from cell state differences and that the contribution of transcriptional bursting is relatively minimal.

Using single-cell technologies to map the human immune system - implications for nephrology

Advances in single-cell technologies are transforming our understanding of cellular identity. For instance, the application of single-cell RNA sequencing and mass cytometry technologies to the study of immune cell populations in blood, secondary lymphoid organs and the renal tract is helping researchers to map the complex immune landscape within the kidney, define cell ontogeny and understand the relationship of kidney-resident immune cells with their circulating counterparts. These studies also provide insights into the interactions of immune cell populations with neighbouring epithelial and endothelial cells in health, and across a range of kidney diseases and cancer. These data have translational potential and will aid the identification of drug targets and enable better prediction of off-target effects. The application of single-cell technologies to clinical renal biopsy samples, or even cells within urine, will improve diagnostic accuracy and assist with personalized prognostication for patients with various kidney diseases. A comparison of immune cell types in peripheral blood and secondary lymphoid organs in healthy individuals and in patients with systemic autoimmune diseases that affect the kidney will also help to unravel the mechanisms that underpin the breakdown in self-tolerance and propagation of autoimmune responses. Together, these immune cell atlases have the potential to transform nephrology.

Integration of eQTL and a Single-Cell Atlas in the Human Eye Identifies Causal Genes for Age-Related Macular Degeneration

Age-related macular degeneration (AMD) is a leading cause of vision loss. To better understand disease pathogenesis and identify causal genes in GWAS loci for AMD risk, we present a comprehensive database of human retina and retinal pigment epithelium (RPE). Our database comprises macular and non-macular RNA sequencing (RNA-seq) profiles from 129 donors, a genome-wide expression quantitative trait loci (eQTL) dataset that includes macula-specific retina and RPE/choroid, and single-nucleus RNA-seq (NucSeq) from human retina and RPE with subtype resolution from more than 100,000 cells. Using NucSeq, we find enriched expression of AMD candidate genes in RPE cells. We identify 15 putative causal genes for AMD on the basis of co-localization of genetic association signals for AMD risk and eye eQTL, including the genes TSPAN10 and TRPM1. These results demonstrate the value of our human eye database for elucidating genetic pathways and potential therapeutic targets for ocular diseases.

RNA sequencing by direct tagmentation of RNA/DNA hybrids

RNA sequencing is widely used to measure gene expression in biomedical research; therefore, improvements in the simplicity and accuracy of the technology are desirable. All existing RNA sequencing methods rely on the conversion of RNA into double-stranded DNA through reverse transcription followed by second-strand synthesis. The latter step requires additional enzymes and purification, and introduces sequence-dependent bias. Here, we show that Tn5 transposase, which randomly binds and cuts double-stranded DNA, can directly fragment and prime the RNA/DNA heteroduplexes generated by reverse transcription. The primed fragments are then subject to PCR amplification. This provides an approach for simple and accurate RNA characterization and quantification.

Transcriptome profiling by RNA sequencing (RNA-seq) has been widely used to characterize cellular status, but it relies on second-strand complementary DNA (cDNA) synthesis to generate initial material for library preparation. Here we use bacterial transposase Tn5, which has been increasingly used in various high-throughput DNA analyses, to construct RNA-seq libraries without second-strand synthesis. We show that Tn5 transposome can randomly bind RNA/DNA heteroduplexes and add sequencing adapters onto RNA directly after reverse transcription. This method, Sequencing HEteRo RNA-DNA-hYbrid (SHERRY), is versatile and scalable. SHERRY accepts a wide range of starting materials, from bulk RNA to single cells. SHERRY offers a greatly simplified protocol and produces results with higher reproducibility and GC uniformity compared with prevailing RNA-seq methods.

Benchmarking principal component analysis for large-scale single-cell RNA-sequencing

BACKGROUND

Principal component analysis (PCA) is an essential method for analyzing single-cell RNA-seq (scRNA-seq) datasets, but for large-scale scRNA-seq datasets, computation time is long and consumes large amounts of memory.

RESULTS

In this work, we review the existing fast and memory-efficient PCA algorithms and implementations and evaluate their practical application to large-scale scRNA-seq datasets. Our benchmark shows that some PCA algorithms based on Krylov subspace and randomized singular value decomposition are fast, memory-efficient, and more accurate than the other algorithms.

CONCLUSION

We develop a guideline to select an appropriate PCA implementation based on the differences in the computational environment of users and developers.

A benchmark of batch-effect correction methods for single-cell RNA sequencing data

BACKGROUND

Large-scale single-cell transcriptomic datasets generated using different technologies contain batch-specific systematic variations that present a challenge to batch-effect removal and data integration. With continued growth expected in scRNA-seq data, achieving effective batch integration with available computational resources is crucial. Here, we perform an in-depth benchmark study on available batch correction methods to determine the most suitable method for batch-effect removal.

RESULTS

We compare 14 methods in terms of computational runtime, the ability to handle large datasets, and batch-effect correction efficacy while preserving cell type purity. Five scenarios are designed for the study: identical cell types with different technologies, non-identical cell types, multiple batches, big data, and simulated data. Performance is evaluated using four benchmarking metrics including kBET, LISI, ASW, and ARI. We also investigate the use of batch-corrected data to study differential gene expression.

CONCLUSION

Based on our results, Harmony, LIGER, and Seurat 3 are the recommended methods for batch integration. Due to its significantly shorter runtime, Harmony is recommended as the first method to try, with the other methods as viable alternatives.

Tumor Functional Heterogeneity Unraveled by scRNA-seq Technologies

Effective cancer treatment has been precluded by the presence of various forms of intratumoral complexity that drive treatment resistance and metastasis. Recent single-cell sequencing technologies are significantly facilitating the characterization of tumor internal architecture during disease progression. New applications and advances occurring at a fast pace predict an imminent broad application of these technologies in many research areas. As occurred with next-generation sequencing (NGS) technologies, once applied to clinical samples across tumor types, single-cell sequencing technologies could trigger an exponential increase in knowledge of the molecular pathways involved in cancer progression and contribute to the improvement of cancer treatment.

Next-Generation Genotoxicology: Using Modern Sequencing Technologies to Assess Somatic Mutagenesis and Cancer Risk

Mutations have a profound effect on human health, particularly through an increased risk of carcinogenesis and genetic disease. The strong correlation between mutagenesis and carcinogenesis has been a driving force behind genotoxicity research for more than 50 years. The stochastic and infrequent nature of mutagenesis makes it challenging to observe and to study. Indeed, decades have been spent developing increasingly sophisticated assays and methods to study these low-frequency genetic errors, in hopes of better predicting which chemicals may be carcinogens, understanding their mode of action, and informing guidelines to prevent undue human exposure. While effective, widely used genetic selection-based technologies have a number of limitations that have hampered major advancements in the field of genotoxicity. Emerging new tools, in the form of enhanced next-generation sequencing platforms and methods, are changing this paradigm. In this review, we discuss rapidly evolving sequencing tools and technologies, such as error-corrected sequencing and single cell analysis, which we anticipate will fundamentally reshape the field. In addition, we consider a variety emerging applications for these new technologies, including the detection of DNA adducts, inference of mutational processes based on genomic site and local sequence contexts, and evaluation of genome engineering fidelity, as well as other cutting-edge challenges for the next 50 years of environmental and molecular mutagenesis research. Environ. Mol. Mutagen. 61:135-151, 2020. © 2019 The Authors. Environmental and Molecular Mutagenesis published by Wiley Periodicals, Inc. on behalf of Environmental Mutagen Society.

Single-Cell Transcriptional Profiling of the Intestinal Epithelium

Emerging single-cell technologies, like single-cell RNA sequencing (scRNA-seq), enable the study of heterogeneous biological systems at cellular resolution. By profiling the set of expressed transcripts in each cell, single-cell transcriptomics has allowed for the cataloging of the cellular constituents of multiple organs and tissues, both in health and disease. In addition, these technologies have provided mechanistic insights into cellular function, cell state transitions, developmental trajectories and lineage relationships, as well as helped to dissect complex, population-level responses to environmental perturbations. scRNA-seq is particularly useful for characterizing the intestinal epithelium because it is a dynamic, rapidly self-renewing tissue comprised of more than a dozen specialized cell types. Here we discuss the fundamentals of single-cell transcriptomics of the murine small intestinal epithelium. We review the principles of proper experimental design and provide methods for the dissociation of the small intestinal epithelium into single cells followed by fluorescence-activated cell sorting (FACS) and for scRNA-seq using the 10× Genomics Chromium platform.

scRNA-seq in medulloblastoma shows cellular heterogeneity and lineage expansion support resistance to SHH inhibitor therapy

Targeting oncogenic pathways holds promise for brain tumor treatment, but inhibition of Sonic Hedgehog (SHH) signaling has failed in SHH-driven medulloblastoma. Cellular diversity within tumors and reduced lineage commitment can undermine targeted therapy by increasing the probability of treatment-resistant populations. Using single-cell RNA-seq and lineage tracing, we analyzed cellular diversity in medulloblastomas in transgenic, medulloblastoma-prone mice, and responses to the SHH-pathway inhibitor vismodegib. In untreated tumors, we find expected stromal cells and tumor-derived cells showing either a spectrum of neural progenitor-differentiation states or glial and stem cell markers. Vismodegib reduces the proliferative population and increases differentiation. However, specific cell types in vismodegib-treated tumors remain proliferative, showing either persistent SHH-pathway activation or stem cell characteristics. Our data show that even in tumors with a single pathway-activating mutation, diverse mechanisms drive tumor growth. This diversity confers early resistance to targeted inhibitor therapy, demonstrating the need to target multiple pathways simultaneously.

Although the hedgehog (HH) pathway is known to be deregulated in medulloblastoma, inhibitors of the pathway have shown disappointing clinical benefit. Using single-cell sequencing in a mouse model of the disease, the authors show that the response to the HH pathway inhibitor vismodegib is cell-type specific.

SMARTer single cell total RNA sequencing

Single cell RNA sequencing methods have been increasingly used to understand cellular heterogeneity. Nevertheless, most of these methods suffer from one or more limitations, such as focusing only on polyadenylated RNA, sequencing of only the 3' end of the transcript, an exuberant fraction of reads mapping to ribosomal RNA, and the unstranded nature of the sequencing data. Here, we developed a novel single cell strand-specific total RNA library preparation method addressing all the aforementioned shortcomings. Our method was validated on a microfluidics system using three different cancer cell lines undergoing a chemical or genetic perturbation and on two other cancer cell lines sorted in microplates. We demonstrate that our total RNA-seq method detects an equal or higher number of genes compared to classic polyA[+] RNA-seq, including novel and non-polyadenylated genes. The obtained RNA expression patterns also recapitulate the expected biological signal. Inherent to total RNA-seq, our method is also able to detect circular RNAs. Taken together, SMARTer single cell total RNA sequencing is very well suited for any single cell sequencing experiment in which transcript level information is needed beyond polyadenylated genes.

Single-Cell RNA Sequencing of the Cardiovascular System: New Looks for Old Diseases

Cardiovascular disease encompasses a wide range of conditions, resulting in the highest number of deaths worldwide. The underlying pathologies surrounding cardiovascular disease include a vast and complicated network of both cellular and molecular mechanisms. Unique phenotypic alterations in specific cell types, visualized as varying RNA expression-levels (both coding and non-coding), have been identified as crucial factors in the pathology underlying conditions such as heart failure and atherosclerosis. Recent advances in single-cell RNA sequencing (scRNA-seq) have elucidated a new realm of cell subpopulations and transcriptional variations that are associated with normal and pathological physiology in a wide variety of diseases. This breakthrough in the phenotypical understanding of our cells has brought novel insight into cardiovascular basic science. scRNA-seq allows for separation of widely distinct cell subpopulations which were, until recently, simply averaged together with bulk-tissue RNA-seq. scRNA-seq has been used to identify novel cell types in the heart and vasculature that could be implicated in a variety of disease pathologies. Furthermore, scRNA-seq has been able to identify significant heterogeneity of phenotypes within individual cell subtype populations. The ability to characterize single cells based on transcriptional phenotypes allows researchers the ability to map development of cells and identify changes in specific subpopulations due to diseases at a very high throughput. This review looks at recent scRNA-seq studies of various aspects of the cardiovascular system and discusses their potential value to our understanding of the cardiovascular system and pathology.

Molecular profiling of single neurons of known identity in two ganglia from the crab Cancer borealis

Single-cell transcriptional profiling has become a widespread tool in cell identification, particularly in the nervous system, based on the notion that genomic information determines cell identity. However, many cell-type classification studies are unconstrained by other cellular attributes (e.g., morphology, physiology). Here, we systematically test how accurately transcriptional profiling can assign cell identity to well-studied anatomically and functionally identified neurons in 2 small neuronal networks. While these neurons clearly possess distinct patterns of gene expression across cell types, their expression profiles are not sufficient to unambiguously confirm their identity. We suggest that true cell identity can only be determined by combining gene expression data with other cellular attributes such as innervation pattern, morphology, or physiology.

Understanding circuit organization depends on identification of cell types. Recent advances in transcriptional profiling methods have enabled classification of cell types by their gene expression. While exceptionally powerful and high throughput, the ground-truth validation of these methods is difficult: If cell type is unknown, how does one assess whether a given analysis accurately captures neuronal identity? To shed light on the capabilities and limitations of solely using transcriptional profiling for cell-type classification, we performed 2 forms of transcriptional profiling—RNA-seq and quantitative RT-PCR, in single, unambiguously identified neurons from 2 small crustacean neuronal networks: The stomatogastric and cardiac ganglia. We then combined our knowledge of cell type with unbiased clustering analyses and supervised machine learning to determine how accurately functionally defined neuron types can be classified by expression profile alone. The results demonstrate that expression profile is able to capture neuronal identity most accurately when combined with multimodal information that allows for post hoc grouping, so analysis can proceed from a supervised perspective. Solely unsupervised clustering can lead to misidentification and an inability to distinguish between 2 or more cell types. Therefore, this study supports the general utility of cell identification by transcriptional profiling, but adds a caution: It is difficult or impossible to know under what conditions transcriptional profiling alone is capable of assigning cell identity. Only by combining multiple modalities of information such as physiology, morphology, or innervation target can neuronal identity be unambiguously determined.

Single-cell sequencing and its applications in head and neck cancer

Head and neck squamous cell carcinoma (HNSCC), like many tumors, is characterized by significant intra-tumoral heterogeneity, namely transcriptional, genetic, and epigenetic differences that define distinct cellular subpopulations. While it has been established that intra-tumoral heterogeneity may have prognostic significance in HNSCC, we are only beginning to describe and define such heterogeneity at a cellular resolution. Recent advances in single-cell sequencing technologies have been critical in this regard, opening new avenues in our understanding of more nuanced tumor biology by identifying distinct cellular subpopulations, dissecting signaling within the tumor microenvironment, and characterizing cellular genomic mutations and copy number aberrations. The combined effect of these insights is likely to be robust and meaningful changes in existing diagnostic and treatment algorithms through the application of novel biomarkers as well as targeted therapeutics. Here, we review single-cell technological and computational advances at the genomic, transcriptomic, and epigenomic levels, and discuss their applications in cancer research and clinical practice, with a specific focus on HNSCC.

High-throughput sequencing of the transcriptome and chromatin accessibility in the same cell

Single-cell RNA sequencing can reveal the transcriptional state of cells, yet provides little insight into the upstream regulatory landscape associated with open or accessible chromatin regions. Joint profiling of accessible chromatin and RNA within the same cells would permit direct matching of transcriptional regulation to its outputs. Here, we describe droplet-based single-nucleus chromatin accessibility and mRNA expression sequencing (SNARE-seq), a method that can link a cell's transcriptome with its accessible chromatin for sequencing at scale. Specifically, accessible sites are captured by Tn5 transposase in permeabilized nuclei to permit, within many droplets in parallel, DNA barcode tagging together with the mRNA molecules from the same cells. To demonstrate the utility of SNARE-seq, we generated joint profiles of 5,081 and 10,309 cells from neonatal and adult mouse cerebral cortices, respectively. We reconstructed the transcriptome and epigenetic landscapes of major and rare cell types, uncovered lineage-specific accessible sites, especially for low-abundance cells, and connected the dynamics of promoter accessibility with transcription level during neurogenesis.

Single-Cell Analysis Using Droplet Microfluidics

Droplet microfluidics has revolutionized the study of single cells. The ability to compartmentalize cells within picoliter droplets in microfluidic devices has opened up a wide range of strategies to extract information at the genomic, transcriptomic, proteomic, or metabolomic level from large numbers of individual cells. Studying the different molecular landscapes at single-cell resolution has provided the authors with a detailed picture of intracellular heterogeneity and the resulting changes in cellular phenotypes. In addition, these technologies have aided in the discovery of rare cells in tumors or in the immune system, and left the authors with a deeper understanding of the fundamental biological processes that determine cell fate. This review aims to provide a detailed overview of the various droplet microfluidic strategies reported in the literature, taking into account the sometimes subtle differences in workflow or reagents that enable or improve certain protocols. Specifically, approaches to targeted- and whole-genome analysis, as well as whole-transcriptome profiling techniques, are reviewed. In addition, an up-to-date overview of new methods to characterize and quantify single-cell protein levels, and of developments to screen secreted molecules such as antibodies, cytokines, or metabolites at the single-cell level, is provided.

Using single-cell transcriptomics to understand functional states and interactions in microbial eukaryotes

Understanding the diversity and evolution of eukaryotic microorganisms remains one of the major challenges of modern biology. In recent years, we have advanced in the discovery and phylogenetic placement of new eukaryotic species and lineages, which in turn completely transformed our view on the eukaryotic tree of life. But we remain ignorant of the life cycles, physiology and cellular states of most of these microbial eukaryotes, as well as of their interactions with other organisms. Here, we discuss how high-throughput genome-wide gene expression analysis of eukaryotic single cells can shed light on protist biology. First, we review different single-cell transcriptomics methodologies with particular focus on microbial eukaryote applications. Then, we discuss single-cell gene expression analysis of protists in culture and what can be learnt from these approaches. Finally, we envision the application of single-cell transcriptomics to protist communities to interrogate not only community components, but also the gene expression signatures of distinct cellular and physiological states, as well as the transcriptional dynamics of interspecific interactions. Overall, we argue that single-cell transcriptomics can significantly contribute to our understanding of the biology of microbial eukaryotes. This article is part of a discussion meeting issue 'Single cell ecology'.

Establishing neuronal diversity in the spinal cord: a time and a place

The vertebrate spinal cord comprises multiple functionally distinct neuronal cell types arranged in characteristic positions. During development, these different types of neurons differentiate from transcriptionally distinct neural progenitors that are arrayed in discrete domains along the dorsal-ventral and anterior-posterior axes of the embryonic spinal cord. This organization arises in response to morphogen gradients acting upstream of a gene regulatory network, the architecture of which determines the spatial and temporal pattern of gene expression. In recent years, substantial progress has been made in deciphering the regulatory network that underlies the specification of distinct progenitor and neuronal cell identities. In this Review, we outline how distinct neuronal cell identities are established in response to spatial and temporal patterning systems, and outline novel experimental approaches to study the emergence and function of neuronal diversity in the spinal cord.

Advances and applications of environmental stress adaptation research

Evolution has produced animals that survive extreme fluctuations in environmental conditions including freezing temperatures, anoxia, desiccating conditions, and prolonged periods without food. For example, the wood frog survives whole-body freezing every winter, arresting all gross physiological functions, but recovers functions upon thawing in the spring. Likewise, many small mammals hibernate for months at a time with minimal metabolic activity, organ perfusion, and movement, yet do not suffer significant muscle atrophy upon arousal. These conditions and the biochemical adaptations employed to deal with them can be viewed as Nature's answer to problems that humans wish to answer, particularly in a biomedical context. This review focuses on recent advances in the field of animal environmental stress adaptation, starting with an emphasis on new areas of research such as epigenetics and microRNA. We then examine new and emerging technologies such as genome editing, novel sequencing applications, and single cell analysis and how these can push us closer to a deeper understanding of biochemical adaptation. Next, evaluate the potential contributions of new high-throughput technologies (e.g. next-generation sequencing, mass spectrometry proteomics) to better understanding the adaptations that support these extreme phenotypes. Concluding, we examine some of the human applications that can be gained from understanding the principles of biochemical adaptation including organ preservation and treatments for conditions such as ischemic stroke and muscle disuse atrophy.

Copolymerization of single-cell nucleic acids into balls of acrylamide gel

We show the use of 5′-Acrydite oligonucleotides to copolymerize single-cell DNA or RNA into balls of acrylamide gel (BAGs). Combining this step with split-and-pool techniques for creating barcodes yields a method with advantages in cost and scalability, depth of coverage, ease of operation, minimal cross-contamination, and efficient use of samples. We perform DNA copy number profiling on mixtures of cell lines, nuclei from frozen prostate tumors, and biopsy washes. As applied to RNA, the method has high capture efficiency of transcripts and sufficient consistency to clearly distinguish the expression patterns of cell lines and individual nuclei from neurons dissected from the mouse brain. By using varietal tags (UMIs) to achieve sequence error correction, we show extremely low levels of cross-contamination by tracking source-specific SNVs. The method is readily modifiable, and we will discuss its adaptability and diverse applications.

The in vivo endothelial cell translatome is highly heterogeneous across vascular beds

Endothelial cells (ECs) are highly specialized across vascular beds. However, given their interspersed anatomic distribution, comprehensive characterization of the molecular basis for this heterogeneity in vivo has been limited. By applying endothelial-specific translating ribosome affinity purification (EC-TRAP) combined with high-throughput RNA sequencing analysis, we identified pan EC-enriched genes and tissue-specific EC transcripts, which include both established markers and genes previously unappreciated for their presence in ECs. In addition, EC-TRAP limits changes in gene expression after EC isolation and in vitro expansion, as well as rapid vascular bed-specific shifts in EC gene expression profiles as a result of the enzymatic tissue dissociation required to generate single-cell suspensions for fluorescence-activated cell sorting or single-cell RNA sequencing analysis. Comparison of our EC-TRAP with published single-cell RNA sequencing data further demonstrates considerably greater sensitivity of EC-TRAP for the detection of low abundant transcripts. Application of EC-TRAP to examine the in vivo host response to lipopolysaccharide (LPS) revealed the induction of gene expression programs associated with a native defense response, with marked differences across vascular beds. Furthermore, comparative analysis of whole-tissue and TRAP-selected mRNAs identified LPS-induced differences that would not have been detected by whole-tissue analysis alone. Together, these data provide a resource for the analysis of EC-specific gene expression programs across heterogeneous vascular beds under both physiologic and pathologic conditions.

Pou2f2 Regulates the Distribution of Dorsal Interneurons in the Mouse Developing Spinal Cord

Spinal dorsal interneurons, which are generated during embryonic development, relay and process sensory inputs from the periphery to the central nervous system. Proper integration of these cells into neuronal circuitry depends on their correct positioning within the spinal parenchyma. Molecular cues that control neuronal migration have been extensively characterized but the genetic programs that regulate their production remain poorly investigated. Onecut (OC) transcription factors have been shown to control the migration of the dorsal interneurons (dINs) during spinal cord development. Here, we report that the OC factors moderate the expression of Pou2f2, a transcription factor essential for B-cell differentiation, in spinal dINs. Overexpression or inactivation of Pou2f2 leads to alterations in the differentiation of dI2, dI3 and Phox2a-positive dI5 populations and to defects in the distribution of dI2-dI6 interneurons. Thus, an OC-Pou2f2 genetic cascade regulates adequate diversification and distribution of dINs during embryonic development.

An ultra high-throughput method for single-cell joint analysis of open chromatin and transcriptome

Simultaneous profiling of transcriptome and chromatin accessibility within single cells is a powerful approach to dissect gene regulatory programs in complex tissues. However, current tools are limited by modest throughput. We now describe an ultra high-throughput method, Paired-seq, for parallel analysis of transcriptome and accessible chromatin in millions of single cells. We demonstrate the utility of Paired-seq for analyzing the dynamic and cell-type-specific gene regulatory programs in complex tissues by applying it to mouse adult cerebral cortex and fetal forebrain. The joint profiles of a large number of single cells allowed us to deconvolute the transcriptome and open chromatin landscapes in the major cell types within these brain tissues, infer putative target genes of candidate enhancers, and reconstruct the trajectory of cellular lineages within the developing forebrain.

Copper and the brain noradrenergic system

Copper (Cu) plays an essential role in the development and function of the brain. In humans, genetic disorders of Cu metabolism may cause either severe Cu deficiency (Menkes disease) or excessive Cu accumulation (Wilson disease) in the brain tissue. In either case, the loss of Cu homeostasis results in catecholamine misbalance, abnormal myelination of neurons, loss of normal brain architecture, and a spectrum of neurologic and/or psychiatric manifestations. Several metabolic processes have been identified as particularly sensitive to Cu dis-homeostasis. This review focuses on the role of Cu in noradrenergic neurons and summarizes the current knowledge of mechanisms that maintain Cu homeostasis in these cells. The impact of Cu misbalance on catecholamine metabolism and functioning of noradrenergic system is discussed.

Single-Cell RNA Sequencing of Plant-Associated Bacterial Communities

Plants in soil are not solitary, hence continually interact with and obtain benefits from a community of microbes ("microbiome"). The meta-functional output from the microbiome results from complex interactions among the different community members with distinct taxonomic identities and metabolic capacities. Particularly, the bacterial communities of the root surface are spatially organized structures composed of root-attached biofilms and planktonic cells arranged in complex layers. With the distinct but coordinated roles among the different member cells, bacterial communities resemble properties of a multicellular organism. High throughput sequencing technologies have allowed rapid and large-scale analysis of taxonomic composition and metabolic capacities of bacterial communities. However, these methods are generally unable to reconstruct the assembly of these communities, or how the gene expression patterns in individual cells/species are coordinated within these communities. Single-cell transcriptomes of community members can identify how gene expression patterns vary among members of the community, including differences among different cells of the same species. This information can be used to classify cells based on functional gene expression patterns, and predict the spatial organization of the community. Here we discuss strategies for the isolation of single bacterial cells, mRNA enrichment, library construction, and analysis and interpretation of the resulting single-cell RNA-Seq datasets. Unraveling regulatory and metabolic processes at the single cell level is expected to yield an unprecedented discovery of mechanisms involved in bacterial recruitment, attachment, assembly, organization of the community, or in the specific interactions among the different members of these communities.

Selective demethylation of two CpG sites causes postnatal activation of the Dao gene and consequent removal of D-serine within the mouse cerebellum

BACKGROUND

Programmed epigenetic modifications occurring at early postnatal brain developmental stages may have a long-lasting impact on brain function and complex behavior throughout life. Notably, it is now emerging that several genes that undergo perinatal changes in DNA methylation are associated with neuropsychiatric disorders. In this context, we envisaged that epigenetic modifications during the perinatal period may potentially drive essential changes in the genes regulating brain levels of critical neuromodulators such as D-serine and D-aspartate. Dysfunction of this fine regulation may contribute to the genesis of schizophrenia or other mental disorders, in which altered levels of D-amino acids are found. We recently demonstrated that Ddo, the D-aspartate degradation gene, is actively demethylated to ultimately reduce D-aspartate levels. However, the role of epigenetics as a mechanism driving the regulation of appropriate D-ser levels during brain development has been poorly investigated to date.

METHODS

We performed comprehensive ultradeep DNA methylation and hydroxymethylation profiling along with mRNA expression and HPLC-based D-amino acids level analyses of genes controlling the mammalian brain levels of D-serine and D-aspartate. DNA methylation changes occurring in specific cerebellar cell types were also investigated. We conducted high coverage targeted bisulfite sequencing by next-generation sequencing and single-molecule bioinformatic analysis.

RESULTS

We report consistent spatiotemporal modifications occurring at the Dao gene during neonatal development in a specific brain region (the cerebellum) and within specific cell types (astrocytes) for the first time. Dynamic demethylation at two specific CpG sites located just downstream of the transcription start site was sufficient to strongly activate the Dao gene, ultimately promoting the complete physiological degradation of cerebellar D-serine a few days after mouse birth. High amount of 5'-hydroxymethylcytosine, exclusively detected at relevant CpG sites, strongly evoked the occurrence of an active demethylation process.

CONCLUSION

The present investigation demonstrates that robust and selective demethylation of two CpG sites is associated with postnatal activation of the Dao gene and consequent removal of D-serine within the mouse cerebellum. A single-molecule methylation approach applied at the Dao locus promises to identify different cell-type compositions and functions in different brain areas and developmental stages.

Single-cell transcriptomes and whole-brain projections of serotonin neurons in the mouse dorsal and median raphe nuclei

Serotonin neurons of the dorsal and median raphe nuclei (DR, MR) collectively innervate the entire forebrain and midbrain, modulating diverse physiology and behavior. To gain a fundamental understanding of their molecular heterogeneity, we used plate-based single-cell RNA-sequencing to generate a comprehensive dataset comprising eleven transcriptomically distinct serotonin neuron clusters. Systematic in situ hybridization mapped specific clusters to the principal DR, caudal DR, or MR. These transcriptomic clusters differentially express a rich repertoire of neuropeptides, receptors, ion channels, and transcription factors. We generated novel intersectional viral-genetic tools to access specific subpopulations. Whole-brain axonal projection mapping revealed that DR serotonin neurons co-expressing vesicular glutamate transporter-3 preferentially innervate the cortex, whereas those co-expressing thyrotropin-releasing hormone innervate subcortical regions in particular the hypothalamus. Reconstruction of 50 individual DR serotonin neurons revealed diverse and segregated axonal projection patterns at the single-cell level. Together, these results provide a molecular foundation of the heterogenous serotonin neuronal phenotypes.

Application of Computational Biology to Decode Brain Transcriptomes

The rapid development of high-throughput sequencing technologies has generated massive valuable brain transcriptome atlases, providing great opportunities for systematically investigating gene expression characteristics across various brain regions throughout a series of developmental stages. Recent studies have revealed that the transcriptional architecture is the key to interpreting the molecular mechanisms of brain complexity. However, our knowledge of brain transcriptional characteristics remains very limited. With the immense efforts to generate high-quality brain transcriptome atlases, new computational approaches to analyze these high-dimensional multivariate data are greatly needed. In this review, we summarize some public resources for brain transcriptome atlases and discuss the general computational pipelines that are commonly used in this field, which would aid in making new discoveries in brain development and disorders.

Systematic Identification of Cell-Cell Communication Networks in the Developing Brain

Since the generation of cell-type specific knockout models, the importance of inter-cellular communication between neural, vascular, and microglial cells during neural development has been increasingly appreciated. However, the extent of communication between these major cell populations remains to be systematically mapped. Here, we describe EMBRACE (embryonic brain cell extraction using FACS), a method to simultaneously isolate neural, mural, endothelial, and microglial cells to more than 94% purity in ∼4 h. Utilizing EMBRACE we isolate, transcriptionally analyze, and build a cell-cell communication map of the developing mouse brain. We identify 1,710 unique ligand-receptor interactions between neural, endothelial, mural, and microglial cells in silico and experimentally confirm the APOE-LDLR, APOE-LRP1, VTN-KDR, and LAMA4-ITGB1 interactions in the E14.5 brain. We provide our data via the searchable “Brain interactome explorer”, available at https://mpi-ie.shinyapps.io/braininteractomeexplorer/. Together, this study provides a comprehensive map that reveals the richness of communication within the developing brain.

•

Isolation of embryonic neural, mural, endothelial, and microglial cells to >94% purity

•

Transcriptome analyses of neural, vascular, and microglial cells from E14.5 brain

•

Generation of inter-cellular communication network with 1,710 unique interactions

•

Established “Brain interactome explorer,” a searchable cell communication database

Association of Myoinositol Transporters with Schizophrenia and Bipolar Disorder: Evidence from Human and Animal Studies

Evidence from animal and human studies has linked myo-inositol (MI) with the pathophysiology and/or treatment of psychiatric disorders such as schizophrenia and bipolar disorder. However, there is still controversy surrounding the definitive role of MI in these disorders. Given that brain MI is differentially regulated by three transporters - SMIT1, SMIT2 and/or HMIT (encoded by the genes: SLC5A3, SLC5A11, and SLC2A13, respectively) - we used available datasets to describe the distribution in mouse and human brain of the different MI transporters and to examine changes in mRNA expression of these transporters in patients with schizophrenia and bipolar disorder. We found a differential distribution of the mRNA of each of the three MI transporters in both human and mouse brain regions. Interestingly, while individual neurons express SMIT1 and HMIT, non-neuronal cells express SMIT2, thus partially accounting for different uptake levels of MI and concordance to downstream second messenger signaling pathways. We also found that the expression of MI transporters is significantly changed in schizophrenia and bipolar disorder in a diagnostic-, brain region- and subtype-specific manner. We then examined the effects of germline deletion in mice of Slc5a3 on behavioral phenotypes related to schizophrenia and bipolar disorder. This gene deletion produces behavioral deficits that mirror some specific symptoms of schizophrenia and bipolar disorder. Finally, chronic administration of MI was able to reverse particular, but not all, behavioral deficits in Slc5a3 knockout mice; MI itself induced some behavioral deficits. Our data support a strong correlation between the expression of MI transporters and schizophrenia and bipolar disorder, and suggest that brain region-specific aberration of one or more of these transporters determines the partial behavioral phenotypes and/or symptomatic pattern of these disorders.

Heterogeneity and plasticity in healthy and atherosclerotic vasculature explored by single-cell sequencing

Cellular characteristics and their adjustment to a state of disease have become more evident due to recent advances in imaging, fluorescent reporter mice, and whole genome RNA sequencing. The uncovered cellular heterogeneity and/or plasticity potentially complicates experimental studies and clinical applications, as markers derived from whole tissue 'bulk' sequencing is unable to yield a subtype transcriptome and specific markers. Here, we propose definitions on heterogeneity and plasticity, discuss current knowledge thereof in the vasculature and how this may be improved by single-cell sequencing (SCS). SCS is emerging as an emerging technique, enabling researchers to investigate different cell populations in more depth than ever before. Cell selection methods, e.g. flow assisted cell sorting, and the quantity of cells can influence the choice of SCS method. Smart-Seq2 offers sequencing of the complete mRNA molecule on a low quantity of cells, while Drop-seq is possible on large numbers of cells on a more superficial level. SCS has given more insight in heterogeneity in healthy vasculature, where it revealed that zonation is crucial in gene expression profiles among the anatomical axis. In diseased vasculature, this heterogeneity seems even more prominent with discovery of new immune subsets in atherosclerosis as proof. Vascular smooth muscle cells and mesenchymal cells also share these plastic characteristics with the ability to up-regulate markers linked to stem cells, such as Sca-1 or CD34. Current SCS studies show some limitations to the number of replicates, quantity of cells used, or the loss of spatial information. Bioinformatical tools could give some more insight in current datasets, making use of pseudo-time analysis or RNA velocity to investigate cell differentiation or polarization. In this review, we discuss the use of SCS in unravelling heterogeneity in the vasculature, its current limitations and promising future applications.

Single-cell transcriptomic profiling of the aging mouse brain

The mammalian brain is complex, with multiple cell types performing a variety of diverse functions, but exactly how each cell type is affected in aging remains largely unknown. Here we performed a single-cell transcriptomic analysis of young and old mouse brains. We provide comprehensive datasets of aging-related genes, pathways and ligand-receptor interactions in nearly all brain cell types. Our analysis identified gene signatures that vary in a coordinated manner across cell types and gene sets that are regulated in a cell-type specific manner, even at times in opposite directions. These data reveal that aging, rather than inducing a universal program, drives a distinct transcriptional course in each cell population, and they highlight key molecular processes, including ribosome biogenesis, underlying brain aging. Overall, these large-scale datasets (accessible online at https://portals.broadinstitute.org/single_cell/study/aging-mouse-brain ) provide a resource for the neuroscience community that will facilitate additional discoveries directed towards understanding and modifying the aging process.

A scalable platform for the development of cell-type-specific viral drivers

Enhancers are the primary DNA regulatory elements that confer cell type specificity of gene expression. Recent studies characterizing individual enhancers have revealed their potential to direct heterologous gene expression in a highly cell-type-specific manner. However, it has not yet been possible to systematically identify and test the function of enhancers for each of the many cell types in an organism. We have developed PESCA, a scalable and generalizable method that leverages ATAC- and single-cell RNA-sequencing protocols, to characterize cell-type-specific enhancers that should enable genetic access and perturbation of gene function across mammalian cell types. Focusing on the highly heterogeneous mammalian cerebral cortex, we apply PESCA to find enhancers and generate viral reagents capable of accessing and manipulating a subset of somatostatin-expressing cortical interneurons with high specificity. This study demonstrates the utility of this platform for developing new cell-type-specific viral reagents, with significant implications for both basic and translational research.

scds: computational annotation of doublets in single-cell RNA sequencing data

MOTIVATION

Single-cell RNA sequencing (scRNA-seq) technologies enable the study of transcriptional heterogeneity at the resolution of individual cells and have an increasing impact on biomedical research. However, it is known that these methods sometimes wrongly consider two or more cells as single cells, and that a number of so-called doublets is present in the output of such experiments. Treating doublets as single cells in downstream analyses can severely bias a study's conclusions, and therefore computational strategies for the identification of doublets are needed.

RESULTS

With scds, we propose two new approaches for in silico doublet identification: Co-expression based doublet scoring (cxds) and binary classification based doublet scoring (bcds). The co-expression based approach, cxds, utilizes binarized (absence/presence) gene expression data and, employing a binomial model for the co-expression of pairs of genes, yields interpretable doublet annotations. bcds, on the other hand, uses a binary classification approach to discriminate artificial doublets from original data. We apply our methods and existing computational doublet identification approaches to four datasets with experimental doublet annotations and find that our methods perform at least as well as the state of the art, at comparably little computational cost. We observe appreciable differences between methods and across datasets and that no approach dominates all others. In summary, scds presents a scalable, competitive approach that allows for doublet annotation of datasets with thousands of cells in a matter of seconds.

AVAILABILITY AND IMPLEMENTATION

scds is implemented as a Bioconductor R package (doi: 10.18129/B9.bioc.scds).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Supervised classification enables rapid annotation of cell atlases

Single-cell molecular profiling technologies are gaining rapid traction, but the manual process by which resulting cell types are typically annotated is labor intensive and rate-limiting. We describe Garnett, a tool for rapidly annotating cell types in single-cell transcriptional profiling and single-cell chromatin accessibility datasets, based on an interpretable, hierarchical markup language of cell type-specific genes. Garnett successfully classifies cell types in tissue and whole organism datasets, as well as across species.

Redefining the Etiologic Landscape of Cerebellar Malformations

Cerebellar malformations are diverse congenital anomalies frequently associated with developmental disability. Although genetic and prenatal non-genetic causes have been described, no systematic analysis has been performed. Here, we present a large-exome sequencing study of Dandy-Walker malformation (DWM) and cerebellar hypoplasia (CBLH). We performed exome sequencing in 282 individuals from 100 families with DWM or CBLH, and we established a molecular diagnosis in 36 of 100 families, with a significantly higher yield for CBLH (51%) than for DWM (16%). The 41 variants impact 27 neurodevelopmental-disorder-associated genes, thus demonstrating that CBLH and DWM are often features of monogenic neurodevelopmental disorders. Though only seven monogenic causes (19%) were identified in more than one individual, neuroimaging review of 131 additional individuals confirmed cerebellar abnormalities in 23 of 27 genetic disorders (85%). Prenatal risk factors were frequently found among individuals without a genetic diagnosis (30 of 64 individuals [47%]). Single-cell RNA sequencing of prenatal human cerebellar tissue revealed gene enrichment in neuronal and vascular cell types; this suggests that defective vasculogenesis may disrupt cerebellar development. Further, de novo gain-of-function variants in PDGFRB, a tyrosine kinase receptor essential for vascular progenitor signaling, were associated with CBLH, and this discovery links genetic and non-genetic etiologies. Our results suggest that genetic defects impact specific cerebellar cell types and implicate abnormal vascular development as a mechanism for cerebellar malformations. We also confirmed a major contribution for non-genetic prenatal factors in individuals with cerebellar abnormalities, substantially influencing diagnostic evaluation and counseling regarding recurrence risk and prognosis.

Emerging Roles of Synapse Organizers in the Regulation of Critical Periods

Experience remodels cortical connectivity during developmental windows called critical periods. Experience-dependent regulation of synaptic strength during these periods establishes circuit functions that are stabilized as critical period plasticity wanes. These processes have been extensively studied in the developing visual cortex, where critical period opening and closure are orchestrated by the assembly, maturation, and strengthening of distinct synapse types. The synaptic specificity of these processes points towards the involvement of distinct molecular pathways. Attractive candidates are pre- and postsynaptic transmembrane proteins that form adhesive complexes across the synaptic cleft. These synapse-organizing proteins control synapse development and maintenance and modulate structural and functional properties of synapses. Recent evidence suggests that they have pivotal roles in the onset and closure of the critical period for vision. In this review, we describe roles of synapse-organizing adhesion molecules in the regulation of visual critical period plasticity and we discuss the potential they offer to restore circuit functions in amblyopia and other neurodevelopmental disorders.

Single-cell technologies in reproductive immunology

The maternal-fetal interface represents a unique immune privileged site that maintains the ability to defend against pathogens while orchestrating the necessary tissue remodeling required for placentation. The recent discovery of novel cellular families (innate lymphoid cells, tissue-resident NK cells) suggests that our understanding of the decidual immunome is incomplete. To understand this complex milieu, new technological developments allow reproductive immunologists to collect increasingly complex data at a cellular resolution. Polychromatic flow cytometry allows for greater resolution in the identification of novel cell types by surface and intracellular protein. Single-cell RNA-seq coupled with microfluidics allows for efficient cellular transcriptomics. The extreme dimensionality and size of data sets generated, however, requires the application of novel computational approaches for unbiased analysis. There are now multiple dimensionality reduction (tSNE, SPADE) and visualization tools (SPICE) that allow researchers to efficiently analyze flow cytometry data. Development of computational tools has also been extended to RNA-seq data (including scRNA-seq), which requires specific analytical tools. Here, we provide an overview and a brief primer for the reproductive immunology community on data acquisition and computational tools for the analysis of complex flow cytometry and RNA-seq data.