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.

Linking mosquito surveillance to dengue fever through Bayesian mechanistic modeling.

Our ability to effectively prevent the transmission of the dengue virus through targeted control of its vector, Aedes aegypti, depends critically on our understanding of the link between mosquito abundance and human disease risk. Mosquito and clinical surveillance data are widely collected, but linking them requires a modeling framework that accounts for the complex non-linear mechanisms involved in transmission. Most critical are the bottleneck in transmission imposed by mosquito lifespan relative to the virus' extrinsic incubation period, and the dynamics of human immunity. We developed a differential equation model of dengue transmission and embedded it in a Bayesian hierarchical framework that allowed us to estimate latent time series of mosquito demographic rates from mosquito trap counts and dengue case reports from the city of Vitória, Brazil. We used the fitted model to explore how the timing of a pulse of adult mosquito control influences its effect on the human disease burden in the following year. We found that control was generally more effective when implemented in periods of relatively low mosquito mortality (when mosquito abundance was also generally low). In particular, control implemented in early September (week 34 of the year) produced the largest reduction in predicted human case reports over the following year. This highlights the potential long-term utility of broad, off-peak-season mosquito control in addition to existing, locally targeted within-season efforts. Further, uncertainty in the effectiveness of control interventions was driven largely by posterior variation in the average mosquito mortality rate (closely tied to total mosquito abundance) with lower mosquito mortality generating systems more vulnerable to control. Broadly, these correlations suggest that mosquito control is most effective in situations in which transmission is already limited by mosquito abundance.

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.

Projected wetland densities under climate change: habitat loss but little geographic shift in conservation strategy.

Climate change poses major challenges for conservation and management because it alters the area, quality, and spatial distribution of habitat for natural populations. To assess species' vulnerability to climate change and target ongoing conservation investments, researchers and managers often consider the effects of projected changes in climate and land use on future habitat availability and quality and the uncertainty associated with these projections. Here, we draw on tools from hydrology and climate science to project the impact of climate change on the density of wetlands in the Prairie Pothole Region of the USA, a critical area for breeding waterfowl and other wetland-dependent species. We evaluate the potential for a trade-off in the value of conservation investments under current and future climatic conditions and consider the joint effects of climate and land use. We use an integrated set of hydrological and climatological projections that provide physically based measures of water balance under historical and projected future climatic conditions. In addition, we use historical projections derived from ten general circulation models (GCMs) as a baseline from which to assess climate change impacts, rather than historical climate data. This method isolates the impact of greenhouse gas emissions and ensures that modeling errors are incorporated into the baseline rather than attributed to climate change. Our work shows that, on average, densities of wetlands (here defined as wetland basins holding water) are projected to decline across the U.S. Prairie Pothole Region, but that GCMs differ in both the magnitude and the direction of projected impacts. However, we found little evidence for a shift in the locations expected to provide the highest wetland densities under current vs. projected climatic conditions. This result was robust to the inclusion of projected changes in land use under climate change. We suggest that targeting conservation towards wetland complexes containing both small and relatively large wetland basins, which is an ongoing conservation strategy, may also act to hedge against uncertainty in the effects of climate change.

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.

Bayesian Modeling of Prion Disease Dynamics in Mule Deer Using Population Monitoring and Capture-Recapture Data.

Epidemics of chronic wasting disease (CWD) of North American Cervidae have potential to harm ecosystems and economies. We studied a migratory population of mule deer (Odocoileus hemionus) affected by CWD for at least three decades using a Bayesian framework to integrate matrix population and disease models with long-term monitoring data and detailed process-level studies. We hypothesized CWD prevalence would be stable or increase between two observation periods during the late 1990s and after 2010, with higher CWD prevalence making deer population decline more likely. The weight of evidence suggested a reduction in the CWD outbreak over time, perhaps in response to intervening harvest-mediated population reductions. Disease effects on deer population growth under current conditions were subtle with a 72% chance that CWD depressed population growth. With CWD, we forecasted a growth rate near one and largely stable deer population. Disease effects appear to be moderated by timing of infection, prolonged disease course, and locally variable infection. Long-term outcomes will depend heavily on whether current conditions hold and high prevalence remains a localized phenomenon.

AGE AND REPEATED BIOPSY INFLUENCE ANTEMORTEM PRP(CWD) TESTING IN MULE DEER (ODOCOILEUS HEMIONUS) IN COLORADO, USA.

Biopsy of rectal mucosa-associated lymphoid tissue provides a useful, but imperfect, live-animal test for chronic wasting disease (CWD) in mule deer (Odocoileus hemionus). It is difficult and expensive to complete these tests on free-ranging animals, and wildlife health managers will benefit from methods that can accommodate test results of varying quality. To this end, we developed a hierarchical Bayesian model to estimate the probability that an individual is infected based on test results. Our model was estimated with the use of data on 210 adult female mule deer repeatedly tested during 2010-14. The ability to identify infected individuals correctly declined with age and may have been influenced by repeated biopsy. Fewer isolated lymphoid follicles (where PrP(CWD) accumulates) were obtained in biopsies of older deer and the proportion of follicles showing PrP(CWD) was reduced. A deer's genotype in the prion gene (PRNP) also influenced detection. At least five follicles were needed in a biopsy to assure a 95% accurate test in PRNP genotype 225SS deer.

Anticipating the prevalence of avian influenza subtypes H9 and H5 in live-bird markets.

An ability to forecast the prevalence of specific subtypes of avian influenza viruses (AIV) in live-bird markets would facilitate greatly the implementation of preventative measures designed to minimize poultry losses and human exposure. The minimum requirement for developing predictive quantitative tools is surveillance data of AIV prevalence sampled frequently over several years. Recently, a 4-year time series of monthly sampling of hemagglutinin subtypes 1-13 in ducks, chickens and quail in live-bird markets in southern China has become available. We used these data to investigate whether a simple statistical model, based solely on historical data (variables such as the number of positive samples in host X of subtype Y time t months ago), could accurately predict prevalence of H5 and H9 subtypes in chickens. We also examined the role of ducks and quail in predicting prevalence in chickens within the market setting because between-species transmission is thought to occur within markets but has not been measured. Our best statistical models performed remarkably well at predicting future prevalence (pseudo-R(2) = 0.57 for H9 and 0.49 for H5), especially considering the multi-host, multi-subtype nature of AIVs. We did not find prevalence of H5/H9 in ducks or quail to be predictors of prevalence in chickens within the Chinese markets. Our results suggest surveillance protocols that could enable more accurate and timely predictive statistical models. We also discuss which data should be collected to allow the development of mechanistic models.

Bayesian multimodel inference for geostatistical regression models.

The problem of simultaneous covariate selection and parameter inference for spatial regression models is considered. Previous research has shown that failure to take spatial correlation into account can influence the outcome of standard model selection methods. A Markov chain Monte Carlo (MCMC) method is investigated for the calculation of parameter estimates and posterior model probabilities for spatial regression models. The method can accommodate normal and non-normal response data and a large number of covariates. Thus the method is very flexible and can be used to fit spatial linear models, spatial linear mixed models, and spatial generalized linear mixed models (GLMMs). The Bayesian MCMC method also allows a priori unequal weighting of covariates, which is not possible with many model selection methods such as Akaike's information criterion (AIC). The proposed method is demonstrated on two data sets. The first is the whiptail lizard data set which has been previously analyzed by other researchers investigating model selection methods. Our results confirmed the previous analysis suggesting that sandy soil and ant abundance were strongly associated with lizard abundance. The second data set concerned pollution tolerant fish abundance in relation to several environmental factors. Results indicate that abundance is positively related to Strahler stream order and a habitat quality index. Abundance is negatively related to percent watershed disturbance.

Spatial variability of western bean cutworm (Lepidoptera: Noctuidae) pheromone trap captures in sprinkler irrigated corn in eastern Colorado.

Strategies for controlling pests are an integral part of any agricultural management plan. Most field crops, such as wheat (Triticum spp.) and corn (Zea mays L.) are managed as if they are homogeneous units. However, pests within fields are rarely homogenous. Development of plans that use targeted pest control tactics requires knowledge of the ecological drivers of the pest species. That is, by understanding the spatio-temporal factors influencing pest populations, we can develop management strategy to prevent or reduce pest damage. This study was conducted to quantify variables influencing the spatial variability of adult male western bean cutworm, Striacosta albicosta (Smith). Striacosta albicosta moths were collected in pheromone traps in two center pivot, irrigated corn fields near Wiggins, CO. We hypothesized that moth abundance would be influenced by the distance from the edge of the field, distance to nearest alternative corn crop and affected by anisotropic effects, such as prevailing wind direction. Greater trap catches of S. albicosta in each of the fields were found with increased proximity to the edge of the field, if the nearest neighboring crop was corn. Prevailing wind direction and directional effects were found to influence abundance. Results serve as a first step toward building a precision pest management system for controlling S. albicosta.

A structured and dynamic framework to advance traits-based theory and prediction in ecology.

Predicting changes in community composition and ecosystem function in a rapidly changing world is a major research challenge in ecology. Traits-based approaches have elicited much recent interest, yet individual studies are not advancing a more general, predictive ecology. Significant progress will be facilitated by adopting a coherent theoretical framework comprised of three elements: an underlying trait distribution, a performance filter defining the fitness of traits in different environments, and a dynamic projection of the performance filter along some environmental gradient. This framework allows changes in the trait distribution and associated modifications to community composition or ecosystem function to be predicted across time or space. The structure and dynamics of the performance filter specify two key criteria by which we judge appropriate quantitative methods for testing traits-based hypotheses. Bayesian multilevel models, dynamical systems models and hybrid approaches meet both these criteria and have the potential to meaningfully advance traits-based ecology.

The role of geographic information systems in wildlife landscape epidemiology: models of chronic wasting disease in Colorado mule deer.

The authors present findings from two landscape epidemiology studies of chronic wasting disease (CWD) in northern Colorado mule deer (Odocoileus hemionus). First, the effects of human land use on disease prevalence were explored by formulating a set of models estimating CWD prevalence in relation to differences in human land use, sex and geographic location. Prevalence was higher in developed areas and among male deer suggesting that anthropogenic influences (changes in land use), differences in exposure risk between sexes and landscape-scaled heterogeneity are associated with CWD prevalence. The second study focused on identifying scales of mule deer movement and mixing that had the greatest influence on the spatial pattern of CWD in north-central Colorado. The authors hypothesised that three scales of mixing - individual, winter subpopulation and summer subpopulation - might control spatial variation in disease prevalence. A fully Bayesian hierarchical model was developed to compare the strength of evidence for each mixing scale. Strong evidence was found indicating that the finest mixing scale corresponded best to the observed spatial distribution of CWD prevalence. This analysis demonstrates how information on the scales of spatial processes that generate observed patterns can be used to gain insight into the epidemiology of wildlife diseases when process data are sparse or unavailable.

Linking chronic wasting disease to mule deer movement scales: a hierarchical Bayesian approach.

Observed spatial patterns in natural systems may result from processes acting across multiple spatial and temporal scales. Although spatially explicit data on processes that generate ecological patterns, such as the distribution of disease over a landscape, are frequently unavailable, information about the scales over which processes operate can be used to understand the link between pattern and process. Our goal was to identify scales of mule deer (Odocoileus hemionus) movement and mixing that exerted the greatest influence on the spatial pattern of chronic wasting disease (CWD) in northcentral Colorado, USA. We hypothesized that three scales of mixing (individual, winter subpopulation, or summer subpopulation) might control spatial variation in disease prevalence. We developed a fully Bayesian hierarchical model to compare the strength of evidence for each mixing scale. We found strong evidence that the finest mixing scale corresponded best to the spatial distribution of CWD infection. There was also evidence that land ownership and habitat use play a role in exacerbating the disease, along with the known effects of sex and age. Our analysis demonstrates how information on the scales of spatial processes that generate observed patterns can be used to gain insight when process data are sparse or unavailable.

Model selection for geostatistical models.

We consider the problem of model selection for geospatial data. Spatial correlation is often ignored in the selection of explanatory variables, and this can influence model selection results. For example, the importance of particular explanatory variables may not be apparent when spatial correlation is ignored. To address this problem, we consider the Akaike Information Criterion (AIC) as applied to a geostatistical model. We offer a heuristic derivation of the AIC in this context and provide simulation results that show that using AIC for a geostatistical model is superior to the often-used traditional approach of ignoring spatial correlation in the selection of explanatory variables. These ideas are further demonstrated via a model for lizard abundance. We also apply the principle of minimum description length (MDL) to variable selection for the geostatistical model. The effect of sampling design on the selection of explanatory covariates is also explored. R software to implement the geostatistical model selection methods described in this paper is available in the Supplement.

Autoregressive models for capture-recapture data: a Bayesian approach.

In this article, we incorporate an autoregressive time-series framework into models for animal survival using capture-recapture data. Researchers modeling animal survival probabilities as the realization of a random process have typically considered survival to be independent from one time period to the next. This may not be realistic for some populations. Using a Gibbs sampling approach, we can estimate covariate coefficients and autoregressive parameters for survival models. The procedure is illustrated with a waterfowl band recovery dataset for northern pintails (Anas acuta). The analysis shows that the second lag autoregressive coefficient is significantly less than 0, suggesting that there is a triennial relationship between survival probabilities and emphasizing that modeling survival rates as independent random variables may be unrealistic in some cases. Software to implement the methodology is available at no charge on the Internet.