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.

Unidirectional trans-Atlantic gene flow and a mixed spawning area shape the genetic connectivity of Atlantic bluefin tuna.

The commercially important Atlantic bluefin tuna (Thunnus thynnus), a large migratory fish, has experienced notable recovery aided by accurate resource assessment and effective fisheries management efforts. Traditionally, this species has been perceived as consisting of eastern and western populations, spawning respectively in the Mediterranean Sea and the Gulf of Mexico, with mixing occurring throughout the Atlantic. However, recent studies have challenged this assumption by revealing weak genetic differentiation and identifying a previously unknown spawning ground in the Slope Sea used by Atlantic bluefin tuna of uncertain origin. To further understand the current and past population structure and connectivity of Atlantic bluefin tuna, we have assembled a unique dataset including thousands of genome-wide single-nucleotide polymorphisms (SNPs) from 500 larvae, young of the year and spawning adult samples covering the three spawning grounds and including individuals of other Thunnus species. Our analyses support two weakly differentiated but demographically connected ancestral populations that interbreed in the Slope Sea. Moreover, we also identified signatures of introgression from albacore (Thunnus alalunga) into the Atlantic bluefin tuna genome, exhibiting varied frequencies across spawning areas, indicating strong gene flow from the Mediterranean Sea towards the Slope Sea. We hypothesize that the observed genetic differentiation may be attributed to increased gene flow caused by a recent intensification of westward migration by the eastern population, which could have implications for the genetic diversity and conservation of western populations. Future conservation efforts should consider these findings to address potential genetic homogenization in the species.

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.

Exploring uncharted territory: new frontiers in environmental DNA for tropical fisheries management.

Tropical ecosystems host a significant share of global fish diversity contributing substantially to the global fisheries sector. Yet their sustainable management is challenging due to their complexity, diverse life history traits of tropical fishes, and varied fishing techniques involved. Traditional monitoring techniques are often costly, labour-intensive, and/or difficult to apply in inaccessible sites. These limitations call for the adoption of innovative, sensitive, and cost-effective monitoring solutions, especially in a scenario of climate change. Environmental DNA (eDNA) emerges as a potential game changer for biodiversity monitoring and conservation, especially in aquatic ecosystems. However, its utility in tropical settings remains underexplored, primarily due to a series of challenges, including the need for a comprehensive barcode reference library, an understanding of eDNA behaviour in tropical aquatic environments, standardized procedures, and supportive biomonitoring policies. Despite these challenges, the potential of eDNA for sensitive species detection across varied habitats is evident, and its global use is accelerating in biodiversity conservation efforts. This review takes an in-depth look at the current state and prospects of eDNA-based monitoring in tropical fisheries management research. Additionally, a SWOT analysis is used to underscore the opportunities and threats, with the aim of bridging the knowledge gaps and guiding the more extensive and effective use of eDNA-based monitoring in tropical fisheries management. Although the discussion applies worldwide, some specific experiences and insights from Indian tropical fisheries are shared to illustrate the practical application and challenges of employing eDNA in a tropical context.

Innovative and practical tools for monitoring and assessing biodiversity status and impacts of multiple human pressures in marine systems.

Human activities at sea can produce pressures and cumulative effects on ecosystem components that need to be monitored and assessed in a cost-effective manner. Five Horizon European projects have joined forces to collaboratively increase our knowledge and skills to monitor and assess the ocean in an innovative way, assisting managers and policy-makers in taking decisions to maintain sustainable activities at sea. Here, we present and discuss the status of some methods revised during a summer school, aiming at better management of coasts and seas. We include novel methods to monitor the coastal and ocean waters (e.g. environmental DNA, drones, imaging and artificial intelligence, climate modelling and spatial planning) and innovative tools to assess the status (e.g. cumulative impacts assessment, multiple pressures, Nested Environmental status Assessment Tool (NEAT), ecosystem services assessment or a new unifying approach). As a concluding remark, some of the most important challenges ahead are assessing the pros and cons of novel methods, comparing them with benchmark technologies and integrating these into long-standing time series for data continuity. This requires transition periods and careful planning, which can be covered through an intense collaboration of current and future European projects on marine biodiversity and ecosystem health.

Trade-offs between reducing complex terminology and producing accurate interpretations from environmental DNA: Comment on "Environmental DNA: What's behind the term?" by Pawlowski et al., (2020).

In a recent paper, "Environmental DNA: What's behind the term? Clarifying the terminology and recommendations for its future use in biomonitoring," Pawlowski et al. argue that the term eDNA should be used to refer to the pool of DNA isolated from environmental samples, as opposed to only extra-organismal DNA from macro-organisms. We agree with this view. However, we are concerned that their proposed two-level terminology specifying sampling environment and targeted taxa is overly simplistic and might hinder rather than improve clear communication about environmental DNA and its use in biomonitoring. This terminology is based on categories that are often difficult to assign and uninformative, and it overlooks a fundamental distinction within eDNA: the type of DNA (organismal or extra-organismal) from which ecological interpretations are derived.

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.

Environmental DNA Metabarcoding: A Promising Tool for Ballast Water Monitoring.

Nonindigenous species are introduced worldwide with ballast water (BW). To prevent further introductions, oceanic BW exchange and BW treatment systems are utilized, but their performance needs to be evaluated. To that aim, characterizing BW communities is essential but usually relies on exhaustive sampling and morphological taxonomic identification, which does not always allow fine-scale taxonomic resolution. Through the analysis of BW samples from 11 vessels arriving to the Chesapeake Bay (USA), we evaluated the potential of environmental DNA (eDNA) metabarcoding for BW monitoring by assessing whether the impact of BW management type could be identified, analyzing the influence of BW sampling access locations on communities, and comparing the accuracy of eDNA for taxonomic assignment and identification of nonindigenous taxa. We found that (1) different sampling access locations of the same tank resulted in different communities, (2) communities from treated and exchanged BW differ, (3) signals of source port and of ocean exchange are observed, (4) eDNA metabarcoding results in more diversity than morphological taxonomy, and (5) the nonindigenous copepod Oithona davisae, not reported before in the Chesapeake Bay, is detected. Overall, this study highlights the potential of eDNA metabarcoding for BW monitoring, but more comprehensive sampling will be needed to optimize the approach.

RAD-seq derived genome-wide nuclear markers resolve the phylogeny of tunas.

Although species from the genus Thunnus include some of the most commercially important and most severely overexploited fishes, the phylogeny of this genus is still unresolved, hampering evolutionary and traceability studies that could help improve conservation and management strategies for these species. Previous attempts based on mitochondrial and nuclear markers were unsuccessful in inferring a congruent and reliable phylogeny, probably due to mitochondrial introgression events and lack of enough phylogenetically informative markers. Here we infer the first genome-wide nuclear marker-based phylogeny of tunas using restriction site associated DNA sequencing (RAD-seq) data. Our results, derived from phylogenomic inferences obtained from 128 nucleotide matrices constructed using alternative data assembly procedures, support a single Thunnus evolutionary history that challenges previous assumptions based on morphological and molecular data.

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.

The RNA-binding protein human antigen R controls global changes in gene expression during Schwann cell development.

An important prerequisite to myelination in peripheral nerves is the establishment of one-to-one relationships between axons and Schwann cells. This patterning event depends on immature Schwann cell proliferation, apoptosis, and morphogenesis, which are governed by coordinated changes in gene expression. Here, we found that the RNA-binding protein human antigen R (HuR) was highly expressed in immature Schwann cells, where genome-wide identification of its target mRNAs in vivo in mouse sciatic nerves using ribonomics showed an enrichment of functionally related genes regulating these processes. HuR coordinately regulated expression of several genes to promote proliferation, apoptosis, and morphogenesis in rat Schwann cells, in response to NRG1, TGFß, and laminins, three major signals implicated in this patterning event. Strikingly, HuR also binds to several mRNAs encoding myelination-related proteins but, contrary to its typical function, negatively regulated their expression, likely to prevent ectopic myelination during development. These functions of HuR correlated with its abundance and subcellular localization, which were regulated by different signals in Schwann cells.

miRanalyzer: an update on the detection and analysis of microRNAs in high-throughput sequencing experiments.

We present a new version of miRanalyzer, a web server and stand-alone tool for the detection of known and prediction of new microRNAs in high-throughput sequencing experiments. The new version has been notably improved regarding speed, scope and available features. Alignments are now based on the ultrafast short-read aligner Bowtie (granting also colour space support, allowing mismatches and improving speed) and 31 genomes, including 6 plant genomes, can now be analysed (previous version contained only 7). Differences between plant and animal microRNAs have been taken into account for the prediction models and differential expression of both, known and predicted microRNAs, between two conditions can be calculated. Additionally, consensus sequences of predicted mature and precursor microRNAs can be obtained from multiple samples, which increases the reliability of the predicted microRNAs. Finally, a stand-alone version of the miRanalyzer that is based on a local and easily customized database is also available; this allows the user to have more control on certain parameters as well as to use specific data such as unpublished assemblies or other libraries that are not available in the web server. miRanalyzer is available at http://bioinfo2.ugr.es/miRanalyzer/miRanalyzer.php.

Comparison of methods to detect copy number alterations in cancer using simulated and real genotyping data.

BACKGROUND: The detection of genomic copy number alterations (CNA) in cancer based on SNP arrays requires methods that take into account tumour specific factors such as normal cell contamination and tumour heterogeneity. A number of tools have been recently developed but their performance needs yet to be thoroughly assessed. To this aim, a comprehensive model that integrates the factors of normal cell contamination and intra-tumour heterogeneity and that can be translated to synthetic data on which to perform benchmarks is indispensable. RESULTS: We propose such model and implement it in an R package called CnaGen to synthetically generate a wide range of alterations under different normal cell contamination levels. Six recently published methods for CNA and loss of heterozygosity (LOH) detection on tumour samples were assessed on this synthetic data and on a dilution series of a breast cancer cell-line: ASCAT, GAP, GenoCNA, GPHMM, MixHMM and OncoSNP. We report the recall rates in terms of normal cell contamination levels and alteration characteristics: length, copy number and LOH state, as well as the false discovery rate distribution for each copy number under different normal cell contamination levels.Assessed methods are in general better at detecting alterations with low copy number and under a little normal cell contamination levels. All methods except GPHMM, which failed to recognize the alteration pattern in the cell-line samples, provided similar results for the synthetic and cell-line sample sets. MixHMM and GenoCNA are the poorliest performing methods, while GAP generally performed better. This supports the viability of approaches other than the common hidden Markov model (HMM)-based. CONCLUSIONS: We devised and implemented a comprehensive model to generate data that simulate tumoural samples genotyped using SNP arrays. The validity of the model is supported by the similarity of the results obtained with synthetic and real data. Based on these results and on the software implementation of the methods, we recommend GAP for advanced users and GPHMM for a fully driven analysis.

Whole transcriptome analysis of a reversible neurodegenerative process in Drosophila reveals potential neuroprotective genes.

BACKGROUND: Neurodegenerative diseases are progressive and irreversible and they can be initiated by mutations in specific genes. Spalt-like genes (Sall) encode transcription factors expressed in the central nervous system. In humans, SALL mutations are associated with hereditary syndromes characterized by mental retardation, sensorineural deafness and motoneuron problems, among others. Drosophila sall mutants exhibit severe neurodegeneration of the central nervous system at embryonic stage 16, which surprisingly reverts later in development at embryonic stage 17, suggesting a potential to recover from neurodegeneration. We hypothesize that this recovery is mediated by a reorganization of the transcriptome counteracting SALL lost. To identify genes associated to neurodegeneration and neuroprotection, we used mRNA-Seq to compare the transcriptome of Drosophila sall mutant and wild type embryos from neurodegeneration and reversal stages. RESULTS: Neurodegeneration stage is associated with transcriptional changes in 220 genes, of which only 5% were already described as relevant for neurodegeneration. Genes related to the groups of Redox, Lifespan/Aging and Mitochondrial diseases are significantly represented at this stage. By contrast, neurodegeneration reversal stage is associated with significant changes in 480 genes, including 424 not previously associated with neuroprotection. Immune response and Salt stress are the most represented groups at this stage. CONCLUSIONS: We identify new genes associated to neurodegeneration and neuroprotection by using an mRNA-Seq approach. The strong homology between Drosophila and human genes raises the possibility to unveil novel genes involved in neurodegeneration and neuroprotection also in humans.

A high density SNP genotyping approach within the 19q13 chromosome region identifies an association of a CNOT3 polymorphism with ankylosing spondylitis.

OBJECTIVE: To identify genomic variants in the 19q13 chromosome region associated with ankylosing spondylitis (AS) in human leucocyte antigen (HLA)-B27-positive populations. METHODS: High-throughput genotyping of 1536 haplotype-tag single nucleotide polymorphisms (SNPs) was performed in 249 patients with AS and 302 healthy controls. Some of the identified associations were validated by genotyping four SNPs in two additional cohorts consisting of 412 cases/301 controls and 144 cases/203 controls. All individuals selected (both cases and controls) were HLA-B27-positive. RESULTS: Two markers in two different genes (CNOT3 and LAIR2) showed significant association (p<10(-3)) with AS. In addition, sliding windows analysis showed association of groups of adjacent SNPs in regions located around CNOT3 (Chr19: 59347459-59356564, p=2.43 × 10(-4) to 6.54 × 10(-4)). The associations were validated by genotyping four SNPs from regions located near LAIR2 and CNOT3 genes (rs1055234, rs8111398, rs2287828 and rs4591276) in two additional cohorts. The CNOT3 polymorphism (rs1055234) remained associated with AS (combined p=9.73 × 10(-6)). One SNP, located downstream of KIR3DL1, was detected which, tested in combination with HLA-Bw4I80, was associated with AS. CONCLUSION: A novel significant association was detected between SNP rs1055234 and AS susceptibility.

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.

Phylogenomic evidence for separate acquisition of plastids in cryptophytes, haptophytes, and stramenopiles.

According to the chromalveolate hypothesis (Cavalier-Smith T. 1999. Principles of protein and lipid targeting in secondary symbiogenesis: euglenoid, dinoflagellate, and sporozoan plastid origins and the eukaryote family tree. J Eukaryot Microbiol 46:347-366), the four eukaryotic groups with chlorophyll c-containing plastids originate from a single photosynthetic ancestor, which acquired its plastids by secondary endosymbiosis with a red alga. So far, molecular phylogenies have failed to either support or disprove this view. Here, we devise a phylogenomic falsification of the chromalveolate hypothesis that estimates signal strength across the three genomic compartments: If the four chlorophyll c-containing lineages indeed derive from a single photosynthetic ancestor, then similar amounts of plastid, mitochondrial, and nuclear sequences should allow to recover their monophyly. Our results refute this prediction, with statistical support levels too different to be explained by evolutionary rate variation, phylogenetic artifacts, or endosymbiotic gene transfer. Therefore, we reject the chromalveolate hypothesis as falsified in favor of more complex evolutionary scenarios involving multiple higher order eukaryote-eukaryote endosymbioses.

A review of the fisheries, life history and stock structure of tropical tuna (skipjack Katsuwonus pelamis, yellowfin Thunnus albacares and bigeye Thunnus obesus) in the Indian Ocean.

Skipjack (Katsuwonus pelamis), yellowfin (Thunnus albacares) and bigeye (Thunnus obesus) tuna are the target species of tropical tuna fisheries in the Indian Ocean, with high commercial value in the international market. High fishing pressure over the past three decades has raised concerns about their sustainability. Understanding life history strategies and stock structure is essential to determine species resilience and how they might respond to exploitation. Here we provide a comprehensive review of available knowledge on the biology, ecology, and stock structure of tropical tuna species in the Indian Ocean. We describe the characteristics of Indian Ocean tropical tuna fisheries and synthesize skipjack, yellowfin, and bigeye tuna key life history attributes such as biogeography, trophic ecology, growth, and reproductive biology. In addition, we evaluate the available literature about their stock structure using different approaches such as analysis of fisheries data, genetic markers, otolith microchemistry and tagging, among others. Based on this review, we conclude that there is a clear lack of ocean basin-scale studies on skipjack, yellowfin and bigeye tuna life history, and that regional stock structure studies indicate that the panmictic population assumption of these stocks should be investigated further. Finally, we identify specific knowledge gaps that should be addressed with priority to ensure a sustainable and effective management of these species.

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.