Keratinocytes (KCs) are the most frequent cells in the epidermis, and they are often isolated and cultured in vitro to study the molecular biology of the skin. Cultured primary cells and various immortalized cells have been frequently used as skin models but their comparability to intact skin has been questioned. Moreover, when analyzing KC transcriptomes, fluctuation of polyA+ RNA content during the KCs' lifecycle has been omitted.

The mammalian cerebral cortex supports cognitive functions such as sensorimotor integration, memory, and social behaviors. Normal brain function relies on a diverse set of differentiated cell types, including neurons, glia, and vasculature. Here, we have used large-scale single-cell RNA sequencing (RNA-seq) to classify cells in the mouse somatosensory cortex and hippocampal CA1 region. We found 47 molecularly distinct subclasses, comprising all known major cell types in the cortex. We identified numerous marker genes, which allowed alignment with known cell types, morphology, and location. We found a layer I interneuron expressing Pax6 and a distinct postmitotic oligodendrocyte subclass marked by Itpr2. Across the diversity of cortical cell types, transcription factors formed a complex, layered regulatory code, suggesting a mechanism for the maintenance of adult cell type identity.

Cell-free DNA has been used for fetal rhesus factor and sex determination, fetal aneuploidy screening, cancer diagnostics and monitoring, and other applications. However current methods of using cell free DNA require amplification, which leads to allelic dropout and bias especially when starting with small amounts of DNA. Here we describe an amplification-free method for sequencing of cell-free DNA, even from low levels of starting material. We evaluated this method in the context of prenatal diagnosis of fetal aneuploidy and compared it with a PCR-based library preparation method as well as a recently described method using unique molecular identifiers (UMI). All methods performed well, however coverage was increased by the amplification-free method and GC-induced bias was reduced by both the amplification-free method and the UMI method. Future diagnostic applications including whole genome sequencing of cell-free DNA will benefit from amplification-free sequencing.

Evidence for distinct human cancer stem cells (CSCs) remains contentious and the degree to which different cancer cells contribute to propagating malignancies in patients remains unexplored. In low- to intermediate-risk myelodysplastic syndromes (MDS), we establish the existence of rare multipotent MDS stem cells (MDS-SCs), and their hierarchical relationship to lineage-restricted MDS progenitors. All identified somatically acquired genetic lesions were backtracked to distinct MDS-SCs, establishing their distinct MDS-propagating function in vivo. In isolated del(5q)-MDS, acquisition of del(5q) preceded diverse recurrent driver mutations. Sequential analysis in del(5q)-MDS revealed genetic evolution in MDS-SCs and MDS-progenitors prior to leukemic transformation. These findings provide definitive evidence for rare human MDS-SCs in vivo, with extensive implications for the targeting of the cells required and sufficient for MDS-propagation.

PAIRED (PRD)-like homeobox genes belong to a class of predicted transcription factor genes. Several of these PRD-like homeobox genes have been predicted in silico from genomic sequence but until recently had no evidence of transcript expression. We found recently that nine PRD-like homeobox genes, ARGFX, CPHX1, CPHX2, DPRX, DUXA, DUXB, NOBOX, TPRX1 and TPRX2, were expressed in human preimplantation embryos. In the current study we characterized these PRD-like homeobox genes in depth and studied their functions as transcription factors. We cloned multiple transcript variants from human embryos and showed that the expression of these genes is specific to embryos and pluripotent stem cells. Overexpression of the genes in human embryonic stem cells confirmed their roles as transcription factors as either activators (CPHX1, CPHX2, ARGFX) or repressors (DPRX, DUXA, TPRX2) with distinct targets that could be explained by the amino acid sequence in homeodomain. Some PRD-like homeodomain transcription factors had high concordance of target genes and showed enrichment for both developmentally important gene sets and a 36 bp DNA recognition motif implicated in Embryo Genome Activation (EGA). Our data implicate a role for these previously uncharacterized PRD-like homeodomain proteins in the regulation of human embryo genome activation and preimplantation embryo development.

We have developed a high-throughput amplification method for generating robust gene expression profiles using single cell or low RNA inputs.

Tumor heterogeneity in breast cancer tumors is today widely recognized. Most of the available knowledge in genetic variation however, relates to the primary tumor while metastatic lesions are much less studied. Many studies have revealed marked alterations of standard prognostic and predictive factors during tumor progression. Characterization of paired primary- and metastatic tissues should therefore be fundamental in order to understand mechanisms of tumor progression, clonal relationship to tumor evolution as well as the therapeutic aspects of systemic disease.

Developmental cell death plays an important role in the construction of functional neural circuits. In vertebrates, the canonical view proposes a selection of the surviving neurons through stochastic competition for target-derived neurotrophic signals, implying an equal potential for neurons to compete. Here we show an alternative cell fitness selection of neurons that is defined by a specific neuronal heterogeneity code. Proprioceptive sensory neurons that will undergo cell death and those that will survive exhibit different molecular signatures that are regulated by retinoic acid and transcription factors, and are independent of the target and neurotrophins. These molecular features are genetically encoded, representing two distinct subgroups of neurons with contrasted functional maturation states and survival outcome. Thus, in this model, a heterogeneous code of intrinsic cell fitness in neighboring neurons provides differential competitive advantage resulting in the selection of cells with higher capacity to survive and functionally integrate into neural networks.

Single-cell RNA-seq has become routine for discovering cell types and revealing cellular diversity, but archived human brain samples still pose a challenge to current high-throughput platforms. We present STRT-seq-2i, an addressable 9600-microwell array platform, combining sampling by limiting dilution or FACS, with imaging and high throughput at competitive cost. We applied the platform to fresh single mouse cortical cells and to frozen post-mortem human cortical nuclei, matching the performance of a previous lower-throughput platform while retaining a high degree of flexibility, potentially also for other high-throughput applications.

The adult human brain comprises more than a thousand distinct neuronal and glial cell types, a diversity that emerges during early brain development. To reveal the precise sequence of events during early brain development, we used single-cell RNA sequencing and spatial transcriptomics and uncovered cell states and trajectories in human brains at 5 to 14 postconceptional weeks (pcw). We identified 12 major classes that are organized as ~600 distinct cell states, which map to precise spatial anatomical domains at 5 pcw. We described detailed differentiation trajectories of the human forebrain and midbrain and found a large number of region-specific glioblasts that mature into distinct pre-astrocytes and pre-oligodendrocyte precursor cells. Our findings reveal the establishment of cell types during the first trimester of human brain development.

Traditionally, neuroscientists have defined the identity of neurons by the cells' location, morphology, connectivity and excitability. However, the direct relationship between these parameters and the molecular phenotypes has remained largely unexplored. Here, we present a method for obtaining full transcriptome data from single neocortical pyramidal cells and interneurons after whole-cell patch-clamp recordings in mouse brain slices. In our approach, termed Patch-seq, a patch-clamp stimulus protocol is followed by the aspiration of the entire somatic compartment into the recording pipette, reverse transcription of RNA including addition of unique molecular identifiers, cDNA amplification, Illumina library preparation and sequencing. We show that Patch-seq reveals a close link between electrophysiological characteristics, responses to acute chemical challenges and RNA expression of neurotransmitter receptors and channels. Moreover, it distinguishes neuronal subpopulations that correspond to both well-established and, to our knowledge, hitherto undescribed neuronal subtypes. Our findings demonstrate the ability of Patch-seq to precisely map neuronal subtypes and predict their network contributions in the brain.

The formation of functional connectivity in the nervous system is governed by axon guidance that instructs nerve growth and branching during development, implying a similarity between neuronal subtypes in terms of nerve extension. We demonstrate the molecular mechanism of another layer of complexity in vertebrates by defining a transcriptional program underlying growth differences between positionally different neurons. The rate of axon extension of the early subset of embryonic dorsal root ganglion sensory neurons is encoded in neurons at different axial levels. This code is determined by a segmental pattern of axial levels of Runx family transcription factor Runx3. Runx3 in turn determines transcription levels of genes encoding cytoskeletal proteins involved in axon extension, including Rock1 and Rock2 which have ongoing activities determining axon growth in early sensory neurons and blocking Rock activity reverses axon extension deficits of Runx3(-/-) neurons. Thus, Runx3 acts to regulate positional differences in axon extension properties apparently without affecting nerve guidance and branching, a principle that could be relevant to other parts of the nervous system.

The cell cycle progression in mouse embryonic stem cells (mESCs) is controlled by ion fluxes that alter cell volume [1]. This suggests that ion fluxes might control dynamic changes in morphology over the cell cycle, such as rounding up of the cell at mitosis. However, specific channels regulating such dynamic changes and the possible interactions with actomyosin complex have not been clearly identified. Following RNAseq transcriptome analysis of cell cycle sorted mESCs, we found that expression of the K(+) ion channel Erg1 peaked in G1 cell cycle phase, which was confirmed by immunostaining. Inhibition of Erg channel activity caused loss of G1 phase cells via non-apoptotic cell death. Cells first lost the ability of membrane blebbing, a typical feature of cultured embryonic stem cells. Continued Erg inhibition further increased cell volume and the cell eventually ruptured. In addition, atomic force measurements on live cells revealed a decreased cortical stiffness after treatment, suggesting alterations in actomyosin organization. When the intracellular osmotic pressure was experimentally decreased by hypertonic solution or block of K(+) ion import via the Na, K-ATPase, cell viability was restored and cells acquired normal volume and blebbing activity. Our results suggest that Erg channels have a critical function in K(+) ion homeostasis of mESCs over the cell cycle, and that cell death following Erg inhibition is a consequence of the inability to regulate cell volume.

The generation of induced pluripotent stem (iPS) cells presents a challenge to normal developmental processes. The low efficiency and heterogeneity of most methods have hindered understanding of the precise molecular mechanisms promoting, and roadblocks preventing, efficient reprogramming. Although several intermediate populations have been described, it has proved difficult to characterize the rare, asynchronous transition from these intermediate stages to iPS cells. The rapid expansion of minor reprogrammed cells in the heterogeneous population can also obscure investigation of relevant transition processes. Understanding the biological mechanisms essential for successful iPS cell generation requires both accurate capture of cells undergoing the reprogramming process and identification of the associated global gene expression changes. Here we demonstrate that in mouse embryonic fibroblasts, reprogramming follows an orderly sequence of stage transitions, marked by changes in the cell-surface markers CD44 and ICAM1, and a Nanog-enhanced green fluorescent protein (Nanog-eGFP) reporter. RNA-sequencing analysis of these populations demonstrates two waves of pluripotency gene upregulation, and unexpectedly, transient upregulation of several epidermis-related genes, demonstrating that reprogramming is not simply the reversal of the normal developmental processes. This novel high-resolution analysis enables the construction of a detailed reprogramming route map, and the improved understanding of the reprogramming process will lead to new reprogramming strategies.

The dynamics and interactions between stem cell pools in the hair follicle (HF), sebaceous gland (SG), and interfollicular epidermis (IFE) of murine skin are still poorly understood. In this study, we used multicolor lineage tracing to mark Lgr6⁺ -expressing basal cells in the HF isthmus, SG, and IFE.We show that these Lgr6⁺ cells constitute long-term self-renewing populations within each compartment in adult skin. Quantitative analysis of clonal dynamics revealed that the Lgr6⁺ progenitor cells compete neutrally in the IFE, isthmus, and SG, indicating population asymmetry as the underlying mode of tissue renewal. Transcriptional profiling of Lgr6⁺ and Lgr6⁺ cells did not reveal a distinct Lgr6⁺ -associated gene expression signature, raising the question of whether Lgr6⁺ expression requires extrinsic niche signals. Our results elucidate the interrelation and behavior of Lgr6⁺ populations in the IFE, HF, and SG and suggest population asymmetry as a common mechanism for homeostasis in several epithelial skin compartments.

The cell cycle coordinates core functions such as replication and cell division. However, cell-cycle-regulated transcription in the control of non-core functions, such as cell identity maintenance through specific transcription factors (TFs) and signalling pathways remains unclear. Here, we provide a resource consisting of mapped transcriptomes in unsynchronized HeLa and U2OS cancer cells sorted for cell cycle phase by Fucci reporter expression. We developed a novel algorithm for data analysis that enables efficient visualization and data comparisons and identified cell cycle synchronization of Notch signalling and TFs associated with development. Furthermore, the cell cycle synchronizes with the circadian clock, providing a possible link between developmental transcriptional networks and the cell cycle. In conclusion we find that cell cycle synchronized transcriptional patterns are temporally compartmentalized and more complex than previously anticipated, involving genes, which control cell identity and development.

Significant heterogeneities in gene expression among individual cells are typically interrogated using single whole cell approaches. However, tissues that have highly interconnected processes, such as in the brain, present unique challenges. Single-nucleus RNA sequencing (SNS) has emerged as an alternative method of assessing a cell's transcriptome through the use of isolated nuclei. However, studies directly comparing expression data between nuclei and whole cells are lacking. Here, we have characterized nuclear and whole cell transcriptomes in mouse single neurons and provided a normalization strategy to reduce method-specific differences related to the length of genic regions. We confirmed a high concordance between nuclear and whole cell transcriptomes in the expression of cell type and metabolic modeling markers, but less so for a subset of genes associated with mitochondrial respiration. Therefore, our results indicate that single-nucleus transcriptome sequencing provides an effective means to profile cell type expression dynamics in previously inaccessible tissues.

Although cardinal cortical interneuron identity is established upon cell-cycle exit, it remains unclear whether specific interneuron subtypes are pre-established, and if so, how their identity is maintained prior to circuit integration. We conditionally removed Sox6 (Sox6-cKO) in migrating somatostatin (Sst+) interneurons and assessed the effects on their mature identity. In adolescent mice, five of eight molecular Sst+ subtypes were nearly absent in the Sox6-cKO cortex without a reduction in cell number. Sox6-cKO cells displayed electrophysiological maturity and expressed genes enriched within the broad class of Sst+ interneurons. Furthermore, we could infer subtype identity prior to cortical integration (embryonic day 18.5), suggesting that the loss in subtype was due to disrupted subtype maintenance. Conversely, Sox6 removal at postnatal day 7 did not disrupt marker expression in the mature cortex. Therefore, Sox6 is necessary during migration for maintenance of Sst+ subtype identity, indicating that subtype maintenance requires active transcriptional programs.

The primary motor cortex (M1) is essential for voluntary fine-motor control and is functionally conserved across mammals1. Here, using high-throughput transcriptomic and epigenomic profiling of more than 450,000 single nuclei in humans, marmoset monkeys and mice, we demonstrate a broadly conserved cellular makeup of this region, with similarities that mirror evolutionary distance and are consistent between the transcriptome and epigenome. The core conserved molecular identities of neuronal and non-neuronal cell types allow us to generate a cross-species consensus classification of cell types, and to infer conserved properties of cell types across species. Despite the overall conservation, however, many species-dependent specializations are apparent, including differences in cell-type proportions, gene expression, DNA methylation and chromatin state. Few cell-type marker genes are conserved across species, revealing a short list of candidate genes and regulatory mechanisms that are responsible for conserved features of homologous cell types, such as the GABAergic chandelier cells. This consensus transcriptomic classification allows us to use patch-seq (a combination of whole-cell patch-clamp recordings, RNA sequencing and morphological characterization) to identify corticospinal Betz cells from layer 5 in non-human primates and humans, and to characterize their highly specialized physiology and anatomy. These findings highlight the robust molecular underpinnings of cell-type diversity in M1 across mammals, and point to the genes and regulatory pathways responsible for the functional identity of cell types and their species-specific adaptations.

Understanding any brain circuit will require a categorization of its constituent neurons. In hippocampal area CA1, at least 23 classes of GABAergic neuron have been proposed to date. However, this list may be incomplete; additionally, it is unclear whether discrete classes are sufficient to describe the diversity of cortical inhibitory neurons or whether continuous modes of variability are also required. We studied the transcriptomes of 3,663 CA1 inhibitory cells, revealing 10 major GABAergic groups that divided into 49 fine-scale clusters. All previously described and several novel cell classes were identified, with three previously described classes unexpectedly found to be identical. A division into discrete classes, however, was not sufficient to describe the diversity of these cells, as continuous variation also occurred between and within classes. Latent factor analysis revealed that a single continuous variable could predict the expression levels of several genes, which correlated similarly with it across multiple cell types. Analysis of the genes correlating with this variable suggested it reflects a range from metabolically highly active faster-spiking cells that proximally target pyramidal cells to slower-spiking cells targeting distal dendrites or interneurons. These results elucidate the complexity of inhibitory neurons in one of the simplest cortical structures and show that characterizing these cells requires continuous modes of variation as well as discrete cell classes.