Search tips
Search criteria 


Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
Arch Virol. Author manuscript; available in PMC 2010 April 12.
Published in final edited form as:
PMCID: PMC2852682

Lack of heterogeneity of HPV16 E7 sequence compared with HPV31 and HPV73 may be related to its unique carcinogenic properties


To assess the role of human papillomavirus virus (HPV) genetics in cervical lesions, we sequenced the E7 gene of HPV16, 31, or 73 from singly infected women who (1) cleared the infection quickly, (2) had type-specific persistent infection, or (3) progressed to CIN2 or worse lesions. Four of the 296 HPV16 E7 nucleotides were variable, compared with 7 of 296 for HPV31 E7 and 4 of 296 for HPV73 E7. While most of the polymorphisms in HPV31 and -73 resulted in non-synonymous amino acid changes, the polymorphisms in the HPV16 E7 resulted in synonymous changes. The lack of heterogeneity of HPV16 E7 suggests high evolutionary purifying selection that might be related to the unique carcinogenicity of HPV16.


Persistent infections with carcinogenic human papillomavirus (HPV) in a high-risk clade of the alpha genus cause virtually all cervical cancers [1, 2]. There are 10–15 carcinogenic HPV genotypes; however, there are substantial differences even among closely related HPV types in their ability to induce malignant transformation.

We are studying how HPV genetic variability contributes to HPV natural history, to gain insight into carcinogenic mechanisms using a biological system that is simpler than human genomic variability. This brief communication examines one of HPV’s early genes, E7, which is an important oncoprotein because it abrogates the function of pRB and is persistently expressed in HPV-induced cancers [3].

In particular, we compared HPV16, 31, and 73, which are genetically closely related types that manifest vast differences in their natural histories and risks for cervical cancer. HPV 16 (alpha-9 species) is by far the most carcinogenic HPV type and the main target of HPV vaccines. It is also the most prevalent carcinogenic type and most likely to persist. It causes half of cervical cancers worldwide and an even higher fraction of extra-cervical HPV-induced cancers (e.g., anal and oropharyngeal) [2]. By comparison, HPV31 (another type in alpha-9 closely related to HPV16) [4] is less prevalent and much less carcinogenic. As an even more striking contrast, HPV73 (alpha-11, a closely related species), while relatively common, is only equivocally carcinogenic [1]. The role of viral genetics in producing these pathologic differences is largely unexplored.

Previous research has shown that the E7 gene of HPV16 is relatively stable, with few nucleotide changes [58], though there is not a lot of research comparing E7 of different HPV genotypes to each other. Our objective in this study was to identify genetic sequence differences in the E7 gene of HPV16, 31 and 73 that could explain the differences in natural history between types.


Study design

This study was a cross-sectional study nested within the Guanacaste Cohort, which is a large population-based study of HPV natural history in a cohort of ~10,000 women aged 18–99 in Costa Rica [9, 10]. Cervical cells collected as part of the parent study (Guanacaste Cohort) were tested for the presence of the HPV genotype at each study visit. In addition, at each visit, a pelvic exam was conducted, and women with suspicious lesions on cytology were further evaluated by colposcopy and biopsy when indicated. We randomly selected 60 women with single HPV16 infection, 54 women with single HPV31 infection, and 20 women with single HPV73 infection who either (1) cleared the infection (n = 20 for HPV16, n = 12 for HPV31, and n = 16 for HPV73) within an average period of less than 1 year, (2) had persistent infection, i.e., the same HPV present at two consecutive visits at least 1 year apart (n = 20 for HPV16; n = 22 for HPV31; and n = 4 for HPV73), or (3) progressed to cervical intra-epithelial neoplasia grade 2 or more severe lesions (CIN2+) (n = 20 of HPV16; n = 13 for HPV31; none of the HPV73 progressed to CIN3+ in the Guanacaste Cohort).

HPV detection

Cervicovaginal samples obtained from participants were stored in standard transport medium until used for polymerase chain reaction based HPV DNA testing using the L1 MY09/11 consensus primers. Dot blot analysis was used for HPV typing as described previously [1, 9].

E7 sequencing

A nested PCR was used to amplify the E7 open reading frame from the original clinical samples using type-specific primers (primer sequences available upon request). PCR-positive samples were sequenced using the inner PCR primer set used for the final round of amplification. The sequences were analyzed using Chromas 2 and compiled using SeqMan within the DNASTAR package. ClustalW as implemented within the BioEdit software package [11] was used to align the sequences to their respective prototypes downloaded from GenBank (HPV16 NC_001526; HPV31 NC_001527; HPV73 NC_006165). Contingency tables showing the association of each viral E7 single nucleotide polymorphism (SNP) by case status are presented.


Table 1 shows a list of the samples available for this project. Samples from eight of the HPV16 samples, and two of the HPV73 samples did not yield amplification products.

Table 1
Study population by case status

Tables 2, ,3,3, and and44 contain data on the nucleotide and amino acid changes observed in the E7 gene of HPV16, 31, and 73, respectively. Four nucleotide base changes (out of 296, 1.4%) were observed in HPV16 E7 compared with 7 changes in HPV31 E7 (7/296, 2.4%) and 4 (4/296, 1.4%) in HPV73 E7. These slight numerical differences are made striking by the following observation: all 4 base pair changes in HPV16 resulted in synonymous changes in the amino acid sequence of HPV16 E7, while most of the base pair changes in HPV31 and 73 resulted in non-synonymous variations, as shown in Tables 3 and and4.4. Thus, at the amino acid level, HPV16 was less variable than HPV31 or HPV73.

Table 2
Nucleotide and amino acid sequence variations in HPV16 E7
Table 3
Nucleotide and amino acid sequence variations in HPV31 E7
Table 4
Nucleotide and amino acid sequence variations in HPV73 E7

There was an insertion of six nucleotide base pairs in one instance in HPV31 E7, which was not observed in HPV16 or -73 E7.

Table 5 shows the results of the viral single nucleotide polymorphisms by case status. We did not observe any effect by case status, although the study was not able to detect differences less than an odds ratio of 5.

Table 5
Mutations by case status


We observed a highly conserved E7 region for HPV16, but not for HPV31 and 73, especially at the level of amino acids, in agreement with previous research [58]. The few nucleotide variations in HPV16 E7 resulted in synonymous amino acid sequence changes, while most of the variations in HPV31 and 73 resulted in non-synonymous changes. Although we did not find the previously described HPV16 E7 nucleotide change in position 647 [58], three of the nucleotide changes we report were also reported by others and also resulted in synonymous amino acid changes [58]. We present two explanations for the lack of non-synonymous changes in HPV16 E7 in this study. First, the HPV16 E7 sequence represents the optimal solution for its function, and any mutants get outcompeted by the prototype. Second, the protein fold of E7 does not allow for mutations, suggesting that HPV16 E7 is at a position in the “fold space” where a single amino acid change unfolds the protein. Our study was not designed to differentiate between these possibilities. It is also possible that our findings were due to chance and hence should be corroborated in other settings.

This study had several unique characteristics. The multidisciplinary approach used coupled laboratory data with a rigorous epidemiologic study of a population-based cohort that is unique for its long and intensive follow-up. Even in this rigorous epidemiologic study, with careful selection of cases, we were unable to detect any associations in relation to case–control status, underscoring the need for further research in this area. This research complements the research community’s ongoing efforts to identify HPV cofactors needed for malignant transformation by studying the role of viral genetic polymorphisms in cancer etiology and contributes to the overall understanding of HPV E7 protein interaction with host cellular factors in cervical cancer development. In summary, the lack of heterogeneity of HPV16 E7 suggests evolutionary purifying selection, which may be related to the unique carcinogenic properties of HPV16. Moreover, the E7 variations are unlikely to explain differences in pathologic potential in this study.

Contributor Information

Mahboobeh Safaeian, Division of Cancer Epidemiology and Genetics, National Cancer Institute, Bethesda, USA ; vog.hin.liam@mnaieafas.

Koenraad van Doorslaer, Albert Einstein College of Medicine, Bronx, USA.

Mark Schiffman, Division of Cancer Epidemiology and Genetics, National Cancer Institute, Bethesda, USA.

Zigui Chen, Albert Einstein College of Medicine, Bronx, USA.

Ana Cecilia Rodriguez, Proyecto Epidemiológico Guanacaste, Fundación INCIENSA, San José, Costa Rica.

Rolando Herrero, Proyecto Epidemiológico Guanacaste, Fundación INCIENSA, San José, Costa Rica.

Allan Hildesheim, Division of Cancer Epidemiology and Genetics, National Cancer Institute, Bethesda, USA.

Robert D. Burk, Albert Einstein College of Medicine, Bronx, USA.


1. Schiffman M, Herrero R, Desalle R, Hildesheim A, Wacholder S, Rodriguez AC, Bratti MC, Sherman ME, Morales J, Guillen D, Alfaro M, Hutchinson M, et al. The carcinogenicity of human papillomavirus types reflects viral evolution. Virology. 2005;337(1):76–84. [PubMed]
2. Bouvard V, Baan R, Straif K, Grosse Y, Secretan B, El GF, Benbrahim-Tallaa L, Guha N, Freeman C, Galichet L, Cogliano V. A review of human carcinogens—part B: biological agents. Lancet Oncol. 2009;10(4):321–322. [PubMed]
3. Stanley MA. Human papillomavirus and cervical carcinogenesis. Best Pract Res Clin Obstet Gynaecol. 2001;15(5):663–676. [PubMed]
4. de Villiers EM, Fauquet C, Broker TR, Bernard HU, zur HH. Classification of papillomaviruses. Virology. 2004;324(1):17–27. [PubMed]
5. Pillai MR, Hariharan R, Babu JM, Lakshmi S, Chiplunkar SV, Patkar M, Tongaonkar H, Dinshaw K, Jayshree RS, Reddy BK, Siddiqui M, Roychoudury S, et al. Molecular variants of HPV-16 associated with cervical cancer in Indian population. Int J Cancer. 2009;125(1):91–103. [PubMed]
6. Fujinaga Y, Okazawa K, Nishikawa A, Yamakawa Y, Fukushima M, Kato I, Fujinaga K. Sequence variation of human papillomavirus type 16 E7 in preinvasive and invasive cervical neoplasias. Virus Genes. 1994;9(1):85–92. [PubMed]
7. Song YS, Kee SH, Kim JW, Park NH, Kang SB, Chang WH, Lee HP. Major sequence variants in E7 gene of human papillomavirus type 16 from cervical cancerous and noncancerous lesions of Korean women. Gynecol Oncol. 1997;66(2):275–281. [PubMed]
8. Eschle D, Durst M, ter MJ, Luande J, Eberhardt HC, Pawlita M, Gissmann L. Geographical dependence of sequence variation in the E7 gene of human papillomavirus type 16. J Gen Virol. 1992;73(Pt 7) 32-1829. [PubMed]
9. Herrero R, Schiffman MH, Bratti C, Hildesheim A, Balmaceda I, Sherman ME, Greenberg M, Cardenas F, Gomez V, Helgesen K, Morales J, Hutchinson M, et al. Design and methods of a population-based natural history study of cervical neoplasia in a rural province of Costa Rica: the Guanacaste Project. Rev Panam Salud Publica. 1997;1(5):362–375. [PubMed]
10. Herrero R, Hildesheim A, Bratti C, Sherman ME, Hutchinson M, Morales J, Balmaceda I, Greenberg MD, Alfaro M, Burk RD, Wacholder S, Plummer M, et al. Population-based study of human papillomavirus infection and cervical neoplasia in rural Costa Rica. J Natl Cancer Inst. 2000;92(6):464–474. [PubMed]
11. Hall TA. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. 1999;41:95–98.