Quanfang Chen
Xiaoning Zhong- Department of Respiratory Medicine, First Affiliated Hospital of Guangxi Medical University, Nanning, China
Objectives: Present study aims to identify the essential mRNAs responsible for the development of brain neurovascular-related metastases (BNM) among lung adenocarcinoma (LUAD) patients. Further, we attempted to predict brain metastases more accurately and prevent their development in LUAD patients.
Methods: Transcriptome data analysis was used to identify differentially expressed mRNAs (DEMs) associated with brain metastasis, and thereby the ferroptosis index (FPI) is calculated using a computational model. Meanwhile, the DEmRNAs linked with FPI, and brain metastasis were derived by the intersection of these two groups of DEMs. We also constructed a ceRNA network containing these DEmRNAs, identifying the HCP5 /hsa-miR-17-5p/HOXA7 axis for analysis. Further, a clinical cohort was employed to validate the regulatory roles of molecules involved in the ceRNA regulatory axis.
Results: Here we report the development of a ceRNA network based on BNM-associated DEMs and FPI-associated DEmRNAs which includes three core miRNAs (hsa-miR-338-3p, hsa-miR-429, and hsa-miR-17-5p), three mRNAs (HOXA7, TBX5, and TCF21), and five lncRNAs (HCP5, LINC00460, TP53TG1). Using gene set enrichment analysis (GSEA) and survival analysis, the potential axis of HCP5 /hsa-miR-17-5p/HOXA7 was further investigated. It is found that HOXA7 and ferroptosis index are positively correlated while inhibiting tumor brain metastasis. It may be that HCP5 binds competitively with miR-17-5p and upregulates HOXA7 to increase iron death limiting brain cancer metastases
Conclusions: The expression of both HOXA7 and HCP5 is positively correlated with FPI, indicating a possible link between ferroptosis and BNM. According to the results of our study, the ferroptosis-related ceRNA HCP5 /hsa-miR-17-5p/HOXA7 axis may contribute to the development of BNM in LUAD patients.
Introduction
Lung adenocarcinoma (LUAD) is the most common histological subtype of non-small cell lung cancer (NSCLC) which frequently occurs in peripheral lung tissue while accounting for approximately 70% of NSCLC and 40% of lung tumors (1–4). There has been an increase in LUAD morbidity and mortality rates worldwide in recent years (5–7). LUAD patients have an extremely low average survival rate (8, 9). In some studies, it has been estimated that there are 90,000 LAUD cases and around 8,00,000 LUAD deaths worldwide each year (10, 11). The 5-year survival rate for LUAD is <20% even though treatment options have evolved in recent times (12, 13). Such a low survival rate for LUAD patients may be attributed to its high infiltration and metastasis rates (14). Metastasis of LUAD often occurs in the brain, where the incidence and prognosis are poor, which poses a significant risk to the patients' health (15, 16). This suggests that brain metastases are a significant cause of treatment failure and death in patients with LUAD. Hence, to improve LUAD patients' prognosis, it is critical to identify biomarkers associated with brain metastases (17).
It has been demonstrated that ferroptosis plays a significant regulatory role in the treatment of many diseases, including cardiovascular, kidney, and oncological disorders (18–21). A ferroptosis process has been observed to influence tumor growth in oncological diseases and its role is expected to be exploited as an anti-cancer therapy target (22, 23). Ferroptosis is associated with the development of NSCLC and inducing ferroptosis may improve the therapeutic potential of NSCLC (24, 25). Some studies suggest that ferroptosis in lung cancer is regulated by USP35, which is thought to affect the growth and progression of the disease and may represent a new therapeutic target (22). Furthermore, ferroptosis appears to enhance the antitumor effects of conventional chemotherapy and radiotherapy (26–28). However, an exact mechanism of action for ferroptosis in LUAD brain metastases is still unknown.
Additionally, the ceRNA network is involved in regulating multiple tumor types (29–32). Several studies have shown that the ceRNA network has an association with lung cancer prognosis and NSCLC cell proliferation (33, 34). There was previously a finding that lincRNA00494 as ceRNA suppresses the proliferation of non-small cell lung cancer cells by regulating SRCIN1 expression (33). Studies have also shown that ceRNA regulatory networks can be used to construct models for assessing the prognosis of patients with malignancies (35). Thus, ferroptosis-related ceRNA may be a good candidate for the identification of new biomarkers in lung cancer diagnosis and prognosis (36).
The major focus of our study was to identify the major mRNAs involved in the development of BNM among LUAD patients and to explore the potential of the ceRNA regulation network for identifying new biomarkers. In addition, this research also intended to uncover any possible therapeutic target to monitor and improve the risk of BNM in LUAD patients with the possibility to extend their survival period.
Methods
Data Download
Due to the lack of BNM information in The Cancer Genome Atlas (TCGA, accessible at https://portal.gdc.cancer.gov/) project, we downloaded the transcriptome profile of LUAD brain neurovascular-related metastases (BNMs) from GSE141685 in the Gene Expression Omnibus (GEO) database (37). This BNM-LUAD dataset is containing bulky RNA sequencing (RNA-seq) data from surgically resected brain metastases of 6 LUAD patients. To investigate the potential molecules and pathways involved in LUAD brain metastasis, we have compared the BNM-LUAD bulk RNA-seq data with nM-LUAD (LUAD with no metastasis) bulk RNA-seq data of 384 patients from the TCGA database. From the TCGA database, we also downloaded miRNA sequences and clinical characterization data on nM-LUAD. In addition, the raw transcript RNA-seq data from the BNM-LUAD dataset were stored in Sequence Read Archive (SRA) file format and downloaded locally for subsequent RNA-seq analysis.
Data Preparation and Recruitment of Patients
The Kallisto software was used to quantify transcript abundance on bulk RNA-seq data (38). The transcripts per million (TPM) units of gene expression are included in the quantitative results of the transcriptome data and the log (TPM+1) for gene expression is calculated. In the next step, we applied the “ComBat” function from the limma package to remove the batch effect and pooling of data (39). To validate the findings of this study, an independent cohort of LUAD patients was also utilized. The Institutional Ethics Review Board of The First Affiliated Hospital of Guangxi Medical University has approved this study and informed written consent was also obtained from all the patients. We randomly collected and then examined 29 tumor tissue samples from the First Affiliated Hospital of Guangxi Medical University where this cohort was enrolled as a clinical cohort. A combination of diagnostic imaging and histology is used to better understand the clinical diagnosis. To obtain tumor tissue samples, surgical excision was used, and excess tumor tissues were collected only after a confirmed pathological diagnosis.
Identification of FPI-Related Genes
In the TCGA nM-LUAD dataset, we used the GSVA package to calculate the enrichment scores of ferroptosis-related gene sets (including ACSL4, ALOX15, COQ10A, COQ10B, FDFT1, GPX4, HMGCR, SLC3A2, SLC7A11, NFE2L2, NOX1, NOX3, NOX4, and NOX5). This followed the ferroptosis index (FPI) calculation by applying the standardized method as described previously (40, 41). These FPI scores are used to quantify the likelihood of ferroptosis of cells within a tumor or tissue and based on median FPI scores, nM-LUAD patients are divided into high and low FPI groups. We next performed differential expression analyses for lncRNAs, miRNAs, and mRNAs between the low and high FPI groups with the limma package (p < 0.05). These differentially expressed mRNAs (DEMs) were defined as ferroptosis-related or FPI-related mRNAs (including FPI-related DEmRNAs, DEmiRNAs, and DElncRNAs). As described previously by utilizing the clusterProfiler R package the FPI-associated DEmRNAs are used for functional pathway analysis including the Kyoto Encyclopedia of Genes and Genomes (KEGG) and gene ontology (GO) (42–44).
Identification of Brian Metastasis-Related Genes
In this study, a significance level of p < 0.05 was set as the cut-off criterion. We have also performed a differential expression analysis between the BNM and no metastasis groups using the “limma” package. The calculated DEMs were named BNM-related mRNAs which are associated with the development of BNM in LUAD.
Identification of Genes Potentially Regulated by Identified FPI-Related miRNAs
We performed intersection analysis studies while including FPI-related DEmRNAs and BNM-related DEmRNAs that may play a crucial role in ferroptosis and brain metastases. Firstly, the miRcode database confirmed a relationship between delncRNAs and demiRNAs (45, 46). After that, we used miRWalk (http://mirwalk.umm.uni-heidelberg.de/) to identify shared DEmiRNAs and DEmRNAs between FPI and BNM (47). As a result, ceRNA networks were constructed based on DElncRNA-DEmiRNA interaction pairs and DEmiRNA-DEmiRNA interaction pairs.
Identification of Genes as the Target of FPI-Related miRNAs Associates to FPI, Brain Metastasis, and FPI-Related miRNAs
The correlation analysis was used to screen possible lncRNA-miRNA and miRNA-mRNA interaction pairs further. A Cox regression analysis was conducted on the molecules in these interaction pairs to predict the probability of survival of LUAD patients after 1, 3, and 6 years. To quantify and predict the survival rates of LUAD patients, the time-dependent receiver operating characteristic curve (ROC) and area under the curve (AUC) were used. We used the “survival” R package to analyze survival data and the “survminer” R package to visualize the data. A Pearson's method and linear regression analysis were then performed to determine the correlation between the variables. The “ggplot2” package was used to plot the heat map and scatter plots of the statistical tests.
Gene Set Enrichment Analysis
GSEA (Gene Set Enrichment Analysis) is a bioinformatics technique for exploring a specific set of functional genes (48, 49). Two groups are generated depending on the expression level of the target gene and critical pathway analysis was then conducted using GSEA. We used the reference dataset “c2.cp.v7.2.symbols.gmt [Curated]” for the enrichment analysis. All adjusted p-values < 0.05 were considered significantly enriched.
Quantitative Reverse Transcription PCR (qRT-PCR)
To quantify the abundance of mRNAs and miRNAs in samples, the total RNA was extracted from cells using Trizol reagent (Invitrogen) and quantified using Nanodrop (Thermo Scientific, Waltham, USA). The expression of miRNAs was determined by stem-loop qPCR (Taq-Man) as mentioned previously as a reference. We conducted a qPCR on the cDNA template with TB GreenTM Premix Ex TaqTM II (Takara; RR820A) to determine the levels of mRNA. In addition, qPCR primer sets for miRNAs as well as mRNAs were obtained from RiboBio while mRNA qPCR primers have been synthesized by Sangon (Shanghai, China), as described in previous studies (50). For the relative quantification of gene expression levels, the CT method was performed in triplicates while GAPDH or U6 are included as internal controls.
Statistical Methods
Statistical analyses were performed using R software (version 4.0.2). The Kaplan-Meier method was used, and survival curves were plotted for survival analyses. Unless otherwise stated, p < 0.05 was found to be significant.
Results
DEGs Related to Ferroptosis Index in LUAD
The flow chart of this study is shown in Figure 1. We calculated the probability of ferroptosis based on the core mRNAs associated with ferroptosis and it was termed the ferroptosis index (FPI) (40). According to their median FPI values, 384 LUAD patients without metastasis (nM-LUAD) in the TCGA-LUAD project were classified into high or low FPI groups. We then identified 5,958 DEmRNAs, 2,380 lncRNAs, and 68 DEmiRNAs to analyze the differences between the low and high FPI groups (Figures 2A–C). It has been proposed that these dysregulated mRNAs are associated with ferroptosis and are known as FPI-related mRNAs. The DEmRNAs that were upregulated in the low FPI group were enriched in RNA splicing, stereocilium, presynaptic active zone, and RNA polymerase complex (Figure 2D). Conversely, pathways that were downregulated in the low FPI group included collagen–containing extracellular matrix, positive regulators of cell adhesion, focal adhesion, and cell adhesion molecule binding (Figure 2E). The analysis described above can identify the mRNAs associated with FPI in LUAD, as well as their possible functions.
Figure 1. The research strategy for this study.
Figure 2. FPI-related DEGs and functional pathway analysis. Volcano plots showed FPI-related (A) DEmRNAs, (B) DEmiRNAs, and (C) DElncRNAs. The plots show the predicted RNAs in ceRNA networks. Red dots indicate upregulated genes in low FPI groups. The low FPI fraction group (D) upregulated (E) downregulated GO and KEGG pathways.
Identification of Shared DEmRNAs and Construction of ceRNA Network
To identify DEmRNAs shared with both FPI and BNM, the differential analysis between BNM-LUAD and nM-LUAD revealed 154 BNM-related DEmRNAs. Then, the intersection analysis of BNM-related DEmRNAs and FPI-related DEmRNAs identified 46 downregulated 44 mRNAs that are common to both brain metastases and ferroptosis (Figures 3A,B). BNM-LUAD downregulated 44 out of the 46 shared mRNAs. Downregulated mRNAs were mostly found to be enriched in leukocyte-mediated cytotoxicity MHC protein complexes, segment specification, HMG box domain binding, allograft rejection, and graft-vs.-host diseases (Figure 3C). Immunologically-related pathways may likely be downregulated in patients with LUAD, contributing to brain metastases.
Figure 3. An intersection analysis of FPI- and BM-associated DEmRNAs, and the construction of ceRNA networks. (A) Venn plot showed 46 shared DEmRNAs among FPI- and BNM-associated DEmRNAs. In BNM and the low FPI group, 44 of these shared DEmRNAs were downregulated. (B) Expression of genes shared by brain metastases and non-metastases. (C) The functional pathways of the shared DEmRNAs downregulated in BNM. (D) Construction of ceRNA networks based on the shared DEmRNAs.
Based on these 44 DEmRNAs [RT1], we constructed a ceRNA network (Figure 3D). The DElncRNAs and DEmiRNAs in the network were derived from DEMs related to FPI. The ceRNA network consists of three core miRNAs (hsa-miR-338-3p, hsa-miR-429, and hsa-miR-17-5p), three mRNAs (HOXA7, TBX5, and TCF21), and five lncRNAs (HCP5, LINC00460, TP53TG1, RGS5, and MAGI2-AS3).
Identification of Potential HCP5/hsa-miR-17-5p/HOXA7 Pathways From the ceRNA Network
Based on the constructed ceRNA network, we performed correlation analysis to identify potential pathways. An analysis of the correlation between these 11 dysregulated molecules is depicted in Figure 4A. Significant correlations (p < 0.01) were observed between the mRNA and miRNA groups. The heat map in Figure 4B shows the expression of these 11 dysregulated molecules in LUAD. The correlation studies identified the HCP5 /hsa-miR-17-5p/HOXA7 axes associated with ferroptosis and BNM from the ceRNA network (Figures 4C,D). There was a negative correlation between the expression of HCP5 and hsa-miR-17-5p (Pearson r = −0.13, p 0.05). Also, the expression of hsa-miR-17-5p and HOXA7 exhibited a negative correlation (Pearson r = −0.17, p 0.01). Additionally, we found a positive correlation between FPI and HXOA7 and HCP5 (p < 0.01), as well as a negative correlation between FPI and hsa-miR-17-5p (p < 0.05). HCP5 overexpression was also linked to a superior primary therapy outcome, suggesting that it may improve the prognosis of LUAD patients (Table 1). Thus, ceRNA HCP5 /hsa-miR-17-5p/HOXA7 may play a role in the development of LUAD.
Figure 4. Identification of potential HCP5 /hsa-miR-17-5p/HOXA7 pathways from the ceRNA network. (A) Correlation analysis of each element in the ceRNA network. (B) FPI and the expression of individual elements in the HCP5 /hsa-miR-17-5p/HOXA7 pathway have statistically significant connections. (C,D) Ferroptosis is significantly linked to each molecule in the HCP5 /hsa-miR-17-5p/HOXA7 pathway; (C) shows the correlation coefficients and (D) shows the statistical differences.
Table 1. The relation between HCP5 expression levels and LUAD clinicopathological features.
Correlation of HCP5, hsa-miR-17-5p, and HOXA7 With LUAD Prognosis and FPI
To further evaluate the effect of HCP5 /hsa-miR-17-5p/HOXA7 in LUAD, ROC curves and analysis of differential expression were used to examine the relationship between HCP5, hsa-miR-17-5p, and HOXA7 and LUAD prognosis and FPI. In ROC analysis, HCP5 (AUC = 0.596, 95% CI: 0.534-0.657), HOXA7 (AUC = 0.617, 95% CI: 0.556–0.678) and hsa-miR-17-5p (AUC = 0.575, 95% CI: 0.512–0.637) significantly predicted low FPI scores (Figures 5A–C). In addition, time-dependent ROC curves indicated that HCP5, HSA-miR-17-5p, and HOXA7 might be associated with the survival of LUAD patients (Figures 5D,E). In addition, there was a statistically significant difference in the expression levels of HCP5, HOXA7, and hsa-miR-17-5p between the high and low FPI groups (Figure 5G). The potential HCP5 /hsa-miR-17-5p/HOXA7 axis is therefore thought to be related to ferroptosis.
Figure 5. Correlation between HCP5, hsa-miR-17-5p, and HOXA7 with LUAD prognosis and FPI. (A–C) ROC curves predicting high or low FPI scores for (A) HCP5, (B) HOXA7, (C) hsa-miR-17-5p. (D–F) ROC curves for predicting 1-, 3-, and 6-year survival rates. (G) Comparison of the expression of HCP, hsa-miR-17-5p, and HOXA7 in high and low FPI groups. *P < 0.05, **P < 0.01, and ***P < 0.001.
Clinical Correlation Analysis of HCP5, hsa-miR-17-5p, and HOXA7
To explain the potential effectiveness of the HCP5 /hsa-miR-17-5p/HOXA7 axis in LUAD at the clinical level, we performed qRT-PCR using available tumor samples from 29 LUAD patients. The expression of HCP5 and HOXA7 correlated positively in clinical samples (Figure 6A; Pearson r = 0.45, p = 0.014), while hsa-miR-17-5p correlated negatively with HCP5 (Figure 6B; Pearson r = −0.39, p = 0.037) as well as HOXA7 (Figure 6C; Pearson r = −0.548, p = 0.002). Furthermore, HOXA7 expression and hsa-miR-17-5p expression were negatively related. In this study, lncRNA HCP5 was found to act as an “RNA sponge” to sequester hsa-miR-17-5p, thereby decreasing the effects of hsa-miR-17-5p on the target HOXA7 mRNA.
Figure 6. Validation of the HCP5 /hsa-miR-17-5p/HOXA7 pathway using GSEA and correlation analysis. (A) A positive correlation between the expression of HOXA7 and HCP5 in clinical cohort samples (Pearson r = 0.45, p = 0.014). (B) The expression of hsa-miR-17-5p and HCP5 was negatively correlated (r = −0.39, p = 0.037). (C) The expression of HOXA7 and hsa-miR-17-5p was negatively correlated (Pearson r = −0.548, p = 0.002). (D,E) GSEA demonstrates pathways involving HOXA7 and HCP5.
Functional Analysis of HCP5 and HOXA7
Using GSEA, we determined the relationship between pathway activation and inhibition when HCP5 and HOXA7 expression levels were high and low, respectively. The GSEA indicated that the tumor samples with elevated HOXA7 expression downregulated miRNA targets in ECM and membrane receptors, as well as metastatic brain tumor pathways, while the apoptosis pathway was upregulated (Figure 6D). Additionally, the tumor samples with elevated HCP5 expression downregulated miRNA targets in ECM and membrane receptors, metastatic brain tumors, as well as and miRNA regulation of the p53 pathway in prostate cancer (Figure 6E). Therefore, increased levels of HOXA7 and HCP5 were associated with the downregulation of tumor brain metastasis. Additionally, the expression of HOXA7 and HCP5 was positively correlated with FPI, indicating a strong relationship between ferroptosis and BNM.
Discussion
According to this study, BNM and ferroptosis are linked in LUAD. It has been observed that the ferroptosis-related ceRNA HCP5/hsa-miR-17-5p/HOXA7 axis may play an important role in the development of BNM in LUAD.
We examined the BNM-related DEMs in LUAD with and without metastatic tumors. In addition, we investigated DEmRNAs associated with FPI in LUAD patients without metastatic tumors between high and low FPI. Furthermore, the intersection studies identified DEmRNAs shared between BMogenesis and ferroptosis. Based on intersection analysis, a ceRNA network associated with BNM and ferroptosis was constructed. The HCP5/hsa-miR-17-5p/HOXA7 axis is thought to play a central role in the ceRNA regulatory network of LUAD. In a clinical cohort, we validated the HCP5 /hsa-miR-17-5p/HOXA7 axis and confirmed our bioinformatic findings.
Previous studies suggest that HOXA7 participates in the pathogenesis of LUAD, providing a new target for the early treatment of LUAD (51). In brain tumors, HOXA7 is involved in glioma progression and affects patient prognosis (52). However, there is a lack of research on the role of HOXA7 in the brain metastasis of lung cancer. Several studies have demonstrated that HOXA7 is involved in tumor metastasis, suggesting that HOXA7 expression is different in different types of cancer.
HOXA1 has diverse effects on tumor progression depending on the type of malignancy. By downregulating HOXA7 and promoting cell migration and invasion, miR-196a plays a pro-cancer role in colorectal cancer (53). It has been shown that HOXA7 expression is elevated in metastatic hepatocellular carcinoma, enhancing the metastasis of hepatocellular carcinoma, and is associated with a poor prognosis for patients (54). Also, HOXA7 is upregulated in the ceRNA network in oral squamous cell carcinoma, increasing the invasion and migration of oral cancer cells (55). Our study is the first to investigate the role of HOXA7 in brain metastasis in LUAD. In patients with brain metastases, HOXA7 expression has been observed to be downregulated. Downregulation of HOXA7 can lead to tumors becoming more likely to metastasize to the brain, which may be associated with a poor prognosis in patients who have brain metastases. Our findings showed a positive correlation between HOXA7 and the ferroptosis index, suggesting HOXA7 may promote ferroptosis. Some studies have revealed that negative regulation of ferroptosis can lead to increased resistance to chemotherapy in lung cancer patients with brain metastases (56). In addition, the promotion of ferroptosis inhibited brain metastasis in HER2-positive breast cancer (57). Thus, we hypothesized that HOXA7 downregulation promotes LUAD brain metastasis, which may be related to ferroptosis induced by HOXA7.
Li et al. demonstrated that miR-17-5p is highly expressed in LUAD and is associated with patient prognosis as a novel marker for clinical diagnosis of NSCLC (58, 59). Moreover, miR-17-5p is a direct target of metastasis-associated LUAD transcript 1 (MALAT1) (60). Additionally, miR-17-5p has been shown to directly target MALT1 (metastasis-associated LUAD transcript 1). A study has shown that increased MALAT1 expression promotes brain metastasis in lung cancer and correlates with patient survival (61). Additionally, miR-17-5p is also involved in brain metastasis of other cancers (62). It is therefore thought that miR-17-5p promotes LUAD brain metastasis, and this may be due to the inhibition of ferroptosis. Even though no studies have indicated a link between miR-17-5p and ferroptosis, some studies have found that miR-17-5p can resist lung cancer cell apoptosis (63). In the present study, miR-17-5p has been observed to negatively regulate HOXA7, thus possibly promoting LUAD brain metastasis. According to our correlation analysis, miR-17-5p promotes LUAD brain metastasis by inhibiting ferroptosis and negatively targeting HOXA7.
HCP5 affects tumor progression differently depending on the type of malignancies. Several studies have explored the relationship between HCP5 and metastasis and prognosis in cancer (64). It was reported by Jiang et al. that HCP5 was highly expressed in tumor tissues of LUAD patients and contributed to epithelial-mesenchymal transition (EMT) of LUAD cells, tumor growth, and metastasis, and was positively correlated with poor patient prognosis (65). Brain metastasis associated with lung cancer is also attributed to EMT (61, 66). We observed that patients with high HCP5 expression had a lower prognosis than those with low HCP5. We found, however, that HCP5 may regulate HOXA7 by binding to miR-17-5p and inhibiting LUAD brain metastasis by adsorbing miR-17-5p, indicating a different mechanism of action for HCP5 to affect LUAD brain metastasis. Furthermore, our study revealed a positive correlation between HCP5 and ferroptosis index, suggesting that HCP5 may promote ferroptosis. Cancer metastasis is promoted by HCP5 through EMT in several tumor diseases (65, 67). Ferroptosis is enhanced by the induction of EMT in cancer cells (68). We, therefore, proposed that HCP5 promotes ferroptosis and inhibits LUAD brain metastasis by upregulating HOXA7 through adsorption of miR-17-5p.
We can speculate that dysregulation of HCP5 expression in LUAD patients may lead to tumor brain metastasis by downregulating HOXA7 by affecting miR-17-5p. Therefore, a decrease in HOXA7 expression may promote the development of brain metastases. In the present study, HOXA7 levels were significantly lower in the BNM group than in the metastasis-free group. According to GSEA, elevated HOXA7 and HCP5 were linked to a decreased incidence of brain metastasis in tumors and may affect brain metastasis by regulating cell membrane surface receptors. In addition, the expression of both HOXA7 and HCP5 was positively correlated with FPI, suggesting a strong correlation between ferroptosis and tumor brain metastasis. The study demands further investigation to explore the relationship between brain tumors and ferroptosis.
The study was limited by the small number of subjects, including clinical study, and the lack of in vitro cell-based validation. We plan to further validate the ceRNA HCP5 /hsa-miR-17-5p/HOXA7 axis' effect on LUAD BNM from basic experiments and clinical samples in our future studies.
Conclusion
In conclusion, bioinformatics and clinical studies have demonstrated a correlation between the expression of HOXA7 and HCP5. It may be possible to link ferroptosis to tumor neurovascular brain metastasis from correlation studies. Moreover, our study demonstrated that the ferroptosis-related ceRNA HCP5 /hsa-miR-17-5p/HOXA7 axis may play an essential role in LUAD BNM.
Data Availability Statement
Publicly available datasets were analyzed in this study. This data can be found here: National Institute of Health (NIH), National Cancer Institute (NCI), The Cancer Genome Atlas (TCGA) Program, https://portal.gdc.cancer.gov/ and TCGA-LUAD and National Center for Biotechnology Information (NCBI), Gene Expression Omnibus (GEO), https://www.ncbi.nlm.nih.gov/geo/, GSE141685.
Ethics Statement
Ethical review and approval was not required for the study on human participants in accordance with the local legislation and institutional requirements. Written informed consent from the patients/participants or patients/participants' legal guardian/next of kin was not required to participate in this study in accordance with the national legislation and the institutional requirements.
Author Contributions
QC, QP, and HG performed the data collection and analyses. QC and XZ analyzed and interpreted the results. QC, YW, and XZ drafted and reviewed the manuscript. All authors read and approved the final manuscript.
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.
Acknowledgments
We thank the reviewers for their outstanding contribution to the improvement of this manuscript. We thank Bullet Edits Limited for the linguistic editing and proofreading of the manuscript.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fneur.2022.878947/full#supplementary-material
Abbreviations
TCGA, the cancer genome atlas; GEO, gene expression omnibus (GEO); LUAD, lung adenocarcinoma; BM, brain metastasis; nM-LUAD, LUAD without metastasis; ceRNA, competing endogenous RNAs; DE-mRNAs, differentially expressed mRNAs; DE-miRNAs, differentially expressed microRNAs; DE-lncRNAs, differentially expressed lncRNAs; FPI, ferroptosis Potential Index; GO, gene ontology; GSEA, gene set enrichment analysis; KEGG, kyoto encyclopedia of genes and genomes; K-M, kruskal-wallis; OS, overall survival; qRT-PCR, quantitative reverse transcription PCR.
References
1. Motono N, Funasaki A, Sekimura A, Usuda K, Uramoto H, et al. Prognostic value of epidermal growth factor receptor mutations and histologic subtypes with lung adenocarcinoma. Med Oncol. (2018) 35:22. doi: 10.1007/s12032-018-1082-y
2. Dong Y, Yang W, Wang J, Zhao J, Qiang Y, Zhao Z, et al. MLW-gcForest: a multi-weighted gcForest model towards the staging of lung adenocarcinoma based on multi-modal genetic data. BMC Bioinformatics. (2019) 20:578. doi: 10.1186/s12859-019-3172-z
3. Zhang W Shen Y Feng Feng G: Predicting the survival of patients with lung adenocarcinoma using a four-gene prognosis risk model. Oncol Lett. (2019). doi: 10.3892/ol.2019.10366
4. Du J, Zhang G, Qiu H, Yu H, Yuan W. The novel circular RNA circ-CAMK2A enhances lung adenocarcinoma metastasis by regulating the miR-615-5p/fibronectin 1 pathway. Cell Mol Biol Lett. (2019) 24:72. doi: 10.1186/s11658-019-0198-1
5. Jemal A, Siegel R, Ward E, Hao Y, Xu J, Thun MJ. Cancer statistics, 2017. CA: A Cancer J Clin. (2017) 67:7–30. doi: 10.3322/caac.21387
6. Jemal, Bray F, Center MM, Ferlay J, Ward E, Forman. Global cancer statistics. CA: A Cancer J Clin. (2011) 61:69–90. doi: 10.3322/caac.20107
7. Sun F, Huang Y, Yang X, Zhan C, Xi J, Lin Z, et al. Solid component ratio influences prognosis of GGO-featured IA stage invasive lung adenocarcinoma. Cancer Imaging. (2020) 20:87. doi: 10.1186/s40644-020-00363-6
8. Hirsch FR, Scagliotti GV, Mulshine JL, Kwon R, Curran WJ Jr, Wu YL, et al. Lung cancer: current therapies and new targeted treatments. The Lancet. (2017) 389:299–311. doi: 10.1016/S0140-6736(16)30958-8
9. Yu Z, Zhu X, Li Y, Liang M, Liu M, Liu Z, et al. Circ-HMGA2 (hsa_circ_0027446) promotes the metastasis and epithelial-mesenchymal transition of lung adenocarcinoma cells through the miR-1236-3p/ZEB1 axis. Cell Death Dis. (2021) 12:313. doi: 10.1038/s41419-021-03601-2
10. Ferlay J, Soerjomataram I, Dikshit R, Eser S, Mathers C, Rebelo M, et al. Cancer incidence and mortality worldwide: sources, methods and major patterns in Globocan 2012: Globocan (2012). Int. J. Cancer. (2015) 136:E359–E386. doi: 10.1002/ijc.29210
11. Zheng H, Huang Q, Huang S, Yang X, Zhu T, Wang W, et al. Senescence inducer Shikonin ROS-dependently suppressed lung cancer progression. Front. Pharmacol. (2018) 9:519. doi: 10.3389/fphar.2018.00519
12. Chen Z, Fillmore CM, Hammerman PS, Kim CF, Wong KK. Non-small-cell lung cancers: a heterogeneous set of diseases. Nat Rev Cancer. (2014) 14:535–46. doi: 10.1038/nrc3775
13. Walter JE, Heuvelmans MA, de Jong PA, Vliegenthart R, van Ooijen PM, Peters RB, et al. Occurrence and lung cancer probability of new solid nodules at incidence screening with low-dose CT: analysis of data from the randomised, controlled NELSON trial. Lancet Oncol. (2016) 17:907–16. doi: 10.1016/S1470-2045(16)30069-9
14. Walter JE, Heuvelmans MA, de Jong PA, Vliegenthart R, van Ooijen PM, Peters RB, et al. Comprehensive profiling reveals distinct microenvironment and metabolism characterization of lung adenocarcinoma. Front. Genet. (2021) 12:619821. doi: 10.3389/fgene.2021.619821
15. Ning M, Chunhua M, Rong J, Yuan L, Jinduo L, Bin W, et al. Diagnostic value of circulating tumor cells in cerebrospinal fluid. Open Medicine. (2016) 11:21–4. doi: 10.1515/med-2016-0005
16. Ning M, Chunhua M, Rong J, Yuan L, Jinduo L, Bin W, et al. Clinical efficacy and safety evaluation of pemetrexed combined with radiotherapy in treatment of patients with lung adenocarcinoma brain metastasis. Oncol Lett. (2019). doi: 10.3892/ol.2019.9894
17. Kawabe T, Phi JH, Yamamoto M, Kim DG, Barfod BE, Urakawa Y. Treatment of brain metastasis from lung cancer. In: Kim, D.G. and Lunsford, L.D. (eds.) Progress in Neurol Surg. (2012):148–155. KARGER, Basel. doi: 10.1159/000331188
18. Fang X, Wang H, Han D, Xie E, Yang X, Wei J, et al. Ferroptosis as a target for protection against cardiomyopathy. Proc Natl Acad Sci USA. (2019) 116:2672–80. doi: 10.1073/pnas.1821022116
19. Friedmann Angeli JP, Schneider M, Proneth B, Tyurina YY, Tyurin VA, Hammond VJ, et al. Inactivation of the ferroptosis regulator Gpx4 triggers acute renal failure in mice. Nat Cell Biol. (2014) 16:1180–91. doi: 10.1038/ncb3064
20. Liang C, Zhang X, Yang M, Dong X. Recent progress in ferroptosis inducers for cancer therapy. Adv Mater. (2019) 31:1904197. doi: 10.1002/adma.201904197
21. Qiu C, Zhang X, Huang B, Wang S, Zhou W, Li C, et al. Disulfiram, a ferroptosis inducer, triggers lysosomal membrane permeabilization by up-regulating ROS in glioblastoma. OTT. (2020) 13:10631–40. doi: 10.2147/OTT.S272312
22. Tang Z, Jiang W, Mao M, Zhao J, Chen J, Cheng N. Deubiquitinase USP35 modulates ferroptosis in lung cancer via targeting ferroportin. Clin Translat Med. (2021) 11:4. doi: 10.1002/ctm2.390
23. Hassannia B, Vandenabeele P, Berghe TV. Targeting ferroptosis to iron out cancer. Cancer Cell. (2019) 35:830–49. doi: 10.1016/j.ccell.2019.04.002
24. Li G, Yang J, Zhao G, Shen Z, Yang K, Tian L, et al. Dysregulation of ferroptosis may involve in the development of non-small-cell lung cancer in Xuanwei area. J Cell Mol Med. (2021) 25:2872–84. doi: 10.1111/jcmm.16318
25. Tang X, Ding H, Liang M, Chen X, Yan Y, Wan N, et al. Curcumin induces ferroptosis in non-small-cell lung cancer via activating autophagy. Thorac Cancer. (2021) 12:1219–1230. doi: 10.1111/1759-7714.13904
26. Lu B, Chen XB, Ying MD, He QJ, Cao J, Yang B. The role of ferroptosis in cancer development and treatment response. Front Pharmacol. (2018) 8:992.. doi: 10.3389/fphar.2017.00992
27. Lei G, Zhang Y, Koppula P, Liu X, Zhang J, Lin SH, et al. The role of ferroptosis in ionizing radiation-induced cell death and tumor suppression. Cell Res. (2020) 30:146–162. doi: 10.1038/s41422-019-0263-3
28. Ye LF, Chaudhary KR, Zandkarimi F, Harken AD, Kinslow CJ, Upadhyayula PS, et al. Higgins, Radiation-induced lipid peroxidation triggers ferroptosis and synergizes with ferroptosis inducers. ACS Chem Biol. 15:469–84. (2020). doi: 10.1021/acschembio.9b00939
29. Hu H, Liu S, Chu A, Chen J, Xing C, Jing J. Comprehensive analysis of ceRNA network of ERCC4 in colorectal cancer. PeerJ. (2021) 9:e12647. doi: 10.7717/peerj.12647
30. Jia X, Meng W, Zhang L, Jia Y, Shi Y, Tong Z. Construction of differentially expressed Her-2 related lncRNA-mRNA-miRNA ceRNA network in Her-2 positive breast cancer. Transl Cancer Res TCR. (2020) 9:2527–33. doi: 10.21037/tcr.2020.03.34
31. Wang Y, Zhang F, Wu D, Wang Q, Nie L, Yu J. A novel circ_0099999/miR-330-5p/FSCN1 ceRNA crosstalk in pancreatic cancer. Autoimmunity. (2021) 54:471–82. doi: 10.1080/08916934.2021.1963958
32. Chu A, Liu J, Yuan Y, Gong Y. Comprehensive analysis of aberrantly expressed ceRNA network in gastric cancer with and without H.pylori infection. J. Cancer. (2019) 10:853–63. doi: 10.7150/jca.27803
33. Dong J, Li B, Lin D, Lu D, Liu C, Lu X, et al. LincRNA00494 suppresses non-small cell lung cancer cell proliferation by regulating SRCIN1 expression as a ceRNA. Front Oncol. 10:79 (2020). doi: 10.3389/fonc.2020.00079
34. Cao Q, Dong Z, Liu S, An G, Yan B, Lei L. Construction of a metastasis-associated ceRNA network reveals a prognostic signature in lung cancer. Cancer Cell Int. (2020) 20:208. doi: 10.1186/s12935-020-01295-8
35. Xu J, Ling T, Dai S, Han S, Ding K. Constructing the ceRNA regulatory network and combining immune cells to evaluate prognosis of colon cancer patients. Front Cell Dev Biol. (2021) 9:686844. doi: 10.3389/fcell.2021.686844
36. Zhao Y, Ma S, Cui Z, Li S, Chen Y, Yin Y, et al. The relationship between LncRNAs and lung adenocarcinoma as well as their ceRNA network. (2021) CBM. 31:165–176. doi: 10.3233/CBM-203078
37. Shih DJ, Nayyar N, Bihun I, Dagogo-Jack I, Gill CM, Aquilanti E, et al. Genomic characterization of human brain metastases identifies drivers of metastatic lung adenocarcinoma. Nat Genet. (2020) 52:371–377. doi: 10.1038/s41588-020-0592-7
38. Bray NL, Pimentel H, Melsted P, Pachter L. Near-optimal probabilistic RNA-seq quantification. Nat Biotechnol. (2016) 34:525–527. doi: 10.1038/nbt.3519
39. Ritchie ME, Phipson B, Wu DI, Hu Y, Law CW, Shi W, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. (2015) 43:e47. doi: 10.1093/nar/gkv007
40. Liu Z, Zhao Q, Zuo ZX, Yuan SQ, Yu K, Zhang Q, et al. Systematic analysis of the aberrances and functional implications of ferroptosis in cancer. iScience. (2020) 23:101302. doi: 10.1016/j.isci.2020.101302
41. Chen C, Zhao J, Liu JN, Sun C. Mechanism and role of the neuropeptide LGI1 Receptor ADAM23 in regulating biomarkers of ferroptosis and progression of esophageal cancer. Dis Markers. (2021) 2021:9227897. doi: 10.1155/2021/9227897
42. Yu G, Wang LG, Han Y, He QY. Clusterprofiler: an R package for comparing biological themes among gene clusters. OMICS: A J Integrative Biol. (2012) 16:284–287. doi: 10.1089/omi.2011.0118
43. Sun C, Ma S, Chen Y, Kim NH, Kailas S, Wang Y, et al. Diagnostic value, prognostic value, and immune infiltration of LOX family members in liver cancer: bioinformatic analysis. Front. Oncol. (2022) 12:843880. doi: 10.3389/fonc.2022.843880
44. Lu L, Kang X, Yi B, Jiang C, Yan X, Chen B, et al. Exploring the mechanism of Yiqi Qingre Ziyin method in regulating neuropeptide expression for the treatment of atrophic rhinitis. Disease Markers. (2022) 2022:1–12. doi: 10.1155/2022/4416637
45. Jeggari A, Marks DS, Larsson E. miRcode: a map of putative microRNA target sites in the long non-coding transcriptome. Bioinformatics. (2012) 28:2062–2063. doi: 10.1093/bioinformatics/bts344
46. Huang J, Liang X, Cai. A potential cerna network for neurological damage in preterm infants. BioMed Research International. (2021) 2021:2628824. doi: 10.1155/2021/2628824
47. Sticht C, De La Torre C, Parveen A, Gretz N. miRWalk: An online resource for prediction of microRNA binding sites. PLoS ONE. (2011) 3:e0206239. doi: 10.1371/journal.pone.0206239
48. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA et al. Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Procee Nat Academy Sci. (2005) 102:15545–50. doi: 10.1073/pnas.0506580102
49. Chen Y, Sun Y, Xu Y, Lin WW, Luo Z, Han Z, et al. Single-cell integration analysis of heterotopic ossification and fibrocartilage developmental lineage: endoplasmic reticulum stress effector xbp1 transcriptionally regulates the notch signaling pathway to mediate fibrocartilage Differentiation.”Oxidative Medicine and Cellular Longevity. (2021) 2021:1–29. doi: 10.1155/2021/7663366
50. Luo Z, Sun Y, Qi B, Lin J, Chen Y, Xu Y, et al. Human bone marrow mesenchymal stem cell-derived extracellular vesicles inhibit shoulder stiffness via let-7a/Tgfbr1 axis. Bioact Mat. (2022) 17:S2452199X22000226. doi: 10.1016/j.bioactmat.2022.01.016
51. Kang JU. Characterization of amplification patterns and target genes on the short arm of chromosome 7 in early-stage lung adenocarcinoma. Mol Med Reports. (2013) 8:1373–1378. doi: 10.3892/mmr.2013.1686
52. Li B, McCrudden CM, Yuen HF, Xi X, Lyu P, Chan KW, et al. CD133 in brain tumor: the prognostic factor. Oncotarget. (2017) 8:11144–59. doi: 10.18632/oncotarget.14406
53. Schimanski CC, Frerichs K, Rahman F, Berger M, Lang H, Galle PR, et al. High miR-196a levels promote the oncogenic phenotype of colorectal cancer cells. WJG. (2009) 15:2089. doi: 10.3748/wjg.15.2089
54. Tang B, Qi G, Sun X, Tang F, Yuan S, Wang Z, et al. HOXA7 plays a critical role in metastasis of liver cancer associated with activation of Snail. Mol Cancer. (2016) 15:57. doi: 10.1186/s12943-016-0540-4
55. Chen H, Zhang Y, Wu K, Qiu X. CircVAPA promotes the proliferation, migration and invasion of oral cancer cells through the miR-132/HOXA7 axis. J Int Med Res. (2021) 49:030006052110132. doi: 10.1177/03000605211013207
56. Liu W, Zhou Y, Duan W, Song J, Wei S, Xia S et al. Glutathione peroxidase 4-dependent glutathione high-consumption drives acquired platinum chemoresistance in lung cancer-derived brain metastasis. Clin Translat Med. (2021) 11:517. doi: 10.1002/ctm2.517
57. Nagpal A, Redvers RP, Ling X, Ayton S, Fuentes M, Tavancheh E, et al. Neoadjuvant neratinib promotes ferroptosis and inhibits brain metastasis in a novel syngeneic model of spontaneous HER2+ve breast cancer metastasis. Breast Cancer Res. (2019) 21:94. doi: 10.1186/s13058-019-1177-1
58. Li J, Li Z, Zhao S, Song Y, Si L, Wang X. Identification key genes, key miRNAs and key transcription factors of lung adenocarcinoma. J Thorac Dis. (2020) 12:1917–1933. doi: 10.21037/jtd-19-4168
59. Zhang Y, Zhang Y, Yin Y, Li S. Detection of circulating exosomal miR-17-5p serves as a novel non-invasive diagnostic marker for non-small cell lung cancer patients. Pathol - Res Prac. (2019) 215:152466. doi: 10.1016/j.prp.2019.152466
60. Wei S, Wang K, Huang X, Tang W, Zhao Z, Zhao Z. Knockdown of the lncRNA MALAT1 alleviates lipopolysaccharide-induced A549 cell injury by targeting the miR-17-5p/FOXA1 axis. Mol Med Report. (2019) 20:2021–9. doi: 10.3892/mmr.2019.10392
61. Shen L, Chen L, Wang Y, Jiang X, Xia H, Zhuang Z. Long noncoding RNA MALAT1 promotes brain metastasis by inducing epithelial-mesenchymal transition in lung cancer. J Neurooncol. (2015) 121:101–108. doi: 10.1007/s11060-014-1613-0
62. Global Analysis of miRNA–mRNA interaction network in breast cancer with brain metastasis. AR. (2017) 37, no. 8,. doi: 10.21873/anticanres.11841
63. Matsubara H, Takeuchi T, Nishikawa E, Yanagisawa K, Hayashita Y, Ebi H, et al. Apoptosis induction by antisense oligonucleotides against miR-17-5p and miR-20a in lung cancers overexpressing miR-17-92. Oncogene. (2007) 26:6099–105. doi: 10.1038/sj.onc.1210425
64. Hu SP, Ge MX, Gao L, Jiang M, Hu KW. LncRNA HCP5 as a potential therapeutic target and prognostic biomarker for various cancers: a meta-analysis and bioinformatics analysis. Cancer Cell Int. (2021) 21:686. doi: 10.1186/s12935-021-02404-x
65. Jiang L, Wang R, Fang LI, Ge X, Chen L, Zhou M, et al. HCP5 is a SMAD3-responsive long non-coding RNA that promotes lung adenocarcinoma metastasis via miR-203/SNAI axis. Theranostics. (2019) 9:2460–74. doi: 10.7150/thno.31097
66. Chen JH, Zhou LY, Xu S, Zheng YL, Wan YF, Hu CP. Overexpression of lncRNA HOXA11-AS promotes cell epithelial–mesenchymal transition by repressing miR-200b in non-small cell lung cancer. Cancer Cell Int. (2017) 17:64. doi: 10.1186/s12935-017-0433-7
67. Gao M, Liu L, Yang Y, Li M, Ma Q, Chang Z. LncRNA HCP5 induces gastric cancer cell proliferation, invasion, and EMT processes through the miR-186-5p/WNT5A Axis Under Hypoxia. Front.Cell Dev. Biol. (2021) 9:663654. doi: 10.3389/fcell.2021.663654
Keywords: HCP5, hsa-miR-17-5p, HOXA7, lung adenocarcinoma, ferroptosis, neurovascular, brain metastasis, biomarker
Citation: Chen Q, Pan Q, Gao H, Wang Y and Zhong X (2022) miR-17-5p/HOXA7 Is a Potential Driver for Brain Metastasis of Lung Adenocarcinoma Related to Ferroptosis Revealed by Bioinformatic Analysis. Front. Neurol. 13:878947. doi: 10.3389/fneur.2022.878947
Received: 18 February 2022; Accepted: 03 May 2022;
Published: 25 May 2022.
Edited by:
Yuzhen Xu, Tongji University, ChinaReviewed by:
Peng-fei Liu, First Affiliated Hospital of Harbin Medical University, ChinaTingyuan Lang, Chongqing University, China
Copyright © 2022 Chen, Pan, Gao, Wang and Zhong. 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: Xiaoning Zhong, uhk3807@163.com