The HPV E2-Host Protein-Protein Interactions: A Complex Hijacking of the Cellular Network.

Over 100 genotypes of human papillomaviruses (HPVs) have been identified as being responsible for unapparent infections or for lesions ranging from benign skin or genital warts to cancer. The pathogenesis of HPV results from complex relationships between viral and host factors, driven in particular by the interplay between the host proteome and the early viral proteins. The E2 protein regulates the transcription, the replication as well as the mitotic segregation of the viral genome through the recruitment of host cell factors to the HPV regulatory region. It is thereby a pivotal factor for the productive viral life cycle and for viral persistence, a major risk factor for cancer development. In addition, the E2 proteins have been shown to engage numerous interactions through which they play important roles in modulating the host cell. Such E2 activities are probably contributing to create cell conditions appropriate for the successive stages of the viral life cycle, and some of these activities have been demonstrated only for the oncogenic high-risk HPV. The recent mapping of E2-host protein-protein interactions with 12 genotypes representative of HPV diversity has shed some light on the large complexity of the host cell hijacking and on its diversity according to viral genotypes. This article reviews the functions of E2 as they emerge from the E2/host proteome interplay, taking into account the large-scale comparative interactomic study.

E2 PROTEIN IN THE HPV LIFE CYCLE AND PATHOGENESIS

The E2 protein is pivotal for both the productive cycle and persistent infections by human papillomaviruses (HPVs). The productive cycle causes benign hyperproliferative lesions (genital or cutaneous warts). However, some lesions have the potential to progress toward malignancies in case of unresolved infection with high-risk HPV (HR-HPV). Carcinogenic conversion is a complex process occurring over long periods of time, thus requiring long-lasting infections. Indeed, persistence of HR-HPV is considered as a major risk factor for cancer development.

The E2 proteins have been mostly studied for the genital α-type HPV, more particularly HR-HPV 16 and 18, which are associated with about 70% of HPV-induced genital cancers and therefore represent major health concerns. Less information is available about E2 of low-risk genital HPV (LR-HPV) or of cutaneous HPVs. Yet, the complete deciphering of E2 activities in the productive life cycle is still required to find a way to eradicate genital and cutaneous warts. In addition, decoding the roles of E2 in viral persistence is crucial to improve the comprehension of the HR-HPV associated carcinogenic conversion, which arises from persistently infected cells. It also may help to understand the pathology of some cutaneous HPVs which persist without generating lesions (asymptomatic infections) in the general population, and trigger productive infection only in immunocompromised individuals.

The E2 proteins contain 350 to 500 amino acids consisting of two conserved modular domains, a 200 amino acids N-terminal domain referred as TAD (TransActivation Domain), and a 100 amino acid C-terminal DNA Binding Domain (DBD). These two domains are separated by an unstructured Hinge region, not conserved and of variable length among the different HPV genotypes. The crystal structure of the two conserved domains has been determined for several genital α-HPV E2 proteins [1, 2], but no information is available concerning the structure of E2 from the cutaneous HPVs.

The DNA binding domain of E2 exhibits a typical three-dimensional structure, which has homology to other viral factors and is considered prototypic of some viral regulation factors [3]. This domain contains both a strong dimerization interface and amino acids making contacts with DNA located along an α-helix [1]. Indeed, E2 binds as a dimer to specific palindromic DNA sequences ACCGN4 CGGT designated E2 binding sites (E2BS). At least four of these E2BS are located within the regulatory region (LCR for Long Control Region) of all HPV genomes [4] (Fig. ). Through the binding to the regulatory region, the E2 protein recruits at the viral genome a range of host cell factors necessary to support its transcription, replication, and mitotic segregation. By binding to the same E2BS, E2 ensures the coordinated regulation of viral DNA transcription and replication, thereby appearing as a pivotal factor of the early/late switch of the productive viral life cycle. An effect on mRNA splicing has been detected as well, suggesting that E2 controls viral gene expression through both transcription regulation and mRNA processing.

Several groups have demonstrated a number of additional functions that are independent of E2 binding to the viral genome, such as NFκB activation [5], induction of apoptosis (for review [6]) or regulation of the host cell cycle [7]. These “autonomous” activities directly impact on the host cells and are likely to be required for the implementation of the productive viral life cycle within stratified epithelia. For example, the stimulation of TNF-induced NFκB activation by E2 was proposed to play a role in the commitment of infected keratinocytes to activation of the differentiation program, which is necessary for the productive phase of viral life cycle [5]. Some of these autonomous E2 functions are specific to HR-HPV, leading to the notion that the E2 proteins have specially developed activities in correlation with the oncogenic power of HPV. For example, a proapoptotic activity is specifically associated with the E2 proteins of HR-HPV and not with the LR-HPV [6, 8, 9]. It was hypothesized that the HR-HPV E2 proteins evolved apoptotic activities to counterbalance the transforming properties of the E6/E7 oncoproteins, which would otherwise be detrimental for progression of viral life cycle. In addition, the HR-HPV E2 proteins cause a G2/M cell cycle arrest, which can induce genomic instability, recognized as a precursor of carcinogenesis [7]. In the case of cutaneous HPV8, E2 was shown to be sufficient to induce the development of skin cancer in transgenic mice [10]. The correlation between E2 and modulation of mechanisms involved in cell transformation has raised the hypothesis that E2 may have a direct role in the early steps of carcinogenic conversion, and thus contribute to the oncogenic potential of HR-HPV in addition to E6 and E7. In late stages of carcinogenesis, the genital HPV genome is integrated into the host chromatin, with the preferential disruption of the E2 gene [11, 12]. The E6/E7 genes are subsequently upregulated and trigger key events for cell transformation. During the advanced steps of malignant progression, the pro-apoptotic activity of HR-HPV E2 proteins may favor the selective inactivation of E2-expressing cells, thereby promoting the clonal selection of E6/E7 over-expressing cells. The role of E2 in the different stages of HPV-associated carcinogenesis has been extensively discussed recently [13].

THE E2 INTERACTOME

Like other viral proteins, E2 exerts its functions by targeting different cell factors. Indeed, almost all of the activities associated with E2 have been linked to direct interaction with host cell factors, either to support viral genome expression, segregation and replication, or to provide conductive cell conditions.

To date, most of the interaction data concerning the E2 proteins arose from separated studies of E2-associated functions. In total, we referenced 57 cellular proteins binding to one or several E2 proteins by database mining (VirHostNet, VirusMINT and PubMed; Table ). Given the major public health concern of cervical cancers, prior studies have primarily focused on E2 proteins from the mucosal HR-HPV HPV16 and HPV18, creating a bias toward the most studied E2 proteins. Interactions were detected with the three domains of E2, although most of them concerned the E2 TAD (Fig. ). To get further insights into the E2-host interplay, we merged all known E2 interactions to build a network, and conducted an analysis based on the GO (Gene Ontology) terms associated with the E2-targeted cellular proteins. For that aim, we used the DAVID database [14], which generates a functional annotation of the data and allows the clustering of similar annotations into functional groups. 18 clusters have been identified within the list of published E2 interacting partners. Fig. () represents the parameters of each cluster: enrichment score, p-value and prevalence (percentage of GO-terms included in the cluster). By combining these parameters, four main functional families emerged presenting overall the highest criteria (high enrichment score, low p-value and high prevalence): transcription regulation, RNA processing, apoptosis and modulation of the cell cycle (Fig. ). This analysis provides insight into the range of E2 functions arising from E2-host interplay, but because of the limited number of genotypes represented, it gives only restricted information regarding HPV pathogenesis.

In order to have a more global outlook of E2 interactions, we recently performed a comparative interactomic assay with the E2 proteins from 12 genotypes representative of HPV diversity in terms of tropism (cutaneous and mucosal) and pathogenic potential (high risk or low risk) and spanning over three HPV clades (α, α and μ HPV) [15]. It consisted of a large-scale identification of E2-interaction partners by yeast two-hybrid, followed by validation of the interactions in mammalian cells with the 12 E2 proteins. 98 cellular proteins were found to interact with at least one of the E2 proteins, of which 76 were new cellular E2 partners. This approach provided the first overview of E2 activities across multiple genotypes. We will discuss below the E2 functions in relation to E2-host protein-protein interactions, in light with the recently performed large-scale interactomic study.

E2 and Transcription Regulation

The most represented category of E2 partners emerging from known interactions is Transcription, which appears in several clusters, involving both many GO term associations and numerous E2 partners (Fig. ). Not surprisingly, such prevailing targeting of cellular transcriptional machinery was reproduced in the interactome with the 12 E2 proteins, substantiating the prominent role of E2 as a transcriptional regulator.

E2 regulates HPV gene transcription either positively or negatively. The regulation of the early promoter transcription depends on the relative position of E2BS along the regulatory region. Within the genital HPV group, two E2BS overlap sequences of the early promoter. In such configuration, E2 acts as a repressor: by binding to these sites, E2 interferes with the binding of the cellular transcription factors TBP and SP1 by steric hindrance [16-18]. Conversely, in the LCR of cutaneous HPVs the E2BS are located farther from the early transcription initiation start, which leads to an overall transcriptional activation by E2 [19]. In addition, the stability of E2/DNA complexes differ according to the sequence of E2BS [16, 20, 21], and this has been proposed to modulate the effect of E2 on early promoter transcription [22].

The early promoter repression is not only due to steric hindrance, but is also the outcome of E2 binding to cellular proteins involved in the recruitment of repressive complexes on the LCR. For example, interaction of E2 with transcriptional regulatory proteins as the BRD4 protein, the histone demethylase SMCX, and the EP400 component of histone acetyltransferase complex have been shown to be required for efficient repression of the HPV18 LCR [23]. In line with these observations, we found in the comparative E2 interactome analysis a high proportion of targeted factors categorized as transcription repressors (Fig. ), highlighting the importance of repression in E2-mediated transcriptional regulation [15].

In addition, the E2 proteins are assumed to activate transcription from the late HPV promoter, which depends on epithelial differentiation and is not easily accessible experimentally. Consequently, little is known about the regulation of late genes by E2 in the context of viral life cycle. The transcriptional activation capacities of E2 were mostly studied using synthetic promoters containing multimerized E2BS upstream of a minimal promoter, or in in vitro transcription assays. It was determined that the HPV E2 proteins can activate promoters over large distances [24], and exhibit different transcriptional activation capacities. Notably the high risk E2 proteins were found to be intrinsically more potent transcriptional activators than the low risk E2 proteins [25]. These differences might result from different interplays between the E2 proteins and the host cell machinery.

Interactions detected with E2 over a wide series of studies identified a large spectrum of transcriptional targets, ranging from basic factors to chromatin-related factors (Fig. ), suggesting that the E2 proteins are able to interfere at multiple levels of transcriptional processes.

Proteins with transcriptional activation properties have been shown to act, at least partially, by manipulating stability and assembly of preinitiation complexes [26]. TFIID is a pivotal component of the cellular transcription initiation complex, and is composed of TBP and TBP-associated factors (TAF) [27]. Regulatory, sequence-specific transcription factors bound in promoter regions allow efficient recruitment of the preinitiation complex through co-activators, mediating protein-protein interactions between regulatory transcription factors and the basal transcription machinery. As depicted in Fig. (), several interactions have been detected between E2 and both basal (TAF, TBP, GTF2B) or regulatory (SP1, HOXC9, NR4A1, C/EBP) transcription factors, as well as with co-activators (TMF), indicating that transcriptional properties rely, at least partly, on modulation of preinitiation complex formation on specific promoters.

Other targets of E2 act at the chromatin level to regulate transcription. Indeed, the best-characterized partner necessary for E2-dependent transcription is BRD4, which binds to acetylated histones and stimulates RNA polymerase II-dependent transcription by recruiting distinct transcriptional regulators [28]. Mutation of the amino acids responsible for the interaction with BRD4 (notably I73 and R37 for HPV16 E2 [29]) strongly affects the transcriptional capacities of E2, making BRD4 an essential mediator of E2’s transcriptional properties [28].

In addition, the E2 proteins are able to modulate transcription through modification of the histone code, by targeting histone-modification factors. Notably, several histone acetyltransferases PCAF, EP300, EP400 or CBP were shown to interact with E2 [30-32]. Histone acetyltransferases function within macromolecular complexes, recruited to their target promoters by interacting with sequence specific transcription factors. Histone acetylation both loosens nucleosome-DNA interactions and promotes interaction of modified histones with other transcription activators proteins.

Lastly, the HPV E2 proteins interact with members of chromatin remodeling complexes, involved in the deposition or displacement of nucleosomes. Two members of chromatin remodeling complexes were known E2 targets, NAP1L1 and hSNF5 [33, 34]. NAP1L1 is primarily involved in replication-coupled nucleosome assembly by mediating the incorporation of histones H2A-H2B dimers in nucleosomes. Its direct binding to E2 from genotypes 5, 8 and 18 has been shown to enhance their transcriptional activation capacities [34], though this interaction might primarily impact on the replication activating functions of E2. On the other hand, the SWI/SNF chromatin remodeling complex was recently shown to enhance HPV18 E2-dependant transcription through direct binding of the core component hSNF5 to E2 [33]. Such interactions might counteract the repressive effect of chromatin through the recruitment of chromatin remodeling complexes to the LCR, in order to regulate transcription.

Numerous interactions thus emerged from studies performed over time, giving clues to some characteristics of E2 transcriptional functions, particularly the broad range of targeted transcriptional processes. However, as discussed earlier, it is difficult to determine to what extent these characteristics are conserved, given that they were described only with a restricted number of HPV genotypes.

The recent comparative mapping of the E2-host protein-protein interactions allowed the identification of 26 transcription regulators targeted by E2, showing a partial overlap with known E2 partners consisting of 8 proteins. It provides an experimental assessment of the interplay between E2 and the transcriptional machinery across multiple genotypes. We were surprised that in fact, if consistently predominant for all genotypes tested, the targeting of cellular transcription factors was quite heterogeneous among different E2 proteins (Fig. , ) [15]. Panel B in Fig. () shows the positive interactions engaged by the different E2 proteins and transcription factors, in correlation with the associated regulatory functions provided by their associated GO terms. Overall, our study corroborated the targeting of a large spectrum of transcription processes by all E2 proteins, which therefore likely represents a common trait of E2 transcriptional function. The range of E2 transcriptional targets encompasses transcription factors and cofactors, associated with transcriptional activation, repression, or possibly with both depending on the promoter context. This analysis highlights a prominent targeting of activators, representing more than half of the partners for all E2 proteins. As noted earlier, E2 also interacted with cellular proteins dedicated to transcriptional repression, and this targeting was more pronounced for the cutaneous β-type E2 proteins than for the α-type E2 proteins.

10 transcription factors are targeted by almost all E2 proteins, and BRD4 was consistently found among the strongest interactors [15]. These common targets may represent a “transcription core”, essential for the transcriptional functions of all E2 proteins. Of the common targets, only the basic transcription factor TBP binds to the HPV regulatory region. The other sequence-specific regulatory transcription factors of this core are TP53, MGA and NR4A1. TP53 primarily activates the transcription of cell cycle arrest or cell death genes [35]. MGA functions either as a repressor or as an activator of MYC-MAX target genes through binding to core E-box DNA sequences upon heterodimerization with the MAX transcription factor [36]. NRA41 (NUR77) is an Orphan nuclear receptor that may act by binding to NGFI-B response element in promoters of cellular genes [37]. Given that none of these factors harbor specific binding sites in HPV LCR, their interactions with the E2 proteins may primarily impact on the transcription of cellular genes rather than of viral genes. It thus emerges that part of E2 transcriptional functions might be involved in the modulation of expression of a common set of cellular genes. Several E2 proteins were indeed shown to affect the transcription of cellular genes [38-42], which could contribute to create cell conditions appropriate for the successive stages of the viral life cycle.

Other factors of the transcription core are co-activators such as MYST2, NMI or the GPS2 co-repressor. GPS2 is a subunit of the NCoR co-repressor complex [43] and its binding to E2 may be involved in the repression of HPV LCR, or affect the expression of NCoR1-repressed genes. NMI increases coactivator protein recruitment to sequence-specific transcription factors [44] and thus may be involved in E2-dependent transcriptional activation. MYST2 belongs to the histone H4-specific acetyltransferase complex HBO1, which acts as a coactivator of TP53-dependent transcription but has also been shown to specifically repress Androgen-Receptor (AR)-mediated transcription [45]. It is thought to regulate DNA replication through chromatin acetylation, and therefore its interaction with E2 may be involved in viral DNA replication, in addition to regulation of transcription.

The two remaining commonly targeted factors are the RSF1 and BAZ1A components of chromatin-remodeling complexes. Interestingly, both factors are part of chromatin remodeling complexes containing the SNF2H helicase. SNF2H is the catalytic subunit of four chromatin-remodeling complexes: CHRAC, RSF, ACF/WCRF, and WICH, which all contain other subunits regulating the specificity or catalytic activity of SNF2H. BAZ1A and RSF1 are accessory, non-catalytic subunits of three of the SNF2H-related complexes: ACF, CHRAC and RSF (Fig. ). The RSF chromatin remodeling complex is primarily involved in transcription regulation by mediating nucleosome deposition and positioning along regular arrays. It can either facilitate transcription or direct transcriptional repression by interacting with transcription factors, as was shown for the HBV pX protein [46]. Its interaction with a wide range of E2 proteins suggests that it may have a general role for E2 transcription functions, potentially both in repression and activation. RSF possesses not only chromatin remodeling activity but also chromatin assembly activity in vitro, and it has been proposed that RSF plays a role in silent chromatin formation [47]. Its recruitment by E2 could thus be involved in long-term silencing of HPV genome associated with viral persistence. Targeting of BAZ1A could have implications in E2-mediated transcriptional repression as well. Indeed, BAZ1A has been shown to interact with NCoR1 [48, 49]. E2 binding to both GPS2 and BAZ1A thus highlights a double targeting of the NCoR complex, which underlines the importance of NCoR-mediated repression for E2 transcriptional functions during infection. In addition, SNF2H forms the ACF complex with BAZ1A, which generates evenly spaced nucleosomes, creating a repressive state of chromatin [50, 51]. It is also involved in DNA replication, in particular it facilitates replication of heterochromatin regions [49]. Together with POLE3 and CHRAC subunits, BAZ1A and SNF2H assemble in the CHRAC complex, which also sets up repressive chromatin, and is involved in nucleosome assembly upon DNA replication [52, 53].

Overall, we propose that factors of the transcription core are instrumental for the transcriptional functions of all E2 proteins, and essential for common E2-dependent transcription processes as the timely regulation of viral gene expression in the course of the viral life cycle.

The remaining non-shared interactions are likely to reflect a degree of specificity in E2 intrinsic transcriptional properties, and may point to mechanistic variability in E2 transcriptional functions. For example, the comparative interactomic study reveals that 16E2 is the only one to bind the general transcription factor TFIIB (GTF2B), while it binds less efficiently to BRD4 when compared to the other E2 proteins [15] (Fig. ). The 16E2 protein has been shown in in-vitro studies to exhibit especially potent transcriptional activation properties, and this was proposed to rely on a more efficient recruitment of basal transcriptional factors than other E2 proteins [22]. The comparative interaction mapping substantiates this hypothesis and suggests it could be due to GTF2B binding. By contrast, BRD4 would not be expected to significantly contribute to the strength of HPV16 E2 transcriptional activity.

All other targeting specificity is associated with cutaneous HPV. Interestingly, most of the cutaneous-specific factors (ENO2, SFRS2, TOB) mediate transcription repression. Their preferred interaction with cutaneous HPV E2 proteins possibly underlies the decreased efficiency of E2-dependent transcription that was observed with synthetic promoters [15]. Conversely, the low level of transactivation is not related to reduced binding to BRD4, since the β-type cutaneous E2 proteins bind to high levels to BRD4. The activation properties of the cutaneous E2 proteins may thus result from a combined recruitment of repressors together with BRD4 and other activators.

The variety of cellular factors targeted by E2 provides an experimental appraisal of the mechanistic diversity of E2 transcriptional activity. It highlights an intricate interplay between the E2 proteins and the host cell transcriptional machinery, and allows the identification of a core set of common targets.

E2 and Viral DNA Replication

By binding to the LCR, E2 also activates the replication of viral DNA. Indeed, the viral origin of replication lies within the promoter-proximal region and contains E2BS as well as a binding site for the viral helicase E1. E2 both binds at high affinity to its cognate sites and to E1 through the TAD, consequently promoting the loading of E1 to the origin of replication ([54] for review). E2 also acts on viral DNA replication by recruiting the host cell replication factors TOP1 or RPA through direct interaction. Only few additional replication factors were identified in the comparative interactome, like POLDIP2, interacting with DNA polymerase δ and PCNA, or ORC3L, a protein involved in the assembly of the pre-replication complex.

E2 and Segregation of the Viral Genome

E2 is also critical for viral genome partitioning [4]. As such, E2 is required for the maintenance of viral episome in the basal dividing cells of the epithelium in the early steps of the viral life cycle. It is also essential to maintain a reservoir of HPV genomes during long-term infections, and is thus a crucial determinant of viral persistence. E2 acts as a bridge between viral episomes and host chromosomes, through binding to viral genome and protein-protein interactions with different factors of the mitotic apparatus. Initially, the targeting of E2 to mitotic chromatin has been demonstrated to rely on the binding of E2 TAD with the double bromodomain protein BRD4 [55]. It was later shown that, while the binding to BRD4 is conserved for all HPV E2 proteins and is consistently required for their transcriptional activation capacities, its role in the targeting of E2 to mitotic chromatin is variable [4]. In fact, several other factors of the mitotic apparatus have since been proposed as relevant for E2-mediated HPV genome partitioning, as the ChLR1 helicase or MKlp2 kinesin [4, 56, 57]. We could confirm the binding of MKlp2 to almost all the tested E2 proteins, and also identified the mitotic spindle-associated factor CEP350, which could be relevant for partitioning of the viral genome.

E2 and RNA Processing

It has been also shown that E2 could assist the splicing of mRNA. This is achieved for the cutaneous E2 proteins through interactions with SR factors, known to be key regulators of RNA processing mechanisms [58]. The mRNA processing machinery functions along with transcription, and some transcriptional activators are also able to interact with essential splicing factors [59]. RNA processing is a highly complex process, involving the formation of the spliceosome composed of mRNA substrates, snRNA (small nuclear RNA), snRNP (small ribonucleoproteins), and associated with splicing regulatory proteins, notably SR proteins.

The relation between E2 and RNA processing, contrary to other E2 functions, has first been uncovered with cutaneous HPVs. β-type E2 proteins harbor a long hinge enriched in arginine, serine and glycine residues. Arginine/serine (RS) dipeptide repeats are the characteristics of the SR superfamily of proteins, primarily involved in processing of mRNA. These observations led to the hypothesis that E2 proteins could share similar functions as cellular SR proteins or interact with them [58]. Accordingly, the HPV5 E2 protein was found to interact with a number of SR proteins [58, 60]. Among others, 5E2 interacts with SFRS1, best known as ASF/SF2 [60], an essential sequence specific splicing factor, as well as with SFSR2 (SC35) and SFRS7 (9G8). By binding to both its cognate DNA elements and to SR proteins, HPV5 E2 was shown to regulate the splicing of transcripts near the promoter, probably by increasing the local concentration of SR factors [58]. The mucosal α-type HPV E2 proteins appear to modulate mRNA splicing as well, in another way since 16E2 protein activates the expression of splicing factors such as ASF/SF2 (SFRS2) [40]. This function appears to be critical given that viral genomes expressing an E2 protein defective for SFRS1 upregulation show reduce levels of viral RNA [40] (reviewed in [61]). Also, 16E2 was recently shown to inhibit the polyadenylation of viral mRNA [62]; further substantiating the participation of E2 to viral mRNA processing.

Regulation of RNA by E2 might be even broader, since E2 interacts with cellular factors regulating RNA processing other than SR proteins. Indeed, SMN1, an RNA binding protein involved in the assembly of snRNP [63], interacts with mucosal HPVE2 proteins and enhances E2-mediated transcriptional activation [64]. For the cutaneous HPV5, E2 binds TNPO3 [60], a protein that mediates the nuclear import of the splicing factors SFRS1 and SFRS2 by recognizing phosphorylated RS domains [65]. These interactions suggest that the regulation of RNA processing by the E2 proteins is extensive and thus must be pivotal for regulation of HPV gene expression.

In the comparative interactomic study conducted with the 12 E2 proteins, we corroborated interactions of E2 proteins with SFRS1, 2 and 7 proteins of the SR family [15]. The interaction data set provides evidence that the targeting of SR splicing factors is not restricted to the β-type E2 proteins, but also extends to the α-type E2 proteins, and therefore is conserved among all HPV genotypes. Interaction with SFRS7 was the strongest for both the cutaneous and the mucosal E2 protein, while SFRS2 only interacted with the E2 proteins from β/μ cutaneous HPV, suggesting that variations may occur in the regulation of mRNA splicing according to HPV tropism. The α-type HPV E2 proteins exhibit greatly reduced interaction levels, in accordance with the presence of short R-alternating sequences in their hinge regions. This indicates that E2 proteins encoded by the genital α-HPVs are able to interact with SR proteins but less efficiently than those from the cutaneous HPVs, which might explain why it had not been detected. A conserved role of E2 in the regulation of RNA splicing through interaction with SR-splicing factors nevertheless emerges, probably involved in the synthesis of viral mRNA, as initially proposed for HPV5 E2 [58]. It substantiates a pivotal role of E2 in the timely regulation of mRNA production through alternative splicing, enabling the sequential expression of viral proteins along the viral life cycle [40]. SRFS1 and 7 are the strongest of all the interactors of cutaneous E2 proteins, suggesting that binding to splicing factors is prevailing for the viral gene expression by these E2. For the α-type E2 proteins, the interaction with splicing factors is not predominant and might thus be less important. Altogether, the E2 proteins nevertheless emerge as primary factors to orchestrate viral regulation, both at the level of transcription and of mRNA splicing.

E2 and Apoptosis

The pro-apoptotic activity of E2 is one of the first described E2 functions independent of its binding to the viral genome (reviewed in [6]). A striking aspect of the involvement of E2 in apoptosis is that it was demonstrated to be specific for the HR-HPV E2 proteins [8, 9], first supporting the idea that E2 proteins might have developed specific activities that correlate with the HPV oncogenic power. The two main pro-apoptotic pathways are dependent on caspase activation. On one hand, the intrinsic pathway can be induced by p53 and involves mitochondrial dysfunctions leading to the activation of caspase 9. On the other hand, the extrinsic pathway is dependent on death receptors signaling and activation of caspase 8. Both pathways eventually merge in the activation of effector caspases, which leads to subsequent cell death. Several interactions between HR-HPV E2 and effectors of both the extrinsic and the intrinsic pathways have been identified over the years, notably with p53, CASP8 and C-FLAR [66-68]. Overall, the interplay between HR-HPV E2 proteins and the apoptotic machinery is complex, as reflected by the variety of interactions identified (see [13]). Other known interactions may have implications for E2-induced apoptosis. For example, SMN1 was shown by Strasswimmer and colleagues to interact with E2 [64]. Mutation in this protein results in a neurodegenerative disease called SMA (Spinal Muscular Atrophy) and recent evidence demonstrated that knock down of SMN1 increases neuronal cell death which might be counteracted by over-expression of the anti-apoptotic Bcl-xL [69]. We speculate that binding of SMN1 to E2 might thus favor apoptosis. A yeast-two hybrid screen also identified NR4A1 as an interacting partner of HPV16 E2 [70]. This protein is a member of the nuclear-receptor superfamily and plays a role in regulating cell growth and apoptosis [71]. NR4A1 triggers cytochrome c release from the mitochondria [72] both in the context of cancer cells and of virus-induced apoptosis [73, 74] by interacting with Bcl2 [75]. This interaction could thus modulate the intrinsic apoptotic pathway and this reinforces the notion that E2 uses various redundant approaches to regulate cell death. One of the final steps of apoptosis following caspase activation, is the alteration of DNA structures (strand breaks, oxidized bases). TOP1 is known to be involved in modification of the topology of supercoiled DNA by decreasing the torsional stress on DNA [76]. TOP1 has recently been linked to apoptosis since down regulation of its expression by RNA interference reduces apoptotic DNA fragmentation [77]. The interaction with E2 [78] might stabilize TOP1 and favor apoptosis completion, suggesting that E2 could impact on multiple steps of apoptosis.

In the large-scale study conducted in our lab, regulation of apoptosis emerged as a functional family targeted by E2, which indicates that the 12 studied E2 proteins have an intrinsic capacity to interact with apoptosis regulators. Thus, the regulation of cell death and survival pathways is not restricted to the mucosal HR-HPV but extends to the LR-mucosal as well as to the cutaneous HPVs. Our data highlight the targeting of both positive and negative regulators of apoptosis, suggesting that the E2 proteins may interfere oppositely with cell death pathways at different steps of the viral life cycle. However, for some of the targets, the effect on cell death can depend on the context (level of damage for TP53) or the protein isoform (BCL2L1), so that the final outcome of the interaction may vary according to cellular conditions, and is difficult to anticipate from this complex picture.

Among the 12 proteins of this family, three were bound by all E2 proteins (CASP8, TAX1BP1 and TP53). For TP53 and CASP8, the binding of E2 may not have similar functional consequences according to HPV genotype. The detection of an interaction between p53 and all the E2 proteins was unexpected based upon previous results describing an interaction of p53 only with the HR-mucosal E2 proteins [9]. From our recent data, we infer that E2 binding to p53 is actually conserved among HPV genotypes, but might have a different impact on p53 function and induce apoptosis only in the case of HR-HPV E2 proteins. For the LR-HPV E2 proteins, binding to p53 may alter other activities, perhaps the induction of a cell cycle arrest. This notion is further supported by the detection of interactions between Caspase 8 and all the E2 proteins tested. Yet, it was shown that among the mucosal HPV, only HR-HPVE2 proteins are able to induce apoptosis through direct binding to caspase 8 [8, 66]. E2-mediated caspase 8 activation is linked to the accumulation of E2 protein in the cytoplasm due to an active nucleo-cytoplasmic shuttling, which does not occur in the case of LR-mucosal HPV [8]. However, the comparative interaction mapping provides evidence for a conserved E2 binding to caspase 8 across all HPV genotypes, suggesting that the outcome of the interaction between the LR-HPV E2 proteins and caspase 8 might be different from that of HR-HPV E2 proteins. It could notably play a role in keratinocyte differentiation, since a role of caspase 8 in skin homeostasis has been recently recognized [79]. For both p53 and caspase 8, the functional impact of interactions with the non-apoptotic E2 proteins could be the opposite of interactions with the proapoptotic E2 proteins, i.e. interfering with death induction. This binding could thus be differentially involved in viral pathogenesis according to HPV oncogenicity. Another hypothesis is that these interactions take place at different times of the HPV life cycle. It can be envisioned that the LR-HPV E2 induces apoptosis only at late stages of the viral life cycle in the upper layers of infected epithelium, in which case it would have escaped detection in cell culture conditions. Indeed, the LR-HPV E2 proteins could trigger cell death only after nuclear membrane breakdown which occurs during keratinocyte differentiation, and would allow the accumulation of E2-Caspase 8 complexes in the cytoplasm enabling caspase activation. Such delayed apoptotic induction could be involved in common viral processes such as viral particles release.

Overall, the functional targeting of apoptosis regulatory factors shows that the E2 proteins actively take part in the manipulation of cell death or survival pathways, tightly orchestrated by the virus throughout infection.

E2 and the Proteasome Ubiquitin System

Controlling or rewiring the ubiquitin-proteasome pathway by viral proteins is likely to have a great impact on the host cell. The consequences could be multiple: degradation of antiviral molecules, and diversion of ubiquitin ligase substrate specificity. There are two types of E3 ubiquitin ligases, the HECT (Homologous to E6AP C-terminus) proteins with an intrinsic catalytic activity, and the Cullin-based complexes composed of a cullin, a RING finger protein, and one or more substrate specificity adapters [80-82].

Several interactions have been shown between the E2 proteins and elements of the ubiquitin proteasome system over the past years. E2 proteins from HR-HPV16 and 18 interact with CDC20 and FZR1 (also known as cdh1), two substrate recognition subunits of the ubiquitin ligase complex APC/C (Anaphase Promoting Complex, Fig. ), while in contrast the LR-HPV6 and 11 E2 proteins are not able to engage in such interaction [7]. The APC/C complex is responsible for anaphase transition through the degradation of specific substrates such as cyclin B. The binding of HR-HPVE2 proteins to CDC20 and FZR1 has been shown to inhibit APC/C-mediated cyclin B degradation, which leads to a metaphase arrest and genomic instability [7]. Another APC/C substrate is SKP2 [83], a substrate recognition subunit of the ubiquitin ligase complex SCF (Fig. ), and its stabilization by E2 activates this E3 ligase complex. Since SCF is responsible for the degradation of the cell cycle negative regulators p21 and p27, this activation results in an enhanced G1 to S transition [84]. Therefore it appears that identified interactions between the HR-HPVE2 proteins and APC/C ubiquitin ligase complex strongly impact on cell cycle. In addition it has been demonstrated by Bellanger et al. that the complex SCF SKP2 can degrade HPV18 E2 protein itself, suggesting a complex and controlled feedback mechanism [85].

Few interactions with proteins involved in ubiquitination mechanisms have also been shown to be involved in the regulation of E2 stability itself. The E2 proteins have short half-lives, and this rapid turnover has been shown to be the result of ubiquitin-mediated proteasome degradation targeting of the N-terminal domain [8, 86-88]. Identification of the ubiquitin ligases responsible for the degradation of some E2 emerged only recently: HPV18 E2 can be degraded by the ubiquitin ligase complex SCFSKP2, and this degradation is mediated by direct interaction with the adaptor SKP2 [85] while for HPV16, it was demonstrated that a cullin 3-based complex (BRC3) mediates ubiquitination and subsequent degradation [88]. It was hypothesized that 16E2 is recruited to this complex through BTB domain proteins, which are substrate adaptors of cullin 3-based complexes (Fig. ) [88, 89]. In line with this hypothesis, two BTB domain–containing proteins, BTBD1 and BTBD2, were identified as partners of HPV16E2 protein in a yeast two-hybrid screen [90]. The process of E2 protein degradation is even more complex, since several cellular proteins have been shown to modulate in E2 stability as TAX1BP1 or BRD4 [87, 91, 92]. Moreover, since HPVs infect different tissue niches, it is unlikely that all E2 proteins are degraded in the same way, and to our knowledge, there are no data available on the stability of E2 proteins from cutaneous β- or μ-HPV types.

Previous studies thus converge on the targeting by the HR-HPV E2 of cullin-based E3 ligase complexes mediated by interactions with substrate adaptors that ensure the specificity of the protein to be degraded. The comparative interaction mapping reinforces the notion of a preferred targeting of E3 ubiquitin ligases by the E2 proteins. Indeed, of 16 factors of the ubiquitin-proteasome pathway found to interact with E2, eight are part of E3 ubiquitin ligase complexes. In line with previously published works, most of them correspond to substrate-adaptors of cullin-based complexes, of which BTBD1 and SPOP are BTB-proteins forming complexes with cullin3, FBX022 is a F-box protein subunit of cullin1-based complex, CDC20 and CDH1 are part of the APC/C. The HR-specificity of E2 binding to CDC20 was recovered among the mucosal HPV, but this binding was also detected with cutaneous E2 proteins, suggesting that the involvement of this interaction in HPV pathogenesis differs according to HPV tropism. Of note, the SPOP adapter binds all E2 proteins, with reduced efficiency for HPV16 E2. Conversely, the HPV16 E2 protein was the only one to bind another adapter of cullin3-complexes, BTBD1, suggesting that the targeting of E2 to cullin3-based complex is conserved but could be mediated through diverse interactions with adaptors [15].

In addition, two HECT domain family proteins HUWE1 and WWP2 were identified as novel E2 partners, even further broadening the potential impact of E2 through interaction with ubiquitin ligases. This is compatible with the notion that E2 could not only modify the action of ubiquitin ligases toward their natural substrates, as was shown for cyclin B or SKP2, but also could divert their substrate specificity toward new substrates, the E2 proteins then acting as novel substrate adaptors. In that case, the targeting of ubiquitin ligase complexes by E2 would have a similar outcome than the binding of E6 to E6AP, which induces the hijacking of E6AP activity toward ubiquitylation of p53 [93]. Such potential E2-dependent substrates are expected to bind to E2 only when in complex with the Ubiquitin ligase and are therefore not accessible through binary protein-protein interaction mapping.

Intracellular Transport

Among known E2-interacting proteins, several are assigned with GO-terms related to intracellular protein transport, such as KIF20A, AP3D1, TNPO3 or TMF1. AP3D1 is part of the complex AP-3, which regulates sorting of membrane proteins [94], and KIF20A acts as a motor for the transport of vesicle and has been shown to associate with a protein from the Rab family, RAB6A [95]. TMF1, an interactor of HPV16 E2, was also demonstrated to be involved in Rab-6 dependent transport processes [96]. Other members of the Rab family have been implicated in the entry pathway of certain types of HPV such HPV31 [97]. Additionally, E2 was shown to interact with TNPO3, involved in the entry steps of other viruses such as HIV-1 [98]. However, the targeted factors are not numerous enough to represent an enriched functional family based on the criteria of GO-term analysis (Fig. ). The large-scale interactomic strategy performed with the 12 E2 proteins identified additional proteins involved in intracellular trafficking, which now emerges as a functional family targeted by E2 [15]. Indeed, 15 cellular E2 targets were classified in intracellular transport family, among which 13 are newly identified E2 partners [15]. Notably, the targeted proteins are more concentrated on vesicular transport, affecting dynamics and maintenance of intracellular membranous organelles (Clathrin, Rab-family protein, lysosomal transport, vesicle transport to and from the golgi apparatus). Even if interactions of the E2 proteins with such factors is surprising when considering that E2 is primarily a nuclear protein, it might be hypothesized that they play a role in the translocation of HPV episomes to the nucleus. Indeed, the regulation of vesicle transport is important for entry processes during HPV infection, and we noticed that most of E2 targets overlap with the pathways of HPV entry in the cell as depicted in Fig. (). E2 is able to link the viral DNA and it interacts with the minor capsid protein L2 [99], which is pivotal for the trafficking of viral genomes from the capsid to the nucleus [100]. It is thus tempting to speculate that E2 participates in the delivery and transport of viral genomes through the cytoplasm by interacting with the identified proteins, and in collaboration with L2. This would imply that E2 is encapsidated in the virion bound to the viral genome, and that it plays a role at very early steps of viral infection. Pseudovirions systems have been developed and they indicate that E2 is not necessary for proper delivery of DNA materials to the nucleus, in contrast to L2 [101]. However, only sparse information is available on the involvement of E2 in pseudovirions infectivity, and these studies were only done with BPV1 pseudvirions. In the BPV1 system, a study reported that E2 enhanced encapsidation of full-length viral DNA and may be packaged with the pseudovirion [102]. These findings, however, were not corroborated in another study showing that E2 expression does not alter BPV1 pseudovirions production and infectivity [101]. In light of the emergence of the functional targeting of intracellular transport pathways, we feel that this issue should be re-evaluated with HPV pseudovirions. We wonder whether the use of over-optimized systems for pseudovirion production (codon optimized ORF) could hide a potential involvement of E2 under normal conditions. Also, it should be taken into consideration that the pseudovirion systems allow the efficient delivery only for small DNA molecule. A potential role of E2 could be uncovered with viral genomic DNA containing E2 binding sites, thereby mimicking authentic HPV genome encapsidation and translocation.

Overall, this compilation of observations raises the intriguing possibility that E2 could take an active part in very early steps of HPV infection, which we feel is worth investigating.

E2 and Keratinocyte Migration/Differentiation

Several lines of evidence point to a modulation of cellular genes involved in keratinocyte migration and differentiation by the E2 proteins. Indeed, E2 represses the transcription of ITGB4, a major integrin for keratinocyte attachment to the matrix, and leads to the detachment of keratinocyte from the underlying matrix [38]. E2 also activates MMP9, a protein involved in cell migration, invasion and metastasis [41]. In addition, E2 was shown to alter the transcription of differentiation markers through cooperation with cellular transcription factors such as C/EBP [103]. Accordingly, in a microarray study, E2 was shown to modify the expression profile of cellular genes, among which many were involved in cell differentiation [104]. It was also observed that E2 induces phenotypical changes typical of terminally differentiated cells, with increased expression of differentiation markers [105]. E2 therefore promotes differentiation by acting as a transcription factor which modulates the expression of a set of cellular genes.

More surprisingly, the regulation of cell differentiation by E2 might be mediated through direct protein-protein interaction. Indeed, it was recently published that E2 stimulates TNF-induced NFκB activation through direct interaction with TRAF5, an intermediate of the NFκB signaling pathway [5]. This could be involved in cell commitment to differentiation, since activation of NFκB by TNF induces an epithelial differentiation program [106]. Other known cellular proteins interacting with E2 have been implicated in cell differentiation or migration, such as CCHCR1, which affects the balance between proliferation and differentiation in keratinocytes [107, 108], or GNB2L1 reported to be involved in migration of carcinoma cells [109]. The interaction of these proteins with E2 in the context of infection might impact on their stability, localization and function, and thus promote the detachment of keratinocytes from the basal lamina and further induction of differentiation.

In the comparative E2 interactomics study, we detected a targeting concentrated around proteins involved in cellular adhesion. Interactions were identified with PVRL1, a protein involved in cell-to-cell junctions [110], PTK2B, a member of the focal adhesion kinase family [111] or FN1, the fibronectin responsible for the attachment to the extracellular matrix [112]. Overall, this reinforces the idea that E2 might impact on the cohesion of the infected epithelium.

We also uncovered a direct interaction of all E2 proteins with ITGB4, which as discussed before, is also transcriptionally repressed by E2 [38]. Such redundant targeting may highlight the importance of inactivating ITGB4 for progression of the viral life cycle. It might provide an early signal to induce differentiation, which requires migration of keratinocyte from basal to upper layers of the epithelia. Interestingly, ITGB4 interaction with E2 is conserved across all HPV genotypes tested, and ITGB4 was also identified as a target of E7 [138]. Theses observations reinforce the notion that the targeting of ITGB4 might be pivotal, and that its sequential regulation by E2 and E7 potentially essential for progression of the viral life cycle.

Overall, these observations point to a potential redundant effect of E2 on keratinocyte differentiation and migration through modulation of gene expression as well as direct targeting of cell-cell and cell-matrix contacts.

CONCLUDING REMARKS

Viral proteins are platforms of interaction that are known to orchestrate host proteins localization and degradation, control their expression levels and post-translational modifications, and rewire signal transduction in the infected cell [113, 114]. Recently, through improvement of detection technologies [115, 116], significant efforts have been made to generate comprehensive protein-protein interaction network maps in order to improve the understanding of virus-host relationship [117-119]. Additional steps should be focused on the development of integrative comparative interactomic approaches in order to decipher the interaction profiles of viral proteins from different pathogenic viruses. It should give insight into how genomic variations drive phenotypic divergence, but also on how viral proteins diversely hijack host cellular structures and pathways.

In the case of HPVs, a great divergence exists between infections with the different viral strains, ranging from variations in the tropisms (oral, genital, cutaneous epithelia) to variations in the pathogenic potential. The variability of HPV-associated lesions likely results from a diverse interplay between viral and host proteins. However, only recent studies have started to encompass large-scale analyses, which will reveal their wide impact on the comprehension of HPV pathogenesis.

Functional analysis of the E2 interactome provides a complete picture of the range of E2 functions, viewed through its interactions with the host proteome. It highlights a prevailing impact of E2 on the cellular processes of transcription, RNA processing, apoptosis, ubiquitination, and intracellular transport. The emergence of this latter unexpected functional family raises the provocative hypothesis that E2 may be involved in early steps of viral infection. In addition, a constant targeting of positive and negative regulators of the same cell processes emerges, indicating a duality in E2 functions, potentially driving opposing activities during the different steps of the viral life cycle. We did not identify any specific cellular functional targeting according to tropism or pathogenic power, suggesting that modulation of these cell processes takes part in the general regulation of the viral life cycle of all HPVs. However, this common functional targeting proved to be achieved through diversified patterns of interactions. Such variations in the interactions profiles most probably drive the contribution of E2 proteins to HPV pathogenesis.

Furthermore, the assessment of previously known interactions with the 12 E2 indicates that interaction data cannot be generalized from functional studies addressing only a subset of E2 proteins. For example, the binding of GTF2B had been detected for both BPV1 E2 and 16 E2, and was assumed to be conserved and involved in the activity of all E2 proteins. However, binding to this factor turned out to be specific for 16E2 and as such may account for some particular aspect of its transcriptional function as discussed earlier. By contrast, other interactions, as with TP53 or caspase 8, had been detected in connection to E2 functions specific for a subset of genotypes, but arise in the comparative interaction mapping as conserved across all HPV genotypes. These observations underline the importance of combining functional studies and interaction mapping to decipher E2 activities in relation with HPV pathogenesis.

Studying viral-host interactions has become an important challenge since targeting protein-protein interaction is now considered as an attractive approach for the development of new therapeutic strategies [120]. The use of small peptide based strategies to counteract specific interaction is studied for long [121]. But the ultimate option to develop an efficient treatment would be to identify a pan-HPV inhibitor. To that end, efforts should be put into the expansion of systematic and comparative approaches to identify protein-protein interactions. The targeting of E2 interactions would be a valuable therapeutic strategy, since E2 is an essential factor of the life cycle of all HPVs, is expressed early in the infection process, and is required for viral persistence.

Table 1.

Literature Curated Interactions. List of Interactions Found for the HPV E2 Proteins in the VirHostNet, virusMINT and PubMed Databases

ENSembl ID	Official Gene Symbol	Name	Synonyms	HPV Type	Refs.
ENSG00000065000	AP3D1	adaptor-related protein complex 3, delta 1 subunit	ADTD	HPV16	[123]
ENSG00000012048	BRCA1	breast cancer 1, early onset	BRCC1, PPP1R53, RNF53	HPV18	[124]
ENSG00000141867	BRD4	bromodomain containing 4	CAP, HUNK1, HUNKI, MCAP	HPV16,18,31,11,6,8,1	[55, 125, 126]
ENSG00000064726	BTBD1	BTB (POZ) domain containing 1		HPV16	[90]
ENSG00000133243	BTBD2	BTB (POZ) domain containing 2		HPV16	[90]
ENSG00000108561	C1QBP	complement component 1, q subcomponent binding protein	gC1Q-R, gC1qR, HABP1, p32, SF2p32	HPV5	[60]
ENSG00000064012	CASP8	caspase 8, apoptosis-related cysteine peptidase	Casp-8, FLICE, MACH, MCH5	HPV16,18	[66]
ENSG00000206355	CCHCR1	coiled-coil alpha-helical rod protein 1	C6orf18, HCR	HPV16	[90]
ENSG00000117399	CDC20	cell division cycle 20 homolog	CDC20A, p55CDC	HPV16,18	[7]
ENSG00000245848	CEBPA	CCAAT/enhancer binding protein (C/EBP), alpha	C/EBP-alpha, CEBP	HPV16,18,8	[103]
ENSG00000172216	CEBPB	CCAAT/enhancer binding protein (C/EBP), beta	C/EBP-beta, CRP2, IL6DBP, LAP, NFIL6, TCF5	HPV16,18,8	[103]
ENSG00000003402	CFLAR	CASP8 and FADD-like apoptosis regulator	c-FLIP, CASH, CASP8AP1, Casper, CLARP, FLAME, FLIP, I-FLICE, MRIT	HPV16	[68]
ENSG00000160917	CPSF4	cleavage and polyadenylation specific factor 4	CPSF30, NAR	HPV16	[62]
ENSG00000005339	CREBBP	CREB binding protein	CBP, KAT3A, RSTS, RTS	HPV18	[31]
ENSG00000036257	CUL3	cullin 3		HPV16	[88]
ENSG00000013573	DDX11	DEAD/H (Asp-Glu-Ala-Asp/His) box polypeptide 11	CHL1, CHLR1, KRG2	HPV16,11	[56]
ENSG00000242372	EIF6	eukaryotic translation initiation factor 6	b(2)gcn, EIF3A, ITGB4BP, p27BBP	HPV16	[123]
ENSG00000100393	EP300	E1A binding protein p300	KAT3B, p300	HPV18,5,8	[32]
ENSG00000183495	EP400	E1A binding protein p400	CAGH32, DKFZP434I225, KIAA1498, KIAA1818, P400, TNRC12	HPV16	[23]
ENSG00000105325	FZR1	fizzy/cell division cycle 20 related 1	CDC20C, CDH1, FZR, FZR2, HCDH, HCDH1, KIAA1242	HPV16,18	[7]
ENSG00000204628	GNB2L1	guanine nucleotide binding protein (G protein), beta polypeptide 2-like 1	Gnb2-rs1, H12.3, RACK1	HPV16	[123]
ENSG00000132522	GPS2	G protein pathway suppressor 2		HPV16	[127]
ENSG00000137947	GTF2B	general transcription factor IIB	TFIIB	HPV16	[128]
ENSG00000180806	HoxC9	homeobox C9	HOX3, HOX3B	HPV16	[90]
ENSG00000114166	KAT2B	K(lysine) acetyltransferase 2B	GCN5, GCN5L, P/CAF, PCAF	HPV16,18,11,6	[30]
ENSG00000126012	KDM5C	lysine (K)-specific demethylase 5C	DXS1272E, JARID1C, SMCX, XE169	HPV16	[23]
ENSG00000112984	KIF20A	kinesin family member 20A	RAB6KIFL	HPV16,11	[57]
ENSG00000135679	MDM2	Mdm2 p53 binding protein homolog	HDM2, HDMX, MGC5370	HPV16	[129]
ENSG00000187109	NAP1L1	nucleosome assembly protein 1-like 1	MGC23410, MGC8688, NAP1, NAP1L, NRP	HPV18,5,8	[34]
ENSG00000123358	NR4A1	nuclear receptor subfamily 4, group A, member 1	GFRP1, HMR, N10, NAK-1, NGFIB, NUR77, TR3	HPV16	[70]
ENSG00000104904	OAZ1	ornithine decarboxylase antizyme 1	AZI, MGC138338, OAZ	HPV16	[123]
ENSG00000143799	PARP1	poly (ADP-ribose) polymerase 1	ADPRT, PARP, PPOL	HPV18	[130]
ENSG00000166851	PLK1	polo-like kinase 1	PLK	HPV5	[131]
ENSG00000132963	POMP	proteasome maturation protein	C13orf12, HSPC014, UMP1	HPV16	[123]
ENSG00000214517	PPME1	protein phosphatase methylesterase 1	PME-1	HPV16	[123]
ENSG00000168066	SF1	splicing factor 1	ZFM1, ZNF162	HPV16	[90]
ENSG00000136450	SFRS1	serine/arginine-rich splicing factor 1	ASF, MGC5228, SF2, SF2p33, SFRS1, SRp30a	HPV5	[60]
ENSG00000161547	SFRS2	serine/arginine-rich splicing factor 2	PR264, SC-35, SC35, SFRS2, SFRS2A	HPV5	[60]
ENSG00000115875	SFRS7	serine/arginine-rich splicing factor 7	9G8, AAG3, HSSG1, RBM37, SFRS7, ZCCHC20, ZCRB2	HPV5	[60]
ENSG00000145604	SKP2	S-phase kinase-associated protein 2 (p45)	FBL1, FBXL1	HPV18	[85]
ENSG00000099956	SMARCB1	SWI/SNF related, matrix associated, actin dependent regulator of chromatin, subfamily b, member 1	BAF47, hSNFS, Ini1, RDT, Sfh1p, SNF5L1, Snr1	HPV18	[33]
ENSG00000172062	SMN1	survival of motor neuron 1, telomeric	BCD541, SMA, SMA1, SMA2, SMA3, SMA@, SMNT	HPV16,18,11	[64]
ENSG00000104852	SNRNP70	small nuclear ribonucleoprotein 70kDa (U1)	RNPU1Z, RPU1, Snp1, SNRP70, U1-70K	HPV5	[58]
ENSG00000185591	SP1	Sp1 transcription factor		HPV18,8	[39. 132]
ENSG00000147133	TAF1	TAF1 RNA polymerase II, TATA box binding protein (TBP)-associated factor, 250kDa	BA2R, CCG1, CCGS, DYT3, DYT3/TAF1, KAT4, NSCL2, TAF2A, TAFII250	HPV16,18	[133]
ENSG00000106290	TAF6	TAF6 RNA polymerase II, TATA box binding protein (TBP)-associated factor, 80kDa	MGC:8964, TAF2E, TAFII70, TAFII80, TAFII85	HPV18	[134]
ENSG00000178913	TAF7	TAF7 RNA polymerase II, TATA box binding protein (TBP)-associated factor, 55kDa	TAF2F, TAFII55	HPV5	[135]
ENSG00000106052	TAX1BP1	Tax1 (human T-cell leukemia virus type I) binding protein 1	CALCOCO3, TXBP151	HPV16,18	[91]
ENSG00000112592	TBP	TATA box binding protein	GTF2D1, SCA17, TFIID	HPV16,8	[135]
ENSG00000144747	TMF1	TATA element modulatory factor 1	ARA160, TMF	HPV16	[90]
ENSG00000064419	TNPO3	transportin 3	IPO12, MTR10A, TRN-SR, TRN-SR2	HPV5	[60]
ENSG00000198900	TOP1	topoisomerase (DNA) I		HPV16	[78]
ENSG00000163781	TOPBP1	topoisomerase (DNA) II binding protein 1	KIAA0259, TOP2BP1	HPV16	[136]
ENSG00000141510	TP53	tumor protein p53	LFS1, p53	HPV16,18	[9,137]
ENSG00000136527	TRA2B	transformer 2 beta homolog	Htra2-beta, SFRS10	HPV5	[60]
ENSG00000082512	TRAF5	TNF receptor-associated factor 5	RNF84	HPV1,3,5,6,8,9,11,16,18,33,39	[5]
ENSG00000175104	TRAF6	TNF receptor-associated factor 6, E3 ubiquitin protein ligase	RNF85	HPV1,3,5,6,8,9,11,16,18,33,39	[5]