ABSTRACT
Human papillomaviruses (HPVs) infect keratinocytes of stratified epithelia. Long-term persistence of infection is a critical risk factor for the development of HPV-induced malignancies. Through the actions of its oncogenes, HPV evades host immune responses to facilitate its productive life cycle. In this work, we discovered a previously unknown function of the HPV16 E5 oncoprotein in the suppression of interferon (IFN) responses. This suppression is focused on keratinocyte-specific IFN-κ and is mediated through E5-induced changes in growth factor signaling pathways, as identified through phosphoproteomics analysis. The loss of E5 in keratinocytes maintaining the complete HPV16 genome results in the derepression of IFNK transcription and subsequent JAK/STAT-dependent upregulation of several IFN-stimulated genes (ISGs) at both the mRNA and protein levels. We also established a link between the loss of E5 and the subsequent loss of genome maintenance and stability, resulting in increased genome integration.
IMPORTANCE Persistent human papillomavirus infections can cause a variety of significant cancers. The ability of HPV to persist depends on evasion of the host immune system. In this study, we show that the HPV16 E5 protein can suppress an important aspect of the host immune response. In addition, we find that the E5 protein is important for helping the virus avoid integration into the host genome, which is a frequent step along the pathway to cancer development.
INTRODUCTION
Human papillomavirus (HPV) infections typically cause benign lesions in cutaneous or mucosal squamous epithelia (1), but persistent infection can result in transition from a benign lesion to malignancy (2–4). The mechanisms involved in HPV persistence are largely unknown, but it is known that the virus must remain undetected by the host immune system in order to persist (1, 5, 6). This is accomplished in part by restricted expression of the major immunogenic proteins, the late capsid proteins, to the upper layers of the epithelium with low levels of immune surveillance (1, 7–9). Further, HPV does not cause viremia, cell lysis, or induction of inflammatory responses (10). In addition to life cycle organization, HPV actively suppresses both innate and adaptive immune responses. The HPV E5, E6, and E7 oncogenes mediate a wide array of anti-inflammatory and immune-evasive functions that render HPV infections and HPV-induced lesions profoundly immunosuppressive (11, 12). The importance of the immune response in HPV infection is indicated by the increased risk of HPV infection and invasive cervical cancer in HIV-infected women and other immunosuppressed patients (13–19).
HPV oncogenes contribute to persistent infection in many ways, most prominently by activating, inhibiting, or modifying host gene expression patterns, including innate immune response genes (20). The least understood HPV oncogene, E5, encodes an 83-amino acid transmembrane protein that has been shown to aid the late stages of the viral life cycle through unknown mechanisms (21, 22). E5 has been reported to regulate multiple growth factor signaling pathways. It is best known for enhancing epidermal growth factor receptor (EGFR) activation by increasing levels of EGFR at the cell surface (23–27). E5 increases c-Met levels and activation in an EGFR-dependent manner (28). In addition, overexpressed E5 has been shown to suppress transforming growth factor beta (TGF-β) signaling by reducing levels of TGF-β receptor II (TGFBR2) transcripts (29).
The interferon (IFN) response is a key component of host innate antiviral immunity. Most type I IFNs are induced through the binding of viral products to pattern recognition receptors (PRRs), leading to activation of interferon response factors (IRFs) and NF-κB to drive the synthesis of IFN molecules (30–32). After secretion, type I IFNs bind to the type I IFN receptor (IFNAR), which induces the phosphorylation and dimerization of signal transducer and activator of transcription 1 (STAT1) and STAT2, which then form a complex and translocate to the nucleus to induce transcription of IFN-stimulated genes (ISGs) (33). ISGs are biochemically diverse but generally function to limit viral replication or spread.
Keratinocytes, targets of HPV, are sometimes designated “immune sentinels” because of the wide array of innate immune sensors and cytokines that they can produce (34–37). Even under resting conditions, keratinocytes secrete an array of cytokines to sustain immune surveillance (38–41). One of these is IFN-κ, which is constitutively expressed in unstimulated keratinocytes (unlike IFN-α/β) (42). The constitutive expression of IFN-κ in resting keratinocytes indicates that IFN-κ expression is subject to fundamentally different regulation than other type I IFNs, such as IFN-α and IFN-β. IFN-κ has been shown to signal through the type I IFN receptor and to promote the upregulation of ISGs that can then promote antiviral immunity (42). To evade IFN responses, HPV oncoproteins repress both basal and agonist-induced transcription of IFNs and ISGs (43–51). IFN-κ, in particular, seems to be a target of HPV both in vivo (52) and in cell culture (11, 53, 54).
Cross talk between IFNs and growth factor signaling pathways has been previously reported. EGFR has been reported to inhibit IFN and ISG expression, including that of IFNK, through IRF1 inhibition (55–62). TGF-β has also been shown to cross talk with IFN pathway regulators (63–67). Our laboratory has recently reported that, in HPV16+ keratinocytes, TGF-β signaling can induce IFNK and ISG expression in keratinocytes harboring episomal HPV16 by reversing methylation of the IFNK promoter (65). This finding indicates that growth factor signaling is critical for regulating the levels of IFNK transcription in HPV16+ cells. As E5 is the viral oncogene most associated with the regulation of host growth factor signaling pathways, in this study, we explored the hypothesis that E5 may play a role in the regulation of IFN-κ in cells containing HPV16. We found that, although IFNK is poorly responsive to PRR agonists, levels of IFNK transcripts are suppressed in the presence of E5. We found that suppression of IFN-κ by E5 involves both the EGFR and TGFBR2 signaling pathways. Furthermore, HPV16 genomes lacking E5 tend to integrate in culture over time at a higher rate than the wild type, suggesting that suppression of IFN-κ by E5 may be required for long-term maintenance of viral genomes, revealing a novel function of E5 in the viral life cycle.
RESULTS
IFNK transcripts are upregulated in response to growth factor treatment.Since human keratinocytes are the targets of HPV, we used primary human foreskin keratinocytes (HFKs) stably maintaining episomally replicating HPV16 genomes (HPV16+ cells) in our studies (65, 68). Unlike IFN-α/β, which require stimulation in order to be expressed at appreciable levels, IFNK is constitutively expressed in resting keratinocytes (50, 53, 65), so it is an important barrier for a pathogen that preferentially targets keratinocytes. Although IFNK transcripts do increase in response to classical Toll-like receptor (TLR) agonists, such as poly(I·C), the induction is weak compared to that of PRR-responsive IFNB (data not shown). We and others have previously reported that IFNK transcription is downregulated upon differentiation (50, 65). We also previously reported, in agreement with others, that IFNK transcript and protein levels are suppressed in HPV-infected cells compared to uninfected HFKs but that this downregulation could be reversed upon treatment with TGF-β (50, 53, 65). These findings show that IFNK is not simply constitutive but can be regulated in response to cellular conditions. We sought to understand the signaling pathways that regulate IFNK in uninfected or HPV16+ cells.
We have previously shown that the growth factor TGF-β could upregulate IFNK in HPV16+ cells (65). We and others have shown that HPV can increase the levels or function of the receptor tyrosine kinases EGFR and c-Met (the hepatocyte growth factor [HGF] receptor) (24, 25, 28). Therefore, we tested whether EGF and HGF could regulate IFNK transcripts in these cells. Values were normalized to untreated controls to determine the effect of each growth factor. Treatment of HPV16+ cells with TGF-β resulted in upregulation of IFNK in cells containing HPV16, as reported previously (Fig. 1) (65). Upregulation was also observed in response to HGF (Fig. 1). EGF could induce IFNK transcripts significantly, but less markedly than TGF-β1 (Fig. 1). These observations are consistent with the premise that IFNK gene expression is not dictated by pathogen recognition, as other type I IFNs are, but rather responds to growth factors present in the epithelium.
IFNK transcripts are readily induced by growth factors. HPV16+ cells and E5 Stop cells were treated with the indicated growth factors for 24 h. Levels of IFNK mRNAs were measured using RT-qPCR. *, P < 0.05; **, P < 0.01; ***, P < 0.001. The error bars indicate standard error of the mean.
IFNK transcripts are suppressed by E5.E5 is the main HPV gene product known to regulate growth factor signaling (69). Because E5 can regulate all three of the pathways shown in Fig. 1 (TGF-β, EGF, and HGF) (23–29), we sought to determine whether E5 can affect IFNK transcription in response to these factors. In order to address the role of E5 in the HPV life cycle, we created keratinocyte-derived cell lines maintaining a mutant HPV16 genome with a stop codon in the E5 open reading frame (ORF) (E5 Stop cells) (28). These mutant viral genomes can immortalize HFKs and maintain themselves as episomes with normal levels of E6/E7 expression (28) but are genetically unable to express E5. We found that the response of IFNK mRNA to TGF-β and HGF was similar to that seen in wild-type HPV16+ cells but that the response of IFNK to EGF in these cells was significantly reduced (Fig. 1), consistent with lower levels of EGFR activity in the cells due to the absence of E5 (28).
To better understand the role of E5 in the regulation of IFNK, we compared the total levels of IFN transcripts in HFKs, HPV16+ cells, and E5 Stop cells. Because the cells were unstimulated, basal levels of IFNA and IFNB transcripts were very low (CT values of 30 to 33). IFNA levels were reduced by the presence of HPV16, and IFNB levels were unchanged; neither differed in the absence of E5 (Fig. 2a). Next, we confirmed previous data showing that HPV16 can inhibit IFNK transcription compared to HFKs. As previously reported, mRNA levels of IFNK were reduced by the presence of HPV16 (65), but this suppression was not as efficient in E5 Stop cells, indicating that E5 contributes to suppression of IFNK by HPV16 (Fig. 2b). Failure of E5 Stop cells to suppress IFN-κ was observed at the protein level, as well (Fig. 2c and d).
IFN-κ is suppressed in an E5-dependent manner. HPV16-containing and E5 Stop cells were isolated, and their RNA or protein was harvested. (a) Levels of IFNA and IFNB mRNAs were measured by RT-qPCR. (b and c) IFNK mRNA (b) and IFN-κ protein (c) levels were measured by RT-qPCR and immunoblotting, respectively. Levels of E7 were also measured by immunoblotting. (d) Quantification of immunoblots. n > 5 for all PCRs, and n > 3 for Western blots. *, P < 0.05; **, P < 0.01; ***, P < 0.001. The error bars indicate standard error of the mean.
These experiments indicate that E5 contributes to IFNK suppression in HPV16+ cells. To determine whether E5 alone is sufficient for IFNK suppression, we examined IFNK transcript levels in h-Tert-immortalized normal oral keratinocytes (NOKs) (70) that were stably transduced with retroviruses expressing HPV16 E5 alone. As E6 is the HPV oncogene previously shown to inhibit IFNK (50, 53), NOKs expressing E6 alone were included as a positive control. Comparing immortalized NOKs with primary HFKs revealed no significant differences in IFNK levels (reference 65 and data not shown), indicating that immortalization by itself does not affect IFNK expression. We observed a reduction in IFNK transcripts in cells expressing E6, consistent with previous reports that E6 can inhibit IFNK transcription (50, 53) (Fig. 3a). Importantly, IFNK transcript levels were reduced significantly in cells expressing E5, at an efficiency at least equal to that of E6 (Fig. 3a). To confirm these results, we transduced primary HFKs with retroviruses that express either E6/E7 alone or E6/E7 along with E5 expressed from the same transcript with an internal ribosome entry site (IRES). Following immortalization, IFNK mRNA levels in these cells were measured, and we observed that E5 suppresses IFNK suppression, even in the presence of E6 (Fig. 3b). These data provide evidence that E5 can suppress IFNK, a novel function of E5.
HPV16 E5 is sufficient to suppress IFN-κ expression. (a) NOKs that were previously hTert immortalized were transduced with retroviruses containing an empty vector (LXSN-Empty), a vector expressing the HPV16 E6 protein alone (LXSN-16E6), or a vector expressing the wild-type HPV16 E5 protein (LXSN-16E5). Following selection and outgrowth, IFNK levels were measured using RT-qPCR. (b) Primary HFKs were transduced with retroviruses expressing either E6/E7 alone or E6/E7 IRES E5. Following selection, the immortalized colonies were pooled, and IFNK levels were measured using RT-qPCR. ***, P < 0.001. The error bars indicate standard error of the mean.
ISG transcripts are suppressed by E5.Because IFN-κ is constitutively expressed in keratinocytes, it can drive basal expression of ISGs (65). Since IFNK transcript and protein levels are higher in E5 Stop cells than inHPV16+ cells, we tested whether ISG levels were also increased. Levels of transcripts for several ISGs, including some that have been shown to inhibit the HPV life cycle (71–73), were suppressed as expected in HPV16+ cells but increased in E5 Stop cells (Fig. 4a). This increase was also observed at the protein level for two of these factors (Fig. 4b). Both isoforms of OAS2 (p69/p71) (74) were usually observed, but differential effects of E5 were not consistently seen. To determine whether E5 alone is sufficient to suppress ISG expression, ISG transcript levels from E5-expressing NOKs were examined by reverse transcription-quantitative PCR (RT-qPCR). Significant reductions were observed in the levels of these ISG mRNAs in E5-expressing NOKs versus NOKs containing the empty vector, similar to what was seen with E6 expression alone (Fig. 4c). Similarly, expression of E5, along with E6/E7, in HFKs resulted in additional inhibition of ISGs, showing that E5 can contribute to ISG suppression even in the presence of E6 and E7 (Fig. 4d). Since IFNK is the only IFN differentially expressed in these cells (Fig. 2), we concluded that increased IFNK transcripts in E5 Stop cells were responsible for the increased levels of ISGs. Further, we concluded that IFN-κ-driven ISG expression was inhibited by E5.
Suppression of ISGs in HPV16+ cells is E5 dependent. (a) Levels of the indicated ISG mRNAs in HFK, HPV16+, and E5 Stop cells were measured using RT-qPCR. (b) Immunoblots showing the levels of OAS2 protein (top) and IFIT1 (bottom) in HFK, HPV16, and E5 Stop cells. (c) Levels of the indicated ISG mRNAs in NOKs transduced with empty or E6- or E5-expressing retroviruses were measured by RT-qPCR. (d) Levels of the indicated ISG mRNAs in HFKs immortalized with retroviruses expressing E6/E7 alone or E6/E7 IRES E5 were measured using RT-qPCR. The statistics in panels c and d are relative to HFK. **, P < 0.01; ***, P < 0.001. The error bars indicate standard error of the mean.
Increased ISG levels suggested that the IFNAR/JAK/STAT pathway could be chronically activated in E5 Stop cells. Indeed, levels of both phosphorylated and total STAT1 were higher in E5 Stop cells than in wild-type HPV16+ cells, consistent both with increased JAK signaling and with STAT1 itself being an ISG (Fig. 5a and b). The increase in STAT1 protein levels is also seen at the mRNA level, with higher STAT1 mRNA levels expressed in E5 Stop cells (Fig. 5c). Additionally, levels of STAT1 mRNA were reduced in NOKs expressing E5 (Fig. 5d). These results indicate that E5 can suppress STAT1 activation and expression. Activation of STAT1 suggests that the high levels of ISGs in E5 Stop cells are due to IFN-κ activating the classical JAK/STAT1 pathway (42). If so, treatment with the JAK1/2 inhibitor ruxolitinib should reduce ISG levels in these cells. We treated HFKs, HPV16+ cells, and E5 Stop cells with ruxolitinib for 24 h. Values were normalized to those of dimethyl sulfoxide (DMSO) treatment for each cell line to reveal the effect of drug treatment. Ruxolitinib treatment resulted in a significant reduction of MX1, IFIT1-IFI56, OAS2, and IFI16 transcript levels in all three cell types (Fig. 6a to d), suggesting that ISG levels are governed by classical JAK/STAT regulation in these cells regardless of the presence of E5. The efficacy of ruxolitinib treatment was confirmed by its ability to inhibit STAT1 phosphorylation (Y701) in response to IFN-β treatment (Fig. 6e). The data also show that chronic ISG expression in HFKs is governed by the JAK/STAT1 pathway, as would be expected downstream of IFN-κ signaling.
STAT1 levels and signaling are suppressed by E5. (a and b) Levels of IFN-κ in HPV16+ and E5 Stop cells were measured by immunoblotting (a) and quantified (b). (c and d) IFNK mRNA levels in HFK, HPV16+, and E5 Stop cells (c) or NOKs transduced with empty vector or E5 (d) were measured using RT-qPCR. *, P < 0.05; **, P < 0.01. The error bars indicate standard error of the mean.
ISG mRNA levels are regulated by JAK/STAT signaling. (a to d) HFK, HPV16+, and E5 Stop cells were treated with DMSO (control) or 10 μM ruxolitinib (Rux) for 24 h, followed by RT-qPCR for the MX1 (a), IFIT1 (b), IFI6 (c), or OAS2 (d) transcripts. (e) The efficacy of Rux treatment was confirmed by its ability to inhibit STAT1 (pY701) in HPV16 and E5 Stop cells treated with IFN-β. *, P < 0.05; **, P < 0.01; ***, P < 0.001. The error bars indicate standard error of the mean.
Phosphoproteomics indicates that E5 regulates EGFR, TGFBR2, and STAT signaling.Having shown that IFNK gene expression is regulated by growth factors and that E5 can inhibit IFNK gene expression and the resulting IFN response, we sought to understand how alteration in growth factor signaling by E5 can affect IFNK levels. In order to build a more complete picture of the various signaling pathways and interactions occurring in cells containing HPV16 with or without E5 expression, we used an antibody proteomics array (Full Moon Biosystems). The array measures the phosphorylated and total levels of various proteins involved in a variety of pathways, including growth factor signaling, cell cycle regulation, and innate immunity. Protein harvested from uninfected HFKs, episomal HPV16+ cells, and episomal E5 Stop cells was processed according to the Full Moon Biosystems protocol and analyzed using the Database for Annotation, Visualization, and Integrated Discovery (DAVID) (75, 76). Comparisons between signal intensities allowed us to determine the fold differences of various phosphorylated proteins and total proteins between the samples and the signaling pathways regulated by HPV16 in an E5-dependent manner in the context of the complete viral genome and physiological host cell type.
HPV16+ cells exhibit increased phosphorylation and total protein levels of 142 proteins that are upregulated at least 2-fold versus those in HFKs. There were only 8 proteins that were decreased in HPV16+ cells versus HFKs. Some of the factors of interest are listed in Table 1, (the full, noncurated table is available in File S1 in the supplemental material), and Table 2 lists the pathways impacted in HPV16+ cells compared to HFKs according to DAVID analysis. Factors associated with a variety of growth factor-related signaling pathways were affected by the presence of HPV16, including the HGF-Met receptor pathway and Ras signaling (28, 77–81). Phosphorylated species of MEK1 and other mitogen-activated protein kinases (MAPKs), such as Erk, as well as the increased activation of EGFR signaling, have previously been shown to be targets of the HPV16 E5 protein (23, 28, 77, 81, 82). These results are consistent with and add to studies of the HPV18 E5 protein in the context of the complete viral genome (82). Both inhibitors of κB alpha and beta (IκB-alpha and IκB-beta) were increased by HPV16. We also observed the decrease of the Rel and p100/p52 NF-κB subunits by HPV16, which is also consistent with suppression of NF-κB activity by HPV oncogenes reported in some previous studies (83–92). STAT1 Y701 phosphorylation was also reduced, further confirming the observations shown in Fig. 5.
Curated list of factors altered by the presence of HPV16 cells versus HFKs
Pathways affected in HPV16+ cells versus HFKs
We compared E5 Stop cells to HPV16+ cells, and 94 total proteins were differentially regulated, with 19 upregulated and 75 downregulated (see File S1). STAT1 Y701 phosphorylation was upregulated, consistent with our observations shown in Fig. 5 (Table 3). We also observed the increased activation of other STATs, including STAT5A and STAT6, which are associated with growth factor signaling and Th2 polarization, respectively (93, 94). Total levels of the Rel NF-κB subunit were increased (Table 3), indicating the potential for NF-κB activity to be increased in E5 Stop cells. The predominant gene ontology pathways associated with downregulated factors were growth factor signaling pathways, including ErbB, Ras, and phosphatidylinositol 3-kinase [PI3K]), MAPK, and vascular endothelial growth factor receptor (VEGFR) (Table 4). Reduced growth factor signaling in E5 Stop cells would be expected based on experiments with overexpressed E5, but these experiments confirm that E5 regulates these pathways in keratinocytes in the context of the complete viral genome.
Curated list of factors altered by the presence of E5 Stop versus HPV16 cells
Pathways affected in E5 Stop versus HPV16+ cells
To determine if the differential regulation of any factors was E5 dependent, we compared the analysis of the HPV16+ cells versus HFKs to the analysis of the E5 Stop cells versus HPV16+ cells and noted factors that were regulated in opposite directions. Table 5 lists the factors either up- or downregulated by HPV16 in an E5-dependent manner. The most notable factors regulated in an E5-dependent manner included STAT1 Y701 phosphorylation, Rel, and activation of various MEKs and MAP kinases. Phosphorylation of p73, a protein that can regulate host cell differentiation via unknown mechanisms (95, 96), was strongly upregulated in an E5-dependent manner.
Factors regulated in an E5-dependent manner
JAK/STAT signaling is not responsible for increased IFNK in E5 Stop cells.The phosphoproteomics array, coupled with our previous data, showed that the EGFR/MAPK and IFNAR/STAT pathways are deregulated in E5 Stop cells compared to wild-type HPV16+ cells. We next sought to determine the impacts of these pathways on IFNK expression. We showed previously that ISGs are dependent on IFNAR/JAK/STAT signaling (Fig. 6), but to determine whether IFNK is similarly dependent, we treated HFK, HPV16+, and E5 Stop cells with ruxolitinib for 24 h. We also treated them with IFN-β to confirm that the inhibitor was working and to determine whether IFNK could respond to other type I IFNs. The RT-qPCR results were normalized to the DMSO treatment for each cell line to reveal the effect of drug treatment. We observed that all three cell types could increase IFNK transcripts a modest 2- to 2.5-fold in response to IFN-β (Fig. 7a). Thus, the presence of the virus does not affect the essential functions of the JAK pathway with respect to IFNK levels. Ruxolitinib alone did not suppress or only weakly suppressed IFNK mRNA. There was modest but consistent downregulation in HPV16+ cells treated with ruxolitinib and IFN-β together. We do not understand the basis for this effect. Knocking down JAK1 resulted in only a modest reduction of IFNK transcripts in HPV16 cells, but not E5 Stop cells (Fig. 7b and c). Together, these data indicate that signaling through the JAK pathway is not a major player in controlling the levels of IFNK.
Basal IFNK is modestly affected by JAK/STAT signaling. (a) HFK, HPV16+, and E5 Stop cells were treated with IFN-β with or without 10 μM ruxolitinib (Rux) for 24 h. IFNK transcript levels were measured using RT-qPCR. (b) The efficacy of JAK1-targeting siRNA (siJAK1) versus nontargeting control (NTC) was tested in HPV16-containing cells by Western blotting. (c) The levels of IFN-κ transcripts in cells transfected with siJAK1 were measured by RT-qPCR. NS, not significant; *, P < 0.05; ***, P < 0.001. The error bars indicate standard error of the mean.
Impact of EGFR signaling on IFNK.We found previously that IFNK was upregulated in response to EGF in HPV16+ and E5 Stop cells (Fig. 1). As MAPK signaling pathways were reduced in E5 Stop cells compared to wild-type HPV16+ cells in our array experiments, we used inhibitors to target EGFR to determine the effect on IFNK transcripts. Cells were treated with AG1478 (a specific EGFR inhibitor) for 24 h. The RT-qPCR results were normalized to the DMSO treatment for each cell line to reveal the effect of drug treatment. In HPV16+ cells treated with AG1478, IFNK transcript levels were increased, consistent with a suppressive effect of EGFR on IFNK that was previously reported by others (Fig. 8a) (55, 97). Upregulation was also seen in HFKs treated with AG1478. Interestingly, in E5 Stop cells, no such increase in IFNK transcripts was observed (Fig. 8a), indicating that EGFR kinase activity does not suppress IFNK in E5 Stop cells. Consistent with previous findings (28), Western blotting revealed that total EGFR levels in the absence of E5 were at or below the levels seen in HFKs, in contrast to the high levels observed in HPV16+ cells (Fig. 8b). Knocking down EGFR with small interfering RNAs (siRNAs) resulted in increased IFNK transcripts in HPV16 cells, but not in E5 Stop cells (Fig. 8c and d). These data indicate that EGFR plays an important role in the ability of E5 to suppress IFNK in HPV-containing cells. It may also play a role in suppressing IFNK in uninfected cells. Treatment of HPV16+ and E5 Stop cells with the Met inhibitor SU11274 resulted in no appreciable change in IFNK transcript levels (data not shown).
IFNK mRNA levels are suppressed by EGFR kinase activity. (a) HFK, HPV16+, and E5 Stop cells were treated for 24 h with 1 μM AG1478. IFNK transcript levels were then measured using RT-qPCR. (b) Levels of EGFR in the indicated cell types were measured by Western blotting (left) and quantified (right). (c) The efficacy of EGFR-targeting siRNA (siEGFR) versus NTC was tested in HPV16-containing cells by Western blotting. (d) The levels of IFN-κ transcripts in cells transfected with siEGFR were measured by RT-qPCR. NS, not significant; *, P < 0.05; **, P < 0.01. The error bars indicate standard error of the mean.
IFN signaling in E5 Stop cells is increased in a TGF-β-dependent manner.We have previously observed that TGF-β1 treatment of HPV16 cells promotes the upregulation of IFNK transcripts (65) (Fig. 1a). We also observed alterations in both canonical and noncanonical TGF-β pathway signaling intermediates in the phosphoproteomics array, consistent with prior observations that HPV16 E5 can suppress the TGF-β signaling pathway by downregulating the TGFBR2 receptor (29). These observations suggested that suppression of the TGF-β pathway may be important for the ability of E5 to suppress IFNK. To determine if TGF-β signaling contributes to IFNK transcript levels, we treated cells with the drug SB431542, a TGFBR2-specific inhibitor. We observed significant reductions in IFNK transcript levels in HFKs, E5 Stop cells, and HPV16+ cells (Fig. 9a). A similar reduction in IFNK transcripts was seen when we knocked down TGFBR2 with siRNA (Fig. 9b and c). These findings indicate that TGF-β signaling is critical for IFNK expression in keratinocytes regardless of the presence of the virus. Analysis of E5 Stop cells showed increased phospho-SMAD3 levels (Fig. 10) and TGFBR2 protein levels (Fig. 10b) versus wild-type HPV16+ cells. These increases were consistent and significant (Fig. 10c), showing that the TGF-β signaling pathway is increased in E5 Stop cells versus wild-type HPV16+ cells. Analysis of HFKs and HPV16 and E5 Stop cells for TGFB1, TGFB2, and TGFBR2 transcript levels showed significant suppression of TGFB1 and an increase of TGFB2 in HPV16+ cells versus HFKs, as previously reported (65) (Fig. 10d). Significantly, both the suppression of TGFB1 and activation of TGFB2 by HPV were E5 dependent. Patterns of TGFBR2 mRNA levels reflected those of TGFB1 and were also E5 dependent, as previously reported (29). NOKs expressing HPV16 E5 also exhibited significantly reduced TGFBR2 transcript levels (Fig. 10e), showing that E5 is sufficient to reduce expression of the gene in the absence of the other viral oncogenes. These data show that TGF-β-induced signaling is suppressed in an E5-dependent manner in cells containing the complete HPV16 genome.
TGF-β signaling is crucial for IFN-κ in keratinocytes. (a) HFK, HPV16+, and E5 Stop cells were treated with 5 or 10 μM SB31542 (SB) for 24 h. IFNK mRNA levels were measured using RT-qPCR. (b) The efficacy of TGFBR2-targeting siRNA (siTGFBR2) versus NTC was tested in HPV16-containing cells by Western blotting. (c) The levels of IFN-κ transcripts in cells transfected with siTGFBR2 were measured by RT-qPCR. **, P < 0.01; ***, P < 0.001. The error bars indicate standard error of the mean.
TGF-β signaling is increased in E5 Stop cells. (a and b) Western blots showing the levels of SMAD3/phosphorylated SMAD3 (a) and TGFBR2 (b). (c) Quantification of results of immunoblotting experiments. (d) mRNAs for TGFB1, TGFB2, and TGFBR2 in HFK, HPV16+, and E5 Stop cells were measured by RT-qPCR. (e) Levels of TGFBR2 mRNA in hTert-immortalized NOKs stably transduced with an empty vector or a vector expressing 16E5 were measured using RT-qPCR. The blots shown are representative. *, P < 0.05; **, P < 0.01; ***, P < 0.001. The error bars indicate standard error of the mean.
E5 does not regulate IFNK promoter methylation.Since HPV16 E6 represses IFNK transcription by inducing promoter methylation and TGF-β increases IFNK transcripts, at least in part, by reversing that methylation (53, 65) (Fig. 11a), we examined whether higher levels of IFNK transcripts in E5 Stop cells were due to reduced levels of promoter methylation. Using bisulfite-sequencing analysis on DNA derived from wild-type HPV16+ cells or E5 Stop cells, we confirmed that CpGs in the IFNK promoter are mostly methylated in wild-type HPV16+ cells (53, 65). We did not observe a significant difference in IFNK promoter methylation in E5 Stop cells compared to wild-type HPV16+ cells, as measured by bisulfite sequencing (Fig. 11b). This finding suggests that higher levels of IFNK expression in E5 Stop cells are not due to demethylation of the promoter. It also suggests that processes other than methylation can regulate IFNK levels in HPV16+ cells.
E5 does not substantially affect IFNK promoter methylation. (a) Diagram of the IFNK promoter showing the locations of the TATA box and five CpGs. (b) Bisulfite-sequencing analysis of the IFNK promoter in wild-type HPV16+ cells and E5 Stop cells. Following bisulfite modification of total DNA, the IFNK promoter was isolated by PCR and inserted into a cloning vector, and individual clones were sequenced. Each row corresponds to a separate clone, and each column corresponds to one of the five CpG sites in the promoter. Five clones each from 3 donor backgrounds were analyzed in the experiment. (c) MeDIP analysis of the CDKN2B (p15) promoter (Pro) and IFNK promoter DNA from HPV16+ cells and E5 Stop cells. DNA was subjected to immunoprecipitation with 0.25 μg/μl 5-methylcytosine-specific (5-hm-C) or preimmune (PI) serum. The immunoprecipitants were subjected to qPCR using primers specific for the p15 or IFNK promoter. n = 9. **, P < 0.01. The error bars indicate standard error of the mean.
Our previous work showed that the IFNK promoter is actively demethylated through a TDG-dependent mechanism (65). Hydroxymethylcytosine (5-hmC) is an intermediate in the TET-TDG-mediated active demethylation pathway (98, 99). Bisulfite sequencing is unable to distinguish between methylated and hydroxymethylated cytosines (100), but unlike methylated C, hydroxymethylation can still permit nucleosomes to be in an open and transcriptionally active state due to the weakening of interactions between DNA and histones (101). Therefore, it is possible that the IFNK promoter is hydroxymethylated in E5 Stop cells. We used methylated DNA immunoprecipitation (MeDIP) to test this idea. Total DNA from wild-type HPV16+ cells or E5 Stop cells was immunoprecipitated using an antibody specific for 5-hmC. Following purification, the DNA was subjected to qPCR analysis to determine enrichment compared to preimmune serum. Thillainadesan et al. showed, using MeDIP, that TGF-β signaling induces 5-hmC conversion at the CDKN2A (p15) promoter during active demethylation in keratinocytes (102), so the p15 promoter was used as a positive control for the presence of hydroxymethylation. As shown in Fig. 11c, the p15 promoter showed enrichment for 5-hmC in E5 Stop cells compared to wild-type HPV16+ cells, as expected. A modest but statistically significant increase in hydroxymethylation at the IFNK promoter was also observed in E5 Stop cells compared to wild-type HPV16+ cells (Fig. 11c). Whether this modest increase is sufficient to explain the increased IFNK transcripts found in E5 Stop cells remains to be determined.
E5 Stop cells have higher levels of viral genome integration over time.Type I IFN signaling can promote episomal loss and integration of the HPV16 genome (103, 104). Because IFN signaling is chronically activated in E5 Stop cells, we predicted that long-term maintenance of viral episomes would be disrupted. Southern blots are routinely and periodically performed on HPV-containing cell lines in the laboratory to ensure that cells used in experiments contain episomal rather than integrated HPV16 genomes (Fig. 12a). Following a review and tabulation of these Southern blots, we found that the frequency of integration in E5 Stop samples was significantly higher than in wild-type HPV16 cells (Table 6). Furthermore, when we considered only the samples that contained episomal viral DNA, quantification of band intensities showed that E5 Stop cells had a significantly lower copy number than wild-type HPV16 cells (Fig. 12b). A low copy number accompanied by a high rate of integration is consistent with the known response of HPV to type I IFN (103–108).
(a) Representative Southern blot indicating the supercoiled (episomal) and nicked circular genome bands. (b) Quantification of episomal viral DNA in wild-type HPV16-containing cells or E5 Stop cells, with HPV16 set to 1. (c) Levels of IFNK mRNA were measured in HFKs and both integrated (Int) and episomal (Epi) HPV16+ and E5 Stop cells (Fig. 1). (d) Model: IFN-κ drives a host response that leads to integration of the viral episome. IFN-κ is positively regulated by the TGF-β pathway and negatively regulated by the EGFR pathway. E5 helps promote long-term episomal maintenance by inhibiting TGF-β signaling and enhancing EGFR signaling, thus reducing levels of IFN-κ. NS, not significant; **, P < 0.01; ***, P < 0.001. The error bars indicate standard error of the mean.
HPV genome status in HPV16+ and E5 Stop cells
To determine how genome integration correlates with IFNK expression, IFNK transcript levels were measured in E5 Stop cell RNA harvested from both episomal and integrated samples. Interestingly, while both integrated and episomal E5 Stop samples showed higher IFNK transcripts than wild-type HPV16+ cells (Fig. 12c), integrated E5 Stop samples showed significantly higher IFNK than E5 Stop cells with episomal genomes (Fig. 12c). This finding is consistent with a model in which the increased IFNK transcripts observed in E5 Stop cells could be driving the increased integration rates of the cells (Fig. 12d). Further, the significant upregulation of IFNK in episomal E5 Stop cells suggests that the loss of E5-induced suppression of IFNK occurs prior to, and may be driving, genome integration. To address the question of whether integration itself may drive IFNK upregulation, we analyzed integrated wild-type HPV16-containing cells and observed no significant change in IFNK transcript levels versus episomal HPV16+ samples (Fig. 12c). This finding confirms that integration per se is not associated with the upregulation of IFNK but that the lack of E5 causes upregulated IFNK and increased IFNK-stimulated signaling.
DISCUSSION
Our findings support a model (Fig. 12d) in which the E5 oncoprotein suppresses TGF-β signaling through inhibition of the TGFBR2 complex and stimulates the MAPK pathway through enhancement of EGFR signaling. Together, these activities prevent induction of IFNK transcription. In the absence of E5, EGFR signaling is reduced and TGFBR2 signaling is increased, which results in increased production of IFN-κ. IFN-κ activates the JAK/STAT signaling pathway, which results in higher levels of ISG expression. Chronic ISG expression results in HPV genome instability and integration (104–108). Together, these findings demonstrate a previously unrecognized role for E5 in the viral life cycle—as an innate immune evasion factor. This role of E5 may explain why the acute phenotypes of E5 mutants in culture systems have generally been mild, despite E5 being conserved among high-risk HPVs (21, 22): the virus can immortalize, maintain itself episomally, and replicate to a reasonable degree in the absence of E5, but it cannot evade the accumulating effects of the innate immune response over time.
This study adds to a growing body of literature showing a multitude of mechanisms encoded by HPVs to evade the IFN response (11, 20, 109). Which particular ISGs are most important for HPV to evade is a question whose answer remains only partially understood. In particular, which are responsible for driving HPV integration and which particular downstream mechanisms drive growth factor-dependent IFNK regulation remain to be determined, although ISGs that are capable of inhibiting HPV have been described (72, 73, 110). The transcription factors involved in controlling IFNK promoter activity have yet to be identified. Our findings suggest that the IFNK promoter utilizes transcription factor complexes other than canonical IRFs to control its expression, in contrast to other type I IFNs. We have found potential IRF1 and AP-1 binding sites in the IFNK promoter through computational studies, but they have not yet been confirmed (data not shown). Further work will be needed to investigate what transcription factors are involved in IFNK expression and activity.
Since we observed increased IFNK transcripts and protein in E5 Stop cells (Fig. 4), we also examined the levels of IFN-α/β, which were found to be equivalent to those in wild-type HPV16+ cells (Fig. 2a). We previously found that increased IFN-κ resulted in increased IFN-α/β transcription (65). E5 Stop cells have levels of E6 and E7 transcripts and protein equivalent to those of wild-type HPV-containing cells (28). Therefore, the activities of E6 and E7, which have been previously reported to inhibit the activities of IRFs and IFN-α/β production (43, 111–115), are sufficient to prevent increased IFN-α/β production in E5 Stop cells. Regardless, we were able to detect increased ISG transcripts and STAT1 phosphorylation in E5 Stop cells (Fig. 4), indicating that IFN-κ upregulation can drive a downstream type I IFN response in HPV-containing keratinocytes, even in the absence of IFN-α/β and in the presence of the other viral oncogenes. Many attempts to knock down IFN-κ have failed (data not shown), so we have not formally proven that high ISG levels are due to IFN-κ, but since IFN-κ is the only type I IFN that is upregulated in E5 Stop cells, we feel that it is reasonable to conclude that IFN-κ is responsible.
Most of what is known about the effects of E5 on cellular pathways has been determined using overexpression systems. We have now confirmed that many of those pathways are in fact regulated by E5 in keratinocytes containing the complete HPV16 genome. We used phosphoproteomics to create a more detailed picture of the cellular effects of E5 than was previously available. Our data reveal that growth factor signaling plays an important role in the regulation of IFNK by HPV16. Some of the findings seem paradoxical. First, the treatment of HPV16+ cells with EGF promoted the upregulation of IFNK transcripts (Fig. 1), but treatment with an EGFR inhibitor also increased IFNK transcripts (Fig. 8). This paradoxical finding may be due to the complexity of the EGFR family and its ligands. EGF-induced EGFR activation drives the EGFR to form heterodimers with the HER2 (ErbB2) receptor (116–118), another member of the EGFR/ErbB family of receptors; however, treatment with AG1478 stimulates the formation of EGFR-EGFR homodimers that are functionally inactive (119). Furthermore, other agonists for EGFR in addition to EGF may result in functionally different consequences (116). Regulation of other ErbB family members by E5 has been reported but is not fully understood (23, 120). Our results also show that E5 Stop cells do not have the high levels of EGFR protein seen in wild-type HPV16+ cells (Fig. 8b), consistent with our previous report but in contrast to cells containing HPV18 (28, 82). In many (but not all) E5 Stop cell lines, EGFR levels are even below that seen in HFKs (Fig. 8b), which may explain why inhibition of EGFR has no effect on IFNK levels in E5 Stop cells, as it does in HFKs (Fig. 8a). These findings suggest that the effects of E5 may differ between HPV types and that E5 may have both qualitative and quantitative effects on EGFR signaling. Much remains to be done to understand how E5 impacts the full EGFR family.
Induction of IFNK by TGF-β1 is associated with promoter demethylation (65). Since E5 Stop cells have increased levels of TGF-β signaling (Fig. 10), we expected that high IFNK transcript levels in E5 Stop cells would be due to low levels of promoter methylation. In contrast to our expectations, the IFNK promoter was methylated in E5 Stop cells (Fig. 11). This finding is reminiscent of our previous finding that TGF-β2 (in contrast to TGF-β1) could increase IFNK levels and IFN signaling without inducing IFNK promoter demethylation (65). We believe there are two possible explanations for these data. Loss of E5 expression could enable recruitment of an unknown transcription factor complex to the promoter that can stimulate IFNK expression despite the presence of promoter methylation. One possibility is the SMAD2/3 complex, which is activated by TGFBR2. Alternatively, the demethylation pathway may actually be active at the IFNK promoter but arrested at the hydroxymethylation step. Our data are consistent with at least some increased level of hydroxymethylation at the IFNK promoter in E5 Stop cells (Fig. 11c). Hydroxymethylation is increasingly recognized as a stable modification of DNA and not just as a transient intermediate of demethylation (121). Whether either or both of these mechanisms is responsible for increased IFNK transcription in the absence of E5 remains to be determined.
For over 30 years, researchers have observed that the portion of the HPV genome containing the E5 gene is frequently lost during the integration of the HPV genome in cervical cancer (122–124). Previous studies investigating the integration of high-risk HPV types in vitro have found that episomal loss of HPV16 and the emergence of clones with integrated HPV16 DNA are associated with transient activation of antiviral response genes inducible by type I IFN (104, 107). Whether this transient response is a cause or a consequence of integration is not clear, but long-term IFN treatment is also associated with accumulation of integrated HPV in culture (103, 108). Integration may be due to the ability of IFN to reduce episome numbers, which allows integrant-bearing cells to dominate the population (103), although another study showed that IFN could induce de novo integration (108). Our data indicate that the loss of E5 results in a TGF-β-mediated increase in IFN signaling. This idea is supported by the fact that the emergence of integrated clones from a previously episomal population correlates not only with an increase in IFN signaling, but also with increased TGF-β signaling (104). Although the authors attributed the consequences of increased TGF-β signaling to cell cycle progression (104), our work supports a model in which TGF-β signaling increases IFN signaling through IFN-κ (65), thereby promoting integration.
Studies investigating integration have generally considered the loss of E2 to be the main driver of cancer progression following integration (125–129). We suggest that, in addition, the loss of E5 expression may have functional consequences as either a cause or a consequence of viral genome integration. One possible scenario notes that the region of the viral genome containing E5 is 180° from the viral replication origin, so it would be the point in theta replication where the replication forks would meet and may be more prone to damage or deletion than other regions of the genome. Loss of E5 would result in failure to suppress TGF-β-mediated IFN signaling and promote integration in a single cell. Since IFNs are secreted factors, an IFN response in a single cell would promote an IFN response in the surrounding population of infected cells, resulting in further episomal loss, driving the selection of integrant-bearing cells until only integrant-bearing cells remained. This idea is supported by in vitro data. Pett et al. reported that increased IFN signaling is observed only when there is a mixed population of episomal and integrated cells; in purely episomal or integrated populations, the IFN response is suppressed (104). The IFN-mediated selection of HPV+ integrants is also likely a very rapid process, since Herdman et al. reported that IFN-β treatment could induce the presence of a previously undetectable integrant in a single passage (103). Therefore, IFN-mediated selection of integrants is likely a rapid process that would be undetectable in the clinic. Our findings suggest that the loss of E5 during integration may have previously uncharacterized functional consequences.
MATERIALS AND METHODS
Cell culture and cell lines used.HFKs were isolated from discarded and deidentified neonatal foreskins; HFKs containing HPV16 genomes (strain W12) were created by transfection and selection, as previously described (68). HPV16 genomes containing a translational stop codon at amino acid 11 in the E5 open reading frame (E5 Stop) were created as previously described (28). The presence of the mutation in E5 Stop cells was confirmed by PCR followed by sequencing (28). HFKs and keratinocyte-derived cell lines were cultivated in E medium with 5% fetal bovine serum (FBS) in the presence of mitomycin C-treated NIH-3T3 J2 fibroblast feeders (68, 130). Cell lines derived from at least three donors were used in separate experiments, and data for all the figures were compiled from at least three individual experiments. Low-passage (less than 10) HFKs and low-passage HFKs containing HPV16 or E5 Stop genomes (less than 12 passages postselection) were used in each experiment. The episomal or integrated status of the viral DNA was confirmed by Southern blotting on total cellular DNA, as previously described (65, 130). Southern blots in which the supercoiled band was visible were scored as episomal; the others were scored as integrated. The intensity of the supercoiled band was measured using ImageJ.
To generate retrovirus stocks, PT67 cells, grown in Dulbecco’s modified Eagle medium (DMEM) (Gibco) supplemented with 5% bovine growth serum (BGS) (HyClone), were seeded at a concentration of 3 million per 10-cm dish. The following day, the cells were transfected overnight with a mixture containing DMEM plus BGS, polyethylenimine (PEI), Opti-MEM medium (Gibco), 1 μg each of the pGag/Pol and pEnv packaging plasmids, and 3 μg of the retroviral plasmid to be packaged. The following day, the medium was aspirated, and fresh DMEM plus BGS was added. After a further 48 h, the supernatant was collected, and centrifuged at 1,000 rpm for 5 min. The debris-free supernatant was then filtered through a 33-mm sterile 0.22-μm syringe filter (Fisher; 09-720-004) and then flash frozen in a dry ice and ethanol bath before being stored at −80°C until it was used for infections.
To create the LXSN-wtE5 retrovirus vector, wild-type HPV16 E5 was cloned from the pUC-HPV16 (strain W12) plasmid using primers (wt16E5 5′ and wt16E5 3′ [Table 7]) that added 5′ EcoRI and 3′ BamHI restriction enzyme sites via PCR. The PCR product was then purified using a Qiagen gel extraction kit and ligated into an EcoRI- and BamHI-double-digested empty LXSN plasmid. After ligation and plasmid propagation, the sequence was verified by DNA sequencing.
Oligonucleotides used in this study
hTert-immortalized NOKs were gifts from the laboratory of Rona Scott, who received them from Karl Munger (70). The NOKs were grown in keratinocyte serum-free medium (ThermoFisher) with EGF and bovine pituitary extract added. The NOK lines were seeded overnight at a density of 3 million cells per 10-cm dish and then infected with one of several LXSN retroviral stocks containing either an empty control or the HPV16-E6, HPV16-E6/7, or HPV16-wtE5 ORF. Viral stocks were mixed with Polybrene (8 μg/ml) and medium and then added to the appropriate plate for overnight infection. Following infection, the cultures underwent selection with 250 μg/ml gentamycin (G418; Corning) for 2 weeks, followed by the outgrowth of infected cells. The cell lines were then used for various experiments, and the levels of the viral transcripts were measured by RT-qPCR to confirm expression of the desired ORF (data not shown).
A pLXSN-based retrovirus vector expressing E6/E7 and E5 was created by synthesis of a Gene Block (IDT). The sequence included the HPV16 E6/E7 sequence and E5 sequence (both strain W12) separated by the encephalomyocarditis virus (ECMV) IRES and a Kozak sequence and flanked by EcoRI and BamHI sites. DNA was rehydrated, digested with EcoRI and BamHI, and ligated into the pLXSN vector. Clones were verified by sequencing. Retrovirus particles were generated and used to infect HFK using the protocol described above. Following selection, the immortalized clones were expanded, pooled, and used for analysis. The levels of the viral transcripts were determined by RT-qPCR to confirm expression of the desired ORF (data not shown). The cell lines were maintained in E medium with 5% FBS with feeders.
The drugs used in this study were as follows: poly(I·C) (Sigma; P9582), recombinant IFN-β (PBL Assay Science; 11415-1), recombinant TGF-β1 (R&D Systems; 240-B-002), recombinant HGF (Gibco; PHG0324), recombinant EGF (Sigma; E9644), AG1478 (Sigma; 658552-5MG), SU11274 (Sigma; S9820-5MG), ruxolitinib (Selleckchem; S1378), and SB431542 (Fisher; 16-414-0). The drugs were reconstituted according to the manufacturers’ instructions and added to monolayer medium at the following concentrations: 25 U/ml (IFN-β), 5 ng/ml (TGF-β1), 100 ng/ml (EGF), 40 ng/ml (HGF), 1 μM (AG1478), 10 μM (ruxolitinib), 1 μM (SU11274), and either 5 μM or 10 μM (SB431542).
RNA extraction, RT-qPCR, and immunoblotting.Total RNA was isolated using RNA-STAT 60 (TelTest, Inc.), digested with RNase-free DNase (Promega), phenol-chloroform extracted, and reverse transcribed using qScript (Quanta) as described previously (65). Quantitative PCR was performed using PerfeCTa SYBR Green SuperMix ROX (Quanta) on an Applied Biosystems StepOne Plus real-time PCR machine using the primers shown in Table 7. Protein lysates were prepared by adding 1× Cell Signaling lysis buffer, supplemented with 1 mM phenylmethylsulfonyl fluoride (PMSF), to cells and incubating them on ice for 5 min, followed by scraping, brief sonication, and clarification by centrifugation. SDS-PAGE and Western blotting were performed essentially as described previously (131), except SDS was excluded from the TGFBR2 gels. Blocking and antibody dilution were performed using Li-Cor Odyssey blocking buffer containing 0.1% Tween 20, and images were acquired using a Li-Cor Odyssey near-infrared imaging system. The antibodies used included IFN-κ (1:500; Abnova; H00056832-M01), phospho-STAT1 (Y701) (1:1,000; Cell Signaling; 7649 P), total STAT1 (1:1,000; Cell Signaling; 9172S), TGFBR2 (1:750; Fisher; PA5-35076), phospho-Smad3 (Ser423/425) (1:2,000; Abcam; ab52903), total JAK1 (1:500; Cell Signaling; 3332), total EGFR (1:2,500; Abcam; ab32562), total Smad2/3 (1:1,000; Cell Signaling; 5678), GAPDH (glyceraldehyde-3-phosphate dehydrogenase) (1:2,500; Santa Cruz; sc-47724), tubulin (1:1,000; NeoMarkers; MS581P), OAS2 (1:1,000; Santa Cruz; sc-271117), IFIT1 (1:1,000; Proteintech; 23247-1-AP), and HPV16-E7 (1:250; Cervimax; vs 13005). Images included in the figures are unsaturated.
siRNA transfection.The DharmaFeCT siRNA protocol was followed. Briefly, half a million keratinocytes with or without an HPV genome were plated and incubated with E medium plus 5% fetal bovine serum (FBS) for 12 h on 6-well plates. siRNA was diluted to 35 nM in serum-free medium (Opti-MEM; Fisher). DharmaFeCT transfection reagent was separately diluted in Opti-MEM. The siRNA and transfection reagents were incubated at room temperature for 5 min, followed by combining the solutions and incubating them at room temperature for 20 min; 1,600 μl of E medium plus 5% FBS and 200 μl of the siRNA mixture were simultaneously added to the cells. Ninety-six hours after transfection, the cells were harvested for RNA or Western analysis as described above.
Bisulfate analysis and MeDIP.Total DNAs were isolated using phenol-chloroform extraction as described previously (21). Bisulfite treatment was performed with Zymo EZ DNA Direct kit methylation, followed by PCR, cloning, and sequencing as previously described (65). Bisulfite-specific primers previously established (the same PCR conditions described in reference 53) were used to amplify the IFN-κ promoter with AmpliTaq Gold 360 Master Mix, as previously performed (65).
For MeDIP analysis, we used a previously described protocol (102); 10 μg total cellular DNA was sonicated using a Bioruptor (Diagenode) to yield 200- to 600-bp fragments, and fragmentation was confirmed using agarose gel electrophoresis. The sonicated DNA (500 ng) was diluted in Tris-EDTA (TE) buffer and denatured at 95°C for 10 min, followed by chilling on ice. Immunoprecipitation (IP) buffer was added to a final concentration (1×) of 10 mM sodium phosphate, pH 7.0, 0.14 M NaCl, 0.05% Triton X-100, followed by control nonspecific polyclonal rabbit serum or 0.25 μg/μl 5-hmC-specific antibody (Active Motif; 39769). Samples were incubated for 2 h at 4°C. DNA-antibody complexes were captured using magnetic protein G beads and incubation with rotation for 2 h at 4°C, followed by three washes in 1× IP buffer for 5 min each time at 4°C. DNA was eluted in elution buffer (50 mM Tris-HCl, pH 8.0, 10 mM EDTA, 0.5% SDS, 0.25 mg/ml proteinase K) for 3 h at 50°C and then purified using an UltraClean PCR cleanup kit (MoBio). qPCR was used to quantify the yield of immunoprecipitated DNA using the primers listed in Table 7.
Proteomics array.The Phospho-Explorer array (Full Moon Biosystems; catalog no. PEX100) was used according to the manufacturer’s instructions. Uninfected HFKs, episomal HPV16-containing cells, and episomal E5 Stop cells were seeded at a concentration of 3 million cells in a 10-cm dish overnight. The following day, two plates for each cell line were washed 3 times in cold phosphate-buffered saline (PBS), followed by harvesting protein using the Full Moon Biosystems kit extraction buffer with PMSF added to inhibit proteases. The protein samples were biotinylated and incubated with the array slides. A GenePix Pro 6.0 array scanner was used to scan the array chips. Two independent experiments were performed, with different donor backgrounds in each. Spot intensities were normalized to the negative controls and then to the average intensity of the GAPDH spots for each array. Normalized intensities were used to calculate ratios between the various treatment groups. P values were calculated based on the log2 value of the intensity for each spot. A >2-fold change in the ratio with a significant P value (<0.05) was used as a cutoff. For each group comparison, we used DAVID to perform an in-depth statistical analysis of protein lists to determine their gene ontology (75, 76).
Statistics.The significance of wild-type versus E5 Stop integration rates was determined using Fisher’s exact test. Significance in other experiments was calculated using Welch’s unequal variances t test. RT-qPCR analysis included at least 3 technical replicates of ≥3 biological experiments, and each immunoblot densitometry analysis included at least 1 technical replicate of ≥6 biological experiments. Western quantification included ≥4 biological replicates.
ACKNOWLEDGMENTS
We thank the members of the Center for Molecular and Tumor Virology for helpful discussions. We thank Rona Scott for the gift of NOKs, Arjun Yusufji for technical assistance, and Liudmila Chesnokova for assistance in processing the phosphoproteomics array.
This work was supported by grants from the National Institute of Allergy and Infectious Diseases (R01AI118904); the National Institute of General Medical Sciences (P30GM110703); the National Heart, Blood, and Lung Institute (HL141155); and the Feist-Weiller Cancer Center.
The content is solely our responsibility and does not necessarily represent the official views of the National Institutes of Health or the Feist-Weiller Cancer Center.
FOOTNOTES
- Received 18 September 2019.
- Accepted 23 October 2019.
- Accepted manuscript posted online 30 October 2019.
Supplemental material is available online only.
- Copyright © 2020 American Society for Microbiology.
REFERENCES
- 1.↵
- 2.↵
- 3.↵
- 4.↵
- 5.↵
- 6.↵
- 7.↵
- 8.↵
- 9.↵
- 10.↵
- 11.↵
- 12.↵
- 13.↵
- 14.↵
- 15.↵
- 16.↵
- 17.↵
- 18.↵
- 19.↵
- 20.↵
- 21.↵
- 22.↵
- 23.↵
- 24.↵
- 25.↵
- 26.↵
- 27.↵
- 28.↵
- 29.↵
- 30.↵
- 31.↵
- 32.↵
- 33.↵
- 34.↵
- 35.↵
- 36.↵
- 37.↵
- 38.↵
- 39.↵
- 40.↵
- 41.↵
- 42.↵
- 43.↵
- 44.↵
- 45.↵
- 46.↵
- 47.↵
- 48.↵
- 49.↵
- 50.↵
- 51.↵
- 52.↵
- 53.↵
- 54.↵
- 55.↵
- 56.↵
- 57.↵
- 58.↵
- 59.↵
- 60.↵
- 61.↵
- 62.↵
- 63.↵
- 64.↵
- 65.↵
- 66.↵
- 67.↵
- 68.↵
- 69.↵
- 70.↵
- 71.↵
- 72.↵
- 73.↵
- 74.↵
- 75.↵
- 76.↵
- 77.↵
- 78.↵
- 79.↵
- 80.↵
- 81.↵
- 82.↵
- 83.↵
- 84.↵
- 85.↵
- 86.↵
- 87.↵
- 88.↵
- 89.↵
- 90.↵
- 91.↵
- 92.↵
- 93.↵
- 94.↵
- 95.↵
- 96.↵
- 97.↵
- 98.↵
- 99.↵
- 100.↵
- 101.↵
- 102.↵
- 103.↵
- 104.↵
- 105.↵
- 106.↵
- 107.↵
- 108.↵
- 109.↵
- 110.↵
- 111.↵
- 112.↵
- 113.↵
- 114.↵
- 115.↵
- 116.↵
- 117.↵
- 118.↵
- 119.↵
- 120.↵
- 121.↵
- 122.↵
- 123.↵
- 124.↵
- 125.↵
- 126.↵
- 127.↵
- 128.↵
- 129.↵
- 130.↵
- 131.↵