Neutrophil-related genes predict prognosis and contribute to immunosuppression in acute myeloid leukemia

doi:10.21203/rs.3.rs-4853209/v1

Download PDF

Article

Neutrophil-related genes predict prognosis and contribute to immunosuppression in acute myeloid leukemia

https://doi.org/10.21203/rs.3.rs-4853209/v1

This work is licensed under a CC BY 4.0 License

You are reading this latest preprint version

Neutrophils, the predominant myeloid cells in human blood, have been implicated as pivotal contributors to cancer development. However, the relationship between neutrophils and acute myeloid leukemia (AML) prognosis remains unclear due to their great plasticity and diversity. In this study, we conducted a comprehensive analysis of gene expression data and clinical records from AML patients sourced from TCGA, GEO, and OHSU databases. We established a neutrophil-based prognostic model incorporating five neutrophil-related genes (CSF3R、BRAF、FFAR2、CD300A and CD37) and the predictive value of the model was validated in both internal and external validation cohorts. Univariate and Multivariate Cox regression analyses further demonstrated that the model remained an independent prognostic factor for overall survival, and a nomogram was constructed to facilitate its clinical application. Notably, CD37 was identified as a crucial neutrophil-related gene associated with adverse AML prognosis. Elevated CD37 expression was indicative of dysregulated cell proliferation, immunosuppressive tumor microenvironment and T cell dysfunction. Overexpression of CD37 was also correlated with increased expression of PD1, CTLA4, CD86, and LAG3. Furthermore, CD37 could serve as a valuable predictor for immunotherapy and represents a druggable target in AML. In summary, our study advances a novel prognostic model based on neutrophil-related gene expression, with CD37 emerging as a promising biomarker for anti-leukemic response prediction. This insight holds promise for personalized prognosis estimation and informed clinical decision-making in AML management.

Biological sciences/Cancer

Biological sciences/Immunology

Health sciences/Biomarkers

Health sciences/Diseases

Health sciences/Oncology

Health sciences/Risk factors

neutrophil

CD37

AML

immunotherapy

tumor microenvironment

Acute myeloid leukemia (AML) is a highly aggressive cancer characterized by a block in myeloid differentiation that leads to the uncontrolled proliferation of myeloblasts^[1]. It stands as the most prevalent leukemia in adults, with an age-adjusted incidence rate of 4.3 cases per 100,000 individuals annually^[2]. Despite recent advances in targeted therapies, AML remains a formidable disease with a dismal prognosis^[3]. Therefore, there is an urgent need for novel therapeutic targets, and one of the promising approaches is cancer immunotherapy. Indeed, targeting an immune checkpoint regulator in combination with chemotherapeutics has been proven to be more effective than the chemotherapeutic regimen alone^{[4, 5]}. With immense success in the treatment of various solid tumors and lymphomas, immune checkpoint inhibitors (ICIs) are now considered as a potential immunotherapeutic option for AML patients^[6]. Encouraging results from clinical trials using ICIs have demonstrated improved response rates and overall survival (OS) in patients with relapsed or refractory AML^[7]. Thus, effective indicators and prognostic biomarkers must be explored to identify AML patients who will benefit from ICIs therapies, ultimately facilitating risk stratification and personalized treatment approaches.

Neutrophils are key mediators of chronic inflammation and innate inflammatory response, as well as important regulators of cancer development^[8]. Neutrophils are abundant and heterogeneous leukocytes in human peripheral blood functioning as the early responder to infection or injury^[9]. Unlike terminally differentiated cell populations, tumor-associated neutrophils are susceptible to environmental factors and exhibit great plasticity and diversity, thereby promoting or interfering with tumor growth^[10]. While neutrophils can potentially inhibit tumor progression by generating anti-tumor factors, prevailing evidence suggests they more commonly act as facilitators of cancer progression and metastasis, by orchestrating various processes such as tumor survival, migration, immune evasion, and angiogenesis^[11–15]. Notably, neutrophils have been implicated in influencing cancer prognosis and responses to diverse anticancer treatments, from chemotherapy to immunotherapy. In numerous malignancies, tumor-associated neutrophil infiltration, peripheral blood neutrophil to lymphocyte ratios (NLR), and neutrophil-based transcriptional profiles are correlated with unfavorable clinical outcomes and reduced drug sensitivity^[16–18]. However, the prognostic role of neutrophils in AML remains elusive, nor is the predictive value of neutrophils in specific therapeutic regimens.

In this study, we conducted a comprehensive analysis of gene expression profiles within the TCGA, GEO, and OHSU databases to assess the influence of neutrophils on AML prognosis. A prognostic model based on neutrophil-associated genes was developed and validated. The relationship between neutrophil-related genes and immune dysregulation, as well as drug sensitivity was further investigated to uncover potential indicators of immunotherapy response and identify novel therapeutic targets for AML treatment.

2.1 Data source and Processing

The Cancer Genome Atlas (TCGA) is a comprehensive database sponsored by the government of the United States, collecting more than 11,000 cases across 33 tumor types. The Gene Expression Omnibus (GEO) is the largest public resource for gene expression data, storing chips, second-generation sequencing, and other high-throughput sequencing data. Beat AML 1.0 was developed from an agreement 6 years ago between Oregon Health & Science University (OHSU) and the Leukemia and Lymphoma Society to develop a consortium of 11 academic medical center partners. OHSU data comes from the Beat AML 1.0 program. We downloaded clinical information and transcriptome data of 151 AML patients from TCGA for our analysis. Overall, 151 patients were randomly divided into two data sets in a 7:3 ratio by the R package “Caret”: training set (n = 107) and internal validation set (n = 44). AML patients from GEO datasets, including GSE10358 (n = 178), GSE14468 (n = 262), GSE37642-GPL570 (n = 136), GSE71014 (n = 104, CN-AML), and OHSU (n = 576) was downloaded as an external validation set.

2.2 Source of neutrophil-related gene

Neutrophil-related genes were retrieved from published studies, which can be classified into three gene sets: genes involved in neutrophil function, specificity, and infiltration, respectively. For genes involved in neutrophil function, 55 genes were extracted from the observation by Rincón et al^[19]. They were vital for the function of cell surface receptors, reactive oxygen species, guanine nucleotide exchange factors, and adhesion receptors in neutrophils. For genes involved in neutrophil specificity, 31 genes specifically expressed in neutrophils were obtained from research defining the immune landscape of human colon cancer^[20]. Finally, 84 genes associated with neutrophil infiltration were obtained from the LM22, a leukocyte gene signature matrix applied in the CIBERSORT algorithm, which was widely used for the estimation of immune cell abundance from tissue expression profiles^[21]. After removing the duplicated genes in the three gene sets, a total of 148 genes were identified as neutrophil-related genes for subsequent analyses.

2.3 Construction and validation of the risk model

To select prognostic genes associated with neutrophils, univariate COX regression was performed on AML patients using the R package "survival", and these genes could be used as candidates for the construction of prognostic models. The forest map of prognostic neutrophil-related genes was constructed using the R package “forestplot”. The prognostic genes were then identified by Least Absolute Shrinkage and Selection Operator (LASSO) regression analysis via the R packages “glmnet” and a 5-gene risk model was established. The function “predict” in the R package “glmnet” was used to calculate the risk score of individual patients. Patients were divided into high-risk and low-risk groups based on median risk scores. The prognostic value of the 5-gene risk models was verified by the R packages "survival" and "survminer" using Kaplan-Meier curves and the log-ranking method, and the sensitivity of the model was assessed by receiver operating characteristic (ROC) curve analysis of the R packages "ROCR" and "rms".

2.4 Nomogram construction

To determine the survival probabilities predicted by the risk model and to facilitate its clinical application, a nomogram was drawn using the R package "rms".

2.5 Survival analysis

Kapan-Meier plots were drawn to elucidate the effect of each prognostic gene on the overall survival (OS) of AML patients. The Log-rank test was used to evaluate their relationships.

2.6 Analysis of differentially expressed genes

The R package “limma” was used to select differentially expressed genes (DEGs). DEGs meeting the following criteria were considered significant: adjusted p-values < 0.05 and | log2 fold change (FC)|\(\:\ge\:\)1.

2.7 Pathway enrichment analysis

Gene Ontology (GO) and Kyoto Gene Genome Encyclopedia (KEGG) pathway analyses were performed via the R package “clusterprofiler”, with a cutoff value of FDR < 0.05. To study the differences in biological processes between the high and low expression groups of CD37, we applied gene set variation analysis (GSVA) with the R package “GSVA”.

2.8 Immune cell infiltration

To assess the composition of immune cells in the tumor microenvironment, we used the ESTIMATE algorithm to evaluate the immune score, tumor purity, and stromal score of tumor samples. CIBERSORT is an analytical tool from the Alizadeh Lab and Newman Lab to impute gene expression profiles and provide an estimation of the abundances of member cell types in a mixed cell population, using gene expression data. xCell is a recently published method based on ssGSEA that estimates the abundance scores of 64 immune cell types, including adaptive and innate immune cells, hematopoietic progenitors, epithelial cells, and extracellular matrix cells. CIBERSORT and xCell algorithm was applied to assess the proportion of immune cell infiltration. The relative abundance of 24 tumor-infiltrating immune cells was obtained by single sample gene set enrichment analysis (ssGSEA) using the R package "GSVA". We then obtained gene sets of CD4 + T cells, naive CD8 + T cells, and naive CD4 + T cells from the Molecular Signatures Database (MSigDB) and calculated their enrichment using the GSEA algorithm.

2.9 Statistical analysis

The T test was used to determine the difference between two groups. The "survminer" package divided patients into the high- and low-risk score groups or high- and low-CD37 groups based on median. The log-rank test was used to determine the P values between the groups in the Kaplan–Meier survival analysis. The univariate and multivariate Cox regression analyses were used to identify the prognostic factors. A ROC curve analysis was used to determine the specificity and sensitivity of the related metrics. R version 4.2.3 software was used for statistical processing. Visualization of data was performed with R software and GraphPad Prism V.9.4.1. P < 0.05 was considered significant.

3.1 Establishment of the neutrophil-based prognostic model

To evaluate the impact of neutrophils on the clinical outcome of AML patients, we performed a univariate Cox regression analysis using the TCGA dataset. The results revealed eight neutrophil-related genes associated with prognosis, including six genes detrimental to OS (LST1, ITGAX, FFAR2, CD300A, ITGAL, CD37), and two genes beneficial to OS (CSF3R, BRAF) (Fig. 1A). Subsequently, we employed LASSO regression analysis on these prognostic genes to construct a neutrophil-based prognostic model. By determining the active coefficients, the model achieved optimal performance with a specified \(\:\lambda\:\) value of 0.0313 and five key prognostic genes (CSF3R, BRAF, FFAR2, CD300A, and CD37) were incorporated (Figs. 1B-D). The individual-level risk score determined by this prognostic model was utilized for risk stratification and prediction of prognosis.

3.2 Validation of the neutrophil-based prognostic model

The neutrophil-based prognostic model was used to calculate the risk score for each patient in the training cohort form TCGA dataset. According to the median risk scores of 0.97 in the training dataset, AML patients were divided into high-risk and low-risk groups. Survival analysis showed that patients in the high-risk group exhibited significantly shorter OS than patients in the low-risk group (P < 0.0001, Fig. 2A). Then a ROC curve was drawn to assess the predictive performance of the prognostic model. As shown in Fig. 2B, the area under the curve (AUC) of 1, 3, and 5 years was 0.750, 0.726, and 0.821, respectively. A similar trend was displayed in the internal validation dataset that patients with high-risk scores demonstrated inferior clinical outcomes (P = 0.044, Fig. 2C), with the AUC values of 0.733, 0.739, and 0.637 at 1, 3, and 5 years (Figs. 2D).

Furthermore, the analysis of external validation datasets from the GSE10358 cohort (n = 178, Figure S1A-B), the GSE14468 cohort (n = 262, Figure S1C-D), the GSE71014 cohort (n = 104, Figure S1E-F), the GSE37642-GPL570 (n = 136, Figure S1G-H), corroborated these results, implying an adverse prognostic impact of the neutrophil-based model. In the merged GEO datasets (n = 680), OS of the high-risk group was significantly shorter than that of the low-risk group (P = 0.02) and the AUC of 1, 3, and 5 years was 0.601, 0.646 and 0.637 (Figs. 2E-F). Additionally, we confirmed our findings using the OHSU dataset (n = 576), where the high-risk group illustrated significantly lower OS compared to the low-risk group (P = 0.0023) and the AUC of 1, 3, and 5 years was 0.593, 0.554 and 0.523 (Figs. 2G-H). Together, through rigorous validation in both internal and external cohorts, our developed risk-score model demonstrated reliable and stable prognostic predictive value.

3.3 Construction of predictive nomogram

To assess whether the neutrophil-based prognostic model impacted OS independent of recognized prognostic indicators for AML, we carried out univariate and multivariate Cox regression analyses. In multivariate models for the TCGA cohort, including the variables with P value less than 0.05 under univariate analysis (Fig. 3A), the risk score remained an independent prognostic indicator for OS (P = 0.015), along with age (P < 0.001), WBC (P = 0.040), cytogenetics risk (P = 0.023), FLT3 mutation (P = 0.003), RUNX1 mutation (P = 0.006) and TP53 mutation (P < 0.001) (Fig. 3B). Subsequently, we developed a nomogram to predict 1-, 3- and 5-year OS in AML patients for the convenience in clinical utility. Based on the results of the multivariate Cox analysis, the risk score, age, WBC, cytogenetics risk, FLT3 mutation, RUNX1 and TP53 mutation were included in the nomogram. The total point scores for all parameters were matched with the corresponding survival time scales (Fig. 3C). Collectively, we validated the remarkable predictive capability of our prognostic signatures and underscored its substantial potential for clinical applicability from multiple perspectives.

3.4 Correlation of neutrophil-related genes with prognosis in AML

With the aim of exploring the function of neutrophil-related genes incorporated into the prognostic model, survival analysis and ROC curve were performed (Fig. 4A-B, Figure S2). The Kaplan-Meier curves indicated that four of these genes (CSF3R, BRAF, CD300A, CD37) were involved in the prognosis of AML, while FFAR2 was not significant. Among these genes, CD37 emerged as the strongest predictor of survival (p = 0.0026), with the AUC values of 0.698, 0.626, and 0.680 at 1, 3, and 5 years (Fig. 4A). The prognostic value of CD37 was further validated using GEO datasets from the GSE10358 cohort (p = 0.00094, Fig. 4B), the GSE14468 cohort (p = 0.041, Fig. 4C) and the GSE71014 cohort (p = 0.0087, Fig. 4D), suggesting that CD37 is the critical neutrophil-related signatures linked to an adverse prognosis of AML.

Several gene expression-based prognostic models for AML have been recently emerged, with the LI24 and LSC17 models showing superior prognostic capabilities and the ability to improve risk stratification^{[22, 23]}. Therefore, it is interesting to explore the predictive power of CD37 expression within the context of these novel schemes. We implemented both models in TCGA database and patients were stratified into high- and low-risk groups as shown in Figs. 4E and F. When applied to each risk group stratified by LI24, CD37 expression remained a significant predictor for OS in the high-risk group (Fig. 4E). Similarly, within the high-risk categories defined by the LSC17 model, CD37 status was able to distinguish survival outcomes (Fig. 4F). In summary, these results indicate that CD37 could serve as a valuable candidate for refining existing classification schemes, potentially improving patient management and treatment stratification in AML.

3.5 Association of CD37 expression with immunity and cellular structure

To gain insight into the involvement of CD37 in AML development, patients in the TCGA dataset were dichotomized at the median expression of CD37. The principal component analysis showed that the gene expression profiles were significantly different between the CD37 high and CD37 low groups (Fig. 5A). The “limma” package was then used in the R environment to identify DEGs between the two groups. A total of 902 DEGs, including 584 up-regulated genes and 318 down-regulated genes, were selected with a cutoff P-value of 0.05 and |fold- change| > 1 (Fig. 5B).

Functional enrichment analysis was performed using the GO and KEGG databases to elucidate the biological functions of DEGs. GO enrichment analysis revealed that upregulated DEGs were primarily associated with immune reaction, including leukocyte mediated immunity, regulation of immune effector process, negative regulation of immune system process, negative regulation of T cell proliferation and regulation of leukocyte mediated immunity (Fig. 5C). Meanwhile, downregulated DEGs were enriched in cellular structure, such as synapse assembly, extracellular matrix organization, extracellular structure organization, and cell junction assembly (Fig. 5D). KEGG enrichment analysis demonstrated that upregulated DEGs were predominantly enriched in antigen processing and presentation, as well as the B cell receptor signaling pathway (Fig. 5E), while downregulated DEGs were mainly associated with the cAMP signaling pathway and GABAergic synapse (Fig. 5F). Furthermore, GSEA analysis yielded similar conclusions to the GO and KEGG analysis. High CD37 expression was closely correlated with immune reaction, including adaptive immune response, leukocyte mediated immunity, lymphocyte mediated immunity, T cell activation, immune receptor activity and antigen processing and presentation (Fig. 5G, I, J). Additionally, high CD37 expression was associated with gene signatures related to cellular structures, such as cell surface, external side of the plasma membrane, membrane side, and vacuole (Fig. 5H). These results allow us to reasonably infer that CD37 is involved in immune regulation of AML.

3.6 CD37-shaped immune disorders context

To further elucidate the role of CD37 in the immune microenvironment of AML, we performed ESTIMATE analysis to assess the level of immune infiltration. The results showed that the CD37 high group had higher immune scores and overall ESTIMATE scores compared to the CD37 low group (Fig. 6A). Consistently, xCell analysis demonstrated that the higher CD37 expression levels, the more immunescore fractions (Fig. 6B). Similar results were observed when the immune fractions were calculated by ESTIMATE and xCell analysis in another two distinct cohorts (GSE10358 and GSE14468, Figure S3A-D). Since cancers with higher human leukocyte antigen (HLA) gene expression were often more immunologically active ^[24], we thereafter investigated the correlation of CD37 and the HLA markers. Indeed, CD37 expression was positively associated with the majority of these HLA genes (Fig. 6C). These findings provided insight into the increased immune cell infiltration of high CD37 expressers, which warrants further exploration.

ssGSEA analysis was then applied to calculate the proportion of each type of immune cell within the tumor microenvironment (TME). As shown in Fig. 6D, neutrophils were significantly enriched in AML patients with a CD37 high score compared to those with a CD37 low score. Furthermore, the preferential enrichment of neutrophils in the CD37 high group was also found in GSE10358 and GSE14468 datasets (Figure S4A-B). In the CD37 high group, myeloid-derived suppressor cells (MDSCs), NK cells, NKT cells and Th1 cells were increased (Fig. 6D), whereas Th2 cells and the ratio of CD4⁺ T /CD8⁺ T cells were decreased (Fig. 6E), indicating the presence of an immunosuppressive environment. This was also true for the higher fractions of immunosuppressive components like MDSC, NK cells and NKT cells in the CD37 high group within two independent AML cohorts (Figure S4C-D). By CIBERSORT analysis, we also observed a significant enrichment of Treg cells and activated NK cells as well as a concomitant reduction in the ratio of CD4⁺ T /CD8⁺ T cells in the CD37 high group, as compared to the CD37 low group (Figure S5A-B). These data suggest that CD37 plays a significant role in shaping the highly disturbed immune microenvironment of AML.

Considering the crucial roles of immune checkpoints in the TME, it is therefore of great interest to evaluate the relationship between CD37 and these immune signatures. Our findings revealed robust positive correlations between CD37 expression and immune checkpoint molecules such as LAG3, LGALS9, CD40, and TNFRPSF8 (Fig. 6F). Moreover, high CD37 expressers exhibited significant increases in those key targets relevant to AML, including PD1, CTLA4, CD86, and LAG3 (Fig. 6G).

3.7 Predictive value of CD37 in responsiveness to immunotherapies

Given the immunophenotypic discrepancy between the CD37 high and CD37 low group, it is reasonable to speculate that CD37 holds promise in predicting the immunotherapy response. Subsequent unsupervised group analysis revealed significant distinctions in the immune profiles of these two groups (Fig. 7A). Further correlation analysis indicated strong associations of CD37 expression with immune cell infiltration, particularly highlighting central memory CD4⁺ T cell, NK cell and MDSCs (Fig. 7B). Additionally, we explored the relationship between CD37 expression and the functional status of immune cells. GSEA analysis showed that CD4⁺ T cells, naïve T cells, naïve CD8⁺ T cells, naïve B cells, naïve CD4⁺ T cells, monocytes and neutrophils gene sets were significantly enriched in the CD37 high group (Fig. 7C). This suggests that CD37 may play a role in coordinating the aberrant function of immune cells.

Recently, a gene expression-based coring system called Tumor Immune Dysfunction and Exclusion (TIDE) has demonstrated outstanding performance in predicting clinical responses to immunotherapy ^[25]. We therefore assessed the relationship between CD37 and the expression signatures associated with T-cell dysfunction and T-cell exclusion. Surprisingly, we discovered a negative correlation of CD37 with T cell exclusion signatures, including MDSCs, M2 subtype of tumor-associated macrophages (TAMs), Exclusion and the TIDE score, whereas a positive correlation with the T-cell dysfunction score, IFNG, and Merck18 signatures (Fig. 7D). This indicates that CD37 molecules might contribute to immune evasion through the induction of T-cell dysfunction. Based on these results, we further explored whether CD37 could serve as a predictive marker for immunotherapy response in AML. In the TCGA AML cohort, we observed significantly higher expression levels of CD37 in predicted responders (Fig. 7E), suggesting CD37 expression could indeed be a valuable predictor for immunotherapy in AML.

Additionally, drug sensitivity was predicted utilizing mRNA expression profiles and anti-tumor drug activity data of CD37 in the GDSC database through the oncoPredict algorithm. As illustrated in Fig. 7F, a total of 146 drugs were found to interact with CD37, indicating its potential as a druggable target. While CD37 is closely associated with immune regulation, small molecular inhibitors such as BCL2 inhibitors could potentially target CD37. These findings implied that treatment targeting a crucial neutrophil-related gene CD37 might improve the TME of AML patients.

The current risk stratification of AML patients relies on chromosomal variants and genetic mutations, categorizing them into high risk, intermediate risk, and low risk. However, these conventional prediction tools have mainly focused on the intrinsic characteristics of cancer, ignoring the contribution of TME to cancer development^[26]. Recently, the involvement of tumor-associated neutrophils, a key component of the TME, has become increasingly prominent in cancer progression. Tumor-infiltrating neutrophils exhibit pro-carcinogenic effects and predict poor prognosis in bladder urothelial cancer^[27]. Furthermore, peripheral NLR has emerged as a potential prognostic factor in various solid tumors^{[28] [29] [30] [31]}, as well as in certain hematological malignancies such as non-Hodgkin lymphoma and multiple myeloma^{[32, 33]}. Additionally, neutrophils have been implicated in influencing the efficacy of immunotherapy and mediating resistance to radiation therapy^[34]. In this study, we investigated the effect of neutrophil-related genes on AML prognosis and explored a potential prediction system. Importantly, CD37 was identified as a critical neutrophil-related gene associated with adverse AML prognosis.

Advancements in next-generation sequencing (NGS) technology have revolutionized our ability to access the genomic profiles of cancer, aiding in precise cancer classification and prognosis prediction. In this study, we focused on investigating the prognostic role of neutrophil-related genes in AML. Our analysis revealed that eight neutrophil-related genes—CSF3R, LST1, ITGAX, BRAF, FFAR2, CD300A, ITGAL, and CD37—exhibited strong correlations with the prognosis of AML. By integrating five of these genes, we constructed an effective prognostic model. Notably, a significant prevalence of acquired CSF3R mutations has been observed in patients with severe congenital neutropenia during the pre-leukemia stage, progressing to overt AML or myelodysplastic syndrome (MDS)^[35]. Intriguingly, there is a high incidence of transformation to MDS or AML in patients who harbor acquired CSF3R mutations, suggesting the importance of CSF3R in the pathogenesis and prognosis of AML^[36]. BRAF gene, belonging to the RAF kinase family, plays pivotal roles in transmitting growth signals in physiological pathways^{[37, 38]}. As a critical constituent of the MAPK pathway, BRAF has been detected as driver mutation in several tumor types, including melanoma, non-small cell lung cancer (NSCLC), and anaplastic thyroid cancer (ATC). Consequently, BRAF has emerged as an appealing target for therapeutic inhibition^[39].

Activation of FFAR2 by microbiota-derived metabolites has been shown to decrease the proliferation of leukemic cells in vitro^[40]. Moreover, FFAR2 itself impacts leukemia cell growth in vivo^[41]. Deletion of FFAR2 compromises the immunosuppressive capabilities of MDSCs on T cells within the tumor microenvironment^[42]. CD300A, also recognized as CMRF-35 or IRp60, belongs to the CD300 cell surface molecule superfamily^[43], which plays a critical role in modulating immune function and contributes to the host response against various diseases, including infectious diseases^{[44, 45]}, cancer^[46], and allergy^{[47, 48]}. The expression levels of CD300A were associated with risk stratification and the clinical relevance of AML. Elevated CD300A expression may serve as an independent adverse prognostic indicator for OS and relapse-free survival (RFS) in AML^[49]. Furthermore, CD300A is involved in regulating neutrophil recruitment, IL-1β production, and participates in the processes of neutrophil apoptosis and efferocytosis^[50].

CD37, belonging to the transmembrane 4 superfamily (TM4SF), is a tetraubiquitin protein prominently expressed on B cell surfaces. CD37 promotes neutrophil adhesion and recruitment by enhancing cytoskeletal function downstream of integrin-mediated adhesion processes^[51]. This characteristic makes CD37 an attractive molecular target for immunotherapy against B-cell lymphomas and leukemias^[52]. AML patients with high CD37 expression were shown to have shorter OS and disease-free survival (DFS)^[53]. While previous studies primarily analyzed CD37 survival data from TCGA datasets, our study extended this analysis by incorporating three additional independent AML datasets and two AML prognostic models. Our investigation provided a comprehensive evaluation of the prognostic relevance of CD37 in AML, confirming its association with poorer outcomes. The consistent findings across multiple datasets and prognostic models underscore CD37's potential as a promising biomarker for assessing cancer risk in AML patients. However, the precise mechanism of CD37 in AML remains to be fully elucidated.

CD37 plays a pivotal role in cancer immunity. Our study unveiled a significant correlation between CD37 expression and gene signatures associated with immune regulation in AML. These signatures encompassed various aspects such as leukocyte-mediated immunity, regulation of immune effector processes, negative regulation of the immune system and inhibition of T cell proliferation. Furthermore, the multifaceted involvement of CD37 protein in diverse biological processes including cell adhesion, motility, differentiation, proliferation, metastasis, growth, survival, trafficking, intercellular communication via exosomes, and immune responses is widely acknowledged^[54]. Additionally, CD37 interacts with key proteins such as integrins, immune receptors, and signaling molecules, which serve as regulators for leukocyte activation, motility, and antigen presentation^[55].

Emerging studies suggest the pivotal role of CD37 in regulating B-cell survival and shaping immune evasion. Specifically, tumors with high expression of CD37 demonstrate increased infiltration by various immune cell subsets, compared to those with lower CD37 expression levels. It has been known that AML patients commonly present with dysfunctional T cells and NK cells at diagnosis^{[56, 57]}. In our investigation, we observed an augmented presence of immunosuppressive cells such as MDSC, Treg cells and NK cells in the high CD37 expression group. Moreover, the elevated ratio of CD4⁺ T cells to CD8⁺ T cells suggests an immunosuppressive leukemic microenvironment. This milieu facilitates AML blast evasion from immune surveillance and collaborates to promote disease progression, which might explain the poorer outcomes observed in individuals with high CD37 expression. Furthermore, CD37 may play a role in orchestrating the dysfunction of antitumor immune cells. Our analysis also reveals enrichment of gene sets associated with CD4⁺ T cells, naive T cells, naive CD8⁺ T cells, naive B cells and naive CD4⁺ T cells in the high CD37 expression group. However, further investigations are warranted to elucidate the precise mechanisms by which CD37 influences the development and activation of immune cells.

ICIs therapies, which rejuvenate the effective antitumor immune response mediated by T cells, have revolutionized cancer treatment. Our research indicates a positive association of CD37 expression with T-cell dysfunction scores in AML. High CD37 expression predicts a more favorable response to immunotherapy. Additionally, CD37 correlates positively with inhibitory immune checkpoints like PD1, CTLA4, CD86, and LAG3, which are often overexpressed in AML cells, presenting promising therapeutic targets. Furthermore, our study reveals a positive correlation between CD37 and MHC class II molecules. The complexes of tetraspanins MHC class II molecules exhibit enrichment with the CD86 co-stimulatory molecule and the HLA-DM peptide editor^[58]. We propose that CD37 regulates peptide presentation on MHC molecules, potentially influencing T-cell activation. Collectively, these findings suggest CD37 as a potential biomarker for AML immunotherapy.

In summary, we developed a prognostic model for AML incorporating five neutrophil-related genes, demonstrating its potential independent prognostic significance. Moreover, we identified a robust link between CD37 expression and immune response in AML. However, the validation of our findings is constrained by insufficient data. Increased clinical samples and prospective studies are necessary to confirm the predictive efficacy of the neutrophil-based model in real-world clinical scenarios. The underlying mechanism of neutrophil-related genes within tumor microenvironment remains ambiguous and requires further exploration in future studies.

Conflicts of Interests

The authors declare no conflict of interest.

Ethical Statement

All procedures performed in studies involving human participants were in accordance with the ethical standards of the institutional and/or national research committee and with the 1964 Helsinki Declaration and its later amendments or comparable ethical standards.

Human Ethics and Consent to Participate declarations: not applicable.

Conflict of Interest:

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Funding

This work was supported by grants from the Distinguished Taishan Scholars in Climbing Plan (tspd20210321), the Distinguished Taishan Scholars Plan (NO. tstp20230653), the National Natural Science Foundation of China (82070160, 82370165), the Fundamental Research Funds for the Central Universities (2022JC012), the Independently Cultivate Innovative Teams of Jinan, Shandong Province (2021GXRC050).

Author Contribution

Conceptualization, DZ and JZ; methodology, DZ; software, DZ; validation, DZ; formal analysis, DZ; YL; TL; XL; investigation, DZ.; resources, DZ; data curation, DZ.; writing—original draft preparation, DZ.; writing—review and editing, JZ.; visualization, DZ.; supervision, JZ.; project administration, JZ.; funding acquisition, JZ. All authors have read and agreed to the published version of the manuscript.

Data Availability

The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.

JORDAN C T. Unique molecular and cellular features of acute myelogenous leukemia stem cells [J]. Leukemia, 2002, 16(4): 559-62.
SHALLIS R M, WANG R, DAVIDOFF A, et al. Epidemiology of acute myeloid leukemia: Recent progress and enduring challenges [J]. Blood Rev, 2019, 36: 70-87.
SHAH A, ANDERSSON T M, RACHET B, et al. Survival and cure of acute myeloid leukaemia in England, 1971-2006: a population-based study [J]. Br J Haematol, 2013, 162(4): 509-16.
YANG H, ROSOVE M H, FIGLIN R A. Tumor lysis syndrome occurring after the administration of rituximab in lymphoproliferative disorders: high-grade non-Hodgkin's lymphoma and chronic lymphocytic leukemia [J]. Am J Hematol, 1999, 62(4): 247-50.
HEINHUIS K M, ROS W, KOK M, et al. Enhancing antitumor response by combining immune checkpoint inhibitors with chemotherapy in solid tumors [J]. Ann Oncol, 2019, 30(2): 219-35.
!!! INVALID CITATION !!! [6].
DAVER N, GARCIA-MANERO G, BASU S, et al. Efficacy, Safety, and Biomarkers of Response to Azacitidine and Nivolumab in Relapsed/Refractory Acute Myeloid Leukemia: A Nonrandomized, Open-Label, Phase II Study [J]. Cancer Discov, 2019, 9(3): 370-83.
HEDRICK C C, MALANCHI I. Neutrophils in cancer: heterogeneous and multifaceted [J]. Nat Rev Immunol, 2022, 22(3): 173-87.
XIAO Y, CONG M, LI J, et al. Cathepsin C promotes breast cancer lung metastasis by modulating neutrophil infiltration and neutrophil extracellular trap formation [J]. Cancer Cell, 2021, 39(3): 423-37.e7.
JAILLON S, PONZETTA A, DI MITRI D, et al. Neutrophil diversity and plasticity in tumour progression and therapy [J]. Nat Rev Cancer, 2020, 20(9): 485-503.
ACHARYYA S, OSKARSSON T, VANHARANTA S, et al. A CXCL1 paracrine network links cancer chemoresistance and metastasis [J]. Cell, 2012, 150(1): 165-78.
COFFELT S B, KERSTEN K, DOORNEBAL C W, et al. IL-17-producing γδ T cells and neutrophils conspire to promote breast cancer metastasis [J]. Nature, 2015, 522(7556): 345-8.
QIAN B Z, LI J, ZHANG H, et al. CCL2 recruits inflammatory monocytes to facilitate breast-tumour metastasis [J]. Nature, 2011, 475(7355): 222-5.
WCULEK S K, MALANCHI I. Neutrophils support lung colonization of metastasis-initiating breast cancer cells [J]. Nature, 2015, 528(7582): 413-7.
ZHUANG X, ZHANG H, LI X, et al. Differential effects on lung and bone metastasis of breast cancer by Wnt signalling inhibitor DKK1 [J]. Nat Cell Biol, 2017, 19(10): 1274-85.
TEMPLETON A J, MCNAMARA M G, ŠERUGA B, et al. Prognostic role of neutrophil-to-lymphocyte ratio in solid tumors: a systematic review and meta-analysis [J]. J Natl Cancer Inst, 2014, 106(6): dju124.
GENTLES A J, NEWMAN A M, LIU C L, et al. The prognostic landscape of genes and infiltrating immune cells across human cancers [J]. Nat Med, 2015, 21(8): 938-45.
SHAUL M E, FRIDLENDER Z G. Tumour-associated neutrophils in patients with cancer [J]. Nat Rev Clin Oncol, 2019, 16(10): 601-20.
RINCÓN E, ROCHA-GREGG B L, COLLINS S R. A map of gene expression in neutrophil-like cell lines [J]. BMC Genomics, 2018, 19(1): 573.
BINDEA G, MLECNIK B, TOSOLINI M, et al. Spatiotemporal dynamics of intratumoral immune cells reveal the immune landscape in human cancer [J]. Immunity, 2013, 39(4): 782-95.
NEWMAN A M, LIU C L, GREEN M R, et al. Robust enumeration of cell subsets from tissue expression profiles [J]. Nat Methods, 2015, 12(5): 453-7.
LI Z, HEROLD T, HE C, et al. Identification of a 24-gene prognostic signature that improves the European LeukemiaNet risk classification of acute myeloid leukemia: an international collaborative study [J]. J Clin Oncol, 2013, 31(9): 1172-81.
NG S W, MITCHELL A, KENNEDY J A, et al. A 17-gene stemness score for rapid determination of risk in acute leukaemia [J]. Nature, 2016, 540(7633): 433-7.
SCHAAFSMA E, FUGLE C M, WANG X, et al. Pan-cancer association of HLA gene expression with cancer prognosis and immunotherapy efficacy [J]. Br J Cancer, 2021, 125(3): 422-32.
JIANG P, GU S, PAN D, et al. Signatures of T cell dysfunction and exclusion predict cancer immunotherapy response [J]. Nat Med, 2018, 24(10): 1550-8.
KLUTH L A, BLACK P C, BOCHNER B H, et al. Prognostic and Prediction Tools in Bladder Cancer: A Comprehensive Review of the Literature [J]. Eur Urol, 2015, 68(2): 238-53.
LIU K, ZHAO K, WANG L, et al. The prognostic values of tumor-infiltrating neutrophils, lymphocytes and neutrophil/lymphocyte rates in bladder urothelial cancer [J]. Pathol Res Pract, 2018, 214(8): 1074-80.
OKAMURA Y, SUGIURA T, ITO T, et al. Neutrophil to lymphocyte ratio as an indicator of the malignant behaviour of hepatocellular carcinoma [J]. Br J Surg, 2016, 103(7): 891-8.
GU X B, TIAN T, TIAN X J, et al. Prognostic significance of neutrophil-to-lymphocyte ratio in non-small cell lung cancer: a meta-analysis [J]. Sci Rep, 2015, 5: 12493.
SUN J, CHEN X, GAO P, et al. Can the Neutrophil to Lymphocyte Ratio Be Used to Determine Gastric Cancer Treatment Outcomes? A Systematic Review and Meta-Analysis [J]. Dis Markers, 2016, 2016: 7862469.
HU K, LOU L, YE J, et al. Prognostic role of the neutrophil-lymphocyte ratio in renal cell carcinoma: a meta-analysis [J]. BMJ Open, 2015, 5(4): e006404.
WANG J, ZHOU X, LIU Y, et al. Prognostic significance of neutrophil-to-lymphocyte ratio in diffuse large B-cell lymphoma: A meta-analysis [J]. PLoS One, 2017, 12(4): e0176008.
ZHANG Q, YANG Q, WENG Y, et al. Neutrophil-to-lymphocyte ratio correlates with prognosis and response to chemotherapy in patients with non-M3 de novo acute myeloid leukemia [J]. Transl Cancer Res, 2021, 10(2): 1013-24.
KENNEDY A, SAHU K K, CERNY J. Role of Immunomodulation of BCG Therapy on AML Remission [J]. Int Med Case Rep J, 2021, 14: 115-9.
KLIMIANKOU M, UENALAN M, KANDABARAU S, et al. Ultra-Sensitive CSF3R Deep Sequencing in Patients With Severe Congenital Neutropenia [J]. Front Immunol, 2019, 10: 116.
GERMESHAUSEN M, BALLMAIER M, WELTE K. Incidence of CSF3R mutations in severe congenital neutropenia and relevance for leukemogenesis: results of a long-term survey [J]. Blood, 2006, 109(1): 93-9.
ULLAH R, YIN Q, SNELL A H, et al. RAF-MEK-ERK pathway in cancer evolution and treatment [J]. Semin Cancer Biol, 2022, 85: 123-54.
FABBRO D, COWAN-JACOB S W, MOEBITZ H. Ten things you should know about protein kinases: IUPHAR Review 14 [J]. Br J Pharmacol, 2015, 172(11): 2675-700.
SUBBIAH V, BAIK C, KIRKWOOD J M. Clinical Development of BRAF plus MEK Inhibitor Combinations [J]. Trends Cancer, 2020, 6(9): 797-810.
BINDELS L B, DEWULF E M, DELZENNE N M. GPR43/FFA2: physiopathological relevance and therapeutic prospects [J]. Trends Pharmacol Sci, 2013, 34(4): 226-32.
BINDELS L B, PORPORATO P E, DUCASTEL S, et al. Ffar2 expression regulates leukaemic cell growth in vivo [J]. Br J Cancer, 2017, 117(9): 1336-40.
ZHAO Z, QIN J, QIAN Y, et al. FFAR2 expressing myeloid-derived suppressor cells drive cancer immunoevasion [J]. J Hematol Oncol, 2024, 17(1): 9.
BORREGO F. The CD300 molecules: an emerging family of regulators of the immune system [J]. Blood, 2013, 121(11): 1951-60.
SILVA R, MOIR S, KARDAVA L, et al. CD300a is expressed on human B cells, modulates BCR-mediated signaling, and its expression is down-regulated in HIV infection [J]. Blood, 2011, 117(22): 5870-80.
VITALLÉ J, TARANCÓN-DÍEZ L, JIMÉNEZ-LEON M R, et al. CD300a identifies a CD4+ memory T cell subset with a higher susceptibility to HIV-1 infection [J]. Aids, 2020, 34(8): 1249-52.
TANG Z, CAI H, WANG R, et al. Overexpression of CD300A inhibits progression of NSCLC through downregulating Wnt/β-catenin pathway [J]. Onco Targets Ther, 2018, 11: 8875-83.
GAUR P, RAHIMLI ALEKBERLI F, KARRA L, et al. Dexamethasone and CD300a activation display additive inhibitory effect on human and murine mast cell functions [J]. Clin Exp Allergy, 2021, 51(10): 1383-6.
KARRA L, GANGWAR R S, PUZZOVIO P G, et al. CD300a expression is modulated in atopic dermatitis and could influence the inflammatory response [J]. Allergy, 2019, 74(7): 1377-80.
ZHUANG H, LI F, SI T, et al. High Expression of CD300A Predicts Poor Survival in Acute Myeloid Leukemia [J]. Acta Haematol, 2023, 146(3): 196-205.
VALIATE B V S, QUEIROZ-JUNIOR C M, LEVI-SCHAFFER F, et al. CD300a contributes to the resolution of articular inflammation triggered by MSU crystals by controlling neutrophil apoptosis [J]. Immunology, 2021, 164(2): 305-17.
BINDELS L B, DEWULF E M, DELZENNE N M. GPR43/FFA2: physiopathological relevance and therapeutic prospects [J]. Trends Pharmacol Sci, 2013, 34(4): 226-32.
BOBROWICZ M, KUBACZ M, SLUSARCZYK A, et al. CD37 in B Cell Derived Tumors-More than Just a Docking Point for Monoclonal Antibodies [J]. Int J Mol Sci, 2020, 21(24).
ZHANG Q, HAN Q, ZI J, et al. CD37 high expression as a potential biomarker and association with poor outcome in acute myeloid leukemia [J]. Biosci Rep, 2020, 40(5).
PAYANDEH Z, NOORI E, KHALESI B, et al. Anti-CD37 targeted immunotherapy of B-Cell malignancies [J]. Biotechnol Lett, 2018, 40(11-12): 1459-66.
ROSENWALD A, WRIGHT G, CHAN W C, et al. The use of molecular profiling to predict survival after chemotherapy for diffuse large-B-cell lymphoma [J]. N Engl J Med, 2002, 346(25): 1937-47.
LE DIEU R, TAUSSIG D C, RAMSAY A G, et al. Peripheral blood T cells in acute myeloid leukemia (AML) patients at diagnosis have abnormal phenotype and genotype and form defective immune synapses with AML blasts [J]. Blood, 2009, 114(18): 3909-16.
COSTELLO R T, SIVORI S, MARCENARO E, et al. Defective expression and function of natural killer cell-triggering receptors in patients with acute myeloid leukemia [J]. Blood, 2002, 99(10): 3661-7.
VEENBERGEN S, VAN SPRIEL A B. Tetraspanins in the immune response against cancer [J]. Immunol Lett, 2011, 138(2): 129-36.

No competing interests reported.

Download PDF

Reviewers agreed at journal
05 Nov, 2024
Reviewers invited by journal
30 Aug, 2024
Editor assigned by journal
30 Aug, 2024
Editor invited by journal
22 Aug, 2024
Submission checks completed at journal
19 Aug, 2024
First submitted to journal
03 Aug, 2024

You are reading this latest preprint version

Neutrophil-related genes predict prognosis and contribute to immunosuppression in acute myeloid leukemia

Status:

Version 1

Abstract

Figures

1 Introduction

2 Materials and methods

2.1 Data source and Processing

2.2 Source of neutrophil-related gene

2.3 Construction and validation of the risk model

2.4 Nomogram construction

2.5 Survival analysis

2.6 Analysis of differentially expressed genes

2.7 Pathway enrichment analysis

2.8 Immune cell infiltration

2.9 Statistical analysis

3 Results

3.1 Establishment of the neutrophil-based prognostic model

3.2 Validation of the neutrophil-based prognostic model

3.3 Construction of predictive nomogram

3.4 Correlation of neutrophil-related genes with prognosis in AML

3.5 Association of CD37 expression with immunity and cellular structure

3.6 CD37-shaped immune disorders context

3.7 Predictive value of CD37 in responsiveness to immunotherapies

4 Discussion

Declarations

Conflicts of Interests

Ethical Statement

Conflict of Interest:

Publisher’s Note

Funding

Author Contribution

Data Availability

References

Additional Declarations

Supplementary Files

Status:

Version 1