ORIGINAL REPORTS
Cancer Prevention and Control
Article Tools
OPTIONS & TOOLS
COMPANION ARTICLES
ARTICLE CITATION
DOI: 10.1200/GO.20.00286 JCO Global Oncology no. 6 (2020) 1276-1281. Published online August 12, 2020.
PMID: 32783640
Prevalence and Subtype Distribution of High-Risk Human Papillomavirus Among Women Presenting for Cervical Cancer Screening at Karanda Mission Hospital
2 University of Toronto, Toronto, Ontario, Canada
3 Beth Israel Deaconess Medical Center, Boston, MA
4 Massachusetts General Hospital, Boston, MA
5 African Institute of Biomedical Science and Technologies Laboratory, Harare, Zimbabwe
6 Beth Israel Deaconess Medical Center, Boston, MA
P.T., R.P., and D.P. are co-first authors of this article.
High-risk human papillomaviruses (hrHPV) are the primary cause of cervical cancer. Human papillomavirus (HPV) vaccination is expected to prevent cervical cancers caused by the HPV types included in vaccines and possibly by cross-protection from other types. This study sought to determine the hrHPV type distribution in women at a rural Zimbabwe hospital.
We implemented a cross-sectional study at the Karanda Mission Hospital. Using the Visual Inspection with Acetic Acid Cervicography technique, clinicians collected cervical swabs from 400 women presenting for screening for cervical cancer. Samples were initially analyzed by Cepheid GeneXpert; candidate hrHPV genotypes were further characterized using the Anyplex II HPV28 Detection Kit.
Twenty-one percent of the 400 women were positive for a high-risk genotype when using the GeneXpert analyzer; 17% were positive when using the multiplex analysis. Almost two thirds of the hrHPV women had a single DNA type identified, whereas one third had multiple genotypes, ranging from 2 to 5. hrHPV was observed more frequently in HIV-positive than in HIV-negative women (27% v 15%). Of the 113 isolates obtained, 77% were hrHPV genotypes not included in the bivalent or quadrivalent vaccines, and 47% represented DNA types not covered in the nonavalent vaccine. Forty-seven percent of the women with hrHPV harbored a single genotype that was not covered by the nonavalent vaccine.
A large fraction of hrHPV isolates from women participating in a cervical cancer screening program in northern Zimbabwe are DNA types not covered by the bivalent, quadrivalent, or nonavalent vaccines. These findings suggest the importance of characterizing the hrHPV DNA types isolated from cervical neoplasia in this population and determining whether cross-immunization against these genotypes develops after administration of the vaccines in current use.
Globally, cervical cancer is the fourth most common cancer and the fourth-leading cause of cancer death in women. 1 Women in low- and middle-income countries are disproportionately affected, accounting for 85% of all cervical cancer–related deaths. 2 In these countries, poor survival rates are associated with poor infrastructure, lack of trained health care providers and available treatment options, and high costs of treatment. 3-5
-
Key Objective
-
What fraction of high-risk human papillomavirus (hrHPV) isolates from women participating in a cervical cancer screening program in Zimbabwe represent DNA types not covered by the bivalent, quadrivalent, or nonavalent vaccines currently in use?
-
Knowledge Generated
-
Seventy-seven percent of the hrHPV genotypes were not covered by the bivalent or quadrivalent vaccines, and 47% of the hrHPV genotypes were not covered by the nonavalent vaccines. Forty-seven percent of women with hrHPV harbored a single genotype that was not covered by the nonavalent vaccine.
-
Relevance
-
The possibility that oncogenic genotypes are not being neutralized by the vaccines in use presents a significant theoretical challenge to the effectiveness of a vaccination program; more research is required to determine the extent to which these vaccines generate cross-protection against these specific hrHPV subtypes.
Zimbabwe ranks fourth in age-standardized cervical cancer incidence rates. 6 Cervical cancer is the most common cancer among women in Zimbabwe. 7 The estimated incidence is 52.1 per 100,000 women, and the age-adjusted mortality rate is 43.1% 7 ; actual rates are likely much higher than these, given underreporting in rural areas. 8,9
Human papillomavirus (HPV) is one of the most prevalent sexually transmitted infections worldwide. High-risk or oncogenic strains are associated with the development of premalignant and malignant epithelial lesions of the female and male anogenital tract. 10,11 More than 95% of cervical cancers are caused by infection of the cervical epithelium with high-risk human papillomavirus (hrHPV). 12 HPV is also considered a comorbid opportunistic infection in the setting of HIV infection. 13 Despite successful antiretroviral therapy, some authors have described rising HPV infections among HIV-positive individuals. 14
Of central importance to Zimbabwe’s cancer strategy is the prevention and early detection of this disease. Supported financially by Gavi, the Vaccine Alliance, the Ministry of Health and Child Care Zimbabwe sponsored pilot projects that delivered HPV vaccinations to girls 10-14 years of age, using the bivalent HPV vaccine. 15 The quadrivalent HPV vaccine has been used in other settings. 16 The bivalent vaccine confers protection against HPV 16 and 18, and the quadrivalent against HPV 6, 11, 16, and 18 (the first 2 are understood to be nononcogenic). HPV types 16 and 18 are reported to be responsible for approximately 70% of cervical cancers in women. 12,17 These vaccines have been shown to be highly effective in preventing high-grade cervical lesions. 18 A nonavalent vaccine is also available and covers an additional 5 oncogenic subtypes of HPV, including HPV 31, 33, 45, 52, and 58, and has been associated with additional reductions in the incidence of cervical intraepithelial neoplasia. 19
The prevalence of hrHPV genotypes seems to vary geographically. A study consisting of HIV-negative women in Zimbabwe showed that HPV 58 was more common than HPV 16 or 18, 20 whereas a study of women with invasive cervical cancer in Harare showed that HPV 16 was the most prevalent high-risk genotype, followed by HPV 33. 21 In a study conducted in rural Zimbabwe, HPV 35 was found to be the most common genotype, followed by HPV 33 and HPV 58. 22 The possibility of oncogene genotypes incompletely or not at all neutralized by the bivalent, quadrivalent, or nonavalent HPV vaccines presents a significant theoretical challenge to the effectiveness of a vaccination program. 9
In a coordinated effort with the nation’s cancer control strategy, the Karanda Mission Hospital (KMH) initiated a screening program using Visual Inspection with Acetic acid and Cinematography (VIAC). The current study was launched to establish the prevalence of hrHPV in women presenting for screening, to characterize the distribution of genotypes in rural northern Zimbabwe, and to identify and quantify hrHPV DNA types not represented in the vaccines being deployed throughout the nation.
We conducted a prospective cross-sectional study of women presenting for cervical cancer screening (VIAC) to KMH, a community hospital located in the Mt Darwin district of rural Zimbabwe. The study was approved by the Joint Research Ethics Committee for the University of Zimbabwe and the Medical Research Council for Zimbabwe (Prevalence and Subtype Distribution of Cervical High-Risk Human Papillomavirus Among Women Presenting for Cervical Cancer Screening at Karanda Mission Hospital, protocol proposal submitted to Medical Research Council Zimbabwe). Participants were recruited to a total of 400 participants.
This sample size was chosen on the basis of the assumptions that, in the Mt Darwin region, the HPV prevalence in HIV-negative and HIV-positive individuals is 5% and 20%, respectively. The HIV prevalence was estimated at 20% (the current HIV prevalence rate among VIAC patients at KMH).
The study consisted of sexually active female patients 30 to 65 years of age, living in the Mt Darwin district, and presenting to KMH for screening as part of the VIAC program. Eligible participants were invited to participate in the study after they were given a thorough description of the study’s goals, and informed consent was obtained from all participants. Patients were excluded if they had a history of cervical cancer or precancerous lesions, had a history of a hysterectomy, or had never engaged in sexual intercourse. If HIV status was not documented or was unknown, an HIV counselor offered rapid HIV testing and counseling. Nurses used a questionnaire to gather data regarding the participants’ age, parity, contraception use, and HIV status.
Nurses were trained in obtaining a sample for HPV analysis. Cervical brushes were used to collect specimens that were then dropped into a collection vial with 3 mL ThinPrep PreservCyt liquid media. Specimens were tested in an onsite laboratory using the GeneXpert HPV assay. Positive specimens were transported to the African Institute of Biomedical Science and Technologies Laboratory in Harare, where DNA extraction and subtyping were performed using multiplex analysis. Briefly, DNA was extracted using the QIAamp MiniElute Virus Spin kit per the manufacturer’s instructions. HPV genotyping was performed using the Anyplex II HPV28 Detection kit (Seegene, Seoul, Republic of Korea) on a BioRad CFX-96 real-time thermocycler per the manufacturer’s instructions. Twenty-eight HPV genotypes (19 high risk and 9 low risk) were tested in each sample. The 19 high-risk genotypes tested were HPV 16, 18, 26, 31, 33, 35, 39, 45, 51, 52, 53, 56, 58, 59, 66, 68, 69, 73, and 82. The 9 low-risk genotypes tested were HPV 6, 11, 40, 42, 43, 44, 54, 61, and 70. Results of both VIAC and HPV testing were communicated to participants.
Statistical analysis was performed using SAS 9.4. Type-specific hrHPV prevalence was calculated with an exact binomial distribution–based 95% CI. Overall hrHPV prevalence was determined in the HIV-positive and HIV-negative groups.
A total of 400 women were recruited into the study and elected to undergo screening for cervical cancer with VIAC at the KMH. Table 1 summarizes the baseline characteristics of the participants.
|
Among the entire cohort of women, 21% (83 of 400) were positive for a high-risk genotype when using the GeneXpert analyzer. Subsequently, 17% (67 of 400) were positive when subjected to multiplex analysis. The data that follow are derived from the 67 samples analyzed by the multiplex technique.
Sixty-four percent of the hrHPV-positive women had a single virus type isolated, whereas approximately one third were infected by multiple DNA types ranging from 2 to 5 per sample (Fig 1). A total of 113 hrHPV isolates were analyzed from the cohort.
Among the cohort of sampled women, 18.5% (74 of 400) were HIV positive. Of the 67 hrHPV-positive individuals, 18 (27%) were HIV positive. The distribution of DNA types in the total population of hrHPV-positive women and the distribution in the hrHPV HIV-positive women are demonstrated in Figure 2.

FIG 2 Distribution (ordinate) of DNA types (abscissa) in the total population of high-risk human papillomavirus (hrHPV)–positive women and the distribution in the hrHPV HIV-positive population: hrHPV in HIV-negative and HIV-positive women (blue); hrHPV in HIV-positive women (red). HPV, human papillomavirus.
Among the 113 isolates obtained from 67 women, 87 (77%) represented hrHPV DNA types that were not covered in the quadrivalent vaccine (types 6, 11, 16, and 18), as listed in Table 2. Of the 113 isolates, 53 (47%) represented DNA types not covered in the nonavalent vaccine (types 6, 11, 16, 18, 31, 33, 45, 52, and 58), as listed in Table 3.
|
|
Among the genotypes not included in the 9-valent vaccine, this study identified DNA types 35, 39, 51, 53, 56, 66, 68, 69, 73. Fourteen (47%) of the 30 women harboring hrHPV not covered by the 9-valent vaccine demonstrated a single DNA type as the only virus identified.
The results of our study demonstrate a high hrHPV prevalence in rural Zimbabwe, as well as a high proportion of hrHPV genotypes not covered by the quadrivalent or nonavalent vaccines. In our study, hrHPV prevalence by the GeneXpert HPV assay was 21%. A recent community-based cross-sectional study in rural Zimbabwe using GeneXpert HPV found a similar HPV prevalence of 17%. 23 These rates are higher than those previously seen in the United States and Europe. 24
We found that hrHPV infection was more common among HIV-positive women than among HIV-negative women (28% of HIV-positive women v 15% of HIV-negative women), which is consistent with the findings of prior studies. 14,25,26 HIV-positive women in sub-Saharan Africa have a higher prevalence of hrHPV infection and a greater proportion of infections with multiple hrHPV genotypes. 27,28 These findings emphasize the importance of screening for cervical cancer in all women, with particular urgency for those who are HIV positive.
Our data also highlighted the representation of hrHPV genotypes other than HPV 16 and 18 in the Mt Darwin district in rural Zimbabwe. These include HPV 35, 45, 52, 58, and 68. Fitzpatrick et al 23 also found high levels of these types in the Hurungwe district of Zimbabwe. Supporting this point is a prior study that demonstrated high levels of types other than HPV 16 and 18 in rural Zimbabwe. 22 Of the 113 hrHPV types identified among women in our study, 15 were HPV 35 and 7 were HPV 68. Both HPV 35 and HPV 68 are not included in any hrHPV vaccine. In the aggregate, 77% of the identified hrHPV DNA types are not covered by the quadrivalent vaccine, and 47% are not covered by the nonavalent vaccine. Considering the frequency with which women are infected with non–vaccine-covered genotypes, it is important to determine whether the identified DNA types are oncogenic in the population under study.
The patients in this study were women who presented for cervical cancer screening. Among these, 15 of the 67 hrHPV-positive women were VIAC positive or harbored lesions specific for cancer. Of interest, 5 women (33%) harbored a genotype not covered by the available vaccines. This suggests the importance of additional research to determine the HPV genotypes in preneoplastic and frank invasive cancer in this population. This information, if consistent with the findings in this study, will encourage the development of vaccines likely to provide optimal protection.
It is important to consider the possibility of cross-immunization, the process by which antibodies to the hrHPV virus-like particles used in the vaccine are able to neutralize virions of related HPV types and thereby prevent disease. 28,29 Harper et al 30 demonstrated partial cross-protection of the bivalent HPV vaccine against HPV 31 and 45. Bivalent vaccination programs have been shown to decrease the population prevalence of HPV genotypes not covered by the vaccine; however, the same study identified low levels of cross-protection with high-risk DNA types other than HPV 31, 33, and 45. 31 Prior studies from Brown et al 32 and Wheeler et al 29 showed low-to-moderate cross-protection against hrHPV types related to HPV 16 and 18. Herrero 33 analyzed the major findings in these 2 studies and noted that the only significant cross-protection against cervical intraepithelial neoplasia 2-3 or adenocarcinoma in situ lesions was protection against HPV 31, which is most related to HPV 16. Furthermore, Wheeler et al 29 performed an intention-to-treat analysis of women who received at least 1 dose of vaccine and returned for follow-up care, regardless of infection status at the start of the study. Among this “catch-up vaccination” population, there was no significant protection against advanced lesions with any combination of HPV types. 33 Data from 2 randomized controlled trials showed no statistically significant cross-protection against persistent infection with HPV 52 or HPV 58. 34,35
The Karanda and Chidamoyo studies indicated that a substantial number of DNA types are not known to be neutralized by any of the vaccines in current use. Cross-protection against the uncovered genotypes identified in rural Zimbabwean women in this study cohort has not been investigated. Additional research is needed to elucidate the frequency with which uncovered high-risk genotypes contribute to the development of cervical intraepithelial neoplasia or cancer, as well as the extent to which vaccines generate cross-protection against these hrHPV subtypes.
Conception and design: Paul Thistle, Hang Lee, Lowell E. Schnipper
Financial support: Lowell E. Schnipper
Collection and assembly of data: Paul Thistle, Rabea Parpia, Justen Manasa, Lowell E. Schnipper
Data analysis and interpretation: All authors
Manuscript writing: All authors
Final approval of manuscript: All authors
Accountable for all aspects of the work: All authors
The following represents disclosure information provided by authors of this manuscript. All relationships are considered compensated unless otherwise noted. Relationships are self-held unless noted. I = Immediate Family Member, Inst = My Institution. Relationships may not relate to the subject matter of this manuscript. For more information about ASCO's conflict of interest policy, please refer to www.asco.org/rwc or ascopubs.org/go/site/misc/authors.html.
Open Payments is a public database containing information reported by companies about payments made to US-licensed physicians (Open Payments).
Research Funding: Gilead Sciences
Patents, Royalties, Other Intellectual Property: I am the co-editor-in-chief of Oncology UpToDate and receive royalties from the company for this responsibility
No other potential conflicts of interest were reported.
ACKNOWLEDGMENT
We acknowledge the excellent contributions made by the following individuals in facilitating the care of patients participating in this clinical study: Z. Chikove (laboratory scientist), T. Murato (laboratory technician), E. Muchenje (facilitator), L. Hungwe (VIAC nurse), L. Chirindo (VIAC nurse), and C. Tafirenyika (VIAC assistant).
| 1. | Arbyn M , Weiderpass E , Bruni L , et al : Estimates of incidence and mortality of cervical cancer in 2018: A worldwide analysis . Lancet Glob Health 8 : e191 - e203 , 2020 Crossref, Medline, Google Scholar |
| 2. | Randall TC , Ghebre R : Challenges in prevention and care delivery for women with cervical cancer in sub-Saharan Africa . Front Oncol 6 : 160 , 2016 Crossref, Medline, Google Scholar |
| 3. | Ajenifuja KO , Gage JC , Adepiti AC , et al : A population-based study of visual inspection with acetic acid (VIA) for cervical screening in rural Nigeria . Int J Gynecol Cancer 23 : 507 - 512 , 2013 Medline, Google Scholar |
| 4. | LaVigne AW , Triedman SA , Randall TC , et al : Cervical cancer in low and middle income countries: Addressing barriers to radiotherapy delivery . Gynecol Oncol Rep 22 : 16 - 20 , 2017 Medline, Google Scholar |
| 5. | Sankaranarayanan R , Ferlay J : Worldwide burden of gynaecological cancer: The size of the problem . Best Pract Res Clin Obstet Gynaecol 20 : 207 - 225 , 2006 Crossref, Medline, Google Scholar |
| 6. | World Cancer Research Fund : Cervical cancer statistics . https://www.wcrf.org/dietandcancer/cancer-trends/cervical-cancer-statistics Google Scholar |
| 7. | Chokunonga E , Borok MZ , Chirenje ZM , et al : Zimbabwe National Cancer Registry: 2014 Annual Report—Pattern of cancer in Zimbabwe. 2015 . www.uicc.org/membership/zimbabwe-national-cancer-registry Google Scholar |
| 8. | Kuguyo O , Matimba A , Tsikai N , et al : Cervical cancer in Zimbabwe: A situation analysis . Pan Afr Med J 27 : 215 , 2017 Medline, Google Scholar |
| 9. | Chin’ombe N , Sebata NL , Ruhanya V , et al : Human papillomavirus genotypes in cervical cancer and vaccination challenges in Zimbabwe . Infect Agent Cancer 9 : 16 , 2014 Medline, Google Scholar |
| 10. | Bosch FX , Lorincz A , Muñoz N , et al : The causal relation between human papillomavirus and cervical cancer . J Clin Pathol 55 : 244 - 265 , 2002 Crossref, Medline, Google Scholar |
| 11. | Palefsky JM , Rubin M : The epidemiology of anal human papillomavirus and related neoplasia . Obstet Gynecol Clin North Am 36 : 187 - 200 , 2009 Medline, Google Scholar |
| 12. | Muñoz N , Bosch FX , de Sanjosé S , et al : Epidemiologic classification of human papillomavirus types associated with cervical cancer . N Engl J Med 348 : 518 - 527 , 2003 Crossref, Medline, Google Scholar |
| 13. | Palefsky J : Human papillomavirus-related disease in people with HIV . Curr Opin HIV AIDS 4 : 52 - 56 , 2009 Crossref, Medline, Google Scholar |
| 14. | Reusser NM , Downing C , Guidry J , et al : HPV carcinomas in immunocompromised patients . J Clin Med 4 : 260 - 281 , 2015 Medline, Google Scholar |
| 15. | Hidle A , Gwati G , Abimbola T , et al : Cost of a human papillomavirus vaccination project, Zimbabwe . Bull World Health Organ 96 : 834 - 842 , 2018 Medline, Google Scholar |
| 16. | Money DM , Moses E , Blitz S , et al : HIV viral suppression results in higher antibody responses in HIV-positive women vaccinated with the quadrivalent human papillomavirus vaccine . Vaccine 34 : 4799 - 4806 , 2016 Medline, Google Scholar |
| 17. | Bailey HH , Chuang LT , duPont NC , et al : American Society of Clinical Oncology statement: Human papillomavirus vaccination for cancer prevention . J Clin Oncol 34 : 1803 - 1812 , 2016 Link, Google Scholar |
| 18. | FUTURE II Study Group : Quadrivalent vaccine against human papillomavirus to prevent high-grade cervical lesions . N Engl J Med 356 : 1915 - 1927 , 2007 Crossref, Medline, Google Scholar |
| 19. | Yang DY , Bracken K : Update on the new 9-valent vaccine for human papillomavirus prevention . Can Fam Physician 62 : 399 - 402 , 2016 Medline, Google Scholar |
| 20. | Fukuchi E , Sawaya GF , Chirenje M , et al : Cervical human papillomavirus incidence and persistence in a cohort of HIV-negative women in Zimbabwe . Sex Transm Dis 36 : 305 - 311 , 2009 Medline, Google Scholar |
| 21. | Stanczuk GA , Kay P , Sibanda E , et al : Typing of human papillomavirus in Zimbabwean patients with invasive cancer of the uterine cervix . Acta Obstet Gynecol Scand 82 : 762 - 766 , 2003 Medline, Google Scholar |
| 22. | Baay MF , Kjetland EF , Ndhlovu PD , et al : Human papillomavirus in a rural community in Zimbabwe: The impact of HIV co-infection on HPV genotype distribution . J Med Virol 73 : 481 - 485 , 2004 Medline, Google Scholar |
| 23. | Fitzpatrick MB , Dube Mandishora RS , Katzenstein DA , et al : hrHPV prevalence and type distribution in rural Zimbabwe: A community-based self-collection study using near-point-of-care GeneXpert HPV testing . Int J Infect Dis 82 : 21 - 29 , 2019 Medline, Google Scholar |
| 24. | Franceschi S , Herrero R , Clifford GM , et al : Variations in the age-specific curves of human papillomavirus prevalence in women worldwide . Int J Cancer 119 : 2677 - 2684 , 2006 Medline, Google Scholar |
| 25. | Sun XW , Kuhn L , Ellerbrock TV , et al : Human papillomavirus infection in women infected with the human immunodeficiency virus . N Engl J Med 337 : 1343 - 1349 , 1997 Crossref, Medline, Google Scholar |
| 26. | Clifford GM , Gonçalves MA , Franceschi S : Human papillomavirus types among women infected with HIV: A meta-analysis . AIDS 20 : 2337 - 2344 , 2006 Crossref, Medline, Google Scholar |
| 27. | McDonald AC , Tergas AI , Kuhn L , et al : Distribution of human papillomavirus genotypes among HIV-positive and HIV-negative women in Cape Town, South Africa . Front Oncol 4 : 48 , 2014 Medline, Google Scholar |
| 28. | Paz-Zulueta M , Álvarez-Paredes L , Rodríguez Díaz JC , et al : Prevalence of high-risk HPV genotypes, categorised by their quadrivalent and nine-valent HPV vaccination coverage, and the genotype association with high-grade lesions . BMC Cancer 18 : 112 , 2018 Medline, Google Scholar |
| 29. | Wheeler CM , Kjaer SK , Sigurdsson K , et al : The impact of quadrivalent human papillomavirus (HPV; types 6, 11, 16, and 18) L1 virus-like particle vaccine on infection and disease due to oncogenic nonvaccine HPV types in sexually active women aged 16-26 years . J Infect Dis 199 : 936 - 944 , 2009 Medline, Google Scholar |
| 30. | Harper DM , Franco EL , Wheeler CM , et al : Sustained efficacy up to 4.5 years of a bivalent L1 virus-like particle vaccine against human papillomavirus types 16 and 18: Follow-up from a randomised control trial . Lancet 367 : 1247 - 1255 , 2006 Crossref, Medline, Google Scholar |
| 31. | Kavanagh K , Pollock KG , Cuschieri K , et al : Changes in the prevalence of human papillomavirus following a national bivalent human papillomavirus vaccination programme in Scotland: A 7-year cross-sectional study . Lancet Infect Dis 17 : 1293 - 1302 , 2017 Crossref, Medline, Google Scholar |
| 32. | Brown DR , Kjaer SK , Sigurdsson K , et al : The impact of quadrivalent human papillomavirus (HPV; types 6, 11, 16, and 18) L1 virus-like particle vaccine on infection and disease due to oncogenic nonvaccine HPV types in generally HPV-naive women aged 16-26 years . J Infect Dis 199 : 926 - 935 , 2009 Crossref, Medline, Google Scholar |
| 33. | Herrero R : Human papillomavirus (HPV) vaccines: Limited cross-protection against additional HPV types . J Infect Dis 199 : 919 - 922 , 2009 Medline, Google Scholar |
| 34. | De Vincenzo R , Ricci C , Conte C , et al : HPV vaccine cross-protection: Highlights on additional clinical benefit . Gynecol Oncol 130 : 642 - 651 , 2013 Medline, Google Scholar |
| 35. | Lu B , Kumar A , Castellsagué X , et al : Efficacy and safety of prophylactic vaccines against cervical HPV infection and diseases among women: A systematic review & meta-analysis . BMC Infect Dis 11 : 13 , 2011 Crossref, Medline, Google Scholar |




