A comprehensive single cell transcriptional landscape of human hematopoietic progenitors

Dynamics of human hematopoietic stem and progenitor cell differentiation to the erythroid lineage

Erythropoiesis, the development of erythrocytes from hematopoietic stem cells, occurs through four phases: erythroid progenitor (EP) development, early erythropoiesis, terminal erythroid differentiation (TED), and maturation. According to the classical model that is based on immunophenotypic profiles of cell populations, each of these phases comprises multiple differentiation states that arise in a hierarchical manner. After segregation of lymphoid potential, erythroid priming begins during progenitor development and progresses through progenitor cell types that have multilineage potential. Complete separation of the erythroid lineage is achieved during early erythropoiesis with the formation of unipotent EPs: burst-forming unit-erythroid and colony-forming unit-erythroid. These erythroid-committed progenitors undergo TED and maturation, which involves expulsion of the nucleus and remodeling to form functional biconcave, hemoglobin-filled erythrocytes. In the last decade or so, many studies employing advanced techniques such as single-cell RNA-sequencing (scRNA-seq) as well as the conventional methods, including colony-forming cell assays and immunophenotyping, have revealed heterogeneity within the stem, progenitor, and erythroblast stages, and uncovered alternate paths for segregation of erythroid lineage potential. In this review, we provide an in-depth account of immunophenotypic profiles of all cell types within erythropoiesis, highlight studies that demonstrate heterogeneous erythroid stages, and describe deviations to the classical model of erythropoiesis. Overall, although scRNA-seq approaches have provided new insights, flow cytometry remains relevant and is the primary method for validation of novel immunophenotypes.

Proteomic and transcriptomic screening demonstrates increased mast cell-derived CCL23 in systemic mastocytosis

BACKGROUND

Systemic mastocytosis (SM) is a heterogeneous group of mast cell-driven diseases diagnosed by bone marrow sampling. However, there are a limited number of available blood disease biomarkers.

OBJECTIVE

Our aim was to identify mast cell-derived proteins that could potentially serve as blood biomarkers for indolent and advanced forms of SM.

METHODS

We performed a plasma proteomics screening coupled with single-cell transcriptomic analysis in SM patients and healthy subjects.

RESULTS

Plasma proteomics screening identified 19 proteins upregulated in indolent disease compared to healthy, and 16 proteins in advanced disease compared to indolent. Among these, 5 proteins, CCL19, CCL23, CXCL13, IL-10, and IL-12Rβ1, were higher in indolent relative to healthy and in advanced disease compared to indolent. Single-cell RNA sequencing demonstrated that CCL23, IL-10, and IL-6 were selectively produced by mast cells. Notably, plasma CCL23 levels correlated positively with known markers of SM disease severity, namely tryptase levels, percentage bone marrow mast cell infiltration, and IL-6.

CONCLUSION

CCL23 is produced predominantly by mast cells in SM, and CCL23 plasma levels are associated with disease severity, correlating positively with established markers of disease burden, thus suggesting that CCL23 is a specific SM biomarker. In addition, the combination of CCL19, CCL23, CXCL13, IL-10, and IL-12Rβ1 may be useful for defining disease stage.

Post-Transplant Administration of G-CSF Impedes Engraftment of Gene Edited Human Hematopoietic Stem Cells by Exacerbating the p53-Mediated DNA Damage Response

Granulocyte colony stimulating factor (G-CSF) is commonly used as adjunct treatment to hasten recovery from neutropenia following chemotherapy and autologous transplantation of hematopoietic stem and progenitor cells (HSPCs) for malignant disorders. However, the utility of G-CSF administration after ex vivo gene therapy procedures targeting human HSPCs has not been thoroughly evaluated. Here, we provide evidence that post-transplant administration of G-CSF impedes engraftment of CRISPR-Cas9 gene edited human HSPCs in xenograft models. G-CSF acts by exacerbating the p53-mediated DNA damage response triggered by Cas9- mediated DNA double-stranded breaks. Transient p53 inhibition in culture attenuates the negative impact of G-CSF on gene edited HSPC function. In contrast, post-transplant administration of G-CSF does not impair the repopulating properties of unmanipulated human HSPCs or HSPCs genetically engineered by transduction with lentiviral vectors. The potential for post-transplant G-CSF administration to aggravate HSPC toxicity associated with CRISPR-Cas9 gene editing should be considered in the design of ex vivo autologous HSPC gene editing clinical trials.

Dissection and Integration of Bursty Transcriptional Dynamics for Complex Systems

RNA velocity estimation is a potentially powerful tool to reveal the directionality of transcriptional changes in single-cell RNA-seq data, but it lacks accuracy, absent advanced metabolic labeling techniques. We developed a novel approach, TopicVelo, that disentangles simultaneous, yet distinct, dynamics by using a probabilistic topic model, a highly interpretable form of latent space factorization, to infer cells and genes associated with individual processes, thereby capturing cellular pluripotency or multifaceted functionality. Focusing on process-associated cells and genes enables accurate estimation of process-specific velocities via a master equation for a transcriptional burst model accounting for intrinsic stochasticity. The method obtains a global transition matrix by leveraging cell topic weights to integrate process-specific signals. In challenging systems, this method accurately recovers complex transitions and terminal states, while our novel use of first-passage time analysis provides insights into transient transitions. These results expand the limits of RNA velocity, empowering future studies of cell fate and functional responses.

Always stressed but never exhausted: how stem cells in myeloid neoplasms avoid extinction in inflammatory conditions

Chronic or recurrent episodes of acute inflammation cause attrition of normal hematopoietic stem cells (HSCs) that can lead to hematopoietic failure but they drive progression in myeloid malignancies and their precursor clonal hematopoiesis. Mechanistic parallels exist between hematopoiesis in chronic inflammation and the continuously increased proliferation of myeloid malignancies, particularly myeloproliferative neoplasms (MPNs). The ability to enter dormancy, a state of deep quiescence characterized by low oxidative phosphorylation, low glycolysis, reduced protein synthesis, and increased autophagy is central to the preservation of long-term HSCs and likely MPN SCs. The metabolic features of dormancy resemble those of diapause, a state of arrested embryonic development triggered by adverse environmental conditions. To outcompete their normal counterparts in the inflammatory MPN environment, MPN SCs co-opt mechanisms used by HSCs to avoid exhaustion, including signal attenuation by negative regulators, insulation from activating cytokine signals, anti-inflammatory signaling, and epigenetic reprogramming. We propose that new therapeutic strategies may be derived from conceptualizing myeloid malignancies as an ecosystem out of balance, in which residual normal and malignant hematopoietic cells interact in multiple ways, only few of which have been characterized in detail. Disrupting MPN SC insulation to overcome dormancy, interfering with aberrant cytokine circuits that favor MPN cells, and directly boosting residual normal HSCs are potential strategies to tip the balance in favor of normal hematopoiesis. Although eradicating the malignant cell clones remains the goal of therapy, rebalancing the ecosystem may be a more attainable objective in the short term.

Hematopoietic stem cells and betaherpesvirus latency

The human betaherpesviruses including human cytomegalovirus (HCMV), human herpesvirus (HHV)-6a and HHV-6b, and HHV-7 infect and establish latency in CD34+ hematopoietic stem and progenitor cells (HPCs). The diverse repertoire of HPCs in humans and the complex interactions between these viruses and host HPCs regulate the viral lifecycle, including latency. Precise manipulation of host and viral factors contribute to preferential maintenance of the viral genome, increased host cell survival, and specific manipulation of the cellular environment including suppression of neighboring cells and immune control. The dynamic control of these processes by the virus regulate inter- and intra-host signals critical to the establishment of chronic infection. Regulation occurs through direct viral protein interactions and cellular signaling, miRNA regulation, and viral mimics of cellular receptors and ligands, all leading to control of cell proliferation, survival, and differentiation. Hematopoietic stem cells have unique biological properties and the tandem control of virus and host make this a unique environment for chronic herpesvirus infection in the bone marrow. This review highlights the elegant complexities of the betaherpesvirus latency and HPC virus-host interactions.

Systems biology approaches to unravel lymphocyte subsets and function

Single-cell technologies have revealed the extensive heterogeneity and complexity of the immune system. Systems biology approaches in immunology have taken advantage of the high-parameter, high-throughput data and analyzed immune cell types in a 'bottom-up' data-driven method. This approach has discovered previously unrecognized cell types and functions. Especially for human immunology, in which experimental manipulations are challenging, systems approach has become a successful means to investigate physiologically relevant contexts. This review focuses on the recent findings in lymphocyte biology, from their development, differentiation into subsets, and heterogeneity in their functions, enabled by these systems approaches. Furthermore, we review examples of the application of findings from systems approach studies and discuss how now to leave the rich dataset in the curse of high dimensionality.

Innovations in hematopoietic stem-cell mobilization: a review of the novel CXCR4 inhibitor motixafortide

Hematopoietic stem-cell transplantation (HCT) and stem-cell-based gene therapies rely on the ability to collect sufficient CD34+ hematopoietic stem and progenitor cells (HSPCs), typically via peripheral blood mobilization. Commonly used HSPC mobilization regimens include single-agent granulocyte colony-stimulating factor (G-CSF), plerixafor, chemotherapy, or a combination of these agents. These regimens, however, frequently require multiple days of injections and leukapheresis procedures to collect adequate HSPCs for HCT (minimum = >2 × 106 CD34+ cells/kg; optimal = 5-6 × 106 CD34+ cells/kg). In addition, these regimens frequently yield suboptimal CD34+ HSPC numbers for HSPC-based gene-edited therapies, given the significantly higher HSPC number needed for successful gene-editing and manufacturing. Meanwhile, G-CSF is associated with common adverse events such as bone pain as well as an increased risk of rare but potentially life-threatening splenic rupture. Moreover, G-CSF is unsafe in patients with sickle-cell disease, a key patient population that may benefit from autologous HSPC-based gene-edited therapies, where it has been associated with unacceptable rates of serious vaso-occlusive and thrombotic events. Motixafortide is a novel CXCR4 inhibitor with extended in vivo activity (>48 h) that has been shown in preclinical and clinical trials to rapidly mobilize robust numbers of HSPCs in preparation for HCT, while preferentially mobilizing increased numbers of more primitive HSPCs by immunophenotyping and single-cell RNA expression profiling. In this review, we present a history of stem-cell mobilization and update of recent innovations in novel mobilization strategies with a specific focus on the development of motixafortide, a long-acting CXCR4 inhibitor, as a novel HSPC mobilizing agent.

Single cell transcriptomics clarifies the basophil differentiation trajectory and identifies pre-basophils upstream of mature basophils

Basophils are the rarest granulocytes and are recognized as critical cells for type 2 immune responses. However, their differentiation pathway remains to be fully elucidated. Here, we assess the ontogenetic trajectory of basophils by single-cell RNA sequence analysis. Combined with flow cytometric and functional analyses, we identify c-Kit^-CLEC12A^hi pre-basophils located downstream of pre-basophil and mast cell progenitors (pre-BMPs) and upstream of CLEC12A^lo mature basophils. The transcriptomic analysis predicts that the pre-basophil population includes previously-defined basophil progenitor (BaP)-like cells in terms of gene expression profile. Pre-basophils are highly proliferative and respond better to non-IgE stimuli but less to antigen plus IgE stimulation than do mature basophils. Although pre-basophils usually remain in the bone marrow, they emerge in helminth-infected tissues, probably through IL-3-mediated inhibition of their retention in the bone marrow. Thus, the present study identifies pre-basophils that bridge the gap between pre-BMPs and mature basophils during basophil ontogeny.

Genetic engineering meets hematopoietic stem cell biology for next-generation gene therapy

The growing clinical success of hematopoietic stem/progenitor cell (HSPC) gene therapy (GT) relies on the development of viral vectors as portable "Trojan horses" for safe and efficient gene transfer. The recent advent of novel technologies enabling site-specific gene editing is broadening the scope and means of GT, paving the way to more precise genetic engineering and expanding the spectrum of diseases amenable to HSPC-GT. Here, we provide an overview of state-of-the-art and prospective developments of the HSPC-GT field, highlighting how advances in biological characterization and manipulation of HSPCs will enable the design of the next generation of these transforming therapeutics.

Basophils beyond allergic and parasitic diseases

Basophils bind IgE via FcεRI-αβγ2, which they uniquely share only with mast cells. In doing so, they can rapidly release mediators that are hallmark of allergic disease. This fundamental similarity, along with some morphological features shared by the two cell types, has long brought into question the biological significance that basophils mediate beyond that of mast cells. Unlike mast cells, which mature and reside in tissues, basophils are released into circulation from the bone marrow (constituting 1% of leukocytes), only to infiltrate tissues under specific inflammatory conditions. Evidence is emerging that basophils mediate non-redundant roles in allergic disease and, unsuspectingly, are implicated in a variety of other pathologies [e.g., myocardial infarction, autoimmunity, chronic obstructive pulmonary disease, fibrosis, cancer, etc.]. Recent findings strengthen the notion that these cells mediate protection from parasitic infections, whereas related studies implicate basophils promoting wound healing. Central to these functions is the substantial evidence that human and mouse basophils are increasingly implicated as important sources of IL-4 and IL-13. Nonetheless, much remains unclear regarding the role of basophils in pathology vs. homeostasis. In this review, we discuss the dichotomous (protective and/or harmful) roles of basophils in a wide spectrum of non-allergic disorders.

DRAG in situ barcoding reveals an increased number of HSPCs contributing to myelopoiesis with age

Ageing is associated with changes in the cellular composition of the immune system. During ageing, hematopoietic stem and progenitor cells (HSPCs) that produce immune cells are thought to decline in their regenerative capacity. However, HSPC function has been mostly assessed using transplantation assays, and it remains unclear how HSPCs age in the native bone marrow niche. To address this issue, we present an in situ single cell lineage tracing technology to quantify the clonal composition and cell production of single cells in their native niche. Our results demonstrate that a pool of HSPCs with unequal output maintains myelopoiesis through overlapping waves of cell production throughout adult life. During ageing, the increased frequency of myeloid cells is explained by greater numbers of HSPCs contributing to myelopoiesis rather than the increased myeloid output of individual HSPCs. Strikingly, the myeloid output of HSPCs remains constant over time despite accumulating significant transcriptomic changes throughout adulthood. Together, these results show that, unlike emergency myelopoiesis post-transplantation, aged HSPCs in their native microenvironment do not functionally decline in their regenerative capacity.

Histone H2A deubiquitinases in the transcriptional programs of development and hematopoiesis: a consolidated analysis

Monoubiquitinated lysine 119 of histone H2A (H2AK119ub) is a highly abundant epigenetic mark, associated with gene repression and deposited on chromatin by the polycomb repressor complex 1 (PRC1), which is an essential regulator of diverse transcriptional programs in mammalian development and tissue homeostasis. While multiple deubiquitinases (DUBs) with catalytic activity for H2AK119ub (H2A-DUBs) have been identified, we lack systematic analyses of their roles and cross-talk in transcriptional regulation. Here, we address H2A-DUB functions in epigenetic regulation of mammalian development and tissue maintenance by conducting a meta-analysis of 248 genomics datasets from 32 independent studies, focusing on the mouse model and covering embryonic stem cells (ESCs), hematopoietic, and immune cell lineages. This covers all the publicly available datasets that map genomic H2A-DUB binding and H2AK119ub distributions (ChIP-Seq), and all datasets assessing dysregulation in gene expression in the relevant H2A-DUB knockout models (RNA-Seq). Many accessory datasets for PRC1-2 and DUB-interacting proteins are also analyzed and interpreted, as well as further data assessing chromatin accessibility (ATAC-Seq) and transcriptional activity (RNA-seq). We report co-localization in the binding of H2A-DUBs BAP1, USP16, and to a lesser extent others that is conserved across different cell-types, and also the enrichment of antagonistic PRC1-2 protein complexes at the same genomic locations. Such conserved sites enriched for the H2A-DUBs and PRC1-2 are proximal to transcriptionally active genes that engage in housekeeping cellular functions. Nevertheless, they exhibit H2AK119ub levels significantly above the genomic average that can undergo further increase with H2A-DUB knockout. This indicates a cooperation between H2A-DUBs and PRC1-2 in the modulation of housekeeping transcriptional programs, conserved across many cell types, likely operating through their antagonistic effects on H2AK119ub and the regulation of local H2AK119ub turnover. Our study further highlights existing knowledge gaps and discusses important directions for future work.

Motixafortide and G-CSF to mobilize hematopoietic stem cells for autologous transplantation in multiple myeloma: a randomized phase 3 trial

Autologous hematopoietic stem cell transplantation (ASCT) improves survival in multiple myeloma (MM). However, many individuals are unable to collect optimal CD34+ hematopoietic stem and progenitor cell (HSPC) numbers with granulocyte colony-stimulating factor (G-CSF) mobilization. Motixafortide is a novel cyclic-peptide CXCR4 inhibitor with extended in vivo activity. The GENESIS trial was a prospective, phase 3, double-blind, placebo-controlled, multicenter study with the objective of assessing the superiority of motixafortide + G-CSF over placebo + G-CSF to mobilize HSPCs for ASCT in MM. The primary endpoint was the proportion of patients collecting ≥6 × 106 CD34+ cells kg-1 within two apheresis procedures; the secondary endpoint was to achieve this goal in one apheresis. A total of 122 adult patients with MM undergoing ASCT were enrolled at 18 sites across five countries and randomized (2:1) to motixafortide + G-CSF or placebo + G-CSF for HSPC mobilization. Motixafortide + G-CSF enabled 92.5% to successfully meet the primary endpoint versus 26.2% with placebo + G-CSF (odds ratio (OR) 53.3, 95% confidence interval (CI) 14.12-201.33, P < 0.0001). Motixafortide + G-CSF also enabled 88.8% to meet the secondary endpoint versus 9.5% with placebo + G-CSF (OR 118.0, 95% CI 25.36-549.35, P < 0.0001). Motixafortide + G-CSF was safe and well tolerated, with the most common treatment-emergent adverse events observed being transient, grade 1/2 injection site reactions (pain, 50%; erythema, 27.5%; pruritis, 21.3%). In conclusion, motixafortide + G-CSF mobilized significantly greater CD34+ HSPC numbers within two apheresis procedures versus placebo + G-CSF while preferentially mobilizing increased numbers of immunophenotypically and transcriptionally primitive HSPCs. Trial Registration: ClinicalTrials.gov , NCT03246529.

Hematopoietic Stem Cells and the Immune System in Development and Aging

Hematopoietic stem cells (HSCs) support haematopoiesis throughout life and give rise to the whole variety of cells of the immune system. Developing in the early embryo, passing through the precursor stage, and maturing into the first HSCs, they undergo a fairly large number of divisions while maintaining a high regenerative potential due to high repair activity. This potential is greatly reduced in adult HSCs. They go into a state of dormancy and anaerobic metabolism to maintain their stemness throughout life. However, with age, changes occur in the pool of HSCs that negatively affect haematopoiesis and the effectiveness of immunity. Niche aging and accumulation of mutations with age reduces the ability of HSCs to self-renew and changes their differentiation potential. This is accompanied by a decrease in clonal diversity and a disturbance of lymphopoiesis (decrease in the formation of naive T- and B-cells) and the predominance of myeloid haematopoiesis. Aging also affects mature cells, regardless of HSC, therefore, phagocytic activity and the intensity of the oxidative burst decrease, and the efficiency of processing and presentation of antigens by myeloid cells is impaired. Aging cells of innate and adaptive immunity produce factors that form a chronic inflammatory background. All these processes have a serious negative impact on the protective properties of the immune system, increasing inflammation, the risk of developing autoimmune, oncological, and cardiovascular diseases with age. Understanding the mechanisms of reducing the regenerative potential in a comparative analysis of embryonic and aging HSCs, the features of inflammatory aging will allow us to get closer to deciphering the programs for the development, aging, regeneration and rejuvenation of HSCs and the immune system.

Single-cell dissection of human hematopoietic reconstitution after allogeneic hematopoietic stem cell transplantation

Hematopoietic stem cell transplantation is an effective regenerative therapy for many malignant, inherited, or autoimmune diseases. However, our understanding of reconstituted hematopoiesis in transplant patients remains limited. Here, we uncover the reconstitution dynamics of human allogeneic hematopoietic stem and progenitor cells (HSPCs) at single-cell resolution after transplantation. Transplanted HSPCs underwent rapid and measurable changes during the first 30 days after transplantation, characterized by a strong proliferative response on the first day. Transcriptomic analysis of HSPCs enabled us to observe that immunoregulatory neutrophil progenitors expressing high levels of the S100A gene family were enriched in granulocyte colony-stimulating factor-mobilized peripheral blood stem cells. Transplant recipients who developed acute graft-versus-host disease (aGVHD) infused fewer S100A^high immunoregulatory neutrophil progenitors, immunophenotyped as Lin^-CD34⁺CD66b⁺CD177⁺, than those who did not develop aGVHD. Therefore, our study provides insights into the regenerative process of transplanted HSPCs in human patients and identifies a potential criterion for identifying patients at high risk for developing aGVHD early after transplant.

Low Expression of BRCA1 as a Potential Relapse Predictor in B-Cell Acute Lymphoblastic Leukemia

B-cell acute lymphoblastic leukemia (B-ALL) is the most common childhood hematological malignancy worldwide. Treatment outcomes have improved dramatically in recent years; despite this, relapse is still a problem, and the potential molecular explanation for this remains an important field of study. We performed microarray and single-cell RNA-Seq data mining, and we selected significant data with a P -value<0.05. We validated BRCA1 gene expression by means of quantitative (reverse transcription-polymerase chain reaction.) We performed statistical analysis and considered a P -value<0.05 significant. We identified the overexpression of breast cancer 1, early onset (BRCA1; P -value=2.52 -134 ), by means of microarray analysis. Moreover, the normal distribution of BRCA1 expression in healthy bone marrow. In addition, we confirmed the increases in BRCA1 expression using real-time (reverse transcription-polymerase chain reaction and determined that it was significantly reduced in patients with relapse ( P -values=0.026). Finally, we identified that the expression of the BRCA1 gene could predict early relapse ( P -values=0.01). We determined that low expression of BRCA1 was associated with B-cell acute lymphoblastic leukemia relapse and could be a potential molecular prognostic marker.

Multi-omic single-cell velocity models epigenome-transcriptome interactions and improves cell fate prediction

Multi-omic single-cell datasets, in which multiple molecular modalities are profiled within the same cell, offer an opportunity to understand the temporal relationship between epigenome and transcriptome. To realize this potential, we developed MultiVelo, a differential equation model of gene expression that extends the RNA velocity framework to incorporate epigenomic data. MultiVelo uses a probabilistic latent variable model to estimate the switch time and rate parameters of chromatin accessibility and gene expression and improves the accuracy of cell fate prediction compared to velocity estimates from RNA only. Application to multi-omic single-cell datasets from brain, skin and blood cells reveals two distinct classes of genes distinguished by whether chromatin closes before or after transcription ceases. We also find four types of cell states: two states in which epigenome and transcriptome are coupled and two distinct decoupled states. Finally, we identify time lags between transcription factor expression and binding site accessibility and between disease-associated SNP accessibility and expression of the linked genes. MultiVelo is available on PyPI, Bioconda and GitHub ( https://github.com/welch-lab/MultiVelo ).

C1Q labels a highly aggressive macrophage-like leukemia population indicating extramedullary infiltration and relapse

Extramedullary infiltration (EMI) is a concomitant manifestation that may indicate poor outcome of acute myeloid leukemia (AML). The underlying mechanism remains poorly understood and therapeutic options are limited. Here, we employed single-cell RNA sequencing on bone marrow (BM) and EMI samples from a patient with AML presenting pervasive leukemia cutis. A complement C1Q+ macrophage-like leukemia subset, which was enriched within cutis and existed in BM before EMI manifestations, was identified and further verified in multiple patients with AML. Genomic and transcriptional profiling disclosed mutation and gene expression signatures of patients with EMI that expressed high levels of C1Q. RNA sequencing and quantitative proteomic analysis revealed expression dynamics of C1Q from primary to relapse. Univariate and multivariate analysis demonstrated adverse prognosis significance of C1Q expression. Mechanistically, C1Q expression, which was modulated by transcription factor MAF BZIP transcription factor B, endowed leukemia cells with tissue infiltration ability, which could establish prominent cutaneous or gastrointestinal EMI nodules in patient-derived xenograft and cell line-derived xenograft models. Fibroblasts attracted migration of the C1Q+ leukemia cells through C1Q-globular C1Q receptor recognition and subsequent stimulation of transforming growth factor β1. This cell-to-cell communication also contributed to survival of C1Q+ leukemia cells under chemotherapy stress. Thus, C1Q served as a marker for AML with adverse prognosis, orchestrating cancer infiltration pathways through communicating with fibroblasts and represents a compelling therapeutic target for EMI.

A quantitative evaluation of topological motifs and their coupling in gene circuit state distributions

One of the major challenges in biology is to understand how gene interactions collaborate to determine overall functions of biological systems. Here, we present a new computational framework that enables systematic, high-throughput, and quantitative evaluation of how small transcriptional regulatory circuit motifs, and their coupling, contribute to functions of a dynamical biological system. We illustrate how this approach can be applied to identify four-node gene circuits, circuit motifs, and motif coupling responsible for various gene expression state distributions, including those derived from single-cell RNA sequencing data. We also identify seven major classes of four-node circuits from clustering analysis of state distributions. The method is applied to establish phenomenological models of gene circuits driving human neuron differentiation, revealing important biologically relevant regulatory interactions. Our study will shed light on a better understanding of gene regulatory mechanisms in creating and maintaining cellular states.

Deep learning models will shape the future of stem cell research

Our ability to understand and control stem cell biology is being augmented by developments on two fronts, our ability to collect more data describing cell state and our capability to comprehend these data using deep learning models. Here we consider the impact deep learning will have in the future of stem cell research. We explore the importance of generating data suitable for these methods, the requirement for close collaboration between experimental and computational researchers, and the challenges we face to do this fairly and effectively. Achieving this will ensure that the resulting deep learning models are biologically meaningful and computationally tractable.

Integrated single-cell transcriptome analysis of CD34 + enriched leukemic stem cells revealed intra- and inter-patient transcriptional heterogeneity in pediatric acute myeloid leukemia

To gain insights into the idiosyncrasies of CD34 + enriched leukemic stem cells, we investigated the nature and extent of transcriptional heterogeneity by single-cell sequencing in pediatric AML. Whole transcriptome analysis of 28,029 AML single cells was performed using the nanowell cartridge-based barcoding technology. Integrated transcriptional analysis identified unique leukemic stem cell clusters of each patient and intra-patient heterogeneity was revealed by multiple LSC-enriched clusters differing in their cell cycle processes and BCL2 expression. All LSC-enriched clusters exhibited gene expression profile of dormancy and self-renewal. Upregulation of genes involved in non-coding RNA processing and ribonucleoprotein assembly were observed in LSC-enriched clusters relative to HSC. The genes involved in regulation of apoptotic processes, response to cytokine stimulus, and negative regulation of transcription were upregulated in LSC-enriched clusters as compared to the blasts. Validation of top altered genes in LSC-enriched clusters confirmed upregulation of TCF7L2, JUP, ARHGAP25, LPAR6, and PRDX1 genes, and serine/threonine kinases (STK24, STK26). Upregulation of LPAR6 showed trend towards MRD positive status (Odds ratio = 0.126; 95% CI = 0.0144-1.10; p = 0.067) and increased expression of STK26 significantly correlated with higher RFS (HR = 0.231; 95% CI = 0.0506-1.052; p = 0.04). Our findings addressed the inter- and intra-patient diversity within AML LSC and potential signaling and chemoresistance-associated targets that warrant investigation in larger cohort that may guide precision medicine in the near future.

Stemness-related genes revealed by single-cell profiling of naïve and stimulated human CD34+ cells from CB and mPB

BACKGROUND

Hematopoietic stem cells (HSCs) from different sources show varied repopulating capacity, and HSCs lose their stemness after long-time ex vivo culture. A deep understanding of these phenomena may provide helpful insights for HSCs.

METHODS

Here, we applied single-cell RNA-seq (scRNA-seq) to analyse the naïve and stimulated human CD34+ cells from cord blood (CB) and mobilised peripheral blood (mPB).

RESULTS

We collected over 16 000 high-quality single-cell data to construct a comprehensive inference map and characterised the HSCs under a quiescent state on the hierarchy top. Then, we compared HSCs in CB with those in mPB and HSCs of naïve samples to those of cultured samples, and identified stemness-related genes (SRGs) associated with cell source (CS-SRGs) and culture time (CT-SRGs), respectively. Interestingly, CS-SRGs and CT-SRGs share genes enriched in the signalling pathways such as mRNA catabolic process, translational initiation, ribonucleoprotein complex biogenesis and cotranslational protein targeting to membrane, suggesting dynamic protein translation and processing may be a common requirement for stemness maintenance. Meanwhile, CT-SRGs are enriched in pathways involved in glucocorticoid and corticosteroid response that affect HSCs homing and engraftment. In contrast, CS-SRGs specifically contain genes related to purine and ATP metabolic process, which is crucial for HSC homeostasis in the stress settings. Particularly, when CT-SRGs are used as reference genes for the construction of the development trajectory of CD34+ cells, lymphoid and myeloid lineages are clearly separated after HSCs/MPPs. Finally, we presented an application through a small-scale drug screening using Connectivity Map (CMap) against CT-SRGs. A small molecule, cucurbitacin I, was found to efficiently expand HSCs ex vivo while maintaining its stemness.

CONCLUSIONS

Our findings provide new perspectives for understanding HSCs, and the strategy to identify candidate molecules through SRGs may be applicable to study other stem cells.

Transplanting rejuvenated blood stem cells extends lifespan of aged immunocompromised mice

One goal of regenerative medicine is to rejuvenate tissues and extend lifespan by restoring the function of endogenous aged stem cells. However, evidence that somatic stem cells can be targeted in vivo to extend lifespan is still lacking. Here, we demonstrate that after a short systemic treatment with a specific inhibitor of the small RhoGTPase Cdc42 (CASIN), transplanting aged hematopoietic stem cells (HSCs) from treated mice is sufficient to extend the healthspan and lifespan of aged immunocompromised mice without additional treatment. In detail, we show that systemic CASIN treatment improves strength and endurance of aged mice by increasing the myogenic regenerative potential of aged skeletal muscle stem cells. Further, we show that CASIN modifies niche localization and H4K16ac polarity of HSCs in vivo. Single-cell profiling reveals changes in HSC transcriptome, which underlie enhanced lymphoid and regenerative capacity in serial transplantation assays. Overall, we provide proof-of-concept evidence that a short systemic treatment to decrease Cdc42 activity improves the regenerative capacity of different endogenous aged stem cells in vivo, and that rejuvenated HSCs exert a broad systemic effect sufficient to extend murine health- and lifespan.

Single cell multi-omic reference atlases of non-human primate immune tissues reveals CD102 as a biomarker for long-lived plasma cells

In response to infection or immunization, antibodies are produced that provide protection against re-exposure with the same pathogen. These antibodies can persist at high titers for decades and are maintained by bone marrow-resident long-lived plasma cells (LLPC). However, the durability of antibody responses to immunization varies amongst vaccines. It is unknown what factors contribute to the differential longevity of serum antibody responses and whether heterogeneity in LLPC contributes to this phenomenon. While LLPC differentiation has been studied extensively in mice, little is known about this population in humans or non-human primates (NHP). Here, we use multi-omic single-cell profiling to identify and characterize the LLPC compartment in NHP. We identify LLPC biomarkers including the marker CD102 and show that CD102 in combination with CD31 identifies LLPC in NHP bone marrow. Additionally, we find that CD102 is expressed by LLPC in mouse and humans. These results further our understanding of the LLPC compartment in NHP, identify biomarkers of LLPC, and provide tissue-specific single cell references for future studies.

Maternal Western-style diet remodels the transcriptional landscape of fetal hematopoietic stem and progenitor cells in rhesus macaques

Maternal obesity adversely impacts the in utero metabolic environment, but its effect on fetal hematopoiesis remains incompletely understood. During late development, the fetal bone marrow (FBM) becomes the major site where macrophages and B lymphocytes are produced via differentiation of hematopoietic stem and progenitor cells (HSPCs). Here, we analyzed the transcriptional landscape of FBM HSPCs at single-cell resolution in fetal macaques exposed to a maternal high-fat Western-style diet (WSD) or a low-fat control diet. We demonstrate that maternal WSD induces a proinflammatory response in FBM HSPCs and fetal macrophages. In addition, maternal WSD consumption suppresses the expression of B cell development genes and decreases the frequency of FBM B cells. Finally, maternal WSD leads to poor engraftment of fetal HSPCs in nonlethally irradiated immunodeficient NOD/SCID/IL2rγ^-/- mice. Collectively, these data demonstrate for the first time that maternal WSD impairs fetal HSPC differentiation and function in a translationally relevant nonhuman primate model.

scCapsNet-mask: an updated version of scCapsNet with extended applicability in functional analysis related to scRNA-seq data

BACKGROUND

With the rapid accumulation of scRNA-seq data, more and more automatic cell type identification methods have been developed, especially those based on deep learning. Although these methods have reached relatively high prediction accuracy, many issues still exist. One is the interpretability. The second is how to deal with the non-standard test samples that are not encountered in the training process.

RESULTS

Here we introduce scCapsNet-mask, an updated version of scCapsNet. The scCapsNet-mask provides a reasonable solution to the issues of interpretability and non-standard test samples. Firstly, the scCapsNet-mask utilizes a mask to ease the task of model interpretation in the original scCapsNet. The results show that scCapsNet-mask could constrain the coupling coefficients, and make a one-to-one correspondence between the primary capsules and type capsules. Secondly, the scCapsNet-mask can process non-standard samples more reasonably. In one example, the scCapsNet-mask was trained on the committed cells, and then tested on less differentiated cells as the non-standard samples. It could not only estimate the lineage bias of less differentiated cells, but also distinguish the development stages more accurately than traditional machine learning models. Therefore, the pseudo-temporal order of cells for each lineage could be established. Following these pseudo-temporal order, lineage specific genes exhibit a gradual increase expression pattern and stem cell associated genes exhibit a gradual decrease expression pattern. In another example, the scCapsNet-mask was trained on scRNA-seq data, and then used to assign cell type in spatial transcriptomics that may contain non-standard sample of doublets. The results show that the scCapsNet-mask not only restored the spatial map but also identified several non-standard samples of doublet.

CONCLUSIONS

The scCapsNet-mask offers a suitable solution to the challenge of interpretability and non-standard test samples. By adding a mask, it has the advantages of automatic processing and easy interpretation compared with the original scCapsNet. In addition, the scCapsNet-mask could more accurately reflect the composition of non-standard test samples than traditional machine learning methods. Therefore, it can extend its applicability in functional analysis, such as fate bias prediction in less differentiated cells and cell type assignment in spatial transcriptomics.

Current knowledge on the early stages of human neutropoiesis

Polymorphonuclear neutrophils are no longer considered as a homogeneous population of terminally differentiated and short-lived cells that belong to the innate immune system only. In fact, data from the past decades have uncovered that neutrophils exhibit large phenotypic heterogeneity and functional versatility that render them more plastic than previously thought. Hence, their precise role as effector cells in inflammation, in immune response and in other pathophysiological processes, including tumors, needs to be better evaluated. In such a complex scenario, common knowledge of the differentiation of neutrophils in bone marrow refers to lineage precursors, starting from the still poorly defined myeloblasts, and proceeding sequentially to promyelocytes, myelocytes, metamyelocytes, band cells, segmented neutrophils, and mature neutrophils, with each progenitor stage being more mature and better characterized. Thanks to the development and utilization of cutting-edge technologies, novel information about neutrophil precursors at stages earlier than the promyelocytes, hence closer to the hematopoietic stem cells, is emerging. Accordingly, this review discusses the main findings related to the very early precursors of human neutrophils and provides our perspectives on human neutropoiesis.

E-Cadherin Expression Distinguishes Mouse from Human Hematopoiesis in the Basophil and Erythroid Lineages

E-cadherin is a key regulator of epithelial cell-cell adhesion, the loss of which accelerates tumor growth and invasion. E-cadherin is also expressed in hematopoietic cells as well as epithelia. The function of hematopoietic E-cadherin is, however, mostly elusive. In this study, we explored the validity of mouse models to functionally investigate the role of hematopoietic E-cadherin in human hematopoiesis. We generated a hematopoietic-specific E-cadherin knockout mouse model. In mice, hematopoietic E-cadherin is predominantly expressed within the basophil lineage, the expression of which is dispensable for the generation of basophils. However, neither E-cadherin mRNA nor protein were detected in human basophils. In contrast, human hematopoietic E-cadherin marks the erythroid lineage. E-cadherin expression in hematopoiesis thereby revealed striking evolutionary differences between the basophil and erythroid cell lineage in humans and mice. This is remarkable as E-cadherin expression in epithelia is highly conserved among vertebrates including humans and mice. Our study therefore revealed that the mouse does not represent a suitable model to study the function of E-cadherin in human hematopoiesis and an alternative means to study the role of E-cadherin in human erythropoiesis needs to be developed.

Raman spectroscopy differ leukemic cells from their healthy counterparts and screen biomarkers in acute leukemia

Precision medicine is important in the treatment of acute leukemia (AL). The target therapies of AL provide an opportunity to reduce the mortality of AL. How AL cells differ from their healthy counterparts is the basis for the development of therapies and the outcome of AL patients. Therefore, a label-free and noninvasive single-cell Raman platform was used to characterize cell molecular profiles and found potential biomarkers from three healthy people and twelve AL patients with more than 90% accuracy. We analyzed myeloblasts, abnormal promyelocytes, monoblasts and B-ALL cells respectively, compared with their healthy counterparts, which could be distinguished by their intrinsic phenotypic Raman spectra using orthogonal partial least squares discriminate analysis (OPLS-DA). Most importantly, we selected statistically significant markers of the four leukemia models. Further analysis of leukemic granulocytes, we found that a combination of the 1003, 1341 and 1579 cm^-1 Raman peaks could discriminate myeloblasts and abnormal promyelocytes from normal granulocytes. The assignments of 1579 cm^-1 gave us a clue to find potential important variables myeloperoxidase related with AL diagnosis. Our study demonstrates the capability of the Raman platform to characterize leukemia cells with non-invasively probing metabolites. The biomarker we identified could be extensible to other blood cells and potentially have a high impact on leukemia therapy.

Multi-omics profiling identifies C1QA/B+ macrophages with multiple immune checkpoints associated with esophageal squamous cell carcinoma (ESCC) liver metastasis

Background

Esophageal squamous cell carcinoma (ESCC) is a highly lethal malignant tumor lacking effective treatments; 20% of ESCC patients develop liver metastasis with an extremely short survival time of ≈5 months. The tumor microenvironment (TME) plays a crucial role in tumor homeostasis, but the relationship between the ESCC TME and liver metastasis is still unknown.

Methods

To identify potential cell populations contributing to ESCC liver metastasis, single-cell RNA (scRNA) sequencing data were analyzed to identify the major cell populations within the TME. Each of the major cell populations was re-clustered to define detailed cell subsets. Thereafter, the gene set variation analysis (GSVA) score was calculated for the bulk RNA-seq data based on the gene signatures of each cell subset. The relationship between the GSVA score of each cellular subset and clinical outcome was further analyzed to identify the cellular subset associated with ESCC liver metastasis, which was validated by multiplex immunohistochemistry.

Results

C1QA/B⁺ tumor-associated macrophages (TAMs) acted as the central regulator of the ESCC TME, closely associated with several key cell subsets. Several immune checkpoints, including CD40, CD47 and LGALS9, were all positively expressed in C1QA/B⁺ macrophages, which may exert central regulatory control of immune evasion by ESCC via these immune checkpoints expressions.

Conclusions

Our results comprehensively revealed the landscape of tumor-infiltrating immune cells associated with ESCC prognosis and metastasis, and suggest a novel strategy for developing immunotherapies for ESCC liver metastasis by targeting C1QA/B⁺ TAMs.

CPHEN-013: Comprehensive phenotyping of hematopoietic stem and progenitor cells in the human fetal liver

Hematopoietic stem cells (HSCs) reside at the top of the hematopoietic hierarchy and can give rise to all the mature blood cell types in our body, while at the same time maintaining a pool of HSCs through self-renewing divisions. This potential is reflected in their functional definition as cells that are capable of long-term multi-lineage engraftment upon transplantation. While all HSCs meet these criteria, subtle differences exist between developmentally different populations of these cells. Here we present a comprehensive overview of traditional and more recently described markers for phenotyping HSCs and their downstream progeny. To address the need to assess the growing number of surface molecules expressed in various HSC-enriched fractions at different developmental stages, we have developed an extensive multi-parameter spectral flow cytometry panel to phenotype hematopoietic stem and multipotent progenitor cells (HSC/MPPs) throughout development. In this study we then employ this panel to comprehensively profile the HSC compartment in the human fetal liver (FL), which is endowed with superior engraftment potential compared to postnatal sources. Spectral cytometry lends an improved resolution of marker expression to our comprehensive approach, allowing to extract combinatorial expression signatures of several relevant HSC/MPP markers to precisely characterize the HSC/MPP fraction in a variety of tissues.

Role of TET dioxygenases in the regulation of both normal and pathological hematopoiesis

The family of ten-eleven translocation dioxygenases (TETs) consists of TET1, TET2, and TET3. Although all TETs are expressed in hematopoietic tissues, only TET2 is commonly found to be mutated in age-related clonal hematopoiesis and hematopoietic malignancies. TET2 mutation causes abnormal epigenetic landscape changes and results in multiple stages of lineage commitment/differentiation defects as well as genetic instability in hematopoietic stem/progenitor cells (HSPCs). TET2 mutations are founder mutations (first hits) in approximately 40–50% of cases of TET2-mutant (TET2^MT) hematopoietic malignancies and are later hits in the remaining cases. In both situations, TET2^MT collaborates with co-occurring mutations to promote malignant transformation. In TET2^MT tumor cells, TET1 and TET3 partially compensate for TET2 activity and contribute to the pathogenesis of TET2^MT hematopoietic malignancies. Here we summarize the most recent research on TETs in regulating of both normal and pathogenic hematopoiesis. We review the concomitant mutations and aberrant signals in TET2^MT malignancies. We also discuss the molecular mechanisms by which concomitant mutations and aberrant signals determine lineage commitment in HSPCs and the identity of hematopoietic malignancies. Finally, we discuss potential strategies to treat TET2^MT hematopoietic malignancies, including reverting the methylation state of TET2 target genes and targeting the concomitant mutations and aberrant signals.

Single-cell multi-omics of human clonal hematopoiesis reveals that DNMT3A R882 mutations perturb early progenitor states through selective hypomethylation

Somatic mutations in cancer genes have been detected in clonal expansions across healthy human tissue, including in clonal hematopoiesis. However, because mutated and wild-type cells are admixed, we have limited ability to link genotypes with phenotypes. To overcome this limitation, we leveraged multi-modality single-cell sequencing, capturing genotype, transcriptomes and methylomes in progenitors from individuals with DNMT3A R882 mutated clonal hematopoiesis. DNMT3A mutations result in myeloid over lymphoid bias, and an expansion of immature myeloid progenitors primed toward megakaryocytic-erythroid fate, with dysregulated expression of lineage and leukemia stem cell markers. Mutated DNMT3A leads to preferential hypomethylation of polycomb repressive complex 2 targets and a specific CpG flanking motif. Notably, the hypomethylation motif is enriched in binding motifs of key hematopoietic transcription factors, serving as a potential mechanistic link between DNMT3A mutations and aberrant transcriptional phenotypes. Thus, single-cell multi-omics paves the road to defining the downstream consequences of mutations that drive clonal mosaicism.

Murine fetal bone marrow does not support functional hematopoietic stem and progenitor cells until birth

While adult bone marrow (BM) hematopoietic stem and progenitor cells (HSPCs) and their extrinsic regulation is well studied, little is known about the composition, function, and extrinsic regulation of the first HSPCs to enter the BM during development. Here, we functionally interrogate murine BM HSPCs from E15.5 through P0. Our work reveals that fetal BM HSPCs are present by E15.5, but distinct from the HSPC pool seen in fetal liver, both phenotypically and functionally, until near birth. We also generate a transcriptional atlas of perinatal BM HSPCs and the BM niche in mice across ontogeny, revealing that fetal BM lacks HSPCs with robust intrinsic stem cell programs, as well as niche cells supportive of HSPCs. In contrast, stem cell programs are preserved in neonatal BM HSPCs, which reside in a niche expressing HSC supportive factors distinct from those seen in adults. Collectively, our results provide important insights into the factors shaping hematopoiesis during this understudied window of hematopoietic development.

Relatively little is known about the first hematopoietic stem and progenitor cells to arrive in the fetal bone marrow. Here they characterize the frequency, function, and molecular identity of fetal BM HSPCs and their bone marrow niche, and show that most BM HSPCs have little hematopoietic function until birth.

Flow Cytometric Identification of Hematopoietic and Leukemic Blast Cells for Tailored Clinical Follow-Up of Acute Myeloid Leukemia

Acute myeloid leukemia (AML) is a myeloid malignancy that is characterized by the accumulation of leukemic blast cells, which originate from hematopoietic stem cells that have undergone leukemic transformation and/or are more mature progenitors that have gained stemness features. Currently, no consensus exists for the flow cytometric identification of normal blast cells and their leukemic counterparts by their antigenic expression profile. Differentiating between the benign cells and the malignant cells is crucial for the further deployment of immunophenotype panels for the clinical follow-up of AML patients. This review provides an overview of immunophenotypic markers that allow the identification of leukemic blast cells in the bone marrow with multiparameter flow cytometry. This technique allows the identification of hematopoietic blast cells at the level of maturing cells by their antigen expression profile. While aberrant antigen expression of a single immunophenotypic marker cell cannot be utilized in order to differentiate leukemic blast cells from normal blast cells, combinations of multiple immunophenotypic markers can enable the distinction of normal and leukemic blast cells. The identification of these markers has provided new perspectives for tailored clinical follow-up, including therapy management, diagnostics, and prognostic purposes. The immunophenotypic marker panels, however, should be developed by carefully considering the variable antigen marker expression profile of individual patients.

PLAG1 dampens protein synthesis to promote human hematopoietic stem cell self-renewal

Hematopoietic stem cell (HSC) dormancy is understood as supportive of HSC function and its long-term integrity. Although regulation of stress responses incurred as a result of HSC activation is recognized as important in maintaining stem cell function, little is understood of the preventive machinery present in human HSCs that may serve to resist their activation and promote HSC self-renewal. We demonstrate that the transcription factor PLAG1 is essential for long-term HSC function and, when overexpressed, endows a 15.6-fold enhancement in the frequency of functional HSCs in stimulatory conditions. Genome-wide measures of chromatin occupancy and PLAG1-directed gene expression changes combined with functional measures reveal that PLAG1 dampens protein synthesis, restrains cell growth and division, and enhances survival, with the primitive cell advantages it imparts being attenuated by addition of the potent translation activator, c-MYC. We find PLAG1 capitalizes on multiple regulatory factors to ensure protective diminished protein synthesis including 4EBP1 and translation-targeting miR-127 and does so independently of stress response signaling. Overall, our study identifies PLAG1 as an enforcer of human HSC dormancy and self-renewal through its highly context-specific regulation of protein biosynthesis and classifies PLAG1 among a rare set of bona fide regulators of messenger RNA translation in these cells. Our findings showcase the importance of regulated translation control underlying human HSC physiology, its dysregulation under activating demands, and the potential if its targeting for therapeutic benefit.

Clonal Hematopoiesis: Role in Hematologic and Non-Hematologic Malignancies

Hematopoietic stem cells (HSCs) ensure the coordinated and balanced production of all hematopoietic cell types throughout life. Aging is associated with a gradual decline of the self-renewal and regenerative potential of HSCs and with the development of clonal hematopoiesis. Clonal hematopoiesis of indeterminate potential (CHIP) defines the clonal expansion of genetically variant hematopoietic cells bearing one or more gene mutations and/or structural variants (such as copy number alterations). CHIP increases exponentially with age and is associated with cancers, including hematologic neoplasia, cardiovascular and other diseases. The presence of CHIP consistently increases the risk of hematologic malignancy, particularly in individuals who have CHIP in association with peripheral blood cytopenia.

A bifurcation concept for B-lymphoid/plasmacytoid dendritic cells with largely fluctuating transcriptome dynamics

Hematopoiesis was considered a hierarchical stepwise process but was revised to a continuous process following single-cell RNA sequencing. However, the uncertainty or fluctuation of single-cell transcriptome dynamics during differentiation was not considered, and the dendritic cell (DC) pathway in the lymphoid context remains unclear. Here, we identify human B-plasmacytoid DC (pDC) bifurcation as large fluctuating transcriptome dynamics in the putative B/NK progenitor region by dry and wet methods. By converting splicing kinetics into diffusion dynamics in a deep generative model, our original computational methodology reveals strong fluctuation at B/pDC bifurcation in IL-7Rα⁺ regions, and LFA-1 fluctuates positively in the pDC direction at the bifurcation. These expectancies are validated by the presence of B/pDC progenitors in the IL-7Rα⁺ fraction and preferential expression of LFA-1 in pDC-biased progenitors with a niche-like culture system. We provide a model of fluctuation-based differentiation, which reconciles continuous and discrete models and is applicable to other developmental systems.

Single-cell transcriptomics reveals the identity and regulators of human mast cell progenitors

Mast cell accumulation is a hallmark of a number of diseases, including allergic asthma and systemic mastocytosis. Immunoglobulin E-mediated crosslinking of the FcεRI receptors causes mast cell activation and contributes to disease pathogenesis. The mast cell lineage is one of the least studied among the hematopoietic cell lineages, and controversies remain about whether FcεRI expression appears during the mast cell progenitor stage or during terminal mast cell maturation. Here, we used single-cell transcriptomics analysis to reveal a temporal association between the appearance of FcεRI and the mast cell gene signature in CD34+ hematopoietic progenitors in adult peripheral blood. In agreement with these data, the FcεRI+ hematopoietic progenitors formed morphologically, phenotypically, and functionally mature mast cells in long-term culture assays. Single-cell transcriptomics analysis further revealed the expression patterns of prospective cytokine receptors regulating development of mast cell progenitors. Culture assays showed that interleukin-3 (IL-3) and IL-5 promoted disparate effects on progenitor cell proliferation and survival, respectively, whereas IL-33 caused robust FcεRI downregulation. Taken together, we showed that FcεRI expression appears at the progenitor stage of mast cell differentiation in peripheral blood. We also showed that external stimuli regulate FcεRI expression of mast cell progenitors, providing a possible explanation for the variable FcεRI expression levels during mast cell development.

Genome-wide analyses of 200,453 individuals yield new insights into the causes and consequences of clonal hematopoiesis

Clonal hematopoiesis (CH), the clonal expansion of a blood stem cell and its progeny driven by somatic driver mutations, affects over a third of people, yet remains poorly understood. Here we analyze genetic data from 200,453 UK Biobank participants to map the landscape of inherited predisposition to CH, increasing the number of germline associations with CH in European-ancestry populations from 4 to 14. Genes at new loci implicate DNA damage repair (PARP1, ATM, CHEK2), hematopoietic stem cell migration/homing (CD164) and myeloid oncogenesis (SETBP1). Several associations were CH-subtype-specific including variants at TCL1A and CD164 that had opposite associations with DNMT3A- versus TET2-mutant CH, the two most common CH subtypes, proposing key roles for these two loci in CH development. Mendelian randomization analyses showed that smoking and longer leukocyte telomere length are causal risk factors for CH and that genetic predisposition to CH increases risks of myeloproliferative neoplasia, nonhematological malignancies, atrial fibrillation and blood epigenetic ageing.

ACUTE AND SUSTAINED ALTERATIONS TO THE BONE MARROW IMMUNE MICROENVIRONMENT FOLLOWING POLYMICROBIAL INFECTION

ABSTRACT

Sepsis is a highly prevalent cause of death in intensive care units. Characterized by severe immune cell derangements, sepsis is often associated with multiorgan dysfunction. For many sepsis survivors, these deficits can persist long after clinical resolution of the underlying infection. Although many studies report on the impact of sepsis on individual immune cell subtypes, a comprehensive analysis of sepsis-induced alterations within and across the immune cell landscape is lacking. In this study, we used single-cell RNA sequencing to assess sepsis-associated transcriptional changes in immune cells isolated from bone marrow at single-cell resolution. We used a high-survival fecal-induced peritonitis sepsis model using Friend leukemia virus B mice. Single-cell RNA sequencing classified 3402 single cells from control subjects into 14 clusters representing long-term hematopoietic stem cell (HSC), short-term HSC, basophil, dendritic cell, eosinophil, erythroblast, erythrocyte, macrophage, neutrophil, natural killer cell, plasma cell, plasmacytoid dendritic cell, pre-B cell, and T memory cell lineages. One day following experimentally induced sepsis, cell type compositions shifted significantly and included notable decreases in HSC and myeloid cell abundance. In addition to proportional cell composition changes, acute sepsis induced significant transcriptional alterations in most immune cell types analyzed-changes that failed to completely resolve 1 month after sepsis. Taken together, we report widespread and persistent transcriptional changes in diverse immune cells in response to polymicrobial infection. This study will serve as a valuable resource for future work investigating acute and/or long-term sepsis-associated immune cell derangements.

ORIGINS: a protein network-based approach to quantify cell pluripotency from scRNA-seq data

Trajectory inference is a common application of scRNA-seq data. However, it is often necessary to previously determine the origin of the trajectories, the stem or progenitor cells. In this work, we propose a computational tool to quantify pluripotency from single cell transcriptomics data. This approach uses the protein-protein interaction (PPI) network associated with the differentiation process as a scaffold and the gene expression matrix to calculate a score that we call differentiation activity. This score reflects how active the differentiation network is in each cell. We benchmark the performance of our algorithm with two previously published tools, LandSCENT (Chen et al., 2019) and CytoTRACE (Gulati et al., 2020), for four healthy human data sets: breast, colon, hematopoietic and lung. We show that our algorithm is more efficient than LandSCENT and requires less RAM memory than the other programs. We also illustrate a complete workflow from the count matrix to trajectory inference using the breast data set.•

ORIGINS is a methodology to quantify pluripotency from scRNA-seq data implemented as a freely available R package.

•

ORIGINS uses the protein-protein interaction network associated with differentiation and the data set expression matrix to calculate a score (differentiation activity) that quantifies pluripotency for each cell.

scDART: integrating unmatched scRNA-seq and scATAC-seq data and learning cross-modality relationship simultaneously

It is a challenging task to integrate scRNA-seq and scATAC-seq data obtained from different batches. Existing methods tend to use a pre-defined gene activity matrix to convert the scATAC-seq data into scRNA-seq data. The pre-defined gene activity matrix is often of low quality and does not reflect the dataset-specific relationship between the two data modalities. We propose scDART, a deep learning framework that integrates scRNA-seq and scATAC-seq data and learns cross-modalities relationships simultaneously. Specifically, the design of scDART allows it to preserve cell trajectories in continuous cell populations and can be applied to trajectory inference on integrated data.

Fbxw11 impairs the repopulation capacity of hematopoietic stem/progenitor cells

BACKGROUND

The ubiquitin-proteasome system plays important roles in maintaining the self-renewal and differentiation of stem and progenitor cells through highly ordered degradation of cellular proteins. Fbxw11, an E3 ligase, participates in many important biological processes by targeting a broad range of proteins. However, its roles in hematopoietic stem/progenitor cells (HSPCs) have not been established.

METHODS

In this study, the effects of Fbxw11 on HSPCs were studied in vitro and in vivo by an overexpression strategy. Real-time PCR was performed to detect the expression of Fbxw11 in hematopoietic subpopulations. Colony-forming assays were performed to evaluate the in vitro function of Fbxw11 on HSPCs. Hoechst 33342 and Ki67 staining was performed to determine the cell-cycle distribution of HSPCs. Competitive transplantation experiments were used to evaluate the effect of Fbxw11 on the reconstitution potential of HSPCs. Single-cell RNA sequencing (scRNA-seq) was employed to reveal the transcriptomic alterations in HSPCs.

RESULTS

The expression of Fbxw11 was higher in Lin^-c-Kit⁺Sca-1⁺ (LSK) cells and myeloid progenitors than in lymphoid progenitors. Fbxw11 played negative roles in colony-forming and quiescence maintenance of HSPCs in vitro. Furthermore, serial competitive transplantation experiments revealed that Fbxw11 impaired the repopulation capacity of HSPCs. The proportion of granulocytes (Gr-1⁺CD11b⁺) in the differentiated mature cells was significantly higher than that in the control group, T cells and B cells were lower. Moreover, scRNA-seq revealed seven cell clusters in HSPCs. In addition, Fbxw11 downregulated the expression of Cebpa, Myc and Arid5b, which are significant regulators of HSPC activity, in most cell clusters.

CONCLUSION

Our data demonstrate that Fbxw11 plays a negative role in the maintenance of HSPCs in vitro and repopulation capacity in vivo. Our data also provide valuable transcriptome references for HSPCs in homeostasis.

Unique molecular and functional features of extramedullary hematopoietic stem and progenitor cell reservoirs in humans

Rare hematopoietic stem and progenitor cell (HSPC) pools outside the bone marrow (BM) contribute to blood production in stress and disease but remain ill-defined. Although nonmobilized peripheral blood (PB) is routinely sampled for clinical management, the diagnosis and monitoring potential of PB HSPCs remain untapped, as no healthy PB HSPC baseline has been reported. Here we comprehensively delineate human extramedullary HSPC compartments comparing spleen, PB, and mobilized PB to BM using single-cell RNA-sequencing and/or functional assays. We uncovered HSPC features shared by extramedullary tissues and others unique to PB. First, in contrast to actively dividing BM HSPCs, we found no evidence of substantial ongoing hematopoiesis in extramedullary tissues at steady state but report increased splenic HSPC proliferative output during stress erythropoiesis. Second, extramedullary hematopoietic stem cells/multipotent progenitors (HSCs/MPPs) from spleen, PB, and mobilized PB share a common transcriptional signature and increased abundance of lineage-primed subsets compared with BM. Third, healthy PB HSPCs display a unique bias toward erythroid-megakaryocytic differentiation. At the HSC/MPP level, this is functionally imparted by a subset of phenotypic CD71+ HSCs/MPPs, exclusively producing erythrocytes and megakaryocytes, highly abundant in PB but rare in other adult tissues. Finally, the unique erythroid-megakaryocytic-skewing of PB is perturbed with age in essential thrombocythemia and β-thalassemia. Collectively, we identify extramedullary lineage-primed HSPC reservoirs that are nonproliferative in situ and report involvement of splenic HSPCs during demand-adapted hematopoiesis. Our data also establish aberrant composition and function of circulating HSPCs as potential clinical indicators of BM dysfunction.

Frequent HLA-DR loss on hematopoietic stem progenitor cells in patients with cyclosporine-dependent aplastic anemia carrying HLA-DR15

To determine whether antigen presentation by HLA-DR on hematopoietic stem progenitor cells (HSPCs) is involved in the development of acquired aplastic anemia (AA), we studied the HLA-DR expression on CD45dimCD34+CD38+ cells in the peripheral blood of 61 AA patients including 23 patients possessing HLA-class I allele-lacking (HLA-class I[-]) leukocytes. HLA-DR-lacking (DR[-]) cells accounted for 13.0-57.1% of the total HSPCs in seven (11.5%) patients with HLA-DR15 who did not possess HLA-class I(-) leukocytes. The incubation of sorted DR(-) HSPCs in the presence of IFN-γ for 72 h resulted in the full restoration of the DR expression. A comparison of the transcriptome profile between DR(-) and DR(+) HSPCs revealed the lower expression of immune response-related genes including co-stimulatory molecules (e.g., CD48, CD74, and CD86) in DR(-) cells, which was not evident in HLA-class I(-) HSPCs. DR(-) cells were exclusively detected in GPI(+) HSPCs in four patients whose HSPCs could be analyzed separately for GPI(+) and GPI(-) HSPCs. These findings suggest that CD4+ T cells specific to antigens presented by HLA-DR15 on HSPCs may contribute to the development of AA as well as the immune escape of GPI(-) HSPCs in a distinct way from CD8+ T cells recognizing HLA-class I-restricted antigens.

Erythroid Cell Research: 3D Chromatin, Transcription Factors and Beyond

Studies of the regulatory networks and signals controlling erythropoiesis have brought important insights in several research fields of biology and have been a rich source of discoveries with far-reaching implications beyond erythroid cells biology. The aim of this review is to highlight key recent discoveries and show how studies of erythroid cells bring forward novel concepts and refine current models related to genome and 3D chromatin organization, signaling and disease, with broad interest in life sciences.

Inference of cell state transitions and cell fate plasticity from single-cell with MARGARET

Despite recent advances in inferring cellular dynamics using single-cell RNA-seq data, existing trajectory inference (TI) methods face difficulty in accurately reconstructing the cell-state manifold and cell-fate plasticity for complex topologies. Here, we present MARGARET (https://github.com/Zafar-Lab/Margaret) for inferring single-cell trajectory and fate mapping for diverse dynamic cellular processes. MARGARET reconstructs complex trajectory topologies using a deep unsupervised metric learning and a graph-partitioning approach based on a novel connectivity measure, automatically detects terminal cell states, and generalizes the quantification of fate plasticity for complex topologies. On a diverse benchmark consisting of synthetic and real datasets, MARGARET outperformed state-of-the-art methods in recovering global topology and cell pseudotime ordering. For human hematopoiesis, MARGARET accurately identified all major lineages and associated gene expression trends and helped identify transitional progenitors associated with key branching events. For embryoid body differentiation, MARGARET identified novel transitional populations that were validated by bulk sequencing and functionally characterized different precursor populations in the mesoderm lineage. For colon differentiation, MARGARET characterized the lineage for BEST4/OTOP2 cells and the heterogeneity in goblet cell lineage in the colon under normal and inflamed ulcerative colitis conditions. Finally, we demonstrated that MARGARET can scale to large scRNA-seq datasets consisting of ∼ millions of cells.

hECA: The cell-centric assembly of a cell atlas

The accumulation of massive single-cell omics data provides growing resources for building biomolecular atlases of all cells of human organs or the whole body. The true assembly of a cell atlas should be cell-centric rather than file-centric. We developed a unified informatics framework for seamless cell-centric data assembly and built the human Ensemble Cell Atlas (hECA) from scattered data. hECA v1.0 assembled 1,093,299 labeled human cells from 116 published datasets, covering 38 organs and 11 systems. We invented three new methods of atlas applications based on the cell-centric assembly: “in data” cell sorting for targeted data retrieval with customizable logic expressions, “quantitative portraiture” for multi-view representations of biological entities, and customizable reference creation for generating references for automatic annotations. Case studies on agile construction of user-defined sub-atlases and “in data” investigation of CAR-T off-targets in multiple organs showed the great potential enabled by the cell-centric ensemble atlas.

•

A unified informatics framework for seamless cell-centric assembly of massive single-cell data

•

Built the general-purpose human Ensemble Cell Atlas (hECA) V1.0 from scattered data

•

Three new methods of applications enabling “in data” cell experiments and portraiture

•

Case studies of agile atlas reconstruction and target therapies side-effect discovery