The genetic drift of human papillomavirus type 16 is a means of reconstructing prehistoric viral spread and the movement of ancient human populations.

We have investigated the diversity of a hypervariable segment of the human papillomavirus type 16 (HPV-16) genome among 301 virus isolates that were collected from 25 different ethnic groups and geographic locations. Altogether, we distinguished 48 different variants that had diversified from one another along five phylogenetic branches. Variants from two of these branches were nearly completely confined to Africa. Variants from a third branch were the only variants identified in Europeans but occurred at lower frequency in all other ethnic groups. A fourth branch was specific for Japanese and Chinese isolates. A small fraction of all isolates from Asia and from indigenous as well as immigrant populations in the Americas formed a fifth branch. Important patterns of HPV-16 phylogeny suggested coevolution of the virus with people of the three major human races, namely, Africans, Caucasians, and East Asians. But several minor patterns are indicative of smaller bottlenecks of viral evolution and spread, which may correlate with the migration of ethnic groups in prehistoric times. The colonization of the Americas by Europeans and Africans is reflected in the composition of their HPV-16 variants. We discuss arguments that today's HPV-16 genomes represent a degree of diversity that evolved over a large time span, probably exceeding 200,000 years, from a precursor genome that may have originated in Africa. The identification of molecular variants is a powerful epidemiological and phylogenetic tool for revealing the ancient spread of papillomaviruses, whose trace through the world has not yet been completely lost.     

Evolution of human papillomavirus type 18: an ancient phylogenetic root in Africa and intratype diversity reflect coevolution with human ethnic groups.

Papillomaviruses are an ideal model system for the study of DNA virus evolution. On several levels, phylogenetic trees of papillomaviruses reflect the relationship of their hosts. Papillomaviruses isolated from remotely related vertebrates form major branches. One branch of human papillomaviruses (HPVs) includes an ape and two monkey papillomaviruses, possibly because the diversification of the viruses predated the separation of the infected-primate taxa. This hypothesis predicts that the root of the evolution of some if not all HPV types should point to Africa, since humans evolved from nonhuman primates in this continent. We tested this hypothesis and compared the genomic sequences of HPV type 18 (HPV-18) isolates from four continents. Diversity within HPV-18 correlates with patterns of the evolution and spread of Homo sapiens: HPV-18 variants, just like HPV-16 variants, are specific for the major human races, with maximal diversity in Africa. Outgroup rooting of the HPV-18 tree against HPV-45, which is closely related to HPV-18, identifies African HPV-18 variants at the root of the tree. The identification of an African HPV-45 isolate further reduces the evolutionary distance between HPV-18 and HPV-45. HPV-18 variants from Amazonian Indians are the closest relatives to those from Japanese and Chinese patients and suggest that a single point mutation in the phylogenetically evaluated genomic segment represents at least 12,000 years of evolution. We estimate that diversity within HPV-18 and probably within other HPV types evolved over a period of more than 200,000 years and that diversity between HPV types evolved over several million years.     

Prevalence of human papillomavirus type 16 variants in the Federal District, Central Brazil

We report the prevalence of human papillomavirus type 16 (HPV-16) variants in women with cervical lesions from the Federal District, Central Brazil. We analyzed 34 HPV-16 samples, identifying the sequence variations of E6 and L1 genes and correlating variant frequency with disease status. The most prevalent HPV-16 variant was the European (50%), followed by Asian-American (41.2%), African-1 (5.9%), and African-2 (2.9%). European and non-European variants appeared in equal frequencies among the cytological types of lesions - atypical squamous or glandular cells of undetermined significance, cytological alterations suggesting HPV infection, cervical intraepithelial neoplasias, squamous cell carcinoma, and adenocarcinoma.     

Intranasal vaccination with recombinant adeno-associated virus type 5 against human papillomavirus type 16 L1

Adeno-associated viruses (AAV) have been developed and evaluated as recombinant vectors for gene therapy in many preclinical studies, as well as in clinical trials. However, only a few approaches have used recombinant AAV (rAAV) to deliver vaccine antigens. We generated an rAAV encoding the major capsid protein L1 (L1h) from the human papillomavirus type 16 (HPV16), aiming to develop a prophylactic vaccine against HPV16 infections, which are the major cause of cervical cancer in women worldwide. A single dose of rAAV5 L1h administered intranasally was sufficient to induce high titers of L1-specific serum antibodies, as well as mucosal antibodies in vaginal washes. Seroconversion was maintained for at least 1 year. In addition, a cellular immune response was still detectable 60 weeks after immunization. Furthermore, lyophilized rAAV5 L1h successfully evoked a systemic and mucosal immune response in mice. These data clearly show the efficacy of a single-dose intranasal immunization against HPV16 based on the recombinant rAAV5L1h vector without the need of an adjuvant.     

Molecular characterization of biologically diverse envelope variants of human immunodeficiency virus type 1 derived from an individual.

The envelope genes of six viruses derived from a single sampling from an individual chronically infected with human immunodeficiency virus type 1 (RJS-4) have been analyzed. Here we present the nucleotide and predicted amino acid sequences of these variants and show a correlation between biological properties and disturbance of the envelope reading frame.     

Identification of a transforming gene of human papillomavirus type 16.

Previously, we observed sequential two-step alteration, growth stimulation, and progression to a more malignant state in NIH 3T3 cells transfected by human papillomavirus type 16 (HPV-16) DNA. In this study, we prepared a cDNA library from RNA extracted from cells transfected with the HPV-16 DNA and isolated cDNA clones which had growth-stimulating activity. Analysis of these cDNA clones indicated that the E7 open reading frame alone is responsible for inducing both steps of this cell transformation.     

Characterization of a human papillomavirus type 16 variant-dependent neutralizing epitope.

We have determined that three type-specific and conformationally dependent monoclonal antibodies, H16.E70, H16.U4, and H16.V5, neutralize pseudotype human papillomavirus type 16 (HPV16) virions in vitro. H16.U4 and H16.V5 neutralized pseudotype virions derived from the German HPV16 variant 114K and the Zairian variant Z-1194 with equal efficiency. In contrast, neutralization of Z-1194 pseudotype virions by H16.E70 was two orders of magnitude weaker than neutralization of 114K pseudotype virions. This difference correlated with enzyme-linked immunosorbent assay reactivity of H16.E70 to L1 virus-like particles of the two variants. A substitution at residue 282 of L1 was responsible for this differential reactivity, suggesting that this residue constitutes part of the H16.E70 epitope.     

Mutational analysis of human papillomavirus type 16 E7 functions.

The human papillomavirus type 16 E7 gene encodes a nuclear oncoprotein (98 amino acids [AAs] long) consisting of three regions: regions 1 (AAs 1 to 20) and 2 (AAs 21 to 40), which show high homology to the sequences of conserved domains 1 and 2, respectively, of adenovirus E1A; and region 3 (AAs 41 to 98) containing two metal-binding motifs Cys-X-X-Cys (AAs 58 and 91 to 94). We constructed AA deletion (substitution) mutants and single-AA substitution mutants of E7 placed under the control of the simian virus 40 promoter and examined their biological functions. Stable expression of E7 protein in monkey COS-1 cells required almost the entire length of E7 and was markedly lowered by the mutations in region 3. Transactivation of the adenovirus E2 promoter in monkey CV-1 cells was lowered by the mutations. It was abolished by changing Cys-24 to Gly and markedly decreased by a mutation at His-2 or at the metal-binding motifs in region 3. Focal transformation of rat 3Y1 cells by E7 was eliminated by changing His-2 to Asp or Cys-24 to Gly and was greatly impaired by changing Cys-61 or Cys-94 to Gly. The transforming function survived mutations at Leu-13 and Cys-68 and deletion of Asp-Ser-Ser (AAs 30 to 32). The data suggest that regions 1 to 3 are required for its functions and that the meta-binding motifs in region 3 are required to maintain a stable or functional structure of the E7 protein.     

1917 avian influenza virus sequences suggest that the 1918 pandemic virus did not acquire its hemagglutinin directly from birds

Wild waterfowl captured between 1915 and 1919 were tested for influenza A virus RNA. One bird, captured in 1917, was infected with a virus of the same hemagglutinin (HA) subtype as that of the 1918 pandemic virus. The 1917 HA is more closely related to that of modern avian viruses than it is to that of the pandemic virus, suggesting (i) that there was little drift in avian sequences over the past 85 years and (ii) that the 1918 pandemic virus did not acquire its HA directly from a bird.     

Prevalence of human papillomavirus (HPV) type 16 variants and rare HPV types in the central Amazon region

Infection by human papillomavirus (HPV) is one of the primary causes of mortality by cancer in northern Brazil. Sexually active women from Manaus, Amazonas, without cytological alterations and women with pre-malignant and malignant cytological alterations were examined for HPV virus, identified via PCR and sequencing. The target region for this study was part of the L1 capsid gene of HPV. Twenty-three samples that were PCR-positive were sequenced. Analysis of 336 bp demonstrated a high incidence of high-risk HPV types in the population of Manaus, identified as HPVs 16, 33, 58, 66, 68. HPV type 16 was the most prevalent, presenting two variants similar to the Asian-American (AA) and East-Asian type (As) variants. A rare HPV type 13 related to "Heck's disease" was also detected. This preliminary provides important information about the HPV circulating in Amazonas State.     

Tumorigenic transformation of murine keratinocytes by the E5 genes of bovine papillomavirus type 1 and human papillomavirus type 16.

To examine the biological properties of the bovine papillomavirus type 1 (BPV) and human papillomavirus type 16 (HPV16) E5 genes, each was cloned separately into a retroviral expression vector and helper-free recombinant viruses were generated in packaging cell lines. The BPV E5 retroviruses efficiently caused morphologic and tumorigenic transformation of cultured lines of murine fibroblasts, whereas the HPV16 E5 viruses were inactive in these assays. In contrast, infection of the p117 established line of murine epidermal keratinocytes with either the BPV or the HPV16 E5 retrovirus resulted in the generation of tumorigenic cells. Pam212 murine keratinocytes were also transformed to tumorigenicity by the HPV16 E5 gene but not by the gene carrying a frameshift mutation. These results establish that the HPV16 E5 gene is a transforming gene in cells related to its normal host epithelial cells.     

Modelling of the human papillomavirus type 16 E5 protein

Cytochrome (cyt) c forms complexes, undergoes a conformational change and becomes partly reduced at interaction with membrane anchored alkaline phosphatase (AP), a glycoprotein which is released into the body fluid in forms differing in hydrophobicity. The proportion of products formed in the mixtures depends on pH, ionic strength, temperature and the buffer composition. The reaction terminates in an equilibrium between cyt c(FeII) and other cyt c conformers. Optimal conditions for the rate of the reaction are 100 mM glycine/NaOH, pH 9.7-9.9, at which 68-74% of cyt c is found in the reduced state. The interaction affects compactness of the haem cleft as shown by changes induced in CD spectra of the Soret region and changes in optical characteristics of phenylalanine, tyrosine and tryptophan residues. Differential scanning calorimetry of AP+cyt c mixtures revealed a creation of at least two types of complexes. A complex formed by non-coulombic binding prevails at substoichiometric AP/cyt c ratios, at higher ratios more electrostatic attraction is involved and at 1:1 molar ratio an apparent complexity of binding forces occurs. The rapid phase of the cyt c(FeII) formation depends on the presence of the hydrophobic alkylacylphosphoinositol (glycosylphosphatidylinositol) moiety, the protein part of the enzyme participates in an electrostatic and much slower phase of cyt c(FeII) creation. The results show that non-coulombic interaction may participate at interaction of cyt c with cellular proteins.     

Worldwide genomic diversity of the high-risk human papillomavirus types 31, 35, 52, and 58, four close relatives of human papillomavirus type 16

Among the more than one hundred formally described human papillomavirus (HPV) types, 18 are referred to as high-risk HPV types due to their association with anogenital cancer. Despite pathogenic similarities, these types form three remotely related taxonomic groups. One of these groups is called HPV species 9 and is formed by HPV-16, the most common and best-studied type, together with HPV-31, -33, -35, -52, -58, and -67. Previous worldwide comparisons of HPV-16 samples showed about 2% nucleotide diversity between isolates, which were subsequently termed variants. The distribution of divergent variants has been found to correlate frequently with the geographic origin and the ethnicity of the infected patients and led to the concept of unique African, European, Asian, and Native American HPV-16 variants. In the current study, we address the question of whether geography and ethnicity also correlate with sequence variations found for HPV-31, -35, -52, and -58. This was done by sequencing the long control region in samples derived from Europe, Asia, and Africa, and from immigrant populations in North and South America. We observed maximal divergence between any two variants within each of these four HPV types ranging from 1.8 to 3.6% based on nucleotide exchanges and, occasionally, on insertions and deletions. Similar to the case with HPV-16, these mutations are not random but indicate a relationship between the variants in form of phylogenetic trees. An interesting example is presented by a 16-bp insert in select variants of HPV-35, which appears to have given rise to additional variants by nucleotide exchanges within the insert. All trees showed distinct phylogenetic topologies, ranging from dichotomic branching in the case of HPV-31 to star phylogenies of the other three types. No clear similarities between these types or between these types and HPV-16 exist. While variant branches in some types were specific for Europe, Africa, or East Asia, none of the four trees reflected human evolution and spread to the extent illustrated by HPV-16. One possible explanation is that the rare HPV types that we studied spread and thereby diversified more slowly than the more abundant HPV-16 and may have established much of today's variant diversity already before the worldwide spread of humans 100,000 years ago. Most variants had prototypic amino acid sequences within the E6 oncoprotein and a segment of the L1 capsid protein. Some had one, two, or three amino acid substitutions in these regions, which might indicate biological and pathogenic diversity between the variants of each HPV type.     

Interaction of the human papillomavirus type 16 E6 oncoprotein with wild-type and mutant human p53 proteins.

The E6 oncoproteins encoded by the cancer-associated human papillomaviruses (HPVs) can associate with and promote the degradation of wild-type p53 in vitro. To gain further insight into this process, the ability of HPV-16 E6 to complex with and promote the degradation of mutant forms of p53 was studied. A correlation between binding and the targeted degradation of p53 was established. Mutant p53 proteins that bound HPV-16 E6 were targeted for degradation, whereas those that did not complex HPV-16 E6 were not degraded. Since the HPV-16 E6-promoted degradation involves the ubiquitin-dependent proteolysis pathway, specific mutations were made in the amino terminus of p53 to examine whether the E6 targeted degradation involved the N-end rule pathway. No requirement for destabilizing amino acids at the N terminus of p53 was found, nor was evidence found that HPV-16 E6 could provide this determinant in trans, indicating that the N-terminal rule pathway is not involved in the E6-promoted degradation of p53.     

Effect of human papillomavirus type 16 oncogenes on MAP kinase activity.

The mitogen-activated protein (MAP) kinase signal transduction pathway is an intracellular signaling cascade which mediates cellular responses to growth and differentiation factors. The MAP kinase pathway can be activated by a wide range of stimuli dependent on the cell types, and this is normally a transient response. Oncogenes such as ras, src, raf, and mos have been proposed to transform cells in part by prolonging the activated stage of components within this signaling pathway. The human papillomavirus (HPV) oncogenes E6 and E7 play an essential role in the in vitro transformation of primary human keratinocytes and rodent cells. The HPV type 16 E5 gene has also been shown to have weak transforming activity and may enhance the epidermal growth factor (EGF)-mediated signal transduction to the nucleus. In the present study, we have investigated the effects of the oncogenic HPV type 16 E5, E6, and E7 genes on the induction of the MAP kinase signaling pathway. The E5 gene induced an increase in the MAP kinase activity both in the absence and in the presence of EGF. In comparison, the E6 and E7 oncoproteins do not alter the MAP kinase activity or prolong the MAP kinase activity induced with EGF. These findings suggest that E5 may function, at least in part, to enhance the cell response through the MAP kinase pathway. However, the transforming activity of E6 and E7 is not associated with alterations in the MAP kinase pathway. These findings are consistent with E5 enhancing the response to growth factor stimulation.     

Identification of the E5 open reading frame of human papillomavirus type 16.

Sequencing of the E5 open reading frame (ORF) of human papillomavirus type 16 revealed an additional nucleotide, a thymidine residue, at position 3903 compared with the original sequence (Seedorf et al., Virology 145:181-185, 1985). The additional T had two effects; first, in reading frame 2, in which the original E5 ORF was predicted, the additional T changed the reading frame downstream of position 3903 to create an ORF, which we designated E5, that terminated at position 4018 and potentially encoded a 52-amino-acid polypeptide. Secondly, in reading frame 3, a new ORF was created (positions 3807 to 4097), which we propose is the authentic papillomavirus type 16 E5 ORF. It contained a methionine residue and encoded an additional 82 amino acids. Both ORFs have been cloned into bacterial expression vectors (pATH), and the fusion proteins have been used to generate polyclonal antibodies in rabbits.