Whole genome sequencing reveals variation in immune and nervous systems associated with harbor seal susceptibility to phocine distemper virus

Anthropogenic pressures have increased the emergence and spread of infectious disease in wildlife in recent decades, with cascading impacts on ecosystem dynamics and elevated risk to public health. Understanding factors that contribute to disease severity, and population- and community-level impacts, is critical to predict and mitigate infectious disease outbreaks. We aimed to investigate how host genetic factors contribute to harbor seal susceptibility to phocine distemper virus (PDV), a virus that has caused variable levels of mortality among pinniped populations across the North Atlantic. We conducted whole-genome sequencing of harbor seals that were stranded during and following a PDV outbreak that occurred in the Northeast United States in 2018 and used genome-wide association and FST-based approaches to identify genetic variants putatively associated with survival. A subset of harbor seals (10 cases, 10 controls) was sequenced at high coverage (~58x) to generate a reference haplotype panel and impute genotypes in remaining samples (37 cases, 45 controls) that were sequenced at a lower coverage (~4x). Upon investigation of 148 644 filtered single nucleotide polymorphisms (SNPs), we identified 33 candidate genes related to innate and adaptive immunity and 37 genes related to the nervous system that were associated with harbor seal survival. These results provide further evidence for the importance of immunogenetic variation and highlight additional key pathways and biological processes that may confer resistance to PDV. Genome-wide approaches, which can reveal the complexity of genetic variation underpinning disease-associated traits, are critical tools to help us understand and potentially predict the impact of future disease outbreaks.