A Semiparametric Bayesian Approach to Epidemics, with Application to the Spread of the Coronavirus MERS in South Korea in 2015

We consider incomplete observations of stochastic processes governing the spread of infectious diseases through finite populations by way of contact. We propose a flexible semiparametric modeling framework with at least three advantages. First, it enables researchers to study the structure of a population contact network and its impact on the spread of infectious diseases. Second, it can accommodate short- and long-tailed degree distributions and detect potential superspreaders, who represent an important public health concern. Third, it addresses the important issue of incomplete data. Starting from first principles, we show when the incomplete-data generating process is ignorable for the purpose of Bayesian inference for the parameters of the population model. We demonstrate the semiparametric modeling framework by simulations and an application to the partially observed MERS epidemic in South Korea in 2015. We conclude with an extended discussion of open questions and directions for future research.

Transmission dynamics of an antimicrobial resistant Campylobacter jejuni lineage in New Zealand's commercial poultry network

Understanding the relative contribution of different between-farm transmission pathways is essential in guiding recommendations for mitigating disease spread. This study investigated the association between contact pathways linking poultry farms in New Zealand and the genetic relatedness of antimicrobial resistant Campylobacter jejuni Sequence Type 6964 (ST-6964), with the aim of identifying the most likely contact pathways that contributed to its rapid spread across the industry. Whole-genome sequencing was performed on 167C. jejuni ST-6964 isolates sampled from across 30 New Zealand commercial poultry enterprises. The genetic relatedness between isolates was determined using whole genome multilocus sequence typing (wgMLST). Permutational multivariate analysis of variance and distance-based linear models were used to explore the strength of the relationship between pairwise genetic associations among the C. jejuni isolates and each of several pairwise distance matrices, indicating either the geographical distance between farms or the network distance of transportation vehicles. Overall, a significant association was found between the pairwise genetic relatedness of the C. jejuni isolates and the parent company, the road distance and the network distance of transporting feed vehicles. This result suggests that the transportation of feed within the commercial poultry industry as well as other local contacts between flocks, such as the movements of personnel, may have played a significant role in the spread of C. jejuni. However, further information on the historical contact patterns between farms is needed to fully characterise the risk of these pathways and to understand how they could be targeted to reduce the spread of C. jejuni.

Using multiple data sources to explore disease transmission risk between commercial poultry, backyard poultry, and wild birds in New Zealand

The movements of backyard poultry and wild bird populations are known to pose a disease risk to the commercial poultry industry. However, it is often difficult to estimate this risk due to the lack of accurate data on the numbers, locations, and movement patterns of these populations. The main aim of this study was to evaluate the use of three different data sources when investigating disease transmission risk between poultry populations in New Zealand including (1) cross-sectional survey data looking at the movement of goods and services within the commercial poultry industry, (2) backyard poultry sales data from the online auction site TradeMe®, and (3) citizen science data from the wild bird monitoring project eBird. The cross-sectional survey data and backyard poultry sales data were transformed into network graphs showing the connectivity of commercial and backyard poultry producers across different geographical regions. The backyard poultry network was also used to parameterise a Susceptible-Infectious (SI) simulation model to explore the behaviour of potential disease outbreaks. The citizen science data was used to create an additional map showing the spatial distribution of wild bird observations across New Zealand. To explore the potential for diseases to spread between each population, maps were combined into bivariate choropleth maps showing the overlap between movements within the commercial poultry industry, backyard poultry trades and, wild bird observations. Network analysis revealed that the commercial poultry network was highly connected with geographical clustering around the urban centres of Auckland, New Plymouth and Christchurch. The backyard poultry network was also a highly active trade network and displayed similar geographic clustering to the commercial network. In the disease simulation models, the high connectivity resulted in all suburbs becoming infected in 96.4 % of the SI simulations. Analysis of the eBird data included reports of over 80 species; the majority of which were identified as coastal seabirds or wading birds that showed little overlap with either backyard or commercial poultry. Overall, our study findings highlight how the spatial patterns of trading activity within the commercial poultry industry, alongside the movement of backyard poultry and wild birds, have the potential to contribute significantly to the spread of diseases between these populations. However, it is clear that in order to fully understand this risk landscape, further data integration is needed; including the use of additional datasets that have further information on critical variables such as environmental factors.

A Mathematical Framework for Predicting Lifestyles of Viral Pathogens

Despite being similar in structure, functioning, and size, viral pathogens enjoy very different, usually well-defined ways of life. They occupy their hosts for a few days (influenza), for a few weeks (measles), or even lifelong (HCV), which manifests in acute or chronic infections. The various transmission routes (airborne, via direct physical contact, etc.), degrees of infectiousness (referring to the viral load required for transmission), antigenic variation/immune escape and virulence define further aspects of pathogenic lifestyles. To survive, pathogens must infect new hosts; the success determines their fitness. Infection happens with a certain likelihood during contact of hosts, where contact can also be mediated by vectors. Besides structural aspects of the host-contact network, three parameters appear to be key: the contact rate and the infectiousness during contact, which encode the mode of transmission, and third the immunity of susceptible hosts. On these grounds, what can be said about the reproductive success of viral pathogens? This is the biological question addressed in this paper. The answer extends earlier results of the author and makes explicit connection to another basic work on the evolution of pathogens. A mathematical framework is presented that models intra- and inter-host dynamics in a minimalistic but unified fashion covering a broad spectrum of viral pathogens, including those that cause flu-like infections, childhood diseases, and sexually transmitted infections. These pathogens turn out as local maxima of numerically simulated fitness landscapes. The models involve differential and integral equations, agent-based simulation, networks, and probability.

Transmission on empirical dynamic contact networks is influenced by data processing decisions

Dynamic contact data can be used to inform disease transmission models, providing insight into the dynamics of infectious diseases. Such data often requires extensive processing for use in models or analysis. Therefore, processing decisions can potentially influence the topology of the contact network and the simulated disease transmission dynamics on the network. In this study, we examine how four processing decisions, including temporal sampling window (TSW), spatial threshold of contact (SpTh), minimum contact duration (MCD), and temporal aggregation (daily or hourly) influence the information content of contact data (indicated by changes in entropy) as well as disease transmission model dynamics. We found that changes made to information content by processing decisions translated to significant impacts to the transmission dynamics of disease models using the contact data. In particular, we found that SpTh had the largest independent influence on information content, and that some output metrics (R₀, time to peak infection) were more sensitive to changes in information than others (epidemic extent). These findings suggest that insights gained from transmission modeling using dynamic contact data can be influenced by processing decisions alone, emphasizing the need to carefully consideration them prior to using contact-based models to conduct analyses, compare different datasets, or inform policy decisions.

The intrinsic vulnerability of networks to epidemics

•

We searched for properties making a network intrinsically vulnerable to epidemics.

•

We conducted simulations on both modelled and real-world contact networks.

•

Network properties may affect outbreak magnitude more than pathogen features.

•

We show how structural properties can be used to infer relative network vulnerability.

Contact networks are convenient models to investigate epidemics, with nodes and links representing potential hosts and infection pathways, respectively. The outcomes of outbreak simulations on networks are driven both by the underlying epidemic model, and by the networks’ structural properties, so that the same pathogen can generate different epidemic dynamics on different networks. Here we ask whether there are general properties that make a contact network intrinsically vulnerable to epidemics (that is, regardless of specific epidemiological parameters). By conducting simulations on a large set of modelled networks, we show that, when a broad range of network topologies is taken into account, the effect of specific network properties on outbreak magnitude is stronger than that of fundamental pathogen features such as transmission rate, infection duration, and immunization ability. Then, by focusing on a large set of real world networks of the same type (potential contacts between field voles, Microtus agrestis), we showed how network structure can be used to accurately assess the relative, intrinsic vulnerability of networks towards a specific pathogen, even when those have limited topological variability. These results have profound implications for how we prevent disease outbreaks; in many real world situations, the topology of host contact networks can be described and used to infer intrinsic vulnerability. Such an approach can increase preparedness and inform preventive measures against emerging diseases for which limited epidemiological information is available, enabling the identification of priority targets before an epidemic event.

Impact of degree truncation on the spread of a contagious process on networks

Understanding how person-to-person contagious processes spread through a population requires accurate information on connections between population members. However, such connectivity data, when collected via interview, is often incomplete due to partial recall, respondent fatigue or study design, e.g., fixed choice designs (FCD) truncate out-degree by limiting the number of contacts each respondent can report. Past research has shown how FCD truncation affects network properties, but its implications for predicted speed and size of spreading processes remain largely unexplored. To study the impact of degree truncation on predictions of spreading process outcomes, we generated collections of synthetic networks containing specific properties (degree distribution, degree-assortativity, clustering), and also used empirical social network data from 75 villages in Karnataka, India. We simulated FCD using various truncation thresholds and ran a susceptible-infectious-recovered (SIR) process on each network. We found that spreading processes propagated on truncated networks resulted in slower and smaller epidemics, with a sudden decrease in prediction accuracy at a level of truncation that varied by network type. Our results have implications beyond FCD to truncation due to any limited sampling from a larger network. We conclude that knowledge of network structure is important for understanding the accuracy of predictions of process spread on degree truncated networks.

Contact among healthcare workers in the hospital setting: developing the evidence base for innovative approaches to infection control

Background

Nosocomial, or healthcare-associated infections (HAI), exact a high medical and financial toll on patients, healthcare workers, caretakers, and the health system. Interpersonal contact patterns play a large role in infectious disease spread, but little is known about the relationship between health care workers’ (HCW) movements and contact patterns within a heath care facility and HAI. Quantitatively capturing these patterns will aid in understanding the dynamics of HAI and may lead to more targeted and effective control strategies in the hospital setting.

Methods

Staff at 3 urban university-based tertiary care hospitals in Canada completed a detailed questionnaire on demographics, interpersonal contacts, in-hospital movement, and infection prevention and control practices. Staff were divided into categories of administrative/support, nurses, physicians, and “Other HCWs” - a fourth distinct category, which excludes physicians and nurses. Using quantitative network modeling tools, we constructed the resulting HCW “co-location network” to illustrate contacts among different occupations and with locations in hospital settings.

Results

Among 3048 respondents (response rate 38%) an average of 3.79, 3.69 and 3.88 floors were visited by each HCW each week in the 3 hospitals, with a standard deviation of 2.63, 1.74 and 2.08, respectively. Physicians reported the highest rate of direct patient contacts (> 20 patients/day) but the lowest rate of contacts with other HCWs; nurses had the most extended (> 20 min) periods of direct patient contact. “Other HCWs” had the most direct daily contact with all other HCWs. Physicians also reported significantly more locations visited per week than nurses, other HCW, or administrators; nurses visited the fewest. Public spaces such as the cafeteria had the most staff visits per week, but the least mean hours spent per visit. Inpatient settings had significantly more HCW interactions per week than outpatient settings.

Conclusions

HCW contact patterns and spatial movement demonstrate significant heterogeneity by occupation. Control strategies that address this diversity among health care workers may be more effective than “one-strategy-fits-all” HAI prevention and control programs.

Electronic supplementary material

The online version of this article (10.1186/s12879-018-3093-x) contains supplementary material, which is available to authorized users.

Contacts between health care workers and patients in a short-stay geriatric unit during the peak of a seasonal influenza epidemic compared with a nonepidemic period

BACKGROUND

Patterns of contacts between health care workers and patients during seasonal epidemics are unknown. Our study objective was to compare the number and duration of contacts between health care workers and patients during a nonepidemic period versus a community influenza epidemic, and to identify supercontactors.

METHODS

Our observational study was conducted in a short-stay geriatric unit of a university hospital. Contacts between individuals were recorded by active radio frequency identification devices. Contact patterns were compared between 2 periods according to contact number and duration. Each care period lasted 5 days in the nonepidemic and influenza epidemic periods.

RESULTS

The study included 21 medical doctors, 43 nurses, and 56 patients. In total, 3,200 contacts (61.4%; 152,700 seconds) were recorded during the nonepidemic period, with 2,013 contacts (38.6%; 92,740 seconds) in the epidemic period (P = .007). More cumulative contacts occurred during the nonepidemic period between nurses and patients (n = 2,638 [82%] vs n = 1,599 [79%]), but not between patients (n = 56 [18%] vs n = 414 [21%]). Contact duration between nurses and patients lasted longer during the nonepidemic period (P = .04). During the epidemic period, 6 nurses (15%) considered to be supercontactors accounted for 44.3% of the total number of contacts with patients.

CONCLUSION

The pattern of contacts between individuals differed according to the presence or not of the community influenza peak that might have influenced the risk of nosocomial influenza.

Eight challenges for network epidemic models

Networks offer a fertile framework for studying the spread of infection in human and animal populations. However, owing to the inherent high-dimensionality of networks themselves, modelling transmission through networks is mathematically and computationally challenging. Even the simplest network epidemic models present unanswered questions. Attempts to improve the practical usefulness of network models by including realistic features of contact networks and of host-pathogen biology (e.g. waning immunity) have made some progress, but robust analytical results remain scarce. A more general theory is needed to understand the impact of network structure on the dynamics and control of infection. Here we identify a set of challenges that provide scope for active research in the field of network epidemic models.