Inferring possible population divergence in Espeletia pycnophylla (Asteraceae) through morphometric and paleogeographic approaches.
The geographic distribution area of E. pycnophylla is actually very fragmented and there are many isolated populations on the Southern Andes of Colombia (Cuatrecasas 1986, Solarte et al. 2007), one in the northern Ecuador and a small population in the central Ecuador (Jorgensen & Leon-Yanes 1999). This isolation is a possible response of the fragmentation of the Paramo continuity during the last 14 000 years. This kind of phenomena induces such inter-population divergence in plant populations due to restricted gene flow, contrasting selecting pressures and derive, the first steps in the microevolutive diversification of plants (Loveless & Hamrick 1984).
To measure this type of divergence some approaches have arisen, such as morphology, citotaxonomy and molecular markers, but other kind of morphometric techniques such as the Elliptic Fourier Analysis of leaf shape has also prove to be a very effective phenotypic based alternative to indirectly estimate genetic variation among and within plant populations, allowing many applications in the study of population evolution (Lonn & Prentice 1995, Olsson et al. 2000, White et al. 1998). This technique accounts for the variation in the shape per se of biological structures (Rohlf & Archie 1984), eliminates the effects of size and recovers a high amount of hereditary variation (Molau & Prentice 1992, Jensen et al. 2002, Rumpunen & Bartish 2002).
Our goal was to test the hypothesis that three geographically isolated populations of E. pycnophylla shows phenotypic divergence correlated with their inter-population spatial distances, as a possible result of their geographic fragmentation, by combining: 1) the variation of leaf shape using the Elliptic Fourier Analysis; 2) a paleogeographic reconstruction of the distribution area of E. pycnophylla in the Southern Colombian Andes during the LateGlacial, using the model of van der Hammen (1974) and 3) an actual estimate of the distribution area of this species in the same locations. We chose to reconstruct the Late-Glacial distribution of E. pycnophylla because according to paleoecological estimations, it was the time when the Paramo vegetation reached the greatest distribution area and was much more continuous than today.
MATERIALS AND METHODS
The species: E. pycnophylla is a long lived giant rosette endemic from the Paramos of Southern Colombia and Northern and Central Ecuador. This species constitute one of the most conspicuous vegetation components of these high altitude ecosystems (Solarte et al. 2007). The main diagnostic characteristics are its sessile and highly pubescent leaves with lanceolate to elliptic shape (Cuatrecasas 1986). The heights of some adult individuals may reach 5m and the stalk diameters range from 6 to 8cm (Laegard 1992). The species are mainly reproduced by a highly auto-incompatible cross pollination (Gonzales et al. 2009) and its seeds lack from pappus, what implies a great restriction for the dispersion over large distances. Moreover, there are no reports of ornitochory, the seeds fall to the ground when the inflorescences are dry enough after the fructification stage, and the seedlings grow up near to the mother plant (Cuatrecasas 1986), producing a clumped pattern of distribution (Benavides et al 2008).
According to Cuatrecasas (1986), E. pycnophylla shows some geographic variation classified in three infraspecific levels: E. pycnophylla var. galerana, E. pycnophylla spp. angelensis and E. pycnophylla spp. llanganatensis (on the central cordillera of Ecuador) but no further evidence has confirmed these groups. Moreover, there are still many isolated populations that remain unstudied and some of them have also such internal patchiness. Here, we considered as a population a group of individuals with a delimited distribution area, spatially isolated from others and clearly identifiable via cartography, according to the classification of Paramo areas in Solarte et al. (2007).
Populations and sampling: We selected three geographically isolated populations of E. pycnophylla in the Southern Colombian Andes, each one being characterized by contrasting climate regimes and geomorphologic characteristics (Solarte et al. 2007). In addition to the spatial distances, the Guaitara river canyon (450m of maximum depth) (Fig. 4) act as a geographical barrier that strengthen the isolation by distance. Within each population we established four sampling areas in each of four orientations (N, S, W, E) and three elevations, which were proportional to the extension of the whole population. Within each sampling area we randomly selected 20 adult individuals (70-100cm tall) and from each individual we took two completely developed and healthy leaves with no damage like hervivory or necrosis (Fig 1). This sampling methodology led us a sample size of 160 leaves per population. Table 1 show some geographic features of each studied population and the definitive sample sizes used for the analyses.
Morphometric analyses: All the collected leaves were pressed and imaged on their abaxial surface with a Cannon[R] PowerShot[R] digital camera at 7.1 megapixels. The images were taken at black and white with standardized conditions (distance, position and angle) upon a white background to increase the image contrast for an optimal acquirement of outlines. We used the software TpsDig to digitalize the outlines and transform those to x,y coordinates. Based on these coordinates, we calculated the Kulh & Giardina's (1982) Elliptical Fourier Descriptors (EFD's) of each leaf outline with the program EFAWin (Isaev 1995). Ten harmonics were enough satisfactory to successfully describe the shape of all leaves. A normalization of the EFD's was applied to maintain invariant the size, location, rotation and starting point of each image, and thereby, recover only the shape per se of the structures (Rohlf & Archie 1984). Additionally, three linear morphometry (LM) variables (total leaf length, maximum leaf width and length-wide proportion) were measured to compare the effectiveness of leaf shape versus leaf size to discriminate between populations.
Statistical analyses: Within and between groups variance was accounted for understanding the morphometric structuring of populations. Both the normalized EFD's and the LM of leaves were used in canonical variates analysis (CVA) with the software SPSS 15.0, excepting the variable length-width proportion, which was analyzed with a Kruskall-Wallis test due to its departures from normality.
Since every harmonic is described by two sine and cosine components, yielding to four coefficients per harmonic (Kulh & Giardina 1982), 40 normalized EFD's were the variables used for shape analysis. Wilk's Lambda values ([LAMBDA]) (being 1 and 0 the minimum and maximum discriminant powers respectively) and the canonical correlation coefficients (being 1 and 0 the maximum and minimum discriminant powers respectively) were used as indicators of the discriminating power of each canonical variate (CV). Post hoc tests of Bonferroni with Hotelling probabilities were applied to inquire differences between population pairs when the CVA showed overall significant discriminations among groups.
[FIGURE 1 OMITTED]
Mantel tests (Mantel 1967) were preformed with the software Past 1.94b to prove if the dissimilarity matrices of geographic and morphometric distances among the three populations were linearly correlated. To build these matrices, the inter-population dissimilarities for both the normalized EFD's and the LM variables were calculated as euclidean distances between the centroids of each pair of populations using Past 1.94b. The chosen alpha for all the analyses was 0.05.
Geographical and paleogeographical analysis: The reconstruction of the distribution area of E. pycnophylla for the Late-Glacial of Pleistocene (14 000 years ago) was made through a cartographic analysis upon the van der Hammen's model of historic altitudinal shifts of Paramo ecosystems in the Colombian Andes (van der Hammen 1974). Based on paleoecological data from fossil pollen grains obtained from different sediment cores, the van der Hammen's model explains how the Paramo ecosystems in Colombia reached their maximum geographic extension in the Late-Glacial, covering an altitudinal range from 2 000 [+ or -] 100m to 3 800 [+ or -] 100m. After that, the increasing global temperatures that characterized the start of the Holocene produced the uplift of the Paramo to higher elevations. This has been happening as a gradual process until nowadays, when the Paramo sensu stricto covers an average altitudinal range from 3500 [+ or -] 100m to 4500 [+ or -] 100m.
We obtained the level curves of the Southern Colombian Andes at every 100m and applied the altimetry estimates of Paramo limits suggested by the van der Hammen's model for the Late-Glacial period (inferior limit: 2 000 [+ or -] 100m; superior limit: 3 800 [+ or -] 100m). We used the SRTM (Shuttle Radar Topography Mission) digital elevation model generated by NASA (National Aeronautics and Spaces Administration) and NGA (National Geospatial Intelligence Agency) at a resolution of 30m. The original image was initially treated via Principal Component Analysis in order to recover the most of the possible information and to facilitate its interpretation in terms of the original variability (Chuvieco 1996). For estimating the current distribution areas of E. pycnophylla we used vegetation coverage images of the Southern Colombian Andes, Landstat TM, October 14 2002 (path 9 row 59, 8 bands, BSQ format, II quadrat, with selection of training fields on screen), which allowed us to cartographically identify the Espeletia pycnophylla-Calamagrostis effusa and Espeletia pycnophylla-Blechnum loxense associations. Field work made by Solarte et al (2007) demonstrates that E. pycnophylla only occurs in these two associations, which are clearly identifiable via Landstat satellite imagery.
Finally, we analyzed the possible relations between the spatial characteristics of the distribution areas by obtaining the respective thematic maps trough a GIS in ArcGis 9.2 platform for both the past and present periods.
Linear Morphometry: Significant differences were found in leaf size among populations. The CVA generated two CV's (Table 2); the first explained the 99.03% of the original data variance and had a positive perfect correlation with the Total Leaf Length (r=1). The Wilk's Lambda ([LAMBDA]) and the canonical correlation values for this function showed highly significant discrimination among populations CCVP from LOP-GVP (p<0.001). The second CV correlated positively with the Maximum Leaf Width (r=0.98), but this variable didn't show significant differences among populations. Fig. 2 shows these results graphically. Likewise, the Length-Width Proportion of leaves were statistically significant to separate populations CCVP from LOP-GVP (p<0.001). These results revealed that according to leaf size, populations GVP and LOP are a relatively homogeneous group with long leaves and a high Length-Width Proportion, whilst CCVP is a contrasting population with short leaves and low Length-Width Proportion (Table 3).
[FIGURE 2 OMITTED]
Leaf shape morphometry with EFD's: Significant differences were found in leaf shape among the three populations. The CVA generated two significant CV's (Table 2), which together explained the 100% of the original EFD's variance. The first CV recovered the 76.2% of the variation and significantly differentiated the three populations (p<0.001). However, this function showed that the CCVP population differs to a greater extent from the rest. The second function recovered the 23.8% of the total variation and mainly differentiated LOP from GVP. Fig. 3 shows these results graphically.
Relation between geographical distances and morphometric differences: High and significant positive linear correlations were found between geographic and leaf size (LM) (r=0.90, p<0.05) and geographic leaf shape (EFD's) (r=0.99) morphometric matrices, according to the Mantel correlation test. This showed that the more geographical distance between any pair of populations, the more morphometric divergence between them. However, this correlation was stronger with the EFD's. This last proves that leaf shape variation among populations is more linearly sensible to the geographic distances than leaf size.
[FIGURE 3 OMITTED]
Geographic and paleogeographic analyses: The map of the paleogeographic reconstruction of the distribution area of E. pycnophylla 14.000 years ago (Fig. 4A) showed a large and uninterrupted central population with some small marginal populations at the south west and north of the distribution, covering a total estimated area of 8 351.4 [km.sup.2]. This ancient population was first split into two parts; the CCVP and a population formed by the conjunction of the current GVP and LOP areas. After that, the fragmentation of the last two populations occurred. The map of the current distribution (Fig 4B) shows a high degree of spatial disjunction, which covers a total estimated area of 669.90 [km.sup.2]. These results demonstrates that the 92% of the geographic distribution area of E. pycnophylla has not only fragmented, but conspicuously reduced since the Late-Glacial through all the Holocene, resulting in a set of geographically isolated populations.
Variation in leaf shape and size: We think that the results of this research support the hypothesis of an allopatric population divergence of E. pycnophylla in the Paramos of the Southern Colombian Andes. The EFD's of leaf shape had greater statistical effectiveness and larger resolution than LM of leaf size in distinguishing inter-population differences. This is in agreement with a general pattern shown by similar studies (Lonn & Prentice 1995, McLellan & Endler 1998, Jensen et al. 2002, Rumpunen & Bartish 2002). However, the results of each individual method deserve a particular explanation. The differences in leaf size could be related to plasticity or to different environmental pressures influencing each population. The significantly smaller leaf size in population CCvP may be a possible adaptation to the sandy soils, lower precipitation rates and lower water availability that govern its habitat (Solarte et al 2007). This is a typical response of plants to reduce the potential area for evapotranspiration and tissue dehydration.
[FIGURE 4 OMITTED]
On the other side, although the differences in leaf shape shown by the EFD's cannot be attributed to conspicuous biologic features or adaptations, making difficult its interpretation, Rohlf & Archie (1984) states that it indeed allows detecting differences in the overall shape per se of biological structures, which usually contain a hereditary control (Dickinson et al. 1987). In fact, some researchers have proved that whilst the LM accounts for a high percentage of environmental and developmental variation and a small percentage of genetic variation, the EFD's applied to leaf shapes do exactly the opposite, serving as very accurate estimators of genetic variation within and between plant populations. For example, Molau & Prentice (1992) and Lonn & Prentice (1995) found that for three Saxifraga species and Hippocrepis emerus respectively, the EFD's only accounted for a very small percentage of variation explained by environmental plasticity. However, the measurement of genetic variation trough molecular markers is still needed to prove if this applies to E. pycnophylla.
Identifying the mechanisms behind the inter-population differences in leaf shape is beyond the scope of this work because we still need to apply genetic measurements and other tools to solve this question. However, we suggest as a hypothesis that the restrictions in seed dispersion and pollination over long distances among isolated populations could be responsible factors.
Relation between geographic distances and morphometric divergence: The significant correlation found between morphometric and geographic distances matrices are also supported by the geographic and paleogeographic analyses. In case of having an important genetic control of leaf shape, this result could indicate that the processes involved with gene flow among the populations of E. pycnophylla are proportionally related with their geographical distances. These kinds of tendencies have been detected in other plant species with truncated gene flow due to restrictions in pollination (Lonn & Prentice 1995, Brunell & Whitkus 1993). In the case of E. pycnophylla, this result may have historic and present implications, which respectively are: 1) Closer populations were joint in a more immediate past than distant ones did, directly or indirectly by intermediate populations; 2) Closer populations as GPV and LOP, which have a minimal linear distance of 12km, have a higher probability of undergoing inbreeding contact than distant ones. This is because the spatial distances between them may not represent striking restrictions for pollination vectors. In a combined perspective of both hypotheses, current nearer populations may undergo gene flow although they were split in the past.
?Evidence of a vicariance event?: Nowadays, an increasing number of palaeoecological evidences shows that the periodic phenomenon of climate change during the Pleistocene (glacial-interglacial cycles) and their repercussion on the dynamics of the spatial distribution of Paramos in the northern Andes were a strong force that influenced many episodes of extinction, speciation, range contractions and expansions for many groups of plants, implying continuous events of vicariance for the biota (Simpson 1971, Van der Hammen 1974, Simpson & Todzia 1990). Applying these ideas to the case of E. pycnophylla and based on our geographic and paleogeographic modeling, we suggest that the three studied populations are the result of a vicariance scenario.
Although our data of leaf shape variation, geographic and paleogeographic analysis of the distribution area of E. pycnophylla supports the mentioned explanations, a greater survey is still needed to test if this pattern is consistent when adding variation from marginal populations and from different locations within each population. Due to the difficult access and mobilization in the Paramos, we only sampled one localized area in each of the three populations; this may have induced such pseudoreplication and makes necessary the addition of new sampling areas in contrasting microhabitats to average unwanted factors. But on the whole, the findings of this study offer a good opportunity to keep searching for biogeographic processes in response to climate change and geographical dynamics in plant species of Paramo ecosystems.
We are grateful to Belisario Cepeda and Maria Elena Solarte for contributing the seminal ideas of the project, to Jorge Bedoya and the Reserva Natural Pueblo Viejo for their authorizations and logistic support during field work; to German Narvaez and the TERRA group who supplied clue information for the geographical analyses; to Jorge Benavides for his useful revision and to IDEA WILD and the VIPRI-UDENAR which provided the financial support.
Benavides, I.F., D.L. Burbano., S.M. Urbano. & M.E. Solarte. 2007. Efecto del gradiente altitudinal sobre aspectos autoecologicos de Espeletia pycnophylla ssp. angelensis Cuatrec. (Asteraceae) en el Paramo El Infiernillo (Narino), Colombia. Actual. Biol. 29: 41-53.
Brunell, M.S. & R. Whitkus. 1993. Patterns of morphological variation in Eriastrum densifolium (Polemoniaceae). Supplement to Am. J. Bot. 80:134.
Carlquist, S. 1974. Island Biology. Columbia University, New York, USA.
Chuvieco, E. 1996. Fundamentos de Teledeteccion Espacial. RIALP, Madrid, Espana.
Cuatrecasas, J. 1986. Speciation and radiation of the Espeletiinae in the Andes, p 276-303. In F. Vuilleumier & M. Monasterio. High altitude tropical biogeography. Oxford University, New York, USA.
Dickinson, T.A., W.H Parker & R.E. Strauss. 1987. Another approach to leaf shape comparisons. Taxon 36: 1-20.
Gonzales, M.S., S.M Urbano & M. Pianda. 2009. Patrones de biologia reproductiva de Espeletia pycnophylla, Puya clava-herculis y Puya cryptantha en dos areas paramunas con diferencias ambientales y de intervencion antropica en el departamento de Narino. Pasto, Narino, Colombia. Memorias del II Congreso Mundial de Paramos, Paramundi. Loja, Ecuador.
Isaev, M.A. & L.N. Denisova. 1995. The computer programs for shape analysis of plant leaves, proceedings of the Mathematics. Computer. Education. International conference, Pushchino, Russia.
Jensen, R.J., K.M Ciofani & L.C. Miramontes. 2002. Lines, outlines, and landmarks: morphometric analyzes of leaves of Acer rubrum, Acer saccharinum (Aceraceae) and their hybrid. Taxon 51: 475-492.
Jorgensen, P.M. & S. Leon-Yanes. 1999. Catalogue of the Vascular Plants of Ecuador. Missouri Botanical Garden, St. Louis, Missouri, USA.
Kuhl, F.P. & C.R. Giardina. 1982. Elliptic Fourier features of a closed contour. Comput. Graph. Image Process. 18: 236-258.
Laegaard, S. 1992. Influence of fire on the grass Paramo vegetation of Ecuador. In H. Balslev & J. Luteyn. Paramo: An Andean Ecosystem under Human Influence. Academic, London, Enlgand.
Lonn, M. & H. Prentice. 1995. The structure of allozyme and leaf shape variation in isolated, range-margin populations of the shrub Hippocrepis emerus (Leguminosae). Ecography 18: 276-285.
Loveless, M.D. & J.L. Hamrick. 1984. Ecological determinants of genetic structure in plant populations. Ann. Rev. Ecol. Syst. 15: 65-95.
Mantel, N. 1967. The detection of disease clustering and a generalized regression approach. Cancer Res. 27: 209-220.
McLellan, T. & J.A. Endler. 1998. The relative success of some methods for measuring and describing the shape of complex objects. Syst. Biol. 47: 264-281.
Molau, U. & H.C. Prentice. 1992. Reproductive system and population structure in three arctic Saxifraga Species. J. Ecol. 80: 149-161.
Monasterio, M. & L. Sarmiento. 1991. Adaptive radiation of Espeletia in the cold Andean tropics. Trends Ecol. Evol. 6: 387-391.
Olsson, A., H. Nybom & H.C. Prentice. 2002. Relationships between Nordic Dogroses (Rosa L. sect. Caninae, Rosaceae) Assessed by RAPDs and Elliptic Fourier Analysis of Leaflet Shape. System. Bot. 25: 511-251.
Rauscher, J.T. 2000. Molecular systematics of the Espeletia complex: evidence from nrITS sequence on the evolution of an Andean adaptative radiation. Ph.D. Thesis, Washington University, St Louis, Missouri, USA.
Rohlf, F.J. & J.W. Archie. 1984. A comparison of Fourier methods for the description of wing shape in mosquitoes (Diptera: Culicidae). Syst. Zool. 33: 302-317.
Rumpunen, K. & I.V. Bartish. 2002. Comparison of differetntiation estimates based on morphometric and molecular data, exemplified by various leaf shape descriptors and RAPD's in the genus Chaenomeles (Rosaceae). Taxon 51: 69-82.
Simpson, B. 1971. Pleistocene changes in the fauna and flora of South America. Science 173: 771-780.
Simpson, B. & C. Todzia. 1990. Patterns and processes in the development of the high Andean flora. Am. J. Bot. 77: 1419-1432.
Solarte, M.A., G. Narvaez, G. Rivas, A.E. Bacca, D. Munoz, J.J. Calderon, C. Torres, V. Figueroa, J. Rengifo, P Martinez, M.T. Davila, B. Cepeda & G. Castillo 2007. Proyecto Estado del arte de la Informacion Biofisica y Socioeconomica de los Paramos de Narino. Tomo I. Grupo de investigacion en Biologia de Paramos y Ecosistemas Andinos, Universidad de Narino -- Corporacion autonoma regional de Narino CORPONARINO, San Juan de Pasto, Narino, Colombia.
van der Hammen, T. & A.M. Cleef. 1986. Development of the high Andean Paramo flora and vegetation, p. 153-201. In F. Vuilleumier & M. Monasterio. Hight altitude tropical biogeography. Oxford University, New York, USA.
van der Hammen, T. 1974. The Pleistocene changes of vegetation and climate in tropical South America. J. Biogeogr. 1: 3-26.
Van der Hammen, T. 1985. The Plio-Pleistocene climatic record of the tropical Andes. J. Geol. Soc. Lond. 142: 483-489.
White, R.J., Prentice, H.C. & Venvijst, T. 1988. Automated image acquisition and morphometric description. Can. J. Bot. 14: 612-623.
Felipe Benavides (1,2), Jorge Burbano (2), Diana Burbano (2), Rodrigo Prieto (3) & Carlos Torres (2,4)
(1.) Grupo de investigacion Geografia Fisica y Problemas Ambientales TERRA Departamento de Geografia, Universidad de Narino, Pasto, Colombia; email@example.com
(2.) Grupo de investigacion Biologia de Paramos y Ecosistemas Andinos, Departamento de Biologia, Universidad de Narino, Pasto, Colombia; firstname.lastname@example.org, email@example.com
(3.) Grupo de investigacion Citogenetica Filogenia y Evolucion de Poblaciones. Universidad del Tolima, Ibague, Colombia; firstname.lastname@example.org
(4.) Grupo de investigacion Cuencas Hidrograficas y Ecosistemas Tropicales - Departamento de Geografia, Universidad de Narino, Pasto, Colombia; email@example.com
Received 16-IX-2009. Corrected 25-VI-2010. Accepted 22-VII-2010.
TABLE 1 Geographic characteristics of the studied populations Population Geographic coordinates CCVP 1[degrees] 02' 32.57" N, 77[degrees] 46'52.8" W GVP 1[degrees] 10' 56.09" N, 77[degrees] 23' 14.4" W LOP 0[degrees] 51" 28.85" N, 77[degrees] 24' 58.18" W Geographic distances between populations (km) Population Altitude (m.a.s.l) n CCVP GVP LOP CCVP 3 322 160 0 GVP 3 500 156 48.21 0 LOP 3 310 160 45.7 35 0 n = sample size. Four leaves were discarded during the statistical analysis because they were outliers outside the 95% confidence limits. CCVP: Chiles Cumbal volcanoes Population GVP: Galeras volcano Population. The respective research authorization for this population was obtained according to the permission DEC309/2000 issued by Parques Nacionales Naturales de Colombia. LOP: Las Ovejas Population The geographic distances were measured on a map as straight lines between the centers of the sampling areas of each pair of populations. TABLE 2 Discrimination statistics for the comparison of the leaf size and leaf shape among the three populations of Espeletia pycnophylla using Linear Morphometry Leaf size (LM) Percentage Wilks Lambda Canonical Function Eigenvalue of variance ([LAMBDA]) 1 1.71 99.03% 0.36 2 0.02 0.97% 0.98 Canonical Function Canonical correlation p 1 0.79 <0.001 2 0.14 non significant Leaf shape (EFD's) Percentage Wilks Lambda Canonical Function Eigenvalue of variance ([LAMBDA]) 1 5.25 76.2% 0.06 2 1.63 23.8% 0.37 Canonical Function Canonical correlation p 1 0.917 <0.001 2 0.788 <0.001 TABLE 3 Average [+ or -] Standard Error (95%) for the linear morphometry variables Population TLL (cm) MLW (cm) GVP 37.98 [+ or -] 4.20 9.12 [+ or -] 1.26 LOP 36.81 [+ or -] 6.25 8.60 [+ or -] 1.17 CCVP 21.61 [+ or -] 2.70 8.48 [+ or -] 1.13 Population LWP GVP 4.22 [+ or -] 0.64 LOP 4.35 [+ or -] 0.94 CCVP 2.56 [+ or -] 0.27 TLL = Total Leaf Length. MLW = Maximum Leaf Width. LWP = Length-Width Proportion.
|Printer friendly Cite/link Email Feedback|
|Author:||Benavides, Felipe; Burbano, Jorge; Burbano, Diana; Prieto, Rodrigo; Torres, Carlos|
|Publication:||Revista de Biologia Tropical|
|Date:||Dec 1, 2010|
|Previous Article:||Caracterizacion espacio-temporal del habitat y presencia de Dermatemys mawii (testudines: dermatemydidae) en la cuenca del Grijalva-Usumacinta,...|
|Next Article:||Young restored forests increase seedling recruitment in abandoned pastures in the Southern Atlantic rainforest.|