Genomic basis for early-life mortality in sharpsnout seabream.

Mortality at early life stages of fishes is common in nature and can be shaped by stochastic and selective processes. Selective mortality has rarely been assessed in natural conditions but can now be studied by combining genomic data with information on different life stages that realates to fitness. Here we investigate selective mortality between settlers and six-month survivors of the sharpsnout seabream by genotype-phenotype/environmental association studies in three localities along a geographic gradient. We gathered information on 105 individuals at 85,031 SNPs, obtained from individual based 2b-RAD libraries, as well as 9 phenotypic and environmental variables derived from individual otolith readings. We found common signals across localities for potential selection, such as lower survival rates for individuals hatching earlier, growing faster and experiencing higher temperatures during their planktonic phase. We identified 122 loci with parallel significant association to phenotypic and environmental variables. Importantly, one of these loci mapped to the exonic region of the il20rb, a gene involved in immune response, in the phylogenetically closest reference genome, showing parallel frequency changes in non-synonymous mutations in the three studied populations. Further temporal assessments are needed to understand how polymorphisms that are key to selective mortality are maintained.

Temperature reduces fish dispersal as larvae grow faster to their settlement size.

As species struggle to cope with rising ocean temperatures, temperate marine assemblages are facing major reorganization. Many benthic species have a brief but critical period dispersing through the plankton, when they are particularly susceptible to variations in temperature. Impacts of rising temperatures can thus ripple through the population with community-wide consequences. However, responses are highly species-specific, making it difficult to discern assemblage-wide patterns in the life histories of different fish species. Here, we evaluate the responses to temperature in the early life histories of several fish species using otolith reconstructive techniques. We also assess the consequences of future warming scenarios to this assemblage. We sampled recent settlers of nine common species across a temperature gradient in the Mediterranean Sea and obtained environmental data for each individual. Using otolith microstructure, we measured early life traits including pelagic larval duration (PLD), growth rate, settlement size, hatching and settlement dates. We used a GLM framework to examine how environmental variables influenced early life-history parameters. We show that increasing temperature results in considerable reduction in the dispersal potential of temperate fish. We find a nearly universal, assemblage-wide decline in pelagic larval duration (PLD) of between 10% and 25%. This was because, with increasing temperature, larvae grew quicker to their settlement size. Settlement size itself was less affected by temperature and appears to be an ontogenetically fixed process. Given current estimates of ocean warming, there could be an assemblage-wide reduction in larval dispersal of up to 50 km across the Mediterranean, reducing connectivity and potentially isolating populations as waters warm.

Warming intensifies the interaction between the temperate seagrass Posidonia oceanica and its dominant fish herbivore Sarpa salpa.

Apart from directly influencing individual life histories of species, climate change is altering key biotic interactions as well, causing community processes to unravel. With rising temperatures, disruptions to producer-consumer relationships can have major knock-on effects, particularly when the producer is a habitat-forming species. We studied how sea surface temperature (SST) modifies multiple pathways influencing the interaction between the foundational seagrass species, Posidonia oceanica, and its main consumer, the fish Sarpa salpa in the Mediterranean Sea. We used a combination of a field-based temperature gradient approaches and experimental manipulations to assess the effect of temperature on seagrass performance (growth) and fish early life history (larval development) as well as on the interaction itself (seagrass palatability and fish foraging activity). Within the range of temperatures assessed, S. salpa larvae grew slightly faster at warmer conditions but maintained their settlement size, resulting in a relatively small reduction in pelagic larval duration (PLD) and potentially reducing dispersion. Under warmer conditions (>24 °C), P. oceanica reduced its growth rate considerably and seemed to display fewer deterring mechanisms as indicated by a disproportionate consumption in choice experiments. However, our field-based observations along the temperature gradient showed no change in fish foraging time, or in other aspects of feeding behaviour. As oceans warm, our results indicate that, while S. salpa may show little change in early life history, its preference towards P. oceanica might increase, which, together with reduced seagrass growth, could considerably intensify the strength of herbivory. It is unclear if P. oceanica meadows can sustain such an intensification, but it will clearly add to the raft of pressures this threatened ecosystem already faces from global and local environmental change.

Individual-based population genomics reveal different drivers of adaptation in sympatric fish.

Connectivity and local adaptation are two contrasting evolutionary forces highly influencing population structure. To evaluate the impact of early-life traits and environmental conditions on genetic structuring and adaptation, we studied two sympatric fish species in the Western Mediterranean Sea: Symphodus tinca and S. ocellatus. We followed an individual-based approach and measured early-life history traits from otolith readings, gathered information on environmental variables and obtained genome-wide markers from genotyping-by-sequencing (GBS). The two species presented contrasting population structure across the same geographic gradient, with high and significant population differentiation in S. ocellatus, mostly determined by oceanographic fronts, and low differentiation and no front effect in S. tinca. Despite their different levels of genetic differentiation, we identified in both species candidate regions for local adaptation by combining outlier analysis with environmental and phenotypic association analyses. Most candidate loci were associated to temperature and productivity in S. ocellatus and to temperature and turbulence in S. tinca suggesting that different drivers may determine genomic diversity and differentiation in each species. Globally, our study highlights that individual-based approach combining genomic, environmental and phenotypic information is key to identify signals of selection and the processes mediating them.

Helping decision making for reliable and cost-effective 2b-RAD sequencing and genotyping analyses in non-model species.

High-throughput sequencing has revolutionized population and conservation genetics. RAD sequencing methods, such as 2b-RAD, can be used on species lacking a reference genome. However, transferring protocols across taxa can potentially lead to poor results. We tested two different IIB enzymes (AlfI and CspCI) on two species with different genome sizes (the loggerhead turtle Caretta caretta and the sharpsnout seabream Diplodus puntazzo) to build a set of guidelines to improve 2b-RAD protocols on non-model organisms while optimising costs. Good results were obtained even with degraded samples, showing the value of 2b-RAD in studies with poor DNA quality. However, library quality was found to be a critical parameter on the number of reads and loci obtained for genotyping. Resampling analyses with different number of reads per individual showed a trade-off between number of loci and number of reads per sample. The resulting accumulation curves can be used as a tool to calculate the number of sequences per individual needed to reach a mean depth ≥20 reads to acquire good genotyping results. Finally, we demonstrated that selective-base ligation does not affect genomic differentiation between individuals, indicating that this technique can be used in species with large genome sizes to adjust the number of loci to the study scope, to reduce sequencing costs and to maintain suitable sequencing depth for a reliable genotyping without compromising the results. Here, we provide a set of guidelines to improve 2b-RAD protocols on non-model organisms with different genome sizes, helping decision-making for a reliable and cost-effective genotyping.