- Research
- Open access
- Published:
Cuproptosis-related ceRNA axis triggers cell proliferation and cell cycle through CBX2 in lung adenocarcinoma
BMC Pulmonary Medicine volume 24, Article number: 85 (2024)
Abstract
Background
Lung adenocarcinoma (LUAD) has high morbidity and mortality. Despite substantial advances in treatment, the prognosis of patients with LUAD remains unfavorable. The ceRNA axis has been reported to play an important role in the pathogenesis of LUAD. In addition, cuproptosis is considered an important factor in tumorigenesis. The expression of CBX2 has been associated with the development of multiple tumors, including LUAD. However, the precise molecular mechanisms through which the cuproptosis-related ceRNA network regulates CBX2 remain unclear.
Methods
The DEGs between tumor and normal samples of LUAD were identified in TCGA database. The “ConsensusClusterPlus” R package was used to perform consensus clustering based on the mRNA expression matrix and cuproptosis-related gene expression profile. Then, LASSO-COX regression analysis was performed to identify potential prognostic biomarkers associated with cuproptosis, and the ceRNA network was constructed. Finally, the mechanisms of ceRNA in LUAD was studied by cell experiments.
Results
In this study, the AC144450.1/miR-424-5p axis was found to promote the progression of LUAD by acting on CBX2. The expression of AC144450.1 and miR-424-5p can be altered to regulate CBX2 and is correlated with cell proliferation and cell cycle of LUAD. Mechanistically, AC144450.1 affects the expression of CBX2 by acting as the ceRNA of miR-424-5p. In addition, a cuproptosis-related model were constructed in this study to predict the prognosis of LUAD.
Conclusions
This study is the first to demonstrate that the AC144450.1/miR-424-5p/CBX2 axis is involved in LUAD progression and may serve as a novel target for its diagnosis and treatment.
Background
Lung adenocarcinoma (LUAD), particularly in advanced stages, is highly malignant and has an unfavorable prognosis [1]. Despite substantial advancements in treatment, the 5-year survival rate of patients with LUAD is 12–15% [2]. Currently available diagnostic methods, such as imaging, biopsy, and molecular testing, have limitations, including the high likelihood of false-negative results and the inability to detect early-stage LUAD [3]. In contrast, novel diagnostic techniques such as next-generation sequencing and companion diagnostics have demonstrated great potential in identifying genomic alterations and biomarkers that can facilitate early detection and prompt personalized treatment [4]. Additionally, novel techniques for detecting early-stage cancer and understanding its progression are under development [5]. Traditional therapeutic strategies such as surgery, chemotherapy, and radiation therapy have disadvantages, including side effects, drug resistance, and the inability to treat advanced-stage LUAD [6]. However, novel therapies such as immunotherapy and targeted therapy, have shown beneficial effects in the treatment of LUAD. Therefore, developing more effective diagnostic and treatment methods is necessary for improving the prognosis of patients with LUAD.
MicroRNAs (miRNAs) modulate gene expression by binding to the 3´-untranslated region (UTR) of target mRNAs. They are small non-coding RNA (ncRNAs) that can either degrade the target mRNA or inhibit translation [7]. miRNAs play a crucial role in regulating the cell cycle, metastasis, angiogenesis, metabolism, and apoptosis, thus influencing tumorigenesis [8]. In lung cancer, miRNAs act as both tumor suppressors and oncomirs, regulating the expression of target mRNAs to influence tumor cell proliferation, angiogenesis, and metastasis [9]. miR-424-5p, a distinct miRNA variant, is involved in regulating gene expression and the progression of several cancer types [10]. For instance, it can hinder the proliferative, migratory, and invasive capabilities of lung cancer cells [11]. Long non-coding RNAs (lncRNAs) are RNA transcripts that exceed 200 nucleotides in length and lack protein-coding ability but play a pivotal role in gene regulation [12]. lncRNAs can interact with DNA, RNA, and proteins to regulate chromatin structure, transcription, post-transcriptional processes, and translation [13]. They act as molecular decoys for miRNAs and regulate splicing, mRNA decay, translation, and protein stability [14]. In addition, they are involved in the development of various types of cancers, including LUAD [15]. For example, the lncRNA MALAT1 induces tumor cell proliferation and metastasis by regulating downstream genes in LUAD [16], and the lncRNA HOTAIR induces tumor cell proliferation and invasion by regulating downstream genes in LUAD [17]. lncRNAs can be used as biomarkers for the diagnosis and treatment of LUAD. The lncRNA AC144450.1, also known as ENSG00000228613, is involved in various biological processes, including cancer development [18]. However, further studies are warranted to examine the precise role of AC144450.1 in biological processes and assess its therapeutic potential. Competitive endogenous RNAs (ceRNAs) is a interaction between RNAs crucial in developing neoplasms [19]. It involves lncRNA acting as an endogenous sponge to effect mRNA expression via sinking miRNA [20]. ceRNAs have been reported to play a vital role in the initiation and development of LUAD [21]. lncRNAs and miRNAs that are inversely related to each other may serve as biomarkers for the diagnosis, prognosis, and treatment of cancer [22].
The presence of excessive copper can induce cell death. The mechanisms underlying copper-induced cell death remained elusive until the identification of cuproptosis [23]. Cuproptosis is a programmed cell death mode induced by copper accumulation and proteotoxic stress [24]. Cuproptosis has been associated with various cancers, including LUAD [25], and cellular metabolism. Certain cancer types often exhibit elevated levels of aerobic respiration [23]. Numerous studies have investigated the relationship between cuproptosis-related genes and diverse tumor characteristics [26]. Some cancers, such as breast cancer and leukemia, have been associated with cuproptosis [27]. Studies have reported the development of risk models based on cuproptosis-related genes to predict prognosis and assess immune cell infiltration in LUAD [28, 29]. In addition, studies have revealed cuproptosis-related genes that can serve as therapeutic targets for LUAD [30].
CBX2 is a gene responsible for encoding an element of the polycomb multiprotein complex that is essential for preserving the transcriptionally repressive state of numerous genes [31]. CBX2 functions as both a transcriptional repressor and activator, thus playing an important role in regulating gene expression, cell differentiation, and cell fate commitment [32]. It has been demonstrated to be upregulated in various cancers, such as breast and colorectal cancers [33]. It is markedly upregulated in LUAD and is associated with an unfavorable prognosis [34]. In addition, it promotes tumor growth and metastasis by regulating gene expression in LUAD [35]. Altogether, CBX2 has emerged as a prospective therapeutic target for LUAD [35]. However, whether ncRNAs regulate CBX2 in LUAD and the underlying regulatory mechanisms remain unclear. Therefore, we attempt to investigate the pontential mechanisms of the cuproptosis-related ceRNA network regulating CBX2.
In this study, AC144450.1 was identified as a ceRNA of miR-424-5p and was found to regulate the cuproptosis-related gene CBX2 in LUAD. Initially, we investigated the regulatory effects of AC144450.1, miR-424-5p, and CBX2 on the malignant progression of LUAD in TCGA database. Subsequently, we predicted and validated the interaction between AC144450.1 and miR-424-5p, demonstrating the role of AC144450.1 as a molecular sponge that regulates CBX2. Mechanistically, the function of miR-424-5p in the malignant progression of LUAD relies on CBX2 to some extent. Altogether, this study introduces a novel and effective approach to diagnosing and treating LUAD.
Methods
Data collection
The RNA-seq data and clinical information of patients with LUAD were extracted from TCGA database. The miRNA microarray dataset GSE62182 was downloaded from the GEO database. Patients with other malignancies and those with incomplete clinical information were excluded.
Differential gene expression analysis
The “limma” package was used to identify differentially expressed genes (DEGs) between lung cancer samples and corresponding adjacent normal tissue samples. Significant DEGs were defined by adjusted P-values of < 0.05 and|fold change (FC)| values of > 1.0.
Functional enrichment analysis
The “clusterProfiler” software package was used to implement GO functional annotation and KEGG pathway enrichment analyses. Significantly enriched pathways were defined by false discovery rates (FDRs) < 0.05. Bubble maps were drawn to visualize the results of enrichment analyses, whereas a chord diagram was drawn using the “circlize” R package to visualize mRNAs in each pathway.
Consistent cluster analysis
The “ConsensusClusterPlus” R package was used to perform consensus clustering based on the mRNA expression matrix and cuproptosis-related gene expression profile. All samples were divided into k groups (k = 2–9). The optimal number of clusters was determined using two methods: the consensus matrix and the cumulative distribution function (CDF). In addition, Kaplan-Meier curves were plotted to compare survival outcomes between subgroups.
Construction of the lncRNA–microRNA–mRNA network
The StarBase and TargetScan databases were used to predict interactions between lncRNAs and mRNAs based on sequence complementarity. The miRDB, miRTarBase, and TargetScan databases were used to predict mRNAs regulated by miRNAs. These databases use algorithms to predict potential miRNA targets based on sequence complementarity and experimental evidence. Finally, a lncRNA-miRNA-mRNA regulatory network based on ceRNAs was constructed by integrating the predicted lncRNA-miRNA and miRNA-mRNA pairs, and the Cytoscape (version 3.7.2) software was used to visualize the ceRNA network.
Analysis of prognostic biomarkers associated with cuproptosis
LASSO-COX regression analysis was performed to identify potential prognostic biomarkers associated with cuproptosis. The risk score of each patient was calculated based on the prognostic model, and the patients were divided into high- and low-risk groups. The lower the risk score, the lower the recurrence risk. Kaplan-Meier curves were plotted to compare survival between the two risk groups. In addition, receiver operating characteristic (ROC) curves were plotted using the “survival ROC” R package.
Cell culture
The human lung cancer cell line A549 and HCC827 used in this study were purchased from the National Authenticated Cell Culture Collection in Shanghai, China. A549 cells were cultured in DMEM supplemented with 10% fetal bovine serum and 1% penicillin–streptomycin (to prevent bacterial contamination) (Gibco, USA). HCC827 were cultured in RPMI-1640 supplemented with 10% fetal bovine serum and 1% penicillin–streptomycin (to prevent bacterial contamination) (Gibco, USA).
Cell transfection
Transfection was performed using Lipofectamine 2000 (Invitrogen) according to the manufacturer’s instructions. Cells were seeded in 6-well plates at an optimal density and transfected with plasmids using Lipofectamine 2000 and Opti-MEM after 24 h. Empty plasmids and plasmids encoding an miR-424-5p mimic and an miR-424-5p inhibitor were procured from GenePharma (Shanghai, China). After 48 h of transfection, the cells were harvested for subsequent experiments.
RNA extraction and qRT-PCR
Total RNA was extracted from cells and reverse transcribed into cDNA using a reverse transcription kit. The synthesized cDNA was amplified via qPCR using the SYBR Premix Ex Taq kit. The relative RNA expression was quantified using the 2-ΔΔCt method. All reagents used for qRT-PCR were purchased from TaKaRa, Japan.
Protein extraction and western blotting
Transfected cells were harvested and lysed in RIPA buffer (Beyotime Biotechnology, China) containing a protease inhibitor mixture. The extracted proteins (in equal amounts) were separated via SDS-PAGE and transferred to PVDF membranes, followed by incubation with primary and secondary antibodies.
CCK8 assay
To evaluate cell viability, cells were seeded in a 96-well plate and cultured until they reached 70–80% confluence. Thereafter, 10 μL of the CCK-8 reagent was added to each well, and the cells were incubated for 2 h. Absorbance was measured at 450 nm using a microplate reader, and the percentage of cell viability was calculated relative to the control group.
Cell cycle
Cells were harvested after transfection, washed twice with precooled PBS, and digested with pancreatin without EDTA. The cells were fixed with 70% ethanol overnight, washed twice with PBS, and stained with 500 μL of PI/RNaSeA solution for 30 min at room temperature. Subsequently, the cell cycle was analyzed on an FACSCalibur flow cytometer (BD).
Luciferase assay
The promoter region of the target gene was cloned into a pGL3-based vector upstream of the firefly luciferase gene. The 3´-untranslated region (UTR) of the target gene downstream of the Renilla luciferase gene was cloned into the pRL-TK vector. Cells were seeded in 24-well plates and cultured until they reached 70–80% confluence. The cells were co-transfected with reporter plasmids and either wild-type (WT) or mutant (MUT) lncRNA AC144450.1 using transfection reagents. Transfected cells were harvested after 24–48 h and washed with PBS. Cell lysates were collected, and the activity of firefly and Renilla luciferases was measured on a photometer.
Statistical analysis
The GraphPad Prism 8 software was used for statistical analysis, and all data were expressed as the mean ± standard deviation. The two-tailed Student’s t-test was used to compare the data of two groups. Survival curves were plotted using the Kaplan-Meier plotter and compared using the log-rank (Mantel-Cox) test. Statistical significance was denoted as follows: *, P < 0.05; **, P < 0.01; ***, P < 0.001.
Results
Identification of lncRNAs, miRNAs, and mRNAs related to the progression of lung cancer
Differentially expressed lncRNAs and mRNAs were identified between 58 matched tumor and normal samples in TCGA-LUAD dataset. A volcano map was constructed to identify significant DEGs based on P-values of < 0.05 and log2FC values of ≥ 1 or ≤ -1 after adjusting for multiple comparisons. The expression of both mRNAs and lncRNAs was significantly different between tumor and normal tissues. Specifically, 1980 upregulated and 2362 downregulated mRNAs as well as 1184 upregulated and 1476 downregulated lncRNAs were identified in lung cancer tissues (Fig. 1A-B and Supplementary Fig. 1A-C). Furthermore, differentially expressed miRNAs were screened in the GSE62182 dataset using GEO2R. A total of 24 miRNAs were upregulated and 89 miRNAs were downregulated in lung cancer tissues. These DEGs were visualized on heatmaps generated using the ggplot2 R package. Tumor and non-tumor tissues were separately clustered based on their respective lncRNA, miRNA, and mRNA expression patterns.
Pathway enrichment analysis of genes related to the progression of lung cancer
Gene set enrichment analysis (GSEA) was performed to examine the significance of the identified DEGs in lung cancer. KEGG analysis revealed that the DEGs were mainly enriched in pathways related to cell proliferation and DNA replication (Fig. 2A and Supplementary Fig. 2A-C). Furthermore, GO analysis revealed that the DEGs were primarily involved in cell proliferation, cell cycle transition, activation of cell cycle-related signaling pathways, and DNA replication (Fig. 2B and Supplementary Fig. 2D-F). Altogether, the results indicated that genes regulating the malignant progression of lung cancer were involved in cell proliferation and cell cycle.
Consensus cluster analysis of cuproptosis-related genes
To examine the potential role of cuproptosis in lung cancer, we performed cluster analysis to evaluate the expression level of cuproptosis of lung cancer in TCGA database, The results revealed that there are three groups of differential expression levels of cuproptosis (Fig. 3A and Supplementary Fig. 3A-B). The group 2 and 3 were choosed for further differential gene analysis according to the survival result (Fig. 3B). A total of 1325 DEGs were identified, including 611 upregulated and 714 downregulated genes (Fig. 3C-D and Supplementary Fig. 3C). GSEA revealed that these DEGs were primarily involved in pathways related to cell proliferation, cell cycle, and other cellular processes (Fig. 3E and Supplementary Fig. 3D).
lncRNA AC144450.1/miR-424-5p/CBX2 may be an important ceRNA axis in the regulation of lung cancer
To gain deeper insights into the regulatory role of lncRNAs in mediating miRNA-mRNA interactions in lung cancer, a ceRNA network was constructed including the differentially expressed lncRNAs, miRNAs, and mRNAs involved in regulating the malignant progression of lung cancer. According to the ceRNA hypothesis, a positive regulatory relationship exists between lncRNAs and mRNAs, whereas a negative regulatory relationship exists between miRNAs and mRNAs. Eventually, 5 lncRNAs, 6 miRNAs, 18 mRNAs, and 34 interaction axes were identified, which were visualized using the Cytoscape (version 3.9.0) software (Fig. 4A and Supplementary Fig. 4). AC144450.1 and miR-424-5p were identified as the lncRNA and miRNA, respectively, involved in the ceRNA axis regulating CBX2 (Fig. 4B). The previous studies found that CBX2 acts as a transcription factor regulating numerous cell functions, which means CBX2 exerts important role in cancer development. Therefore, the AC144450.1/miR-424-5p/CBX2 axis was selected for further experiments.
lncRNA AC144450.1/miR-424-5p/CBX2 axis may serve as a prognostic biomarker for lung cancer
Previous studies have demonstrated that AC144450.1, miR-424-5p, and CBX2 are potential prognostic biomarkers for lung cancer. Because the malignant progression of lung cancer is the result of the collective regulatory effects of multiple gene networks, using a single gene as a diagnostic marker has limitations. To investigate the potential prognostic value of the identified biomarkers, LASSO regression analysis (Fig. 5A-E) was performed to construct a cuproptosis-related prognostic signature incorporating the three biomarkers. According to this signature, the risk scores of patients with LUAD were calculated using the following formula: LINC01991 * (-0.124) + FSCN1 * 0.0007 + EGLN3 * 0.0014 + FOSL13 * 0.0006. The patients were divided into high- and low-risk groups based on the median risk score. The risk scores were significantly associated with survival in both groups. The prognosis of high-risk patients was poorer than that of low-risk patients, and the AUC value of the polygenic model was higher than that of the single-gene model (Fig. 5F and J, Supplementary Fig. 5A-D). Subsequently, the high predictive performance of the prognostic model was verified in the validation set (Fig. 5G-I). Altogether, the prognostic model constructed based on AC144450.1, miR-424-5p, and CBX2 detected lung cancer and predicted the prognosis more accurately.
lncRNA AC144450.1/miR-424-5p/CBX2 regulates the malignant progression of lung cancer
According to the abovementioned results, the AC144450.1/miR-424-5p/CBX2 axis functions as the ceRNA network that regulates the occurrence and development of lung cancer. However, the precise role of AC144450.1, miR-424-5p, and CBX2 in the malignant progression of lung cancer remains unclear. Therefore, we divided the samples into high- and low-expression groups based on the median expression of the three factors and performed differential expression and enrichment analyses for each group. The results revealed that miRNA- and mRNA-regulated genes were enriched in pathways associated with cell proliferation and cell cycle (Fig. 6A and Supplementary Fig. 6A-B). Besides, we analysed the differentially expressed levels of the three genes in EGFR wild type and mutation type groups in TCGA LUAD patients, the results showed that CBX2 is significantly higher expression comparing with the EGFR mutation group (Supplementary Fig. 6C). Next, the expression of AC144450.1, miR-424-5p, and CBX2 was individually knocked down in lung cancer cells to examine the regulatory effects of the ceRNA network on cell proliferation and cell cycle. The results of CCK8 assay, colony formation assay, and cell cycle analysis revealed that the knockdown of each gene led to a significant reduction in cell proliferation and altered the cell cycle of lung cancer cells (Fig. 6B-D and Supplementary Fig. 7A-C). These results were consistent with those of enrichment analyses. Altogether, the results revealed that AC144450.1, miR-424-5p, and CBX2 individually regulated the malignant progression of lung cancer.
lncRNA AC144450.1 and mir-424-5p regulate the CBX2 protein in lung cancer
To examine the relationship among AC144450.1, miR-424-5p, and CBX2 proteins, we used antisense oligonucleotides (ASOs) to knock down AC144450.1 and assessed changes in the expression of miR-424-5p, CBX2 using western blotting. The results revealed that AC144450.1 knockdown decreased the expression of AC144450.1, resulting in subsequent downregulation of CBX2 and upregulation of miR-424-5p (Fig. 7A and Supplementary Fig. 8A). Furthermore, miR-424-5p mimic and inhibitor were used to overexpress and inhibit it, respectively, and changes in the expression of AC144450.1 and CBX2 protein were detected. The results revealed that the signaling pathway regulated by AC144450.1 and CBX2 was negatively correlated with the expression of miR-424-5p (Fig. 7B-C and Supplementary Fig. 8B-C). Altogether, these results indicate that AC144450.1 and miR-424-5p regulate the CBX2 protein in lung cancer, which is consistent with the ceRNA hypothesis.
AC144450.1 as the ceRNA of miR-424-5p regulates the CBX2 protein and affects the malignant progression of lung cancer
A complementary mutation was induced in AC144450.1 and the binding site of miR-424-5p was predicted to assess whether AC144450.1 is the ceRNA of miR-424-5p (Fig. 8A). Dual-luciferase reporter assay was used to detect changes in luciferase activity caused by changes in miR-424-5p levels before and after AC144450.1 mutation induced. The results revealed that miR-424-5p regulatedthe luciferase activity of WT AC144450.1 but not that of mutant AC144450.1 (Fig. 8B and Supplementary Fig. 8D). Thereafter, we predicted the enrichment of genes co-regulated by two or more genes in the AC144450.1/miR-424-5p/CBX2 axis on the lung cancer malignant progression pathway (Fig. 8C). The intersection genes of miR-424-5p and CBX2 were significantly enriched in the proliferation-related pathway (Fig. 8D). To validate this finding, CCK8 and clone formation assays were performed to examine cell proliferation after inhibiting miR-424-5p in shNC or shCBX2 A549 cells. The results showed that the proliferation was activated by miR-424-5p inhibitor in shNC group. However, when CBX2 was knocked down with shRNAs, the effects of the miR-424-5p inhibitor on proliferation were rescued (Fig. 8E-F and Supplementary Fig. 8E-F). Altogether, these results suggest that miR-424-5p acts as a regulator of CBX2 and affects the malignant behavior of lung cancer.
Discussion
LUAD is a highly malignant tumor with a poor prognosis. The development of more effective diagnostic and treatment strategies for LUAD is a major focus of research at present. miRNAs are small ncRNAs that play an important role in tumor progression by regulating the cell cycle, metastasis, angiogenesis, and metabolism. They regulate gene expression and are involved in the development, progression, metastasis, diagnosis, and prognosis of cancer. For instance, the downregulation of miR-7, which has an aberrant expression pattern in lung cancer, has been demonstrated to inhibit the mRNA expression of EGFR and Raf1 [36]. Similar to miRNAs, lncRNAs play an important role in LUAD. They regulate the proliferative, invasive, and migratory capabilities of LUAD cells by interacting with miRNAs, WNT pathway, and other mechanisms [37]. The abnormal expression of lncRNAs has been associated with metastasis, advanced pathological stages, and an unfavorable prognosis among patients with LUAD. For example, dysregulated lncRNAs such as ZNF503-AS1, CYP4F26P, and RP11-108M12.3 have been associated with a poor prognosis in LUAD [38]. Therefore, lncRNAs may serve as potential targets for developing novel therapeutic strategies for LUAD. ceRNAs represent a unique mechanism of gene regulation and play an important role in the physiology and pathogenesis of diseases, including cancer. ceRNA networks of mRNAs, lncRNAs, and miRNAs have been identified in LUAD [15]. These networks regulate mRNA expression and interactions among the marker genes of tumor-associated macrophages. In this study, we identified a ceRNA regulatory network in LUAD and verified that AC144450.1, a ceRNA of miR-424-5P, regulates the malignant progression of LUAD cells. The findings of this study add to the existing knowledge of the ceRNA regulatory network in LUAD.
Cuproptosis is a copper-dependent cell death mechanism that differs from other existing cell death mechanisms. It is closely associated with cellular metabolism. Specific cancer types often exhibit elevated aerobic respiration. Studies have extensively investigated the relationship between cuproptosis-related genes and diverse tumor characteristics in multiple cancers. These studies have revealed the significant role of cuproptosis-related genes in signaling pathways associated with tumorigenesis [23]. For treating early lung cancer, cuproptosis could potentially be a valuable therapeutic approach [39]. Research into the role of cuproptosis in the treatment of early-stage lung cancer has focused on understanding the mechanisms of this form of cell death and developing prognostic models based on cuproptosis-related biomarkers [40]. Due to the significant impact of cuproptosis-related genes and pathways on mitochondrial function, researchers have focused on analyzing the prognostic value through bioinformatics and constructing prognostic models to achieve early diagnosis at the gene level [41]. A study explored the role of 16-cuproptosis-related long non-coding RNAs in predicting the clinical prognosis of patients with lung adenocarcinoma [42]. Further research on cuproptosis holds promise in enhancing our understanding of cancer and guiding clinical interventions. A notable application of cuproptosis in tumor diagnosis and treatment is the construction of risk or prognostic models [43]. The relationship between cuproptosis and LUAD has attracted substantial interest, and researchers have used online databases to identify genes involved in this relationship [44]. LUAD subtypes associated with cuproptosis have been identified and associated with clinicopathological characteristics, prognosis, and biological pathways [45]. In a study, weighted gene co-expression network analysis was used to establish a cuproptosis-based risk model for predicting the prognosis of LUAD [46]. Cuproptosis-related genes play an essential role in diagnosing and treating LUAD. In this study, the LUAD samples were uniformly clustered based on the expression of 13 cuproptosis-related genes. Patients with LUAD were divided into two groups based on the levels of cuproptosis, and differential expression analysis was performed to identify significant cuproptosis-related genes associated with the prognosis of LUAD. Altogether, the results of this study provide a potential therapeutic target for LUAD.
The CBX2 gene is responsible for maintaining the transcriptionally repressive state of numerous genes, which is crucial for regulating transcription and negatively modulating cell differentiation. CBX2 regulates the expression of various genes involved in embryonic development, sexual determination, stem cell differentiation, and primary sexual characteristics. A recent study revealed that CBX2 inhibits cancer cell proliferation and invasion by regulating the Akt/PI3K pathway [47]. Knockdown of CBX2 leads to a significant reduction in tumor cell proliferation and promotes anoikis, which increases the sensitivity of tumor cells to pods [48]. The abnormally high expression of CBX2 has been closely associated with the degree of malignancy, TNM stage, lymph node metastasis, and poor patient prognosis. Additionally, CBX2 may also influence the sensitivity of cancer cells to chemotherapeutic drugs [49]. In a study, in vitro and in vivo experiments demonstrated that knockdown of CBX2 remarkably prevented cell growth and metastasis in LUAD [35]. In addition, a similar study reported that concurrent high expression of CBX2 and EZH2 promoted the growth and metastasis of LUAD [50]. Altogether, these studies suggest that CBX2 functions as an oncogene in LUAD and serves as a promising therapeutic target for the disease. Furthermore, the involvement of CBX2 in ceRNA networks has been reported in diverse cancer types [51]. The ceRNA regulatory network involving CBX2, PBK, and AP002478.1 influences the progression of hepatocellular carcinoma [52]. CBX2 may contribute to cancer progression through ceRNA networks, suggesting that it is a promising therapeutic target for cancer. However, no studies have reported the involvement of CBX2 in ceRNA networks associated with LUAD. In this study, we established a ceRNA regulatory network involving CBX2, AC144450.1, and miR-424-5p. This network was found to be associated with the progression of LUAD. These findings provide a theoretical basis for the development of CBX2-targeted diagnostic and therapeutic strategies for LUAD.
Conclusions
In conclusion, this study revealed that the lncRNA AC144450.1/miR-424-5p/CBX2 axis is a novel ceRNA network that regulates the malignant progression of LUAD. AC144450.1 functions as a molecular sponge of miR-424-5p and competitively regulates the expression of CBX2. In addition, the regulatory effects of miR-424-5p on the malignant progression of LUAD rely on CBX2 to a certain extent. These results provide novel insights into the complex regulatory mechanisms underlying the pathogenesis of LUAD and introduce novel targets for the diagnosis and treatment of the disease.
Data availability
The data generated or analysed during this study are included in this published article and any related materials are available from the corresponding author on reasonable request.
Abbreviations
- CBX2:
-
Chromobox protein homolog 2
- ceRNAs:
-
Competitive endogenous RNAs
- GSEA:
-
Gene set enrichment analysis
- lncRNAs:
-
Long non-coding RNAs
- LUAD:
-
Lung adenocarcinoma
- miRNAs:
-
MicroRNAs
- ncRNAs:
-
Non-coding RNA
- ROC:
-
Receiver operating characteristic
- UTR:
-
Untranslated region
- WT:
-
Wild-type
References
Rudin CM, Brambilla E, Faivre-Finn C, Sage JJNRDP. Small-cell lung cancer. 2021;7:3.
Myers DJ, Wallen JM. Lung adenocarcinoma. StatPearls. StatPearls Publishing; 2022.
Kazdal D, Hofman V, Christopoulos P, Ilié M, Stenzinger A, Hofman PJG, Chromosomes et al. Fusion-positive non‐small cell lung carcinoma: biological principles, clinical practice, and diagnostic implications. 2022;61:244–60.
Bennett CW, Berchem G, Kim YJ, El-Khoury VJO. Cell-free DNA and next-generation sequencing in the service of personalized medicine for lung cancer. 2016;7:71013.
Crosby D, Bhatia S, Brindle KM, Coussens LM, Dive C, Emberton M et al. Early detection of cancer. 2022;375:eaay9040.
Nooreldeen R, Bach HJI. Current and future development in lung cancer diagnosis. 2021;22:8661.
Statello L, Guo C-J, Chen L-L. Huarte MJNrMcb. Gene regulation by long non-coding RNAs and its biological functions. 2021;22:96–118.
Chen Y, Gao D-Y, Huang LJA. In vivo delivery of miRNAs for cancer therapy: challenges and strategies. 2015;81:128–41.
Lee SS, Cheah YKJJ. The interplay between microRNAs and cellular components of tumour microenvironment (TME) on non-small-cell lung cancer (NSCLC) progression. 2019;2019.
Guo C, Ju Q-q, Zhang C-x, Gong M, Li Z-l. Gao Y-yJBc. Overexpression of HOXA10 is associated with unfavorable prognosis of acute myeloid leukemia. 2020;20:1–12.
Zhang M, Zeng J, Zhao Z, Liu ZJM. Loss of MiR-424‐3p, not miR‐424‐5p, confers chemoresistance through targeting YAP1 in non‐small cell lung cancer. 2017;56:821–32.
Zhang X, Wang W, Zhu W, Dong J, Cheng Y, Yin Z et al. Mechanisms and functions of long non-coding RNAs at multiple regulatory levels. 2019;20:5573.
Briata P, Gherzi RJN-cR. Long non-coding RNA-ribonucleoprotein networks in the post-transcriptional control of gene expression. 2020;6:40.
Li Z, Zhao W, Wang M, Zhou XJGEPC. The role of long noncoding RNAs in gene expression regulation. 2019:1–17.
Wu X, Sui Z, Zhang H, Wang Y, Yu ZJFO. Integrated Analysis of lncRNA–Mediated ceRNA network in lung adenocarcinoma. 2020;10:554759.
Tiansheng G, Junming H, Xiaoyun W, Peixi C, Shaoshan D, Qianping CJCJ. lncRNA metastasis-associated lung adenocarcinoma transcript 1 promotes proliferation and invasion of non-small cell lung cancer cells via down-regulating miR-202 expression. 2020;22:375.
Liu W, Yin N-C, Liu H, Nan K-JJG. Cav-1 promote lung cancer cell proliferation and invasion through lncRNA HOTAIR. 2018;641:335–40.
Zhu H, Zhou Y, Wang Q, Yang X, Ding C, Xiong YJB et al. Long non-coding RNA LALTOP promotes non-small cell lung cancer progression by stabilizing topoisomerase IIα mRNA. 2021;574:56–62.
Bai Q, Pan Z, Nabi G, Rashid F, Liu Y, Khan SJAJCR. Emerging role of competing endogenous RNA and associated noncoding RNAs in thyroid cancer. 2022;12:961.
Zhang Z, Cheng J, Wu Y, Qiu J, Sun Y, Tong XJM. LncRNA HOTAIR controls the expression of Rab22a by sponging miR-373 in ovarian cancer. 2016;14:2465–72.
Sui J, Li Y-H, Zhang Y-Q, Li C-Y, Shen X, Yao W-Z et al. Integrated analysis of long non-coding RNA–associated ceRNA network reveals potential lncRNA biomarkers in human lung adenocarcinoma. 2016;49:2023–36.
Vecera M, Sana J, Lipina R, Smrcka M, Slaby OJI. Long non-coding RNAs in gliomas: from molecular pathology to diagnostic biomarkers and therapeutic targets. 2018;19:2754.
Xie J, Yang Y, Gao Y, He JJMC. Cuproptosis: mechanisms and links with cancers. 2023;22:46.
Chi H, Peng G, Wang R, Yang F, Xie X, Zhang J et al. Cuprotosis programmed-cell-death-related lncRNA signature predicts prognosis and immune landscape in PAAD patients. 2022;11:3436.
Liu T, Cai L, Hua H, Jiang X, Xu X, Zhang T et al. Cuprotosis patterns are associated with tumor mutation burden and immune landscape in lung adenocarcinoma2022; 2022.
Sha S, Si L, Wu X, Chen Y, Xiong H, Xu Y et al. Prognostic analysis of cuproptosis-related gene in triple-negative breast cancer. 2022;13:922780.
Zhao Q, Qi TJFO. The implications and prospect of cuproptosis-related genes and copper transporters in cancer progression. 2023;13:1117164.
Ma S, Zhu J, Wang M, Zhu J, Wang W, Xiong Y et al. A cuproptosis-related long non-coding RNA signature to predict the prognosis and immune microenvironment characterization for lung adenocarcinoma. 2022;11:2079.
Zhang P, Pei S, Liu J, Zhang X, Feng Y, Gong Z et al. Cuproptosis-related lncRNA signatures: Predicting prognosis and evaluating the tumor immune microenvironment in lung adenocarcinoma. 2023;12:1088931.
Zhang H, Shi Y, Yi Q, Wang C, Xia Q, Zhang Y et al. A novel defined cuproptosis-related gene signature for predicting the prognosis of lung adenocarcinoma. 2022;13:975185.
Kim JJ, Kingston REJNRG. Context-specific polycomb mechanisms in development. 2022;23:680–95.
Morey L, Santanach A, Blanco E, Aloia L, Nora EP, Bruneau BG et al. Polycomb regulates mesoderm cell fate-specification in embryonic stem cells through activation and repression mechanisms. 2015;17:300–15.
Kiani S, Akhavan-Niaki H, Fattahi S, Kavoosian S, Jelodar NB, Bagheri N et al. Purified sulforaphane from broccoli (Brassica oleracea var. italica) leads to alterations of CDX1 and CDX2 expression and changes in miR-9 and miR-326 levels in human gastric cancer cells. 2018;678:115–23.
Zheng S, Lv P, Su J, Miao K, Xu H, Li MJA. Overexpression of CBX2 in breast cancer promotes tumor progression through the PI3K/AKT signaling pathway. 2019;11:1668.
Hu F-F, Chen H, Duan Y, Lan B, Liu C-J, Hu H et al. CBX2 and EZH2 cooperatively promote the growth and metastasis of lung adenocarcinoma. 2022;27:670–84.
Weng W, Wei Q, Toden S, Yoshida K, Nagasaka T, Fujiwara T et al. Circular RNA ciRS-7—a promising prognostic biomarker and a potential therapeutic target in colorectal cancer. 2017;23:3918–28.
Chen W, Zhao W, Zhang L, Wang L, Wang J, Wan Z et al. MALAT1-miR-101-SOX9 feedback loop modulates the chemo-resistance of lung cancer cell to DDP via wnt signaling pathway. 2017;8:94317.
Dong X, Shao C, Tu J, Chen D, Xu S. Construction and validation of a prognostic signature based on anoikis-related lncRNAs in lung Adenocarcinoma. 2023.
Wang M, Zheng L, Ma S, Lin R, Li J, Yang S. Cuproptosis: emerging biomarkers and potential therapeutics in cancers. Front Oncol. 2023;13:1288504.
Ma S, Zhu J, Wang M, Zhu J, Wang W, Xiong Y, et al. A cuproptosis-related long non-coding RNA signature to predict the prognosis and immune microenvironment characterization for lung adenocarcinoma. Transl Lung Cancer Res. 2022;11:2079–93.
Li Q, Wang T, Zhou Y, Shi J. Cuproptosis in lung cancer: mechanisms and therapeutic potential. Mol Cell Biochem. 2023.
Wang F, Lin H, Su Q, Li C. Cuproptosis-related lncRNA predict prognosis and immune response of lung adenocarcinoma. World J Surg Oncol. 2022;20:275.
Chen Y, Chen X, Wang XJFC, Biology D. Identification of a prognostic model using cuproptosis-related genes in uveal melanoma. 2022;10:973073.
Jiang G, Song C, Wang X, Xu Y, Li H, He Z et al. The multi-omics analysis identifies a novel cuproptosis-anoikis-related gene signature in prognosis and immune infiltration characterization of lung adenocarcinoma. 2023;9.
Wang Y, Zhang C, Ji C, Jin W, He X, Yu S et al. Molecular subtypes based on cuproptosis-related genes and immune profiles in lung adenocarcinoma. 2022;13:1006938.
Ling X, Zhang L, Fang C, Liang H, Zhu J, Ma JJTLCR. Development of a cuproptosis-related signature for prognosis prediction in lung adenocarcinoma based on WGCNA. 2023;12:754.
Wang L, Ren B, Zhuang H, Zhong Y, Nan YJTiCR. Treatment. CBX2 induces glioma cell proliferation and invasion through the Akt/PI3K pathway. 2021;20:15330338211045831.
Wheeler LJ, Watson ZL, Qamar L, Yamamoto TM, Post MD, Berning AA et al. CBX2 identified as driver of anoikis escape and dissemination in high grade serous ovarian cancer. 2018;7:92.
Chen WY, Zhang XY, Liu T, Liu Y, Zhao YS, Pang DJO. Chromobox homolog 2 protein: a novel biomarker for predicting prognosis and taxol sensitivity in patients with breast cancer. 2017;13:1149–56.
Entezari M, Taheriazam A, Paskeh MDA, Sabouni E, Zandieh MA, Aboutalebi M et al. The pharmacological and biological importance of EZH2 signaling in lung cancer. 2023;160:114313.
Yan Y, Lu Y, Mao K, Zhang M, Liu H, Zhou Q et al. Identification and validation of a prognostic four-genes signature for hepatocellular carcinoma: integrated ceRNA network analysis. 2019;13:618–30.
Gu X, Li H, Sha L, Zhao WJBRI. Construction and comprehensive analyses of a competing endogenous RNA network in tumor-node-metastasis stage I hepatocellular carcinoma. 2020;2020.
Acknowledgements
The authors would like to thank all members of the Liu laboratory.
Funding
This work was supported by the National Natural Science Foundation of China for Young Scholars grants 82203317 (to Quanxing Liu); Medical and Scientific Research Project of Science and Health Joint Department of Chongqing, Grant/Award Number: 2020MSXM031 (to Quanxing Liu); Chongqing Technology Innovation and Application Development Key Projects, Grant/Award Number: CSTB2022TIAD-CUX0019; Clinical Research Foundation of Army Medical University, Grant/Award Number: 2018JSLC033; 2018XLC2014; 2018XLC3035; National Natural Science Foundation of China, Grant/Award Number: 81972190; 81702247.
Author information
Authors and Affiliations
Contributions
JW collected and analyzed the data. GF and CL performed the cell experiments. LC was in charge of data checking. QL designed and supervised the whole study, and was a major contributor in writing the manuscript. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Below is the link to the electronic supplementary material.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Wu, J., Fu, G., Luo, C. et al. Cuproptosis-related ceRNA axis triggers cell proliferation and cell cycle through CBX2 in lung adenocarcinoma. BMC Pulm Med 24, 85 (2024). https://doi.org/10.1186/s12890-024-02887-0
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12890-024-02887-0