1 Albert Einstein College of Medicine, Department of Epidemiology and Population Health, 1300 Morris Park Avenue, Bronx, NY 10461, USA
2 University of California, San Francisco, San Francisco, CA, USA
3 Maimonides Medical Center, Brooklyn, NY, USA
4 Southern Illinois University School of Medicine, Springfield, IL, USA
5 The Henry M. Jackson Foundation, Rockville, MD, USA
6 University of Southern California, Los Angeles, CA, USA
7 Johns Hopkins Bloomberg School of Public Health, Baltimore, MD, USA
8 National Institute of Child Health and Human Development, Bethesda, MD, USA
Correspondence
Nicolas F. Schlecht
nschlech{at}aecom.yu.edu
![]() |
ABSTRACT |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
INTRODUCTION |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Several epidemiological studies have found that certain HPV16 and 18 variants may be more strongly associated with cervical cancer (Burk et al., 2003; Hildesheim et al., 2001
; Ho et al., 1993
; Londesborough et al., 1996
; Nindl et al., 1999
; Villa et al., 2000
; Xi et al., 1995
, 1997
; Zehbe et al., 1998
, 2001b
). By convention, variants are defined as HPVs that differ by <5 % of the nucleotide sequence in non-coding regions and <2 % of nucleotides in coding regions (Bernard et al., 1994
). HPVs with genotype differences greater than this are considered to be different HPV types. HPV16 and 18 variants are typically classified into European, African, Asian-American, North-American or Amerindian lineages according to their population of origin (Villa et al., 2000
; Yamada et al., 1995
).
The most consistent observation has been a higher prevalence of non-European variants of HPV16 in cervical cancer and high-grade cervical lesions than would be expected based on their relative prevalence in normal tissues. Taken as a whole, recent studies suggest that there is a two- to threefold greater association of non-European variants with high-grade cervical lesions than is found for European variants of HPV16 (Hildesheim et al., 2001; Xi et al., 2002
).
There have been far fewer studies of HPV18, but the limited data available suggest that, as for HPV16, non-European variants of HPV18 may be more common than expected in cancer specimens (Burk et al., 2003; Hecht et al., 1995
; Lizano et al., 1997
; Terry et al., 1997
) and high-grade cervical lesions (Villa et al., 2000
). Confirmation of these results in large, well-designed studies is greatly needed, since the number of HPV18-positive specimens genotyped in earlier investigations was small, and the algorithms for classifying the non-European variant groups varied (Burk et al., 2003
; Villa et al., 2000
).
Since almost all studies to date have been cross-sectional or case-control investigations, very little is known regarding the natural history of HPV16 and 18 variants to help explain the possible differences in their associations with cervical cancer. These differences could involve events early in the course of the multistage process of cervical tumorigenesis, such as differences in viral persistence and the risk of developing cervical lesions. Alternatively, differences in the risk of cancer between HPV16 and 18 variants could relate to later events, such as risk of transition from high-grade lesions to cancer. Additionally, highly oncogenic HPV16 and 18 variants might be expected to have a particularly aggressive natural history in the presence of diminished host immune status, such as in HIV-positive women (Chaturvedi et al., 2004; Icenogle et al., 1992
; Mayrand et al., 2000
; Perez-Gallego et al., 2001
; Xi et al., 1998
). Thus, HIV-positive women represent a unique opportunity to study the effects of host immune status on variant type and behaviour. To date, however, there have been only a few small studies in HIV-positive women, each involving fewer than 40 HPV16- or 18-positive women. These small studies have reported conflicting results regarding the effects of HIV coinfection on the distribution of HPV variants (Chaturvedi et al., 2004
; Icenogle et al., 1992
; Perez-Gallego et al., 2001
), and none have prospectively assessed the effects of variant lineage on the natural history of infection.
Therefore, in a large prospective cohort of HIV-positive and -negative women, the Women's Interagency HIV Study (WIHS), we determined the distribution of HPV16 and 18 variant lineages, the associations of variant lineage with race, HIV coinfection and immune status, and whether variant lineages differed in their persistence and/or their associations with squamous intraepithelial lesions (SIL).
![]() |
METHODS |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Clinical laboratory data.
Pap smears were examined at a single laboratory and categorized using the 1991 Bethesda system for cytological diagnosis. Two independent cytotechnologists examined all Pap smears. Any Pap smear considered abnormal by either cytotechnologist, as well as 10 % of all negative smears, were examined by a cytopathologist. Women were referred for colposcopy for any abnormal cytology result. Treatment decisions were individualized, but the study-wide protocol recommended excision or ablation for women with known or suspected high-grade cervical disease. T-cell subsets were determined by immunofluorescence using flow cytometry in laboratories participating in the AIDS Clinical Trials Quality Assurance Program. Quantification of HIV-1 RNA in plasma was performed using the isothermal nucleic acid sequence-based amplification (NASBA/Nuclisens) method (Organon Teknika) in laboratories participating in the National Institutes of Health (NIH)/ National Institute of Allergy and Infectious Disease (NIAID) Virology Quality Assurance Laboratory proficiency testing program. The lower limit of quantification through September 1997 was 4000 copies ml1 using a 0·1 ml sample; from October 1997 to December 1998 the lower limit was 400 copies ml1 using 0·2 ml; after January 1999 the lower limit was 80 copies ml1 using 1·0 ml.
Detection of HPV DNA.
Exfoliated cervicovaginal cells were tested for the presence of HPV DNA, following digestion with proteinase K, using a well-established PCR protocol that amplifies a highly conserved 450 bp segment in the L1 viral gene (flanked by primers MY09/MY11/HMB01) (Palefsky et al., 1999). In brief, following proteinase K digestion, 210 µl of each cell digest was used in reactions containing 10 mM Tris/HCl, 50 mM KCl, 4 mM MgCl2, 200 µM of each deoxyribonucleotide triphosphate, 2·5 U Taq DNA polymerase and 0·5 µM of each primer. There were 35 amplification cycles (95 °C for 20 s, 55 °C for 30 s and 72 °C for 30 s, with a 5 min extension period at 72 °C on the last cycle). Primer set PC04/GH20, which amplifies a 268 bp cellular
-globin DNA fragment, was included in each assay as an internal control. The amplification products were probed for the presence of HPV DNA using a generic probe mixture, as well as for HPV DNA on a type-specific basis using filters individually hybridized with biotinylated oligonucleotide probes for HPV types 6, 11, 13, 16, 18, 26, 31, 32, 33, 34, 35, 39, 40, 42, 45, 51, 52, 53, 54, 55, 56, 57, 58, 59, 61, 62, 64, 66, 67, 68, 69, 70, 71 (AE8), 72, 73 (PAP238A), 81 (AE7), 83 (PAP291), 82 (W13B/AE2), 84 (PAP155), 85 (AE5), 89 (AE6), AE9 and AE10.
HPV16/18 variant analysis.
All specimens found to be positive for either HPV16 or 18 from any of the first four visits were assessed for the variant present. Variant classification was determined by sequencing of the E6 regions of HPV16 and 18. The E6 region for HPV16 was amplified using a previously described protocol and primer set (Wheeler et al., 1997). The primers used in the direct sequencing reactions were specific for HPV16 or 18 and did not amplify other HPV types, thus reducing ambiguity due to mixed infections with other HPV types. The HPV18 E6 region was amplified by single-tube nested PCR using the outer primers AGTAACCGAAAACGGTCGGG (forward) and CGGGCTGGTAAATGTTGATGA (reverse), and inner primers CGGTGTATATAAAAGATG (forward) and TGCTCGTGACATAGAA (reverse). PCR products were then isolated by gel electrophoresis using the Qiagen Gel extraction kit and sequenced on both strands for the entire E6 regions of the HPV16 and 18 genomes. The upstream regulatory region (URR) was sequenced in seven HPV16-positive specimens that could not be sequenced for E6. The URR was similarly amplified by nested PCR and isolated by gel electrophoresis using a separate protocol (Burk et al., 2003
). In total, 243 HPV16- and 225 HPV18-positive samples identified by L1 PCR were tested with a success rate for amplification of the E6 region of approximately 75 %, which did not vary by variant group, HPV viral load (as determined by MY09/11 PCR dot-blot intensity) or CD4+ T-cell count.
For efficiency, we determined the variant present in only the first HPV16- or 18-positive CVL in each series of sequential positive specimens (e.g. the CVL at visit 3, for a patient who was HPV16-positive at visits 3, 4 and 5). To confirm that this approach did not result in misclassification of the variant present at later visits (e.g. at visits 4 and 5), we determined the variant present in the very last HPV16/18-positive specimen of each series in a random sample of 20 % of subjects (n=82). Only in one instance did we detect a change in variant lineage classification (HPV18 African to European). In an additional subsample of 22 women with HPV16 infection that persisted for three or more visits, we also sequenced the HPV16 present at each and every visit, and found all samples to contain the same variant and sequence at each visit. Based on the consistency of these results, we deemed it appropriate to assume that the same variant was present in sequential visits when HPV16 or 18 infections persisted.
E6 variant classification.
Variants were grouped according to lineages categorized as previously reported (Ho et al., 1993; Ong et al., 1993
; Yamada et al., 1995
). HPV16 variant lineages were designated European (E), African-1 (Af-1), African-2 (Af-2), North-American (NA) or Asian-American (AA) (Fig. 1
a). HPV18 variant lineages were likewise classified into E, Af-1, Af-2 or Asian/Amerindian (As/Ai) lineages. By convention, for HPV16, the E variant containing a thymine (T) at nucleotide position 350 was designated the primary reference sequence to which all other HPV16 variants were compared (Yamada et al., 1995
), whereas convention designates the As/Ai variant as the HPV18 primary reference sequence (Ong et al., 1993
). Our principal comparison, set a priori, was between E variants of HPV16 or 18 and non-E lineages, since, as discussed above, these were the groupings most consistently found to have differences in their associations with cervical cancer.
|
To study the effects of HPV16 and 18 variant lineage on the persistence of infection we used multivariable Cox proportional hazards regression. Specifically, persistence was measured using time-to-clearance of an HPV16 or 18 infection that was newly detected (studying incidently detected HPV makes it possible to estimate time of onset), with clearance defined using a stringent criterion of at least two sequential negative HPV-PCR results (to minimize concerns that a single false-negative result might affect the findings). The mid-point calendar date between two consecutive visits was used to estimate the time of event. Because this analysis was restricted to incident infections the data were more limited. Therefore, we used non-oncogenic HPV types as the reference group, with HPV16 and 18 variant lineage included as a variable in each multivariable Cox-model; that is, we calculated the hazard ratio (HR) for time-to-clearance for each variant lineage relative to that of all non-oncogenic HPV types, which includes a broad category of HPV types with many incident and clearance events so as to provide a stable statistical base for estimation of HRs. This approach also made it possible to stratify the baseline hazard by HPV-type viral load (determined in MY09/11 PCR), HIV serostatus and CD4+ T-cell count detected at time of infection, to best control for these important factors in the analysis. Further adjustment for age, ethnic group and other potential confounders (see above) was achieved by including these variables in the multivariable Cox models. The Wald test was used to measure the statistical significance of differences in the HRs between variant groups. Proportional hazards assumptions were checked and no violations were found.
Lastly, to study the cross-sectional relationship of HPV16 and 18 variant lineage with SIL, we used multivariable logistic regression controlling for age, HPV-viral load, CD4+ T-cell count, HIV serostatus, ethnic group and whether the HPV infection was prevalent or incident. Restricting observations to incidently detected infections alone did not change the associations. Statistical analyses were performed using Stata-8 and SAS-9.
![]() |
RESULTS |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Fig. 1(a) presents the distribution of HPV16 variants observed in our study and their lineage classification. E variants were the most common HPV16 group detected (64 %, n=124), followed by Af (28 %, n=55) and AA/NA (8 %, n=16) variants. Of the E variants detected, most (n=55) were the primary reference sequence (28 %); 31 (16 %) were similar to the primary reference sequence except for a single nucleotide difference at position 350, a G (guanine) instead of a T (thymine), designated E-350G, which can result in a leucine-to-valine amino-acid change; and 19 (10 %) differed from these two major E subgroups by a single nucleotide. An additional 13 (7 %) were considered Asian variants because of a G-to-T substitution at position 178, though still part of the E lineage. Among the common non-E variants detected, 37 (19 %) were categorized as the Af-1 subgroup reference, 8 (4 %) as the Af-2 subgroup reference and 12 (6 %) as the AA subgroup reference; the remainder had small differences (mostly single nucleotide changes) from these three major sequences. Notably, 11 (6 %) HIV-positive subjects appeared to have coinfection with more than one HPV16 variant, as determined by the detection of two different nucleotides at a single position. Instances of HPV16 coinfection were also more likely to occur in subjects with poor immune status (i.e. seven subjects had a CD4+ T-cell count of <200 cells mm3 and an HIV-RNA level of >20 000 copies ml1). Seven additional HPV16-positive specimens that could not be amplified by E6 PCR were instead amplified and sequenced within the URR (not shown) and found to include six E and one AA isolate. We also identified two novel HPV16 variant sequences: an NA variant with C (cytosine) instead of T at position 109, and an As variant with G instead of T at position 183. Interestingly, the latter substitution also likely results in a non-synonomous (non-silent) amino-acid change: isoleucine to arginine.
HPV18 was less heterogeneous than HPV16, with a small group of variants accounting for most HPV18-positive specimens. In addition, unlike HPV16, the most common HPV18 variants were Af (48 %, n=77), followed by E (38 %, n=62) and As/Ai (14 %, n=23) (Fig. 1b), and the vast majority of the variants contained the exact subgroup reference sequence for each lineage. The exception was Af-2, which had a common variant (n=14) containing a T instead of a G at position 342. Coinfection with more than one variant was rare in HPV18-positive women, being detected in just 2/162 (1·2 %) of HPV-positive women. Three E samples were also positive for the T18-7 variant (see Fig. 1b
) (de Boer et al., 2005
), but with additional (silent) substitutions at positions 149 and 377.
Interestingly, we found that all of the major HPV16 and 18 variant groups detected in our population, and even their subgroups, could be distinguished based on the nucleotides found at just a small number of key positions. With respect to HPV16, these positions were 131G, 132C/T, 143/5GT, 178G, 350G and 532G (as shown in Fig. 1a). For example, Asian variants, which are a subgroup of the E lineage, could be distinguished from the E primary reference sequence solely by the presence of a G at position 178. Non-E variants, including AA, NA and Af lineages, could be distinguished from the E primary reference sequence by a T at position 145 (or alternatively at 335). Within the Af lineage, the two main subgroups, Af-1 and Af-2, could be distinguished from the E primary reference sequence by a C-to-G substitution at position 143 and/or either G-to-C (Af-1) or G-to-T (Af-2) substitutions at position 132. An analogous algorithm could be derived to distinguish the four major HPV18 variant groups detected in our population, based on departures from the As/Ai primary reference sequence at three positions in the E6 region: 317C, 548G and 549A. All HPV18 E variants carried a C-to-A (adenine) substitution at position 549, whereas Af variants carried an additional A-to-G substitution at position 548. Af variants could be further divided into Af-1 and Af-2 subgroups, based on a T-to-C substitution at position 317.
Table 1 shows the distribution of HPV16 and 18 variants by race. Caucasians and Hispanics/Others had similar distributions of HPV16 variants; the majority of Caucasians (81 %) and those classified as Hispanic/Other (81 %) were positive for HPV16 E variants, and had similar distributions of E variant subgroups (not shown) as well as a similar prevalence of Af and AA/NA variants. African-Americans (including Hispanics of African-American origin), in contrast, had a higher prevalence of HPV16 Af variants (42 %) than either Caucasians (14 %) (P=0·001) or those classified as Hispanic/Other (9 %) (P<0·001). Consistent with this, African-Americans also had a significantly higher prevalence of HPV18 Af variants (60 %) than Caucasians (13 %; P<0·001). However, the distribution of HPV18 variants in women classified as Hispanic/Other differed from Caucasians. The prevalence of HPV18 Af (38 %), As/Ai (15 %) and E (48 %) in the Hispanic/Other group was intermediate between African-Americans and Caucasians. Adjustment for potential confounders, including HIV serostatus, CD4+ T-cell count, number of sexual partners, HPV-viral load and study recruitment site did not alter the strong associations between race and the distribution of HPV16 and 18 variants. The distribution of variants between the different racial groups varied little across the different recruitment sites (Brooklyn, Bronx, Washington DC, Chicago, Los Angeles and San Francisco) or across the other demographic characteristics listed previously. Similar findings (not shown) were observed in analyses that excluded Other ethnicities from the Hispanic/Other group (i.e. a group composed of only Hispanics) and compared HPV16 and 18 variant distributions by race.
|
|
|
|
![]() |
DISCUSSION |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
A number of different HPV16 and 18 variants were detected in this population of women at high risk of HPV infection. Most of these variants were categorized as having one of a limited number of common sequences that have been previously identified (Ong et al., 1993; Yamada et al., 1995
). There was a low prevalence of HPV16 AA variants compared to E and Af variants (Yamada et al., 1997
), but this was not unusual given the small number of Hispanic subjects found in our population (Berumen et al., 2001
; Casas et al., 1999
; Chaturvedi et al., 2004
; del Refugio Gonzalez-Losa et al., 2004
; Picconi et al., 2003
; Villa et al., 2000
). More surprising was the detection of a T-to-G substitution at position 183 in 12 HPV16 E variants which has not been reported among E variants except in a small number of cases of in situ or invasive cervical cancer (Burk et al., 2003
; Matsumoto et al., 2000
; Zehbe et al., 1998
). Additionally, we detected two novel HPV16 variant sequences: an NA variant with a T-to-C substitution at 109, and an As variant with a T-to-G substitution at 183, the latter of which also likely results in a non-synonymous amino-acid change. Furthermore, three HPV18 E samples were positive for the T18-7 variant that has been identified previously (de Boer et al., 2005
; Ong et al., 1993
) but with additional (silent) substitutions at positions 149 and 377.
Of special note, we found that a small number of key nucleotide positions could be used to distinguish all of the major variant lineages and perhaps even their subgroups. With respect to HPV16, these positions were 131/2, 143/5, 178, 350 and 532, and for HPV18, the relevant nucleotide positions were 317, 548 and 549. In addition to these, a polymorphic site within the HPV16 E6 gene (at position 276) has also been identified in an E subgroup variant found in Asian and European populations (de Boer et al., 2004; Hu et al., 2001a
, b
; Kammer et al., 2002
; Matsumoto et al., 2000
; Zehbe et al., 2001a
). If confirmed by other studies, it may be possible to exploit these key nucleotide positions to develop simplified high-throughput methods for identifying the variants present in cervical specimens (Gemignani et al., 2004
; Wheeler et al., 1997
).
The greatest predictor of variant lineage was race. In particular, Af variants of HPV16 and 18 were significantly more common in African-Americans than in Caucasians. While this may initially seem an obvious finding, it is of scientific interest to note that even outside Africa, Af variants of HPV16 and 18 commonly infect women of African descent. The reason for this is not concretely known. Although it is likely to reflect race preference in sexual mixing behaviours (i.e. women of African descent may be more likely to have sexual partners of African descent), immunogenetic (e.g. human leukocyte antigen haplotype) and other factors may also play a role (Giannoudis & Herrington, 2001; Terry et al., 1997
; van Duin et al., 2000
; Zehbe et al., 2001a
).
We observed no association between the HPV16 and 18 variant lineage and HIV coinfection or markers of immune status (i.e. CD4+ T-cell count and HIV-RNA level) among HIV-positive women. Prior data from the literature are limited. Among the few studies that included HIV-positive women, two small studies (n<25 in both) reported non-significant increases in detection of non-E variants for HPV16 (Icenogle et al., 1992; Perez-Gallego et al., 2001
), and one showed a significantly greater prevalence of 350G variants (Chaturvedi et al., 2004
). While the authors of the latter study controlled for age, no adjustments were made for immune status (i.e. CD4+ T-cell count and HIV-RNA level) or race, and 35 % (n=33) of the women had SIL, of which two-thirds were high-grade. An earlier study of HPV16 variants in anal specimens from HIV-positive and -negative men and women in San Francisco, which involved a subset of women from the WIHS and men seen in the same clinical centre, found no association between HPV16 non-E variants and HIV serostatus (Da Costa et al., 2002
). Interestingly, we observed a slightly higher frequency of coinfection with more than one HPV16 variant than in some previous studies (Wheeler et al., 1997
), which in our population was largely restricted to HIV-positive women with poor immune status (i.e. CD4+ T-cell count <200 cells mm3 and HIV-RNA level >20 000 copies ml1). If mixed HPV-variant infections are in fact higher among immune-compromised HIV-infected individuals, this observation may potentially support a model for immunological protection. However, such infections represented only a fraction of HPV-positive subjects in our study (less than 6 %), and further investigation into their potential association with immune status and effect on HPV persistence is needed. Taken as a whole, there are currently little data to suggest that the distribution of cervical HPV16 or 18 variants is affected by HIV coinfection or host immune status.
There are also scant prior data regarding the association of HPV16 and 18 variant lineage with the natural history of infection in either HIV-positive or -negative women. We detected no significant differences in time-to-clearance among HPV16 or 18 variant lineages. While we also unexpectedly found non-oncogenic HPV types tended to persist longer than either HPV16 or 18 types, the relative time-to-clearance of oncogenic and non-oncogenic HPV has varied between studies (Franco et al., 1999; Giuliano et al., 2002
; Liaw et al., 2001
; Richardson et al., 2003
). Our finding is not unique (Ahdieh et al., 2001
; Moscicki et al., 2004
) and may relate to subset analysis (HPV16 and 18 alone) (Rousseau et al., 2001
). In HIV-positive adolescents, Moscicki et al. (2004)
found no differences in persistence of HPV16/18 and non-oncogenic types. The few published studies that have compared the persistence of HPV16 between variant lineages report conflicting findings (Bontkes et al., 1998
; Londesborough et al., 1996
; Villa et al., 2000
; Xi et al., 2002
). A common limitation to these typically small studies was that they combined both prevalent and incident infections (Bontkes et al., 1998
; Londesborough et al., 1996
; Villa et al., 2000
), even though prevalent HPV infections are more likely to persist than are incidently detected HPV infections, and the time of onset cannot be estimated with prevalent infections (Ho et al., 1998
). Consistent with our results, an earlier study by Xi et al. (2002)
, which also measured persistence of incident HPV16 infections (n=62), found no differences in persistence between reference-like variants (including E variants) and non-reference-like (non-E-like) variants.
The current study is to our knowledge the first to assess the persistence of different HPV18 variant lineages, so there are no relevant data for comparison. We observed a higher risk of persistence (i.e. lower clearance rate) of HPV18 As/Ai variants compared to E variants that approximated statistical significance. However, there is no a priori biologic reason we are aware of to have anticipated this result, and the finding should be considered cautiously until confirmed by other independent studies.
Reports of a stronger association of cervical cancer with non-E than E variants of HPV16 is the main reason for the current interest in studying HPV16 lineages (Berumen et al., 2001; Burk et al., 2003
; Hildesheim et al., 2001
; Londesborough et al., 1996
; Zehbe et al., 1998
). However, an association between HPV16 variant lineage and low-grade lesions has not been established. In fact, a recent prospective study found no association between HPV16 variant lineage and risk of any SIL [relative risk=1·0, 95 % confidence interval (CI) 0·42·4], whereas the incidence of cervical intraepithelial neoplasia-2/3 (CIN-2/3) was increased 3·5-fold (95 % CI 1·011·8) in the presence of a non-E versus E variant (Xi et al., 2002
). Consistent with these data, we found no significant differences in the association of SIL (mainly low-grade SIL) with HPV16 variant lineage.
For HPV18, we observed a significantly reduced odds of SIL associated with Af relative to E variants. Unfortunately there are few prior studies of HPV18 variant lineage and cervical neoplasia for purposes of comparison. Villa et al. (2000) reported an increased odds of SIL for non-E HPV18 variants relative to E variants, although they combined these with non-E HPV16 variants in their analyses (i.e. they did not study HPV18 separately), and only one subject harboured an HPV18 Af variant. Burk et al. (2003)
found no association between HPV18 variant lineage and squamous cell carcinoma, but the investigators combined Af and E lineages in their analysis. De Boer et al. (2005)
did not detect any Af variants in adenocarcinoma samples collected from women from three international populations (Surinamese, Dutch and Indonesian), but observed a higher prevalence of Ai/AA variants in adenocarcinoma compared to squamous cell carcinoma. It would therefore seem that individual HPV18 variant lineages might carry differential risks for cervical cancers. Whether the HPV18 Af lineage is associated with a lower risk of cervical neoplasia, as observed in the current study, will need to be confirmed in independent populations.
This investigation has several important limitations that should be noted. First, even though this study represents the largest of its kind in HIV-positive women, we were often limited in our inability to study variants except in a grouped fashion (i.e. by lineage), rather than by specific variant subgroup. Secondly, 75 % of, rather than all, specimens positive for HPV16 or 18 by MY09/MY11 PCR were successfully amplified and sequenced for E6. The E6 amplification rates did not vary by HPV viral load estimated by dot-blot intensity in the MY09/MY11 PCR assay, or by immune status or other parameters, but still we cannot be certain how the excluded samples might have affected our results. Lastly, cytology was used to assess clinical end points. Although both cytology and histology have similar levels of inter-rater agreement, histological confirmation of the presence and grade of neoplasia would have added certainty to the diagnoses.
To summarize, we detected a broad range of HPV16 and 18 variants in our population of HIV-positive and -negative women, and found that a small number of key nucleotide positions could be used to distinguish all of the major variant groups and perhaps even their subgroups. Whether this approach to differentiating HPV16 and 18 variants can be generalized needs to be examined in independent populations. HIV coinfection and immune status (i.e. CD4+ T-cell count and HIV-RNA level) were not associated with the distribution of variants, but we observed a strong association of HPV16 and 18 variant lineage with race. In particular, the data showed that even outside Africa, Af variants of HPV16 and 18 are more common in women of African descent, most likely due to race preference in sexual mixing behaviours and/or immunogenetic factors. We did not, however, find evidence to suggest that non-E HPV16 variant lineage, which is reported to be high-risk for cervical cancer, was associated with either longer time-to-clearance of infection or greater risk of SIL (primarily low-grade SIL). Similarly, no association of HPV18 lineage with time-to-clearance of infection, and only an inverse association between HPV18 Af variants and SIL was observed. Thus, if it is correct that non-E HPV16 or 18 variant lineage is a risk factor for cancer, our data, as well as the results of a previous prospective study of HPV16, collectively suggest that the effect is most likely at a later stage in the multistage process of cervical tumorigenesis (e.g. a greater risk of progression following establishment of cervical neoplasia).
![]() |
ACKNOWLEDGEMENTS |
---|
![]() |
REFERENCES |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Barkan, S. E., Melnick, S. L., Preston-Martin, S. & 7 other authors (1998). The Women's Interagency HIV Study. WIHS Collaborative Study Group. Epidemiology 9, 117125.[CrossRef][Medline]
Bernard, H. U., Chan, S. Y., Manos, M. M., Ong, C. K., Villa, L. L., Delius, H., Peyton, C. L., Bauer, H. M. & Wheeler, C. M. (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, 10771085.[Medline]
Berumen, J., Ordonez, R. M., Lazcano, E. & 7 other authors (2001). Asian-American variants of human papillomavirus 16 and risk for cervical cancer: a case-control study. J Natl Cancer Inst 93, 13251330.
Bontkes, H. J., van Duin, M., de Gruijl, T. D. & 11 other authors (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, 166171.[CrossRef][Medline]
Bosch, F. X., Manos, M. M., Munoz, N. & 7 other authors (1995). Prevalence of human papillomavirus in cervical cancer: a worldwide perspective. International biological study on cervical cancer (IBSCC) Study Group. J Natl Cancer Inst 87, 796802.
Burk, R. D., Terai, M., Gravitt, P. E. & 10 other authors (2003). Distribution of human papillomavirus types 16 and 18 variants in squamous cell carcinomas and adenocarcinomas of the cervix. Cancer Res 63, 72157220.
Casas, L., Galvan, S. C., Ordonez, R. M., Lopez, N., Guido, M. & Berumen, J. (1999). Asian-american variants of human papillomavirus type 16 have extensive mutations in the E2 gene and are highly amplified in cervical carcinomas. Int J Cancer 83, 449455.[CrossRef][Medline]
Chaturvedi, A. K., Brinkman, J. A., Gaffga, A. M., Dumestre, J., Clark, R. A., Braly, P. S., Dunlap, K., Kissinger, P. J. & Hagensee, M. E. (2004). Distribution of human papillomavirus type 16 variants in human immunodeficiency virus type 1-positive and -negative women. J Gen Virol 85, 12371241.
Da Costa, M. M., Hogeboom, C. J., Holly, E. A. & Palefsky, J. M. (2002). Increased risk of high-grade anal neoplasia associated with a human papillomavirus type 16 E6 sequence variant. J Infect Dis 185, 12291237.[CrossRef][Medline]
de Boer, M. A., Peters, L. A., Aziz, M. F., Siregar, B., Cornain, S., Vrede, M. A., Jordanova, E. S., Kolkman-Uljee, S. & Fleuren, G. J. (2004). Human papillomavirus type 16 E6, E7, and L1 variants in cervical cancer in Indonesia, Suriname, and The Netherlands. Gynecol Oncol 94, 488494.[CrossRef][Medline]
de Boer, M. A., Peters, L. A., Aziz, M. F., Siregar, B., Cornain, S., Vrede, M. A., Jordanova, E. S. & Fleuren, G. J. (2005). Human papillomavirus type 18 variants: histopathology and E6/E7 polymorphisms in three countries. Int J Cancer 114, 422425.[CrossRef][Medline]
del Refugio Gonzalez-Losa, M., Laviada Mier y Teran, M. A., Puerto-Solis, M. & Garcia-Carranca, A. (2004). Molecular variants of HPV type 16 E6 among Mexican women with LSIL and invasive cancer. J Clin Virol 29, 9598.[CrossRef][Medline]
Franco, E. L., Villa, L. L., Sobrinho, J. P., Prado, J. M., Rousseau, M. C., Desy, M. & Rohan, T. E. (1999). Epidemiology of acquisition and clearance of cervical human papillomavirus infection in women from a high-risk area for cervical cancer. J Infect Dis 180, 14151423.[CrossRef][Medline]
Gemignani, F., Landi, S., Chabrier, A., Smet, A., Zehbe, I., Canzian, F. & Tommasino, M. (2004). Generation of a DNA microarray for determination of E6 natural variants of human papillomavirus type 16. J Virol Methods 119, 95102.[CrossRef][Medline]
Giannoudis, A. & Herrington, C. S. (2001). Human papillomavirus variants and squamous neoplasia of the cervix. J Pathol 193, 295302.[CrossRef][Medline]
Giuliano, A. R., Harris, R., Sedjo, R. L. & 7 other authors (2002). Incidence, prevalence, and clearance of type-specific human papillomavirus infections: The Young Women's Health Study. J Infect Dis 186, 462469.[CrossRef][Medline]
Hecht, J. L., Kadish, A. S., Jiang, G. & Burk, R. D. (1995). Genetic characterization of the human papillomavirus (HPV) 18 E2 gene in clinical specimens suggests the presence of a subtype with decreased oncogenic potential. Int J Cancer 60, 369376.[Medline]
Hildesheim, A., Schiffman, M., Bromley, C. & 12 other authors (2001). Human papillomavirus type 16 variants and risk of cervical cancer. J Natl Cancer Inst 93, 315318.
Ho, L., Chan, S. Y., Burk, R. D. & 7 other authors (1993). The genetic drift of human papillomavirus type 16 is a means of reconstructing prehistoric viral spread and the movement of ancient human populations. J Virol 67, 64136423.
Ho, G. Y., Bierman, R., Beardsley, L., Chang, C. J. & Burk, R. D. (1998). Natural history of cervicovaginal papillomavirus infection in young women. N Engl J Med 338, 423428.
Hu, X., Pang, T., Guo, Z., Mazurenko, N., Kisseljov, F., Ponten, J. & Nister, M. (2001a). HPV16 E6 gene variations in invasive cervical squamous cell carcinoma and cancer in situ from Russian patients. Br J Cancer 84, 791795.[CrossRef][Medline]
Hu, X., Pang, T., Guo, Z., Ponten, J., Nister, M. & Bernard Afink, G. (2001b). Oncogene lineages of human papillomavirus type 16 E6, E7 and E5 in preinvasive and invasive cervical squamous cell carcinoma. J Pathol 195, 307311.[CrossRef][Medline]
Icenogle, J. P., Laga, M., Miller, D., Manoka, A. T., Tucker, R. A. & Reeves, W. C. (1992). Genotypes and sequence variants of human papillomavirus DNAs from human immunodeficiency virus type 1-infected women with cervical intraepithelial neoplasia. J Infect Dis 166, 12101216.[Medline]
Kammer, C., Tommasino, M., Syrjanen, S., Delius, H., Hebling, U., Warthorst, U., Pfister, H. & Zehbe, I. (2002). Variants of the long control region and the E6 oncogene in European human papillomavirus type 16 isolates: implications for cervical disease. Br J Cancer 86, 269273.[CrossRef][Medline]
Liaw, K. L., Hildesheim, A., Burk, R. D. & 9 other authors (2001). A prospective study of human papillomavirus (HPV) type 16 DNA detection by polymerase chain reaction and its association with acquisition and persistence of other HPV types. J Infect Dis 183, 815.[CrossRef][Medline]
Lizano, M., Berumen, J., Guido, M. C., Casas, L. & Garcia-Carranca, A. (1997). Association between human papillomavirus type 18 variants and histopathology of cervical cancer. J Natl Cancer Inst 89, 12271231.
Londesborough, P., Ho, L., Terry, G., Cuzick, J., Wheeler, C. & Singer, A. (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, 364368.[CrossRef][Medline]
Matsumoto, K., Yoshikawa, H., Nakagawa, S. & 8 other authors (2000). Enhanced oncogenicity of human papillomavirus type 16 (HPV16) variants in Japanese population. Cancer Lett 156, 159165.[CrossRef][Medline]
Mayrand, M. H., Coutlee, F., Hankins, C., Lapointe, N., Forest, P., de Ladurantaye, M. & Roger, M. (2000). Detection of human papillomavirus type 16 DNA in consecutive genital samples does not always represent persistent infection as determined by molecular variant analysis. J Clin Microbiol 38, 33883393.
Moscicki, A. B., Ellenberg, J. H., Farhat, S. & Xu, J. (2004). Persistence of human papillomavirus infection in HIV-infected and -uninfected adolescent girls: risk factors and differences, by phylogenetic type. J Infect Dis 190, 3745.[CrossRef][Medline]
Nindl, I., Rindfleisch, K., Lotz, B., Schneider, A. & Durst, M. (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, 203207.[CrossRef][Medline]
Ong, C. K., Chan, S. Y., Campo, M. S. & 7 other authors (1993). Evolution of human papillomavirus type 18: an ancient phylogenetic root in Africa and intratype diversity reflect coevolution with human ethnic groups. J Virol 67, 64246431.
Palefsky, J. M., Minkoff, H., Kalish, L. A. & 7 other authors (1999). Cervicovaginal human papillomavirus infection in human immunodeficiency virus-1 (HIV)-positive and high-risk HIV-negative women. J Natl Cancer Inst 91, 226236.
Perez-Gallego, L., Moreno-Bueno, G., Sarrio, D., Suarez, A., Gamallo, C. & Palacios, J. (2001). Human papillomavirus-16 E6 variants in cervical squamous intraepithelial lesions from HIV-negative and HIV-positive women. Am J Clin Pathol 116, 143148.[CrossRef][Medline]
Picconi, M. A., Alonio, L. V., Sichero, L., Mbayed, V., Villa, L. L., Gronda, J., Campos, R. & Teyssie, A. (2003). Human papillomavirus type-16 variants in Quechua aboriginals from Argentina. J Med Virol 69, 546552.[CrossRef][Medline]
Richardson, H., Kelsall, G., Tellier, P., Voyer, H., Abrahamowicz, M., Ferenczy, A., Coutlee, F. & Franco, E. L. (2003). The natural history of type-specific human papillomavirus infections in female university students. Cancer Epidemiol Biomarkers Prev 12, 485490.
Rousseau, M. C., Pereira, J. S., Prado, J. C., Villa, L. L., Rohan, T. E. & Franco, E. L. (2001). Cervical coinfection with human papillomavirus (HPV) types as a predictor of acquisition and persistence of HPV infection. J Infect Dis 184, 15081517.[CrossRef][Medline]
Strickler, H. D., Palefsky, J. M., Shah, K. V. & 16 other authors (2003). Human papillomavirus type 16 and immune status in human immunodeficiency virus-seropositive women. J Natl Cancer Inst 95, 10621071.
Terry, G., Ho, L. & Cuzick, J. (1997). Analysis of E2 amino acid variants of human papillomavirus types 16 and 18 and their associations with lesion grade and HLA DR/DQ type. Int J Cancer 73, 651655.[CrossRef][Medline]
van Duin, M., Snijders, P. J., Vossen, M. T., Klaassen, E., Voorhorst, F., Verheijen, R. H., Helmerhorst, T. J., Meijer, C. J. & Walboomers, J. M. (2000). Analysis of human papillomavirus type 16 E6 variants in relation to p53 codon 72 polymorphism genotypes in cervical carcinogenesis. J Gen Virol 81, 317325.
Villa, L. L., Sichero, L., Rahal, P., Caballero, O., Ferenczy, A., Rohan, T. & Franco, E. L. (2000). Molecular variants of human papillomavirus types 16 and 18 preferentially associated with cervical neoplasia. J Gen Virol 81, 29592968.
Walboomers, J. M., Jacobs, M. V., Manos, M. M. & 7 other authors (1999). Human papillomavirus is a necessary cause of invasive cervical cancer worldwide. J Pathol 189, 1219.[CrossRef][Medline]
Wheeler, C. M., Yamada, T., Hildesheim, A. & Jenison, S. A. (1997). Human papillomavirus type 16 sequence variants: identification by E6 and L1 lineage-specific hybridization. J Clin Microbiol 35, 1119.[Abstract]
Xi, L. F., Demers, G. W., Koutsky, L. A., Kiviat, N. B., Kuypers, J., Watts, D. H., Holmes, K. K. & Galloway, D. A. (1995). Analysis of human papillomavirus type 16 variants indicates establishment of persistent infection. J Infect Dis 172, 747755.[Medline]
Xi, L. F., Koutsky, L. A., Galloway, D. A., Kuypers, J., Hughes, J. P., Wheeler, C. M., Holmes, K. K. & Kiviat, N. B. (1997). Genomic variation of human papillomavirus type 16 and risk for high grade cervical intraepithelial neoplasia. J Natl Cancer Inst 89, 796802.
Xi, L. F., Critchlow, C. W., Wheeler, C. M. & 9 other authors (1998). Risk of anal carcinoma in situ in relation to human papillomavirus type 16 variants. Cancer Res 58, 38393844.[Abstract]
Xi, L. F., Carter, J. J., Galloway, D. A., Kuypers, J., Hughes, J. P., Lee, S. K., Adam, D. E., Kiviat, N. B. & Koutsky, L. A. (2002). Acquisition and natural history of human papillomavirus type 16 variant infection among a cohort of female university students. Cancer Epidemiol Biomarkers Prev 11, 343351.
Yamada, T., Wheeler, C. M., Halpern, A. L., Stewart, A. C., Hildesheim, A. & Jenison, S. A. (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, 77437753.[Abstract]
Yamada, T., Manos, M. M., Peto, J., Greer, C. E., Munoz, N., Bosch, F. X. & Wheeler, C. M. (1997). Human papillomavirus type 16 sequence variation in cervical cancers: a worldwide perspective. J Virol 71, 24632472.[Abstract]
Zehbe, I., Wilander, E., Delius, H. & Tommasino, M. (1998). Human papillomavirus 16 E6 variants are more prevalent in invasive cervical carcinoma than the prototype. Cancer Res 58, 829833.[Abstract]
Zehbe, I., Tachezy, R., Mytilineos, J. & 7 other authors (2001a). Human papillomavirus 16 E6 polymorphisms in cervical lesions from different European populations and their correlation with human leukocyte antigen class II haplotypes. Int J Cancer 94, 711716.[CrossRef][Medline]
Zehbe, I., Voglino, G., Wilander, E., Delius, H., Marongiu, A., Edler, L., Klimek, F., Andersson, S. & Tommasino, M. (2001b). p53 codon 72 polymorphism and various human papillomavirus 16 E6 genotypes are risk factors for cervical cancer development. Cancer Res 61, 608611.
Received 24 March 2005;
accepted 16 June 2005.
HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
INT J SYST EVOL MICROBIOL | MICROBIOLOGY | J GEN VIROL |
J MED MICROBIOL | ALL SGM JOURNALS |