In the present cross-sectional study, the prevalence of HPV genotypes in the biological samples of a population from the southwest of Iran (Fars province) was assessed. The results showed a positivity rate of 29.9% in which the high-risk genotypes comprised the most genotypes in the studied population. To determine the prevalence of HPV genotypes, several studies were performed in different parts of Iran. Of them, a previous study from the northeast of Iran (2013 to 2018) reported a positivity rate of 48.4% for HPV genotypes [30]. Another study reported a prevalence of 49.5% on the samples which were from different provinces of Iran, from April 2011 to April 2016 [21]. In addition, a recent study from the capital city of Iran, Tehran, reported a positivity rate of 53% on the samples collected from 2017 to 2021 [28]. Also, another study on 12076 biological samples reported a positivity rate of 38.68% which was carried out from 2016 to 2018 [5]. In our neighboring country Iraq, with whom we share a lot of religious ties, the positive rate was reported as 20.9% on samples collected between January 2018 and September 2020. [9]. Another study reported an HPV-positive rate of 17.96% in Iraqi women aged between 15 to 50 years [44]. The inconsistency between the positive rates in our results and the others may be due to the application of different methods (PCR, Hybridization) and variations in the commercial kits capable of detecting different numbers of genotypes. In addition, community awareness of HPV infections and its possible complications as well as vaccine availability and more intensive preventive behaviors may also be the reason for the decline in HPV transmission.
Distributions of HPV genotypes revealed that the samples with positive high-risk genotypes were more frequent than positive low-risk genotypes. Of the positive samples, the most frequent genotypes were HPV6 and HPV16, respectively. In concordance with our study, genotypes 6 and 16 were also the most low- and high-risk genotypes in the other studies reported from Iran [16, 27, 30].
Geographically, comparing the distributions of HPV in our result with different continents and countries, the highest prevalence has been reported in the Asian, African, and European regions, respectively. Also, the prevalence of HPV infection in the European areas revealed a lower rate of positivity in the developed country in comparison with developing countries [3, 32, 38]. Like ours, a global study reported that high-risk genotypes are more frequent than low-risk ones [6, 12, 17]. Furthermore, in agreement with our results, a study from China, where we have a lot of travel and business dealings, reported a higher frequency of high-risk genotypes (19.2%) compared with low-risk ones (6.4%) [40].
Remarkably, we observed that genotypes 53, 66, 51, and 39 were more frequent high-risk genotypes after HPV16, respectively. However, other Iranian studies reported different high-risk genotypes following HPV16 such as genotypes 52 [16], 66 [14], 51 [30], and 68 [31]. Given that the available 9-valent HPV vaccine is protective against high-risk genotypes 16, 52, 31, 18, 58, 45, and 33 as well as low-risk genotypes 6 and 11, and based on the other Iranian study and ours, it seems that some more frequent HPV high-risk genotypes such as 53, 66, 51 and 39 among Iranian population were not covered by the currently available vaccine. In addition, in the present study, we observed that genotype 44 which is the second most frequent low-risk genotype is not included in the vaccine. Exploring the distribution of HPV genotypes in other countries revealed that in China, genotypes 52, 58, 16, 51, and 56 were more frequent, respectively [45]. Another study from Beijing, China, reported that genotypes 52, 58, 16, 39, and 51 were the most common high-risk genotypes [39]. A study from Weifang reported that genotypes 16, 52, 58, 53, and 68 were more frequent. A report from Maputo city, Mozambique, showed that genotypes 52, 35, 16, 53, 58, and 51 were the high-risk genotypes [24]. Researchers from Shanghai, China, reported that high-risk genotypes 52, 58, 16, 53, 51, and 81 were more prevalent [15] and A study from the Metropolitan Area of Naples, Italy, reported that genotypes 16, 31, 18, and 51 were the most frequent types [15]. Accordingly, it seems that the other more frequent high-risk genotypes with worldwide distribution should be added to the new vaccine formulation to significantly decline the outcomes of HPV infections in the world especially in our country, Iran.
Investigating the association of HPV genotypes with age range showed that a higher frequency of HPV infection was observed in age groups less than 30 years old, while, in the individuals with age > 60 years old the prevalence rate of HPV infection statistically reduced. In agreement with our study, another study from Iran and other countries reported the same association between HPV prevalence and age ranges. Almost in all studies, the highest frequency was observed in the ages < 30 years old [8, 16, 26, 36]. CDC recommends two schedules for HPV vaccination based on the age of vaccinated individuals. A two-dose series (0, 6–12 months) for most persons who initiate vaccination at ages 9 through 14 years and a three-dose series (0, 1–2, 6 months) for persons who initiate vaccination at ages 15 through 45 years, and also for immunocompromised patients [1]. So, it seems that based on the CDC recommendation and the prevalence rate of HPV infection, HPV vaccination for individuals less than 30 years old could be most effective in preventing the high-risk age group from future cancer-associated HPV.
A previous study reported that in 20–50% of HPV-infected women, especially young women, co-infection with more than one HPV type is common. Infection with multiple HPV types is often considered a risk factor leading to the development of cervical cancer [10]. In agreement with previous studies, we also observed 35% co-infection in our study groups. Among them, co-infection with 10 genotypes were also detected in two samples and co-infection with very high-risk genotypes 16 and 18 were determined in one percent of positive samples. Given the frequency of co-infection in our and the previous study and according to that more co-infection was observed between genotype 6, the most frequent genotype, and the high-risk genotypes 16, 66, and 53, respectively, it seems that vaccination with a vaccine containing more than 2 (Cervarix) and 4 (Gardasil, 4vHPV) genotypes like 9-valent is essential for the future vaccination strategy. In a Chinese large study, co-infection of genotypes 16 and 18 was similar to ours and they also reported that co-infection of HPV16 with other high-risk genotypes was 23.54% [42]. Among the multi-genotype HPV infection and disease outcomes, it has been reported that the risk of high-grade squamous epithelial lesions increases with the number of HPV types (Odds ratio for single, 2 to 3 and 4 to 6 HPV types were 41.5, 91.7 and 424, respectively), compared with HPV negative samples [35]. A recent study reported that TLR4 rs10759931 is protective against multiple high-risk HPV infections and in contrast, haplotype ACAC (at the loci rs7873784, rs4986791, rs4986790 and rs4986791) in the TLR-4 as pathogen recognition receptor (PRR) on the innate immune cells increases the risk of infection with multiple high-risk genotypes [25]. So it seems that gene variants associated with immune responses can act as predisposing factors in the multiple HPV genotype infection. Our study limitations were that we selected the study population from individuals without and with symptoms referring to a subspecialty clinic for routine checkups which imply that they may be a high-risk group instead of the normal population that can affect the prevalence of HPV in the community. In addition, there was no information about pop smear results and clinical manifestations in the HPV-infected individuals to compare with HPV-negative ones. Furthermore, the situation of cervical cancer in the patients was not clear. As the positive rate of HPV infection was higher (29.9%) in our study compared with the results obtained from other population-based studies (9.73% in Ethiopia [33], 13.22% in China [19], and 13.6% in Canary Island [2], to determine the HPV distributions reflecting HPV prevalence, future studies with higher sample size collected from healthy individuals in the different parts of Iran is warranted.