Optimally capturing latency dynamics in models of tuberculosis transmission

Abstract

Although different structures are used in modern tuberculosis (TB) models to simulate TB latency, it remains unclear whether they are all capable of reproducing the particular activation dynamics empirically observed. We aimed to determine which of these structures replicate the dynamics of progression accurately. We reviewed 88 TB-modelling articles and classified them according to the latency structure employed. We then fitted these different models to the activation dynamics observed from 1352 infected contacts diagnosed in Victoria (Australia) and Amsterdam (Netherlands) to obtain parameter estimates. Six different model structures were identified, of which only those incorporating two latency compartments were capable of reproducing the activation dynamics empirically observed. We found important differences in parameter estimates by age. We also observed marked differences between our estimates and the parameter values used in many previous models. In particular, when two successive latency phases are considered, the first period should have a duration that is much shorter than that used in previous studies. In conclusion, structures incorporating two latency compartments and age-stratification should be employed to accurately replicate the dynamics of TB latency. We provide a catalogue of parameter values and an approach to parameter estimation from empiric data for calibration of future TB-models.

Keywords: Mathematical modelling; Parameter estimation; Risk of disease activation; Tuberculosis latency.