Unlocking pathogen genotyping information for public health by mathematical modeling

Disease modeling for public health: added value, challenges, and institutional constraints

Public health policymakers face increasingly complex questions and decisions and need to deal with an increasing quantity of data and information. For policy advisors to make use of scientific evidence and to assess available intervention options effectively and therefore indirectly for those deciding on and implementing public health policies, mathematical modeling has proven to be a useful tool. In some areas, the use of mathematical modeling for public health policy support has become standard practice at various levels of decision-making. To make use of this tool effectively within public health organizations, it is necessary to provide good infrastructure and ensure close collaboration between modelers and policymakers. Based on experience from a national public health institute, we discuss the strategic requirements for good modeling practice for public health. For modeling to be of maximal value for a public health institute, the organization and budgeting of mathematical modeling should be transparent, and a long-term strategy for how to position and develop mathematical modeling should be in place.

Prospects of elimination of HIV with test-and-treat strategy

Recently, there has been much debate about the prospects of eliminating HIV from high endemic countries by a test-and-treat strategy. This strategy entails regular HIV testing in the entire population and starting antiretroviral treatment immediately in all who are found to be HIV infected. We present the concept of the elimination threshold and investigate under what conditions of treatment uptake and dropout elimination of HIV is feasible. We used a deterministic model incorporating an accurate description of disease progression and variable infectivity. We derived explicit expressions for the basic reproduction number and the elimination threshold. Using estimates of exponential growth rates of HIV during the initial phase of epidemics, we investigated for which populations elimination is within reach. The concept of the elimination threshold allows an assessment of the prospects of elimination of HIV from information in the early phase of the epidemic. The relative elimination threshold quantifies prospects of elimination independently of the details of the transmission dynamics. Elimination of HIV by test-and-treat is only feasible for populations with very low reproduction numbers or if the reproduction number is lowered significantly as a result of additional interventions. Allowing low infectiousness during primary infection, the likelihood of elimination becomes somewhat higher. The elimination threshold is a powerful tool for assessing prospects of elimination from available data on epidemic growth rates of HIV. Empirical estimates of the epidemic growth rate from phylogenetic studies were used to assess the potential for elimination in specific populations.

How the dynamics and structure of sexual contact networks shape pathogen phylogenies

The characteristics of the host contact network over which a pathogen is transmitted affect both epidemic spread and the projected effectiveness of control strategies. Given the importance of understanding these contact networks, it is unfortunate that they are very difficult to measure directly. This challenge has led to an interest in methods to infer information about host contact networks from pathogen phylogenies, because in shaping a pathogen's opportunities for reproduction, contact networks also shape pathogen evolution. Host networks influence pathogen phylogenies both directly, through governing opportunities for evolution, and indirectly by changing the prevalence and incidence. Here, we aim to separate these two effects by comparing pathogen evolution on different host networks that share similar epidemic trajectories. This approach allows use to examine the direct effects of network structure on pathogen phylogenies, largely controlling for confounding differences arising from population dynamics. We find that networks with more heterogeneous degree distributions yield pathogen phylogenies with more variable cluster numbers, smaller mean cluster sizes, shorter mean branch lengths, and somewhat higher tree imbalance than networks with relatively homogeneous degree distributions. However, in particular for dynamic networks, we find that these direct effects are relatively modest. These findings suggest that the role of the epidemic trajectory, the dynamics of the network and the inherent variability of metrics such as cluster size must each be taken into account when trying to use pathogen phylogenies to understand characteristics about the underlying host contact network.

Complex population dynamics and the coalescent under neutrality

Estimates of the coalescent effective population size N(e) can be poorly correlated with the true population size. The relationship between N(e) and the population size is sensitive to the way in which birth and death rates vary over time. The problem of inference is exacerbated when the mechanisms underlying population dynamics are complex and depend on many parameters. In instances where nonparametric estimators of N(e) such as the skyline struggle to reproduce the correct demographic history, model-based estimators that can draw on prior information about population size and growth rates may be more efficient. A coalescent model is developed for a large class of populations such that the demographic history is described by a deterministic nonlinear dynamical system of arbitrary dimension. This class of demographic model differs from those typically used in population genetics. Birth and death rates are not fixed, and no assumptions are made regarding the fraction of the population sampled. Furthermore, the population may be structured in such a way that gene copies reproduce both within and across demes. For this large class of models, it is shown how to derive the rate of coalescence, as well as the likelihood of a gene genealogy with heterochronous sampling and labeled taxa, and how to simulate a coalescent tree conditional on a complex demographic history. This theoretical framework encapsulates many of the models used by ecologists and epidemiologists and should facilitate the integration of population genetics with the study of mathematical population dynamics.