Divergent sensory and immune gene evolution in sea turtles with contrasting demographic and life histories

Sea turtles represent an ancient lineage of marine vertebrates that evolved from terrestrial ancestors over 100 MYA, yet the genomic basis of the unique physiological and ecological traits enabling these species to thrive in diverse marine habitats remains largely unknown. Additionally, many populations have drastically declined due to anthropogenic activities over the past two centuries, and their recovery is a high global conservation priority. We generated and analyzed high-quality reference genomes for the leatherback (Dermochelys coriacea) and green (Chelonia mydas) turtles, representing the two extant sea turtle families. These genomes are highly syntenic and homologous, but localized regions of non-collinearity were associated with higher copy numbers of immune, zinc-finger, and olfactory receptor (OR) genes in green turtles, with ORs related to waterborne odorants greatly expanded in green turtles. Our findings suggest that divergent evolution of these key gene families may underlie immunological and sensory adaptations assisting navigation, occupancy of neritic versus pelagic environments, and diet specialization. Reduced collinearity was especially prevalent in microchromosomes, with greater gene content, heterozygosity, and genetic distances between species, supporting their critical role in vertebrate evolutionary adaptation. Finally, diversity and demographic histories starkly contrasted between species, indicating that leatherback turtles have had a low yet stable effective population size, exhibit extremely low diversity compared to other reptiles, and harbor a higher genetic load compared to green turtles, reinforcing concern over their persistence under future climate scenarios. These genomes provide invaluable resources for advancing our understanding of evolution and conservation best practices in an imperiled vertebrate lineage. Statement of significance Sea turtle populations have undergone recent global declines. We analyzed de novo assembled genomes for both extant sea turtle families through the Vertebrate Genomes Project to inform their conservation and evolutionary biology. These highly conserved genomes were differentiated by localized gene-rich regions of divergence, particularly within microchromosomes, suggesting that these genomic elements play key functional roles in the evolution of sea turtles and possibly other vertebrates. We further demonstrate that dissimilar evolutionary histories impact standing genomic diversity and genetic load, and are critical to consider when using these metrics to assess adaptive potential and extinction risk. Our results also demonstrate how reference genome quality impacts inferences of comparative and conservation genomics analyses that need to be considered in their application.

[1]  Christopher C. Kyriazis,et al.  The critically endangered vaquita is not doomed to extinction by inbreeding depression , 2022, Science.

[2]  Anna E. Savage,et al.  Adaptive evolution of major histocompatibility complex class I immune genes and disease associations in coastal juvenile sea turtles , 2022, Royal Society Open Science.

[3]  Robert M. Waterhouse,et al.  The era of reference genomes in conservation genomics. , 2022, Trends in ecology & evolution.

[4]  A. Weeks,et al.  Conservation genetics as a management tool: The five best-supported paradigms to assist the management of threatened species , 2021, Proceedings of the National Academy of Sciences.

[5]  J. DeWoody,et al.  Life-history traits and habitat availability shape genomic diversity in birds: implications for conservation , 2021, Proceedings of the Royal Society B.

[6]  Alonzo Alfaro-Núñez,et al.  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 , 2021, Heredity.

[7]  Huanming Yang,et al.  Comparative genomics provides insights into the aquatic adaptations of mammals , 2021, Proceedings of the National Academy of Sciences.

[8]  Fergal J. Martin,et al.  Population genomics of the critically endangered kākāpō , 2021, Cell genomics.

[9]  K. Bjorndal,et al.  Divergence and hybridization in sea turtles: Inferences from genome data show evidence of ancient gene flow between species , 2021, Molecular ecology.

[10]  J. Graves,et al.  Microchromosomes are building blocks of bird, reptile, and mammal chromosomes , 2021, Proceedings of the National Academy of Sciences.

[11]  P. Hedrick,et al.  The crucial role of genome-wide genetic variation in conservation , 2021, Proceedings of the National Academy of Sciences.

[12]  J. DeWoody,et al.  The long‐standing significance of genetic diversity in conservation , 2021, Molecular ecology.

[13]  Q. Cronk,et al.  Whole‐genome sequencing and genome regions of special interest: Lessons from major histocompatibility complex, sex determination, and plant self‐incompatibility , 2021, Molecular ecology.

[14]  U. Ramakrishnan,et al.  Genomic evidence for inbreeding depression and purging of deleterious genetic variation in Indian tigers , 2021, Proceedings of the National Academy of Sciences.

[15]  Yoshiaki Tanaka,et al.  Morphological features of the nasal cavities of hawksbill, olive ridley, and black sea turtles: Comparative studies with green, loggerhead and leatherback sea turtles , 2021, PloS one.

[16]  A. Prasad,et al.  Evaluating the role of reference-genome phylogenetic distance on evolutionary inference , 2021, bioRxiv.

[17]  H. Shaffer,et al.  A global phylogeny of turtles reveals a burst of climate-associated diversification on continental margins , 2021, Proceedings of the National Academy of Sciences.

[18]  A. Mooers,et al.  Do We Need to Identify Adaptive Genetic Variation When Prioritizing Populations for Conservation? , 2021, Conservation Genetics.

[19]  D. Kondoh,et al.  The nasal cavity in sea turtles: adaptation to olfaction and seawater flow , 2021, Cell and tissue research.

[20]  Joanna L. Kelley,et al.  The changing face of genome assemblies: Guidance on achieving high‐quality reference genomes , 2020, Molecular ecology resources.

[21]  Anushya Muruganujan,et al.  PANTHER version 16: a revised family classification, tree-based classification tool, enhancer regions and extensive API , 2020, Nucleic Acids Res..

[22]  Silvio C. E. Tosatto,et al.  The InterPro protein families and domains database: 20 years on , 2020, Nucleic Acids Res..

[23]  L. M. Zimmerman,et al.  The reptilian perspective on vertebrate immunity: 10 years of progress , 2020, Journal of Experimental Biology.

[24]  Voichita D. Marinescu,et al.  Progressive Cactus is a multiple-genome aligner for the thousand-genome era , 2020, Nature.

[25]  M. Fabbri,et al.  Olfactory receptor gene evolution is unusually rapid across Tetrapoda and outpaces chemosensory phenotypic change. , 2020, Current zoology.

[26]  Ryan K. Schott,et al.  The tuatara genome reveals ancient features of amniote evolution , 2020, Nature.

[27]  Y. Yamaguchi,et al.  Behavioral effects of scents from male mature Rathke glands on juvenile green sea turtles (Chelonia mydas) , 2020, The Journal of veterinary medical science.

[28]  Graham M. Hughes,et al.  Six reference-quality genomes reveal evolution of bat adaptations , 2020, Nature.

[29]  M. Waters,et al.  Rapid range shifts and megafaunal extinctions associated with late Pleistocene climate change , 2020, Nature Communications.

[30]  Y. Yamaguchi,et al.  Computed tomographic analysis of internal structures within the nasal cavities of green, loggerhead and leatherback sea turtles , 2020, Anatomical record.

[31]  A. Phillippy,et al.  Reference genome and demographic history of the most endangered marine mammal, the vaquita , 2020, bioRxiv.

[32]  Sergey Koren,et al.  Towards complete and error-free genome assemblies of all vertebrate species , 2020, Nature.

[33]  R. W. Baird,et al.  Runs of homozygosity in killer whale genomes provide a global record of demographic histories , 2020, bioRxiv.

[34]  W. Funk,et al.  Population genomics for wildlife conservation and management , 2020, Molecular ecology.

[35]  Jeffrey C. Mangel,et al.  Enhanced, coordinated conservation efforts required to avoid extinction of critically endangered Eastern Pacific leatherback turtles , 2020, Scientific Reports.

[36]  J. Wyneken,et al.  Feeding behavior and visual field differences in loggerhead and leatherback sea turtles may explain differences in longline fisheries interactions , 2020 .

[37]  Alexander Suh,et al.  Identifying the causes and consequences of assembly gaps using a multiplatform genome assembly of a bird‐of‐paradise , 2019, bioRxiv.

[38]  Andrew G. Clark,et al.  RepeatModeler2: automated genomic discovery of transposable element families , 2019, bioRxiv.

[39]  S. Kelly,et al.  OrthoFinder: phylogenetic orthology inference for comparative genomics , 2019, Genome Biology.

[40]  Parice A. Brandies,et al.  The Value of Reference Genomes in the Conservation of Threatened Species , 2019, Genes.

[41]  Jingqi Zhou,et al.  Convergent degeneration of olfactory receptor gene repertoires in marine mammals , 2019, BMC Genomics.

[42]  Christopher C. Kyriazis,et al.  Strongly deleterious mutations are a primary determinant of extinction risk due to inbreeding depression , 2019, bioRxiv.

[43]  Bent Petersen,et al.  Narwhal Genome Reveals Long-Term Low Genetic Diversity despite Current Large Abundance Size , 2019, iScience.

[44]  J. Weishampel,et al.  Effects of future sea level rise on coastal habitat , 2019, The Journal of Wildlife Management.

[45]  Kelly R. Stewart,et al.  A versatile Rapture (RAD‐Capture) platform for genotyping marine turtles , 2019, Molecular ecology resources.

[46]  K. Lohmann,et al.  There and back again: natal homing by magnetic navigation in sea turtles and salmon , 2019, Journal of Experimental Biology.

[47]  Nicholas J. Davison,et al.  Killer whale genomes reveal a complex history of recurrent admixture and vicariance , 2019, bioRxiv.

[48]  K. White,et al.  Giant tortoise genomes provide insights into longevity and age-related disease , 2018, Nature Ecology & Evolution.

[49]  R. Wayne,et al.  Genomic signatures of extensive inbreeding in Isle Royale wolves, a population on the threshold of extinction , 2018, Science Advances.

[50]  B. Shapiro,et al.  Conservation of biodiversity in the genomics era , 2018, Genome Biology.

[51]  J. DeWoody,et al.  Runs of homozygosity have utility in mammalian conservation and evolutionary studies , 2018, Conservation Genetics.

[52]  T. Zerjal,et al.  Cancer- and behavior-related genes are targeted by selection in the Tasmanian devil (Sarcophilus harrisii) , 2018, PloS one.

[53]  B. Godley,et al.  Mediterranean sea turtles: current knowledge and priorities for conservation and research , 2018 .

[54]  B. Godley,et al.  Mitochondrial DNA short tandem repeats unveil hidden population structuring and migration routes of an endangered marine turtle , 2018 .

[55]  Graham J. Etherington,et al.  Adaptation and conservation insights from the koala genome , 2018, Nature Genetics.

[56]  S. Johnsen,et al.  Green sea turtle (Chelonia mydas) population history indicates important demographic changes near the mid-Pleistocene transition , 2018, Marine Biology.

[57]  Christophe Klopp,et al.  D-GENIES: dot plot large genomes in an interactive, efficient and simple way , 2018, PeerJ.

[58]  M. Hofreiter,et al.  Extended and Continuous Decline in Effective Population Size Results in Low Genomic Diversity in the World’s Rarest Hyena Species, the Brown Hyena , 2018, Molecular biology and evolution.

[59]  Peter K. Joshi,et al.  Runs of homozygosity: windows into population history and trait architecture , 2018, Nature Reviews Genetics.

[60]  P. Dutton,et al.  Environmental Warming and Feminization of One of the Largest Sea Turtle Populations in the World , 2018, Current Biology.

[61]  M. Read,et al.  A review of fibropapillomatosis in Green turtles (Chelonia mydas). , 2016, Veterinary journal.

[62]  Diego Ortega-Del Vecchyo,et al.  Genomic Flatlining in the Endangered Island Fox , 2016, Current Biology.

[63]  Nathan F. Putman,et al.  Multi-Modal Homing in Sea Turtles: Modeling Dual Use of Geomagnetic and Chemical Cues in Island-Finding , 2016, Front. Behav. Neurosci..

[64]  H. Ellegren,et al.  PSMC analysis of effective population sizes in molecular ecology and its application to black‐and‐white Ficedula flycatchers , 2016, Molecular ecology.

[65]  D. Ray,et al.  Contrasting Patterns of Evolutionary Diversification in the Olfactory Repertoires of Reptile and Bird Genomes , 2016, Genome biology and evolution.

[66]  M. Diekhans,et al.  Genomic legacy of the African cheetah, Acinonyx jubatus , 2015, Genome Biology.

[67]  Wen J. Li,et al.  Reference sequence (RefSeq) database at NCBI: current status, taxonomic expansion, and functional annotation , 2015, Nucleic Acids Res..

[68]  S. Kelly,et al.  OrthoFinder: solving fundamental biases in whole genome comparisons dramatically improves orthogroup inference accuracy , 2015, Genome Biology.

[69]  Agostinho Antunes,et al.  Olfactory Receptor Subgenomes Linked with Broad Ecological Adaptations in Sauropsida. , 2015, Molecular biology and evolution.

[70]  Luis E Escobar,et al.  A global map of suitability for coastal Vibrio cholerae under current and future climate conditions. , 2015, Acta tropica.

[71]  Peter H. Sudmant,et al.  Mountain gorilla genomes reveal the impact of long-term population decline and inbreeding , 2015, Science.

[72]  G. Litman,et al.  Massive expansion and functional divergence of innate immune genes in a protostome , 2015, Scientific Reports.

[73]  L. Duret,et al.  Comparative population genomics in animals uncovers the determinants of genetic diversity , 2014, Nature.

[74]  E. Teeling,et al.  How and why should we implement genomics into conservation? , 2014, Evolutionary applications.

[75]  B. Godley,et al.  Pan-Atlantic analysis of the overlap of a highly migratory species, the leatherback turtle, with pelagic longline fisheries , 2014, Proceedings of the Royal Society B: Biological Sciences.

[76]  Melissa J. Landrum,et al.  RefSeq: an update on mammalian reference sequences , 2013, Nucleic Acids Res..

[77]  S. Pan,et al.  Genome analysis and signature discovery for diving and sensory properties of the endangered Chinese alligator , 2013, Cell Research.

[78]  K. Lohmann,et al.  Detection of coastal mud odors by loggerhead sea turtles: a possible mechanism for sensing nearby land , 2013 .

[79]  D. Pike Climate influences the global distribution of sea turtle nesting , 2013 .

[80]  Bronwen L. Aken,et al.  The draft genomes of soft–shell turtle and green sea turtle yield insights into the development and evolution of the turtle–specific body plan , 2013, Nature Genetics.

[81]  David Haussler,et al.  HAL: a hierarchical format for storing and analyzing multiple genome alignments , 2013, Bioinform..

[82]  Heng Li Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM , 2013, 1303.3997.

[83]  I. Hanski,et al.  High genetic load in an old isolated butterfly population , 2012, Proceedings of the National Academy of Sciences.

[84]  Y. Levy,et al.  Mitochondrial DNA STR analysis as a tool for studying the green sea turtle (Chelonia mydas) populations: the Mediterranean Sea case study. , 2012, Marine genomics.

[85]  Xiongfei Zhang,et al.  Limitations of the rhesus macaque draft genome assembly and annotation , 2012, BMC Genomics.

[86]  J. A. Sánchez,et al.  Heterozygosity-fitness correlations in the gilthead sea bream Sparus aurata using microsatellite loci from unknown and gene-rich genomic locations. , 2011, Journal of fish biology.

[87]  Heng Li,et al.  A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data , 2011, Bioinform..

[88]  Roderic B. Mast,et al.  Global Conservation Priorities for Marine Turtles , 2011, PloS one.

[89]  R. Ohlemüller,et al.  Rapid Range Shifts of Species Associated with High Levels of Climate Warming , 2011, Science.

[90]  R. Durbin,et al.  Inference of human population history from individual whole-genome sequences. , 2011, Nature.

[91]  M. DePristo,et al.  The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. , 2010, Genome research.

[92]  Menna E. Jones,et al.  MHC gene copy number variation in Tasmanian devils: implications for the spread of a contagious cancer , 2010, Proceedings of the Royal Society B: Biological Sciences.

[93]  David Haussler,et al.  Cactus Graphs for Genome Comparisons , 2010, RECOMB.

[94]  Hans H. Cheng,et al.  Comparison of linkage disequilibrium and haplotype diversity on macro- and microchromosomes in chicken , 2009, BMC Genetics.

[95]  K. Lohmann,et al.  Perception of airborne odors by loggerhead sea turtles , 2009, Journal of Experimental Biology.

[96]  M. Chaloupka,et al.  Rise and Fall over 26 Years of a Marine Epizootic in Hawaiian Green Sea Turtles , 2009, Journal of wildlife diseases.

[97]  Gonçalo R. Abecasis,et al.  The Sequence Alignment/Map format and SAMtools , 2009, Bioinform..

[98]  Matthew H. Godfrey,et al.  Climate change and marine turtles , 2009 .

[99]  Hanyi Zhuang,et al.  Odor coding by a mammalian receptor repertoire , 2009, Neuroscience Research.

[100]  N. Mrosovsky,et al.  Leatherback turtles: the menace of plastic. , 2009, Marine pollution bulletin.

[101]  C. Limpus,et al.  Ontogenetic changes in diet and habitat use in green sea turtle (Chelonia mydas) life history , 2008 .

[102]  Manuel A. R. Ferreira,et al.  PLINK: a tool set for whole-genome association and population-based linkage analyses. , 2007, American journal of human genetics.

[103]  R. Reina,et al.  Reassessment of the Leatherback Turtle (Dermochelys coriacea) Nesting Population at Parque Nacional Marino Las Baulas, Costa Rica: Effects of Conservation Efforts , 2007 .

[104]  A. Barragán,et al.  Conservation and Biology of the Leatherback Turtle in the Mexican Pacific , 2007 .

[105]  J. Jackson,et al.  Conservation implications of historic sea turtle nesting beach loss , 2006 .

[106]  M. Nei,et al.  Evolutionary dynamics of olfactory and other chemosensory receptor genes in vertebrates , 2006, Journal of Human Genetics.

[107]  J. DeWoody,et al.  On the estimation of genome-wide heterozygosity using molecular markers. , 2005, The Journal of heredity.

[108]  Nansheng Chen,et al.  Using RepeatMasker to Identify Repetitive Elements in Genomic Sequences , 2009, Current protocols in bioinformatics.

[109]  G. Hewitt Genetic consequences of climatic oscillations in the Quaternary. , 2004, Philosophical transactions of the Royal Society of London. Series B, Biological sciences.

[110]  Craig Moritz,et al.  Strategies to protect biological diversity and the evolutionary processes that sustain it. , 2002, Systematic biology.

[111]  K. Bjorndal,et al.  Historical Overfishing and the Recent Collapse of Coastal Ecosystems , 2001, Science.

[112]  P. Dutton,et al.  Global phylogeography of the leatherback turtle ( Dermochelys coriacea ) , 1999 .

[113]  R. Hirayama Oldest known sea turtle , 1998, Nature.

[114]  J. Davenport Temperature and the life-history strategies of sea turtles , 1997 .

[115]  A. Rodionov [Micro vs. macro: structural-functional organization of avian micro- and macrochromosomes]. , 1996, Genetika.

[116]  J. Mcvey,et al.  Olfactory-Based Orientation in Artificially Imprinted Sea Turtles , 1984, Science.

[117]  N. Mrosovsky,et al.  Body Temperature of Dermochelys coriacea: Warm Turtle from Cold Water , 1972, Science.

[118]  D. Ehrenfeld,et al.  CHEMORECEPTION IN THE MIGRATORY SEA TURTLE, CHELONIA MYDAS , 1972 .

[119]  Mosè Manni,et al.  BUSCO: Assessing Genome Assembly and Annotation Completeness. , 2019, Methods in molecular biology.

[120]  M. Salmon,et al.  Role of chemical and visual cues in food recognition by leatherback posthatchlings (Dermochelys coriacea L). , 2003, Zoology.

[121]  C. L. Yntema,et al.  Incubation temperature and sex ratio in hatchling loggerhead turtles: a preliminary report , 1979 .