ORIGINAL RESEARCH article

Front. Immunol., 12 January 2026

Sec. Cancer Immunity and Immunotherapy

Volume 16 - 2025 | https://doi.org/10.3389/fimmu.2025.1709417

Identification and validation of HOXB3 hypomethylation as a novel prognostically epigenetic biomarker in acute myeloid leukemia

  • 1. Department of Hematology, The Affiliated People’s Hospital of Jiangsu University, Zhenjiang Clinical Research Center of Hematology, Zhenjiang, Jiangsu, China

  • 2. Institute of Hematology, Jiangsu University, Zhenjiang, Jiangsu, China

  • 3. The Key Lab of Precision Diagnosis and Treatment in Hematologic Malignancies of Zhenjiang City, Zhenjiang, Jiangsu, China

  • 4. Department of Radiation Oncology, The Affiliated People’s Hospital of Jiangsu University, Zhenjiang, Jiangsu, China

Abstract

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% (). Genetic alterations and epigenetic modifications lead to abrogated differentiation, blocked apoptosis, and uncontrolled proliferation of hematopoietic stem cells (). 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 (, ). 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) ().

Homeobox (HOX) genes belong to a family of homeodomain containing transcription factors that are mainly involved in early development and hematopoiesis (). 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 (). HOXA family genes, particularly HOXA9 have been thoroughly studied in AML (). Aberrant methylation was preferentially found in HOXA gene sequences, including HOXA1, HOXA4, HOXA5, HOXA9, HOXA10, and HOXA11 which affect their expression (). Moreover, our previous study also identified HOXA9 hypomethylation in HOXA family as an epigenetic biomarker for predicting prognosis and guiding treatment choice in AML (). Among HOXB family, although some studies have suggested that several HOXB genes are upregulated in certain types of AML (–), 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 (). Clinical data and laboratory findings of AML patients were downloaded from the cBioPortal (https://www.cbioportal.org/) (). 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 (, ). 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 ().

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 (, ). 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 (). 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 (). The Relative HOXB3 methylation levels were calculated using the 2-ΔΔCT method.

Bioinformatics analysis

The detailed bioinformatics analysis referred to our previous article ().

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

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

VariablesTotal AMLNon-M3 AML
Univariate analysisMultivariate analysisUnivariate analysisMultivariate analysis
HR (95% CI)PHR (95% CI)PHR (95% CI)PHR (95% CI)P
Overall Survival
HOXB2 methylation0.869 (0.615-1.229)0.4281.003 (0.703-1.432)0.986
HOXB3 methylation0.613 (0.432-0.869)0.0060.719 (0.487-1.064)0.0990.620 (0.433-0.890)0.0090.620 (0.433-0.890)0.009
HOXB4 methylation0.929 (0.658-1.313)0.6760.999 (0.700-1.425)0.995
HOXB5 methylation0.603 (0.425-0.857)0.0050.700 (0.473-1.037)0.0750.746 (0.521-1.070)0.1120.965 (0.626-1.487)0.872
HOXB6 methylation0.639 (0.451-0.906)0.0120.821 (0.542-1.244)0.3520.706 (0.493-1.012)0.0580.838 (0.563-1.247)0.384
HOXB7 methylation0.813 (0.575-1.149)0.2410.844 (0.591-1.206)0.352
HOXB8 methylation0.992 (0.702-1.402)0.9661.138 (0.797-1.624)0.477
HOXB9 methylation0.943 (0.667-1.334)0.7410.817 (0.731-1.488)0.817
Disease-free Survival
HOXB2 methylation0.855 (0.605-1.208)0.3750.983 (0.689-1.403)0.925
HOXB3 methylation0.608 (0.429-0.802)0.0050.698 (0.477-1.022)0.0640.621 (0.433-0.890)0.0090.621 (0.433-0.890)0.009
HOXB4 methylation0.939 (0.665-1.326)0.7201.000 (0.701-1.427)0.998
HOXB5 methylation0.613 (0.432-0.869)0.0060.708 (0.483-1.036)0.0760.774 (0.541-1.109)0.1620.992 (0.649-1.516)0.970
HOXB6 methylation0.633 (0.447-0.898)0.0100.803 (0.532-1.211)0.2950.713 (0.498-1.022)0.0650.837 (0.566-1.236)0.371
HOXB7 methylation0.771 (0.546-1.090)0.1421.124 (0.736-1.717)0.5890.811 (0.568-1.158)0.249
HOXB8 methylation0.980 (0.694-1.385)0.9101.127 (0.790-1.609)0.509
HOXB9 methylation0.991 (0.702-1.400)0.9591.085 (0.761-1.548)0.652

Univariate and multivariate Cox regression analysis of HOXB genes methylation for overall survival and leukemia-free survival in AML patients.

HR, hazard ratio; CI, confidence interval; Multivariate analysis includes variables with P < 0.200 in univariate analysis.

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

Patient’s parametersTotalHypo- (n=97)Hyper- (n=97)P value
Sex, male/female104/9046/5158/390.113
Median age, years (range)57 (18-88)60 (21-88)55 (18-83)0.100
Median WBC, ×109/L (range)16.6 (0.4-298.4)19.9 (0.5-298.4)16.0 (0.4-223.8)0.071
Median PB blasts, % (range)36 (0-98)22 (0-98)40.50 (0-97)0.522
Median BM blasts, % (range)73.5 (30-100)75 (30-100)71.5 (30-100)0.283
FAB classifications0.009
 M0195140.051
 M14421230.864
 M24224180.384
 M3186120.215
 M44122190.725
 M5221750.011
 M6321>0.999
 M73030.246
 No data2020.497
Karyotypes<0.001
 normal916328<0.001
 t(15;17)166100.435
 t(8;21)7250.444
 inv(16)12012<0.001
 +8945>0.999
 del(5)101>0.999
 -7/del(7)844>0.999
 11q234130.621
 others144100.163
 complex2811170.307
 No data422>0.999
Risks (cytogenetic)<0.001
 Good358270.001
 Intermediate1127042<0.001
 Poor4317260.166
 No data422>0.999
Risks (molecular)<0.001
 Good37829<0.001
 Intermediate1036835<0.001
 Poor5120310.102
 No data312>0.999
Gene mutations
FLT3 (+/-)55/139 (28.4%)36/61 (37.1%)19/78 (19.6%)0.010
NPM1 (+/-)53/141 (27.3%)47/50 (48.5%)6/91 (6.2%)<0.001
DNMT3A (+/-)47/147 (24.2%)40/57 (41.2%)7/90 (7.2%)<0.001
IDH2 (+/-)18/176 (9.3%)8/89 (8.2%)10/87 (10.3%)0.805
IDH1 (+/-)19/175 (9.8%)11/86 (11.3%)8/89 (8.3%)0.630
TET2 (+/-)16/178 (8.2%)10/87 (10.3%)6/91 (6.2%)0.435
RUNX1 (+/-)15/179 (8.2%)4/93 (4.1%)11/86 (11.3%)0.104
TP53 (+/-)16/178 (8.2%)7/90 (7.2%)9/88 (9.3%)0.795
NRAS (+/-)15/179 (8.2%)7/90 (7.2%)8/89 (8.2%)>0.999
CEBPA (+/-)13/181 (6.7%)5/92 (5.2%)8/89 (8.2%)0.568
WT1 (+/-)11/183 (5.7%)5/92 (5.2%)6/91 (6.2%)>0.999
PTPN11 (+/-)9/185 (4.6%)6/91 (6.2%)3/94 (3.1%)0.497
KIT (+/-)8/186 (4.1%)3/94 (3.1%)5/92 (5.2%)0.721
U2AF1 (+/-)8/186 (4.1%)3/94 (3.1%)5/92 (5.2%)0.721
KRAS (+/-)6/188 (3.1%)4/93 (4.1%)2/95 (2.1%)0.683

Correlation of clinic-pathologic characteristics in AML patients with HOXB3 hypo- and hyper-methylation.

AML, acute myeloid leukemia; WBC, white blood cells; PB, peripheral blood; BM, bone marrow; FAB, French-American-British classification. Cytogenetic and molecular risk classifications are based on the 2017 European LeukemiaNet (ELN) classification.

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 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

VariablesTotal AMLNon-M3 AML
Univariate analysisMultivariate analysisUnivariate analysisMultivariate analysis
HR (95% CI)PHR (95% CI)PHR (95% CI)PHR (95% CI)P
Overall survival
Age1.039 (1.026-1.053)<0.0011.026 (1.012-1.040)<0.0011.032 (1.018-1.046)<0.0011.020 (1.005-1.035)0.007
WBC1.003 (0.999-1.006)0.1191.002 (0.999-1.006)0.1501.001 (0.997-1.005)0.558
Molecular risks1.870 (1.474-2.371)<0.0012.246 (1.704-2.962)<0.0011.701 (1.303-2.222)<0.0012.173 (1.610-2.934)<0.001
Treatment regimen0.540 (0.380-0.767)0.0010.479 (0.319-0.720)<0.0010.446 (0.311-0.639)<0.0010.432 (0.286-0.654)<0.001
FLT3 mutation1.271 (0.870-1.858)0.2151.352 (0.915-2.000)0.1301.423 (0.932-2.174)0.102
NPM1 mutation1.220 (0.836-1.780)0.3021.063 (0.725-1.557)0.756
DNMT3A mutation1.615 (1.100-2.373)0.0151.153 (0.751-1.771)0.5141.431 (0.971-2.109)0.0701.147 (0.749-1.756)0.528
HOXB3 methylation0.613 (0.432-0.869)0.0060.664 (0.467-0.945)0.0230.620 (0.433-0.890)0.0090.656 (0.455-0.946)0.024
Disease-free survival
Age1.034 (1.021-1.047)<0.0011.024 (1.010-1.038)0.0011.027 (1.013-1.040)<0.0011.021 (1.007-1.036)0.004
WBC1.003 (1.000-1.006)0.0941.002 (0.999-1.006)0.1391.001 (0.998-1.005)0.422
Molecular risks1.796 (1.424-2.265)<0.0012.040 (1.563-2.661)<0.0011.625 (1.253-2.109)<0.0012.005 (1.503-2.675)<0.001
Treatment regimen0.607 (0.427-0.862)0.0050.565 (0.380-0.839)0.0050.522 (0.365-0.748)<0.0010.532 (0.357-0.792)0.002
FLT3 mutation1.252 (0.857-1.829)0.2451.343 (0.909-1.983)0.1391.442 (0.944-2.201)0.090
NPM1 mutation1.263 (0.865-1.843)0.2271.107 (0.755-1.622)0.604
DNMT3A mutation1.524 (1.040-2.235)0.0311.059 (0.691-1.624)0.7921.360 (0.924-2.001)0.1191.059 (0.692-1.620)0.793
HOXB3 methylation0.608 (0.429-0.862)0.0050.627 (1.563-2.661)0.0090.621 (0.433-0.890)0.0090.671 (0.460-0.980)0.039

Univariate and multivariate Cox regression analysis of potential risk factors for overall survival and leukemia free survival in AML patients.

HR, hazard ratio; CI, confidence interval; WBC, white blood cells. Variables including age (continuous variables), WBC (continuous variables), treatment regimen (with transplantation vs. without transplantation), molecular risks (good, intermediate, poor and unknown; classified by the 2017 European LeukemiaNet classification), FLT3/NPM1/DNMT3A mutation (wild type vs. mutant) and HOXB3 methylation (hyper- vs. hypo-). Multivariate analysis includes variables with P < 0.200 in univariate analysis.

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 (–, ). Moreover, positively correlated DEmiRs (top 5) such as miR-10, miR-196a and miR-1 were reported to have pro-leukemic effects (–), whereas negatively correlated DEmiRs (top 5) such as miR-193b and miR-379 were reported to have anti-leukemic effects in AML (, ). Moreover, Gene Ontology (GO) analysis demonstrated that these DEGs were involved in multiple biological processes such as development (Figure 3d).

Figure 3

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

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 (). 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 (). 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 (). 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 (). 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 (). Moreover, HOX overexpression may be correlated with a distinct biologic subset of AML with intermediate cytogenetics (). 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 (–). Moreover, several studies also addressed the closely associations of aberrant HOXB expression/methylation with NPM1, FLT3, and DNMT3A mutations in AML (, , ). 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 (). 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 (–). Consistently, Chen et al. found that HOXB5 is closely associated with the expression of other HOX family genes (), 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.

Statements

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

Summary

Keywords

AML, HOXB, HoxB3, hypomethylation, prognosis, regulation

Citation

Zhang T, Chang R, Xie F, Xu Z, Chu M, Wu X, Qian J and Zhou J (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

Revised

11 December 2025

Accepted

22 December 2025

Published

12 January 2026

Volume

16 - 2025

Edited by

Ananya Datta Mitra, UC Davis Health, United States

Reviewed by

Irma Olarte, Hospital General de México Dr. Eduardo Liceaga, Mexico

Alireza Ghezavati, Quest Diagnostics, United States

Updates

Copyright

*Correspondence: Jing-dong Zhou,

† These authors have contributed equally to this work

Disclaimer

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics