Human Papillomavirus 16 (HPV16) E2 Repression of TWIST1 Transcription Is a Potential Mediator of HPV16 Cancer Outcomes

HPV16-positive cancers have a better clinical outcome that their non-HPV anatomical counterparts. Furthermore, the presence of HPV16 E2 RNA predicts a better outcome for HPV16-positive tumors; the reasons for this are not known.

HPV infects basal epithelial cells, and upon nuclear entry, a host of cellular factors activate transcription from the viral long control region (LCR) (6). The resultant viral transcript is processed into individual gene RNAs that are then translated. The viral oncoproteins E6 and E7 target several cellular proteins and disrupt their functions, including the tumor suppressors p53 and pRb, respectively (7,8). Both p53 and pRb are transcription factors; therefore, the presence of E6 and E7 in cells results in a disruption of host gene transcription that contributes to the oncogenic properties of HPV. HPV uses two proteins to regulate replication of the viral genome. The E2 protein forms homodimers via a carboxyl terminus domain and binds to four 12-bp palindromic sequences within the viral LCR. Three of these surround the viral origin of replication (9), and via the amino terminus of the E2 protein, the viral helicase E1 is recruited to the A/T-rich viral origin of replication (10). E1 forms a dihexameric helicase that then replicates the viral genome in association with host polymerases (11)(12)(13)(14). Following infection, the virus replicates to around 50 copies per cell to establish the infection. This copy number is maintained as the infected cell migrates through the epithelium before amplifying in the upper layers; the L1/L2 structural proteins are then expressed and the viral genomes encapsulated, resulting in viral particles that egress from the upper layers of the epithelium (15).
As well as acting as a replication factor, the E2 protein can regulate transcription. Where E2 sites are present upstream from a heterologous promoter such as herpes simplex virus 1 thymidine kinase (tk) promoter, E2 can activate transcription (16,17). In addition, overexpression of E2 can repress transcription from the viral LCR (18)(19)(20). Given the ability of E2 to act as a transcription factor, the ability of E2 to regulate host gene transcription has been studied. E2 can regulate transcription via AP1 (21)(22)(23)(24)(25), nuclear receptors (26), and C/EBP (27). Transient E2 overexpression identified global gene changes induced by E2 (28)(29)(30). To gain a greater understanding of E2 regulation of host gene transcription by physiologically tolerated levels of E2, we generated stable cell lines expressing E2 and identified E2-induced host gene expression changes. This was originally done in U2OS cells (31,32). However, we wished to develop a more physiologically relevant model, and to do this, we used N/Tert-1 cells, foreskin keratinocytes immortalized by telomerase. We generated N/Tert-1 cell lines stably expressing HPV16 E2. For comparison, we also prepared N/Tert-1 cells that contained the entire HPV16 genome and have previously demonstrated that N/Tert-1 cells support late stages of the HPV16 life cycle, making it an appropriate model for the study of HPV16 (33). Transcriptome sequence (RNA-seq) analysis demonstrated a significant overlap between genes regulated by E2 and those regulated by the entire HPV16 genome (34). Many innate immune genes were repressed by HPV16 E2 and the entire genome, and these genes are also repressed by E6 and E7 expression. We recently demonstrated that one of these genes, SAMHD1, is a restriction factor for HPV16, as it controls the viral life cycle in the differentiating epithelium (35). However, we wished to identify a gene that was only regulated by E2, not by E6/E7. Our RNA-seq analysis predicted that TWIST1 was transcriptionally repressed by E2 and the entire HPV16 genome in N/Tert-1 cells; in addition, TWIST1 was also downregulated in HPV16-positive versus HPV16negative head and neck cancer (34).
TWIST1 is a basic helix-loop-helix transcription factor critical for promoting epithelial-to-mesenchymal transition (EMT) and embryogenesis (36,37). EMT is a critical process in cancer, as this trait is associated with high-grade malignancy and resistance to chemotherapeutic agents (38)(39)(40)(41). EMT is an epigenetic process that proceeds independently from DNA mutations (42). In addition, EMT promotes immune escape of cancer cells (43).
The potential repression of TWIST1 expression by HPV16 E2 is intriguing, as it could play a role in dictating therapeutic outcomes. Several reports have demonstrated that the expression of E2 in HPV16-positive tumors predicts improved survival, and repression of TWIST1 would correlate with this improved survival (44)(45)(46). Here, we demonstrate that TWIST1 is transcriptionally repressed by E2 in N/Tert-1 cells and that TWIST1 is downregulated in HPV-positive head and neck cancer versus HPV negative. The mechanism of E2 repression is not due to DNA methylation but involves direct binding of E2 to the TWIST1 promoter, repressing transcription, and inducing repressive histone markers. We demonstrate that TWIST1 target genes are also downregulated in E2expressing cells and that wound healing is compromised. As EMT is a process that occurs during wound healing, the observed change in wound healing suggests that E2 contributes to EMT suppression in N/Tert-1 cells (47). Two HPV16-positive head and neck cancer cell lines were also studied, one that has episomal viral genomes and therefore expresses E2, and the other with an integrated HPV16 genome that has lost E2 expression. In the episomal genome-containing cell line, the presence of E2 resulted in decreased levels of TWIST1 compared with the integrated non-E2-expressing cell line. Finally, neither E6 nor E7 is able to regulate the expression of TWIST1. Overall, our results support the idea that E2 represses TWIST1 expression during the HPV16 life cycle and that this downregulation persists into HPV16-positive tumors. TWIST1 repression would promote a better patient outcome, a hallmark of E2 expression (44)(45)(46).

RESULTS
TWIST1 is transcriptionally repressed by E2, but not E6/E7. Our RNA-seq analysis of N/Tert-11E2 (expressing HPV16 E2 only) and N/Tert-11HPV16 (containing the entire HPV16 genome) demonstrated that there was a highly statistically significant overlap between the genes regulated by E2 and the entire HPV16 genome (34). In addition, using our analysis of The Cancer Genome Atlas (TCGA) data, we observed a significant downregulation of TWIST1 in HPV16-positive versus negative head and neck cancers (HNSCCs) (33,34). Figure 1A summarizes the expression of TWIST1 in both episomal and integrated HPV16-positive versus HPV16-negative HNSCCs. Because integration disrupts E2 gene expression, integration and episomal groups were determined depending on their E2 expression (46,48). TWIST1 mRNA expression data were obtained from 528 TCGA tumors using the cBio Cancer Genomic Portal (49,50). We have previously characterized HPV16 status and viral integration in these samples (46). TWIST1 mRNA expression was then compiled and reported by HPV16 status (Fig. 1A). We found that TWIST1 mRNA is expressed at statistically significantly lower levels in episomal HPV16-positive tumors than those that were either integrated or HPV16 negative. There was no statistical difference in TWIST1 expression between integrated and HPV-negative tumors. Figure 1B summarizes our previous data from our RNA-seq analysis in N/Tert-1 cells. TWIST1 was found to be significantly downregulated in N/Tert-1 cells expressing HPV16 E2 or the entire viral episome compared to parental control cells (33,34). Figure 2A validates the downregulation of TWIST1 RNA expression in N/Tert-1 cells expressing E2 or the entire HPV16 genome (compare lanes 2 and 3, respectively, with lane 1). We have also previously reported the overexpression of HPV16 E6 and E7 in the same cell background, and these oncogenes did not alter the expression of TWIST1 (lane 4). This is in contrast with our prior study of innate immune response gene regulation by HPV16, which demonstrated that E2, E6, and E7 can all repress these genes (34). Therefore, TWIST1 is the first gene we have determined to be likely exclusively regulated by E2 during the HPV16 life cycle. To confirm that the RNA expression was reflected at the functional protein level, Western blotting was carried out for TWIST1 (Fig. 2B). There is a clear downregulation of TWIST1 protein expression in cells expressing E2 or HPV16 (compare lanes 2 and 3, respectively, with lane 1) but not in cells expressing E6/E7 (compare lane 4 with lane 1). This was repeated another two times and the results quantified (Fig. 2C); both E2 and HPV16 induce a significant reduction in TWIST1 protein levels, while E6/E7 do not. To confirm that TWIST1 function is also downregulated in the presence of E2, we monitored expression of TWIST1 target genes vimentin and N-cadherin in N/Tert-11Vec, N/Tert-11E2, N/Tert-11HPV16, and N/Tert-11E6/E7 cells (Fig. 2D). The expression levels of these TWIST1 target genes are reflective of the TWIST1 levels in the cell. This demonstrates that there is a functional loss of TWIST1 in E2 and HPV16 expressing N/Tert-1.
E2 binds to the TWIST1 promoter region and directly represses transcription. Our previous work demonstrated that E2 represses the transcription of innate immune genes by regulating the DNA methylation of the corresponding promoters (34). To determine whether E2 is also regulating TWIST1 expression via DNA methylation of the promoter, we treated the cells with 1 mM decitabine, a DNA methylase inhibitor that relieves E2-mediated repression of innate immune genes (Fig. 3). Figure 3A demonstrates that the drug has worked, as there is an increase in both MX1 and IFIT1 when decitabine is added, as we observed previously (34). The increase occurs in the vector control cells (compare lanes 3 and 4 versus 1 and 2), but the relief of repression is greater in the presence of E2 (compare lanes 7 and 8 versus 5 and 6) and HPV16 (compare lanes 9 and 10 versus 11 and 12); these results duplicate those we observed previously (34). However, the TWIST1 levels were not altered in any of the cell lines following treatment with decitabine. Therefore, TWIST1 is regulated by E2 differently from innate immune response genes and is independent of DNA methylation.
We next investigated whether E2 could directly repress transcription from the TWIST1 promoter. A construct containing the TWIST1 promoter upstream from the luciferase gene (pTWIST-luc) was cotransfected with E2 into N/Tert-1 cells to determine whether E2 can repress expression directly from the TWIST1 promoter (Fig. 4A). The expression of E2 resulted in a ;10-fold reduction in luciferase activity, demonstrating that E2 directly represses this promoter. To confirm that this was not due to E2 toxicity in the transiently transfected cells, we included transcriptional activation studies on  (46,48). Four hundred sixty-six samples were HPV negative. Forty samples retained HPV16 episome status, and 17 samples had integrated HPV16 genomes. Firehose Legacy TWIST1 mRNA expression data were obtained using cBioPortal. TWIST1 mRNA between HPV16 status was then compared using two-tailed Student's t test. Vertical axis in log 10 . (B) N/Tert-1 RNA-sequencing results adapted from reference 34. TWIST1 expression was compared to parental N/Tert-1 cells with no HPV16 or HPV16 E2. **, P , 0.05 using Bonferroni correction when applicable.
ptk6E2-luc, a plasmid we have shown is responsive to E2. Figure 4B demonstrates that cotransfection with E2 increases transcription from ptk6E2-luc, demonstrating that there is no overall toxicity induced by E2 in the N/Tert-1 cells. To determine whether E2 can bind directly to the TWIST1 promoter, we carried out a chromatin immunoprecipitation assay (ChIP) using N/Tert-11Vec, N/Tert-11E2, and N/Tert-11HPV16 cells using a sheep E2 antibody as previously described (51,52) (Fig. 4C). There was a significant increase in signal in N/Tert-11E2 and N/Tert-11HPV16 compared with the signal obtained with N/Tert-11Vec (compare lanes 2 and 3, respectively, with lane 1), demonstrating that E2 binds to the TWIST1 promoter region when overexpressed but also in the context of the entire HPV16 genome. We next investigated whether the levels of repressive chromatin markers implicated in TWIST1 regulation are changed in the presence of E2. H3K9me2 is a repressive chromatin marker involved in the regulation of TWIST1 expression (53,54), and the levels of this marker on the TWIST1 promoter are increased in N/Tert-11E2 and N/Tert-11HPV16 compared with N/Tert-11Vec ( Fig. 4D; compare lanes 2 and 3, respectively, with lane 1). These results suggest that E2 interacts with the TWIST1 promoter directly and modifies the local epigenetic environment around the promoter, leading to a reduced level of transcription.
E2 expression reduces the rate of wound healing in N/Tert-1 cells. The repression of TWIST1 expression by E2 suggested suppression of the EMT phenotype in  4). GAPDH is shown as an internal control. Western blotting was visualized using a Li-Cor system. (C) Western blotting was quantitated, and TWIST1 protein expression was calculated relative to vector using ImageJ. (D) Repression of TWIST1 by E2 and HPV16 leads to reduction in EMT marker expression. N/Tert-1 cells were harvested for RNA and processed for cDNA. qRT-PCR was performed for TWIST1 target genes CDH2 and VIM, which encode the proteins N-cadherin and vimentin, respectively. Results are expressed as fold change from that observed in the vector control N/Tert-1 cells (lane 1). Data in panels A, C, and D represent the averages of at least 3 independent experiments, and error bars indicate standard error of the mean. *, P , 0.05. N/Tert-1 cells. EMT is not a defined status but is a spectrum of phenotypes that range from epithelial through to fully mesenchymal (47). The N/Tert-11E2 cells do not look appreciably different from N/Tert-1 control cells; therefore, we propose that there is a subtle influence of E2 on the EMT status of these cells. During wound healing, there is an EMT transition of the wounded epithelia cells that promotes migration and eventual wound closure (55). Our previous studies on U2OS cells demonstrated that the expression of E2 slowed the ability of these cells to close wounds in monolayer cells (32); therefore, we carried out "scratch" assays with our N/Tert-1 cells in order to investigate wound healing. Figure 5A shows images of the results from the woundhealing experiment. Twenty hours following wound induction, the wound is almost completely healed in N/Tert-11Vec cells (top 3 panels). However, with N/Tert-11E2 and N/Tert-11HPV16, there is a failure to close the wound after 20 h (second and third from top panels, respectively). To determine whether the oncogenes E6/E7 play a role in regulating wound healing, we also determined wound closure in N/Tert-11E6/E7 cells (characterized and described in reference 34). The presence of the viral oncogenes made no difference to the wound healing (compare the bottom three panels with the top three). Therefore, the failure to close the wound is reflective of TWIST1 levels and E2 expression. This assay was repeated several times and the data quantified (Fig. 5B). There is a significant delay in wound healing 12 and 20 h after the initial "scratch" in both the N/Tert-11E2 and N/Tert-11HPV16 cell lines compared with N/Tert-11Vec cells.
E2 and TWIST1 levels inversely correlate in HPV16-positive head and neck cancer cell lines. Figure 1 demonstrates that there is, on average, less expression of TWIST1 in HPV-positive head and neck cancers than in HPV-negative head and neck cancers. Our results suggest that the expression of E2 is influential on the expression of TWIST1 within the HPV-positive group. Previous studies suggested that UMSCC104 contained episomal HPV16 genomes (and therefore retained E2 expression), while UMSCC47 had integrated viral genomes (and therefore had lost E2 expression). We confirmed that UMSCC104 expressed E2 and that UMSCC47 did not (Fig. 6A). We next investigated the TWIST1 levels in these cell lines. TWIST1 RNA expression was significantly downregulated in UMSCC104 compared with UMSCC47 (Fig. 6B), and the levels of RNA expression were reflected in a downregulation of TWIST1 protein levels (Fig. 6C). Western blotting was repeated several times and quantitated; there is significantly less TWIST1 protein expressed in UMSCC47 than in UMSCC104. The results in these two HPV16-positive head and neck cancer cell lines support the model of E2 repressing TWIST1 gene expression.

DISCUSSION
The HPV E2 protein plays multiple crucial roles during the viral life cycle. It is essential for replication of the viral genome, can regulate the segregation of the viral genome into daughter cells, and has transcriptional properties with the potential to regulate transcription from the viral and human genomes (9, 56). Regulation of host gene expression by E2, and its importance in the viral life cycle, is relatively understudied compared with its role in viral replication and segregation. Our recent work demonstrates that E2 regulates host gene transcription that is relevant during the viral life cycle (34). That report demonstrated an overlap in the ability of E2 and E6/E7 to repress the expression of innate immune genes. In this current report, we extend our observations on the transcriptional regulatory functions of E2. We demonstrate that E2 can repress TWIST1 gene expression and that this repression is via a distinct mechanism compared to innate immune gene repression. E6 and E7 did not affect the levels of TWIST1 RNA and protein; therefore, it is likely that TWIST1 expression is regulated by E2 during the HPV16 life cycle. For innate immune genes, methylation of the gene promoters played a significant role in repressing expression, but for TWIST1, DNA methylation plays no role. Rather, E2 binds to the TWIST1 promoter and actively represses transcription from this region. E2 also increases H3K9me2 at the promoter, a hallmark of transcriptionally repressed genes previously implicated in regulation of TWIST1 transcription (53). To our knowledge, this is the first time E2 has been demonstrated to directly bind to a host gene promoter and induce repressive epigenetic markers. The repression of TWIST1 by E2 resulted in downregulation of TWIST1 target genes vimentin and N-cadherin. During wound healing, there is an EMT-like transition that promotes wound closure, and wound healing in E2-expressing N/Tert-1 is compromised compared to parental cells (55).
Repression of TWIST1 expression by E2 may play a role in HPV cancer outcomes. Overexpression of TWIST1 is associated with poorer overall survival in head and neck cancer (57), while E2 expression is associated with improved outcomes in HPV-positive head and neck cancer (44)(45)(46). We demonstrate here that TWIST1 is downregulated in HPV16-positive tumors that retain E2 expression compared with HPV16 tumors that have no E2 expression or are HPV negative. We therefore highlight this correlation between E2 expression and TWIST1 repression, which may contribute to the better clinical outcomes of HPV16-positive tumors that retain E2 expression. We demonstrate, in two HPV16-positive head and neck cancer cell lines (one that is E2 positive, one that is negative), that E2 expression inversely correlates with TWIST1 expression, which correlates with our observations in N/Tert-1 cells and TCGA data sets. While these lines are not isogenic, they provide a model for studying the repression of TWIST1 in cancer cell lines that is potentially mediated by expression of the E2 protein. Not all E2-positive head and neck tumors have downregulation of TWIST1 compared with E2-negative tumors as evidenced by TCGA data, but there is a significant trend for downregulation. Not only has elevated TWIST1 been associated with poorer survival in several cancers, including head and neck, but TWIST1 protein attenuates the response to chemotherapeutic drugs, which provides a rationale for the worse clinical outcomes in high TWIST1 expression patients (47,55).
The mechanism of E2 repression of TWIST1 expression remains to be fully elucidated. The TWIST1 promoter is methylated in cancer, although this did not correlate with low TWIST1 RNA or protein levels (58). Methylation does not play a role in the mechanism that E2 uses to repress TWIST1 levels, as treatment of cells with decitabine did not relieve TWIST1 repression; therefore, E2 uses multiple mechanisms to regulate host gene transcription. We demonstrate here for the first time that E2 binds to and represses transcription from the TWIST1 promoter and modifies chromatin around the promoter start site by inducing elevated levels of H3K9me2, a repressive marker. The transcription factor SP1 is involved in the basal transcription levels from the TWIST1 promoter, and displacement of this factor from the promoter repressed TWIST1 transcription (59). The E2 protein has been shown to displace SP1 from HPV LCRs, resulting in transcriptional repression, and may act similarly on the TWIST1 promoter (60,61). This active repression may be required to block TWIST1 expression, as STAT3 is activated by HPV in cervical cancer cells (62), and active STAT3 is an activator of TWIST1 expression, promoting EMT (63). We demonstrate here that there is an increase in the level of the repressive marker H3K9me2 at the TWIST1 promoter in the presence of E2. The chromatin repressor NuRD complex can be recruited to the TWIST1 promoter via DOC1, resulting in transcriptional repression via induction of a nucleosome on the TWIST1 promoter. It is feasible that E2 recruits a repressor complex to the TWIST1 promoter in order to regulate transcriptional repression via nucleosome assembly (64). E2 functionally interacts with BRD4, and this interaction is involved in regulating E2 repression of the viral LCR. Again, such a mechanism may be used by E2 to contribute to the repression of the TWIST1 promoter (65). Future studies will characterize the mechanism(s) that E2 uses to repress the TWIST1 promoter. Enhancing our understanding of this mechanism is important, as E2 offers a model system for repressing TWIST1 transcription. For example, if H3K9me2 methylation of the TWIST1 promoter is a major contributor to repression by E2, then induction of this modification via drug treatment could be an opportunity to promote a better response of high TWIST1-expressing tumors to chemotherapeutic agents. In addition, if TWIST1 repression is an important mechanism promoting the HPV16 life cycle in epithelial cells, then reversing this repression offers an opportunity to disrupt this process and alleviate HPV16 infections and disease. Another question is, why does HPV16 repress the expression of TWIST1? While TWIST1 expression is important in mouse embryogenesis, demonstrating an important role in cellular proliferation and differentiation, there is less known about the effect of TWIST1 on epithelial cell differentiation (37). HPV manipulates epithelial differentiation in order to create an environment that supports viral replication (66). While TWIST1 is clearly a promoter of EMT, the effect of this protein on the epithelial differentiation pathway is less clear and is worthy of future study. Perhaps the repression of TWIST1 plays a role in promoting an environment that allows HPV16 replication in the infected epithelium. TWIST1 and HPV16 E2 interact with the NF-κB pathway, and downregulation of TWIST1 may promote the ability of E2 to regulate NF-κB, which could be important for the HPV16 life cycle (67,68). TWIST1 has been shown to downregulate expression of C/EBPa, a protein that E2 functionally interacts with to regulate transcription (27,69). Therefore, downregulation of TWIST1 expression may prevent disruption of E2's ability to regulate host gene transcription via C/EBPa during the viral life cycle.
Recently, it has been demonstrated that C/EBPa is a crucial factor for epithelial maintenance and prevents EMT; thus, E2 downregulation of TWIST1 may enhance the ability of C/EBPa to carry out this function (70). This would also promote the enhanced epithelial status of the E2-positive cells.
In conclusion, this report demonstrates that HPV16 E2 downregulates the expression of TWIST1. This occurs in cells that contain the entire HPV16 genome and also occurs in HPV16-positive cancers that are E2 positive. Given the important role that TWIST1 plays in EMT, cancer outcomes, and response to chemotherapeutic agents, a fuller understanding of the interaction of HPV16 with TWIST1 is warranted.

MATERIALS AND METHODS
Differential expression in TCGA. Head and neck squamous cell carcinoma (HNSCC) TWIST1 mRNA expression data were obtained from The Cancer Genome Atlas using cBio Cancer Genomics Portal (49,50). HNSCC tumor samples (Firehose Legacy) were analyzed for HPV status and viral genome integration as previously described utilizing the HPV16 E2:E7 ratio of mRNA sequencing reads (46,71). Samples without available TWIST1 mRNA expression data were omitted. Data for 528 HNSC samples were available. Of these, 522 were identified with available TWIST1 mRNA expression and HPV status, which was correlated and reported using R. Statistical analysis was performed using two-way Student's t test with Bonferroni correction for two comparisons.
SYBR green reverse transcription-quantitative PCR. RNA was isolated using the SV Total RNA isolation system (Promega) following the manufacturer's instructions. We reverse transcribed 2 mg of RNA into cDNA using the high-capacity reverse transcription kit (Applied Biosystems). cDNA and relevant primers were added to PowerUp SYBR green master mix (Applied Biosystems), and real-time PCR was performed using 7500 Fast real-time PCR system. The primer sequences utilized were as follows: Decitabine treatment. N/Tert-1 cells were plated at a density of 1.5 Â 10 5 in 6-well plates (60 mm 2 / well). The following day, cells were treated with 1 mM decitabine or 1 mM DMSO for 72 h as previously described (34). Afterward, the cells were harvested and processed for reverse transcription-quantitative PCR (qRT-PCR) as described above.
Immunoblotting. N/Tert-1, UMSCC47, and UMSCC104 cells were trypsinized, washed with 1Â phosphate-buffered saline (PBS), and resuspended in 2Â pellet volume protein lysis buffer (0.5% Nonidet P-40, 50 mM Tris [pH 7.8], and150 mM NaCl) supplemented with protease inhibitor (Roche Molecular Biochemicals) and phosphatase inhibitor cocktail (Sigma). The cell suspension was lysed for 30 min on ice and then centrifuged for 20 min at 184,000 relative centrifugal force (rcf) at 4°C. Protein concentration was determined using the Bio-Rad protein estimation assay. We boiled 25-mg protein samples in equal volume 2Â Laemmli sample buffer (Bio-Rad). Samples were run down a Novex 4 to 12% Tris-glycine gel (Invitrogen) and transferred onto a nitrocellulose membrane (Bio-Rad) at 30 V overnight using the wet blot method. Membranes were blocked with Odyssey (PBS) blocking buffer (diluted 1:1 with 1Â PBS) at room temperature for 1 h and probed with indicated primary antibody diluted in Odyssey blocking buffer. Membranes were then washed with PBS supplemented with 0.1% Tween (PBS-Tween) and probed with the indicated Odyssey secondary antibody (goat anti-mouse IRdye 800CW or goat anti-rabbit IRdye 680CW) diluted in Odyssey blocking buffer at 1:10,000. Membranes were washed and underwent infrared scanning using the Odyssey CLx Li-Cor imaging system. Immunoblots were quantified using ImageJ. The following primary antibodies were used for immunoblotting: HPV16 E2 (TVG 261) at 1:1,000 dilution from Abcam; glyceraldehyde-3-phosphate dehydrogenase (GAPDH) (catalog no. sc-47724) and p53 (catalog no. sc-47698) at 1:1,000 dilution from Santa Cruz Biotechnology, and TWIST1 (catalog no. 25465-1-AP) at 1:300 from Proteintech.
ChIP assay. N/Tert-1 cells were plated at a density of 2 Â 10 6 in 150-mm 2 plates. The following day, the cells were harvested via scraping and processed for chromatin as previously described (72). Chromatin concentration was determined with a NanoDrop spectrophotometer. Approximately 100 mg of chromatin was used per antibody experiment. The following antibodies were used for ChIP: 2 ml sheep anti-HPV16 E2 (amino acids 1 to 201) prepared and purified by Dundee Cell Products, United Kingdom; 2 mg rabbit anti-histone H3K9me2 (Abcam; catalog no. ab1220). Chromatin was then processed for qPCR. ChIP DNA primers include Twist1 forward, 59-TCAGGCCAATGACACTGCT-39, and Twist1 reverse, 59-GACGGTGTGGATGGCCCCGA-39.
Transcription assay. We plated out 5 Â 10 5 N/Tert-1 cells on 100-mm 2 plates, and they were transfected 24 h later with either 1 mg HPV16 E2 plasmid and 1 mg pTWIST1-Luc/ptk6E2-luc, or 1 mg pTWIST1-Luc/ptk6E2-luc alone using Lipofectamine 3000 (Thermo Fisher Scientific) according to the manufacturer's instructions as previously described (73). pTWIST1-luc contains the human promoter and has been described previously (74,75). ptk6E2-luc has been described (16,17). Briefly, the cells were harvested 72 h posttransfection utilizing the Promega reporter lysis buffer and analyzed for luciferase using the Promega luciferase assay system. Concentrations were normalized to protein levels, as measured by the Bio-Rad protein estimation assay mentioned above. Relative fluorescence units were measured using the BioTek Synergy H1 hybrid reader.
Wound healing assay. We plated 2.5 Â 10 5 N/Tert-1 cells at the center of 6-well plates (60 mm 2 / well). The cells were left to grow to confluence. Afterward, the monolayer was scratched using a 1,000ml pipette tip, creating an ;1-mm wound. Wounds were imaged at 0-, 12-, and 20-h intervals (Zeiss; Axiovert 200 M microscope). Multiple images were taken randomly along the wounds, and measurements were taken from the leading cell edge. Wounds were measured using ImageJ software, and wound healing was calculated as a ratio of the 0-h time point.
Statistical analyses. The standard error was calculated from no less than three independent experiments. The only exception was in the TWIST1-Luc transcription assay where statistical power was achieved after 2 replicates. Significance was determined using two-tailed Student's t test. Bonferroni correction for significance was utilized when indicated.