Tarik Gheit- Infections and Cancer Biology Group, International Agency for Research on Cancer (IARC), Lyon, France
Papillomaviridae is a family of small non-enveloped icosahedral viruses with double-stranded circular DNA. More than 200 different human papillomaviruses (HPVs) have been listed so far. Based on epidemiological data, a subgroup of alphapapillomaviruses (alpha HPVs) was referred to as high-risk (HR) HPV types. HR HPVs are the etiological agents of anogenital cancer and a subset of head and neck cancers. The cutaneous HPV types, mainly from beta and gamma genera, are widely present on the surface of the skin in the general population. However, there is growing evidence of an etiological role of betapapillomaviruses (beta HPVs) in non-melanoma skin cancer (NMSC), together with ultraviolet (UV) radiation. Studies performed on mucosal HR HPV types, such as 16 and 18, showed that both oncoproteins E6 and E7 play a key role in cervical cancer by altering pathways involved in the host immune response to establish a persistent infection and by promoting cellular transformation. Continuous expression of E6 and E7 of mucosal HR HPV types is essential to initiate and to maintain the cellular transformation process, whereas expression of E6 and E7 of cutaneous HPV types is not required for the maintenance of the skin cancer phenotype. Beta HPV types appear to play a role in the initiation of skin carcinogenesis, by exacerbating the accumulation of UV radiation-induced DNA breaks and somatic mutations (the hit-and-run mechanism), and they would therefore act as facilitators rather than direct actors in NMSC. In this review, the natural history of HPV infection and the transforming properties of various HPV genera will be described, with a particular focus on describing the state of knowledge about the role of cutaneous HPV types in NMSC.
Introduction
Papillomaviridae is a family of small non-enveloped icosahedral viruses with double-stranded circular DNA, which range in length from 5,748 bp for Sparus aurata papillomavirus 1 (SaPV1) to 8,607 bp for canine papillomavirus type 1 (CPV1). Papillomaviruses (PVs) infect basal keratinocytes of the mucosal and cutaneous epithelia of both animals (reptiles, birds, marsupials, and others) and humans (1). In 2016, the first characterization of a PV in fish (SaPV1) rendered this family of viruses much older than expected, with an emergence 450 million years ago (2). PVs are considered to be host-restricted; however, in rare cases cross-species transmission may also occur (3, 4).
Based on the nucleotide sequences of the major capsid protein L1, the PV study group within the International Committee on Taxonomy of Viruses (http://ictv.global/report/papillomaviridae) has classified the PVs into 53 genera, among which only 5 genera include PVs that infect humans (HPVs). To be classified as a novel PV type, the nucleotide sequence of L1 must share < 90% similarity with other PVs (5). During the past decade, there has been an exponential increase in the identification of new betapapillomaviruses (beta HPVs) and gammapapillomaviruses (gamma HPVs) that have been discovered with the advent of new technologies such as next-generation sequencing. More than 200 different HPV types have been listed by the International HPV Reference Center (www.hpvcenter.se) (6), and this number continues to expand (7).
HPV types are organized into five major genera: alpha, beta, gamma, mu, and nu (5). The genus gamma includes the majority of the known HPVs, with 99 types, followed by the genera alpha (n = 65) and beta (n = 54). The genera mu and nu include only 3 and 1 types, respectively (www.hpvcenter.se, on 2019-01-30). A recent study reported the identification of the complete genome of 83 previously unknown HPV types, among which 69 were classified as gamma (8).
Based on epidemiological and biological data, a subgroup of 12 mucosal HPV types classified in the genus Alphapapillomavirus (alpha HPV), i.e., HPV16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, and 59, referred to as high-risk (HR) HPV types, has been classified as carcinogenic (IARC Group 1), and eight other HPV types, i.e., HPV26, 53, 66, 67, 68, 70, 73, and 82, have been classified as probably or possibly carcinogenic (IARC Groups 2A and 2B) (9, 10). The HR HPV types are the etiological agents of several cancers, such as those of the cervix, vagina, vulva, anus, penis, and a subset of head and neck cancers (HNCs), particularly oropharyngeal cancer (11). In 2012, 15.4% of cancers worldwide were attributed to carcinogenic infections (12). Together, the HPV-related cancers represent 630,000 new cancer cases per year, which account for ~30% of all cancers induced by infectious agents (10). Epidemiological and biological studies showed that HPV16 is the most oncogenic of the HR HPVs. The genus alpha also includes low-risk (LR) HPV types that cause benign lesions. HPV6 and 11, the most studied LR HPV types, induce benign genital warts or condylomata acuminata (13). These HPV types are also found in recurrent respiratory papillomatosis (14), in children. The genus alpha also includes a few cutaneous HPV types (HPV2, 3, 7, 10, 27, 28, and 57), which cause common and plantar warts (15–18).
The genera beta, gamma, mu, and nu contain HPVs that infect cutaneous epithelia. The genus alpha also contains some HPV types, such as HPV2, 3, and 10, that colonize the human skin. The mu and nu genera include only a few HPV types, whereas the genus beta includes more than 54 HPV types, which are subdivided into five species (β1–5), and the genus gamma includes 98 types subdivided into 27 species (6). In the immunocompetent population, these so-called cutaneous HPV types induce asymptomatic chronic infections. However, members of the genera gamma, mu, and nu can induce benign skin lesions, e.g., cutaneous papillomas or warts (5). In addition, there is growing evidence of an etiological role of beta HPVs in non-melanoma skin cancer (NMSC), most likely in association with ultraviolet (UV) radiation [reviewed in (19–21)].
In this review, we describe the natural history of mucosal and cutaneous HPV infections, and discuss the transforming properties of a subset of them that have been shown or are suspected to play a causative role in various human cancers. In particular, we focus on describing the state of knowledge about the role of cutaneous HPV types in NMSC, because many existing reviews already cover all the aspects of mucosal HPV type infection and its role in transformation.
Genomic Organization of HPV and Viral Gene Products
The HPV genome is organized into three major regions: (i) a long control region (LCR), also called the upstream regulatory region (URR), located between the L1 and E6 open reading frames (ORFs), which contains the early promoter and regulatory element involved in viral DNA replication and transcription; (ii) the early region, which encodes the E1, E2, E4, E5, E6, and E7 proteins involved in viral gene expression, replication, and survival; and (iii) the late region, which encodes the structural proteins (L1 and L2). This organization is shared among all alpha HPVs (Figure 1). However, only four ORFs (those of E1, E2, L1, and L2) are necessary to fulfill the requirements to ensure the viral replication and shedding of the virus, and are present in all known PVs (22). Certain HPV genera and types lack an ORF. For example, the E5 ORF is lacking from HPV types that belong to the genera beta, gamma, and mu, and both the E5 and E6 ORFs are lacking from three gamma HPV types (HPV101, 103, and 108) (23, 24). LR HPV types (i.e. HPV6 and 11) encode for two E5-like proteins E5γ and E5δ (25). Beta and gamma HPV types have the particularity of harboring a shorter LCR compared with members of other genera (Figure 1).
Figure 1. Genome organization of alpha and beta HPV types. The circular dsDNA genome of high-risk alpha and beta HPV types is represented by a black circle. The positions of the early and late genes are represented by blue and orange, respectively. Both represented PV genera have the potential to express an E8∧E2C transcript that encodes for a protein that includes an E8 domain fused to the hinge and DNA-binding domains of E2. LCR, long control region.
All PVs also have the potential to express an E8∧E2C transcript (or equivalent), encoding for a protein that includes the E8 domain fused to the hinge and DNA-binding domains of E2 (26). This protein acts as a transcriptional repressor, and represses E1/E2-dependent replication of the viral origin (27–29).
Natural History of Mucosal HPV Infections
Mucosal HPV infections occur during the first sexual exposures in early adulthood (30), although non-sexual infection may also be possible (31). Although a majority of sexually active women will acquire a genital HPV infection, most (>90%) cervical HPV infections are resolved by the host immune system within 1–2 years (32, 33) and give rise only to asymptomatic infections. However, a minority of HPV infections become persistent. The risk of developing epithelial cell abnormalities and cancer is then increased (32). Persistent infection may be explained by several factors. Alcohol consumption and the HR HPV load synergistically increased the risk of persistent HR HPV infection in women (34). Cigarette smoking alone also plays an important role in the acquisition of persistent HPV infections in women, by decreasing the probability of clearing oncogenic infections (35). Smoking is also a risk factor for a persistent oral HPV infection (36). Host genetic risk factors may also predispose an individual to persistent HPV infection and developing cervical cancer, as supported by the high rate of heritability of cervical cancer (37–39). Defective immune response due to genetic variations (e.g., inflammasome genetics) has been associated with virus persistence and progression to cervical cancer (39–41). Human leukocyte antigen (HLA) genes also play a key role in the persistence of HPV infections and progression to cervical cancer, depending on their ability to bind to HPV antigens (39, 42).
Many natural HPV intra-type variants that differ from the prototypes have been identified, and some of them have been associated with the persistence and clinical outcome of cervical HPV infections (43–46). Several studies showed that EUR-350 T and EUR-350 G, which are HPV16 E6 variants, could influence viral persistence (47, 48). Recently, Mirabello et al. (49) evaluated a large collection of 5,570 HPV16-infected case–control samples using HPV whole-genome sequencing, and showed that E7 variation from the prototype sequence greatly decreases the risk of invasive cancer. That study also highlighted the presence of rarely detected variants of HPV16 that are consistent with the antiviral activity of the human APOBEC3 cytidine deaminase family. Moreover, APOBEC3A is stabilized in human keratinocytes by HR HPV E7, but not by E7 of LR HPV types (50).
The HR mucosal HPV genome integrates into host chromosomes in anogenital and oropharyngeal cancers, and it is considered an important driver of carcinogenesis.
There is not yet consensus on whether this is an early or a late event in carcinogenesis (51, 52). The integration occurs frequently within the E1 or E2 proteins, leading to dysregulation of expression of the oncoproteins E6 and E7 due to loss of the E2 repression function (53). The stability of E6 and E7 mRNAs is increased when HPV16 DNA is integrated into the human genome (54).
The frequency of viral integration increases with the severity of the cervical precancerous lesions. Integration of HPV into the host genome is reported in a majority of cervical cancers (83%) (55), and HPV18 has been found integrated into almost all cervical cancers (56, 57). HPV16 can be present as episomal or integrated forms, or both (58).
The DNA integration event could contribute to the HPV-driven malignant transformation of cervical cells, by affecting the expression of key host genes near the integration site (59). In addition, the existence of long-distance regulation among the integrated HPV fragment, the MYC gene, and the 8q24.22 region has been reported (60).
Although it was commonly accepted that integration occurred randomly into the host genome, several studies have shown that non-random HPV integration occurs. Integration of HPV16 sequences appears to take place preferentially in common fragile sites in human chromosomes (61, 62), in transcriptionally active regions (63), or near microRNAs (64).
Carcinogenesis can also be enhanced by epigenetic events. In cervical cancer, HPV16 E2 binding sites (E2 BS1 and E2 BS2), which play a role in transcriptional repression of E6 and E7 oncoproteins, are heavily methylated. The binding of E2 protein on these sites is then prevented, which leads to the upregulation of E6 and E7 oncoproteins (65). Host or viral epigenetic changes may be enhanced by the integration of HPV. A recent study showed that L1 gene methylation increases significantly according to the grade of the cervical lesion, and a high methylation rate of this gene correlates with the physical status of HPV integration (66, 67).
Data on HPV integration in HNC are still sparse and contradictory. Parfenov et al. (68) reported that 25 out of 35 HNCs showed integration of HPV into the human genome. Other studies reported HPV integration rates of 39% in oral squamous cell carcinoma (SCC) (69) and 43% in HNC (70), which is lower compared with HPV-positive cervical cancer. Episomal, integrated, and mixed forms of HPV are found in HPV-positive HNC. Moreover, a recent study, based on analysis of the Cancer Genome Atlas sequencing data, reported for the first time the presence of viral–human hybrid episomes in HPV-positive HNC (71).
Natural History of Cutaneous HPV Infections and Clinical Implications
Cutaneous HPV types are ubiquitous and are widespread in the general population. Up to 90% of healthy individuals tested positive for beta HPV types (72–74). The infection occurs in young children through skin-to-skin contact (75, 76). The cutaneous HPVs have been proposed to infect hair follicle stem cells of healthy individuals, where they constitute a reservoir of persistent infection (77, 78). DNA from cutaneous HPVs is frequently detected in hair bulbs, independently of the anatomical region where the hairs were plucked (79). Schmitt et al. (80) showed that keratinocyte stem cells (KSCs) in rabbit hair follicles co-localize with the primary target cells of cottontail rabbit papillomavirus (CRPV). They showed presence of CRPV early transcripts in clonogenic cells of the hair follicles soon after infection, suggesting that hair follicle stem cells are the initial target cells for CRPV (80), which was shown to induce cutaneous warts and carcinomas (81).
Epidemiological studies in humans showed that the cutaneous HPV population present in hair follicles mirrors the HPV prevalence in the skin of the same individual, which makes eyebrow hair an excellent sampling method to characterize the individual cutaneous HPV population (79, 82). Studies on intra-familial transmission showed that similar spectra of beta HPV types are present within members of the same family. Babies and their parents share some of the beta HPV types (83), which can persist for many years on healthy skin (84, 85). Transmission of beta and gamma HPVs has been demonstrated in couples (86).
Exposure to cutaneous HPV is common. Serological studies measuring antibodies against the type-specific L1 major capsid protein of HPV showed that 52% of the Dutch population and 67% of the Italian population were exposed to beta HPV infections (74). A recent seroprevalence study, based on ten beta HPV types, showed that 39% of healthy men were seropositive for at least one beta HPV type (87). Seroconversion for beta and gamma HPV types appears to be slow and to increase with age (88). Low viral loads in immunocompetent individuals (89) and continuous renewal of the infected keratinocytes may explain why only half of the infected individuals develop antibodies against cutaneous HPV types (90, 91).
Cutaneous HPVs are highly prevalent in the rare hereditary disease epidermodysplasia verruciformis (EV), which tends to progress to cutaneous SCC (cSCC), frequently located at sun-exposed anatomical sites (92). HPV5 and 8, two members of species β1, were isolated from cSCC in patients with EV (93). These two types were classified as possibly carcinogenic to humans (IARC Group 2B) (10).
In organ transplant recipients, a 65–250-fold increased risk of developing SCC compared with the general population was reported (94, 95). In addition, in HIV-positive individuals, several studies reported at least a 2-fold increased risk of cSCC compared with HIV-uninfected people (96, 97). The correlation between the immunodeficiency state and an increased risk of developing NMSC suggested a possible role of infectious agents (98), such as cutaneous HPV.
This observation is corroborated by growing evidence showing an association between cutaneous infection and the risk of developing NMSC in immunocompetent individuals under certain conditions (e.g., UV radiation exposure) [reviewed in (21, 99)], in particular in cSCC (100–102).
Epidemiological and biological data suggested that beta HPV types, and species β1 and β2 in particular, may be linked to cSCC development in immunocompetent individuals (99, 100) [reviewed in (21)].
Interestingly, the cutaneous HPV viral load was higher in actinic keratosis, which is considered to be the precursor lesion of cSCC, compared with cSCC (< 1 copy per cell), suggesting a possible role of cutaneous HPV types in the initiation of skin carcinogenesis but not in the maintenance of the cancer phenotype, by exacerbating the accumulation of UV radiation-induced DNA breaks and somatic mutations (the hit-and-run mechanism) (103, 104). The absence of HPV mRNA (105) and the lack of evidence for the integration event in skin tumors support this scenario.
The Mastomys natalensis papillomavirus (MnPV)-infected rodent model Mastomys coucha was previously described and used to evaluate the role of PVs in NMSC. In these mice, which can be naturally infected by MnPV, skin lesions such as papillomas or keratoacanthomas can be induced (106). Using this model, Hasche et al. (107) clearly demonstrated cooperation between UVB radiation exposure and MnPV infection in the first step of NMSC initiation. MnPV-positive animals chronically exposed to UVB radiation developed lesions significantly more frequently than MnPV-negative animals (107). Moreover, the authors showed higher MnPV viral loads in well-differentiated and keratinizing SCCs (KSCCs) compared with normal skin. Of note, these KSCC lesions share histological similarities to human SCC pre-malignant lesions, such as actinic keratosis, that normally show higher cutaneous HPV DNA loads than SCC, as discussed above. In addition, non-keratinizing SCCs (nKSCCs) occurring in mice UV-irradiated at higher doses showed low levels of (or no) viral DNA and often harbored p53 mutations, again mirroring the evolution of SCC in humans. The loss of viral DNA and viral gene expression in nKSCC lesions is most likely due to their undifferentiated state, which no longer sustains the infection. However, despite the lack of viral DNA, animals with nKSCCs developed antibody response to MnPV capsids, which clearly highlights past exposure to the virus, similar to the situation observed in SCC cases (108). Therefore, the Mastomys model shows parallels with a natural infection by cutaneous HPV and provides a good model to study the association between cutaneous HPV infection, UV radiation and SCC. Most importantly, this model provides further evidence for a hit-and-run mechanism.
Many independent studies reported the presence of beta HPV types at different anatomical regions other than the skin, such as the oral mucosal epithelium, genital sites, and the anal canal (109–112), or investigated the role of cutaneous HPV in various malignancies (e.g., male external lesions, breast cancers, salivary gland tumors, esophageal cancer) other than skin cancer, without showing any association (113–116). However, a recent study reported an increased risk of HNC in individuals who tested positive for HPV5 of species β1, as well as other HPV types of species γ11 and γ12 (117). The same trend was reported in another study, in which β1 HPV5 and β2 HPV122 were significantly associated with HNC (118).
Whereas, beta HPVs are known to infect cutaneous tissues, epidemiological data showed that β3 HPV types are also present in the mucosal epithelium, suggesting a dual tropism (111, 119). Moreover, this species shares biological similarities with mucosal HR HPV types such as HPV16 in in vitro and in vivo experimental models (120, 121). A study showed that β3 HPV49 transgenic mice, after 4-nitroquinoline 1-oxide treatment, were prone to develop upper digestive tract tumors (121).
Epidemiological and biological data on the role of gamma HPV types are sparse, and do not support an etiological role in NMSC. However, this genus includes an increasing number of members that may deserve more investigation. Using deep sequencing, a member of species γ27 (HPV197) has been isolated exclusively from skin cancers. However, additional studies are required to demonstrate an etiological link between this type and skin cancer (122, 123).
Entry and Life Cycle of HPV
The viral particle reaches the basal layer of the epithelium via micro-wounds or micro-fissures, or via hair follicles (80, 124, 125), to infect basal keratinocytes or stem epithelial cells. In the cervix, HR HPV may directly reach the single layer of cuboidal epithelial cells from the squamocolumnar junction located between the endocervix and the ectocervix (126); this anatomical region is therefore more sensitive to infection. In the oral cavity, HPV infection is believed to occur in tonsillar crypts, where a comparable cellular structure exists (127).
The infectious entry of the virus into the cell is intimately linked to proliferating cells of the regenerating basal epithelium. Broniarczyk et al. (128) showed that the viral particle can remain infectious after several weeks on the surface of senescent cells, which are resistant to HPV infection; however, reactivation of the cell cycle by p53 siRNA led to the entry of the virus (128).
To initiate the infection, the virus binds to the heparan sulfate chains of proteoglycans (HSPGs) located on the cell membrane or on the extracellular matrix (ECM) through the major capsid protein L1, leading to cyclophilin (CyP) B-mediated conformational changes of the capsid structure, which expose the minor capsid protein L2 on the surface of the viral particle (Figure 2). The cleavage of L1 in the extracellular space by a serine protease, kallikrein-8, appears to be crucial for the efficient externalization of L2 (129). On the ECM, the virus can also bind to a transient binding receptor, laminin-322 (formerly called laminin-5) (130).
Figure 2. Entry of HPV. The virus binds to the heparan sulfate chains of proteoglycans (HSPGs) located on the cell membrane or on the extracellular matrix through the major capsid protein L1, which leads to the conformational change of the capsid and the subsequent L2 externalization on the surface of the viral particle. Next, the N-terminus of the minor capsid protein L2 is cleaved at a conserved furin cleavage consensus site, which facilitates the binding of HPV to an unidentified secondary receptor, allowing subsequent internalization that is actin-dependent. The L2–viral DNA complex and a portion of L1 are delivered to the trans-Golgi network, through retromer- and retriever-mediated vesicular trafficking. Inside the late endosomes, the capsid of the viral particle is disassembled in a low-pH environment, leading to the dissociation of the major capsid protein L1 from the minor capsid L2, mediated by host-cell cyclophilins. The L2–viral DNA complex is then delivered into the nucleus to the nuclear domain (ND10), during cell mitosis.
Next, the N-terminus of the minor capsid protein L2 is cleaved with a CyP-independent mechanism at a conserved furin cleavage consensus site to expose the L2 amino acids 17–36 (the RG1 neutralizing epitope), which appears to be important for the interaction of the capsid with an unidentified secondary receptor (131, 132). The binding of HPV to the secondary receptor allows subsequent internalization.
Using HPV16 pseudovirion (PsV) particles covalently labeled with some fluorophores Schelhaas et al. (133) showed that HPV16 PsVs moved along the outside cell membrane in correspondence to the actin protrusions and were transported inside the cells by actin retrograde flow. Upon inhibition of the actin flow by blebbistatin, the virus particles could not be actively transported along actin-rich protrusions, which resulted in reduced efficiency but not in complete abolishment of the infection (133). In a more recent work, using different compound inhibitors and siRNAs, Schelhaas et al. (134) showed that the entry of HPV16 into HeLA and HaCaT cells was actin-dependent and clathrin-, AP2-, caveolin-, flotillin-, and dynamin-2-independent (134). Although the exact mechanism of membrane protrusion/vesicle formation during virus endocytosis remains unclear, Schelhaas's study clearly showed that actin polymerization and depolymerization are crucial for HPV16 endocytosis, in particular for scission of endocytic vesicles (134). Once inside the cell, the virus needs to reach the replication machinery within the nucleus (Figure 2). The virus is delivered to an early endosome or macropinosome-like endosome that matures into a late endosome [reviewed in (135)]. Inside the late endosomes, the capsid of the viral particle is disassembled in a low-pH environment, leading to the dissociation of the major capsid protein L1 from the minor capsid L2, mediated by host-cell CyPs (136). However, a recent study showed that residual conformationally intact L1 protein remains in complex with the viral genome (137). The L2–viral DNA complex traffics to the trans-Golgi network, mediated by the retromer complex (138). Retromer is a heteropentameric complex composed of dimer of sorting nexins (SNX1, SNX2, SNX5, and SNX6) and a vacuolar protein sorting (Vps) trimer containing Vps26, Vps29, and Vps35, which play a central role in the retrieval of several different cargo proteins from the endosome to the trans-Golgi network. The trafficking also involves retriever, a multi-protein complex analogous to retromer (139), and involves the cellular adaptor protein SNX17. This complex is associated with endosomes and is essential for the recycling of multiple cargoes. A proteomic approach has shown that SNX17 interacts with HPV16 L2 (140). The observation that the SNX17 binding site is conserved among different PV genera suggests that it has an important role in regulating the PV life cycle and replication. The SNX17/L2 interaction is indeed required for HPV16 infection and appears to play key roles at different stages of the infection, being important for effective capsid disassembly and L1 dissociation. The interaction of L2 with SNX17 also contributes to the maintenance of viral capsids in the late endosomal compartment, by protecting them from lysosomal degradation and allowing L2–DNA complexes to egress from the endosomes (141). A related protein, SNX27, is part of the retromer complex and also interacts with L2, but through its PDZ domain. SNX27 appears to be important for virion trafficking and to cooperate with SNX17 in this event (142).
L2 is required for efficient trafficking of the viral genome to the nucleus (138) and remains associated with the viral DNA when the transport vesicles deliver the L2–viral DNA complex inside the nucleus [reviewed in (143)], most likely together with residual L1.
Delivery of the L2–viral DNA complex into the nucleus requires cell cycle progression into early mitosis (144). Nuclear envelope disruption is also required for nuclear import of the L2–viral DNA complex (145) (Figure 2).
Inside the nucleus, L2 mediates the viral genome delivery to the nuclear domain (ND10) (146, 147), also known as promyelocytic leukemia (PML) bodies, where early viral transcription and replication will take place, as for many DNA viruses. ND10 structures are host restriction factors that limit viral infection as part of the intrinsic defense against viral infection (148). Stepp et al. (149) showed that Sp100 nuclear antigen, one of the components of ND10, represses transcription, replication, and establishment of incoming HPV DNA in the early stages of infection. Upon initial infection, the minor capsid protein L2 leads to the alteration of the composition of ND10 protein, leading to the release or degradation of Sp100 (150). The Daxx protein, another component of ND10, is also recruited (150). It was previously reported that this protein modulates the early gene expression and the transient replication of HPV genomes in U2OS cells (151).
The replication cycle of HPV is linked to the differentiation of the infected epithelium, and starts with a first step, called “establishment replication,” which consists of maintaining a constant number of episomal copies (50–100 per cell) (Figure 3) [reviewed in (152)]. Viral DNA replication relies on the host DNA replication machinery and is supported by the early viral proteins E1 and E2 (153). Upon nuclear entry, viral DNA replication is initiated by the binding of E2 on specific sites located on the LCR, which is required for the recruitment and binding of E1 helicase to the viral origin of replication (154). The cellular proteins TopBP1 and Brd4 are involved in the initiation of the HPV16 E1/E2-mediated DNA replication (155). Limited viral genome amplification may be obtained by the interaction of NCOR/SMRT repressor complexes with E8∧E2C proteins, which inhibits viral replication (156, 157). After this initial step, which aimed to generate a low copy number of genomes, the maintenance phase is initiated (Figure 3) [reviewed in (152)]. This phase consists of creating the conditions to maintain a constant number of viral genomes in the nuclei of undifferentiated basal cells as an extrachromosomal genome to create a persistent infection. The viral genomes need to be correctly segregated during cell division through attachment/tethering of the virus genome to the host-cell chromosome. The transactivation domain at the N-terminal part of E2, with the bromodomain protein Brd4, interacts with the host mitotic chromosomes. Concomitantly, the E2 DNA-binding/dimerization domain binds to E2 binding sites in the LCR of the viral genome (158–162). Brd4 interacts with the E2 transactivation domain of most PVs (163). ChIR1, an ATP-dependent DNA helicase, seems to play a key role by regulating the chromatin association of HPV16 E2, and in maintaining the episomal form of HPV16 (164–166). The structural maintenance of chromosomes (SMC) proteins SMC5 and SMC6 may also play a role in viral genome maintenance by interacting with E2 (167). In addition to E1 and E2, the oncoproteins E6 and/or E7 are required for stable episomal maintenance of HPV (168, 169).
Figure 3. Life cycle of HPV. The first step of the replication cycle of HPV, called “establishment replication,” consists of maintaining a constant number of episomal copies. Viral DNA replication relies on the host DNA replication machinery and is supported by the early viral proteins E1 and E2. After this initial step, the maintenance phase is initiated. This phase consists of creating the conditions to maintain a constant number of viral genomes in the nuclei of undifferentiated basal cells as an extrachromosomal genome to create a persistent infection. In addition to E1 and E2, the oncoproteins E6 and/or E7 are required for stable episomal maintenance of HPV. Upon cell differentiation in stratified epithelium, vegetative or productive viral replication is initiated, with the subsequent production of progeny virions. Here, the oncoproteins E6 and E7 expressed at relatively low levels in differentiated cells play a key role by inactivating tumor suppressor proteins. The activation in differentiated epithelial cells of the late promoter, located in the E7 region, leads to the production of the high levels of E1 and E2 viral proteins required to ensure viral DNA replication. E4 and E5 contribute to efficient productive replication. The capsid proteins L1 and L2 are expressed from the late promoter, and are involved in the encapsidation of newly replicated genomes, resulting in virion release in the superficial layers during desquamation. E4 also plays a role at this step of the viral life cycle, by interacting with the keratin network. The right side of the schematic representation illustrates that a minority of HPV infections become persistent, thus increasing the risk of keratinocyte transformation, through integration of viral DNA into the human genome, which leads to the loss of E2 repression function and subsequent E6 and E7 upregulation.
Upon cell differentiation in stratified epithelium, vegetative or productive viral replication is initiated, with the subsequent production of progeny virions (170, 171), [reviewed in (152)] (Figure 3). Uninfected cells leave the basal layer for terminal differentiation, are withdrawn from the cell cycle, and stop the replication of the DNA. Because the amplification of the viral genome requires cellular conditions that enable cell proliferation and thus DNA replication, the virus has developed strategies to prevent cell-cycle arrest and apoptosis signals. Here, the oncoproteins E6 and E7 expressed at relatively low levels in differentiated cells play a key role by inactivating tumor suppressor proteins (e.g., p53, retinoblastoma protein [pRb]) and activating signal transduction, to ensure that the infected cells remain active and progress to the S phase. The activation in differentiated epithelial cells of the late promoter (P670 for HPV16, P811 for HPV18, and P742 for HPV31), located in the E7 region, leads to the production of the high levels of E1 and E2 viral proteins required to ensure viral DNA replication. E4 and E5 contribute to efficient productive replication (172, 173). The capsid proteins L1 and L2 are expressed from the late promoter, and are involved in the encapsidation of newly replicated genomes, resulting in virion release in the superficial layers during desquamation (174) (Figure 3). E4 also plays a role at this step of the viral life cycle, by interacting with the keratin network (175). A more complete description of the function of each of the viral proteins is provided in Table 1.
Table 1. Main features and functions of the early and late gene products from alpha and beta HPV types.
Epigenetic regulation of HPV transcription plays an important role in the virus life cycle. HPV regulatory epigenetic mechanisms have recently been partly elucidated by Pentland et al. (176). The key players in this regulation are the chromatin-organizing CCCTC-binding factor (CTCF) and Yin Yang 1 (YY1), binding, respectively, to the E2 ORF and the viral LCR. CTCF and YY1 are involved in the formation of a loop allowing E2 and LCR interaction and leading to the formation of epigenetically repressed chromatin in the HPV18 genome. This event leads to attenuated expression of oncoproteins in undifferentiated cells. CTCF plays a key role in this event; in fact, the HPV18 ΔCTCF genome, carrying a mutation in the CTCF binding site, shows higher H3K4me3 and reduced H3K27me3 marks in the viral enhancer and in the downstream early promoter. This is coherent with reduced recruitment of polycomb repressor complex 2 (PCR2), enrichment of RNA pol II binding, and a consequent increase in the transcriptional activity of the viral early promoter. Pentland et al. (176) also showed that both CTCF and YY1 are necessary to repress viral transcription. During cell differentiation, a drop in YY1 expression levels leads to loss of chromatin loop formation and to increased expression of the oncoproteins consecutive to epigenetic depression of the viral genome (176).
Transforming Activities of PVs
Studies performed on mucosal HR HPV types, such as HPV16 and 18, showed that both oncoproteins E6 and E7 play a key role in cervical cancer by altering pathways involved in the host immune response to establish a persistent infection and promote cellular transformation. The E6 (150 amino acids) and E7 (100 amino acids) oncoproteins play a central role in carcinogenesis by interacting with a large number of cellular proteins involved in key cellular events, such as the cell cycle and apoptosis control. In addition to E6 and E7, HR HPV types encode a small hydrophobic oncoprotein of < 90 amino acids, E5, which also appears to play a role in HPV-induced carcinogenesis. Continuous expression of the E6 and E7 oncoproteins is essential to initiate and to maintain the cellular transformation process. Indeed, the use of different strategies to inhibit the function of E6 and E7 in HPV-positive cancer cells resulted in cell growth arrest by apoptosis or senescence (177, 178). The main targets of the E6 and E7 oncoproteins are p53 and pRb, respectively; p53 regulates DNA damage response and apoptosis, and pRb tightly controls the cell cycle. Although E6 and E7 of mucosal HR HPV types were extensively studied, functional data on E6 and E7 of beta HPV types are sparse. Moreover, the absence of cutaneous E6 and E7 transcripts in skin tumors suggests an alternative mechanism for promoting cancer development (hit-and-run mechanism). The beta and gamma HPVs lack the E5 ORF.
Transforming Activities of E6
The E6 oncoprotein of HPV16 is a small protein of ~150 amino acids, made up of two zinc-like fingers joined by an interdomain linker (179). E6 of HR HPV types has the ability to bind to the LXXLL peptide motif of cellular proteins. Among them, the E6-associated protein (E6AP) is an E3 ubiquitin ligase that targets proteins for ubiquitination and degradation by the proteasome (180). The best-characterized E6/E6AP interaction is with the p53 tumor suppressor (181). In normal cells, p53 plays a key role to safeguard the integrity of the cellular genome, avoiding the proliferation of cells with damaged DNA. Therefore, upon DNA damage insults, p53 will either induce cell-cycle arrest at the G1 phase to allow DNA repair or, if the DNA damage is too extensive, activate apoptotic pathways. E6/E6AP interaction leads to conformational changes of E6 that allow the association with the p53 pro-apoptotic tumor suppressor (182), leading to its ubiquitin-mediated proteasome degradation and, as a consequence, to accumulation of DNA damage and to genomic instability (183).
Although E6 of the β3 HPV49 type also showed the ability to degrade p53 by an E6AP-dependent mechanism, E6s of HPV species β2 do not share this property and have evolved alternative mechanisms to counteract p53 functions [reviewed in (19)]. For example, E6 of HPV23 interferes with the ability of the homeodomain-interacting protein kinase 2 (HIPK2) to phosphorylate p53 on serine 46 to activate its apoptotic function upon UVB irradiation (184).
In addition to inducing its degradation, E6 of the mucosal HR HPV types can inactivate p53 and abolish its transactivation in vivo (185), by targeting its transcriptional coactivator CBP/p300. This event leads to the displacement of p53 from CBP (186) and to the inhibition of p300-dependent p53 acetylation (187).
As for HR HPV types, the E6 oncoproteins of β1 HPV types (HPV5, 8) also have the ability to interact with p300 with high efficiency (188). HPV5 and 8 E6s inhibit the association of AKT with the p300 C-terminus that is needed to ensure p300 stability, thus leading to its proteasomal degradation (188). The degradation of p300 results in decreased protein levels of ATR, a PI3 kinase family member, which plays a key role in UV radiation damage signaling. This event in turn reduces ATR's ability to protect cells against UVB radiation-induced damage (189). The reduced ATR level results in decreased phosphorylation and subsequent attenuated accumulation of p53 in cells in response to UVB radiation exposure. In addition, in HPV5, 8, and 38 E6-expressing cells, arrest in the G1 phase of the cell cycle is prevented, which hampers DNA repair. Both events lead to increased persistence of thymine dimers and UVB radiation-induced double-strand breaks in these cells (189). More recently, Hufbauer et al. (190) showed a lack of phosphorylation of ATM, ATR, and Chk1 in HPV8 E6 expressing-monolayer and organotypic cultures leading to the impairment of DNA damage sensing and repair. Consistent with impaired activation of the ATM/ATR pathway, immunohistochemical analysis revealed the presence of thymine dimers in UVB radiation-treated E6-expressing cells (190), further confirming a role of beta HPV E6 proteins in facilitating UV radiation-induced DNA damage.
E6 of HPV8 may also play a role in skin carcinogenesis by down-regulating CCAAT/enhancer-binding protein α (C/EBPα), a transcription factor with a role in cell differentiation. MicroRNA-203 (miR-203), previously shown to be an important regulator of epidermal proliferation and differentiation, is induced by C/EBPα. In HPV8 E6 expressing keratinocytes, downregulation of C/EBPα via E6-mediated p300 degradation led to loss of miR-203 expression and induction of the ΔNP63α protein, altering the proliferation abilities of the infected epithelial cells (191).
Similar to E6s of HPV species β1, E6 of HPV38 (species β2) was able to prevent p300-mediated acetylation of p53, an event that resulted in its inhibition and thus led to the immortalization of primary keratinocytes (192), indicating that the ability of E6 proteins of HR cutaneous HPVs to induce changes in p53 acetylation status is crucial for the transforming ability of these viruses (192).
E6 has additional cellular targets that play an important role in apoptosis (e.g., Bak, survivin, TNF R1, FADD, procaspase 8) (193–197).
The anti-apoptotic protein Bak, a member of the Bcl-2 family of proteins, plays an important role in regulating the apoptotic process by forming pores with Bax to permeabilize the outer mitochondrial membrane. Bak has been shown to be targeted for degradation by different HPVs from the alpha and beta genera, an event that contributes to preventing UV radiation-induced apoptosis. The ability of the E6 proteins of beta HPVs (e.g., HPV5, 8, or 38) to target Bak for degradation, is important for promoting the survival of DNA-damaged cells and therefore for the progression of NMSC. Bak interacts with the ubiquitin ligase E6AP; however, E6 proteins of HPV5 or 8 that do not interact with E6AP (198) can still induce Bak degradation. Moreover, HPV5 E6 mediates Bak proteolysis when expressed in E6AP-null mouse embryo fibroblast (MEFs) (199), indicating that HPV5 degrades Bak via an E6AP-independent mechanism. Holloway et al. (200) have shown that HERC1 ubiquitin ligase is required for HPV5 E6-mediated Bak degradation (200). In response to UV irradiation, Bak is dephosphorylated on K113 and interacts with HERC1. The association of HERC1 and Bak is dependent upon E6 expression. In fact, HERC1 is recruited by HPV5 E6 to degrade Bak upon UV radiation-induced damage. The degradation of Bak by E6 prevents the release of proapoptotic factors including apoptosis-inducing factor (AIF) from the mitochondria into the nucleus, an event that prevents cells from undergoing UV radiation-induced apoptosis (201).
Those experiments highlight the ability of beta HPVs to enhance UV radiation-induced DNA damage, supporting a role of beta HPVs at the early stage of skin carcinogenesis.
Beta HPVs display additional features for increasing the oncogenic potential of UV radiation exposure by attenuating DNA damage repair. For instance, the binding and destabilizing of p300 by beta HPV5 and 8 E6s also lead to lower levels of BRCA1 and BRCA2, critical to the homology-dependent repair of double-strand breaks (202). Viarisio et al. (104, 121) reported a strong reduction of IL-18 production in HPV38 transgenic mice exposed to UV radiation (104, 121). IL-18 reduces UV radiation-induced DNA damage by the induction of DNA repair (203), suggesting that by inducing a reduction of IL-18 secretion in the mouse skin, HPV38 early proteins may repress the UV radiation-induced inflammasome response, favoring DNA damage accumulation and therefore promoting the carcinogenic process (104).
All HR HPV E6 proteins contain C-terminal PDZ domain-binding motifs (PBMs) that enable interaction with several cellular PDZ domain-containing proteins involved in the regulation of different processes, such as cell polarity, migration, and attachment (204). The interaction of HR HPV E6 with the PDZ proteins leads to their proteasome-mediated degradation (205). The PBM is absent from E6s of LR HPV types, suggesting a key role for this domain in HR HPV-induced carcinogenesis.
A phosphor-acceptor site for protein kinase A (PKA) is located within the PBM. Phosphorylation at this site negatively regulates the PDZ domain-binding activity. The PKA recognition sequence is highly conserved among E6 of HR HPV types and is not present in LR HPV E6 oncoproteins, further suggesting that the PDZ domain-binding activity is crucial in HR HPV-mediated transformation (206, 207). In addition, phosphorylation at this site is required for HPV18 E6 to interact with 14-3-3ξ, contributing to the stabilization of the oncoprotein (208).
Different PDZ domain-containing proteins have been identified as targets of E6 of HR HPV types: DLG1, a human homolog of Drosophila discs large 1, and hScrib, a homolog of the Drosophila scribble protein (209–211), which play an important role in the polarity of the epithelial cells and have been identified as tumor suppressors (212); the molecular scaffolds membrane-associated guanylate kinase homology proteins, such as MAGI-1 (213), which is involved in the modulation of epithelial cell adhesion and tight junction integrity, acts as a tumor suppressor through the stabilization of PTEN (214, 215); and another HPV16 PDZ domain-containing protein, CAL, which is involved in the regulation of intracellular vesicular trafficking (216). The multi-PDZ domain protein MUPP1, which negatively regulates cellular proliferation, is targeted for degradation by HPV18 E6 and possibly by other HR HPV E6 oncoproteins, thus contributing to the transforming activities of E6 oncoproteins (217). The PDZ domain-containing proteins PSD95, TIP-2/GIPC, and NHERF1, which are involved in different signaling pathways, are also targeted by E6 (218–220).
In addition to playing a role in many cellular processes, the PBM of E6 proteins also regulates the HPV life cycle; mutations of this motif were shown to lead to lower levels of episomal HPV genome and cell growth (221). HPV8 E6, which lacks the PBM, can target the PDZ domain-containing protein syntenin 2, by reducing its expression level in a MAML1-independent manner, by an epigenetic mechanism. This event plays a key role in hampering differentiation of HPV8-expressing keratinocytes (222) and provides a further example of how skin and mucosal HPV types evolved different mechanisms to alter normal cell homeostasis and induce cellular transformation.
Using a proteomic approach, Thomas et al. (223) compared the ability of E6 PBMs of HPV types with different oncogenic potentials to interact with PDZ substrates (223). This comprehensive analysis showed that the number of PDZ substrates varies according to the HPV type. Most of the PDZ substrates were recognized by HR types (i.e., HPV16 and 18), whereas HPV66, which is a possible HR HPV type, recognizes only a limited number of substrates. Of note, whereas PDZ substrates such as DLG1 are common targets of the E6 proteins irrespective of their oncogenic potential, the ability of E6 PBMs to interact with hScrib, a component of the Scribble polarity complex, is correlated with the transforming ability of the studied HPV types (223).
Telomerase activity also plays an important role during cervical carcinogenesis, because it increases with the grade of the cervical lesions (224). The ability to activate the promoter of hTERT (human telomerase reverse transcriptase), the catalytic subunit of the telomerase, by different HPVs was assessed using a luciferase-based assay. The results of that study showed that several HR HPV types significantly upregulated hTERT promoter activity (225). The induction of hTERT requires the E6/E6AP complex to interact with the heterodimer c-Myc/Max (226–229). Moreover, E6 recruits histone acetyltransferase (HAT) at the hTERT promoter, which facilitates the access of transcriptional activators to key regulatory sites through the opening of the chromatin (230).
NFX1-91, a constitutive hTERT repressor, is targeted by E6/E6AP for ubiquitin-dependent degradation (230). The induction of hTERT expression by E6/E6AP is also dependent on the activation of NFX1-123, which constitutes a level of post-transcriptional regulation (231). Together, these events lead to extending the life span of primary human keratinocytes, and to their immortalization. Proteomics studies comparing the cellular protein-binding partners of the E6s of HPV types of different genera have shown that whereas all the studied alpha HPV E6s have a high affinity for E6AP and HERC2 (232), beta and gamma HPV E6s specifically interact with the cellular transcriptional co-activator MAML1 (233). The latter interaction is required to inhibit Notch1, a key player in skin differentiation and in the activation of transcription of the cell-cycle inhibitor p21. Interaction of beta HPV E6 with the Notch pathway appears to be important to promote the transformation process of the infected keratinocytes (234, 235). The different affinities of E6s of different genera for specific cellular targets could be a consequence of their tropism and co-evolution with different type of tissues (236).
Transforming Activities of E7
The E7 oncoprotein is made up of ~100 amino acids and has a C-terminal zinc-binding domain. HPV E7 presents three conserved regions, homologous to adenovirus E1A: CR1, CR2 in the N-terminal part, and CR3 in the C-terminal part. Moreover, HPV E7 comprises two nuclear localization sequences (NLSs) and one nuclear export sequence (NES), which enable E7 to be located in both the nuclear and cytoplasmic compartments, where E7 has different functions (237).
The CR1 domain has an UBR4/p600 binding site; the 600-kDa pRb-associated factor p600 is required for membrane morphogenesis, anchorage-independent growth, and cellular transformation (238, 239). Moreover, the interaction between BPV-1 E7 and p600 inhibits apoptosis, contributing to viral-induced transformation (240). UBR4/p600 ubiquitin ligase also appears to be involved in E7-mediated proteasomal degradation of the potential tumor suppressor PTPN14 (241), identified as a target of HR HPV E7 in cervical tumors (241, 242). The ability of HPV16 E7 to bind to p600 correlates with its capacity to transform cells (238). Similarly to E6, E7 binds to CBP/p300 and to its associated factor (pCAF) by interacting with the CBP TAZ2 domain (also known as CH3). The E7 binding site responsible for its interaction with the CBP TAZ2 domain overlaps with the LXCXE motif, which is crucial for E7 to bind to pRb. Formation of the ternary complex of E7/pRb and pCAF leads to the acetylation of pRb, an event that appears to play an important role in HPV-mediated transformation, probably by stabilizing the E7/pRb interaction (243, 244).
The CR2 domain comprises serine residues at positions 31 and 32 that are susceptible to phosphorylation by casein kinase II (CKII). Phosphorylation at these sites contributes to modulating some of the E7 functions. For example, CKII phosphorylation of E7 increases the binding affinity of HPV16 E7 for the TATA box-binding protein (TBP) (245), and it is required for efficient transformation by E7 (246). Interestingly, an additional phosphorylation site at serine position 29 exists in a natural E7 variant (N29S) and leads to increased levels of phosphorylation by CKII, which increases the interaction of E7 with TBP and pRb, and its transforming activity in primary rodent cells (247).
CR2 also comprises a conserved LXCXE motif located at residues 22–26 in HPV16 E7. This motif is essential for the interaction with pRb and related pocket proteins (p107 and p130) (248). An intact LXCXE motif is needed to induce cellular DNA synthesis and mediates cell transformation of immortalized NIH3T3 cells. Residues in the CR3 domain of E7 also seem to be required for its interaction with pRb (249). pRb and related pocket proteins (p107 and p130) play a key role in the regulation of the cell cycle by binding to and inhibiting E2F transcription factors (E2Fs), which regulates the expression of genes required for DNA synthesis and cell division. HR HPV E7 oncoproteins are able to target pRb and other members of the retinoblastoma family for proteasome-mediated degradation. This event requires the recruitment of the cullin-2 ubiquitin ligase complex, which binds to the CR1 region (250). Upon pRb interaction with HPV E7, E2F is released, with subsequent transcriptional induction of G1–S cell-cycle genes, such as cyclin A and cyclin E. The latter products are required for the activation of the cyclin-dependent kinase (CDK) complexes, which force the entry of HPV-infected cells into the S phase (251). This event leads to unscheduled growth and accumulation of genomic instability. LR HPV E7s bind to but are not efficient in targeting pRb for degradation (252). The ability of E7 to interact with pRb is not sufficient for HPV-induced transformation, which, as mentioned above, requires the interaction of E7 with additional cellular proteins (253).
The functionally and structurally pRb-related proteins, p107 and p130, are also targeted by the HPV E7 via the LXCXE motif, and both have the ability to bind to and inactivate E2Fs. These proteins are important negative regulators of the cell cycle, and they are believed to also act as tumor suppressors. However, unlike pRb, their role as tumor suppressors is still controversial, because related genetic alterations are not frequent in human cancers. Both HR and LR HPV E7s target p130 for proteasomal degradation, most likely by different mechanisms, thus highlighting an important role of this pocket protein in the HPV life cycle (254).
Analysis of HPV E7-interacting proteins showed that HPV16 E7 interacts with E2F6, which has the particularity of lacking the binding domain for pRb and acts as a transcriptional repressor. The association of HPV16 E7 with E2F6 abrogates its repressive activity on target genes, which results in increased DNA replication of HPV16-infected cells (255).
Using an organotypic culture of skin, Akgül et al. (256) showed that beta HPV8 E7-transduced keratinocytes displayed an invasive phenotype, in addition to losing their polarity (256). HPV8 E7-transduced keratinocytes also showed enhanced production of extracellular proteinases (MMP-1, MMP-8, and MT1-MMP) that play a role in the degradation of components of the basement membrane and the ECM during cell invasion. In another study, the same group used a transgenic mouse model expressing the E7 oncoprotein of HPV8 in the epidermis to investigate the molecular basis for HPV-induced invasion of skin keratinocytes. HPV8 E7 expression in mouse skin led to a reduction of E-cadherin and a subsequent upregulation of N-cadherin, an event that causes epithelial–mesenchymal transition. E7-positive keratinocytes also showed increased fibronectin expression and secretion. In addition, HPV8 E7-expressing cell surface showed enhanced levels of the integrin α3 chain, which plays a crucial role at the pre-invasive stage (257). E7 of beta-HPVs may play an important role in promoting an invasive phenotype in keratinocytes and may provide additional evidence to support a role of cutaneous HPVs in skin carcinogenesis.
HPV8 also has the ability to maintain cells in a “stemness-like” phenotype. Using in vitro colony formation and tumor sphere assays, Hufbauer et al. (258) showed that beta HPVs are able to increase the number of stem cell-like cells, as determined by measuring the number of HPV8 E2-, E6-, and E7-expressing cells harboring high levels of epithelial stemness cell surface markers, such as CD44 and epithelial cell adhesion molecule (EpCAM). EpCAM induction was also observed in organotypic skin cultures of HPV5, 8, and 16 E7-expressing cells, as well as in skin lesions of patients with EV (258). E7 expression also led to a reduction of the epithelial differentiation marker calgranulin B. Together, these results highlight a possible role of beta HPVs in skin carcinogenesis by increasing the number of stem cell-like cells present during early carcinogenesis and delaying cell differentiation.
Lanfredini et al. (259) evaluated KSC populations located in different regions of the hair follicle (259). In that study, using HPV8 transgenic (HPV8tg) mice that express the entire HPV8 early region, the authors identified KSCs from the upper part of the junctional zone as the site from which the expansion to the interfollicular epidermis takes place during the initial step of skin carcinogenesis. This population of stem cells express leucine-rich repeats and immunoglobulin-like domains protein 1 (Lrig-1) and are characterized by nuclear staining for the stemness marker p63. Skin accumulation of the p63 isoform, ΔNp63, which is induced by the expression of HPV8 early proteins, promotes the proliferation and expansion of Lrig1-positive KSCs. In healthy skin, the ΔNP63α transcript is directly targeted by miR-203, an event that plays a major role in the regulation of epidermal proliferation and differentiation, leading to “stemness” loss (260, 261). Of note, the skin lesions of patients with EV show strong p63 immunostaining associated with very low levels of miR-203.
In addition to features cited above, E7 can also deregulate the cell cycle by binding and abrogating the inhibitory effect of CDK inhibitors like p21Cip1 and p27Kip1. This interaction abolishes their inhibitory effect on the cell cycle and highlights a different mechanism by which E7 favors cell proliferation (262–264).
E6 and E7 of HPV38 and 49, which can efficiently immortalize human primary keratinocytes, do not have the ability to target pRb for degradation when expressed in these cells (120). However, similarly to E7 of LR mucosal types HPV38 and 49 E7s can efficiently bind to pRb. In contrast, beta HPV8 E7 binds to pRb to a lesser extent compared with HPV16 E7, but it is still able to repress pRb levels when expressed in keratinocytes (265, 266).
Studies from our laboratory have shown that different beta HPV types (e.g., HPV38) are able to trigger pRb hyperphosphorylation, an event that leads to release of E2F and to the transcriptional activation of E2F target genes [reviewed in (19)]. In addition, E7 of HPV38 shows the ability to counteract p53-mediated apoptosis by inducing accumulation of the p73 isoform, ΔNp73, by two different mechanisms: (i) by inducing the recruitment of a post-translationally modified form of p53 to the promoter of ΔNp73, which activates its transcription (267), and (ii) by mediating ΔNp73 protein stabilization with an IκB kinase beta (IKKβ)-mediated mechanism (268). The latter studies highlight that cutaneous HPVs have evolved different mechanisms than HR mucosal HPVs, to target key cellular proteins (pRb and p53) and promote cellular transformation.
Transforming Activities of E5
Although the transforming activities of E6 and E7 oncoproteins are well characterized, the role of the E5 oncoprotein is still poorly understood (269). E5 is a hydrophobic membrane-associated protein that localizes to the endoplasmic reticulum, Golgi apparatus, and plasma membrane (270, 271). E5 forms a hexameric viroporin complex, which modulates ion homeostasis (269). E5 viroporin is required for the hyperactivation of mitotic signaling, and it plays a key role during the life cycle of HR HPV types. Recombinant HPV16 E5 membrane channel activity is sensitive to compounds such as adamantine and rimantadine. Inhibition of E5 viroporin activity in primary human keratinocytes harboring HPV18 genome leads to decreased ERK-MAPK phosphorylation and cyclin B1 levels, and to an increase in the expression of differentiation marker (272). E5 viroporin activity is therefore critical for maintaining mitogenic signaling and delaying expression of differentiation marker during the productive stages of the HPV18 life cycle; therefore, targeting viroporin could pave the way for an effective antiviral strategy.
The expression of E5 increased the immortalization effects of HPV16 E6 and E7 (273). E5 is involved in the transformation of primary human keratinocytes by potentiating epidermal growth factor receptor (EGFR) signaling (274). In 16E5 transgenic mice, EGFR is required for epithelial hyperplasia (275). The activation of EGFR leads to induction of the expression of cyclooxygenase-2 (COX-2), which plays an important role in cervical carcinogenesis, through nuclear factor-kappa B (NF-κB) and AP-1 (276). Via stimulation of EGFR, HPV16 E5 induces Met, a growth factor receptor involved in tumor cell motility and cancer metastasis; moreover, Met signaling is required for proper differentiation-dependent viral gene expression (277).
HPV16 E5 protein enhances, in cooperation with EGFR signaling, the down-regulation of tumor suppressor p27Kip1 levels, which promotes cell-cycle progression (278). The level of another CDK inhibitor, p21, is reduced in cells that over-express E5 (279). Both events lead to an increase in cell proliferation and in the percentage of cells in the S phase.
HPV16 E5 also contributes to carcinogenesis by inhibiting apoptosis by multiple molecular strategies. The effect of HPV16 E5 on apoptosis was assessed using raft cultures in which E5 prevented FasL- or TRAIL-mediated apoptosis (280). A recent study reported that HPV16 E5 specifically down-regulates miR-196a in cervical cancer cells. The same study showed that miR-196a expression affects cell proliferation, cell growth, and apoptosis, thus confirming the role of HPV16 E5 in cervical cancer development (281). Another anti-apoptotic strategy consists of decreasing pro-apoptotic Bak and Bax levels, and increasing the expression of the anti-apoptotic Bcl-2. HPV16 E5 inhibits apoptosis of cervical cancer cells by stimulating the ubiquitin-proteasome-mediated degradation of Bax protein (282). It has also been reported that HPV16 E5 protects E5-expressing foreskin keratinocytes from UVB radiation-induced apoptosis (283).
All these data indicate that expression of E5 contributes to HPV-driven cellular transformation. However, unlike E6 and E7, E5 is not required for the maintenance of the oncogenic phenotype, because the E5 ORF is frequently disrupted during HPV integration, thus rather highlighting a role of E5 at the initial step of carcinogenesis (284).
HPV types from different genera (beta, gamma, and mu) lack the E5 ORF, suggesting that the corresponding functions are performed by other viral proteins.
E5 forms a complex with EVER1 and EVER2, proteins encoded by the TMC6 and TMC8 genes, respectively, (92). Most of the patients with EV who are susceptible to persistent beta HPV infections carry mutations in these genes. The EVER proteins interact with the zinc transporter ZnT-1 in order to regulate intracellular zinc homeostasis. Zinc sequestration by the ZnT-1/EVER complex may inhibit the replication of EV-associated HPV types. In contrast, mutations in EVER proteins lead to the disruption of the ZnT-1/EVER complex, and an increase in the level of zinc in the cytoplasm, resulting in increased replication of EV-associated HPV types (285), [reviewed in (286)]. A recent study reported that patients with EV can also display the typical EV phenotype, bearing non-synonymous calcium- and integrin-binding protein 1 (CIB1) mutations and lacking TMC6 and TMC8 mutations (287). CIB1 does not seem to interact with ZnT-1, but this protein forms a complex with EVER1 and EVER2. The CIB1-EVER1-EVER2 complex appears to act as a restriction factor for infection by cutaneous HPVs that do not express E5. In patients with EV, impairment of CIB1, EVER1, or EVER2 allows replication of these HPV types and development of EV lesions on the skin (287).
Conclusive Remarks, Challenges, and Open Questions
In the past decades, epidemiological and biological studies have clearly demonstrated the role of mucosal HR HPV types in human carcinogenesis. While co-evolving with the human host, HPVs have developed strategies to evade the immune system and create the optimal conditions to persist in the host for many years. During HPV chronic infection, the E6 and E7 oncoproteins interact with different cellular proteins to prevent apoptosis, counteract the cellular senescence program, and promote unscheduled cell growth. These activities of HR HPV E6 and E7 are part of a viral strategy to complete the viral life cycle and generate progeny. The expression of the oncoproteins is kept under the tight control of E2. Loss of E2 repressive functions, by viral integration or by other mechanisms, which are discussed in this review, leads to constitutive high expression of E6 and E7, a key event in HPV-mediated transformation (Figure 3). Although the mechanisms leading to cervical cancer have been well-characterized, there is still a gap in the understanding of the pathogenesis of other HPV-related malignancies, e.g., oropharyngeal cancer. The proportion of HPV-driven oropharyngeal cancers has been increasing over the past decades in many parts of the world, mainly in Europe and North America (288–290). The natural history of HPV infection in the oral cavity and oropharynx is poorly studied and needs further investigation.
Of note, recent studies identified HPV16 E6 antibody as a good diagnostic and prognostic marker of oropharyngeal cancer. HPV16 E6 seropositivity was detected more than 10 years before the diagnosis (291, 292).
Moreover, although the available prophylactic vaccines have demonstrated their effectiveness in the context of prevention of anogenital cancers (293), very limited epidemiological data are available on the efficacy of the HPV vaccine on oral HPV infection (294). This deserves further investigation.
Prophylactic vaccines show a good efficiency in covering targeted HPVs, but they are inefficient in existing infections. Therefore, therapeutic vaccines are needed to fill this gap. Several studies have reported the development of promising therapeutic vaccines. Their use in the near future could be expected; however, their cost could be high. This limitation could be overcome by the use of inexpensive technologies (i.e., plant-made therapeutics) (295–299).
One of the most recurrent questions regarding the prophylactic vaccine is whether HPV vaccination will result in the occupation of a newly vacant ecological niche by non-targeted HPV types. This hypothesis has been referred as “type replacement. ” To date, the available data do not support this hypothesis. However, as reported in recent studies, a few types (i.e., HPV51 and 52) may need to be monitored after vaccination (300, 301).
The bacterial microbiome composition may also play a role in the outcome of mucosal HPV infection. A chronic mycoplasma infection was shown to promote cervical dysplasia induced by HPV (302). Therefore, identifying the bacterial microbiome could allow the identification of patients at high risk for developing cervical cancer.
Our knowledge of the biology of beta and gamma HPV types is still uncomplete. For example, the expression patterns of cutaneous HPV types and associated splicing patterns are largely unknown, except for HPV5 and 8 (https://pave.niaid.nih.gov/). There are also important gaps in the understanding of the life cycle of cutaneous HPV types, which also deserve further research.
In February, 2009, an IARC Monograph Working Group met at IARC to assess the carcinogenicity of biological agents including beta and gamma HPV types. HPV5 and 8 were classified as possibly carcinogenic to humans in patients with EV (IARC Group 2B), whereas the other cutaneous HPV types were each determined to be “not classifiable as to its carcinogenicity to humans” (IARC Group 3), on the basis of insufficient epidemiological and biological data available at that time (10). In fact, although several epidemiological and biological studies support a role of β1 and β2 HPV types in NMSC, the lack of cutaneous HPV E6 and E7 expression in SCC has raised the question of whether the viral proteins are required in the carcinogenesis process.
Since that 2009 meeting, strong epidemiological and experimental data in favor of a role of cutaneous HPV types as cofactors in skin carcinogenesis have been published, as described in this review. Indeed, recent findings from our group and others have shown that unlike the expression of oncoproteins of HR HPV types, cutaneous HPV E6 and E7 expression appears to be required only at the initial step of skin carcinogenesis by exacerbating the deleterious effects of UV radiation (104). Over the past decade, studies using experimental animals (HPV38 transgenic mice and the MnPV-infected rodent model Mastomys coucha) have provided strong in vivo evidence for a synergism between UV radiation and cutaneous HPV types in the development of cSCC (104, 107). Most importantly, both models provided further evidence for a hit-and-run mechanism. However, additional studies are still necessary to validate this hypothesis in humans.
In addition to these experimental data, a recent prospective study of two organ transplant recipient cohorts showed that having five or more different beta HPV types in eyebrow hair, and a high beta HPV viral load, was associated with cSCC carcinogenesis, providing strong evidence for a role of beta HPVs in this cancer, whereas no evidence was found for basal cell carcinoma (303). The development of vaccination strategies that target cutaneous HPV infections would be of great benefit to these patients (304, 305).
Another open question regards the tropism of certain HPVs. As previously specified in this review, a dual tropism of HPV species β3 has been proposed. Moreover, this species displays biological similarities with mucosal HR HPV types (111, 120). Our recent data highlighted that E6 and E7 of HPV49 share some features with HR mucosal types, and efficiently transform epithelial cells in in vitro and in vivo models (120, 121). In addition, in line with these biological data, epidemiological studies have shown β3 HPVs and other cutaneous types in the oral cavity mucosa. Beta and gamma HPV types have also been isolated from the nasopharynx and the nasal mucosa (119, 306). Moreover, a recent study based on a prospective design reported that beta or gamma HPV types were associated with the incidence of HNC (117). Those findings warrant further studies to better understand the role of those so-called “cutaneous HPV types” in mucosotropic cancers (111, 119).
A large number of gamma HPV types have been isolated; however, their biological activities are poorly studied and deserve more investigation. Finally, some of the known but not yet biologically characterized beta or gamma HPV types may display transforming activity and also deserve attention. In addition to these not-yet-studied types, one may also consider the large number of novel HPVs found in recent years, as well as the rapid evolution of viral detection and screening techniques and large data set analysis pipelines.
For all these reasons, we believe that in the coming years the field of HPV research will continue to be a source of surprises and excitement for the scientific community.
Author Contributions
The author confirms being the sole contributor of this work and has approved it for publication. TG wrote the manuscript and designed the figures.
Disclaimer
The author is staff member of the International Agency for Research on Cancer / World Health Organization. The author alone is responsible for the views expressed in this article and he does not necessarily represent the decisions, policy or views of the International Agency for Research on Cancer/World Health Organization.
Conflict of Interest Statement
The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
I would like to thank Dr. Massimo Tommasino and Dr. Rosita Accardi for helpful discussions. I am grateful to Dr. Karen Müller and Ms. Jessica Cox for editing, and to Ms. Nicole Suty for her help with preparation of this manuscript.
References
1. Bravo IG, de Sanjose S, Gottschling M. The clinical importance of understanding the evolution of papillomaviruses. Trends Microbiol. (2010) 18:432–8. doi: 10.1016/j.tim.2010.07.008
2. Lopez-Bueno A, Mavian C, Labella AM, Castro D, Borrego JJ, Alcami A, et al. Concurrence of iridovirus, polyomavirus, and a unique member of a new group of fish papillomaviruses in lymphocystis disease-affected gilthead sea bream. J Virol. (2016) 90:8768–79. doi: 10.1128/JVI.01369-16
3. Bravo IG, Felez-Sanchez M. Papillomaviruses: Viral evolution, cancer and evolutionary medicine. Evol Med Public Health. (2015) 2015:32–51. doi: 10.1093/emph/eov003
4. Gottschling M, Goker M, Stamatakis A, Bininda-Emonds OR, Nindl I, Bravo IG. Quantifying the phylodynamic forces driving papillomavirus evolution. Mol Biol Evol. (2011) 28:2101–13. doi: 10.1093/molbev/msr030
5. de Villiers EM, Fauquet C, Broker TR, Bernard HU, zur Hausen H. Classification of papillomaviruses. Virology. (2004) 324:17–27. doi: 10.1016/j.virol.2004.03.033
6. Muhr LSA, Eklund C, Dillner J. Towards quality and order in human papillomavirus research. Virology. (2018) 519:74–6. doi: 10.1016/j.virol.2018.04.003
7. Brancaccio RN, Robitaille A, Dutta S, Cuenin C, Santare D, Skenders G, et al. Generation of a novel next-generation sequencing-based method for the isolation of new human papillomavirus types. Virology. (2018) 520:1–10. doi: 10.1016/j.virol.2018.04.017
8. Pastrana DV, Peretti A, Welch NL, Borgogna C, Olivero C, Badolato R, et al. Metagenomic Discovery of 83 new human papillomavirus types in patients with immunodeficiency. mSphere. (2018) 3:e00645–18. doi: 10.1128/mSphereDirect.00645-18
9. Halec G, Schmitt M, Dondog B, Sharkhuu E, Wentzensen N, Gheit T, et al. Biological activity of probable/possible high-risk human papillomavirus types in cervical cancer. Int J Cancer. (2013) 132:63–71. doi: 10.1002/ijc.27605
10. Bouvard V, Baan R, Straif K, Grosse Y, Secretan B, El Ghissassi F, et al. A review of human carcinogens–Part B: biological agents. Lancet Oncol. (2009) 10:321–2. doi: 10.1016/S1470-2045(09)70096-8
11. Haedicke J, Iftner T. Human papillomaviruses and cancer. Radiother Oncol. (2013) 108:397–402. doi: 10.1016/j.radonc.2013.06.004
12. Plummer M, de Martel C, Vignat J, Ferlay J, Bray F, Franceschi S. Global burden of cancers attributable to infections in 2012: a synthetic analysis. Lancet Global Health. (2016) 4:e609–16. doi: 10.1016/S2214-109X(16)30143-7
13. Ball SL, Winder DM, Vaughan K, Hanna N, Levy J, Sterling JC, et al. Analyses of human papillomavirus genotypes and viral loads in anogenital warts. J Med Virol. (2011) 83:1345–50. doi: 10.1002/jmv.22111
14. Derkay CS, Wiatrak B. Recurrent respiratory papillomatosis: a review. Laryngoscope. (2008) 118:1236–47. doi: 10.1097/MLG.0b013e31816a7135
15. McLaughlin-Drubin ME, Meyers J, Munger K. Cancer associated human papillomaviruses. Curr Opin Virol. (2012) 2:459–66. doi: 10.1016/j.coviro.2012.05.004
16. King CM, Johnston JS, Ofili K, Tam M, Palefsky J, Da Costa M, et al. Human papillomavirus types 2, 27, and 57 Identified in plantar verrucae from HIV-positive and HIV-negative individuals. J Am Pediatric Med Assoc. (2014) 104:141–6. doi: 10.7547/0003-0538-104.2.141
17. Hogendoorn GK, Bruggink SC, de Koning MNC, Eekhof JAH, Hermans KE, Rissmann R, et al. Morphological characteristics and human papillomavirus genotype predict the treatment response in cutaneous warts. Br J Dermatol. (2018) 178:253–60. doi: 10.1111/bjd.15758
18. Bruggink SC, de Koning MN, Gussekloo J, Egberts PF, Ter Schegget J, Feltkamp MC, et al. Cutaneous wart-associated HPV types: prevalence and relation with patient characteristics. J Clin Virol. (2012) 55:250–5. doi: 10.1016/j.jcv.2012.07.014
19. Accardi R, Gheit T. Cutaneous HPV and skin cancer. Presse Med. (2014) 43:e435–43. doi: 10.1016/j.lpm.2014.08.008
20. Viarisio D, Gissmann L, Tommasino M. Human papillomaviruses and carcinogenesis: well-established and novel models. Curr Opin Virol. (2017) 26:56–62. doi: 10.1016/j.coviro.2017.07.014
21. Tommasino M. The biology of beta human papillomaviruses. Virus Res. (2017) 231:128–38. doi: 10.1016/j.virusres.2016.11.013
22. Van Doorslaer K, McBride AA. Molecular archeological evidence in support of the repeated loss of a papillomavirus gene. Sci Rep. (2016) 6:33028. doi: 10.1038/srep33028
23. Chen Z, Schiffman M, Herrero R, Desalle R, Burk RD. Human papillomavirus (HPV) types 101 and 103 isolated from cervicovaginal cells lack an E6 open reading frame (ORF) and are related to gamma-papillomaviruses. Virology. (2007) 360:447–53. doi: 10.1016/j.virol.2006.10.022
24. Nobre RJ, Herraez-Hernandez E, Fei JW, Langbein L, Kaden S, Grone HJ, et al. E7 oncoprotein of novel human papillomavirus type 108 lacking the E6 gene induces dysplasia in organotypic keratinocyte cultures. J Virol. (2009) 83:2907–16. doi: 10.1128/JVI.02490-08
25. Bravo IG, Alonso A. Mucosal human papillomaviruses encode four different E5 proteins whose chemistry and phylogeny correlate with malignant or benign growth. J Virol. (2004) 78:13613–26. doi: 10.1128/JVI.78.24.13613-13626.2004
26. McBride AA. The papillomavirus E2 proteins. Virology. (2013) 445:57–79. doi: 10.1016/j.virol.2013.06.006
27. Zobel T, Iftner T, Stubenrauch F. The papillomavirus E8-E2C protein represses DNA replication from extrachromosomal origins. Mol Cell Biol. (2003) 23:8352–62. doi: 10.1128/MCB.23.22.8352-8362.2003
28. Lace MJ, Anson JR, Thomas GS, Turek LP, Haugen TH. The E8–E2 gene product of human papillomavirus type 16 represses early transcription and replication but is dispensable for viral plasmid persistence in keratinocytes. J Virol. (2008) 82:10841–53. doi: 10.1128/JVI.01481-08
29. Ammermann I, Bruckner M, Matthes F, Iftner T, Stubenrauch F. Inhibition of transcription and DNA replication by the papillomavirus E8-E2C protein is mediated by interaction with corepressor molecules. J Virol. (2008) 82:5127–36. doi: 10.1128/JVI.02647-07
30. Burchell AN, Winer RL, de Sanjose S, Franco EL. Chapter 6: Epidemiology and transmission dynamics of genital HPV infection. Vaccine. (2006) 24 (Suppl. 3):S3/52–61. doi: 10.1016/j.vaccine.2006.05.031
31. Ryndock EJ, Meyers C. A risk for non-sexual transmission of human papillomavirus? Exp Rev Anti-Infect Ther. (2014) 12:1165–70. doi: 10.1586/14787210.2014.959497
32. Schiffman M, Castle PE, Jeronimo J, Rodriguez AC, Wacholder S. Human papillomavirus and cervical cancer. Lancet. (2007) 370:890–907. doi: 10.1016/S0140-6736(07)61416-0
33. Maucort-Boulch D, Plummer M, Castle PE, Demuth F, Safaeian M, Wheeler CM, et al. Predictors of human papillomavirus persistence among women with equivocal or mildly abnormal cytology. Int J Cancer. (2010) 126:684–91. doi: 10.1002/ijc.24752
34. Oh HY, Seo SS, Kim MK, Lee DO, Chung YK, Lim MC, et al. Synergistic effect of viral load and alcohol consumption on the risk of persistent high-risk human papillomavirus infection. PLoS ONE. (2014) 9:e104374. doi: 10.1371/journal.pone.0104374
35. Giuliano AR, Sedjo RL, Roe DJ, Harri R, Baldwi S, Papenfuss MR, et al. Clearance of oncogenic human papillomavirus (HPV) infection: effect of smoking (United States). Cancer Causes Control. (2002) 13:839–46. doi: 10.1023/A:1020668232219
36. Haukioja A, Asunta M, Soderling E, Syrjanen S. Persistent oral human papillomavirus infection is associated with smoking and elevated salivary immunoglobulin G concentration. J Clin Virol. (2014) 61:101–6. doi: 10.1016/j.jcv.2014.06.012
37. Magnusson PK, Sparen P, Gyllensten UB. Genetic link to cervical tumours. Nature. (1999) 400:29–30. doi: 10.1038/21801
38. Hemminki K, Dong C, Vaittinen P. Familial risks in cervical cancer: is there a hereditary component? Int J Cancer. (1999) 82:775–81. doi: 10.1002/(SICI)1097-0215(19990909)82:6<775::AID-IJC1>3.0.CO;2-V
39. Chattopadhyay K. A comprehensive review on host genetic susceptibility to human papillomavirus infection and progression to cervical cancer. Indian J Human Genet. (2011) 17:132–44. doi: 10.4103/0971-6866.92087
40. Wang SS, Gonzalez P, Yu K, Porras C, Li Q, Safaeian M, et al. Common genetic variants and risk for HPV persistence and progression to cervical cancer. PLoS ONE. (2010) 5:e8667. doi: 10.1371/journal.pone.0008667
41. Pontillo A, Bricher P, Leal VN, Lima S, Souza PR, Crovella S. Role of inflammasome genetics in susceptibility to HPV infection and cervical cancer development. J Med Virol. (2016) 88:1646–51. doi: 10.1002/jmv.24514
42. Bernal-Silva S, Granados J, Gorodezky C, Alaez C, Flores-Aguilar H, Cerda-Flores RM, et al. HLA-DRB1 Class II antigen level alleles are associated with persistent HPV infection in Mexican women; a pilot study. Infect Agent Cancer. (2013) 8:31. doi: 10.1186/1750-9378-8-31
43. Chen AA, Heideman DA, Boon D, Gheit T, Snijders PJ, Tommasino M, et al. Human papillomavirus 45 genetic variation and cervical cancer risk worldwide. J Virol. (2014) 88:4514–21. doi: 10.1128/JVI.03534-13
44. Chen AA, Heideman DA, Boon D, Chen Z, Burk RD, De Vuyst H, et al. Human papillomavirus 33 worldwide genetic variation and associated risk of cervical cancer. Virology. (2014) 448:356–62. doi: 10.1016/j.virol.2013.10.033
45. Cornet I, Gheit T, Iannacone MR, Vignat J, Sylla BS, Del Mistro A, et al. HPV16 genetic variation and the development of cervical cancer worldwide. Br J Cancer. (2013) 108:240–4. doi: 10.1038/bjc.2012.508
46. Chen AA, Gheit T, Franceschi S, Tommasino M, Clifford GM, Group IHVS. Human Papillomavirus 18 genetic variation and cervical cancer risk worldwide. J Virol. (2015) 89:10680–7. doi: 10.1128/JVI.01747-15
47. Gheit T, Cornet I, Clifford GM, Iftner T, Munk C, Tommasino M, et al. Risks for persistence and progression by human papillomavirus type 16 variant lineages among a population-based sample of Danish women. Cancer Epidemiol Biomarkers Prev. (2011) 20:1315–21. doi: 10.1158/1055-9965.EPI-10-1187
48. Cornet I, Gheit T, Clifford GM, Combes JD, Dalstein V, Franceschi S, et al. Human papillomavirus type 16 E6 variants in France and risk of viral persistence. Infect Agent Cancer. (2013) 8:4. doi: 10.1186/1750-9378-8-4
49. Mirabello L, Yeager M, Yu K, Clifford GM, Xiao Y, Zhu B, et al. HPV16 E7 genetic conservation is critical to carcinogenesis. Cell. (2017) 170:1164–74 e6. doi: 10.1016/j.cell.2017.08.001
50. Westrich JA, Warren CJ, Klausner MJ, Guo K, Liu CW, Santiago ML, et al. Human Papillomavirus 16 E7 Stabilizes APOBEC3A protein by inhibiting Cullin 2-dependent protein degradation. J Virol. (2018) 92:e01318–17. doi: 10.1128/JVI.01318-17
51. zur Hausen H. Papillomaviruses causing cancer: evasion from host-cell control in early events in carcinogenesis. J Natl Cancer Inst. (2000) 92:690–8. doi: 10.1093/jnci/92.9.690
52. Peitsaro P, Johansson B, Syrjanen S. Integrated human papillomavirus type 16 is frequently found in cervical cancer precursors as demonstrated by a novel quantitative real-time PCR technique. J Clin Microbiol. (2002) 40:886–91. doi: 10.1128/JCM.40.3.886-891.2002
53. Collins SI, Constandinou-Williams C, Wen K, Young LS, Roberts S, Murray PG, et al. Disruption of the E2 gene is a common and early event in the natural history of cervical human papillomavirus infection: a longitudinal cohort study. Cancer Res. (2009) 69:3828–32. doi: 10.1158/0008-5472.CAN-08-3099
54. Jeon S, Lambert PF. Integration of human papillomavirus type 16 DNA into the human genome leads to increased stability of E6 and E7 mRNAs: implications for cervical carcinogenesis. Proc Natl Acad Sci USA. (1995) 92:1654–8. doi: 10.1073/pnas.92.5.1654
55. Cancer Genome Atlas Research N, Albert Einstein College of M, Analytical Biological S, Barretos Cancer H Baylor College of M, Beckman Research Institute of C,ity of H et al. Integrated genomic and molecular characterization of cervical cancer. Nature. (2017) 543:378–84. doi: 10.1038/nature21386
56. Xu B, Chotewutmontri S, Wolf S, Klos U, Schmitz M, Durst M, et al. Multiplex Identification of human papillomavirus 16 DNA integration sites in cervical carcinomas. PLoS ONE. (2013) 8:e66693. doi: 10.1371/journal.pone.0066693
57. Gray E, Pett MR, Ward D, Winder DM, Stanley MA, Roberts I, et al. In vitro progression of human papillomavirus 16 episome-associated cervical neoplasia displays fundamental similarities to integrant-associated carcinogenesis. Cancer Res. (2010) 70:4081–91. doi: 10.1158/0008-5472.CAN-09-3335
58. Hong D, Liu J, Hu Y, Lu X, Li B, Li Y, et al. Viral E6 is overexpressed via high viral load in invasive cervical cancer with episomal HPV16. BMC Cancer. (2017) 17:136. doi: 10.1186/s12885-017-3124-9
59. Zhang R, Shen C, Zhao L, Wang J, McCrae M, Chen X, et al. Dysregulation of host cellular genes targeted by human papillomavirus (HPV) integration contributes to HPV-related cervical carcinogenesis. Int J Cancer. (2016) 138:1163–74. doi: 10.1002/ijc.29872
60. Shen C, Liu Y, Shi S, Zhang R, Zhang T, Xu Q, et al. Long-distance interaction of the integrated HPV fragment with MYC gene and 8q24.22 region upregulating the allele-specific MYC expression in HeLa cells. Int J Cancer. (2017) 141:540–8. doi: 10.1002/ijc.30763
61. Glover TW, Berger C, Coyle J, Echo B. DNA polymerase alpha inhibition by aphidicolin induces gaps and breaks at common fragile sites in human chromosomes. Hum Genet. (1984) 67:136–42. doi: 10.1007/BF00272988
62. Thorland EC, Myers SL, Gostout BS, Smith DI. Common fragile sites are preferential targets for HPV16 integrations in cervical tumors. Oncogene. (2003) 22:1225–37. doi: 10.1038/sj.onc.1206170
63. Christiansen IK, Sandve GK, Schmitz M, Durst M, Hovig E. Transcriptionally active regions are the preferred targets for chromosomal HPV integration in cervical carcinogenesis. PLoS ONE. (2015) 10:e0119566. doi: 10.1371/journal.pone.0119566
64. Schmitz M, Driesch C, Jansen L, Runnebaum IB, Durst M. Non-random integration of the HPV genome in cervical cancer. PLoS ONE. (2012) 7:e39632. doi: 10.1371/journal.pone.0039632
65. Leung TW, Liu SS, Leung RC, Chu MM, Cheung AN, Ngan HY. HPV 16 E2 binding sites 1 and 2 become more methylated than E2 binding site 4 during cervical carcinogenesis. J Med Virol. (2015) 87:1022–33. doi: 10.1002/jmv.24129
66. Oka N, Kajita M, Nishimura R, Ohbayashi C, Sudo T. L1 gene methylation in high-risk human papillomaviruses for the prognosis of cervical intraepithelial neoplasia. Int J Gynecol Cancer. (2013) 23:235–43. doi: 10.1097/IGC.0b013e31827da1f6
67. Torres-Rojas FI, Alarcon-Romero LDC, Leyva-Vazquez MA, Ortiz-Ortiz J, Mendoza-Catalan MA, Hernandez-Sotelo D, et al. Methylation of the L1 gene and integration of human papillomavirus 16 and 18 in cervical carcinoma and premalignant lesions. Oncol Lett. (2018) 15:2278–86. doi: 10.3892/ol.2017.7596
68. Parfenov M, Pedamallu CS, Gehlenborg N, Freeman SS, Danilova L, Bristow CA, et al. Characterization of HPV and host genome interactions in primary head and neck cancers. Proc Natl Acad Sci USA. (2014) 111:15544–9. doi: 10.1073/pnas.1416074111
69. Olthof NC, Speel EJ, Kolligs J, Haesevoets A, Henfling M, Ramaekers FC, et al. Comprehensive analysis of HPV16 integration in OSCC reveals no significant impact of physical status on viral oncogene and virally disrupted human gene expression. PLoS ONE. (2014) 9:e88718. doi: 10.1371/journal.pone.0088718
70. Vojtechova Z, Sabol I, Salakova M, Turek L, Grega M, Smahelova J, et al. Analysis of the integration of human papillomaviruses in head and neck tumours in relation to patients' prognosis. Int J Cancer. (2016) 138:386–95. doi: 10.1002/ijc.29712
71. Nulton TJ, Olex AL, Dozmorov M, Morgan IM, Windle B. Analysis of the cancer genome Atlas sequencing data reveals novel properties of the human papillomavirus 16 genome in head and neck squamous cell carcinoma. Oncotarget. (2017) 8:17684–99. doi: 10.18632/oncotarget.15179
72. Antonsson A, Forslund O, Ekberg H, Sterner G, Hansson BG. The ubiquity and impressive genomic diversity of human skin papillomaviruses suggest a commensalic nature of these viruses. J Virol. (2000) 74:11636–41. doi: 10.1128/JVI.74.24.11636-11641.2000
73. de Koning MN, Weissenborn SJ, Abeni D, Bouwes Bavinck JN, Euvrard S, Green AC, et al. Prevalence and associated factors of betapapillomavirus infections in individuals without cutaneous squamous cell carcinoma. J Gen Virol. (2009) 90:1611–21. doi: 10.1099/vir.0.010017-0
74. Bouwes Bavinck JN, Neale RE, Abeni D, Euvrard S, Green AC, Harwood CA, et al. Multicenter study of the association between betapapillomavirus infection and cutaneous squamous cell carcinoma. Cancer Res. (2010) 70:9777–86. doi: 10.1158/0008-5472.CAN-10-0352
75. Hsu JY, Chen AC, Keleher A, McMillan NA, Antonsson A. Shared and persistent asymptomatic cutaneous human papillomavirus infections in healthy skin. J Med Virol. (2009) 81:1444–9. doi: 10.1002/jmv.21529
76. Antonsson A, Karanfilovska S, Lindqvist PG, Hansson BG. General acquisition of human papillomavirus infections of skin occurs in early infancy. J Clin Microbiol. (2003) 41:2509–14. doi: 10.1128/JCM.41.6.2509-2514.2003
77. Plasmeijer EI, Struijk L, Bouwes Bavinck JN, Feltkamp MC. Epidemiology of cutaneous human papillomavirus infections. Cancer Treat Res. (2009) 146:143–57. doi: 10.1007/978-0-387-78574-5_13
78. Egawa N, Egawa K, Griffin H, Doorbar J. Human papillomaviruses; epithelial tropisms, and the development of neoplasia. Viruses. (2015) 7:3863–90. doi: 10.3390/v7072802
79. Kohler A, Forschner T, Meyer T, Ulrich C, Gottschling M, Stockfleth E, et al. Multifocal distribution of cutaneous human papillomavirus types in hairs from different skin areas. Br J Dermatol. (2007) 156:1078–80. doi: 10.1111/j.1365-2133.2007.07809.x
80. Schmitt A, Rochat A, Zeltner R, Borenstein L, Barrandon Y, Wettstein FO, et al. The primary target cells of the high-risk cottontail rabbit papillomavirus colocalize with hair follicle stem cells. J Virol. (1996) 70:1912–22.
81. Rous P, Beard JW. The progression to carcinoma of virus-induced rabbit papillomas (Shope). J Exp Med. (1935) 62:523–48. doi: 10.1084/jem.62.4.523
82. Plasmeijer EI, Neale RE, Buettner PG, de Koning MN, Ter Schegget J, Quint WG, et al. Betapapillomavirus infection profiles in tissue sets from cutaneous squamous cell-carcinoma patients. Int J Cancer. (2010) 126:2614–21. doi: 10.1002/ijc.24991
83. Weissenborn SJ, De Koning MN, Wieland U, Quint WG, Pfister HJ. Intrafamilial transmission and family-specific spectra of cutaneous betapapillomaviruses. J Virol. (2009) 83:811–6. doi: 10.1128/JVI.01338-08
84. Hazard K, Karlsson A, Andersson K, Ekberg H, Dillner J, Forslund O. Cutaneous human papillomaviruses persist on healthy skin. J Invest Dermatol. (2007) 127:116–9. doi: 10.1038/sj.jid.5700570
85. de Koning MN, Struijk L, Bavinck JN, Kleter B, ter Schegget J, Quint WG, et al. Betapapillomaviruses frequently persist in the skin of healthy individuals. J Gen Virol. (2007) 88:1489–95. doi: 10.1099/vir.0.82732-0
86. Moscicki AB, Ma Y, Gheit T, McKay-Chopin S, Farhat S, Widdice LE, et al. Prevalence and transmission of beta and gamma human papillomavirus in heterosexual couples. Open Forum Infect Dis. (2017) 4:ofw216. doi: 10.1093/ofid/ofw216
87. Rahman S, Pierce Campbell CM, Waterboer T, Rollison DE, Ingles DJ, Torres BN, et al. Seroprevalence of cutaneous human papillomaviruses (HPVs) among men in the multinational HPV Infection in Men study. J Gen Virol. (2016) 97:3291–301. doi: 10.1099/jgv.0.000620
88. Michael KM, Waterboer T, Pfister H, Gariglio M, Majewski S, Favre M, et al. Seroreactivity of 38 human papillomavirus types in epidermodysplasia verruciformis patients, relatives, and controls. J Invest Dermatol. (2010) 130:841–8. doi: 10.1038/jid.2009.356
89. Weissenborn S, Neale RE, Waterboer T, Abeni D, Bavinck JN, Green AC, et al. Beta-papillomavirus DNA loads in hair follicles of immunocompetent people and organ transplant recipients. Med Microbiol Immunol. (2012) 201:117–25. doi: 10.1007/s00430-011-0212-3
90. Dillner J. The serological response to papillomaviruses. Semin Cancer Biol. (1999) 9:423–30. doi: 10.1006/scbi.1999.0146
91. Michael KM, Waterboer T, Sehr P, Rother A, Reidel U, Boeing H, et al. Seroprevalence of 34 human papillomavirus types in the German general population. PLoS Pathog. (2008) 4:e1000091. doi: 10.1371/journal.ppat.1000091
92. Orth G. Host defenses against human papillomaviruses: lessons from epidermodysplasia verruciformis. Curr Topics Microbiol Immunol. (2008) 321:59–83. doi: 10.1007/978-3-540-75203-5_3
93. Orth G, Jablonska S, Favre M, Croissant O, Jarzabek-Chorzelska M, Rzesa G. Characterization of two types of human papillomaviruses in lesions of epidermodysplasia verruciformis. Proc Natl Acad Sci USA. (1978) 75:1537–41. doi: 10.1073/pnas.75.3.1537
94. O'Reilly Zwald F, Brown M. Skin cancer in solid organ transplant recipients: advances in therapy and management: part I. epidemiology of skin cancer in solid organ transplant recipients. J Am Acad Dermatol. (2011) 65:253–61. doi: 10.1016/j.jaad.2010.11.062
95. Nindl I, Rosl F. Molecular concepts of virus infections causing skin cancer in organ transplant recipients. Am J Transl. (2008) 8:2199–204. doi: 10.1111/j.1600-6143.2008.02392.x
96. Silverberg MJ, Leyden W, Warton EM, Quesenberry CP Jr, Engels EA, Asgari MM. HIV infection status, immunodeficiency, and the incidence of non-melanoma skin cancer. J Natl Cancer Inst. (2013) 105:350–60. doi: 10.1093/jnci/djs529
97. Asgari MM, Ray GT, Quesenberry CP Jr, Katz KA, Silverberg MJ. Association of multiple primary skin cancers with human immunodeficiency virus infection, CD4 count, and viral load. JAMA Dermatol. (2017) 153:892–6. doi: 10.1001/jamadermatol.2017.1716
98. Grulich AE, van Leeuwen MT, Falster MO, Vajdic CM. Incidence of cancers in people with HIV/AIDS compared with immunosuppressed transplant recipients: a meta-analysis. Lancet. (2007) 370:59–67. doi: 10.1016/S0140-6736(07)61050-2
99. Howley PM, Pfister HJ. Beta genus papillomaviruses and skin cancer. Virology. (2015) 479–480:290–6. doi: 10.1016/j.virol.2015.02.004
100. Iannacone MR, Gheit T, Pfister H, Giuliano AR, Messina JL, Fenske NA, et al. Case-control study of genus-beta human papillomaviruses in plucked eyebrow hairs and cutaneous squamous cell carcinoma. Int J Cancer. (2014) 134:2231–44. doi: 10.1002/ijc.28552
101. Iannacone MR, Gheit T, Waterboer T, Giuliano AR, Messina JL, Fenske NA, et al. Case-control study of cutaneous human papillomaviruses in squamous cell carcinoma of the skin. Cancer Epidemiol Biomarkers Prev. (2012) 21:1303–13. doi: 10.1158/1055-9965.EPI-12-0032
102. Faust H, Andersson K, Luostarinen T, Gislefoss RE, Dillner J. Cutaneous human papillomaviruses and squamous cell carcinoma of the skin: nested case-control study. Cancer Epidemiol Biomarkers Prev. (2016) 25:721–4. doi: 10.1158/1055-9965.EPI-15-1290
103. Weissenborn SJ, Nindl I, Purdie K, Harwood C, Proby C, Breuer J, et al. Human papillomavirus-DNA loads in actinic keratoses exceed those in non-melanoma skin cancers. J Invest Dermatol. (2005) 125:93–7. doi: 10.1111/j.0022-202X.2005.23733.x
104. Viarisio D, Muller-Decker K, Accardi R, Robitaille A, Durst M, Beer K, et al. Beta HPV38 oncoproteins act with a hit-and-run mechanism in ultraviolet radiation-induced skin carcinogenesis in mice. PLoS Pathog. (2018) 14:e1006783. doi: 10.1371/journal.ppat.1006783
105. Arron ST, Ruby JG, Dybbro E, Ganem D, Derisi JL. Transcriptome sequencing demonstrates that human papillomavirus is not active in cutaneous squamous cell carcinoma. J Invest Dermatol. (2011) 131:1745–53. doi: 10.1038/jid.2011.91
106. Nafz J, Kohler A, Ohnesorge M, Nindl I, Stockfleth E, Rosl F. Persistence of Mastomys natalensis papillomavirus in multiple organs identifies novel targets for infection. J Gen Virol. (2007) 88:2670–8. doi: 10.1099/vir.0.82955-0
107. Hasche D, Stephan S, Braspenning-Wesch I, Mikulec J, Niebler M, Grone HJ, et al. The interplay of UV and cutaneous papillomavirus infection in skin cancer development. PLoS Pathog. (2017) 13:e1006723. doi: 10.1371/journal.ppat.1006723
108. Andersson K, Waterboer T, Kirnbauer R, Slupetzky K, Iftner T, de Villiers EM, et al. Seroreactivity to cutaneous human papillomaviruses among patients with nonmelanoma skin cancer or benign skin lesions. Cancer Epidemiol Biomarkers Prev. (2008) 17:189–95. doi: 10.1158/1055-9965.EPI-07-0405
109. Wong MCS, Vlantis AC, Liang M, Wong PY, Ho WCS, Boon SS, et al. Prevalence and epidemiologic profile of oral infection with alpha, beta, and gamma papillomaviruses in an asian Chinese population. J Infect Dis. (2018) 218:388–97. doi: 10.1093/infdis/jiy160
110. Nunes EM, Lopez RVM, Sudenga SL, Gheit T, Tommasino M, Baggio ML, et al. Concordance of Beta-papillomavirus across anogenital and oral anatomic sites of men: the HIM Study. Virology. (2017) 510:55–9. doi: 10.1016/j.virol.2017.07.006
111. Hampras SS, Rollison DE, Giuliano AR, McKay-Chopin S, Minoni L, Sereday K, et al. Prevalence and concordance of cutaneous beta human papillomavirus infection at mucosal and cutaneous sites. J Infect Dis. (2017) 216:92–6. doi: 10.1093/infdis/jix245
112. Smelov V, Muwonge R, Sokolova O, McKay-Chopin S, Eklund C, Komyakov B, et al. Beta and gamma human papillomaviruses in anal and genital sites among men: prevalence and determinants. Sci Rep. (2018) 8:8241. doi: 10.1038/s41598-018-26589-w
113. Pierce Campbell CM, Messina JL, Stoler MH, Jukic DM, Tommasino M, Gheit T, et al. Cutaneous human papillomavirus types detected on the surface of male external genital lesions: a case series within the HPV Infection in Men Study. J Clin Virol. (2013) 58:652–9. doi: 10.1016/j.jcv.2013.10.011
114. ElAmrani A, Gheit T, Benhessou M, McKay-Chopin S, Attaleb M, Sahraoui S, et al. Prevalence of mucosal and cutaneous human papillomavirus in Moroccan breast cancer. Papillomavirus Res. (2018) 5:150–5. doi: 10.1016/j.pvr.2018.04.003
115. Chen AA, Gheit T, Stellin M, Lupato V, Spinato G, Fuson R, et al. Oncogenic DNA viruses found in salivary gland tumors. Oral Oncol. (2017) 75:106–10. doi: 10.1016/j.oraloncology.2017.11.005
116. Agalliu I, Chen Z, Wang T, Hayes RB, Freedman ND, Gapstur SM, et al. Oral alpha, beta, and gamma HPV types and risk of incident esophageal cancer. Cancer Epidemiol Biomarkers Prev. (2018) 27:1168–75. doi: 10.1158/1055-9965.EPI-18-0287
117. Agalliu I, Gapstur S, Chen Z, Wang T, Anderson RL, Teras L, et al. Associations of oral alpha-, beta-, and gamma-human papillomavirus types with risk of incident head and neck cancer. JAMA Oncol. (2016) 2:599–606. doi: 10.1001/jamaoncol.2015.5504
118. Sabol I, Smahelova J, Klozar J, Mravak-Stipetic M, Gheit T, Tommasino M, et al. Beta-HPV types in patients with head and neck pathology and in healthy subjects. J Clin Virol. (2016) 82:159–65. doi: 10.1016/j.jcv.2016.07.019
119. Forslund O, Johansson H, Madsen KG, Kofoed K. The nasal mucosa contains a large spectrum of human papillomavirus types from the Betapapillomavirus and Gammapapillomavirus genera. J Infect Dis. (2013) 208:1335–41. doi: 10.1093/infdis/jit326
120. Cornet I, Bouvard V, Campo MS, Thomas M, Banks L, Gissmann L, et al. Comparative analysis of transforming properties of E6 and E7 from different beta human papillomavirus types. J Virol. (2012) 86:2366–70. doi: 10.1128/JVI.06579-11
121. Viarisio D, Muller-Decker K, Zanna P, Kloz U, Aengeneyndt B, Accardi R, et al. Novel ss-HPV49 transgenic mouse model of upper digestive tract cancer. Cancer Res. (2016) 76:4216–25. doi: 10.1158/0008-5472.CAN-16-0370
122. Grace M, Munger K. Proteomic analysis of the gamma human papillomavirus type 197 E6 and E7 associated cellular proteins. Virology. (2017) 500:71–81. doi: 10.1016/j.virol.2016.10.010
123. Arroyo Muhr LS, Hultin E, Bzhalava D, Eklund C, Lagheden C, Ekstrom J, et al. Human papillomavirus type 197 is commonly present in skin tumors. Int J Cancer. (2015) 136:2546–55. doi: 10.1002/ijc.29325
124. Kines RC, Thompson CD, Lowy DR, Schiller JT, Day PM. The initial steps leading to papillomavirus infection occur on the basement membrane prior to cell surface binding. Proc Natl Acad Sci USA. (2009) 106:20458–63. doi: 10.1073/pnas.0908502106
125. Schiller JT, Day PM, Kines RC. Current understanding of the mechanism of HPV infection. Gynecol Oncol. (2010) 118(Suppl. 1):S12–7. doi: 10.1016/j.ygyno.2010.04.004
126. Herfs M, Yamamoto Y, Laury A, Wang X, Nucci MR, McLaughlin-Drubin ME, et al. A discrete population of squamocolumnar junction cells implicated in the pathogenesis of cervical cancer. Proc Natl Acad Sci USA. (2012) 109:10516–21. doi: 10.1073/pnas.1202684109
127. Rautava J, Syrjanen S. Biology of human papillomavirus infections in head and neck carcinogenesis. Head Neck Pathol. (2012) 6 (Suppl. 1):S3–15. doi: 10.1007/s12105-012-0367-2
128. Broniarczyk J, Ring N, Massimi P, Giacca M, Banks L. HPV-16 virions can remain infectious for 2 weeks on senescent cells but require cell cycle re-activation to allow virus entry. Sci Rep. (2018) 8:811. doi: 10.1038/s41598-017-18809-6
129. Cerqueira C, Samperio Ventayol P, Vogeley C, Schelhaas M. Kallikrein-8 Proteolytically processes human papillomaviruses in the extracellular space to facilitate entry into host cells. J Virol. (2015) 89:7038–52. doi: 10.1128/JVI.00234-15
130. Culp TD, Budgeon LR, Marinkovich MP, Meneguzzi G, Christensen ND. Keratinocyte-secreted laminin 5 can function as a transient receptor for human papillomaviruses by binding virions and transferring them to adjacent cells. J Virol. (2006) 80:8940–50. doi: 10.1128/JVI.00724-06
131. Bronnimann MP, Calton CM, Chiquette SF, Li S, Lu M, Chapman JA, et al. Furin cleavage of l2 during papillomavirus infection: minimal dependence on cyclophilins. J Virol. (2016) 90:6224–34. doi: 10.1128/JVI.00038-16
132. Rubio I, Seitz H, Canali E, Sehr P, Bolchi A, Tommasino M, et al. The N-terminal region of the human papillomavirus L2 protein contains overlapping binding sites for neutralizing, cross-neutralizing and non-neutralizing antibodies. Virology. (2011) 409:348–59. doi: 10.1016/j.virol.2010.10.017
133. Schelhaas M, Ewers H, Rajamaki ML, Day PM, Schiller JT, Helenius A. Human papillomavirus type 16 entry: retrograde cell surface transport along actin-rich protrusions. PLoS Pathog. (2008) 4:e1000148. doi: 10.1371/journal.ppat.1000148
134. Schelhaas M, Shah B, Holzer M, Blattmann P, Kuhling L, Day PM, et al. Entry of human papillomavirus type 16 by actin-dependent, clathrin- and lipid raft-independent endocytosis. PLoS Pathog. (2012) 8:e1002657. doi: 10.1371/journal.ppat.1002657
135. Siddiqa A, Broniarczyk J, Banks L. Papillomaviruses and endocytic trafficking. Int J Mol Sci. (2018) 19:E2619 . doi: 10.3390/ijms19092619
136. Bienkowska-Haba M, Williams C, Kim SM, Garcea RL, Sapp M. Cyclophilins facilitate dissociation of the human papillomavirus type 16 capsid protein L1 from the L2/DNA complex following virus entry. J Virol. (2012) 86:9875–87. doi: 10.1128/JVI.00980-12
137. DiGiuseppe S, Bienkowska-Haba M, Guion LGM, Keiffer TR, Sapp M. Human papillomavirus major capsid protein L1 remains associated with the incoming viral genome throughout the entry process. J Virol. (2017) 91:e00537–17. doi: 10.1128/JVI.00537-17
138. Popa A, Zhang W, Harrison MS, Goodner K, Kazakov T, Goodwin EC, et al. Direct binding of retromer to human papillomavirus type 16 minor capsid protein L2 mediates endosome exit during viral infection. PLoS Pathog. (2015) 11:e1004699. doi: 10.1371/journal.ppat.1004699
139. McNally KE, Faulkner R, Steinberg F, Gallon M, Ghai R, Pim D, et al. Retriever is a multiprotein complex for retromer-independent endosomal cargo recycling. Nat Cell Biol. (2017) 19:1214–25. doi: 10.1038/ncb3610
140. Bergant Marusic M, Ozbun MA, Campos SK, Myers MP, Banks L. Human papillomavirus L2 facilitates viral escape from late endosomes via sorting nexin 17. Traffic. (2012) 13:455–67. doi: 10.1111/j.1600-0854.2011.01320.x
141. Bergant M, Peternel S, Pim D, Broniarczyk J, Banks L. Characterizing the spatio-temporal role of sorting nexin 17 in human papillomavirus trafficking. J Gen Virol. (2017) 98:715–25. doi: 10.1099/jgv.0.000734
142. Pim D, Broniarczyk J, Bergant M, Playford MP, Banks L. A Novel PDZ domain interaction mediates the binding between human papillomavirus 16 L2 and sorting nexin 27 and modulates virion trafficking. J Virol. (2015) 89:10145–55. doi: 10.1128/JVI.01499-15
143. DiGiuseppe S, Bienkowska-Haba M, Sapp M. Human papillomavirus entry: hiding in a bubble. J Virol. (2016) 90:8032–5. doi: 10.1128/JVI.01065-16
144. Pyeon D, Pearce SM, Lank SM, Ahlquist P, Lambert PF. Establishment of human papillomavirus infection requires cell cycle progression. PLoS Pathog. (2009) 5:e1000318. doi: 10.1371/journal.ppat.1000318
145. Aydin I, Weber S, Snijder B, Samperio Ventayol P, Kuhbacher A, Becker M, et al. Large scale RNAi reveals the requirement of nuclear envelope breakdown for nuclear import of human papillomaviruses. PLoS Pathog. (2014) 10:e1004162. doi: 10.1371/journal.ppat.1004162
146. Ascoli CA, Maul GG. Identification of a novel nuclear domain. J Cell Biol. (1991) 112:785–95. doi: 10.1083/jcb.112.5.785
147. Day PM, Baker CC, Lowy DR, Schiller JT. Establishment of papillomavirus infection is enhanced by promyelocytic leukemia protein (PML) expression. Proc Natl Acad Sci USA. (2004) 101:14252–7. doi: 10.1073/pnas.0404229101
148. Tavalai N, Adler M, Scherer M, Riedl Y, Stamminger T. Evidence for a dual antiviral role of the major nuclear domain 10 component Sp100 during the immediate-early and late phases of the human cytomegalovirus replication cycle. J Virol. (2011) 85:9447–58. doi: 10.1128/JVI.00870-11
149. Stepp WH, Meyers JM, McBride AA. Sp100 provides intrinsic immunity against human papillomavirus infection. mBio. (2013) 4:e00845–13. doi: 10.1128/mBio.00845-13
150. Florin L, Schafer F, Sotlar K, Streeck RE, Sapp M. Reorganization of nuclear domain 10 induced by papillomavirus capsid protein l2. Virology. (2002) 295:97–107. doi: 10.1006/viro.2002.1360
151. Kivipold P, Vosa L, Ustav M, Kurg R. DAXX modulates human papillomavirus early gene expression and genome replication in U2OS cells. Virol J. (2015) 12:104. doi: 10.1186/s12985-015-0335-z
152. Doorbar J, Quint W, Banks L, Bravo IG, Stoler M, Broker TR, et al. The biology and life-cycle of human papillomaviruses. Vaccine. (2012) 30 (Suppl. 5):F55–70. doi: 10.1016/j.vaccine.2012.06.083
153. McKinney CC, Kim MJ, Chen D, McBride AA. Brd4 activates early viral transcription upon human papillomavirus 18 infection of primary keratinocytes. mBio. (2016) 7:e01644–16. doi: 10.1128/mBio.01644-16
154. Sanders CM, Stenlund A. Recruitment and loading of the E1 initiator protein: an ATP-dependent process catalysed by a transcription factor. EMBO J. (1998) 17:7044–55. doi: 10.1093/emboj/17.23.7044
155. Gauson EJ, Donaldson MM, Dornan ES, Wang X, Bristol M, Bodily JM, et al. Evidence supporting a role for TopBP1 and Brd4 in the initiation but not continuation of human papillomavirus 16 E1/E2-mediated DNA replication. J Virol. (2015) 89:4980–91. doi: 10.1128/JVI.00335-15
156. Dreer M, Fertey J, van de Poel S, Straub E, Madlung J, Macek B, et al. Interaction of NCOR/SMRT repressor complexes with papillomavirus E8∧E2C proteins inhibits viral replication. PLoS Pathog. (2016) 12:e1005556. doi: 10.1371/journal.ppat.1005556
157. Dreer M, van de Poel S, Stubenrauch F. Control of viral replication and transcription by the papillomavirus E8∧E2 protein. Virus Res. (2017) 231:96–102. doi: 10.1016/j.virusres.2016.11.005
158. Bastien N, McBride AA. Interaction of the papillomavirus E2 protein with mitotic chromosomes. Virology. (2000) 270:124–34. doi: 10.1006/viro.2000.0265
159. Jang MK, Shen K, McBride AA. Papillomavirus genomes associate with BRD4 to replicate at fragile sites in the host genome. PLoS Pathog. (2014) 10:e1004117. doi: 10.1371/journal.ppat.1004117
160. McBride AA, Jang MK. Current understanding of the role of the Brd4 protein in the papillomavirus lifecycle. Viruses. (2013) 5:1374–94. doi: 10.3390/v5061374
161. Abroi A, Ilves I, Kivi S, Ustav M. Analysis of chromatin attachment and partitioning functions of bovine papillomavirus type 1 E2 protein. J Virol. (2004) 78:2100–13. doi: 10.1128/JVI.78.4.2100-2113.2004
162. Rogers A, Waltke M, Angeletti PC. Evolutionary variation of papillomavirus E2 protein and E2 binding sites. Virol J. (2011) 8:379. doi: 10.1186/1743-422X-8-379
163. McPhillips MG, Oliveira JG, Spindler JE, Mitra R, McBride AA. Brd4 is required for e2-mediated transcriptional activation but not genome partitioning of all papillomaviruses. J Virol. (2006) 80:9530–43. doi: 10.1128/JVI.01105-06
164. Parish JL, Bean AM, Park RB, Androphy EJ. ChlR1 is required for loading papillomavirus E2 onto mitotic chromosomes and viral genome maintenance. Mol Cell. (2006) 24:867–76. doi: 10.1016/j.molcel.2006.11.005
165. Parish JL, Rosa J, Wang X, Lahti JM, Doxsey SJ, Androphy EJ. The DNA helicase ChlR1 is required for sister chromatid cohesion in mammalian cells. J Cell Sci. (2006) 119:4857–65. doi: 10.1242/jcs.03262
166. Harris L, McFarlane-Majeed L, Campos-Leon K, Roberts S, Parish JL. The cellular DNA helicase ChlR1 regulates chromatin and nuclear matrix attachment of the human papillomavirus 16 e2 protein and high-copy-number viral genome establishment. J Virol. (2017) 91:16. doi: 10.1128/JVI.01853-16
167. Bentley P, Tan MJA, McBride AA, White EA, Howley PM. The SMC5/6 complex interacts with the papillomavirus E2 protein and influences maintenance of viral episomal DNA. J Virol. (2018) 92:e00356–18. doi: 10.1128/JVI.00356-18
168. Thomas JT, Hubert WG, Ruesch MN, Laimins LA. Human papillomavirus type 31 oncoproteins E6 and E7 are required for the maintenance of episomes during the viral life cycle in normal human keratinocytes. Proc Natl Acad Sci USA. (1999) 96:8449–54. doi: 10.1073/pnas.96.15.8449
169. Lorenz LD, Rivera Cardona J, Lambert PF. Inactivation of p53 rescues the maintenance of high risk HPV DNA genomes deficient in expression of E6. PLoS Pathog. (2013) 9:e1003717. doi: 10.1371/journal.ppat.1003717
170. Moody CA, Laimins LA. Human papillomavirus oncoproteins: pathways to transformation. Nat Rev Cancer. (2010) 10:550–60. doi: 10.1038/nrc2886
171. Maglennon GA, McIntosh P, Doorbar J. Persistence of viral DNA in the epithelial basal layer suggests a model for papillomavirus latency following immune regression. Virology. (2011) 414:153–63. doi: 10.1016/j.virol.2011.03.019
172. Hummel M, Hudson JB, Laimins LA. Differentiation-induced and constitutive transcription of human papillomavirus type 31b in cell lines containing viral episomes. J Virol. (1992) 66:6070–80.
173. Fehrmann F, Klumpp DJ, Laimins LA. Human papillomavirus type 31 E5 protein supports cell cycle progression and activates late viral functions upon epithelial differentiation. J Virol. (2003) 77:2819–31. doi: 10.1128/JVI.77.5.2819-2831.2003
174. Longworth MS, Laimins LA. Pathogenesis of human papillomaviruses in differentiating epithelia. Microbiol Mol Biol Rev. (2004) 68:362–72. doi: 10.1128/MMBR.68.2.362-372.2004
175. Khan J, Davy CE, McIntosh PB, Jackson DJ, Hinz S, Wang Q, et al. Role of calpain in the formation of human papillomavirus type 16 E1∧E4 amyloid fibers and reorganization of the keratin network. J Virol. (2011) 85:9984–97. doi: 10.1128/JVI.02158-10
176. Pentland I, Campos-Leon K, Cotic M, Davies KJ, Wood CD, Groves IJ, et al. Disruption of CTCF-YY1-dependent looping of the human papillomavirus genome activates differentiation-induced viral oncogene transcription. PLoS Biol. (2018) 16:e2005752. doi: 10.1371/journal.pbio.2005752
177. Butz K, Ristriani T, Hengstermann A, Denk C, Scheffner M, Hoppe-Seyler F. siRNA targeting of the viral E6 oncogene efficiently kills human papillomavirus-positive cancer cells. Oncogene. (2003) 22:5938–45. doi: 10.1038/sj.onc.1206894
178. Fujii T, Saito M, Iwasaki E, Ochiya T, Takei Y, Hayashi S, et al. Intratumor injection of small interfering RNA-targeting human papillomavirus 18 E6 and E7 successfully inhibits the growth of cervical cancer. Int J Oncol. (2006) 29:541–8. doi: 10.3892/ijo.29.3.541
179. Howie HL, Katzenellenbogen RA, Galloway DA. Papillomavirus E6 proteins. Virology. (2009) 384:324–34. doi: 10.1016/j.virol.2008.11.017
180. Scheffner M, Huibregtse JM, Vierstra RD, Howley PM. The HPV-16 E6 and E6-AP complex functions as a ubiquitin-protein ligase in the ubiquitination of p53. Cell. (1993) 75:495–505. doi: 10.1016/0092-8674(93)90384-3
181. Scheffner M, Werness BA, Huibregtse JM, Levine AJ, Howley PM. The E6 oncoprotein encoded by human papillomavirus types 16 and 18 promotes the degradation of p53. Cell. (1990) 63:1129–36. doi: 10.1016/0092-8674(90)90409-8
182. Huibregtse JM, Scheffner M, Howley PM. Localization of the E6-AP regions that direct human papillomavirus E6 binding, association with p53, and ubiquitination of associated proteins. Mol Cell Biol. (1993) 13:4918–27. doi: 10.1128/MCB.13.8.4918
183. Martinez-Zapien D, Ruiz FX, Poirson J, Mitschler A, Ramirez J, Forster A, et al. Structure of the E6/E6AP/p53 complex required for HPV-mediated degradation of p53. Nature. (2016) 529:541–5. doi: 10.1038/nature16481
184. Muschik D, Braspenning-Wesch I, Stockfleth E, Rosl F, Hofmann TG, Nindl I. Cutaneous HPV23 E6 prevents p53 phosphorylation through interaction with HIPK2. PLoS ONE. (2011) 6:e27655. doi: 10.1371/journal.pone.0027655
185. Pim D, Storey A, Thomas M, Massimi P, Banks L. Mutational analysis of HPV-18 E6 identifies domains required for p53 degradation in vitro, abolition of p53 transactivation in vivo and immortalisation of primary BMK cells. Oncogene. (1994) 9:1869–76.
186. Zimmermann H, Degenkolbe R, Bernard HU, O'Connor MJ. The human papillomavirus type 16 E6 oncoprotein can down-regulate p53 activity by targeting the transcriptional coactivator CBP/p300. J Virol. (1999) 73:6209–19.
187. Thomas MC, Chiang CM. E6 oncoprotein represses p53-dependent gene activation via inhibition of protein acetylation independently of inducing p53 degradation. Mol Cell. (2005) 17:251–64. doi: 10.1016/j.molcel.2004.12.016
188. Howie HL, Koop JI, Weese J, Robinson K, Wipf G, Kim L, et al. Beta-HPV 5 and 8 E6 promote p300 degradation by blocking AKT/p300 association. PLoS Pathog. (2011) 7:e1002211. doi: 10.1371/journal.ppat.1002211
189. Wallace NA, Robinson K, Howie HL, Galloway DA. HPV 5 and 8 E6 abrogate ATR activity resulting in increased persistence of UVB induced DNA damage. PLoS Pathog. (2012) 8:e1002807. doi: 10.1371/journal.ppat.1002807
190. Hufbauer M, Cooke J, van der Horst GT, Pfister H, Storey A, Akgul B. Human papillomavirus mediated inhibition of DNA damage sensing and repair drives skin carcinogenesis. Mol Cancer. (2015) 14:183. doi: 10.1186/s12943-015-0453-7
191. Marthaler AM, Podgorska M, Feld P, Fingerle A, Knerr-Rupp K, Grasser F, et al. Identification of C/EBPalpha as a novel target of the HPV8 E6 protein regulating miR-203 in human keratinocytes. PLoS Pathog. (2017) 13:e1006406. doi: 10.1371/journal.ppat.1006406
192. Muench P, Probst S, Schuetz J, Leiprecht N, Busch M, Wesselborg S, et al. Cutaneous papillomavirus E6 proteins must interact with p300 and block p53-mediated apoptosis for cellular immortalization and tumorigenesis. Cancer Res. (2010) 70:6913–24. doi: 10.1158/0008-5472.CAN-10-1307
193. Thomas M, Banks L. Human papillomavirus (HPV) E6 interactions with Bak are conserved amongst E6 proteins from high and low risk HPV types. J Gen Virol. (1999) 80:1513–7. doi: 10.1099/0022-1317-80-6-1513
194. Borbely AA, Murvai M, Konya J, Beck Z, Gergely L, Li F, et al. Effects of human papillomavirus type 16 oncoproteins on survivin gene expression. J Gen Virol. (2006) 87:287–94. doi: 10.1099/vir.0.81067-0
195. Filippova M, Song H, Connolly JL, Dermody TS, Duerksen-Hughes PJ. The human papillomavirus 16 E6 protein binds to tumor necrosis factor (TNF) R1 and protects cells from TNF-induced apoptosis. J Biol Chem. (2002) 277:21730–9. doi: 10.1074/jbc.M200113200
196. Filippova M, Parkhurst L, Duerksen-Hughes PJ. The human papillomavirus 16 E6 protein binds to Fas-associated death domain and protects cells from Fas-triggered apoptosis. J Biol Chem. (2004) 279:25729–44. doi: 10.1074/jbc.M401172200
197. Filippova M, Johnson MM, Bautista M, Filippov V, Fodor N, Tungteakkhun SS, et al. The large and small isoforms of human papillomavirus type 16 E6 bind to and differentially affect procaspase 8 stability and activity. J Virol. (2007) 81:4116–29. doi: 10.1128/JVI.01924-06
198. Bedard KM, Underbrink MP, Howie HL, Galloway DA. The E6 oncoproteins from human betapapillomaviruses differentially activate telomerase through an E6AP-dependent mechanism and prolong the lifespan of primary keratinocytes. J Virol. (2008) 82:3894–902. doi: 10.1128/JVI.01818-07
199. Simmonds M, Storey A. Identification of the regions of the HPV 5 E6 protein involved in Bak degradation and inhibition of apoptosis. Int J Cancer. (2008) 123:2260–6. doi: 10.1002/ijc.23815
200. Holloway A, Simmonds M, Azad A, Fox JL, Storey A. Resistance to UV-induced apoptosis by beta-HPV5 E6 involves targeting of activated BAK for proteolysis by recruitment of the HERC1 ubiquitin ligase. Int J Cancer. (2015) 136:2831–43. doi: 10.1002/ijc.29350
201. Leverrier S, Bergamaschi D, Ghali L, Ola A, Warnes G, Akgul B, et al. Role of HPV E6 proteins in preventing UVB-induced release of pro-apoptotic factors from the mitochondria. Apoptosis. (2007) 12:549–60. doi: 10.1007/s10495-006-0004-1
202. Wallace NA, Robinson K, Howie HL, Galloway DA. beta-HPV 5 and 8 E6 disrupt homology dependent double strand break repair by attenuating BRCA1 and BRCA2 expression and foci formation. PLoS Pathog. (2015) 11:e1004687. doi: 10.1371/journal.ppat.1004687
203. Schwarz T, Schwarz A. DNA repair and cytokine responses. J Invest Dermatol Symp Proc. (2009) 14:63–6. doi: 10.1038/jidsymp.2009.3
204. Javier RT. Cell polarity proteins: common targets for tumorigenic human viruses. Oncogene. (2008) 27:7031–46. doi: 10.1038/onc.2008.352
205. Songyang Z, Fanning AS, Fu C, Xu J, Marfatia SM, Chishti AH, et al. Recognition of unique carboxyl-terminal motifs by distinct PDZ domains. Science. (1997) 275:73–7. doi: 10.1126/science.275.5296.73
206. Kuhne C, Gardiol D, Guarnaccia C, Amenitsch H, Banks L. Differential regulation of human papillomavirus E6 by protein kinase A: conditional degradation of human discs large protein by oncogenic E6. Oncogene. (2000) 19:5884–91. doi: 10.1038/sj.onc.1203988
207. Zhang Y, Dasgupta J, Ma RZ, Banks L, Thomas M, Chen XS. Structures of a human papillomavirus (HPV) E6 polypeptide bound to MAGUK proteins: mechanisms of targeting tumor suppressors by a high-risk HPV oncoprotein. J Virol. (2007) 81:3618–26. doi: 10.1128/JVI.02044-06
208. Boon SS, Banks L. High-risk human papillomavirus E6 oncoproteins interact with 14-3-3zeta in a PDZ binding motif-dependent manner. J Virol. (2013) 87:1586–95. doi: 10.1128/JVI.02074-12
209. Nakagawa S, Huibregtse JM. Human scribble (Vartul) is targeted for ubiquitin-mediated degradation by the high-risk papillomavirus E6 proteins and the E6AP ubiquitin-protein ligase. Mol Cell Biol. (2000) 20:8244–53. doi: 10.1128/MCB.20.21.8244-8253.2000
210. Gardiol D, Kuhne C, Glaunsinger B, Lee SS, Javier R, Banks L. Oncogenic human papillomavirus E6 proteins target the discs large tumour suppressor for proteasome-mediated degradation. Oncogene. (1999) 18:5487–96. doi: 10.1038/sj.onc.1202920
211. Kranjec C, Tomaic V, Massimi P, Nicolaides L, Doorbar J, Banks L. The high-risk HPV E6 target scribble (hScrib) is required for HPV E6 expression in cervical tumour-derived cell lines. Papillomavirus Res. (2016) 2:70–7. doi: 10.1016/j.pvr.2016.04.001
212. Kapil S, Sharma BK, Patil M, Elattar S, Yuan J, Hou SX, et al. The cell polarity protein Scrib functions as a tumor suppressor in liver cancer. Oncotarget. (2017) 8:26515–31. doi: 10.18632/oncotarget.15713
213. Kranjec C, Banks L. A systematic analysis of human papillomavirus (HPV) E6 PDZ substrates identifies MAGI-1 as a major target of HPV type 16 (HPV-16) and HPV-18 whose loss accompanies disruption of tight junctions. J Virol. (2011) 85:1757–64. doi: 10.1128/JVI.01756-10
214. Dobrosotskaya IY, James GL. MAGI-1 interacts with beta-catenin and is associated with cell-cell adhesion structures. Biochem Biophys Res Commun. (2000) 270:903–9. doi: 10.1006/bbrc.2000.2471
215. Kotelevets L, van Hengel J, Bruyneel E, Mareel M, van Roy F, Chastre E. Implication of the MAGI-1b/PTEN signalosome in stabilization of adherens junctions and suppression of invasiveness. FASEB J. (2005) 19:115–7. doi: 10.1096/fj.04-1942fje
216. Jeong KW, Kim HZ, Kim S, Kim YS, Choe J. Human papillomavirus type 16 E6 protein interacts with cystic fibrosis transmembrane regulator-associated ligand and promotes E6-associated protein-mediated ubiquitination and proteasomal degradation. Oncogene. (2007) 26:487–99. doi: 10.1038/sj.onc.1209837
217. Lee SS, Glaunsinger B, Mantovani F, Banks L, Javier RT. Multi-PDZ domain protein MUPP1 is a cellular target for both adenovirus E4-ORF1 and high-risk papillomavirus type 18 E6 oncoproteins. J Virol. (2000) 74:9680–93. doi: 10.1128/JVI.74.20.9680-9693.2000
218. Handa K, Yugawa T, Narisawa-Saito M, Ohno S, Fujita M, Kiyono T. E6AP-dependent degradation of DLG4/PSD95 by high-risk human papillomavirus type 18 E6 protein. J Virol. (2007) 81:1379–89. doi: 10.1128/JVI.01712-06
219. Favre-Bonvin A, Reynaud C, Kretz-Remy C, Jalinot P. Human papillomavirus type 18 E6 protein binds the cellular PDZ protein TIP-2/GIPC, which is involved in transforming growth factor beta signaling and triggers its degradation by the proteasome. J Virol. (2005) 79:4229–37. doi: 10.1128/JVI.79.7.4229-4237.2005
220. Accardi R, Rubino R, Scalise M, Gheit T, Shahzad N, Thomas M, et al. E6 and E7 from human papillomavirus type 16 cooperate to target the PDZ protein Na/H exchange regulatory factor 1. J Virol. (2011) 85:8208–16. doi: 10.1128/JVI.00114-11
221. Lee C, Laimins LA. Role of the PDZ domain-binding motif of the oncoprotein E6 in the pathogenesis of human papillomavirus type 31. J Virol. (2004) 78:12366–77. doi: 10.1128/JVI.78.22.12366-12377.2004
222. Marx B, Miller-Lazic D, Doorbar J, Majewski S, Hofmann K, Hufbauer M, et al. HPV8-E6 Interferes with Syntenin-2 expression through deregulation of differentiation, methylation and phosphatidylinositide-kinase dependent mechanisms. Front Microbiol. (2017) 8:1724. doi: 10.3389/fmicb.2017.01724
223. Thomas M, Myers MP, Massimi P, Guarnaccia C, Banks L. Analysis of multiple HPV E6 PDZ interactions defines type-specific PDZ fingerprints that predict oncogenic potential. PLoS Pathog. (2016) 12:e1005766. doi: 10.1371/journal.ppat.1005766
224. Branca M, Giorgi C, Ciotti M, Santini D, Di Bonito L, Costa S, et al. Upregulation of telomerase (hTERT) is related to the grade of cervical intraepithelial neoplasia, but is not an independent predictor of high-risk human papillomavirus, virus persistence, or disease outcome in cervical cancer. Diag Cytopathol. (2006) 34:739–48. doi: 10.1002/dc.20554
225. Van Doorslaer K, Burk RD. Association between hTERT activation by HPV E6 proteins and oncogenic risk. Virology. (2012) 433:216–9. doi: 10.1016/j.virol.2012.08.006
226. Cohen SB, Graham ME, Lovrecz GO, Bache N, Robinson PJ, Reddel RR. Protein composition of catalytically active human telomerase from immortal cells. Science. (2007) 315:1850–3. doi: 10.1126/science.1138596
227. Wolyniec K, Shortt J, de Stanchina E, Levav-Cohen Y, Alsheich-Bartok O, Louria-Hayon I, et al. E6AP ubiquitin ligase regulates PML-induced senescence in Myc-driven lymphomagenesis. Blood. (2012) 120:822–32. doi: 10.1182/blood-2011-10-387647
228. James MA, Lee JH, Klingelhutz AJ. HPV16-E6 associated hTERT promoter acetylation is E6AP dependent, increased in later passage cells and enhanced by loss of p300. Int J Cancer. (2006) 119:1878–85. doi: 10.1002/ijc.22064
229. Lebel R, McDuff FO, Lavigne P, Grandbois M. Direct visualization of the binding of c-Myc/Max heterodimeric b-HLH-LZ to E-box sequences on the hTERT promoter. Biochemistry. (2007) 46:10279–86. doi: 10.1021/bi700076m
230. Xu M, Luo W, Elzi DJ, Grandori C, Galloway DA. NFX1 interacts with mSin3A/histone deacetylase to repress hTERT transcription in keratinocytes. Mol Cell Biol. (2008) 28:4819–28. doi: 10.1128/MCB.01969-07
231. Katzenellenbogen RA, Egelkrout EM, Vliet-Gregg P, Gewin LC, Gafken PR, Galloway DA. NFX1-123 and poly(A) binding proteins synergistically augment activation of telomerase in human papillomavirus type 16 E6-expressing cells. J Virol. (2007) 81:3786–96. doi: 10.1128/JVI.02007-06
232. Mortensen F, Schneider D, Barbic T, Sladewska-Marquardt A, Kuhnle S, Marx A, et al. Role of ubiquitin and the HPV E6 oncoprotein in E6AP-mediated ubiquitination. Proc Natl Acad Sci USA. (2015) 112:9872–7. doi: 10.1073/pnas.1505923112
233. Brimer N, Drews CM, Vande Pol SB. Association of papillomavirus E6 proteins with either MAML1 or E6AP clusters E6 proteins by structure, function, and evolutionary relatedness. PLoS Pathog. (2017) 13:e1006781. doi: 10.1371/journal.ppat.1006781
234. Brimer N, Lyons C, Wallberg AE, Vande Pol SB. Cutaneous papillomavirus E6 oncoproteins associate with MAML1 to repress transactivation and NOTCH signaling. Oncogene. (2012) 31:4639–46. doi: 10.1038/onc.2011.589
235. Meyers JM, Spangle JM, Munger K. The human papillomavirus type 8 E6 protein interferes with NOTCH activation during keratinocyte differentiation. J Virol. (2013) 87:4762–7. doi: 10.1128/JVI.02527-12
236. White EA, Howley PM. Proteomic approaches to the study of papillomavirus-host interactions. Virology. (2013) 435:57–69. doi: 10.1016/j.virol.2012.09.046
237. Knapp AA, McManus PM, Bockstall K, Moroianu J. Identification of the nuclear localization and export signals of high risk HPV16 E7 oncoprotein. Virology. (2009) 383:60–8. doi: 10.1016/j.virol.2008.09.037
238. Huh KW, DeMasi J, Ogawa H, Nakatani Y, Howley PM, Munger K. Association of the human papillomavirus type 16 E7 oncoprotein with the 600-kDa retinoblastoma protein-associated factor, p600. Proc Natl Acad Sci USA. (2005) 102:11492–7. doi: 10.1073/pnas.0505337102
239. Nakatani Y, Konishi H, Vassilev A, Kurooka H, Ishiguro K, Sawada J, et al. p600, a unique protein required for membrane morphogenesis and cell survival. Proc Natl Acad Sci USA. (2005) 102:15093–8. doi: 10.1073/pnas.0507458102
240. DeMasi J, Chao MC, Kumar AS, Howley PM. Bovine papillomavirus E7 oncoprotein inhibits anoikis. J Virol. (2007) 81:9419–25. doi: 10.1128/JVI.00422-07
241. Szalmas A, Tomaic V, Basukala O, Massimi P, Mittal S, Konya J, et al. The PTPN14 tumor suppressor is a degradation target of human papillomavirus E7. J Virol. (2017) 91:e00057–17. doi: 10.1128/JVI.00057-17
242. White EA, Munger K, Howley PM. High-risk human papillomavirus E7 proteins target PTPN14 for degradation. mBio. (2016) 7:16. doi: 10.1128/mBio.01530-16
243. Avvakumov N, Torchia J, Mymryk JS. Interaction of the HPV E7 proteins with the pCAF acetyltransferase. Oncogene. (2003) 22:3833–41. doi: 10.1038/sj.onc.1206562
244. Jansma AL, Martinez-Yamout MA, Liao R, Sun P, Dyson HJ, Wright PE. The high-risk HPV16 E7 oncoprotein mediates interaction between the transcriptional coactivator CBP and the retinoblastoma protein pRb. J Mol Biol. (2014) 426:4030–48. doi: 10.1016/j.jmb.2014.10.021
245. Massimi P, Pim D, Storey A, Banks L. HPV-16 E7 and adenovirus E1a complex formation with TATA box binding protein is enhanced by casein kinase II phosphorylation. Oncogene. (1996) 12:2325–30.
246. Barbosa MS, Edmonds C, Fisher C, Schiller JT, Lowy DR, Vousden KH. The region of the HPV E7 oncoprotein homologous to adenovirus E1a and Sv40 large T antigen contains separate domains for Rb binding and casein kinase II phosphorylation. EMBO J. (1990) 9:153–60. doi: 10.1002/j.1460-2075.1990.tb08091.x
247. Zine El Abidine A, Tomaic V, Bel Haj Rhouma R, Massimi P, Guizani I, Boubaker S, et al. A naturally occurring variant of HPV-16 E7 exerts increased transforming activity through acquisition of an additional phospho-acceptor site. Virology. (2017) 500:218–25. doi: 10.1016/j.virol.2016.10.023
248. Munger K, Werness BA, Dyson N, Phelps WC, Harlow E, Howley PM. Complex formation of human papillomavirus E7 proteins with the retinoblastoma tumor suppressor gene product. EMBO J. (1989) 8:4099–105. doi: 10.1002/j.1460-2075.1989.tb08594.x
249. Todorovic B, Hung K, Massimi P, Avvakumov N, Dick FA, Shaw GS, et al. Conserved region 3 of human papillomavirus 16 E7 contributes to deregulation of the retinoblastoma tumor suppressor. J Virol. (2012) 86:13313–23. doi: 10.1128/JVI.01637-12
250. Huh K, Zhou X, Hayakawa H, Cho JY, Libermann TA, Jin J, et al. Human papillomavirus type 16 E7 oncoprotein associates with the cullin 2 ubiquitin ligase complex, which contributes to degradation of the retinoblastoma tumor suppressor. J Virol. (2007) 81:9737–47. doi: 10.1128/JVI.00881-07
251. Chellappan S, Kraus VB, Kroger B, Munger K, Howley PM, Phelps WC, et al. Adenovirus E1A, simian virus 40 tumor antigen, and human papillomavirus E7 protein share the capacity to disrupt the interaction between transcription factor E2F and the retinoblastoma gene product. Proc Natl Acad Sci USA. (1992) 89:4549–53. doi: 10.1073/pnas.89.10.4549
252. Roman A, Munger K. The papillomavirus E7 proteins. Virology. (2013) 445:138–68. doi: 10.1016/j.virol.2013.04.013
253. Banks L, Edmonds C, Vousden KH. Ability of the HPV16 E7 protein to bind RB and induce DNA synthesis is not sufficient for efficient transforming activity in NIH3T3 cells. Oncogene. (1990) 5:1383–9.
254. Barrow-Laing L, Chen W, Roman A. Low- and high-risk human papillomavirus E7 proteins regulate p130 differently. Virology. (2010) 400:233–9. doi: 10.1016/j.virol.2010.01.034
255. McLaughlin-Drubin ME, Huh KW, Munger K. Human papillomavirus type 16 E7 oncoprotein associates with E2F6. J Virol. (2008) 82:8695–705. doi: 10.1128/JVI.00579-08
256. Akgul B, Garcia-Escudero R, Ghali L, Pfister HJ, Fuchs PG, Navsaria H, et al. The E7 protein of cutaneous human papillomavirus type 8 causes invasion of human keratinocytes into the dermis in organotypic cultures of skin. Cancer Res. (2005) 65:2216–23. doi: 10.1158/0008-5472.CAN-04-1952
257. Heuser S, Hufbauer M, Steiger J, Marshall J, Sterner-Kock A, Mauch C, et al. The fibronectin/alpha3beta1 integrin axis serves as molecular basis for keratinocyte invasion induced by betaHPV. Oncogene. (2016) 35:4529–39. doi: 10.1038/onc.2015.512
258. Hufbauer M, Biddle A, Borgogna C, Gariglio M, Doorbar J, Storey A, et al. Expression of betapapillomavirus oncogenes increases the number of keratinocytes with stem cell-like properties. J Virol. (2013) 87:12158–65. doi: 10.1128/JVI.01510-13
259. Lanfredini S, Olivero C, Borgogna C, Calati F, Powell K, Davies KJ, et al. HPV8 field cancerization in a transgenic mouse model is due to Lrig1+ keratinocyte stem cell expansion. J Invest Dermatol. (2017) 137:2208–16. doi: 10.1016/j.jid.2017.04.039
260. Lena AM, Shalom-Feuerstein R, Rivetti di Val Cervo P, Aberdam D, Knight RA, Melino G, et al. miR-203 represses 'stemness' by repressing DeltaNp63. Cell Death Different. (2008) 15:1187–95. doi: 10.1038/cdd.2008.69
261. Yi R, Poy MN, Stoffel M, Fuchs E. A skin microRNA promotes differentiation by repressing 'stemness'. Nature. (2008) 452:225–9. doi: 10.1038/nature06642
262. Shin MK, Balsitis S, Brake T, Lambert PF. Human papillomavirus E7 oncoprotein overrides the tumor suppressor activity of p21Cip1 in cervical carcinogenesis. Cancer Res. (2009) 69:5656–63. doi: 10.1158/0008-5472.CAN-08-3711
263. Funk JO, Waga S, Harry JB, Espling E, Stillman B, Galloway DA. Inhibition of CDK activity and PCNA-dependent DNA replication by p21 is blocked by interaction with the HPV-16 E7 oncoprotein. Genes Dev. (1997) 11:2090–100. doi: 10.1101/gad.11.16.2090
264. Jones DL, Alani RM, Munger K. The human papillomavirus E7 oncoprotein can uncouple cellular differentiation and proliferation in human keratinocytes by abrogating p21Cip1-mediated inhibition of cdk2. Genes Dev. (1997) 11:2101–11. doi: 10.1101/gad.11.16.2101
265. Akgul B, Ghali L, Davies D, Pfister H, Leigh IM, Storey A. HPV8 early genes modulate differentiation and cell cycle of primary human adult keratinocytes. Exp Dermatol. (2007) 16:590–9. doi: 10.1111/j.1600-0625.2007.00569.x
266. Schmitt A, Harry JB, Rapp B, Wettstein FO, Iftner T. Comparison of the properties of the E6 and E7 genes of low- and high-risk cutaneous papillomaviruses reveals strongly transforming and high Rb-binding activity for the E7 protein of the low-risk human papillomavirus type 1. J Virol. (1994) 68:7051–9.
267. Accardi R, Dong W, Smet A, Cui R, Hautefeuille A, Gabet AS, et al. Skin human papillomavirus type 38 alters p53 functions by accumulation of deltaNp73. EMBO Rep. (2006) 7:334–40. doi: 10.1038/sj.embor.7400615
268. Accardi R, Scalise M, Gheit T, Hussain I, Yue J, Carreira C, et al. IkappaB kinase beta promotes cell survival by antagonizing p53 functions through DeltaNp73alpha phosphorylation and stabilization. Mol Cell Biol. (2011) 31:2210–26. doi: 10.1128/MCB.00964-10
269. Wetherill LF, Holmes KK, Verow M, Muller M, Howell G, Harris M, et al. High-risk human papillomavirus E5 oncoprotein displays channel-forming activity sensitive to small-molecule inhibitors. J Virol. (2012) 86:5341–51. doi: 10.1128/JVI.06243-11
270. Oetke C, Auvinen E, Pawlita M, Alonso A. Human papillomavirus type 16 E5 protein localizes to the Golgi apparatus but does not grossly affect cellular glycosylation. Arch Virol. (2000) 145:2183–91. doi: 10.1007/s007050070048
271. Hu L, Ceresa BP. Characterization of the plasma membrane localization and orientation of HPV16 E5 for cell-cell fusion. Virology. (2009) 393:135–43. doi: 10.1016/j.virol.2009.07.034
272. Wetherill LF, Wasson CW, Swinscoe G, Kealy D, Foster R, Griffin S, et al. Alkyl-imino sugars inhibit the pro-oncogenic ion channel function of human papillomavirus (HPV) E5. Antiviral Res. (2018) 158:113–21. doi: 10.1016/j.antiviral.2018.08.005
273. Stoppler MC, Straight SW, Tsao G, Schlegel R, McCance DJ. The E5 gene of HPV-16 enhances keratinocyte immortalization by full-length DNA. Virology. (1996) 223:251–4. doi: 10.1006/viro.1996.0475
274. Barbaresi S, Cortese MS, Quinn J, Ashrafi GH, Graham SV, Campo MS. Effects of human papillomavirus type 16 E5 deletion mutants on epithelial morphology: functional characterization of each transmembrane domain. J Gen Virol. (2010) 91:521–30. doi: 10.1099/vir.0.016295-0
275. Genther Williams SM, Disbrow GL, Schlegel R, Lee D, Threadgill DW, Lambert PF. Requirement of epidermal growth factor receptor for hyperplasia induced by E5, a high-risk human papillomavirus oncogene. Cancer Res. (2005) 65:6534–42. doi: 10.1158/0008-5472.CAN-05-0083
276. Kim SH, Oh JM, No JH, Bang YJ, Juhnn YS, Song YS. Involvement of NF-kappaB and AP-1 in COX-2 upregulation by human papillomavirus 16 E5 oncoprotein. Carcinogenesis. (2009) 30:753–7. doi: 10.1093/carcin/bgp066
277. Scott ML, Coleman DT, Kelly KC, Carroll JL, Woodby B, Songock WK, et al. Human papillomavirus type 16 E5-mediated upregulation of Met in human keratinocytes. Virology. (2018) 519:1–11. doi: 10.1016/j.virol.2018.03.021
278. Pedroza-Saavedra A, Lam EW, Esquivel-Guadarrama F, Gutierrez-Xicotencatl L. The human papillomavirus type 16 E5 oncoprotein synergizes with EGF-receptor signaling to enhance cell cycle progression and the down-regulation of p27(Kip1). Virology. (2010) 400:44–52. doi: 10.1016/j.virol.2010.01.009
279. Liao S, Deng D, Hu X, Wang W, Li L, Li W, et al. HPV16/18 E5, a promising candidate for cervical cancer vaccines, affects SCPs, cell proliferation and cell cycle, and forms a potential network with E6 and E7. Int J Mol Med. (2013) 31:120–8. doi: 10.3892/ijmm.2012.1168
280. Kabsch K, Mossadegh N, Kohl A, Komposch G, Schenkel J, Alonso A, et al. The HPV-16 E5 protein inhibits TRAIL- and FasL-mediated apoptosis in human keratinocyte raft cultures. Intervirology. (2004) 47:48–56. doi: 10.1159/000076642
281. Liu C, Lin J, Li L, Zhang Y, Chen W, Cao Z, et al. HPV16 early gene E5 specifically reduces miRNA-196a in cervical cancer cells. Sci Rep. (2015) 5:7653. doi: 10.1038/srep07653
282. Oh JM, Kim SH, Cho EA, Song YS, Kim WH, Juhnn YS. Human papillomavirus type 16 E5 protein inhibits hydrogen-peroxide-induced apoptosis by stimulating ubiquitin-proteasome-mediated degradation of Bax in human cervical cancer cells. Carcinogenesis. (2010) 31:402–10. doi: 10.1093/carcin/bgp318
283. Zhang B, Spandau DF, Roman A. E5 protein of human papillomavirus type 16 protects human foreskin keratinocytes from UV B-irradiation-induced apoptosis. J Virol. (2002) 76:220–31. doi: 10.1128/JVI.76.1.220-231.2002
284. Valle GF, Banks L. The human papillomavirus (HPV)-6 and HPV-16 E5 proteins co-operate with HPV-16 E7 in the transformation of primary rodent cells. J Gen Virol. (1995) 76:1239–45. doi: 10.1099/0022-1317-76-5-1239
285. Lazarczyk M, Pons C, Mendoza JA, Cassonnet P, Jacob Y, Favre M. Regulation of cellular zinc balance as a potential mechanism of EVER-mediated protection against pathogenesis by cutaneous oncogenic human papillomaviruses. J Exp Med. (2008) 205:35–42. doi: 10.1084/jem.20071311
286. Kalinska-Bienias A, Kowalewski C, Majewski S. The EVER genes - the genetic etiology of carcinogenesis in epidermodysplasia verruciformis and a possible role in non-epidermodysplasia verruciformis patients. Postepy Dermatol Alergol. (2016) 33:75–80. doi: 10.5114/ada.2016.59145
287. de Jong SJ, Crequer A, Matos I, Hum D, Gunasekharan V, Lorenzo L, et al. The human CIB1-EVER1-EVER2 complex governs keratinocyte-intrinsic immunity to beta-papillomaviruses. J Exp Med. (2018) 215:2289–310. doi: 10.1084/jem.20170308
288. Chaturvedi AK, Engels EA, Anderson WF, Gillison ML. Incidence trends for human papillomavirus-related and -unrelated oral squamous cell carcinomas in the United States. J Clin Oncol. (2008) 26:612–9. doi: 10.1200/JCO.2007.14.1713
289. Nasman A, Attner P, Hammarstedt L, Du J, Eriksson M, Giraud G, et al. Incidence of human papillomavirus (HPV) positive tonsillar carcinoma in Stockholm, Sweden: an epidemic of viral-induced carcinoma? Int J Cancer. (2009) 125:362–6. doi: 10.1002/ijc.24339
290. Hong AM, Grulich AE, Jones D, Lee CS, Garland SM, Dobbins TA, et al. Squamous cell carcinoma of the oropharynx in Australian males induced by human papillomavirus vaccine targets. Vaccine. (2010) 28:3269–72. doi: 10.1016/j.vaccine.2010.02.098
291. Kreimer AR, Johansson M, Waterboer T, Kaaks R, Chang-Claude J, Drogen D, et al. Evaluation of human papillomavirus antibodies and risk of subsequent head and neck cancer. J Clin Oncol. (2013) 31:2708–15. doi: 10.1200/JCO.2012.47.2738
292. Kreimer AR, Johansson M, Yanik EL, Katki HA, Check DP, Lang Kuhs KA, et al. Kinetics of the human papillomavirus type 16 E6 antibody response prior to oropharyngeal cancer. J Natl Cancer Inst. (2017) 109:djx005 . doi: 10.1093/jnci/djx005
293. Harper DM, DeMars LR. HPV vaccines - a review of the first decade. Gynecol Oncol. (2017) 146:196–204. doi: 10.1016/j.ygyno.2017.04.004
294. Guo T, Eisele DW, Fakhry C. The potential impact of prophylactic human papillomavirus vaccination on oropharyngeal cancer. Cancer. (2016) 122:2313–23. doi: 10.1002/cncr.29992
295. Vici P, Pizzuti L, Mariani L, Zampa G, Santini D, Di Lauro L, et al. Targeting immune response with therapeutic vaccines in premalignant lesions and cervical cancer: hope or reality from clinical studies. Exp Rev Vaccines. (2016) 15:1327–36. doi: 10.1080/14760584.2016.1176533
296. Yang A, Farmer E, Wu TC, Hung CF. Perspectives for therapeutic HPV vaccine development. J Biomed Sci. (2016) 23:75. doi: 10.1186/s12929-016-0293-9
297. Kim HJ, Kim HJ. Current status and future prospects for human papillomavirus vaccines. Arch Pharmacal Res. (2017) 40:1050–63. doi: 10.1007/s12272-017-0952-8
298. Hefferon K. Plant virus expression vectors: a powerhouse for global health. Biomedicines. (2017) 5:E44 . doi: 10.3390/biomedicines5030044
299. Chabeda A, Yanez RJR, Lamprecht R, Meyers AE, Rybicki EP, Hitzeroth II. Therapeutic vaccines for high-risk HPV-associated diseases. Papillomavirus Res. (2018) 5:46–58. doi: 10.1016/j.pvr.2017.12.006
300. Gray P, Luostarinen T, Vanska S, Eriksson T, Lagheden C, Man I, et al. Occurrence of human papillomavirus (HPV) type replacement by sexual risk-taking behaviour group: Post-hoc analysis of a community randomized clinical trial up to 9 years after vaccination (IV). Int J Cancer. (2019). doi: 10.1002/ijc.32189. [Epub ahead of print].
301. Gray P, Palmroth J, Luostarinen T, Apter D, Dubin G, Garnett G, et al. Evaluation of HPV type-replacement in unvaccinated and vaccinated adolescent females-Post-hoc analysis of a community-randomized clinical trial (II). Int J Cancer. (2018) 142:2491–500. doi: 10.1002/ijc.31281
302. Klein C, Gonzalez D, Samwel K, Kahesa C, Mwaiselage J, Aluthge N, et al. Relationship between the cervical microbiome, HIV status, and precancerous lesions. mBio. (2019) 10 10:e02785–18. doi: 10.1128/mBio.02785-18
303. Bouwes Bavinck JN, Feltkamp MCW, Green AC, Fiocco M, Euvrard S, Harwood CA, et al. Human papillomavirus and posttransplantation cutaneous squamous cell carcinoma: a multicenter, prospective cohort study. Am J Transl. (2018) 18:1220–30. doi: 10.1111/ajt.14537
304. Huber B, Schellenbacher C, Shafti-Keramat S, Jindra C, Christensen N, Kirnbauer R. Chimeric L2-based virus-like particle (VLP) vaccines targeting cutaneous human papillomaviruses (HPV). PLoS ONE. (2017) 12:e0169533. doi: 10.1371/journal.pone.0169533
305. Vinzon SE, Braspenning-Wesch I, Muller M, Geissler EK, Nindl I, Grone HJ, et al. Protective vaccination against papillomavirus-induced skin tumors under immunocompetent and immunosuppressive conditions: a preclinical study using a natural outbred animal model. PLoS Pathog. (2014) 10:e1003924. doi: 10.1371/journal.ppat.1003924
306. Ostrbenk A, Kocjan BJ, Hosnjak L, Li J, Deng Q, Sterbenc A, et al. Identification of a Novel human papillomavirus, type HPV199, isolated from a nasopharynx and anal canal, and complete genomic characterization of papillomavirus species gamma-12. PLoS ONE. (2015) 10:e0138628. doi: 10.1371/journal.pone.0138628
Keywords: Papillomavirus, infection and cancer, transformation, anogenital cancer, skin cancer
Citation: Gheit T (2019) Mucosal and Cutaneous Human Papillomavirus Infections and Cancer Biology. Front. Oncol. 9:355. doi: 10.3389/fonc.2019.00355
Received: 10 January 2019; Accepted: 17 April 2019;
Published: 08 May 2019.
Edited by:
Mauro G. Tognon, University of Ferrara, ItalyReviewed by:
Baki Akgül, Universität zu Köln, GermanyLawrence Banks, International Centre for Genetic Engineering and Biotechnology, Italy
Copyright © 2019 Gheit. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Tarik Gheit, Z2hlaXR0QGlhcmMuZnI=