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

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

Mitochondrial genomes and genetic structure of the Kemp's ridley sea turtle (Lepidochelys kempii)

The Kemp's ridley (Lepidochelys kempii) is the world's most endangered sea turtle species and is primarily distributed in the Gulf of Mexico. In the United States, South Padre Island, Texas serves as a key nesting ground for the species. Genetic studies of the Kemp's ridley have been used to aid in conservation and management practices, with the mitochondrial control region as the most commonly used marker due to its perceived hypervariability and ease of sequencing. However, with the advent of next generation sequencing technology, targeting complete mitochondrial genomes is now feasible. Here, we describe a more complete mitochondrial genome for the Kemp's ridley than has been previously published in literature and demonstrate a cost-effective and efficient method for obtaining complete mitochondrial genomes from sea turtles. We compare the genetic diversity and taxonomic resolution obtained from whole mitochondrial genomes to that obtained from the mitochondrial control region alone. We compare current genetic diversity with previous records. Furthermore, we evaluate the genetic structure between the breeding stock in South Padre Island and that of deceased Kemp's ridleys recovered on the Northern coast of the Gulf of Mexico after the 2010 BP Deepwater Horizon oil spill, and of Kemp's ridleys stranded on the East Coast of the United States. Our results show that complete mitochondrial genomes provide greater resolution than the control region alone. They also show that the genetic diversity of the Kemp's ridley has remained stable, despite large population declines, and that the genetic makeup of deceased turtles stranded after the Deepwater Horizon oil spill is indistinguishable from the breeding stock in South Padre Island, Texas.

OPEN DATA BADGE

This article has earned an Open Data Badge for making publicly available the digitally-shareable data necessary to reproduce the reported results. The data is available at https://www.ncbi.nlm.nih.gov/genbank/.

Identifying genetic lineages through shape: An example in a cosmopolitan marine turtle species using geometric morphometrics

The green turtle (Chelonia mydas) is a globally distributed marine species whose evolutionary history has been molded by geological events and oceanographic and climate changes. Divergence between Atlantic and Pacific clades has been associated with the uplift of the Panama Isthmus, and inside the Pacific region, a biogeographic barrier located west of Hawaii has restricted the gene flow between Central/Eastern and Western Pacific populations. We investigated the carapace shape of C. mydas from individuals of Atlantic, Eastern Pacific, and Western Pacific genetic lineages using geometric morphometrics to evaluate congruence between external morphology and species' phylogeography. Furthermore, we assessed the variation of carapace shape according to foraging grounds. Three morphologically distinctive groups were observed which aligned with predictions based on the species' lineages, suggesting a substantial genetic influence on carapace shape. Based on the relationship between this trait and genetic lineages, we propose the existence of at least three distinct morphotypes of C. mydas. Well-defined groups in some foraging grounds (Galapagos, Costa Rica and New Zealand) may suggest that ecological or environmental conditions in these sites could also be influencing carapace shape in C. mydas. Geometric morphometrics is a suitable tool to differentiate genetic lineages in this cosmopolitan marine species. Consequently, this study opens new possibilities to explore and test ecological and evolutionary hypotheses in species with wide morphological variation and broad geographic distribution range.

Mixed stock analysis of juvenile green turtles aggregating at two foraging grounds in Fiji reveals major contribution from the American Samoa Management Unit

In this study we assessed the breeding population, or Management Unit (MU), origin of green turtles (Chelonia mydas) present at Yadua Island and Makogai Island foraging grounds in Fiji, central South Pacific. Based on analysis of mitochondrial (mt) DNA sequences from 150 immature green turtles caught during surveys carried out in 2015–2016, we identified a total of 18 haplotypes, the most common being CmP22.1 (44%) which is a primary haplotype characterizing the American Samoa breeding population. Results of a Bayesian mixed-stock analysis reveals that the two foraging grounds are used by green turtles from the American Samoa MU (72%, Credible Interval (CI): 56–87%), New Caledonia MU (17%, CI: 6–26%) and French Polynesia MU (7%, CI: 0–23%). The prominence of the contribution we found from the American Samoa MU compared to that of French Polynesia, both which have historic telemetry and tagging data showing connectivity with Fijian foraging areas, may reflect the current relative abundance of these two nesting populations and draws attention to a need to update population surveys and identify any significant nesting in Fiji that may have been overlooked.

First Record of Fibropapillomatosis in a Green Turtle Chelonia mydas from the Baja California Peninsula

Fibropapillomatosis (FP) is characterized by multiple fibroepithelial tumors in all parts of the skin and has been reported in sea turtles worldwide. Clinically infected individuals are often emaciated and anemic. In Mexico, however, there are few records of this disease. In this study of green turtles Chelonia mydas in Laguna San Ignacio in Baja California Sur (BCS), we noted one juvenile with multifocal fibropapilloma lesions on the external upper surface of its eyes and hind flippers. Light microscopy revealed hyperkeratosis, epidermal hyperplasia, dermal papillary projections, and fibroblast proliferation. Electron microscopy revealed viral particles. Biopsies of normal skin were done to determine the origin of the turtle through genetic analysis. Its mitochondrial DNA matched that of a haplotype (CMP2) from a Hawaiian green turtle population. Finding FP in a turtle captured in BCS elucidates the need for further monitoring along the west coast of Mexico. Further investigation should include testing tumors to detect and characterize any chelonid herpesviruses and explore any association with FP and other diseases that pose a health risk to other sea turtle species. Received March 26, 2016; accepted August 3, 2016.

The developmental biogeography of hawksbill sea turtles in the North Pacific

High seas oceanic ecosystems are considered important habitat for juvenile sea turtles, yet much remains cryptic about this important life‐history period. Recent progress on climate and fishery impacts in these so‐called lost years is promising, but the developmental biogeography of hawksbill sea turtles (Eretmochelys imbricata) has not been widely described in the Pacific Ocean. This knowledge gap limits the effectiveness of conservation management for this globally endangered species. We address this with 30 years of stranding observations, 20 years of bycatch records, and recent simulations of natal dispersal trajectories in the Hawaiian Archipelago. We synthesize the analyses of these data in the context of direct empirical observations, anecdotal sightings, and historical commercial harvests from the insular Pacific. We find hawksbills 0–4 years of age, measuring 8–34 cm straight carapace length, are found predominantly in the coastal pelagic waters of Hawaii. Unlike other species, we find no direct evidence of a prolonged presence in oceanic habitats, yet satellite tracks of passive drifters (simulating natal dispersal) and our small sample sizes suggest that an oceanic phase for hawksbills cannot be dismissed. Importantly, despite over 600 million hooks deployed and nearly 6000 turtle interactions, longline fisheries have never recorded a single hawksbill take. We address whether the patterns we observe are due to population size and gear selectivity. Although most sea turtle species demonstrate clear patterns of oceanic development, hawksbills in the North Pacific may by contrast occupy a variety of ecosystems including coastal pelagic waters and shallow reefs in remote atolls. This focuses attention on hazards in these ecosystems – entanglement and ingestion of marine debris – and perhaps away from longline bycatch and decadal climate regimes that affect sea turtle development in oceanic regions.

Population structure and phylogeography reveal pathways of colonization by a migratory marine reptile (Chelonia mydas) in the central and eastern Pacific

Climate, behavior, ecology, and oceanography shape patterns of biodiversity in marine faunas in the absence of obvious geographic barriers. Marine turtles are an example of highly migratory creatures with deep evolutionary lineages and complex life histories that span both terrestrial and marine environments. Previous studies have focused on the deep isolation of evolutionary lineages (>3 mya) through vicariance; however, little attention has been given to the pathways of colonization of the eastern Pacific and the processes that have shaped diversity within the most recent evolutionary time. We sequenced 770 bp of the mtDNA control region to examine the stock structure and phylogeography of 545 green turtles from eight different rookeries in the central and eastern Pacific. We found significant differentiation between the geographically separated nesting populations and identified five distinct stocks (F_ST = 0.08–0.44, P < 0.005). Central and eastern Pacific Chelonia mydas form a monophyletic group containing 3 subclades, with Hawaii more closely related to the eastern Pacific than western Pacific populations. The split between sampled central/eastern and western Pacific haplotypes was estimated at around 0.34 mya, suggesting that the Pacific region west of Hawaii has been a more formidable barrier to gene flow in C. mydas than the East Pacific Barrier. Our results suggest that the eastern Pacific was colonized from the western Pacific via the Central North Pacific and that the Revillagigedos Islands provided a stepping-stone for radiation of green turtles from the Hawaiian Archipelago to the eastern Pacific. Our results fit with a broader paradigm that has been described for marine biodiversity, where oceanic islands, such as Hawaii and Revillagigedo, rather than being peripheral evolutionary “graveyards”, serve as sources and recipients of diversity and provide a mechanism for further radiation.

Predicting connectivity of green turtles at Palmyra Atoll, central Pacific: a focus on mtDNA and dispersal modelling

Population connectivity and spatial distribution are fundamentally related to ecology, evolution and behaviour. Here, we combined powerful genetic analysis with simulations of particle dispersal in a high-resolution ocean circulation model to investigate the distribution of green turtles foraging at the remote Palmyra Atoll National Wildlife Refuge, central Pacific. We analysed mitochondrial sequences from turtles (n = 349) collected there over 5 years (2008-2012). Genetic analysis assigned natal origins almost exclusively (approx. 97%) to the West Central and South Central Pacific combined Regional Management Units. Further, our modelling results indicated that turtles could potentially drift from rookeries to Palmyra Atoll via surface currents along a near-Equatorial swathe traversing the Pacific. Comparing findings from genetics and modelling highlighted the complex impacts of ocean currents and behaviour on natal origins. Although the Palmyra feeding ground was highly differentiated genetically from others in the Indo-Pacific, there was no significant differentiation among years, sexes or stage-classes at the Refuge. Understanding the distribution of this foraging population advances knowledge of green turtles and contributes to effective conservation planning for this threatened species.

Reproductive biology and genetic diversity of the green turtle (Chelonia mydas) in Vamizi island, Mozambique

Introduction

Vamizi, an Island located in the Western Indian Ocean, is visited by a small and not fully characterized green turtle (Chelonia mydas (L.)) population. This population is threatened by natural hazards and several human activities, which are used to identify conservation priorities for marine turtles.

It was our aim to contribute to the knowledge of marine turtles that nest in Vamizi, with respect to its regional management, and to an area that may possibly be included on the UNESCO World Heritage List due to its potential Outstanding Universal Value.

Case description

Here, we evaluate the nesting parameters (incubation period, clutch size, hatching and emergence successes rates) and patterns over an 8-year (2003 – 2010) conservation program. We also present the results of genetic diversity based on the analysis of approximately an 850 pb fragment of the mitochondrial DNA control region.

Discussion and evaluation

We found that Vamizi beaches host a small number of nesting females, approximately 52 per year, but these have shown a reduction in their length. High hatching success (88.5 ± SD 17.2%, N = 649), emergence success rates (84.5 ± SD 20.4%, N = 649) were observed, and genetic diversity (N = 135), with 11 haplotypes found (7 new). It was also observed, in the later years of this study, a reduction in the incubation period, a dislocation of the nesting peak activity and an increase in the number of flooded nests and an increase of the number of nests in areas with lower human activity.

Conclusions

Some resilience and behavioral plasticity seems to occur regarding human territory occupancy and climate changes. However, regardless of the results, aspects like what seems to be the reduction of some cohorts, the number of flooded nests and the diminishing of the incubation period (East and South facing beaches), show that conservation efforts have to be improved.

Electronic supplementary material

The online version of this article (doi:10.1186/2193-1801-3-540) contains supplementary material, which is available to authorized users.

Characterization of polymorphic microsatellite markers for the green turtle (Chelonia mydas)

We describe primers and polymerase chain reaction conditions to amplify 12 microsatellite loci from the green turtle (Chelonia mydas), including one dinucleotide, four trinucleotide and seven tetranucleotide loci. The primers were tested on 78 individuals from a Pacific population nesting in the Hawaiian Islands. The primer pairs developed in this study yielded an average of 8.33 alleles per locus (range of 3-15 alleles), an average observed heterozygosity of 0.668 (range 0.309-0.910), and an average polymorphic information content of 0.647 (range 0.287-0.894).

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.

Comparative phylogeny and historical perspectives on population genetics of the Pacific hawksbill (Eretmochelys imbricata) and green turtles (Chelonia mydas), inferred from feeding populations in the Yaeyama Islands, Japan

Mitochondrial DNA sequence polymorphisms and patterns of genetic diversity represent the genealogy and relative impacts of historical, geographic, and demographic events on populations. In this study, historical patterns of population dynamics and differentiation in hawksbill (Eretmochelys imbricata) and green turtles (Chelonia mydas) in the Pacific were estimated from feeding populations in the Yaeyama Islands, Japan. Phylogenetic relationships of the haplotypes indicated that hawksbill and green turtles in the Pacific probably underwent very similar patterns and processes of population dynamics over the last million years, with population subdivision during the early Pleistocene and population expansion after the last glacial maximum. These significant contemporary historical events were suggested to have been caused by climatic and sea-level fluctuations. On the other hand, comparing our results to long-term population dynamics in the Atlantic, population subdivisions during the early Pleistocene were specific to Pacific hawksbill and green turtles. Therefore, regional differences in historical population dynamics are suggested. Despite limited sampling locations, these results are the first step in estimating the historical trends in Pacific sea turtles by using phylogenetics and population genetics.