Abstract
Introduction
Acute myeloid leukemia (AML) is a blood cancer characterized by excessive growth of immature myeloid cells and makes up 15% to 20% of childhood acute leukemia. 1 Despite intensified chemotherapeutic regimens and enhanced supportive care measures, and even though allogeneic stem cell transplantation is considered effective for pediatric AML, the overall survival rate for children is still below 70%. 2 In addition, after initial remission, the recurrence rate remains at 15% to 20%. 3 Prognostic biomarkers have the potential to provide valuable insights into the possible outcomes of leukemia, assisting in patient stratification, targeted therapy, and disease monitoring in clinical practice, improving patient prognoses. 4 There is a strong need to identify biomarkers and develop novel biological therapeutic targets to increase pediatric AML survival rates and decrease treatment cytotoxicity.
MicroRNAs (miRNAs) are a class of small noncoding RNA molecules that typically range from 19 to 25 nucleotides in length; they regulate gene expression at the posttranscriptional stage, 5 play a vital role in AML tumorigenesis, and influence hematopoietic cell differentiation and proliferation. 6 It is becoming increasingly clear that miRNAs are involved in leukemia diagnosis, therapeutic response, and prognosis, suggesting their key role in the progression of leukemia.7–9 miR-133a is a muscle-specific miRNA; recent studies have shown that it plays an important role in organ fibrosis, inflammatory responses, and cancer development.10–12 Yang et al found that miRNA-133a is an independent prognostic marker in patients with osteosarcoma and is correlated with cumulative survival. 13 Flavia et al found that miR-133a was expressed at low levels in K562 cells, but was undetectable in HL-60 cells. 14 Interestingly, miR-133a-5p silenced TIM-3 expression, a known biomarker of AML leukemic stem cells. 15 miR-135a is associated with various types of cancer, playing a dual role as a tumor suppressor in some cases while exerting oncogenic functions in others. Low miR-135a levels are associated with a higher risk of AML relapse. 16
Our previous study found that induction of chemotherapeutics such as all-trans-retinoic acid (ATRA) could decrease the activity of homeobox transcription factor caudal type homeobox 2 (CDX2), leading to reduced transcription of miR-125b and enhanced expression of core binding factor β (CBFβ); this promotes myeloid cell differentiation, counteracting leukemogenesis. 17 However, the mechanisms by which ATRA regulates the abnormal expression of CDX2 remain unclear. Our previous miRNA profiling analysis showed that miR-133a and miR-135a were downregulated at diagnosis in pediatric patients with AML compared to healthy controls. 18 Previous studies have also found that ATRA promotes the expression of miR-133a and miR-135a.19,20 These findings raises the possibility that miR-133a and miR-135a might play an intermediary role between ATRA and CDX2 with involvement in the development or pathogenesis of pediatric AML.
In this study, we examined the roles of miR-133a and miR-135a in the ATRA/CDX2/miR-125b/CBFβ pathway and investigated the relationship between miR-133a/135a and complete survival in pediatric patients with AML.
Materials and Methods
Patients and Samples
This study enrolled 59 pediatric patients, including five AML subtypes (AML-M1 to AML-M5), along with 9 healthy controls. Written informed consent was obtained from the parents or guardians of all participants. Diagnoses were established based on clinical features and morphological criteria following the French (American) classification. The biological and clinical characteristics of the patients are presented in Table 1. Bone marrow samples were obtained between January 2007 and April 2014 at the First Affiliated Hospital of Sun Yat-Sen University. This study was approved by the Ethics Committee of the First Affiliated Hospital of Sun Yat-Sen University (No. [2030] 070). This study conformed to the REMARK guidelines. 21
Characteristics of Patients with AML.
AML: acute myeloid leukemia.
Note. aWith at least 3 years of clinical follow-up.
Cell Culture and RNA/Protein Isolation
NB4 and HL-60 (human myeloid leukemia cell lines) were cultured in RPMI-1640 (HyClone, UT, USA). Both cultures were supplemented with 10% fetal bovine serum (HyClone) and incubated in a 5% CO2 atmosphere at 37 °C. Cell lines were treated with 1 µM ATRA at a density of 1 × 106 cells per well. Total RNA and proteins were extracted using Tirol (Invitrogen) according to the manufacturer's instructions.
Quantitative Real-Time PCR Analysis
Quantitative reverse transcription PCR (qRT-PCR) was performed to detect mature miRNAs and CDX2 mRNA. A total of 68 bone marrow samples were used to determine the expression of miR-133a and miR-135a by qRT-PCR, including nine controls and 59 patients with AML (38 at diagnosis, 13 in CR, and 8 in relapse, all of which were not in pairs). Total RNA was quantified using SYBR Premix Ex Taq II (TaKaRa) and the ABI Step One Detection System (Applied Biosystems). Detection of mature miR-133a and miR-135a in total extracted RNA was performed using the Hairpin-itTM miRNA Real-time PCR Quantization Kit (GenePharma, Shanghai, China), which included a stem-loop-like RT primer and miRNA-specific primers. As an internal control, U6 small RNA was amplified using the same method. mRNA quantification was performed using endogenous β-actin mRNA as an internal control within the same system. Each assay was performed in triplicate.
Cell Transfection and Vector Constructs
The following materials were acquired from GenePharma (Shanghai, China): miR-133a, miR-135a mimics, and scrambled oligonucleotides (miR-NC; negative control); H2O (denoted as miR-mock or anti-mock); miR-133a and miR-135a antisense sequence complementary to mature miR-133a and miR-135a (anti-miR-133a/anti-miR-135a); and miR antisense NC (anti-NC; negative control for miR-133a/miR-135a antisense). HEK-293 T cells were cultured in 24-well plates at a density of 1 × 105 cells per well in a complete growth medium (0.5 mL) and transfected with mimics or miR-NC using lipofectamine 2000 (Invitrogen). NB4 cells were cultured using similar methods and electrotransfected with miRNA mimics or antisense miRNAs at a concentration of 100 nM. Transfected cells were harvested at 48 h posttransfection for western blotting. The fragments of the 3`-UTR of CDX2 mRNA containing the wild-type and mutant binding sites for miR-133a and miR-135a were synthesized by Songon (Shanghai, China) and then inserted into the psiCHECK-2 vector (Promega) for functional analysis. Wild-type and mutant insertions were confirmed using DNA sequencing.
Target Prediction and Luciferase Reporter Assay
To identify possible candidate miRNAs targeting CDX2, we employed PICTAR (http://pictar.mdc-berlin.de/) and TargetScan (http://www.targetscan.org/) as prediction algorithms. Transient transfections were performed in HEK-293 T cells, and cotransfection of CDX2 constructs (psi-Check2-CDX2-WT/Mut), psiCheck control, miR-133a or miR-135a, miR-133a, or miR-135a negative control, and the Renilla luciferase plasmid pRL-TK was performed using Lipofectamine 2000 (Invitrogen). A dual-luciferase reporter assay system (Promega) was used to assess luciferase activity. Relative activity was computed by normalizing the firefly luciferase activity as to the Renilla luciferase activity as an internal control.
Western Blot Assay
Cells were collected 48 h after transfections and lysed with protease and phosphatase inhibitors in RIPA buffer (Pierce, Rockford, IL, USA). The same amount of proteins (2 µg) and loading control (Fermentas) were separated by 10% SDS-PAGE and subsequently transferred onto a polyvinylidene fluoride membrane (Milli-pore, Beijing, China). The membranes were incubated with antibodies against CDX2 (1:5000, BioGenex, SanRamon, CA, USA) or β-tubulin (1:1000, Cell Signaling Technology, MA, USA) at 4 ˚C overnight. After 1 h of incubation with an HRP-conjugated secondary antibody (antirabbit, Sigma-Aldrich), an ECL Plus immunoblotting detection system (GE Healthcare Biosciences) was used to detect protein levels.
Cell Proliferation Assays
A Cell Counting Kit-8 (CCK-8) was used to perform a cell proliferation assay (Dojindo Molecular Technologies, Shanghai, China). Cells were transfected by electroporation with 100 nM miR-133a, miR-135a mimics, 200 nM anti-miR-133a, and anti-miR-135a antisense along with the relative doses of miR-NC and anti-NC. After transfection, 1 × 104 cells per well were seeded into 96-well plates and the CCK-8 assay was performed at 24, 48, 72, and 96 h to measure cell proliferation.
Statistical Analysis
Statistical analyses were conducted using SPSS (Version 15.0; Chicago, IL, USA). Luciferase and qRT-PCR data are presented as the mean ± standard deviation (sd) from three independent experiments. Unless otherwise indicated, differences among groups were assessed using Student's t-test. Statistical significance was set at P < .05.
Results
Differing Expression of CDX2 in ALL and AML, Subtypes of AML, and Primary Cells From Patients
According to the database, CDX2 was more highly expressed in AML than in ALL, particularly in M6 cells (Figure 1A-1B). To confirm that CDX2 acts as an oncogene in AML, we examined the expression of CDX2 in the primary cells of patients. As shown in Figure 1C, CDX2 was expressed at low levels in patients in complete remission compared with samples taken at diagnosis and relapse.

Differences in expression of CDX2 in ALL and AML, subtypes of AML, and primary cells from patients. (A) Comparison of CDX2 expression between patients with ALL and patients with AML according to the database. (B) CDX2 was highly expressed in AML-M6 compared with other AML subtypes. (C) The expression level of CDX2 in different disease stages in pediatric patients with AML was quantified by qRT-PCR. CDX2 expression was lower in CR patients compared with patients at diagnosis or with relapse. (*P < .05, **P < .01, ***P < .001).
miR-133a and miR-135a are Predicted to Interact with CDX2
To identify potential molecules targeting CDX2 involved in ATRA-mediated differentiation, we searched miRNA databases for several candidates (Figure 2A-2B). Combining previous research and our experimental results (Figure 2C), miR-133a and miR-135a were computationally predicted to possibly interact with a presumed binding site located in the 3'-UTR of CDX2 (Figure 2D-2E). Luciferase reporter assays were performed to confirm these findings.

miR-133a and miR-135a are predicted to interact with CDX2. (A/B) TargetScan and PICTAR were used to search for candidate miRNAs. (C) HL-60 cells were treated with ATRA for 48 h, and only miR-133a and miR-135a were significantly expressed. (D/E) miR-133a and miR-135a can target the CDX2 3′UTR. Sequences of the miR-133a and miR-135a binding sites in human CDX2 genomic regions. (F/G) Luciferase reporter assays were employed to confirm the interaction of miR-133a and miR-135a and the CDX2 3′ UTR. Firefly luciferase activity was normalized to Renilla luciferase activity, and the obtained results were presented with reference to the control (NC).
As shown in Figure 2F-2G, when the wild-type of CDX2 was cotransfected with miR-133a and miR-135a mimics into HEK-293 T cells, both miRNAs significantly repressed luciferase activity compared to miR-NC by 42.7% and 44.1%, respectively. However, with mutant CDX2, the deletion of miR-133a and miR-135a seed recognition sequences did not reduce luciferase activity.
miR-133a and miR-135a Downregulate CDX2 Expression Through Translational Inhibition
Further verification experiments were performed to determine the potential effects of miR-133a and miR-135a on CDX2 expression. Overexpression and knockdown of miR-133a and miR-135a in NB4 and HL-60 cells, respectively, indicated that the expression of CDX2 protein was significantly reduced in miR-133a- and miR-135a-transfected cells compared with cells transfected with miR-NC and miR-mock (Figure 3B). Transfecting cells with anti-miR-133a and anti-miR-135a upregulated CDX2 protein expression in comparison with anti-NC or anti-mock (Figure 3A).

miR-133a and miR-135a downregulate CDX2 expression through translational repression. (A/B) The expression level of CDX2 was analyzed by western blotting after the overexpression of miR-133a and miR-135a mimics or the knockdown of endogenous miR-133a and miR-135a using miR-133a and miR-135a antisense in NB4 and HL-60 cells.
miR-133a and miR-135a Inhibit the Proliferation of Two AML Cell Lines
To further demonstrate the potential role of miR-133a and miR-135a in myeloid cell proliferation, we performed a CCK-8 assay on the NB4 and HL-60 cell lines when miR-133a and miR-135a were overexpressed or downregulated, respectively. The results indicated that the survival rates of NB4 and HL-60 cells transfected with anti-133a and anti-135a antibodies were significantly increased (Figure 4A-4B). Cell survival rates in the miR-133a and miR-135a groups were significantly lower than those in the anti-NC and mock groups (Figure 4C-4D).

miR-133a and miR-135a inhibit the proliferation of two AML cell lines. (A/B) After transfection with anti-133a and anti-135a, the proliferation of NB4 and HL-60 cells was significantly promoted. (C/D) Overexpression of miR-133a and miR-135a inhibited proliferation in NB4 and HL-60 cells.
miR-133a and miR-135a Are Associated With Different Disease Stages, Subtypes, and Survival
We also investigated the expression of miR-133a and miR-135a at different disease stages and in different subtypes. As shown in Figure 5A-5B, miR-133a and miR-135a were more highly expressed in patients in CR than in patients at diagnosis or with relapse. Figure 5C-5D further demonstrates the significant differences between CR and non-CR. In addition, the two miRNAs had high expression levels in AML subtypes M1-3. Lower expression was observed in other subtypes, including AML M4-5 (Figure 5E-5F).

miR-133a and miR-135a are correlated with different pediatric AML disease stages, subtypes, and survival. (A/B) The expression of miR-133a and miR-135a was measured by RT-PCR in nine controls and 59 patients with AML. Both were significantly upregulated in CR patients compared to diagnosis and relapse. (C/D) Comparison of the expression of miR-133a and miR-135a between CR and non-CR patients by qRT-PCR. (E/F) MiR-133a and miR-135a were highly expressed in M1-M3 patients in comparison with other subtypes of pediatric AML. (G/H) Patients with elevated expression levels of miR-133a and miR-135a exhibited significantly improved leukemia-free survival (LFS) compared to those with lower expression. (I/J) The expression of miR-133a and miR-135a was not associated with gender.
A positive correlation was found between elevated expression of miR-133a and miR-135a and improved leukemia-free survival rates during the 5-year follow-up period (Figure 5G-5H). There was no significant difference in expression between male and female mice (Figure 5I-5J).
Discussion
Relapse and failure to respond to first line treatments remain a barrier to survival in children suffering from AML. 22 It is of great importance to clarify the mechanism of leukemogenesis in AML and identify prognostic factors to design novel targeted therapies, predict clinical outcomes, and guide risk-adjusted clinical management.23,24 In recent years, considerable evidence has emerged indicating that miRNAs serve as both oncogenes or tumor suppressors, acting as potential novel biomarkers for leukemia diagnosis, treatment, and prognosis determination.25,26 Our previous studies demonstrated aberrant expression of miR-133a and miR-135a in leukemia cells 18 ; however, few studies have investigated the expression and function of miR-133a and miR-135a in pediatric AML.
The caudal-type homeobox genes are homologous to the caudal-type genes of Drosophila and include CDX1, CDX2, and CDX4. CDX2 is located on chromosome 13q12-13 and consists of three exons and two introns, encoding a protein of 311 amino acids. 27 The CDX2 protein binds to a corresponding DNA region through a helix-loop helix and regulates DNA expression as a transcription factor. 28 CDX2 is primarily recognized as a developmental regulator of early embryogenesis and hematopoiesis.29–31 In normal hematopoietic stem cells (HSC), CDX2 is not expressed; however, it is highly expressed in patients with AML, advanced chronic myeloid leukemia, and high-risk myelodysplastic syndrome.32–34 Aberrant expression of CDX2 has been implicated in promoting the transformation of HSCs into leukemia stem cells, and it is believed that CDX2 might be necessary for leukemia growth. 35 In the present study, we observed aberrant CDX2 expression in pediatric AML, similar to that reported in ectopic expression in adult leukemia. We also identified differential CDX2 expression patterns in different disease stages and subtypes of pediatric AML.
MiR-133a is a confirmed tumor suppressor that has low expression levels in various cancers, such as ovarian, colorectal, bladder, breast, prostate, and esophageal cancers.36–43 miR-135a has also been shown to be a novel tumor regulator; in classical Hodgkin's lymphoma and gastric cancer, miR-135a inhibits tumor cell proliferation by targeting JAK2 to suppress STAT3 activation, decreasing cyclin D1 and Bcl-xL expression.44,45 Additionally, in vitro experiments demonstrated that miR-135a exerted a substantial inhibitory effect on the migratory and invasive capabilities of lung cancer cells. 46 However, in bladder cancer, miR-135a promotes cell proliferation by targeting PHLPP2 and FOXO1, which function as onco-miR. 47 In this study, we found that miR-133a and miR-135a exhibited promising anticancer effects in pediatric AML by directly targeting CDX2, indicating that these miRNAs may be novel therapeutic targets for pediatric AML. We established a correlation between the ATRA-miRNA-CDX2 pathway for the first time, providing fundamental theoretical support for further clarifying the pathogenesis of leukemia and elucidating the mechanism of drug action.
Homeobox (HOX) transcript antisense RNA (HOTAIR) is a lncRNA located on chromosome 12. It can form competitive endogenous RNAs with miRNAs and participate in various biological processes, including proliferation, differentiation, and apoptosis. 48 Liu et al reported that HOTAIR regulates HER2 expression by sponging miR-331-3p in gastric cancer. 49 In addition, Nashwa et al found a significant difference in HOTAIR, c-kit, and miR-193a expression between patients with AML that had 12p abnormalities and those without. Patients with low HOTAIR and c-kit expression exhibited significantly better survival. 50 In contrast, Hu et al found that HOTAIR was significantly downregulated in patients with AML compared to the iron-deficiency anemia control group and that HOTAIR promoted ATRA-induced differentiation by regulating the cell cycle regulator p21 through miR-17-5p. 51 In addition, CDX2 was identified as a transcription factor for HOTAIR in hepatocellular carcinoma by using ChIPBase. 52 Considering this, the potential ATRA/HOTAIR/miRNA/CDX2 pathway deserves further exploration.
We also investigated the relationship between miR-133a and miR-135a expression in different disease stages, subtypes, and survival rates. The results showed that miR-133a and miR-135a were more highly expressed in patients in complete remission than those at diagnosis and with relapse, as well as in patients with better leukemia-free survival. Together, these findings collectively indicate that increased expression levels of miR-133a and miR-135a are associated with a more favorable prognosis and enhanced survival outcomes in pediatric AML. However, this study has some limitations. The mechanism by which ATRA enhances miRNA133a/135a expression remains unclear, and further exploration is required to elucidate the mechanisms underlying leukemogenesis as well as the development of drug resistance and relapse. Additionally, considerable efforts are necessary to discover novel therapeutics, enhance prognosis, and improve overall quality of life.
Conclusion
In conclusion, our study showed that ATRA promoted the expression of miR-133a and miR-135a in AML cells, inhibiting CDX2 expression to inhibit AML cell proliferation. Additionally, high expression of miR-133a and miR-135a were commonly seen in patients in CR and better survival. These findings indicate that miR-133a and miR-135a may serve as potential prognostic markers and novel therapeutic targets in pediatric AML.
Footnotes
Authors’ Note
Yu-Cai Cheng, Zhong Fan, and Cong Liang, contributed equally to this article.
Data Availability
The datasets generated during the study are available from the corresponding author.
Declaration of Conflicting Interests
The authors declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Ethical Statement
The study was approved by the Ethics Committee of the First Affiliated Hospital of Sun Yat-Sen University. The patients and guardians have signed the informed consent form.
Funding
The authors disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by the Medical Science and Technology Research Foundation of Guangdong Province, China, the Shenzhen Fundamental Research Program (Grant Nos. A2021328, A2022443, and JCYJ20190809142619278).
