Increasing marine trophic web knowledge through DNA analyses of fish stomach content: a step towards an ecosystem-based approach to fisheries research.

Multispecies and ecosystem models, which are key for the implementation of ecosystem-based approaches to fisheries management, require extensive data on the trophic interactions between marine organisms, including changes over time. DNA metabarcoding, by allowing the simultaneous taxonomic identification of the community present in hundreds of samples, could be used for speeding up large-scale stomach content data collection. Yet, for DNA metabarcoding to be routinely implemented, technical challenges should be addressed, such as the potentially complicated sampling logistics, the detection of a high proportion of predator DNA, and the inability to provide reliable abundance estimations. Here, we present a DNA metabarcoding assay developed to examine the diet of five commercially important fish, which can be feasibly incorporated into routinary samplings. The method is devised to speed up the analysis process by avoiding the stomach dissection and content extraction steps, while preventing the amplification of predator DNA by using blocking primers. Tested in mock samples and in real stomach samples, the method has proven effective and shows great effectiveness discerning diet variations due to predator ecology or prey availability. Additionally, by applying our protocol to mackerel stomachs previously analyzed by visual inspection, we showcase how DNA metabarcoding could complement visually based data by detecting overlooked prey by the visual approach. We finally discuss how DNA metabarcoding-based data can contribute to trophic data collection. Our work reinforces the potential of DNA metabarcoding for the study and monitoring of fish trophic interactions and provides a basis for its incorporation into routine monitoring programs, which will be critical for the implementation of ecosystem-based approaches to fisheries management.

Shrinking body size of European anchovy in the Bay of Biscay.

Decreased body size is often cited as a major response to ocean warming. Available evidence, however, questions the actual emergence of shrinking trends and the prevalence of temperature-driven changes in size over alternative drivers. In marine fish, changes in food availability or fluctuations in abundance, including those due to size-selective fishing, provide compelling mechanisms to explain changes in body size. Here, based on three decades of scientific survey data (1990-2021), we report a decline in the average body size-length and weight-of anchovy, Engraulis encrasicolus L., in the Bay of Biscay. Shrinking was evident in all age classes, from juveniles to adults. Allometric adjustment indicated slightly more pronounced declines in weight than in total length, which is consistent with a change toward a slender body shape. Trends in adult weight were nonlinear, with rates accelerating to an average decline of up to 25% decade-1 during the last two decades. We found a strong association between higher anchovy abundance and reduced juvenile size. The effect of density dependence was less clear later in life, and temperature became the best predictor of declines in adult size. Theoretical analyses based on a strategic model further suggested that observed patterns are consistent with a simultaneous, opposing effect of rising temperatures on accelerating early growth and decreasing adult size as predicted by the temperature-size rule. Macroecological assessment of ecogeographical-Bergmann's and James'-rules in anchovy size suggested that the observed decline largely exceeds intraspecific variation and might be the result of selection. Limitations inherent in the observational nature of the study recommend caution and a continued assessment and exploration of alternative drivers. Additional evidence of a climate-driven regime shift in the region suggests, however, that shrinking anchovy sizes may signal a long-lasting change in the structure and functioning of the Bay of Biscay ecosystem.

Climate regime shifts and biodiversity redistribution in the Bay of Biscay.

Global ocean warming, wave extreme events, and accelerating sea-level rise are challenges that coastal communities must address to anticipate damages in coming decades. The objective of this study is to undertake a time-series analysis of climate change (CC) indicators within the Bay of Biscay, including the Basque coast. We used an integrated and flexible methodology, based on Generalized Additive Mixed Models, to detect trends on 19 indicators (including marine physics, chemistry, atmosphere, hydrology, geomorphology, biodiversity, and commercial species). The results of 87 long-term time series analysed (~512,000 observations), in the last four decades, indicate four groups of climate regime shifts: 1) A gradual shift associated with CC starting in the 1980s, with a warming of the sea surface down to 100 m depth in the bay (0.10-0.25 °C per decade), increase in air temperature and insolation. This warming may have impacted on benthic community redistribution in the Basque coast, favouring warm-water species relative to cold-water species. Weight at age for anchovy and sardine decreased in the last two decades. 2) Deepening of the winter mixed layer depth in the south-eastern bay that probably led to increases in nutrients, surface oxygen, and chlorophyll concentration. Current increases on chlorophyll and zooplankton (i.e., copepods) biomass are contrary to those expected under CC scenarios in the region. 3) Sea-level rise (1.5-3.5 cm per decade since 1990s), associated with CC. 4) Increase of extreme wave height events of 16.8 cm per decade in the south-eastern bay, probably related to stormy conditions in the last decade, with impacts on beach erosion. Estimating accurate rates of sea warming, sea-level rise, extreme events, and foreseeing the future pathways of marine productivity, are key to define the best adaptation measures to minimize negative CC impacts in the region.

Marine water environmental DNA metabarcoding provides a comprehensive fish diversity assessment and reveals spatial patterns in a large oceanic area.

Current methods for monitoring marine fish (including bony fishes and elasmobranchs) diversity mostly rely on trawling surveys, which are invasive, costly, and time-consuming. Moreover, these methods are selective, targeting a subset of species at the time, and can be inaccessible to certain areas. Here, we used environmental DNA (eDNA), the DNA present in the water column as part of shed cells, tissues, or mucus, to provide comprehensive information about fish diversity in a large marine area. Further, eDNA results were compared to the fish diversity obtained in pelagic trawls. A total of 44 5 L-water samples were collected onboard a wide-scale oceanographic survey covering about 120,000 square kilometers in Northeast Atlantic Ocean. A short region of the 12S rRNA gene was amplified and sequenced through metabarcoding generating almost 3.5 million quality-filtered reads. Trawl and eDNA samples resulted in the same most abundant species (European anchovy, European pilchard, Atlantic mackerel, and blue whiting), but eDNA metabarcoding resulted in more detected bony fish and elasmobranch species (116) than trawling (16). Although an overall correlation between fishes biomass and number of reads was observed, some species deviated from the common trend, which could be explained by inherent biases of each of the methods. Species distribution patterns inferred from eDNA metabarcoding data coincided with current ecological knowledge of the species, suggesting that eDNA has the potential to draw sound ecological conclusions that can contribute to fish surveillance programs. Our results support eDNA metabarcoding for broad-scale marine fish diversity monitoring in the context of Directives such as the Common Fisheries Policy or the Marine Strategy Framework Directive.

Including ecosystem descriptors in current fishery data collection programmes to advance towards a holistic monitoring: Seabird abundance attending demersal trawlers.

To develop an ecosystem-based approach to fisheries management a holistic perspective is necessary that goes beyond target species management to preserve ecosystem functioning and, therefore, secure future food availability. To achieve these objectives, current fishery data collection programmes should widen their objectives to include the collection of ecosystem descriptors to effectively take advantage of funding and human resources in relation to fisheries monitoring already in place. From an ecological perspective, fishing discards are food subsidies unnaturally available that can profoundly impact the life history traits and population dynamics of seabirds, as well as community structure. In 2015, we took advantage of the Data Collection Framework (DCF) programme, monitoring the Basque trawling fleet, to monitor seabird abundance associated with trawlers as an additional task to be performed by the observers. The main objectives were (1) to develop a standard protocol from an interdisciplinary expert committee, (2) to obtain baseline information of seabird association with trawlers to be able to track changes and (3) to understand seabird abundance in relation to discard facilitation processes based on environmental and fishing parameters. Based on the developed standard protocol, more than 21 species of seabirds were recorded in 241 fishing hauls in the Bay of Biscay (ICES areas 8abd). This biogeographic area is an important migratory flyway and wintering area, where maximum number of seabirds attending trawlers were recorded between December and April (2015-2019). Based on Generalised Linear Models, seabird abundance was higher in multiple situations: harsher environmental conditions, less available discards, shallower areas, increased period of discard availability, higher number of trawlers simultaneously operating in the same fishing ground, in addition to the months of February and April. The core discarding ground was located in the outer French shelf between the latitudinal range of 44.5°N and 46°N, especially for otter trawl fishing. Our approach puts into value the seabird counts that the observers of the DCF can perform systematically, collecting relevant information on the effect of trawling on other biodiversity components such as seabirds. This information will be critical to respond to the application of the reform of the Common Fisheries Policy regarding the effect of the Landing Obligation that seeks fishing sustainability.

BioTIME: A database of biodiversity time series for the Anthropocene.

MOTIVATION: The BioTIME database contains raw data on species identities and abundances in ecological assemblages through time. These data enable users to calculate temporal trends in biodiversity within and amongst assemblages using a broad range of metrics. BioTIME is being developed as a community-led open-source database of biodiversity time series. Our goal is to accelerate and facilitate quantitative analysis of temporal patterns of biodiversity in the Anthropocene. MAIN TYPES OF VARIABLES INCLUDED: The database contains 8,777,413 species abundance records, from assemblages consistently sampled for a minimum of 2 years, which need not necessarily be consecutive. In addition, the database contains metadata relating to sampling methodology and contextual information about each record. SPATIAL LOCATION AND GRAIN: BioTIME is a global database of 547,161 unique sampling locations spanning the marine, freshwater and terrestrial realms. Grain size varies across datasets from 0.0000000158 km2 (158 cm2) to 100 km2 (1,000,000,000,000 cm2). TIME PERIOD AND GRAIN: BioTIME records span from 1874 to 2016. The minimal temporal grain across all datasets in BioTIME is a year. MAJOR TAXA AND LEVEL OF MEASUREMENT: BioTIME includes data from 44,440 species across the plant and animal kingdoms, ranging from plants, plankton and terrestrial invertebrates to small and large vertebrates. SOFTWARE FORMAT: .csv and .SQL.

Population structure of Atlantic mackerel inferred from RAD-seq-derived SNP markers: effects of sequence clustering parameters and hierarchical SNP selection.

Restriction-site-associated DNA sequencing (RAD-seq) and related methods are revolutionizing the field of population genomics in nonmodel organisms as they allow generating an unprecedented number of single nucleotide polymorphisms (SNPs) even when no genomic information is available. Yet, RAD-seq data analyses rely on assumptions on nature and number of nucleotide variants present in a single locus, the choice of which may lead to an under- or overestimated number of SNPs and/or to incorrectly called genotypes. Using the Atlantic mackerel (Scomber scombrus L.) and a close relative, the Atlantic chub mackerel (Scomber colias), as case study, here we explore the sensitivity of population structure inferences to two crucial aspects in RAD-seq data analysis: the maximum number of mismatches allowed to merge reads into a locus and the relatedness of the individuals used for genotype calling and SNP selection. Our study resolves the population structure of the Atlantic mackerel, but, most importantly, provides insights into the effects of alternative RAD-seq data analysis strategies on population structure inferences that are directly applicable to other species.