Introduction

Molecular Medicine Reports

1791-2997 1791-3004

D.A. Spandidos

10.3892/mmr.2023.12933

MMR-27-2-12933

Articles

Evaluation of the proteomic landscape of HPV E7‑induced alterations in human keratinocytes reveal therapeutically relevant pathways for cervical cancer

Gandhi

Sivasangkary

1 Mohamad Razif

Muhammad Fazril

1 Othman

Shatrah

1 2 Chakraborty

Sajib

3 Nor Rashid

Nurshamimi

1 2

1Department of Molecular Medicine, Faculty of Medicine, University of Malaya, 50603 Kuala Lumpur, Malaysia 2Drug Design and Development Research Group, University of Malaya, 50603 Kuala Lumpur, Malaysia 3Translational System Biology Laboratory, Department of Biochemistry and Molecular Biology, University of Dhaka, Dhaka 1000, People's Republic of Bangladesh

Correspondence to: Dr Nurshamimi Nor Rashid, Department of Molecular Medicine, Faculty of Medicine, University of Malaya, Block F, Jalan Professor Diraja Ungku Aziz, 50603 Kuala Lumpur, Malaysia, E-mail: nurshamimi@um.edu.my

02 2023

09 01 2023

27 2

25052022 18102022

2023

The lack of specific and accurate therapeutic targets poses a challenge in the treatment of cervical cancer (CC). Global proteomics has the potential to characterize the underlying and intricate molecular mechanisms that drive the identification of therapeutic candidates for CC in an unbiased manner. The present study assessed human papillomavirus (HPV)-induced proteomic alterations to identify key cancer hallmark pathways and protein-protein interaction (PPI) networks, which offered the opportunity to evaluate the possibility of using these for targeted therapy in CC. Comparative proteomic profiling of HPV-transfected (HPV16/18 E7), HPV-transformed (CaSki and HeLa) and normal human keratinocyte (HaCaT) cells was performed using the liquid chromatography-tandem mass spectrometry (LC-MS/MS) technique. Both label-free quantification and differential expression analysis were performed to assess differentially regulated proteins in HPV-transformed and -transfected cells. The present study demonstrated that protein expression was upregulated in HPV-transfected cells compared with in HPV-transformed cells. This was probably due to the ectopic expression of E7 protein in the former cell type, in contrast to its constitutive expression in the latter cell type. Subsequent pathway visualization and network construction demonstrated that the upregulated proteins in HPV16/18 E7-transfected cells were predominantly associated with a diverse array of cancer hallmarks, including the mTORC1 signaling pathway, MYC targets V1, hypoxia and glycolysis. Among the various proteins present in the cancer hallmark enrichment pathways, phosphoglycerate kinase 1 (PGK1) was present across all pathways. Therefore, PGK1 may be considered as a potential biomarker. PPI analysis demonstrated a direct interaction between p130 and polyubiquitin B, which may lead to the degradation of p130 via the ubiquitin-proteasome proteolytic pathway. In summary, elucidation of the key signaling pathways in HPV16/18-transfected and -transformed cells may aid in the design of novel therapeutic strategies for clinical application such as targeted therapy and immunotherapy against cervical cancer.

human papillomavirus human keratinocytes E7 proteomics bioinformatics analysis

Malaysian Ministry of Higher Education through the Fundamental Research Grant Scheme

FRGS/1/2016/SKK08/UM/02/16 (FP015-2016)

The present study was supported by the Malaysian Ministry of Higher Education through the Fundamental Research Grant Scheme [grant no. FRGS/1/2016/SKK08/UM/02/16 (FP015-2016)].

Introduction

Human papillomavirus (HPV) infections are strongly associated with cervical cancer (CC), and with the incidence of oropharyngeal cancer (1) and anogenital cancer (2). HPV infections are also typically linked to skin or mucosal lesions, such as warts. HPV affects both men and women; however, the disease burden is commonly observed among women due to their high susceptibility for cervical infections (3). Approximately 90% of CC cases are caused by one or more HPV types, particularly HPV16 (detected in ~50% of cases) and HPV18 (detected in 10–15% of cases) (4). HPV16 and −18 are known as high-risk HPV (HR-HPV) subtypes, alongside another 15 types, including HPV31, −33-35, −39, −45, −51, −52, −56, −58, −59, −68, −73 and −82 (5). Usually, HR-HPV infections persist for 1–2 years. Within this period, the virion replicates for months in host tissue with regular transformation, along with active suppression of both innate and adaptive immune responses (6). HPV relies greatly on its potential to control both viral and host gene expression (7). Viral gene transcript regulation is directly related to keratinocyte differentiation in the host (8). Previous studies have reported that HPV-infected keratinocytes exhibit significantly reduced expression of numerous inflammatory mediators (9,10).

HPV is able to promote immune evasion via the expression of oncogenic proteins, which are responsible for the modulation of several immune mechanisms, including antigen presentation and inflammatory pathways (11). HPV E6 and E7 oncoproteins support viral oncogenicity as they are expressed concurrently with whole tumor progression and deregulate host cell gene expression (12). Both E6 and E7 proteins exhibit numerous functions; however, they are mainly involved in the inactivation of p53 and retinoblastoma protein (pRb) (13). The inactivation of p53 and pRb can result in the loss of control over DNA damage repair and cell cycle regulation. The E7 oncoprotein binds to and induces degradation of pocket proteins (such as pRb, p107 and p130), which leads to the release of E2F transcription factor family proteins and the expression of S-phase genes (14). However, the suppression of pRb family proteins is not entirely influenced by E7, which suggests that other mechanisms are involved. Furthermore, the degradation of p107 and p130 by E7 is also related to the suppression of DREAM complex-mediated genes (13). Various proteins across different pathways, including activator protein-1 (15), HIF1 (16), krüppel-like factor 4 (17), p73 (18) and other host cell factors have been identified as E7 targets. However, their functional interactions that encourage virus replication and carcinogenesis are not yet fully understood (7).

Understanding the pathogenic mechanisms and signaling pathways associated with HR-HPV, along with the identification of differentially expressed protein-protein interactions (PPIs), could aid in the design of novel therapeutic strategies that target HPV-associated CC. Therefore, the present study evaluated changes in protein expression levels between normal human keratinocytes (HaCaT) transfected with recombinant HPV16/18 E7 and HPV-transformed cells (Caski and HeLa) to identify proteins that may contribute to cancer hallmark enrichment. Deciphering alterations in cellular protein interactions and their associated pathways during the viral life cycle is essential to understand the evolution of HPV E7 function in cervical carcinogenesis.

Materials and methods <sec> <title>Cell lines

Human HaCaT keratinocyte cells (Addexbio Technologies), HPV16-positive CC CaSki (passage 9) and HPV18-positive CC HeLa (passage 11) cell lines (both from the Laboratory of Virology, Department of Medical Microbiology, University of Malaya, Kuala Lumpur, Malaysia), were used in the present study. The HaCaT cells were authenticated using STR profiling. CaSki cells carry integrated viral HPV16 DNA. HeLa cells contain integrated viral HPV18 DNA. All cell lines were maintained in T25 flasks with 5 ml Dulbecco's Modified Eagle's Medium (Gibco; Thermo Fisher Scientific, Inc.), 10% heat-inactivated fetal bovine serum (Gibco; Thermo Fisher Scientific, Inc.) and 1% penicillin/streptomycin (Gibco; Thermo Fisher Scientific, Inc.). Cells were maintained under standard incubation conditions of 10% CO₂ with 95% humidity at 37°C.

Transfection of HPV16/HPV18 E7 recombinant plasmids in human keratinocytes

Transfection of HaCaT cells was performed according to the manufacturer's protocol. Briefly, FuGENE^® HD transfection buffer (Promega Corporation) was added to 2 µg of each plasmid [pMSCVpuro; Addgene (Fig. S1) (control); pMSCVpuro-HPV16E7; and pMSCVpuro-HPV18E7], which were diluted in 100 µl serum-free Opti-MEM™ Reduced Serum medium (Gibco; Thermo Fisher Scientific, Inc.). Subsequently, 100 µl transfection mixture containing 2 µg of the nucleic acid was added to HaCaT cells (3×10⁵) grown in a 6-well plate and incubated at 37°C for 48 h. The cells were then washed and maintained using 0.5 µg/ml puromycin (Gibco; Thermo Fisher Scientific, Inc.) for 24 h post-transfection and incubated for a further 48 h to establish stable transformants before harvesting. Transfection efficiency was assessed using reverse transcription-quantitative PCR (RT-qPCR) and the results are presented in Fig. S2.

RNA extraction and RT-qPCR

Total RNA was extracted using the FavorPrep™ Blood/Cultured Cell Total RNA Mini Kit (Favorgen). The cDNA was prepared using the RevertAid™ First Strand cDNA Synthesis Kit (Thermo Fisher Scientific, Inc.) according to the manufacturer's protocol. qPCR was performed using Absolute SYBR Green ROX (Thermo Fisher Scientific, Inc.) in triplicate using an ABI PRISM™ 7900HT sequence detector (Applied Biosystems; Thermo Fisher Scientific, Inc.) (19). The primers used for the qPCR reactions are presented in Table SI. PCR was performed under the following thermocycling conditions: 95°C for 15 min, followed by 40 cycles at 95°C for 15 sec, 58°C for 30 sec and 72°C for 30 sec. Optimization was performed for each primer set and the relative quantification was performed with normalization against β-actin (20).

Protein extraction and quantification

Proteins were extracted using RIPA lysis buffer [NaCl, 150 mM; 1% Triton X-100; Tris-HCl, 50mM (pH, 8.0); 0.5% sodium deoxycholate; 0.1% SDS; Halt™ Protease Inhibitor Cocktail (100X; Thermo Fisher Scientific, Inc.)]. Cells were washed twice with ice-cold PBS and mixed with ice-cold RIPA lysis buffer. The cell lysate was centrifuged at 13,000 × g for 20 min at 4°C. The supernatant, containing soluble proteins, was collected and stored at −80°C until further analysis. The concentration of the extracted proteins was determined using the Bradford assay (21). Bovine serum albumin (MilliporeSigma) was used as a protein standard to construct the calibration curve.

In-solution tryptic digestion and desalting

A total of ~500 µg protein from each sample were digested in 50 mM ammonium bicarbonate digestion buffer together with 100 mM DTT reducing buffer and incubated at 95°C for 5 min. The samples were cooled to ambient temperature and alkylated in 100 mM iodoacetamide buffer for 20 min in the dark. Trypsin is a serine protease that is highly specific, digesting protein into peptides by cutting at the carboxyl side of arginine and lysine residues. Samples were then incubated with trypsin (Promega Corp.) at 37°C for 3 h and another 1 µl trypsin was added for overnight incubation at 30°C. The digested peptides were desalted, lyophilized and stored at −80°C until further experimentation (22).

Protein identification using nano-electrospray ionization-liquid chromatography-tandem mass spectrometry (nano-ESI-LC-MS/MS)

Protein identification using nano-ESI-LC-MS/MS was performed with certain modifications (23). Trypsin-digested peptides (500 µg) were loaded onto a 300 Å, SB-C18, 160 nl enrichment column and 75 µm × 150 mm analytical column (cat. no. G4240-62010; Agilent Technologies, Inc.) with a flow rate of 4 µl/min using a capillary pump and 0.5 µl/min using an Agilent 1200 nano pump. The eluted peptides were subjected to nano-ESI MS/MS using an Agilent 1200 Series HPLC-Chip/MS System Interface, coupled with the Agilent 6550 Q-TOF LC/MS system (Agilent Technologies, Inc.). Injection volume was adjusted to 1 µl/sample. The mobile phases were 0.1% formic acid in water (solution A) and 90% acetonitrile with 0.1% formic acid (solution B) for 4 min and 70% solution B for 3 min, using an Agilent 1200 Series nanoflow LC pump. Ion polarity was set to positive ionization mode, the drying gas flow rate was 5.0 µl/min and the temperature was fixed at 325°C. The fragmentor and capillary voltage were set at 360 and 1,900 V, respectively. The spectra were acquired in MS/MS mode with an MS scan range of 110–3,000 m/z and MS/MS scan range of 50–3,000 m/z. Precursor charge selection was set as a double-, triple- or more than triple-charged state, with the exclusion of precursors 1,221.9906 m/z (z=1) set as reference ions. Data were extracted using MH⁺ (positive ion) mass range between 50–3,200 Da and processed using the Agilent Spectrum Mill MS Proteomics Workbench software packages version B.04.00 (Agilent Technologies, Inc.). The raw data files obtained from the LC-MS/MS were processed using PEAKS DB analysis to evaluate the relative protein abundance.

De novo identification of peptides using PEAKS Studio 7.0

PEAKS DB is a proteomic software package for MS/MS designed for peptide sequencing, protein identification and quantification (24). Peptide identification was performed using automated de novo sequencing in PEAKS Studio 7.0 (Bioinformatics Solution, Inc.) (25). Proteins were identified from HPV-transfected and -transformed cells using the Uniprot Homo sapiens, type Swissprot and Trembl databases, processed with PEAKS 7.0 (Bioinformatics Solution, Inc.) and carbamidomethylation was set as a fixed modification. Subsequently, high-confidence proteins were identified by setting a false discovery rate (FDR) threshold of 1%, unique peptide ≥1 and −10lgP>20.

Label-free quantification (LFQ)

Label-free mass spectrometry-based quantitative approaches provide powerful, fast and low-cost tools for analyzing protein changes in complex biological samples in several large-scale biomarker discovery studies (26). However, to avoid any variation errors in the performance of LC and MS, a carefully controlled normalization step is required. LFQ was performed using PEAKS studio 7.0 (Bioinformatics Solution, Inc.) the abundance of proteins was calculated using normalized spectral protein intensity (LFQ intensity), in which proteins were quantified by comparing the number of identified MS/MS spectra from the same protein in each of the multiple LC-MS/MS data sets. It is possible that an increase in protein abundance results in a higher number of proteolytic peptides and vice versa. In turn, a larger number of proteolytic peptides leads to higher protein sequence coverage and increased number of both identified unique peptides and total MS/MS spectra (spectral count) for each protein (27). HaCaT-pMSCV puro was used as a control to normalize the LFQ of fold changes for significant proteins. The pheatmap package was used to present differentially expressed proteins with hierarchical clustering.

Data acquisition and statistical analysis

Three biological replicates for each cell line, including the control, underwent LC-MS/MS. P-values were obtained using the Benjamini-Hochberg FDR. Statistical analysis was performed using GraphPad Prism version 5.0.0 for Windows (GraphPad Software, Inc.) and the data are presented as the mean ± standard deviation (SD). A comparison between the transfected samples and control was conducted using Student'st-test. and P<0.05 was considered to indicate a statistically significant difference.

Bioinformatics analysis <sec> <title>UpSet plot analysis

An UpSet plot is commonly applied to visualize a dataset with more than three overlapping patterns. It is similar to a Venn diagram but more comprehensive (28). An UpSet plot (https://upset.app/) for HaCaT control, HaCaT-16 E7, HaCaT-18 E7, CaSki and HeLa datasets was generated to identify unique and overlapping proteins across all five cell lines.

Volcano plot analysis

The identified proteins were further analyzed using Perseus (version 1.5.3.0) (29) for differential expression analysis. Firstly, the reverse, site-only and contaminant peptides were removed from the dataset, with missing values input using a normal distribution. This analysis identified significant protein changes using Student's t-test with an FDR<0.05 between HPV-transfected and HPV-transformed cells. The data were presented using a volcano plot, where the x-axis represented the fold change and the y-axis presented the -logP-value.

Molecular signature database (MSigDB) analysis

The MSigDB (30) is one of the largest and most popular repositories of gene sets for use with the Gene Set Enrichment Analysis (GSEA) software tool. The GSEA tool was used to characterize protein expression levels in HaCaT control and HPV-transfected cells with the accession no. MSV000090470 (massive.ucsd.edu). The cancer hallmarks for upregulated proteins in HPV-transfected cells were also identified (31). The P-value summary was used to determine the FDR according to the method described by Benjamini (32). The proteins were identified for each cancer hallmark, sorted by their P-value and FDR values. The top-scoring protein with a summary FDR<0.01 was considered to encompass the definitive cancer hallmark set.

PPI network

Cytoscape 3.8.2 is a freely available platform for network visualization and analysis (Cytoscape_v3.8.2). PPI network construction in Cytoscape requires each protein in the input file to have the same identifiers. Therefore, UniProt ID mapping with Cytoscape was applied to standardize the identifiers. All upregulated proteins in HPV-transfected cells were selected to create a PPI network. The interactors in the network were carefully evaluated to identify potential interactions between the nodes.

Search Tools for the Retrieval of Interacting Genes (STRING) functional enrichment analysis

Gene Ontology (GO) enrichment analysis was applied to identify the function of differentially expressed proteins based on three main aspects as follows: i) Molecular function (MF); ii) cellular component (CC); and iii) biological process (BP) in which they were mutually involved. The list of identified proteins was evaluated using version 10 of the STRING database (http://string10.embl.de) to predict the identified proteins. The interaction score was set to a high confidence level of 0.700. GO terms with FDR<0.01 were defined as the enriched terms for the differentially expressed proteins; Homo sapiens was selected as the organism.

Results <sec> <title>Protein identification in HPV-transfected and -transformed cells

Protein profiling analysis was performed using total protein fractions from each cell type to identify differentially expressed proteins in HPV-transfected human keratinocytes (HaCaT-HPV16/18 E7) and HPV-transformed cells (CaSki and HeLa), with empty vector pMSCV-puro human keratinocytes (HaCaT) used as a control. After quality control assessment of each protein, including average mass, number of peptides and number of uniquely expressed peptides, high confidence proteins (−10lgP>20 considered significant) were identified. A total of 41 proteins were identified in HaCaT cells (Table SII), 29 proteins were identified in HaCaT-pMSCVpuro (control) cells (Table SIII), 85 proteins were identified in HaCaT-HPV16 E7 cells (Table SIV) and 104 proteins were identified in HaCaT-HPV18 E7 cells (Table SV). For HPV-transformed cells, 45 and 39 proteins were identified in CaSki (Table SVI) and HeLa (Table SVII) cells, respectively. The results of the present study demonstrated that glucose-6-phosphate isomerase (GPI) and fructose-bisphosphate aldolase (ALDOA) were the two common proteins among the cell lines. GAPDH, tubulin α-1A chain (TUBA1A), tubulin α chain, heat shock 70kDa protein 8 isoform 2 variant (HSPA8), profilin (PFN1), keratin type I cytoskeletal 17, annexin 5, transketolase, annexin 2 (ANXA2), keratin type II cytoskeletal 8 (KRT8) and peroxiredoxin-1 (PRDX1) were identified in CaSki and HeLa cells. Tubulin α-1B chain and elongation factor 2 (EEF2) were common proteins identified in HaCaT-HPV16/18 E7, CaSki and HeLa cells.

LFQ and heat map clustering analysis between HPV-transfected and -transformed cells

The high-confidence proteins (P=0.01) were identified across HPV-transfected and -transformed cells were subjected to LFQ to assess the differentially expressed proteins with significant fold changes (Table SVIII). HaCaT control cells were used for normalization. A total of 62 differentially expressed proteins were identified. A heat map of differentially expressed proteins was generated with legend color bar, column and row annotations (Fig. 1). All proteins were subsequently grouped based on the hierarchical clustering between HPV-transfected and -transformed cells. The upregulated proteins were clustered at the top of the heat map (blue), whereas the downregulated proteins were clustered at the bottom (red). The results of the present study identified 37 upregulated and 13 downregulated proteins. The hierarchical clustering (column) of HPV-transfected cells demonstrated more proteins in common between HaCaT-HPV16 E7 and HaCaT-HPV18 E7 cells compared with those identified in both HPV-transformed cells. Furthermore, hierarchical clustering (row) resulted in two distinct clusters, separated by upregulated and downregulated proteins based on assessed fold change. Most proteins in HPV-transfected cells were upregulated. EIF4A1, PFN1, HSP90AB1, MYL6, EEF1D, ACTN1, TUBA1A, ANXA1, ANXA2, EEF2, HSPA8, SERPINB5, PPIA, ENO1, ALSOA, TUBB4B, HSPD1, CFL1, GAPDH, NME2, PPA1, MYH9, DOPNI1, TXN, UBB, PKM2, PGK1, HSPA1A, TUBB, S100A10, GSTP1, S100-P, PRDX1, AHNAK, SET, CCT6A and HEL-S-48 were upregulated in HPV-transfected cells compared with in both HaCaT (control) and HPV-transformed cells (Cluster 1). However, Cluster 2 demonstrated a mixed pattern with both upregulated and downregulated proteins. The upregulated proteins in Cluster 2 were HNRNPA1, PKM, DDX21, HNRNPA2B1, NCL, HNRNPD, HNRNPU, RPLP2, HNRNPF, HMGB1, HIST1H1D, HNRNPK and TPI1. The downregulated proteins were HSP70-5, ERBB2, GPI, LDHA, SFN, KRT7, EZR, HSPB1, S100A2, EEF1A1, PGAM1 and EL52.

Volcano plot of proteins with significantly increased expression levels in HPV-transfected and -transformed cells

Proteins in the aforementioned heat map analysis (Fig. 1) were further analyzed using a volcano plot to assess the proteins with the most significantly increased protein expression levels in HPV-transfected and -transformed cells (Fig. 2). The volcano plot presented fold-change differences and protein expression level distribution in HPV-transfected cells vs. HaCaT cells, and in HPV-transformed cells vs. HaCaT cells. Proteins were presented in graphs according to fold change (difference) and significance (−log10 P-value). The 15 upregulated proteins mainly observed in HPV-transfected cells were GAPDH, EEF1D, PRDX1, TUBA1A, AHNAK, EEF2, ANXA2, ANXA1, ACTN1, MYL6, SERPINB5, HSPA8, EIF4A1, HSP90AB1 and PFN1. Among these proteins, EIF4A1, HSP90AB1, HSPA8 and SERPINB5 were presented toward the top right of the plot, which indicated high statistical significance and fold change. EIF4A1 serves a crucial role in the transformation and progression of various types of cancer as it is part of the EIF4F complex that controls initiation rates of pro-oncogenic mRNAs involving the PI3K/Akt/mTOR signaling pathway (33). HSP90AB1, a chaperone protein, is often upregulated in cancerous cells and it is able to stabilize their protein functions with activating mutation (34). HSPA8 is a member of the HSP70 family, which collectively functions as a buffering system for cellular stress which is required for cancer cell survival (35). Fold-change differences in the expression levels of the remaining proteins were not statistically significant.

UpSet intersection plot of HPV-transfected and -transformed cells

The total number of quantified proteins in HaCaT control, HPV-transfected (HaCaT-HPV16 E7 and HaCaT-HPV18 E7) and HPV-transformed (CaSki and HeLa) cells were further evaluated using UpSet plot analysis (Fig. 3). The UpSet intersection plot was presented as a matrix layout to visualize overlaps and differences between qualified proteins across all five cell lines. Dark circles in the matrix indicate sets that are part of the intersection. A high number of linkages indicate a marked association of proteins between the respective cell lines (36). A total of 18 protein linkages were demonstrated between HaCaT-HPV16 E7 and HaCaT-HPV18 E7 cells, the most in the entire analysis (GAPDH, PKM, PPIA, HSPD1, ENO1, TUBB4B, EEF2, ANXA2, ANXA1, TUBA1A, ALDOA, NME2, EEFID, ACTN1, MYL6, CFL1, HSPA8 and HNRNPA1). The second-highest number of protein linkages (n=17) was demonstrated across all five cell lines (GAPDH, PPIA, HSPD1, ENO1, TUBB4B, EEF2, ANXA2, ANXA1, TUBA1A, ALDOA, NME2, EEFID, ACTN1, MYL6, CFL1, HSPA8 and HNRNPA1). Only two protein linkages, KRT7 and EZR, were demonstrated between HeLa and CaSki cells. The lowest number of protein linkages (n=1) was demonstrated in three different groups of linkages as follows: i) HeLa, CaSki, HaCaT control and HaCaT-HPV18 E7; ii) HeLa, CaSki, HaCaT control and HaCaT-HPV16 E7; and iii) HeLa, CaSki, HaCaT-HPV16 E7 and HaCaT-HPV18 E7. The representative chromatograms presented the intensity count for GAPDH, ENO1, ALDOA and HSPA8 proteins with the time taken to pass through the column assessed using the LC-MS/MS analysis of the HPV-transfected (HaCaT-HPV16 E7 and HaCaT-HPV18 E7) and HPV-transformed (CaSki and HeLa) cells (Fig. 4).

Cancer hallmark enrichment of upregulated proteins in HPV-transfected cells

Cancer hallmark enrichment was analyzed based on the P-value and FDR values of the aforementioned upregulated proteins in HaCaT-HPV16/18 E7 cells. The analysis demonstrated that these differentially expressed proteins were significantly enriched across nine pathways (Fig. 5A). A total of 17 of these proteins were involved in the top four enriched pathways: MTORC1 signaling, glycolysis, hypoxia and MYC target VI (Fig. 5B). A total of 10 proteins (PGK1, ALDOA, ENO1, PPIA, GAPDH, HSPD1, PRDX1, PPA1, CCT6A and TXN) were involved with the activation of the mTORC1 complex. A total of seven proteins (PGK1, ALDOA, GAPDH, PRDX1, HSP90AB1, SET and PFN1) were involved with the MYC target VI pathway. Another seven proteins were demonstrated to be involved in cellular response to hypoxia (PGK1, ALDOA, ENO1, PPIA, HSPD1, ANXA2 and EIF4A1). Furthermore, PGK1, ALDOA, ENO1, PPIA, GAPDH, PKM and MYH9 proteins were involved in glycolysis. Notably, PGK1, ALDOA and ENO1 were involved in all four pathways. Other noteworthy cancer hallmark enrichment pathways were apical junction, reactive oxygen species (ROS) pathway, P13K/AKT/mTOR signaling pathway, unfolded protein response and fatty acid metabolism. A total of three proteins (PGK1, ANXA2 and ACTN1) were demonstrated to be engaged in the apical junction pathway. PGK1, PPA1 and PKM were linked to the ROS pathway. PGK1, ACTN1 and EEF2 were linked with the PI3K/AKT/mTOR signaling pathway. PGK1, HSP90AB1 and S100A1 were involved with the unfolded protein response, a cellular stress response related to the endoplasmic reticulum. Finally, PGK1, ENO1 and HSPA1A were associated with the metabolism of fatty acids. Overall, the results demonstrated that PGK1 expression was upregulated and was involved in all cancer hallmark pathways.

Interactions of upregulated proteins in HPV-transfected cells and HPV-transformed cells with p130, involucrin (IVL) and keratin 10

The upregulated proteins in HPV-transfected and -transformed cells that interacted with p130, IVL and keratin 10 were S100A10, EEF1D, ANXA2, GSTP1, TUBB4B, EEF2, AHNAK, CFL1, TUBA1A, HSP90AB1, ACTN1, PFN1, HSPA1A, HSPA8, PKM, MYH9, SET, NME2, ANXA1, UBB, EIF4A1, ENO1, PGK1, ALDOA, TXN, PRDX1, PPIA, PPA1, CCT6A, GAPDH, HSPD1, TUBB, MYL6 and SERPINB5.

The PPI network for upregulated proteins in HaCaT-HPV16 E7 and HaCaT-HPV18 E7 cells were compared with the target protein network, p130 (also known as RBL2), IVL and keratin 10. p130 was chosen based on its host protein localization, whereas keratin 10 and IVL were selected as protein markers for cellular differentiation. Based on the heat map analysis, PPIs between HPV-transfected and -transformed cells, with their respective log₂ fold changes, were compared. The PPI network demonstrated that p130 only interacted with polyubiquitin-B protein (UBB), which suggested that UBB may be subjected to p130-mediated degradation (Fig. 6). Furthermore, the PPI network also demonstrated that IVL interacted with a calcium-binding protein, S100A10.

GO analysis

All 62 proteins from the LFQ analysis were subjected to functional classification annotation. GO analysis was performed to generate classification clusters in the categories Biological Process (BP), Cellular Component (CC) and Molecular Function (MF). The top upregulated and downregulated differentially expressed proteins linked to BP (Fig. 7A and B), MF (Fig. 7C and D) and CC (Fig. 7E and F) were presented. The subcategories and the number of proteins across 20 GO terms were also labelled. For example, in the BP category, ‘cellular process’ and ‘regulation of biological process’ both exceeded 30 upregulated proteins (Fig. 7A), whereas the ‘localization’ subcategory had eight downregulated proteins (Fig. 7B). In Fig. 7C, top 20 GO terms in the category MF for the upregulated proteins are provided. A high protein count was found in binding and protein binding categories with 34 and 28 proteins, respectively. The downregulated proteins in MF were enriched in 13 GO terms. Protein and enzyme binding categories had a high protein count with 12 and 11 proteins, respectively (Fig. 7D). The top 20 GO terms corresponding to the CC category for the upregulated proteins are provided in Fig. 7E. Cell and cytoplasm had the highest enrichment with 38 proteins, followed by intracellular organelles, cytosol and membrane-bounded organelle with 30, 31 and 32 proteins, respectively. The downregulated proteins were enriched in 12 GO terms corresponding to CC with the cytosol exhibiting the highest enrichment with 12 proteins, followed by cytoplasmic vesicle with 7 proteins and myelin sheath with 5 proteins (Fig. 7F).

The number of proteins in each subcategory was determined; however, the value must be compared with the actual number of occurrences and the expected number of occurrences for each category to draw a reasonable conclusion.

Discussion

HPV E7 was the first oncogene of all the HPV oncogenes to be identified (15). E7 serves a crucial role in driving cells towards cancer and it may trigger cancer characteristics during the process of viral genome replication. Therefore, HPV E7 gene manipulation may be an effective therapy in CC. Furthermore, transfection of keratinocytes with the HPV E7 oncogene is a suitable model to study the oncogenic changes induced by HPV infection.

In the present study, proteomics technology was used to evaluate the protein expression profiles of HR-HPV E7 types (HPV16 and 18) in HPV-transfected and -transformed cells. HaCaT cells were transfected with recombinant HPV 16/18 E7 and HPV-transformed cells CaSki and HeLa were used to represent native infections with HPV16 and 18, respectively. Label-free proteomics was performed to profile their protein contents. The results of the present study demonstrated that GPI and ALDOA were common proteins identified in all cell lines, with these proteins being involved in the catalytic activity of glycolysis (37,38). Glycolysis in tumor cells provides energy to support both rapid proliferation and increased metabolic requirements for macromolecule synthesis. GPI is a dimeric enzyme that acts in the second step of glycolysis, where it catalyzes the conversion of glucose-6-phosphate to fructose-6-phosphate (37). Furthermore, it also functions as a cytokine/growth factor induced by c-Myc and HIF-1 (39), and is markedly expressed in numerous types of cancer, such as bladder, colon, stomach, kidney, lung and ovarian cancer, and lymphoma (40). ALDOA functions as a key enzyme catalyzing the reversible reaction of fructose-1,6-bisphosphate to glyceraldehyde-3-phosphate and dihydroxyacetone phosphate in glycolysis. It also serves an essential role in ATP biosynthesis, is ubiquitous in all organs or cells and is upregulated in numerous types of cancer, including cervical adenocarcinoma (41).

The heat map generated based on the LFQ analysis performed in the present study demonstrated that 37 proteins (Cluster 1) were upregulated in HPV-transfected cells compared with in control and HPV-transformed cells. These cluster 1 proteins were further analyzed using volcano plot analysis, which demonstrated that GAPDH, EEF1D, PRDX1, TUBA1A, AHNAK, EEF2, ANXA2, ANXA1, ACTN1, MYL6, SERPINB5, HSPA8, EIF4A1, HSP90AB1 and PFN1 proteins were significantly upregulated in transfected cells. A previous study reported that TUBA1A, PFN1 and ANXA2 proteins were upregulated in cervical squamous cell carcinoma, and HSPA8 and KRT8 were downregulated (42). EEF1D has previously been reported to be upregulated in various types of cancer, such as colorectal carcinoma, esophageal carcinoma, glioblastoma, glioma, liver, lymphoma, medulloblastoma, melanoma, oral squamous cell carcinoma, osteosarcoma, prostate and papillary renal cell carcinoma, thus indicating that the protein may act as a promoter for cell proliferation and tumor growth (43). AHNAK acts as a tumor suppressor via activation of the TGFβ/Smad3 signaling cascade, which arrests the cell cycle in the G₀/G₁ phase and downregulates c-Myc expression during cell growth (44). It is closely associated with metastasis of aggressive tumors. A previous study reported unregulated ANXA1 protein expression in cells transfected with HPV16 E6/E7, which suggested its involvement in HPV-mediated carcinogenesis (45). Furthermore, SERPINB5 is differentially expressed in different cancer types. It is upregulated in gastric adenocarcinoma, and breast, colon, gallbladder, ovarian and pancreatic cancer, but downregulated in prostate and gastric cancer (46). SERPINB5 exhibits tumor-suppressive properties due to its nuclear localization, binding to chromatin and inhibiting cancer cell metastasis (47). Overexpression of EIF4A1 has been reported to promote CC progression due to its ATP-dependent RNA helicase activity in mRNA translation of oncoproteins involved in cell apoptosis and proliferation (48). HSP90AB1 is required for cancer cell invasion and migration, and its protein expression levels have been reported to be upregulated in HPV-transfected cells (49).

There were four significant cancer hallmark pathways identified through cancer hallmark enrichment analysis, the mTORC1 signaling pathway, MYC target VI, hypoxia and glycolysis. The mTORC1 signaling pathway was the prominent hallmark in the present study. In cancer, cells often use the mTOR signaling pathway as a mechanism to enhance their proliferation. mTOR is a Ser/Thr kinase that performs various functions, and is associated with growth (increase in cell mass and size), proliferation, survival, autophagy, metabolism and cytoskeletal organization (50). mTOR activity has been reported to be dysregulated in numerous types of cancer as it serves a vital role in the autophagy of tumor cells. Notably, the inhibitory effect of rapamycin on mTOR activity may increase cell autophagic flux and result in decreased tumor growth. Furthermore, previous studies have demonstrated that rapamycin may promote the formation of autophagosomes and induce autophagosome-lysosome fusion (51,52). Ji and Zheng (53) reported that the mTOR signaling pathway activated cervical carcinoma, and mTOR-specific small interfering RNA was revealed to effectively suppress HeLa cell proliferation via inhibiting the cell cycle and increasing apoptosis, which is similar to the mechanism of action of rapamycin. Rapamycin is a highly specific inhibitor of mTOR that has been used to impede cell proliferation (54). HPV-transfected cells exhibit highly reduced pRB protein expression levels as a consequence of functional E6 and E7 expression (55). It was previously reported that rapamycin resistance in the proliferation of human keratinocytes expressing HPV16 was associated with the ability of E7 to induce pRb degradation. HPV16 E7 was also reported to have conferred resistance to the anti-proliferative effect of rapamycin. This was also associated with the integrity of the LxCxE motif, which has been reported to affect rapamycin resistance (56).

The second significant hallmark was the MYC target VI pathway. MYC is a transcription factor that regulates multiple human genes that promote cell proliferation (57). It also affects apoptosis via alterations to the pro-and anti-apoptotic members of the BCL-2 family, activates telomerase and regulates the expression of vascular endothelial growth factor, which is associated with angiogenesis (58). These downstream targets make MYC one of the most influential oncogenes. Both mTORC1 signaling and MYC target VI (59) are involved in cell proliferation.

The third cancer hallmark was the hypoxia signaling pathway, which is governed by HIF stabilization. While adapting to hypoxia, tumor cells can become more aggressive and become resistant to therapeutics. Hypoxia induces changes to gene expression and the subsequent proteome changes can have notable effects on various functions, which may negatively affect patient prognosis (60). Notably, slowly dividing cells in hypoxic regions are able to escape most cytotoxic drugs, as these treatments target rapidly dividing cells. Cancer stem cells may also be present in poorly hypoxic regions, thus ensuring epithelial-to-mesenchymal transition (61). Tumor cell survival under hypoxia or substances which block HIF-1α and HIF-2α-linked signaling pathways makes tumor cells adapt to hypoxia. They do so by inducing metabolic reprogramming, improving the survival of tumor cells, and supporting both angiogenesis and metastasis (62,63). Previous studies have reported that high glucose concentrations (25 mM), which are a common occurrence in the blood of patients with uncontrolled diabetes, may efficiently counteract hypoxic E6/E7 repression. The basis of glucose-linked effects on gene expression is complex, and may involve epigenetic mechanisms and specific transcription factors, such as MondoA/ChREBP-Mlx, NF-κB, c-Myc and SP1 (64). It is necessary to study how E6/E7 repression under hypoxia influences viral antigen presentation on HPV-positive cancer cells, thereby assisting their escape from host immune defense mechanisms.

Glycolysis was the final cancer hallmark addressed in the present study. Viral proteins regulate the cell cycle via interacting with the tumor suppressor proteins p53 and pRB. HPV E6/E7, as well as E5 and E2, favor the Warburg effect and can contribute to radioresistance and chemoresistance, supporting glycolytic enzyme activities, Krebs cycle and respiratory chain inhibition (65). These processes lead to the accelerated production of ATP, which may satisfy the energy demands of cancer cells during proliferation. In this manner, HPV proteins may promote cancer hallmarks; however, it is also possible that during early HPV infection, the Warburg effect may aid efficient viral replication (66).

The present study demonstrated that PGK1 was present in all cancer hallmark enrichment pathways, which indicated it as a potential biomarker. Previous studies reported that the protein expression levels of PGK1 were elevated in breast cancer (67), astrocytoma (68), metastatic colon cancer (69) and pancreatic ductal adenocarcinoma (70). Furthermore, its mRNA expression levels have been shown to be elevated in gastric cancer (71). PGK1 is an essential enzyme in aerobic glycolysis, which catalyzes the reversible transfer of a phosphate group from 1,3-bisphosphoglycerate to ADP, thus producing 3-phosphoglycerate and ATP. PGK1 can affect the function of some transcription factors, such as β-catenin, a tumor-associated oncoprotein (72). PGK1 is the upstream regulator of β-catenin, which affects tumor growth, proliferation, invasion, metastasis, angiogenesis and drug resistance (69,73,74). Furthermore, PGK1 serves an important role in the tumor occurrence and progression not only as a metabolic enzyme, but also as a protein kinase. Mitochondrial PGK1 activates pyruvate dehydrogenase kinase isoenzyme 1, through which tumor cells are able to inhibit mitochondrial pyruvate metabolism and promote the Warburg effect (75). Abnormal expression of PGK1 has been detected in tumor tissues, and also in peripheral blood and saliva samples of patients (76). Therefore, PGK1 may be considered a potential target for tumor therapy and may become a popular molecule in tumor therapy research. However, the role of PGK1 in different tumors may vary according to its tissue specificity and associated level of expression. Furthermore, the development of therapeutic drugs targeting PGK1 according to its function is also an important consideration. Therefore, PGK1 has a broad research potential in cancer, especially as a therapeutic target for cervical cancer.

The PPI network analysis performed in the present study demonstrated that p130 interacted exclusively with UBB, which is subjected to p130-mediated degradation via proteasomal degradation (77). The UBB protein is part of the ubiquitin-proteasome system (UPS) that is associated with the degradation of various intracellular proteins in eukaryotic cells (78). It is also involved in cellular signaling pathways, such as cell cycle control, cell survival, proliferation, transcription, DNA repair, apoptosis, cellular metabolism, membrane trafficking and ubiquitination, which are vital for the immune response (79). E7 inactivates most cellular substrates via the interaction of UPS components, leading to their degradation at the proteasome (80). The present study demonstrated that IVL interacted with a calcium-binding protein, S100A10. This protein appears as a small dimeric helix-loop-helix tightly associated with ANXA2 and S100A10 tends to be degraded in the absence of ANXA2. It is involved in the intracellular post-entry trafficking of several membrane-bound proteins (81). It has also been reported to mediate the migration of macrophages to the tumor site (82).

The present study also investigated the functional classification of differentially expressed proteins identified using LFQ analysis to create three distinct clusters: BP, CC and MF. The present study was able to distinguish proteins according to the main clusters affected by HPV infection. Due to its interactions with a wide range of proteins, the effects of E7 were demonstrated across numerous cellular processes, including viral replication, transformation, cell cycle and cell death (83,84). The active site for binding of tumor suppressors with HPV E7 oncoprotein is within conserved region 2, encompassing the LxCxE motif. This motif is responsible for binding with cellular targets. Conserved regions 2 and 3 of HPV are responsible for the degradation of tumor suppressors that ultimately lead to inhibition of cell cycle arrest (85).

In conclusion, developments in HPV research have continued to identify numerous mechanisms that can be exploited by the virus to overcome cellular growth controls. At present, development of accurate therapeutic targets for CC treatment still poses a challenge. The results of the present study demonstrated the importance of elucidating the involvement of multiple pathways perturbed by HPV infection. Furthermore, it was demonstrated that PGK1 was present across all cancer hallmark enrichment pathways. Further investigations are required to identify other PGK1 functions, which could further evaluate its potential as a therapeutic target in CC.

Supplementary Material

Supporting Data

Acknowledgements

Not applicable.

Availability of data and materials

The data generated in the present study maybe obtained from the MassIVE database (massive.ucsd.edu) under accession no. MSV000079070. The direct URL is ftp://massive.ucsd.edu/MSV000090470/. The datasets were summarized in Table SI, Table SII, Table SIII, Table SIV, Table SV, Table SVI. The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Authors' contributions

SG, NFMR, SO, SC and NNR contributed to the conception and design of the present study. Material preparation, data collection and analysis were performed by SG. NNR and MFMR confirm the authenticity of all the raw data. The first draft of the manuscript was written by SG. NNR, MFMR and SO contributed to acquisition of data. SC assisted with data analysis and interpretation. NNR, MFMR, SO and SC critically reviewed the manuscript. All authors have read and approved the final manuscript.

Ethics approval and consent to participate

Not applicable.

Patient consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

References 1

Berman

Schiller

Human papillomavirus in cervical cancer and oropharyngeal cancer: One cause, two diseases

Cancer123221922292017

10.1002/cncr.30588

28346680

Sand

Munk

Jensen

Svahn

Frederiksen

Kjaer

Long-Term risk for noncervical anogenital cancer in women with previously diagnosed high-grade cervical intraepithelial neoplasia: A danish nationwide cohort study

Cancer Epidemiol Biomarkers Prev25109010972016

10.1158/1055-9965.EPI-15-1291

27358257

Fuller

Hinyard

Factors associated with HPV vaccination in young males

J Community Health42112711322017

10.1007/s10900-017-0361-4

28432549

Haverkos

O'Mara

Co-carcinogenesis: Human papillomaviruses, coal tar derivatives, and squamous cell cervical cancer

Front Microbiol822532017

10.3389/fmicb.2017.02253

29180993

Muñoz

Bosch

De Sanjosé

Herrero

Castellsagué

Shah

Snijders

Meijer

International Agency for Research on Cancer Multicenter Cervical Cancer Study Group

Epidemiologic classification of human papillomavirus types associated with cervical cancer

N Engl J Med3485185272003

10.1056/NEJMoa021641

12571259

Della Fera

Warburton

Coursey

Khurana

McBride

Persistent human papillomavirus infection

Viruses133212021

10.3390/v13020321

33672465

Songock

Kim

Bodily

The human papillomavirus E7 oncoprotein as a regulator of transcription

Virus Res23156752017

10.1016/j.virusres.2016.10.017

27818212

Graham

Faizo

AAA

Control of human papillomavirus gene expression by alternative splicing

Virus Res23183952017

10.1016/j.virusres.2016.11.016

27867028

Hong

Akhavan

Lee

Kim

Kang

Park

Shin

Proinflammatory cytokine TNFα promotes HPV-associated oral carcinogenesis by increasing cancer stemness

Int J Oral Sci1232020

10.1038/s41368-019-0069-7

31911577

Richards

Wasson

Watherston

Doble

Eric Blair

Wittmann

Macdonald

The human papillomavirus (HPV) E7 protein antagonises an Imiquimod-induced inflammatory pathway in primary human keratinocytes

Sci Rep5129222015

10.1038/srep12922

26268216

de Freitas

de Oliveira

THA

Barros

JrVenuti

hrHPV E5 oncoprotein: Immune evasion and related immunotherapies

J Exp Clin Cancer Res36712017

10.1186/s13046-017-0541-1

28545552

Yeo-The

NSL

Ito

Jha

High-risk human papillomaviral oncogenes E6 and E7 target key cellular pathways to achieve oncogenesis

Int J Mol Sci1917062018

10.3390/ijms19061706

Fischer

Uxa

Stanko

Magin

Engeland

Human papilloma virus E7 oncoprotein abrogates the p53-p21-DREAM pathway

Sci Rep726032017

10.1038/s41598-017-02831-9

28572607

Bienkowska-Haba

Luszczek

Zwolinska

Scott

Sapp

Genome-Wide transcriptome analysis of human papillomavirus 16-infected primary keratinocytes reveals subtle perturbations mostly due to E7 protein expression

J Virol94e01360192020

10.1128/JVI.01360-19

31748387

Pal

Kundu

Human papillomavirus E6 and E7: The cervical cancer hallmarks and targets for therapy

Front Microbiol1031162020

10.3389/fmicb.2019.03116

32038557

Cuninghame

Jackson

Zehbe

Hypoxia-inducible factor 1 and its role in viral carcinogenesis

Virology4564573703832014

24698149

Gunasekharan

Andrade

Laimins

Post-Transcriptional regulation of KLF4 by high-risk human papillomaviruses is necessary for the differentiation-dependent viral life cycle

PLoS Pathog12e10057472016

10.1371/journal.ppat.1005747

27386862

Sen

Ganguly

Modulation of DNA methylation by human papillomavirus E6 and E7 oncoproteins in cervical cancer

Oncol Lett1511222018

29285184

Gandhi

Nor Rashid

Mohamad Razif

Othman

Proteasomal degradation of p130 facilitate cell cycle deregulation and impairment of cellular differentiation in high-risk Human Papillomavirus 16 and 18 E7 transfected cells

Mol Biol Rep48512151332021

10.1007/s11033-021-06509-4

34169395

Schmittgen

Livak

Analyzing real-time PCR data by the comparative C(T) method

Nat Protoc3110111082008

10.1038/nprot.2008.73

18546601

Olson

BJSC

Assays for determination of protein concentration

Curr Protoc Pharmacol73A.3A.1A.3A.322016

10.1002/cpph.3

27248579

Tan

Revisiting Notechis scutatus venom: On shotgun proteomics and neutralization by the ‘bivalent’ Sea Snake Antivenom

J Proteomics14433382016

10.1016/j.jprot.2016.06.004

27282922

Zainal Abidin

Rajadurai

Chowdhury

Ahmad Rusmili

Othman

Naidu

Proteomic characterization and comparison of malaysian tropidolaemus wagleri and cryptelytrops purpureomaculatus venom using shotgun-proteomics

Toxins (Basel)82992016

10.3390/toxins8100299

27763534

Zhang

Xin

Shan

Chen

Xie

Yuen

Zhang

Lajoie

PEAKS DB: de novo sequencing assisted database search for sensitive and accurate peptide identification

Mol Cell Proteomics11M111.0105872012

10.1074/mcp.M111.010587

Zhang

Hendrie

Liang

Doherty-Kirby

Lajoie

PEAKS: Powerful software for peptide de novo sequencing by tandem mass spectrometry

Rapid Commun Mass Spectrom17233723422003

10.1002/rcm.1196

14558135

Levin

Schwarz

Wang

Leweke

Bahn

Label-free LC-MS/MS quantitative proteomics for large-scale biomarker discovery in complex samples

J Sep Sci30219822032007

10.1002/jssc.200700189

17668910

Washburn

Wolters

Yates

Large-scale analysis of the yeast proteome by multidimensional protein identification technology

Nat Biotechnol192422472001

10.1038/85686

11231557

Alsallakh

Aigner

Miksch

Hauser

Radial sets: Interactive visual analysis of large overlapping sets

IEEE Trans Vis Comput Graph19249625052013

10.1109/TVCG.2013.184

24051816

Tyanova

Temu

Sinitcyn

Carlson

Hein

Geiger

Mann

Cox

The Perseus computational platform for comprehensive analysis of (prote)omics data

Nat Methods137317402016

10.1038/nmeth.3901

27348712

Liberzon

Birger

Thorvaldsdottir

Ghandi

Mesirov

Tamayo

The molecular signatures database (MSigDB) hallmark gene set collection

Cell Syst14174252015

10.1016/j.cels.2015.12.004

26771021

Subramanian

Tamayo

Mootha

Mukherjee

Ebert

Gillette

Paulovich

Pomeroy

Golub

Lander

Mesirov

Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles

Proc Natl Acad Sci USA10215545155502005

10.1073/pnas.0506580102

16199517

Benjamini

Hochberg

Controlling the false discovery rate: A practical and powerful approach to multiple testing

J R Statist Soc B572893001995

Lin

Zhang

Cai

Liang

Chen

Qin

Han

Yuan

Liu

Systematic analysis of gene expression alteration and co-expression network of eukaryotic initiation factor 4A-3 in cancer

J Cancer9456845772018

10.7150/jca.27655

30588240

Zeng

Wang

Hsp90 up-regulates PD-L1 to promote HPV-positive cervical cancer via HER2/PI3K/AKT pathway

Mol Med271302021

10.1186/s10020-021-00384-2

34666670

Shan

Zhou

Zhang

Identification of HSPA8 as a candidate biomarker for endometrial carcinoma by using iTRAQ-based proteomic analysis

Onco Targets Ther9216921792016

27110132

Potriquet

Laohaviroj

Bethony

Mulvenna

A modified FASP protocol for high-throughput preparation of protein samples for mass spectrometry

PLoS One12e01759672017

10.1371/journal.pone.0175967

28750034

Fermo

Vercellati

Marcello

Zaninoni

Aytac

Cetin

Capolsini

Casale

Paci

Zanella

Clinical and molecular spectrum of glucose-6-phosphate isomerase deficiency. report of 12 new cases

Front Physiol104672019

10.3389/fphys.2019.00467

31133865

Kawai

Uemura

Munakata

Takahashi

Haraguchi

Nishimura

Hata

Matsuda

Ikenaga

Murata

Fructose-bisphosphate aldolase A is a key regulator of hypoxic adaptation in colorectal cancer cells and involved in treatment resistance and poor prognosis

Int J Oncol505255342017

10.3892/ijo.2016.3814

28000858

Ždralević

Marchiq

de Padua

MMC

Parks

Pouysségur

Metabolic plasiticy in cancers-distinct role of glycolytic enzymes GPI, LDHs or membrane transporters MCTs

Front Oncol73132017

10.3389/fonc.2017.00313

29326883

Han

Deng

Sun

Luo

Liang

Zhang

Peng

Tian

GPI is a prognostic biomarker and correlates with immune infiltrates in lung adenocarcinoma

Front Oncol117526422021

10.3389/fonc.2021.752642

34912709

Saito

Takasawa

Aoyama

Akimoto

Ota

Magara

Murata

Hirohashi

Hasegawa

Aldolase A promotes epithelial-mesenchymal transition to increase malignant potentials of cervical adenocarcinoma

Cancer Sci111307130812020

10.1111/cas.14524

32530543

Saritha

Abdul

Arun

Meenakshi

Analysis of differentially expressed proteins in the exfoliated cells of normal and squamous cell carcinoma of the uterine cervix to define candidate markers for cervical cancer

Int J Bioc Biotechnol56266362016

Peng

Gao

Chen

Shen

Han

Chen

Lin

Chen

Kang

The role of EEF1D in disease pathogenesis: A narrative review

Ann Transl Med916002021

10.21037/atm-21-5025

34790806

Sohn

Shin

Yoo

Goh

Lee

Bae

Ahnak promotes tumor metastasis through transforming growth factor-β-mediated epithelial-mesenchymal transition

Sci Rep8143792018

10.1038/s41598-018-32796-2

30258109

Calmon

Sichero

Boccardo

Villa

Rahal

HPV16 E6 regulates annexin 1 (ANXA1) protein expression in cervical carcinoma cell lines

Virology49635412016

10.1016/j.virol.2016.05.016

27240147

Manawapat-Klopfer

Thomsen

Martus

Munk

Russ

Gmuender

Frederiksen

Haedicke-Jarboui

Stubenrauch

Kjaer

Iftner

TMEM45A, SERPINB5 and p16INK4A transcript levels are predictive for development of high-grade cervical lesions

Am J Cancer Res6152415362016

27508094

Chang

Liu

Ragunanan

Shiue

SERPINB5 Expression: Association with CCRT Response and Prognostic Value in Rectal Cancer

Int J Med Sci153763842018

10.7150/ijms.22823

29511373

Liang

Zhou

Chen

Wen

Downregulation of eukaryotic initiation factor 4A1 improves radiosensitivity by delaying DNA double strand break repair in cervical cancer

Oncol Lett14697669822017

29163714

Wang

Deng

Mou

Liu

Wang

Hsp90ab1 stabilizes LRP5 to promote epithelial-mesenchymal transition via activating of AKT and Wnt/β-catenin signaling pathways in gastric cancer progression

Oncogene38148915072019

10.1038/s41388-018-0532-5

30305727

Saxton

Sabatini

mTOR signaling in growth, metabolism, and disease

Cell1689609762017

10.1016/j.cell.2017.02.004

28283069

Rezazadeh

Norooznezhad

Mansouri

Jahani

Mostafaie

Mohammadi

Modarressi

Rapamycin reduces cervical cancer cells viability in hypoxic condition: Investigation of the role of autophagy and apoptosis

Onco Targets Ther13423942472020

10.2147/OTT.S249985

32547058

Dossou

Basu

The emerging roles of mTORC1 in macromanaging autophagy

Cancers (Basel)1114222019

10.3390/cancers11101422

31554253

Zheng

Activation of mTOR signaling pathway contributes to survival of cervical cancer cells

Gynecol Oncol1171031082010

10.1016/j.ygyno.2009.12.020

20102778

Lamming

Inhibition of the mechanistic target of rapamycin (mTOR)-rapamycin and beyond

Cold Spring Harb Perspect Med6a0259242016

10.1101/cshperspect.a025924

27048303

Yim

Park

The role of HPV E6 and E7 oncoproteins in HPV-associated cervical carcinogenesis

Cancer Res Treat373193242005

10.4143/crt.2005.37.6.319

19956366

Rabachini

Boccardo

Andrade

Perez

Nonogaki

Cuccovia

Villa

HPV-16 E7 expression up-regulates phospholipase D activity and promotes rapamycin resistance in a pRB-dependent manner

BMC Cancer184852018

10.1186/s12885-018-4392-8

29703186

Chanvorachote

Sriratanasak

Nonpanya

C-myc contributes to malignancy of lung Cancer: A potential anticancer drug target

Anticancer Res406096182020

10.21873/anticanres.13990

32014901

Chen

Olopade

MYC and breast cancer

Genes Cancer16296402010

10.1177/1947601910378691

21779462

Meyer

Penn

Reflecting on 25 years with MYC

Nat Rev Cancer89769902008

10.1038/nrc2231

19029958

Roma-Rodrigues

Mendes

Baptista

Fernandes

Targeting tumor microenvironment for cancer therapy

Int J Mol Sci208402019

10.3390/ijms20040840

30781344

Shibue

Weinberg

EMT, CSCs, and drug resistance: The mechanistic link and clinical implications

Nat Rev Clin Oncol146116292017

10.1038/nrclinonc.2017.44

28397828

Balamurugan

HIF-1 at the crossroads of hypoxia, inflammation, and cancer

Int J Cancer138105810662016

10.1002/ijc.29519

25784597

Semenza

Hypoxia-inducible factors: Mediators of cancer progression and targets for cancer therapy

Trends Pharmacol Sci332072142012

10.1016/j.tips.2012.01.005

22398146

Havula

Hietakangas

Glucose sensing by ChREBP/MondoA-Mlx transcription factors

Semin Cell Dev Biol236406472012

10.1016/j.semcdb.2012.02.007

22406740

Martínez-Ramírez

Carrillo-García

Contreras-Paredes

Ortiz-Sánchez

Cruz-Gregorio

Lizano

Regulation of cellular metabolism by high-risk human papillomaviruses

Int J Mol Sci1918392018

10.3390/ijms19071839

29932118

Medda

Duca

Chiocca

Human papillomavirus and cellular pathways: Hits and targets

Pathogens102622021

10.3390/pathogens10030262

33668730

Zhang

Tai

Wong

Chiu

Sethi

Koay

Proteomic study reveals that proteins involved in metabolic and detoxification pathways are highly expressed in HER-2/neu-positive breast cancer

Mol Cell Proteomics4168616962005

10.1074/mcp.M400221-MCP200

16048908

Yan

Yang

Xiao

Zou

Zhang

Liu

Over-expression of cofilin-1 and phosphoglycerate kinase 1 in astrocytomas involved in pathogenesis of radioresistance

CNS Neurosci Ther187297362012

10.1111/j.1755-5949.2012.00353.x

22742733

Ahmad

Glatzle

Bajaeifer

Bühler

Lehmann

Königsrainer

Vollmer

Sipos

Ahmad

Northoff

Phosphoglycerate kinase 1 as a promoter of metastasis in colon cancer

Int J Oncol435865902013

10.3892/ijo.2013.1971

23727790

Hwang

Liang

Chien

Overexpression and elevated serum levels of phosphoglycerate kinase 1 in pancreatic ductal adenocarcinoma

Proteomics6225922722006

10.1002/pmic.200500345

16493704

Zieker

Königsrainer

Tritschler

Löffler

Beckert

Traub

Nieselt

Bühler

Weller

Gaedcke

Phosphoglycerate kinase 1 a promoting enzyme for peritoneal dissemination in gastric cancer

Int J Cancer126151315202010

19688824

Rojas-Pirela

Andrade-Alviárez

Rojas

Kemmerling

Cáceres

Michels

Concepción

Quiñones

Phosphoglycerate kinase: Structural aspects and functions, with special emphasis on the enzyme from Kinetoplastea

Open Biol102003022020

10.1098/rsob.200302

33234025

Lowy

Clements

Bishop

Kong

Bonney

Sisco

Aronow

Fenoglio-Preiser

Groden

β-Catenin/Wnt signaling regulates expression of the membrane type 3 matrix metalloproteinase in gastric cancer

Cancer Res66473447412006

10.1158/0008-5472.CAN-05-4268

16651426

Yamada

Takaoka

Naishiro

Hayashi

Maruyama

Maesawa

Ochiai

Hirohashi

Transactivation of the multidrug resistance 1 gene by T-cell factor 4/beta-catenin complex in early colorectal carcinogenesis

Cancer Res60476147662000

10987283

Jiang

Meisenhelder

Yang

Hawke

Zheng

Xia

Aldape

Hunter

Mitochondria-translocated PGK1 functions as a protein kinase to coordinate glycolysis and the TCA cycle in tumorigenesis

Mol Cell617057192016

10.1016/j.molcel.2016.02.009

26942675

Non-metabolic functions of glycolytic enzymes in tumorigenesis

Oncogene36262926362017

10.1038/onc.2016.410

27797379

Kalejta

Shenk

Proteasome-dependent, ubiquitin-independent degradation of the Rb family of tumor suppressors by the human cytomegalovirus pp71 protein

Proc Natl Acad Sci USA100326332682003

10.1073/pnas.0538058100

12626766

Schwartz

Ciechanover

Targeting proteins for destruction by the ubiquitin system: Implications for human pathobiology

Annu Rev Pharmacol Toxicol4973962009

10.1146/annurev.pharmtox.051208.165340

18834306

Đukić

Lulić

Thomas

Skelin

Bennett Saidu

Grce

Banks

Tomaić

HPV oncoproteins and the ubiquitin proteasome system: A signature of malignancy?

Pathogens91332020

10.3390/pathogens9020133

32085533

Gupta

Singh

Varshney

Khan

Delineating crosstalk mechanisms of the ubiquitin proteasome system that regulate apoptosis

Front Cell Dev Biol6112018

10.3389/fcell.2018.00011

29479529

Taylor

Skeate

Kast

Annexin A2 in virus infection

Front Microbiol929542018

10.3389/fmicb.2018.02954

30568638

Xia

Braunstein

Toomey

Zhong

Rao

S100 proteins as an important regulator of macrophage inflammation

Front Immunol819082018

10.3389/fimmu.2017.01908

29379499

McLaughlin-Drubin

Huh

Munger

Human papillomavirus type 16 E7 oncoprotein associates with E2F6

J Virol82869587052008

10.1128/JVI.00579-08

18579589

McLaughlin-Drubin

Meyers

Munger

Cancer associated human papillomaviruses

Curr Opin Virol24594662012

10.1016/j.coviro.2012.05.004

22658985

Aarthy

Singh

Interpretations on the interaction between protein tyrosine phosphatase and E7 oncoproteins of high and low-risk HPV: A computational perception

ACS Omega616472164872021

10.1021/acsomega.1c01619

34235319

Figure 1.

Heat map presentation of the relative protein expression levels (log₂FC) for 62 proteins compared with HaCaT cells. The columns represent the HPV-transformed and -transfected cell lines, and rows indicate the proteins. The color-coded intensities represent the expression level of each protein. The upregulated and downregulated proteins are presented in blue and red, respectively. FC, fold change; HPV, human papillomavirus.

Figure 2.

Volcano plot of the 62 differentially expressed proteins in human papillomavirus-transfected and -transformed cells vs. the HaCaT control. The y-axis presents the mean expression value [-log(P-value)] and the x-axis presents the FC difference. Red dots represent significantly differentially expressed proteins (P<0.05). Grey dots represent differentially expressed proteins that did not reach statistical significance (P>0.05). FC, fold change.

Figure 3.

UpSet intersection plot based on protein linkages between HaCaT (control), HPV-transfected (HaCaT-HPV16 E7 and HaCaT-HPV18 E7) and HPV-transformed (CaSki and HeLa) cells. Horizontal bars represent the number of protein linkages identified in each cell line. Vertical bars represent the number of protein linkages present between the cell lines. Dots indicate cell lines in which the proteins were identified. HPV human papillomavirus.

Figure 4.

Chromatograms of the proteins corresponding to the liquid chromatography-tandem mass spectrometry analysis of the HPV-transfected (HaCaT-HPV16 E7 and HaCaT-HPV18 E7) and HPV-transformed (CaSki and HeLa) cells. The x-axis represents time (min) while the y-axis represents intensity count. ALDOA, fructose-bisphosphate aldolase; HSPA8, heat shock 70kDa protein 8 isoform 2 variant; HPV human papillomavirus.

Figure 5.

Enrichment analysis delineates cancer hallmarks associated with human papillomavirus transfection in HaCaT-16/18 E7. (A) Cancer hallmark enrichment pathways based on upregulated proteins in human papillomavirus-transfected HaCaT-16/18 E7 cells. Dots represent the protein count in each enrichment pathway. The color represents the logP-value of the most significant proteins in the enrichment pathways. (B) Top four cancer hallmark enrichment pathways and their associated proteins. Blue indicates proteins involved in the enrichment pathways and grey indicates proteins not involved in the enrichment pathways. ROS, reactive oxygen species.

Figure 6.

Protein-protein interaction network of upregulated proteins in human papillomavirus-transfected cells. FC, fold change.

Figure 7.

Gene Ontology analysis of (A) upregulated and (B) downregulated differentially expressed proteins in biological process. Gene Ontology analysis of (C) upregulated and (D) downregulated differentially expressed proteins in molecular function. Gene Ontology analysis of (E) upregulated and (F) downregulated differentially expressed proteins in cellular component.