Background:

Cervical cancer screening with high-risk human papillomavirus (HrHPV) testing is being introduced. Most HrHPV infections are transient, requiring triage tests to identify individuals at highest risk for progression to cervical cancer. Head-to-head comparisons of available strategies for screening and triage are needed. Endometrial and ovarian cancers could be amenable to similar testing.

Methods:

Between 2016 and 2020, discarded cervical cancer screening specimens from women ages 25 to 65 undergoing screening at Kaiser Permanente Northern California were collected. Specimens were aliquoted, stabilized, and stored frozen. Human papillomavirus (HPV), cytology, and histopathology results as well as demographic and cofactor information were obtained from electronic medical records (EMR). Follow-up collection of specimens was conducted for 2 years, and EMR-based data collection was planned for 5 years.

Results:

Collection of enrollment and follow-up specimens is complete, and EMR-based follow-up data collection is ongoing. At baseline, specimens were collected from 54,957 HPV-positive, 10,215 HPV-negative/Pap-positive, and 12,748 HPV-negative/Pap-negative women. Clinical history prior to baseline was available for 72.6% of individuals, of which 53.9% were undergoing routine screening, 8.6% recently had an abnormal screen, 30.3% had previous colposcopy, and 7.2% had previous treatment. As of February 2021, 55.7% had one or more colposcopies, yielding 5,563 cervical intraepithelial neoplasia grade 2 (CIN2), 2,756 cervical intraepithelial neoplasia grade 3 (CIN3), and 146 cancer histopathology diagnoses.

Conclusions:

This robust population-based cohort study represents all stages of cervical cancer screening, management, and posttreatment follow-up.

Impact:

The IRIS study is a unique and highly relevant resource allowing for natural history studies and rigorous evaluation of candidate HrHPV screening and triage markers, while permitting studies of biomarkers associated with other gynecologic cancers.

The recognition that persistent infections with carcinogenic HPV were a necessary cause of almost all cervical cancers (1–3) and its immediate precursor abnormalities (precancer) has led to introduction of human papillomavirus (HPV) testing into cervical cancer screening (4–7). Across primary HPV screening strategies the question is how best to manage HPV-positive individuals, most of whom have transient infections that will not progress to cancer. We cannot yet accurately predict which individuals testing HPV-positive will fail to control their infections and eventually develop precancer and, if left untreated, possibly cancer. Yet, it is not feasible to send all people testing HPV-positive to colposcopy-guided biopsy to search for precancer. Therefore, screening strategies that utilize high-risk HPV (HrHPV) tests require additional triage tests to determine who among those screening HrHPV-positive need to go to colposcopy. Clinical scenarios postcolposcopy and posttreatment also require management decisions that incorporate HrHPV test and triage.

In primary HPV screening, currently recommended triage strategies for individuals testing HPV-positive include HPV16/18 genotyping and/or cytology, approaches with limitations including lower sensitivity of cytology and the reduced utility of HPV16/18 typing in HPV-vaccinated populations (8, 9). New robust and efficient biomarkers are required to triage HPV-positive women. Cytology-based p16/Ki-67 dual stain is reproducible with higher sensitivity and specificity than cytology, resulting in improved risk stratification that reduces colposcopy referral while finding more precancer (10). Because it provides excellent long-term risk stratification compared with cytology, dual stain–negative, HPV-positive women might be released to extended surveillance intervals (11). An automated algorithm using machine learning for cytology can match manual performance and provide a quantitative result, allowing for evaluation of dual stain at higher thresholds that confer a risk similar to high grade squamous intraepithelial lesion (HSIL) cytology (12). Another promising biomarker is HPV and host DNA methylation which is associated with increased disease severity (13–15). The performance of these novel biomarkers in clinical scenarios such as postcolposcopy and posttreatment has not been established.

In existing management guidelines, HPV16 is considered uniquely carcinogenic and HPV18 and HPV45 are associated with higher risks, especially for adenocarcinomas; all warranting immediate colposcopy. Differentiating between other HPV16-related types HPV31, HPV33, HPV35, HPV52, and HPV58 that pose substantially higher risk than the lowest-risk group consisting of HPV39, HPV51, HPV56, HPV59, and HPV68 may provide additional risk stratification (16).

The abundance of new test options demands well-powered head—to–head evaluations in a well-defined screening population to evaluate their performance relative to each other and to identify optimal management strategies. Such evaluations are required for screening and management guidelines that provide specific recommendations on how new technologies should be incorporated in clinical management (6, 8). To determine appropriate recommendations for a given test result, absolute risks of cervical intraepithelial neoplasia grade 3 or worse (CIN3+) are required at multiple follow-up time points (17, 18). Beyond screening, novel biomarkers have the potential to improve management of individuals who have abnormal screening results that do not trigger immediate colposcopy referral, individuals with colposcopy findings that do not require treatment, and individuals previously treated (8).

To address this, we created a biorepository of cervical specimens from a large population to evaluate the performance of biomarkers in cervical cancer screening and management while also allowing for evaluation of the impact of HPV vaccination on screening and management performance. Cervical biobank studies storing residual liquid-based cytology and HPV testing samples are a robust resource for natural history studies and longitudinal analyses of novel biomarkers (19–21). The Improving Risk Informed HPV Screening (IRIS) study includes stored cervical specimens from screening and follow-up linked with relevant cofactors, clinical history, and follow-up information. Specimens from the IRIS cohort will enable well-powered etiologic analyses of host and viral somatic changes, HPV genome variation, association of host ancestry, cervicovaginal microbiome, and HPV cofactors with the progression from HPV infection to precancer. In addition, the IRIS study will allow for the evaluation of biomarkers associated with other gynecologic cancers, including endometrial and ovarian cancers in prospectively collected cervical specimens from a large population.

Study population

Kaiser Permanente Northern California (KPNC) is a large integrated healthcare delivery system providing cervical cancer screening to individuals age 25 and older with concurrent HPV and Pap testing (“cotesting”) every 3 years. The KPNC Regional Laboratory processes almost 300,000 pairs of cytology and HPV tests each year. The representative collection of baseline specimens occurred between 2016 and 2018 from a consecutive sample of individuals age 25 to 65 undergoing cotesting at KPNC. After baseline discard specimens were collected, members were mailed an opt-out consent letter informing about the study and providing the opportunity to opt-out of participation. Women could respond by phoning a toll-free number or returning a prestamped, preaddressed letter. Specimens for women who chose to opt out or had returned mail were destroyed and their specimen test results and personal health information are not used and removed from the study dataset.

Specimens were selected according to cotest result, with an objective of sampling almost all HPV-positive individuals and a random sample of individuals testing negative for both HPV and cytology, to allow weighting results back to the full KPNC population (Table 1). The time periods of collection varied by sampling strata: HPV-positive, HPV-negative/cytology-positive, and HPV-negative/cytology-negative (referred to as “cotest negative”). As presented in Table 1, HPV-positive and HPV-negative/Pap-positive specimens were collected during the entire period (April 2016–March 2018 and January 2016–August 2018) while cotest-negative specimens were sampled early in the study to maximize follow-up time during the study period (March 2016–January 2017).

Table 1.

IRIS cohort to full KPNC population age 25+ screened between January 2016 and August 2018.

Baseline cotest resultUnique individuals with one or more cotest results between January 2016 and August 2018aSampling time periods for each cotest groupUnique individuals with cotest result during sampling time periodNumber of individuals sampled (percent of unique individuals with cotest during sampling period)Opt-out received (percent of individuals that participated)Final individuals in IRIS cohort
Total 897,680  305,083 81,348 (26.7%) 3,428 (4.2%) 77,920 
HPV/Papb 809,199 Mar 2016–Jan 2017c 228,749 13,422 (5.9%) 674 (5.0%) 12,748 
HPV/Pap+b 13,512 Jan 2016–Aug 2018 13,437 10,791 (80.3%) 576 (5.3%) 10,215 
HPV+b 74,969 Apr 2016–Mar 2018 62,897 57,135 (90.8%) 2,178 (3.8%) 54,957 
Baseline cotest resultUnique individuals with one or more cotest results between January 2016 and August 2018aSampling time periods for each cotest groupUnique individuals with cotest result during sampling time periodNumber of individuals sampled (percent of unique individuals with cotest during sampling period)Opt-out received (percent of individuals that participated)Final individuals in IRIS cohort
Total 897,680  305,083 81,348 (26.7%) 3,428 (4.2%) 77,920 
HPV/Papb 809,199 Mar 2016–Jan 2017c 228,749 13,422 (5.9%) 674 (5.0%) 12,748 
HPV/Pap+b 13,512 Jan 2016–Aug 2018 13,437 10,791 (80.3%) 576 (5.3%) 10,215 
HPV+b 74,969 Apr 2016–Mar 2018 62,897 57,135 (90.8%) 2,178 (3.8%) 54,957 

aNumber used for weighting to full KPNC population with one or more cotests between January 2016 and August 2018.

bResults for the first cotest collected between January 2016 and August 2018 are assigned to each person screened. HPV/Pap+ includes HPV-negative/ASC-US cytology. Table excludes 29 individuals with missing HPV result and 365 individuals with missing Pap result. We also exclude 2,195 women with samples collected outside of the sampling time window and 670 women who were <25 years of age.

cHPV/Pap was not sampled during July/August 2016.

The details of specimen collection, processing, aliquoting, shipping, and storage are described separately (22). Briefly, Surepath (Becton Dickinson) and standard transport medium (STM; Qiagen) specimens collected for cytology and HPV testing were stabilized as needed and transferred into storage tubes as previously described (21). Immediate cryopreservation was not possible in the clinical setting, but specimens were frozen as soon as possible, allowing for compatibility with future analysis of RNA and DNA protein that would allow a range of assays from RNA expression, methylation, and microbiome. They were then shipped to the NCI Biorepository for long-term storage. From a stratified random sample of HrHPV-negative and -positive specimens, a total of 48,794 liquid-based cytology slides were created and stained with p16/ki-67 dual stain between June 2016 and March 2018 (22).

Follow-up HPV and cytology specimens for IRIS individuals were collected through November 2020. Deidentified HPV, cytology, and histopathology results as well as demographic and cofactor information for individuals were obtained from electronic medical records (EMR) for baseline and will continue to be collected for an additional 5 years or more. The IRIS study is an NCI-funded public resource. Study data will be made available for research purposes in accordance with institutional review boards upon request. The number of unique individuals in the KPNC screening program with cotest results during the study period was also collected. The study was approved by ethical review boards at Kaiser Permanente and the NCI.

Analysis

The number of follow-up cotest specimens collected as of November 2020 were summarized by time of baseline collection and baseline cotest result. Each individual's status in the cervical cancer screening and management process was classified hierarchically using their clinical history prior to the baseline cotest. Those with no previous information available were categorized as an “unknown clinical history.” Those with any history of excisional treatment were considered “previously treated.” Individuals with a record of colposcopy prior to baseline were categorized as “postcolposcopy” and further classified as “recent” or “distant” depending on whether the colposcopy was conducted within 42 months prior to baseline cotest. Those whose previous screening result was abnormal [HPV-positive or cytology low-grade squamous intraepithelial lesion (LSIL) or worse] were categorized as “previous screen abnormal” while those whose previous screen was normal [HPV-negative and cytology atypical squamous cells of undetermined significance (ASC-US) or negative) were categorized as “routine screening.” Based upon the 3-year screening intervals at KPNC, routine screening was divided by time of screening visit into “early” (<=30 months), “timely” (over 30 and less than or equal to 42 months), or “late” (>42 months). This clinical history among women participating in the IRIS study was compared with the entire KPNC screening population. Baseline cotest results and follow-up activities were summarized by categories of history of previous screening and treatment. Follow-up activities included repeat cotests, colposcopy visits, and histology results from biopsy and excisional procedures. The current membership in status was also summarized.

Statistical considerations

The planned number of individuals enrolled in IRIS was 50,000 HPV-positive individuals, 10,000 HPV-negative/cytology-positive individuals, and 10,000 women negative for both HPV and cytology with the expectation that approximately 2,500, 50, and 5 precancers (CIN3+) would be detected in each group, respectively, most of which would be diagnosed within 3 to 4 years after baseline.

This sample size was based upon calculations for statistical power to evaluate novel biomarkers for primary screening and triage of HPV-positive results (23). The required precision depended upon the anticipated performance characteristics of individual risk markers or combined strategies. The precision of absolute risk estimates would be high overall, for example for evaluation of primary screening, a 3-year CIN3 risk of 5.5%, 1.4%, and 0.45% would have a SE of 0.22%, 0.16%, and 0.06%, respectively. For triage, an immediate CIN3 risk of 19.7% and 1.4% would have a SE of 0.7% and 0.11%, respectively.

Between January 2016 and August 2018, 897,680 KPNC individuals age 25 and older had a cotest at KPNC (Table 1). Baseline specimens from individuals with cotest negative, HPV-positive/cytology-positive, and HPV-positive results were collected at different, overlapping time periods. During these time periods, 305,083 individuals had a cotest across the three strata (Table 1). Overall, the IRIS study collected baseline specimens for 81,348 individuals– representing 26.7% of the 305,083 individuals with a cotest during the sampling time periods. The proportion of the individuals in the sampling strata with a specimen collected varied by cotest result: 5.9% of women testing cotest-negative, 80.3% of women testing HPV-negative/Pap-positive, and 90.8% of women testing HPV-positive were sampled, respectively. Opt-out requests were received for 3,428 (4.2%) individuals and those with HPV-positive baseline results were less likely to opt out than other sampling groups (3.8% vs. 5.3% HPV-negative/Pap-positive, and 5.0% HPV-positive; P < 0.001). The final analytic cohort included 77,920 individuals.

Demographic characteristics of individuals are presented in Table 2. Those with unknown screening history were younger and therefore less likely to have been a KPNC member for 5 years or more (mean age of individuals with 5 or more years membership was 42.2 years). Individuals with previous abnormal screen were also more likely to be younger (likely because they had not had previous colposcopy or treatment).

Table 2.

Demographic characteristics of individuals in IRIS cohort by history of screening and treatment.

Routine screeningPrevious screen abnormalPostcolposcopyPosttreatmentUnknown screening history
n%n%n%n%n%
Total30,5201004,85410017,1301004,09510021,321100
Age 
 Mean (SD) 40.7 (11.7)  36.1 (11.0)  41.0 (11.7)  44.2 (12.3)  35.6 (10.5)  
 Median (IQR) 39 (30–50)  32 (28–42)  38 (31–50)  42 (34–54)  32 (27–41)  
 25–29 6,662 21.8 1,732 35.7 3,082 18.0 414 10.1 8,341 39.1 
 30–39 9,183 30.1 1,739 35.8 6,211 36.3 1,355 33.1 6,958 32.6 
 40–49 6,740 22.1 654 13.5 3,419 20.0 946 23.1 3,088 14.5 
 50–59 5,495 18.0 457 9.4 2,873 16.8 804 19.6 2,143 10.1 
 60–65 2,289 7.5 228 4.7 1,169 6.8 350 8.6 718 3.4 
 66+ 151 0.5 44 0.9 376 2.2 226 5.5 73 0.3 
Race/ethnicity 
 White 13,512 44.3 2,143 44.2 7,653 44.7 1,928 47.1 9,061 42.5 
 Hispanic 7,190 23.6 1,082 22.3 3,988 23.3 876 21.4 4,610 21.6 
 African American 2,380 7.8 484 10.0 1,453 8.5 354 8.6 1,265 5.9 
 Asian/Pacific islander 5,624 18.4 833 17.2 3,071 17.9 728 17.8 4,631 21.7 
 Multiracial/other 1,814 5.9 312 6.4 965 5.6 209 5.1 1,754 8.2 
Month of previous membership 
 <2 years 1,545 5.1 1,020 21.0 1,334 7.8 216 5.3 18,100 84.9 
 2–4.99 years 7,740 25.4 1,372 28.3 3,859 22.5 702 17.1 2,233 10.5 
 5–7.99 years 6,596 21.6 880 18.1 3,636 21.2 847 20.7 416 2.0 
 8+ years 14,593 47.8 1,569 32.3 8,277 48.3 2,327 56.8 327 1.5 
 Missing 46 0.2 13 0.3 24 0.1 0.1 245 1.1 
Body mass index 
 Under/normal weight 11,006 36.1 2,006 41.3 6,515 38.0 1,494 36.5 9,643 45.2 
 Overweight 8,120 26.6 1,194 24.6 4,559 26.6 1,154 28.2 5,094 23.9 
 Obese 6,690 21.9 1,021 21.0 3,688 21.5 918 22.4 3,811 17.9 
 Very obese 1,887 6.2 278 5.7 972 5.7 225 5.5 966 4.5 
 Unknown 2,817 9.2 355 7.3 1,396 8.2 304 7.4 1,807 8.5 
Routine screeningPrevious screen abnormalPostcolposcopyPosttreatmentUnknown screening history
n%n%n%n%n%
Total30,5201004,85410017,1301004,09510021,321100
Age 
 Mean (SD) 40.7 (11.7)  36.1 (11.0)  41.0 (11.7)  44.2 (12.3)  35.6 (10.5)  
 Median (IQR) 39 (30–50)  32 (28–42)  38 (31–50)  42 (34–54)  32 (27–41)  
 25–29 6,662 21.8 1,732 35.7 3,082 18.0 414 10.1 8,341 39.1 
 30–39 9,183 30.1 1,739 35.8 6,211 36.3 1,355 33.1 6,958 32.6 
 40–49 6,740 22.1 654 13.5 3,419 20.0 946 23.1 3,088 14.5 
 50–59 5,495 18.0 457 9.4 2,873 16.8 804 19.6 2,143 10.1 
 60–65 2,289 7.5 228 4.7 1,169 6.8 350 8.6 718 3.4 
 66+ 151 0.5 44 0.9 376 2.2 226 5.5 73 0.3 
Race/ethnicity 
 White 13,512 44.3 2,143 44.2 7,653 44.7 1,928 47.1 9,061 42.5 
 Hispanic 7,190 23.6 1,082 22.3 3,988 23.3 876 21.4 4,610 21.6 
 African American 2,380 7.8 484 10.0 1,453 8.5 354 8.6 1,265 5.9 
 Asian/Pacific islander 5,624 18.4 833 17.2 3,071 17.9 728 17.8 4,631 21.7 
 Multiracial/other 1,814 5.9 312 6.4 965 5.6 209 5.1 1,754 8.2 
Month of previous membership 
 <2 years 1,545 5.1 1,020 21.0 1,334 7.8 216 5.3 18,100 84.9 
 2–4.99 years 7,740 25.4 1,372 28.3 3,859 22.5 702 17.1 2,233 10.5 
 5–7.99 years 6,596 21.6 880 18.1 3,636 21.2 847 20.7 416 2.0 
 8+ years 14,593 47.8 1,569 32.3 8,277 48.3 2,327 56.8 327 1.5 
 Missing 46 0.2 13 0.3 24 0.1 0.1 245 1.1 
Body mass index 
 Under/normal weight 11,006 36.1 2,006 41.3 6,515 38.0 1,494 36.5 9,643 45.2 
 Overweight 8,120 26.6 1,194 24.6 4,559 26.6 1,154 28.2 5,094 23.9 
 Obese 6,690 21.9 1,021 21.0 3,688 21.5 918 22.4 3,811 17.9 
 Very obese 1,887 6.2 278 5.7 972 5.7 225 5.5 966 4.5 
 Unknown 2,817 9.2 355 7.3 1,396 8.2 304 7.4 1,807 8.5 

The IRIS cohort consisted of women undergoing cotesting for a variety of indications (Table 3): about one quarter (27.4%) had no previous clinical history at KPNC. 39.2% were routine screening visits, 6.2% had a previous abnormal screening test, 22.0% had a previous colposcopy, and 5.3% had been previously treated. Most routine screening visits were conducted around 3 years (30–42 months) after their previous screen. IRIS individuals with previous colposcopy were more likely to have had the procedure in the past 42 months versus later (14.4% vs. 7.6% of total cohort). According to the study design, individuals selected for IRIS were more likely to have previous treatment, colposcopy, and abnormal screening compared with the overall population. By study design, individuals with abnormal test results were more likely to be selected for IRIS and therefore they were less likely to be in routine screening than the overall population with a cotest at KPNC (39.2% vs. 58.4%). The 30,520 individuals in routine screening represent approximately 58.4% of the full routine screening population.

Table 3.

History of screening and treatment by baseline cotest result among individuals in IRIS cohort.

Baseline cotest result
High-grade+LSILASC-USNILM
Participant clinical historyaIRIS cohort% of IRIS cohort% of KPNC populationa,bHPV+HPVHPV+HPVHPV+HPVHPV+HPV
Total 77,920 100% 100% 3,449 1,031 11,678 1,272 15,984 7,912 23,846 12,748 
Routine screening 30,520 39.2% 58.4% 722 427 3,418 506 4,824 3,894 8,964 7,765 
(previous screen negative)    20.9% 41.4% 29.3% 39.8% 30.2% 49.2% 37.6% 60.9% 
 Early (≤30 months) 4,473 5.7% 14.6% 87 54 600 66 782 568 1,371 945 
    2.5% 5.2% 5.1% 5.2% 4.9% 7.2% 5.7% 7.4% 
 Timely (≤42 months) 17,630 22.6% 37.7% 349 267 1,801 316 2,646 2,316 4,941 4,994 
    10.1% 25.9% 15.4% 24.8% 16.6% 29.3% 20.7% 39.2% 
 Late (>42 months) 8,417 10.8% 7.1% 286 106 1,017 124 1,396 1,010 2,652 1,826 
    8.3% 10.3% 8.7% 9.7% 8.7% 12.8% 11.1% 14.3% 
Previous screen abnormal 4,854 6.2% 3.3% 258 16 812 56 1,227 229 2,052 204 
    7.5% 1.6% 7.0% 4.4% 7.7% 2.9% 8.6% 1.6% 
Postcolposcopy 17,130 22.0% 10.3% 892 245 3,302 305 4,465 1,551 4,985 1,385 
    25.9% 23.8% 28.3% 24.0% 27.9% 19.6% 20.9% 10.9% 
 Recent (<42 months) 11,225 14.4% 3.2% 659 151 2,545 182 3,191 748 3,255 494 
    19.1% 14.6% 21.8% 14.3% 20.0% 9.5% 13.7% 3.9% 
 Distant (42+ months) 5,905 7.6% 7.1% 233 94 757 123 1,274 803 1,730 891 
    6.8% 9.1% 6.5% 9.7% 8.0% 10.1% 7.3% 7.0% 
Posttreatment 4,095 5.3% 2.5% 392 130 631 83 807 453 1,133 466 
    11.4% 12.6% 5.4% 6.5% 5.0% 5.7% 4.8% 3.7% 
Unknown clinical history 21,321 27.4% 25.5% 1,185 213 3,515 322 4,661 1,785 6,712 2,928 
    34.4% 20.7% 30.1% 25.3% 29.2% 22.6% 28.1% 23.0% 
Baseline cotest result
High-grade+LSILASC-USNILM
Participant clinical historyaIRIS cohort% of IRIS cohort% of KPNC populationa,bHPV+HPVHPV+HPVHPV+HPVHPV+HPV
Total 77,920 100% 100% 3,449 1,031 11,678 1,272 15,984 7,912 23,846 12,748 
Routine screening 30,520 39.2% 58.4% 722 427 3,418 506 4,824 3,894 8,964 7,765 
(previous screen negative)    20.9% 41.4% 29.3% 39.8% 30.2% 49.2% 37.6% 60.9% 
 Early (≤30 months) 4,473 5.7% 14.6% 87 54 600 66 782 568 1,371 945 
    2.5% 5.2% 5.1% 5.2% 4.9% 7.2% 5.7% 7.4% 
 Timely (≤42 months) 17,630 22.6% 37.7% 349 267 1,801 316 2,646 2,316 4,941 4,994 
    10.1% 25.9% 15.4% 24.8% 16.6% 29.3% 20.7% 39.2% 
 Late (>42 months) 8,417 10.8% 7.1% 286 106 1,017 124 1,396 1,010 2,652 1,826 
    8.3% 10.3% 8.7% 9.7% 8.7% 12.8% 11.1% 14.3% 
Previous screen abnormal 4,854 6.2% 3.3% 258 16 812 56 1,227 229 2,052 204 
    7.5% 1.6% 7.0% 4.4% 7.7% 2.9% 8.6% 1.6% 
Postcolposcopy 17,130 22.0% 10.3% 892 245 3,302 305 4,465 1,551 4,985 1,385 
    25.9% 23.8% 28.3% 24.0% 27.9% 19.6% 20.9% 10.9% 
 Recent (<42 months) 11,225 14.4% 3.2% 659 151 2,545 182 3,191 748 3,255 494 
    19.1% 14.6% 21.8% 14.3% 20.0% 9.5% 13.7% 3.9% 
 Distant (42+ months) 5,905 7.6% 7.1% 233 94 757 123 1,274 803 1,730 891 
    6.8% 9.1% 6.5% 9.7% 8.0% 10.1% 7.3% 7.0% 
Posttreatment 4,095 5.3% 2.5% 392 130 631 83 807 453 1,133 466 
    11.4% 12.6% 5.4% 6.5% 5.0% 5.7% 4.8% 3.7% 
Unknown clinical history 21,321 27.4% 25.5% 1,185 213 3,515 322 4,661 1,785 6,712 2,928 
    34.4% 20.7% 30.1% 25.3% 29.2% 22.6% 28.1% 23.0% 

Note: High-grade+ includes atypical glandular cells of undetermined significance (AGUS), atypical squamous cells, cannot rule out high-grade squamous intraepithelial lesion (ASC-H), HSILs, and carcinoma.

Abbreviation: NILM, negative for intraepithelial lesion or malignancy.

aClinical history at first cotest during period between January 2016 and August 2018.

bDistribution at KPNC for individuals with one or more cotest results between January 2016 and August 2018.

Baseline cotest results varied by clinical history with negative baseline screening results more common among women in routine screening versus previous colposcopy or treatment (Table 3). Individuals with a high-grade baseline cytology result were more likely to have no prior screening record or previous treatment. Due to the sampling scheme, IRIS oversampled individuals with HPV-negative/HSIL or worse cytology results and they were more likely to have been previously treated.

Follow-up activities stratified by category of previous clinical history are presented in Table 4. As of February 2020, 75.3% of IRIS individuals were still active KPNC individuals. Individuals with unknown clinical history were less likely to be active members (63.9%), compared with those who had previous treatment (83.3%) or colposcopy (80.5%). As of November 2020, most individuals (83.7%) had either a follow-up cotest and/or a colposcopy. Specifically, 75.9% had a follow-up cotest (99.0% of which were collected for IRIS) and 55.7% had colposcopy either at baseline or during follow-up. These procedures were more common when the individual had previous abnormal screening results, colposcopy, or treatment. Through February 2021, 5,563 cervical intraepithelial neoplasia grade 2 (CIN2), 2,756 CIN3 (including 219 adenocarcinoma in situ), and 146 cervical cancers (including 69 adenocarcinoma) were identified among IRIS individuals so far, with passive follow-up continuing. Among individuals with CIN2+ histologic diagnoses, CIN2 was more common among individuals undergoing timely or early routine screening (226 of 303, 74.6%) and (835 of 1,095, 76.3%) while individuals with previous treatment had a greater proportion of CIN3 (346 of 806, 42.9%) and cancer (27 of 807, 3.3%) diagnoses.

Table 4.

History of previous screening and treatment and follow-up among individuals in IRIS study.

Active KPNC membersbNumber (%) with follow-up cotestscHistology resultsd
Participant's clinical historyaIRIS cohortN (%)12 or moreNumber (%) with 1 or more colposcopy visitsdNumber with cotest, colpo, or histology follow-upCIN2CIN3Cancer
Total 77,920 58,659 35,498 23,610 43,401 65,199 5,563 2,756 146 
  (75.3%) (45.6%) (30.3%) (55.7%) (83.7%)    
Routine screening 30,520 24,372 15,336 7,710 12,571 24,470 1,520 560 23 
(previous screen negative)  (79.9%) (50.3%) (25.3%) (41.2%) (80.2%)    
 Early (≤30 months) 4,473 3,419 2,118 1,391 2,033 3,735 226 72 
  (76.4%) (47.4%) (31.1%) (45.5%) (83.5%)    
 Timely (≤42 months) 17,630 14,589 9,335 4,501 7,032 14,487 835 251 
  (82.8%) (53.0%) (25.5%) (39.9%) (82.2%)    
 Late (>42 months) 8,417 6,364 3,883 1,818 3,506 6,248 459 237 
  (75.6%) (46.1%) (21.6%) (41.7%) (74.2%)    
 Previous screen abnormal 4,854 3,462 1,955 1,966 3,946 4,453 495 272 
  (71.3%) (40.3%) (40.5%) (81.3%) (91.7%)    
 Postcolposcopy 17,130 13,793 7,389 6,887 12,659 15,717 1,558 663 31 
  (80.5%) (43.1%) (40.2%) (73.9%) (91.8%)    
 Recent (<42 months) 11,225 8,905 4,496 5,081 9,371 10,620 1,204 486 19 
  (79.3%) (40.1%) (45.3%) (83.5%) (94.6%)    
 Distant (42+ months) 5,905 4,888 2,893 1,806 3,288 5,097 354 177 12 
  (82.8%) (49.0%) (30.6%) (55.7%) (86.3%)    
Posttreatment 4,095 3,409 1,543 1,761 2,875 3,768 434 346 27 
  (83.3%) (37.7%) (43.0%) (70.2%) (92.0%)    
Unknown clinical history 21,321 13,623 9,275 5,286 11,350 16,791 1,556 915 56 
  (63.9%) (43.5%) (24.8%) (53.2%) (78.8%)    
Active KPNC membersbNumber (%) with follow-up cotestscHistology resultsd
Participant's clinical historyaIRIS cohortN (%)12 or moreNumber (%) with 1 or more colposcopy visitsdNumber with cotest, colpo, or histology follow-upCIN2CIN3Cancer
Total 77,920 58,659 35,498 23,610 43,401 65,199 5,563 2,756 146 
  (75.3%) (45.6%) (30.3%) (55.7%) (83.7%)    
Routine screening 30,520 24,372 15,336 7,710 12,571 24,470 1,520 560 23 
(previous screen negative)  (79.9%) (50.3%) (25.3%) (41.2%) (80.2%)    
 Early (≤30 months) 4,473 3,419 2,118 1,391 2,033 3,735 226 72 
  (76.4%) (47.4%) (31.1%) (45.5%) (83.5%)    
 Timely (≤42 months) 17,630 14,589 9,335 4,501 7,032 14,487 835 251 
  (82.8%) (53.0%) (25.5%) (39.9%) (82.2%)    
 Late (>42 months) 8,417 6,364 3,883 1,818 3,506 6,248 459 237 
  (75.6%) (46.1%) (21.6%) (41.7%) (74.2%)    
 Previous screen abnormal 4,854 3,462 1,955 1,966 3,946 4,453 495 272 
  (71.3%) (40.3%) (40.5%) (81.3%) (91.7%)    
 Postcolposcopy 17,130 13,793 7,389 6,887 12,659 15,717 1,558 663 31 
  (80.5%) (43.1%) (40.2%) (73.9%) (91.8%)    
 Recent (<42 months) 11,225 8,905 4,496 5,081 9,371 10,620 1,204 486 19 
  (79.3%) (40.1%) (45.3%) (83.5%) (94.6%)    
 Distant (42+ months) 5,905 4,888 2,893 1,806 3,288 5,097 354 177 12 
  (82.8%) (49.0%) (30.6%) (55.7%) (86.3%)    
Posttreatment 4,095 3,409 1,543 1,761 2,875 3,768 434 346 27 
  (83.3%) (37.7%) (43.0%) (70.2%) (92.0%)    
Unknown clinical history 21,321 13,623 9,275 5,286 11,350 16,791 1,556 915 56 
  (63.9%) (43.5%) (24.8%) (53.2%) (78.8%)    

aStatus at first cotest during period between January 2016 and August 2018.

bActive KPNC member as of February 2020.

cCotests through November 2020.

dColposcopy visits and histology results through February 2021.

With over 77,000 individuals included and 2,500+ precancers (CIN3) detected in the first 2 to 4 years of follow-up, the IRIS study is well-positioned to rapidly evaluate novel biomarkers in a variety of clinical scenarios to inform the changing landscape of cervical screening and management. The IRIS study will also permit studies of HPV carcinogenesis. To date, most of the women in the IRIS cohort remain KPNC members and follow-up information will be obtained from EMR records for at least 5 more years to identify outcomes as they occur over time during follow-up clinical care.

KPNC has a very large screening program that tested almost 900,000 members during a 20-month period between 2016 and 2018. By embedding the study within the regional laboratory of a large integrated health delivery system, baseline and follow-up discard specimens were systematically identified, retrieved, and processed, and stored without disrupting regular lab processes. Access to EMRs allowing for ascertainment of deidentified relevant clinical history, test results, as well as current and long-term outcomes creates additional efficiency for this research. Within the first 2 to 4 years of follow-up, our target numbers of precancer (CIN3) cases were reached for individuals with HPV-positive results at baseline (n = 2,701) and HPV-negative/Pap-positive results at baseline (n = 53). In the lowest-risk sampling stratum (HPV-negative/Pap-negative result at baseline), two of the anticipated five CIN3 cases had been detected. Case ascertainment is continuing in subsequent screening and management rounds. Ascertainment of outcomes among cotest-negatives is expected to take longer compared with HPV-positives, since clinical management is based on these results.

Because exact sampling fractions are available, analyses from the IRIS study can be weighted back to the entire KPNC screening population. The membership represents a well-screened population; their demographic characteristics align with those of the U.S. census-enumerated population in the Bay Area Metropolitan Statistical Area, with the exception of lowest and highest incomes (24). While the KPNC screening populations represents a high-resource setting, the risk estimates calculated within strata of cytology and HPV result from KPNC have shown to be portable to some other U.S. settings (18).

In the United States, the management of abnormal cervical cancer screening results are guided by the American Society for Colposcopy and Cervical Pathology (ASCCP) Consensus Guidelines. The 2019 ASCCP Risk-Based Consensus Management Guidelines were redesigned to current test results and history of cervical cancer screening tests and treatments for precancer. While KPNC represents an optimal scenario with regard to having prior history available, widely applicable estimates can be generated by collapsing across subsets of the population with and without known prior history. Importantly, the framework utilizes clinical history and risk calculations given test results to incorporate new technologies without the need for frequent interim guidance or full consensus conferences. Outside of IRIS, the full KPNC screening population is a unique large screening cohort that provides data for analyses pertaining to HPV testing, cytology, and other cofactors (17, 25–31).The KPNC screening population includes 1.5 million individuals aged 25 to 65 years routinely screened and provides immediate and 5-year risks of CIN3+ for combinations of current test results paired with history of screening test and colposcopy/biopsy result (17). Because the IRIS cohort is a subcohort within this larger effort, the results from IRIS will be directly applicable to the risk estimates that underlie current screening and management guidelines in the United States and will allow rapid update of guidelines when new technologies receive regulatory approval. As new biomarkers are developed, the large IRIS study is well poised to evaluate assays and provide risk calculations for many clinical scenarios, allowing direct head—to–head comparisons of screening, triage, and management tests which can accelerate moving new promising assays into clinical practice.

While the FDA continues to evaluate, approve, and regulate new technologies for management of abnormal cervical cancer screening results, the FDA frequently refers to clinical guidelines such as the 2019 ASCCP Consensus Guidelines, and its continuation, the Enduring Guidelines effort, to provide specific recommendations on how new technologies should be incorporated in patient management. Similarly, international guidelines are using risk-based approaches to evaluate new technologies and make practice recommendations (32). IRIS allows for evaluation of both common testing scenarios such as routine screening and unknown screening history as well as unique clinical scenarios such as postcolposcopy, posttreatment, and previous abnormal screening result in the absence of colposcopy. Using a risk-based approach, follow-up data from IRIS can be used to evaluate novel biomarkers to inform future clinical practice guidelines.

Future analyses will consider the performance of primary screening and triage markers among vaccinated individuals, using the KPNC vaccination registry, and simulated vaccinated cohorts based on models that include HPV typing. Future etiologic studies will evaluate HPV and host somatic changes, epigenetics, and the cervicovaginal microbiome among other markers (13, 33–37). Prior to IRIS, the HPV Persistence and Progression (PaP) Cohort was a very large “convenience sample” created between 2007 and 2010 by storing residual STM specimens from among approximately 55,000 mainly HPV-positive women who were cotested or triaged at KPNC. Three thousand, four hundred and ninety-nine women from the PaP Cohort were also enrolled in IRIS, allowing for even long-term natural history studies.

The IRIS study collected specimens and data from individuals managed according to KPNC clinical practice guidelines. Therefore, management and disease ascertainment were based on screening and management test results. For example, screen-positive individuals are more likely to have repeat sampling than screen-negative individuals and risk analyses require sophisticated statistical methodologies such as prevalence-incidence mixture models (18, 38–40). Furthermore, all assays and natural history studies must be conducted from one of the two clinical collection mediums, STM and SurePath, used at KPNC and collected in IRIS. Although clinical follow-up data is available for 83.7% of the IRIS cohort, 24.7% were no longer KPNC members, preventing the acquisition of long-term outcomes for those individuals.

Beyond cervical endpoints, the IRIS study also is well-placed to study markers for endometrial and gynecologic cancers because cells from the genital tract, including the fallopian tubes and the endometrium shed and can be detected in cervical samples (41, 42). Given the large population of individuals with repeat cervical samples, the IRIS study is amenable to evaluate biomarkers for gynecologic cancers such as somatic mutations and methylation for carcinogenesis and early detection studies. In summary, IRIS is a large, unique cohort that allows for natural history studies and head—to–head evaluation of established assays and new biomarkers for cervical cancer screening and management.

T. Raine-Bennett reports other support from NCI of NIH during the conduct of the study. No disclosures were reported by the other authors.

The content is solely the responsibility of the authors and does not necessarily represent the official views of the NIH.

J.C. Gage: Conceptualization, formal analysis, supervision, validation, investigation, methodology, writing–original draft, project administration, writing–review and editing. T. Raine-Bennett: Conceptualization, data curation, supervision, methodology, project administration, writing–review and editing. M. Schiffman: Conceptualization, formal analysis, investigation, visualization, methodology, writing–review and editing. M.A. Clarke: Formal analysis, validation, investigation, writing–review and editing. L.C. Cheung: Data curation, software, formal analysis, validation, writing–review and editing. N.E. Poitras: Conceptualization, data curation, software, formal analysis, supervision, validation, investigation, methodology, writing–review and editing. N.E. Varnado: Data curation, supervision, investigation, visualization, methodology, project administration. H.A. Katki: Conceptualization, data curation, formal analysis, validation, investigation, visualization, writing–review and editing. P.E. Castle: Conceptualization, data curation, software, investigation, visualization, writing–review and editing. B. Befano: Data curation, software, formal analysis, supervision, validation, investigation, visualization, methodology, writing–review and editing. M. Chandra: Data curation, software, validation, investigation. G. Rydzak: Data curation, software, formal analysis, validation, writing–review and editing. T. Lorey: Conceptualization, resources, supervision, validation, project administration, writing–review and editing. N. Wentzensen: Conceptualization, data curation, formal analysis, supervision, funding acquisition, validation, investigation, methodology, project administration, writing–review and editing.

The authors acknowledge the extensive efforts of laboratory study staff, including Stephanie Phelps, Oscar Lee, Kennedy Hardemion, Howard Moffett, and Chinar Sheth. This paper is dedicated to the memory of Barbara Fetterman, who coordinated the KP Regional Lab through 2017 and coordinated the KP-NCI collaborations. Her long-term determination, tenacity, and commitment were critical for success of these projects. Research reported in this publication was supported by the Intramural research award of the NIH. This publication was supported by grant number T32CA09168 from the NIH.

The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

1.
Bosch
FX
,
Lorincz
A
,
Munoz
N
,
Meijer
CJ
,
Shah
KV
. 
The causal relation between human papillomavirus and cervical cancer
.
J Clin Pathol
2002
;
55
:
244
65
.
2.
Schiffman
M
,
Castle
PE
,
Jeronimo
J
,
Rodriguez
AC
,
Wacholder
S
. 
Human papillomavirus and cervical cancer
.
Lancet
2007
;
370
:
890
907
.
3.
Bouvard
V
,
Baan
R
,
Straif
K
,
Grosse
Y
,
Secretan
B
,
El Ghissassi
F
, et al
A review of human carcinogens–Part B: biological agents
.
Lancet Oncol
2009
;
10
:
321
2
.
4.
Ronco
G
,
Dillner
J
,
Elfstrom
KM
,
Tunesi
S
,
Snijders
PJ
,
Arbyn
M
, et al
Efficacy of HPV-based screening for prevention of invasive cervical cancer: follow-up of four European randomised controlled trials
.
Lancet
2014
;
383
:
524
32
.
5.
Dillner
J
,
Rebolj
M
,
Birembaut
P
,
Petry
KU
,
Szarewski
A
,
Munk
C
, et al
Long term predictive values of cytology and human papillomavirus testing in cervical cancer screening: joint European cohort study
.
BMJ
2008
;
337
:
a1754
.
6.
Fontham
ETH
,
Wolf
AMD
,
Church
TR
,
Etzioni
R
,
Flowers
CR
,
Herzig
A
, et al
Cervical cancer screening for individuals at average risk: 2020 guideline update from the American Cancer Society
.
CA Cancer J Clin
2020
;
70
:
321
46
.
7.
Kyrgiou
M
,
Arbyn
M
,
Bergeron
C
,
Bosch
FX
,
Dillner
J
,
Jit
M
, et al
Cervical screening: ESGO-EFC position paper of the European Society of Gynaecologic Oncology (ESGO) and the European Federation of Colposcopy (EFC)
.
Br J Cancer
2020
;
123
:
510
7
.
8.
Perkins
RB
,
Guido
RS
,
Castle
PE
,
Chelmow
D
,
Einstein
MH
,
Garcia
F
, et al
2019 ASCCP risk-based management consensus guidelines for abnormal cervical cancer screening tests and cancer precursors
.
J Low Genit Tract Dis
2020
;
24
:
102
31
.
9.
Wheeler
CM
,
Hunt
WC
,
Joste
NE
,
Key
CR
,
Quint
WG
,
Castle
PE
. 
Human papillomavirus genotype distributions: implications for vaccination and cancer screening in the United States
.
J Natl Cancer Inst
2009
;
101
:
475
87
.
10.
Wentzensen
N
,
Fetterman
B
,
Castle
PE
,
Schiffman
M
,
Wood
SN
,
Stiemerling
E
, et al
p16/Ki-67 dual stain cytology for detection of cervical precancer in HPV-positive women
.
J Natl Cancer Inst
2015
;
107
:
djv257
.
11.
Clarke
MA
,
Cheung
LC
,
Castle
PE
,
Schiffman
M
,
Tokugawa
D
,
Poitras
N
, et al
Five-year risk of cervical precancer following p16/Ki-67 dual-stain triage of HPV-positive women
.
JAMA Oncol
2019
;
5
:
181
6
.
12.
Wentzensen
N
,
Lahrmann
B
,
Clarke
MA
,
Kinney
W
,
Tokugawa
D
,
Poitras
N
, et al
Accuracy and efficiency of deep-learning-based automation of dual stain cytology in cervical cancer screening
.
J Natl Cancer Inst
2021
;
113
:
72
9
.
13.
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
.
Clin Cancer Res
2018
;
24
:
2194
202
.
14.
Kelly
H
,
Benavente
Y
,
Pavon
MA
,
De Sanjose
S
,
Mayaud
P
,
Lorincz
AT
. 
Performance of DNA methylation assays for detection of high-grade cervical intraepithelial neoplasia (CIN2+): a systematic review and meta-analysis
.
Br J Cancer
2019
;
121
:
954
65
.
15.
Bowden
SJ
,
Kalliala
I
,
Veroniki
AA
,
Arbyn
M
,
Mitra
A
,
Lathouras
K
, et al
The use of human papillomavirus DNA methylation in cervical intraepithelial neoplasia: A systematic review and meta-analysis
.
EBioMedicine
2019
;
50
:
246
59
.
16.
Demarco
M
,
Egemen
D
,
Raine-Bennett
TR
,
Cheung
LC
,
Befano
B
,
Poitras
NE
, et al
A study of partial human papillomavirus genotyping in support of the 2019 ASCCP risk-based management consensus guidelines
.
J Low Genit Tract Dis
2020
;
24
:
144
7
.
17.
Egemen
D
,
Cheung
LC
,
Chen
X
,
Demarco
M
,
Perkins
RB
,
Kinney
W
, et al
Risk estimates supporting the 2019 ASCCP risk-based management consensus guidelines
.
J Low Genit Tract Dis
2020
;
24
:
132
43
.
18.
Cheung
LC
,
Egemen
D
,
Chen
X
,
Katki
HA
,
Demarco
M
,
Wiser
AL
, et al
2019 
ASCCP risk-based management consensus guidelines: methods for risk estimation, recommended management, and validation
.
J Low Genit Tract Dis
2020
;
24
:
90
101
.
19.
Arbyn
M
,
Andersson
K
,
Bergeron
C
,
Bogers
JP
,
von Knebel-Doebertitz
M
,
Dillner
J
. 
Cervical cytology biobanks as a resource for molecular epidemiology
.
Methods Mol Biol
2011
;
675
:
279
98
.
20.
Perskvist
N
,
Norlin
L
,
Dillner
J
. 
The process of moving from a regionally based cervical cytology biobank to a national infrastructure
.
Biopreserv Biobank
2015
;
13
:
94
7
.
21.
LaMere
BJ
,
Kornegay
J
,
Fetterman
B
,
Sadorra
M
,
Shieh
J
,
Castle
PE
. 
Human papillomavirus genotyping after denaturation of specimens for Hybrid Capture 2 testing: feasibility study for the HPV persistence and progression cohort
.
J Virol Methods
2007
;
146
:
80
5
.
22.
Raine-Bennett
T
,
Gage
JC
,
Poitras
N
,
Chandra
M
,
Varnado
N
,
Befano
B
, et al
Development of a large biorepository of cervical specimens for the improving risk informed HPV screening study (IRIS)
.
J Clin Virol
2021
;
145
:
105014
.
23.
Wentzensen
N
,
Wacholder
S
. 
From differences in means between cases and controls to risk stratification: a business plan for biomarker development
.
Cancer Discov
2013
;
3
:
148
57
.
24.
Krieger
N
. 
Overcoming the absence of socioeconomic data in medical records: validation and application of a census-based methodology
.
Am J Public Health
1992
;
82
:
703
10
.
25.
Katki
HA
,
Kinney
WK
,
Fetterman
B
,
Lorey
T
,
Poitras
NE
,
Cheung
L
, et al
Cervical cancer risk for women undergoing concurrent testing for human papillomavirus and cervical cytology: a population-based study in routine clinical practice
.
Lancet Oncol
2011
;
12
:
663
72
.
26.
Landy
R
,
Schiffman
M
,
Sasieni
PD
,
Cheung
LC
,
Katki
HA
,
Rydzak
G
, et al
Absolute risks of cervical precancer among women who fulfill exiting guidelines based on HPV and cytology cotesting
.
Int J Cancer
2020
;
146
:
617
26
.
27.
Castle
PE
,
Kinney
WK
,
Xue
X
,
Cheung
LC
,
Gage
JC
,
Poitras
NE
, et al
Role of screening history in clinical meaning and optimal management of positive cervical screening results
.
J Natl Cancer Inst
2019
;
111
:
820
7
.
28.
Castle
PE
,
Kinney
WK
,
Xue
X
,
Cheung
LC
,
Gage
JC
,
Zhao
FH
, et al
Effect of several negative rounds of human papillomavirus and cytology co-testing on safety against cervical cancer: an observational cohort study
.
Ann Intern Med
2018
;
168
:
20
9
.
29.
Schiffman
M
,
Kinney
WK
,
Cheung
LC
,
Gage
JC
,
Fetterman
B
,
Poitras
NE
, et al
Relative performance of HPV and cytology components of cotesting in cervical screening
.
J Natl Cancer Inst
2018
;
110
:
501
8
.
30.
Rendle
KA
,
Schiffman
M
,
Cheung
LC
,
Kinney
WK
,
Fetterman
B
,
Poitras
NE
, et al
Adherence patterns to extended cervical screening intervals in women undergoing human papillomavirus (HPV) and cytology cotesting
.
Prev Med
2018
;
109
:
44
50
.
31.
Clarke
MA
,
Fetterman
B
,
Cheung
LC
,
Wentzensen
N
,
Gage
JC
,
Katki
HA
, et al
Epidemiologic evidence that excess body weight increases risk of cervical cancer by decreased detection of precancer
.
J Clin Oncol
2018
;
36
:
1184
91
.
32.
Arbyn
M
,
Rezhake
R
,
Yuill
S
,
Canfell
K
. 
Triage of HPV-positive women in Norway using cytology, HPV16/18 genotyping and HPV persistence
.
Br J Cancer
2020
;
122
:
1577
9
.
33.
Pinheiro
M
,
Gage
JC
,
Clifford
GM
,
Demarco
M
,
Cheung
LC
,
Chen
Z
, et al
Association of HPV35 with cervical carcinogenesis among women of African ancestry: Evidence of viral-host interaction with implications for disease intervention
.
Int J Cancer
2020
;
147
:
2677
86
.
34.
Zhu
B
,
Xiao
Y
,
Yeager
M
,
Clifford
G
,
Wentzensen
N
,
Cullen
M
, et al
Mutations in the HPV16 genome induced by APOBEC3 are associated with viral clearance
.
Nat Commun
2020
;
11
:
886
.
35.
Mirabello
L
,
Yeager
M
,
Yu
K
,
Clifford
GM
,
Xiao
Y
,
Zhu
B
, et al
HPV16 E7 genetic conservation is critical to carcinogenesis
.
Cell
2017
;
170
:
1164
74
.
36.
Mirabello
L
,
Frimer
M
,
Harari
A
,
McAndrew
T
,
Smith
B
,
Chen
Z
, et al
HPV16 methyl-haplotypes determined by a novel next-generation sequencing method are associated with cervical precancer
.
Int J Cancer
2015
;
136
:
E146
53
.
37.
Usyk
M
,
Zolnik
CP
,
Castle
PE
,
Porras
C
,
Herrero
R
,
Gradissimo
A
, et al
Cervicovaginal microbiome and natural history of HPV in a longitudinal study
.
PLoS Pathog
2020
;
16
:
e1008376
.
38.
Landy
R
,
Cheung
LC
,
Schiffman
M
,
Gage
JC
,
Hyun
N
,
Wentzensen
N
, et al
Challenges in risk estimation using routinely collected clinical data: The example of estimating cervical cancer risks from electronic health-records
.
Prev Med
2018
;
111
:
429
35
.
39.
Cheung
LC
,
Pan
Q
,
Hyun
N
,
Schiffman
M
,
Fetterman
B
,
Castle
PE
, et al
Mixture models for undiagnosed prevalent disease and interval-censored incident disease: applications to a cohort assembled from electronic health records
.
Stat Med
2017
;
36
:
3583
95
.
40.
Hyun
N
,
Cheung
LC
,
Pan
Q
,
Schiffman
M
,
Katki
HA
. 
Flexible risk prediction models for left or interval-censored data from electronic health records
.
Ann Appl Stat
2017
;
11
:
1063
84
.
41.
Kinde
I
,
Bettegowda
C
,
Wang
Y
,
Wu
J
,
Agrawal
N
,
Shih Ie
M
, et al
Evaluation of DNA from the Papanicolaou test to detect ovarian and endometrial cancers
.
Sci Transl Med
2013
;
5
:
167ra4
.
42.
Wang
Y
,
Li
L
,
Douville
C
,
Cohen
JD
,
Yen
TT
,
Kinde
I
, et al
Evaluation of liquid from the Papanicolaou test and other liquid biopsies for the detection of endometrial and ovarian cancers
.
Sci Transl Med
2018
;
10
:
eaap8793
.
This work is licensed under a Creative Commons Attribution-NonCommercial-NoDerivs International 4.0 License.