Prevalence and Genotype Distribution of High-Risk Human Papillomavirus Infection in Women with Abnormal Cervical Cytology: A Population-Based Study in Shanxi Province, China.

PURPOSE: High-risk human papillomavirus (HR-HPV) infection is widely known as the major cause of cervical cancer and there are notable differences in HR-HPV prevalence and genotype distribution in different populations. Women with abnormal cervical cytology are at increased risk of cervical cancer; however, the genotype distribution of HR-HPV in women with abnormal cervical cytology remains unclear.
METHODS: A total of 2,300 women with abnormal cervical cytology (from 39,988 women completing a baseline survey in a cohort established during June 2014 to December 2014) were enrolled in this study. All participants gave informed consent and completed a questionnaire about characteristics related to HPV infection. HPV genotypes were identified using flow-through hybridization, and cytology was assessed by the ThinPrep cytological test. Data were analyzed using SPSS 22.0 for Windows.
RESULTS: The overall prevalence of HR-HPV in the 2,300 women with abnormal cervical cytology was 32%, with single and multiple HR-HPV infections making up 70.2% and 29.8%, respectively. The top-five HR-HPV genotypes were HPV16 (13.5%), HPV58 (5.7%), HPV52 (4.9%), HPV53 (2.5%), and HPV51 (2.3%). The prevalence of HR-HPV in atypical squamous cells of undetermined significance, low-grade squamous intraepithelial lesions, and high-grade squamous intraepithelial lesions or higher was 30.8%, 36.5%, and 54.9%, respectively, showing an increasing trend with severity of cervical cytology (χ 2 trend=13.952, p<0.001). The prevalence of HPV16 and HPV33 increased significantly with the degree of cytological abnormality. HR-HPV infection risk was statistically higher in women aged 35-45 years, with low education, infrequent bathing, multiple gravidity, multiple parity, history of gynecological diseases, and premenopause.
CONCLUSION: HR-HPV infection in women with abnormal cervical cytology was 32%, and the top-five HR-HPV genotypes were HPV16, HPV58, HPV52, HPV53, and HPV51 in Shanxi Province, China. These results shed light on demographic and behavioral characteristics related to HR-HPV infection in women with abnormal cervical cytology and provide an insight for the development of HPV vaccines.

Introduction

Human papillomavirus (HPV) infection has been identified as a definite human carcinogen,1 and accounts for >50% of infection-linked cancers in females.2 Cervical cancer ranks as the fourth-commonest cancer among women worldwide, with an estimated 570,000 new cases and 311,000 deaths in 2018.3 Annually, there are 102,000 new cases and 31,000 deaths in China, and the incidence of cervical cancer shows a substantially increasing trend.4 Shanxi Province is a high-incidence area of cervical cancer in China.5 High-risk HPV (HR-HPV) has been identified as an etiological factor in cervical cancer and precancerous lesions.6

To date, >100 HPV genotypes have been identified,7 of which approximately 40 are associated with genital tract infection.8 Among the numerous members of the HPV family, 15 HR-HPV genotypes (types 16, 18, 31, 33, 35, 39, 45, 51, 52, 53, 56, 58, 59, 66, and 68) are etiologically associated with >98% of cervical cancers. Effective prophylactic vaccines against the most important carcinogenic HPV types are available, but they do not cover the main types of HR-HPV in different countries and regions.9 Moreover, prevalence and genotype distribution of HPV vary substantially by ethnicity, demographic and behavioral characteristics, and health status.10–12 However, few studies have been reported on women with abnormal cervical cytology.

Cervical cytology testing is used as a diagnostic aid for earlier detection of cervical cancer.13 Studies have found that incidence of cervical intraepithelial neoplasia (CIN) 2–3 and cervical cancer in women with smears revealing ASC-H (atypical squamous cells — cannot exclude high-grade squamous intraepithelial lesions) was 56.5% and 9.3%, respectively,14 and the incidence of CIN 2–3 and cervical cancer in women with high-grade squamous intraepithelial lesion (HSIL) smears was 62.9% and 25.8%.15 HPV-positive women with atypical squamous cells of undetermined significance (ASC-US)/low-grade squamous intraepithelial lesion (LSIL) cytology had a 5–year CIN3+ risk of 22.2%, while for those with normal cytology, this risk was 7.9%.16 These studies clearly demonstrated that women with abnormal cervical cytology had high rates of significant cervical lesions, particularly invasive cervical disease. Many studies have suggested that HR-HPV infection has a high probability of transferring persistent infection and subsequently leading to CIN and cervical cancer17 and that abnormal cervical cytology represents an early lesion of CIN and cervical cancer. It is obviously important to understand the major distribution characteristics and genotypes of HR-HPV in women with abnormal cervical cytology, in order to implement targeted measures to control the persistence of HR-HPV infection and reduce the risk of malignant progression in high-risk women.

The prevalence of HPV infection in women with abnormal cervical cytology is 47.0%–72.8%.18,19 HPV prevalence increases with increasing severity of cervical lesions: from 12% in normal cytology to 89% in invasive cervical cancer.20 The prevalence of HR-HPV in women with abnormal cervical cytology is 92.4%, significantly higher than those with normal cytology (71%).21 In Shanxi Province in China, the prevalence of HR-HPV in normal pathology, CIN1, and CIN2+ is 11.1%, 72.5%, and 91.4%, respectively.22 However, the distribution of HR-HPV infection varies in different countries and regions, and demographics as well as behavioral characteristics related to HR-HPV infection in women with cervical abnormalities are poorly understood. The primary aim of this study was to describe the prevalence and genotype distribution of HR-HPV infection in women with abnormal cervical cytology in Shanxi to provide a foundation for preventing and controlling cervical lesions.

Methods

Study Population

This study was based on data obtained from a baseline survey of a cohort established during June 2014 to December 2014 in Shanxi. The present survey was conducted in two counties (Jiexiu and Yangqv) of Shanxi. A woman was considered eligible to enter the study if she had resided in Shanxi for at least 1 year, had current or past sexual activity, was not pregnant, and had no history of cervical cancer, precancerous lesions, cervical treatments (eg, loop electrosurgical excision procedure, conization, and adnexectomy), or injected HPV vaccine and agreed to participate in the present study. A total of 39,988 women aged 19–65 years were included. All participants completed a demographic and behavioral characteristic–related questionnaire, medical examination, and ThinPrep cytological test (TCT). Of these, 37,219 were negative for intraepithelial lesion or malignancy and 2,769 showed abnormal cervical cytology: 2,305 with ASC-US, 82 with ASC-H, 316 with LSIL, 54 with HSIL, ten with squamous cervical carcinoma, and two with atypical glandular cells. Following that, 2,769 participants with abnormal cervical cytology were referred for HPV testing. Of those, 68 refused, 387 received inadequate investigation, eight had not fully completed medical examinations, and six specimens were not suitable for detection. Finally, 2,300 women were included in the present study (Figure 1). Informed consent was signed by all participants, and the study was approved by the Institutional Review Board of Shanxi Medical University (2013–003).

Figure 1

Flowchart of participants in the study. ASC-US+, atypical squamous cells of undetermined significance and above.

Data and Sample Collection

Details of demographic characteristics, history of gynecological diseases (including vaginitis, pelvic inflammation disease, uterine fibroids, hyperplasia endometria, and ovarian tumors), lifestyle habits, and sexual and reproductive history were collected by trained interviewers using a structured questionnaire. A sterile cotton swab was used to wipe a secretion specimen from the cervix, then a TCT-specimen brush was inserted into the cervical canal to a depth of approximately 1 cm to the junction of the squamous epithelium and columnar epithelium at the outer cervix and rotated clockwise for three to five circles. Next, the TCT-specimen brush was put into cell-preservation solution. The specimens were then subjected to TCT by an experienced cytologist. Afterward, cervical specimens were collected from the cervix using an HPV-testing brush and plugged into a vial containing transport medium. Eventually, samples were processed for HPV genotyping within 2 days.

ThinPrep Cytological Test

Exfoliated cervical cells were obtained from each participant using a TCT-specimen brush. All liquid-based cytology specimens were independently diagnosed by two experienced cytopathologists. If the diagnosis differed between the two cytopathologists, the sample was reviewed by a third cytopathologist and consensus obtained. Cytology results were diagnosed according to the 2001 Bethesda system and classified as negative for intraepithelial lesion or malignancy, ASC-US, ASC-H, LSIL, HSIL, squamous cervical carcinoma, and atypical glandular cells. A diagnosis of ASC-US or higher was considered cytological abnormality.

HR-HPV Detection

HPV detection and genotyping were tested using a HybriMax HPV GenoArray kit (Hybribio Biotechnology, Chaozhou, China) according to the manufacturer’s instructions. In sum, 21 types of HPV — 15 HR-HPV types (16, 18, 31, 33, 35, 39, 45, 51, 52, 53, 56, 58, 59, 66 and 68) and six low-risk HPV types (6, 11, 42, 43, 44, and CP8304) — were identified by flow-through hybridization using a TC96/G/H6 HPV DNA–amplification analyzer and an HMM2 fast nucleic acid molecule–hybridization instrument (HybriBio). Briefly, PCR was performed in a 25 μL reaction mixture containing 5 μL extracted DNA, 0.75 μL DNA Taq polymerase and 19.25 μL PCR-mix solution containing MY09/11 primer system. The PCR protocol was denaturion at 95°C for 9 minutes, followed by 40 cycles of 20 seconds at 95°C, 30 seconds at 55°C, and 30 seconds at 72°C, and finally an extension at 72°C for 5 minutes. Positive and negative controls were run in each PCR process. HPV-genotype results were determined by the position of the HPV-genotype probes on the microarray chip. With addition of NBT/BCIP solution to display the results, a positive result was indicated by a clearly visible indigo dot. Multiple dots indicated multiple infections. Participants with any of the HR-HPV–positive genotypes were regarded as HR-HPV infected. Single HR-HPV infection was defined as only one genotype of HR-HPV in a sample being detected, while multiple HR-HPV infections were two or more HR-HPV types being detected.

Statistical Analyses

Statistical analyses were performed using SPSS 22.0 (IBM, Armonk, NY, USA) for Windows. Numerical data were examined by χ2 tests, and ORs and 95% CIs were calculated using unconditional logistic regression. Graphics were constructed using GraphPad Prism 7.0 (GraphPad Software). All p-values were two-sided, and statistical significance was defined as p<0.05.

Results

Demographic Characteristics of Participants

The demographic characteristics of the study population are summarized in Table 1. Ultimately, a total of 2,300 females aged 19–65 (50±13) years were evaluated. Most patients (63.6%) had received education to junior high school and below, and only 838 (36.4%) to senior high school and above. Nearly 1,525 (66.3%) were housewives or farmers. Only 96 (4.2%) respondents were divorced, widowed, or separated. In this survey, 1877(81.6%) participants’ yearly income were <¥10,000.

Table 1

Demographic Characteristics of the 2,300 Participants

	n	%
Age (years)
<35	183	8.0
35–44	497	21.6
45–54	867	37.7
≥55	753	32.7
Education
Junior high school and below	1,462	63.6
Senior high school and above	838	36.4
Marital status
Married/de facto	2,204	95.8
Divorced/separated/widowed	96	4.2
Occupation
Farmer/housewife	1,525	66.3
Others	775	33.7
Average yearly income (¥)
<3,000	525	22.8
3,000–6,999	789	34.3
7,000–9,999	563	24.5
≥10,000	423	18.4

Distribution of HR-HPV Genotypes

A total of 736 of the 2,300 (32%) women had HR-HPV infection, with the top-five types being HPV16 (13.5%), HPV58 (5.7%), HPV52 (4.9%), HPV53 (2.5%), and HPV51 (2.3%). Prevalence of single and multiple-type HR-HPV infection was 22.5% (517 of 2,300) and 9.5% (219 of 2,300), with corresponding proportions of 70.2% (517 of 736) and 29.8% (219 of 736), respectively, in infected women. Among those with multiple infections, infection of two genotypes was most commonly observed — 76.7% (168 of 219). In addition, we found that HPV16, HPV58, and HPV39 occurred mainly as a single infection in all populations. The details are shown in Figure 2.

Figure 2

Genotype distribution of HR-HPV infection in women with abnormal cervical cytology.

Demographic and Behavioral Characteristics Related to HR-HPV Infection

Age-specific prevalence of HR-HPV infection exhibited a peak of 46.3% (230 of 497) in the age-group 35–45 years and was lowest (23.9%, 180 of 753) in those aged 55–65 years. Age-group 35–45 years, low education, less frequent bathing, multiple gravidity (number of pregnancies), multiple parity (number of deliveries), history of gynecological diseases, and premenopause placed women at high risk of HR-HPV infection (Table 2), mainly for HPV16, HPV58, and HPV52. Of note, we found high infection rates for HPV16 and HPV51 in those <35 years old. No significant association was found between occupation, marital status, average annual income, or condom use and HR-HPV infection.

Table 2

Associations Between Demographic/Behavioral Characteristics and HR-HPV Infection in Women with Abnormal Cervical Cytology

	n	HR-HPV Infection, n (%)	Wald χ2	p	OR (95% CI)
Age (years)
<35	183	65 (35.5)			1.000
35–44	497	230 (46.3)	6.257	0.012	1.564 (1.102–2.220)
45–54	867	261 (30.1)	2.064	0.151	0.782 (0.559–1.094)
≥55	753	180 (23.9)	10.123	0.001	0.570 (0.403–0.806)
Education
Junior high school and below	1,462	492 (33.7)			1.000
Senior high school and above	838	244 (29.1)	5.028	0.025	0.810 (0.674–0.974)
Occupation
Farmer/housewife	1,525	486 (31.9)			1.000
Others	775	250 (32.3)	0.036	0.850	1.018 (0.846–1.225)
Marital status
Married/de facto	2,204	699 (31.7)			1.000
Divorced/separated/widowed	96	37 (38.5)	1.957	0.162	1.350 (0.887–2.056)
Average yearly income (¥)
<3,000	525	182 (34.7)			1.000
3,000–6,999	789	251 (31.8)	1.162	0.281	0.879 (0.696–1.111)
7,000–9,999	563	175 (31.1)	1.581	0.209	0.850 (0.660–1.095)
≥10,000	423	128 (30.3)	2.056	0.151	0.818 (0.621–1.076)
Frequency of bathing
More than once a week	1,570	421 (26.8)			1.000
Once a –week to once a month	135	58 (43.0)	52.794	<0.001	2.075 (1.704–2.527)
Less than once a month	595	257 (43.2)	15.515	<0.001	2.056 (1.436–2.942)
Gravidity
<2	675	197 (29.2)			1.000
2–3	1,390	449 (32.3)	2.052	0.152	1.158 (0.947–1.415)
>3	235	90 (38.3)	6.660	0.010	1.506 (1.104–2.055)
Parity
<2	602	171 (28.4)			1.000
2	1,059	336 (31.7)	1.996	0.158	1.171 (0.941–1.459)
≥3	639	229 (35.8)	7.812	0.005	1.408 (1.108–1.789)
Menopause
No	1,124	472 (42.0)			1.000
Yes	1,176	264 (22.4)	98.422	<0.001	0.400 (0.334–0.479)
Condom use
No	2,178	703 (32.3)			1.000
Yes	122	33 (27.0)	1.445	0.229	0.778 (0.517–1.172)
History of gynecological diseases
No	2,021	582 (28.8)			1.000
Yes	279	154 (55.2)	73.386	<0.001	3.046 (2.361–3.930)

Distribution of HR-HPV Genotypes in Various Abnormal Cervical Cytology Groups

Among the 2,300 women with abnormal cervical cytology, ASC-US, LSIL, and HSIL+ accounted for 86.5%, 11.3%, and 2.2%, respectively. The prevalence of HR-HPV in ASC-US, LSIL, and HSIL+ was 30.8% (613 of 1,989), 36.5% (95 of 260), and 54.9% (28 of 51), respectively, showing an increasing trend with severity of cervical cytology (χ2trend=13.952, p<0.001), particularly HPV16 and HPV33. The top-five HR-HPV types were HPV16 (13.1%), HPV58 (5.5%), HPV52 (4.9%), HPV53 (2.6%), and HPV51 (2.4%) in ASC-US, HPV16 (12.3%), HPV58 (8.1%), HPV52 (5.8%), HPV33 (3.1%), and HPV18 (2.7%) in LSIL, and HPV16 (33.3%), HPV33 (7.8%), HPV31 (3.9%), HPV56 (3.9%), and HPV58 (2.0%) in HSIL+ (Table 3).

Table 3

HR-HPV Prevalence and Genotype Distribution Related to Cervical Cytology Status

Genotype	Abnormal Cervical Cytology			χ2	p	χ2trend	P
	ASC-US (n=1,989)	LSIL (n=260)	HSIL+ (n=51)
HR-HPV	613 (30.8)	95 (36.5)	28 (54.9)	16.028	<0.001	13.952	<0.001
HPV16	261 (13.1)	32 (12.3)	17 (33.3)	17.762	<0.001	6.236	0.013
HPV58	110 (5.5)	21 (8.1)	1 (2.0)	3.709	0.149	0.225	0.670
HPV52	97 (4.9)	15 (5.8)	1 (2.0)	1.044	0.597	0.032	0.910
HPV53	52 (2.6)	5 (1.9)	1 (2.0)	0.272	0.882	0.453	0.539
HPV51	48 (2.4)	6 (2.3)	0	0.506	0.820	0.669	0.514
HPV33	37 (1.9)	8 (3.1)	4 (7.8)	8.016	0.014	8.090	0.008
HPV31	37 (1.9)	6 (2.3)	2 (3.9)	1.952	0.363	1.090	0.369
HPV39	38 (1.9)	4 (1.5)	1 (2.0)	0.266	0.929	0.079	0.859
HPV56	37 (1.9)	5 (1.9)	2 (3.9)	1.665	0.389	0.563	0.464
HPV66	37 (1.9)	4 (1.5)	1 (2.0)	0.268	0.854	0.051	0.860
HPV68	34 (1.7)	4 (1.5)	0	0.116	1	0.593	0.563
HPV18	28 (1.4)	7 (2.7)	1 (2.0)	2.980	0.173	1.772	0.223
HPV59	18 (0.9)	1 (0.4)	0	1.800	0.380	0.006	1
HPV35	16 (0.8)	0	0	1.592	0.472	2.254	0.149
HPV45	3 (0.2)	0	0	1.136	1	0.420	1

Note: HSIL+, diagnosis of HSIL and/or higher.

Discussion

In the present study, HR-HPV prevalence was 32% in women with abnormal cervical cytology, higher than population-based prevalence observed in Shanxi (15.2%).22 Castellsagué et al23 reported that HPV16, HPV18, HPV31, HPV58, and HPV52 were the top-five HR-HPV types in both normal and abnormal cervical cytology, according to estimates from 193 countries. Hernan et al24 found that HPV16, HPV53, HPV52, HPV58, and HPV59 were the top-five genotypes in women with abnormal cervical cytology from Bogotá, Colombia. Our findings indicated that the top-five HR-HPV genotypes were HPV16, HPV58, HPV52, HPV53, and HPV51 in women with abnormal cervical cytology in Shanxi, while HPV18 was in 12th position. The distribution of HR-HPV genotypes varies greatly worldwide, and these differences might be related to the complex geographical and biological interplay among different HPV genotypes and host immunogenetic factors.

Our results showed that the proportions of single and multiple HR-HPV infections were 70.2% and 29.8%, respectively. HPV16, HPV58, and HPV39 occur mainly as a single infection in all populations. A study conducted in Cyprus also showed that single HR-HPV infection accounted for 77.2% and multiple HR-HPV infections 22.8% in women with both abnormal cervical cytology and HR-HPV.19 There are currently three prophylactic vaccines targeting various HPV types: the bivalent vaccine (HPV16/18), the quadrivalent vaccine (HPV16/18/6/11) and the nine-valent vaccine (HPV16/18/6/11/31/33/45/52/58),25 but they do not cover the most frequent HR-HPV genotypes in Shanxi, such as HPV53 and HPV51, from our present study. Therefore, further work will be needed to develop next-generation vaccines that can target more HR-HPV types.

Li et al26 reported that HR-HPV prevalence presented two age peaks in the general Chinese population: one at age 15–24 years and the other in women aged 35–49 years. HR-HPV infection is a sexually transmitted disease, and the higher rate of infection in these women might point to more frequent sexual intercourse. As such, cervical cells might be more vulnerable to damage and HR-HPV infection. In our study, the prevalence of HR-HPV was much higher in those aged 35–45 years abnormal cervical cytology. These results were inconsistent with those for all age-groups, suggesting that women aged 35–45 years should be regarded as a high-risk group and a focus of prevention among women with abnormal cervical cytology. Additionally, we observed that HPV16, HPV58, and HPV52 presented the highest prevalence in women aged 35–45 years, suggesting that HPV16, HPV58, and HPV52 infection in this region should be given attention. Moreover, the nine-valent vaccine appears to be more valuable than the quadrivalent vaccine in women with abnormal cervical cytology.

Remarkably, it should be noted that the HR-HPV infection rate decreased with increasing education, the same as in a study conducted in Yunnan Province, China.27 Women with junior high school education and below had limited knowledge about cervical cancer and its prevention. It will be necessary to improve the education of women to prevent HR-HPV infection and cervical lesions. Our findings provide evidence that women bathing less frequently were more vulnerable to HR-HPV infection. Bathing infrequently results in increasing viral or bacterial reproduction, which can be curbed by frequently cleaning the lower genital tract.28

We found that the prevalence of HR-HPV significantly increased with increasing gravidity and parity, which might because of the hormonal changes of pregnancy. The female sex hormones — estrogen and progesterone — might influence susceptibility to HR-HPV infection. Changes in susceptibility to HR-HPV might also be the result of hormone-induced changes in host adaptive immunoresponse.29 Our previous studies have shown that menopause was a protective factor in HR-HPV infection, which could be ascribed to the decreased estrogen level of postmenopausal women.30–32 These results indicated that HPV DNA testing and effective HPV vaccines should be implemented in premenopausal women. Studies have reported that women with gynecological diseases were prone to HR-HPV infection, due to unstable sexual hormone levels, poor function of the immune system, lower resistance to HR-HPV, and weakviral clearance.33,34

Our findings indicated that HPV16 is the most frequently identified type in women with ASC-US, LSIL, and HSIL+ in Shanxi, confirming that HPV16 is the main HR-HPV type associated with cervical lesions. Furthermore, we found that infection rates of HPV16 and HPV33 increased significantly with degree of cytological abnormality. Our analyses suggest that HPV16 and HPV33 infection may be biomarkers to predict cervical lesion progression in women with abnormal cervical cytology.

Conclusion

HR-HPV prevalence was 32% and the top-five HR-HPV genotypes were HPV16, HPV58, HPV52, HPV53, and HPV51 in women with abnormal cervical cytology in Shanxi. Age 35–45 years, low education, less frequent bathing, multiple gravidity, multiple parity, history of gynecological diseases, and premenopause placed women at high risk of HR-HPV infection. Our results provide insight into HR-HPV infection–related demographic and behavioral characteristics in women with abnormal cervical cytology, and provide evidence for establishing strategies to control HR-HPV infection and reducing the risk of malignant progression in high-risk women. Genotype distribution in the region may point to a promising future for HPV-vaccine research and development. Based on these important leads, prospective cohort studies are needed to provide more powerful evidence.