Global population genetic structure and male-mediated gene flow in the green sea turtle (Chelonia mydas): analysis of microsatellite loci

Demographic changes in Pleistocene sea turtles were driven by past sea level fluctuations affecting feeding habitat availability

Pleistocene environmental changes are generally assumed to have dramatically affected species’ demography via changes in habitat availability, but this is challenging to investigate due to our limited knowledge of how Pleistocene ecosystems changed through time. Here, we tracked changes in shallow marine habitat availability resulting from Pleistocene sea level fluctuations throughout the last glacial cycle (120–14 thousand years ago; kya) and assessed correlations with past changes in genetic diversity inferred from genome‐wide SNPs, obtained via ddRAD sequencing, in Caribbean hawksbill turtles, which feed in coral reefs commonly found in shallow tropical waters. We found sea level regression resulted in an average 75% reduction in shallow marine habitat availability during the last glacial cycle. Changes in shallow marine habitat availability correlated strongly with past changes in hawksbill turtle genetic diversity, which gradually declined to ~1/4th of present‐day levels during the Last Glacial Maximum (LGM; 26–19 kya). Shallow marine habitat availability and genetic diversity rapidly increased after the LGM, signifying a population expansion in response to warming environmental conditions. Our results suggest a positive correlation between Pleistocene environmental changes, habitat availability and species’ demography, and that demographic changes in hawksbill turtles were potentially driven by feeding habitat availability. However, we also identified challenges associated with disentangling the potential environmental drivers of past demographic changes, which highlights the need for integrative approaches. Our conclusions underline the role of habitat availability on species’ demography and biodiversity, and that the consequences of ongoing habitat loss should not be underestimated.

The population genomic structure of green turtles (Chelonia mydas) suggests a warm-water corridor for tropical marine fauna between the Atlantic and Indian oceans during the last interglacial

The occasional westward transport of warm water of the Agulhas Current, "Agulhas leakage", around southern Africa has been suggested to facilitate tropical marine connectivity between the Atlantic and Indian oceans, but the "Agulhas leakage" hypothesis does not explain the signatures of eastward gene flow observed in many tropical marine fauna. We investigated an alternative hypothesis: the establishment of a warm-water corridor during comparatively warm interglacial periods. The "warm-water corridor" hypothesis was investigated by studying the population genomic structure of Atlantic and Southwest Indian Ocean green turtles (N = 27) using 12,035 genome-wide single nucleotide polymorphisms (SNPs) obtained via ddRAD sequencing. Model-based and multivariate clustering suggested a hierarchical population structure with two main Atlantic and Southwest Indian Ocean clusters, and a Caribbean and East Atlantic sub-cluster nested within the Atlantic cluster. Coalescent-based model selection supported a model where Southwest Indian Ocean and Caribbean populations diverged from the East Atlantic population during the transition from the last interglacial period (130-115 thousand years ago; kya) to the last glacial period (115-90 kya). The onset of the last glaciation appeared to isolate Atlantic and Southwest Indian Ocean green turtles into three refugia, which subsequently came into secondary contact in the Caribbean and Southwest Indian Ocean when global temperatures increased after the Last Glacial Maximum. Our findings support the establishment of a warm-water corridor facilitating tropical marine connectivity between the Atlantic and Southwest Indian Ocean during warm interglacials.

First Evidence of Chelonid Herpesvirus 5 (ChHV5) Infection in Green Turtles (Chelonia mydas) from Sabah, Borneo

Fibropapillomatosis (FP) of sea turtles is characterised by cutaneous tumours and is associated with Chelonid herpesvirus 5 (ChHV5), an alphaherpesvirus from the family Herpesviridae. Here, we provide the first evidence of ChHV5-associated FP in endangered Green turtles (Chelonia mydas) from Sabah, which is located at the northern region of Malaysian Borneo. The aims of our study were firstly, to determine the presence of ChHV5 in both tumour exhibiting and tumour-free turtles using molecular techniques and secondly, to determine the phylogeography of ChHV5 in Sabah. We also aim to provide evidence of ChHV5 infection through histopathological examinations. A total of 115 Green turtles were sampled from Mabul Island, Sabah. We observed three Green turtles that exhibited FP tumours and were positive for ChHV5. In addition, six clinically healthy turtles (with no presence of tumours) were also positive for the virus based on Polymerase Chain Reaction of three viral genes (Capsid protein gene UL18, Glycoprotein H gene UL22, and Glycoprotein B gene UL27). The prevalence of the ChHV5 was 5.22% in asymptomatic Green turtles. Epidermal intranuclear inclusions were identified in tumour lesions upon histopathological examination. In addition, phylogenetic analyses of the UL18, UL22, UL27, and UL30 gene sequences showed a worldwide distribution of the ChHV5 strain with no clear distinction based on geographical location suggesting an interoceanic connection and movement of the sea turtles. Thus, the emergence of ChHV5 in Green turtles in the waters of Sabah could indicate a possible threat to sea turtle populations in the future and requires further monitoring of the populations along the Bornean coast.

Green, yellow or black? Genetic differentiation and adaptation signatures in a highly migratory marine turtle

Marine species may exhibit genetic structure accompanied by phenotypic differentiation related to adaptation despite their high mobility. Two shape-based morphotypes have been identified for the green turtle (Chelonia mydas) in the Pacific Ocean: the south-central/western or yellow turtle and north-central/eastern or black turtle. The genetic differentiation between these morphotypes and the adaptation of the black turtle to environmentally contrasting conditions of the eastern Pacific region has remained a mystery for decades. Here we addressed both questions using a reduced-representation genome approach (Dartseq; 9473 neutral SNPs) and identifying candidate outlier loci (67 outlier SNPs) of biological relevance between shape-based morphotypes from eight Pacific foraging grounds (n = 158). Our results support genetic divergence between morphotypes, probably arising from strong natal homing behaviour. Genes and enriched biological functions linked to thermoregulation, hypoxia, melanism, morphogenesis, osmoregulation, diet and reproduction were found to be outliers for differentiation, providing evidence for adaptation of C. mydas to the eastern Pacific region and suggesting independent evolutionary trajectories of the shape-based morphotypes. Our findings support the evolutionary distinctness of the enigmatic black turtle and contribute to the adaptive research and conservation genomics of a long-lived and highly mobile vertebrate.

Large-scale connectivity, cryptic population structure, and relatedness in Eastern Pacific Olive ridley sea turtles (Lepidochelys olivacea)

Endangered species are grouped into genetically discrete populations to direct conservation efforts. Mitochondrial control region (mtCR) haplotypes are used to elucidate deep divergences between populations, as compared to nuclear microsatellites that can detect recent structuring. When prior populations are unknown, it is useful to subject microsatellite data to clustering and/or ordination population inference. Olive ridley sea turtles (Lepidochelys olivacea) are the most abundant sea turtle, yet few studies have characterized olive ridley population structure. Recently, clustering results of olive ridleys in the Eastern Tropical Pacific Ocean suggested weak structuring (F ST = 0.02) between Mexico and Central America. We analyzed mtCR haplotypes, new microsatellite genotypes from Costa Rica, and preexisting microsatellite genotypes from olive ridleys across the Eastern Tropical Pacific, to further explore population structuring in this region. We subjected inferred populations to multiple analyses to explore the mechanisms behind their structuring. We found 10 mtCR haplotypes from 60 turtles nesting at three sites in Costa Rica, but did not detect divergence between Costa Rican sites, or between Central America and Mexico. In Costa Rica, clustering suggested one population with no structuring, but ordination suggested four cryptic clusters with moderate structuring (F ST = 0.08, p < .001). Across the Eastern Tropical Pacific, ordination suggested nine cryptic clusters with moderate structuring (F ST = 0.103, p < .001) that largely corresponded to Mexican and Central American populations. All ordination clusters displayed significant internal relatedness relative to global relatedness (p < .001) and contained numerous sibling pairs. This suggests that broadly dispersed family lineages have proliferated in Eastern Tropical Pacific olive ridleys and corroborates previous work showing basin-wide connectivity and shallow population structure in this region. The existence of broadly dispersed kin in Eastern Tropical Pacific olive ridleys has implications for management of olive ridleys in this region, and adds to our understanding of sea turtle ecology and life history, particularly in light of the natal-homing paradigm.

Genomics overrules mitochondrial DNA, siding with morphology on a controversial case of species delimitation

Species delimitation is a major quest in biology and is essential for adequate management of the organismal diversity. A challenging example comprises the fish species of red snappers in the Western Atlantic. Red snappers have been traditionally recognized as two separate species based on morphology: Lutjanus campechanus (northern red snapper) and L. purpureus (southern red snapper). Recent genetic studies using mitochondrial markers, however, failed to delineate these nominal species, leading to the current lumping of the northern and southern populations into a single species ( L. campechanus). This decision carries broad implications for conservation and management as red snappers have been commercially over-exploited across the Western Atlantic and are currently listed as vulnerable. To address this conflict, we examine genome-wide data collected throughout the range of the two species. Population genomics, phylogenetic and coalescent analyses favour the existence of two independent evolutionary lineages, a result that confirms the morphology-based delimitation scenario in agreement with conventional taxonomy. Despite finding evidence of introgression in geographically neighbouring populations in northern South America, our genomic analyses strongly support isolation and differentiation of these species, suggesting that the northern and southern red snappers should be treated as distinct taxonomic entities.

Climate change resilience of a globally important sea turtle nesting population

Few studies have looked into climate change resilience of populations of wild animals. We use a model higher vertebrate, the green sea turtle, as its life history is fundamentally affected by climatic conditions, including temperature-dependent sex determination and obligate use of beaches subject to sea level rise (SLR). We use empirical data from a globally important population in West Africa to assess resistance to climate change within a quantitative framework. We project 200 years of primary sex ratios (1900-2100) and create a digital elevation model of the nesting beach to estimate impacts of projected SLR. Primary sex ratio is currently almost balanced, with 52% of hatchlings produced being female. Under IPCC models, we predict: (a) an increase in the proportion of females by 2100 to 76%-93%, but cooler temperatures, both at the end of the nesting season and in shaded areas, will guarantee male hatchling production; (b) IPCC SLR scenarios will lead to 33.4%-43.0% loss of the current nesting area; (c) climate change will contribute to population growth through population feminization, with 32%-64% more nesting females expected by 2120; (d) as incubation temperatures approach lethal levels, however, the population will cease growing and start to decline. Taken together with other factors (degree of foraging plasticity, rookery size and trajectory, and prevailing threats), this nesting population should resist climate change until 2100, and the availability of spatial and temporal microrefugia indicates potential for resilience to predicted impacts, through the evolution of nest site selection or changes in nesting phenology. This represents the most comprehensive assessment to date of climate change resilience of a marine reptile using the most up-to-date IPCC models, appraising the impacts of temperature and SLR, integrated with additional ecological and demographic parameters. We suggest this as a framework for other populations, species and taxa.

Genomic evolution, recombination, and inter-strain diversity of chelonid alphaherpesvirus 5 from Florida and Hawaii green sea turtles with fibropapillomatosis

Chelonid alphaherpesvirus 5 (ChHV5) is a herpesvirus associated with fibropapillomatosis (FP) in sea turtles worldwide. Single-locus typing has previously shown differentiation between Atlantic and Pacific strains of this virus, with low variation within each geographic clade. However, a lack of multi-locus genomic sequence data hinders understanding of the rate and mechanisms of ChHV5 evolutionary divergence, as well as how these genomic changes may contribute to differences in disease manifestation. To assess genomic variation in ChHV5 among five Hawaii and three Florida green sea turtles, we used high-throughput short-read sequencing of long-range PCR products amplified from tumor tissue using primers designed from the single available ChHV5 reference genome from a Hawaii green sea turtle. This strategy recovered sequence data from both geographic regions for approximately 75% of the predicted ChHV5 coding sequences. The average nucleotide divergence between geographic populations was 1.5%; most of the substitutions were fixed differences between regions. Protein divergence was generally low (average 0.08%), and ranged between 0 and 5.3%. Several atypical genes originally identified and annotated in the reference genome were confirmed in ChHV5 genomes from both geographic locations. Unambiguous recombination events between geographic regions were identified, and clustering of private alleles suggests the prevalence of recombination in the evolutionary history of ChHV5. This study significantly increased the amount of sequence data available from ChHV5 strains, enabling informed selection of loci for future population genetic and natural history studies, and suggesting the (possibly latent) co-infection of individuals by well-differentiated geographic variants.

Population genetic structure between Yap and Palau for the coral Acropora hyacinthus

Information on connectivity is becoming increasingly in demand as marine protected areas are being designed as an integral part of a network to protect marine resources at the ecosystem level. Larval dispersal and population structure, however, remain very difficult to assess. Here, we tested the predictions of a detailed oceanographic connectivity model of larval dispersal and coral recruitment within Palau and between Palau and Yap, which was developed to support the review of the existing network of marine protected areas in Palau. We used high throughput microsatellite genotyping of the coral Acropora hyacinthus to characterize population genetic structure. Pairwise F' ST values between Palau and Yap (0.10), Palau and Ngulu (0.09) and Yap and Ngulu (0.09) were all significant and similar to pairwise F' ST values of sites within Palau (0.02-0.12) and within Yap (0.02-0.09) highlighting structure at island scale and indicating that recruitment may be even more localized than previously anticipated. A bottleneck test did not reveal any signs of a founder effect between Yap and Palau. Overall, the data supports the idea that recovery of A. hyacinthus in Palau did not come exclusively from a single source but most likely came from a combination of areas, including sites within Palau. In light of these results there seems to be very little connectivity around the barrier reef and management recommendation would be to increase the number or the size of MPAs within Palau.

Mitochondrial divergence between slow- and fast-aging garter snakes

Mitochondrial function has long been hypothesized to be intimately involved in aging processes--either directly through declining efficiency of mitochondrial respiration and ATP production with advancing age, or indirectly, e.g., through increased mitochondrial production of damaging free radicals with age. Yet we lack a comprehensive understanding of the evolution of mitochondrial genotypes and phenotypes across diverse animal models, particularly in species that have extremely labile physiology. Here, we measure mitochondrial genome-types and transcription in ecotypes of garter snakes (Thamnophis elegans) that are adapted to disparate habitats and have diverged in aging rates and lifespans despite residing in close proximity. Using two RNA-seq datasets, we (1) reconstruct the garter snake mitochondrial genome sequence and bioinformatically identify regulatory elements, (2) test for divergence of mitochondrial gene expression between the ecotypes and in response to heat stress, and (3) test for sequence divergence in mitochondrial protein-coding regions in these slow-aging (SA) and fast-aging (FA) naturally occurring ecotypes. At the nucleotide sequence level, we confirmed two (duplicated) mitochondrial control regions one of which contains a glucocorticoid response element (GRE). Gene expression of protein-coding genes was higher in FA snakes relative to SA snakes for most genes, but was neither affected by heat stress nor an interaction between heat stress and ecotype. SA and FA ecotypes had unique mitochondrial haplotypes with amino acid substitutions in both CYTB and ND5. The CYTB amino acid change (Isoleucine → Threonine) was highly segregated between ecotypes. This divergence of mitochondrial haplotypes between SA and FA snakes contrasts with nuclear gene-flow estimates, but correlates with previously reported divergence in mitochondrial function (mitochondrial oxygen consumption, ATP production, and reactive oxygen species consequences).

How do climate-linked sex ratios and dispersal limit range boundaries?

BACKGROUND

Geographic ranges of ectotherms such as reptiles may be determined strongly by abiotic factors owing to causal links between ambient temperature, juvenile survival and individual sex (male or female). Unfortunately, we know little of how these factors interact with dispersal among populations across a species range. We used a simulation model to examine the effects of dispersal, temperature-dependent juvenile survival and sex determining mechanism (temperature-dependent sex determination (TSD) and genotypic sex determination (GSD)) and their interactions, on range limits in populations extending across a continuous range of air temperatures. In particular, we examined the relative importance of these parameters for population persistence to recommend targets for future empirical research.

RESULTS

Dispersal influenced the range limits of species with TSD to a greater extent than in GSD species. Whereas male dispersal led to expanded species ranges across warm (female-producing) climates, female dispersal led to expanded ranges across cool (male-producing) climates. Two-sex dispersal eliminated the influence of biased sex ratios on ranges.

CONCLUSION

The results highlight the importance of the demographic parameter of sex ratio in determining population persistence and species range limits.

DNA fingerprinting in zoology: past, present, future

In 1962, Thomas Kuhn famously argued that the progress of scientific knowledge results from periodic 'paradigm shifts' during a period of crisis in which new ideas dramatically change the status quo. Although this is generally true, Alec Jeffreys' identification of hypervariable repeat motifs in the human beta-globin gene, and the subsequent development of a technology known now as 'DNA fingerprinting', also resulted in a dramatic shift in the life sciences, particularly in ecology, evolutionary biology, and forensics. The variation Jeffreys recognized has been used to identify individuals from tissue samples of not just humans, but also of many animal species. In addition, the technology has been used to determine the sex of individuals, as well as paternity/maternity and close kinship. We review a broad range of such studies involving a wide diversity of animal species. For individual researchers, Jeffreys' invention resulted in many ecologists and evolutionary biologists being given the opportunity to develop skills in molecular biology to augment their whole organism focus. Few developments in science, even among the subsequent genome discoveries of the 21st century, have the same wide-reaching significance. Even the later development of PCR-based genotyping of individuals using microsatellite repeats sequences, and their use in determining multiple paternity, is conceptually rooted in Alec Jeffreys' pioneering work.

Finding the 'lost years' in green turtles: insights from ocean circulation models and genetic analysis

Organismal movement is an essential component of ecological processes and connectivity among ecosystems. However, estimating connectivity and identifying corridors of movement are challenging in oceanic organisms such as young turtles that disperse into the open sea and remain largely unobserved during a period known as 'the lost years'. Using predictions of transport within an ocean circulation model and data from published genetic analysis, we present to our knowledge, the first basin-scale hypothesis of distribution and connectivity among major rookeries and foraging grounds (FGs) of green turtles (Chelonia mydas) during their 'lost years'. Simulations indicate that transatlantic dispersal is likely to be common and that recurrent connectivity between the southwestern Indian Ocean and the South Atlantic is possible. The predicted distribution of pelagic juvenile turtles suggests that many 'lost years hotspots' are presently unstudied and located outside protected areas. These models, therefore, provide new information on possible dispersal pathways that link nesting beaches with FGs. These pathways may be of exceptional conservation concern owing to their importance for sea turtles during a critical developmental period.

The interplay of homing and dispersal in green turtles: a focus on the southwestern atlantic

Current understanding of spatial ecology is insufficient in many threatened marine species, failing to provide a solid basis for conservation and management. To address this issue for globally endangered green turtles, we investigated their population distribution by sequencing a mitochondrial control region segment from the Rocas Atoll courtship area (n = 30 males) and four feeding grounds (FGs) in Brazil (n = 397), and compared our findings to published data (n (nesting) = 1205; n (feeding) = 1587). At Rocas Atoll, the first Atlantic courtship area sequenced to date, we found males were differentiated from local juveniles but not from nesting females. In combination with tag data, this indicates possible male philopatry. The most common haplotypes detected at the study sites were CMA-08 and CMA-05, and significant temporal variation was not revealed. Although feeding grounds were differentiated overall, intra-regional structure was less pronounced. Ascension was the primary natal source of the study FGs, with Surinam and Trindade as secondary sources. The study clarified the primary connectivity between Trindade and Brazil. Possible linkages to African populations were considered, but there was insufficient resolution to conclusively determine this connection. The distribution of FG haplotype lineages was nonrandom and indicative of regional clustering. The study investigated impacts of population size, geographic distance, ocean currents, and juvenile natal homing on connectivity, addressed calls for increased genetic sampling in the southwestern Atlantic, and provided data important for conservation of globally endangered green turtles.

Genetic structure and cryptic diversity of Onychodactylus japonicus (Amphibia, Caudata, Hynobiidae) in northeastern Honshu, Japan, as revealed by allozymic analysis

We conducted a comprehensive allozymic analysis of 393 specimens of the Japanese clawed salamander, Onychodactylus japonicus, from 33 populations of northeastern Honshu, Japan. As a result, these populations exhibited extensive geographic genetic differentiation, and four major genetic groups (N-Tohoku, S-Tohoku, Tsukuba, and SW-Honshu groups) were consistently recognized. Of these, the Tsukuba group was geographically isolated from all the others, whereas the N- and S-Tohoku groups, and the S-Tohoku and SW-Honshu groups, respectively, were nearly parapatric, without distinct geographic barriers. The magnitude of genetic distances between the four groups, except for between the N- and S-Tohoku groups, was as large as that normally found among different hynobiid species. A structure analysis detected no admixture of the N- and S-Tohoku groups, whereas few hybrids were found between the S-Tohoku and SW-Honshu around their contact zone. However, genetic exchange between these parapatric groups appeared to be infrequent, suggesting the presence of some isolation mechanisms between them. Within each group, only the S-Tohoku group exhibited an extensive level of population genetic structure that roughly distinguishes the eastern, central, and northwestern subgroups, indicating the complexity of the phylogeographic traits of this group. These results strongly suggest that populations of O. japonicus from northeastern Japan encompass several cryptic species.

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.

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

The Mediterranean population of the green sea turtle Chelonia mydas is critically endangered. Genetic analysis of this population using the ordinary haplotyping system, based on sequence analysis of a segment of the mitochondrial DNA (mtDNA) D-loop (control region), revealed very little variation. The most common haplotype, CM-A13, was observed in all but three individuals in hundreds of samples in previous studies. In search for a more informative marker we sequenced the 3' of the mitochondrial control region which contains an AT-rich microsatellite. We found a unique pattern that consists of four AT short tandem repeats (STRs) with varying copy numbers. This allowed us to construct a new haplotyping system composed of four different STR sizes for each mtDNA sequence. Our new mitochondrial STR (mtSTR) haplotyping approach revealed 33 different haplotypes within the nesting and stranded sea turtles along the Mediterranean Israeli seashore. The Israeli coast nesting females had 10 different haplotypes that can be used for monitoring and conservation purposes. The mtSTR haplotyping system can clearly assist in fingerprinting of individual turtles. Moreover, it can be used for estimating phylogenetic distances within populations. This case study shows that the mtSTR haplotyping is applicable for the study of global green sea turtle populations and could also be considered as markers of genetic variability in other species.

Phylogenetic and morphologic analyses of a coastal fish reveals a marine biogeographic break of terrestrial origin in the southern Caribbean

Background

Marine allopatric speciation involves interplay between intrinsic organismal properties and extrinsic factors. However, the relative contribution of each depends on the taxon under study and its geographic context. Utilizing sea catfishes in the Cathorops mapale species group, this study tests the hypothesis that both reproductive strategies conferring limited dispersal opportunities and an apparent geomorphologic barrier in the Southern Caribbean have promoted speciation in this group from a little studied area of the world.

Methodology/Principal Findings

Mitochondrial gene sequences were obtained from representatives of the Cathorops mapale species group across its distributional range from Colombia to Venezuela. Morphometric and meristic analyses were also done to assess morphologic variation. Along a ∼2000 km transect, two major lineages, Cathorops sp. and C. mapale, were identified by levels of genetic differentiation, phylogenetic reconstructions, and morphological analyses. The lineages are separated by ∼150 km at the Santa Marta Massif (SMM) in Colombia. The northward displacement of the SMM into the Caribbean in the early Pleistocene altered the geomorphology of the continental margin, ultimately disrupting the natural habitat of C. mapale. The estimated ∼0.86 my divergence of the lineages from a common ancestor coincides with the timing of the SMM displacement at ∼0.78 my.

Main Conclusions/Significance

Results presented here support the hypothesis that organismal properties as well as extrinsic factors lead to diversification of the Cathorops mapale group along the northern coast of South America. While a lack of pelagic larval stages and ecological specialization are forces impacting this process, the identification of the SMM as contributing to allopatric speciation in marine organisms adds to the list of recognized barriers in the Caribbean. Comparative examination of additional Southern Caribbean taxa, particularly those with varying life history traits and dispersal capabilities, will determine the extent by which the SMM has influenced marine phylogeography in the region.

Investigating potential for depensation in marine turtles: how low can you go?

Where mechanisms inherent within the biology of a species affect individual fitness at low density, demographic-scale depensation may occur, hastening further decline and leading ultimately to population extirpation and species extinction. Reduction in fertility at low population densities has been identified in marine and terrestrial species. Using data on hatch success and hatchling-emergence success as proxies for fertilization success, we conducted a global meta-analysis of data from breeding aggregations of green turtles (Chelonia mydas) and loggerhead turtles (Caretta caretta). We found that there has been no reduction in fertility in small nesting aggregations in either of these species worldwide. We considered mechanisms within the mating strategies and reproductive biology of marine turtles that may allow for novel genetic input and facilitate enhanced gene flow among rookeries. Behavioral reproductive mechanisms, such as natal philopatry and polyandry, may mitigate potential impacts of depensation and contribute to the resilience of these species.

Cytonuclear patterns of genetic diversity and the intricate evolutionary history of the inland silverside (Menidia beryllina)

DNA sequence variation at a mitochondrial and a nuclear intron locus was surveyed within and among multiple populations of the inland silverside (Menidia beryllina) from the southeastern United States and revealed discordant phylogenetic patterns but similar patterns of population genetic variation across nuclear and mitochondrial loci. Mitochondrial variation was geographically structured, with strongly supported monophyletic assemblages among Gulf of Mexico population samples and a close association of the St John's River (SJ) population with these same samples. Nuclear alleles were not strongly structured geographically, with little support for monophyly within or across basins. Conversely, population genetic parameters indicate that the bulk of genetic diversity for both genomes resides within and among Gulf of Mexico populations and that diversity within the Atlantic is largely restricted to the SJ population. The contrast in genetic variation and population phylogenies appears to be a function of historical demographic processes, most likely directed by fluctuating geomorphology of the Florida peninsula in response to North American glaciation cycles.

Species identification using the cytochrome b gene of commercial turtle shells

Turtle shells and their gelled products are familiar in some countries as foods, tonics and medicines. These shells may come from endangered and protected species, requiring the identification of the species present to enforce national and international legislation. We report on the design of five combinations of primer pairs for the identification of turtle shells and shell fragments used as ornaments, food products and medicines. The types of samples used are those encountered frequently and will typically contain highly degraded DNA. The success rate for species identification using the test described is dependent upon the choice of primer sets used and the length of the expected amplification product. Gelled products were simulated by the process of decoction for up to 12 h, after which all the turtle species could be identified from the liquid samples. This study establishes a method for the identification of commercial turtle shells and illustrates a simulated case using gelled products.

Population genetic structure of Earth's largest fish, the whale shark (Rhincodon typus)

Large pelagic vertebrates pose special conservation challenges because their movements generally exceed the boundaries of any single jurisdiction. To assess the population structure of whale sharks (Rhincodon typus), we sequenced complete mitochondrial DNA control regions from individuals collected across a global distribution. We observed 51 single site polymorphisms and 8 regions with indels comprising 44 haplotypes in 70 individuals, with high haplotype (h = 0.974 +/- 0.008) and nucleotide diversity (pi = 0.011 +/- 0.006). The control region has the largest length variation yet reported for an elasmobranch (1143-1332 bp). Phylogenetic analyses reveal no geographical clustering of lineages and the most common haplotype was distributed globally. The absence of population structure across the Indian and Pacific basins indicates that oceanic expanses and land barriers in Southeast Asia are not impediments to whale shark dispersal. We did, however, find significant haplotype frequency differences (AMOVA, Phi(ST) = 0.107, P < 0.001) principally between the Atlantic and Indo-Pacific populations. In contrast to other recent surveys of globally distributed sharks, we find much less population subdivision and no evidence for cryptic evolutionary partitions. Discovery of the mating and pupping areas of whale sharks is key to further population genetic studies. The global pattern of shared haplotypes in whale sharks provides a compelling argument for development of broad international approaches for management and conservation of Earth's largest fish.

Dispersal, recruitment and migratory behaviour in a hawksbill sea turtle aggregation

We investigated the dispersal, recruitment and migratory behaviour of the hawksbill sea turtle (Eretmochelys imbricata), among different life-history stages and demographic segments of the large hawksbill turtle aggregation at Mona Island, Puerto Rico. There were significant differences in both mitochondrial DNA (mtDNA) haplotype diversity and haplotype frequencies among the adult males, females and juveniles examined, but little evidence for temporal heterogeneity within these same groups sampled across years. Consistent with previous studies and the hypothesis of strong natal homing, there were striking mtDNA haplotype differences between nesting females on Mona Island and nesting females in other major Caribbean rookeries. Breeding males also showed strong, albeit weaker, genetic evidence of natal homing. Overall, Bayesian mixed-stock analysis suggests that Mona Island was the natal rookery for 79% (65-94%) of males in the aggregation. In contrast, the Mona Island rookery accounted for only a small subset of the new juvenile recruits to the foraging grounds or in the population of older juvenile hawksbills turtles on Mona. Instead, both new recruits and the older juvenile hawksbill turtles on Mona more likely recruited from other Caribbean rookeries, suggesting that a mechanism besides natal homing must be influencing recruitment to feeding habitats. The difference in the apparent degree of natal homing behaviour among the different life-history stages of hawksbill turtles at Mona Island underscores the complexity of the species' life-history dynamics and highlights the need for both local and regional conservation efforts.

Population genetics and phylogeography of sea turtles

The seven species of sea turtles occupy a diversity of niches, and have a history tracing back over 100 million years, yet all share basic life-history features, including exceptional navigation skills and periodic migrations from feeding to breeding habitats. Here, we review the biogeographic, behavioural, and ecological factors that shape the distribution of genetic diversity in sea turtles. Natal homing, wherein turtles return to their region of origin for mating and nesting, has been demonstrated with mtDNA sequences. These maternally inherited markers show strong population structure among nesting colonies while nuclear loci reveal a contrasting pattern of male-mediated gene flow, a phenomenon termed 'complex population structure'. Mixed-stock analyses indicate that multiple nesting colonies can contribute to feeding aggregates, such that exploitation of turtles in these habitats can reduce breeding populations across the region. The mtDNA data also demonstrate migrations across entire ocean basins, some of the longest movements of marine vertebrates. Multiple paternity occurs at reported rates of 0-100%, and can vary by as much as 9-100% within species. Hybridization in almost every combination among members of the Cheloniidae has been documented but the frequency and ultimate ramifications of hybridization are not clear. The global phylogeography of sea turtles reveals a gradient based on habitat preference and thermal regime. The cold-tolerant leatherback turtle (Dermochelys coriacea) shows no evolutionary partitions between Indo-Pacific and Atlantic populations, while the tropical green (Chelonia mydas), hawksbill (Eretmochelys imbricata), and ridleys (Lepidochelys olivacea vs. L. kempi) have ancient separations between oceans. Ridleys and loggerhead (Caretta caretta) also show more recent colonization between ocean basins, probably mediated by warm-water gyres that occasionally traverse the frigid upwelling zone in southern Africa. These rare events may be sufficient to prevent allopatric speciation under contemporary geographic and climatic conditions. Genetic studies have advanced our understanding of marine turtle biology and evolution, but significant gaps persist and provide challenges for the next generation of sea turtle geneticists.

Genetic differentiation across a latitudinal gradient in two co-occurring butterfly species: revealing population differences in a context of climate change

Genetic differentiation within a species' range is determined by natural selection, genetic drift, and gene flow. Selection and drift enhance genetic differences if populations are sufficiently isolated, while gene flow precludes differentiation and local adaptation. Over large geographical areas, these processes can create a variety of scenarios, ranging from admixture to a high degree of population differentiation. Genetic differences among populations may signal functional differences within a species' range, potentially leading to population or ecotype-specific responses to global change. We investigated differentiation within the geographical range of two butterfly species along a broad latitudinal gradient. This gradient is the primary axis of climatic variation, and many ecologists expect populations at the poleward edge of this gradient to expand under climate change. Our study species inhabit a shared ecosystem and differ in body size and resource specialization; both also find their poleward range limit on an island. We find evidence for divergence of peripheral populations from the core in both taxa, suggesting the potential for genetic distinctiveness at the leading edge of climate change. We also find differences between the species in the extent of peripheral differentiation with the smaller and more specialized species showing greater population divergence (microsatellites and mtDNA) and reduced gene flow (mtDNA). Finally, gene flow estimates in both species differed strongly between two marker types. These findings suggest caution in assuming that populations are invariant across latitude and thus will respond as a single ecotype to climatic change.

Female philopatry in a heterogeneous environment: ordinary conditions leading to extraordinary ESS sex ratios

Background

We use a simulation-based model to study the impact of female philopatry and heterogeneity of habitat quality on the evolution of primary sex ratio.

Results

We show that these conditions may lead to strongly biased ESS habitat-dependent sex ratios, under two kinds of density-dependent population regulation. ESS sex ratios are always biased towards females in good habitats, towards males in poor habitats, and are generally equilibrated considering the whole population. Noticeably, the predicted bias of sex ratio usually increases with decreasing female philopatry.

Conclusion

The selection forces responsible for these results are fully described. This study provides a new perspective on the evolutionary significance of temperature sex determination. We discuss the case of turtles by comparing our theoretical results with field observations.

Phylogeography of the green turtle, Chelonia mydas, in the Southwest Indian Ocean

Patterns of mitochondrial DNA (mtDNA) variation were used to analyse the population genetic structure of southwestern Indian Ocean green turtle (Chelonia mydas) populations. Analysis of sequence variation over 396 bp of the mtDNA control region revealed seven haplotypes among 288 individuals from 10 nesting sites in the Southwest Indian Ocean. This is the first time that Atlantic Ocean haplotypes have been recorded among any Indo-Pacific nesting populations. Previous studies indicated that the Cape of Good Hope was a major biogeographical barrier between the Atlantic and Indian Oceans because evidence for gene flow in the last 1.5 million years has yet to emerge. This study, by sampling localities adjacent to this barrier, demonstrates that recent gene flow has occurred from the Atlantic Ocean into the Indian Ocean via the Cape of Good Hope. We also found compelling genetic evidence that green turtles nesting at the rookeries of the South Mozambique Channel (SMC) and those nesting in the North Mozambique Channel (NMC) belong to separate genetic stocks. Furthermore, the SMC could be subdivided in two different genetic stocks, one in Europa and the other one in Juan de Nova. We suggest that this particular genetic pattern along the Mozambique Channel is attributable to a recent colonization from the Atlantic Ocean and is maintained by oceanic conditions in the northern and southern Mozambique Channel that influence early stages in the green turtle life cycle.

Estimating population structure under nonequilibrium conditions in a conservation context: continent-wide population genetics of the giant Amazon river turtle, Podocnemis expansa (Chelonia; Podocnemididae)

Giant Amazon river turtles, Podocnemis expansa, are indigenous to the Amazon, Orinoco, and Essequibo River basins, and are distributed across nearly the entire width of the South American continent. Although once common, their large size, high fecundity, and gregarious nesting, made P. expansa especially vulnerable to over-harvesting for eggs and meat. Populations have been severely reduced or extirpated in many areas throughout its range, and the species is now regulated under Appendix II of the Convention on International Trade in Endangered Species. Here, we analyse data from mitochondrial DNA sequence and multiple nuclear microsatellite markers with an array of complementary analytical methods. Results show that concordance from multiple data sets and analyses can provide a strong signal of population genetic structure that can be used to guide management. The general lack of phylogeographic structure but large differences in allele and haplotype frequencies among river basins is consistent with fragmented populations and female natal-river homing. Overall, the DNA data show that P. expansa populations lack a long history of genetic differentiation, but that each major tributary currently forms a semi-isolated reproductive population and should be managed accordingly.

Phylogeographic history and gene flow among giant Galápagos tortoises on southern Isabela Island

Volcanic islands represent excellent models with which to study the effect of vicariance on colonization and dispersal, particularly when the evolution of genetic diversity mirrors the sequence of geological events that led to island formation. Phylogeographic inference, however, can be particularly challenging for recent dispersal events within islands, where the antagonistic effects of land bridge formation and vicariance can affect movements of organisms with limited dispersal ability. We investigated levels of genetic divergence and recovered signatures of dispersal events for 631 Galápagos giant tortoises across the volcanoes of Sierra Negra and Cerro Azul on the island of Isabela. These volcanoes are among the most recent formations in the Galápagos (<0.7 million years), and previous studies based on genetic and morphological data could not recover a consistent pattern of lineage sorting. We integrated nested clade analysis of mitochondrial DNA control region sequences, to infer historical patterns of colonization, and a novel Bayesian multilocus genotyping method for recovering evidence of recent migration across volcanoes using eleven microsatellite loci. These genetic studies illuminate taxonomic distinctions as well as provide guidance to possible repatriation programs aimed at countering the rapid population declines of these spectacular animals.