Abstract:Introduction Cancer immunotherapy and targeted therapy have yielded impressive clinical efficacy in acute B lymphoblastic leukemia (B-ALL). Despite the initial high complete remission (CR) rate following first-line therapy, treatment refractoriness and disease relapse remain are correlated with dismal survival. By the time the malignant cells generate, they are accompanied by a rich network of stromal cells and cytokines in bone marrow (BM). This tumor microenvironment (TME) represents an important feature of the biology of B-ALL but also shapes the clinical behavior of the disease. It remains to be confirmed whether the cellular composition and transcriptional heterogeneity impacts the clinical effects of B-ALL. Herein, we analyzed the immune cell infiltration features and related marker genes for B-ALL based on single cell RNA sequencing (scRNA-seq) data, which would be of significance for the development of novel immunotherapies. Methods ScRNA-seq data of 11373 BM cells from 3 B-ALL patients were obtained from the Gene Expression Omnibus (GEO, GSE153358). After quality control and data normalization, cell filtration and marker genes extraction were performed by the Seurat package. Principal component analysis (PCA) and t-distributed stochastic neighbor embedding (t-SNE) were then applied to cluster cells, following with cell types' definition and gene expression profiles in total subsets. Cell clustering was demonstrated using t-SNE-1 and t-SNE-2. In order to determine the cellular characteristics of TME cells mainly mediated by STING pathway, dendritic cell (DC) and B cell were extracted and further plotted gene expression including immunosuppressive molecules and STING pathway, respectively. The pseudo-time analysis was finally performed by Monocle package to display B cell development trajectory and gene expressions over time. Results 9 cell subsets in B-ALL BM, including naïve and memory CD4+ T cell, CD14+ monocyte, B cell, CD8+ T cell, FCGRA3+monocyte, nature killer (NK) cell, DC, and platelet, were identified based on t-SNE analysis (Fig.1A). Top 10 marker genes in each cell cluster were presented in heat map (Fig.1B). Through analyzing differentially expressed genes, we found that BM B cells hardly expressed PD-L1 (CD274), but partially carried TMEM173 (STING), NFKB1 (NF-κB) and GSDMD. In addition, immune cells in BM TME broadly distributed and highly expressed STING and NF-κB, indicating the potential response to type I innate immune response and higher sensitivity to STING agonists than PD-1 antibody (Fig.2A and B). Previous studies had revealed that STING pathway participated in the activation of DC following with production of type I IFNs. We further isolated DC from 9 subsets and profiled the gene expression features. T-SNE analysis revealed 3 subtypes of DC in BM, marker genes comparison further identified as monocyte derived DC, CD1C-CD14-DC and myeloid conventional DC. Cyclic GMP-AMP synthase (cGAS), STING and NF-κB were highly expressed in each type in compared with PDCD1 (PD-1) and highest existed on myeloid conventional DC (Fig.3). We then explored B cell subsets to determine whether STING pathway could induce cell pyroptosis in BM B cells. Different subgroup of B cells shared similar marker genes, companying with higher expression of NF-κB and GDSMD (Fig.4). Furthermore, pseudo-time analysis plotted the development trajectory of malignant B cells. The results showed that GSDMD gradually increased along with cell development, suggesting that STING agonist would be sensitive to mature B cells (Fig.5). Subsets analysis shown that anti-tumor immune response of DC and pyroptosis of B cell might be triggered through STING pathway activation. Conclusion Our study profiled for the first time the expression of STING pathway in BM DC and B cell from B-ALL patients based on single-cell transcriptome. Combination of STING agonist and conventional immunotherapy had been shown prospects in antitumor therapy. STING agonist is expected to be an adjuvant drug for B-ALL immunotherapies in the future. Keywords: Single-cell RNA sequencing; acute B lymphoblastic leukemia; STING; PD-1; immunotherapy. Disclosures No relevant conflicts of interest to declare.

Multi-Faceted Effects of ST6Gal1 Expression on Precursor B-Lineage Acute Lymphoblastic Leukemia

ST6GALNAC4 Promotes Hepatocellular Carcinogenesis by Inducing Abnormal Glycosylation.

Differential Expression of ST6GAL1 in the Tumor Progression of Colorectal Cancer.

[Effect of B4GALT1 on Proliferation of Its Co-cultured Human Acute Myeloid Leukemia Cell Line in Bone Marrow Stromal Cells].

Single-Cell Transcriptome Reveals the Role of Sting Pathway in Acute B Lymphoblastic Leukemia Bone Marrow Cells

The Sialyltransferase ST3GAL6 Influences Homing and Survival in Multiple Myeloma.

Remodeling of marrow hematopoietic stem and progenitor cells by non-self ST6Gal-1 sialyltransferase

Evaluation of ST6Gal1 expression and clinicopathological significance in human glioma

Hypersialylation of Tumor Cells Promotes Pancreatic Cancer Progression

ST6GAL1 inhibits metastasis of hepatocellular carcinoma via modulating sialylation of MCAM on cell surface

TSLP-Induced Alterations of Multiple Signaling Pathways in Primary CRLF2 B-ALL Xenografts

ST6GAL2 Downregulation Inhibits Cell Adhesion and Invasion and is Associated with Improved Patient Survival in Breast Cancer.

The β-galactoside α2,6-sialyltranferase 1 (ST6GAL1) inhibits the colorectal cancer metastasis by stabilizing intercellular adhesion molecule-1 via sialylation.

Remodelling of the glycome of B-cell precursor acute lymphoblastic leukemia cells developing drug-tolerance

Single-cell RNA-seq analysis revealed the stemness of a specific cluster of B cells in acute lymphoblastic leukemia progression

ST6GAL1 sialyltransferase promotes acinar to ductal metaplasia and pancreatic cancer progression

Identification of ST3GAL6-AS1 Interactive Protein in Multiple Myeloma

The pro-inflammatory cytokines IL-1β and IL-6 promote upregulation of the ST6GAL1 sialyltransferase in pancreatic cancer cells

Potential molecular mechanisms of ETV6-RUNX1-positive B progenitor cell cluster in acute lymphoblastic leukemia revealed by single-cell RNA sequencing

ASSESSING THE IMPACT OF A B-ALL-ASSOCIATED GERMLINE ETV6 VARIANT ON MURINE HEMATOPOIESIS AND STEM CELL FUNCTION