Next Article in Journal
Effect of Local Winds on Salinity Intrusion in the Columbia River Estuary
Next Article in Special Issue
Seasonal Changes in Upper Thermal Tolerances of Freshwater Thai Fishes
Previous Article in Journal
An Investigation on the Effect of Leakages on the Water Quality Parameters in Distribution Networks
Previous Article in Special Issue
Environmentally Realistic Waterborne Atrazine Exposure Affects Behavior in Poecilia latipinna
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Past and Present Environmental Factors Differentially Influence Genetic and Morphological Traits of Italian Barbels (Pisces: Cyprinidae)

1
Department of Theoretical and Applied Sciences, University of Insubria, 21100 Varese, Italy
2
CNR Water Research Institute, 28922 Verbania, Italy
3
Department of Chemistry, Biology and Biotechnologies, University of Perugia, 06123 Perugia, Italy
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Water 2023, 15(2), 325; https://doi.org/10.3390/w15020325
Submission received: 12 December 2022 / Revised: 6 January 2023 / Accepted: 10 January 2023 / Published: 12 January 2023

Abstract

:
Local adaptation and phenotypic plasticity can lead to environment-related morphological and genetic variations in freshwater fish. Studying the responses of fish to environmental changes is crucial to understand their vulnerability to human-induced changes. Here, we used a latitudinal gradient as a proxy for past and present environmental factors and tested its influences on both genetic and morphological patterns. We selected as a suitable biogeographic model, the barbels, which inhabit 17 Adriatic basins of the central-southern Italian Peninsula, and explored association among attributes from genetic, morphological, and environmental analyses. The analysis of the mitochondrial DNA control region evidenced a southward significant increase in the number of private haplotypes, supporting the isolation of the southernmost populations related to the Mio-Pleistocene events. In contrast, morphology was mainly affected by changes in the present environmental conditions. Particularly, the number of scales and fish coloration were clearly associated to latitude, and thus thermal and hydrological conditions. Other morphometric and functional traits varied under the selective pressure of other environmental factors like elevation and distance from headwater. These results highlight the sensitivity of barbels to climate changes, which can serve as a basis for future eco-evolutionary and conservation studies.

Graphical Abstract

1. Introduction

In lotic ecosystems, environmental conditions (e.g., hydrological regime, channel and substrate characteristics, and water physico-chemical parameters) play an active role in modelling fish morphology relevant to motion and trophic adaption [1,2,3,4]. Different availability of habitats and resources may determine different ecological niches for fish populations inhabiting distinct rivers. According to the niche variation hypothesis [5], niche width variation at the population level, mainly resulting from the increased variation at the between-individual level, is accompanied by greater morphological and genetic variations [6,7,8,9]. Phenotypic variations in fish may result in fixed traits (i.e., traits that do not change as fast as the environment [10]), which arise from adaptation following natural selection or other evolutionary forces, or in plastic response to changing environments (phenotypic plasticity) [11] to the point that both populations of fish belonging to the same species or to different species may display both morphological and genetic differences reflecting the environmental conditions of the rivers they inhabit [12,13,14,15].
Global changes, particularly the current climate change, are modifying both river habitats and resources [16,17], and could thus lead to changes in the ecological niches of fish populations. Both air temperature increases and a decrease in low-flow values are causing a relevant increment of water temperature in many rivers [18,19,20,21,22], which is one of the, if not the, most important drivers of phenotypic variations in fish [23,24,25]. Since temperature and other environmental factors such as hydrology, photoperiod, length of growing season, and productivity vary with latitude [26], the study of morphological and genetic differences of fish populations along latitudinal gradients may provide important indications of the effects of current climate changes.
Variation in meristic traits in fish is mainly expected as a consequence of the current temperature increase and of changes in other environmental parameters such as salinity and dissolved oxygen [27]. In contrast, genetic diversity, particularly in neutral (i.e., not under selection) DNA regions, such as those commonly used for population genetic studies, typically reflect historical demographic evolution due to past climatic and geological changes e.g., [28]. Nevertheless, both genetic and phenotypic variabilities can also be altered by recent anthropogenic activities such as the introduction of exotic species (e.g., following hybridization [29]) and river damming [30,31,32]. Thus, separating genetic and phenotypic variations induced by local adaptation developed during historical evolution from the variations deriving by relatively recent events, provides valuable insights for both eco-evolutionary processes and for a conservation perspective [33].
Despite its relatively small area (approximately 300,000 km2), the Italian Peninsula is characterized by a high environmental heterogeneity due to the complex structure of its orography (40% of its area is mountainous), and relevant biogeographical and geological events [34]. The wide latitudinal extension of the peninsula (approximately 10 degrees) places it between different climatic zones (continental, temperate, and Mediterranean zones), which affect the hydrography and the hydrological regime of Italian watercourses [35]. In particular, watercourses on the Adriatic side of central-southern Italy are influenced by the presence of the Apennines, are relatively limited in length (<100 km), and flow directly into the sea. Moving southward, these rivers have a hydrological regime increasingly prone to intermittency [36], with climate change also contributing to water scarcity via altered runoff [37]. Changes in conditions of the river environment along the latitudinal gradient may have influenced both morphological and genetic traits of fish, thus making Italy an ideal geographical model to investigate the effects of changes in temperature and hydrological conditions on freshwater fish.
Barbels (genus Barbus Daudin, 1805) are among the most widespread primary freshwater fishes inhabiting different geographic districts of Europe including the Italian Peninsula [38]. They consist of a diverse group of medium to large size cyprinids with a high level of endemism and varied ecology [38,39,40]. Four of the six native barbel species are fluviolacustrine inhabiting middle and lower reaches of the rivers: Padanian barbel Barbus plebejus Bonaparte, 1839, Tiber barbel Barbus tyberinus Bonaparte, 1839, Barbus samniticus Lorenzoni et al., 2021, and Barbus fucini Costa, 1853. These species are vicariant in different ichthyogeographic districts: Padano-Venetian, Tuscany-Latium, northern Apulia-Campania, and southern Apulia-Campania, respectively [38,41]. Specifically, three species inhabit the basins flowing into the Adriatic Sea: B. plebejus is present from the Krka River in Croatia to the Tronto River in Italy including the Po River basin (i.e., the largest Italian basin, with a 71,000 km2 area), B. samniticus is present from the Vomano River to the Biferno River, and B. fucini was detected in the Fortore and Ofanto basins [38]. Speciation in the mentioned districts and basins started following the Messinian “Lago Mare” phase of the Mediterranean Sea (6.1–5.1 MYA–Miocene), and species distribution was further influenced by Pleistocenic geomorphological events. After the last glacial maximum (Würm glacial phase 15,000–18,000 years ago), when the paleo-Po basin reached the current Vomano basin and biological homogenization occurred, the independent basins flowing into the middle-Adriatic Sea started to isolate from the south to the north due to marine regression [42]. Only a few characteristics allow the discrimination between the four vicariant Italian fluviolacustrine barbels [38,43,44,45]. Also, the trait mainly examined to discriminate among B. plebejus and B. tyberinus in the field, i.e., the number (and therefore size) of scales along the lateral line [46,47], can lead to erroneous species attribution [43,45] since it is likely influenced by environmental conditions like for other meristic traits [38].
Here we used fluviolacustrine barbel populations distributed along the Italian Peninsula to disentangle the influence on morphology and genetic diversity of Mio-Pleistocene isolation of watersheds compared to recent habitat features. This was done by comparing genetic and phenotypic variability along a latitudinal gradient that can be seen as a proxy of changes in river environmental conditions. Thus, we present a study on both the genetic and morphological characterization of barbel populations inhabiting 17 rivers distributed from north to south along the Adriatic side of central-southern Italy (Figure 1) including three out of four Italian ichthyogeographic districts (Padano-Venetian, northern and southern Apulia-Campania). We posited that: i) the genetic variability of the studied populations is related to the geomorphological evolution of the river basins they inhabit with no relation with current environmental conditions, and ii) both intra- and inter-specific morphological variations are related to local environmental conditions, representing a response to the present habitat features potentially varying due to climate change. For this purpose, we employed different approaches including morphological and molecular analyses, and environmental investigation.

2. Materials and Methods

2.1. Study Area and Barbel Samples

In previous works [41,45], we genetically and morphologically characterized 24 barbel populations from sampling sites within Adriatic basins in central-southern Italy. Since some populations were affected by the invasion of the exotic fluviolacustrine European barbel Barbus barbus L. hybridizing with the native Italian species [48,49,50,51], thus eroding their genotypes and altering their morphologies [44,51], we tried to find additional purebred populations in vicinity to increase the dataset size. Thus, in August 2020 we sampled barbel populations (approximately 20–30 individuals per population) from six Italian rivers (Borello, Marecchia, Foglia, Salinello, Tordino, and Trigno) located on the Adriatic side of central-southern Italy (Figure 1 and Table 1).
Following their collection by electric fishing, fish were anaesthetized (MS-222) and photographed for morphological analysis from their left side using a Nikon D300 camera (24–85 mm lens) positioned by means of a tripod on a table. Scales were removed from the left side of each fish and stored in 70% ethanol for age determination. Scalimetry was performed by two operators independently [52] under a stereomicroscope coupled with a camera with images stored within an archive built using the image analysis system IAS 2000 (QEA’s IASLab software). Determination done by two independent readers allowed to control for potential errors in age determination by applying a quality control approach for which the independent determinations were cross-checked and in case of mismatches, the scale was redetermined or discarded. Although age determination tends to overestimate fish ages [53], this approach allowed us to select fish of relatively similar age (see Section 2.3). A caudal-fin clip was taken, and samples were stored in 100% ethanol for genetic analysis. Following their recovery, all fish were released at their capture sites.

2.2. Genetic Data

Total genomic DNA was extracted from selected individuals using a proteinase K digestion, followed by sodium chloride extraction and ethanol precipitation (i.e., salting out method [54]). First, for excluding populations invaded by B. barbus, the mitochondrial cytochrome b of a few individuals (4–10) per population was amplified through PCR using primers L15267 and H16461 [55], and a fragment of 1121 bp was analyzed. Then, a fragment (867 bp) of the mitochondrial control region (D-loop) locus from 13–21 individuals belonging to purebred populations was amplified using D-loops × F and D-loopd × R [56,57] primers pair (Table 1; Borello and Marecchia populations were excluded from further analysis due to the presence of B. barbus).
Amplifications were performed for both loci using the Q5 High-Fidelity Master Mix (New England Biolabs Inc., Ipswich, Massachusetts, MA, USA) in a 10 μL reaction volume containing 10 ng of template DNA and 0.25 μM of each primer. The PCR profile was set with an initial denaturation at 98 °C for 3 min followed by 35 cycles of 30 s at 98 °C, 90 s at 55 °C, and 90 s at 72 °C, with a final extension step at 72 °C for 10 min. Amplicons were then purified using EuroSAP PCR Enzymatic Clean-up kit (EuroClone s.p.a, Pero, Milan) and sequenced with single pass by Macrogen Europe B. V. on a 3730XL DNA sequencer. Sequences were aligned using CLUSTAL-W [58], as implemented in Bioedit software [59], and further checked manually to eliminate ambiguities. All the haplotypes detected in this study were deposited in GenBank (ON661665-ON661675 for the D-loop mtDNA). Species attribution was based on the genetic results presented in this and previous studies [41,45].

2.3. Morphological Data

Morphological data were collected from photos taken from adult fish of similar age (2+, 3+, and 4+) to reduce potential allometric bias (Table 1). The sex of the fish was not determined due to the non-lethal nature of our sampling. Thus, potential sources of variation due to sexual dimorphisms were not controlled even though these should be limited to the anal fin [38,47]. For geometric morphometrics, 33 landmarks (Figure 2a) were captured from each photo using the R [60] Geomorph function “digitize2d” [61].
For common morphological analysis, a total of 40 traits (25 morphometric and 15 meristic traits) were analyzed. In addition to the four meristics showed in Figure 2b, coloration (presence or absence of dots on body–Db, scales–Ds, dorsal fin–Ddf, anal fin–Daf, and caudal fin–Dcf, spots on body–Sb, external pigmentation on scales–EPs, and color of ventral–Cvf, anal–Caf, dorsal–Cdf, and caudal–Ccf–fin) was examined. Fin color was determined independently by two observers following the descriptions reported in previous studies [39,46,47]. Morphometric measurements were standardized to the overall mean standard length [62] to further reduce the effects of size and allometry.
For the analysis of functional traits, 10 ratios of morphological measurements mainly associated with food acquisition and locomotion were taken into consideration [63] (Table 2 and Figure 2b). The knowledge of the relationship between functional morphology and environmental gradients allows for better predictions of the organisms’ responses to environmental change [64,65].

2.4. Environmental Data

At each site in Table 1, the percentage coverages of different riverbed substrates, different mesohabitats, shading, plants, and refuges were detected as well as the average depth and width (Table S1). Moreover, geological elements and land use were calculated at basin scale (Table S1) using QGIS [74], and geolithological and land use maps [75,76]. The climate characterizing the different watersheds was determined by downloading data on temperature and precipitation (Table S1) recorded at the meteorological station closest to the barbel sampling sites from SCIA [77]; the average values for the period 1998–2018 were considered for the analysis. In order to consider the water chemistry and benthos status of barbel habitats sensu Water Framework Directive (2000/60/EC), the mean values of LIMeco and STAR_ICMi indices [78,79] relative to the river reaches near the sampling sites (Table S1) were collected from reports by regional environmental protection agencies for the period 2015–2017 [80,81,82,83]. The total river length and the quotient of average flow at the river outlet and total basin area (Q/A) were also retrieved [84,85,86,87] and the latter was used as indicator of the river hydrological regime (Table S1).

2.5. Data Analysis

2.5.1. Genetic Data

Number of haplotypes (n) and polymorphic sites (S) were calculated for each barbel population, and both nucleotide (π) and haplotype (H) diversity indices were estimated using DnaSP v6.12.01 [88]. Moreover, the relative frequency of population-specific haplotypes (PSh) was used as a proxy of genetic isolation. To visualize the geographical distribution of haplotypes and relative frequencies, a median-joining network was created through the software PopArt [89].
Linear regression of all genetic parameters with latitude was tested using XLSTAT2014 software, and only significant regressions (p ≤ 0.05) were reported in this study.

2.5.2. Morphological Data

Generalized Procrustes analysis was performed to remove non-shape variations, occurring through variation in position, orientation, and size, as implemented in MorphoJ software [90]. Through the same software, shape variations within the dataset were analyzed by canonical variate analysis (CVA) using populations as a priori groups. Mahalanobis distances (MDs) between populations in the CVA were calculated using permutation tests (10,000 replicates) and then MDs were represented graphically by non-metric multidimensional scaling (NMDS).
Linear regression of NMDS axes and all morphological and functional traits with latitude was tested using XLSTAT2014 software, and only significant regressions (p ≤ 0.05) were reported. Regressions were examined using both average values and all measured values of the mentioned traits. For meristic traits, the relative frequency of each population was considered.

2.5.3. Environmental Data

Correlations among latitude and the other environmental factors considered in this study were determined using Pearson product moment correlation coefficient. Then, to detect the influence of all environmental factors on morphological and genetic traits of barbel populations, a redundancy analysis (RDA) was carried out. Environmental data were adopted as explanatory variables, and morphometric, functional, and genetic traits as response variables. To reduce the number of variables, principal component analysis (PCA) was performed with data on basin geology and land use, substrate and mesohabitat of sampling sites, climate, and morphometric traits, and the first component (also the second component for morphometric traits) was used in the RDA. All variables were standardized prior to the RDA. The analyses were performed using XLSTAT2014 software.

3. Results

3.1. Genetics of Barbel Populations along a Latitudinal Gradient

Among the six new populations, only four (i.e., Foglia, Salinello, Tordino, and Trigno) were purebred and thus were retained in this study along with other 13 purebred populations from previous works (Table 1). The median-joining network built on the 17 populations identified three main haplotype groups (Figure 3) separated by 16 to 22 mutation steps and corresponding to the native species B. plebejus, B. samniticus, and B. fucini. Among the newly analyzed populations, barbels found in Foglia, Salinello, and Tordino rivers belonged to B. plebejus species, while Trigno population consisted of B. samniticus specimens (Table 1). Barbus plebejus is thus present from Foglia (S1) to Tordino (S11), B. samniticus from Vomano (S12) to Trigno (S15), and B. fucini in Tappino (S16) and Ofanto (S17) rivers. From B. plebejus populations, 34 different haplotypes were detected, 13 from B. samniticus populations and six from B. fucini populations, for a total of 53 haplotypes. The number of haplotypes per population ranged from two to 11 (Table 3). The Barbus plebejus haplogroup presented a dumbbell pattern with the two most frequent haplotypes separated by two mutation steps. Several B. plebejus haplotypes were shared across different populations with only two populations, namely Tenna (S7) and Aso (S8), characterized, respectively, by two exclusive haplotypes (Figure 3). The Barbus samniticus network was characterized by a heterogeneous pattern where the two most frequent haplotypes were connected by up to 12 mutation steps, with Vomano (S12) including four population-specific haplotypes. Haplotypes detected from B. fucini were all population-specific haplogroups where four mutation steps were postulated between the major types (Figure 3).
The total number of haplotypes (n) detected in each barbel population significantly increased with latitude while the relative frequency of population-specific haplotypes (PSh) significantly decreased (Figure 4). No significant associations with the latitudinal gradient were observed for the number of polymorphic sites (S), and haplotype (H) and nucleotide (π) diversity indices. However, accounting for the two diversity indices (Table 3), some considerations could be done: (i) populations from S1 to S4 displayed the highest values of both H (0.68–0.87 on average) and π (0.002 on average), indicating populations with a good genetic structure compared to the other populations, and (ii) populations S7, S12, and S13, displayed the lowest values of both H (0.10–0.30 on average) and π (0.0001–0.0005 on average), probably indicating a recent marked demographic contraction.

3.2. Morphology of Barbel Populations along a Latitudinal Gradient

Only a few morphological differences were detected along the CVA axes (Figure 5a). These differences were mainly associated to mouth orientation, head length, and the caudal fin shape. Specimens belonging to S15 (Trigno) population were the most morphologically different from the others (Figure 5b), i.e., the average Mahalanobis distance (MD) from the other populations was the highest (9.9). The maximum MD value (12.6) was detected among S15 and S13 (Giardino) populations despite belonging to the same species (B. samniticus) and located at close latitudes. No significant relationships were observed among the axes of NMDS displaying MDs and latitude.
However, considering morphological traits singularly, some significant linear regressions with latitude were found (Figure 6). On average, the maximum caudal fin depth (CFD) was negatively related to latitude, while the meristic traits associated to the number of scales (i.e., NSLL, NSALL, and NSULL) were positively related to latitude. Also, the relative frequency of individuals having an orange caudal fin (Ccf, the others had an orange/grey caudal fin) increased with latitude while the relative frequency of individuals having dots on scales (Ds) decreased with latitude. If all measured values instead of the average value of each population were considered, other significant linear regressions were detected but R2 values were low (0.02–0.08, Figure S1). In this case, head depth (HD), minimum caudal peduncle depth (CPD), caudal fin base length (CFBL), and anal fin length (AFL) decreased with increasing latitude while height of the third dorsal fin branched ray (HDOR3), dorsal fin insertion length (DFIL), and anal fin insertion length (AFIL) showed an opposite relationship.
No significant linear regressions were detected among average values of functional traits and latitude. Also, in this case, accounting for all measured values, significant positive relations were found with relative eye size (RES) and maximum body length (BLmax), while significant negative relations were detected with body lateral shape (BLS) and caudal peduncle throttling (CPT) (Figure S1).
Considering only the 11 populations belonging to B. plebejus species (intra-species variation), the relationships among latitude and NSLL (R2 = 0.72, p < 0.001, average values), NSALL (R2 = 0.56, p = 0.008, average values), NSULL (R2 = 0.07, p < 0.001, all values), HD (R2 = 0.13, p < 0.001, all values), CPD (R2 = 0.55, p = 0.009, average values), CFBL (R2 = 0.37, p = 0.047, average values), AFL (R2 = 0.13, p < 0.001, all values), AFIL (R2 = 0.06, p = 0.002, all values), RES (R2 = 0.55, p = 0.009, average values), and BLS (R2 = 0.12, p < 0.001, all values) were confirmed. Only AFIL and BLS showed opposite relationships between the two datasets.

3.3. Environmental Factors Influencing Barbel Populations

Among the environmental factors, latitude was positively correlated with precipitation (r = 0.71, p = 0.001) and the average flow/basin area (Q/A) ratio (r = 0.47, p = 0.054).
As shown in the RDA (Figure 7), other environmental factors besides latitude explained the morphological variation of barbel populations. The overall variation explained by both RDA axes was 64%.
Among the explanatory variables, RDA axis 1 was mainly negatively associated to Q/A (score −0.56) and elevation (score −0.42), while positively associated to distance from headwater (score 0.41). Among the response variables, this first axis was positively associated to the first component of the PCA on morphometric traits (i.e., longer pectoral, anal, dorsal, and caudal fins; larger head depth and length; longer jaws; and lower position of mouth, eye, and pectoral fin; Figure S2, score 0.94), body lateral shape (BLS, score 0.74), pectoral fin size (PFS, score 0.74), and relative maxillary length (RML, score 0.71), and was negatively associated to NMDS2 (Figure 5b, score −0.88), oral gape position (OGP, score −0.83), pectoral fin vertical position (PFVP, score −0.83), and eye vertical position (EVP, score −0.75).
RDA axis 2 was mainly negatively associated to elevation (score −0.53), to the first axis of the PCA on substrate (i.e., dominance of larger substrate, Figure S2, score −0.52), and to average depth (score −0.48); these are the explanatory variables. It was also negatively associated to the relative eye size (RES, score −0.80), which is the response variable. Moreover, this second axis was positively associated to the second component of the PCA on morphometric traits (i.e., higher length of fin insertions and ventral fin, and higher body depth, Figure S2, score 0.81).
Both the genetic diversity indices were weakly associated to both RDA axes (H scores: 0.44 with F1 and −0.22 with F2; π scores: 0.26 with F1 and −0.32 with F2).

4. Discussion

The primary goal of this study was to investigate the relationships between genetic and morphological variations, detected across 17 Italian barbel populations, and changes in past and present environmental conditions, mainly represented by the latitudinal gradient of approximately 2.5 degrees of the study area. Our previous studies [38,41] have highlighted that Mio-Pleistocenic geomorphological events determined fragmented geographical isolation of the 17 investigated rivers located along the Adriatic side of the Italian Peninsula, leading to barbel speciation. Moreover, a possible cline in meristic traits causing taxonomic identification errors was already pointed out [43,45]. However, the existence of a latitudinal gradient at both morphological and genetic level in Italian barbels was only hypothesized and never specifically tested.
Our results showed that the relative frequency of population-specific haplotypes increased southward, whilst the total number of haplotypes decreased, supporting the isolation of southern barbel populations. These results highlight that the modern-day distribution of freshwater fish throughout multiple rivers reflects the ancient connectivity of their drainages and the subsequent geomorphological history [91], partly supporting our initial hypothesis where a lack of a relationship between genetic variability and local environmental conditions was posited. In contrast, the overall genetic variability expressed by nucleotide and haplotype diversity indices did not offer any clear pattern from the latitudinal point of view. Within each species, we found populations with relative higher haplotype and nucleotide diversity, and populations with intermediate to low levels of one or both genetic diversity indices, independently from the latitudinal positioning of the inhabited river. The RDA analysis revealed some weak correlations between both nucleotide and haplotype diversity indices and some environmental characteristics, some of which are consistent with previous studies, which found correlations between genetic diversity and ecosystem size (i.e., river length in the RDA analysis [32]), habitat features (i.e., refuges and geology in this study [28]), and hydrological regime (represented in this study by the Q/A [92]). Although these results may suggest an influence of local environmental conditions on genetic variability, the RDA scores of both the genetic indices were too low to support a strong relationship between genetic variability and the environmental parameters. Several factors can determine the lack or non-detection of a relationship with local environmental features [28], especially in case of mitochondrial DNA markers. Mitochondrial DNA markers are less sensitive to recent environmental changes than DNA markers of, or those linked to, adaptive loci. This has prevented us to further speculate on the role of environment in shaping the genetic diversity of Italian barbels, which is thus more likely to reflect the past geomorphological evolution of the river catchments, at least at the DNA markers examined.
The main morphological differences previously detected among the four Italian fluviolacustrine barbel species (B. tyberinus, B. plebejus, B. samniticus, and B. fucini) were mainly associated to the number of scales on the lateral line and the shape of the caudal fin; the former trait mainly differentiates B. plebejus from the other three species while the latter mainly differentiate B. tyberinus [38]. Barbus tyberinus was, however, excluded from the present study because of its absence in the Adriatic rivers and because of the difficulties in finding purebred populations. Our study confirmed the previous findings and detected a significant increase of some meristics with latitude. These results agree with the Jordan’s rule [27,65,93,94,95], an ecomorphological rule that postulates an inverse relationship between meristic traits and latitude (particularly temperature) in fish. Additionally, we found that individuals at lower latitudes displayed different coloration, supporting our second hypothesis for which morphological traits are related to current environmental conditions rather than past evolution and species attribution. Skin color variations in teleosts play an important role in life history and behavior, such as camouflage, photoprotection, or antipredation, and are influenced by environmental changes [96,97,98,99]. These variations can be classified into morphological and physiological: the former include variations in skin pigment cell numbers, as well as changes in density, morphology, and distribution [100], and take days or weeks to complete, with a long-term stimulus following comprehensive adaptation to certain conditions such as background pattern, diet, and UV light [101,102,103,104]; the latter are triggered by short-term stimuli that facilitate the dispersion and translocation of skin chromatophores from the periphery to the center [105]. The darker livery of barbels at decreasing latitudes may thus be related to a darker background or to differences in diet or UV light. In southern rivers, this could be also associated to lower flows and water depths, favoring camouflage with the riverbed substrate [106].
The only morphometric trait having a relevant association to latitude was larger caudal fin at lower latitudes. This result was not confirmed by the analysis of the partial dataset which included only B. plebejus populations (intra-species variation) but was supported by negative relationships with latitude of other caudal fin measurements. The other consistent result on morphological traits analysis was the increase of visual acuity (i.e., relative eye size) and the decrease of head size at increasing latitudes. However, considering the whole dataset, an increase of the maximum body length at increasing latitudes was also detected, and this agrees with the temperature-size rule, which postulates that ectotherms living in warmer conditions (i.e., lower latitudes) grow faster as juveniles but attain smaller adult body sizes [107], likely due to growth and reproduction trade-offs across temperatures [108]. The relationship between body size and latitude can be controversial in fish, being different from species to species [109]. Nevertheless, our results were consistent with a previous report dealing with 24 B. barbus populations [110].
Our results thus demonstrate that meristic traits, fish body shape, and some morphometric traits of barbels may be sensitive to climate changes since they have been responsive to an overall limited latitudinal cline. This is an important result given that these traits are not widely studied in this context. Despite several studies assessing the relationship between climate change and phenotypic traits exists (reviewed by Crozier and Hutchings [111]), they focus only on certain “climate-sensitive” traits (e.g., growth, survival, migration, and spawn).
Besides the relationships with latitude, we also found significant relationships between fish morphology and other environmental parameters, which strengthen the potential sensitivity of morphological traits to current climate change. The result of the RDA highlighted an evident role of the ratio between average flow and basin area, elevation, and distance from headwater in shaping fish morphology. Barbel populations were sampled at different elevations (from 28 to 398 m asl) and distance from headwater (from 0.5 to 93 km). Particularly, fin and head size increased at decreasing elevation (in the lower river reaches), while the position of the eyes, mouth, and pectoral fin were lower, indicating a change of the fish position in the water column moving from the headwater to the river outlet. Evidence of the influence of the elevation gradient on morphological traits is widely available in literature for other fish species [112,113].
However, questions remain as to whether the observed morphological variations are a result of phenotypic plasticity (not yet fixed genotypes) or local adaptations (fixed genotypes) and how fast these variations appear across generations. Future experimental studies may shed light on this since understanding these aspects as well as their mechanistic basis is crucial to fully comprehend the vulnerability of fish populations to climate change [114] and other anthropogenic stressors [115].

5. Conclusions

This study proves that the environmental factors, represented by a latitudinal gradient, have an influence on both genetic and morphological traits of three species of freshwater fish belonging to the Barbus genus. While genetic variations are mainly related to latitude due to the historical evolution of the landmass and climatic Mio-Pleistocenic events, morphological variations are mainly affected by changes in current environmental conditions. Particularly, meristic traits such as the number of scales and fish coloration are clearly associated to latitude, and thus should not be used for taxonomical identification in the field. Besides latitude, other environmental factors displayed an important influence on fish morphology, which then appear to be potentially sensitive to climate and human-mediated changes. This should be better understood to unravel the vulnerability or adaptability of barbels to these stressors.
Studies like this one provide important information on the biology of freshwater species, information that is also necessary for their conservation and for the correct management of the rivers they inhabit. The management should take into consideration the current and forecasted alterations of both hydrological and thermal regimes due to the synergistic effects of anthropogenic pressure and climate change.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/w15020325/s1, Table S1: List and brief description of the environmental variables considered in this study. Figure S1: Significant linear regressions of morphological and functional traits with latitude. HD = head depth, CPD = minimum caudal peduncle depth, CFBL = caudal fin base length, HDOR3 = height of the third dorsal fin branched ray, DFIL = dorsal fin insertion length, AFL = anal fin length, AFIL = anal fin insertion length, RES = relative eye size, BLS = body lateral shape, CPT = caudal peduncle throttling, BLmax = maximum body length. Figure S2: Principal component analysis of some environmental characteristics of the investigated areas (i.e., geology and land use of basins, substrate and mesohabitat of sampling sites, and climate), and of morphometric traits (see abbreviations in Figure 2).

Author Contributions

Conceptualization, S.Q., V.D.S., M.L., A.C., S.Z. and I.V.; formal analysis, S.Q. and V.D.S.; investigation, M.L., A.C., S.Q. and V.D.S.; data curation, A.C., S.Q., V.D.S. and I.V.; writing—original draft preparation, S.Q. and V.D.S.; writing—review and editing, S.Q., V.D.S., M.L., A.C., S.Z. and I.V.; visualization, S.Q. and V.D.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The genetic dataset generated during the current study is available in GenBank (ON661665-ON661675 for the D-loop mtDNA) while the morphological and environmental datasets are available from the corresponding authors upon reasonable request.

Acknowledgments

The authors thank the fisheries departments for the local authorization to promote research activities in the field. Furthermore, we are particularly thankful to all the students involved in laboratory work and analyses, to Pietro Volta who provided valuable suggestions on a previous version of the manuscript, and to Alex Ceriani for English editing.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sagnes, P.; Statzner, B. Hydrodynamic abilities of riverine fish: A functional link between morphology and velocity use. Aquat. Living Resour. 2009, 22, 79–91. [Google Scholar] [CrossRef] [Green Version]
  2. Franssen, N.R.; Stewart, L.K.; Schaefer, J.F. Morphological divergence and flow induced phenotypic plasticity in a native fish from anthropogenically altered stream habitats. Ecol. Evol. 2013, 3, 4648–4657. [Google Scholar] [CrossRef] [PubMed]
  3. Samways, K.M.; Leavitt, P.R.; Magnan, P.; Rodríguez, M.A.; Peres-Neto, P.R. Convergent polymorphism between stream and lake habitats: The case of brook char. Can. J. Fishe. Aquat. Sci. 2015, 72, 1406–1414. [Google Scholar] [CrossRef]
  4. Thacker, C.E.; Gkenas, C. Morphometric convergence among European sand gobies in freshwater (Gobiiformes: Gobionellidae). Ecol. Evol. 2019, 9, 8087–8103. [Google Scholar] [CrossRef] [PubMed]
  5. Van Valen, L. Morphological variation and width of ecological niche. Am. Nat. 1965, 99, 377–390. [Google Scholar] [CrossRef]
  6. Bolnick, D.; Svanbäck, R.; Araújo, M.; Persson, L. Comparative support for the niche variation hypothesis that more generalized populations also are more heterogeneous. Proc. Natl. Acad. Sci. USA 2007, 104, 10075–10079. [Google Scholar] [CrossRef] [Green Version]
  7. Snowberg, L.K.; Hendrix, K.M.; Bolnick, D.I. Covarying variances: More morphologically variable populations also exhibit more diet variation. Oecologia 2015, 178, 89–101. [Google Scholar] [CrossRef]
  8. Li, W.; Zhai, D.; Wang, C.; Gao, X.; Liu, H.; Cao, W. Relationships Among Trophic Niche Width, Morphological Variation, and Genetic Diversity of Hemiculter leucisculus in China. Front. Ecol. Evol. 2021, 9, 691218. [Google Scholar] [CrossRef]
  9. Yedier, S.; Bostanci, D. Intra-and Inter-species Discriminations: Scorpaena species from the Black Sea, Aegean Sea, Sea of Marmara and Mediterranean Sea. Sci. Mar. 2021, 85, 197–209. [Google Scholar] [CrossRef]
  10. Scheiner, S.M. Genetics and evolution of phenotypic plasticity. Annu. Rev. Ecol. Syst. 1993, 24, 35–68. [Google Scholar] [CrossRef]
  11. Price, T.D.; Qvarnstrom, A.; Irwin, D.E. The role of phenotypic plasticity in driving genetic evolution. Proc. Biol. Sci. 2003, 270, 1433–1440. [Google Scholar] [CrossRef] [Green Version]
  12. Bower, L.M.; Saenz, D.E.; Winemiller, K.O. Widespread convergence in stream fishes. Biol. J. Linn. Soc. 2021, 133, 863–879. [Google Scholar] [CrossRef]
  13. Kirk, M.A.; Rahel, F.J.; Laughlin, D.C. Environmental filters of freshwater fish community assembly along elevation and latitudinal gradients. Glob. Ecol. Biogeogr. 2022, 31, 470–485. [Google Scholar] [CrossRef]
  14. Bower, L.M.; Winemiller, K.O. Fish assemblage convergence along stream environmental gradients: An intercontinental analysis. Ecography 2019, 42, 1691–1702. [Google Scholar] [CrossRef] [Green Version]
  15. Higham, T.E.; Ferry, L.A.; Schmitz, L.; Irschick, D.J.; Starko, S.; Anderson, P.S.; Bergmann, P.J.; Jamniczky, H.A.; Monteiro, L.R.; Navon, D.; et al. Linking ecomechanical models and functional traits to understand phenotypic diversity. Trends Ecol. Evol. 2021, 36, 860–873. [Google Scholar] [CrossRef]
  16. Viganò, G.; Confortola, G.; Fornaroli, R.; Cabrini, R.; Canobbio, S.; Mezzanotte, V.; Bocchiola, D. Effects of future climate change on a river habitat in an Italian alpine catchment. J. Hydrol. Eng. 2016, 21, 04015063. [Google Scholar] [CrossRef]
  17. Fuso, F.; Stucchi, L.; Bonacina, L.; Fornaroli, R.; Bocchiola, D. Evaluation of water temperature under changing climate and its effect on river habitat in a regulated Alpine catchment. J. Hydrol. 2022, 616, 128816. [Google Scholar] [CrossRef]
  18. Mantua, N.; Tohver, I.; Hamlet, A. Climate change impacts on streamflow extremes and summertime stream temperature and their possible consequences for freshwater salmon habitat in Washington State. Clim. Chang. 2010, 102, 187–223. [Google Scholar] [CrossRef]
  19. Isaak, D.J.; Wollrab, S.; Horan, D.; Chandler, G. Climate change effects on stream and river temperatures across the northwest US from 1980-2009 and implications for salmonid fishes. Clim. Chang. 2012, 113, 499–524. [Google Scholar] [CrossRef] [Green Version]
  20. Schneider, C.; Laizé, C.L.R.; Acreman, M.C.; Florke, M. How will climate change modify river flow regimes in Europe? Hydrol. Earth Syst. Sci. 2013, 17, 325–339. [Google Scholar] [CrossRef]
  21. van Vliet, M.T.; Franssen, W.H.; Yearsley, J.R.; Ludwig, F.; Haddeland, I.; Lettenmaier, D.P.; Kabat, P. Global river discharge and water temperature under climate change. Glob. Environ. Chang. 2013, 23, 450–464. [Google Scholar] [CrossRef]
  22. Ye, F.; Kameyama, S. Long-term nationwide spatiotemporal changes of freshwater temperature in Japan during 1982–2016. J. Environ. Manage. 2021, 281, 111866. [Google Scholar] [CrossRef] [PubMed]
  23. Benejam, L.; Alcaraz, C.; Sasal, P.; Simon-Levert, G.; García-Berthou, E. Life history and parasites of the invasive mosquitofish (Gambusia holbrooki) along a latitudinal gradient. Biol. Invasions 2009, 11, 2265–2277. [Google Scholar] [CrossRef]
  24. Parra, I.; Nicola, G.G.; Vøllestad, L.A.; Elvira, B.; Almodóvar, A. Latitude and altitude differentially shape life history trajectories between the sexes in non-anadromous brown trout. Evol. Ecol. 2014, 28, 707–720. [Google Scholar] [CrossRef]
  25. Riesch, R.; Martin, R.A.; Diamond, S.E.; Jourdan, J.; Plath, M.; Brian Langerhans, R. Thermal regime drives a latitudinal gradient in morphology and life history in a livebearing fish. Biol. J. Linn. Soc. 2018, 125, 126–141. [Google Scholar] [CrossRef] [Green Version]
  26. Rohde, K. Latitudinal gradients in species diversity: The search for the primary cause. Oikos 1992, 65, 514–527. [Google Scholar] [CrossRef] [Green Version]
  27. Barlow, G.W. Causes and significance of morphological variation in fishes. Syst. Zool. 1961, 10, 105–117. [Google Scholar] [CrossRef] [Green Version]
  28. Osborne, M.J.; Perkin, J.S.; Gido, K.B.; Turner, T.F. Comparative riverscape genetics reveals reservoirs of genetic diversity for conservation and restoration of Great Plains fishes. Mol. Ecol. 2014, 23, 5663–5679. [Google Scholar] [CrossRef] [Green Version]
  29. De Santis, V.; Quadroni, S.; Britton, R.J.; Carosi, A.; Gutmann Roberts, C.; Lorenzoni, M.; Crosa, G.; Zaccara, S. Biological and trophic consequences of genetic introgression between endemic and invasive Barbus fishes. Biol. Invasions 2021, 23, 3351–3368. [Google Scholar] [CrossRef]
  30. Prunier, J.G.; Dubut, V.; Loot, G.; Tudesque, L.; Blanchet, S. The relative contribution of river network structure and anthropogenic stressors to spatial patterns of genetic diversity in two freshwater fishes: A multiple-stressors approach. Freshwater Biol. 2018, 63, 6–21. [Google Scholar] [CrossRef]
  31. Radojković, N.; Marinović, Z.; Milošković, A.; Radenković, M.; Đuretanović, S.; Lujić, J.; Simić, V. Effects of stream damming on morphological variability of fish: Case study on large spot barbel Barbus balcanicus. Turk. J. Fish. Aquat. Sci. 2018, 19, 231–239. [Google Scholar] [CrossRef]
  32. Machado, C.B.; Braga-Silva, A.; Freitas, P.D.; Galetti Jr, P.M. Damming shapes genetic patterns and may affect the persistence of freshwater fish populations. Freshwater Biol. 2022, 67, 603–618. [Google Scholar] [CrossRef]
  33. Berbel-Filho, W.M.; Martinez, P.A.; Ramos, T.; Torres, R.A.; Lima, S.M. Inter-and intra-basin phenotypic variation in two riverine cichlids from northeastern Brazil: Potential eco-evolutionary damages of Sao Francisco interbasin water transfer. Hydrobiologia 2016, 766, 43–56. [Google Scholar] [CrossRef]
  34. Schmitt, T.; Fritz, U.; Delfino, M.; Ulrich, W.; Habel, J.C. Biogeography of Italy revisited: Genetic lineages confirm major phylogeographic patterns and a pre-Pleistocene origin of its biota. Front. Zool. 2021, 18, 34. [Google Scholar] [CrossRef]
  35. Gumiero, B.; Maiolini, B.; Comiti, F.; Moroni, F. The Italian Rivers. In Rivers of Europe, 2nd ed.; Tockner, K., Zarfl, C., Robinson, C.T., Eds.; Elsevier: Amsterdam, The Netherlands, 2022; pp. 657–685. [Google Scholar] [CrossRef]
  36. Skoulikidis, N.T.; Sabater, S.; Datry, T.; Morais, M.M.; Buffagni, A.; Dörflinger, G.; Zogaris, S.; Sánchez-Montoya, M.D.M.; Bonada, N.; Kalogianni, E.; et al. Non-perennial Mediterranean rivers in Europe: Status, pressures, and challenges for research and management. Sci. Total Environ. 2017, 577, 1–18. [Google Scholar] [CrossRef]
  37. Darvini, G.; Memmola, F. Assessment of the impact of climate variability and human activities on the runoff in five catchments of the Adriatic Coast of south-central Italy. J. Hydrol. Reg. Stud. 2020, 31, 100712. [Google Scholar] [CrossRef]
  38. Lorenzoni, M.; Carosi, A.; Quadroni, S.; De Santis, V.; Vanetti, I.; Delmastro, G.B.; Zaccara, S. Cryptic diversity within endemic Italian barbels: Revalidation and description of new Barbus species (Teleostei: Cyprinidae). J. Fish Biol. 2021, 98, 1433–1449. [Google Scholar] [CrossRef]
  39. Kottelat, M.; Freyhof, J. Handbook of European Freshwater Fishes; Publications Kottelat: Cornol, Switzerland; Berlin, Germany, 2007. [Google Scholar]
  40. Buonerba, L.; Zaccara, S.; Delmastro, G.B.; Lorenzoni, M.; Salzburger, W.; Gante, H.F. 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–129. [Google Scholar] [CrossRef]
  41. Zaccara, S.; Quadroni, S.; De Santis, V.; Vanetti, I.; Carosi, A.; Britton, R.; Lorenzoni, M. Genetic and morphological analysis reveal a complex biogeographic pattern in the endemic barbel populations of the southern Italian peninsula. Ecol. Evol. 2019, 9, 10185–10197. [Google Scholar] [CrossRef]
  42. Bianco, P.G. Potential role of the palaeohistory of the Mediterranean and Paratethys basins on the early dispersal of Euro-Mediterranean freshwater fishes. Ichthyol. Explor. Fres. 1990, 1, 167–184. [Google Scholar]
  43. Livi, S.; De Innocentiis, S.; Longobardi, A.; Cataudella, S.; Tancioni, L.; Rampacci, M.; Marino, G. Genetic structure of Barbus spp. populations in the Marche Region of central Italy and its relevance to conservation actions. J. Fish Biol. 2013, 82, 806–826. [Google Scholar] [CrossRef] [PubMed]
  44. Geiger, M.F.; Schreiner, C.; Delmastro, G.B.; 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, 1038–1055. [Google Scholar] [CrossRef] [PubMed]
  45. Zaccara, S.; Quadroni, S.; Vanetti, I.; Carosi, A.; La Porta, G.; Crosa, G.; Britton, R.; Lorenzoni, M. Morphologic and genetic variability in the Barbus fishes (Teleostei, Cyprinidae) of Central Italy. Zool. Scr. 2019, 48, 289–301. [Google Scholar] [CrossRef]
  46. Bianco, P.G. A revision of the Italian Barbus species (Cypriniformes: Cyprinidae). Ichthyol. Explor. Fres. 1995, 6, 305–324. [Google Scholar]
  47. Lorenzoni, M.; Carosi, A.; Angeli, V.; Bicchi, A.; Pedicillo, G.; Viali, P. Individuazione e Riconoscimento dei Barbi Autoctoni nel Bacino del Fiume Paglia. Provincia di Terni-Assessorato alla Programmazione Faunistica; Arti grafiche Tiezzi: Terni, Italy, 2006. [Google Scholar]
  48. Meraner, A.; Venturi, A.; Ficetola, G.F.; Rossi, S.; Candiotto, A.; Gandolfi, A. Massive invasion of alien Barbus barbus and introgressive hybridization with endemic Barbus plebejus in Northern Italy: Where, how and why? Mol. Ecol. 2013, 22, 5295–5312. [Google Scholar] [CrossRef]
  49. Carosi, A.; Ghetti, L.; La Porta, G.; Lorenzoni, M. Ecological effects of the European barbel Barbus barbus (L., 1758) (Cyprinidae) invasion on native barbel populations in the Tiber River basin (Italy). Eur. Zool. J. 2017, 84, 420–435. [Google Scholar] [CrossRef] [Green Version]
  50. Piccoli, F.; Burgazzi, G.; Laini, A.; Ferrari, C.; Voccia, A.; Filonzi, L.; Bolpagni, R.; Nonnis Marzano, F. Barbel species arrangement in a regional Natura 2000 network (Emilia Romagna, northern Italy): An altitudinal perspective. J. Limnol. 2017, 76, 140–147. [Google Scholar] [CrossRef] [Green Version]
  51. Zaccara, S.; Quadroni, S.; De Santis, V.; Vanetti, I.; Carosi, A.; Crosa, G.; Britton, R.J.; Lorenzoni, M. Genetic and phenotypic displacement of an endemic Barbus complex by invasive European barbel Barbus barbus in Central Italy. Biol. Invasions 2020, 23, 521–535. [Google Scholar] [CrossRef]
  52. Bagenal, T. Methods for Assessment of Fish Production in Fresh Waters; Blackwell Scientific Publications: Oxford, UK, 1978. [Google Scholar]
  53. Vilizzi, L.; Copp, G.H.; Britton, J.R. Age and growth of European barbel Barbus barbus (Cyprinidae) in the small, mesotrophic River Lee and relative to other populations in England. Knowl. Managt. Aquatic Ecosyst. 2013, 409, 9. [Google Scholar] [CrossRef] [Green Version]
  54. Aljanabi, S.M.; Martinez, I. Universal and rapid salt-extraction of high quality genomic DNA for PCR-based techniques. Nucleic Acids Res. 1997, 25, 4692–4693. [Google Scholar] [CrossRef]
  55. Briolay, J.; Galtier, N.; Brito, R.M.; Bouvet, Y. Molecular phylogeny of Cyprinidae inferred from cytochrome b DNA sequences. Mol. Phylogenet. Evol. 1998, 9, 100–108. [Google Scholar] [CrossRef]
  56. Rossi, G.; Zuffi, G.; Gandolfi, G.; Marchi, A.; Rinaldi, M.; Valli, M. Analisi della Distribuzione delle Specie del Genere Barbus Cuvier, 1871 nei Bacini Idrografici della Regione Abruzzo; Dipartimento di Scienze Biologiche, Geologiche e Ambientali dell’Università di Bologna: Bologna, Italy, 2013. [Google Scholar]
  57. Antognazza, C.M.; Andreou, D.; Zaccara, S.; Britton, R.J. Loss of genetic integrity and biological invasions result from stocking and introductions of Barbus barbus: Insights from rivers in England. Ecol. Evol. 2016, 6, 1280–1292. [Google Scholar] [CrossRef]
  58. Thompson, J.D.; Higgins, D.G.; Gibson, T.J. CLUSTAL W: Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22, 4673–4680. [Google Scholar] [CrossRef] [Green Version]
  59. Hall, T.A. BioEdit: A user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp. Ser. 1999, 41, 95–98. [Google Scholar]
  60. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2019; Available online: https://www.R-project.org/ (accessed on 6 May 2021).
  61. Adams, D.C.; Collyer, M.L.; Kaliontzopoulou, A.; Geomorph: Software for Geometric Morphometric Analysis. R Package Version 3.0.6. 2018. Available online: https://cran.r-project.org/package=geomorph (accessed on 6 May 2021).
  62. Beacham, T.D. Meristic and morphometric variation in pink salmon (Oncorhynchus gorbuscha) in southern British Columbia and Puget Sound. Can. J. Zool. 1985, 63, 366–372. [Google Scholar] [CrossRef] [Green Version]
  63. Xiang, T.; Dong, X.; Grenouillet, G. Ecological and biological traits of non-native freshwater fish species differentiate them from native species in China. Ecol. Indic. 2021, 131, 108218. [Google Scholar] [CrossRef]
  64. Olden, J.D.; Jackson, D.A.; Peres-Neto, P.R. Predictive models of fish species distributions: A note on proper validation and chance predictions. T. Am. Fish. Soc. 2002, 131, 329–336. [Google Scholar] [CrossRef]
  65. Ha, L.M.; Iguchi, K.I. Geographical continuity and discontinuity in the meristic characteristics of ayus of the southern subspecies Plecoglossus altivelis ryukyuensis. Ichthyol. Res. 2021, 68, 177–181. [Google Scholar] [CrossRef]
  66. Toussaint, A.; Charpin, N.; Brosse, S.; Villéger, S. Global functional diversity of freshwater fish is concentrated in the Neotropics while functional vulnerability is widespread. Sci. Rep. 2016, 6, 22125. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Reecht, Y.; Rochet, M.-J.; Trenkel, V.M.; Jennings, S.; Pinnegar, J.K. Use of morphological characteristics to define functional groups of predatory fishes in the Celtic Sea. J. Fish Biol. 2013, 83, 355–377. [Google Scholar] [CrossRef] [Green Version]
  68. Winemiller, K.O. Ecomorphological diversification in lowland freshwater fish assemblages from five biotic regions. Ecol. Monogr. 1991, 61, 343–365. [Google Scholar] [CrossRef]
  69. Boyle, K.S.; Horn, M.H. Comparison of feeding guild structure and ecomorphology of intertidal fish assemblages from central California and central Chile. Mar. Ecol. Prog. Ser. 2006, 319, 65–84. [Google Scholar] [CrossRef]
  70. Dumay, O.; Tari, P.S.; Tomasini, J.A.; Mouillot, D. Functional groups of lagoon fish species in Languedoc Roussillon, southern France. J. Fish Biol. 2004, 64, 970–983. [Google Scholar] [CrossRef]
  71. Lefcheck, J.S.; Buchheister, A.; Laumann, K.M.; Stratton, M.A.; Sobocinski, K.L.; Chak, S.T.C.; Clardy, T.R.; Reynolds, P.L.; Latour, R.J.; Duffy, J.E. Dimensions of biodiversity in Chesapeake Bay demersal fishes: Patterns and drivers through space and time. Ecosphere 2014, 5, art14. [Google Scholar] [CrossRef]
  72. Fulton, C.J.; Bellwood, D.R.; Wainwright, P.C. The relationship between swimming ability and habitat use in wrasses (Labridae). Mar. Biol. 2001, 139, 25–33. [Google Scholar] [CrossRef]
  73. Webb, P.W. Body form, locomotion and foraging in aquatic vertebrates. Am. Zool. 1984, 24, 107–120. [Google Scholar] [CrossRef]
  74. QGIS Development Team. QGIS Geographic Information System. QGIS Association. 2022. Available online: http://www.qgis.or (accessed on 4 June 2021).
  75. Available online: http://wms.pcn.minambiente.it/ogc?map=/ms_ogc/WMS_v1.3/Vettoriali/Carta_geolitologica.map (accessed on 4 June 2021).
  76. Available online: http://wms.pcn.minambiente.it/ogc?map=/ms_ogc/WMS_v1.3/Vettoriali/Corine_Land_Cover2012.map (accessed on 4 June 2021).
  77. Available online: http://www.scia.isprambiente.it (accessed on 4 October 2021).
  78. DECRETO 8 Novembre 2010, n. 260. Regolamento Recante i Criteri Tecnici per la Classificazione dello Stato dei Corpi Idrici Superficiali, per la Modifica delle Norme Tecniche del Decreto Legislativo 3 Aprile 2006, n. 152, Recante Norme in Materia Ambientale, Predisposto ai Sensi dell’articolo 75, Comma 3, del Medesimo Decreto Legislativo. (11G0035) (GU Serie Generale n.30 del 07-02-2011-Suppl. Ordinario n. 31). Available online: https://www.gazzettaufficiale.it/eli/id/2011/02/07/011G0035/sg (accessed on 9 July 2021).
  79. Buffagni, A.; Erba, S. Intercalibrazione e classificazione di qualità ecologica dei fiumi per la 2000/60/EC (WFD): l’indice STAR-ICMi. IRSA-CNR Not. Dei Metod. Anal. 2007, 1, 94–100. [Google Scholar]
  80. Available online: https://www.arpa.marche.it/images/PUBBLICAZIONI/BACINI_2015-2017/Relazione%20fiumi%202015-2017_senza_accorp_lug2019_revfeb2020.pdf (accessed on 6 October 2021).
  81. Available online: https://www.artaabruzzo.it/acque-superficiali.php?id_page=1 (accessed on 6 October 2021).
  82. Available online: https://www.arpamolise.it/Acque/Super/rapporto_2018.pdf (accessed on 6 October 2021).
  83. Available online: https://www.arpa.puglia.it/pagina3381_stato-ecologico-delle-acque-superficiali-interne.html (accessed on 6 October 2021).
  84. Available online:. Available online: https://www.regione.marche.it/Regione-Utile/Ambiente/Tutela-delle-acque/PTA#Documentazione (accessed on 8 October 2021).
  85. Available online: https://www.regione.abruzzo.it/content/piano-tutela-delle-acque (accessed on 8 October 2021).
  86. Available online: https://www.regione.molise.it/flex/cm/pages/ServeBLOB.php/L/IT/IDPagina/13780 (accessed on 8 October 2021).
  87. Available online: http://www.adb.basilicata.it/adb/pStralcio/distretto_bas.asp (accessed on 8 October 2021).
  88. Rozas, J.; Ferrer-Mata, A.; Sanchez-DelBarrio, J.C.; Guirao-Rico, S.; Librado, P.; Ramos-Onsins, S.E.; Sánchez-Gracia, A. DnaSP 6: DNA sequence polymorphism analysis of large data sets. Mol. Biol. Evol. 2017, 34, 3299–3302. [Google Scholar] [CrossRef]
  89. Leigh, J.W.; Bryant, D. Data from: PopART: Full-feature software for haplotype network construction. Methods Ecol. Evol. 2015, 6, 1110–1116. [Google Scholar] [CrossRef]
  90. Klingenberg, C.P. MorphoJ: An integrated software package for geometric morphometrics. Mol. Ecol. Resour. 2011, 11, 353–357. [Google Scholar] [CrossRef]
  91. Won, H.; Jeon, H.B.; Suk, H.Y. Evidence of an ancient connectivity and biogeodispersal of a bitterling species, Rhodeus notatus, across the Korean Peninsula. Sci. Rep. 2020, 10, 1011. [Google Scholar] [CrossRef] [Green Version]
  92. Gido, K.B.; Whitney, J.E.; Perkin, J.S.; Turner, T.F. Fragmentation, connectivity and fish species persistence in freshwater ecosystems. In Conservation of Freshwater Fishes; Closs, G.P., Krkosek, M., Olden, J.D., Eds.; Cambridge University Press: Cambridge, UK, 2016; pp. 292–323. [Google Scholar]
  93. Jordan, D.S. The Days of a Man: Being Memories of a Naturalist, Teacher, and Minor Prophet of Democracy (Vol. 2); World Book Company: Chicago, IL, USA, 1922. [Google Scholar]
  94. Hubbs, C.L. Variations in the number of vertebrae and other meristic characters of fishes correlated with the temperature of water during development. Am. Nat. 1922, 56, 360–372. [Google Scholar] [CrossRef] [Green Version]
  95. McDowall, R. Jordan’s and other ecogeographical rules, and the vertebral number in fishes. J. Biogeogr. 2008, 35, 501–508. [Google Scholar] [CrossRef]
  96. Cal, L.; Suarez-Bregua, P.; Cerdá-Reverter, J.M.; Braasch, I.; Rotllant, J. Fish pigmentation and the melanocortin system. Comp. Biochem. Phys. A 2017, 211, 26–33. [Google Scholar] [CrossRef] [Green Version]
  97. Nüsslein-Volhard, C.; Singh, A.P. How fish color their skin: A paradigm for development and evolution of adult patterns: Multipotency, plasticity, and cell competition regulate proliferation and spreading of pigment cells in zebrafish coloration. BioEssays 2017, 39, 1600231. [Google Scholar] [CrossRef] [Green Version]
  98. Mizusawa, K.; Yamamura, Y.; Kasagi, S.; Cerdá-Reverter, J.M.; Takahashi, A. Expression of genes for melanotropic peptides and their receptors for morphological color change in goldfish Carassius auratus. Gen. Comp. Endocr. 2018, 264, 138–150. [Google Scholar] [CrossRef]
  99. Luo, M.; Lu, G.; Yin, H.; Wang, L.; Atuganile, M.; Dong, Z. Fish pigmentation and coloration: Molecular mechanisms and aquaculture perspectives. Rev. Aquacult. 2021, 13, 2395–2412. [Google Scholar] [CrossRef]
  100. Chavin, W. Fundamental aspects of morphological melanin color changes in vertebrate skin. Am. Zool. 1969, 9, 505–520. [Google Scholar] [CrossRef]
  101. Sugimoto, M. Morphological color changes in fish: Regulation of pigment cell density and morphology. Microsc. Res. Techniq. 2002, 58, 496–503. [Google Scholar] [CrossRef] [PubMed]
  102. Henning, F.; Jones, J.C.; Franchini, P.; Meyer, A. Transcriptomics of morphological color change in polychromatic Midas cichlids. BMC Genom. 2013, 14, 171. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  103. Rashidian, G.; Rainis, S.; Prokić, M.D.; Faggio, C. Effects of different levels of carotenoids and light sources on swordtail fish (Xiphophorus helleri) growth, survival rate and reproductive parameters. Nat. Prod. Res. 2020, 35, 3675–3686. [Google Scholar] [CrossRef] [PubMed]
  104. Wang, L.; Luo, M.; Yin, H.; Zhu, W.; Fu, J.; Dong, Z. Effects of background adaptation on the skin color of Malaysian red tilapia. Aquaculture 2020, 521, 735061. [Google Scholar] [CrossRef]
  105. Aspengren, S.; Sköld, H.; Wallin, M. Different strategies for color change. Cell. Mol. Life Sci. 2009, 66, 187–191. [Google Scholar] [CrossRef]
  106. Donnelly, W.A.; Dill, L.M. Evidence for crypsis in coho salmon, Oncorhynchus kisutch (Walbaum), parr: Substrate colour preference and achromatic reflectance. J. Fish Biol. 1984, 25, 183–195. [Google Scholar] [CrossRef]
  107. Atkinson, D. Temperature and organism size –a biological law for ectotherms. Adv. Ecol. Res. 1994, 25, 1–58. [Google Scholar] [CrossRef]
  108. Wootton, H.F.; Morrongiello, J.R.; Schmitt, T.; Audzijonyte, A. Smaller adult fish size in warmer water is not explained by elevated metabolism. Ecol. Lett. 2022, 25, 1177–1188. [Google Scholar] [CrossRef]
  109. Rypel, A.L. The cold-water connection: Bergmann’s rule in north American freshwater fishes. Am. Nat. 2014, 183, 147–156. [Google Scholar] [CrossRef]
  110. Blanck, A.; Lamouroux, N. Large-scale intraspecific variation in life-history traits of European freshwater fish. J. Biogeogr. 2007, 34, 862–875. [Google Scholar] [CrossRef]
  111. Crozier, L.G.; Hutchings, J.A. Plastic and evolutionary responses to climate change in fish. Evol. Appl. 2014, 7, 68–87. [Google Scholar] [CrossRef]
  112. Rahel, F.J.; Hubert, W.A. Fish assemblages and habitat gradients in a Rocky Mountain-Great Plains stream: Biotic zonation and additive patterns of community change. Trans. Am. Fish. Soc. 1991, 120, 319–332. [Google Scholar] [CrossRef]
  113. Olinger, C.T.; Peoples, B.K.; Frimpong, E.A. Reproductive life history of Heterandria bimaculata (Heckel, 1848) (Poeciliinae: Poeciliidae) in the Honduran interior highlands: Trait variation along an elevational gradient. Neotrop. Ichthyol. 2016, 14, e150050. [Google Scholar] [CrossRef] [Green Version]
  114. Michel, M.J.; Chien, H.; Beachum, C.E.; Bennett, M.G.; Knouft, J.H. Climate change, hydrology, and fish morphology: Predictions using phenotype-environment associations. Climatic Chang. 2017, 140, 563–576. [Google Scholar] [CrossRef]
  115. Rudolfsen, T.; Watkinson, D.A.; Poesch, M. Morphological divergence of the threatened Rocky Mountain sculpin (Cottus sp.) is driven by biogeography and flow regime: Implications for mitigating altered flow regime to freshwater fishes. Aquat. Conserv. 2018, 28, 78–86. [Google Scholar] [CrossRef]
Figure 1. Map of Italy with the new sampling sites (red circles) and the locations of the other purebred barbel populations considered in the present study (black circles). The dashed red line indicates the watershed.
Figure 1. Map of Italy with the new sampling sites (red circles) and the locations of the other purebred barbel populations considered in the present study (black circles). The dashed red line indicates the watershed.
Water 15 00325 g001
Figure 2. (a) Position of the 33 landmarks used for body shape analysis: (1) anterior tip of snout, (2, 3) anterior and posterior end of the eye, (4) orthogonal projection on the dorsal profile of the eye center, (5) lateral projection of the eye center on the insertion of the operculum, (6) intersection of the operculum at the lateral profile, (7, 8) ventral and dorsal end of gills, (9) anterior insertion of pectoral fin, (10, 11) orthogonal projections on the dorsal and ventral profile of the anterior insertion of pectoral fin, (12,13) anterior and posterior insertion of dorsal fin, (14) insertion of pelvic fin, (15, 16) posterior and anterior insertion of anal fin, (17, 18) anterior attachment of dorsal and ventral membrane of caudal fin, (19) base of middle caudal rays, (20, 21) orthogonal projections on the dorsal and ventral profile of the base of middle caudal rays, (22) fork, (23, 24) orthogonal projections on the dorsal and ventral profile of fork, (25, 26) end of the upper and lower lobe of caudal fin, (27, 28) lateral projections of anterior attachment of dorsal and ventral membrane of caudal fin, (29) opening of mouth, (30) posterior end of jaw, (31) posterior-most edge of opercle, (32) supraoccipital, posteromedial tip (visible as an indentation in dorsal surface), (33) posterior nare, posterior margin. (b) The 25 morphometric traits (BLmax, maximum body length; BL, body (standard) length; MOD, mouth-operculum distance; POD, pre-orbital distance; HD, head depth; ED, eye diameter; EH, eye position; MO, oral gape position; MJL, maxillary jaw length; LLL, lower jaw length; BD, maximum body depth (i.e., maximum height); HDOR1, height of the dorsal fin ossified ray; HDOR3, height of the third dorsal fin branched ray; DFIL, dorsal fin insertion length; AFL, anal fin length; AFIL, anal fin insertion length; PFL, pectoral fin length; PFI, pectoral fin position; VFL, ventral fin length; CPD, minimum caudal peduncle depth (i.e., minimum height); CFD, maximum caudal fin depth; CFBL, caudal fin base length; FD, fork depth; LLL, length of the lower lobe of caudal fin; ULL, length of the upper lobe of caudal fin) and 4 meristic traits (NDBR, the number of dorsal fin branched rays; NSLL, the number of scales on the lateral line, and on rows above–NSALL–and under–NSULL–the lateral line) considered for morphological analysis.
Figure 2. (a) Position of the 33 landmarks used for body shape analysis: (1) anterior tip of snout, (2, 3) anterior and posterior end of the eye, (4) orthogonal projection on the dorsal profile of the eye center, (5) lateral projection of the eye center on the insertion of the operculum, (6) intersection of the operculum at the lateral profile, (7, 8) ventral and dorsal end of gills, (9) anterior insertion of pectoral fin, (10, 11) orthogonal projections on the dorsal and ventral profile of the anterior insertion of pectoral fin, (12,13) anterior and posterior insertion of dorsal fin, (14) insertion of pelvic fin, (15, 16) posterior and anterior insertion of anal fin, (17, 18) anterior attachment of dorsal and ventral membrane of caudal fin, (19) base of middle caudal rays, (20, 21) orthogonal projections on the dorsal and ventral profile of the base of middle caudal rays, (22) fork, (23, 24) orthogonal projections on the dorsal and ventral profile of fork, (25, 26) end of the upper and lower lobe of caudal fin, (27, 28) lateral projections of anterior attachment of dorsal and ventral membrane of caudal fin, (29) opening of mouth, (30) posterior end of jaw, (31) posterior-most edge of opercle, (32) supraoccipital, posteromedial tip (visible as an indentation in dorsal surface), (33) posterior nare, posterior margin. (b) The 25 morphometric traits (BLmax, maximum body length; BL, body (standard) length; MOD, mouth-operculum distance; POD, pre-orbital distance; HD, head depth; ED, eye diameter; EH, eye position; MO, oral gape position; MJL, maxillary jaw length; LLL, lower jaw length; BD, maximum body depth (i.e., maximum height); HDOR1, height of the dorsal fin ossified ray; HDOR3, height of the third dorsal fin branched ray; DFIL, dorsal fin insertion length; AFL, anal fin length; AFIL, anal fin insertion length; PFL, pectoral fin length; PFI, pectoral fin position; VFL, ventral fin length; CPD, minimum caudal peduncle depth (i.e., minimum height); CFD, maximum caudal fin depth; CFBL, caudal fin base length; FD, fork depth; LLL, length of the lower lobe of caudal fin; ULL, length of the upper lobe of caudal fin) and 4 meristic traits (NDBR, the number of dorsal fin branched rays; NSLL, the number of scales on the lateral line, and on rows above–NSALL–and under–NSULL–the lateral line) considered for morphological analysis.
Water 15 00325 g002
Figure 3. D-loop median-joining network. Circle size is proportional to the observed haplotype frequencies and black points and dashes represent unobserved haplotypes and potential intermediates. Dotted circles show species boundaries (B. plebejus, B. samniticus, and B. fucini).
Figure 3. D-loop median-joining network. Circle size is proportional to the observed haplotype frequencies and black points and dashes represent unobserved haplotypes and potential intermediates. Dotted circles show species boundaries (B. plebejus, B. samniticus, and B. fucini).
Water 15 00325 g003
Figure 4. Linear regressions among latitude and the total number of haplotypes (n haplotypes) or the relative frequency of population-specific haplotypes (PSh) detected in each barbel population.
Figure 4. Linear regressions among latitude and the total number of haplotypes (n haplotypes) or the relative frequency of population-specific haplotypes (PSh) detected in each barbel population.
Water 15 00325 g004
Figure 5. (a) Output of the canonical variate analysis (CVA) on fish body shape carried out between the 17 barbel populations (see Table 1); 90% confidence ellipses were drawn for each population. Wireframe graphs indicate the shape changes along each axis (from grey to dashed black line). (b) Non-parametric multidimensional scaling (NMDS) showing the Mahalanobis distances (MDs) among the populations detected through the CVA.
Figure 5. (a) Output of the canonical variate analysis (CVA) on fish body shape carried out between the 17 barbel populations (see Table 1); 90% confidence ellipses were drawn for each population. Wireframe graphs indicate the shape changes along each axis (from grey to dashed black line). (b) Non-parametric multidimensional scaling (NMDS) showing the Mahalanobis distances (MDs) among the populations detected through the CVA.
Water 15 00325 g005aWater 15 00325 g005b
Figure 6. Significant linear regressions among morphological traits (average values) and latitude. CFD = maximum caudal fin depth, Ccf = color of caudal fin (calculated as relative frequency of orange color), Ds = presence of dots on scales (calculated as relative frequency of individuals having dots), NSLL = the number of scales on the lateral line, NSALL = the number of scales on rows above the lateral line, NSULL = the number of scales on rows under the lateral line.
Figure 6. Significant linear regressions among morphological traits (average values) and latitude. CFD = maximum caudal fin depth, Ccf = color of caudal fin (calculated as relative frequency of orange color), Ds = presence of dots on scales (calculated as relative frequency of individuals having dots), NSLL = the number of scales on the lateral line, NSALL = the number of scales on rows above the lateral line, NSULL = the number of scales on rows under the lateral line.
Water 15 00325 g006
Figure 7. Redundancy analysis (RDA) performed using environmental characteristics as explanatory variables (red circles, for some variables, the first axis of the PCA shown in Figure S2 was used) and morphometric (NMDS axes in Figure 4 and PCA axes in Figure S2), functional (see abbreviations in Table 2) and genetic (H = haplotype diversity, π = nucleotide diversity) traits as response variables (green–morphological or blue–genetic–squares). Different sites/populations are indicated by black circles.
Figure 7. Redundancy analysis (RDA) performed using environmental characteristics as explanatory variables (red circles, for some variables, the first axis of the PCA shown in Figure S2 was used) and morphometric (NMDS axes in Figure 4 and PCA axes in Figure S2), functional (see abbreviations in Table 2) and genetic (H = haplotype diversity, π = nucleotide diversity) traits as response variables (green–morphological or blue–genetic–squares). Different sites/populations are indicated by black circles.
Water 15 00325 g007
Table 1. Main characteristics of the sampling sites, and genetic attribution of barbel populations. ID = site/population identification code, Ng = number of samples used for genetic (D-loop) analysis, Nm = number of samples used for morphological analysis.
Table 1. Main characteristics of the sampling sites, and genetic attribution of barbel populations. ID = site/population identification code, Ng = number of samples used for genetic (D-loop) analysis, Nm = number of samples used for morphological analysis.
WatercourseBasinLatitude (N)Longitude (E)Elevation (m)Distance from headwater (km)SpeciesIDNgNmReferences
BorelloSavio44°3′8.1″12°10′47″7730No purebred---This study
MarecchiaMarecchia43°49′35.2″12°15′2.4″33527No purebred---This study
FogliaFoglia43°45′1.2″12°27′13.3″30025B. plebejusS11323This study
BossoMetauro43°30′49.8″12°33′42.8″39810B. plebejusS24310[45]
BuranoMetauro43°28′ 48.3″12°37′16.4″34114B. plebejusS31410[45]
EsinoEsino43°24′24.4″12°58′47″23719B. plebejusS42610[45]
PotenzaPotenza43°16′59.8″13°20′5.1″15045B. plebejusS52510[45]
ChientiChienti43°9′ 35.4″13°14′37″26840B. plebejusS63410[45]
TennaTenna43°1′10.6″13° 25′ 34″36314B. plebejusS71910[45]
AsoAso43°1′25.5″13°37′20.1″18538B. plebejusS82710[45]
TrontoTronto42°46′12.2″13°24′29.4″37039B. plebejusS92810[45]
SalinelloSalinello42°47′17″13°53′28.2″2839B. plebejusS102129This study
TordinoTordino42°42′8″13°49′49.9″9744B. plebejusS111931This study
VomanoVomano42°35′21.6″13°40′51.4″20333B. samniticusS122510[45]
GiardinoAterno-Pescara42°10′25.8″13°49′51.4″2470.5B. samniticusS132410[41]
SangroSangro42°3′10.7″14°20′51.3″11693B. samniticusS142310[41]
TrignoTrigno41° 54′ 14.2″14°38′43.7″12161B. samniticusS151823This study
TappinoFortore41°33′13.2″14°52′33.9″23930B. fuciniS162710[41]
OfantoOfanto41°4′26.1″15°32′46.4″21974B. fuciniS172010[41]
Table 2. List of ten functional traits associated with food acquisition and locomotion (slightly modified from Xiang et al. [63], see Figure 2b).
Table 2. List of ten functional traits associated with food acquisition and locomotion (slightly modified from Xiang et al. [63], see Figure 2b).
Functional TraitCodeFormulaPotential Link with Fish FunctionsReferences
Maximum body lengthBLmaxBLmaxSize is linked to metabolism, trophic impacts, locomotion ability, nutrient cycling[66]
Body elongationBEBL/BDHydrodynamism[67]
Eye vertical positionEVPEH/BDPosition of fish and/or of its prey in the water column[68]
Relative eye sizeRESED/HDVisual acuity[69]
Oral gape positionOGPMO/BDFeeding position in the water column[70,71]
Relative maxillary lengthRMLMJL/HDSize of mouth and strength of jaw[66]
Body lateral shapeBLSHD/BDHydrodynamism and head size[66]
Pectoral fin vertical positionPFVPPFI/BDPectoral fin use for swimming[70]
Pectoral fin sizePFSPFL/BLPectoral fin use for swimming[72]
Caudal peduncle throttlingCPTCFD/CPDCaudal propulsion efficiency through reduction of drag[73]
Table 3. Total number of haplotyopes (n), relative frequency of population-specific haplotypes (PSh), number of polymorphic sites (S), and haplotype (H) and nucleotide (π) diversity of the 17 barbel populations investigated along the Adriatic side.
Table 3. Total number of haplotyopes (n), relative frequency of population-specific haplotypes (PSh), number of polymorphic sites (S), and haplotype (H) and nucleotide (π) diversity of the 17 barbel populations investigated along the Adriatic side.
SpeciesWatercourseID nPShSπH
B. plebejusFogliaS1 70.380.002 ± 0.0010.79 ± 0.11
BossoS2 70.460.002 ± 0.0010.87 ± 0.06
BuranoS3 40.350.002 ± 0.0010.68 ± 0.07
EsinoS4 110.5110.002 ± 0.0010.85 ± 0.05
PotenzaS5 40.030.0006 ± 0.00060.46 ± 0.10
ChientiS6 40.340.001 ± 0.0010.6 ± 0.07
TennaS7 21.010.0001 ± 0.00020.10 ± 0.09
AsoS8 31.040.001 ± 0.0010.55 ± 0.08
TrontoS9 40.540.001 ± 0.0010.54 ± 0.09
SalinelloS10 40.560.001 ± 0.0010.61 ± 0.09
TordinoS11 60.550.001 ± 0.0010.70 ± 0.10
B. samniticusVomanoS12 51.060.0005 ± 0.00050.30 ± 0.12
GiardinoS13 30.730.0003 ± 0.00040.16 ± 0.10
SangroS14 20.530.002 ± 0.0010.52 ± 0.03
TrignoS15 50.840.001 ± 0.0010.76 ± 0.06
B. fuciniTappinoS16 31.040.001 ± 0.00080.32 ± 0.10
OfantoS17 31.030.0007 ± 0.00060.43 ± 0.12
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Quadroni, S.; De Santis, V.; Carosi, A.; Vanetti, I.; Zaccara, S.; Lorenzoni, M. Past and Present Environmental Factors Differentially Influence Genetic and Morphological Traits of Italian Barbels (Pisces: Cyprinidae). Water 2023, 15, 325. https://doi.org/10.3390/w15020325

AMA Style

Quadroni S, De Santis V, Carosi A, Vanetti I, Zaccara S, Lorenzoni M. Past and Present Environmental Factors Differentially Influence Genetic and Morphological Traits of Italian Barbels (Pisces: Cyprinidae). Water. 2023; 15(2):325. https://doi.org/10.3390/w15020325

Chicago/Turabian Style

Quadroni, Silvia, Vanessa De Santis, Antonella Carosi, Isabella Vanetti, Serena Zaccara, and Massimo Lorenzoni. 2023. "Past and Present Environmental Factors Differentially Influence Genetic and Morphological Traits of Italian Barbels (Pisces: Cyprinidae)" Water 15, no. 2: 325. https://doi.org/10.3390/w15020325

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop