Open Access

Co-infection with HPV types from the same species provides natural cross-protection from progression to cervical cancer

Infectious Agents and Cancer20149:26

Received: 4 April 2014

Accepted: 22 July 2014

Published: 12 August 2014



The worldwide administration of bivalent and quadrivalent HPV vaccines has resulted in cross-protection against non-vaccine HPV types. Infection with multiple HPV types may offer similar cross-protection in the natural setting. We hypothesized that infections with two or more HPV types from the same species, and independently, infections with two or more HPV types from different species, associate with protection from high-grade lesions.


We recruited a cohort of 94 HIV, HPV-positive women from Botswana, with Grade 2 or higher cervical intraepithelial neoplasia. Infections with 2 or more HPV types from a single species associated with reduced lesion severity in univariate analysis (OR = 0.41, 95% CI 0.18-0.97, p = 0.042), when adjusted for the presence of HPV 16 or 18 types (OR = 0.41, 95% CI 0.17-1.00, p = 0.049), or all high-risk HPV type infections (OR = 0.38, 95% CI 0.16-0.90, p = 0.028). Infections with 2 or more HPV types from different species did not associate (OR = 0.68, 95% CI 0.25-1.81, p = 0.435).


Our findings show that co-infections with genetically similar HPV types reduce the likelihood of progression to high-grade lesions in HIV positive women, an effect not observed in co-infections with taxonomically different HPV types. This observation is possibly caused by an immune cross-protection through a similar mechanism to that observed after HPV vaccination.


HPV Taxonomy Cervical cancer Coinfection Immune cross-protection Vaccine


Cervical cancer is the third most common cancer in women worldwide and the number one cause of cancer-related mortality in Sub-Saharan Africa [1]. Infection with an oncogenic human papillomavirus (HPV) type is necessary, but not sufficient, for progression to cervical cancer [2]. Most HPV types do not cause cancer, and position on the viral phylogeny affects disease risk [3, 4]. HPV phylogeny is based on genetic distance between isolates, with the genus serving as the broadest taxon. Within a genus, similar viruses are classified into species, and viruses within a species are subdivided into types [5, 6]. Only the alpha genus of HPV is associated with cervical cancer [7]. The carcinogenic HPV types fall into 5 species, α5, α6, α7, α9 and α11 [3] and are termed high-risk species. Two HPV types, 16 and 18, are responsible for approximately 70% of cervical cancer cases worldwide; with the majority of the remainder ascribed to 11 other types [8].

Vaccination against HPV was a major advance in cervical cancer prevention, offering protection against a few disease causing types. Currently, 2 prophylactic vaccines are approved for human use. The tetravalent vaccine covers four HPV types, offering protection against cervical cancer (types 16 and 18) as well as genital warts (types 6 and 11). The bivalent vaccine covers types 16 and 18, targeting only cervical cancer [5, 6, 9]. Understanding the effects of HPV co-infections is essential in determining the long term effects of the vaccines. Availability of such results can better inform the projections on the effect of broad vaccine implementation, and the subsequent prevalence of cervical cancer. Available data suggest that both vaccines have a variable level of cross-protection against types closely related HPV types [10]. It is hypothesized that this protection can be attributed to the genetic similarity between the vaccine types and those against which cross-protection is afforded, as the types for which cross-protection was observed are in the same species as HPV 16 and 18.

In this study we determined whether taxonomic relationships of co-infecting HPV types associated with progression to advanced cervical lesions (CIN3+). We tested two independent hypotheses: 1) infections with two or more HPV types from the same species associate with progression to high-grade lesions; or 2) infections with two or more HPV types from different species associate with progression to high-grade lesions.


We recruited HIV-positive women living in Botswana who were diagnosed with a Grade 2 or higher cervical intraepithelial neoplasia (CIN2 or above) on colposcopy. All patients provided informed consent.

Specimen collection and processing have been previously described [11]. Briefly, Pap smear swabs were lysed in 500 μL of lysis buffer (Roche MagNA Pure LC DNA Isolation Kit) for 30 minutes at room temperature. DNA was extracted, amplified and analyzed following the manufacturer’s specifications (Roche Linear Array® HPV Genotyping test). This method can detect 13 high-risk HPV genotypes (genotypes 16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, 59, and 68) and 24 low-risk HPV genotypes (genotypes 6, 11, 26, 40, 42, 53, 54, 55, 61, 62, 64, 66, 67, 69, 70, 71, 72, 73, 81, 82, 83, 84, IS39, and CP6108).

HPV types were assigned to species using conventional taxonomic criteria [5, 6]. Briefly, using the L1 ORF of papillomaviruses, viruses with <60% sequence homology in this region are considered to be in a different genus. Sequence homology between 60 and 70% constitutes different species, while different HPV types within a species have between 70 and 89% nucleotide identity [5, 6]. Only HPV types in the alpha-papillomavirus genus were considered in our study.

Univariate logistic regression was used to determine association between CIN3 lesions (as opposed to CIN2) and three categories of HPV co-infection, separately. The HPV co-infection categories were: 1) two or more HPV types from two or more species, 2) two or more HPV types from a single species, 3) two or more HPV types from a single high-risk species (α5, α6, α7, α9 or α11 [3]).

We also tested the association between CIN3 lesions and variables previously associated with CIN progression; namely: baseline CD4 count, patient age at the time of screening, patient age at the time of first intercourse, and lifetime number of sexual partners [12]. Univariate logistic regression was performed between presence of CIN3 and each variable. Univariate analyses were used to separately model the three categories of patient HPV as the dependent variable and patient CD4 counts either at presentation or the last available CD4 value as the independent variables.

Multivariate logistic regression models of CIN3 were tested, separately using the 3 HPV co-infection categories above, along with other variables statistically significant in univariate analyses. Presence of HPV 16 or 18, or the presence of a high-risk HPV of any type was also included in the model.

Study population

94 women with CIN2, CIN3, and squamous cell carcinoma made up the final study population. 41 of the patients had CIN2, 52 had CIN3, and one had squamous cell carcinoma (SCC). The squamous cell carcinoma patient was categorized as CIN3 in the logistic regression analyses (ordinal analyses with a separate SCC category yielded comparable results, Additional file 1: Table S1). The mean age of study participants was 36.78 and the median number of HPV type co-infections per patient was 3 (Table 1).

Table 1

Summary statistics of study population


Mean (st dev)

Median (IQR)


36.78 (5.90)

36 (33-39)

HPV infections per person

3.93 (2.62)

3 (2-5)

Age at first intercourse*

18.13 (2.32)

18 (16-20)

Number of lifetime sexual partners*

5.86 (5.26)

4 (3-6)

CD4 count during first visit**

254.85 (207.29)

200 (109-365)

CD4 count during last visit

452.43 (208.12)

426 (285-578)

Data available for 63 patients* and for 61 patients**.

The classification of HPV types into species is presented in Additional file 1: Table S2 [13]. 79 patients were infected with multiple HPV types. 72 patients had two or more HPV types from different species, 53 patients had two or more HPV types from a single species, and 40 patients had two or more HPV types from a single high-risk species. We found co-infections with only single HPV types of multiple species in 26 patients. Seven patients were infected with multiple HPV types from one species exclusively, all of which were high-risk species. 46 patients had infections with both HPV types in multiple species and multiple HPV types in the same species. 90 patients were infected with a high-risk HPV type, while 50 patients had an infection with either HPV 16 or 18.

Univariate analyses

Infection with two or more HPV types from a single species protected against progression to higher risk lesions (OR = 0.41, 95% CI 0.18-0.97, p = 0.042; Table 2). Infection with two or more HPV from a single high-risk species was also protective (OR = 0.37, 95% CI 0.16-0.86, p = 0.021). Infection with two or more HPV types from different species did not associate with progression (OR = 0.68, 95% CI 0.25-1.81, p = 0.435).
Table 2

Univariate and multivariate logistic regression results modeling CIN3 presence/absence or HPV co-infections



Odds ratio

95% CI

p value

Presence of CIN3

Multiple types single species




Presence of CIN3

Multiple types different species




Presence of CIN3

Multiple types single high-risk species




Presence of CIN3

HPV 16 or 18




Presence of CIN3





Presence of CIN3

Sexual Partners*




Presence of CIN3

Age at First Sex*




Presence of CIN3

CD4 count at presentationϮ




Presence of CIN3

Last available CD4 count




Presence of CIN3

Multiple types from single species; high-risk HPVǂ




Presence of CIN3

Multiple types from single species; HPV 16 or 18




Presence of CIN3

Multiple types from different species; high-risk HPVǂ




Presence of CIN3

Multiple types from different species; HPV 16 or 18




Presence of CIN3

Multiple types from single high-risk species; high-risk HPVǂ




Presence of CIN3

Multiple types from single high-risk species; HPV 16 or 18




2+ HPV infections

CD4 count at presentation




2+ HPV infections

Last available CD4 count




Multiple types single species

CD4 count at presentation




Multiple types single species

Last available CD4 count




Multiple types single high-risk species

CD4 count at presentation




Multiple types single high-risk species

Last available CD4 count




Multiple types different species

CD4 count at presentation




Multiple types different species

Last available CD4 count




Data available for 63 patients* and for 61 patientsϮ. 90 patients had high risk HPV typesǂ.

Univariate regression of CIN3 with a patient’s CD4 count at first presentation, last available CD4 count, patient’s age, number of lifetime sexual partners, and most of the individual high-risk HPV types were not statistically significant (Table 2, Additional file 1: Table S3). Patients with HPV types 39 were always diagnosed with CIN3, and those with HPV 26, 40, or 73 always had CIN2 (Additional file 1: Table S3). HPV 51 also associated with CIN3 (OR = 0.16, 95% CI 0.03-0.81, p = 0.027) (Additional file 1: Table S3). In our sample, infection with either HPV 16 or 18 was not associated with progression to CIN3 (p = 0.62). Initial and last available CD4 counts were not associated with the presence of any of the HPV co-infection categories (Table 2).

Multivariate analyses

In multivariate analysis, the association of CIN3 with two or more HPV type infections from multiple species was not significant when adjusting for high-risk HPV type infections (OR = 0.61, 95% CI 0.22-1.63, p = 0.322; Table 2). Association of CIN3 with two or more HPV type infections from a single species was significant when adjusting for high-risk HPV type infections (single species co-infection OR = 0.38, 95% CI 0.16-0.90, p = 0.028; Table 2). Infections with two or more HPV types from a single high-risk species were significant in multivariate analyses adjusted for all high-risk HPV (OR = 0.42, 95% CI 0.18-0.97, p = 0.044; Table 2).


Our results indicate that the presence of infections with two or more HPV types from the same species protects from progression to high-grade CIN. While our data did not longitudinally ascertain the sequence of infection events, it is likely that some portion of our patients were infected with a low-risk HPV type prior to the acquisition of a high-risk type from the same species. This sequence of events could boost the immune system against the high-risk types and provide cross-protection via a mechanism similar to that afforded for HPV 31 and 45 infections when using vaccines [14].

Cross-reactive immune response based on taxonomic placement has been previously considered among HPV vaccines [15]. Immune responses to natural infection have been shown to be weaker than those elicited by the vaccines [14]. However, the persistence of HPV infections coupled with the immunodominant nature of the response, generating antibodies to only a subset of antigens, makes natural cross-protection an important factor to consider [16]. Cross-reactivity could offer a degree of protection capable of preventing or delaying carcinogenesis to only a subset of similar type infections (same species) but not strong enough to clear the infection. This mechanism could account for the effects observed in our study. This conclusion is further supported by the result that co-infections with HPV types from different species did not associate with protection for progressing to CIN3, while univariate analyses of same species infections in only the five high-risk species provided stronger associations and effect sizes than same species infections in general (Table 2, Additional file 1: Table S1). Hence sequence similarity among co-infections increases protection.

HIV and resultant low CD4 counts have been previously associated with increased risk of single and multiple HPV infections [17]. However, in our study, HIV-associated immunosuppression was not predictive of multiple HPV infections (Table 2). We observed no effect of CD4 counts on the likelihood of acquisition of co-infections with multiple HPV species, regardless of category (Table 2). This leads us to suggest that our findings are applicable beyond the niche of HIV + patients.

Our data, observed in an unvaccinated, natural setting, indicate that protection from any high-risk species type is likely to confer some degree of protection against all types in that species, supporting the argument that vaccines would do better if all high-risk species were covered. Therefore, widespread use of HPV vaccines can indirectly affect the prevalence of some non-vaccine types. Since the bivalent vaccine covers species α7 and α9, and the quadrivalent adds protection to α10, the burden of HPV carcinogenesis in a fully vaccinated population might shift to the high-risk species that are not currently covered by either vaccine; namely α5, α6, and α11. Adding virus-like particles from at least one HPV type in each of these species should therefore be a priority for future vaccine development.

This study was reviewed and approved by the Institutional Review Boards of the University of Pennsylvania and the Botswana Ministry of Health.



RSS was supported by Public Health Service award T32 GM07347 from the National Institute of General Medical Studies for the Vanderbilt Medical-Scientist Training Program. SMW and RSS were partially supported by NIH grant P20 GM103534. NMZ was supported in part by NIH grants R01AI097045 and P30AI45008 (Penn Center for AIDS Research). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Authors’ Affiliations

Center for Human Genetics Research, Department of Molecular Physiology and Biophysics, Vanderbilt University Medical Center
Department of Genetics, Geisel School of Medicine, Dartmouth College
Division of Infectious Disease, University of Pennsylvania
Botswana-University of Pennsylvania Partnership
Department of Medicine, University of Botswana
Center for Human Genetics Research, Vanderbilt University


  1. Ferlay J, Shin HR, Bray F, Forman D, Mathers C, Parkin DM: Estimates of worldwide burden of cancer in 2008: GLOBOCAN 2008. Int J Canc J Int Canc. 2010, 127 (12): 2893-2917.View ArticleGoogle Scholar
  2. Walboomers JM, Jacobs MV, Manos MM, Bosch FX, Kummer JA, Shah KV, Snijders PJ, Peto J, Meijer CJ, Munoz N: Human papillomavirus is a necessary cause of invasive cervical cancer worldwide. J Pathol. 1999, 189 (1): 12-19.PubMedView ArticleGoogle Scholar
  3. Schiffman M, Rodriguez AC, Chen Z, Wacholder S, Herrero R, Hildesheim A, Desalle R, Befano B, Yu K, Safaeian M: A population-based prospective study of carcinogenic human papillomavirus variant lineages, viral persistence, and cervical neoplasia. Cancer Res. 2010, 70 (8): 3159-3169.PubMedPubMed CentralView ArticleGoogle Scholar
  4. Burk RD, Chen Z, Van Doorslaer K: Human papillomaviruses: genetic basis of carcinogenicity. Publ Health Genomics. 2009, 12 (5–6): 281-290.View ArticleGoogle Scholar
  5. de Villiers EM, Fauquet C, Broker TR, Bernard HU, zur Hausen H: Classification of papillomaviruses. Virology. 2004, 324 (1): 17-27.PubMedView ArticleGoogle Scholar
  6. Bernard HU, Burk RD, Chen Z, van Doorslaer K, zur Hausen H, de Villiers EM: Classification of papillomaviruses (PVs) based on 189 PV types and proposal of taxonomic amendments. Virology. 2010, 401 (1): 70-79.PubMedPubMed CentralView ArticleGoogle Scholar
  7. Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52: 696-704.PubMedView ArticleGoogle Scholar
  8. Doorslaer KV: Evolution of the Papillomaviridae. Virology. 2013, 445: 11-20.PubMedView ArticleGoogle Scholar
  9. Schiffman M, Clifford G, Buonaguro FM: Classification of weakly carcinogenic human papillomavirus types: addressing the limits of epidemiology at the borderline. Infect Agents Canc. 2009, 4: 8-View ArticleGoogle Scholar
  10. Clifford G, Franceschi S, Diaz M, Munoz N, Villa LL: Chapter 3: HPV type distribution in women with and without cervical neoplastic diseases. Vaccine. 2006, 24 (Suppl 3): 26-34.View ArticleGoogle Scholar
  11. Barr E, Tamms G: Quadrivalent human papillomavirus vaccine. Clin Infect Dis: Offic Publ Infect Dis Soc Am. 2007, 45 (5): 607-609.Google Scholar
  12. Jaquet A, Horo A, Ekouevi DK, Toure B, Coffie PA, Effi B, Lenaud S, Messou E, Minga A, Sasco AJ: Risk factors for cervical intraepithelial neoplasia in HIV-infected women on antiretroviral treatment in Côte d'Ivoire. West Africa PLoS One. 2014, 9 (3): e90625-PubMedView ArticleGoogle Scholar
  13. Keam SJ, Harper DM: Human papillomavirus types 16 and 18 vaccine (recombinant, AS04 adjuvanted, adsorbed) [Cervarix]. Drugs. 2008, 68 (3): 359-372.PubMedView ArticleGoogle Scholar
  14. Mariani L, Venuti A: HPV vaccine: an overview of immune response, clinical protection, and new approaches for the future. J Transl Med. 2010, 8: 105-PubMedPubMed CentralView ArticleGoogle Scholar
  15. Malagon T, Drolet M, Boily MC, Franco EL, Jit M, Brisson J, Brisson M: Cross-protective efficacy of two human papillomavirus vaccines: a systematic review and meta-analysis. Lancet Infect Dis. 2012, 12 (10): 781-789.PubMedView ArticleGoogle Scholar
  16. Siddiqui MA, Perry CM: Human papillomavirus quadrivalent (types 6, 11, 16, 18) recombinant vaccine (Gardasil). Drugs. 2006, 66 (9): 1263-1271. discussion 1272-1263PubMedView ArticleGoogle Scholar
  17. Camargo M, Soto-De Leon SC, Munoz M, Sanchez R, Peňa-Herrera D, Pineda-Peňa AC, Sussmann O, Paez C, Perez-Prados A, Patarroyo ME: Human papillomavirus detection in women with and without human immunodeficiency virus infection in Colombia. BMC Cancer. 2014, 14 (1): 451-PubMedPubMed CentralView ArticleGoogle Scholar


© Sobota et al.; licensee BioMed Central Ltd. 2014

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.