Bursts of vertex activation and epidemics in evolving networks

Structural inequalities exacerbate infection disparities

During the COVID-19 pandemic, the world witnessed a disproportionate infection rate among marginalized and low-income groups. Despite empirical evidence suggesting that structural inequalities in society contribute to health disparities, there has been little attempt to offer a computational and theoretical explanation to establish its plausibility and quantitative impact. Here, we focus on two aspects of structural inequalities: wealth inequality and social segregation. Our computational model demonstrates that (a) due to the inequality in self-quarantine ability, the infection gap widens between the low-income and high-income groups, and the overall infected cases increase, (b) social segregation between different socioeconomic status (SES) groups intensifies the disease spreading rates, and (c) the second wave of infection can emerge due to a false sense of safety among the medium and high SES groups. By performing two data-driven analyses, one on the empirical network and economic data of 404 metropolitan areas of the United States and one on the daily Covid-19 data of the City of Chicago, we verify that higher segregation leads to an increase in the overall infection cases and higher infection inequality across different ethnic/socioeconomic groups. These findings together demonstrate that reducing structural inequalities not only helps decrease health disparities but also reduces the spread of infectious diseases overall.

Temporal scaling theory for bursty time series with clusters of arbitrarily many events

Long-term temporal correlations in time series in a form of an event sequence have been characterized using an autocorrelation function that often shows a power-law decaying behavior. Such scaling behavior has been mainly accounted for by the heavy-tailed distribution of interevent times, i.e., the time interval between two consecutive events. Yet, little is known about how correlations between consecutive interevent times systematically affect the decaying behavior of the autocorrelation function. Empirical distributions of the burst size, which is the number of events in a cluster of events occurring in a short time window, often show heavy tails, implying that arbitrarily many consecutive interevent times may be correlated with each other. In the present study, we propose a model for generating a time series with arbitrary functional forms of interevent time and burst size distributions. Then, we analytically derive the autocorrelation function for the model time series. In particular, by assuming that the interevent time and burst size are power-law distributed, we derive scaling relations between power-law exponents of the autocorrelation function decay, interevent time distribution, and burst size distribution. These analytical results are confirmed by numerical simulations. Our approach helps to rigorously and analytically understand the effects of correlations between arbitrarily many consecutive interevent times on the decaying behavior of the autocorrelation function.

Bursts of communication increase opinion diversity in the temporal Deffuant model

Human interactions create social networks forming the backbone of societies. Individuals adjust their opinions by exchanging information through social interactions. Two recurrent questions are whether social structures promote opinion polarisation or consensus and whether polarisation can be avoided, particularly on social media. In this paper, we hypothesise that not only network structure but also the timings of social interactions regulate the emergence of opinion clusters. We devise a temporal version of the Deffuant opinion model where pairwise social interactions follow temporal patterns. Individuals may self-organise into a multi-partisan society due to network clustering promoting the reinforcement of local opinions. Burstiness has a similar effect and is alone sufficient to refrain the population from consensus and polarisation by also promoting the reinforcement of local opinions. The diversity of opinions in socially clustered networks thus increases with burstiness, particularly, and counter-intuitively, when individuals have low tolerance and prefer to adjust to similar peers. The emergent opinion landscape is well-balanced regarding groups' size, with relatively short differences between groups, and a small fraction of extremists. We argue that polarisation is more likely to emerge in social media than offline social networks because of the relatively low social clustering observed online, despite the observed online burstiness being sufficient to promote more diversity than would be expected offline. Increasing the variance of burst activation times, e.g. by being less active on social media, could be a venue to reduce polarisation. Furthermore, strengthening online social networks by increasing social redundancy, i.e. triangles, may also promote diversity.

The use of networks in spatial and temporal computational models for outbreak spread in epidemiology: A systematic review

OBJECTIVES

To examine recent literature in order to present a comprehensive overview of the current trends as regards the computational models used to represent the propagation of an infectious outbreak in a population, paying particular attention to those that represent network-based transmission.

METHODS

a systematic review was conducted following the PRISMA (Preferred Reporting Items for Systematic Reviews and Meta-Analyses) guidelines. Papers published in English between 2010 and September 2021 were sought in the ACM Digital Library, IEEE Xplore, PubMed and Scopus databases.

RESULTS

Upon considering their titles and abstracts, 832 papers were obtained, of which 192 were selected for a full content-body check. Of these, 112 studies were eventually deemed suitable for quantitative and qualitative analysis. Emphasis was placed on the spatial and temporal scales studied, the use of networks or graphs, and the granularity of the data used to evaluate the models. The models principally used to represent the spreading of outbreaks have been stochastic (55.36%), while the type of networks most frequently used are relationship networks (32.14%). The most common spatial dimension used is a region (19.64%) and the most used unit of time is a day (28.57%). Synthetic data as opposed to an external source were used in 51.79% of the papers. With regard to the granularity of the data sources, aggregated data such as censuses or transportation surveys are the most common.

CONCLUSION

We identified a growing interest in the use of networks to represent disease transmission. We detected that research is focused on only certain combinations of the computational model, type of network (in both the expressive and the structural sense) and spatial scale, while the search for other interesting combinations has been left for the future.

Routing Strategies for Isochronal-Evolution Random Matching Network

In order to abstract away a network model from some real-world networks, such as navigation satellite networks and mobile call networks, we proposed an Isochronal-Evolution Random Matching Network (IERMN) model. An IERMN is a dynamic network that evolves isochronally and has a collection of edges that are pairwise disjoint at any point in time. We then investigated the traffic dynamics in IERMNs whose main research topic is packet transmission. When a vertex of an IERMN plans a path for a packet, it is permitted to delay the sending of the packet to make the path shorter. We designed a routing decision-making algorithm for vertices based on replanning. Since the IERMN has a specific topology, we developed two suitable routing strategies: the Least Delay Path with Minimum Hop (LDPMH) routing strategy and the Least Hop Path with Minimum Delay (LHPMD) routing strategy. An LDPMH is planned by a binary search tree and an LHPMD is planned by an ordered tree. The simulation results show that the LHPMD routing strategy outperformed the LDPMH routing strategy in terms of the critical packet generation rate, number of delivered packets, packet delivery ratio, and average posterior path lengths.

Impact of human contact patterns on epidemic spreading in time-varying networks

Human contact behaviors involve both dormant and active processes. The dormant (active) process goes from the disappearance (creation) to the creation (disappearance) of an edge. The dormant (active) time is the elapsed time since the edge became dormant (active). Many empirical studies have revealed that dormant and active times in human contact behaviors tend to show a long-tailed distribution. Previous researches focused on the impact of the dormant process on spreading dynamics. However, the epidemic spreading happens on the active process. This raises the question of how the active process affects epidemic spreading in complex networks. Here, we propose a novel time-varying network model in which the distributions of both the dormant time and active time of edges are adjustable. We develop a pairwise approximation method to describe the spreading dynamical processes in the time-varying networks. Through extensive numerical simulations, we find that the epidemic threshold is proportional to the mean dormant time and inversely proportional to the mean active time. The attack rate decreases with the increase of mean dormant time and increases with the increase of mean active time. It is worth noting that the epidemic threshold and the attack rate (e.g., the infected density in the steady state) are independent of the heterogeneities of the dormant time distribution and the active time distribution. Increasing the heterogeneity of the dormant time distribution accelerates epidemic spreading while increasing the heterogeneity of the active time distribution slows it down.

The interplay between disease spreading and awareness diffusion in multiplex networks with activity-driven structure

The interplay between disease and awareness has been extensively studied in static networks. However, most networks in reality will evolve over time. Based on this, we propose a novel epidemiological model in multiplex networks. In this model, the disease spreading layer is a time-varying network generated by the activity-driven model, while the awareness diffusion layer is a static network, and the heterogeneity of individual infection and recovery ability is considered. First, we extend the microscopic Markov chain approach to analytically obtain the epidemic threshold of the model. Then, we simulate the spread of disease and find that stronger heterogeneity in the individual activities of a physical layer can promote disease spreading, while stronger heterogeneity of the virtual layer network will hinder the spread of disease. Interestingly, we find that when the individual infection ability follows Gaussian distribution, the heterogeneity of infection ability has little effect on the spread of disease, but it will significantly affect the epidemic threshold when the individual infection ability follows power-law distribution. Finally, we find the emergence of a metacritical point where the diffusion of awareness is able to control the onset of the epidemics. Our research could cast some light on exploring the dynamics of epidemic spreading in time-varying multiplex networks.

The dynamical formation of ephemeral groups on networks and their effects on epidemics spreading

In network models of propagation processes, the individual, microscopic level perspective is the norm, with aggregations studied as possible outcomes. On the contrary, we adopted a mesoscale perspective with groups as the core element and in this sense we present a novel agent-group dynamic model of propagation in networks. In particular, we focus on ephemeral groups that dynamically form, create new links, and dissolve. The experiments simulated 160 model configurations and produced results describing cases of consecutive and non-consecutive dynamic grouping, bounded or unbounded in the number of repetitions. Results revealed the existence of complex dynamics and multiple behaviors. An efficiency metric is introduced to compare the different cases. A Null Model analysis disclosed a pattern in the difference between the group and random models, varying with the size of groups. Our findings indicate that a mesoscopic construct like the ephemeral group, based on assumptions about social behavior and absent any microscopic level change, could produce and describe complex propagation dynamics. A conclusion is that agent-group dynamic models may represent a powerful approach for modelers and a promising new direction for future research in models of coevolution between propagation and behavior in society.

Social fluidity mobilizes contagion in human and animal populations

Humans and other group-living animals tend to distribute their social effort disproportionately. Individuals predominantly interact with a small number of close companions while maintaining weaker social bonds with less familiar group members. By incorporating this behavior into a mathematical model, we find that a single parameter, which we refer to as social fluidity, controls the rate of social mixing within the group. Large values of social fluidity correspond to gregarious behavior, whereas small values signify the existence of persistent bonds between individuals. We compare the social fluidity of 13 species by applying the model to empirical human and animal social interaction data. To investigate how social behavior influences the likelihood of an epidemic outbreak, we derive an analytical expression of the relationship between social fluidity and the basic reproductive number of an infectious disease. For species that form more stable social bonds, the model describes frequency-dependent transmission that is sensitive to changes in social fluidity. As social fluidity increases, animal-disease systems become increasingly density-dependent. Finally, we demonstrate that social fluidity is a stronger predictor of disease outcomes than both group size and connectivity, and it provides an integrated framework for both density-dependent and frequency-dependent transmission.

Impact of temporal correlations on high risk outbreaks of independent and cooperative SIR dynamics

We first propose a quantitative approach to detect high risk outbreaks of independent and coinfective SIR dynamics on three empirical networks: a school, a conference and a hospital contact network. This measurement is based on the k-means clustering method and identifies proper samples for calculating the mean outbreak size and the outbreak probability. Then we systematically study the impact of different temporal correlations on high risk outbreaks over the original and differently shuffled counterparts of each network. We observe that, on the one hand, in the coinfection process, randomization of the sequence of the events increases the mean outbreak size of high-risk cases. On the other hand, these correlations do not have a consistent effect on the independent infection dynamics, and can either decrease or increase this mean. Randomization of the daily pattern correlations has no strong impact on the size of the outbreak in either the coinfection or the independent spreading cases. We also observe that an increase in the mean outbreak size does not always coincide with an increase in the outbreak probability; therefore, we argue that merely considering the mean outbreak size of all realizations may lead us into falsely estimating the outbreak risks. Our results suggest that some sort of contact randomization in the organizational level in schools, events or hospitals might help to suppress the spreading dynamics while the risk of an outbreak is high.

Spectral Analysis of Epidemic Thresholds of Temporal Networks

Many complex systems can be modeled as temporal networks with time-evolving connections. The influence of their characteristics on epidemic spreading is analyzed in a susceptible-infected-susceptible epidemic model illustrated by the discrete-time Markov chain approach. We develop the analytical epidemic thresholds in terms of the spectral radius of weighted adjacency matrix by averaging temporal networks, e.g., periodic, nonperiodic Markovian networks, and a special nonperiodic non-Markovian network (the link activation network) in time. We discuss the impacts of statistical characteristics, e.g., bursts and duration heterogeneity, as well as time-reversed characteristic on epidemic thresholds. We confirm the tightness of the proposed epidemic thresholds with numerical simulations on seven artificial and empirical temporal networks and show that the epidemic threshold of our theory is more precise than those of previous studies.

Spreading of computer viruses on time-varying networks

Social networks are the prime channel for the spreading of computer viruses. Yet the study of their propagation neglects the temporal nature of social interactions and the heterogeneity of users' susceptibility. Here, we introduce a theoretical framework that captures both properties. We study two realistic types of viruses propagating on temporal networks featuring Q categories of susceptibility and derive analytically the invasion threshold. We found that the temporal coupling of categories might increase the fragility of the system to cyber threats. Our results show that networks' dynamics and their interplay with users' features are crucial for the spreading of computer viruses.

Information diffusion backbones in temporal networks

Progress has been made in understanding how temporal network features affect the percentage of nodes reached by an information diffusion process. In this work, we explore further: which node pairs are likely to contribute to the actual diffusion of information, i.e., appear in a diffusion trajectory? How is this likelihood related to the local temporal connection features of the node pair? Such deep understanding of the role of node pairs is crucial to tackle challenging optimization problems such as which kind of node pairs or temporal contacts should be stimulated in order to maximize the prevalence of information spreading. We start by using Susceptible-Infected (SI) model, in which an infected (information possessing) node could spread the information to a susceptible node with a given infection probability β whenever a contact happens between the two nodes, as the information diffusion process. We consider a large number of real-world temporal networks. First, we propose the construction of an information diffusion backbone GB(β) for a SI spreading process with an infection probability β on a temporal network. The backbone is a weighted network where the weight of each node pair indicates how likely the node pair appears in a diffusion trajectory starting from an arbitrary node. Second, we investigate the relation between the backbones with different infection probabilities on a temporal network. We find that the backbone topology obtained for low and high infection probabilities approach the backbone GB(β → 0) and GB(β = 1), respectively. The backbone GB(β → 0) equals the integrated weighted network, where the weight of a node pair counts the total number of contacts in between. Finally, we explore node pairs with what local connection features tend to appear in GB(β = 1), thus actually contribute to the global information diffusion. We discover that a local connection feature among many other features we proposed, could well identify the (high-weight) links in GB(β = 1). This local feature encodes the time that each contact occurs, pointing out the importance of temporal features in determining the role of node pairs in a dynamic process.

Epidemic spreading on time-varying multiplex networks

Social interactions are stratified in multiple contexts and are subject to complex temporal dynamics. The systematic study of these two features of social systems has started only very recently, mainly thanks to the development of multiplex and time-varying networks. However, these two advancements have progressed almost in parallel with very little overlap. Thus, the interplay between multiplexity and the temporal nature of connectivity patterns is poorly understood. Here, we aim to tackle this limitation by introducing a time-varying model of multiplex networks. We are interested in characterizing how these two properties affect contagion processes. To this end, we study susceptible-infected-susceptible epidemic models unfolding at comparable timescale with respect to the evolution of the multiplex network. We study both analytically and numerically the epidemic threshold as a function of the multiplexity and the features of each layer. We found that higher values of multiplexity significantly reduce the epidemic threshold especially when the temporal activation patterns of nodes present on multiple layers are positively correlated. Furthermore, when the average connectivity across layers is very different, the contagion dynamics is driven by the features of the more densely connected layer. Here, the epidemic threshold is equivalent to that of a single layered graph and the impact of the disease, in the layer driving the contagion, is independent of the multiplexity. However, this is not the case in the other layers where the spreading dynamics is sharply influenced by it. The results presented provide another step towards the characterization of the properties of real networks and their effects on contagion phenomena.

Correlated bursts in temporal networks slow down spreading

Spreading dynamics has been considered to take place in temporal networks, where temporal interaction patterns between nodes show non-Poissonian bursty nature. The effects of inhomogeneous interevent times (IETs) on the spreading have been extensively studied in recent years, yet little is known about the effects of correlations between IETs on the spreading. In order to investigate those effects, we study two-step deterministic susceptible-infected (SI) and probabilistic SI dynamics when the interaction patterns are modeled by inhomogeneous and correlated IETs, i.e., correlated bursts. By analyzing the transmission time statistics in a single-link setup and by simulating the spreading in Bethe lattices and random graphs, we conclude that the positive correlation between IETs slows down the spreading. We also argue that the shortest transmission time from one infected node to its susceptible neighbors can successfully explain our numerical results.

Exponential random graph models to evaluate the movement of backyard chickens after the avian influenza crisis in 2004-2005, Thailand

The aim of this study was to use exponential random graph models (ERGMs) to explain networks of movement of backyard chickens in provinces which had been hotspots for avian influenza outbreaks in Thailand during 2004-2005. We used structured questionnaires to collect data for the period January to December 2009 from participants who were involved in the backyard chicken farming network in three avian influenza hotspots (Ratchaburi, Suphan Buri, and Nakhon Pathom provinces) in Thailand. From 557 questionnaires, we identified nodes, points of entry and exit from nodes, and activities relating to backyard chicken farming and movement of chickens, and generated ERGMs based on non-festive periods (Model 1) and the Chinese New Year period (Model 2). In Model 1, k-star (the central node is connected to k other nodes) connections were predominant (P <  0.001). In Model 2, the frequency of movement increased by 10.62 times, k-star connections were still predominant (P <  0.001), and the model was scale-free. Hubs were formed from owners/observers in the arenas/training fields, farmers who raised chickens for consumption only, and traders. In conclusion, our models indicated that, if avian influenza was introduced during non-festive periods, the authorities would need to regularly restrict the movement of chickens. However, during high-frequency periods of movement of backyard chickens, authorities would also need to focus on the network hubs. Our research can be used by the relevant authorities to improve control measures and reduce the risk or lessen the magnitude of disease spread during an avian influenza epidemic.

Automated monitoring of behavior reveals bursty interaction patterns and rapid spreading dynamics in honeybee social networks

Social networks mediate the spread of information and disease. The dynamics of spreading depends, among other factors, on the distribution of times between successive contacts in the network. Heavy-tailed (bursty) time distributions are characteristic of human communication networks, including face-to-face contacts and electronic communication via mobile phone calls, email, and internet communities. Burstiness has been cited as a possible cause for slow spreading in these networks relative to a randomized reference network. However, it is not known whether burstiness is an epiphenomenon of human-specific patterns of communication. Moreover, theory predicts that fast, bursty communication networks should also exist. Here, we present a high-throughput technology for automated monitoring of social interactions of individual honeybees and the analysis of a rich and detailed dataset consisting of more than 1.2 million interactions in five honeybee colonies. We find that bees, like humans, also interact in bursts but that spreading is significantly faster than in a randomized reference network and remains so even after an experimental demographic perturbation. Thus, while burstiness may be an intrinsic property of social interactions, it does not always inhibit spreading in real-world communication networks. We anticipate that these results will inform future models of large-scale social organization and information and disease transmission, and may impact health management of threatened honeybee populations.

Epidemic Threshold in Continuous-Time Evolving Networks

Current understanding of the critical outbreak condition on temporal networks relies on approximations (time scale separation, discretization) that may bias the results. We propose a theoretical framework to compute the epidemic threshold in continuous time through the infection propagator approach. We introduce the weak commutation condition allowing the interpretation of annealed networks, activity-driven networks, and time scale separation into one formalism. Our work provides a coherent connection between discrete and continuous time representations applicable to realistic scenarios.

Lifetime-preserving reference models for characterizing spreading dynamics on temporal networks

To study how a certain network feature affects processes occurring on a temporal network, one often compares properties of the original network against those of a randomized reference model that lacks the feature in question. The randomly permuted times (PT) reference model is widely used to probe how temporal features affect spreading dynamics on temporal networks. However, PT implicitly assumes that edges and nodes are continuously active during the network sampling period - an assumption that does not always hold in real networks. We systematically analyze a recently-proposed restriction of PT that preserves node lifetimes (PTN), and a similar restriction (PTE) that also preserves edge lifetimes. We use PT, PTN, and PTE to characterize spreading dynamics on (i) synthetic networks with heterogeneous edge lifespans and tunable burstiness, and (ii) four real-world networks, including two in which nodes enter and leave the network dynamically. We find that predictions of spreading speed can change considerably with the choice of reference model. Moreover, the degree of disparity in the predictions reflects the extent of node/edge turnover, highlighting the importance of using lifetime-preserving reference models when nodes or edges are not continuously present in the network.

Random walks on activity-driven networks with attractiveness

Virtually all real-world networks are dynamical entities. In social networks, the propensity of nodes to engage in social interactions (activity) and their chances to be selected by active nodes (attractiveness) are heterogeneously distributed. Here, we present a time-varying network model where each node and the dynamical formation of ties are characterized by these two features. We study how these properties affect random-walk processes unfolding on the network when the time scales describing the process and the network evolution are comparable. We derive analytical solutions for the stationary state and the mean first-passage time of the process, and we study cases informed by empirical observations of social networks. Our work shows that previously disregarded properties of real social systems, such as heterogeneous distributions of activity and attractiveness as well as the correlations between them, substantially affect the dynamical process unfolding on the network.

Burstiness and tie activation strategies in time-varying social networks

The recent developments in the field of social networks shifted the focus from static to dynamical representations, calling for new methods for their analysis and modelling. Observations in real social systems identified two main mechanisms that play a primary role in networks' evolution and influence ongoing spreading processes: the strategies individuals adopt when selecting between new or old social ties, and the bursty nature of the social activity setting the pace of these choices. We introduce a time-varying network model accounting both for ties selection and burstiness and we analytically study its phase diagram. The interplay of the two effects is non trivial and, interestingly, the effects of burstiness might be suppressed in regimes where individuals exhibit a strong preference towards previously activated ties. The results are tested against numerical simulations and compared with two empirical datasets with very good agreement. Consequently, the framework provides a principled method to classify the temporal features of real networks, and thus yields new insights to elucidate the effects of social dynamics on spreading processes.

Prediction of Pig Trade Movements in Different European Production Systems Using Exponential Random Graph Models

In most European countries, data regarding movements of live animals are routinely collected and can greatly aid predictive epidemic modeling. However, the use of complete movements’ dataset to conduct policy-relevant predictions has been so far limited by the massive amount of data that have to be processed (e.g., in intensive commercial systems) or the restricted availability of timely and updated records on animal movements (e.g., in areas where small-scale or extensive production is predominant). The aim of this study was to use exponential random graph models (ERGMs) to reproduce, understand, and predict pig trade networks in different European production systems. Three trade networks were built by aggregating movements of pig batches among premises (farms and trade operators) over 2011 in Bulgaria, Extremadura (Spain), and Côtes-d’Armor (France), where small-scale, extensive, and intensive pig production are predominant, respectively. Three ERGMs were fitted to each network with various demographic and geographic attributes of the nodes as well as six internal network configurations. Several statistical and graphical diagnostic methods were applied to assess the goodness of fit of the models. For all systems, both exogenous (attribute-based) and endogenous (network-based) processes appeared to govern the structure of pig trade network, and neither alone were capable of capturing all aspects of the network structure. Geographic mixing patterns strongly structured pig trade organization in the small-scale production system, whereas belonging to the same company or keeping pigs in the same housing system appeared to be key drivers of pig trade, in intensive and extensive production systems, respectively. Heterogeneous mixing between types of production also explained a part of network structure, whichever production system considered. Limited information is thus needed to capture most of the global structure of pig trade networks. Such findings will be useful to simplify trade networks analysis and better inform European policy makers on risk-based and more cost-effective prevention and control against swine diseases such as African swine fever, classical swine fever, or porcine reproductive and respiratory syndrome.

Recalibrating disease parameters for increasing realism in modeling epidemics in closed settings

Background

The homogeneous mixing assumption is widely adopted in epidemic modelling for its parsimony and represents the building block of more complex approaches, including very detailed agent-based models. The latter assume homogeneous mixing within schools, workplaces and households, mostly for the lack of detailed information on human contact behaviour within these settings. The recent data availability on high-resolution face-to-face interactions makes it now possible to assess the goodness of this simplified scheme in reproducing relevant aspects of the infection dynamics.

Methods

We consider empirical contact networks gathered in different contexts, as well as synthetic data obtained through realistic models of contacts in structured populations. We perform stochastic spreading simulations on these contact networks and in populations of the same size under a homogeneous mixing hypothesis. We adjust the epidemiological parameters of the latter in order to fit the prevalence curve of the contact epidemic model. We quantify the agreement by comparing epidemic peak times, peak values, and epidemic sizes.

Results

Good approximations of the peak times and peak values are obtained with the homogeneous mixing approach, with a median relative difference smaller than 20 % in all cases investigated. Accuracy in reproducing the peak time depends on the setting under study, while for the peak value it is independent of the setting. Recalibration is found to be linear in the epidemic parameters used in the contact data simulations, showing changes across empirical settings but robustness across groups and population sizes.

Conclusions

An adequate rescaling of the epidemiological parameters can yield a good agreement between the epidemic curves obtained with a real contact network and a homogeneous mixing approach in a population of the same size. The use of such recalibrated homogeneous mixing approximations would enhance the accuracy and realism of agent-based simulations and limit the intrinsic biases of the homogeneous mixing.

Electronic supplementary material

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

Individual-based approach to epidemic processes on arbitrary dynamic contact networks

The dynamics of contact networks and epidemics of infectious diseases often occur on comparable time scales. Ignoring one of these time scales may provide an incomplete understanding of the population dynamics of the infection process. We develop an individual-based approximation for the susceptible-infected-recovered epidemic model applicable to arbitrary dynamic networks. Our framework provides, at the individual-level, the probability flow over time associated with the infection dynamics. This computationally efficient framework discards the correlation between the states of different nodes, yet provides accurate results in approximating direct numerical simulations. It naturally captures the temporal heterogeneities and correlations of contact sequences, fundamental ingredients regulating the timing and size of an epidemic outbreak, and the number of secondary infections. The high accuracy of our approximation further allows us to detect the index individual of an epidemic outbreak in real-life network data.

Information content of contact-pattern representations and predictability of epidemic outbreaks

To understand the contact patterns of a population--who is in contact with whom, and when the contacts happen--is crucial for modeling outbreaks of infectious disease. Traditional theoretical epidemiology assumes that any individual can meet any with equal probability. A more modern approach, network epidemiology, assumes people are connected into a static network over which the disease spreads. Newer yet, temporal network epidemiology, includes the time in the contact representations. In this paper, we investigate the effect of these successive inclusions of more information. Using empirical proximity data, we study both outbreak sizes from unknown sources, and from known states of ongoing outbreaks. In the first case, there are large differences going from a fully mixed simulation to a network, and from a network to a temporal network. In the second case, differences are smaller. We interpret these observations in terms of the temporal network structure of the data sets. For example, a fast overturn of nodes and links seem to make the temporal information more important.

Dynamics of history-dependent epidemics in temporal networks

The structural properties of temporal networks often influence the dynamical processes that occur on these networks, e.g., bursty interaction patterns have been shown to slow down epidemics. In this paper, we investigate the effect of link lifetimes on the spread of history-dependent epidemics. We formulate an analytically tractable activity-driven temporal network model that explicitly incorporates link lifetimes. For Markovian link lifetimes, we use mean-field analysis for computing the epidemic threshold, while the effect of non-Markovian link lifetimes is studied using simulations. Furthermore, we also study the effect of negative correlation between the number of links spawned by an individual and the lifetimes of those links. Such negative correlations may arise due to the finite cognitive capacity of the individuals. Our investigations reveal that heavy-tailed link lifetimes slow down the epidemic, while negative correlations can reduce epidemic prevalence. We believe that our results help shed light on the role of link lifetimes in modulating diffusion processes on temporal networks.

Diffusion on networked systems is a question of time or structure

Network science investigates the architecture of complex systems to understand their functional and dynamical properties. Structural patterns such as communities shape diffusive processes on networks. However, these results hold under the strong assumption that networks are static entities where temporal aspects can be neglected. Here we propose a generalized formalism for linear dynamics on complex networks, able to incorporate statistical properties of the timings at which events occur. We show that the diffusion dynamics is affected by the network community structure and by the temporal properties of waiting times between events. We identify the main mechanism--network structure, burstiness or fat tails of waiting times--determining the relaxation times of stochastic processes on temporal networks, in the absence of temporal-structure correlations. We identify situations when fine-scale structure can be discarded from the description of the dynamics or, conversely, when a fully detailed model is required due to temporal heterogeneities.

The basic reproduction number as a predictor for epidemic outbreaks in temporal networks

The basic reproduction number R₀—the number of individuals directly infected by an infectious person in an otherwise susceptible population—is arguably the most widely used estimator of how severe an epidemic outbreak can be. This severity can be more directly measured as the fraction of people infected once the outbreak is over, Ω. In traditional mathematical epidemiology and common formulations of static network epidemiology, there is a deterministic relationship between R₀ and Ω. However, if one considers disease spreading on a temporal contact network—where one knows when contacts happen, not only between whom—then larger R₀ does not necessarily imply larger Ω. In this paper, we numerically investigate the relationship between R₀ and Ω for a set of empirical temporal networks of human contacts. Among 31 explanatory descriptors of temporal network structure, we identify those that make R₀ an imperfect predictor of Ω. We find that descriptors related to both temporal and topological aspects affect the relationship between R₀ and Ω, but in different ways.

Predicting epidemic risk from past temporal contact data

Understanding how epidemics spread in a system is a crucial step to prevent and control outbreaks, with broad implications on the system’s functioning, health, and associated costs. This can be achieved by identifying the elements at higher risk of infection and implementing targeted surveillance and control measures. One important ingredient to consider is the pattern of disease-transmission contacts among the elements, however lack of data or delays in providing updated records may hinder its use, especially for time-varying patterns. Here we explore to what extent it is possible to use past temporal data of a system’s pattern of contacts to predict the risk of infection of its elements during an emerging outbreak, in absence of updated data. We focus on two real-world temporal systems; a livestock displacements trade network among animal holdings, and a network of sexual encounters in high-end prostitution. We define the node’s loyalty as a local measure of its tendency to maintain contacts with the same elements over time, and uncover important non-trivial correlations with the node’s epidemic risk. We show that a risk assessment analysis incorporating this knowledge and based on past structural and temporal pattern properties provides accurate predictions for both systems. Its generalizability is tested by introducing a theoretical model for generating synthetic temporal networks. High accuracy of our predictions is recovered across different settings, while the amount of possible predictions is system-specific. The proposed method can provide crucial information for the setup of targeted intervention strategies.

Following the emergence of a transmissible disease epidemic, interventions and resources need to be prioritized to efficiently control its spread. While the knowledge of the pattern of disease-transmission contacts among hosts would be ideal for this task, the continuously changing nature of such pattern makes its use less practical in real public health emergencies (or otherwise highly resource-demanding when possible). We show that in such situations critical knowledge to assess the real-time risk of infection can be extracted from past temporal contact data. An index expressing the conservation of contacts over time is proposed as an effective tool to prioritize interventions, and its efficiency is tested considering real data on livestock movements and on human sexual encounters.

Steady state and mean recurrence time for random walks on stochastic temporal networks

Random walks are basic diffusion processes on networks and have applications in, for example, searching, navigation, ranking, and community detection. Recent recognition of the importance of temporal aspects on networks spurred studies of random walks on temporal networks. Here we theoretically study two types of event-driven random walks on a stochastic temporal network model that produces arbitrary distributions of interevent times. In the so-called active random walk, the interevent time is reinitialized on all links upon each movement of the walker. In the so-called passive random walk, the interevent time is reinitialized only on the link that has been used the last time, and it is a type of correlated random walk. We find that the steady state is always the uniform density for the passive random walk. In contrast, for the active random walk, it increases or decreases with the node's degree depending on the distribution of interevent times. The mean recurrence time of a node is inversely proportional to the degree for both active and passive random walks. Furthermore, the mean recurrence time does or does not depend on the distribution of interevent times for the active and passive random walks, respectively.

Evolutionary dynamics of time-resolved social interactions

Cooperation among unrelated individuals is frequently observed in social groups when their members combine efforts and resources to obtain a shared benefit that is unachievable by an individual alone. However, understanding why cooperation arises despite the natural tendency of individuals toward selfish behavior is still an open problem and represents one of the most fascinating challenges in evolutionary dynamics. Recently, the structural characterization of the networks in which social interactions take place has shed some light on the mechanisms by which cooperative behavior emerges and eventually overcomes the natural temptation to defect. In particular, it has been found that the heterogeneity in the number of social ties and the presence of tightly knit communities lead to a significant increase in cooperation as compared with the unstructured and homogeneous connection patterns considered in classical evolutionary dynamics. Here, we investigate the role of social-ties dynamics for the emergence of cooperation in a family of social dilemmas. Social interactions are in fact intrinsically dynamic, fluctuating, and intermittent over time, and they can be represented by time-varying networks. By considering two experimental data sets of human interactions with detailed time information, we show that the temporal dynamics of social ties has a dramatic impact on the evolution of cooperation: the dynamics of pairwise interactions favors selfish behavior.

Causality-driven slow-down and speed-up of diffusion in non-Markovian temporal networks

Recent research has highlighted limitations of studying complex systems with time-varying topologies from the perspective of static, time-aggregated networks. Non-Markovian characteristics resulting from the ordering of interactions in temporal networks were identified as one important mechanism that alters causality and affects dynamical processes. So far, an analytical explanation for this phenomenon and for the significant variations observed across different systems is missing. Here we introduce a methodology that allows to analytically predict causality-driven changes of diffusion speed in non-Markovian temporal networks. Validating our predictions in six data sets we show that compared with the time-aggregated network, non-Markovian characteristics can lead to both a slow-down or speed-up of diffusion, which can even outweigh the decelerating effect of community structures in the static topology. Thus, non-Markovian properties of temporal networks constitute an important additional dimension of complexity in time-varying complex systems.

Temporal stability of network partitions

We present a method to find the best temporal partition at any time scale and rank the relevance of partitions found at different time scales. This method is based on random walkers coevolving with the network and as such constitutes a generalization of partition stability to the case of temporal networks. We show that, when applied to a toy model and real data sets, temporal stability uncovers structures that are persistent over meaningful time scales as well as important isolated events, making it an effective tool to study both abrupt changes and gradual evolution of a network mesoscopic structures.

Interplay of network dynamics and heterogeneity of ties on spreading dynamics

The structure of a network dramatically affects the spreading phenomena unfolding upon it. The contact distribution of the nodes has long been recognized as the key ingredient in influencing the outbreak events. However, limited knowledge is currently available on the role of the weight of the edges on the persistence of a pathogen. At the same time, recent works showed a strong influence of temporal network dynamics on disease spreading. In this work we provide an analytical understanding, corroborated by numerical simulations, about the conditions for infected stable state in weighted networks. In particular, we reveal the role of heterogeneity of edge weights and of the dynamic assignment of weights on the ties in the network in driving the spread of the epidemic. In this context we show that when weights are dynamically assigned to ties in the network, a heterogeneous distribution is able to hamper the diffusion of the disease, contrary to what happens when weights are fixed in time.

Birth and death of links control disease spreading in empirical contact networks

We investigate what structural aspects of a collection of twelve empirical temporal networks of human contacts are important to disease spreading. We scan the entire parameter spaces of the two canonical models of infectious disease epidemiology—the Susceptible-Infectious-Susceptible (SIS) and Susceptible-Infectious-Removed (SIR) models. The results from these simulations are compared to reference data where we eliminate structures in the interevent intervals, the time to the first contact in the data, or the time from the last contact to the end of the sampling. The picture we find is that the birth and death of links, and the total number of contacts over a link, are essential to predict outbreaks. On the other hand, the exact times of contacts between the beginning and end, or the interevent interval distribution, do not matter much. In other words, a simplified picture of these empirical data sets that suffices for epidemiological purposes is that links are born, is active with some intensity, and die.

Addressing population heterogeneity and distribution in epidemics models using a cellular automata approach

BACKGROUND

The spread of an infectious disease is determined by biological and social factors. Models based on cellular automata are adequate to describe such natural systems consisting of a massive collection of simple interacting objects. They characterize the time evolution of the global system as the emergent behaviour resulting from the interaction of the objects, whose behaviour is defined through a set of simple rules that encode the individual behaviour and the transmission dynamic.

METHODS

An epidemic is characterized trough an individual-based-model built upon cellular automata. In the proposed model, each individual of the population is represented by a cell of the automata. This way of modeling an epidemic situation allows to individually define the characteristic of each individual, establish different scenarios and implement control strategies.

RESULTS

A cellular automata model to study the time evolution of a heterogeneous populations through the various stages of disease was proposed, allowing the inclusion of individual heterogeneity, geographical characteristics and social factors that determine the dynamic of the desease. Different assumptions made to built the classical model were evaluated, leading to following results: i) for low contact rate (like in quarantine process or low density population areas) the number of infective individuals is lower than other areas where the contact rate is higher, and ii) for different initial spacial distributions of infected individuals different epidemic dynamics are obtained due to its influence on the transition rate and the reproductive ratio of disease.

CONCLUSIONS

The contact rate and spatial distributions have a central role in the spread of a disease. For low density populations the spread is very low and the number of infected individuals is lower than in highly populated areas. The spacial distribution of the population and the disease focus as well as the geographical characteristic of the area play a central role in the dynamics of the desease.

Toward an endgame: finding and engaging people unaware of their HIV-1 infection in treatment and prevention

Epidemic modeling suggests that a major scale-up in HIV treatment could have a dramatic impact on HIV incidence. This has led both researchers and policymakers to set a goal of an "AIDS-Free Generation." One of the greatest obstacles to achieving this objective is the number of people with undiagnosed HIV infection. Despite recent innovations, new research strategies are needed to identify, engage, and successfully treat people who are unaware of their infection.

Epidemiologically optimal static networks from temporal network data

One of network epidemiology's central assumptions is that the contact structure over which infectious diseases propagate can be represented as a static network. However, contacts are highly dynamic, changing at many time scales. In this paper, we investigate conceptually simple methods to construct static graphs for network epidemiology from temporal contact data. We evaluate these methods on empirical and synthetic model data. For almost all our cases, the network representation that captures most relevant information is a so-called exponential-threshold network. In these, each contact contributes with a weight decreasing exponentially with time, and there is an edge between a pair of vertices if the weight between them exceeds a threshold. Networks of aggregated contacts over an optimally chosen time window perform almost as good as the exponential-threshold networks. On the other hand, networks of accumulated contacts over the entire sampling time, and networks of concurrent partnerships, perform worse. We discuss these observations in the context of the temporal and topological structure of the data sets.

To understand how diseases spread in a population, it is important to study the network of people in contact. Many methods to model epidemic outbreaks make the assumption that one can treat this network as static. In reality, we know that contact patterns between people change in time, and old contacts are soon irrelevant—it does not matter that we know Marie Antoinette's lovers to understand the HIV epidemic. This paper investigates methods for constructing networks of people that are as relevant as possible for disease spreading. The most promising method we call exponential-threshold network works by letting contacts contribute less, the further from the beginning of an outbreak they take place. We investigate the methods both on artificial models of the contact patterns and empirical data. Except searching for the optimal network representation, we also investigate how the structure of the original data set affects the performance of the representations.

Flow motifs reveal limitations of the static framework to represent human interactions

Networks are commonly used to define underlying interaction structures where infections, information, or other quantities may spread. Although the standard approach has been to aggregate all links into a static structure, some studies have shown that the time order in which the links are established may alter the dynamics of spreading. In this paper, we study the impact of the time ordering in the limits of flow on various empirical temporal networks. By using a random walk dynamics, we estimate the flow on links and convert the original undirected network (temporal and static) into a directed flow network. We then introduce the concept of flow motifs and quantify the divergence in the representativity of motifs when using the temporal and static frameworks. We find that the regularity of contacts and persistence of vertices (common in email communication and face-to-face interactions) result on little differences in the limits of flow for both frameworks. On the other hand, in the case of communication within a dating site and of a sexual network, the flow between vertices changes significantly in the temporal framework such that the static approximation poorly represents the structure of contacts. We have also observed that cliques with 3 and 4 vertices containing only low-flow links are more represented than the same cliques with all high-flow links. The representativity of these low-flow cliques is higher in the temporal framework. Our results suggest that the flow between vertices connected in cliques depend on the topological context in which they are placed and in the time sequence in which the links are established. The structure of the clique alone does not completely characterize the potential of flow between the vertices.