Skip to main content

ORIGINAL RESEARCH article

Front. Microbiol., 04 January 2023
Sec. Virology

Prevalence of human papillomavirus genotypes and related cervical morphological results in southern Hunan Province of China, 2018–2020: Baseline measures at a tertiary institution prior to mass human papillomavirus vaccination

Zhihua LanZhihua Lan1Jing ZhangJing Zhang1Hongtao LiHongtao Li1Rongfang HeRongfang He1Qiang ZhaoQiang Zhao1Fang Yang
Fang Yang2*
  • 1Department of Pathology, The First Affiliated Hospital, Hengyang Medical School of University of South China, Hengyang, China
  • 2Department of Anorectal Surgery of Traditional Chinese Medicine, The First Affiliated Hospital, University of South China, Hengyang, China

Background: Human papillomavirus (HPV) infection is closely correlated with cervical lesions. However, the HPV prevalence varies among different districts. This retrospective study investigated the HPV genotype distribution and its relationship with cervical lesions in southern Hunan Province.

Methods: The database at our Pathology Department was searched for HPV and morphological results during 2018–2020 were reviewed. A total of 49,955 gynaecological inpatients and outpatients, each of whom underwent HPV testing based on the amplification of L1 sequence and reverse dot blot hybridization, were included in this study. Available cytology and relevant histological examination results were reviewed. Enrolled cases were categorized into seven groups based on their age. Household registry and educational level were evaluated.

Results: Seven thousand two hundred eighty-six females were positive for HPV and the overall HPV positivity rate was 14.59%. The top five most prevalent HPV genotypes were HPV52, 16, 58, 53 and 51 (22.98%, 17.54%, 14.29%, 7.47%, and 5.70%, respectively). The HPV prevalence curve specific to the seven age groups showed a bimodal distribution. High school education or blow and rural residence served as risk factors for HPV infection. Single infection was the main type of HPV infection, and multiple infections occurred in 21.92% of the infected women. Of 3,148 cases had cytological results, 1,149 (36.50%) had abnormal cytological abnormalities. Among 2,833 participants with histological examination, 2000 (70.60%) had cervical abnormalities. Cytological and pathological abnormalities were mainly associated with infection with HPV16, 52 and 58. Further analysis found that HPV16, 58, 52 and 33 were the main genotypes associated with high-grade squamous lesions (HSIL) and that HPV16, 31, 33 and 58 were independent risk factors for HSIL (odd ratio [OR] = 1.70, 1.99, 2.59, 2.29; 95% CI = 1.41–2.03, 1.17–3.41, 1.88–3.59, 1.85–2.82; all p < 0.05). HPV16, 58, 52, 18 and 33 were the most frequent genotypes detected in squamous cell carcinoma (SCC) patients, and HPV16 and 18 were independent risk factors for cervical carcinomas (OR = 6.72;95% CI = 5.48–8.25; p < 0.001; OR = 1.67, 95% CI = 1.22–2.30; p = 0.001).

Conclusion: This retrospective study demonstrated the prevalence and the distribution characteristics of HPV infection and its correlation with cervical lesions in southern Hunan Province. The comprehensive results of this survey can guide HPV vaccine research to protect against some common genotypes in China.

Introduction

Cervical cancer, the world’s fourth most common cancer, is becoming a concerning public health problem and accounts for 40.1% of new global cases of cancer in eastern Asia (Sung et al., 2021). China, the largest developing country, reported 109,741 new cases of cervical cancer and 59,060 cervical cancer-related deaths in 2020 (Cao et al., 2021).

Human papillomavirus (HPV) is a small, nonenveloped double-stranded DNA virus. More than 200 genotypes have been identified, of which approximately 40 genotypes infect the genital tract (Bzhalava et al., 2015; Luo et al., 2020). Most HPV infections are transient and can be eradicated by the immune system. Persistent infection induces cervical lesions. Based on pathogenicity, HPVs are divided into a low-risk category, which includes HPVs that can lead to cervical low-grade squamous intraepithelial lesions (LSILs), and a high-risk category, which includes HPVs that are strong carcinogenic factors that can induce cervical high-grade squamous intraepithelial lesions (HSIL) and invasive squamous cell carcinoma (SCC). More than 90% of HSIL and SCCs are HPV-associated. Therefore, HPV vaccination has been provided to prevent cervical lesions.

Accordingly, the 2-valent vaccine against HPV16/18 and 4-valent vaccine against HPV6/11/16/18 provide 70% protection against cervical cancer. The protection rate increases to 90% in persons who receive the 9-valent vaccine targeting HPV6/11/16/18/31/33/45/52/58 (Yan et al., 2021). However, vaccine appointments are not easy, especially for the 9-valent vaccine, and there are strict age limitations in China. It is worth noting that HPV vaccines play a prophylactic but not therapeutic role in cervical epithelial lesions and that their adverse effects should not be ignored (Morimoto et al., 2015). Additionally, the geographical distributions of HPV genotypes vary from region to region. Although HPV16 is generally considered highly carcinogenic and the most prevalent genotype responsible for more than 50% of cervical cancers in most regions of the world (Bletsa et al., 2021), some eastern Asian countries, such as Korea, Japan and China, have demonstrated higher prevalence rates of HPV58 and HPV52 in women (Asato et al., 2004; Nah et al., 2017; Zhang et al., 2018).

Therefore, investigating and understanding the prevalence of HPV genotypes and their related cervical lesions in a specific region is of great importance for formulating public health strategies and guiding the application of the HPV vaccine. The first study on the prevalence of HPV in the Hengyang district demonstrated that the infection rate between April 2010 and March 2012 was 22.6% (Li et al., 2013). During 2012–2018, the overall positive rate was 18.71% (Luo et al., 2021). However, the correlation between HPV prevalence and cervical morphology in Hunan Province had not been intensively studied. The gynaecology department of the First Affiliated Hospital of University of South China, the largest medical institution in southern Hunan Province, performs HPV genotype detection in approximately 20,000 people per year. This retrospective investigation reviewed HPV detection results, the prevalence of genotypes and the characteristics of HPV-associated morphologic abnormalities to provide useful references for the detection, prevention and control of related cervical lesions in southern Hunan Province.

Materials and methods

Population and criteria

The database of our pathology department was searched for HPV test and related morphological examination results from HPV-positive patients from January 2018 to December 2020. We limited the interval between HPV testing and the receipt of associated morphological results (if available) in the same patient to within 180 days to ensure their correlation when analysing their relationship. The inclusion criteria were as follows: (1) history of sexual activity at any age; (2) had not been vaccinated against HPV; and (3) the first positive result from patients with several HPV test. The exclusion criteria included the following: (1) repeated results from the same patients; (2) prior physical, chemo- or radiation therapy for cervical lesions; or (3) history of conization or hysterectomy. Correlations with cytological and histological examination results after positive HPV test results were analysed.

Cervical sample collection

Cervical samples were collected by a gynaecologist via two sampling brushes. The brush for the HPV test was submerged in 2 ml 0.9% saline for subsequent analysis within 24 h. The sample for the cytology test was placed into cell storage solution (Guangzhou Anbiping Medical Company Technology Co., Ltd.) for cytological testing.

Human papillomavirus DNA genotyping

Human papillomavirus testing was carried out using the amplification of L1 sequence and reverse dot blot hybridization. HPV DNA was extracted from sampled exfoliated cervical cells and hybridized with the target HPV sequence on the hybridization chip, which included 18 high-risk HPV types (16, 18, 31, 33, 35, 45, 52, 66, 39, 58, 26, 51, 56, 59, 68, 82, 53, and 73) and 10 low-risk HPV types (6, 11, 40, 42, 44, 61, 83, 55, 43, and 81). Final result interpretation was complemented after chip washing and visualization. Chips that showed interactions with cells from the HPV-positive sample exhibited one or more blue spots at the corresponding position referring to the specific HPV genotype. The above procedures were performed with a 28 HPV GenoArray Diagnostic Kit (Guangzhou Anbiping Medical Company Technology Co., Ltd.) according to the manufacturer’s instructions.

ThinPrep cytology test

The Sedimentation Cell Prep Plus LBC (liquid-based cytology) Processor under the liquid-based preparation (LBP) system (LBP-2601, Guangzhou Anbiping Medical Company Technology Co., Ltd.) was used. Cells from the exfoliative sample were automatically sedimented onto a glass slide, forming a diagnostic area of 13 mm in diameter. Cytology evaluation was performed in conformity with the Bethesda 2014 criteria to identify negative intraepithelial lesion or malignancy (NILM); atypical squamous cells of undetermined significance (ASC-US); atypical squamous cells, cannot exclude high-grade squamous intraepithelial lesion (ASC-H); atypical glandular cells-not otherwise specified (AGC-NOS); atypical glandular cell-favour neoplasia (AGC-FN); adenocarcinoma (AC); low-grade squamous intraepithelial lesion (LSIL); high-grade squamous intraepithelial lesions (HSIL); and SCC.

Histological diagnosis

Cervical biopsy and other histological examination results were evaluated by two senior pathologists. Any cases with varying results were submitted to another senior doctor and confirmed by their discussion. The diagnosis was based on the 2020 World Health Organization (WHO) (Fifth Edition) classification criteria, including no neoplastic lesions, LSIL (low-grade dysplasia corresponding to cervical intraepithelial neoplasia grade I (CIN1)), HSIL (high-grade dysplasia corresponding to CIN2 and CIN3), SCC, AC (AC in situ and invasive AC) and adenosquamous carcinoma (ASC). If two or more pathological examination results were available for the same patient, only the most serious result was included.

Statistical analysis

Human papillomavirus, cytological, histological and clinical data were analysed in Excel 2020 and R software (x64 4.1.2).

The HPV positivity rate, single genotype infection rate and composition of multiple genotype infections were evaluated. Analysis was performed on the basic population information, as well as household registration and educational level. The HPV-positive cases were also stratified by age group (<20, 20–29, 30–39, 40–49, 50–59, 60–69, and ≥70 years) and genotype for comparison. Correlations of different genotypes with cytological categories (NILM, ASC-US, ASC-H, LSIL, HSIL, SCC, AGC-NOS and AGC-FN) and histological categories (negative, LSIL, HSIL, SCC and AC) were calculated.

The odds ratios and relative 95% confidence intervals for each HPV genotype in HSIL and cervical cancers were calculated and compared to those of other genotypes. p < 0.05 was considered statistically significant.

Results

Gross human papillomavirus prevalence

Data from 49,955 HPV screening tests performed between 2018 and 2020 were included. 26,268 cases were from rural areas and the other 23,687 were from urban regions. The age of all the patients ranged from 14 to 91 years, with an average of 41.9 ± 11.1 years. A total of 7,286 patients, accounting for 14.59% of the total population, were HPV positive. Among them, 5,507 (75.58, 95% CI, 74.58–76.57%) were infected with high-risk genotypes, and 1779 (24.42, 95% CI, 23.43–25.42%) were infected with only low-risk genotypes. The positive rate in rural women were 16.36%, which was significantly higher than 12.62% in urban females (p < 0.05). Additionally, the infective rate in cases with college degree or above were 11.87%, which was significantly lower than 14.99% in females without college education (p < 0.05).

Human papillomavirus prevalence by age group

The HPV prevalence curve specific to the seven age groups showed a bimodal distribution (Figure 1). The HPV infection rate of the <20 years group was 27.32% (95% CI, 23.11–31.84%) and had the highest peak. In the 20–29, 30–39 and 40–49 years age groups, the infection rates decreased were 13.14% (95% CI, 12.35–13.96%), 13.04% (95% CI, 12.48–13.61%) and 13.26% (95% CI, 12.75–13.79%), respectively. With increasing age, the prevalence rate increased to 17.25% (95% CI, 16.49–18.03%) in the 50–59 years group, with a second peak (22.93, 95% CI, 21.27–24.64%) in the 60–69 years group. The ≥70 years group accounted for 19.73% (95% CI, 16.40–23.41%) of the infected population. Regardless of age group, the positive patients were mainly infected with high-risk HPV genotypes (Table 1). A total of 22.6% (95% CI, 15.33–31.35%) of the positive patients in the <20 years group had only low-risk HPV genotypes. The rates in the 20–29, 30–39, 40–49, 50–59, 60–69 and ≥70 years groups were 15.71% (95% CI, 13.40–18.25%), 14.99% (95% CI, 13.36–16.73%), 16.09% (95% CI, 14.58–17.69%), 13.23% (95% CI, 11.61–14.98%), 9.30% (95% CI, 7.03–12.02%) and 9.71% (95% CI, 4.75–17.13%), respectively.

FIGURE 1
www.frontiersin.org

Figure 1. The prevalence of HPV infection in different age groups from 2018 to 2020.

TABLE 1
www.frontiersin.org

Table 1. HPV prevalence in different age groups.

Of all the HPV genotypes, HPV52 was the most prevalent genotype found in 1674 patients, accounting for 22.98% (95% CI, 22.01, 23.96%) of infections. However, in the groups aged <20 and 60–69 years, HPV16 ranked first (Table 2; Figure 2A). Among the low-risk genotypes, HPV81 had the highest infection rate in the 40–49, 50–59, 60–69 and ≥70 years age groups. In addition, HPV06 in the groups aged <20 and 20–29 years and HPV44 in the group aged 30–39 years had the highest infection rates among the low-risk genotypes (Figure 2B).

TABLE 2
www.frontiersin.org

Table 2. Distribution of HPV genotypes in different age groups.

FIGURE 2
www.frontiersin.org

Figure 2. The distribution of HPV genotypes in different age groups. Each bar represents the number of infections due to corresponding HPV genotypes. Black, yellow, pink, red, blue, green and violet indicate the numbers of infections in the <20, 20–29, 30–39, 40–49, 50–59, 60–69 and >70 years groups for the corresponding genotypes, respectively. (A) High-risk HPV infection. (B) Low-risk HPV infection.

Distributions of single and multiple HPV infections in different genotype and age groups

Regarding HPV26, 45, 56, 59, 82, 06, 40, 42 43 and 83, combinations with other genotypes were more common than single infections (Figure 3). However, single HPV infection was the most prevalent pattern in the different age groups and occurred in 5686 (78.04%; 95% CI, 77.07–78.99%) patients (Table 3). The most common multiple-infection situations were dual infections (16.50%; 95% CI, 15.62–17.34), followed by 3-strain infections (3.94%; 95% CI, 3.50–4.41%), 4-strain infections (1.04%; 95% CI, 0.82–1.30%), 5-strain infections (0.26%; 95% CI, 0.16–0.41%), and six-strain or more infections (0.25%; 95% CI, 0.15–0.39%; Figure 4A). The top 10 dual HPV combinations were as follows: 16/52, 16/58, 52/58, 52/53, 51/52, 18/52, 16/53, 52/81, 51/58 and 52/68 (Figure 4B). HPV81 was the only low-risk genotype among these 10 combinations. The multiple HPV infection rates in different age groups were different. The highest rate was in the <20 years group and reached 43.48% (95% CI, 34.26–53.04%; p = 0.19). As age increased, the rate decreased to 24.89% (95% CI, 22.10–27.84%), 19.02% (95% CI, 17.22–20.91%) and 18.53% (95% CI, 16.93–20.22%) in the 20–29, 30–39 and 40–49 years groups, respectively. However, the elderly groups (50–59, 60–69 and ≥70 years) accounted for 23.39% (95% CI, 21.25–25.44%), 30.77% (95% CI, 26.96–34.78%) and 28.16% (95% CI, 19.73–37.87%) of multiple HPV infections, respectively. The distribution also showed a bimodal pattern.

FIGURE 3
www.frontiersin.org

Figure 3. Distribution of HPV genotypes during 2018–2020. Each bar represents the number of infections due to corresponding HPV genotypes. Red indicates the number of cases with only single-type infections for the corresponding HPV genotype, and yellow indicates the number of cases due to multiple infections. (A) High-risk HPV infection. (B) Low-risk HPV infection.

TABLE 3
www.frontiersin.org

Table 3. Distribution of single-type and multiple-type infection among different age groups.

FIGURE 4
www.frontiersin.org

Figure 4. (A) Distribution of multitype HPV infection. Each part represents the number of cases due to multiple HPV infections. Orange, red, olivine, green, blue and violet indicate the number of dual, 3-strain, 4-strain, 5-strain, and six-strain or more infections, respectively. (B) The green bar represents the number of the top 10 dual HPV infections.

Distribution of HPV genotypes according to different cytological results

To further analyse the 7,286 HPV-positive samples, cytological results from 3,148 patients in our database were analysed. The distribution of liquid-based cytological examination results of different HPV genotypes is shown in Table 4 and Figure 5. Among such patients, 1,149 (36.50%; 95% CI, 34.81–38.21%) had abnormal cytological abnormalities and 1999 (63.50%; 95% CI, 61.79–65.19%) had negative results. HPV52 was the most common genotype, accounting for 25.26% (95% CI, 23.37–27.23%) of NILM, 26.13% (95% CI, 21.99–30.60) of ASC-US and 23.32% (95% CI, 18.25–29.03) of LSIL. However, HPV16 was the leading genotype found in 44.36% (95% CI, 35.75–53.22%; p = 0.22) of ASC-H, 46.86% (95% CI, 41.27–50.50%; p = 0.29) of HSIL and 71.43% (95% CI, 41.90–90.61%; p = 0.18) of SCC, followed by HPV58. Among patients without high-risk HPV infection, HPV81 was the main genotype, accounting for 15 (3.54%; 95% CI, 2.01–5.81%) cases of ASC-US, accompanied by HPV61, accounting for 4 (1.58%; 95% CI, 0.43–4.00%) cases of LSIL.

TABLE 4
www.frontiersin.org

Table 4. Distribution of HPV genotypes in different cytology results.

FIGURE 5
www.frontiersin.org

Figure 5. HPV genotype infection rates in the NILM, ASC-US, LSIL, ASC-H, HSIL and SCC groups. Each bar represents the number of infections due to corresponding HPV genotypes. Green, violet, pink, blue, red and yellow represent the numbers of NILM, ASC-US, LSIL, ASC-H, HSIL and SCC cases attributable to the corresponding genotype, respectively. (A) High-risk HPV genotypes with cytological results. (B) Low-risk HPV genotypes with cytological results.

Distribution of HPV genotypes according to different pathological results

A total of 2,833 participants were referred for histological examination. Among them, 2,000 patients had pathological abnormalities (70.60, 95% CI, 68.88–72.27%). The distribution of pathological results of cases due to different HPV genotypes is shown in Table 5 and Figure 6. Of the 760 cases of LSIL, HPV52 was the dominant genotype (28.68%; 95% CI, 25.49–32.04%), followed by HPV16 (15.00%; 95% CI, 12.54–17.74%), HPV58 (14.47%; 95% CI, 12.05–17.18%) and HPV53 (9.74%; 95% CI, 7.72–12.07%). Of the 714 cases of HSIL, HPV16 was the dominant genotype (39.22%; 95% CI, 35.62–42.91%), followed by HPV58 (25.77%; 95% CI, 22.60–29.14%), HPV52 (21.71%; 95% CI, 18.74–24.92%) and HPV33 (10.08%; 95% CI, 7.97–12.53%). Of the 481 cases of SCC, HPV16 was the dominant genotype (67.98%; 95% CI, 63.61–72.13%), followed by HPV58 (11.02%; 95% CI, 8.36–14.16%), HPV52 (8.32%; 95% CI, 6.01–11.15%) and HPV18 (7.28%; 95% CI, 5.12–9.97%). Of the 44 cases of AC, HPV18 (50.00%; 95% CI, 34.56–65.43%; p = 1.00) and HPV16 (43.18%; 95% CI, 28.35–58.97; p = 0.45) were the main genotypes. As shown in Table 6, patients infected with HPV16, 31, 33 and 58 had increased risks of HSIL (OR = 1.70, 95% CI = 1.41–2.03; OR = 1.99, 95% CI = 1.17–3.41; OR = 2.59, 95% CI = 1.88–3.59; OR = 2.29, 95% CI = 1.85–2.82). Patients with HPV16 and 18 were at higher risk of cervical carcinoma (including SCC, AC and ASC; OR = 6.72, 95% CI = 5.48–8.25; OR = 1.67, 95% CI = 1.22–2.30).

TABLE 5
www.frontiersin.org

Table 5. Distribution of HPV genotypes in different pathological result.

FIGURE 6
www.frontiersin.org

Figure 6. HPV genotype infection rates in the negative, LSIL, HSIL and SCC groups. Each bar represents the number of infections due to the corresponding HPV genotypes. Grey, violet, red and black represent negative results, LSIL, HSIL and SCC for the corresponding genotype, respectively. (A) High-risk HPV genotypes with pathological results. (B) Low-risk HPV genotypes with pathological results.

TABLE 6
www.frontiersin.org

Table 6. Odds ratios of high-risk HPV genotypes in HSIL and cervical cancers.

Discussion

Action toward achieving the global elimination of cervical cancer has been proposed by the World Health Organization (WHO) since 2018 (Kumar et al., 2019). Free screening for cervical cancer and breast cancer has been offered to rural women in China since 2009. At the beginning of 2022, urban women were also enrolled in the programme. As a sensitive method, HPV testing is usually used in combination with cytology for screening for cervical lesions (Hamers et al., 2022).

The fact that the HPV prevalence varies considerably on the basis of geography and population has prompted many investigations into regional epidemical strategies. Disparities also exist within countries.

The overall HPV infection rate in our study was 14.59%, which was lower than those in Guizhou (16.95), Sichuan (23.84%), Fujian (20.57%), Guangdong (19.81%), Shanghai (17.92%), Shandong (28.4%), Liaoning (16.1%), Henan (19.7%), and Jiangsu (26.92%) but higher than those in Xinjiang (14.02%) and Yunnan (12.9%; Li et al., 2016, 2020; Sun et al., 2017; Zhao et al., 2018; Chen et al., 2019; Jiang et al., 2019; Wang et al., 2019, 2022; Zhang et al., 2019, 2020; Luo et al., 2020). Compared to previous study, the HPV positive rate decreased significantly, which may be due to increasing health awareness led to increase participation in cervical screen (Luo et al., 2021). The distribution of the age-specific HPV infection rate presented a bimodal pattern, which was in accordance with the results of the majority of previous studies. The <20 years age group had the highest infection rate, possibly because new sexual behaviour increases the risk of virus exposure and an undeveloped immune system has an insufficient capability for virus clearance (Bergqvist et al., 2021; Yan et al., 2021). The 60–69 years age group had the second highest HPV infection rate. The reason could be a decreased ability to clear recent infection due to age-related immune senescence and persistent infections in older women from earlier exposures (Kang et al., 2014). Notably, the numbers of participants aged <20 and ≥60 were significantly less than those in the other groups. In fact, only a small minority of women under the age of 20 reported sexual behaviour. Fewer elderly women aged 60 years or older received HPV testing in the clinic, possibly because of reluctance due to an insufficient understanding or a lack of symptoms (He and He, 2020). Consistent with previous study (Yang et al., 2022), the current study shows that women with college education or from urban areas had lower risk of HPV infection. Higher education and convenient urban medical service may improve the health attention and practice of the population.

Single-type HPV infection and high-risk HPV infection in our study were the most common types of infections. The most common genotype worldwide is HPV16 (Crow, 2012), similar to that in northernmost China (Jiang et al., 2019; Zhang et al., 2020; Wang et al., 2022). However, the most prevalent genotype in southernmost China is HPV52 (Li et al., 2016, 2020; Chen et al., 2019; Zhang et al., 2019; Luo et al., 2020). In the south of Hunan province, the most common genotype has been changed from HPV 16 to HPV 52 (Li et al., 2013; Luo et al., 2021). Regardless of whether HPV16 or 52 ranks first, the top three genotypes in China are HPV16, 52 and 58 (Li et al., 2016, 2020; Sun et al., 2017; Zhang et al., 2018, 2019, 2020; Chen et al., 2019; Jiang et al., 2019; Wang et al., 2019, 2022; Luo et al., 2020). As the second most carcinogenic genotype in the world, HPV18 ranks only fourth in some Chinese areas or is not even in the top five in other provinces (Li et al., 2016; Zhang et al., 2018). Moreover, HPV53 is the fourth or fifth most prevalent genotype in most areas of China (Sun et al., 2017; Chen et al., 2019, 2021; Jiang et al., 2019; Wang et al., 2019, 2022; Zhang et al., 2019, 2020; Li et al., 2020; Luo et al., 2020), and is the most prevalent genotype among healthy women in other Asian countries (Ouh et al., 2018). Our study also demonstrated that HPV81 was the most common low-risk type, with the same overwhelming advantage in other Chinese provinces (Zhao et al., 2018; Chen et al., 2019; Zhang et al., 2019; Luo et al., 2020; Yan et al., 2021).

The multiple-infection rate in HPV-positive patients in the different age groups also showed a bimodal distribution. The peak occurred in the <20 years group, and the second highest rate was observed in the 60–69 years group, similar to the distribution of the age-specific HPV infection rate and likely explained by the same reasons (Li et al., 2020). Patients with multiple HPV infections had a longer duration of HPV infection and an increased risk of cervical lesions compared to patients with a single HPV infection (Kim et al., 2021). High-risk HPVs were also the predominant genotypes in coinfections. Among the top 10 dual infections, HPV82 was the only low-risk genotype. However, the oncogenic potential of each HPV genotype under coinfection conditions could not be assessed accurately because interactions and competition between the genotypes were unclear (Bernard et al., 2013).

Cytology screening in patients with a high-risk HPV genotype is the conventional method (Rebolj et al., 2022). The distribution of different cervical cytology results in each HPV genotype group, especially the high-risk genotypes, is helpful in evaluation of the risk of morphological changes. Our investigation showed that HPV16, 52 and 58 accounted for the top three genotypes resulting in cervical cytopathologic abnormalities, which was in accordance with data from other studies (Luo et al., 2020; Yan et al., 2021). Compared to the other genotypes, HPV16 was the major species in the ASC-H, HSIL and SCC cytology samples. Similar to other observations, HPV52 was the dominant genotype in cases with negative cytology and ASC-US (Li et al., 2020; Yan et al., 2021). Some HSIL were related to uncommon genotypes, such as HPV31 and HPV33. Approximately 16.67% of 126 patients with HPV33 had HSIL. In contrast, HSIL occurred in only 18 (approximately 9.52%) HPV18-positive cases. Colposcopy is recommended for women with HPV16/18, even if NILM is absent according to the American Society for Colposcopy and Cervical Pathology (ASCCP; Perkins et al., 2021). Based on our study, it is also appropriate to perform colposcopy screening in women with HPV52/58/31/33 in China.

The distribution of different pathological results for each genotype was similar to the distribution of cytopathological results. HPV16, 58, and 52 were also the major genotypes in women with histological abnormalities. Relative risk analysis demonstrated that women infected with HPV16, 31, 33, and 58 had an increased risk of HSIL compared with the other genotypes, which supports the above recommendation for additional colposcopy examinations in the population in China (Long et al., 2018; Yan et al., 2021). In our study, HPV16 and 18 were also the top two oncogenic genotypes associated with cervical cancer, consistent with the results of a previous study (Yan et al., 2021).

According to the guidelines from the American Cancer Society, HPV testing in combination with cytology or cytology alone should be phased out once full access to primary HPV examination for cervical cancer screening is available without barriers (Fontham et al., 2020). Therefore, HPV testing and prevention will play a more important role in epidemiological management in different regions. HPV16/18 were detected in 75.89% of cervical cancer samples in our study, suggesting that approximately 70% of cervical carcinomas could be prevented by bivalent or quadrivalent HPV vaccination (Gonzalez-Bosquet et al., 2020). Moreover, only 24 cases, accounting for 4.51% of cervical cancer patients, were found to be caused by high-risk genotypes not covered by the 9-valent vaccine. China and other 193 countries made the commitment to make sure that 90% of girls will be fully vaccinated by the age of 15 years in 2030 (Xia et al., 2020). However, the application of this vaccine is limited by low availability of imported vaccine and a long vaccination process (Yan et al., 2021). Although China has not integrated HPV vaccination into State immune programme, some provinces such as Fujian, Zhejiang and Guangdong have implemented free administration of domestic bivalent vaccine for girls under 15 years old since 2022. As a populous developing country, China’s national campaign on HPV vaccination is restricted by shortage of vaccine supply and high price of imported vaccine. Therefore, the development of sufficient and inexpensive domestic preventive HPV vaccines covering HPV16, 18, 31, 33, 39, 51, 52, 53, and 58 should be accelerated to achieved the goal.

The major limitation of our study was the fact that participants were enrolled from a single institution and not by randomized sampling of the population; therefore, the results might not be representative of all women in southern Hunan Province. In addition, some HPV-positive patients were lost to follow-up (e.g., attended other hospitals for morphological examination) and thus were not included in the morphological study, which resulted in a selection bias.

Conclusion

This study analysed the prevalence of HPV and the age-specific HPV genotype distribution during 2017–2020 in southern Hunan Province. Moreover, the relationships between HPV genotypes and cervical morphological lesions were also analysed. Overall, HPV prevalence was 14.59% showing a decreased trend compared to previous study in the region. Bimodal distribution of age-specific HPV infection rate demonstrated females aged ≤20 or >60 years were susceptible to HPV. Significantly, HPV 52, 16 and 58 were main genotypes occurred in cervical abnormalities. Additionally, Infection of HPV 31, 33, 39, 51, and 53 should not be ignored in the follow-up management. Our study provides epidemiological data for the development of domestic HPV vaccines and future national campaign on HPV vaccination.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

The studies involving human participants were reviewed and approved by the Ethics Committee of the First Affiliated Hospital of University of South China. Written informed consent from the participants’ legal guardian/next of kin was not required to participate in this study in accordance with the national legislation and the institutional requirements.

Author contributions

ZL designed the study and drafted the manuscript. HL performed the HPV analyses. JZ, RH, and QZ collected and reviewed the data. FY performed the statistical analysis and reviewed the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by The First Affiliated Hospital, Hengyang Medical School, University of South China.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Abbreviations

HPV, human papillomavirus; LBC, liquid-based cytology; LBP, liquid-based preparation; NILM, intraepithelial lesion or malignancy; ASC-US, atypical squamous cells of undetermined significance; ASC-H, atypical squamous cells, not excluding high-grade squamous intraepithelial lesion; AGC-NOS, atypical glandular cell-not otherwise specified; AGC-FN, atypical glandular cell-favour neoplasia; AC, adenocarcinoma; LSIL, low-grade squamous intraepithelial lesion; HSIL, high-grade squamous intraepithelial lesion; CIN, cervical intraepithelial neoplasia; CI, confidence interval; WHO, World Health Organization.

References

Asato, T., Maehama, T., Nagai, Y., Kanazawa, K., Uezato, H., and Kariya, K. (2004). A large case-control study of cervical cancer risk associated with human papillomavirus infection in Japan, by nucleotide sequencing-based genotyping. J. Infect. Dis. 189, 1829–1832. doi: 10.1086/382896

PubMed Abstract | CrossRef Full Text | Google Scholar

Bergqvist, L., Kalliala, I., Aro, K., Auvinen, E., Jakobsson, M., Kiviharju, M., et al. (2021). Distribution of HPV genotypes differs depending on Behavioural factors among young women. Microorganisms. 9:750. doi: 10.3390/microorganisms9040750

PubMed Abstract | CrossRef Full Text | Google Scholar

Bernard, E., Pons-Salort, M., Favre, M., Heard, I., Delarocque-Astagneau, E., Guillemot, D., et al. (2013). Comparing human papillomavirus prevalences in women with normal cytology or invasive cervical cancer to rank genotypes according to their oncogenic potential: a meta-analysis of observational studies. BMC Infect. Dis. 13:73. doi: 10.1186/1471-2334-13-373

PubMed Abstract | CrossRef Full Text | Google Scholar

Bletsa, G., Zagouri, F., Amoutzias, G. D., Nikolaidis, M., Zografos, E., Markoulatos, P., et al. (2021). Genetic variability of the HPV16 early genes and LCR. Present and future perspectives. Expert. Rev. Mol. Med. 23:e19. doi: 10.1017/erm.2021.18

PubMed Abstract | CrossRef Full Text | Google Scholar

Bzhalava, D., Eklund, C., and Dillner, J. (2015). International standardization and classification of human papillomavirus types. Virology 476, 341–344. doi: 10.1016/j.virol.2014.12.028

PubMed Abstract | CrossRef Full Text | Google Scholar

Cao, W., Chen, H. D., Yu, Y. W., Li, N., and Chen, W. Q. (2021). Changing profiles of cancer burden worldwide and in China: a secondary analysis of the global cancer statistics 2020. Chin. Med. J. (Engl) 134, 783–791. doi: 10.1097/CM9.0000000000001474

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, R., Fu, Y., You, B., Li, Y., Yao, Y., Wang, X., et al. (2021). Clinical characteristics of single human papillomavirus 53 infection: a retrospective study of 419 cases. BMC Infect. Dis. 21:1158. doi: 10.1186/s12879-021-06853-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Z., Li, Q., Huang, Q., Liu, H., Jiang, H., Chen, Z., et al. (2019). Characteristics of human papillomaviruses distribution in Guizhou Province, China. Virol. J. 16:123. doi: 10.1186/s12985-019-1239-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Crow, J. M. (2012). HPV: the global burden. Nature 488, S2–S3. doi: 10.1038/488S2a

CrossRef Full Text | Google Scholar

Fontham, E., Wolf, A., Church, T. R., Etzioni, R., Flowers, C. R., Herzig, A., et al. (2020). Cervical cancer screening for individuals at average risk: 2020 guideline update from the American Cancer Society. CA Cancer J. Clin. 70, 321–346. doi: 10.3322/caac.21628

PubMed Abstract | CrossRef Full Text | Google Scholar

Gonzalez-Bosquet, E., Gibert, M., Serra, M., Hernandez-Saborit, A., and Gonzalez-Fernandez, A. (2020). Candidate HPV genotypes not included in the 9-valent vaccine for prevention of CIN 2-3. Int. J. Gynecol. Cancer 30, 954–958. doi: 10.1136/ijgc-2019-001069

CrossRef Full Text | Google Scholar

Hamers, F. F., Poullié, A. I., and Arbyn, M. (2022). Updated evidence-based recommendations for cervical cancer screening in France. Eur. J. Cancer Prev. 31, 279–286. doi: 10.1097/CEJ.0000000000000701

PubMed Abstract | CrossRef Full Text | Google Scholar

He, L., and He, J. (2020). Attitudes towards HPV self-sampling among women in Chengdu, China: a cross-sectional survey. J. Med. Screen. 27, 201–206. doi: 10.1177/0969141319895543

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, L., Tian, X., Peng, D., Zhang, L., Xie, F., Bi, C., et al. (2019). HPV prevalence and genotype distribution among women in Shandong Province, China: analysis of 94,489 HPV genotyping results from Shandong's largest independent pathology laboratory. PLoS One 14:e0210311. doi: 10.1371/journal.pone.0210311

PubMed Abstract | CrossRef Full Text | Google Scholar

Kang, L. N., Castle, P. E., Zhao, F. H., Jeronimo, J., Chen, F., Bansil, P., et al. (2014). A prospective study of age trends of high-risk human papillomavirus infection in rural China. BMC Infect. Dis. 14:96. doi: 10.1186/1471-2334-14-96

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, M., Park, N. J., Jeong, J. Y., and Park, J. Y. (2021). Multiple human papilloma virus (HPV) infections are associated with HSIL and persistent HPV infection status in Korean patients. Viruses 13:1342. doi: 10.3390/v13071342

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumar, P., Gupta, S., Das, A. M., and Das, B. C. (2019). Towards global elimination of cervical cancer in all groups of women. Lancet Oncol. 20:e237. doi: 10.1016/S1470-2045(19)30170-6

CrossRef Full Text | Google Scholar

Li, H., Li, P., Huang, L., Sun, L., Ren, H., and Li, P. (2020). Prevalence characteristics of cervical human papillomavirus (HPV) infection in the Zhoupu District, Shanghai City, China. Virol. J. 17:84. doi: 10.1186/s12985-020-01352-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Z., Liu, F., Cheng, S., Shi, L., Yan, Z., Yang, J., et al. (2016). Prevalence of HPV infection among 28,457 Chinese women in Yunnan Province, Southwest China. Sci. Rep. 6:21039. doi: 10.1038/srep21039

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, H., Zhang, J., Chen, Z., Zhou, B., and Tan, Y. (2013). Prevalence of human papillomavirus genotypes among women in Hunan Province. Eur. J. Obstet. Gynecol. Reprod. Biol. X 170, 202–205. doi: 10.1016/j.ejogrb.2013.05.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Long, W., Yang, Z., Li, X., Chen, M., Liu, J., Zhang, Y., et al. (2018). HPV-16, HPV-58, and HPV-33 are the most carcinogenic HPV genotypes in southwestern China and their viral loads are associated with severity of premalignant lesions in the cervix. Virol.J. 15:94. doi: 10.1186/s12985-018-1003-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Luo, L. P., He, P., Liu, Q. T., Jiang, Y. H., Zhang, Y. N., Li, Q. Z., et al. (2021). Prevalence and genotype distribution of HPV infection among 214,715 women from southern China, 2012-2018: baseline measures prior to mass HPV vaccination. BMC Infect. Dis. 21:328. doi: 10.1186/s12879-021-06019-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Luo, Q., Jiang, N., Wu, Q., Wang, J., and Zhong, J. (2020). Prevalence and genotype distribution of HPV and cervical pathological results in Sichuan Province, China: a three years surveys prior to mass HPV vaccination. Virol. J. 17:100. doi: 10.1186/s12985-020-01366-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Morimoto, A., Ueda, Y., Egawa-Takata, T., Yagi, A., Terai, Y., Ohmichi, M., et al. (2015). Effect on HPV vaccination in Japan resulting from news report of adverse events and suspension of governmental recommendation for HPV vaccination. Int. J. Clin. Oncol. 20, 549–555. doi: 10.1007/s10147-014-0723-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Nah, E. H., Cho, S., Kim, S., and Cho, H. I. (2017). Human papillomavirus genotype distribution among 18,815 women in 13 Korean cities and relationship with cervical cytology findings. Ann. Lab. Med. 37, 426–433. doi: 10.3343/alm.2017.37.5.426

PubMed Abstract | CrossRef Full Text | Google Scholar

Ouh, Y. T., Min, K. J., Cho, H. W., Ki, M., Oh, J. K., Shin, S. Y., et al. (2018). Prevalence of human papillomavirus genotypes and precancerous cervical lesions in a screening population in the Republic of Korea, 2014-2016. J. Gynecol. Oncol. 29:e14. doi: 10.3802/jgo.2018.29.e14

PubMed Abstract | CrossRef Full Text | Google Scholar

Perkins, R. B., Guido, R. L., Saraiya, M., Sawaya, G. F., Wentzensen, N., Schiffman, M., et al. (2021). Summary of current guidelines for cervical cancer screening and Management of Abnormal Test Results: 2016-2020. J Womens Health (Larchmt) 30, 5–13. doi: 10.1089/jwh.2020.8918

PubMed Abstract | CrossRef Full Text | Google Scholar

Rebolj, M., Mathews, C. S., and Denton, K., HPV Pilot Steering Group (2022). Cytology interpretation after a change to HPV testing in primary cervical screening: observational study from the English pilot. Cancer Cytopathol. 130, 531–541. doi: 10.1002/cncy.22572

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, P., Song, Y., Ruan, G., Mao, X., Kang, Y., Dong, B., et al. (2017). Clinical validation of the PCR-reverse dot blot human papillomavirus genotyping test in cervical lesions from Chinese women in the Fujian province: a hospital-based population study. J. Gynecol. Oncol. 28:e50. doi: 10.3802/jgo.2017.28.e50

PubMed Abstract | CrossRef Full Text | Google Scholar

Sung, H., Ferlay, J., Siegel, R. L., Laversanne, M., Soerjomataram, I., Jemal, A., et al. (2021). Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J. Clin. 71, 209–249. doi: 10.3322/caac.21660

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., Song, Y., Wei, X., Wang, G., Sun, R., Wang, M., et al. (2022). Prevalence and distribution of human papillomavirus genotypes among women attending gynecology clinics in northern Henan Province of China. Virol. J. 19:6. doi: 10.1186/s12885-020-06935-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, J., Tang, D., Wang, K., Wang, J., Zhang, Z., Chen, Y., et al. (2019). HPV genotype prevalence and distribution during 2009-2018 in Xinjiang, China: baseline surveys prior to mass HPV vaccination. BMC Womens Health 19:90. doi: 10.1186/s12905-019-0785-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Xia, C. F., Qiao, Y. L., Zhang, Y., and Zhao, F. H. (2020). WHO's global strategy of cervical cancer elimination and the challenges and initiatives in China. Zhonghua Yi Xue Za Zhi 100, 3484–3488. doi: 10.3760/cma.j.cn112137-20200909-02606

PubMed Abstract | CrossRef Full Text | Google Scholar

Yan, X., Shen, L., Xiao, Y., Wang, Q., Li, F., and Qian, Y. (2021). Prevalence, characteristics, and distribution of HPV genotypes in women from Zhejiang Province, 2016-2020. Virol. J. 18:208. doi: 10.1186/s12985-021-01676-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, D., Zhang, J., Cui, X., Ma, J., Wang, C., and Piao, H. (2022). Status and epidemiological characteristics of high-risk human papillomavirus infection in multiple centers in Shenyang. Front. Microbiol. 13:985561. doi: 10.3389/fmicb.2022.985561

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, C., Cheng, W., Liu, Q., Guan, Q., and Zhang, Q. (2019). Distribution of human papillomavirus infection: a population-based study of cervical samples from Jiangsu Province. Virol. J. 16:67. doi: 10.1186/s12985-019-1175-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, C., Huang, C., Zheng, X., and Pan, D. (2018). Prevalence of human papillomavirus among Wenzhou women diagnosed with cervical intraepithelial neoplasia and cervical cancer. Infect. Agent. Cancer. 13:37. doi: 10.1186/s13027-018-0211-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, J., Zhang, D., Yang, Z., Wang, X., and Wang, D. (2020). The role of human papillomavirus genotyping for detecting high-grade intraepithelial neoplasia or cancer in HPV-positive women with normal cytology: a study from a hospital in northeastern China. BMC Cancer 20:443. doi: 10.1186/s12885-020-06935-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, P., Liu, S., Zhong, Z., Hou, J., Lin, L., Weng, R., et al. (2018). Prevalence and genotype distribution of human papillomavirus infection among women in northeastern Guangdong Province of China. BMC Infect. Dis. 18:204. doi: 10.1186/s12879-018-3105-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: human papillomavirus, genotypes, morphological examination, cervical intraepithelial neoplasm, cervical cytology, histopathology

Citation: Lan Z, Zhang J, Li H, He R, Zhao Q and Yang F (2023) Prevalence of human papillomavirus genotypes and related cervical morphological results in southern Hunan Province of China, 2018–2020: Baseline measures at a tertiary institution prior to mass human papillomavirus vaccination. Front. Microbiol. 13:1094560. doi: 10.3389/fmicb.2022.1094560

Received: 10 November 2022; Accepted: 08 December 2022;
Published: 04 January 2023.

Edited by:

Zigui Chen, The Chinese University of Hong Kong, China

Reviewed by:

Siaw Shi Boon, The Chinese University of Hong Kong, China
Leabaneng Tawe, University of Botswana, Botswana

Copyright © 2023 Lan, Zhang, Li, He, Zhao and Yang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Fang Yang, ✉ yangfangsina@163.com

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.