European minnows through time: museum collections aid genetic assessment of species introductions in freshwater fishes (Cyprinidae: Phoxinus species complex)

A chromosome-level, haplotype-resolved genome assembly and annotation for the Eurasian minnow (Leuciscidae: Phoxinus phoxinus) provide evidence of haplotype diversity

BACKGROUND

In this study, we present an in-depth analysis of the Eurasian minnow (Phoxinus phoxinus) genome, highlighting its genetic diversity, structural variations, and evolutionary adaptations. We generated an annotated haplotype-phased, chromosome-level genome assembly (2n = 50) by integrating high-fidelity (HiFi) long reads and chromosome conformation capture data (Hi-C).

RESULTS

We achieved a haploid size of 940 megabase pairs (Mbp) for haplome 1 and 929 Mbp for haplome 2 with high scaffold N50 values of 36.4 Mb and 36.6 Mb and BUSCO scores of 96.9% and 97.2%, respectively, indicating a highly complete genome assembly. We detected notable heterozygosity (1.43%) and a high repeat content (approximately 54%), primarily consisting of DNA transposons, which contribute to genome rearrangements and variations. We found substantial structural variations within the genome, including insertions, deletions, inversions, and translocations. These variations affect genes enriched in functions such as dephosphorylation, developmental pigmentation, phagocytosis, immunity, and stress response. In the annotation of protein-coding genes, 30,980 messenger RNAs and 23,497 protein-coding genes were identified with a high completeness score, which further underpins the high contiguity of our genome assemblies. We performed a gene family evolution analysis by comparing our proteome to 10 other teleost species, which identified immune system gene families that prioritize histone-based disease prevention over NB-LRR-related-based immune responses. Additionally, demographic analysis indicates historical fluctuations in the effective population size of P. phoxinus, likely correlating with past climatic changes.

CONCLUSIONS

This annotated, phased reference genome provides a crucial resource for resolving the taxonomic complexity within the genus Phoxinus and highlights the importance of haplotype-phased assemblies in understanding haplotype diversity in species characterized by high heterozygosity.

Unclear host taxonomy hinders parasite studies: An up-to-date checklist of the protozoan and metazoan parasites of Phoxinus minnows (Teleostei: Leuciscidae)

Phoxinus minnows are small freshwater fish found throughout Eurasia. The taxonomy of this genus is still under investigation, and new species are regularly described. They are frequently introduced outside their native range due to their use as live bait for angling, and such fish introductions can have diverse impacts on the recipient areas. Before the rise in interest regarding the diversity of this genus in the past 15 years, studies carried out on their biology and ecology reported their findings for the Eurasian minnow Phoxinus phoxinus Linnaeus, 1758, which was then considered ubiquitous. A parasites checklist for Phoxinus minnows has yet to be published, and the recent progress on the taxonomy of this genus has enabled us to propose a reassignment of parasite records to their proper host. The most thoroughly studied potential species were the North European species Phoxinus sp. 7 sensu Palandačić et al. 2017, Phoxinus isetensis, and P. phoxinus. We were able to highlight a gap in knowledge for numerous species that have received little-to-no attention in terms of parasitology. Systematic molecular identification of Phoxinus minnows should be conducted when studying their parasites, as a reliable identification of the host is vital in parasitology research. Future research will certainly reveal a species-specific composition of the parasitofauna occurring in Phoxinus, especially among widespread parasites such as Gyrodactylus and Diplostomum. Such specificities could even become tools for assessing the introduction routes of Phoxinus species.

Unraveling the molecular mechanisms of fish physiological response to freshwater salinization: A comparative multi-tissue transcriptomic study in a river polluted by potash mining

Freshwater salinization is an escalating global environmental issue that threatens freshwater biodiversity, including fish populations. This study aims to uncover the molecular basis of salinity physiological responses in a non-native minnow species (Phoxinus septimaniae x P. dragarum) exposed to saline effluents from potash mines in the Llobregat River, Barcelona, Spain. Employing high-throughput mRNA sequencing and differential gene expression analyses, brain, gills, and liver tissues collected from fish at two stations (upstream and downstream of saline effluent discharge) were examined. Salinization markedly influenced global gene expression profiles, with the brain exhibiting the most differentially expressed genes, emphasizing its unique sensitivity to salinity fluctuations. Pathway analyses revealed the expected enrichment of ion transport and osmoregulation pathways across all tissues. Furthermore, tissue-specific pathways associated with stress, reproduction, growth, immune responses, methylation, and neurological development were identified in the context of salinization. Rigorous validation of RNA-seq data through quantitative PCR (qPCR) underscored the robustness and consistency of our findings across platforms. This investigation unveils intricate molecular mechanisms steering salinity physiological response in non-native minnows confronting diverse environmental stressors. This comprehensive analysis sheds light on the underlying genetic and physiological mechanisms governing fish physiological response in salinity-stressed environments, offering essential knowledge for the conservation and management of freshwater ecosystems facing salinization.

An integrative approach highlights the discrepancy in the genetic, phenotypic, and presumptive taxonomic structure of Phoxinus (Actinopterygii, Leuciscidae, Phoxininae) in Bulgaria

The present drainage network of Bulgaria is the result of a complex Neogene and Quaternary evolution. Karst, which has developed on 23% of the territory, further complicates the hydrological pattern. Fresh waters of Bulgaria drain into the Black Sea and the Aegean Sea basins and can be roughly divided into the Danube (Middle and Lower Danube), non-Danube Black Sea, East Aegean, and West Aegean hydrological regions. Phoxinus, a small leuciscid fish, has a mosaic distribution in all four of these regions, inhabiting small mountainous and semi-mountainous streams. Based on morphology, it was identified as three species, Phoxinus phoxinus in the Danube, Phoxinus strandjae in the non-Danube, and Phoxinus strymonicus in West Aegean region. Later, molecular data revealed Phoxinus csikii and Phoxinus lumaireul in the Middle Danube and P. csikii in the Lower Danube. Phoxinus has been the focus of many studies, showing a high molecular and morphological diversity, which is not entirely consistent with previous morphology-only-based taxonomic concepts. In this study, molecular (a mitochondrial marker and a nuclear marker) and morphological data from both historical and recently sampled collections were analysed to assess the applicability of the integrative approach in Phoxinus. The results showed a significant influence of the complex paleo- and recent hydrology on the currently observed genetic structure of the considered populations and species. Furthermore, the study also demonstrated a strong influence of phenotypic plasticity on the morphological analysis of Phoxinus and the lack of a clear differentiation between P. csikii and P. strandjae. A barcoded specimen was designated as neotype to fix the species named P. strandjae in the current taxonomic concept. Finally, a significant discordance between genetically delimited clades and phenotypic groups did not allow a proper delineation of the species distributed in Bulgaria, demonstrating that more molecular markers are needed for further taxonomic study of the Phoxinus complex.

An annotated catalogue of selected historical type specimens, including genetic data, housed in the Natural History Museum Vienna

Museum collections are an important source for resolving taxonomic issues and species delimitation. Type specimens as name-bearing specimens, traditionally used in morphology-based taxonomy, are, due to the progress in historical DNA methodology, increasingly used in molecular taxonomic studies. Museum collections are subject to constant deterioration and major disasters. The digitisation of collections offers a partial solution to these problems and makes museum collections more accessible to the wider scientific community. The Extended Specimen Approach (ESA) is a method of digitisation that goes beyond the physical specimen to include the historical information stored in the collection. The collections of the Natural History Museum Vienna represent one of the largest non-university research centres in Europe and, due to their size and numerous type specimens, are frequently used for taxonomic studies by visiting and resident scientists. Recently, a version of ESA was presented in the common catalogue of the Fish and Evertebrata Varia collections and extended to include genetic information on type specimens in a case study of a torpedo ray. Here the case study was extended to a heterogeneous selection of historical type series from different collections with the type locality of Vienna. The goal was to apply the ESA, including genetic data on a selected set of type material: three parasitic worms, three myriapods, two insects, twelve fishes, and one bird species. Five hundred digital items (photographs, X-rays, scans) were produced, and genetic analysis was successful in eleven of the 21 type series. In one case a complete mitochondrial genome was assembled, and in another case ten short fragments (100-230 bp) of the cytochrome oxidase I gene were amplified and sequenced. For five type series, genetic analysis confirmed their taxonomic status as previously recognised synonyms, and for one the analysis supported its status as a distinct species. For two species, genetic information was provided for the first time. This catalogue thus demonstrates the usefulness of ESA in providing digitised data of types that can be easily made available to scientists worldwide for further study.

Phoxinus abanticus, a new species from the Lake Abant drainage in Türkiye (Teleostei: Leuciscidae)

Phoxinus abanticus, a new species, is described from the Lake Abant basin. It is distinguished from Phoxinus species in Türkiye and adjacent waters by the presence of fewer lateral line scales (60-69, vs. 75-91 in Phoxinus colchicus, 75-90 in Phoxinus strandjae); a deeper caudal peduncle (caudal peduncle depth: 1.8-2.3 times in length, vs. 2.4-2.9 in P. colchicus; 2.5-3.2 in P. strandjae); the absence of scales in the breast of males (vs. present); and ventral body reddish in nuptial colouration pattern for male (vs. brackish). The new species, P. abanticus, is also distinguished from its closest relative, P. strandjae, by a minimum of 3.40% genetic distance in the mtDNA cytochrome b (cyt b) gene.

A new set of microsatellite markers for Phoxinus lumaireul senso lato, Phoxinus marsilii and Phoxinus krkae for population and molecular taxonomic studies

Minnows of the genus Phoxinus are common and an often highly abundant fish species in Palearctic freshwater habitats. Phoxinus species have a complex evolutionary history, phylogenetic relationships are not well understood and there are a number of unresolved taxonomic problems. There are currently 23 different mitochondrial genetic lineages identified in the genus Phoxinus, 13 of which are recognized as valid species. The taxonomic status of these lineages requires resolution, including the degree to which they can interbreed. Suitable nuclear molecular markers for studies of population divergence and interbreeding between morphotypes and mitochondrial lineages are lacking for Phoxinus species. Therefore, the authors developed a set of microsatellite markers using genomic information from Phoxinus lumaireul and tested their suitability for this and two related species, Phoxinus krkae and Phoxinus marsilii. Out of 16 microsatellite candidate loci isolated, 12 were found to be in Hardy-Weinberg equilibrium when tested on two P. lumaireul senso lato populations. Seven loci amplified across the three species, enabling the study of intraspecific genetic diversity and population structure within P. marsilii and P. krkae. The markers were able to clearly resolve differences among the three tested species, including the recently described P. krkae, and are therefore suitable for the detection of introgression and hybridization among populations consisting of mixtures of two or more of P. lumaireul s. l., P. marsilii and P. krkae.

The minnow Phoxinus lumaireul (Leuciscidae) shifts the Adriatic-Black Sea basin divide in the north-western Dinaric Karst region

Karst landscapes are characterized by intermittent and sinking streams. The most common method used to study underground hydrological connections in karst is tracing tests. However, a more biologically oriented approach has been suggested: analysis of the genetic structure of aquatic organisms. Biological tracers can be sought among trogloxenes, that is, surface species that occasionally enter caves and groundwater. One such example is the fish genus Phoxinus, which exhibits high genetic diversity and complex phylogeography in the Balkan Peninsula. In the north-western Dinaric Karst, the complex hydrological network was digitalized in 2020. Contemporaneously, Phoxinus lumaireul populations in the Slovenian Dinaric Karst were intensively sampled and analysed for fragments of two mitochondrial genes and one nuclear gene. The derived phylogeographic structure and data on hydrological connections were compared to evaluate support for three alternative scenarios: The genetic structure (1) is a consequence of the ongoing geneflow through underground connections, (2) reflects a previous hydrological network or (3) is an outcome of anthropogenic translocations. The results suggest that the first two scenarios seem to have played a major role, while the third has not had profound effects on the genetic composition. Comparison between the genetic structure of Slovenian Dinaric Karst sampling sites and that of hydrologically isolated reference sampling sites indicated a greater genetic connectivity in the former. Moreover, the range of Adriatic (1a) and Black Sea (1c) haplotypes does not correspond to the Adriatic-Black Sea basin divide but is shifted northwards.

No strict host specificity: Brain metacercariae Diplostomum petromyzifluviatilis Müller (Diesing, 1850) are conspecific with Diplostomum sp. Lineage 4 of Blasco-Costa et al. (2014)

Metacercariae of Diplostomum petromyzifluviatilis (Digenea, Diplostomidae) from the brain of European river lamprey Lampetra fluviatilis from the Baltic Sea basin and Arctic lamprey Lethenteron camtschaticum from the White Sea basin were studied with the use of genetic and morphological methods. Phylogenetic analysis based on cox1 marker showed that the parasites of both lamprey species were conspecific with Diplostomum sp. Lineage 4 of Blasco-Costa et al. (2014). The name Diplostomum petromyzifluviatilis Müller (Diesing, 1850) has historical precedence as a species described from the brain of lampreys and should be used in genus nomenclature. There were no morphological qualitative differences between the metacercariae from the two lamprey species but those from L. fluviatilis were larger than those from L. camtschaticum. We expanded the data on the second intermediate hosts and the localization of D. petromyzifluviatilis, showing that its metacercariae occur not only in the brain of lampreys but also in the brain and the retina of three-spined stickleback Gasterosteus aculeatus and the vitreous humour of the perch Perca fluviatilis across the European part of the Palearctic.

A comprehensive DNA barcode inventory of Austria's fish species

Austria is inhabited by more than 80 species of native and non-native freshwater fishes. Despite considerable knowledge about Austrian fish species, the latest Red List of threatened species dates back 15 years and a systematic genetic inventory of Austria's fish species does not exist. To fulfill this deficit, we employed DNA barcoding to generate an up-to-date and comprehensive genetic reference database for Austrian fish species. In total, 639 newly generated cytochrome c oxidase subunit 1 (COI) sequences were added to the 377 existing records from the BOLD data base, to compile a near complete reference dataset. Standard sequence similarity analyses resulted in 83 distinct clusters almost perfectly reflecting the expected number of species in Austria. Mean intraspecific distances of 0.22% were significantly lower than distances to closest relatives, resulting in a pronounced barcoding gap and unique Barcode Index Numbers (BINs) for most of the species. Four cases of BIN sharing were detected, pointing to hybridization and/or recent divergence, whereas in Phoxinus spp., Gobio spp. and Barbatula barbatula intraspecific splits, multiple BINs and consequently cryptic diversity were observed. The overall high identification success and clear genetic separation of most of the species confirms the applicability and accuracy of genetic methods for bio-surveillance. Furthermore, the new DNA barcoding data pinpoints cases of taxonomic uncertainty, which need to be addressed in further detail, to more precisely assort genetic lineages and their local distribution ranges in a new National Red-List.

Reconnecting research and natural history museums in Italy and the need of a national collection biorepository

In Italy, differently from other countries, a national museum of natural history is not present. This absence is due, among other reasons, to its historical political fragmentation up to 1870, which led to the establishment of medium-sized museums, mostly managed by local administrations or universities. Moreover, a change of paradigm in biological research, at the beginning of the 20^th century, contributed to privilege experimental studies in universities and facilitated the dismissal of descriptive and exploratory biology, which formed the basis of the taxonomic research carried out by natural history museums. Consequently, only a few museums have a provision of curatorial staff, space and material resources adequate to maintain their original mission of discovering the natural world, by conducting a regular research activity accompanied by field campaigns. The creation of a national research centre for the study of biodiversity, facilitating interconnections among the existing natural history museums could be a solution and is here supported, together with a centralised biorepository to host collections and vouchers, to the benefit of current and future taxonomic research and environmental conservation. Such an institution should find place and realisation within the recently proposed National Biodiversity Future Center (NBFC) planned within the National Plan of Recovery and Resilience (PNRR). Pending upon the creation of this new national centre, a network among the existing museums should coordinate their activities.

Molecular analysis reveals multiple native and alien Phoxinus species (Leusciscidae) in the Netherlands and Belgium

Fresh waters are among the most endangered ecosystems, one of the problems being the lack of data on biodiversity. In the center of the missing knowledge are cryptic species, two (or more) species classified as a single one due to their (seemingly) indistinguishable morphology. Lack of research and stabilizing selection are reflected in the cryptic diversity of the genus Phoxinus (Leusciscidae), the studies of which have intensified over the last two decades and reveal undetected taxonomic complexity. Moreover, some of the Phoxinus lineages act as invasive species, while others are endangered by their alien counterparts. Minnows have been intentionally (as food for predatory fish species) or unintentionally (with other fries) stocked causing hybridisation zones in Norway, Portugal, Spain, France, Italy, Germany and Austria. Given that genetic identity and lineage assignment of Phoxinus from Belgium and the Netherlands have not been researched, the goal of the study was to examine available samples from known localities in the area in order to infer- whether they are native or not. For this purpose, the barcoding region cytochrome oxidase I, another mitochondrial gene cytochrome b, a nuclear recombination activating gene 1 and a combination of these markers from a wider neighboring region were analyzed. The study found four different Phoxinus species/lineages occurring in Belgium and the Netherlands: P. phoxinus, P. csikii, P. septimaniae and genetic lineage 11 (possibly P.cf. morella). While the first seem to be native, the other three were probably introduced.

Infirm effect of phylogeny on morphometric features in a cryptic Gobio species complex

Several recent notes prove that taxonomic relations of close relative animal groups (species complexes or cryptic species) can be revealed by the combined use of genetic and morphologic methodologies. At the same time scarce information can be found about how phylogeny, population origin, and sexual dimorphism affect the morphometric features of these species. In our present work, we performed simultaneous phylogenetic and morphological studies on the taxonomically still questionable Carpathian stream dwelling gudgeons (Cyprinidae, Gobio) by using two different methodologies (distance based and geometric morphometry). Our results were in correspondence with the previous findings, showing the presence of three phylogenetically more or less distinct groups in the area. The results of the whole-body geometric and the traditional, distance-based morphometry reflected the extent of phylogenetic differences. While the results of geometric scale morphometry did not correspond with the genetic subdivisions. Results of three way PERMANOVA analyses showed that the phylogenetic effects on morphometry is less considerable as the population origin or the sexual dimorphism at these cyprinid taxa. Our investigation contributed to the better understanding of the taxonomy of fish stocks in the Carpathian Basin, and to their conservation, but additional investigations will be needed to clarify the exact taxonomic position of the gudgeons (’Gobio sp1’) dominating the eastern part of the studied drainage.

Increased spatial resolution of sampling in the Carpathian basin helps to understand the phylogeny of central European stream-dwelling gudgeons

Background

Phylogenetic studies of widespread European fish species often do not completely cover their entire distribution area, and some areas are often excluded from analyses than others. For example, Carpathian stocks are often omitted from these surveys or are under-represented in the samples. However, this area served as an extra-Mediterranean refugia for many species; therefore, it is assumed that fish stocks here may show special phylogenetic features. For this reason, increased spatial resolution of sampling, namely revealing genetic information from unexamined Carpathian areas within the range of doubtful taxa, may help us better understand their phylogenetic features. To test this hypothesis, a phylogenetic investigation using a partial mtCR sequence data was conducted on 56 stream-dwelling freshwater fish (Gobio spp.) individuals collected from 11 rivers of the data-deficient Southeastern Carpathian area. Moreover, we revieved the available phylogenetic data of Middle-Danubian stream-dwelling gudgeon lineages to delineate their distribution in the area.

Results

Seven out of the nine detected haplotypes were newly described, suggesting the studied area hosts distinct and diverse Gobio stocks. Two valid species (G. obtusirostris, G. gobio), and a haplogroup with doubtful phylogenetic position” G. sp. 1" were detected in the area, showing a specific spatial distribution pattern. The distribution of the detected lineages in the Middle-Danubian area correspond with recent and paleo hydrogeographic features, at the same time mainly on their bordering areas show considerable overlap.

Conclusions

Despite the relatively limited geographic range of the study, our results provide important information which contributes to a better understanding of the phylogenetic, taxonomic and distribution features of Central European gudgeons. The genetically confirmed distribution data of the indicated lineages corresponds well with the recent and near-recent hydrogeographic features of the area, and assumes several hybrid zones in the Carpathian Basin. Additionally, the results show that the middle and lower Danubian watershed cannot be excluded from the range of G. gobio. Moreover, the” G. sp. 1", is slightly differentiated but phylogenetically distinct entity, and is the only Gobio taxa thus far detected in the middle and lower Tisza-basin. However, further investigations are necessary to clarify the taxonomic position of this group.

Adding DNA barcoding to stream monitoring protocols - What's the additional value and congruence between morphological and molecular identification approaches?

Although aquatic macroinvertebrates and freshwater fishes are important indicators for freshwater quality assessments, the morphological identification to species-level is often impossible and thus especially in many invertebrate taxa not mandatory during Water Framework Directive monitoring, a pragmatism that potentially leads to information loss. Here, we focus on the freshwater fauna of the River Sieg (Germany) to test congruence and additional value in taxa detection and taxonomic resolution of DNA barcoding vs. morphology-based identification in monitoring routines. Prior generated morphological identifications of juvenile fishes and aquatic macroinvertebrates were directly compared to species assignments using the identification engine of the Barcode of Life Data System. In 18% of the invertebrates morphology allowed only assignments to higher systematic entities, but DNA barcoding lead to species-level assignment. Dissimilarities between the two approaches occurred in 7% of the invertebrates and in 1% of the fishes. The 18 fish species were assigned to 20 molecular barcode index numbers, the 104 aquatic invertebrate taxa to 113 molecular entities. Although the cost-benefit analysis of both methods showed that DNA barcoding is still more expensive (5.30–8.60€ per sample) and time consuming (12.5h), the results emphasize the potential to increase taxonomic resolution and gain a more complete profile of biodiversity, especially in invertebrates. The provided reference DNA barcodes help building the foundation for metabarcoding approaches, which provide faster sample processing and more cost-efficient ecological status determination.

Species composition of introduced and natural minnow populations of the Phoxinus cryptic complex in the westernmost part of the Po River Basin (north Italy)

Invasive alien species are a major driver of biodiversity loss, with their impacts potentially more intense when complexes of cryptic species are involved. In freshwaters, the anthropogenic manipulation of fish communities has resulted in altered fish communities, and in Europe has increased the complexity of Phoxinus species assemblages. Here, we investigated the Phoxinus communities of the westernmost part of the Po river basin, where adjacent freshwater ecosystems (Alpine high-altitude lakes and lowland streams) are representative of different management strategies (i.e. manipulated fish communities via stocking in Alpine lakes vs. natural populations in streams). We tested the genetic composition of the cryptic Phoxinus populations inhabiting these waters, as the species are morphologically indistinct. Sequences of the mitochondrial cytochrome oxidase I (COI) were obtained from 239 specimens, with the results indicating that 17 Alpine high-altitude lakes are now populated by a complex of Phoxinus species, comprising P. septimaniae (native to the Mediterranean area of France), P. csikii (native to the Central Balkans) and P. lumaireul (native to the North Adriatic Sea basins). Their introduction resulted from their use as angling live baits. Minnow populations in lowland streams were primarily comprised of native P. lumaireul, with only a single P. csikii specimen detected. While nuclear sequences of the recombination activating gene 1 (RAG1) marker were not useful for tracking the presence of alien alleles in these stream populations, the COI data emphasised the importance of using molecular tools to investigate cryptic species complexes that have been modified by anthropogenic activities.