Visual, spectral, and microchemical quantification of crystalline anomalies in otoliths of wild and cultured delta smelt.

Developmental abnormalities in otoliths can impact growth and survival in teleost fishes. Here, we quantified the frequency and severity of developmental anomalies in otoliths of delta smelt (Hypomesus transpacificus), a critically endangered estuarine fish that is endemic to the San Francisco Estuary. Left-right asymmetry and anomalous crystalline polymorphs (i.e., vaterite) were quantified and compared between wild and cultured populations using digital image analysis. Visual estimates of vaterite were validated using X-ray diffraction, Raman spectroscopy, laser ablation ICPMS, and electron probe microanalysis. Results indicated that cultured delta smelt were 80 times more likely to contain a vateritic otolith and 18 times more likely to contain relatively large (≥ 15%) amounts of vaterite. Similarly, cultured fish exhibited 30% greater asymmetry than wild fish. These results indicate that cultured delta smelt exhibit a significantly higher frequency of vestibular abnormalities which are known to reduce fitness and survival. Such hatchery effects on otolith development could have important implications for captive culture practices and the supplementation of wild fish populations with cultured individuals.

Polygenic discrimination of migratory phenotypes in an estuarine forage fish.

Migration is a complex phenotypic trait with some species containing migratory and nonmigratory individuals. Such life history variation may be attributed in part to plasticity, epigenetics, or genetics. Although considered semianadromous, recent studies using otolith geochemistry have revealed life history variation within the critically endangered Delta Smelt. Broadly categorizable as migratory or freshwater residents, we examined Restriction site Associated DNA sequencing data to test for a relationship between genetic variation and migratory behaviors. As previously shown, we found no evidence for neutral population genetic structure within Delta Smelt; however, we found significant evidence for associations between genetic variants and life history phenotypes. Furthermore, discriminant analysis of principal components, hierarchical clustering, and machine learning resulted in accurate assignment of fish into the freshwater resident or migratory classes based on their genotypes. These results suggest the presence of adaptive genetic variants relating to life history variation within a panmictic population. Mechanisms that may lead to this observation are genotype dependent habitat choice and spatially variable selection, both of which could operate each generation and are not exclusive. Given that the population of cultured Delta Smelt are being used as a refugial population for conservation, as a supply for wild population supplementation, and currently represent the majority of all living individuals of this species, we recommend that the hatchery management strategy consider the frequencies of life history-associated alleles and how to maintain this important aspect of Delta Smelt biological variation while under captive propagation.

Experimental validation of otolith-based age and growth reconstructions across multiple life stages of a critically endangered estuarine fish.

BACKGROUND: The application of otolith-based tools to inform the management and conservation of fishes first requires taxon- and stage-specific validation. The Delta Smelt (Hypomesus transpacificus), a critically endangered estuarine fish that is endemic to the upper San Francisco Estuary (SFE), California, United States, serves as a key indicator species in the SFE; thus, understanding this species' vital rates and population dynamics is valuable for assessing the overall health of the estuary. Otolith-based tools have been developed and applied across multiple life stages of Delta Smelt to reconstruct age structure, growth, phenology, and migration. However, key methodological assumptions have yet to be validated, thus limiting confidence in otolith-derived metrics that are important for informing major water management decisions in the SFE. METHODS: Using known-age cultured Delta Smelt and multiple independent otolith analysts, we examined otolith formation, otolith-somatic proportionality, aging accuracy and precision, left-right symmetry, and the effects of image magnification for larval, juvenile, and adult Delta Smelt. RESULTS: Overall, otolith size varied linearly with fish size (from 10-60 mm), explaining 99% of the variation in fish length, despite a unique slope for larvae < 10 mm. Otolith-somatic proportionality was similar among wild and cultured specimens. Aging precision among independent analysts was 98% and aging accuracy relative to known ages was 96%, with age estimates exhibiting negligible differences among left and right otoliths. Though error generally increased with age, percent error decreased from 0-30 days-post-hatch, with precision remaining relatively high (≥ 95%) thereafter. Increased magnification (400×) further improved aging accuracy for the oldest, slowest-growing individuals. Together, these results indicate that otolith-based techniques provide reliable age and growth reconstructions for larval, juvenile, and adult Delta Smelt. Such experimental assessments across multiple developmental stages are key steps toward assessing confidence in otolith-derived metrics that are often used to assess the dynamics of wild fish populations.

Biogeochemical processes create distinct isotopic fingerprints to track floodplain rearing of juvenile salmon.

Floodplains represent critical nursery habitats for a variety of fish species due to their highly productive food webs, yet few tools exist to quantify the extent to which these habitats contribute to ecosystem-level production. Here we conducted a large-scale field experiment to characterize differences in food web composition and stable isotopes (δ¹³C, δ¹5N, δ³4S) for salmon rearing on a large floodplain and adjacent river in the Central Valley, California, USA. The study covered variable hydrologic conditions including flooding (1999, 2017), average (2016), and drought (2012-2015). In addition, we determined incorporation rates and tissue fractionation between prey and muscle from fish held in enclosed locations (experimental fields, cages) at weekly intervals. Finally, we measured δ³4S in otoliths to test if these archival biominerals could be used to reconstruct floodplain use. Floodplain-reared salmon had a different diet composition and lower δ13C and δ³4S (δ¹³C = -33.02±2.66‰, δ³4S = -3.47±2.28‰; mean±1SD) compared to fish in the adjacent river (δ¹³C = -28.37±1.84‰, δ³4S = +2.23±2.25‰). These isotopic differences between habitats persisted across years of extreme droughts and floods. Despite the different diet composition, δ¹5N values from prey items on the floodplain (δ¹5N = 7.19±1.22‰) and river (δ¹5N = 7.25±1.46‰) were similar, suggesting similar trophic levels. The food web differences in δ13C and δ³4S between habitats were also reflected in salmon muscle tissue, reaching equilibrium between 24-30 days (2014, δ¹³C = -30.74±0.73‰, δ³4S = -4.6±0.68‰; 2016, δ¹³C = -34.74 ±0.49‰, δ³4S = -5.18±0.46‰). δ³4S measured in sequential growth bands in otoliths recorded a weekly time-series of shifting diet inputs, with the outermost layers recording time spent on the floodplain (δ³4S = -5.60±0.16‰) and river (δ³4S = 3.73±0.98‰). Our results suggest that δ¹³C and δ³4S can be used to differentiate floodplain and river rearing habitats used by native fishes, such as Chinook Salmon, across different hydrologic conditions and tissues. Together these stable isotope analyses provide a toolset to quantify the role of floodplains as fish habitats.

Complex life histories discovered in a critically endangered fish.

Effective conservation of endangered species requires knowledge of the full range of life-history strategies used to maximize population resilience within a stochastic and ever-changing environment. California's endemic Delta Smelt (Hypomesus transpacificus) is rapidly approaching extinction in the San Francisco Estuary, placing it in the crossfire between human and environmental uses of limited freshwater resources. Though managed as a semi-anadromous species, recent studies have challenged this lifecycle model for Delta Smelt, suggesting the species is an estuarine resident with several localized "hot-spots" of abundance. Using laser-ablation otolith strontium isotope microchemistry, we discovered three distinct life-history phenotypes including freshwater resident (FWR), brackish-water resident (BWR), and semi-anadromous (SA) fish. We further refined life-history phenotypes using an unsupervised algorithm and hierarchical clustering and found that in the last resilient year-class, the FWR (12%) and BWR (7%) comprised a small portion of the population, while the majority of fish were SA (81%). Furthermore, the semi-anadromous fish could be clustered into at least four additional life-history phenotypes that varied by natal origin, dispersal age and adult salinity history. These diverse life-history strategies should be incorporated into future conservation and management efforts aimed at preventing the extinction of Delta Smelt in the wild.

Rocks, teeth, and tools: New insights into early Neanderthal mobility strategies in South-Eastern France from lithic reconstructions and strontium isotope analysis.

Neanderthals had complex land use patterns, adapting to diversified landscapes and climates. Over the past decade, considerable progress has been made in reconstructing the chronology, land use and subsistence patterns, and occupation types of sites in the Rhône Valley, southeast France. In this study, Neanderthal mobility at the site of Payre is investigated by combining information from lithic procurement analysis ("chaîne evolutive" and "chaîne opératoire" concepts) and strontium isotope analysis of teeth (childhood foraging area), from two units (F and G). Both units date to the transition from Marine Isotope Stage (MIS) 8 to MIS 7, and show similar environmental conditions, but represent contrasting occupation durations. Level Gb (unit G) represents a long-term year-round use, in contrast to short-term seasonal use of the cave in level Fb (unit F). For both levels, lithic material and food were generally collected from a local to semi-local region. However, in level Gb, lithic materials were mainly collected from colluviums and food collected in the valley, whereas in level Fb, lithic procurement focused primarily on alluvial deposits and food was collected from higher elevation plateaus. These procurement or exchange patterns might be related to flint availability, knapping advantages of alluvial flint or occupation duration. The site of Payre is located in a flint rich circulation corridor and the movement of groups or exchanges between groups were organized along a north-south axis on the plateaus or towards the east following the river. The ridges were widely used as they are rich in flint, whereas the Rhône Valley is not an important source of lithic raw materials. Compared to other western European Middle Palaeolithic sites, these results indicate that procurement strategies have a moderate link with occupation types and duration, and with lithic technology. The Sr isotope ratios broadly match the proposed foraging areas, with the Rhône Valley being predominantly used in unit G and the ridges and limestone plateaus in unit F. While lithic reconstructions and childhood foraging are not directly related this suggests that the three analysed Neanderthals spend their childhood in the same general area and supports the idea of mobile Neanderthals in the Rhône Valley and neighbouring higher elevation plateaus. The combination of reconstructing lithic raw material sources, provisioning strategies, and strontium isotope analyses provides new details on how Neanderthals at Payre practised land use and mobility in the Early Middle Palaeolithic.

Calibrating temperature reconstructions from fish otolith oxygen isotope analysis for California's critically endangered Delta Smelt.

RATIONALE: Oxygen isotope ratios (δ18 O values) of fish otoliths (ear bones) are valuable geochemical tracers of water conditions and thermal life history. Delta Smelt (Hypomesus transpacificus) are osmerid forage fish endemic to the San Francisco Estuary, California, USA, that are on the verge of extinction. These fish exhibit a complex life history that allows them to survive in a dynamic estuarine environment; however, a rapidly warming climate threatens this thermally sensitive species. Here we quantify the accuracy and precision of using δ18 O values in otoliths to reconstruct the thermal life histories of Delta Smelt. METHODS: Delta Smelt were reared for 360 days using three different water sources with different ambient δ18 Owater values (-8.75‰, -5.28‰, and -4.06‰) and different water temperatures (16.4°C, 16.7°C, 18.7°C, and 20.5°C). Samples were collected after 170 days (n = 28) and 360 days (n = 14) post-hatch. In situ δ18 O values were measured from the core of the otolith to the dorsal edge using secondary ion mass spectrometry (SIMS) to reconstruct temporally resolved thermal life histories. RESULTS: The δ18 Ootolith values for Delta Smelt varied as a linear inverse function of water temperature: 1000 ln α = 18.39 (±0.43, 1SE)(103 TK-1 ) - 34.56 (±1.49, 1SE) and δ18 Ootolith(VPDB) - δ18 Owater (VPDB) = 31.34(±0.09, 1SE) - 0.19(±0.01, 1SE) × T ° C. When the ambient δ18 Owater value is known, this species-specific temperature-dependent oxygen isotope fractionation model facilitated the accurate (0.25°C) and precise (±0.37°C, 2σ) reconstruction of the water temperature experienced by the fish. In contrast, the use of existing general fractionation equations resulted in inaccurate temperature reconstructions. CONCLUSIONS: The species-specific δ18 Ootolith fractionation equation allowed for accurate and precise reconstructions of water temperatures experienced by Delta Smelt. Characterization of ambient δ18 Owater values remains a critical next step for reconstructing thermal life histories of wild Delta Smelt. This tool will provide new insights into habitat utilization, potential thermal refugia, and resilience to future warming for this critically endangered fish.

IsoFishR: An application for reproducible data reduction and analysis of strontium isotope ratios (87Sr/86Sr) obtained via laser-ablation MC-ICP-MS.

The IsoFishR application is a data reduction and analysis tool for laser-ablation strontium isotope data, following common best practices and providing reliable and reproducible results. Strontium isotope ratios (87Sr/86Sr) are a powerful geochemical tracer commonly applied in a wide range of scientific fields and laser-ablation inductively coupled mass spectrometry is considered the method of choice to obtain spatially resolved 87Sr/86Sr isotope ratios from a variety of sample materials. However, data reduction and analyses methods are variable between different research groups and research communities limiting reproducibility between studies. IsoFishR provides a platform to standardize these methods and can be used for both spot and time-resolved line transects. Furthermore, it provides advanced data analysis tools and filters for outlier removal, noise reduction, and visualization of time resolved data. The application can be downloaded from GitHub (https://github.com/MalteWillmes/IsoFishR) and the source code is available, encouraging future development and evolution of this software.

A bioavailable strontium isoscape for Western Europe: A machine learning approach.

Strontium isotope ratios (87Sr/86Sr) are gaining considerable interest as a geolocation tool and are now widely applied in archaeology, ecology, and forensic research. However, their application for provenance requires the development of baseline models predicting surficial 87Sr/86Sr variations ("isoscapes"). A variety of empirically-based and process-based models have been proposed to build terrestrial 87Sr/86Sr isoscapes but, in their current forms, those models are not mature enough to be integrated with continuous-probability surface models used in geographic assignment. In this study, we aim to overcome those limitations and to predict 87Sr/86Sr variations across Western Europe by combining process-based models and a series of remote-sensing geospatial products into a regression framework. We find that random forest regression significantly outperforms other commonly used regression and interpolation methods, and efficiently predicts the multi-scale patterning of 87Sr/86Sr variations by accounting for geological, geomorphological and atmospheric controls. Random forest regression also provides an easily interpretable and flexible framework to integrate different types of environmental auxiliary variables required to model the multi-scale patterning of 87Sr/86Sr variability. The method is transferable to different scales and resolutions and can be applied to the large collection of geospatial data available at local and global levels. The isoscape generated in this study provides the most accurate 87Sr/86Sr predictions in bioavailable strontium for Western Europe (R2 = 0.58 and RMSE = 0.0023) to date, as well as a conservative estimate of spatial uncertainty by applying quantile regression forest. We anticipate that the method presented in this study combined with the growing numbers of bioavailable 87Sr/86Sr data and satellite geospatial products will extend the applicability of the 87Sr/86Sr geo-profiling tool in provenance applications.