Patterns of host-parasite coinvasion promote enemy release and specialist parasite spillover.

Species invasion and redistribution, driven by climate change and other anthropogenic influences, alter global biodiversity patterns and disrupt ecosystems. As host species move, they can bring their associated parasites with them, potentially infecting resident species, or leave their parasites behind, enhancing their competitive ability in their new ranges. General rules to predict why invading hosts will retain some parasites but not others are relatively unexplored, and the potential predictors are numerous, ranging from parasite life history to host community composition. In this study, we focus on the parasite retention process during host invasion. We used the Global Mammal Parasite Database to identify terrestrial mammal hosts sampled for parasites in both native and non-native ranges. We then selected predictors likely to play a role in parasite retention, such as parasite type, parasite specialism, species composition of the invaded community, and the invading host's phylogenetic or trait-based similarity to the new community. We modelled parasite retention using boosted regression trees, with a suite of 25 predictors describing parasite and host community traits. We further tested the generality of our predictions by cross-validating models on data for other hosts and invasion locations. Our results show that parasite retention is nonrandom and predictable across hosts and invasions. It is broadly shaped by parasite type and parasite specialism, with more specialist parasites that infect many closely related hosts more likely to be retained. This trend is pronounced across parasite types; helminths, however, show a more uniform likelihood of retention regardless of specificity. Overall, we see that most parasites are not retained (11% retained), meaning many invasive species may benefit from enemy release. However, species redistribution does have the potential to spread parasites, and this also has great relevance to understanding conservation implications of species invasions. We see that specialist parasites are most likely to coinvade with their hosts, which suggests that species closely related to the invasive hosts are most likely to be affected by parasite spillover.

Vector communities under global change may exacerbate and redistribute infectious disease risk.

Vector-borne parasites may be transmitted by multiple vector species, resulting in an increased risk of transmission, potentially at larger spatial scales compared to any single vector species. Additionally, the different abilities of patchily distributed vector species to acquire and transmit parasites will lead to varying degrees of transmission risk. Investigation of how vector community composition and parasite transmission change over space due to variation in environmental conditions may help to explain current patterns in diseases but also informs our understanding of how patterns will change under climate and land-use change. We developed a novel statistical approach using a multi-year, spatially extensive case study involving a vector-borne virus affecting white-tailed deer transmitted by Culicoides midges. We characterized the structure of vector communities, established the ecological gradient controlling change in structure, and related the ecology and structure to the amount of disease reporting observed in host populations. We found that vector species largely occur and replace each other as groups, rather than individual species. Moreover, community structure is primarily controlled by temperature ranges, with certain communities being consistently associated with high levels of disease reporting. These communities are essentially composed of species previously undocumented as potential vectors, whereas communities containing putative vector species were largely associated with low levels, or even absence, of disease reporting. We contend that the application of metacommunity ecology to vector-borne infectious disease ecology can greatly aid the identification of transmission hotspots and an understanding of the ecological drivers of parasite transmission risk both now and in the future.

Host and parasite traits predict cross-species parasite acquisition by introduced mammals.

Species invasions and range shifts can lead to novel host-parasite communities, but we lack general rules on which new associations are likely to form. While many studies examine parasite sharing among host species, the directionality of transmission is typically overlooked, impeding our ability to derive principles of parasite acquisition. Consequently, we analysed parasite records from the non-native ranges of 11 carnivore and ungulate species. Using boosted regression trees, we modelled parasite acquisition within each zoogeographic realm of a focal host's non-native range, using a suite of predictors characterizing the parasites themselves and the host community in which they live. We found that higher parasite prevalence among established hosts increases the likelihood of acquisition, particularly for generalist parasites. Non-native host species are also more likely to acquire parasites from established host species to which they are closely related; however, the acquisition of several parasite groups is biased to phylogenetically specialist parasites, indicating potential costs of parasite generalism. Statistical models incorporating these features provide an accurate prediction of parasite acquisition, indicating that measurable host and parasite traits can be used to estimate the likelihood of new host-parasite associations forming. This work provides general rules to help anticipate novel host-parasite associations created by climate change and other anthropogenic influences.

The statistics of epidemic transitions.

Emerging and re-emerging pathogens exhibit very complex dynamics, are hard to model and difficult to predict. Their dynamics might appear intractable. However, new statistical approaches-rooted in dynamical systems and the theory of stochastic processes-have yielded insight into the dynamics of emerging and re-emerging pathogens. We argue that these approaches may lead to new methods for predicting epidemics. This perspective views pathogen emergence and re-emergence as a "critical transition," and uses the concept of noisy dynamic bifurcation to understand the relationship between the system observables and the distance to this transition. Because the system dynamics exhibit characteristic fluctuations in response to perturbations for a system in the vicinity of a critical point, we propose this information may be harnessed to develop early warning signals. Specifically, the motion of perturbations slows as the system approaches the transition.

Vector-borne parasite invasion in communities across space and time.

While vector-borne parasite transmission often operates via generalist-feeding vectors facilitating cross-species transmission in host communities, theory describing the relationship between host species diversity and parasite invasion in these systems is underdeveloped. Host community composition and abundance vary across space and time, generating opportunities for parasite invasion. To explore how host community variation can modify parasite invasion potential, we develop a model for vector-borne parasite transmission dynamics that includes a host community of arbitrary richness and species' abundance. To compare invasion potential across communities, we calculate the community basic reproductive ratio of the parasite. We compare communities comprising a set of host species to their subsets, which allows for flexible scenario building including the introduction of novel host species and species loss. We allow vector abundance to scale with, or be independent of, community size, capturing regulation by feeding opportunities and non-host effects such as limited oviposition sites. Motivated by equivocal data relating host species competency to abundance, we characterize plausible host communities via phenomenological relationships between host species abundance and competency. We identify an underappreciated mechanism whereby changes to communities simultaneously alter average competency and the vector to host ratio and demonstrate that the interaction can profoundly influence invasion potential.

Anticipating epidemic transitions with imperfect data.

Epidemic transitions are an important feature of infectious disease systems. As the transmissibility of a pathogen increases, the dynamics of disease spread shifts from limited stuttering chains of transmission to potentially large scale outbreaks. One proposed method to anticipate this transition are early-warning signals (EWS), summary statistics which undergo characteristic changes as the transition is approached. Although theoretically predicted, their mathematical basis does not take into account the nature of epidemiological data, which are typically aggregated into periodic case reports and subject to reporting error. The viability of EWS for epidemic transitions therefore remains uncertain. Here we demonstrate that most EWS can predict emergence even when calculated from imperfect data. We quantify performance using the area under the curve (AUC) statistic, a measure of how well an EWS distinguishes between numerical simulations of an emerging disease and one which is stationary. Values of the AUC statistic are compared across a range of different reporting scenarios. We find that different EWS respond to imperfect data differently. The mean, variance and first differenced variance all perform well unless reporting error is highly overdispersed. The autocorrelation, autocovariance and decay time perform well provided that the aggregation period of the data is larger than the serial interval and reporting error is not highly overdispersed. The coefficient of variation, skewness and kurtosis are found to be unreliable indicators of emergence. Overall, we find that seven of ten EWS considered perform well for most realistic reporting scenarios. We conclude that imperfect epidemiological data is not a barrier to using EWS for many potentially emerging diseases.

Phylogenetic aggregation increases zoonotic potential of mammalian viruses.

While many viruses of wild mammals are capable of infecting humans, our understanding of zoonotic potential is incomplete. Viruses vary in their degree of generalism, characterized by the phylogenetic relationships of their hosts. Among the dimensions of this phylogenetic landscape, phylogenetic aggregation, which is largely overlooked in studies of parasite host range, emerges in this study as a key predictor of zoonotic status of viruses. Plausibly, viruses that exhibit aggregation, typified by discrete clusters of related host species, may (i) have been able to close the phylogenetic distance to humans, (ii) have subsequently acquired an epidemiologically relevant host and (iii) exhibit relatively high fitness in realized host communities, which are frequently phylogenetically aggregated. These mechanisms associated with phylogenetic aggregation may help explain why correlated fundamental traits, such as the ability of viruses to replicate in the cytoplasm, are associated with zoonoses.

Ebola cases and health system demand in Liberia.

In 2014, a major epidemic of human Ebola virus disease emerged in West Africa, where human-to-human transmission has now been sustained for greater than 12 months. In the summer of 2014, there was great uncertainty about the answers to several key policy questions concerning the path to containment. What is the relative importance of nosocomial transmission compared with community-acquired infection? How much must hospital capacity increase to provide care for the anticipated patient burden? To which interventions will Ebola transmission be most responsive? What must be done to achieve containment? In recent years, epidemic models have been used to guide public health interventions. But, model-based policy relies on high quality causal understanding of transmission, including the availability of appropriate dynamic transmission models and reliable reporting about the sequence of case incidence for model fitting, which were lacking for this epidemic. To investigate the range of potential transmission scenarios, we developed a multi-type branching process model that incorporates key heterogeneities and time-varying parameters to reflect changing human behavior and deliberate interventions in Liberia. Ensembles of this model were evaluated at a set of parameters that were both epidemiologically plausible and capable of reproducing the observed trajectory. Results of this model suggested that epidemic outcome would depend on both hospital capacity and individual behavior. Simulations suggested that if hospital capacity was not increased, then transmission might outpace the rate of isolation and the ability to provide care for the ill, infectious, and dying. Similarly, the model suggested that containment would require individuals to adopt behaviors that increase the rates of case identification and isolation and secure burial of the deceased. As of mid-October, it was unclear that this epidemic would be contained even by 99% hospitalization at the planned hospital capacity. A new version of the model, updated to reflect information collected during October and November 2014, predicts a significantly more constrained set of possible futures. This model suggests that epidemic outcome still depends very heavily on individual behavior. Particularly, if future patient hospitalization rates return to background levels (estimated to be around 70%), then transmission is predicted to remain just below the critical point around Reff = 1. At the higher hospitalization rate of 85%, this model predicts near complete elimination in March to June, 2015.

Predicting cryptic links in host-parasite networks.

Networks are a way to represent interactions among one (e.g., social networks) or more (e.g., plant-pollinator networks) classes of nodes. The ability to predict likely, but unobserved, interactions has generated a great deal of interest, and is sometimes referred to as the link prediction problem. However, most studies of link prediction have focused on social networks, and have assumed a completely censused network. In biological networks, it is unlikely that all interactions are censused, and ignoring incomplete detection of interactions may lead to biased or incorrect conclusions. Previous attempts to predict network interactions have relied on known properties of network structure, making the approach sensitive to observation errors. This is an obvious shortcoming, as networks are dynamic, and sometimes not well sampled, leading to incomplete detection of links. Here, we develop an algorithm to predict missing links based on conditional probability estimation and associated, node-level features. We validate this algorithm on simulated data, and then apply it to a desert small mammal host-parasite network. Our approach achieves high accuracy on simulated and observed data, providing a simple method to accurately predict missing links in networks without relying on prior knowledge about network structure.

Spatiotemporal Fluctuations and Triggers of Ebola Virus Spillover.

Because the natural reservoir of Ebola virus remains unclear and disease outbreaks in humans have occurred only sporadically over a large region, forecasting when and where Ebola spillovers are most likely to occur constitutes a continuing and urgent public health challenge. We developed a statistical modeling approach that associates 37 human or great ape Ebola spillovers since 1982 with spatiotemporally dynamic covariates including vegetative cover, human population size, and absolute and relative rainfall over 3 decades across sub-Saharan Africa. Our model (area under the curve 0.80 on test data) shows that spillover intensity is highest during transitions between wet and dry seasons; overall, high seasonal intensity occurs over much of tropical Africa; and spillover intensity is greatest at high (>1,000/km2) and very low (<100/km2) human population densities compared with intermediate levels. These results suggest strong seasonality in Ebola spillover from wild reservoirs and indicate particular times and regions for targeted surveillance.

Estimating parasite host range.

Estimating the number of host species that a parasite can infect (i.e. host range) provides key insights into the evolution of host specialism and is a central concept in disease ecology. Host range is rarely estimated in real systems, however, because variation in species relative abundance and the detection of rare species makes it challenging to confidently estimate host range. We applied a non-parametric richness indicator to estimate host range in simulated and empirical data, allowing us to assess the influence of sampling heterogeneity and data completeness. After validating our method on simulated data, we estimated parasite host range for a sparsely sampled global parasite occurrence database (Global Mammal Parasite Database) and a repeatedly sampled set of parasites of small mammals from New Mexico (Sevilleta Long Term Ecological Research Program). Estimation accuracy varied strongly with parasite taxonomy, number of parasite occurrence records, and the shape of host species-abundance distribution (i.e. the dominance and rareness of species in the host community). Our findings suggest that between 20% and 40% of parasite host ranges are currently unknown, highlighting a major gap in our understanding of parasite specificity, host-parasite network structure, and parasite burdens.

Global Mammal Parasite Database version 2.0.

Illuminating the ecological and evolutionary dynamics of parasites is one of the most pressing issues facing modern science, and is critical for basic science, the global economy, and human health. Extremely important to this effort are data on the disease-causing organisms of wild animal hosts (including viruses, bacteria, protozoa, helminths, arthropods, and fungi). Here we present an updated version of the Global Mammal Parasite Database, a database of the parasites of wild ungulates (artiodactyls and perissodactyls), carnivores, and primates, and make it available for download as complete flat files. The updated database has more than 24,000 entries in the main data file alone, representing data from over 2700 literature sources. We include data on sampling method and sample sizes when reported, as well as both "reported" and "corrected" (i.e., standardized) binomials for each host and parasite species. Also included are current higher taxonomies and data on transmission modes used by the majority of species of parasites in the database. In the associated metadata we describe the methods used to identify sources and extract data from the primary literature, how entries were checked for errors, methods used to georeference entries, and how host and parasite taxonomies were standardized across the database. We also provide definitions of the data fields in each of the four files that users can download.

Predictability of helminth parasite host range using information on geography, host traits and parasite community structure.

Host-parasite associations are complex interactions dependent on aspects of hosts (e.g. traits, phylogeny or coevolutionary history), parasites (e.g. traits and parasite interactions) and geography (e.g. latitude). Predicting the permissive host set or the subset of the host community that a parasite can infect is a central goal of parasite ecology. Here we develop models that accurately predict the permissive host set of 562 helminth parasites in five different parasite taxonomic groups. We developed predictive models using host traits, host taxonomy, geographic covariates, and parasite community composition, finding that models trained on parasite community variables were more accurate than any other covariate group, even though parasite community covariates only captured a quarter of the variance in parasite community composition. This suggests that it is possible to predict the permissive host set for a given parasite, and that parasite community structure is an important predictor, potentially because parasite communities are interacting non-random assemblages.

The potential for sexual transmission to compromise control of Ebola virus outbreaks.

Recent evidence suggests that sexual contact may give rise to transmission of Ebola virus long after infection has been cleared from blood. We develop a simple mathematical model that incorporates contact transmission and sexual transmission parametrized from data relating to the 2013-2015 West African Ebola epidemic. The model explores scenarios where contact transmission is reduced following infection events, capturing behaviour change, and quantifies how these actions reducing transmission may be compromised by sexual transmission in terms of increasing likelihood, size and duration of outbreaks. We characterize the extent to which sexual transmission operates in terms of the probability of initial infection resolving to sexual infectiousness and the sexual transmission rate, and relate these parameters to the overall case burden. We find that sexual transmission can have large effects on epidemic dynamics (increasing attack ratios from 25% in scenarios without sexual transmission but with contact-transmission-reducing behaviour, up to 80% in equivalent scenarios with sexual transmission).

Vector species richness increases haemorrhagic disease prevalence through functional diversity modulating the duration of seasonal transmission.

Although many parasites are transmitted between hosts by a suite of arthropod vectors, the impact of vector biodiversity on parasite transmission is poorly understood. Positive relationships between host infection prevalence and vector species richness (SR) may operate through multiple mechanisms, including (i) increased vector abundance, (ii) a sampling effect in which species of high vectorial capacity are more likely to occur in species-rich communities, and (iii) functional diversity whereby communities comprised species with distinct phenologies may extend the duration of seasonal transmission. Teasing such mechanisms apart is impeded by a lack of appropriate data, yet could highlight a neglected role for functional diversity in parasite transmission. We used statistical modelling of extensive host, vector and microparasite data to test the hypothesis that functional diversity leading to longer seasonal transmission explained variable levels of disease in a wildlife population. We additionally developed a simple transmission model to guide our expectation of how an increased transmission season translates to infection prevalence. Our study demonstrates that vector SR is associated with increased levels of disease reporting, but not via increases in vector abundance or via a sampling effect. Rather, the relationship operates by extending the length of seasonal transmission, in line with theoretical predictions.

Interspecific Contact and Competition May Affect the Strength and Direction of Disease-Diversity Relationships for Directly Transmitted Microparasites.

The frequency of opportunities for transmission is key to the severity of directly transmitted disease outbreaks in multihost communities. Transmission opportunities for generalist microparasites often arise from competitive and trophic interactions. Additionally, contact heterogeneities within and between species either hinder or promote transmission. General theory incorporating competition and contact heterogeneities for disease-diversity relationships is underdeveloped. Here, we present a formal framework to explore disease-diversity relationships for directly transmitted parasites that infect multiple host species, including influenza viruses, rabies virus, distemper viruses, and hantaviruses. We explicitly include host regulation via intra- and interspecific competition, where the latter can be dependent on or independent of interspecific contact rates (covering resource utilization overlap, habitat selection preferences, and temporal niche partitioning). We examine how these factors interact with frequency- and density-dependent transmission along with traits of the hosts in the assemblage, culminating in the derivation of a relationship describing the propensity for parasite fitness to decrease in species assemblages relative to that in single-host species. This relationship reveals that increases in biodiversity do not necessarily suppress frequency-dependent parasite transmission and that regulation of hosts via interspecific competition does not always lead to a reduction in parasite fitness. Our approach explicitly shows that species identity and ecological interactions between hosts together determine microparasite transmission outcomes in multispecies communities.

Infectious disease transmission and behavioural allometry in wild mammals.

Animals' social and movement behaviours can impact the transmission dynamics of infectious diseases, especially for pathogens transmitted through close contact between hosts or through contact with infectious stages in the environment. Estimating pathogen transmission rates and R0 from natural systems can be challenging. Because host behavioural traits that underlie the transmission process vary predictably with body size, one of the best-studied traits among animals, body size might therefore also predict variation in parasite transmission dynamics. Here, we examine how two host behaviours, social group living and the intensity of habitat use, scale allometrically using comparative data from wild primate, carnivore and ungulate species. We use these empirical relationships to parameterize classical compartment models for infectious micro- and macroparasitic diseases, and examine how the risk of pathogen invasion changes as a function of host behaviour and body size. We then test model predictions using comparative data on parasite prevalence and richness from wild mammals. We report a general pattern suggesting that smaller-bodied mammal species utilizing home ranges more intensively experience greater risk for invasion by environmentally transmitted macroparasites. Conversely, larger-bodied hosts exhibiting a high degree of social group living could be more readily invaded by directly transmitted microparasites. These trends were supported through comparison of micro- and macroparasite species richness across a large number of carnivore, primate and ungulate species, but empirical data on carnivore macroparasite prevalence showed mixed results. Collectively, our study demonstrates that combining host behavioural traits with dynamical models of infectious disease scaled against host body size can generate testable predictions for variation in parasite risk across species; a similar approach might be useful in future work focused on predicting parasite distributions in local host communities.

Refugia and the evolutionary epidemiology of drug resistance.

Drug resistance is a long-standing economic, veterinary and human health concern in human and animal populations. Efficacy of prophylactic drug treatments targeting a particular pathogen is often short-lived, as drug-resistant pathogens evolve and reach high frequency in a treated population. Methods to combat drug resistance are usually costly, including use of multiple drugs that are applied jointly or sequentially, or development of novel classes of drugs. Alternatively, there is growing interest in exploiting untreated host populations, refugia, for the management of drug resistance. Refugia do not experience selection for resistance, and serve as a reservoir for native, drug-susceptible pathogens. The force of infection from refugia may dilute the frequency of resistant pathogens in the treated population, potentially at an acceptable cost in terms of overall disease burden. We examine this concept using a simple mathematical model that captures the core mechanisms of transmission and selection common to many host-pathogen systems. We identify the roles of selection and gene flow in determining the utility of refugia.