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.

Prototypical pacemaker neurons interact with the resident microbiota.

Pacemaker neurons exert control over neuronal circuit function by their intrinsic ability to generate rhythmic bursts of action potential. Recent work has identified rhythmic gut contractions in human, mice, and hydra to be dependent on both neurons and the resident microbiota. However, little is known about the evolutionary origin of these neurons and their interaction with microbes. In this study, we identified and functionally characterized prototypical ANO/SCN/TRPM ion channel-expressing pacemaker cells in the basal metazoan Hydra by using a combination of single-cell transcriptomics, immunochemistry, and functional experiments. Unexpectedly, these prototypical pacemaker neurons express a rich set of immune-related genes mediating their interaction with the microbial environment. Furthermore, functional experiments gave a strong support to a model of the evolutionary emergence of pacemaker cells as neurons using components of innate immunity to interact with the microbial environment and ion channels to generate rhythmic contractions.

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.

Spatially resolved multiomics on the neuronal effects induced by spaceflight in mice.

Impairment of the central nervous system (CNS) poses a significant health risk for astronauts during long-duration space missions. In this study, we employed an innovative approach by integrating single-cell multiomics (transcriptomics and chromatin accessibility) with spatial transcriptomics to elucidate the impact of spaceflight on the mouse brain in female mice. Our comparative analysis between ground control and spaceflight-exposed animals revealed significant alterations in essential brain processes including neurogenesis, synaptogenesis and synaptic transmission, particularly affecting the cortex, hippocampus, striatum and neuroendocrine structures. Additionally, we observed astrocyte activation and signs of immune dysfunction. At the pathway level, some spaceflight-induced changes in the brain exhibit similarities with neurodegenerative disorders, marked by oxidative stress and protein misfolding. Our integrated spatial multiomics approach serves as a stepping stone towards understanding spaceflight-induced CNS impairments at the level of individual brain regions and cell types, and provides a basis for comparison in future spaceflight studies. For broader scientific impact, all datasets from this study are available through an interactive data portal, as well as the National Aeronautics and Space Administration (NASA) Open Science Data Repository (OSDR).

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.

High cardiomyocyte diversity in human early prenatal heart development.

Cardiomyocytes play key roles during cardiogenesis, but have poorly understood features, especially in prenatal stages. Here, we characterized human prenatal cardiomyocytes, 6.5-7 weeks post-conception, by integrating single-cell RNA sequencing, spatial transcriptomics, and ligand-receptor interaction information. Using a computational workflow developed to dissect cell type heterogeneity, localize cell types, and explore their molecular interactions, we identified eight types of developing cardiomyocyte, more than double compared to the ones identified in the Human Developmental Cell Atlas. These have high variability in cell cycle activity, mitochondrial content, and connexin gene expression, and are differentially distributed in the ventricles, including outflow tract, and atria, including sinoatrial node. Moreover, cardiomyocyte ligand-receptor crosstalk is mainly with non-cardiomyocyte cell types, encompassing cardiogenesis-related pathways. Thus, early prenatal human cardiomyocytes are highly heterogeneous and develop unique location-dependent properties, with complex ligand-receptor crosstalk. Further elucidation of their developmental dynamics may give rise to new therapies.

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.

A small noncoding RNA links ribosome recovery and translation control to dedifferentiation during salamander limb regeneration.

Building a blastema from the stump is a key step of salamander limb regeneration. Stump-derived cells temporarily suspend their identity as they contribute to the blastema by a process generally referred to as dedifferentiation. Here, we provide evidence for a mechanism that involves an active inhibition of protein synthesis during blastema formation and growth. Relieving this inhibition results in a higher number of cycling cells and enhances the pace of limb regeneration. By small RNA profiling and fate mapping of skeletal muscle progeny as a cellular model for dedifferentiation, we find that the downregulation of miR-10b-5p is critical for rebooting the translation machinery. miR-10b-5p targets ribosomal mRNAs, and its artificial upregulation causes decreased blastema cell proliferation, reduction in transcripts that encode ribosomal subunits, diminished nascent protein synthesis, and retardation of limb regeneration. Taken together, our data identify a link between miRNA regulation, ribosome biogenesis, and protein synthesis during newt limb regeneration.

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.