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.

NLRP1 inflammasome activation induces pyroptosis of hematopoietic progenitor cells.

Cytopenias are key prognostic indicators of life-threatening infection, contributing to immunosuppression and mortality. Here we define a role for Caspase-1-dependent death, known as pyroptosis, in infection-induced cytopenias by studying inflammasome activation in hematopoietic progenitor cells. The NLRP1a inflammasome is expressed in hematopoietic progenitor cells and its activation triggers their pyroptotic death. Active NLRP1a induced a lethal systemic inflammatory disease that was driven by Caspase-1 and IL-1ß but was independent of apoptosis-associated speck-like protein containing a CARD (ASC) and ameliorated by IL-18. Surprisingly, in the absence of IL-1ß-driven inflammation, active NLRP1a triggered pyroptosis of hematopoietic progenitor cells resulting in leukopenia at steady state. During periods of hematopoietic stress induced by chemotherapy or lymphocytic choriomeningitis virus (LCMV) infection, active NLRP1a caused prolonged cytopenia, bone marrow hypoplasia, and immunosuppression. Conversely, NLRP1-deficient mice showed enhanced recovery from chemotherapy and LCMV infection, demonstrating that NLRP1 acts as a cellular sentinel to alert Caspase-1 to hematopoietic and infectious stress.

Haemopedia RNA-seq: a database of gene expression during haematopoiesis in mice and humans.

During haematopoiesis, haematopoietic stem cells differentiate into restricted potential progenitors before maturing into the many lineages required for oxygen transport, wound healing and immune response. We have updated Haemopedia, a database of gene-expression profiles from a broad spectrum of haematopoietic cells, to include RNA-seq gene-expression data from both mice and humans. The Haemopedia RNA-seq data set covers a wide range of lineages and progenitors, with 57 mouse blood cell types (flow sorted populations from healthy mice) and 12 human blood cell types. This data set has been made accessible for exploration and analysis, to researchers and clinicians with limited bioinformatics experience, on our online portal Haemosphere: https://www.haemosphere.org. Haemosphere also includes nine other publicly available high-quality data sets relevant to haematopoiesis. We have added the ability to compare gene expression across data sets and species by curating data sets with shared lineage designations or to view expression gene vs gene, with all plots available for download by the user.

Plasmacytoid dendritic cell heterogeneity is defined by CXCL10 expression following TLR7 stimulation.

Plasmacytoid dendritic cells (pDCs) play a critical role in bridging the innate and adaptive immune systems. pDCs are specialized type I interferon (IFN) producers, which has implicated them as initiators of autoimmune pathogenesis. However, little is known about the downstream effectors of type I IFN signaling that amplify autoimmune responses. Here, we have used a chemokine reporter mouse to determine the CXCR3 ligand responses in DCs subsets. Following TLR7 stimulation, conventional type 1 and type 2 DCs (cDC1 and cDC2, respectively) uniformly upregulate CXCL10. By contrast, the proportion of chemokine positive pDCs was significantly less, and stable CXCL10+ and CXCL10- populations could be distinguished. CXCL9 expression was induced in all cDC1s, in half of the cDC2 but not by pDCs. The requirement for IFNAR signaling for chemokine reporter expression was interrogated by receptor blocking and deficiency and shown to be critical for CXCR3 ligand expression in Flt3-ligand-derived DCs. Chemokine-producing potential was not concordant with the previously identified markers of pDC heterogeneity. Finally, we show that CXCL10+ and CXCL10- populations are transcriptionally distinct, expressing unique transcriptional regulators, IFN signaling molecules, chemokines, cytokines, and cell surface markers. This work highlights CXCL10 as a downstream effector of type I IFN signaling and suggests a division of labor in pDCs subtypes that likely impacts their function as effectors of viral responses and as drivers of inflammation.

Early lineage priming by trisomy of Erg leads to myeloproliferation in a Down syndrome model.

Down syndrome (DS), with trisomy of chromosome 21 (HSA21), is the commonest human aneuploidy. Pre-leukemic myeloproliferative changes in DS foetal livers precede the acquisition of GATA1 mutations, transient myeloproliferative disorder (DS-TMD) and acute megakaryocytic leukemia (DS-AMKL). Trisomy of the Erg gene is required for myeloproliferation in the Ts(1716)65Dn DS mouse model. We demonstrate here that genetic changes specifically attributable to trisomy of Erg lead to lineage priming of primitive and early multipotential progenitor cells in Ts(1716)65Dn mice, excess megakaryocyte-erythroid progenitors, and malignant myeloproliferation. Gene expression changes dependent on trisomy of Erg in Ts(1716)65Dn multilineage progenitor cells were correlated with those associated with trisomy of HSA21 in human DS hematopoietic stem and primitive progenitor cells. These data suggest a role for ERG as a regulator of hematopoietic lineage potential, and that trisomy of ERG in the context of DS foetal liver hemopoiesis drives the pre-leukemic changes that predispose to subsequent DS-TMD and DS-AMKL.

Mpl expression on megakaryocytes and platelets is dispensable for thrombopoiesis but essential to prevent myeloproliferation.

Thrombopoietin (TPO) acting via its receptor, the cellular homologue of the myeloproliferative leukemia virus oncogene (Mpl), is the major cytokine regulator of platelet number. To precisely define the role of specific hematopoietic cells in TPO-dependent hematopoiesis, we generated mice that express the Mpl receptor normally on stem/progenitor cells but lack expression on megakaryocytes and platelets (Mpl(PF4cre/PF4cre)). Mpl(PF4cre/PF4cre) mice displayed profound megakaryocytosis and thrombocytosis with a remarkable expansion of megakaryocyte-committed and multipotential progenitor cells, the latter displaying biological responses and a gene expression signature indicative of chronic TPO overstimulation as the underlying causative mechanism, despite a normal circulating TPO level. Thus, TPO signaling in megakaryocytes is dispensable for platelet production; its key role in control of platelet number is via generation and stimulation of the bipotential megakaryocyte precursors. Nevertheless, Mpl expression on megakaryocytes and platelets is essential to prevent megakaryocytosis and myeloproliferation by restricting the amount of TPO available to stimulate the production of megakaryocytes from the progenitor cell pool.

Concordant mast cell and basophil production by individual hematopoietic blast colony-forming cells.

Previous studies have shown that mouse bone marrow cells can produce mast cells when stimulated in vitro by stem cell factor (SCF) and interleukin-3 (IL-3). Experiments to define the marrow cells able to generate mast cells showed that the most active subpopulations were the Kit(+) Sca1(-) progenitor cell fraction and the more ancestral Kit(+) Sca1(+) blast colony-forming cell fraction. In clonal cultures, up to 64% of blast colony-forming cells were able to generate mast cells when stimulated by SCF and IL-3, and, of these, the most active were those in the CD34(-) Flt3R(-) long-term repopulating cell fraction. Basophils, identified by the monoclonal antibody mMCP-8 to mouse mast cell serine protease-8, were also produced by 50% of blast colony-forming cells with a strong concordance in the production of both cell types by individual blast colony-forming cells. Enriched populations of marrow-derived basophils were shown to generate variable numbers of mast cells after a further incubation with SCF and IL-3. The data extend the repertoire of lineage-committed cells able to be produced by multipotential hematopoietic blast colony-forming cells and show that basophils and mast cells can have common ancestral cells and that basophils can probably generate mast cells at least under defined in vitro conditions.

Plant phenological responses to a long-term experimental extension of growing season and soil warming in the tussock tundra of Alaska.

Climate warming is strongly altering the timing of season initiation and season length in the Arctic. Phenological activities are among the most sensitive plant responses to climate change and have important effects at all levels within the ecosystem. We tested the effects of two experimental treatments, extended growing season via snow removal and extended growing season combined with soil warming, on plant phenology in tussock tundra in Alaska from 1995 through 2003. We specifically monitored the responses of eight species, representing four growth forms: (i) graminoids (Carex bigellowii and Eriophorum vaginatum); (ii) evergreen shrubs (Ledum palustre, Cassiope tetragona, and Vaccinium vitis-idaea); (iii) deciduous shrubs (Betula nana and Salix pulchra); and (iv) forbs (Polygonum bistorta). Our study answered three questions: (i) Do experimental treatments affect the timing of leaf bud break, flowering, and leaf senescence? (ii) Are responses to treatments species-specific and growth form-specific? and (iii) Which environmental factors best predict timing of phenophases? Treatment significantly affected the timing of all three phenophases, although the two experimental treatments did not differ from each other. While phenological events began earlier in the experimental plots relative to the controls, duration of phenophases did not increase. The evergreen shrub, Cassiope tetragona, did not respond to either experimental treatment. While the other species did respond to experimental treatments, the total active period for these species did not increase relative to the control. Air temperature was consistently the best predictor of phenology. Our results imply that some evergreen shrubs (i.e., C. tetragona) will not capitalize on earlier favorable growing conditions, putting them at a competitive disadvantage relative to phenotypically plastic deciduous shrubs. Our findings also suggest that an early onset of the growing season as a result of decreased snow cover will not necessarily result in greater tundra productivity.

Reduced lymphocyte longevity and homeostatic proliferation in lamin B receptor-deficient mice results in profound and progressive lymphopenia.

The lamin B receptor (LBR) is a highly unusual inner nuclear membrane protein with multiple functions. Reduced levels are associated with decreased neutrophil lobularity, whereas complete absence of LBR results in severe skeletal dysplasia and in utero/perinatal lethality. We describe a mouse pedigree, Lym3, with normal bone marrow and thymic development but profound and progressive lymphopenia particularly within the T cell compartment. This defect arises from a point mutation within the Lbr gene with only trace mutant protein detectable in homozygotes, albeit sufficient for normal development. Reduced T cell homeostatic proliferative potential and life span in vivo were found to contribute to lymphopenia. To investigate the role of LBR in gene silencing in hematopoietic cells, we examined gene expression in wild-type and mutant lymph node CD8 T cells and bone marrow neutrophils. Although LBR deficiency had a very mild impact on gene expression overall, for common genes differentially expressed in both LBR-deficient CD8 T cells and neutrophils, gene upregulation prevailed, supporting a role for LBR in their suppression. In summary, this study demonstrates that LBR deficiency affects not only nuclear architecture but also proliferation, cell viability, and gene expression of hematopoietic cells.

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.

Regulation of hematopoietic stem cells by their mature progeny.

Thrombopoietin (TPO), acting through its receptor Mpl, has two major physiological roles: ensuring production of sufficient platelets via stimulation of megakaryocyte production and maintaining hematopoietic stem cell (HSC) quiescence. Mpl also controls circulating TPO concentration via receptor-mediated internalization and degradation. Here, we demonstrate that the megakaryocytosis and increased platelet mass in mice with mutations in the Myb or p300 genes causes reduced circulating TPO concentration and TPO starvation of the stem-cell compartment, which is exacerbated because these cells additionally exhibit impaired responsiveness to TPO. HSCs from Myb(Plt4/Plt4) mice show altered expression of TPO-responsive genes and, like HSCs from Tpo and Mpl mutant mice, exhibit increased cycling and a decline in the number of HSCs with age. These studies suggest that disorders of platelet number can have profound effects on the HSC compartment via effects on the feedback regulation of circulating TPO concentration.

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.

Estimating the proportion of microarray probes expressed in an RNA sample.

A fundamental question in microarray analysis is the estimation of the number of expressed probes in different RNA samples. Negative control probes available in the latest microarray platforms, such as Illumina whole genome expression BeadChips, provide a unique opportunity to estimate the number of expressed probes without setting a threshold. A novel algorithm was proposed in this study to estimate the number of expressed probes in an RNA sample by utilizing these negative controls to measure background noise. The performance of the algorithm was demonstrated by comparing different generations of Illumina BeadChips, comparing the set of probes targeting well-characterized RefSeq NM transcripts with other probes on the array and comparing pure samples with heterogenous samples. Furthermore, hematopoietic stem cells were found to have a larger transcriptome than progenitor cells. Aire knockout medullary thymic epithelial cells were shown to have significantly less expressed probes than matched wild-type cells.

Dual requirement for the ETS transcription factors Fli-1 and Erg in hematopoietic stem cells and the megakaryocyte lineage.

Fli-1 and Erg are closely related members of the Ets family of transcription factors. Both genes are translocated in human cancers, including Ewing's sarcoma, leukemia, and in the case of Erg, more than half of all prostate cancers. Although evidence from mice and humans suggests that Fli-1 is required for megakaryopoiesis, and that Erg is required for normal adult hematopoietic stem cell (HSC) regulation, their precise physiological roles remain to be defined. To elucidate the relationship between Fli-1 and Erg in hematopoiesis, we conducted an analysis of mice carrying mutations in both genes. Our results demonstrate that there is a profound genetic interaction between Fli-1 and Erg. Double heterozygotes displayed phenotypes more dramatic than single heterozygotes: severe thrombocytopenia, with a significant deficit in megakaryocyte numbers and evidence of megakaryocyte dysmorphogenesis, and loss of HSCs accompanied by a reduction in the number of committed hematopoietic progenitor cells. These results illustrate an indispensable requirement for both Fli-1 and Erg in normal HSC and megakaryocyte homeostasis, and suggest these transcription factors may coregulate common target genes.

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.

Fine mapping of Leishmania major susceptibility Locus lmr2 and evidence of a role for Fli1 in disease and wound healing.

Genetic linkage studies of the host response to Leishmania major, the causative agent of cutaneous leishmaniasis, have identified significant genetic complexity in humans and mice. In the mouse model, multiple loci have been implicated in susceptibility to infection, but to date, the genes underlying these loci have not been identified. We now describe the contribution of a novel candidate gene, Fli1, to both L. major resistance and enhanced wound healing. We have previously mapped the L. major response locus, lmr2, to proximal chromosome 9 in a genetic cross between the resistant C57BL/6 strain and the susceptible BALB/c strain. We now show that the presence of the resistant C57BL/6 lmr2 allele in susceptible BALB/c mice confers an enhanced L. major resistance and wound healing phenotype. Fine mapping of the lmr2 locus permitted the localization of the lmr2 quantitative trait locus to a 5-Mb interval comprising 21 genes, of which microarray analysis was able to identify differential expression in 1 gene-Fli1. Analysis of Fli1 expression in wounded and L. major-infected skin and naïve and infected lymph nodes validated the importance of Fli1 in lesion resolution and wound healing and identified 3 polymorphisms in the Fli1 promoter, among which a GA repeat element may be the important contributor.

Interconversion between Tumorigenic and Differentiated States in Acute Myeloid Leukemia.

Tumors are composed of phenotypically heterogeneous cancer cells that often resemble various differentiation states of their lineage of origin. Within this hierarchy, it is thought that an immature subpopulation of tumor-propagating cancer stem cells (CSCs) differentiates into non-tumorigenic progeny, providing a rationale for therapeutic strategies that specifically eradicate CSCs or induce their differentiation. The clinical success of these approaches depends on CSC differentiation being unidirectional rather than reversible, yet this question remains unresolved even in prototypically hierarchical malignancies, such as acute myeloid leukemia (AML). Here, we show in murine and human models of AML that, upon perturbation of endogenous expression of the lineage-determining transcription factor PU.1 or withdrawal of established differentiation therapies, some mature leukemia cells can de-differentiate and reacquire clonogenic and leukemogenic properties. Our results reveal plasticity of CSC maturation in AML, highlighting the need to therapeutically eradicate cancer cells across a range of differentiation states.

Transcriptional profiling of eosinophil subsets in interleukin-5 transgenic mice.

Eosinophils are important in fighting parasitic infections and are implicated in the pathogenesis of asthma and allergy. IL-5 is a critical regulator of eosinophil development, controlling proliferation, differentiation, and maturation of the lineage. Mice that constitutively express IL-5 have in excess of 10-fold more eosinophils in the hematopoietic organs than their wild type (WT) counterparts. We have identified that much of this expansion is in a population of Siglec-F high eosinophils, which are rare in WT mice. In this study, we assessed transcription in myeloid progenitors, eosinophil precursors, and Siglec-F medium and Siglec-F high eosinophils from IL-5 transgenic mice and in doing so have created a useful resource for eosinophil biologists. We have then utilized these populations to construct an eosinophil trajectory based on gene expression and to identify gene sets that are associated with eosinophil lineage progression. Cell cycle genes were significantly associated with the trajectory, and we experimentally demonstrate an increasing trend toward quiescence along the trajectory. Additionally, we found gene expression changes associated with constitutive IL-5 signaling in eosinophil progenitors, many of which were not observed in eosinophils.