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. Supplementary Information: The online version contains supplementary material available at 10.1007/s10592-022-01465-3.

Quantifying local ecological knowledge to model historical abundance of long-lived, heavily-exploited fauna.

Deriving robust historical population trends for long-lived species subject to human exploitation is challenging in scenarios where long-term scientific data are scarce or unavailable, as often occurs for species affected by small-scale fisheries and subsistence hunting. The importance of Local Ecological Knowledge (LEK) in data-poor scenarios is increasingly recognized in conservation, both in terms of uncovering historical trends and for engaging community stewardship of historic information. Building on previous work in marine historical ecology and local ecological knowledge, we propose a mixed socio-ecological framework to reliably document and quantify LEK to reconstruct historical population trends. Our method can be adapted by interdisciplinary teams to study various long-lived taxa with a history of human use. We demonstrate the validity of our approach by reconstructing long-term abundance data for the heavily-exploited East Pacific green turtle (Chelonia mydas) in Baja California, Mexico, which was driven to near extinction by a largely unregulated fishery from the early 1950s to the 1980s. No scientific baseline abundance data were available for this time-frame because recent biological surveys started in 1995 after all green turtle fisheries in the area were closed. To fill this data gap, we documented LEK among local fishers using ethnographic methods and obtained verified, qualitative data to understand the socio-environmental complexity of the green turtle fishery. We then established an iterative framework to synthesize and quantify LEK using generalized linear models (GLMs) and nonlinear regression (NLR) to generate a standardized, LEK-derived catch-per-unit-effort (CPUE) time-series. CPUE is an index of abundance that is compatible with contemporary scientific survey data. We confirmed the accuracy of LEK-derived CPUE estimates via comparisons with fisheries statistics available for 1962-1982. We then modeled LEK-derived abundance trends prior to 1995 using NLR. Our model established baseline abundance and described historical declines, revealing that the most critical (exponential) decline occurred between 1960 and 1980. This robust integration of LEK data with ecological science is of critical value for conservation and management, as it contributes to a holistic view of a species' historic and contemporary conservation status.

Recent advances on the estimation of the thermal reaction norm for sex ratios.

Temperature-dependent sex determination, or TSD, is a widespread phenomenon in reptiles. The shape of the relationship between constant incubation temperature and sex ratio defines the TSD pattern. The TSD pattern is considered a life-history parameter important for conservation because the wider the range of temperatures producing both sexes, the more resilient the species is to climate change impacts. We review the different published equations and methodologies that have been used to model TSD patterns. We describe a new flexible model that allows for an asymmetrical pattern around the pivotal temperature, which is the constant temperature producing both sexes in equal proportions. We show that Metropolis-Hastings with Markov chain produced by a Monte Carlo process has many advantages compared to maximum likelihood and is preferred. Finally, we apply the models to results from incubation experiments using eggs from the marine turtle Lepidochelys olivacea originating in Northeast Indian, East Pacific, and West Atlantic Regional Management Units (RMUs) and find large differences in pivotal temperatures but not in transitional ranges of temperatures.

Deepwater Horizon oil spill impacts on sea turtles could span the Atlantic.

We investigated the extent that the 2010 Deepwater Horizon oil spill potentially affected oceanic-stage sea turtles from populations across the Atlantic. Within an ocean-circulation model, particles were backtracked from the Gulf of Mexico spill site to determine the probability of young turtles arriving in this area from major nesting beaches. The abundance of turtles in the vicinity of the oil spill was derived by forward-tracking particles from focal beaches and integrating population size, oceanic-stage duration and stage-specific survival rates. Simulations indicated that 321 401 (66 199-397 864) green (Chelonia mydas), loggerhead (Caretta caretta) and Kemp's ridley (Lepidochelys kempii) turtles were likely within the spill site. These predictions compared favourably with estimates from in-water observations recently made available to the public (though our initial predictions for Kemp's ridley were substantially lower than in-water estimates, better agreement was obtained with modifications to mimic behaviour of young Kemp's ridley turtles in the northern Gulf). Simulations predicted 75.2% (71.9-76.3%) of turtles came from Mexico, 14.8% (11-18%) from Costa Rica, 5.9% (4.8-7.9%) from countries in northern South America, 3.4% (2.4-3.5%) from the United States and 1.6% (0.6-2.0%) from West African countries. Thus, the spill's impacts may extend far beyond the current focus on the northern Gulf of Mexico.

Does polyandry really pay off? The effects of multiple mating and number of fathers on morphological traits and survival in clutches of nesting green turtles at Tortuguero.

Despite the long debate of whether or not multiple mating benefits the offspring, studies still show contradictory results. Multiple mating takes time and energy. Thus, if females fertilize their eggs with a single mating, why to mate more than once? We investigated and inferred paternal identity and number of sires in 12 clutches (240 hatchlings) of green turtles (Chelonia mydas) nests at Tortuguero, Costa Rica. Paternal alleles were inferred through comparison of maternal and hatchling genotypes, and indicated multiple paternity in at least 11 of the clutches (92%). The inferred average number of fathers was three (ranging from 1 to 5). Moreover, regression analyses were used to investigate for correlation of inferred clutch paternity with morphological traits of hatchlings fitness (emergence success, length, weight and crawling speed), the size of the mother, and an environmental variable (incubation temperature). We suggest and propose two different comparative approaches for evaluating morphological traits and clutch paternity, in order to infer greater offspring survival. First, clutches coded by the exact number of fathers and second by the exact paternal contribution (fathers who gives greater proportion of the offspring per nest). We found significant differences (P < 0.05) in clutches coded by the exact number of fathers for all morphological traits. A general tendency of higher values in offspring sired by two to three fathers was observed for the length and weight traits. However, emergence success and crawling speed showed different trends which unable us to reach any further conclusion. The second approach analysing the paternal contribution showed no significant difference (P > 0.05) for any of the traits. We conclude that multiple paternity does not provide any extra benefit in the morphological fitness traits or the survival of the offspring, when analysed following the proposed comparative statistical methods.

First approximation to congenital malformation rates in embryos and hatchlings of sea turtles.

BACKGROUND: Congenital malformations in sea turtles have been considered sporadical. Research carried out in the Mexican Pacific revealed high levels of congenital malformations in the olive ridley, but little or no information is available for other species. We present results from analyses of external congenital malformations in olive ridley, green, and hawskbill sea turtles from Mexican rookeries on the Pacific coast and Gulf of Mexico. METHODS: We examined 150 green and hawksbill nests and 209 olive ridley nests during the 2010 and 2012 nesting seasons, respectively. Olive ridley eggs were transferred to a hatchery and incubated in styrofoam boxes. Nests from the other two species were left in situ. Number of eggs, live and dead hatchlings, and eggs with or without embryonic development were registered. Malformation frequency was evaluated with indices of prevalence and severity. RESULTS: Mortality levels, prevalence and severity were higher in olive ridley than in hawksbill and green sea turtles. Sixty-three types of congenital malformations were observed in embryos, and dead or live hatchlings. Of these, 38 are new reports; 35 for wild sea turtles, three for vertebrates. Thirty-one types were found in hawksbill, 23 in green, and 59 in olive ridley. The head region showed a higher number of malformation types. Malformation levels in the olive ridley were higher than previously reported. CONCLUSION: Olive ridleys seem more prone to the occurrence of congenital malformations than the other two species. Whether the observed malformation levels are normal or represent a health problem cannot be currently ascertained without long-term assessments.

Geographic patterns of genetic variation in a broadly distributed marine vertebrate: new insights into loggerhead turtle stock structure from expanded mitochondrial DNA sequences.

Previous genetic studies have demonstrated that natal homing shapes the stock structure of marine turtle nesting populations. However, widespread sharing of common haplotypes based on short segments of the mitochondrial control region often limits resolution of the demographic connectivity of populations. Recent studies employing longer control region sequences to resolve haplotype sharing have focused on regional assessments of genetic structure and phylogeography. Here we synthesize available control region sequences for loggerhead turtles from the Mediterranean Sea, Atlantic, and western Indian Ocean basins. These data represent six of the nine globally significant regional management units (RMUs) for the species and include novel sequence data from Brazil, Cape Verde, South Africa and Oman. Genetic tests of differentiation among 42 rookeries represented by short sequences (380 bp haplotypes from 3,486 samples) and 40 rookeries represented by long sequences (∼800 bp haplotypes from 3,434 samples) supported the distinction of the six RMUs analyzed as well as recognition of at least 18 demographically independent management units (MUs) with respect to female natal homing. A total of 59 haplotypes were resolved. These haplotypes belonged to two highly divergent global lineages, with haplogroup I represented primarily by CC-A1, CC-A4, and CC-A11 variants and haplogroup II represented by CC-A2 and derived variants. Geographic distribution patterns of haplogroup II haplotypes and the nested position of CC-A11.6 from Oman among the Atlantic haplotypes invoke recent colonization of the Indian Ocean from the Atlantic for both global lineages. The haplotypes we confirmed for western Indian Ocean RMUs allow reinterpretation of previous mixed stock analysis and further suggest that contemporary migratory connectivity between the Indian and Atlantic Oceans occurs on a broader scale than previously hypothesized. This study represents a valuable model for conducting comprehensive international cooperative data management and research in marine ecology.

Nesting ecology of Chelonia mydas Testudines: Cheloniidae on the Guanahacabibes Peninsula, Cuba

The nesting colony of green sea turtles Chelonia mydas at Guanahacabibes Peninsula Biosphere Reserve and National Park is one of the largest in the Cuban archipelago; however, little information about its nesting ecology is available. Temporal and spatial variation in nesting and reproductive success as well as morphometric characteristics of gravid females were used to ecologically characterize this colony. Nine beaches of the Southernmost coast of Guanahacabibes Peninsula were monitored for 14 years 1998-2012 to determine green turtle nesting activity, from May to September peak nesting season in this area. Beach dimensions were measured to determine nest density using the length and the area. Afterward the beaches were divided in two categories, index and secondary. Females were measured and tagged to compare new tagged females 823 with returning tagged females 140. Remigration interval was also determined. Temporal variation was identified as the annual number of nesting emergences and oviposits per female, with apparent peaks in reproductive activity on a biennial cycle in the first six years followed by periods of annual increase in nest number 2003-2008 and periods of decreasing number of nests 2010-2012. We also found intra-seasonal variation with the highest nesting activity in July, particularly in the second half of the month. The peak emergence time was 22:00-02:00hr. In terms of spatial variation, smaller beaches had the highest nest density and nesting was more frequent 6-9m from the high tide line, where hatchling production was maximized although hatchling success was high on average, above 80. Morphometric analysis of females was made and newly tagged turtles were smaller on average than remigrants. Our results are only a first attempt at characterizing Guanahacabibes populations but have great value for establishing conservation priorities within the context of national management plans, and for efficient monitoring and protection of nesting beaches.

La colonia de Tortuga verde Chelonia mydas que anida en la Península de Guanahacabibeses es una de las más grandes del archipiélago cubano pero existe poca información disponible sobre la misma. Por ello utilizamos la variación espacial y temporal de la anidación así como el éxito reproductivo y las características morfométricas de las hembras grávidas para realizar la primera caracterización ecológica de esta colonia. Se estudiaron nueve playas durante 14 años 1998-2012 desde mayo hasta septiembre pico de anidación. Estas se dividieron en dos categorías: índices y secundarias. Las hembras fueron medidas y marcadas para comparar las nuevas hembras marcadas 823 con las remigrantes 140. Se detectó una variación temporal en la anidación a tres niveles: interanual, dentro de la temporada y durante la noche. Se observa un ciclo bienal en los primeros seis años, luego un incremento anual 2003-2008 seguido por un período de decrecimiento en las anidaciones 2010-2012. Dentro de la temporada se encontró un pico en el mes de julio aunque cambia de la primera quincena a la segunda en años diferentes y la hora de emergencia preferentemente de 22:00-02:00hr. En la variación espacial, las playas pequeñas tuvieron mayor densidad de anidación, especialmente entre los 6-9m hasta la línea media de marea alta, donde se maximiza la producción de neonatos, aunque el éxito de emergencia de los neonatos se mantuvo como promedio por encima del 80 en toda la playa. Morfométricamente las hembras remigrantes resultaron mayores como promedio que las marcadas por primera vez. Estos resultados son una primera aproximación para la caracterización de la colonia de Guanahacabibes pero tienen gran valor para establecer prioridades de conservación dentro del contexto de los planes de manejo nacionales así como para garantizar un monitoreo eficiente y la protección de las playas de anidación.

Nesting ecology of Chelonia mydas (Testudines: Cheloniidae) on the Guanahacabibes Peninsula, Cuba.

The nesting colony of green sea turtles (Chelonia mydas) at Guanahacabibes Peninsula Biosphere Reserve and National Park is one of the largest in the Cuban archipelago; however, little information about its nesting ecology is available. Temporal and spatial variation in nesting and reproductive success as well as morphometric characteristics of gravid females were used to ecologically characterize this colony. Nine beaches of the Southernmost coast of Guanahacabibes Peninsula were monitored for 14 years (1998-2012) to determine green turtle nesting activity, from May to September (peak nesting season in this area). Beach dimensions were measured to determine nest density using the length and the area. Afterward the beaches were divided in two categories, index and secondary. Females were measured and tagged to compare new tagged females (823) with returning tagged females (140). Remigration interval was also determined. Temporal variation was identified as the annual number of nesting emergences and oviposits per female, with apparent peaks in reproductive activity on a biennial cycle in the first six years followed by periods of annual increase in nest number (2003-2008) and periods of decreasing number of nests (2010-2012). We also found intra-seasonal variation with the highest nesting activity in July, particularly in the second half of the month. The peak emergence time was 22:00-02:00 hr. In terms of spatial variation, smaller beaches had the highest nest density and nesting was more frequent 6-9m from the high tide line, where hatchling production was maximized although hatchling success was high on average, above 80%. Morphometric analysis of females was made and newly tagged turtles were smaller on average than remigrants. Our results are only a first attempt at characterizing Guanahacabibes' populations but have great value for establishing conservation priorities within the context of national management plans, and for efficient monitoring and protection of nesting beaches.

Re-examination of population structure and phylogeography of hawksbill turtles in the wider Caribbean using longer mtDNA sequences.

Management of the critically endangered hawksbill turtle in the Wider Caribbean (WC) has been hampered by knowledge gaps regarding stock structure. We carried out a comprehensive stock structure re-assessment of 11 WC hawksbill rookeries using longer mtDNA sequences, larger sample sizes (N = 647), and additional rookeries compared to previous surveys. Additional variation detected by 740 bp sequences between populations allowed us to differentiate populations such as Barbados-Windward and Guadeloupe (F (st) = 0.683, P < 0.05) that appeared genetically indistinguishable based on shorter 380 bp sequences. POWSIM analysis showed that longer sequences improved power to detect population structure and that when N < 30, increasing the variation detected was as effective in increasing power as increasing sample size. Geographic patterns of genetic variation suggest a model of periodic long-distance colonization coupled with region-wide dispersal and subsequent secondary contact within the WC. Mismatch analysis results for individual clades suggest a general population expansion in the WC following a historic bottleneck about 100 000-300 000 years ago. We estimated an effective female population size (N (ef)) of 6000-9000 for the WC, similar to the current estimated numbers of breeding females, highlighting the importance of these regional rookeries to maintaining genetic diversity in hawksbills. Our results provide a basis for standardizing future work to 740 bp sequence reads and establish a more complete baseline for determining stock boundaries in this migratory marine species. Finally, our findings illustrate the value of maintaining an archive of specimens for re-analysis as new markers become available.

Multiple distant origins for green sea turtles aggregating off Gorgona Island in the Colombian eastern Pacific.

Mitochondrial DNA analyses have been useful for resolving maternal lineages and migratory behavior to foraging grounds (FG) in sea turtles. However, little is known about source rookeries and haplotype composition of foraging green turtle aggregations in the southeastern Pacific. We used mitochondrial DNA control region sequences to identify the haplotype composition of 55 green turtles, Chelonia mydas, captured in foraging grounds of Gorgona National Park in the Colombian Pacific. Amplified fragments of the control region (457 bp) revealed the presence of seven haplotypes, with haplotype (h) and nucleotide (π) diversities of h = 0.300±0.080 and π = 0.009±0.005 respectively. The most common haplotype was CMP4 observed in 83% of individuals, followed by CMP22 (5%). The genetic composition of the Gorgona foraging population primarily comprised haplotypes that have been found at eastern Pacific rookeries including Mexico and the Galapagos, as well as haplotypes of unknown stock origin that likely originated from more distant western Pacific rookeries. Mixed stock analysis suggests that the Gorgona FG population is comprised mostly of animals from the Galapagos rookery (80%). Lagrangian drifter data showed that movement of turtles along the eastern Pacific coast and eastward from distant western and central Pacific sites was possible through passive drift. Our results highlight the importance of this protected area for conservation management of green turtles recruited from distant sites along the eastern Pacific Ocean.

Global conservation priorities for marine turtles.

Where conservation resources are limited and conservation targets are diverse, robust yet flexible priority-setting frameworks are vital. Priority-setting is especially important for geographically widespread species with distinct populations subject to multiple threats that operate on different spatial and temporal scales. Marine turtles are widely distributed and exhibit intra-specific variations in population sizes and trends, as well as reproduction and morphology. However, current global extinction risk assessment frameworks do not assess conservation status of spatially and biologically distinct marine turtle Regional Management Units (RMUs), and thus do not capture variations in population trends, impacts of threats, or necessary conservation actions across individual populations. To address this issue, we developed a new assessment framework that allowed us to evaluate, compare and organize marine turtle RMUs according to status and threats criteria. Because conservation priorities can vary widely (i.e. from avoiding imminent extinction to maintaining long-term monitoring efforts) we developed a "conservation priorities portfolio" system using categories of paired risk and threats scores for all RMUs (n = 58). We performed these assessments and rankings globally, by species, by ocean basin, and by recognized geopolitical bodies to identify patterns in risk, threats, and data gaps at different scales. This process resulted in characterization of risk and threats to all marine turtle RMUs, including identification of the world's 11 most endangered marine turtle RMUs based on highest risk and threats scores. This system also highlighted important gaps in available information that is crucial for accurate conservation assessments. Overall, this priority-setting framework can provide guidance for research and conservation priorities at multiple relevant scales, and should serve as a model for conservation status assessments and priority-setting for widespread, long-lived taxa.

Regional management units for marine turtles: a novel framework for prioritizing conservation and research across multiple scales.

BACKGROUND: Resolving threats to widely distributed marine megafauna requires definition of the geographic distributions of both the threats as well as the population unit(s) of interest. In turn, because individual threats can operate on varying spatial scales, their impacts can affect different segments of a population of the same species. Therefore, integration of multiple tools and techniques--including site-based monitoring, genetic analyses, mark-recapture studies and telemetry--can facilitate robust definitions of population segments at multiple biological and spatial scales to address different management and research challenges. METHODOLOGY/PRINCIPAL FINDINGS: To address these issues for marine turtles, we collated all available studies on marine turtle biogeography, including nesting sites, population abundances and trends, population genetics, and satellite telemetry. We georeferenced this information to generate separate layers for nesting sites, genetic stocks, and core distributions of population segments of all marine turtle species. We then spatially integrated this information from fine- to coarse-spatial scales to develop nested envelope models, or Regional Management Units (RMUs), for marine turtles globally. CONCLUSIONS/SIGNIFICANCE: The RMU framework is a solution to the challenge of how to organize marine turtles into units of protection above the level of nesting populations, but below the level of species, within regional entities that might be on independent evolutionary trajectories. Among many potential applications, RMUs provide a framework for identifying data gaps, assessing high diversity areas for multiple species and genetic stocks, and evaluating conservation status of marine turtles. Furthermore, RMUs allow for identification of geographic barriers to gene flow, and can provide valuable guidance to marine spatial planning initiatives that integrate spatial distributions of protected species and human activities. In addition, the RMU framework--including maps and supporting metadata--will be an iterative, user-driven tool made publicly available in an online application for comments, improvements, download and analysis.