Discrepancies between observed data and predictions from mathematical modelling of the impact of screening interventions on Chlamydia trachomatis prevalence

Mathematical modelling studies of C. trachomatis transmission predict that interventions to screen and treat chlamydia infection will reduce prevalence to a greater degree than that observed in empirical population-based studies. We investigated two factors that might explain this discrepancy: partial immunity after natural infection clearance and differential screening coverage according to infection risk. We used four variants of a compartmental model for heterosexual C. trachomatis transmission, parameterized using data from England about sexual behaviour, C. trachomatis testing, diagnosis and prevalence, and Markov Chain Monte Carlo methods for statistical inference. In our baseline scenario, a model in which partial immunity follows natural infection clearance and the proportion of tests done in chlamydia-infected people decreases over time fitted the data best. The model predicts that partial immunity reduced susceptibility to reinfection by 68% (95% Bayesian credible interval 46–87%). The estimated screening rate was 4.3 (2.2–6.6) times higher for infected than for uninfected women in 2000, decreasing to 2.1 (1.4–2.9) in 2011. Despite incorporation of these factors, the model still predicted a marked decline in C. trachomatis prevalence. To reduce the gap between modelling and data, advances are needed in knowledge about factors influencing the coverage of chlamydia screening, the immunology of C. trachomatis and changes in C. trachomatis prevalence at the population level.