Species and population specific gene expression in blood transcriptomes of marine turtles

Testosterone identifies hatchling sex for Mojave desert tortoises (Gopherus agassizii)

The threatened Mojave desert tortoise (Gopherus agassizii) exhibits temperature-dependent sex determination, and individuals appear externally sexually monomorphic until sexual maturity. A non-surgical sex identification method that is suitable for a single in situ encounter with hatchlings is essential for minimizing handling of wild animals. We tested (1) whether plasma testosterone quantified by enzyme-linked immunosorbent assay differentiated males from females in 0-3 month old captive hatchlings, and (2) whether an injection of follicle-stimulating hormone (FSH) differentially elevates testosterone in male hatchlings to aid in identifying sex. We validated sex by ceolioscopic (laparoscopic) surgery. We then fit the testosterone concentrations to lognormal distributions and identified the concentration below which individuals are more likely female, and above which individuals are more likely male. Using a parametric bootstrapping procedure, we estimated a 0.01-0.04% misidentification rate for naïve testosterone samples, and a 1.26-1.39% misidentification rate for challenged (post-FSH injection) testosterone samples. Quantification of plasma testosterone concentration from small volume (0.1 mL) blood samples appears to be a viable, highly accurate method to identify sex of 0-3 month old hatchlings and could be a valuable tool for conservation measures and investigation of trends and variation in sex ratios for in situ wild nests.

Population transcriptogenomics highlights impaired metabolism and small population sizes in tree frogs living in the Chernobyl Exclusion Zone

BACKGROUND

Individual functional modifications shape the ability of wildlife populations to cope with anthropogenic environmental changes. But instead of adaptive response, human-altered environments can generate a succession of deleterious functional changes leading to the extinction of the population. To study how persistent anthropogenic changes impacted local species' population status, we characterised population structure, genetic diversity and individual response of gene expression in the tree frog Hyla orientalis along a gradient of radioactive contamination around the Chernobyl nuclear power plant.

RESULTS

We detected lower effective population size in populations most exposed to ionizing radiation in the Chernobyl Exclusion Zone that is not compensated by migrations from surrounding areas. We also highlighted a decreased body condition of frogs living in the most contaminated area, a distinctive transcriptomics signature and stop-gained mutations in genes involved in energy metabolism. While the association with dose will remain correlational until further experiments, a body of evidence suggests the direct or indirect involvement of radiation exposure in these changes.

CONCLUSIONS

Despite ongoing migration and lower total dose rates absorbed than at the time of the accident, our results demonstrate that Hyla orientalis specimens living in the Chernobyl Exclusion Zone are still undergoing deleterious changes, emphasizing the long-term impacts of the nuclear disaster.

Immunity in Sea Turtles: Review of a Host-Pathogen Arms Race Millions of Years in the Running

The immune system of sea turtles is not completely understood. Sea turtles (as reptiles) bridge a unique evolutionary gap, being ectothermic vertebrates like fish and amphibians and amniotes like birds and mammals. Turtles are ectotherms; thus, their immune system is influenced by environmental conditions like temperature and season. We aim to review the turtle immune system and note what studies have investigated sea turtles and the effect of the environment on the immune response. Turtles rely heavily on the nonspecific innate response rather than the specific adaptive response. Turtles' innate immune effectors include antimicrobial peptides, complement, and nonspecific leukocytes. The antiviral defense is understudied in terms of the diversity of pathogen receptors and interferon function. Turtles also mount adaptive responses to pathogens. Lymphoid structures responsible for lymphocyte activation and maturation are either missing in reptiles or function is affected by season. Turtles are a marker of health for their marine environment, and their immune system is commonly dysregulated because of disease or contaminants. Fibropapillomatosis (FP) is a tumorous disease that afflicts sea turtles and is thought to be caused by a virus and an environmental factor. We aim, by exploring the current understanding of the immune system in turtles, to aid the investigation of environmental factors that contribute to the pathogenesis of this disease and provide options for immunotherapy.

Blood-based gene expression as non-lethal tool for inferring salinity-habitat history of European eel (Anguilla anguilla)

The European eel is a facultative catadromous species, meaning that it can skip the freshwater phase or move between marine and freshwater habitats during its continental life stage. Otolith microchemistry, used to determine the habitat use of eel or its salinity history, requires the sacrifice of animals. In this context, blood-based gene expression may represent a non-lethal alternative. In this work, we tested the ability of blood transcriptional profiling to identify the different salinity-habitat histories of European eel. Eels collected from different locations in Norway were classified through otolith microchemistry as freshwater residents (FWR), seawater residents (SWR) or inter-habitat shifters (IHS). We detected 3451 differentially expressed genes from blood by comparing FWR and SWR groups, and then used that subset of genes in a machine learning approach (i.e., random forest) to the extended FWR, SWR, and IHS group. Random forest correctly classified 100% of FWR and SWR and 83% of the IHS using a minimum of 30 genes. The implementation of this non-lethal approach may replace otolith-based microchemistry analysis for the general assessment of life-history tactics in European eels. Overall, this approach is promising for the replacement or reduction of other lethal analyses in determining certain fish traits.

Using transcriptomics to predict and visualize disease status in bighorn sheep (Ovis canadensis)

Increasing risk of pathogen spillover coupled with overall declines in wildlife population abundance in the Anthropocene make infectious disease a relevant concern for species conservation worldwide. While emerging molecular tools could improve our diagnostic capabilities and give insight into mechanisms underlying wildlife disease risk, they have rarely been applied in practice. Here, employing a previously reported gene transcription panel of common immune markers to track physiological changes, we present a detailed analysis over the course of both acute and chronic infection in one wildlife species where disease plays a critical role in conservation, bighorn sheep (Ovis canadensis). Differential gene transcription patterns distinguished between infection statuses over the course of acute infection and differential correlation (DC) analyses identified clear changes in gene co-transcription patterns over the early stages of infection, with transcription of four genes-TGFb, AHR, IL1b and MX1-continuing to increase even as transcription of other immune-associated genes waned. In a separate analysis, we considered the capacity of the same gene transcription panel to aid in differentiating between chronically infected animals and animals in other disease states outside of acute disease events (an immediate priority for wildlife management in this system). We found that this transcription panel was capable of accurately identifying chronically infected animals in the test dataset, though additional data will be required to determine how far this ability extends. Taken together, our results showcase the successful proof of concept and breadth of potential utilities that gene transcription might provide to wildlife disease management, from direct insight into mechanisms associated with differential disease response to improved diagnostic capacity in the field.

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.

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.