Skip to main content

Methylation estimates the risk of precancer in HPV-infected women with discrepant results between cytology and HPV16/18 genotyping



Vigilant management of women with high-risk human papillomavirus (hrHPV) is necessary in cancer screening programs. To this end, we evaluated the performance of S5 (targeting DNA methylation in HPV16, HPV18, HPV31, HPV33, and human gene EPB41L3) to predict cervical intraepithelial neoplasia grade 2 or higher (CIN2+) in a sample of hrHPV-infected women referred to colposcopy in the FRIDA Study, a large screening trial in Mexico. A nested case-control sample with women referred to colposcopy either by atypical squamous cells of undetermined significance or higher (ASCUS+) in cytology and/or positive for HPV types 16 or 18 was tested by S5. Seventy-nine cases of CIN2+ were age-matched to 237 controls without a diagnosis of CIN2+ (<CIN2). DNA from exfoliated cervical cells was bisulfite converted and PCR amplified for S5 targets, and methylation was quantified at specific cytosines by pyrosequencing.


The S5 classifier separated women with CIN2+ from <CIN2 with a highly significant area under the curve (AUC) of 0.75 (95% CI 0.69–0.82), while AUC for CIN3+ was 0.81 (95% CI 0.74–0.89). To optimize sensitivity and specificity for Mexico, an alternative S5 cutoff of 3.7 was implemented to account for overall higher methylation seen in our already triaged women. All three invasive cancers were detected by methylation or HPV16/18 but none by cytology. Sensitivity of S5 for CIN2+ was 62% (95% CI 50.4–72.7%), specificity was 73% (95% CI 66.9–78.5%), and adjusted PPV was 15.1% (95% CI 12.0–18.3%). In contrast, the crude sensitivity of HPV16/18 detection and cytology were 63.3% (95% CI 51.7–73.9%) and 57.0% (95% CI 45.3–68.1%) respectively; specificity was 29.1% (95% CI 23.4–35.3%) and 62.4% (95% CI 55.9–68.6%) respectively, while adjusted PPV was 6.4% (95% CI 4.9–8.1%) and 10.5% (95% CI 8.0–13.1%), respectively. Methylation testing could reduce colposcopy referrals by 30 to 50% with virtually no loss of sensitivity for CIN2+ and CIN3+.


S5 testing on hrHPV-positive women significantly increased diagnostic information compared to triage by HPV16/18 plus cytology and appears to have clinical utility as an additional test to substantially lessen burdens on colposcopy.

Trial registration

The FRIDA Study is registered in, number NCT02510027.


Cervical cancer is the fourth leading cause of cancer-related death in women worldwide, and developing regions carry the greatest burden of this disease. In Mexico, cervical cancer is ranked second as a cause of death from malignancy in women [1,2,3,4]. The recognition of persistent infection with high-risk human papillomavirus (hrHPV) as a major cause of cervical cancer [5,6,7,8] has favored the development of new technologies for hrHPV detection. The introduction of these tests as part of primary screening for cervical cancer in various countries has greatly improved the detection of precancerous lesions [9,10,11,12]. Some new testing procedures represent a breakthrough in screening [9, 13, 14]; however, hrHPV testing programs must incorporate triage before colposcopy referral because of the low specificity of hrHPV assays [15, 16]. The use of triage tests should enable the relatively accurate identification of the small proportion of hrHPV-positive women who need to go to colposcopy given their increased risk of developing cancer. Triage can reduce over-diagnosis and over-treatment of clinically non-relevant hrHPV infections [3, 17,18,19]. Various types of triage tests have been proposed, including cytology, p16/Ki67 immunocytochemistry, and HPV16/18 detection, but these methods still have important limitations such as subjectivity of interpretation, low sensitivity, and low positive predictive value (PPV) [15, 16, 20].

DNA methylation, the enzymatic addition of a methyl group to the carbon at position 5 of the cytosine ring in a cytosine-phosphate-guanine (CpG) site, has been shown to play an important role in the development of the carcinogenic process [21,22,23,24,25]. Many studies have found a strong association between the methylation of host and viral genome with the development of CIN2, CIN3, and cancer [22, 23, 26,27,28,29]. Indeed, various quantitative combination methylation assays are currently under evaluation with the most common host genes being EPB41L3, MAL, CADM, FAM19A4, and MIR124, as well as CpG sites in the late regions of various HPV genomes [24, 30,31,32,33,34,35]. According to Clarke et al., methylation across all 12 carcinogenic HPV types is associated with precancer (CIN3+), providing a better sensitivity than cytology as triage test [36]. Other methylation marker panels designed from methylation-specific real-time PCR assays are now being validated with promising results [37], in particular, the QIAsure methylation test, which is a multiplex real-time methylation-specific PCR [38], and the GynTect assay based on the detection of DNA methylation of human marker gene regions [37,38,39]. In our previous studies, we have consistently found that methylation assays targeting the most carcinogenic HPV types help to distinguish progressing infections and more severe lesions [22, 27, 32]. Indeed, we have found that the performance of methylation to detect precancer improves as more HPV viral types and human genes are added [30, 40]. In particular, the development of a methylation score—called S5—combining DNA methylation levels of HPV16, HPV18, HPV31, HPV33, and the human tumor suppressor gene EPB41L3 has arisen as a promising biomarker panel that can be used to triage hrHPV positive women. At 90% sensitivity, S5 provided a specificity of 49% to identify or predict the risk of a high-grade cervical lesion in a UK colposcopy referral population [31]; in other words, one half of the women sent for colposcopy (regarded as an invasive examination) could have been spared costs, a biopsy, and anxiety while being reassured of highly sensitive detectability if they had a precancer. The S5 classifier was developed in a colposcopy study [31] and validated in screening studies in the UK and Canada [28, 32]. The aim of our study was to extend validation of the diagnostic performance of S5 and to find new applications in geographically and ethnically different women undergoing triage in a large population-based HPV screening cohort in Mexico. In this study, we explored the utility of methylation testing in women with discrepant results between cervical cytology and HPV16/18 genotyping. We give the term “already-triaged” to the women in our study to avoid possible confusion over the role of an additional methylation test.


The 79 cases and 237 controls were composed mostly of women aged between 30 and 49 (92.4%), with a median age of 37 years. Figure 1 shows the consort diagram of our nested case-control study: 29,759 cervical specimens were tested for hrHPV (Cobas 4800). The grouped prevalence of 13 types of hrHPV plus HPV66 (now regarded as a low-risk HPV) was 10.8%. From 3228 hrHPV-positive women, 863 (26.7%) were positive on initial triage (HPV16/18+ and/or ASCUS+); 672 of these women were referred to colposcopy, and 561 attended for colposcopy prior to our nested case-control study cutoff date of July 2015; the colposcoped women were the subjects of our methylation study.

Fig. 1
figure 1

Consort diagram of FRIDA nested case-control triage study showing the numbers of women in each step. Triage positive women included HPV16/18 positive (<ASCUS) (n = 508), ASCUS+ (HPV16/18 negative) (n = 277), and HPV16/18 positive and ASCUS+ (n = 78). *Prior to July 2015, 672 out of 863 triage-positive women attended the colposcopy evaluation. Then, 561 women out of 672 who underwent colposcopy had histology results recorded from April 2013 to the time of the study cutoff date on July 15. These 561 women are represented in our sampling frame, from which we selected all the CIN2+ cases (79) and a random selection of three controls per case matched by age (1 CIN1 and 2 NEG). In total, 316 participants were included in our methylation analysis: 79 CIN2+ cases, 79 CIN1, and 158 NEG. The cases included all the CIN2+ detected until July 2015. The three controls per case were randomly selected and matched by age from women with histological diagnoses of CIN1 or less. The remaining 245 triage-positive women with histology results of CIN1 and negative would have been controls but were not selected by the sampling method as we already had adequate power for the study. Abbreviations: hrHPV, high-risk human papillomavirus; HPV16/18, human papillomavirus type 16 or type 18; ASCUS+, atypical squamous cells of undetermined significance or worse; NEG, histologically negative; CIN, cervical intraepithelial neoplasia (of grades 1, 2, and 3)

There were no significant differences with respect to demographics and lifestyle characteristics between cases and controls (Additional file 1). Twenty-seven percent of cases and 30% of controls began their sexual life before the age of 16. The majority of cases and 51% of controls reported having had one lifetime sexual partner. The proportion of women who had used hormonal contraceptives for more than 12 months was the same (12.7%) in both groups. Most women in both groups reported not having used a condom during the last 12 months (75% cases and 74% controls). 48.1% of cases and 36.7% of controls had four or more pregnancies before FRIDA enrollment. Overall, 90% of the study population had never smoked (89.8% cases and 93.3% controls).

Focusing on the HPV results, 60.8% of cases and 46.8% of controls were HPV16-positive (p = 0.038). However, only 6.3% of cases as compared to 26.6% of controls were HPV18-positive (p = 0.000). With respect to other hrHPV types, 24.1% of cases and 19.4% of controls were positive for HPV31 (p = 0.422). Similarly, 6.3% of cases and 5.9% of controls were HPV33-positive (p = 1.00). Fifty-five percent of both cases and controls were positive for other types of hrHPV plus HPV66. With respect to cytology, 57% of cases and 38% of controls showed ASCUS+ (p = 0.004) (Table 1).

Table 1 hrHPV type-specific prevalence and cytological evaluation in cases and controls

The S5 classifier showed a highly significant increase proportional to the severity of lesions (Cuzick test for trend, p < 0.0001). Median methylation was 1.3 in histopathologically negative samples (NEG), 1.7 in CIN1, 4.4 in CIN2, 6.5 in CIN3, and 19.2 in cervical cancer (CC) (Additional file 2). Figure 2 shows the distribution of the S5 classifier by histopathological diagnosis. The Mann-Whitney U test revealed highly significant methylation differences in the following pairwise comparisons: NEG vs CIN2 (p = 0.01), NEG vs. CIN3 (p < 0.001), NEG vs CC (p < 0.001), CIN1 vs CIN2 (p = 0.03), CIN1 vs. CIN3 (p < 0.001), CIN1 vs CC (p = 0.004), CIN2 vs CC (p = 0.007), and CIN3 vs CC (p = 0.02).

Fig. 2
figure 2

Comparison of S5 methylation classifier in histologically negative (NEG), CIN1, CIN2, CIN3, and cervical cancer cases (CC). The S5 classifier was significantly different between the following group comparisons: NEG vs CIN2 (p = 0.01), NEG vs CIN3 (p < 0.001), NEG vs CC (p < 0.001), CIN1 vs CIN2 (p = 0.03), CIN1 vs CIN3 (p < 0.001), CIN1 vs CC (p = 0.004), CIN2 vs CC (p = 0.007), and CIN3 vs CC (p = 0.02). Other comparisons were not significant (NEG vs CIN1 and CIN2 vs CIN3). Abbreviations: NEG, histologically negative; CIN, cervical intraepithelial neoplasia (of grades 1, 2, and 3); CC, cervical cancer. The top of box represents the upper quartile (p75), bottom the lower quartile (p25), and the line the median (p50). The upper whisker extends to the largest point of the inter-quartile range from the upper quartile. The lower whisker extends to the smallest point of the inter-quartile range from the lower quartile. The outliers are plotted as individual points for each lesion grade. The Cuzick test for trend was highly significant (p < 0.001)

ROC analysis of the S5 classifier for detecting CIN2+ and CIN3+ gave areas under the curve (AUC) of 0.75 (95% CI 0.69–0.82) and 0.81 (95% CI 0.74–0.89) respectively (Fig. 3). A new S5 cutoff value of 3.7 was selected for better discriminating CIN2+ lesions from <CIN2 diagnoses in our already-triaged Mexican women with abnormal cytology and/or HPV16/18 infection (Fig. 3). This cutoff yielded a relative sensitivity and specificity of 62.0% (95% CI 50.4–72.7) and 73.0% (95% CI 66.9–78.5) respectively for CIN2+ (Table 2). The same cutoff for discriminating CIN3+ from <CIN3 produced a sensitivity of 70.3% (95% CI 56.7–81.8) and a specificity of 76.6% (95% CI 70.7–81.9). In addition, we calculated the performance of the S5 classifier with a cutoff of 0.8 which was previously validated in the UK and Canadian screening populations. We observed a better sensitivity for CIN2+ of 86.1% (95% CI 76.5–92.8) (p = 0.0005), but a worse specificity of 40.1% (95% CI 33.8–46.6) (p < 0.0001) at the 0.8 S5 cutoff.

Fig. 3
figure 3

Receiver operator characteristic (ROC) and area under the curve (AUC) of S5 methylation for detecting CIN2+ or CIN3+. The blue diamond denotes the sensitivity and specificity of S5 at a cutoff of 3.7 for CIN2+. The red circle denotes the sensitivity and specificity at the S5 cutoff of 0.8 predefined for use in the UK for CIN2+. The cutoff for cytology alone was ASCUS+. Abbreviations: ASCUS+, atypical squamous cells of undetermined significance or worse

Table 2 Performance of the S5 methylation classifier, cytology, and HPV16/18 genotyping for detecting CIN2+ and CIN3+

Of the 316 triage-positive samples, 218 were positive for HPV16/18 and 134 were ASCUS+ (32 were HSIL). The unadjusted sensitivity of HPV16/18 genotyping for detecting CIN2+ was 63.3% (95% CI 51.7–73.9%), with a specificity of 29.1% (95% CI 23.4–35.3, Table 2). Cytology alone with a cutoff of ASCUS+ had an unadjusted sensitivity of 57.0 (95% CI 45.3–68.1) and an unadjusted specificity of 62.4% (95% CI 55.9–68.8), while cytology with a cutoff of HSIL+ had an unadjusted sensitivity and specificity of 27.8 (95% CI 18.3–39.1) and 95.8% (95% CI 92.8–98.0) respectively.

We investigated the positivity of different triage tests for NEG, CIN1, CIN2, CIN3, and CC histopathological diagnoses (Table 3). In the context of cervical cancer screening in Mexico, a greater specificity would reduce the number of women called for colposcopy (Table 3). For example, HPV16/18 genotyping was positive in 72.2% of histologically negative women as opposed to 27.2% positive with S5 classifier at a cutoff of 3.7. All three cancers were identified by HPV16/18 detection and by the S5 classifier, but all were missed by cytology.

Table 3 Positivity of different triage tests for histologically negative (NEG), CIN1, CIN2, CIN3, and CC

We subsequently explored the clinical utility of the S5 methylation classifier as a second or reflex triage test (Fig. 4). The women in our study had increased risk of CIN2+ whether they were positive for either or both ASCUS and HPV16/18; therefore, we defined three risk groups of women based on the outcomes of the latter two tests. The first group consisted of HPV16/18-positive and ASCUS-positive women. We consider that these women should be referred to colposcopy based on the first triage alone. A second group represented a triage discrepancy consisting of women negative for HPV16/18 but positive for an ASCUS+ result. In our study, these women were also called for colposcopy, but in retrospect, some of them were considered to have gone unnecessarily. The use of S5 methylation as a second triage test in this group would have reduced colposcopy referral by 50% for a CIN2+ endpoint (p = 0.0000) or 43% for a CIN3+ endpoint (p = 0.0000). In the third group, women were HPV16/18 positive but had normal cytology. All these women were also called for colposcopy. Here again, the use of S5 methylation would have reduced colposcopy referral by 30% (CIN2+ endpoint, p = 0.0000) (Fig. 4a) or 28% (CIN3+ endpoint, p = 0.0000) (Fig. 4b). In both groups 2 and 3, the sensitivity of S5 as the second triage test was statistically no different than in a scenario where all women were sent to colposcopy (Table 3 and Fig. 4). This means that using S5 as a reflex triage test would reduce the number of women unnecessarily sent to colposcopy by 30 to 50% while maintaining the same sensitivity to detect CIN2+ (or CIN3+) as a triage based on a combination of HPV16/18 and cytology.

Fig. 4
figure 4

Benefit of using the S5 classifier as a second triage test for colposcopy referral for a) CIN2+ and b) CIN3+ endpoints. S5 helped reduce unnecessary colposcopy referrals in both the “HPV16/18 neg, ASCUS+” and “HPV16/18 pos, normal cytology” groups. Under the current Mexican algorithm, all HPV16/18 positive and/or ASCUS positive are referred to colposcopy, but we propose to use S5 as second triage test for the discrepant triage groups (“HPV16/18 neg, ASCUS+” and “HPV16/18 pos, normal cytology”) which then decreases the false-positive rate. We therefore defined the following hypothesis to test S5 benefits. In group 1, women were HPV16/18 positive and ASCUS positive and must be referred to colposcopy without any additional triage procedure. In group 2, women were HPV16/18 negative and ASCUS positive and were called for colposcopy, but if we had used S5 methylation as a second triage test, it would have reduced by 50% (CIN2 endpoint) or 43% (CIN3 endpoint) the number of false positives referred to colposcopy. Using S5 methylation as a second triage for women in group 3, who were HPV16/18 positive with normal cytology, would have reduced unnecessary colposcopy referrals by 30% in CIN2+ and 28% in CIN3+. Abbreviations: ASCUS+, atypical squamous cells of undetermined significance or worse; CIN, cervical intraepithelial neoplasia (of grades 1, 2, and 3). The frequency shows the absolute number of women in each group


There is increasing evidence that cervical precancer arises primarily from persistent hrHPV infection in concert with reasonably well-understood abnormal epigenetic events. Host cells can respond to HPV infection with a cascade of defense mechanisms such as induction of interferon pathways, increased APOBEC-related DNA editing, and increased methylation. Such activity can get out of control and result in genome integrity problems, which in pluripotent cells have dangerous effects. Since the defense mechanisms are unable to subdue persistent HPV infections in a timely manner, they may instead participate in concert with continued over-expression of the HPV E6 and E7 oncoproteins as drivers towards cancer [28, 29, 41].

DNA methylation is a particularly tractable feature of the epigenetic machinery because patterns of CpG methylation are measurable as signatures and are stably retained during mitotic cell division. This feature encouraged us to explore the clinical utility of the S5 DNA methylation classifier as a triage test in our study of women with abnormal cytology and/or HPV16 and HPV18 infections. An important finding of our study is that the S5 classifier, at a modified cutoff of 3.7 (more favorable for women already triaged by cytology or HPV16/18 testing), was able to detect 62.0% of CIN2+, 70.3% of CIN3+, and all cancers in our study. In particular, S5 was able to provide important additional information for women with discrepant results between HPV16/18 and cytology. We believe that a combined triage approach that is cost-effective, with good specificity and a sensitivity of greater than 60% for CIN2/3, is generally safe for use in Mexico and also in women residing in low- and middle-income countries (LMIC), because of two important considerations: firstly, it is not crucial to detect all precancers in the first screen because most of these lesions regress or progress quite slowly; however, cancers are not reversible, and in our study, all cancers were detected by S5. Furthermore, in two other studies, S5 was also shown to have a ~ 100% sensitivity for cervical cancer [28, 32]. Secondly, additional intensive follow-up of all triage test-negative women is a standard recommendation, meaning that most if not all cancers in development can be caught at an early stage when cure rates are high (> 95%) as long as the women alerted to follow-up adhere to recommendations. A pressing priority in Mexico and indeed in most LMIC is to not overwhelm colposcopy clinics. A workload increase in these clinics could lead to increasingly poor outcomes and potentially to increases in cancer due to hurried and inadequate assessment of women, who may be given false-negative diagnoses or who may be lost to follow-up [42, 43]. Our results support a proposal to create enhanced risk scores that combine methylation of target genes with current cytology and HPV16/18 triage testing options for hrHPV screen-positive women.

Another major finding of our study is confirmation of the diagnostic utility of S5 seen in two earlier validation studies, a large screening study in the UK and a randomized control trial in Canada [28, 32]. In the hrHPV-based screening population from the UK, using a cutoff of 0.8 for S5, the AUC obtained for CIN2+ was 0.78 (95% CI 0.69–0.88), similar to the AUC of 0.75 (95% CI 0.69–0.8) seen in our FRIDA population. The 0.8 cutoff used in the UK population gave a higher sensitivity but a lower triage specificity compared with our S5 cutoff of 3.7. In previous studies of S5 performance, the sensitivity for cancer at the suggested triage cutoff of 0.8 was 100%. In our study, we showed that even at a cutoff of 3.7, all cancers were detected. Similar data have also been shown in other studies [28, 32], which indicates that few if any cancers would be missed when employing DNA methylation signatures for triage, even with a rather high S5 cutoff of 3.7 or perhaps higher.

Most lesions resulting from an infection with hrHPV eventually regress and only a very small number (< 5%) ever go on to invasive cancer in the absence of screening [3, 6]. To decrease costs and to avoid a large amount of unnecessary treatment, it is of great importance to identify hrHPV-positive women who are most likely to develop truly precancerous high grade CIN and to clearly separate them from women with CIN2 and CIN3 who will not develop cervical cancer at all or at least not over several decades [9, 14, 44]. Previous studies have proposed classifiers based on the methylation of HPV types (classifiers S1, S2 [22], and S4 [30]) and or human genes [40].

The S5 classifier was positive in 27% of our selected controls, but it is worthwhile to recall that these were women with either HPV16/18 infection or an abnormal cytology who were not subsequently diagnosed with CIN2+. Previous studies have observed that DNA methylation can predict the persistence of HPV infection of different viral types [22,23,24, 45]. Thus, it can be suggested that women, for whom a first colposcopy evaluation detected no significant lesions, could be carefully evaluated in follow-up visits to document quantitative changes in methylation consistent with the clearance or progression of disease.

Our results are applicable to methylation triage of women who are discrepantly positive for HPV16/18 and/or cytology ASCUS+; however, we were not able to assess if S5 could be used for triage of all hrHPV-positive women. The use of S5 could reduce colposcopy rates by up to 50% without a loss of sensitivity to detect cancer, which would be a big saving for any LMIC, not just because of lower financial costs but because of overall lower clinical burden. Furthermore, as shown by our results on women with concordant positive results for both HPV16/18 and cytology, the S5 test may also be useful for women who had a normal colposcopy. In this situation, a negative methylation result (S5 = 0) would provide reassurance against missed disease, especially if there is a high risk of loss to follow-up. The new cutoff proposed and different values of high methylation need to be validated in the setting of LMIC to estimate the clinical utility in populations with a similar context of income and disease prevention strategies to estimate the risk of cervical cancer.

The use of our particular triage population underestimates the number of potential cases of CIN2+ attributable to non-HPV16/18 hrHPV types; however, many of these women may be identified on follow-up. We consider that an introduction of S5 DNA methylation testing for routine use in self-validated expert clinical diagnostics laboratories is a reasonable course and, in some cases, could be undertaken quite quickly, being advantageous to further triage of women with HPV16/18 and/or abnormal cytology results. Successful implementation would allow colposcopy clinics to see only the highest-risk women and thereby attain higher quality and better manage overtreatment. Excessive burdening of clinical services and overtreatment of women are negative features of medical practice today due to an increasing flood of new women identified in expanded national HPV screening programs [46]. Based on practical experience from our laboratory, the cost of an S5 DNA methylation test is not more than a routine HPV screening and reflex genotyping. Thus, the S5 test could be applied to a minority of the women, and incremental costs could be balanced against the expected 30 to 50% reduction in colposcopy referral costs. Also, of note, the three cancers were detected by both HPV16/18 testing and the S5 test, but all of them were missed by cytology (Table 3). If confirmed by larger studies, this observation may allow for a generalized rapid reduction in cervical cancer in LMIC, which have been particularly poorly served by cytology screening over the past 50 years. Expected large-scale cervical cancer incidence reductions from implementation of the prophylactic HPV vaccine are still many decades away, and it is prudent to take more preventative actions in the meantime.

Our study has some important strengths including that it is the largest screening multi-triage study of cervical cancer in an LMIC. All primary clinical and laboratory procedures (screening, triage, colposcopy, and histopathological assessments) were done completely within Mexico. We selected all available cases and three matching controls, the latter specimens being selected at the same time and diagnosed in the same corresponding months as the cases. The accuracy of the colposcopic evaluation in FRIDA was very high and was not typical of Mexico, LMICs, or even developed countries such as the USA and Europe because it was performed by specially trained experts and also did not depend only on the visual expertise of the colposcopists. At least one biopsy was taken from each quadrant of the cervix in all women to reduce verification bias. Additionally, the histopathological evaluation and diagnostic confirmation of the biopsies and/or endocervical samples of all women included in our study were reviewed by a panel of pathologists, including an international expert. Moreover, all methylation measurements were done by blinded laboratory staff, including a scientist from Mexico (RH), with the aim of avoiding differential information bias.


The S5 methylation classifier provides good sensitivity, specificity, and positive predictive value as a triage test for colposcopy referral in Mexican women with HPV16/18+ and/or ASCUS+, especially if CIN2/3 with low methylation do not progress, which is a subject of ongoing studies. We believe that our study and the existing literature on this topic support the use of methylation biomarkers in the triage of hrHPV-positive women within cervical cancer prevention programs in both highly developed settings and in LMICs. S5 methylation triage can reduce the numbers of needed colposcopies by an additional 30 to 50% as compared to triage by cytology and HPV16/18 genotyping while detecting all cancers. Such steps would greatly reduce existing costs and allow healthcare staff to better focus on the women at real risk of cervical cancer.


Study population

We performed a nested case-control analysis in a subset of hrHPV-positive women participating in the FRIDA screening study, all of whom underwent colposcopy (Fig. 1). Details of the study methods have been described previously [47]. Briefly, FRIDA is a population-based study that recruited 36,212 women from a target population of close to 60,000 women aged 30–64 years, between April 2013 and February 2016. Women were users of public health services, who attended the Cervical Cancer Screening Program in 100 primary health care facilities in the state of Tlaxcala, Mexico. The main aim of FRIDA was to evaluate the performance of different triage tests for hrHPV-positive women. Women who were pregnant at the time of recruitment or had a prior hysterectomy were excluded. Informed consent was required for women to participate in the study. The population sampling for this nested case-control study was selected from a pool of women attending colposcopy because they were either singly positive or positive for both cytology and HPV16/18 between April 2013 and July 2015, when 29,759 women had been enrolled. Women who were infected by HPV types other than HPV16 or HPV18 who did not have cytological abnormalities were excluded for this analysis.

FRIDA is registered in, number NCT02510027, and was approved by the Institutional Review Boards (IRBs) of the participating institutions: National Public Health Institute (INSP) [1094], Tlaxcala State Ministry of Health [SS.DECI-OI-13/12], and the Mexican regulatory agency COFEPRIS [CAS/OR/01/CAS/123300410C0044-3578/2012].

Selection of cases and controls

The nested case-control triage study considered hrHPV-positive women who attended for colposcopy as a result of a positive test by HPV16/18 (n = 586) or a cytological diagnosis of atypical squamous cells of undetermined significance or worse (ASCUS+) (n = 277). HPV testing for 13 high-risk HPV types and HPV66, followed by HPV16/18 genotyping, were performed by the Cobas® 4800 DNA test (Roche, Pleasanton, CA, USA) as described below. Endocervical curettage was performed for all women referred for colposcopic evaluation, during which one biopsy was also collected from what the colposcopist suspected to be the most abnormal zone of each quadrant. Evaluation of all histological samples was conducted by a standardized panel of pathologists. If two pathologists agreed in their diagnosis, this was the final conclusion. If the diagnoses were discordant, an additional interpretation was made by a third expert pathologist. In this case, the assigned diagnosis was based on the agreement of two out of three pathologists or, if totally discordant, on the result given by the expert pathologist.

Cases and controls were selected based on the histological status. All women diagnosed with CIN2, CIN3, or invasive cervical cancer were defined as cases. For each case identified, we selected three age-matched controls: one control diagnosed with low-grade cervical intraepithelial neoplasia (CIN1) and two controls with no evidence of any cervical intraepithelial lesion (NEG). Controls were selected randomly within a 2-month time frame of the case diagnosis.

Collection and shipment of cervical samples

DNA methylation assays were performed on an aliquot of cervical samples obtained from the screening visit. In this visit, a first cervical sample was collected in a SurePath (Becton Dickinson, Sparks, MD, USA) vial and used for the cytology procedure. Cervical samples for HPV and methylation testing were collected as a second sample and placed in a vial containing ThinPrep® preservative (Hologic, Inc., Bedford, MA). Both samples were collected using a Cervex-Brush® (Rovers®) in the screening visit. The samples were temporarily stored at room temperature at the health center until they were delivered to HPV laboratory facilities where aliquots were made and kept frozen at − 70 °C. The samples were shipped on dry ice to the Molecular Epidemiology Laboratory of the Wolfson Institute of Preventive Medicine (Queen Mary University of London, UK) for analysis.

HPV detection

The samples were tested for hrHPV using the Cobas® 4800 HPV test (Roche Molecular Systems, Pleasanton, CA). This system extracted DNA and amplified by PCR using specific primers to human beta globin and hrHPV. HPV16 and HPV18 are identified individually and other HPV types (16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, 59, 66, and 68) were detected as a pool.

Since the Cobas® assay does not give individual typing of HPV31 and HPV33, the prevalence of these two HPV types was estimated by the detection of amplicons specific for these two HPV types in the PCR reactions (HPV31L1 and HPV33L2) leading up to the methylation assays.

Methylation assays

For methylation assays, the DNA was extracted from cervical samples using a QIAamp mini kit (Qiagen, Hilden, Germany) following the manufacturer’s instructions. Subsequently, the DNA was quantified by UV absorption using a NanoDrop spectrophotometer. An aliquot of 500 ng of DNA underwent bisulfite conversion in which unmethylated cytosines were converted to uracil using the EZ DNA methylation kit (Zymo Research, Irvine, CA) following the manufacturer’s instructions. This kit is designed to reduce degradation and minimize the loss of DNA. The DNA was eluted in 20 μl of Elution Buffer.

After DNA conversion, the PyroMark PCR Kit (Qiagen) was used to amplify the following CpG sites: HPV16 L1 (6367 and 6389), HPV16 L2 (4238, 4247, 4259, 4268, and 4275), HPV18 L2 (4256, 4261, 4265, 4269, 4275, and 4282), HPV31 L1 (6352, 6364), HPV33 L2 (5557, 5560, 5566, and 5572), and human gene EPB41L3 (425, 427, and 438). PCR products were pyrosequenced using a PyroMark Q96 ID instrument (Qiagen), and the proportions of cytosine and thymine were quantified for each CpG site. The laboratory was blinded to the cytology and histology results. The procedures and quality control of these experiments have been reported previously [23, 27].

Data analysis

Socio-demographic, lifestyle, and sexual behavior of cases and controls were summarized as means or proportions. We computed the specific prevalence of HPV types (16, 18, 31, and 33) within our study population. Differences across all variables between cases and controls were tested using a two-sided, difference of proportions by a binomial test or difference of medians by a Mann-Whitney U test. Our primary goal was to assess the clinical performance of the S5 classifier. This score was calculated using the average methylation values of all the target regions using the following standardized equation [31]:

S5 = EPB41L3*(30.9) + HPV16L1.3*(13.7) + HPV16L2*(4.3) + HPV18L2*(8.4) + HPV31L1*(22.4) + HPV33L2*(20.3)

We created boxplots to illustrate the distribution of the S5 classifier according to the histopathological diagnosis of the lesions (NEG, CIN1, CIN2, CIN3, and invasive cervical cancer (CC)). We used the Mann-Whitney U test for comparing S5 differences between different disease categories and the Cuzick test for trend to determine if methylation increased significantly as a function of greater histology result.

In order to maximize the sensitivity and specificity of the S5 score after the triage test (HPV16/18 and cytology) we calculated receiver operating characteristic curves (ROC) to estimate areas under the curve (AUC) and to obtain a new cutoff value that best discriminated hrHPV-positive women with CIN2+ (cases) from controls who did not have histologic abnormalities, in order to maximize the sum of sensitivity and specificity of the S5 score among our already-triaged samples. Based on this new cutoff value, we set a disease endpoint of CIN2+ to estimate sensitivity, specificity, positive predictive value (PPV), and negative predictive value (NPV); the PPVs were also adjusted for CIN2+ and CIN3+ prevalence (10.0% and 6.8% respectively) using the following formula: PPV = (Sn × Pr)/((Sn × Pr) + (1 − Sp) × (1 − Pr)), where Sn is sensitivity, Sp is specificity, and Pr is the CIN2+ or CIN3+ prevalence. A ROC curve was also computed for CIN3+ versus NEG/CIN1. We compared the performance of S5 with HPV16/18 detection and cytology to identify CIN2+ using crude estimators because it was not possible to adjust for the bias created by not referring all hrHPV positive women to colposcopy; in effect, HPV16/18 diagnostic performance was assessed versus a comparator composed of HPV16/18 plus cytology ASCUS+; similarly, ASCUS+ cytology performance was compared to the same composite comparator. A second analysis was done using the established 0.8 cutoff validated earlier in a UK screening population for detection of CIN2+ and CIN3+. All p values were estimated as two sided, with a confidence interval of 95%. Statistical analyses were performed using STATA software v.14.0 (StataCorp LP, College Station, TX).

Availability of data and materials

The full datasets generated from and/or analyzed during the current study are not publicly available due to the need to protect participant confidentiality; however, selected anonymized data can be made available by the FRIDA steering committee. Inquiries should be communicated to the corresponding author who will work with other authorized members of the FRIDA team to consider all sufficiently specified and reasonable requests.



Atypical squamous cells of undetermined significance


Area under the curve


Cervical cancer


Confidence interval


Cervical intraepithelial neoplasia grade 1


Cervical intraepithelial neoplasia grade 2


Cervical intraepithelial neoplasia grade 2 or higher


Cervical intraepithelial neoplasia grade 3


Cytosine-phosphate-Guanine site


High-risk human papillomavirus


High-grade squamous intraepithelial lesions


Late region 1


Late region 2


Low-grade squamous intraepithelial lesions


Negative for cervical intraepithelial neoplasia


Negative predictive value


Positive predictive value


Receiver operating characteristic


  1. IARC. GLOBOCAN2018 [Internet]. [cited 2018 Dec 22]. Available from:

  2. Popat K, McQueen K, Feeley TW. The global burden of cancer. Best practice & research Clinical anaesthesiology. 2013;27(4).

  3. Schiffman M, Castle PE, Jeronimo J, Rodriguez AC, Wacholder S. Human papillomavirus and cervical cancer. Lancet. 2007;370(9590).

    Article  CAS  Google Scholar 

  4. Lazcano-Ponce E, Lorincz AT, Salmeron J, Fernandez I, Cruz A, Hernandez 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(10).

    Article  Google Scholar 

  5. Walboomers JM, Jacobs MV, Manos MM, Bosch FX, Kummer JA, Shah KV, et al. Human papillomavirus is a necessary cause of invasive cervical cancer worldwide. Jounal Pathology. 1999;189(1).

  6. Franco EL, Duarte-Franco E, Ferenczy A. Cervical cancer: epidemiology, prevention and the role of human papillomavirus infection. Cmaj. 2001;164:7.

    Google Scholar 

  7. Munoz N, Bosch FX, de Sanjose S, Herrero R, Castellsague X, Shah KV, et al. Epidemiologic classification of human papillomavirus types associated with cervical cancer. New England Journal Medicine. 2003;348(6).

    Article  Google Scholar 

  8. Lorincz AT, Reid R, Jenson AB, Greenberg MD, Lancaster WKR. Human papillomavirus infection of the cervix: relative risk associations of 15 common anogenital types. Obstet Gynecol. 1992;9:3.

    Google Scholar 

  9. Ratnam S, Franco EL, Ferenczy A. Human papillomavirus testing for primary screening of cervical cancer precursors. Cancer Epidemiol Biomarkers Prev. 2000;9:9.

    Google Scholar 

  10. Campos NG, Maza M, Alfaro K, Gage JC, Castle PE, Felix JC, et al. The comparative and cost-effectiveness of HPV-based cervical cancer screening algorithms in El Salvador. International journal of cancer. 2015;137(4).

    Article  CAS  Google Scholar 

  11. Ronco G, Dillner J, Elfström KM, Tunesi S, Snijders PJF, Arbyn M, et al. Efficacy of HPV-based screening for prevention of invasive cervical cancer: follow-up of four European randomised controlled trials. The Lancet. 2014;383(9916).

    Article  Google Scholar 

  12. Castle PE, Schiffman M, Wheeler CM, Wentzensen N, Gravitt PE. Impact of improved classification on the association of human papillomavirus with cervical precancer. American journal of epidemiology. 2010;171(2).

    Article  Google Scholar 

  13. Schiffman M, Wheeler CM, Castle PE. Human papillomavirus DNA remains detectable longer than related cervical cytologic abnormalities. Journal of Infect Dissease. 2002;186(8).

    Article  CAS  Google Scholar 

  14. Tota JE, Chevarie-Davis M, Richardson LA, Devries M, Franco EL. Epidemiology and burden of HPV infection and related diseases: implications for prevention strategies. Prev Med. 2011;53:1.

    Article  Google Scholar 

  15. Arbyn M, Sasieni P, Meijer CJLM, Clavel C, Koliopoulos G, Dillner J. Chapter 9: Clinical applications of HPV testing: a summary of meta-analyses. Vaccine. 2006;24(Suppl. 3):78–89.

    Article  Google Scholar 

  16. Rao A, Sandri MT, Sideri M, Young S, Sharma A, Behrens C. Comparison of hybrid capture 2 high-risk HPV results in the low positive range with cobas® HPV test results from the ATHENA study. Journal of Clinical Virology. 2013;58(1).

    Article  Google Scholar 

  17. Fields RP. Human papillomavirus DNA versus papanicolaou screening tests for cervical cancer. New England Journal Medicine. 2008;358(6).

  18. Flores YN, Bishai DM, Lorincz A, Shah KV, Lazcano-Ponce E, Hernandez M, et al. HPV testing for cervical cancer screening appears more cost-effective than Papanicolau cytology in Mexico. Cancer Causes Control. 2011;22(2).

    Article  Google Scholar 

  19. Mayrand MH, Duarte-Franco E, Rodrigues I, Walter SD, Hanley J, Ferenczy A, et al. Human papillomavirus DNA versus Papanicolaou screening tests for cervical cancer. New England Journal Medicine. 2007;357(16).

    Article  CAS  Google Scholar 

  20. Huh WK, Williams E, Huang J, Bramley T, Poulios N. Cost effectiveness of human papillomavirus-16/18 genotyping in cervical cancer screening. Applied Health Economics and Health Policy. 2015;13(1).

    Article  Google Scholar 

  21. Lai HC, Lin YW, Huang TH, Yan P, Huang RL, Wang HC, et al. Identification of novel DNA methylation markers in cervical cancer. International Journal of Cancer. 2008;123(1):161–167.

    Article  CAS  Google Scholar 

  22. Lorincz AT, Brentnall AR, Vasiljevic N, Scibior-Bentkowska D, Castanon A, Fiander A, et al. HPV16 L1 and L2 DNA methylation predicts high-grade cervical intraepithelial neoplasia in women with mildly abnormal cervical cytology. International Journal of Cancer. 2013;133(3).

    Article  CAS  Google Scholar 

  23. Mirabello L, Schiffman M, Ghosh A, Rodriguez AC, Vasiljevic N, Wentzensen N, et al. Elevated methylation of HPV16 DNA is associated with the development of high grade cervical intraepithelial neoplasia. International Journal of Cancer. 2013;132(6).

    Article  Google Scholar 

  24. Lorincz AT. Cancer diagnostic classifiers based on quantitative DNA methylation. Expert review of molecular diagnostics. 2014;14(3).

    Article  CAS  Google Scholar 

  25. Denis H, Ndlovu MN, Fuks F. Regulation of mammalian DNA methyltransferases: a route to new mechanisms. EMBO Rep. 2011;12(7).

    Article  CAS  Google Scholar 

  26. Louvanto K, Franco EL, Ramanakumar AV, Vasiljevi N, Scibior-Bentkowska D, Koushik A, et al. Methylation of viral and host genes and severity of cervical lesions associated with human papillomavirus type 16. International Journal of Cancer. 2015;136:6.

    Article  Google Scholar 

  27. Vasiljević N, Scibior-Bentkowska D, Brentnall A, Cuzick J, Lorincz A. A comparison of methylation levels in HPV18, HPV31 and HPV33 genomes reveals similar associations with cervical precancers. Journal of Clinical Virology. 2014;59(3).

    Article  Google Scholar 

  28. Cook DA, Krajden M, Brentnall AR, Gondara L, Chan T, Law JH, et al. Evaluation of a validated methylation triage signature for human papillomavirus positive women in the HPV FOCAL cervical cancer screening trial. International Journal of Cancer. 2018.

  29. Nedjai B, Reuter C, Ahmad A, Banwait R, Warman R, Carton J, et al. Molecular progression to cervical precancer, epigenetic switch or sequential model? International Journal of Cancer. 2018;143(7).

    Article  CAS  Google Scholar 

  30. Brentnall AR, Vasiljevic̈ N, Scibior-Bentkowska D, Cadman L, Austin J, Szarewski A, et al. A DNA methylation classifier of cervical precancer based on human papillomavirus and human genes. International Journal of Cancer. 2014;135(6).

    Article  CAS  Google Scholar 

  31. Brentnall AR, Vasiljevic N, Scibior-Bentkowska D, Cadman L, Austin J, Cuzick J, et al. HPV33 DNA methylation measurement improves cervical pre-cancer risk estimation of an HPV16, HPV18, HPV31 and EPB41L3 methylation classifier. Cancer Biomarkers. 2015;15(5).

    Article  CAS  Google Scholar 

  32. Lorincz AT, Brentnall AR, Scibior-Bentkowska D, Reuter C, Banwait R, Cadman L, et al. Validation of a DNA methylation HPV triage classifier in a screening sample. International Journal of Cancer. 2016;138(11).

    Article  CAS  Google Scholar 

  33. Van BR, Van Der MJ, Snijders PJF, Rodriquez-manfredi A, Harmsel B, Van Den MHAM, et al. CADM1 and MAL methylation status in cervical scrapes is representative of the most severe underlying lesion in women with multiple cervical biopsies. International Journal of Cancer. 2016;138.

  34. How Kit A, Nielsen HM, Tost J. DNA methylation based biomarkers: practical considerations and applications. Biochimie. 2012;94(11).

    Article  CAS  Google Scholar 

  35. Clarke MA, Luhn P, Gage JC, Bodelon C, Dunn ST, Walker J, et al. Discovery and validation of candidate host DNA methylation markers for detection of cervical precancer and cancer. International Journal of Cancer. 2017;710.

  36. Clarke MA, Gradissimo A, Schiffman M, Lam J, Sollecito CC, Fetterman B, et al. Human papillomavirus DNA methylation as a biomarker for cervical precancer: consistency across 12 genotypes and potential impact on management of HPV-positive women. Clinical Cancer Research. 2018;24(9).

    Article  CAS  Google Scholar 

  37. Schmitz M, Eichelkraut K, Schmidt D, Zeiser I, Hilal Z, Tettenborn Z, et al. Performance of a DNA methylation marker panel using liquid-based cervical scrapes to detect cervical cancer and its precancerous stages. BMC Cancer. 2018;18(1).

  38. Floore A, Hesselink A, Oštrbenk A, Alcaniz E, Rothe B, Pedersen H, et al. Intra- and inter-laboratory agreement of the FAM19A4/mir124-2 methylation test: results from an international study. 2019;33(e22854).

  39. Schmitz M, Wunsch K, Hoyer H, Scheungraber C, Runnebaum IB, Hansel A, et al. Performance of a methylation specific real-time PCR assay as a triage test for HPV-positive women. Clinical epigenetics. 2017;9.

  40. Vasiljević N, Scibior-Bentkowska D, Brentnall AR, Cuzick J, Lorincz AT. Credentialing of DNA methylation assays for human genes as diagnostic biomarkers of cervical intraepithelial neoplasia in high-risk HPV positive women. Gynecologic Oncology. 2014;132(3).

    Article  Google Scholar 

  41. Lorincz AT. Virtues and weaknesses of DNA methylation as a test for cervical cancer. Acta Cytologica. 2016.

  42. Lazcano-Ponce E, Lorincz A, Torres L, Salmerón J, Cruz A, Rojas R. Specimen self-collection and HPV DNA screening in a pilot study of 100,242 women. International Journal of Cancer. 2014;116.

  43. Beal CM, Salmerón J, Flores YN, Torres L, Granados-garcía V, Dugan E, et al. Cost analysis of different cervical cancer screening strategies in Mexico. Salud Publica de Mexico. 2014;56(5).

    Article  Google Scholar 

  44. de Sanjose S, Quint WG, Alemany L, Geraets DT, Klaustermeier JE, Lloveras B, et al. Human papillomavirus genotype attribution in invasive cervical cancer: a retrospective cross-sectional worldwide study. The lancet oncology [Internet]. 2010 Nov [cited 2014 Jul 23];11(11). Available from:

  45. Wentzensen N, Sun C, Ghosh A, Kinney W, Mirabello L, Wacholder S, et al. Methylation of HPV18, HPV31, and HPV45 genomes and cervical intraepithelial neoplasia grade 3. J Natl Cancer Inst. 2012;104(22).

    Article  CAS  Google Scholar 

  46. Snijders P, Cubie H, Lorincz A, Ogilvie G, Franceschi S, Arbyn M, et al. Eurogin roadmap 2017: triage strategies for the management of HPV-positive women in cervical screening programs. International Journal of Cancer. 2018;143(4).

  47. Torres-Ibarra L, Lazcano-Ponce E, Franco EL, Cuzick J, Hernández-Ávila M, Lorincz A, et al. Triage strategies in cervical cancer detection in Mexico: methods of the FRIDA study. Salud Publica de Mexico. 2016;58(2).

    Article  Google Scholar 

Download references


We would like to thank all of the health staff in Tlaxcala for their enthusiastic work as well as the HPV laboratory at INSP and the molecular laboratory at QMUL for valuable support. The authors would also like to extend special acknowledgments to the FRIDA study population for their generous participation as well as the Tlaxcala FRIDA field team for collecting the data and supervising field procedures. The study sponsors had no role in the study design, in the collection, analysis, and interpretation of data, in the writing of the manuscript, or in the decision to submit the manuscript for publication.

*FRIDA Study Group: Attila Lorincz, Cosette Wheeler, Patti Gravitt, Eduardo Lazcano, Leticia Torres, Leith León, Paula Ramírez, Berenice Rivera, Eduardo L. Franco, Jack Cuzick, Pablo Méndez, Jorge Salmerón, Mauricio Hernández, Anna Barbara Moscicki, Yvonne Flores, Enrique Carmona, Kathleen M. Schmeler, David Bishai, Pilar Hernández, Rubi Hernández, Indira Mendiola.


This study was supported by the Instituto Nacional de Salud Publica of México, the Coordinación de Investigación en Salud del Instituto Mexicano del Seguro Social, the Secretaria de Salud Tlaxcala, the Instituto Nacional de las Mujeres, and the Consejo Nacional de Ciencia y Tecnología. Additional support was provided by Cancer Research UK (grant number C569/A10404) for the research done at the Molecular Epidemiology Laboratory, Queen Mary University of London.

The study sponsors did not play any role in designing the study, collecting, analyzing, or interpreting the data, writing the report, or submitting this paper for publication.

Author information

Authors and Affiliations




RHL, AL, LTI, CR, DSB, and JS were responsible for the study design and had full access to all the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis. JS, BR, PR, IMP, LLM, ELP, and LTI were the site investigators; they contributed to study implementation and study coordination. BR, PR, LLM, RHL, RW, and JS were responsible for the data collection and study supervision. RHL, AL, LTI, CR, DSB, BN, LLM, ELP, JC, and JS were responsible for drafting the manuscript and contributed to the data interpretation. RHL, AL, PR, LTI, and JC performed and reviewed the statistical analyses. All authors provided critical revision of content and approved the final version before submission.

Corresponding author

Correspondence to Attila T. Lorincz.

Ethics declarations

Ethics approval and consent to participate

The FRIDA study was reviewed and approved by the research and ethics committees of the Instituto Nacional de Salud Pública (INSP) (1094), Mexican regulatory agency COFEPRIS (CAS/OR/01/CAS/123300410C0044-3578/2012), and Tlaxcala State Ministry of Health (SSET) (SS.DECI-OI-13/12). The purpose and procedures of the study were explained to all participants, and informed verbal consent was obtained. The FRIDA Study design and management were supported by a group of principal investigators, steering committee, and subcommittees including molecular triage, quality control of cytology and pathology, and safety and case management.

Competing interests

JS and ELP have received funding for research projects from Qiagen, Roche, Beckton Dickinson, Gen-Probe, DICIPA, and Arbor Vita. The other authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Hernández-López, R., Lorincz, A.T., Torres-Ibarra, L. et al. Methylation estimates the risk of precancer in HPV-infected women with discrepant results between cytology and HPV16/18 genotyping. Clin Epigenet 11, 140 (2019).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: