Mitochondrial DNA STR analysis as a tool for studying the green sea turtle (Chelonia mydas) populations: the Mediterranean Sea case study

Conservation Genetics of the Critically Endangered Southern River Terrapin (Batagur affinis) in Malaysia: Genetic Diversity and Novel Subspecies Distribution Ranges

A population genetics study was carried out on the Southern River terrapin (Batagur affinis) from four places in Peninsular Malaysia: Pasir Gajah, Kemaman (KE), Terengganu; Bukit Pinang (BP), Kedah; Bota Kanan (BK), Perak; and Bukit Paloh, Kuala Berang (KB), Terengganu. The goal of this study is to identify genetic differences in two subspecies of B. affinis in Malaysia. No previous reports were available on the genetic diversity, phylogenetic relationships and matrilineal hereditary structure of these terrapin populations in Malaysia. The sequencing identified 46 single nucleotide polymorphisms that defined six mitochondrial haplotypes in the Southern River terrapins. Tajima's D test and Fu's Fs neutrality tests were performed to evaluate the signatures of recent historical demographic events. Based on the tests, the B. affinis edwardmolli was newly subspecies identified in the west coast-northern region of Kedah state. In addition, the B. affinis edwardmolli in Bukit Paloh, Kuala Berang (KB), Terengganu (Population 4), was shown to have a single maternal lineage compared to other populations. Low genetic diversity, but significant genetic differences, were detected among the studied Southern River terrapin populations.

Detection of high heteroplasmy in complete loggerhead and hawksbill sea turtles mitochondrial genomes using RNAseq

Sea turtle populations around the world face rapid decline due to the effect of anthropogenic and environmental factors. Among the affected populations are those of hawksbill turtles (Eretmochelys imbricata) and loggerhead turtles (Caretta caretta), which is why a greater effort is currently being made in their monitoring and tracing. The intragenic degree of heteroplasmic mutations, commonly associated with diseases of variable symptoms, has not been analyzed in these species. In this study, heteroplasmy in the complete mitogenome (mtDNA) of three loggerhead turtles and one hawksbill turtle was identified from data obtained by RNAseq. Individuals Cc3, Ei1, Cc1 and Cc2 presented 0.3, 1.7, 1.8 and 7.1% of heteroplasmic mutations in all their mtDNA, respectively. The protein-coding genes that presented the highest percentage of heteroplasmy were ND4 and ND5 in individual Cc2 with 16 and 38.6%, respectively. Of the tRNA genes, only tRNA^Tyr was heteroplasmic in the four individuals with 5.63% (Cc1), 25.35% (Ei1 and Cc2) and 49.3% (Cc3). In this study, we identified the critical sites of heteroplasmy in each individual and the genetic variability of their mitogenomes. The data obtained represents the baseline for future projects that evaluate the population status of these species.

Structural Variation of the Turtle Mitochondrial Control Region

The present study describes the most comprehensive comparison of turtle mtD-loop regions to date. The primary structure was compared from DNA sequences accessed from GenBank from 48 species in 13 families of extant turtles, and secondary structures of the mtD-loop region were inferred from thermal stabilities, using the program Mfold, for each superfamiliy of turtles. Both primary and secondary structures were found to be highly variable across the order. The Cryptodira showed conservation in the primary structure at conserved sequence blocks (CSBs), but the Pleurodira displayed limited conservation of primary structural characters, other than the coreTAS, a binding site for the helicase TWINKLE, which was highly conserved in the Central and Right Domains across the order. No secondary structure was associated with a TAS, but an AT-rich fold (secondary structure) near the 3' terminus of the mtD-loop region was detected in all turtle superfamilies. Mapping of character states of structural features of the mtD-loop region revealed that most character states were autapomorphies and inferred a number of homoplasies. The Left Domain of turtles, containing no highly conserved structural elements, likely does not serve a functional role; therefore, the Central Domain in turtles is likely equivalent to the Left Domain of mammals. The AT-rich secondary structural element near the 3' terminus of the mtD-loop region may be conserved across turtles because of a functional role, perhaps containing the Light Strand Promotor, or perhaps interacting with the TWINKLE-coreTAS complex in the Central and Right Domains to regulate mtDNA replication and transcription.

Seascape Genetics and the Spatial Ecology of Juvenile Green Turtles

Understanding how ocean currents impact the distribution and connectivity of marine species, provides vital information for the effective conservation management of migratory marine animals. Here, we used a combination of molecular genetics and ocean drift simulations to investigate the spatial ecology of juvenile green turtle (Chelonia mydas) developmental habitats, and assess the role of ocean currents in driving the dispersal of green turtle hatchlings. We analyzed mitochondrial (mt)DNA sequenced from 358 juvenile green turtles, and from eight developmental areas located throughout the Southwest Indian Ocean (SWIO). A mixed stock analysis (MSA) was applied to estimate the level of connectivity between developmental sites and published genetic data from 38 known genetic stocks. The MSA showed that the juvenile turtles at all sites originated almost exclusively from the three known SWIO stocks, with a clear shift in stock contributions between sites in the South and Central Areas. The results from the genetic analysis could largely be explained by regional current patterns, as shown by the results of passive numerical drift simulations linking breeding sites to developmental areas utilized by juvenile green turtles. Integrating genetic and oceanographic data helps researchers to better understand how marine species interact with ocean currents at different stages of their lifecycle, and provides the scientific basis for effective conservation management.

Population-specific signatures of intra-individual mitochondrial DNA heteroplasmy and their potential evolutionary advantages

Heteroplasmy is the existence of more than one mitochondrial DNA (mtDNA) variant within a cell. The evolutionary mechanisms of heteroplasmy are not fully understood, despite being a very common phenomenon. Here we combined heteroplasmy measurements using high throughput sequencing on green turtles (Chelonia mydas) with simulations to understand how heteroplasmy modulates population diversity across generations and under different demographic scenarios. We found heteroplasmy to be widespread in all individuals analysed, with consistent signal in individuals across time and tissue. Significant shifts in haplotype composition were found from mother to offspring, signalling the effect of the cellular bottleneck during oogenesis as included in the model. Our model of mtDNA inheritance indicated that heteroplasmy favoured the increase of population diversity through time and buffered against population bottlenecks, thus indicating the importance of this phenomenon in species with reduced population sizes and frequent population bottlenecks like marine turtles. Individuals with recent haplotypes showed higher levels of heteroplasmy than the individuals with ancient haplotypes, suggesting a potential advantage of maintaining established copies when new mutations arise. We recommend using heteroplasmy through high throughput sequencing in marine turtles, as well as other wildlife populations, for diversity assessment, population genetics, and mixed stock analysis.

Population recovery changes population composition at a major southern Caribbean juvenile developmental habitat for the green turtle, Chelonia mydas

Understanding the population composition and dynamics of migratory megafauna at key developmental habitats is critical for conservation and management. The present study investigated whether differential recovery of Caribbean green turtle (Chelonia mydas) rookeries influenced population composition at a major juvenile feeding ground in the southern Caribbean (Lac Bay, Bonaire, Caribbean Netherlands) using genetic and demographic analyses. Genetic divergence indicated a strong temporal shift in population composition between 2006-2007 and 2015-2016 (ϕ_ST = 0.101, P < 0.001). Juvenile recruitment (<75.0 cm straight carapace length; SCL) from the north-western Caribbean increased from 12% to 38% while recruitment from the eastern Caribbean region decreased from 46% to 20% between 2006-2007 and 2015-2016. Furthermore, the product of the population growth rate and adult female abundance was a significant predictor for population composition in 2015-2016. Our results may reflect early warning signals of declining reproductive output at eastern Caribbean rookeries, potential displacement effects of smaller rookeries by larger rookeries, and advocate for genetic monitoring as a useful method for monitoring trends in juvenile megafauna. Furthermore, these findings underline the need for adequate conservation of juvenile developmental habitats and a deeper understanding of the interactions between megafaunal population dynamics in different habitats.

How many came home? Evaluating ex situ conservation of green turtles in the Cayman Islands

Ex situ management is an important conservation tool that allows the preservation of biological diversity outside natural habitats while supporting survival in the wild. Captive breeding followed by re-introduction is a possible approach for endangered species conservation and preservation of genetic variability. The Cayman Turtle Centre Ltd was established in 1968 to market green turtle (Chelonia mydas) meat and other products and replenish wild populations, thought to be locally extirpated, through captive breeding. We evaluated the effects of this re-introduction programmme using molecular markers (13 microsatellites, 800-bp D-loop and simple tandem repeat mitochondrial DNA sequences) from captive breeders (N = 257) and wild nesting females (N = 57) (sampling period: 2013-2015). We divided the captive breeders into three groups: founders (from the original stock), and then two subdivisions of F₁ individuals corresponding to two different management strategies, cohort 1995 ("C1995") and multicohort F₁ ("MCF1"). Loss of genetic variability and increased relatedness was observed in the captive stock over time. We found no significant differences in diversity among captive and wild groups, and similar or higher levels of haplotype variability when compared to other natural populations. Using parentage and sibship assignment, we determined that 90% of the wild individuals were related to the captive stock. Our results suggest a strong impact of the re-introduction programmme on the present recovery of the wild green turtle population nesting in the Cayman Islands. Moreover, genetic relatedness analyses of captive populations are necessary to improve future management actions to maintain genetic diversity in the long term and avoid inbreeding depression.

Characterization of a subtropical hawksbill sea turtle (Eretmocheyles imbricata) assemblage utilizing shallow water natural and artificial habitats in the Florida Keys

In order to provide information to better inform management decisions and direct further research, vessel-based visual transects, snorkel transects, and in-water capture techniques were used to characterize hawksbill sea turtles in the shallow marine habitats of a Marine Protected Area (MPA), the Key West National Wildlife Refuge in the Florida Keys. Hawksbills were found in hardbottom and seagrass dominated habitats throughout the Refuge, and on man-made rubble structures in the Northwest Channel near Cottrell Key. Hawksbills captured (N = 82) were exclusively juveniles and subadults with a straight standard carapace length (SSCL) ranging from 21.4 to 69.0cm with a mean of 44.1 cm (SD = 10.8). Somatic growth rates were calculated from 15 recaptured turtles with periods at large ranging from 51 to 1188 days. Mean SSCL growth rate was 7.7 cm/year (SD = 4.6). Juvenile hawksbills (<50 cm SSCL) showed a significantly higher growth rate (9.2 cm/year, SD = 4.5, N = 11) than subadult hawksbills (50–70 cm SSCL, 3.6 cm/year, SD = 0.9, N = 4). Analysis of 740 base pair mitochondrial control region sequences from 50 sampled turtles yielded 12 haplotypes. Haplotype frequencies were significantly different compared to four other Caribbean juvenile foraging aggregations, including one off the Atlantic coast of Florida. Many-to-one mixed stock analysis indicated Mexico as the primary source of juveniles in the region and also suggested that the Refuge may serve as important developmental habitat for the Cuban nesting aggregation. Serum testosterone radioimmunoassay results from 33 individuals indicated a female biased sex ratio of 3.3 females: 1 male for hawksbills in the Refuge. This assemblage of hawksbills is near the northern limit of the species range, and is one of only two such assemblages described in the waters of the continental United States. Since this assemblage resides in an MPA with intensive human use, basic information on the assemblage is vital to resource managers charged with conservation and species protection in the MPA.

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.