Introduction
Simuliids, known generally as “black flies” in English and “jejenes” in Spanish, are widely distributed in streams and rivers worldwide from sea level to altitudes of 5 000 m (Coscaron & Coscaron-Arias 1995, Adler et al. 2004, Adler & Crosskey 2017). A total of 2 204 extant species of black flies are known to science (Adler & Crosskey 2017) of which 74 have been recorded in Colombia (Moncada-Álvarez et al. 2017). Black flies are often the dominant form of secondary productivity in boreal lotic ecosystems (Malmqvist et al. 2004, Currie & Adler 2008). Due to their activity as filter feeders, these aquatic dipterans play a fundamental role in the transfer of energy from fine particulate organic matter top redators, represented by other macroinvertebrates and fish. They also are important in the supply of larger particles via their feces to a broad spectrum of deposit-feeding detritivores (Wallace & Webster 1996). Simuliid larvae typically are associated with well-oxygenated waters, although certain species can be found at low concentrations of dissolved oxygen (Welch & Tindell 1992). Consequently, the immature stages can be useful bioindicators of water quality (Martínez & Portillo 1999, Carle et al. 2015). However, few ecological studies of black flies have addressed taxonomic resolution to the species level (McCreadie & Adler 2008), and this has hindered a more accurate use of these organisms as bioindicators.
Most ecological studies of black flies have been conducted in Nearctic and Palearctic regions (Figueiró & Gil-Azevedo 2010, Figueiró et al. 2014). In the Neotropical region, ecological studies have mostly been carried out in Brazil, especially at low elevations (Hamada & McCreadie 1999, Hamada et al. 2002, McCreadie et al. 2005, Figueiró et al. 2006, 2008, 2012, Landeiro et al. 2009, Dos Santos et al. 2010, Couceiro et al. 2014, Docile et al. 2015), with some studies in Venezuela (Grillet & Barrera 1997). In Ecuadorian highland Andean streams the simuliids are one of the four most abundant families during both the wet and dry season (Jacobsen et al. 1997, Jacobsen & Encalada 1998), and one of the most abundant groups when riparian coverings are forests and pastures (Iñiguez-Armijos et al. 2016). These studies and other published on macroinvertebrates in Neotropical high mountain rivers only refer to the simuliids as a family without sufficient ecological information on the species.
In Colombia most simuliid literature is taxonomic (Adler & Crosskey 2017). Early ecological notes were recorded for Colombian black fly s pecies in the Sierra Nevada de Santa Marta (Lewis & Lee-Potter 1964). Bueno et al. (1979) published some ecological annotations for Simulium (Hemicnetha) muiscorium. An analysis of historical ecology of the Central and Eastern Cordilleras in Colombia revealed an elevational species gradient in which the lower areas correspond to the oldest regions with the most primitive species (Muñoz & Miranda-Esquivel 1997). Reports of the presence and abundance of black flies in the Fucha River (a stream of the Eastern Hills of Bogotá) were based on limited environmental observations (Restrepo & Rincón 2009), while a few ecological notes were presented for oviposition behavior of Simulium dinellii (Mantilla et al. 2015). Some studies have addressed the medical importance of simuliids as vectors of onchocerciasis in Colombia (Tidwell et al. 1980, Corredor et al. 1998), the transmission of which is currently eradicated from Colombia (WHO 2013). Other reports have dealt with the role of these dipterans in bird diseases (Lotta et al. 2016).
Simuliid species of the temperate zones, as well as species from lower elevations in Neotropical regions are different from those reported in the Colombian Andean region (Moncada-Álvarez et al. 2017). Therefore, one would expect that on a species level, ecological responses of high mountain river taxa in Colombia are different. A study of bioindication in Brazil showed that Simuliidae species were present in streams with intermediate urban pollution impacts, but absent in heavily impacted sites (Docile etal. 2015). These authors suggested that some species were associated with more pristine breeding sites, while others might be good bioindicators of moderately impacted streams. However, the Brazilian species are different from those found in Colombian waters, so it is necessary to study the native Colombian black fly fauna to know their responses to environmental variables.
The majority of the watercourses in the upper basin of the Bogotá River are used for agricultural, domestic, and industrial purposes. The Arzobispo, Chorro de Padilla (also known as San Francisco), La Vieja, and El Delirio (also known as La Osa) rivers have their upriver areas protected by the Empresa de Acueducto y Alcantarillado de Bogotá (EAAB), and these rivers eventually feed into the Bogotá River. Since, small, highland streams and rivers can influence limnological processes in lower reaches of the Andean basins (Allan et al. 2006), a long-term biomonitoring and environmental assessment of these watercourses is required to account for the impact imposed by anthropic practices. Given the potential value of black flies as bioindicators of water quality, knowledge about their biodiversity and relative abundance could facilitate the conservation of riverine ecosystems, essential for the maintenance of the urban lotic environment and water reserves for cities like Bogotá. Therefore, in the present study the relationship between water stream physical and chemical factors and the presence of black fly species was investigated in four streams of the Bogotá River watershed. This allows to explore the capability of simuliids as indicators of water quality.
Methods
Study area
The study area lies within the watershed of the Bogotá River, which has four sub-basins: Teusacá, Juan Amarillo (Salitre), Fucha, and Tunjuelito. The streams and rivers of these sub-basins run through the Eastern Hills around Bogota before entering the urban area, which lies in an extensive, high (La Osa) (Fig. 1). The list of environmental features of each sampled river is presented in Table. 1.
Sampling
Each river was sampled four times, twice during the rainy season (March and April) and twice during the beginning of the dry season (May and June) in 2012. A total of 16 samples per river were taken during the study. In each river a sampling section of 30 m including a combination of pools and riffles was randomly established. The same sampling section of each river was visited during the dry and rainy seasons. In each of these sections two people collected simuliid immature states during two hours of effort for a total of four hours per sampling. Abundance data were reported as number of individuals captured during these four hours. Forceps were used to manually collect simulliid larvae and pupae from available substrates (stones, rocks, pieces of wood, leaves, cloths, plastics) along the established sampling section of each river. No adults were collected because they were absent during the sampling periods.
Specimens were fixed in 70 % alcohol and transported to the Medical Entomology Laboratory at Universidad Nacional de Colombia (LEMUN) for identification using taxonomic keys by Muñoz et al. (1984), Muñoz & Coscarón (1999), Hamada & Grillet (2001), Shelley et al. (2004), and Coscarón & Coscarón-Arias (2007). Species identification was completed shortly after collection. The following physicochemical variables were measured in situ with a water sample from the stream surface: water temperature and pH were recorded with a potentiometer Schott® HandyLab; and nitrates and dissolved oxygen were assessed with Aqua-Merck® kits. In addition, water velocity was determined following the floating-object method reported by Hauer & Lamberti (2007).
Data analysis
The following ecological indices were calculated: the Shannon-Weaver diversity (natural logarithm), Simpson dominance, and Pielou’s equitability (i.e. evenness) indices. To examine the differences of the physical, chemical and biotic variables between the rivers, the nonparametric variance analysis of Kruskal-Wallis was applied to the data. Subsequently, post-hoc pairwise comparisons, based on the Tukey-Kramer test, were performed to establish the set of data that produced the differences. A principal component analysis (PCA) was then performed to inspect the relation among physicochemical variables coupled to a detrended correspondence analysis (DCA) to explore the spatial-temporal organization of species. Two nonmetric multidimensional scaling (NMDS) analyses were performed, one for species and another for environmental variables. For all these ordination procedures, data were square root-transformed due to small sample sizes and to reduce statistical biases. Non-parametric Spearman correlations were calculated to examine the statistical relationship between abiotic factors and taxon abundance using data of the 16 sampling cases. The similarity of the sampling sites was determined based on simuliid composition, using a Bray-Curtis similarity cluster analysis. All statistical analyses were conducted with the PAST program, version 3.18 (Hammer et al. 2001), and species accumulation curves were computed with the EstimateS program (Colwell & Elsensohn 2014), using the nonparametric estimators Chao1 (abundance) and Chao2 (incidence).
Results
A total of 1 492 larvae and pupae of five taxa across five species/species complexes were collected (Fig. 2): Simulium muiscorum (Bueno, Moncada & Muñoz), Simulium furcillatum (Wygodzinsky & Coscaron), Simulium ignescens (Wygodzinsky), Gigantodax cervicornis (Wygodzinsky), and the Gigantodax ortizi complex. The G. ortizi complex as treated here includes at least the following species: Gigantodax chilensis (Philippi), Gigantodax paramorum (Wygodzinsky & Coscaron), and Gigantodax multituberculatus (Wygodzinsky & Coscaron). Diagnostic characters for these three species of Gigantodax are equivocal, hence, until their taxonomy is clarified through molecular studies or capture of adult forms these three species were recognized as a single taxonomic entity.
Simulium muiscorum was the most abundant taxon in all rivers (68.6 % of the total relative abundance) excluding La Vieja, where it was absent in all sampled months but May. Gigantodax cervicornis was the least abundant species (0.13 %) and was found only in La Vieja in June (Fig. 2). According to the Kruskal-Wallis analysis the abundances of S. muiscorum and the G. ortizi complex were significantly different between rivers (H = 10.3, p = 0.02, H = 9.7, p = 0.02). The Tukey-Kramer test indicated that La Vieja River was the stream that accounted for this difference. There were no significant differences in the abundances of simuliid species between climatic seasons (H = 1.9, p = 0.2).
The obtained species accumulation curves, based on 16 samples and 5 taxa, were asymptotic for the Chao 1 and 2 estimators (Fig. 3), showing that all species of black flies present probably were collected. Thus, the estimators suggested that the sampling program had 100 % efficiency, both in numbers of individuals of each taxon (Chao 1) and as species richness (Chao 2).
The four rivers were acidic with high dissolved oxygen levels (> 5.9 mg L-1), cool temperatures (< 12 °C), moderately high mean stream velocities (1.3 - 1.4 m s-1), and appreciable nitrate concentrations (> 13.8 mg L-1) (Table 2). Dissolved oxygen and pH were higher in Chorro de Padilla. Arzobispo had larger amounts of nitrates and, together with La Vieja, had a slightly faster water velocity. El Delirio had the lowest pH and nitrates. Abiotic variables had values with no significant differences between the rivers, but during the dry period temperature was statistically higher (H = 5.1, p = 0.02) (11.5 °C in dryness vs. 10.4 ° C in rains) and water velocity faster (H = 6.1, p = 0.01) (1.65 m s-1 in dryness vs. 1.07 m s-1 in rains).
PCA analysis on abiotic variables showed that principal component (PC) axis 1 explained 41.98 % of the total variability and PC axis 2 explained 27.53 % (Fig. 4). The cumulative explanation for the two components was near 70 % of all variation among physicochemical data. The PC axis 1 was associated mainly with pH and dissolved oxygen and the PC axis 2 with nitrates. Arzobispo river in May and June was related with stream velocity and temperature, parameters that had high scores with axis 1 (0.41 and 0.46, respectively). Oxygen and pH were associated with El Delirio river in June (scores with axis 1: 0.55 and 0.57 respectively) and nitrates with La Vieja and Chorro de Padilla rivers in May (score with axis 2: 0.6). The other sites and samples did not show a clear relation with the physicochemical variables.
The DCA (Fig. 5) with a high eigenvalue of 0.532 on the first axis and a minor value of 0.147 on the second axis, indicated that La Vieja was characterized by the presence of the G. ortizi complex and that Arzobispo and Chorro de Padilla shared a predominance of S. muiscorum. The NMDS analysis for species with a stress value of 0.08, indicated that the composition of black flies in Arzobispo and Chorro de Padilla was similar (Fig. 6A). In contrast, the simuliid fauna of La Vieja was distinct, and in El Delirio simuliids also showed some degree of segregation (except in June). The NMDS analysis for environmental variables (with a stress value of 0.03) split data into two groups. One such group encompassed all April samplings with a handful of samplings across sites and months (Fig. 6B, right hand side cluster). The other group included all March samples in El Delirio and La Vieja; May samples in Chorro de Padilla, La Vieja and Arzobispo; and June samples in Arzobispo (Fig. 6B, left hand side cluster). All sites falling into the second group featured the highest levels of nitrates.
The abundance of S. muiscorum was negatively correlated with dissolved oxygen (r = -0.6, p = 0.02, n = 16), suggesting an increase in abundance in less oxygenated waters. In contrast, S. ignescens abundance correlated positively with oxygen (r = 0.6, p = 0.02, n = 16); S. ignescens tended to be slightly more abundant in more oxygenated waters. S. muiscorum also correlated negatively both with temperature (r = -0.5, p = 0.04, n = 16) and stream velocity (r = -0.6, p = 0.02, n = 16). Abundance of S. furcilatum was directly correlated with stream velocity (r = 0.5, p = 0.04, n = 16), but the correlation between abundance of the entire simuliid community and water velocity was negative (r = -0.8, p < 0.05, n = 16). The mean species number per river on each sampling date was 1.8 - 2.5. The low equitability (< 0.5) of the sampled rivers confirmed the dominance of some taxa (S. muiscorum and G. ortizi complex) over others. Shannon-Weaver diversity index had in general low values and was highest in El Delirio (0.86 nats) and lowest in Chorro de Padilla (0.14 nats) (Fig. 7A). In El Delirio, Simpson’s dominance index was low, and equitability was high, whereas in La Vieja, the converse was observed (Figs. 7B and 7C).
As indicated by the cluster in Fig. 8, Arzobispo and Chorro de Padilla were virtually identical in their simuliid composition; the composition of El Delirio resembled that of Arzobispo and Chorro de Padilla; and that of La Vieja was the most different.
Discussion
Ecological studies of black flies in much of the Neotropical Region are scarce, often addressing these dipterans at the family level and occasionally at the genus level (Prat et al. 2009). Ecological responses, however, are typically species-specific or even population-specific. The family Simuliidae is well known for its cryptic species (Adler et al. 2010, Conflitti et al. 2015). The inability to reliably distinguish nominal species of the G. ortizi complex potentially limits not only high values of diversity, but also it limits the possibility of establishing the responses of each species to environmental variations. We suggest the capture of adults or raising pupae to improve taxonomic identification. Other tools, such as molecular and cytogenetic techniques, and obtaining pharate adults should be used to correctly identify the species from the G. ortizi complex.
The rivers of the Eastern Hills of Bogotá have few simuliid species (a total of 5), with Simulium species being more abundant in the central (Arzobispo) and southern (Chorro de Padilla, El Delirio) rivers, and the Gigantodax species predominating to the north (La Vieja river). Also, a low number of simuliid taxa (7) have been recorded in watercourses of the municipality of La Calera (Muñoz et al. 1984) near the rivers of the current study. For other territories of Colombia (Departments of Cundinamarca and Tolima, located in the central Andean region of the country), Muñoz & Miranda-Esquivel (1997) recorded 38 species (15 species of Gigantodax and 23 species of Simulium). The mean number (1.8 - 2.5) of simuliid species in each of our rivers is lower than the typical number (2.3 - 4.7) in Nearctic and Neotropical streams and rivers (McCreadie et al. 2005). This lower number of species in the rivers of the Eastern Hills around Bogotá may be due mainly to the fact that the G. ortizi complex includes at least three cryptic species. It is probable that once the taxonomic difficulties of this complex are solved the species richness in the area will increase.
Dissolved oxygen, pH, and temperature are physicochemical parameters to which aquatic organisms are sensitive (Roldan & Ramírez 2008). Changes in these parameters in the four rivers of this study are mainly influenced by climatic variation. No industrial pollutants or household wastes are released into these rivers. Rain, however, can carry eroded sediments into the rivers, increasing substances such as nitrates (Roldán & Ramírez 2008). The high prevalence of Gigantodax in La Vieja could indicate that species of this taxon have some preference for the slightly more acidic waters of this river. This assumption, however, requires further testing, both in the field and in the laboratory. Also, it is likely that the characteristics of the substrate in La Vieja (cobbles and boulders) favor the establishment of species of this genus. This hypothesis, like the previous assumption, should be studied in future research.
Given the appropriate range of physicochemical factors, biotic factors (e.g., competition) and stochastic factors can also influence species assemblage at a given site (McCreadie et al. 1997). One of the factors influencing black fly diversity is elevation, which is associated with temperature (Srisuka etal. 2015, Ya’cob etal. 2016). The maximum difference in elevation among sampling sites was only 266 m, and it was coupled to a temperature difference of ca. 2 °C. This temperature difference was able to differentiate samples in May and June, during the transition from the rainy to the dry season, when waters were slightly warmer (as the Kruskal-Wallis analysis seems to confirm).
Water velocity influences the abundance of Diptera in Andean rivers (Donato 2008), and often it is correlated with simuliid distribution at various spatial scales (Orth & Maughan 1983, Morin et al. 1986, McCreadie & Colbo 1993). Although stream velocity was uniform across the rivers in this study (Table 2), the differences reflected in the standard deviation were sufficient to cause a reduction in the abundance of the total simuliid community when stream velocity was higher. Remarkably, S. furcilatum seemed to prefer swifter waters, which means that a greater speed seems to favor its abundance. Either way, water velocity can have an important effect on the abundance of these dipterans. S. muiscorum, which predominated in three of the four rivers (Chorro de Padilla, Arzobispo, and El Delirio), tends to have a clumped larval distribution. The individuals’ aggregations might allow the larvae to optimize filtering of particulate matter when water velocity increases (Eymann 1991), and perhaps a clumped population offers some protection from predators. However, the negative relationship between the total simuliid abundance and water velocity could mean that not all species are favored by the increase in water velocity. The influence of other physical and chemical variables was weak or not correlated with species abundance, perhaps due to the moderate degree of uniformity in the environmental conditions of all rivers. As previously mentioned, temporal variations in environmental conditions were restricted to temperature and stream velocity and did not significantly influence short-term faunal composition or species distribution.
Diversity indices were low. This is partly due to cryptic speciation in the Simuliidae (Adler et al. 2010; Conflitti et al. 2015). An example of this phenomenon is the G. ortizii complex, which could include three species of this genus in the studied rivers. Another reason for the low diversity values is that only Simuliid taxa were considered. The inclusion of other invertebrate groups would have increased diversity index values. All sampled streams had low equitability (that is, the distribution of individuals between populations was uneven), a clear indication of the trend of dominance of some taxa (mainly S. muiscorum). Research on the Simuliidae by Dos Santos etal. (2010) shows that the dominance in a given region (like the Guaravera river, Paraná State, Brazil) is common. In that basin the authors observed that some simuliids species were much more dominant than others. The Simpson index for El Delirio revealed that in this river species predominance was lowest and at the same time, it was the river with greater equitability, possibly because this is a second-order river (Table. 1), namely, with contributions from other streams, that could favor a slightly greater diversity. When observing the physicochemical variables of this river (Table 2), it is feasible to think that its colder, more uniform temperatures (from May-June), more oxygenated and nitrate-low environment could lead to a richer and more diverse assemblage of black flies. Furthermore, it can be hypothesized that the greater richness and diversity observed in El Delirio river is also a reflection of the intense competition for niches and habitats taking place in this river, and this prevents dominance of a single species.
The four rivers are geographically close to one another, and all are located in the same cool temperature zone, with similar elevation, flow, and slope. The rivers in the southern and central sectors of the study area (Arzobispo, Chorro de Padilla, and El Delirio) are more similar to one-another than to La Vieja, which is located in the north section, where the G. ortizi complex predominated. The clustering of Arzobispo and Chorro de Padilla probably reflects their physical proximity (2.4 k m). Despite the apparent environmental homogeneity in the studied rivers, the current results reveal that there is some variation in the black fly species assemblages among s i tes. Whether this is significant or not, this variation is clearly shown by the differences between black fly species correlations to different habitats or environmental variables. Thus, some simuliids species (e.g. S. ignescens) could be suitable indicators of good water quality, and some taxa (e.g. S. muiscorum) might be more tolerant to lower oxygen levels in slightly polluted rivers.
In conclusion, the results of this study show that the species of simuliids in the rivers of the Eastern Hills around Bogotá can be potential indicators of the environmental conditions of these lotic ecosystems. In this way, the presence of different black fly taxa of could indicate the conditions of each river. The present study is preliminary and exploratory, so it would be necessary in future studies to increase the number of river ecosystems and deepen the biota-environmental relationships, including many more variables, especially relating to those of sanitation, such as coliform bacteria, chemical and biological oxygen demands, phosphorus and electrical conductivity, among others.