HOXA10‑AS: A novel oncogenic long non‑coding RNA in glioma

The tip of the iceberg-The roles of long noncoding RNAs in acute myeloid leukemia

Long noncoding RNAs (lncRNAs) are traditionally defined as RNA transcripts longer than 200 nucleotides that have no protein coding potential. LncRNAs have been identified to be dysregulated in various types of cancer, including the deadly hematopoietic cancer-acute myeloid leukemia (AML). Currently, survival rates for AML have reached a plateau necessitating new therapeutic targets and biomarkers to improve treatment options and survival from the disease. Therefore, the identification of lncRNAs as novel biomarkers and therapeutic targets for AML has major benefits. In this review, we assess the key studies which have recently identified lncRNAs as important molecules in AML and summarize the current knowledge of lncRNAs in AML. We delve into examples of the specific roles of lncRNA action in AML such as driving proliferation, differentiation block and therapy resistance as well as their function as tumor suppressors and utility as biomarkers. This article is categorized under: RNA in Disease and Development > RNA in Disease.

Comparative computational analysis to distinguish mesenchymal stem cells from fibroblasts

Introduction

Mesenchymal stem cells (MSCs) are considered to be the most promising stem cell type for cell-based therapies in regenerative medicine. Based on their potential to home to diseased body sites following a therapeutically application, these cells could (i) differentiate then into organ-specific cell types to locally restore injured cells or, most prominently, (ii) foster tissue regeneration including immune modulations more indirectly by secretion of protective growth factors and cytokines. As tissue-resident stem cells of mesenchymal origin, these cells are morphologically and even molecularly- at least concerning the classical marker genes- indistinguishable from similar lineage cells, particularly fibroblasts.

Methods

Here we used microarray-based gene expression and global DNA methylation analyses as well as accompanying computational tools in order to specify differences between MSCs and fibroblasts, to further unravel potential identity genes and to highlight MSC signaling pathways with regard to their trophic and immunosuppressive action.

Results

We identified 1352 differentially expressed genes, of which in the MSCs there is a strong signature for e.g., KRAS signaling, known to play essential role in stemness maintenance, regulation of coagulation and complement being decisive for resolving inflammatory processes, as well as of wound healing particularly important for their regenerative capacity. Genes upregulated in fibroblasts addressed predominately transcription and biosynthetic processes and mapped morphological features of the tissue. Concerning the cellular identity, we specified the already known HOX code for MSCs, established a potential HOX code for fibroblasts, and linked certain HOX genes to functional cell-type-specific properties. Accompanied methylation profiles revealed numerous regions, especially in HOX genes, being differentially methylated, which might provide additional biomarker potential.

Discussion

Conclusively, transcriptomic together with epigenetic signatures can be successfully be used for the definition (cellular identity) of MSCs versus fibroblasts as well as for the determination of the superior functional properties of MSCs, such as their immunomodulatory potential.

Delineating the glioblastoma stemness by genes involved in cytoskeletal rearrangements and metabolic alterations

Literature data on glioblastoma ongoingly underline the link between metabolism and cancer stemness, the latter is one responsible for potentiating the resistance to treatment, inter alia due to increased invasiveness. In recent years, glioblastoma stemness research has bashfully introduced a key aspect of cytoskeletal rearrangements, whereas the impact of the cytoskeleton on invasiveness is well known. Although non-stem glioblastoma cells are less invasive than glioblastoma stem cells (GSCs), these cells also acquire stemness with greater ease if characterized as invasive cells and not tumor core cells. This suggests that glioblastoma stemness should be further investigated for any phenomena related to the cytoskeleton and metabolism, as they may provide new invasion-related insights. Previously, we proved that interplay between metabolism and cytoskeleton existed in glioblastoma. Despite searching for cytoskeleton-related processes in which the investigated genes might have been involved, not only did we stumble across the relation to metabolism but also reported genes that were found to be implicated in stemness. Thus, dedicated research on these genes in GSCs seems justifiable and might reveal novel directions and/or biomarkers that could be utilized in the future. Herein, we review the previously identified cytoskeleton/metabolism-related genes through the prism of glioblastoma stemness.

HTR1D functions as a key target of HOXA10-AS/miR-340-3p axis to promote the malignant outcome of pancreatic cancer via PI3K-AKT signaling pathway

Competing endogenous RNAs (ceRNAs) are a newly discovered class of molecular regulators involved in many diseases, especially tumors. Therefore, exploration of the potential ceRNA regulatory network regarding the occurrence and development of pancreatic cancer will provide a new theoretical basis for its diagnosis and treatment. Based on the above background, we applied a bioinformatics approach to mine the public database The Cancer Genome Atlas (TCGA) and performed a series of subsequent molecular biology assays to confirm the hypothesis that HOXA10-AS/ miR-340-3p/HTR1D axis could modulate the malignant progression of pancreatic cancer. Here, our present study demonstrated that the expression level of HTR1D, positively correlated with the level of lncRNA HOXA10-AS and negatively associated with the level of miR-340-3p, was significantly increased in pancreatic cancer cell lines (PCs) compared with that in normal HPDE6-C7 cells. Knocking down HTR1D obviously inhibited the proliferation and migration of PCs and promoted apoptosis by upregulating p-AKT. Elevated miR-340-3p blocked the progression of pancreatic cancer by downregulating HTR1D. Lessened level of lncRNA HOXA10-AS reduced the sponging of miR-340-3p, resulting in an increase of miR-340-3p and a subsequent decrease of HTR1D to ultimately suppress the malignant biological behaviors of cancer. These data illustrated that the HOXA10-AS/miR-340-3p/HTR1D ceRNA axis acted a crucial part in the malignant biological behavior of pancreatic cancer in an AKT-dependent manner.

+HOXA10-AS Promotes Malignant Phenotypes of Gastric Cancer via Upregulating HOXA10

Objective

To study the role of long noncoding RNA HOXA10-AS in gastric cancer (GC) and its underlying mechanism which is one of the most common and fetal malignancies. Long noncoding RNA HOXA10-AS is highly expressed and acts in an oncogenic role in cancers. However, its roles in GC are still unknown.

Methods

The expression of HOXA10-AS and HOXA10 in GC tissues from the TCGA database was analyzed. Western blot and qRT-PCR assays were applied to examine the expression of HOXA10-AS and HOXA10. Cell proliferation was evaluated with CCK-8 and EdU incorporation assays. Cell apoptosis was analyzed by flow cytometry. Migratory and invasive capacities were evaluated with wound healing and transwell assays.

Results

HOXA10-AS and HOXA10 were upregulated in GC, and their expressions were positively correlated. Knockdown of HOXA10-AS inhibited HOXA10 expression in GC cells. Furthermore, knockdown of HOXA10-AS restrained GC cell proliferation, migration, and invasion but promoted apoptosis. In addition, overexpression of HOXA10-AS promoted malignant phenotypes of GC cells, but all these effects could be reversed by knockdown of HOXA10.

Conclusion

HOXA10-AS promoted GC cell proliferation, migration and invasion and enhanced apoptosis via upregulating HOXA10. Our study implies a novel regulatory mechanism of malignant phenotypes and provides potential therapeutic targets for GC.

LncRNA HOXA10-AS Activated by E2F1 Facilitates Proliferation and Migration of Nasopharyngeal Carcinoma Cells Through Sponging miR-582-3p to Upregulate RAB31

Nasopharyngeal carcinoma (NPC) is a kind of head and neck cancer with a characteristic regional distribution. Increasing evidence has illustrated that long noncoding RNAs (lncRNAs) exert the regulatory function in tumor development. Nevertheless, the specific functions of lncRNA HOXA10 antisense RNA (HOXA10-AS) in NPC remain to be clarified. In this research, quantitative reverse transcription polymerase chain reaction detected HOXA10-AS expression in NPC cells. Cell counting kit-8, colony formation, and transwell assays were utilized to measure the proliferation and migration of NPC cells. Moreover, mechanism assays detected the interaction of different genes. Briefly, HOXA10-AS was highly expressed in NPC cells. HOXA10-AS down-regulation restrained NPC cell proliferation and migration. Further, HOXA10-AS could bind to miR-582-3p by acting as a competing endogenous RNA. Besides, Ras-related protein Rab-31 (RAB31) was proven as the target gene of miR-582-3p. Additionally, E2F transcription factor 1 (E2F1) acted as a transcription factor to activate HOXA10-AS expression. In the final rescue assays, we observed that the effect of HOXA10-AS depletion on NPC cell growth could be fully reversed by RAB31 overexpression or miR-582-3p inhibition. In short, our research proved that HOXA10-AS activated by E2F1 facilitated proliferation and migration of NPC cells through sponging miR-582-3p to upregulate RAB31.

A Pan-Cancer Landscape of HOX-Related lncRNAs and Their Association With Prognosis and Tumor Microenvironment

The highly conserved homology cassette family (HOX) as well as 18 referenced long non-coding antisense transcripts (HOXATs) play vital roles in the development of some cancers. Nevertheless, their expression patterns as well as their association with cancer prognosis and the tumor microenvironment (TME) in pan-cancers are still unclear. Here, based on public databases, the expression levels of HOXATs, their prognostic potentials, and correlation with tumor mutation burden (TMB), immune cell infiltration, immune subtype, immune response-related genes, and stemness scores corresponding to 33 tumor types were analyzed systematically using R language. The results of the analysis indicated that different cancer tissues show different HOXAT expression profiles. Further, HOXAT expression showed association with cancer prognosis and immune and stemness regulation. Gene set enrichment analysis also demonstrated that HOXATs participate in cancer- and immune-related pathways, and based on their expression levels, HOTAIRM1 and HOXB-AS1 showed potential involvement in oncogenesis as well as possible involvement in immune regulation across a variety of cancer types. Further investigation also confirmed a significantly higher expression of HOXB-AS1 in GBM than in lower grade glioma tissues. Importantly, in vitro cell function experiments indicated that HOXB-AS1 supports cancer stem cell and plays a fundamental role in glioma metastasis. In conclusion, our results provide valuable resources that can guide the investigation of the mechanisms related to the role of HOXATs in cancers as well as therapeutic analysis in this regard.

Construction and comprehensive analysis of a competitive endogenous RNA network to reveal potential biomarkers for the malignant differentiation of glioma

BACKGROUND

Long noncoding RNAs (lncRNAs) can act as microRNA (miRNA) sponges to regulate protein-coding gene expression; therefore, lncRNAs are considered major components of the competitive endogenous RNA (ceRNA) network and have attracted growing attention. This study explored the regulatory mechanisms and functional roles of lncRNAs as ceRNAs in the malignant differentiation of low-grade glioma (LGG) to glioblastoma (GBM) and their potential impact on the prognosis of patients with GBM.

METHODS

LncRNA and messenger RNA (mRNA) data were extracted from the Cancer Genome Atlas (TCGA) database from 156 GBM samples and 529 LGG samples. Separately, the miRNA expression data were downloaded from the Gene Expression Omnibus database, with the GSE112009 dataset containing miRNA expression data from 10 GBM samples and 15 LGG samples. Weighted gene coexpression network analysis was performed to screen the glioma grade-related lncRNAs. Then, a ceRNA network was established. The database for annotation, visualization, and integrated discovery was adopted to conduct functional enrichment analysis based on 57 upregulated differentially expressed mRNAs in the ceRNA network. Finally, Kaplan-Meier curves were created for the survival analysis of 13 hub lncRNA by combining the clinical data of GBM patients in TCGA.

RESULTS

A ceRNA network including 16 lncRNAs, 18 miRNAs, and 78 mRNAs specific to the malignant differentiation of LGG to GBM was established. The 57 upregulated differentially expressed mRNAs in the ceRNA network were significantly enriched in 35 gene ontology terms and 5 pathways. The survival analysis showed that 2 lncRNAs (LINC00261 and HOXA10-AS) were prognostic biomarkers for patients with GBM in TCGA.

CONCLUSION

The proposed ceRNA network may help elucidate the regulatory mechanism by which lncRNAs function as ceRNAs and contribute to the malignant differentiation of LGG to GBM. Importantly, the candidate lncRNAs, miRNAs, and mRNAs involved in the ceRNA network can be further evaluated as potential therapeutic targets and prognostic biomarkers for GBM.

Long non-coding RNA PSMA3-AS1 promotes glioma progression through modulating the miR-411-3p/HOXA10 pathway

BACKGROUND

Glioma is a common type of brain tumor and is classified as low and high grades according to morphology and molecules. Growing evidence has proved that long non-coding RNAs (lncRNAs) play pivotal roles in numerous tumors or diseases including glioma. Proteasome 20S subunit alpha 3 antisense RNA 1 (PSMA3-AS1), as a member of lncRNAs, has been disclosed to play a tumor-promoting role in cancer progression. However, the role of PSMA3-AS1 in glioma remains unknown. Therefore, we concentrated on researching the regulatory mechanism of PSMA3-AS1 in glioma.

METHODS

PSMA3-AS1 expression was detected using RT-qPCR. Functional assays were performed to measure the effects of PSMA3-AS1 on glioma progression. After that, ENCORI ( http://starbase.sysu.edu.cn/ ) database was used to predict potential genes that could bind to PSMA3-AS1, and miR-411-3p was chosen for further studies. The interaction among PSMA3-AS1, miR-411-3p and homeobox A10 (HOXA10) were confirmed through mechanism assays.

RESULTS

PSMA3-AS1 was verified to be up-regulated in glioma cells and promote glioma progression. Furthermore, PSMA3-AS1 could act as a competitive endogenous RNA (ceRNA) for miR-411-3p to regulate HOXA10 and thus affecting glioma progression.

CONCLUSION

PSMA3-AS1 stimulated glioma progression via the miR-411-3p/HOXA10 pathway, which might offer a novel insight for the therapy and treatment of glioma.

There and Back Again: Hox Clusters Use Both DNA Strands

Bilaterian animals operate the clusters of Hox genes through a rich repertoire of diverse mechanisms. In this review, we will summarize and analyze the accumulated data concerning long non-coding RNAs (lncRNAs) that are transcribed from sense (coding) DNA strands of Hox clusters. It was shown that antisense regulatory RNAs control the work of Hox genes in cis and trans, participate in the establishment and maintenance of the epigenetic code of Hox loci, and can even serve as a source of regulatory peptides that switch cellular energetic metabolism. Moreover, these molecules can be considered as a force that consolidates the cluster into a single whole. We will discuss the examples of antisense transcription of Hox genes in well-studied systems (cell cultures, morphogenesis of vertebrates) and bear upon some interesting examples of antisense Hox RNAs in non-model Protostomia.

PLEK2, RRM2, GCSH: A Novel WWOX-Dependent Biomarker Triad of Glioblastoma at the Crossroads of Cytoskeleton Reorganization and Metabolism Alterations

Glioblastoma is one of the deadliest human cancers. Its malignancy depends on cytoskeleton reorganization, which is related to, e.g., epithelial-to-mesenchymal transition and metastasis. The malignant phenotype of glioblastoma is also affected by the WWOX gene, which is lost in nearly a quarter of gliomas. Although the role of WWOX in the cytoskeleton rearrangement has been found in neural progenitor cells, its function as a modulator of cytoskeleton in gliomas was not investigated. Therefore, this study aimed to investigate the role of WWOX and its collaborators in cytoskeleton dynamics of glioblastoma. Methodology on RNA-seq data integrated the use of databases, bioinformatics tools, web-based platforms, and machine learning algorithm, and the obtained results were validated through microarray data. PLEK2, RRM2, and GCSH were the most relevant WWOX-dependent genes that could serve as novel biomarkers. Other genes important in the context of cytoskeleton (BMP4, CCL11, CUX2, DUSP7, FAM92B, GRIN2B, HOXA1, HOXA10, KIF20A, NF2, SPOCK1, TTR, UHRF1, and WT1), metabolism (MTHFD2), or correlation with WWOX (COL3A1, KIF20A, RNF141, and RXRG) were also discovered. For the first time, we propose that changes in WWOX expression dictate a myriad of alterations that affect both glioblastoma cytoskeleton and metabolism, rendering new therapeutic possibilities.

Comprehensive analysis of LncRNAs expression profiles in an in vitro model of steatosis treated with Exendin-4

Background and aims

The hallmark of non-alcoholic fatty liver disease (NAFLD) is the excessive hepatic lipid accumulation. Currently, no pharmacotherapy exists for NAFLD. However, the glucagon-like peptide-1 receptor agonists have recently emerged as potential therapeutics. Here, we sought to identify the long non-coding RNAs (LncRNAs) associated with the steatosis improvement induced by the GLP-1R agonist Exendin-4 (Ex-4) in vitro.

Methods

Steatosis was induced in HepG2 cells with oleic acid. The transcriptomic profiling was performed using total RNA extracted from untreated, steatotic, and Ex-4-treated steatotic cells. We validated a subset of differentially expressed LncRNAs with qRT-PCR and identified the most significantly enriched cellular functions associated with the relevant LncRNAs.

Results

We confirm that Ex-4 improves steatosis in HepG2 cells. We found 379 and 180 differentially expressed LncRNAs between untreated and steatotic cells and between steatotic and Ex-4-treated steatotic cells, respectively. Interestingly, 22 upregulated LncRNAs in steatotic cells became downregulated with Ex-4 exposure, while 50 downregulated LncRNAs in steatotic cells became upregulated in the presence of Ex-4. Although some LncRNAs, such as MALAT1, H19, and NEAT1, were previously associated with NAFLD, the association of others with steatosis and the positive effect of Ex-4 is being reported for the first time. Functional enrichment analysis identified many critical pathways, including fatty acid and pyruvate metabolism, and insulin, PPAR, Wnt, TGF-β, mTOR, VEGF, NOD-like, and Toll-like receptors signaling pathways.

Conclusion

Our results suggest that LncRNAs may play essential roles in the mechanisms underlying steatosis improvement in response to GLP-1R agonists and warrant further functional studies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12967-021-02885-4.

Silencing lncRNA HOXA10-AS decreases cell proliferation of oral cancer and HOXA10-antisense RNA can serve as a novel prognostic predictor

Objective

Long noncoding (lnc)RNAs regulate multiple biological processes including cancer. Oral squamous cell carcinoma (OSCC) is a common malignancy with poor prognosis. We aimed to identify the function of lncRNA HOXA10 antisense RNA (HOXA10-AS) and its clinical significance.

Methods

We used differential expression analysis to identify aberrantly expressed lncRNAs associated with OSCC. We identified key genes related to HOXA10-AS and their biological functions using bioinformatics tools and functional enrichment analyses. We predicted the function of HOXA10-AS using gene set enrichment and variation analyses and analyzed proliferation markers at the mRNA and protein levels. Finally, we silenced HOXA10-AS using antisense oligonucleotide and assessed proliferation ability using a cell counting kit (CCK8) and clone formation assays.

Results

In total, 506 aberrantly expressed lncRNAs were identified. HOXA10-AS was identified as a risk factor for OSCC and its expression was positively associated with tumor grade. We identified hub genes involved in regulating proliferation and predicted that HOXA10-AS is associated with an active cell cycle and increased proliferation. Silencing HOXA10-AS decreased proliferation in OSCC cell lines.

Conclusions

HOXA10-AS is involved in cell proliferation and silencing it decreases proliferation. Thus, HOXA10-AS could serve as prognostic biomarker and therapeutic target for OSCC.

FAM225B Is a Prognostic lncRNA for Patients with Recurrent Glioblastoma

Objective

The overall survival of patients with recurrent glioblastoma (rGBM) is quite different, so clinical outcome prediction is necessary to guide personalized clinical treatment for patients with rGBM. The expression level of lncRNA FAM225B was analyzed to determine its prognostic value in rGBMs.

Methods

We collected 109 samples of Chinese Glioma Genome Atlas (CGGA) RNA sequencing dataset and divided into training set and validation set. Then, we analyzed the expression of FAM225B, clinical characteristics, and overall survival (OS) information. Kaplan-Meier survival analysis was used to estimate the OS distributions. The prognostic value of FAM225B in rGBMs was tested by univariate and multivariate Cox regression analyses. Moreover, we analyzed the biological processes and signaling pathways of FAM225B.

Results

We found that FAM225B was upregulated in rGBMs (P = 0.0009). The expression of FAM225B increased with the grades of gliomas (P < 0.0001). The OS of rGBMs in the low-expression group was significantly longer than that in the high-expression group (P = 0.0041). Similar result was found in the training set (P = 0.0340) and verified in the validation set (P = 0.0292). In multivariate Cox regression analysis, FAM225B was identified to be an independent prognostic factor for rGBMs (P = 0.003). Biological process and KEGG pathway analyses implied FAM225B mainly played a functional role on transcription, regulation of transcription, cell migration, focal adhesion, etc.

Conclusions

FAM225B is expected to be as a new prognostic biomarker for the identification of rGBM patients with poor outcome. And our study provided a potential therapeutic target for rGBMs.

Transcriptional landscape of the embryonic chicken Müllerian duct

Background

Müllerian ducts are paired embryonic tubes that give rise to the female reproductive tract in vertebrates. Many disorders of female reproduction can be attributed to anomalies of Müllerian duct development. However, the molecular genetics of Müllerian duct formation is poorly understood and most disorders of duct development have unknown etiology. In this study, we describe for the first time the transcriptional landscape of the embryonic Müllerian duct, using the chicken embryo as a model system. RNA sequencing was conducted at 1 day intervals during duct formation to identify developmentally-regulated genes, validated by in situ hybridization.

Results

This analysis detected hundreds of genes specifically up-regulated during duct morphogenesis. Gene ontology and pathway analysis revealed enrichment for developmental pathways associated with cell adhesion, cell migration and proliferation, ERK and WNT signaling, and, interestingly, axonal guidance. The latter included factors linked to neuronal cell migration or axonal outgrowth, such as Ephrin B2, netrin receptor, SLIT1 and class A semaphorins. A number of transcriptional modules were identified that centred around key hub genes specifying matrix-associated signaling factors; SPOCK1, HTRA3 and ADGRD1. Several novel regulators of the WNT and TFG-β signaling pathway were identified in Müllerian ducts, including APCDD1 and DKK1, BMP3 and TGFBI. A number of novel transcription factors were also identified, including OSR1, FOXE1, PRICKLE1, TSHZ3 and SMARCA2. In addition, over 100 long non-coding RNAs (lncRNAs) were expressed during duct formation.

Conclusions

This study provides a rich resource of new candidate genes for Müllerian duct development and its disorders. It also sheds light on the molecular pathways engaged during tubulogenesis, a fundamental process in embryonic development.

Identification of methylation-driven genes related to the prognosis of papillary renal cell carcinoma: a study based on The Cancer Genome Atlas

Background

Aberrant DNA methylation patterns are involved in the pathogenesis of papillary renal cell carcinoma (pRCC). This study aimed to investigate the potential of methylation-driven genes as biomarkers in determining the prognosis of pRCC by bioinformatics analysis.

Methods

DNA methylation and transcriptome profiling data were downloaded from The Cancer Genome Atlas database. Methylation-driven genes (MDGs) were obtained using MethylMix R package. A Cox regression model was used to screen for pRCC prognosis-related MDGs, and a linear risk model based on MDG methylation profiles was constructed. A combined methylation and gene expression survival analysis was performed to further explore the prognostic value of MDGs independently.

Results

A total of 31 MDGs were obtained. Univariate and multivariate Cox regression analysis identified eight genes (CASP1, CD68, HOXD3, HHLA2, HOXD9, HOXA10-AS, TMEM71, and PLA2G16), which were used to construct a predictive model associated with overall survival in pRCC patients. Combined DNA methylation and gene expression survival analysis revealed that C19orf33, GGT6, GIPC2, HHLA2, HOXD3, HSD17B14, PLA2G16, and TMEM71 were significantly associated with patients’ survival.

Conclusion

Through the analysis of MDGs in pRCC, this study identified potential biomarkers for precision treatment and prognosis prediction, and provided the basis for future research into the molecular mechanism of pRCC.

Quantitative Analysis of Differential Expression of HOX Genes in Multiple Cancers

Transcription factors encoded by Homeobox (HOX) genes play numerous key functions during early embryonic development and differentiation. Multiple reports have shown that mis-regulation of HOX gene expression plays key roles in the development of cancers. Their expression levels in cancers tend to differ based on tissue and tumor type. Here, we performed a comprehensive analysis comparing HOX gene expression in different cancer types, obtained from The Cancer Genome Atlas (TCGA), with matched healthy tissues, obtained from Genotype-Tissue Expression (GTEx). We identified and quantified differential expression patterns that confirmed previously identified expression changes and highlighted new differential expression signatures. We discovered differential expression patterns that are in line with patient survival data. This comprehensive and quantitative analysis provides a global picture of HOX genes' differential expression patterns in different cancer types.