Host mobility drives pathogen competition in spatially structured populations

Effect of Human Behavior on the Evolution of Viral Strains During an Epidemic

It is well known in the literature that human behavior can change as a reaction to disease observed in others, and that such behavioral changes can be an important factor in the spread of an epidemic. It has been noted that human behavioral traits in disease avoidance are under selection in the presence of infectious diseases. Here, we explore a complementary trend: the pathogen itself might experience a force of selection to become less “visible,” or less “symptomatic,” in the presence of such human behavioral trends. Using a stochastic SIR agent-based model, we investigated the co-evolution of two viral strains with cross-immunity, where the resident strain is symptomatic while the mutant strain is asymptomatic. We assumed that individuals exercised self-regulated social distancing (SD) behavior if one of their neighbors was infected with a symptomatic strain. We observed that the proportion of asymptomatic carriers increased over time with a stronger effect corresponding to higher levels of self-regulated SD. Adding mandated SD made the effect more significant, while the existence of a time-delay between the onset of infection and the change of behavior reduced the advantage of the asymptomatic strain. These results were consistent under random geometric networks, scale-free networks, and a synthetic network that represented the social behavior of the residents of New Orleans.

Activity-driven network modeling and control of the spread of two concurrent epidemic strains

The emergency generated by the current COVID-19 pandemic has claimed millions of lives worldwide. There have been multiple waves across the globe that emerged as a result of new variants, due to arising from unavoidable mutations. The existing network toolbox to study epidemic spreading cannot be readily adapted to the study of multiple, coexisting strains. In this context, particularly lacking are models that could elucidate re-infection with the same strain or a different strain-phenomena that we are seeing experiencing more and more with COVID-19. Here, we establish a novel mathematical model to study the simultaneous spreading of two strains over a class of temporal networks. We build on the classical susceptible-exposed-infectious-removed model, by incorporating additional states that account for infections and re-infections with multiple strains. The temporal network is based on the activity-driven network paradigm, which has emerged as a model of choice to study dynamic processes that unfold at a time scale comparable to the network evolution. We draw analytical insight from the dynamics of the stochastic network systems through a mean-field approach, which allows for characterizing the onset of different behavioral phenotypes (non-epidemic, epidemic, and endemic). To demonstrate the practical use of the model, we examine an intermittent stay-at-home containment strategy, in which a fraction of the population is randomly required to isolate for a fixed period of time.

Diffusion dynamics of competing information on networks

Information diffusion on social networks has been described as a collective outcome of threshold behaviors in the framework of threshold models. However, since the existing models do not take into account individuals' optimization problems, it remains an open question what dynamics emerge in the diffusion process when individuals face multiple (and possibly incompatible) information sources. Here, we develop a microfounded general threshold model that enables us to analyze the collective dynamics of individual behavior in the propagation of multiple information sources. The analysis reveals that the virality of competing information sources is fundamentally indeterminate. When individuals maximize coordination with neighbors, the diffusion process is described as a saddle path, thereby leading to unpredictable symmetry breaking. When individuals' choices are irreversible, there is a continuum of stable equilibria where a certain degree of social polarization takes place by chance.

Reconstructing contact network structure and cross-immunity patterns from multiple infection histories

Interactions within a population shape the spread of infectious diseases but contact patterns between individuals are difficult to access. We hypothesised that key properties of these patterns can be inferred from multiple infection data in longitudinal follow-ups. We developed a simulator for epidemics with multiple infections on networks and analysed the resulting individual infection time series by introducing similarity metrics between hosts based on their multiple infection histories. We find that, depending on infection multiplicity and network sampling, multiple infection summary statistics can recover network properties such as degree distribution. Furthermore, we show that by mining simulation outputs for multiple infection patterns, one can detect immunological interference between pathogens (i.e. the fact that past infections in a host condition future probability of infection). The combination of individual-based simulations and analysis of multiple infection histories opens promising perspectives to infer and validate transmission networks and immunological interference for infectious diseases from longitudinal cohort data.

Connectivity, reproduction number, and mobility interact to determine communities' epidemiological superspreader potential in a metapopulation network

Disease epidemic outbreaks on human metapopulation networks are often driven by a small number of superspreader nodes, which are primarily responsible for spreading the disease throughout the network. Superspreader nodes typically are characterized either by their locations within the network, by their degree of connectivity and centrality, or by their habitat suitability for the disease, described by their reproduction number (R). Here we introduce a model that considers simultaneously the effects of network properties and R on superspreaders, as opposed to previous research which considered each factor separately. This type of model is applicable to diseases for which habitat suitability varies by climate or land cover, and for direct transmitted diseases for which population density and mitigation practices influences R. We present analytical models that quantify the superspreader capacity of a population node by two measures: probability-dependent superspreader capacity, the expected number of neighboring nodes to which the node in consideration will randomly spread the disease per epidemic generation, and time-dependent superspreader capacity, the rate at which the node spreads the disease to each of its neighbors. We validate our analytical models with a Monte Carlo analysis of repeated stochastic Susceptible-Infected-Recovered (SIR) simulations on randomly generated human population networks, and we use a random forest statistical model to relate superspreader risk to connectivity, R, centrality, clustering, and diffusion. We demonstrate that either degree of connectivity or R above a certain threshold are sufficient conditions for a node to have a moderate superspreader risk factor, but both are necessary for a node to have a high-risk factor. The statistical model presented in this article can be used to predict the location of superspreader events in future epidemics, and to predict the effectiveness of mitigation strategies that seek to reduce the value of R, alter host movements, or both.

Mixing protocols in the public goods game

If interaction partners in social dilemma games are not selected randomly from the population but are instead determined by a network of contacts, it has far reaching consequences for the evolutionary dynamics. Selecting partners randomly leads to a well-mixed population, where pattern formation is essentially impossible. This rules out important mechanisms that can facilitate cooperation, most notably network reciprocity. In contrast, if interactions are determined by a lattice or a network, then the population is said to be structured, where cooperators can form compact clusters that protect them from invading defectors. Between these two extremes, however, there is ample middle ground that can be brought about by the consideration of temporal networks, mobility, or other coevolutionary processes. The question that we here seek to answer is, when does mixing on a lattice actually lead to well-mixed conditions? To that effect, we use the public goods game on a square lattice, and we consider nearest-neighbor and random mixing with different frequencies, as well as a mix of both mixing protocols. Not surprisingly, we find that nearest-neighbor mixing requires a higher frequency than random mixing to arrive at the well-mixed limit. The differences between the two mixing protocols are most expressed at intermediate mixing frequencies, whilst at very low and very high mixing frequencies the two almost converge. We also find a near universal exponential growth of the average size of cooperator clusters as their fraction increases from zero to one, regardless of whether this increase is due to increasing the multiplication factor of the public goods, decreasing the frequency of mixing, or gradually shifting the mixing from random to nearest neighbors.

Modelling the impact of testing, contact tracing and household quarantine on second waves of COVID-19

While severe social-distancing measures have proven effective in slowing the coronavirus disease 2019 (COVID-19) pandemic, second-wave scenarios are likely to emerge as restrictions are lifted. Here we integrate anonymized, geolocalized mobility data with census and demographic data to build a detailed agent-based model of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) transmission in the Boston metropolitan area. We find that a period of strict social distancing followed by a robust level of testing, contact-tracing and household quarantine could keep the disease within the capacity of the healthcare system while enabling the reopening of economic activities. Our results show that a response system based on enhanced testing and contact tracing can have a major role in relaxing social-distancing interventions in the absence of herd immunity against SARS-CoV-2.

Modeling the impact of social distancing, testing, contact tracing and household quarantine on second-wave scenarios of the COVID-19 epidemic

The new coronavirus disease 2019 (COVID-19) has required the implementation of severe mobility restrictions and social distancing measures worldwide. While these measures have been proven effective in abating the epidemic in several countries, it is important to estimate the effectiveness of testing and tracing strategies to avoid a potential second wave of the COVID-19 epidemic. We integrate highly detailed (anonymized, privacy-enhanced) mobility data from mobile devices, with census and demographic data to build a detailed agent-based model to describe the transmission dynamics of SARS-CoV-2 in the Boston metropolitan area. We find that enforcing strict social distancing followed by a policy based on a robust level of testing, contact-tracing and household quarantine, could keep the disease at a level that does not exceed the capacity of the health care system. Assuming the identification of 50% of the symptomatic infections, and the tracing of 40% of their contacts and households, which corresponds to about 9% of individuals quarantined, the ensuing reduction in transmission allows the reopening of economic activities while attaining a manageable impact on the health care system. Our results show that a response system based on enhanced testing and contact tracing can play a major role in relaxing social distancing interventions in the absence of herd immunity against SARS-CoV-2.

Interplay between competitive and cooperative interactions in a three-player pathogen system

In ecological systems, heterogeneous interactions between pathogens take place simultaneously. This occurs, for instance, when two pathogens cooperate, while at the same time, multiple strains of these pathogens co-circulate and compete. Notable examples include the cooperation of human immunodeficiency virus with antibiotic-resistant and susceptible strains of tuberculosis or some respiratory infections with Streptococcus pneumoniae strains. Models focusing on competition or cooperation separately fail to describe how these concurrent interactions shape the epidemiology of such diseases. We studied this problem considering two cooperating pathogens, where one pathogen is further structured in two strains. The spreading follows a susceptible-infected-susceptible process and the strains differ in transmissibility and extent of cooperation with the other pathogen. We combined a mean-field stability analysis with stochastic simulations on networks considering both well-mixed and structured populations. We observed the emergence of a complex phase diagram, where the conditions for the less transmissible, but more cooperative strain to dominate are non-trivial, e.g. non-monotonic boundaries and bistability. Coupled with community structure, the presence of the cooperative pathogen enables the coexistence between strains by breaking the spatial symmetry and dynamically creating different ecological niches. These results shed light on ecological mechanisms that may impact the epidemiology of diseases of public health concern.

Coevolution spreading in complex networks

The propagations of diseases, behaviors and information in real systems are rarely independent of each other, but they are coevolving with strong interactions. To uncover the dynamical mechanisms, the evolving spatiotemporal patterns and critical phenomena of networked coevolution spreading are extremely important, which provide theoretical foundations for us to control epidemic spreading, predict collective behaviors in social systems, and so on. The coevolution spreading dynamics in complex networks has thus attracted much attention in many disciplines. In this review, we introduce recent progress in the study of coevolution spreading dynamics, emphasizing the contributions from the perspectives of statistical mechanics and network science. The theoretical methods, critical phenomena, phase transitions, interacting mechanisms, and effects of network topology for four representative types of coevolution spreading mechanisms, including the coevolution of biological contagions, social contagions, epidemic-awareness, and epidemic-resources, are presented in detail, and the challenges in this field as well as open issues for future studies are also discussed.

Host contact dynamics shapes richness and dominance of pathogen strains

The interaction among multiple microbial strains affects the spread of infectious diseases and the efficacy of interventions. Genomic tools have made it increasingly easy to observe pathogenic strains diversity, but the best interpretation of such diversity has remained difficult because of relationships with host and environmental factors. Here, we focus on host-to-host contact behavior and study how it changes populations of pathogens in a minimal model of multi-strain interaction. We simulated a population of identical strains competing by mutual exclusion and spreading on a dynamical network of hosts according to a stochastic susceptible-infectious-susceptible model. We computed ecological indicators of diversity and dominance in strain populations for a collection of networks illustrating various properties found in real-world examples. Heterogeneities in the number of contacts among hosts were found to reduce diversity and increase dominance by making the repartition of strains among infected hosts more uneven, while strong community structure among hosts increased strain diversity. We found that the introduction of strains associated with hosts entering and leaving the system led to the highest pathogenic richness at intermediate turnover levels. These results were finally illustrated using the spread of Staphylococcus aureus in a long-term health-care facility where close proximity interactions and strain carriage were collected simultaneously. We found that network structural and temporal properties could account for a large part of the variability observed in strain diversity. These results show how stochasticity and network structure affect the population ecology of pathogens and warn against interpreting observations as unambiguous evidence of epidemiological differences between strains.

Endemicity and prevalence of multipartite viruses under heterogeneous between-host transmission

Multipartite viruses replicate through a puzzling evolutionary strategy. Their genome is segmented into two or more parts, and encapsidated in separate particles that appear to propagate independently. Completing the replication cycle, however, requires the full genome, so that a systemic infection of a host requires the concurrent presence of several particles. This represents an apparent evolutionary drawback of multipartitism, while its advantages remain unclear. A transition from monopartite to multipartite viral forms has been described in vitro under conditions of high multiplicity of infection, suggesting that cooperation between defective mutants is a plausible evolutionary pathway towards multipartitism. However, it is unknown how the putative advantages that multipartitism might enjoy at the microscopic level affect its epidemiology, or if an explicit advantange is needed to explain its ecological persistence. In order to disentangle which mechanisms might contribute to the rise and fixation of multipartitism, we here investigate the interaction between viral spreading dynamics and host population structure. We set up a compartmental model of the spread of a virus in its different forms and explore its epidemiology using both analytical and numerical techniques. We uncover that the impact of host contact structure on spreading dynamics entails a rich phenomenology of ecological relationships that includes cooperation, competition, and commensality. Furthermore, we find out that multipartitism might rise to fixation even in the absence of explicit microscopic advantages. Multipartitism allows the virus to colonize environments that could not be invaded by the monopartite form, while homogeneous contacts between hosts facilitate its spread. We conjecture that these features might have led to an increase in the diversity and prevalence of multipartite viral forms concomitantly with the expansion of agricultural practices.

Persistent spatial patterns of interacting contagions

The spread of infectious diseases, rumors, fashions, and innovations are complex contagion processes, embedded in network and spatial contexts. While the studies in the former context are intensively expanded, the latter remains largely unexplored. In this paper, we investigate the pattern formation of an interacting contagion, where two infections, A and B, interact with each other and diffuse simultaneously in space. The contagion process for each follows the classical susceptible-infected-susceptible kinetics, and their interaction introduces a potential change in the secondary infection propensity compared to the baseline reproduction number R_{0}. We show that the nontrivial spatial infection patterns arise when the susceptible individuals move faster than the infected and the interaction between the two infections is neither too competitive nor too cooperative. Interestingly, the system exhibits pattern hysteresis phenomena, i.e., quite different parameter regions for patterns exist in the direction of increasing or decreasing R_{0}. Decreasing R_{0} reveals remarkable enhancement in contagion prevalence, meaning that the eradication becomes difficult compared to the single-infection or coinfection without space. Linearization analysis supports our observations, and we have identified the required elements and dynamical mechanism, which suggests that these patterns are essentially Turing patterns. Our work thus reveals new complexities in interacting contagions and paves the way for further investigation because of its relevance to both biological and social contexts.

Limited role of spatial self-structuring in emergent trade-offs during pathogen evolution

Pathogen transmission and virulence are main evolutionary variables broadly assumed to be linked through trade-offs. In well-mixed populations, these trade-offs are often ascribed to physiological restrictions, while populations with spatial self-structuring might evolve emergent trade-offs. Here, we reexamine a spatially-explicit, SIR model of the latter kind proposed by Ballegooijen and Boerlijst with the aim of characterising the mechanisms causing the emergence of the trade-off and its structural robustness. Using invadability criteria, we establish the conditions under which an evolutionary feedback between transmission and virulence mediated by pattern formation can poise the system to a critical boundary separating a disordered state (without emergent trade-off) from a self-structured phase (where the trade-off emerges), and analytically calculate the functional shape of the boundary in a certain approximation. Beyond evolutionary parameters, the success of an invasion depends on the size and spatial structure of the invading and invaded populations. Spatial self-structuring is often destroyed when hosts are mobile, changing the evolutionary dynamics to those of a well-mixed population. In a metapopulation scenario, the systematic extinction of the pathogen in the disordered phase may counteract the disruptive effect of host mobility, favour pattern formation and therefore recover the emergent trade-off.

Diversity of hysteresis in a fully cooperative coinfection model

We propose a fully cooperative coinfection model in which singly infected individuals are more likely to acquire a second disease than susceptible ones and doubly infected individuals are also assumed to be more contagious than singly infected ones. The dynamics of such a fully cooperative coinfection model is investigated through the well-mixed approach. In particular, discontinuous outbreak transitions from the disease free state or the low prevalence state to the high prevalence state can be separately observed as a disease transmission rate crosses a threshold α_o from the below when the epidemic is still in the early stages. Moreover, discontinuous eradications from the high prevalence state to the low prevalence or disease free state are also separately seen as the transmission rate reaches a threshold α_e(<α_o) from the above when the outbreak occurs. Such phenomena constitute three types of hysteresis, where only one type has been identified before. Complete characterization of these three types of hysteresis in terms of parameters measuring the uniformity of the model is both analytically and numerically provided.

Moment closure of infectious diseases model on heterogeneous metapopulation network

The global transmission of infectious diseases poses huge threats to human. Traditional heterogeneous mean-field models on metapopulation networks ignore the heterogeneity of individuals who are in different disease states in subpopulations with the same degree, resulting in inaccuracy in predicting the spread of disease. In this paper, we take heterogeneity of susceptible and infectious individuals in subpopulations with the same degree into account, and propose a deterministic unclosed general model according to Markov process on metapopulation networks to curve the global transmission of diseases precisely. Then we make the general model closed by putting forward two common assumptions: a two-dimensional constant distribution and a two-dimensional log-normal distribution, where the former is equivalent to the heterogeneous mean-field model, and the latter is a system of weighted ordinary differential equations. Further we make a stability analysis for two closed models and illustrate the results by numerical simulations. Next, we conduct a series of numerical simulations and stochastic simulations. Results indicate that our general model extends and optimizes the mean-field model. Finally, we investigate the impacts of total mobility rate on disease transmission and find that timely and comprehensive travel restriction in the early stage is an effective prevention and control of infectious diseases.

Spatial evolutionary epidemiology of spreading epidemics

Most spatial models of host-parasite interactions either neglect the possibility of pathogen evolution or consider that this process is slow enough for epidemiological dynamics to reach an equilibrium on a fast timescale. Here, we propose a novel approach to jointly model the epidemiological and evolutionary dynamics of spatially structured host and pathogen populations. Starting from a multi-strain epidemiological model, we use a combination of spatial moment equations and quantitative genetics to analyse the dynamics of mean transmission and virulence in the population. A key insight of our approach is that, even in the absence of long-term evolutionary consequences, spatial structure can affect the short-term evolution of pathogens because of the build-up of spatial differentiation in mean virulence. We show that spatial differentiation is driven by a balance between epidemiological and genetic effects, and this quantity is related to the effect of kin competition discussed in previous studies of parasite evolution in spatially structured host populations. Our analysis can be used to understand and predict the transient evolutionary dynamics of pathogens and the emergence of spatial patterns of phenotypic variation.

Dynamic, spatial models of parasite transmission in wildlife: Their structure, applications and remaining challenges

Individual differences in contact rate can arise from host, group and landscape heterogeneity and can result in different patterns of spatial spread for diseases in wildlife populations with concomitant implications for disease control in wildlife of conservation concern, livestock and humans. While dynamic disease models can provide a better understanding of the drivers of spatial spread, the effects of landscape heterogeneity have only been modelled in a few well-studied wildlife systems such as rabies and bovine tuberculosis. Such spatial models tend to be either purely theoretical with intrinsic limiting assumptions or individual-based models that are often highly species- and system-specific, limiting the breadth of their utility. Our goal was to review studies that have utilized dynamic, spatial models to answer questions about pathogen transmission in wildlife and identify key gaps in the literature. We begin by providing an overview of the main types of dynamic, spatial models (e.g., metapopulation, network, lattice, cellular automata, individual-based and continuous-space) and their relation to each other. We investigate different types of ecological questions that these models have been used to explore: pathogen invasion dynamics and range expansion, spatial heterogeneity and pathogen persistence, the implications of management and intervention strategies and the role of evolution in host-pathogen dynamics. We reviewed 168 studies that consider pathogen transmission in free-ranging wildlife and classify them by the model type employed, the focal host-pathogen system, and their overall research themes and motivation. We observed a significant focus on mammalian hosts, a few well-studied or purely theoretical pathogen systems, and a lack of studies occurring at the wildlife-public health or wildlife-livestock interfaces. Finally, we discuss challenges and future directions in the context of unprecedented human-mediated environmental change. Spatial models may provide new insights into understanding, for example, how global warming and habitat disturbance contribute to disease maintenance and emergence. Moving forward, better integration of dynamic, spatial disease models with approaches from movement ecology, landscape genetics/genomics and ecoimmunology may provide new avenues for investigation and aid in the control of zoonotic and emerging infectious diseases.

Human mobility networks and persistence of rapidly mutating pathogens

Rapidly mutating pathogens may be able to persist in the population and reach an endemic equilibrium by escaping hosts' acquired immunity. For such diseases, multiple biological, environmental and population-level mechanisms determine the dynamics of the outbreak, including pathogen's epidemiological traits (e.g. transmissibility, infectious period and duration of immunity), seasonality, interaction with other circulating strains and hosts' mixing and spatial fragmentation. Here, we study a susceptible-infected-recovered-susceptible model on a metapopulation where individuals are distributed in sub-populations connected via a network of mobility flows. Through extensive numerical simulations, we explore the phase space of pathogen's persistence and map the dynamical regimes of the pathogen following emergence. Our results show that spatial fragmentation and mobility play a key role in the persistence of the disease whose maximum is reached at intermediate mobility values. We describe the occurrence of different phenomena including local extinction and emergence of epidemic waves, and assess the conditions for large-scale spreading. Findings are highlighted in reference to previous studies and to real scenarios. Our work uncovers the crucial role of hosts' mobility on the ecological dynamics of rapidly mutating pathogens, opening the path for further studies on disease ecology in the presence of a complex and heterogeneous environment.

Disease Localization in Multilayer Networks

We present a continuous formulation of epidemic spreading on multilayer networks using a tensorial representation, extending the models of monoplex networks to this context. We derive analytical expressions for the epidemic threshold of the susceptible-infected-susceptible (SIS) and susceptible-infected-recovered dynamics, as well as upper and lower bounds for the disease prevalence in the steady state for the SIS scenario. Using the quasistationary state method, we numerically show the existence of disease localization and the emergence of two or more susceptibility peaks, which are characterized analytically and numerically through the inverse participation ratio. At variance with what is observed in single-layer networks, we show that disease localization takes place on the layers and not on the nodes of a given layer. Furthermore, when mapping the critical dynamics to an eigenvalue problem, we observe a characteristic transition in the eigenvalue spectra of the supra-contact tensor as a function of the ratio of two spreading rates: If the rate at which the disease spreads within a layer is comparable to the spreading rate across layers, the individual spectra of each layer merge with the coupling between layers. Finally, we report on an interesting phenomenon, the barrier effect; i.e., for a three-layer configuration, when the layer with the lowest eigenvalue is located at the center of the line, it can effectively act as a barrier to the disease. The formalism introduced here provides a unifying mathematical approach to disease contagion in multiplex systems, opening new possibilities for the study of spreading processes.

Networks are all around. They describe the flow of information, the movement of people and goods via multiple modes of transportation, and the spread of disease across interconnected populations. Traditionally, networks have been studied as if they were a single layer, which flattens out hierarchies such as social circles. Multilayer networks, which consider each of those circles as a layer, are more accurate descriptions of real-world networks and their use can have deep implications for understanding the dynamics of the system. Using the spread of disease as a model, we have developed a mathematical framework that accounts for the multilayer structure, and we have identified several behaviors that emerge from this analysis.

The framework relies on tensors, mathematical objects that allow us to represent multidimensional data in a compact way. Through mathematical analysis and numerical simulations, we find a number of interesting features such as the existence of multiple epidemic thresholds and transmission rates beyond which the number of individuals that catch a disease is non-negligible. We also show the existence of disease localization, a scenario in which the disease cannot escape a layer and jump to another.

Our work provides a unifying mathematical approach to studying disease transmission among multilayered populations. There are still many aspects to investigate such as how to use these results to help contain an epidemic as well as how the picture changes in more complex scenarios. Disease-like models can also be used to explore other networks such as the propagation of information.

Competing opinions and stubborness: Connecting models to data

We introduce a general contagionlike model for competing opinions that includes dynamic resistance to alternative opinions. We show that this model can describe candidate vote distributions, spatial vote correlations, and a slow approach to opinion consensus with sensible parameter values. These empirical properties of large group dynamics, previously understood using distinct models, may be different aspects of human behavior that can be captured by a more unified model, such as the one introduced in this paper.

The impact of aggregating serogroups in dynamic models of Neisseria meningitidis transmission

Background

Neisseria meningitidis (Nm) is a pathogen of multiple serogroups that is highly prevalent in many populations. Serogroups associated with invasive meningococcal disease (IMD) in Canada, for example, include A, B, C, W-135, X and Y. IMD is a rare but serious outcome of Nm infection, and can be prevented with vaccines that target certain serogroups. This has stimulated the development of dynamic models to evaluate vaccine impact. However, these models typically aggregate the various Nm serogroups into a small number of combined groups, instead of modelling each serogroup individually. The impact of aggregation on dynamic Nm model predictions is poorly understood. Our objective was to explore the impact of aggregation on dynamic model predictions.

Methods

We developed two age-structured agent-based models--a 2-strain model and a 4-strain model--to simulate vaccination programs in the Canadian setting. The 2-strain model was used to explore two different groupings: C, versus all other serogroups combined; and B, versus all other serogroups combined. The 4-strain model used the four groupings: C, B, Neisseria lactamica, versus all other serogroups combined. We compared the predicted impact of monovalent C vaccine, quadrivalent ACWY vaccine (MCV-4), and monovalent B vaccine (4CMenB) on the prevalence of serogroup carriage under these different models.

Results

The 2-strain and 4-strain models predicted similar overall impacts of vaccines on carriage prevalence, especially with respect to the vaccine-targeted serogroups. However, there were some significant quantitative and qualitative differences. Declines in vaccine-targeted serogroups were more rapid in the 2-strain model than the 4-strain model, for both the C and the 4CMenB vaccines. Sustained oscillations, and evidence for multiple attractors (i.e., different types of dynamics for the same model parameters but different initial conditions), occurred in the 4-strain model but not the 2-strain model. Strain replacement was also more pronounced in the 4-strain model, on account of the 4-strain model spreading prevalence more thinly across groups and thus enhancing competitive interactions.

Conclusions

Simplifying assumptions like aggregation of serogroups can have significant impacts on dynamic model predictions. Modellers should carefully weigh the advantages and disadvantages of aggregation when formulating models for multi-strain pathogens.

Electronic supplementary material

The online version of this article (doi:10.1186/s12879-015-1015-8) contains supplementary material, which is available to authorized users.

How demography-driven evolving networks impact epidemic transmission between communities

In this paper, we develop a complex network susceptible-infected-susceptible (SIS) model to investigate the impact of demographic factors on disease spreads. We carefully capture the transmission by short-time travelers, by assuming the susceptibles randomly travel to another community, stay for a daily time scale, and return. We calculate the basic reproductive number R0 and analyze the relevant stability of the equilibria (disease-free equilibrium and endemic equilibrium) of the model by applying limiting system theory and comparison principle. The results reveal that the disease-free equilibrium is globally asymptotically stable given R0<1, whereas the condition R0>1 leads to a globally asymptotically stable endemic equilibrium. Our numerical simulations show that demographic factors, such as birth, immigration, and short-time travels, play important roles in epidemic propagation from one community to another. Moreover, we quantitatively demonstrate how the distribution of individual's network degree would affect the result of disease transmission.

Global epidemic invasion thresholds in directed cattle subpopulation networks having source, sink, and transit nodes

Through the characterization of a metapopulation cattle disease model on a directed network having source, transit, and sink nodes, we derive two global epidemic invasion thresholds. The first threshold defines the conditions necessary for an epidemic to successfully spread at the global scale. The second threshold defines the criteria that permit an epidemic to move out of the giant strongly connected component and to invade the populations of the sink nodes. As each sink node represents a final waypoint for cattle before slaughter, the existence of an epidemic among the sink nodes is a serious threat to food security. We find that the relationship between these two thresholds depends on the relative proportions of transit and sink nodes in the system and the distributions of the in-degrees of both node types. These analytic results are verified through numerical realizations of the metapopulation cattle model.

Characterising two-pathogen competition in spatially structured environments

Different pathogens spreading in the same host population often generate complex co-circulation dynamics because of the many possible interactions between the pathogens and the host immune system, the host life cycle, and the space structure of the population. Here we focus on the competition between two acute infections and we address the role of host mobility and cross-immunity in shaping possible dominance/co-dominance regimes. Host mobility is modelled as a network of traveling flows connecting nodes of a metapopulation, and the two-pathogen dynamics is simulated with a stochastic mechanistic approach. Results depict a complex scenario where, according to the relation among the epidemiological parameters of the two pathogens, mobility can either be non-influential for the competition dynamics or play a critical role in selecting the dominant pathogen. The characterisation of the parameter space can be explained in terms of the trade-off between pathogen's spreading velocity and its ability to diffuse in a sparse environment. Variations in the cross-immunity level induce a transition between presence and absence of competition. The present study disentangles the role of the relevant biological and ecological factors in the competition dynamics, and provides relevant insights into the spatial ecology of infectious diseases.

Epidemic outbreaks in two-scale community networks

We consider a model for the diffusion of epidemics in a population that is partitioned into local communities. In particular, assuming a mean-field approximation, we analyze a continuous-time susceptible-infected-susceptible (SIS) model that has appeared recently in the literature. The probability by which an individual infects individuals in its own community is different from the probability of infecting individuals in other communities. The aim of the model, compared to the standard, nonclustered one, is to provide a compact description for the presence of communities of local infection where the epidemic process is faster compared to the rate at which it spreads across communities. Ultimately, it provides a tool to express the probability of epidemic outbreaks in the form of a metastable infection probability. In the proposed model, the spatial structure of the network is encoded by the adjacency matrix of clusters, i.e., the connections between local communities, and by the vector of the sizes of local communities. Thus, the existence of a nontrivial metastable occupancy probability is determined by an epidemic threshold which depends on the clusters' size and on the intercommunity network structure.

Competitive epidemic spreading over arbitrary multilayer networks

This study extends the Susceptible-Infected-Susceptible (SIS) epidemic model for single-virus propagation over an arbitrary graph to an Susceptible-Infected by virus 1-Susceptible-Infected by virus 2-Susceptible (SI_{1}SI_{2}S) epidemic model of two exclusive, competitive viruses over a two-layer network with generic structure, where network layers represent the distinct transmission routes of the viruses. We find analytical expressions determining extinction, coexistence, and absolute dominance of the viruses after we introduce the concepts of survival threshold and absolute-dominance threshold. The main outcome of our analysis is the discovery and proof of a region for long-term coexistence of competitive viruses in nontrivial multilayer networks. We show coexistence is impossible if network layers are identical yet possible if network layers are distinct. Not only do we rigorously prove a region of coexistence, but we can quantitate it via interrelation of central nodes across the network layers. Little to no overlapping of the layers' central nodes is the key determinant of coexistence. For example, we show both analytically and numerically that positive correlation of network layers makes it difficult for a virus to survive, while in a network with negatively correlated layers, survival is easier, but total removal of the other virus is more difficult.

Domination-time dynamics in susceptible-infected-susceptible virus competition on networks

When two viruses compete for healthy nodes in a simple network and both spreading rates are above the epidemic threshold, only one virus will survive. However, if we prevent the viruses from dying out, rich dynamics emerge. When both viruses are identical, one virus always dominates the other, but the dominating and dominated virus alternate. We show in the complete graph that the domination time depends on the total number of infected nodes at the beginning of the domination period and, moreover, that the distribution of the domination time decays exponentially yet slowly. When the viruses differ moderately in strength and/or speed the weaker and/or slower virus can still dominate the other but for a short time. Interestingly, depending on the number of infected nodes at the start of a domination period, being quicker can be a disadvantage.

Metapopulation epidemic models with heterogeneous mixing and travel behaviour

Background

Determining the pandemic potential of an emerging infectious disease and how it depends on the various epidemic and population aspects is critical for the preparation of an adequate response aimed at its control. The complex interplay between population movements in space and non-homogeneous mixing patterns have so far hindered the fundamental understanding of the conditions for spatial invasion through a general theoretical framework. To address this issue, we present an analytical modelling approach taking into account such interplay under general conditions of mobility and interactions, in the simplifying assumption of two population classes.

Methods

We describe a spatially structured population with non-homogeneous mixing and travel behaviour through a multi-host stochastic epidemic metapopulation model. Different population partitions, mixing patterns and mobility structures are considered, along with a specific application for the study of the role of age partition in the early spread of the 2009 H1N1 pandemic influenza.

Results

We provide a complete mathematical formulation of the model and derive a semi-analytical expression of the threshold condition for global invasion of an emerging infectious disease in the metapopulation system. A rich solution space is found that depends on the social partition of the population, the pattern of contacts across groups and their relative social activity, the travel attitude of each class, and the topological and traffic features of the mobility network. Reducing the activity of the less social group and reducing the cross-group mixing are predicted to be the most efficient strategies for controlling the pandemic potential in the case the less active group constitutes the majority of travellers. If instead traveling is dominated by the more social class, our model predicts the existence of an optimal across-groups mixing that maximises the pandemic potential of the disease, whereas the impact of variations in the activity of each group is less important.

Conclusions

The proposed modelling approach introduces a theoretical framework for the study of infectious diseases spread in a population with two layers of heterogeneity relevant for the local transmission and the spatial propagation of the disease. It can be used for pandemic preparedness studies to identify adequate interventions and quantitatively estimate the corresponding required effort, as well as in an emerging epidemic situation to assess the pandemic potential of the pathogen from population and early outbreak data.