Incorporating distance metrics and temporal trends to refine mixed stock analysis

The distribution of marine organisms is shaped by geographic distance and oceanographic features like currents. Among migratory species, individuals from multiple populations may share feeding habitats seasonally or across life stages. Here, we introduce a modification for many-to-many mixed stock models to include distance between breeding and foraging sites as an ecological covariate and evaluate how the composition of green turtle, Chelonia mydas, juvenile mixed stock aggregations changed in response to population growth over time. Our modified many-to-many model is more informative and generally tightens credible intervals over models that do not incorporate distance. Moreover, we identified a decrease in genetic diversity in a Florida nesting site and two juvenile aggregations. Mixed stock aggregations in central Florida have changed from multiple sources to fewer dominant source populations over the past ~ 20 years. We demonstrate that shifts in contributions from source populations to mixed stock aggregations are likely associated with nesting population growth. Furthermore, our results highlight the importance of long-term monitoring and the need for periodical reassessment of reproductive populations and juvenile aggregations. Understanding how mixed stock aggregations change over time and how different life stages are connected is fundamental for the development of successful conservation plans for imperiled species.

Adaptive evolution of major histocompatibility complex class I immune genes and disease associations in coastal juvenile sea turtles

Characterizing polymorphism at the major histocompatibility complex (MHC) genes is key to understanding the vertebrate immune response to disease. Despite being globally afflicted by the infectious tumour disease fibropapillomatosis (FP), immunogenetic variation in sea turtles is minimally explored. We sequenced the α 1 peptide-binding region of MHC class I genes (162 bp) from 268 juvenile green (Chelonia mydas) and 88 loggerhead (Caretta caretta) sea turtles in Florida, USA. We recovered extensive variation (116 alleles) and trans-species polymorphism. Supertyping analysis uncovered three functional MHC supertypes corresponding to the three well-supported clades in the phylogeny. We found significant evidence of positive selection at seven amino acid sites in the class I exon. Random forest modelling and risk ratio analysis of Ch. mydas alleles uncovered one allele weakly associated with smooth FP tumour texture, which may be associated with disease outcome. Our study represents the first characterization of MHC class I diversity in Ch. mydas and the largest sample of sea turtles used to date in any study of adaptive genetic variation, revealing tremendous genetic variation and high adaptive potential to viral pathogen threats. The novel associations we identified between MHC diversity and FP outcomes in sea turtles further highlight the importance of evaluating genetic predictors of disease, including MHC and other functional markers.

Gene expression changes with tumor disease and leech parasitism in the juvenile green sea turtle skin transcriptome

Emerging infectious diseases are a major threat to biodiversity in the 21st century. Fibropapillomatosis (FP) is an epithelial tumor disease that affects immature and adult marine turtles worldwide, particularly green turtles (Chelonia mydas). We know little about the host factors contributing to FP susceptibility, in part because transcriptomic studies that compare transcript expression in turtles with and without FP are lacking. Here, we performed RNA-Seq on healthy skin tissue from immature C. mydas in the Indian River Lagoon, Florida, USA, comparing turtles (1) with and without FP and (2) with and without leech parasites, a putative vector of FP. We assembled a de novo C. mydas skin transcriptome to identify transcripts with significant differential expression (DE) across FP and leech categories. Significant DE transcripts were found across FP and leech comparisons, including 10 of the same transcripts with DE across both comparisons. Leech-positive individuals significantly upregulated different immune and viral interaction transcripts than did leech-negative individuals, including viral interaction transcripts associated with herpesvirus interactions. This finding strengthens the role of marine leeches as mechanical vectors of Chelonid herpesvirus 5 (ChHV5) which has been implicated as a causative agent of FP. FP-positive turtles upregulated several tumor progression and suppression transcripts relative to FP-negative turtles, which had no significant DE tumor progression transcripts. FP-positive turtles also upregulated significantly more protein interaction transcripts than FP-negative turtles. DE transcripts across leech comparisons showed no functional enrichment, whereas DE transcripts across FP comparisons showed some GO terms were enriched in FP-positive and FP negative turtles. Notably, only FP-negative turtles were enriched for GO terms involved in acquired and inflammatory immune gene regulation. Overall, our DE transcripts included several candidate genes that may play important roles in C. mydas resistance to or recovery from FP, highlighting that transcriptomics provides a promising venue to understand this impactful disease. Continued investigation of C. mydas responses to FP and leech affliction is imperative for species persistence and the conservation of marine ecosystems worldwide due to the essential role of sea turtles in ecosystem function and stability.