Genetic structure of harbour porpoise Phocoena phocoena populations in the northwest Atlantic based on mitochondrial and nuclear markers

Using eDNA to Supplement Population Genetic Analyses for Cryptic Marine Species: Identifying Population Boundaries for Alaska Harbour Porpoises

Isolation by distance and biogeographical boundaries define patterns of population genetic structure for harbour porpoise along the Pacific coast from California to British Columbia. Until recently, inadequate sample sizes in many regions constrained efforts to characterise population genetic structure throughout the coastal waters of Alaska. Here, tissue samples from beachcast strandings and fisheries bycatch were supplemented with targeted environmental DNA (eDNA) samples in key regions of Alaska coastal and inland waters. Using a geographically explicit, hierarchical approach, we examined the genetic structure of Alaska harbour porpoises, using both mitochondrial DNA (mtDNA) sequence data and multilocus SNP genotypes. Despite a lack of evidence of genetic differentiation from nuclear SNP loci, patterns of relatedness and genetic differentiation from mtDNA suggest natal philopatry at multiple geographic scales, with limited gene flow among sites possibly mediated by male dispersal. A priori clustering of sampled areas at an intermediate scale (eastern and western Bering Sea, Gulf of Alaska and Southeast Alaska) best explained the genetic variance (12.37%) among regions. In addition, mtDNA differentiation between the Gulf of Alaska and eastern Bering Sea, and among regions within the Gulf of Alaska, indicated significant genetic structuring of harbour porpoise populations in Southeast Alaska. The targeted collection of eDNA samples from strata within Southeast Alaska was key for elevating the statistical power of our mtDNA dataset, and findings indicate limited dispersal between neighbouring strata within coastal and inland waters. These results provide evidence supporting a population boundary within the currently recognised Southeast Alaska Stock. Together, these findings will prove useful for ongoing management efforts to reduce fisheries conflict and conserve genetic diversity in this iconic coastal species.

Are There Barriers Separating the Pink River Dolphin Populations (Inia boliviensis, Iniidae, Cetacea) within the Mamoré-Iténez River Basins (Bolivia)? An Analysis of Its Genetic Structure by Means of Mitochondrial and Nuclear DNA Markers

The pink river dolphin, or bufeo, is one of the dolphins which lives in the rivers of the Orinoco and Amazon basins in South America. The Bolivian bufeo population is considered a differentiated species (Inia boliviensis) from the Amazon and Orinoco species (Inia geoffrensis). Until now, no study has completed an extensive population genetics analysis of the bufeo in Bolivian rivers. We analyzed 82 bufeos from different rivers from the Mamoré and Iténez (Guaporé) river basins for the mt control region (CR), nuclear microsatellites, and DQB-1 gene sequences to determine if the inner rapids of these Bolivian river basins have some influence on the genetic structure of this species. The first relevant result was that the genetic diversity for CR, and the microsatellites were substantially lower in the Bolivian bufeos than in the dolphins studied in other areas of the Amazon and Orinoco. However, the DQB-1 gene sequences yielded similar genetic diversity to those found in other areas. The second relevant result is the existence of some significant genetic heterogeneity among the bufeo populations within Bolivia, although in a small degree, but this differentiation is independent of the inner rapids of the Bolivian rivers we sampled. The third relevant result was the existence of significant isolation by distance for the CR, but not for microsatellites and DQB-1 gene sequences. This was related to differential gene flow capacity of females (philopatric) and males (less philopatric and more migrants) and, possibly, to different selective patterns affecting the molecular markers studied. The fourth relevant result was related to diverse demographic changes of these bufeos. At least two or three bottleneck events and one or two population expansions have occurred in the Bolivian bufeo population. The major part of these events occurred during the Pleistocene.

No leading-edge effect in North Atlantic harbor porpoises: Evolutionary and conservation implications

Understanding species responses to past environmental changes can help forecast how they will cope with ongoing climate changes. Harbor porpoises are widely distributed in the North Atlantic and were deeply impacted by the Pleistocene changes with the split of three subspecies. Despite major impacts of fisheries on natural populations, little is known about population connectivity and dispersal, how they reacted to the Pleistocene changes, and how they will evolve in the future. Here, we used phylogenetics, population genetics, and predictive habitat modeling to investigate population structure and phylogeographic history of the North Atlantic porpoises. A total of 925 porpoises were characterized at 10 microsatellite loci and one quarter of the mitogenome (mtDNA). A highly divergent mtDNA lineage was uncovered in one porpoise off Western Greenland, suggesting that a cryptic group may occur and could belong to a recently discovered mesopelagic ecotype off Greenland. Aside from it and the southern subspecies, spatial genetic variation showed that porpoises from both sides of the North Atlantic form a continuous system belonging to the same subspecies (Phocoena phocoena phocoena). Yet, we identified important departures from random mating and restricted dispersal forming a highly significant isolation by distance (IBD) at both mtDNA and nuclear markers. A ten times stronger IBD at mtDNA compared with nuclear loci supported previous evidence of female philopatry. Together with the lack of spatial trends in genetic diversity, this IBD suggests that migration-drift equilibrium has been reached, erasing any genetic signal of a leading-edge effect that accompanied the predicted recolonization of the northern habitats freed from Pleistocene ice. These results illuminate the processes shaping porpoise population structure and provide a framework for designing conservation strategies and forecasting future population evolution.

Population structure in a continuously distributed coastal marine species, the harbor porpoise, based on microhaplotypes derived from poor-quality samples

Harbor porpoise in the North Pacific are found in coastal waters from southern California to Japan, but population structure is poorly known outside of a few local areas. We used multiplexed amplicon sequencing of 292 loci and genotyped clusters of single nucleotide polymoirphisms as microhaplotypes (N = 271 samples) in addition to mitochondrial (mtDNA) sequence data (N = 413 samples) to examine the genetic structure from samples collected along the Pacific coast and inland waterways from California to southern British Columbia. We confirmed an overall pattern of strong isolation-by-distance, suggesting that individual dispersal is restricted. We also found evidence of regions where genetic differences are larger than expected based on geographical distance alone, implying current or historical barriers to gene flow. In particular, the southernmost population in California is genetically distinct (F_ST  = 0.02 [microhaplotypes]; 0.31 [mtDNA]), with both reduced genetic variability and high frequency of an otherwise rare mtDNA haplotype. At the northern end of our study range, we found significant genetic differentiation of samples from the Strait of Georgia, previously identified as a potential biogeographical boundary or secondary contact zone between harbor porpoise populations. Association of microhaplotypes with remotely sensed environmental variables indicated potential local adaptation, especially at the southern end of the species' range. These results inform conservation and management for this nearshore species, illustrate the value of genomic methods for detecting patterns of genetic structure within a continuously distributed marine species, and highlight the power of microhaplotype genotyping for detecting genetic structure in harbor porpoises despite reliance on poor-quality samples.

Kin relationships in cultural species of the marine realm: case study of a matrilineal social group of sperm whales off Mauritius island, Indian Ocean

Understanding the organization and dynamics of social groups of marine mammals through the study of kin relationships is particularly challenging. Here, we studied a stable social group of sperm whales off Mauritius, using underwater observations, individual-specific identification, non-invasive sampling and genetic analyses based on mitochondrial sequencing and microsatellite profiling. Twenty-four sperm whales were sampled between 2017 and 2019. All individuals except one adult female shared the same mitochondrial DNA (mtDNA) haplotype-one that is rare in the western Indian Ocean-thus confirming with near certainty the matrilineality of the group. All probable first- and second-degree kin relationships were depicted in the sperm whale social group: 13 first-degree and 27 second-degree relationships were identified. Notably, we highlight the likely case of an unrelated female having been integrated into a social unit, in that she presented a distinct mtDNA haplotype and no close relationships with any members of the group. Investigating the possible matrilineality of sperm whale cultural units (i.e. vocal clans) is the next step in our research programme to elucidate and better apprehend the complex organization of sperm whale social groups.

Identifying mechanisms of genetic differentiation among populations in vagile species: historical factors dominate genetic differentiation in seabirds

Elucidating the factors underlying the origin and maintenance of genetic variation among populations is crucial for our understanding of their ecology and evolution, and also to help identify conservation priorities. While intrinsic movement has been hypothesized as the major determinant of population genetic structuring in abundant vagile species, growing evidence indicates that vagility does not always predict genetic differentiation. However, identifying the determinants of genetic structuring can be challenging, and these are largely unknown for most vagile species. Although, in principle, levels of gene flow can be inferred from neutral allele frequency divergence among populations, underlying assumptions may be unrealistic. Moreover, molecular studies have suggested that contemporary gene flow has often not overridden historical influences on population genetic structure, which indicates potential inadequacies of any interpretations that fail to consider the influence of history in shaping that structure. This exhaustive review of the theoretical and empirical literature investigates the determinants of population genetic differentiation using seabirds as a model system for vagile taxa. Seabirds provide a tractable group within which to identify the determinants of genetic differentiation, given their widespread distribution in marine habitats and an abundance of ecological and genetic studies conducted on this group. Herein we evaluate mitochondrial DNA (mtDNA) variation in 73 seabird species. Lack of mutation-drift equilibrium observed in 19% of species coincided with lower estimates of genetic differentiation, suggesting that dynamic demographic histories can often lead to erroneous interpretations of contemporary gene flow, even in vagile species. Presence of land across the species sampling range, or sampling of breeding colonies representing ice-free Pleistocene refuge zones, appear to be associated with genetic differentiation in Tropical and Southern Temperate species, respectively, indicating that long-term barriers and persistence of populations are important for their genetic structuring. Conversely, biotic factors commonly considered to influence population genetic structure, such as spatial segregation during foraging, were inconsistently associated with population genetic differentiation. In light of these results, we recommend that genetic studies should consider potential historical events when identifying determinants of genetic differentiation among populations to avoid overestimating the role of contemporary factors, even for highly vagile taxa.

Global phylogeography and genetic diversity of the long-finned pilot whale Globicephala melas, with new data from the southeastern Pacific

The matrilineal long-finned pilot whale presents an antitropical distribution and is divided into two subspecies, one in the temperate seas of the Southern Hemisphere and the other restricted to the North Atlantic and Mediterranean. Until now, population genetic and phylogeographic studies have included localities of most of its Northern Hemisphere distribution, while only the southwestern Pacific has been sampled in the Southern Hemisphere. We add new genetic data from the southeastern Pacific to the published sequences. Low mitochondrial and nuclear diversity was encountered in this new area, as previously reported for other localities. Four haplotypes were found with only one new for the species. Fifteen haplotypes were detected in the global dataset, underlining the species’ low diversity. As previously reported, the subspecies shared two haplotypes and presented a strong phylogeographic structure. The extant distribution of this species has been related to dispersal events during the Last Glacial Maximum. Using the genetic data and Approximate Bayesian Calculations, this study supports this historical biogeographic scenario. From a taxonomic perspective, even if genetic analyses do not support the subspecies category, this study endorses the incipient divergence process between hemispheres, thus maintaining their status and addressing them as Demographically Independent Populations is recommended.

Population Structure of Riverine and Coastal Dolphins Sotalia fluviatilis and Sotalia guianensis: Patterns of Nuclear and Mitochondrial Diversity and Implications for Conservation

Coastal and freshwater cetaceans are particularly vulnerable due to their proximity to human activity, localized distributions, and small home ranges. These species include Sotalia guianensis, found in the Atlantic and Caribbean coastal areas of central and South America, and Sotalia fluviatilis, distributed in the Amazon River and tributaries. We investigated the population structure and genetic diversity of these 2 species by analyses of mtDNA control region and 8-10 microsatellite loci. MtDNA analyses revealed strong regional structuring for S. guianensis (i.e., Colombian Caribbean vs. Brazilian Coast, FST = 0.807, ΦST = 0.878, P < 0.001) especially north and south of the Amazon River mouth. For S. fluviatilis, population structuring was detected between the western and eastern Amazon (i.e., Colombian Amazon vs. Brazilian Amazon, FST = 0.085, ΦST = 0.277, P < 0.001). Haplotype and nucleotide diversity were higher for S. fluviatilis. Population differentiation was supported by analysis of the microsatellite loci (S. guianensis, northern South America vs. southern South America FST = 0.275, Jost's D = 0.476, P < 0.001; S. fluviatilis, western and eastern Amazon FST = 0.197, Jost's D = 0.364, P < 0.001). Most estimated migration rates in both species overlapped with zero, suggesting no measurable migration between most of the sampling locations. However, for S. guianensis, there was measurable migration in neighboring sampling locations. These results indicate that the small home ranges of these species may act to restrict gene flow between populations separated by relatively short distances, increasing the risk of extirpation of some localized populations in the future if existing threats are not minimized.

Water, water everywhere: environmental DNA can unlock population structure in elusive marine species

Determining management units for natural populations is critical for effective conservation and management. However, collecting the requisite tissue samples for population genetic analyses remains the primary limiting factor for a number of marine species. The harbour porpoise (Phocoena phocoena), one of the smallest cetaceans in the Northern Hemisphere, is a primary example. These elusive, highly mobile small animals confound traditional approaches of collecting tissue samples for genetic analyses, yet their nearshore habitat makes them highly vulnerable to fisheries by-catch and the effects of habitat degradation. By exploiting the naturally shed cellular material in seawater and the power of next-generation sequencing, we develop a novel approach for generating population-specific mitochondrial sequence data from environmental DNA (eDNA) using surface seawater samples. Indications of significant genetic differentiation within a currently recognized management stock highlights the need for dedicated eDNA sampling throughout the population's range in southeast Alaska. This indirect sampling tactic for characterizing stock structure of small and endangered marine mammals has the potential to revolutionize population assessment for otherwise inaccessible marine taxa.

Genetic footprint of population fragmentation and contemporary collapse in a freshwater cetacean

Understanding demographic trends and patterns of gene flow in an endangered species is crucial for devising conservation strategies. Here, we examined the extent of population structure and recent evolution of the critically endangered Yangtze finless porpoise (Neophocaena asiaeorientalis asiaeorientalis). By analysing genetic variation at the mitochondrial and nuclear microsatellite loci for 148 individuals, we identified three populations along the Yangtze River, each one connected to a group of admixed ancestry. Each population displayed extremely low genetic diversity, consistent with extremely small effective size (≤92 individuals). Habitat degradation and distribution gaps correlated with highly asymmetric gene-flow that was inefficient in maintaining connectivity between populations. Genetic inferences of historical demography revealed that the populations in the Yangtze descended from a small number of founders colonizing the river from the sea during the last Ice Age. The colonization was followed by a rapid population split during the last millennium predating the Chinese Modern Economy Development. However, genetic diversity showed a clear footprint of population contraction over the last 50 years leaving only ~2% of the pre-collapsed size, consistent with the population collapses reported from field studies. This genetic perspective provides background information for devising mitigation strategies to prevent this species from extinction.

Does silvoagropecuary landscape fragmentation affect the genetic diversity of the sigmodontine rodent Oligoryzomys longicaudatus?

BACKGROUND

Fragmentation of native forests is a highly visible result of human land-use throughout the world. In this study, we evaluated the effects of landscape fragmentation and matrix features on the genetic diversity and structure of Oligoryzomys longicaudatus, the natural reservoir of Hantavirus in southern South America. We focused our work in the Valdivian Rainforest where human activities have produced strong change of natural habitats, with an important number of human cases of Hantavirus.

METHODS

We sampled specimens of O. longicaudatus from five native forest patches surrounded by silvoagropecuary matrix from Panguipulli, Los Rios Region, Chile. Using the hypervariable domain I (mtDNA), we characterized the genetic diversity and evaluated the effect of fragmentation and landscape matrix on the genetic structure of O. longicaudatus. For the latter, we used three approaches: (i) Isolation by Distance (IBD) as null model, (ii) Least-cost Path (LCP) where genetic distances between patch pairs increase with cost-weighted distances, and (iii) Isolation by Resistance (IBR) where the resistance distance is the average number of steps that is needed to commute between the patches during a random walk.

RESULTS

We found low values of nucleotide diversity (π) for the five patches surveyed, ranging from 0.012 to 0.015, revealing that the 73 sampled specimens of this study belong to two populations but with low values of genetic distance (γST ) ranging from 0.022 to 0.099. Likewise, we found that there are no significant associations between genetic distance and geographic distance for IBD and IBR. However, we found for the LCP approach, a significant positive relationship (r = 0.737, p = 0.05), with shortest least-cost paths traced through native forest and arborescent shrublands.

DISCUSSION

In this work we found that, at this reduced geographical scale, Oligoryzomys longicaudatus shows genetic signs of fragmentation. In addition, we found that connectivity between full growth native forest remnants is mediated by the presence of dense shrublands and native forest corridors. In this sense, our results are important because they show how native forest patches and associated routes act as source of vector species in silvoagropecuary landscape, increasing the infection risk on human population. This study is the first approach to understand the epidemiological spatial context of silvoagropecuary risk of Hantavirus emergence. Further studies are needed to elucidate the effects of landscape fragmentation in order to generate new predictive models based on vector intrinsic attributes and landscape features.

Mixing of porpoise ecotypes in southwestern UK waters revealed by genetic profiling

Contact zones between ecotypes are windows for understanding how species may react to climate changes. Here, we analysed the fine-scale genetic and morphological variation in harbour porpoises (Phocoena phocoena) around the UK by genotyping 591 stranded animals at nine microsatellite loci. The data were integrated with a prior study to map at high resolution the contact zone between two previously identified ecotypes meeting in the northern Bay of Biscay. Clustering and spatial analyses revealed that UK porpoises are derived from two genetic pools with porpoises from the southwestern UK being genetically differentiated, and having larger body sizes compared to those of other UK areas. Southwestern UK porpoises showed admixed ancestry between southern and northern ecotypes with a contact zone extending from the northern Bay of Biscay to the Celtic Sea and Channel. Around the UK, ancestry blends from one genetic group to the other along a southwest--northeast axis, correlating with body size variation, consistent with previously reported morphological differences between the two ecotypes. We also detected isolation by distance among juveniles but not in adults, suggesting that stranded juveniles display reduced intergenerational dispersal. The fine-scale structure of this admixture zone raises the question of how it will respond to future climate change and provides a reference point for further study.

Population structure and cryptic speciation in bonnethead sharks Sphyrna tiburo in the south-eastern U.S.A. and Caribbean

Population structure and lineage diversification within a small, non-dispersive hammerhead shark species, the bonnethead shark Sphyrna tiburo, was assessed. Sphyrna tiburo is currently described as one continuously distributed species along the Atlantic continental margins of North, Central and South America, but recent genetic analysis of an insular population (Trinidad) suggests the possibility of cryptic speciation. To address this issue S. tiburo were sampled at six sites along c. 6200 km of continuous, continental coastline and from one island location (Grand Bahama) across a discontinuity in their distribution (the Straits of Florida), in order to test if they constitute a single lineage over this distribution. A total of 1030 bp of the mitochondrial control region (CR) was obtained for 239 S. tiburo, revealing 73 distinct haplotypes, high nucleotide diversity (0·01089) and a pair of highly divergent lineages estimated to have separated 3·61-5·62 million years ago. Mitochondrial cytochrome oxidase I and nuclear internal transcribed spacer loci show the same pattern. Divergence is similar within S. tiburo to that observed between established elasmobranch sister species, providing further evidence of cryptic speciation. A global AMOVA based on CR confirms that genetic diversity is primarily partitioned among populations (Φ_ST = 0·828, P < 0·001) because the divergent lineages are almost perfectly segregated between Belize and North America-The Bahamas. An AMOVA consisting solely of the North American and Bahamian samples is also significantly different from zero (Φ_ST = 0·088, P < 0·001) and pairwise F_ST is significantly different between all sites. These findings suggest that S. tiburo comprises a species complex and supports previous research indicating fine population structure, which has implications for fisheries management and biodiversity conservation.

Historical dimensions of population structure in a continuously distributed marine species: The case of the endemic Chilean dolphin

The complementarity of historical and contemporary processes contributes to understanding the genetic structure of continuously distributed marine species with high dispersal capabilities. Cephalorhynchus eutropia, has a continuous coastal distribution with strong genetic differentiation identified by nuclear DNA markers. We explored the historical dimension of this genetic differentiation between northern and southern populations to evaluate phylogeographic structure. Additionally, we conducted mtDNA and microsatellite analyses to detect past and recent demographic changes. The southern population was characterized by lower genetic diversity with a signal of population expansion, likely associated with ice retreat and habitat extension after the Last Glacial Maximum (LGM). In contrast, structure within the northern population was more consistent with stable historical population size. Approximate Bayesian Computation analyses suggested that during the LGM, C. eutropia persisted in the northern area; while the south was colonized by dispersal ~11,000 years ago followed by population expansion. This study shows that Chilean dolphin population structure is consistent with predictions from the Expansion-Contraction biogeographic model, with a poleward post-glacial shift revealed in current genetic structure. The results also confirm the validity of the population units previously identified, demonstrating their historical origin and highlighting the utility of integrating genetic markers with different temporal scale resolutions.

Parentage-Based Group Composition and Dispersal Pattern Studies of the Yangtze Finless Porpoise Population in Poyang Lake

Social behaviors are poorly known for the critically endangered Yangtze finless porpoise (YFP, Neophocaena asiaeorientalis asiaeorientalis). Here, group composition and dispersal patterns of the YFP population living in the Poyang Lake were studied by parentage-based pedigree analyses using 21 microsatellite loci and a 597 bp segment of the mitochondrial DNA control region. In this study, 21 potential mother-offspring pairs and six potential father-offspring pairs (including two potential parents-offspring pairs) were determined, among which 12 natural mother-offspring groups and a maternal group of three generations were found. No genetically-determined fathers were found associated with their offspring. This study also found that maternally related porpoises at the reproductive state tend to group together. This suggest maternal relationship and reproductive state may be factors for grouping in the YFP population. In natural mother-offspring groups, male offspring were all younger than two years old, which suggest male offspring may leave their mothers at approximately two years of age, or at least they were not in tight association with their mothers as they may have been under two years old. However, female offspring can stay longer with their mothers and can reproduce in the natal group.

Is there a distinct harbor porpoise subpopulation in the Marmara Sea?

Genetic population structure of geographically isolated endangered Black Sea harbor porpoise (Phocoena phocoena relicta) is little known in Turkish waters, especially in the Turkish Straits System (TSS- Marmara Sea, Bosphorus and Dardanelles), which connects the Black Sea and the Aegean Sea. Mitochondrial DNA sequences of 70 new individuals sampled in the Turkish Black Sea, TSS and Aegean Sea, revealed five new haplotypes from the Black Sea. The findings support the idea that harbor porpoises from the Black Sea dispersed into the Aegean through the TSS. Considering signatures of population expansion, all subpopulations showed a signature of population expansion. The network data and the Фst calculations indicated that the Marmara Sea subpopulation was significantly differentiated from all of the other subpopulations, and supports the notion of its isolated. The finding of a potential management unit (MU) within an already heavily impacted subpopulation as a whole suggests that the individuals of P. p. relicta inhabiting the Marmara Sea require a very rigorous conservation strategy to ensure the survival of this subpopulation, represented by its unique haplotype.

Population Genetics of Franciscana Dolphins (Pontoporia blainvillei): Introducing a New Population from the Southern Edge of Their Distribution

Due to anthropogenic factors, the franciscana dolphin, Pontoporia blainvillei, is the most threatened small cetacean on the Atlantic coast of South America. Four Franciscana Management Areas have been proposed: Espiritu Santo to Rio de Janeiro (FMA I), São Paulo to Santa Catarina (FMA II), Rio Grande do Sul to Uruguay (FMA III), and Argentina (FMA IV). Further genetic studies distinguished additional populations within these FMAs. We analyzed the population structure, phylogeography, and demographic history in the southernmost portion of the species range. From the analysis of mitochondrial DNA control region sequences, 5 novel haplotypes were found, totalizing 60 haplotypes for the entire distribution range. The haplotype network did not show an apparent phylogeographical signal for the southern FMAs. Two populations were identified: Monte Hermoso (MH) and Necochea (NC)+Claromecó (CL)+Río Negro (RN). The low levels of genetic variability, the relative constant size over time, and the low levels of gene flow may indicate that MH has been colonized by a few maternal lineages and became isolated from geographically close populations. The apparent increase in NC+CL+RN size would be consistent with the higher genetic variability found, since genetic diversity is generally higher in older and expanding populations. Additionally, RN may have experienced a recent split from CL and NC; current high levels of gene flow may be occurring between the latter ones. FMA IV would comprise four franciscana dolphin populations: Samborombón West+Samborombón South, Cabo San Antonio+Buenos Aires East, NC+CL+Buenos Aires Southwest+RN and MH. Results achieved in this study need to be taken into account in order to ensure the long-term survival of the species.

Microsatellite markers reveal strong genetic structure in the endemic Chilean dolphin

Understanding genetic differentiation and speciation processes in marine species with high dispersal capabilities is challenging. The Chilean dolphin, Cephalorhynchus eutropia, is the only endemic cetacean of Chile and is found in two different coastal habitats: a northern habitat with exposed coastlines, bays and estuaries from Valparaíso (33°02′S) to Chiloé (42°00′S), and a southern habitat with highly fragmented inshore coastline, channels and fjords between Chiloé and Navarino Island (55°14′S). With the aim of evaluating the potential existence of conservation units for this species, we analyzed the genetic diversity and population structure of the Chilean dolphin along its entire range. We genotyped 21 dinucleotide microsatellites for 53 skin samples collected between 1998 and 2012 (swab: n = 8, biopsy: n = 38, entanglement n = 7). Bayesian clustering and spatial model analyses identified two genetically distinct populations corresponding to the northern and southern habitats. Genetic diversity levels were similar in the two populations (He: 0.42 v/s 0.45 for southern and northern populations, respectively), while effective size population was higher in the southern area (Ne: 101 v/s 39). Genetic differentiation between these two populations was high and significant (F_ST = 0.15 and R_ST = 0.19), indicating little or no current gene flow. Because of the absence of evident geographical barriers between the northern and southern populations, we propose that genetic differentiation may reflect ecological adaptation to the different habitat conditions and resource uses. Therefore, the two genetic populations of this endemic and Near Threatened species should be considered as different conservation units with independent management strategies.

Multiple populations of pantropical spotted dolphins in Hawaiian waters

Understanding gene flow and dispersal patterns is important for predicting effects of natural events and anthropogenic activities on animal populations. In Hawaii, most species of odontocetes are managed as single populations. Recent exceptions include false killer whales, spinner dolphins, and common bottlenose dolphins, for which studies have shown fidelity to individual islands or groups of islands. Our study focused on pantropical spotted dolphins. We analyzed mitochondrial control region and 11 microsatellite loci from 101 individuals from 4 areas: Hawaii, Maui/Lanai, Oahu, and Kauai/Niihau. We examined F ST, F' ST, R ST, Jost's D, and ΦST and used TESS to estimate number of populations and assignment probabilities. Our results support genetic differentiation among Hawaii, Maui/Lanai, and Oahu and suggest that pantropical spotted dolphins near Kauai/Niihau are likely transient and in low numbers. Between island regions, F ST for microsatellites ranged from 0.016 to 0.045 and for mtDNA, from 0.011 to 0.282. F ' ST, ranged from 0.098 to 0.262 for microsatellites and 0.019 to 0.415 for mtDNA. R ST and ΦST showed similar results to F ST for microsatellites and mtDNA respectively, and Jost's D fell between F ST and F ' ST. TESS supported 3 populations, and greatest mean assignment probability by island region ranged from 0.50 to 0.72. The private alleles method indicated migration rates among regions from 1.49 to 3.45, and effective population size of the island of Hawaii was estimated to be 220. There was no strong evidence to support sex-biased dispersal or group fidelity. Considering this study in the larger context of other odontocete population studies and studies of connectivity, we suggest genetic differentiation may be mediated by behavior adapted to differing habitat types and niches.

Genetic diversity and population structure of the critically endangered Yangtze finless porpoise (Neophocaena asiaeorientalis asiaeorientalis) as revealed by mitochondrial and microsatellite DNA

Ecological surveys have indicated that the population of the critically endangered Yangtze finless porpoise (YFP, Neophocaena asiaeorientalis asiaeorientalis) is becoming increasingly small and fragmented, and will be at high risk of extinction in the near future. Genetic conservation of this population will be an important component of the long-term conservation effort. We used a 597 base pair mitochondrial DNA (mtDNA) control region and 11 microsatellite loci to analyze the genetic diversity and population structure of the YFP. The analysis of both mtDNA and microsatellite loci suggested that the genetic diversity of the YFP will possibly decrease in the future if the population keeps declining at a rapid rate, even though these two types of markers revealed different levels of genetic diversity. In addition, mtDNA revealed strong genetic differentiation between one local population, Xingchang–Shishou (XCSS), and the other five downstream local populations; furthermore, microsatellite DNA unveiled fine but significant genetic differentiation between three of the local populations (not only XCSS but also Poyang Lake (PY) and Tongling (TL)) and the other local populations. With an increasing number of distribution gaps appearing in the Yangtze main steam, the genetic differentiation of local populations will likely intensify in the future. The YFP is becoming a genetically fragmented population. Therefore, we recommend attention should be paid to the genetic conservation of the YFP.

Patterns of population structure for inshore bottlenose dolphins along the eastern United States

Globally distributed, the bottlenose dolphin (Tursiops truncatus) is found in a range of offshore and coastal habitats. Using 15 microsatellite loci and mtDNA control region sequences, we investigated patterns of genetic differentiation among putative populations along the eastern US shoreline (the Indian River Lagoon, Florida, and Charleston Harbor, South Carolina) (microsatellite analyses: n = 125, mtDNA analyses: n = 132). We further utilized the mtDNA to compare these populations with those from the Northwest Atlantic, Gulf of Mexico, and Caribbean. Results showed strong differentiation among inshore, alongshore, and offshore habitats (ФST = 0.744). In addition, Bayesian clustering analyses revealed the presence of 2 genetic clusters (populations) within the 250 km Indian River Lagoon. Habitat heterogeneity is likely an important force diversifying bottlenose dolphin populations through its influence on social behavior and foraging strategy. We propose that the spatial pattern of genetic variation within the lagoon reflects both its steep longitudinal transition of climate and also its historical discontinuity and recent connection as part of Intracoastal Waterway development. These findings have important management implications as they emphasize the role of habitat and the consequence of its modification in shaping bottlenose dolphin population structure and highlight the possibility of multiple management units existing in discrete inshore habitats along the entire eastern US shoreline.

Low diversity in the mitogenome of sperm whales revealed by next-generation sequencing

Large population sizes and global distributions generally associate with high mitochondrial DNA control region (CR) diversity. The sperm whale (Physeter macrocephalus) is an exception, showing low CR diversity relative to other cetaceans; however, diversity levels throughout the remainder of the sperm whale mitogenome are unknown. We sequenced 20 mitogenomes from 17 sperm whales representative of worldwide diversity using Next Generation Sequencing (NGS) technologies (Illumina GAIIx, Roche 454 GS Junior). Resequencing of three individuals with both NGS platforms and partial Sanger sequencing showed low discrepancy rates (454-Illumina: 0.0071%; Sanger-Illumina: 0.0034%; and Sanger-454: 0.0023%) confirming suitability of both NGS platforms for investigating low mitogenomic diversity. Using the 17 sperm whale mitogenomes in a phylogenetic reconstruction with 41 other species, including 11 new dolphin mitogenomes, we tested two hypotheses for the low CR diversity. First, the hypothesis that CR-specific constraints have reduced diversity solely in the CR was rejected as diversity was low throughout the mitogenome, not just in the CR (overall diversity π = 0.096%; protein-coding 3rd codon = 0.22%; CR = 0.35%), and CR phylogenetic signal was congruent with protein-coding regions. Second, the hypothesis that slow substitution rates reduced diversity throughout the sperm whale mitogenome was rejected as sperm whales had significantly higher rates of CR evolution and no evidence of slow coding region evolution relative to other cetaceans. The estimated time to most recent common ancestor for sperm whale mitogenomes was 72,800 to 137,400 years ago (95% highest probability density interval), consistent with previous hypotheses of a bottleneck or selective sweep as likely causes of low mitogenome diversity.

A European melting pot of harbour porpoise in the French Atlantic coasts inferred from mitochondrial and nuclear data

Field surveys have reported a global shift in harbour porpoise distribution in European waters during the last 15 years, including a return to the Atlantic coasts of France. In this study, we analyzed genetic polymorphisms at a fragment of the mitochondrial control region (mtDNA CR) and 7 nuclear microsatellite loci, for 52 animals stranded and by-caught between 2000 and 2010 along the Atlantic coasts of France. The analysis of nuclear and mitochondrial loci provided contrasting results. The mtDNA revealed two genetically distinct groups, one closely related to the Iberian and African harbour porpoises, and the second related to individuals from the more northern waters of Europe. In contrast, nuclear polymorphisms did not display such a distinction. Nuclear markers suggested that harbour porpoises behaved as a randomly mating population along the Atlantic coasts of France. The difference between the two kinds of markers can be explained by differences in their mode of inheritance, the mtDNA being maternally inherited in contrast to nuclear loci that are bi-parentally inherited. Our results provide evidence that a major proportion of the animals we sampled are admixed individuals from the two genetically distinct populations previously identified along the Iberian coasts and in the North East Atlantic. The French Atlantic coasts are clearly the place where these two previously separated populations of harbour porpoises are now admixing. The present shifts in distribution of harbour porpoises along this coast is likely caused by habitat changes that will need to be further studied.

Mitogenomic sequences better resolve stock structure of southern Greater Caribbean green turtle rookeries

Analyses of mitochondrial control region polymorphisms have supported the presence of several demographically independent green turtle (Chelonia mydas) rookeries in the Greater Caribbean region. However, extensive sharing of common haplotypes based on 490-bp control region sequences confounds assessment of the scale of natal homing and population structure among regional rookeries. We screened the majority of the mitochondrial genomes of 20 green turtles carrying the common haplotype CM-A5 and representing the rookeries of Buck Island, St. Croix, United States Virgin Islands (USVI); Aves Island, Venezuela; Galibi, Suriname; and Tortuguero, Costa Rica. Five single-nucleotide polymorphisms (SNPs) were identified that subdivided CM-A5 among regions. Mitogenomic pairwise φ(ST) values of eastern Caribbean rookery comparisons were markedly lower than the respective pairwise F(ST) values. This discrepancy results from the presence of haplotypes representing two divergent lineages in each rookery, highlighting the importance of choosing the appropriate test statistic for addressing the study question. Haplotype frequency differentiation supports demographic independence of Aves Island and Suriname, emphasizing the need to recognize the smaller Aves rookery as a distinct management unit. Aves Island and Buck Island rookeries shared mitogenomic haplotypes; however, frequency divergence suggests that the Buck Island rookery is sufficiently demographically isolated to warrant management unit status for the USVI rookeries. Given that haplotype sharing among rookeries is common in marine turtles with cosmopolitan distributions, mitogenomic sequencing may enhance inferences of population structure and phylogeography, as well as improve the resolution of mixed stock analyses aimed at estimating natal origins of foraging turtles.

Temporal trends and spatial differences of perfluoroalkylated substances in livers of harbor porpoise (Phocoena phocoena) populations from Northern Europe, 1991-2008

Long-term temporal trends (1991-2008) and spatial differences of perfluoroalkylated substances (PFASs) were investigated in harbor porpoise (Phocoena phocoena) liver samples of juvenile females from the Baltic and North Sea. Additionally, spatial differences between the populations in the Baltic Sea and Atlantic Ocean (i.e. Iceland and Norway) and the influence of the body mass, age and sex on the PFAS concentrations were examined. Perfluorooctane sulfonate (PFOS) was the predominant compound with a concentration range of 160-2425 ng/g wet weight (ww), followed by perfluorooctane sulfonamide (PFOSA, 1-237 ng/g ww) and perfluoroundecanoate (PFUnA, 3-124 ng/g ww). In terms of temporal trends, perfluoroalkyl sulfonates (PFSAs) and PFOSA concentrations decreased over time, while, conversely, the C(9)-C(13) perfluoroalkyl carboxylate (PFCA) concentrations increased. Spatial distribution of the contaminant concentrations showed consistently higher concentrations in the Baltic Sea and lowest concentrations in the Icelandic population of the Atlantic Ocean.

Population genetic structure and social kin associations of franciscana dolphin, Pontoporia blainvillei

We investigated population and social structure of the franciscana dolphin, Pontoporia blainvillei, an endemic and the most endangered cetacean of the southwestern Atlantic Ocean. We analyzed samples from the Rio de la Plata estuary obtained in Uruguayan waters and from the Atlantic Ocean obtained in both Uruguayan and Brazilian waters. Mitochondrial and microsatellite DNA markers were used to study differentiation between the estuary and the ocean and the association between kinship and social group structure. Although multilocus analyses suggested that franciscanas are structured into 2 subpopulations (K = 2, divergence among clusters: F(ST) = 0.06, P = 0.002; R(ST) = 0.3, P = 0.001), mitochondrial markers did not support such divergence (F(ST) = 0.02, P = 0.12; Ф(ST) = 0.06, P =0.06). However, these units are not entirely segregated geographically. Regarding social structure, some groups are composed by first-order related individuals (R ≥ 0.5, P < 0.5). Overall, the data suggest that matrilines could be the social unit in this species. We argue that the divergence found could be associated to local adaptation and social structure, resulting from either feature leading to a recent divergence or reflecting equilibrium between local differentiation and gene flow. This evidence supports considering franciscanas from the Rio de la Plata estuary a discrete management unit.

Mitochondrial control region haplotypes of the South American sea lion Otaria flavescens (Shaw, 1800)

The South American sea lion, Otaria flavescens, is widely distributed along the Pacific and Atlantic coasts of South America. However, along the Brazilian coast, there are only two nonbreeding sites for the species (Refúgio de Vida Silvestre da Ilha dos Lobos and Refúgio de Vida Silvestre do Molhe Leste da Barra do Rio Grande), both in Southern Brazil. In this region, the species is continuously under the effect of anthropic activities, mainly those related to environmental contamination with organic and inorganic chemicals and fishery interactions. This paper reports, for the first time, the genetic diversity of O. flavescens found along the Southern Brazilian coast. A 287-bp fragment of the mitochondrial DNA control region (D-loop) was analyzed. Seven novel haplotypes were found in 56 individuals (OFA1-OFA7), with OFA1 being the most frequent (47.54%). Nucleotide diversity was moderate (π = 0.62%) and haplotype diversity was relatively low (67%). Furthermore, the median joining network analysis indicated that Brazilian haplotypes formed a reciprocal monophyletic clade when compared to the haplotypes from the Peruvian population on the Pacific coast. These two populations do not share haplotypes and may have become isolated some time back. Further genetic studies covering the entire species distribution are necessary to better understand the biological implications of the results reported here for the management and conservation of South American sea lions.

Isolation by environmental distance in mobile marine species: molecular ecology of franciscana dolphins at their southern range

The assessment of population structure is a valuable tool for studying the ecology of endangered species and drafting conservation strategies. As we enhance our understanding about the structuring of natural populations, it becomes important that we also understand the processes behind these patterns. However, there are few rigorous assessments of the influence of environmental factors on genetic patterns in mobile marine species. Given their dispersal capabilities and localized habitat preferences, coastal cetaceans are adequate study species for evaluating environmental effects on marine population structure. The franciscana dolphin, a rare coastal cetacean endemic to the Western South Atlantic, was studied to examine these issues. We analysed genetic data from the mitochondrial DNA and 12 microsatellite markers for 275 franciscana samples utilizing frequency-based, maximum-likelihood and Bayesian algorithms to assess population structure and migration patterns. This information was combined with 10 years of remote sensing environmental data (chlorophyll concentration, water turbidity and surface temperature). Our analyses show the occurrence of genetically isolated populations within Argentina, in areas that are environmentally distinct. Combined evidence of genetic and environmental structure suggests that isolation by distance and a process here termed isolation by environmental distance can explain the observed correlations. Our approach elucidated important ecological and conservation aspects of franciscana dolphins, and has the potential to increase our understanding of ecological processes influencing genetic patterns in other marine species.