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.

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.

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.

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.

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.

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.

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).

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

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⁻⁵ 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⁻⁴ to 3.96 × 10⁻⁷). 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⁻⁵ 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.

[Knowledge and attitudes of practice of health workers in Bizkaia health centers toward diet and cholesterol]

OBJECTIVE

Analyze the extent to which criteria and practice of the Sanitary Personnel approximate the recommendations of the Consensus on the Control of Hypercholesterolemia in Spain, find out how extensive their dietary knowledge in and how skillful they are to interview and counsel a patient with hyperlipidemia.

DESIGN

Transversal study with simple aleatory sampling using a questionnaire. PARTICIPANTS AND LOCATION: 101 physicians and 95 nurses working in Primary Health Care Teams in Bizkaia, during the monthly of June 1991; deputies for periods inferior to 3 months were excluded.

MEASUREMENTS AND MAIN RESULTS

82% of physicians (P) counselled dietetic treatment for cholesterolemic levels coinciding with those set in the Consensus and 70% demanded excessive analytical controls for hypercholesterolemic patients under treatment already normalized; 37% (P) and 12% of nurses (N) knew that the polyunsaturated/saturated fat ratio is the main dietetic factor; various errors were detected in questions relating to olive oil, unsaturated fat in fish and fibre values of legume; accurate responses between 13-54% (N) and 23-77% (P).

CONCLUSIONS

The "mediterranean diet" is not sufficiently known; our ability for the dietetic interview and personalized counsel needs to be improved. In the light of the above and taking into accounted the gaps and deviations existing between our practice and the Consensus and its applications, we consider a better diffusion of the cited document is necessary among Primary Health Care Teams.