Sardines at a junction: Seascape genomics reveals ecological and oceanographic drivers of variation in the NW Mediterranean Sea.

By evaluating genetic variation across the entire genome, one can address existing questions in a novel way while raising new ones. The latter includes how different local environments influence adaptive and neutral genomic variation within and among populations, providing insights into local adaptation of natural populations and their responses to global change. Here, under a seascape genomic approach, ddRAD data of 4609 single nucleotide polymorphisms (SNPs) from 398 sardines (Sardina pilchardus) collected in 11 Mediterranean and one Atlantic site were generated. These were used along with oceanographic and ecological information to detect signals of adaptive divergence with gene flow across environmental gradients. The studied sardines constitute two clusters (FST  = 0.07), a pattern attributed to outlier loci, highlighting putative local adaptation. The trend in the number of days with sea surface temperature above 19°C, a critical threshold for successful sardine spawning, was crucial at all levels of population structuring with implications on the species' key biological processes. Outliers link candidate SNPs to the region's environmental heterogeneity. Our findings provide evidence for a dynamic equilibrium in which population structure is maintained by physical and ecological factors under the opposing influences of migration and selection. This dynamic in a natural system warrants continuous monitoring under a seascape genomic approach that might benefit from a temporal and more detailed spatial dimension. Our results may contribute to complementary studies aimed at providing deeper insights into the mechanistic processes underlying population structuring. Those are key to understanding and predicting future changes and responses of this highly exploited species in the face of climate change.

Organophosphate ester plasticizers in edible fish from the Mediterranean Sea: Marine pollution and human exposure.

Concentrations of organophosphate esters (OPEs) plasticizers were analysed in the present study. Fifty-five fish samples belonging to three highly commercial species, European sardine (Sardina pilchardus), European anchovy (Engraulis encrasicolus), and European hake (Merluccius merluccius), were taken from the Western Mediterranean Sea. OPEs were detected in all individuals, except for two hake samples, with concentrations between 0.38 and 73.4 ng/g wet weight (ww). Sardines presented the highest mean value with 20.5 ± 20.1 ng/g ww, followed by anchovies with 14.1 ± 8.91 ng/g ww and hake with 2.48 ± 1.76 ng/g ww. The lowest OPE concentrations found in hake, which is a partial predator of anchovy and sardine, and the higher δ15N values (as a proxy of trophic position), may indicate the absence of OPEs biomagnification. Eleven out of thirteen tested OPEs compounds were detected, being diphenyl cresyl phosphate (DCP) one of the most frequently detected in all the species. The highest concentration values were obtained for tris(1,3-dichloro-2-propyl) phosphate (TDClPP), trihexyl phosphate (THP), and tris(2-butoxyethyl) phosphate (TBOEP), for sardines, anchovies, and hakes, respectively. The human health risk associated with the consumption of these fish species showing that their individual consumption would not pose a considerable threat to public health regarding OPE intake.

Trophic niche overlap between round sardinella (Sardinella aurita) and sympatric pelagic fish species in the Western Mediterranean.

The northward expansion of round sardinella (Sardinella aurita) in the Mediterranean Sea, together with declines and fluctuations in biomass and landings of European sardine (Sardina pilchardus) and anchovy (Engraulis encrasicolus) observed in recent decades, may suggest potential inter-specific competition in the pelagic domain. The coexistence of sympatric zooplanktivorous fish species might therefore be exposed in part to trophic niche overlap and competition for food. Combining visual diet characterization under the microscope with DNA metabarcoding from stomach contents of fish collected in spring results show that predation on relatively large krill is equally important for sardinella than for the other two niche overlapping species. Furthermore, an important overlap is found in their isotopic niche, especially with anchovy, using nitrogen (δ15N) and carbon (δ13C) stable isotopes in muscle tissue. In fact, the three fish species are able to feed effectively in the whole prey size spectrum available during the sampled season, from the smallest diatoms and copepods to the larger prey (i.e., decapods and euphausiids), including fish larvae. Moreover, effective predation upon other large prey like siphonophores, which is observed only when multi-proxy analyses in stomach contents are applied, might also be relevant in the diet of sardinella. The overlapping diet composition in spring, together with the effective use of food resource by sardinella, can be of special interest in potential future scenarios with warmer water temperature leading to lower zooplankton and/or higher jellyfish availability, where sardinella may take advantage over other species due to its feeding plasticity.

Energy content of anchovy and sardine using surrogate calorimetry methods.

European anchovy (Engraulis encrasicolus) and sardine (Sardina pilchardus) are crucial species for the marine ecosystem of the Northwestern Mediterranean Sea. They account for a high percentage of fish landings and they represent an important economic income for the fishery sector. Concerns over their stock status are rising in recent years as biomass, growth, reproductive capacity, and body condition of both species are declining, with latitudinal variations. Therefore, there is an urgent need for a body condition monitoring scheme. Energy storage variability has important implications for both fish recruitment and population structure. Direct condition indices, such as energy density (ED) with bomb calorimetry, are highly reliable for measuring the energy content, but time-consuming. Alternatively, fatmeter analysis and relative condition index (Kn) have been proposed as effective indirect methods. The aim of this study is to test the application of fatmeter as a surrogate of bomb calorimetry to infer the energy content of sardine and anchovy. To validate its use, fatmeter values were compared with both ED and Kn values. Individuals of both species were sampled monthly for a year in order to assess seasonal variations in energy content. Our results highlight that fatmeter measurements are strongly correlated with calorimetry ED for sardine, while a weaker but significant correlation was found for anchovy. The observed differences between the two species are related to their breeding strategies. Based on this study, Kn cannot be considered a good proxy of the energy density of sardine, in particular during the resting period. By contrast, fatmeter analysis appears to be a faster and suitable method to evaluate the energy content of both species routinely. In addition, we provide a linear model to infer ED from fatmeter values for both small pelagic fish. Eventually, these findings could be used to implement body condition monitoring protocols and boost continuous large-scale monitoring.

Latitudinal changes in the trophic structure of benthic coastal food webs along the Antarctic Peninsula.

Stable isotopes of C and N have been used to assess the effect of the duration of the sea ice season on the structure of benthic, Antarctic coastal food webs. Samples of suspended particulate organic matter, macroalgae and macroinvertebrates were collected at five subtidal rocky sites across a latitudinal gradient along the western Antarctic Peninsula and among the South Shetland Islands. We tested the hypotheses that trophic positions of omnivores decrease, and food web structure becomes more redundant at higher latitudes. A latitudinal shift in the isotope baseline was detected for both δ13C and δ15N, but the trophic positions of macroinvertebrates and their relative positions along the δ13C axis and were basically constant across sites, even after rescaling stable isotope ratios to account for shifts in the baseline. Although the northernmost and southernmost study sites differed significantly in most of the metrics of the food web structure, changes with latitude and the duration of the sea ice season were non-monotonic. Highly productive phytoplankton blooms during the ice-free season at Esperanza Bay and Rothera Point may explain the observed pattern, as they result in a massive supply of planktonic organic matter to the detritus bank in the seabed and, hence, shorten the C range and increase trophic redundancy. If this hypothesis is correct, the intensity of the summer phytoplankton bloom can be as relevant for the structure of the benthic marine food web as the duration of the sea ice season.

Main drivers of spatial change in the biomass of commercial species between summer and winter in the NW Mediterranean Sea.

There is a general lack of information related to the spatial structure and functioning of marine ecosystems considering seasonality. Here, we modeled the biomass distribution of eight commercial marine species in the northwestern Mediterranean Sea during winter and summer. We hypothesised that the seasonal differences of the water column and the spatial heterogeneity of oceanographic conditions in the study area could result in seasonal variations on the species biomass distributions. We employed a Bayesian hierarchical species distribution modelling approach (B-SDM) with data from two experimental trawl surveys to analyse which are the significant drivers in each season. Our results showed that bathymetry, temperature and fishing patterns are important variables explaining the species spatial biomass distributions. Furthermore, we found seasonal differentiation in the spatial distribution of biomass for all the studied species. Our results provide essential knowledge about the seasonal distributions of key species in the Mediterranean Sea, with important management implications.

A review of the combined effects of climate change and other local human stressors on the marine environment.

Climate change (CC) is a key, global driver of change of marine ecosystems. At local and regional scales, other local human stressors (LS) can interact with CC and modify its effects on marine ecosystems. Understanding the response of the marine environment to the combined effects of CC and LS is crucial to inform marine ecosystem-based management and planning, yet our knowledge of the potential effects of such interactions is fragmented. At a global scale, we explored how cumulative effect assessments (CEAs) have addressed CC in the marine realm and discuss progress and shortcomings of current approaches. For this we conducted a systematic review on how CEAs investigated at different levels of biological organization ecological responses, functional aspects, and the combined effect of CC and HS. Globally, the effects of 52 LS and of 27 CC-related stressors on the marine environment have been studied in combination, such as industrial fisheries with change in temperature, or sea level rise with artisanal fisheries, marine litter, change in sediment load and introduced alien species. CC generally intensified the effects of LS at species level. At trophic groups and ecosystem levels, the effects of CC either intensified or mitigated the effects of other HS depending on the trophic groups or the environmental conditions involved, thus suggesting that the combined effects of CC and LS are context-dependent and vary among and within ecosystems. Our results highlight that large-scale assessments on the spatial interaction and combined effects of CC and LS remain limited. More importantly, our results strengthen the urgent need of CEAs to capture local-scale effects of stressors that can exacerbate climate-induced changes. Ultimately, this will allow identifying management measures that aid counteracting CC effects at relevant scales.

A trophic latitudinal gradient revealed in anchovy and sardine from the Western Mediterranean Sea using a multi-proxy approach.

This work combines state-of-the-art methods (DNA metabarcoding) with classic approaches (visual stomach content characterization and stable isotope analyses of nitrogen (δ15N) and carbon (δ13C)) to investigate the trophic ecology of anchovy (Engraulis encrasicolus) and sardine (Sardina pilchardus) at high taxonomic and spatial resolution in the Western Mediterranean Sea. Gut contents observed are in accordance with the dietary plasticity generally described for anchovy and sardine, suggesting a diet related to the opportunistic ingestion of available prey in a certain area and/or time. Genetic tools also showed modest inter-specific differences regarding ingested species. However, inter-specific and intra-specific differences in ingested prey frequencies and prey biomass reflected a latitudinal signal that could indicate a more effective predation on large prey like krill by anchovy versus sardine, as well as a generalized higher large prey ingestion by both species southwards. In fact, both species presented lower δ15N in the northernmost area. This latitudinal gradient indicates changes in the trophic ecology of anchovy and sardine that coincide with previously described better biological conditions for fish in the southern part of the study area as well as higher landings of both species in recent years.

Ingestion of microplastics and occurrence of parasite association in Mediterranean anchovy and sardine.

We quantified the incidence of microplastics in the gut contents of the European sardine (Sardina pilchardus) and anchovy (Engraulis encrasicolus) in the Northwestern Mediterranean Sea and tested which variables influence this abundance, including the prevalence of parasites (i.e., trematoda larvae and nematodes). We detected a 58% occurrence of microplastics ingestion in sardines and a 60% in anchovies. With respect to sardines, the individuals with lower body conditions were found to have the highest microplastics ingestion probabilities, whereas in anchovies such probabilities were observed in individuals with higher gonadosomatic indices and smaller size. The areas with the highest microplastics ingestion probabilities were the Gulf of Alicante for sardines and the Gulf of Lion - Ebro Delta for anchovies. Both species showed a positive relationship between parasites and microplastics ingestion. These results highlight that both parasitism and ingestion of microplastics are concerns for the health of marine stocks and human consumers.