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.

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.

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.

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.

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.

A high-throughput chromatin immunoprecipitation approach reveals principles of dynamic gene regulation in mammals.

Understanding the principles governing mammalian gene regulation has been hampered by the difficulty in measuring in vivo binding dynamics of large numbers of transcription factors (TF) to DNA. Here, we develop a high-throughput Chromatin ImmunoPrecipitation (HT-ChIP) method to systematically map protein-DNA interactions. HT-ChIP was applied to define the dynamics of DNA binding by 25 TFs and 4 chromatin marks at 4 time-points following pathogen stimulus of dendritic cells. Analyzing over 180,000 TF-DNA interactions we find that TFs vary substantially in their temporal binding landscapes. This data suggests a model for transcription regulation whereby TF networks are hierarchically organized into cell differentiation factors, factors that bind targets prior to stimulus to prime them for induction, and factors that regulate specific gene programs. Overlaying HT-ChIP data on gene-expression dynamics shows that many TF-DNA interactions are established prior to the stimuli, predominantly at immediate-early genes, and identified specific TF ensembles that coordinately regulate gene-induction.

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.

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.

The H3K27 demethylase Utx regulates somatic and germ cell epigenetic reprogramming.

Induced pluripotent stem cells (iPSCs) can be derived from somatic cells by ectopic expression of different transcription factors, classically Oct4 (also known as Pou5f1), Sox2, Klf4 and Myc (abbreviated as OSKM). This process is accompanied by genome-wide epigenetic changes, but how these chromatin modifications are biochemically determined requires further investigation. Here we show in mice and humans that the histone H3 methylated Lys 27 (H3K27) demethylase Utx (also known as Kdm6a) regulates the efficient induction, rather than maintenance, of pluripotency. Murine embryonic stem cells lacking Utx can execute lineage commitment and contribute to adult chimaeric animals; however, somatic cells lacking Utx fail to robustly reprogram back to the ground state of pluripotency. Utx directly partners with OSK reprogramming factors and uses its histone demethylase catalytic activity to facilitate iPSC formation. Genomic analysis indicates that Utx depletion results in aberrant dynamics of H3K27me3 repressive chromatin demethylation in somatic cells undergoing reprogramming. The latter directly hampers the derepression of potent pluripotency promoting gene modules (including Sall1, Sall4 and Utf1), which can cooperatively substitute for exogenous OSK supplementation in iPSC formation. Remarkably, Utx safeguards the timely execution of H3K27me3 demethylation observed in embryonic day 10.5-11 primordial germ cells (PGCs), and Utx-deficient PGCs show cell-autonomous aberrant epigenetic reprogramming dynamics during their embryonic maturation in vivo. Subsequently, this disrupts PGC development by embryonic day 12.5, and leads to diminished germline transmission in mouse chimaeras generated from Utx-knockout pluripotent cells. Thus, we identify Utx as a novel mediator with distinct functions during the re-establishment of pluripotency and germ cell development. Furthermore, our findings highlight the principle that molecular regulators mediating loss of repressive chromatin during in vivo germ cell reprogramming can be co-opted during in vitro reprogramming towards ground state pluripotency.

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.

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.

A critical role for ceramide synthase 2 in liver homeostasis: II. insights into molecular changes leading to hepatopathy.

We have generated a mouse that cannot synthesize very long acyl chain (C22-C24) ceramides (Pewzner-Jung, Y., Park, H., Laviad, E. L., Silva, L. C., Lahiri, S., Stiban, J., Erez-Roman, R., Brugger, B., Sachsenheimer, T., Wieland, F. T., Prieto, M., Merrill, A. H., and Futerman, A. H. (2010) J. Biol. Chem. 285, 10902-10910) due to ablation of ceramide synthase 2 (CerS2). As a result, significant changes were observed in the sphingolipid profile of livers from these mice, including elevated C16-ceramide and sphinganine levels. We now examine the functional consequences of these changes. CerS2 null mice develop severe nonzonal hepatopathy from about 30 days of age, the age at which CerS2 expression peaks in wild type mice, and display increased rates of hepatocyte apoptosis and proliferation. In older mice there is extensive and pronounced hepatocellular anisocytosis with widespread formation of nodules of regenerative hepatocellular hyperplasia. Progressive hepatomegaly and noninvasive hepatocellular carcinoma are also seen from approximately 10 months of age. Even though CerS2 is found at equally high mRNA levels in kidney and liver, there are no changes in renal function and no pathological changes in the kidney. High throughput analysis of RNA expression in liver revealed up-regulation of genes associated with cell cycle regulation, protein transport, cell-cell interactions and apoptosis, and down-regulation of genes associated with intermediary metabolism, such as lipid and steroid metabolism, adipocyte signaling, and amino acid metabolism. In addition, levels of the cell cycle regulator, the cyclin dependent-kinase inhibitor p21(WAF1/CIP1), were highly elevated, which occurs by at least two mechanisms, one of which may involve p53. We propose a functional rationale for the synthesis of sphingolipids with very long acyl chains in liver homeostasis and in cell physiology.

Fine mapping of AHI1 as a schizophrenia susceptibility gene: from association to evolutionary evidence.

In previous studies, we identified a locus for schizophrenia on 6q23.3 and proposed the Abelson helper integration site 1 (AHI1) as the candidate gene. AHI1 is expressed in the brain and plays a key role in neurodevelopment, is involved in Joubert syndrome, and has been recently associated with autism. The neurodevelopmental role of AHI1 fits with etiological hypotheses of schizophrenia. To definitively confirm our hypothesis, we searched for associations using a dense map of the region. Our strongest findings lay within the AHI1 gene: single-nucleotide polymorphisms rs11154801 and rs7759971 showed significant associations (P=6.23E-06; P=0.84E-06) and haplotypes gave P values in the 10E-8 to 10E-10 range. The second highest significant region maps close to AHI1 and includes the intergenic region between BC040979 and PDE7B (rs2038549 at P=9.70E-06 and rs1475069 at P=6.97E-06), and PDE7B and MAP7. Using a sample of Palestinian Arab families to confirm these findings, we found isolated signals. While these results did not retain their significance after correction for multiple testing, the joint analysis across the 2 samples supports the role of AHI1, despite the presence of heterogeneity. Given the hypothesis of positive selection of schizophrenia genes, we resequenced a 11 kb region within AHI1 in ethnically defined populations and found evidence for a selective sweep. Network analysis indicates 2 haplotype clades, with schizophrenia-susceptibility haplotypes clustering within the major clade. In conclusion, our data support the role of AHI1 as a susceptibility gene for schizophrenia and confirm it has been subjected to positive selection, also shedding light on new possible candidate genes, MAP7 and PDE7B.

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.

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.