Assessing recovery of spectacled eiders using a Bayesian decision analysis.

Assessing species status and making classification decisions under the Endangered Species Act is a critical step towards effective species conservation. However, classification decisions are liable to two errors: i) failing to classify a species as threatened or endangered that should be classified (underprotection), or ii) classifying a species as threatened or endangered when it is not warranted (overprotection). Recent surveys indicate threatened spectacled eider populations are increasing in western Alaska, prompting the U.S. Fish and Wildlife Service to reconsider the federal listing status. There are multiple criteria set for assessing spectacled eider status, and here we focus on the abundance and decision analysis criteria. We estimated population metrics using state-space models for Alaskan breeding populations of spectacled eiders. We projected abundance over 50 years using posterior estimates of abundance and process variation to estimate the probability of quasi-extinction. The decision analysis maps the risk of quasi-extinction to the loss associated with making a misclassification error (i.e., underprotection) through a loss function. Our results indicate that the Yukon Kuskokwim Delta breeding population in western Alaska has met the recovery criteria but the Arctic Coastal Plain population in northern Alaska has not. The methods employed here provide an example of accounting for uncertainty and incorporating value judgements in such a way that the decision-makers may understand the risk of committing a misclassification error. Incorporating the abundance threshold and decision analysis in the reclassification criteria greatly increases the transparency and defensibility of the classification decision, a critical aspect for making effective decisions about species management and conservation.

Incomplete host immunity favors the evolution of virulence in an emergent pathogen.

Immune memory evolved to protect hosts from reinfection, but incomplete responses that allow future reinfection may inadvertently select for more-harmful pathogens. We present empirical and modeling evidence that incomplete immunity promotes the evolution of higher virulence in a natural host-pathogen system. We performed sequential infections of house finches with Mycoplasma gallisepticum strains of various levels of virulence. Virulent bacterial strains generated stronger host protection against reinfection than less virulent strains and thus excluded less virulent strains from infecting previously exposed hosts. In a two-strain model, the resulting fitness advantage selected for an almost twofold increase in pathogen virulence. Thus, the same immune systems that protect hosts from infection can concomitantly drive the evolution of more-harmful pathogens in nature.

Evolution of pathogen virulence across space during an epidemic.

We explore pathogen virulence evolution during the spatial expansion of an infectious disease epidemic in the presence of a novel host movement trade-off, using a simple, spatially explicit mathematical model. This work is motivated by empirical observations of the Mycoplasma gallisepticum invasion into North American house finch (Haemorhous mexicanus) populations; however, our results likely have important applications to other emerging infectious diseases in mobile hosts. We assume that infection reduces host movement and survival and that across pathogen strains the severity of these reductions increases with pathogen infectiousness. Assuming these trade-offs between pathogen virulence (host mortality), pathogen transmission, and host movement, we find that pathogen virulence levels near the epidemic front (that maximize wave speed) are lower than those that have a short-term growth rate advantage or that ultimately prevail (i.e., are evolutionarily stable) near the epicenter and where infection becomes endemic (i.e., that maximize the pathogen basic reproductive ratio). We predict that, under these trade-offs, less virulent pathogen strains will dominate the periphery of an epidemic and that more virulent strains will increase in frequency after invasion where disease is endemic. These results have important implications for observing and interpreting spatiotemporal epidemic data and may help explain transient virulence dynamics of emerging infectious diseases.

Parallel patterns of increased virulence in a recently emerged wildlife pathogen.

The evolution of higher virulence during disease emergence has been predicted by theoretical models, but empirical studies of short-term virulence evolution following pathogen emergence remain rare. Here we examine patterns of short-term virulence evolution using archived isolates of the bacterium Mycoplasma gallisepticum collected during sequential emergence events in two geographically distinct populations of the host, the North American house finch (Haemorhous [formerly Carpodacus] mexicanus). We present results from two complementary experiments, one that examines the trend in pathogen virulence in eastern North American isolates over the course of the eastern epidemic (1994-2008), and the other a parallel experiment on Pacific coast isolates of the pathogen collected after M. gallisepticum established itself in western North American house finch populations (2006-2010). Consistent with theoretical expectations regarding short-term or dynamic evolution of virulence, we show rapid increases in pathogen virulence on both coasts following the pathogen's establishment in each host population. We also find evidence for positive genetic covariation between virulence and pathogen load, a proxy for transmission potential, among isolates of M. gallisepticum. As predicted by theory, indirect selection for increased transmission likely drove the evolutionary increase in virulence in both geographic locations. Our results provide one of the first empirical examples of rapid changes in virulence following pathogen emergence, and both the detected pattern and mechanism of positive genetic covariation between virulence and pathogen load are consistent with theoretical expectations. Our study provides unique empirical insight into the dynamics of short-term virulence evolution that are likely to operate in other emerging pathogens of wildlife and humans.

Evolution of virulence in heterogeneous host communities under multiple trade-offs.

Many pathogens and parasites are transmitted through hosts that differ in species, sex, genotype, or immune status. In addition, virulence (here defined as disease-induced mortality) and transmission can vary during the infectious period within hosts of different state. Most models of virulence evolution assume that transmission and virulence are constant over the infectious period and that the host population is homogenous. Here, we examine a multispecies susceptible-infected-recovered (SIR) model where transmission occurs within and between species, and transmission and virulence varied during the infectious period. This allows us to understand virulence evolution in a broader range of situations that characterize many emerging diseases. Because emerging pathogens are by definition new to their host populations, they should be expected to rapidly adapt after emergence. We illustrate these evolutionary effects using the framework of adaptive dynamics to examine how virulence evolves after emergence in response to the relative strength of selection on pathogen fitness and mutational variance for virulence. We illustrate the role of evolution by simulating adaptive walks to an evolutionarily stable virulence. We found that the magnitude of between-species transmission and the relative timing of transmission and mortality across species were of primary importance for determining the evolutionarily stable virulence.

Pathogenicity and immunogenicity of three Mycoplasma gallisepticum isolates in house finches (Carpodacus mexicanus).

Mycoplasma gallisepticum (MG) has become a common cause of conjunctivitis in free-living house finches (Carpodacus mexicanus) since its emergence in the early 1990s. To date, temporal and spatial genotypic variation in MG has been documented, but phenotypic variation in pathogenicity and immunogenicity has not been examined. House finches were inoculated with MG isolates Virginia (VA)1994, California (CA)2006, or North Carolina (NC)2006, which were cultured from free-living house finches with conjunctivitis in 1994, 2006, and 2006, respectively. Infection with NC2006 resulted in the most severe eye lesions, highest pathogen loads, and highest levels of pathogen-specific lachrymal and serum antibodies. Infection with CA2006 caused the least severe eye lesions, lowest pathogen load, and lowest levels of antibodies. A small number of birds in each group developed protracted, severe disease in spite of robust antibody responses, suggesting that immunopathology may contribute to the lesions. Immunoblot analyses indicated that isolates are antigenically similar; thus, there may be partial cross-protection if a house finch encounters two or more strains of MG throughout the course of its lifetime. This study provides evidence that MG strains or strain variants circulating in house finch populations vary in their ability to cause disease, induce antibody responses, and persist in the host.

Evolution of virulence when transmission occurs before disease.

Most models of virulence evolution assume that transmission and virulence are constant during an infection. In many viral (HIV and influenza), bacterial (TB) and prion (BSE and CWD) systems, disease-induced mortality occurs long after the host becomes infectious. Therefore, we constructed a model with two infected classes that differ in transmission rate and virulence in order to understand how the evolutionarily stable strategy (ESS) depends on the relative difference in transmission and virulence between classes, on the transition rate between classes and on the recovery rate from the second class. We find that ESS virulence decreases when expressed early in the infection or when transmission occurs late in an infection. When virulence occurred relatively equally in each class and there was disease recovery, ESS virulence increased with increased transition rate. In contrast, ESS virulence first increased and then decreased with transition rate when there was little virulence early in the infection and a rapid recovery rate. This model predicts that ESS virulence is highly dependent on the timing of transmission and pathology after infection; thus, pathogen evolution may either increase or decrease virulence after emergence in a new host.

Spatial and temporal patterns of chronic wasting disease: fine-scale mapping of a wildlife epidemic in Wisconsin.

Emerging infectious diseases threaten wildlife populations and human health. Understanding the spatial distributions of these new diseases is important for disease management and policy makers; however, the data are complicated by heterogeneities across host classes, sampling variance, sampling biases, and the space-time epidemic process. Ignoring these issues can lead to false conclusions or obscure important patterns in the data, such as spatial variation in disease prevalence. Here, we applied hierarchical Bayesian disease mapping methods to account for risk factors and to estimate spatial and temporal patterns of infection by chronic wasting disease (CWD) in white-tailed deer (Odocoileus virginianus) of Wisconsin, U.S.A. We found significant heterogeneities for infection due to age, sex, and spatial location. Infection probability increased with age for all young deer, increased with age faster for young males, and then declined for some older animals, as expected from disease-associated mortality and age-related changes in infection risk. We found that disease prevalence was clustered in a central location, as expected under a simple spatial epidemic process where disease prevalence should increase with time and expand spatially. However, we could not detect any consistent temporal or spatiotemporal trends in CWD prevalence. Estimates of the temporal trend indicated that prevalence may have decreased or increased with nearly equal posterior probability, and the model without temporal or spatiotemporal effects was nearly equivalent to models with these effects based on deviance information criteria. For maximum interpretability of the role of location as a disease risk factor, we used the technique of direct standardization for prevalence mapping, which we develop and describe. These mapping results allow disease management actions to be employed with reference to the estimated spatial distribution of the disease and to those host classes most at risk. Future wildlife epidemiology studies should employ hierarchical Bayesian methods to smooth estimated quantities across space and time, account for heterogeneities, and then report disease rates based on an appropriate standardization.

Host culling as an adaptive management tool for chronic wasting disease in white-tailed deer: a modelling study.

Emerging wildlife diseases pose a significant threat to natural and human systems. Because of real or perceived risks of delayed actions, disease management strategies such as culling are often implemented before thorough scientific knowledge of disease dynamics is available. Adaptive management is a valuable approach in addressing the uncertainty and complexity associated with wildlife disease problems and can be facilitated by using a formal model.We developed a multi-state computer simulation model using age, sex, infection-stage, and seasonality as a tool for scientific learning and managing chronic wasting disease (CWD) in white-tailed deer Odocoileus virginianus. Our matrix model used disease transmission parameters based on data collected through disease management activities. We used this model to evaluate management issues on density- (DD) and frequency-dependent (FD) transmission, time since disease introduction, and deer culling on the demographics, epizootiology, and management of CWD.Both DD and FD models fit the Wisconsin data for a harvested white-tailed deer population, but FD was slightly better. Time since disease introduction was estimated as 36 (95% CI, 24-50) and 188 (41->200) years for DD and FD transmission, respectively. Deer harvest using intermediate to high non-selective rates can be used to reduce uncertainty between DD and FD transmission and improve our prediction of long-term epidemic patterns and host population impacts. A higher harvest rate allows earlier detection of these differences, but substantially reduces deer abundance.Results showed that CWD has spread slowly within Wisconsin deer populations, and therefore, epidemics and disease management are expected to last for decades. Non-hunted deer populations can develop and sustain a high level of infection, generating a substantial risk of disease spread. In contrast, CWD prevalence remains lower in hunted deer populations, but at a higher prevalence the disease competes with recreational hunting to reduce deer abundance.Synthesis and applications. Uncertainty about density- or frequency-dependent transmission hinders predictions about the long-term impacts of chronic wasting disease on cervid populations and the development of appropriate management strategies. An adaptive management strategy using computer modelling coupled with experimental management and monitoring can be used to test model predictions, identify the likely mode of disease transmission, and evaluate the risks of alternative management responses.

Immune response to sympatric and allopatric parasites in a snail-trematode interaction.

BACKGROUND: The outcome of parasite exposure depends on the (1) genetic specificity of the interaction, (2) induction of host defenses, and (3) parasite counter defenses. We studied both the genetic specificity for infection and the specificity for the host-defense response in a snail-trematode interaction (Potamopyrgus antipodarum-Microphallus sp.) by conducting a reciprocal cross-infection experiment between two populations of host and parasite. RESULTS: We found that infection was greater in sympatric host-parasite combinations. We also found that the host-defense response (hemocyte concentration) was induced by parasite exposure, but the response did not increase with increased parasite dose nor did it depend on parasite source, host source, or host-parasite combination. CONCLUSION: The results are consistent with a genetically specific host-parasite interaction, but inconsistent with a general arms-race type interaction where allocation to defense is the main determinant of host resistance.

Host sex and local adaptation by parasites in a snail-trematode interaction.

One of the leading theories for the evolutionary stability of sex in eukaryotes relies on parasite-mediated selection against locally common host genotypes (the Red Queen hypothesis). As such, parasites would be expected to be better at infecting sympatric host populations than allopatric host populations. Here we examined all published and unpublished infection experiments on a snail-trematode system (Potamopyrgus antipodarum and Microphallus sp., respectively). A meta-analysis demonstrated significant local adaptation by the parasite, and a variance components analysis showed that the variance due to the host-parasite interaction far exceeded the variance due to the main effects of host source and parasite source. The meta-analysis also indicated that asexual host populations were more resistant to allopatric sources of parasites than were (mostly) sexual host populations, but we found no significant differences among parasite populations in the strength of local adaptation. This result suggests that triploid asexual snails are more resistant to remote sources of parasites, but the parasite has, through coevolution, overcome the difference. Finally, we found that the degree of local adaptation did not depend on the genetic distance among host populations. Taken together, the results demonstrate that the parasites are adapted, on average, to infecting their local host populations and suggest that they may be a factor in selecting against common host genotypes in natural populations.