Population-based study on the prevalence of and risk factors for human papillomavirus infection in Qujing of Yunnan province, Southwest China
- Sun Lu-lu1,
- Jin Qiong1,
- Li Hui2,
- Zhou Xian-rong3,
- Song Zhi-qin4,
- Cheng Xue-mei1,
- Tao Tao1,
- Liang Bing1,
- Xu Lin5,
- Wang Yu-rong6,
- Zhen Yan7,
- He Ji-wen8 and
- Shen Keng1Email author
© Lu-lu et al.; licensee BioMed Central Ltd. 2012
Received: 21 April 2011
Accepted: 2 August 2012
Published: 8 August 2012
Human papillomavirus (HPV) infection causes cervical cancer and premalignant lesions of the cervix. Prevalence of HPV infection and HPV genotypes vary among different regions. However there is no data on the prevalence of HPV infection and HPV genotypes from southwest China. This study was undertaken to determine the prevalence of and risk factors for HR-HPV infection in Qujing of Yunnan province, southwest China to provide comprehensive baseline data for future screening strategies.
A sample of 5936 women was chosen by the multi-stage stratified cluster sampling method with selection probabilities proportional to size (PPS). An epidemiological questionnaire was conducted via a face-to-face interview and cervical specimens were taken for HPV DNA testing by Digene Hybrid Capture 2 (HC2) test. HPV Genotyping Reverse Hybridization Test was used for HPV genotyping. Proportions were compared by Chi-squared tests, and logistic regression was utilized to evaluate risk factors.
The median age was 38 years and the inter-quartile range was from 31 years to 47 years. 97.3% of the study population was Han nationality. Overall prevalence of HR-HPV infection was 8.3% (494/5936) and bimodal age distribution of HPV infection was observed. The five most prevalent HR-HPV genotypes were HPV-16(3.4%), HPV-56(1.7%), HPV-58(1.4%), HPV-33(1.2%) and HPV-52(0.88%). Multiple HPV infections were identified in 50.5% (208/412) of the positive genotyping specimens. Multivariate logistic regression model indicated that parity (OR = 1.35, 95% CI: 1.18-1.53, p < 0.0001) was a risk factor for HR-HPV infection, and age of 50–65 years (OR = 0.60, 95% CI: 0.45-0.80, p = 0.0005), being married or in stable relationship (OR = 0.55, 95% CI: 0.31-0.96, p = 0.035) were protective factors.
This study provided baseline data on HR-HPV prevalence in the general female population in Qujing of Yunnan province, southwest China. The finding of multiple HPV infections and bimodal age distribution revealed that HPV screening is necessary for perimenopausal women in future.
KeywordsHuman papillomavirus Genotype Cervix Epidemiology China
Polymerase chain reaction
Cervical cancer is the second most common cancer in women worldwide.  More than 100 HPV types have been identified and nearly 40 types are known to infect the genital tract . HPV 16 and 18 have been shown to be causative in the etiology of cervical cancer . Women with persistent infection with oncogenic HPV subtypes may develop cervical cancer.
There are about 20 high-risk subtypes of HPV associated with cervical cancer and precancerous lesions.  Overall, HPV-16 and HPV-18 account for approximately 70% of all cervical cancer diagnosed worldwide each year . HPV genotypes differ greatly in their geographic distribution [5, 6]. However, there is no data on the prevalence of HPV infection from southwest China. Qujing located at an altitude of 2000 meters in southwest China. As the second largest city of Yunnan province, the population is 6.16 million, and there are various ethnic groups. This study seeks to describe the prevalence of HPV infection, HPV genotype distribution as well as risk factors for HR-HPV infection in the general female population of Qujing, Yunnan province, southwest China to provide comprehensive baseline data for assessing its potential impact on cervical pathology and future screening strategies of cervical cancer in China.
Study population and survey sampling design
n0: estimated sample size;
π: expected prevalence of HPV infection, 5% was used;
δ: allowing within 15% of the range for a prevalence of 5%, δ = 0.15π;
α = 0.05, u α = 1.96.
In the expectation of a response rate of 80% and expanded sample size of 150%, the sample size was adjusted using n = 1.5n0: 0.80 = 6083. The estimated sample size for the survey was 6100 rounded to the nearest hundred. Non-contactables were defined as individuals who were not contactable on three occasions, and refusals defined as individuals who declined to participate in the study. And women who had a hysterectomy, were menstruating or pregnant at the time of the study were excluded. Sample size supplement is not necessary if non-contactables and refusals account for less than 20%.
Multi-stage stratified cluster sampling method
Central city zone
Community units / Township
Community Residents’ Committee
Verbal informed consent was obtained before the questionnaire interview. Participants were asked to complete a questionnaire via a face-to-face interview, which was designed to elicit socio-demographic status, sexual behavior, history of sexually transmitted infection (STI), parity, and other probable risk factors for cervical cancer. The questionnaire was piloted on 85 outpatients from gynecological clinic in Peking Union Medical College Hospital (PUMCH). The internal consistency reliability (Cronbach’s α) for the full scale was 0.825. Furthermore, face validity and content validity of the questionnaire were both confirmed by epidemiologists and experts from Department of Obstetrics and Gynecology, PUMCH. The results indicated that the questionnaire had good reliability and validity.
Cervical specimen and HR-HPV detection
Informed consent was obtained before the collection of the cervical specimens. Samples of exfoliated cervical cells were collected with DNAPap cervical sampler brush (Digene) during gynecological examinations. The sampler was inserted 1–1.5 cm into the endocervical canal and rotated 3–5 full turns in counterclockwise direction. The cervical sampler was then placed into transport medium and stored at 4°C. HR-HPV detection was carried out using the Digene Hybrid Capture 2 ® (HC2) test according to the manufacture’s instructions, as described previously. The results are given as a relative light unit (RLU) ratio and a positive result was defined as RLU/Co ≥ 1.00 (5000 copies of HPV genome).
Genotyping was carried out by HPV Genotyping Reverse Hybridization Test RUO detection kit. Phenol-chloroform was used to extract HPV DNA from the residual of HC2 positive samples. And the human β-globin gene was used as internal control (IC) to verify the quality for polymerase chain reaction (PCR). DNA of specimens with negative internal control amplification result were re-extracted and re-amplified. GP5+/6+ consensus primers were used according to the manufacture’s instructions .
HPV Genotyping Reverse Hybridization Test identifies 18 high risk or probably high risk HPV genotypes. Besides 13 HR-HPV genotypes covered by HC2 (including HPV-16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, 59 and 68), there are 5 other HPV genotypes that are not included in the current HC2 probe cocktail (including HPV-26, 53, 66, 73 and 82) .
The final results were detected by colorimetric change on the membrane strip immortalized with type-specific oligonucleotides probes. All procedures were carried out according to the manufacture's instructions. HR-HPV positive cases were defined as those with single-type (a single HR type) infection and those with multiple-type (two or more HR types) infection. Type-specific prevalence included that in both single-type and multiple-type infection .
Quality assurance and quality control
Several efforts were made to standardize all procedures in each study site. Manuals were written with instructions on how to conduct the fieldwork, interviewing, collecting cervical samples, processing and storing the samples collected. The training program for doctors and crew of local stuff was held a week before investigation. Data were key-entered twice by two experienced typists. Independent databases were formed and were compared with each other. Any mismatch was corrected after inter-comparison of the databases and a proportion of 10% of the records in the final database were drawn up randomly for re-correction manually.
For the laboratory tests, HPV DNA testing by the HC2 was conducted using the same commercial kit and the same machine in local lab of Qujing. Cervical samples were transported to central lab in Shenzhen for HPV genotyping thereafter. All of the laboratory technicians were blinded to the clinical profile of the clients. HPV genotyping were independent of HPV DNA tests.
Ethical approval was granted by the Research Ethics Committee of Peking Union Medical College Hospital, Peking Union Medical College & the Chinese Academy of Medical Sciences.
Descriptive statistics (medians and proportions) were utilized to characterize the variables. Proportions were compared using Chi-squared tests. Odds ratios (ORs) with 95% confidence intervals (CIs) were calculated using unconditional logistic regression (Mantel-Haenszel analysis) to estimate independent risk factors for HR-HPV infection. Age and variables associated with HR-HPV infection on univariate logistic regression model were included as candidates in multivariate model. Parity, age of 50–65 years and marital status were included in the final model by stepwise selection.
Data analysis was conducted using the SAS System for Windows (SAS Institute Inc., USA). A level of 0.05 was chosen to indicate statistical significance.
Qujing is the second largest city of Yunnan province in southwestern China, with an average altitude of 2000 m. It has a population of 6 million and has eight ethnic groups. Because of the limited medical resources, cervical cancer screening is not widely utilized, especially in rural areas of Qujing. Of the 6116 women who participated in the study, 180 were excluded because of inadequate cytology specimen, negative β-globin or inadequate questionnaire information, leaving a total of 5936 women included in the statistical analyses. The median age was 38 years and the inter-quartile range was from 31 years to 47 years. 97.3% of the study population was Han nationality.
HPV prevalence for specific types, both overall and by age group (n = 5936) *
Prevalence by age group,%
18-24 (n = 374)
25-29 (n = 856)
30-34 (n = 977)
35-39 (n = 1035)
40-44 (n = 813)
45-49 (n = 781)
50-54 (n = 447)
55-59 (n = 383)
60-65 (n = 270)
Any HR types
Risk profile for HR-HPV infection
Major risk factors for HR-HPV infection adjusted for age
HPV DNA positive No. (%)
The Han nationality
Income per year
Married or in stable relationship
Age at first intercourse
Age of menarche
Oral contraceptive pills
Life time sexual partners
Husband’s life time sexual partners
Summary of the final multivariate model (Stepwise Selection, slentry = 0.5, slstay = 0.05)
1.182 ~ 1.534
50-65 years old
0.446 ~ 0.797
Being married or in stable relationship
0.312 ~ 0.958
Differences in HR-HPV prevalence are seen between countries. Previous studies have documented that HR-HPV prevalence in women attending routine cervical screening program in Europe was 15.6%. It was 13-18% in Latin America, 17% in India and 26% in Nigeria , which were much higher than the prevalence of this study in China. To date, there has been no population-based data in the prevalence of HR-HPV infection in southwestern China. The present study demonstrated that HR-HPV infection in the general female population in Qujing is 8.3%. Compared with other studies, HR-HPV infection in our study population is lower than that in central China of Shanxi (12.2%), northeast China of Shenyang (11.7%), southeast China of Zhejiang (10.2%), and southeastern coast China of Shenzhen (13.5%) and Taiwan (32.4%) [12–17]. However some studies that reported even lower HR-HPV infection are the ones in Beijing (5.8%) and our previous study in Tibet (7.1%) [7, 8]. HR-HPV infection rate may vary among different regions because of different culture and different life styles. One interpretation attributed to the relatively lower prevalence of HPV infection in this study may be that Qujing located in a close mountainous area keeping with the traditional Chinese social habits and life styles. Stable marital status and sexual partner may provide a protection for HPV infection [18–22].
The persistence and clearance of HPV infection result in the actual HPV infection rate. [23, 24] Bimodal age distribution of HR-HPV is also observed in this study. The first peak of HR-HPV infection is obvious in women aged 18–24 years, which may be due to primary exposure to HPV after sexual initiation and lack of adaptive immune responses in young women. [25–27] The second peak of HR-HPV infection is observed in women aged 55–59 years which should be deserved more attention. It is assumed that the immunologic and physiologic deregulation caused by hormone fluctuations may explain the high infection rate around menopausal women [28, 29]. Nevertheless, the positive HR-HPV infection around menopausal period women may imply the persistent HR-HPV infection [28, 30, 31]. Therefore, HPV detection is clinically valuable for perimenopausal women in cervical cancer screening program. For perimenopausal women with positive HR-HPV DNA test result, regular follow-up and standard management should be considered, including repeated cytology and colposcopy if ASCUS or more severe cytological findings were reported, because they have higher risk for the development of cervical cancer . However, because of regional economy imbalance in China, HPV testing is currently not available in some regions, especially in west China.
The five most common HPV types in this study were HPV-16, -56, -58, -33 and −52 respectively. Some studies had shown that HPV-52 and −58 were more predominant and overrepresented in cervical cancer cases in Asia [5, 12, 13, 15, 17, 32]. Data of this study revealed that HPV-16,-56, -58, -33 and −52 may play a role in the etiopathology of cervical cancer. Meanwhile, the high prevalence of multiple HPV infection in this study is also detected in young women with sexual activity as well as in women around menopausal period. Some researches have reported that women with multiple HPV infection were at increased risk of precancerous lesions. [33, 34] Our previous study in Tibet revealed that significant difference was evident when comparing the percentage of abnormal cytological results in multiple HPV infection with single HPV infection . Therefore, the role of multiple HPV infection in etiology of cervical cancer should be further investigated, especially for the research on HPV vaccine.
HPV vaccines have shown type-restricted prophylactic efficacy for genital lesions. However, it does not have clinical significance for women who have been infected with HPV. As a result, HPV vaccine program will be of greater benefit to younger women. We assume that the HPV prophylactic vaccines including HPV-52 and −58 may offer higher protection for young women in China and other Asian countries.
This study is the first report on HR-HPV prevalence, HPV genotypes as well as the risk factors for HR-HPV infection in the general female population in southwestern China. The baseline data of this study will contribute to making strategies for women’s health care and implementation of cervical cancer screening program in China.
We sincerely thank all the staff and participants who took part in this study. This work was funded by National Science and Technology Infrastructure Program “The National Key Technologies R&D Program of China” (grant 2008BAI57B01) and China Women Development Foundation “Healthy Mother Express” Research Fund.
- Parkin DM, Bray F, Ferlay J, Pisani P: Global cancer statistics, 2002. CA –A Cancer Journal for Clinicians 2005, 55:74–108.PubMedView Article
- Munoz N, Castellsague X, de Gonzalez AB, Gissmann L: Chapter 1: HPV in the etiology of human cancer. Vaccine 2006,24(Suppl 3):1–10.View Article
- Walboomers JM, Jacobs MV, Manos MM, Bosch FX, Kummer JA, Shah KV, Snijders PJ, Peto J, Meijer CJ, Muñoz N: Human papillomavirus is a necessary cause of invasive cervical cancer worldwide. Journal of Pathology 1999, 189:12–19.PubMedView Article
- Smith JS, Lindsay L, Hoots B, Keys J, Franceschi S, Winer R: Clifford GM Human papillomavirus type distribution in invasive cervical cancer and high-grade cervical lesions: a meta-analysis update. Int J Cancer 2007, 121:621–632.PubMedView Article
- Clifford GM, Gallus S, Herrero R, Munoz N, Snijders PJ, Vaccarella S, Anh PT, Ferreccio C, Hieu NT, Matos E, Molano M, Rajkumar R, Ronco G, de Sanjose S, Shin HR, Sukvirach S, Thomas JO, Tunsakul S, Meijer CJ, Franceschi S: Worldwide distribution of human papillomavirus types in cytologically normal women in the International Agency for Research on Cancer HPV prevalence surveys: a pooled analysis. Lancet 2005, 366:991–998.PubMedView Article
- Li Y, Wang Y, Jia C, Ma Y, Lan Y, Wang S: Detection of human papillomavirus genotypes with liquid bead microarray in cervical lesions of northern Chinese patients. Cancer Genet Cytogenet 2008, 182:12–17.PubMedView Article
- Zhao R, Zhang WY, Wu MH, Zhang SW, Pan J, Zhu L, Zhang YP, Li H, Gu YS, Liu XZ: Human papillomavirus infection in Beijing, People’s Republic of China: a population-based study. Br J Cancer 2009, 101:1635–1640.PubMedView Article
- Jin Q, Shen K, Li H, Zhou XR, Huang HF, LENG JH: Age-specific prevalence of human papillomavirus by grade of cervical cytology in Tibetan women. Chin Med J 2010, 123:2004–2011.PubMed
- Manos MM, Ting Y, Wright DK, Lewis AJ, Broker TR, Wolinsky SM: Use of polymerase chain reaction amplification for the detection of genital human papillomaviruses. Cancer Cells 1989, 7:209–214.
- Bouvard V, Baan R, Straif K, Grosse Y, Secretan B, El Ghissassi F, Benbrahim-Tallaa L, Guha N, Freeman C, Galichet L, Cogliano V: A review of human carcinogens–Part B: biological agents. Lancet Oncol 2009, 10:321–322.PubMedView Article
- de Villiers EM, Fauquet C, Broker TR, Bernard HU, zur Hausen H: Classification of papillomaviruses. Virology 2004, 324:17–27.PubMedView Article
- Dai M, Bao YP, Li N, Clifford GM, Vaccarella S, Snijders PJ, Huang RD, Sun LX, Meijer CJ, Qiao YL, Franceschi S: Human papillomavirus infection in Shanxi Province, People’s Republic of China: a population-based study. Br J Cancer 2006, 95:96–101.PubMedView Article
- Li LK, Dai M, Clifford GM, Yao WQ, Arslan A, Li N, Shi JF, Snijders PJ, Meijer CJ, Qiao YL, Franceschi S: Human papillomavirus infection in Shenyang City, People’s Republic of China: A population-based study. Br J Cancer 2006, 95:1593–1597.PubMedView Article
- Ye J, Cheng X, Chen X, Ye F, Lü W, Xie X: Prevalence and risk profile of cervical human papillomavirus infection in Zhejiang Province, southeast China: a population-based study. Virology Journal 2010, 7:66.PubMedView Article
- Wu RF, Dai M, Qiao YL, Clifford GM, Liu ZH, Arslan A, Li N, Shi JF, Snijders PJ, Meijer CJ, Franceschi S: Human papillomavirus infection in women in Shenzhen City, People’s Republic of China, a population typical of recent Chinese urbanization. Int J Cancer 2007, 121:1306–1311.PubMedView Article
- Tsao KC, Huang CG, Kuo YB, Chang TC, Sun CF, Chang CA, Yang SL, Chan EC: Prevalence of Human Papillomavirus Genotypes in Northern Taiwanese Women. J Med Virol 2010, 82:1739–1745.PubMedView Article
- Bao YP, Li N, Smith JS, Qiao YL: Human papillomavirus type distribution in women from Asia: a meta-analysis. Int J Gynecol Cancer 2008, 18:71–79.PubMedView Article
- Deacon JM, Evans CD, Yule R, Desai M, Binns W, Taylor C, Peto J: Sexual behaviour and smoking as determinants of cervical HPV infection and of CIN3 among those infected: a case–control study nested within the Manchester cohort. Br J Cancer 2000, 83:1565–1572.PubMedView Article
- Moscicki AB, Hills N, Shiboski S, Powell K, Jay N, Hanson E, Miller S, Clayton L, Farhat S, Broering J, Darragh T, Palefsky J: Risks for incident human papillomavirus infection and low-grade squamous intraepithelial lesion development in young females. Jama 2001, 285:2995–3002.PubMedView Article
- La Vecchia C, Franceschi S, Decarli A, Fasoli M, Gentile A, Parazzini F, Regallo M: Sexual factors, venereal diseases, and the risk of intraepithelial and invasive cervical neoplasia. Cancer 1986, 58:935–941.PubMedView Article
- Winer RL, Lee SK, Hughes JP, Adam DE, Kiviat NB, Koutsky LA: Genital human papillomavirus infection: incidence and risk factors in a cohort of female university students. Am J Epidemiol 2003, 157:218–226.PubMedView Article
- Chan PK, Ho WC, Wong MC, Chang AR, Chor JS, Yu MY: Epidemiologic risk profile of infection with different groups of human papillomaviruses. J Med Virol 2009, 81:1635–1644.PubMedView Article
- Brown DR, Shew ML, Qadadri B, Neptune N, Vargas M, Tu W, Juliar BE, Breen TE, Fortenberry JD: A longitudinal study of genital human papillomavirus infection in a cohort of closely followed adolescent women. J Infect Dis 2005, 191:182–192.PubMedView Article
- Goodman MT, Shvetsov YB, McDuffie K, Wilkens LR, Zhu X, Thompson PJ, Ning L, Killeen J, Kamemoto L, Hernandez BY: Prevalence, acquisition, and clearance of cervical human papillomavirus infection among women with normal cytology: Hawaii Human Papillomavirus Cohort Study. Cancer Res 2008, 68:8813–8824.PubMedView Article
- Winer RL, Feng Q, Hughes JP, O’Reilly S, Kiviat NB, Koutsky LA: Risk of female human papillomavirus acquisition associated with first male sex partner. J Infect Dis 2008, 197:279–282.PubMedView Article
- Moscicki AB, Shiboski S, Broering J, Powell K, Clayton L, Jay N, Darragh TM, Brescia R, Kanowitz S, Miller SB, Stone J, Hanson E, Palefsky J: The natural history of human papillomavirus infection as measured by repeated DNA testing in adolescent and young women. J Pediatr 1998, 132:277–284.PubMedView Article
- Sherman ME, Schiffman M, Cox JT: Effects of age and human papilloma viral load on colposcopy triage: data from the randomized Atypical Squamous Cells of Undetermined Significance/Low-Grade Squamous Intraepithelial Lesion Triage Study (ALTS). J Natl Cancer Inst 2002, 94:102–107.PubMedView Article
- Althoff KN, Paul P, Burke AE, Viscidi R, Sangaramoorthy M, Gravitt PE: Correlates of cervicovaginal human papillomavirus detection in perimenopausal women. J Womens Health (Larchmt) 2009, 18:1341–1346.View Article
- Jacobs MV, Walboomers JM, Snijders PJ, Voorhorst FJ, Verheijen RH, Fransen-Daalmeijer N, Meijer CJ: Distribution of 37 mucosotropic HPV types in women with cytologically normal cervical smears: the age-related patterns for high-risk and low-risk types. Int J Cancer 2000, 87:221–227.PubMedView Article
- Castle PE, Schiffman M, Herrero R, Hildesheim A, Rodriguez AC, Bratti MC, Sherman ME, Wacholder S, Tarone R, Burk RD: A prospective study of age trends in cervical human papillomavirus acquisition and persistence in Guanacaste, Costa Rica. J Infect Dis 2005, 191:1808–1816.PubMedView Article
- Franco EL, Villa LL, Sobrinho JP, Prado JM, Rousseau MC, Desy M, Rohan TE: Epidemiology of acquisition and clearance of cervical human papillomavirus infection in women from a high-risk area for cervical cancer. J Infect Dis 1999, 180:1415–1143.PubMedView Article
- Lin H, Ma YY, Moh JS, Ou YC, Shen SY, ChangChien CC: High prevalence of genital human papillomavirus type 52 and 58 infection in women attending gynecologic practitioners in South Taiwan. Gynecol Oncol 2006, 101:40–45.PubMedView Article
- Herrero R, Castle PE, Schiffman M, Bratti MC, Hildesheim A, Morales J, Alfaro M, Sherman ME, Wacholder S, Chen S, Rodriguez AC, Burk RD: Epidemiologic profile of type-specific human papillomavirus infection and cervical neoplasia in Guanacaste, Costa Rica. J Infect Dis 2005, 191:1796–1807.PubMedView Article
- An HJ, Cho NH, Lee SY, Kim IH, Lee C, Kim SJ, Mun MS, Kim SH, Jeong JK: Correlation of cervical carcinoma and precancerous lesions with human papillomavirus (HPV) genotypes detected with the HPV DNA chip microarray method. Cancer 2003, 97:1672–1680.PubMedView Article
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 (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.