A model for leveraging animal movement to understand spatio-temporal disease dynamics.

The ongoing explosion of fine-resolution movement data in animal systems provides a unique opportunity to empirically quantify spatial, temporal and individual variation in transmission risk and improve our ability to forecast disease outbreaks. However, we lack a generalizable model that can leverage movement data to quantify transmission risk and how it affects pathogen invasion and persistence on heterogeneous landscapes. We developed a flexible model 'Movement-driven modelling of spatio-temporal infection risk' (MoveSTIR) that leverages diverse data on animal movement to derive metrics of direct and indirect contact by decomposing transmission into constituent processes of contact formation and duration and pathogen deposition and acquisition. We use MoveSTIR to demonstrate that ignoring fine-scale animal movements on actual landscapes can mis-characterize transmission risk and epidemiological dynamics. MoveSTIR unifies previous work on epidemiological contact networks and can address applied and theoretical questions at the nexus of movement and disease ecology.

Perceived Stress and Molecular Bacterial Vaginosis in the National Institutes of Health Longitudinal Study of Vaginal Flora.

Vaginal microbiota provide the first line of defense against urogenital infections primarily through protective actions of Lactobacillus species Perceived stress increases susceptibility to infection through several mechanisms, including suppression of immune function. We investigated whether stress was associated with deleterious changes to vaginal bacterial composition in a subsample of 572 women in the Longitudinal Study of Vaginal Flora, sampled from 1999 through 2002. Using Cox proportional hazards models, both unadjusted and adjusted for sociodemographic factors and sexual behaviors, we found that participants who exhibited a 5-unit-increase in Cohen's Perceived Stress Scale had greater risk (adjusted hazard ratio (HR) = 1.40, 95% confidence interval (CI): 1.13, 1.74) of developing molecular bacterial vaginosis (BV), a state with low Lactobacillus abundance and diverse anaerobic bacteria. A 5-unit increase in stress score was also associated with greater risks of transitioning from the L. iners-dominated community state type (26% higher) to molecular-BV (adjusted HR = 1.26, 95% CI: 1.01, 1.56) or maintaining molecular-BV from baseline (adjusted HR = 1.23, 95% CI: 1.01, 1.47). Inversely, women with baseline molecular-BV reporting a 5-unit stress increase were less likely to transition to microbiota dominated by L. crispatus, L. gasseri, or L. jensenii (adjusted HR = 0.81, 95% CI: 0.68, 0.99). These findings suggest that psychosocial stress is associated with vaginal microbiota composition, inviting a more mechanistic exploration of the relationship between psychosocial stress and molecular-BV.

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.

Effects of social structure and management on risk of disease establishment in wild pigs.

Contact heterogeneity among hosts determines invasion and spreading dynamics of infectious disease, thus its characterization is essential for identifying effective disease control strategies. Yet, little is known about the factors shaping contact networks in many wildlife species and how wildlife management actions might affect contact networks. Wild pigs in North America are an invasive, socially structured species that pose a health concern for domestic swine given their ability to transmit numerous devastating diseases such as African swine fever (ASF). Using proximity loggers and GPS data from 48 wild pigs in Florida and South Carolina, USA, we employed a probabilistic framework to estimate weighted contact networks. We determined the effects of sex, social group and spatial distribution (monthly home-range overlap and distance) on wild pig contact. We also estimated the impacts of management-induced perturbations on contact and inferred their effects on ASF establishment in wild pigs with simulation. Social group membership was the primary factor influencing contacts. Between-group contacts depended primarily on space use characteristics, with fewer contacts among groups separated by >2 km and no contacts among groups >4 km apart within a month. Modelling ASF dynamics on the contact network demonstrated that indirect contacts resulting from baiting (a typical method of attracting wild pigs or game species to a site to enhance recreational hunting) increased the risk of disease establishment by ~33% relative to direct contact. Low-intensity population reduction (<5.9% of the population) had no detectable impact on contact structure but reduced predicted ASF establishment risk relative to no population reduction. We demonstrate an approach for understanding the relative role of spatial, social and individual-level characteristics in shaping contact networks and predicting their effects on disease establishment risk, thus providing insight for optimizing disease control in spatially and socially structured wildlife species.

Mixed ancestry from wild and domestic lineages contributes to the rapid expansion of invasive feral swine.

Invasive alien species are a significant threat to both economic and ecological systems. Identifying the processes that give rise to invasive populations is essential for implementing effective control strategies. We conducted an ancestry analysis of invasive feral swine (Sus scrofa, Linnaeus, 1758), a highly destructive ungulate that is widely distributed throughout the contiguous United States, to describe introduction pathways, sources of newly emergent populations and processes contributing to an ongoing invasion. Comparisons of high-density single nucleotide polymorphism genotypes for 6,566 invasive feral swine to a comprehensive reference set of S. scrofa revealed that the vast majority of feral swine were of mixed ancestry, with dominant genetic associations to Western heritage breeds of domestic pig and European populations of wild boar. Further, the rapid expansion of invasive feral swine over the past 30 years was attributable to secondary introductions from established populations of admixed ancestry as opposed to direct introductions of domestic breeds or wild boar. Spatially widespread genetic associations of invasive feral swine to European wild boar deviated strongly from historical S. scrofa introduction pressure, which was largely restricted to domestic pigs with infrequent, localized wild boar releases. The deviation between historical introduction pressure and contemporary genetic ancestry suggests wild boar-hybridization may contribute to differential fitness in the environment and heightened invasive potential for individuals of admixed domestic pig-wild boar ancestry.

Optimal spatial prioritization of control resources for elimination of invasive species under demographic uncertainty.

Populations of invasive species often spread heterogeneously across a landscape, consisting of local populations that cluster in space but are connected by dispersal. A fundamental dilemma for invasive species control is how to optimally allocate limited fiscal resources across local populations. Theoretical work based on perfect knowledge of demographic connectivity suggests that targeting local populations from which migrants originate (sources) can be optimal. However, demographic processes such as abundance and dispersal can be highly uncertain, and the relationship between local population density and damage costs (damage function) is rarely known. We used a metapopulation model to understand how budget and uncertainty in abundance, connectivity, and the damage function, together impact return on investment (ROI) for optimal control strategies. Budget, observational uncertainty, and the damage function had strong effects on the optimal resource allocation strategy. Uncertainty in dispersal probability was the least important determinant of ROI. The damage function determined which resource prioritization strategy was optimal when connectivity was symmetric but not when it was asymmetric. When connectivity was asymmetric, prioritizing source populations had a higher ROI than allocating effort equally across local populations, regardless of the damage function, but uncertainty in connectivity structure and abundance reduced ROI of the optimal prioritization strategy by 57% on average depending on the control budget. With low budgets (monthly removal rate of 6.7% of population), there was little advantage to prioritizing resources, especially when connectivity was high or symmetric, and observational uncertainty had only minor effects on ROI. Allotting funding for improved monitoring appeared to be most important when budgets were moderate (monthly removal of 13-20% of the population). Our result showed that multiple sources of observational uncertainty should be considered concurrently for optimizing ROI. Accurate estimates of connectivity direction and abundance were more important than accurate estimates of dispersal rates. Developing cost-effective surveillance methods to reduce observational uncertainties, and quantitative frameworks for determining how resources should be spatially apportioned to multiple monitoring and control activities are important and challenging future directions for optimizing ROI for invasive species control programs.

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.

"One Health" or Three? Publication Silos Among the One Health Disciplines.

The One Health initiative is a global effort fostering interdisciplinary collaborations to address challenges in human, animal, and environmental health. While One Health has received considerable press, its benefits remain unclear because its effects have not been quantitatively described. We systematically surveyed the published literature and used social network analysis to measure interdisciplinarity in One Health studies constructing dynamic pathogen transmission models. The number of publications fulfilling our search criteria increased by 14.6% per year, which is faster than growth rates for life sciences as a whole and for most biology subdisciplines. Surveyed publications clustered into three communities: one used by ecologists, one used by veterinarians, and a third diverse-authorship community used by population biologists, mathematicians, epidemiologists, and experts in human health. Overlap between these communities increased through time in terms of author number, diversity of co-author affiliations, and diversity of citations. However, communities continue to differ in the systems studied, questions asked, and methods employed. While the infectious disease research community has made significant progress toward integrating its participating disciplines, some segregation--especially along the veterinary/ecological research interface--remains.

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.

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.

Inferring invasive species abundance using removal data from management actions.

Evaluation of the progress of management programs for invasive species is crucial for demonstrating impacts to stakeholders and strategic planning of resource allocation. Estimates of abundance before and after management activities can serve as a useful metric of population management programs. However, many methods of estimating population size are too labor intensive and costly to implement, posing restrictive levels of burden on operational programs. Removal models are a reliable method for estimating abundance before and after management using data from the removal activities exclusively, thus requiring no work in addition to management. We developed a Bayesian hierarchical model to estimate abundance from removal data accounting for varying levels of effort, and used simulations to assess the conditions under which reliable population estimates are obtained. We applied this model to estimate site-specific abundance of an invasive species, feral swine (Sus scrofa), using removal data from aerial gunning in 59 site/time-frame combinations (480-19,600 acres) throughout Oklahoma and Texas, USA. Simulations showed that abundance estimates were generally accurate when effective removal rates (removal rate accounting for total effort) were above 0.40. However, when abundances were small (<50) the effective removal rate needed to accurately estimates abundances was considerably higher (0.70). Based on our post-validation method, 78% of our site/time frame estimates were accurate. To use this modeling framework it is important to have multiple removals (more than three) within a time frame during which demographic changes are minimized (i.e., a closed population; ≤3 months for feral swine). Our results show that the probability of accurately estimating abundance from this model improves with increased sampling effort (8+ flight hours across the 3-month window is best) and increased removal rate. Based on the inverse relationship between inaccurate abundances and inaccurate removal rates, we suggest auxiliary information that could be collected and included in the model as covariates (e.g., habitat effects, differences between pilots) to improve accuracy of removal rates and hence abundance estimates.

African swine fever at the wildlife-livestock interface: challenges for management and outbreak response within invasive wild pigs in the United States.

African swine fever (ASF) causes significant morbidity and mortality in both domestic and wild suids (Sus scrofa), and disease outbreaks convey profound economic costs to impacted industries due to death loss, the cost of culling exposed/infected animals as the primary disease control measure, and trade restrictions. The co-occurrence of domestic and wild suids significantly complicates ASF management given the potential for wild populations to serve as persistent sources for spillover. We describe the unique threat of African swine fever virus (ASFV) introduction to the United States from epidemiological and ecological perspectives with a specific focus on disease management at the wild-domestic swine interface. The introduction of ASF into domestic herds would require a response focused on containment, culling, and contact tracing. However, detection of ASF among invasive wild pigs would require a far more complex and intensive response given the challenges of detection, containment, and ultimately elimination among wild populations. We describe the state of the science available to inform preparations for an ASF response among invasive wild pigs, describe knowledge gaps and the associated studies needed to fill those gaps, and call for an integrated approach for preparedness that incorporates the best available science and acknowledges sociological attributes and the policy context needed for an integrated disease response.

Deriving spatially explicit direct and indirect interaction networks from animal movement data.

Quantifying spatiotemporally explicit interactions within animal populations facilitates the understanding of social structure and its relationship with ecological processes. Data from animal tracking technologies (Global Positioning Systems ["GPS"]) can circumvent longstanding challenges in the estimation of spatiotemporally explicit interactions, but the discrete nature and coarse temporal resolution of data mean that ephemeral interactions that occur between consecutive GPS locations go undetected. Here, we developed a method to quantify individual and spatial patterns of interaction using continuous-time movement models (CTMMs) fit to GPS tracking data. We first applied CTMMs to infer the full movement trajectories at an arbitrarily fine temporal scale before estimating interactions, thus allowing inference of interactions occurring between observed GPS locations. Our framework then infers indirect interactions-individuals occurring at the same location, but at different times-while allowing the identification of indirect interactions to vary with ecological context based on CTMM outputs. We assessed the performance of our new method using simulations and illustrated its implementation by deriving disease-relevant interaction networks for two behaviorally differentiated species, wild pigs (Sus scrofa) that can host African Swine Fever and mule deer (Odocoileus hemionus) that can host chronic wasting disease. Simulations showed that interactions derived from observed GPS data can be substantially underestimated when temporal resolution of movement data exceeds 30-min intervals. Empirical application suggested that underestimation occurred in both interaction rates and their spatial distributions. CTMM-Interaction method, which can introduce uncertainties, recovered majority of true interactions. Our method leverages advances in movement ecology to quantify fine-scale spatiotemporal interactions between individuals from lower temporal resolution GPS data. It can be leveraged to infer dynamic social networks, transmission potential in disease systems, consumer-resource interactions, information sharing, and beyond. The method also sets the stage for future predictive models linking observed spatiotemporal interaction patterns to environmental drivers.

Habitat quality influences trade-offs in animal movement along the exploration-exploitation continuum.

To successfully establish itself in a novel environment, an animal must make an inherent trade-off between knowledge accumulation and exploitation of knowledge gained (i.e., the exploration-exploitation dilemma). To evaluate how habitat quality affects the spatio-temporal scale of switching between exploration and exploitation during home range establishment, we conducted experimental trials comparing resource selection and space-use of translocated animals to those of reference individuals using reciprocal translocations between habitat types of differing quality. We selected wild pigs (Sus scrofa) as a model species to investigate hypotheses related to the movement behavior of translocated individuals because they are globally distributed large mammals that are often translocated within their introduced range to facilitate recreational hunting. Individuals translocated to higher quality habitat (i.e. higher proportions of bottomland hardwood habitats) exhibited smaller exploratory movements and began exploiting resources more quickly than those introduced to lower quality areas, although those in lower-quality areas demonstrated an increased rate of selection for preferred habitat as they gained knowledge of the landscape. Our data demonstrate that habitat quality mediates the spatial and temporal scale at which animals respond behaviorally to novel environments, and how these processes may determine the success of population establishment.

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.

Role of social structure in establishment of an invasive large mammal after translocation.

BACKGROUND: Data on the movement behavior of translocated wild pigs is needed to develop appropriate response strategies for containing and eliminating new source populations following translocation events. We conducted experimental trials to compare the home range establishment and space-use metrics, including the number of days and distance traveled before becoming range residents, for wild pigs translocated with their social group and individually. RESULTS: We found wild pigs translocated with their social group made less extensive movements away from the release location and established a stable home range ~5 days faster than those translocated individually. We also examined how habitat quality impacted the home range sizes of translocated wild pigs and found wild pigs maintained larger ranges in areas with higher proportion of low-quality habitat. CONCLUSION: Collectively, our findings suggest translocations of invasive wild pigs have a greater probability of establishing a viable population near the release site when habitat quality is high and when released with members of their social unit compared to individuals moved independent of their social group or to low-quality habitat. However, all wild pigs translocated in our study made extensive movements from their release location, highlighting the potential for single translocation events of either individuals or groups to have far-reaching consequences within a much broader landscape beyond the location where they are released. These results highlight the challenges associated with containing populations in areas where illegal introduction of wild pigs occurs, and the need for rapid response once releases are identified. © 2023 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Individual-level patterns of resource selection do not predict hotspots of contact.

Contact among animals is crucial for various ecological processes, including social behaviors, disease transmission, and predator-prey interactions. However, the distribution of contact events across time and space is heterogeneous, influenced by environmental factors and biological purposes. Previous studies have assumed that areas with abundant resources and preferred habitats attract more individuals and, therefore, lead to more contact. To examine the accuracy of this assumption, we used a use-available framework to compare landscape factors influencing the location of contacts between wild pigs (Sus scrofa) in two study areas in Florida and Texas (USA) from those influencing non-contact space use. We employed a contact-resource selection function (RSF) model, where contact locations were defined as used points and locations without contact as available points. By comparing outputs from this contact RSF with a general, population-level RSF, we assessed the factors driving both habitat selection and contact. We found that the landscape predictors (e.g., wetland, linear features, and food resources) played different roles in habitat selection from contact processes for wild pigs in both study areas. This indicated that pigs interacted with their landscapes differently when choosing habitats compared to when they encountered other individuals. Consequently, relying solely on the spatial overlap of individual or population-level RSF models may lead to a misleading understanding of contact-related ecology. Our findings challenge prevailing assumptions about contact and introduce innovative approaches to better understand the ecological drivers of spatially explicit contact. By accurately predicting the spatial distribution of contact events, we can enhance our understanding of contact based ecological processes and their spatial dynamics.

Genomic regions associated with pseudorabies virus infection status in naturally infected feral swine (Sus scrofa).

Pseudorabies virus (PRV)-the causative agent of Aujeszky's disease-was eliminated from commercial pig production herds in the United States (US) in 2004; however, PRV remains endemic among invasive feral swine (Sus scrofa). The circulation of PRV among abundant, widespread feral swine populations poses a sustained risk for disease spillover to production herds. Risk-based surveillance has been successfully implemented for PRV in feral swine populations in the US. However, understanding the role of host genetics in infection status may offer new insights into the epidemiology and disease dynamics of PRV that can be applied to management strategies. Genetic mechanisms underlying host susceptibility to PRV are relatively unknown; therefore, we sought to identify genomic regions associated with PRV infection status among naturally infected feral swine using genome-wide association studies (GWAS) and gene set enrichment analysis of single nucleotide polymorphism data (GSEA-SNP). Paired serological and genotypic data were collected from 6,081 feral swine distributed across the invaded range within the contiguous US. Three complementary study populations were developed for GWAS: 1) comprehensive population consisting of feral swine throughout the invaded range within the contiguous US; 2) population of feral swine under high, but temporally variable PRV infection pressure; and 3) population of feral swine under temporally stable, high PRV infection pressure. We identified one intronic SNP associated with PRV infection status within candidate gene AKAP6 on autosome 7. Various gene sets linked to metabolic pathways were enriched in the GSEA-SNP. Ultimately, improving disease surveillance efforts in feral swine will be critical to further understanding of the role host genetics play in PRV infection status, helping secure the health of commercial pork production.