Sea Turtle Population Genomic Discovery: Global and Locus-Specific Signatures of Polymorphism, Selection, and Adaptive Potential

The first high-quality chromosome-level genome of Eretmochelys imbricata using HiFi and Hi-C data

Eretmochelys imbricata, a critically endangered sea turtle inhabiting tropical oceans and protected across the world, had an unknown genome sequence until now. In this study, we used HiFi reads and Hi-C technology to assemble a high-quality, chromosome-level genome of E. imbricata. The genome size was 2,138.26 Mb, with contig N50 length of 123.49 Mb and scaffold N50 of 137.21 Mb. Approximately 97.52% of the genome sequence was anchored to 28 chromosomes. A total of 20,206 protein-coding genes were predicted. We also analyzed the evolutionary relationships, gene family expansions, and positive selection of E. imbricata. Our results revealed that E. imbricata diverged from Chelonia mydas 38 million years ago and had enriched olfactory receptors and aging-related genes. Our genome will be useful for studying E. imbricata and its conservation.

Tissue and Temperature-Specific RNA-Seq Analysis Reveals Genomic Versatility and Adaptive Potential in Wild Sea Turtle Hatchlings (Caretta caretta)

Simple Summary

Digital transcriptomics is rapidly emerging as a powerful new technology for modelling the environmental dynamics of the adaptive landscape in diverse lineages. This is particularly valuable in taxa such as turtles and tortoises (order Testudines) which contain a large fraction of endangered species at risk due to anthropogenic impacts on the environment, including pollution, overharvest, habitat degradation, and climate change. Sea turtles (family Cheloniidae) in particular invite a genomics-enabled approach to investigating their remarkable portfolio of adaptive evolution. Our de novo transcriptome assemblies and measurements of tissue- and temperature- specific global gene expression in the loggerhead sea turtle (Caretta caretta) reveal the genomic basis for potential resiliency in this endangered flagship species, and are crucial to future management and conservation strategies with attention to changing climates. We summarize the interactions among differentially expressed genes by producing network visualizations, and highlight the shared biological pathways related to development, migration, immunity, and longevity reported in the avian and reptilian literature. Our original results for loggerhead sea turtles provide a large, new comparative genomic resource for the investigation of genotype–phenotype relationships in amniotes.

Abstract

Background: Digital transcriptomics is rapidly emerging as a powerful new technology for modelling the environmental dynamics of the adaptive landscape in diverse lineages. This is particularly valuable in taxa such as turtles and tortoises (order Testudines) which contain a large fraction of endangered species at risk due to anthropogenic impacts on the environment, including pollution, overharvest, habitat degradation, and climate change. Sea turtles (family Cheloniidae) in particular invite a genomics-enabled approach to investigating their remarkable portfolio of adaptive evolution. The sex of the endangered loggerhead sea turtle (Caretta caretta) is subject to temperature-dependent sex determination (TSD), a mechanism by which exposure to temperatures during embryonic development irreversibly determines sex. Higher temperatures produce mainly female turtles and lower temperatures produce mainly male turtles. Incubation temperature can have long term effects on the immunity, migratory ability, and ultimately longevity of hatchlings. We perform RNA-seq differential expression analysis to investigate tissue- and temperature-specific gene expression within brain (n = 7) and gonadal (n = 4) tissue of male and female loggerhead hatchlings. Results: We assemble tissue- and temperature-specific transcriptomes and identify differentially expressed genes relevant to sexual development and life history traits of broad adaptive interest to turtles and other amniotic species. We summarize interactions among differentially expressed genes by producing network visualizations, and highlight shared biological pathways related to migration, immunity, and longevity reported in the avian and reptile literature. Conclusions: The measurement of tissue- and temperature-specific global gene expression of an endangered, flagship species such as the loggerhead sea turtle (Caretta caretta) reveals the genomic basis for potential resiliency and is crucial to future management and conservation strategies with attention to changing climates. Brain and gonadal tissue collected from experimentally reared loggerhead male and female hatchlings comprise an exceedingly rare dataset that permits the identification of genes enriched in functions related to sexual development, immunity, longevity, and migratory behavior and will serve as a large, new genomic resource for the investigation of genotype–phenotype relationships in amniotes.

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.

Population Viability of Sea Turtles in the Context of Global Warming

Sea turtles present a model for the potential impacts of climate change on imperiled species, with projected warming generating concern about their persistence. Various sea turtle life-history traits are affected by temperature; most strikingly, warmer egg incubation temperatures cause female-biased sex ratios and higher embryo mortality. Predictions of sea turtle resilience to climate change are often focused on how resulting male limitation or reduced offspring production may affect populations. In the present article, by reviewing research on sea turtles, we provide an overview of how temperature impacts on incubating eggs may cascade through life history to ultimately affect population viability. We explore how sex-specific patterns in survival and breeding periodicity determine the differences among offspring, adult, and operational sex ratios. We then discuss the implications of skewed sex ratios for male-limited reproduction, consider the negative correlation between sex ratio skew and genetic diversity, and examine consequences for adaptive potential. Our synthesis underscores the importance of considering the effects of climate throughout the life history of any species. Lethal effects (e.g., embryo mortality) are relatively direct impacts, but sublethal effects at immature life-history stages may not alter population growth rates until cohorts reach reproductive maturity. This leaves a lag during which some species transition through several stages subject to distinct biological circumstances and climate impacts. These perspectives will help managers conceptualize the drivers of emergent population dynamics and identify existing knowledge gaps under different scenarios of predicted environmental change.

Helping decision making for reliable and cost-effective 2b-RAD sequencing and genotyping analyses in non-model species

High-throughput sequencing has revolutionized population and conservation genetics. RAD sequencing methods, such as 2b-RAD, can be used on species lacking a reference genome. However, transferring protocols across taxa can potentially lead to poor results. We tested two different IIB enzymes (AlfI and CspCI) on two species with different genome sizes (the loggerhead turtle Caretta caretta and the sharpsnout seabream Diplodus puntazzo) to build a set of guidelines to improve 2b-RAD protocols on non-model organisms while optimising costs. Good results were obtained even with degraded samples, showing the value of 2b-RAD in studies with poor DNA quality. However, library quality was found to be a critical parameter on the number of reads and loci obtained for genotyping. Resampling analyses with different number of reads per individual showed a trade-off between number of loci and number of reads per sample. The resulting accumulation curves can be used as a tool to calculate the number of sequences per individual needed to reach a mean depth ≥20 reads to acquire good genotyping results. Finally, we demonstrated that selective-base ligation does not affect genomic differentiation between individuals, indicating that this technique can be used in species with large genome sizes to adjust the number of loci to the study scope, to reduce sequencing costs and to maintain suitable sequencing depth for a reliable genotyping without compromising the results. Here, we provide a set of guidelines to improve 2b-RAD protocols on non-model organisms with different genome sizes, helping decision-making for a reliable and cost-effective genotyping.