Adaptive immunity-related gene expression profile is correlated with clinical phenotype in patients with acute myeloid leukemia.

BACKGROUND: Acute myeloid leukemia (AML) is a common and lethal hematopoietic malignancy that is highly dependent on the immune microenvironment. However, light has yet to be shed on the landscape of adaptive immunity-related genes. This work aimed to uncover the novel molecular events in AML and potential therapeutic strategies for AML treatment.
METHODS: For the current research, the transcriptional information of 732 genes that participate in adaptive immunity was collected from 173 patients with AML, and the patients were grouped into different cohorts based on the different expression patterns. The correlations between gene expression and clinical characteristics, including prognosis, were studied.
RESULTS: According to the notably different expressions of adaptive immunity-related genes, the 173 patients were divided into 2 clusters and 3 subclusters. No significant differences in overall survival (OS) or progression-free survival (PFS) were detected between the clusters or subclusters. There were obvious discrepancies found in age, peripheral blood (PB) blast percentage, and French-American-British (FAB) classification between each cluster or subcluster. The patients in cluster 1 were older and more of them had M5 type; the patients in cluster 2 were younger and more of them had M2 type. Further, 81 genes were significantly correlated with age and 101 genes were significantly correlated with PB blast percentage. Comparison of the prognosis between each FAB type revealed that patients with M3 type displayed the most favorable OS and PFS. Among the differentially expressed genes (DEGs), CLEC2B expression was much lower in M2 patients than in patients with other types (P<0.001), and its high expression indicated a worse outcome (12.4 vs. 46.5 months of OS).
CONCLUSIONS: This study has uncovered the expression profile of adaptive immunity-related genes in AML. The different gene expression patterns are not associated with survival, but are significantly correlated the FAB types. CLEC2B expression is low in patients with M2 type and is negatively correlated with prognosis, thus revealing a potential therapeutic target for AML. 2021 Annals of Translational Medicine. All rights reserved.

Introduction

Acute myeloid leukemia (AML) is a genetically heterogeneous disease which is characterized by abnormal hematopoiesis caused by the accumulation and expansion of immature myeloid cells in the bone marrow and peripheral blood (PB) (1). Although chemotherapy can prolong the prognosis of AML, most patients suffer relapse (2). Allogeneic transplantation is an effective therapy for AML, but relapse after transplantation occurs commonly, leading to a dismal outcome (3). Progress in the genetic investigation of AML has seen researchers excavate the etiology of leukemogenesis, which has promoted the development of potential novel targeted therapies (4). Hematopoietic cell-specific immune molecules, like CD33, have been found to be promising and ideal targets for AML treatment (5,6). Previous studies have revealed that the aberrant expression of immunity-related genes is closely associated with clinical outcome in patients with AML (7,8). Therefore, a comprehensive understanding of immune-related genes and molecular genetic characteristics is essential for diagnosis and outcome prediction in AML.

In multiple tumors, the immune system can identify and eliminate cancerous cells in a process termed cancer immunosurveillance. However, some studies have evidenced that the immune system can also promote tumor progression (9). In humans, innate immunity constitutes the primary facility of host defense against various anomalies, and it also guides the adaptive immune response (10). AML is generally recognized as an immunoresponsive malignancy and remains the most common indication for receiving allogeneic hematopoietic stem cell transplantation (11). Previous studies have shown that the suppression of immunity-related antigens in AML cells is correlated with its relapse (12,13). Previous studies showed that immunophenotype is not only useful in the initial diagnosis but also in monitoring and determining prognosis of the disease through minimal residual disease (MRD) testing (14). In a study by Wang et al., the results showed that the abnormal antigen expression of AML is tightly linked with karyotype abnormality. Detection of immunophenotype may help to diagnose and classify AML (15). However, it is still unclear how the expression profile of adaptive immune-related genes affects leukemia development or prognosis.

The current study investigated the expression characteristics and landscape of adaptive immunity-related genes in a previously reported AML population. The subgroups with specific expression patterns were identified and their correlation with survival was assessed. Additionally, we studied the associations of adaptive immunity-related gene expression with clinical features, and evaluated the prognostic roles of individual genes after the detection of differentially expressed genes (DEGs). This work aimed to uncover the novel molecular events in AML and potential therapeutic strategies for AML treatment. We present the following article in accordance with the REMARK reporting checklist (available at http://dx.doi.org/10.21037/atm-21-2720).

Methods

Patients and RNA-seq

Normalized RNA sequencing data of 200 patients with AML and their corresponding clinical records were obtained from The Cancer Genome Atlas (TCGA) database (cbioportal.org) (16). Detailed clinical characteristics recorded including patients’ age, sex, tumor grade, pathological information, and laboratory test results were collected from the TGCA database. For all 200 patients, AML was diagnosed and validated according to the results of pathological examination. The level of gene expression was expressed as the z-score of messenger RNA (mRNA) and comparisons were conducted between each participant. Due to the lack of related transcriptional information, 27 cases were ruled out of the final analysis, and finally, 173 eligible cases were enrolled into the study. These datasets were publicly available, and ethical approval was granted prior to their use. The study was conducted in accordance with the Declaration of Helsinki (as revised in 2013).

Bioinformatics

Similar to previous studies, the present study focused on genes related to adaptive immunity which are statistically annotated in the Kyoto Encyclopedia of Genes and Genomes (KEGG) database (kegg.jp/hsa04514) (17). The profile of 141 genes involved in cell adhesion in AML were assessed, and 8 of them were excluded due to a lack of expression information. To distinguish samples based on their gene expression profile, a cluster analysis of the genetic expression of integral gene sets was performed. Cases with similar gene expression patterns among the entire study population were identified. The transcription levels of related genes were expressed as mRNA z-scores and clustered using a hierarchical clustering algorithm via a Stanford program, as previously described (18). A tumor stage-specific cluster heat map and pattern were generated with Java TreeView (jtreeview.sourceforge.net) (19) and GraphPad Prism (Version 8.0, GraphPad Software, Inc., San Diego, CA, USA).

Prognostic relevance analyses

The prognostic roles of chromatin remodeling-related genes were investigated through comparison of their survivals between groups with different gene expression levels. The 4 major outcomes included overall survival (OS), progression-free survival (PFS), disease-free survival (DFS), and disease-specific survival (DSS). These outcomes were accessed using GraphPad Prism (Version 8.0, GraphPad Software, Inc.). The survival rates of different clusters were compared to reveal the relationship between gene expression levels and prognosis. Additionally, analysis of the difference in OS between cohorts with low or high expression levels of specific genes was conducted in GraphPad Prism (Version 8.0, GraphPad Software, Inc.).

Statistical analysis

SPSS 24.0 (IBM, Inc., Armonk, NY, USA) was used to carry out statistical analyses. Continuous variables were expressed as mean ± standard deviation (SD). Categorical variables were expressed as numbers and compared using the χ2 test or Fisher’s exact test. Differences in gene expression levels between clusters were detected by analysis of variance (ANOVA). Correlations between variables were determined through regression analyses. Survival curves were plotted and compared between groups using the log-rank test in GraphPad Prism (Version 8.0, GraphPad Software, Inc.). P<0.05 was considered to indicate statistical significance.

Results

The expression profile of adaptive immunity-related genes is significantly associated with clinical characteristics of AML

According to hierarchical clustering, patients with AML were grouped into cohorts on the basis of having similar adaptive immunity-related gene expression profile. There were 2 clusters and 3 subclusters identified for all 173 patients, of which cluster 1 and subcluster 1 contained the same patients. The prognosis of each cohort was compared, and no significant difference was apparent in OS or PFS (P>0.05). Meanwhile, cluster_2 and subcluster_3 showed a slightly worse PFS than the other cohorts, although the difference was not statistically significant (P>0.05) (, ). Comparison of patients’ clinical features revealed that significant differences existed in age, PB blast percentage, and French-American-British (FAB) classification distribution among both the clusters and subclusters (P<0.05). Cluster_1/subcluster_1 had a significantly higher average age, a lower PB blast percentage, and more patients with M4 or M5 type ().

Figure 1

The expression profile of adaptive immunity-related genes is significantly associated with AML clinical characteristics. An insignificant difference in overall survival (OS) and progression-free survival (PFS) between the clusters (A) and subclusters (B) (P>0.05).

Table 1

Survival differences in the cohorts

Identified cohort	OS	Hazard ratio (95% CI of ratio)	P value	PFS	Hazard ratio (95% CI of ratio)	P value
Cluster_1	14.5	1.176 (0.757–1.827)	0.449	13.4	1.355 (0.806–2.277)	0.214
Cluster_2	20.5			19.3
Subcluster_1	14.5		0.687	13.4		0.461
Subcluster_2	16.4			14.2
Subcluster_3	22.3			19.3

OS, overall survival; PFS, progression-free survival.

Table 2

Comparison of clinical features between each cluster or subcluster

Clinical features	Cluster_1 (n=44)	Cluster_2 (n=129)	P value	Subcluster_1 (n=44)	Subcluster_2 (n=37)	Subcluster_3 (n=92)	P value
Age	60.89±13.41	53.29±16.59	0.007	60.89±13.41	55.05±13.92	52.59±17.57	0.019
Sex
Female	24	57	0.294	24	13	44	0.210
Male	20	72		20	24	48
BM blast percentage	68.64	69.24	0.857	68.64	67.89	69.78	0.866
WBC	45.18±38.08	33.72±48.30	0.155	45.18±38.08	35.12±44.66	33.15±49.92	0.356
PB blast percentage	26.93	44.02	0.002	26.93	35.67	47.36	0.002
FAB type
M0	0	16	0.000	0	7	9	0.000
M1	7	37		7	4	33
M2	4	34		4	9	25
M3	0	16		0	1	15
M4	17	17		17	12	5
M5	16	2		16	1	1
M6	0	2		0	1	1
M7	0	3		0	1	2
NC	0	2		0	1	1

BM, bone marrow; WBC, white blood cell; PB, peripheral blood; FAB, French-American-British.

Significant associations exist between adaptive immunity-related gene expression and age and FAB types in AML

Additionally, we detected the specific differences between each identified cohort. Patients in cluster 1/subcluster 1 were older than patients in the other subgroups (). The dominant FAB types in cluster 1 were M4 (38.64%) and M5 (36.36%), while the dominant FAB types in cluster 2 were M1 (28.68%) and M2 (26.36%) (). Compared to subcluster 1, more patients with M4 type and M2 type were included in subcluster 2, and more with M1 type and M2 type were included in subcluster 3 ().

Figure 2

The subgroups according to adaptive immunity-related gene expression are significantly associated with age and FAB type in AML. (A,B) Patients in cluster 1 and subcluster 1 were older than patients in other subgroups; (C) the dominant FAB types in cluster 1 were M4 (38.64%) and M5 (36.36%), while in cluster 2, the dominant FAB types were M1 (28.68%) and M2 (26.36%); (D) compared to subcluster 1, more patients with M4 type and M2 type were included in subcluster 2, and more patients with M1 type and M2 type were included in subcluster 3. *, P<0.05; **, P<0.01. FAB, French-American-British; AML, acute myeloid leukemia.

The expression of adaptive immunity-related genes is significantly correlated with age change and PB blast percentage

We found that cluster_1 had a lower PB blast percentage (26.93 vs. 44.02, P=0.002) and a higher mean age (60.89 vs. 53.29, P=0.006) than the cluster_2 (). To assess the age-associated molecules, the correlation of adaptive immunity-related gene expression with age was evaluated. There were 81 genes markedly correlated with age, including 61 positively correlated genes and 20 negatively correlated genes (). Among these age-associated molecules, F-box protein 32 (FBXO32), cytotoxic T-lymphocyte-associated protein 4 (CTLA4), and killer cell lectin-like receptor F1 (KLRF1) had the strongest positive correlations, and ORAI calcium release-activated calcium modulator 2 (ORAI2), cell division cycle 23 (CDC23), and cluster of differentiation 96 (CD96) had the strongest negative correlations. Assessment of the relationship of adaptive immunity-related genes with PB blast percentage showed that there were 101 genes notably correlated with PB blast percentage, including 6 positively correlated genes and 95 negatively correlated genes (). Among these PB blast percentage-associated molecules, zinc and ring finger 1 (ZNRF1), ubiquitin-conjugating enzyme E2 Z (UBE2Z), and cluster of differentiation 99 (CD99) had the strongest positive correlations, and tubulin alpha 4A (TUBA4A), cathepsin H (CTSH), and protein tyrosine phosphatase receptor type J (PTPRJ) had the strongest negative correlations.

Figure 3

Adaptive immunity-related gene expression was significantly correlated with age and peripheral blood blast percentage. (A) There were 81 genes markedly correlated with age, including 61 positively correlated genes and 20 negatively correlated genes; (B) there were 101 genes markedly correlated with PB blast percentage, including 6 positively correlated genes and 95 negatively correlated genes. PB, peripheral blood.

The expression of adaptive immunity-related genes is correlated with AML FAB types

Since the distributions of FAB subtypes differed significantly for each gene expression pattern, we next studied the correlation of adaptive immunity-related gene expression with FAB types. We compared OS and PFS between patients with each FAB type and uncovered those patients with M3 type showed the most favorable outcome (). Assessment of DEGs between M3 type and the other types showed that 306 genes were significantly varied (P <0.05) and functionally related (). KEGG functional enrichment analysis of these genes revealed that phagosome function, antigen processing and presentation, and ubiquitin-mediated proteolysis were the 3 most significant functional sets among DEGs in patients with M3 type (). Among the DEGs, C-type lectin domain family 2 member B (CLEC2B) showed a significantly low expression level in M3 type and a high expression in the other types (). Moreover, high expression of CLEC2B indicated an unfavorable OS of AML (12.40 vs. 46.50 months, P=0.003) ().

Figure 4

Expression of adaptive immunity-related genes was correlated with AML FAB types. (A) Patients with M3 type showed the most favorable outcome; (B) 306 genes were significantly varied and functionally related; (C) phagosome function, antigen processing and presentation, and ubiquitin-mediated proteolysis were the 3 most significant functional sets among DEGs in the M3 type; (D) CLEC2B showed a significantly low expression level in M3 type but was highly expressed in the other types. (E) high expression of CLEC2B was significantly correlated to an unfavorable overall survival of AML. ****, P<0.0001. FAB, French-American-British; AML, Acute myeloid leukemia; DEGs, differentially expressed genes.

Discussion

Variation in adaptive immunity-related gene expression is significantly associated with clinical characteristics and relapse in patients with AML. Therefore, illuminating the landscape of adaptive immunity-related gene expression can greatly facilitate the understanding of AML. For example, some studies showed that CD25 is an independent prognostic factor in AML patients and alternative therapies for CD25-positive AML patients are needed (20). In the present work, we assessed an expression profile containing 732 genes and grouped 173 patients into different cohorts according to the expression pattern. Comparison of outcomes among the patients showed that the identified cluster_1 (or subcluster_1) had a shorter OS (14.5 vs. 20.5 months) and PFS (13.4 vs. 19.3 months). Nevertheless, no statistically significant difference was detected (P>0.05). This negative result might be attributable to the limited sample size. Regarding clinical characteristics, patients in cluster_1 were older on average (60.89 vs. 53.29 years), which supported cluster_1 having a worse prognosis. PB blast percentage and blood marrow blast percentage are considered to be important prognostic factors for AML, along with the relative ratio (21,22). According to our data, no difference was detected in the BM blast percentage between each cluster or subcluster. However, cluster_1 was associated with a remarkably lower PB blast percentage than the other clusters.

To analyze the correlations of adaptive immunity-related gene expression with age and PB blast percentage, we performed correlation analyses. There were 81 adaptive immunity-related genes that were obviously correlated with age including 61 positively correlated genes and 20 negatively correlated genes. A relationship between differences in adaptive immunity-related genes and older age was revealed. There were 101 genes that were correlated with PB blast percentage. Interestingly, among these PB blast percentage-associated genes, only 6 genes were positively correlated with increased PB blast percentage. The expression of most PB blast percentage-associated genes was negatively correlated with increased PB blast percentage . These results suggest that expression differences in adaptive immunity-related genes are closely associated with age and PB blast percentage in AML.

We uncovered that the distribution of FAB subtypes differed significantly between the cohorts identified based on adaptive immunity-related gene expression. In cluster 1, the most widely distributed FAB types were M4 and M5, while in cluster 2, M1 and M2 were dominantly distributed. These data indicated a significant correlation between adaptive immunity and FAB subtypes. Previous studies have shown that gene aberrations were closely related to the distribution of FAB types and associated with AML prognosis (23). In the current work, we found that patients with M3 type exhibited a prolonged OS and PFS compared to patients with non-M3 types, which is consistent with the results of a previous report (24). Additional assessment of the difference in adaptive immunity-related gene expression between M3 type and non-M3 types indicated that CLEC2B expression differed significantly. In patients with the M3 type, CLEC2B expression was low, while in patients with non-M3 types, it was significantly elevated. Moreover, the high expression of CLEC2B was associated with a worse OS in patients with AML. CLEC-2 is activated by the transmembrane protein podoplanin, which is found outside of the vasculature and is upregulated in development, inflammation, and cancer, and there is also evidence for additional ligands (25). In melanoma, CLEC2B was reported to be highly expressed in melanoma tissues and associated with a poor prognosis (26). Together, these results suggest that CLEC2B has the potential to serve as a prognostic marker for predicting the survival of AML.

In conclusion, our study has demonstrated that adaptive immunity-related genes showed significant expression differences in AML, and the expression patterns were closely associated with age, FAB type, and PB blast percentage in patients with AML. CLEC2B was specifically lowly expressed in patients with M3 type, and its high expression predicted a poor outcome; therefore, it might be pursued as a prognostic factor and treatment target. As for limitation, the results of our present study need to be further investigated in clinical trials in the future.

The article’s supplementary files as