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.

Variation in the levels of anisakid infection in the European anchovy Engraulis encrasicolus (Linnaeus) from the Bay of Biscay during the period 2000-2023 (ICES Subarea 8).

The European anchovy Engraulis encrasicolus is one of the most important commercial species in the Bay of Biscay (ICES Subarea 8), and our analysis focused on the analysis of the temporal mean abundance, prevalence, and intensity of Anisakis spp. larvae species in anchovies from ICES Subarea 8 in the years 2000, 2001, 2014-2016, and 2019-2023. Prevalence in adult individuals of anchovy was only 1% in 2000 but increased to 90% in 2014. Since 2015, the prevalence has decreased, and the number of individuals affected in 2023 accounted for 17.6%. The mean abundance showed a similar trend, with a peak of 3.79 nematodes/anchovy in 2014, falling to 0.21 in 2023. The species A. simplex sensu stricto and A. pegreffii were identified by PCR/SANGER sequencing and PCR/RLFP techniques in 2019 and 2020. Anisakis simplex (s.s.) was the most abundant species and, according to the results returned by these two techniques, it accounted for an average of 62.4% and 52.1% of total nematodes in 2019 and 2020, respectively. The results of studies monitoring infection levels in anchovies showed that the mean abundance and prevalence changed over the course of the study period and that the proportion of different species of Anisakis is also subject to variation from year to year.

Global mesozooplankton communities show lower connectivity in deep oceanic layers.

Mesozooplankton is a key component of the ocean, regulating global processes such as the carbon pump, and ensuring energy transfer from lower to higher trophic levels. Yet, knowledge on mesozooplankton diversity, distribution and connectivity at global scale is still fragmented. To fill this gap, we applied DNA metabarcoding to mesozooplankton samples collected during the Malaspina-2010 circumnavigation expedition across the Atlantic, Indian and Pacific oceans from the surface to bathypelagic depths. We highlight the still scarce knowledge on global mesozooplankton diversity and identify the Indian Ocean and the deep sea as the oceanic regions with the highest proportion of hidden diversity. We report no consistent alpha-diversity patterns for mesozooplankton at a global scale, neither across vertical nor horizontal gradients. However, beta-diversity analysis suggests horizontal and vertical structuring of mesozooplankton communities mostly attributed to turnover and reveals an increase in mesozooplankton beta-diversity with depth, indicating reduced connectivity at deeper layers. Additionally, we identify a water mass type-mediated structuring of mesozooplankton bathypelagic communities instead of an oceanic basin-mediated as observed at upper layers. This suggests limited dispersal at deep ocean layers, most likely due to weaker currents and lower mixing of water mass types, thus reinforcing the importance of oceanic currents and barriers to dispersal in shaping global plankton communities.

Revisiting the mercury cycle in marine sediments: A potential multifaceted role for Desulfobacterota.

Marine sediments impacted by urban and industrial pollutants are typically exposed to reducing conditions and represent major reservoirs of toxic mercury species. Mercury methylation mediated by anaerobic microorganisms is favored under such conditions, yet little is known about potential microbial mechanisms for mercury detoxification. We used culture-independent (metagenomics, metabarcoding) and culture-dependent approaches in anoxic marine sediments to identify microbial indicators of mercury pollution and analyze the distribution of genes involved in mercury reduction (merA) and demethylation (merB). While none of the isolates featured merB genes, 52 isolates, predominantly affiliated with Gammaproteobacteria, were merA positive. In contrast, merA genes detected in metagenomes were assigned to different phyla, including Desulfobacterota, Actinomycetota, Gemmatimonadota, Nitrospirota, and Pseudomonadota. This indicates a widespread capacity for mercury reduction in anoxic sediment microbiomes. Notably, merA genes were predominately identified in Desulfobacterota, a phylum previously associated only with mercury methylation. Marker genes involved in the latter process (hgcAB) were also mainly assigned to Desulfobacterota, implying a potential central and multifaceted role of this phylum in the mercury cycle. Network analysis revealed that Desulfobacterota were associated with anaerobic fermenters, methanogens and sulfur-oxidizers, indicating potential interactions between key players of the carbon, sulfur and mercury cycling in anoxic marine sediments.

Analysis of potential drivers of spatial and temporal changes in anisakid larvae infection levels in European hake, Merluccius merluccius (L.), from the North-East Atlantic fishing grounds.

We analysed the spatial and temporal variability of Anisakis larvae infection in hake (Merluccius merluccius) from the North-East Atlantic from 1998 to 2020 and the potential drivers (i.e., environmental and host abundance) of such variation. The results showed that hake from separate sea areas in the North Atlantic have marked differences in temporal abundance levels. Hake larger than 60 cm were all parasitized in all ICES (International Council for the Exploration of the Sea) subareas 6, 7, and 8. The belly flaps were the most parasitized parts of the flesh, accounting for 92% of the total. Individuals of Anisakis simplex, Anisakis pegreffii, Anisakis spp. and a hybrid of Anisakis simplex × pegreffii were genetically identified, and Anisakis simplex as the most abundant (88-100%). An ecological niche model of Anisakis occurrence in fishes in the NE Atlantic was built to define the thermal optimum and environmental ranges for salinity, depth, chlorophyll concentration, and diffuse attenuation. The temporal variability of anisakid infection in fishes in the last two decades indicated an increase in the NE Atlantic at an annual rate of 31.7 nematodes per total number of specimens examined per year. This rise in infection levels could be triggered by the increase in intermediate host fish stocks, especially hake in the area.

Infection Rate in Seabasses Fed with Viscera Parasitised by Anisakid Larvae.

PURPOSE: It has been suggested that the removal of infected viscera on board is responsible for the high prevalence of anisakid larvae present in wild fish species. The aim of this work is to assess the re-infection capacity of anisakid larvae in European seabasses, a natural host species for the parasite by feeding with pieces of parasitised hake liver under controlled experimental conditions. METHODS: To prove this potential link between manipulation and re-infestation, 50 farmed seabasses free of anisakid nematodes were fed with fresh hake liver pieces naturally infested with anisakid larvae. RESULTS: After digestion periods from 4 to 21 days, the seabasses showed a prevalence of Anisakis of 6%, and a low retention rate of 0.11 larvae/seabass after four days' digestion, and 0.0021 after 21 day digestion. Two nematodes were found in the intestine and in the visceral cavity, and 13 Anisakis were found partially digested in the stomach of one same individual after 4 day digestion. Results showed that only a small number of Anisakis ingested with the viscera were able to reinfect the seabasses, as most of the larvae seemed to be quickly digested or defecated. CONCLUSION: it seems that the availability of larvae that could re-enter the life cycle and re-infect a fish after the removal and discarding the infected viscera on board could be much less important than commonly believed.

Evidence of stock connectivity, hybridization, and misidentification in white anglerfish supports the need of a genetics-informed fisheries management framework.

Understanding population connectivity within a species as well as potential interactions with its close relatives is crucial to define management units and to derive efficient management actions. However, although genetics can reveal mismatches between biological and management units and other relevant but hidden information such as species misidentification or hybridization, the uptake of genetic methods by the fisheries management process is far from having been consolidated. Here, we have assessed the power of genetics to better understand the population connectivity of white (Lophius piscatorius) and its interaction with its sister species, the black anglerfish (Lophius budegassa). Our analyses, based on thousands of genome-wide single nucleotide polymorphisms, show three findings that are crucial for white anglerfish management. We found (i) that white anglerfish is likely composed of a single panmictic population throughout the Northeast Atlantic, challenging the three-stock based management, (ii) that a fraction of specimens classified as white anglerfish using morphological characteristics are genetically identified as black anglerfish (L. budegassa), and iii) that the two Lophius species naturally hybridize leading to a population of hybrids of up to 20% in certain areas. Our results set the basics for a genetics-informed white anglerfish assessment framework that accounts for stock connectivity, revises and establishes new diagnostic characters for Lophius species identification, and evaluates the effect of hybrids in the current and future assessments of the white anglerfish. Furthermore, our study contributes to provide additional evidence of the potentially negative consequences of ignoring genetic data for assessing fisheries resources.

A microbial mandala for environmental monitoring: Predicting multiple impacts on estuarine prokaryote communities of the Bay of Biscay.

Routine monitoring of benthic biodiversity is critical for managing and understanding the anthropogenic impacts on marine, transitional and freshwater ecosystems. However, traditional reliance on morphological identification generally makes it cost-prohibitive to increase the scale of monitoring programmes. Metabarcoding of environmental DNA has clear potential to overcome many of the problems associated with traditional monitoring, with prokaryotes and other microorganisms showing particular promise as bioindicators. However, due to the limited knowledge regarding the ecological roles and responses of environmental microorganisms to different types of pressure, the use of de novo approaches is necessary. Here, we use two such approaches for the prediction of multiple impacts present in estuaries and coastal areas of the Bay of Biscay based on microbial communities. The first (Random Forests) is a machine learning method while the second (Threshold Indicator Taxa Analysis and quantile regression splines) is based on de novo identification of bioindicators. Our results show that both methods overlap considerably in the indicator taxa identified, but less for sequence variants. Both methods also perform well in spite of the complexity of the studied ecosystem, providing predictive models with strong correlation to reference values and fair to good agreement with ecological status groups. The ability to predict several specific types of pressure is especially appealing. The cross-validated models and biotic indices developed can be directly applied to predict the environmental status of estuaries in the same geographical region, although more work is needed to evaluate and improve them for use in new regions or habitats.

Pan-regional marine benthic cryptobiome biodiversity patterns revealed by metabarcoding Autonomous Reef Monitoring Structures.

Autonomous Reef Monitoring Structures (ARMS) have been applied worldwide to characterize the critical yet frequently overlooked biodiversity patterns of marine benthic organisms. In order to disentangle the relevance of environmental factors in benthic patterns, here, through standardized metabarcoding protocols, we analyse sessile and mobile (<2 mm) organisms collected using ARMS deployed across six regions with different environmental conditions (3 sites × 3 replicates per region): Baltic, Western Mediterranean, Adriatic, Black and Red Seas, and the Bay of Biscay. A total of 27,473 Amplicon Sequence Variants (ASVs) were observed ranging from 1,404 in the Black Sea to 9,958 in the Red Sea. No ASVs were shared among all regions. The highest number of shared ASVs was between the Western Mediterranean and the Adriatic Sea (116) and Bay of Biscay (115). Relatively high numbers of ASVs (103), mostly associated with the genus Amphibalanus, were also shared between the lower salinity seas (Baltic and Black Seas). We found that compositional differences in spatial patterns of rocky-shore benthos are determined slightly more by dispersal limitation than environmental filtering. Dispersal limitation was similar between sessile and mobile groups, while the sessile group had a larger environmental niche breadth than the mobile group. Further, our study can provide a foundation for future evaluations of biodiversity patterns in the cryptobiome, which can contribute up to 70% of the local biodiversity.

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.

Novel Vibrio spp. Strains Producing Omega-3 Fatty Acids Isolated from Coastal Seawater.

Omega-3 long-chain polyunsaturated fatty acids (LC-PUFAs), such as eicosapentaenoic acid (EPA) (20:5n-3) and docosahexaenoic acid (DHA) (22:6n-3), are considered essential for human health. Microorganisms are the primary producers of omega-3 fatty acids in marine ecosystems, representing a sustainable source of these lipids, as an alternative to the fish industry. Some marine bacteria can produce LC-PUFAs de novo via the Polyunsaturated Fatty Acid (Pfa) synthase/ Polyketide Synthase (PKS) pathway, which does not require desaturation and elongation of saturated fatty acids. Cultivation-independent surveys have revealed that the diversity of microorganisms harboring a molecular marker of the pfa gene cluster (i.e., pfaA-KS domain) is high and their potential distribution in marine systems is widespread, from surface seawater to sediments. However, the isolation of PUFA producers from marine waters has been typically restricted to deep or cold environments. Here, we report a phenotypic and genotypic screening for the identification of omega-3 fatty acid producers in free-living bacterial strains isolated from 5, 500, and 1000 m deep coastal seawater from the Bay of Biscay (Spain). We further measured EPA production in pelagic Vibrio sp. strains collected at the three different depths. Vibrio sp. EPA-producers and non-producers were simultaneously isolated from the same water samples and shared a high percentage of identity in their 16S rRNA genes, supporting the view that the pfa gene cluster can be horizontally transferred. Within a cluster of EPA-producers, we found intraspecific variation in the levels of EPA synthesis for isolates harboring different genetic variants of the pfaA-KS domain. The maximum production of EPA was found in a Vibrio sp. strain isolated from a 1000 m depth (average 4.29% ± 1.07 of total fatty acids at 10 °C, without any optimization of culturing conditions).

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.

Solute carrier family 2 member 1 is involved in the development of nonalcoholic fatty liver disease.

UNLABELLED: Susceptibility to develop nonalcoholic fatty liver disease (NAFLD) has genetic bases, but the associated variants are uncertain. The aim of the present study was to identify genetic variants that could help to prognose and further understand the genetics and development of NAFLD. Allele frequencies of 3,072 single-nucleotide polymorphisms (SNPs) in 92 genes were characterized in 69 NAFLD patients and 217 healthy individuals. The markers that showed significant allele-frequency differences in the pilot groups were subsequently studied in 451 NAFLD patients and 304 healthy controls. Besides this, 4,414 type 2 diabetes mellitus (T2DM) cases and 4,567 controls were genotyped. Liver expression of the associated gene was measured and the effect of its potential role was studied by silencing the gene in vitro. Whole genome expression, oxidative stress (OS), and the consequences of oleic acid (OA)-enriched medium on lipid accumulation in siSLC2A1-THLE2 cells were studied by gene-expression analysis, dihydroethidium staining, BODIPY, and quantification of intracellular triglyceride content, respectively. Several SNPs of SLC2A1 (solute carrier family 2 [facilitated glucose transporter] member 1) showed association with NAFLD, but not with T2DM, being the haplotype containing the minor allele of SLC2A1 sequence related to the susceptibility to develop NAFLD. Gene-expression analysis demonstrated a significant down-regulation of SLC2A1 in NAFLD livers. Enrichment functional analyses of transcriptome profiles drove us to demonstrate that in vitro silencing of SLC2A1 induces an increased OS activity and a higher lipid accumulation under OA treatment. CONCLUSIONS: Genetic variants of SLC2A1 are associated with NAFLD, and in vitro down-regulation of this gene promotes lipid accumulation. Moreover, the oxidative response detected in siSLC2A1-THLE2 cells corroborated the antioxidant properties previously related to this gene and linked the most representative clinical characteristics of NAFLD patients: oxidative injury and increased lipid storage.

Fine mapping of a major histocompatibility complex in ankylosing spondylitis: association of the HLA-DPA1 and HLA-DPB1 regions.

OBJECTIVE: To investigate the potential association of major histocompatibility complex (MHC) markers other than HLA-B27 with ankylosing spondylitis (AS). METHODS: A total of 603 patients with AS and 542 healthy control subjects, all of whom were HLA-B27 positive, were selected for this study based on clinical criteria. First, high-density genotyping across the MHC region (2,360 single-nucleotide polymorphisms [SNPs]) was performed in a cohort of 191 patients and 241 control subjects. After a fine-mapping study, 5 SNPs from the HLA-DPA1/DPB1 region were validated in a second cohort of 412 patients with AS and 301 healthy control subjects. RESULTS: Seventeen SNPs located within or near the HLA-DPA1 and HLA-DPB1 loci showed association with AS (P = 1.38 × 10⁻5 to 0.05). In addition, multimarker tests, both linkage disequilibrium and sliding windows, showed association of some groups of adjacent SNPs within the HLA-DPA1/DPB1 region with AS (P = 1.0 × 10⁻4 to 3.96 × 10⁻7). We validated the association by genotyping 5 SNPs from the DPA1/DPB1 region in an additional cohort and obtained P values from 6.42 × 10⁻5 to 0.01 in the analysis of the combined cohorts. Subtyping analysis of HLA-DPA1 and HLA-DPB1 showed that HLA-DPA1*01:03, A1*02:01, and B1*13:01 were the subtypes most susceptible to AS. CONCLUSION: HLA markers and linkage disequilibrium blocks near HLA-DPA1 and HLA-DPB1 are statistically associated with AS. We identified a region located around the HLA-DPA1 and HLA-DPB1 loci associated with AS, another region within the MHC that is different from HLA-B27.