Land-use changes affect the functional structure of stream fish assemblages in the Brazilian Savanna

We investigated the mechanisms involved in the relationship between landuse changes and aquatic biodiversity, using stream fish assemblages of the Brazilian Savanna (i.e., Cerrado) as a study model. We tested the prediction that landscape degradation would decrease environmental heterogeneity and change predominant physical-habitat types, which in turn would decrease the functional diversity and alter the functional identity of fish assemblages. We sampled fish from 40 streams in the Upper Paraná River basin, and assessed catchment and instream conditions. We then conducted an ecomorphological analysis to functionally characterize all species (36) and quantify different facets of the functional structure of assemblages. We detected multiple pathways of the impacts from landscape changes on the fish assemblages. Catchment degradation reduced the stream-bed complexity and the heterogeneity of canopy shading, decreasing assemblage functional specialization and divergence. Landscape changes also reduced the water volume and the amount of large rocks in streams, resulting in decreased abundances of species with large bodies and with morphological traits that favor swimming in the water column. We conclude that land-use intensification caused significant changes in aquatic biodiversity in the Cerrado, reinforcing the need to pay special attention to this global hotspot.


INTRODUCTION
A primary goal of contemporary ecology is to understand how biotic communities respond to increasing human-induced pressures across ecosystems. Natural landscapes are being strongly modified by land-use changes, a widespread process on Earth, but especially critical for the tropics (Dirzo et al., 2014). Besides holding huge amounts of biodiversity, the tropics are marked by the conflicting relationship between unsustainable economic growth and the maintenance of ecosystem functioning and services (Costanza et al., 1997;Gardner et al., 2013). The Neotropical Savanna (hereafter Cerrado) is a typical example of such relationships. Although it harbors high levels of species richness and endemism, the Cerrado is one of the most threatened vegetation domains in South America; therefore, it is considered a global biodiversity hotspot for conservation (Myers et al., 2000;Mittermeier et al., 2004). Because of favorable climate, relief and soils (Goedert, 1989), the Cerrado has undergone intense agribusiness pressure in recent decades. As a result, more than 80% of the original vegetation has been modified for anthropogenic uses (Silva et al., 2013), but only 7% of the area is legally protected (Soares-Filho et al., 2014). This increasing rate of landscape degradation has led to extensive and intensive impacts on several ecosystems (Santos et al., 2017), including headwater streams and their biota (Carvalho, Tejerina-Garro, 2015;Agra et al., 2021).
Because of the riverscape hierarchical configuration, the replacement of natural vegetation cover by agricultural or urbanized areas in a basin promotes a variety of 3/21 ni.bio.br | scielo.br/ni Ludmilla R.P. Alvarenga, Paulo S. Pompeu, Cecília G. Leal, Robert M. Hughes, Daniela C. Fagundes and Rafael P. Leitão changes in the physical-habitat structure of streams (Wang et al., 2006;Leal et al., 2016;Hughes et al., 2019). For example, catchment deforestation drives soil erosion and sedimentation, and reduces leaf and wood inputs into streams. These processes result in altered substrate types (e.g., from rocky substrate and leaf banks to sand-clay dominated streams), and lower levels of stream-bed complexity and stability, among many other alterations in channel morphology (Leal et al., 2016;Castro et al., 2018;Leitão et al., 2018). Reduced canopy cover means reduced inputs of allochthonous items (e.g., terrestrial insects, fruits and seeds) and increased light penetration, altering energy dynamics for stream food webs (Nakano et al., 1999;Bojsen, Barriga, 2002;Carvalho et al., 2017a). These alterations create two fundamental changes for stream systems: reduced environmental heterogeneity and a shift in predominant habitat types. Consequently, biological assemblages are expected to respond both in terms of diversity and identity/type of species.
The distribution of species over environmental gradients is considered a product of organism-habitat adjustments mediated by their functional traits (e.g., morphology, physiology or behavior). Thus, it is expected that a set of coexisting organisms in a given assemblage reflects, at least in part, the relationship of their traits with the local habitat (Poff, Ward, 1990). Such relationships can be understood both in terms of the variability (i.e., functional diversity) and the predominance of traits (i.e., functional identity) within a biological assemblage (Mouillot et al., 2013). Therefore, in the context of land-use impacts on stream physical habitat, we can expect that reduced environmental heterogeneity would result in lower levels of functional diversity (i.e., Biodiversity-Heterogeneity Hypothesis;MacArthur, MacArthur, 1961). Moreover functional identity may be directly affected by changes in predominant habitat types (i.e., Habitat Templet Hypothesis; Southwood, 1977).
Given the capacity to reveal mechanistic processes driving community assembly (Mouillot et al., 2013), and better performance of functional indexes than pure taxonomic indicators (Teresa, Casatti, 2017), functional approaches have been increasingly applied to understand biodiversity changes over degradation gradients. Particularly for Brazilian headwaters, there has been a remarkable increase in literature focusing on the effects of land use on the functional structure of fish assemblages (Bordignon et al., 2015;Casatti et al., 2015;Roa-Fuentes et al., 2019;Zeni et al., 2020;Teresa et al., 2021). However, some authors have pointed out the context-dependency of this relationship, showing significant differences between study regions resulting from particularities such as geomorphology, landscape degradation level and history (Teresa, Casatti, 2017;Leitão et al., 2018). Additionally, those studies seldom explicitly considered the hierarchical structure of the landscape-habitat-ichthyofauna relationship (but see Dala-Corte et al., 2016 andLeitão et al., 2018). These factors potentially limit our capacity to detect general patterns in the biodiversity response to degradation, and to propose efficient conservation and management strategies for stream fish assemblages. In this context, we conducted a study along a wide gradient of landscape degradation and used path analysis to investigate how land use and cover, mediated by alterations in local physical habitat, affect the functional structure of Cerrado stream fish assemblages. We tested the predictions that landscape degradation would decrease environmental heterogeneity and change predominant physical-habitat types, which in turn would decrease the functional diversity and alter the functional identity of fish assemblages. This index ranges from 0 (absence of the three land-use types) to 400 (catchment fully occupied by urban areas); the higher the CDI value, the greater the level of catchment degradation upstream from each sample site (Ligeiro et al., 2013). The weights established for each type of land use were based on previous studies that considered their severity to the hydromorphological degradation for streams on both temperate (e.g., Böhmer et al., 2004;Feld, 2004;Rawer-Jost et al., 2004) and tropical regions. The latter includes previous assessments conducted for quantifying gradients of catchment anthropogenic disturbance in our study area or in other hydrological units in the Cerrado (e.g., Ligeiro et al., 2013;Carvalho et al., 2017b;Silva et al., 2017;Firmiano et al., 2020;Martins et al., 2020).
Stream physical habitat. Physical habitat measurements were made along the entire 150 m site, which was divided into ten sections by 11 equidistant cross-sectional transects (see Fig. S1; adapted from Peck et al., 2006). At each of the 11 transects, we obtained: the width of the seasonal bed and of the wetted area of the stream, using a metric tape; the canopy coverage in six points, using a convex densiometer; the angles of both margins, using a clinometer; five equidistant measurements of water depth, using a bathymetric ruler, with respective estimates of the embeddedness level, type of substrates (e.g., coarse gravel, fine gravel, sand, leaf bank, roots), and proportion of the transect with fish shelter (e.g., large wood, macrophytes, undercut banks, overhanging vegetation). In secondary transects placed between a pair of contiguous transects, we also measured the wetted width of the stream channel and five equidistant estimates of substrate types. Along the thalweg profile in each of the ten sections, we obtained 10 equidistant longitudinal measurements of water depth, presence of fine sediments (i.e., clay and sand), and type of meso-habitats (e.g., riffles, smooth flow, pool). Along the total area of each section, we also quantified large wood (> 1.5 m long and > 10 cm in diameter) and channel slope (by clinometer).
These measures yielded more than 200 metrics describing several aspects of the stream physical habitat (e.g., channel morphoplogy, substrate, fish shelter; Kaufmann et al., 1999;Hughes, Peck, 2008;Kaufmann, Faustini, 2012). We selected a subset of 13 variables representative of habitat heterogeneity, complexity, stability and type to run our models based on (i) our expert knowledge and (ii) the literature showing which variables can both be affected by land use and affect stream fish assemblage structure (Kaufmann et al., 1999;Junqueira et al., 2016;Castro et al., 2018;Leal et al., 2018;Leitão et al., 2018). Next, we excluded two variables because of their high correlation with others (Pearson correlation > 0.70) (see Fig. S2), reaching a final subset of 11 variables (Tab. 1). The selection of a reduced subset of variables was important to comply with the analytical constraints of the structural equation modeling (see below) in which the number of variables needs to be balanced against the number of samples.  Fish sampling. Fish were collected during daylight hours, using semicircular sieves (80 cm in diameter, 1 mm mesh) and seines (4 m long, 2 m high and 5 mm mesh). The effort was standardized over a period of two hours (12 min per section) of sampling by three people along the 150-m site. In the field, all specimens were anesthetized in a lethal dose of Eugenol's solution, and then fixed in 10% formalin. In the laboratory, fish were preserved in 70% alcohol and identified to species using keys and consultation with specialists. Voucher specimens of all species are deposited in the Fish Collection at the Federal University of Lavras (CIUFLA), Minas Gerais, Brazil (voucher numbers in Tab. S3).
Functional structure of the fish assemblages. To assess the functional structure of the fish assemblages, an ecomorphological analysis was first performed. For this, body mass and morphometric measurements (see Fig. S4) were taken in adult specimens (one to seven per species) and then combined into 12 quantitative ecomorphological traits. Together, these traits represent complementary functional characteristics, such as habitat use, foraging and locomotion (Tab. 2; see Leitão et al., 2018 for further details). The mean values of each trait by species were first scaled (mean = 0 and standard deviation = 1) and used to build a matrix of Euclidean distance between each pair of species. Based on this matrix, we performed a Principal Component Analysis (PCA) to reduce dimensionality and build the functional space with the species pool (i.e., all species sampled in the study). We kept the first three axes of the PCA (64.8% of accumulated variance), which were necessary to accommodate the high quality of the functional space and minimize the loss of information (Villéger et al., 2008;Maire et al., 2015). However, this procedure still led to the exclusion of 10 sites that lacked enough species to calculate the functional indexes in a three-dimensional functional space.
Based on the position of fish species in the three-dimensional functional space, and their relative abundances, five complementary indices were calculated to represent the functional diversity of each assemblage: Functional Richness (FRic), Functional Evenness (FEve), Functional Divergence (FDiv), Functional Originality (FOri), and Functional Specialization (FSpe). FRic is the volume of the functional space filled by the species of an assemblage, representing the combination of the functional traits. FEve measures the regularity of the distance among species and the uniformity of their abundances in the functional space. FDiv indicates how species abundance is distributed in relation to the assemblage centroid in the functional space, being low when more abundant species are closer to the centroid, and high when the most abundant species are at the edges of the functional space (Villéger et al., 2008). FOri quantifies the degree of uniqueness, being expressed by the distance of each species to its closest neighbor in the functional space. FSpe quantifies the degree of specialization of a species, expressed by its distance to the centroid of the functional space, considering all the species pool (Mouillot et al., 2013). FOri and FSpe are initially calculated for each species and then scaled for the assemblage by calculating the mean value of the species in an assemblage. To estimate functional identity, we used the community-weighted mean of a trait (CWM;Lavorel et al., 2008). Because we used an ordination technique (PCA) to construct the functional space (i.e., we did not build the functional space directly with the original ecomorphological traits as the functional dimensions), CWM is expressed as the abundance-weighted average  value for each PC axis (i.e., CWM1, CWM2 and CWM3), and the loading of each trait is considered in the subsequent interpretations regarding the CWM values. All functional indices were calculated using the multidimFD package in R (R Development Core Team, 2019).

Data analysis.
Structural equation modeling (SEM) was performed to investigate the possible ways by which landscape degradation, mediated by changes in stream physical habitat, affect fish assemblage functional structure. This method is capable of simultaneously evaluating several pathways to represent the functioning of the system from relationships between variables at different spatial extents and representing a mechanistic approach to ecological studies (Shipley, 2000). Based on the literature focusing on the relationships between land use, physical habitat and stream biodiversity 9/21 ni.bio.br | scielo.br/ni Ludmilla R.P. Alvarenga, Paulo S. Pompeu, Cecília G. Leal, Robert M. Hughes, Daniela C. Fagundes and Rafael P. Leitão (e.g., De Paula et al., 2011;Casatti et al., 2012;Teresa, Casatti, 2012;Junqueira et al., 2016;Castro et al., 2018;Leitão et al., 2018), we tested two hypothetical models, both with the CDI as a predictor of environmental degradation at the landscape spatial extent (Fig. 3).
The first model presents a structure in which the biological response is the assemblage functional diversity expressed by the five above-mentioned indices (i.e., FRic, FEve, FDiv, FOri and FSpe) and the site predictive variables are metrics that represent stream habitat heterogeneity and environmental stability ( Fig. 3A; see ecological meaning of each physical-habitat metric in Tab. 1). Because of the processes associated with soil erosion, sedimentation, and removal of forest vegetation, we expected that catchment degradation would negatively affect streambed complexity (RP100) and stability (LRBS), meso-habitat diversity (DIV_MH), wood volume (V1W_MSQ), the variation in channel shading by riparian vegetation (VCDENMID), and the amount of natural shelters for fish (XFC_NAT) (Fig. 3A).
The second model presents a structure in which the biological response is the functional identity of the fish assemblages (i.e., CWM1, CWM2 and CWM3), and the predictive variables are the site habitat types ( Fig. 3B; Tab. 1). In this sense, we expected changes in the functional identity in response to landscape degradation mediated by decreased wood volume in the channel (V1W_MSQ) and mean canopy shading (XCDEMID). We also expected that functional identity changes would result from catchment erosion and sedimentation processes that could reduce wetted area (XWXD), the amount of large substrates (PCT_BIGR) and the frequency of fast-flowing meso-habitats (PCT_ FAST), while increasing the amount of fine sediments (PCT_FN) (Fig. 3B). We used Bollen-Stine bootstrap (1000 draws) to evaluate the overall fit of the models, which is an adaptation of the chi-square statistic considered robust to non-normal data distributions (Bollen, Stine, 1992) and measures the correspondence between the model and the observed data structure. Standardized path coefficients that were not statistically significant were retained in the model (i.e., we did not re-specify the structural model a posteriori) but were not included in figures to simplify visualization of the significant paths. Modeling was conducted through use of the lavaan package in R (R Development Core Team, 2019).
Effects of landscape degradation on functional diversity. Catchment degradation was negatively associated with residual pools (RP100), stream-bed stability (LRBS), and canopy cover variability (VCDENMID). Given that RP100 was positively associated with FDiv (+0.40), FEve (+0.32) and FSpe (+0.27), CDI had a negative indirect effect on these three functional diversity metrics through this pathway ( Fig. 4A; see the product of the two significant direct effects linking them, from CDI to RP100 and from RP100 to functional indexes, in Tab. 3). VCDENMID was positively associated with FSpe (+0.38), FRic (+0.36) and FOri (+0.51), meaning a negative indirect effect of CDI on these functional metrics ( Fig. 4A; Tab. 3). On the other hand, VCDENMID was negatively associated with FEve (-0.28) and LRBS was negatively associated with FRic (-0.32) and FOri (-0.36), meaning a positive indirect effect of CDI on functional diversity through these pathways ( Fig. 4A; Tab. 3). Summing up the multiple pathways of land use impacts on functional diversity, there was a stronger total negative effect on FSpe (sum of the products of all paths from CDI to FSpe = -0.23) and FDiv (-0.16), a small negative effect on FEve (-0.03) and FOri (-0.03), and no effect on FRic (Tab. 3). Although FRic was positively influenced by the diversity of meso-habitats, the latter was not associated with the CDI. The number of species strongly and positively affected FRic, but it was not affected by any physical habitat predictor. Large wood volume and natural shelters for fish were positively associated with CDI, but they had no effect on any biodiversity metric (Fig. 4A).
Effect of landscape degradation on functional identity. Landscape degradation at the catchment extent (CDI) also affected aspects of fish assemblage functional identity, mediated by a reduction of the wetted area (XWXD) and the amount of large substrates (PCT_BIGR) in the stream channel (Fig. 4B). Both physical-habitat metrics were positively associated with CWM2, leading to an indirect and total negative effect of CDI on this functional index (sum of the products of all paths from CDI to CWM2 = -0.35; Tab. 3). This means that the greater the CDI score, the less the water volume  The mean canopy shading and the amount of fine sediment were positively associated with, respectively, CWM1 and CWM2, but those physical-habitat metrics were not influenced by catchment degradation (Fig. 4B). CWM3 was not affected by any habitat metric.
12/21 ni.bio.br | scielo.br/ni Landscape degradation and fish functional diversity TABLE 3 | Total and indirect effects of catchment degradation (CDI), mediated by instream physical-habitat metrics, on the functional structure of the stream fish assemblages from the Araguari River basin. FDiv: functional divergence, FOri: functional originality, FSpe: functional specialization, FRic: functional richness, FEve: functional evenness, Nb_sp: number of species, CWM1-3: functional identity. The indirect effect of CDI on each functional index, mediated by instream physical habitat, is calculated by multiplying the standardized path coefficients (taking into account the signal, positive or negative) of the two significant direct effects linking them (arrows in Fig. 4): i) from CDI to a physical-habitat metric; ii) from that physical-habitat metric to the functional index. For instance, the indirect effect of CDI on FDiv via RP100 is -0.16 ((-0.39) * (+0.40)). The total effect of CDI on a given functional index is the sum of all significant indirect effects linking them. For instance, the total effect of CDI on FSpe is -0.23 (via RP100 (-0.11) + via VCDENMID (-0.12)). See graphical representation and the direct-effect paths in Fig. 4.

DISCUSSION
In this study we aimed to understand mechanisms that explain the relationship between land uses (i.e., agriculture, pasture and urbanization) and aquatic biodiversity, using stream fish assemblages of the Cerrado as a model. This relationship was investigated by focusing on how the variability and types of fish functional traits responded to instream habitat changes associated with landscape degradation. We used two independent and complementary frameworks, the first hypothesizing that decreased environmental heterogeneity would decrease functional diversity, and the second hypothesizing that changes in stream habitat types would alter the functional identity of assemblages. We observed multiple pathways and various intensities of impacts of landscape degradation on fish assemblages, in great part corroborating our expectations. Catchment degradation (CDI) reduced the heterogeneity and stability of stream habitats, leading to a decrease in four (FDiv, FSpe, FEve and FOri) out of the five functional diversity indexes tested. Changes in catchment condition also resulted in changes in the habitat types, altering the predominant functional traits (i.e., functional identity) within assemblages. These results were allied to the fact that all the physical habitat metrics that were negatively influenced by catchment degradation (i.e., RP100, LRBS, VCDENMID, WXWD, PCT_BIGR) were those that mediated the paths of land-use effects on the assemblage functional structure. Therefore, we can conclude that land-use intensification resulted in significant changes and losses in aquatic biodiversity in the Cerrado, reinforcing the need to pay special attention to this global hotspot. The negative impacts of landscape degradation on the functional diversity of fish assemblages were mediated by a decrease in residual pool volume. This variable (RP100) is expressed as the relative residual thalweg depth (Kaufmann et al., 1999(Kaufmann et al., , 2009, which contributes to micro-habitat diversity resulting from the variability in depth along the stream (Petty, Grossman, 1996), and is inversely related to silting processes (Kaufmann et al., 1999). Therefore, reduced levels of RP100 are likely explained by soil erosion from watersheds (Allan et al., 1997;Allan, 2004), a typical result of agrobusiness activities in the Cerrado (Casatti et al., 2006;Castro et al., 2018). Given that RP100 is a direct indicator of stream-bed complexity (Kaufmann, Faustini, 2012), its reduction should affect diversity components of the fish assemblages. Indeed, this degradation pathway reduced assemblage functional divergence (FDiv), specialization (FSpe) and evenness (FEve).
Another important pathway of impact resulting from landscape degradation is the reduced variability in channel shading by canopy vegetation (VCDENMID), also an indicator of reduced environmental heterogeneity of these systems (e.g., deforested areas can be completely homogeneous in terms of shading over the channel). The distribution of organisms within streams frequently follows patch dynamics because of the heterogeneous distribution of habitat and food resources (Petty, Grossman, 1996;Leitão et al., 2015). For example, several fishes forage in leaf banks or close to tree roots and branches, which are more frequent where marginal vegetation is denser (Da Silva Gonçalves et al., 2018). On the other hand, periphyton feeders tend to occur in areas with greater light incidence and greater algae production (Bojsen, Barriga, 2002), a condition found even in pristine systems (e.g., wider stream stretch not fully covered by the canopy or running in natural riparian-forest clearings). Therefore, scaling up these specificities to the assemblage level, we can expect a positive relationship between the variation in shading over the stream channel and functional diversity. This expectation was corroborated in this study, given that three different functional diversity metrics (FRic, FSpe and FOri) were positively associated with VCDENMID.
Overall, the most negatively affected metrics of functional diversity were FSpe and FDiv. The first was doubly reduced, by both the decrease in RP100 and VCDENMID in degraded catchments, and the latter was the most negatively affected metric considering a single pathway. In biological and ecological terms, reducing FSpe and FDiv means reducing the level of specialization and the number of individuals with highly specialized ecomorphological traits (i.e., far from the center of gravity of the functional space; Villéger et al., 2008). Assuming the fine-tuned relationship between morphology and function (Winemiller, 1991), the loss of such traits from an assemblage can represent the disruption of important ecological processes within the system (Leitão et al., 2016), such as the use and partitioning of habitat and food, with consequent changes in patterns of species packing (Winemiller, 1991), and several other ecological interactions among coexisting species (Mouillot et al., 2013).
Not only diversity, but fish functional identity was also significantly affected by land use changes. Landscape degradation reduced the wetted area (XWXD) and the percent of large substrates (PCT_BIGR) in the stream channel, which led to a decrease in abundance and occurrence of species with larger bodies and with morphological traits that favor fast swimming and endurance in the water column (e.g., deep caudal fins, pectoral fins positioned lower on the body, and lower values for body transversal surface). Given that XWXD takes into account the channel width and the thalweg depth, it is an indication of the amount of water in the system (Kaufmann et al., 1999;Hughes, Peck, 2008). Large rocks also contribute to the formation of deep mesohabitats, such as pools and backwater, by retaing water and shaping the stream channel (Rincón, 1999). Therefore, reducing XWXD and PCT_BIGR represents the loss of critical habitats for large and mid-water column feeders (e.g., Astyanax and Psalidodon species in the studied systems). On the other hand, shallow streams would favor smaller fishes with high body transversal surface, represented here by the poeciliids, especially the widespread invasive species, Poecilia reticulata Peters, 1859 (see its position in the ecomorphological space, Fig. 5).
These results partially corroborates Carvalho et al. (2017b), who developed fish-based multimetric indices (MMI) to assess the condition of Cerrado streams. Carvalho et al. (2017b) studied our sampled sites in the Araguari River basin and other hydrological units from the Upper Paraná River, and found that higher abundances of characiforms and P. reticulata were indicative of, respectivelly, least-and most-disturbed streams. Besides occupying deeper areas, several stream-dwelling characins are allochthonous feeders, and dependent on least disturbed habitats and adjacent forest (Ferreira, Casatti, 2006). On the other hand, P. reticulata is extremely tolerant to poor water quality (Terra et al., 2013;Carvalho et al., 2019). Therefore, although using completely different aproaches (ecomorphological versus biological/physiological traits), it was possible to identify convergent patters in the relationship between landscape degradation and stream fish assemblages in the Cerrado.
Although the overall results of this study corroborated our initial predictions about the negative relationship between land use changes, stream physical habitat and the 15/21 ni.bio.br | scielo.br/ni Ludmilla R.P. Alvarenga, Paulo S. Pompeu, Cecília G. Leal, Robert M. Hughes, Daniela C. Fagundes and Rafael P. Leitão functional structure of the fish assemblages, we also found unexpected relationships. For instance, although negatively associated with CDI, the relative bed stability of streams (LRBS) reduced FRic and FOri, and VCDENMID reduced FEve. Therefore, summing up all pathways of impacts, the total effect of catchment degradation on each of these three functional diversity metrics was almost null (i.e., negative and positive effects compensated each other by different pathways). Another unexpected result was the positive relationship of CDI with the amount of natural shelter and wood volume within streams. A possible interpretation for this result is the influence of deforestation history in the Cerrado, particularly of the riparian forests. Streams surrounded by more vegetation are likely to have greater wood volume in their channels (Leal et al., 2016). However, in the initial periods following riparian deforestation, high amounts of wood inputs can be expected (De Paula et al., 2011). Also, riparian deforestation leads to increased abundance of aquatic macrophytes, which offer considerable shelter for fishes. Further investigations evaluating landscape degradation history and deforestation timing (e.g., Brejão et al., 2018) may clarify these unexpected relationships.
Tropical landscapes are facing intense alterations, especially resulting from the expansion of agricultural lands over natural ecosystems (Malhi et al., 2014), severely threatening global hotspots of biodiversity such as the Cerrado (Mittermeier et al., 2004). Particularly regarding Brazilian freshwaters, environmental legislation seems insufficient to preserve the ichthyofauna (Dala-Corte et al., 2020;Leal et al., 2020). Additionally, conservation strategies still largely neglect aquatic biota, although it was recently shown that focusing on this component of biodiversity would represent impressive gains with minimal extra financial costs (Leal et al., 2020). Therefore, continued efforts and innovative tools are urgently needed to better understand, control and remediate the human-induced impacts on the functioning of freshwater ecosystems (Leitão et al., 2018). We are in accordance with the growing consensus that one of these tools is to seek for and quantify biological responses taking into account the multiple facets of biodiversity (Mouillot et al., 2013;Teresa, Casatti, 2017). For instance, we did not detect changes in taxonomic diversity (i.e., the number of fish species), despite the fact that different aspects of the functional structure of fish assemblages were clearly responsive to alterations in stream habitat type and heterogeneity. Therefore, we believe that, by demonstrating pathways by which land-use changes affect the functional diversity and identity of stream fish assemblages, this study bring important insights to improve conservation and management strategies to one of the most ecologically and economically important regions of the Neotropics. and Companhia Energética de Minas Gerais (P&D ANEEL/CEMIG-GT487 and P&D ANEEL/CEMIG-GT599 -PROECOS Project). LRPA received CAPES scholarship (88887.199567/2018-00), PSP received CNPq productivity grant (303548/2017-7), CGL received FAPESP funding (2017/25383-0), RMH received a Fulbright Brazil scholarship, and RPL received CNPq funding (436007/2018-5).