Within-Host Bayesian Joint Modeling of Longitudinal and Time-to-Event Data of Leishmania Infection

The host immune system plays a significant role in managing and clearing pathogen material during an infection, but this complex process presents numerous challenges from a modeling perspective. There are many mathematical and statistical models for these kinds of processes that take into account a wide range of events that happen within the host. In this work, we present a Bayesian joint model of longitudinal and time-to-event data of Leishmania infection that considers the interplay between key drivers of the disease process: pathogen load, antibody level, and disease. The longitudinal model also considers approximate inflammatory and regulatory immune factors. In addition to measuring antibody levels produced by the immune system, we adapt data from CD4+ and CD8+ T cell proliferation, and expression of interleukin 10, interferon-gamma, and programmed cell death 1 as inflammatory or regulatory factors mediating the disease process. The model is developed using data collected from a cohort of dogs naturally exposed to Leishmania infantum. The cohort was chosen to start with healthy infected animals, and this is the majority of the data. The model also characterizes the relationship features of the longitudinal outcomes and time of death due to progressive Leishmania infection. In addition to describing the mechanisms causing disease progression and impacting the risk of death, we also present the model’s ability to predict individual trajectories of Canine Leishmaniosis (CanL) progression. The within-host model structure we present here provides a way forward to address vital research questions regarding the understanding progression of complex chronic diseases such as Visceral Leishmaniasis, a parasitic disease causing significant morbidity worldwide.