Recent demographic history and present fine-scale structure in the Northwest Atlantic leatherback (Dermochelys coriacea) turtle population

Natal Origin and Spatiotemporal Distribution of Leatherback Turtle (Dermochelys coriacea) Strandings at a Foraging Hotspot in Temperate Waters of the Southwest Atlantic Ocean

Leatherback turtles migrate long distances between nesting beaches and distant foraging areas worldwide. This study analyzes the genetic diversity, life history stage, spatiotemporal distribution, and associated threats of a foraging aggregation in the Southwest Atlantic Ocean. A total of 242 leatherbacks stranded or bycaught by artisanal fisheries were recorded from 1997 to 2021 in Uruguay, with sizes ranging from 110.0 to 170.0 cm carapace lengths, indicating that the aggregation is composed of large juveniles and adults. Results of Bayesian mixed-stock analysis show that leatherbacks come primarily from the West African rookeries, based on mitochondrial DNA sequences obtained from 59 of the turtles representing seven haplotypes, including a novel one (Dc1.7). The main threat identified in the area is the fisheries bycatch but most of the carcasses observed were badly decomposed. There was significant seasonal and interannual variability in strandings that is likely associated with the availability of prey and the intensity of the fishing effort. Taken together, these findings reinforce the importance of these South American foraging areas for leatherbacks and the need to determine regional habitat use and migratory routes across the broader Atlantic region, in order to develop effective conservation measures to mitigate threats both at nesting beaches and foraging areas.

Global phylogeography of ridley sea turtles (Lepidochelys spp.): evolution, demography, connectivity, and conservation

Globally distributed marine taxa are well suited for investigations of biogeographic impacts on genetic diversity, connectivity, and population demography. The sea turtle genus Lepidochelys includes the wide-ranging and abundant olive ridley (L. olivacea), and the geographically restricted and ‘Critically Endangered’ Kemp’s ridley (L. kempii). To investigate their historical biogeography, we analyzed a large dataset of mitochondrial DNA (mtDNA) sequences from olive (n = 943) and Kemp’s (n = 287) ridleys, and genotyped 15 nuclear microsatellite loci in a global sample of olive ridleys (n = 285). We found that the ridley species split ~ 7.5 million years ago, before the Panama Isthmus closure. The most ancient mitochondrial olive ridley lineage, located in the Indian Ocean, was dated to ~ 2.2 Mya. Both mitochondrial and nuclear markers revealed significant structure for olive ridleys between Atlantic (ATL), East Pacific (EP), and Indo-West Pacific (IWP) areas. However, the divergence of mtDNA clades was very recent (< 1 Mya) with low within- clade diversity, supporting a recurrent extinction-recolonization model for these ocean regions. All data showed that ATL and IWP groups were more closely related than those in the EP, with mtDNA data supporting recent recolonization of the ATL from the IWP. Individual olive ridley dispersal between the ATL, EP, and IN/IWP could be interpreted as more male- than female-biased, and genetic diversity was lowest in the Atlantic Ocean. All populations showed signs of recent expansion, and estimated time frames were concordant with their recent colonization history. Investigating species abundance and distribution changes over time is central to evolutionary biology, and this study provides a historical biogeographic context for marine vertebrate conservation and management.

Divergence and hybridization in sea turtles: Inferences from genome data show evidence of ancient gene flow between species

Reconstructing past events of hybridization and population size changes are required to understand speciation mechanisms and current patterns of genetic diversity, and ultimately contribute to species' conservation. Sea turtles are ancient species currently facing anthropogenic threats including climate change, fisheries, and illegal hunting. Five of the seven extant sea turtle species are known to currently hybridize, especially along the Brazilian coast where some populations can have ~32%-42% of hybrids. Although frequently observed today, it is not clear what role hybridization plays in the evolutionary diversification of this group of reptiles. In this study, we generated whole genome resequencing data of the five globally distributed sea turtle species to estimate a calibrated phylogeny and the population size dynamics, and to understand the role of hybridization in shaping the genomes of these ancient species. Our results reveal discordant species divergence dates between mitochondrial and nuclear genomes, with a high frequency of conflicting trees throughout the nuclear genome suggesting that some sea turtle species frequently hybridized in the past. The reconstruction of the species' demography showed a general decline in effective population sizes with no signs of recovery, except for the leatherback sea turtle. Furthermore, we discuss the influence of reference bias in our estimates. We show long-lasting ancestral gene flow events within Chelonioidea that continued for millions of years after initial divergence. Speciation with gene flow is a common pattern in marine species, and it raises questions whether current hybridization events should be considered as a part of these species' evolutionary history or a conservation issue.

New Genetic Insights About Hybridization and Population Structure of Hawksbill and Loggerhead Turtles From Brazil

An extremely high incidence of hybridization among sea turtles is found along the Brazilian coast. This atypical phenomenon and its impact on sea turtle conservation can be elucidated through research focused on the evolutionary history of sea turtles. We assessed high-quality multilocus haplotypes of 143 samples of the 5 species of sea turtles that occur along the Brazilian coast to investigate the hybridization process and the population structure of hawksbill (Eretmochelys imbricata) and loggerhead turtles (Caretta caretta). The multilocus data were initially used to characterize interspecific hybrids. Introgression (F2 hybrids) was only confirmed in hatchlings of F1 hybrid females (hawksbill × loggerhead), indicating that introgression was either previously overestimated and F2 hybrids may not survive to adulthood, or the first-generation hybrid females nesting in Brazil were born as recent as few decades ago. Phylogenetic analyses using nuclear markers recovered the mtDNA-based Indo-Pacific and Atlantic lineages for hawksbill turtles, demonstrating a deep genetic divergence dating from the early Pliocene. In addition, loggerhead turtles that share a common feeding area and belong to distinct Indo-Pacific and Atlantic mtDNA clades present no clear genetic differentiation at the nuclear level. Finally, our results indicate that hawksbill and loggerhead rookeries along the Brazilian coast are likely connected by male-mediated gene flow.

Survival and breeding interval of an endangered marine vertebrate, the leatherback turtle Dermochelys coriacea, in French Guiana

In the context of global change, endangered species such as sea turtles undergo strong population dynamics changes. Understanding demographic processes inducing such changes is critical for developing appropriate measures for conservation and management. Nesting females of the French Guiana population of leatherback sea turtles Dermochelys coriacea represented 40% of the world’s total in 2001; today, they represent only 10%. Here, we analyze data from the long-term monitoring program of leatherbacks in French Guiana, based on 28 yr of capture-mark-recapture data from 46051 individuals in northwestern French Guiana. We used multievent models (multistate capture-recapture models with state uncertainty) to represent the main peculiarity of the life cycle, intermittent reproduction, and to take advantage of the use of several different types of marks to account for mark loss and incomplete detection. The average annual adult survival probability (mean ± SE) 0.789 ± 0.009 is low compared to published estimates for this species. By combining the estimates of departure and return probabilities, we provide an estimate of the interval among laying seasons equal to 2.777 ± 0.118 yr, which matches previous findings. Taking survival into account, this interval translates into an average of 1.704 ± 0.034 further reproductive seasons for an adult female having just bred, which is remarkably low compared to other long-lived vertebrates. The representation of the life cycle and mark loss in our models provides a sound structure for broader and richer analyses to extract biological information from marked individuals and offers an invaluable source of demographic information on leatherbacks, a species for which only a small segment of the population is accessible to intermittent sampling.

Genetic population structure and demography of an apex predator, the tiger shark Galeocerdo cuvier

Population genetics has been increasingly applied to study large sharks over the last decade. Whilst large shark species are often difficult to study with direct methods, improved knowledge is needed for both population management and conservation, especially for species vulnerable to anthropogenic and climatic impacts. The tiger shark, Galeocerdo cuvier, is an apex predator known to play important direct and indirect roles in tropical and subtropical marine ecosystems. While the global and Indo-West Pacific population genetic structure of this species has recently been investigated, questions remain over population structure and demographic history within the western Indian (WIO) and within the western Pacific Oceans (WPO). To address the knowledge gap in tiger shark regional population structures, the genetic diversity of 286 individuals sampled in seven localities was investigated using 27 microsatellite loci and three mitochondrial genes (CR,COI, and cytb). A weak genetic differentiation was observed between the WIO and the WPO, suggesting high genetic connectivity. This result agrees with previous studies and highlights the importance of the pelagic behavior of this species to ensure gene flow. Using approximate Bayesian computation to couple information from both nuclear and mitochondrial markers, evidence of a recent bottleneck in the Holocene (2,000-3,000 years ago) was found, which is the most probable cause for the low genetic diversity observed. A contemporary effective population size as low as 111 [43,369] was estimated during the bottleneck. Together, these results indicate low genetic diversity that may reflect a vulnerable population sensitive to regional pressures. Conservation measures are thus needed to protect a species that is classified as Near Threatened.

Demographic inference through approximate-Bayesian-computation skyline plots

The skyline plot is a graphical representation of historical effective population sizes as a function of time. Past population sizes for these plots are estimated from genetic data, without a priori assumptions on the mathematical function defining the shape of the demographic trajectory. Because of this flexibility in shape, skyline plots can, in principle, provide realistic descriptions of the complex demographic scenarios that occur in natural populations. Currently, demographic estimates needed for skyline plots are estimated using coalescent samplers or a composite likelihood approach. Here, we provide a way to estimate historical effective population sizes using an Approximate Bayesian Computation (ABC) framework. We assess its performance using simulated and actual microsatellite datasets. Our method correctly retrieves the signal of contracting, constant and expanding populations, although the graphical shape of the plot is not always an accurate representation of the true demographic trajectory, particularly for recent changes in size and contracting populations. Because of the flexibility of ABC, similar approaches can be extended to other types of data, to multiple populations, or to other parameters that can change through time, such as the migration rate.

Full circumpolar migration ensures evolutionary unity in the Emperor penguin

Defining reliable demographic models is essential to understand the threats of ongoing environmental change. Yet, in the most remote and threatened areas, models are often based on the survey of a single population, assuming stationarity and independence in population responses. This is the case for the Emperor penguin Aptenodytes forsteri, a flagship Antarctic species that may be at high risk continent-wide before 2100. Here, using genome-wide data from the whole Antarctic continent, we reveal that this top-predator is organized as one single global population with a shared demography since the late Quaternary. We refute the view of the local population as a relevant demographic unit, and highlight that (i) robust extinction risk estimations are only possible by including dispersal rates and (ii) colony-scaled population size is rather indicative of local stochastic events, whereas the species' response to global environmental change is likely to follow a shared evolutionary trajectory.

Delimiting populations is crucial for conserving threatened species. Using genome-wide data from the whole of Antarctica, Cristofari et al. show that Emperor penguins are organised into a single global population that have shared demography since the late Quarternary.

Mixed-stock analysis in green turtles Chelonia mydas: mtDNA decipher current connections among west Atlantic populations

The green turtle Chelonia mydas undertakes wide-ranging migrations between feeding and nesting sites, resulting in mixing and isolation of genetic stocks. We used mtDNA control region to characterize the genetic composition, population structure, and natal origins of C. mydas in the West Atlantic Ocean, at one feeding ground (State of Rio de Janeiro, Brazil), and three Caribbean nesting grounds (French Guiana, Guadeloupe, and Suriname). The feeding ground presented considerable frequency of common haplotypes from the South Atlantic, whereas the nesting sites presented a major contribution of the most common haplotype from the Caribbean. MSA revealed multiple origins of individuals at the feeding ground, notably from Ascension Island, Guinea Bissau, and French Guiana. This study enables a better understanding of the dispersion patterns and highlights the importance of connecting both nesting and feeding areas. Effective conservation initiatives need to encompass these ecologically and geographically distinct sites as well as those corridors connecting them.

Phylogeography, Genetic Diversity, and Management Units of Hawksbill Turtles in the Indo-Pacific

Hawksbill turtle (Eretmochelys imbricata) populations have experienced global decline because of a history of intense commercial exploitation for shell and stuffed taxidermied whole animals, and harvest for eggs and meat. Improved understanding of genetic diversity and phylogeography is needed to aid conservation. In this study, we analyzed the most geographically comprehensive sample of hawksbill turtles from the Indo-Pacific Ocean, sequencing 766 bp of the mitochondrial control region from 13 locations (plus Aldabra, n = 4) spanning over 13500 km. Our analysis of 492 samples revealed 52 haplotypes distributed in 5 divergent clades. Diversification times differed between the Indo-Pacific and Atlantic lineages and appear to be related to the sea-level changes that occurred during the Last Glacial Maximum. We found signals of demographic expansion only for turtles from the Persian Gulf region, which can be tied to a more recent colonization event. Our analyses revealed evidence of transoceanic migration, including connections between feeding grounds from the Atlantic Ocean and Indo-Pacific rookeries. Hawksbill turtles appear to have a complex pattern of phylogeography, showing a weak isolation by distance and evidence of multiple colonization events. Our novel dataset will allow mixed-stock analyses of hawksbill turtle feeding grounds in the Indo-Pacific by providing baseline data needed for conservation efforts in the region. Eight management units are proposed in our study for the Indo-Pacific region that can be incorporated in conservation plans of this critically endangered species.

Contextualising the Last Survivors: Population Structure of Marine Turtles in the Dominican Republic

Nesting by three species of marine turtles persists in the Dominican Republic, despite historic threats and long-term population decline. We conducted a genetic survey of marine turtles in the Dominican Republic in order to link them with other rookeries around the Caribbean. We sequenced a 740bp fragment of the control region of the mitochondrial DNA of 92 samples from three marine turtle species [hawksbill (n = 48), green (n = 2) and leatherback (n = 42)], and incorporated published data from other nesting populations and foraging grounds. The leatherback turtle (Dermochelys coriacea) in the Dominican Republic appeared to be isolated from Awala-Yalimapo, Cayenne, Trinidad and St. Croix but connected with other Caribbean populations. Two distinct nesting populations of hawksbill turtles (Eremochelys imbricata) were detected in the Dominican Republic and exhibited interesting patterns of connectivity with other nesting sites and juvenile and adult male foraging aggregations. The green sea turtle (Chelonia mydas) has almost been extirpated from the Dominican Republic and limited inference could be made from our samples. Finally, results were compared with Lagrangian drifting buoys and published Lagrangian virtual particles that travelled through the Dominican Republic and Caribbean waters. Conservation implications of sink-source effects or genetic isolation derived from these complex inter-connections are discussed for each species and population.