Transforming properties of the cottontail rabbit papillomavirus oncoproteins Le6 and SE6 and of the E8 protein

The rabbit papillomavirus model: a valuable tool to study viral-host interactions

Cottontail rabbit papillomavirus (CRPV) was the first DNA virus shown to be tumorigenic. The virus has since been renamed and is officially known as Sylvilagus floridanus papillomavirus 1 (SfPV1). Since its inception as a surrogate preclinical model for high-risk human papillomavirus (HPV) infections, the SfPV1/rabbit model has been widely used to study viral-host interactions and has played a pivotal role in the successful development of three prophylactic virus-like particle vaccines. In this review, we will focus on the use of the model to gain a better understanding of viral pathogenesis, gene function and host immune responses to viral infections. We will discuss the application of the model in HPV-associated vaccine testing, in therapeutic vaccine development (using our novel HLA-A2.1 transgenic rabbits) and in the development and validation of novel anti-viral and anti-tumour compounds. Our goal is to demonstrate the role the SfPV1/rabbit model has played, and continues to play, in helping to unravel the intricacies of papillomavirus infections and to develop tools to thwart the disease. This article is part of the theme issue 'Silent cancer agents: multi-disciplinary modelling of human DNA oncoviruses'.

The Papillomavirus Episteme: a major update to the papillomavirus sequence database

The Papillomavirus Episteme (PaVE) is a database of curated papillomavirus genomic sequences, accompanied by web-based sequence analysis tools. This update describes the addition of major new features. The papillomavirus genomes within PaVE have been further annotated, and now includes the major spliced mRNA transcripts. Viral genes and transcripts can be visualized on both linear and circular genome browsers. Evolutionary relationships among PaVE reference protein sequences can be analysed using multiple sequence alignments and phylogenetic trees. To assist in viral discovery, PaVE offers a typing tool; a simplified algorithm to determine whether a newly sequenced virus is novel. PaVE also now contains an image library containing gross clinical and histopathological images of papillomavirus infected lesions. Database URL: https://pave.niaid.nih.gov/.

Cutaneous papillomavirus E6 proteins must interact with p300 and block p53-mediated apoptosis for cellular immortalization and tumorigenesis

The binding of the papillomavirus E6 protein to E6AP and the induction of p53 degradation are common features of high-risk genital human papillomaviruses (HPV); cutaneous HPVs, on the other hand, lack these capacities. Nevertheless, several cutaneous HPV types of the beta-genus, such as HPV38 are associated with tumor formation when combined with genetic predisposition, immunosuppression, or UV exposure. In an animal model system, the cottontail rabbit papillomavirus (CRPV) rapidly induces skin cancer without additional cofactors, and CRPVE6 and E7 immortalize rabbit keratinocytes in vitro. However, CRPVE6 neither interacts with E6AP and p53 nor does it induce p53 degradation. In this study, we show that the interaction of CRPVE6, or HPV38E6, with the histone acetyltransferase p300 is crucial to inhibit the ability of p53 to induce apoptosis. Strikingly, E6 mutants deficient for p300 binding are incapable of preventing p53 acetylation, p53-dependent transcription, and apoptosis induction. Moreover, E6 mutants deficient for p300 binding cannot contribute to HPV38-induced immortalization of human keratinocytes or CRPV-induced tumor formation. Our findings highlight changes in the p53 acetylation status mediated by the viral E6 protein as a crucial requirement in the ability of high-risk cutaneous papillomaviruses to immortalize primary keratinocytes and induce tumors. Cancer Res; 70(17); 6913-24. (c)2010 AACR.

Papillomavirus DNA complementation in vivo

Recent phylogenic studies indicate that DNA recombination could have occurred in ancient papillomavirus types. However, no experimental data are available to demonstrate this event because of the lack of human papillomavirus infection models. We have used the cottontail rabbit papillomavirus (CRPV)/rabbit model to study pathogenesis and immunogenicity of different mutant genomes in vivo. Although the domestic rabbit is not a natural host for CRPV infection, it is possible to initiate infection with naked CRPV DNA cloned into a plasmid and monitor papilloma outgrowth on these animals. Taking advantage of a large panel of mutants based on a CRPV strain (Hershey CRPV), we tested the hypothesis that two non-viable mutant genomes could induce papillomas by either recombination or complementation. We found that co-infection with a dysfunctional mutant with an E2 transactivation domain mutation and another mutant with an E7 ATG knock out generated papillomas in rabbits. DNA extracted from these papillomas contained genotypes from both parental genomes. Three additional pairs of dysfunctional mutants also showed similar results. Individual wild type genes were also shown to rescue the function of corresponding dysfunctional mutants. Therefore, we suggest that complementation occurred between these two non-viable mutant PV genomes in vivo.

The E7 protein of the cottontail rabbit papillomavirus immortalizes normal rabbit keratinocytes and reduces pRb levels, while E6 cooperates in immortalization but neither degrades p53 nor binds E6AP

Human papillomaviruses (HPVs) cause cervical cancer and are associated with the development of non-melanoma skin cancer. A suitable animal model for papillomavirus-associated skin carcinogenesis is the infection of domestic rabbits with the cottontail rabbit papillomavirus (CRPV). As the immortalizing activity of CRPV genes in the natural target cells remains unknown, we investigated the properties of CRPV E6 and E7 in rabbit keratinocytes (RK) and their influence on the cell cycle. Interestingly, CRPV E7 immortalized RK after a cellular crisis but showed no such activity in human keratinocytes. Co-expressed CRPV E6 prevented cellular crisis. The HPV16 or CRPV E7 protein reduced rabbit pRb levels thereby causing rabbit p19(ARF) induction and accumulation of p53 without affecting cellular proliferation. Both CRPV E6 proteins failed to degrade rabbit p53 in vitro or to bind E6AP; however, p53 was still inducible by mitomycin C. In summary, CRPV E7 immortalizes rabbit keratinocytes in a species-specific manner and E6 contributes to immortalization without directly affecting p53.

Genetics of epidermodysplasia verruciformis: Insights into host defense against papillomaviruses

Epidermodysplasia verruciformis (EV) is a rare autosomal recessive genodermatosis associated with a high risk of skin carcinoma. EV results from an abnormal susceptibility to infection by specific human papillomavirus (HPV) genotypes (beta-papillomaviruses) which include the potentially oncogenic HPV5. EV-specific HPVs are considered as harmless for the general population. EV was recently found to be caused by invalidating mutations in two adjacent, related, novel genes, EVER1/TMC6 and EVER2/TMC8. EVER genes encode transmembrane proteins located in the endoplasmic reticulum, which are likely to function as modifiers of ion transporters or channels and to be involved in signal transduction. It was proposed that EV was a primary defect of innate immunity. Our hypothesis is that EVER proteins act as restriction factors for EV-specific HPVs in keratinocytes, and that EV represents a primary deficiency of intrinsic immunity against certain papillomaviruses.

Papillomavirus genome structure, expression, and post-transcriptional regulation

Papillomaviruses are a group of small non-enveloped DNA tumor viruses whose infection usually causes benign epithelial lesions (warts). Certain types of HPVs, such as HPV-16, HPV-18, and HPV-31, have been recognized as causative agents of cervical cancer and anal cancer and their infections, which arise via sexual transmission, are associated with more than 95% of cervical cancer. Papillomaviruses infect keratinocytes in the basal layer of stratified squamous epithelia and replicate in the nucleus of infected keratinocytes in a differentiation-dependent manner. Viral gene expression in infected cells depends on cell differentiation and is tightly regulated at the transcriptional and post-transcriptional levels. A noteworthy feature of all papillomavirus transcripts is that they are transcribed as a bicistronic or polycistronic form containing two or more ORFs and are polyadenylated at either an early or late poly(A) site. In the past ten years, remarkable progress has been made in understanding how this complex viral gene expression is regulated at the level of transcription (such as via DNA methylation) and particularly post-transcription (including RNA splicing, polyadenylation, and translation). Current knowledge of papillomavirus mRNA structure and RNA processing has provided some clues on how to control viral oncogene expression. However, we still have little knowledge about which mRNAs are used to translate each viral protein. Continuing research on post-transcriptional regulation of papillomavirus infection will remain as a future focus to provide more insights into papillomavirus-host interactions, the virus life-cycle, and viral oncogenesis.

Cottontail rabbit papillomavirus E8 protein is essential for wart formation and provides new insights into viral pathogenesis

The cottontail rabbit papillomavirus (CRPV) a and b subtypes display a conserved E8 open reading frame encoding a 50-amino-acid hydrophobic protein, with structural similarities to the E5 transmembrane oncoprotein of genital human PVs (HPVs). CRPV E8 has been reported to play a role in papilloma growth but not to be essential in papilloma formation. Here we report that the knockout of E8 start codon almost prevented wart induction upon biobalistic inoculation of viral DNA onto rabbit skin. The scarce warts induced showed very slow growth, despite sustained expression of E6 and E7 oncogenes. This points to an essential role of E8 in disturbing epidermal homeostasis. Using a yeast two-hybrid screen, we found that E8 interacted with the zinc transporter ZnT1, protocadherin 1 (PCDH1), and AHNAK/desmoyokin, three proteins as yet unrelated to viral pathogenesis or cell transformation. HPV16 E5 also interacted with these proteins in two-hybrid assay. CRPV E8 mainly localized to the Golgi apparatus and the early endosomes of transfected keratinocytes and colocalized with ZnT1, PCDH1, and AHNAK. We showed that ZnT1 and PCDH1 formed a complex and that E8 disrupted this complex. CRPV E8, like HPV16 E5, increased epidermal growth factor (EGF)-dependent extracellular signal-regulated kinase 1/2 (ERK1/2) phosphorylation and both the EGF-dependent and the EGF-independent activity of activating protein-1 (AP-1). Competition experiments with a nonfunctional truncated ZnT1 protein showed that E8-ZnT1 interaction was required for AP-1 activation. Our data identify CRPV E8 as a key player in papilloma induction and unravel novel cellular targets for inducing the proliferation of keratinocytes.

Association of cottontail rabbit papillomavirus E6 oncoproteins with the hDlg/SAP97 tumor suppressor

Papillomaviruses are small DNA viruses that infect epithelial tissues and cause warts. Human papillomavirus (HPV) infection is the primary risk factor for the development of cervical cancer. The E6 and E7 oncogenes are the only genes consistently expressed in HPV-positive cervical cancer cells. Cottontail rabbit papillomavirus (CRPV) induces papillomas and carcinomas on cottontail and domestic rabbits and provides an excellent animal model of HPV infection and vaccine development. CRPV encodes three transforming proteins; LE6, SE6, and E7. Each of these proteins is required for papilloma formation. Like HPV E7, the CRPV E7 protein binds to the tumor suppressor pRB. In contrast, unlike HPV E6, the CRPV E6 proteins do not bind the tumor suppressor p53. Although more than a dozen cellular proteins have been identified as HPV E6 interacting proteins, nothing is known about the cellular interacting proteins of CRPV E6s. Here we describe the association of CRPV E6s with hDlg/SAP97, the mammalian homolog of the Drosophila discs large tumor suppressor protein. HPV E6 has previously shown to bind and target hDlg/SAP97 for degradation. Our results demonstrate that both LE6 and SE6 interact with hDlg/SAP97, although their association does not lead to the degradation of hDlg/SAP97. The PDZ domains of hDlg were shown to be sufficient for interaction with CRPV E6 proteins while the C-terminus of CRPV E6 is essential for the interaction with hDlg. The association of hDlg with SE6 may be important but not sufficient for the transformation of NIH 3T3 cells by SE6. Importantly, a CRPV SE6 mutant defective for papilloma formation did not interact with hDlg. These results suggest that interaction with hDlg/SAP97 plays a role in the biological function of CRPV E6s.

Characterization of three rabbit oral papillomavirus oncogenes

Rabbit oral papillomavirus (ROPV) induces warts in mucosal tissues, and represents a useful model for understanding host-virus interactions that are reflected in mucosal/HPV infections. ROPV induces benign papillomas that regress in 100% of infected rabbits. We previously reported the complete genome sequence of ROPV. However, the oncogenic potential of this virus is unknown because of immunologically mediated regression. The purpose of this study was to characterize the transforming proteins of E6, E7, and E8 genes of ROPV. E6, E7, and E8 genes of ROPV were cloned into the expression vector PCR3. Two hybrid CRPV-ROPV E6 genes were also constructed and tested together with the three wild-type ROPV genes. Each construct was transfected into NIH3T3 cells and stable transfected cell lines were established. Transforming properties of ROPV E6, E7, and E8 were tested via anchorage-independent growth of cells in agar plates and tumor growth in athymic mice. Cells with ROPV E6, E7, or E8 formed colonies in agar and tumors in athymic mice. These observations suggest that ROPV E6, E7, and E8 are oncogenic.

GM-CSF enhances protective immunity to cottontail rabbit papillomavirus E8 genetic vaccination in rabbits

We have reported previously that cottontail rabbit papillomavirus (CRPV) E8 gene immunization induced strong protection against virus challenge. In this study, we primed E8 gene vaccination with mouse granulocyte-macrophage colony-stimulating factor (mGM-CSF), a cytokine that induces differentiation and local recruitment of professional antigen-presenting cells. EIII/JC inbred rabbits were divided into four groups receiving vaccinations with the following constructs: mGM-CSF plus E8, mGM-CSF only, E8 only and vector only. After three immunizations at intervals of 3 weeks, rabbits were challenged with viral DNA at six scarified sites. Papillomas grew on all vaccinated rabbits 4 weeks after inoculation. At week 5, papillomas on four rabbits of mGM-CSF plus E8 and one of E8 only rabbits began to regress. At week 11, all the papillomas on rabbits in the GM-CSF plus E8 vaccination group regressed (regression rate = 100%); regression rates of the mGM-CSF only and E8 only vaccination groups were 50 and 25%, respectively. All papillomas on the vector immunized rabbits remained persistent until the end of the experiment (0%). Antibodies to mGM-CSF were detected in rabbit serum by Western blot. Rabbits vaccinated with E8 plus mGM-CSF or E8 only group had positive Delayed-type hypersensitivity (DTH) skin test to different E8 peptides. These results demonstrated that mGM-CSF could enhance the effects of E8 immunization in rabbits to CRPV infection through cell-mediated immune responses.

Amino acid residues in the carboxy-terminal region of cottontail rabbit papillomavirus E6 influence spontaneous regression of cutaneous papillomas

Previous studies have identified two different strains of cottontail rabbit papillomavirus (CRPV) that differ by approximately 5% in base pair sequence and that perform quite differently when used to challenge New Zealand White (NZW) rabbit skin. One strain caused persistent lesions (progressor strain), and the other induced papillomas that spontaneously regressed (regressor strain) at high frequencies (J. Salmon, M. Nonnenmacher, S. Caze, P. Flamant, O. Croissant, G. Orth, and F. Breitburd, J. Virol. 74:10766-10777, 2000; J. Salmon, N. Ramoz, P. Cassonnet, G. Orth, and F. Breitburd, Virology 235:228-234, 1997). We generated a panel of CRPV genomes that contained chimeric and mutant progressor and regressor strain E6 genes and assessed the outcome upon infection of both outbred and EIII/JC inbred NZW rabbits. The carboxy-terminal 77-amino-acid region of the regressor CRPV strain E6, which contained 15 amino acid residues that are different from those of the equivalent region of the persistent CRPV strain E6, played a dominant role in the conversion of the persistent CRPV strain to one showing high rates of spontaneous regressions. In addition, a single amino acid change (G252E) in the E6 protein of the CRPV progressor strain led to high frequencies of spontaneous regressions in inbred rabbits. These observations imply that small changes in the amino acid sequences of papillomavirus proteins can dramatically impact the outcome of natural host immune responses to these viral infections. The data imply that intrastrain differences between separate isolates of a single papillomavirus type (such as human papillomavirus type 16) may contribute to a collective variability in host immune responses in outbred human populations.

Intracutaneous DNA vaccination with the E8 gene of cottontail rabbit papillomavirus induces protective immunity against virus challenge in rabbits

The cottontail rabbit papillomavirus (CRPV)-rabbit model has been used in several studies for testing prophylactic and therapeutic papillomavirus vaccines. Earlier observations had shown that the CRPV nonstructural genes E1, E2, and E6 induced strong to partial protective immunity against CRPV infection. In this study, we found that CRPV E8 immunization eliminated virus-induced papillomas in EIII/JC inbred rabbits (100%) and provided partial protection (55%) against virus challenge in outbred New Zealand White rabbits. CRPV-E8 is a small open reading frame, coding for a 50-amino-acid protein, that is colinear with the CRPV E6 gene and has features similar to those of the bovine papillomavirus and human papillomavirus E5 genes. Papillomas that grew on E8-vaccinated outbred rabbits were significantly smaller than those on vector-vaccinated rabbits (P < 0.01; t test). Delayed-type hypersensitivity skin tests showed that some of the E8-vaccinated rabbits had positive responses to E8-specific peptides.

Variation in the nucleotide sequence of cottontail rabbit papillomavirus a and b subtypes affects wart regression and malignant transformation and level of viral replication in domestic rabbits

We previously reported the partial characterization of two cottontail rabbit papillomavirus (CRPV) subtypes with strikingly divergent E6 and E7 oncoproteins. We report now the complete nucleotide sequences of these subtypes, referred to as CRPVa4 (7,868 nucleotides) and CRPVb (7,867 nucleotides). The CRPVa4 and CRPVb genomes differed at 238 (3%) nucleotide positions, whereas CRPVa4 and the prototype CRPV differed by only 5 nucleotides. The most variable region (7% nucleotide divergence) included the long regulatory region (LRR) and the E6 and E7 genes. A mutation in the stop codon resulted in an 8-amino-acid-longer CRPVb E4 protein, and a nucleotide deletion reduced the coding capacity of the E5 gene from 101 to 25 amino acids. In domestic rabbits homozygous for a specific haplotype of the DRA and DQA genes of the major histocompatibility complex, warts induced by CRPVb DNA or a chimeric genome containing the CRPVb LRR/E6/E7 region showed an early regression, whereas warts induced by CRPVa4 or a chimeric genome containing the CRPVa4 LRR/E6/E7 region persisted and evolved into carcinomas. In contrast, most CRPVa, CRPVb, and chimeric CRPV DNA-induced warts showed no early regression in rabbits homozygous for another DRA-DQA haplotype. Little, if any, viral replication is usually observed in domestic rabbit warts. When warts induced by CRPVa and CRPVb virions and DNA were compared, the number of cells positive for viral DNA or capsid antigens was found to be greater by 1 order of magnitude for specimens induced by CRPVb. Thus, both sequence variation in the LRR/E6/E7 region and the genetic constitution of the host influence the expression of the oncogenic potential of CRPV. Furthermore, intratype variation may overcome to some extent the host restriction of CRPV replication in domestic rabbits.

Characterization of transformation function of cottontail rabbit papillomavirus E5 and E8 genes

Cottontail rabbit papillomavirus (CRPV) induces rabbit skin papillomas, which progress to invasive carcinoma in some animals. Two early genes, E7 and E6, have been demonstrated previously to be oncogenes. In this study, we identified two additional transforming genes, E8 and E5. Both E8 and E5 stimulated C127 and BALB/c A31 (A31) cell proliferation and affected cell cycle transition. The E8 and E5 transfectants lost cell contact inhibition, reaching a high saturation density when cultured up to 2 weeks. E8-C127 transfectants formed colonies in soft agar in the presence of platelet-derived growth factor (PDGF) while E5-C127 transfectants formed colonies without the requirement for PDGF. E8-C127 transfectants were highly tumorigenic whereas E5-C127 transfectants showed a weak tumorigenicity in nude mice. Both E8 and E5 A31 transfectants failed to form colonies in soft agar even in the presence of platelet-derived growth factor (PDGF) and did not develop tumors in nude mice. These results clearly showed that CRPVE8 and E5 are oncogenes and that the PDGF beta-receptor signaling pathway may be involved in E8-mediated C127 cell transformation. The difference in colony formation in soft agar and tumorigenicity in nude mice between C127 and A31 cell lines indicates that additional alterations in cellular gene expression are needed for E5- and E8-transfected cells to acquire a malignant phenotype.

Molecular targets for human papillomaviruses: prospects for antiviral therapy

A substantial medical need exists for the development of antiviral medicines for the treatment of diseases associated with infection by human papillomaviruses (HPVs). HPVs are associated with various benign and malignant lesions including benign genital condyloma, common skin warts, laryngeal papillomas and anogenital cancer. Since treatment options are limited and typically not very satisfactory, the development of safe and effective antiviral drugs for HPV could have substantial clinical impact. In the last few years, exciting advances have been made in our understanding of papillomavirus replication and the effects that the virus has on growth of the host cell. Although still somewhat rudimentary, techniques have been developed for limited virion production in vitro offering the promise of more rapid advances in the dissection and understanding of the virus life cycle. Of the 8-10 HPV gene products that are made during infection, only one encodes enzymatic activities, the E1 helicase. Successful antiviral therapies have traditionally targeted viral enzymes such as polymerases, kinases and proteases. In contrast, macromolecular interactions which mediate the functions of E6, E7 and E2 are thought to be more difficult targets for small molecule therapy.

The E6 and E7 genes of human papilloma virus-type 16 protect primary astrocyte cultures from injury

Many oncogenes are implicated in the regulation of apoptosis as well as in control of the cell cycle, and several have been shown to protect cells from injury. We tested whether the expression of human papillomavirus type 16 genes E6 and E7 could protect primary astrocytes from injury. Retroviral vectors were used to express E6, E7, or E6E7 in primary murine astrocyte cultures. Astrocytes expressing E6E7 suffered less than half the injury seen in controls after exposure to 400 microM H2O2. When we compared astrocytes expressing only E6 or E7, cells expressing E7 alone were protected to a greater extent and from more severe injury than those expressing E6. E6E7, E6 and E7 all provided protection from 30 h glucose deprivation, but again E7 provided the best protection, reducing injury to less than a third of that seen in controls. Expression of E7 alone decreases vulnerability to both hydrogen peroxide and glucose deprivation injury while E6 alone primarily decreases glucose deprivation injury. These results demonstrate that expression of human papillomavirus oncogenes can reduce the vulnerability of astrocytes to oxidative stress and nutrient deprivation.

A comparative analysis of the interactions of the E6 proteins from cutaneous and genital papillomaviruses with p53 and E6AP in correlation to their transforming potential

Common necessity for all papillomaviruses is to induce DNA synthesis in quiescent cells. This is commonly achieved by the E7 gene product, which interferes with the function of members of the retinoblastoma family controlling transition from the G1-phase to the S-phase of the cell cycle. Uncontrolled entry into S-phase activates, however, negative growth control signals which have to be bypassed to achieve production of progeny viruses. In addition to inherent activities of the E7 protein, high risk genital types encode an E6 protein that overcomes p53-mediated G1-arrest and apoptosis in concert with the cellular factor E6AP by targeting p53 for the enhanced ubiquitin-dependent degradation. The key question, which of these functions of genital E6 and E7 proteins is responsible for the carcinogenic phenotype, is still not completely answered. In contrast to high risk genital types no immortalizing or transforming activities have been found for the E7 proteins of the high risk cutaneous HPV8 and 47. On the other hand the ability of the E6 protein to transform established rodent fibroblasts seems to be a property shared by high risk genital and cutaneous types. To examine the existence of a common E6-mediated transforming pathway for both virus groups we compared the properties of the cutaneous E6 proteins with already known functions of E6 proteins of genital viruses. For this we analyzed the E6 proteins of low nak and high risk cutaneous and genital papillomaviruses with respect to cell transformation, to their abilities to bind, degradate, and influence the activity of human p53, and to bind E6AP. The results of our study demonstrate a clear lack of interaction between the transforming E6 proteins of HPV1 and HPV8 and both cellular proteins p53 and E6AP. In contrast, we found E6AP-independent binding of HPV16 E6 and HPV6 E6 to p53, although both proteins were different in their transforming potential. Of all four proteins investigated, only HPV16 E6 was able to bind to p53 and E6AP and to induce degradation of the p53 protein in the reticulocyte system. When we investigated in frame deletion mutants of the E6 protein of HPV16 for their abilities to bind to p53 or E6AP, degradate, and inhibit the transactivation function of p53 and to transform rodent fibroblasts, no correlation between the different activities could be found. Mutants still able to bind p53 and E6AP lacked transforming ability and other mutants that were transformation-competent were deficient in p53 and E6AP bindings.

A cottontail rabbit papillomavirus strain (CRPVb) with strikingly divergent E6 and E7 oncoproteins: An insight in the evolution of papillomaviruses

We previously observed that warts induced by an isolate of cottontail rabbit papillomavirus (CRPV) showed incomplete instead of systemic regression in some domestic rabbits. We report that the viral isolate contained, as a major component, a CRPV strain (CRPVb) showing an unexpectedly high divergence in the E6 and E7 open reading frames (ORFs), compared to the prototype CRPVa present in the isolate as a minor component. The E6 and E7 oncoproteins of CRPVa and -b disclosed only 87.5% identical amino acids and differed in size by three and two amino acids, respectively. This divergence involved (i) a great number (4.4%) of nucleotide substitutions and a high rate (83.3%) of nonsynonymous mutations; (ii) mutations changing the E6 and E7 stop codons; and (iii) in-frame sequence insertions in the E6 ORF (18 nucleotides) and downstream of the mutated E7 stop codon (6 nucleotides), both likely to result from a duplication of adjacent sequences. These extensive differences could account for distinct biological and antigenic properties. Strikingly, only four (0.8%) amino acids of the L1 major capsid protein were variable. Thus, it seems likely that sequence duplications and mutations affecting stop codons exert a strong selection pressure on the fixation of nonsynonymous mutations and that phylogenetic calculations based only on point mutations may misevaluate the time scale of the evolution of papillomaviruses.