Combinatorial decomposition of an outbreak signature

Abstract

We use mathematically rigorous definitions of epidemiological concepts in order to derive a sequential combinatorial model of disease outbreak decomposition. We define the idea of a population specific ‘disease signature’ and use this in order to decompose and further understand outbreaks as incidents of spatial and temporal spread of disease exposure both in, and across, populations. This allows us to differentiate between different disease spread scenarios with a level of sensitivity that previous models were unable to provide. This perspective leads us to propose a new practical definition for ‘outbreak’. In addition, we are able to use this model to understand, estimate, and, in some cases, correct for, the likely instances of reporting error inherent in disease surveillance.

We demonstrate our model first with a hypothetical outbreak scenario and then in an analysis of suspected outbreaks of waterborne diseases in Massachusetts (MA) in 1995.

Introduction

Mathematical models have long been recognized as useful epidemiological tools. They provide a foundation for quantitative predictions, allow for rigorous testing of hypotheses, and necessitate clear definitions of concepts and parameters. The complicated and diverse array of infectious diseases lead to the generation of generalized models that need to be tailored by the use of carefully generated parameters to provide direct insight into the mechanisms of transmission of a particular pathogen, the rate of infection and likelihood of widespread outbreak given certain circumstances, all of which have direct applications for health care management and disease control. Traditionally, these tailored parameters have been mass action transfer terms governing the respective likelihoods of an individual transitioning from susceptible to infected to recovered over time (SIR models, cf. [1]). In cases involving complex circumstances such as multiple distinct populations or repeated, isolated exposures, these models can become complicated, the determination of appropriate mass action terms for each separate population can be difficult [2] and, in some cases, the use of mass action terms themselves can be inappropriate for the focus of the investigation [3]. Additionally, while these mass action transfer terms are mathematically meaningful, they are clinically difficult to measure, creating a disparity between mathematical elegance and usefulness.

By focusing our models on a narrower set of pathogens, those where the link between exposure and infection is clearly defined (as opposed to diseases where there can be multiple and confounding factors), we are able to use the timing of disease incidence and different etiologies specific to the different affected subpopulations to fully understand the dynamics of disease outbreaks. We here propose a method of sequential combinatorial decomposition to accomplish this narrower focus, allowing us to incorporate an understanding of the different temporal distributions governing the transition from susceptible to infected to recovered associated with each population. This method embodies a compromise between the complexity of individual behavior and the broad-brush assumption of mass action, population averages, and is based on a set of clinically measurable parameters. Our system of choice for this study will be waterborne illness due to the clearly defined direct link between exposure and infection. While we have chosen this system for study here, this method may be applied to any system so long as that link is unambiguously understood. All implications of our model are meant to be representative only of this type of system, though the theory may be generalized to others.

Unlike most SIR models that focus solely or primarily on secondary (human-to-human) transmission, we emphasize primary transmission for waterborne illnesses (i.e. an infection from an environmental contaminant/external point source). This is an important aspect because many waterborne diseases have this sort of external transmission that at least sparks an outbreak. We here present a set of rigorous definitions and operational rules that lead to a natural characterization of disease spread. This provides an outbreak signature decomposition model through heterogeneous populations. Additionally, our model will provide a natural, practical, mathematical definition of an ‘outbreak’. We will present our method by analyzing both data from a simulated scenario and actual data from the suspected cryptosporidiosis and giardiasis outbreaks in Massachusetts, USA during 1995 [4].