Nucleotide sequence and phylogenetic classification of human papillomavirus type 59

Phylogenetic analysis of Alphapapillomavirus based on L1, E6 and E7 regions suggests that carcinogenicity and tissue tropism have appeared multiple times during viral evolution

Members of the Alphapapillomavirus genus are causative agents for cervix cancer and benign lesions in humans. These viruses are classified according to sequence similarities in their L1 region. Yet, viral carcinogenicity has been associated with variations in the proteins encoded by the E6 and E7 genes. In order to relate evolutionary history with origin of carcinogenicity, we performed phylogenetic reconstructions using both nucleotide and predicted amino acid sequences of the L1, E6 and E7 genes. Whilst phylogenetic analysis of L1 reconstructed genus evolutionary history, phylogenies based on E6 and E7 proteins support the idea that mutations at amino acids S/Tx [V/L] (E6) and LxCxE (E7) might be responsible for carcinogenic potential. These findings indicate that virulence within Alphapapillomavirus have appeared multiple times during evolution. Our results reveal that oncogenic potential is not a monophyletic clade-specific adaptation but might be the result of positive selection on random mutations occurring on proteins involved in host infection during viral diversification.

Evolution and taxonomic classification of alphapapillomavirus 7 complete genomes: HPV18, HPV39, HPV45, HPV59, HPV68 and HPV70

BACKGROUND

The species Alphapapillomavirus 7 (alpha-7) contains human papillomavirus genotypes that account for 15% of invasive cervical cancers and are disproportionately associated with adenocarcinoma of the cervix. Complete genome analyses enable identification and nomenclature of variant lineages and sublineages.

METHODS

The URR/E6 region was sequenced to screen for novel variants of HPV18, 39, 45, 59, 68, 70, 85 and 97 from 1147 cervical samples obtained from multiple geographic regions that had previously been shown to contain an alpha-7 HPV isolate. To study viral heterogeneity, the complete 8 kb genome of 128 isolates, including 109 sequenced for this analysis, were annotated and analyzed. Viral evolution was characterized by constructing phylogenic trees using maximum-likelihood and Bayesian algorithms. Global and pairwise alignments were used to calculate total and ORF/region nucleotide differences; lineages and sublineages were assigned using an alphanumeric system. The prototype genome was assigned to the A lineage or A1 sublineage.

RESULTS

The genomic diversity of alpha-7 HPV types ranged from 1.1% to 6.7% nucleotide sequence differences; the extent of genome-genome pairwise intratype heterogeneity was 1.1% for HPV39, 1.3% for HPV59, 1.5% for HPV45, 1.6% for HPV70, 2.1% for HPV18, and 6.7% for HPV68. ME180 (previously a subtype of HPV68) was designated as the representative genome for HPV68 sublineage C1. Each ORF/region differed in sequence diversity, from most variable to least variable: noncoding region 1 (NCR1) / noncoding region 2 (NCR2) > upstream regulatory region (URR) > E6 / E7 > E2 / L2 > E1 / L1.

CONCLUSIONS

These data provide estimates of the maximum viral genomic heterogeneity of alpha-7 HPV type variants. The proposed taxonomic system facilitates the comparison of variants across epidemiological and molecular studies. Sequence diversity, geographic distribution and phylogenetic topology of this clinically important group of HPVs suggest an independent evolutionary history for each type.

Evolutionary dynamics of variant genomes of human papillomavirus types 18, 45, and 97

Human papillomavirus type 18 (HPV18) and HPV45 account for approximately 20% of all cervix cancers. We show that HPV18, HPV45, and the recently discovered HPV97 comprise a clade sharing a most recent common ancestor within HPV alpha7 species. Variant lineages of these HPV types were classified by sequence analysis of the upstream regulatory region/E6 region among cervical samples from a population-based study in Costa Rica, and 27 representative genomes from each major variant lineage were sequenced. Nucleotide variation within HPV18 and HPV45 was 3.82% and 2.39%, respectively, and amino acid variation was 4.73% and 2.87%, respectively. Only 18 nucleotide variations, of which 10 were nonsynonymous, were identified among three HPV97 genomes. Full-genome comparisons revealed maximal diversity between HPV18 African and non-African variants (2.6% dissimilarity), whereas HPV18 Asian-American [E1 (AA)] and European (E2) variants were closely related (less than 0.5% dissimilarity); HPV45 genomes had a maximal difference of 1.6% nucleotides. Using a Bayesian Markov chain Monte Carlo (MCMC) method, the divergence times of HPV18, -45, and -97 from their most recent common ancestors indicated that HPV18 diverged approximately 7.7 million years (Myr) ago, whereas HPV45 and HPV97 split off around 5.7 Myr ago, in a period encompassing the divergence of the great ape species. Variants within the HPV18/45/97 lineages were estimated to have diverged from their common ancestors in the genus Homo within the last 1 Myr (<0.7 Myr). To investigate the molecular basis of HPV18, HPV45, and HPV97 evolution, regression models of codon substitution were used to identify lineages and amino acid sites under selective pressure. The E5 open reading frame (ORF) of HPV18 and the E4 ORFs of HPV18, HPV45, and HPV18/45/97 had nonsynonymous/synonymous substitution rate ratios (d(N)/d(S)) over 1 indicative of positive Darwinian selection. The L1 ORF of HPV18 genomes had an increased proportion of nonsynonymous substitutions (4.93%; average d(N)/d(S) ratio [M3] = 0.3356) compared to HPV45 (1.86%; M3 = 0.1268) and HPV16 (2.26%; M3 = 0.1330) L1 ORFs. In contrast, HPV18 and HPV16 genomes had similar amino acid substitution rates within the E1 ORF (2.89% and 3.24%, respectively), while HPV45 E1 was highly conserved (amino acid substitution rate was 0.77%). These data provide an evolutionary history of this medically important clade of HPVs and identify an unexpected divergence of the L1 gene of HPV18 that may have clinical implications for the long-term use of an L1-virus-like particle-based prophylactic vaccine.

Progression from productive infection to integration and oncogenic transformation in human papillomavirus type 59-immortalized foreskin keratinocytes

Studies of changes in the virus and host cell upon progression from human papillomavirus (HPV) episomal infection to integration are critical to understanding HPV-related malignant transformation. However, there exist only a few in vitro models of both productive HPV infection and neoplastic progression on the same host background. We recently described a unique foreskin keratinocyte cell line (ERIN 59) that contains HPV 59 (a close relative of HPV 18). Early passages of ERIN 59 cells (passages 9-13) contained approximately 50 copies of episomes/cell, were feeder cell-dependent, and could be induced to differentiate and produce infectious virus in a simple culture system. We now report that late passage cells (passages greater than 50) were morphologically different from early passage cells, were feeder cell independent, and did not differentiate or produce virus. These late passage cells contained HPV in an integrated form. An integration-derived oncogene transcript was expressed in late passage cells. The E2 open reading frame was interrupted in this transcript at nucleotide 3351. Despite a lower viral genome copy number in late passage ERIN 59 cells, expression of E6/E7 oncogene transcripts was similar to early passage cells. We conclude that ERIN 59 cells are a valuable cell line representing a model of progression from HPV 59 episomal infection and virus production to HPV 59 integration and associated oncogenic transformation on the same host background.

Human papillomavirus type 59 immortalized keratinocytes express late viral proteins and infectious virus after calcium stimulation

Human papillomavirus type 59 (HPV 59) is an oncogenic type related to HPV 18. HPV 59 was recently propagated in the athymic mouse xenograft system. A continuous keratinocyte cell line infected with HPV 59 was created from a foreskin xenograft grown in an athymic mouse. Cells were cultured beyond passage 50. The cells were highly pleomorphic, containing numerous abnormally shaped nuclei and mitotic figures. HPV 59 sequences were detected in the cells by DNA in situ hybridization in a diffuse nuclear distribution. Southern blots were consistent with an episomal state of HPV 59 DNA at approximately 50 copies per cell. Analysis of the cells using a PCR/reverse blot strip assay, which amplifies a portion of the L1 open reading frame, was strongly positive. Differentiation of cells in monolayers was induced by growth in F medium containing 2 mM calcium chloride for 10 days. Cells were harvested as a single tissue-like sheet, and histologic analysis revealed a four-to-six cell-thick layer. Transcripts encoding involucrin, a cornified envelope protein, and the E1/E4 and E1/E4/L1 viral transcripts were detected after several days of growth in F medium containing 2 mM calcium chloride. The E1/E4 and L1 proteins were detected by immunohistochemical analysis, and virus particles were seen in electron micrographs in a subset of differentiated cells. An extract of differentiated cells was prepared by vigorous sonication and was used to infect foreskin fragments. These fragments were implanted into athymic mice. HPV 59 was detected in the foreskin xenografts removed 4 months later by DNA in situ hybridization and PCR/reverse blot assay. Thus, the complete viral growth cycle, including production on infectious virus, was demonstrated in the HPV 59 immortalized cells grown in a simple culture system.

Human papillomavirus and cervical cancer

Of the many types of human papillomavirus (HPV), more than 30 infect the genital tract. The association between certain oncogenic (high-risk) strains of HPV and cervical cancer is well established. Although HPV is essential to the transformation of cervical epithelial cells, it is not sufficient, and a variety of cofactors and molecular events influence whether cervical cancer will develop. Early detection and treatment of precancerous lesions can prevent progression to cervical cancer. Identification of precancerous lesions has been primarily by cytologic screening of cervical cells. Cellular abnormalities, however, may be missed or may not be sufficiently distinct, and a portion of patients with borderline or mildly dyskaryotic cytomorphology will have higher-grade disease identified by subsequent colposcopy and biopsy. Sensitive and specific molecular techniques that detect HPV DNA and distinguish high-risk HPV types from low-risk HPV types have been introduced as an adjunct to cytology. Earlier detection of high-risk HPV types may improve triage, treatment, and follow-up in infected patients. Currently, the clearest role for HPV DNA testing is to improve diagnostic accuracy and limit unnecessary colposcopy in patients with borderline or mildly abnormal cytologic test results.

Identification of six putative novel human papillomaviruses (HPV) and characterization of candidate HPV type 87

Six putative novel human papillomavirus (HPV) types were detected by using general primers for a conserved L1 HPV region in patients examined in gynecologic centers. One of the isolates, detected in samples from 4 patients with koilocytic atypia at cervical cytology (3 of whom were also infected with human immunodeficiency virus type 1), was completely sequenced, identified as a new HPV genotype, and designated candidate HPV87 (candHPV87) by the Reference Center for Human Papillomavirus. candHPV87 shows the classic HPV genome organization and the absence of a functional E5 coding region. Phylogenetic analysis documented that the candHPV87 genome clusters within the A3 group of HPVs, together with HPV61, HPV72, HPV83, HPV84 and candHPV86, which have been completely sequenced, and a number of other putative novel genotypes (two of which are described in this work), which have been partially characterized. To address the growth-enhancing potential of candHPV87, the E6 and E7 putative coding regions were cloned and expressed in tissue cultures. The data indicate that both proteins stimulate cell division in tissue cultures more than those of low-risk HPVs, though not as much as those of HPV16. Taken together, the clinical, molecular, and biological data suggest that the novel papillomavirus characterized in the present study is a low- to intermediate-risk HPV.

Relationships between 80 human papillomavirus genotypes and different grades of cervical intraepithelial neoplasia: association and causality

To clarify the causal relationship between human papillomavirus (HPV) and cervical intraepithelial neoplasia (CIN), we analyzed 386 unfixed biopsy specimens by blot hybridization at Tm -40 degrees C, targeting 38 skin and 42 genital HPVs. By the recognition of PstI, BanI, and MspI cleavage patterns, single genital, but no skin's HPVs were identified with more than 10 copies per cell in 354 CIN (88 CIN I, 94 CIN II, and 172 CIN III). HPVs 40, 42, 43, 54, 62, or 71 was found in 10 CIN I, while HPVs 18, 30, 39, 51, 56, 59, 66, 68, 69, or 82 was found in 35 CIN I, 20 CIN II, or 8 CIN III. On the other hand, HPVs 16, 31, 33, 35, 52, 58, or 67 was identified in 43 CIN I, 74 CIN II, or 164 CIN III. The results are strongly indicative that most genital HPVs have potency to induce CIN I; however, HPV 16 and its closely related types are able to efficiently induce CIN III. We discuss the definition of causal HPV for CIN with regard to viral prevalence and viral load.

Propagation of human papillomavirus type 59 in the athymic mouse xenograft system

Studies of human papillomavirus (HPV) infection are hampered by the lack of an adequate culture system. The athymic mouse xenograft system permits propagation of HPV, but only a few HPV types have been grown in this manner. To produce an oncogenic type for studies of HPV pathogenesis, a condylomata acuminata lesion from an immunosuppressed patient was used to prepare an infectious extract. The patient's lesion was shown by PCR analysis to contain abundant HPV 59 (an oncogenic type) and a lesser amount of HPV 11, a nononcogenic type. The extract was used to infect human foreskin tissue which was subsequently implanted into athymic mice. Characterization of implants recovered after 3-4 months of growth revealed the presence of HPV 59 exclusively. A second extract was prepared from one of these implants and used in an additional experiment to demonstrate the passage of HPV 59. Compared to the histopathologic changes induced by the prototypic nononcogenic HPV 11, infection with HPV 59 caused a higher degree of basal cell crowding, less acanthosis, minimal papillomatosis and less pronounced koilocytosis.

Ubiquitous presence of E6 and E7 transcripts in human papillomavirus-positive cervical carcinomas regardless of its type

The presence of human papillomavirus (HPV) DNA in almost all of the cervical carcinomas is one of the most compelling evidence for the viral carcinogenesis. HPVs are thought to induce cervical carcinoma most likely through the expression of E6 and E7 genes presumably by inactivating the tumor suppressor proteins, p53 and pRb, respectively. Thus far, the presence of HPV E6 and E7 transcripts have been identified only in cervical carcinoma-derived cell lines harboring type 16 or 18, and in a limited number of cervical neoplasia specimens positive for type 16, 18, 33 or 51. To see whether the expression of E6 and E7 genes is an essential finding in HPV-positive cervical carcinoma and cervical intraepithelial neoplasia (CIN), we constructed a reverse transcription-polymerase chain reaction (RT-PCR) assay using a pair of consensus primers in the E6 and E7 regions. Using the assay, E6 transcripts (full-length E6/E7 transcripts) and E7 transcripts (spliced E6/E7 transcripts, E6* mRNA) were identified in 97% (30/31) and 100% (all 31) of cervical carcinomas and in 100% (all 23) and 74% (17/23) of CINs, respectively. This assay also revealed unknown splice donor and acceptor sites of E6* mRNA of less frequent HPV types 31, 35, 52, 56, 58 and 59 based on sequence analyses of the PCR products. Thus, the present study demonstrates that E6 and E7 transcripts of HPV exist in virtually all HPV-positive cervical neoplasia specimens except for the absence of E7 transcripts in some of CINs.

Complete nucleotide sequence, genomic organization and phylogenetic analysis of a novel genital human papillomavirus type, HLT7474-S

A novel human papillomavirus (HPV) type, HLT7474-S, was isolated from a cervical scraping of a female sex worker with a wart virus infection. The complete DNA sequence of 7812 bp was derived from four overlapping PCR products and authenticated by RFLP analysis. The L1 gene exhibited 78% identity to those of its most closely related known HPV types in group A7, comprising HPV types 18, 39, 45, 59, 68 and 70. The genomic organization and phylogenetic analysis of HLT7474-S and group A7 HPVs reiterated their relatedness. Of significance were the strong sequence similarity, phylogenetic relationship and conservation of critical motifs between the transforming E6 and E7 of HLT7474-S and E6 of HPV-18 and E7 of HPV-59, respectively. These features clearly suggest that HLT7474-S is a high-risk genital HPV isolate, closely related to HPV-18 and other members of the A7 group of genital HPVs.

The genomes of three of four novel HPV types, defined by differences of their L1 genes, show high conservation of the E7 gene and the URR

The DNA genomes of four new human papillomaviruses, HPV 75, HPV 76, HPV 77, and HPV 80, have been cloned, sequenced, and characterized. HPV 75, HPV 76 (both HPV 49-related), and HPV 77 (HPV 29-related) were isolated from benign cutaneous warts and HPV 80 (HPV 15-related) from histologically normal skin. HPV 77 has also been demonstrated in dysplastic warts and squamous cell carcinomas of the skin. The sequence data presented in this study led to a proposed modification of the definition of a new HPV type. The high degree of DNA sequence similarity between the E7 ORF of HPV 77 and HPV 29 (97.7%), as opposed to the E6 (82.8%) and L1 (85.3%) ORFs, might suggest conservation of a specific function or a possible recombinational event. Only the E6 and L1 ORFs of HPV 75 and HPV 76 have a similarity lower than 90%, whereas the DNA sequences of their upstream regulatory regions (URRs) share a similarity of 93%. The E7, E1, and E4 ORFs, as well as the URR of HPV 15 and HPV 80, share sequence similarities higher than 90%. Such a divergence in the similarity between different segments of the virus genomes of closely related HPV types has not been noted to date. A detailed comparative sequence analysis was performed. HPV 75, HPV 76, and HPV 80 revealed features characteristic of truly cutaneous HPV types, whereas HPV 77 shared several characteristics with the mucosal HPV types, some of which may have functional consequences.

Identification of cis-regulatory elements in the upstream regulatory region of human papillomavirus type 59

Human papillomavirus type 59 (HPV-59) was cloned from a vulvar intraepithelial neoplasia and the complete nucleotide sequence was determined. This virus is closely related to HPV-18 and -39 (60% homology in nucleotide sequence) and is grouped with the genital HPV types. In the present paper, we demonstrate that the HPV-59 E2 transactivator represses its E6 promoter-mediated transcription. We have also analyzed cis-regulatory elements in the upstream regulatory region (URR) of HPV-59 using chloramphenicol acetyl transferase assays as well as electrophoresis mobility shift assays (EMSA). The results allow for a subdivision of the HPV-59 URR into three regions of activity: distal (nt 7149-7493), central (nt 7493-7742), and proximal (nt 7742-7748). In particular, the 250 bp (nt 7493-7742) of the central region plays an important role as a constitutive enhancer element for the maximal transcription of the E6 promoter. Our results suggest that the transcription factors AP1, Oct1, SP1 and unidentified factors bind to the HPV-59 E6 promoter region, whereas NF1, GRE and TFIID fail to bind despite the presence of putative binding sites in the DNA sequence.

Transforming activities of human papillomavirus type 59 E5, E6 and E7 open reading frames in mouse C127 cells

The DNA sequence from a human papillomavirus type 59 (HPV 59) has been recently determined. The HPV 59 genome consists of 7896 nucleotides (nt). A comparative analysis of this sequence with the sequences of other HPVs revealed the closest homology to HPV 18 (71%). To test the transforming activities of HPV 59 DNA and its gene products, several plasmids expressing HPV 59 open reading frames (ORF) were constructed. The E5, E6, and E7 ORFs of HPV 59 were inserted into pRc/CMV vector containing a promoter of cytomegalovirus to test the transforming activities of these ORFs. When these DNAs were transfected into mouse C127 cells, all three ORFs were independently able to transform C127 cells in the presence of G418, although the full length HPV 59 DNAs failed to induce the focus-formation. The E7 ORF showed the strongest transforming activity and the E5 ORF exhibited the weakest transforming activity. Cell lines transformed by E5, E6, and E7 ORFs were established and they grew anchorage-independently. The presence of HPV 59 ORF DNA was confirmed by polymerase chain reaction (PCR) and Southern blot analysis in HPV 59 ORFs-transformed cell lines.

Human papillomavirus type 70 genome cloned from overlapping PCR products: complete nucleotide sequence and genomic organization

The genome of human papillomavirus (HPV) type 70 (HPV 70), isolated from a cervical condyloma, was obtained by cloning overlapping PCR products. By automated DNA sequence analysis, the genome was found to consist of 7,905 bp with a G + C content of 40%. The genomic organization showed the characteristic features shared by other sequenced HPVs. Nucleotide sequence comparison with previously known HPV types demonstrated the closest homology with HPV 68 (82%), HPV 39 (82%), HPV 18 (70%), HPV 45 (70%), and HPV 59 (70%). Comparison with seven other partially sequenced HPV 70 isolates showed homologies of between 100 and 99.5%. Cloning of overlapping PCR products and automated DNA sequence analysis was found to be a feasible method of obtaining full-length sequences of HPVs.

Two novel genital human papillomavirus (HPV) types, HPV68 and HPV70, related to the potentially oncogenic HPV39

The genomes of two novel human papillomavirus (HPV) types, HPV68 and HPV70, were cloned from a low-grade cervical intraepithelial neoplasia and a vulvar papilloma, respectively, and partially sequenced. Both types are related to HPV39, a potentially oncogenic virus. HPV68 and HPV70 were also detected in genital intraepithelial neoplasia from three patients and one patient, respectively. Comparison with sequence data in the literature indicates that the subgenomic ME180-HPV DNA fragment, cloned from a carcinoma cell line, corresponds to an HPV68 subtype and that several HPV DNA fragments amplified by PCR from genital neoplasia represent worldwide distributed variants of HPV68 and HPV70.

Evolution of the feline-subgroup parvoviruses and the control of canine host range in vivo

A related group of parvoviruses infects members of many different carnivore families. Some of those viruses differ in host range or antigenic properties, but the true relationships are poorly understood. We examined 24 VP1/VP2 and 8 NS1 gene sequences from various parvovirus isolates to determine the phylogenetic relationships between viruses isolated from cats, dogs, Asiatic raccoon dogs, mink, raccoons, and foxes. There were about 1.3% pairwise sequence differences between the VP1/VP2 genes of viruses collected up to four decades apart. Viruses from cats, mink, foxes, and raccoons were not distinguished from each other phylogenetically, but the canine or Asiatic raccoon dog isolates formed a distinct clade. Characteristic antigenic, tissue culture host range, and other properties of the canine isolates have previously been shown to be determined by differences in the VP1/VP2 gene, and we show here that there are at least 10 nucleotide sequence differences which distinguish all canine isolates from any other virus. The VP1/VP2 gene sequences grouped roughly according to the time of virus isolation, and there were similar rates of sequence divergence among the canine isolates and those from the other species. A smaller number of differences were present in the NS1 gene sequences, but a similar phylogeny was revealed. Inoculation of mutants of a feline virus isolate into dogs showed that three or four CPV-specific differences in the VP1/VP2 gene controlled the in vivo canine host range.