Demographic History and Adaptive Evolution of Indo-Pacific Bottlenose Dolphins (Tursiops aduncus) in Western Australia

Demographic processes can substantially affect a species' response to changing ecological conditions, necessitating the combined consideration of genetic responses to environmental variables and neutral genetic variation. Using a seascape genomics approach combined with population demographic modelling, we explored the interplay of demographic and environmental factors that shaped the current population structure in Indo-Pacific bottlenose dolphins (Tursiops aduncus) along the Western Australian coastline. We combined large-scale environmental data gathered via remote sensing with RADseq genomic data from 133 individuals at 19 sampling sites. Using population genetic and outlier detection analyses, we identified three distinct genetic clusters, coinciding with tropical, subtropical and temperate provincial bioregions. In contrast to previous studies, our demographic models indicated that populations occupying the paleo-shoreline split into two demographically independent lineages before the last glacial maximum (LGM). A subsequent split after the LGM 12-15 kya gave rise to the Shark Bay population, thereby creating the three currently observed clusters. Although multi-locus heterozygosity declined from north to south, dolphins from the southernmost cluster inhabiting temperate waters had higher heterozygosity in potentially adaptive loci compared to dolphins from subtropical and tropical waters. These findings suggest ongoing adaptation to cold-temperate waters in the southernmost cluster, possibly linked to distinct selective pressures between the different bioregions. Our study demonstrated that in the marine realm, without apparent physical boundaries, only a combined approach can fully elucidate the intricate environmental and genetic interactions shaping the evolutionary trajectory of marine mammals.

Inferring the origin of new D-loop haplotypes of loggerhead sea turtles (Testudinata: Cheloniidae) from the Southwest Atlantic lineage

The populations of the loggerhead turtles, Caretta caretta, present four main D-loop mitochondrial haplogroups that are distributed across the Indo-Pacific, Mediterranean, and Atlantic oceans. The Southwestern Atlantic (SWA) is one of the Regional Management Units (RMUs) of loggerheads, characterized by unique haplotypes, high nest density, and distinct life history traits. Detecting new D-loop haplogroups is important, particularly endemic ones, as they can enhance our understanding of their life history within the RMUs and contribute to the resolution of mixed stock analysis. In this study, we conducted a series of phylogenetic delimitation and network analyses to identify, validate, and infer the origin of four new D-loop haplotypes detected in the loggerhead populations from the SWA. Our findings demonstrate that these new D-loop haplotypes are valid and unique to the SWA lineage, potentially aiding in the delimitation of individuals' origins and the inference of their lineage.

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.

Mitochondrial genomes and genetic structure of the Kemp's ridley sea turtle (Lepidochelys kempii)

The Kemp's ridley (Lepidochelys kempii) is the world's most endangered sea turtle species and is primarily distributed in the Gulf of Mexico. In the United States, South Padre Island, Texas serves as a key nesting ground for the species. Genetic studies of the Kemp's ridley have been used to aid in conservation and management practices, with the mitochondrial control region as the most commonly used marker due to its perceived hypervariability and ease of sequencing. However, with the advent of next generation sequencing technology, targeting complete mitochondrial genomes is now feasible. Here, we describe a more complete mitochondrial genome for the Kemp's ridley than has been previously published in literature and demonstrate a cost-effective and efficient method for obtaining complete mitochondrial genomes from sea turtles. We compare the genetic diversity and taxonomic resolution obtained from whole mitochondrial genomes to that obtained from the mitochondrial control region alone. We compare current genetic diversity with previous records. Furthermore, we evaluate the genetic structure between the breeding stock in South Padre Island and that of deceased Kemp's ridleys recovered on the Northern coast of the Gulf of Mexico after the 2010 BP Deepwater Horizon oil spill, and of Kemp's ridleys stranded on the East Coast of the United States. Our results show that complete mitochondrial genomes provide greater resolution than the control region alone. They also show that the genetic diversity of the Kemp's ridley has remained stable, despite large population declines, and that the genetic makeup of deceased turtles stranded after the Deepwater Horizon oil spill is indistinguishable from the breeding stock in South Padre Island, Texas.

OPEN DATA BADGE

This article has earned an Open Data Badge for making publicly available the digitally-shareable data necessary to reproduce the reported results. The data is available at https://www.ncbi.nlm.nih.gov/genbank/.

Intrapopulation variability in the timing of ontogenetic habitat shifts in sea turtles revealed using δ15 N values from bone growth rings

Determining location and timing of ontogenetic shifts in the habitat use of highly migratory species, along with possible intrapopulation variation in these shifts, is essential for understanding mechanisms driving alternate life histories and assessing overall population trends. Measuring variations in multi-year habitat-use patterns is especially difficult for remote oceanic species. To investigate the potential for differential habitat use among migratory marine vertebrates, we measured the naturally occurring stable nitrogen isotope (δ¹⁵ N) patterns that differentiate distinct ocean regions to create a 'regional isotope characterization', analysed the δ¹⁵ N values from annual bone growth layer rings from dead-stranded animals, and then combined the bone and regional isotope data to track individual animal movement patterns over multiple years. We used humeri from juvenile North Pacific loggerhead turtles (Caretta caretta), animals that undergo long migrations across the North Pacific Ocean (NPO), using multiple discrete regions as they develop to adulthood. Typical of many migratory marine species, ontogenetic changes in habitat use throughout their decades-long juvenile stage is poorly understood, but each potential habitat has unique foraging opportunities and spatially explicit natural and anthropogenic threats that could affect key life-history parameters. We found a bimodal size/age distribution in the timing that juveniles underwent an ontogenetic habitat shift from the oceanic central North Pacific (CNP) to the neritic east Pacific region near the Baja California Peninsula (BCP) (42·7 ± 7·2 vs. 68·3 ± 3·4 cm carapace length, 7·5 ± 2·7 vs. 15·6 ± 1·7 years). Important to the survival of this population, these disparate habitats differ considerably in their food availability, energy requirements and threats, and these differences can influence life-history parameters such as growth, survival and future fecundity. This is the first evidence of alternative ontogenetic shifts and habitat-use patterns for juveniles foraging in the eastern NPO. We combine two techniques, skeletochronology and stable isotope analysis, to reconstruct multi-year habitat-use patterns of a remote migratory species, linked to estimated ages and body sizes of individuals, to reveal variable ontogeny during the juvenile life stage that could drive alternate life histories and that has the potential to illuminate the migration patterns for other species with accretionary tissues.