Host adaptation to novel pathogen introduction: Predicting conditions that promote evolutionary rescue.

Novel pathogen introduction can have drastic consequences for naive host populations, and outcomes can be difficult to predict. Evolutionary rescue (ER) provides a foundation for understanding whether hosts are driven to extinction or survive via adaptation. Currently, patterns of host population dynamics alongside evidence of adaptation are used to infer ER. However, the gap between established ER theory and complexity inherent in natural systems makes interpreting empirical patterns difficult because they can be confounded with ecological drivers of survival under current theory. To bridge this gap, we expand ER theory to include biological selective agents, such as pathogens. We find birth processes to be more important than previously theorised in determining ER potential. We employ a novel framework evaluating ER potential within natural systems and gain ability to identify system characteristics that make ER possible. Identifying these characteristics allows a shift from retrospective observation to a predictive mindset, and our findings suggest that ER occurrence may be more limited than previously thought. We use the plague system of Yersinia pestis infecting Cynomys ludovicianus (black-tailed prairie dogs) and Spermophilus beecheyi (California ground squirrels) as a case study.

Realistic assumptions about spatial locations and clustering of premises matter for models of foot-and-mouth disease spread in the United States.

Spatially explicit livestock disease models require demographic data for individual farms or premises. In the U.S., demographic data are only available aggregated at county or coarser scales, so disease models must rely on assumptions about how individual premises are distributed within counties. Here, we addressed the importance of realistic assumptions for this purpose. We compared modeling of foot and mouth disease (FMD) outbreaks using simple randomization of locations to premises configurations predicted by the Farm Location and Agricultural Production Simulator (FLAPS), which infers location based on features such as topography, land-cover, climate, and roads. We focused on three premises-level Susceptible-Exposed-Infectious-Removed models available from the literature, all using the same kernel approach but with different parameterizations and functional forms. By computing the basic reproductive number of the infection (R0) for both FLAPS and randomized configurations, we investigated how spatial locations and clustering of premises affects outbreak predictions. Further, we performed stochastic simulations to evaluate if identified differences were consistent for later stages of an outbreak. Using Ripley's K to quantify clustering, we found that FLAPS configurations were substantially more clustered at the scales relevant for the implemented models, leading to a higher frequency of nearby premises compared to randomized configurations. As a result, R0 was typically higher in FLAPS configurations, and the simulation study corroborated the pattern for later stages of outbreaks. Further, both R0 and simulations exhibited substantial spatial heterogeneity in terms of differences between configurations. Thus, using realistic assumptions when de-aggregating locations based on available data can have a pronounced effect on epidemiological predictions, affecting if, where, and to what extent FMD may invade the population. We conclude that methods such as FLAPS should be preferred over randomization approaches.

Continental-scale dynamics of avian influenza in U.S. waterfowl are driven by demography, migration, and temperature.

Emerging diseases of wildlife origin are increasingly spilling over into humans and domestic animals. Surveillance and risk assessments for transmission between these populations are informed by a mechanistic understanding of the pathogens in wildlife reservoirs. For avian influenza viruses (AIV), much observational and experimental work in wildlife has been conducted at local scales, yet fully understanding their spread and distribution requires assessing the mechanisms acting at both local, (e.g., intrinsic epidemic dynamics), and continental scales, (e.g., long-distance migration). Here, we combined a large, continental-scale data set on low pathogenic, Type A AIV in the United States with a novel network-based application of bird banding/recovery data to investigate the migration-based drivers of AIV and their relative importance compared to well-characterized local drivers (e.g., demography, environmental persistence). We compared among regression models reflecting hypothesized ecological processes and evaluated their ability to predict AIV in space and time using within and out-of-sample validation. We found that predictors of AIV were associated with multiple mechanisms at local and continental scales. Hypotheses characterizing local epidemic dynamics were strongly supported, with age, the age-specific aggregation of migratory birds in an area and temperature being the best predictors of infection. Hypotheses defining larger, network-based features of the migration processes, such as clustering or between-cluster mixing explained less variation but were also supported. Therefore, our results support a role for local processes in driving the continental distribution of AIV.

Predicting functional responses in agro-ecosystems from animal movement data to improve management of invasive pests.

Functional responses describe how changing resource availability affects consumer resource use, thus providing a mechanistic approach to prediction of the invasibility and potential damage of invasive alien species (IAS). However, functional responses can be context dependent, varying with resource characteristics and availability, consumer attributes, and environmental variables. Identifying context dependencies can allow invasion and damage risk to be predicted across different ecoregions. Understanding how ecological factors shape the functional response in agro-ecosystems can improve predictions of hotspots of highest impact and inform strategies to mitigate damage across locations with varying crop types and availability. We linked heterogeneous movement data across different agro-ecosystems to predict ecologically driven variability in the functional responses. We applied our approach to wild pigs (Sus scrofa), one of the most successful and detrimental IAS worldwide where agricultural resource depredation is an important driver of spread and establishment. We used continental-scale movement data within agro-ecosystems to quantify the functional response of agricultural resources relative to availability of crops and natural forage. We hypothesized that wild pigs would selectively use crops more often when natural forage resources were low. We also examined how individual attributes such as sex, crop type, and resource stimulus such as distance to crops altered the magnitude of the functional response. There was a strong agricultural functional response where crop use was an accelerating function of crop availability at low density (Type III) and was highly context dependent. As hypothesized, there was a reduced response of crop use with increasing crop availability when non-agricultural resources were more available, emphasizing that crop damage levels are likely to be highly heterogeneous depending on surrounding natural resources and temporal availability of crops. We found significant effects of crop type and sex, with males spending 20% more time and visiting crops 58% more often than females, and both sexes showing different functional responses depending on crop type. Our application demonstrates how commonly collected animal movement data can be used to understand context dependencies in resource use to improve our understanding of pest foraging behavior, with implications for prioritizing spatiotemporal hotspots of potential economic loss in agro-ecosystems.

Need for speed: An optimized gridding approach for spatially explicit disease simulations.

Numerical models for simulating outbreaks of infectious diseases are powerful tools for informing surveillance and control strategy decisions. However, large-scale spatially explicit models can be limited by the amount of computational resources they require, which poses a problem when multiple scenarios need to be explored to provide policy recommendations. We introduce an easily implemented method that can reduce computation time in a standard Susceptible-Exposed-Infectious-Removed (SEIR) model without introducing any further approximations or truncations. It is based on a hierarchical infection process that operates on entire groups of spatially related nodes (cells in a grid) in order to efficiently filter out large volumes of susceptible nodes that would otherwise have required expensive calculations. After the filtering of the cells, only a subset of the nodes that were originally at risk are then evaluated for actual infection. The increase in efficiency is sensitive to the exact configuration of the grid, and we describe a simple method to find an estimate of the optimal configuration of a given landscape as well as a method to partition the landscape into a grid configuration. To investigate its efficiency, we compare the introduced methods to other algorithms and evaluate computation time, focusing on simulated outbreaks of foot-and-mouth disease (FMD) on the farm population of the USA, the UK and Sweden, as well as on three randomly generated populations with varying degree of clustering. The introduced method provided up to 500 times faster calculations than pairwise computation, and consistently performed as well or better than other available methods. This enables large scale, spatially explicit simulations such as for the entire continental USA without sacrificing realism or predictive power.

Inferring infection hazard in wildlife populations by linking data across individual and population scales.

Our ability to infer unobservable disease-dynamic processes such as force of infection (infection hazard for susceptible hosts) has transformed our understanding of disease transmission mechanisms and capacity to predict disease dynamics. Conventional methods for inferring FOI estimate a time-averaged value and are based on population-level processes. Because many pathogens exhibit epidemic cycling and FOI is the result of processes acting across the scales of individuals and populations, a flexible framework that extends to epidemic dynamics and links within-host processes to FOI is needed. Specifically, within-host antibody kinetics in wildlife hosts can be short-lived and produce patterns that are repeatable across individuals, suggesting individual-level antibody concentrations could be used to infer time since infection and hence FOI. Using simulations and case studies (influenza A in lesser snow geese and Yersinia pestis in coyotes), we argue that with careful experimental and surveillance design, the population-level FOI signal can be recovered from individual-level antibody kinetics, despite substantial individual-level variation. In addition to improving inference, the cross-scale quantitative antibody approach we describe can reveal insights into drivers of individual-based variation in disease response, and the role of poorly understood processes such as secondary infections, in population-level dynamics of disease.

Invasion of two tick-borne diseases across New England: harnessing human surveillance data to capture underlying ecological invasion processes.

Modelling the spatial spread of vector-borne zoonotic pathogens maintained in enzootic transmission cycles remains a major challenge. The best available spatio-temporal data on pathogen spread often take the form of human disease surveillance data. By applying a classic ecological approach-occupancy modelling-to an epidemiological question of disease spread, we used surveillance data to examine the latent ecological invasion of tick-borne pathogens. Over the last half-century, previously undescribed tick-borne pathogens including the agents of Lyme disease and human babesiosis have rapidly spread across the northeast United States. Despite their epidemiological importance, the mechanisms of tick-borne pathogen invasion and drivers underlying the distinct invasion trajectories of the co-vectored pathogens remain unresolved. Our approach allowed us to estimate the unobserved ecological processes underlying pathogen spread while accounting for imperfect detection of human cases. Our model predicts that tick-borne diseases spread in a diffusion-like manner with occasional long-distance dispersal and that babesiosis spread exhibits strong dependence on Lyme disease.

Ecological Traits Driving the Outbreaks and Emergence of Zoonotic Pathogens.

Infectious diseases that are transmitted from wildlife hosts to humans, such as the Ebola virus and MERS virus, can be difficult to understand because the pathogens emerge from complex multifaceted ecological interactions. We use a wildlife-pathogen system-prairie dogs (Cynomys ludovicianus) and the plague bacterium (Yersinia pestis)-to describe aspects of disease ecology that apply to many cases of emerging infectious disease. We show that the monitoring and surveillance of hosts and vectors during the buildup to disease outbreaks are crucial for understanding pathogen-transmission dynamics and that a community-ecology framework is important to identify reservoir hosts. Incorporating multidisciplinary approaches and frameworks may improve wildlife-pathogen surveillance and our understanding of seemingly sporadic and rare pathogen outbreaks.

Identification of migratory bird flyways in North America using community detection on biological networks.

Migratory behavior of waterfowl populations in North America has traditionally been broadly characterized by four north-south flyways, and these flyways have been central to the management of waterfowl populations for more than 80 yr. However, previous flyway characterizations are not easily updated with current bird movement data and fail to provide assessments of the importance of specific geographical regions to the identification of flyways. Here, we developed a network model of migratory movement for four waterfowl species, Mallard (Anas platyrhnchos), Northern Pintail (A. acuta), American Green-winged Teal (A. carolinensis), and Canada Goose (Branta canadensis), in North America, using bird band and recovery data. We then identified migratory flyways using a community detection algorithm and characterized the importance of smaller geographic regions in identifying flyways using a novel metric, the consolidation factor. We identified four main flyways for Mallards, Northern Pintails, and American Green-winged Teal, with the flyway identification in Canada Geese exhibiting higher complexity. For Mallards, flyways were relatively consistent through time. However, consolidation factors revealed that for Mallards and Green-winged Teal, the presumptive Mississippi flyway was potentially a zone of high mixing between other flyways. Our results demonstrate that the network approach provides a robust method for flyway identification that is widely applicable given the relatively minimal data requirements and is easily updated with future movement data to reflect changes in flyway definitions and management goals.

Network analysis of host-virus communities in bats and rodents reveals determinants of cross-species transmission.

Bats are natural reservoirs of several important emerging viruses. Cross-species transmission appears to be quite common among bats, which may contribute to their unique reservoir potential. Therefore, understanding the importance of bats as reservoirs requires examining them in a community context rather than concentrating on individual species. Here, we use a network approach to identify ecological and biological correlates of cross-species virus transmission in bats and rodents, another important host group. We show that given our current knowledge the bat viral sharing network is more connected than the rodent network, suggesting viruses may pass more easily between bat species. We identify host traits associated with important reservoir species: gregarious bats are more likely to share more viruses and bats which migrate regionally are important for spreading viruses through the network. We identify multiple communities of viral sharing within bats and rodents and highlight potential species traits that can help guide studies of novel pathogen emergence.

Host and viral ecology determine bat rabies seasonality and maintenance.

Rabies is an acute viral infection that is typically fatal. Most rabies modeling has focused on disease dynamics and control within terrestrial mammals (e.g., raccoons and foxes). As such, rabies in bats has been largely neglected until recently. Because bats have been implicated as natural reservoirs for several emerging zoonotic viruses, including SARS-like corona viruses, henipaviruses, and lyssaviruses, understanding how pathogens are maintained within a population becomes vital. Unfortunately, little is known about maintenance mechanisms for any pathogen in bat populations. We present a mathematical model parameterized with unique data from an extensive study of rabies in a Colorado population of big brown bats (Eptesicus fuscus) to elucidate general maintenance mechanisms. We propose that life history patterns of many species of temperate-zone bats, coupled with sufficiently long incubation periods, allows for rabies virus maintenance. Seasonal variability in bat mortality rates, specifically low mortality during hibernation, allows long-term bat population viability. Within viable bat populations, sufficiently long incubation periods allow enough infected individuals to enter hibernation and survive until the following year, and hence avoid an epizootic fadeout of rabies virus. We hypothesize that the slowing effects of hibernation on metabolic and viral activity maintains infected individuals and their pathogens until susceptibles from the annual birth pulse become infected and continue the cycle. This research provides a context to explore similar host ecology and viral dynamics that may explain seasonal patterns and maintenance of other bat-borne diseases.

A comparison of bats and rodents as reservoirs of zoonotic viruses: are bats special?

Bats are the natural reservoirs of a number of high-impact viral zoonoses. We present a quantitative analysis to address the hypothesis that bats are unique in their propensity to host zoonotic viruses based on a comparison with rodents, another important host order. We found that bats indeed host more zoonotic viruses per species than rodents, and we identified life-history and ecological factors that promote zoonotic viral richness. More zoonotic viruses are hosted by species whose distributions overlap with a greater number of other species in the same taxonomic order (sympatry). Specifically in bats, there was evidence for increased zoonotic viral richness in species with smaller litters (one young), greater longevity and more litters per year. Furthermore, our results point to a new hypothesis to explain in part why bats host more zoonotic viruses per species: the stronger effect of sympatry in bats and more viruses shared between bat species suggests that interspecific transmission is more prevalent among bats than among rodents. Although bats host more zoonotic viruses per species, the total number of zoonotic viruses identified in bats (61) was lower than in rodents (68), a result of there being approximately twice the number of rodent species as bat species. Therefore, rodents should still be a serious concern as reservoirs of emerging viruses. These findings shed light on disease emergence and perpetuation mechanisms and may help lead to a predictive framework for identifying future emerging infectious virus reservoirs.

Persistence of black-tailed prairie-dog populations affected by plague in northern Colorado, USA.

The spatial distribution of prairie dog (Cynomys ludovicianus) colonies in North America has changed from large, contiguous populations to small, isolated colonies in metapopulations. One factor responsible for this drastic change in prairie-dog population structure is plague (caused by the bacterium Yersinia pestis). We fit stochastic patch occupancy models to 20 years of prairie-dog colony occupancy data from two discrete metapopulations (west and east) in the Pawnee National Grassland in Colorado, USA, that differ in connectivity among suitable habitat patches. We conducted model selection between two hypothesized modes of plague movement: independent of prairie-dog dispersal (colony-area) vs. plague movement consistent with prairie-dog dispersal (connectivity to extinct colonies). The best model, which fit the data well (area under the curve [AUC]: 0.94 west area; 0.79 east area), revealed that over time the proportion of extant colonies was better explained by colony size than by connectivity to extinct (plagued) colonies. The idea that prairie dogs are not likely to be the main vector that spreads Y. pestis across the landscape is supported by the observation that colony extinctions are primarily caused by plague, prairie-dog dispersal is short range, and connectivity to extinct colonies was not selected as a factor in the models. We also conducted simulations with the best model to examine long-term patterns of colony occupancy and persistence of prairie-dog metapopulations. In the case where the metapopulations persist, our model predicted that the western metapopulation would have a colony occupancy rate approximately 2.5 times higher than that of the eastern metapopulation (-50% occupied colonies vs. 20%) in 50 years, but that the western metapopulation has -80% chance of extinction in 100 years while the eastern metapopulation has a less than 25% chance. Extinction probability of individual colonies depended on the frequency with which colonies of the same size class occurred in the metapopulation. Thus, the long-term persistence of prairie-dog metapopulations depended on specific details of the metapopulation.

Minimizing the threat of pandemic emergence from avian influenza in poultry systems.

BACKGROUND: Live-animal markets are a culturally important feature of meat distribution chains in many populations, yet they provide an opportunity for the maintenance and transmission of potentially emergent zoonotic pathogens. The ongoing human outbreak of avian H7N9 in China highlights the need for increased surveillance and control in these live-bird markets (LBMs). DISCUSSION: Closure of retail markets in affected areas rapidly decreased human cases to rare, sporadic occurrence, but little attention has been paid thus far to the role of upstream elements of the poultry distribution chain such as wholesale markets. This could partly explain why transmission in poultry populations has not been eliminated more broadly. We present surveillance data from both wholesale live-bird markets (wLBMs) and rLBMs in Shantou, China (from 2004-2006), and call on disease-dynamic theory to illustrate why closing rLBMs has only minor effects on the overall volume of transmission. We show that the length of time birds stay in rLBMs can severely limit transmission there, but that the system-wide effect may be reduced substantially by high levels of transmission upstream of retail markets. SUMMARY: Management plans that minimize transmission throughout the entire poultry supply chain are essential for minimizing exposure to the public. These include reducing stay-time of birds in markets to 1 day, standardizing poultry supply chains to limit transmission in pre-retail settings, and monitoring strains with epidemiological traits that pose a high risk of emergence. These actions will further limit human exposure to extant viruses and reduce the likelihood of the emergence of novel strains by decreasing the overall volume of transmission.

An exploration of within-herd dynamics of a transboundary livestock disease: A foot and mouth disease case study.

Transboundary livestock diseases are a high priority for policy makers because of the serious economic burdens associated with infection. In order to make well informed preparedness and response plans, policy makers often utilize mathematical models to understand possible outcomes of different control strategies and outbreak scenarios. Many of these models focus on the transmission between herds and the overall trajectory of the outbreak. While the course of infection within herds has not been the focus of the majority of models, a thorough understanding of within-herd dynamics can provide valuable insight into a disease system by providing information on herd-level biological properties of the infection, which can be used to inform decision making in both endemic and outbreak settings and to inform larger between-herd models. In this study, we develop three stochastic simulation models to study within-herd foot and mouth disease dynamics and the implications of different empirical data-based assumptions about the timing of the onset of infectiousness and clinical signs. We also study the influence of herd size and the proportion of the herd that is initially infected on the outcome of the infection. We find that increasing herd size increases the duration of infectiousness and that the size of the herd plays a more significant role in determining this duration than the number of initially infected cattle in that herd. We also find that the assumptions made regarding the onset of infectiousness and clinical signs, which are based on contradictory empirical findings, can result in the predictions about when infection would be detectable differing by several days. Therefore, the disease progression used to characterize the course of infection in a single bovine host could have significant implications for determining when herds can be detected and subsequently controlled; the timing of which could influence the overall predicted trajectory of outbreaks.

Environmental transmission of influenza A virus in mallards.

IMPORTANCE: Wild birds are the natural reservoir hosts of influenza A viruses. Highly pathogenic strains of influenza A viruses pose risks to wild birds, poultry, and human health. Thus, understanding how these viruses are transmitted between birds is critical. We conducted an experiment where we experimentally infected mallards which are ducks that are commonly exposed to influenza viruses. We exposed several contact ducks to the experimentally infected duck to estimate the probability that a contact duck would become infected from either exposure to the virus shed directly from the infected duck or shared water contaminated with the virus from the infected duck. We found that environmental transmission from contaminated water best predicted the probability of transmission to naïve contact ducks, relatively low levels of virus in the water were sufficient to cause infection, and the probability of a naïve duck becoming infected varied over time.

Geo-temporal patterns to design cost-effective interventions for zoonotic diseases -the case of brucellosis in the country of Georgia.

Introduction: Control of zoonosis can benefit from geo-referenced procedures. Focusing on brucellosis, here the ability of two methods to distinguish disease dissemination patterns and promote cost-effective interventions was compared. Method: Geographical data on bovine, ovine and human brucellosis reported in the country of Georgia between 2014 and 2019 were investigated with (i) the Hot Spot (HS) analysis and (ii) a bio-geographical (BG) alternative. Results: More than one fourth of all sites reported cases affecting two or more species. While ruminant cases displayed different patterns over time, most human cases described similar geo-temporal features, which were associated with the route used by migrant shepherds. Other human cases showed heterogeneous patterns. The BG approach identified small areas with a case density twice as high as the HS method. The BG method also identified, in 2018, a 2.6-2.99 higher case density in zoonotic (human and non-human) sites than in non-zoonotic sites (which only reported cases affecting a single species) -a finding that, if corroborated, could support cost-effective policy-making. Discussion: Three dissemination hypotheses were supported by the data: (i) human cases induced by sheep-related contacts; (ii) human cases probably mediated by contaminated milk or meat; and (iii) cattle and sheep that infected one another. This proof-of-concept provided a preliminary validation for a method that may support cost-effective interventions oriented to control zoonoses. To expand these findings, additional studies on zoonosis-related decision-making are recommended.

Model-guided fieldwork: practical guidelines for multidisciplinary research on wildlife ecological and epidemiological dynamics.

Infectious disease ecology has recently raised its public profile beyond the scientific community due to the major threats that wildlife infections pose to biological conservation, animal welfare, human health and food security. As we start unravelling the full extent of emerging infectious diseases, there is an urgent need to facilitate multidisciplinary research in this area. Even though research in ecology has always had a strong theoretical component, cultural and technical hurdles often hamper direct collaboration between theoreticians and empiricists. Building upon our collective experience of multidisciplinary research and teaching in this area, we propose practical guidelines to help with effective integration among mathematical modelling, fieldwork and laboratory work. Modelling tools can be used at all steps of a field-based research programme, from the formulation of working hypotheses to field study design and data analysis. We illustrate our model-guided fieldwork framework with two case studies we have been conducting on wildlife infectious diseases: plague transmission in prairie dogs and lyssavirus dynamics in American and African bats. These demonstrate that mechanistic models, if properly integrated in research programmes, can provide a framework for holistic approaches to complex biological systems.

Modeling nation-wide U.S. swine movement networks at the resolution of the individual premises.

The spread of infectious livestock diseases is a major cause for concern in modern agricultural systems. In the dynamics of the transmission of such diseases, movements of livestock between herds play an important role. When constructing mathematical models used for activities such as forecasting epidemic development, evaluating mitigation strategies, or determining important targets for disease surveillance, including between-premises shipments is often a necessity. In the United States (U.S.), livestock shipment data is not routinely collected, and when it is, it is not readily available and mostly concerned with between-state shipments. To bridge this gap in knowledge and provide insight into the complete livestock shipment network structure, we have developed the U.S. Animal Movement Model (USAMM). Previously, USAMM has only existed for cattle shipments, but here we present a version for domestic swine. This new version of USAMM consists of a Bayesian model fit to premises demography, county-level livestock industry variables, and two limited data sets of between-state swine movements. The model scales up the data to simulate nation-wide networks of both within- and between-state shipments at the level of individual premises. Here we describe this shipment model in detail and subsequently explore its usefulness with a rudimentary predictive model of the prevalence of porcine epidemic diarrhea virus (PEDv) across the U.S. Additionally, in order to promote further research on livestock disease and other topics involving the movements of swine in the U.S., we also make 250 synthetic premises-level swine shipment networks with complete coverage of the entire conterminous U.S. freely available to the research community as a useful surrogate for the absent shipment data.