ARTICULO ORIGINAL
Use of epidermal characters as bioindicators of environmental pollution
Uso de caracteres epidérmicos como bioindicadores en contaminación ambiental
Mirta Olga Arriaga, Mónica L. Stampacchio, María G. Fernández Pepi, Patricia Eleonora Perelman y Ana María Faggi
Museo Argentino de Ciencias Naturales (Argentine Museum of Natural Sciences), Av. A. Gallardo 470, C1405DJR, Buenos Aires, Argentina. Fax: (54) (011) 4982-6595, ext. 144-145 <marriaga@macn.gov.ar > <mgfernandezpepi@ macn.gov.ar>
Abstract
The present study explores the use of micromorphological leaf epidermal characters of Ficus benjamina (evergreen) and Fraxinus pennsylvanica (deciduous) as indicators of pollution in the Metropolitan Area of Buenos Aires. Mature leaves of both species were sampled along an urbanperiurban gradient, with a rural area as control for epidermal characters. Parameters examined from abaxial and upper epidermis were: (1) length and width of stomata, (2) number of stomata per area unit, (3) length and width of epidermal cells, (4) number of epidermal cells per area unit. Data were analysed with basic statistical and multivariate techniques. Principal Component Analysis (PCA) clearly grouped foliar samples according to three zones: urban, periurban and rural areas with different pollution values. No correlations were observed between epidermal character variations and climate variables, air pollutants or C, N, P and metals concentration in soil. Only one positive correlation (p = 0.006) was found between stomata length in leaves from Ficus and iron concentration in soil.
Key words: Bioindicators; Urban-periurban gradient; Seasonal variation; Areal variations
Resumen
En el presente trabajo se explora el uso de los caracteres micromorfológicos de epidermis de hoja de Ficus benjamina (especie perenne) y Fraxinus pennsylvanica (especie caduca) como indicadores de contaminación en el área metropolitana de Buenos Aires. A lo largo de un gradiente urbano-periurbano, con un área rural como control, se colectaron hojas de ambas especies. Se examinaron parámetros de epidermis abaxial y adaxial: largo y ancho de estomas, número de estomas por unidad de área, largo y ancho de células epidérmicas, número de células epidérmicas por unidad de área. Los datos fueron analizados con técnicas básicas de estadística y análisis multivariado. El análisis de componentes principales (PCA) claramente agrupó las muestras foliares en tres zonas: urbana, periurbana y rural, con diferentes valores de contaminación. No se observaron correlaciones entre las variaciones de los caracteres epidérmicos y las variables clima, contaminación ambiental o C, N, P y concentración de metales en el suelo. Se encontró una única correlación positiva (p = 0.006) entre el largo de estomas en hojas de Ficus sp. y la concentración de Fe en el suelo.
Palabras clave: Bioindicador; Gradiente urbano-periurbano; Variación estacional; Variación areal
Introduction
Environmental pollution is increasingly higher due to human activities and results in the accumulation of heavy metals that cause various types of problems to living beings (Rai et al., 2005). Man’s daily activities cause pollution of different kinds and degrees; indeed, air, soil and water get contaminated as a result of industrial and agricultural activities, transport, etc.
In big cities, air pollution levels are ground of concern. Monitoring programs and search for metals in environmental samples have become widely established as direct and indirect measurement methods. Biological monitoring constitutes an economic alternative for atmospheric pollution studies (Bellis et al., 2003; Castro & Faggi 2008; De Temmerman et al., 2004; Dion et al., 1993; Jasan et al., 2004; Panichev& McCrindle, 2004; Pignata et al., 2002; Walkenhorst et al., 1993). Plants have been observed to be far more sensitive to pollution than animals and man and are therefore used as indicators (Rani et al., 2006). Studies have recorded changes in plants due to all kinds of environmental pollutants, and most of these works refer to physiological alterations (Kabata- Pendias & Pendias, 2001; Linster, 1991; Ochiai, 1987; Patra & Sharma, 2000). Among passive monitor’s assays, the multielement analysis of tree bark from urban trees (Castro et al., 2008; Faggi et al., 2008; Jasan et al., 2004; Perelman et al., 2006; Steubing et al., 2002) has been recently introduced to detect and measure the elements of air pollution retained in tree bark.
The degree of impact in plants depends on pollutants concentration, location of entry into plant and species under consideration. Different species may present varying sensitivity/tolerance levels to different contaminating agents (Srivastava, 1999). Alderdice (1967) set two categories of toxic effects: acute toxicity (due to high pollution levels during short periods of time, usually lethal) and chronic toxicity (due to low pollution levels over long periods of time, lethal or sublethal).
Toxic effects produced by pollutants on soil organisms and plants have usually been studied in laboratories under controlled conditions, but field studies have been scarce (Patra & Sharma, 2000). Stomatal and epidermal cell size, lower frequency, thickening of cell wall, epicuticular wax deposition alterations and chlorosis are among the structural modifications in leaves deriving from pollution (Mhatre & Chaphekar, 1985; Mohapatra et al., 1991; Rao & Dubey 1991; Setia et al., 1994; Srivastava, 1999). In Brazil, Sant’Anna-Santos et al., (2008) worked with Genipa americana L. (Rubiaceae) in order to characterize injuries on leaf structure and micromorphology, and to evaluate its degree of susceptibility to simulated acid rain. They found necrotic interveinal spots on leaf blade, with plasmolized guard cells and cuticle rupture. Their results highlighted the relevance of anatomical data for precocious diagnosis of injury and to determine species sensitivity to acid rain. Working with Celosia cristata samples subjected to gaseous and particulated pollution deriving from heavy traffic in Indian routes, Srivastava (1999) observed that general plant growth was affected with severe distortions in foliar epidermal characters. Pointing out the importance of cuticle and epidermal features in the determination of tolerance/sensitivity of each species to environmental pollutants.
Taking into account these characteristics, some authors consider foliar epidemis as a bioindicator of environmental quality (Masuch et al., 1992; Alves et al., 2008; Balasooriya et al., 2009). A multielement analysis was carried out in the metropolitan area of Buenos Aires, using Fraxinus pennsylvanica tree bark as a bioindicator of air pollution (Castro et al, 2008; Perelman et al., 2006). Results confirmed the existence of a distribution pattern along an urban-periurban gradient (central, residential, periurban). Simultaneously, basic properties and heavy metal levels were analysed in soil from the basis of the same trees (López et al., 2006). Results were then compared with those from studies performed in the city of Mendoza (Faggi et al., 2008). Later, Fujiwara et al. (2006) correlated the concentration of elements in particulate matter (PM) from street dust with the concentration of pollutants retained in tree bark and with the concentration of pollutants in air in the same sites. They concluded that tree barks are more sensitive as passive bioindicators than street dust.
Our objective was to evaluate the possibility of using epidermis as a pollution indicator in trees from Buenos Aires metropolitan urban and periurban areas, given the ease and low cost of the technique. To recognise areas of different pollution levels and the effect produced by the same level of pollutants over a short time period (deciduous species) and over a long time period (evergreen species).
Materials and Methods
Selection of species
A deciduous and an evergreen tree were selected, namely Fraxinus pennsylvanica Marshall, and Ficus benjamina L, respectively. These species were selected given their high frequency along the streets of the metropolitan area of Buenos Aires, and because of the different life span of their leaves. If morphological or epidermal modifications were observed, their life span could allow us to determine whether they responded to a stimulus over a short period of time (stationary) or to the cumulative action of pollutants over long periods of time. Perelman et al. (2006) conducted sampling for tree bark at the same time as that reported for soil by López (2006).
Sample sites
Leaf samples were collected from urban and periurban areas characterized by different contamination levels. Three different sites belong to the urban area: (Ur 1) Constitución, affected by heavy car traffic and a railway terminal, (Ur 2) Flores, a residential neighbourhood affected by car traffic, and (Ur3) downtown, with dense traffic incidence. A public park affected by the international airport and its access highway (Figure 1) conform the periurban area Ezeiza (Per). Chascomús (Buenos Aires Province), a rural environment (Rur), was selected as the “noncontaminated zone” where leaves from the same species under study were used as “blank samples”, to provide a baseline to test foliar anatomic characters.
Figure 1. Location of sites sampled. References: (Ur) urban area: (Ur1) Constitución, (Ur2) Flores, (Ur3) downtown; (Per) periurban area; (Rur) rural area.
Figura 1. Ubicación de las zonas de muestreo. Referencias: (Ur) área urbana: (Ur1) Constitución, (Ur2) Flores, (Ur3) centro; (Per) área periurbana; (Rur) área rural.
Collection of Samples
Ten to twenty leaves were taken from the external part of the canopy, nearly 2 m high. Samples were taken from 10 trees per site every season, at the same time as sampling for tree bark (Perelman et al., 2006) and soil for pollution studies (López et al., 2006) which documented air pollution and field conditions.
Specimens analysed
A total of 18 Ficus specimens (8 from spring and 10 from autumn), and of 16 Fraxinus (16 from Spring and 13 from autumn) were analysed. Samples were kept in open paper bags.
Preparation of leaves
In Fraxinus, a compound odd pinnate leaf species, the bigger foliole from the apical pair was selected. Epidermal peels from this foliole and from Ficus leaves were examined under optical microscopy (MO) and scanning electronic microscopy (SEM). Traditional preparation techniques were utilized (Arriaga, 1986; D’Ambroggio de Argüeso, 1986).
Anatomical data registration
The macroscopic characteristics analysed were leaf colour, stains on leaves, leaf/foliole contour, apex, nervation, chlorosis and necrosis. The microscopic characteristics analysed in upper and abaxial epidermis were: epidermal cells length and width (CLUE, CWUE, CLAE, CWAE); stomata cells length and width (SLUE, SWUE, SLAE, SWAE); stomata number per area unit (considering units of 19.110 μm2) (SNUE, SNAE); epidermic cell number per transect (CNUE, CNAE, considering number of complete cells placed on 10 transects 155 μm long, randomised take on each slide).
Climate variables considered and concentration of elements in soil
In order to evaluate the eventual incidence of climate conditions on the variations of the epidermis characteristics considered, information was requested from the Servicio Meteorológico Nacional (National Meteorology Service) on temperature, effective heliophany and precipitation corresponding to the time when samples were taken. The concentration of C, N, P and metals present as soil cations (López et al., 2006), and metals retained in barks (Perelman et al., 2006) deriving from air pollution in the same locations were contrasted to evaluate their influence on the anatomical variations registered.
Data analysis
Data from observations were recorded on a matrix, which was then statistically analysed in order to find interspecific, zonal and/or seasonal differences and similarities. Mean and standard deviation values were recorded on the matrix. PCA (principal component analysis) was used as ordination method to describe the morphological variability among samples. Data were transformed by the Pearson moment-product correlation coefficient on a similarity matrix used to obtain the principal components (Sneath & Sokal, 1973; James & McCulloch, 1990). Invariant characters were previously removed from the analysis. Each individual was considered an operational taxonomic unit (OTU).
This ordination method was used to identify the main groups according to locality of collection. PCA was repeated several times to detect similar groups in each locality. Hence, once a group of OTUs related to a given locality was recognised, it was then separated and reanalysed alone. This procedure was repeated until well-differentiated groups of OTUs were determined. To describe the variability among groups, the characters that contributed most to the variability of the first five PCA components (r > 0.45) were analysed, their axes being selected according to group representation.
Discriminant analysis (DA) (Snealth& Sokal, 1973; Affifi & Clark, 1984) was applied in order to verify the consistency of groups resulting from PCA, on the basis of correlated characters. The empirical method (Affifi & Clark, 1984) was used to estimate goodness of fit of the classification. Based on Mahalanobis distances of each specimen from the different centroids of each group, the respective subsequent classification probabilities were computed. Thus, the proportion of individuals correctly classified for each group was calculated. Numerical analyses were carried out using NTSYS-pc program, version 2.02h (Rohlf, 1992).
Mean and standard deviations of quantitative characters, and the mode of discrete variables were calculated to describe differences among the most similar groups of species. One-way analysis of variance (for continuous characters; ANOVA) and non-parametric Kruskal- Wallis one-way analysis (for discrete characters) were performed in order to independently analyse morphological differences among groups. The goodness of fit of a normal distribution was assessed with the Kolmogorov-Smirnov one-sample test for all continuous characters. Variance homogeneity amongst groups was verified by Bartlett’s test.
Correlations were also carried out according to anatomical characters which presented variability in each species between either localities or times of collection and environmental variables (weather conditions and amount of chemicals). Results were compared with those obtained from measurements of soil contamination, metal retention in cortex and PM, from the same locations (Lopez et al., 2006; Perelman et al., 2006). Correlations were performed considering the characteristics and concentrations of cations, both independently as well as each character contrasted against all metals together, should there be any synergistic effect. STATISTICA 99 Edition (Stat Soft, Inc., 1995) was used for these analyses.
Results
Principal Component Analysis (PCA) clearly grouped foliar samples according to three zones: urban (Ur), periurban (Per) and rural (Rur), for each specie in each period (Figure 2).
Figure 2. Epidermal characters analised with multivariate techniques (PCA) References: Ur (urban area), Per (periurban area), Rur (rural area).
Figura 2. Gráfico de de caracteres epidérmicos analizados con técnicas de analisis multivariados (PCA). Referencias: Ur (area Urbana), Per (area periurbana), Rur (area rural)
Discriminant Analysis (DA) corroborated PCA’s classification. DA showed a clear discriminate pattern taking into account that 100% of specimens from each site were correctly classified into the groups established.
Fraxinus
In autumn, Fraxinus samples presented significant differences between rural and periurban species regarding length of abaxial epidermal stomata, while urban stomata were the longest. Samples from the urban area presented meaningful differences in stomata width of the abaxial epidermis with respect to samples from the other areas. Indeed, rural stomata were narrower than urban stomata, but wider than periurban stomata. Upper epidermal cells from the rural area were narrower than those from urban and periurban areas. Abaxial epidermal cells from the periurban zone showed meaningful differences in width with respect to the other two zones. They were narrower than those from the urban and rural areas (Table 1).
Table 1. Fraxinus characters (significant differences between areas * p < 0.05, ** p < 0.001 and between periods in the same area (# p < 0.05 in all cases)
Tabla 1. Caracteres de Fraxinus (diferencias significativas entre áreas * p < 0.05, ** p < 0.001 y entre períodos en la misma área (# p< 0.05 en todos los casos)
Fraxinus spring samples showed meaningful differences in length and width of abaxial epidermal stomata, in abaxial epidermis cell width, and in upper epidermis cell length. Rural and also urban and periurban areas presented meaningful differences regarding length and width of abaxial epidermal stomata. Rural stomata from upper epidermal cells were the shortest (Table 1). A comparison of Fraxinus characters between seasons per area revealed meaningful differences in urban and periurban samples, but not in rural samples.
In the latter, stomata in abaxial epidermis were shorter in autumn than in spring samples, and epidermal cell number was greater in autumn than in spring in both epidermis (Table 1). In periurban samples, meaningful differences correspond to stomatal number, length and width and cell number, all from the abaxial epidermis. In autumn, stomata and cell number from this epidermis were greater, while in spring, stomata in this epidermis were longer and wider (Table 1).
Ficus
In autumn leaves we found meaningful differences between samples from rural and urban areas regarding stomata width and cell length from abaxial epidermis. Both characters in samples from urban sites were bigger. We also found meaningful differences between samples from rural and periurban sites regarding the number of abaxial epidermal cells. Periurban specimens showed greater number of cells (Table 2).
Table 2. Ficus characters (* significant differences between areas and # between periods in the same area, p < 0.05 in all cases)
Tabla 2. Caracteres de Ficus (*diferencias significativas entre areas y # entre periodos en la misma área, p < 0.05 en todos los casos)
In spring leaves we only found meaningful differences regarding the number of stomata cells from abaxial epidermis in Ficus specimens from urban and periurban areas. These cells were more numerous in periurban areas. Regarding cell length from the abaxial epidermis, we found meaningful differences between samples from rural and periurban areas. Cells from periurban areas were longer (Table 2).
When we compared Ficus samples from spring and autumn, only specimens from the periurban area presented meaningful differences just in length and width of upper epidermal cells: autumn samples were shorter and wider than those from Spring (Table 2).
Correlation with climate variables
No correlations were found in any of the two seasons, between variations of leaf characters measured in the two species studied, and the climate variations considered (mean temperature, mean precipitation and heliophany).
Correlation with elements from soil or retained in tree bark
We contrasted data of character variations measured in this study with data on C, N, P and metals concentrations in soil (López et al., 2006) or elements from air retained in Fraxinus bark (from Perelman et al., 2006). The results revealed only one positive correlation (p = 0.006) between stomata length in Ficus leaves and soil Fe concentration.
Discussion
In order to evaluate morphological alterations produced by different pollutants, Sharma & Tyree (1973) compared leaf size from several populations of Liquidambar styraciflua L. (Hamamelidaceae) and verified diminished foliar length in specimens living next to a primary source of pollution, especially of particulate matter. Mandre et al. (1995) worked in the vicinity of Kunda (Estonia) with Pinus leaves from an area near a source of alkaline dust (cement industry), and observed diminished culm in leaves and in total tree length as symptoms of pollution.
However, we observed no alterations in shape, size or colour, nor the presence of chlorosis or necrosis in leaves from any of the areas sampled.
Alves et al. (2008) performed studies in foliar anatomy of Eugenia (Myrtaceae) in Sao Paulo, a large Brazilian metropolis. They compared individuals from heavily contaminated areas (urban zone: Dos Bandeirantes Ave.) and individuals from less contaminated areas (rural zone: Jarinu) and found that epidermis in leaves from urban zones showed statistically meaningful differences in stomatic density (higher density) than that in leaves from rural areas. Higher stomata density was also recorded by other authors in plants from urban zones subjected to air pollutants (Masuch et al., 1992; Päänkkönen et al., 1995 y 1997).
Working on Ligustrum lucidum Aiton f. (Oleaceae) in Córdoba city (Argentina), Bruno et al. (2007) also found higher stomata frequency in plants due to high pollution, presumably as a compensation for the negative effect of contaminants on photosynthesis.
Balasooriya et al. (2009), in the surroundings of Gent (Belgium) used Taraxacun officinalis (Asteraceae) in order to assess the potential of anatomical leaf characteristics (specifically leaf area, stomata density and cell size, among others) as effective parameters to biomonitor urban habitat quality. They monitored 4 different types of lands: U (urban), UB (urban green and forests), H (harbour and industrial areas) and P (pasture and grass lands), and found higher stomata density and bigger stomata in areas with better air quality.
Rani et al (2006) studied the effect of railway engine emission on the epidermal characteristics of some field plants growing along the railway track of Hapur Junction (India) affected by the emission of gases from diesel engines. Testing Calotropis procera, Croton bonplandianum and Cannabis sativa, they found that all three followed the same change pattern in leaf morphology due to environmental stress. Stomata and epidermal cells frequency declined in plants growing in close vicinity of railway tracks. This decline was greater regarding stomata frequency than regarding epidermal cells, and resulted in a meaningful fall in stomata index, accompanied by increased size of epidermal cells. Salgare and Thorat (1990) presented similar results on some trees from Andheri (West) Bombay.
Perhaps lower stomata frequency and smaller size allowed less gaseous exchange between environment and plant, thus protecting the latter from toxic inhalation. Since differentiation of stomata mother cells require the division of epidermal cells, decrease in stomata frequency should normally be accompanied by increased size of epidermal cells (Sant’Anna-Santos et al., 2008).
We also observed differences in length and frequency of stomata and epidermal cells between samples from polluted and non-polluted areas. Although in our study those differences were statistically not meaningful, they would indicate a clear tendency for highly polluted areas to produce in spring and autumn longer and wider stomata in both epidermises of the two species studied. Stomata number in spring and autumn Ficus samples was lower in highly polluted areas than in less polluted ones. In Fraxinus, stomata number resulted higher in most polluted areas in spring and autumn.
After contrasting differences, no correlation was found in our work with pollutants retained in tree barks or soil from the same sites, except for Ficus (evergreen), where we found correlation between high concentration of iron in soil and stomata length.
Differences between our results and those from other authors (regarding the correlation between foliar morphology and pollutants) may be due to different pollution levels, not as high in our city as in the other big cities studied; different sensitivity level of our selected species, or else the period of time in which the species were subjected to pollution. In the case of Fraxinus (deciduous species), the period of time during which leaves suffered pollution might not have been long enough to express alterations in epidermal characters.
Nevertheless, the differences found through PCA analysis clearly determined specific groups belonging to areas with different pollution values within the metropolitan area of the city of Buenos Aires
Conclusions
No differences in macroscopical characters but some in microscopical ones were detected between seasons or among areas in the species studied.
No correlation was observed between variations of the epidermal characters and the climate variables (mean temperature, mean precipitation and heliophany) or air pollutant values considered. Regarding the relation between variations in epidermal characters and soil
concentration of heavy metals, a positive correlation was found (p = 0.006) between stomata length of leaves in Ficus and iron concentration in soil. Differences in PCA analysis clearly determine groups belonging to areas with different pollution values within the metropolitan area of Buenos Aires.
Aknowledgements
We are indebted to the Servicio Meteorológico Nacional (National Meteorology Service) for the meteorological data in this study.
Bibliography
1. AFFIFI, A. A. & V. CLARK, 1984. Computer- aided multivariate analysis. Van Nostrand Reinhold. New York.
2. ALDERDICE, D. F., 1967. The detection and measurement of water pollution: biological assays. Canada Fishing Reports 9: 33- 39.
3. ALVES, E. S., F. TRESMONDI & E. L. LONGUI, 2008. Análise estrutural de folhas de Eugenia uniflora L. (Myrtaceae) coletadas em ambientes rural e urbano, SP, Brasil. Acta Botanica Brasilica 22(1): 241-248.
4. ARRIAGA, M. O., 1986. Metodología adaptada al estudio de hábitos alimentarios en insectos herbívoros. Comunicaciones del Museo Argentino de Ciencias Naturales II(15): 103-111.
5. BALASOORIYA, B.L.W.K., R. SAMSON, F. MBIKWA, U. W. A. VITHARANA, P. BOECKX & M. VAN MEIRVENNE, 2009. Biomonitoring of urban habitat quality by anatomical and chemical leaf characteristics. Environmental and Experimental Botany 65(2-3): 386-394.
6. BELLIS, D., K. SATAKE, K. TSUNODA& C. W. MCLEOD, 2003. Environmental monitoring of historical change in arsenic deposition with tree bark pockets. Journal of Environmental Monitoring 5: 671-674.
7. BRUNO, G., L. STIEFKENS, M. HADID, I. LISCOVSKY, M. T. COSA& N. DOTTORI, 2007. Efecto de la contaminación ambiental en la anatomía de la hoja de Ligustrum lucidum (Oleaceae). Boletín de la Sociedad Argentina de Botánica 42(3-4): 231-236.
8. CASTRO, M. A. & M. A. FAGGI (Comp.), 2008. Bioindicadores de la contaminación ambiental. Editorial Dunken, pp. 159.
9. CASTRO, M. A., L. E. NAVARRO, P. PERELMAN & A. FAGGI, 2008. Análisis multielemental de cortezas de árboles de Fraxinus pennsylvanica a lo largo de un gradiente urbano-periurbano en el área metropolitana de Buenos Aires, Argentina. Estimación de la polución atmosférica utilizando la técnica ICP-OES. In: Castro, M. A. & M. A. Faggi (Comp.), Bioindicadores de la contaminación ambiental. Editorial Dunken, pp. 159.
10. D’AMBROGGIO DE ARGÜESO, A. 1986. Manual de técnicas en histología vegetal. Ed. Hemisferio Sur. Buenos Aires.
11. DE TEMMERMAN, L., J. N. B. BELL, J. P. GARREC, A. KLUMPP, G. M. H. KRAUSE & A. E. G. TONNEIJCK, 2004. Biomonitoring of air pollution with plants. A consideration for the future. In: Klumpp, A., W. Ansel, & G. Klumpp (Eds.), Urban air pollution, Bioindicators and Environmental Awareness (pp. 337- 373). Cuvillier Verlag, Gottingen.
12. DION, M., S. LORANGER, G. KENNEDY, F. GOURCHESNE & J. ZAYED, 1993. Evaluation of black spruce (Picea mariana) as a bioindicator of aluminion contamination. Water, Air and Soil pollution 71: 29-41.
13. FAGGI, A., P. PERELMAN, M. A. CASTRO, & E. MARTINEZ-CARRETERO, 2008. Caso Mendoza, Area Metropolitana de Mendoza. In: Castro, M. A. & M. A. Faggi (Comp.), Bioindicadores de la contaminación ambiental (pp. 41-49). Editorial Dunken.
14. FUJIWARA, F., M. DOS SANTOS, J. MARRERO, G. POLLA, D. GÓMEZ, L. DAWIDOWSKI & P. SMICHOWSKI, 2006. Fractionation of eleven elements by chemical bonding from airborne particulate matter collected in an industrial city of Argentina. Journal of Environmental Monitoring 8: 1-10.
15. JAMES, F. C. & C. E. MCCULLOCH, 1990. Multivariate analisys in ecology and systematic: panacea or Pandora’s box? Annual Review of Ecology and Systematic 21: 129-166.
16. JASAN, R., T. VERBURG, T. WOLTERBEEK, R. PLÁ & M. PIGNATA, 2004. On the use of the lichen Ramalina celastri (Spreng.) Krong. & Swincs. as an indicator of atmospheric pollution in the province of Córdoba, Argentina, considering both lichen physiological parameters and element concentrations. Journal of Radioanalytical and Nuclear Chemistry 259: 93.
17. KABATA-PENDIAS, A. & H. PENDIAS, 2001. Trace elements in soils and plants. CRC Press, Florida, USA.
18. LINSTER, M., 1991. The impact of sewage sludge on agriculture. Towards Sustainable Agricultural Development 11: 320- 336.
19. LÓPEZ, S. C., P. PERELMAN, M. RIVARA, M. A. CASTRO & A. FAGGI, 2006. Soil characteristics and metal concentration along an urbanisation gradient in Buenos Aires, Argentina. Multequina 15: 69-80.
20. MANDRE, M., J. RAUK & K. OTS, 1995. Needle and shoot growth. In: Mandre, M. (Ed.), Dust pollution and forest ecosystems. A study of conifers in an alkalized environment. Publication 3 (pp. 103- 111). Tallinn: Institute of Ecology.
21. MASUCH, G., J. T . FRANZ, H. G. KICINSKI & A. KETTRUP, 1992. Histological and biochemical differences of lightly and severely injured spruce needles of two stands in northrhine Westphalia. Environmental and Experimental Botany 32: 163-182.
22. MHATRE, G. N. & S. B. CHAPHEKAR, 1985. The effect of the mercury on some aquatic plants. Environmental Pollution Ser. A. 39: 207-216.
23. MOHAPATRA, A. & A. K. PANIGRAHI, 1991. Effect of mercuric chloride on the pigment content of a mulberry plant. Pollution Research 10(3): 123-133.
24. OCHIAI, E. I., 1987. General Principles of Biochemistry of the Elements. Plenum Press, New York.
25. PÄÄNKKÖNEN, E., T. HOLOPAINEN& L. KÄRENLAMPI, 1995. Ageing-related anatomical and ultraestructural changes in leaves of Birch (Betula pendula Roth.) clones as affected by low ozone exposure. Annals of Botany 75: 285-294.
26. PÄÄNKKÖNEN, E., T. HOLOPAINEN& L. KÄRENLAMPI, 1997. Differences in growth, leaf senescence and injury, and stomatal density in Birch (Betula pendula Roth.) in relation to ambient levels of ozone in Finland. Environmental Pollution 96: 117-127.
27. PANICHEV, N. & R. I. MCCRINDLE, 2004. The application of bioindicators for the assessment of air pollution. Journal of Environmental Monitoring 6: 121-123.
28. PATRA, M. & A. SHARMA, 2000. Mercury toxicity in plants. Botanical Review 66: 379-422.
29. PERELMAN, P. E., M. A. CASTRO, L. E. NAVARRO, M. RECCHI, M. ARRIAGA, S. LÓPEZ, E. MARTÍNEZ CARRETERO & A. FAGGI, 2006. Análisis elemental de cortezas de fresno (Fraxinus pennsylvanica) a lo largo de un gradiente urbano-periurbano en la metrópolis de Buenos Aires. Revista del Museo Argentino de Ciencias Naturales, n.s., 8(2): 231-236.
30. PIGNATA, M. L., G. L. GUDIÑO, E. D. WANNAZ, R. R. PLA, C. M. GONZÁLEZ, H. A. CARRERAS& L. ORELLANA, 2002. Atmospheric quality and distribution of heavy metals in Argentina employing Tillandsia capillaris as a biomonitoring. Environmental Pollution 120: 59-68.
31. RAI, V., S. KHATOON, S. S. BISHT& S. MEHROTRA, 2005. Effect of cadmium on growth, ultramorphology of leaf and secondary metabolites of Phyllanthus amarus Schum. and Thonn. Chemosphere 61(11): 1644-50.
32. RANI, M., N. PAL & R. K. SHARMA, 2006. Effect of railway engines emission on the micromorphology of some field plants. Journal of Environmental Biology 27(2): 373-376.
33. RAO, M. V. & P. S. DUBEY, 1991. Detoxifying mechanism: Probable role in determining the plant response to SO2 under different light intensities. Indian Journal of Environment and Toxicology 1(1): 39-45.
34. ROHLF, F. J., 1992. NTSYS-Pc numerical taxonomy and multivariate analysis system. Applied Biostatistics Inc. New York.
35. SALGARE, S. A. & V. B. THORAT, 1990. Effect of auto-exhaust pollution at Andheri (west) Bombay on the micromorphology of some trees. Journal of Ecobiology 2(4): 267-272.
36. SANT’ANNA-SANTOS, B. F., L. CAMPOS DA SILVA, A. A. ACEVEDO & R. AGUILAR, 2008. Effects of simulated acid rain on leaf anatomy and micromorphology of Genipa americana L. (Rubiaceae). Brazilian Archives of Biology and Technology 49: 313-321.
37. SHARMA, G. K. & J. TYREE, J., 1973. Geographic leaf cuticular and gross morphological variations in Liquidambar styraciflua L. and their possible relationship to environmental pollution. Botanical Gazette 134: 179-184.
38. SETIA, P., N. CHILDS, E. WAILES& J. LIVEZEY, 1994. The U. S. Rice Industry. Agricultural Economic Report No. 700, U.S. Department of Agriculture Economic Research 8 Service, Washington, D.C.
39.SNEATH, P. H. A. & R. R. SOKAL, 1973. Numerical taxonomy. The principles and practice of numerical classification. Ed. W. H. Freeman. San Francisco, CA.
40. SRIVASTAVA, K., 1999. Celosia as indicator of automobile pollution. In XVI International Botanical Congress, Abstracts (pp. 2580). St. Louis, USA.
41. STATSOFT, INC., 1995. Statistica for Windows [Computer program manual]. Tulsa, OK, USA.
42. STEUBING, L., R. GODOY & M. ALBERDI, 2002. Métodos de ecología vegetal 2002. Editorial Universitaria, Santiago de Chile, 345 pp.
43. WALKENHORST, A., J. HEGEMEYER& W. BRECKLE, 1993. Passive monitoring of airborne pollutants, particularly trace metals, with tree bark. In: Market, B. (Ed.), Plants as biomonitors. Indicators for heavy metals in the terrestrial environment (pp. 523-540). VCH, Weinheim, New York, Basel, Cambridge.
Recibido 05/2014
Aceptado 07/2014