Background: We used self-administered vaginal swabs to assess the incidence and clearance of carcinogenic human papillomavirus (HPV) infections in rural Rakai, Uganda.

Methods: Women provided self-administered vaginal swab at annual home-based visits. Type-specific carcinogenic HPV incidence and clearance and risk factors were assessed.

Results: Carcinogenic HPV incidence was 17.3 per 100 person-years among HIV-positive women compared with 7.0 per 100 person-years among HIV-negative women (P < 0.001). HPV-51 had the highest incidence followed by HPV-16 (1.8 per 100 and 1.5 per 100 person-years, respectively). In multivariate model, HIV-positive women were twice as likely to have incident infection compared with HIV-negative women. Younger women were at higher risk for incident infection, as were women with higher lifetime and recent sexual partners, and high perception of AIDS. Married women were less likely to have incident infection. Approximately half of all carcinogenic HPV infections cleared over the study follow-up of 3 years. HPV-31, HPV-35, and HPV-16 had the lowest clearance (16.7%, 27.9%, and 38.3%, respectively). In multivariate model, HIV-positive, women over 30 years with higher HPV viral, burden and more lifetime sex partners were less likely to clear infections.

Conclusions: Self-collected vaginal swabs provide accurate HPV exposure assessment for studying HPV exposure and epidemiology and can be an important tool for research in populations unwilling to undergo pelvic exam. (Cancer Epidemiol Biomarkers Prev 2008;17(6):1300–7)

Human papillomavirus (HPV) is the etiologic agent for cervical cancer and the most common viral sexually transmitted infection. HPV infections represent a major public health burden worldwide as a cause of cervical cancer in women, genital warts in men and women, other anogenital cancers, and some nasopharyngeal cancers. Acquisition of HPV is especially common among young sexually active adolescents with 3-year cumulative incidence estimated at more than 40% (1-3). Although HPV point prevalence is very high, most clear spontaneously, especially among young women (4, 5); however, persistent infection with high-risk HPV (HR-HPV) types is considered the crucial step in the development of cervical cancer (6).

As a consequence of cervical cytology screening, cervical cancer incidence and mortality has been drastically reduced in developed counties, whereas such a decline has not occurred in less developed countries where 80% of cervical cancers occur. Beyond lack of, or inadequate coverage of screening programs (7, 8) in developing countries, women rarely present for a routine pelvic exam if they are asymptomatic for disease. Hence, to increase screening coverage in settings with low compliance for pelvic exam–based screening, use of self-collected vaginal swabs has been widely advocated and evaluated in the literature (9-23), showing good agreement between self-collected vaginal swabs and clinician-collected cervical swabs for HPV DNA detection. This allows self-collection of vaginal samples to be used as a screening tool for HPV. Furthermore, the possibility of primary prevention of HPV infection and thus cervical cancer via HPV vaccination has been realized. Information on the epidemiology of HPV in understudied regions is fundamental to planning and evaluating future vaccine implementation programs. Finally, the prevention of HPV and cervical cancer is of particular urgency in Africa where there are an estimated 17 million women living with HIV (24).

We previously showed comparable HPV DNA detection from self-collected vaginal and physician-collected cervical swabs in this population (25). To investigate the utility of self-collected vaginal swabs beyond screening, in this analysis, our objective was to investigate the epidemiology (incidence and clearance) of HPV using self-collected vaginal swabs from a general rural African population. If our findings, based on self-collected vaginal samples, are similar to findings from other studies of the natural history of HPV, which relied on clinician-obtained cervical samples, self-sampling could represent an acceptable and valid alternative for epidemiologic studies in the developing world where low compliance with pelvic exams may lead to self-selection of participants and a reluctance to provide repeat samples for research studies. Further, self-sampling may provide an efficient way to monitor vaccine effectiveness in these high-risk populations.

Study Population

Participants are from the ongoing Rakai Community Cohort Study described in detail elsewhere (26). Briefly, the Rakai Community Cohort Study conducts annual home-based surveys, with comprehensive interviewer-administered questionnaires on sociodemographic, behavioral, and health characteristics, on an average population of 12,000 adults ages 15 to 49 years. Biological samples are collected following the interview. Venous blood is collected for HIV-1, syphilis, and herpes simplex type 2 serology. Starting in 1998, self-collected vaginal swabs were obtained during routine follow-up visits. Consenting women were asked to provide a self-collected vaginal sample at home, using the Digene sampler kit, with a Dacron swab. Women were asked to squat and insert a sterile 20-cm Dacron swab into the vagina up to the vault and to rotate the swab three times in the vaginal vault. The swab was then placed in 1 mL standard transport medium (Digene) and kept on ice until transported back to the field laboratories and frozen at -80°C. All samples were periodically shipped to Johns Hopkins University and stored at -80°C until the time of assay.

Although majority (>85%) of women provided self-collected vaginal swabs at each Rakai follow-up visit, for this longitudinal study, our analytic sample was composed of women who provided three to four self-collected swabs at consecutive follow-up visits.

The study was approved by the institutional review boards in Uganda and at Johns Hopkins University and Columbia University.

Laboratory Analyses

HIV Serology. The HIV status of each participant was determined using two enzyme immunoassays (Vironostika HIV-1; Organon Teknika and Cambridge Biotech), and discordant enzyme immunoassay results or new seroconversions were confirmed by Western blot (HIV-1 Western Blot; Bio-Merieux-Vitek).

HPV Detection and Genotyping. HPV status was determined using a two-stage approach. First, using hybrid capture 2 (HC2; Digene), all samples were screened for the presence of 13 carcinogenic HPV genotypes most clinically relevant to cervical premalignant and malignant lesions. Subsequently, all HC2 HPV-positive samples and 10% random sample of HC2 HPV-negative samples were genotyped using the PGMY09/11 L1 consensus primer PCR to amplify HPV DNA and reverse hybridization technique (Roche Molecular Systems) to detect 37 HPV genotypes (27, 28).

HPV DNA Testing by HC2. HC2 is a Food and Drug Administration approved, commercially available HPV test that collectively detects 13 carcinogenic HPV types (HPV 16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, 59, and 68) without distinguishing the HPV genotype. HC2 is a signal amplification assay that uses antibody capture of HPV DNA and RNA probe hybrids and chemiluminescent signal detection. The HC2 assay was done according to the manufacturer's instructions on a 100 μL aliquot of each specimen, a strategy that preserved sufficient sample volume for the additional PCR-based analyses.

HPV DNA Testing by Roche-PCR Line Blot Assay: DNA Extraction. HPV genotyping was done using the reverse line-blot assay as described previously (28). Briefly, 50 μL standard transport medium sample was denatured with 10× digestion buffer for 1 h at 65°C followed by heat denaturation at 95°C for 10 min. DNA was precipitated with ethanol and ammonium acetate at -20°C overnight. After centrifugation at 21,000 × g for 30 min at 4°C, DNA pellet was dried and resuspended in 25 μL Tris-EDTA and stored at -20°C (13).

DNA (5 μL) was amplified using the PGMY09/11 system that coamplifies HPV genotypes and human β-globin internal control target in a single reaction (13, 27, 28). Genotypes were discriminated by hybridizing the biotinylated PCR products to a probe array (RMS Linear Array) developed by Roche Molecular Systems, as described earlier, which can differentiate between 37 HPV types (27).

Appropriate quality-control samples were imbedded during the DNA purification, PCR, and genotype discrimination steps.

Statistical Analyses

Incidence of HR-HPV Infection. Genotype-specific HPV incident infection was defined as a positive test result for that genotype not present at study baseline. Subsequently, the 13 HR-HPV (16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, 59, and 68) genotypes most associated with cervical cancer were combined into one outcome variable of “incidence of new carcinogenic HPV” infection. Person-time methods were used to calculate the type-specific incidence density (29).

Univariate and multivariate Poisson regression was used to evaluate the effect of exposure variables on the risk of new HR-HPV infection. Incidence of any of the 13 HR-HPV was measured as the number of new infections in a given time interval divided by the population at risk during the interval. The log of the HPV incidence rate is assumed to be linearly associated with the risk factors, and the model parameters, after exponentiation, can be interpreted as relative incidence, which is similar to the relative risk. Behavioral and demographic factors measured at every visit were entered into the models as time-dependent covariates. When subjects missed study visits, the exposures were treated as missing. The independent contribution of each exposure variable, after adjustment for the effects of the other covariates, was expressed as the adjusted relative incidence rate or rate ratio (RR).

Clearance of HR-HPV Infection. Subjects included in the clearance analysis had to have been positive at least once for any of the 13 HR-HPV genotypes during their follow-up visits.

Methods for multiple failure-time data were used to evaluate the effect of covariates on the risk of type-specific HR-HPV clearance. The unit of analysis for the carcinogenic HPV clearance is the HPV infection. Hence, individuals infected with multiple HPV genotypes could clear one but not another type, and one or more clearance events could occur for the same individual. We first assigned clearance for each HPV genotype and then constructed an overall variable for combined clearance. We assumed clearance events to be unordered, which means that clearing one HPV type was independent on clearing another HPV type. We then used survival analysis methods to evaluate factors associated with type-specific HPV clearance. Time was measured from the first time a woman tested positive for a specific HPV type, and the event, type-specific HPV clearance, was defined as the first time a woman tested negative following a positive measurement. Those who never cleared the specific HPV-type infection were administratively censored at the end of the study follow-up. Because it is not reasonable to assume that clearance of each HPV type should have the same baseline hazard, the Cox model was stratified by each HPV type. Relative hazards <1 indicate lower clearance.

Incidence of Carcinogenic HPV Infection

One thousand seventy-nine women contributed 2,954 visits for the incidence analysis. The median number of visits per person was 3, with an interquartile range of 2 to 4 visits. There were 24 women with only 1 visit who were excluded from this analysis, leaving 1,055 women in the analytic sample.

When we compared the analytic sample with all women in the Rakai cohort (n = 6,520) at the HPV study baseline visit, we observed significant differences in certain characteristics suggesting that the analytic sample had a lower risk profile for HPV infection. Women in the analytic sample were slightly younger (median age, 28 years; interquartile range, 23-36) compared with those not in the analytic sample (median age, 29 years; interquartile range, 22-40; P = 0.005), less likely to be HIV positive compared with others (14% versus 21%; P < 0.0001), more likely to be currently married (75% versus 63%; P < 0.0001), and with fewer lifetime sexual partners (P < 0.0001; data not shown). Thus, the analytic sample had a lower risk profile than the general Rakai population.

Table 1 shows the univariate associations with new HR-HPV infections. Incidence was constant over time of observation, indicating that there were no temporal trends over the study duration. HR-HPV incidence was highest among younger women and decreased with age. There were no associations with education. Married women had a lower incidence of HR-HPV infection compared with single women. Incidence was lower among gravid women compared with nulligravid women. Incidence was lower among never-pregnant women; however, this was due to confounding by single marital status.

Table 1.

Factors associated with HR-HPV incidence: crude and adjusted rates

Person-yearsNew infectionsIncidence per 100 person-yearsCrude RR (95% CI)Adjusted RR (95% CI)
Demographics      
    Visit number      
        1 862 82 9.51  
        2 1,021 77 7.54 0.79 (0.58-1.08)  
        3 170 19 11.16 1.17 (0.71-1.93)  
    Age (y)      
        15-19 64 14.09 4.20 (1.88-9.70) 6.23 (2.39-16.25) 
        20-24 449 47 10.46 3.12 (1.72-5.66) 4.00 (2.07-7.71) 
        25-29 508 45 8.86 2.64 (1.45-4.81) 2.61 (1.36-5.03) 
        30-34 286 34 11.88 3.54 (1.90-6.60) 3.04 (1.55-5.95) 
        35-39 329 29 8.80 2.62 (1.39-4.96) 2.48 (1.26-4.90) 
        ≥40 417 14 3.36 
    Education      
        0 211 15 7.10  
        Primary 1,403 117 8.34 1.17 (0.69-2.01)  
        Secondary 353 38 10.78 1.51 (0.83-2.75)  
        Technical/professional 86 9.31 1.31 (0.56-3.09)  
    Currently married      
        No 521 58 11.13 
        Yes 1,533 120 7.83 0.70 (0.51-0.96) 0.62 (0.43-0.90) 
    Ever pregnant      
        No 46 17.47  
        Yes 2,008 170 8.47 0.48 (0.24-0.98)  
    Gravidity      
        0 44 18.28  
        1-5 1,161 113 9.74 0.53 (0.26-1.09)  
        ≥6 845 55 6.51 0.36 (0.17-0.75)  
Health status      
    HIV status      
        Negative 1,647 115 7.00 
        Positive 307 53 17.29 2.48 (1.79-3.43) 2.04 (1.43-2.91) 
        Missing 100 10 9.98 —  
    STD symptoms since last visit      
        0 968 73 7.55  
        1 385 38 9.87 1.31 (0.88-1.94)  
        ≥2 701 67 9.55 1.27 (0.91-1.76)  
    Current STD symptoms      
        0 1,622 141 8.69  
        1 239 18 7.53 0.87 (0.53-1.42)  
        ≥2 193 19 9.86 1.13 (0.70-1.83)  
    Health symptoms      
        0 680 55 8.09  
        1 978 79 8.07 1.00 (0.71-1.41)  
        ≥2 396 44 11.12 1.38 (0.93-2.04)  
Behaviors      
    Age at first sex      
        5-16 502 47 9.36 1.03 (0.73-1.45)  
        No sex/≥17 1,097 100 9.11  
    Lifetime sex partners      
        0/1 571 29 5.08 
        2 614 58 9.45 1.86 (1.19-2.91) 1.96 (1.21-3.18) 
        ≥3 870 91 10.46 2.06 (1.36-3.13) 2.13 (1.31-3.47) 
    Sex partners since last visit      
        0/1 258 14 5.43 
        >2 1,796 164 9.13 1.68 (0.97-2.90) 2.34 (1.25-4.38) 
    Perceived AIDS risk: partner      
        None/do not know 1,408 100 7.16 
        High 639 78 12.20 1.72 (1.28-2.31) 1.52 (1.11-2.08) 
    Perceived AIDS risk: self      
        None 1,297 92 7.09  
        Do not know 247 20 8.10 1.14 (0.70-1.85)  
        High 510 66 12.94 1.82 (1.33-2.50)  
    Self and partner drinking      
        Neither drank 704 62 8.81  
        Partner only 642 53 8.26 0.94 (0.65-1.35)  
        Self only 85 8.27 0.94 (0.43-2.05)  
        Both drank 608 55 9.05 1.03 (0.71-1.48)  
Person-yearsNew infectionsIncidence per 100 person-yearsCrude RR (95% CI)Adjusted RR (95% CI)
Demographics      
    Visit number      
        1 862 82 9.51  
        2 1,021 77 7.54 0.79 (0.58-1.08)  
        3 170 19 11.16 1.17 (0.71-1.93)  
    Age (y)      
        15-19 64 14.09 4.20 (1.88-9.70) 6.23 (2.39-16.25) 
        20-24 449 47 10.46 3.12 (1.72-5.66) 4.00 (2.07-7.71) 
        25-29 508 45 8.86 2.64 (1.45-4.81) 2.61 (1.36-5.03) 
        30-34 286 34 11.88 3.54 (1.90-6.60) 3.04 (1.55-5.95) 
        35-39 329 29 8.80 2.62 (1.39-4.96) 2.48 (1.26-4.90) 
        ≥40 417 14 3.36 
    Education      
        0 211 15 7.10  
        Primary 1,403 117 8.34 1.17 (0.69-2.01)  
        Secondary 353 38 10.78 1.51 (0.83-2.75)  
        Technical/professional 86 9.31 1.31 (0.56-3.09)  
    Currently married      
        No 521 58 11.13 
        Yes 1,533 120 7.83 0.70 (0.51-0.96) 0.62 (0.43-0.90) 
    Ever pregnant      
        No 46 17.47  
        Yes 2,008 170 8.47 0.48 (0.24-0.98)  
    Gravidity      
        0 44 18.28  
        1-5 1,161 113 9.74 0.53 (0.26-1.09)  
        ≥6 845 55 6.51 0.36 (0.17-0.75)  
Health status      
    HIV status      
        Negative 1,647 115 7.00 
        Positive 307 53 17.29 2.48 (1.79-3.43) 2.04 (1.43-2.91) 
        Missing 100 10 9.98 —  
    STD symptoms since last visit      
        0 968 73 7.55  
        1 385 38 9.87 1.31 (0.88-1.94)  
        ≥2 701 67 9.55 1.27 (0.91-1.76)  
    Current STD symptoms      
        0 1,622 141 8.69  
        1 239 18 7.53 0.87 (0.53-1.42)  
        ≥2 193 19 9.86 1.13 (0.70-1.83)  
    Health symptoms      
        0 680 55 8.09  
        1 978 79 8.07 1.00 (0.71-1.41)  
        ≥2 396 44 11.12 1.38 (0.93-2.04)  
Behaviors      
    Age at first sex      
        5-16 502 47 9.36 1.03 (0.73-1.45)  
        No sex/≥17 1,097 100 9.11  
    Lifetime sex partners      
        0/1 571 29 5.08 
        2 614 58 9.45 1.86 (1.19-2.91) 1.96 (1.21-3.18) 
        ≥3 870 91 10.46 2.06 (1.36-3.13) 2.13 (1.31-3.47) 
    Sex partners since last visit      
        0/1 258 14 5.43 
        >2 1,796 164 9.13 1.68 (0.97-2.90) 2.34 (1.25-4.38) 
    Perceived AIDS risk: partner      
        None/do not know 1,408 100 7.16 
        High 639 78 12.20 1.72 (1.28-2.31) 1.52 (1.11-2.08) 
    Perceived AIDS risk: self      
        None 1,297 92 7.09  
        Do not know 247 20 8.10 1.14 (0.70-1.85)  
        High 510 66 12.94 1.82 (1.33-2.50)  
    Self and partner drinking      
        Neither drank 704 62 8.81  
        Partner only 642 53 8.26 0.94 (0.65-1.35)  
        Self only 85 8.27 0.94 (0.43-2.05)  
        Both drank 608 55 9.05 1.03 (0.71-1.48)  

HIV-positive women were at a significantly elevated risk of new HR-HPV infection compared with HIV-negative women [RR, 2.48; 95% confidence interval (95% CI), 1.79-3.43]. Self-reported numbers of current sexually transmitted disease (STD) symptoms or self-reported general health symptoms were not associated with new HR-HPV infection (Table 1).

Age at first sex was not associated with new HR-HPV infection (Table 1), whereas HR-HPV incidence was significantly elevated among women who reported 2 and ≥3 lifetime sexual partners [RR (95% CI), 1.86 (1.19-2.91) and 2.06 (1.36-3.13), respectively] compared with women with 0 or 1 lifetime partners. Incidence was marginally elevated among women with ≥2 sex partners in the past year. Women who perceived themselves or their partner to be at high risk for HIV/AIDS infection were more likely to have new HR-HPV infection. Drinking alcohol was not associated with incident HR-HPV.

Risk factors associated with new HR-HPV infections in the univariate analyses were evaluated in multivariate models (Table 1). After adjustment, risk of incident HR-HPV was significantly increased with younger age, HIV infection, unmarried women, 2 and ≥3 lifetime sex partners, and ≥2 partners in the past year. Women who perceived their recent partner to be at high risk for HIV/AIDS infection were more likely to have new HR-HPV infection.

Table 2 shows incidence rates for the 13 HR-HPV individually and combined. HPV-51 had the highest incidence rate followed by HPV-16, whereas HPV-31 and HPV-35 had the lowest incidence rates. A total of 178 new HR-HPV infections were observed during 2,054 person-years of observation, reflecting an incidence rate of 8.7 per 100 person-years. Stratifying by HIV serostatus, among 9 of the 13 HR-HPV types, the risk of a new infection was significantly greater among HIV-positive compared with HIV-negative women (types 16, 18, 35, 45, 51, 56, 58, 59, and 68). The highest incidence RR comparing HIV-positive with HIV-negative women was observed for HPV-58. Overall, the risk of a new HPV infection was significantly higher among HIV-positive women (17.3 per 100 person-years) compared with HIV-negative women (7.0 per 100 person-years; P < 0.0001).

Table 2.

Type-specific HR-HPV incidence, overall and HIV stratified

HPV typeOverall
HIV negative
HIV positive
HIV+/HIV-
Incident casesPerson-years at riskRate per 100 person-yearsIncident casesPerson-years at riskRate per 100 person-yearsIncident casesPerson-years at riskRate per 100 person-yearsCrude incidence RR (95% CI)Adjusted incidence RR (95% CI)*
16 32 2,106 1.5 20 1,689 1.2 10 312 3.2 2.71 (1.13-6.06) 3.09 (1.39-6.85) 
18 18 2,131 0.9 11 1,692 0.7 333 2.1 3.23 (1.06-9.13) 3.19 (1.17-8.71) 
31 10 2,157 0.5 1,712 0.5 336 0.3 0.60 (0.01-4.09) 0.56 (0.07-4.60) 
33 17 2,124 0.8 14 1,691 0.8 327 0.6 0.74 (0.08-3.22) 0.51 (0.12-2.30) 
35 13 2,126 0.6 1,702 0.4 315 1.9 4.63 (1.29-16.09) 3.88 (1.25-12.03) 
39 18 2,141 0.9 12 1,698 0.7 336 1.2 1.68 (0.40-5.56) 2.01 (0.60-6.68) 
45 23 2,128 1.1 13 1,694 0.8 328 2.7 3.58 (1.35-9.04) 3.41 (1.40-8.29) 
51 37 2,100 1.8 20 1,678 1.2 16 315 5.1 4.26 (2.06-8.66) 3.58 (1.80-7.11) 
52 22 2,090 1.1 16 1,678 1.0 307 1.3 1.47 (0.33-4.24) 1.17 (0.38-3.60) 
56 20 2,131 0.9 12 1,699 0.7 325 1.9 2.61 (0.80-7.53) 2.80 (1.00-7.86) 
58 14 2,125 0.7 1,688 0.4 329 2.1 5.99 (1.72-21.56) 6.15 (1.94-19.55) 
59 22 2,147 1.0 12 1,703 0.7 338 2.4 3.36 (1.19-8.94) 3.16 (1.24-8.09) 
68 18 2,128 0.9 10 1,684 0.6 336 2.1 3.51 (1.13-10.21) 2.99 (1.10-8.20) 
Any 13 HPV 178 2,054 8.7 115 1,647 7.0 53 307 17.3 2.47 (1.75-3.45) 2.30 (1.57-3.37) 
HPV typeOverall
HIV negative
HIV positive
HIV+/HIV-
Incident casesPerson-years at riskRate per 100 person-yearsIncident casesPerson-years at riskRate per 100 person-yearsIncident casesPerson-years at riskRate per 100 person-yearsCrude incidence RR (95% CI)Adjusted incidence RR (95% CI)*
16 32 2,106 1.5 20 1,689 1.2 10 312 3.2 2.71 (1.13-6.06) 3.09 (1.39-6.85) 
18 18 2,131 0.9 11 1,692 0.7 333 2.1 3.23 (1.06-9.13) 3.19 (1.17-8.71) 
31 10 2,157 0.5 1,712 0.5 336 0.3 0.60 (0.01-4.09) 0.56 (0.07-4.60) 
33 17 2,124 0.8 14 1,691 0.8 327 0.6 0.74 (0.08-3.22) 0.51 (0.12-2.30) 
35 13 2,126 0.6 1,702 0.4 315 1.9 4.63 (1.29-16.09) 3.88 (1.25-12.03) 
39 18 2,141 0.9 12 1,698 0.7 336 1.2 1.68 (0.40-5.56) 2.01 (0.60-6.68) 
45 23 2,128 1.1 13 1,694 0.8 328 2.7 3.58 (1.35-9.04) 3.41 (1.40-8.29) 
51 37 2,100 1.8 20 1,678 1.2 16 315 5.1 4.26 (2.06-8.66) 3.58 (1.80-7.11) 
52 22 2,090 1.1 16 1,678 1.0 307 1.3 1.47 (0.33-4.24) 1.17 (0.38-3.60) 
56 20 2,131 0.9 12 1,699 0.7 325 1.9 2.61 (0.80-7.53) 2.80 (1.00-7.86) 
58 14 2,125 0.7 1,688 0.4 329 2.1 5.99 (1.72-21.56) 6.15 (1.94-19.55) 
59 22 2,147 1.0 12 1,703 0.7 338 2.4 3.36 (1.19-8.94) 3.16 (1.24-8.09) 
68 18 2,128 0.9 10 1,684 0.6 336 2.1 3.51 (1.13-10.21) 2.99 (1.10-8.20) 
Any 13 HPV 178 2,054 8.7 115 1,647 7.0 53 307 17.3 2.47 (1.75-3.45) 2.30 (1.57-3.37) 

NOTE: Slight discrepancy in incident case numbers and person-time is due to visits with missing HIV results.

*

Adjusted for lifetime number of partners and age.

Clearance of HPV Infections

One hundred ninety-five HPV-positive women contributed 389 observations (the unit of analysis is type-specific HPV infections) and 465 person-times for analysis of viral clearance. Two hundred forty-nine of the infections cleared over the follow-up period (53.55%).

Table 3 presents findings of factors associated with HPV clearance. Women ages >30 years were less likely to clear infections compared with women ages 15 to 29 years. Women with higher education were more likely to clear infection compared with those with no education, as were currently married women. Ever-pregnancy and gravidity were not associated with clearance in univariate analyses.

Table 3.

Factors associated with HR-HPV clearance: crude and adjusted rates

No. subjectsTime at riskClearance rateCrude clearance rate (95% CI)Adjusted clearance rate (95% CI)
Demographics      
    Age (y)      
        15-29 219 250 61.60 
        ≥30 170 215 44.18 0.77 (0.63-0.96) 0.84 (0.68-1.02) 
    Education      
        None 57 74 43.24 
        Primary 257 306 52.94 1.19 (0.87-1.64) 1.31 (0.92-1.87) 
        Secondary and more 75 85 64.71 1.49 (1.07-2.08) 1.58 (1.08-2.30) 
    Currently married      
        No 156 198 44.94  
        Yes 233 267 59.93 1.31 (1.04-1.65)  
    Ever pregnant      
        No 14 16 56.25  
        Yes 375 449 53.45 0.84 (0.62-1.13)  
    Gravidity      
        1-5 251 296 55.74  
        ≥6 138 169 49.70 0.97 (0.77-1.24)  
Health status      
    HIV status      
        Negative 194 219 66.21 
        Positive 175 224 39.73 0.70 (0.56-0.87) 0.75 (0.59-0.96) 
    STD symptoms since last visit      
        0 137 166 51.81  
        1 51 61 55.74 1.04 (0.79-1.36)  
        ≥2 201 238 54.20 1.05 (0.86-1.29)  
    Current STD symptoms      
        0 246 288 56.25  
        1 63 78 50.00 0.93 (0.71-1.21)  
        ≥2 80 99 48.48 0.91 (0.69-1.18)  
    Health symptoms      
        0 109 129 56.59  
        1 145 172 55.23 1.03 (0.82-1.28)  
        ≥2 135 164 49.39 0.93 (0.73-1.18)  
    Multiple HPV infections      
        No 122 141 59.57  
        Yes 267 324 50.93 0.90 (0.73-1.11)  
    HC2 relative light unit categories      
        >0-10 161 190 61.05 
        ≥10-100 141 162 52.47 0.76 (0.53-1.08) 0.86 (0.70-1.07) 
        ≥100 87 113 42.48 0.68 (0.44-1.06) 0.76 (0.58-0.99) 
Behaviors      
    Age at first sex      
        5-15 199 233 55.79 1.05 (0.84-1.31)  
        ≥16 185 227 50.22  
    Lifetime no. sex partners      
        1 68 76 68.42 
        2 101 124 50.81 0.77 (0.60-1.00) 0.74 (0.57-0.97) 
        ≥3 220 265 50.57 0.79 (0.63-0.99) 0.92 (0.72-1.17) 
    Sex partners since last visit      
        0 44 59 42.37  
        1 323 379 55.94 1.18 (0.87-1.59)  
        ≥2 22 27 44.44 0.88 (0.46-1.66)  
    Perceived AIDS risk: partner      
        None 66 77 59.74  
        Do not know 158 186 53.76 0.95 (0.74-1.23)  
        High 164 201 51.24 0.87 (0.67-1.11)  
    Perceived AIDS risk: self      
        None 170 191 59.69  
        Do not know 68 81 49.38 0.90 (0.69-1.17)  
        High 151 193 49.22 0.89 (0.72-1.07)  
    Self and partner drinking      
        Neither drank 138 161 57.14  
        Partner only 125 154 53.25 0.92 (0.71-1.19)  
        Self only 17 20 55.00 0.95 (0.63-1.45)  
        Both drank 108 129 49.61 0.91 (0.72-1.16)  
No. subjectsTime at riskClearance rateCrude clearance rate (95% CI)Adjusted clearance rate (95% CI)
Demographics      
    Age (y)      
        15-29 219 250 61.60 
        ≥30 170 215 44.18 0.77 (0.63-0.96) 0.84 (0.68-1.02) 
    Education      
        None 57 74 43.24 
        Primary 257 306 52.94 1.19 (0.87-1.64) 1.31 (0.92-1.87) 
        Secondary and more 75 85 64.71 1.49 (1.07-2.08) 1.58 (1.08-2.30) 
    Currently married      
        No 156 198 44.94  
        Yes 233 267 59.93 1.31 (1.04-1.65)  
    Ever pregnant      
        No 14 16 56.25  
        Yes 375 449 53.45 0.84 (0.62-1.13)  
    Gravidity      
        1-5 251 296 55.74  
        ≥6 138 169 49.70 0.97 (0.77-1.24)  
Health status      
    HIV status      
        Negative 194 219 66.21 
        Positive 175 224 39.73 0.70 (0.56-0.87) 0.75 (0.59-0.96) 
    STD symptoms since last visit      
        0 137 166 51.81  
        1 51 61 55.74 1.04 (0.79-1.36)  
        ≥2 201 238 54.20 1.05 (0.86-1.29)  
    Current STD symptoms      
        0 246 288 56.25  
        1 63 78 50.00 0.93 (0.71-1.21)  
        ≥2 80 99 48.48 0.91 (0.69-1.18)  
    Health symptoms      
        0 109 129 56.59  
        1 145 172 55.23 1.03 (0.82-1.28)  
        ≥2 135 164 49.39 0.93 (0.73-1.18)  
    Multiple HPV infections      
        No 122 141 59.57  
        Yes 267 324 50.93 0.90 (0.73-1.11)  
    HC2 relative light unit categories      
        >0-10 161 190 61.05 
        ≥10-100 141 162 52.47 0.76 (0.53-1.08) 0.86 (0.70-1.07) 
        ≥100 87 113 42.48 0.68 (0.44-1.06) 0.76 (0.58-0.99) 
Behaviors      
    Age at first sex      
        5-15 199 233 55.79 1.05 (0.84-1.31)  
        ≥16 185 227 50.22  
    Lifetime no. sex partners      
        1 68 76 68.42 
        2 101 124 50.81 0.77 (0.60-1.00) 0.74 (0.57-0.97) 
        ≥3 220 265 50.57 0.79 (0.63-0.99) 0.92 (0.72-1.17) 
    Sex partners since last visit      
        0 44 59 42.37  
        1 323 379 55.94 1.18 (0.87-1.59)  
        ≥2 22 27 44.44 0.88 (0.46-1.66)  
    Perceived AIDS risk: partner      
        None 66 77 59.74  
        Do not know 158 186 53.76 0.95 (0.74-1.23)  
        High 164 201 51.24 0.87 (0.67-1.11)  
    Perceived AIDS risk: self      
        None 170 191 59.69  
        Do not know 68 81 49.38 0.90 (0.69-1.17)  
        High 151 193 49.22 0.89 (0.72-1.07)  
    Self and partner drinking      
        Neither drank 138 161 57.14  
        Partner only 125 154 53.25 0.92 (0.71-1.19)  
        Self only 17 20 55.00 0.95 (0.63-1.45)  
        Both drank 108 129 49.61 0.91 (0.72-1.16)  

NOTE: Hazard rates <1 indicate lower clearance rate.

HIV-positive women were significantly less likely to clear infection compared with HIV-negative women (RR, 0.70; 95% CI, 0.56-0.87). Clearance was not associated with symptoms suggestive of STD, AIDS, opportunistic infections, or infection with multiple HPV types (Table 3). Women with greater HPV shedding as indicated by the higher relative light units (based on HC2 assays) were less likely to clear their infection compared with those in the lowest relative light units.

Clearance rates were lower in women who reported 2 or ≥3 lifetime partners compared with women with 1 lifetime sex partner (RR, 0.77 and 0.79, respectively). There were no associations with age at first sex, recent sex partners, AIDS perception, or drinking.

In multivariate analysis (Table 3), HIV-positive women were less likely to clear their infection (adjusted RR, 0.75; 95% CI, 0.59-0.96) compared with HIV-negative women, as were women ages >30 years compared with women ages ≤29 years (adjusted RR, 0.84; 95% CI, 0.68-1.02). Clearance was also reduced in women reporting ≥2 lifetime partners and among women with higher HPV viral burden. Women with secondary or higher education were more likely to clear their infection. Stratified analysis by HIV status showed that among HIV-positive women, higher HPV viral burden, and higher number of lifetime sexual partners and age (marginal effect) were associated with lower clearance. However, among HIV-negative women, age was associated with a marginal decrease in clearance, but education was associated with higher clearance (data not shown).

Table 4 presents clearance rates for the 13 different HR-HPV types. Highest clearance rates were observed for HPV-18 and HPV-68 followed by HPV-39. Lowest clearance occurred for HPV-31 and HPV-35 followed by HPV-16. Stratifying by HIV status, for almost all viral types (with the exception of HPV-56), we observed lower clearance rate among HIV-positive compared with HIV-negative women. HPV-33 achieved statistical significance after adjusting for age and lifetime number of sexual partners.

Table 4.

Clearance rate for the 13 HR-HPV types (unit of analysis is infections)

HPV typeOverall
HIV negative
HIV positive
HIV−/HIV+
No. clearedPerson-years at riskClearance rateNo. clearedPerson-yearsClearance rateNo. clearedPerson-yearsClearance rateCrude clearance (95% CI)Adjusted clearance (95% CI)*
16 23 60 38.3 10 19 52.6 10 38 26.3 1.74 (0.81-3.73) 1.61 (0.74-3.51) 
18 22 27 81.5 14 16 87.5 10 80.0 1.10 (0.75-1.62) 1.22 (0.85-1.77) 
31 12 16.7 50.0 10 10.0 4.00 (0.41-39.24) 2.57 (0.60-11.00) 
33 19 42 45.2 14 22 63.6 16 25.0 2.32 (0.79-6.78) 3.96 (1.88-8.31) 
35 12 43 27.9 11 45.4 32 21.9 2.13 (0.67-6.77) 2.60 (0.91-7.47) 
39 17 26 65.4 11 15 73.3 44.4 1.34 (0.71-2.51) 1.44 (0.69-2.99) 
45 20 31 64.5 12 16 75 11 45.4 1.88 (0.86-3.97) 1.33 (0.57-3.08) 
51 36 59 61.0 21 33 63.6 14 24 58.3 1.13 (0.69-1.82) 1.11 (0.66-1.88) 
52 35 64 54.7 20 29 69.0 14 34 41.2 1.68 (0.90-3.13) 1.70 (0.90-3.24) 
56 20 31 64.5 14 50.0 11 15 73.3 0.67 (0.33-1.36) 0.37 (0.18-0.76) 
58 18 30 60.0 11 15 73.3 14 42.9 1.57 (0.77-3.22) 1.76 (0.91-3.39) 
59 14 42.9 57.1 33.3 1.67 (0.43-6.52) 1.00 (0.29-3.43) 
68 19 26 73.1 15 20 75.0 60.0 1.23 (0.55-2.72) 1.41 (0.66-3.01) 
Any 13 HPV 249 465 53.5 145 219 66.2 89 224 39.7 1.42 (1.14-1.78) 1.40 (1.10-1.77) 
HPV typeOverall
HIV negative
HIV positive
HIV−/HIV+
No. clearedPerson-years at riskClearance rateNo. clearedPerson-yearsClearance rateNo. clearedPerson-yearsClearance rateCrude clearance (95% CI)Adjusted clearance (95% CI)*
16 23 60 38.3 10 19 52.6 10 38 26.3 1.74 (0.81-3.73) 1.61 (0.74-3.51) 
18 22 27 81.5 14 16 87.5 10 80.0 1.10 (0.75-1.62) 1.22 (0.85-1.77) 
31 12 16.7 50.0 10 10.0 4.00 (0.41-39.24) 2.57 (0.60-11.00) 
33 19 42 45.2 14 22 63.6 16 25.0 2.32 (0.79-6.78) 3.96 (1.88-8.31) 
35 12 43 27.9 11 45.4 32 21.9 2.13 (0.67-6.77) 2.60 (0.91-7.47) 
39 17 26 65.4 11 15 73.3 44.4 1.34 (0.71-2.51) 1.44 (0.69-2.99) 
45 20 31 64.5 12 16 75 11 45.4 1.88 (0.86-3.97) 1.33 (0.57-3.08) 
51 36 59 61.0 21 33 63.6 14 24 58.3 1.13 (0.69-1.82) 1.11 (0.66-1.88) 
52 35 64 54.7 20 29 69.0 14 34 41.2 1.68 (0.90-3.13) 1.70 (0.90-3.24) 
56 20 31 64.5 14 50.0 11 15 73.3 0.67 (0.33-1.36) 0.37 (0.18-0.76) 
58 18 30 60.0 11 15 73.3 14 42.9 1.57 (0.77-3.22) 1.76 (0.91-3.39) 
59 14 42.9 57.1 33.3 1.67 (0.43-6.52) 1.00 (0.29-3.43) 
68 19 26 73.1 15 20 75.0 60.0 1.23 (0.55-2.72) 1.41 (0.66-3.01) 
Any 13 HPV 249 465 53.5 145 219 66.2 89 224 39.7 1.42 (1.14-1.78) 1.40 (1.10-1.77) 
*

Adjusted for lifetime number of partners and age.

Data on the epidemiology of HPV from resource poor countries are sparse, mainly due to lack of longitudinal studies. A constraint on performing these studies has been recruiting and retaining participants to studies, which require pelvic examination for collection of cervical specimens. This analysis based on self-collected vaginal swabs showed that HR-HPV incidence and clearance from a rural community in Uganda mirrors studies of the epidemiology of HPV based on cervical samples. These findings have important implications for future longitudinal studies in settings where resource constraints limit pelvic exams and, most importantly, where women are reluctant to participate in studies that rely on pelvic exams for data collection (25, 30-32).

We were comfortable using self-collected vaginal swab to determine the epidemiology of HPV, because in an earlier paper we validated HPV DNA detection from self-collected vaginal swabs to physician-collected cervical swabs showing very good agreement between vaginal and cervical DNA detection (25).

Similar to other studies, we found younger age and higher number of recent and lifetime sex partners to be significant risk factors for new HPV infection (1, 33). Younger age reflects onset of sexual activity as well as limited acquired immunity. The association of recent numbers of sexual partners with incident detection of HR-HPV suggests that part of the new infections are due to sexual behavior; however, the association of higher numbers of lifetime sexual partners with incident infection suggests that part of the new infections may be due to reactivated past infections. Married women had lower risk of new HPV infection compared with single women and women with high perception of their partner's AIDS risk were more likely to be at risk of new HPV infection. This may be a marker of the participant or her partner's sexual behaviors that were not directly measured or were misreported.

HIV-positive women were twice as likely to have a new HR-HPV infection in accord with published studies from the United States and Europe (34, 35). With the exception of HPV-31 and HPV-33, the incident RR for the other 11 HR-HPV ranged from 1.5 to 6 times higher among HIV-positive women compared with HIV-negative women (HPV-52 and HPV-58, respectively). The incidence RR for HPV-16 was 3.09 among HIV-positive compared with HIV-negative women, similar to findings by Strickler et al. (36), suggesting that HPV-16 may be weakly associated with immune status in HIV-positive women compared with other types and be better at evading the immune system, which can lead to higher persistence and its stronger association with cervical cancer.

HPV-31, HPV-35, and HPV-16 had the lowest clearance rates as reported in other studies (37, 38). Similar to other studies (35, 34, 39), we observed lower HPV clearance rates among HIV-positive women compared with HIV-negative women. We found HPV viral burden to be associated with lower clearance, and this was highly significant among HIV-positive women. This may be an indication of the damage by HIV to the immune system and its inability to perform normally. Furthermore, higher HPV viral load is also associated with worse cytopathology and hence could be an indication of higher-grade cervical lesions among HIV-positive women. However, although we did not have pathology on all women, among those with colposcopy (and biopsy if it was indicated), we observed more lesions among HIV-positive women than HIV-negative women. Additionally, we found that age affected HPV clearance: women ages >30 years were less likely to clear an infection compared with women ages <30 years, consistent with findings of other studies (1, 4, 34). We also found that women with more lifetime sex partners were less likely to clear their HPV infection. A similar finding was observed by Shew et al. (40) among adolescents. This is suggestive of multiple past exposures or reactivated infections acquired earlier, which may be less prone to clearance.

Interpretations of our results should take into account some limitations. Our application of the HC2 assay to all samples, followed by genotyping of all HC2-positive samples and 10% randomly selected HC2-negative samples, assumed that HC2 detected all 13 HPV types in its probe with reasonable accuracy. We found 8% HC2-negative samples were PCR positive, which could have resulted in an underestimation of the true incidence and clearance of HR-HPV. Differential virus quantity used in HC2 and line-blot assays or variations in the laboratory techniques could be possible explanation for the discordant results between HC2 and line-blot assay. To minimize differences due to laboratory practices and ensure quality control, one person (M.S.) performed all laboratory assays, and negative control cell lines were included for DNA isolation and positive and negative controls were included in the PCR and genotyping procedures. Reasons for differences between HC2 and line-blot warrant further investigation.

Further, all factors found to be associated with HPV infection in other studies were not measured in this study. We could not assess the effect of CD4 cell counts and HIV viral load on incidence and clearance of HPV. Additionally, because no women were using oral contraceptives or smoking cigarettes in this sample, we could not assess their effects on incidence or clearance.

Specific to the incidence and clearance analysis, some women may have acquired or lost their infection in between-study visits, which could lead to underestimation of incidence or overestimation of duration of infection. With regard to the latter, our objective was not to report on the duration of infection but rather to characterize women who cleared their infection.

Because the epidemiology of HPV as measured by incidence and clearance in this representative rural Ugandan population was similar to studies that relied on cervical sampling, recruitment and retention of women for natural history studies of HPV may be increased by self-collection as an alternate to clinician-collected cervical specimens. In settings where obtaining the necessary data for the natural history of HPV poses a challenge for researchers, self-collection can provide a reliable method for obtaining samples and may be adopted in trials to increase compliance and reduce costs associated with conducting pelvic exams. Moreover, determining the effectiveness of vaccines in resource poor and HIV endemic countries is dependent on knowledge of natural history and epidemiology of HPV and can provide an important tool for future research or vaccine evaluation.

No potential conflicts of interest were disclosed.

Grant support: Fogarty HIV-Associated Malignancies, New and Minority Investigator Award, Fogarty AIDS Training and Research Program, Johns Hopkins University and NIH/National Institute of Allergy and Infectious Diseases STI training grant T32AI50056 (M. Safaeian).

Note: Presented at 23rd International Papillomavirus Conference, Prague, Czech Republic, September 2006.

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.

We thank the Roche Diagnostics for providing the reagents and line blots for the genotyping assay and Digene for providing the HC2 probe B kits at reduced costs.

1
Ho GY, Bierman R, Beardsley L, Chang CJ, Burk RD. Natural history of cervicovaginal papillomavirus infection in young women.
N Engl J Med
1998
;
338
:
423
–8.
2
Winer RL, Lee SK, Hughes JP, Adam DE, Kiviat NB, Koutsky LA. Genital human papillomavirus infection: incidence and risk factors in a cohort of female university students.
Am J Epidemiol
2003
;
157
:
218
–26.
3
Woodman CB, Collins S, Winter H, et al. Natural history of cervical human papillomavirus infection in young women: a longitudinal cohort study.
Lancet
2001
;
357
:
1831
–6.
4
Hildesheim A, Schiffman MH, Gravitt PE, et al. Persistence of type-specific human papillomavirus infection among cytologically normal women.
J Infect Dis
1994
;
169
:
235
–40.
5
Molano M, Van den BA, Plummer M, et al. Determinants of clearance of human papillomavirus infections in Colombian women with normal cytology: a population-based, 5-year follow-up study.
Am J Epidemiol
2003
;
158
:
486
–94.
6
Schiffman M. Integration of human papillomavirus vaccination, cytology, and human papillomavirus testing.
Cancer
2007
;
111
:
145
–53.
7
Sankaranarayanan R, Black RJ, Swaminathan R, Parkin DM. An overview of cancer survival in developing countries. IARC Sci Publ 135-173; 1998.
8
Gordis L. The epidemiologic approach to evaluation of screening programs. Epidemiology. W.B. Saunders Co.; 2000. p. 258–76.
9
Baldwin S, Santos C, Mendez BE, et al. Comparison of type-specific human papillomavirus data from self and clinician directed sampling.
Gynecol Oncol
2005
;
97
:
612
–7.
10
Belinson JL, Qiao YL, Pretorius RG, et al. Shanxi Province cervical cancer screening study. II. Self-sampling for high-risk human papillomavirus compared to direct sampling for human papillomavirus and liquid based cervical cytology.
Int J Gynecol Cancer
2003
;
13
:
819
–26.
11
Dannecker C, Siebert U, Thaler CJ, Kiermeir D, Hepp H, Hillemanns P. Primary cervical cancer screening by self-sampling of human papillomavirus DNA in internal medicine outpatient clinics.
Ann Oncol
2004
;
15
:
863
–9.
12
Garcia F, Barker B, Santos C, et al. Cross-sectional study of patient- and physician-collected cervical cytology and human papillomavirus.
Obstet Gynecol
2003
;
102
:
266
–72.
13
Gravitt PE, Lacey JV, Jr., Brinton LA, et al. Evaluation of self-collected cervicovaginal cell samples for human papillomavirus testing by polymerase chain reaction.
Cancer Epidemiol Biomarkers Prev
2001
;
10
:
95
–100.
14
Harper DM, Noll WW, Belloni DR, Cole BF. Randomized clinical trial of PCR-determined human papillomavirus detection methods: self-sampling versus clinician-directed-biologic concordance and women's preferences.
Am J Obstet Gynecol
2002
;
186
:
365
–73.
15
Hillemanns P, Kimmig R, Huttemann U, Dannecker C, Thaler CJ. Screening for cervical neoplasia by self-assessment for human papillomavirus DNA.
Lancet
1999
;
354
:
1970
.
16
Kahn JA, Slap GB, Huang B, et al. Comparison of adolescent and young adult self-collected and clinician-collected samples for human papillomavirus.
Obstet Gynecol
2004
;
103
:
952
–9.
17
Lack N, West B, Jeffries D, et al. Comparison of non-invasive sampling methods for detection of HPV in rural African women.
Sex Transm Infect
2005
;
81
:
239
–41.
18
Lorenzato FR, Singer A, Ho L, et al. Human papillomavirus detection for cervical cancer prevention with polymerase chain reaction in self-collected samples.
Am J Obstet Gynecol
2002
;
186
:
962
–8.
19
Palmisano ME, Gaffga AM, Daigle J, et al. Detection of human papillomavirus DNA in self-administered vaginal swabs as compared to cervical swabs.
Int J STD AIDS
2003
;
14
:
560
–7.
20
Petignat P, Hankins C, Walmsley S, et al. Self-sampling is associated with increased detection of human papillomavirus DNA in the genital tract of HIV-seropositive women.
Clin Infect Dis
2005
;
41
:
527
–34.
21
Salmeron J, Lazcano-Ponce E, Lorincz A, et al. Comparison of HPV-based assays with Papanicolaou smears for cervical cancer screening in Morelos State, Mexico.
Cancer Causes Control
2003
;
14
:
505
–12.
22
Sellors JW, Lorincz AT, Mahony JB, et al. Comparison of self-collected vaginal, vulvar and urine samples with physician-collected cervical samples for human papillomavirus testing to detect high-grade squamous intraepithelial lesions.
CMAJ
2000
;
163
:
513
–8.
23
Wright TC, Jr., Denny L, Kuhn L, Pollack A, Lorincz A. HPV DNA testing of self-collected vaginal samples compared with cytologic screening to detect cervical cancer.
JAMA
2000
;
283
:
81
–6.
24
2006 Report on the global AIDS epidemic. 2006. Internet communication.
25
Safaeian M, Kiddugavu M, Gravitt PE, et al. Comparability of self-collected vaginal swabs and physician-collected cervical swabs for detection of human papillomavirus infections in Rakai, Uganda.
Sex Transm Dis
2007
;
34
:
429
–36.
26
Wawer MJ, Gray RH, Sewankambo NK, et al. A randomized, community trial of intensive sexually transmitted disease control for AIDS prevention, Rakai, Uganda.
AIDS
1998
;
12
:
1211
–25.
27
Gravitt PE, Peyton CL, Alessi TQ, et al. Improved amplification of genital human papillomaviruses.
J Clin Microbiol
2000
;
38
:
357
–61.
28
Gravitt PE, Peyton CL, Apple RJ, Wheeler CM. Genotyping of 27 human papillomavirus types by using L1 consensus PCR products by a single-hybridization, reverse line blot detection method.
J Clin Microbiol
1998
;
36
:
3020
–7.
29
Breslow NE, Day NE, editors. The design and analysis of cohort studies. Lyon: IARC; 1987.
30
Brown CL, Screening patterns for cervical cancer: how best to reach the unscreened population.
J Natl Cancer Inst Monogr
1996
;
21
:
7
–11.
31
Byrd TL, Chavez R, Wilson KM. Barriers and facilitators of cervical cancer screening among Hispanic women.
Ethn Dis
2007
;
17
:
129
–34.
32
Mutyaba T, Faxelid E, Mirembe F, Weiderpass E. Influences on uptake of reproductive health services in Nsangi community of Uganda and their implications for cervical cancer screening.
Reprod Health
2007
;
4
:
4
.
33
Franco EL, Villa LL, Sobrinho JP, et al. Epidemiology of acquisition and clearance of cervical human papillomavirus infection in women from a high-risk area for cervical cancer.
J Infect Dis
1999
;
180
:
1415
–23.
34
Ahdieh L, Klein RS, Burk R, et al. Prevalence, incidence, and type-specific persistence of human papillomavirus in human immunodeficiency virus (HIV)-positive and HIV-negative women.
J Infect Dis
2001
;
184
:
682
–90.
35
Sun XW, Kuhn L, Ellerbrock TV, Chiasson MA, Bush TJ, Wright TC, Jr. Human papillomavirus infection in women infected with the human immunodeficiency virus.
N Engl J Med
1997
;
337
:
1343
–9.
36
Strickler HD, Palefsky JM, Shah KV, et al. Human papillomavirus type 16 and immune status in human immunodeficiency virus-seropositive women.
J Natl Cancer Inst
2003
;
95
:
1062
–71.
37
Munoz N, Mendez F, Posso H, et al. Incidence, duration, and determinants of cervical human papillomavirus infection in a cohort of Colombian women with normal cytological results.
J Infect Dis
2004
;
190
:
2077
–87.
38
Richardson H, Kelsall G, Tellier P, et al. The natural history of type-specific human papillomavirus infections in female university students.
Cancer Epidemiol Biomarkers Prev
2003
;
12
:
485
–90.
39
Strickler HD, Burk RD, Fazzari M, et al. Natural history and possible reactivation of human papillomavirus in human immunodeficiency virus-positive women.
J Natl Cancer Inst
2005
;
97
:
577
–86.
40
Shew ML, Fortenberry JD, Tu W, et al. Association of condom use, sexual behaviors, and sexually transmitted infections with the duration of genital human papillomavirus infection among adolescent women.
Arch Pediatr Adolesc Med
2006
;
160
:
151
–6.