Predicting and controlling infectious disease epidemics using temporal networks

An algorithm to build synthetic temporal contact networks based on close-proximity interactions data

Small populations (e.g., hospitals, schools or workplaces) are characterised by high contact heterogeneity and stochasticity affecting pathogen transmission dynamics. Empirical individual contact data provide unprecedented information to characterize such heterogeneity and are increasingly available, but are usually collected over a limited period, and can suffer from observation bias. We propose an algorithm to stochastically reconstruct realistic temporal networks from individual contact data in healthcare settings (HCS) and test this approach using real data previously collected in a long-term care facility (LTCF). Our algorithm generates full networks from recorded close-proximity interactions, using hourly inter-individual contact rates and information on individuals' wards, the categories of staff involved in contacts, and the frequency of recurring contacts. It also provides data augmentation by reconstructing contacts for days when some individuals are present in the HCS without having contacts recorded in the empirical data. Recording bias is formalized through an observation model, to allow direct comparison between the augmented and observed networks. We validate our algorithm using data collected during the i-Bird study, and compare the empirical and reconstructed networks. The algorithm was substantially more accurate to reproduce network characteristics than random graphs. The reconstructed networks reproduced well the assortativity by ward (first-third quartiles observed: 0.54-0.64; synthetic: 0.52-0.64) and the hourly staff and patient contact patterns. Importantly, the observed temporal correlation was also well reproduced (0.39-0.50 vs 0.37-0.44), indicating that our algorithm could recreate a realistic temporal structure. The algorithm consistently recreated unobserved contacts to generate full reconstructed networks for the LTCF. To conclude, we propose an approach to generate realistic temporal contact networks and reconstruct unobserved contacts from summary statistics computed using individual-level interaction networks. This could be applied and extended to generate contact networks to other HCS using limited empirical data, to subsequently inform individual-based epidemic models.

Analysing dynamic work systems using DynEAST: a demonstration of concept

The capability of current Ergonomics methods to capture dynamism is limited, stifling our understanding of work-as-done, distributed situational awareness and organisational drift. This paper provides a demonstration of concept of DynEAST; an extension of the EAST framework underpinned by principles from Dynamic Network Analysis, to capture elements of dynamism within work systems. The DynEAST concept is applied to a railway maintenance case study. Case study findings demonstrate how DynEAST outputs can be used to advance our understanding of the aforementioned phenomena and better equip practitioners for current and future Ergonomics challenges.Practitioner summary: This paper introduces the DynEAST method. DynEAST enables HF/E practitioners to model and analyse dynamic features of complex work systems. The development of DynEAST is timely due to the concurrent proliferation of increasingly complex sociotechnical systems and stagnation of HF/E methods development; particularly those able to model systemic dynamism. Abbreviations: DynEAST: dynamic event analysis of systemic teamwork; EAST: dynamic event analysis of systemic teamwork; HF/E: human factors and ergonomics; HF: human factors; DNA: dynamic network analysis; HTA: hierarchal task analysis; CWA: cognitive work analysis; CAST: causal analysis based on system theory; STAMP: system theoretic accident model and processes; FRAM: functional resonance analysis method; SNA: social network analysis; DSA: distributed situational awareness; PPO: possession protection officer; PO: protection officer; RTS: railway track signals; LPA: local possession authority; SMEs: subject matter experts.

Risk-aware temporal cascade reconstruction to detect asymptomatic cases

This paper studies the problem of detecting asymptomatic cases in a temporal contact network in which multiple outbreaks have occurred. We show that the key to detecting asymptomatic cases well is taking into account both individual risk and the likelihood of disease-flow along edges. We consider both aspects by formulating the asymptomatic case detection problem as a directed prize-collecting Steiner tree (Directed PCST) problem. We present an approximation-preserving reduction from this problem to the directed Steiner tree problem and obtain scalable algorithms for the Directed PCST problem on instances with more than 1.5M edges obtained from both synthetic and fine-grained hospital data. On synthetic data, we demonstrate that our detection methods significantly outperform various baselines (with a gain of 3.6 × ). We apply our method to the infectious disease prediction task by using an additional feature set that captures exposure to detected asymptomatic cases and show that our method outperforms all baselines. We further use our method to detect infection sources ("patient zero") of outbreaks that outperform baselines. We also demonstrate that the solutions returned by our approach are clinically meaningful by presenting case studies.

Hidden dependence of spreading vulnerability on topological complexity

Many dynamical phenomena in complex systems concern spreading that plays out on top of networks with changing architecture over time-commonly known as temporal networks. A complex system's proneness to facilitate spreading phenomena, which we abbreviate as its "spreading vulnerability," is often surmised to be related to the topology of the temporal network featured by the system. Yet, cleanly extracting spreading vulnerability of a complex system directly from the topological information of the temporal network remains a challenge. Here, using data from a diverse set of real-world complex systems, we develop the "entropy of temporal entanglement" as a quantity to measure topological complexities of temporal networks. We show that this parameter-free quantity naturally allows for topological comparisons across vastly different complex systems. Importantly, by simulating three different types of stochastic dynamical processes playing out on top of temporal networks, we demonstrate that the entropy of temporal entanglement serves as a quantitative embodiment of the systems' spreading vulnerability, irrespective of the details of the processes. In being able to do so, i.e., in being able to quantitatively extract a complex system's proneness to facilitate spreading phenomena from topology, this entropic measure opens itself for applications in a wide variety of natural, social, biological, and engineered systems.

Epidemic models characterize seizure propagation and the effects of epilepsy surgery in individualized brain networks based on MEG and invasive EEG recordings

Epilepsy surgery is the treatment of choice for drug-resistant epilepsy patients. However, seizure-freedom is currently achieved in only 2/3 of the patients after surgery. In this study we have developed an individualized computational model based on MEG brain networks to explore seizure propagation and the efficacy of different virtual resections. Eventually, the goal is to obtain individualized models to optimize resection strategy and outcome. We have modelled seizure propagation as an epidemic process using the susceptible-infected (SI) model on individual brain networks derived from presurgical MEG. We included 10 patients who had received epilepsy surgery and for whom the surgery outcome at least one year after surgery was known. The model parameters were tuned in in order to reproduce the patient-specific seizure propagation patterns as recorded with invasive EEG. We defined a personalized search algorithm that combined structural and dynamical information to find resections that maximally decreased seizure propagation for a given resection size. The optimal resection for each patient was defined as the smallest resection leading to at least a 90% reduction in seizure propagation. The individualized model reproduced the basic aspects of seizure propagation for 9 out of 10 patients when using the resection area as the origin of epidemic spreading, and for 10 out of 10 patients with an alternative definition of the seed region. We found that, for 7 patients, the optimal resection was smaller than the resection area, and for 4 patients we also found that a resection smaller than the resection area could lead to a 100% decrease in propagation. Moreover, for two cases these alternative resections included nodes outside the resection area. Epidemic spreading models fitted with patient specific data can capture the fundamental aspects of clinically observed seizure propagation, and can be used to test virtual resections in silico. Combined with optimization algorithms, smaller or alternative resection strategies, that are individually targeted for each patient, can be determined with the ultimate goal to improve surgery outcome. MEG-based networks can provide a good approximation of structural connectivity for computational models of seizure propagation, and facilitate their clinical use.

Quantifying agent impacts on contact sequences in social interactions

Human social behavior plays a crucial role in how pathogens like SARS-CoV-2 or fake news spread in a population. Social interactions determine the contact network among individuals, while spreading, requiring individual-to-individual transmission, takes place on top of the network. Studying the topological aspects of a contact network, therefore, not only has the potential of leading to valuable insights into how the behavior of individuals impacts spreading phenomena, but it may also open up possibilities for devising effective behavioral interventions. Because of the temporal nature of interactions—since the topology of the network, containing who is in contact with whom, when, for how long, and in which precise sequence, varies (rapidly) in time—analyzing them requires developing network methods and metrics that respect temporal variability, in contrast to those developed for static (i.e., time-invariant) networks. Here, by means of event mapping, we propose a method to quantify how quickly agents mingle by transforming temporal network data of agent contacts. We define a novel measure called contact sequence centrality, which quantifies the impact of an individual on the contact sequences, reflecting the individual’s behavioral potential for spreading. Comparing contact sequence centrality across agents allows for ranking the impact of agents and identifying potential ‘behavioral super-spreaders’. The method is applied to social interaction data collected at an art fair in Amsterdam. We relate the measure to the existing network metrics, both temporal and static, and find that (mostly at longer time scales) traditional metrics lose their resemblance to contact sequence centrality. Our work highlights the importance of accounting for the sequential nature of contacts when analyzing social interactions.

An integrated simulation framework for the prevention and mitigation of pandemics caused by airborne pathogens

In this work, we developed an integrated simulation framework for pandemic prevention and mitigation of pandemics caused by airborne pathogens, incorporating three sub-models, namely the spatial model, the mobility model, and the propagation model, to create a realistic simulation environment for the evaluation of the effectiveness of different countermeasures on the epidemic dynamics. The spatial model converts images of real cities obtained from Google Maps into undirected weighted graphs that capture the spatial arrangement of the streets utilized next for the mobility of individuals. The mobility model implements a stochastic agent-based approach, developed to assign specific routes to individuals moving in the city, through the use of stochastic processes, utilizing the weights of the underlying graph to deploy shortest path algorithms. The propagation model implements both the epidemiological model and the physical substance of the transmission of an airborne pathogen (in our approach, we investigate the transmission parameters of SARS-CoV-2). The deployment of a set of countermeasures was investigated in reducing the spread of the pathogen, where, through a series of repetitive simulation experiments, we evaluated the effectiveness of each countermeasure in pandemic prevention.

A pilot study using financial transactions’ spatial information to define high-risk neighborhoods and distribution pattern of COVID-19

Background

Development of surveillance systems based on big data sources with spatial information is necessitated more than ever during this pandemic. Here, we present our pilot results of a new technique for the incorporation of spatial information of transactions and a vital registry of COVID-19 to evaluate the disease spread.

Methods

We merged two databases of laboratory-confirmed national COVID-19 registry of Iran and financial transactions of point-of-sale devices from February to March 2020 as our training data sources. Spatial information was used for the visualization of maps and movements of sick individuals. We used the point-of-sale devices-related guild to check for the dynamics of financial transactions and effectiveness of quarantines.

Findings

In the study period, 174,428 confirmed cases were in the COVID-19 registry with accompanying transactions information. In total, 13,924,982 financial transactions were performed by them, with a mean of 1.2 per day for each person. All guilds had a decreasing pattern of “risky” transactions except for grocery stores and pharmacies. The latter showed a decreasing pattern by impose of lockdowns. Different cities were the hotspot of disease transmission as many “high-risk” transactions were performed in them, among which Tehran (mainly its central neighborhoods) and southern cities of Lake Urmia predominated. Lockdowns indicated that the disease gradually became less transmissible.

Interpretation

Financial transactions can be readily used for epidemics surveillance. Semi real-time results of such iterations can be informative for policy makers, guild owners, and general population to prepare safer commuting and merchandise spaces.

Revealing mechanisms of infectious disease spread through empirical contact networks

The spread of pathogens fundamentally depends on the underlying contacts between individuals. Modeling the dynamics of infectious disease spread through contact networks, however, can be challenging due to limited knowledge of how an infectious disease spreads and its transmission rate. We developed a novel statistical tool, INoDS (Identifying contact Networks of infectious Disease Spread) that estimates the transmission rate of an infectious disease outbreak, establishes epidemiological relevance of a contact network in explaining the observed pattern of infectious disease spread and enables model comparison between different contact network hypotheses. We show that our tool is robust to incomplete data and can be easily applied to datasets where infection timings of individuals are unknown. We tested the reliability of INoDS using simulation experiments of disease spread on a synthetic contact network and find that it is robust to incomplete data and is reliable under different settings of network dynamics and disease contagiousness compared with previous approaches. We demonstrate the applicability of our method in two host-pathogen systems: Crithidia bombi in bumblebee colonies and Salmonella in wild Australian sleepy lizard populations. INoDS thus provides a novel and reliable statistical tool for identifying transmission pathways of infectious disease spread. In addition, application of INoDS extends to understanding the spread of novel or emerging infectious disease, an alternative approach to laboratory transmission experiments, and overcoming common data-collection constraints.

Network models are widely used to understand and predict infectious disease spread in human and animal populations. However, the choice of network model often relies on subjective expert knowledge or disease transmission experiments that are time-consuming and difficult to perform. We developed a novel tool, called INoDS (Identifying contact Networks of infectious Disease Spread), that uses robust statistical approach to establish relevance of a network model in explaining transmission pathways of an infectious disease outbreak. We used computer simulations and real-world dataset to test the accuracy of our tool and robustness to missing data. We found that INoDS is robust to common data-collection constraints, broadly applicable and accurate compared to current approaches. The tool that we have developed can therefore provide immediate epidemiological insights in the event of an epidemic outbreak, and can be used to improve targeted disease control.

Epidemic dynamics on metapopulation networks with node2vec mobility

Metapopulation models have been a powerful tool for both theorizing and simulating epidemic dynamics. In a metapopulation model, one considers a network composed of subpopulations and their pairwise connections, and individuals are assumed to migrate from one subpopulation to another obeying a given mobility rule. While how different mobility rules affect epidemic dynamics in metapopulation models has been studied, there have been relatively few efforts on comparison of the effects of simple (i.e., unbiased) random walks and more complex mobility rules. Here we study a susceptible-infectious-susceptible (SIS) dynamics in a metapopulation model in which individuals obey a parametric second-order random-walk mobility rule called the node2vec. We map the second-order mobility rule of the node2vec to a first-order random walk in a network whose each node is a directed edge connecting a pair of subpopulations and then derive the epidemic threshold. For various networks, we find that the epidemic threshold is large (therefore, epidemic spreading tends to be suppressed) when the individuals infrequently backtrack or infrequently visit the common neighbors of the currently visited and the last visited subpopulations than when the individuals obey the simple random walk. The amount of change in the epidemic threshold induced by the node2vec mobility is in general not as large as, but is sometimes comparable with, the one induced by the change in the diffusion rate for individuals.

Analyzing spatial mobility patterns with time-varying graphical lasso: Application to COVID-19 spread

This work applies the time-varying graphical lasso (TVGL) method, an extension of the traditional graphical lasso approach, to address learning time-varying graphs from spatiotemporal measurements. Given georeferenced data, the TVGL method can estimate a time-varying network where an edge represents a partial correlation between two nodes. To achieve this, we use a COVID-19 data set from the Argentine province of Chaco. As an application, we use the estimated network to study the impact of COVID-19 confinement measures and evaluate whether the measures produced the expected result.

Concurrency measures in the era of temporal network epidemiology: a review

Diseases spread over temporal networks of interaction events between individuals. Structures of these temporal networks hold the keys to understanding epidemic propagation. One early concept of the literature to aid in discussing these structures is concurrency-quantifying individuals' tendency to form time-overlapping 'partnerships'. Although conflicting evaluations and an overabundance of operational definitions have marred the history of concurrency, it remains important, especially in the area of sexually transmitted infections. Today, much of theoretical epidemiology uses more direct models of contact patterns, and there is an emerging body of literature trying to connect methods to the concurrency literature. In this review, we will cover the development of the concept of concurrency and these new approaches.

Fast and principled simulations of the SIR model on temporal networks

The Susceptible-Infectious-Recovered (SIR) model is the canonical model of epidemics of infections that make people immune upon recovery. Many of the open questions in computational epidemiology concern the underlying contact structure's impact on models like the SIR model. Temporal networks constitute a theoretical framework capable of encoding structures both in the networks of who could infect whom and when these contacts happen. In this article, we discuss the detailed assumptions behind such simulations-how to make them comparable with analytically tractable formulations of the SIR model, and at the same time, as realistic as possible. We also present a highly optimized, open-source code for this purpose and discuss all steps needed to make the program as fast as possible.

The Effectiveness of Interactive Visualization Techniques for Time Navigation of Dynamic Graphs on Large Displays

Dynamic networks can be challenging to analyze visually, especially if they span a large time range during which new nodes and edges can appear and disappear. Although it is straightforward to provide interfaces for visualization that represent multiple states of the network (i.e., multiple timeslices) either simultaneously (e.g., through small multiples) or interactively (e.g., through interactive animation), these interfaces might not support tasks in which disjoint timeslices need to be compared. Since these tasks are key for understanding the dynamic aspects of the network, understanding which interactive visualizations best support these tasks is important. We present the results of a series of laboratory experiments comparing two traditional approaches (small multiples and interactive animation), with a more recent approach based on interactive timeslicing. The tasks were performed on a large display through a touch interface. Participants completed 24 trials of three tasks with all techniques. The results show that interactive timeslicing brings benefit when comparing distant points in time, but less benefits when analyzing contiguous intervals of time.

Epidemic spreading in a weighted pig trade network

The inclusion of edge weights can add valuable insights in the spreading processes within trade networks and may identify factors influencing the final epidemic size. The aim of the study was to evaluate the effect of different network versions on the outcome of an epidemiological model. The weighted network versions included the number of trade contacts (A), the sum of delivered livestock (B) and the mean number of delivered livestock per trade contact (C). Furthermore, other factors, e.g. transmission probability and farm type of primary outbreak, were tested for their impact on the final epidemic size. From 2013-2014, data from a pig trade network in Northern Germany was recorded containing 678 farms connected by 1,018 directed edges. An epidemiological model was implemented considering a higher probability of disease spread for edges with a higher weight for each of the combinations between network version and transmission probability. Only transmission routes following the network structure were considered for disease transmission. The outcome of the epidemiological model (number of infected farms) was tested with a generalized linear mixed model including the fixed effects network version (unweighted, A, B, C), transmission probability and farm type of primary outbreak (breeding farm, farrowing farm, finishing farm, farrow-to-finishing farm, unknown) as well as all twofold interactions. The results revealed that all fixed effects as well as all twofold interactions were significant (p ≤ 0.05), i.e. in the following only the impact of the interactions on the number of infected farms can be interpreted. Network versions B and C showed in all combinations the highest number of infected farms independent of the underlying transmission probability. The unweighted network and network version A showed a significant increase of infected farms with increasing transmission probability. All interactions including the farm type of primary outbreak revealed a significant higher number of infected farms for farm types located at the beginning of the production chain, e.g. breeding farms. These farm types reached also more other farms in 1-4 steps compared to farm types located near to the end of the production chain. The inclusion of edge weights has a significant effect on the outcome of epidemiological models and dependent on the chosen edge weight the results need to be interpreted accordingly.

Stay with your community: Bridges between clusters trigger expansion of COVID-19

In this study, the spread of virus infection was simulated using artificial human networks. Here, real-space urban life was modeled as a modified scale-free network with constraints. To date, the scale-free network has been adopted for modeling online communities in several studies. However, in the present study, it has been modified to represent the social behaviors of people where the generated communities are restricted and reflect spatiotemporal constraints in real life. Furthermore, the networks have been extended by introducing multiple cliques in the initial step of network construction and enabling people to contact hidden (zero-degree) as well as popular (large-degree) people. Consequently, four findings and a policy proposal were obtained. First, "second waves" were observed in some cases of the simulations even without external influence or constraints on people's social contacts or the releasing of the constraints. These waves tend to be lower than the first wave and occur in "fresh" clusters, that is, via the infection of people who are connected in the network but have not been infected previously. This implies that the bridge between infected and fresh clusters may trigger a new spread of the virus. Second, if the network changes its structure on the way of infection spread or after its suppression, a second wave larger than the first can occur. Third, the peak height in the time series of the number of infected cases depends on the difference between the upper bound of the number of people each member actually meets and the number of people they choose to meet during the period of infection spread. This tendency is observed for the two kinds of artificial networks introduced here and implies the impact of bridges between communities on the virus spreading. Fourth, the release of a previously imposed constraint may trigger a second wave higher than the peak of the time series without introducing any constraint so far previously, if the release is introduced at a time close to the peak. Thus, overall, both the government and individuals should be careful in returning to society where people enjoy free inter-community contact.

Generative models of simultaneously heavy-tailed distributions of interevent times on nodes and edges

Intervals between discrete events representing human activities, as well as other types of events, often obey heavy-tailed distributions, and their impacts on collective dynamics on networks such as contagion processes have been intensively studied. The literature supports that such heavy-tailed distributions are present for interevent times associated with both individual nodes and individual edges in networks. However, the simultaneous presence of heavy-tailed distributions of interevent times for nodes and edges is a nontrivial phenomenon, and its origin has been elusive. In the present study, we propose a generative model and its variants to explain this phenomenon. We assume that each node independently transits between a high-activity and low-activity state according to a continuous-time two-state Markov process and that, for the main model, events on an edge occur at a high rate if and only if both end nodes of the edge are in the high-activity state. In other words, two nodes interact frequently only when both nodes prefer to interact with others. The model produces distributions of interevent times for both individual nodes and edges that resemble heavy-tailed distributions across some scales. It also produces positive correlation in consecutive interevent times, which is another stylized observation for empirical data of human activity. We expect that our modeling framework provides a useful benchmark for investigating dynamics on temporal networks driven by non-Poissonian event sequences.

Fast Response to Superspreading: Uncertainty and Complexity in the Context of COVID-19

Although the first coronavirus disease 2019 (COVID-19) wave has peaked with the second wave underway, the world is still struggling to manage potential systemic risks and unpredictability of the pandemic. A particular challenge is the "superspreading" of the virus, which starts abruptly, is difficult to predict, and can quickly escalate into medical and socio-economic emergencies that contribute to long-lasting crises challenging our current ways of life. In these uncertain times, organizations and societies worldwide are faced with the need to develop appropriate strategies and intervention portfolios that require fast understanding of the complex interdependencies in our world and rapid, flexible action to contain the spread of the virus as quickly as possible, thus preventing further disastrous consequences of the pandemic. We integrate perspectives from systems sciences, epidemiology, biology, social networks, and organizational research in the context of the superspreading phenomenon to understand the complex system of COVID-19 pandemic and develop suggestions for interventions aimed at rapid responses.

Topological epidemic model: Theoretical insight into underlying networks

Although there are various models of epidemic diseases, there are a few individual-based models that can guide susceptible individuals on how they should behave in a pandemic without its appropriate treatment. Such a model would be ideal for the current coronavirus disease 2019 (COVID-19) pandemic. Thus, here, we propose a topological model of an epidemic disease, which can take into account various types of interventions through a time-dependent contact network. Based on this model, we show that there is a maximum allowed number of persons one can see each day for each person so that we can suppress the epidemic spread. Reducing the number of persons to see for the hub persons is a key countermeasure for the current COVID-19 pandemic.

Active and inactive quarantine in epidemic spreading on adaptive activity-driven networks

We consider an epidemic process on adaptive activity-driven temporal networks, with adaptive behavior modeled as a change in activity and attractiveness due to infection. By using a mean-field approach, we derive an analytical estimate of the epidemic threshold for susceptible-infected-susceptible (SIS) and susceptible-infected-recovered (SIR) epidemic models for a general adaptive strategy, which strongly depends on the correlations between activity and attractiveness in the susceptible and infected states. We focus on strong social distancing, implementing two types of quarantine inspired by recent real case studies: an active quarantine, in which the population compensates the loss of links rewiring the ineffective connections towards nonquarantining nodes, and an inactive quarantine, in which the links with quarantined nodes are not rewired. Both strategies feature the same epidemic threshold but they strongly differ in the dynamics of the active phase. We show that the active quarantine is extremely less effective in reducing the impact of the epidemic in the active phase compared to the inactive one and that in the SIR model a late adoption of measures requires inactive quarantine to reach containment.

Relevance of temporal cores for epidemic spread in temporal networks

Temporal networks are widely used to represent a vast diversity of systems, including in particular social interactions, and the spreading processes unfolding on top of them. The identification of structures playing important roles in such processes remains largely an open question, despite recent progresses in the case of static networks. Here, we consider as candidate structures the recently introduced concept of span-cores: the span-cores decompose a temporal network into subgraphs of controlled duration and increasing connectivity, generalizing the core-decomposition of static graphs. To assess the relevance of such structures, we explore the effectiveness of strategies aimed either at containing or maximizing the impact of a spread, based respectively on removing span-cores of high cohesiveness or duration to decrease the epidemic risk, or on seeding the process from such structures. The effectiveness of such strategies is assessed in a variety of empirical data sets and compared to baselines that use only static information on the centrality of nodes and static concepts of coreness, as well as to a baseline based on a temporal centrality measure. Our results show that the most stable and cohesive temporal cores play indeed an important role in epidemic processes on temporal networks, and that their nodes are likely to include influential spreaders.

Infectious Disease Research in the Era of Big Data

Infectious disease research spans scales from the molecular to the global—from specific mechanisms of pathogen drug resistance, virulence, and replication to the movement of people, animals, and pathogens around the world. All of these research areas have been impacted by the recent growth of large-scale data sources and data analytics. Some of these advances rely on data or analytic methods that are common to most biomedical data science, while others leverage the unique nature of infectious disease, namely its communicability. This review outlines major research progress in the past few years and highlights some remaining opportunities, focusing on data or methodological approaches particular to infectious disease.

Illustration of Different Disease Transmission Routes in a Pig Trade Network by Monopartite and Bipartite Representation

Simple Summary

Besides direct animal movements between farms; indirect transmission routes of pathogens can have an immense impact on network structure and disease spread in animal trade networks. This study integrated these indirect transmission routes between farms via transport companies or feed supply as bipartite networks; which were compared to the monopartite animal movements network representing the direct transmission route. Both bipartite networks were projected on farm level to enable a comparison to the monopartite network. The number of edges increased immensely from the monopartite animal movements network to both projected networks. Thus, farms can be highly connected over indirect connections, although they are not directly trading animals. The ranking of the animals according to their centrality parameters, indicating their importance for the network, showed moderate correlations only between the animal movements and the transportation network. The epidemiological models based on the different network representations revealed significantly more infected farms for the networks including indirect transmission routes compared to the direct animal movements. Indirect transmission routes had an immense impact on the outcome of centrality parameters, as well as on the spreading process within the network. This knowledge is needed to understand disease spread and to establish reliable prevention and control measurements.

Abstract

Besides the direct transport of animals, also indirect transmission routes, e.g., contact via contaminated vehicles, have to be considered. In this study, the transmission routes of a German pig trade network were illustrated as a monopartite animal movements network and two bipartite networks including information of the transport company and the feed producer which were projected on farm level (n = 866) to enable a comparison. The networks were investigated with the help of network analysis and formed the basis for epidemiological models to evaluate the impact of different transmission routes on network structure as well as on potential epidemic sizes. The number of edges increased immensely from the monopartite animal movements network to both projected networks. The median centrality parameters revealed clear differences between the three representations. Furthermore, moderate correlation coefficients ranging from 0.55 to 0.68 between the centrality values of the animal movements network and the projected transportation network were obtained. The epidemiological models revealed significantly more infected farms for both projected networks (70% to 100%) compared to the animal movements network (1%). The inclusion of indirect transmission routes had an immense impact on the outcome of centrality parameters as well as on the results of the epidemiological models.

Dynamic contact networks of patients and MRSA spread in hospitals

Methicillin-resistant Staphylococcus aureus (MRSA) is a difficult-to-treat infection. Increasing efforts have been taken to mitigate the epidemics and to avoid potential outbreaks in low endemic settings. Understanding the population dynamics of MRSA is essential to identify the causal mechanisms driving the epidemics and to generalise conclusions to different contexts. Previous studies neglected the temporal structure of contacts between patients and assumed homogeneous behaviour. We developed a high-resolution data-driven contact network model of interactions between 743,182 patients in 485 hospitals during 3,059 days to reproduce the exact contact sequences of the hospital population. Our model captures the exact spatial and temporal human contact behaviour and the dynamics of referrals within and between wards and hospitals at a large scale, revealing highly heterogeneous contact and mobility patterns of individual patients. A simulation exercise of epidemic spread shows that heterogeneous contacts cause the emergence of super-spreader patients, slower than exponential polynomial growth of the prevalence, and fast epidemic spread between wards and hospitals. In our simulated scenarios, screening upon hospital admittance is potentially more effective than reducing infection probability to reduce the final outbreak size. Our findings are useful to understand not only MRSA spread but also other hospital-acquired infections.

Containing misinformation spreading in temporal social networks

Many researchers from a variety of fields, including computer science, network science, and mathematics, have focused on how to contain the outbreaks of Internet misinformation that threaten social systems and undermine societal health. Most research on this topic treats the connections among individuals as static, but these connections change in time, and thus social networks are also temporal networks. Currently, there is no theoretical approach to the problem of containing misinformation outbreaks in temporal networks. We thus propose a misinformation spreading model for temporal networks and describe it using a new theoretical approach. We propose a heuristic-containing (HC) strategy based on optimizing the final outbreak size that outperforms simplified strategies such as those that are random-containing and targeted-containing. We verify the effectiveness of our HC strategy on both artificial and real-world networks by performing extensive numerical simulations and theoretical analyses. We find that the HC strategy dramatically increases the outbreak threshold and decreases the final outbreak threshold.

Fitness preferential attachment as a driving mechanism in bitcoin transaction network

Bitcoin is the earliest cryptocurrency and among the most successful ones to date. Recently, its dynamical evolution has attracted the attention of the research community due to its completeness and richness in historical records. In this paper, we focus on the detailed evolution of bitcoin trading with the aim of elucidating the mechanism that drives the formation of the bitcoin transaction network. Our empirical investigation reveals that although the temporal properties of the transaction network possesses scale-free degree distribution like many other networks, its formation mechanism is different from the commonly assumed models of degree preferential attachment or wealth preferential attachment. By defining the fitness value of each node as the ability of the node to attract new connections, we have instead uncovered that the observed scale-free degree distribution results from the intrinsic fitness of each node following a power-law distribution. Our finding thus suggests that the "good-get-richer" rather than the "rich-get-richer" paradigm operates within the bitcoin ecosystem. Based on these findings, we propose a model that captures the temporal generative process by means of a fitness preferential attachment and data-driven birth/death mechanism. Our proposed model is able to produce structural properties in good agreement with those obtained from the empirical bitcoin network.

Spreading dynamics in a cattle trade network: Size, speed, typical profile and consequences on epidemic control strategies

Infections can spread among livestock notably because infected animals can be brought to uncontaminated holdings, therefore exposing a new group of susceptible animals to the disease. As a consequence, the structure and dynamics of animal trade networks is a major focus of interest to control zoonosis. We investigate the impact of the chronology of animal trades on the dynamics of the process. Precisely, in the context of a basic SI model spreading, we measure on the French database of bovine transfers to what extent a snapshot-based analysis of the cattle trade networks overestimates the epidemic risks. We bring into light that an analysis taking into account the chronology of interactions would give a much more accurate assessment of both the size and speed of the process. For this purpose, we model data as a temporal network that we analyze using the link stream formalism in order to mix structural and temporal aspects. We also show that in this dataset, a basic SI spreading comes down in most cases to a simple two-phases scenario: a waiting period, with few contacts and low activity, followed by a linear growth of the number of infected holdings. Using this portrait of the spreading process, we identify efficient strategies to control a potential outbreak, based on the identification of specific elements of the link stream which have a higher probability to be involved in a spreading process.

The pig transport network in Switzerland: Structure, patterns, and implications for the transmission of infectious diseases between animal holdings

The topology of animal transport networks contributes substantially to how fast and to what extent a disease can transmit between animal holdings. Therefore, public authorities in many countries mandate livestock holdings to report all movements of animals. However, the reported data often does not contain information about the exact sequence of transports, making it impossible to assess the effect of truck sharing and truck contamination on disease transmission. The aim of this study was to analyze the topology of the Swiss pig transport network by means of social network analysis and to assess the implications for disease transmission between animal holdings. In particular, we studied how additional information about transport sequences changes the topology of the contact network. The study is based on the official animal movement database in Switzerland and a sample of transport data from one transport company. The results show that the Swiss pig transport network is highly fragmented, which mitigates the risk of a large-scale disease outbreak. By considering the time sequence of transports, we found that even in the worst case, only 0.34% of all farm-pairs were connected within one month. However, both network connectivity and individual connectedness of farms increased if truck sharing and especially truck contamination were considered. Therefore, the extent to which a disease may be transmitted between animal holdings may be underestimated if we only consider data from the official animal movement database. Our results highlight the need for a comprehensive analysis of contacts between farms that includes indirect contacts due to truck sharing and contamination. As the nature of animal transport networks is inherently temporal, we strongly suggest the use of temporal network measures in order to evaluate individual and overall risk of disease transmission through animal transportation.

Particle velocity controls phase transitions in contagion dynamics

Interactions often require the proximity between particles. The movement of particles, thus, drives the change of the neighbors which are located in their proximity, leading to a sequence of interactions. In pathogenic contagion, infections occur through proximal interactions, but at the same time, the movement facilitates the co-location of different strains. We analyze how the particle velocity impacts on the phase transitions on the contagion process of both a single infection and two cooperative infections. First, we identify an optimal velocity (close to half of the interaction range normalized by the recovery time) associated with the largest epidemic threshold, such that decreasing the velocity below the optimal value leads to larger outbreaks. Second, in the cooperative case, the system displays a continuous transition for low velocities, which becomes discontinuous for velocities of the order of three times the optimal velocity. Finally, we describe these characteristic regimes and explain the mechanisms driving the dynamics.

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

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

Objective measures for sentinel surveillance in network epidemiology

Assume one has the capability of determining whether a node in a network is infectious or not by probing it. Then problem of optimizing sentinel surveillance in networks is to identify the nodes to probe such that an emerging disease outbreak can be discovered early or reliably. Whether the emphasis should be on early or reliable detection depends on the scenario in question. We investigate three objective measures from the literature quantifying the performance of nodes in sentinel surveillance: the time to detection or extinction, the time to detection, and the frequency of detection. As a basis for the comparison, we use the susceptible-infectious-recovered model on static and temporal networks of human contacts. We show that, for some regions of parameter space, the three objective measures can rank the nodes very differently. This means sentinel surveillance is a class of problems, and solutions need to chose an objective measure for the particular scenario in question. As opposed to other problems in network epidemiology, we draw similar conclusions from the static and temporal networks. Furthermore, we do not find one type of network structure that predicts the objective measures, i.e., that depends both on the data set and the SIR parameter values.

Speeding up non-Markovian first-passage percolation with a few extra edges

One model of real-life spreading processes is that of first-passage percolation (also called the SI model) on random graphs. Social interactions often follow bursty patterns, which are usually modelled with independent and identically distributed heavy-tailed passage times on edges. On the other hand, random graphs are often locally tree-like, and spreading on trees with leaves might be very slow due to bottleneck edges with huge passage times. Here we consider the SI model with passage times following a power-law distribution ℙ(ξ>t)∼t-α with infinite mean. For any finite connected graph G with a root s, we find the largest number of vertices κ(G,s) that are infected in finite expected time, and prove that for every k≤κ(G,s), the expected time to infect k vertices is at most O(k1/α). Then we show that adding a single edge from s to a random vertex in a random tree 𝒯 typically increases κ(𝒯,s) from a bounded variable to a fraction of the size of 𝒯, thus severely accelerating the process. We examine this acceleration effect on some natural models of random graphs: critical Galton--Watson trees conditioned to be large, uniform spanning trees of the complete graph, and on the largest cluster of near-critical Erdős‒Rényi graphs. In particular, at the upper end of the critical window, the process is already much faster than exactly at criticality.

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.

The reachability of contagion in temporal contact networks: how disease latency can exploit the rhythm of human behavior

Background

The symptoms of many infectious diseases influence their host to withdraw from social activity limiting their potential to spread. Successful transmission therefore requires the onset of infectiousness to coincide with a time when the host is socially active. Since social activity and infectiousness are both temporal phenomena, we hypothesize that diseases are most pervasive when these two processes are synchronized.

Methods

We consider disease dynamics that incorporate behavioral responses that effectively shorten the infectious period of the pathogen. Using data collected from face-to-face social interactions and synthetic contact networks constructed from empirical demographic data, we measure the reachability of this disease model and perform disease simulations over a range of latent period durations.

Results

We find that maximum transmission risk results when the disease latent period (and thus the generation time) are synchronized with human circadian rhythms of 24 h, and minimum transmission risk when latent periods are out of phase with circadian rhythms by 12 h. The effect of this synchronization is present for a range of disease models with realistic disease parameters and host behavioral responses.

Conclusions

The reproductive potential of pathogens is linked inextricably to the host social behavior required for transmission. We propose that future work should consider contact periodicity in models of disease dynamics, and suggest the possibility that disease control strategies may be designed to optimize against the effects of synchronization.

Electronic supplementary material

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

Natural Disasters and Cholera Outbreaks: Current Understanding and Future Outlook

PURPOSE OF REVIEW

Diarrheal diseases remain a serious global public health threat, especially for those populations lacking access to safe water and sanitation infrastructure. Although association of several diarrheal diseases, e.g., cholera, shigellosis, etc., with climatic processes has been documented, the global human population remains at heightened risk of outbreak of diseases after natural disasters, such as earthquakes, floods, or droughts. In this review, cholera was selected as a signature diarrheal disease and the role of natural disasters in triggering and transmitting cholera was analyzed.

RECENT FINDINGS

Key observations include identification of an inherent feedback loop that includes societal structure, prevailing climatic processes, and spatio-temporal seasonal variability of natural disasters. Data obtained from satellite-based remote sensing are concluded to have application, although limited, in predicting risks of a cholera outbreak(s). We argue that with the advent of new high spectral and spatial resolution data, earth observation systems should be seamlessly integrated in a decision support mechanism to be mobilize resources when a region suffers a natural disaster. A framework is proposed that can be used to assess the impact of natural disasters with response to outbreak of cholera, providing assessment of short- and long-term influence of climatic processes on disease outbreaks.

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.

Epidemic spreading in modular time-varying networks

We investigate the effects of modular and temporal connectivity patterns on epidemic spreading. To this end, we introduce and analytically characterise a model of time-varying networks with tunable modularity. Within this framework, we study the epidemic size of Susceptible-Infected-Recovered, SIR, models and the epidemic threshold of Susceptible-Infected-Susceptible, SIS, models. Interestingly, we find that while the presence of tightly connected clusters inhibits SIR processes, it speeds up SIS phenomena. In this case, we observe that modular structures induce a reduction of the threshold with respect to time-varying networks without communities. We confirm the theoretical results by means of extensive numerical simulations both on synthetic graphs as well as on a real modular and temporal network.

Sampling of temporal networks: Methods and biases

Temporal networks have been increasingly used to model a diversity of systems that evolve in time; for example, human contact structures over which dynamic processes such as epidemics take place. A fundamental aspect of real-life networks is that they are sampled within temporal and spatial frames. Furthermore, one might wish to subsample networks to reduce their size for better visualization or to perform computationally intensive simulations. The sampling method may affect the network structure and thus caution is necessary to generalize results based on samples. In this paper, we study four sampling strategies applied to a variety of real-life temporal networks. We quantify the biases generated by each sampling strategy on a number of relevant statistics such as link activity, temporal paths and epidemic spread. We find that some biases are common in a variety of networks and statistics, but one strategy, uniform sampling of nodes, shows improved performance in most scenarios. Given the particularities of temporal network data and the variety of network structures, we recommend that the choice of sampling methods be problem oriented to minimize the potential biases for the specific research questions on hand. Our results help researchers to better design network data collection protocols and to understand the limitations of sampled temporal network data.

Cost-efficient vaccination protocols for network epidemiology

We investigate methods to vaccinate contact networks—i.e. removing nodes in such a way that disease spreading is hindered as much as possible—with respect to their cost-efficiency. Any real implementation of such protocols would come with costs related both to the vaccination itself, and gathering of information about the network. Disregarding this, we argue, would lead to erroneous evaluation of vaccination protocols. We use the susceptible-infected-recovered model—the generic model for diseases making patients immune upon recovery—as our disease-spreading scenario, and analyze outbreaks on both empirical and model networks. For different relative costs, different protocols dominate. For high vaccination costs and low costs of gathering information, the so-called acquaintance vaccination is the most cost efficient. For other parameter values, protocols designed for query-efficient identification of the network’s largest degrees are most efficient.

Finding methods to identify important spreaders—and consequently protocols to identify individuals to vaccinate in targeted vaccination campaigns—is one of the most important topics of network theory. Earlier studies typically make some assumption about what information is available about the contact network that the disease spreads over. Then they try to optimize an objective function—either the average outbreak size in disease simulations, or (simpler) the size of the largest connected component. For public-health practitioners, gathering the network information cannot be detached from the decision process—their cost function includes the costs for both the vaccination itself and mapping of the network. This is the first paper to evaluate the cost efficiency of vaccination protocols—a problem that is much more relevant and not so much more complicated, than the oversimplified objective functions optimized in previous studies. We find a “no-free lunch” situation, where different protocols proposed in the past are most efficient at different cost scenarios. However, some methods are never cost efficient due to the amount of information they need. What protocol that is the best depends on network structure in a non-trivial way. We use both analytical and simulation techniques to reach these conclusions.

Concurrency-Induced Transitions in Epidemic Dynamics on Temporal Networks

Social contact networks underlying epidemic processes in humans and animals are highly dynamic. The spreading of infections on such temporal networks can differ dramatically from spreading on static networks. We theoretically investigate the effects of concurrency, the number of neighbors that a node has at a given time point, on the epidemic threshold in the stochastic susceptible-infected-susceptible dynamics on temporal network models. We show that network dynamics can suppress epidemics (i.e., yield a higher epidemic threshold) when the node's concurrency is low, but can also enhance epidemics when the concurrency is high. We analytically determine different phases of this concurrency-induced transition, and confirm our results with numerical simulations.

Interplay between epidemic spread and information propagation on metapopulation networks

The spread of an infectious disease has been widely found to evolve with the propagation of information. Many seminal works have demonstrated the impact of information propagation on the epidemic spreading, assuming that individuals are static and no mobility is involved. Inspired by the recent observation of diverse mobility patterns, we incorporate the information propagation into a metapopulation model based on the mobility patterns and contagion process, which significantly alters the epidemic threshold. In more details, we find that both the information efficiency and the mobility patterns have essential impacts on the epidemic spread. We obtain different scenarios leading to the mitigation of the outbreak by appropriately integrating the mobility patterns and the information efficiency as well. The inclusion of the impacts of the information propagation into the epidemiological model is expected to provide an support to public health implications for the suppression of epidemics.

SCOUT: simultaneous time segmentation and community detection in dynamic networks

Many evolving complex real-world systems can be modeled via dynamic networks. An important problem in dynamic network research is community detection, which finds groups of topologically related nodes. Typically, this problem is approached by assuming either that each time point has a distinct community organization or that all time points share a single community organization. The reality likely lies between these two extremes. To find the compromise, we consider community detection in the context of the problem of segment detection, which identifies contiguous time periods with consistent network structure. Consequently, we formulate a combined problem of segment community detection (SCD), which simultaneously partitions the network into contiguous time segments with consistent community organization and finds this community organization for each segment. To solve SCD, we introduce SCOUT, an optimization framework that explicitly considers both segmentation quality and partition quality. SCOUT addresses limitations of existing methods that can be adapted to solve SCD, which consider only one of segmentation quality or partition quality. In a thorough evaluation, SCOUT outperforms the existing methods in terms of both accuracy and computational complexity. We apply SCOUT to biological network data to study human aging.

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.

Data-driven network modelling of disease transmission using complete population movement data: spread of VTEC O157 in Swedish cattle

European Union legislation requires member states to keep national databases of all bovine animals. This allows for disease spread models that includes the time-varying contact network and population demographic. However, performing data-driven simulations with a high degree of detail are computationally challenging. We have developed an efficient and flexible discrete-event simulator SimInf for stochastic disease spread modelling that divides work among multiple processors to accelerate the computations. The model integrates disease dynamics as continuous-time Markov chains and livestock data as events. In this study, all Swedish livestock data (births, movements and slaughter) from July 1st 2005 to December 31st 2013 were included in the simulations. Verotoxigenic Escherichia coli O157:H7 (VTEC O157) are capable of causing serious illness in humans. Cattle are considered to be the main reservoir of the bacteria. A better understanding of the epidemiology in the cattle population is necessary to be able to design and deploy targeted measures to reduce the VTEC O157 prevalence and, subsequently, human exposure. To explore the spread of VTEC O157 in the entire Swedish cattle population during the period under study, a within- and between-herd disease spread model was used. Real livestock data was incorporated to model demographics of the population. Cattle were moved between herds according to real movement data. The results showed that the spatial pattern in prevalence may be due to regional differences in livestock movements. However, the movements, births and slaughter of cattle could not explain the temporal pattern of VTEC O157 prevalence in cattle, despite their inherently distinct seasonality.

Temporal network structures controlling disease spreading

We investigate disease spreading on eight empirical data sets of human contacts (mostly proximity networks recording who is close to whom, at what time). We compare three levels of representations of these data sets: temporal networks, static networks, and a fully connected topology. We notice that the difference between the static and fully connected networks-with respect to time to extinction and average outbreak size-is smaller than between the temporal and static topologies. This suggests that, for these data sets, temporal structures influence disease spreading more than static-network structures. To explain the details in the differences between the representations, we use 32 network measures. This study concurs that long-time temporal structures, like the turnover of nodes and links, are the most important for the spreading dynamics.

Temporal correlation coefficient for directed networks

Previous studies dealing with network theory focused mainly on the static aggregation of edges over specific time window lengths. Thus, most of the dynamic information gets lost. To assess the quality of such a static aggregation the temporal correlation coefficient can be calculated. It measures the overall possibility for an edge to persist between two consecutive snapshots. Up to now, this measure is only defined for undirected networks. Therefore, we introduce the adaption of the temporal correlation coefficient to directed networks. This new methodology enables the distinction between ingoing and outgoing edges. Besides a small example network presenting the single calculation steps, we also calculated the proposed measurements for a real pig trade network to emphasize the importance of considering the edge direction. The farm types at the beginning of the pork supply chain showed clearly higher values for the outgoing temporal correlation coefficient compared to the farm types at the end of the pork supply chain. These farm types showed higher values for the ingoing temporal correlation coefficient. The temporal correlation coefficient is a valuable tool to understand the structural dynamics of these systems, as it assesses the consistency of the edge configuration. The adaption of this measure for directed networks may help to preserve meaningful additional information about the investigated network that might get lost if the edge directions are ignored.