Files > Volume 6 > Vol 6 No 4 2021 > Vol 6 No 2 2021
INVESTIGATION / RESEARCH
Size variation of morphological traits in Bosmina freyi and its relation with environmental variables in a tropical eutrophic reservoir
*Yury Catalina López-Cardona1, Edison Parra-García12, Jaime Palacio-Baena3, Silvia Lucía Villabona-González4
Available from: http://dx.doi.org/10.21931/RB/2021.06.02.16
We assessed the size variation of morphological traits in Bosmina freyi regarding changes in environmental variables, the biomass of invertebrate predators, and algal food availability in two depths of the photic zone, the riverine zone, and near the dam zone (lacustrine zone) in The Riogrande II reservoir. In 200 individuals of B. freyi, using the software TpsDig2 we measured the body size, mucron and antennule lengths, and the antennule aperture percentage. Using the Mann-Whitney U test, we assessed the differences between these traits considering the zones and the photic depths; however, we used a canonical discriminant analysis with morphologic traits and environmental variables. Measured morphological traits showed a heterogeneous distribution between sampled zones and depths (p < 0.05). The highest values mucron and antennule lengths and the smallest antennule aperture angle were observed on small body size individuals, associated with physical, chemical, and biological characteristics in the riverine zone and the subsurface. Size structure distribution in B. freyi was related to changes in water temperature, trophic state, depredation, availability, and quality of food, of which implications related to the zooplankton community structure, predator-prey relations, and energy flow in the reservoir.
Keywords: Cladocera, predation, size structure, trophic state.
Cladocera is a zooplankton assemblage that occupies an intermediate position in the aquatic food chain 1 because they ingest protozoa 2, algae, bacteria, and detritus lentic systems, and simultaneously they constitute one of the chosen prey for invertebrates and vertebrates 3.
Some species of this group have a wide variation in size structure due to genetic factors 4 and environmental factors as temperature 5, food availability 6, pH 7, and acidification1.
Bosmina genera is a group of small cladoceran filter feeders 8 with variations in body size, shape, and length of antennules and mucron 1,5,9.
According to Masson S (2004) 10 in Canadian lakes, increasing chlorophyll-a concentrations favors small individuals' predominance. Other authors detail that some changes in Bosmina size are because the predation pressure by Copepoda, Rotifera 11, 12 and macroinvertebrates of Chaoborus genera 1, 8, 9, 13. 4 indicate that predation by invertebrates induces an increase of size and the development of defense structures as the elongation of mucron and antennules in Bosmina.
Bosmina freyi14 contributes more than 60% of the total zooplankton biomass in The Riogrande II reservoir 15 due to its broad trophic plasticity 16 and resistance to high turbidity conditions17. This species was reported for the first time in 2004 in South-America 18, for Colombia in 2010 (Elmoor-Loureiro, com.pers.), and consequently, there is a lack of information about the structure of morphological traits.
We hypothesized that if the size distribution of morphological traits of Bosmina freyi in the vertical and longitudinal axis in The Riogrande II reservoir changes concerning environmental variables and the presence of invertebrate predators, so in the reservoir zone with the highest trophic state and warmer water, individuals of small size will dominate and in zones with high invertebrate predators biomass, individuals will have a mucron and antennule with a longer length and an antennule closed-angle.
The Riogrande II reservoir is located in the northeastern of Colombia in The Río Grande II basin at 2200 m.a.s.l (6º 33' 50"- 6º 28' 07" N and 75º 32' 30"- 75º 26' 10" O). Has a total volume up to landfill level (2270 m.a.s.l) of 240 million of m3 and an average area of 1214 ha. This waterbody receives water from the Grande and Chico rivers and the Las Ánimas creek 19. This system is mainly used for hydroelectric generation and potable water supply 20. The Riogrande II reservoir has higher turbidity and a greater degree of eutrophication at riverine zones, mainly in the Chico River zone, due to sediment washing process in the basin and the domestic, industrial, agricultural, and livestock wastewater discharge from the closer cities (Entrerríos and San Pedro de los Milagros) located in the two principal rivers 21.
This research was conducted as part of the “Estudio de la problemática ambiental de tres embalses de Empresas Públicas de Medellín para la gestión integral y adecuada del recurso hídrico”, we carried out five sampled areas in The Riogrande II reservoir, in the subsurface and the photic zone edge (Figure 1) in September 2011. The highest dominance of Bosmina freyi happened in this period 15.
Figure 1. Study area and location of sampling stations
Physical, chemical and biologic variables data information (chlorophyll-a and plankton biomass) were took from the “Estudio de la problemática ambiental de tres embalses de Empresas Públicas de Medellín para la gestión integral y adecuada del recurso hídrico” 15, 22, 23, 24.
Methods for zooplankton sample collection and the estimation of Carlson index of trophic state (TSI) modified by 25 are described in 14.
We took 20 specimens of B. freyi (n = 200 individuals), and we used a photographic camera adapted to an optical microscope and were individually photographed to determine morphological traits using TpsDig2 (Figure 2). Lengths were measured in µm and the antennule angle in the opening percentage, assuming that a right angle (90°) has an aperture of 100%.
Figure 2. Body length (BL), mucron length (ML), antennule total length (AL) and antennule angle (AA) of Bosmina freyi (Modified from 5)
We used a Mann–Whitney U test because data did not fulfill the normality of residuals and homogeneity of variances. The last is to assess the significance of the differences between individuals' morphological traits from the zones and depths.
We also made a Generalized Discriminant Analysis using the measured morphological traits and the environmental variables. We calculated the scores and canonical vectors for graphic representation of the essential canonical discriminant functions in terms of explained variance, allowing a simple interpretation in the canonical space of differentiation among morphological traits of B. freyi and the environmental matrix 26. Data analyses were made using R Wizard V 2.3 Software 27.
The riverine zone had the highest temperature degree, greater eutrophication, and higher phytoplankton biomass (chlorophyll-a), rotifers (Asplanchna girodi De Geurne, 1888), and copepods predators (Mesocyclops longisetus Thiébaud, 1912; Metacyclops leptopus Kiefer, 1927; and M. mendocinus Wierzejski, 1892) biomass. Meanwhile, the lacustrine zone recorded the highest transparency, dissolved oxygen concentration, and Cryptomonas sp. biomass (Table 1).
B. freyi body length varied over the range of 202.64 and 457.02 µm (average = 294.7 µm), mucron length between 14.87 and 29.94 µm (average = 22.47 µm), antennule length between 94.21 and 170.25 µm (average = 134.11 µm) and the antennule opening angle ranged between 72.48 and 88.21% (average = 81.05%).
According to the Mann-Whitney U test, mucron and antennule lengths and the opening angle were significantly different among sampling zones and depths (Figure 3, a & b). Meanwhile, body length was significantly different among depths; it exhibited a significantly higher value in the specimens from the photic zone boundary (Figure 3 b). While mucron and antennule lengths were higher in individuals from the riverine zone and the subsurface concerning those from the lacustrine zone and the photic zone boundary, antennule aperture showed an opposite behavior.
Figure 3. Morphological traits distribution BL, ML, AL, and AA into two zones (a) and depths (b) in the photic zone of The Riogrande II reservoir, September 2011.
Table 1. Physical, chemical, and biological variables in The Riogrande II reservoir's photic zone, September 2011. C.V (%) = Coefficient of Variation.
A variance of 100% explains canonical discriminant analysis among zones and depths in the first axis. In both cases, classification and cross-validation values were 100%, clearly discriminating between riverine zones and lacustrine zone, as well as the subsurface and the photic zone boundary.
While individuals of B. freyi with higher mucron and antennule lengths are related with the riverine zone (warmer waters, as well as a greater degree of eutrophication and higher phytoplankton and predatory copepods biomass), in the lacustrine zone (more oxygenated and clear waters), specimens showed a tendency to present higher body size and an antennule more separated from the body (Figure 4 a & 5).
Depths were discriminated mainly by the presence of individuals with a higher mucron and antennule lengths in the subsurface, related to a higher temperature and greater dissolved oxygen availability. By contrast, individuals with higher body size and antennule aperture were mainly associated with the photic zone boundary, with higher predatory copepods and rotifers biomass, as well as phytoplankton and Cryptomonas sp. biomass (Figure 4 b & 5).
Figure 4. Canonical Discriminant Analysis applied to B. freyi morphological traits, environmental variables and phytoplankton, and zooplankton biomass. Linear regression model for zones (a) and depths (b) in the photic zone of The Riogrande II reservoir, September 2011.
Distribution map of size variation of morphological traits in Bosmina freyi
Figure 5. Distribution map of size variation of morphological traits in Bosmina freyi in The Riogrande II reservoir.
Melão, M. G. (1999)28 mentioned that temperature increase is a factor that directly affect organism's life cycle and population attributes, causing a time decrease in the organism's growth and an increase in population growth rates, which favors large populations with small size individuals 29, such as we observed in this study in the riverine zone and in the reservoir subsurface, where occurred the highest number of individuals with youth traits and higher mucron and antennule lengths 30.
Higher values of Carlson Index (TSI) and chlorophyll-a concentration in the riverine zone also favor small individuals due to a continuous energy source to maintain high metabolic rates 31. So, B. freyi population size structure in the reservoir seems to behave similarly to the general pattern of cladoceran size structure, in which smaller species are correlated to conditions of a more significant trophic state 32.
It has been reported that biotic variables have an indirect influence on mucron and antennule sizes and a direct effect on spatial and temporal population's dynamic, affecting the presence of predators and establishing food restrictions 33. This situation implies changes in size structure and the lengths of species attributes 1, as evidenced by this study.
The highest biomass of predatory copepods and rotifers in the riverine zone matched with has been reported in some studies about the influence of predation by invertebrates on morphological traits of Bosmina genera 9, 34, 35. Significantly 36, claim that copepods can consume Bosminidae species with sizes less than 400 µm. Nevertheless, the highest mucron and antennule size could reduce predation rates of small individuals of B. freyi by copepods, as it happens with other cladoceran species 37.
Antennule aperture may also be a defense mechanism against predators, and it depends on the predator attack mechanism. Some authors such as 38, 39 consider that a curved antennule (lower aperture) as that of the individuals detected in the riverine zone, protects the ventral part of cladocerans from some predators of genus Mesocyclops and especially from Asplanchna girodi, a rotifer species that dominated the biomass of rotifers in The Riogrande II reservoir during the study made by 15.
The highest biomass of the smaller organisms of B. freyi, of the predatory copepods and rotifers, may also be due to the highest primary production in the riverine zone 40 because the highest availability of resources may determine its presence in the longitudinal axis of the reservoir through the tributaries, especially the Chico River 39.
In the vertical axis, the size distribution of B. freyi structures could be attributed to vertical migration, common in the zooplanktonic community, particularly in Bosmina genera 41. In The Riogrande II reservoir 15 found a gradient in the vertical distribution of the zooplankton biomass. While smaller species seem to prefer superficial layers of the reservoir, the largest ones tend to be located in the deepest layers of the photic zone, probably as a strategy to avoid predation and selective removal of bigger species subsurface. This could resemble the vertical distribution of sizes in B. freyi population in which the individuals of larger size and smaller antennule and mucron length were mainly concentrated in the deepest layer of the photic zone.
Predatory copepods and rotifers agglomerated in this layer, and probably, as the individuals of B. freyi with a larger size (higher than 400 µm), they avoid predation/or are selectively removed from the subsurface 15.
The highest abundance of more significant individuals and smaller mucron and antennule lengths tended to coincide with the increase of Cryptomonas sp biomass in the zone near the dam (lacustrine zone), the photic zone boundary. The last implies that maybe these genera of algae constitute an essential resource for adults of B. freyi 30 in the reservoir 24 said that nanoplanktonic and mixotrophic species of Cryptophyta are a persistent component of zooplankton diet from The Riogrande II reservoir; apparently, it favors its growth and reproduction rates because it provides enough organic carbon, nitrates, and phosphates. Furthermore, it is very feasible that B. freyi complements its diet with detritus 42 and that alimentary strategy becomes more critical on small individuals of B. freyi together with the dissolved substances and organisms of picoplankton ingestion (< 2 µm).
Results showed that the distribution in the size structure in B. freyi is related to changes in biotic and abiotic variables. So, high temperatures and a greater nutrients concentration favored the growth of small individuals' populations, mainly in the reservoir's riverine zones.
Small sizes individuals seem to be protected from predation by other zooplankters as rotifers and copepods by a high mucron and antennule length, as well as by a closed antennule. By contrast, individuals from B. freyi with a bigger size and adult features just like predatory rotifers and copepods seem to avoid predation pressure, probably from visual predators as vertebrates, remaining mainly in the deepest zone of the photic layer. Furthermore, these cladocerans coexist with the highest densities of Cryptomonas sp., one of the most nutritious resources for Bosmina in the reservoir.
Consequently, the distribution of the size structure of B. freyi affects the structure of the zooplanktonic community, the predator-prey trophic relationships, and the flow of energy and matter in the reservoir. Nevertheless, factors do not involve here, as competence and seasonality may also be determinants and affect the population attributes of B. freyi.
Thanks to the program “Estudio de la problemática ambiental de tres embalses de Empresas Públicas de Medellín para la gestión integral y adecuada del recurso hídrico”, for the provision of information; and to the research groups Limnobase-Biotamar and GAIA from the Universidad de Antioquia for their support. Thanks to Isabel Cristina Gil for the paper translation.
1. Labaj, A.L., Korosi, J.B., Kurek, J., Jeziorski, A., Keller, W.B. and Smol, J.P. (2016). Response of Bosmina size structure to the acidification and recovery of lakes near Sudburdy, Canada. Journal of. Limnology, 75(s2): 22–29. http://doi.org/10.4081/jlimnol.2016.1183
2. Demott, W.R. (1982). Feeding selectivities and relative ingestion rates of Daphnia and Bosmina. Limnology Oceanography, 27(3): 518–527.
3. De Bernardi, R. and Peters, R. (1987). Why Daphnia? In Peters, R. and De Bernardi, R., ed. Daphnia. Memorie Dell’Istituto Italiano Di Idrobiologia, 45, pp. 1-9.
4. Black, R.W. (1980). The Nature and Causes of Cyclomorphosis in a Species of the Bosmina longirostris Complex. Ecology, 61(5): 1122–1132. http://www.jstor.org/stable/1936832.
5. Kappes, H. and Sinsch, U. (2002). Morphological variation in Bosmina longirostris (OF Muller, 1785) (Crustacea: Cladocera): consequence of cyclomorphosis or indication of cryptic species? Journal of Zoological Systematics and Evolutionary Research, 40, 113–122.
6. Kappes, H. and Sinsch, U. (2002b). Temperature and predator induced phenotypic plasticity in Bosmina cornuta and B. pellucida (Crustacea: Cladocera): Freshwater Biology, 47, 1944–1955.
7. Korosi, J.B., Paterson, A.M., Desellas, A.M. and Smol, J.P. (2008). Linking mean body size of pelagic Cladocera to environmental variables in Precambrian Shield lakes: a paleolimnological approach. Journal of Limnology, 67, 22-34.
8. Korosi, J.B., Kurek, J. and Smol, J. (2013). A review on utilizing Bosmina size structure archived in lake sediments to infer historic shifts in predation regimes. Journal of Plankton Research, 35 (2): 444–460. http://doi.org/10.1093/plankt/fbt007
9. Chang, K. and Hanazato, T. (2003). Seasonal and reciprocal succession and cyclomorphosis of two Bosmina species (Cladocera, Crustacea) coexisting in a lake : their relationship with invertebrate predators. Journal of Plankton Research, 25 (2): 141–150. Retrieved from: http://plankt.oxfordjournals.org/
10. Masson, S., Pinel-Alloul, B. and Dutilleul, P. (2004). Spatial heterogeneity of zooplankton biomass and size structure in southern Québec lakes: variation among lakes and within lake among epi-, meta- and hypolimnion strata. Journal of Plankton Research, 26, 1441–14458.
11. Guiset, A. (1977). Stomach contents in Asplanchna and Ploesoma. Archiv für Hydrobiologie- BeiheftErgebnisse der Limnologie, 8, 126-129.
12. Salt, G.W. (1977). An analysis of the diets of five sympatric species of Asplachna. Archiv für Hydrobiologie–BeiheftErgebnisse der Limnologie, 8, 123-125.
13. Alexander, M. and Hotchkiss, S. (2010). Bosmina remains in lake sediment as indicators of zooplankton community composition. Journal of Paleolimnology, 43, 51–59. http://doi.org/10.1007/s10933-009-9312-0
14. De Melo, R. D., and Hebert, P. D. (1994). A taxonomic reevaluation of North American Bosminidae. Canadian Journal of Zoology, 72(10), 1808-1825.
15. Villabona-González, S.L., Ramírez-Restrepo, J.J., Palacio-Baena, J.A. and Costa-Bonecker, C. (2015). Respuesta de la biomasa zooplanctónica a los gradientes de estado trófico y precipitación de un embalse tropical. Revista de la Academia Colombiana de Ciencias Exactas, Físicas y Naturales, 39(152): 374-388.
16. Kerner, M., Ertl, S. and Spitzy, A. (2004). Trophic diversity within the planktonic food web of the Elbe Estuary determined on isolated individual species by 13C analysis. Journal of Plankton Research, 26(9): 1039-1048.
17. Popp, A., Hoagland, K. and Hergenrader, G. (1996). Zooplankton community response to reservoir aging. Hydrobiologia, 339, 13-21.
18. Elmoor-Loureiro, L.M.A., Mendonça-Galvão, L. and Padovesi-Fonseca, C. (2004). New cladoceran records from lake Paranoá, central Brazil. Brazilian Journal of Biology, 64(3A): 415-422.
19. Villabona-González, S.L., Buitrago-Amariles, R.F., Ramírez-Restrepo, J.J. and Palacio-Baena, J.A. (2014). Biomasa de rotíferos de dos embalses con diferentes estados tróficos (Antioquia, Colombia) y su relación con algunas variables limnológicas. Actualidades Biológicas, 36(101): 149-162.
20. Empresas Públicas De Medellín. 1994. Proyecto de aprovechamiento múltiple de Riogrande. Declaración de impacto ambiental., Medellín, Colombia.
21. Zabala, A.M. (2013). Evaluación del estado trófico del embalse Riogrande II (Antioquia, Colombia) a través del comportamiento espacial y temporal de los nutrientes. Trabajo de grado, Universidad de Antioquia, Medellín, Colombia, 112 p.
22. Palacio-Baena, J., Herrera-Loaiza, N., López-Muñoz, M., Palacio-Betancourt, H. and Rodríguez, M. (2013). Limnoecología de los embalses Riogrande II, La Fe y Porce II. In: Palacio J., ed. Informe técnico embalse Riogrande II “Estudio de la problemática ambiental de tres embalses de Empresas Públicas de Medellín para la gestión integral y adecuada del recurso hídrico”. Universidad de Antioquia, Universidad Nacional sede Medellín, Empresas Públicas de Medellín, pp. 444-528.
23. Toro, M., Moreno, A., Chalarca, D. and Grajales, H. (2013). Dinámica ambiental de los nutrientes (especies químicas del nitrógeno, fósforo y sílice) en los embalses La Fe, Riogrande II y Porce II a través de modelos conceptuales y matemáticos, como herramienta de gestión de la calidad de agua en los tres sistemas. In: Palacio J. ed. Informe técnico embalse Riogrande II “Estudio de la problemática ambiental de tres embalses de Empresas Públicas de Medellín para la gestión integral y adecuada del recurso hídrico”. Universidad de Antioquia, Universidad Nacional sede Medellín, Empresas Públicas de Medellín, pp. 242-297.
24. López-Muñoz, M., Ramírez-Restrepo, J., Palacio-Baena, J., Echenique, R., Mattos-Bicudo, C. and Parra-García, E. (2016). Biomasa del fitoplancton eucariota y su disponibilidad para la red trófica del embalse Riogrande II (Antioquia, Colombia): Revista de la Academia Colombiana de Ciencias Exactas, Físicas y Naturales, 40 (155): 244-253.
25. Toledo, A.P., Talarico, M., Chinez, S.J. and Agudo, E.G. (1983). A aplicação de modelos simplificados para a avaliação de procesos de eutrofização em lagos e reservatórios tropicais. XIX Congreso interamericano de Engenharia e ambiental. Camboriú, Brasil.
26. Friendly, M. (2007). HE plots for multivariate general linear models. Journal of Computational and Graphical Statistics, 16, 421-444.
27. Guisande, C., Heine, J., González-Dacosta, J. and García-Roselló, E. (2014) Rwizard Software. University of Vigo. Vigo, Spain.
28. Melão, M. G. (1999). Desenvolvimento e aspectos reprodutivos de cladóceros e copépodos de águas continentais brasileiras. In: Pompêo, M., ed. Perspectivas da Limnologia no Brasil. São Luís: Gráfica e Editora União, pp. 45–57.
29. Gillooly, J. (2000). Effect of body size and temperature on generation time in zooplankton. Journal of Plankton Research, 22, 241–251.
30. Lagergren, R., Svensson, J.E. and Stenson, J.A. (2007). Models of ontogenetic allometry in cladoceran morphology studies. Hidrobiologia, 594, 109–116. http://doi.org/10.1007/s10750-007-9085-2
31. Hall, D.J., Threlkeld, S.T., Burns, C.W. and Crowley, P.H. (1976). The size-efficiency hypothesis and the size structure of zooplankton communities. Annual Review of Ecology and Systematics, 7, 177-208.
32. Corgosinho, P. and Pinto-Coelho, R. (2006). Zooplankton biomass, abundance and allometric patterns along an eutrophic gradient at Furnas Reservoir (Minas Gerais , Brazil): Acta Limnologica Brasiliensia, 182, 213–224.
33. Hart, R. and Bychek, E. (2011). Body size in freshwater planktonic crustaceans: an overview of extrinsic determinants and modifying influences of biotic interactions. Hidrobiologia, 668, 61–108. http://doi.org/10.1007/s10750-0100400-y
34. Conde-Porcuna, J.M., Morales-Baquero, R. and Cruz-Pizarro, L. (1993). Effectiveness of the caudal spine as a defense mechanism in Keratella cochlearis. Hydrobiologia, 255/256, 283-287
35. Kappes, H., Mechenich, C. and Sinsch, U. (2000). Long-term dynamics of Asplanchna priodonta in LakeWindsborn with comments on the diet. Hydrobiologia, 432, 91–100.
36. Wong, C.K. and Sprules, W.G. (1985). Size-selective feeding by the predatory copepod Epischura lacustris Forbes. Canadian Journal of Fisheries and Aquatic Sciences, 42, 189-193.
37. Branstrator, D.K. (1998). Predicting diet composition from body length in the zooplankton predator Leptodora kindtii. Limnology and Oceanography, 43, 530–535.
38. Sakamoto, M. and Hanazato, T. (2008). Antennule shape and body size of Bosmina: key factors determining its vulnerability to predacious Copepoda. Limnology, 9, 27-34. http://doi.org/10.1007/s10201-007-0231-3
39. Giraldo, M. (2010). Composición de la dieta del rotífero Asplanchna girodi (De Guerne, 1888) y su variación temporal y espacial en el embalse Riogrande II (Antioquia, Colombia): Trabajo de grado, Universidad de Antioquia, Medellín, Colombia, 80 p.
40. Bonecker, C., Azevedo, F. and Simões, N. (2012). Zooplankton body-size structure and biomass in tropical floodplain lakes: relationship with planktivorous fishes. Acta Limnologica Brasiliensia, 23(3): 217–228. http://doi.org/10.1590/S2179975X2012005000005
41. Gavilán-Díaz, R. (1990). Flutuações Nictemerais dos Fatores Ecológicos na Represa de Barra Bonita-Médio Tietê-SP. Dissertação de Mestre, Universidade Federal de São Carlos, Brasil, 157 p.
42. Parra-García, E. A., Rivera-Parra, N., Picazo, A., and Camacho, A. (2020): Detrital food chain as a possible mechanism to support the trophic structure of the planktonic community in the photic zone of a tropical reservoir. Limnetica, 39(1): 511-524.
Received: 14 November 2020
Accepted: 14 February 2021
*Yury Catalina López-Cardona1, Edison Parra-García12, Jaime Palacio-Baena3, Silvia Lucía Villabona-González4
1Grupo de Investigación en Limnología Básica y Experimental y Biología y Taxonomía Marina, Bióloga, Universidad de Antioquia, Medellín, Colombia. orcid.org/0000-0001-8030-7632
2Grupo de Investigación en Limnología Básica y Experimental y Biología y Taxonomía Marina, Candidato a Dr., Universidad de Antioquia, Medellín, Colombia. orcid.org/0000-0003-4957-7970
Grupo de Investigación en Gestión y Modelación Ambiental, Universidad de Antioquia, Medellín, Colombia.
4Grupo de Investigación en Limnología y Recursos Hídricos, Profesor Asociado departamento de Ciencias Exactas y Naturales, Universidad Católica de Oriente, Rionegro, Colombia. orcid.org/0000-0002-6136-4495
*Autor corresponsal: Universidad de Antioquia, Instituto de Biología, programa de Biología, Calle 67 No. 53-108, PBX: +(57)(4) 219 8332, Medellín, Antioquia, Colombia. Email: [email protected]