# The association of HPV-16 seropositivity and natural immunity to reinfection: insights from compartmental models

- Igor A Korostil
^{1}Email author, - Suzanne M Garland
^{2, 3, 4}, - Matthew G Law
^{1}and - David G Regan
^{1}

**13**:83

https://doi.org/10.1186/1471-2334-13-83

© Korostil et al.; licensee BioMed Central Ltd. 2013

**Received: **28 June 2012

**Accepted: **25 January 2013

**Published: **13 February 2013

## Abstract

### Background

Seroreactivity, processes of seroconversion and seroreversion, in the context of HPV infection has been investigated in numerous studies. However, the data resulting from these studies are usually not accounted for in mathematical transmission models of various HPV types due to gaps in our understanding of the nature of seroreactivity and its implications for HPV natural history.

### Methods

In this study we selected a number of simple but plausible compartmental transmission models of HPV-16, differing in assumptions regarding the relation between seropositivity and immunity, and attempted to calibrate them to Australian HPV seroprevalence data for females and males, as well as DNA prevalence data for females, using a Bayesian model comparison procedure. We ranked the models according to both their simplicity and ability to be fitted to the data.

### Results

Our results demonstrate that models with seroreversion where seropositivity indicates only a partial or very short-term full protection against re-infection generate age-specific HPV DNA prevalence most consistent with the observed data when compared with other models.

### Conclusions

Models supporting the notion that seropositive individuals are fully immune to reinfection demonstrated consistently inferior fits to the data than other models making no such assumption.

### Keywords

Seroreactivity Compartmental model HPV Natural immunity## Background

Genital human papillomaviruses (HPV) are viral sexually transmitted infections (STIs) with around 40 types having tropism for the anogenital region. High-risk (oncogenic) HPV types 16 and 18 are more virulent than others and associated with about 70-76% of cervical cancers [1, 2]. Of these two, type 16 is more prevalent [3, 4] and responsible for the highest proportion of cervical cancers (> 50%) [1, 5]. To prevent the spread of cervical cancer and other HPV associated diseases such as genital warts caused by the low-risk (nononcogenic) HPV types 6 and 11, many developed countries are now implementing comprehensive vaccination programs [6–8] utilising highly effective prophylactic HPV vaccines Cervarix® (bivalent vaccine; protects against HPV types 16 and 18) or Gardasil® (quadrivalent vaccine; protects against HPV types 6, 11, 16 and 18).

Because cancer generally develops long after initial infection with HPV, the actual impact of vaccination programs for cancer prevention will not be known for decades after these programs have commenced. Mathematical models have therefore been commonly employed to predict the potential population-level impact of vaccination under different vaccination scenarios and assumptions regarding vaccine properties.

Mathematical transmission models can be constructed in a number of ways but deterministic compartmental models are commonly used due to their relative simplicity and tractability [9–14]. A typical compartmental model is described by a nonlinear system of ordinary differential equations (ODEs) governing changes in the number (or proportions) of individuals in pre-specified subgroups of the modelled population over time.

A necessary element of every modelling study is model calibration. Calibration is performed by adjustment of parameter values to ensure that the model predictions, which are intrinsically uncertain, are consistent with available real-life data. The accuracy of the calibration process can be iteratively improved as more data become available. Often HPV models are calibrated to HPV incidence or prevalence data collected in a particular country or jurisdiction. In view of increasing availability of data related to seroreactivity (production of antibodies in response to infection, known as seroconversion, and their decay, or seroreversion), it is timely to investigate the present possibilities to use them for model calibration. A number of studies (for example, [15–18]) report seroprevalence by age or other characteristics, estimated times from HPV DNA detection to detection of seropositivity, and rates of seroconversion or seroreversion. With the exception of a single modelling study we are aware of [9], these data have not been considered in the development and calibration of transmission because the relationship between seropositivity and immunity is not well understood. However, we believe that it is worthwhile to investigate a few possible associations that may exist between seroreactivity and HPV transmission. The key association that we focus on in this study is that between seropositivity and natural immunity developed after resolution of an HPV infection.

In this study, we develop eight compartmental models based on types SIS (Susceptible-Infected-Susceptible), SIR (Susceptible-Infected-Recovered) and SIRS (Susceptible-Infected-Recovered-Susceptible) [19, 20], which incorporate different assumptions regarding the relation between seropositivity and immunity. Our models aim to obey the principle of parsimony, which loosely states that among competing hypotheses the simplest one should be selected. In practice, this means that if the literature presents several conflicting views on an aspect of HPV transmission, we prefer to adopt the one which is described by fewer parameters. By minimizing the number of parameters we increase their explanatory power. We rank the models in terms of their simplicity combined with ability to be fitted to Australian HPV-16 DNA prevalence and seroprevalence data.

## Methods

### Modelled population

Since we intended to calibrate our models to Australian data, it was important to ensure that the population we modelled was a reasonably accurate representation of the sexually active heterosexual Australian population. We defined the modelled population as a set of non-overlapping groups of individuals stratified by gender, age, sexual activity and infection state. In compartmental models (sometimes referred to as population-based models, in contrast with individual-based models), each of these groups (“compartments”) is assumed to be large enough to behave independently of individual stochastic effects. Throughout this paper, when we refer to an “individual” from a particular compartment, we actually mean a descriptor representing the whole population in that compartment, whose attributes are averaged attributes of that population. The age structure of the population was represented by 48 one-year age groups in the range 12 to 59 years of age. This was motivated by the following factors: 1) to model HPV we need to model only the sexually active Australian population, which excludes those younger than a certain age: in our models individuals start sexual activity at 15, but we also included 12–14 year olds, to allow for possible extension of the model should sexual behaviour data for this age group become available; 2) the sexual behaviour data we used [21] do not cover individuals older than 60, and no alternative data were available.

Sexual activity was described by four groups defined by the annual number of new sexual partners. These groups are numbered 1 to 4 in order of increasing activity and contain 60%, 27%, 11% and 2% of the modelled population, respectively [22]. Infection states were as follows: susceptible (S; an individual is susceptible if he or she is at risk of infection); infected (I; an individual who is currently infected; infected individuals are assumed to be both infectious and DNA positive); and recovered (R; an individual in this infection state has resolved his or her infection, tests HPV DNA negative, and is fully immune to re-infection). For compartmental models, the total size of the population is not important for analysis, but was set arbitrarily to 100,000. The population was assumed to be closed, i.e. immigrants, emigrants and temporary visitors were not considered (whether this plays an important role for HPV transmission in Australia has not been established but determining this is outside the scope of the present study). Mortality was also not considered for the following reasons: 1) we did not model progression to cancer and there is no mortality directly associated with HPV infection; 2) the age-specific mortality profile for the age range of the modelled population is fairly flat such that deaths from other causes can be ignored [23]. While births are not modelled explicitly, the modelled population is replenished with 12 year-olds at an annual rate equal to that at which exit from the sexually active population occurs at age 60.

### Sexual mixing

Sexual behaviour in the Australian population is described in our models by means of a mixing matrix which quantifies the rate of new partner acquisition by males and females based on their age and level of sexual activity. Our implementation of the mixing matrix is as previously employed and described in [9, 11], and is based on a formulation developed by Garnett and Anderson [24]. The parameters of the matrix are specified as in [22] based on an analysis of data from the Australian Study of Health and Relationships (ASHR) [21].

### Natural history of HPV-16

**Implementation of seroreactivity in the compared models, their DIC scores and ranking**

Model | Effect of seropositivity on transmission | Seroconversion before clearance | Seroreversion | DIC score | Rank | ||
---|---|---|---|---|---|---|---|

No reinfection while seropositive | Reduced risk of re-infection while seropositive | Risk of re-infection while seropositive is unchanged | |||||

SIS | – | ✓ | – | – | – | −101.3 | 3 |

SIS | – | ✓ | – | – | ✓ | −108.8 | 1 |

SIR | ✓ | – | – | – | – | −28.6 | 8 |

SIR | ✓ | – | – | ✓ | – | −31.2 | 7 |

SIRS | ✓ | – | – | – | ✓ | −40.4 | 6 |

SIRS | ✓ | – | – | ✓ | ✓ | −45.7 | 5 |

SIRS | – | – | ✓ | – | – | −93.4 | 4 |

SIRS | – | – | ✓ | – | ✓ | −105.6 | 2 |

_{1}does not incorporate seroreversion (the decay or loss of antibodies detectable by current assays in an individual): seropositive individuals remain seropositive for life. On the other hand, in SIS

_{2}they are assumed to be losing antibodies at a constant rate while in the state S+, so there is a steady migration of individuals from S+ to S- as they lose whatever degree of immunity they had.

**Model parameters and their prior distributions**

Parameter description | Symbol | Prior | Source |
---|---|---|---|

Per-partnership probability of transmission from female to male used to calculate the force of infection |
| U(0.10-1.00) | [25] |

Per-partnership probability of transmission from male to female used to calculate the force of infection |
| U(0.10-1.00) | [25] |

Average duration of infection for males | T | U(0.60,1.70) | |

Average duration of infection for females | T | U(0.75,1.50) | |

Average rate of loss of immunity for males; defined as 1/T | r | U(0.01,0.33) (SIRS | [31] (SIRS |

Average rate of loss of immunity for females; defined as 1/Tim,f i.e. the inverse of the of natural immunity for females | r | U(0.01,0.33) (SIRS | [31] (SIRS |

Probability of seroconversion for males | p | U(0.01,0.30) | [15] |

Probability of seroconversion for females | p | U(0.40,0.70) | [32] |

Average rate of seroreversion for males | r | U(0.01,0.10) | [33] |

Average rate of seroreversion for females | r | U(0.10,1.00) | |

Average degree of immunity for seropositive males | s | U(0.00,1.00) | Not available |

Average degree of immunity for seropositive females | s | U(0.10,1.00) | |

Average time to seroconversion for males (a proportion of T | T | U(0.50,0.95) | [15] |

Average time to conversion for females (a proportion of T | T | U(0.50,0.95) | |

Degree of assortativity by age group |
| U(0.10,0.90) | Not available |

Degree of assortativity by sexual activity group |
| U(0.10,0.90) | Not available |

_{1}). In SIR

_{2}, individuals can test seropositive and still remain infected for some time.

Model SIRS_{1} is essentially SIR_{1} with waning immunity. Note that we assume no difference in the rates of loss of immunity between seropositive (R+) and seronegative (R-) immune individuals, and losing immunity is equivalent to losing seropositive status. Similarly, SIRS_{2} is an extension of SIR_{2}. Model SIRS_{3}, just like SIR_{1} or SIRS_{1}, assumes that clearance and seroconversion are synchronous, but seropositivity is not an indication of immunity. Consequently, seropositive individuals lose immunity at the same rate as seronegative ones and then can become infected while testing seropositive. Seropositive status is life-long.

Finally, SIRS_{4} is SIRS_{3} with seroreversion. Both susceptible and immune individuals who seroconverted due to previous infection are losing antibodies at a constant rate which is different for males and females.

Note that in order to limit complexity, we chose not to model the scenario whereby an individual can serorevert while infected, since in the infected state the level of antibodies can be assumed to be high - there is, however, no evidence to convincingly support this hypothesis.

Model parameters are gender specific, which allows for possible differences in HPV-16 natural history between females and males. The ordinary differential equations describing the models included in this comparison are provided in the Additional file 1: Technical Appendix.

### Model comparison and calibration

According to the Bayesian approach we adopted, it is necessary to formulate our beliefs about each of the model parameters, before any data produced by the models have been observed, as probability distributions. These distributions are known as prior distributions or “priors”. The priors we used in this study are given in Table 2 and discussed in the Additional file 1: Technical Appendix. We applied a forward simulation procedure with adaptive Markov chain Monte Carlo (as described in [38]) to calibrate each of the models under consideration to Australian HPV-16 seroprevalence data for males and females [17] and DNA prevalence data for females aged 15–39 [39]. During the calibration procedure, quantities required to apply the deviance information criterion (DIC) [40] were calculated. We chose DIC as a means to quantify the parsimony of the models and the goodness of fit of the outcomes they produced to real data after careful consideration of a number of other statistics, such as Bayes factor and Akaike information criterion (see [41, 42] or [43] for an extensive overview of approaches to model selection) because it is simple to implement given the samples generated by a Markov chain Monte Carlo simulation. It is pertinent to note that fitting was done to all data at once and the results we obtained would not necessarily coincide with the results produced via fitting to only some of the data (for example, only HPV DNA prevalence). Another important detail is that we did not fix the two sexual mixing parameters *ε*
_{
a
} and *ε*
_{
r
} describing assortativity by age and sexual activity group (see the Additional file 1: Technical Appendix), instead treating them as ordinary parameters with the assigned uniform prior distributions U(0.1,0.9). This was motivated by the fact that these parameters are very difficult to derive for a particular population based on data from currently available sexual behaviour surveys. Consequently, sexual behaviour was not enforced to be identical in all models.

## Results

The DIC values, calculated for each tested model, are presented in Table 1. There is no rigorous specification of what magnitude of difference in DIC scores indicates a strong preference in favour of a model with the lower score, but some authors recommend using a difference of 10, as a rule of thumb [44]. By this rule of thumb, all models in which seropositivity is associated with full immunity (SIR_{1}, SIRS_{1}, SIR_{2}, SIRS_{2}) are clearly inferior to the other models. Allowing seroconversion prior to clearance of infection in SIR_{1}, and SIRS_{1} (which turn them into SIR_{2} and SIRS_{2}, respectively) somewhat improves their scores, but these are still not competitive. Calibration plots for all models can be found in the Additional file 1: Technical Appendix, along with comprehensive descriptions of the posterior distributions for model parameters.

Here we would like to briefly comment on some of the inferred parameter values for the two “best” models SIS_{2} and SIRS_{4}. Firstly, we observe that in SIS_{2} the per-partnership transmission probability from male to female (βf, posterior median 0.806 and the 95% Highest Posterior Density (HPD) interval, i.e. the shortest interval in parameter space which contains 95% of the distribution, (0.514-0.999)) is higher than that from female to male (βm , posterior median 0.59, 95% HPD interval 0.248-0.961)). This is also the case for SIRS_{4}, where the posterior mean for βf is 0.885 against 0.695 for βm. These values are consistent with the values predicted in other modelling studies: for example, β (assumed to be the same for female to male and male to female) was estimated at 0.8 (median) with the 95% posterior interval (0.6, 0.99) in [6], at 0.6 in [9] and 0.4 in [7]. There was, however, a study which reported higher rates of female-to-male relative to male-to-female transmission [45]. Posteriors for the average durations of HPV-16 infection are left-skewed for both genders, with median at 1.367 (years) for males and 1.30 for females in SIS_{2} and 1.48 (males) and 1.367 (females) in SIRS_{4}. In SIS_{2} the probability of seroconversion for males pm is low (median at 0.135), and for females (pf) it is not higher than the values reported in literature. In particular, its posterior median is at 0.494 while the 95% HPD interval is (0.4-0.654), which is in agreement with 0.5-0.6 suggested in [32]. Somewhat higher pm and pf were observed for SIRS_{4}. The inferred values for the degree of immunity for males do not let us make any meaningful conclusions regarding whether or not males are protected, because sm appears to have little influence on the model performance, which is evident from its nearly flat posterior and 95% HPD interval (0.001-0.903). In contrast, the degree of protection for females, sf, has a non-flat posterior, and its 95% HPD interval (0.100-0.810) suggests that we can at least be reasonably confident that it is certainly not complete and does not exceed 0.81, which is an important implication. Another modelling study [14], where degrees of natural immunity were introduced in a manner similar to ours, estimated them at 0.5 for both genders. Rates of loss of immunity (SIRS_{4}) were high for both males and females, 95% HPD interval for males is (0.365-0.999) and for females (0.403-1.0). These indicate very short average durations of natural immunity, namely, 1–2.74 years for males and 1–2.48 years for females. Finally, the rates of seroreversion under SIS_{2} are low but higher for females than for males (median 0.08 against 0.03). Under SIRS_{4} these are very similar (median 0.079 for females and 0.03 for males).

## Discussion

The results we obtained show that models assuming that seropositive individuals are fully and permanently protected from reinfection with HPV-16 are clearly inferior to the other models making no such strong assumptions. This conclusion is based on DIC scores. It is important to realise that DIC does not detect a ‘correct’ model in terms of HPV-16 transmission mechanism. Instead, it provides a quantitative model ranking which discourages complexity and is based on the ability of models under consideration (among which the ‘correct’ model may not even be present) to be fitted to the data. Hence, if a simpler model can be calibrated to the data at least as well as a more complex model, it will get a better DIC score. To receive a better DIC ranking, a more complex model would have to justify its complexity by producing a notably better fit than its simpler competitors. To further clarify the context in which our results should be viewed, we mention that our results can be meaningfully interpreted only if we completely rely on the available data – should these be extended or replaced, our results would inevitably change too. Another important aspect is that the DIC ranking factors in how well the models can be fitted to all data at once, for both males and females. If we, for instance, restricted ourselves to only calibrating the models to HPV seroprevalence, the resulting model ranking would likely be different.

As is evident from Table 1, our ‘best’ model is SIS_{2}, closely followed by SIRS_{4}. The difference in DIC scores between the two models is not substantial and hence does not imply that SIS_{2} is clearly preferable. We should note that the reason why SIS_{2} outscored SIS_{1} is inclusion of seroreversion. Indeed, it is the only difference between the models. Seroreversion in SIS_{2} is implemented with the help of two additional parameters (rsr,m and rsr,f), as compared with SIS_{1}, and nonetheless, it improved the fit substantially enough to overcome penalisation for extra parameters and get ahead of SIS_{1} by 8.5 points. The benefits of seroreversion in SIS_{1} are predictable since without it, SIS_{1} can not capture declining seroprevalence in older females. For the same reason, SIRS_{4} provided a significant improvement over SIRS_{3}. We see that seroreversion in SIS and SIRS models is crucial in terms of improving the fit to data, even though the rate of seroreversion is low.

Although the highest ranking models SIS_{2} and SIRS_{4} have different structures, as we mentioned in Results, the fitted durations of full natural immunity in SIRS_{4} are very short. Hence, this model is approaching a limit case when it almost becomes SIS_{2} (see Figure 1 and Figure 3).

It is our view, given what is currently known about immunity (in particular, the reported association between seropositivity and reductions in the number of incident infections in seropositive individuals [34, 35]), that the protection mechanism assumed in SIS_{2} may be a more realistic representation of naturally acquired protective immunity than a short but full immunity as in SIRS_{4}.

It is important to note that nearly all information available regarding the possible association of seropositivity with protective immunity has come from studies of females. The only study of males in this context that we are aware of [16] suggests that for males seropositivity is possible without any immunity. No substantiated inferences in regard to the existence of protective immunity in males resulted from our study: SIS_{2} was not sensitive to variations in the degree of immunity in males. To increase sensitivity, the amount of data used for model specification and calibration and/or their accuracy should be increased, which we expect to happen in future, when, for example, HPV DNA seroprevalence data for males become available.

Our models have a number of limitations. In particular, we assumed the duration of immunity to be the same for all ages, which is unlikely to be true in reality, and the probability of seroconversion to be independent of an individual’s age though there is some evidence to the contrary [34]. Additionally, compartmental models are inherently biased in certain respects. Because they assume a sexual contact is effectively instantaneous, to achieve better fit to real data, compartmental models need to compensate for a somewhat lowered level of sexual activity by maintaining higher probabilities of transmission and longer durations of infection (see [46] for detailed discussion). It is also important to remember that considerable uncertainty remains in our understanding of HPV natural history which influence our specification of priors for model parameters. Also, reliability of data obtained from sexual behavior surveys may be arguable. Finally, the results of this study rely on the data we calibrated our models to, which had their own limitations (see [17, 39] for discussion). Perhaps, the most evident limitation is that HPV-16 prevalence data only covered women aged 15–39.

## Conclusions

In conclusion, the models which provided the optimal combination of parsimony and goodness of fit to the currently available Australian data are these where seropositivity indicates only a partial (or very short full) immunity against re-infection and seroreversion is assumed to be taking place. Future studies will no doubt provide greater insight into the nature of acquired immunity and its association with seropositivity, enabling us to build more accurate models.

## Declarations

### Acknowledgements

We thank Julia Brotherton for her valuable help with the HPV-16 prevalence data extraction from the WHINURS study database.

This work was supported by an Australian Research Council Linkage Project [LP0883831] and a National Health and Medical Research Council Program Grant [568971]. The Kirby Institute receives funding from the Australian Government Department of Health and Ageing. The views expressed in this publication do not necessarily represent the position of the Australian Government. The Kirby Institute is affiliated with the Faculty of Medicine, University of New South Wales.

## Authors’ Affiliations

## References

- Munoz N, Bosch FX, de Sanjosé S, Herrero R, Castellsagué X, Shah KV, Snijders PJF, Meijer CJLM, Research on Cancer Multicenter Cervical Cancer Study Group IA: Epidemiologic classification of human papillomavirus types associated with cervical cancer. N Engl J Med. 2003, 348 (6): 518-527. 10.1056/NEJMoa021641.View ArticlePubMedGoogle Scholar
- Ni Li SF, Rebecca H-J, Snijders PJF, Clifford GM: Human papillomavirus type distribution in 30,848 invasive cervical cancers worldwide: variation by geographical region, histological type and year of publication. Int J Cancer. 2011, 128: 927-935. 10.1002/ijc.25396.View ArticlePubMedGoogle Scholar
- Jamison JH, Kaplan DW, Hamman R, Eagar R, Beach R, Douglas JM: Spectrum of genital human papillomavirus infection in a female adolescent population. Sex Transm Dis. 1995, 22 (4): 236-243. 10.1097/00007435-199507000-00007.View ArticlePubMedGoogle Scholar
- Kitchener HC, Almonte M, Wheeler P, Desai M, Gilham C, Bailey A, Sargent A, Peto J, Group ARTISTICTS: HPV testing in routine cervical screening: cross sectional data from the ARTISTIC trial. Br J Cancer. 2006, 95 (1): 56-61. 10.1038/sj.bjc.6603210.View ArticlePubMedPubMed CentralGoogle Scholar
- Clifford GM, Smith JS, Plummer M, Muñoz N, Franceschi S: Human papillomavirus types in invasive cervical cancer worldwide: a meta-analysis. Br J Cancer. 2003, 88 (1): 63-73. 10.1038/sj.bjc.6600688.View ArticlePubMedPubMed CentralGoogle Scholar
- Bogaards JA, Xiridou M, Coupe VM, Meijer CJ, Wallinga J, Berkhof J: Model-based estimation of viral transmissibility and infection-induced resistance from the age-dependent prevalence of infection for 14 high-risk types of human papillomavirus. Am J Epidemiol. 2010, 171 (7): 817-825. 10.1093/aje/kwp466.View ArticlePubMedGoogle Scholar
- Baussano I, Ronco G, Segnan N, French K, Vineis P, Garnett GP: HPV-16 infection and cervical cancer: modeling the influence of duration of infection and precancerous lesions. Epidemics. 2010, 2 (1): 21-28. 10.1016/j.epidem.2010.02.002.View ArticlePubMedGoogle Scholar
- Burchell AN, Richardson H, Mahmud SM, Trottier H, Tellier PP, Hanley J, Coutlee F, Franco EL: Modeling the sexual transmissibility of human papillomavirus infection using stochastic computer simulation and empirical data from a cohort study of young women in Montreal, Canada. Am J Epidemiol. 2006, 163 (6): 534-543. 10.1093/aje/kwj077.View ArticlePubMedGoogle Scholar
- Barnabas RV, Laukkanen P, Koskela P, Kontula O, Lehtinen M, Garnett GP: Epidemiology of HPV 16 and cervical cancer in Finland and the potential impact of vaccination: mathematical modelling analyses. PLoS Med. 2006, 3 (5): e138-10.1371/journal.pmed.0030138.View ArticlePubMedPubMed CentralGoogle Scholar
- Choi YH, Jit M, Gay N, Cox A, Garnett GP, Edmunds WJ: Transmission dynamic modelling of the impact of human papillomavirus vaccination in the United Kingdom. Vaccine. 2010, 28 (24): 4091-4102. 10.1016/j.vaccine.2009.09.125.View ArticlePubMedGoogle Scholar
- Elbasha EH, Dasbach EJ, Insinga RP: Model for assessing human papillomavirus vaccination strategies. Emerg Infect Dis. 2007, 13 (1): 28-41. 10.3201/eid1301.060438.View ArticlePubMedPubMed CentralGoogle Scholar
- Elbasha EH, Dasbach EJ, Insinga RP: A multi-type HPV transmission model. Bull Math Biol. 2008, 70 (8): 2126-2176. 10.1007/s11538-008-9338-x.View ArticlePubMedGoogle Scholar
- Hughes JP, Garnett GP, Koutsky L: The theoretical population-level impact of a prophylactic human papilloma virus vaccine. Epidemiology. 2002, 13 (6): 631-639. 10.1097/00001648-200211000-00006.View ArticlePubMedGoogle Scholar
- Kim JJ, Andres-Beck B, Goldie SJ: The value of including boys in an HPV vaccination programme: a cost-effectiveness analysis in a low-resource setting. Br J Cancer. 2007, 97 (9): 1322-1328. 10.1038/sj.bjc.6604023.View ArticlePubMedPubMed CentralGoogle Scholar
- Edelstein ZR, Carter JJ, Garg R, Winer RL, Feng Q, Galloway DA, Koutsky LA: Serum antibody response following genital alpha9 human papillomavirus infection in young men. J Infect Dis. 2011, 204 (2): 209-216. 10.1093/infdis/jir242.View ArticlePubMedPubMed CentralGoogle Scholar
- Lu B, Viscidi RP, Wu Y, Lee J-H, Nyitray AG, Villa LL, Lazcano-Ponce E, da Silva RJC, Baggio ML, Quiterio M: Prevalent serum antibody is not a marker of immune protection against acquisition of oncogenic HPV16 in men. Cancer Res. 2011, 72: 3-Google Scholar
- Newall AT, Brotherton JML, Quinn HE, McIntyre PB, Backhouse J, Gilbert L, Esser MT, Erick J, Bryan J, Formica N: Population seroprevalence of human papillomavirus types 6, 11, 16, and 18 in men, women, and children in Australia. Clin Infect Dis. 2008, 46 (11): 1647-1655. 10.1086/587895.View ArticlePubMedGoogle Scholar
- Syrjänen S, Waterboer T, Sarkola M, Michael K, Rintala M, Syrjänen K, Grenman S, Pawlita M: Dynamics of human papillomavirus serology in women followed up for 36 months after pregnancy. J Gen Virol. 2009, 90 (Pt 6): 1515-1526.View ArticlePubMedGoogle Scholar
- Keeling MJ, Rohani P: Modeling Infectious Diseases in Humans and Animals. 2008, Princeton, USA: Princeton University PressGoogle Scholar
- Vynnycky E, White RG: An Introduction to Infectious Disease Modelling. 2010, New York, USA: Oxford University PressGoogle Scholar
- Smith AMA, Rissel CE, Richters J, Grulich AE, de Visser RO: Sex in Australia: The rationale and methods of the Australian Study of Health and Relationships. Aust N Z J Public Health. 2003, 27: 106-117. 10.1111/j.1467-842X.2003.tb00797.x.View ArticlePubMedGoogle Scholar
- Regan DG, Philp DJ, Hocking JS, Law MG: Modelling the population-level impact of vaccination on the transmission of human papillomavirus type 16 in Australia. Sex Health. 2007, 4 (3): 147-163. 10.1071/SH07042.View ArticlePubMedGoogle Scholar
- Australian Bureau of Statistics; 3302.0 - Deaths, Australia. 2011, http://www.abs.gov.au,
- Garnett GP, Anderson RM: Balancing sexual partnerships in an age and activity stratified model of HIV transmission in heterosexual populations. IMA J Math Appl Med Biol. 1994, 11 (3): 161-192. 10.1093/imammb/11.3.161.View ArticlePubMedGoogle Scholar
- Burchell AN, Coutlée F, Tellier P-P, Hanley J, Franco EL: Genital transmission of human papillomavirus in recently formed heterosexual couples. J Infect Dis. 2011, 204 (11): 1723-1729. 10.1093/infdis/jir644.View ArticlePubMedPubMed CentralGoogle Scholar
- Giuliano AR, Lee J-H, Fulp W, Villa LL, Lazcano E, Papenfuss MR, Abrahamsen M, Salmeron J, Anic GM, Rollison DE: Incidence and clearance of genital human papillomavirus infection in men (HIM): a cohort study. Lancet. 2011, 377 (9769): 932-940. 10.1016/S0140-6736(10)62342-2.View ArticlePubMedPubMed CentralGoogle Scholar
- Giuliano AR, Lu B, Nielson CM, Flores R, Papenfuss MR, Lee J-H, Abrahamsen M, Harris RB: Age-specific prevalence, incidence, and duration of human papillomavirus infections in a cohort of 290 US men. J Infect Dis. 2008, 198 (6): 827-835. 10.1086/591095.View ArticlePubMedGoogle Scholar
- Giuliano AR, Harris R, Sedjo RL, Baldwin S, Roe D, Papenfuss MR, Abrahamsen M, Inserra P, Olvera S, Hatch K: Incidence, prevalence, and clearance of type-specific human papillomavirus infections: The Young Women’s Health Study. J Infect Dis. 2002, 186 (4): 462-469. 10.1086/341782.View ArticlePubMedGoogle Scholar
- Richardson H, Kelsall G, Tellier P, Voyer H, Abrahamowicz M, Ferenczy A, Coutlée F, Franco EL: The natural history of type-specific human papillomavirus infections in female university students. Cancer Epidemiol Biomarkers Prev. 2003, 12 (6): 485-490.PubMedGoogle Scholar
- Trottier H, Mahmud S, Prado JCM, Sobrinho JS, Costa MC, Rohan TE, Villa LL, Franco EL: Type-specific duration of human papillomavirus infection: implications for human papillomavirus screening and vaccination. J Infect Dis. 2008, 197 (10): 1436-1447. 10.1086/587698.View ArticlePubMedGoogle Scholar
- Shah KV, Viscidi RP, Alberg AJ, Helzlsouer KJ, Comstock GW: Antibodies to human papillomavirus 16 and subsequent in situ or invasive cancer of the cervix. Cancer Epidemiol Biomarkers Prev. 1997, 6 (4): 233-237.PubMedGoogle Scholar
- Carter JJ, Koutsky LA, Hughes JP, Lee SK, Kuypers J, Kiviat N, Galloway DA: Comparison of human papillomavirus types 16, 18, and 6 capsid antibody responses following incident infection. J Infect Dis. 2000, 181 (6): 1911-1919. 10.1086/315498.View ArticlePubMedGoogle Scholar
- Thompson DL, Douglas JM, Foster M, Hagensee ME, Diguiseppi C, Barón AE, Cameron JE, Spencer TC, Zenilman J, Malotte CK: Seroepidemiology of infection with human papillomavirus 16, in men and women attending sexually transmitted disease clinics in the United States. J Infect Dis. 2004, 190 (9): 1563-1574. 10.1086/423817.View ArticlePubMedGoogle Scholar
- Wang SS, Schiffman M, Herrero R, Carreon J, Hildesheim A, Rodriguez AC, Bratti MC, Sherman ME, Morales J, Guillen D: Determinants of human papillomavirus 16 serological conversion and persistence in a population-based cohort of 10 000 women in Costa Rica. Br J Cancer. 2004, 91 (7): 1269-1274. 10.1038/sj.bjc.6602088.View ArticlePubMedPubMed CentralGoogle Scholar
- Wentzensen N, Rodriguez AC, Viscidi R, Herrero R, Hildesheim A, Ghosh A, Morales J, Wacholder S, Guillen D, Alfaro M: A competitive serological assay shows naturally acquired immunity to human papillomavirus infections in the Guanacaste Natural History Study. J Infect Dis. 2011, 204 (1): 94-102. 10.1093/infdis/jir209.View ArticlePubMedPubMed CentralGoogle Scholar
- Ho GYF, Studentsov Y, Hall CB, Bierman R, Beardsley L, Lempa M, Burk RD: Risk factors for subsequent cervicovaginal human papillomavirus (HPV) infection and the protective role of antibodies to HPV-16 virus-like particles. J Infect Dis. 2002, 186 (6): 737-742. 10.1086/342972.View ArticlePubMedGoogle Scholar
- Safaeian M, Porras C, Schiffman M, Rodriguez AC, Wacholder S, Gonzalez P, Quint W, van Doorn L-J, Sherman ME, Xhenseval V: Epidemiological study of anti-HPV16/18 seropositivity and subsequent risk of HPV16 and −18 infections. J Natl Cancer Inst. 2010, 102 (21): 1653-1662. 10.1093/jnci/djq384.View ArticlePubMedPubMed CentralGoogle Scholar
- Korostil IA, Peters GW, Cornebise J, Regan DG: Adaptive Markov chain Monte Carlo forward projection for statistical analysis in epidemic modelling of human papillomavirus. Stat Med. 2012, 10.1002/sim.5590.Google Scholar
- Garland SM, Brotherton JML, Condon JR, McIntyre PB, Stevens MP, Smith DW, Tabrizi SN, study group WHINURS: Human papillomavirus prevalence among indigenous and non-indigenous Australian women prior to a national HPV vaccination program. BMC Med. 2011, 9: 104-10.1186/1741-7015-9-104.View ArticlePubMedPubMed CentralGoogle Scholar
- Gelman A, Carlin JB, Stern HS, Rubin DB: Bayesian data analysis. 2004, Chapman & Hall: CRC, 2Google Scholar
- Claeskens G, Hjort NL: Model Selection and Model Averaging. 2008, New York, USA: Cambridge University PressView ArticleGoogle Scholar
- Congdon P: Applied Bayesian Modelling. 2003, UK: John Wiley & Sons, LtdView ArticleGoogle Scholar
- Burnham KP, Anderson DR: Model Selection and Multi-Model Inference: A Practical Information-Theoretic Approach. 2010, New York, USA: Springer-Verlag, 2Google Scholar
- Spiegelhalter DJ, Best NG, Carlin BP, van der Linde A: Bayesian measures of model complexity and fit. Journal of Royal Statistical Society B. 2002, 64: 583-639. 10.1111/1467-9868.00353.View ArticleGoogle Scholar
- Hernandez BY, Wilkens LR, Zhu X, Thompson P, McDuffie K, Shvetsov YB, Kamemoto LE, Killeen J, Ning L, Goodman MT: Transmission of human papillomavirus in heterosexual couples. Emerg Infect Dis. 2008, 14 (6): 888-894.View ArticlePubMedPubMed CentralGoogle Scholar
- Ferguson NM, Garnett GP: More Realistic Models of Sexually Transmitted Disease Transmission Dynamics: Sexual Partnership Networks, Pair Models, and Moment Closure. Sexuall. 2000, 27: 600-609.Google Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2334/13/83/prepub

### Pre-publication history

## Copyright

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.