Previous Article | Next Article 
Journal of Clinical Microbiology, March 2008, p. 1060-1066, Vol. 46, No. 3
0095-1137/08/$08.00+0 doi:10.1128/JCM.02202-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.
Human Papillomavirus Type 16 Variant Analysis of E6, E7, and L1 Genes and Long Control Region in Biopsy Samples from Cervical Cancer Patients in North India
Shailja Pande,1
Neeraj Jain,1
Bhupesh K. Prusty,1
Suresh Bhambhani,2
Sanjay Gupta,2
Rajyashri Sharma,3
Swaraj Batra,4 and
Bhudev C. Das1*
Division of Molecular Oncology,1
Division of Cytopathology, Institute of Cytology and Preventive Oncology (ICMR), I-7, Sector 39, Noida 201301, India,2
Department of Obstetrics and Gynecology, JNMCH, Aligarh, India,3
Department of Obstetrics and Gynecology, LNJP, New Delhi, India4
Received 14 November 2007/
Accepted 2 January 2008

ABSTRACT
High-risk human papillomaviruses (HPVs), particularly HPV types
16 and 18 (HPV-16 and HPV-18, respectively), play a cardinal
role in the etiology of cervical cancer. The most prevalent
type, HPV-16, shows intratypic sequence variants that are known
to differ in oncogenic potential and geographic distribution.
This study was designed to analyze sequence variations in E6,
E7, and L1 genes and the LCR (for long control region) of HPV-16
in cervical cancer patients to identify the most prevalent and
novel HPV-16 variants and to correlate them with the severity
of the disease. Cervical biopsies from 60 HPV-16-positive cancer
cases were analyzed by PCR and DNA sequencing. The most frequently
observed variations were T350G (100%) in E6, T789C (87.5%) in
E7, A6695C (54.5%) in L1, and G7521A (91.1%) in the LCR. In
addition, only one novel variant (T527A) in E6 and four new
variants each in L1 (A6667C, A6691G, C6906T, and A6924C) and
in the LCR (C13T, A7636C, C7678T, and G7799A) were identified.
While E7 was found to be highly conserved, the variant 350G
of E6 was the most prevalent in all of the histopathological
grades. The majority of LCR variants were found at the YY1 transcription
factor binding sites. Interestingly, a complete absence of the
Asian lineage and a high prevalence of European lineages in
E6, E7, L1, and the LCR (85, 86.7, 67.7, and 63.3%, respectively)
indicate a possible epidemiological linkage between Europe and
India with regard to the dissemination of HPV-16 infections
in India.

INTRODUCTION
Cervical cancer is the major cancer in Indian women and a leading
cause of cancer deaths. Every year, more than 130,000 new cases
and about 70,000 deaths are recorded. The persistent infection
by specific types of high-risk human papillomaviruses (HR-HPVs)
is essential for the progression of cervical lesions (
6,
11,
31,
60), and women who are infected with HR-HPVs are likely
to develop cancer (
3,
4,
27,
40). Various studies have demonstrated
that more than 70% of invasive cervical cancers harbor HPV type
16 (HPV-16) and HPV-18 (
17,
31), and the products of viral transforming
E6 and E7 genes have been shown to contribute to tumorigenesis
by functionally inactivating two important cellular tumor suppressor
proteins, p53 and retinoblastoma (
5,
14,
30,
55).
More than 100 types of HPV are known, but only about 30 types are associated with anogenital cancer. According to the Papillomavirus Nomenclature Committee, a new HPV type is defined by a nucleotide sequence variation of more than 10% compared to that of other known HPV types in the E6, E7, and L1 open reading frames. Those differing by 2 to 10% are referred to as subtypes, whereas intratype variants may vary by up to 2% in the coding region and 5% in the noncoding region compared to that of the prototype (2, 10).
On the basis of sequence variations in E6, L1, L2, and the long control region (LCR), HPV-16 variants have been identified and grouped into six distinct phylogenetic branches: E (European), AA (Asian-American), Af1 (African 1), Af2 (African 2), As (Asian), and NA (North American) (54, 56, 57). These variants have been found to show different geographic distributions, with various oncogenic potentials. A number of sequence variations have been reported for HPV-16 E6, E7, and L1 genes as well as in the LCR in cervical cancer (33, 39, 48, 55, 56). Studies also have shown that specific intratype variants may influence the persistence of HPV infection and the progression of precursor lesions to cancer (27, 58, 59). HPV variants also may affect virus assembly, immunologic responses, pathogenicity, p53 degradation, immortalization activity, and the regulation of transcription (15, 18, 24, 25, 37, 51). These variations immediately affect the sensitivity and specificity of different PCR-based genotype diagnostic methods.
Of the two HPV-16 oncogenes E6 and E7, E6 has been found to show more variations than E7, which is relatively conserved (48, 50, 56, 58, 59). The analysis of the L1 gene, which codes for the viral major capsid protein, is of immense importance because of its high diagnostic value. The range of intratype variation observed in this region allows the distinction and assessment of known and novel HPV types (46). This is also an important target for the development of HPV vaccines. Very few reports are available on HPV variants from this region (9, 56). The LCR, which regulates viral transcription, has been found to be the most variable region of the HPV-16 genome (41, 53, 56). Several authors have evaluated the association of specific HPV-16 variants with viral persistence and the development of cervical intraepithelial neoplasia lesions (3, 27, 34, 50, 52).
In India, as high as 98% of the cervical carcinoma cases are found to harbor HPV infection, and the most prevalent (
80%) type is HPV-16 (7, 8, 22). Although the prevalence of HPV and cervical cancer in India is the highest in the world, there is not much information on HPV variants from different regions of the Indian subcontinent. In the present study, we have examined the sequence variations in E6, E7, and L1 genes and the LCR of the most prevalent HR-HPV type, HPV-16, and correlated them with the age of the women, histopathologic grades of tumors, and oncogenic potential.

MATERIALS AND METHODS
Study population and specimen collection.
Out of a total of 90 tumor samples screened by standard procedures
(
20,
38), 60 HPV-16-positive cervical cancer cases were recruited
for the variant analysis in the present study. The tumor samples
were from biopsies collected for the routine diagnosis of cervical
cancer from the Cancer Clinic of the Department of Gynecology
and Obstetrics, Lok Nayak Hospital, New Delhi, India, and Jawaharlal
Nehru Medical College, Aligarh, India. Fresh tumor biopsies
were collected in chilled phosphate-buffered saline (PBS) and
stored at –70°C in a deep freezer for further processing
later. Informed consent was obtained from all patients, and
the study was approved by the institutional ethical committee.
The diagnosis of histopathological grades was done independently
by two pathologists. Discrepancies were resolved by a third
pathologist to arrive at a final decision.
DNA extraction and typing of HPV.
High-molecular-weight genomic DNA from cervical biopsies was isolated by the standard method of proteinase K digestion and phenol-chloroform extraction that is routinely followed in our laboratory (8, 20).
The detection of HPV and the typing of HR-HPV-16 and HR-HPV-18 were carried out using consensus and type-specific primers described earlier (20, 38). The amplification of the β-globin gene served as the internal control. PCR was performed using the in-house PCR protocol that is routinely followed in our laboratory (7, 8, 17). Briefly, the method involved a 25-µl reaction mix containing 100 to 200 ng DNA, 10 mM Tris-Cl (pH 8.4), 50 mM KCl, 1.5 mM MgCl2, 12.5 µM of each deoxynucleoside triphosphate (dATP, dCTP, dGTP, and dTTP), 5 pmol of each oligonucleotide primer, and 0.5 U Taq DNA polymerase (Applied Biosystems). The temperature profile used for amplification was an initial denaturation at 95°C for 5 min, followed by 30 cycles with denaturation at 95°C for 30 s, annealing at 55°C for 30 s, and extension at 72°C for 30 s, which was extended for 4 min in the final cycle. The oligonucleotide primers were synthesized in an automated Applied Biosystems DNA synthesizer (model 381A; Foster City, CA) and purified by high-pressure liquid chromatography.
Variant analysis of E6, E7, and L1 genes and the LCR of HPV-16 by PCR and direct sequencing.
HPV-16 E6-, E7-, and LCR-specific PCRs were performed with the following specific primers: for HPV-16 E6 (nucleotide [nt] 83 to 559) with an amplicon size of 476 bp, 5'-GAA ACC GGT TAG TAT AAA AGC AGA C-3' and 5'-AGC TGG GTT TCT CTA CGT GTT CT-3'; for E7 (nt 562 to 858) with an amplicon size of 296 bp, 5'-CCA TAA TAT AAG GGG TCG GTG GA-3' and 5'-TTT TTC CAC TAA CAG CCT CTA CAT-3'; for the LCR (nt 7437 to 7456 and nt 119 to 147) with an amplicon size of 617 bp, 5'-CCA TTT TGT AGC TTC AAC CCG-3' and 5'-AAG TGT GGT AAC TTT CTG GGT CGC TCC TG-3'(54). The L1 gene was amplified partially (450 bp) using MY 11/MY 09, the consensus primer, with sequences 5'-GCM CAG GGW CAT AAY AAT GG-3' and 5'-CGT CCM ARR GGA WAC TGA-3', where M = A and C, W = A and T, Y = C and T, and R = A and G (38).
For variant analysis, the PCR products were directly sequenced on an automated DNA sequencer (310 ABI Prism genetic analyzer; Applied Biosystems) according to the manufacturer's protocol. PCR products first were purified using an ammonium acetate-ethanol precipitation method to remove unused deoxynucleoside triphosphates and primers and then cycle sequenced using the BigDye Terminator sequencing ready reaction mix (Applied Biosystems) on a Gene Amp PCR 9700 (Applied Biosystems). The raw data were collected using ABI Prism 310 collection software, analyzed using sequencing analysis software (version 3.4.1) on a Macintosh operating system (version 9.1), and compared to the sequence of GenBank/EMBL/DDBJ sequence no. U89348. All of the samples were reverse sequenced to corroborate the findings.
Identification of variants.
Variants were identified using the prototype sequence (HPV-16R), which belongs to the European lineage, as the standard for comparisons and nucleotide position numbering (32), and variants were classified into the six major branches European (E), Asian-American (AA), Asian (As), African 1 (Af1), African 2 (Af2), and North American (NA), as described by Yamada et al. (56).
Statistical analysis.
The statistical analysis was performed using the
2 test to determine the association between HPV-16 variants and the age and histopathologic status of the cancer patients. Fisher's exact test was employed for comparisons of small numbers. A P value of less than 0.05 was considered statistically significant.

RESULTS
Out of 90 cervical tumor samples screened with L1 consensus
primers, 65 (72.2%) were positive for HPV, and further genotyping
with type-specific primers revealed 92.3% (
n = 60) positivity
for HPV-16 and 7.7% (
n = 5) positivity for HPV-18. All 60 patients
recruited were exclusively infected with HPV-16, and none of
them showed the presence of HPV-18. The patients were 30 to
80 years old, with a mean age of 51.3 years (standard deviation,
12.8). The patients were divided into two age groups: 30 to
50 years and more than 50 years. All of the cases were of squamous
cell carcinoma (SCC); 23 (38.3%) were well-differentiated (WDSCC),
35 (58.3%) were moderately differentiated (MDSCC), and 2 (3.4%)
were poorly differentiated squamous cell carcinomas (PDSCC).
The E6, E7, and L1 genes and the LCR were analyzed for nucleotide
variations and were compared to the prototype sequence HPV-16R
of the E lineage.
Association between HPV variants and histopathological grades.
For the E6 gene, variants were found to be more common than the prototype in both WDSCC (68.2%) and MDSCC (65.8%) cases (the MDSCC cases include two PDSCC cases), while the prototypes of E7 (86.4 and 86.8%, respectively) and L1 (60.0 and 66.7%, respectively) (Table 1) were more common than variants. The frequency of LCR variants also was found to be higher in WDSCC and MDSCC (81.8 and 71.1%, respectively) than that of the prototype (Table 1).
E6 sequence variations.
Nucleotide sequences of the complete E6 open reading frame in
all 60 cancer patients screened were compared to the HPV-16
reference sequence. Forty-two (70%) patients showed seven different
types of E6 gene mutations, of which four, G145T, C335T, T350G,
and T527A, were missense mutations: glutamine

histidine (Q14H),
histidine

tyrosine (H78Y), leucine

valine (L83V), and serine

threonine
(S142T), respectively. The remaining three, T286A, A289G, and
A532G, led to silent mutations at codons 61, 62, and 143 for
amino acids alanine, valine, and serine, respectively. The most
frequently observed mutation was the European variant T350G
(Fig.
1a), which was detected in all 42 (100%) patients either
alone (34 [81%]) or along with other variants (8 [19%]) (Fig.
2). Additional E6 gene variations were observed only in those
cases that showed the 350G variation, and no case with the 350T
sequence had any nucleotide change. The majority of cases (85%)
were found to be of E lineage, 10% belonged to branch NA, 3.3%
to AA, and none to As, Af1, or Af2 (see Table
3). Interestingly,
a novel mutation, T527A, was recorded along with 350G (Fig.
2) for the first time.
E7 sequence variations.
The E7 gene revealed only five nucleotide variations in 8 (13.3%)
patients, whereas the other 52 (86.7%) cases showed the prototype
sequence. All of the variants were silent nucleotide alterations;
one was a transversion, and the other four were transition substitutions.
The silent mutations detected were G666A, T732C, T789C, T795G,
and T828C for the amino acids glutamic acid (E), phenylalanine
(F), isoleucine (Il), threonine (T), and isoleucine (Il) at
codons 35, 57, 76, 78, and 89, respectively. E7 sequencing data
indicated that 86.7% belonged to branch E (prototype), 11.6%
to Af2, 1.7% to As, and none to NA, AA, or Af1 (see Table
3).
L1 sequence variations.
A partial sequence of the L1 gene was screened using a 450-bp amplicon for samples from 31 cervical cancer patients. A total of 13 nucleotide variations were detected in 11 (35.5%) patients, while the majority (20 [64.5%]) showed the prototype sequence. Of 13 variants, 5 (38.5%) led to missense mutations and 8 (61.5%) led to silent mutations.
Of four variants found independently (Fig. 2), two were silent mutations, A6667C and A6691G, coding for serine at codons 343 and 351, whereas A6803T and A6964C led to the missense mutations threonine
serine (T389S) and lysine
asparagine (K442N). The most frequently detected L1 variation, A6695C (Fig. 1c), was found in six (54.5%) cases. Interestingly, four novel L1 variants consisting of two silent mutations, at A6667C and A6691G, and two missense mutations, serine
phenylalanine (S423F) at C6906T and glutamine
proline (Q429P) at A6924C, were observed (Table 2). The novel L1 variants that exhibited silent mutations had only prototype sequences at E6 and E7 gene regions and the LCR (Fig. 2). The L1 sequence data showed that 67.7% belonged to branch E, 3.2% to branch E (G131), 12.9% to AA, 6.5% to NA, and none to As, Af1, or Af2; 9.7% were new variants (Table 3).
View this table:
[in this window]
[in a new window]
|
TABLE 2. Novel variants of HPV-16 E6, L1, and the LCR according to nucleotide position, variant, altered amino acid, and mutation
|
LCR sequence variations.
The LCR of HPV-16 showed the highest number of nucleotide variations
in various known and unknown regulatory binding sites of a number
of cellular and/or viral transcription/regulatory factors. A
total of 45 (75%) variant cases that were analyzed showed 17
different point mutations (nine transitions and eight transversions).
Of these, seven (G7521A, A7636C, C7689A, T7714G, C7792T, G7826A,
and A7839C) were detected in the known binding sites for various
transcription factors, such as YY1, AP1, TEF1, NF1, and Oct-1.
The most commonly observed LCR variation was the transition
substitution G7521A (Fig.
1d), which was detected in 41 (91.1%)
patients at one of the many YY1 binding sites. This mutation
occurred alone in as many as 22 (53.7%) cases (Fig.
2).
Another frequently observed LCR variation was a transversion mutation, T7714G, detected in nine (20%) cases, mostly along with G7521A, which is located in the NF1 transcription factor binding site. Two other LCR variants, T7743G and C7764T, were detected upstream of the NF1 binding site in 11.1 and 22.2% of cases, respectively. We report here four novel LCR variants, C13T, A7636C, C7678T, and G7799A (Table 2), of which three were transition mutations and one (A7636C) was a transversion mutation located at the AP1 binding site.
A phylogenetic classification of LCR variants was found to exhibit the same pattern as that observed for E6, E7, and L1, showing that a majority (63.3%) belonged to branch E (P), 13.3% to AA, 11.6% to Af1, 3.3% to Af2, and none to NA or As. Four (8.3%) novel LCR variants were observed for the first time from India (Table 3).

DISCUSSION
HPV variant data are important in developing HPV diagnostics,
vaccines, and other therapeutic approaches to control virus-induced
diseases. HPV-16 variants have been shown to have different
biological as well as biochemical effects, resulting in altered
oncogenic potentials (
1,
47,
59). The oncogenicity of distinct
HPV variants also may differ between geographical regions because
of differences in the population related to the distribution
of HLA alleles (
28).
Nucleotide sequence variations observed in E6, E7, and L1 genes and the LCR of HPV-16 in 60 cervical carcinoma cases showed a high prevalence of 350G E6 variants (the prototype is 350T) in all of the histopathologic grades (WDSCC, MDSCC, and PDSCC). This is in good agreement with previous reports (1, 59), which demonstrated an increased frequency of 350G E6 variants in cervical cancer leading to an amino acid change of leucine to valine at position 83 (L83V), and this has been associated with the progression of cervical lesions (59). In contrast, several authors failed to find a significant role of the E6 variants in cervical cancer (19, 33, 50). A study of the prevalence of the E6 gene mutations in a subset of south Indian cervical cancer patients found only 20% of cancer patients with 350G variants, which is in sharp contrast to the present study, which shows that 78% of the patients had 350G variants. Our results, however, are in good agreement with those of a study by Sathish et al. (39) from the southeastern part of India.
We have observed that nucleotide positions 83 to 144, 146 to 285, and 351 to 526 are the most conserved regions of HPV-16 E6, which seems to be important for silencing E6 expression using short interfering RNA and/or ribozyme treatment. Also, the amino acids located in these regions may play an important immunogenic role in new vaccine strategies (21). The amino acid substitutions analyzed in this study are Q14H, H78Y, and L83V, and they also were observed by Stöppler and his group (47). Most interestingly, 100% distribution (42/42) of the HPV-16 E6 E variant 350G was detected; it was found either alone (81%; 34/42) or in combination with additional E6 gene mutations. The E6 variations G145T and C335T that were detected with 350G also have been observed in AA, NA, Af1, and Af2 lineages. A novel E6 gene mutation, T527A (S142T), that replaced serine with threonine was observed in only one (2.4%) case. In sharp contrast to a report from China that showed a high prevalence of the HPV-16 As lineage (55), our results as well as the reports of other authors from India (36, 39) suggest the complete absence of the As lineage in India.
HPV-16 E7 gene mutation has been reported to be rare in various geographic and ethnic populations around the world (16, 36). However, several studies, particularly from Japan, Korea, and Indonesia, have reported a high frequency of E7 gene mutation (65 to 75%) in cervical cancer patients (9, 45). Interestingly, a recent report from China showed that 100% of cervical cancer cases had nucleotide variations in the HPV-16 E7 region (55). In contrast, we observed a highly conserved E7 region; as high as 86.7% of patients showed no sequence variation, and only 13.3% showed variation in the form of silent mutations, thus having no effect on proteins. Because of the conserved E7, it can be far more easily targeted in order to silence its effects at the RNA or protein level. The most frequently reported E7 variation, A647G, also has been found to be completely absent in India. Any change in the two transforming genes E6 and E7 may lead to altered biological function and oncogenicity of the proteins encoded, which can affect the natural history of HPV infection (13).
Examinations of the L1 gene also did not reveal any significant variations, and 64.5% (20/31) of the samples showed the prototype sequence. Very few reports are available on the HPV-16 L1 variants (9, 53, 54, 56), and none are from south Asia. Interestingly, the common L1 variant T6862C observed in all previous studies was not detected; instead, four novel nucleotide variations, two silent and two missense mutations, at codons 343, 351, 423, and 429 were detected (Table 2). This is indicative of the hypervariability of the L1 gene in this region. The importance of the most frequently observed L1 gene variation, A6695 (T353P), lies in the fact that a polar uncharged amino acid (threonine) is replaced by a nonpolar aliphatic amino acid that may have an effect on the structure or function of the L1 protein, which may play an important role in immune recognition and vaccine development strategies. Amino acid changes among molecular variants also have been shown to affect the efficiency of HPV-16 L1 proteins to self assemble into virus-like particles (25, 49). This variability could lead to conformational changes within epitopes relevant for viral neutralization (44).
The LCR, which contains the upstream regulatory region, is the binding site of various cellular and viral transcription factors that either activate or suppress the p97 promoter activity, which regulates the transcription of various HPV-16 genes, especially the E6 and E7 oncogenes. It has been reported to be the most variable segment of the HPV-16 genome in different populations (23, 26, 53, 56). We observed a similar pattern, as 75% (45/60) of patients had 17 different LCR variants, of which four (C13T, A7636C, C7678T, and G7799A) were found to be novel types, including an A7636C substitution that was detected in the AP1 binding site (4.4%). Since LCR is the major site for the transcriptional control of the virus, studies are under way to correlate the observed alterations with the biological/biochemical properties and different clinicopathological features of the disease. The number of nucleotide alterations was highest for the LCR, followed by the L1, E6, and E7 genes, and this is in good agreement with results from a previous report (56). The most frequently observed mutation in the LCR was G7521A, which is located at the YY1 binding site, was found in as many as 91% of cases, and causes the repression of HPV transcription. This mutation has been found to be uniformly distributed in the majority of cervical cancer patients throughout the world (26, 41, 53, 56) but is not found in asymptomatic carriers (41). Mutations in the YY1 binding sites have been shown to promote the p97 promoter activity by three- to sixfold (12, 29). It also has been shown to quench AP1 activity, thereby repressing HPV-16 transcription (35). The role of mutation at the YY1 binding site in cancer cases assumes a greater importance, as this mutation is absent from noncancerous lesions and asymptomatic carriers (42). The presence of mutations in the binding sites of various transcription factors or in their close vicinities is intriguing, since these factors enhance the transcription activity of HPV-16 oncogenes. A low frequency (4.4%) of mutation found in the Oct-1 binding site appears to be important, since Oct-1 is an important part of the basal transcriptional machinery and is known to down-regulate HPV expression (43).
Most interestingly, when the distribution of HPV-16 variants with respect to E6, E7, and L1 genes and the LCR was compared to that of different lineages, an extremely high prevalence of the E lineage (85, 86.7, 67.7, and 63.3%, respectively) and the complete absence of the As lineage (36, 39) were observed. The occurrence of AA and NA classes in this study is in accordance with other studies from the southern part of India (56). It is important that HPV variants are correlated with the severity of lesions through follow-up data to elucidate the biological significance of HPV-16 variants during cervical carcinogenesis, thus delineating the utility of such correlations in developing reliable diagnostics and effective therapeutics, including vaccines against HPV.

ACKNOWLEDGMENTS
Shailja Pande thanks the Department of Science and Technology,
Government of India, New Delhi, for providing financial support
under the women scientists scheme.
We thank L. Satyanarayana, Suresh Hedau, and Alok C. Bharti of the ICPO for their help with the manuscript.

FOOTNOTES
* Corresponding author. Mailing address: Division of Molecular Oncology, Institute of Cytology and Preventive Oncology, I-7, Sector 39, Noida 201301, India. Phone: 91-120-2578837. Fax: 91-120-2579473. E-mail:
bcdas48{at}hotmail.com 
Published ahead of print on 16 January 2008. 

REFERENCES
1 - Andersson, S., M. Alemi, E. Rylander, A. Strand, B. Larsson, J. Sallstrom, and E. Wilander. 2000. Uneven distribution of HPV 16 E6 prototype and variant (L83V) oncoprotein in cervical neoplastic lesions. Br. J. Cancer 83:307-310.[CrossRef][Medline]
2 - Bernard, H. U., S. Y. Chan, M. M. Manos, C. K. Ong, L. L. Villa, H. Delius, C. L. Peyton, H. M. Bauer, and C. M. Wheeler. 1994. Identification and assessment of known and novel human papillomaviruses by polymerase chain reaction amplification, restriction fragment length polymorphisms, nucleotide sequence, and phylogenetic algorithms. J. Infect. Dis. 170:1077-1085.[Medline]
3 - Bontkes, H. J., M. van Duin, T. D. de Gruijl, M. F. Duggan-Keen, J. M. Walboomers, M. J. Stukart, R. H. Verheijen, T. J. Helmerhorst, C. J. Meijer, R. J. Scheper, F. R. Stevens, P. A. Dyer, P. Sinnott, and P. L. Stern. 1998. HPV 16 infection and progression of cervical intra-epithelial neoplasia: analysis of HLA polymorphism and HPV 16 E6 sequence variants. Int. J. Cancer 78:166-171.[CrossRef][Medline]
4 - Castle, P. E., M. Schiffman, R. Herrero, A. Hildesheim, A. C. Rodriguez, M. C. Bratti, M. E. Sherman, S. Wacholder, R. Tarone, and R. D. Burk. 2005. A prospective study of age trends in cervical human papillomavirus acquisition and persistence in Guanacaste, Costa Rica. J. Infect. Dis. 191:1808-1816.[CrossRef][Medline]
5 - Crook, T., J. A. Tidy, and K. H. Vousden. 1991. Degradation of p53 can be targeted by HPV E6 sequences distinct from those required for p53 binding and trans-activation. Cell 67:547-556.[CrossRef][Medline]
6 - Cuzick J., G. Terry, L. Ho, T. Hollingworth, and M. Anderson. 1992. Human papillomavirus type 16 in cervical smears as predictor of high-grade cervical intraepithelial neoplasia. Lancet 18:959-960.
7 - Das, B. C., J. K. Sharma, V. Gopalakrishna, and U. K. Luthra. 1992. Analysis by polymerase chain reaction of the physical state of human papillomavirus type 16 DNA in cervical preneoplastic and neoplastic lesions. J. Gen. Virol. 73:2327-2336.[Abstract/Free Full Text]
8 - Das, B. C., J. K. Sharma, V. Gopalkrishna, D. K. Das, V. Singh, L. Gissmann, H. zur Hausen, and U. K. Luthra. 1992. A high frequency of human papillomavirus DNA sequences in cervical carcinomas of Indian women as revealed by Southern blot hybridization and polymerase chain reaction. J. Med. Virol. 36:239-245.[CrossRef][Medline]
9 - de Boer, M. A., L. A. Peters, M. F. Aziz, B. Siregar, S. Cornain, M. A. Vrede, E. S. Jordanova, S. Kolkman-Uljee, and G. J. Fleuren. 2004. Human papillomavirus type 16 E6, E7, and L1 variants in cervical cancer in Indonesia, Suriname, and The Netherlands. Gynecol. Oncol. 94:488-494.[CrossRef][Medline]
10 - de Villiers, E. M. 1994. Human pathogenic papillomavirus types: an update. Curr. Top. Microbiol. Immunol. 186:1-12.[Medline]
11 - de Villiers, E. M., C. Fauquet, T. R. Broker, H. U. Bernard, and H. zur Hausen. 2004. Classification of papillomaviruses. Virology 324:17-27.[CrossRef][Medline]
12 - Dong, X. P., and H. Pfister. 1999. Overlapping YY1- and aberrant SP1-binding sites proximal to the early promoter of human papillomavirus type 16. J. Gen. Virol. 80:2097-2101.[Abstract/Free Full Text]
13 - Duensing, S., and K. Munger. 2002. The human papillomavirus type 16 E6 and E7 oncoproteins independently induce numerical and structural chromosome instability. Cancer Res. 62:7075-7082.[Abstract/Free Full Text]
14 - Dyson, N., P. M. Howley, K. Münger, and E. Harlow. 1989. The human papilloma virus-16 E7 oncoprotein is able to bind to the retinoblastoma gene product. Science 243:934-937.[Abstract/Free Full Text]
15 - Ellis, J. R., P. J. Keating, J. Baird, E. F. Hounsell, D. V. Renouf, M. Rowe, D. Hopkins, M. F. Duggan-Keen, J. S. Bartholomew, L. S. Young, et al. 1995. The association of an HPV16 oncogene variant with HLA-B7 has implications for vaccine design in cervical cancer. Nat. Med. 1:464-470.[CrossRef][Medline]
16 - Eschle, D., M. Durst, J. ter Meulen, J. Luande, H. C. Eberhardt, M. Pawlita, and L. Gissmann. 1992. Geographical dependence of sequence variation in the E7 gene of human papillomavirus type 16. J. Gen. Virol. 73:1829-1832.[Abstract/Free Full Text]
17 - Gupta, A., R. Arora, S. Gupta, B. K. Prusty, U. Kailash, S. Batra, and B. C. Das. 2006. Human papillomavirus DNA in urine samples of women with or without cervical cancer and their male partners compared with simultaneously collected cervical/penile smear or biopsy specimens. J. Clin. Virol. 37:190-194.[CrossRef][Medline]
18 - Hildesheim, A., M. Schiffman, C. Bromley, S. Wacholder, R. Herrero, A. Rodriguez, M. C. Bratti, M. E. Sherman, U. Scarpidis, Q. Q. Lin, M. Terai, R. L. Bromley, K. Buetow, R. J. Apple, and R. D. Burk. 2001. Human papillomavirus type 16 variants and risk of cervical cancer. J. Natl. Cancer Inst. 93:315-318.[Free Full Text]
19 - Hu, X., T. Pang, Z. Guo, N. Mazurenko, F. Kisseljov, J. Ponten, and M. Nister. 2001. HPV16 E6 gene variations in invasive cervical squamous cell carcinoma and cancer in situ from Russian patients. Br. J. Cancer 84:791-795.[CrossRef][Medline]
20 - Jain, N., V. Singh, S. Hedau, S. Kumar, M. K. Daga, R. Dewan, N. S. Murthy, S. A. Husain, and B. C. Das. 2005. Infection of human papillomavirus type 18 and p53 codon 72 polymorphism in lung cancer patients from India. Chest 128:3999-4007.[CrossRef][Medline]
21 - Jiang, M., and J. Milner. 2002. Selective silencing of viral gene expression in HPV-positive human cervical carcinoma cells treated with siRNA, a primer of RNA interference. Oncogene 21:6041-6048.[CrossRef][Medline]
22 - Kailash, U., C. C. Soundararajan, R. Lakshmy, R. Arora, S. Vivekanandhan, and B. C. Das. 2006. Telomerase activity as an adjunct to high-risk human papillomavirus types 16 and 18 and cytology screening in cervical cancer. Br. J. Cancer 95:1250-1257.[CrossRef][Medline]
23 - Kämmer, C., U. Warthorst, N. Torrez-Martinez, C. M. Wheeler, and H. Pfister. 2000. Sequence analysis of the long control region of human papillomavirus type 16 variants and functional consequences for P97 promoter activity. J. Gen. Virol. 81:1975-1981.[Abstract/Free Full Text]
24 - Kast, W. M., R. M. Brandt, J. Sidney, J. W. Drijfhout, R. T. Kubo, H. M. Grey, C. J. Melief, and A. Sette. 1994. Role of HLA-A motifs in identification of potential CTL epitopes in human papillomavirus type 16 E6 and E7 proteins. J. Immunol. 152:3904-3912.[Abstract]
25 - Kirnbauer, R., J. Taub, H. Greenstone, R. Roden, M. Durst, L. Gissmann, D. R. Lowy, and J. T. Schiller. 1993. Efficient self-assembly of human papillomavirus type 16 L1 and L1-L2 into virus-like particles. J. Virol. 67:6929-6936.[Abstract/Free Full Text]
26 - Kurvinen, K., M. Yliskoski, S. Saarikoski, K. Syrjanen, and S. Syrjanen. 2000. Variants of the long control region of human papillomavirus type 16. Eur. J. Cancer 36:1402-1410.[CrossRef][Medline]
27 - Londesborough, P., L. Ho, G. Terry, J. Cuzick, C. Wheeler, and A. Singer. 1996. Human papillomavirus genotype as a predictor of persistence and development of high-grade lesions in women with minor cervical abnormalities. Int. J. Cancer 69:364-368.[CrossRef][Medline]
28 - Matsumoto, K., H. Yoshikawa, S. Nakagawa, X. Tang, T. Yasugi, K. Kawana, S. Sekiya, Y. Hirai, I. Kukimoto, T. Kanda, and Y. Taketani. 2000. Enhanced oncogenicity of human papillomavirus type 16 (HPV16) variants in Japanese population. Cancer Lett. 156:159-165.[CrossRef][Medline]
29 - May, M., X. P. Dong, E. Beyer-Finkler, F. Stubenrauch, P. G. Fuchs, and H. Pfister. 1994. The E6/E7 promoter of extrachromosomal HPV16 DNA in cervical cancers escapes from cellular repression by mutation of target sequences for YY1. EMBO J. 13:1460-1466.[Medline]
30 - Münger, K., and P. M. Howley. 2002. Human papillomavirus immortalization and transformation functions. Virus Res. 89:213-228.[CrossRef][Medline]
31 - Muñoz, N., F. X. Bosch, S. de Sanjose, R. Herrero, X. Castellsague, K. V. Shah, P. J. Snijders, and C. J. Meijer. 2003. Epidemiologic classification of human papillomavirus types associated with cervical cancer. N. Engl. J. Med. 348:518-527.[Abstract/Free Full Text]
32 - Myers, G., H. U. Bernard, H. Delius, C. Baker, J. Icenogle, A. Halpern, and C. Wheeler. 1995. Human papillomavirus. A compilation of analysis of nucleic acid and amino acid sequences. Publication LA-UR 95-3675. Los Alamos National Laboratory, Los Alamos, NM.
33 - Nindl, I., K. Rindfleisch, B. Lotz, A. Schneider, and M. Durst. 1999. Uniform distribution of HPV 16 E6 and E7 variants in patients with normal histology, cervical intra-epithelial neoplasia and cervical cancer. Int. J. Cancer 82:203-207.[CrossRef][Medline]
34 - Nobbenhuis, M. A., J. M. Walboomers, T. J. Helmerhorst, L. Rozendaal, A. J. Remmink, E. K. Risse, H. C. van der Linden, F. J. Voorhorst, P. Kenemans, and C. J. Meijer. 1999. Relation of human papillomavirus status to cervical lesions and consequences for cervical-cancer screening: a prospective study. Lancet 354:20-25.[CrossRef][Medline]
35 - O'Connor, M., and H. U. Bernard. 1995. Oct-1 activates the epithelial-specific enhancer of human papillomavirus type 16 via a synergistic interaction with NFI at a conserved composite regulatory element. Virology 207:77-88.[CrossRef][Medline]
36 - Radhakrishna Pillai, M., S. Sreevidya, B. H. Pollock, P. G. Jayaprakash, and B. Herman. 2002. Human papillomavirus type 16 E6 and E7 gene variations in Indian cervical cancer. Gynecol. Oncol. 87:268-273.[CrossRef][Medline]
37 - Rastogi, T., A. Hildesheim, and R. Sinha. 2004. Opportunities for cancer epidemiology in developing countries. Nat. Rev. Cancer 4:909-917.[CrossRef][Medline]
38 - Resnick, R. M., M. T. Cornelissen, D. K. Wright, G. H. Eichinger, H. S. Fox, J. ter Schegget, and M. M. Manos. 1990. Detection and typing of human papillomavirus in archival cervical cancer specimens by DNA amplification with consensus primers. J. Natl. Cancer Inst. 82:1477-1484.[Abstract/Free Full Text]
39 - Sathish, N., P. Abraham, A. Peedicayil, G. Sridharan, and G. Chandy. 2005. HPV 16 E6 sequence variations in Indian patients with cervical neoplasia. Cancer Lett. 229:93-99.[CrossRef][Medline]
40 - Schlecht, N. F., R. W. Platt, E. Duarte-Franco, M. C. Costa, J. P. Sobrinho, J. C. Prado, A. Ferenczy, T. E. Rohan, L. L. Villa, and E. L. Franco. 2003. Human papillomavirus infection and time to progression and regression of cervical intraepithelial neoplasia. J. Natl. Cancer Inst. 95:1336-1343.[Abstract/Free Full Text]
41 - Schmidt, M., W. Kedzia, and A. Gozdzicka-Jozefiak. 2001. Intratype HPV16 sequence variation within LCR of isolates from asymptomatic carriers and cervical cancers. J. Clin. Virol. 23:65-77.[CrossRef][Medline]
42 - Seedorf, K., G. Krammer, M. Durst, S. Suhai, and W. G. Rowekamp. 1985. Human papillomavirus type 16 DNA sequence. Virology 145:181-185.[CrossRef][Medline]
43 - Sibbet, G. J., S. Cuthill, and M. S. Campo. 1995. The enhancer in the long control region of human papillomavirus type 16 is up-regulated by PEF-1 and down-regulated by Oct-1. J. Virol. 69:4006-4011.[Abstract]
44 - Sichero, L., S. Ferreira, H. Trottier, E. Duarte-Franco, A. Ferenczy, E. L. Franco, and L. L. Villa. 2007. High grade cervical lesions are caused preferentially by non-European variants of HPVs 16 and 18. Int. J. Cancer 120:1763-1768.[CrossRef][Medline]
45 - Song, Y. S., S. H. Kee, J. W. Kim, N. H. Park, S. B. Kang, W. H. Chang, and H. P. Lee. 1997. Major sequence variants in E7 gene of human papillomavirus type 16 from cervical cancerous and noncancerous lesions of Korean women. Gynecol. Oncol. 66:275-281.[CrossRef][Medline]
46 - Stewart, A. C., A. M. Eriksson, M. M. Manos, N. Munoz, F. X. Bosch, J. Peto, and C. M. Wheeler. 1996. Intratype variation in 12 human papillomavirus types: a worldwide perspective. J. Virol. 70:3127-3136.[Abstract]
47 - Stöppler, M. C., K. Ching, H. Stoppler, K. Clancy, R. Schlegel, and J. Icenogle. 1996. Natural variants of the human papillomavirus type 16 E6 protein differ in their abilities to alter keratinocyte differentiation and to induce p53 degradation. J. Virol. 70:6987-6993.[Abstract/Free Full Text]
48 - Tornesello, M. L., M. L. Duraturo, I. Salatiello, L. Buonaguro, S. Losito, G. Botti, G. Stellato, S. Greggi, R. Piccoli, S. Pilotti, B. Stefanon, G. De Palo, S. Franceschi, and F. M. Buonaguro. 2004. Analysis of human papillomavirus type-16 variants in Italian women with cervical intraepithelial neoplasia and cervical cancer. J. Med. Virol. 74:117-126.[CrossRef][Medline]
49 - Touze, A., S. El Mehdaoui, P. Y. Sizaret, C. Mougin, N. Munoz, and P. Coursaget. 1998. The L1 major capsid protein of human papillomavirus type 16 variants affects yield of virus-like particles produced in an insect cell expression system. J. Clin. Microbiol. 36:2046-2051.[Abstract/Free Full Text]
50 - van Duin, M., P. J. Snijders, M. T. Vossen, E. Klaassen, F. Voorhorst, R. H. Verheijen, T. J. Helmerhorst, C. J. Meijer, and J. M. Walboomers. 2000. Analysis of human papillomavirus type 16 E6 variants in relation to p53 codon 72 polymorphism genotypes in cervical carcinogenesis. J. Gen. Virol. 81:317-325.[Abstract/Free Full Text]
51 - Veress, G., M. Murvai, K. Szarka, A. Juhasz, J. Konya, and L. Gergely. 2001. Transcriptional activity of human papillomavirus type 16 variants having deletions in the long control region. Eur. J. Cancer 37:1946-1952.[CrossRef][Medline]
52 - Villa, L. L., P. Rahal, and E. L. Franco. 1997. Molecular variant analysis as a tool in natural history studies of human papillomavirus infection and cervical neoplasia, p. 379-385. In E. L. Franco and J. Monsonego (ed.), New developments in cervical cancer screening and prevention. Blackwell Science, Oxford, United Kingdom.
53 - Villa, L. L., L. Sichero, P. Rahal, O. Caballero, A. Ferenczy, T. Rohan, and E. L. Franco. 2000. Molecular variants of human papillomavirus types 16 and 18 preferentially associated with cervical neoplasia. J. Gen. Virol. 81:2959-2968.[Abstract/Free Full Text]
54 - Wheeler, C. M., T. Yamada, A. Hildesheim, and S. A. Jenison. 1997. Human papillomavirus type 16 sequence variants: identification by E6 and L1 lineage-specific hybridization. J. Clin. Microbiol. 35:11-19.[Abstract]
55 - Wu, Y., Y. Chen, L. Li, G. Yu, Y. He, and Y. Zhang. 2006. Analysis of mutations in the E6/E7 oncogenes and L1 gene of human papillomavirus 16 cervical cancer isolates from China. J. Gen. Virol. 87:1181-1188.[Abstract/Free Full Text]
56 - Yamada, T., M. M. Manos, J. Peto, C. E. Greer, N. Munoz, F. X. Bosch, and C. M. Wheeler. 1997. Human papillomavirus type 16 sequence variation in cervical cancers: a worldwide perspective. J. Virol. 71:2463-2472.[Abstract]
57 - Yamada, T., C. M. Wheeler, A. L. Halpern, A. C. Stewart, A. Hildesheim, and S. A. Jenison. 1995. Human papillomavirus type 16 variant lineages in United States populations characterized by nucleotide sequence analysis of the E6, L2, and L1 coding segments. J. Virol. 69:7743-7753.[Abstract]
58 - Zehbe, I., G. Voglino, H. Delius, E. Wilander, and M. Tommasino. 1998. Risk of cervical cancer and geographical variations of human papillomavirus 16 E6 polymorphisms. Lancet 352:1441-1442.[CrossRef][Medline]
59 - Zehbe, I., E. Wilander, H. Delius, and M. Tommasino. 1998. Human papillomavirus 16 E6 variants are more prevalent in invasive cervical carcinoma than the prototype. Cancer Res. 58:829-833.[Abstract/Free Full Text]
60 - zur Hausen, H. 2002. Papillomaviruses and cancer: from basic studies to clinical application. Nat. Rev. Cancer 2:342-350.[CrossRef][Medline]
Journal of Clinical Microbiology, March 2008, p. 1060-1066, Vol. 46, No. 3
0095-1137/08/$08.00+0 doi:10.1128/JCM.02202-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.