CRISTIAN ROMÁN-P.
CÉSAR ROMÁN-VALENCIA
DONALD C. TAPHORN
Universidad del Valle, Facultad de Ciencias Exactas y Naturales, Departamento de Biología, Cali, Colombia. cromanpa94@gmail.com
Universidad del Quindío, Laboratorio de Ictiología, Apartado 2639, Armenia, Quindío, Colombia. ceroman@uniquindio.edu.co
1822 North Charles Street, Belleville, Illinois, 62221, USA. taphorn@gmail.com
ABSTRACT
Hemibrycon brevispini is a Neotropical characid fish endemic in La Venada Creek, a headwater tributary of the Quindío River of Colombia ( Cauca River drainage). It is mainly a diurnal insectivore with a diet dominated by benthic dipterans (Chironomidae, Simuliidae, Psychodidae, Culicidae, Calliphoridae, Dixidae and Muscidae), hymenopterans (Formicidae and Vespidae) and ephemeropterans (Baetidae), as well as allochthonous prey and items eaten accidentally. Microhabitats of mountain streams with lower water velocity tend to have more riparian vegetation and the associated terrestrial arthropods that are consumed by H. brevispini. It has three peaks in reproduction: December, April and August. Average fecundity was 776 mature oocytes per female.
Key words. Diet, reproduction, conservation, Neotropical fishes, natural history.
RESUMEN
Hemibrycon brevispini es un pez carácido neotropical endémico de la quebrada La Venada, un afluente del río Quindío en Colombia (cuenca del río Cauca). Esta especie es predominantemente insectívora diurna, con una dieta dominada por dípteros (Chironomidae, Simuliidae, Psychodidae, Culicidae, Calliphoridae, Dixidae y Muscidae), himenópteros (Formicidae y Vespidae) y efemerópteros (Baetidae), además de presas de origen alóctono y otras definidas como consumo accidental. Los hábitats con baja velocidad de agua sustentan mayor vegetación ribereña asociada a artrópodos terrestres, consumidos por H. brevispini. Su reproducción tiene tres picos: diciembre, abril y agosto. Su fecundidad promedio es de 776 oocitos por hembra.
Palabras clave. Dieta, reproducción, conservación, peces neotropicales, historia natural.
Recibido: 07/10/2013
Aceptado: 29/10/2014
INTRODUCTION
The genus Hemibrycon consists of fishes characterized by the presence of more than four teeth on the maxilla (in adults) (Eigenmann 1927, Román-Valencia et al. 2013). A phylogenetic analysis of Hemibrycon determined its monophyly based on four synapomorphies: ectopterygoids with widened ventral anterior projection, four to six times wider than posterior part; a red spot present in life on ventral margin of caudal peduncle; a postero-ventral projection on the pterotic and first infraorbitals gradually decreasing in width from posterior tip and located near posterior part of antorbital (Arcila-Mesa 2008). Hemibrycon brevispini Román-Valencia & Arcila-Mesa was described from La Venada and Quebrada Negra Creeks tributaries of the Santo Domingo River, Quindío River basin, upper Cauca, Andes of Colombia (Román-Valencia & Arcila-Mesa 2009).
Fourteen species of Hemibrycon have been described from the Cauca-Magdalena River Basin in Colombia, but there are few studies of their ecology that provide baseline information to determine their conservation status or provide guidelines for the management of many species that have relatively small geographical ranges and populations. In fact, habitat, extensive diet and reproductive data are only available for two Hemibrycon species from the Magdalena-Cauca River Basin: H. boquiae (Román-Valencia et al. 2008) and H. quindos (Román-Valencia & Botero 2006), but short notes have been published about of H. brevispini (Román-Valencia & Arcila-Mesa 2009), H. antioquiae, H. fasciatus and H. cardalensis (Román-Valencia et al. 2013), H. cairoense (Román-Valencia & Arcila-Mesa 2009), H. paez, H. raqueliae, H. virolinica and H. yacopiae (Román-Valencia & Arcila-Mesa 2010), H. palomae (Román-Valencia et al. 2010) and H rafaelense (Román-Valencia & Arcila-Mesa 2008). H. brevispini is endemic to La Venada and La Negra Creeks that are both tributaries of the Quindío River, in the upper Cauca River drainage (Román-Valencia & Arcila-Mesa 2009). Aspects of H. brevispini diet, reproduction and habitat were analyzed order to provide baseline information useful for conservation and management efforts of this endemic species and its habitat, especially considering the increasing impacts of hydropower and mining development in the Colombian Andes.
MATERIALS AND METHODS
Data collection and study area description. Sampling sites are distributed along the entire length of La Venada Creek, from its origin to its mouth where it discharges in to Quebrada Negra Creek. Fishes were collected in the middle and lower reaches of La Venada Creek (4o 26 47.4 N & 75 o 40 44.3 W, 1661 m.a.s.l. and 4o 26 54.9 N & 75o 40 48.8 W, 1307 m.a.s.l.), a tributary of Quebrada La Negra Creek, which in turn is a tributary of the Santo Domingo/Quindío/upper Cauca River system in the Andes of Colombia. Thus, La Venada Creek is a primary or secondary stream, in the hierarchical classification system of Allan (1995). Fish were captured on two days of each month from July 2011 to November 2012, using a 2 x 0.5 m seine net, with 5 mm mesh and a 2 m cast net with a 10 mm mesh between 0900 and 1300 hr, sampling much of the creek. This period included dry seasons (June-August and January-February) and wet seasons (March-May and September-December) (Fig. 1).
A total of 122 specimens of H. brevispini were collected and placed on ice to decrease the rate of enzymatic digestion of stomach contents, as recommended by Bowen (1996). Samples were dissected the same day as collection at the Ichthyology Laboratory of the Universidad of Quindío, Armenia , Colombia (IUQ). An incision was made along the ventrum, and the stomach, intestine and gonads were extracted. After dissection, specimens were fixed in formalin (10%) for 15 days and then in 70% ethanol and deposited in the fish collection (IUQ). Gonads were weighed using an analytical balance (Adventurer-Ohaus H226) with 0.0001 g precision and subsequently preserved in 70% ethanol. Mitutoyo digital calipers with 0.01 mm precision were used to measure standard and total length of the fish, stomach length, stomach width and intestine length.
Habitat measurements. A pH meter (Hanna HI 921ON) was used to measure pH, air temperature and water temperature. A digital oxygen meter (OX1196) was used to measure dissolved oxygen and saturation. Geographical coordinates were recorded from a GPS unit (Garmin eTrex 10).
Diet. Prey items in fish stomach contents were identified to the lowest taxonomic resolution possible (order, family, genus) (Borror et al. 1992, Roldàn 1996). Stomach contents were analyzed using numeric and frequency methods (Hyslop 1980, Hynes 1950) and a volumetric method (Capitoli, 1992, Pedley & Jones 1978). The Index of Relative Importance (IRI) (Oda & Parrish, 1981), proposed by Pinkas et al. (1971) was used to determine the importance of each food item.
IRI = %Fo (%N + %V)
Where, %V= percent volume, % Fo= observed frequency percent and %N= proportion of food type.
Prey item identification was done for individuals in different states of digestion or based on diagnostic structures, however unidentified remains or parts of organisms were not treated as items; these are listed, but not included in the statistical diet analyses. The emptiness coefficient (V) (Hyslop 1980) was also calculated to reveal the months included on the feeding period of the species,
V = n x 100
N
Where, n= number of empty stomachs, N= total number of stomachs examined.
A centered principal component analysis was made using stomach contents abundances (%N) and the proportion of the total volume of each prey item (called %PCA after Billy et al. [2000]). In this method, total row (of diet items encountered in an individual stomach) is equal to 1. For the analysis, the families were grouped within orders to enable a better explanation of the diet. This analysis reveals the pattern of dispersion of prey found in stomachs. Representation on two axes (Gabriel 1981) was done with the first factorial plane; the position of each item is equivalent to the position of a stomach containing 100% of the prey species. Each stomach is at the centroid of the prey items, with each prey species being given a weight equal to its proportion in the stomach. Computations and graphical displays were performed with the ADE-4 package (Thioulouse et al. 1997) running in R-software (R Development Core Team 2013). Correlation analysis between the %PCA scores and standard length was performed to examine patterns of prey consumption among fishes of different sizes.
Abundance distribution normality of food items was evaluated using the Kolmogorov-Smirnov test with a 5% of significance (α= 0.05). Based on this a Kruskal-Wallis was done to evaluate differences in prey abundance (numbers of individuals) among seasons, sex and maturation stage. A correlation analysis based on log-transformed data was applied to examine relationships among the variables: standard length (SL), total length (TL), intestine length (IL) it is in mm, stomach width (STW) in g, stomach length (STL) in mm, stomach weight (STWE) in g, gonad weight (GW) in g, fecundity (FE) in number of oocytes and total weight (TWE) in g. The Past 2.11 (Hammer et al. 2001) and R-software (R Development Core Team 2013) programs were used for statistical analyses.
ANOSIM was used to test dietary preferences among food items. The similarity matrix was generated with the transformed data (log[x+1]) for consumed prey using the Bray-Curtis similarity and the observed relationships were compared based on 9999 permutations. The value of R lies between -1 and 1, where 0 indicates that low and high similarities are perfectly mixed, thus there is no preference in dietary items.
Condition factor (K). A condition factor (K) was calculated and used to evaluate the populations condition (Wootton 1992, Vazzoler 1996, Bagenal & Tesch 1978),
K = Wt x 100
Lsb
Where, Wt (g) = total weight, Ls (mm) = standard length and b, la relation length weight.
Reproduction. To determine the reproductive season, temporal variation in the gonadosomatic index (GSI) was evaluated. GSI was equal to
GSI = Wo x 100
Wc
where Wc= Wt-Wo, and Wo (g) = gonad weight, Wt (g) = total weight, and Wc (g) = body weight (Vazzoler 1996). Spawning seasons were identified as peaks in mean GSI. Size at sexual maturity was determined using the graphic method of Sokal & Rohlf (1995) that identifies the maturation size as that for which 50% of the population is reproducing. Sex ratio was evaluated using chi-squared (X2) and proportion of males and females.
Fecundity. Fecundity was determined using the dry subsample method (Ricker, 1971), and absolute fecundity (Fa) was calculated using only mature females according to the formula
Fa = ∑ no
No
Where no= number of oocytes per female, and No= total number of females. Oocyte diameter was measured using millimetric graph paper by counting the number of oocytes fitting into 10 mm of the line on the paper and dividing by 10, and later calculating the average number of oocytes on a one-dimensional space for each ovary with oocytes.
RESULTS
Habitat. La Venada Creek is a primary stream in both its highest, and lower sections, with a width of 2- 3 m, and a depth 0.5- 1 m during both the rainy and dry seasons. Substrate is mostly rocky, with some sand and decomposing vegetation. For most of the length studied shore vegetation is not natural, consisting of white ginger (Hedychium coronarium), bamboo (Guadua angustifolia), coffee trees (Coffea arabica) and banana plantations (Musa spp.).
Ambient and surface water temperature ranged on average from 18˚C to 20.0°C during low-water season, and from 16.2°C to 20.5°C during the rainy season. Oxygen saturation was around 84% and dissolved oxygen was 6.0mg/L during the dry season, but values were in average higher during the rainy season: 92% and 7.9mg/L respectively. In the dry season pH was near 7.6 and during the wet season 7.1 (Table 1).
Hemibrycon brevispini was found in the middle and lower reaches of La Venada Creek along with Carlastyanax aurocaudatus (Eigenmann 1913), Astroblepus cf. cyclopus Humboldt 1805, Brycon henni Eigenmann 1913, Bryconamericus caucanus Eigenmann 1913, Cetopsorhamdia boquillae Eigenmann 1922, Chaetostoma cf. fischeri Steindachner 1879, Parodon caliensis Boulenger 1895, Poecilia caucana Steindachner 1880, Trichomycterus caliense Eigenmann 1918 and T. chapmani Eigenmann 1918.
Digestive tract morphology. The stomach of H. brevispini is longer (mean=12.7 mm, S.D. = 2.79) than wide (mean 7 mm, S.D. = 1.95) and is located in the anterior portion of the coelomic cavity, sometimes thickly covered with fat. Two pyloric caecae are present on the anterior part of the stomach. A significant, positive correlation was found between the intestine length and total length (r=0.67; p=0.017) as well as the intestine length and standard length (r=0.67; p=0,017; Fig. 6).
Diet. The feeding activity of this species is constant throughout the year (emptiness coefficient [V] = 0.82%). Prey abundance did not show normal distribution (p>0.05). The Kruskal-Wallis analysis revealed no difference in the abundance of the diet items between males vs. females or immature vs. adults (KW P= 0.25; df= 1; KW P= 0.31, df= 1 respectively). This was also the case when diet abundances were compared between wet and dry seasons (KW P= 0.98; df= 1). Thus season, sex and maturity are not determinant factors of the prey abundance in the diet of H. brevispini, and hence it appears that the items found in its diet are in constant supply all year round. Stones, feathers and nematodes found in the digestive tract of some individuals were considered occasional and accidental.
Stomach contents analysis revealed 41 total prey categories consumed by H. brevispini (Table 2, Fig. 2), with Diptera being the most frequently consumed (%N= 11.44; %FO= 16.96; %V= 15.01; IRI= 448.39), followed by Hymenoptera: Formicidae (%N= 15.84; %FO= 13.84; %V= 13.59; IRI= 407.21) and Ephemeroptera: Baetidae (%N= 16.09; %FO= 9.86; %V= 13.44; IRI= 291.19). Occasionally ingested organisms (classified as such based on their low consumption frequency) included Coleoptera, Araneae, Dyctioptera, and other allochthonous material. The degree of digestion of items found in stomach contents and the hour of capture allow us to infer that feeding is diurnal in H. brevispini, when they take mostly benthic organisms, some arthropods from the water column or that have fallen into the water from shoreline vegetation. If this species were a nocturnal feeder, the samples made during the day would not have found identifiable stomach contents that for the most part showed little effects of digestion (soft tissue still present) and indication that prey ingestion had occurred shortly before capture.
In the normalized principal component analysis (%PCA) (Fig. 2) of prey item abundances components one (17.18%), two (13.62%) and three (9.33%) explain only 40.14% of total variance. This is a consequence of the elevated heterogeneity of volume of prey consumed and large number of different prey items consumed by different individuals of H. brevispini. The principal component analysis recovered only a small percentage of this variability, leaving 59.06% of the variance unexplained. Thus, only a preliminary approximation of the trophic characteristics of this species is possible based on PCA, but it is evident that widely different amounts of a wide variety of prey ítems are eaten. The items Coleoptera, bird feathers, gravel, Diptera, Nematoda and Hemiptera accounted for much of the variation in individual diets, whereas Ephemeroptera, Hymenoptera and others items accounted for little variation and were distributed near the sample centroid. Diptera were very abundant in only a few stomachs and, overall, were not as common as the other items located near the centroid. Most of the individual stomachs in Figure 1 grouped near the origin are indicating that individuals employ a foraging strategy that exploits both dominant and rare prey items. The analysis indicates a broad trophic niche and relatively low between-individual variation in diets. Correlation analysis did not reveal an association between the first axis scores (% PCA1) and fish size (r=0.12, p=0.615), indicating no prey preferences in relation to fish length. ANOSIM revealed non-significant differences between the items found among individuals (R= -0.45, P=1).
Condition factor (K). In adults (female), the lowest K values were obtained from July 2011 (dry season), November-December 2011 and October 2012, these coincide with the wet season when fewer prey items were found in stomachs; low K values were obtained for adult males in September 2011 (wet season), February 2012 and July 2012 both of which coincide with the dry season. Maximum K values were obtained for females in October 2011, February and June 2012, during the wet and dry seasons, respectively; maximum K values for males were observed in July 2011 (dry season), October 2011 and 2012, both during the wet season. The lowest and maximum values in condition factor K contrast with the gonadosomatic index value GSI (Fig. 4) and cannot be interpreted as related to gonadal development. In immature specimens, variation was observed in condition factor values with the notable increase occurring between May and November of 2012 and July 2011, which corresponds to the wet and dry season; the variation in K values is considerable in both immature and adults. Comparison of maximum and minimum values by sex for both adults and immatures does not reveal great differences, but lowest values are from April 2012 and highest just one month later in May of 2012.
Reproduction. The gonosomatic index (GSI) (Fig. 4) showed high variation, up to one order of magnitude, among months of this study. Hemibrycon brevispini reproduction has three peaks during the year, with two large peaks in GSI in December 2011 and August 2012, and another smaller peak in April 2012.The two larger peaks coincide with the transition from wet to dry season (December) and from dry to wet season (August), and the smaller peak (April) is during the rainy season (Fig. 1). Females had higher GSI values than males throughout the year.
The high number of males present in the population is remarkable: 65.2% were males and 34.7% of the population individuals were females, giving a sex ratio of 1.9 with a predominance of males during the entire study period; significant differences exist as a result (X2= 9.26, df= 1, p= 0.05). For females, the size at sexual maturity is 76.3 mm SL and for males 68.7 mm SL (Fig. 5). Moreover, the size difference between the sexes is statistically significant (KW, P< 0.05, df= 1).
Fecundity. Average fecundity was 776 oocytes, and the mean diameter of mature oocytes was 0.85 mm (S.D. = 0.26). A non-significant low correlation value was found between fecundity and SL (r= 0.1, P> 0.05, Fig. 6). The mean weight of an oocyte was 3.4 x 10- 4 g (S.D. =1,37 x 10-6). Total body weight was significantly and positively correlated with gonad weight (r= 0.66, p= 0.019).
DISCUSSION
In streams of the upper Cauca River drainage, two species of Hemibrycon feed heavily on benthic insects such as Ephemeroptera, Odonata, and Trichoptera (Román-Valencia & Botero 2006, Román-Valencia et al. 2008), and this was also found for H. brevispini. Similar diets have been reported for other characid genera of the upper Cauca River drainage such as Creagrutus brevipinnis (Román-Valencia 1998), Roeboides dayi (Román-Valencia et al. 2003), Argopleura magdalenensis (Román-Valencia & Perdomo 2004), Carlastyanax aurocaudatus (Román-Valencia & Ruiz 2005), and Bryconamericus caucanus (Román-Valencia & Muñoz 2001a, Román-Valencia et al. 2008). It is commonly accepted that immature aquatic stages of insects are an abundant alimentary resource in Neotropical montane streams; however H. brevispini also consumed large amounts of ants (terrestrial Hymenoptera).
Pools and other habitats with low current velocities (0.23 to 0.67 m/s, mean= 0.35 m/s) in La Venada Creek inhabited by H. brevispini tended to have more riparian vegetation that probably supports ants and other terrestrial arthropods. Field observations of riparian vegetation in La Venada Creek and other similar Andean streams indicate that areas with lower water velocity are often wider than swift-water reaches and support denser riparian vegetation that in turn offers refuge and food for fishes.
Among allochthonous items found in their diet (14 of 42 categories), ants (Hymenoptera: Formicidae) were much more important than other categories (e.g., Vespidae, Diptera: Muscidae, Heteroptera, Auchenorrhynca, Chrysomelidae, Ptilodactylidae, Lepidoptera, Miriapoda: Diplopoda, Arachnida: Araneae, Seeds, Vegetative tissue, Feather, Pteridophyta, Dyctioptera).
Similar to the observed diet of H. brevispini, Román-Valencia et al. (2008) reported that Hemibrycon boquiae consumed a large proportion of Diptera larvae (Chironomidae, Simuliidae, Tipulidae, Ceratopogonidae and Muscidae): to evaluate why some items but not others are found in stomach contents it would be necessary to conduct prey preference and relative abundance studies, not stomach content occurrence as we present here. However, the ecological characteristics of Neotropical dipterans coincides with the preponderance of exploitation of these prey by H. brevispini. Diptera have pupae and larvae with aquatic or semi-aquatic habitats, in both running and quiet waters (Foote 1987, Brown 2001, Merritt et al. 2003, Courtney & Merritt 2008, Courtney et al. 2008). Among Diptera, Nematoceran families (especially Tipuloidea and Chironomoidea) are a preponderant component of aquatic communities, frequently eaten by primary consumers. Armitage et al. (1995) mention that Chironomidae are one of the most abundant families present in freshwater habitats, making them prone to capture.
We also found nematodes in H. brevispini, and interpret from their intact state that they are parasites and not prey; which coincides with findings reported for Bryconamericus caucanus (Román-Valencia & Muñoz 2001a). The relationship of intestine and body length is a general diet indicator. Herbivores have relatively long intestines compared to carnivores and omnivores (Wootton 1992, Kramer & Bryant 1995). The relatively short intestine length of H. brevispini indicates carnivory. Hemibrycon brevispini consumed mostly benthic arthropods (35 of the 45 food categories), and plant material was rare in stomachs.
The spawning period for H. brevispini differs greatly from patterns reported for other species in the genus. Hemibrycon boquiae spawns from July to September, the transition from dry to wet season (Román-Valencia et al. 2008). Hemibrycon quindos spawns from March to September in both wet and dry seasons (Román-Valencia & Botero 2006) (Fig. 1). Variation observed in the gonosomatic index of H. brevispini; could be caused by a seasonal life history (Winemiller 1989, 1992) influenced by a seasonal increase in food availability.
The low abundance of females during the sampling period could be due to bias imposed by males in sex ratios, since according some authors (Daiber 1977, Morse 1981, Petrie 1983, Clutton-Brock 1988, Goto et al. 1999, Goto et al. 2000, Trivers 1972) this could be a consequence of protandria, in which later maturation of females and dominance by males during juvenile life phases causes differential early mortality. Temperature affects are often cited as a principal cause of sex ration bias, but, the differential reproductive success between sexes or sampling bias have also been indicated as possible explanations. So our results may be a consequence of protandria and later maturation in females. Size at sexual maturity was larger for females than that of males. This may be a strategy of higher initial investment in somatic body weight which allows an increased investment in gonads at a later time. This trade-off strategy implies an extension of the time transcurred before entering the reproductive life phase, which in turn may permit higher mortality before reproduction. Although fitness in one hand is augmented by morphological characteristics, it is at the same time reduced by the inverse relationship between development time and the possibility of death.
Winemiller (1992) noted that medium-sized characids usually have a periodic life strategy characterized by late maturity, large numbers of eggs and low survival. Mean fecundity (776 oocytes) for H. brevispini is high when compared to other species of the genus: H. boquiae (376 oocytes) (Román-Valencia et al. 2008) and H. quindos (445 oocytes) (Román-Valencia & Botero 2006). When compared to other characids in the area, H. brevispini has lower fecundity than Bryconamericus caucanus (3759) (Román-Valencia & Muñoz 2001a), B. galvisi (1391 oocytes) (Román-Valencia & Muñoz 2001b), and Creagrutus brevipinnis (613 oocytes) (Román-Valencia 1998), but it has higher fecundity than Carlastyanax aurocaudatus (181 oocytes) (Román-Valencia & Ruiz-C. 2005).
The periodic life strategy is also evident in the larger average size of females with respect to males, a tendency shared with other species of characids common in the area. A larger size at sexual maturity was found for H. brevispini ( 76.3 mm SL for females and 68.7 mm SL for males) when compared to other species of the genus: H. boquiae ( 65 mm SL for females and 45 mm SL for males) (Román-Valencia et al. 2008) and H. quindos ( 53 mm SL for females and 50 mm SL for males) (Román-Valencia 1998); and when compared to other characids in the area: C. brevispinnis ( 40 mm SL) (Román-Valencia 1998), Bryconamericus caucanus ( 50 mm SL for female and 40 mm SL for male) (Román-Valencia et al. 2008), Carlastyanax aurocaudatus ( 35 mm SL for female and 40 mm SL for male) (Román-Valencia & Ruiz-C. 2005), except for B. galvisi (57.5- 89.9 mm SL for female and 61.3- 81.1 mm SL for males) (Román-Valencia & Muñoz 2001b).
Although we found no statistically significant seasonal differences in water quality parameters measured (pH, temperature, dissolved oxygen) between the wet and dry seasons (see also other data and parameters measured in Román-Valencia et al. 2005), this does not mean that changes in these parameters do not affect Hemibrycon brevispini. Physico-chemical parameters for habitats of H. boquiae (Román-Valencia et al. 2008), H. quindos (Román- Valencia & Botero 2006) and Hyphessobrycon poecilioides (García-Alzate & Román-Valencia 2008) also showed few significant seasonal differences.
Although physico-chemical parameters measured do not yet indicate decreased water quality in La Venada Creek (see also other data in Román-Valencia et al. 2005), threats not evaluated here, occurring near the study site from mining and current cultivated areas make this endemic species future uncertain. In this study area, other endemic species of Hemibrycon are present: H. boquiae, H. palomae and H. quindos, and are threatened in similar ways as H. brevispini. Studies of trophic and reproductive ecology of this species will provide a useful baseline for future impact studies.
ACKNOWLEDGMENTS
We thank K. Winemiller and three anonymous reviewers for valuable criticism and suggestions on the manuscript. We thank Biology students of the University of Quindío, Armenia (IUQ) for help with sampling.
REFERENCES
1. Allan, J.D. 1995. Stream ecology: structure and function of running waters. Chapman & Hall, New York, USA.
2. Arcila-Mesa, D.K. 2008. Análisis filogenético y biogeográfico de las especies pertenecientes al género Hemibrycon (Ostariophysi, Characiformes). Trabajo de grado, Programa de Biología, Universidad del Quindio, Armenia, 240p.
3. Armitage, P.D., P.S. Cranston & L.C.V. Pinder (eds). 1995. The Chironomidae: Biology and Ecology of Non-Biting Midges. Chapman and Hall, London. XII+572 pp.
4. Bagenal, T.B. & F.W. Tesch. 1978. Age and growth. In: T. Begenal (ed.). Methods for assessment of fish production in fresh waters, 3rd Edn. IBP Handbook No. 3, Blackwell Science Publications, Oxford: 101-136.
5. Billy, V. de C., S. Dolec & D. Chessel.2000. Biplot presentation of diet composition data: an alternative for fish stomach contents analysis. Journal Fish Biology 56:961-973
6. Borror, D.J., C.A. Triplehorn & N.F. Johnson. 1992. An Introduction to the study of Insects. Six Edition. Saunders College Publishing. U.S.A. 187-202.
7. Bowen, S.H. 1996. Quantitative description of the diet. In : Fisheries techniques p. 513-522. In: B. R. Murphy & D. W. Willis, editors. 2nd edition. American Fisheries Society, Bethesda, Maryland.
8. Brown, B.V. 2001. Flies, gnats, and mosquitoes. Pp. 815-826. In: S.A. Levin (ed). Encylopedia of Biodiversity. Academic Press, London.
9. Capitoli, R.R. 1992. Métodos para estimar volúmenes do conteudo alimentar de peixes e macroinvertebrados. Atlántica, Rio Grande 4:117-120
10. Clutton-Brock, T.H. 1988. Reproductive success. pp 472-485. In: T. H. Clutton-Brock (ed.). Reproductive success. University of Chicago Press, Chicago.
11. Courtney, G.W. & R.W. Merritt. 2008. Aquatic Diptera: Part one: larvae of aquatic Diptera. Pp. 687-722. In: R.W. Merritt, K.W. Cummins & M.B. Berg (eds). An Introduction to the Aquatic Insects of North America. Fourth Edition. Kendall/Hunt Publishing, Dubuque, Iowa. 1158 pp.
12. Courtney, G.W., R.W. Merritt, K.W. Cummins, M.B. Berg, D.W. Webb, & B.A. Foote. 2008. Ecological and distributional data for larval aquatic Diptera. Pp. 747771. In: R.W. Merritt, K.W. Cummins, & M.B. Berg (eds). An Introduction to the Aquatic Insects of North America. Fourth Edition. Kendall/Hunt Publishing, Dubuque, Iowa. 1158 pp.
13. Daiber, F.C. 1977. Salt-marsh animals: distributions related to tidal flooding, salinity and vegetation. In: Ecosystems of the world. I wet coastal ecosystems, ed. V. J. Chapman. Elsevier Scientific Publishing Company, Amsterdam, pp. 79-108.
14. Eigenmann, C.H. 1927. The American Characidae. Memoirs Museum Comparative Zoology XLIII: 311-428.
15. Foote, B.A. (coordinator). 1987. Order Diptera. Pp. 690-915. In: F. W. Stehr (ed). Immature Insects. Volume 2. Kendall/Hunt Publishing, Dubuque, Iowa. 975 pp.
16. Gabriel, K.R. 1981. Biplot display of multivariate matrices for inspection of data and diagnosis. In: Interpreting multivariate data, p.147-174. Barnett V Ed. Wiley, New York.
17. García-Alzate, C. A. & C. Román-Valencia. 2008. Biología alimentaria y reproductiva de Hyphessobrycon poecilioides (Pisces: Characidae) en la cuenca del río La Vieja, Alto Cauca, Colombia Revista Museo Argentino Ciencias Naturales, n.s. 10: 17-27.
18. Goto, R, T. Kayaba, S. Adachi & K. Yamauchi. 2000. Effects of temperature on sex determination in the marbled sole (Limanda yokohamae). Fisheries Science 66: 400.
19. Goto, R., T. Mori, K. Kawamata, T. Matsubara, S. Mizuno, S. Adachi & K. Yamauchi. 1999. Effects of temperature on gonadal sex determination in barfin flounder (Verasper moseri). Fisheries Science 65: 884-887.
20. Hammer, Ø., D.A.T. Harper & P.D. Ryan. 2001. PAST: Paleontological statistics software package for education and data analysis. Palaeontologia Electronica 4(1): 9p.
21. Hynes, H.B.N. 1950. The food of fresh water Sticklebacks (Gasterosteus aculeatus and Pygosteus pungitius), with a review of methods used in studies of the food of fiches. Journal Animal Ecology 19:36-58.
22. Hyslop, E.J. 1980. Stomach contents analysis a review and methods and their application. Journal Fish Biology 17: 411-429.
23. Kramer, D. & M. Bryant. 1995. Intestine length in the fishes of a tropical stream: 2. Relationships to diet: the long and short of a convoluted issue. Environmental Biology Fishes 42: 129-141.
24. Merritt, R.W., G.W. Courtney & J.B. Keiper. 2003. Diptera. Pp. 324-340. In: V.H. Resh & R.T. Cardé (eds). Encyclopedia of Insects. Academic Press, London.
25. Morse, W.W. 1981. Reproduction of the summer flounder, (Paralichthys dentatus)(L.) Journal Fish Biology 19: 189-203.
26. Oda, D.K. & J.D. Parrish. 1981. Ecology of commercial snappers and groupers introduced to Hawaiian reefs. Proc. Fourth Intern. Coral reef Sym. 1: 59-67
27. Pedley, R.B. & J.W. Jones. 1978. The comparative feeding behaviour of Brown trout Salmo trutta L. and Atlántico salmon Salmo salar L. In: Llyndwytwch, Wales . Journal Fish Biology 12: 253-256.
28. Petrie, M. 1983. Mate choice in role-reversed species. Pag.:167-179. In: P. Bateson (ed.) Mate choice. Cambridge University Press, Cambridge.
29. Pinkas, L., M.S. Oliphant & L.R. Iverson. 1971. Food habits of albacore, bluefin tuna, and bonito in California waters. Fishery Bulletin. 152: 1-105.
30. R Development Core Team 2013. R: A language and environment for statistical computing. R foundation for statistical computing, Vienna, Austria .
31. Ricker, E. 1971. Methods for assessment of fish production in freshwater IBP. Handbook 3. Blackwell Sci. Publ., Oxford and Edinburgh. F. A.Davis, Philadelphia. xiv + 348 p.
32. Roldán, P.G. 1996. Guía para el estudio de los macro-invertebrados acuáticos del Departamento de Antioquía. FEN-Colombia, Colciencias-Universidad de Antioquia, Medellín. 217p.
33. Román-Valencia, C. 1998. Alimentación y reproducción de Creagrutus brevipinnis (Pisces: Characidae) en alto Cauca, Colombia. Revista de Biología Tropical 46: 783-789.
34. Román-Valencia, C. & A. Muñoz. 2001a. Ecologìa trófica y reproductiva de Bryconamericus caucanus (Pisces: Characidae). Bolletino Museum Regionalli Science Naturalli, Torino 18: 459-467.
35. Román-Valencia, C. & A. Muñoz, A. 2001b. Alimentación y reproducción de Bryconamericus galvisi (Pisces: Characidae) en Alto Putumayo, Amazonia Colombiana. Boletin Ecotropica: Ecosistemas Tropicales 35: 37-50.
36. Román-Valencia, C., A. Botero & R. Ruiz-C. 2003. Trophic and reproductive ecology of Roeboides dayi (Teleostei: Characidae) from upper Rio Cauca, Colombia. Bolletino Museum Regionalli Science Naturalli, Torino 20: 487-496.
37. Román-Valencia, C. & A. Perdomo. 2004. Ecología trófica y reproductiva de Argopleura magdalenensis (Pisces: Characidae) en cuenca alta de los ríos Cauca y Magdalena, Colombia. Revista Museo Argentino de Ciencias Naturales n. s. 6:175-182.
38. Román-Valencia, C., J.G. Cadavid, J.A. Vanegas & D.K. Arcila. 2005. Análisis de algunas variables físicas, químicas y biológicas en tres quebradas de la cuenca Alta del Rio Cauca, Colombia. Revista de Investigaciones, Universidad del Quindío, 15: 83-93
39. Román-Valencia, C. & A. Botero. 2006. Trophic and reproductive ecology of a species of Hemibrycon (Pisces: Characidae) in Tinajas creek, Quindio River drainage,upper Cauca basin, Colombia . Revista Museo Argentino de Ciencias Naturales n.s. 8:1-8.
40. Román-Valencia, C, R.I. Ruiz-C. & A. Giraldo. 2008. Dieta y reproducción de dos especies sintópicas: Hemibrycon boquiae y Bryconamericus caucanus (Pisces: Characidae) en la quebrada Boquía, Río Quindío, Alto Cauca, Colombia. Revista Museo Argentino de Ciencias Naturales n.s. 10: 55-62.
41. Román-Valencia, C. & D.K. Arcila-Mesa. 2008. Hemibrycon rafaelense n.sp. (Characiformes, Characidae), a new species from the upper Cauca River, with keys to Colombian species. Animal Biodiversity and Conservation 31.1:67-75.
42. Román-Valencia, C. & D.K. Arcila-Mesa. 2009.Two new species of Hemibrycon (Characiformes, Characidae) from the Magdalena River, Colombia. Animal Biodiversity and Conservation 32(2): 77-87.
43. Román-Valencia, C. & D.K. Arcila-Mesa. 2010. Five new species of Hemibrycon (Characiformes: Characidae) from the Rio Magdalena basin, Colombia . Revista de Biología Tropical 58: 339-356.
44. Román-Valencia, C., C.A. García-Alzate, R.I. Ruiz-C. & D.C. Taphorn. 2010. A new species of Hemibrycon (Teleostei: Characiformes: Characidae) from the Roble River, Alto Cauca, Colombia, with a key to species known from the Magdalena-Cauca River Basin. Vertebrate Zoology 60:99-105.
45. Román-Valencia, C., R.I. Ruiz-C., D.C. Taphorn, N. Mancera-Rodríguez & C.A. García-Alzate. 2013. Three new species of Hemibrycon (Characiformes: Characidae) from the Magdalena River Basin, Colombia. Revista de Biología Tropical 61: 1365-1387.
46. Sokal, R.R. & Rohlf. 1995. F.J.Biometry. Third edition. W. H. Freeman and Co. New York.
47. Thioulouse, J., D. Chessel, S. Doledec & J.M. Olivier. 1997. ADE-4: a multivariate analysis and graphical display software. Stat Comput 7: 45-83.
48. Trivers, R.L. 1972. Parental investment and sexual selection. Pages 136- 179 in B. G. Campbell, ed. sexual selection and the descent of man, Aldine, Chicago.
49. Vazzoler, A.E. de A. 1996. Biología de reprodução de peixes teleósteos: teoría y práctica. Universidad Estatal do Brasil. Maringa, EDUEM; São Paulo: SBI. 169p.
50. Winemiller, K.O. 1989. Patterns of variation in life history among South American fishes in seasonal environments. Oecologia 81: 225-241.
51. Winemiller, K.O. & K.A. Rose. 1992. Patterns of life-history diversification in North American fishes: implications for population regulation. Canada Journal Fishery Aquatic Sciences 49: 2196-2218.