Methylation of High-Risk Human Papillomavirus Genomes Are Associated with Cervical Precancer in HIV-Positive Women

Human papillomavirus (HPV) infection is a common sexually transmitted infection (STI) and most sexually active women will acquire at least one HPV infection at some point during their lives (1). Persistent infection by a high-risk HPV (HR-HPV) type is the cause of nearly all cervical cancers (2). However, only a small number of HR-HPV infections result in cervical precancerous lesions, of which a third go on to become cancers (1, 3). Worldwide, cervical cancer is the fourth most common cancer and cause of cancer-related deaths in women (4), and there is especially high incidence of cervical precancer and cancer among women with HIV (5).

Cervical cancer screening has largely depended on cytology and more recently HPV DNA testing (6). However, their combined (cotest) positive predictive value (PPV) remains moderate (7). The majority of women who are cotest positive do not have clinically relevant cervical disease upon colposcopy and biopsy. Therefore, biomarkers capable of improving cervical cancer screening specificity and PPV are greatly needed, especially as several generations of women are too old to receive HPV vaccination, both in high-resource and limited-resource settings.

Epigenetic changes of the HPV viral genome have been identified as a promising biomarker for cervical precancer and cancer (8). Viral genome methylation can influence protein expression, replication, and susceptibility to immune surveillance, cell proliferation, and risk of virus-related disease. DNA methylation occurs preferentially in regions where a cytosine occurs next to a guanine nucleotide linked to one another by a phosphate bond (termed a CpG site). As we recently reviewed (8), there are no classical CpG islands within the HPV genome, but there are regions of high CpG density, and conservation of CpG sites across HPV types suggest these sites are biologically important (8–11). A correlation has been made between high levels of CpG methylation within the HPV genome and cervical precancer development (7, 12), especially methylation levels in the HPV capsid L1 and L2 genes (13).

Despite the high rates of HPV and cervical disease in HIV-positive women, recent prospective studies have shown that HPV cotesting in HIV-positive women has high sensitivity and negative predictive values (NPV), as in HIV-negative women; strongly discriminating between those at high and low-risk of current as well as subsequent cervical precancer (14). Therefore, recent cervical cancer screening guidelines in the United States recommend the use of HR-HPV cotesting in women infected with HIV (15); albeit, as in the general population, the specificity and PPV of cotesting in HIV-positive women are moderate. Given the much higher rate of positive cotests in HIV-infected women than in the general population, the need for biomarkers capable of improving the accuracy of cervical cancer screening in HIV-positive women is especially great.

However, there are little data regarding HPV CpG methylation and precancer in HIV-positive women. This issue is compounded by the diversity of HR-HPV types detected in HIV-positive precancer and cancer cases (16), as it requires methylation assays for a wide range of HR-HPV types to appropriately characterize the relation of HPV CpG methylation with precancer. Therefore, we developed high-throughput next-generation sequencing (NGS) HPV CpG methylation assays able to assess HR-HPV types 16/18/31/33/35/39/45/51/52/56/58/59. The current investigation was designed as a proof-of-concept study of these new assays, and to determine whether there are differences in HPV CpG methylation levels between HIV-positive women with and without cervical precancer. We also evaluated whether HPV CpG methylation could be used in precancer cases with multiple concurrent HR-HPV infections to identify the type causally related to the lesion. That is, we hypothesized that one HR-HPV type would have markedly higher methylation than other types detected, consistent with there being a single HR-HPV type causally related to each lesion.

The Women's Interagency HIV Study (WIHS) is a large ongoing cohort of HIV-positive (n = 2,793) and HIV-negative (n = 975) women enrolled from six clinical sites (Brooklyn, NY; Bronx, NY; Washington, DC; Chicago, IL; San Francisco, CA; Los Angeles, CA), as described previously (14, 17). These women are followed every 6 months with questionnaires, general physical and pelvic examinations, including Papanicolaou (Pap) smear and cervicovaginal lavage (CVL) for HPV testing. Written informed consent was obtained from all participants, and the WIHS study was approved by each local institutional review board.

The current case–control study was performed only in HIV-infected women. As CIN3 is the established precancerous lesion, cases were limited to CIN3 (N = 72, CIN3) selected from women with available HPV DNA results. Controls had no CIN2 or greater at any time in the WIHS and were initially matched 1:1 to cases based on age, CD4 count, and HPV type. The HPV type present was defined as the most oncogenic type present in each precancer case defined hierarchically (HPV16 > HPV18 > other HPV types) and selected so that cases would include N = 20 positive for HPV16, and N = 25 other alpha-9 HPV species (i.e., HPV31/33/35/52/58); a total of N = 25 alpha-7 species (i.e., HPV18/39/45/59); and as many alpha-5/alpha-6 (i.e., HPV51 and 56) as available. Those precancers and controls with the fewest HPV types present were prioritized for inclusion. Overall, 1 fewer HPV16 (N = 19), and 1 greater alpha-7 (N = 26) than planned were available for the current pilot study. Too few alpha-5/alpha-6 were available to analyze individually. Furthermore, approximately 32% (23/72) of cases and 13% (9/72) controls had multiple HR-HPV infections (i.e., ≥2 HR-HPV types), and all data were used to maximize the efficiency of this proof-of-concept investigation (see Statistical analysis).

The HPV DNA detection methods have been described in detail elsewhere (16, 18, 19). Briefly, HPV DNA was detected using a well-established MY09/11 primer system PCR assay and probed for more than 40 individual HPV types using filters hybridized with type-specific biotinylated oligonucleotides. Consistent with recommendations from the International Association for Research on Cancer (IARC), HPV types 16/18/31/33/35/39/45/51/52/56/58/59 were considered oncogenic (20, 21).

Isolated cervicovaginal DNA was treated with bisulfite using the Qiagen Epitect Bisulfite Kit (Qiagen) as recommended by the manufacturer. Negative and positive controls, DNase-free water and HPV-positive DNA, respectively, were processed as described for the clinical samples. SiHa (ATCC HTB-35) and HeLa (ATCC CCL-2) cells were obtained from ATCC in January 2017, shipped frozen to Albert Einstein College of Medicine (Bronx, NY), and maintained for 20 passages in DMEM supplemented with 10% FBS. The cells tested negative for Mycoplasma in March 2017 using the e-Myco Mycoplasma PCR Detection Kit 2.0 (iNtRON Biotechnology) prior to DNA extraction using QIAamp DNA Blood Mini kit (Qiagen). The cell lines were authenticated in May 2017 by mapping HPV integration sites described previously (22, 23) using a deep-sequencing NGS methodology. Following bisulfite treatment and subsequent PCR, unmethylated cytosines (C's) were converted to uridine (U's) and then thymidine (T's), whereas methylated C's remain unmodified (8).

Up to three regions within the viral genome (E2E4/L1/L2) of HR-HPV types 16/18/31/33/35/39/45/51/52/56/58/59) were selected for quantitative site-specific CpG methylation using NGS assays. Methylation-specific primers (Integrated DNA Technologies) for bisulfite-converted DNA were designed using MethPrimer (24; http://www.urogene.org/methprimer/index1.html). Primers were designed to amplify fragments comprising at least 3 CpG sites per selected region. Each forward primer contained a unique 12 base pair (bp) Golay barcode attached to the 5′ end of the forward primer sequence (ref. 7; see Bioinformatics). Reverse primers were also labeled with a 12 bp barcode to distinguish each sample in downstream bioinformatics analyses while increasing the number of multiplexing combinations (7).

Barcoded PCR was performed using the Qiagen PyroMark PCR Kit, allowing enhanced specificity of primer annealing to increase amplification efficiency in bisulfite-converted DNA sequences. Each PCR reaction contained 1 μL of bisulfite converted DNA, 2.5 mmol/L MgCl₂, 1× PyroMark Master Mix, 5 μmol/L of forward and reverse primers in a total volume of 25 μL. The conditions and annealing temperatures for each assay were as follows: initial denaturation step at 95°C for 15 minutes followed by 50 cycles: 30 sec at 94°C; 30 sec at the optimized primer-specific annealing temperature (56°C–58°C); 30 sec at 72°C, and a final extension for 10 minutes at 72°C.

PCR products were evaluated by 3% gel electrophoresis stained with ethidium bromide to confirm the expected amplicon sizes (range = 148–292 bp). Samples with and without visible bands in the gel were pooled together as part of each assay using 2 μL or 20 μL, respectively, (i.e., independent pools per each PCR assay) to adjust for disproportionate numbers of reads from samples with a visible DNA band. PCR amplicons from the pooled samples were purified by gel excision using Qiagen QIAquick Gel Extraction Kit according to manufacturer's guidelines, and then used to create libraries for NGS.

Libraries were constructed using KAPA LTP Library Preparation Kit for Illumina platforms (Kapa Biosystems) according to the manufacturer's instructions. Libraries were sequenced using 300 bp paired-end reads on an Illumina MiSeq at the Albert Einstein College of Medicine Epigenomics Core Facility (Bronx, NY). Overall, percent methylation of 117 CpG sites was evaluated for the selected HR-HPV types, including those reported previously (13, 25, 26).

Illumina sequencing data files (13.4 M reads) were first filtered for low-quality reads using Prinseq v0.20.4 (27) for a minimum average Phred score of 25 (11.4 M reads; ref. 28). After preprocessing of the Illumina sequencing files, an in-house pipeline was used to extract CpG methylation percentage information by demultiplexing the reads using the novobarcode computer package (Novocraft Technologies Sdn Bhd), based on the 12 bp Golay forward and reverse barcodes uniquely assigned to each sample as described above (7.1 M reads). Bismark v0.15.3 and Bowtie v2.2.3 (29, 30) with samtools v1.2 were used for paired-end mapping to HR-HPV reference genomes bisulfite converted in silico to determine the methylation state of CpG and non-CpG cytosines. On average, the number of reads per fragment was 17,713 (range 1,409–120,461) with the exception of HPV58, with average read counts less than 5,000. Custom R v3.2.1 (31) scripts were incorporated in the pipeline to process the Bismark output by: (i) producing a methylation pattern for every sequenced molecule, indicating whether cytosines were methylated (+), unmethylated (−), missing (o) or in disagreement for paired-end reads (x); (ii) generating counts of each unique pattern per sample; and (iii) calculating site-specific CpG methylation percentages for each assayed cytosine per sample by comparing the ratio of methylated C's to the total number of methylated and unmethylated (C+T) at each CpG site screened. The ratio of “C/(C + T)” indicates the proportion of methylated cytosines at each CpG site for the assayed sample.

Reproducibility of intraplate duplicate PCR-positive samples for each assay was evaluated on the basis of Pearson and Spearman correlation coefficients and all were R² > 0.7. Individual site-specific CpG methylation levels were analyzed by HR-HPV type using median and range within each PCR fragment assayed using the Mann–Whitney U test. Standard logistic regression (rather than conditional logistic regression) was performed to determine case–control differences in percent methylation distributions for every CpG site by defining a tertile cut-off (T3 vs. T1 and T2, i.e., high methylation vs. low methylation) based on data in controls. This approach made it possible to incorporate data for multiple different HPV types from cases and controls with ≥2 HPV types, and thereby maximize statistical power in this proof-of-concept pilot study, which would not have been possible with conditional logistic regression models. Note, however, that for HPV type–specific analyses each case and control only contributed a single observation (i.e., for that one HPV type). Thus, statistical methods to address multiple observations per subject were not applicable.

Conversely, for analyses involving multiple homologous CpG sites, we used logistic regression models that incorporated a generalized estimating equation (GEE), as women could contribute data to more than one HPV type even if from the same phylogenetic group, resulting in multiple observations, and all models were adjusted for the matching variables, age and CD4 count. To conduct these analyses, tertiles were separately determined for each HPV type and CpG site using the methylation level of controls to define cut-off values. The data for all HPV types and CpG sites within a given grouping of homologous sites were then combined and the GEE models were fitted to estimate the ORs, 95% confidence intervals, and P values. Thus, the primary a priori hypothesis formally tested using multivariate models in this study was that phylogenetically related HPV types with homologous CpG sites would, as a group, show greater methylation at selected CpG sites among cervical precancer cases than controls; an approach which maximized statistical power. SAS 9.3 (SAS Institute Inc) was used to conduct all analyses. All statistical tests were two-sided, with P < 0.05 used to define statistical significance.

We hypothesized that in precancer cases, one HR-HPV type would have markedly higher methylation than other types detected, consistent with there being a single HR-HPV type causally related to each lesion, whereas among controls (in the absence of disease) we predicted there would be no standout HR-HPV type with high methylation levels relative to the other HPV types detected. Therefore, we assessed the heterogeneity in methylation levels between HR-HPV types within individuals and then summary differences between cases and controls. Briefly, for each case and control with multiple HR-HPV types, a Z-score was calculated as “Z = (m - μ_control)/σ_control” for each type, where μ and σ were respectively the mean and SD of the methylation values. The Z-scores measured how high the methylation level was relative to the mean level across all control samples with that HPV type for a given CpG site/HR-HPV type. A heatmap of these values among cases and controls with multiple HR-HPV types was created using a custom python script incorporating matplotlib and seaborn software packages (32, 33), and then used to visualize the heterogeneity of methylation levels. We hypothesized that cases but not controls would have one HR-HPV type with markedly higher methylation than the other types detected. Therefore, we used the Kruskal–Wallis test to compare the average differences in Z-scores among cases relative to that in controls. The expectation was that cases would have significantly greater differences than controls, as cases would more often have a single HR-HPV type with much higher methylation levels than all other HR-HPV detected.

Table 1 shows selected characteristics of cases and matched controls at the time of case diagnosis. These subjects were matched by age (median = 38 and 37 years, for cases and controls, respectively), race/ethnicity (e.g., 61% and 63% Black), current CD4 count (median = 237 and 278), and HR-HPV type (e.g., HPV16 was detected in 26% of cases and 26% of controls). Among the other variables assessed, the only differences that reached statistical significance were parity (median = 3 vs. 2 live births).

To evaluate the association of HPV CpG methylation with CIN3, site-specific CpG methylation levels were determined at a total of 117 CpG sites across each of 12 HR-HPV types (Fig. 1; Supplementary Table S1). For HPV16, cases were significantly more likely to have higher median CpG methylation levels at each CpG site within the L1 fragment (all sites P ≤ 0.01). The ORs comparing the highest to the lowest two tertiles (T3 vs. T1 and T2) for the four L1 CpG sites ranged from 4.69 (95% CI: 1.20–18.4; P = 0.03) to 8.12 (95% CI: 1.87–35.2; P = 0.005). In contrast, the L2 CpG methylation levels for HPV16 were not significantly associated with precancer (Supplementary Table S1).

Other alpha-9 HR-HPV types (HPV31/35/52/58) also showed significant case–control differences at specific CpG sites (Supplementary Table S1). For example, there were five HPV31 L2 CpG sites with median methylation levels positively associated with case status (P ≤ 0.01), and based on tertiles, the strongest association was OR = 28.0 (95% CI: 2.82–278; P = 0.004; see Supplementary Table S1). HPV35 L2 had five CpG sites with significant differences in median methylation levels between cases and controls, and by tertiles the strongest association was OR = 22.0, 95% CI: 1.54–314 (P = 0.02). Similarly to our observations in HPV16 L1 region, HPV35 L2/L1 region showed significant associations in median levels between cases and controls (Supplementary Table S1). HPV58 had several CpG sites in the E2/E4 region (rather than L1/L2) that showed significant positive associations of methylation levels with case status by the Mann–Whitney U test (Supplementary Table S1).

Similarly, for alpha-7 HPV types (HPV18/39/45) there were significant case–control differences at specific CpG sites. HPV18 had increased methylation at two L1 CpG sites (e.g., HPV18-L1-7062: P = 0.02; OR = 7.87, 95% CI: 1.10–56.1, P = 0.04), and HPV39 methylation was positively associated with case status at a CpG site located at the transition between the L1 and L2 genes (CpG#: HPV39-L2L1-5731). For HPV45, there were an insufficient number of cases and controls to adequately analyze, although the results were in the expected direction (Supplementary Table S1).

We scanned CpG locations within L1 and L2 of the alpha-9 species group, and in the alpha-7 species group, to identify homologous sites that could be combined for analysis. Figure 2A shows the fragment alignments for all 12 HR-HPV types surveyed and identifies multiple homologous CpG sites present in some or all of the genomes (vertical dark bars). We posited that phylogenetically related HR-HPV types should have similar case–control associations at these homologous CpG sites. Tertile methylation levels (T3 vs. T1 and T2) were used to categorize these data, rather than absolute values (which may vary between HR-HPV types and CpG sites), as discussed in the Statistical analysis. The alpha-9 species analyses showed significant results for the combined data from HPV16L1, HPV31L2, HPV35L2L1, and similar results were observed for other HPV types (see Fig. 2B). For example, significant case–control differences were observed across the homologous alpha-9 CpG sites HPV16-L1-5608/HPV31-L1-5521/HPV35-L2L1-5570 (OR = 7.28; 95% CI: 2.75–19.3; P = 0.0001), and alpha-7 CpG sites HPV18-L1-7062/HPV45-L1-7066 (OR = 6.94; 95% CI: 1.23–39.3; P = 0.03).

We addressed the hypothesis that, for cases but not controls, one HR-HPV type would have markedly higher methylation than the other types detected, consistent with there being a single HR-HPV type causally related to each lesion. First, a Z-score was determined for each CpG site in each HR-HPV in all women with multiple HR-HPV types. A heatmap of Z-score levels (Fig. 3) was then generated, and confirmed that most precancer cases had a single HR-HPV with much higher methylation levels than all other HR-HPV detected, whereas this was not true for controls. To formally quantify this case–control difference, for each individual woman, we calculated the absolute (subtracted) difference between the highest Z-score and the Z-score for each of the other HR-HPVs (using the mean of these subtracted differences to represent that woman). To summarize these results, we then compared the average of the individual differences for cases and controls using the Kruskal–Wallis test. As predicted, the average difference among cases was significantly greater than for controls (see Fig. 3 insert; P = 3.041e⁻⁵).

This study found significant differences in HR-HPV CpG methylation levels between HIV-positive women with and without cervical precancer. Prior studies in HIV-negative populations also found that HR-HPV methylation levels are strongly associated with detection of cervical precancer and cancer, and might be a promising biomarker with potential to improve the specificity and PPV of primary HR-HPV testing or of cotesting (HR-HPV testing in combination with Pap tests; refs. 7, 8, 12, 25, 34–40). Improving the accuracy of cervical cancer screening is important as the majority of HIV-positive and HIV-negative women who are HR-HPV or cotest positive do not have clinically relevant cervical disease upon colposcopy and biopsy. Because a diversity of HR-HPV causes cervical precancer in HIV-positive women, we developed NGS methylation assays for 12 different HR-HPVs. The development and use of these NGS assays contributes to the clinical relevance of this proof-of-concept study, as these high-throughput assays allow specimens from a large number of women to be concurrently assessed, which will be important for use in real-world cervical cancer screening programs (7). In addition, the current data provide interesting preliminary evidence that HR-HPV CpG methylation levels might be useful in determining the causal HPV-type in precancer and cancer cases who test positive for multiple HR-HPV types. From a biological perspective, furthermore, the study also provided evidence that homologous CpG sites across phylogenetically related HR-HPV types share similar associations with cervical disease, and probably epigenetic mechanisms of carcinogenesis.

More specifically, we observed several strong cervical precancer case–control associations in HIV-positive women for alpha-9 HR-HPV types, especially HPV16, 31, and 35, as well as moderate associations for HPV58. The data were also suggestive of case–control differences related to methylation of HPV18 and other alpha-7 HR-HPV types, HPV39 and 45. We then combined tertile data across homologous L1/L2 CpG sites of related HR-HPVs, and found case–control differences across homologous alpha-9 CpG sites, and alpha-7 CpG sites, indicating a similar role in cervical disease for related HR-HPV, though the mechanisms underlying these relationships remain an important area of investigation.

Multiple concurrent HR-HPV infections are common in HIV-positive women, but have not been extensively studied in the context of type-specific CpG methylation levels. Because several sources of evidence have shown that each CIN3 and cancer lesions are caused by a single HR-HPV type (25, 40), we used a novel approach to determine whether HPV CpG methylation patterns could be used to identify the HR-HPV causally associated with CIN3 in HIV-positive cases with multiple concurrent HPV. First we evaluated each HR-HPV genome for the CpG methylation site most strongly associated with precancer based on Z-scores. Then we tested and confirmed the hypothesis that one HR-HPV type would have markedly higher methylation than all other types detected, consistent with there being a single HR-HPV type causally related to each lesion. While this was not true in every precancer case with multiple HR-HPV types, it was observed in most. In fact, for a few cases it was not possible to ascertain the causal HR-HPV type using the Z-score even though it was more accurate than comparing the heatmap color coding. This could possibly be due to the development of a recent cervical precancerous lesion that might still accumulate increasing viral methylation over time (i.e., for cases where the Z-score was not elevated above control levels). These findings now warrant validation in studies that involve HPV PCR testing of microdissected histologic specimens to confirm the “causal” HPV type inferred from HPV CpG methylation levels of exfoliated cervical cells. This was not possible in the current study as we did not have histologic specimens from WIHS women. If confirmed, though, the findings may be particularly relevant for use in limited resource settings, where inexpensive high-throughput molecular HPV assays are increasingly being used, as recommended by the World Health Organization (41). Persistence of the same HR-HPV following precancer treatment is an important risk factor for disease recurrence (42–44).

Although this study had several strengths, including a well-characterized cohort of HIV-positive women with detailed cervical HPV, cytology, and histologic data, there were also limitations that warrant discussion. In particular, as a proof-of-concept study, sample size and power were inadequate to conduct subgroup analyses, some effect estimates were based on few cases, and we did not correct for multiple comparisons. The main a priori hypothesis formally tested using multivariate models was that phylogenetically related HPV types with homologous CpG sites would, as a group, show greater methylation at selected CpG sites among cervical precancer cases than controls. By prioritizing selection of cases and controls with few HPV types present for inclusion in this pilot study, our design under-represented the true proportion with multiple HPV types in the WIHS cohort, which will need to be addressed in future larger and more comprehensive studies. On the other hand, a number of a priori predicted associations were highly significant, and most additional results were in keeping with our initial hypotheses. The cross-sectional design was another major limitation and, therefore, the case–control associations we observed may represent reverse causality. That said, even if the HPV methylation–precancer relationships we characterized are not causal associations, our findings and assays could still be clinically important as biomarkers for secondary cancer prevention in cervical cancer screening (15, 16), whereas causality will need to be studied in appropriate prospective studies.

Overall, this study presents novel initial data showing that cervical precancer in HIV-positive women is associated with high levels of methylation at CpG sites in the L1/L2 region. In combination with prior data in the general population, these findings raise the possibility that with continued improvements the high-throughput NGS methylation assays we utilized could help improve the accuracy of HR-HPV testing and cotesting in cervical cancer screening. These assays might even be useful in determining the HPV-type causally associated with precancer in those cases with multiple HR-HPV. However, future laboratory studies are warranted to understand the mechanisms underlying the relation of L1/L2 methylation and cervical disease development. Larger, more comprehensive prospective molecular epidemiologic studies will be needed to understand temporality between methylation and cervical disease, to determine whether early risk stratification is possible, and understand the possible role of L1/L2 methylation in the natural history of HR-HPV and cervical disease progression.

H.D. Strickler reports receiving other commercial research support from Roche, BD Diagnositics, and Arbor Vita. No potential conflicts of interest were disclosed by the other authors.

The contents of this publication are solely the responsibility of the authors and do not represent the official views of the NIH.

Conception and design: A. Gradissimo, H.D. Strickler, R.D. Burk

Development of methodology: A. Gradissimo, X. Xie, R.D. Burk

Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): J. Palefsky, I.-E. Eltoum, M.A. Fischl, K. Anastos, H. Minkoff, L.C. Flowers, N.A. Hessol, H.D. Strickler, R.D. Burk

Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): A. Gradissimo, J. Lam, J.D. Attonito, X. Xie, S. Shrestha, H.D. Strickler, R.D. Burk

Writing, review, and/or revision of the manuscript: A. Gradissimo, J.D. Attonito, J. Palefsky, L.S. Massad, X. Xie, I.-E. Eltoum, L. Rahangdale, M.A. Fischl, K. Anastos, H. Minkoff, G. D'Souza, L.C. Flowers, C. Colie, S. Shrestha, N.A. Hessol, H.D. Strickler, R.D. Burk

Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): A. Gradissimo, X. Xie, H. Minkoff, R.D. Burk

Study supervision: M.A. Fischl, H.D. Strickler, R.D. Burk

Other (interpretation of the data): M.A. Fischl

Other (PI of project that collected the samples and data): H.D. Strickler

The authors would like to acknowledge the contribution of the Epigenomics Core Facility at the Albert Einstein College of Medicine for providing assistance with the Illumina sequencing and helpful discussions on this work. HPV testing and the current analyses were supported by the NCI (grant numbers R01-CA-085178 and R01-CA-174634, to H.D. Strickler; P30-CA-013330). Data in this article were collected by the Women's Interagency HIV Study (WIHS). WIHS (principal investigators): UAB-MS WIHS (to M. Saag, M.C. Kempf, and D. Konkle-Parker), U01-AI-103401; Atlanta WIHS (to I. Ofotokun and G. Wingood), U01-AI-103408; Bronx WIHS (to K. Anastos), U01-AI-035004; Brooklyn WIHS (to H. Minkoff and D. Gustafson), U01-AI-031834; Chicago WIHS (to M. Cohen and A. French), U01-AI-034993; Metropolitan Washington WIHS (to S. Kassaye), U01-AI-034994; Miami WIHS (to M.A. Fischl and L. Metsch), U01-AI-103397; UNC WIHS (to A. Adimora), U01-AI-103390; Connie Wofsy Women's HIV Study, Northern California (to R. Greenblatt, B. Aouizerat, and P. Tien), U01-AI-034989; WIHS Data Management and Analysis Center (to S. Gange and E. Golub), U01-AI-042590; Southern California WIHS (to J. Milam), U01-HD-032632 (WIHS I – WIHS IV). The WIHS is funded primarily by the National Institute of Allergy and Infectious Diseases (NIAID; grant number P30-AI-050410, to R. Swanstrom), with additional cofunding from the Eunice Kennedy Shriver National Institute of Child Health and Human Development (NICHD), the NCI, the National Institute on Drug Abuse (NIDA), and the National Institute on Mental Health (NIMH). Targeted supplemental funding for specific projects is also provided by the National Institute of Dental and Craniofacial Research (NIDCR), the National Institute on Alcohol Abuse and Alcoholism (NIAAA), the National Institute on Deafness and other Communication Disorders (NIDCD), and the NIH Office of Research on Women's Health. WIHS data collection is also supported by UL1-TR000004 (UCSF CTSA) and UL1-TR000454 (Atlanta CTSA).

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.