Skip to main content

Prognostic and immune infiltrative biomarkers of CENPO in pan–cancer and its relationship with lung adenocarcinoma cell proliferation and metastasis

Abstract

Background

The centromere protein O (CENPO) is an important member of the centromere protein family. However, the role of CENPO in pan–cancer and immune infiltration has not been reported. Here, we investigated the role of CENPO in pan–cancer and further validated its role in lung adenocarcinoma (LUAD) by in vitro experiments.

Method

The UCSC Xena database and The Cancer Genome Atlas (TCGA)–LUAD data were used to assess the expression levels of CENPO. The potential value of CENPO as a diagnostic and prognostic biomarker for pan–cancer was evaluated using TCGA data and the GEPIA database. The -expression profiles of LUAD patients and the corresponding clinical data were downloaded for correlation analysis. The role of CENPO in immune infiltration was investigated using the UCSC Xena database. Subsequently, qRT–PCR was performed to detect the expression of CENPO. Cell proliferation, migration, and invasion were determined using CCK–8, wound–healing assay, and transwell assay, respectively.

Results

CENPO is highly expressed in most cancers, and the upregulation of CENPO is associated with poor prognosis in many cancers. CENPO expression correlates with age, TNM stage, N stage, T stage, and receipt of radiotherapy in LUAD patients, and LUAD patients with high CENPO expression have poorer overall survival (OS) and disease–free survival (DFS). In addition, CENPO expression is associated with immune cell infiltration and immune checkpoint inhibitors. Moreover, the expression of CENPO was closely related to the expression of tumor mutational load and microsatellite instability. In vitro experiments showed that CENPO expression was increased in LUAD cell lines and that knockdown of CENPO significantly inhibited the proliferation, cell invasion, and migration ability of LUAD cells.

Conclusion

CENPO may be a potential pan–cancer biomarker and oncogene, especially in LUAD. In addition, CENPO is associated with immune cell infiltration and may serve as a new molecular therapeutic target and effective prognostic marker for LUAD.

Peer Review reports

Introduction

Cancer is a leading cause of death, with 19.3 million new cancer cases and 10 million cancer deaths worldwide in 2020 [1]. The burden of cancer morbidity and mortality is rapidly escalating [2]. Despite the expanding number of treatments, including improvements in surgical modalities, immunotherapy, targeted therapies, and combination therapies [3], mortality from liver and lung cancers continues to increase yearly [4]. With the rapid development of next–generation sequencing, pan–cancer research has been widely used to identify tumor molecular markers and signaling pathways for a more comprehensive and in–depth understanding of the molecular mechanisms of tumor development [5,6,7,8]. Therefore, searching for novel diagnostic and prognostic biomarkers for cancer can help improve the prognosis of cancer patients.

Centromere protein O (CENPO), a member of the centromere protein (CENP) family, plays an important role in meiosis and mitosis, especially in chromosome segregation and normal cell division processes [9]. Many CENP family members have been identified, including CENPA, CENPL, CENPI and CENPU, and CENPO, which are associated with poor prognosis [6, 9,10,11,12]. CENPO–like proteins are associated with the prevention of premature sister chromatid separation and are associated with the cell cycle [13]. Recent studies have demonstrated that breast cancer patients with high CENPO expression have a poor prognosis and are independent factors affecting the distant recurrence–free survival (DRFS) of breast cancer patients [14]. High CENPO expression promotes the proliferation of gastric cancer cells and is positively associated with poor prognosis in gastric cancer [12]. CENPO is highly expressed in CRC (Colorectal Cancer) tissues and cells as positively associated with the malignancy of CRC degree and positively correlated with tumor progression possibly through EMT and PI3K/AKT signaling pathways [15]. CENPO upregulation promoted bladder cancer progression [16]. However, an investigation into the role of CENPO in pan–cancer is still lacking. We found that CENPO was significantly upregulated in multiple cancers and correlated with lung adenocarcinoma (LUAD) staging, overall survival (OS), and DFS. Thus, we chose LUAD to validate our bioinformatics results.

This study assessed CENPO expression in various human tumors using The Cancer Genome Atlas (TCGA) RNA–seq. The prognostic impact of aberrant CENPO expression on pan–cancer was investigated using the GEPIA database, and the impact of CENPO expression on the prognosis of LUAD patients was analyzed in conjunction with clinical data. Functional enrichment and pathway analysis were then performed, and correlations with immune cell infiltration, tumor mutational load (TMB), microsatellite instability (MSI), and gene methylation were investigated to understand further the biological mechanisms of CENPO in the pathogenesis of LUAD. Subsequently, bioinformatics results were validated by in vitro experiments. We preliminarily explored the potential mechanism of action of CENPO in LUAD, and CENPO could be used as an immunotherapeutic target and prognostic marker for pan–cancer, especially LUAD.

Methods

Material and methods

CENPO mRNA expression analysis

The UCSC Xena database GDC–TCGA data were used to analyze the differential mRNA expression of CENPO in pan–cancer and the relationship between CENPO, and cancer stage was analyzed by the “Stage plots” module of R. Lung adenocarcinoma RNA–seq sequencing profiles were obtained from the UCSC Xena database TCGA–LUAD data (517 LUAD cases and 58 normal cases), (https://xenabrowser.net/datepages/). We first mined the TCGA–LUAD dataset for differentially expressed genes (DEGs) in LUAD samples vs normal samples by the R limma package and then presented DEGs as a volcano map using Volcano map [17] (the criteria for determining differential gene expression were P < 0.05 and logFC absolute value > 1). Cases with insufficient or missing data were removed from subsequent data processing.

Diagnostic value and prognostic analysis of CENPO

The diagnostic value of CENPO in pan–cancer was determined by plotting ROC curves from the GDC–TCGA data of the UCSC Xena database. AUC > 0.8 was thought as a high diagnostic value. In addition, the potential prognostic value of CENPO in pan–cancer was analyzed by overall survival (OS) and disease–free survival (DFS) from the GEPIA database (http://gepia.cancer-pku.cn/). The GEPIA database was divided into high and low expression groups based on the median CENPO expression in 33 cancers. KM survival curve analysis was performed based on the overall survival and disease-free survival of patients in the above groups, respectively, and risk ratios were calculated based on the Cox PH model.

Methylation analysis and tumor mutation load

To screen methylation–regulated genes, we analyzed the correlation between CENPO expression and methylation in pan–cancer using GDC–TCGA data from the UCSC Xena database. Tumor mutational load (TMB) and microsatellite instability (MSI) analyses were performed on 33 tumors in the UCSC Xena database GDC–TCGA data via the “fmsb” package of R4.2.0.3.4. Next, the data were visualized using the radar chart package, and the obtained TMB correlation coefficients were all between 0.8 and -0.8, and the obtained MSI correlation coefficients were all between 0.4 and -0.4.

Correlation of CENPO expression with clinical factors of LUAD and prognostic analysis

The correlation between CENPO expression and clinical characteristics, including gender, age, TNM stage, T stage, N stage, distant metastasis, and history of radiation therapy, were analyzed in 517 patients with LUAD in the TCGA database. To further evaluate the prognostic value of CENPO in patients with LUAD, univariate and multifactorial COX proportional risk regression analyses were performed to investigate the effects of CENPO and various clinical characteristics on OS in patients with LUAD, and HR and 95% confidence intervals were evaluated. In addition, we constructed column line plots based on CENPO expression values and pathological staging to facilitate the application of CENPO in the clinical prognostic assessment of LUAD. Finally, we assessed the predictive accuracy of the column line plots by calibration curves.

Functional enrichment and immune infiltration analysis

Functional enrichment analysis of DEGs was performed by clusterProfiler in R software for Gene Ontology (GO) with the Kyoto Encyclopedia of Genes and Genomes (KEGG) [18, 19]. GO is commonly used to annotate genes and analyze their biological processes, which helps to reveal the functional mechanisms behind the patterns observed in transcriptomics, genomics, and proteomics data, which are classified into cellular components (CC), molecular functions (MF), and biological processes (BP). Furthermore, KEGG is a reference knowledge base for the biological interpretation of large–scale molecular datasets. From this, key points of pathway mechanisms involved in differential genes in living systems can be obtained. Stromal cell score (StromalScore), immune cell score (ImmuneScore), and combined stromal and immune cell score (ESTIMATEScore) were performed by R4.2.0 for each tumor sample of TCGA–LUAD data in the UCSC Xena database, and the expression of CENPO was analyzed for its correlation with the tumor microenvironment correlation, immune cells, immune checkpoints.

Drug sensitivity analysis

The list of cgp2016ExprRma drugs was obtained by the R package, CENPO was divided into high and low expression groups by median, and CENPO expression in LUAD was analyzed with a drug sensitivity analysis (IC50). For filtered conditions, the P–value was less than 0.001.

Cell culture condition

The human normal lung epithelial cell line Beas–2b, LUAD cell lines XWLC–05, H838, A549, SPCA–1, and H1299 (Cell Bank of Chinese Academy of Sciences, Kunming, China) were placed in an incubator at 37 °C containing 5% CO2 using Roswell Park Memorial Institute (RPMI)– 1640 medium (GIBCO) and 10% fetal bovine serum (FBS, GIBCO).

RNA interference and quantitative real–time PCR

GEMA chemically synthesized CENPO–specific small interfering RNAs (CENPO siRNAs), negative control siRNAs, and three RNA interference target sequences were designed (Table 1). Subsequently, cells were inoculated in 6–well plates and transfected when cell fusion reached approximately 50%. CENPO expression plasmids were transfected into cells using Lipofectamine 2000 (Invitrogen) transfection reagent according to the manufacturer's instructions. After 6 h of transfection, the cells were harvested 48 h later, and qRT–PCR verified the transfection efficiency. According to the kit instructions, RNA was extracted from A549 and H1299 cells 48 h after transfection using a Trizol reagent (Invitrogen, Carlsbad, CA, USA). An iScript cDNA Synthesis Kit (Bio–Rad, USA) was used to reverse transcribe the RNA to cDNA, followed by iTaq Universal SYBR (Bio–Rad, USA). Finally, iTaq Universal SYBR (Bio–Rad, USA) was used for qRT–PCR. The relative mRNA expression of CENPO was quantified using the cycling threshold (Ct) value and normalized using the 2−ΔΔCt method. The primer sequences are summarized in Table 1. GAPDH was used as a reference control.

Table 1 The sequences of primers for qRT-PCR

Cell counting kit–8 assay and colony formation assay

Cell proliferation was assessed using the CCK–8 assay. A549 and H1299 cells were digested after 48 h transfection and inoculated in 96–well plates at a density of 2,000 cells/well. After 1, 2, 3, 4, and 5 days of incubation, 10 μl CCK8 solution was added sequentially to each well, and the absorbance (OD) was measured at a wavelength of 450 nm using an enzyme marker after 2 h. The 5–day data were calculated, and cell proliferation curves were plotted. The experiment was repeated three times under the same conditions. A549 and H1299 cells were inoculated in 6–well plates at 500 cells/well and incubated in an incubator at 37 °C in humid air containing 5% CO2 for 14 days until the colonies were visible to the naked eye, followed by washing the colonies twice with PBS and staining with 0.25% crystal violet for 1 h. The colonies were counted and statistically analyzed using ImageJ software.

Cell cycle test

Cells were collected 48 h after transfection and each tube contained 1 × 106 cells. Using the cell cycle assay kit from biosharp, each sample was treated with 500 µl of staining buffer, 25 μl of propidium iodide staining solution (20 ×), and 10 µl of RNase A (50 ×) staining for 30 min. Immediately thereafter, the samples were assayed by flow cytometry.

Wound–healing assay and Transwell invasion assay

The si–transfected (siNC, siCENPO) A549 and H1299 cells were inoculated in 6–well plates at a density of 1 × 106 cells/well. The cells were placed under a microscope for image acquisition at 0 h, 24 h, and 48 h. First, the distance (μm) between scratches at different time points was quantified using Image J software (National Institutes of Health) to calculate the healing area (initial scratch width– final scratch width)/initial scratch width × 100%). Next, A549 and H1299 cells were placed in a Transwell chamber (24–well, 8 mm well) (Corning) at a density of 3 × 104 cells/well and incubated for 24 h at 37 °C in a 5% CO2 incubator. The inner chamber contains 100 μl of cell suspension, and the outer chamber was 600 μl of 1640 medium containing 20% FBS. Afterward, the non–invasive cells in the upper chamber were removed Next, the cells attached to the polycarbonate membrane were fixed with 4% pre–cooled paraformaldehyde for 30 min and stained with 0.1% crystal violet for 20 min at room temperature. Afterward, the cells were placed under a 100 × microscope, and images were captured from five randomly selected fields of view.

Statistical analysis

Data were analyzed and plotted using SPSS 25.0 (IBM, Armonk, NY, USA) and GraphPad Prism 8.0 software (GraphPad Software Inc., San Diego, CA, USA), and analyzed using R language. The results represent at least three replicate experiments and are expressed as Mean \(\pm\) SD. Differences between two groups were compared using an unpaired t–test, and data from more than three groups were analyzed using one–way ANOVA. The Wilcoxon or Kruskal–Wallis test was employed to contrast the data from different groups. P < 0.05 was considered statistically significant (*P < 0.05; **P < 0.01; ***P < 0.001; ****P < 0.0001).

Results

CENPO is upregulated in multiple cancers

We analyzed CENPO expression in pan–cancer using the UCSC Xena database GDC–TCGA data. The results obtained indicate that CENPO is upregulated in most cancers (Fig. 1A), including urothelial carcinoma of the bladder (BLCA), breast cancer invasive carcinoma (BRCA), cervical cancer (CESC), bile duct cancer (CHOL), colon adenocarcinoma (COAD), esophageal cancer (ESCA), glioblastoma multiforme (GBM), head and neck squamous cell carcinoma (HNSC), hepatocellular carcinoma (LIHC), lung adenocarcinoma (LUAD), lung squamous cell carcinoma (LUSC), prostate adenocarcinoma (PRAD), rectal adenocarcinoma (READ), sarcomatoid lung cancer (SARC), gastric adenocarcinoma (STAD), and uterine somatic endometrial cancer (UCEC). However, CENPO was downregulated in renal cell carcinoma (KICH), and thyroid cancer (THCA) (Fig. 1A). In addition, CENPO expression is correlated with the staging of BRCA, KIRC, LIHC, LUAD, LUSC, and THCA (Fig. 1B). Based on the median expression of CENPO in 517 LUAD patients divided into high and low expression groups, differential analysis was performed to obtain 960 differential genes (DEGs) (Fig. 1C). Analysis of differential genes in 58 pairs of LUAD in paired samples with complete data of cancer and paracancerous tissue produced a volcano plot showing 4,220 DEGs (Fig. 1D). From the analysis of CENPO expression in 58 pairs of paired samples with complete data of cancer versus paraneoplastic tissues in LUAD, it could be concluded that the expression of CENPO was significantly higher in tumor tissues of LUAD patients than in normal tissues (Fig. 1E).

Fig. 1
figure 1

CENPO expression profile in pan–cancer. A CENPO expression in TIMER 2.0 database. B Correlation between CENPO expression and pan–cancer stage. C 960 DEGs were associated between LUAD and CENPO. D 58 pairs of paired LUAD tissue DEGs. E 58 pairs of paired LUAD tissues with CENPO differential expression

Diagnostic and prognostic value of CENPO in pan–cancer

We analyzed the diagnostic value of CENPO in various cancers using ROC curves. ROC diagnostic curves were drawn based on the median CENPO expression of 33 cancers, which were divided into high and low expression groups. Figure 2(A–P) shows that CENPO had high diagnostic value in STAD, SARC, READ, LUSA, UCEC, LUAD, GBM, HNSC, LIHC, KICH, BRCA, COAD, ESCA, and BLCA (AUC: 0.8–1), and AUC = 0.638 in THCA, AUC = 0.704 in PRAD. These results suggest that CENPO has good diagnostic value for 14 cancer types and acceptable diagnostic value for THCA and PRAD. KM survival curves showed that increased expression of CENPO was linked to poor OS in UVM (P = 0.0064), KIRC (P = 0.046), SKCM (P = 0.00041), SARC (P = 0.0053), LIHC (P = 0. 0015), LUAD (P = 0.0088), MESO (P = 5.8e–05), ACC (P = 0.0012), KICH (P = 0.0011), LGG (P = 8.4e–06), and PCPG (P = 0.026) (Fig. 3A–L). In contrast, increased expression of CENPO was associated with better OS for CHYM (P = 0.045) and READ (P = 0.046; Fig. 3H, M). DFS results showed that increased expression of CENPO was associated with poor prognosis in KICH (P = 0.03), LUAD (P = 0.014), LIHC (P = 0.0081), ACC (P = 0.00033), SARC (P = 0. 0093), SKCM (P = 0.042), UVM (P = 0.011), PAAD (P = 0.0079), BLCA (P = 0.017), and LGG (P = 0.02; Fig. 4A–J). In contrast, increased expression of CENPO was associated with better DFS in PEAD (P = 0.0093) (Fig. 4K). We downloaded 33 pan-cancer data from the TCGA database and performed CENPO differential expression analysis and survival analysis, taking the intersection of cancers with differential CENPO expression and cancers with significantly associated OS to obtain four cancers, KICH, LIHC, READ, and SARC. The clinicopathological characteristics of these four cancers and the expression of CENPO were then subjected to COX regression analysis, which showed that CENPO and Stage were independent risk factors for KICH. CENPO, age, and M were independent risk factors for LIHC. CENPO, age, and N stage were independent risk factors for READ. CENPO and age were independent risk factors for SARC (Supplementary Fig. 1).

Fig. 2
figure 2

Diagnostic value of CENPO in pan–cancer analysis. AP ROC curves of CENPL in THCA, STAD, SARC, READ, PRAD, LUSA, UCEC, LUAD, GBM, HNSC, LIHC, KICH, BRCA, COAD, ESCA, and BLCA

Fig. 3
figure 3

Prognostic value of CENPO in pan–cancer analysis (AM) Kaplan–Meier analysis of the association between CENPO expression and OS. TPM: Transcripts Per Kilobase of exon model per Million mapped reads

Fig. 4
figure 4

Prognostic value of CENPO in pan–cancer analysis. AKaplan–Meier analysis of the association between CENPO expression and DFS

Methylation analysis and tumor mutation load

Methylation is a novel programmed inflammatory cell death type closely associated with tumor development [20]. To further investigate the mechanisms underlying the aberrant expression of CENPO in various cancers, we analyzed the relationship between CENPO expression and genomic methylation. The results showed that CENPO was positively associated with methylation in UCEC, THYM, THCA, and BRCA (Fig. 5A–D). Conversely, CENPO was negatively correlated with methylation in PCPG, PRAD, SKCM, TGCT, LUAD, LIHC, CESC, KICH, and DLBC (Fig. 5E–M). These findings suggest that methylation changes may affect the expression of CENPO in pan–cancer. Tumor mutational load (TMB) and microsatellite instability (MSI) influence tumorigenesis and progression. Therefore, the correlation between CENPO expression and TMB or MSI in 33 common cancer types was analyzed. CENPO levels were significantly correlated with TMB in ACC, UCES, THYM, THCA, STAD, SKCM, SARC, READ, PRAD, PAAD, LUSC, LUAD, LGG, COAD, CHOL, BRCA and BLCA (Fig. 5N). On the other hand, UVM, UCEC, SARC, LUSC, LUAD, LIHC, KIRP, ESCA, DLBC, COAD, CESC, BLCA, and CENPO were significantly correlated with MSI (Fig. 5O).

Fig. 5
figure 5

Methylation analysis and tumor mutation load. Correlation analysis of CENPO with methylation in pan–cancer. AD CENPO expression is positively correlated with methylation. EM CENPO expression is negatively correlated with methylation. Correlation of CENPO expression with TMB (N). Correlation of CENPO expression with MSI (O)

Analysis of the clinical correlation of CENPO in LUAD

The correlation between CENPO expression and clinical characteristics of 517 LUAD patients in the TCGA database was analyzed. CENPO expression was correlated with gender, age, T stage, N stage, TNM stage, and whether or not radiotherapy was administered, and not with recurrence or distant metastasis (Table 2). These results suggest that CENPO expression may be related to the malignant pathological progression of LUAD. Meanwhile, to further evaluate the prognostic value of CENPO in LUAD, univariate and multifactorial COX regression, and ROC curve analyses were performed. COX regression analysis showed that CENPO expression was associated with an increased risk of death in patients with LUAD (P = 0.021, HR:1.237, 95% CI 1.033–1.481) (Fig. 6A), and multifactorial COX regression analysis showed that TNM stage and the presence or absence of radiation therapy were independent risk factors (P < 0.05) (Fig. 6B). In addition, to facilitate the application of CENPO in clinical evaluation, we constructed column line graphs (Fig. 6C) based on the expression and pathological stage of CENPO. Finally, we evaluated the accuracy of the prognostic assessment model for LUAD patients after 1, 3, and 5 years by calibration curves (Fig. 6D). The column line graph had relatively good accuracy.

Table 2 Correlation of clinicopathological features with CENPO expression in LUAD
Fig. 6
figure 6

Prognostic analysis of CENPO in LUAD (AB) The independent prognostic significance of CENPO in LUAD was analyzed by univariate and multivariate COX regression. C Columnar line graphs based on CENPO expression and pathological staging. D Corrected analysis plot of column line graphs

Functional enrichment and tumor microenvironment

To further explore the function of CENPO in LUAD, we performed GO and KEGG functional enrichment analysis, as shown in Supplementary Table 2, yielding 970 GO items (766 BP items, 90 CC items, 114 MF items) and 35 KEGG pathways. The top 10 items under each classification are included in the bubble plots (Fig. 7A, B). GO analysis showed that these molecules were mainly involved in organelle fission, nuclear division, chromosome segregation, mitotic nuclear division, nuclear chromosome segregation, etc. Cellular component enrichment showed that these molecules were mainly enriched in microtubule, chromosomal region, spindle, and collagen–containing extracellular matrix, etc. Molecular functional enrichment analysis showed that these molecules were mainly associated with tubulin binding, microtubule binding, ATP hydrolysis activity, cytoskeletal motor activity, and microtubule motor activity. KEGG analysis showed that CENPO was associated with Cell cycle, Oocyte meiosis, protein digestion and absorption, Pancreatic secretion, Bile secretion, Progesterone–mediated oocyte maturation, Complement and coagulation cascades, Drug metabolism–cytochrome P450, ECM–receptor interaction, Salivary secretion, p53 signaling pathway, Homologous recombination, Fanconi anemia and other signaling pathways. Tumor microenvironment score (TME score) showed that StromalScore, ImmuneScore, and ESTIMATEScore were higher in the CENPO low expression group than in the CENPO high expression group (Fig. 7C). The results of correlation analysis between CENPO expression and immune cells showed that CENPO in LUAD was significantly and positively correlated with memory–activated CD4 T cells, Macrophages M0, Macrophages M1, T cells CD8, Neutrophils, B cells naive, and resting NK cells. Significant negative correlations were found with Mast cells resting, Dendritic cells resting, B cells memory, Monocytes, and memory resting CD4 T cells (Fig. 7D). In addition, we evaluated the correlation between CENPO expression and immune checkpoints, and showed that CENPO was significantly positively correlated with immune checkpoints such as CD276, CD274, TNFSF4, TNFRSF9, and LAG3. Meanwhile, CENPO was negatively correlated with immune checkpoints such as LGALS9, TNFRSF14, CD40LG, and HHLA2 (Fig. 7E).

Fig. 7
figure 7

Functional enrichment and immune infiltration analysis. AB GO and KEGG functional enrichment analysis of the molecules interacted with CENPO. C Correlation analysis of CENPO high and low expression groups with TME scores. D Correlation of CENPO expression with immune cell infiltration. E Correlation analysis of CENPO expression with immune checkpoint inhibitors

Drug sensitivity analysis

One of the challenges of precision medicine is identifying biomarkers that predict cancer cells' response to drug therapy. Large–scale cancer pharmacogenomic studies offer the possibility to develop new markers [21]. Our screening results showed that the IC50 of 114 drugs correlated with the expression of CENPO (P < 0.001), 15 of which were clinically used lung cancer therapeutics (Fig. 8A–O). These data suggest that CENPO can be used as a biomarker for predicting LUAD drug therapy.

Fig. 8
figure 8

Drug sensitivity analysis. AO CENPO expression in LUAD with drug sensitivity analysis (IC50)

CENPO knockdown inhibits LUAD cell proliferation

To investigate the cellular function of CENPO in LUAD, we examined the mRNA expression levels of CENPO in lung normal epithelial cells BEAS–2b and five LUAD cell lines XWLC–05, H838, A549, SPCA–1, and H1299 by qRT–PCR. CENPO was significantly upregulated in the LUAD cell lines (Fig. 9A). A549 and H1299 cell lines were selected for CENPO knockdown for functional studies. Three siRNAs targeting CENPO were used for the knockdown, and the most efficient siRNA (CENPO–1) was selected for the next experiments (Fig. 9B, C). We investigated the effect of CENPO knockdown on cell proliferation by CCK–8 and clone formation assays. We showed that the downregulation of CENPO significantly reduced the proliferation viability of A549 and H1299 cells (Fig. 9E, F) and the number of colonies formed by the cells (Fig. 9D).

Fig. 9
figure 9

Validation of CENPO expression in LUAD. A CENPO mRNA expression in LUAD cell lines and normal lung epithelial cells. BC Knockdown efficiency of CENPO mRNA in A549 and H1299 cells. D Clone formation assay. EF CCK–8 assay to detect the proliferation of A549 and H1299 cells

Downregulation of CENPO inhibits the migration and invasion of LUAD cells

After cell scratching experiments, we found that the healing ability of A549 and H1299 cells was significantly reduced after the knockdown of CENPO (Fig. 10A, B). Furthermore, the Transwell invasion assay was also performed, and the results showed that the number of H1299 and A549 cells invaded was significantly reduced after the knockdown of CENPO (Fig. 10C). Based on the above results, we inferred that the knockdown of CENPO could inhibit migration and invasion of LUAD cell lines A549 and H1299. Furthermore, WB results showed that the knockdown of CENPO resulted in decreased expression of CDK2 and CyclinB1 in A549 and H1299 cells (Fig. 10D). Flow cytometry results showed that H1299 exhibited a consistent trend of change with A549 cells, although there was no statistically significant difference. It indicates that the knockdown of CENPO significantly induced G0/G1 phase arrest in A549 and H1299 cells (Fig. 10E–F).

Fig. 10
figure 10

Effect of CENPO knockdown on the invasion and migration of A549 and H1299 cells. AC The effect of CENPO knockdown on migration and invasion of A549 and H1299 cells was assessed using a wound–healing assay and transwell assay. D Western blot detected the expression of CDK2 and CyclinB1 in A549 and H1299 cells. ECell cycle distribution was analyzed by Propidium iodide staining 48 h after A549 and H1299 cells were transfected with Si-CENPA or CENPA plasmids.The results represent experiments that were repeated at least three times. The results represent experiments that were repeated at least three times. Data are expressed as mean ± standard deviation

Discussion

Mitotic granules are an important component of chromosome segregation during meiosis and mitosis in normal cells. Their abnormal function can lead to mitotic arrest and altered chromosome numbers [22]. It has been reported that Centromere proteins control chromosome segregation and that abnormalities in Centromere proteins may lead to chromosomal instability in tumors, which in turn lead to tumor progression and drug resistance [23]. CENPO is associated with the development and progress of various tumors, including LUAD. These data suggest the potential of CENPO as a prognostic biomarker.

In this study, we verified that CENPO was significantly highly expressed in most cancer tissues and correlated with the TNM stage by the TCGA database. ROC curve analysis indicated that CENPO might serve as a pan–cancer diagnostic biomarker. Meanwhile, KM survival analysis showed that increased expression of CENPO was associated with poorer OS of UVM, KIRC, SKCM, SARC, LIHC, LUAD, MESO, ACC, KICH, LGG, and PCPG, whereas the increased expression of CENPO was associated with KICH, LUAD, LUSC, SARC, SKCM, UVM, ACC, BLCA LGG, LIHC, and PAAD poorer DFS. These findings suggest that CENPO may play a significant role in tumorigenesis and prognosis, and that abnormal CENPO may be a potential prognostic predictive molecule for cancer. Further analysis revealed that CENPO expression in LUAD was significantly correlated with gender, TNM stage, N stage, T stage, and whether radiotherapy was administered. Columnar plots were constructed to further evaluate the role of CENPO in LUAD prognosis.

Epigenetic mechanisms regulate the expression of immunosynaptic proteins by tumor cells, thereby suppressing immune function [24]. DNA methylation is a key epigenetic mechanism of immune regulation, and cancer cells often use epigenetic dysregulation to silence tumor suppressors or activate oncogenes [25]. Our results showed that CENPO expression was negatively correlated with methylation in LUAD, suggesting that CENPO may be involved in epigenetic modifications to regulate mitophagy assembly and mitophagy function [26]. The tumor microenvironment is a complex survival environment for tumor cells, mainly characterized by immune cells, mesenchymal cells, blood vessels, and an extracellular matrix [27]. In recent years, there has been increasing evidence that the tumor microenvironment and the level of immune cell infiltration are closely related to the growth and progression of tumor cells [28,29,30,31]. Although previous studies have demonstrated the prognostic importance of CENPO in various types of cancer, there are few reports on whether aberrant CENPO expression influences immune cell infiltration. Immune cells of the tumor microenvironment are associated with the efficacy of immunotherapy, including CD4 + T cells [32] CD8 + T cells [33], NK cells [34], and macrophage cells [35]. Our results show for the first time that CENPO expression in LUAD is significant positive associated with the infiltration of memory–activated CD4 T cells, M0 macrophages, M1 macrophages, T cells CD8, Neutrophils, naive B cells, and NK cell resting; our results also show a significant negative correlation with the resting Mast cells, resting Dendritic cells, memory B cells, Monocytes, and memory resting CD4 T cells. Macrophages have been reported to promote tumor progression through tissue remodeling and angiogenesis [36]. High infiltration of M1 macrophages significantly reduced survival time in prostate cancer patients [29]. The M1 macrophage content was significantly higher in the metastatic group than in the non–metastatic group of medulloblastoma [37]. Our results also showed that increased CENPO expression was associated with poor LUAD survival.

The level of immune cell infiltration in TME is critical for tumor prognosis [38]. An immune–based assay called “Immunoscore” was defined to quantify T–cell infiltration in situ and proved superior to the TNM classification of colorectal cancer patients [39]. Furthermore, ImmuneScore has been reported to predict that survival in CRC patients is superior to microsatellite instability [40]. Our results showed that StromalScore, ImmuneScore, and ESTIMATEScore were significantly lower in the CENPO high expression group, which may indicate that increased CENPO expression is associated with disturbed immune status in LUAD patients. Immune checkpoints are involved in physiological immune responses, and immune checkpoints and their ligands are upregulated in the tumor microenvironment of various malignancies to induce antitumor immune responses [41]. The clinical application of immune checkpoint inhibitors has improved the prognosis of many patients with advanced tumors [42]; we showed that CENPO was positively correlated with immune checkpoint CD276, TNFSF4, TNFSF9, CD247, and LAG3. Furthermore, it was negatively correlated with immune checkpoints LGALS9, CD44, TNFRSF14, CD40LG, and HHLA2. These results suggest that CENPO plays an important immunomodulatory role and may serve as a potential immunotherapeutic relevant biomarker. Previous studies reported that the CENP family is associated with immune function [6, 43], and our bioinformatics results also demonstrated that CENPO is associated with immune regulation. We then demonstrated by in vitro experiments that CENPO was highly expressed in A549 and H1299 cell lines and that the knockdown of CENPO significantly inhibited the proliferative activity, migration, and invasive ability of A549 and H1299 cells, suggesting that CENPO may affect LUAD metastasis. Studies have reported that the knockdown of CENPO inhibited GC cell growth, induced apoptosis, and decreased the expression of ATM, cyclin D1, and c–Jun [12]. Knockdown of CENPO in CRC downregulated the expression of N–cadherin, Vimentin, Snail, and CCND1, inhibited cell proliferation, and attenuated migration and invasion [15]. This is consistent with our Bioinformatics Analysis, the functional enrichment results showed significant enrichment of cell cycle, P53 pathway, and extracellular matrix. Thus, our study suggests that the knockdown of CENPO significantly decrease the expression of CDK2 and CyclinB1, and induce G0/G1 arrest in A549 and H1299 cells. We speculate that CENPO may promote LUAD progression by regulating the cell cycle. These results suggest that CENPO may be a potential molecular target for anti–tumor therapy in LUAD, and the specific mechanism of CENPO in regulating the development of LUAD will be further investigated in the future.

Availability of data and materials

All data generated or analyzed during this study are included in this published article (and its Supplementary Information files).

The datasets generated during and/or analyzed during the current study are available in the UCSC Xena(https://xenabrowser.net/datepages/) and GEPIA repository GEPIA (Gene Expression Profiling Interactive Analysis) (cancer-pku.cn).

Abbreviations

CENPO:

Centromere protein O

CRC:

Colorectal Cancer

DEG:

Differentially expressed genes

DFS:

Disease–Free Survival

GEO:

Gene Expression Omnibus

GO:

Gene Ontology

GSEA:

Gene Set Enrichment Analysis

KEGG:

Kyoto Encyclopedia of Genes and Genomes

LUAD:

Lung adenocarcinoma

MSI:

Microsatellite instability

OS:

Overall survival

TNM:

Tumor node metastasis classification

TME:

Tumor microenvironment

TMB:

Tumor mutation burden

TCGA:

The Cancer Genome Atlas

References

  1. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, Bray F. Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J Clin. 2021;71(3):209–49.

    Article  PubMed  Google Scholar 

  2. Dyba T, Randi G, Bray F, Martos C, Giusti F, Nicholson N, Gavin A, Flego M, Neamtiu L, Dimitrova N, et al. The European cancer burden in 2020: Incidence and mortality estimates for 40 countries and 25 major cancers. Eur J Cancer. 2021;157:308–47.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Durante M, Reppingen N, Held K. Immunologically augmented cancer treatment using modern radiotherapy. Trends Mol Med. 2013;19(9):565–82.

    Article  PubMed  CAS  Google Scholar 

  4. Hashim D, Boffetta P, La Vecchia C, Rota M, Bertuccio P, Malvezzi M, Negri E. The global decrease in cancer mortality: trends and disparities. Ann Oncol. 2016;27(5):926–33.

    Article  PubMed  CAS  Google Scholar 

  5. Li W, Liu J. The Prognostic and Immunotherapeutic Significance of AHSA1 in Pan-Cancer, and Its Relationship With the Proliferation and Metastasis of Hepatocellular Carcinoma. Front Immunol. 2022;13:845585.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  6. Feng Z, Chen Y, Cai C, Tan J, Liu P, Chen Y, Shen H, Zeng S, Han Y. Pan-Cancer and Single-Cell Analysis Reveals CENPL as a Cancer Prognosis and Immune Infiltration-Related Biomarker. Front Immunol. 2022;13:916594.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  7. Pan-cancer analysis of whole genomes. Nature. 2020;578(7793):82–93.

    Article  Google Scholar 

  8. Weinstein J, Collisson E, Mills G, Shaw K, Ozenberger B, Ellrott K, Shmulevich I, Sander C, Stuart J. The Cancer Genome Atlas Pan-Cancer analysis project. Nat Genet. 2013;45(10):1113–20.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Sharma A, Dimitrov S, Hamiche A, Van Dyck E. Centromeric and ectopic assembly of CENP-A chromatin in health and cancer: old marks and new tracks. Nucleic Acids Res. 2019;47(3):1051–69.

    Article  PubMed  CAS  Google Scholar 

  10. Ding N, Li R, Shi W, He C. CENPI is overexpressed in colorectal cancer and regulates cell migration and invasion. Gene. 2018;674:80–6.

    Article  PubMed  CAS  Google Scholar 

  11. Liu Y, Yao Y, Liao B, Zhang H, Yang Z, Xia P, Jiang X, Ma W, Wu X, Mei C, et al. A positive feedback loop of CENPU/E2F6/E2F1 facilitates proliferation and metastasis via ubiquitination of E2F6 in hepatocellular carcinoma. Int J Biol Sci. 2022;18(10):4071–87.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  12. Cao Y, Xiong J, Li Z, Zhang G, Tu Y, Wang L, Jie Z. CENPO expression regulates gastric cancer cell proliferation and is associated with poor patient prognosis. Mol Med Rep. 2019;20(4):3661–70.

    PubMed  PubMed Central  CAS  Google Scholar 

  13. Hori T, Okada M, Maenaka K, Fukagawa T. CENP-O class proteins form a stable complex and are required for proper kinetochore function. Mol Biol Cell. 2008;19(3):843–54.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  14. Zhang S, Xie Y, Tian T, Yang Q, Zhou Y, Qiu J, Xu L, Wen N, Lv Q, Du Z. High expression levels of centromere protein A plus upregulation of the phosphatidylinositol 3-kinase/Akt/mammalian target of rapamycin signaling pathway affect chemotherapy response and prognosis in patients with breast cancer. Oncol Lett. 2021;21(5):410.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  15. Liu Z, Chen C, Yan M, Zeng X, Zhang Y, Lai D. CENPO regulated proliferation and apoptosis of colorectal cancer in a p53-dependent manner. Discov Oncol. 2022;13(1):8.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Liu Y, Xiong S, Liu S, Chen J, Yang H, Liu G, Li G. Analysis of Gene Expression in Bladder Cancer: Possible Involvement of Mitosis and Complement and Coagulation Cascades Signaling Pathway. J Comput Biol. 2020;27(6):987–98.

    Article  PubMed  CAS  Google Scholar 

  17. Ritchie M, Phipson B, Wu D, Hu Y, Law C, Shi W, Smyth G. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Wu T, Hu E, Xu S, Chen M, Guo P, Dai Z, Feng T, Zhou L, Tang W, Zhan L, et al. clusterProfiler 4.0: A universal enrichment tool for interpreting omics data. Innovation (Camb). 2021;2(3):100141.

    PubMed  CAS  Google Scholar 

  19. Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  20. Park C, Jeong K, Park J, Jung S, Bae J, Kim K, Ock C, Kim M, Keam B, Kim T, et al. Pan-cancer methylation analysis reveals an inverse correlation of tumor immunogenicity with methylation aberrancy. Cancer Immunol Immunother. 2021;70(6):1605–17.

    Article  PubMed  CAS  Google Scholar 

  21. Ba-Alawi W, Nair S, Li B, Mammoliti A, Smirnov P, Mer A, Penn L, Haibe-Kains B. Bimodal Gene Expression in Patients with Cancer Provides Interpretable Biomarkers for Drug Sensitivity. Can Res. 2022;82(13):2378–87.

    Article  CAS  Google Scholar 

  22. Liu X, Wang H, Zhao G. Centromere Protein A Goes Far Beyond the Centromere in Cancers. Mol Cancer Res. 2022;20(1):3–10.

    Article  PubMed  CAS  Google Scholar 

  23. Slee R, Steiner C, Herbert B, Vance G, Hickey R, Schwarz T, Christan S, Radovich M, Schneider B, Schindelhauer D, et al. Cancer-associated alteration of pericentromeric heterochromatin may contribute to chromosome instability. Oncogene. 2012;31(27):3244–53.

    Article  PubMed  CAS  Google Scholar 

  24. Berglund A, Mills M, Putney R, Hamaidi I, Mulé J, Kim S. Methylation of immune synapse genes modulates tumor immunogenicity. J Clin Investig. 2020;130(2):974–80.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  25. Chen C, Yin N, Yin B, Lu Q. DNA methylation in thoracic neoplasms. Cancer Lett. 2011;301(1):7–16.

    Article  PubMed  CAS  Google Scholar 

  26. Samel A, Cuomo A, Bonaldi T, Ehrenhofer-Murray A. Methylation of CenH3 arginine 37 regulates kinetochore integrity and chromosome segregation. Proc Natl Acad Sci USA. 2012;109(23):9029–34.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  27. Anderson N, Simon M. The tumor microenvironment. Current biology : CB. 2020;30(16):R921–5.

    Article  PubMed  CAS  Google Scholar 

  28. Zhang S, Zhang E, Long J, Hu Z, Peng J, Liu L, Tang F, Li L, Ouyang Y, Zeng Z. Immune infiltration in renal cell carcinoma. Cancer Sci. 2019;110(5):1564–72.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  29. Wu Z, Chen H, Luo W, Zhang H, Li G, Zeng F, Deng F. The Landscape of Immune Cells Infiltrating in Prostate Cancer. Front Oncol. 2020;10:517637.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Xu C, Sui S, Shang Y, Yu Z, Han J, Zhang G, Ntim M, Hu M, Gong P, Chen H, et al. The landscape of immune cell infiltration and its clinical implications of pancreatic ductal adenocarcinoma. J Adv Res. 2020;24:139–48.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  31. Brassart-Pasco S, Brézillon S, Brassart B, Ramont L, Oudart J, Monboisse J. Tumor Microenvironment: Extracellular Matrix Alterations Influence Tumor Progression. Front Oncol. 2020;10:397.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Borst J, Ahrends T, Bąbała N, Melief C, Kastenmüller W. CD4 T cell help in cancer immunology and immunotherapy. Nat Rev Immunol. 2018;18(10):635–47.

    Article  PubMed  CAS  Google Scholar 

  33. Raskov H, Orhan A, Christensen J, Gögenur I. Cytotoxic CD8 T cells in cancer and cancer immunotherapy. Br J Cancer. 2021;124(2):359–67.

    Article  PubMed  CAS  Google Scholar 

  34. Bi J, Tian Z. NK Cell Dysfunction and Checkpoint Immunotherapy. Front Immunol. 1999;2019:10.

    Google Scholar 

  35. Tan Y, Wang M, Zhang Y, Ge S, Zhong F, Xia G, Sun C. Tumor-Associated Macrophages: A Potential Target for Cancer Therapy. Front Oncol. 2021;11:693517.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  36. Schmieder A, Michel J, Schönhaar K, Goerdt S, Schledzewski K. Differentiation and gene expression profile of tumor-associated macrophages. Semin Cancer Biol. 2012;22(4):289–97.

    Article  PubMed  CAS  Google Scholar 

  37. Zhang J, Yuan X, Wang Y, Liu J, Li Z, Li S, Liu Y, Gong X, Sun Y, Wu W, et al. Tumor-Associated Macrophages Correlate With Prognosis in Medulloblastoma. Front Oncol. 2022;12:893132.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  38. Xu H, Xu Q, Yin L. Prognostic value of tumor immune cell infiltration patterns in colon adenocarcinoma based on systematic bioinformatics analysis. Cancer Cell Int. 2021;21(1):344.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  39. Bruni D, Angell H, Galon J. The immune contexture and Immunoscore in cancer prognosis and therapeutic efficacy. Nat Rev Cancer. 2020;20(11):662–80.

    Article  PubMed  CAS  Google Scholar 

  40. Mlecnik B, Bindea G, Angell H, Maby P, Angelova M, Tougeron D, Church S, Lafontaine L, Fischer M, Fredriksen T, et al. Integrative Analyses of Colorectal Cancer Show Immunoscore Is a Stronger Predictor of Patient Survival Than Microsatellite Instability. Immunity. 2016;44(3):698–711.

    Article  PubMed  CAS  Google Scholar 

  41. Toor S, Sasidharan Nair V, Decock J, Elkord E. Immune checkpoints in the tumor microenvironment. Semin Cancer Biol. 2020;65:1–12.

    Article  PubMed  CAS  Google Scholar 

  42. Kocikowski M, Dziubek K, Parys M: Hyperprogression Under Immune Checkpoint-Based Immunotherapy-Current Understanding, The Role of PD-1/PD-L1 Tumour-Intrinsic Signalling, Future Directions and a Potential Large Animal Model. Cancers 2020, 12(4).

  43. Shao C, Wang Y, Duan H, Ding P, Zhang Y, Ning J, Han J, Jiang T, Yan X. CENP-UMitosis-related gene as a potential biomarker in malignancy. Ann Transl Med. 2021;9(23):1744.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

Download references

Acknowledgements

We acknowledge UCSC Xena and GEPIA for the data sources.

Funding

This work was supported by the National Natural Science Foundation of China (No. 81560380).

Author information

Authors and Affiliations

Authors

Contributions

YC conceived and designed the study. YBW performed the experiments and drafted the manuscript. DWY completed the data analysis. YL, FHL, WBS, and HL supervised the study and prepared figures. LT, ZLF, YW, FL, and YL revised the manuscript and prepared tables.

Corresponding author

Correspondence to Yan Chen.

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.

Supplementary Information

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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wang, Y., Ye, D., Li, Y. et al. Prognostic and immune infiltrative biomarkers of CENPO in pan–cancer and its relationship with lung adenocarcinoma cell proliferation and metastasis. BMC Cancer 23, 735 (2023). https://doi.org/10.1186/s12885-023-11233-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12885-023-11233-2