Hierarchical population structure and habitat differences in a highly mobile marine species: the Atlantic spotted dolphin

Genetic analysis of harvest samples reveals population structure in a highly mobile generalist carnivore

Delineating wildlife population boundaries is important for effective population monitoring and management. The bobcat (Lynx rufus) is a highly mobile generalist carnivore that is ecologically and economically important. We sampled 1225 bobcats harvested in South Dakota, USA (2014-2019), of which 878 were retained to assess genetic diversity and infer population genetic structure using 17 microsatellite loci. We assigned individuals to genetic clusters (K) using spatial and nonspatial Bayesian clustering algorithms and quantified differentiation (F ST and G ST ″ ) among clusters. We found support for population genetic structure at K = 2 and K = 4, with pairwise F ST and G ST ″ values indicating weak to moderate differentiation, respectively, among clusters. For K = 2, eastern and western clusters aligned closely with historical bobcat management units and were consistent with a longitudinal suture zone for bobcats previously identified in the Great Plains. We did not observe patterns of population genetic structure aligning with major rivers or highways. Genetic divergence observed at K = 4 aligned roughly with ecoregion breaks and may be associated with environmental gradients, but additional sampling with more precise locational data may be necessary to validate these patterns. Our findings reveal that cryptic population structure may occur in highly mobile and broadly distributed generalist carnivores, highlighting the importance of considering population structure when establishing population monitoring programs or harvest regulations. Our study further demonstrates that for elusive furbearers, harvest can provide an efficient, broad-scale sampling approach for genetic population assessments.

Seascape Genetics of the Atlantic Spotted Dolphin (Stenella frontalis) Based on Mitochondrial DNA

The Atlantic spotted dolphin (Stenella frontalis) is endemic to tropical, subtropical, and warm temperate waters of the Atlantic Ocean. Throughout its distribution, both geographic distance and environmental variation may contribute to population structure of the species. In this study, we follow a seascape genetics approach to investigate population differentiation of Atlantic spotted dolphins based on a large worldwide dataset and the relationship with marine environmental variables. The results revealed that the Atlantic spotted dolphin exhibits population genetic structure across its distribution based on mitochondrial DNA control region (mtDNA-CR) data. Analyses based on the contemporary landscape suggested, at both the individual and population level, that the population genetic structure is consistent with the isolation-by-distance model. However, because geography and environmental matrices were correlated, and because in some, but not all analyses, we found a significant effect for the environment, we cannot rule out the addition contribution of environmental factors in structuring genetic variation. Future analyses based on nuclear data are needed to evaluate whether local processes, such as social structure and some level of philopatry within populations, may be contributing to the associations among genetic structure, geographic, and environmental distance.

Sea surface temperature, rather than land mass or geographic distance, may drive genetic differentiation in a species complex of highly dispersive seabirds

Seabirds, particularly Procellariiformes, are highly mobile organisms with a great capacity for long dispersal, though simultaneously showing high philopatry, two conflicting life-history traits that may lead to contrasted patterns of genetic population structure. Landmasses were suggested to explain differentiation patterns observed in seabirds, but philopatry, isolation by distance, segregation between breeding and nonbreeding zones, and oceanographic conditions (sea surface temperatures) may also contribute to differentiation patterns. To our knowledge, no study has simultaneously contrasted the multiple factors contributing to the diversification of seabird species, especially in the gray zone of speciation. We conducted a multilocus phylogeographic study on a widespread seabird species complex, the little shearwater complex, showing highly homogeneous morphology, which led to considerable taxonomic debate. We sequenced three mitochondrial and six nuclear markers on all extant populations from the Atlantic (lherminieri) and Indian Oceans (bailloni), that is, five nominal lineages from 13 populations, along with one population from the eastern Pacific Ocean (representing the dichrous lineage). We found sharp differentiation among populations separated by the African continent with both mitochondrial and nuclear markers, while only mitochondrial markers allowed characterizing the five nominal lineages. No differentiation could be detected within these five lineages, questioning the strong level of philopatry showed by these shearwaters. Finally, we propose that Atlantic populations likely originated from the Indian Ocean. Within the Atlantic, a stepping-stone process accounts for the current distribution. Based on our divergence time estimates, we suggest that the observed pattern of differentiation mostly resulted from historical and current variation in sea surface temperatures.

Fine-scale seascape genomics of an exploited marine species, the common cockle Cerastoderma edule, using a multimodelling approach

Population dynamics of marine species that are sessile as adults are driven by oceanographic dispersal of larvae from spawning to nursery grounds. This is mediated by life-history traits such as the timing and frequency of spawning, larval behaviour and duration, and settlement success. Here, we use 1725 single nucleotide polymorphisms (SNPs) to study the fine-scale spatial genetic structure in the commercially important cockle species Cerastoderma edule and compare it to environmental variables and current-mediated larval dispersal within a modelling framework. Hydrodynamic modelling employing the NEMO Atlantic Margin Model (AMM15) was used to simulate larval transport and estimate connectivity between populations during spawning months (April-September), factoring in larval duration and interannual variability of ocean currents. Results at neutral loci reveal the existence of three separate genetic clusters (mean F ST = 0.021) within a relatively fine spatial scale in the north-west Atlantic. Environmental association analysis indicates that oceanographic currents and geographic proximity explain over 20% of the variance observed at neutral loci, while genetic variance (71%) at outlier loci was explained by sea surface temperature extremes. These results fill an important knowledge gap in the management of a commercially important and overexploited species, bringing us closer to understanding the role of larval dispersal in connecting populations at a fine geographic scale.

Integrating morphological and genetic data at different spatial scales in a cosmopolitan marine turtle species: challenges for management and conservation

Patterns of genetic structure in highly mobile marine vertebrates may be accompanied by phenotypic variation. Most studies in marine turtles focused on population genetic structure have been performed at rookeries. We studied whether genetic and morphological variation of the endangered green turtle (Chelonia mydas) is consistent geographically, focusing on foraging grounds. An association between population genetic structure and body shape variation at broad (inter-lineage) and fine (foraging grounds) scales was predicted and analysed using mitochondrial DNA and geometric morphometrics. Although genetic and phenotypic differentiation patterns were congruent between lineages, no fine-scale association was found, suggesting adaptive divergence. Connectivity among Pacific foraging grounds found here suggests that temperatures of ocean surface currents may influence the genetic structure of C. mydas on a broad scale. Our results suggest that vicariance, dispersal, life-history traits and ecological conditions operating in foraging grounds have shaped the intraspecific morphology and genetic diversity of this species. Considering a range of geographic and temporal scales is useful when management strategies are required for cosmopolitan species. Integrating morphological and genetic tools at different spatial scales, conservation management is proposed based on protection of neutral and adaptive diversity. This approach opens new questions and challenges, especially regarding conservation genetics in cosmopolitan species.

Comparison of persistent organic pollutants (POPs) between small cetaceans in coastal and estuarine waters of the northern Gulf of Mexico

Small cetaceans continue to be exposed to elevated levels of persistent organic pollutants (POPs). The goals of this study were to use data from remote biopsy sampling and photographic-identification to compare POP concentrations between small cetacean stocks in the northern Gulf of Mexico. During 2015-2017, 74 remote biopsies were collected in St. Andrew Bay and adjacent coastal waters from two species: common bottlenose dolphin (Tursiops truncatus) (N = 28, ♀; N = 42, ♂) and Atlantic spotted dolphin (Stenella frontalis) (N = 2, ♀; N = 2, ♂). Common bottlenose dolphin POP concentrations were significantly higher in St. Andrew Bay than coastal waters. Male St. Andrew Bay dolphins had the highest Σ DDT (dichlorodiphenyl-dichloroethane) levels measured in the southeastern U.S. (67 μg/g, 50-89 μg/g; geometric mean and 95% CI) and showed a significant negative relationship between Σ DDT and sighting distance from a St. Andrew Bay point source.

Population genetic structure of Guiana dolphin (Sotalia guianensis) from the southwestern Atlantic coast of Brazil

Sotalia guianensis is a small dolphin that is vulnerable to anthropogenic impacts. Along the Brazilian Atlantic coast, this species is threatened with extinction. A prioritized action plan for conservation strategies relies on increased knowledge of the population. The scarcity of studies about genetic diversity and assessments of population structure for this animal have precluded effective action in the region. Here, we assessed, for the first time, the genetic differentiation at 14 microsatellite loci in 90 S. guianensis specimens stranded on the southeastern Atlantic coast of the State of Espírito Santo, Brazil. We estimated population parameters and structure, measured the significance of global gametic disequilibrium and the intensity of non-random multiallelic interallelic associations and constructed a provisional synteny map using Bos taurus, the closest terrestrial mammal with a reference genome available. All microsatellite loci were polymorphic, with at least three and a maximum of ten alleles each. Allele frequencies ranged from 0.01 to 0.97. Observed heterozygosity ranged from 0.061 to 0.701. The mean inbreeding coefficient was 0.103. Three loci were in Hardy-Weinberg disequilibrium even when missing genotypes were inferred. Although 77 of the 91 possible two-locus associations were in global gametic equilibrium, we unveiled 13 statistically significant, sign-based, non-random multiallelic interallelic associations in 10 two-locus combinations with either coupling (D' values ranging from 0.782 to 0.353) or repulsion (D' values -0.517 to -1.000) forces. Most of the interallelic associations did not involve the major alleles. Thus, for either physically or non-physically linked loci, measuring the intensity of non-random interallelic associations is important for defining the evolutionary forces at equilibrium. We uncovered a small degree of genetic differentiation (F_ST = 0.010; P-value = 0.463) with a hierarchical clustering into one segment containing members from the southern and northern coastal regions. The data thus support the scenario of little genetic structure in the population of S. guianensis in this geographic area.

Spatially Explicit Analysis of Genome-Wide SNPs Detects Subtle Population Structure in a Mobile Marine Mammal, the Harbor Porpoise

The population structure of the highly mobile marine mammal, the harbor porpoise (Phocoena phocoena), in the Atlantic shelf waters follows a pattern of significant isolation-by-distance. The population structure of harbor porpoises from the Baltic Sea, which is connected with the North Sea through a series of basins separated by shallow underwater ridges, however, is more complex. Here, we investigated the population differentiation of harbor porpoises in European Seas with a special focus on the Baltic Sea and adjacent waters, using a population genomics approach. We used 2872 single nucleotide polymorphisms (SNPs), derived from double digest restriction-site associated DNA sequencing (ddRAD-seq), as well as 13 microsatellite loci and mitochondrial haplotypes for the same set of individuals. Spatial principal components analysis (sPCA), and Bayesian clustering on a subset of SNPs suggest three main groupings at the level of all studied regions: the Black Sea, the North Atlantic, and the Baltic Sea. Furthermore, we observed a distinct separation of the North Sea harbor porpoises from the Baltic Sea populations, and identified splits between porpoise populations within the Baltic Sea. We observed a notable distinction between the Belt Sea and the Inner Baltic Sea sub-regions. Improved delineation of harbor porpoise population assignments for the Baltic based on genomic evidence is important for conservation management of this endangered cetacean in threatened habitats, particularly in the Baltic Sea proper. In addition, we show that SNPs outperform microsatellite markers and demonstrate the utility of RAD-tags from a relatively small, opportunistically sampled cetacean sample set for population diversity and divergence analysis.

Genomic Methods Take the Plunge: Recent Advances in High-Throughput Sequencing of Marine Mammals

The dramatic increase in the application of genomic techniques to non-model organisms (NMOs) over the past decade has yielded numerous valuable contributions to evolutionary biology and ecology, many of which would not have been possible with traditional genetic markers. We review this recent progression with a particular focus on genomic studies of marine mammals, a group of taxa that represent key macroevolutionary transitions from terrestrial to marine environments and for which available genomic resources have recently undergone notable rapid growth. Genomic studies of NMOs utilize an expanding range of approaches, including whole genome sequencing, restriction site-associated DNA sequencing, array-based sequencing of single nucleotide polymorphisms and target sequence probes (e.g., exomes), and transcriptome sequencing. These approaches generate different types and quantities of data, and many can be applied with limited or no prior genomic resources, thus overcoming one traditional limitation of research on NMOs. Within marine mammals, such studies have thus far yielded significant contributions to the fields of phylogenomics and comparative genomics, as well as enabled investigations of fitness, demography, and population structure. Here we review the primary options for generating genomic data, introduce several emerging techniques, and discuss the suitability of each approach for different applications in the study of NMOs.

Contrasting patterns of population connectivity between regions in a commercially important mollusc Haliotis rubra: integrating population genetics, genomics and marine LiDAR data

Estimating contemporary genetic structure and population connectivity in marine species is challenging, often compromised by genetic markers that lack adequate sensitivity, and unstructured sampling regimes. We show how these limitations can be overcome via the integration of modern genotyping methods and sampling designs guided by LiDAR and SONAR data sets. Here we explore patterns of gene flow and local genetic structure in a commercially harvested abalone species (Haliotis rubra) from southeastern Australia, where the viability of fishing stocks is believed to be dictated by recruitment from local sources. Using a panel of microsatellite and genomewide SNP markers, we compare allele frequencies across a replicated hierarchical sampling area guided by bathymetric LiDAR imagery. Results indicate high levels of gene flow and no significant genetic structure within or between benthic reef habitats across 1400 km of coastline. These findings differ to those reported for other regions of the fishery indicating that larval supply is likely to be spatially variable, with implications for management and long-term recovery from stock depletion. The study highlights the utility of suitably designed genetic markers and spatially informed sampling strategies for gaining insights into recruitment patterns in benthic marine species, assisting in conservation planning and sustainable management of fisheries.