Underexpression of HOXA11 Is Associated with Treatment Resistance and Poor Prognosis in Glioblastoma

An atlas of causal and mechanistic drivers of interpatient heterogeneity in glioma

Grade IV glioma, formerly known as glioblastoma multiforme (GBM) is the most aggressive and lethal type of brain tumor, and its treatment remains challenging in part due to extensive interpatient heterogeneity in disease driving mechanisms and lack of prognostic and predictive biomarkers. Using mechanistic inference of node-edge relationship (MINER), we have analyzed multiomics profiles from 516 patients and constructed an atlas of causal and mechanistic drivers of interpatient heterogeneity in GBM (gbmMINER). The atlas has delineated how 30 driver mutations act in a combinatorial scheme to causally influence a network of regulators (306 transcription factors and 73 miRNAs) of 179 transcriptional "programs", influencing disease progression in patients across 23 disease states. Through extensive testing on independent patient cohorts, we share evidence that a machine learning model trained on activity profiles of programs within gbmMINER significantly augments risk stratification, identifying patients who are super-responders to standard of care and those that would benefit from 2 nd line treatments. In addition to providing mechanistic hypotheses regarding disease prognosis, the activity of programs containing targets of 2 nd line treatments accurately predicted efficacy of 28 drugs in killing glioma stem-like cells from 43 patients. Our findings demonstrate that interpatient heterogeneity manifests from differential activities of transcriptional programs, providing actionable strategies for mechanistically characterizing GBM from a systems perspective and developing better prognostic and predictive biomarkers for personalized medicine.

Intratumoral IL-12 delivery via mesenchymal stem cells combined with PD-1 blockade leads to long-term antitumor immunity in a mouse glioblastoma model

BACKGROUND

Although PD-1 blockade is effective for treating several types of cancer, the efficacy of this agent in glioblastoma is largely limited. To overcome non-responders and the immunosuppressive tumor microenvironment, combinational immunotherapeutic strategies with anti-PD-1 need to be considered. Here, we developed IL-12-secreting mesenchymal stem cells (MSC_IL-12) with glioblastoma tropism and evaluated the therapeutic effects of anti-PD-1, MSC_IL-12, and their combination against glioblastoma.

METHODS

Therapeutic responses were evaluated using an immunocompetent mouse orthotopic model. Tumor-infiltrating lymphocytes (TILs) were analyzed using immunofluorescent imaging. Single-cell transcriptome was obtained from mouse brains after treatments.

RESULTS

Anti-PD-1 and MSC_IL-12 showed complete tumor remission in 25.0% (4/16) and 23.1% (3/13) of glioblastoma-implanted mice, respectively, and their combination yielded synergistic antitumor efficacy indicated by 50.0% (6/12) of complete tumor remission. Analyses of TILs revealed that anti-PD-1 increased CD8+ T cells, while MSC_IL-12 led to infiltration of CD4+ T cells and NK cells. Both therapies reduced frequencies of Tregs. All these aspects observed in each monotherapy group were superimposed in the combination group. Notably, no tumor growth was observed upon rechallenge in cured mice, indicating long-term immunity against glioblastoma provoked by the therapies. Single-cell RNA-seq data confirmed these results and revealed that the combined treatment led to immune-favorable tumor microenvironment-CD4+, CD8+ T cells, effector memory T cells, and activated microglia were increased, whereas exhausted T cells, Tregs, and M2 polarized microglia were reduced.

CONCLUSION

Anti-PD-1 and MSC_IL-12 monotherapies show long-term therapeutic responses, and their combination further enhances antitumor efficacy against glioblastoma via inducing immune-favorable tumor microenvironment.

Development and validation of a predictive model in diagnosis and prognosis of primary glioblastoma patients based on Homeobox A family

BACKGROUND

Homeobox A (HOXA) family is involved in the development of malignancies as either tumor suppressors or oncogenes. However, their roles in glioblastoma (GBM) and clinical significance have not been fully elucidated.

METHODS

HOXA mutation and expressions in pan-cancers were investigated using GSCA and Oncomine, which in GBM were validated by cBioPortal, Chinese Glioma Genome Atlas (CGGA), and The Cancer Genome Atlas (TCGA) datasets. Kaplan-Meier analyses were conducted to determine prognostic values of HOXAs at genetic and mRNA levels. Diagnostic roles of HOXAs in tumor classification were explored by GlioVis and R software. Independent prognostic HOXAs were identified using Cox survival analyses, the least absolute shrinkage and selection operator (LASSO) regression, quantitative real-time PCR, and immunohistochemical staining. A HOXAs-based nomogram survival prediction model was developed and evaluated using Kaplan-Meier analysis, time-dependent Area Under Curve, calibration plots, and Decision Curve Analysis in training and validation cohorts.

RESULTS

HOXAs were highly mutated and overexpressed in pan-cancers, especially in CGGA and TCGA GBM datasets. Genetic alteration and mRNA expression of HOXAs were both found to be prognostic. Specific HOXAs could distinguish IDH mutation (HOXA1-7, HOXA9, HOXA13) and molecular GBM subtypes (HOXA1-2, HOXA9-11, HOXA13). HOXA1/2/3/10 were confirmed to be independent prognostic members, with high expressions validated in clinical GBM tissues. The HOXAs-based nomogram model exhibited good prediction performance and net benefits for patients in training and validation cohorts.

CONCLUSION

HOXA family has diagnostic values, and the HOXAs-based nomogram model is effective in survival prediction, providing a novel approach to support the treatment of GBM patients.

Molecular implications of HOX genes targeting multiple signaling pathways in cancer

Homeobox (HOX) genes encode highly conserved homeotic transcription factors that play a crucial role in organogenesis and tissue homeostasis. Their deregulation impacts the function of several regulatory molecules contributing to tumor initiation and progression. A functional bridge exists between altered gene expression of individual HOX genes and tumorigenesis. This review focuses on how deregulation in the HOX-associated signaling pathways contributes to the metastatic progression in cancer. We discuss their functional significance, clinical implications and ascertain their role as a diagnostic and prognostic biomarker in the various cancer types. Besides, the mechanism of understanding the theoretical underpinning that affects HOX-mediated therapy resistance in cancers has been outlined. The knowledge gained shall pave the way for newer insights into the treatment of cancer.

HOXA11-AS aggravates microglia-induced neuroinflammation after traumatic brain injury

Long noncoding RNAs (lncRNAs) participate in many pathophysiological processes after traumatic brain injury by mediating neuroinflammation and apoptosis. Homeobox A11 antisense RNA (HOXA11-AS) is a member of the lncRNA family that has been reported to participate in many inflammatory reactions; however, its role in traumatic brain injury remains unclear. In this study, we established rat models of traumatic brain injury using a weight-drop hitting device and injected LV-HOXA11-AS into the right lateral ventricle 2 weeks before modeling. The results revealed that overexpression of HOXA11-AS aggravated neurological deficits in traumatic brain injury rats, increased brain edema and apoptosis, promoted the secretion of proinflammatory factors interleukin-1β, interleukin-6, and tumor necrosis factor α, and promoted the activation of astrocytes and microglia. Microglia were treated with 100 ng/mL lipopolysaccharide for 24 hours to establish in vitro cell models, and then transfected with pcDNA-HOXA11-AS, miR-124-3p mimic, or sh-MDK. The results revealed that HOXA11-AS inhibited miR-124-3p expression and boosted MDK expression and TLR4-nuclear factor-κB pathway activation. Furthermore, lipopolysaccharide enhanced potent microglia-induced inflammatory responses in astrocytes. Forced overexpression of miR-124-3p or downregulating MDK repressed microglial activation and the inflammatory response of astrocytes. However, the miR-124-3p-mediated anti-inflammatory effects were reversed by HOXA11-AS. These findings suggest that HOXA11-AS can aggravate neuroinflammation after traumatic brain injury by modulating the miR-124-3p-MDK axis. This study was approved by the Animal Protection and Use Committee of Southwest Medical University (approval No. SMU-2019-042) on February 4, 2019.

Long Non-Coding RNAs in Diagnosis, Treatment, Prognosis, and Progression of Glioma: A State-of-the-Art Review

Glioma is the most common malignant central nervous system tumor with significant mortality and morbidity. Despite considerable advances, the exact molecular pathways involved in tumor progression are not fully elucidated, and patients commonly face a poor prognosis. Long non-coding RNAs (lncRNAs) have recently drawn extra attention for their potential roles in different types of cancer as well as non-malignant diseases. More than 200 lncRNAs have been reported to be associated with glioma. We aimed to assess the roles of the most investigated lncRNAs in different stages of tumor progression and the mediating molecular pathways in addition to their clinical applications. lncRNAs are involved in different stages of tumor formation, invasion, and progression, including regulating the cell cycle, apoptosis, autophagy, epithelial-to-mesenchymal transition, tumor stemness, angiogenesis, the integrity of the blood-tumor-brain barrier, tumor metabolism, and immunological responses. The well-known oncogenic lncRNAs, which are upregulated in glioma, are H19, HOTAIR, PVT1, UCA1, XIST, CRNDE, FOXD2-AS1, ANRIL, HOXA11-AS, TP73-AS1, and DANCR. On the other hand, MEG3, GAS5, CCASC2, and TUSC7 are tumor suppressor lncRNAs, which are downregulated. While most studies reported oncogenic effects for MALAT1, TUG1, and NEAT1, there are some controversies regarding these lncRNAs. Expression levels of lncRNAs can be associated with tumor grade, survival, treatment response (chemotherapy drugs or radiotherapy), and overall prognosis. Moreover, circulatory levels of lncRNAs, such as MALAT1, H19, HOTAIR, NEAT1, TUG1, GAS5, LINK-A, and TUSC7, can provide non-invasive diagnostic and prognostic tools. Modulation of expression of lncRNAs using antisense oligonucleotides can lead to novel therapeutics. Notably, a profound understanding of the underlying molecular pathways involved in the function of lncRNAs is required to develop novel therapeutic targets. More investigations with large sample sizes and increased focus on in-vivo models are required to expand our understanding of the potential roles and application of lncRNAs in glioma.

HOXA11 plays critical roles in disease progression and response to cytarabine in AML

Lysine methyltransferase 2A (KMT2A, also known as MLL) translocations (MLL-t) are frequently associated with mutations in RAS pathway genes in acute myeloid leukemia (AML). Previous findings with a mouse model showed that cooperation of MLL/AF10 with tyrosine-protein phosphatase non-receptor type 11 (PTPN11)^G503A accelerated leukemia development, but increased cytarabine (Ara-C) sensitivity of leukemia cells. To identify the genes responsible for reduced survival and Ara-C resistance, transcriptomic profiling between six pairs of mouse MLL/AF10(OM-LZ) leukemia cells harboring activating and wild-type KRAS or PTPN11 was compared. A total of 23 differentially expressed genes (DEGs) with >1.5-fold-change between the paired cell lines were identified. The Gene Ontology (GO) terms overrepresented in these 23 DEGs included ‘immune system process’, ‘actin filament binding’, ‘cellular response to interferon-alpha’ and ‘sequence-specific DNA’. Among the four genes (Hoxa11, PR domain zinc finger protein 5, Iroquois-class homeodomain protein IRX-5 and homeobox protein PKNOX2) mapped to the GO term ‘sequence-specific DNA’, HOXA11 upregulation was associated with AML harboring MLL-t and RAS signaling mutations based on a meta-analysis using data deposited in Oncomine™ and analysis of the clinical samples in the present study. Microarray data revealed that only Hoxa11 was upregulated in those cells harboring activating PTPN11. Functional studies of Hoxa11 knockdown or overexpression in MLL/AF10(OM-LZ) cells revealed that Hoxa11 expression levels were associated with survival in vivo and Ara-C sensitivity/apoptosis in vitro. In addition, Hoxa11 regulated the expression of the apoptosis-related genes, NF-κB inhibitor α, transcription factor p65 and transformation-related protein p53. Furthermore, the results of a meta-analysis using Heuser's AML dataset supported the finding that chemotherapy responders have higher expression levels of HOXA11. These results indicated that the expression of HOXA11 increased cell apoptosis and predicted an improved response to Ara-C in AML.

Downregulation of HOXA11 enhances endometrial cancer malignancy and cisplatin resistance via activating PTEN/AKT signaling pathway

PURPOSE

Endometrial cancer is the most common malignant tumor of female genital system worldwide. Homeobox A11 (HOXA11) is an evolutionarily conserved Homeobox gene closely implicated in carcinogenesis. However, the mechanisms of HOXA11 in the progression and cisplatin resistance of endometrial cancer remain unclear.

METHODS

The expression of HOXA11 was analyzed based on 548 endometrial cancer and 35 control tissues from The Cancer Genome Atlas (TCGA) database. Transwell assay was performed to investigate the effect of HOXA11 on endometrial cell migration and invasion. TUNEL staining was carried out to assay the role of HOXA11 in endometrial cell apoptosis. Western blot was employed to detect the protein levels of B cell lymphoma-2 (Bcl-2), Bcl-2 associated X (Bax), cleaved caspase-3, matrix metalloproteinase-2/9 (MMP/9), phosphatase and tensin homolog (PTEN), protein kinase B (AKT) and p-AKT.

RESULTS

TCGA data showed that HOXA11 expression was significantly down-regulated in endometrial cancer tissue samples. The overexpression of HOXA11 promoted the apoptosis, but inhibited the proliferation, migration and invasion of endometrial cancer cells. HOXA11 knockdown with small interfering RNA (siRNA) considerably repressed cell apoptosis, while promoted cell proliferation, migration, and invasion through PTEN/AKT signaling pathway. Interestingly, HOXA11 was lowly expressed in Ishikawa cells treated with cisplatin. In addition, HOXA11 knockdown increased the resistance of endometrial cancer to cisplatin through activating PTEN/AKT signaling pathway.

CONCLUSION

Low HOXA11 expression may promote the proliferation, migration, invasion of endometrial cancer cells, and increase their resistance to cisplatin through activating PTEN/AKT pathway.

HOX gene cluster (de)regulation in brain: from neurodevelopment to malignant glial tumours

HOX genes encode a family of evolutionarily conserved homeodomain transcription factors that are crucial both during development and adult life. In humans, 39 HOX genes are arranged in four clusters (HOXA, B, C, and D) in chromosomes 7, 17, 12, and 2, respectively. During embryonic development, particular epigenetic states accompany their expression along the anterior-posterior body axis. This tightly regulated temporal-spatial expression pattern reflects their relative chromosomal localization, and is critical for normal embryonic brain development when HOX genes are mainly expressed in the hindbrain and mostly absent in the forebrain region. Epigenetic marks, mostly polycomb-associated, are dynamically regulated at HOX loci and regulatory regions to ensure the finely tuned HOX activation and repression, highlighting a crucial epigenetic plasticity necessary for homeostatic development. HOX genes are essentially absent in healthy adult brain, whereas they are detected in malignant brain tumours, namely gliomas, where HOX genes display critical roles by regulating several hallmarks of cancer. Here, we review the major mechanisms involved in HOX genes (de)regulation in the brain, from embryonic to adult stages, in physiological and oncologic conditions. We focus particularly on the emerging causes of HOX gene deregulation in glioma, as well as on their functional and clinical implications.

N6-Methylandenosine-Related lncRNAs Are Potential Biomarkers for Predicting the Overall Survival of Lower-Grade Glioma Patients

The prognostic value of N6-methylandenosine-related long non-coding RNAs (m6A-related lncRNAs) was investigated in 646 lower-grade glioma (LGG) samples from The Cancer Genome Atlas (TCGA) and the Chinese Glioma Genome Atlas (CGGA) datasets. We implemented Pearson correlation analysis to explore the m6A-related lncRNAs, and then univariate Cox regression analysis was performed to screen their prognostic roles in LGG patients. Twenty-four prognostic m6A-related lncRNAs were identified as prognostic lncRNAs and they were inputted in a least absolute shrinkage and selection operator (LASSO) Cox regression to establish a m6A-related lncRNA prognostic signature (m6A-LPS, including 9 m6A-related prognostic lncRNAs) in the TCGA dataset. Corresponding risk scores of patients were calculated and divided LGG patients into low- and high-risk subgroups by the median value of risk scores in each dataset. The m6A-LPS was validated in the CGGA dataset and it showed a robust prognostic ability in the stratification analysis. Principal component analysis showed that the low- and high-risk subgroups had distinct m6A status. Enrichment analysis indicated that malignancy-associated biological processes, pathways and hallmarks were more common in the high-risk subgroup. Moreover, we constructed a nomogram (based on m6A-LPS, age and World Health Organization grade) that had a strong ability to forecast the overall survival (OS) of the LGG patients in both datasets. We also establish a competing endogenous RNA (ceRNA) network based on seven of the twenty-four m6A-related lncRNAs. Besides, we also detected five m6A-related lncRNA expression levels in 22 clinical samples using quantitative real-time polymerase chain reaction assay.

Context-dependent HOX transcription factor function in health and disease

During animal development, HOX transcription factors determine the fate of developing tissues to generate diverse organs and appendages. The power of these proteins is striking: mis-expressing a HOX protein causes homeotic transformation of one body part into another. During development, HOX proteins interpret their cellular context through protein interactions, alternative splicing, and post-translational modifications to regulate cell proliferation, cell death, cell migration, cell differentiation, and angiogenesis. Although mutation and/or mis-expression of HOX proteins during development can be lethal, changes in HOX proteins that do not pattern vital organs can result in survivable malformations. In adults, mutation and/or mis-expression of HOX proteins disrupts their gene regulatory networks, deregulating cell behaviors and leading to arthritis and cancer. On the molecular level, HOX proteins are composed of DNA binding homeodomain, and large regions of unstructured, or intrinsically disordered, protein sequence. The primary roles of HOX proteins in arthritis and cancer suggest that mutations associated with these diseases in both the structured and disordered regions of HOX proteins can have substantial functional effects. These insights lead to new questions critical for understanding and manipulating HOX function in physiological and pathological conditions.

Methylation in HOX Clusters and Its Applications in Cancer Therapy

HOX genes are commonly known for their role in embryonic development, defining the positional identity of most structures along the anterior–posterior axis. In postembryonic life, HOX gene aberrant expression can affect several processes involved in tumorigenesis such as proliferation, apoptosis, migration and invasion. Epigenetic modifications are implicated in gene expression deregulation, and it is accepted that methylation events affecting HOX gene expression play crucial roles in tumorigenesis. In fact, specific methylation profiles in the HOX gene sequence or in HOX-associated histones are recognized as potential biomarkers in several cancers, helping in the prediction of disease outcomes and adding information for decisions regarding the patient’s treatment. The methylation of some HOX genes can be associated with chemotherapy resistance, and its identification may suggest the use of other treatment options. The use of epigenetic drugs affecting generalized or specific DNA methylation profiles, an approach that now deserves much attention, seems likely to be a promising weapon in cancer therapy in the near future. In this review, we summarize these topics, focusing particularly on how the regulation of epigenetic processes may be used in cancer therapy.

Expression Profile and Prognostic Values of HOXA Family Members in Laryngeal Squamous Cell Cancer

The homeobox A cluster (HOXA) gene family, comprising 11 members, is involved in a wide spectrum of biological functions in human cancers. However, there is little research on the expression profile and prognostic values of HOXA genes in laryngeal squamous cell cancer (LSCC). Based on updated public resources and integrative bioinformatics analysis, we assessed the expression profile and prognostic values of the HOXA family members. Expression and methylation data on HOXA family members were obtained from The Cancer Genome Atlas (TCGA). The prognostic values of HOXA members and clinical features were identified. A gene set enrichment analysis (GSEA) was conducted to explore the mechanism underlying the involvement of HOXA members in LSCC. The associations between tumor immune infiltrating cells (TIICs) and the HOXA family members were evaluated using the Tumor Immune Estimation Resource (TIMER) database. HOXA2 and HOXA4 were downregulated and HOXA7 and HOXA9–13 were upregulated in LSCC. Upregulation of HOXA10, HOXA11, and HOXA13, along with two clinical characteristics (M stage and gender), were associated with a poor LSCC prognosis based on the results of univariate and multivariate Cox proportional hazards regression analyses. Although there were no significant correlations between TIICs and HOXA members, the GSEA results indicated that HOXA members participate in multiple biological processes underlying tumorigenesis. This study comprehensively analyzed the HOXA members, providing insights for further investigation of the HOXA family members as potential targets in LSCC.

Long Non-Coding RNA HOXA11-AS Promotes Non-Small Cell Lung Cancer Tumorigenesis Through microRNA-148a-3p/DNMT1 Regulatory Axis

OBJECTIVE

Our present study aimed to further investigate the molecular basis of long non-coding RNA homeobox A11 antisense (HOXA11-AS) in the tumorigenesis of non-small cell lung cancer (NSCLC).

METHODS

HOXA11-AS, microRNA-148a-3p (miR-148a-3p), and DNA methyltransferase 1 (DNMT1) mRNA levels were measured by RT-qPCR assay. DNMT1 protein level was determined by Western blot assay. Cell proliferative capacity and apoptotic rate were determined by CCK-8 assay and flow cytometry analysis, respectively. The relationships of HOXA11-AS, miR-148a-3p, and DNMT1 were tested through bioinformatics analysis, luciferase assay, and RNA pull down assay. Mouse xenograft models of NSCLC were established to examine the biological function of HOXA11-AS in vivo.

RESULTS

HOXA11-AS expression was notably upregulated and miR-148a-3p expression was conspicuously downregulated in NSCLC tissues and cells. HOXA11-AS knockdown curbed NSCLC cell proliferation and promoted cell apoptosis through directly increasing miR-148a-3p expression. Moreover, miR-148a-3p overexpression suppressed NSCLC cell proliferation and induced cell apoptosis. HOXA11-AS functioned as a competing endogenous RNA (ceRNA) of miR-148a-3p to increase DNMT1 expression in NSCLC cells. And, DNMT1 upregulation weakened the influence of HOXA11-AS1 loss on NSCLC cell proliferation and apoptosis. Additionally, HOXA11-AS knockdown suppressed NSCLC xenograft growth by upregulating miR-148a-3p and downregulating DNMT1 in vivo.

CONCLUSION

HOXA11-AS facilitated NSCLC tumorigenesis through miR-148a-3p/DNMT1 axis in vitro and in vivo, deepening our understanding of the molecular basis of HOXA11-AS in the development of NSCLC.

A self-enforcing HOXA11/Stat3 feedback loop promotes stemness properties and peritoneal metastasis in gastric cancer cells

Rationale: Peritoneal metastasis is one of the most common and life-threatening metastases in gastric cancer patients. The disseminated gastric cancer cells forming peritoneal metastasis exhibit a variety of characteristics that contrast with those of adjacent epithelial cell of gastric mucosa and even primary gastric cancer cells. We hypothesized that the gene expression profiles of peritoneal foci could reveal the identities of genes that might function as metastatic activator.

Methods: In this study, we show, using in vitro, in vivo, in silico and gastric cancer tissues studies in humans and mice, that Homoebox A11 (HOXA11) potently promote peritoneal metastasis of gastric cancer cells.

Results: Its mechanism of action involves alternation of cancer stemness and subsequently enhancement of the adhesion, migration and invasion and anti-apoptosis. This is achieved, mainly, through formation of a positive feedback loop between HOXA11 and Stat3, which is involved in the stimulation of Stat3 signaling pathway.

Conclusions: These observations uncover a novel peritoneal metastatic activator and demonstrate the association between HOXA11, Stat3 and cancer stemness of gastric cancer cells, thereby revealing a previously undescribed mechanism of peritoneal metastasis.

MicroRNA-128/homeobox B8 axis regulates ovarian cancer cell progression

MicroRNA-128 (miR-128) has been found to be dysregulated and might function as a tumour suppressor in various cancers, including ovarian cancer. However, the underlying mechanism of miR-128 in ovarian cancer has not been fully understood. The miR-128 and homeobox B8 (HOXB8) levels in clinical samples and cultured cell lines were measured using qRT-PCR and/or Western blot analysis. Cell proliferation was assessed using Cell Counting Kit-8 assay. Cell apoptosis was determined using flow cytometry. The association between miR-128 and HOXB8 was confirmed using dual-luciferase reporter assay. Results showed that decreased miR-128 expression and increased HOXB8 expression were observed in ovarian cancer tissues and cell lines. Transfection with miR-128 mimics suppressed the cell proliferation and enhanced paclitaxel sensitivity in ovarian cancer cell lines. miR-128 directly targeted HOXB8 in ovarian cancer cell lines. Knockdown of HOXB8 abolished the effects of miR-128 inhibitor on ovarian cancer cell proliferation and paclitaxel sensitivity. Summarily, miR-128 displayed a tumour suppressor role in ovarian cancer via targeting HOXB8. It is supposed that miR-128 might be effective for targeting therapy for ovarian cancer.

DNA methylation changes following DNA damage in prostate cancer cells

Many cancer therapies operate by inducing double-strand breaks (DSBs) in cancer cells, however treatment-resistant cells rapidly initiate mechanisms to repair damage enabling survival. While the DNA repair mechanisms responsible for cancer cell survival following DNA damaging treatments are becoming better understood, less is known about the role of the epigenome in this process. Using prostate cancer cell lines with differing sensitivities to radiation treatment, we analysed the DNA methylation profiles prior to and following a single dose of radiotherapy (RT) using the Illumina Infinium HumanMethylation450 BeadChip platform. DSB formation and repair, in the absence and presence of the DNA hypomethylating agent, 5-azacytidine (5-AzaC), were also investigated using γH2A.X immunofluorescence staining. Here we demonstrate that DNA methylation is generally stable following a single dose of RT; however, a small number of CpG sites are stably altered up to 14 d following exposure. While the radioresistant and radiosensitive cells displayed distinct basal DNA methylation profiles, their susceptibility to DNA damage appeared similar demonstrating that basal DNA methylation has a limited influence on DSB induction at the regions examined. Recovery from DSB induction was also similar between these cells. Treatment with 5-AzaC did not sensitize resistant cells to DNA damage, but rather delayed recruitment of phosphorylated BRCA1 (S1423) and repair of DSBs. These results highlight that stable epigenetic changes are possible following a single dose of RT and may have significant clinical implications for cancer treatment involving recurrent or fractionated dosing regimens.

Long noncoding RNA HOXA11-AS promotes gastric cancer cell proliferation and invasion via SRSF1 and functions as a biomarker in gastric cancer

BACKGROUND

Gastric cancer (GC) is the fourth most frequent malignancy all over the world. The diagnosis of GC is challenging and the prognosis of GC is very unfavorable. Accumulating evidence reveals that serum long noncoding RNAs (lncRNAs) can function as biomarkers in various types of cancers, including GC.

AIM

To explore the level and molecular mechanism of the lncRNA HOXA11-AS in GC and the diagnostic and prognostic significance of serum HOXA11-AS in GC.

METHODS

HOXA11-AS levels in GC tissue, cell lines, and serum samples were measured. The correlation between HOXA11-AS expression and clinicopathological characteristics was analyzed. The role of HOXA11-AS in the diagnosis and prognosis of GC was evaluated. Cell function assays were performed for exploration of the roles of HOXA11-AS in GC cells. Moreover, Western blot was performed to explore the target regulated by HOXA11-AS in GC cells.

RESULTS

Up-regulation of HOXA11-AS was found in GC tissues, cell lines, and serum samples. In GC patients, decreased serum HOXA11-AS levels were negatively related with tumor size, TNM stage, and lymph node metastasis. The area under the receiver operating characteristic curve of serum HOXA11-AS in the diagnosis of GC was 0.924 (95%CI: 0.881-0.967; sensitivity, 0.787; specificity 0.978). Results of the Kaplan-Meier survival curves suggested the GC patients with a lower HOXA11-AS level having a better overall survival rate. HOXA11-AS promoted GC cell proliferation and invasion. SRSF1 may be the target regulated by HOXA11-AS in GC cells.

CONCLUSION

HOXA11-AS promotes GC cell proliferation and invasion via SRSF1 and may function as a promising marker in GC.

An inverse interaction between HOXA11 and HOXA11-AS is associated with cisplatin resistance in lung adenocarcinoma

HOXA11, which is a member of the homeobox (HOX) gene family, and its natural antisense transcript (NAT) HOXA11-AS have been reported to be closely related to the development of lung cancer. We aimed to investigate their specific roles in cisplatin (DDP) resistance in lung adenocarcinoma (LUAD). First, we found that HOXA11 is hypermethylated and significantly downregulated in a DDP-resistant A549 cell line (A549/DDP) and LUAD tissues, while the HOXA11-AS expression level is elevated. Although HOXA11 and HOXA11-AS mRNA overlap in the 5'-untranslated region (5' UTR) and share two CpG islands, DNA methylation only regulates the expression of HOXA11. Then, we found that HOXA11 and HOXA11-AS have an inverse interaction by transfecting their siRNAs and overexpression vectors into A549 and A549/DDP cells. A dual-luciferase reporter assay further confirmed that the overlapping 5'UTR is essential for the bidirectional regulation between HOXA11 and HOXA11-AS. Functional analysis showed that knockdown of HOXA11 expression in A549 cells induced DDP resistance and activated Akt/β-catenin signaling, while overexpression of HOXA11 in A549/DDP cells increased DDP sensitivity and inhibited Akt/β-catenin signaling. Moreover, HOXA11-AS knockdown in A549 cells increased DDP sensitivity and inhibited Akt/β-catenin signaling, while the overexpression of HOXA11-AS in A549/DDP cells induced DDP resistance and activated Akt/β-catenin signaling. In conclusion, our study demonstrates that the inverse interaction between HOXA11 and HOXA11-AS promotes DDP resistance in LUAD.

Targeting HOX/PBX dimer formation as a potential therapeutic option in esophageal squamous cell carcinoma

Homeobox genes are known to be classic examples of the intimate relationship between embryogenesis and tumorigenesis, which are a family of transcriptional factors involved in determining cell identity during early development, and also dysregulated in many malignancies. Previously, HOXB7, HOXC6 and HOXC8 were found abnormally upregulated in esophageal squamous cell carcinoma (ESCC) tissues compared with normal mucosa and seen as poor prognostic predictors for ESCC patients, and were shown to promote cell proliferation and anti‐apoptosis in ESCC cells. These three HOX members have a high level of functional redundancy, making it difficult to target a single HOX gene. The aim of the present study was to explore whether ESCC cells are sensitive to HXR9 disrupting the interaction between multiple HOX proteins and their cofactor PBX, which is required for HOX functions. ESCC cell lines (KYSE70, KYSE150, KYSE450) were treated with HXR9 or CXR9, and coimmunoprecipitation and immunofluorescent colocalization were carried out to observe HOX/PBX dimer formation. To further investigate whether HXR9 disrupts the HOX pro‐oncogenic function, CCK‐8 assay and colony formation assay were carried out. Apoptosis was assessed by flow cytometry, and tumor growth in vivo was investigated in a xenograft model. RNA‐seq was used to study the transcriptome of HXR9‐treated cells. Results showed that HXR9 blocked HOX/PBX interaction, leading to subsequent transcription alteration of their potential target genes, which are involved in JAK‐signal transducer and activator of transcription (STAT) activation and apoptosis inducement. Meanwhile, HXR9 showed an antitumor phenotype, such as inhibiting cell proliferation, inducing cell apoptosis and significantly retarding tumor growth. Therefore, it is suggested that targeting HOX/PBX may be a novel effective treatment for ESCC.

Upregulation of HOXA11 during the progression of lung adenocarcinoma detected via multiple approaches

The altered expression of homeobox (HOX)A11 has been observed in various malignant tumor types, but it has remained to be determined in human lung adenocarcinoma (LUAD). In the present study, the expression of HOXA11 in LUAD and the potential associated mechanisms were assessed. Data from The Cancer Genome Atlas and Oncomine microarrays were gathered and in‑house polymerase chain reaction data were produced to investigate the altered expression of HOXA11 in LUAD and its association with various clinicopathological characteristics. Genes co‑expressed with HOXA11 were also identified by searching the cBioPortal and Multi Experiment Matrix databases, and performing a bioinformatics analysis, through which the potential molecular mechanisms of HOXA11 in LUAD were explored. The data analyses indicated that HOXA11 was overexpressed in the LUAD samples, and together with its co‑expressed genes, it was indicated to participate in various key signaling pathways, including the focal adhesion, extracellular matrix‑receptor interaction, axon guidance and small cell lung cancer signaling pathways. Furthermore, collagen type III α 1 chain (COL3A1), ephrin B2 (EFNB2), integrin subunit α 8 (ITGA8) and syndecan 2 (SDC2) were confirmed to be differentially expressed in LUAD vs. normal controls at the mRNA and protein level. Of note, LUAD patients with low expression of HOXA11 and ITGB1 had better overall survival rates. The present study indicated that HOXA11 may function as an oncogene in LUAD, and HOXA11 protein probably combines with ITGB1, COL3A1, EFNB2, ITGA8 and SDC2 to have a role in the focal adhesion pathway.

HOXA11 antisense long noncoding RNA (HOXA11-AS): A promising lncRNA in human cancers

The cancers are the leading cause of disease‐related deaths worldwide with a high risk of morbidity and mortality. Long noncoding RNAs (lncRNAs) play a critical role in a wide range of biological processes, including tumorigenesis. HOXA11‐AS (NCRNA00076), the antisense strands of HOXA11 gene, was initially revealed in a mouse embryonic cDNA library in 2009 and it was a fairly novel lncRNA. This review summarized the advanced research progression concerning the expression and role of HOXA11‐AS in different human malignancies. The expression of HOXA11‐AS is aberrantly altered in many cancers, either as a tumor suppressor or as a tumor accelerator. The different underlying mechanism of HOXA11‐AS in different cancers (including, nonsmall cell lung cancers, osteosarcoma, uveal melanoma, glioma, hepatocellular carcinoma, gastric cancer, breast cancer, cervical cancer, ovarian cancer, colorectal cancer, ovarian cancer, and glioblastoma) was also detailed. These findings lead us to conclude that HOXA11‐AS participate in the complex network of cancers and plays an important role in the tumorigenesis and progression. Functional HOXA11‐AS could be a promising biomarker for early detection as well as prognosis evaluation in cancer patients. Future HOXA11‐AS‐targeted intervention may become a valuable novel therapeutic tool, improving the clinical management of cancers.

Evaluation of the HOXA11 level in patients with lung squamous cancer and insights into potential molecular pathways via bioinformatics analysis

Background

This study was carried out to discover the underlying role that HOXA11 plays in lung squamous cancer (LUSC) and uncover the potential corresponding molecular mechanisms and functions of HOXA11-related genes.

Methods

Twenty-three clinical paired LUSC and non-LUSC samples were utilized to examine the level of HOXA11 using quantitative real-time polymerase chain reaction (qRT-PCR). The clinical significance of HOXA11 was systematically analyzed based on 475 LUSC and 18 non-cancerous adjacent tissues from The Cancer Genome Atlas (TCGA) database. A total of 102 LUSC tissues and 121 non-cancerous tissues were available from Oncomine to explore the expressing profiles of HOXA11 in LUSC. A meta-analysis was carried out to further assess the differential expression of HOXA11 in LUSC, including in-house qRT-PCR data, expressing data extracted from TCGA and Oncomine databases. Moreover, the enrichment analysis and potential pathway annotations of HOXA11 in LUSC were accomplished via Gene Oncology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). The expression of hub genes and according correlations with HOXA11 were assessed to further explore the biological role of HOXA11 in LUSC.

Results

HOXA11 expression in LUSC had a tendency to be upregulated in comparison to adjacent non-cancerous tissues by qRT-PCR. TCGA data displayed that HOXA11 was remarkably over-expressed in LUSC compared with that in non-LUSC samples, and the area under curves (AUC) was 0.955 (P < 0.001). A total of 1523 co-expressed genes were sifted for further analysis. The most significant term enriched in the KEGG pathway was focal adhesion. Among the six hub genes of HOXA11, including PARVA, ILK, COL4A1, COL4A2, ITGB1, and ITGA5, five (with the exception of COL4A1) were significantly decreased compared with the normal lung tissues. Moreover, the expression of ILK was negatively related to HOXA11 (r = − 0.141, P = 0.002).

Conclusion

High HOXA11 expression may lead to carcinogenesis and the development of LUSC. Furthermore, co-expressed genes might affect the prognosis of LUSC.

Identification of radiation responsive genes and transcriptome profiling via complete RNA sequencing in a stable radioresistant U87 glioblastoma model

The absence of major progress in the treatment of glioblastoma (GBM) is partly attributable to our poor understanding of both GBM tumor biology and the acquirement of treatment resistance in recurrent GBMs. Recurrent GBMs are characterized by their resistance to radiation. In this study, we used an established stable U87 radioresistant GBM model and total RNA sequencing to shed light on global mRNA expression changes following irradiation. We identified many genes, the expressions of which were altered in our radioresistant GBM model, that have never before been reported to be associated with the development of radioresistant GBM and should be concertedly further investigated to understand their roles in radioresistance. These genes were enriched in various biological processes such as inflammatory response, cell migration, positive regulation of epithelial to mesenchymal transition, angiogenesis, apoptosis, positive regulation of T-cell migration, positive regulation of macrophage chemotaxis, T-cell antigen processing and presentation, and microglial cell activation involved in immune response genes. These findings furnish crucial information for elucidating the molecular mechanisms associated with radioresistance in GBM. Therapeutically, with the global alterations of multiple biological pathways observed in irradiated GBM cells, an effective GBM therapy may require a cocktail carrying multiple agents targeting multiple implicated pathways in order to have a chance at making a substantial impact on improving the overall GBM survival.

Identification of skin-related lncRNAs as potential biomarkers that involved in Wnt pathways in keloids

The long non-coding RNAs (lncRNAs) regulating encoding transcripts/genes involved in Wnt signalling pathway in keloids is largely unclear. We used a pathway-focused lncRNA microarray to detect the differentiated expression profiles of both lncRNAs and genes involved in Wnt pathway, thus a total of 116 Wnt-targeted genes and 69 Wnt-related lncRNAs aberrantly expressed in keloids were initially identified. A stepwise bioinformatics was further performed to find skin-related lncRNA/gene pairs in Wnt pathway in keloids. Firstly, an lncRNA/gene co-expression network with clustered functional modules was constructed; simultaneously, 114 Wnt-genes regarding to dermis were online enriched using Phenotype Enrichment. Secondly, 17 skin-related keloid-aberrant Wnt-genes were acquired by overlapping the 114 skin-related Wnt-genes with the 116 keloid-aberrant Wnt-genes. Thirdly, after co-expression coefficient of each lncRNA/gene profile being ranked respectively, 11 top co-expressed lncRNAs characterized with the highest co-expression coefficients to the 17 genes were identified. Fourthly, seven of the 11 top co-expressed lncRNAs exhibiting array-detected aberrant expression in keloids, together with their 12 most interactive Wnt-genes, were selected to undergo in-pair intracellularly quantitative PCR validation in keloids. As a result, four lncRNAs including CACNA1G-AS1, HOXA11-AS, LINC00312 and RP11-91I11.1 with their six paired Wnt-genes undergoing both array-and-qPCR as well as lncRNA-and-gene double validation were finally identified as skin-related lncRNA/gene pairs that involved in Wnt signalling pathway in keloids. In conclusion, in-depth exploration on these easily-accessible lncRNAs in keloids might aid to find the novel target on how to maintain highly recurrent tumours benign via Wnt-involved network regulation.

Functional analysis of HOXA10 and HOXB4 in human medulloblastoma cell lines

Medulloblastoma (MB) is a malignant childhood brain tumor which at molecular level is classified into at least four major subtypes: WNT, SHH, group C and group D differing in response to treatment. Previous studies have associated changes in expression levels and activation of certain HOX genes with MB development. In the present study, we investigate the role of HOX genes in two attributes acquired by tumor cells: migration and proliferation potential, as well as, in vivo tumorigenic potential. We analyzed UW402, UW473, DAOY and ONS-76 human pediatric MB cell lines and cerebellum primary cultures. Two-color microarray-based gene expression analysis was used to identify differentially expressed HOX genes. Among the various HOX genes significantly overexpressed in DAOY and ONS-76 cell lines compared to UW402 and UW473 cell lines, HOXA10 and HOXB4 were selected for further analysis. The expression levels of these HOX genes were validated by real-time PCR. A mouse model was used to study the effect of the HOXA10 and HOXB4 genes on the in vivo tumorigenic potential and the in vitro proliferative and migration potential of MB cell lines. Our results show that the inhibition of HOXA10 in DAOY cell line led to increased in vitro cell migration while in vitro cell proliferation or in vivo tumorigenic potential were unaffected. We also observed that induced expression of HOXB4 in the UW473 cell line significantly reduced in vitro cell proliferation and migration capability of UW473 cells with no effect on the in vivo tumorigenicity. This suggests that HOXA10 plays a role in migration events and the HOXB4 gene is involved in proliferation and migration processes of medulloblastoma cells, however, it appears that these genes are not essential for the tumorigenic process of these cells.

Long non‑coding RNA HOXA11‑AS promotes cell proliferation and metastasis in human breast cancer

Breast cancer is one of the most frequently occurring malignancies in female cancers worldwide, however, its detailed mechanism of tumorigenesis remains to be elucidated. Long non-coding RNAs (LncRNAs) have previously been demonstrated to be important in multiple cancers, including breast cancer. The present study aimed to elucidate the molecular mechanism of the effects of the novel Lnc RNA HOXA11‑AS, on cell proliferation and metastasis in breast cancer. The data revealed that the relative transcript level of HOXA11‑AS was upregulated in vivo and in vitro in models of breast cancer. Knockdown of HOXA11‑AS in MDA‑MB‑231 and MDA‑MB‑436 breast cancer cell lines inhibited the formation of cell colonies and arrested the cell cycle at the G0/G1 phase. Depletion of HOXA11‑AS using two specific short interfering (si)RNAs against HOXA11‑AS (siHOXA11‑AS‑1 and siHOXA11‑AS‑2) additionally suppressed the cell proliferative rate. Furthermore, transwell assays and wound‑healing analysis revealed that siRNA transfection inhibited cell migration and invasion by ~50% in the two cell lines. The results of the present study demonstrated the oncogenic role of HOXA11‑AS in breast cancer, providing novel clues for the future clinical diagnosis and treatment of early stage breast cancer patients.

Clinical Significance and Effect of lncRNA HOXA11-AS in NSCLC: A Study Based on Bioinformatics, In Vitro and in Vivo Verification

HOXA11 antisense RNA (HOXA11-AS) has been shown to be involved in tumorigenesis and development of different cancers. However, the role of HOXA11-AS in non-small cell lung cancer (NSCLC) remains unclear. In this study, we firstly explored and confirmed the expression of HOXA11-AS in NSCLC tissues and cells. Cytometry, CCK-8, cell scratch, migration, Matrigel invasion and flow cytometry assays were performed to determine the biological impact of HOXA11-AS in vitro. Furthermore, a chick embryo chorioallantoic membrane (CAM) model of NSCLC was constructed to explore the effect of HOXA11-AS on tumorigenicity and angiogenesis in vivo. Additionally, bioinformatics analyses were performed to investigate the prospective pathways of HOXA11-AS co-expressed genes. As results, HOXA11-AS was markedly highly expressed in NSCLC tissues and cells. Furthermore, the proliferation, migration, invasion, tumorigenic and angiogenic ability of NSCLC cells were all inhibited and apoptosis was induced after HOXA11-AS knock-down. HOXA11-AS RNAi also led to cell cycle arrest on G0/G1 or G2/M phase. In addition, the non-small cell lung cancer pathway might be involved in regulating the co-expressed genes of HOXA11-AS in NSCLC. These results indicate that HOXA11-AS plays pivotal roles in NSCLC and it can become a novel therapeutic direction for treating NSCLC.

Microarray-based analysis and clinical validation identify ubiquitin-conjugating enzyme E2E1 (UBE2E1) as a prognostic factor in acute myeloid leukemia

Background

Previous research suggested that single gene expression might be correlated with acute myeloid leukemia (AML) survival. Therefore, we conducted a systematical analysis for AML prognostic gene expressions.

Methods

We performed a microarray-based analysis for correlations between gene expression and adult AML overall survival (OS) using datasets GSE12417 and GSE8970. Positive findings were validated in an independent cohort of 50 newly diagnosed, non-acute promyelocytic leukemia (APL) AML patients by quantitative RT-PCR and survival analysis.

Results

Microarray-based analysis suggested that expression of eight genes was each associated with 1-year and 3-year AML OS in both GSE12417 and GSE8970 datasets (p < 0.05). Next, we validated our findings in an independent cohort of AML samples collected in our hospital. We found that ubiquitin-conjugating enzyme E2E1 (UBE2E1) expression was adversely correlated with AML survival (p = 0.04). Multivariable analysis showed that UBE2E1^high patients had a significant shorter OS and shorter progression-free survival after adjusting other known prognostic factors (p = 0.03). At last, we found that UBE2E1 expression was negatively correlated with patients’ response to induction chemotherapy (p < 0.05).

Conclusions

In summary, we demonstrated that UBE2E1 expression was a novel prognostic factor in adult, non-APL AML patients.

Electronic supplementary material

The online version of this article (doi:10.1186/s13045-016-0356-0) contains supplementary material, which is available to authorized users.