- 1The First Clinical Medical College of Gansu University of Chinese Medicine (Gansu Provincial Hospital), Lanzhou, China
- 2Department of Orthopedics, Gansu Provincial Hospital, Lanzhou, China
Background: Dysregulation of long non-coding (lncRNA) has been reported in various solid tumors. HOXA cluster antisense RNA 2 (HOXA-AS2) is a newly identified lncRNA with abnormal expression in several human malignancies. However, its prognostic value remains controversial. This meta-analysis synthesized available data to clarify the association between HOXA-AS2 expression levels and clinical prognosis in multiple cancers.
Methods: Four public databases (Embase, PubMed, Web of Science, The Cochrane Library) were used to identify eligible studies. Hazard ratios (HRs) and odds ratios (ORs) with their 95% confidence intervals (CIs) were combined to assess the correlation of HOXA-AS2 expression with survival outcomes and clinicopathological features of cancer patients. Publication bias was measured using Begg’s funnel plot and Egger’s regression test, and the stability of the combined results was measured using sensitivity analysis. Additionally, multiple public databases were screened and extracted to validate the results of this meta-analysis.
Results: The study included 20 studies, containing 1331 patients. The meta-analysis showed that the overexpression of HOXA-AS2 was associated with poor overall survival (HR = 2.06, 95% CI 1.58–2.69, p < 0.001). In addition, the high expression of HOXA-AS2 could forecast advanced tumor stage (OR = 3.89, 95% CI 2.90–5.21, p < 0.001), earlier lymph node metastasis (OR = 3.48, 95% CI 2.29–5.29, p < 0.001), larger tumor size (OR = 2.36, 95% CI 1.52–3.66, p < 0.001) and earlier distant metastasis (OR = 3.54, 95% CI 2.00–6.28, p < 0.001). However, other clinicopathological features, including age (OR = 1.09, 95% CI 0.86–1.38, p = 0.467), gender (OR = 0.92, 95% CI 0.72–1.18, p = 0.496), depth of invasion (OR = 2.13, 95% CI 0.77–5.90, p = 0.146) and differentiation (OR = 1.02, 95% CI 0.65–1.59, p = 0.945) were not significantly different from HOXA-AS2 expression.
Conclusion: Our study showed that the overexpression of HOXA-AS2 was related to poor overall survival and clinicopathological features. HOXA-AS2 may serve as a potential prognostic indicator and therapeutic target for tumor treatment.
Introduction
Cancer is the second greatest cause of death in most parts of the world, and it has become the most common public pathological condition on the planet (Chu et al., 2020). Traditional cancer treatments, such as surgery, adjuvant medical treatment, and actinotherapy have improved dramatically over the last century (Zhang et al., 2019). Despite this, 5-year cancer survival rates remain poor, particularly for patients with advanced tumor stage or metastasis (Li et al., 2019). One of the most significant causes is the lack of a good biomarker for detecting cancer early and predicting the clinical outcome of cancer patients (Ye et al., 2019). The significance of biomarkers in cancer has garnered increased attention in recent years across various fields, and they are thought to play critical roles in effectively screening or diagnosing cancer (Tang et al., 2020).
Long noncoding RNAs (lncRNAs) are RNA molecules that are longer than 200 nucleotides and cannot code for proteins (Zhou et al., 2019). A huge number of lncRNAs are produced during the active transcription of the human genome (Morlando and Fatica, 2018). One of the functions of lncRNAs in vivo is as tumor suppressors or oncogenes (Xu et al., 2021). Increasing evidence suggest that lncRNAs play a synergistic role in tumorigenesis or tumor suppression and that aberrant lncRNA expression is linked to cell proliferation, growth, and metastasis (Wang Y. et al., 2020). The development of RNA-targeted therapies has presented possibility of lncRNA-guided cancer therapy (Bhan et al., 2017). The inhibition of lncRNA function by RNA depletion and the removal of lncRNA exons encoding essential functional domains using splice-switching oligonucleotides may be the mechanism for targeting lncRNAs for cancer therapy (Kole et al., 2012; Schmitt and Chang, 2016). Therefore, functional lncRNA can be used as a biomarker for cancer diagnosis and for predicting treatment outcome and patient prognosis (Wang J. et al., 2020).
LncRNA HOXA cluster antisense RNA 2 (HOXA-AS2) is located on chromosome 7p15.2, a 1048-bp lncRNA, between the HOXA3 and HOXA4 genes of the HOXA cluster (Liu et al., 2019). Previous studies found that HOXA-AS2 was up-regulated in certain cancers. The increased expression of HOXA-AS2 typically predicts poor prognosis for patients with several cancers including cervical cancer (CC) (Chen and He, 2021), oral squamous cell carcinoma (OSCC) (Chen et al., 2021), lung cancer (LC) (Li and Jiang, 2017; Cui et al., 2019; Liu et al., 2019), colorectal cancer (CRC) (Li et al., 2016; Ding et al., 2017), breast cancer (BC) (Fang et al., 2017), thyroid cancer (TC) (Xia et al., 2018; Jiang et al., 2019), hepatocellular carcinoma (HCC) (Wang et al., 2016; Zhang et al., 2018; Lu et al., 2020), acute myeloid leukemia (AML) (Qu et al., 2020), bladder cancer (Wang F. et al., 2019), osteosarcoma (OSA) (Wang et al., 2018; Wang L. et al., 2019), lower-grade glioma (LGG) (Wu et al., 2019), prostate cancer (PCa) (Xiao and Song, 2020), gastric cancer (GC) (Xie et al., 2015). A high level of HOXA-AS2 expression is associated with poor overall survival (OS) and clinicopathological characteristics such as differentiation, tumor node metastasis (TNM) stage, lymph node metastasis (LNM). However, it is not clear the prognostic value of HOXA-AS2, as most of the published studies were based on a small group of patients. We explored the prognostic value of HOXA-AS2 in pan-cancer for the first time using meta-analysis. Furthermore, we further validated and explored the prognostic value of HOXA-AS2 in multiple databases through bioinformatics analysis, and explored HOXA-AS2-related genes and potential pathways. Also, the role of HOXA-AS2 in tumor immunity was investigated to identify the potential of HOXA-AS2 as a novel tumor marker and therapeutic target.
Materials and methods
Registration
The study was registered in the International Platform of Registered Systematic Review and Meta-Analysis Protocols (the registration number is: CRD42021292257). Because the present study was a systematic review and meta-analysis, Institutional Review Board (IRB) approval was not required.
Search strategy
Quality meta-analysis guidelines were followed to search for and find related papers in the Embase, PubMed, Web of Science, and The Cochrane Library. Key terms include the following: “HOXA-AS2” “long noncoding RNA HOXA-AS2” “lncRNA HOXA-AS2” “HOXA cluster antisense RNA 2” “HOXA3as” “neoplasm” “cancer” “malignancy” “neoplasia” “melanoma” “tumor” “sarcoma” “carcinoma” or “adenoma”. These terms were used to maximize the likelihood of finding a relevant article. The literature search included articles revealed as of 15 November 2021. A manual search of the reference lists of the retrieved literature was performed to confirm the eligible studies included. Any conflicts between the inclusion and exclusion clauses were resolved through group discussion.
Participants, interventions, and comparators
Studies that complied with the following criteria were eventually included: The inclusion criteria were: (a) the use of real-time quantitative polymerase chain reaction (RT-qPCR) analysis to determine the expression of HOXA-AS2 in neoplastic tissues; (b) patients diagnosed with cancer, and the study described a link between HOXA-AS2 and survival data or clinicopathology; (c) the patients were divided into two groups according to the expression level of HOXA-AS2, and (d) the quantitative hazard ratios (HRs) of OS could be extracted from the text or survival curve. The exclusion criteria were: (a) studies not related to tumors or HOXA-AS2; (b) duplicate publications; (c) reviews, conference abstracts, or case reports, and (d) studies that lacked relevant data.
Data extraction
Two researchers extracted information from each study, and any disagreement was resolved by discussing it with a third author. We obtained the following data and information from every study: (a) first author, (b) publication year, (c) country of origin, (d) cancer type, (e) number of samples, (f) HOXA-AS2 expression detection technique, (g) cut-off value, (h) sample size with high and low HOXA-AS2 expression, (i) HRs and 95% confidence intervals (CIs) for OS, (j) clinicopathologic parameters, and (k) follow-up times. OS data was directly obtained or extracted from the Kaplan-Meier (KM) curves using Engauge Digitizer version 4.1 software and the HRs and 95% CIs were computed.
Quality assessment
Two reviewers extracted information individually based on the inclusion and exclusion criteria. Some disagreements were resolved in consultation with a third reviewer. The quality of the studies was assessed using the Newcastle–Ottawa scale (NOS). The scale uses nine elements to judge a study, and a score of one is satisfied for an exact item. Total scores range from 0 to 9. A NOS score of ≥ 7 represents high-quality analysis results.
Statistical analysis
All statistical analyses were conducted using Stata software (version 12.0). The correlation of HOXA-AS2 expression with survival and clinicopathological features of tumor patients was assessed using HRs and odds ratios (ORs) with their 95% CIs, respectively. The chi-squared test and I2 statistic were preferred to determine the heterogeneity between studies. If there is strong heterogeneity (PQ < 0.1, I2 > 50%), we considered the random-effect model was applied, and the fixed-effect model was applied otherwise. All results are shown as Forest plots. Egger’s test and Begg’s funnel plot were used to evaluate publication bias, and sensitivity analysis was conducted to evaluate the robustness of the results.
Public data and tools
HOXA-AS2 expression levels in tumors and normal tissues of different solid tumors were analyzed by the Gene Expression Profiling Interactive Analysis (GEPIA, http://gepia.cancer-pku.cn) online database (based on TCGA and GTEx databases) (cutoff, p < 0.01). The survival outcomes were then verified by plotting the correlation between HOXA-AS2 expression and OS as a KM curve. Moreover, we further explored the prognostic value of HOXA-AS2 in various cancers using the Biomarker Exploration of Solid Tumors (BEST, https://rookieutopia.com) online tool. Additionally, we further explored the correlation between HOXA-AS2 and drug response through the CellMiner database (Reinhold et al., 2012) using the R-package “readxl”, “impute” and “limma” options.
Correlation of HOXA-AS2 expression with tumor immunity
Initially, we analyzed the relationship between HOXA-AS2 expression and the level of immune cell infiltration (ICI) in various cancers based on the R packages “ggExtra”, “ggpubr” and “ggplot2” via the CIBERSORT tool. Next, stromal and immune scores were calculated for each tumor sample using the ESTIMATE algorithm. The correlation between HOXA-AS2 expression and tumor microenvironment (TME) was assessed by the R package “ESTIMATE” and “limma”. In addition, the relation of HOXA-AS2 expression and tumor mutational load (TMB), tumor microsatellite instability (MSI) and immune checkpoint genes were further evaluated. TMB scores were calculated by Perl scripts, and MSI scores were determined from TCGA database mutation data. The results were visualized using the R package “RColorBrewer” and “reshape2”.
Analysis of HOXA-AS2-related genes and construction of signaling pathway network
To further investigate the value of HOXA-AS2, we obtained related genes from the MEM-Multi Experiment Matrix database. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were carried out. The top findings with a p-value of less than 0.05 were deemed significant. Finally, we used Cytoscape software to create a signal pathway network.
Results
Identification of articles
A total of 157 records were found in four electronic databases (Embase = 57, PubMed = 47, Web of Science = 53, The Cochrane Library = 0). Ninety-five duplicate articles were deleted using Endnote X9 software. After screening the titles and abstracts, 34 articles were excluded because they were not associated with the review topic or as a result of reviews, meta-analysis, letters, or expert opinions. Hence, a full-text examination was conducted for 28 articles. One article was excluded because we were unable to extract data. Two papers were eliminated because they were reviews, and five articles were eliminated because they were cell-based studies. In the end, 20 articles were included in the final meta-analysis (Figure 1).
Characteristics of the included articles
All selected articles were published between 2015 and 2021 and included 1331 patients, all of whom were from China. The smallest sample size was 27, and the largest sample size was 128. Among the twenty studies, one focused on LGG (Wu et al., 2019), one on CC (Chen and He, 2021), one on OSCC (Chen et al., 2021), three on LC (Li and Jiang, 2017; Cui et al., 2019; Liu et al., 2019), two on CRC (Li et al., 2016; Ding et al., 2017), one on BC (Fang et al., 2017), two on TC (Xia et al., 2018; Jiang et al., 2019), three on HCC (Wang et al., 2016; Zhang et al., 2018; Lu et al., 2020), two on OSA (Wang et al., 2018; Wang L. et al., 2019), one on PCa (Xiao and Song, 2020), one on GC (Xie et al., 2015), one on bladder cancer (Wang F. et al., 2019), and one on AML (Qu et al., 2020). The expression of the indicated genes in cancer tissues was measured by RT-qPCR. All eligible studies were dichotomized into low and high HOXA-AS2 expression groups based on a cut-off value. The follow-up time ranged from 25 to 120 months. All included studies were cohort studies, 65% (13/20) of which reported OS rates. The main characteristics of the eligible studies are shown in Table 1.
Association of HOXA-AS2 with OS
A total of 834 patients in 13 studies reported a link between HOXA-AS2 expression and OS. Since there was no significant heterogeneity between the studies, a fixed-effects model was used to calculate the HR and 95% CI. The pooled HR was 2.06, which indicated that HOXA-AS2 overexpression predicted poor OS in these patients with neoplasms (Figure 2A). Furthermore, KM survival analysis was applied to determine OS in different subgroups of patients according to tumor type (digestive system, respiratory system, others) (Figure 2B), sample size (n ≥ 60 or n < 60) (Figure 2C), follow-up time (≥60 months or < 60 months) (Figure 2D), and NOS score (NOS scores ≥ 7 or < 7) (Figure 2E). As depicted in Table 2, higher HOXA-AS2 expression levels were significantly associated with worse OS.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/944278/fgene-13-944278-HTML/image_m/fgene-13-944278-g002.jpg)
FIGURE 2. Relationship between HOXA-AS2 expression and overall survival. (A) Forest plots for association of HOXA-AS2 expression with overall survival. (B) Subgroup analysis stratified by cancer type. (C) Subgroup analysis stratified by sample size. (D) Subgroup analysis stratified by follow-up time. (E) Subgroup analysis stratified by NOS score.
Association between HOXA-AS2 and clinicopathologic parameters
Correlations between HOXA-AS2 expression and the clinicopathological features of the patients are shown in Table 3. The meta-analysis results showed that higher HOXA-AS2 expression levels tended to be significantly associated with advanced tumor stage (OR = 3.89, 95% CI 2.90–5.21, p < 0.001) (Figure 3C), earlier LNM (OR = 3.48, 95% CI 2.29–5.29, p < 0.001) (Figure 3D), larger tumor size (OR = 2.36, 95% CI 1.52–3.66, p < 0.001) (Figure 3E) and earlier distant metastasis (OR = 3.54, 95% CI 2.00–6.28, p < 0.001) (Figure 3H). However, age (OR = 1.09, 95% CI 0.86–1.38, p = 0.467) (Figure 3A), gender (OR = 0.92, 95% CI 0.72–1.18, p = 0.496) (Figure 3B), depth of invasion (OR = 2.13, 95% CI 0.77–5.90, p = 0.146) (Figure 3G), and differentiation (OR = 1.02, 95% CI 0.65–1.59, p = 0.945) (Figure 3F), had no significant link with increased HOXA-AS2 expression levels.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/944278/fgene-13-944278-HTML/image_m/fgene-13-944278-g003.jpg)
FIGURE 3. Forest plots for association of HOXA-AS2 expression with clinicopathological features. (A) Age. (B) Gender. (C) Tumor stage. (D) Lymph node metastasis. (E) Tumor size. (F) Differentiation. (G) Depth of invasion. (H) Distant metastasis.
Publication bias and sensitivity analysis
We employ Begg’s funnel plot and Egger’s regression test to identify publication bias of OS. The shape of Begg’s funnel was essentially symmetrical, with no visible asymmetry (Figure 4A), and Egger’s regression analysis did not reveal the presence of publication bias (Pr > |t| = 0.971). We ran a sensitivity analysis by eliminating one qualified study to analyze the influence of a single study on the result. According to the analysis, the results were not significantly influenced (Figure 4B). This verifies the reliability of the meta-analysis conclusions.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/944278/fgene-13-944278-HTML/image_m/fgene-13-944278-g004.jpg)
FIGURE 4. (A) Begg’s funnel plot of HOXA-AS2 for overall survival. (B) A sensitivity analysis of pooled HR for overall survival.
Validation of HOXA-AS2 expression in public databases
We used the TCGA dataset to analyze the degree of HOXA AS2 expression in various tumors to further corroborate our results. HOXA-AS2 was aberrantly expressed in glioblastoma multiforme (GMB), acute myeloid leukemia (LAML), pancreatic adenocarcinoma (PAAD), and thymoma (THYM), compared to normal controls (Supplementary Figure S1A). A violin plot revealed that the degree of HOXA-AS2 expression in human cancer was highly related to the clinical stage (Supplementary Figure S1B). We used GEPIA to create survival graphs by combining HOXA-AS2 expression data with the OS data of patients with malignancies in the entire TCGA dataset, which included 9491patients separated into high (4741) and low (4750) groups of HOXA-AS2 expression based on median levels. The results showed that increased HOXA-AS2 expression predicted poor OS, confirming the meta-analysis results (Supplementary Figure S1C). Additionally, we explored the link of HOXA-AS2 expression and tumor prognosis using Cox regression model through the BEST online tool. The findings revealed that there was a significant correlation between HOXA-AS2 expression and the prognosis of GBM, stomach adenocarcinoma (STAD), LGG, adrenocortical carcinoma (ACC), skin cutaneous melanoma (SKCM), breast invasive carcinoma (BRCA), lung adenocarcinoma (LUAD), bladder urothelial carcinoma (BLCA), sarcoma (SARC), and CRC in at least two datasets (Figure 5).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/944278/fgene-13-944278-HTML/image_m/fgene-13-944278-g005.jpg)
FIGURE 5. The correlation between HOXA-AS2 expression and survival in different tumor types was analyzed in BEST. Cox regression analysis of (A) GBM, (B) STAD, (C) LGG, (D) ACC, (E) SKCM, (F) BRCA, (G) LUAD, (H) BLCA, (I) SARC, (J) CRC. *p < 0.05, **p < 0.01, and ***p < 0.001.
HOXA-AS2 and drug response
To further explore the significance of HOXA-AS2 in guiding cancer treatment, we analyzed the relationship between HOXA-AS2 and drug response. The results revealed that patients with high HOXA-AS2 expression had a better drug response to XL−147, Cpd−401, cordycepin, fenretinide, estramustine, and arsenic trioxide. In contrast, AS−703569, ENMD−2076, SB−1317, benzaldehyde (BEN), staurosporine, aminoflavone, amonafide, midostaurin, sapitinib, and KW−2449 had a better drug response in patients with low HOXA-AS2 expression (Figure 6).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/944278/fgene-13-944278-HTML/image_m/fgene-13-944278-g006.jpg)
FIGURE 6. An illustration of the association between HOXA-AS2 expression and expected medication response.
Correlation analysis of HOXA-AS2 expression with tumor immunity
Correlation analysis between HOXA-AS2 expression and ICI levels identified remarkable correlations between HOXA-AS2 expression and ICI levels in BRCA (n = 7), kidney renal papillary cell carcinoma (KIRP) (n = 6), kidney renal clear cell carcinoma (KIRC) (n = 5), head and Neck squamous cell carcinoma (HNSC) (n = 4), STAD (n = 4), thyroid carcinoma (THCA) (n = 3), testicular germ cell tumors (TGCT) (n = 3), BLCA (n = 3), LGG (n = 3), THYM (n = 2), ovarian serous cystadenocarcinoma (OV) (n = 2), esophageal carcinoma (ESCA) (n = 2), prostate adenocarcinoma (PRAD) (n = 1), mesothelioma (MESO) (n = 1), lung squamous cell carcinoma (LUSC) (n = 1), LAML (n = 1), and cervical squamous cell carcinoma and endocervical adenocarcinoma (CESC) (n=1). Detailed information on the subpopulations of infiltrating immune cells in various cancer types is illustrated in Figure 7. HOXA-AS2 expression was negatively correlated with the levels of infiltrating M0 macrophages in ESCA, NHSC, CESC, THCA, BRCA, LUSC, and BLCA (Figure 7A). Similarly, HOXA-AS2 expression was negatively associated with the levels of infiltrating neutrophils in STAD, HNSC, KIRC, and BRCA (Figure 7B). In regard to monocytes, their infiltration levels were negatively correlated with the HOXA-AS2 expression in LGG (Figure 7C). HOXA-AS2 expression was also negatively related to the levels of infiltrating M1 macrophages in BLCA, and KIRP (Figure 7D). Moreover, HOXA-AS2 expression was negatively correlated with the levels of infiltrating M2 macrophages in BRCA, but positively associated with TGCT (Figure 7E). HOXA-AS2 expression was positively associated with the levels of infiltrating naive B cells in STAD, OV, and BRCA, but positively associated with KIRP and TGCT (Figure 7F). HOXA-AS2 expression was negatively correlated with the infiltrating levels of activated CD4 memory T cells in BLCA, KIRC, KIRP, STAD, TGCT, and THYM (Figure 7G). Furthermore, HOXA-AS2 expression presented a positive relationship with the levels of infiltrating plasma cells in BRCA (Figure 7H). The levels of infiltrating CD8 T cells were positively associated with HOXA-AS2 expression in PRAD, HNSC, and BRCA, but negatively related in KIRP (Figure 7I). A positive association with infiltrating follicular helper T cells was identified in KIRC (Figure 7J). The levels of infiltrating resting CD4 memory T cells were positively associated with HOXA-AS2 expression in KIRP, and LGG (Figure 7K). HOXA-AS2 expression was negatively associated with the levels of infiltrating activated dendritic cells in MESO, HNSC, KIRC, ESCA, and THCA (Figure 7L). Moreover, HOXA-AS2 expression was positively associated with the levels of infiltrating resting DCs in BRCA (Figure 7M). The levels of infiltrating activated mast cells were negatively associated with HOXA-AS2 expression in LGG (Figure 7N). In contrast, the levels of infiltrating resting mast cells were positively correlated with HOXA-AS2 expression in THYM, STAD, and KIRP, but negatively related in KIRC and LAML (Figure 7O). The levels of infiltrating memory B cells were positively associated with HOXA-AS2 expression in THCA, but negatively associated in OV (Figure 7P).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/944278/fgene-13-944278-HTML/image_m/fgene-13-944278-g007.jpg)
FIGURE 7. Correlation between HOXA-AS2 gene expression and the level of immune cell infiltration in pan-cancerous tissues. HOXA-AS2 expression significantly correlated with infiltrating levels of M0 macrophages in BLCA, BRCA, CESC, ESCA, NHSC, LUSC, and THCA (A), neutrophils in BRCA, HNSC, KIRC, and STAD (B), monocytes in LGG (C), M1 macrophages in BLCA and KIRP (D), M2 macrophages in BRCA and TGCT (E), naive B cells in BRCA, KIRP, OV, STAD, and TGCT (F), CD4 memory T cells in BLCA, KIRC, KIRP, STAD, TGCT, and THYM (G), plasma cells in BRCA (H), CD8 T cells in BRCA, HNSC, KIRP, and PRAD (I), follicular helper T cells in KIRC (J), resting CD4 memory T cells in KIRP and LGG (K), activated dendritic cells in ESCA, HNSC, KIRC, MESO, and THCA (L), resting DCs in BRCA (M), activated mast cells in LGG (N), resting mast cells in KIRC, KIRP, LAML, STAD, and THYM (O), memory B cells in OV and THCA (P).
Furthermore, to explore its relationship with the TME, we analyzed the association of HOXA-AS2 expression with stromal and immune scores. Our findings showed that HOXA-AS2 expression correlated with the stromal scores of 13 cancers, the top 6 tumors were BRCA, pheochromocytoma and paraganglioma (PCPG), THCA, LGG, TGCT, and LUAD (Figure 8A); and with the immune scores of 11 cancers, the top 6 tumors were PRAD, LGG, BRCA, THCA, PCPG, and LUAD (Figure 8B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/944278/fgene-13-944278-HTML/image_m/fgene-13-944278-g008.jpg)
FIGURE 8. Correlations between HOXA-AS2 expression and immunity, including stromal score, immune score, immune checkpoint genes, TMB and MSI in cancers. (A) Correlation of HOXA-AS2 expression with the stromal score in pan-cancer. (B) Correlation of HOXA-AS2 expression with the immune score in pan-cancer. (C) Correlation of HOXA-AS2 expression and immune checkpoint genes. (D) The radar chart displays the correlation of TMB with HOXA-AS2 expression. The red curve indicates the correlation coefficient, and the blue value indicates the range. (E) The radar chart displays the correlation of MSI with HOXA-AS2 expression. The blue curve indicates the correlation coefficient, and the green value indicates the range. *p < 0.05, **p < 0.01, and ***p < 0.001.
Moreover, the correlation of HOXA-AS2 expression with immune checkpoint genes showed that CD44, CD40, VSIR, LGALS9, TNFRSF14 and TNFRSF25 were significantly associated with HOXA-AS2 expression in several cancers, especially in LGG, PRAD, LUAD, THCA, and HNSC (Figure 8C). We evaluated the relation between HOXA-AS2 expression and TMB/MSI as well, and found that there was a significant positive correlation between its expression and TMB in LGG, THYM, KIRC, and HNSC, while there was a significant negative correlation in STAD, BRCA, colon adenocarcinoma (COAD), PRAD, SKCM, BLCA, uterine carcinosarcoma (UCS), and PCPG (Figure 8D). Also, there was a significant positive correlation between HOXA-AS2 expression and MSI in kidney chromophobe (KICH) and BRCA, while there was a significant negative correlation in COAD and STAD (Figure 8E).
Analysis of HOXA-AS2-related genes
The top 150 genes were screened co-expressed by HOXA-AS2 using the MEM-Multi Experiment Matrix database. HOXA2, HOXA5, and HOXA3 were the top three co-expressed genes ranked by p-value and interrelated with HOXA-AS2 expression (Supplementary Figure S2). In addition, GO and KEGG pathway analyses were performed to explore the underlying molecular mechanisms (Figure 9, Table 4), and a signaling pathway network was constructed using Cytoscape software (Figure 10).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/944278/fgene-13-944278-HTML/image_m/fgene-13-944278-g009.jpg)
FIGURE 9. GO terms and the KEGG pathway. (A) GO enrichment of target genes in biological process ontolagy (p < 0.05).(B) GO enichment of target genes in cellular component ontology (p < 0.05). (C) GO enrichment of target genes in molecular function ontology (p < 0.05).(D) The top 8 pathways related to the differentially expressed genes by the KEGG database analysis. BP, biological process; CC, cellular component GO, gene ontology analysis; KEGG, Kyoto Encyclopedia of Genes and Genomes; MF, molecular function.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/944278/fgene-13-944278-HTML/image_m/fgene-13-944278-g010.jpg)
FIGURE 10. Differentially expressed gene interaction network analysis. Green nodes represent target genes and sky blue nodes represent the related pathway. As indicated in red, HOXA-AS2 localized at the center of the network.
Discussion
According to global cancer statistics 2020, there were approximately 19.3 million new cancer cases and 10.0 million cancer deaths worldwide (Sung et al., 2021). Despite the variety of treatments available today, cancer has a high rate of recurrence and death, resulting in increased costs and poor patient prognosis (Huang et al., 2020). For most types of cancer, early detection and treatment improve prognoses. LncRNA was recently found to have huge clinical value in the early diagnosis and novel treatment of patients with cancer (Dunn et al., 2010). Numerous studies showed that lncRNAs were involved in a variety of physiological and pathological processes and had important effects on tumorigenesis and tumor growth (Renganathan and Felley-Bosco, 2017). For instance, prostate cancer-associated transcript 6 (PCAT6) was significantly increased in various cancers. The overexpression of PCAT6 was closely correlated with OS, TNM stage, distant metastasis, LNM, tumor size, and the degree of differentiation in cancer patients (Shi et al., 2021), which may be a new cancer-related biomarker. In recent years, numerous studies found that the HOXA-AS2 was overexpressed in a variety of tumor types, such as CC (Chen and He, 2021), OSCC (Chen et al., 2021), PCa (Xiao and Song, 2020) and so on. The function of HOXA-AS2 in many malignancies has yet to be realized. Thus, we conducted the current meta-analysis to assess the predictive and clinical importance of HOXA-AS2 aberrant expression in cancer patients.
Our meta-analysis found that HOXA-AS2 overexpression was related to a lower chance of survival. We also examined the connection between HOXA-AS2 overexpression and several clinicopathological features. HOXA-AS2 overexpression was associated with tumor stage, LNM, tumor size and distant metastasis. However, there were no significant associations with age, gender, differentiation, or depth of invasion. Furthermore, we performed further analysis of the prognostic role of HOXA-AS2 in various cancers using several public databases. Among them, Cox regression analysis indicated that HOXA-AS2 had a better prognostic value in GBM, STAD, LGG, ACC, SKCM, BRCA, LUAD, BLCA, SARC, and CRC. These results suggest that HOXA-AS2 is a predictor of poor prognosis in cancer patients.
Our findings demonstrated a strong correlation between HOXA-AS2 expression and immunity in multiple tumors. Tumor-infiltrating immune cells play an irreplaceable role in tumor development. Previous research has shown that HOXA-AS2 could influence glioma progression by regulating Treg cell proliferation and immune tolerance (Zhong et al., 2022). We found that HOXA-AS2 was associated with multiple infiltrating immune cells in a variety of tumors. However, the regulatory mechanism of HOXA-AS2 on ICI still remains to be further confirmed by abundant experiments. The interaction between TME and tumor cells is decisive for tumor survival and progression. Immune cells and stromal cells are key components of the TME (Xiao and Yu, 2021). We found that HOXA-AS2 expression correlated with immune cell scores in 11 tumors and with stromal cell scores in 13 tumors. It indicates that HOXA-AS2 has an essential role in the TME. TMB and MSI are well-directed for tumor immunotherapy. Our findings indicated that HOXA-AS2 expression correlated with TMB in 12 cancers and with MSI in 4 cancers. In summary, HOXA-AS2 may further influence the prognosis of cancer patients via modulation of tumor immunity.
Although HOXA-AS2 was demonstrated to be a major predictive factor for patients with various malignancies in several studies, the basic principle of how HOXA-AS2 caused cancer remains unknown. According to the results of this study, the overexpression of HOXA-AS2 can significantly aid cancer growth and metastasis. In contrast, inhibited HOXA-AS2 expression, significantly reduced cell proliferation, migration, and invasion, as well as the carcinogenesis process. In PCa, HOXA-AS2 exhibited a negative connection with miR-509–3p. The inhibition of HOXA-AS2 prevented PCa cells from proliferating and migrating (Xiao and Song, 2020). This suggested that HOXA-AS2 could be used as a therapeutic target to treat PCa. In addition, by decreasing miR-520c-3p expression, HOXA-AS2 enhanced the growth and spread of HCC (Wang et al., 2016). HOXA-AS2 was elevated in OSCC tissues and increased OSCC cell proliferation by sponging miR-567/CDK8 (Chen et al., 2021). Table 5 highlights the association between HOXA-AS2 and malignancies to investigate functionally associated genes.
We used the MEM-Multi Experiment Matrix database to predict target genes and perform the signaling pathway analysis of HOXA-AS2 to further investigate its value. HOXA2, HOXA5, and HOXA3, all of which play important roles in cancer, were strongly associated with HOXA-AS2 expression in our study. Following that, we conducted GO analysis, which indicated that the sequence-specific DNA binding, extracellular exosomes, and angiogenesis of HOXA-AS2 were all significantly related. HOXA AS2 was highly connected to cancer-associated pathways in KEGG analysis.
Nevertheless, this meta-analysis had several limitations. First, some HRs and the corresponding 95% CIs were extracted from KM curves. Second, the qualifying studies were all performed in China, so it is unclear whether the results can be generalized to patients in other countries. To address this limitation, we validated the correlation between HOXA-AS2 expression and prognosis of cancer patients in public databases. Third, the included studies were inconsistent in dividing expression according to cut-off values. Additionally, only a few trials were included, and some cancer types had very low sample sizes. Thus, more clinical investigations should be conducted to assess the potential prognostic role of HOXA-AS2 expression in cancer types that were not included.
Conclusion
In summary, this meta-analysis found that HOXA-AS2 overexpression was linked to the poor prognosis of cancer patients and could be used as a new prognostic biomarker and therapeutic target for various malignancies. The predictive usefulness of HOXA-AS2 in tumors has to be confirmed in more studies, including other cancer types.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Materials, further inquiries can be directed to the corresponding authors.
Author contributions
XL and WL conceived this study. FZ, MG, and XP were responsible for the collection, extraction, and analysis of the data. FZ, GZ, and HZ was responsible for writing the paper. DZ, XX, and FC performed the quality evaluation and completed data analysis. XL and WL reviewed the paper. All authors read and approved the final manuscript.
Funding
This work was financially supported by Natural Science Foundation of Gansu Province (No. 20JR10RA365); Key R and D Program of Gansu Province-Social Development (No. 20YF8FA092).
Acknowledgments
We would like to thank all the researchers and study participants for their contributions.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2022.944278/full#supplementary-material
References
Bhan, A., Soleimani, M., and Mandal, S. S. (2017). Long noncoding RNA and cancer: A new paradigm. Cancer Res. 77 (15), 3965–3981. doi:10.1158/0008-5472.Can-16-2634
Chen, R., and He, P. (2021). Long noncoding RNA HOXA-AS2 accelerates cervical cancer by the miR-509-3p/BTN3A1 axis. J. Pharm. Pharmacol. 73 (10), 1387–1396. doi:10.1093/jpp/rgab090
Chen, R., Wang, X., Zhou, S., and Zeng, Z. (2021). LncRNA HOXA-AS2 promotes tumor progression by suppressing miR-567 expression in oral squamous cell carcinoma. Cancer Manag. Res. 13, 5443–5455. doi:10.2147/cmar.S305946
Chu, H. Y., Chen, Y. J., Hsu, C. J., Liu, Y. W., Chiou, J. F., Lu, L. S., et al. (2020). Physical cues in the microenvironment regulate stemness-dependent homing of breast cancer cells. Cancers (Basel) 12 (8), E2176. doi:10.3390/cancers12082176
Cui, T. J., Lin, G. S., Dai, Y. M., Zheng, J. P., Chen, Z., Chen, Q., et al. (2019). LncRNA HOXA-AS2 regulates microRNA-216a-5p to promote malignant progression of non-small cell lung cancer. Eur. Rev. Med. Pharmacol. Sci. 23 (3), 264–273. doi:10.26355/eurrev_201908_18656
Ding, J., Xie, M., Lian, Y., Zhu, Y., Peng, P., Wang, J., et al. (2017). Long noncoding RNA HOXA-AS2 represses P21 and KLF2 expression transcription by binding with EZH2, LSD1 in colorectal cancer. Oncogenesis 6 (1), e288. doi:10.1038/oncsis.2016.84
Dunn, B. K., Wagner, P. D., Anderson, D., and Greenwald, P. (2010). Molecular markers for early detection. Semin. Oncol. 37 (3), 224–242. doi:10.1053/j.seminoncol.2010.05.007
Fang, Y., Wang, J., Wu, F., Song, Y., Zhao, S., and Zhang, Q. (2017). Long non-coding RNA HOXA-AS2 promotes proliferation and invasion of breast cancer by acting as a miR-520c-3p sponge. Oncotarget 8 (28), 46090–46103. doi:10.18632/oncotarget.17552
Huang, S., Yang, J., Fong, S., and Zhao, Q. (2020). Artificial intelligence in cancer diagnosis and prognosis: Opportunities and challenges. Cancer Lett. 471, 61–71. doi:10.1016/j.canlet.2019.12.007
Jiang, L., Wu, Z., Meng, X., Chu, X., Huang, H., and Xu, C. (2019). LncRNA HOXA-AS2 facilitates tumorigenesis and progression of papillary thyroid cancer by modulating the miR-15a-5p/HOXA3 Axis. Hum. Gene Ther. 30 (5), 618–631. doi:10.1089/hum.2018.109
Kole, R., Krainer, A. R., and Altman, S. (2012). RNA therapeutics: Beyond RNA interference and antisense oligonucleotides. Nat. Rev. Drug Discov. 11 (2), 125–140. doi:10.1038/nrd3625
Li, Q., Dai, Y., Wang, F., and Hou, S. (2016). Differentially expressed long non-coding RNAs and the prognostic potential in colorectal cancer. Neoplasma 63 (6), 977–983. doi:10.4149/neo_2016_617
Li, W., Li, J., Mu, H., Guo, M., and Deng, H. (2019). MiR-503 suppresses cell proliferation and invasion of gastric cancer by targeting HMGA2 and inactivating WNT signaling pathway. Cancer Cell. Int. 19, 164. doi:10.1186/s12935-019-0875-1
Li, Y., and Jiang, H. (2017). Up-regulation of long non-coding RNA HOXA-AS2 in non-small cell lung cancer is associated with worse survival outcome. Int. J. Clin. Exp. Pathol. 10 (9), 9690–9696.
Liu, Y., Lin, X., Zhou, S., Zhang, P., Shao, G., and Yang, Z. (2019). Long noncoding RNA HOXA-AS2 promotes non-small cell lung cancer progression by regulating miR-520a-3p. Biosci. Rep. 39 (5), BSR20190283. doi:10.1042/bsr20190283
Lu, Q., Gao, J., Tang, S., Li, Z., Wang, X., Deng, C., et al. (2020). Integrated RNA sequencing and single-cell mass cytometry reveal a novel role of LncRNA HOXA-AS2 in tumorigenesis and stemness of hepatocellular carcinoma. Onco. Targets. Ther. 13, 10901–10916. doi:10.2147/OTT.S272717
Morlando, M., and Fatica, A. (2018). Alteration of epigenetic regulation by long noncoding RNAs in cancer. Int. J. Mol. Sci. 19 (2), E570. doi:10.3390/ijms19020570
Qu, Y., Wang, Y., Wang, P., Lin, N., Yan, X., and Li, Y. (2020). Overexpression of long noncoding RNA HOXA-AS2 predicts an adverse prognosis and promotes tumorigenesis via SOX4/PI3K/AKT pathway in acute myeloid leukemia. Cell. Biol. Int. 44 (8), 1745–1759. doi:10.1002/cbin.11370
Reinhold, W. C., Sunshine, M., Liu, H., Varma, S., Kohn, K. W., Morris, J., et al. (2012). CellMiner: A web-based suite of genomic and pharmacologic tools to explore transcript and drug patterns in the NCI-60 cell line set. Cancer Res. 72 (14), 3499–3511. doi:10.1158/0008-5472.Can-12-1370
Renganathan, A., and Felley-Bosco, E. (2017). Long noncoding RNAs in cancer and therapeutic potential. Adv. Exp. Med. Biol. 1008, 199–222. doi:10.1007/978-981-10-5203-3_7
Schmitt, A. M., and Chang, H. Y. (2016). Long noncoding RNAs in cancer pathways. Cancer Cell. 29 (4), 452–463. doi:10.1016/j.ccell.2016.03.010
Shi, S. B., Cheng, Q. H., Gong, S. Y., Lu, T. T., Guo, S. F., Song, S. M., et al. (2021). PCAT6 may be a new prognostic biomarker in various cancers: A meta-analysis and bioinformatics analysis. Cancer Cell. Int. 21 (1), 370. doi:10.1186/s12935-021-02079-4
Sung, H., Ferlay, J., Siegel, R. L., Laversanne, M., Soerjomataram, I., Jemal, A., et al. (2021). Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. Ca. Cancer J. Clin. 71 (3), 209–249. doi:10.3322/caac.21660
Tang, L., Chen, Y. L., Tang, X., Wei, D., Xu, X. Y., and Yan, F. (2020). Long noncoding RNA DCST1-AS1 promotes cell proliferation and metastasis in triple-negative breast cancer by forming a positive regulatory loop with miR-873-5p and MYC. J. Cancer 11 (2), 311–323. doi:10.7150/jca.33982
Wang, F., Wu, D., Chen, J., Chen, S., He, F., Fu, H., et al. (2019a). Long non-coding RNA HOXA-AS2 promotes the migration, invasion and stemness of bladder cancer via regulating miR-125b/Smad2 axis. Exp. Cell. Res. 375 (1), 1–10. doi:10.1016/j.yexcr.2018.11.005
Wang, F., Yang, H., Deng, Z., Su, Y., Fang, Q., and Yin, Z. (2016). HOX antisense lincRNA HOXA-AS2 promotes tumorigenesis of hepatocellular carcinoma. Cell. Physiol. biochem. 40 (1-2), 287–296. doi:10.1159/000452545
Wang, J., Wang, W., Tang, Q., Lu, L., Luo, Z., Li, W., et al. (2020a). Long non-coding RNA lnc-GNAT1-1 suppresses liver cancer progression via modulation of epithelial-mesenchymal transition. Front. Genet. 11, 1029. doi:10.3389/fgene.2020.01029
Wang, L., Wang, L., and Zhang, X. (2019b). Knockdown of lncRNA HOXA-AS2 inhibits viability, migration and invasion of osteosarcoma cells by miR-124-3p/E2F3. Onco. Targets. Ther. 12, 10851–10861. doi:10.2147/ott.S220072
Wang, Y., Yang, C., Liu, X., Zheng, J., Zhang, F., Wang, D., et al. (2020b). Transcription factor AP-4 (TFAP4)-upstream ORF coding 66 aa inhibits the malignant behaviors of glioma cells by suppressing the TFAP4/long noncoding RNA 00520/microRNA-520f-3p feedback loop. Cancer Sci. 111 (3), 891–906. doi:10.1111/cas.14308
Wang, Y., Zhang, R., Cheng, G., Xu, R., and Han, X. (2018). Long non-coding RNA HOXA-AS2 promotes migration and invasion by acting as a ceRNA of miR-520c-3p in osteosarcoma cells. Cell. Cycle 17 (13), 1637–1648. doi:10.1080/15384101.2018.1489174
Wu, L., Zhu, X., Song, Z., Chen, D., Guo, M., Liang, J., et al. (2019). Long non-coding RNA HOXA-AS2 enhances the malignant biological behaviors in glioma by epigenetically regulating RND3 expression. Onco. Targets. Ther. 12, 9407–9419. doi:10.2147/ott.S225678
Xia, F., Chen, Y., Jiang, B., Du, X., Peng, Y., Wang, W., et al. (2018). Long noncoding RNA HOXA-AS2 promotes papillary thyroid cancer progression by regulating miR-520c-3p/S100A4 pathway. Cell. Physiol. biochem. 50 (5), 1659–1672. doi:10.1159/000494786
Xiao, S., and Song, B. (2020). LncRNA HOXA-AS2 promotes the progression of prostate cancer via targeting miR-509-3p/PBX3 axis. Biosci. Rep. 40 (8), BSR20193287. doi:10.1042/bsr20193287
Xiao, Y., and Yu, D. (2021). Tumor microenvironment as a therapeutic target in cancer. Pharmacol. Ther. 221, 107753. doi:10.1016/j.pharmthera.2020.107753
Xie, M., Sun, M., Zhu, Y. N., Xia, R., Liu, Y. W., Ding, J., et al. (2015). Long noncoding RNA HOXA-AS2 promotes gastric cancer proliferation by epigenetically silencing P21/PLK3/DDIT3 expression. Oncotarget 6 (32), 33587–33601. doi:10.18632/oncotarget.5599
Xu, X., Zhong, Z., Shao, Y., and Yi, Y. (2021). Prognostic value of MEG3 and its correlation with immune infiltrates in gliomas. Front. Genet. 12, 679097. doi:10.3389/fgene.2021.679097
Ye, J., Sun, H., Feng, Z., Zhang, Q., Xia, Y., Ji, Y., et al. (2019). Prognostic significance of LncRNA GHET1 expression in various cancers: A systematic review and meta-analysis. Biosci. Rep. 39 (10), BSR20190608. doi:10.1042/bsr20190608
Zhang, C., Ren, X., He, J., Wang, W., Tu, C., and Li, Z. (2019). The prognostic value of long noncoding RNA SNHG16 on clinical outcomes in human cancers: A systematic review and meta-analysis. Cancer Cell. Int. 19, 261. doi:10.1186/s12935-019-0971-2
Zhang, Y., Xu, J., Zhang, S., An, J., Zhang, J., Huang, J., et al. (2018). HOXA-AS2 promotes proliferation and induces epithelial-mesenchymal transition via the miR-520c-3p/GPC3 Axis in hepatocellular carcinoma. Cell. Physiol. biochem. 50 (6), 2124–2138. doi:10.1159/000495056
Zhong, C., Tao, B., Li, X., Xiang, W., Peng, L., Peng, T., et al. (2022). HOXA-AS2 contributes to regulatory T cell proliferation and immune tolerance in glioma through the miR-302a/KDM2A/JAG1 axis. Cell. Death Dis. 13 (2), 160. doi:10.1038/s41419-021-04471-4
Keywords: lncRNA, HOXA-AS2, cancers, prognosis, meta-analysis, bioinformatics analysis
Citation: Zhang F, Zhang G, Zhang H, Pu X, Chi F, Zhang D, Xin X, Gao M, Luo W and Li X (2022) HOXA-AS2 may be a potential prognostic biomarker in human cancers: A meta-analysis and bioinformatics analysis. Front. Genet. 13:944278. doi: 10.3389/fgene.2022.944278
Received: 15 May 2022; Accepted: 27 October 2022;
Published: 10 November 2022.
Edited by:
Yilin Zhang, The University of Chicago, United StatesReviewed by:
Qi Shengcai, Shanghai Stomatology Prevention Hospital, ChinaXiaoqi Yang, Huazhong University of Science and Technology, China
Copyright © 2022 Zhang, Zhang, Zhang, Pu, Chi, Zhang, Xin, Gao, Luo and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xingyong Li, gssylxy@163.com; Wenyuan Luo, cf15104026385@163.com
†These authors share first authorship