- 1Department of Hematology, The Affiliated People’s Hospital of Jiangsu University, Zhenjiang Clinical Research Center of Hematology, Zhenjiang, Jiangsu, China
- 2Institute of Hematology, Jiangsu University, Zhenjiang, Jiangsu, China
- 3The Key Lab of Precision Diagnosis and Treatment in Hematologic Malignancies of Zhenjiang City, Zhenjiang, Jiangsu, China
- 4Department of Radiation Oncology, The Affiliated People’s Hospital of Jiangsu University, Zhenjiang, Jiangsu, China
Background: Aberrant expression of Homeobox (HOX) genes has been observed in acute myeloid leukemia (AML), but their epigenetic regulatory mechanisms remain largely elusive. Previously, we identified HOXA9 hypomethylation, among HOXA family genes as an epigenetic biomarker for predicting clinical outcomes and guiding treatment choices in AML. Herein, we further investigated the methylation of HOXB family members in AML and determined its clinical implications.
Methods: We first systematically analyzed the association of HOXB methylation with expression and clinical outcomes in AML from The Cancer Genome Atlas (TCGA) database. Next, the candidate prognosis-related gene HOXB3 was selected for clinical relevance analysis and further verified in another independent cohort from our hospital.
Results: Hypomethylation of HOXB3, which was negatively associated with its expression, was correlated with adverse prognosis among HOXB family genes in AML from TCGA datasets. Clinically, AML patients with HOXB3 hypomethylation had unique clinical subtypes and cytogenetic/molecular patterns, including FAB-M5, normal karyotype, cytogenetic/molecular-intermediate risks, and mutations in FLT3-ITD, NPM1 and DNMT3A. Despite these associations, HOXB3 hypomethylation was an independent prognostic biomarker for AML. Bioinformatics analysis demonstrated the association of HOXB3 hypomethylation with several leukemia-related genes (HOXB family genes, miR-10, miR-196a, miR-1, miR-193b and miR-379) in AML. Subsequently, we further validated aberrant HOXB3 hypomethylation and its epigenetic regulatory role in AML.
Conclusions: HOXB3 hypomethylation, which is associated with HOXB3 overexpression, is a frequent event in AML. AML with HOXB3 hypomethylation usually has unique genetic patterns such as a normal karyotype, cytogenetic/molecular-intermediate risk, and mutations in FLT3-ITD, NPM1 and DNMT3A. Despite these associations, HOXB3 hypomethylation may serve as an independent prognostic biomarker for AML.
Introduction
Acute myeloid leukemia (AML) is a heterogeneous aggressive blood cancer with an unsatisfactory prognosis and 5-year overall survival (OS) rate of less than 50% (1). Genetic alterations and epigenetic modifications lead to abrogated differentiation, blocked apoptosis, and uncontrolled proliferation of hematopoietic stem cells (1). Over the past two decades, we have gained deeper insights into the roles of cytogenetics and molecular biology in the pathogenesis and progression of AML (1, 2). In recent years, epigenetic modifications, such as aberrant DNA methylation, a major biological process in silencing the transcription of tumor suppressor genes, have contributed to the pathophysiology of myeloid malignancies and hypomethylating agents are commonly used to treat AML and myelodysplastic syndromes/neoplasm (MDS) (3).
Homeobox (HOX) genes belong to a family of homeodomain containing transcription factors that are mainly involved in early development and hematopoiesis (4). A total of 39 members have been identified in the human HOX gene family, which is divided into four families: HOXA, HOXB, HOXC, and HOXD (4). HOXA family genes, particularly HOXA9 have been thoroughly studied in AML (5). Aberrant methylation was preferentially found in HOXA gene sequences, including HOXA1, HOXA4, HOXA5, HOXA9, HOXA10, and HOXA11 which affect their expression (6). Moreover, our previous study also identified HOXA9 hypomethylation in HOXA family as an epigenetic biomarker for predicting prognosis and guiding treatment choice in AML (7). Among HOXB family, although some studies have suggested that several HOXB genes are upregulated in certain types of AML (8–10), their epigenetic dysregulation in AML has not yet been elucidated. Following our previous study, we systematically investigated methylation of HOXB family members in AML and determined their clinical significance.
Materials and methods
Public datasets
This study included a cohort of 200 AML patients from The Cancer Genome Atlas (TCGA) public databases (11). Clinical data and laboratory findings of AML patients were downloaded from the cBioPortal (https://www.cbioportal.org/) (12). HOXB family gene expression data was available for 173 patients, whereas HOXB family gene methylation data was available for 194 patients (170 overlapping patients with both HOXB expression and methylation).
Patients
A total of 54 de novo AML patients were also included in this research as the validation cohort. Clinico-pathological features of the patients were presented in Supplementary Table S1. After informed consents were obtained from all participants, bone marrow (BM) samples were collected and were prepared for BM mononuclear cells (BMMNCs) as reported (13, 14). Twenty healthy individuals served as the controls.
Cell lines and cell culture
Human leukemic cell lines (KASUMI-1, OCI-AML-2, OCI-AML-3, THP-1, MV4-11, SKM-1, NOMO-1, HEL and DAMI) and BM stoma cell line (HS-5) were used in this study. Cell culture was performed as previous described (7).
Hypomethylation agent and demethylation studies
The KASUMI-1 and OCI-AML-2 cells in 4 mL (5×105 cells/mL) were treated with the hypomethylation agent 5-aza-2’-deoxycytidine (5-aza-dC) at a final concentration of 0 μM (control), 2 μM, 4 μM, and 10 μM for 3 days. All the cell lines were cultured until they are harvested for RNA and DNA extraction.
RNA isolation, reverse transcription and RQ-PCR
Total RNA was isolated from cells using TRIzol reagent, followed by reverse transcription for cDNA synthesis as previously reported (13, 14). Real-time quantitative PCR (RQ-PCR) was conducted to examine the expression levels of HOXB3 by the usage of TB Green® Premix Ex Taq™ II (Tli RNaseH Plus) (RR820A, TaKaRa, Tokyo, Japan). The primers used for the detection of HOXB3 expression were 5’-CTTGGACCGGCTGTTGG-3’ (Forward) and 5’-TTGTCGTAGTAGGTGGCTTT-3’ (Reverse). The RQ-PCR conditions were as follows: 95°C for 5 min, 40 cycles for 10 s at 95°C, 30 s at 56°C, 1 min at 72°C. The reference gene GAPDH was used as our previous literature (7). The relative HOXB3 expression levels were measured using 2-ΔΔCt method.
DNA isolation, bisulfite modification, and RQ-MSP
Genomic DNA was extracted from the cells using a Puregene Blood Kit (158023, QIAGEN, Duesseldorf, Germany), followed by bisulfite modification with the usage of EZ DNA Methylation-Gold™ Kit (D5002, ZYMO RESEARCH, Orange County, California). Real-time quantitative methylation-specific PCR (RQ-MSP) was performed to detect the methylation level of HOXB3 by the usage of TB Green® Premix Ex Taq™ II (Tli RNaseH Plus) (RR820A, TaKaRa, Tokyo, Japan). The primers utilized for the detection of HOXB3 methylation were 5’-TTTCGGATCGTTTATACGC-3’ (Methylation Forward) and 5’-CACTTCATACGCCGATTCTA-3’ (Methylation Reverse) as well as 5’-GTTTTTGGATTGTTTATATGT-3’ (Unmethylation Forward) and 5’-CCACTTCATACACCAATTCTAA-3’ (Unmethylation Reverse). The conditions of RQ-MSP program were 95°C for 5 min, 40 cycles for 10 s at 95°C, 30 s at 56°C (Methylation primers) or 54 °C (Unmethylation primers), 30 s at 72°C, 30s at 75 °C. The reference gene ALU was detected as our previous literature (7). The Relative HOXB3 methylation levels were calculated using the 2-ΔΔCT method.
Bioinformatics analysis
The detailed bioinformatics analysis referred to our previous article (7).
Statistical analysis
SPSS 22.0 and GraphPad 8.0 software were used for the statistical analysis. The Mann-Whitney U-test/Kruskal-Wallis test/Student’s T-test and Pearson’s χ2-test/Fisher’s exact test were used to compare the continuous and categorical variables, respectively. Kaplan-Meier method (log-rank test) and Cox regression analysis (proportional hazards) were performed to determine the prognostic impact of HOXB family genes methylation on OS and disease-free survival (DFS). The receiver operating characteristic (ROC) curve and area under the ROC curve (AUC) were conducted to evaluate the ability of HOXB3 methylation to distinguish patients with AML from controls. A two-sided P value less than 0.05 was considered statistically significant.
Results
Identification of HOXB3 hypomethylation among HOXB family genes associated with prognosis with expression in AML
To identify the crucial epigenetic markers in HOXB family in AML, we first analyzed the prognostic value of HOXB members methylation with available data (HOXB2/HOXB3/HOXB4/HOXB5/HOXB6/HOXB7/HOXB8/HOXB9) using TCGA databases. The patients were divided into two groups (hyper- and hypo- methylation) based on the median methylation level of each gene. In total AML patients, hypo-methylation of HOXB3, HOXB5 and HOXB6 was markedly correlated with shorter OS (P = 0.006, 0.004 and 0.011, respectively) and DFS (P = 0.005, 0.005, and 0.009, respectively), whereas the other members did not affect survival (both P>0.05, Figure 1a). As patients with AML-M3 have extremely favorable outcomes, we further explored the prognostic impact of these members in non-M3 AML. Notably, only HOXB3 hypo-methylation had a significant prognostic effect on OS and LFS (both P = 0.009, Figure 1a).
Figure 1. Identification of HOXB family genes methylation correlated with prognosis and expression in AML. (a) The statistical value of the prognostic significance of each HOXB family gene methylation on survival in AML patients (determined by Kaplan-Meier analysis). Moreover, the survival curve of the impact of HOXB3, HOXB5 and HOXB6 methylation on overall/disease-free survival in total/non-M3 AML patients with significant difference were also presented. Total_AML_OS: The prognostic effect of each HOXB family gene methylation on overall survival among total AML; Total_AML_DFS: The prognostic effect of each HOXB family gene methylation on disease free survival among total AML; nonM3_AML_OS: The prognostic effect of each HOXB family gene methylation on overall survival among non-M3 AML; nonM3_AML_DFS: The prognostic effect of each HOXB family gene methylation on disease free survival among non-M3 AML. NS, no significant; *P < 0.05; **P < 0.01. (b) The statistical value of the correlation analysis between each HOXB gene expression and methylation in AML patients (determined by Spearman test). The depth of color differences represents correlation coefficients. *P < 0.05; **P < 0.01; ***P < 0.001.
We further tested the prognostic effect of HOXB expression on OS and LFS in AML patients using Cox regression analysis. In all AML patients, none of the HOXB members had prognostic effects on OS or DFS (all P>0.05, Table 1). However, among non-M3 AML patients, HOXB3 hypo-methylation significantly affected the OS and DFS (both P = 0.009, Table 1).
Table 1. Univariate and multivariate Cox regression analysis of HOXB genes methylation for overall survival and leukemia-free survival in AML patients.
Promoter-associated DNA methylation is a major mechanism regulating gene expression. We analyzed the association between HOXB methylation and expression in patients with AML. A significant negative correlation was observed for HOX2/HOXB3/HOXB4/HOXB5/HOXB6/HOXB7 but not for HOXB8/HOXB9 (Figure 1b).
Taken together, these results suggest that HOXB3 hypo-methylation may be the most valuable prognostic biomarker among the HOXB genes in AML, especially in non-M3 AML, and was selected for further analysis.
HOXB3 hypomethylation associated with specific subtypes of AML
To evaluate the clinical implications of HOXB3 methylation in AML, we compared the clinico-pathological characteristics between HOXB3 hypo- and hyper-methylation groups. There were no significant differences between the two groups in sex, age, white blood cell count, and peripheral blood (PB)/BM blast ratio (all P>0.05, Table 2). However, marked differences were observed between the two groups in the distribution of the French-American-British (FAB) classification (P = 0.009), karyotypic classification (P < 0.001), and molecular/cytogenetic risk classifications (all P < 0.001) (Table 2). Hypomethylation of HOXB3 was significantly occurred in FAB-M5 (P = 0.011, Table 2). Among the karyotypic classification, HOXB3 hypomethylation was more likely to occur in normal karyotype (P < 0.001), and rarely occurred in inv(16) (P < 0.001) (Table 2). As expected, hypomethylation of HOXB3 was occurred markedly in intermediate karyotype and rarely occurred in good karyotype among both cytogenetic and molecular risk classifications (all P ≤ 0.001, Table 2). Among common gene mutations, hypomethylation of HOXB3 was associated with FLT3 (P = 0.010), NPM1 (P < 0.001), DNMT3A mutations (P < 0.001) (Table 2).
Table 2. Correlation of clinic-pathologic characteristics in AML patients with HOXB3 hypo- and hyper-methylation.
Based on the dramatical clinical relations, we further investigated the HOXB3 methylation level in AML patients among different FAB subtypes, karyotypic classifications, and molecular/cytogenetic risks. Expectedly, significant differences were confirmed in the FAB subtypes (P = 0.004), karyotypic classifications (P < 0.001), and molecular/cytogenetic risks (both P < 0.001) (Figures 2a–d). Next, we explored HOXB3 methylation levels in AML patients with and without NPM1/FLT3/DNMT3A/RUNX1 gene mutations. Similar results were observed in NPM1 (P < 0.001), FLT3 (P < 0.001), DNMT3A (P < 0.001) mutation and RUNX1 wild-type (P = 0.009) (Figures 2e–h).
Figure 2. The associations of HOXB3 methylation with specific subtypes of AML. (a) HOXB3 methylation level among different French-American-British (FAB) subtypes of AML. (b) HOXB3 methylation level among different cytogenetic risk subgroups of AML. (c) HOXB3 methylation level among different karyotypic subgroups of AML. (d) HOXB3 methylation level among different molecular risk subgroups of AML. (e) HOXB3 methylation level in groups of AML patients with and without NPM1 mutations. (f) HOXB3 methylation level in groups of AML patients with and without FLT3 mutations. (g) HOXB3 methylation level in groups of AML patients with and without DNMT3A mutations. (h) HOXB3 methylation level in groups of AML patients with and without RUNX1 mutations. The differences between groups were determined by Mann-Whitney U test/Kruskal-Wallis test.
The independent prognostic value of HOXB3 hypomethylation in AML
Since a close correlation was found between HOXB3 hypo-methylation and common prognostic factor karyotypic classifications, we further performed multivariate Cox regression to confirm the independent prognostic impact of HOXB3 hypo-methylation in AML after adjusting for prognosis-related factors. Cox regression multivariate analysis indicated that HOXB3 expression was an independent risk factor affecting OS and DFS in both total AML (P = 0.023 and 0.009, respectively) and non-M3 AML (P = 0.024 and 0.039, respectively) patients (Table 3).
Table 3. Univariate and multivariate Cox regression analysis of potential risk factors for overall survival and leukemia free survival in AML patients.
Molecular alterations correlated with HOXB3 dysregulation in AML
To investigate the biological network correlated with HOXB3 dysregulation in AML, we studied the transcriptomes of AML samples with low and high HOXB3 expression from TCGA database. A total of 784 differentially expressed genes (DEGs, including mRNAs and long noncoding RNAs) (208 positively correlated and 576 negatively correlated) and 81 differentially expressed miRNAs (DEmiRs) (8 positively correlated and 73 negatively correlated) were identified between the two groups (Figures 3a–c; Supplementary Table S2). The positively correlated DEGs (Top 10) were HOXB family members, which were reported to have pro-leukemic effects and were associated with clinical outcomes of AML (8–10, 15). Moreover, positively correlated DEmiRs (top 5) such as miR-10, miR-196a and miR-1 were reported to have pro-leukemic effects (16–18), whereas negatively correlated DEmiRs (top 5) such as miR-193b and miR-379 were reported to have anti-leukemic effects in AML (19, 20). Moreover, Gene Ontology (GO) analysis demonstrated that these DEGs were involved in multiple biological processes such as development (Figure 3d).
Figure 3. Molecular signatures associated with HOXB3 expression in AML. (a) Expression heatmap of differentially expressed genes including mRNAs/lncRNAs (DEGs) between low and high HOXB3 expression groups in AML (|log2 FC|>1.5, FDR<0.05 and P < 0.05). (b) Volcano plot of DEGs between low and high HOXB3 expression groups in AML. (c) Expression heatmap of differentially expressed miRNAs between low and high HOXB3 expression groups in AML (|log2 FC|>1.5, FDR<0.05 and P < 0.05). (d) Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of above DEGs.
Validation of HOXB3 hypomethylation and its epigenetically regulatory role in AML
To validate epigenetic dysregulation of HOXB3 during leukemogenesis, we first determined HOXB3 methylation levels in de novo AML and controls from the validation cohort. HOXB3 methylation levels in AML patients were markedly reduced compared with those in the controls (P = 0.047, Figure 4a). ROC curve analysis indicated that HOXB3 methylation might act as a potential biomarker for distinguishing AML patients from controls, with an AUC of 0.651 (95% CI: 0.532-0.770, P = 0.047) (Figure 4b). Next, we detected the methylation and expression of HOXB3 in nine AML cell lines and one BM stoma cell line. KASUMI-1 and OCI-AML-2 showed the highest HOXB3 methylation levels with a relatively low HOXB3 expression levels, whereas OCI-AML-3 and SKM-1 presented the lowest HOXB3 methylation levels with a relatively high HOXB3 expression levels (Figure 4c-d). HOXB3 expression showed a nearly negative correlation with HOXB3 methylation in nine AML cell lines and one BM stoma cell line (R=-0.612, P = 0.067), indicating that HOXB3 methylation plays a vital role in regulating HOXB3 expression. To confirm the epigenetic regulatory role, KASUMI-1 and OCI-AML-2 cell lines were treated with the hypomethylation agent 5-aza-dC. Expectedly, HOXB3 expression levels were markedly increased after 5-aza-dC treatment (Figures 4e-f).
Figure 4. Validation of aberrant HOXB3 hypomethylation and its epigenetically regulatory role in AML. (a) HOXB3 methylation in AML patients from our hospital. (b) Receiver operating characteristic curve analysis of HOXB3 methylation in distinguishing AML patients from controls. (c) HOXB3 methylation in AML cell lines. (d) HOXB3 expression in AML cell lines. (e) HOXB3 expression in KASUMI-1 cell-line after 5-aza-dC treatment at different dose. (f) HOXB3 expression in OCI-AML-2 cell-line after 5-aza-dC treatment at different dose. The differences between groups were determined by Mann-Whitney U test and Student T test. *P < 0.05; **P < 0.01; ***P < 0.001; ns: no significant.
Discussion
Accumulating evidence has demonstrated aberrant HOX gene expression in diverse blood cancers, especially AML, but the epigenetic regulatory mechanism in AML remains largely elusive (5). Previously, we systematically identified the prognostic value of HOXA family gene methylation in AML, and revealed HOXA9 hypomethylation is an epigenetic biomarker for predicting clinical outcomes and guiding treatment choice in AML (7). Herein, we further determined HOXB family genes methylation in AML, and revealed HOXB3 hypomethylation as an epigenetic biomarker independently affects clinical outcome in AML. Similarly, Lindblad et al. reported that HOXB2 and HOXB3 expressions are independent prognostic markers in AML (9). As we have also demonstrated the epigenetic regulatory role of HOXB3 methylation in AML, we propose that the prognostic value and functional role of HOXB3 hypomethylation in AML mainly depends on its expression. In accordance with our study, Qu et al. showed that decreased methylation of HOXB3 and HOXB4 is associated with increased expression of both HOXB genes specific to mid-risk AML (21). These results suggest HOXB3 methylation and expression are potential biomarkers for AML in clinical use.
Distinct patterns of HOX family gene expression have been described in defined cytogenetic and molecular subsets of patients with AML. Drabkin et al. reported that favorable and poor prognostic chromosomal rearrangements in AML were associated with distinct levels of HOX expression (22). Moreover, HOX overexpression may be correlated with a distinct biologic subset of AML with intermediate cytogenetics (23). Consistently, we observed a marked associations of HOXB3 hypomethylation with FAB subtypes (M5), different karyotypes (normal karyotype), and molecular/cytogenetic risk classifications (intermediate risk). Moreover, we also found that HOXB3 hypomethylation correlated with NPM1, FLT3, and DNMT3A mutations. Similarly, previous studies have indicated that the high expression of several HOX genes is associated with NPM1, FLT3 and DNMT3A mutations (23–26). Moreover, several studies also addressed the closely associations of aberrant HOXB expression/methylation with NPM1, FLT3, and DNMT3A mutations in AML (10, 21, 27). These results indicated that HOXB3 is correlated with certain genetic contexts in AML and may be a common downstream gene with potential targets in these specific AML subtypes.
HOX family members have been identified as transcription factors that activate their own expression or cross-regulate their family members or their cofactors (4). In our study, we also observed that HOXB3 is closely related to the expression of other HOX family genes enriched in the same AML-related function. Moreover, several leukemia-related miRNAs, such as miR-10, miR-196a, miR-1, miR-193b and miR-379, have also been found to be associated with HOXB3 dysregulation (16–20). Consistently, Chen et al. found that HOXB5 is closely associated with the expression of other HOX family genes (10), which supports our results.
Although our study yielded significant results, it has some limitations in the current research. First, the clinical implications of HOXB3 hypomethylation in AML were drawn from public databases and validated using only a small cohort of AML patients from our hospital. These results need to be verified in a larger cohort of AML patients. Second, few functional and mechanical studies have determined the underlying role of HOXB3 hypomethylation in AML. Future studies should incorporate these experiments in order to further corroborate our research findings.
Conclusion
Taken together, HOXB3 hypomethylation, which is associated with HOXB3 overexpression, is a frequent event in AML. AML with HOXB3 hypomethylation usually has unique genetic patterns such as a normal karyotype, cytogenetic/molecular-intermediate risk, and mutations in FLT3-ITD, NPM1 and DNMT3A. Despite these associations, HOXB3 hypomethylation may serve as an independent prognostic biomarker for AML.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.
Ethics statement
The studies involving humans were approved by Ethics Committee of the Affiliated People’s Hospital of Jiangsu University. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
JZ: Conceptualization, Funding acquisition, Methodology, Validation, Writing – review & editing. TZ: Data curation, Formal Analysis, Funding acquisition, Investigation, Writing – original draft. RC: Investigation, Methodology, Writing – original draft. FX: Investigation, Methodology, Writing – original draft. ZX: Formal Analysis, Software, Writing – original draft. MC: Methodology, Writing – original draft. XW: Methodology, Writing – original draft. JQ: Data curation, Funding acquisition, Resources, Writing – original draft.
Funding
The author(s) declared that financial support was received for this work and/or its publication. The work was supported by National Natural Science Foundation of China (82300164, 82270179), Natural Science Foundation of Jiangsu Province (BK20221287, BK20230296), Research Project of Jiangsu Commission of Health (M2022123), Social Development Foundation of Zhenjiang (SH2025041).
Conflict of interest
The authors declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declared that generative AI was not used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2025.1709417/full#supplementary-material
Abbreviations
AML, acute myeloid leukemia; OS, overall survival; MDS, myelodysplastic syndromes; HOX, homeobox; MDS, myelodysplastic syndromes/neoplasm; TCGA, The Cancer Genome Atlas; BM, bone marrow; BMMNCs, BM mononuclear cells; 5-aza-dC, 5-aza-2’-deoxycytidine; RQ-PCR, real-time quantitative PCR; RQ-MSP, Real-time quantitative methylation-specific PCR; DFS, disease-free survival; ROC; receiver operating characteristic; AUC, area under the ROC curve; PB, peripheral blood; FAB, French-American-British; DEGs, differentially expressed genes (including mRNAs and lncRNAs); DemiRs, differentially expressed miRNAs; WHO, World Health Organization; GO, Gene Ontology.
References
1. DiNardo CD, Erba HP, Freeman SD, and Wei AH. Acute myeloid leukaemia. Lancet. (2023) 401:2073–86. doi: 10.1016/S0140-6736(23)00108-3
2. Snaith O, Poveda-Rogers C, Laczko D, Yang G, and Morrissette JJD. Cytogenetics and genomics of acute myeloid leukemia. Best Pract Res Clin Haematol. (2024) 37:101533. doi: 10.1016/j.beha.2023.101533
3. Estey EH. Epigenetics in clinical practice: the examples of azacitidine and decitabine in myelodysplasia and acute myeloid leukemia. Leukemia. (2013) 27:1803–12. doi: 10.1038/leu.2013.173
4. Alharbi RA, Pettengell R, Pandha HS, and Morgan R. The role of HOX genes in normal hematopoiesis and acute leukemia. Leukemia. (2013) 27:1000–8. doi: 10.1038/leu.2012.356
5. Collins CT and Hess JL. Role of HOXA9 in leukemia: dysregulation, cofactors and essential targets. Oncogene. (2016) 35:1090–8. doi: 10.1038/onc.2015.174
6. Roux B, Picou F, Debeissat C, Koubi M, Gallay N, Hirsch P, et al. Aberrant DNA methylation impacts HOX genes expression in bone marrow mesenchymal stromal cells of myelodysplastic syndromes and de novo acute myeloid leukemia. Cancer Gene Ther. (2022) 29:1263–75. doi: 10.1038/s41417-022-00441-w
7. Xie F, Zhang TJ, Zhang XL, Xu ZJ, Qiao L, Wang Y, et al. Identification of HOXA9 methylation as an epigenetic biomarker predicting prognosis and guiding treatment choice in acute myeloid leukemia. BMC Cancer. (2025) 25:215. doi: 10.1186/s12885-025-13633-y
8. Umeda S, Yamamoto K, Murayama T, Hidaka M, Kurata M, Ohshima T, et al. Prognostic significance of HOXB4 in de novo acute myeloid leukemia. Hematology. (2012) 17:125–31. doi: 10.1179/102453312X13376952196250
9. Lindblad O, Chougule RA, Moharram SA, Kabir NN, Sun J, Kazi JU, et al. The role of HOXB2 and HOXB3 in acute myeloid leukemia. Biochem Biophys Res Commun. (2015) 467:742–7. doi: 10.1016/j.bbrc.2015.10.071
10. Chen M, Qu Y, Yue P, and Yan X. The prognostic value and function of HOXB5 in acute myeloid leukemia. Front Genet. (2021) 12:678368. doi: 10.3389/fgene.2021.678368
11. Cancer Genome Atlas Research Network, Ley TJ, Miller C, Ding L, Raphael BJ, Mungall AJ, et al. Genomic and epigenomic landscapes of adult de novo acute myeloid leukemia. N Engl J Med. (2013) 368:2059–74. doi: 10.1056/NEJMoa1301689
12. Cerami E, Gao J, Dogrusoz U, Gross BE, Sumer SO, Aksoy BA, et al. The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. Cancer Discov. (2012) 2:401–4. doi: 10.1158/2159-8290.CD-12-0095
13. Zhou JD, Zhao YJ, Leng JY, Gu Y, Xu ZJ, Ma JC, et al. DNA methylation-mediated differential expression of DLX4 isoforms has opposing roles in leukemogenesis. Cell Mol Biol Lett. (2022) 27:59. doi: 10.1186/s11658-022-00358-0
14. Zhang TJ, Xu ZJ, Wen XM, Gu Y, Ma JC, Yuan Q, et al. SLIT2 promoter hypermethylation-mediated SLIT2-IT1/miR-218 repression drives leukemogenesis and predicts adverse prognosis in myelodysplastic neoplasm. Leukemia. (2022) 36:2488–98. doi: 10.1038/s41375-022-01659-1
15. Kusakabe M, Sun AC, Tyshchenko K, Wong R, Nanda A, Shanna C, et al. Synthetic modeling reveals HOXB genes are critical for the initiation and maintenance of human leukemia. Nat Commun. (2019) 10:2913. doi: 10.1038/s41467-019-10510-8
16. Havelange V, Ranganathan P, Geyer S, Nicolet D, Huang X, Yu X, et al. Implications of the miR-10 family in chemotherapy response of NPM1-mutated AML. Blood. (2014) 123:2412–5. doi: 10.1182/blood-2013-10-532374
17. Fan B, Wang L, Hu T, Zheng L, and Wang J. Exosomal miR-196a-5p secreted by bone marrow mesenchymal stem cells inhibits ferroptosis and promotes drug resistance of acute myeloid leukemia. Antioxid Redox Signal. (2025) 42:933–53. doi: 10.1089/ars.2024.0882
18. Ghazaryan A, Wallace JA, Tang WW, Barba C, Lee SH, Bauer KM, et al. miRNA-1 promotes acute myeloid leukemia cell pathogenesis through metabolic regulation. Front Genet. (2023) 14:1192799. doi: 10.3389/fgene.2023.1192799
19. Bhayadia R, Krowiorz K, Haetscher N, Jammal R, Emmrich S, Obulkasim A, et al. Endogenous tumor suppressor microRNA-193b: therapeutic and prognostic value in acute myeloid leukemia. J Clin Oncol. (2018) 36:1007–16. doi: 10.1200/JCO.2017.75.2204
20. Wu H, Zhao L, Guo H, Xie Y, Hu J, and Tan X. miR-379-5p inhibited the proliferation of acute myeloid leukemia cells through negative regulation of YBX1. Turk J Haematol. (2025) 42:92–9. doi: 10.4274/tjh.galenos.2025.2024.0424
21. Qu X, Davison J, Du L, Storer B, Stirewalt DL, Heimfeld S, et al. Identification of differentially methylated markers among cytogenetic risk groups of acute myeloid leukemia. Epigenetics. (2015) 10:526–35. doi: 10.1080/15592294.2015.1048060
22. Drabkin HA, Parsy C, Ferguson K, Guilhot F, Lacotte L, Roy L, et al. Quantitative HOX expression in chromosomally defined subsets of acute myelogenous leukemia. Leukemia. (2002) 16:186–95. doi: 10.1038/sj.leu.2402354
23. Roche J, Zeng C, Barón A, Gadgil S, Gemmill RM, Tigaud I, et al. Hox expression in AML identifies a distinct subset of patients with intermediate cytogenetics. Leukemia. (2004) 18:1059–63. doi: 10.1038/sj.leu.2403366
24. Andreeff M, Ruvolo V, Gadgil S, Zeng C, Coombes K, Chen W, et al. HOX expression patterns identify a common signature for favorable AML. Leukemia. (2008) 22:2041–7. doi: 10.1038/leu.2008.198
25. Brunetti L, Gundry MC, Sorcini D, Guzman AG, Huang YH, Ramabadran R, et al. Mutant NPM1 maintains the leukemic state through HOX expression. Cancer Cell. (2018) 34:499–512.e9. doi: 10.1016/j.ccell.2018.08.005
26. Tan YT, Sun Y, Zhu SH, Ye L, Zhao CJ, Zhao WL, et al. Deregulation of HOX genes by DNMT3A and MLL mutations converges on BMI1. Leukemia. (2016) 30:1609–12. doi: 10.1038/leu.2016.15
27. Liu HC, Shih LY, May Chen MJ, Wang CC, Yeh TC, Lin TH, et al. Expression of HOXB genes is significantly different in acute myeloid leukemia with a partial tandem duplication of MLL vs. a MLL translocation: a cross-laboratory study. Cancer Genet. (2011) 204:252–9. doi: 10.1016/j.cancergen.2011.02.003
Keywords: AML, HOXB, HoxB3, hypomethylation, prognosis, regulation
Citation: Zhang T-j, Chang R, Xie F, Xu Z-j, Chu M-q, Wu X-c, Qian J and Zhou J-d (2026) Identification and validation of HOXB3 hypomethylation as a novel prognostically epigenetic biomarker in acute myeloid leukemia. Front. Immunol. 16:1709417. doi: 10.3389/fimmu.2025.1709417
Received: 20 September 2025; Accepted: 22 December 2025; Revised: 11 December 2025;
Published: 12 January 2026.
Edited by:
Ananya Datta Mitra, UC Davis Health, United StatesReviewed by:
Irma Olarte, Hospital General de México Dr. Eduardo Liceaga, MexicoAlireza Ghezavati, Quest Diagnostics, United States
Copyright © 2026 Zhang, Chang, Xie, Xu, Chu, Wu, Qian and Zhou. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jing-dong Zhou, emhvdWppbmdkb25nQHVqcy5lZHUuY24=
†These authors have contributed equally to this work
Ran Chang1,2,3,4†