Key issues in assessing threats to sea turtles: knowledge gaps and future directions

Sea turtles are an iconic group of marine megafauna that have been exposed to multiple anthropogenic threats across their different life stages, especially in the past decades. This has resulted in population declines, and consequently many sea turtle populations are now classified as threatened or endangered globally. Although some populations of sea turtles worldwide are showing early signs of recovery, many still face fundamental threats. This is problematic since sea turtles have important ecological roles. To encourage informed conservation planning and direct future research, we surveyed experts to identify the key contemporary threats (climate change, direct take, fisheries, pollution, disease, predation, and coastal and marine development) faced by sea turtles. Using the survey results and current literature, we also outline knowledge gaps in our understanding of the impact of these threats and how targeted future research, often involving emerging technologies, could close those gaps

Data Development and Analysis Pathways for Marine Mammals and Turtles: Creating a User Interface

A major obstacle in genetic research is developing streamlined methods for analyzing large amounts of data. The statistical computer programming language R provides users with the ability to develop packages containing specific functions in order to create more accessible data analysis pipelines. However, writing code in R can still be intimidating to those with little to no coding experience. Fortunately, the R package shiny provides a framework for developing web applications based on R functions. Using shiny, we developed a user-friendly web application containing functions of the R package strataG. The strataG package contains several functions for summarizing genetic data and analyzing population structure. Researchers who do not know how to use R or may not feel comfortable with using strataG on the R command line can now use a browser-based graphical user interface to load their raw data, select and run analyses, and save the outputs generated. User-friendly data analysis tools such as ours will assist researchers as the field of biology continues to increasingly demand that scientists have advanced computing skills. The user interface of this web application is presented here for “Data Development and Analysis Pathways for Marine Mammals and Turtles.” The server will be presented by Warren Asfazadour, et al. separately in the same series

A Scientist's Guide to Achieving Broader Impacts through K-12 STEM Collaboration.

The National Science Foundation and other funding agencies are increasingly requiring broader impacts in grant applications to encourage US scientists to contribute to science education and society. Concurrently, national science education standards are using more inquiry-based learning (IBL) to increase students' capacity for abstract, conceptual thinking applicable to real-world problems. Scientists are particularly well suited to engage in broader impacts via science inquiry outreach, because scientific research is inherently an inquiry-based process. We provide a practical guide to help scientists overcome obstacles that inhibit their engagement in K-12 IBL outreach and to attain the accrued benefits. Strategies to overcome these challenges include scaling outreach projects to the time available, building collaborations in which scientists' research overlaps with curriculum, employing backward planning to target specific learning objectives, encouraging scientists to share their passion, as well as their expertise with students, and transforming institutional incentives to support scientists engaging in educational outreach

Environment and phenology shape local adaptation in thermal performance

Populations within species often exhibit variation in traits that reflect local adaptation and further shape existing adaptive potential for species to respond to climate change. However, our mechanistic understanding of how the environment shapes trait variation remains poor. Here, we used common garden experiments to quantify thermal performance in eight populations of the marine snail Urosalpinx cinerea across thermal gradients on the Atlantic and the Pacific coasts of North America. We then evaluated the relationship between thermal performance and environmental metrics derived from time-series data. Our results reveal a novel pattern of ‘mixed’ trait performance adaptation, where thermal optima were positively correlated with spawning temperature (cogradient variation), while maximum trait performance was negatively correlated with season length (countergradient variation). This counterintuitive pattern probably arises because of phenological shifts in the spawning season, whereby ‘cold’ populations delay spawning until later in the year when temperatures are warmer compared to ‘warm’ populations that spawn earlier in the year when temperatures are cooler. Our results show that variation in thermal performance can be shaped by multiple facets of the environment and are linked to organismal phenology and natural history. Understanding the impacts of climate change on organisms, therefore, requires the knowledge of how climate change will alter different aspects of the thermal environment

Diminished warming tolerance and plasticity in low-latitude populations of a marine gastropod

Models of species response to climate change often assume that physiological traits are invariant across populations. Neglecting potential intraspecific variation may overlook the possibility that some populations are more resilient or susceptible than others, creating inaccurate predictions of climate impacts. In addition, phenotypic plasticity can contribute to trait variation and may mediate sensitivity to climate. Quantifying such forms of intraspecific variation can improve our understanding of how climate can affect ecologically important species, such as invasive predators. Here, we quantified thermal performance (tolerance, acclimation capacity, developmental traits) across seven populations of the predatory marine snail (Urosalpinx cinerea) from native Atlantic and non-native Pacific coast populations in the USA. Using common garden experiments, we assessed the effects of source population and developmental acclimation on thermal tolerance and developmental traits of F1 snails. We then estimated climate sensitivity by calculating warming tolerance (thermal tolerance − habitat temperature), using field environmental data. We report that low-latitude populations had greater thermal tolerance than their high latitude counterparts. However, these same low-latitude populations exhibited decreased thermal tolerance when exposed to environmentally realistic higher acclimation temperatures. Low-latitude native populations had the greatest climate sensitivity (habitat temperatures near thermal limits). In contrast, invasive Pacific snails had the lowest climate sensitivity, suggesting that these populations are likely to persist and drive negative impacts on native biodiversity. Developmental rate significantly increased in embryos sourced from populations with greater habitat temperature but had variable effects on clutch size and hatching success. Thus, warming can produce widely divergent responses within the same species, resulting in enhanced impacts in the non-native range and extirpation in the native range. Broadly, our results highlight how intraspecific variation can alter management decisions, as this may clarify whether management efforts should be focused on many or only a few populations

Sublethal salinity stress contributes to habitat limitation in an endangered estuarine fish.

As global change alters multiple environmental conditions, predicting species' responses can be challenging without understanding how each environmental factor influences organismal performance. Approaches quantifying mechanistic relationships can greatly complement correlative field data, strengthening our abilities to forecast global change impacts. Substantial salinity increases are projected in the San Francisco Estuary, California, due to anthropogenic water diversion and climatic changes, where the critically endangered delta smelt (Hypomesus transpacificus) largely occurs in a low-salinity zone (LSZ), despite their ability to tolerate a much broader salinity range. In this study, we combined molecular and organismal measures to quantify the physiological mechanisms and sublethal responses involved in coping with salinity changes. Delta smelt utilize a suite of conserved molecular mechanisms to rapidly adjust their osmoregulatory physiology in response to salinity changes in estuarine environments. However, these responses can be energetically expensive, and delta smelt body condition was reduced at high salinities. Thus, acclimating to salinities outside the LSZ could impose energetic costs that constrain delta smelt's ability to exploit these habitats. By integrating data across biological levels, we provide key insight into the mechanistic relationships contributing to phenotypic plasticity and distribution limitations and advance the understanding of the molecular osmoregulatory responses in nonmodel estuarine fishes

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. 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 noncollinearity 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 with other reptiles, and harbor a higher genetic load compared with 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

Assessments at multiple levels of biological organization allow for an integrative determination of physiological tolerances to turbidity in an endangered fish species.

Turbidity can influence trophic levels by altering species composition and can potentially affect fish feeding strategies and predator-prey interactions. The estuarine turbidity maximum, described as an area of increased suspended particles, phytoplankton and zooplankton, generally represents a zone with higher turbidity and enhanced food sources important for successful feeding and growth in many fish species. The delta smelt (Hypomesus transpacificus) is an endangered, pelagic fish species endemic to the San Francisco Estuary and Sacramento-San Joaquin River Delta, USA, where it is associated with turbid waters. Turbidity is known to play an important role for the completion of the species' life cycle; however, turbidity ranges in the Delta are broad, and specific requirements for this fish species are still unknown. To evaluate turbidity requirements for early life stages, late-larval delta smelt were maintained at environmentally relevant turbidity levels ranging from 5 to 250 nephelometric turbidity units (NTU) for 24 h, after which a combination of physiological endpoints (molecular biomarkers and cortisol), behavioural indices (feeding) and whole-organism measures (survival) were determined. All endpoints delivered consistent results and identified turbidities between 25 and 80 NTU as preferential. Delta smelt survival rates were highest between 12 and 80 NTU and feeding rates were highest between 25 and 80 NTU. Cortisol levels indicated minimal stress between 35 and 80 NTU and were elevated at low turbidities (5, 12 and 25 NTU). Expression of stress-related genes indicated significant responses for gst, hsp70 and glut2 in high turbidities (250 NTU), and principal component analysis on all measured genes revealed a clustering of 25, 35, 50 and 80 NTU separating the medium-turbidity treatments from low- and high-turbidity treatments. Taken together, these data demonstrate that turbidity levels that are either too low or too high affect delta smelt physiological performance, causing significant effects on overall stress, food intake and mortality. They also highlight the need for turbidity to be considered in habitat and water management decisions

Plasma proteomics of green turtles (\u3cem\u3eChelonia mydas\u3c/em\u3e) reveals pathway shifts and potential biomarker candidates associated with health and disease

Evaluating sea turtle health can be challenging due to an incomplete understanding of pathophysiologic responses in these species. Proteome characterization of clinical plasma samples can provide insights into disease progression and prospective biomarker targets. A TMT-10-plex-LC–MS/MS platform was used to characterize the plasma proteome of five, juvenile, green turtles (Chelonia mydas) and compare qualitative and quantitative protein changes during moribund and recovered states. The 10 plasma samples yielded a total of 670 unique proteins. Using ≥1.2-fold change in protein abundance as a benchmark for physiologic upregulation or downregulation, 233 (34.8%) were differentially regulated in at least one turtle between moribund and recovered states. Forty-six proteins (6.9%) were differentially regulated in all five turtles with two proteins (0.3%) demonstrating a statistically significant change. A principle component analysis showed protein abundance loosely clustered between moribund samples or recovered samples and for turtles that presented with trauma (n = 3) or as intestinal floaters (n = 2). Gene Ontology terms demonstrated that moribund samples were represented by a higher number of proteins associated with blood coagulation, adaptive immune responses and acute phase response, while recovered turtle samples included a relatively higher number of proteins associated with metabolic processes and response to nutrients. Abundance levels of 48 proteins (7.2%) in moribund samples significantly correlated with total protein, albumin and/or globulin levels quantified by biochemical analysis. Differentially regulated proteins identified with immunologic and physiologic functions are discussed for their possible role in the green turtle pathophysiologic response and for their potential use as diagnostic biomarkers. These findings enhance our ability to interpret sea turtle health and further progress conservation, research and rehabilitation programs for these ecologically important species

Large-Scale Patterns of Green Turtle Trophic Ecology in the Eastern Pacific Oceans

Trophic position and niche width are fundamental components of a species\u27 ecology, reflecting resource use, and influencing key demographic parameters such as somatic growth, maturation, and survival. Concepts about a species\u27 trophic niche space have important implications for local management and habitat protection, and can shed light about resilience to changing climate for species occurring over broad spatial scales. For elusive marine animals such as sea turtles, trophic niche is challenging to study, and researchers often rely on other metrics, such as isotopic niche, as a proxy. Here, stable isotope analysis (delta C-13 and delta N-15 values) was conducted on bulk skin tissue of 718 green turtles (Chelonia mydas) distributed among 16 foraging areas in the eastern Pacific from the USA to Chile, a range spanning similar to 10,000 km. Compound-specific nitrogen isotope analysis of amino acids (CSIA-AA) was applied to 21 turtles among seven sites. Isotopic niche space was determined via Bayesian ellipse area (BEA) and convex hull area (CHA) analyses of bulk isotope values, which were also used along with amino acid delta N-15 values to determine trophic position (TP). Substantial variability in bulk tissue delta C-13 and delta N-15 values was found within and among sites, and amino acid delta N-15 values confirmed this was largely due to spatial differences in baseline nitrogen isotopic compositions, but also to a lesser extent from TP differences among the green turtle foraging populations. Isotope niche space varied among sites, influenced by the diversity of prey types and relative input of terrestrial- vs. marine-derived nutrients; BEAs were the most suitable measurement of isotopic niche space due to the larger influence of outlying values with the CHA approach. Amino acid isotope-derived TP estimates that accounted for local habitat conditions (e.g., mixed seagrass/macroalgae diet) performed the best among several approaches; TP ranged from 2.3 to 3.6, which indicates an omnivorous diet for most populations. In addition to providing additional spatial resolution for delta C-13 and delta N-15 isoscapes in the eastern Pacific, especially in coastal habitats, this study further establishes CSIA-AA as an effective tool to study the trophic ecology of sea turtles across a variety of food webs and habitats