TERRESTRIAL PATHOGEN POLLUTANT, TOXOPLASMA GONDII, THREATENS HAWAIIAN MONK SEALS (NEOMONACHUS SCHAUINSLANDI) FOLLOWING HEAVY RUNOFF EVENTS.

Toxoplasmosis is a major threat to Hawaiian monk seals (Neomonachus schauinslandi) in the main Hawaiian Islands where seal habitat overlaps with substantial human and domestic cat populations. As the definitive hosts, members of the Felidae are the sole sources contaminating the environment with infectious oocysts; these oocysts can be transported into the marine environment, thereby threatening marine mammals. To understand environmental factors influencing Hawaiian monk seal exposure to Toxoplasma gondii, we examined monk seal strandings from toxoplasmosis in relationship to location and rainfall patterns throughout the main Hawaiian Islands. Using a case-control study design, we compared mortalities due to toxoplasmosis (cases) with those from other causes (controls). We found that cases were up to 25 times more likely than controls to occur after heavy runoff events. The greatest odds ratio was observed when rainfall occurred 3 wk before strandings, potentially indicating important timelines in the disease process. Our results suggest that heavy rainfall frequently delivers sufficient numbers of oocysts to infect Hawaiian monk seals. With infectious doses of as low as a single oocyst, any contaminated runoff constitutes a risk to Hawaii's endangered monk seal.

ESTABLISHING HEMATOLOGY AND SERUM CHEMISTRY REFERENCE INTERVALS FOR WILD HAWAIIAN MONK SEALS ( NEOMONACHUS SCHAUINSLANDI).

Hematology and serum chemistry reference intervals have been previously established for the endangered Hawaiian monk seal ( Neomonachus schauinslandi) as an imperative measure for health assessments. Monitoring the health of the wild population depends upon reference intervals that are context specific; hence we developed reference intervals from fresh samples, as opposed to frozen, from wild monk seals. This study builds on the number of parameters from previous efforts by using samples collected between 2004 and 2015 from wild monk seals. Blood samples were analyzed by a single veterinary diagnostic laboratory within 24 hr of collection from apparently healthy, wild seals during research activities. Reference intervals were determined based on the analytical steps outlined by the American Society for Veterinary Clinical Pathology. These comprehensive hematology and serum chemistry reference intervals enable more consistent and systematic interpretation of results, which will guide individual and population-level health assessment and decision-making research and recovery activities.

Survey for Placental Disease and Reproductive Pathogens in the Endangered Hawaiian Monk Seal ( Neomonachus schauinslandi).

There is considerable temporal and spatial variability in the reproductive rates of Hawaiian monk seals (HMS; Neomonachus schauinslandi). Poor reproductive performance limits the recovery of this endangered species; however, causal factors are not fully understood. There is serologic evidence that HMS are exposed to pathogens that can impact reproductive success, but the prevalence of placental infections in HMS has not been evaluated. Placental tissues ( n=50), including tissues from 25% of known HMS births, were opportunistically collected in 2011 from six Northwestern Hawaiian Islands and three main Hawaiian Islands. Reproductive histories of the sampled females were representative of the breeding population, as determined through comparisons in age of primiparity and mature reproductive rate. Placental tissues were examined histologically and screened by PCR for Coxiella burnetii, Brucella spp., Chlamydia spp., Leptospira spp., herpesviruses, and Toxoplasma gondii. There was no histologic evidence of placental pathology, and molecular analyses were negative. These negative results can be used to estimate pathogen prevalence in the nonsampled population. For an approximate population size of 1,300 HMS, we can estimate with 99% confidence that the prevalence of each pathogen tested is 9% or less. This is low relative to other pinnipeds and indicates that factors other than reproductive pathology, such as resource limitation, may drive variability in HMS reproductive rates. Further investigation into the cumulative impacts of resource limitation and other stressors on HMS reproduction is warranted.

Model recommendations meet management reality: implementation and evaluation of a network-informed vaccination effort for endangered Hawaiian monk seals.

Where disease threatens endangered wildlife populations, substantial resources are required for management actions such as vaccination. While network models provide a promising tool for identifying key spreaders and prioritizing efforts to maximize efficiency, population-scale vaccination remains rare, providing few opportunities to evaluate performance of model-informed strategies under realistic scenarios. Because the endangered Hawaiian monk seal could be heavily impacted by disease threats such as morbillivirus, we implemented a prophylactic vaccination programme. We used contact networks to prioritize vaccinating animals with high contact rates. We used dynamic network models to simulate morbillivirus outbreaks under real and idealized vaccination scenarios. We then evaluated the efficacy of model recommendations in this real-world vaccination project. We found that deviating from the model recommendations decreased the efficiency; requiring 44% more vaccinations to achieve a given decrease in outbreak size. However, we gained protection more quickly by vaccinating available animals rather than waiting to encounter priority seals. This work demonstrates the value of network models, but also makes trade-offs clear. If vaccines were limited but time was ample, vaccinating only priority animals would maximize herd protection. However, where time is the limiting factor, vaccinating additional lower-priority animals could more quickly protect the population.

MODELING A MORBILLIVIRUS OUTBREAK IN HAWAIIAN MONK SEALS (NEOMONACHUS SCHAUINSLANDI) TO AID IN THE DESIGN OF MITIGATION PROGRAMS.

We developed a stochastic susceptible-exposed-infectious-removed (SEIR) model to simulate a range of plausible morbillivirus outbreak scenarios in a randomly mixing population of 170 endangered Hawaiian monk seals (Neomonachus schauinslandi). We then modeled realistic vaccination and quarantine measures to determine the potential efficacy of such mitigation efforts. Morbillivirus outbreaks represent substantial risk to monk seals-91% of simulated baseline outbreaks grew (R0>1), and in one-third of the scenarios all, or nearly all, individuals were infected. Simulated vaccination efforts in response to an outbreak were not effective in substantially reducing infections, largely because of the prolonged interval between vaccination and immunity. Prophylactic vaccination, in contrast, could be an effective tool for preventing outbreaks. Herd immunity is practically achievable because of the small sizes of monk seal populations and the animals' accessibility on shore. Adding realistic spatial structure to the model, as informed by movement of seals tracked in the main Hawaiian Islands with the use of telemetry, greatly reduced the simulated impact of outbreaks (≤10 seals were infected in 62% of spatially structured simulations). Although response vaccination remained relatively ineffective, spatial segregation allowed herd immunity to be achieved through prophylactic vaccination with less effort. In a randomly mixing population of 170 seals, 86% would need to be vaccinated to achieve herd immunity in 95% of simulated outbreaks, compared to only approximately 60% in three spatially segregated subgroups with the same combined abundance. Simulations indicate that quarantining a modest number (up to 20) of ill seals has the potential to extinguish even fast-growing outbreaks rapidly. The efficacy of quarantine, however, is highly dependent upon rapid detection and response. We conclude that prophylactic vaccination combined with a quarantine program supported by vigilant surveillance and rapid, reliable diagnosis could greatly mitigate the threat of a morbillivirus outbreak in Hawaiian monk seals.

ESTIMATING CONTACT RATES OF HAWAIIAN MONK SEALS (NEOMONACHUS SCHAUINSLANDI) USING SOCIAL NETWORK ANALYSIS.

Understanding disease transmission dynamics, which are in part mediated by rates and patterns of social contact, is fundamental to predicting the likelihood, rate of spread, impacts, and mitigation of disease outbreaks in wildlife populations. Contact rates, which are important parameters required for epidemiologic models, are difficult to estimate. The endangered Hawaiian monk seal (Neomonachus schauinslandi) may be particularly vulnerable to morbillivirus outbreaks, due to its low abundance, lack of genetic diversity, and history of isolation from mammalian diseases. Morbillivirus epizootics have had devastating effects on other seal populations. We constructed social networks based on visual observations of individually identifiable monk seals associating onshore to estimate contact rates, assuming random mixing, and also to investigate contact patterns of different age and sex classes. Contact rates estimated from two island populations in 4 yr were remarkably similar, indicating any two individuals have about a one in 1,000 chance of making contact on any given day. Further, contact patterns within and among age and sex classes were statistically different from random. The methods we used could be broadly applied to empirically derive contact rates using association data. These rates are critical for epidemiologic modelling to simulate wildlife disease outbreaks and to inform science-based prevention and mitigation programs.

Disease risk in a dynamic environment: the spread of tick-borne pathogens in Minnesota, USA.

As humans and climate change alter the landscape, novel disease risk scenarios emerge. Understanding the complexities of pathogen emergence and subsequent spread as shaped by landscape heterogeneity is crucial to understanding disease emergence, pinpointing high-risk areas, and mitigating emerging disease threats in a dynamic environment. Tick-borne diseases present an important public health concern and incidence of many of these diseases are increasing in the United States. The complex epidemiology of tick-borne diseases includes strong ties with environmental factors that influence host availability, vector abundance, and pathogen transmission. Here, we used 16 years of case data from the Minnesota Department of Health to report spatial and temporal trends in Lyme disease (LD), human anaplasmosis, and babesiosis. We then used a spatial regression framework to evaluate the impact of landscape and climate factors on the spread of LD. Finally, we use the fitted model, and landscape and climate datasets projected under varying climate change scenarios, to predict future changes in tick-borne pathogen risk. Both forested habitat and temperature were important drivers of LD spread in Minnesota. Dramatic changes in future temperature regimes and forest communities predict rising risk of tick-borne disease.

Evaluating spatial overlap and relatedness of white-tailed deer in a chronic wasting disease management zone.

Wildlife disease transmission, at a local scale, can occur from interactions between infected and susceptible conspecifics or from a contaminated environment. Thus, the degree of spatial overlap and rate of contact among deer is likely to impact both direct and indirect transmission of infectious diseases such chronic wasting disease (CWD) or bovine tuberculosis. We identified a strong relationship between degree of spatial overlap (volume of intersection) and genetic relatedness for female white-tailed deer in Wisconsin's area of highest CWD prevalence. We used volume of intersection as a surrogate for contact rates between deer and concluded that related deer are more likely to have contact, which may drive disease transmission dynamics. In addition, we found that age of deer influences overlap, with fawns exhibiting the highest degree of overlap with other deer. Our results further support the finding that female social groups have higher contact among related deer which can result in transmission of infectious diseases. We suggest that control of large social groups comprised of closely related deer may be an effective strategy in slowing the transmission of infectious pathogens, and CWD in particular.

The walk is never random: subtle landscape effects shape gene flow in a continuous white-tailed deer population in the Midwestern United States.

One of the pervasive challenges in landscape genetics is detecting gene flow patterns within continuous populations of highly mobile wildlife. Understanding population genetic structure within a continuous population can give insights into social structure, movement across the landscape and contact between populations, which influence ecological interactions, reproductive dynamics or pathogen transmission. We investigated the genetic structure of a large population of deer spanning the area of Wisconsin and Illinois, USA, affected by chronic wasting disease. We combined multiscale investigation, landscape genetic techniques and spatial statistical modelling to address the complex questions of landscape factors influencing population structure. We sampled over 2000 deer and used spatial autocorrelation and a spatial principal components analysis to describe the population genetic structure. We evaluated landscape effects on this pattern using a spatial autoregressive model within a model selection framework to test alternative hypotheses about gene flow. We found high levels of genetic connectivity, with gradients of variation across the large continuous population of white-tailed deer. At the fine scale, spatial clustering of related animals was correlated with the amount and arrangement of forested habitat. At the broader scale, impediments to dispersal were important to shaping genetic connectivity within the population. We found significant barrier effects of individual state and interstate highways and rivers. Our results offer an important understanding of deer biology and movement that will help inform the management of this species in an area where overabundance and disease spread are primary concerns.

Emerging prion disease drives host selection in a wildlife population.

Infectious diseases are increasingly recognized as an important force driving population dynamics, conservation biology, and natural selection in wildlife populations. Infectious agents have been implicated in the decline of small or endangered populations and may act to constrain population size, distribution, growth rates, or migration patterns. Further, diseases may provide selective pressures that shape the genetic diversity of populations or species. Thus, understanding disease dynamics and selective pressures from pathogens is crucial to understanding population processes, managing wildlife diseases, and conserving biological diversity. There is ample evidence that variation in the prion protein gene (PRNP) impacts host susceptibility to prion diseases. Still, little is known about how genetic differences might influence natural selection within wildlife populations. Here we link genetic variation with differential susceptibility of white-tailed deer to chronic wasting disease (CWD), with implications for fitness and disease-driven genetic selection. We developed a single nucleotide polymorphism (SNP) assay to efficiently genotype deer at the locus of interest (in the 96th codon of the PRNP gene). Then, using a Bayesian modeling approach, we found that the more susceptible genotype had over four times greater risk of CWD infection; and, once infected, deer with the resistant genotype survived 49% longer (8.25 more months). We used these epidemiological parameters in a multi-stage population matrix model to evaluate relative fitness based on genotype-specific population growth rates. The differences in disease infection and mortality rates allowed genetically resistant deer to achieve higher population growth and obtain a long-term fitness advantage, which translated into a selection coefficient of over 1% favoring the CWD-resistant genotype. This selective pressure suggests that the resistant allele could become dominant in the population within an evolutionarily short time frame. Our work provides a rare example of a quantifiable disease-driven selection process in a wildlife population, demonstrating the potential for infectious diseases to alter host populations. This will have direct bearing on the epidemiology, dynamics, and future trends in CWD transmission and spread. Understanding genotype-specific epidemiology will improve predictive models and inform management strategies for CWD-affected cervid populations.

The role of genetics in chronic wasting disease of North American cervids.

Chronic wasting disease (CWD) is a major concern for the management of North American cervid populations. This fatal prion disease has led to declines in populations which have high CWD prevalence and areas with both high and low infection rates have experienced economic losses in wildlife recreation and fears of potential spill-over into livestock or humans. Research from human and veterinary medicine has established that the prion protein gene (Prnp) encodes the protein responsible for transmissible spongiform encephalopathies (TSEs). Polymorphisms in the Prnp gene can lead to different prion forms that moderate individual susceptibility to and progression of TSE infection. Prnp genes have been sequenced in a number of cervid species including those currently infected by CWD (elk, mule deer, white-tailed deer, moose) and those for which susceptibility is not yet determined (caribou, fallow deer, sika deer). Over thousands of sequences examined, the Prnp gene is remarkably conserved within the family Cervidae; only 16 amino acid polymorphisms have been reported within the 256 amino acid open reading frame in the third exon of the Prnp gene. Some of these polymorphisms have been associated with lower rates of CWD infection and slower progression of clinical CWD. Here we review the body of research on Prnp genetics of North American cervids. Specifically, we focus on known polymorphisms in the Prnp gene, observed genotypic differences in CWD infection rates and clinical progression, mechanisms for genetic TSE resistance related to both the cervid host and the prion agent and potential for natural selection for CWD-resistance. We also identify gaps in our knowledge that require future research.

Diversity and distribution of white-tailed deer mtDNA lineages in chronic wasting disease (CWD) outbreak areas in southern Wisconsin, USA.

Chronic wasting disease (CWD) is a transmissible spongiform encephalopathy affecting North American cervids. Because it is uniformly fatal, the disease is a major concern in the management of white-tailed deer populations. Management programs to control CWD require improved knowledge of deer interaction, movement, and population connectivity that could influence disease transmission and spread. Genetic methods were employed to evaluate connectivity among populations in the CWD management zone of southern Wisconsin. A 576-base-pair region of the mitochondrial DNA of 359 white-tailed deer from 12 sample populations was analyzed. Fifty-eight variable sites were detected within the sequence, defining 43 haplotypes. While most sample populations displayed similar levels of haplotype diversity, individual haplotypes were clustered on the landscape. Spatial clusters of different haplotypes were apparent in distinct ecoregions surrounding CWD outbreak areas. The spatial distribution of mtDNA haplotypes suggests that clustering of the deer matrilineal groups and population connectivity are associated with broad-scale geographic landscape features. These landscape characteristics may also influence the contact rates between groups and therefore the potential spread of CWD; this may be especially true of local disease spread between female social groups. Our results suggest that optimal CWD management needs to be tailored to fit gender-specific dispersal behaviors and regional differences in deer population connectivity. This information will help wildlife managers design surveillance and monitoring efforts based on population interactions and potential deer movement among CWD-affected and unaffected areas.

Oral rabies vaccine (ORV) bait uptake by captive striped skunks.

Aerial delivery of oral rabies vaccine (ORV) baits has proven effective in large-scale efforts to immunize wildlife against rabies, and in North America this strategy currently is being used to immunize foxes (Urocyon cinereoargenteus and Vulpes vulpes), raccoons (Procyon lotor), and coyotes (Canis latrans). Skunks are also a major reservoir and vector of rabies, but at present oral vaccines for use in skunks are not licensed. Furthermore, given differences in morphology (smaller jaws) and behavior (food handling and consumption), it is unknown if baits currently used in ORV campaigns would be effective for skunks. Because oral vaccine delivery is contingent upon puncture of the vaccine container (VC), baits need to be sufficiently attractive to elicit selection and consumption. Manipulation of the bait to facilitate vaccine ingestion by the target species is a critical element for an effective ORV bait. The objectives of this study were to assess manipulation and consumption of current ORV baits by striped skunks (Mephitis mephitis). We conducted four independent trials with penned animals and various baits to assess bait selection frequency, VC puncture frequency, and consumption. Video recorded trials were used to assess attractiveness of baits and consumption behavior of skunks. Bait characteristics, such as texture, size, and flavor influenced selection and consumption. Fish and chicken flavors were preferred and vaccine containers within selected baits were likely to be punctured. Vaccine ingestion seemed more likely if VCs were directly coated with the bait matrix. To make baits attractive to skunks and to ensure puncture of the VC, modifications to current baits should consider a smaller size, a meat-flavored matrix, a slightly pressurized VC, and a direct coating of matrix on the VC.