Home | About | Journals | Submit | Contact Us | Français |

**|**PLoS One**|**v.7(1); 2012**|**PMC3265500

Formats

Article sections

Authors

Related links

PLoS One. 2012; 7(1): e30736.

Published online 2012 January 24. doi: 10.1371/journal.pone.0030736

PMCID: PMC3265500

Julia Kravchenko,^{1,}^{*} Igor Akushevich,^{2} Staci L. Sudenga,^{3} Craig M. Wilson,^{3} Emily B. Levitan,^{3} and Sadeep Shrestha^{3,}^{*}

Beatriz G. J. Grinsztejn, Editor^{}

Instituto de Pesquisa Clínica Evandro Chagas/Fundação Oswaldo Cruz, Brazil

Conceived and designed the experiments: JK IA SS. Performed the experiments: JK IA SS. Analyzed the data: JK IA SS SLS. Contributed reagents/materials/analysis tools: IA EBL CMW. Wrote the paper: JK IA SS. Provided input in developments of models: EBL. Provided critical reviews of the manuscript: SLS EBL CMW.

Received 2011 July 12; Accepted 2011 December 22.

Copyright Kravchenko et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

This article has been cited by other articles in PMC.

HIV-1-positive patients clear the human papillomavirus (HPV) infection less frequently than HIV-1-negative. Datasets for estimating HPV clearance probability often have irregular measurements of HPV status and risk factors. A new transitional probability-based model for estimation of probability of HPV clearance was developed to fully incorporate information on HIV-1-related clinical data, such as CD4 counts, HIV-1 viral load (VL), highly active antiretroviral therapy (HAART), and risk factors (measured quarterly), and HPV infection status (measured at 6-month intervals).

Data from 266 HIV-1-positive and 134 at-risk HIV-1-negative adolescent females from the Reaching for Excellence in Adolescent Care and Health (REACH) cohort were used in this study. First, the associations were evaluated using the Cox proportional hazard model, and the variables that demonstrated significant effects on HPV clearance were included in transitional probability models. The new model established the efficacy of CD4 cell counts as a main clearance predictor for all type-specific HPV phylogenetic groups. The 3-month probability of HPV clearance in HIV-1-infected patients significantly increased with increasing CD4 counts for HPV16/16-like (p<0.001), HPV18/18-like (p<0.001), HPV56/56-like (p=0.05), and low-risk HPV (p<0.001) phylogenetic groups, with the lowest probability found for HPV16/16-like infections (21.60±1.81% at CD4 level 200 cells/mm^{3}, p<0.05; and 28.03±1.47% at CD4 level 500 cells/mm^{3}). HIV-1 VL was a significant predictor for clearance of low-risk HPV infections (p<0.05). HAART (with protease inhibitor) was significant predictor of probability of HPV16 clearance (p<0.05). HPV16/16-like and HPV18/18-like groups showed heterogeneity (p<0.05) in terms of how CD4 counts, HIV VL, and HAART affected probability of clearance of each HPV infection.

This new model predicts the 3-month probability of HPV infection clearance based on CD4 cell counts and other HIV-1-related clinical measurements.

HIV-1-positive women clear HPV infections 4–10 times more slowly than HIV-1-negative, and HIV-1-infected patients with CD4+ T-lymphocytes cell count (CD4) <200 cells/mm^{3} show the slowest clearance [1], [2]. Understanding the role of immunosuppression in risk of persistence of sexually transmitted human papillomavirus (HPV) infection, a main risk factor for cervical intraepithelial neoplasia and a central etiologic agent of cervical cancer, and clarifying how this risk is modified by other factors (such as co-infections, antiretroviral therapy, and behavioral factors) is important for optimization of follow-up strategy [3]. Information from longitudinal studies about factors that affect the probability of type-specific HPV clearance can be used to estimate the impact of cervical cancer interventions. However, this is often problematic because HPV persistence is loosely defined as detection of the same HPV type at two or more subsequent visits, ranging from 2 months to 7 years [4], [5], [6], [7], [8], [9], and its probability depends on intervals between the tests. Further, the analysis is complicated by the possibility of co-infections with multiple HPV types as well as by varying length of intervals between missing measurements of HPV status. Analytical methods that can fully utilize real-life heterogeneous data, specifically, clinical data that are unevenly scheduled (for example, at 6 months for HPV and at 3-month intervals for HIV clinical data) could be useful in studies of the factors having an impact on probability of clearance of HPV infection.

Here, we describe a transition probability model for studying the relationship between immune status (based on CD4 cell count) and probability of HPV clearance in HIV-infected patients. The Reaching for Excellence in Adolescent Care and Health (REACH) dataset used in this study planned for measurement of HIV-1 status every 3 months and HPV status every 6 months. As in other cohorts, this study had missing and irregular visit measurements: only 82% of time intervals between measurements of HPV statuses were performed as scheduled (i.e., every 6 months), while other HPV tests were done at 3-, 9-, 12-, or other-month intervals. Analyzing these data with standard techniques would require multiple assumptions about the definition of HPV clearance and censoring time for each individual for time-dependent predictors; as a result, part of the dataset would not be utilized. We have developed transition probability-based models and applied an HPV/HIV co-infected cohort to estimate 3-month HPV clearance probabilities while maximizing all available data for both HIV and HPV in the estimations.

Data on 266 HIV-1-positive and 134 high-risk HIV-1-negative adolescent females from the REACH cohort were analyzed. The REACH study design and methods of quarterly follow up with HIV-1 testing, immunophenotyping, HIV-1 RNA viral load (VL) and collection of biological specimens, demographic and behavioral factors, and other clinical data, along with incidence and prevalence of HPV infections, have been previously described [2], [10], [11]. Briefly, between 1996 and 2000, adolescents aged 12–19 years who were HIV-1-positive and comparable at-risk HIV-1-negative persons were recruited into a longitudinal study at 15 clinical sites in the United States. HIV-1-related clinical data and risk factors were measured every 3 months. At enrollment and every 6 months thereafter, cervical lavage samples were tested for HPV infection by MY09/MY11/HMB01-based PCR and for 30 HPV type-specific probes with a chemiluminescent dot blot procedure [12]. PCR-based HPV data were classified as follows: negative; positive for specific HPV types; or “positive, type unknown” (when the sample was positive for the generic probe but not for specific HPV type). PCR amplification of a human β-globin gene segment was used for internal DNA quality control, and samples negative for this assay were excluded from the analysis. For certain types of analyses, HPVs were categorized according to phylogenetic patterns [2], [13] into: 1) 16/16-like (16, 31, 52, 58, 67); 2) 18/18-like (18, 39, 45, 59, 26, 51); 3) 56/56-like (56, 53, 66); and 4) low-risk (6, 11, 42, 44, 54, 40, 13, 32, 62, 72, 2, 57, 55). The highly active antiretroviral therapy (HAART) was defined as a combination of two nucleoside reverse transcriptase inhibitors and either a protease inhibitor (PI) or a non-nucleoside reverse transcriptase inhibitor, or a zidovudine/lamivudine combination regimen plus another antiretroviral drug. Data on antiretroviral therapy were obtained through interviews and chart reviews for current prescriptions, and adherence data were obtained through interviews as previously described [14].

All the participants of the study provided written informed consent in the parent study at the Adolescent Medicine HIV/AIDS Research Network for the REACH project, and the UAB Institutional Review Board approved this sub-study. The parent study and this sub-study conformed to the procedures for informed consent (parental permission was obtained wherever required) approved by institutional review boards at all sponsoring organizations and to human-experimentation guidelines set forth by the United States Department of Health and Human Services.

*Cox proportional hazard* regression was used to test which covariates have significant effects on probability of HPV clearance for each of four phylogenetic HPV groups. First, an univariable regression analysis was used for a wide spectrum of variables including demographic characteristics, clinical exams, antiretroviral therapy, behavioral factors, and coinfections (as listed in Table 1). From this, a subset of predictors with significant effects was selected and considered in multivariable analysis. For this analysis, SAS PROC PHREG (Cary, NC) was used.

The model developed in the present paper falls under the general category of transitional Markov models and is referred to as a *transitional probability-based model*. These kinds of approaches, also known as regressive or conditional models, are used in epidemiology for analyses of dependent binary observations [15], [16], [17]. The model allows for estimating the probability of changing HPV status in patients with specific HPV type, i.e., a chance to clear HPV infection during the follow-up period. To account for HIV-1-related clinical data such as CD4 counts, HIV-1 VL, and other HIV-1-related factors, which were measured every 3 months, HPV infection status was reconstructed for the same time interval considering transition probabilities. The transition probabilities were referred to as , where was initial state of type-specific HPV infection ( corresponded to absence, and to presence of type-specific HPV infection), and and corresponded to the status of type-specific HPV infection at the end of a 3-month period. Vector denoted the set of most influential predictors of HPV clearance probability, such as CD4 count, HIV-1 VL, HAART, and HPV type [2], [18]. For example, the probability corresponded to the situations where HPV infection observed at a recent visit was cleared in 3 months.

With this model, the information about missed measurements at odd visits can be reconstructed using the previous and forthcoming measurements (see Figure 1A). The transition probability between two subsequent visits with measured HPV status could be presented as , where *i* and *j* describe the HPV status at the first and third visits, respectively. The status at a second visit is unknown, and, therefore, the sum over two possible intermediate states contributes to the observed transition probability between two subsequent even visits. Parameters and in the above formula denote the sets of respective predictors for transitions between first-to-second and second-to-third visits, respectively. The likelihood function is the product over all transfers with known HPV status. If data are taken exactly according to the cohort's measurement design, the likelihood is

(1)

Here, covers all individuals in the dataset, —all transitions between states with measured HPV virus status represented by indices and in individual_{,} and and are the vectors of predictors measured at the beginning of time period of respective transition. The dependence of transition probabilities on predictors are modeled in logistic regression style as:

(2)

where intercepts and refer to logarithms of odds of changing the type-specific HPV status for zeroth values of predictors, and parameters and describe the effects of respective predictors [19].

If the number of missed HPV status is 0 (zeroth) or varies from two to five, the likelihood function could be generalized by summing over all intermediate states. For example, for transitions with three missed HPV statuses (see Figure 1B), the contribution to the likelihood function could be presented as:

(3)

where , , and are unmeasured (0 or 1) HPV status in three intermediate states. The set of observed transfers with measured HPV statuses for a specific type of HPV virus or for HPV group is the input dataset for likelihood maximization. Note that CD4 count and HIV-1 VL also could have missing measurements; however, the appearance of these missing variables is not related to the study design, so we assume that they are missing at random, and any standard approach for filling missing data at random can be applied (e.g., imputation with the mean conditional on observed values of other variables or simply linear interpolation). Because in the REACH cohort, the fractions of missing values of predictors of HPV clearance were relatively small (about 2% for CD4 count in HIV-1-positive patients, about 5% for HIV-1 VL, and about 1% for HAART), no notable impact of specific scheme of their filling was expected. The reported results were obtained using a linear interpolation to fill missing values in the predictors (but not HPV) status.

To test the different hypotheses regarding the possible effects of different potential predictors on probability of HPV clearance, a two-stage approach was designed. The two-stage approach complemented the advantages of methods used at each stage such as nonparametric estimates of hazard ratios in the Cox model and opportunities for evaluating 3-month probabilities for models (1) – (3). The comparison of the results obtained from both approaches allows for validating the properties of the new method. At the first stage, the effects were evaluated using the Cox proportional hazard model, and then the variables that demonstrated significant effects on HPV clearance in the Cox analysis were included in the logistic type models (substantively specified below). Note, the results can be obtained only under specific assumptions that are necessary to identify the times of HPV incidence and clearance/censoring. Moreover, the specific assumptions also are required to decide how to treat missing data on unknown HPV status. There is no consensus in literature about the choice of specific assumptions in these situations. Based on the published studies, our set of assumptions included: 1) time of HPV incidence was 90 days before the first exam with positive HPV status; 2) patient was removed when a time period between measurements exceeded 450 days; 3) when no time period between measurements exceeded 450 days, the clearance time was the time of the first (of the two subsequent) negative exams, ignoring any missing visits; 4) when the last exam was positive, the censoring time was the time of this exam plus 180 days; and 5) when the last exam was negative (after the positive), the censoring time was the time of this exam.

The approach used at the second stage (i.e., the formalism based on eq. (1–3)) does not require these assumptions. All logistic-type models described the probability of a 3-month HPV clearance for different predictors using the base equation 2; their parameters were estimated using the likelihood described in equation 1 and generalized in the style of equation 3. The CD4 count parameter was selected to be tested by the model as a main predictor of HPV clearance in HIV-1-positive patients based on multiple studies [2], [20], [21]—the base model (M1) had it as a main parameter and an argument of the exponent for the model in equation 2 was . Then, additional models were developed, extending the base model as follows:

M2 —includes the effect of the presence of HIV-1 infection on intercept of , where is the binary variable characterizing the HIV-1 seropositivity () or HIV-1 seronegativity (). While the base model M1 is applied separately for HIV-1-positive and HIV-1-negative patients, the M2 as well as the M3 models are designed to investigate whether and how specifically the effect of CD4 counts on HPV clearance probability is influenced by the presence of HIV-1 infection, which may have an effect on HPV clearance beyond the CD4 counts effect (e.g., this effect could be further compared with the patients with other immunodeficiencies such as those with organ transplants or inherited immune disorders).

M3 —extends the M2 model by including the possible interference of the effects of HIV-1 seropositivity and CD4 count on HPV clearance . This model evaluates the modifying effect of CD4 on HPV clearance by the presence of HIV-1-infection.

M4 —in HIV-1-infected patients, it introduces CD4 count as a squared parameter in (the exponent argument is ) to test the assumption about non-linearity (squared CD4) in interrelations between CD4 count and HPV clearance probability in HIV-1-positive patients.

M5 —in HIV-1-positive patients, the effect of CD4 count is described as a piecewise-linear interpolation with an arbitrary set , of nodes (i.e., CD4 count values at which the linear functions are joined):

where is the indicator function and ; parameters referring to the logarithms of odds of changing the HPV status for . This model tests whether different shapes of interrelations between CD4 counts and probability of HPV clearance are possible in HIV-1-positive patients at a very low (<200 cells/mm^{3}), low (200–499 cells/mm^{3}), and normal range of CD4 counts.

M6 —includes the logarithm of HIV-1 VL () to evaluate the effect of HIV-1 VL on probability of HPV clearance in HIV-1-positive patients controlling the level of CD4 count.

M7 —investigates the effects of HAART that include PIs on the intercept , where is a binary variable of HAART; when HAART (with PI) was applied at the time of the visit. This model tests whether HAART with PIs may have an additional effect on interrelations between CD4 count and probability of HPV clearance in HIV-1-infected patients.

The nonlinear optimization techniques as implemented in PROC NLP in SAS 9.2 (Cary, NC) were used for the likelihood maximization in all these models to estimate a 3-month probability of HPV clearance. For the majority of calculations, the intercept and the effects of the predictors were chosen as model parameters and estimated with the standard errors (SEs). Because of the functional relation between parameter *β* (which describes the effects of predictor of clearance) and the clearance probability for a certain value of predictor, the latter can be used as a model parameter instead of *β*. Estimation of this model using the Proc NLP allows for evaluating its standard error.

The main characteristics of studied patients are listed in Table 1. The average age of HIV-1-infected adolescent females were 16.8 (±1.1) years old, and 78.6% of them were African Americans. On average, they had about 11 HIV-1 status-related visits/examinations per patient (quarterly); however, HPV data were collected only biannually, so during half of these visits. Table 2 provides detailed descriptions of incident and prevalent HPV infections: among oncogenic HPVs, HPV16, 31, 52, and 67 (HPV16/16-like group), HPV59 and 26 (HPV18/18-like group), and HPV56 and 53 (HPV56/56-like group) were more often (p<0.05) registered in HIV-1-infected than HIV-1-negative patients.

The results of the Cox proportional hazard analysis for HIV-1-infected adolescent females are shown in Table 3. No significant HRs were obtained for such parameters as having <6 or >15 lifetime sexual partners, ever being a smoker, or being infected with *Trichomonas vaginalis* or *Neisseria gonorrhea*. CD4 count was a consistent predictor for clearance of HPV16/16-like, 18/18-like, and low-risk groups. For these HPV groups, a significantly higher probability of HPV clearance was at CD4 levels higher than 500 cells/mm^{3} (compared with those at <500 cells/mm^{3}); the difference was especially pronounced for CD4 categorical cutout at 200 cells/mm^{3}. For HIV-1 VL, significant effect on HPV clearance (for all except HPV56/56-like group) was observed only without CD4 count as a second parameter (in the univariable analysis). While being on HAART was a significant predictor of HPV16/16-like clearance only, being on HAART with PI significantly increased probability of clearance of 16/16-like, 18/18-like, and low-risk HPV infection. Also, being infected with *Chlamydia trachomatis* could be a positive predictor for low-risk HPV clearance.

Sensitivity analysis was performed to check the stability of HRs estimated by the Cox model. In many cases, the HRs essentially shifted when assumptions changed: e.g., the effect of CD4 being >200 cells/mm^{3} changed for HPV16/16-like infection from 1.68 to 1.86, and from 2.53 to 3.07 (still remaining significant) when 365 days were used in assumptions #2 and #3 instead of 450 days. Another example is the changes of the HAART effect from 1.77 to 2.03 for HPV16/16-like infection, which remained highly significant while using 270 days instead of 450 in these assumptions.

The results obtained from the basic model and from the models describing the effects of HIV-1 VL and HAART (with PIs) are shown in Table 4, and the results obtained from all models (M1–M7) are presented in Table S1. The probability of HPV clearance in HIV-1-infected patients increased with increasing CD4 level for all HPV groups: parameter *β _{11}* (the log odds ratio, describing the effect of CD4 count on HPV clearance) was 1.15 (0.27) for HPV16/16-like (p<0.001), 1.58 (0.36) for HPV18/18-like (p<0.001), 0.72 (0.38) for HPV56/56-like (p=0.059), and 1.5 (0.41) for low-risk HPVs (p<0.001). HPV16/16-like infection was least likely to clear at low CD4 cell count (<200 cells/mm

When two oncogenic groups of HPV infections—HPV16/16-like and HPV18/18-like—were analyzed for each HPV type separately, both groups showed heterogeneity in terms of how probability of type-specific HPV clearance was affected by CD4 counts, HIV VL, and HAART with PI (see Table 6): while lower clearance probability was registered for HPV16/16-like than for the 18/18-like group (p<0.05). Interestingly, HPV16 and HPV18 alone had an equal chance to be cleared at all CD4 levels examined: i.e., 20.29(3.76)% and 18.16(4.04)% at CD4 level 200 cells/mm^{3}, 26.63(3.19)% and 28.66(4.0)% at CD4 level 400 cells/mm^{3}, and 34.12(3.49)% and 42.12(6.24)% at CD4 level 600 cells/mm^{3}, for HPV16 and HPV18, respectively). The effect of HIV-1 VL was significant for clearance of HPV58 (16-like group) and HPV59 (18-like group), and the effect of HAART (with PIs) was significant for HPV16 clearance. In average, HPV67 had a higher probability to be cleared than the other 16/16-like HPV types: e.g., 64.11 (15.72)% vs. 25.76 (1.48)% at CD4 level 400 cells/mm^{3} (p<0.05). Likewise, HPV18 had the lowest probability to be cleared (18.16 (4.04)%), and HPV26 the highest (42.12 (8.39)%) than other 18/18-like HPV types at CD4 level 200 cells/mm^{3} compared to the group at average (27.40 (2.38)%), p<0.05.

The longitudinal studies of HPV infection are important for determining the covariates and outcomes associated with HPV persistence, which leads to the development of cancer. Traditionally, the rate of HPV clearance are usually compared in HIV-1-positive patient subgroups based on baseline CD4 counts (such as <200, 200–500, and >500 cells/mm^{3}) using Kaplan-Meier curves and Cox proportional hazards models [1], [2], [22], [23]. The new method developed in this study is based on the logistic-type model and allows for prediction of future HPV status, conditional on its current status and the measurements of factors that are potential predictors of an HPV clearance event; i.e., CD4 count was considered to be a main predictor with HIV-1 VL and HAART (with PI) as additional predictors. The model parameters were estimated by maximizing the likelihood function constructed as the product over transfers with known HPV statuses (measured every 6 months) and HIV-1-related covariates (measured every 3 months). Similar to Kong et al. [17], our methodology is based on conditional probabilities that take into account multiple correlations between individual outcomes measured longitudinally. However, the developed method extends the opportunities of approaches by Kong as well as several other researchers [1], [2], [24], [25] by reconstructing between-the-measurements HPV statuses (i.e., presence or absence of HPV infection). This approach allows for inclusion of the whole longitudinal dataset, thus increasing the accuracy of prediction of probability of HPV clearance without making multiple assumptions about how the time of incidence, clearance, and censoring events could be reconstructed (as it is required for Kaplan-Meier and Cox analyses). The shorter intervals accessible in our method allow for taking into account the dynamics of potential predictors, which could change quickly (such as CD4 count, HIV-1 VL, and HAART regimen). This model allows for calculating the clearance probability with subsequent confirmation in another 3 months——by transferring the probabilities such as , thus corresponding to the routine definition of HPV clearance when the absence of HPV type-specific infection is required for two subsequent visits. Opposite to the Cox model, in which HRs are estimated for time-dependent covariates, the developed approach allows us to estimate the transition probability and evaluate its standard errors. Since the developed model provides the hazard function for probability of HPV infection clearance, the respective survival function and characteristics of time-to-clearance distribution also can be evaluated: e.g., time to clearance (in months) could be estimated as , and a median of HPV clearance time as . The approaches utilizing the generalized estimating equation (GEE)—they take into account the mutual correlations in clearance of different HPV types and modeling mixed effects, allowing individuals to have their own characteristics (i.e., distributed in a population)—could be used to further enrich the developed base model; Xue et al. [26] recently reviewed the series of approaches that can be used for similar generalizations.

In both approaches, a transitional probability-based model and Cox regression model, CD4 count was a significant predictor of clearance of all phylogenetic HPV groups in HIV-1-infected adolescent females; also, certain effects of HAART (with PI) on clearance of HPV16/16-like and HPV18/18-like infections were observed. However, while in the Cox model, being HIV-1-infected had a minor effect only on HPV56/56-like clearance, in the transitional probability model, this factor was a significant predictor of clearance of HPV16/16-like, HPV56/56-like, and low-risk HPVs.

In immunodeficient patients, the mechanisms by which immune deficiency increases the risk of persistence of HPV infection are still poorly understood: some alterations in dendritic antigen-presenting cells, Langerhans cells, and macrophages function, as well as a deficient cytotoxic lymphocyte response to E6 and E7 proteins, might be the contributing factors [27], [28]. The results obtained from our study about the role of CD4 in HPV clearance corroborate previous reports from the REACH cohort, as well as several other studies on adult HIV-infected females [1], [2], [29], [30], [31], [32], [33], [34]. However, there is no agreement about the role CD4 play in clearance of individual types of HPVs; e.g., in several studies on both HIV-1-negative and HIV-1-positive females, it has been shown that HPV16 infection has a lower probability of clearance than other HPV types, possibly due to its greater ability to escape immunologic surveillance [35], [36], [37], while other studies did not demonstrate such a difference [38]. In our study, a lower clearance probability was registered for the HPV16/16-like than for the 18/18-like group, while HPV16 had an equal-with-HPV18 chance to be cleared at both pathologic and normal CD4 counts. The observed heterogeneity of phylogenetic groups of HPV infection in terms of a probability of HPV clearance may depend not only on CD4 counts and other predictors measured at current time (such as HIV-1 VL and HAART with PI), but also from the history of HPV type-specific infection (e.g., from the time since HPV acquisition, which is an unobserved variable), depending on a prevalent or incident type-specific HPV infection. Further analysis could be performed using non-Markov approaches to model unobserved time since HPV acquisition.

HIV infection, independent of CD4 count, has also been suggested to be a predictor of persistence of HPV infection in HIV-1-positive women. This may imply an alternate mechanism besides CD4, e.g., via alteration of the cytokine response to HPV infection in the cervical mucus [2], [39], [40], [41], [42]. In our study, being HIV-1-positive affected the probability of clearance of HPV16/16-like, 18/18-like, and low-risk HPVs. In the REACH cohort–based study by Moscicki et al. [2], when only subjects with normal CD4 counts (i.e., ≥500 cells/mm^{3}) were considered, the multivariable regression showed high significance of HIV status as an independent predictor of HPV clearance event (HR=1.60, p=0.012).

Currently, prognostic importance of high HIV-1 VL for HPV clearance is not clear, but it likely increases the risk of persistence of HPV infection at low CD4 cell counts [32], [43]. In our study, HIV-1 VL could affect the clearance of low-risk HPVs and certain oncogenic HPVs (e.g., HPV58 and HPV59). The apparent impact of HAART on HPV incidence, clearance, and persistence also is not clear [11], [23], [32], [44], [45]. In our study, when HAART was analyzed taking into account its PI component, a significant effect was observed for HPV16, and minor effects were observed for HPV16/16-like, HPV18/18-like, and HPV59. In vitro studies have shown that specific PIs inhibit the ability of HPV16 E6 to degrade p53 and selectively kill E-6-dependent cervical carcinoma cells [46]. Previous crossover analyses in REACH suggested no significant effect of HAART on HPV clearance [11]; however, the effect of PI was not examined, as it is incorporated in this new method. These results require further investigation with longer follow up and more detailed analysis of HAART/PI history and dose/exposure.

The observation on *C. trachomatis* increasing probability of clearance of low-risk HPV falls in with the results from animal studies about potential role of interferon-γ as local “protector” against other (i.e., non-*Chlamydia*) infections [47], [48]. Oncogenic HPV types could be strong enough to avoid this mechanism; recently, it has been speculated that *C. trachomatis* could have effect on oncogenic HPV types [49].

The results of this study have several limitations. While the prevalence and incidence of HPV infections among HIV-1-positive adolescents in the REACH study is high [11], some of the associations may have been limited by the relatively smaller sample size of HIV-1-negative individuals along with the lower HPV infection rate. Due to the populations served at the REACH recruitment sites, young African-American women were a significant proportion of the population; therefore, the results may not be fully generalized to other populations. Also, the interrelations described in this study were obtained on a cohort of young adolescent girls with relatively short histories of HIV-1 infections, who are generally healthy and whose immune response to the infection may differ from older women; for example, it has been shown in several studies that older age was associated with higher risk of HPV persistence in both HIV-infected and HIV-uninfected women [24], [50]. Regarding the approach, the simple version of the model was intentionally chosen as a base model, resulting in some limitations; e.g., there was no distinction between the effects of incident infection and re-infection, no correlations between clearance of distinct HPV types in one individual were modeled, and the time after the incidence was not explicitly represented. Due to the two-step design of the study, some variables which were statistically insignificant were not included into the second step of the analysis thus potentially compromising the robustness of the model. Nevertheless, the limitations can be overcome by the extensions of the proposed approach using approaches specifically developed for analyses of HPV clearance [26] and those that were successfully used in other related research areas, e.g., g-formula [16] or a (binomial) stochastic process model [51], [52], [53].

In summary, our new model estimates a probability for HPV clearance of type-specific HPV groups at a 3-month period by coordinating uneven time scales of measurements on biannual HPV status and other quarterly HIV-1-related clinical data and risk factors.

HPV type-specific characteristics obtained from the M1 – M7 models tested in the REACH cohort of HIV-1-positive and HIV-1-negative adolescent girls enrolled in 1996-2000.

(DOC)

Click here for additional data file.^{(91K, doc)}

The authors thank the investigators and staff [listed in *J Adolesc Health* 2001;29 (suppl): 5 to 6] of the Adolescent Medicine HIV/AIDS Research Network (1994–2001) and the youth who participated in the REACH project for their valuable contributions.

**Competing Interests: **The authors have declared that no competing interests exist.

**Funding: **This work was supported in part by a career developmental award (SS) from the SPORE in Cervical Cancer (5P50CA098252-05). EBL was supported by the Agency for Healthcare Research and Quality (K12HS019465). SLS was supported by the Cancer Prevention and Control Training Program (5R25CA047888). The REACH study was supported by the National Institute of Child Health and Human Development (U01-HD32830). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

1. Ahdieh L, Munoz A, Vlahov D, Trimble CL, Timpson LA, et al. Cervical neoplasia and repeated positivity of human papillomavirus infection in human immunodeficiency virus-seropositive and -seronegative women. Am J Epidemiol. 2000;151:1148–1157. [PubMed]

2. Moscicki AB, Ellenberg JH, Farhat S, Xu J. Persistence of human papillomavirus infection in HIV-infected and -uninfected adolescent girls: risk factors and differences, by phylogenetic type. J Infect Dis. 2004;190:37–45. [PubMed]

3. Schiffman MH, Castle P. Epidemiologic studies of a necessary causal risk factor: human papillomavirus infection and cervical neoplasia. J Natl Cancer Inst. 2003;95:E2. [PubMed]

4. Woodman CB, Collins SI, Young LS. The natural history of cervical HPV infection: unresolved issues. Nat Rev Cancer. 2007;7:11–22. [PubMed]

5. Ho GY, Burk RD, Klein S, Kadish AS, Chang CJ, et al. Persistent genital human papillomavirus infection as a risk factor for persistent cervical dysplasia. J Natl Cancer Inst. 1995;87:1365–1371. [PubMed]

6. Kulmala SM, Shabalova IP, Petrovitchev N, Syrjanen KJ, Gyllensten UB, et al. Type-specific persistence of high-risk human papillomavirus infections in the New Independent States of the former Soviet Union Cohort Study. Cancer Epidemiol Biomarkers Prev. 2007;16:17–22. [PubMed]

7. Kjaer SK, van den Brule AJ, Paull G, Svare EI, Sherman ME, et al. Type specific persistence of high risk human papillomavirus (HPV) as indicator of high grade cervical squamous intraepithelial lesions in young women: population based prospective follow up study. Bmj. 2002;325:572. [PMC free article] [PubMed]

8. Nobbenhuis MA, Helmerhorst TJ, van den Brule AJ, Rozendaal L, Voorhorst FJ, et al. Cytological regression and clearance of high-risk human papillomavirus in women with an abnormal cervical smear. Lancet. 2001;358:1782–1783. [PubMed]

9. Kotloff KL, Wasserman SS, Russ K, Shapiro S, Daniel R, et al. Detection of genital human papillomavirus and associated cytological abnormalities among college women. Sex Transm Dis. 1998;25:243–250. [PubMed]

10. Rogers AS, Futterman DK, Moscicki AB, Wilson CM, Ellenberg J, et al. The REACH Project of the Adolescent Medicine HIV/AIDS Research Network: design, methods, and selected characteristics of participants. J Adolesc Health. 1998;22:300–311. [PubMed]

11. Shrestha S, Sudenga SL, Smith JS, Bachmann LH, Wilson CM, et al. The impact of highly active antiretroviral therapy on prevalence and incidence of cervical human papillomavirus infections in HIV-positive adolescents. BMC Infect Dis. 2010;10:295. [PMC free article] [PubMed]

12. Moscicki AB, Ellenberg JH, Vermund SH, Holland CA, Darragh T, et al. Prevalence of and risks for cervical human papillomavirus infection and squamous intraepithelial lesions in adolescent girls: impact of infection with human immunodeficiency virus. Arch Pediatr Adolesc Med. 2000;154:127–134. [PubMed]

13. Van Ranst M, Kaplan JB, Burk RD. Phylogenetic classification of human papillomaviruses: correlation with clinical manifestations. J Gen Virol. 1992;73(Pt 10):2653–2660. [PubMed]

14. Murphy DA, Wilson CM, Durako SJ, Muenz LR, Belzer M. Antiretroviral medication adherence among the REACH HIV-infected adolescent cohort in the USA. AIDS Care. 2001;13:27–40. [PubMed]

15. Bonney GE. Logistic regression for dependent binary observations. Biometrics. 1987;43:951–973. [PubMed]

16. Taubman SL, Robins JM, Mittleman MA, Hernan MA. Intervening on risk factors for coronary heart disease: an application of the parametric g-formula. Int J Epidemiol. 2009;38:1599–1611. [PMC free article] [PubMed]

17. Kong X, Archer KJ, Moulton LH, Gray RH, Wang MC. Parametric frailty models for clustered data with arbitrary censoring: application to effect of male circumcision on HPV clearance. BMC Medical Research Methodology. 2010;10:40. [PMC free article] [PubMed]

18. Lehtinen M, Luukkaala T, Wallin KL, Paavonen J, Thoresen S, et al. Human papillomavirus infection, risk for subsequent development of cervical neoplasia and associated population attributable fraction. J Clin Virol. 2001;22:117–124. [PubMed]

19. Hosmer DW, Lemeshow S. New York: Wiley; 2000. Applied Logistic Regression.

20. Harris TG, Burk RD, Palefsky JM, Massad LS, Bang JY, et al. Incidence of cervical squamous intraepithelial lesions associated with HIV serostatus, CD4 cell counts, and human papillomavirus test results. Jama. 2005;293:1471–1476. [PubMed]

21. Clark RA, Brandon W, Dumestre J, Pindaro C. Clinical manifestations of infection with the human immunodeficiency virus in women in Louisiana. Clin Infect Dis. 1993;17:165–172. [PubMed]

22. Shrestha S, Wang C, Aissani B, Wilson CM, Tang J, et al. Interleukin-10 gene (IL10) Polymorphisms and Human Papillomavirus Clearance among Immunosuppressed Adolescents. Cancer Epidemiology Biomarkers and Prevention. 2007;16:1–7. [PubMed]

23. Khanna N. HAART use in women with HIV and influence on cervical intraepithelial neoplasia: a clinical opinion. J Low Genit Tract Dis. 2002;6:111–115. [PubMed]

24. Ahdieh L, Klein RS, Burk R, Cu-Uvin S, Schuman P, 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–690. [PubMed]

25. Williams DA. The analysis of binary responses from toxicological experiments involving reproduction and teratogenicity. Biometrics. 1975;31:949–952. [PubMed]

26. Xue X, Gange SJ, Zhong Y, Burk RD, Minkoff H, et al. Marginal and mixed-effects models in the analysis of human papillomavirus natural history data. Cancer Epidemiol Biomarkers Prev. 19:159–169. [PMC free article] [PubMed]

27. Nakagawa M, Stites DP, Farhat S, Sisler JR, Moss B, et al. Cytotoxic T lymphocyte responses to E6 and E7 proteins of human papillomavirus type 16: relationship to cervical intraepithelial neoplasia. The Journal of infectious diseases. 1997;175:927–931. [PubMed]

28. Sobhani I, Walker F, Roudot-Thoraval F, Abramowitz L, Johanet H, et al. Anal carcinoma: incidence and effect of cumulative infections. Aids. 2004;18:1561–1569. [PubMed]

29. Sun XW, Kuhn L, Ellerbrock TV, Chiasson MA, Bush TJ, et al. Human papillomavirus infection in women infected with the human immunodeficiency virus. N Engl J Med. 1997;337:1343–1349. [PubMed]

30. Moscicki AB, Ellenberg JH, Crowley-Nowick P, Darragh TM, Xu J, et al. Risk of high-grade squamous intraepithelial lesion in HIV-infected adolescents. J Infect Dis. 2004;190:1413–1421. [PubMed]

31. Vermund SH, Kelley KF, Klein RS, Feingold AR, Schreiber K, et al. High risk of human papillomavirus infection and cervical squamous intraepithelial lesions among women with symptomatic human immunodeficiency virus infection. Am J Obstet Gynecol. 1991;165:392–400. [PubMed]

32. Denny L, Boa R, Williamson AL, Allan B, Hardie D, et al. Human papillomavirus infection and cervical disease in human immunodeficiency virus-1-infected women. Obstet Gynecol. 2008;111:1380–1387. [PubMed]

33. Sun XW, Ellerbrock TV, Lungu O, Chiasson MA, Bush TJ, et al. Human papillomavirus infection in human immunodeficiency virus-seropositive women. Obstet Gynecol. 1995;85:680–686. [PubMed]

34. Fennema JS, van Ameijden EJ, Coutinho RA, van den Hoek AA. HIV, sexually transmitted diseases and gynaecologic disorders in women: increased risk for genital herpes and warts among HIV-infected prostitutes in Amsterdam. Aids. 1995;9:1071–1078. [PubMed]

35. Molano M, Van den Brule A, Plummer M, Weiderpass E, Posso H, 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–494. [PubMed]

36. Strickler HD, Palefsky JM, Shah KV, Anastos K, Klein RS, et al. Human papillomavirus type 16 and immune status in human immunodeficiency virus-seropositive women. J Natl Cancer Inst. 2003;95:1062–1071. [PubMed]

37. Schiffman M, Herrero R, Desalle R, Hildesheim A, Wacholder S, et al. The carcinogenicity of human papillomavirus types reflects viral evolution. Virology. 2005;337:76–84. [PubMed]

38. Franco EL, Villa LL, Sobrinho JP, Prado JM, Rousseau MC, 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–1423. [PubMed]

39. Koshiol JE, Schroeder JC, Jamieson DJ, Marshall SW, Duerr A, et al. Time to clearance of human papillomavirus infection by type and human immunodeficiency virus serostatus. Int J Cancer. 2006;119:1623–1629. [PubMed]

40. Rowhani-Rahbar A, Hawes SE, Sow PS, Toure P, Feng Q, et al. The impact of HIV status and type on the clearance of human papillomavirus infection among Senegalese women. J Infect Dis. 2007;196:887–894. [PubMed]

41. Crowley-Nowick PA, Ellenberg JH, Vermund SH, Douglas SD, Holland CA, et al. Cytokine profile in genital tract secretions from female adolescents: impact of human immunodeficiency virus, human papillomavirus, and other sexually transmitted pathogens. J Infect Dis. 2000;181:939–945. [PubMed]

42. Silverberg MJ, Schneider MF, Silver B, Anastos KM, Burk RD, et al. Serological detection of human papillomavirus type 16 infection in human immunodeficiency virus (HIV)-positive and high-risk HIV-negative women. Clinical and vaccine immunology : CVI. 2006;13:511–519. [PMC free article] [PubMed]

43. Lorincz AT, Castle PE, Sherman ME, Scott DR, Glass AG, et al. Viral load of human papillomavirus and risk of CIN3 or cervical cancer. Lancet. 2002;360:228–229. [PubMed]

44. Lillo FB, Ferrari D, Veglia F, Origoni M, Grasso MA, et al. Human papillomavirus infection and associated cervical disease in human immunodeficiency virus-infected women: effect of highly active antiretroviral therapy. J Infect Dis. 2001;184:547–551. [PubMed]

45. Rosenblatt A, Guidi HGC. K; 2009. Human papillomavirus infection in HIV-infected individuals Human Papillomavirus: Springer-Verlag Berlin and Heidelberg GmbH & Co. pp. 163–177.

46. Hampson L, Kitchener HC, Hampson IN. Specific HIV protease inhibitors inhibit the ability of HPV16 E6 to degrade p53 and selectively kill E6-dependent cervical carcinoma cells in vitro. Antiviral Therapy. 2006;11:813–825. [PubMed]

47. Johansson M, Schon K, Ward M, Lycke N. Studies in knockout mice reveal that anti-chlamydial protection requires TH1 cells producing IFN-gamma: is this true for humans? Scandinavian Journal of Immunology. 1997;46:546–552. [PubMed]

48. Su H, Caldwell HD. CD4+ T cells play a significant role in adoptive immunity to Chlamydia trachomatis infection of the mouse genital tract. Infection and Immunity. 1995;63:3302–3308. [PMC free article] [PubMed]

49. Koskela P, Anttila T, Bjorge T, Brunsvig A, Dillner J, et al. Chlamydia trachomatis infection as a risk factor for invasive cervical cancer. International journal of cancer Journal international du cancer. 2000;85:35–39. [PubMed]

50. 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–428. [PubMed]

51. Yashin AI, Arbeev KG, Akushevich I, Kulminski A, Akushevich L, et al. Stochastic model for analysis of longitudinal data on aging and mortality. Mathematical Biosciences. 2007;208:538–551. [PMC free article] [PubMed]

52. Yashin AI, Manton KG. Effects of unobserved and partially covariate processes on system failure: A review of Models and estimation strategies. . Statistical Science. 1997;12:20–34.

53. Akushevich I, Kulminski A, Manton K. Life Tables with Covariates: Dynamic Model for Nonlinear Analysis of Longitudinal Data. Mathematical Population Studies. 2005;12:51–80.

Articles from PLoS ONE are provided here courtesy of **Public Library of Science**

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's national science library in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |