HOXA9 promotes MYC-mediated leukemogenesis by maintaining gene expression for multiple anti-apoptotic pathways

Bahcc1 is critical for the aberrant epigenetic program in a mouse model of MLL-ENL-mediated leukemia

ABSTRACT

In leukemogenesis, genotoxic stress in hematopoietic stem and progenitor cells (HSPCs) drives individual context-dependent programs of malignant transformation. In light of the various differentiation stages of HSPCs based on a recently revised definition using CD150/CD48, our analyses showed that a subpopulation of long-term repopulating HSCs was most susceptible to MLL-ENL-mediated transformation. An analysis of the molecular mechanism identified Bromo-adjacent homology domain and coiled-coil containing 1 (Bahcc1), which encodes a reader molecule of trimethylated histone H3 lysine 27 (H3K27me3), as a candidate gene involved in distinct susceptibility to leukemic transformation. Interestingly, Bahcc1 was previously reported to be highly expressed in acute myeloid leukemia (AML) with an unfavorable prognosis, including some cases of MLL-rearranged AML. We found that MLL-ENL upregulated Bahcc1 through binding to its promoter, and that Bahcc1 was involved in MLL-ENL-mediated immortalization at least partly through repression of H3K27me3-marked Cdkn1c. Analyses using bone marrow transplantation in mice showed that depletion of Bahcc1 suppressed the leukemogenic activity of MLL-ENL. In a public database, high BAHCC1 expression was found to be associated with a poor prognosis in pediatric AML, in which BAHCC1 expression was significantly lower in MLL-AF9-AML than in other MLL-fusion-AML. These findings shed light on the distinct immortalization potential of HSPCs and suggest a novel MLL-fusion-Bahcc1 axis, which may lead to development of molecular targeted therapy against MLL-fusion-mediated leukemia.

Exploring the expression and clinical significance of the miR-140-3p-HOXA9 axis in colorectal cancer

PURPOSE

This study aims to investigate the expression patterns and clinical significance of miR-140-3p and homeobox A9 (HOXA9) in colorectal cancer (CRC) selected by bioinformatic study, while elucidating their potential interplay.

METHODS

The microRNA expression profiles of paired colorectal cancer and matched normal tissues were retrieved from the Gene Expression Omnibus Database. Differentially expressed microRNAs and microRNA candidates were filtered and subjected to further analysis. Clinicopathological data, along with paraffin-embedded samples of colorectal tumor tissues were collected to facilitate comprehensive analysis. Expression levels of miR-140-3p and HOXA9 were quantified using qRT-PCR and immunohistochemistry. Survival rates were determined using the Kaplan-Meier method, and the COX regression model was utilized to identify independent prognostic factors that impact the overall prognosis.

RESULTS

MiR-140-3p was significantly downregulated in colorectal tumors compared to normal tissue, and HOXA9 was identified as a previously unreported potential downstream target. HOXA9 expression was elevated in tumors compared to normal tissues. Reduced miR-140-3p expression was associated with lymph node metastasis, while high HOXA9 expression correlated with both lymph node metastasis and lympho-vascular invasion. Patients with low miR-140-3p and high HOXA9 expression had a poorer prognosis. HOXA9 was identified as an independent risk factor for CRC patient survival.

CONCLUSION

The miR-140-3p-HOXA9 signaling disruption is closely linked to lymph node metastasis and unfavorable prognosis in CRC. This axis shows promise as a clinical biomarker for predicting the CRC patient survival and a potential therapeutic target.

Epigenomic mapping reveals distinct B cell acute lymphoblastic leukemia chromatin architectures and regulators

B cell lineage acute lymphoblastic leukemia (B-ALL) is composed of diverse molecular subtypes, and while transcriptional and DNA methylation profiling has been extensively examined, the chromatin landscape is not well characterized for many subtypes. We therefore mapped chromatin accessibility using ATAC-seq in primary B-ALL cells from 156 patients spanning ten molecular subtypes and present this dataset as a resource. Differential chromatin accessibility and transcription factor (TF) footprint profiling were employed and identified B-ALL cell of origin, TF-target gene interactions enriched in B-ALL, and key TFs associated with accessible chromatin sites preferentially active in B-ALL. We further identified over 20% of accessible chromatin sites exhibiting strong subtype enrichment and candidate TFs that maintain subtype-specific chromatin architectures. Over 9,000 genetic variants were uncovered, contributing to variability in chromatin accessibility among patient samples. Our data suggest that distinct chromatin architectures are driven by diverse TFs and inherited genetic variants that promote unique gene-regulatory networks.

HOXA9 transcription factor is a double-edged sword: from development to cancer progression

The HOXA9 transcription factor serves as a molecular orchestrator in cancer stemness, epithelial-mesenchymal transition (EMT), metastasis, and generation of the tumor microenvironment in hematological and solid malignancies. However, the multiple modes of regulation, multifaceted functions, and context-dependent interactions responsible for the dual role of HOXA9 as an oncogene or tumor suppressor in cancer remain obscure. Hence, unravelling its molecular complexities, binding partners, and interacting signaling molecules enables us to comprehend HOXA9-mediated transcriptional programs and molecular crosstalk. However, it is imperative to understand its central role in fundamental biological processes such as embryogenesis, foetus implantation, hematopoiesis, endothelial cell proliferation, and tissue homeostasis before designing targeted therapies. Indeed, it presents an enormous challenge for clinicians to selectively target its oncogenic functions or restore tumor-suppressive role without altering normal cellular functions. In addition to its implications in cancer, the present review also focuses on the clinical applications of HOXA9 in recurrence and drug resistance, which may provide a broader understanding beyond oncology, open new avenues for clinicians for accurate diagnoses, and develop personalized treatment strategies. Furthermore, we have also discussed the existing therapeutic options and accompanying challenges in HOXA9-targeted therapies in different cancer types.

Graphical Abstract

SET-PP2A complex as a new therapeutic target in KMT2A (MLL) rearranged AML

KMT2A-rearranged (KMT2A-R) is an aggressive and chemo-refractory acute leukemia which mostly affects children. Transcriptomics-based characterization and chemical interrogation identified kinases as key drivers of survival and drug resistance in KMT2A-R leukemia. In contrast, the contribution and regulation of phosphatases is unknown. In this study we uncover the essential role and underlying mechanisms of SET, the endogenous inhibitor of Ser/Thr phosphatase PP2A, in KMT2A-R-leukemia. Investigation of SET expression in acute myeloid leukemia (AML) samples demonstrated that SET is overexpressed, and elevated expression of SET is correlated with poor prognosis and with the expression of MEIS and HOXA genes in AML patients. Silencing SET specifically abolished the clonogenic ability of KMT2A-R leukemic cells and the transcription of KMT2A targets genes HOXA9 and HOXA10. Subsequent mechanistic investigations showed that SET interacts with both KMT2A wild type and fusion proteins, and it is recruited to the HOXA10 promoter. Pharmacological inhibition of SET by FTY720 disrupted SET-PP2A interaction leading to cell cycle arrest and increased sensitivity to chemotherapy in KMT2A-R-leukemic models. Phospho-proteomic analyses revealed that FTY720 reduced the activity of kinases regulated by PP2A, including ERK1, GSK3β, AURB and PLK1 and led to suppression of MYC, supporting the hypothesis of a feedback loop among PP2A, AURB, PLK1, MYC, and SET. Our findings illustrate that SET is a novel player in KMT2A-R leukemia and they provide evidence that SET antagonism could serve as a novel strategy to treat this aggressive leukemia.

Epigenomic mapping in B-cell acute lymphoblastic leukemia identifies transcriptional regulators and noncoding variants promoting distinct chromatin architectures

B-cell lineage acute lymphoblastic leukemia (B-ALL) is comprised of diverse molecular subtypes and while transcriptional and DNA methylation profiling of B-ALL subtypes has been extensively examined, the accompanying chromatin landscape is not well characterized for many subtypes. We therefore mapped chromatin accessibility using ATAC-seq for 10 B-ALL molecular subtypes in primary ALL cells from 154 patients. Comparisons with B-cell progenitors identified candidate B-ALL cell-of-origin and AP-1-associated cis-regulatory rewiring in B-ALL. Cis-regulatory rewiring promoted B-ALL-specific gene regulatory networks impacting oncogenic signaling pathways that perturb normal B-cell development. We also identified that over 20% of B-ALL accessible chromatin sites exhibit strong subtype enrichment, with transcription factor (TF) footprint profiling identifying candidate TFs that maintain subtype-specific chromatin architectures. Over 9000 inherited genetic variants were further uncovered that contribute to variability in chromatin accessibility among individual patient samples. Overall, our data suggest that distinct chromatin architectures are driven by diverse TFs and inherited genetic variants which promote unique gene regulatory networks that contribute to transcriptional differences among B-ALL subtypes.

Endemic Burkitt lymphoma avatar mouse models for exploring inter-patient tumor variation and testing targeted therapies

Endemic Burkitt lymphoma (BL) is a childhood cancer in sub-Saharan Africa characterized by Epstein-Barr virus and malaria-associated aberrant B-cell activation and MYC chromosomal translocation. Survival rates hover at 50% after conventional chemotherapies; therefore, clinically relevant models are necessary to test additional therapies. Hence, we established five patient-derived BL tumor cell lines and corresponding NSG-BL avatar mouse models. Transcriptomics confirmed that our BL lines maintained fidelity from patient tumors to NSG-BL tumors. However, we found significant variation in tumor growth and survival among NSG-BL avatars and in Epstein-Barr virus protein expression patterns. We tested rituximab responsiveness and found one NSG-BL model exhibiting direct sensitivity, characterized by apoptotic gene expression counterbalanced by unfolded protein response and mTOR pro-survival pathways. In rituximab-unresponsive tumors, we observed an IFN-α signature confirmed by the expression of IRF7 and ISG15. Our results demonstrate significant inter-patient tumor variation and heterogeneity, and that contemporary patient-derived BL cell lines and NSG-BL avatars are feasible tools to guide new therapeutic strategies and improve outcomes for these children.

A pro B cell population forms the apex of the leukemic hierarchy in Hoxa9/Meis1-dependent AML

Relapse is a major challenge to therapeutic success in acute myeloid leukemia (AML) and can be partly associated with heterogeneous leukemic stem cell (LSC) properties. In the murine Hoxa9/Meis1-dependent (H9M) AML model, LSC potential lies in three defined immunophenotypes, including Lin^-cKit⁺ progenitor cells (Lin^-), Gr1⁺CD11b⁺cKit⁺ myeloid cells, and lymphoid cells (Lym⁺). Previous reports demonstrated their interconversion and distinct drug sensitivities. In contrast, we here show that H9M AML is hierarchically organized. We, therefore, tracked the developmental potential of LSC phenotypes. This unexpectedly revealed a substantial fraction of Lin^- LSCs that failed to regenerate Lym⁺ LSCs, and that harbored reduced leukemogenic potential. However, Lin^- LSCs capable of producing Lym⁺ LSCs as well as Lym⁺ LSCs triggered rapid disease development suggestive of their high relapse-driving potential. Transcriptional analyses revealed that B lymphoid master regulators, including Sox4 and Bach2, correlated with Lym⁺ LSC development and presumably aggressive disease. Lentiviral overexpression of Sox4 and Bach2 induced dedifferentiation of H9M cells towards a lineage-negative state in vitro as the first step of lineage conversion. This work suggests that the potency to initiate a partial B lymphoid primed transcriptional program as present in infant AML correlates with aggressive disease and governs the H9M LSC hierarchy.

Role of HOXA9 in solid tumors: mechanistic insights and therapeutic potential

HOXA9 functioning as a transcription factor is one of the members of HOX gene family, which governs multiple cellular activities by facilitating cellular signal transduction. In addition to be a driver in AML which has been widely studied, the role of HOXA9 in solid tumor progression has also received increasing attention in recent years, where the aberrant expression of HOXA9 is closely associated with the prognosis of patient. This review details the signaling pathways, binding partners, post-transcriptional regulation of HOXA9, and possible inhibitors of HOXA9 in solid tumors, which provides a reference basis for further study on the role of HOXA9 in solid tumors.

Methylome Profiling of PD-L1-Expressing Glioblastomas Shows Enrichment of Post-Transcriptional and RNA-Associated Gene Regulation

Glioblastomas are the most frequent primary brain tumors in adults. They show highly malignant behavior and devastating outcomes. Since there are still no targeted therapies available, median survival remains in the range of 12 to 15 months for glioblastoma patients. Programmed Cell Death Ligand 1 (PD-L1) is a promising novel candidate in precision medicine. Here, we performed integrated epigenome-wide methylation profiling of 866,895 methylation-specific sites in 20 glioblastoma samples comparing PD-L1 high- (i.e., TPS (tumor proportion score) > 30%) and PD-L1 low-expressing glioblastomas (i.e., TPS < 10%). We found 12,597 significantly differentially methylated CpGs (DMCG) (Δβ ≥ 0.1 and p-value < 0.05) in PD-L1 high- compared with PD-L1 low-expressing glioblastomas. These DMCGs were annotated to 2546 tiling regions, 139 promoters, 107 genes, and 107 CpG islands. PD-L1 high-expressing glioblastomas showed hypomethylation in 68% of all DMCGs. Interestingly, the list of the top 100 significantly differentially methylated genes showed the enrichment of regulatory RNAs with 19 DMCGs in miRNA, snoRNAs, lincRNAs, and asRNAs. Gene Ontology analysis showed the enrichment of post-transcriptional and RNA-associated pathways in the hypermethylated gene regions. In summary, dissecting the methylomes depending on PD-L1 status revealed significant alterations in RNA regulation and novel molecular targets in glioblastomas.

Comprehensive analysis of a homeobox family gene signature in clear cell renal cell carcinoma with regard to prognosis and immune significance

The homeobox (HOX) family genes have been linked to multiple types of tumors, while their effect on malignant behaviors of clear cell renal cell carcinoma (ccRCC) and clinical significance remains largely unknown. Here, we comprehensively analyzed the expression profiles and prognostic value of HOX genes in ccRCC using datasets from The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC) databases. We developed a prognostic signature comprising eight HOX genes (HOXB1, HOXA7, HOXB5, HOXD8, HOXD9, HOXB9, HOXA9, and HOXA11) for overall survival prediction in ccRCC and it allowed patients to be subdivided into high- and low-risk groups. Kaplan-Meier survival analysis in all the internal and external cohorts revealed significant difference in clinical outcome of patients in different risk groups, indicating the satisfactory predictive power of the signature. Additionally, we constructed a prognostic nomogram by integrating signature-derived risk score and clinical factors such as gender, age, T and M status, which might be helpful for clinical decision-making and designing tailored management schedules. Immunological analysis revealed that the regulatory T cells (Tregs) infiltrated differently between the two subgroups in both TCGA and ICGC cohorts. ssGSEA method showed that the enrichment scores for mast cells were significantly lower in high-risk group compared with the low-risk group, which was consistent in both TCGA and ICGC cohorts. As for the related immune function, the enrichment scores of APC co-inhibition, para-inflammation, and type II IFN response were consistently lower in high-risk group in both cohorts. Of the eight HOX genes, the mRNA and protein levels of HOXD8 were downregulated in ccRCC than that in normal tissues, and decreased expression of HOXD8 was associated with increased tumor grade and stage, and lymph node metastasis. Survival analysis revealed that lower expression of HOXD8 predicted worse overall survival in ccRCC. In conclusion, our HOX gene-based signature was a favorable indicator to predict the prognosis of ccRCC cases and associated with immune cell infiltration. HOXD8 might be a tumor suppressor gene in ccRCC and a potential predictor of tumor progression.

Downregulation of hsa-miR-548d-3p and overexpression of HOXA9 in diffuse large B-cell lymphoma patients and the risk of R-CHOP chemotherapy resistance and disease progression

INTRODUCTION

Routine categorization of DLBCL patients into GCB and non-GCB groups by Hans' criteria could not accurately predict chemotherapy resistance and disease progression in patients treated with standard R-CHOP therapy. There is a need to identify better biomarker predictors to enhance assisted selection of chemotherapy regimens for DLBCL patients.

AIM OF THE STUDY

To identify dysregulated miRNAs and mRNAs that are predictive of resistance to R-CHOP chemotherapy or disease progression in patients with DLBCL.

METHODS

miRNA and mRNA profiling were performed on archival FFPE samples of the DLBCL patients. miRabel and miRNet bioinformatic tools were applied to determine experimental validated miRNA-mRNA target interaction. The significance of the genomic predictive values was assessed using adjusted odds ratios (AOR) and 95% confidence intervals (CI).

RESULTS

19/36 were R-CHOP therapy-resistant whilst 17/36 were R-CHOP therapy-sensitive. Ten dysregulated miRNAs and 12 dysregulated mRNAs were identified in therapy-resistant DLBCL patients. These dysregulated miRNAs and mRNA cause therapy resistance and disease progression in DLBCL patients, most likely via upregulation of the anti-apoptotic protein bcl2, activation of the JAK/STAT signalling pathway and dysregulation of p53 pathway. Downregulation of hsa-miR-548d-3p and overexpression of HOXA9 mRNA were significantly associated with therapy resistance and disease progression in DLBCL patients [hsa-miR-548d-3p AOR: 0.258, 95%CI: 0.097-0.684, p = 0.006].

CONCLUSION

DLBCL patients with downregulation of hsa-miR-548d-3p and overexpression of HOXA9 mRNA are more likely to experience R-CHOP therapy resistance and disease progression.