Six years genotype distribution of Human Papillomavirus in Calabria Region, Southern Italy: a retrospective study

Background Although analysis of the Human papillomavirus (HPV) genotype spread in a particular area has a crucial impact on public health and prevention programmes, there is a lack of epidemiological data regarding HPV in the Calabria region of Italy. We therefore update information on HPV age/genotype distribution by retrospectively analysing a cohort of women, with and without cervical lesions, living in Calabria, who underwent HPV DNA testing; moreover, we also evaluated HPV age/genotype distribution in a subset of patients with cervical lesions. Methods Cervical scrape specimens obtained from 9590 women (age range 20–75 years) from January 2010 to December 2015 were tested for HPV DNA. Viral types were genotyped by Linear Array HPV Genotyping® test (Roche, USA) at the Clinical Microbiology Operative Unit of six hospitals located in four provinces of the Calabria region. Cervical scrape specimens were also used to perform Pap smears for cytological analysis in a subset of 405 women; cytological classification of the samples was performed according to the Bethesda classification system. Results A total of 2974 women (31%) (C.I. 95% 30.09–31.94) were found to be HPV DNA positive for at least one (57.3%) or several (42.7%) HPV genotypes. Of single genotype HPV infections, 46.5% and 36.4 % were classed as high-risk (HR, Group 1) and low-risk (LR, Group 3) respectively, while 16.9% were classed as probably/possibly carcinogenic and 0.2% undetermined risk. Stratified by age, total HPV distribution, showed the highest prevalence within the range 30–39 years (37.2%), while single genotype infection distribution displayed a peak in women from the age range 20–29 years (37.5%). The most common high-risk HPV type was HPV 16 (19.1%), followed by HPV 31 (9.1%). Conclusions We provide epidemiological data on HPV age/genotype distribution in women living in the Calabria region with or without cytological abnormalities, further to the enhancement of HPV screening/prevention programmes for the local population.


Background
Human papillomavirus (HPV) infection of the anogenital tract is the most common sexually transmitted infection; it has a broad range of clinical manifestations, ranging from subclinical and self-limited to persistent and associated with malignant progression. [1].
The mucosal HPV genotypes, grouped into the genus alpha papillomavirus, have been classified into high-risk (HR HPV) and low-risk types (LR HPV), according to their malignant potential [2]. HPV16 and HPV18 are known to be powerful carcinogens, but an International Agency for Research on Cancer (IARC) working group has classified a total of twelve, evolutionarily linked HPV genotypes (16,18,31,33,35,39,45,51,52,56,58, and 59) as presenting a high risk of carcinogenesis to humans (Group I). Based upon their phylogenetic relatedness, other types were classified as possible (Group 2B) or probable carcinogens (Group 2A) [3].
HR HPV types are a well-established cause of cervical cancer, and there is a growing body of data on the role of HR HPV in other anogenital cancers (anal, vulval, vaginal, penile), as well and head and neck cancers [4,5].
In Italy, cervical cancer is the third most common female cancer in women aged 15 to 44 years, and, according to the Ministry of Health, each year about 3,200 new cases of cervical cancer are diagnosed with 1,500 women dying from this cancer [6]. Epidemiological studies, carried out in different Italian regions [7][8][9], on the prevalence of HPV in women with normal cervix or cytological abnormalities have confirmed HPV 16 as the most frequent high-risk genotype identified in all neoplastic and non-neoplastic conditions.
However, several studies have reported heterogeneous HPV-type distributions in different geographical areas [10,11]. This depends on a complex interaction between viral molecular characteristics and host immunity, as well as the sexual behaviours and age of the target population. It is therefore essential for the distribution of type-specific HPV prevalence to be analysed in each geographical area in order for efficacious public awareness and prevention programmes to be developed.
Indeed, cervical cancer control consists of primary prevention by vaccination, and secondary prevention through screening to diagnose and treat precancerous cervical lesions. However, adherence to screening in the Calabria region is lower than the national standard (respectively 33% and 39%, steps 2010-12) [12], and therefore a regionwide public information campaign is urgently required.
With this in mind, we set out to collect epidemiological data on the prevalence of HPV infection and age/ genotype distribution in women who underwent HPV DNA testing in several centres across the region. Moreover, we also evaluated HPV age/genotype distribution in a subset of patients with cervical lesions.
For all samples DNA extraction was performed by an automated method using EasyMag (bioMérieux); eluted DNA was used for Polymerase Chain Reaction (PCR) amplification of a 450-bp fragment from the L1 HPV region, using primer PGMY09/11, by Linear Array HPV Genotyping Test (LA) (Roche, USA) [13]. An additional primer pair targets the human β-globin gene to provide a control for cell adequacy, extraction and amplification. HPV genotypes were detected by hybridization using a reverse line blot system with type-specific probes for simultaneous detection of 37 HPV genotypes. Negative and positive controls were provided with the kit and used in each test.
The study was designed as a retrospective analysis. The design was approved by the Catanzaro University Hospital Ethics Committee in compliance with the Declaration of Helsinki.

Statistical analysis
Statistical analysis was carried out on GraphPad Prism version 6.07 to establish the percentage values and their confidence intervals (CI 95%). The odds ratio (OD) value was determined using the online calculator at: https://www.medcalc.org/calc/odds_ratio.php. [14]. The OD value was used to establish the association between ASC-US, L-SIL/NILM and specific HPV groups. In all findings a p value <0.05 was considered statistically significant.
Concerning the prevalence of infections with the most representative HPV genotypes, we observed possibly carcinogenic single genotype infections accounting for 16% of overall HPV infections, and, in particular, possibly high-risk (pHR) HPV 53 exhibited a prevalence of 12.3% in multiple infections (Table 2). In single HPV infections, on the other hand, the four most common HR HPV types were HPV 16, 31, 51 and 58, with a prevalence of 16.4%, 6.8%, 4.2% and 3.7%, respectively. Stratification of the most common HR HPV types by age revealed a peak prevalence for HPV 16 and HPV 31 in the 30-39 year age group; HPV 51 was mainly detected in women aged 20-29 while HPV 58 was mostly found in the 30-39 and 40-50 age groups (Fig. 1). As a whole, these findings support the necessity for vaccination of young women in order to prevent HR HPV infections.

HPV infections and cervical lesions
In a subset of 405 women with or without cytologically demonstrated cervical lesions, analysed between January to December 2015, we also evaluated HPV status/age distribution according to cervical lesion. Accordingly, women (median age 39.  When the women's age was matched with cervical lesions and HPV status, we found that the mean age of women with negative concordant status for HPV DNA plus cytological analysis was 39.3, whereas the mean age of HPV DNA-positive women with negative cytology was 32.3 years. The mean age of women with both Pap  (Table 3). Further analysis showed that high-risk Group 1 prevalence increased in women with ASC-US (9.6%) and L-SIL (13.1%), as compared with women displaying normal cytology (4.9%), while HR HPV Group 2A/2B frequencies in women with ASC-US, L-SIL and normal cytology were 4.4%, 4.7% and 1.2% respectively. In addition, HR Group 2A/2B types showed a stronger association with L-SIL (OR 11.7, 95% CI 4.1-33.2) than HR Group 1 (OR 8.1, 95% CI 4.3-15.1) and LR Group 3 types (OR 4.7, 95% CI 2.3-9.8).

Discussion
To counter the lack of epidemiological data on HPV in Calabria (to our knowledge this is the first of such studies), we gathered data on the prevalence of HPV infection and age/genotype distribution in women attending six different hospitals in the region. This multicentre analysis revealed a total prevalence of HPV infection of 31.0 %, which, as per the available literature on the subject, varied according to age [15]. In particular, we detected the highest prevalence of single infections (37.5%) in women aged 20-29. This is in line with the results of several studies that have shown a greater prevalence of HPV in younger women, followed by a reduction of HPV infections with increasing age [16,17]. That being said, we also found a peak of total, as well a multiple, infections in the 30-39-year age group (37.2%), which may reflect local sexual behavioural trends (e.g. number of partners).
Regarding the different HR HPV genotypes in the tested population, we detected several that were most represented in both single and multiple infections, namely HPV 16, 31, 51 and 58. HPV 18, a type strongly associated with cervical cancer, was detected in 3.6% of women. Our results are similar to those recently reported by Veldhuijzen et al., in which two European female populations were tested in HPV screening trials [18]. Comparable results were also obtained in Sardinia, a major Italian island, during a study on the female population exhibiting cervical cytological abnormalities [19]. Also, our data are consistent with those previously reported on HPV types prevalence and distribution in Italian population [8,9].
The high prevalence of HPV 16 in both single and multiple infections, mainly in 30-39 age group, confirms the widespread distribution of this virus. It is also important to note that HPV 18, the second most common high-risk HPV types in cervical cancer, was rather high in the 30-39 age group, even though its prevalence was very low (3.6%) across the overall population [20][21][22]. Another important finding was the distribution of HPV types not covered by current vaccines. Among these, HPV 51, 53, 73 and 66 were observed to have a high prevalence among the women analysed in our sample. As a whole, these observations may be of interest to those working on vaccination and local HPV prevention schemes.
In particular, the coverage of the HPV vaccine schedule in the Calabria region is very low (ranging from 46.36% to 70.88 % of girls born in the period from 1997 to 2003) [23]; increasing HPV-vaccination coverage, both across the local population and in terms of HPV targeting, therefore represents a major public health issue. In Table 3 HPV distribution across different cervical cytological statuses. Prevalence of HPV positivity in single and multiple infections, according to cytological cervical lesion this regard, our findings underscore the importance of keeping up-to-date information on HPV distribution, both to monitor vaccine efficacy and the spread of the disease.

Conclusions
Our epidemiological data on the HPV age/genotype distribution in the Calabria region of Italy indicate the need to improve the local HPV screening plan, whose efficacy should be closely monitored in order to maximise efficacy while minimising the cost to the local health service.