Prevalence and Genotype Distribution of Human Papillomavirus Among Healthy Females in Beijing, China, 2016-2019.

OBJECTIVE: Human papillomavirus (HPV) infection, especially with high-risk HPV (HR-HPV) genotypes, is closely associated with cervical cancer. This study aimed to observe the epidemiological characteristics of HPV infection among healthy women in Beijing, China.
MATERIALS AND METHODS: Cervical specimens were collected from 29,436 healthy women, who underwent health check-ups in Peking Union Medical College Hospital between 2016 and 2019. A commercial kit was used for the detection of 15 HR-HPV and two low-risk HPV (LR-HPV) genotypes.
RESULTS: A total of 3586 (12.18%) participants tested positive for HPV, 3467 of which were infected with HR-HPVs. The most prevalent genotypes were HPV52, 58, 16, 51, and 56. Moreover, while infection with a single genotype (9.84%) was more prevalent, HPV16+52 was the most common combination in those infected with multiple HPVs. Furthermore, the highest infection rate among age groups was in women aged <25 years (20.92%). No significant difference in the prevalence was observed from 2016 to 2019. However, HPV incidence in Beijing was significantly different than that in all other areas in China, except for Zhengzhou (p < 0.05).
CONCLUSION: Our findings could serve as potential reference for better understanding of the epidemiological characteristics of HPV infection in Beijing.

Introduction

Cervical cancer is one of the most common malignancies in women. In 2018, it was ranked fourth in terms of incidence and mortality in females.1 Globally, roughly 570,000 new cases and 311,000 deaths worldwide were associated with this disease in 2018, one-fifth of which were reported in China.2 Moreover, the global average age of women diagnosed with cervical cancer is 53 years, ranging from 44 to 68 years, and the global average age at death from cervical cancer is 59 years, ranging from 45 to 76 years.3 Hence, cervical cancer remains a major public health problem affecting females, especially those in China.

Human papillomavirus (HPV) has been detected in almost all cases of cervical cancers.4 At present, more than 100 HPV genotypes are known, and these are classified into high-risk HPV (HR-HPV), intermediate risk HPV (IR-HPV) and low-risk HPV (LR-HPV), based on their ability to induce premalignant and malignant transformation.5,6 Over 100 cervical cancer genotypes have been recognized, and at least 15 of those have been classified as HR-HPV.7 Persistent infection with HR-HPV is the main risk factor of cervical cancer, whereas infection with LR-HPV is related with genital warts.8 Two HR-HPV genotypes, HPV16 and HPV18, account for roughly 70% of all cervical cancer cases.9 Other genotypes associated with cervical adenocarcinomas include HPV31, 33, 35, 39, 45, 51, 52, 56, 58, 59, and 68.10

HPV vaccines have been developed as a means of preventing cervical cancer and its precancerous lesions and have been effective in preventing and controlling HPV infection in some developed countries.11–13 However, because of differences in demographic, geographic, and ethnic factors, the prevalence and genotype distribution of HPV also varies across countries and regions worldwide,14–16 and this might lead to different distribution levels of cervical cancer subtypes and precancerous cervical lesions.17,18 Hence, it is very important to fully investigate the prevalence and genotype distribution of HPV in order to develop proper monitoring, and prevention and control tactics against the virus.

In this study, data from healthy females in Beijing area were analyzed, and the prevalence and distribution of HPV were evaluated in order to provide a scientific basis for the development and evaluation of the effectiveness of HPV vaccines and HPV screening in this region.

Materials and Methods

Study Population

This was a retrospective study to determine the prevalence of HPV infection. We performed a cross-sectional analysis, wherein we collected data of 29,436 healthy females aged 18–86 years old who visited the Peking Union Medical College Hospital and received a routine medical examination between 2016 and 2019. To be included in this study, the healthy women: (1) were aged 18 years or more; (2) were a permanent resident of the local area; (3) were mentally and physically healthy; (4) were not pregnant; and (5) had not undergone a total hysterectomy or cervicectomy. All enrolled patients underwent a routine gynecological check; their cervical specimens were examined for HPV DNA genotyping.

Ethics Statement

The study adhered to strict the principles of the Declaration of Helsinki and was approved by the Ethics Committee of Peking Union Medical College Hospital of the Chinese Academy of Medical Sciences (No. S-K1573). Informed consent was taken from all participants who joined the study.

Sample Collection

Sample collection was performed by a gynecologist. Cervical exfoliated cell specimens were collected using a disposable cervical brush. Thereafter, the brush was transferred to a sampling tube containing a cell preservation liquid. All the samples were stored at 4 °C and then transported to the molecular biology laboratory; subsequent tests were performed within 48 h.

HPV DNA Testing

HPV DNA was extracted using an Autrax automatic nucleic acid extraction workstation (Liferiver Biotechnology Limited Corp. Shanghai, China). HPV DNA was amplified using a Roche LightCycler 480 instrument (Roche, USA) according to the manufacturer’s instructions. HPV genotyping was performed using two HPV genotyping kits for 17 HPV types (Liferiver Biotechnology Limited Corp. Shanghai, China), wherein one kit included nine HR-HPV types (HPV16, 18, 39, 45, 56, 59, 66, 68 and 82) and six IR-HPV types (HPV31, 33, 35, 51, 52, and 58), and another included two LR-HPV types (HPV6 and 11). These assay kits rely on a TaqMan-based real-time fluorescent quantitative polymerase chain reaction (PCR) method. The amplification procedure of the HR-HPV and IR-HPV genotyping kit was performed as follows: 2 min of denaturation at 94 °C, followed by 40 cycles of denaturation at 93 °C for 10 s, 31 s of annealing and elongation at 62 °C. Single-point fluorescence was detected at 62 °C. The lowest limit of detection (LOD) was 1 × 104 copies/mL. The LR-HPV genotyping kit comprised specific primers and probes designed according to the conserved regions of the HPV6 and HPV11 genomes, and was used to detect the double positive sample of HPV6 and HPV11. The amplification procedure was as follows: 2 min of pre-denaturation at 37 °C, then 2 min of denaturation at 94 °C, followed by 40 cycles of denaturation at 93 °C for 15 s, 60 s of annealing and elongation at 60 °C. Single-point fluorescence was detected at 62 °C. The lowest LOD was 1 × 103 copies/mL. Internal quality control and external quality assessment (negative control and positive control) were performed to ascertain the reliability of the results. Internal quality control was used to ensure that the cervical epithelial cells had been collected.

Analysis of Distribution and Prevalence

Genotypic data from 2016 to 2019 were retroactively collected from our hospital database and subsequently analyzed to perform a comparison between HPV infection prevalence in Beijing and other areas of China, for which some recently published regional studies based on a larger population sample were collected.19–30 The other regional studies were categorized into four macro-geographical regions: east (Shanghai, Zhejiang, Jiangsu, and Jiangxi), west (Xinjiang and Sichuan), south (Wuhan, Guangdong, and Yunnan), and north (Zhengzhou, Shaanxi, and Shandong).

Statistical Analysis

Participant age data were presented as the mean ± standard deviation. All participants were divided into 10 groups according to their age during HPV testing, with each group covering a 5-year range. Statistical analyses were performed using SPSS version 20.0 (IBM Corp., Armonk, NY, USA). Differences between groups were compared by chi-squared test (χ2). For all statistical analyses, p-values were two-sided and p < 0.05 was considered statistically significant.

Results

Overall and Type-Specific HPV Prevalence

This study analyzed data of HPV tests of 29,436 healthy women collected between 2016 and 2019. The age of the women ranged from 18 to 86 years, and the mean age was 44.95 ± 11.38 years. Of the 29,436 subjects, 3586 women were positive for HPV, with an overall HPV infection rate of 12.18%. The constituent ratio of the age groups of HPV-positive patients is shown in Figure 1, which indicates that the groups with the highest percentages of healthy women were of age groups 45–49 years (15.28%, 548/3586) and 50–54 years (14.33%, 514/3586).

Figure 1

Constituent ratio of different age groups of HPV infection.

HR-HPV was detected in 3467 individuals, with an infection rate of 11.78% (3467/29,436). The most prevalent genotypes were HPV52, with a frequency of 3.04% (895/29,436), followed by HPV58, 16, 51 and 56, with frequencies of 1.84% (541/29,436), 1.82% (535/29,436), 1.18% (346/29,436) and 1.06% (311/29,436), respectively (Figure 2).

Figure 2

Distribution of HPV genotypes among healthy women.

The prevalence of infection with HPV16/18, HPV6/11/16/18, and HPV16/18/52/58 was 20.8% (720/3586), 23.01% (825/3586), and 52.43% (1880/3586), respectively, in the HPV-positive patients.

Distribution Characteristics of Single and Multiple Infections

The analysis of the characteristics of HPV infections across all the samples showed that infection with single and multiple HPVs were at 9.84% (2897/29436) and 2.34% (689/29436), respectively. Among those infected with multiple HPVs, double infection was the most prevalent, and affected 547 females (Table 1). Further analysis of double infections revealed ten common combinations of HPV subtypes (Table 2), of which, HPV16+52 was the most common subtype, accounting for 5.08% (35/689) of double infections.

Table 1

Frequency and Prevalence of Single and Multiple HPV Infection Among Healthy Women

Genotype	Positive Number	Prevalence (%)
Single	2897	9.84
Double	547	1.86
Triple	104	0.35
Quadruple	30	0.10
Quintet	7	0.02
Sextuple	1	0.01

Table 2

The Top Ten Combinations of Genotypes in HPV Double Infections

Combination of Types	Positive Number	Prevalence (%)
16+52	35	5.08
39+68	26	3.77
52+58	23	3.34
51+52	20	2.90
51+58	16	2.32
52+56	16	2.32
16+58	15	2.18
66+56	12	1.74
39+52	12	1.74
56+58	12	1.74

Age-Specific Prevalence of HPV Infection

HPV prevalence in different age groups is shown in Figure 3A. The results indicate that there was a significant difference in the infection rates among the different age groups (p < 0.05). Moreover, healthy women who were younger than 25 years had the highest HPV prevalence (20.92%, 64/306), whereas healthy women who were over 65 years had the lowest prevalence (10.97%, 154/1404).

Figure 3

Prevalence of HPV grouped by age. (A) Prevalence of HPV infection in different age groups; (B) Prevalence of single and multiple infection of HPV in different age groups; (C) Prevalence of HPV16, 51, 52, 56, and 58 in different age groups.

Significant differences were also found among the ten age groups in terms of single and multiple infections (p < 0.05). The trend was similar to the results presented earlier, wherein the prevalence of a single infection was much higher than that of a multiple infection, and the rate of single and multiple infections in healthy women aged less than 25 years old was significantly higher than that in other groups (p < 0.05) (Figure 3B).

Figure 3C reveals the variation tendency of the five most common HPV subtypes (HPV16, 51, 52, 56 and 58) among different age groups in healthy females. There were significant differences in the distribution of these five subtypes among the ten age groups (p < 0.05). In this study, the prevalence of HPV52 in each age group was much higher than that of the other four subtypes. Moreover, the positive rate of HPV16 was higher than that of HPV58 in all age groups, except for those aged 40–64 years.

Prevalence of HPV Grouped by Year

From 2016 to 2019, the annual prevalence of HPV infection among the healthy females was not significantly different (p > 0.05). In the past 4 years, cases with a single infection were consistently higher than those with multiple infections, and this trend was similar to that of overall infection discussed earlier (Table 3). Infections with each subtype from 2016 to 2019 are summarized in Table 4, and significant differences in distribution were seen for subtypes of HPV39, 58, 59, 68, 6 and 11 during this period (p < 0.05). Differences among the annual HPV-positive rates of healthy women aged 30–34 years were significant (p < 0.05) during the study period (Table 5).

Table 3

Comparison of Annual Total, Single, and Multiple Infection Rates of HPV from 2016 to 2019

Year	N	Total Infection, n (%)	Single, n (%)	Multiple, n (%)
2016	6398	796(12.44)	637(9.96)	159(2.49)
2017	8511	1016(11.94)	821(9.65)	195(2.29)
2018	7756	962(12.40)	787(10.15)	175(2.36)
2019	6771	812(11.99)	652(9.63)	160(2.36)
P-value		0.691	0.655	0.818

Table 4

Comparison the Prevalence of Different HPV Subtypes from 2016 to 2019

HPV subtype	2016, n (%)	2017, n (%)	2018, n (%)	2019, n (%)	P-value
HPV16	125(1.95%)	143(1.68%)	157(2.02%)	110(1.62%)	0.184
HPV18	56(0.88%)	51(0.60%)	52(0.67%)	40(0.59%)	0.153
HPV31	46(0.72%)	72(0.85%)	65(0.84%)	42(0.62%)	0.344
HPV33	33(0.52%)	36(0.42%)	31(0.40%)	32(0.47%)	0.735
HPV35	28(0.44%)	45(0.53%)	43(0.55%)	38(0.56%)	0.744
HPV39	69(1.08%)	67(0.79%)	81(1.04%)	48(0.71%)	0.044
HPV45	12(0.19%)	16(0.19%)	15(0.19%)	14(0.21%)	0.993
HPV51	72(1.13%)	85(1.00%)	91(1.17%)	98(1.45%)	0.081
HPV52	181(2.83%)	271(3.18%)	249(3.21%)	194(2.87%)	0.387
HPV56	71(1.11%)	80(0.94%)	91(1.17%)	69(1.02%)	0.497
HPV58	143(2.24%)	164(1.93%)	121(1.56%)	113(1.67%)	0.016
HPV59	34(0.53%)	40(0.47%)	66(0.85%)	46(0.68%)	0.013
HPV66	51(0.80%)	75(0.88%)	49(0.63%)	54(0.80%)	0.336
HPV68	39(0.61%)	44(0.52%)	56(0.72%)	74(1.09%)	0.000
HPV82	16(0.25%)	17(0.20%)	11(0.14%)	12(0.18%)	0.519
HPV (6+11)	21(0.33%)	48(0.56%)	18(0.23%)	32(0.47%)	0.005

Table 5

Comparison of HPV Infection Rate of Different Ages from 2016 to 2019

Age Groups (Years)	2016, n (%)	2017, n (%)	2018, n (%)	2019, n (%)	P-value
<25	13(22.41)	22(23.91)	17(20.00)	12(16.90)	0.726
25–29	62(17.87)	84(13.35)	108(16.22)	73(12.87)	0.096
30–34	107(15.71)	126(11.99)	110(10.61)	118(12.67)	0.018
35–39	96(10.67)	126(10.99)	110(9.80)	131(12.54)	0.234
40–44	125(11.94)	127(10.37)	117(11.93)	85(10.20)	0.429
45–49	119(10.59)	163(11.98)	152(13.15)	114(12.51)	0.286
50–54	118(10.92)	172(12.31)	120(11.58)	104(11.91)	0.753
55–59	72(13.82)	88(12.12)	114(14.92)	89(12.52)	0.376
60–64	52(14.40)	69(13.48)	67(12.52)	50(11.16)	0.541
≥65	32(11.51)	39(10.51)	47(12.57)	36(9.45)	0.564

Region-Specific Prevalence of HPV Infection

The prevalence of HPV infection in Beijing and in 12 different areas of China was compared separately (Table 6). The positive rate of HPV obtained in our study (Beijing) was significantly different from those in other areas (p < 0.05), except for Zhengzhou. Moreover, the highest rate of HPV infection was in Shaanxi Province (30.21%, 4263/14,111),29 whereas the lowest was in Zhengzhou City (12.09%, 6020/49,793),20 with an incidence rate similar to that observed in our study. The three most common genotypes of HPV infections, namely HPV 16, 52 and 58, were consistent among these regions, and the infection rates of these three subtypes among different regions are listed in Table 7.

Table 6

Comparison of HPV Infection Rates in Different Areas in China

Regions	Positive Number	Prevalence (%)	P-value
Beijing City	3586	12.18
Shanghai City19	6360	18.95	<0.001
Zhengzhou City20	6020	12.09	0.7
Wuhan City21	2436	17.68	<0.001
Xinjiang Province22	5287	14.02	<0.001
Jiangxi Province23	16,065	22.49	<0.001
Sichuan Province24	2641	24.11	<0.001
Yunnan Province25	3681	12.94	<0.001
Zhejiang Province26	197,367	20.54	<0.001
Guangdong Province27	5690	19.81	<0.001
Jiangsu Province28	16,775	26.92	<0.001
Shaanxi Province29	4263	30.21	<0.001
Shandong Province30	26,839	28.40	<0.001

Table 7

Prevalence of the Three Most Common HPV Genotypes of in Different Areas in China

Regions	N	Prevalence of Genotypes (%)
		HPV16	HPV52	HPV58
Beijing City	29,436	1.82	3.04	1.84
Shanghai City19	33,562	3.34	2.48	2.63
Zhengzhou City20	49,793	2.32	2.41	1.59
Wuhan City21	13,775	2.56	4.23	2.37
Xinjiang Province22	37,722	3.79	2.47	1.76
Jiangxi Province23	71,435	2.60	2.23	1.76
Sichuan Province24	10,953	3.33	3.74	3.21
Yunnan Province25	28,457	1.71	2.08	1.04
Zhejiang Province26	961,029	3.62	3.36	3.24
Guangdong Province27	28,730	5.01	2.63	1.91
Jiangsu Province28	62,317	5.06	5.09	3.14
Shaanxi Province29	14,111	10.71	5.24	5.43
Shandong Province30	94,489	5.77	5.14	3.46

Discussion

Our cross-sectional study evaluated the prevalence and genotype distribution of HPV in healthy women in China, with the aim of better understanding the epidemiological characteristics of HPV infection in Beijing. In this study, the prevalence of HPV infection among apparently healthy women in Beijing was 12.17%. In contrast, a recent research in Beijing has shown that 21.06% women were infected with HPV.31 The difference may have been mainly due to differences in the targeted populations screened, since the former study was based on all women who visited the hospital and received an HPV test, whereas our study was based on healthy women who received a routine medical examination.

In agreement with previous studies, we found that HPV16, 52, and 58 were the three most common subtypes in most regions in China.32,33 However, this result was distinct from that of a meta-analysis involving 1 million women from five continents, which indicated that the three most commonly identified HR-HPV genotypes worldwide were HPV16, 18, and 52.34 Although the majority of the 29,436 healthy females in this study were infected with a single HPV genotype, the possibility of multiple-type infections cannot be ignored. Multiple studies suggest that infections with multiple HPV genotypes might increase the length of persistent infection, and possibly further increase the risk of developing cervical cancer.35,36 Moreover, the three most common combinations of genotypes in multiple HPV infections were HPV16+52, HPV39+68, and HPV52+58. Interestingly, while HPV39 and 68 were not part of the three most commonly detected HPV genotypes, they were part of the common combinations of genotypes in multiple infections. This may be linked to the tendency of HPV68 infection rate to increase yearly; HPV39 infection was ranked sixth in prevalence in this study. However, this finding should be confirmed in larger cohorts.

In our study, HPV infection rates were different across different age groups, and this was also observed in previous reports.34,37 Our results showed that the prevalence of HPV significantly peaked below the age of 25 years, then showed a clearer trend of moderating in the middle age groups, and was a fluctuating variation of two little peaks after the age of 45. In terms of multiple infections, there was a gradual upward trend in the over 45 years old age groups. The higher HPV prevalence in young women might be related to more frequent sexual intercourse and inadequate awareness of safety.38,39 Moreover, they were likely more sensitive to HPV infections after the beginning of sexual behavior, as a result of immature immune protection; however, majority of infections are temporary.40 When women reach menopause, fluctuations in hormones cause a series of physiological and psychological changes, and immune dysfunction. It is possible that this places them in a more vulnerable category for HPV, owing to their diminishing ability of eliminating and inhibiting the virus.41,42 Therefore, HPV screening should be taken much more seriously for women over the age of 45. In some studies, researchers observed a bimodal curve (or an U‐shape curve) in the age‐specific HPV distribution.22,24,43 However, the phenomenon of the curve was not obvious in our study, and whether or not it exists in healthy women remains to be seen in follow-up studies.

While variations in the incidence of HPV infection were observed from 2016 to 2019, the results of this study showed no significant difference in the total, single, and multiple infections observed within four years, and suggested that the positive rate of HPV may tend to stabilize in this population. However, significant changes in the prevalence of some genotypes were observed during this period. Specifically, the prevalence of HPV58 tended to decrease, whereas that of HPV68 tended to increase and even ranked fourth in 2019. As an important part of common combinations of genotypes, it is also a noteworthy point. For future studies, we will continue to collect data, and observe the trend of positive rates of HPV and different genotypes over the years.

There is a great deal of research indicating that HPV infections are significantly region-specific.33,44,45 Compared with similar studies in other parts of China, the prevalence observed in our study was close to what has been reported in Zhengzhou, however, it was lower than that reported in the other 11 regions. While HPV16, 52 and 58 were consistently the three most frequently observed subtypes across regions, substantial differences in their prevalence in different regions were observed. This difference in the prevalence and genotype distribution of HPV might have been due to the different geographical and economic conditions, and lifestyles in the regions considered. Apart from that, the diversity in the study populations might have also contributed to the observed difference. Hence, a multicenter study should be developed to confirm the differences further.

Three types of vaccines, namely, the bivalent (HPV16/18), quadrivalent (HPV6/11/16/18), and 9-valent HPV vaccine (HPV6/11/16/18/31/33/45/52/58), have been approved by the National Medical Products Administration in China. The 9-valent HPV vaccine covers a broader swath and has a stronger protective effect against HPV. However, it has not been used extensively in mainland China, probably because it has a strict age limit that covers only young women (16–26 years old) and its price is higher than that of others. A study on patients with cervical cancer in West China found that the five most common subtypes of HPV were HPV16, 18, 58, 53, and 33.46 In our study, two of these genotypes (HPV16 and 58) also had higher prevalence in healthy women. In addition, in multiple HPV infections, HPV16+52 was the most common combination observed. Furthermore, the positive rate of HPV16/18/52/58 accounted for more than half of all HPV infections in our study. In the light of the above findings, it is necessary to make a few adjustments in the HPV genotypes being targeted by future vaccines. For example, an HPV vaccine that covers the four additional HR-HPV types (HPV16/18/52/58) should be developed and included in the approved vaccines, as this may be more conducive to the implementation of the HPV vaccination program in China.

Conclusions

In summary, this retrospective study found that the five most prevalent HPV genotypes among healthy women in Beijing were HPV52, 58, 16, 51, and 56. A significant prevalence peak of HPV appeared in younger women aged less than 25 years. Moreover, the age-specific and region-specific prevalence of HPV infection and distribution of genotypes in healthy women were also evaluated. To the best of our knowledge, this is the first study to evaluate HPV infection among healthy women in Beijing. Our findings can be considered as benefit for better understanding the epidemiological characteristics of HPV infection in Beijing. However, further large-scale and multicenter studies are needed to fully understand the prevalence and genotype distribution of HPV in healthy Chinese women. The limitations of the study must be acknowledged. There was a lack in analyzing the potential risk factors for HPV infection to safeguard patient privacy; the correlation between HPV genotypes and cervical cytology grade could not been elucidated. In future studies, we aim to expand the collection of sample information and analyze the correlation between cervical cytology grade and other potential risk factors such as smoking, drinking habit, education, occupation, marital status, age at sexual debut, contraceptive use, and number of sex partners in the last six months.