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.

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.

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.

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.

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.

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.

Smart-seq2 for sensitive full-length transcriptome profiling in single cells.

Single-cell gene expression analyses hold promise for characterizing cellular heterogeneity, but current methods compromise on either the coverage, the sensitivity or the throughput. Here, we introduce Smart-seq2 with improved reverse transcription, template switching and preamplification to increase both yield and length of cDNA libraries generated from individual cells. Smart-seq2 transcriptome libraries have improved detection, coverage, bias and accuracy compared to Smart-seq libraries and are generated with off-the-shelf reagents at lower cost.

A reference transcriptome and inferred proteome for the salamander Notophthalmus viridescens.

Salamanders have a remarkable capacity to regenerate complex tissues, such as limbs and brain, and are therefore an important comparative model system for regenerative medicine. Despite these unique properties among adult vertebrates, the genomic information for amphibians in general, and salamanders in particular, is scarce. Here, we used massive parallel sequencing to reconstruct a de novo reference transcriptome of the red spotted newt (Notophthalmus viridescens) containing 118,893 transcripts with a N50 length of 2016 nts. Comparisons to other vertebrates revealed a newt transcriptome that is comparable in size and characteristics to well-annotated vertebrate transcriptomes. Identification of putative open reading frames (ORFs) enabled us to infer a comprehensive proteome, including the annotation of 19,903 newt proteins. We used the identified domain architectures (DAs) to assign ORFs phylogenetic positions, which also revealed putative salamander specific proteins. The reference transcriptome and inferred proteome of the red spotted newt will facilitate the use of systematic genomic technologies for regeneration studies in salamanders and enable evolutionary analyses of vertebrate regeneration at the molecular level.

Arrangements in the modular evolution of proteins.

It has been known for the last couple of decades that proteins evolve partly through rearrangements of larger fragments, typically domains. These units are considered the basic modules of protein structure, evolution and function. In the last few years, the analysis of protein-domain rearrangements has provided us with functional and evolutionary insights and has aided improved functional predictions and domain assignments to previously uncharacterised genes and proteins. Although some mechanisms that govern modular rearrangements of protein domains have been uncovered, such as the addition or deletion of a single N- or C-terminal domain, much is still unknown about the genetics behind these arrangements.

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.

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.

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.

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.

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.

Quantitative assessment of the structural bias in protein-protein interaction assays.

With recent publications of several large-scale protein-protein interaction (PPI) studies, the realization of the full yeast interaction network is getting closer. Here, we have analysed several yeast protein interaction datasets to understand their strengths and weaknesses. In particular, we investigate the effect of experimental biases on some of the protein properties suggested to be enriched in highly connected proteins. Finally, we use support vector machines (SVM) to assess the contribution of these properties to protein interactivity. We find that protein abundance is the most important factor for detecting interactions in tandem affinity purifications (TAP), while it is of less importance for Yeast Two Hybrid (Y2H) screens. Consequently, sequence conservation and/or essentiality of hubs may be related to their high abundance. Further, proteins with disordered structure are over-represented in Y2H screens and in one, but not the other, large-scale TAP assay. Hence, disordered regions may be important both in transient interactions and interactions in complexes. Finally, a few domain families seem to be responsible for a large part of all interactions. Most importantly, we show that there are method-specific biases in PPI experiments. Thus, care should be taken before drawing strong conclusions based on a single dataset.