Molecular and spatial landmarks of early mouse skin development.

A wealth of specialized cell populations within the skin facilitates its hair-producing, protective, sensory, and thermoregulatory functions. How the vast cell-type diversity and tissue architecture develops is largely unexplored. Here, with single-cell transcriptomics, spatial cell-type assignment, and cell-lineage tracing, we deconstruct early embryonic mouse skin during the key transitions from seemingly uniform developmental precursor states to a multilayered, multilineage epithelium, and complex dermal identity. We identify the spatiotemporal emergence of hair-follicle-inducing, muscle-supportive, and fascia-forming fibroblasts. We also demonstrate the formation of the panniculus carnosus muscle (PCM), sprouting blood vessels without pericyte coverage, and the earliest residence of mast and dendritic immune cells in skin. Finally, we identify an unexpected epithelial heterogeneity within the early single-layered epidermis and a signaling-rich periderm layer. Overall, this cellular and molecular blueprint of early skin development-which can be explored at https://kasperlab.org/tools-establishes histological landmarks and highlights unprecedented dynamic interactions among skin cells.

Single-Cell RNA Analysis Reveals Cell-Intrinsic Functions of CAR T Cells Correlating with Response in a Phase II Study of Lymphoma Patients.

PURPOSE: Although CD19 chimeric antigen receptor T cells (CAR-T) therapy has shown remarkable success in B-cell malignancies, a substantial fraction of patients do not obtain a long-term clinical response. This could be influenced by the quality of the individual CAR-T infusion product. To shed some light on this, clinical outcome was correlated to characteristics of CAR-T infusion products. PATIENTS AND METHODS: In this phase II study, patients with B-cell lymphoma (n = 23) or leukemia (n = 1) received one or two infusions of third-generation CD19-directed CAR-Ts (2 × 108/m2). The clinical trial was registered at clinicaltrials.gov: NCT03068416. We investigated the transcriptional profile of individual CD19 CAR-T infusion products using targeted single-cell RNA sequencing and multicolor flow cytometry. RESULTS: Two CAR-T infusions were not better than one in the settings used in this study. As for the CAR-T infusion products, we found that effector-like CD8+CAR-Ts with a high polyfunctionality, high cytotoxic and cytokine production profile, and low dysfunctional signature were associated with clinical response. An extended ex vivo expansion time during CAR-T manufacturing negatively influenced the proportion of effector CD8+CAR-Ts in the infusion product. CONCLUSIONS: We identified cell-intrinsic characteristics of effector CD8+CAR-Ts correlating with response that could be used as an indicator for clinical outcome. The results in the study also serve as a guide to CAR-T manufacturing practices.

Single-cell multi-omics of human preimplantation embryos shows susceptibility to glucocorticoids.

The preconceptual, intrauterine, and early life environments can have a profound and long-lasting impact on the developmental trajectories and health outcomes of the offspring. Given the relatively low success rates of assisted reproductive technologies (ART; ∼25%), additives and adjuvants, such as glucocorticoids, are used to improve the success rate. Considering the dynamic developmental events that occur during this window, these exposures may alter blastocyst formation at a molecular level, and as such, affect not only the viability of the embryo and the ability of the blastocyst to implant, but also the developmental trajectory of the first three cell lineages, ultimately influencing the physiology of the embryo. In this study, we present a comprehensive single-cell transcriptome, methylome, and small RNA atlas in the day 7 human embryo. We show that, despite no change in morphology and developmental features, preimplantation glucocorticoid exposure reprograms the molecular profile of the TE lineage, and these changes are associated with an altered metabolic and inflammatory response. Our data also suggest that glucocorticoids can precociously mature the TE sublineages, supported by the presence of extravillous trophoblast markers in the polar sublineage and presence of X Chromosome dosage compensation. Further, we have elucidated that epigenetic regulation-DNA methylation and microRNAs (miRNAs)-likely underlies the transcriptional changes observed. This study suggests that exposures to exogenous compounds during preimplantation may unintentionally reprogram the human embryo, possibly leading to suboptimal development and longer-term health outcomes.

Limited access to antigen drives generation of early B cell memory while restraining the plasmablast response.

Cell fate decisions during early B cell activation determine the outcome of responses to pathogens and vaccines. We examined the early B cell response to T-dependent antigen in mice by single-cell RNA sequencing. Early after immunization, a homogeneous population of activated precursors (APs) gave rise to a transient wave of plasmablasts (PBs), followed a day later by the emergence of germinal center B cells (GCBCs). Most APs rapidly exited the cell cycle, giving rise to non-GC-derived early memory B cells (eMBCs) that retained an AP-like transcriptional profile. Rapid decline of antigen availability controlled these events; provision of excess antigen precluded cell cycle exit and induced a new wave of PBs. Fate mapping revealed a prominent contribution of eMBCs to the MBC pool. Quiescent cells with an MBC phenotype dominated the early response to immunization in primates. A reservoir of APs/eMBCs may enable rapid readjustment of the immune response when failure to contain a threat is manifested by increased antigen availability.

Characterization of the nuclear and cytosolic transcriptomes in human brain tissue reveals new insights into the subcellular distribution of RNA transcripts.

Transcriptome analysis has mainly relied on analyzing RNA sequencing data from whole cells, overlooking the impact of subcellular RNA localization and its influence on our understanding of gene function, and interpretation of gene expression signatures in cells. Here, we separated cytosolic and nuclear RNA from human fetal and adult brain samples and performed a comprehensive analysis of cytosolic and nuclear transcriptomes. There are significant differences in RNA expression for protein-coding and lncRNA genes between cytosol and nucleus. We show that transcripts encoding the nuclear-encoded mitochondrial proteins are significantly enriched in the cytosol compared to the rest of protein-coding genes. Differential expression analysis between fetal and adult frontal cortex show that results obtained from the cytosolic RNA differ from results using nuclear RNA both at the level of transcript types and the number of differentially expressed genes. Our data provide a resource for the subcellular localization of thousands of RNA transcripts in the human brain and highlight differences in using the cytosolic or the nuclear transcriptomes for expression analysis.

Tissue-specific transcriptional imprinting and heterogeneity in human innate lymphoid cells revealed by full-length single-cell RNA-sequencing.

The impact of the microenvironment on innate lymphoid cell (ILC)-mediated immunity in humans remains largely unknown. Here we used full-length Smart-seq2 single-cell RNA-sequencing to unravel tissue-specific transcriptional profiles and heterogeneity of CD127+ ILCs across four human tissues. Correlation analysis identified gene modules characterizing the migratory properties of tonsil and blood ILCs, and signatures of tissue-residency, activation and modified metabolism in colon and lung ILCs. Trajectory analysis revealed potential differentiation pathways from circulating and tissue-resident naïve ILCs to a spectrum of mature ILC subsets. In the lung we identified both CRTH2+ and CRTH2- ILC2 with lung-specific signatures, which could be recapitulated by alarmin-exposure of circulating ILC2. Finally, we describe unique TCR-V(D)J-rearrangement patterns of blood ILC1-like cells, revealing a subset of potentially immature ILCs with TCR-δ rearrangement. Our study provides a useful resource for in-depth understanding of ILC-mediated immunity in humans, with implications for disease.

Author Correction: Single cell transcriptomics identifies stem cell-derived graft composition in a model of Parkinson's disease.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Spatio-molecular domains identified in the mouse subthalamic nucleus and neighboring glutamatergic and GABAergic brain structures.

The subthalamic nucleus (STN) is crucial for normal motor, limbic and associative function. STN dysregulation is correlated with several brain disorders, including Parkinson's disease and obsessive compulsive disorder (OCD), for which high-frequency stimulation of the STN is increasing as therapy. However, clinical progress is hampered by poor knowledge of the anatomical-functional organization of the STN. Today, experimental mouse genetics provides outstanding capacity for functional decoding, provided selective promoters are available. Here, we implemented single-nuclei RNA sequencing (snRNASeq) of the mouse STN followed through with histological analysis of 16 candidate genes of interest. Our results demonstrate that the mouse STN is composed of at least four spatio-molecularly defined domains, each distinguished by defined sets of promoter activities. Further, molecular profiles dissociate the STN from the adjoining para-STN (PSTN) and neighboring structures of the hypothalamus, mammillary nuclei and zona incerta. Enhanced knowledge of STN´s internal organization should prove useful towards genetics-based functional decoding of this clinically relevant brain structure.

Single cell transcriptomics identifies stem cell-derived graft composition in a model of Parkinson's disease.

Cell replacement is a long-standing and realistic goal for the treatment of Parkinson's disease (PD). Cells for transplantation can be obtained from fetal brain tissue or from stem cells. However, after transplantation, dopamine (DA) neurons are seen to be a minor component of grafts, and it has remained difficult to determine the identity of other cell types. Here, we report analysis by single-cell RNA sequencing (scRNA-seq) combined with comprehensive histological analyses to characterize intracerebral grafts from human embryonic stem cells (hESCs) and fetal tissue after functional maturation in a pre-clinical rat PD model. We show that neurons and astrocytes are major components in both fetal and stem cell-derived grafts. Additionally, we identify a cell type closely resembling a class of recently identified perivascular-like cells in stem cell-derived grafts. Thus, this study uncovers previously unknown cellular diversity in a clinically relevant cell replacement PD model.

Conbase: a software for unsupervised discovery of clonal somatic mutations in single cells through read phasing.

Accurate variant calling and genotyping represent major limiting factors for downstream applications of single-cell genomics. Here, we report Conbase for the identification of somatic mutations in single-cell DNA sequencing data. Conbase leverages phased read data from multiple samples in a dataset to achieve increased confidence in somatic variant calls and genotype predictions. Comparing the performance of Conbase to three other methods, we find that Conbase performs best in terms of false discovery rate and specificity and provides superior robustness on simulated data, in vitro expanded fibroblasts and clonal lymphocyte populations isolated directly from a healthy human donor.

Single-cell RNA sequencing reveals midbrain dopamine neuron diversity emerging during mouse brain development.

Midbrain dopamine (mDA) neurons constitute a heterogenous group of cells that have been intensely studied, not least because their degeneration causes major symptoms in Parkinson's disease. Understanding the diversity of mDA neurons - previously well characterized anatomically - requires a systematic molecular classification at the genome-wide gene expression level. Here, we use single cell RNA sequencing of isolated mouse neurons expressing the transcription factor Pitx3, a marker for mDA neurons. Analyses include cells isolated during development up until adulthood and the results are validated by histological characterization of newly identified markers. This identifies seven neuron subgroups divided in two major branches of developing Pitx3-expressing neurons. Five of them express dopaminergic markers, while two express glutamatergic and GABAergic markers, respectively. Analysis also indicate evolutionary conservation of diversity in humans. This comprehensive molecular characterization will provide a valuable resource for elucidating mDA neuron subgroup development and function in the mammalian brain.

Author Correction: A comprehensive map coupling histone modifications with gene regulation in adult dopaminergic and serotonergic neurons.

In the originally published version of this Article, financial support was not fully acknowledged. The PDF and HTML versions of the Article have now been corrected to include support to Thomas Perlmann provided by Knut and Alice Wallenberg Foundation (grant 2013.0075) and Swedish Research Council (VR; grant 2016-02506).

Sequentially acting SOX proteins orchestrate astrocyte- and oligodendrocyte-specific gene expression.

SOX transcription factors have important roles during astrocyte and oligodendrocyte development, but how glial genes are specified and activated in a sub-lineage-specific fashion remains unknown. Here, we define glial-specific gene expression in the developing spinal cord using single-cell RNA-sequencing. Moreover, by ChIP-seq analyses we show that these glial gene sets are extensively preselected already in multipotent neural precursor cells through prebinding by SOX3. In the subsequent lineage-restricted glial precursor cells, astrocyte genes become additionally targeted by SOX9 at DNA regions strongly enriched for Nfi binding motifs. Oligodendrocyte genes instead are prebound by SOX9 only, at sites which during oligodendrocyte maturation are targeted by SOX10. Interestingly, reporter gene assays and functional studies in the spinal cord reveal that SOX3 binding represses the synergistic activation of astrocyte genes by SOX9 and NFIA, whereas oligodendrocyte genes are activated in a combinatorial manner by SOX9 and SOX10. These genome-wide studies demonstrate how sequentially expressed SOX proteins act on lineage-specific regulatory DNA elements to coordinate glial gene expression both in a temporal and in a sub-lineage-specific fashion.

Exploring parasite heterogeneity using single-cell RNA-seq reveals a gene signature among sexual stage Plasmodium falciparum parasites.

The malaria parasite has a complex lifecycle, including several events of differentiation and stage progression, while actively evading immunity in both its mosquito and human hosts. Important parasite gene expression and regulation during these events remain hidden in rare populations of cells. Here, we combine a capillary-based platform for cell isolation with single-cell RNA-sequencing to transcriptionally profile 165 single infected red blood cells (iRBCs) during the intra-erythrocytic developmental cycle (IDC). Unbiased analyses of single-cell data grouped the cells into eight transcriptional states during IDC. Interestingly, we uncovered a gene signature from the single iRBC analyses that can successfully discriminate between developing asexual and sexual stage parasites at cellular resolution, and we verify five, previously undefined, gametocyte stage specific genes. Moreover, we show the capacity of detecting expressed genes from the variable gene families in single parasites, despite the sparse nature of data. In total, the single parasite transcriptomics holds promise for molecular dissection of rare parasite phenotypes throughout the malaria lifecycle.

A comprehensive map coupling histone modifications with gene regulation in adult dopaminergic and serotonergic neurons.

The brain is composed of hundreds of different neuronal subtypes, which largely retain their identity throughout the lifespan of the organism. The mechanisms governing this stability are not fully understood, partly due to the diversity and limited size of clinically relevant neuronal populations, which constitute a technical challenge for analysis. Here, using a strategy that allows for ChIP-seq combined with RNA-seq in small neuronal populations in vivo, we present a comparative analysis of permissive and repressive histone modifications in adult midbrain dopaminergic neurons, raphe nuclei serotonergic neurons, and embryonic neural progenitors. Furthermore, we utilize the map generated by our analysis to show that the transcriptional response of midbrain dopaminergic neurons following 6-OHDA or methamphetamine injection is characterized by increased expression of genes with promoters dually marked by H3K4me3/H3K27me3. Our study provides an in vivo genome-wide analysis of permissive/repressive histone modifications coupled to gene expression in these rare neuronal subtypes.

Prostaglandin E2 suppresses human group 2 innate lymphoid cell function.

BACKGROUND: Group 2 innate lymphoid cells (ILC2s) are involved in the initial phase of type 2 inflammation and can amplify allergic immune responses by orchestrating other type 2 immune cells. Prostaglandin (PG) E2 is a bioactive lipid that plays protective roles in the lung, particularly during allergic inflammation. OBJECTIVE: We set out to investigate how PGE2 regulates human ILC2 function. METHODS: The effects of PGE2 on human ILC2 proliferation and intracellular cytokine and transcription factor expression were assessed by means of flow cytometry. Cytokine production was measured by using ELISA, and real-time quantitative PCR was performed to detect PGE2 receptor expression. RESULTS: PGE2 inhibited GATA-3 expression, as well as production of the type 2 cytokines IL-5 and IL-13, from human tonsillar and blood ILC2s in response to stimulation with a combination of IL-25, IL-33, thymic stromal lymphopoietin, and IL-2. Furthermore, PGE2 downregulated the expression of IL-2 receptor α (CD25). In line with this observation, PGE2 decreased ILC2 proliferation. These effects were mediated by the combined action of E-type prostanoid receptor (EP) 2 and EP4 receptors, which were specifically expressed on ILC2s. CONCLUSION: Our findings reveal that PGE2 limits ILC2 activation and propose that selective EP2 and EP4 receptor agonists might serve as a promising therapeutic approach in treating allergic diseases by suppressing ILC2 function.

Neuropilin-1 Is Expressed on Lymphoid Tissue Residing LTi-like Group 3 Innate Lymphoid Cells and Associated with Ectopic Lymphoid Aggregates.

Here, we characterize a subset of ILC3s that express Neuropilin1 (NRP1) and are present in lymphoid tissues, but not in the peripheral blood or skin. NRP1+ group 3 innate lymphoid cells (ILC3s) display in vitro lymphoid tissue inducer (LTi) activity. In agreement with this, NRP1+ ILC3s are mainly located in proximity to high endothelial venules (HEVs) and express cell surface molecules involved in lymphocyte migration in secondary lymphoid tissues via HEVs. NRP1 was also expressed on mouse fetal LTi cells, indicating that NRP1 is a conserved marker for LTi cells. Human NRP1+ ILC3s are primed cells because they express CD45RO and produce higher amounts of cytokines than NRP1- cells, which express CD45RA. The NRP1 ligand vascular endothelial growth factor A (VEGF-A) served as a chemotactic factor for NRP1+ ILC3s. NRP1+ ILC3s are present in lung tissues from smokers and patients with chronic obstructive pulmonary disease, suggesting a role in angiogenesis and/or the initiation of ectopic pulmonary lymphoid aggregates.

Single-Cell Analysis Reveals a Close Relationship between Differentiating Dopamine and Subthalamic Nucleus Neuronal Lineages.

Stem cell engineering and grafting of mesencephalic dopamine (mesDA) neurons is a promising strategy for brain repair in Parkinson's disease (PD). Refinement of differentiation protocols to optimize this approach will require deeper understanding of mesDA neuron development. Here, we studied this process using transcriptome-wide single-cell RNA sequencing of mouse neural progenitors expressing the mesDA neuron determinant Lmx1a. This approach resolved the differentiation of mesDA and neighboring neuronal lineages and revealed a remarkably close relationship between developing mesDA and subthalamic nucleus (STN) neurons, while also highlighting a distinct transcription factor set that can distinguish between them. While previous hESC mesDA differentiation protocols have relied on markers that are shared between the two lineages, we found that application of these highlighted markers can help to refine current stem cell engineering protocols, increasing the proportion of appropriately patterned mesDA progenitors. Our results, therefore, have important implications for cell replacement therapy in PD.

The heterogeneity of human CD127(+) innate lymphoid cells revealed by single-cell RNA sequencing.

Innate lymphoid cells (ILCs) are increasingly appreciated as important participants in homeostasis and inflammation. Substantial plasticity and heterogeneity among ILC populations have been reported. Here we have delineated the heterogeneity of human ILCs through single-cell RNA sequencing of several hundreds of individual tonsil CD127(+) ILCs and natural killer (NK) cells. Unbiased transcriptional clustering revealed four distinct populations, corresponding to ILC1 cells, ILC2 cells, ILC3 cells and NK cells, with their respective transcriptomes recapitulating known as well as unknown transcriptional profiles. The single-cell resolution additionally divulged three transcriptionally and functionally diverse subpopulations of ILC3 cells. Our systematic comparison of single-cell transcriptional variation within and between ILC populations provides new insight into ILC biology during homeostasis, with additional implications for dysregulation of the immune system.

Tn5 transposase and tagmentation procedures for massively scaled sequencing projects.

Massively parallel DNA sequencing of thousands of samples in a single machine-run is now possible, but the preparation of the individual sequencing libraries is expensive and time-consuming. Tagmentation-based library construction, using the Tn5 transposase, is efficient for generating sequencing libraries but currently relies on undisclosed reagents, which severely limits development of novel applications and the execution of large-scale projects. Here, we present simple and robust procedures for Tn5 transposase production and optimized reaction conditions for tagmentation-based sequencing library construction. We further show how molecular crowding agents both modulate library lengths and enable efficient tagmentation from subpicogram amounts of cDNA. The comparison of single-cell RNA-sequencing libraries generated using produced and commercial Tn5 demonstrated equal performances in terms of gene detection and library characteristics. Finally, because naked Tn5 can be annealed to any oligonucleotide of choice, for example, molecular barcodes in single-cell assays or methylated oligonucleotides for bisulfite sequencing, custom Tn5 production and tagmentation enable innovation in sequencing-based applications.