Genome-wide association study of an unusual dolphin mortality event reveals candidate genes for susceptibility and resistance to cetacean morbillivirus

Genomic Divergence and the Evolution of Ecotypes in Bottlenose Dolphins (Genus Tursiops)

Climatic changes have caused major environmental restructuring throughout the world's oceans. Marine organisms have responded to novel conditions through various biological systems, including genomic adaptation. Growing accessibility of next-generation DNA sequencing methods to study nonmodel species has recently allowed genomic changes underlying environmental adaptations to be investigated. This study used double-digest restriction-site associated DNA (ddRAD) sequence data to investigate the genomic basis of ecotype formation across currently recognized species and subspecies of bottlenose dolphins (genus Tursiops) in the Southern Hemisphere. Subspecies-level genomic divergence was confirmed between the offshore common bottlenose dolphin (T. truncatus truncatus) and the inshore Lahille's bottlenose dolphin (T. t. gephyreus) from the southwestern Atlantic Ocean (SWAO). Similarly, subspecies-level divergence is suggested between inshore (eastern Australia) Indo-Pacific bottlenose dolphin (T. aduncus) and the proposed Burrunan dolphin (T. australis) from southern Australia. Inshore bottlenose dolphin lineages generally had lower genomic diversity than offshore lineages, a pattern particularly evident for T. t. gephyreus, which showed exceptionally low diversity. Genomic regions associated with cardiovascular, musculoskeletal, and energy production systems appear to have undergone repeated adaptive evolution in inshore lineages across the Southern Hemisphere. We hypothesize that comparable selective pressures in the inshore environment drove similar adaptive responses in each lineage, supporting parallel evolution of inshore bottlenose dolphins. With climate change altering marine ecosystems worldwide, it is crucial to gain an understanding of the adaptive capacity of local species and populations. Our study provides insights into key adaptive pathways that may be important for the long-term survival of cetaceans and other organisms in a changing marine environment.

A targeted approach to investigating immune genes of an iconic Australian marsupial

Disease is a contributing factor to the decline of wildlife populations across the globe. Koalas, iconic yet declining Australian marsupials, are predominantly impacted by two pathogens, Chlamydia and koala retrovirus. Chlamydia is an obligate intracellular bacterium and one of the most widespread sexually transmitted infections in humans worldwide. In koalas, Chlamydia infections can present as asymptomatic or can cause a range of ocular and urogenital disease signs, such as conjunctivitis, cystitis and infertility. In this study, we looked at differences in response to Chlamydia in two northern populations of koalas using a targeted gene sequencing of 1209 immune genes in addition to genome‐wide reduced representation data. We identified two MHC Class I genes associated with Chlamydia disease progression as well as 25 single nucleotide polymorphisms across 17 genes that were associated with resolution of Chlamydia infection. These genes are involved in the innate immune response (TLR5) and defence (TLR5, IFNγ, SERPINE1, STAT2 and STX4). This study deepens our understanding of the role that genetics plays in disease progression in koalas and leads into future work that will use whole genome resequencing of a larger sample set to investigate in greater detail regions identified in this study. Elucidation of the role of host genetics in disease progression and resolution in koalas will directly contribute to better design of Chlamydia vaccines and management of koala populations which have recently been listed as “endangered.”

Riverscape Genomics Clarifies Neutral and Adaptive Evolution in an Amazonian Characin Fish (Triportheus albus)

Understanding the role of natural selection in the evolution of wild populations is challenging due to the spatial complexity of natural systems. The richest diversity of freshwater fishes in the world is found in the Amazon Basin, a system where marked hydrochemical differences exist at the interface of major rivers with distinct “water colors” (i.e., black, white, and clear water). We hypothesize that divergent natural selection associated with these “aquatic ecotones” influences population-level adaptive divergence in the non-migratory Amazonian fish fauna. This hypothesis was tested using a landscape genomics framework to compare the relative contribution of environmental and spatial factors to the evolutionary divergence of the Amazonian characin fish Triportheus albus. The framework was based on spatial data, in situ hydrochemical measurements, and 15,251 filtered SNPs (single nucleotide polymorphisms) for T. albus sampled from three major Amazonian rivers. Gradient Forest, redundancy analysis (RDA) and BayPass analyses were used to test for signals of natural selection, and model-based and model-free approaches were used to evaluate neutral population differentiation. After controlling for a signal of neutral hierarchical structure which was consistent with the expectations for a dendritic system, variation in turbidity and pH were key factors contributing to adaptive divergence. Variation in genes involved in acid-sensitive ion transport pathways and light-sensitive photoreceptor pathways was strongly associated with pH and turbidity variability. This study improves our understanding of how natural selection and neutral evolution impact on the distribution of aquatic biodiversity from the understudied and ecologically complex Amazonia.

Seascape genomics of coastal bottlenose dolphins along strong gradients of temperature and salinity

Heterogeneous seascapes and strong environmental gradients in coastal waters are expected to influence adaptive divergence, particularly in species with large population sizes where selection is expected to be highly efficient. However, these influences might also extend to species characterized by strong social structure, natal philopatry and small home ranges. We implemented a seascape genomic study to test this hypothesis in Indo-Pacific bottlenose dolphins (Tursiops aduncus) distributed along the environmentally heterogeneous coast of southern Australia. The datasets included oceanographic and environmental variables thought to be good predictors of local adaptation in dolphins and 8,081 filtered single nucleotide polymorphisms (SNPs) genotyped for individuals sampled from seven different bioregions. From a neutral perspective, population structure and connectivity of the dolphins were generally influenced by habitat type and social structuring. Genotype-environment association analysis identified 241 candidate adaptive loci and revealed that sea surface temperature and salinity gradients influenced adaptive divergence in these animals at both large- (1,000s km) and fine-scales (<100 km). Enrichment analysis and annotation of candidate genes revealed functions related to sodium-activated ion transport, kidney development, adipogenesis and thermogenesis. The findings of spatial adaptive divergence and inferences of putative physiological adaptations challenge previous suggestions that marine megafauna is most likely to be affected by environmental and climatic changes via indirect, trophic effects. Our work contributes to conservation management of coastal bottlenose dolphins subjected to anthropogenic disturbance and to efforts of clarifying how seascape heterogeneity influences adaptive diversity and evolution in small cetaceans.

Molecular signatures in cetacean morbillivirus and host species proteomes: Unveiling the evolutionary dynamics of an enigmatic pathogen?

Cetacean morbillivirus (CeMV) infects marine mammals often causing a fatal respiratory and neurological disease. Recently, CeMV has expanded its geographic and host species range, with cases being reported worldwide among dolphins, whales, seals, and other aquatic mammalian species, and therefore has emerged as the most threatening nonanthropogenic factor affecting marine mammal's health and conservation. Extensive research efforts have aimed to understand CeMV epidemiology and ecology, however, the molecular mechanisms underlying its transmission and pathogenesis are still poorly understood. In particular, the field suffers from a knowledge gap on the structural and functional properties of CeMV proteins and their host interactors. Nevertheless, the body of scientific literature produced in recent years has inaugurated new investigational trends, driving future directions in CeMV molecular research. In this mini-review, the most recent literature has been summarized in the context of such research trends, and categorized into four priority research topics, such as (1) the interaction between CeMV glycoprotein and its host cell receptors across several species; (2) the CeMV molecular determinants responsible for different disease phenotype; (3) the host molecular determinants responsible for differential susceptibility to CeMV infection; (4) the CeMV molecular determinants responsible for difference virulence among circulating CeMV strains. Arguably, these are the most urgent topics that need to be investigated and that most promisingly will help to shed light on the details of CeMV evolutionary dynamics in the immediate future.

Whole genomes reveal multiple candidate genes and pathways involved in the immune response of dolphins to a highly infectious virus

Wildlife species are challenged by various infectious diseases that act as important demographic drivers of populations and have become a great conservation concern particularly under growing environmental changes. The new era of whole genome sequencing provides new opportunities and avenues to explore the role of genetic variants in the plasticity of immune responses, particularly in non-model systems. Cetacean morbillivirus (CeMV) has emerged as a major viral threat to cetacean populations worldwide, contributing to the death of thousands of individuals of multiple dolphin and whale species. To understand the genomic basis of immune responses to CeMV, we generated and analysed whole genomes of 53 Indo-Pacific bottlenose dolphins (Tursiops aduncus) exposed to Australia's largest known CeMV-related mortality event that killed at least 50 dolphins from three different species. The genomic data set consisted of 10,168,981 SNPs anchored onto 23 chromosome-length scaffolds and 77 short scaffolds. Whole genome analysis indicated that levels of inbreeding in the dolphin population did not influence the outcome of an individual. Allele frequency estimates between survivors and nonsurvivors of the outbreak revealed 15,769 candidate SNPs, of which 689 were annotated to 295 protein coding genes. These included 50 genes with functions related to innate and adaptive immune responses, and cytokine signalling pathways and genes thought to be involved in immune responses to other morbilliviruses. Our study characterised genomic regions and pathways that may contribute to CeMV immune responses in dolphins. This represents a stride towards clarifying the complex interactions of the cetacean immune system and emphasises the value of whole genome data sets in understanding genetic elements that are essential for species conservation, including disease susceptibility and adaptation.

Genomic association with pathogen carriage in bighorn sheep (Ovis canadensis)

Genetic composition can influence host susceptibility to, and transmission of, pathogens, with potential population-level consequences. In bighorn sheep (Ovis canadensis), pneumonia epidemics caused by Mycoplasma ovipneumoniae have been associated with severe population declines and limited recovery across North America. Adult survivors either clear the infection or act as carriers that continually shed M. ovipneumoniae and expose their susceptible offspring, resulting in high rates of lamb mortality for years following the outbreak event. Here, we investigated the influence of genomic composition on persistent carriage of M. ovipneumoniae in a well-studied bighorn sheep herd in the Wallowa Mountains of Oregon, USA. Using 10,605 SNPs generated using RADseq technology for 25 female bighorn sheep, we assessed genomic diversity metrics and employed family-based genome-wide association methodologies to understand variant association and genetic architecture underlying chronic carriage. We observed no differences among genome-wide diversity metrics (heterozygosity and allelic richness) between groups. However, we identified two variant loci of interest and seven associated candidate genes, which may influence carriage status. Further, we found that the SNP panel explained ~55% of the phenotypic variance (SNP-based heritability) for M. ovipneumoniae carriage, though there was considerable uncertainty in these estimates. While small sample sizes limit conclusions drawn here, our study represents one of the first to assess the genomic factors influencing chronic carriage of a pathogen in a wild population and lays a foundation for understanding genomic influence on pathogen persistence in bighorn sheep and other wildlife populations. Future research should incorporate additional individuals as well as distinct herds to further explore the genomic basis of chronic carriage.

Sarcoptic mange severity is associated with reduced genomic variation and evidence of selection in Yellowstone National Park wolves (Canis lupus)

Population genetic theory posits that molecular variation buffers against disease risk. Although this "monoculture effect" is well supported in agricultural settings, its applicability to wildlife populations remains in question. In the present study, we examined the genomics underlying individual-level disease severity and population-level consequences of sarcoptic mange infection in a wild population of canids. Using gray wolves (Canis lupus) reintroduced to Yellowstone National Park (YNP) as our focal system, we leveraged 25 years of observational data and biobanked blood and tissue to genotype 76,859 loci in over 400 wolves. At the individual level, we reported an inverse relationship between host genomic variation and infection severity. We additionally identified 410 loci significantly associated with mange severity, with annotations related to inflammation, immunity, and skin barrier integrity and disorders. We contextualized results within environmental, demographic, and behavioral variables, and confirmed that genetic variation was predictive of infection severity. At the population level, we reported decreased genome-wide variation since the initial gray wolf reintroduction event and identified evidence of selection acting against alleles associated with mange infection severity. We concluded that genomic variation plays an important role in disease severity in YNP wolves. This role scales from individual to population levels, and includes patterns of genome-wide variation in support of the monoculture effect and specific loci associated with the complex mange phenotype. Results yielded system-specific insights, while also highlighting the relevance of genomic analyses to wildlife disease ecology, evolution, and conservation.