Abstract
Exotic species invasions often result in native biodiversity loss, i.e. a lower taxonomic diversity, but current knowledge on invasions effects underlined a potential increase of functional diversity. We thus explored the connections between functional diversity and exotic species invasions, while accounting for their environmental drivers, using a fine-resolution large dataset of Mediterranean stream fish communities. While functional diversity of native and exotic species responded similarly to most environmental constraints, we found significant differences in the effects of altitude and in the different ranking of constraints. These differences suggest that invasion dynamics could play a role in overriding some major environmental drivers. Our results also showed that a lower diversity of ecological traits in communities (about half of less disturbed communities) corresponded to a high invasion degree, and that the exotic component of communities had typically less diverse ecological traits than the native one, even when accounting for stream order and species richness. Overall, our results suggest that possible outcomes of severe exotic species invasions could include a reduced functional diversity of invaded communities, but analyzing data with finer ecological, temporal and spatial resolutions would be needed to pinpoint the causal relationship between invasions and functional diversity.
Similar content being viewed by others
Introduction
Biodiversity is a key element of ecosystem functioning and characterizes its resilience to different pressures1, but suffers from a general worldwide decline2. Among the causes of global biodiversity loss, exotic invasions are often placed at the top3,4, but mechanisms and impact of species invasions might vary in different ecosystems, taxa and spatial scales. While a lot of attention has been devoted to the consequences of exotic invasions on species taxonomic diversity at different geographical scales5, hardly any regional extinctions have been recorded in aquatic taxa such as fish (with some notable exceptions, e.g.6,7), but extirpation, species substitution and decrease of native biomass have all been reported as a local result of exotic invasions8,9,10,11. The detection of exotic invasions‘ effects could thus depend on the spatial and biological detail level of the data used, so that some effects could be potentially overlooked or misinterpreted if the data is insufficiently detailed.
For example, fish introductions could enrich functional diversity at the regional scale, because, from a presence/absence perspective, they can increase regional species richness and consequently increase the number of functional traits present in that area12. Functional diversity describes the distinctive assemblage of morphological, biochemical, physiological, structural, phenological or behavioral traits that characterizes living communities, and that is highly coupled to environmental conditions13. Ecofunctional diversity is a subset of functional diversity, focusing on the combination of ecological traits in communities14, but is referred to simply as functional diversity hereafter, because it is a more common term. Ecological traits in a community are usually selected by factors such as habitat diversity, geography, land use or water chemistry15, but established exotic species are often generalists16 and invasion dynamics (e.g. human-aided dispersion) could partly override habitat selectivity.
Competition with introduced exotic species could act as an additional filter for native species, sometimes stronger than environmental gradients11. Furthermore, widespread and severe biological invasions could result in the taxonomic homogenization of invaded communities17,18,19, with few exotic species dominating heavily invaded areas, a pattern which has been detected at least in plants20. Most likely, taxonomic homogenization would also affect functional diversity of invaded areas, as it leads to communities with a lower number of species and a lower variety of traits, but it remains unclear whether this is a significant element of exotic invasions.
Exotic species are a main cause of the loss of biological diversity in the Mediterranean region21,22, particularly in freshwater habitats23. The functional structure of Mediterranean freshwater fish communities is relatively uncharted, and an ecological trait characterization of native and exotic fish species has only been recently defined for some areas24,25,26. This finally enables further research on the linkages between biological invasions, functional diversity and the environment.
We used a spatially-broad yet very detailed dataset, comprising several river basins in northern Italy at a late invasion stage (i.e. >30 years after major invasions), as a test case to explore the outcomes of exotic species invasion in freshwater fish communities and to investigate the relationships between environmental factors, invasions and functional diversity. We used boosted regression tree (BRT) analysis to test the hypothesis (H1) that functional diversity of exotic and native species would respond differently to environmental variables, as invasion dynamics could temporarily override habitat selectivity. We then used spatial and regression analyses to test whether (H2) the overall functional diversity of communities would be negatively or positively affected by different degrees of exotic invasions. Our results would ultimately reveal the connections between functional diversity and exotic species invasions, on the background of habitat filtering.
Results
BRT analysis showed that environmental variables typically had different magnitude and direction of relative influence on exotic and native species functional diversity (Fig. 1). While altitude was one of the most significant variables negatively affecting native species functional diversity (and positively that of exotic ones), high temperature, low salinity or high turbidity were clearly linked to higher exotic species functional diversity (Fig. 1). Forest cover was associated to lower functional diversity of both native and exotic species, while brackish water was associated to a higher functional diversity of native species (Fig. 1).
Several sites showed minimum levels of exotic invasion (invasion degree ≤ 10% for 140 sites, 41.9% of the total) and most of these sites hosted completely native communities (126 sites, 37.6% of the total). However, the majority of sites were invaded (209 sites, 62.4% of the total) and the invasion degree was relatively severe in most of them (invasion degree ≥ 50% for 134 sites, 40% of the total), including some sites where the community exclusively comprised exotic species (10 sites, 3% of the total). Our spatial analysis underlined that the most severely invaded sites were mostly located on the lower stretches of most watercourses examined (Fig. 2a). Moreover, in invaded communities, there was a clear spatial overlap between the most invaded areas and the areas where functional diversity was lowest (Fig. 2b).
In invaded sites, functional diversity was also clearly negatively linked to the invasion degree (Fig. 3). At the highest invasion degrees, a decrease of nearly 50% of functional diversity, compared with less disturbed communities, could be observed.
Low-order streams (in the uplands) showed lower values of species richness and functional diversity than high-order streams (in the lowlands) (Fig. 4a,b). Species richness was low in upland streams, but these habitats showed a larger variation in functional diversity than lowland streams. Canals, in the lowlands, had lower richness and functional diversity values than natural rivers in the same area (Fig. 4a,b). The difference between native and exotic functional diversity was consistently lower than what suggested by their relative richness (Fig. 4c,d), with native species usually showing higher values of functional diversity relative to their richness. Exotic species showed lower functional diversity values than native species in all natural rivers, which was most evident in higher order streams, but the pattern was reversed in canals (Fig. 4d).
Both native and exotic functional diversity in communities were linked to species richness and tended to saturate with growing richness (Fig. 5), albeit with different rates and asymptotes. However, accounting for species richness, community functional diversity for exotic species was generally lower than that of native ones (Fig. 5).
Discussion
Functional diversity of native and exotic species, as expressed by the FDis metric, responded similarly to most environmental constraints. However, confirming our first hypothesis (H1), we found significant differences in the effects of altitude (positive effect for native and negative effect for exotic species) and in the different ranking of constraints. This suggests that invasion dynamics could play a role in overriding some major environmental drivers of functional diversity. Our results also showed that a high invasion degree corresponded to a lower functional diversity of fish communities, confirming our second hypothesis (H2). Exotic species typically had a lower functional diversity than native species, even when accounting for stream order and species richness, which could be at the root of the observed lower functional diversity of highly invaded communities. Our results thus suggest that possible outcomes of severe exotic species invasions could include a reduced functional diversity of freshwater fish communities, and that analyzing data with finer ecological and spatial resolutions could help further investigate these effects.
The different ranking of environmental drivers for exotic and native species might indicate that there are some differences in the relevance of factors shaping their respective functional diversity patterns. For example, our BRT analysis highlighted the role of native marine species entering brackish habitats in contributing to the fish functional diversity in these areas. Our analysis results also suggested that geographical factors like altitude could play a major role to control invasion dynamics, as altitude appeared to be the main factor of divergence for native and exotic functional diversity. In addition to temperature, altitude is linked to a number of other factors, including habitat fragmentation but also habitat quality. Mountain streams typically have lower habitat complexity and resources availability but higher habitat quality than rivers in the lowlands27,28. However, upland streams can be more fragmented due to the presence of dispersal barriers29,30.
The high degree of invasion in the lower stretches of all basins also implies that lowland areas constitute a massive reservoir of exotic species, which could exert a significant propagule pressure towards higher altitudes31,32. This pressure is held in check, at least in some cases, by impassable barriers to migration and by unfavorable ecological conditions in upstream areas33. Habitat fragmentation and other consequences of water abstraction for human use are a worldwide issue often deeply interlinked with exotic species invasions34,35, freshwater communities functional diversity36, and homogenization processes37. Unfortunately, habitat connectivity restoration near heavily invaded areas may not be sufficient to restore native biodiversity, as it could favor exotic invasions in upstream areas. Furthermore, some fish species can colonize upstream areas through human-mediated transport even if migration barriers are not removed38,39. Our analysis did not show sharp regional invasion gradients in several non-fragmented rivers, albeit the spatial resolution of our data should have been sufficient to underline habitat fragmentation at that scale. Overall, exotic species invasion was still halted in upland areas, and we can only speculate that this could be attributed to natural gradients (e.g. temperature and habitat factors).
Perhaps the most immediate outcome of our work is a visualization of the degree of exotic freshwater fish invasion on a broad geographical scale, which highlighted the alarming decline in native species richness and abundance at most sites in the lowlands of our focus area. A high degree of invasion corresponded also to a lower functional diversity of invaded communities. Canals were among the most heavily invaded areas and a good example of this mechanism, with lower overall functional diversity than rivers in the same geographical area. This could be the result of simplified habitats, but also of the native species extinctions in this area, where only few natives remained at the time of sampling9,10,11, but still showing higher functional diversity for the few remaining native species than for the more abundant exotic ones. Lowland rivers with high-stream order class had generally more heterogeneous habitats than canals, offering more spatial and trophic niche options to species and thus presumably allowing for a higher overall functional diversity, yet showed a similar functional diversity to that detected in canals. This pattern could potentially be caused by the presence of exotic species, which have led to local homogenizations of the fish fauna in lowland rivers and canals that are connected to them. Conversely, stream order 2 sites had a very low invasion degree (at the extreme left of Fig. 3), yet the functional diversity of these communities was not much higher than more invaded sites. This perhaps indicates that invasion in these areas is not at the initial stages, and has caused a very moderate decline (if any) of native species, thus resulting in communities where the moderate loss is compensated by a corresponding increase in functional diversity provided by exotic species. Invasion in these areas probably did not progress further because of environmental constraints (e.g. water temperature or habitat) or competition with native species11.
Although there are only few documented cases of extinctions in freshwater fish6,7,8, previous studies in this area have highlighted how freshwater fish invasions could substitute native with exotic species, resulting in a decrease of native taxonomic diversity at the local level, especially when abundances are taken into account9,10. Taxonomic and functional substitutions have been previously advocated as non-detrimental or even beneficial for the environment and biodiversity40,41. However, from a functional perspective, our data show that exotic species generally had a lower functional diversity than native species across similar habitats, which was evident in higher order streams and in canals. As a consequence, substitution of native species with exotic ones could decrease the overall functional diversity of the community.
It has been previously suggested that the availability of habitat and trophic niches in the environment might be the ultimate factor that defines a limit to the possible expression of functional diversity. In this framework, if more species are added to the system functional diversity should not be positively affected, but rather show an asymptotic saturation of the available functional space42. While our results clearly show this asymptotic saturation, they also show that, in invaded communities, the functional diversity of exotic species was generally lower than native ones, perhaps because fewer exotic species with a limited number of traits are introduced. This could be at the root of the decrease observed in overall functional diversity of the highly invaded communities (>50% invasion degree), as the exotic component becomes predominant and drives the overall functional diversity of the community. This also suggests that a relatively smaller number of species with similar traits could be leading the invasion in the lowlands, and that high degrees of invasion might coincide with higher faunal homogenization i.e. Figure 3 of the present work,11,19,43. Our view is also supported by the relatively higher functional diversity of native species in the lowlands, despite their richness not being higher than exotic ones44. A possible explanation is the longer evolution history in the area, which allowed natives to be more functionally diverse and match the ecological niches in their native environment more closely45, albeit no studies have dealt with this aspect in fish communities, so far. Thus, our results clearly indicate that the overall functional diversity of highly invaded communities might be diminished, resulting in an overall heavier loss of diversity and adding a new dimension to the effects of exotic invasions, which should be of great interest to conservation science and management but needs to be further investigated. Our results disagree with the outcomes of the study by12, conducted at the macro-regional scale. This difference could stem from our different approach to the same problem: we analyzed a set of sites at a late-invasion stage rather than using a pre/post invasion approach, but also used species abundance measures and a finer spatial scale, which were advocated by12. The divergence in our outcomes indicates that different spatial and biological resolutions, and different study setups, deeply affect the results of the analysis (see e.g.46). It is also possible that the number and quality of functional traits investigated could further modulate the ultimate outcome of the analysis, suggesting that the role of these factors in affecting the results would need to be tested by future studies.
Given our results, we encourage further studies that deal with functional diversity and, more in general, a wider use of functional approaches before a consensus is formed on the ultimate effects of exotic invasions on diversity. Spatially- and temporally-detailed measures of species abundances would be crucial in these studies, as subtler changes in abundances could be detected before regional extinctions occur. Future studies should investigate in further detail whether a true functional loss has occurred, using a wider geographical scale and comparing areas with different invasion degrees (or theoretical reference communities). Furthermore, large-scale patterns of community diversity could be explored using measures of relative functional diversity, to express the amount of functional diversity expressed by an average species unit. We also encourage studies that explore these mechanisms in different taxa and, whenever possible, compare results across taxa. Ultimately, other aspects of functional diversity should also be investigated, as they might provide useful insights on the final outcomes of exotic and native species interactions. These are all necessary elements to fully understand the ecological consequences of functional diversity loss and its significance for ecosystem functioning at large.
Materials and Methods
Study area
Northern Italy is a dwelling area for more than 17 million humans, with consequent agricultural activities and livestock farming. This area has a Mediterranean continental climate, with an annual average precipitation of 1036 mm and a mean temperature of 12 °C. The largest river basin in Italy, the Po River basin (71,000 km2), is included in this area and we focused our investigation on the Po River itself (in all its course), the Oglio River (one of the most important hydrographic left tributaries of the Po River) and the hydrographic right tributaries in the Emilia-Romagna region. We also analyzed two additional groups of rivers outside of this basin: the Brenta River (north-east of the Po River basin) and the Romagna rivers basin (south of the Po River basin). Overall, a total of 335 sampling sites in 105 watercourses were included in this study, covering a wide range of freshwater habitats, different altitudinal zones and environmental conditions (see Supplementary Fig. 1).
Organic material originating from villages, small towns and livestock farms is the main source of river pollution for rivers in the uplands. Conversely, urbanization and intensive agriculture, causing high nutrient loads and consequent eutrophication, are the main factors affecting lowland rivers10. A network of drainage canals was established in the lowlands to support agricultural irrigation around the 19th century, with hydrological management aimed at supplying both irrigation and drainage needs10,47. The study area was in a late invasion stage11 at the time of sampling, since loss of native species and exotic invasion occurred already prior to 199710, before the data analyzed here were collected.
Data collection
Fish data were collected within monitoring programs of the Emilia-Romagna region48, the Padova Province49, the Po River50 and the Oglio River51 over a relatively long-term period (1999–2010). Despite this, community turnover was not a relevant factor in our study, because fish communities are typically stable over such timescales and, despite the ample interval, most of the data were collected over a limited timeframe44,52. Fish sampling was performed between spring and autumn by electrofishing, combined with nets in sites of higher water depth and conductivity (e.g. lower stretches of the rivers), further details on fish sampling procedures are described in8,11.
Fish species were classified according to53, taking into account recent taxonomic determinations and common names as listed in FishBase54. Each species was categorized as native or exotic: a species was considered as native when naturally present in a specific basin and as exotic when human-introduced, irrespective of the time elapsed since the introduction. Hybrid specimens or uncertain species were excluded from this study in order to avoid taxonomic asymmetries.
Abundance of each species was expressed with Moyle classes55 ranging from 1 (lower abundance, 1–2 individuals per site) to 5 (higher abundance, more than 50 individuals per site). Unfortunately, classes of numerical abundance tend to overestimate the ecological significance of small-bodied species and underestimate that of large-bodied ones, but direct measures of body-mass were not taken at the time of sampling. Attempting to overcome this, a weight was assigned to each species based on their average adult size in the area (1 = small body up to ∼150 g; 2 = medium body ∼150–400 g; 3 = large body over ∼400 g, derived from FishBase and unpublished data) and multiplied by Moyle abundance classes, in order to obtain a body-mass-corrected abundance, hereafter defined simply as abundance24.
Water physicochemical sampling was performed with standard methods, in temporal and spatial proximity to the fish sampling, by different Regional Environmental Protection Agencies (ARPAs, in Italian) for the Po, the Brenta and the rivers in the Emilia-Romagna region. The Oglio River Water Authority carried out the water sampling in the Oglio River. The geographical position and the elevation of each site were recorded. Eight physicochemical variables were monitored: water temperature, electrical conductivity, chemical oxygen demand, biological oxygen demand, total suspended solids, total phosphorus, ammonia (NH4+) and nitrate nitrogen (NO3−).
Land cover data were obtained from the CORINE database (2012, https://www.eea.europa.eu/data-and-maps/data/copernicus-land-monitoring-service-corine). In the lowlands, where estimation of watershed areas is more difficult due to low slopes and human-regulated flow, the land cover of the whole river basin or of the administrative province was used. CORINE land cover classes were merged in five categories based on the main land use in order to better describe the study area: urban use, agricultural use, forest, other natural area, freshwater and brackish water. Land cover was expressed as the share of each of these categories in the watershed of each site.
Fish functional traits
In order to investigate the functional composition of fish communities, five different ecological functions were examined: feeding, reproduction, migration, tolerance and habitat use. Within these ecological functions, all fish species were classified in guilds, each representing an ecofunctional trait (Supplementary Table 1, see also Noble, et al.56). Ecological functions, guilds and classification for most species in this study were taken from24, where all available information was used to identify appropriate guilds for each species. The same methodology was applied to classify euryhaline species that were not included in previous work (see Supplementary Table 2). A total of 59 fish species were sampled in the study area; of these 37 were native and 22 were exotic species, relatively to the national territory.
Influence of environmental variables on functional diversity
Fish functional diversity was investigated through the functional dispersion metric FDis,57. FDis was calculated using one matrix with species abundance and another matrix of functional traits of fish communities through the dbFD function58. FDis calculates the spread of traits in multidimensional space by measuring the distance of each species from a centroid weighted by the species abundance.
As a result, FDis measures the relative diversity of functional traits in a community, communities with larger FDis values are more diverse in species traits combinations, whereas low FDis values refer to communities with more traits in common. The advantages in choosing the FDis metric include that it should theoretically not be overly affected by species richness, can be computed from any distance or dissimilarity measures and take into account any number and type of traits (also qualitative traits, as in this study), is not strongly influenced by outliers and accounts for species abundances57. FDis was calculated for both the whole fish community and for exotic and native species separately, using species abundance as weights to define the relative representation of each trait in the community.
The relative influence of geographical and land use features, as well as water physicochemical variables, on functional diversity was investigated with a machine learning method: the boosted regression trees (BRT)59. BRT analysis is an efficient method to describe any non-linear relationships between variables (e.g. thresholds) and incorporate interactions between variables. Compared to traditional regression methods, BRT analysis combines a large number of simple tree models using the boosting technique to improve the predictive performance. BRT analysis was applied to native and exotic species separately, to investigate differences between the factors affecting the two categories. The relative influence (positive or negative) of each variable was determined by the prevalent direction of its effects on the functional dispersion of species, as no clear unimodal trends were found in our data.
Guild and species abundance matrices for each sampling site were Hellinger transformed60,61 to standardize variations among both species and community size, respectively. Environmental variables expressed as percentages were arcsine transformed, while others were log-transformed. All statistical analyses were performed in R software version 3.4.362. FDis was calculated through the homonymous R package58. BRT was performed with the ‘gbm’ R package63, using standard values (Gaussian distribution, bag fraction of 0.75 and shrinkage of 0.001).
Impacts of exotic invasions on functional diversity
The degree of exotic invasion in each site was estimated through the abundance of exotic versus native species in each community (% of exotic abundance). This measure was then spatially analyzed for all watercourses within each basin in the study area, through linear kriging64.
We also evaluated the spatial overlap and correlation between invasion degree and functional diversity (as expressed through the FDis metric calculated using ecofunctional traits,65). First, we used linear kriging to represent functional diversity, based on the whole community, for invaded sites (i.e. where both native and exotic species were present, 198 out of 335 sites, excluding an invaded site at an altitude >400 m a.s.l.). We selected these sites because we could clearly evaluate invasion effects only where both exotic and native species are present, but also, more importantly, because it allowed us to focus our analysis on a geographically uniform area (the lowlands, below 400 m a.s.l.) and thus avoid potentially-confounding geographical factors (i.e. altitude, as species-poor communities at high altitudes typically have low functional diversity, even if not invaded). This area also had a relatively homogeneous fish community at baseline conditions, so that we could consider the functional diversity of least invaded communities as a close proxy to the baseline conditions of reference communities. Following this, we used the same sites to also explore the variations in overall functional diversity along the invasion degree gradient (between 2.3% and 96.5%).
To investigate the impact of exotic invasions on functional diversity we also examined the variations in native and exotic species richness for each stream order, alongside the respective functional diversity of each site, trying to detect patterns indicative of invasion effects. Stream order for each site was calculated using a Digital Elevation Model (DEM) (http://www.sinanet.isprambiente.it/it/sia-ispra/download-mais/dem20/view). Flow direction and accumulation, as well as the watershed of each sampling site, were calculated based on the DEM layer. For the entire river network generated by flow accumulation, stream order was derived with the Strahler method66. This procedure was reliable for upland streams, but it was less so in the lowland, possibly due to the fact that flow direction and magnitude in the lowlands are not always natural because of human intervention. The stream order was thus manually checked and revised when necessary in lowland rivers and streams. Rivers were grouped into four classes based on stream order: class 1 (Strahler stream order 1 and 2), class 2 (stream order 3 and 4), class 3 (stream order 5 and 6) and class 4 (stream order > 6). As drainage and irrigation canals could not be assigned into any natural class, a separate class called “Canals” was created. Canals are man-made environments, usually characterized by low habitat heterogeneity and controlled hydrology, located in the lowlands southwest of the Po River, near its delta.
We further examined the relationship between native and exotic richness and their relative functional diversity. We used non-linear regressions, evaluated through a combination of Akaike’s Information Criterion (AIC) weights67 and R2 values, to estimate best-fitting curves for functional diversity distributions versus species richness and invasion degree. Non-linear best-fitting regressions were estimated using the Curve Expert Professional 2.6 software68 and the results of the five best regression models for each relationship were reported in Supplementary Table 3. The spatial analyses were performed with ArcGIS software69, using its Kriging Tool and Hydrology Spatial Analyst Tool.
Data availability
Underlying data for this paper is made available also for review as Supplementary Material and will be made publicly available through OpenScienceFramework should the manuscript be accepted for publication.
References
Hooper, D. U. et al. Effects of biodiversity on ecosystem functioning: a consensus of current knowledge. Ecological monographs 75, 3–35 (2005).
Wiens, J. J. Climate-related local extinctions are already widespread among plant and animal species. PLoS Biology 14, e2001104 (2016).
Leprieur, F., Beauchard, O., Blanchet, S., Oberdorff, T. & Brosse, S. Fish Invasions in the World’s River Systems: When Natural Processes Are Blurred by Human Activities. PLOS Biology 6, e28, https://doi.org/10.1371/journal.pbio.0060028 (2008).
Sala, O. E. et al. Global biodiversity scenarios for the year 2100. science 287, 1770–1774 (2000).
Sax, D. F. & Gaines, S. D. Species diversity: from global decreases to local increases. Trends in Ecology & Evolution 18, 561–566 (2003).
Townsend, C. R. Invasion biology and ecological impacts of brown trout Salmo trutta in New Zealand. Biological Conservation 78, 13–22 (1996).
Kitchell, J. F., Schindler, D. E., Ogutu-Ohwayo, R. & Reinthal, P. N. The Nile perch in Lake Victoria: interactions between predation and fisheries. Ecological Applications 7, 653–664 (1997).
Lanzoni, M., Milardi, M., Aschonitis, V., Fano, E. A. & Castaldelli, G. A regional fish inventory of inland waters in Northern Italy reveals the presence of fully exotic fish communities. The European Zoological Journal 85, 1–7, https://doi.org/10.1080/24750263.2017.1415384 (2018).
Milardi, M., Lanzoni, M., Gavioli, A., Fano, E. A. & Castaldelli, G. Long-term fish monitoring underlines a rising tide of temperature tolerant, rheophilic, benthivore and generalist exotics, irrespective of hydrological conditions. Journal of Limnology. https://doi.org/10.4081/jlimnol.2018.1745 (2018).
Castaldelli, G. et al. Introduction of exotic fish species and decline of native species in the lower Po basin, north-eastern Italy. Aquatic Conservation: Marine and Freshwater Ecosystems 23, 405–417, https://doi.org/10.1002/aqc.2345 (2013).
Milardi, M. et al. Run to the hills: exotic fish invasions and water quality degradation drive native fish to higher altitudes. Science of The Total Environment 624, 1325–1335, https://doi.org/10.1016/j.scitotenv.2017.12.237 (2018).
Toussaint, A. et al. Non-native species led to marked shifts in functional diversity of the world freshwater fish faunas. Ecology letters (2018).
Cadotte, M. W., Carscadden, K. & Mirotchnick, N. Beyond species: functional diversity and the maintenance of ecological processes and services. Journal of applied ecology 48, 1079–1087 (2011).
Poff, N. L. & Ward, J. Physical habitat template of lotic systems: recovery in the context of historical pattern of spatiotemporal heterogeneity. Environmental management 14, 629 (1990).
Comte, L., Cucherousset, J., Boulêtreau, S. & Olden, J. D. Resource partitioning and functional diversity of worldwide freshwater fish communities. Ecosphere 7, e01356 (2016).
Ruesink, J. L. Global analysis of factors affecting the outcome of freshwater fish introductions. Conservation Biology 19, 1883–1893 (2005).
Olden, J. D., Poff, N. L., Douglas, M. R., Douglas, M. E. & Fausch, K. D. Ecological and evolutionary consequences of biotic homogenization. Trends in ecology & evolution 19, 18–24 (2004).
McKinney, M. L. & Lockwood, J. L. Biotic homogenization: a few winners replacing many losers in the next mass extinction. Trends in ecology & evolution 14, 450–453 (1999).
Leprieur, F., Beauchard, O., Hugueny, B., Grenouillet, G. & Brosse, S. Null model of biotic homogenization: a test with the European freshwater fish fauna. Divers. Distrib. 14, 291–300 (2008).
Van Kleunen, M., Weber, E. & Fischer, M. A meta-analysis of trait differences between invasive and non-invasive plant species. Ecology letters 13, 235–245 (2010).
Cuttelod, A., García, N., Malak, D. A., Temple, H. J. & Katariya, V. The Mediterranean: a biodiversity hotspot under threat. Wildlife in a Changing World–an analysis of the 2008 IUCN Red List of Threatened Species 89 (2009).
Médail, F. & Quézel, P. Biodiversity hotspots in the Mediterranean Basin: setting global conservation priorities. Conservation biology 13, 1510–1513 (1999).
Crivelli, A. Are fish introductions a threat to endemic freshwater fishes in the northern Mediterranean region? Biological Conservation 72, 311–319 (1995).
Milardi, M. & Castaldelli, G. A novel approach to an ecofunctional fish index for Mediterranean countries. Ecological Indicators 89, 376–385, https://doi.org/10.1016/j.ecolind.2018.02.022 (2018).
Colin, N., Villéger, S., Wilkes, M., de Sostoa, A. & Maceda-Veiga, A. Functional diversity measures revealed impacts of non-native species and habitat degradation on species-poor freshwater fish assemblages. Science of the Total Environment 625, 861–871 (2018).
Villéger, S., Grenouillet, G. & Brosse, S. Decomposing functional β‐diversity reveals that low functional β-diversity is driven by low functional turnover in E uropean fish assemblages. Global Ecology and Biogeography 22, 671–681 (2013).
Bouska, K. L. Discontinuities and functional resilience of large river fish assemblages. Ecosphere 9, e02351 (2018).
Aschonitis, V. et al. Proposing priorities of intervention for the recovery of native fish populations using hierarchical ranking of environmental and exotic species impact. J. Environ. Manage. 210, 36–50 (2018).
Dynesius, M. & Nilsson, C. Fragmentation and flow regulation of river systems in the northern third of the world. Science 266, 753–762 (1994).
Tockner, K. & Stanford, J. A. Riverine flood plains: present state and future trends. Environ. Conserv. 29, 308–330 (2002).
Lockwood, J. L., Cassey, P. & Blackburn, T. The role of propagule pressure in explaining species invasions. Trends in Ecology & Evolution 20, 223–228 (2005).
Simberloff, D. The role of propagule pressure in biological invasions. Annual Review of Ecology, Evolution, and Systematics 40, 81–102 (2009).
Gavioli, A. et al. Exotic species, rather than low flow, negatively affect native fish in the Oglio River, Northern Italy. River Research and Applications, 1–11, https://doi.org/10.1002/rra.3324 (2018).
Bunn, S. E. & Arthington, A. H. Basic principles and ecological consequences of altered flow regimes for aquatic biodiversity. Environmental management 30, 492–507 (2002).
Johnson, P. T., Olden, J. D. & Vander Zanden, M. J. Dam invaders: impoundments facilitate biological invasions into freshwaters. Frontiers in Ecology and the Environment 6, 357–363 (2008).
Oliveira, A. G., Baumgartner, M. T., Gomes, L. C., Dias, R. M. & Agostinho, A. A. Long‐term effects of flow regulation by dams simplify fish functional diversity. Freshwater Biology 63, 293–305 (2018).
Poff, N. L., Olden, J. D., Merritt, D. M. & Pepin, D. M. Homogenization of regional river dynamics by dams and global biodiversity implications. Proceedings of the National Academy of Sciences 104, 5732–5737 (2007).
Hulme, P. E. Trade, transport and trouble: managing invasive species pathways in an era of globalization. Journal of applied ecology 46, 10–18 (2009).
Olden, J. D., Kennard, M. J. & Pusey, B. J. Species invasions and the changing biogeography of Australian freshwater fishes. Global Ecology and Biogeography 17, 25–37 (2008).
Schlaepfer, M. A., Sax, D. F. & Olden, J. D. The potential conservation value of non-native species. Conservation Biology 25, 428–437 (2011).
Griffiths, C. J. & Harris, S. Prevention of secondary extinctions through taxon substitution. Conservation Biology 24, 645–646 (2010).
Villéger, S., Miranda, J. R., Hernández, D. F. & Mouillot, D. Contrasting changes in taxonomic vs. functional diversity of tropical fish communities after habitat degradation. Ecological Applications 20, 1512–1522 (2010).
Clavel, J., Julliard, R. & Devictor, V. Worldwide decline of specialist species: toward a global functional homogenization? Frontiers in Ecology and the Environment 9, 222–228 (2011).
Gavioli, A., Milardi, M., Castaldelli, G., Fano, E. A. & Soininen, J. Diversity patterns of native and exotic fish species suggest homogenization processes, but partly fail to highlight extinction threats. Divers. Distrib. 25, 983–994 (2019).
Marks, C. O. & Lechowicz, M. J. Alternative designs and the evolution of functional diversity. The American Naturalist 167, 55–66 (2005).
Shuai, F., Lek, S., Li, X. & Zhao, T. Biological invasions undermine the functional diversity of fish community in a large subtropical river. Biol. Invasions 20, 2981–2996 (2018).
Milardi, M., Chapman, D., Lanzoni, M., Long, J. M. & Castaldelli, G. First evidence of bighead carp wild recruitment in Western Europe, and its relation to hydrology and temperature. PloS one 12, e0189517, https://doi.org/10.1371/journal.pone.0189517 (2017).
Regione Emilia Romagna. (eds Giuseppe Castaldelli & Remigio Rossi) (Assessorato Attività Produttive, Sviluppo Economico e Piano Telematico., Bologna, 2008).
Provincia di Padova. Carta ittica: aggiornamento dello stato delle conoscenze sui popolamenti ittici della Provincia di Padova (2010).
Autorita’ di bacino del Po. Monitoraggio dell’ittiofauna e redazione della Carta ittica del Fiume Po. Qualità dell’ittiofauna e del macrobenthos del fiume Po. (2008).
Consorzio dell’Oglio. Progetto di Sperimentazione dei rilasci per il deflusso minimo vitale nel fiume Oglio sub lacuale. Relazione conclusiva 2009–2015. (2016).
Korhonen, J. J., Soininen, J. & Hillebrand, H. A quantitative analysis of temporal turnover in aquatic species assemblages across ecosystems. Ecology 91, 508–517 (2010).
Kottelat, M. & Freyhof, J. Handbook of European freshwater fishes. (Publications Kottelat, 2007).
Froese, R. & Pauly, D. FishBase. Version 03/2017. World Wide Web electronic publication, Froese, R. & Pauly, D. (eds.). Available at, www. fishbase. org. Accessed April 2017 (2017).
Moyle, P. B. & Nichols, R. D. Ecology of some native and introduced fishes of the Sierra Nevada foothills in central California. Copeia, 478–490 (1973).
Noble, R., Cowx, I., Goffaux, D. & Kestemont, P. Assessing the health of European rivers using functional ecological guilds of fish communities: standardising species classification and approaches to metric selection. Fisheries Management and Ecology 14, 381–392 (2007).
Laliberté, E. & Legendre, P. A distance-based framework for measuring functional diversity from multiple traits. Ecology 91, 299–305 (2010).
Laliberté, E., Legendre, P., Shipley, B. & Laliberté, M. E. Package ‘FD’. Version 1, 12 (2014).
Elith, J., Leathwick, J. R. & Hastie, T. A working guide to boosted regression trees. Journal of Animal Ecology 77, 802–813 (2008).
Legendre, P. & De Cáceres, M. Beta diversity as the variance of community data: dissimilarity coefficients and partitioning. Ecology letters 16, 951–963 (2013).
Legendre, P. & Gallagher, E. D. Ecologically meaningful transformations for ordination of species data. Oecologia 129, 271–280 (2001).
R Core Team. R: A language and environment for statistical computing. Vienna, Austria (2017).
Hickey, J. et al. gbm: Generalized boosted regression models. R Package Version 2 (2016).
Papritz, A. & Stein, A. Spatial prediction by linear kriging, in Spatial statistics for remote sensing 83–113 (Springer, 1999).
Milardi, M. et al. The role of species introduction in modifying the functional diversity of native communities. Science of the Total Environment V., 699, p. 134364 (2020).
Strahler, A. N. Quantitative analysis of watershed geomorphology. Eos, Transactions American Geophysical Union 38, 913–920 (1957).
Wagenmakers, E.-J. & Farrell, S. AIC model selection using Akaike weights. Psychonomic Bulletin & Review 11, 192–196 (2004).
Curve Expert Professional 2.6, a comprehensive data analysis software for Windows (Hixson, Tennessee, 2017).
ArcGIS Desktop (Environmental Systems Research Institute, Redlands, CA, 2011).
Acknowledgements
We thank LL.D. V.E. Manduca and Dr. M. Rizzoli of the Fisheries Bureau of the Emilia-Romagna Region for providing the Fish Inventories data in the context of a long-term research collaboration. The Oglio River Water Authority (Consorzio dell’Oglio, in Italian) is also acknowledged for providing fish and water quality data for the Oglio River. We also thank Dr R. Spaggiari and Dr S. Franceschini of the Emilia-Romagna Region Environmental Protection Agency (ARPA-EMR), the Piemonte Region Environmental Protection Agency (ARPA-Piemonte) and the Veneto Region Environmental Protection Agency (ARPAV) for providing the water quality database. We also thank Jyrki Lappalainen, for his kind help with statistical analyses.
Author information
Authors and Affiliations
Contributions
M.M. and A.G. conceived the idea and designed the methodology. A.G. collected and analyzed the data. M.M. led the writing of the manuscript. M.M., A.G., J.S. and G.C. contributed critically to the drafts and gave final approval for publication.
Corresponding author
Ethics declarations
Competing interests
No specific funding was received for this work and the authors declare no competing interests (financial or non-financial).
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
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 license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license 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 license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Milardi, M., Gavioli, A., Soininen, J. et al. Exotic species invasions undermine regional functional diversity of freshwater fish. Sci Rep 9, 17921 (2019). https://doi.org/10.1038/s41598-019-54210-1
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-019-54210-1
This article is cited by
-
Response of functional diversity of fish communities to habitat alterations in small lowland rivers
Biodiversity and Conservation (2024)
-
Trophic niche partitioning among non-native fish species coexisting in a Colombian high Andean reservoir
Biological Invasions (2024)
-
Long-term impact of an invasive predator on the diversity of fish assemblages in a neotropical reservoir
Biological Invasions (2024)
-
Natural and anthropogenic factors drive large-scale freshwater fish invasions
Scientific Reports (2022)
-
Spatial pattern of freshwater habitats and their prioritization using additive partitions of beta diversity of inhabitant piscine assemblages in the Terai–Dooars ecoregion of Eastern Himalayas
Limnology (2022)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.