Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Factors associated with high-risk human papillomavirus infection and high-grade cervical neoplasia: A population-based study in Paraguay

  • Elena Kasamatsu,

    Roles Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Supervision, Writing – original draft, Writing – review & editing

    Affiliation Instituto de Investigaciones en Ciencias de la Salud, Universidad Nacional de Asunción (IICS-UNA), San Lorenzo, Paraguay

  • María Isabel Rodríguez Riveros,

    Roles Data curation, Formal analysis, Investigation, Methodology, Supervision

    Affiliation Instituto de Investigaciones en Ciencias de la Salud, Universidad Nacional de Asunción (IICS-UNA), San Lorenzo, Paraguay

  • Ana María Soilan,

    Roles Investigation, Methodology

    Affiliation Hospital Materno Infantil de San Lorenzo, Ministerio de Salud Pública y Bienestar Social, San Lorenzo, Paraguay

  • Marina Ortega,

    Roles Investigation, Methodology

    Affiliation Hospital Nacional, Ministerio de Salud Pública y Bienestar Social, Itauguá, Paraguay

  • Pamela Mongelós,

    Roles Investigation, Methodology

    Affiliation Instituto de Investigaciones en Ciencias de la Salud, Universidad Nacional de Asunción (IICS-UNA), San Lorenzo, Paraguay

  • Malvina Páez,

    Roles Investigation, Methodology

    Affiliation Instituto de Investigaciones en Ciencias de la Salud, Universidad Nacional de Asunción (IICS-UNA), San Lorenzo, Paraguay

  • Amalia Castro,

    Roles Investigation, Methodology

    Affiliation Instituto de Investigaciones en Ciencias de la Salud, Universidad Nacional de Asunción (IICS-UNA), San Lorenzo, Paraguay

  • Carmen Cristaldo,

    Roles Investigation, Methodology

    Affiliation Hospital Nacional, Ministerio de Salud Pública y Bienestar Social, Itauguá, Paraguay

  • Fátima Romina Báez,

    Roles Investigation, Methodology

    Affiliation Instituto de Investigaciones en Ciencias de la Salud, Universidad Nacional de Asunción (IICS-UNA), San Lorenzo, Paraguay

  • Claudia Carolina Centurión,

    Roles Investigation, Methodology

    Affiliation Instituto de Investigaciones en Ciencias de la Salud, Universidad Nacional de Asunción (IICS-UNA), San Lorenzo, Paraguay

  • Jaime Vester,

    Roles Formal analysis, Methodology

    Affiliation Instituto de Investigaciones en Ciencias de la Salud, Universidad Nacional de Asunción (IICS-UNA), San Lorenzo, Paraguay

  • Hernán Barrios,

    Roles Investigation, Methodology

    Affiliation Instituto de Investigaciones en Ciencias de la Salud, Universidad Nacional de Asunción (IICS-UNA), San Lorenzo, Paraguay

  • Griselda Villalba,

    Roles Investigation, Methodology

    Affiliation Hospital Materno Infantil de San Lorenzo, Ministerio de Salud Pública y Bienestar Social, San Lorenzo, Paraguay

  • María Luisa Amarilla,

    Roles Investigation, Methodology

    Affiliation Hospital Nacional, Ministerio de Salud Pública y Bienestar Social, Itauguá, Paraguay

  • Graciela Giménez,

    Roles Investigation, Methodology

    Affiliation Instituto de Investigaciones en Ciencias de la Salud, Universidad Nacional de Asunción (IICS-UNA), San Lorenzo, Paraguay

  • Elodie Caubere,

    Roles Investigation, Methodology

    Affiliation International Agency for Research on Cancer, IARC/WHO, Lyon, France

  • María de la Luz Hernández,

    Roles Data curation, Methodology

    Affiliation International Agency for Research on Cancer, IARC/WHO, Lyon, France

  • Armando Baena,

    Roles Formal analysis, Writing – review & editing

    Affiliation International Agency for Research on Cancer, IARC/WHO, Lyon, France

  • Maribel Almonte,

    Roles Conceptualization, Data curation, Formal analysis, Methodology, Validation, Writing – review & editing

    Affiliation International Agency for Research on Cancer, IARC/WHO, Lyon, France

  • Rolando Herrero,

    Roles Conceptualization, Formal analysis, Funding acquisition, Methodology, Validation, Writing – review & editing

    Affiliation International Agency for Research on Cancer, IARC/WHO, Lyon, France

  • Laura Patricia Mendoza ,

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Supervision, Validation, Writing – original draft, Writing – review & editing

    lauramendozatorres@gmail.com

    Affiliation Instituto de Investigaciones en Ciencias de la Salud, Universidad Nacional de Asunción (IICS-UNA), San Lorenzo, Paraguay

  •  [ ... ],
  • for the ESTAMPA Paraguay Center study group

    Membership of the ESTAMPA Paraguay Center study group is provided in the Acknowledgments.

  • [ view all ]
  • [ view less ]

Abstract

Background

Cervical cancer (CC) is one of the leading causes of cancer mortality among women from Paraguay, with high incidence and mortality rates (31.2 and 16 per 100 000 women, respectively). Although the risk factors associated with high-risk human papillomavirus (hrHPV) infection and preneoplastic cervical lesions are widely studied, population-based characteristics of particular settings may influence the feasibility of HPV-based CC screening implementation. This study aimed to explore factors associated with hrHPV infection and high-grade cervical neoplasia in hrHPV-positive (hrHPV+) women from Paraguay.

Methods

A total of 5677 women aged 30–64 years from the Central Department of Paraguay were screened with HPV test (Hybrid Capture 2) and Pap smear. Sociodemographic and risk factor interviews were conducted. hrHPV+ women were referred to colposcopy and women with an abnormal colposcopy had a biopsy taken. The outcomes recorded were the hrHPV status and the presence of high-grade cervical intraepithelial neoplasia or worse (CIN2+) among hrHPV+ women. Associations were investigated using multivariate logistic regressions.

Results

hrHPV prevalence was 13.8% (95%CI 13.0–14.8). This value decreased with the age of women (p-trend<0.001) and increased with the lifetime number of sexual partners (p-trend<0.001) and number of previous female partners of their current male partner if women had had one lifetime sexual partner (p-trend<0.001), increasing from 3.06 (95%CI 0.073–20.9) if partners had had one previous female partner to 9.19 (95%CI 2.36–61.1) if they had had eight or more. In hrHPV+ women, CIN2+ prevalence was 10.7% (95%CI 8.58–13.2) and increased with time since the last Pap smear (p-trend<0.001) and with the increasing number of pregnancies (p-trend = 0.05).

Conclusion

In these settings, the sexual behavior of women and their male partners is associated with hrHPV infection. In hrHPV+ women, underscreening practices and multiple pregnancies are associated with CIN2+. This knowledge can contribute to public health policies for CC prevention and control in Paraguay.

Introduction

Cervical cancer (CC) ranks globally as the fourth most frequently diagnosed cancer and the fourth leading cause of cancer death in women, with an estimated 570 000 new cases and 311 000 deaths in 2018 [1] and more than 85% of the cases occurring in developing countries. In Paraguay, the burden of CC is a significant health problem with high incidence and mortality rates of 31.2 and 16 per 100 000 women, respectively [1].

It has been recognized that the main cause of CC and precancerous lesions is high-risk human papillomavirus (hrHPV) [2,3]. A worldwide study has shown that the most common HPV types associated with CC are HPVs 16, 18, 31 33, 35, 45, 52 and 58 [4], which account for 91% of the cases. In Paraguay, the most prevalent types in CC are HPVs 16, 18 and 45 (78.5%) [5,6], whereas the most prevalent ones in high-grade squamous intraepithelial lesions diagnosed by histology or cytology are HPVs 16, 18 and 33 [7]. HPV is one of the most common sexually transmitted infections. Most of the infections are transient and clear within two years without any treatment. Persistent infection with oncogenic hrHPV is considered as a necessary cause but not sufficient to develop cervical precancerous lesions and cancer [8,9]. Other cofactors such as multiparity, long-term hormonal contraceptive use, tobacco smoking and HIV infections, influence the risk of progression to high-grade squamous intraepithelial lesion or cancer in hrHPV infected women [1012].

In low-and middle-income countries, cytology-based programs are difficult to implement and are characterized by low coverage and unsatisfactory follow-up. It has been demonstrated that HPV testing is substantially more sensitive than cytology for detection of high-grade cervical lesions and, although less specific, it has large potential to reduce CC rates [1319]; thus, this method has been recommended by WHO as the primary cervical screening method to implement worldwide [20]. The national cytological screening program established in Paraguay more than 40 years ago has not been successful in controlling CC because of the low sensitivity of cytology, unsatisfactory follow-up of women with cytological abnormalities and the limited coverage reported in Paraguayan women (less than 20%) [21]. Currently the use of HPV testing is limited to private patients and to one of the reference public gynecology health centers but is not yet included in the national screening program for cervical cancer prevention and control.

This study is part of a large international project entitled “Multicentric study of cervical cancer screening and triage with human papillomavirus testing—ESTAMPA study” (NCT01881659) coordinated by the Prevention and Implementation Group of the International Agency for Research on Cancer/World Health Organization (IARC/WHO) and being conducted in nine countries of Latin America (12 centers), including Paraguay, since 2014. The aim of the ESTAMPA study is to investigate the performance of emerging CC screening and triage techniques in women aged 30–64 years old and to evaluate the feasibility of different approaches for the implementation of organized HPV-based screening programs. The target number of women to be recruited in the different research centers across Latin America is around 50 000.

In Paraguay, HPV types associated with CC and preneoplastic lesions have been previously studied but there are no population-based data or information on HPV cofactors. Therefore, within the framework of the ESTAMPA study, the aim of this study was to explore the prevalence and factors associated with hrHPV infection and CIN2+ in hrHPV+ women in the Central Department of Paraguay. The information collected may influence the feasibility of HPV-based cervical cancer screening implementation and could contribute to public health policies for CC prevention and control in Paraguay.

Materials and methods

Study design

This was a cross-sectional screening study, in which hrHPV testing was carried out in 5677 women aged 30–64 years old recruited in two districts of the Central Department of Paraguay during 2014–2018. All hrHPV+ women were referred to colposcopy, histological diagnosis and treatment as needed. The first outcome recorded was the hrHPV status of each woman, whereas the second one was the presence of high-grade cervical intraepithelial neoplasia or worse (CIN2+), including CIN2, CIN3 and CC among hrHPV+ women. To assure high-quality data and comparability between districts, standardized study forms and procedures, including monitoring and quality control were applied, based on the ESTAMPA protocol. Additionally, study investigators and coordinators received specific training in good clinical practice, field work, pathology, colposcopy, molecular biology and statistics.

The Central Department of Paraguay is made up of 19 districts and includes the city of Asunción, which is the capital of the country. According to data from the Dirección General de Estadística, Encuestas y Censos (DGEEC), in 2018 approximately 426 739 women aged 30 to 64 years were living in the Central Department. [22]

The two districts mentioned above were selected for convenience from a well-defined catchment area of the Central Department with a baseline population of approximately 92 804 women aged 30 to 64 years. These districts included women of medium/low socioeconomic level and were close to primary health care centers as well as close to a referral hospital to facilitate the referral of hrHPV+ women to colposcopy visits and for treatment in case of CIN2+ (Fig 1). One of the selected districts was San Lorenzo, distant 14 km away from Asunción, the capital city of Paraguay where, with the collaboration of the Expanded Program on Immunization (EPI) the study group organized the census described below. In San Lorenzo, the referral Hospital Materno-Infantil de San Lorenzo, dependent on the Ministry of Public Health and Social Welfare of Paraguay (MSPyBS) collaborated with this study. The other district was Itauguá, 30km away from Asunción, including the urbanization developed by the National Housing Council of Paraguay (CONAVI) under the coverage of the referral Hospital Nacional de Itauguá, MSPyBS.

thumbnail
Fig 1. Map showing the political division of Paraguay, the location of Central Department and Itauguá and San Lorenzo Distrits.

https://doi.org/10.1371/journal.pone.0218016.g001

Information booklets were prepared and a door-to-door census was carried out, at which women or household members were informed about cervical screening, HPV testing, and the study. Potential eligible women aged 30 to 64 years old were invited to attend screening at the local clinic mainly by telephone calls. Women attending the clinic were then informed of the study and asked to sign an informed consent. A total of 691 out of 10 773 potentially eligible women were ineligible due to exclusion criteria such as hysterectomy, pregnancy/postpartum, not available for pelvic examination, pre-cancer treatment, no previous sexual intercourse, history of CC or serious pre-existing medical conditions (e.g. advanced cancer, terminal renal failure). Out of 10 082 potentially eligible women that fulfilled the inclusion criteria, 4984 accepted to participate in the study. Other 693 eligible women from the Central Department were also enrolled by invitation (poster, talks) at the Universidad Nacional de Asunción, Paraguay. Thus, a total of 5677 women were included in the study (Fig 2).

thumbnail
Fig 2. Flow diagram of participants.

*Added participants enrolled by invitations (posters, talks) at the Universidad Nacional de Asunción, Paraguay. †Women included in the analysis for the first outcome: hrHPV status and associated risk factors. ‡Women included in the analysis for the second outcome: presence of CIN2+ and associated risk factors.

https://doi.org/10.1371/journal.pone.0218016.g002

Eligible women were invited to participate in the study and scheduled to visit the nearest health center, through repeated phone calls. At the health center, after initial information and orientation about the study and the signature of the informed consent, trained staff administered a socio-demographic and risk factor questionnaires (Supporting information S1S6 Files) to obtain data on years of education, year of last Pap smear, age at first sexual intercourse, lifetime number of sexual partners, number of previous female partners of their current male partner, number of pregnancies, use of contraceptive methods, cigarette use and years and lifetime number of cigarettes smoked. Then, trained nurses performed a pelvic exam and cervical sample collection for Pap smear and detection of hrHPV. The samples collected were labeled and sent to the laboratory of the Instituto de Investigaciones en Ciencias de la Salud, Universidad Nacional de Asunción (IICS-UNA) for processing. hrHPV+ women were adequately informed and invited to have a colposcopy done. hrHPV-negative (hrHPV-) women returned to regular screening with Pap smear, offered free of charge by the MSPyBS.

hrHPV detection

Cervical samples for hrHPV testing were collected in PreservCyt medium and 4 ml of each residual PreservCyt sample was centrifuged and resuspended in 150 μl of a mixture of specimen transport medium and denaturation buffer provided in the kit. A 75 μl aliquot of each resuspended sample was used for each Hybrid Capture 2 assay, a sandwich capture molecular hybridization assay that uses chemiluminescence detection. All assays were performed according to the manufacturer’s protocol, with the results reported as a ratio of relative light units (RLU/CO). Samples with an RLU/CO ratio >1.0 were considered positive for hrHPV.

Colposcopy, histopathological study and clinical management

All hrHPV+ women were referred to colposcopy for diagnosis and treatment by two colposcopists trained in the study protocol. In women with abnormal colposcopic findings, 2–3 biopsies from the aceto-white lesions observed were taken, fixed in 10% neutral buffered formalin and transported to a pathology laboratory for histological procedures. Histologically, cervical lesions were classified as normal, CIN1, CIN2, CIN3, and CC. The final histopathological diagnosis was based on the highest grade of the lesion observed in biopsy or in the histological specimen of the large loop excision of the transformation zone (LLETZ) of each woman. All histology samples were read by the local pathologist experienced in CC pathology trained in the study protocol. Women with CIN2+ were treated with LLETZ and those with CC were referred to a tertiary care center. Women with negative colposcopy or histology less than CIN2 (<CIN2) were invited to a follow-up hrHPV screening at 18 months, whereas those with positive colposcopy or CIN2+ had a final colposcopy. Women with histological CIN2+ were treated adequately and hrHPV- women or with <CIN2 are given clear follow-up indications to continue with regular screening within the local health care system.

Statistical analysis

The main outcomes were the hrHPV infection status of each woman and the CIN2+ lesions present in hrHPV+ women. Bivariate analyses between sociodemographic characteristics and known risk factors (independent variables) for hrHPV infection and CIN2+ were conducted using Chi-square or Mann-Whitney tests for categorical or numerical variables, respectively. Age-adjusted and multivariate-adjusted odds ratios (OR) with their corresponding 95% confidence intervals (CI) were estimated from logistic regressions to explore associations between independent variables and outcomes. Independent variables regarding well-established risk factors for both hrHPV infection (sexual behavior variables) [23,24] and CIN2+ (number of pregnancies, cigarette consumption and contraceptive use) [25,26] were selected for the corresponding multivariate analyses adjusted by age, education level and screening history (time since the last cytology). Herein we included age of first intercourse, use of condom and lifetime number of sexual partners for the hrHPV multivariate analysis, and number of pregnancies, years of cigarette use and years of hormonal contraceptive use for the CIN2+ multivariate analysis. Additionally, we evaluated whether the sexual history of the current partner of women with only one lifetime sexual partner was associated with hrHPV infection. Therefore, an additional analysis among women with only one lifetime sexual partner was conducted. Analyses were performed in R (R Core Team (2018). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. URL https://www.R-project.org/) and all tests were two-sided with a significance level of 0.05.

Ethical considerations

The protocol was approved by the Ethics Committee of the International Agency for Research on Cancer, World Health Organization (IARC/WHO), Pan American Health Organization (PAHO) and by IICS/UNA.

Results

Characteristics of the study population

A total of 5677 of the 10 082 30-64-year-old eligible women and 693 women invited by the UNA were enrolled. Among them, 4354 women were from the San Lorenzo district and 630 women from the Itauguá district. The distribution of socio-demographic, reproductive and sexual behavioral characteristics of the study population is presented in Table 1. The mean/SD age was 44.5/9.25. A total of 5496 (96.9%) reported having had Pap smears and almost half of these women (48.1%) had the last Pap smear in the previous 2 years. However, 1375 (24.2%) had the last Pap smear >5 years or never before the interview; 4246 (74.8%) had high school education or less, 2978 (52.5%) had their first intercourse at 17–20 years, 3925 (69.3%) had ≥2 lifetime sexual partners 438 (25.2%) out of the 1739 women with one lifetime sexual partner answered about the number of previous female partners of their current male partner, 2428 (42.8%) had 2–3 pregnancies and 4844 (85.3%) never smoked. Regarding contraceptive methods, 3145 (55.5%) never used condoms and 3850 (67.8%) used hormonal contraceptives, among which 2643 (46.6%) reported using them for less than 10 years.

thumbnail
Table 1. Characteristics and prevalence of hrHPV in 30-69-year-old women from the Central Department Paraguay.

https://doi.org/10.1371/journal.pone.0218016.t001

Prevalence of hrHPV infection and CIN2+ in hrHPV+ women

The overall hrHPV prevalence was 13.8% (95%CI 13.0–14.8). The highest prevalence of hrHPV infection was observed among women aged 30–34 years (21.12%, 95%CI 18.5–24.0) and decreased with increasing age to 11.8% (95%CI 9.49–14.6) in 55-59-year-old women (Table 1).

Among the 694 hrHPV+ women with histopathological diagnosis, 74 (10.7%, 95%CI 8.58–13.2) had CIN2+. The mean/SD age of CIN2+ women was 40.2/7.42, with higher prevalences in women aged 35–39 years (14.1%, 95%CI 9.31–20.8), 40–44 years (13.6%, 95%CI 8.27–21.5) and 50–54 years (13.7%, 95%CI 7.61–23.4) (Table 2). The histopathological diagnosis of CIN2+ showed 8 cases of CIN2, 55 CIN 3 and 11 CC.

thumbnail
Table 2. Characteristics and prevalence of CIN2+ in 30-69-year-old hrHPV+ women with histopathological diagnosis from the Central Department, Paraguay.

https://doi.org/10.1371/journal.pone.0218016.t002

Risk factors for hrHPV infection and CIN2+ among hrHPV+ women

Risk factors for hrHPV infection were investigated in 5677 women, whereas risk factors for CIN2+ were analyzed in 694 hrHPV+ women (Fig 2). Table 3 shows the risk factors associated with hrHPV infection. The multivariable adjusted analysis showed that the risk of hrHPV infection decreased with the age of the woman (p-trend<0.001) and increased with the number of lifetime sexual partners, compared ≥6 partners with one partner (OR = 1.81, 95%CI 1.37–2.38, p-trend<0.001). Among the 1739 women with one lifetime partner the risk of hrHPV infection increased with the number of previous female partners of their current male partner (p-trend<0.001).

thumbnail
Table 3. Odds ratios of hrHPV associated with sociodemographic characteristics and selected risk factors among 30-69-year-old women from the Central Department, Paraguay.

https://doi.org/10.1371/journal.pone.0218016.t003

Comparison of women with one lifetime sexual partner whose partners had not had previous female partners showed that the adjusted OR of having hrHPV increased from 3.06 (95%CI 0.73–20.9) if partners had had one previous female partner to 9.19 (95%CI 2.36–61.1) if they had had ≥8 (Table 4).

thumbnail
Table 4. Odds ratios of hrHPV infection associated with sociodemographic characteristics and selected risk factors among 30-69-year-old women with one lifetime sexual partner from the Central Department, Paraguay.

https://doi.org/10.1371/journal.pone.0218016.t004

The time of the last Pap smear, education level, age of first intercourse, having had different sexual partners during the last year, number of pregnancies, cigarette use, condom use and hormonal contraceptive use were not significantly associated with hrHPV infection.

Among hrHPV+ women, the risk of CIN2+ increased with the time of the last Pap smear (p-trend<0.001). As compared with women who had had a Pap smear in the last two years, those having had a Pap smear within 2–5 years were at increased risk (OR = 1.25, 95%CI 0.62–2.44), with a further increase in risk of CIN2+ for women with more than 5 years since the last Pap smear (OR = 3.00, 95%CI 1.67–5.47). In addition, the risk of CIN2+ was higher among women with ≥4 pregnancies (OR = 2.72, 95%CI 1.11–7.45) than among those with 0–1 pregnancy (p-trend = 0.05). The age, education level, years of hormonal contraceptive use and years of cigarette use were not significantly associated with the risk of CIN2+ (Table 5).

thumbnail
Table 5. Odds ratios of CIN2+ associated with sociodemographic characteristics and selected risk factors among 30-69-year-old hrHPV+ women from the Central Department, Paraguay.

https://doi.org/10.1371/journal.pone.0218016.t005

Discussion

This is the first population-based study of hrHPV infection and CC precursors using hrHPV testing in Paraguay, which assesses the prevalence and risk factors of hrHPV infection and CIN2+ in hrHPV+ women living in two districts of the Central Department. The results showed that the prevalence of hrHPV infection (13.8%) was consistent with that observed in Paraguayan women without cytological abnormalities in five indigenous communities of Paraguayan Chaco (16.1%) and those coming from reference health centers in the Central Department, of Paraguay (13.5%). [7,27]. The higher hrHPV prevalence observed in indigenous women could be due to the fact that the median age was lower (30) than that of women enrolled in the present study (44.5).

The prevalence of HPV differs worldwide, with higher prevalence values in Sub-Saharan African regions (24%), Latin America and the Caribbean (16.1%), Eastern Europe (14.2%) and Southeastern Asia (14%) [28]. The hrHPV prevalence found in the present study in Paraguayan women is comparable to that found in more recent results from Latin America and the Caribbean centers participating in the ESTAMPA study, using comparable methodologies (Almonte, personal communication).

This study showed decreasing hrHPV prevalence with increasing age, consistent with previous studies [2831]. A non-significant small second peak was observed in postmenopausal women aged 60–64. This second peak at perimenopausal or postmenopausal ages has been previously described in Latin American women and could be related to a cohort effect due to higher HPV exposure of older women when they were young, to the reactivation of latent HPV infection by decreased immune response or hormonal changes or to new HPV infection with an increase of low risk HPV [29, 3235].

In the present study, we also explored the risk factors for hrHPV infection and found a strong association with the lifetime number of sexual partners. In Latin America and the Caribbean, the acquisition of HPV infection is associated with an increased number of lifetime sexual partners as well as with the number of recent sexual partners [5,16,29,33,36]. A pooled analysis of the IARC HPV surveys of 11 337 women from Vietnam, Thailand, Korea, Nigeria, Spain, Argentina, Chile, Mexico and Colombia, confirmed that the lifetime number of sexual partners was the most important risk factor for HPV infection and that the heterogeneity in the ORs for HPV positivity associated with multiple sexual partners is probably due to the fact that some women underreport their lifetime number of sexual partners, which is more severe in the traditional Asian Society than in Europe, Africa or Latin America [23]. In addition, a clear increase in risk was observed with the number of previous female partners of the current male partners of women with one lifetime sexual partner. The strongest association was observed for partners with eight or more additional female partners (OR = 9.19, 95%CI 2.36–61.1; p-trend<0.001) compared to women who reported that their partners had had only one lifetime sexual partner. From our point of view, considering the cultural characteristics of Paraguayan women, including their language (guaraní language), well trained study field workers made all efforts and collect as adequate as possible information about the number of lifetime sexual partners, different sexual partners during the last year and number of previous female partners of their husband or current male partner in women with one lifetime sexual partner. However, we can not exclude potential reporting bias, even if the woman assures to know the sexual behavior of her partner.

In hrHPV+ women, we observed a strong association between CIN2+ and the time to their last Pap smear, particularly in those who had had their last Pap smear more than 5 years before the interview (OR = 3.00, 95% CI 1.67–5.47; p-trend<0.001). A previous study on CC in Paraguay, showed that never having had a Pap smear was one of the main risk factors for CC (OR = 26.7, 95%CI 2.4–296.9) [5]. This finding is important for CC prevention and control program in Paraguay, in which the cytology-based national screening program is opportunistic, with limited coverage (20%) [21], although Pap testing has been free of charge for many years. Also, this knowledge can contribute to supporting the incorporation of hrHPV testing, that allows extension of screening intervals as the primary screening method in our country.

Our results also showed that having had a large number of pregnancies was associated with higher risk of CIN2+, with an almost three-fold increase for women with more than four pregnancies compared with those with 0–1 pregnancies. Several studies have reported the association of high parity with the risk of CIN2+, mainly in younger women [10,11,3739], probably because of the prolonged period of exposure of the transformation zone to hrHPV and/or other cofactors in highly parous women [38]. The increased levels of estrogen and progesterone during pregnancy are probably responsible for the alterations in the transformation zone, potentially inducing a reduced immune response to HPV infection and influencing risk of persistence or progression [39]. A large collaborative epidemiological study on CC and hormonal contraceptives concluded that the risk of CC in current hormonal contraceptive users increases with the increasing duration of use [40]. Another study has suggested that neither high parity nor long-term use of hormonal contraceptives are associated with HPV prevalence, but may be involved in the transition from HPV infection to neoplastic cervical lesions [25].

Two large multicentric studies coordinated by the IARC have demonstrated that the education level is inversely and consistently associated with CC risk, but not with HPV prevalence [41]. It has also been found that the association between low socioeconomic status and level of education (especially no schooling) increases the risk for precancerous cervical lesions and invasive cancer in women with hrHPV [37]. In the present study, the socioesconomic status was not analyzed and the education level was not associated with CIN2+.

Regarding the knowledge, attitudes and practice related to HPV and CC, in a study conducted among women aged 30 years or more of a marginal riverside neighborhood of Asunción, Paraguay, only 10% of the women answered that they knew about HPV and 27% responded that Pap testing was associated with CC prevention [42]. Education, information and orientation about HPV infection, risk factors and also on the methods for the prevention of CC are necessary, to improve the strategy of prevention and control of CC. Furthermore, based on the association of the sexual behavior of the male partner observed in this study, we consider it appropriate to promote the participation of husbands or male partners of women in the information, orientation and education programs about prevention of HPV infection and the risk of hrHPVassociated genital cancer.

One of the strengths of our study was that the analysis of risk factors for CIN2+ was restricted to hrHPV+ women, a fact that allowed eliminating potential confounding by factors associated with the acquisition of infection. Another strength is that a standard protocols for data and specimen collection, laboratory testing and quality control, as well as a clinical management of HPV positive women that was validated by an international group of experts. Future analyses will provide additional detail, including HPV typing, multiple HPV infections and other sexually transmitted co-infecions, socio-economic factors, environmental and cultural factors, access to health centers and others.

In conclusion, this study provides the first population-based data on the risk factors for hrHPV infection and preneoplastic cervical lesions in an unvaccinated and not previously screened with HPV testing Paraguayan population. The sexual behavior of male partners was found to be associated with hrHPV infection, whereas underscreening practice and high number of pregnancies were found to be associated with CIN2+. Greater efforts to promote CC screening in our country are needed, particularly including educational campaigns on HPV infection and CC prevention targeting not only women but the whole population.

Supporting information

S3 File. Socio-demographic form in english.

https://doi.org/10.1371/journal.pone.0218016.s003

(PDF)

S4 File. Socio-demographic form in spanish.

https://doi.org/10.1371/journal.pone.0218016.s004

(PDF)

S5 File. Risk factors questionnaire form in english.

https://doi.org/10.1371/journal.pone.0218016.s005

(PDF)

S6 File. Risk factors questionnaire form in spanish.

https://doi.org/10.1371/journal.pone.0218016.s006

(PDF)

Acknowledgments

The authors thank the International Agency for Research on Cancer (IARC), the Ministry of Health and Science of Paraguay (Hospital Nacional de Itaugua, Hospital Materno Infantil de San Lorenzo and Primary Health Care Units), and Instituto de Investigaciones en Ciencias de la Salud for all their support.

ESTAMPA Paraguay Center study group are listed below. Laura Patricia Mendoza, Elena Kasamatsu, María Isabel Rodríguez Riveros, Ana María Soilan, Marina Ortega, Pamela Mongelós, Malvina Páez, Amalia Castro, Carmen Cristaldo, Fátima Romina Báez, Claudia Carolina Centurión, Jaime Vester, Hernán Barrios, Griselda Villalba, María Luisa Amarilla, Graciela Giménez.

References

  1. 1. Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin.2018; 68:394–424. pmid:30207593
  2. 2. Walboomers JM, Jacobs MV, Manos MM, Bosch FX, Kummer A, Shah KV, et al. Human papillomavirus is a necessary cause of invasive cancer worldwide. J Pathol 1999; 189:12–9. pmid:10451482
  3. 3. Muñoz N, Bosch FX, de Sanjose S, Herrero R, Castellsagué X, Shah KV, et al. Epidemiologic classification of human papillomavirus type associated with cervical cancer. N Engl J Med 2003; 348:518–27. pmid:12571259
  4. 4. de Sanjose S, Quint WG, Alemany L, Geraets DT, Klaustermeier JE, Lloveras B, et al. Human papillomavirus genotype attribution in invasive cervical cancer: a retrospectrive cross-sectional worldwide study. Lancet Oncol 2010; 11:1048–56. pmid:20952254
  5. 5. Rolon PA, Smith JS, Muñoz N, Klug SJ, Herrero R, Bosch X, et al. Human papillomavirus infection and invasive cancer in Paraguay. Int J Cancer 2000; 85:486–91. pmid:10699919
  6. 6. Kasamatsu E, Cubilla AL, Alemany L, Chaux A, Tous S, Mendoza L, et al. Type-specific human papillomavirus distribution in invasive cervical carcinomas in Paraguay: a study of 432 cases. J Med Virol 2012; 84:1628–35. pmid:22930512
  7. 7. Mendoza LP, Arbiza J, Páez M, Kasamatsu E, Castro A, Giménez G, et al. Ditribution of human papillomavirus genotypes in paraguayan women according to severity of cervical lesión. J Med Virol 2011; 83:1351–57. pmid:21618554
  8. 8. Kjær SK, Frederiksen K, Munk C, Iftner T. Long-term absolute risk of cervical intraepithelial neoplasia grade 3 or worse following human papillomavirus infection: role of persistence. J Natl Cancer Inst. 2010; 102 (19):1478–88. pmid:20841605
  9. 9. Rositch A, Koshiol J, Hudgens M, Razzaghi H, Backes DM, Pimienta JM, et al. Patterns of persistent genital human papillomavirus infection among women worldwide: a literature review and meta-analysis. Int J Cancer. 2013; 133(6): 1271–85. pmid:22961444
  10. 10. Hildesheim A, Herrero R, Castle PE, Wacholder S, Bratti MC, Sherman ME, et al. HPV co-factors related to the development of cervical cancer: results from a population-based study in Costa Rica. Br J Cancer 2001; 84(9):1219–26. pmid:11336474
  11. 11. Muñoz N, Franceschi S, Bosetti C, Moreno V, Herrero R, Smith JS, et al. Role of parity and human papillomavirus in cervical cancer: the IARC multicentric case-control study. Lancet. 2002, 359:1093–101. S0140-6736(02)08151-5. pmid:11943256
  12. 12. Almonte M, Albero G, Molano M, Carcamo C, García PJ, Pérez G. Risk factors for human papillomavirus exposure and co-factors for cervical cancer in Latin America and the Caribbean. Vaccine. 2008; 19; 26 Suppl 11:L16–36.
  13. 13. Muñoz N, Franco E, Herrero R, Andrus JK, de Quadros C, Goldie SJ, et al. Recommendations for cervical cancer prevention in Latin America and the Caribbean. Vaccine 2008; 26S:L96–L107.
  14. 14. Murillo R, Almonte M, Pereira A, Ferrer E, Gamboa OA, Jerónimo J, et al. Cervical cancer screening programs in Latin America and the Caribbean. Vaccine 2008; 26S:L37–L48.
  15. 15. Sankaranarayanan R, Nene BM, Shastri S, Jayant K, Muwonge R, Budukh A, et al. HPV screening for cervical cancer in rural India. N Engl J 2009; 360:1385–94.
  16. 16. Almonte M, Murillo R, Sanchez GI, Jerónimo J, Salmerón J, Ferreccio C, et al. Nuevos paradigmas y desafíos en la prevención y control del cancer de cuello uterino en America Latina. Salud Pública Mex 2010,52:544–49. pmid:21271014
  17. 17. Lazcano-Ponce E, Lorincz AT, Salmerón J, Fernández I, Cruz A, Hernández P, et al. A pilot study of HPV DNA and cytology testing in 50,159 women in the routine Mexican Social Security Program. Cancer Causes Control 2010, 21:1693–700. pmid:20617376
  18. 18. Arbyn M, Ronco G, Anttila A, Meijer C, Poljak M, Ogilvie G, Koliopoulus G, et al. Evicence regarding human papillomavirus testing in secondary prevention of cervical cancer. Vaccine 2012; 30. suppl 5: F88–F99.
  19. 19. Castle PE, de Sanjose S, Quiao Y. Belinson JL, Lazcano-Ponce E, Kinney W. Introduction of human papillomavirus DNA screening in the world: 15 years of experience. Vaccine 2012, 305:F117–F122.
  20. 20. WHO guidelines for screening and treatment of precancerous lesions for cervical cancer prevention, 2013.
  21. 21. Ministerio de Salud Pública y Bienestar Social, Paraguay. OPS/OMS. Manual nacional de normas y procedimientos para la prevención y el control de cáncer de cuello uterino. 2010.
  22. 22. DGEEC. Proyección de la población por sexo y edad, según distrito para el 2018. pp. 452; 456; 466.
  23. 23. Vaccarella S, Franceschi S, Herrero R, Muñoz N, Snijders PJ, Clifford GM, et al. Sexual behavior, condom use, and human papillomavirus: pooled analysis of the IARC human papillomavirus prevalence surveys. Cancer Epidemiol Biomarkers Rev. 2006; 15(2):326–33.
  24. 24. Roura E, Ifter T, Vidart JA, Kjaer SK, Bosch FX, Muñoz N, et al. Predictors of human papillomavirus infection in women undergoing routine cervical cancer screening in Spain: the CLEOPATRE study. BMC Infect Dis 2012; 26: 145.
  25. 25. Vaccarella S, Herrero R, Dai M, Snijders PJ, Meijer CJ, Thomas O, et al. Reproductive factors, oral contraceptive use, and human papillomavirus infection: pooled analysis of the IARC HPV prevalence surveys, Cancer Epidemiol. Biomark. Prev 2006: 15:2148–53.
  26. 26. Castellsagué X, Muñoz N. Chapter 3: Cofactors in human papillomavirus carcinogenesis-role of parity, oral contraceptives, and tobacco smoking. J Natl Cancer Inst Monogr 2003 (31):20–8 pmid:12807941
  27. 27. Mendoza L, Mongelós P, Páez M, Castro A, Rodríguez-Riveros MI, Giménez G, et al. Human papillomavirus and other genital infections in indigenous women from Paraguay: a cross-sectional analytical study. BMC Infectious Dis 2013, 13:531–39.
  28. 28. Bruni L, Diaz M, Castellxagué X, Ferrer E, Bosch FX, de Sanjose S. Cervical human papillomavirus prevalence in 5 continents: meta-analysis of 1 millon women with normal cytology. JID 2010; 202:1789–99. pmid:21067372
  29. 29. Lazcano-Ponce E, Herrero R, Muñoz N, Cruz A, Shah KV, Alonso P, Hernández P, et al. Epidemiology of HPV infection among Mexican women with normal cervical cytology. Int J Cancer 2001; 912:412–20.
  30. 30. de Sanjose S, Diaz M, Castellsagué X, Clifford G, Bruni L, Muñoz N, et al. Worldwide prevalence and genotype distribution of cervical human papillomavirus DNA in women with normal cytology: a meta-analysis. Lancet Infect Dis 2007, 7:453–59. pmid:17597569
  31. 31. Ginindza TG, Diamini X, Almonte M, Herrero R, Jolly PE, Tsoka-Gwegweni J, et al. Prevalence of and associated risk factors for high risk human papillomavirus among sexually active women, Swaziland. Plos ONE 2017; 12(1): e0170189. pmid:28114325
  32. 32. Molano M, Posso H, Weiderpass E, van den Brule AJC, Ronderos M, Franceschi S, et al. Prevalence and determinants of HPV infection among Colombian women with normal cytology. Br. J Cancer, 2002:87:324–33. pmid:12177803
  33. 33. Ferrecio C, Prado RB, Luzoro AV, Ampuero S, Snijders PJ, Meijer CJ, et al. Population-based prevalence and age distribution of human papillomavirus among women in Santiago, Chile. Cancer Epidemiol Biomark Prev 2004,13:2271–76.
  34. 34. Herrero R, Hildesheim A, Bratti C, Sherman ME, Hutchinson M, Morales J, et al. Population-based study of human papillomavirus infection and cervical neoplasia in rural Costa Rica. J Nat Cancer Inst 2000, 92:464–74. pmid:10716964
  35. 35. Herrero R, Castle PE, Schiffman M, Bratti MC, Hildesheim A, Morales J, et al. Epidemiologic profile of type-specific human papillomavirus infection and cervical neoplasia in Guanacaste, Costa Rica. J Infect Dis 2005, 191:1796–807. pmid:15871111
  36. 36. Ortiz AP, Romaguera J, Pérez CM, González D, Muñoz C, González L, et al. Prevalence, genotyping and correlates of anogenital HPV infection in a population-based sample of women in Puerto Rico. Papillomavirus Res 2016, 2:89–96. pmid:29074191
  37. 37. Almonte M, Ferreccio C, González M, Delgado JM, Buckley H, Luciani S, et al. Risk factors for high-risk human papillomavirus infection and cofactors for high-grade cervical disease in Peru. Int J Gynecol Cancer 2011; 21 (9): 1654–63. pmid:21892094
  38. 38. Muwonge R, Ngo Mbus L, Ngoma T, Gombe Mbalawa Ch, Dolo A, Ganda Manuel M, et al. Socio-demographic and reproductive determinants of cervical neoplasia in seven sub-Sahara African countries. Cancer Causes Control 2016; 27 (12) 1437–46. pmid:27822586
  39. 39. Roura E, Travier N, Waterboer T, de Sanjose S, Bosch FX, Pawlita M et al. The influence of hormonal factors on the risk of developing cervical cancer and pre-cancer: results from the EPIC Cohort. Plos ONE 2016,25:11 (1) e0147029. pmid:26808155
  40. 40. International Collaboration of Epidemiological Studies of Cervical Cancer. Cervical Cancer and hormonal contraceptives: collaborative reanalysis of individual data for 16 573women with cervical cancer and 35 509 women without cervical cancer from 24 epidemiological studies. Int J Cancer 2007; 119(5): 1108–24.
  41. 41. Franceschi S, Plummer M, Clifford G, de Sanjose S, Bosch X, Herrero R, et al. Differences in the risk of cervical cancer in human papillomavirus infection by educational level. Br J Cancer 2009; 101:865–70. pmid:19654578
  42. 42. Paez M, Rodríguez-Riveros MI, Kasamatsu E, Castro A, Orué E, Lampert N, et al. Knowledge, attitudes and practice on human papillomavirus (HPV) and cervix cancer in women who were 30 years old or more, from a riverside neighborhood of Asunción (Bañado Sur).2012. Rev Univ Ind Santander Salud 2016; 48 (1):37–44.