The Human Cell Atlas bone marrow single-cell interactive web portal

Integrated single cell analysis shows chronic alcohol drinking disrupts monocyte differentiation in the bone marrow

Chronic heavy alcohol drinking (CHD) rewires monocytes and macrophages toward heightened inflammatory states with compromised antimicrobial defenses that persist after 1-month abstinence. To determine whether these changes are mediated through alterations in the bone marrow niche, we profiled monocytes and hematopoietic stem cell progenitors (HSCPs) from CHD rhesus macaques using a combination of functional assays and single cell genomics. CHD resulted in transcriptional profiles consistent with increased activation and inflammation within bone marrow resident monocytes and macrophages. Furthermore, CHD resulted in transcriptional signatures associated with increased oxidative and cellular stress in HSCP. Differentiation of HSCP in vitro revealed skewing toward monocytes expressing "neutrophil-like" markers with greater inflammatory responses to bacterial agonists. Further analyses of HSCPs showed broad epigenetic changes that were in line with exacerbated inflammatory responses within monocytes and their progenitors. In summary, CHD alters HSCPs in the bone marrow leading to the production of monocytes poised to generate dysregulated hyper-inflammatory responses.

Single cell genomics as a transformative approach for aquaculture research and innovation

Single cell genomics encompasses a suite of rapidly maturing technologies that measure the molecular profiles of individual cells within target samples. These approaches provide a large up-step in biological information compared to long-established 'bulk' methods that profile the average molecular profiles of all cells in a sample, and have led to transformative advances in understanding of cellular biology, particularly in humans and model organisms. The application of single cell genomics is fast expanding to non-model taxa, including aquaculture species, where numerous research applications are underway with many more envisaged. In this review, we highlight the potential transformative applications of single cell genomics in aquaculture research, considering barriers and potential solutions to the broad uptake of these technologies. Focusing on single cell transcriptomics, we outline considerations for experimental design, including the essential requirement to obtain high quality cells/nuclei for sequencing in ectothermic aquatic species. We further outline data analysis and bioinformatics considerations, tailored to studies with the under-characterized genomes of aquaculture species, where our knowledge of cellular heterogeneity and cell marker genes is immature. Overall, this review offers a useful source of knowledge for researchers aiming to apply single cell genomics to address biological challenges faced by the global aquaculture sector though an improved understanding of cell biology.

Single-cell RNA sequencing of human non-hematopoietic bone marrow cells reveals a unique set of inter-species conserved biomarkers for native mesenchymal stromal cells

BACKGROUND

Native bone marrow (BM) mesenchymal stem/stromal cells (BM-MSCs) participate in generating and shaping the skeleton and BM throughout the lifespan. Moreover, BM-MSCs regulate hematopoiesis by contributing to the hematopoietic stem cell niche in providing critical cytokines, chemokines and extracellular matrix components. However, BM-MSCs contain a heterogeneous cell population that remains ill-defined. Although studies on the taxonomy of native BM-MSCs in mice have just started to emerge, the taxonomy of native human BM-MSCs remains unelucidated.

METHODS

By using single-cell RNA sequencing (scRNA-seq), we aimed to define a proper taxonomy for native human BM non-hematopoietic subsets including endothelial cells (ECs) and mural cells (MCs) but with a focal point on MSCs. To this end, transcriptomic scRNA-seq data were generated from 5 distinct BM donors and were analyzed together with other transcriptomic data and with computational biology analyses at different levels to identify, characterize and classify distinct native cell subsets with relevant biomarkers.

RESULTS

We could ascribe novel specific biomarkers to ECs, MCs and MSCs. Unlike ECs and MCs, MSCs exhibited an adipogenic transcriptomic pattern while co-expressing genes related to hematopoiesis support and multilineage commitment potential. Furthermore, by a comparative analysis of scRNA-seq of BM cells from humans and mice, we identified core genes conserved in both species. Notably, we identified MARCKS, CXCL12, PDGFRA, and LEPR together with adipogenic factors as archetypal biomarkers of native MSCs within BM. In addition, our data suggest some complex gene nodes regulating critical biological functions of native BM-MSCs together with a preferential commitment toward an adipocyte lineage.

CONCLUSIONS

Overall, our taxonomy for native BM non-hematopoietic compartment provides an explicit depiction of gene expression in human ECs, MCs and MSCs at single-cell resolution. This analysis helps enhance our understanding of the phenotype and the complexity of biological functions of native human BM-MSCs.

Early activation of inflammatory pathways in UBA1-mutated hematopoietic stem and progenitor cells in VEXAS

VEXAS (vacuoles, E1 enzyme, X-linked, autoinflammatory, somatic) syndrome is a pleiotropic, severe autoinflammatory disease caused by somatic mutations in the ubiquitin-like modifier activating enzyme 1 (UBA1) gene. To elucidate VEXAS pathophysiology, we performed transcriptome sequencing of single bone marrow mononuclear cells and hematopoietic stem and progenitor cells (HSPCs) from VEXAS patients. HSPCs are biased toward myeloid (granulocytic) differentiation, and against lymphoid differentiation in VEXAS. Activation of multiple inflammatory pathways (interferons and tumor necrosis factor alpha) occurs ontogenically early in primitive hematopoietic cells and particularly in the myeloid lineage in VEXAS, and inflammation is prominent in UBA1-mutated cells. Dysregulation in protein degradation likely leads to higher stress response in VEXAS HSPCs, which positively correlates with inflammation. TCR usage is restricted and there are increased cytotoxicity and IFN-γ signaling in T cells. In VEXAS syndrome, both aberrant inflammation and myeloid predominance appear intrinsic to hematopoietic stem cells mutated in UBA1.

Single-cell chromatin accessibility profiling of acute myeloid leukemia reveals heterogeneous lineage composition upon therapy-resistance

Acute myeloid leukemia (AML) is a heterogeneous disease characterized by high rate of therapy resistance. Since the cell of origin can impact response to therapy, it is crucial to understand the lineage composition of AML cells at time of therapy resistance. Here we leverage single-cell chromatin accessibility profiling of 22 AML bone marrow aspirates from eight patients at time of therapy resistance and following subsequent therapy to characterize their lineage landscape. Our findings reveal a complex lineage architecture of therapy-resistant AML cells that are primed for stem and progenitor lineages and spanning quiescent, activated and late stem cell/progenitor states. Remarkably, therapy-resistant AML cells are also composed of cells primed for differentiated myeloid, erythroid and even lymphoid lineages. The heterogeneous lineage composition persists following subsequent therapy, with early progenitor-driven features marking unfavorable prognosis in The Cancer Genome Atlas AML cohort. Pseudotime analysis further confirms the vast degree of heterogeneity driven by the dynamic changes in chromatin accessibility. Our findings suggest that therapy-resistant AML cells are characterized not only by stem and progenitor states, but also by a continuum of differentiated cellular lineages. The heterogeneity in lineages likely contributes to their therapy resistance by harboring different degrees of lineage-specific susceptibilities to therapy.

Profiling of miRNAs and their interfering targets in peripheral blood mononuclear cells from patients with chronic myeloid leukaemia

Introduction

MicroRNAs may be implicated in the acquisition of drug resistance in chronic myeloid leukemia as they regulate the expression of not only BCR-ABL1 but also genes associated with the activation of drug transfer proteins or essential signaling pathways.

Methods

To understand the impact of specifically expressed miRNAs in chronic myeloid leukemia and their target genes, we collected peripheral blood mononuclear cells (PBMC) from patients diagnosed with chronic myeloid leukemia (CML) and healthy donors to determine whole miRNA expression by small RNA sequencing and screened out 31 differentially expressed microRNAs (DE-miRNAs) with high expression. With the utilization of miRNA set enrichment analysis tools, we present here a comprehensive analysis of the relevance of DE-miRNAs to disease and biological function. Furthermore, the literature-based miRNA-target gene database was used to analyze the overall target genes of the DE-miRNAs and to define their associated biological responses. We further integrated DE-miRNA target genes to identify CML miRNA targeted gene signature singscore (CMTGSS) and used gene-set enrichment analysis (GSEA) to analyze the correlation between CMTGSS and Hallmark gene-sets in PBMC samples from clinical CML patients. Finally, the association of CMTGSS stratification with multiple CML cell lineage gene sets was validated in PBMC samples from CML patients using GSEA.

Results

Although individual miRNAs have been reported to have varying degrees of impact on CML, overall, our results show that abnormally upregulated miRNAs are associated with apoptosis and aberrantly downregulated miRNAs are associated with cell cycle. The clinical database shows that our defined DE-miRNAs are associated with the prognosis of CML patients. CMTGSS-based stratification analysis presented a tendency for miRNAs to affect cell differentiation in the blood microenvironment.

Conclusion

Collectively, this study defined differentially expressed miRNAs by miRNA sequencing from clinical samples and comprehensively analyzed the biological functions of the differential miRNAs in association with the target genes. The analysis of the enrichment of specific myeloid differentiated cells and immune cells also suggests the magnitude and potential targets of differentially expressed miRNAs in the clinical setting. It helps us to make links between the different results obtained from the multi-faceted studies to provide more potential research directions.

Transcriptomic classes of BCR-ABL1 lymphoblastic leukemia

In BCR-ABL1 lymphoblastic leukemia, treatment heterogeneity to tyrosine kinase inhibitors (TKIs), especially in the absence of kinase domain mutations in BCR-ABL1, is poorly understood. Through deep molecular profiling, we uncovered three transcriptomic subtypes of BCR-ABL1 lymphoblastic leukemia, each representing a maturation arrest at a stage of B-cell progenitor differentiation. An earlier arrest was associated with lineage promiscuity, treatment refractoriness and poor patient outcomes. A later arrest was associated with lineage fidelity, durable leukemia remissions and improved patient outcomes. Each maturation arrest was marked by specific genomic events that control different transition points in B-cell development. Interestingly, these events were absent in BCR-ABL1⁺ preleukemic stem cells isolated from patients regardless of subtype, which supports that transcriptomic phenotypes are determined downstream of the leukemia-initialing event. Overall, our data indicate that treatment response and TKI efficacy are unexpected outcomes of the differentiation stage at which this leukemia transforms.

Pan-cancer T cell atlas links a cellular stress response state to immunotherapy resistance

Tumor-infiltrating T cells offer a promising avenue for cancer treatment, yet their states remain to be fully characterized. Here we present a single-cell atlas of T cells from 308,048 transcriptomes across 16 cancer types, uncovering previously undescribed T cell states and heterogeneous subpopulations of follicular helper, regulatory and proliferative T cells. We identified a unique stress response state, TSTR, characterized by heat shock gene expression. TSTR cells are detectable in situ in the tumor microenvironment across various cancer types, mostly within lymphocyte aggregates or potential tertiary lymphoid structures in tumor beds or surrounding tumor edges. T cell states/compositions correlated with genomic, pathological and clinical features in 375 patients from 23 cohorts, including 171 patients who received immune checkpoint blockade therapy. We also found significantly upregulated heat shock gene expression in intratumoral CD4/CD8+ cells following immune checkpoint blockade treatment, particularly in nonresponsive tumors, suggesting a potential role of TSTR cells in immunotherapy resistance. Our well-annotated T cell reference maps, web portal and automatic alignment/annotation tool could provide valuable resources for T cell therapy optimization and biomarker discovery.

Neutrophils and emergency granulopoiesis drive immune suppression and an extreme response endotype during sepsis

Sepsis arises from diverse and incompletely understood dysregulated host response processes following infection that leads to life-threatening organ dysfunction. Here we showed that neutrophils and emergency granulopoiesis drove a maladaptive response during sepsis. We generated a whole-blood single-cell multiomic atlas (272,993 cells, n = 39 individuals) of the sepsis immune response that identified populations of immunosuppressive mature and immature neutrophils. In co-culture, CD66b⁺ sepsis neutrophils inhibited proliferation and activation of CD4⁺ T cells. Single-cell multiomic mapping of circulating hematopoietic stem and progenitor cells (HSPCs) (29,366 cells, n = 27) indicated altered granulopoiesis in patients with sepsis. These features were enriched in a patient subset with poor outcome and a specific sepsis response signature that displayed higher frequencies of IL1R2⁺ immature neutrophils, epigenetic and transcriptomic signatures of emergency granulopoiesis in HSPCs and STAT3-mediated gene regulation across different infectious etiologies and syndromes. Our findings offer potential therapeutic targets and opportunities for stratified medicine in severe infection.

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.

Erythroid/megakaryocytic differentiation confers BCL-XL dependency and venetoclax resistance in acute myeloid leukemia

Myeloid neoplasms with erythroid or megakaryocytic differentiation include pure erythroid leukemia, myelodysplastic syndrome with erythroid features, and acute megakaryoblastic leukemia (FAB M7) and are characterized by poor prognosis and limited treatment options. Here, we investigate the drug sensitivity landscape of these rare malignancies. We show that acute myeloid leukemia (AML) cells with erythroid or megakaryocytic differentiation depend on the antiapoptotic protein B-cell lymphoma (BCL)-XL, rather than BCL-2, using combined ex vivo drug sensitivity testing, genetic perturbation, and transcriptomic profiling. High-throughput screening of >500 compounds identified the BCL-XL-selective inhibitor A-1331852 and navitoclax as highly effective against erythroid/megakaryoblastic leukemia cell lines. In contrast, these AML subtypes were resistant to the BCL-2 inhibitor venetoclax, which is used clinically in the treatment of AML. Consistently, genome-scale CRISPR-Cas9 and RNAi screening data demonstrated the striking essentiality of BCL-XL-encoding BCL2L1 but not BCL2 or MCL1, for the survival of erythroid/megakaryoblastic leukemia cell lines. Single-cell and bulk transcriptomics of patient samples with erythroid and megakaryoblastic leukemias identified high BCL2L1 expression compared with other subtypes of AML and other hematological malignancies, where BCL2 and MCL1 were more prominent. BCL-XL inhibition effectively killed blasts in samples from patients with AML with erythroid or megakaryocytic differentiation ex vivo and reduced tumor burden in a mouse erythroleukemia xenograft model. Combining the BCL-XL inhibitor with the JAK inhibitor ruxolitinib showed synergistic and durable responses in cell lines. Our results suggest targeting BCL-XL as a potential therapy option in erythroid/megakaryoblastic leukemias and highlight an AML subgroup with potentially reduced sensitivity to venetoclax-based treatments.

Integrated single cell analysis shows chronic alcohol drinking disrupts monocyte differentiation in the bone marrow niche

Chronic alcohol drinking rewires circulating monocytes and tissue-resident macrophages towards heightened inflammatory states with compromised anti-microbial defenses. As these effects remain consistent in short-lived monocytes after a 1-month abstinence period it is unclear whether these changes are restricted to the periphery or mediated through alterations in the progenitor niche. To test this hypothesis, we profiled monocytes/macrophages and hematopoietic stem cell progenitors (HSCP) of the bone marrow compartment from rhesus macaques after 12 months of ethanol consumption using a combination of functional assays and single cell genomics. Bone marrow-resident monocytes/macrophages from ethanol-consuming animals exhibited heightened inflammation. Differentiation of HSCP in vitro revealed skewing towards monocytes expressing neutrophil-like markers with heightened inflammatory responses to bacterial agonists. Single cell transcriptional analysis of HSCPs showed reduced proliferation but increased inflammatory markers in mature myeloid progenitors. We observed transcriptional signatures associated with increased oxidative and cellular stress as well as oxidative phosphorylation in immature and mature myeloid progenitors. Single cell analysis of the chromatin landscape showed altered drivers of differentiation in monocytes and progenitors. Collectively, these data indicate that chronic ethanol drinking results in remodeling of the transcriptional and epigenetic landscapes of the bone marrow compartment leading to altered functions in the periphery.

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.

Establishing a human bone marrow single cell reference atlas to study ageing and diseases

Introduction

Ageing in the human bone marrow is associated with immune function decline that results in the elderly being vulnerable to illnesses. A comprehensive healthy bone marrow consensus atlas can serve as a reference to study the immunological changes associated with ageing, and to identify and study abnormal cell states.

Methods

We collected publicly available single cell transcriptomic data of 145 healthy samples encompassing a wide spectrum of ages ranging from 2 to 84 years old to construct our human bone marrow atlas. The final atlas has 673,750 cells and 54 annotated cell types.

Results

We first characterised the changes in cell population sizes with respect to age and the corresponding changes in gene expression and pathways. Overall, we found significant age-associated changes in the lymphoid lineage cells. The naïve CD8+ T cell population showed significant shrinkage with ageing while the effector/memory CD4+ T cells increased in proportion. We also found an age-correlated decline in the common lymphoid progenitor population, in line with the commonly observed myeloid skew in haematopoiesis among the elderly. We then employed our cell type-specific ageing gene signatures to develop a machine learning model that predicts the biological age of bone marrow samples, which we then applied to healthy individuals and those with blood diseases. Finally, we demonstrated how to identify abnormal cell states by mapping disease samples onto the atlas. We accurately identified abnormal plasma cells and erythroblasts in multiple myeloma samples, and abnormal cells in acute myeloid leukaemia samples.

Discussion

The bone marrow is the site of haematopoiesis, a highly important bodily process. We believe that our healthy bone marrow atlas is a valuable reference for studying bone marrow processes and bone marrow-related diseases. It can be mined for novel discoveries, as well as serve as a reference scaffold for mapping samples to identify and investigate abnormal cells.

Characterization of bone marrow heterogeneity in NK-AML (M4/M5) based on single-cell RNA sequencing

Normal karyotype acute myeloid leukemia (NK-AML) is a heterogeneous hematological malignancy that contains a minor population of self-renewing leukemia stem cells (LSCs), which complicate efforts to achieve long-term survival. We performed single-cell RNA sequencing to profile 39,288 cells from 6 bone marrow (BM) aspirates including 5 NK-AML (M4/M5) patients and 1 healthy donor. The single-cell transcriptome atlas and gene expression characteristics of each cell population in NK-AML (M4/M5) and healthy BM were obtained. In addition, we identified a distinct LSC-like cluster with possible biomarkers in NK-AML (M4/M5) and verified 6 genes using qRT‒PCR and bioinformatic analyses. In conclusion, we utilized single-cell technologies to provide an atlas of NK-AML (M4/M5) cell heterogeneity, composition, and biomarkers with implications for precision medicine and targeted therapies.

From helping to regulating – A transcriptomic profile of Ifng+ Il10+ Il21+ Cd4+ Th1 cells indicates their role in regulating inflammation during experimental trypanosomosis

Introduction

Trypanosoma evansi parasite infections cause a chronic animal wasting disease called Surra, and cases of atypical Human Trypanosomosis (aHT). In experimental models, T. evansi infections are hallmarked by the early onset of excessive inflammation. Therefore, balancing the production of inflammatory cytokines by anti-inflammatory IL-10 is crucial for prolonged survival.

Methods

To improve the understanding of trypanosomosis induced immunopathology, we used scRNA-seq data from an experimental chronic T. evansi infection mouse model, resembling natural infection in terms of disease characteristics.

Results and discussion

For the first time, obtained results allowed to assess the transcriptomic profile and heterogeneity of splenic CD4+ T cell subsets, during a trypanosome infection. Here, the predominant subpopulation of T cells was represented by Tbx21(T-bet)+Ccr5+ Id2+ type 1 helper T cells (Th1), followed by Icos+ Cxcr5+Follicular T helper cells (Tfh) and very minor fraction of Il2ra(CD25)+Foxp3+ regulatory T cells (Tregs). Interestingly, the profile of Th1 cells shows that besides Ifng, these cells express high levels of Il10 and Il21, coding for anti-inflammatory and immunoregulatory cytokines. This coincides with the elevated expression of key genes involved in IL-10 and IL-21 secretion pathway such as Stat1 and Stat3, as well as the transcriptional factors Prdm1 (Blimp 1), and Maf (c-Maf). In contrast, there is virtually no IL-10 transcription detected in the Treg population. Finally, differential gene expression and gene ontology analysis of infection-induced Ifng+ Il10+ Il21+ Th1 cells highlights their suppressive function on T cell activation, differentiation and INF-γ production itself. This indicates that during trypanosome infections, the Ifng+ Il10+ Il21+ Th1 cells, rather than Tregs, assume an immune regulatory role that is needed for dampening inflammation.

Decision level integration of unimodal and multimodal single cell data with scTriangulate

Decisively delineating cell identities from uni- and multimodal single-cell datasets is complicated by diverse modalities, clustering methods, and reference atlases. We describe scTriangulate, a computational framework to mix-and-match multiple clustering results, modalities, associated algorithms, and resolutions to achieve an optimal solution. Rather than ensemble approaches which select the "consensus", scTriangulate picks the most stable solution through coalitional iteration. When evaluated on diverse multimodal technologies, scTriangulate outperforms alternative approaches to identify high-confidence cell-populations and modality-specific subtypes. Unlike existing integration strategies that rely on modality-specific joint embedding or geometric graphs, scTriangulate makes no assumption about the distributions of raw underlying values. As a result, this approach can solve unprecedented integration challenges, including the ability to automate reference cell-atlas construction, resolve clonal architecture within molecularly defined cell-populations and subdivide clusters to discover splicing-defined disease subtypes. scTriangulate is a flexible strategy for unified integration of single-cell or multimodal clustering solutions, from nearly unlimited sources.

Perspectives on single-nucleus RNA sequencing in different cell types and tissues

Single-cell RNA sequencing has become a powerful and essential tool for delineating cellular diversity in normal tissues and alterations in disease states. For certain cell types and conditions, there are difficulties in isolating intact cells for transcriptome profiling due to their fragility, large size, tight interconnections, and other factors. Single-nucleus RNA sequencing (snRNA-seq) is an alternative or complementary approach for cells that are difficult to isolate. In this review, we will provide an overview of the experimental and analysis steps of snRNA-seq to understand the methods and characteristics of general and tissue-specific snRNA-seq data. Knowing the advantages and limitations of snRNA-seq will increase its use and improve the biological interpretation of the data generated using this technique.

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.

Temporal molecular program of human hematopoietic stem and progenitor cells after birth

Hematopoietic stem and progenitor cells (HSPCs) give rise to the blood system and maintain hematopoiesis throughout the human lifespan. Here, we report a transcriptional census of human bone-marrow-derived HSPCs from the neonate, infant, child, adult, and aging stages, showing two subpopulations of multipotent progenitors separated by CD52 expression. From birth to the adult stage, stem and multipotent progenitors shared similar transcriptional alterations, and erythroid potential was enhanced after the infant stage. By integrating transcriptome, chromatin accessibility, and functional data, we further showed that aging hematopoietic stem cells (HSCs) exhibited a bias toward megakaryocytic differentiation. Finally, in comparison with the HSCs from the cord blood, neonate bone-marrow-derived HSCs were more quiescent and had higher long-term regeneration capability and durable self-renewal. Taken together, this work provides an integral transcriptome landscape of HSPCs and identifies their dynamics in post-natal steady-state hemopoiesis, thereby helping explore hematopoiesis in development and diseases.

Dissecting CD8+ T cell pathology of severe SARS-CoV-2 infection by single-cell immunoprofiling

Introduction

SARS-CoV-2 infection results in varying disease severity, ranging from asymptomatic infection to severe illness. A detailed understanding of the immune response to SARS-CoV-2 is critical to unravel the causative factors underlying differences in disease severity and to develop optimal vaccines against new SARS-CoV-2 variants.

Methods

We combined single-cell RNA and T cell receptor sequencing with CITE-seq antibodies to characterize the CD8⁺ T cell response to SARS-CoV-2 infection at high resolution and compared responses between mild and severe COVID-19.

Results

We observed increased CD8⁺ T cell exhaustion in severe SARS-CoV-2 infection and identified a population of NK-like, terminally differentiated CD8⁺ effector T cells characterized by expression of FCGR3A (encoding CD16). Further characterization of NK-like CD8⁺ T cells revealed heterogeneity among CD16⁺ NK-like CD8⁺ T cells and profound differences in cytotoxicity, exhaustion, and NK-like differentiation between mild and severe disease conditions.

Discussion

We propose a model in which differences in the surrounding inflammatory milieu lead to crucial differences in NK-like differentiation of CD8⁺ effector T cells, ultimately resulting in the appearance of NK-like CD8⁺ T cell populations of different functionality and pathogenicity. Our in-depth characterization of the CD8⁺ T cell-mediated response to SARS-CoV-2 infection provides a basis for further investigation of the importance of NK-like CD8⁺ T cells in COVID-19 severity.

The Value of a Comprehensive Genomic Evaluation in Prenatal Diagnosis of Genetic Diseases: A Retrospective Study

Currently, there are still many challenges in prenatal diagnosis, such as limited or uncertain fetal phenotyping, variant interpretation, and rapid turnaround times. The aim of this study was to illustrate the value of a comprehensive genomic evaluation in prenatal diagnosis. We retrospectively reviewed 20 fetuses with clinically significant copy number variants (CNVs) detected by chromosomal microarray analysis (CMA) and no further exome sequencing testing in our tertiary center between 2019 and 2020. The residual DNA from the prenatal cases was used for the parallel implementation of CNV sequencing (CNV-seq) and trio-based clinical exome sequencing (trio-CES). CMA revealed 26 clinically significant CNVs (18 deletions and eight duplications) in 20 fetuses, in which five fetuses had two or more CNVs. There were eight fetuses with pathogenic CNVs (e.g., del 1p36), nine fetuses with likely pathogenic CNVs (e.g., dup 22q11.21), and three fetuses with variants of unknown significance (VOUS, e.g., dup 1q21.1q21.2). Trio-CES identified four fetuses with likely pathogenic mutations (SNV/InDels). Of note, a fetus was detected with a maternally inherited hemizygous variant in the SLX4 gene due to a 16p13.3 deletion on the paternal chromosome. The sizes of CNVs detected by CNV-seq were slightly larger than that of the SNP array, and four cases with mosaic CNVs were all identified by CNV-seq. In conclusion, microdeletion/duplication syndromes and monogenic disorders may co-exist in a subject, and CNV deletion may contribute to uncovering additional recessive disease alleles. The application of a comprehensive genomic evaluation (CNVs and SNV/InDels) has great value in the prenatal diagnosis arena. CNV-seq based on NGS technology is a reliable and a cost-effective technique for identifying CNVs.

CAMML with the Integration of Marker Proteins (ChIMP)

MOTIVATION

Cell typing is a critical task in the analysis of single-cell data, particularly when studying complex diseased tissues. Unfortunately, the sparsity and noise of single-cell data make accurate cell typing of individual cells difficult. To address these challenges, we previously developed the CAMML method for multi-label cell typing of single-cell RNA-sequencing (scRNA-seq) data. CAMML uses weighted gene sets to score each profiled cell for multiple potential cell types. While CAMML outperforms other scRNA-seq cell typing techniques, it only leverages transcriptomic data so cannot take advantage of newer multi-omic single-cell assays that jointly profile gene expression and protein abundance (e.g. joint scRNA-seq/CITE-seq).

RESULTS

We developed the CAMML with the Integration of Marker Proteins (ChIMP) method to support multi-label cell typing of individual cells jointly profiled via scRNA-seq and CITE-seq. ChIMP combines cell type scores computed on scRNA-seq data via the CAMML approach with discretized CITE-seq measurements for cell type marker proteins. The multi-omic cell type scores generated by ChIMP allow researchers to more precisely and conservatively cell type joint scRNA-seq/CITE-seq data.

AVAILABILITY AND IMPLEMENTATION

An implementation of this work is available on CRAN at https://cran.r-project.org/web/packages/CAMML/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Comprehensive Transcriptomic Analysis of VISTA in Acute Myeloid Leukemia: Insights into Its Prognostic Value

The V-domain Ig suppressor of T-cell activation (VISTA) has been recognized as a critical negative regulator of antitumor immune response and is gaining growing interest as a potential pharmacological target in immunotherapy. This molecule is highly expressed in hematopoietic stem cells and myeloid compartment, and it has been found upmodulated in acute myeloid leukemia (AML). However, VISTA-associated immune features are relatively unexplored in myeloid malignancies. Herein, we aimed to explore whether this immune checkpoint regulator could play a role in the generation of an immune escape environment in AML patients. We characterized VISTA mRNA expression levels in leukemia cell lines and in large publicly available cohorts of specimens from bone marrow of healthy individuals and AML patients at diagnosis by deploying bulk and single-cell RNA sequencing. We also defined the correlations with leukemia-associated burden using results of whole-exome sequencing of AML samples at disease onset. We showed that VISTA expression linearly increased across the myeloid differentiation tree in normal hematopoiesis. Accordingly, its transcript was highly enriched in AML cell lines as well as in AML patients at diagnosis presenting with myelomonocytic and monocytic differentiation. A strong correlation was seen with NPM1 mutations regardless of the presence of FLT3 lesions. Furthermore, VISTA expression levels at baseline correlated with disease recurrence in patients with normal karyotype and NPM1 mutations, a subgroup traditionally considered as favorable according to current diagnostic schemes. Indeed, when compared to patients with long-term remission (>5 years after standard chemotherapy regimens), cases relapsing within 2 years from diagnosis had increased VISTA expression in both leukemia and T cells. Our results suggest a rationale for developing VISTA-targeted therapeutic strategies to treat molecularly defined subgroups of AML patients to prevent disease recurrence and treatment resistance.

Advances in Single-Cell Sequencing Technology and Its Application in Poultry Science

Single-cell sequencing (SCS) uses a single cell as the research material and involves three dimensions: genes, phenotypes and cell biological mechanisms. This type of research can locate target cells, analyze the dynamic changes in the target cells and the relationships between the cells, and pinpoint the molecular mechanism of cell formation. Currently, a common problem faced by animal husbandry scientists is how to apply existing science and technology to promote the production of high-quality livestock and poultry products and to breed livestock for disease resistance; this is also a bottleneck for the sustainable development of animal husbandry. In recent years, although SCS technology has been successfully applied in the fields of medicine and bioscience, its application in poultry science has been rarely reported. With the sustainable development of science and technology and the poultry industry, SCS technology has great potential in the application of poultry science (or animal husbandry). Therefore, it is necessary to review the innovation of SCS technology and its application in poultry science. This article summarizes the current main technical methods of SCS and its application in poultry, which can provide potential references for its future applications in precision breeding, disease prevention and control, immunity, and cell identification.

Single-cell immune profiling reveals thymus-seeding populations, T cell commitment, and multilineage development in the human thymus

T cell development in the mouse thymus has been studied extensively, but less is known regarding T cell development in the human thymus. We used a combination of single-cell techniques and functional assays to perform deep immune profiling of human T cell development, focusing on the initial stages of prelineage commitment. We identified three thymus-seeding progenitor populations that also have counterparts in the bone marrow. In addition, we found that the human thymus physiologically supports the development of monocytes, dendritic cells, and NK cells, as well as limited development of B cells. These results are an important step toward monitoring and guiding regenerative therapies in patients after hematopoietic stem cell transplantation.

Targeting the EIF2AK1 Signaling Pathway Rescues Red Blood Cell Production in SF3B1-Mutant Myelodysplastic Syndromes With Ringed Sideroblasts

SF3B1 mutations, which occur in 20% of patients with myelodysplastic syndromes (MDS), are the hallmarks of a specific MDS subtype, MDS with ringed sideroblasts (MDS-RS), which is characterized by the accumulation of erythroid precursors in the bone marrow and primarily affects the elderly population. Here, using single-cell technologies and functional validation studies of primary SF3B1-mutant MDS-RS samples, we show that SF3B1 mutations lead to the activation of the EIF2AK1 pathway in response to heme deficiency and that targeting this pathway rescues aberrant erythroid differentiation and enables the red blood cell maturation of MDS-RS erythroblasts. These data support the development of EIF2AK1 inhibitors to overcome transfusion dependency in patients with SF3B1-mutant MDS-RS with impaired red blood cell production.

SIGNIFICANCE

MDS-RS are characterized by significant anemia. Patients with MDS-RS die from a shortage of red blood cells and the side effects of iron overload due to their constant need for transfusions. Our study has implications for the development of therapies to achieve long-lasting hematologic responses. This article is highlighted in the In This Issue feature, p. 476.

Interleukin-1 contributes to clonal expansion and progression of bone marrow fibrosis in JAK2V617F-induced myeloproliferative neoplasm

Chronic inflammation is frequently associated with myeloproliferative neoplasms (MPN), but the role of inflammation in the pathogenesis of MPN remains unclear. Expression of the proinflammatory cytokine interleukin-1 (IL-1) is elevated in patients with MPN as well as in Jak2V617F knock-in mice. Here, we show that genetic deletion of IL-1 receptor 1 (IL-1R1) normalizes peripheral blood counts, reduces splenomegaly and ameliorates bone marrow fibrosis in homozygous Jak2V617F mouse model of myelofibrosis. Deletion of IL-1R1 also significantly reduces Jak2V617F mutant hematopoietic stem/progenitor cells. Exogenous administration of IL-1β enhances myeloid cell expansion and accelerates the development of bone marrow fibrosis in heterozygous Jak2V617F mice. Furthermore, treatment with anti-IL-1R1 antibodies significantly reduces leukocytosis and splenomegaly, and ameliorates bone marrow fibrosis in homozygous Jak2V617F mice. Collectively, these results suggest that IL-1 signaling plays a pathogenic role in MPN disease progression, and targeting of IL-1R1 could be a useful strategy for the treatment of myelofibrosis.

Clonal hematopoiesis of indeterminate potential, DNA methylation, and risk for coronary artery disease

Age-related changes to the genome-wide DNA methylation (DNAm) pattern observed in blood are well-documented. Clonal hematopoiesis of indeterminate potential (CHIP), characterized by the age-related acquisition and expansion of leukemogenic mutations in hematopoietic stem cells (HSCs), is associated with blood cancer and coronary artery disease (CAD). Epigenetic regulators DNMT3A and TET2 are the two most frequently mutated CHIP genes. Here, we present results from an epigenome-wide association study for CHIP in 582 Cardiovascular Health Study (CHS) participants, with replication in 2655 Atherosclerosis Risk in Communities (ARIC) Study participants. We show that DNMT3A and TET2 CHIP have distinct and directionally opposing genome-wide DNAm association patterns consistent with their regulatory roles, albeit both promoting self-renewal of HSCs. Mendelian randomization analyses indicate that a subset of DNAm alterations associated with these two leading CHIP genes may promote the risk for CAD.

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.

Regulatory Programs of B-cell Activation and Germinal Center Reaction Allow B-ALL Escape from CD19 CAR T-cell Therapy

Chimeric antigen receptor (CAR) T-cell therapy has led to tremendous successes in the treatment of B-cell malignancies. However, a large fraction of treated patients relapse, often with disease expressing reduced levels of the target antigen. Here, we report that exposing CD19+ B-cell acute lymphoblastic leukemia (B-ALL) cells to CD19 CAR T cells reduced CD19 expression within hours. Initially, CD19 CAR T cells caused clustering of CD19 at the T cell-leukemia cell interface followed by CD19 internalization and decreased CD19 surface expression on the B-ALL cells. CD19 expression was then repressed by transcriptional rewiring. Using single-cell RNA sequencing and single-cell assay for transposase-accessible chromatin using sequencing, we demonstrated that a subset of refractory CD19low cells sustained decreased CD19 expression through transcriptional programs of physiologic B-cell activation and germinal center reaction. Inhibiting B-cell activation programs with the Bruton's tyrosine kinase inhibitor ibrutinib increased the cytotoxicity of CD19 CAR T cells without affecting CAR T-cell viability. These results demonstrate transcriptional plasticity as an underlying mechanism of escape from CAR T cells and highlight the importance of combining CAR T-cell therapy with targeted therapies that aim to overcome this plasticity. See related Spotlight by Zhao and Melenhorst, p. 1040.

HMGA2 expression defines a subset of human AML with immature transcriptional signature and vulnerability to G2/M inhibition

High-mobility group AT-hook 2 (HMGA2) is a nonhistone chromatin-binding protein that is normally expressed in stem cells of various tissues and aberrantly detected in several tumor types. We recently observed that one-fourth of human acute myeloid leukemia (AML) specimens express HMGA2, which associates with a very poor prognosis. We present results indicating that HMGA2+ AMLs share a distinct transcriptional signature representing an immature phenotype. Using single-cell analyses, we showed that HMGA2 is expressed in CD34+ subsets of stem cells and early progenitors, whether normal or derived from AML specimens. Of interest, we found that one of the strongest gene expression signatures associated with HMGA2 in AML is the upregulation of G2/M checkpoint genes. Whole-genome CRISPR/Cas9 screening in HMGA2 overexpressing cells further revealed a synthetic lethal interaction with several G2/M checkpoint genes. Accordingly, small molecules that target G2/M proteins were preferentially active in vitro and in vivo on HMGA2+ AML specimens. Together, our findings suggest that HMGA2 is a key functional determinant in AML and is associated with stem cell features, G2/M status, and related drug sensitivity.

CST6 suppresses osteolytic bone disease in multiple myeloma by blocking osteoclast differentiation

Osteolytic bone disease is a hallmark of multiple myeloma (MM). A significant fraction (~20%) of MM patients do not develop osteolytic lesions (OL). The molecular basis for the absence of bone disease in MM is not understood. We combined PET-CT and gene expression profiling (GEP) of purified bone marrow (BM) CD138+ MM cells from 512 newly diagnosed MM patients to reveal that elevated expression of cystatin M/E (CST6) was significantly associated with the absence of OL in MM. An enzyme-linked immunosorbent assay revealed a strong correlation between CST6 levels in BM serum/plasma and CST6 mRNA expression. Both recombinant CST6 protein and BM serum from patients with high CST6 significantly inhibited the activity of the osteoclast-specific protease cathepsin K, and blocked osteoclast differentiation and function. Recombinant CST6 inhibited bone destruction in ex vivo and in vivo myeloma models. Single cell RNA-sequencing identified that CST6 attenuates polarization of monocytes to osteoclast precursors. Furthermore, CST6 protein blocks osteoclast differentiation by suppressing cathepsin-mediated cleavage of NF-κB/p100 and TRAF3 following RANKL stimulation. Secretion by MM cells of CST6, an inhibitor of osteoclast differentiation and function, suppresses osteolytic bone disease in MM and probably other diseases associated with osteoclast-mediated bone loss.

Single-cell RNA-seq of primary bone marrow neutrophils from female and male adult mice

Widespread sex-dimorphism is observed in the mammalian immune system. Consistently, studies have reported sex differences in the transcriptome of immune cells at the bulk level, including neutrophils. Neutrophils are the most abundant cell type in human blood, and they are key components of the innate immune system as they form a first line of defense against pathogens. Neutrophils are produced in the bone marrow, and differentiation and maturation produce distinct neutrophil subpopulations. Thus, single-cell resolution studies are crucial to decipher the biological significance of neutrophil heterogeneity. However, since neutrophils are very RNA-poor, single-cell profiling of these cells has been technically challenging. Here, we generated a single-cell RNA-seq dataset of primary neutrophils from adult female and male mouse bone marrow. After stringent quality control, we found that previously characterized neutrophil subpopulations can be detected in both sexes. Additionally, we confirmed that canonical sex-linked markers are differentially expressed between female and male cells across neutrophil subpopulations. This dataset provides a groundwork for comparative studies on the lifelong transcriptional sexual dimorphism of neutrophils.

Synchrotron-based characterization of arthroprosthetic CoCrMo particles in human bone marrow

Particles released from cobalt-chromium-molybdenum (CoCrMo) alloys are considered common elicitors of chronic inflammatory adverse effects. There is a lack of data demonstrating particle numbers, size distribution and elemental composition of bone marrow resident particles which would allow for implementation of clinically relevant test strategies in bone marrow models at different degrees of exposure. The aim of this study was to investigate metal particle exposure in human periprosthetic bone marrow of three types of arthroplasty implants. Periprosthetic bone marrow sections from eight patients exposed to CoCrMo particles were analyzed via spatially resolved and synchrotron-based nanoscopic X-ray fluorescence imaging. These analyses revealed lognormal particle size distribution patterns predominantly towards the nanoscale. Analyses of particle numbers and normalization to bone marrow volume and bone marrow cell number indicated particle concentrations of up to 1 × 10¹¹ particles/ml bone marrow or 2 × 10⁴ particles/bone marrow cell, respectively. Analyses of elemental ratios of CoCrMo particles showed that particularly the particles' Co content depends on particle size. The obtained data point towards Co release from arthroprosthetic particles in the course of dealloying and degradation processes of larger particles within periprosthetic bone marrow. This is the first study providing data based on metal particle analyses to be used for future in vitro and in vivo studies of possible toxic effects in human bone marrow following exposure to arthroprosthetic CoCrMo particles of different concentration, size, and elemental composition. Graphical abstract.

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.

Pollock: fishing for cell states

Motivation

The use of single-cell methods is expanding at an ever-increasing rate. While there are established algorithms that address cell classification, they are limited in terms of cross platform compatibility, reliance on the availability of a reference dataset and classification interpretability. Here, we introduce Pollock, a suite of algorithms for cell type identification that is compatible with popular single-cell methods and analysis platforms, provides a set of pretrained human cancer reference models, and reports interpretability scores that identify the genes that drive cell type classifications.

Results

Pollock performs comparably to existing classification methods, while offering easily deployable pretrained classification models across a wide variety of tissue and data types. Additionally, it demonstrates utility in immune pan-cancer analysis.

Availability and implementation

Source code and documentation are available at https://github.com/ding-lab/pollock. Pretrained models and datasets are available for download at https://zenodo.org/record/5895221.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

Antigen presentation safeguards the integrity of the hematopoietic stem cell pool

Hematopoietic stem and progenitor cells (HSPCs) are responsible for the production of blood and immune cells. Throughout life, HSPCs acquire oncogenic aberrations that can cause hematological cancers. Although molecular programs maintaining stem cell integrity have been identified, safety mechanisms eliminating malignant HSPCs from the stem cell pool remain poorly characterized. Here, we show that HSPCs constitutively present antigens via major histocompatibility complex class II. The presentation of immunogenic antigens, as occurring during malignant transformation, triggers bidirectional interactions between HSPCs and antigen-specific CD4+ T cells, causing stem cell proliferation, differentiation, and specific exhaustion of aberrant HSPCs. This immunosurveillance mechanism effectively eliminates transformed HSPCs from the hematopoietic system, thereby preventing leukemia onset. Together, our data reveal a bidirectional interaction between HSPCs and CD4+ T cells, demonstrating that HSPCs are not only passive receivers of immunological signals but also actively engage in adaptive immune responses to safeguard the integrity of the stem cell pool.

The journey of neutropoiesis: how complex landscapes in bone marrow guide continuous neutrophil lineage determination

Neutrophils are the most abundant white blood cell, and they differentiate in homeostasis in the bone marrow from hematopoietic stem cells (HSCs) via multiple intermediate progenitor cells into mature cells that enter the circulation. Recent findings support a continuous model of differentiation in the bone marrow of heterogeneous HSCs and progenitor populations. Cell fate decisions at the levels of proliferation and differentiation are enforced through expression of lineage-determining transcription factors and their interactions, which are influenced by intrinsic (intracellular) and extrinsic (extracellular) mechanisms. Neutrophil homeostasis is subjected to positive-feedback loops, stemming from the gut microbiome, as well as negative-feedback loops resulting from the clearance of apoptotic neutrophils by mature macrophages. Finally, the cellular kinetics regarding the replenishing of the mature neutrophil pool is discussed in light of recent contradictory data.

Recent updates on the biological basis of heterogeneity in bone marrow stromal cells/skeletal stem cells

Based on studies over the last several decades, the self-renewing skeletal lineages derived from bone marrow stroma could be an ideal source for skeletal tissue engineering. However, the markers for osteogenic precursors; i.e., bone marrowderived skeletal stem cells (SSCs), in association with other cells of the marrow stroma (bone marrow stromal cells, BMSCs) and their heterogeneous nature both in vivo and in vitro remain to be clarified. This review aims to highlight: i) the importance of distinguishing BMSCs/SSCs from other "mesenchymal stem/stromal cells", and ii) factors that are responsible for their heterogeneity, and how these factors impact on the differentiation potential of SSCs towards bone. The prospective role of SSC enrichment, their expansion and its impact on SSC phenotype is explored. Emphasis has also been given to emerging single cell RNA sequencing approaches in scrutinizing the unique population of SSCs within the BMSC population, along with their committed progeny. Understanding the factors involved in heterogeneity may help researchers to improvise their strategies to isolate, characterize and adopt best culture practices and source identification to develop standard operating protocols for developing reproducible stem cells grafts. However, more scientific understanding of the molecular basis of heterogeneity is warranted that may be obtained from the robust high-throughput functional transcriptomics of single cells or clonal populations.

Hematopoietic differentiation is characterized by a transient peak of entropy at a single-cell level

Background

Mature blood cells arise from hematopoietic stem cells in the bone marrow by a process of differentiation along one of several different lineage trajectories. This is often represented as a series of discrete steps of increasing progenitor cell commitment to a given lineage, but as for differentiation in general, whether the process is instructive or stochastic remains controversial. Here, we examine this question by analyzing single-cell transcriptomic data from human bone marrow cells, assessing cell-to-cell variability along the trajectories of hematopoietic differentiation into four different types of mature blood cells. The instructive model predicts that cells will be following the same sequence of instructions and that there will be minimal variability of gene expression between them throughout the process, while the stochastic model predicts a role for cell-to-cell variability when lineage commitments are being made.

Results

Applying Shannon entropy to measure cell-to-cell variability among human hematopoietic bone marrow cells at the same stage of differentiation, we observed a transient peak of gene expression variability occurring at characteristic points in all hematopoietic differentiation pathways. Strikingly, the genes whose cell-to-cell variation of expression fluctuated the most over the course of a given differentiation trajectory are pathway-specific genes, whereas genes which showed the greatest variation of mean expression are common to all pathways. Finally, we showed that the level of cell-to-cell variation is increased in the most immature compartment of hematopoiesis in myelodysplastic syndromes.

Conclusions

These data suggest that human hematopoietic differentiation could be better conceptualized as a dynamical stochastic process with a transient stage of cellular indetermination, and strongly support the stochastic view of differentiation. They also highlight the need to consider the role of stochastic gene expression in complex physiological processes and pathologies such as cancers, paving the way for possible noise-based therapies through epigenetic regulation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01264-9.

Of vascular defense, hemostasis, cancer, and platelet biology: an evolutionary perspective

We have established considerable expertise in studying the role of platelets in cancer biology. From this expertise, we were keen to recognize the numerous venous-, arterial-, microvascular-, and macrovascular thrombotic events and immunologic disorders are caused by severe, acute-respiratory-syndrome coronavirus 2 (SARS-CoV-2) infections. With this offering, we explore the evolutionary connections that place platelets at the center of hemostasis, immunity, and adaptive phylogeny. Coevolutionary changes have also occurred in vertebrate viruses and their vertebrate hosts that reflect their respective evolutionary interactions. As mammals adapted from aquatic to terrestrial life and the heavy blood loss associated with placentalization-based live birth, platelets evolved phylogenetically from thrombocytes toward higher megakaryocyte-blebbing-based production rates and the lack of nuclei. With no nuclei and robust RNA synthesis, this adaptation may have influenced viral replication to become less efficient after virus particles are engulfed. Human platelets express numerous receptors that bind viral particles, which developed from archetypal origins to initiate aggregation and exocytic-release of thrombo-, immuno-, angiogenic-, growth-, and repair-stimulatory granule contents. Whether by direct, evolutionary, selective pressure, or not, these responses may help to contain virus spread, attract immune cells for eradication, and stimulate angiogenesis, growth, and wound repair after viral damage. Because mammalian and marsupial platelets became smaller and more plate-like their biophysical properties improved in function, which facilitated distribution near vessel walls in fluid-shear fields. This adaptation increased the probability that platelets could then interact with and engulf shedding virus particles. Platelets also generate circulating microvesicles that increase membrane surface-area encounters and mark viral targets. In order to match virus-production rates, billions of platelets are generated and turned over per day to continually provide active defenses and adaptation to suppress the spectrum of evolving threats like SARS-CoV-2.

Elucidating minimal residual disease of paediatric B-cell acute lymphoblastic leukaemia by single-cell analysis

Minimal residual disease that persists after chemotherapy is the most valuable prognostic marker for haematological malignancies and solid cancers. Unfortunately, our understanding of the resistance elicited in minimal residual disease is limited due to the rarity and heterogeneity of the residual cells. Here we generated 161,986 single-cell transcriptomes to analyse the dynamic changes of B-cell acute lymphoblastic leukaemia (B-ALL) at diagnosis, residual and relapse by combining single-cell RNA sequencing and B-cell-receptor sequencing. In contrast to those at diagnosis, the leukaemic cells at relapse tended to shift to poorly differentiated states, whereas the changes in the residual cells were more complicated. Differential analyses highlighted the activation of the hypoxia pathway in residual cells, resistant clones and B-ALL with MLL rearrangement. Both in vitro and in vivo models demonstrated that inhibition of the hypoxia pathway sensitized leukaemic cells to chemotherapy. This single-cell analysis of minimal residual disease opens up an avenue for the identification of potent treatment opportunities for B-ALL.

Single-Cell Technologies to Decipher the Immune Microenvironment in Myeloid Neoplasms: Perspectives and Opportunities

Myeloid neoplasms (MN) are heterogeneous clonal disorders arising from the expansion of hematopoietic stem and progenitor cells. In parallel with genetic and epigenetic dynamics, the immune system plays a critical role in modulating tumorigenesis, evolution and therapeutic resistance at the various stages of disease progression. Single-cell technologies represent powerful tools to assess the cellular composition of the complex tumor ecosystem and its immune environment, to dissect interactions between neoplastic and non-neoplastic components, and to decipher their functional heterogeneity and plasticity. In addition, recent progress in multi-omics approaches provide an unprecedented opportunity to study multiple molecular layers (DNA, RNA, proteins) at the level of single-cell or single cellular clones during disease evolution or in response to therapy. Applying single-cell technologies to MN holds the promise to uncover novel cell subsets or phenotypic states and highlight the connections between clonal evolution and immune escape, which is crucial to fully understand disease progression and therapeutic resistance. This review provides a perspective on the various opportunities and challenges in the field, focusing on key questions in MN research and discussing their translational value, particularly for the development of more efficient immunotherapies.

Diversity of human NK cell developmental pathways defined by single-cell analyses

Human natural killer (NK) and innate lymphoid cells (ILCs) include diverse specialized phenotypic and functional subsets that reflect their roles as innate immune effector cells present in tissue and circulation. In recent years, significant advances have been made in better defining their tissue resident phenotypes, developmental pathways, and phenotypic plasticity. Here we offer a brief review of new insights into human NK cell diversity specifically defined by next generation sequencing and single-cell transcriptomic studies and integrate these into our current models of human NK cell developmental trajectories and mature subsets. These studies highlight both a deeper understanding of innate lymphoid cell differentiation and homeostasis and underscore critical questions that remain outstanding in the field.

Mapping atopic dermatitis and anti-IL-22 response signatures to type 2-low severe neutrophilic asthma

BACKGROUND

Transcriptomic changes in patients who respond clinically to biological therapies may identify responses in other tissues or diseases.

OBJECTIVE

We sought to determine whether a disease signature identified in atopic dermatitis (AD) is seen in adults with severe asthma and whether a transcriptomic signature for patients with AD who respond clinically to anti-IL-22 (fezakinumab [FZ]) is enriched in severe asthma.

METHODS

An AD disease signature was obtained from analysis of differentially expressed genes between AD lesional and nonlesional skin biopsies. Differentially expressed genes from lesional skin from therapeutic superresponders before and after 12 weeks of FZ treatment defined the FZ-response signature. Gene set variation analysis was used to produce enrichment scores of AD and FZ-response signatures in the Unbiased Biomarkers for the Prediction of Respiratory Disease Outcomes asthma cohort.

RESULTS

The AD disease signature (112 upregulated genes) encompassing inflammatory, T-cell, TH2, and TH17/TH22 pathways was enriched in the blood and sputum of patients with asthma with increasing severity. Patients with asthma with sputum neutrophilia and mixed granulocyte phenotypes were the most enriched (P < .05). The FZ-response signature (296 downregulated genes) was enriched in asthmatic blood (P < .05) and particularly in neutrophilic and mixed granulocytic sputum (P < .05). These data were confirmed in sputum of the Airway Disease Endotyping for Personalized Therapeutics cohort. IL-22 mRNA across tissues did not correlate with FZ-response enrichment scores, but this response signature correlated with TH22/IL-22 pathways.

CONCLUSIONS

The FZ-response signature in AD identifies severe neutrophilic asthmatic patients as potential responders to FZ therapy. This approach will help identify patients for future asthma clinical trials of drugs used successfully in other chronic diseases.

CAMML: Multi-Label Immune Cell-Typing and Stemness Analysis for Single-Cell RNA-sequencing

Inferring the cell types in single-cell RNA-sequencing (scRNA-seq) data is of particular importance for understanding the potential cellular mechanisms and phenotypes occurring in complex tissues, such as the tumor-immune microenvironment (TME). The sparsity and noise of scRNA-seq data, combined with the fact that immune cell types often occur on a continuum, make cell typing of TME scRNA-seq data a significant challenge. Several single-label cell typing methods have been put forth to address the limitations of noise and sparsity, but accounting for the often overlapped spectrum of cell types in the immune TME remains an obstacle. To address this, we developed a new scRNA-seq cell-typing method, Cell-typing using variance Adjusted Mahalanobis distances with Multi-Labeling (CAMML). CAMML leverages cell type-specific weighted gene sets to score every cell in a dataset for every potential cell type. This allows cells to be labelled either by their highest scoring cell type as a single label classification or based on a score cut-off to give multi-label classification. For single-label cell typing, CAMML performance is comparable to existing cell typing methods, SingleR and Garnett. For scenarios where cells may exhibit features of multiple cell types (e.g., undifferentiated cells), the multi-label classification supported by CAMML offers important benefits relative to the current state-of-the-art methods. By integrating data across studies, omics platforms, and species, CAMML serves as a robust and adaptable method for overcoming the challenges of scRNA-seq analysis.

Genetic and Transcriptional Contributions to Relapse in Normal Karyotype Acute Myeloid Leukemia

To better understand clonal and transcriptional adaptations after relapse in patients with acute myeloid leukemia (AML), we collected presentation and relapse samples from six normal karyotype AML cases. We performed enhanced whole-genome sequencing to characterize clonal evolution, and deep-coverage single-cell RNA sequencing on the same samples, which yielded 142,642 high-quality cells for analysis. Identifying expressed mutations in individual cells enabled us to discriminate between normal and AML cells, to identify coordinated changes in the genome and transcriptome, and to identify subclone-specific cell states. We quantified the coevolution of genetic and transcriptional heterogeneity during AML progression, and found that transcriptional changes were significantly correlated with genetic changes. However, transcriptional adaptation sometimes occurred independently, suggesting that clonal evolution does not represent all relevant biological changes. In three cases, we identified cells at diagnosis that likely seeded the relapse. Finally, these data revealed a conserved relapse-enriched leukemic cell state bearing markers of stemness, quiescence, and adhesion. SIGNIFICANCE: These data enabled us to identify a relapse-enriched leukemic cell state with distinct transcriptional properties. Detailed case-by-case analyses elucidated the complex ways in which the AML genome, transcriptome, and immune microenvironment interact to evade chemotherapy. These analyses provide a blueprint for evaluating these factors in larger cohorts.This article is highlighted in the In This Issue feature, p. 1.