Spreading Processes in Multiplex Metapopulations Containing Different Mobility Networks

Dynamics of a susceptible-infected-recovered model on complex networks with interregional migration

We present a susceptible-infected-recovered model based on a dynamic flow network that describes the epidemic process on complex metapopulation networks. This model views population regions as interconnected nodes and describes the evolution of each region using a system of differential equations. The next-generation matrix method is used to derive the global basic reproduction number for three cases: a general network with homogeneous infection rates in all regions, a fully connected network, and a star network with heterogeneous infection and recovery rates. For the homogeneous case, we show that this global basic reproduction number is independent of the migration rates between regions. However, the rate of convergence of each region to an equilibrium state exhibits a much larger variance in random (Erdős-Rényi) networks compared to small-scale (Barabási-Albert) networks. For the general heterogeneous case, we report interesting results, namely that the global basic reproduction number decays exponentially with respect to the smallest nonzero Laplacian eigenvalue (algebraic connectivity). Furthermore, we demonstrate both analytically and numerically that as the network's algebraic connectivity increases, either by increasing the average node degree of each region or the global migration rate, the global basic reproduction number decreases and converges to the ratio of the average local infection rate to the average local recovery rate, meaning that the lower bound of the global basic reproduction rate does not equal the mean of local basic reproduction rates.

On epidemic spreading in metapopulation networks with time-varying contact patterns

Considering that people may change their face-to-face communication patterns with others depending on the season, we propose an epidemic model that incorporates a time-varying contact rate on a metapopulation network and its second-neighbor network. To describe the time-varying contact mode, we utilize a switched system and define two forms of the basic reproduction number corresponding to two different restrictions. We provide the theoretical proof for the stability of the disease-free equilibrium and confirm periodic stability conditions using simulations. The simulation results reveal that as the period of the switched system lengthens, the amplitude of the final infected density increases; however, the peak infected density within a specific period remains relatively unchanged. Interestingly, as the basic reproduction number grows, the amplitude of the final infected density within a period gradually rises to its maximum and then declines. Moreover, the contact rate that occupies a longer duration within a single period has a more significant influence on epidemic spreading. As the values of different contact rates progressively increase, the recovery rate, natural birth rate, and natural death rate all decrease, leading to a larger final infection density.

Dimension reduction in higher-order contagious phenomena

We investigate epidemic spreading in a deterministic susceptible-infected-susceptible model on uncorrelated heterogeneous networks with higher-order interactions. We provide a recipe for the construction of one-dimensional reduced model (resilience function) of the N-dimensional susceptible-infected-susceptible dynamics in the presence of higher-order interactions. Utilizing this reduction process, we are able to capture the microscopic and macroscopic behavior of infectious networks. We find that the microscopic state of nodes (fraction of stable healthy individual of each node) inversely scales with their degree, and it becomes diminished due to the presence of higher-order interactions. In this case, we analytically obtain that the macroscopic state of the system (fraction of infectious or healthy population) undergoes abrupt transition. Additionally, we quantify the network's resilience, i.e., how the topological changes affect the stable infected population. Finally, we provide an alternative framework of dimension reduction based on the spectral analysis of the network, which can identify the critical onset of the disease in the presence or absence of higher-order interactions. Both reduction methods can be extended for a large class of dynamical models.

Identifying partial topology of simplicial complexes

This paper investigates partial topology identification of simplicial complexes based on adaptive synchronization. For the nodes of interest, the interactions that they participate in can be accurately reconstructed by designing adaptive controllers and parameter estimators. Particularly, not only pairwise interactions but a higher-order structure can be effectively recovered by our method. Moreover, a new linear independence condition with a rigorous definition is established for parameter estimators to converge asymptotically to the true values. Numerical simulations on a general two-dimensional simplicial complex as well as a real-world structure are provided to show the validity of the result and discuss the influence of different parameters on the identification process.

Intervention of resource allocation strategies on spatial spread of epidemics

Medical resources are crucial in mitigating epidemics, especially during pandemics such as the ongoing COVID-19. Thereby, reasonable resource deployment inevitably plays a significant role in suppressing the epidemic under limited resources. When an epidemic breaks out, people can produce resources for self-protection or donate resources to help others for treatment. That is, the exchange of resources also affects the transmission between individuals, thus, altering the epidemic dynamics. To understand factors on resource deployment and the interplay between resource and transmission we construct a metapopulation network model with resource allocation. Our results indicate actively or promptly donating resources is not helpful to suppress the epidemic under both homogeneous population distribution and heterogeneous population distribution. Besides, strengthening the speed of resources production can significantly increase the recovery rate so that they reduce the final outbreak size. These results may provide policy guidance toward epidemic containment.

A metapopulation approach to identify targets for Wolbachia-based dengue control

Over the last decade, the release of Wolbachia-infected Aedes aegypti into the natural habitat of this mosquito species has become the most sustainable and long-lasting technique to prevent and control vector-borne diseases, such as dengue, zika, or chikungunya. However, the limited resources to generate such mosquitoes and their effective distribution in large areas dominated by the Aedes aegypti vector represent a challenge for policymakers. Here, we introduce a mathematical framework for the spread of dengue in which competition between wild and Wolbachia-infected mosquitoes, the cross-contagion patterns between humans and vectors, the heterogeneous distribution of the human population in different areas, and the mobility flows between them are combined. Our framework allows us to identify the most effective areas for the release of Wolbachia-infected mosquitoes to achieve a large decrease in the global dengue prevalence.

Metapopulation epidemic models with a universal mobility pattern on interconnected networks

The global pandemic of the coronavirus disease 2019 (COVID-19) exemplifies the influence of human mobility on epidemic spreading. A framework called the movement-interaction-return (MIR) model is a model to study the impact of human mobility on epidemic spreading. In this paper, we investigate epidemic spreading in interconnected metapopulation networks. Specifically, we incorporate the human mobility pattern called the radiation model into the MIR model. As a result, the proposed model is more realistic in comparison to the original MIR model. We use the tensorial framework to develop Markovian equations that describe the dynamics of the proposed model on interconnected metapopulation networks. Then we derive the corresponding epidemic thresholds by converting tensors into matrices. Comprehensive numerical simulations confirm our analysis.

Contagion-diffusion processes with recurrent mobility patterns of distinguishable agents

The analysis of contagion-diffusion processes in metapopulations is a powerful theoretical tool to study how mobility influences the spread of communicable diseases. Nevertheless, many metapopulation approaches use indistinguishable agents to alleviate analytical difficulties. Here, we address the impact that recurrent mobility patterns, and the spatial distribution of distinguishable agents, have on the unfolding of epidemics in large urban areas. We incorporate the distinguishable nature of agents regarding both their residence and their usual destination. The proposed model allows both a fast computation of the spatiotemporal pattern of the epidemic trajectory and the analytical calculation of the epidemic threshold. This threshold is found as the spectral radius of a mixing matrix encapsulating the residential distribution and the specific commuting patterns of agents. We prove that the simplification of indistinguishable individuals overestimates the value of the epidemic threshold.

Identifying influential subpopulations in metapopulation epidemic models using message-passing theory

Identifying influential subpopulations in metapopulation epidemic models has far-reaching potential implications for surveillance and intervention policies of a global pandemic. However, there is a lack of methods to determine influential nodes in metapopulation models based on a rigorous mathematical background. In this study, we derive the message-passing theory for metapopulation modeling and propose a method to determine influential spreaders. Based on our analysis, we identify the most dangerous city as a potential seed of a pandemic when applied to real-world data. Moreover, we particularly assess the relative importance of various sources of heterogeneity at the subpopulation level, e.g., the number of connections and mobility patterns, to determine properties of spreading processes. We validate our theory with extensive numerical simulations on empirical and synthetic networks considering various mobility and transmission probabilities. We confirm that our theory can accurately predict influential subpopulations in metapopulation models.

Modelling the geographical spread of HIV among MSM in Guangdong, China: a metapopulation model considering the impact of pre-exposure prophylaxis

Men who have sex with men (MSM) make up the majority of new human immunodeficiency virus (HIV) diagnoses among young people in China. Understanding HIV transmission dynamics among the MSM population is, therefore, crucial for the control and prevention of HIV infections, especially for some newly reported genotypes of HIV. This study presents a metapopulation model considering the impact of pre-exposure prophylaxis (PrEP) to investigate the geographical spread of a hypothetically new genotype of HIV among MSM in Guangdong, China. We use multiple data sources to construct this model to characterize the behavioural dynamics underlying the spread of HIV within and between 21 prefecture-level cities (i.e. Guangzhou, Shenzhen, Foshan, etc.) in Guangdong province: the online social network via a gay social networking app, the offline human mobility network via the Baidu mobility website, and self-reported sexual behaviours among MSM. Results show that PrEP initiation exponentially delays the occurrence of the virus for the rest of the cities transmitted from the initial outbreak city; hubs on the movement network, such as Guangzhou, Shenzhen, and Foshan are at a higher risk of 'earliest' exposure to the new HIV genotype; most cities acquire the virus directly from the initial outbreak city while others acquire the virus from cities that are not initial outbreak locations and have relatively high betweenness centralities, such as Guangzhou, Shenzhen and Shantou. This study provides insights in predicting the geographical spread of a new genotype of HIV among an MSM population from different regions and assessing the importance of prefecture-level cities in the control and prevention of HIV in Guangdong province. This article is part of the theme issue 'Data science approach to infectious disease surveillance'.

Phase transitions in virology

Viruses have established relationships with almost every other living organism on Earth and at all levels of biological organization: from other viruses up to entire ecosystems. In most cases, they peacefully coexist with their hosts, but in most relevant cases, they parasitize them and induce diseases and pandemics, such as the AIDS and the most recent avian influenza and COVID-19 pandemic events, causing a huge impact on health, society, and economy. Viruses play an essential role in shaping the eco-evolutionary dynamics of their hosts, and have been also involved in some of the major evolutionary innovations either by working as vectors of genetic information or by being themselves coopted by the host into their genomes. Viruses can be studied at different levels of biological organization, from the molecular mechanisms of genome replication, gene expression and encapsidation, to global pandemics. All these levels are different and yet connected through the presence of threshold conditions allowing for the formation of a capsid, the loss of genetic information or epidemic spreading. These thresholds, as occurs with temperature separating phases in a liquid, define sharp qualitative types of behaviour. Thesephase transitionsare very well known in physics. They have been studied by means of simple, but powerful models able to capture their essential properties, allowing us to better understand them. Can the physics of phase transitions be an inspiration for our understanding of viral dynamics at different scales? Here we review well-known mathematical models of transition phenomena in virology. We suggest that the advantages of abstract, simplified pictures used in physics are also the key to properly understanding the origins and evolution of complexity in viruses. By means of several examples, we explore this multilevel landscape and how minimal models provide deep insights into a diverse array of problems. The relevance of these transitions in connecting dynamical patterns across scales and their evolutionary and clinical implications are outlined.

Synchronizability of double-layer dumbbell networks

Synchronization of multiplex networks has been a topical issue in network science. Dumbbell networks are very typical structures in complex networks which are distinguished from both regular star networks and general community structures, whereas the synchronous dynamics of a double-layer dumbbell network relies on the interlink patterns between layers. In this paper, two kinds of double-layer dumbbell networks are defined according to different interlayer coupling patterns: one with the single-link coupling pattern between layers and the other with the two-link coupling pattern between layers. Furthermore, the largest and smallest nonzero eigenvalues of the Laplacian matrix are calculated analytically and numerically for the single-link coupling pattern and also obtained numerically for the two-link coupling pattern so as to characterize the synchronizability of double-layer dumbbell networks. It is shown that interlayer coupling patterns have a significant impact on the synchronizability of multiplex systems. Finally, a numerical example is provided to verify the effectiveness of theoretical analysis. Our findings can facilitate company managers to select optimal interlayer coupling patterns and to assign proper parameters in terms of improving the efficiency and reducing losses of the whole team.

An overview of epidemic models with phase transitions to absorbing states running on top of complex networks

Dynamical systems running on the top of complex networks have been extensively investigated for decades. But this topic still remains among the most relevant issues in complex network theory due to its range of applicability. The contact process (CP) and the susceptible-infected-susceptible (SIS) model are used quite often to describe epidemic dynamics. Despite their simplicity, these models are robust to predict the kernel of real situations. In this work, we review concisely both processes that are well-known and very applied examples of models that exhibit absorbing-state phase transitions. In the epidemic scenario, individuals can be infected or susceptible. A phase transition between a disease-free (absorbing) state and an active stationary phase (where a fraction of the population is infected) are separated by an epidemic threshold. For the SIS model, the central issue is to determine this epidemic threshold on heterogeneous networks. For the CP model, the main interest is to relate critical exponents with statistical properties of the network.

Modeling partial lockdowns in multiplex networks using partition strategies

National stay-at-home orders, or lockdowns, were imposed in several countries to drastically reduce the social interactions mainly responsible for the transmission of the SARS-CoV-2 virus. Despite being essential to slow down the COVID-19 pandemic, these containment measures are associated with an economic burden. In this work, we propose a network approach to model the implementation of a partial lockdown, breaking the society into disconnected components, or partitions. Our model is composed by two main ingredients: a multiplex network representing human contacts within different contexts, formed by a Household layer, a Work layer, and a Social layer including generic social interactions, and a Susceptible-Infected-Recovered process that mimics the epidemic spreading. We compare different partition strategies, with a twofold aim: reducing the epidemic outbreak and minimizing the economic cost associated to the partial lockdown. We also show that the inclusion of unconstrained social interactions dramatically increases the epidemic spreading, while different kinds of restrictions on social interactions help in keeping the benefices of the network partition.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s41109-021-00366-7.

Network structure-based interventions on spatial spread of epidemics in metapopulation networks

Mathematical modeling of epidemics is fundamental to understand the mechanism of the disease outbreak and provides helpful indications for effectiveness of interventions for policy makers. The metapopulation network model has been used in the analysis of epidemic dynamics by taking individual migration between patches into account. However, so far, most of the previous studies unrealistically assume that transmission rates within patches are the same, neglecting the nonuniformity of intervention measures in hindering epidemics. Here, based on the assumption that interventions deployed in a patch depend on its population size or economic level, which have shown a positive correlation with the patch's degree in networks, we propose a metapopulation network model to explore a network structure-based intervention strategy, aiming at understanding the interplay between intervention strategy and other factors including mobility patterns, initial population, as well as the network structure. Our results demonstrate that interventions to patches with different intensity are able to suppress the epidemic spreading in terms of both the epidemic threshold and the final epidemic size. Specifically, the intervention strategy targeting the patches with high degree is able to efficiently suppress epidemics. In addition, a detrimental effect is also observed depending on the interplay between the intervention measures and the initial population distribution. Our study opens a path for understanding epidemic dynamics and provides helpful insights into the implementation of countermeasures for the control of epidemics in reality.

Travel restrictions during pandemics: A useful strategy?

Though carrying considerable economic and societal costs, restricting individuals' traveling freedom appears as a logical way to curb the spreading of an epidemic. However, whether, under what conditions, and to what extent travel restrictions actually exert a mitigating effect on epidemic spreading are poorly understood issues. Recent studies have actually suggested the opposite, i.e., that allowing some movements can hinder the propagation of a disease. Here, we explore this topic by modeling the spreading of a generic contagious disease where susceptible, infected, or recovered point-wise individuals are uncorrelated random-walkers evolving within a space comprising two equally sized separated compartments. We evaluate the spreading process under different separation conditions between the two spatial compartments and a forced relocation schedule. Our results confirm that, under certain conditions, allowing individuals to move from regions of high to low infection rates may turn out to have a positive effect on aggregate; such positive effect is nevertheless reduced if a directional flow is allowed. This highlights the importance of considering travel restriction policies alternative to classical ones.

Infection spreading and recovery in a square lattice

We investigate spreading and recovery of disease in a square lattice, and, in particular, emphasize the role of the initial distribution of infected patches in the network on the progression of an endemic and initiation of a recovery process, if any, due to migration of both the susceptible and infected hosts. The disease starts in the lattice with three possible initial distribution patterns of infected and infection-free sites, viz., infected core patches (ICP), infected peripheral patches (IPP), and randomly distributed infected patches (RDIP). Our results show that infection spreads monotonically in the lattice with increasing migration without showing any sign of recovery in the ICP case. In the IPP case, it follows a similar monotonic progression with increasing migration; however, a self-organized healing process starts for higher migration, leading the lattice to full recovery at a critical rate of migration. Encouragingly, for the initial RDIP arrangement, chances of recovery are much higher with a lower rate of critical migration. An eigenvalue-based semianalytical study is made to determine the critical migration rate for realizing a stable infection-free lattice. The initial fraction of infected patches and the force of infection play significant roles in the self-organized recovery. They follow an exponential law, for the RDIP case, that governs the recovery process. For the frustrating case of ICP arrangement, we propose a random rewiring of links in the lattice allowing long-distance migratory paths that effectively initiate a recovery process. Global prevalence of infection thereby declines and progressively improves with the rewiring probability that follows a power law with the critical migration and leads to the birth of emergent infection-free networks.

Lessons from being challenged by COVID-19

We present results of different approaches to model the evolution of the COVID-19 epidemic in Argentina, with a special focus on the megacity conformed by the city of Buenos Aires and its metropolitan area, including a total of 41 districts with over 13 million inhabitants. We first highlight the relevance of interpreting the early stage of the epidemic in light of incoming infectious travelers from abroad. Next, we critically evaluate certain proposed solutions to contain the epidemic based on instantaneous modifications of the reproductive number. Finally, we build increasingly complex and realistic models, ranging from simple homogeneous models used to estimate local reproduction numbers, to fully coupled inhomogeneous (deterministic or stochastic) models incorporating mobility estimates from cell phone location data. The models are capable of producing forecasts highly consistent with the official number of cases with minimal parameter fitting and fine-tuning. We discuss the strengths and limitations of the proposed models, focusing on the validity of different necessary first approximations, and caution future modeling efforts to exercise great care in the interpretation of long-term forecasts, and in the adoption of non-pharmaceutical interventions backed by numerical simulations.

Spreading of two interacting diseases in multiplex networks

We consider the interacting processes between two diseases on multiplex networks, where each node can be infected by two interacting diseases with general interacting schemes. A discrete-time individual-based probability model is rigorously derived. By the bifurcation analysis of the equilibrium, we analyze the outbreak condition of one disease. The theoretical predictions are in good agreement with discrete-time stochastic simulations on scale-free networks. Furthermore, we discuss the influence of network overlap and dynamical parameters on the epidemic dynamical behaviors. The simulation results show that the network overlap has almost no effect on both epidemic threshold and prevalence. We also find that the epidemic threshold of one disease does not depend on all system parameters. Our method offers an analytical framework for the spreading dynamics of multiple processes in multiplex networks.

Intervention threshold for epidemic control in susceptible-infected-recovered metapopulation models

Metapopulation epidemic models describe epidemic dynamics in networks of spatially distant patches connected via pathways for migration of individuals. In the present study, we deal with a susceptible-infected-recovered (SIR) metapopulation model where the epidemic process in each patch is represented by an SIR model and the mobility of individuals is assumed to be a homogeneous diffusion. We consider two types of patches including high-risk and low-risk ones under the assumption that a local patch is changed from a high-risk one to a low-risk one by an intervention. We theoretically analyze the intervention threshold which indicates the critical fraction of low-risk patches for preventing a global epidemic outbreak. We show that an intervention targeted to high-degree patches is more effective for epidemic control than a random intervention. The theoretical results are validated by Monte Carlo simulations for synthetic and realistic scale-free patch networks. The theoretical results also reveal that the intervention threshold depends on the human mobility network and the mobility rate. Our approach is useful for exploring better local interventions aimed at containment of epidemics.

Complex Networks: a Mini-review

Network analysis is a powerful tool that provides us a fruitful framework to describe phenomena related to social, technological, and many other real-world complex systems. In this paper, we present a brief review about complex networks including fundamental quantities, examples of network models, and the essential role of network topology in the investigation of dynamical processes as epidemics, rumor spreading, and synchronization. A quite of advances have been provided in this field, and many other authors also review the main contributions in this area over the years. However, we show an overview from a different perspective. Our aim is to provide basic information to a broad audience and more detailed references for those who would like to learn deeper the topic.

Analytical treatment for cyclic three-state dynamics on static networks

Whenever a dynamical process unfolds on static networks, the dynamical state of any focal individual will be exclusively influenced by directly connected neighbors, rather than by those unconnected ones, hence the arising of the dynamical correlation problem, where mean-field-based methods fail to capture the scenario. The dynamic correlation coupling problem has always been an important and difficult problem in the theoretical field of physics. The explicit analytical expressions and the decoupling methods often play a key role in the development of corresponding field. In this paper, we study the cyclic three-state dynamics on static networks, which include a wide class of dynamical processes, for example, the cyclic Lotka-Volterra model, the directed migration model, the susceptible-infected-recovered-susceptible epidemic model, and the predator-prey with empty sites model. We derive the explicit analytical solutions of the propagating size and the threshold curve surface for the four different dynamics. We compare the results on static networks with those on annealed networks and made an interesting discovery: for the symmetrical dynamical model (the cyclic Lotka-Volterra model and the directed migration model, where the three states are of rotational symmetry), the macroscopic behaviors of the dynamical processes on static networks are the same as those on annealed networks; while the outcomes of the dynamical processes on static networks are different with, and more complicated than, those on annealed networks for asymmetric dynamical model (the susceptible-infected-recovered-susceptible epidemic model and the predator-prey with empty sites model). We also compare the results forecasted by our theoretical method with those by Monte Carlo simulations and find good agreement between the results obtained by the two methods.

Markovian approach to tackle the interaction of simultaneous diseases

The simultaneous emergence of several abrupt disease outbreaks or the extinction of some serotypes of multistrain diseases are fingerprints of the interaction between pathogens spreading within the same population. Here, we propose a general and versatile benchmark to address the unfolding of both cooperative and competitive interacting diseases. We characterize the explosive transitions between the disease-free and the epidemic regimes arising from the cooperation between pathogens and show the critical degree of cooperation needed for the onset of such abrupt transitions. For the competing diseases, we characterize the mutually exclusive case and derive analytically the transition point between the full-dominance phase, in which only one pathogen propagates, and the coexistence regime. Finally, we use this framework to analyze the behavior of the former transition point as the competition between pathogens is relaxed.

Explosive transitions induced by interdependent contagion-consensus dynamics in multiplex networks

We introduce a model to study the interplay between information spreading and opinion formation in social systems. Our framework consists in a two-layer multiplex network where opinion dynamics takes place in one layer, while information spreads on the other one. The two dynamical processes are mutually coupled in such a way that the control parameters governing the dynamics of the node states at one layer depend on the dynamical states at the other layer. In particular, we consider the case in which consensus is favored by the common adoption of information, while information spreading is boosted between agents sharing similar opinions. Numerical simulations of the model point out that, when the coupling between the dynamics of the two layers is strong enough, a double explosive transition, i.e., a discontinuous transition both in consensus dynamics and in information spreading appears. Such explosive transitions lead to bi-stability regions in which the consensus-informed states and the disagreement-uninformed states are both stable solutions of the intertwined dynamics.