Research Paper Volume 16, Issue 9 pp 7752—7773
Long non-coding RNA HOTAIR: from pan-cancer analysis to colorectal cancer-related uridine metabolism
- 1 Department of Gastroenterology, The Affiliated Hospital of Yangzhou University, Yangzhou University, Yangzhou, China
- 2 Department of Oncology, The Second Affiliated Hospital of Nanjing Medical University, Nanjing, China
Received: September 26, 2023 Accepted: March 29, 2024 Published: May 1, 2024
https://doi.org/10.18632/aging.205781How to Cite
Copyright: © 2024 Chen et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Abstract
Long non-coding RNAs (lncRNAs) are involved significantly in the development of human cancers. lncRNA HOTAIR has been reported to play an oncogenic role in many human cancers. Its specific regulatory role is still elusive. And it might have enormous potential to interpret the malignant progression of tumors in a broader perspective, that is, in pan-cancer. We comprehensively investigated the effect of HOTAIR expression on tumor prognosis across human malignancies by analyzing multiple cancer-related databases like The Cancer Genome Atlas (TCGA) and Tumor Immune Estimation Resource (TIMER). Bioinformatics data indicated that HOTAIR was overexpressed in most of these human malignancies and was significantly associated with the prognosis of patients with cancer, especially in colorectal cancer (CRC). Subsequently, this study further clarified the utility of HOTAIR that downregulation of its expression could result in reduced proliferation and invasion of CRC cells. Mechanistically, HOTAIR upregulated the metabolic enzymes UPP1 by recruiting histone methyltransferase EZH2, thereby increasing the tumor progression. Our results highlight the essential role of HOTAIR in pan-cancer and uridine bypass, suggesting that the HOTAIR/EZH2/UPP1 axis might be a novel target for overcoming CRC. We anticipate that the role of HOTAIR in metabolism could be important in the context of CRC and even exploited for therapeutic purposes.
Introduction
Cancer is a major public health problem with high global mortality and is the second leading cause of death worldwide. Despite tremendous breakthroughs in the field of tumors now, most cancer remains incurable [1]. Although the overall mortality rate has been effectively controlled due to the development of pathogenesis, diagnosis and treatment, the understanding of its biological characteristics still does not meet the requirements of translational medicine, with the result that the annual incidence rate continues to increase [2, 3]. Given the increased mortality in cancer patients, there is a need to further explore the mechanisms of cancer to improve prognosis and level of treatment. Through public data analysis and basic experimental exploration, we can now have a more comprehensive understanding of the important role of certain genes in human cancer progression.
Long non-coding RNA (lncRNAs) is an RNA molecule that is longer than 200 nucleotides and does not have the ability to encode proteins, playing a central role in different cancers through various potential molecular mechanism [4]. HOTAIR is located on chromosome 12 and is transcribed from the HOXC locus [5]. According to current research, HOTAIR regulates downstream gene expression by recruitment of chromatin modifiers and either by competitive binding with translation factor at the transcriptional and/or posttranscriptional levels [6, 7]. HOTAIR is an oncogenic RNA whose expression level is associated with many clinical features such as tumor grade and prognosis and has been explored to be oncogenic in colorectal cancer (CRC). HOTAIR also regulates various downstream protein expression via various mechanism axis. Enhancer of zeste homolog 2 (EZH2), a histone methyltransferase, regulate downstream target genes expression via specifically inducing the methylation on Lysine-27 of histone 3 (H3K27me) in multiple cancers [8, 9]. Uridine phosphorylase 1 (UPP1) functions in its phosphorolysis of uridine to uracil and ribose-1-phosphate, which is considered as an indispensable metabolic enzyme in the pyrimidine salvage biological processes and is upregulated in many cancers, including CRC [10]. When glucose is limiting, uridine might serve as alternative sources fulfil energy requirements. Uridine bypass makes CRC cells capable of using uridine and promoting tumor growth [11].
In conclusion, we conducted a comprehensive summary about the biological functions of HOTAIR with pan-cancers. The results show that HOTAIR is highly expressed in most tumors indicating that HOTAIR is linked to cancer development and prognosis. Furthermore, the expression level of HOTAIR also can reflect the immune infiltration level of tumor tissues. In CRC, HOTAIR upregulated the metabolic enzymes UPP1 by recruiting EZH2, thereby promoting the tumor progression, suggesting a novel metabolic axis for CRC therapy.
Materials and Methods
Data collection
We used Gent2 database (http://gent2.appex.kr/gent2/), TIMER database (http://timer.cistrome.org/), Gene Expression Profiling Interactive Analysis (GEPIA) database (http://gepia2.cancer-pku.cn/#index) to compare the expression of HOTAIR between tumor tissues and normal tissues among 33 types of cancers. We obtained the overall survival data to investigate the association between the expression level of HOTAIR and survival status by Kaplan–Meier Plotter across various tumors in the GEPIA database and PrognoScan database. The cutoff value was set as cutoff-high (50%) and cutoff-low (50%) in distinguishing the high or low expression of HOTAIR. cBioPortal tool was used to analyze the condition of alteration frequency, mutated site information and mutation pattern of HOTAIR across all tumors. The infiltration data were used to demonstrate if there was a link between HOTAIR expression and infiltration in TIMER database. We used the TIMER and GEPIA database to explore the co-expression genes of HOTAIR, then performed gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis in HOTAIR and neighborhood genes.
Gent2 database
We used Gent2 database to explore the expression level of HOTAIR in pan-cancer level [12]. The database can reflect the genes expression in pan-cancer.
TIMER database
TIMER2 is a database for analysis of immune infiltrates in various cancer types. TIMER2 provides immune infiltrates’ abundances estimated from the TCGA database to explore the immunological characteristics of tumors. The infiltration data were used to demonstrate if there was a link between the expression level of HOTAIR and infiltration [13].
GEPIA database
GEPIA2 is a tool containing differential expression gene data and clinical data to assess the prognostic value of specific genes. Patient cases were divided into high expression and low expression depending on the gene expression level [14]. Hypothesis testing used the log-rank test.
PrognoScan database survival
We got the overall survival data to investigate the clinical links between HOTAIR expression level and survival status by the PrognoScan database [15].
The cBioPortal database
cBioPortal tool (14) was used to analyze the condition of alteration frequency, mutated site information and mutation pattern of HPOTAIR across TCGA tumors [16].
LinkedOmics database
The LinkedOmics is a database that includes 32 TCGA cancer cohorts [17]. In our study, LinkedOmics database was used to explore genes differentially expressed with HOTAIR in the TCGA COAD cohort. The results were represented by volcano plots and heat maps.
Gene enrichment analysis
We obtained the top 200 genes with a co-expression pattern to HOTAIR in all TCGA tumor tissue data in the “Similar Genes Detection” module of GEPIA2.
Clinical tissue data
Twenty-five pairs of tumor and normal tissue samples were collected from patients with CRC. Tissue samples obtained by surgical treatment are frozen in liquid nitrogen. None of the patients had received specific treatment. CRC diagnosis was confirmed via pathological examination. All the informed consent forms of patients are signed and ethical approval for this experiment was approved by the Nanjing Medical University, Medical Ethics Committee.
Cell culture, transfection and qRT- PCR
The COAD cell lines (SW480, SW620, HCT-116, DLD1, LOVO and HT-29) and the normal cell line (NCM460) were all obtained from Nanjing medical university laboratory. All cell lines were authenticated by STR profiling and were routinely tested for mycoplasma contamination. Cells were cultured in suitable culture media added with 1% penicillin/streptomycin and 10% fetal bovine serum (FBS). siRNAs (HIPPOBIO Biotechnology, Nanjing, China) of HOTAIR transfected in the SW480 and HCT-116 cell lines, followed the manufacturer’s instructions. Gene knockdown was achieved by transfecting cells with specific siRNAs or shRNA using Lipofectamine 3000 (Invitrogen, USA). The total RNA was extracted using TRIzol reagent (Invitrogen, USA). The expression levels of target genes and GAPDH were explored by qRT-PCR utilizing the SYBR qPCR (Vazyme, Nanjing, China). The primers were listed in Supplementary Data 1.
Plasmid construction
Three UPP1 over-expressing plasmid and shRNAs targeting UPP1 were designed and purchased from GeneCopoeia Company (USA). They were transfected into CRC cells with lipo3000 (Invitrogen, USA) following the manufacturer’s manual. The sequences of all shRNAs were summarized in Supplementary Data 1.
Cell proliferation assays
First, SW480 (1500 cells/well) and HCT-116 (1000 cells/well) cells were seeded in 96-well plates after cell transfection. CCK-8 reagent (Vazyme, Nanjing, China) was added into each well and incubated at 37°C for 2 h, the absorbance was measured at 450 nm using a microplate reader. CRC cells were planted into 6-well plates to perform colony formation assays. After a 2-week incubation period, we used 4% paraformaldehyde to fix the cells which were stained with crystal violet.
Transwell assay
CRC cells were seeded in the plate upper chambers. The bottom chamber was added with full medium supplemented with 10% FBS. Following incubation at 24 hours, we used 4% paraformaldehyde fixing the cells on the lower surface and crystal violet staining the cells. Subsequently, cells were photographed under a microscope.
ChIP assay and RIP assay
Chromatin immunoprecipitation (ChIP) assay and RNA immunoprecipitation (RIP) assay was conducted following manufacturer’s instruction (Millipore, USA). The sequences of primers for the UPP1 promoter were summarized in Supplementary Data 1.
Western blot assay
We extracted total protein from CRC cells using RIPA buffer added with protease inhibitor. SDS–PAGE was used to separate proteins which was transferred onto PVDF membranes. Specific antibodies: UPP1 (Abcam, UK; ab128854), EZH2 (Proteintech, China; 21800-1-AP), Actin antibody was used as a control.
Statistical analysis
All experiments were conducted using GraphPad Prism 8, and the results were expressed as means ± SD derived from at least three independent samples. To assess differences, an unpaired two-tailed t-test, as recommended for independent analysis, was employed. Statistical significance was determined by a p-value below 0.05.
Results
Expression levels of HOTAIR in the pan-cancer analysis
We conducted a comparative analysis to assess HOTAIR expression in pan-cancer by comparing it between tumor and normal samples by the Gent2 database (Figure 1A). Compared with normal tissues, HOTAIR expression level was upregulated in various tumors, including colorectal cancer, breast cancer, brain cancer, kidney cancer, lung cancer, pancreatic cancer, skin cancer, stomach cancer and uterus cancer (P < 0.05). In addition, HOTAIR expression level was decreased in skin cancer (P < 0.05). Second, we analyzed the different expression level in HOTAIR using the TIMER database (Figure 1B). Figure 1B showed that overexpression of HOTAIR was observed in 17 types of tumors including colon adenocarcinoma (COAD), cholangiocarcinoma (CHOL), breast invasive carcinoma (BRCA), esophageal carcinoma (ESCA), glioblastoma multiforme (GBM), head and neck squamous cell carcinoma (HNSC), kidney renal clear cell carcinoma (KIRC), kidney renal papillary cell carcinoma (KIRP), liver hepatocellular carcinoma (LIHC), lung adenocarcinoma (LUAD), lung squamous cell carcinoma (LUSC), Pheochromocytoma and Paraganglioma (PCPG), rectum adenocarcinoma (READ), stomach adenocarcinoma (STAD) and thyroid carcinoma (THCA). However, the expression level of HOTAIR was reduced in Kidney Chromophobe (KICH) (all P < 0.05). Furthermore, we utilized the GEPIA dataset to investigate HOTAIR expression across human tumors (Supplementary Figure 1). Our findings revealed that HOTAIR expression level was higher in BRCA compared to neighboring non-cancerous tissues (P < 0.05; Supplementary Figure 2).
Associations between the expression of HOTAIR and clinicopathological features progression
We analyzed the expression level of HOTAIR at different pathology stages via the GEPIA database (Supplementary Figure 3). The HOTAIR expression varied significantly in COAD and READ (P < 0.05). We further explored the correlation between the RNA expression level and patients’ clinicopathological parameters in various cancers. We found a positive correlation between the age, gender, race and tumor stage and expression level of HOTAIR at the same time in ACC, BRCA, COAD, KIRC, MESO and THYM (Figure 2). These results showed that HOTAIR expression level could impact the prognosis in patients. To analyse the diagnostic value of HOTAIR, this study calculated and plotted the diagnostic receiver operating characteristic (ROC). We used the data from TCGA dataset plotting the diagnostic ROC for pan-cancer analysis (Supplementary Figure 4). Results showed that the area under the curve (AUC) of HOTAIR diagnostic ROC was 0.794 in BRCA, 0.829 in KICH, 0.661 in KIRC, 0.705 in LUAD, 0.847 in LUSC, 0.759 in PRAD, 0.934 in STAD, 0.648 in COAD, and 0.846 in PAAD indicating that HOTAIR may be a promising diagnostic biomarker in human cancers. Interestingly, we found in the results that the ROC curve was still well represented in the CRC. This result clearly reconfirmed the role of HOTAIR in clinical practice.
Figure 2. Correlation between HOTAIR expression level and clinicopathological parameters in human cancers including age, gender, race and stage across all tumors in TCGA. The red color indicates a positive correlation (0–1), while the blue color represents a negative correlation (−1–0). The correlation with P-value < 0.05 is considered as statistically significant. Statistically non-significant correlations values are marked with a cross.
Prognostic analysis of HOTAIR in human cancers
We conducted a prognosis analysis using data obtained from PrognoScan and GEPIA to examine the relationship between HOTAIR expression levels and overall survival (OS) across various cancer types. In PrognoScan database, the results revealed that HOTAIR expression associated with prognosis in colorectal and breast cancers. (Supplementary Figure 5A–5D). Furthermore, using the GEPIA database, we found that high HOTAIR expression in adrenocortical carcinoma (ACC), COAD and KIRC was associated with shorter OS and DFS (Figure 3A, 3B), indicating a detrimental effect of HOTAIR overexpression in these cancer types. Conversely, in GBM, mesothelioma (MESO), cervical squamous cell carcinoma and endocervical adenocarcinoma (CESC), LUAD and Pancreatic adenocarcinoma (PAAD), overexpression of HOTAIR was associated with longer OS, suggesting a potentially protective role of HOTAIR in these specific cancer contexts. These findings highlight the complexity and background dependence of the role of HOTAIR in tumorigenesis and progression, whose function exhibits completely opposite properties in different tumor types. The internal mechanisms controlling these different effects are likely to be extremely complex and require further investigation.
Figure 3. Relationship between HOTAIR expression level and patient survival in TCGA tumors. Relationship between HOTAIR gene expression and survival overall (A), disease-free survival (B) was assessed in all TCGA tumors using GEPIA2. The positive results of survival map and Kaplan-Meier curves are listed.
Genetic alterations of HOTAIR
We used the cBioPortal online database to conduct an analysis of genetic alterations in HOTAIR across TCGA types. We explored that HOTAIR was altered in 2.7% patients including amplification and deep deletion (Supplementary Figure 6A). Among the alterations, the most common type was amplification. The most prevalent genetic alteration in HOTAIR was the “amplification” type, observed in nearly all cancer types. (Supplementary Figure 6B). Notably, the highest frequency of HOTAIR alterations (10%) was observed in PAAD patients, primarily in the form of “amplification”. In addition, “deep deletion” of HOTAIR occurred in 1% of patients with BRCA, and mutations in HOTAIR were the deep deletion type in 0.5% of patients with LIHC. Next, we concluded the relationship between the types of the alterations and the levels of HOTAIR mRNA. As shown in Supplementary Figure 6C, among the alterations the level of HOTAIR mRNA increased from the deletion group, amplification group, gain group to the diploid group in a sequential manner.
HOTAIR associated with the CNV events across human cancers
The progression of tumors is inseparable from the immune system, and the rise of immunotherapy is changing the paradigm of cancer treatment. Therefore, to examining the possibility of that HOTAIR might be explored as an immune target in cancer therapy, we took advantage of the TIMER database to investigate the association between the expression level of HOTAIR and immune infiltration by six main immune cells including B cells, CD8+ T cells, CD4+ T cells, macrophages, neutrophils, and dendritic cells (Figure 4). Among different cancers, HOTAIR CNV is linked to the infiltrating level of different types of immune cells: B cells, CD8+ T cells, CD4+ T cells, macrophages, neutrophils and dendritic cells in the BRCA, B cells and neutrophils in COAD, B cells, CD4+ T cells, macrophages and dendritic cells in the LUAD and B cells, CD4+ T cells, neutrophils and dendritic cells in the PAAD. Thus, HOTAIR may affect the immunity of different cancers and do so through potential distinct mechanisms.
Immune infiltration analysis
We explored the correlations between HOTAIR and immune infiltration. As shown in Figure 5A, in GBMLGG, LGG and PRAD, there is a positively correlation between the immune score and HOTAIR expression. Conversely, there is a negative correlation between the immune score and HOTAIR expression in SARC, STES and STAD. These results indicate a tight correlation between HOTAIR expression levels and tumor immune responsiveness in most tumor types (Figure 5B). In COAD, HOTAIR positively correlated with four immune cell types (B cells, neutrophils, CD4+ T cells and macrophages). From the single-cell level, HOTAIR expression level was most highly in malignant cells (Figure 5C). Figure 5D shows the close connection between high HOTAIR expression with relative percentage of cells.
Figure 5. The correlation of HOTAIR expression with immune infiltration levels. (A) The correlation between HOTAIR and immune score. (B) The correlation between HOTAIR and immune cells. (C) Cells from the GSE148673 dataset were mapped on the tSNE plot. (D) tSNE plot illustrating HOTAIR expression profile at the cell level.
Gene enrichment analysis of HOTAIR
To elucidate the functional mechanism of HOTAIR in carcinogenesis, we employed the GEPIA tool to identify genes that exhibited co-expression with HOTAIR. A heatmap is presented to illustrate the positive correlation between HOTAIR and these 5 genes (HOXC10, HOXC11, HOXC12, HOXC13, and NKD2) (Supplementary Figure 7A). HOTAIR expression correlated positively with HOXC10 (homeobox C10) (R = 0.38), HOXC11 (homeobox C11) (R = 0.6), HOXC12 (homeobox C12) (R = 0.66), HOXC13 (homeobox C13) (R = 0.46) and NKD2 (NKD inhibitor of WNT signaling pathway 2) (R = 0.6) (all p < 0.001) (Supplementary Figure 7B). Furthermore, gene pathway enrichment analysis revealed that the biological processes involving the top HOTAIR-associated genes were closely associated with protein modifications-related pathways (Supplementary Figure 7C). The co-expression genes of HOTAIR were analyzed by LinkedOmics database to further explore the biological effects of HOTAIR. As shown in Supplementary Figure 7D, the expression level of HOTAIR affects the expression of a large number of genes. Genes positively (Supplementary Figure 7E) and negatively (Supplementary Figure 7F) both correlated with HOTAIR are shown in the heat map.
LINC00467 promotes CRC cell proliferation and migration
Differential analysis and prognostic analysis showed that HOTAIR played an important biological role in CRC. By analyzing the TCGA dataset, we summarized the correlation between HOTAIR expression levels and clinical indicators such as age and stage (p < 0.05) (Table 1). First, we examined the expression level of HOTAIR in related CRC cell lines by RT-qPCR and HOTAIR expression was highly expressed in CRC cells compared with the normal cell line NCM460 (Figure 6A). Thus, we selected SW480 and HCT116 cell lines for further research. We also confirmed the knockdown efficacy of siRNAs (si1 and si2) targeting HOTAIR through RT-qPCR. Compared with adjacent normal tissues, LINC00467 expression level was upregulated in tumor tissue samples (Figure 6B). The proliferation assay (CCK-8 and colony assays) proved that the low expression of HOTAIR significantly suppressed cell proliferation (Figure 6C, 6D). Our study also demonstrated that HOTAIR siRNA suppressed the cell migration abilities in both SW480 and HCT116 cell lines (Figure 6E).
Table 1. Correlations between HOTAIR expression and clinicopathological characteristics in TCGA-COAD.
Characteristic | Low expression of HOTAIR | High expression of HOTAIR | p |
n | 322 | 322 | |
T stage, n (%) | 0.180 | ||
T1 | 11 (1.7%) | 9 (1.4%) | |
T2 | 65 (10.1%) | 46 (7.2%) | |
T3 | 210 (32.8%) | 226 (35.3%) | |
T4 | 33 (5.1%) | 41 (6.4%) | |
N stage, n (%) | 0.004 | ||
N0 | 201 (31.4%) | 167 (26.1%) | |
N1 | 75 (11.7%) | 78 (12.2%) | |
N2 | 44 (6.9%) | 75 (11.7%) | |
M stage, n (%) | 0.103 | ||
M0 | 240 (42.6%) | 235 (41.7%) | |
M1 | 36 (6.4%) | 53 (9.4%) | |
Age, median (IQR) | 66 (57, 74) | 69 (59, 77) | 0.034 |
Figure 6. Effect of HOTAIR in CRC cell proliferation and migration. (A) Transcription level of HOTAIR in CRC cell lines and NCM460. (B) Transcription level of HOTAIR in HCT116 and SW480 were significantly downregulated by si-HOTAIR transfection, respectively. (C) Relative expression of HOTAIR in CRC tissues was assessed by comparing with normal samples (n = 25). (D) Cell proliferation was assessed by CCK assay. (E) Transwell assay employed to detect the migration ability of HOTAIR knockdown cells. *p < 0.05; **p < 0.01; ***p < 0.001.
HOTAIR affects CRC progression through uridine bypass via EZH2/UPP1 axis
To further explore the mechanism of how HOTAIR regulates colorectal cancer proliferation and migration, we performed an RNA sequencing analysis (Figure 7A). Among the large number of differentially expressed genes, we selected the top 10 genes with the highest expression for experimental verification (Figure 7B). In the qRT-PCR analysis, we found the highest expression of UPP1 gene in the case of HOTAIR knockdown. Then we designed the HOTAIR expression vector and three shRNAs to study its biological role (Figure 7C). Studies have shown that lncRNAs can regulate the downstream genes expression, such as EZH2, Ago2, LSD1, and SUZ12, by interacting with RNA-binding proteins [18, 19]. So, we carried out RIP experiment to verify and screen and the results showed the strongest binding between HOTAIR and EZH2, indicating that HOTAIR interacted with EZH2 specifically (Figure 7D). Furthermore, we observed that UPP1 was down-regulated at both mRNA and protein levels when EZH2 was knocked down, indicating that EZH2 promoted the high expression of UPP1 (Figure 7E, 7F). Chromatin immunoprecipitation (Chip) assay showed that EZH2 directly bound to the UPP1 promoter region and induced H3K27 trimethylation. Meanwhile, knockdown of HOTAIR reduced EZH2 and H3K27 trimethylation levels through the UPP1 promoter (Figure 7G). These results showed that HOTAIR promotes CRC progression through the upregulate of UPP1 via interacting with EZH2. The expression level of UPP1 is significantly increased in CRC samples (Figure 7H). UPP1 is a uridine phosphorylase catalysing the phosphate-dependent catabolism of uridine into R1P and uracil (Figure 7I). Glucose oxidation drives cellular bioenergetics. Recent studies have confirmed that UPP1 liberates uridine-derived ribose to fuel tumor metabolism and thereby support cells proliferation [20]. This phenomenon depends on uridine being present and expression of UPP1 (Figure 7J). Colony-forming assays have also demonstrated that HOTAIR could affect CRC progression through uridine bypass via EZH2/UPP1 axis (Figure 7K).
Figure 7. HOTAIR affects CRC proliferation through uridine bypass via EZH2/UPP1 axis. (A) A hierarchically clustered heatmap of differentially expressed genes in CRC cells after transfection of si-HOTAIR or NC-siRNAs. (B) Ten representative genes expression levels in CRC cells depleted of HOTAIR. (C) The knockdown and over-expressing efficiency of UPP1 were detected by RT-qPCR. (D) RIP assays were performed in CRC cells to show HOTAIR co-immunoprecipitation with EZH2, LSD1, Ago2 and SUZ12. (E, F) EZH2 and UPP1 mRNA and proteins level in CRC cells transfected with si-EZH2 or siRNA-NC by qRT-PCR and western blot analysis. (G) EZH2 occupancy on the UPP1 promoters was upregulated by HOTAIR knockdown by ChIP-qPCR assay. (H) UPP1 was highly expressed in CRC tissues compared with normal tissues through TCGA. (I) Reaction catalysed by UPP1 proteins. (J) Cell growth assays of HCT116 control cells in glucose-free media in the presence of 10 mM of either glucose, uridine or ribose. (K) Colony-forming assays assessed CRC cell proliferation in the presence of si-HOTAIR, sh-UPP1, OE-UPP1, glucose, uridine or ribose.
Discussion
HOTAIR was discovered by Howard Chang’s group that repressors the expression of homeobox gene D cluster (HOXD) by recruiting a transcriptional corepressor [21]. The human HOTAIR gene is located in the intergenic region between HOXC11 and HOXC12 in the HOXC gene cluster on chromosome 12 [22]. Since the first appearance of HOTAIR in the field of oncology, its highly expressed role has been considered as an oncogene [23–26]. For example, HOTAIR contributes to 5FU Resistance through suppressing miR-218 and activating NF-κB/TS signaling in colorectal cancer [27]. Consistent with the results of our study, HOTAIR was highly expressed in cancer tissues and positively correlated with cancer proliferation and migration in multiple tumor cohorts. Further studies with a larger sample size will help to draw reliable conclusions about the expression of HOTAIR in tumor tissues. Although more and more studies have shown that HOTAIR has important clinical implications in a variety of tumors [28–30], Up to now, the role of HOTAIR in pan-cancer has not been reported. Therefore, we need to better understand the role of HOTAIR in pan-cancer and its potential prognostic value, as well as the molecular mechanism of its action.
We found that HOTAIR genes are highly expressed in a total of seventeen cancers and lowly expressed in one cancer by TIMER database. The differences in samples led to the differences in the results of different databases. Sample differences lead to differences in the results of different databases, but as long as there are enough samples, the cancer-promoting role of HOTAIR is still plausible. Despite great progress has been achieved in recent years, the mortality of CRC is still increasing [31]. Effective prognostic prediction is of great significance to improve the survival of patients with CRC. However, to date, prognostic biomarkers remain limited. High expression of HOTAIR in CRC was identified, and a high expression of HOTAIR was associated with poor survival in patients with CRC. Interestingly, the pan-cancer analysis also revealed a markedly poor survival in tumors with high HOTAIR expression, which indicated that HOTAIR might participate in tumorigenesis. Through genetic alteration analysis of HOTAIR, we observed that HOTAIR produced significant alterations in a variety of cancers. Therefore, we can hypothesize that HOTAIR expression and mutation are potential parts of tumor biology.
The relationship between immunity cells and cancer cells is very close and immunotherapy represents the future of clinical cancer treatment [32].
Given the important role of HOTAIR in a variety of tumors, we further explored the correlation between tumor immunity and HOTAIR. Accumulating evidence suggests that HOTAIR is closely related to the proliferation and metastasis of human cancers [33, 34]. However, according to the current research, there are few studies on the relationship between HOTAIR expression level and immune infiltration.
We observed a significant correlation between HOTAIR expression and tumor immune cell infiltration, including B cells, CD8+ T cells, CD4+ T cells, macrophages, neutrophils, and dendritic cells. This suggests that HOTAIR may impact tumor processing and prognosis through its influence on cancer immunity. In this study, we provide the first statistical evidence of a positive association between HOTAIR expression and the infiltration level of CD4+ T cells in COAD. Pre-clinical and clinical investigations have identified intra-tumoral cytotoxic CD4+ T cells capable of directly eliminating cancerous cells [35]. Studies have demonstrated that tumor-derived exosomes facilitate the immunosuppressive function of B-regs by delivering HOTAIR in COAD. Tumor-derived HOTAIR directs B cells towards a regulatory phenotype characterized by programmed cell death-ligand 1 (PDL1) expression in CRC, thereby inducing PDL1-expressing B cells to suppress CD8+T cell activity [36]. The pivotal role of HOTAIR in both immunity and tumorigenesis necessitates further investigations to elucidate its biological functions and underlying mechanisms.
In our study, further investigation confirmed a robust positive correlation between HOTAIR expression and the HOXC cluster genes. The HOTAIR gene is transcribed in an antisense orientation relative to its flanking HOXC11 and HOXC12 genes. Previous studies have demonstrated that increased DNA methylation in an intergenic CpG island located between HOXC12 and HOTAIR is positively associated with HOTAIR expression [37]. Additionally, we observed that HOTAIR and its interacting genes primarily focus on “distal pattern formation”, “histone modification”, “covalent chromatin modification”, etc.
We then validated the function of HOTAIR in COAD and observed significant overexpression of HOTAIR, which was consistent with current research, using bioinformatics analysis. Furthermore, we explored the oncogenic role of HOTAIR through cell experiments and found that it potentiates cancer metastasis and tumor progression by recruiting EZH2 [18]. Similar to previous studies, we identified HOTAIR as an independent prognostic factor in CRC. Specifically, blocking HOTAIR with siRNA significantly suppressed CRC cancer cell growth and migration in vitro. In future studies, we aim to confirm the function of HOTAIR both in vivo and in vitro using more advanced molecular biology techniques. Glucose serves as a source of energy. When glucose is limiting, alternative nutrients such as undine could fuel ATP production and gluconeogenesis [38]. HOTAIR could affects CRC progression through uridine bypass via EZH2/UPP1 axis which upregulated UPP1 transcription level and uridine catabolism in CRC cells. We propose that a new signaling axis controlling UPP1 transcription thus maintaining the homeostasis of energy and biosynthesis.
While our study has improved our understanding of HOTAIR from a pan-cancer perspective, there are still several limitations. Firstly, the sample size of cancer patients in the TCGA database was significantly larger than that of normal patients. Secondly, the absence of clinical factors in the public database, such as specific details regarding medication and/or surgical treatment, also impacts patient prognosis. Additionally, our study does have a limitation concerning the association between HOTAIR and immunotherapy. The involvement of HOTAIR in immunotherapy should be further validated through clinical and cellular experiments, such as co-culturing tumor cells with immune cells using interference for HOTAIR. Lastly, this present study is retrospective, and prospective studies should be conducted to address the limitations inherent to retrospective research design. Despite these limitations, our study provides valuable insights into investigating the function of HOTAIR in cancers and identifies potential targets and prognostic markers for CRC treatment. In conclusion, our pan-cancer analysis offers a comprehensive overview of the oncogenic roles played by HOTAIR across various human cancers. Overexpression of HOTAIR generally indicates poor prognosis for cancer patients. Our study characterizes HOTAIR expression patterns across different cancer types and highlights its potential value as a predictive biomarker while shedding light on further exploration into its prognostic and therapeutic potential. Growing evidence that HOTAIR plays crucial roles in carcinogenesis. In the study, HOTAIR was demonstrated to promote tumorigenesis via recruiting EZH2 in CRC cells, indicating the clinical importance of HOTAIR/EZH2/UPP1 axis as a promising therapeutic target for CRC. Recently, the rapid progress in mRNA vaccine design and delivery technologies has significantly expedited the development and clinical application of mRNA-based cancer vaccines. Future investigations could potentially explore synergistic combinations of HOTAIR cancer vaccines with immunotherapeutic approaches to enhance the efficacy of cancer treatment.
Abbreviations
lncRNAs: long non-coding RNAs; TIMER: Tumor Immune Estimation Resource; GEPIA: Gene Expression Profiling Interactive Analysis; TCGA: The Cancer Genome Atlas; GEO: The Gene Expression Omnibus; KEGG: Kyoto Encyclopedia of Genes and Genomes; GO: gene ontology; GTEx: Genotype-Tissue Expression; OS: overall survival; DFS: disease-free survival; CPTAC: Clinical Proteomics Tumor Analysis Consortium; FBS: fetal bovine serum; BRCA: breast invasive carcinoma; CHOL: cholangiocarcinoma; COAD: colon adenocarcinoma; ESCA: esophageal carcinoma; GBM: glioblastoma multiforme; HNSC: head and neck squamous cell carcinoma; KIRP: kidney renal clear cell carcinoma; KIRC: kidney renal papillary cell carcinoma; LIHC: liver hepatocellular carcinoma; LUAD: lung adenocarcinoma; LUSC: lung squamous cell carcinoma; PCPG: pheochromocytoma and paraganglioma; READ: rectum adenocarcinoma; STAD: stomach adenocarcinoma; THCA: thyroid carcinoma; KICH: kidney chromophobe; MESO: mesothelioma; CESC: cervical squamous cell carcinoma and endocervical adenocarcinoma; PAAD: pancreatic adenocarcinoma; ROC: receiver operating characteristic; AUC: area under the curve; PRC2: polycomb repressive complex 2; HOXD: homeobox gene D cluster; PDL1: programmed cell death-ligand 1.
Author Contributions
All of the authors worked collaboratively on the work presented here. XC, SW, XJ and YD designed the experiments and supervised the study. SW and MZ searched the articles and made figures; XC wrote this manuscript. All authors read and approved the final manuscript.
Conflicts of Interest
The authors declare no conflicts of interest related to this study.
Ethical Statement and Consent
All the informed consent forms of patients are signed and ethical approval for this experiment was given by the Nanjing Medical University, Medical Ethics Committee (No. [2019]-KY-121).
Funding
This work was supported by grants from the National Natural Science Foundation of China (No. 82273084), the Strengthening Health Care via Science and Education Project.
References
- 1. Qi J, Li M, Wang L, Hu Y, Liu W, Long Z, Zhou Z, Yin P, Zhou M. National and subnational trends in cancer burden in China, 2005-20: an analysis of national mortality surveillance data. Lancet Public Health. 2023; 8:e943–55. https://doi.org/10.1016/S2468-2667(23)00211-6 [PubMed]
- 2. Fitzgerald RC, Antoniou AC, Fruk L, Rosenfeld N. The future of early cancer detection. Nat Med. 2022; 28:666–77. https://doi.org/10.1038/s41591-022-01746-x [PubMed]
- 3. Siegel RL, Giaquinto AN, Jemal A. Cancer statistics, 2024. CA Cancer J Clin. 2024; 74:12–49. https://doi.org/10.3322/caac.21820 [PubMed]
- 4. Nojima T, Proudfoot NJ. Mechanisms of lncRNA biogenesis as revealed by nascent transcriptomics. Nat Rev Mol Cell Biol. 2022; 23:389–406. https://doi.org/10.1038/s41580-021-00447-6 [PubMed]
- 5. Woo CJ, Kingston RE. HOTAIR lifts noncoding RNAs to new levels. Cell. 2007; 129:1257–9. https://doi.org/10.1016/j.cell.2007.06.014 [PubMed]
- 6. Marchese FP, Huarte M. Long non-coding RNAs and chromatin modifiers: their place in the epigenetic code. Epigenetics. 2014; 9:21–6. https://doi.org/10.4161/epi.27472 [PubMed]
- 7. Smillie CL, Sirey T, Ponting CP. Complexities of post-transcriptional regulation and the modeling of ceRNA crosstalk. Crit Rev Biochem Mol Biol. 2018; 53:231–45. https://doi.org/10.1080/10409238.2018.1447542 [PubMed]
- 8. Duan R, Du W, Guo W. EZH2: a novel target for cancer treatment. J Hematol Oncol. 2020; 13:104. https://doi.org/10.1186/s13045-020-00937-8 [PubMed]
- 9. Wang J, Yu X, Gong W, Liu X, Park KS, Ma A, Tsai YH, Shen Y, Onikubo T, Pi WC, Allison DF, Liu J, Chen WY, et al. EZH2 noncanonically binds cMyc and p300 through a cryptic transactivation domain to mediate gene activation and promote oncogenesis. Nat Cell Biol. 2022; 24:384–99. https://doi.org/10.1038/s41556-022-00850-x [PubMed]
- 10. Caradoc-Davies TT, Cutfield SM, Lamont IL, Cutfield JF. Crystal structures of Escherichia coli uridine phosphorylase in two native and three complexed forms reveal basis of substrate specificity, induced conformational changes and influence of potassium. J Mol Biol. 2004; 337:337–54. https://doi.org/10.1016/j.jmb.2004.01.039 [PubMed]
- 11. Skinner OS, Blanco-Fernández J, Goodman RP, Kawakami A, Shen H, Kemény LV, Joesch-Cohen L, Rees MG, Roth JA, Fisher DE, Mootha VK, Jourdain AA. Salvage of ribose from uridine or RNA supports glycolysis in nutrient-limited conditions. Nat Metab. 2023; 5:765–76. https://doi.org/10.1038/s42255-023-00774-2 [PubMed]
- 12. Park SJ, Yoon BH, Kim SK, Kim SY. GENT2: an updated gene expression database for normal and tumor tissues. BMC Med Genomics. 2019; 12:101. https://doi.org/10.1186/s12920-019-0514-7 [PubMed]
- 13. Li T, Fan J, Wang B, Traugh N, Chen Q, Liu JS, Li B, Liu XS. TIMER: A Web Server for Comprehensive Analysis of Tumor-Infiltrating Immune Cells. Cancer Res. 2017; 77:e108–10. https://doi.org/10.1158/0008-5472.CAN-17-0307 [PubMed]
- 14. Tang Z, Li C, Kang B, Gao G, Li C, Zhang Z. GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. 2017; 45:W98–102. https://doi.org/10.1093/nar/gkx247 [PubMed]
- 15. Mizuno H, Kitada K, Nakai K, Sarai A. PrognoScan: a new database for meta-analysis of the prognostic value of genes. BMC Med Genomics. 2009; 2:18. https://doi.org/10.1186/1755-8794-2-18 [PubMed]
- 16. Cerami E, Gao J, Dogrusoz U, Gross BE, Sumer SO, Aksoy BA, Jacobsen A, Byrne CJ, Heuer ML, Larsson E, Antipin Y, Reva B, Goldberg AP, et al. The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. Cancer Discov. 2012; 2:401–4. https://doi.org/10.1158/2159-8290.CD-12-0095 [PubMed]
- 17. Vasaikar SV, Straub P, Wang J, Zhang B. LinkedOmics: analyzing multi-omics data within and across 32 cancer types. Nucleic Acids Res. 2018; 46:D956–63. https://doi.org/10.1093/nar/gkx1090 [PubMed]
- 18. Huang KB, Zhang SP, Zhu YJ, Guo CH, Yang M, Liu J, Xia LG, Zhang JF. Hotair mediates tumorigenesis through recruiting EZH2 in colorectal cancer. J Cell Biochem. 2019; 120:6071–7. https://doi.org/10.1002/jcb.27893 [PubMed]
- 19. Chiyomaru T, Fukuhara S, Saini S, Majid S, Deng G, Shahryari V, Chang I, Tanaka Y, Enokida H, Nakagawa M, Dahiya R, Yamamura S. Long non-coding RNA HOTAIR is targeted and regulated by miR-141 in human cancer cells. J Biol Chem. 2014; 289:12550–65. https://doi.org/10.1074/jbc.M113.488593 [PubMed]
- 20. Nwosu ZC, Ward MH, Sajjakulnukit P, Poudel P, Ragulan C, Kasperek S, Radyk M, Sutton D, Menjivar RE, Andren A, Apiz-Saab JJ, Tolstyka Z, Brown K, et al. Uridine-derived ribose fuels glucose-restricted pancreatic cancer. Nature. 2023; 618:151–8. https://doi.org/10.1038/s41586-023-06073-w [PubMed]
- 21. Rinn JL, Kertesz M, Wang JK, Squazzo SL, Xu X, Brugmann SA, Goodnough LH, Helms JA, Farnham PJ, Segal E, Chang HY. Functional demarcation of active and silent chromatin domains in human HOX loci by noncoding RNAs. Cell. 2007; 129:1311–23. https://doi.org/10.1016/j.cell.2007.05.022 [PubMed]
- 22. Tsai MC, Manor O, Wan Y, Mosammaparast N, Wang JK, Lan F, Shi Y, Segal E, Chang HY. Long noncoding RNA as modular scaffold of histone modification complexes. Science. 2010; 329:689–93. https://doi.org/10.1126/science.1192002 [PubMed]
- 23. Yang L, Peng X, Li Y, Zhang X, Ma Y, Wu C, Fan Q, Wei S, Li H, Liu J. Long non-coding RNA HOTAIR promotes exosome secretion by regulating RAB35 and SNAP23 in hepatocellular carcinoma. Mol Cancer. 2019; 18:78. https://doi.org/10.1186/s12943-019-0990-6 [PubMed]
- 24. Zhao W, Geng D, Li S, Chen Z, Sun M. LncRNA HOTAIR influences cell growth, migration, invasion, and apoptosis via the miR-20a-5p/HMGA2 axis in breast cancer. Cancer Med. 2018; 7:842–55. https://doi.org/10.1002/cam4.1353 [PubMed]. Retraction in: Cancer Med. 2023; 12:22437. https://doi.org/10.1002/cam4.6773 [PubMed]
- 25. Liu XH, Sun M, Nie FQ, Ge YB, Zhang EB, Yin DD, Kong R, Xia R, Lu KH, Li JH, De W, Wang KM, Wang ZX. Lnc RNA HOTAIR functions as a competing endogenous RNA to regulate HER2 expression by sponging miR-331-3p in gastric cancer. Mol Cancer. 2014; 13:92. https://doi.org/10.1186/1476-4598-13-92 [PubMed]
- 26. Hu CY, Su BH, Lee YC, Wang CT, Yang ML, Shen WT, Fu JT, Chen SY, Huang WY, Ou CH, Tsai YS, Kuo FC, Shiau AL, et al. Interruption of the long non-coding RNA HOTAIR signaling axis ameliorates chemotherapy-induced cachexia in bladder cancer. J Biomed Sci. 2022; 29:104. https://doi.org/10.1186/s12929-022-00887-y [PubMed]
- 27. Li P, Zhang X, Wang L, Du L, Yang Y, Liu T, Li C, Wang C. lncRNA HOTAIR Contributes to 5FU Resistance through Suppressing miR-218 and Activating NF-κB/TS Signaling in Colorectal Cancer. Mol Ther Nucleic Acids. 2017; 8:356–69. https://doi.org/10.1016/j.omtn.2017.07.007 [PubMed]
- 28. Xie W, Chu M, Song G, Zuo Z, Han Z, Chen C, Li Y, Wang ZW. Emerging roles of long noncoding RNAs in chemoresistance of pancreatic cancer. Semin Cancer Biol. 2022; 83:303–18. https://doi.org/10.1016/j.semcancer.2020.11.004 [PubMed]
- 29. Wang Q, Li X, Ren S, Su C, Li C, Li W, Yu J, Cheng N, Zhou C. HOTAIR induces EGFR-TKIs resistance in non-small cell lung cancer through epithelial-mesenchymal transition. Lung Cancer. 2020; 147:99–105. https://doi.org/10.1016/j.lungcan.2020.06.037 [PubMed]
- 30. Zhang S, Zheng F, Zhang L, Huang Z, Huang X, Pan Z, Chen S, Xu C, Jiang Y, Gu S, Zhao C, Zhang Q, Shi G. LncRNA HOTAIR-mediated MTHFR methylation inhibits 5-fluorouracil sensitivity in esophageal cancer cells. J Exp Clin Cancer Res. 2020; 39:131. https://doi.org/10.1186/s13046-020-01610-1 [PubMed]
- 31. Siegel RL, Miller KD, Wagle NS, Jemal A. Cancer statistics, 2023. CA Cancer J Clin. 2023; 73:17–48. https://doi.org/10.3322/caac.21763 [PubMed]
- 32. Efimova I, Catanzaro E, Van der Meeren L, Turubanova VD, Hammad H, Mishchenko TA, Vedunova MV, Fimognari C, Bachert C, Coppieters F, Lefever S, Skirtach AG, Krysko O, Krysko DV. Vaccination with early ferroptotic cancer cells induces efficient antitumor immunity. J Immunother Cancer. 2020; 8:e001369. https://doi.org/10.1136/jitc-2020-001369 [PubMed]
- 33. Zhang L, Yu Z, Qu Q, Li X, Lu X, Zhang H. Exosomal lncRNA HOTAIR Promotes the Progression and Angiogenesis of Endometriosis via the miR-761/HDAC1 Axis and Activation of STAT3-Mediated Inflammation. Int J Nanomedicine. 2022; 17:1155–70. https://doi.org/10.2147/IJN.S354314 [PubMed]
- 34. Wang Y, Yi K, Liu X, Tan Y, Jin W, Li Y, Zhou J, Wang H, Kang C. HOTAIR Up-Regulation Activates NF-κB to Induce Immunoescape in Gliomas. Front Immunol. 2021; 12:785463. https://doi.org/10.3389/fimmu.2021.785463 [PubMed]
- 35. Oh DY, Fong L. Cytotoxic CD4+ T cells in cancer: Expanding the immune effector toolbox. Immunity. 2021; 54:2701–11. https://doi.org/10.1016/j.immuni.2021.11.015 [PubMed]
- 36. Xie Z, Xia J, Jiao M, Zhao P, Wang Z, Lin S, Xing Y, Li Y, Lu Z, Zhong Z, Miao C, Zhou P, Qian J, et al. Exosomal lncRNA HOTAIR induces PDL1+ B cells to impede anti-tumor immunity in colorectal cancer. Biochem Biophys Res Commun. 2023; 644:112–21. https://doi.org/10.1016/j.bbrc.2023.01.005 [PubMed]
- 37. Lu L, Zhu G, Zhang C, Deng Q, Katsaros D, Mayne ST, Risch HA, Mu L, Canuto EM, Gregori G, Benedetto C, Yu H. Association of large noncoding RNA HOTAIR expression and its downstream intergenic CpG island methylation with survival in breast cancer. Breast Cancer Res Treat. 2012; 136:875–83. https://doi.org/10.1007/s10549-012-2314-z [PubMed]
- 38. Ward MH, Nwosu ZC, Lyssiotis CA. Uridine: as sweet as sugar for some cells? Cell Res. 2023; 33:898–9. https://doi.org/10.1038/s41422-023-00860-w [PubMed]