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

Rolling with the punches: Organism-environment interactions shape spatial pattern of adaptive differentiation in the widespread mantis shrimp Oratosquilla oratoria

Investigating spatial pattern of adaptive variation and its underlying processes can inform the adaptive potential distributed within species ranges, which is increasingly important in the context of a changing climate. A correct interpretation of adaptive variation pattern requires that population history and the ensuing population genetic structure are taken into account. Here we carried out such a study by integrating population genomic analyses, demographic model testing and species distribution modeling to investigate patterns and causes of adaptive differentiation in a widespread mantis shrimp, Oratosquilla oratoria, along a replicated, broad-scale temperature gradient in the northwestern Pacific (NWP). Our results supported a strong hierarchical ecogeographic structure dominated by habitat-linked divergence among O. oratoria populations accompanied with introgressive hybridization. A combined FST outlier and environmental correlation analyses revealed remarkable temperature-associated clines in allele frequency across paired North-South populations on Chinese and Japanese coasts, and identified a suite of loci associated with temperature adaptation. Further demographic model testing revealed the observed clinal variation derived partly from Pleistocene divergence followed by recent secondary contact. More importantly, the likelihood of hybridization is predicted to increase as climate change progresses, which would break barriers to gene flow and enable the spread of adaptive genetic variation. These results support that not only is temperature-driven adaptive differentiation occurs in O. oratoria but that such pattern is likely attributed to ancient adaptive variation, sustained by contemporary ocean conditions and a semi-permeable barrier to gene flow maintained by selection. They moreover provide genomic insights into the distribution of adaptive potential across O. oratoria' s species range. This work can serve as a case study to characterize adaptive diversity of marine species in the NWP by integrating environmental and genetic data at temporal and spatial scales in a population genomic framework, which would improve management and conservation actions under climate change.

Seascapes Shaped the Local Adaptation and Population Structure of South China Coast Yellowfin Seabream (Acanthopagrus latus)

Understanding the genetic composition and regional adaptation of marine species under environmental heterogeneity and fishing pressure is crucial for responsible management. In order to understand the genetic diversity and adaptability of yellowfin seabream (Acanthopagrus latus) along southern China coast, this study was conducted a seascape genome analysis on yellowfin seabream from the ecologically diverse coast, spanning over 1600 km. A total of 92 yellowfin seabream individuals from 15 sites were performed whole-genome resequencing, and 4,383,564 high-quality single nucleotide polymorphisms (SNPs) were called. By conducting a genotype-environment association analysis, 29,951 adaptive and 4,328,299 neutral SNPs were identified. The yellowfin seabream exhibited two distinct population structures, despite high gene flow between sites. The seascape genome analysis revealed that genetic structure was influenced by a variety of factors including salinity gradients, habitat distance, and ocean currents. The frequency of allelic variation at the candidate loci changed with the salinity gradient. Annotation of these loci revealed that most of the genes are associated with osmoregulation, such as kcnab2a, kcnk5a, and slc47a1. These genes are significantly enriched in pathways associated with ion transport including G protein-coupled receptor activity, transmembrane signaling receptor activity, and transporter activity. Overall, our findings provide insights into how seascape heterogeneity affects adaptive evolution, while providing important information for regional management in yellowfin seabream populations.

Local adaptation with gene flow in a highly dispersive shark

Adaptive divergence in response to environmental clines are expected to be common in species occupying heterogeneous environments. Despite numerous advances in techniques appropriate for non-model species, gene-environment association studies in elasmobranchs are still scarce. The bronze whaler or copper shark (Carcharhinus brachyurus) is a large coastal shark with a wide distribution and one of the most exploited elasmobranchs in southern Africa. Here, we assessed the distribution of neutral and adaptive genomic diversity in C. brachyurus across a highly heterogeneous environment in southern Africa based on genome-wide SNPs obtained through a restriction site-associated DNA method (3RAD). A combination of differentiation-based genome-scan (outflank) and genotype-environment analyses (redundancy analysis, latent factor mixed models) identified a total of 234 differentiation-based outlier and candidate SNPs associated with bioclimatic variables. Analysis of 26,299 putatively neutral SNPs revealed moderate and evenly distributed levels of genomic diversity across sites from the east coast of South Africa to Angola. Multivariate and clustering analyses demonstrated a high degree of gene flow with no significant population structuring among or within ocean basins. In contrast, the putatively adaptive SNPs demonstrated the presence of two clusters and deep divergence between Angola and all other individuals from Namibia and South Africa. These results provide evidence for adaptive divergence in response to a heterogeneous seascape in a large, mobile shark despite high levels of gene flow. These results are expected to inform management strategies and policy at the national and regional level for conservation of C. brachyurus populations.

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.

Whole-genome sequencing reveals fine-scale environment-associated divergence near the range limits of a temperate reef fish

Environmental variation is increasingly recognized as an important driver of diversity in marine species despite the lack of physical barriers to dispersal and the presence of pelagic stages in many taxa. A robust understanding of the genomic and ecological processes involved in structuring populations is lacking for most marine species, often hindering management and conservation action. Cunner (Tautogolabrus adspersus) is a temperate reef fish with both pelagic early life-history stages and strong site-associated homing as adults; the species is also of interest for use as a cleaner fish in salmonid aquaculture in Atlantic Canada. We aimed to characterize genomic and geographic differentiation of cunner in the Northwest Atlantic. To achieve this, a chromosome-level genome assembly for cunner was produced and used to characterize spatial population structure throughout Atlantic Canada using whole-genome sequencing. The genome assembly spanned 0.72 Gbp and 24 chromosomes; whole-genome sequencing of 803 individuals from 20 locations from Newfoundland to New Jersey identified approximately 11 million genetic variants. Principal component analysis revealed four regional Atlantic Canadian groups. Pairwise FST and selection scans revealed signals of differentiation and selection at discrete genomic regions, including adjacent peaks on chromosome 10 across multiple pairwise comparisons (i.e. FST 0.5-0.75). Redundancy analysis suggested association of environmental variables related to benthic temperature and oxygen range with genomic structure. Results suggest regional scale diversity in this temperate reef fish and can directly inform the collection and translocation of cunner for aquaculture applications and the conservation of wild populations throughout the Northwest Atlantic.

Genomics reveals a genetically isolated population of the Pacific white-sided dolphin (Lagenorhynchus obliquidens) distributed in the Sea of Japan

The processes by which animals become genetically isolated in an open environment such as the ocean have not yet been fully elucidated. Morphologically different populations of Pacific white-sided dolphin Lagenorhynchus obliquidens are observed sympatrically in the Sea of Japan. However, genetic studies that exclusively used limited mitochondrial loci or microsatellite DNA have failed to demonstrate the existence of genetically distinct populations. Here, to reveal the population structure, we analysed genome-wide population genetic data using single nucleotide polymorphisms (SNPs) gathered in 2018-2021 from all domestic captive individuals in aquaria, the majority of which originated from the wild, as well as from some stranded individuals, together covering a wide area of coastal water around Japan (n = 123). Multiplexed intersimple sequence repeat genotyping-by-sequencing analysis was performed to obtain the SNP data. Principal coordinate analysis and the clustering method structure indicated that two genetically-distinct populations exist, with little interpopulation gene flow revealed. In addition, the genotypic segregation was reflected in differences in external morphotype. Furthermore, a population demographic analysis based on the whole-genome sequences of an individual from each population indicated that sea-level changes during the Last Glacial Period probably led to allopatric divergence of this species in a limited area of the Sea of Japan, with that group subsequently sharing a distribution area with the other population. These findings yield insights into the formation of genetically isolated sympatric populations in the ocean.

Conservation genomics of an endangered arboreal mammal following the 2019-2020 Australian megafire

The impacts of a changing climate threaten species, populations and ecosystems. Despite these significant and large-scale impacts on threatened species, many remain understudied and have little to no genetic information available. The greater glider, Petauroides volans, is an endangered species highly sensitive to the predicted changes in temperature under a changing climate and was recently severely impacted by a megafire natural disaster (85% estimated population loss). Baseline genetic data is essential for conservation management and for detecting detrimental changes in fire-effected populations. We collected genetic samples within 2 years post the 2019-2020 catastrophic Australian bushfires to examine adaptive potential, baseline genetic diversity and population structure, across their southern range in the state of New South Wales. Population genomic analyses were conducted using 8493 genome-wide SNPs for 86 greater glider individuals across 14 geographic locations. Substantial genetic structure was detected across locations, with low genetic diversity and effective population sizes observed in isolated areas. Additionally, we found signals of putative adaptation in response to temperature in greater gliders using a genotype-environment association analysis. These findings have important implications for the management of greater glider populations by identifying at-risk populations and identifying adaptive potential. We demonstrate the importance of baseline genetic information for endangered species as a practical approach to conservation. This is particularly important given the threat that changes in temperatures and megafire events, as predicted under a changing climate, poses for this species.

Molecular Assessments, Statistical Effectiveness Parameters and Genetic Structure of Captive Populations of Tursiops truncatus Using 15 STRs

Simple Summary

The bottlenose dolphins are one of the most used species in entertainment, assisted therapy, education, and research on welfare. However, their maintenance in captivity requires powerful and sensitive tools for preserving their diversity. The number of genetic markers for this purpose remains controversial, restraining the marine species’ genetic diversity determination. We aimed to select 15 hypervariable molecular markers whose statistical parameters were made in 210 captive dolphins from 18 Mexican centers to support their usefulness. The proposed set of markers allowed us to obtain a genetic fingerprint of each dolphin. Additionally, we identified the structure of the captive population, analyzing the groups according to the capture location. Such characterization is key for maintaining the captive species’ biodiversity rates within conservation and reintroduction programs. However, these 15 genetic markers can also be helpful for small- isolated populations, subspecies and other genera of endangered and vulnerable species.

Abstract

Genetic analysis is a conventional way of identifying and monitoring captive and wildlife species. Knowledge of statistical parameters reinforcing their usefulness and effectiveness as powerful tools for preserving diversity is crucial. Although several studies have reported the diversity of cetaceans such as Tursiops truncatus using microsatellites, its informative degree has been poorly reported. Furthermore, the genetic structure of this cetacean has not been fully studied. In the present study, we selected 15 microsatellites with which 210 dolphins were genetically characterized using capillary electrophoresis. The genetic assertiveness of this set of hypervariable markers identified one individual in the range of 6.927e¹³ to 1.806e¹⁶, demonstrating its substantial capability in kinship relationships. The genetic structure of these 210 dolphins was also determined regarding the putative capture origin; a genetic stratification (k = 2) was found. An additional dolphin group of undetermined origin was also characterized to challenge the proficiency of our chosen markers. The set of markers proposed herein could be a helpful tool to guarantee the maintenance of the genetic diversity rates in conservation programs both in Tursiops truncatus and across other odontocetes, Mysticeti and several genera of endangered and vulnerable species.