Identification of COPA as a potential prognostic biomarker and pharmacological intervention target of cervical cancer by quantitative proteomics and experimental verification
Journal of Translational Medicine volume 20, Article number: 18 (2022)
Cervical cancer is the most fatal gynecological carcinoma in the world. It is urgent to explore novel prognostic biomarkers and intervention targets for cervical cancer.
Through integrated quantitative proteomic strategy, we investigated the protein expression profiles of cervical cancer; 28 fresh frozen tissue samples (11 adenocarcinoma (AC), 12 squamous cell carcinoma (SCC) and 5 normal cervixes (HC)) were included in discover cohort; 45 fresh frozen tissue samples (19 AC, 18 SCC and 8 HC) were included in verification cohort; 140 paraffin-embedded tissues samples of cervical cancer (85 AC and 55 SCC) were used for immunohistochemical evaluation (IHC) of coatomer protein subunit alpha (COPA) as a prognostic biomarker for cervical cancer; how deficiency of COPA affects cell viability and tumorigenic ability of cervical cancer cells (SiHa cells and HeLa cells) were evaluated by cell counting kit-8 and clone formation in vitro.
We identified COPA is a potential prognostic biomarker for cervical cancer in quantitative proteomics analysis. By retrospective IHC analysis, we additionally verified the proteomics results and demonstrated moderate or strong IHC staining for COPA is an unfavourable independent prognostic factor for cervical cancer. We also identified COPA is a potential pharmacological intervention target of cervical cancer by a series of in vitro experiments.
This study is the first to demonstrate that COPA may contribute to progression of cervical cancer. It can serve as a potential prognostic biomarker and promising intervention target for cervical cancer.
Cervical cancer is the most fatal gynecological carcinoma worldwide . Approximately 13 million women have been threatened by this malignant tumor. Early diagnosis and accurate prediction of prognosis are essential for guiding the clinical treatment . ThinPrep cytology test (TCT)  and high-risk human papillomavirus (HPV) DNA test  were used for early screening . Astonishingly, however, incidence rate of cervical cancer is rising steadily [6,7,8]. Unfortunately, there has been no significant improvement in treatment strategies for cervical cancer in the last decades. This brings forward the urgent need for further exploring prognostic biomarkers and innovative therapies to improve the prognosis of patients with cervical cancer.
As cancer treatment continues to move toward a focus on precision medicine, omics diagnostic methods are advantageous in early diagnosing and monitoring the progression of the disease. Exploration of biological characterization will facilitate the choice of targeted therapies for cervical cancer. Proteomics has emerged as a powerful tool in identifying potential prognostic biomarkers and pharmacological intervention targets in cancers . Compared to other technologies, Mass spectrometry (MS) has been widely explored in surveillance and precision treatment of various cancer, standing out for its high accuracy, sensitivity, resolution, and throughput. Technological advances have optimized the performance of MS in proteomics and other analysis [10, 11], which opens new possibilities to strengthen surveillance of cancer patients in clinical situations. However, there are limited data regarding cervical cancer [12,13,14,15].
The present study was designed to investigate global protein expression profiles of cervical cancer. We employed an integrated analysis strategy comprising MS-based quantitative proteomics and parallel reaction monitoring (PRM)-based targeted proteomics in fresh frozen tissue samples of cervical cancer and healthy cervix samples. In this study, we focused on coatomer protein subunit alpha COPA, a tumor-promoting gene. COPA is involved in the movement of vesicles within the Golgi and retrograde transport of cargo proteins between endoplasmic reticulum and Golgi . Recently, COPA was reported involves in a novel carcinogenesis of hepatocarcinoma, a virus-related tumor [15, 17]. However, the function of COPA in cervical cancer is unknown. In untargeted (“discovery”) proteomics and targeted (“verification”) proteomics, COPA was consistently identified as a potential prognostic biomarker of cervical cancer. We additionally verified the expression of COPA in a large number of clinical samples of cervical cancer using IHC analysis. Furthermore, we demonstrated moderate (score 2) or strong (score 3) IHC staining for COPA was a unfavorable independent risk factor for cervical cancer. Finally, through a series of in vitro experiments, the notable contribution of COPA to the progression of cervical cancer was verified in two cervical cancer cell lines.
Together, we are the first to report COPA can serve as a potential prognostic biomarker and pharmacological intervention target for cervical cancer.
Materials and methods
Patients and samples
Our study was designed and analyzed according to the flow chart as shown in Fig. 1. This study has been approved by the Biomedical Ethics Committee of Guangdong Provincial People's Hospital (Permission Number: LW-  No.1, KY-Q-2021–077-01) and Zhuhai People's Hospital. All of the patients provided written informed consent for their data to be used for this research. Tissue samples of from patients with cervical cancer were collected from patients who undergone radical hysterectomy and pelvic lymphadenectomy between January 2021 and May 2021 at the Department of Gynecology, Guangdong Provincial People's Hospital and the Department of Gynecology, Zhuhai People's Hospital. Normal cervix samples from patients with leiomyoma undergoing total hysterectomy during the same time period were obtained and used as healthy controls. TCT and HPV DNA tests were negative in all patients in the control group. For untargeted proteomic analysis, 28 fresh frozen tissue samples, including 11 cervical AC and 12 cervical SCC, and 5 normal cervixes were collected. For targeted proteomics analysis, 45 fresh frozen tissue samples, including 19 cervical AC and 18 cervical SCC, and 8 normal cervixes were collected. For retrospective immunohistochemistry (IHC) analysis, 155 archival formalin-fixed and paraffin-embedded (FFPE) samples of cervical cancer were included. We processed the COPA staining data by removing the samples with staining failure due to tissue detachment. Finally, 140 samples were successfully scored for COPA staining. Basic characteristics of included patients with cervical cancer were summarized in Additional file 1: Table S1, which showed no significant differences (p > 0.05) in three cohorts.
Extraction and tryptic digestion of protein
Tissue samples were treated by ground into powder with liquid nitrogen, and then sonicated three times on the ice in lysis buffer (1% protease inhibitor cocktail and 8 M urea). After centrifugation at 12,000×g for 10 min at 4 °C, the supernatant was collected, and then the total protein concentration was determined using a BCA Protein Assay Kit according to manufacturer’s instructions. For digestion, the soluble proteins were treated with dithiothreitol (5 mM) at 56 °C for 30 min, and then they were alkylated using iodoacetamide (11 mM) at room temperature for 15 min in the dark. Then the protein samples were diluted via adding 100 mM TEAB. Eventually, trypsin was added at a trypsin-to-protein mass ratio of 1:50 in the first digestion overnight, and then a trypsin-to-protein mass ratio of 1:100 in a second digestion for 4 h.
Quantitative proteomic analysis
Candidate protein biomarkers of cervical cancer was identified via liquid chromatography-tandem mass spectrometry. The pre-processed peptides were detected using nano spray ion (NSI) source followed by tandem MS/MS in TIMS ToF Pro, which was coupled online to UPLC system. Resulting MS/MS data were then processed applying the MaxQuant version 18.104.22.168. Tandem mass spectra were searched against the UniProt human database , which is a large resource of protein sequences and associated detailed annotation. We specified Trypsin/P as cleavage enzyme which allowed up to two missing cleavages. Mass tolerance of precursor ions was set as 20 ppm and 5 ppm in the first and main search, respectively. For mass tolerance for fragment ions, it was set as 20 ppm. A fixed modification was specified by the carbamidomethyl on Cys, and variable modifications were specified by acetylation modification and oxidation on Met. The false discovery rate (FDR) was set as < 1%, and the minimum score was set as > 40 for modified peptides.
Targeted proteomics analysis
Based on potential biological functions of the identified differentially expressed proteins (DEPs), 20 candidate proteins from the discovery measurements were selected for further PRM analysis. According to the previous method of protein digestion, the protein samples were extracted from 45 fresh frozen tissues. Peptides were detected by tandem MS/MS in Q ExactiveTM Plus with UPLC system. The MS data were processed by Skyline version 3.6.
For the peptide settings, the cleavage enzyme was specified by trypsin [KR/P], and the max missing cleavage was set as two. Peptide length was set as 8–25 and we set carbamidomethyl on Cys and oxidation on Met as the variable modification, with maximum variable modifications as three.
For the transition settings, we set the precursor charges as two and three, ion charges as one and two, and ion types as b, y, and p. We set the product ions from ion three to the last ion, and the ion match tolerance was 0.02 Da.
Retrospective immunohistochemical analysis
Paraffin sections were prepared and IHC staining were performed as described previously . Rabbit monoclonal antibody targeting COPA (ab181224) were obtained from Abcam. The dilution for the primary antibody was 1:100 for COPA. In order to eliminate inter-rater variability and performed blind analysis, the percentage of positive tumor cells and staining intensity of IHC staining were scored by three experienced pathologists from the Department of Pathology in Zhuhai People's Hospital and Department of Pathology in Guangdong Provincial People's Hospital. The results were scored as follows: negative (score 0), weak (score 1), moderate (score 2) and strong (score 3) staining.
Through searching the Human Protein Atlas (HPA) database (http://www.proteinatlas.org/), highly expressed COPA RNA in SiHa cells and HeLa cells was determined (Additional file 1: Fig. S1). Therefore, these two cell lines were used in the subsequent experiments. SiHa cells and HeLa cells were obtained from the American Type Culture Collection (ATCC). SiHa cells were cultured in Minimum Essential Medium (Gibco) supplemented with 10% fetal bovine serum, 100 U/ML penicillin and 100 mg/L streptomycin (Invitrogen). HeLa cells were cultured in Dulbecco's Modified Eagle Medium (Gibco) supplemented with 10% fetal bovine serum, 100 U/ML penicillin and 100 mg/L streptomycin (Invitrogen).
Transfection of small interfering RNA and evaluation biological behaviors of cervical cancer cells
Predesigned small interfering RNA (siRNA) specific for COPA (4390824) and control siRNA (4390846) were obtained from ThermoFisher Scientific. COPA siRNA or control siRNA were introduced into SiHa cells and HeLa cells using the Lipofectamine 3000 reagent (ThermoFisher Scientific) according to the manufacturer’s instructions. To confirm the decrease of COPA protein expression in SiHa cells and HeLa cells with COPA deletion, Western blotting were performed as previously described . The dilutions for the primary antibodies were as follow: 1:1000 for COPA and 1:2000 for β-actin. We observed the impact of COPA on the biological behaviors of SiHa cells and HeLa cells using cell counting kit-8 and colony formation. For determination of cell viability, a cell-counting kit 8 (Dojindo Molecular Technologies, Japan) was used according to the user manual. To detect cell tumorigenicity, clone formation assay was performed as previously described .
Bioinformatics and statistics
To better understand the proteomic variances, DEPs (Differentially expressed proteins) and KEGG analysis between comparison groups were performed using ANOVA test (p < 0.05). We investigated the relationship among the DEPs and visualized with protein–protein interaction (PPI) networks by the STRING webserver (https://string-db.org/)  and Cytoscape software . Kaplan–Meier survival analysis based on TCGA database was executed to explore the relationship between expression level of COPA RNA and prognosis of cervical cancer patients using GEPIA2 software . For IHC analysis, the correlation between the expression level of COPA and the clinicopathological characteristics was estimated using the Mann–Whitney U test for continuous variables and Fisher's exact test or χ2 for categorical variables. Violin plot, heatmaps and stacked column chart were applied to demonstrate the correlation. Kaplan–Meier survival analysis was used to calculate the cumulative probability of overall survival (OS) of patients in IHC cohort, log-rank tests were used to evaluate the differences. A two-step Cox regression analysis was performed using SPSS (version 21.0) to determine the potential prognostic value of COPA in clinical situations. The hazard ratios (HRs) and their 95% confidence intervals (CIs) were also estimated. In vitro experimental phase, data were presented as mean ± standard error. Statistical significances were analyzed with one-way ANOVA using SPSS. A p-value < 0.05 was considered to be statistically significant.
Screening of candidate proteins using untargeted quantitative proteomic
In discovery proteomics analysis, overall 7685 proteins were recognized. Among the identified proteins, 6982, 7178, and 5442 proteins were found in the cervical AC group, cervical SCC group, and healthy control group, respectively (Fig. 2A). We quantified 6602 proteins with one or more unique peptides. Compared with the healthy control group, a total of 1679 quantified proteins were significantly differentially expressed in the cervical AC group, including 1085 upregulated proteins and 594 downregulated proteins. We identified 2083 DEPs between the cervical SCC and healthy control groups, including 1309 upregulated and 774 downregulated proteins. Moreover, we identified 822 DEPs between the cervical AC and cervical SCC groups, including 427 upregulated and 395 downregulated proteins. A two-dimensional principal component analysis was executed to display the abundance variation of proteins within and among the three groups. Based on the log2-ratio of each sample over the average of all samples, we observed that cervical cancer group (cervical AC group and cervical SCC group) was completely separated from the healthy group (Fig. 2B).
We analyzed the DEPs in different comparable groups (AC vs. healthy controls, SCC vs. healthy controls, AC vs. SCC). A total of 85 DEPs overlapped across all the comparable groups. In addition, 369 DEPs were only present in the AC vs. SCC comparable group, 228 DEPs were only displayed in the AC vs. healthy controls comparable group, and 528 DEPs were only present in the SCC vs. healthy controls comparable group (Fig. 2C). The upregulated and downregulated DEPs in the three comparable groups were also analyzed. A total of 22 upregulated DEPs overlapped across all three comparable groups. In addition, 317 upregulated DEPs were only present in the AC vs. SCC comparable group, 163 upregulated DEPs were only displayed in the AC vs. healthy controls comparable group, and 473 upregulated DEPs were only present in the SCC vs. healthy controls comparable group (Fig. 2D). Six downregulated DEPs overlapped across the three comparable groups. In addition, 340 downregulated DEPs were only present in the AC vs. SCC comparable group, 65 downregulated DEPs were only displayed in the AC vs. healthy controls comparable group, and 294 downregulated DEPs were only present in the SCC vs. healthy controls comparable group (Fig. 2E).
Using KEGG database for unbiased clustering and cluster specific enrichment analysis, we found that the enriched pathways in the three groups were different (Fig. 2F). All the pathways were clustered into six clusters. Pathways enriched in the healthy group mainly included cluster5 (complement and coagulation cascades, focal adhesion, proteoglycans in cancer, and phospholipase D signaling pathway) and cluster6 (vascular smooth muscle contraction, cGMP-PKG signaling pathway, gap junction, and apoptosis). We also found that cluster4 (amino sugar and nucleotide sugar metabolism, other glycan degradation, and glycosaminoglycan degradation) was significantly enriched in the AC group, and cluster1 (p53 signaling pathway, proteasome, and ubiquitin mediated proteolysis) was specifically enriched in the SCC group.
Validation of candidate proteins and identification COPA as a potential biomarker using targeted proteomics
In targeted proteomics, 19 proteins were successfully quantified and verified. Among them, seven proteins were significantly up-regulated (p < 0.05) and five proteins were down-regulated (p < 0.05) in cervical SCC versus healthy controls (p < 0.05); Six proteins were significantly up-regulated (p < 0.05) and four proteins were down-regulated (p < 0.05) in cervical AC versus healthy controls; Six verified proteins were significantly up-regulated, while four verified proteins were down-regulated in cervical AC compared to cervical SCC (Table 1). Based on quantification of the following two unique peptides: GITGVDLFGTTDAVVK (Fig. 3A) and CPLSGACYSPEFK (Fig. 3B), COPA was verified as a potential biomarker of cervical cancer via targeted proteomics. COPA was observed to up-regulated significantly in cervical AC group versus healthy control group, in the cervical SCC group versus healthy control group, and in the cervical AC group versus cervical SCC group with a fold change (FC) of 2.8 (p = 0.000253), 2.07 (p = 0.004) and 1.35 (p = 0.0295), respectively (Fig. 3C). Then, further bioinformatics analysis and clinical value analysis of COPA were performed.
Functional analysis of COPA using bioinformatics
To obtain an unbiased overview of COPA function in a biological context, the STRING webserver and Cytoscape software were used to create a PPI network including all DEPs from cervical cancer group (both AC and SCC) vs. healthy controls. A total of 45 DEPs with first neighbour relationships (undirected) to COPA were found (Fig. 3D). The analysis of the KEGG pathway shows that the 46 DEPs were mainly enriched in the following three pathways: SNARE interactions in vesicular transport, legionellosis and endocytosis pathway (Fig. 3E). Based on TCGA database, the KM curve indicated that the prognosis (quartile survival) of patients with high COPA RNA expression were significantly worse (p = 0.042), (Fig. 3F). The hazard ratio (HR) based on the Cox PH model was calculated, high expressed COPA RNA was an independent prognostic factor of cervical cancer with an HR of 2.1 (p = 0.047). Therefore, we can speculate that high expressed COPA protein may also be an unfavorable prognostic factor of patients with cervical cancer. Subsequently, the expression of COPA protein in cervical cancer tissues was additionally verified in IHC cohort.
Correlation analysis between COPA expression and clinicopathological characteristics in IHC cohort
Human Protein Atlas analyses for the COPA (https://www.proteinatlas.org/ENSG00000122218-COPA/cell) predicted that COPA protein would be detected in the nucleoplasm, cytosol, and Golgi apparatus (Additional file 1: Fig. S2). COPA protein was retrospectively evaluated in a larger cohort of patients with cervical cancer using IHC analysis (Fig. 4), 140 patients were successfully scored for COPA staining. There were 137 (97.9%) samples were positively stained for COPA. We further investigated the correlation between COPA staining and clinicopathological characteristics of patients with cervical cancer. The clinicopathological characteristics of patients and the distribution of COPA staining were shown in Table 2. The specific results of the analysis are as follows:
As shown in Fig. 5 and Table 2, the positive rate in cervical AC and cervical SCC was 98.8% (Fig. 5A) and 96.4% (Fig. 5B), respectively. Three (2.1%) sample was negative (score 0) stained for COPA in the cervical cancer sample. Only one (1.2%) tissue sample was negative (score 0) stained for COPA in the cervical AC, and 2 (3.6%) tissue sample was negative (score 0) stained for COPA in the cervical SCC. The expression of COPA in cervical AC was mainly (90.6%) moderate (score 2) and strong (score 3) expression; 8.2% of the patients with cervical AC were weakly positive for COPA. In comparison, moderate (score 2) and strong (score 3) expression of COPA in cervical SCC was less, that is 61.8%; 34.5% of the patients with cervical SCC were weakly positive for COPA.
The frequency distribution of the age of 140 patients with cervical cancer was displayed in Additional file 1: Figure S3. The violin plot indicated the distribution of ages of patients with different expressions of COPA. The lower quartile, median, and upper quartile ages of cervical AC patients with weak (score 1) COPA expression were 38 years, 49 years, and 55 years, respectively. For the cervical AC patients with moderate (score 2) COPA expression, the median age was 47 years, while the lower and upper quartiles were 45 years and 55 years, respectively. Among the 32 cervical AC patients with strong COPA expression (score 3), the lower quartile, median, and upper quartile ages were 46.25 years, 52 years, and 54.75 years, respectively (Fig. 5C). The lower quartile, median, and upper quartile ages of cervical SCC patients with weak (score 1) COPA expression were 39 years, 48 years, and 57 years, respectively. For the cervical SCC patients with moderate (score 2) COPA expression, median age was 56 years, and lower quartile and upper quartile were 44 years and 62 years, respectively. Among the nine cervical SCC patients with strong COPA expression (score 3), the lower quartile, median, and upper quartile ages were 43.5 years, 57 years, and 65 years, respectively (Fig. 5D). We found that the median age of patients with low expression of COPA was similar in the two pathological types of cervical cancer, but the median age of patients with moderate (score 2) and strong (score 3) expression of COPA in cervical AC were relatively younger than that of in cervical SCC.
The FIGO stage of 140 patients ranged from stage IA1 to IIIC1. Most of patients were distributed in stages IB1 to IIA2, with stage IB1 being the most dominant, followed by stage IB2. Negative (score 0) or weak (score 1) expression of COPA in cervical AC patients were all diagnosed as stage IB, including 6 cases of stage IB1 and 2 cases of stage IB2; Even in stage I patients with cervical AC, 86.2% had moderate (score 2) and high expression of COPA. The strong expression of COPA was found in all cervical AC patients diagnosed with stage IIA and above (Fig. 5E). Similarly, cervical SCC patients with negative COPA expression (score 0) were all diagnosed as stage IB2; However, the positive expression of COPA dispersed in cervical SCC were relatively scattered than that in cervical AC; 94.7% patients of cervical SCC with weak COPA expression (score 1) were diagnosed as stage I (including 15.8% for stage IA1, 5.3% for stage IA2, 5.3% for stage IB1, 42.1% for stage IB2, and 26.3% for stage IB3) except one (5.3%) stage IIA1 patient with cervical SCC; patients with moderate COPA expression (score 2) were mainly (56.0%) categorized into stages IB (including 8.0% for stage IB1, 32.0% for stage IB2, and 16.0% for stage IB3), followed by stages IA (20.0%, including 8.0% for stage IA1, 12.0% for stage IA2) and IIA (16.0%, including 4.0% for stage IIA1, 12.0% for stage IIA2); patients of cervical SCC with high COPA expression (score 3) were mainly (88.9%) diagnosed as stage IB (including 22.2% for stage IB1, 33.3% for stage IB2, and 33.3% for stage IB3) except one (11.1%) stage IA1 patient with cervical SCC (Fig. 5F). We found that all tissues of AC (Figs. 5G) and most tissues of SCC (Figs. 5H) with stage II and above tumors were strongly positive stained of COPA (score 3), while a considerable number of tissues from patients with stage I cervical cancer were weakly stained for COPA (score 1). These results demonstrated that COPA protein level was significantly elevated in advanced cervical cancer, which suggested COPA may notable contributions to the progression of cervical cancer. High expressed COPA protein may also be an unfavorable prognostic factor of patients with cervical cancer.
Parametrial invasion (PI) was present in four (2.9%) patients with cervical cancer, the two pathological subtypes accounted for 50% respectively; all cervical AC patients with PI were in the moderate COPA expression (score 2) group (Fig. 6A). PI was also present in two (3.6%) patients with cervical SCC, one patient with weak COPA expression (score 1) and the other with moderate COPA expression (score 2) (Fig. 6B).
Lymph node metastasis
Twenty patients (14.3%) of cervical cancer with lymph node metastasis (LNM) were found by histopathology, most patients (80.0%) were in the moderate COPA expression (score 2) group and the high COPA expression (score 3) group. Eleven (12.9%) patients with cervical AC were presented LNM. Of these, 90.9% patients were in the moderate COPA expression (score 2) group and the high COPA expression (score 3) group, except one (9.1%) patient with weak COPA expression (Fig. 6A). Nine (16.4%) patients with cervical SCC were found LNM. Of these, COPA was moderately (score 2) or highly (score 3) expressed in 6 (66.7%) cases, weakly expressed in 2 (22.2%) cases, and negatively expressed in one (11.1%) case (Fig. 6B).
Patients with cervical cancer were categorized into two categories as follows: tumor size < 2 cm and ≥ 2 cm: 52 (37.1%) patients were classified into the former, 88 (62.9%) patients were classified into the latter. We found tumor size of 48 (56.5%) patients with cervical AC is larger than or equal to 2 cm, moderately (score 2) or highly expressed (score 3) COPA were found in 46 (95.8%) samples. The remarkably elevated COPA was significantly correlated with tumor size ≥ 2 cm in cervical AC (p = 0.04). Besides, 24 cervical AC patients with bulky tumor (tumor size ≥ 4 cm) were all in the moderate (score 2) and the high (score 3) COPA expression group (Fig. 6A). Meanwhile, tumor size of 40 (72.7%) patients with cervical SCC is larger than or equal to 2 cm, moderately (score 2) or highly expressed (score 3) COPA were found in 24 (60.0%) patients; negatively (score 0) or weakly expressed (score 1) COPA were found in 16 (40.0%) cervical SCC patients (Fig. 6B).
Deep stromal invasion
Deep stromal invasion (DSI) was found in eighty-five (60.7%) cervical cancer samples by histopathology. Specifically, DSI was presented in 46 (54.1%) patients with cervical AC and 39 (70.9%) patients with cervical SCC; Forty-three (93.5%) cervical AC patients with DSI were in the moderate (score 2) COPA expression group or high (score 3) COPA expression group (Fig. 6A). Meanwhile, 23 (59.0%) cervical SCC patients with DSI was in the moderate (score 2) COPA expression group or high (score 3) COPA expression group (Fig. 6B).
Lymphatic vascular space involvement
Lymphatic vascular space involvement (LVSI) was found in 38 (27.1%) cervical cancer samples by histopathology. To be specific, LVSI was presented in 15 (17.6%) patients with cervical AC and 23 (41.8%) patients with cervical SCC; All samples (100.0%) of cervical AC with LVSI were positively expressed COPA, and 93.3% of them were in the moderate (score 2) or high (score 3) COPA expression group (Fig. 6A). In the cervical SCC samples with LVSI, 95.7% were positively expressed COPA, and 69.6% was in the moderate (score 2) or high (score 3) COPA expression group, 26.1% was in the weak (score 1) COPA expression group; negative (score 0) COPA expression was found in only one (4.3%) sample of cervical SCC with LVSI (Fig. 6B).
Validation of COPA staining can serve as an independent prognostic factor for cervical cancer
Through a survival analysis, we discovered that the patients with the moderate (score 2) or strong (score 3) COPA staining possessed worse survival than the patients with the negative (score 0) or weak (score 1) COPA staining (p = 0.0033) (Fig. 7A). Area under the curve (AUC) of ROC curve was 0.626 (Fig. 7B). We performed a univariate Cox regression analysis to assess the influences of patients’ clinicopathological parameters and COPA staining on patients’ OS. The results indicated that moderate (score 2) or strong (score 3) staining of COPA, advanced stage and lymph node metastasis were unfavourable prognostic factors of cervical cancer. Then, we performed a multivariate Cox regression analysis, moderate (score 2) or strong (score 3) COPA staining (HR = 8.946; 95% CI: 1.218–65.714, p = 0.031) and advanced stage were independent prognostic factor of cervical cancer (Table 3).
Depletion COPA inhibited the aggressive behaviors of cervical cancer cells in vitro
SiHa cells (cervical SCC) and HeLa cells (cervical AC) exhibit relatively high endogenous COPA levels. To further elucidate the functional effects of COPA in cervical cancer cells, these two cervical cancer cell lines were used in vitro intervention experiments. We observed the effects of COPA on cell viability and tumorigenic ability in the two COPA knockdown cervical cancer cells. Using Specific siRNAs for COPA, we successfully knocked out COPA in SiHa cells and HeLa cells. The protein abundances of COPA in SiHa cells (Fig. 8A) and HeLa cells (Fig. 8B) were descended robustly in response to transcriptional suppression of COPA. COPA silencing significantly inhibited cell viability in SiHa cells (p < 0.05) (Fig. 8C) and in HeLa cells (p < 0.05) (Fig. 8D). In addition, deficiency of COPA dramatically interfered the tumorigenic ability of SiHa cells (p < 0.05) (Fig. 8E) and of SiHa cells (p < 0.05) (Fig. 8F). These observations provided solid evidence that COPA makes notable contributions to progression of cervical cancer, COPA is a novel pharmacological intervention target for cervical cancer. We demonstrated that deficiency of COPA dramatically interfered the aggressive behaviors of cervical cancer cells in vitro.
Protein mass spectrometry can be generally classified into untargeted (“discovery”) proteomics and targeted (“verification”) proteomics . Targeted proteomics has made steady advances and was named Nature Methods’ 2012 Method of the year because of its special suitability for biomarker candidate verification and clinical applications . More recent, high mass accuracy-based parallel reaction monitoring (PRM)-MS shows outstanding performance characteristics . However, previous tissue-based proteomics of cervical cancer were based on untargeted proteomics. The validation data of these studies were antibody-based technologies such as western blotting and immunohistochemistry, but not targeted proteomics [15, 21, 28,29,30,31,32,33,34,35,36,37]. To our knowledge, this is the first study using the integrated analysis strategy of untargeted proteomics and targeted proteomics to investigate protein expression profiles of cervical cancer. COPA was consistently identified as a potential prognostic biomarker of cervical cancer in two stage proteomic. The validation of the proteomic findings in a large number of clinical samples is beneficial to push forward clinic applications. Therefore, additional IHC analysis of COPA in cervical cancer was performed in a large amount of FFPE samples. In IHC cohort, we not only verified the results of proteomics, but also defined the clinical value of COPA staining in cervical cancer.
Currently, surgery remains the primary treatment of cervical cancer. Adjuvant treatments (radiotherapy or chemoradiotherapy) are determined based on histopathological risk factors [38,39,40]. The prognosis of cervical cancer is closely related to the stage and histopathological risk factors. Accurately predict the prognosis would provide new opportunities for optimizing therapeutic protocols. Parametrical invasion, positive surgical margins, and lymph node metastasis are high-risk factors for cervical cancer. Large tumor (tumor size ≥ 4 cm), DSI, and LVSI are considered intermediate-risk factors [38,39,40]. Recent studies reported that cervical cancer patients with tumor size ≥ 2 cm have unfavorable 5-year overall survival . Our results of IHC analysis demonstrated that COPA protein level was significantly elevated in advanced cervical cancer. Notably, the remarkably elevated COPA was significantly correlated with tumor size ≥ 2 cm in cervical AC (p = 0.04). The moderate (score 2) or strong COPA staining is an unfavourable independent prognostic factor of cervical cancer. These findings further demonstrates that COPA can serve as a valuable prognostic biomarker for cervical cancer and is suitable for clinical application.
COPA is part of the coatomer protein complex I, which is involved in the movement of vesicles within the Golgi and retrograde transport of cargo proteins between endoplasmic reticulum and Golgi . Our functional analysis of COPA and its first neighbor DEPs found the mainly enriched pathways were SNARE interactions in vesicular transport, legionellosis and endocytosis pathway. Therefore, we speculate that COPA could serve as a tumor-promoting gene in cervical cancer trough participate in the regulation of the aforementioned pathways. Recent studies have underscored the role of COPA in a novel carcinogenic mechanism. Downregulation of intracellular COPA levels has been shown to dramatically reduce tumorigenic ability of hepatocarcinoma cells  and inhibit the proliferation of prostate cancer cells . Since, our previous analysis showed that the level of COPA protein was significantly increased in advanced cervical cancer, we also speculated that COPA could be an attractive pharmacological intervention target of cervical cancer. One would query what is the evidence to support this speculation? To this end, we explored COPA functional role in two cervical cancer cell lines. Upon the COPA knockdown, cell viability and tumorigenesis ability of SiHa cells and HeLa cells were dramatically reduced. Therefore, COPA is critical in the progression of cervical cancer. Depletion COPA inhibited the aggressive behaviors of cervical cancer cells in vitro. So far, we have identified COPA as a potential prognostic biomarker and a pharmacological intervention target for cervical cancer on multiple levels.
Indubitably, although the results are promising, some limitations must be clarified in this study. First, we retrospectively analysed the correlation between the IHC staining of COPA and clinicopathological characteristics of patients with cervical cancer, its clinical utility would be more convincing using prospective data to assess. Second, the signalling mechanisms of COPA in cervical cancer is interpreted based on bioinformatics analyses, further experimental studies to validate our findings are imperative. Third, although we confirmed in vitro that COPA is an attractive pharmacological intervention target for cervical cancer, in vivo intervention experiments are also warranted in the future.
To summarize, we analyzed protein expression profiles in cervical cancer through integrating proteomic strategy. COPA was identified as a potential prognostic biomarker for cervical cancer, which was additionally verified by IHC staining. Notable, the critical functional role of COPA in cervical cancer was determined via in vitro experiments. All our data consistently suggested that COPA can serve as a potential prognostic biomarker and pharmacological intervention target for cervical cancer. We hope that these findings will offer some useful insights for subsequent studies and clinical practice.
All data that support the findings of this study are available from the corresponding author upon reasonable request.
Ahmad A, Ansari IA. A Comprehensive Review on Crosstalk of Human Papilloma Virus oncoproteins and developmental/self-renewal pathways during the pathogenesis of uterine cervical cancer. Curr Mol Med. 2020;89:234.
Cohen JD, Li L, Wang Y, Thoburn C, Afsari B, Danilova L, Douville C, Javed AA, Wong F, Mattox A. Detection and localization of surgically resectable cancers with a multi-analyte blood test. Science. 2018;359:926–30.
Lu X, Wang Z, Zhang S, Konijnenberg AP, Ouyang Y, Zhao C, Cai Y. Microscopic phase reconstruction of cervical exfoliated cell under partially coherent illumination. J Biophotonics. 2021;14:e202000401.
Zorzi M, Mistro A, Farruggio A, Bartolomeis L, Frayle-Salamanca H, Baboci L, Bertazzo A, Cocco P, Fedato C, Gennaro M, et al. Use of a high-risk human papillomavirus DNA test as the primary test in a cervical cancer screening programme: a population-based cohort study. BJOG. 2013;120:1260–7.
Han R, Shi Q, Wu S, Yin D, Peng M, Dong D, Zheng Y, Guo Y, Zhang R, Hu F. Dissemination of Carbapenemases (KPC, NDM, OXA-48, IMP, and VIM) Among Carbapenem-Resistant Enterobacteriaceae Isolated From Adult and Children Patients in China. Front Cell Infect Microbiol. 2020;10:314.
Islami F, Fedewa SA, Jemal A. Trends in cervical cancer incidence rates by age, race/ethnicity, histological subtype, and stage at diagnosis in the United States. Prev Med. 2019;123:316–23.
Wu C, Zhu X, Kang Y, Cao Y, Lu P, Zhou W, Zhou H, Zhang Y, Song Y. Epidemiology of Humanpapilloma virus infection among women in Fujian China. BMC Public Health. 2017;18:95.
Zhou J, Wu SG, Sun JY, Li FY, Lin HX, Chen QH, He ZY. Comparison of clinical outcomes of squamous cell carcinoma, adenocarcinoma, and adenosquamous carcinoma of the uterine cervix after definitive radiotherapy: a population-based analysis. J Cancer Res Clin Oncol. 2017;143:115–22.
Sun X, Wan JJ, Qian K. Designed microdevices for in vitro diagnostics. Small Methods. 2017;1:1700196.
Cao J, Shi X, Gurav DD, Huang L, Su H, Li K, Niu J, Zhang M, Wang Q, Jiang M. Metabolic fingerprinting on synthetic alloys for medulloblastoma diagnosis and radiotherapy evaluation. Adv Mater. 2020;32:2000906.
Su H, Li X, Huang L, Cao J, Zhang M, Vedarethinam V, Di W, Hu Z, Qian K. Plasmonic alloys reveal a distinct metabolic phenotype of early gastric cancer. Adv Mater. 2021;33:2007978.
Pappa KI, Christou P, Xholi A, Mermelekas G, Kontostathi G, Lygirou V, Makridakis M, Zoidakis J, Anagnou NP. Membrane proteomics of cervical cancer cell lines reveal insights on the process of cervical carcinogenesis. Int J Oncol. 2018;53:2111–22.
Zhang L, Jin J, Zhang L, Hu R, Gao L, Huo X, Liu D, Ma X, Wang C, Han J, et al. Quantitative analysis of differential protein expression in cervical carcinoma cells after zeylenone treatment by stable isotope labeling with amino acids in cell culture. J Proteomics. 2015;126:279–87.
Xu T, Pang Q, Zhou D, Zhang A, Luo S, Wang Y, Yan X. Proteomic investigation into betulinic acid-induced apoptosis of human cervical cancer HeLa cells. PLoS ONE. 2014;9:e105768.
Choi YP, Kang S, Hong S, Xie X, Cho NH. Proteomic analysis of progressive factors in uterine cervical cancer. Proteomics. 2005;5:1481–93.
Arakel EC, Schwappach B. Formation of COPI-coated vesicles at a glance. J Cell Sci. 2018;131:89.
Song Y, An O, Ren X, Chan THM, Tay DJT, Tang SJ, Han J, Hong H, Ng VHE, Ke X, et al. RNA editing mediates the functional switch of COPA in a novel mechanism of hepatocarcinogenesis. J Hepatol. 2021;74:135–47.
UniProt Consortium T. UniProt: the universal protein knowledgebase. Nucleic Acids Res. 2018;46:2699.
Zhou Q, Andersson R, Hu D, Bauden M, Kristl T, Sasor A, Pawłowski K, Pla I, Hilmersson KS, Zhou M, et al. Quantitative proteomics identifies brain acid soluble protein 1 (BASP1) as a prognostic biomarker candidate in pancreatic cancer tissue. EBioMedicine. 2019;43:282–94.
Bao H, Sin TK, Zhang G. Activin A induces tumorigenesis of leiomyoma via regulation of p38β MAPK-mediated signal cascade. Biochem Biophys Res Commun. 2020;529:379–85.
Peng X, Xu X, Wang Y, Hawke DH, Yu S, Han L, Zhou Z, Mojumdar K, Jeong KJ, Labrie M. A-to-I RNA editing contributes to proteomic diversity in cancer. Cancer Cell. 2018;33:817–28.
Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, Simonovic M, Doncheva NT, Morris JH, Bork P, et al. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019;47:D607-d613.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13:2498–504.
Tang Z, Kang B, Li C, Chen T, Zhang Z. GEPIA2: an enhanced web server for large-scale expression profiling and interactive analysis. Nucleic Acids Res. 2019;47:W556-w560.
Sobsey CA, Ibrahim S, Richard VR, Gaspar V, Mitsa G, Lacasse V, Zahedi RP, Batist G, Borchers CH. Targeted and untargeted proteomics approaches in biomarker development. Proteomics. 2020;20:1900029.
Marx V. Targeted proteomics. Nat Methods. 2013;10:19.
Gallien S, Duriez E, Crone C, Kellmann M, Moehring T, Domon B. Targeted proteomic quantification on quadrupole-orbitrap mass spectrometer. Mol Cell Proteomics. 2012;11:1709–23.
Kontostathi G, Zoidakis J, Anagnou NP, Pappa KI, Vlahou A, Makridakis M. Proteomics approaches in cervical cancer: focus on the discovery of biomarkers for diagnosis and drug treatment monitoring. Expert Rev Proteomics. 2016;13:731–45.
Zhao Q, He Y, Wang X-L, Zhang Y-X, Wu Y-M. Differentially expressed proteins among normal cervix, cervical intraepithelial neoplasia and cervical squamous cell carcinoma. Clin Transl Oncol. 2015;17:620–31.
Escobar-Hoyos LF, Yang J, Zhu J, Cavallo J-A, Zhai H, Burke S, Koller A, Chen EI, Shroyer KR. Keratin 17 in premalignant and malignant squamous lesions of the cervix: proteomic discovery and immunohistochemical validation as a diagnostic and prognostic biomarker. Mod Pathol. 2014;27:621–30.
Wang W, Jia H, Huang J, Liang Y, Tan H, Geng H, Guo L, Yao S. Identification of biomarkers for lymph node metastasis in early-stage cervical cancer by tissue-based proteomics. Br J Cancer. 2014;110:1748–58.
Song JY, Bae HS, Koo DH, Lee JK, Jung HH, Lee KW, Lee NW. Candidates for tumor markers of cervical cancer discovered by proteomic analysis. J Korean Med Sci. 2012;27:1479–85.
Fukushima C, Murakami A, Yoshitomi K, Sueoka K, Nawata S, Nakamura K, Sugino N. Comparative proteomic profiling in squamous cell carcinoma of the uterine cervix. Proteomics-Clin Appl. 2011;5:133–40.
Lomnytska MI, Becker S, Hellman K, Hellström AC, Souchelnytskyi S, Mints M, Hellman U, Andersson S, Auer G. Diagnostic protein marker patterns in squamous cervical cancer. Proteomics-Clin Appl. 2010;4:17–31.
Arnouk H, Merkley MA, Podolsky RH, Stöppler H, Santos C, Álvarez M, Mariategui J, Ferris D, Lee JR, Dynan WS. Characterization of molecular markers indicative of cervical cancer progression. Proteomics-Clinical Applications. 2009;3:516–27.
Zhu X, Lv J, Yu L, Zhu X, Wu J, Zou S, Jiang S. Proteomic identification of differentially-expressed proteins in squamous cervical cancer. Gynecol Oncol. 2009;112:248–56.
Bae SM, Lee C-H, Cho YL, Nam KH, Kim YW, Kim CK, Han BD, Lee YJ, Chun HJ, Ahn WS. Two-dimensional gel analysis of protein expression profile in squamous cervical cancer patients. Gynecol Oncol. 2005;99:26–35.
Huang H, Feng YL, Wan T, Zhang YN, Cao XP, Huang YW, Xiong Y, Huang X, Zheng M, Li YF, et al. Effectiveness of sequential chemoradiation vs concurrent chemoradiation or radiation alone in adjuvant treatment after hysterectomy for cervical cancer: The STARS Phase 3 Randomized Clinical Trial. JAMA Oncol. 2021;7:361–9.
Cibula D, Pötter R, Planchamp F, Avall-Lundqvist E, Fischerova D, Haie Meder C, Köhler C, Landoni F, Lax S, Lindegaard JC, et al. The European Society of Gynaecological Oncology/European Society for Radiotherapy and Oncology/European Society of Pathology guidelines for the management of patients with cervical cancer. Radiother Oncol. 2018;127:404–16.
Koh WJ, Abu-Rustum NR, Bean S, Bradley K, Campos SM, Cho KR, Chon HS, Chu C, Clark R, Cohn D, et al. Cervical Cancer, Version 32019, NCCN Clinical Practice Guidelines in Oncology. J Natl Compr Canc Netw. 2019;17:64–84.
Kato T, Takashima A, Kasamatsu T, Nakamura K, Mizusawa J, Nakanishi T, Takeshima N, Kamiura S, Onda T, Sumi T, et al. Clinical tumor diameter and prognosis of patients with FIGO stage IB1 cervical cancer (JCOG0806-A). Gynecol Oncol. 2015;137:34–9.
Iglesias-Gato D, Wikström P, Tyanova S, Lavallee C, Thysell E, Carlsson J, Hägglöf C, Cox J, Andrén O, Stattin P, et al. The Proteome of Primary Prostate Cancer. Eur Urol. 2016;69:942–52.
We thank Jingjie PTM BioLab (Hangzhou) Co. Ltd., Hangzhou, Zhejiang, China, for technical assistance with LC-MS/MS.
Institutional review board
The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Ethics Committee of Guangdong Provincial People's Hospital and Zhuhai People's Hospital (Zhuhai Hospital affiliated with Jinan University) (protocol code LW- No.1, KY-Q-2021–077-01).
This research was funded by State Key Development Program for Basic Research of China, Grant Number 2017YFA0205200; the National Natural Science Foundation of China, Grant Numbers 81901857 and 82002170; the Science and Technology Planning Project of Guangzhou, Guangdong Province, China, Grant Number 201904010026.
Informed consent was obtained from all subjects involved in the study.
The authors declare no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
. Baseline characteristics of included patients with cervical cancer. Table S2. The histological types of included cervical adenocarcinoma. Figure S1. RNA expression of COPA in cell lines based on HPA database. (A) The cell lines ordered by human tissues and organ; (B) The cell lines ordered by descending RNA expression order. Figure S2. HPA (The Human Protein Atlas) analyses for the COPA predicted that the protein was detected in nucleoplasm, cytosol and Golgi apparatus and was predicted to be secreted. Figure S3. Histograms showing the distributions of ages of the 140 patients with cervical cancer in the IHC analysis. Figure S4. The positions of COPA plotted in the volcano plots of cervical cancer (AC and SCC) vs. healthy controls.
About this article
Cite this article
Bao, H., Li, X., Cao, Z. et al. Identification of COPA as a potential prognostic biomarker and pharmacological intervention target of cervical cancer by quantitative proteomics and experimental verification. J Transl Med 20, 18 (2022). https://doi.org/10.1186/s12967-021-03218-1