PD-1 regulates ILC3-driven intestinal immunity and homeostasis.

Interleukin-(IL) 22 production by intestinal group 3 innate lymphoid cells (ILC3) is critical to maintain gut homeostasis. However, IL-22 needs to be tightly controlled; reduced IL-22 expression is associated with intestinal epithelial barrier defect while its overexpression promotes tumor development. Here, using a single-cell ribonucleic acid sequencing approach, we identified a core set of genes associated with increased IL-22 production by ILC3. Among these genes, programmed cell death 1 (PD-1), extensively studied in the context of cancer and chronic infection, was constitutively expressed on a subset of ILC3. These cells, found in the crypt of the small intestine and colon, displayed superior capacity to produce IL-22. PD-1 expression on ILC3 was dependent on the microbiota and was induced during inflammation in response to IL-23 but, conversely, was reduced in the presence of Notch ligand. PD-1+ ILC3 exhibited distinct metabolic activity with increased glycolytic, lipid, and polyamine synthesis associated with augmented proliferation compared with their PD-1- counterparts. Further, PD-1+ ILC3 showed increased expression of mitochondrial antioxidant proteins which enable the cells to maintain their levels of reactive oxygen species. Loss of PD-1 signaling in ILC3 led to reduced IL-22 production in a cell-intrinsic manner. During inflammation, PD-1 expression was increased on natural cytotoxicity receptor (NCR)- ILC3 while deficiency in PD-1 expression resulted in increased susceptibility to experimental colitis and failure to maintain gut barrier integrity. Collectively, our findings uncover a new function of the PD-1 and highlight the role of PD-1 signaling in the maintenance of gut homeostasis mediated by ILC3 in mice.

A risk-reward examination of sample multiplexing reagents for single cell RNA-Seq.

Single-cell RNA sequencing (scRNA-Seq) has emerged as a powerful tool for understanding cellular heterogeneity and function. However the choice of sample multiplexing reagents can impact data quality and experimental outcomes. In this study, we compared various multiplexing reagents, including MULTI-Seq, Hashtag antibody, and CellPlex, across diverse sample types such as human peripheral blood mononuclear cells (PBMCs), mouse embryonic brain and patient-derived xenografts (PDXs). We found that all multiplexing reagents worked well in cell types robust to ex vivo manipulation but suffered from signal-to-noise issues in more delicate sample types. We compared multiple demultiplexing algorithms which differed in performance depending on data quality. We find that minor improvements to laboratory workflows such as titration and rapid processing are critical to optimal performance. We also compared the performance of fixed scRNA-Seq kits and highlight the advantages of the Parse Biosciences kit for fragile samples. Highly multiplexed scRNA-Seq experiments require more sequencing resources, therefore we evaluated CRISPR-based destruction of non-informative genes to enhance sequencing value. Our comprehensive analysis provides insights into the selection of appropriate sample multiplexing reagents and protocols for scRNA-Seq experiments, facilitating more accurate and cost-effective studies.

Late fetal hematopoietic failure results from ZBTB11 deficiency despite abundant HSC specification.

Hematopoiesis produces diverse blood cell lineages to meet the basal needs and sudden demands of injury or infection. A rapid response to such challenges requires the expansion of specific lineages and a prompt return to balanced steady-state levels, necessitating tightly coordinated regulation. Previously we identified a requirement for the zinc finger and broad complex, tramtrak, bric-a-brac domain-containing 11 (ZBTB11) transcription factor in definitive hematopoiesis using a forward genetic screen for zebrafish myeloid mutants. To understand its relevance to mammalian systems, we extended these studies to mice. When Zbtb11 was deleted in the hematopoietic compartment, embryos died at embryonic day (E) 18.5 with hematopoietic failure. Zbtb11 hematopoietic knockout (Zbtb11hKO) hematopoietic stem cells (HSCs) were overabundantly specified from E14.5 to E17.5 compared with those in controls. Overspecification was accompanied by loss of stemness, inability to differentiate into committed progenitors and mature lineages in the fetal liver, failure to seed fetal bone marrow, and total hematopoietic failure. The Zbtb11hKO HSCs did not proliferate in vitro and were constrained in cell cycle progression, demonstrating the cell-intrinsic role of Zbtb11 in proliferation and cell cycle regulation in mammalian HSCs. Single-cell RNA sequencing analysis identified that Zbtb11-deficient HSCs were underrepresented in an erythroid-primed subpopulation and showed downregulation of oxidative phosphorylation pathways and dysregulation of genes associated with the hematopoietic niche. We identified a cell-intrinsic requirement for Zbtb11-mediated gene regulatory networks in sustaining a pool of maturation-capable HSCs and progenitor cells.

A three-stage developmental pathway for human Vγ9Vδ2 T cells within the postnatal thymus.

Vγ9Vδ2 T cells are the largest population of γδ T cells in adults and can play important roles in providing effective immunity against cancer and infection. Many studies have suggested that peripheral Vγ9Vδ2 T cells are derived from the fetal liver and thymus and that the postnatal thymus plays little role in the development of these cells. More recent evidence suggested that these cells may also develop postnatally in the thymus. Here, we used high-dimensional flow cytometry, transcriptomic analysis, functional assays, and precursor-product experiments to define the development pathway of Vγ9Vδ2 T cells in the postnatal thymus. We identify three distinct stages of development for Vγ9Vδ2 T cells in the postnatal thymus that are defined by the progressive acquisition of functional potential and major changes in the expression of transcription factors, chemokines, and other surface markers. Furthermore, our analysis of donor-matched thymus and blood revealed that the molecular requirements for the development of functional Vγ9Vδ2 T cells are delivered predominantly by the postnatal thymus and not in the periphery. Tbet and Eomes, which are required for IFN-γ and TNFα expression, are up-regulated as Vγ9Vδ2 T cells mature in the thymus, and mature thymic Vγ9Vδ2 T cells rapidly express high levels of these cytokines after stimulation. Similarly, the postnatal thymus programs Vγ9Vδ2 T cells to express the cytolytic molecules, perforin, granzyme A, and granzyme K. This study provides a greater understanding of how Vγ9Vδ2 T cells develop in humans and may lead to opportunities to manipulate these cells to treat human diseases.

Zinc Supplementation with or without Additional Micronutrients Does Not Affect Peripheral Blood Gene Expression or Serum Cytokine Level in Bangladeshi Children.

Preventive zinc supplementation provided as a stand-alone dispersible tablet, or via home fortification as multiple micronutrient powders (MNPs), has been considered a potential strategy to prevent zinc deficiency and improve health (including immune) outcomes among children in low- and middle-income countries. However, the impact of zinc supplementation on immune profiles has not been well characterized. We sought to define the effect of zinc supplementation on peripheral blood gene expression and cytokine levels among young children in Dhaka, Bangladesh. In a sub-study of a large randomized, controlled, community-based efficacy trial where children 9-11 months of age received one of the following interventions on a daily basis for 24 weeks: (1) MNPs containing 10 mg of zinc; (2) dispersible tablet containing 10 mg zinc; or (3) placebo powder, we used RNA sequencing to profile the peripheral blood gene expression, as well as highly sensitive multiplex assays to detect cytokine profiles. We profiled samples from 100 children enrolled in the parent trial (zinc MNPs 28, zinc tablets 39, placebo 33). We did not detect an effect from either zinc intervention on differential peripheral blood gene expression at the end of the intervention, or an effect from the intervention on changes in gene expression from baseline. We also did not detect an effect from either intervention on cytokine concentrations. Exploratory analysis did not identify an association between undernutrition (defined as stunting, underweight or wasting) and peripheral blood gene expression. Zinc interventions in children did not produce a gene expression or cytokine signature in the peripheral blood. However, this study demonstrates a proof of principle that sensitive multi-omic techniques can be applied to samples collected in field studies.

The histone acetyltransferase HBO1 promotes efficient tip cell sprouting during angiogenesis.

Blood vessel growth and remodelling are essential during embryonic development and disease pathogenesis. The diversity of endothelial cells (ECs) is transcriptionally evident and ECs undergo dynamic changes in gene expression during vessel growth and remodelling. Here, we investigated the role of the histone acetyltransferase HBO1 (KAT7), which is important for activating genes during development and for histone H3 lysine 14 acetylation (H3K14ac). Loss of HBO1 and H3K14ac impaired developmental sprouting angiogenesis and reduced pathological EC overgrowth in the retinal endothelium. Single-cell RNA sequencing of retinal ECs revealed an increased abundance of tip cells in Hbo1-deficient retinas, which led to EC overcrowding in the retinal sprouting front and prevented efficient tip cell migration. We found that H3K14ac was highly abundant in the endothelial genome in both intra- and intergenic regions, suggesting that HBO1 acts as a genome organiser that promotes efficient tip cell behaviour necessary for sprouting angiogenesis. This article has an associated 'The people behind the papers' interview.

Clonal multi-omics reveals Bcor as a negative regulator of emergency dendritic cell development.

Despite advances in single-cell multi-omics, a single stem or progenitor cell can only be tested once. We developed clonal multi-omics, in which daughters of a clone act as surrogates of the founder, thereby allowing multiple independent assays per clone. With SIS-seq, clonal siblings in parallel "sister" assays are examined either for gene expression by RNA sequencing (RNA-seq) or for fate in culture. We identified, and then validated using CRISPR, genes that controlled fate bias for different dendritic cell (DC) subtypes. This included Bcor as a suppressor of plasmacytoid DC (pDC) and conventional DC type 2 (cDC2) numbers during Flt3 ligand-mediated emergency DC development. We then developed SIS-skew to examine development of wild-type and Bcor-deficient siblings of the same clone in parallel. We found Bcor restricted clonal expansion, especially for cDC2s, and suppressed clonal fate potential, especially for pDCs. Therefore, SIS-seq and SIS-skew can reveal the molecular and cellular mechanisms governing clonal fate.

Single-cell analyses reveal the clonal and molecular aetiology of Flt3L-induced emergency dendritic cell development.

Regulation of haematopoietic stem and progenitor cell (HSPC) fate is crucial during homeostasis and under stress conditions. Here we examine the aetiology of the Flt3 ligand (Flt3L)-mediated increase of type 1 conventional dendritic cells (cDC1s). Using cellular barcoding we demonstrate this occurs through selective clonal expansion of HSPCs that are primed to produce cDC1s and not through activation of cDC1 fate by other HSPCs. In particular, multi/oligo-potent clones selectively amplify their cDC1 output, without compromising the production of other lineages, via a process we term tuning. We then develop Divi-Seq to simultaneously profile the division history, surface phenotype and transcriptome of individual HSPCs. We discover that Flt3L-responsive HSPCs maintain a proliferative 'early progenitor'-like state, leading to the selective expansion of multiple transitional cDC1-primed progenitor stages that are marked by Irf8 expression. These findings define the mechanistic action of Flt3L through clonal tuning, which has important implications for other models of 'emergency' haematopoiesis.

Transcriptomic Profiling of Human Pluripotent Stem Cell-derived Retinal Pigment Epithelium over Time.

Human pluripotent stem cell (hPSC)-derived progenies are immature versions of cells, presenting a potential limitation to the accurate modelling of diseases associated with maturity or age. Hence, it is important to characterise how closely cells used in culture resemble their native counterparts. In order to select appropriate time points of retinal pigment epithelium (RPE) cultures that reflect native counterparts, we characterised the transcriptomic profiles of the hPSC-derived RPE cells from 1- and 12-month cultures. We differentiated the human embryonic stem cell line H9 into RPE cells, performed single-cell RNA-sequencing of a total of 16,576 cells to assess the molecular changes of the RPE cells across these two culture time points. Our results indicate the stability of the RPE transcriptomic signature, with no evidence of an epithelial-mesenchymal transition, and with the maturing populations of the RPE observed with time in culture. Assessment of Gene Ontology pathways revealed that as the cultures age, RPE cells upregulate expression of genes involved in metal binding and antioxidant functions. This might reflect an increased ability to handle oxidative stress as cells mature. Comparison with native human RPE data confirms a maturing transcriptional profile of RPE cells in culture. These results suggest that long-term in vitro culture of RPE cells allows the modelling of specific phenotypes observed in native mature tissues. Our work highlights the transcriptional landscape of hPSC-derived RPE cells as they age in culture, which provides a reference for native and patient samples to be benchmarked against.

A new lymphoid-primed progenitor marked by Dach1 downregulation identified with single cell multi-omics.

A classical view of blood cell development is that multipotent hematopoietic stem and progenitor cells (HSPCs) become lineage-restricted at defined stages. Lin-c-Kit+Sca-1+Flt3+ cells, termed lymphoid-primed multipotent progenitors (LMPPs), have lost megakaryocyte and erythroid potential but are heterogeneous in their fate. Here, through single-cell RNA sequencing, we identify the expression of Dach1 and associated genes in this fraction as being coexpressed with myeloid/stem genes but inversely correlated with lymphoid genes. Through generation of Dach1-GFP reporter mice, we identify a transcriptionally and functionally unique Dach1-GFP- subpopulation within LMPPs with lymphoid potential with low to negligible classic myeloid potential. We term these 'lymphoid-primed progenitors' (LPPs). These findings define an early definitive branch point of lymphoid development in hematopoiesis and a means for prospective isolation of LPPs.

Membrane budding is a major mechanism of in vivo platelet biogenesis.

How platelets are produced by megakaryocytes in vivo remains controversial despite more than a century of investigation. Megakaryocytes readily produce proplatelet structures in vitro; however, visualization of platelet release from proplatelets in vivo has remained elusive. We show that within the native prenatal and adult environments, the frequency and rate of proplatelet formation is incompatible with the physiological demands of platelet replacement. We resolve this inconsistency by performing in-depth analysis of plasma membrane budding, a cellular process that has previously been dismissed as a source of platelet production. Our studies demonstrate that membrane budding results in the sustained release of platelets directly into the peripheral circulation during both fetal and adult life without induction of cell death or proplatelet formation. In support of this model, we demonstrate that in mice deficient for NF-E2 (the thrombopoietic master regulator), the absence of membrane budding correlates with failure of in vivo platelet production. Accordingly, we propose that membrane budding, rather than proplatelet formation, supplies the majority of the platelet biomass.

Benchmarking single cell RNA-sequencing analysis pipelines using mixture control experiments.

Single cell RNA-sequencing (scRNA-seq) technology has undergone rapid development in recent years, leading to an explosion in the number of tailored data analysis methods. However, the current lack of gold-standard benchmark datasets makes it difficult for researchers to systematically compare the performance of the many methods available. Here, we generated a realistic benchmark experiment that included single cells and admixtures of cells or RNA to create 'pseudo cells' from up to five distinct cancer cell lines. In total, 14 datasets were generated using both droplet and plate-based scRNA-seq protocols. We compared 3,913 combinations of data analysis methods for tasks ranging from normalization and imputation to clustering, trajectory analysis and data integration. Evaluation revealed pipelines suited to different types of data for different tasks. Our data and analysis provide a comprehensive framework for benchmarking most common scRNA-seq analysis steps.

scPipe: A flexible R/Bioconductor preprocessing pipeline for single-cell RNA-sequencing data.

Single-cell RNA sequencing (scRNA-seq) technology allows researchers to profile the transcriptomes of thousands of cells simultaneously. Protocols that incorporate both designed and random barcodes have greatly increased the throughput of scRNA-seq, but give rise to a more complex data structure. There is a need for new tools that can handle the various barcoding strategies used by different protocols and exploit this information for quality assessment at the sample-level and provide effective visualization of these results in preparation for higher-level analyses. To this end, we developed scPipe, an R/Bioconductor package that integrates barcode demultiplexing, read alignment, UMI-aware gene-level quantification and quality control of raw sequencing data generated by multiple protocols that include CEL-seq, MARS-seq, Chromium 10X, Drop-seq and Smart-seq. scPipe produces a count matrix that is essential for downstream analysis along with an HTML report that summarises data quality. These results can be used as input for downstream analyses including normalization, visualization and statistical testing. scPipe performs this processing in a few simple R commands, promoting reproducible analysis of single-cell data that is compatible with the emerging suite of open-source scRNA-seq analysis tools available in R/Bioconductor and beyond. The scPipe R package is available for download from https://www.bioconductor.org/packages/scPipe.

Deciphering the Innate Lymphoid Cell Transcriptional Program.

Innate lymphoid cells (ILCs) are enriched at mucosal surfaces, where they provide immune surveillance. All ILC subsets develop from a common progenitor that gives rise to pre-committed progenitors for each of the ILC lineages. Currently, the temporal control of gene expression that guides the emergence of these progenitors is poorly understood. We used global transcriptional mapping to analyze gene expression in different ILC progenitors. We identified PD-1 to be specifically expressed in PLZF+ ILCp and revealed that the timing and order of expression of the transcription factors NFIL3, ID2, and TCF-1 was critical. Importantly, induction of ILC lineage commitment required only transient expression of NFIL3 prior to ID2 and TCF-1 expression. These findings highlight the importance of the temporal program that permits commitment of progenitors to the ILC lineage, and they expand our understanding of the core transcriptional program by identifying potential regulators of ILC development.

Derivation of novel human ground state naive pluripotent stem cells.

Mouse embryonic stem (ES) cells are isolated from the inner cell mass of blastocysts, and can be preserved in vitro in a naive inner-cell-mass-like configuration by providing exogenous stimulation with leukaemia inhibitory factor (LIF) and small molecule inhibition of ERK1/ERK2 and GSK3ß signalling (termed 2i/LIF conditions). Hallmarks of naive pluripotency include driving Oct4 (also known as Pou5f1) transcription by its distal enhancer, retaining a pre-inactivation X chromosome state, and global reduction in DNA methylation and in H3K27me3 repressive chromatin mark deposition on developmental regulatory gene promoters. Upon withdrawal of 2i/LIF, naive mouse ES cells can drift towards a primed pluripotent state resembling that of the post-implantation epiblast. Although human ES cells share several molecular features with naive mouse ES cells, they also share a variety of epigenetic properties with primed murine epiblast stem cells (EpiSCs). These include predominant use of the proximal enhancer element to maintain OCT4 expression, pronounced tendency for X chromosome inactivation in most female human ES cells, increase in DNA methylation and prominent deposition of H3K27me3 and bivalent domain acquisition on lineage regulatory genes. The feasibility of establishing human ground state naive pluripotency in vitro with equivalent molecular and functional features to those characterized in mouse ES cells remains to be defined. Here we establish defined conditions that facilitate the derivation of genetically unmodified human naive pluripotent stem cells from already established primed human ES cells, from somatic cells through induced pluripotent stem (iPS) cell reprogramming or directly from blastocysts. The novel naive pluripotent cells validated herein retain molecular characteristics and functional properties that are highly similar to mouse naive ES cells, and distinct from conventional primed human pluripotent cells. This includes competence in the generation of cross-species chimaeric mouse embryos that underwent organogenesis following microinjection of human naive iPS cells into mouse morulas. Collectively, our findings establish new avenues for regenerative medicine, patient-specific iPS cell disease modelling and the study of early human development in vitro and in vivo.

Deterministic direct reprogramming of somatic cells to pluripotency.

Somatic cells can be inefficiently and stochastically reprogrammed into induced pluripotent stem (iPS) cells by exogenous expression of Oct4 (also called Pou5f1), Sox2, Klf4 and Myc (hereafter referred to as OSKM). The nature of the predominant rate-limiting barrier(s) preventing the majority of cells to successfully and synchronously reprogram remains to be defined. Here we show that depleting Mbd3, a core member of the Mbd3/NuRD (nucleosome remodelling and deacetylation) repressor complex, together with OSKM transduction and reprogramming in naive pluripotency promoting conditions, result in deterministic and synchronized iPS cell reprogramming (near 100% efficiency within seven days from mouse and human cells). Our findings uncover a dichotomous molecular function for the reprogramming factors, serving to reactivate endogenous pluripotency networks while simultaneously directly recruiting the Mbd3/NuRD repressor complex that potently restrains the reactivation of OSKM downstream target genes. Subsequently, the latter interactions, which are largely depleted during early pre-implantation development in vivo, lead to a stochastic and protracted reprogramming trajectory towards pluripotency in vitro. The deterministic reprogramming approach devised here offers a novel platform for the dissection of molecular dynamics leading to establishing pluripotency at unprecedented flexibility and resolution.

High-throughput chromatin immunoprecipitation for genome-wide mapping of in vivo protein-DNA interactions and epigenomic states.

Dynamic protein binding to DNA elements regulates genome function and cell fate. Although methods for mapping in vivo protein-DNA interactions are becoming crucial for every aspect of genomic research, they are laborious and costly, thereby limiting progress. Here we present a protocol for mapping in vivo protein-DNA interactions using a high-throughput chromatin immunoprecipitation (HT-ChIP) approach. By using paramagnetic beads, we streamline the entire ChIP and indexed library construction process: sample transfer and loss is minimized and the need for manually labor-intensive procedures such as washes, gel extraction and DNA precipitation is eliminated. All of this allows for fully automated, cost effective and highly sensitive 96-well ChIP sequencing (ChIP-seq). Sample preparation takes 3 d from cultured cells to pooled libraries. Compared with previous methods, HT-ChIP is more suitable for large-scale in vivo studies, specifically those measuring the dynamics of a large number of different chromatin modifications/transcription factors or multiple perturbations.