Skip to main content

Mitochondrial phylogeny and taxonomic revision of Italian and Slovenian fluvio-lacustrine barbels, Barbus sp. (Cypriniformes, Cyprinidae)

Abstract

Background

Barbels are ray finned cyprinid fishes of the Old-World with partially unresolved, intricate taxonomy. Within the Barbus sensu lato paraphyletic assemblage, Barbus sensu stricto is a monophyletic tetraploid lineage of Europe, northern Africa and Middle East, including two monophyletic sibling genera: Barbus and Luciobarbus. Italy, Slovenia and northern Croatia are natively inhabited by several entities of the genus Barbus, whose relationships and taxonomic ranks are still unclear. Aim of the present work is to focus on phylogeography of Italian and Slovenian barbels, with an appraisal of their current taxonomy.

Results

One hundred fifty specimens were collected in 78 sampling sites from 33 main watersheds, widely distributed along Italian and Slovenian ichthyogeographic districts. We amplified two mitochondrial markers, cytochrome b (cytb) and control region (D-loop), to infer a robust phylogeny for our sample and investigate on species delimitation.

Our results strongly indicate all Italian and Adriatic Slovenian fluvio-lacustrine barbels to be comprised into at least three distinct species. We provide a proposal of taxonomic revision and a list of synonymies for two of them and a new description under the International Code of Zoological Nomenclature rules for the third one.

Conclusions

If nuclear data will confirm our findings, at least three specific entities should be acknowledged across our sampling area. Namely, the three species are (i) Barbus plebejus, in the Padano-Venetian district; (ii) Barbus tyberinus, in the Tuscany-Latium district; (iii) Barbus oscensis Rossi & Plazzi sp. nov., in the Tyrrhenian and southernmost-Adriatic parts of Apulia-Campania district. Finally, we briefly discuss the implications of such a taxonomic scenario on conservation policies.

Background

The taxonomy of barbels (Actinopterygii: Cyprinidae) is a vexing question which remains partially unresolved. Barbus sensu lato is a paraphyletic taxon of the Old World, which has normally three levels of ploidy and is comprised by at least 800 species [1]. Conversely, Barbus sensu stricto (Barbus Cuvier and Cloquet 1816 sensu [2, 3]) is a monophyletic tetraploid genus of Europe, northern Africa and Middle East [4].

Doadrio [5] proposed further morphological subdivision of this taxon into two monophyletic sister subgenera: the nominotypical Barbus, and Luciobarbus Heckel 1843: the former is distributed in the northern part of the distribution range, whereas the latter is distributed in the southern watersheds.

Further molecular analyses confirmed the existence of these two clades within Barbus s. s [6,7,8,9,10,11,12,13].. They are currently recognised as distinct genera [14,15,16,17,18,19,20,21], and a third genus, Aulopyge Heckel, 1841, has been found to be closely related to the ancestor of both and has been included in the group [4, 9]. Aulopyge is a monotypic genus, with A. huegelii living in Croatia and Bosnia-Herzegovina.

The Italian peninsula is currently subdivided into three different districts or regions for what ichthyogeography is concerned. The Padano-Venetian district (PV [22, 23]), alias Padan region [24], corresponds to the Paleo-Po last glacial maximum catchment and includes Adriatic rivers of Slovenia and Croatia as well. Nearly all the remaining Italian mainland is subdivided into a northern Tuscany-Latium district (TL [22, 23]) and a southern Apulia-Campania district (AC [25]); territories (including also Italian islands and Corsica) with poorness or absence of primary freshwater fish species [22, 26] are not classified as districts or regions (Fig. 1).

Fig. 1
figure1

Approximate delimitation of ichthyogeographic districts in the area of interest (redrawn from districts in [22, 26, 27]) and localization of specimen collection. Abbreviations: PV, Padano-Venetian district; TL, Tuscany-Latium district; AC, Apulia-Campania district; DAN, Danubian district, NC, not classified as district. Map was generated by GR using the software QGIS 2.18 (https://qgis.org/it/site/index.html)

Native Italian barbels belong to genus Barbus [5, 14, 28] and are currently enlisted as four different species [15]: Barbus caninus Bonaparte 1839; Barbus balcanicus Kotlík et al. 2002; Barbus plebejus Bonaparte 1839; Barbus tyberinus Bonaparte 1839.

Barbus caninus and B. balcanicus are two rheophilic sibling vicariant species of the PV that also inhabit few TL and many not-Adriatic Dinaric/Balkan drainages, respectively [15, 18, 28,29,30,31]; although human manipulation cannot be excluded, their distribution may be linked to natural river capture phenomena.

As stated in [15, 28, 31, 32], the fluvio-lacustrine B. plebejus and B. tyberinus are vicariant sibling species too: the former native of PV [33], the latter naturally inhabiting TL and part of AC [34]. Other authors [35, 36] claimed these to be two morpho-species sensu Ruse [37], and pointed out that morphological differences between populations could be highly influenced by environmental factors, or caused by hybridization phenomena. Therefore, they do not recognize B. tyberinus as a valid species and include its populations and distribution range in those of B. plebejus.

Further insights did not settle the issue: focusing on intra- versus inter-group phenotypic similarity/dissimilarity, the existence of the two morpho-species was claimed by Lorenzoni et al. [38], but rejected by Livi et al. [39]. Nevertheless, even if the scenario proposed by Livi et al. [39] is correct, the two taxa may well represent geographically isolated, cryptic species sensu Bickford et al. [40].

Within genus Barbus introgressive hybridization has been repeatedly documented [18, 31, 41,42,43,44,45,46], occurring either in natural hybrid zones, (e.g., B. balcanicus × B. plebejus and B. caninus × B. plebejus [18]) or after a human-mediated secondary contact, e.g. between B. plebejus, B. tyberinus and B. barbus (L.) [18, 45, 47], the latter being a native species of the Danubian basin.

In fact, indigenous Italian barbel distribution is altered by anthropic activities [24], and invasive species introduction plays a key role as well [48], so that the existence of introgressive hybridization phenomena could hamper the analysis and disentanglement of morphological and genetic differentiation between the native Barbus species.

Despite that the taxonomic ranking debate is still unresolved, the use of the name B. tyberinus increased over time by leading organizations as the International Union for Conservation of Nature [49] – Comitato Italiano IUCN [50], the Institute for Biodiversity Science and Sustainability [51], and the Species 2000 and ITIS Integrated Taxonomic Information System [52]. Unstable nomenclature could cause wildlife management aberrations, with population of B. tyberinus possibly losing the conservation status they had as a part of B. plebejus (species listed in Annexes II and V of the European Union Habitats Directive and in Appendix III of the Bern Convention).

Recently, Buonerba et al. [18] analysed mitochondrial and nuclear markers in Barbus specimens from PV and TL and demonstrated the genetic distinguishability and close relatedness of B. tyberinus and B. plebejus. The genetic distinguishability of B. tyberinus and B. plebejus was also observed by [53, 54] using fragments of the D-loop control region and cytochrome b gene. Interestingly, they also identified two new, distinct and allopatric Barbus clades. The former was only found in Eastern Central Italy, in the Adriatic watersheds of Vomano [53], Aterno-Pescara, Sangro and Biferno [54]. The latter is distributed in Adriatic watersheds to the south of Biferno river (Fortore and Ofanto) and in southernmost investigated Tyrrhenian watersheds (Liri-Garigliano, Volturno and Sele) [54].

The existence of genetically distinct allopatric taxa at least identifies them as Evolutionary Significant Units [55], but their taxonomic rank remains undefined. Nonetheless, phylogeographic structures and genetic similarity/dissimilarity between and within territories could be a diagnostic key for the definition of ichthyogeographic (or more generally biogeographic) districts [56,57,58]. Therefore, AC can be split into its northern Adriatic part (NAAC) and its Tyrrhenian and southern Adriatic part (TSAAC).

Specifically, in the cited works, nine sampling sites from AC were analysed covering four out of the five main watersheds of NAAC (lacking Trigno river basin), but only five out of the 14 main watersheds of TSAAC (missing, north to south, the Adriatic Candelaro, Carapelle, Cervaro and the Tyrrhenian Savone, Sarno, Tusciano, Alento, Mingardo, Bussento river basins) in which Barbus species are possibly autochthonous according to Bianco [28]. Moreover, two mitochondrial markers were used, but the sequencing of the fragment of the cytochrome b gene was limited to a subsample of 26 specimens showing 26 different D-loop haplotypes [54].

In this work, we analysed phylogenetic relationships within Italian barbels increasing the number of analysed localities, primarily in TSAAC and in the Slovenian part of PV. We present a thorough phylogenetic reconstruction of Italian barbels using two mitochondrial markers for each specimen instead of collapsing to haplotypes. This allowed to assess genetic and morphological differences in terms of similarity/dissimilarity between and within clades. We focused on the hypothesis of the existence of two new species previously undescribed under the rules of the International Code of Zoological Nomenclature (ICZN). In this scenario, the current nomenclature in use for Italian fluvio-lacustrine barbels would be obsolete, since the two putative new species are until now enlisted as part of B. tyberinus.

We then found supporting evidence for the description of a new species corresponding to the TSAAC clade, for which we propose the name B. oscensis, though we advise caution until present data are confirmed by thorough molecular studies using also nuclear markers. Finally, we also review range, description and synonymy for the already described B. plebejus and B. tyberinus and give suggestions for conservation policies of Italian barbels.

Methods

Specimen collection, morphologic analysis, identification and preparation

After acquiring the relevant approvals and permissions to collect animals, we collected 138 specimens in 72 sampling sites from 32 main watersheds (i.e., basins of rivers emptying into a sea), widely distributed along PV, TL and AC districts (29, 31 and 12 sampling sites, respectively; Fig. 1), to which 12 specimens in 6 sampling sites from the Danubian district have been added (DAN; Fig. 1). The distribution of sampling sites per district, main watershed and river is provided in Additional file 1. No animals were excluded from subsequent analyses.

Specimens were caught using electrofishing devices with direct current and electrical settings that minimize possible stresses to the animals [59], made unconscious with the anaesthetic 2-phenoxyethanol 0,5 ml/L, and rapidly identified utilizing meristic and qualitative external characters (not to harm any specimen) as in the dichotomous key of Kottelat and Freyhof [15] for barbels of Apennine Peninsula and Adriatic basin of Slovenia and Croatia, provided as Additional file 2. Fishes were then biopsied regardless of morphological determination and released alive and conscious in the place where caught. Biopsies (clips of caudal fin small enough to avoid a significant effect on specimens’ motility) were stored in 100% ethanol and refrigerated at 4 °C.

Six specimens from a single target sampling site belonging to the new clade described in this work were collected, humanely euthanized and fixed in 10% formalin for museum conservation and conforming to Directive 2010/63/EU (provided that all the permissions requested under the Italian law had been granted). Morphological description of the six specimens was done on the basis of morphometrics measures (eye diameter; preorbital distance; mouth-operculum distance; length of pectoral fin; length of ventral fin; length of anal fin; height of the third dorsal fin ossified ray) and meristic counts (branched rays in the dorsal, anal and pelvic fins; scales on the lateral line, above the lateral line, below the lateral line; number of circumpeduncolar scales; gill rakers in the lower arch and in the upper arch; total gill rakers; pharyngeal teeth in the left and in the right side; serrae on ossified ray of fin), following previous works on this species complex [28, 33, 54, 60].

PCR amplification and sequencing

DNA was extracted from caudal fin clips stored in 100% ethanol using a standard phenol:chloroform protocol [61]. PCR amplification of portions of cytochrome b gene (cytb) and tRNA-Pro/control region (trnP/D-loop) was carried out with GoTaq® Flexi DNA Polymerase (Promega, Madison, WI, USA), as follows: 10 μL 5× Green GoTaq® Flexi Buffer, MgCl2 (3 mM), ddNTPs (800 μM each), primers (500 nM each), 1 U GoTaq® DNA Polymerase, 40 ng template DNA, ddH2O up to 50 μL. Primers were designed for the present study and are BARBUS 8 (5′-GCGCTAGGGAGGAGTTTA-3′) and BARBUS 5 (5’TTTTAACCGAGACCAATGAC-3′) for cytb and dloop sxF (D1) (5′-AAAGCATCGGTCTTGTAATC-3′) and dloop dxR (D2) (5′-GAGTTTTCTAGGACCCATCTTA-3′) for trnP/D-loop. Annealing settings were 60 °C/35″ and 55 °C–58 °C/35″, respectively. PCR results were visualized using a 1% electrophoresis agarose gel stained with ethidium bromide. Amplicons were purified and sequenced at the Macrogen Europe facility (Amsterdam, The Netherlands). Electropherograms were edited using MEGA 7 [62].

Phylogenetic analyses

This paper presents the first phylogenetic reconstruction of Italian barbels using two mitochondrial markers for each specimen. To our knowledge, all the previous papers focusing on Italian barbels (e.g., [18, 45, 53, 54, 63]) presented haplotypes instead of sequences of single individuals. As a consequence, it was not possible to include all the available literature in tree reconstruction, with the only exception of [54], who published the association of 19 D-loop haplotypes with the respective cytb haplotype: these sequences were added to our dataset. Cyprinus carpio (L.) and Luciobarbus graellsii (Steindachner, 1866) sequences were downloaded from GenBank (Accession Numbers DQ868875/JN105352 and JN049525/MG827110, respectively) and used as outgroups.

The T-Coffee [64] algorithm was used for single alignments. Nucleotide sequences of cytb and D-loop were separately aligned through the M-Coffee approach, starting from MAFFT [65] and Muscle [66] libraries. Sites with low or noisy phylogenetic signal were masked using Gblocks 0.91b [67]. The cytb alignment was further subdivided into the three codon positions using a custom-tailored Python script (available from FP upon request), obtaining five datasets: cytb1_1, cytb1_2, cytb1_3, trnP, and dloop.

The datasets were concatenated into the final dataset and PartitionFinder 1.1.0 [68] was used to decide whether to apply a single-partition or a multiple-partition scheme, as well as to select molecular evolution models; the Bayesian Information Criterion and the greedy approach were chosen. Three methods were then selected to reconstruct the phylogenetic tree of Italian barbels.

First, the software RaxML 8.2.11 [69] was used for the Maximum Likelihood (ML) inference under the CAT model. The Best-Known Likelihood Tree (BKLT) was computed and then it was annotated with bootstrap support values and using 1000 bootstrap replicates. Furthermore, Bayesian Inference was also carried out with MrBayes 3.2.7 [70] using two separate runs, four chains, and 10,000,000 generations of MC3, sampling every 100 trees. Convergence between runs and burn-in were estimated looking to four diagnostics: standard deviation of average split frequencies sampled every 1000 generation, Potential Scale Reduction Factor (PSRF [71]), plot of log-likelihoods on the MCMC generations, and minimum Estimated Sample Sizes. Finally, the phylogenetic inference was carried out using IQ-TREE 1.7-beta7 [72] with 1000 ultrafast bootstrap replicates [73]. In the IQ-TREE analysis, substitution models were selected using ModelFinder [74] and the best partitioning scheme was selected with the greedy strategy implemented in ModelFinder [68, 75].

Taxonomic unit definition and networks

Two different barcode gap (see, f.i., [76, 77]) approaches were explored to test taxonomic unit definition among the native Italian barbel clades (B. plebejus, B. tyberinus, Barbus sp. clade 4, and Barbus sp. TSAAC clade). First, we applied the Automatic Barcode Gap Discovery (ABGD [78]), using the K-2-P distance and a transition/transversion ratio set to 2. Prior intraspecific distance ranged from 0.001 and 0.01, with 20 steps. Moreover, we used the Species Delimitation Plugin of Geneious Prime® 2021.0.3 (Biomatters Ltd., Auckland, New Zealand) to compute the Rosenberg’s PAB [79] and the mean probabilities of making a correct identification of an unknown specimen of the focal taxon following [80].

Given the short length of the amplified trnP fragment (76 bp), the minimum spanning network was computed from multiple cytb and D-loop sequence alignments only. Newly produced sequences were added to previous data taken from [18, 45, 53, 54, 63]; the software PopART v 1.7 [81] was used to draw the network, using the statistical parsimony criterion and setting ε = 0. Data from literature were assigned the relative geographical abundancies upon retrieval in the original publication. Unfortunately, the geographical information associated to some Tyrrhenian watershed haplotypes by [54] was missing or misspelt (see Table 1 in [54]), therefore we could not unambiguously assign them to either TSAAC or NAAC district. All sequences used for single-marker minimum spanning networks are listed in Additional file 3. Conversely, data from the present work refer to single specimens, therefore an abundancy of 1 was assigned to the relative ichthyogeographical district (see Additional file 3). Finally, for the combined cytb-dloop network, we retained only those 155 entries for which an association between cytb and dloop was available; the analysis was carried out as above.

Morphological analysis

Data registered for the formal description of the new clade types were added to those of the corresponding SI2 lineage of [54] and compared to other Italian clades from the same work (B barbus, B. plebejus, B. tyberinus, and the SI1 lineage, corresponding to the NAAC lineage). Single measurements were not published in [54], therefore it was not possible to standardize morphometric data using the Beacham [82] formula to reduce the effects of size and allometry. Therefore, statistical analyses (ANOVA and Tukey post-hoc test) were performed exclusively on meristic traits published by [54]. Aggregate data (average value, standard deviation and sample size [54];) were merged with new data and a statistical analysis was carried out using custom scripts available from FP upon request.

Meristic data used for the formal description of the types of the new clade were compared to the values of B. plebejus and B. tyberinus published in the taxonomic revision of Bianco [28, 33, 34], which includes the respective type specimens. Since maximum ranges and usual vales (generically defined) are only reported therein, no quantitative measure of data dispersion was available; hence, no statistical test was performed and only a graphical comparison was carried out.

Nomenclatural acts

The electronic edition of this article conforms to the requirements of the amended International Code of Zoological Nomenclature, and hence the new name contained herein is available under that Code from the electronic edition of this article. The electronic edition of this work was published in a journal with an ISSN, and has been archived and is available from the following digital repositories: PubMed Central, LOCKSS.

Results

Phylogenetic analysis

The final dataset is composed by 171 sequences, including the only 19 sequences that were possible to obtain from [54] and two outgroups: C. carpio and L. graellsii. However, three specimens (GR633Omb11BA, GR570Usi200BA, and GR568Usi200BA) from our sample turned out to be L. graellsii as well. Indeed, this Iberian species is enlisted within Italian non-native barbels [17, 83]: thus, it is unsurprising to sporadically catch L. graellsii specimens among native Italian barbels.

The final cytb alignment is 1141 bp long, while the final trnP/D-loop alignment is 620 bp long (including 76 and 544 bp of trnP and D-loop, respectively). Both markers were amplified from all individuals: new sequences were uploaded to GenBank under Accession Numbers MG495623-MG495922. Specimens vouchers, taxonomy following our phylogenetic reconstruction, locality (country, river, altitude) and accession numbers are provided as Additional file 4.

The strategy selected by PartitionFinder was to keep all markers together in a single partition and the best-fitting molecular evolutionary model was TrN + I + G. Conversely, the ModelFinder algorithm of IQ-TREE selected a partitioning scheme with four partitions: cytb_1 + trnP, cytb_2, cytb_3, and dloop. The best-fitting molecular evolutionary models were K-2-P + G4, TPM3u + F + I, TIM2 + F, and HKY + F + R3, respectively.

The three trees yielded identical topologies, the only differences being in the node support values (Fig. 2): some internal nodes are not statistically supported in some trees, but topologies are never mutually exclusive (Fig. 2 and Additional file 5).

Fig. 2
figure2

Consensus phylogenetic tree of barbels of Italy and Adriatic catchments of Slovenia and Croatia, based on combined cytb/trnP/D-loop data. Cyprinus carpio and Luciobarbus graellsii were used as outgroups. Node support values are shown in the following order: RAxML boostrap proportion (BP)/Bayesian posterior probability (PP)/IQ-TREE ultrafast bootstrap (UF_Boot) support values. For sake of clarity, only internal nodes are shown: nodes with BP < 60, PP < 0.950, and UF-Boot support value< 85 were collapsed. For complete trees, explicit GenBank Accession Numbers, and branch lengths, we refer to Additional file 5. Photograph by AM

Italian fluvio-lacustrine barbels comprise a single monophyletic clade, which is highly supported, with a bootstrap proportion (BP) of 98, a posterior probability (PP) of 1.000 and an ultrafast bootstrap (UF-Boot) support value of 95. Within this clade, four main lineages were statistically evidenced. The first lineage, corresponding to B. tyberinus, almost exclusively encompasses samples from TL district (BP = 100; PP = 1.000, UF-Boot support value = 100). The only exception is represented by the B. tyberinus AQ84Lir330BA specimen caught in the Liri river, which is part of the northernmost watershed (Liri-Garigliano) of the TSAAC district, at the boundary with TL.

The second lineage, corresponding to B. plebejus (BP = 76; PP = 1.000; UF-Boot support value = 88.0), is mainly comprised by samples from the PV district, but few samples were also caught in TL and AC (see Additional files 4 and 5).

A third lineage (BP = 65; PP = 1.000, UF-Boot support value = 97.0) consists exclusively of samples from the TSAAC district and five out nineteen samples from [54] are nested within this clade. A possible fourth lineage (PP = 1.000; UF-Boot support value = 98) was only found in the NAAC district and is comprised by the remaining samples from [54]; this lineage was not supported as monophyletic in the RAxML analysis. Phylogenetic relationships between these four clades are not well resolved: the Bayesian inference suggest the first clade to be the sister group of a monophyletic clade comprised by the remaining three (PP = 0.954; UF-Boot support value = 61).

Taxonomic unit definition and networks

The number of groups identified in the ABGD analysis ranged from 7 to 2 moving from a prior intraspecific divergence (P) equal to 0.001/0.003 through P = 0.01 (Fig. 3). In the former case (low Ps), we retrieved two groups with a single specimen (either RA606Sin150BA or a specimen from [54] corresponding to Accession Numbers MK728798/MK728816), one group with two specimens (PU245Mta595BA and PU246Mta595BA), and four groups corresponding to the four clades of native Italian barbels depicted in Fig. 2. In the latter case (high Ps), a group with a single specimen (RA606Sin150BA, as before) and a group with all remaining specimens were recovered. At an intermediate P = 0.0033, this large assembly was split into PV specimens, one side, and all remaining specimens, the other side. The four clades were consistently recovered as reciprocally monophyletic across the three tree-building methods, with Rosenberg’s PAB always smaller than 3.10 × 10− 11 (Additional file 6); moreover, the mean probabilities of making a correct identification of an unknown specimen of the focal taxon were always greater than 0.90, with the only exception of the NAAC district lineage, that showed slightly smaller values (Additional file 6).

Fig. 3
figure3

Automatic Barcode Gap Discovery. The number of identified groups is plotted over the prior intraspecific divergence used to identify the barcode gap. Major monophyletic groups, corresponding to clades in Fig. 2, are shown close to the respective points; minor groups, comprised either by a single OTU or by two OTUs, are detailed in the text. If possible, the method recursively splits (blue circles) the original partition (green circles), but in most case they overlap (red circles)

The minimum spanning network computed by the concatenated cytb-dloop alignment (155 sequences/1761 bp) identifies the four clades as above (Fig. 4). The cluster of Danubian haplotypes is connected to the B. tyberinus TL clade and then to the B. plebejus PV clade. From PV haplotypes, haplotypes found in TSAAC and NAAC districts are separated by few mutation steps; moreover, haplotypes of ambiguous geographical origin (see below) connect these Southern Italy lineages. However, single-marker networks including a larger amount of available data (205 sequences/1141 bp and 230 sequences/547 bp for cytb and dloop, respectively) yielded slightly different results: Danubian cytb haplotypes are connected to TSAAC haplotypes and then to PV ones, while Danubian D-loop haplotypes are connected to TL haplotypes and then to PV ones (Additional file 7). PV haplotypes are always connected to TL haplotypes, one side, and NAAC ones, the other side (Fig. 4; Additional file 7).

Fig. 4
figure4

Minimum spanning network computed on the combined cytb/D-loop dataset, using published abundance data from [54]. Samples that was not possible to unambiguously assign to an ichthyogeographical district are shaded in red and brown

Within each clade, the mean cytb uncorrected p-distance is one order of magnitude lower than between clades; the four Italian clades show between-groups mean uncorrected p-distance values ranging from 1.55 to 2.22%, while these values rise to 3.86–4.25%, 8.31–8.96%, and 8.38–9.19% when the four Italian clades are compared with B. barbus, B. balcanicus, and B. caninus lineages, respectively (Table 1; see also Additional file 6). Similar results were retrieved from the D-loop fragment, which, as expected, is in some cases more variable: in some clades the intra-group mean uncorrected p-distance is greater than 1% (B. balcanicus and B. caninus). However, the four Italian clades show between-groups mean uncorrected p-distance values ranging from 0.83 to 2.64%, while these values rise to 4.86–5.74%, 7.37–8.07%, and 5.81–6.02% when the four Italian clades are compared with B. barbus, B. balcanicus, and B. caninus lineages, respectively (Table 1).

Table 1 Within- and between-groups uncorrected p-distance (mean ± standard deviation) among the seven Barbus lineages shown in Fig. 2

Morphological analysis

Morphological character values of type specimens of the TSAAC clade are reported in Additional file 8 together with meristic character values of B. plebejus and B. tyberinus in the taxonomic revision of Bianco [28, 33, 34], which includes type specimens.

Graphical morphological comparison between the three clades shows moderate to high overlapping values: the most reliable characters, number of scales on lateral line and around caudal peduncle, only permit discrimination of B. plebejus. The other two taxa remain indistinguishable. Graphical comparisons are provided as Additional file 9.

Meristic character values used to test morphological differences are reported in Table 2 together with ANOVA results. Meristic traits of the type specimens of the Barbus sp. TSAAC lineage here described were merged with those of the corresponding SI2 linage of [54] and tested against the other clades recorded by the same authors (B. barbus, B. plebejus, B. tyberinus, and Barbus sp. NAAC lineage). In all cases, results from ANOVA were highly significant (p < 0.001), meaning that a structure is present in the provided groups.

Table 2 Meristic traits of the four Italian fluvio-lacustrine barbel clades along with Barbus barbus. Data from 6 specimens of the Barbus sp. TSAAC lineage presented in this paper have been added to data published by [54]. Results from the Analysis of Variance (ANOVA) are also shown. The table lists mean ± standard deviation; observed ranges are reported in parentheses

Tukey post-hoc test for pairwise comparisons (provided as Additional file 10), confirmed most differences between the analysed groups as significant, hence evidencing morphological differentiation between them.

Concluding, given the molecular and morphological results presented above, we propose the Barbus sp. TSAAC clade to be a new species: Barbus oscensis Rossi & Plazzi sp. nov.

Discussion

Phylogenetics and systematics of barbels in continental Italy

Our analysis confirms the monophyly and genetic distinguishability of at least three out of the four taxa previously identified [18, 53, 54]: a Padano-Venetian clade, a Tuscany-Latium clade, and a new clade distributed in Tyrrhenian watersheds and southernmost-Adriatic part of Apulia-Campania district, deepening genetic structure knowledge, taxonomic position and geographic distribution of the latter. The reciprocal monophyly of the two established clades and the new clade is repeatedly confirmed in our analysis, and the three groups were not identified only for relatively high values of prior intraspecific divergence, nor were B. tyberinus and B. plebejus (Fig. 3). In a nutshell, this means that, if B. tyberinus and B. plebejus are distinguished as monophyletic entities, the same must hold for the two AC Barbus sp. clades.

Conversely, the fourth lineage identified by [54] is distributed in the northernmost-Adriatic part of Apulia-Campania district, but further insights are required to assess its identity. In fact, some of the original coordinates [54] are lacking or result in unclear spatial reference system, and therefore geographic details assumed as valid here were taken from text and cartography of the cited work, but there is no clear indication of sampling sites. Moreover, in our phylogenetic reconstruction, this clade was not always recovered as monophyletic (Fig. 2; see also Additional file 6) and its relationships with other lineages are not unambiguously resolved and appear to be highly dependent from the dataset (see Figs. 2, 3 and 4; Additional file 7).

As is the case for previous work, our study is based on mitochondrial markers that are commonly employed in freshwater fish systematics. Therefore, we advise caution while interpreting present results, until they are confirmed by future studies using nuclear markers. A thorough genetic approach may shed light on hybridization and introgression phenomena as well, which may be easily overlooked using mitochondrial datasets. Alternative explanations of our results, such as retentions of ancestral polymorphisms and introgressions, were repeatedly observed in Barbinae subfamily [20, 41, 42, 83,84,85,86,87,88].

Nonetheless, our results are supported by available morphological data and our analyses do suggest that the three identified taxa are new species with mutually exclusive distribution. In fact, the two known exceptions can be explained either by human-mediated or natural events.

  1. (i)

    Samples of Padano-Venetian lineage found outside the PV district show a patchwork pattern distribution and have no original haplotype. Moreover, Tuscany-Latium and Apulia-Campania district were historically subject to transplantation of native species from PV district [89]. This well-documented phenomenon (also known as ichthyofaunal “padanization”, from the Padan Plain) is due to governmental institutional restocking plans: public and private ichthyogenic centers located in the PV district were used as main sources of specimens starting from the end of the nineteenth century [23]. The opposite phenomenon, the transplantation of native species from Tuscany-Latium or Apulia-Campania districts to other districts, is not documented.

  2. (ii)

    Three haplotypes of the TL lineage were found in the Liri-Garigliano basin in the present work, while [54] individuated two haplotypes of Tyrrhenian and southern Adriatic parts of Apulia-Campania species. The Liri-Garigliano basin is the northernmost watershed of the Apulia-Campania district at the boundary with the Tuscany-Latium district. The Gari river, which is now part of the Liri-Garigliano basin, was probably flowing in the southern Volturno basin until the growing of Roccamonfina volcano, which took place approximately 630,000 years ago [90]. This may explain the presence of haplotypes from different lineages in this river.

Indeed, each of the three Italian fluvio-lacustrine barbels species recognized here do not trespass the extension of a single Italian mainland ichthyogeographic district sensu [25].

Taxonomy of barbels in continental Italy

The Padano-Venetian lineage fully corresponds to the species B. plebejus, whose first valid description (based on type specimens from the PV district only) is the one of Bonaparte in 1839 [91]. A revision of the synonymies for this species is provided as Appendix 1.

The Tuscany-Latium lineage exactly corresponds to the species originally described as B. fluviatilis tyberinus by Bonaparte in 1839 [91] (based on type specimens sampled in the TL district only), which is an older synonym of B. fucini Costa 1853; this species was previously considered [28, 34] inhabiting the Apulia-Campania district as well. A revision of the synonymies for this species is provided as Appendix 2.

Finally, the Tyrrhenian and southern Adriatic parts of Apulia-Campania lineage has no valid description as a species under the rules of the International Code of Zoological Nomenclature (ICZN). We then propose the name B. oscensis sp. nov. and provide the original description, diagnosis, and pictures hereafter (Fig. 5).

Fig. 5
figure5

Holotype of Barbus oscensis; original pictures. a Lateral view with registration details of the Museum of Zoology of Bologna – MZB. b Outer gill rakers count of the first (anterior-most) gill arch on the right side. c Pharyngeal teeth count. d Dorsal fin rays, with undistinguished serration of the last unbranched ray. e Scales below lateral line. f Scales above lateral line. g Circumpeduncolar scales. h Anal fin rays. i Pelvic fin rays. R, right; L, left; ε, lower limb of the gill arch; δ, upper limb of the gill arch; ₺, distinctly formed rudiment; α, inner row of pharyngeal teeth; β, middle row of pharyngeal teeth; γ, outer row of pharyngeal teeth; ₸, fallen or broken tooth; ζ, last unbranched ray

Original description of Barbus oscensis Rossi, G. & Plazzi, F. sp. nov.

Family Cyprinidae Rafinesque 1815

Subfamily Barbinae Bleeker 1859

Genus Barbus Cuvier and Cloquet 1816

Barbus oscensis Rossi & Plazzi sp. nov. (Fig. 5)

Pro parte synonymy

Barbus fucini Costa 1838 sensu [23]: 463 (partim: TSAAC); Barbus fluviatilis plebejus (non Bonaparte 1839) [92] (partim: TSAAC); Barbus plebejus (non Bonaparte 1839) [93]: 72–77 (partim: TSAAC) [94]: 11 (partim: TSAAC) [95]: 43 (partim: Volturno River near Venafro); [96]: 42–43 (partim: TSAAC); [97]: 40–41 (partim: TSAAC); [98]: 274–275 (partim: TSAAC); [35]: 198–202 (partim: TSAAC); Barbus Barbus plebejus (non Bonaparte 1839) [99]: 172 (partim: TSAAC); [100]: 19 (partim: TSAAC); [101]: 274–275 (partim: TSAAC); Barbus tyberinus (non Bonaparte 1839) [28]: 313–318, figs. 4b, 6–7 (Partim: TSAAC); [102]: 50–51 (Partim: TSAAC); [14]: 234 (Partim: TSAAC); [34]: 427 (Partim: TSAAC).

Type specimens

All specimens were sampled with electrofishing from a single sampling site and conserved in alcohol 70% after formalin fixation in Museum of Zoology of Bologna (MZB).

Holotype

MZB 201001 (GenBank Accession Numbers MG495912 and MG495762 for cytb and D-loop, respectively).

Paratypes

MZB 201002–201,006 (GenBank Accession Numbers MG495913, MG495914, MG495915, MG495917, MG495918 for cytb and MG495763, MG495764, MG495765, MG495767, MG495768 for D-loop).

Type locality

Fiumara di Atella stream (River Ofanto watershed; Adriatic side of Southern Apennine), 340 m above sea level in Basilicata Region at the foot of Vulture Mountain in the vicinity of “il Calvario”.

Distribution

TSAAC

Etymology

Specific name oscensis is derived from Osci: the name of an ancient Italian people inhabiting a territory strongly overlapping the range of the taxon during the Iron Age.

Diagnosis

Tyrrhenian and southern Adriatic parts of Apulia-Campania lineage shows a genetic differentiation from other fluvio-lacustrine barbel groups statistically supported in all the phylogenetic analyses hereby conducted. Moreover, the amount of genetic difference of this clade with the well-established species B. plebejus and B. tyberinus is comparable with the genetic distance between them.

Even from a morphological perspective this clade is statistically distinguishable from other fluvio-lacustrine barbel groups. However, these differences, albeit relevant to support a phylogenetic differentiation, are not reliable enough in terms of field discrimination of the different groups. This is even more true when considering that many of those differences rely (see also [54]) on morphometric traits that can only be identified after image analysis and statistic tests. In this sense, the TSAAC lineage could be considered an almost cryptic species that recently speciated from other fluvio-lacustrine species. Being a new species previously undescribed under International Code of Zoological Nomenclature rules, although yet observed as a proper taxonomic entity in other works [54], we provide a full description under the proposed name of Barbus oscensis sp. nov.

The holotype and paratypes of the species are conserved in the Museum of Zoology of Bologna (MZB).

Morphology

Count and measurement are given in Additional file 8; last unbranched ray of the dorsal fin poorly ossified and with undistinguished serration; serration could be distinguishable in younger specimens as in the other two fluvio-lacustrine subspecies; ossification reduced in the distal portion of the ray.

Remarks

B. oscensis sp. nov is a fluvio-lacustrine barbel species with TSAAC range, parapatric to B. plebejus (PV range) and B. tyberinus (TL range). Conversely, on the basis of morphological analysis [28, 34], B. oscensis was enlisted undistinguished in B. tyberinus species. The TSAAC range identified by the presence of this species include the Tyrrhenian basins from Liri-Garigliano to Bussento and the Adriatic basins from Fortore to Ofanto.

Conservation and management of Italian barbels

The phylogenetic subdivision of B. plebejus into at least three Barbus species supports the partitioning of continental Italy in at least the three districts proposed by Bianco and de Filippo [25] on the basis of district-specific endemisms and subendemisms. In fact, each of the three allopatric fluvio-lacustrine barbels species hereby identified is a district-specific endemism too and an Evolutionary Significant Unit (ESU) by itself.

Results achieved by this study relies on a wide distribution of the sampling sites that, differently from previous works, covers most watersheds of Italian mainland and Slovenia. However, it must be noted that further investigation is needed in order to fully characterize the areal of the newly described species, B. oscensis, and verify the possibly the existence of a fourth species.

Moreover, since our results evidence a mutually exclusive range distribution of fluvio-lacustrine barbels species with morphological character moderately to strongly overlapping, it cannot be neglected that morphologically indistinguishable unknown endemisms may inhabit the few watersheds not investigated here. This is partially confirmed by the distinct lineage previously observed [53, 54], outside the range of our sampling design, whose sequences mostly comprise our Barbus sp. clade 4 (Fig. 2). To date, we regard at the northernmost-Adriatic part of the Apulia-Campania district as a currently unidentified district, possibly inhabited by a further undescribed species. Therefore, further investigations are needed for barbel populations of Abruzzo and Molise Regions, as well as of the Ionic side of the Basilicata Region.

Conclusions and final remarks

The up-to-date taxonomic revision proposed in the present paper has clear effects on management and conservation policies. Historically, freshwater fish repopulations (either for sport fishing or conservation policies) were managed taking into account administrative borders and specific ranks rather than ichthyogeographic districts and ESUs. As fluvio-lacustrine cyprinids, the three species revised in the present study occupy similar habitats in their respective allopatric ranges; they can therefore be considered as vicariants species. Thus, the biology of Italian barbels and the lack of Italian laws dealing with transplantations within national borders led to indiscriminate repopulation and genetic erosion of wildlife Barbus populations – and the situation may get still worse.

When the Habitat Directive of the European Union (Council Directive 92/43/EEC) was written, all these species were considered as B. plebejus; therefore, every Italian fluvio-lacustrine barbel species should inherit the status of B. plebejus and deserves the proper effort in terms on conservation. In fact, each species must be protected as a separate ESU, avoiding transplantation at least between ichthyogeographic districts. On a precautionary principle, transplantation should be avoided even between different main watersheds, so that ESUs from currently unidentified districts are also protected.

Furthermore, since the taxonomy used in wildlife conservation laws (e.g., the Habitat Directive of the European Union or the Bern Convention) cannot be easily updated following scientific evidence, every wildlife conservation law addressed to any of the fluvio-lacustrine Italian barbels should be extended to any (albeit pro parte) possible taxonomic synonym here reported (see Discussion and Appendices 1 and 2). Concluding, since taxonomic revisions increased exponentially after the development of molecular analyses, wildlife legislation must also target ESUs and Management Units (beside species) in order to reach long-term effectiveness.

Availability of data and materials

The datasets generated during and analysed during the current study are available in the NCBI GenBank repository, under Accession Numbers MG495623-MG495922. Specimens vouchers, taxonomy following our phylogenetic reconstruction, locality (country, river, altitude) and accession numbers are provided as S5. A custom R script was developed to separate intra- and inter-clade distances and is freely available for download from GitHub (https://github.com/federicoplazzi/inter-intra). All B. oscensis type specimens are conserved in alcohol 70% after formalin fixation in Museum of Zoology of Bologna (MZB).

Abbreviations

PV:

Padano-Venetian district

TL:

Tuscany-Latium district

AC:

Apulia-Campania district

NAAC:

Northern Adriatic part of Apulia-Campania district

TSAAC:

Tyrrhenian and southern Adriatic part of Apulia-Campania district

DAN:

Danubian district

cytb :

Cytochrome b

trnP :

tRNA-Pro

BP:

Bootstrap proportion

PP:

Posterior probability

UF-Boot:

Ultrafast bootstrap support value

P:

Prior intraspecific distance

References

  1. 1.

    Chenuil A, Galtier N, Berrebi P. A test of the hypothesis of an autopolyploid vs. allopolyploid origin for a tetraploid lineage: application to the genus Barbus (Cyprinidae). Heredity. 1999;82(4):373–80. https://0-doi-org.brum.beds.ac.uk/10.1038/sj.hdy.6884890.

    Article  PubMed  Google Scholar 

  2. 2.

    Howes GJ. The phylogenetic position of the Yugoslavian cyprinid fish genus Aulopyge Heckel, 1841: with an appraisal of the genus Barbus Cuvier and Cloquet, 1816 and the subfamily Cyprininae. Bulletin of the British Museum of Natural History. Zool Series. 1987;52:165–96.

    Google Scholar 

  3. 3.

    Berrebi P, Kottelat M, Skelton P, Rab P. Systematics of Barbus: state of the art and heuristic comments. Folia Zool. 1996;45(Supplement 1):5–12.

    Google Scholar 

  4. 4.

    Berrebi P, Tsigenopoulos CS. Phylogenetic organization of the genus Barbus sensu stricto: a review based on data obtained using molecular markers. In: Banarescu P, Bogutskaya NG, editors. The freshwater fishes of Europe (vol. 5/II: Cyprinidae 2, Part II: Barbus). Wiebelsheim: Aula-Verlag; 2003. p. 11–22.

    Google Scholar 

  5. 5.

    Doadrio I. Phylogenetic relationships and classification of western palaearctic species of the genus Barbus (Osteichthyes, Cyprinidae). Aquat Living Resour. 1990;3(4):265–82. https://0-doi-org.brum.beds.ac.uk/10.1051/alr:1990028.

    Article  Google Scholar 

  6. 6.

    Zardoya R, Doadrio I. Phylogenetic relationships of Iberian cyprinids: systematic and biogeographical implications. Proc R Soc Lond Ser B Biol Sci. 1998;265(1403):1365–72. https://0-doi-org.brum.beds.ac.uk/10.1098/rspb.1998.0443.

    CAS  Article  Google Scholar 

  7. 7.

    Zardoya R, Doadrio I. Molecular evidence on the evolutionary and biogeographical patterns of European cyprinids. J Mol Evol. 1999;49(2):227–37. https://0-doi-org.brum.beds.ac.uk/10.1007/PL00006545.

    CAS  Article  PubMed  Google Scholar 

  8. 8.

    Zardoya R, Economidis PS, Doadrio I. Phylogenetic relationships of Greek Cyprinidae: molecular evidence for at least two origins of the Greek cyprinid fauna. Mol Phylogenet Evol. 1999;13(1):122–31. https://0-doi-org.brum.beds.ac.uk/10.1006/mpev.1999.0630.

    CAS  Article  PubMed  Google Scholar 

  9. 9.

    Tsigenopoulos CS, Berrebi P. Molecular phylogeny of North Mediterranean freshwater barbs (genus Barbus: Cyprinidae) inferred from cytochrome b sequences: biogeographic and systematic implications. Mol Phylogenet Evol. 2000;14(2):165–79. https://0-doi-org.brum.beds.ac.uk/10.1006/mpev.1999.0702.

    CAS  Article  PubMed  Google Scholar 

  10. 10.

    Machordom A, Doadrio I. Evidence of a Cenozoic Betic–Kabilian connection based on freshwater fish phylogeography (Luciobarbus, Cyprinidae). Mol Phylogenet Evol. 2001a;18(2):252–362. https://0-doi-org.brum.beds.ac.uk/10.1006/mpev.2000.0876.

    CAS  Article  PubMed  Google Scholar 

  11. 11.

    Machordom A, Doadrio I. Evolutionary history and speciation modes in the cyprinid genus Barbus. Proc R Soc Lond Ser B Biol Sci. 2001b;268(1473):1297–306. https://0-doi-org.brum.beds.ac.uk/10.1098/rspb.2001.1654.

    CAS  Article  Google Scholar 

  12. 12.

    Durand JD, Tsigenopoulos CS, Ünlü E, Berrebi P. Phylogeny and biogeography of the family Cyprinidae in the Middle East inferred from cytochrome b DNA—evolutionary significance of this region. Mol Phylogenet Evol. 2002;22(1):91–100. https://0-doi-org.brum.beds.ac.uk/10.1006/mpev.2001.1040.

    CAS  Article  PubMed  Google Scholar 

  13. 13.

    Tsigenopoulos CS, Durand JD, Ünlü E, Berrebi P. Rapid radiation of the Mediterranean Luciobarbus species (Cyprinidae) after the Messinian salinity crisis of the Mediterranean Sea, inferred from mitochondrial phylogenetic analysis. Biol J Linn Soc. 2003;80(2):207–22. https://0-doi-org.brum.beds.ac.uk/10.1046/j.1095-8312.2003.00237.x.

    Article  Google Scholar 

  14. 14.

    Bianco PG. Diversity of Barbinae fishes in southern Europe with description of a new genus and a new species (Cyprinidae). Ital J Zool. 1998;65(Supplement):125–36.

    Article  Google Scholar 

  15. 15.

    Kottelat M, Freyhof J. Handbook of European freshwater fishes. Cornol: Kottelat – Berlin: Freyhof; 2007.

    Google Scholar 

  16. 16.

    Clavero M, Hermoso V. Reservoirs promote the taxonomic homogenization of fish communities within river basins. Biodivers Conserv. 2011;20(1):41–57. https://0-doi-org.brum.beds.ac.uk/10.1007/s10531-010-9945-3.

    Article  Google Scholar 

  17. 17.

    Buonerba L, Pompei L, Lorenzoni M. First record of Iberian barbel Luciobarbus graellsii (Steindachner, 1866) in the Tiber River (Central Italy). BioInvasions Rec. 2013;2(4):297–301. https://0-doi-org.brum.beds.ac.uk/10.3391/bir.2013.2.4.06.

    Article  Google Scholar 

  18. 18.

    Buonerba L, Zaccara S, Delmastro GB, Lorenzoni M, Salzburger W, Gante HF. Intrinsic and extrinsic factors act at different spatial and temporal scales to shape population structure, distribution and speciation in Italian Barbus (Osteichthyes: Cyprinidae). Mol Phylogenet Evol. 2015;89:115–29. https://0-doi-org.brum.beds.ac.uk/10.1016/j.ympev.2015.03.024.

    Article  PubMed  Google Scholar 

  19. 19.

    Berrebi P, Chenuil A, Kotlík P, Machordom A, Tsigenopoulos CS. Disentangling the evolutionary history of the genus Barbus sensu lato, a twenty years adventure. In: Alves MJ, Cartaxana A, Correia AM, Lopes LF, editors. Professor Carlos Almaça (1934–2010) – Estado da Arte em Áreas Cientificas que Desenvolveu. Lisboa: Museu Nacional de Historia Natural e da Ciencia; 2014. p. 29–55.

    Google Scholar 

  20. 20.

    Gante HF, Doadrio I, Alves MJ, Dowling TE. Semi-permeable species boundaries in Iberian barbels (Barbus and Luciobarbus, Cyprinidae). BMC Evol Biol. 2015;15:1.

    Article  Google Scholar 

  21. 21.

    Casal-López M, Pere S, Yahyaoui A, Doadrio I. Taxonomic review of the genus Luciobarbus Heckel, 1843 (Actinopterygii, Cyprinidae) from northwestern Morocco with the description of three new species. Graellsia. 2015;71:e027.

    Article  Google Scholar 

  22. 22.

    Bianco PG. L’inquadramento zoogeografico dei pesci d’acqua dolce d’Italia e problemi determinati dalle falsificazioni faunistiche in Biologia e gestione dell’ittiofauna autoctona. In: AIIAD, editors. Atti del II Convegno nazionale dell’Associazione Italiana Ittiologi Acque Dolci AIIAD. Regione Piemonte: AIIAD; 1987. pp. 41–65.

  23. 23.

    Bianco PG. L’ittiofauna continentale dell’Appennino umbro-marchigiano, barriera semipermeabile allo scambio di componenti primarie tra gli opposti versanti dell’Italia centrale. Biogeographia. 1994;17:427–85.

    Google Scholar 

  24. 24.

    Zerunian S. Condannati all’estinzione? Biodiversità, biologia, minacce e strategie di conservazione dei Pesci d’acqua dolce indigeni in Italia. Bologna: Edagricole; 2002.

    Google Scholar 

  25. 25.

    Bianco PG, de Filippo G. (eds). Contributo alla conoscenza della fauna ittica d’acqua dolce in aree protette d’Italia. Res Wildl Conserv. 3. USA: IGF Publications; 2011.

  26. 26.

    Bianco PG. Potential role of the palaeohistory of the Mediterranean and Paratethys basins on the early dispersal of Euro-Mediterranean freshwater fishes. Ichthyol Explor Freshw. 1990;1:167–84.

    Google Scholar 

  27. 27.

    Bianco PG. An update on the status of native and exotic freshwater fishes of Italy. J Appl Ichthyol. 2014;30(1):62–77. https://doi.org/10.1111/jai.12291.

    Article  Google Scholar 

  28. 28.

    Bianco PG. A revision of the Italian Barbus species (Cypriniformes: Cyprinidae). Ichthyol Explor Freshw. 1995;6:305–24.

    Google Scholar 

  29. 29.

    Kotlík P, Berrebi P. Genetic subdivision and biogeography of the Danubian rheophilic barb Barbus petenyi inferred from phylogenetic analysis of mitochondrial DNA variation. Mol Phylogenet Evol. 2002;24(1):10–8. https://0-doi-org.brum.beds.ac.uk/10.1016/S1055-7903(02)00264-6.

    Article  PubMed  Google Scholar 

  30. 30.

    Kotlík P, Tsigenopoulos CS, Ráb P, Berrebi P. Two new Barbus species from the Danube River basin, with redescription of B. petenyi (Teleostei: Cyprinidae). Folia Zool. 2002;51:227–40.

    Google Scholar 

  31. 31.

    Tsigenopoulos CS, Kotlík P, Berrebi P. Biogeography and pattern of gene flow among Barbus species (Teleostei: Cyprinidae) inhabiting the Italian Peninsula and neighbouring Adriatic drainages as revealed by allozyme and mitochondrial sequence data. Biol J Linn Soc. 2002;75:3–99.

    Article  Google Scholar 

  32. 32.

    Ketmaier V, Finamore F, Largiadèr C, Milone M, Bianco PG. Phylogeography of bleaks Alburnus spp. (Cyprinidae) in Italy, based on cytochrome b data. J Fish Biol. 2009;75(5):997–1017. https://0-doi-org.brum.beds.ac.uk/10.1111/j.1095-8649.2009.02357.x.

    CAS  Article  PubMed  Google Scholar 

  33. 33.

    Bianco PG. Barbus plebejus Bonaparte, 1839. In: Banarescu P, Bogutskaya NG, editors. The freshwater fishes of Europe (vol. 5/II: Cyprinidae 2, Part II: Barbus). Wiebelsheim: Aula-Verlag; 2003a. p. 339–64.

    Google Scholar 

  34. 34.

    Bianco PG. Barbus tyberinus Bonaparte, 1839. In: Banarescu P, Bogutskaya NG, editors. The freshwater fishes of Europe (vol. 5/II: Cyprinidae 2, Part II: Barbus). Wiebelsheim: Aula-Verlag; 2003b. p. 427–49.

    Google Scholar 

  35. 35.

    Gandolfi G, Zerunian S, Torricelli P, Marconato A. I Pesci delle acque interne italiane. Roma: Istituto Poligrafico e Zecca dello Stato; 1991.

    Google Scholar 

  36. 36.

    Zerunian S. Pesci delle acque interne d’Italia. Quaderni Conservazione Natura, vol. 20. Roma: Ministero dell’Ambiente e Istituto Nazionale Fauna Selvatica; 2004.

    Google Scholar 

  37. 37.

    Ruse M. Definitions of species in biology. Br J Philos Sci. 1969;20(2):97–119. https://0-doi-org.brum.beds.ac.uk/10.1093/bjps/20.2.97.

    Article  Google Scholar 

  38. 38.

    Lorenzoni M, Carosi A, Angeli V, Bicchi A, Pedicillo G, Viali P. Individuazione e riconoscimento dei barbi autoctoni nel bacino del fiume Paglia. Terni: Arti Grafiche Iezzi; 2006.

    Google Scholar 

  39. 39.

    Livi S, de Innocentiis S, Longobardi A, Cataudella S, Tancioni L, Rampacci M, et al. Genetic structure of Barbus spp. populations in the Marches Region of central Italy and its relevance to conservation actions. J Fish Biol. 2013;82(3):806–26. https://0-doi-org.brum.beds.ac.uk/10.1111/jfb.12021.

    CAS  Article  PubMed  Google Scholar 

  40. 40.

    Bickford D, Lohman DJ, Sodhi NS, Ng PKL, Meier R, Winker K, et al. Cryptic species as a window on diversity and conservation. Trends Ecol Evol. 2007;22(3):148–55. https://0-doi-org.brum.beds.ac.uk/10.1016/j.tree.2006.11.004.

    Article  Google Scholar 

  41. 41.

    Machordom A, Berrebi P, Doadrio I. Spanish barbel hybridization detected using enzymatic markers: Barbus meridionalis Risso × Barbus haasi Mertens (Osteichthyes, Cyprinidae). Aquat Living Resour. 1990;3(4):295–303. https://0-doi-org.brum.beds.ac.uk/10.1051/alr:1990030.

    Article  Google Scholar 

  42. 42.

    Berrebi P, Cattaneo-Berrebi G, Le Brun N. Natural hybridization of two species of tetraploid barbels: Barbus meridionalis and Barbus barbus (Osteichtyes, Cyprinidae) in southern France. Biol J Linn Soc. 1993;48(4):319–33. https://0-doi-org.brum.beds.ac.uk/10.1111/j.1095-8312.1993.tb02094.x.

    Article  Google Scholar 

  43. 43.

    Chenuil A, Crespin L, Pouyaud L, Berrebi P. Autosomal differences between males and females in hybrid zones: a first report from Barbus barbus and Barbus meridionalis (Cyprinidae). Heredity. 2004;93(2):128–34. https://0-doi-org.brum.beds.ac.uk/10.1038/sj.hdy.6800455.

    CAS  Article  PubMed  Google Scholar 

  44. 44.

    Lajbner Z, Šlechtová V, Šlechta V, Švátora M, Berrebi P, Kotlík P. Rare and asymmetrical hybridization of the endemic Barbus carpathicus with its widespread congener B. barbus. J Fish Biol. 2009;74(2):418–36. https://0-doi-org.brum.beds.ac.uk/10.1111/j.1095-8649.2008.02098.x.

    CAS  Article  PubMed  Google Scholar 

  45. 45.

    Meraner A, Venturi A, Ficetola GF, Rossi S, Candiotto A, Gandolfi A. Massive invasion of exotic Barbus barbus and introgressive hybridization with endemic Barbus plebejus in northern Italy: where, how and why? Mol Ecol. 2013;22(21):5295–312. https://0-doi-org.brum.beds.ac.uk/10.1111/mec.12470.

    CAS  Article  PubMed  Google Scholar 

  46. 46.

    Levin BA, Gandlin AA, Simonov ES, Levina MA, Barmintseva AE, Japoshvili B, et al. Phylogeny, phylogeography and hybridization of Caucasian barbels of the genus Barbus (Actinopterygii, Cyprinidae). Mol Phylogenet Evol. 2019;135:31–44. https://0-doi-org.brum.beds.ac.uk/10.1016/j.ympev.2019.02.025.

    Article  PubMed  Google Scholar 

  47. 47.

    Geiger MF, Schreiner C, Delmastro GB, Herder F. Combining geometric morphometrics with molecular genetics to investigate a putative hybrid complex: a case study with barbels Barbus spp. (Teleostei: Cyprinidae). J Fish Biol. 2016;88(3):1038–55. https://0-doi-org.brum.beds.ac.uk/10.1111/jfb.12871.

    CAS  Article  PubMed  Google Scholar 

  48. 48.

    Zerunian S. Problematiche di conservazione dei Pesci d’acqua dolce italiani. Biologia Ambientale. 2007;21:49–55.

    Google Scholar 

  49. 49.

    IUCN. The IUCN red list of threatened species. Version 2016-1: International Union for Conservation of Nature and Natural Resources; 2019. Available from: https://www.iucnredlist.org

    Google Scholar 

  50. 50.

    Bianco PG, Caputo V, Ferrito V, Lorenzoni M, Nonnis Marzano F, Stefani F, et al. Barbus tyberinus. In: Rondinini C, Battistoni A, Peronace V, Teofili C, editors. Lista Rossa IUCN dei Vertebrati Italiani. Roma: Comitato Italiano IUCN e Ministero dell’Ambiente e della Tutela del Territorio e del Mare; 2013.

    Google Scholar 

  51. 51.

    Eschmeyer WN, Fricke R, van der Laan R. Catalog of fishes: genera, species, references: California Academy of Sciences; 2019. Available from: http://researcharchive.calacademy.org/research/ichthyology/catalog/fishcatmain.asp

    Google Scholar 

  52. 52.

    Froese R, Pauly D. FishBase (version Jan 2016). In: Roskov Y, Abucay L, Orrell T, Nicolson D, Flann C, Bailly N, et al., editors. Species 2000 & ITIS catalogue of life, 2016 annual checklist. Leiden: Naturalis; 2016.

    Google Scholar 

  53. 53.

    Zaccara S, Quadroni S, Vanetti I, Carosi A, La Porta G, Crosa G, et al. Morphologic and genetic variability in the Barbus fishes (Teleostei, Cyprinidae) of Central Italy. Zool Scr. 2019a;48(3):289–301. https://0-doi-org.brum.beds.ac.uk/10.1111/zsc.12341.

    Article  Google Scholar 

  54. 54.

    Zaccara S, Quadroni S, De Santis V, Vanetti I, Carosi A, Britton R, et al. Genetic and morphological analyses reveal a complex biogeographic pattern in the endemic barbel populations of the southern Italian peninsula. Ecol Evol. 2019b;9(18):10185–97. https://0-doi-org.brum.beds.ac.uk/10.1002/ece3.5521.

    Article  PubMed  PubMed Central  Google Scholar 

  55. 55.

    Moritz C. Defining evolutionarily significant units for conservation. Trends Ecol Evol. 1994;9:373–4.

    CAS  PubMed  Article  Google Scholar 

  56. 56.

    Zaccara S, Delmastro GB. Tyrrhenian basins of Ligury as a new peri-Mediterranean ichthyogeographic district? Population structure of Telestes muticellus (Osteichthyes, Cyprinidae), a primary freshwater fish. Hydrobiologia. 2009;632(1):285–95. https://0-doi-org.brum.beds.ac.uk/10.1007/s10750-009-9851-4.

    CAS  Article  Google Scholar 

  57. 57.

    Marchetto F, Zaccara S, Muenzel FM, Salzburger W. Phylogeography of the Italian vairone (Telestes muticellus, Bonaparte 1837) inferred by microsatellite markers: evolutionary history of a freshwater fish species with a restricted and fragmented distribution. BMC Evol Biol. 2010;10(1):111. https://0-doi-org.brum.beds.ac.uk/10.1186/1471-2148-10-111.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  58. 58.

    Ketmaier V, Bianco PG. Understanding and conserving genetic diversity in a world dominated by alien introductions and native transfers: the case study of primary and peripheral freshwater fishes in southern Europe. In: Closs GP, Krkosek M, Olden JD, editors. Conservation of freshwater fishes. Cambridge: Cambridge University Press; 2015. p. 506–34.

    Google Scholar 

  59. 59.

    Snyder DE. Electrofishing and its harmful effects on fish. Information and Technology Report USGS/BRD/ITR--2003–0002. Denver: U. S. Government Printing Office; 2003.

    Google Scholar 

  60. 60.

    Holčík J, Banarescu P, Evans D. General introduction to fishes. In: Holčík J, editor. The freshwater fishes of Europe. Vol. 1, Part II. Wiesbaden: Aula Verlag; 1989. p. 18–147.

    Google Scholar 

  61. 61.

    Sambrook J, Russell DW. Purification of nucleic acids by extraction with phenol:chloroform. CSH Protoc. 2006;2006(1):pdb.prot4455. https://0-doi-org.brum.beds.ac.uk/10.1101/pdb.prot4455.

    Article  PubMed  Google Scholar 

  62. 62.

    Kumar S, Stecher G, Tamura K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33(7):1870–4. https://0-doi-org.brum.beds.ac.uk/10.1093/molbev/msw054.

    CAS  Article  PubMed  Google Scholar 

  63. 63.

    Zaccara S, Antognazza CM, Buonerba L, Britton R, Crosa G. Human-mediated contact zone between endemic and invasive Barbus species (Osteichthyes: Cyprinidae) in a regulated lowland river: genetic inferences and conservation implications. Ital J Zool. 2014;81(4):571–83. https://0-doi-org.brum.beds.ac.uk/10.1080/11250003.2014.944225.

    Article  Google Scholar 

  64. 64.

    Notredame C, Higgins DG, Heringa J. T-Coffee: a novel method for fast and accurate multiple sequence alignment. J Mol Biol. 2000;302(1):205–17. https://0-doi-org.brum.beds.ac.uk/10.1006/jmbi.2000.4042.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  65. 65.

    Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30(4):772–80. https://0-doi-org.brum.beds.ac.uk/10.1093/molbev/mst010.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  66. 66.

    Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32(5):1792–7. https://0-doi-org.brum.beds.ac.uk/10.1093/nar/gkh340.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  67. 67.

    Castresana J. Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol Biol Evol. 2000;17(4):540–52. https://0-doi-org.brum.beds.ac.uk/10.1093/oxfordjournals.molbev.a026334.

    CAS  Article  PubMed  Google Scholar 

  68. 68.

    Lanfear R, Calcott B, Ho SYW, Guindon S. PartitionFinder: combined selection of partitioning schemes and substitution models for phylogenetic analyses. Mol Biol Evol. 2012;29(6):1695–701. https://0-doi-org.brum.beds.ac.uk/10.1093/molbev/mss020.

    CAS  Article  Google Scholar 

  69. 69.

    Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014;30(9):1312–3. https://0-doi-org.brum.beds.ac.uk/10.1093/bioinformatics/btu033.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  70. 70.

    Ronquist F, Teslenko M, van der Mark P, Ayres DL, Darling A, Höhna S, et al. MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst Biol. 2012;61(3):539–42. https://0-doi-org.brum.beds.ac.uk/10.1093/sysbio/sys029.

    Article  PubMed  PubMed Central  Google Scholar 

  71. 71.

    Gelman A, Rubin DB. Inference from iterative simulation using multiple sequences. Stat Sci. 1992;7:457–511.

    Google Scholar 

  72. 72.

    Nguyen L-T, Schmidt HA, von Haeseler A, Minh BQ. IQ-TREE: a fast and effective stochastic algorithm for estimating maximum likelihood phylogenies. Mol Biol Evol. 2015;32(1):268–74. https://0-doi-org.brum.beds.ac.uk/10.1093/molbev/msu300.

    CAS  Article  Google Scholar 

  73. 73.

    Hoang DT, Chernomor O, von Haeseler A, Minh BQ, Vinh LS. UFBoot2: improving the ultrafast bootstrap approximation. Mol Biol Evol. 2018;35(2):518–22. https://0-doi-org.brum.beds.ac.uk/10.1093/molbev/msx281.

    CAS  Article  Google Scholar 

  74. 74.

    Kalyaanamoorthy S, Minh BQ, Wong TKF, von Haeseler A, Jermiin LS. ModelFinder: fast model selection for accurate phylogenetic estimates. Nat Methods. 2017;14(6):587–9. https://0-doi-org.brum.beds.ac.uk/10.1038/nmeth.4285.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  75. 75.

    Chernomor O, von Haeseler A, Minh BQ. Terrace aware data structure for phylogenomic inference from supermatrices. Syst Biol. 2016;65(6):997–1008. https://0-doi-org.brum.beds.ac.uk/10.1093/sysbio/syw037.

    Article  PubMed  PubMed Central  Google Scholar 

  76. 76.

    Hebert PDN, Cywinska A, Ball SL, Dewaard JR. Biological identifications through DNA barcodes. Proc R Soc Lond Ser B Biol Sci. 2003;270(1512):313–21. https://0-doi-org.brum.beds.ac.uk/10.1098/rspb.2002.2218.

    CAS  Article  Google Scholar 

  77. 77.

    Kim S, Eo H-S, Koo H, Choi J-K, Kim W. DNA barcode-based molecular identification system for fish species. Mol Cells. 2010;30(6):507–12. https://0-doi-org.brum.beds.ac.uk/10.1007/s10059-010-0148-2.

    CAS  Article  PubMed  Google Scholar 

  78. 78.

    Puillandre N, Lambert A, Brouillet S, Achaz G. ABGD, automatic barcode gap discovery for primary species delimitation. Mol Ecol. 2012;21(8):1864–77. https://0-doi-org.brum.beds.ac.uk/10.1111/j.1365-294X.2011.05239.x.

    CAS  Article  Google Scholar 

  79. 79.

    Rosenberg NA. Statistical tests for taxonomic distinctiveness from observations of monophyly. Evolution. 2007;61(2):317–23. https://0-doi-org.brum.beds.ac.uk/10.1111/j.1558-5646.2007.00023.x.

    Article  PubMed  Google Scholar 

  80. 80.

    Ross HA, Murugan S, Li WLS. Testing the reliability of genetic methods of species identification via simulation. Syst Biol. 2008;57(2):216–30. https://0-doi-org.brum.beds.ac.uk/10.1080/10635150802032990.

    Article  PubMed  Google Scholar 

  81. 81.

    Leigh JW, Bryant D. PopArt: full-feature software for haplotype network construction. Methods Ecol Evol. 2015;6(9):1110–6. https://0-doi-org.brum.beds.ac.uk/10.1111/2041-210X.12410.

    Article  Google Scholar 

  82. 82.

    Beacham TD. Meristic and morphometric variation in pink salmon (Oncorhynchus gorbuscha) in southern British Columbia and Puget Sound. Can J Zool. 1985;63(2):366–72. https://0-doi-org.brum.beds.ac.uk/10.1139/z85-056.

    Article  Google Scholar 

  83. 83.

    Slechtova V, Slechta V, Berrebi P. Diversité génétique du genre Barbus dans les rivières Tchèques et Slovaques: résultats préliminaires. Cah Ethol. 1993;13:207–10.

    Google Scholar 

  84. 84.

    Crespin L, Berrebi P. L'hybridation naturelle entre le barbeau commun et le barbeau méridional en France: compte rendu de dix années de recherche. Bull Fr Pêche Piscic. 1994;334:177–89.

    Article  Google Scholar 

  85. 85.

    Scribner KT, Page KS, Bartron ML. Hybridization in freshwater fishes: a review of case studies and cytonuclear methods of biological inference. Rev Fish Biol Fisher. 2000;10(3):293–323. https://0-doi-org.brum.beds.ac.uk/10.1023/A:1016642723238.

    Article  Google Scholar 

  86. 86.

    Tsigenopoulos CS, Rab P, Naran D, Berrebi P. Multiple origins of polyploidy in the phylogeny of southern African barbs (Cyprinidae) as inferred from mtDNA markers. Heredity. 2002;88(6):466–73. https://0-doi-org.brum.beds.ac.uk/10.1038/sj.hdy.6800080.

    CAS  Article  PubMed  Google Scholar 

  87. 87.

    Freyhof J, Lieckfeldt D, Pitra C, Ludwig A. Molecules and morphology: evidence for introgression of mitochondrial DNA in Dalmatian cyprinids. Mol Phylogenet Evol. 2005;37(2):347–54. https://0-doi-org.brum.beds.ac.uk/10.1016/j.ympev.2005.07.018.

    CAS  Article  PubMed  Google Scholar 

  88. 88.

    Touil A, Casal-Lopez M, Bouhadad R, Doadrio I. Phylogeny and phylogeography of the genus Luciobarbus (Haeckel, 1843) in Algeria inferred from mitochondrial DNA sequence variation. Mitochondrial DNA A DNA Mapp Seq Anal. 2019;30(2):332–44. https://0-doi-org.brum.beds.ac.uk/10.1080/24701394.2018.1526928.

    CAS  Article  PubMed  Google Scholar 

  89. 89.

    Bianco PG, Ketmaier V. Anthropogenic changes in the freshwater fish fauna of Italy, with reference to the central region and Barbus graellsii, a newly established alien species of Iberian origin. J Fish Biol. 2001;59(Supplement A):190–208.

    Article  Google Scholar 

  90. 90.

    Wightman EM. Topographic survey in the Liri Valley, southern Lazio, Italy. Curr Anthropol. 1978;19(2):389–90. https://0-doi-org.brum.beds.ac.uk/10.1086/202087.

    Article  Google Scholar 

  91. 91.

    Bonaparte CL. Iconografia della Fauna Italica, per le quattro classi degli animali vertebrati. Tomo 3, Pesci. Roma: Salviucci; 1832-1841. https://0-doi-org.brum.beds.ac.uk/10.5962/bhl.title.70395.

    Book  Google Scholar 

  92. 92.

    De Filippi F. Pesci finora osservati in Lombardia. Notizie Naturali e Civili su la Lombardia. 1844;1:389–406.

    Google Scholar 

  93. 93.

    Canestrini G. Prospetto critico dei pesci d’acqua dolce d’Italia. Archivio per la Zoologia, l’Anatomia e la Fisiologia. 1866;4:47–187.

    Google Scholar 

  94. 94.

    Canestrini G. Fauna d’Italia–Parte terza. Pesci. Milano: Vallardi; 1874.

    Google Scholar 

  95. 95.

    Giglioli EH. Elenco dei Mammiferi degli Uccelli e dei Rettili ittiofagi appartenenti alla fauna italica e catalogo degli anfibi e dei Pesci italiani. Firenze: Stamperia Reale; 1880.

    Google Scholar 

  96. 96.

    Scotti L. La distribuzione dei pesci d'acqua dolce in Italia. Giornale italiano di pesca e acquicoltura. 1898;2:42–8.

    Google Scholar 

  97. 97.

    Supino F. I pesci d’acqua dolce d’Italia. Milano: Hoepli; 1916.

    Google Scholar 

  98. 98.

    Brunelli G, Chiappi T. I pesci d’acqua dolce. In: Monti R, Brunelli G, Chiappi T, editors. La pesca nei mari e nelle acque interne d’Italia. Notiziario tecnico e legislativo e repertorio della industria e del commercio dei prodotti pescherecci, vol. 2. Roma: Istituto Poligrafico dello Stato; 1931.

    Google Scholar 

  99. 99.

    Koller O. Eine kritische Übersicht über die bisher beschriebenen mittel-und südeuropäischen Arten der Cyprinidengattung Barbus Cuv. Sitzungsberichte der Kaiserlichen Akademie der Wissenschaften– mathematisch-naturwissenschaftliche Classe. 1926;135:167–202.

    Google Scholar 

  100. 100.

    Bini G. I pesci delle acque interne d’Italia. Roma: Garzanti; 1962.

    Google Scholar 

  101. 101.

    Tortonese E. Osteichthyes (Pesci ossei). Parte Prima. Fauna d’Italia, 10. Bologna: Calderini; 1970.

    Google Scholar 

  102. 102.

    Kottelat M. European freshwater fish. Biologia Bratislava. 1997;52(Supplement 5):1–271.

    Google Scholar 

  103. 103.

    Cuvier GLCFD, Valenciennes A. Histoire naturelle des poissons, vol. 16. Paris: Bertrand, Strasbourg: Levrault; 1842.

    Google Scholar 

  104. 104.

    Heckel J, Kner R. Die Süsswasserfische der östreichischen Monarchie. Leipzig: Engelmann; 1858.

    Google Scholar 

  105. 105.

    Nardo GD. Prospetti sistematici degli animali delle province venete e del mare Adriatico e distinzione delle specie in gruppi relativi alla loro geografia fisica ed all'interesse economico statistico che presentano. Venezia: Antonelli; 1860.

  106. 106.

    Steindachner F. Zur Fischfauna von Isonzo. Verhandlungen der Zoologisch-Botanischen Gesellschaft. 1865;15:141–2.

    Google Scholar 

  107. 107.

    Günther AC. Catalogue of the fishes in the British museum, vol. 7. London: British Museum; 1868.

    Google Scholar 

  108. 108.

    Vuković TI, Ivanović B. Slatkovodne ribe Jugoslavije. Sarajevo: Zemaljski muzej Bosne i Hercegovine; 1971.

    Google Scholar 

  109. 109.

    Gridelli E. I pesci d'acqua dolce della Venezia Giulia. Trieste-Udine: Consorzio per la tutela della pesca nella Venezia Giulia; 1935.

    Google Scholar 

  110. 110.

    Karaman M. Süsswasserfische der Türkei. 8 Teil. Revision der Barben Europas, Vorderasiens und Nordafrikas. Arch Fisch. 1971;22:165–74.

    Google Scholar 

  111. 111.

    Bianco PG. I Pesci d’acqua dolce d’Italia: note su un recente contributo. Atti della Società italiana di scienze naturali e del Museo civico di storia naturale di Milano. 1988;129:146–58.

    Google Scholar 

  112. 112.

    Cuvier GLCFD. Le Règne Animal, distribué d’après son organisation, pour servir de base à l’histoi-re naturelle des animaux et d'introduction à l’anatomie comparée. 2nd ed. Paris: Déterville; 1829.

    Google Scholar 

  113. 113.

    Costa OG. Fauna del regno di Napoli, ossia enumerazione di tutti gli animali che abitano le diverse regioni di questo regno e le acque che le bagnano contenente la descrizione de’ nuovi o poco esattamente conosciuti. Pesci. Parti I-III. Napoli: Azzolino; 1830-1857.

    Google Scholar 

  114. 114.

    Vinciguerra D. Guida del Museo di Zoologia della R. Università di Roma–Fauna locale–Specie animali della provincia di Roma esistenti nella nuova collezione. Parte 3. Pesci. Bollettino del Museo di Zoologia della Reale Università di Roma. 1890;1:1–39.

    Google Scholar 

Download references

Acknowledgements

We thank Fran Janžekovič, Samo Podgornik and Danilo Puklavec for specimen collection in Slovenia; Marco Guida and Antonella Giorgio for the support to every step of sampling in Southern Italy; Franco Anacardi, Andrea Lenuzza, Luca Magnani, Simona Piccini, and Marco Rizzoli for contribution in sampling design; Giorgio Marchesan, Giacomo Serafini, Valentina Mingazzini, Michela Zattini and Mariachiara Zanichelli for early molecular analysis; Amir Pourshaban for knowledge sharing on biogeographical and phylogenetic topics. A special thank goes to Matteo Monti for the professional and friendly contribution in anatomical pictures taking. We are grateful to Francesco Zaccanti for the scientific support in every step of the research, to Gaetano Caricato for the location of the sampling site for type collection and to Antonio Bonfitto for advices in taxonomic tasks. We want to especially thank Pier Giorgio Bianco for the long talks and the detailed lectures and, in a nutshell, for his entire work, from which this article took the maximum inspiration. Finally, this work was significantly improved thanks to two anonymous reviewers.

Funding

No funding was received for the present study.

Author information

Affiliations

Authors

Contributions

GR conceived and designed the study. GR, FP and RF supervised it during all stages. GR, GZ, AM, SDB, and MV carried out sampling campaigns. GZ, MV, and PM performed wetlab experiments. FP and GZ analyzed sequences. FP carried out the phylogenetic analysis. AM prepared the samples for museum conservation, carried out morphological analysis, and took the photograph for Fig. 2. GR performed taxonomic analysis and wrote Appendices. GR, FP, GZ, AM, and RF drafted the manuscript; GR, FP, and RF reviewed, revised and finalized the manuscript; all authors have read and approved the final manuscript.

Corresponding author

Correspondence to Federico Plazzi.

Ethics declarations

Ethics approval and consent to participate

All samplings and procedures were carried out in accordance with relevant guidelines and regulations, with special reference to the positive opinion 110717 issued on 05/22/2019 by the Comitato per il Benessere degli Animali of the University of Bologna and to the clearance 698/2019-PR issued on 10/14/2019 by the licensing committee in charge, the Italian Ministero della Salute – Direzione Generale della Sanità Animale e dei Farmaci Veterinari – Ufficio 6, in accordance to the Italian legislative decree 26/2014, art. 31.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

Sampling sites distribution, along with coordinates in EPSG 32632 reference systems.

Additional file 2.

Dichotomous key for barbels of Apennine Peninsula and Adriatic basin of Slovenia and Croatia [15].

Additional file 3.

GenBank Accession Numbers and geographical abundancies of sequences used for single-marker minimum spanning network computations (see Additional file 6). TL, Tuscany-Latium district; PV, Padano-Venetian district; NAAC, northernmost-Adriatic part of Apulia-Campania district; TSAAC, Tyrrhenian and southernmost-Adriatic parts of Apulia-Campania district; DAN, Danubian district; AccNum, GenBank Accession Number.

Additional file 4.

Dataset, taxonomy, river/watershed location, sampling coordinates in EPSG 32632 reference systems and GenBank Accession Numbers.

Additional file 5.

Original phylogenetic trees in NEXUS format as computed by IQ-TREE, MrBayes, and RAxML, respectively. Each OTU is labeled with the resulting clade (see Fig. 2) and with GenBank Accession Numbers (Additional file 4).

Additional file 6.

Tests for chance occurrence of reciprocal monophyly. All analyses were carried out using the Species Delimitation Plugin of Geneious Prime® 2021.0.3.

Additional file 7.

Single-marker minimum spanning networks computed from sequences and geographical abundancies listed in S4. a, cytb; b, D-loop. Asterisks mark samples that was not possible to unambiguously assign to an ichthyogeographical district (see also Fig. 4).

Additional file 8.

Meristic data of Barbus plebejus and B. tyberinus from [28, 33, 34] and morphometric and meristic data of B. oscensis (i.e., the B. tyberinus TSAAC clade).

Additional file 9.

Meristic data of Barbus plebejus (N = 153) and Barbus tyberinus (N = 168) from [28, 33, 34] and Barbus oscensis (i.e., the B. tyberinus TSAAC clade; N = 6) from original counts. Vertical lines: observed range; dithered boxes and horizontal lines: usual values. Characters with no evident variability between subspecies are not shown.

Additional file 10.

Tukey post-hoc test results. The table lists the Q value; the p-value is shown in parentheses. ***, p < 0.001; **, p < 0.01. NAAC, Barbus sp. clade 4; TSAAC, Barbus sp. TSAAC clade; TL, B. tyberinus TL clade; PV, B. plebejus PV clade; DAN, B. barbus clade.

Appendices

Appendix 1

Review of the original description of Barbus plebejus

Family Cyprinidae Rafinesque 1815

Subfamily Barbinae Bleeker 1859

Genus Barbus Cuvier and Cloquet 1816

Barbus plebejus Bonaparte 1839

Synonymy

Barbus plebejus [91]: vol 3, no. 25, pl. 110, fig. 1 (Lake Como); Valenciennes, in [103]: 139, pl. 462 (Po River); [104]: 82, fig. 38 (near Milan and Triest); [105]: 72 (upper Adriatic Province); [106]: 142 (Isonzo-Soča River); [93]: 72 (Partim: PV); [107]: 88 (Lago Maggiore and Dalmatia); [108]: 179 (Isonzo-Soča, Zrmanja and Krka Rivers, Adriatic watershed of former Yugoslavia); [35]: 198–202 (Partim: PV); [23]: 461 (PV); [28]: 307, figs. 2-4 (Partim: Northern Italy, Adriatic Watershed of Central Italy north of Vomano excluded, Istria, Dalmazia north to the river Krka included); [14]: 128 (PV including the Zrmanja and Krka Rivers in Croatia); [33]: 339 (PV); B. eques [91]: vol 3, no. 25, pl. 110, fig. 2 (Partim: Northern Italy); B. fluviatilis non Valenciennes 1842: [92]: vol. 1, 594 (Lombardy); B. barbus plebejus [99]: 172, Tab I, figs. 6,7 (Italy); [109]: 100 (Northern Italy); [101]: X, 274, fig. 11A (Partim: PV); B. plebejus plebejus [110]: 200 (partim: PV); [111]: 154 (Partim: PV).

Type specimens

The type specimens defined by Bianco [33] for B. plebejus within syntypes of Bonaparte’s original collection, sampled in PV.

Lectotype

Lectotype of B. plebejus Bonaparte: Collection of the Academy of Natural Sciences of Philadelphia No. ANSP 6183, from Lake Como.

Paralectotypes

5 paralectotypes of B. plebejus Bonaparte: Collection of the Academy of Natural Sciences of Philadelphia No. ANSP 6184–6188, from Lake Como.

Type locality

Lake Como.

Distribution

PV

Etymology

The specific name plebejus is the Latin word for plebeian, meaning common.

Diagnosis

Due to overlapping, meristic characters have low or no reliability for morphological discrimination from other species; therefore, the diagnosis is mainly based on genetic differences in mitochondrial markers (cytb and D-loop).

Morphology

Count and measurement are given in S8; last unbranched ray of the dorsal fin moderately ossified and with fine serration; ossification reduced in the distal portion of the ray, serration tend to become less evident as the fish becomes older.

Remarks

B. plebejus Bonaparte 1839 is a fluvio-lacustrine barbel species perfectly matching the PV, parapatric to B. tyberinus (TL range) and to the new species proposed in this work (TSAAC range). Originally mentioned (as nomen nudum) under the name of B. plebejus by Cuvier [112], the first valid description of this taxon, B. plebejus Bonaparte, 1839 [91], is exclusively based on type specimens from PV; hence, no revision or emendation is needed.

Appendix 2

Review of the original description of Barbus tyberinus

Family Cyprinidae Rafinesque 1815

Subfamily Barbinae Bleeker 1859

Genus Barbus Cuvier and Cloquet 1816

Barbus tyberinus Bonaparte 1839

Synonymy

Barbus fluviatilis tyberinus [91]: fasc. 25, pl. 110 fig. 3 (Tiber River near Rome); Barbus eques [91]: fasc. 25 pl. 110 fig. 2 (partim: River Arno in Tuscany Region); Valenciennes in [103]: 141 (Florence); Barbus canalii Valenciennes in [103]: 143–144 (Topino River, misspelled Topico); Barbus fucini [113]: folio 79, pl 11, figs 1-7 (lake Fucino) [113]; sensu Bianco [23]: 463 (partim: TL); Barbus fluviatilis plebejus (non Bonaparte 1839) [92] (partim: TL); Barbus plebejus (non Bonaparte 1839), [93]: 72–77 (partim: TL); [94]: 11 (partim: TL); [95]: 43 (partim: Castelnuovo in Garfagnana, Arno River at Florence, Casentino region); [114]: 28 (Rome province); [96]: 42–43 (partim: TL); [97]: 40–41 (partim: TL); [98]: 274–275 (partim: TL); [35]: 198–202 (partim: TL); Barbus Barbus plebejus (non Bonaparte 1839) [99]: 172 (partim: TL); [100]: 19 (partim: TL); [101]: 274–275 (partim: TL); Barbus tyberinus (non Bonaparte 1839) [28]: 313–318, figs. 4b, 6–7 (Partim: TL); [102]: 50–51 (Partim: TL); [14]: 234 (Partim: TL); [34]: 427 (Partim: TL).

Type specimens

The type specimens defined by Bianco [34] within syntypes of Bonaparte’s and Costa’s original collections, sampled in TL.

Lectotype

Lectotype of B. fluviatilis tyberinus: Collection of the Academy of Natural Sciences of Philadelphia No. ANSP 6152 (Bonaparte’s original number 426), from river Tevere. 122 mm SL (166 mm TL); LL 55; 15.5 row of scales above and 9.5 below LL; 26 circumpeduncular scales; D 8; 8 (6 + 2) gill rakers.

Paralectotypes

27 paralectotypes of B. fluviatilis tyberinus: Collection of the Academy of Natural Sciences of Philadelphia No. ANSP 6153–6179, from river Tevere. 50–188 mm SL; 47–61 LL; 12.5–15.5 row of scales above and 7.5–9.5 below LL; 24–28 circumpeduncular scales; D 8; 8–10 (6–8+ 1–3) gill rakers. Lectotype of B. canalii. Muséum national d’Histoire Naturelle in Paris No. MNHN 1412; from river Topino. 146 mm SL; 53 LL; 12.5 row of scales above and 8.5 scales below LL; 24 circumpeduncular scales; 8 (7 + 1) gill rakers. 5 paralectotypes of B. eques Bonaparte: Collection of the Academy of Natural Sciences of Philadelphia No. ANSP 6144, 6146–6149, from river Arno.

Type locality

River Tevere near Rome.

Distribution

TL

Etymology

The specific name tyberinus is derived from the name of the River Tevere, the type locality.

Diagnosis

Same considerations as for B. plebejus.

Morphology

Count and measurement are given in S8; last unbranched ray of the dorsal fin poorly ossified and with fine serration; ossification reduced in the distal portion of the ray, serration tend to become undistinguished as the fish becomes older.

Remarks

B. tyberinus Bonaparte 1839 is a fluvio-lacustrine barbel species perfectly matching the TL, parapatric to B. plebejus (PV range) and to the new species proposed in this work (TSAAC range); the presence of the species in Liri-Garigliano basin, the northernmost watershed of Apulia-Campania district at the boundary with Tuscany-Latium district is likely to be naturally due to the geomorphological phenomena.

Conversely, on the basis of morphological analysis [28, 34] B. tyberinus was considered including every non-PV fluvio-lacustrine populations; therefore, its range was assumed covering also AC.

The first valid description of this taxon, B. fluviatilis tyberinus Bonaparte 1839 is based on type specimens all sampled in TL as individuated by Bianco [28], therefore the original description of Bonaparte must be considered correct; anyway, being not a subspecies of B. fluviatilis, rather a different species, Bianco [28] identified as valid the combination B. tyberinus Bonaparte 1839. Correctly, following article 51.3.2. of ICZN, no parenthesis was added around author’s name and date. Moreover, this combination is also a synonym of B. eques Valenciennes in Cuvier and Valenciennes 1842, of B. canalii Valenciennes in Cuvier and Valenciennes 1842 and B. fucini Costa 1853.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Rossi, G., Plazzi, F., Zuffi, G. et al. Mitochondrial phylogeny and taxonomic revision of Italian and Slovenian fluvio-lacustrine barbels, Barbus sp. (Cypriniformes, Cyprinidae). BMC Zool 6, 8 (2021). https://0-doi-org.brum.beds.ac.uk/10.1186/s40850-021-00073-x

Download citation

Keywords

  • Phylogeography
  • Fluvio-lacustrine barbels
  • Mitochondrial markers
  • Ichthyogeographic districts
  • Taxonomy