Cattle transport network predicts endemic and epidemic foot-and-mouth disease risk on farms in Turkey

The structure of contact networks affects the likelihood of disease spread at the population scale and the risk of infection at any given node. Though this has been well characterized for both theoretical and empirical networks for the spread of epidemics on completely susceptible networks, the long-term impact of network structure on risk of infection with an endemic pathogen, where nodes can be infected more than once, has been less well characterized. Here, we analyze detailed records of the transportation of cattle among farms in Turkey to characterize the global and local attributes of the directed—weighted shipments network between 2007-2012. We then study the correlations between network properties and the likelihood of infection with, or exposure to, foot-and-mouth disease (FMD) over the same time period using recorded outbreaks. The shipments network shows a complex combination of features (local and global) that have not been previously reported in other networks of shipments; i.e. small-worldness, scale-freeness, modular structure, among others. We find that nodes that were either infected or at high risk of infection with FMD (within one link from an infected farm) had disproportionately higher degree, were more central (eigenvector centrality and coreness), and were more likely to be net recipients of shipments compared to those that were always more than 2 links away from an infected farm. High in-degree (i.e. many shipments received) was the best univariate predictor of infection. Low in-coreness (i.e. peripheral nodes) was the best univariate predictor of nodes always more than 2 links away from an infected farm. These results are robust across the three different serotypes of FMD observed in Turkey and during periods of low-endemic prevalence and high-prevalence outbreaks.     

Infectious disease transmission and contact networks in wildlife and livestock

The use of social and contact networks to answer basic and applied questions about infectious disease transmission in wildlife and livestock is receiving increased attention. Through social network analysis, we understand that wild animal and livestock populations, including farmed fish and poultry, often have a heterogeneous contact structure owing to social structure or trade networks. Network modelling is a flexible tool used to capture the heterogeneous contacts of a population in order to test hypotheses about the mechanisms of disease transmission, simulate and predict disease spread, and test disease control strategies. This review highlights how to use animal contact data, including social networks, for network modelling, and emphasizes that researchers should have a pathogen of interest in mind before collecting or using contact data. This paper describes the rising popularity of network approaches for understanding transmission dynamics in wild animal and livestock populations; discusses the common mismatch between contact networks as measured in animal behaviour and relevant parasites to match those networks; and highlights knowledge gaps in how to collect and analyse contact data. Opportunities for the future include increased attention to experiments, pathogen genetic markers and novel computational tools.     

The dynamics of sexual contact networks: effects on disease spread and control

Sexually transmitted pathogens persist in populations despite the availability of biomedical interventions and knowledge of behavioural changes that would reduce individual-level risk. While behavioural risk factors are shared between many sexually transmitted infections, the prevalence of these diseases across different risk groups varies. Understanding this heterogeneity and identifying better control strategies depends on an improved understanding of the complex social contact networks over which pathogens spread. To date, most efforts to study the impact of sexual network structure on disease dynamics have focused on static networks. However, the interaction between the dynamics of partnership formation and dissolution and the dynamics of transmission plays a role, both in restricting the effective network accessible to the pathogen, and in modulating the transmission dynamics. We present a simple method to simulate dynamical networks of sexual partnerships. We inform the model using survey data on sexual attitudes and lifestyles, and investigate how the duration of infectiousness changes the effective contact network over which disease may spread. We then simulate several control strategies: screening, vaccination and behavioural interventions. Previous theory and research has advanced the importance of core groups for spread and control of STD. Our work is consistent with the importance of core groups, but extends this idea to consider how the duration of infectiousness associated with a particular pathogen interacts with host behaviours to define these high risk subpopulations. Characteristics of the parts of the network accessible to the pathogen, which represent the network structure of sexual contacts from the “point of view” of the pathogen, are substantially different from those of the network as a whole. The pathogen itself plays an important role in determining this effective network structure; specifically, we find that if the pathogen’s duration of infectiousness is short, infection is more concentrated in high-activity, high-concurrency individuals even when all other factors are held constant. Widespread screening programmes would be enhanced by follow-up interventions targeting higher-risk individuals, because screening shortens the expected duration of infectiousness and causes a greater relative decrease in prevalence among lower-activity than in higher-activity individuals. Even for pathogens with longer durations of infectiousness, our findings suggest that targeting vaccination and behavioural interventions towards high-activity individuals provides comparable benefits to population-wide interventions.     

How the Dynamics and Structure of Sexual Contact Networks Shape Pathogen Phylogenies

The characteristics of the host contact network over which a pathogen is transmitted affect both epidemic spread and the projected effectiveness of control strategies. Given the importance of understanding these contact networks, it is unfortunate that they are very difficult to measure directly. This challenge has led to an interest in methods to infer information about host contact networks from pathogen phylogenies, because in shaping a pathogen''s opportunities for reproduction, contact networks also shape pathogen evolution. Host networks influence pathogen phylogenies both directly, through governing opportunities for evolution, and indirectly by changing the prevalence and incidence. Here, we aim to separate these two effects by comparing pathogen evolution on different host networks that share similar epidemic trajectories. This approach allows use to examine the direct effects of network structure on pathogen phylogenies, largely controlling for confounding differences arising from population dynamics. We find that networks with more heterogeneous degree distributions yield pathogen phylogenies with more variable cluster numbers, smaller mean cluster sizes, shorter mean branch lengths, and somewhat higher tree imbalance than networks with relatively homogeneous degree distributions. However, in particular for dynamic networks, we find that these direct effects are relatively modest. These findings suggest that the role of the epidemic trajectory, the dynamics of the network and the inherent variability of metrics such as cluster size must each be taken into account when trying to use pathogen phylogenies to understand characteristics about the underlying host contact network.     

Raccoon Social Networks and the Potential for Disease Transmission

Raccoons are an important vector of rabies and other pathogens. The degree to which these pathogens can spread through a raccoon population should be closely linked to association rates between individual raccoons. Most studies of raccoon sociality have found patterns consistent with low levels of social connectivity within populations, thus the likelihood of direct pathogen transmission between raccoons is theoretically low. We used proximity detecting collars and social network metrics to calculate the degree of social connectivity in an urban raccoon population for purposes of estimating potential pathogen spread. In contrast to previous assumptions, raccoon social association networks were highly connected, and all individuals were connected to one large social network during 15 out of 18 months of study. However, these metrics may overestimate the potential for a pathogen to spread through a population, as many of the social connections were based on relatively short contact periods. To more closely reflect varying probabilities of pathogen spread, we censored the raccoon social networks based on the total amount of time spent in close proximity between two individuals per month. As this time criteria for censoring the social networks increased from one to thirty minutes, corresponding measures of network connectivity declined. These findings demonstrate that raccoon populations are much more tightly connected than would have been predicted based on previous studies, but also point out that additional research is needed to calculate more precise transmission probabilities by infected individuals, and determine how disease infection changes normal social behaviors.     

Population genetic structure of a common host predicts the spread of white‐nose syndrome,an emerging infectious disease in bats

Landscape complexity influences patterns of animal dispersal, which in turn may affect both gene flow and the spread of pathogens. White‐nose syndrome (WNS) is an introduced fungal disease that has spread rapidly throughout eastern North America, causing massive mortality in bat populations. We tested for a relationship between the population genetic structure of the most common host, the little brown myotis (Myotis lucifugus), and the geographic spread of WNS to date by evaluating logistic regression models of WNS risk among hibernating colonies in eastern North America. We hypothesized that risk of WNS to susceptible host colonies should increase with both geographic proximity and genetic similarity, reflecting historical connectivity, to infected colonies. Consistent with this hypothesis, inclusion of genetic distance between infected and susceptible colonies significantly improved models of disease spread, capturing heterogeneity in the spatial expansion of WNS despite low levels of genetic differentiation among eastern populations. Expanding our genetic analysis to the continental range of little brown myotis reveals strongly contrasting patterns of population structure between eastern and western North America. Genetic structure increases markedly moving westward into the northern Great Plains, beyond the current distribution of WNS. In western North America, genetic differentiation of geographically proximate populations often exceeds levels observed across the entire eastern region, suggesting infrequent and/or locally restricted dispersal, and thus relatively limited opportunities for pathogen introduction in western North America. Taken together, our analyses suggest a possibly slower future rate of spread of the WNS pathogen, at least as mediated by little brown myotis.     

Contribution of Company Affiliation and Social Contacts to Risk Estimates of Between-Farm Transmission of Avian Influenza

Background

Models of between-farm transmission of pathogens have identified service vehicles and social groups as risk factors mediating the spread of infection. Because of high levels of economic organization in much of the poultry industry, we examined the importance of company affiliation, as distinct from social contacts, in a model of the potential spread of avian influenza among broiler poultry farms in a poultry-dense region in the United States. The contribution of company affiliation to risk of between-farm disease transmission has not been previously studied.

Methodology/Principal Findings

We obtained data on the nature and frequency of business and social contacts through a national survey of broiler poultry growers in the United States. Daily rates of contact were estimated using Monte Carlo analysis. Stochastic modeling techniques were used to estimate the exposure risk posed by a single infectious farm to other farms in the region and relative risk of exposure for farms under different scenarios. The mean daily rate of vehicular contact was 0.82 vehicles/day. The magnitude of exposure risk ranged from <1% to 25% under varying parameters. Risk of between-farm transmission was largely driven by company affiliation, with farms in the same company group as the index farm facing as much as a 5-fold increase in risk compared to farms contracted with different companies. Employment of part-time workers contributed to significant increases in risk in most scenarios, notably for farms who hired day-laborers. Social visits were significantly less important in determining risk.

Conclusions/Significance

Biosecurity interventions should be based on information on industry structure and company affiliation, and include part-time workers as potentially unrecognized sources of viral transmission. Modeling efforts to understand pathogen transmission in the context of industrial food animal production should consider company affiliation in addition to geospatial factors and pathogen characteristics. Restriction of social contacts among farmers may be less useful in reducing between-farm transmission.     

Prevalence and infection intensity of Nosema in honey bee (Apis mellifera L.) colonies in Virginia

Nosema ceranae is a recently described pathogen of Apis mellifera and Apis cerana. Relatively little is known about the distribution or prevalence of N. ceranae in the United States. To determine the prevalence and potential impact of this new pathogen on honey bee colonies in Virginia, over 300 hives were sampled across the state. The samples were analyzed microscopically for Nosema spores and for the presence of the pathogen using real-time PCR. Our studies indicate that N. ceranae is the dominant species in Virginia with an estimated 69.3% of hives infected. Nosema apis infections were only observed at very low levels (2.7%), and occurred only as co-infections with N. ceranae. Traditional diagnoses based on spore counts alone do not provide an accurate indication of colony infections. We found that 51.1% of colonies that did not have spores present in the sample were infected with N. ceranae when analyzed by real-time PCR. In hives that tested positive for N. ceranae, average C_T values were used to diagnose a hive as having a low, moderate, or a heavy infection intensity. Most infected colonies had low-level infections (73%), but 11% of colonies had high levels of infection and 16% had moderate level infections. The prevalence and mean levels of infection were similar in different regions of the state.     

Combining Epidemiological and Genetic Networks Signifies the Importance of Early Treatment in HIV-1 Transmission

Inferring disease transmission networks is important in epidemiology in order to understand and prevent the spread of infectious diseases. Reconstruction of the infection transmission networks requires insight into viral genome data as well as social interactions. For the HIV-1 epidemic, current research either uses genetic information of patients'' virus to infer the past infection events or uses statistics of sexual interactions to model the network structure of viral spreading. Methods for a reliable reconstruction of HIV-1 transmission dynamics, taking into account both molecular and societal data are still lacking. The aim of this study is to combine information from both genetic and epidemiological scales to characterize and analyse a transmission network of the HIV-1 epidemic in central Italy.We introduce a novel filter-reduction method to build a network of HIV infected patients based on their social and treatment information. The network is then combined with a genetic network, to infer a hypothetical infection transmission network. We apply this method to a cohort study of HIV-1 infected patients in central Italy and find that patients who are highly connected in the network have longer untreated infection periods. We also find that the network structures for homosexual males and heterosexual populations are heterogeneous, consisting of a majority of ‘peripheral nodes’ that have only a few sexual interactions and a minority of ‘hub nodes’ that have many sexual interactions. Inferring HIV-1 transmission networks using this novel combined approach reveals remarkable correlations between high out-degree individuals and longer untreated infection periods. These findings signify the importance of early treatment and support the potential benefit of wide population screening, management of early diagnoses and anticipated antiretroviral treatment to prevent viral transmission and spread. The approach presented here for reconstructing HIV-1 transmission networks can have important repercussions in the design of intervention strategies for disease control.     

Initiation of Batrachochytrium dendrobatidis infection in the absence of physical contact with infected hosts – a field study in a high altitude lake

Understanding transmission is a critical prerequisite for predicting disease dynamics and impacts on host populations. It is well established that Batrachochytrium dendrobatidis (Bd), the amphibian fungal pathogen responsible for chytridiomycosis, can be transmitted directly, through physical contact with an infected host. However, indirect pathways of transmission remain poorly investigated. We conducted a five‐week long field infection experiment at a high altitude mountain lake in the French Pyrenees to investigate Bd transmission pathways in larval midwife toads Alytes obstetricans. Uninfected naïve tadpoles were co‐housed either with infected tadpoles (direct and indirect transmission) or with uninfected ones (indirect transmission only). We found that physical contact with an infected host is not necessary for initial infection with Bd and that all tadpoles became infected after only four weeks. However, physical contact with infected tadpoles led to a faster spread within a tadpole group and resulted in higher Bd loads and subsequently higher mortality. Our findings clearly demonstrate that in A. obstetricans, Bd can quickly spread in a population even without physical contact. Our experiment therefore stresses the importance of indirect transmission of Bd zoospores in infected lakes for disease dynamics, especially when a reservoir species such as A. obstetricans is present.     

Epidemic threshold and network structure: The interplay of probability of transmission and of persistence in small-size directed networks

A growing number of studies are investigating the effect of contact structure on the dynamics of epidemics in large-scale complex networks. Whether findings thus obtained apply also to networks of small size, and thus to many real-world biological applications, is still an open question. We use numerical simulations of disease spread in directed networks of 100 individual nodes with a constant number of links. We show that, no matter the type of network structure (local, small-world, random and scale-free), there is a linear threshold determined by the probability of infection transmission between connected nodes and the probability of infection persistence in an infected node. The threshold is significantly lower for scale-free networks compared to local, random and small-world ones only if super-connected nodes have a higher number of links both to and from other nodes. The starting point, the node at which the epidemic starts, does not affect the threshold conditions, but has a marked influence on the final size of the epidemic in all kinds of network. There is evidence that contact structure has an influence on the average final size of an epidemic across all starting nodes, with significantly lower values in scale-free networks at equilibrium. Simulations in scale-free networks show a distinctive time-series pattern, which, if found in a real epidemic, can be used to infer the underlying network structure. The findings have relevance also for meta-population ecology and species conservation.     

Coevolutionary cycling of host sociality and pathogen virulence in contact networks

Infectious diseases may place strong selection on the social organization of animals. Conversely, the structure of social systems can influence the evolutionary trajectories of pathogens. While much attention has focused on the evolution of host sociality or pathogen virulence separately, few studies have looked at their coevolution. Here we use an agent-based simulation to explore host-pathogen coevolution in social contact networks. Our results indicate that under certain conditions, both host sociality and pathogen virulence exhibit continuous cycling. The way pathogens move through the network (e.g., their interhost transmission and probability of superinfection) and the structure of the network can influence the existence and form of cycling.     

Host and pathogen ecology drive the seasonal dynamics of a fungal disease,white-nose syndrome

Seasonal patterns in pathogen transmission can influence the impact of disease on populations and the speed of spatial spread. Increases in host contact rates or births drive seasonal epidemics in some systems, but other factors may occasionally override these influences. White-nose syndrome, caused by the emerging fungal pathogen Pseudogymnoascus destructans, is spreading across North America and threatens several bat species with extinction. We examined patterns and drivers of seasonal transmission of P. destructans by measuring infection prevalence and pathogen loads in six bat species at 30 sites across the eastern United States. Bats became transiently infected in autumn, and transmission spiked in early winter when bats began hibernating. Nearly all bats in six species became infected by late winter when infection intensity peaked. In summer, despite high contact rates and a birth pulse, most bats cleared infections and prevalence dropped to zero. These data suggest the dominant driver of seasonal transmission dynamics was a change in host physiology, specifically hibernation. Our study is the first, to the best of our knowledge, to describe the seasonality of transmission in this emerging wildlife disease. The timing of infection and fungal growth resulted in maximal population impacts, but only moderate rates of spatial spread.     

Brown spider monkeys (Ateles hybridus): a model for differentiating the role of social networks and physical contact on parasite transmission dynamics

Elevated risk of disease transmission is considered a major cost of sociality, although empirical evidence supporting this idea remains scant. Variation in spatial cohesion and the occurrence of social interactions may have profound implications for patterns of interindividual parasite transmission. We used a social network approach to shed light on the importance of different aspects of group-living (i.e. within-group associations versus physical contact) on patterns of parasitism in a neotropical primate, the brown spider monkey (Ateles hybridus), which exhibits a high degree of fission–fusion subgrouping. We used daily subgroup composition records to create a ‘proximity’ network, and built a separate ‘contact’ network using social interactions involving physical contact. In the proximity network, connectivity between individuals was homogeneous, whereas the contact network highlighted high between-individual variation in the extent to which animals had physical contact with others, which correlated with an individual''s age and sex. The gastrointestinal parasite species richness of highly connected individuals was greater than that of less connected individuals in the contact network, but not in the proximity network. Our findings suggest that among brown spider monkeys, physical contact impacts the spread of several common parasites and supports the idea that pathogen transmission is one cost associated with social contact.     

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.     

Coexistence and specialization of pathogen strains on contact networks

The coexistence of different pathogen strains has implications for pathogen variability and disease control and has been explained in a number of different ways. We use contact networks, which represent interactions between individuals through which infection could be transmitted, to investigate strain coexistence. For sexually transmitted diseases the structure of contact networks has received detailed study and has been shown to be a vital determinant of the epidemiological dynamics. By using analytical pairwise models and stochastic simulations, we demonstrate that network structure also has a profound influence on the interaction between pathogen strains. In particular, when the population is serially monogamous, fully cross-reactive strains can coexist, with different strains dominating in network regions with different characteristics. Furthermore, we observe specialization of different strains in different risk groups within the network, suggesting the existence of diverging evolutionary pressures.     

Contact networks structured by sex underpin sex‐specific epidemiology of infection

Contact networks are fundamental to the transmission of infection and host sex often affects the acquisition and progression of infection. However, the epidemiological impacts of sex‐related variation in animal contact networks have rarely been investigated. We test the hypothesis that sex‐biases in infection are related to variation in multilayer contact networks structured by sex in a population of European badgers Meles meles naturally infected with Mycobacterium bovis. Our key results are that male‐male and between‐sex networks are structured at broader spatial scales than female‐female networks and that in male‐male and between‐sex contact networks, but not female‐female networks, there is a significant relationship between infection and contacts with individuals in other groups. These sex differences in social behaviour may underpin male‐biased acquisition of infection and may result in males being responsible for more between‐group transmission. This highlights the importance of sex‐related variation in host behaviour when managing animal diseases.     

Comparison of Contact Patterns Relevant for Transmission of Respiratory Pathogens in Thailand and the Netherlands Using Respondent-Driven Sampling

Understanding infection dynamics of respiratory diseases requires the identification and quantification of behavioural, social and environmental factors that permit the transmission of these infections between humans. Little empirical information is available about contact patterns within real-world social networks, let alone on differences in these contact networks between populations that differ considerably on a socio-cultural level. Here we compared contact network data that were collected in the Netherlands and Thailand using a similar online respondent-driven method. By asking participants to recruit contact persons we studied network links relevant for the transmission of respiratory infections. We studied correlations between recruiter and recruited contacts to investigate mixing patterns in the observed social network components. In both countries, mixing patterns were assortative by demographic variables and random by total numbers of contacts. However, in Thailand participants reported overall more contacts which resulted in higher effective contact rates. Our findings provide new insights on numbers of contacts and mixing patterns in two different populations. These data could be used to improve parameterisation of mathematical models used to design control strategies. Although the spread of infections through populations depends on more factors, found similarities suggest that spread may be similar in the Netherlands and Thailand.     

Spread of plague among black-tailed prairie dogs is associated with colony spatial characteristics

Sylvatic plague (Yersinia pestis) is an exotic pathogen that is highly virulent in black-tailed prairie dogs (Cynomys ludovicianus) and causes widespread colony losses and individual mortality rates >95%. We investigated colony spatial characteristics that may influence inter-colony transmission of plague at 3 prairie dog colony complexes in the Great Plains. The 4 spatial characteristics we considered include: colony size, Euclidean distance to nearest neighboring colony, colony proximity index, and distance to nearest drainage (dispersal) corridor. We used multi-state mark–recapture models to determine the relationship between these colony characteristics and probability of plague transmission among prairie dog colonies. Annual mapping of colonies and mark–recapture analyses of disease dynamics in natural colonies led to 4 main results: 1) plague outbreaks exhibited high spatial and temporal variation, 2) the site of initiation of epizootic plague may have substantially influenced the subsequent inter-colony spread of plague, 3) the long-term effect of plague on individual colonies differed among sites because of how individuals and colonies were distributed, and 4) colony spatial characteristics were related to the probability of infection at all sites although the relative importance and direction of relationships varied among sites. Our findings suggest that conventional prairie dog conservation management strategies, including promoting large, highly connected colonies, may need to be altered in the presence of plague. © 2011 The Wildlife Society     

Covariation between the physiological and behavioral components of pathogen transmission: host heterogeneity determines epidemic outcomes

Although heterogeneity in contact rate, physiology, and behavioral response to infection have all been empirically demonstrated in host–pathogen systems, little is known about how interactions between individual variation in behavior and physiology scale‐up to affect pathogen transmission at a population level. The objective of this study is to evaluate how covariation between the behavioral and physiological components of transmission might affect epidemic outcomes in host populations. We tested the consequences of contact rate covarying with susceptibility, infectiousness, and infection status using an individual‐based, dynamic network model where individuals initiate and terminate contacts with conspecifics based on their behavioral predispositions and their infection status. Our results suggest that both heterogeneity in physiology and subsequent covariation of physiology with contact rate could powerfully influence epidemic dynamics. Overall, we found that 1) individual variability in susceptibility and infectiousness can reduce the expected maximum prevalence and increase epidemic variability; 2) when contact rate and susceptibility or infectiousness negatively covary, it takes substantially longer for epidemics to spread throughout the population, and rates of epidemic spread remained suppressed even for highly transmissible pathogens; and 3) reductions in contact rate resulting from infection‐induced behavioral changes can prevent the pathogen from reaching most of the population. These effects were strongest for theoretical pathogens with lower transmissibility and for populations where the observed variation in contact rate was higher, suggesting that such heterogeneity may be most important for less infectious, more chronic diseases in wildlife. Understanding when and how variability in pathogen transmission should be modelled is a crucial next step for disease ecology.