www.nature.com/scientificreports OPEN Cdc6 contributes to abrogating the G1 checkpoint under hypoxic conditions in HPV E7 expressing Received: 12 January 2017 cells Accepted: 21 April 2017 Published: xx xx xxxx 1 2 2 3,4 1 1 Hanxiang Chen , Qishu Zhang , Lijun Qiao , Xueli Fan , Weifang Zhang , Weiming Zhao & Jason J. Chen The human papillomavirus (HPV) plays a central role in cervical carcinogenesis and its oncogene E7 is essential in this process. We showed here that E7 abrogated the G1 cell cycle checkpoint under hypoxia and analyzed key cell cycle related proteins for their potential role in this process. To further explore the mechanism by which E7 bypasses hypoxia-induced G1 arrest, we applied a proteomic approach and used mass spectrometry to search for proteins that are differentially expressed in E7 expressing cells under hypoxia. Among differentially expressed proteins identified, Cdc6 is a DNA replication initiation factor and exhibits oncogenic activities when overexpressed. We have recently demonstrated that Cdc6 was required for E7-induced re-replication. Significantly, here we showed that Cdc6 played a role in E7-mediated G1 checkpoint abrogation under hypoxic condition, and the function could possibly be independent from its role in DNA replication initiation. This study uncovered a new function of Cdc6 in regulating cell cycle progression and has important implications in HPV-associated cancers. Human papillomaviruses (HPVs) are double-strand, non-enveloped small DNA viruses . HPV is one of the most 2 3, common sexually transmitted infections worldwide . To date, over 170 genotypes of HPV have been identified and can be classified into two major groups: cutaneous and mucosal HPV. Infection by HPV may lead to the formation of warts, benign lesions, cervical and several other cancers. According to the clinical prognosis of the lesions they cause, mucosal (genital) HPV types can be categorized as either “high-risk” or “low-risk” types. Up to 99% of cervical cancers contain high-risk HPV . In addition, HPV has been detected in over 80% of oropharyn- geal cancers . HPV infects the basal layer of cervical epithelium and then relies on the differentiation of the host cell to complete its life cycle. HPV encodes proteins that promote S-phase re-entry in differentiating keratinocytes . Hence, HPV can manipulate the cell cycle by establishing a milieu in the differentiated keratinocytes supportive for viral DNA amplification. Some of these cell cycle alteration activities may be correlated with HPV-associated carcinogenesis. The E6 oncoprotein leads to the rapid ubiquitination and degradation of p53 while E7 binds 9 10, 11 and promotes the degradation of pRb, leading to the release of E2F and uncontrolled cell proliferation . pRb-independent functions of E7 have also been demonstrated . Under normal conditions, DNA damage arrests cells in G1 phase and prevents cells with damaged DNA from multiplying, and allowing the cellular repair sys- tems to fix damaged DNA. E7-expressing cells bypass the G1 arrest induced by DNA damage . The mechanism by which E7 regulates G1 checkpoint has been under extensive study yet is still not fully understood. We have 14, 15 recently shown that Cdk1 and WDHD1 play a key role in G1/S transition in E7-expressing cells . Cell division cycle 6 (Cdc6) is an essential regulator of DNA replication in eukaryotic cells. The well-established function of Cdc6 is to assemble prereplicative complexes (preRCs) at origins of replication during G1 phase . As a key factor for origin licensing, Cdc6 is responsible for the loading of MCM onto the origins of replication and is Department of Pathogenic Biology and Key Laboratory of Infection and Immunity of Shandong Province, Shandong University School of Basic Medical Sciences, Jinan, Shandong, 250012, China. The Cancer Research Center, Shandong University School of Medicine, Jinan, Shandong, 250012, China. Department of Medicine, University of Massachusetts Medical School, Worcester, MA, 01532, USA. Department of Dermatology, Xijing Hospital, Fourth Military Medical University, Xian, Shanxi, 710032, China. Correspondence and requests for materials should be addressed to W.Z. (email: email@example.com) or J.J.C. (email: firstname.lastname@example.org) Scientific Repo R ts | 7: 2927 | DOI:10.1038/s41598-017-03060-w 1 www.nature.com/scientificreports/ essential for the initiation of DNA replication . In G1/S transition, Cdc6 promotes cell cycle progression by acti- 18, 19 vating Cdk2, which is bounded by p21 or p27, in an ATP dependent way . Cdc6 knockdown leads to cell cycle arrest and induces apoptosis . Cdc6 is prone to being overexpressed in most cancer cells because of dysfunction in the pRb-E2F transcriptional pathway . Deregulation of Cdc6 led to the inactivation of the INK4/ARF locus, 22, 23 which encodes three important tumor suppressor genes, p16INK4a, p15INK4b, and p53 activator ARF . Cdc6 has been identified as a biological marker for cervical cancer in early detection . We have recently shown that Cdc6 is up-regulated in E7-expressing cells and plays an important role in E7-mediated re-replication . The microenvironment of a solid tumor is characterized by irregular vascularization, poor nutrient and oxygen supply. The continuously increasing cell number and the demand of O exacerbate the hypoxic stress. Hypoxia inducible factor 1α (HIF-1α) is a central molecule involved in mediating these effects in cancer cells. Of note, in general, human cancers express high levels of HIF-1α not only due to the hypoxic tumor microen- vironment, but also because of the dysregulated signaling pathway for catering and adapting the challenging circumstances. As a transcription factor, HIF-1α regulates multiple genes that involved in energy metabolism, angiogenesis and apoptosis. HIF-1α arrest cell cycle at G1 phase by up-regulating the expression of Cdk inhib- 28, 29 itors p21 or p27 under hypoxia . A non-transcriptional mechanism of HIF-1α arrest of cell cycle was also reported . In cervical cancer, HPV E7 increases HIF-1α mediated transcription by inhibiting the binding of his- tone deacetylases , leading to HIF-1α accumulation and VEGF expression, which may contribute to enhanced 32, 33 angiogenesis . Glioma cells expressing HPV-16 E7 showed a G2/M arrest with concomitant decrease in G1 and S phases subject to hypoxia . The cell cycle profiles in other types of cells expressing HPV E7 under hypoxia remain to be determined. In this study, we demonstrated that E7 abrogated the hypoxia-induced G1 arrest. We then took a proteomic approach to search for proteins that are differentially expressed in E7 expressing cells under hypoxia. Cdc6 was found to be up-regulated in E7 expressing cells under hypoxia. Signic fi antly, we demonstrated that Cdc6 played a role in E7-mediated G1 checkpoint abrogation under hypoxic condition. This study suggested a new function of Cdc6 in regulating cell cycle progression and has important implications in HPV-associated cancers. Materials and Methods Cell culture and reagents. Human retinal pigment epithelium cell line (RPE1) was maintained in a 1:1 dilution of DMEM-Ham’s F-12 medium plus 10% FBS. RPE1 cells stably expressing pBabe or HPV-16 E7 were established by retrovirus-mediated infection using the pBabe-puro-based retroviral construct. Cells were selected with 10.5 μg/ml puromycin for 3 to 6 days. Aer t ft he infected cells were pooled and expanded, they were main- tained in puromycin (6.5 μg/ml) and used within 15 passages. The HPV-16 positive human cervical cancer cell line CaSki was purchased from American Type Culture Collection (ATCC). CaSki cells were cultured in DMEM supplemented with 10% FBS and antibiotics. Deferoxamine (DFO, Sigma) and CoCl (Sigma) were used to mimic hypoxia at a dose of 200 μM and 500 μM, respectively. Cell viability assay. RPE1 vector and E7 cells were seeded in 96-well plates at a density of 7,000 cells per well. After being treated with DFO (from 50 μM to 300 μM) or CoCl (from 0.1 mM to 1 mM) for 72 hours, the cell viability was determined by the Cell Counting Kit-8 (Dojindo) according to the manufacture and read by the Infinite M200 PRO microplate reader (Tecan Group Ltd.) at 450 nm. Flow cytometry. For cell cycle analysis, cells were collected at various time points, fixed with 70% ethanol overnight at 4 °C, resuspended with 50 μg/ml RNaseA plus 50 μg/ml propidium iodide (PI), and then analyzed by fluorescence-activated cell sorting (FACS, Becton Dickinson). Cell cycle data analysis was done using Modfit LT (version 3.10). Each histogram is consisted of the Data Line, Model Fit Line and the model components. The beginning of the S-phase Rectangle is the mean position of the G1 Gaussian and the end is at the mean position 35, 36 of the G2/M Gaussian . For the bromodeoxyuridine (BrdU) labeling experiment, we used a FITC BrdU Flow Kit (BD Biosciences, San Diego, CA) and followed the manufacturer’s instruction. The BrdU data was analyzed by using FCS Express. Immunoblotting. Total cell lysates were prepared in RIPA lysis buffer (Santa Cruz). The protein concentra- tion was measured by the bicinchoninic acid (BCA) protein assay reagent (Invitrogen). Equal amounts of protein from each cell lysate were separated in an SDS polyacrylamide gel (PAGE) and transferred onto a polyvinylidene difluoride (PVDF) membrane. Membranes were reacted with antibodies against pRb (BD Biosciences), HIF-1α (BD Biosciences), Cdc6 (Santa Cruz), Cdk1 (Santa Cruz), Cdk2 (Santa Cruz), p53 (BD Biosciences), p21 (BD Biosciences) and p27 (Santa Cruz). HRP conjugated goat anti-mouse or anti-rabbit were used as secondary anti- bodies. GAPDH (Santa Cruz) or β-tubulin (Sigma) was used to indicate loading amount of total proteins. ImageJ (NIH) was used to quantify gel images. RNA interference (RNAi). Cells were seeded onto a 60 mm dish the day before transfection. Cells were transfected with a final concentration of 20 nM siRNA per target gene using Lipofectamine 2000 transfection reagents according to the manufacturer (Invitrogen). For gene knockdown analysis, cells were harvested 72 hours post-transfection and specific protein levels were analyzed by immunoblot- ting. For cell cycle analysis, twenty-four hours after transfection, cells were treated with DFO (200 μM) for an additional 8 hours. All of the siRNAs were synthesized from Sigma and the oligo-nucleotide sequences were as follows: siCdc6 sense strand: 5 ′ -CUUC C CA C CUUAUA C CA GA dTdT-3 ′ ; siRNA209 (siE6) sense strand: 5′-UCCAUAUGCUGUAUGUGAUdTdT-3′; siRNA198 (siE6E7) sense strand: 5′ -GCACACAC GUAGACAUUC G dTdT-3′. Negative control siRNA (siCon) sense strand: 5′-UUCUCCGAACGUGUCACGUdTdT-3′. Scientific Repo R ts | 7: 2927 | DOI:10.1038/s41598-017-03060-w 2 www.nature.com/scientificreports/ Reverse-transcription PCR and real-time PCR assay. Total RNA was isolated using TRIzol reagent (Invitrogen) and then transcribed to cDNA using PrimeScriptTM RT Reagent KIT with gDNA Eraser (TaKaRa). Real-time quantitative PCR was performed by the Bio-Rad CFX96 Touch Real-Time PCR Detection system. The levels of mRNA were measured by SYBR Green assay using SYBR Premix Ex TaqTM (TaKaRa) according to the manufacturer’s instructions and GAPDH was used as a control house- keeping gene. Melting curve analysis was performed to confirm amplification of specific transcripts. Each reaction was run parallel and in triplicate. The expression levels of transcripts were calculated by the rela- tive quantification (2−ΔΔCt) study method. All primer sequences were listed as follows: Cdc6 (forward, 5′-ACCTATGCAACACTCCCCATT-3′; reverse, 5′-TGGCTAGTTCTCTTTTGCTAGGA-3′), HIF-1α (forward, 5′-CACCACAGGACAGTACAGGAT-3′; reverse, 5′-CGTGCTGAATAATACCACTCACA-3′), Cdk1 (forward, 5′-AAACTACAGGTCAAGTGGTAGCC-3′; reverse, 5′-TCCTGCATAAGCACATCCTGA-3′), Cdk2 (forward, 5′-CCAGGAGTTACTTCTATGCCTGA-3′; reverse, 5′-TTCATCCAGGGGAGGTACAAC-3′), GAPDH (for- ward, 5′-GCACCGTCAAGGCTGAGAAC-3′; reverse, 5′-TGGTGAAGACGCCAGTGGA-3′), HPV-16 E6 (for- ward, 5′-ACAAACCGTTGTGTGATTTGTT-3′; reverse, 5′-CAGTGGCTTTTGACAGTTAATACA-3′), HPV-16 E7 (forward, 5′-GAACCGGACAGAGCCCATTA-3′; reverse, 5′-ACACTTGCAACAAAAGGTTACA-3′). Mass spectrometry. Sample preparation was performed according to filter aided sample preparation (FASP) method . The prepared digests were analyzed by LC-MS/MS using a Thermo EASY-nLC II system (Thermo Scientific, Denmark), which was interfaced to an LTQ Orbitrap Elite/Velos Pro mass spectrometer (Thermo Scientific). LC MS/MS files were exported from XCalibur. Protein identification was subjected to strict data QC filtering within MASCOT and SQUEST. Label-free quantification of identified proteins was performed by Progenesis LC-MS software package from Nonlinear Dynamics. Functional classification and pathway analysis of differentially expressed proteins (cut-off point 1.75) were analyzed by gene ontology (GO) via DAVID Bioinformatics Resources 6.7. For GO analyses, at least two genes were needed to be present in a cluster, and a probability (P) value below 0.05 was considered significant. Statistical analysis. All data are shown as means and standard error of the mean (SEM). The Student t test was used to compare the differences between means. Significance was set at a P value of < 0.05. Results HPV-16 E7 abrogates G1 checkpoint under hypoxia in epithelial cells. Hypoxia has been shown to inhibit cell proliferation in a wide variety of cell types while cancer cells could proliferate to some extent under 38–41 hypoxic condition . Here we assessed the effect of HPV-16 E7 on cell cycle progression under hypoxia. RPE1 cells stably expressing HPV-16 E7 were established . Cell viability was analyzed aer ft treating cells with DFO or CoCl to mimic hypoxic conditions. As shown in Fig. 1a, the cell viability of control cells was attenuated with increasing drug concentration while E7 expressing cells had a significantly higher rate. To examine cell cycle profiles of E7 expressing cells under hypoxia, cells were cultured in hypoxic chamber at 1% O , or treated with DFO and then analyzed by flow cytometry. In hypoxic chamber, 74.5% of vector cells accumulated in G1 phase while less (65.1%) E7 expressing cells arrested at G1 (Fig. 1b). Similar observations were found aer D ft FO treat- ment (72.5% vector cells and 61.2% E7 expressing cells at G1) (Fig. 1b), suggesting that hypoxia induced cell cycle arrest in G1 phase was partly abrogated by HPV-16 E7. To examine the effect of E7 in promoting S-phase entry of cells under hypoxia more directly, we measured BrdU incorporation aer t ft reating E7 expressing cells with DFO. Significantly, much more BrdU incorporation (30.3% vs. 16.1%) was observed in RPE1 E7 cells compared with vector control cells (Fig. 1c). Similar results were observed when cells were treated with CoCl (data not shown). es Th e results demonstrate that E7 facilitates S-phase entry under hypoxia. Hypoxia leads to G1 arrest in HPV-16 E7 knockdown cervical cancer cells. Given the association between HPV-16 E7 and cell cycle progression under hypoxia, we examined the effect by using RNAi technology in HPV-16 positive cervical carcinoma cells. CaSki cells are derived from a cervical cancer patient and contain around 600 integrated HPV-16 genome copies. In high-risk HPV associated cervical cancer cells, E6*I is an extremely efficient spliced transcript and the mRNAs of E6* I are responsible for HPV-16 E7 production . The siRNA targeting E7 potentially impact the expression of E6. Given by this, in order to examine the function of E7 on hypoxia-induced G1 arrest, we utilized siRNA198 and siRNA209 to selectively knockdown HPV-16 E6E7 and E6, respectively . The expressions of E6 and E7 targets, p53 and pRb, were used to evaluate the knockdown efficiency. As shown in Fig. 2a, aer E6 k ft nocked down, the expression of p53 was increased in both CaSki-siE6E7 and CaSki-siE6 cells. Compared with siE6, the expression of pRb was increased in CaSki-siE6E7 cells, indicating the knockdown of E7. We also assessed the expression of E6 and E7 at mRNA level (Fig. 2b). Knockdown E6E7 decreased both E6 and E7 mRNA expressions while knockdown E6 specifically impaired the expression of E6. es Th e results suggest that CaSki-siE6E7 cells could be used as control for CaSki-siE6 cells in studying the function of E7. Notably, E6 knockdown led to increased p53 protein expression, which is consistent with previous obser- vations of p53 upregulation in E7 expressing cells . We next determined the cell cycle profiles of CaSki-siE6E7 and CaSki-siE6 cells under hypoxia. Both E6E7 and E6 knockdown in CaSki cells led to an increased number of cells at the G1 phase (Fig. 2c). Aer b ft eing treated with DFO, 65.1% of CaSki-siE6 cells stayed in G1 phase while 80.3% of CaSki-siE6E7 cells accumulated at G1 (Fig. 2c), indicating that E7 promotes G1/S transition under hypoxia in cervical cancer cells. To demonstrate the role of E7 in promoting S-phase entry of cervical cancer cells more directly, we measured BrdU incorporation upon DFO treatment. We found that knockdown of E7 by siRNAs led to a significant reduc- tion in BrdU incorporation (from 7.8% for CaSki-siE6 to 4.1% for CaSki-siE6E7) ae ft r DFO treatment (Fig. 2d). Scientific Repo R ts | 7: 2927 | DOI:10.1038/s41598-017-03060-w 3 www.nature.com/scientificreports/ Figure 1. HPV-16 E7 abrogates hypoxia-induced G1 arrest. (a) RPE1 E7 cells were treated with hypoxia mimic drug DFO or CoCl at different concentrations for 72 hours, cell viability was determined using CCK8 assay. Data from a representative of 3 experiments are shown. (b) RPE1 vector and RPE1 E7 cells were incubated in hypoxic chamber at 1% oxygen or with 200 μM DFO for 8 hours, respectively. Cells were stained by PI and examined by flow cytometry. Data from a representative of 4 experiments are shown (Upper panel) and are summarized (Lower panel). (c) RPE1 vector and E7 expressing cells were treated with 200 μM DFO for 6 hours and labeled with 20 nM BrdU for 2 additional hours. Cells were stained with anti-BrdU antibody, counterstained with 7-AAD, analyzed by flow cytometry. Data from a representative of 4 experiments are shown (Upper panel) and the percentage of BrdU-positive cells was gated. Data are summarized (Lower panel). Error bars reflected the standard error of the mean. *P < 0.05; **P < 0.01; ***P < 0.001. es Th e results demonstrate an important role of E7 in the G1 cell cycle control and S-phase entry of cervical cancer cells under hypoxia. Level of cell cycle-related proteins in E7 expressing cells under hypoxic condition. To under- stand the mechanism by which HPV-16 E7 abrogates hypoxia-induced G1 arrest, we examined the expressions of cell cycle-related proteins under hypoxic condition. As shown in Fig. 3a, the steady-state level of hypoxia marker HIF-1α was higher in E7 expressing cells compared with the vector control cells and increased aer D ft FO treatment. These results are consistent with a previous report . Both Cdk1 and Cdk2 levels were decreased under hypoxia while their expression levels were still higher in E7 expressing cells until 8 hours post treatment. No significant differences were found for the expressions of Cdk4, Cdk6 and cyclin A aer D ft FO treatment. The level of p53 remained high in DFO-treated E7 expressing cells. In contrast, although p53 expression was low in the vector control cells, it was elevated ae ft r DFO treatment. Both p27 and p21 expressions showed similar patterns until the treatment time reached 16 hours, where they dropped (Fig. 3a). We also examined the expressions of a few selected cell cycle related proteins ae ft r cells were cultured in low oxygen chamber. As shown in Fig. 3b, the steady-state level of HIF-1α was higher in E7 expressing cells and increased aer ce ft lls were cultured with 1% O for 8 hours. The expressions of p53 and p21 did not change much in E7 expressing cells but increased in vector cells aer ft culturing in low oxygen. Both Cdk1 and Cdk2 expressions were decreased aer ft cells were cultured in hypoxic chamber. Next we determined the mRNA levels of HIF-1α, Cdk1 and Cdk2. As shown in Fig. 3c, compared with the vector control cells, the level of HIF-1α was higher in regularly cultured E7 expressing cells. Upon DFO treat- ment, the level of HIF-1α mRNA in E7 expressing cells went up at 8 hours post DFO treatment and dropped at 16 hours post treatment. In the vector control cells, however, the level of HIF-1α mRNA decreased after DFO treatment. A similar pattern of mRNA expression was observed for Cdk2 in both vector and E7 expressing cells (Fig. 3c). The level of Cdk1 mRNA was also decreased after DFO treatment. Nonetheless, it was consistently higher in E7 expressing cells than the vector control cells. In summary, under hypoxic condition, the mRNA levels of HIF-1α, Cdk1 and Cdk2 were consistently higher in E7 expressing cells, suggesting that Cdk1 and Cdk2 may counteract the negative effect of HIF-1α on cell cycle. Scientific Repo R ts | 7: 2927 | DOI:10.1038/s41598-017-03060-w 4 www.nature.com/scientificreports/ Figure 2. Hypoxia leads to G1 arrest in HPV-16 E7 knockdown cervical cancer cells. HPV-16 positive CaSki cells were transfected with siRNAs targeting E7 or E6 (siRNA198 or siRNA209). (a) Forty-eight hours later, the steady-state levels of p53 and pRb were examined by Immunoblotting (full-length blots are presented in Supplementary Figure 1). (b) Thirty-six hours post-transfection, the mRNA levels of HPV-16 E6 and HPV-16 E7 were determined by RT-qPCR. (c) Twenty-four hours post-transfection, cells were treated with 200 μM DFO for 8 hours. Cells were stained with PI and analyzed by flow cytometry. A representative experiment of 3 is shown (Upper panel) and summarized (Lower panel). (d) Twenty-four hours aer t ft ransfection, cells were treated with 200 μM DFO for 6 hours and then labeled with BrdU for 2 additional hours. Cells were stained with anti-BrdU antibody, counterstained with 7-AAD, analyzed by flow cytometry. Data from a representative of 3 experiments are shown (Upper panel) and the percentage of BrdU-positive cells was gated. Mean percentage of BrdU positive CaSki cells are summarized (Lower panel). Error bars reflected the standard error of the mean. *P < 0.05; **P < 0.01; ***P < 0.001. Proteomic analysis of E7 abrogation of G1 checkpoint under hypoxia. To systematically explore the mechanism by which HPV-16 E7 abrogates hypoxia-induced G1 arrest, we took a proteomic approach and performed label-free mass spectrometry to search for proteins that are differentially expressed between E7 expressing cells and vector control cells under hypoxia. One hundred and eighty-two differentially expressed pro- teins were identified at cut-off point 1.75 fold (Supplementary Table 1). Among these proteins, many of them are 46 47 implicated in E7 associated carcinogenesis, such as PAK2 and SPARC . We then conducted functional classifi- cation and pathway analysis. Pathways with altered protein member expression include translation, cell motion, protein complex assembly, DNA replication and oxidation reduction (Fig. 4a). These pathways are consistent with known biological activities of E7, which has been reported to regulate multiple genes to modulate translation , Scientific Repo R ts | 7: 2927 | DOI:10.1038/s41598-017-03060-w 5 www.nature.com/scientificreports/ Figure 3. Expression of cell cycle-related proteins under hypoxic conditions. The steady-state levels of HIF-1α, p53, p21, p27, Cdk1, Cdk2, Cdk4, Cdk6 and cyclinA in RPE1 cells treated with 200 μM DFO (a) or 1% O chamber (b) were examined by immunoblotting (Left panels) (original blots are presented in Supplementary Figure 2) and quantified (Right panels). A representative of 3 independent experiments is shown. (c) The mRNA level of HIF-1α, Cdk1 and Cdk2 were detected by RT-qPCR. GAPDH was used as a control. *P < 0.05; **P < 0.01; ***P < 0.001. 49 50 51 promote cell migration , play a role in macromolecular complex subunit organization , DNA replication and redox balancing . Among the differentially expressed proteins, Cdc6 draws our attention. We have recently demonstrated that Cdc6 is upregulated in HPV E7 expressing cells and plays an important role in E7-induced re-replication . In addition, it has been shown that Cdc6 activates p21 or p27 bound Cdk2-cyclin A/E complexes that may poten- tially contribute to cell cycle progression. Our M.S. study showed that the protein level of Cdc6 in E7 expressing cells was 1.79-fold higher than in the vector control cells under hypoxic condition (Fig. 4b). We then confirmed the protein expression of Cdc6 by Immunoblotting aer ce ft lls were treated in hypoxic chamber (Fig. 4c). In nor- moxia, consistent with what we observed , Cdc6 expression was higher in E7 expressing cells. Under hypoxia, Cdc6 protein level decreased in both vector and E7 expressing cells. However, the steady-state level of Cdc6 was still higher in E7 expressing cells. The mRNA level of Cdc6 was also higher in E7 expressing cells (Fig. 4d). Scientific Repo R ts | 7: 2927 | DOI:10.1038/s41598-017-03060-w 6 www.nature.com/scientificreports/ Figure 4. Label-free mass spectrometry quantification of RPE1 vector and E7 expressing cells under hypoxia. (a) RPE1 vector and E7 expressing cells treated with 200 μM DFO were prepared for mass spectrometry. Functional classification of differentially expressed proteins that were quantified by M.S. (b) Expression level of Cdc6 in vector and E7 expressing cells under hypoxia was determined by M.S. (c) Verification of Cdc6 steady-state level in vector and E7 expressing cells treated in hypoxic chamber at 1% oxygen by immunoblotting (original blots are presented in Supplementary Figure 3). (d) Verification of Cdc6 mRNA level in vector and E7 expressing cells treated with 200 μM DFO by RT-qPCR. Error bars reflected the standard error of the mean. *P < 0.05. Interestingly, in response to hypoxia, the level of Cdc6 mRNA was increased and higher in E7 expressing cells, suggesting that the Cdc6 protein was degraded more rapidly in E7 expressing cells. Cdc6 is important for E7 cells to bypass G1 checkpoint under hypoxia. Since we have observed that Cdc6 was one of the noticeable genes up-regulated in response to hypoxia, we decided to examine its contribution to E7-mediated G1 checkpoint abrogation in response to hypoxia. We took an RNAi approach to knockdown Cdc6 in the cells. At a final concentration of 3 nM siRNA, Cdc6 was reduced to a level close to what was observed in the vector control cells (Fig. 5a). Under such a condition, the normal DNA replication function of Cdc6 should be intact. Accordingly, we knocked down Cdc6 and the cell cycle profile of E7 expressing cells under hypoxia was exam- ined. Notably, the effect of HPV 16 E7 on abrogating G1 checkpoint under hypoxia was abolished aer k ft nocking down Cdc6. Significantly, more cells were at G1 aer ft Cdc6 knockdown (63.9% for negative control siRNA while 74.2% for Cdc6 siRNA) (Fig. 5b). Next we performed BrdU assay to assess S phase entry in E7 expressing cells after Cdc6 was knocked down. As shown in Fig. 5c, the number of cells in S phase was reduced from 19.1% to 5.4% aer C ft dc6 knockdown. These results indicate that Cdc6 is critical for E7 expressing cells in bypassing G1 checkpoint under hypoxic conditions. Discussion In this study, we showed that HIF-1α was up-regulated in HPV-16 E7 expressing cells. While hypoxia results in cell cycle arrest of normal cells at G1 phase of the cell cycle via HIF-1α, HPV-16 E7 expressing cells bypassed this G1 checkpoint. To explore the mechanism by which E7 abrogates hypoxia-induced G1 arrest, we used mass spec- trometry and identified nearly 200 proteins that are differentially expressed in E7 expressing cells under hypoxia. Among the proteins identified, Cdc6 has been implicated in G1 checkpoint regulation and exhibits oncogenic activities. Our studies demonstrated that Cdc6 was required for E7-mediated G1 checkpoint abrogation under hypoxic conditions. This study uncovered a new function of Cdc6 that might be involved in regulating cell cycle progression under hypoxia and has important implications in HPV-associated cancers. One important characteristic of solid tumor is poor oxygen supply. Hypoxia leads to increased production of HIF-1α that may lead to angiogenesis and reduced cell proliferation. In cervical cancer, the high expression of 31–33 HIF-1α leads to increased VEGF expression and angiogenesis . e Th cell cycle profiles in cells expressing HPV E7 under hypoxia were not clear. Glioma cells expressing HPV-16 E7 showed a concomitant decrease in G1 and S phases subject to hypoxia . Here we demonstrated that HPV-16 E7 abrogates G1 checkpoint under hypoxia and rendered cells enter into S-phase. This function of E7 may contribute to cervical carcinogenesis. In this study, we uncovered a new function of Cdc6 controlling G1 checkpoint under hypoxic conditions. We believe the elevated Cdc6 in E7 expressing cells was resulted from the deregulation of the pRb-E2F transcriptional pathway. Despite elevated levels of p53, p21 and p27 were found in E7 expressing cells, they were considered 53, 54 55 dysfunctional . It is believed that DNA replication and cell cycle progression are coupled , we speculated that the DNA replication licensing factors may play a role in cell cycle checkpoint control indirectly. Therefore, the Scientific Repo R ts | 7: 2927 | DOI:10.1038/s41598-017-03060-w 7 www.nature.com/scientificreports/ Figure 5. Cdc6 is important for E7 expressing cells to alleviate G1 arrest under hypoxia. (a) Cdc6 in RPE1 E7 cells was knocked down by siRNA at indicated doses and analyzed by immunoblotting (Upper panel) (original blots are presented in Supplementary Figure 3) and RT-qPCR (Lower panel, 3 nM siRNA was used). (b) RPE1 E7 cells were transfected with either Cdc6 or negative control siRNA. Cells were treated with 200 μM DFO 24 hours post-transfection and cultured for an additional 8 hours, stained with PI and analyzed by flow cytometry. A representative experiment of 3 is shown (Upper panel) and summarized (Lower panel). (c) RPE1 E7 cells were transfected with either Cdc6 or negative control siRNA. Twenty-four hours later, cells were treated with 200 μM DFO for 6 hours and then labeled with BrdU for 2 additional hours. Cells were stained with anti-BrdU antibody, counterstained with 7-AAD, analyzed by flow cytometry. Data from a representative of 3 experiments are shown (Upper panel). Data are summarized (Lower panel). Error bars reflected the standard error of the mean. *P < 0.05; **P < 0.01. observed G1 arrest ae ft r Cdc6 knockdown could be a result of DNA replication initiation issues instead of a G1 checkpoint response. However, the level of Cdc6 after siRNA knockdown was comparable to that of the vector control cells, suggesting that the regular DNA replication process was intact. These results suggest that G1 check- point regulation and DNA replication initiation are separate functions of Cdc6. Although our study demonstrated a role for Cdc6 in bypassing G1 checkpoint in E7 expressing cells, the mechanism remains to be established. Previous studies have demonstrated that Cdc6 activates p21-associated Cdk2 in an ATP-dependent way . Since we have reported that Cdk1 plays a critical role in DNA damage-induced G1 checkpoint abrogation by E7, we speculate that Cdc6 may work through Cdk1 in bypassing the G1 checkpoint. Future study should determine whether Cdk1 indeed contributes to G1 checkpoint abrogation in E7 expressing cells under hypoxia. In summary, HPV-16 E7 contributed to G1 checkpoint abrogation under hypoxia by up-regulating Cdc6. We demonstrated Cdc6 was a key regulator in G1/S transition, which may modulate the interaction between p21 and Cdk1. These results reveal that a new function of Cdc6 might be involved in promoting cell cycle progression under hypoxic condition in HPV-16 E7 expressing cells, and may have important implications in HPV-associated cancers. References 1. zur Hausen, H. Papillomaviruses in the causation of human cancers - a brief historical account. Virology. 384, 260–265 (2009). 2. Baseman, J. G. & Koutsky, L. A. The epidemiology of human papillomavirus infections. Journal of clinical virology: the official publication of the Pan American Society for Clinical Virology. 32(Suppl 1), S16–24 (2005). 3. Van Doorslaer, K. et al. e P Th apillomavirus Episteme: a central resource for papillomavirus sequence data and analysis. Nucleic acids research. 41, D571–578 (2013). 4. de Villiers, E. M. Cross-roads in the classification of papillomaviruses. Virology. 445, 2–10 (2013). 5. Burd, E. M. Human papillomavirus and cervical cancer. Clinical microbiology reviews. 16, 1–17 (2003). 6. Gao, G. et al. A novel RT-PCR method for quantification of human papillomavirus transcripts in archived tissues and its application in oropharyngeal cancer prognosis. International journal of cancer. 132, 882–890 (2013). 7. Banerjee, N. S. et al. Conditionally activated E7 proteins of high-risk and low-risk human papillomaviruses induce S phase in postmitotic, differentiated human keratinocytes. Journal of virology. 80, 6517–6524 (2006). 8. Chen, J. J., Reid, C. E., Band, V. & Androphy, E. J. Interaction of papillomavirus E6 oncoproteins with a putative calcium-binding protein. Science (New York, N.Y.). 269, 529–531 (1995). Scientific Repo R ts | 7: 2927 | DOI:10.1038/s41598-017-03060-w 8 www.nature.com/scientificreports/ 9. Dyson, N., Howley, P. M., Munger, K. & Harlow, E. The human papilloma virus-16 E7 oncoprotein is able to bind to the retinoblastoma gene product. Science (New York, N.Y.). 243, 934–937 (1989). 10. Johnson, D. G., Schwarz, J. K., Cress, W. D. & Nevins, J. R. Expression of transcription factor E2F1 induces quiescent cells to enter S phase. Nature. 365, 349–352 (1993). 11. Chen, J. J. Genomic Instability Induced By Human Papillomavirus Oncogenes. North American journal of medicine & science. 3, 43–47 (2010). 12. Shin, M. K., Sage, J. & Lambert, P. F. Inactivating all three rb family pocket proteins is insufficient to initiate cervical cancer. Cancer research. 72, 5418–5427 (2012). 13. Demers, G. W., Espling, E., Harry, J. B., Etscheid, B. G. & Galloway, D. A. Abrogation of growth arrest signals by human papillomavirus type 16 E7 is mediated by sequences required for transformation. Journal of virology. 70, 6862–6869 (1996). 14. Fan, X. & Chen, J. J. Role of Cdk1 in DNA damage-induced G1 checkpoint abrogation by the human papillomavirus E7 oncogene. Cell cycle (Georgetown, Tex.). 13, 3249–3259 (2014). 15. Zhou, Y. et al. Role of WDHD1 in Human Papillomavirus-Mediated Oncogenesis Identified by Transcriptional Profiling of E7- Expressing Cells. Journal of virology. 90, 6071–6084 (2016). 16. Borlado, L. R. & Mendez, J. CDC6: from DNA replication to cell cycle checkpoints and oncogenesis. Carcinogenesis. 29, 237–243 (2008). 17. Cocker, J. H., Piatti, S., Santocanale, C., Nasmyth, K. & Diffley, J. F. An essential role for the Cdc6 protein in forming the pre- replicative complexes of budding yeast. Nature. 379, 180–182 (1996). 18. Kan, Q., Jinno, S., Yamamoto, H., Kobayashi, K. & Okayama, H. ATP-dependent activation of p21WAF1/CIP1-associated Cdk2 by Cdc6. Proceedings of the National Academy of Sciences of the United States of America. 105, 4757–4762 (2008). 19. Uranbileg, B. et al. Cdc6 protein activates p27KIP1-bound Cdk2 protein only after the bound p27 protein undergoes C-terminal phosphorylation. e J Th ournal of biological chemistry. 287, 6275–6283 (2012). 20. Lau, E., Zhu, C., Abraham, R. T. & Jiang, W. The functional role of Cdc6 in S-G2/M in mammalian cells. EMBO reports. 7, 425–430 (2006). 21. Hateboer, G. et al. Cell cycle-regulated expression of mammalian CDC6 is dependent on E2F. Molecular and cellular biology. 18, 6679–6697 (1998). 22. Gonzalez, S. et al. Oncogenic activity of Cdc6 through repression of the INK4/ARF locus. Nature. 440, 702–706 (2006). 23. Kim, W. Y. & Sharpless, N. E. The regulation of INK4/ARF in cancer and aging. Cell. 127, 265–275 (2006). 24. Bonds, L., Baker, P., Gup, C. & Shroyer, K. R. Immunohistochemical localization of cdc6 in squamous and glandular neoplasia of the uterine cervix. Archives of pathology & laboratory medicine. 126, 1164–1168 (2002). 25. Fan, X., Zhou, Y. & Chen, J. J. Role of Cdc6 in re-replication in cells expressing human papillomavirus E7 oncogene. Carcinogenesis. 37, 799–809 (2016). 26. Zhong, H. et al. Overexpression of hypoxia-inducible factor 1alpha in common human cancers and their metastases. Cancer research. 59, 5830–5835 (1999). 27. Pugh, C. W. & Ratcliffe, P. J. Regulation of angiogenesis by hypoxia: role of the HIF system. Nature medicine. 9, 677–684 (2003). 28. Cho, Y. S. et al. HIF-1alpha controls keratinocyte proliferation by up-regulating p21(WAF1/Cip1). Biochimica et biophysica acta. 1783, 323–333 (2008). 29. Gardner, L. B. et al. Hypoxia inhibits G1/S transition through regulation of p27 expression. e J Th ournal of biological chemistry. 276, 7919–7926 (2001). 30. Hubbi, M. E. et al. A nontranscriptional role for HIF-1alpha as a direct inhibitor of DNA replication. Science signaling. 6, ra10 (2013). 31. Bodily, J. M., Mehta, K. P. & Laimins, L. A. Human papillomavirus E7 enhances hypoxia-inducible factor 1-mediated transcription by inhibiting binding of histone deacetylases. Cancer research. 71, 1187–1195 (2011). 32. Nakamura, M. et al. Hypoxia-specific stabilization of HIF-1alpha by human papillomaviruses. Virology. 387, 442–448 (2009). 33. Rodolico, V. et al. Hypoxia inducible factor-1 alpha expression is increased in infected positive HPV16 DNA oral squamous cell carcinoma and positively associated with HPV16 E7 oncoprotein. Infectious agents and cancer. 6, 18 (2011). 34. Moon, S. U. et al. The expression of human papillomavirus type 16 (HPV16 E7) induces cell cycle arrest and apoptosis in radiation and hypoxia resistant glioblastoma cells. Molecular medicine reports. 4, 1247–1253 (2011). 35. Bagwell, C. B. et al. Optimizing flow cytometric DNA ploidy and S-phase fraction as independent prognostic markers for node- negative breast cancer specimens. Cytometry. 46, 121–135 (2001). 36. Bagwell, C. B. et al. DNA and cell cycle analysis as prognostic indicators in breast tumors revisited. Clinics in laboratory medicine. 21, 875–895, x (2001). 37. Zhou, N., Yuan, S., Wang, R., Zhang, W. & Chen, J. J. Role of dual specificity tyrosine-phosphorylation-regulated kinase 1B (Dyrk1B) in S-phase entry of HPV E7 expressing cells from quiescence. Oncotarget. 6, 30745–30761 (2015). 38. Koshiji, M. et al. HIF-1alpha induces cell cycle arrest by functionally counteracting Myc. e EMB Th O journal. 23, 1949–1956 (2004). 39. Gordan, J. D., Bertout, J. A., Hu, C. J., Diehl, J. A. & Simon, M. C. HIF-2alpha promotes hypoxic cell proliferation by enhancing c-myc transcriptional activity. Cancer cell. 11, 335–347 (2007). 40. Goda, N. et al. Hypoxia-inducible factor 1alpha is essential for cell cycle arrest during hypoxia. Molecular and cellular biology. 23, 359–369 (2003). 41. Gordan, J. D., Thompson, C. B. & Simon, M. C. HIF and c-Myc: sibling rivals for control of cancer cell metabolism and proliferation. Cancer cell. 12, 108–113 (2007). 42. Fan, X., Liu, Y., Heilman, S. A. & Chen, J. J. Human papillomavirus E7 induces rereplication in response to DNA damage. Journal of virology. 87, 1200–1210 (2013). 43. Tang, S., Tao, M., McCoy, J. P. Jr. & Zheng, Z. M. The E7 oncoprotein is translated from spliced E6*I transcripts in high-risk human papillomavirus type 16- or type 18-positive cervical cancer cell lines via translation reinitiation. Journal of virology. 80, 4249–4263 (2006). 44. Tang, S., Tao, M., McCoy, J. P. Jr. & Zheng, Z. M. Short-term induction and long-term suppression of HPV16 oncogene silencing by RNA interference in cervical cancer cells. Oncogene. 25, 2094–2104 (2006). 45. Zhang, W. et al. The human papillomavirus type 58 E7 oncoprotein modulates cell cycle regulatory proteins and abrogates cell cycle checkpoints. Virology. 397, 139–144 (2010). 46. Wu, R., Abramson, A. L., Symons, M. H. & Steinberg, B. M. Pak1 and Pak2 are activated in recurrent respiratory papillomas, contributing to one pathway of Rac1-mediated COX-2 expression. International journal of cancer. 127, 2230–2237 (2010). 47. Yiu, G. K. et al. SPARC (secreted protein acidic and rich in cysteine) induces apoptosis in ovarian cancer cells. e Th American journal of pathology. 159, 609–622 (2001). 48. Songock, W. K., Kim, S. M. & Bodily, J. M. The human papillomavirus E7 oncoprotein as a regulator of transcription. Virus research (2016). 49. Hu, D. et al. HPV-16 E6/E7 promotes cell migration and invasion in cervical cancer via regulating cadherin switch in vitro and in vivo. Archives of gynecology and obstetrics. 292, 1345–1354 (2015). 50. Roman, A. & Munger, K. The papillomavirus E7 proteins. Virology. 445, 138–168 (2013). 51. Munger, K. et al. Mechanisms of human papillomavirus-induced oncogenesis. Journal of virology. 78, 11451–11460 (2004). 52. De Marco, F. Oxidative stress and HPV carcinogenesis. Viruses. 5, 708–731 (2013). Scientific Repo R ts | 7: 2927 | DOI:10.1038/s41598-017-03060-w 9 www.nature.com/scientificreports/ 53. Zerfass-Thome, K. et al. Inactivation of the cdk inhibitor p27KIP1 by the human papillomavirus type 16 E7 oncoprotein. Oncogene. 13, 2323–2330 (1996). 54. Massimi, P. & Banks, L. Repression of p53 transcriptional activity by the HPV E7 proteins. Virology. 227, 255–259 (1997). 55. Kelly, T. J., Nurse, P. & Forsburg, S. L. Coupling DNA replication to the cell cycle. Cold Spring Harbor symposia on quantitative biology. 58, 637–644 (1993). 56. Tada, S. Cdt1 and geminin: role during cell cycle progression and DNA damage in higher eukaryotes. Frontiers in bioscience: a journal and virtual library. 12, 1629–1641 (2007). Acknowledgements We thank members of our laboratories for technical assistance and helpful discussion. We thank Wesley Haynes for English editing. This work, including the efforts of Jason Chen, was funded by National Natural Science Foundation of China (Grant No. 81471944). Weiming Zhao was supported by the National Natural Science Foundation of China (NSFC) (81472552). Weifang Zhang was supported by the National Natural Science Foundation of China (NSFC) (81571986). Author Contributions H.C. performed experiments and drae ft d the manuscripts. Q.Z. collected and analyzed the MS data. L.Q., X.F. and W.Z. assisted in analyzing data. J.C. and W.Z. designed the experiments, reviewed and edited the manuscripts. Additional Information Supplementary information accompanies this paper at doi:10.1038/s41598-017-03060-w Competing Interests: The authors declare that they have no competing interests. Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Cre- ative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not per- mitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/. © The Author(s) 2017 Scientific Repo R ts | 7: 2927 | DOI:10.1038/s41598-017-03060-w 10
Scientific Reports – Springer Journals
Published: Jun 7, 2017
It’s your single place to instantly
discover and read the research
that matters to you.
Enjoy affordable access to
over 18 million articles from more than
15,000 peer-reviewed journals.
All for just $49/month
Query the DeepDyve database, plus search all of PubMed and Google Scholar seamlessly
Save any article or search result from DeepDyve, PubMed, and Google Scholar... all in one place.
Get unlimited, online access to over 18 million full-text articles from more than 15,000 scientific journals.
Read from thousands of the leading scholarly journals from SpringerNature, Elsevier, Wiley-Blackwell, Oxford University Press and more.
All the latest content is available, no embargo periods.
“Hi guys, I cannot tell you how much I love this resource. Incredible. I really believe you've hit the nail on the head with this site in regards to solving the research-purchase issue.”Daniel C.
“Whoa! It’s like Spotify but for academic articles.”@Phil_Robichaud
“I must say, @deepdyve is a fabulous solution to the independent researcher's problem of #access to #information.”@deepthiw
“My last article couldn't be possible without the platform @deepdyve that makes journal papers cheaper.”@JoseServera