Human bone marrow assessment by single-cell RNA sequencing, mass cytometry, and flow cytometry

RAC1 inhibition ameliorates IBSP-induced bone metastasis in lung adenocarcinoma

Macrophage-to-osteoclast differentiation (osteoclastogenesis) plays an essential role in tumor osteolytic bone metastasis (BM), while its specific mechanisms remain largely uncertain in lung adenocarcinoma BM. In this study, we demonstrate that integrin-binding sialoprotein (IBSP), which is highly expressed in the cancer cells from bone metastatic and primary lesions of patients with lung adenocarcinoma, can facilitate BM and directly promote macrophage-to-osteoclast differentiation independent of RANKL/M-CSF. In vivo results further suggest that osteolytic BM in lung cancer specifically relies on IBSP-induced macrophage-to-osteoclast differentiation. Mechanistically, IBSP regulates the Rac family small GTPase 1 (Rac1)-NFAT signaling pathway and mediates the forward shift of macrophage-to-osteoclast differentiation, thereby leading to early osteolysis. Moreover, inhibition of Rac1 by EHT-1864 or azathioprine in mice models can remarkably alleviate IBSP-induced BM of lung cancer. Overall, our study suggests that tumor-secreted IBSP promotes BM by inducing macrophage-to-osteoclast differentiation, with potential as an early diagnostic maker for BM, and Rac1 can be the therapeutic target for IBSP-promoted BM in lung cancer.

Global burden and trends of leukemia attributable to high body mass index risk in adults over the past 30 years

Background

High BMI (Body Mass Index) is a significant factor impacting health, with a clear link to an increased risk of leukemia. Research on this topic is limited. Understanding the epidemiological trends of leukemia attributable to high BMI risk is crucial for disease prevention and patient support.

Methods

We obtained the data from the Global Burden of Disease Study, analyzing the ASR (age-standardized rates), including ASDR (age-standardized death rate) and age-standardized disability-adjusted life years (DALYs) rate, and estimated annual percentage change (EAPC) by gender, age, country, and region from 1990 to 2019.

Results

In 2019, deaths and DALYs have significantly increased to 21.73 thousand and 584.09 thousand. The global age-standardized death and DALYs rates have slightly increased over the past 30 years (EAPCs: 0.34 and 0.29). Among four common leukemia subtypes, only CML (Chronic Myeloid Leukemia) exhibited a significant decrease in ASDR and age-standardized DALYs rate, with EAPC of -1.74 and -1.52. AML (Acute Myeloid Leukemia) showed the most pronounced upward trend in ASDR, with an EAPC of 1.34. These trends vary by gender, age, region, and national economic status. Older people have been at a significantly greater risk. Females globally have borne a higher burden. While males have shown an increasing trend. The regions experiencing the greatest growth in ASR were South Asia. The countries with the largest increases were Equatorial Guinea. However, It is worth noting that there may be variations among specific subtypes of leukemia. Regions with high Socio-demographic Index (SDI) have had the highest ASR, while low-middle SDI regions have shown the greatest increase in these rates. All ASRs values have been positively correlated with SDI, but there has been a turning point in medium to high SDI regions.

Conclusions

Leukemia attributable to high BMI risk is gradually becoming a heavier burden globally. Different subtypes of leukemia have distinct temporal and regional patterns. This study's findings will provide information for analyzing the worldwide disease burden patterns and serve as a basis for disease prevention, developing suitable strategies for the modifiable risk factor.

Single-cell genomics-based immune and disease monitoring in blood malignancies

Achieving long-term disease control using therapeutic immunomodulation is a long-standing concept with a strong tradition in blood malignancies. Besides allogeneic hematopoietic stem cell transplantation that continues to provide potentially curative treatment for otherwise challenging diagnoses, recent years have seen impressive progress in immunotherapies for leukemias and lymphomas with immune checkpoint blockade, bispecific monoclonal antibodies, and CAR T cell therapies. Despite their success, non-response, relapse, and immune toxicities remain frequent, thus prioritizing the elucidation of the underlying mechanisms and identifying predictive biomarkers. The increasing availability of single-cell genomic tools now provides a system's immunology view to resolve the molecular and cellular mechanisms of immunotherapies at unprecedented resolution. Here, we review recent studies that leverage these technological advancements for tracking immune responses, the emergence of immune resistance, and toxicities. As single-cell immune monitoring tools evolve and become more accessible, we expect their wide adoption for routine clinical applications to catalyze more precise therapeutic steering of personal immune responses.

Non-viral DNA delivery and TALEN editing correct the sickle cell mutation in hematopoietic stem cells

Sickle cell disease is a devastating blood disorder that originates from a single point mutation in the HBB gene coding for hemoglobin. Here, we develop a GMP-compatible TALEN-mediated gene editing process enabling efficient HBB correction via a DNA repair template while minimizing risks associated with HBB inactivation. Comparing viral versus non-viral DNA repair template delivery in hematopoietic stem and progenitor cells in vitro, both strategies achieve comparable HBB correction and result in over 50% expression of normal adult hemoglobin in red blood cells without inducing β-thalassemic phenotype. In an immunodeficient female mouse model, transplanted cells edited with the non-viral strategy exhibit higher engraftment and gene correction levels compared to those edited with the viral strategy. Transcriptomic analysis reveals that non-viral DNA repair template delivery mitigates P53-mediated toxicity and preserves high levels of long-term hematopoietic stem cells. This work paves the way for TALEN-based autologous gene therapy for sickle cell disease.

Mapping the cellular biogeography of human bone marrow niches using single-cell transcriptomics and proteomic imaging

Non-hematopoietic cells are essential contributors to hematopoiesis. However, heterogeneity and spatial organization of these cells in human bone marrow remain largely uncharacterized. We used single-cell RNA sequencing (scRNA-seq) to profile 29,325 non-hematopoietic cells and discovered nine transcriptionally distinct subtypes. We simultaneously profiled 53,417 hematopoietic cells and predicted their interactions with non-hematopoietic subsets. We employed co-detection by indexing (CODEX) to spatially profile over 1.2 million cells. We integrated scRNA-seq and CODEX data to link predicted cellular signaling with spatial proximity. Our analysis revealed a hyperoxygenated arterio-endosteal neighborhood for early myelopoiesis, and an adipocytic localization for early hematopoietic stem and progenitor cells (HSPCs). We used our CODEX atlas to annotate new images and uncovered mesenchymal stromal cell (MSC) expansion and spatial neighborhoods co-enriched for leukemic blasts and MSCs in acute myeloid leukemia (AML) patient samples. This spatially resolved, multiomic atlas of human bone marrow provides a reference for investigation of cellular interactions that drive hematopoiesis.

Single-cell RNA sequencing reveals the expansion of circulating tissue-homing B cell subsets in secondary acute dengue viral infection

The roles of antibodies secreted by subsets of B cells in dengue virus (DENV) infection have been extensively studied, yet, the contribution of tissue-homing B cells to antiviral immunity remains unclear. In this study, we performed a comprehensive analysis of B cell subpopulations in peripheral blood samples from DENV-infected patients using single-cell RNA-sequencing (scRNA-seq) datasets and flow cytometry. We showed that plasma cells (PCs) and plasmablasts (PBs) were the predominant B cell populations during the acute phase of secondary natural DENV infection, but not in convalescent phase nor in healthy controls. Interestingly, these cells expressed proliferation, adhesion, and tissue-homing genes, including SELPLG, a homing marker of the skin, the initial infected site of DENV. Flow cytometry analysis confirmed a significant upregulation of cell surface expression of a cutaneous lymphocyte-associated antigen (CLA) encoded by SELPLG in PCs and PBs, compared to naive and memory B cells from the same patients. The analysis of an independent single-cell B-cell receptor sequencing (scBCR-seq) dataset of DENV-infected patients revealed that the peripheral blood PCs and PBs exhibited the highest clonal expansion in secondary DENV infection compared to other B cell subsets. These clonally expanded cells also expressed the highest levels of tissue-homing genes, including SELPLG. In addition, by utilizing a public scRNA-seq dataset of SARS-CoV2 infection, we demonstrated the upregulation of several tissue-homing genes in PCs and PBs. Our study provides evidence for the potential roles of tissue-homing B cell subsets in the context of immune responses against viral infections in humans.

Direct comparison of mass cytometry and single-cell RNA sequencing of human peripheral blood mononuclear cells

Single-cell methods offer a high-resolution approach for characterizing cell populations. Many studies rely on single-cell transcriptomics to draw conclusions regarding cell state and behavior, with the underlying assumption that transcriptomic readouts largely parallel their protein counterparts and subsequent activity. However, the relationship between transcriptomic and proteomic measurements is imprecise, and thus datasets that probe the extent of their concordance will be useful to refine such conclusions. Additionally, novel single-cell analysis tools often lack appropriate gold standard datasets for the purposes of assessment. Integrative (combining the two data modalities) and predictive (using one modality to improve results from the other) approaches in particular, would benefit from transcriptomic and proteomic data from the same sample of cells. For these reasons, we performed single-cell RNA sequencing, mass cytometry, and flow cytometry on a split-sample of human peripheral blood mononuclear cells. We directly compare the proportions of specific cell types resolved by each technique, and further describe the extent to which protein and mRNA measurements correlate within distinct cell types.

A single cell framework identifies functionally and molecularly distinct multipotent progenitors in adult human hematopoiesis

Hematopoietic multipotent progenitors (MPPs) regulate blood cell production to appropriately meet the biological demands of the human body. Human MPPs remain ill-defined whereas mouse MPPs have been well characterized with distinct immunophenotypes and lineage potencies. Using multiomic single cell analyses and complementary functional assays, we identified new human MPPs and oligopotent progenitor populations within Lin-CD34+CD38dim/lo adult bone marrow with distinct biomolecular and functional properties. These populations were prospectively isolated based on expression of CD69, CLL1, and CD2 in addition to classical markers like CD90 and CD45RA. We show that within the canonical Lin-CD34+CD38dim/loCD90CD45RA-MPP population, there is a CD69+ MPP with long-term engraftment and multilineage differentiation potential, a CLL1+ myeloid-biased MPP, and a CLL1-CD69-erythroid-biased MPP. We also show that the canonical Lin-CD34+CD38dim/loCD90-CD45RA+ LMPP population can be separated into a CD2+ LMPP with lymphoid and myeloid potential, a CD2-LMPP with high lymphoid potential, and a CLL1+ GMP with minimal lymphoid potential. We used these new HSPC profiles to study human and mouse bone marrow cells and observe limited cell type specific homology between humans and mice and cell type specific changes associated with aging. By identifying and functionally characterizing new adult MPP sub-populations, we provide an updated reference and framework for future studies in human hematopoiesis.

Clonal hematopoiesis and its impact on the aging osteo-hematopoietic niche

Clonal hematopoiesis (CH) defines a premalignant state predominantly found in older persons that increases the risk of developing hematologic malignancies and age-related inflammatory diseases. However, the risk for malignant transformation or non-malignant disorders is variable and difficult to predict, and defining the clinical relevance of specific candidate driver mutations in individual carriers has proved to be challenging. In addition to the cell-intrinsic mechanisms, mutant cells rely on and alter cell-extrinsic factors from the bone marrow (BM) niche, which complicates the prediction of a mutant cell's fate in a shifting pre-malignant microenvironment. Therefore, identifying the insidious and potentially broad impact of driver mutations on supportive niches and immune function in CH aims to understand the subtle differences that enable driver mutations to yield different clinical outcomes. Here, we review the changes in the aging BM niche and the emerging evidence supporting the concept that CH can progressively alter components of the local BM microenvironment. These alterations may have profound implications for the functionality of the osteo-hematopoietic niche and overall bone health, consequently fostering a conducive environment for the continued development and progression of CH. We also provide an overview of the latest technology developments to study the spatiotemporal dependencies in the CH BM niche, ideally in the context of longitudinal studies following CH over time. Finally, we discuss aspects of CH carrier management in clinical practice, based on work from our group and others.

SuperCellCyto: enabling efficient analysis of large scale cytometry datasets

Advancements in cytometry technologies have enabled quantification of up to 50 proteins across millions of cells at single cell resolution. Analysis of cytometry data routinely involves tasks such as data integration, clustering, and dimensionality reduction. While numerous tools exist, many require extensive run times when processing large cytometry data containing millions of cells. Existing solutions, such as random subsampling, are inadequate as they risk excluding rare cell subsets. To address this, we propose SuperCellCyto, an R package that builds on the SuperCell tool which groups highly similar cells into supercells. SuperCellCyto is available on GitHub ( https://github.com/phipsonlab/SuperCellCyto ) and Zenodo ( https://doi.org/10.5281/zenodo.10521294 ).

Comparative assessment of cytometry by time-of-flight and full spectral flow cytometry based on a 33-color antibody panel

Mass cytometry and full spectrum flow cytometry have recently emerged as new promising single cell proteomic analysis tools that can be exploited to decipher the extensive diversity of immune cell repertoires and their implication in human diseases. In this study, we evaluated the performance of mass cytometry against full spectrum flow cytometry using an identical 33-color antibody panel on four healthy individuals. Our data revealed an overall high concordance in the quantification of major immune cell populations between the two platforms using a semi-automated clustering approach. We further showed a strong correlation of cluster assignment when comparing manual and automated clustering. Both comparisons revealed minor disagreements in the quantification and assignment of rare cell subpopulations. Our study showed that both single cell proteomic technologies generate highly overlapping results and substantiate that the choice of technology is not a primary factor for successful biological assessment of cell profiles but must be considered in a broader design framework of clinical studies.

An immunophenotype-coupled transcriptomic atlas of human hematopoietic progenitors

Analysis of the human hematopoietic progenitor compartment is being transformed by single-cell multimodal approaches. Cellular indexing of transcriptomes and epitopes by sequencing (CITE-seq) enables coupled surface protein and transcriptome profiling, thereby revealing genomic programs underlying progenitor states. To perform CITE-seq systematically on primary human bone marrow cells, we used titrations with 266 CITE-seq antibodies (antibody-derived tags) and machine learning to optimize a panel of 132 antibodies. Multimodal analysis resolved >80 stem, progenitor, immune, stromal and transitional cells defined by distinctive surface markers and transcriptomes. This dataset enables flow cytometry solutions for in silico-predicted cell states and identifies dozens of cell surface markers consistently detected across donors spanning race and sex. Finally, aligning annotations from this atlas, we nominate normal marrow equivalents for acute myeloid leukemia stem cell populations that differ in clinical response. This atlas serves as an advanced digital resource for hematopoietic progenitor analyses in human health and disease.

The RTK-RAS signaling pathway is enriched in patients with rare acute myeloid leukemia harboring t(16;21)(p11;q22)/FUS::ERG

Acute myeloid leukemia (AML) with t(16;21)(p11;q22)/FUS::ERG is a rare AML subtype associated with poor prognosis. However, its clinical and molecular features remain poorly defined. We determined the clinicopathological, genomic, and transcriptomic characteristics and outcomes of patients with AML harboring FUS::ERG at our center. Thirty-six AML patients harboring FUS::ERG were identified, with an incidence rate of 0.3%. These patients were characterized by high lactate dehydrogenase levels (median: 838.5 U/L), elevated bone marrow blast counts (median: 71.5%), and a CD56-positive immunophenotype (94.3%). Notably, we found that RTK-RAS GTPase (RAS) pathway genes, including NRAS (33%) and PTPN11 (24%), were frequently mutated in this subtype. Transcriptome analysis revealed enrichment of the phosphatidylinositol-3-kinase-Akt (PI3K-Akt), mitogen-activated protein kinase (MAPK), and RAS signaling pathways and upregulation of BCL2, the target of venetoclax, in FUS::ERG AML compared to RUNX1::RUNX1T1 AML, a more common AML subtype with good prognosis. The median event-free survival in patients with FUS::ERG AML was 11.9 (95% confidence interval [CI]: 9.0-not available [NA]) months and the median overall survival was 18.2 (95% CI: 12.4-NA) months. Allogeneic hematopoietic stem cell transplantation failed to improve outcomes. Overall, the high incidence of RTK-RAS pathway mutations and high expression of BCL2 may indicate promising therapeutic targets in this high-risk AML subset.

Mapping the Cellular Biogeography of Human Bone Marrow Niches Using Single-Cell Transcriptomics and Proteomic Imaging

The bone marrow is the organ responsible for blood production. Diverse non-hematopoietic cells contribute essentially to hematopoiesis. However, these cells and their spatial organization remain largely uncharacterized as they have been technically challenging to study in humans. Here, we used fresh femoral head samples and performed single-cell RNA sequencing (scRNA-Seq) to profile 29,325 enriched non-hematopoietic bone marrow cells and discover nine transcriptionally distinct subtypes. We next employed CO-detection by inDEXing (CODEX) multiplexed imaging of 18 individuals, including both healthy and acute myeloid leukemia (AML) samples, to spatially profile over one million single cells with a novel 53-antibody panel. We discovered a relatively hyperoxygenated arterio-endosteal niche for early myelopoiesis, and an adipocytic, but not endosteal or perivascular, niche for early hematopoietic stem and progenitor cells. We used our atlas to predict cell type labels in new bone marrow images and used these predictions to uncover mesenchymal stromal cell (MSC) expansion and leukemic blast/MSC-enriched spatial neighborhoods in AML patient samples. Our work represents the first comprehensive, spatially-resolved multiomic atlas of human bone marrow and will serve as a reference for future investigation of cellular interactions that drive hematopoiesis.

Monocytes in allo-HSCT with aged donors secrete IL-1/IL-6/TNF to increase the risk of GVHD and damage the aged HSCs

Aging is considered a critical factor of poor prognosis in allogenic hemopoietic stem cell transplantation (allo-HSCT). To elucidate the underlying mechanisms, we comprehensively reintegrated our clinical data from patients after allo-HSCT and public single-cell transcriptomic profile from post-allo-HSCT and healthy individuals, demonstrating that old donors were more prone to acute GVHD (aGVHD) with pronounced inflammation accumulation and worse overall survival (OS). We also found the presence of inflammation-related CXCL2+ HSC subpopulation during aging with significantly enriched pro-inflammatory pathways. Shifting attention to the HSC microenvironment, we deciphered that IL-1/IL-6 and TRAIL (i.e., TNFSF10) ligand‒receptor pair serves as the crucial bridge between CD14/CD16 monocytes and hematopoietic stem/progenitor cells (HSPCs). The profound upregulation of these signaling pathways during aging finally causes HSC dysfunction and lineage-biased differentiation. Our findings provide the theoretical basis for achieving tailored GVHD management and enhancing allo-HSCT regimens efficacy for aged donors.

Interpreting single-cell messages in normal and aberrant hematopoiesis with the Cell Marker Accordion

Single-cell technologies offer a unique opportunity to explore cellular heterogeneity in hematopoiesis, reveal malignant hematopoietic cells with clinically significant features and measure gene signatures linked to pathological pathways. However, reliable identification of cell types is a crucial bottleneck in single-cell analysis. Available databases contain dissimilar nomenclature and non-concurrent marker sets, leading to inconsistent annotations and poor interpretability. Furthermore, current tools focus mostly on physiological cell types, lacking extensive applicability in disease. We developed the Cell Marker Accordion, a user-friendly platform for the automatic annotation and biological interpretation of single-cell populations based on consistency weighted markers. We validated our approach on peripheral blood and bone marrow single-cell datasets, using surface markers and expert-based annotation as the ground truth. In all cases, we significantly improved the accuracy in identifying cell types with respect to any single source database. Moreover, the Cell Marker Accordion can identify disease-critical cells and pathological processes, extracting potential biomarkers in a wide variety of contexts in human and murine single-cell datasets. It characterizes leukemia stem cell subtypes, including therapy-resistant cells in acute myeloid leukemia patients; it identifies malignant plasma cells in multiple myeloma samples; it dissects cell type alterations in splicing factor-mutant cells from myelodysplastic syndrome patients; it discovers activation of innate immunity pathways in bone marrow from mice treated with METTL3 inhibitors. The breadth of these applications elevates the Cell Marker Accordion as a flexible, faithful and standardized tool to annotate and interpret hematopoietic populations in single-cell datasets focused on the study of hematopoietic development and disease.

Prediction model for EBV infection following HLA haploidentical matched hematopoietic stem cell transplantation

AIMS

Allogeneic hematopoietic stem cell transplantation (allo-HSCT) is an effective treatment for hematological malignancies. However, viral infections, particularly EBV infection, frequently occur following allo-HSCT and can result in multi-tissue and organ damage. Due to the lack of effective antiviral drugs, these infections can even progress to post-transplant lymphoproliferative disorders (PTLD), thereby impacting the prognosis. In light of this, our objective is to develop a prediction model for EBV infection following allo-HSCT.

METHODS

A total of 466 patients who underwent haploidentical hematopoietic stem cell transplantation (haplo-HSCT) between September 2019 and December 2020 were included in this study. The patients were divided into a development cohort and a validation cohort based on the timing of their transplantation. Our aim was to develop and validate a grading scale using these cohorts to predict the risk of EBV infection within the first year after haplo-HSCT. Additionally, single-cell RNA sequencing (sc-RNAseq) data from the bone marrow of healthy donors were utilized to assess the impact of age on immune cells and viral infection.

RESULTS

In the multivariate logistic regression model, four predictors were retained: donor age, female-to-male transplant, graft MNC (mononuclear cell) dose, and CD8 dose. Based on these predictors, an EBV reactivation predicting score system was constructed. The scoring system demonstrated good calibration in both the derivation and validation cohorts, as confirmed by the Hosmer-Lemeshow test (p > 0.05). The scoring system also exhibited favorable discriminative ability, as indicated by the C statistics of 0.72 in the derivation cohort and 0.60 in the validation cohort. Furthermore, the clinical efficacy of the scoring system was evaluated using Kaplan-Meier curves based on risk ratings. The results showed significant differences in EBV reactivation rates between different risk groups, with p-values less than 0.001 in both the derivation and validation cohorts, indicating robust clinical utility. The analysis of sc-RNAseq data from the bone marrow of healthy donors revealed that older age had a profound impact on the quantity and quality of immune subsets. Functional enrichment analysis highlighted that older age was associated with a higher risk of infection. Specifically, CD8 + T cells from older individuals showed enrichment in the pathway of "viral carcinogenesis", while older CD14 + monocytes exhibited enrichment in the pathway of "regulation of viral entry into host cell." These findings suggest that older age may contribute to an increased susceptibility to viral infections, as evidenced by the altered immune profiles observed in the sc-RNAseq data.

CONCLUSION

Overall, these results demonstrate the development and validation of an effective scoring system for predicting EBV reactivation after haplo-HSCT, and provide insights into the impact of age on immune subsets and viral infection susceptibility based on sc-RNAseq analysis of healthy donors' bone marrow.

Cellular hierarchy insights reveal leukemic stem-like cells and early death risk in acute promyelocytic leukemia

Acute promyelocytic leukemia (APL) represents a paradigm for targeted differentiation therapy, with a minority of patients experiencing treatment failure and even early death. We here report a comprehensive single-cell analysis of 16 APL patients, uncovering cellular compositions and their impact on all-trans retinoic acid (ATRA) response in vivo and early death. We unveil a cellular differentiation hierarchy within APL blasts, rooted in leukemic stem-like cells. The oncogenic PML/RARα fusion protein exerts branch-specific regulation in the APL trajectory, including stem-like cells. APL cohort analysis establishes an association of leukemic stemness with elevated white blood cell counts and FLT3-ITD mutations. Furthermore, we construct an APL-specific stemness score, which proves effective in assessing early death risk. Finally, we show that ATRA induces differentiation of primitive blasts and patients with early death exhibit distinct stemness-associated transcriptional programs. Our work provides a thorough survey of APL cellular hierarchies, offering insights into cellular dynamics during targeted therapy.

Single-cell analysis of immune and stroma cell remodeling in clear cell renal cell carcinoma primary tumors and bone metastatic lesions

BACKGROUND

Despite therapeutic advances, once a cancer has metastasized to the bone, it represents a highly morbid and lethal disease. One third of patients with advanced clear cell renal cell carcinoma (ccRCC) present with bone metastasis at the time of diagnosis. However, the bone metastatic niche in humans, including the immune and stromal microenvironments, has not been well-defined, hindering progress towards identification of therapeutic targets.

METHODS

We collected fresh patient samples and performed single-cell transcriptomic profiling of solid metastatic tissue (Bone Met), liquid bone marrow at the vertebral level of spinal cord compression (Involved), and liquid bone marrow from a different vertebral body distant from the tumor site but within the surgical field (Distal), as well as bone marrow from patients undergoing hip replacement surgery (Benign). In addition, we incorporated single-cell data from primary ccRCC tumors (ccRCC Primary) for comparative analysis.

RESULTS

The bone marrow of metastatic patients is immune-suppressive, featuring increased, exhausted CD8 + cytotoxic T cells, T regulatory cells, and tumor-associated macrophages (TAM) with distinct transcriptional states in metastatic lesions. Bone marrow stroma from tumor samples demonstrated a tumor-associated mesenchymal stromal cell population (TA-MSC) that appears to be supportive of epithelial-to mesenchymal transition (EMT), bone remodeling, and a cancer-associated fibroblast (CAFs) phenotype. This stromal subset is associated with poor progression-free and overall survival and also markedly upregulates bone remodeling through the dysregulation of RANK/RANKL/OPG signaling activity in bone cells, ultimately leading to bone resorption.

CONCLUSIONS

These results provide a comprehensive analysis of the bone marrow niche in the setting of human metastatic cancer and highlight potential therapeutic targets for both cell populations and communication channels.

Genomic landscape of patients with germline RUNX1 variants and familial platelet disorder with myeloid malignancy

ABSTRACT

Familial platelet disorder with associated myeloid malignancies (FPDMM) is caused by germline RUNX1 mutations and characterized by thrombocytopenia and increased risk of hematologic malignancies. We recently launched a longitudinal natural history study for patients with FPDMM. Among 27 families with research genomic data by the end of 2021, 26 different germline RUNX1 variants were detected. Besides missense mutations enriched in Runt homology domain and loss-of-function mutations distributed throughout the gene, splice-region mutations and large deletions were detected in 6 and 7 families, respectively. In 25 of 51 (49%) patients without hematologic malignancy, somatic mutations were detected in at least 1 of the clonal hematopoiesis of indeterminate potential (CHIP) genes or acute myeloid leukemia (AML) driver genes. BCOR was the most frequently mutated gene (in 9 patients), and multiple BCOR mutations were identified in 4 patients. Mutations in 6 other CHIP- or AML-driver genes (TET2, DNMT3A, KRAS, LRP1B, IDH1, and KMT2C) were also found in ≥2 patients without hematologic malignancy. Moreover, 3 unrelated patients (1 with myeloid malignancy) carried somatic mutations in NFE2, which regulates erythroid and megakaryocytic differentiation. Sequential sequencing data from 19 patients demonstrated dynamic changes of somatic mutations over time, and stable clones were more frequently found in older adult patients. In summary, there are diverse types of germline RUNX1 mutations and high frequency of somatic mutations related to clonal hematopoiesis in patients with FPDMM. Monitoring changes in somatic mutations and clinical manifestations prospectively may reveal mechanisms for malignant progression and inform clinical management. This trial was registered at www.clinicaltrials.gov as #NCT03854318.

A rigorous benchmarking of alignment-based HLA typing algorithms for RNA-seq data

Accurate identification of human leukocyte antigen (HLA) alleles is essential for various clinical and research applications, such as transplant matching and drug sensitivities. Recent advances in RNA-seq technology have made it possible to impute HLA types from sequencing data, spurring the development of a large number of computational HLA typing tools. However, the relative performance of these tools is unknown, limiting the ability for clinical and biomedical research to make informed choices regarding which tools to use. Here we report the study design of a comprehensive benchmarking of the performance of 12 HLA callers across 682 RNA-seq samples from 8 datasets with molecularly defined gold standard at 5 loci, HLA-A, -B, -C, -DRB1, and -DQB1. For each HLA typing tool, we will comprehensively assess their accuracy, compare default with optimized parameters, and examine for discrepancies in accuracy at the allele and loci levels. We will also evaluate the computational expense of each HLA caller measured in terms of CPU time and RAM. We also plan to evaluate the influence of read length over the HLA region on accuracy for each tool. Most notably, we will examine the performance of HLA callers across European and African groups, to determine discrepancies in accuracy associated with ancestry. We hypothesize that RNA-Seq HLA callers are capable of returning high-quality results, but the tools that offer a good balance between accuracy and computational expensiveness for all ancestry groups are yet to be developed. We believe that our study will provide clinicians and researchers with clear guidance to inform their selection of an appropriate HLA caller.

Single-cell profiling of the microenvironment in human bone metastatic renal cell carcinoma

Bone metastasis is of common occurrence in renal cell carcinoma with poor prognosis, but no optimal treatment approach has been established for bone metastatic renal cell carcinoma. To explore the potential therapeutic targets for bone metastatic renal cell carcinoma, we profile single cell transcriptomes of 6 primary renal cell carcinoma and 9 bone metastatic renal cell carcinoma. We also include scRNA-seq data of early-stage renal cell carcinoma, late-stage renal cell carcinoma, normal kidneys and healthy bone marrow samples in the study to better understand the bone metastasis niche. The molecular properties and dynamic changes of major cell lineages in bone metastatic environment of renal cell carcinoma are characterized. Bone metastatic renal cell carcinoma is associated with multifaceted immune deficiency together with cancer-associated fibroblasts, specifically appearance of macrophages exhibiting malignant and pro-angiogenic features. We also reveal the dominance of immune inhibitory T cells in the bone metastatic renal cell carcinoma which can be partially restored by the treatment. Trajectory analysis showes that myeloid-derived suppressor cells are progenitors of macrophages in the bone metastatic renal cell carcinoma while monocytes are their progenitors in primary tumors and healthy bone marrows. Additionally, the infiltration of immune inhibitory CD47+ T cells is observed in bone metastatic tumors, which may be a result of reduced phagocytosis by SIRPA-expressing macrophages in the bone microenvironment. Together, our results provide a systematic view of various cell types in bone metastatic renal cell carcinoma and suggest avenues for therapeutic solutions.

Precise single-cell transcriptomic mapping of normal and leukemic cell states reveals unconventional lineage priming in acute myeloid leukemia

Initial classification of acute leukemia involves the assignment of blasts to cell states within the hematopoietic hierarchy based on morphological and immunophenotypic features. Yet, these traditional classification approaches lack precision, especially at the level of immature blasts. Single-cell RNA-sequencing (scRNA-seq) enables precise determination of cell state using thousands of markers, thus providing an opportunity to re-examine present-day classification schemes of acute leukemia. Here, we developed a detailed reference map of human bone marrow hematopoiesis from 263,519 single-cell transcriptomes spanning 55 cellular states. Cell state annotations were benchmarked against purified cell populations, and in-depth characterization of gene expression programs underlying hematopoietic differentiation was undertaken. Projection of single-cell transcriptomes from 175 samples spanning acute myeloid leukemia (AML), mixed phenotype acute leukemia (MPAL), and acute erythroid leukemia (AEL) revealed 11 subtypes involving distinct stages of hematopoietic differentiation. These included AML subtypes with notable lymphoid or erythroid lineage priming, challenging traditional diagnostic boundaries between AML, MPAL, and AEL. Quantification of lineage priming in bulk patient cohorts revealed specific genetic alterations associated with this unconventional lineage priming. Integration of transcriptional and genetic information at the single-cell level revealed how genetic subclones can induce lineage restriction, differentiation blocks, or expansion of mature myeloid cells. Furthermore, we demonstrate that distinct cellular hierarchies can co-exist within individual patients, providing insight into AML evolution in response to varying selection pressures. Together, precise mapping of hematopoietic cell states can serve as a foundation for refining disease classification in acute leukemia and understanding response or resistance to emerging therapies.

SINGLE CELL DISSECTION OF DEVELOPMENTAL ORIGINS AND TRANSCRIPTIONAL HETEROGENEITY IN B-CELL ACUTE LYMPHOBLASTIC LEUKEMIA

Sequencing of bulk tumor populations has improved genetic classification and risk assessment of B-ALL, but does not directly examine intratumor heterogeneity or infer leukemia cellular origins. We profiled 89 B-ALL samples by single-cell RNA-seq (scRNA-seq) and compared them to a reference map of normal human B-cell development established using both functional and molecular assays. Intra-sample heterogeneity was driven by cell cycle, metabolism, differentiation, and inflammation transcriptional programs. By inference of B lineage developmental state composition, nearly all samples possessed a high abundance of pro-B cells, with variation between samples mainly driven by sub-populations. However, ZNF384- r and DUX4- r B-ALL showed composition enrichment of hematopoietic stem cells, BCR::ABL1 and KMT2A -r ALL of Early Lymphoid progenitors, MEF2D -r and TCF3::PBX1 of Pre-B cells. Enrichment of Early Lymphoid progenitors correlated with high-risk clinical features. Understanding variation in transcriptional programs and developmental states of B-ALL by scRNA-seq refines existing clinical and genomic classifications and improves prediction of treatment outcome.

scDREAMER for atlas-level integration of single-cell datasets using deep generative model paired with adversarial classifier

Integration of heterogeneous single-cell sequencing datasets generated across multiple tissue locations, time, and conditions is essential for a comprehensive understanding of the cellular states and expression programs underlying complex biological systems. Here, we present scDREAMER ( https://github.com/Zafar-Lab/scDREAMER ), a data-integration framework that employs deep generative models and adversarial training for both unsupervised and supervised (scDREAMER-Sup) integration of multiple batches. Using six real benchmarking datasets, we demonstrate that scDREAMER can overcome critical challenges including skewed cell type distribution among batches, nested batch-effects, large number of batches and conservation of development trajectory across batches. Our experiments also show that scDREAMER and scDREAMER-Sup outperform state-of-the-art unsupervised and supervised integration methods respectively in batch-correction and conservation of biological variation. Using a 1 million cells dataset, we demonstrate that scDREAMER is scalable and can perform atlas-level cross-species (e.g., human and mouse) integration while being faster than other deep-learning-based methods.

Bone marrow vasculature advanced in vitro models for cancer and cardiovascular research

Cardiometabolic diseases and cancer are among the most common diseases worldwide and are a serious concern to the healthcare system. These conditions, apparently distant, share common molecular and cellular determinants, that can represent targets for preventive and therapeutic approaches. The bone marrow plays an important role in this context as it is the main source of cells involved in cardiovascular regeneration, and one of the main sites of liquid and solid tumor metastasis, both characterized by the cellular trafficking across the bone marrow vasculature. The bone marrow vasculature has been widely studied in animal models, however, it is clear the need for human-specific in vitro models, that resemble the bone vasculature lined by endothelial cells to study the molecular mechanisms governing cell trafficking. In this review, we summarized the current knowledge on in vitro models of bone marrow vasculature developed for cardiovascular and cancer research.

Integrated Genomic and Transcriptomic Analysis Improves Disease Classification and Risk Stratification of MDS with Ring Sideroblasts

PURPOSE

Ring sideroblasts (RS) define the low-risk myelodysplastic neoplasm (MDS) subgroup with RS but may also reflect erythroid dysplasia in higher risk myeloid neoplasm. The benign behavior of MDS with RS (MDSRS+) is limited to SF3B1-mutated cases without additional high-risk genetic events, but one third of MDSRS+ carry no SF3B1 mutation, suggesting that different molecular mechanisms may underlie RS formation. We integrated genomic and transcriptomic analyses to evaluate whether transcriptome profiles may improve current risk stratification.

EXPERIMENTAL DESIGN

We studied a prospective cohort of MDSRS+ patients irrespective of World Health Organization (WHO) class with regard to somatic mutations, copy-number alterations, and bone marrow CD34+ cell transcriptomes to assess whether transcriptome profiles add to prognostication and provide input on disease classification.

RESULTS

SF3B1, SRSF2, or TP53 multihit mutations were found in 89% of MDSRS+ cases, and each mutation category was associated with distinct clinical outcome, gene expression, and alternative splicing profiles. Unsupervised clustering analysis identified three clusters with distinct hemopoietic stem and progenitor (HSPC) composition, which only partially overlapped with mutation groups. IPSS-M and the transcriptome-defined proportion of megakaryocyte/erythroid progenitors (MEP) independently predicted survival in multivariable analysis.

CONCLUSIONS

These results provide essential input on the molecular basis of SF3B1-unmutated MDSRS+ and propose HSPC quantification as a prognostic marker in myeloid neoplasms with RS.

ZNF683 marks a CD8+ T cell population associated with anti-tumor immunity following anti-PD-1 therapy for Richter syndrome

Unlike many other hematologic malignancies, Richter syndrome (RS), an aggressive B cell lymphoma originating from indolent chronic lymphocytic leukemia, is responsive to PD-1 blockade. To discover the determinants of response, we analyze single-cell transcriptome data generated from 17 bone marrow samples longitudinally collected from 6 patients with RS. Response is associated with intermediate exhausted CD8 effector/effector memory T cells marked by high expression of the transcription factor ZNF683, determined to be evolving from stem-like memory cells and divergent from terminally exhausted cells. This signature overlaps with that of tumor-infiltrating populations from anti-PD-1 responsive solid tumors. ZNF683 is found to directly target key T cell genes (TCF7, LMO2, CD69) and impact pathways of T cell cytotoxicity and activation. Analysis of pre-treatment peripheral blood from 10 independent patients with RS treated with anti-PD-1, as well as patients with solid tumors treated with anti-PD-1, supports an association of ZNF683high T cells with response.

Evaluation of zero counts to better understand the discrepancies between bulk and single-cell RNA-Seq platforms

Recent advances in sample preparation and sequencing technology have made it possible to profile the transcriptomes of individual cells using single-cell RNA sequencing (scRNA-Seq). Compared to bulk RNA-Seq data, single-cell data often contain a higher percentage of zero reads, mainly due to lower sequencing depth per cell, which affects mostly measurements of low-expression genes. However, discrepancies between platforms are observed regardless of expression level. Using four paired datasets with multiple samples each, we investigated technical and biological factors that can contribute to this expression shift. Using two separate machine learning models we found that, in addition to expression level, RNA integrity, gene or UTR3 length, and the number of transcripts potentially also influence the occurrence of zeros. These findings could enable the development of novel analytical methods for cross-platform expression shift correction. We also identified genes and biological pathways in our diverse datasets that consistently showed differences when assessed at the single cell versus bulk level to assist in interpreting analysis across transcriptomic platforms. At the gene level, 25 genes (0.12%) were found in all datasets as discordant, but at the pathway level, 7 pathways (2.02%) showed shared enrichment in discordant genes.

Multimodal deep learning approaches for single-cell multi-omics data integration

Integrating single-cell multi-omics data is a challenging task that has led to new insights into complex cellular systems. Various computational methods have been proposed to effectively integrate these rapidly accumulating datasets, including deep learning. However, despite the proven success of deep learning in integrating multi-omics data and its better performance over classical computational methods, there has been no systematic study of its application to single-cell multi-omics data integration. To fill this gap, we conducted a literature review to explore the use of multimodal deep learning techniques in single-cell multi-omics data integration, taking into account recent studies from multiple perspectives. Specifically, we first summarized different modalities found in single-cell multi-omics data. We then reviewed current deep learning techniques for processing multimodal data and categorized deep learning-based integration methods for single-cell multi-omics data according to data modality, deep learning architecture, fusion strategy, key tasks and downstream analysis. Finally, we provided insights into using these deep learning models to integrate multi-omics data and better understand single-cell biological mechanisms.

Single-cell reference mapping to construct and extend cell-type hierarchies

Single-cell genomics is now producing an ever-increasing amount of datasets that, when integrated, could provide large-scale reference atlases of tissue in health and disease. Such large-scale atlases increase the scale and generalizability of analyses and enable combining knowledge generated by individual studies. Specifically, individual studies often differ regarding cell annotation terminology and depth, with different groups specializing in different cell type compartments, often using distinct terminology. Understanding how these distinct sets of annotations are related and complement each other would mark a major step towards a consensus-based cell-type annotation reflecting the latest knowledge in the field. Whereas recent computational techniques, referred to as 'reference mapping' methods, facilitate the usage and expansion of existing reference atlases by mapping new datasets (i.e. queries) onto an atlas; a systematic approach towards harmonizing dataset-specific cell-type terminology and annotation depth is still lacking. Here, we present 'treeArches', a framework to automatically build and extend reference atlases while enriching them with an updatable hierarchy of cell-type annotations across different datasets. We demonstrate various use cases for treeArches, from automatically resolving relations between reference and query cell types to identifying unseen cell types absent in the reference, such as disease-associated cell states. We envision treeArches enabling data-driven construction of consensus atlas-level cell-type hierarchies and facilitating efficient usage of reference atlases.

The Immunophenotypic Profile of Healthy Human Bone Marrow

Flow cytometry enables multiparametric characterization of hematopoietic cell immunophenotype. Deviations from normal immunophenotypic patterns comprise a cardinal feature of many hematopoietic neoplasms, underscoring the ongoing essentiality of flow cytometry as a diagnostic tool. However, understanding of aberrant hematopoiesis requires an equal understanding of normal hematopoiesis as a comparator. In this review, we outline key features of healthy adult hematopoiesis and lineage specification as illuminated by flow cytometry and provide diagrams illustrating what a diagnostician may observe in flow cytometric plots. These features provide a profile of baseline hematopoiesis, to which clinical samples with suspected neoplasia may be compared.

Decoding Aging Hallmarks at the Single-Cell Level

Organismal aging exhibits wide-ranging hallmarks in divergent cell types across tissues, organs, and systems. The advancement of single-cell technologies and generation of rich datasets have afforded the scientific community the opportunity to decode these hallmarks of aging at an unprecedented scope and resolution. In this review, we describe the technological advancements and bioinformatic methodologies enabling data interpretation at the cellular level. Then, we outline the application of such technologies for decoding aging hallmarks and potential intervention targets and summarize common themes and context-specific molecular features in representative organ systems across the body. Finally, we provide a brief summary of available databases relevant for aging research and present an outlook on the opportunities in this emerging field.

Quantifying Cell-State Densities in Single-Cell Phenotypic Landscapes using Mellon

Cell-state density characterizes the distribution of cells along phenotypic landscapes and is crucial for unraveling the mechanisms that drive cellular differentiation, regeneration, and disease. Here, we present Mellon, a novel computational algorithm for high-resolution estimation of cell-state densities from single-cell data. We demonstrate Mellon's efficacy by dissecting the density landscape of various differentiating systems, revealing a consistent pattern of high-density regions corresponding to major cell types intertwined with low-density, rare transitory states. Utilizing hematopoietic stem cell fate specification to B-cells as a case study, we present evidence implicating enhancer priming and the activation of master regulators in the emergence of these transitory states. Mellon offers the flexibility to perform temporal interpolation of time-series data, providing a detailed view of cell-state dynamics during the inherently continuous developmental processes. Scalable and adaptable, Mellon facilitates density estimation across various single-cell data modalities, scaling linearly with the number of cells. Our work underscores the importance of cell-state density in understanding the differentiation processes, and the potential of Mellon to provide new insights into the regulatory mechanisms guiding cellular fate decisions.

Single-cell transcriptomics and epigenomics unravel the role of monocytes in neuroblastoma bone marrow metastasis

Metastasis is the major cause of cancer-related deaths. Neuroblastoma (NB), a childhood tumor has been molecularly defined at the primary cancer site, however, the bone marrow (BM) as the metastatic niche of NB is poorly characterized. Here we perform single-cell transcriptomic and epigenomic profiling of BM aspirates from 11 subjects spanning three major NB subtypes and compare these to five age-matched and metastasis-free BM, followed by in-depth single cell analyses of tissue diversity and cell-cell interactions, as well as functional validation. We show that cellular plasticity of NB tumor cells is conserved upon metastasis and tumor cell type composition is NB subtype-dependent. NB cells signal to the BM microenvironment, rewiring via macrophage mgration inhibitory factor and midkine signaling specifically monocytes, which exhibit M1 and M2 features, are marked by activation of pro- and anti-inflammatory programs, and express tumor-promoting factors, reminiscent of tumor-associated macrophages. The interactions and pathways characterized in our study provide the basis for therapeutic approaches that target tumor-to-microenvironment interactions.

Cell-cell interactome of the hematopoietic niche and its changes in acute myeloid leukemia

The bone marrow (BM) is a complex microenvironment, coordinating the production of billions of blood cells every day. Despite its essential role and its relevance to hematopoietic diseases, this environment remains poorly characterized. Here we present a high-resolution characterization of the niche in health and acute myeloid leukemia (AML) by establishing a single-cell gene expression database of 339,381 BM cells. We found significant changes in cell type proportions and gene expression in AML, indicating that the entire niche is disrupted. We then predicted interactions between hematopoietic stem and progenitor cells (HSPCs) and other BM cell types, revealing a remarkable expansion of predicted interactions in AML that promote HSPC-cell adhesion, immunosuppression, and cytokine signaling. In particular, predicted interactions involving transforming growth factor β1 (TGFB1) become widespread, and we show that this can drive AML cell quiescence in vitro. Our results highlight potential mechanisms of enhanced AML-HSPC competitiveness and a skewed microenvironment, fostering AML growth.

Nanobody-tethered transposition enables multifactorial chromatin profiling at single-cell resolution

Chromatin states are functionally defined by a complex combination of histone modifications, transcription factor binding, DNA accessibility and other factors. Current methods for defining chromatin states cannot measure more than one aspect in a single experiment at single-cell resolution. Here we introduce nanobody-tethered transposition followed by sequencing (NTT-seq), an assay capable of measuring the genome-wide presence of up to three histone modifications and protein-DNA binding sites at single-cell resolution. NTT-seq uses recombinant Tn5 transposase fused to a set of secondary nanobodies (nb). Each nb-Tn5 fusion protein specifically binds to different immunoglobulin-G antibodies, enabling a mixture of primary antibodies binding different epitopes to be used in a single experiment. We apply bulk-cell and single-cell NTT-seq to generate high-resolution multimodal maps of chromatin states in cell culture and in human immune cells. We also extend NTT-seq to enable simultaneous profiling of cell surface protein expression and multimodal chromatin states to study cells of the immune system.

CIForm as a Transformer-based model for cell-type annotation of large-scale single-cell RNA-seq data

Single-cell omics technologies have made it possible to analyze the individual cells within a biological sample, providing a more detailed understanding of biological systems. Accurately determining the cell type of each cell is a crucial goal in single-cell RNA-seq (scRNA-seq) analysis. Apart from overcoming the batch effects arising from various factors, single-cell annotation methods also face the challenge of effectively processing large-scale datasets. With the availability of an increase in the scRNA-seq datasets, integrating multiple datasets and addressing batch effects originating from diverse sources are also challenges in cell-type annotation. In this work, to overcome the challenges, we developed a supervised method called CIForm based on the Transformer for cell-type annotation of large-scale scRNA-seq data. To assess the effectiveness and robustness of CIForm, we have compared it with some leading tools on benchmark datasets. Through the systematic comparisons under various cell-type annotation scenarios, we exhibit that the effectiveness of CIForm is particularly pronounced in cell-type annotation. The source code and data are available at https://github.com/zhanglab-wbgcas/CIForm.

Global burden of hematologic malignancies and evolution patterns over the past 30 years

Hematologic malignancies are among the most common cancers, and understanding their incidence and death is crucial for targeting prevention, clinical practice improvement, and research resources appropriately. Here, we investigated detailed information on hematological malignancies for the period 1990-2019 from the Global Burden of Disease study. The age-standardized incidence rate (ASIR), the age-standardized death rate (ASDR), and the corresponding estimated annual percentage changes (EAPC) were calculated to assess temporal trends in 204 countries and territories over the past 30 years. Globally, incident cases of hematologic malignancies have been increasing since 1990, reaching 1343.85 thousand in 2019, but the ASDR for all types of hematologic malignancies has been declining. The ASDR for leukemia, multiple myeloma, non-Hodgkin lymphoma, and Hodgkin lymphoma were 4.26, 1.42, 3.19, and 0.34 per 100,000 population in 2019, respectively, with Hodgkin lymphoma showing the most significant decline. However, the trend varies by gender, age, region, and the country's economic situation. The burden of hematologic malignancies is generally higher in men, and this gender gap decreases after peaking at a given age. The regions with the largest increasing trend in the ASIR of leukemia, multiple myeloma, non-Hodgkin lymphoma, and Hodgkin lymphoma were Central Europe, Eastern Europe, East Asia, and Caribbean, respectively. In addition, the proportion of deaths attributed to high body-mass index continued to rise across regions, especially in regions with high socio-demographic indices (SDI). Meanwhile, the burden of leukemia from occupational exposure to benzene and formaldehyde was more widespread in areas with low SDI. Thus, hematologic malignancies remain the leading cause of the global tumor burden, with growing absolute numbers but sharp among several age-standardized measures over the past three decades. The results of the study will inform analysis of trends in the global burden of disease for specific hematologic malignancies and develop appropriate policies for these modifiable risks.

MASI enables fast model-free standardization and integration of single-cell transcriptomics data

Single-cell transcriptomics datasets from the same anatomical sites generated by different research labs are becoming increasingly common. However, fast and computationally inexpensive tools for standardization of cell-type annotation and data integration are still needed in order to increase research inclusivity. To standardize cell-type annotation and integrate single-cell transcriptomics datasets, we have built a fast model-free integration method, named MASI (Marker-Assisted Standardization and Integration). We benchmark MASI with other well-established methods and demonstrate that MASI outperforms other methods, in terms of integration, annotation, and speed. To harness knowledge from single-cell atlases, we demonstrate three case studies that cover integration across biological conditions, surveyed participants, and research groups, respectively. Finally, we show MASI can annotate approximately one million cells on a personal laptop, making large-scale single-cell data integration more accessible. We envision that MASI can serve as a cheap computational alternative for the single-cell research community.

Immune heterogeneity in cardiovascular diseases from a single-cell perspective

A variety of immune cell subsets occupy different niches in the cardiovascular system, causing changes in the structure and function of the heart and vascular system, and driving the progress of cardiovascular diseases (CVDs). The immune cells infiltrating the injury site are highly diverse and integrate into a broad dynamic immune network that controls the dynamic changes of CVDs. Due to technical limitations, the effects and molecular mechanisms of these dynamic immune networks on CVDs have not been fully revealed. With recent advances in single-cell technologies such as single-cell RNA sequencing, systematic interrogation of the immune cell subsets is feasible and will provide insights into the way we understand the integrative behavior of immune populations. We no longer lightly ignore the role of individual cells, especially certain highly heterogeneous or rare subpopulations. We summarize the phenotypic diversity of immune cell subsets and their significance in three CVDs of atherosclerosis, myocardial ischemia and heart failure. We believe that such a review could enhance our understanding of how immune heterogeneity drives the progression of CVDs, help to elucidate the regulatory roles of immune cell subsets in disease, and thus guide the development of new immunotherapies.

Prediction of HLA genotypes from single-cell transcriptome data

The human leukocyte antigen (HLA) locus plays a central role in adaptive immune function and has significant clinical implications for tissue transplant compatibility and allelic disease associations. Studies using bulk-cell RNA sequencing have demonstrated that HLA transcription may be regulated in an allele-specific manner and single-cell RNA sequencing (scRNA-seq) has the potential to better characterize these expression patterns. However, quantification of allele-specific expression (ASE) for HLA loci requires sample-specific reference genotyping due to extensive polymorphism. While genotype prediction from bulk RNA sequencing is well described, the feasibility of predicting HLA genotypes directly from single-cell data is unknown. Here we evaluate and expand upon several computational HLA genotyping tools by comparing predictions from human single-cell data to gold-standard, molecular genotyping. The highest 2-field accuracy averaged across all loci was 76% by arcasHLA and increased to 86% using a composite model of multiple genotyping tools. We also developed a highly accurate model (AUC 0.93) for predicting HLA-DRB345 copy number in order to improve genotyping accuracy of the HLA-DRB locus. Genotyping accuracy improved with read depth and was reproducible at repeat sampling. Using a metanalytic approach, we also show that HLA genotypes from PHLAT and OptiType can generate ASE ratios that are highly correlated (R2 = 0.8 and 0.94, respectively) with those derived from gold-standard genotyping.

Mechanisms of response and resistance to combined decitabine and ipilimumab for advanced myeloid disease

The challenge of eradicating leukemia in patients with acute myelogenous leukemia (AML) after initial cytoreduction has motivated modern efforts to combine synergistic active modalities including immunotherapy. Recently, the ETCTN/CTEP 10026 study tested the combination of the DNA methyltransferase inhibitor decitabine together with the immune checkpoint inhibitor ipilimumab for AML/myelodysplastic syndrome (MDS) either after allogeneic hematopoietic stem cell transplantation (HSCT) or in the HSCT-naïve setting. Integrative transcriptome-based analysis of 304 961 individual marrow-infiltrating cells for 18 of 48 subjects treated on study revealed the strong association of response with a high baseline ratio of T to AML cells. Clinical responses were predominantly driven by decitabine-induced cytoreduction. Evidence of immune activation was only apparent after ipilimumab exposure, which altered CD4+ T-cell gene expression, in line with ongoing T-cell differentiation and increased frequency of marrow-infiltrating regulatory T cells. For post-HSCT samples, relapse could be attributed to insufficient clearing of malignant clones in progenitor cell populations. In contrast to AML/MDS bone marrow, the transcriptomes of leukemia cutis samples from patients with durable remission after ipilimumab monotherapy showed evidence of increased infiltration with antigen-experienced resident memory T cells and higher expression of CTLA-4 and FOXP3. Altogether, activity of combined decitabine and ipilimumab is impacted by cellular expression states within the microenvironmental niche of leukemic cells. The inadequate elimination of leukemic progenitors mandates urgent development of novel approaches for targeting these cell populations to generate long-lasting responses. This trial was registered at www.clinicaltrials.gov as #NCT02890329.

The shaky foundations of simulating single-cell RNA sequencing data

BACKGROUND

With the emergence of hundreds of single-cell RNA-sequencing (scRNA-seq) datasets, the number of computational tools to analyze aspects of the generated data has grown rapidly. As a result, there is a recurring need to demonstrate whether newly developed methods are truly performant-on their own as well as in comparison to existing tools. Benchmark studies aim to consolidate the space of available methods for a given task and often use simulated data that provide a ground truth for evaluations, thus demanding a high quality standard results credible and transferable to real data.

RESULTS

Here, we evaluated methods for synthetic scRNA-seq data generation in their ability to mimic experimental data. Besides comparing gene- and cell-level quality control summaries in both one- and two-dimensional settings, we further quantified these at the batch- and cluster-level. Secondly, we investigate the effect of simulators on clustering and batch correction method comparisons, and, thirdly, which and to what extent quality control summaries can capture reference-simulation similarity.

CONCLUSIONS

Our results suggest that most simulators are unable to accommodate complex designs without introducing artificial effects, they yield over-optimistic performance of integration and potentially unreliable ranking of clustering methods, and it is generally unknown which summaries are important to ensure effective simulation-based method comparisons.

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.

Clustering-independent estimation of cell abundances in bulk tissues using single-cell RNA-seq data

Single-cell RNA-sequencing has transformed the study of biological tissues by enabling transcriptomic characterizations of their constituent cell states. Computational methods for gene expression deconvolution use this information to infer the cell composition of related tissues profiled at the bulk level. However, current deconvolution methods are restricted to discrete cell types and have limited power to make inferences about continuous cellular processes like cell differentiation or immune cell activation. We present ConDecon, a clustering-independent method for inferring the likelihood for each cell in a single-cell dataset to be present in a bulk tissue. ConDecon represents an improvement in functionality and accuracy with respect to current deconvolution methods. Using ConDecon, we discover the implication of neurodegenerative microglial inflammatory pathways in the mesenchymal transformation of ependymoma, recapitulate spatial patterns of cell differentiation during zebrafish embryogenesis, and make temporal inferences from bulk ATAC-seq data. Overall, ConDecon significantly enhances our understanding of dynamic cellular processes within bulk tissue samples.

Review: cellularity in bone marrow autografts for bone and fracture healing

The use of autografts, as primary cell and tissue source, is the current gold standard approach to treat critical size bone defects and nonunion defects. The unique mixture of the autografts, containing bony compartments and bone marrow (BM), delivers promising results. Although BM mesenchymal stromal cells (BM-MSCs) still represent a major target for various healing approaches in current preclinical research and respective clinical trials, their occurrence in the human BM is typically low. In vitro expansion of this cell type is regulatory challenging as well as time and cost intensive. Compared with marginal percentages of resident BM-MSCs in BM, BM mononuclear cells (BM-MNCs) contained in BM aspirates, concentrates, and bone autografts represent a readily available abundant cell source, applicable within hours during surgical procedures without the need for time-consuming and regulatory challenging cell expansion. This benefit is one reason why autografting has become a clinical standard procedure. However, the exact anatomy and cellularity of BM-MNCs in humans, which is strongly correlated to their unique mode of action and wide application range remains to be elucidated. The aim of this review was to present an overview of the current knowledge on these specific cell types found in human BM, emphasize the contribution of BM-MNCs in bone healing, highlight donor site dependence, and discuss limitations in the current isolation and subsequent characterization procedures. Hereby, the most recent and relevant examples of human BM-MNC cell characterization, flow cytometric analyses, and findings are summarized, with a strong focus on bone therapy.

Biologically informed deep learning to query gene programs in single-cell atlases

The increasing availability of large-scale single-cell atlases has enabled the detailed description of cell states. In parallel, advances in deep learning allow rapid analysis of newly generated query datasets by mapping them into reference atlases. However, existing data transformations learned to map query data are not easily explainable using biologically known concepts such as genes or pathways. Here we propose expiMap, a biologically informed deep-learning architecture that enables single-cell reference mapping. ExpiMap learns to map cells into biologically understandable components representing known 'gene programs'. The activity of each cell for a gene program is learned while simultaneously refining them and learning de novo programs. We show that expiMap compares favourably to existing methods while bringing an additional layer of interpretability to integrative single-cell analysis. Furthermore, we demonstrate its applicability to analyse single-cell perturbation responses in different tissues and species and resolve responses of patients who have coronavirus disease 2019 to different treatments across cell types.