Evaluating the inﬂuence of environmental variables on ﬁsh assemblages along Tropical Andes: considerations from ecology to conservation

The lack of knowledge about freshwater biodiversity in remote areas like montane Andean rivers is noticeable and hinder the development of efﬁcient conservation plans for aquatic ecosystems. Habitat requirements of tropical Andean ﬁshes have been poorly analysed, despite the relevance of these ecological features on the ﬁsh conservation strategies. The present study aimed to examine ﬁsh assemblages along the altitudinal gradient in tropical Andean rivers and their relationship to habitat features to detect the most relevant abiotic factors and their implication on conservation. This study was conducted on a regional scale, along six basins located in tropical Andes (500 to 2,692 m a.s.l.). We studied ﬁsh diversity and distribution and their habitats along these elevation gradients. We found 6,320 specimens from 53 genera. Our analysis for the six most common and widely distributed genera ( Astroblepus , Brycon , Bryconamericus , Chaetostoma , Pimelodella and Trichomycterus ) underlines the relevance of hydromorphological variables on ﬁsh community structure in the Andean montane rivers, whereas human alterations seem to be less signiﬁcant. These ﬁndings could be useful for current management and conservation strategies considering present threats over tropical montane rivers.


Introduction
Fish biodiversity is the highest amongst the vertebrates, with about 35,500 species currently recognised, approximately the same species as that of all non-fish vertebrates combined. Amongst them, almost 18,000 fishes are freshwater-living species (Fricke et al., 2020). Despite this group's relevance in ecological and even economic and social terms (Dudgeon et al., 2006;Hermoso et al., 2018), freshwater fishes are more threatened and less known than other vertebrates (Darwall et al., 2011;Miqueleiz et al., 2020). Freshwater ecosystems are some of the most endangered globally, being lost at an alarming proportion, definitely faster than terrestrial biodiversity (Barrett et al., 2018;Reid et al., 2019).
Amongst freshwater fishes, South America is the most diverse area on earth, with around 40% of freshwater fishes occurring in the Neotropical region (Tisseuil et al., 2013;Miqueleiz et al., 2020), where indeed most of the high-priority freshwater ecoregions are located (Hermoso et al., 2017). Neotropical fishes have been studied mainly in the lowlands of major river systems, whereas montane and foothill tropical Andean fishes have been less studied (Anderson & Maldonado-Ocampo, 2011;Tognelli et al., 2019).
In the Neotropics, the Andean region retains the highest diversity of mountain fishes with almost 700 species, with a high percentage of endemic species Tognelli et al., 2019). The uniqueness of fish diversity increases with altitude in the tropical Andes (Jaramillo-Villa et al., 2010;Oberdorff et al., 2019;Jézéquel et al., 2020), making this region a hotspot on fish biodiversity. In this scenario, principal threats over aquatic ecosystems and tropical Andean fishes are deforestation, water withdrawals and pollution (Araújo-Flores et al., 2021), exotic species introduction, hydropower development, land-use and climate change (Anderson & Maldonado-Ocampo, 2011;Lujan et al., 2013;Pelicice et al., 2017).
Despite the relevance of tropical Andean fishes, limited information exists on the distribution, conservation status, population dynamics, and habitat preferences for most of these species. Ecological and natural history knowledge of many species is scarce; several species are known only from the type locality and original description (Reis et al., 2016;Tognelli et al., 2019). Regarding their conservation status, a group of experts lead by IUCN Survival Commission evaluated 666 Andean tropical fishes (Tognelli et al., 2019). A third of assessed species were categorised as data deficient, and population trends are unknown for 77% of assessed species. Even accounting for this study, additional assessments are justified given patterns of both biodiversity and anthropogenic stressors (Hermoso et al., 2018;Ríos-Touma & Ramírez, 2019). Although general habitat priorities are recognised-i.e. aquatic macrohabitat characteristics-, few or nothing is known about specific habitat requirements. A better understanding of fish habitat relationships is essential for designing conservation strategies and efficiently managing tropical Andean freshwater ecosystems.
Mountain fishes have been traditionally considered particularly sensitive to human impacts in water quality and quantity (Liang et al., 2013;Encalada et al., 2019;Richardson, 2019). With the subsequent exploitation of resources, human populations' growth has led to widespread decreases in the native biodiversity of mountain freshwater ecosystems (Anderson et al., 2011;Richardson, 2019). Regarding to tropical Andes, headwater streams are amongst the most likely to suffer high levels of biological diversity loss in the future, given climate change predictions and current impacts (Herrera-R et al., 2020). The high Andean region, specially the northern inter-Andean valleys, have a long history of high human pressure. They concentrate most of the population of the countries along the mountain range notably affecting water systems (Buytaert & De Bièvre, 2012). It is subjected to severe agricultural land use development and deforestation that alter the landscape's capacity to provide ecological services (Leitão et al., 2018;Encalada et al., 2019). Land-use changes across the catchment and riparian zones are essential considerations in constant efforts to protect headwater streams and conserve their biodiversity. These pressures to freshwater ecosystems may also reduce the availability of ecosystem services, such as food security and nutrition, indispensable to human populations (Anderson et al., 2011;Buytaert & De Bièvre, 2012;Lessmann et al., 2019).
Considering this high vulnerability scenario, to implement efficient management and conservation plans for freshwater systems in the region is critical. To that end, a better understanding of freshwater species ecology and diversity distribution is essential. Patterns in fish assemblage and diversification in such a hard environment as headwaters of Andes mountains arise from the interaction between habitat features and specific evolutionary constrains (Grossman et al., 1998;Carvajal-Quintero et al., 2015), providing essential considerations for riverine ecosystem conservation efforts. Nevertheless, this crucial ecological knowledge connecting environmental characteristics and fish diversity is yet scarce and dim (Alofs et al., 2014;Costa et al., 2016). Studies about habitat requirements show that mountain fish assemblages are mainly conditioned by variables related to water quality -as physicochemical parameters-, stream hydromorphology, and local habitat features (Quist et al., 2004;Li et al., 2011;López-Delgado et al., 2020). Habitat preferences of Andean mountain fishes' studies have been elaborated on in specific rivers and places, with sometimes different and controversial results (Chará et al., 2006;Alexiades & Encalada, 2017). Therefore, we are far from understanding the habitat requirement for the majority of Neotropical mountain fish species. Tropical Andes might be subject to exceptionally high extinction risk, as it represents a hotspot with high threat status, small species ranges, and high data deficiency (Alofs et al., 2014;Pelayo-Villamil et al., 2015;Miqueleiz et al., 2020). Knowledge about the biological and ecological requirements of species extant in this region is essential for their conservation. Studies about the Andean region's fish community structure and habitat suitability are scarce, considering specific taxon (Vélez-Espino, 2003, 2006Chará et al., 2006) or particular basins (Jaramillo-Villa et al., 2010;Tobes et al., 2016b;Miranda et al., 2018). To our knowledge, there are no studies on fish assemblage variation about environmental features considering on a regional scale, including trans and cis-Andean basins, considering jointly different basins and rivers present on Andean mountains.
The goals of this study were (a) to examine the distribution of tropical Andean montane fishes along an altitude gradient in different freshwater ecosystems and (b) to learn which abiotic factors might be most important in determining fish assemblage and populations abundance. Our final aim is to discuss the possible implications of our results on fish conservation.
We hypothesised that specific environmental settings and human pressures strongly condition fish assemblages in the Andean region. These abiotic features can be evaluated according to the influence on fish biodiversity. Anthropogenic impacts are strongly conditioning the presence and abundance of freshwater fishes in the Andean foothill, where constraining environmental features have significantly shaped the biological characteristics and distribution of freshwater species (Lujan et al., 2013;Tobes et al., 2016aTobes et al., , 2016b. We anticipated that several of these abiotic factors would substantially modify the structure of river biodiversity and we aspire to find sound patterns explaining this relation along the Andean basins.

Study area
The present study includes 101 sampling points in six basins located in Colombia, Ecuador and Peru into the Andean mountain range ([ 500 m above sea level,  Table 1). From North to South, the first basin is the Suaza River, belonging to the Magdalena Basin, which drains most Colombian Andes to the Caribbean Sea. Sampling surveys include remote-first points in the core area of the Cueva de Los Guácharos National Park . The Hacha River is in the Colombian Andes too, belongs to the Caquetá River Basin, one of the largest tributaries of the Amazon River in Colombia. This river descends by the piedmont to its confluence with the Orteguaza River in the Amazonian plateau (Tobes et al., 2016b).
In Ecuador, three basins were sampled. The upper valley of the Mira River constitutes the border between Colombia and Ecuador. The river continues in Colombia until the mouth into the Pacific Ocean. Mashpi and Chirape River Basins are located in the north-western Andean ridge of Ecuador (Morabowen et al., 2019). Both rivers drain on the Esmeraldas River into the North Andean Pacific Slopes-Rio Atrato freshwater ecoregion (Abell et al., 2008). Finally, Alto Madre de Dios River is the upper stretch of the Madre de Dios River, in the southwestern Peruvian Amazon, tributary of the Madeira Basin, one of the largest tributaries of the Amazon River and the unique tributary on its right riverbank with Andean influence (Benejam et al., 2018). (pst) and cropland (crp) extend (percent cover), soil erosion rate (ero, km ha-1 year-1) and human footprint index (hft, 9 10). These values were obtained from the HydroATLAS database (Linke et al., 2019) Sampling protocol and data acquisition To estimate fish abundance and biomass, specimens were captured by electrofishing surveys (FEG 1500, EFKO-Elektrofischfanggeräte, Leutkirch, Germany and Hans Grassl Model IG200/ 2D, 300-600 V, 0.2-2 A, Schönau am Königssee, Germany), during the dry season between 2004 and 2017, following a single-run depletion methodology and estimation of fish abundance (Meador et al., 2003). Collected fish were anesthetised and, subsequently, counted, their total length measured and then released after the survey. The fish abundances at each sampling point were expressed as catch per unit effort (CPUE). This CPUE was calculated as the number of fishes captured in a sampling site per unit of time (Meador et al., 2003). The wadeable sampling stretches were continuous in all streams, but we ensured that principal habitat types were sampled. Surveys were performed in accordance with the procedures set out in CEN standard EN 14011 and with the corresponding authorisations from the local administrations in the different sampling campaigns. Voucher specimens of each species were taken to the laboratory; and identified to the species level with the nomenclatural assignments according to Eschmeyer's Catalogue of Fishes (Fricke et al., 2020), or to genus level when taxonomic information was insufficient for reliable identification. Carrillo, T. Pamballo (Universidad Tecnológica Indoamérica, Quito, Ecuador), J. Valdiviezo-Rivera (Instituto Nacional de Biodiversidad, Quito, Ecuador) and the authors of this article. Voucher specimens are deposited in Museo de Zoología de la Universidad Tecnológica Indoamérica, Quito, Ecuador; Instituto Nacional de Biodiversidad, Quito, Ecuador; Museo de la Escuela Politécnica Nacional, Quito, Ecuador; Museo de Historia Natural de la Universidad Nacional San Marcos, Lima, Perú; and Museo de Historia Natural de la Universidad de la Amazonía, Florencia, Colombia.
Abiotic data of the river site were measured after fish sampling. These data were categorised according to five categories. Firstly, physicochemical parameters (water temperature (°C), pH, conductivity (lS cm -1 ), and dissolved oxygen (mg l -1 )) were measured by digital multiparameter probes (Hanna Instruments, HI 98,129 Combo Waterpro, oximeter HI 9146,   Along each transect, we measured width (considering only wetted surface), water depth and water velocity, and determinate if the point was shaded or not, and the type of substrate at 1-m intervals, according to fish habitat assessment protocols (Simonson et al., 1994). Shading (%) was visually estimated to nearest 10%, and considering shade when the sun is at its zenith (90°a bove the horizon). Water velocity was measured with a rheometer (manufactured, boat display SpeedCoach OC 2, Nielsen-Kellerman Co., Boothwyn, PA, USA).
Fourth, three habitat quality indices were measured: the qualitative habitat evaluation index QHEI (Gazendam et al., 2011), the Index of Fluvial Habitat IHF and the Riparian Forest Quality Index QBR adapted to Andean mountains (Acosta et al., 2009). The QHEI jointly considers different habitat parameters related to substrate and embeddedness, river hydromorphology (instream structure, water speed and depth, morphological diversity), riparian vegetation cover and channel alteration. The IHF characterises the physical habitat and evaluates the relationship between habitat heterogeneity and the stream channel's physical variables, which are influenced by hydrology and substrate composition. The IHF considers the substrate diversity and embeddedness, the sedimentation in pools, the frequency of riffles and the diversity of flow velocity and depths. The index also evaluates other elements that contribute to increasing physical habitat diversity, as the presence and dominance of different groups of aquatic vegetation, accumulation of leaves and drifted wood (Acosta et al., 2009).
The QBR index considers aspects of the riparian forest, including canopy cover and structure and channel alterations. The QBR is calculated in the field by a field surveyor (an 'observer') who is familiar with the most common tree and shrub species found in the study areas. The index must be calculated in river or stream lengths of 50 m (upstream areas) or 100 m (middle and lower reaches). Both river banks and analysed upstream or downstream from observer should be considered together.
This index is a summation of four parts: (1) the percentage of riparian cover of all kind of plants except annuals; (2) the structure and rate of forest cover, including shrubs and other low-lying vegetation; (3) the quality of this cover, considering native and non-native tree species and the percentage of gallery-type structure; and (4) the naturalness of the river channel, according to morphological changes for human activity, channel reduction, river linearisation  Fig. 3 Redundancy analysis biplot for the sampling sites scores (grey circles) and fish genus scores (red circles). Abiotic variables (green vectors) are plotted as correlations with fitted site scores, and their length is proportional to the explained variability. WT water temperature (°C), OX dissolved oxygen (mgÁl -1 ), CD conductivity (lSÁcm -1 ), pH, WM mean width (m), DM mean depth (cm), VEL water velocity (mÁs -1 ), FIN fines, GRA gravels, PEB pebbles, BOU boulders, SHA shade (%), IHF Index of Fluvial Habitat, QHEI Qualitative Habitat Evaluation Index, crp cropland extend, hft human footprint index and elimination of meanders. These differences are measured in a simple, quantitative way according methodology explained in Munné et al. (2003). The index varies between 0 and 100 points, and it is the sum of these four scores.
Finally, to examine the effect of human activity, four indicators of impacts were considered (Fig. 1): the cropland and pasture extent, the soil erosion and the human footprint. As percent cover in the total watershed upstream of reach pour point, the extension of pastures and croplands is a good indicator of human development and impact. EarthStat provides a global data set of croplands and pastures circa 2000 by combining agricultural inventory data and satellitederived land cover data (Ramankutty et al., 2008). Soil erosion rates expressed as a mass of soil lost per unit area and time (km ha -1 year -1 ) measure land-use changes and human impact. These rates were estimated using a combination of remote sensing, GIS modelling and census data (Borrelli et al., 2017). Finally, the accumulated Human Footprint index in the total watershed upstream of the reach pour point was considered. The Human Footprint represents the relative human influence. This index considers different human pressures: (1) built environments, (2) population density, (3) electric infrastructure, (4) croplands, (5) pasture lands, (6) roads, (7) railways, and (8) navigable waterways. These human impacts were combined in a unique index, the Human Footprint (Venter et al., 2016). These values were obtained for each point from the HydroATLAS database, a global compilation of hydro-environmental river reach characteristics at 15 arc-second resolution (Linke et al., 2019).

Statistical analysis
Analyses were conducted at two levels. At the basic level, we carry out an exploratory analysis based on non-parametric correlations (Spearman's rho) amongst the abiotic variables and fish abundances of more representative species (present in all the basins and collected in more than 20 places). The relationships between fish assemblages and abiotic explanatory variables were explored at sampling points by redundancy analysis. Redundancy analysis is a constrained ordination method, similar to multiple regressions. It estimates the amount of variation in multiple dependant variables based on multiple predictor variables (Legendre & Legendre, 1998).
To carry out the redundancy analysis, environmental variables, except pH, were log-transformed (ln x ? 1) to improve normality and because they were measured using different units. According to the recommendation of this method, logit transformation was used too for percentage data to analyse proportions in ecological research as an alternative to the arcsine square root transformation (Warton & Hui, 2011). Posteriorly, these variables were submitted to a Spearman correlation to identify multicollinearity. In the case of redundant variables (r s [|0.60|, P \ 0.05), those with higher predictive importance for tropical stream fish assemblages, according to the literature, were retained (Prudente et al., 2017). Data on fish abundance were root transformed.
To reduce the influence of rare species on multivariate analyses, poorly represented species (occurring at only ten sites) were removed from this analysis, and taxa were grouped by genus. Because of the low representativeness of many species, data disaggregation to species level causes a severe loss of information. Grouping species into genera allows a more consistent analysis of the obtained data. The significance level was set at P = 0.05. Statistical analyses were conducted with the PAST 4.05 statistical software (Hammer et al., 2001).

Results
A total of 6,320 specimens were considered in this study from 101 surveys carried out in the Tropical Andes ([ 500 m above sea level). From them, seventy-four species were identified, whereas 2,065 specimens (32.7% of the total catch) belonging to 25 genera were not identified until the species level (Annex II displays the number of specimens collected per sampling points). Because of the high proportion of unidentified specimens until species level, posterior analyses have been performed grouping species on genus level. These species are clustered on 53 genera, distributed between 500 and 2,692 m of elevation (Fig. 2).
Only six genera (Astroblepus, Brycon, Bryconamericus, Chaetostoma, Pimelodella and Trichomycterus) were extant in all the basins and more than 20 surveys, being the climbing catfishes (genus Astroblepus) the most frequent species, present on 79 of 101 sampling points and distributed in all the altitudinal ranges. The representative genera of our analyses show a significant correlation between abiotic variables and fish abundances in different ways, depending on genera ( Table 2). The climbing catfishes show the most distinct pattern compared to other genera. Whereas this genus shows preference to cold waters and narrow and shallow stretches, other genera show a positive correlation with the increase of temperature (Bryconamericus, Brycon and Pimelodella) and morphological variables as width and depth (Chaetostoma and Pimelodella). Estimated fish abundances do not show a significant relationship with substrate types, except for genus Brycon, positively correlated to the presence of gravels and pebbles (Table 2; Fig. 3).
Regarding multicollinearity of abiotic variables, significant Spearman correlation values (r s [|0.60|, P \ 0.05) were observed for dissolved oxygen (OX) in relation to the Riparian Forest Quality Index (QBR) (r s = 0.75); Human Footprint (HFT) in relation to pasture extent (PST) and soil erosion (ERO) (r s = 0.70 and r s = 0.66). The variables OX and HFT were retained in redundancy analysis. Finally, we select 17 environmental variables, including four physicochemical parameters (water temperature, pH, conductivity and dissolved oxygen), four hydromorphological variables (mean width, mean depth, water velocity, and percentage of shade), two habitat quality indices (QHEI and IHF), five substrate types, and two human impact indicators (cropland extent and the human footprint).
The total fraction of the variance explained by the redundancy analysis was 39.4%. The first two canonical axes explained 17.5% and 10.8% of the data variance, respectively. The redundancy analysis's first axis was principally influenced by hydromorphological variables (mean width and mean depth of sampling points) and physicochemical parameters (water temperature, conductivity and dissolved oxygen). In contrast, the water temperature was again the primary determinant of assemblage structure along the second axis, and additionally, other variables as the percentage of bedrock substrate, dissolved oxygen or water velocity (Fig. 3).
Climbing catfishes (genus Astroblepus) are separated from the remaining, and they are closely related to high mountain stretches of the rivers, with lower temperatures and higher environmental quality indices. The pencil catfishes (genus Trichomycterus) appear in a second position regarding altitude (Fig. 2). They are shown with climbing catfishes in the upper river stretches of high mountains. Three genera of characoids (Bryconamericus, Brycon and Pseudochalceus) are the next genus in the altitudinal gradient, related to the presence of gravels and rocky substrates, high pH values, and shadow places. In places with lower elevations, armoured catfishes of the Chaetostoma genus and other genera are related to more profound and wider river stretches (Table 2; Fig. 3).
Regarding the human impact indicators, cropland and pasture extent are negatively correlated with some genera (Chaetostoma and Pimelodella), showing their preference for low-impact places (Table 2). However, climbing and pencil catfishes (Astroblepus and Trichomycterus genus, respectively), the inhabitants of the highest areas (Fig. 2), exhibit a positive correlation with pastures, as well as a positive association with the soil erosion in the case of Astroblepus, and with the human footprint in the case of Trichomycterus. These two genera are related to cropland extent in the redundancy analysis (Fig. 3).

Discussion
Fish beta-diversity in the Andes is higher than for other montane places in the world because of the particular evolutionary process in the Neotropical montane fishes (Schaefer & Arroyave, 2010;Oberdorff et al., 2019;López-Delgado et al., 2020). Ninetynine taxa (seventy-four species, twenty-five undetermined) were registered in this study, in addition to a high percentage of endemic taxa (Annex II). Whereas at headwaters only few but endemic species are present, the number of species grows downstream, according to the ordinary succession of fish species in the Andean elevation gradient (Lujan et al., 2013) (Fig. 2). The increase in habitat diversity, environmental stability and trophic resources from upstream to downstream cause the rise in species, as it has been widely reflected in the literature (Carvajal-Quintero et al., 2015;Benejam et al., 2018).
Ecological attributes determine fish assemblages in freshwater ecosystems and show diverse relationships amongst them. According to previous studies, analysed fish populations are distributed through altitudinal gradient depending on habitat characteristics (Fig. 2), from torrential headwater to moderate current in the intermediate places on the piedmont, between mountain upstream, and lowland reaches (Jaramillo-Villa et al., 2010;Encalada et al., 2019). Regarding substrate preferences, tropical Andean rivers present in general rocky substrates, increasing the diversity of sizes downstream by the erosion effect, offering shelter to mountain fish species and their principal food resources (Anderson & Maldonado-Ocampo, 2011;Lujan et al., 2013;Miranda et al., 2018).
Redundancy analysis presents similar results regarding fish distribution in the altitudinal gradient and their relationships with abiotic variables (Fig. 3). Vector direction and length of environmental variables in this analysis indicate the relative magnitude and influence of a particular variable on the fish assemblages (Borcard et al., 2011). Water temperature is considered one of the most relevant factor conditioning river fish assemblages, widely discussed in the literature (Kriaučiūnien_ e et al., 2019;Moody et al., 2019). In mountain ecosystems, the variation of temperature across the elevation gradient determines more remarkably the local patterns of fish populations (Carvajal-Quintero et al., 2015;Reis et al., 2016). Water temperature directly affects the metabolism, growth and fecundity of fishes, and it is considered one of the most important factors limiting fish assemblages in high elevated places (Tejerina-Garro et al., 2005;Jaramillo-Villa et al., 2010). Mountain ecosystems show low metabolism rates due to cooler habitats, and only fishes adapted to low growth and fecundity rates, principally invertivore species, can survive in the highest elevations (Román-Valencia, 2001;Benejam et al., 2018), whereas fish diversity, richness and abundance increase in lower altitudes because of habitat diversification and the rise in metabolism (Jaramillo-Villa et al. 2010;Carvajal-Quintero et al., 2015).
However, whereas dissolved oxygen concentration should be negatively correlated with water temperature, the highest elevation places ([ 1,000 m a.s.l.) show low levels of dissolved oxygen (Fig. 3) due to low partial pressures of atmospheric oxygen at high mountain altitudes (Lujan et al., 2013). This lack of oxygen and low temperatures in the highest places could determinate the presence of only adapted species because of low metabolism (Moody et al., 2019;Cornejo et al., 2020).
Hydromorphological variables are considered relevant on fish habitat relationships too: steep stretches with the abundance of rocky substrates and speedywaters impede the access of predators and promote the presence of adapted running-water species. Besides, shadow-water shores are used as refuges for juveniles of some species. On the other hand, the existence of pools permits the occurrence of other bigger species (Jaramillo-Villa et al., 2010;Lujan et al., 2013).
Our results highlight the relevance of these hydromorphological variables on fish community structure in the Andean montane rivers; scores of redundancy analysis indicate which environmental variables are significant in determining the fish community composition (Braak & Verdonschot, 1995). Unexpectedly, the direct human impact seems to be less critical (Table 2; Fig. 3). Furthermore, climbing (Astroblepus) and pencil (Trichomycterus) catfishes present a positive correlation with some human impact indicators (Table 2; Fig. 3). These correlations highlight the capacity of these genera to subsist in high impacted places; their abundances are more conditioned by other variables than human activities, as will be discussed later.
Although these high correlations with some types of human alterations on the basin territory of these genera, our studied basins do not have substantial modifications compared to other basins in the Andean region. Obviously, artificial changes as cropland and pasture development or urban expansion and pollution have definitely adverse effects on fish conservation (Fu et al., 2003;Leitão et al., 2018). However, some representative montane species have been localised in highly impacted places in the present study and others. Fish biodiversity seems to be more conditioned by hydromorphological characteristics than anthropic alterations. In a recent study, positive relationships amongst highly modified agricultural streams and abundances of the climbing catfish Astroblepus vaillanti (Regan, 1904) have been found (Alexiades & Encalada, 2017). Authors conclude that A. vaillanti distribution might be driven by an invasive predator displacement, looking for refuges in the altered places. Other researchers found that variables associated mainly with river size and substratum structure were significant determinants of fish assemblage (Chará et al., 2006;Rodrigues-Filho et al., 2017). In addition, the climbing catfish Astroblepus mariae (Fowler, 1919) was detected in high-altered places, where morphological features of the creeks determine fish distribution (Buitrago-Suarez et al., 2015). In other regions, fish community structure is related to water quality and habitat structure variables, but mainly to hydromorphological features (Li et al., 2011;Tobes et al., 2016a).
Climbing catfishes are the most common genus in this study, showing the highest abundance and occurrence values, and reaching the highest elevated places (along with pencil catfishes only for the Alto Madre de Dios Basin, Fig. 2). Astroblepus is the principal representative fish genus of the Andes' high mountains, widely distributed from Panama to Bolivia, being the unique species present at sites close to the headwaters (Schaefer & Arroyave, 2010;Benejam et al., 2018). This genus is particularly adapted to steep-torrential streams with a rapid current, using these stretches as refuges from predators, and using local food resources (Pouilly et al., 2006;Lujan et al., 2013). The knowledge of this genus is scarce in general, especially concerning biogeography and systematics (Schaefer & Arroyave, 2010;DoNascimiento et al., 2017). Regarding conservation status on the Astroblepus species, most of the information available is for only one species, A. ubidiai Pellegrin, 1931, in the Mira River Basin (Vélez-Espino, 2003, 2005, 2006. This species is catalogued as critically endangered by IUCN Red List, considering the total area of occupancy, the subpopulation's fragmentation, and the environmental deterioration (water extraction, pollution, habitat loss, fishing and invasive specie's impact). Although the distribution and population's status of most Astroblepus species is unknown, general threats over A. ubidiai are generalised on the Andes (Arratia, 1983;Román-Valencia, 2001). It is expected similar conservation status than A. ubidiai for other species of this genus (Román-Valencia, 2001;Buitrago-Suarez et al., 2015;Alexiades & Encalada, 2017). Their restricted distributions, limited dispersal capacity, high adaptation to extreme environmental conditions, and the preference of invasive trout for those habitats (Vimos et al., 2015) increase the probability to categorise these fishes as threatened (Anderson & Maldonado-Ocampo, 2011).
As discussed above, the hydromorphological features of Andean mountain streams are especially relevant to the habitat requirements of these fish species (Table 2; Fig. 3). However, habitat alteration has been poorly studied in the Andean montane streams, despite the broad knowledge of these impacts Richardson, 2019). The adverse effects of these alterations on mountain headwaters, as the development of small obstacles, can be similar to those reported for larger dams regarding fish diversity conservation (Alexandre & Almeida, 2010;Tobes et al., 2016a). One of the main consequences of constructing these small obstacles is the hydromorphological modification of stretches, affecting fish communities (Tobes et al., 2016a;Latli et al., 2019). Similarly, other habitat alterations as channelisation or dredging produce significant modifications on fish assemblages (Fischer et al., 2012;Reis, 2013;Chiu & Suen, 2016).
To our knowledge, before this study, there was little published information on the Andean montane fishes and habitat relationships, and most of the research was conducted mainly considering a particular species (Vélez-Espino, 2003, 2006Chará et al., 2006) or certain places and basins (Jaramillo-Villa et al., 2010;Tobes et al., 2016b). With our analysis, regional ecological patterns explaining habitat preference arise, and we through some light on the main drivers for fish distribution and assemblage composition. Whereas taxonomic determination is the first step on species conservation, knowledge of population trends, ecological features, and principal threats are necessary to determine the conservation status of those species (Araújo-Flores et al., 2021). Fish habitat is a key component to delimit these statuses and develop successful conservation plans (Costa & Barletta, 2016;Tognelli et al., 2019).
Finally, some considerations should be discussed regarding the lack of knowledge and the known taxonomic impediment (Giangrande 2003;Donaldson et al., 2016). In this study, taxonomic experts were not able to identify until a species level a third of the studied specimens. Whereas taxonomists' expertise is beyond any doubt, the lack of knowledge explains the incapacity to identify these species (Tognelli et al., 2019;Miqueleiz et al., 2020). For this reason, species are grouped in genera in the analyses. However, different species belonging to one genus can show different habitat preferences (Anderson & Maldonado-Ocampo, 2011). For example, pencil catfishes of genus Trichomycterus include more than 170 species, with current high description rates and complex taxonomic history (Donin et al., 2020). Species of this genus are present in the entire South American continent, from sea level to [ 4,000 m of altitude, using different habitats and resources (Chará et al., 2006). These specific differences should be taken into consideration in the interpretation of our results. Whereas general habitat preferences of mountain species could be considered for most species of a genus (Anderson & Maldonado-Ocampo, 2011), some species present habitat adaptations that could differ from our results.
Most species of studied taxa have complex taxonomic relationships that are currently unclear, and researchers are not in full agreement about their taxonomic status (Lessmann et al., 2016). The uncertainty and complexity of Neotropical freshwater fish's taxonomy (Hubert & Renno, 2006;Olden et al., 2010) highlight the necessity to increase resources of taxonomists with relevant and high-quality skills in this area, as essential knowledge on ecological and conservation biology (Giangrande, 2003;Mace, 2004). Deficient taxonomic data (the Linnean shortfall) is a severe and general impediment on environmental studies in the Neotropical region (Collen et al., 2008;Anderson & Maldonado-Ocampo, 2011).

Conclusions
Our results show that variability in hydromorphological variables has relevance to the habitat requirements of tropical Andean fishes. River morphology configures the distribution and fish assemblage in these mountain ecosystems. Our study underlines the significance of hydromorphological variables on fish community structure in the area, whereas human alterations seem to be less significant. Our results suggested that the conservation of tropical Andean fishes requires maintaining habitat heterogeneity on a regional scale. These findings could assist in local management and conservation activities and be helpful in a global conservation assessment of Andean fishes, considering current threats over hydromorphology and general habitat features of tropical montane rivers. However, more data are needed to determine habitat fish relationships in the montane tropical Andes, allowing better strategies to conserve fish diversity in these aquatic ecosystems. There is a serious need to increase scientific knowledge and effective conservation of this fauna and other biodiversity living in freshwater montane tropical Andes and other upland regions. rivers, provided facilities, knowledge and assistance during fieldtrips in Ecuador. The Asociación para la Conservación de la Cuenca Amazónica (ACCA), kindly provided valuable information, collaboration and lodging during our fieldwork. This study was sponsored by the fellowship program III-B of the Ministry of External Affairs (MAEC-AECID), and the fellowships program for mobility of senior professors and researchers to foreign academic and research centres, Ministry of Education and Vocational Training of the Government of Spain (2019). Blanca Rios-Touma had the support of project AMB.BRT.20.01 (UDLA). The vice-rectory of research of the Universidad de la Amazonia provided logistic and economic resources for the achievement of this work.
Author contributions RM, BR-T, IT, HO, MP-R, SG and JMA-F conceived and designed the analyses, all authors collected the data, HO, IT, SG and JMA-F coordinated the identification of the specimens, RM and IT performed the analyses, RM, BR-T, AF-L, JMA-F and IT wrote the paper.
Funding Open Access funding provided thanks to the CRUE-CSIC agreement with Springer Nature. The Asociación de Amigos-University of Navarra awarded a doctoral grant to Ibon Tobes (2012Tobes ( -2014. This study was sponsored by the fellowship program III-B of the Ministry of External Affairs (MAEC-AECID), and the fellowships program for mobility of senior professors and researchers to foreign academic and research centres, Ministry of Education and Vocational Training of the Government of Spain (2019). Blanca Rios-Touma had the support of project AMB.BRT.20.01 (UDLA). The vice-rectory of research of the Universidad de la Amazonia provided logistic and economic resources for the achievement of this work.
Data availability Data used for this analysis are included on the supplementary material.

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