SciELO - Scientific Electronic Library Online

 
vol.63 issue3Mites associated with two species of the genusOdontotaenius (Coleoptera: Passalidae) in MexicoWater pollution and distribution of the black fly (Diptera: Simuliidae) in the Atlantic Forest, Brazil author indexsubject indexarticles search
Home Pagealphabetic serial listing  

Services on Demand

Journal

Article

Indicators

Related links

Share


Revista de Biología Tropical

On-line version ISSN 0034-7744Print version ISSN 0034-7744

Rev. biol. trop vol.63 n.3 San José Jul./Sep. 2015

 

Invertebrados terrestres

Fluctuating asymmetry in Apis mellifera (Hymenoptera: Apidae) as bioindicator of anthropogenic environments

Asimetría fluctuante en Apis mellifera (Hymenop-tera: Apidae) como bioindicador de ambientes antropo-génicos

Lorena Andrade Nunes 1  

Edilson Divino de Araújo 2  

Luís Carlos Marchini 3  

1Universidade Estadual do Sudoeste da Bahia, Pos-graduado em Genética, Biodiversidade e Conservado. Rua José Moreira Sobrinho, s/n, Jequiezinho, Jequié, BA, 45206-190, Brazil; lorenunes2@gmail.com

2Laboratorio de Genética e Conservado de Recursos Naturais-GECON, Departamento de Biologia Universidade Federal de Sergipe, Campus Universitário José Aloízio de Campos; edaraujo@yahoo.com.br

3Laboratorio de Apicultura Departamento de Entomologia e Acarologia, Escola Superior de Agricultura "Luiz de Queiroz"- ESALQ/USP; lcmarc@usp.br

ABSTRACT

The successful distribution of A. mellifera is due to their ability to adjust to seasonal variations, considerable control over their internal physical environment and exploration of different resources. However, their populations have experienced different forms and levels of environmental pressure. This research aimed to verify the phenotypic plasticity in both size and shape of wings in A. mellifera using fluctuating asymmetry, based on geometric morphometrics from apiaries located in sites with high and low levels of anthropization. We sampled 16 locations throughout all five geographic regions of Brazil. At each site, samples were collected from 20 beehives installed in apiaries: 10 installed near high anthropogenic environments (Cassilàndia - MS, Fortaleza - CE, Maringá - PR, Aquidauana - MS, Rolim de Moura - RO, Riachuelo - SE, Ubirata - PR and Piracicaba - SP), and 10 in sites with low levels of human disturbance (Cassilàndia - MS, Itapiúna CE, Uniao da Vitoria - PR, Aquidauana - MS, Rolim de Moura - RO, Pacatuba - SE, Erval Seco - RS, Rio Claro - SP). A sample of 10 individuals was taken in each hive, totaling 200 per location, for a total of 1 600 individuals. We used fluctuating asymmetry (FA) in size and shape of the forewing through geometric morphometrics. The FA analysis was conducted in order to check bilateral differences. The indexes of size and shape were submitted to analysis of variance (ANOVA), where the characters evaluated were used as factors to verify the size and shape differences. The results indicated an asymmetry on the shape of the wing (P < 0.001) but no asymmetry was observed on wing size. Considering FA as an environmental response and high and low impacted areas as a fixed factor, we observed significant differences (P < 0.05). The results for the wing shape in A. mellifera demonstrated that this feature undergoes more variation during ontogeny compared to the variation in size. We concluded that bee samples collected from colonies with higher levels of human disturbance had higher wing-shape asymmetry; the variation of fluctuating asymmetry in the wing shape of honeybees can be used as an indicator of the degree of environmental anthropization.

Key words: centroid size; forewings; geometric morphometrics; honeybees; shape; anthropization

RESUMEN

La distribución exitosa de A. mellifera se debe a su capacidad para adaptarse a las variaciones estacionales, controlar considerablemente su ambiente físico interno y por la exploración de recursos. Sin embargo, sus poblaciones experimentan diferentes formas y niveles de presión ambiental. Esta investigación evaluó colmenares, ubicadas en entornos con distintos niveles de antropización. Ambientes considerados altamente antropogénicos fueron escogidos: áreas urbanas, agrícolas con aplicaciones de insecticidas, y las industrias con grandes áreas (regiones con altos índices de contaminación). Por otra parte, los entornos considerados con bajos niveles de perturbación humana fueron elegidos: áreas de preservación permanente de los bosques restantes (Mata Atlántica, Planalto, Cerrado), áreas ribereñas o reservas ecológicas. Se muestrearon 16 localidades que abarcan las cinco regiones del Brasil. Estamos utilizando la asimetría fluctuante desde el tamaño y la forma del ala anterior por técnicas de morfometría geométrica. En cada sitio, las muestras se obtuvieron de las 20 colmenas instaladas en los colmenares: 10 instalados cerca de entornos altamente antropogénicos (Cassilándia - MS, Fortaleza - CE, Maringá - PR, Aquidauana - MS, Rolim de Moura - RO, Riachuelo - SE, Ubirata - PR y Piracicaba - SP) y 10 en sitios de bajo nivel de perturbación humana (Cassilándia - MS, ItapiúnaCE, Uniao da Vitória - PR, Aquidauana - MS, Rolim de Moura - RO, Pacatuba- SE, Erval Seco - RS, Rio Claro - SP). Los índices de tamaño y forma fueron sometidos a análisis de varianza (ANOVA), donde se utilizaron los caracteres evaluados como factores, para verificar las diferencias de tamaño y forma. La asimetría de análisis fluctuante (AF) se llevó a cabo con el fin de comprobar las diferencias bilaterales. Los resultados indican la existencia de la asimetría de la forma del ala (P < 0.001), pero no se observó asimetría del tamaño del ala. Considerando AF como respuesta ambiental y áreas de alto y bajo grado de alteración humana como factor fijo, observamos diferencias significativas (P < 0.05). Los resultados, para la forma de ala de la A. mellifera, muestran que esta característica se somete a más variación durante la ontogenia en comparación con la variación en el tamaño. Llegamos a la conclusión de que las colonias de abejas recogidas en ambientes con niveles más altos de la perturbación humana tienen una mayor asimetría en forma de ala, por lo que la asimetría fluctuante en forma de alas de las abejas puede ser utilizada como un indicador del grado de antropización del medio ambiente.

Palabras clave: tamaño de centroide; alas anteriores; mor-fometría geométrica; abejas; forma; antropización

During embryonic development, the organism may face disruption from a genetic or environmental nature that needs to surmount for the appropriate phenotypic expression. This instability can be assessed by random disturbances suffered during development, enabling verification of increased phenotypic variance within one species (Knierim et al, 2007). The intense environmental changes resulting from human activities, such as deforestation and pollution, create increasing concerns about the ecological consequences of these human activities on natural populations (Polak, Opoka, &amp; Cartwright, 2002). Moreover, genetic changes such as mutation, selection, inbreeding, and chromosomal disorder can increase asymmetry (Parsons, 1992).

Changes in the individuals' development, in natural populations, show up as promising bioindicators of environmental quality. Most part of the environment effects on the pheno-typic responses is due to phenotypic plasticity (Scheiner, 1993). Plasticity is essential in adapting the organism to heterogeneous and unstable environments, since it allows the organism dynamic maintenance, although not all plastic responses are adaptive (Gotthard &amp; Nylin, 1995; Via et al., 1995).

The fluctuating asymmetry (FA) is commonly used to evaluate subtle random and non-directional changes between the symmetry planes of the individuals during their onto-genetic development. FA has been studied because they reflect both genetic and environmental stresses (Palmer &amp; Strobeck, 1986; Parsons, 1990).

The occurrence of asymmetric fluctuations is more related to organisms that develop under stress, because they generally have greater difficulty for regulating their development. Therefore, the fluctuating asymmetry can be used as a tool to estimate stress, since it represents a sensitive indicator of organism development, being helpful to elucidate how organisms respond to environmental and genetic changes (Lempa et al., 2000; Silva, Lomonaco, Augusto, &amp; Kerr, 2009).The animals that exhibit bilateral symmetry, as the insects, usually have similar morphological structures on both sides of the body. Since the development of these structures is under genetic control, it is expected that they exhibit morphological symmetry on both sides, as the products of expression of the same genome (Leary &amp; Allendorf, 1989). Bees have a great plasticity to occupy different niches and, on that basis, are able to establish themselves as wild populations in different regions. Since the beginning of the Africanization process of Apis mellifera Linnaeus in the Americas, this species has been the subject of many studies, making it possible to observe its great adaptability and wide distribution, and enabling the maintenance of a diverse gene pool (Nunes, Araujo, Marchini, &amp; Moreti, 2012). The successful distribution of A. mel-lifera is due to their ability to adjust to seasonal variations and apply considerable control over their internal physical environment and by the resource exploration (Oliveira &amp; Cunha, 2005). However, their populations experience different forms and levels of environmental pressure.

This study aimed to verify the phenotypic plasticity in both size and shape of wings in A. mellifera using fluctuating asymmetry based on geometric morphometrics. Furthermore, we evaluated whether adverse conditions related to human actions affect the bilateral symmetry of wings or not.

Materials and Methods

The study was conducted in a two years period, from April 2009 to April 2011. We sampled 16 locations in Brazil, distributed in 13 municipalities (Table 1), covering five geographic regions of Brazil. At each site, samples were collected from the 20 beehives installed in apiaries: 10 installed near highly anthropogenic environments (Cassilandia - MS, Fortaleza -CE, Maringa - PR, Aquidauana - MS, Rolim de Moura - RO, Riachuelo - SE, Ubirata - PR and Piracicaba - SP), and 10 in places considered of low levels of human disturbance (Cassilan-dia - MS, Itapiuna - CE, Uniao da Vitoria - PR, Aquidauana - MS, Rolim de Moura - RO, Pacatuba - SE, ErvalSeco - RS, Rio Claro - SP). A total of 10 individuals were sampled in each hive using falcon tubes, totaling 200 individuals per location, and 1 600 when both high and low impacted areas were considered. According to Palmer (1994), a total of 30 individuals are enough to verify the presence of FA.

The environment was chosen in a way to maximize the degree of conservation and anthropization. Regarding the classification of places with high and low degree of human disturbance, environments considered highly anthropogenic were chosen for major urban, farming areas with insecticide applications, and industries with large areas (regions with high pollution rates). On the other hand, environments considered with low levels of human disturbance were chosen because they are areas of permanent preservation of remaining forests (Atlantic Forest, Planalto, Cerrado), riparian areas or ecological reserves (Table 1).

Table 1: Locations, geographic coordinates, characterization of the sample collection area and kind of impact of Apis mellifera population 

All sampled specimens were stored in ethanol at -20 °C. Later, the forewings of both right and left sides of each individual were removed, with the help of entomological forceps. Subsequently, these structures were placed between two microscopic slides, to keep them flat and prevent distortion, during the image capture process. The images were acquired and digitized with a camera attached to a stereomicroscope using the Leica Application Suite version 3.4.1. From the images of the wings, coordinates of the landmarks were obtained on both sides using software Tps-DIG2 version 1.40 (Rohlf, 1998), as shown in figure 1.

Four hundred individuals from four locations were selected for the asymmetric preliminary analysis, in order to verify the measurement error. Hence, each individual was measured twice. According to Palmer (1994), tests such as this should be used, to verify the measurement error, using at least 30 individuals.

The indexes of size were submitted to analysis of variance (ANOVA), where the characters evaluated were used as factors, to verify the size differences. The fluctuating asymmetry (FA) analysis was conducted in order to check bilateral differences. To calculate FA, both sides of the wing were measured, as well, the mean difference between the right and left sides, i.e., FA = [(E|(R-L)|/n], and compared using the paired T-test, to see whether the distributions have means equal to zero (Palmer &amp; Strobeck, 1986; Swaddle, Witter, &amp; Cuthill, 1994). FA was measured in all individuals for the determination of size and shape.

The univariate analysis of variance, ANOVA, was performed, to verify the asymmetry of the wing centroid size, where the centroid size was treated as an independent variable, the side of the body as a fixed effect, and the individuals as a random effect (Klingenberg &amp; McIntyre, 1998). To obtain the F value for the side effect and individuals, we used the interaction between individual and side as the denominator. To obtain the F value for the interaction between individual and side, the measurement error as denominator was used (Palmer, 1994; Klingenberg &amp; Mclntyre, 1998).

According to the genetic model of Via and Lande (1985), which evaluates quantitatively the phenotypic plasticity, the total phenotypic variation of a character is the sum of its genetic and environmental components, as well the interaction between these factors. Therefore, the occurrence of phenotypic plasticity in size has been estimated by the environmental component of variation (Falconer, 1989).

Figure 1: Forewing of Apis mellifera with anatomical landmarks used for analyzes. 

Procrustes ANOVA analysis was performed to analyze the asymmetry of the wing shape of the populations and the configurations of the species aligned (Klingenberg &amp; Mclntyre, 1998; Palmer &amp; Strobeck, 2003). In a similar way, we analyzed the asymmetry to the centroid size, where the side of the wings was treated as a fixed effect, and the individuals as random effect. Thus, the effect of individual consists in individual variation of the shape, and the side and individual effect corresponding to FA. The term of the interaction between side and individual is the measure of fluctuating asymmetry (small random differences between left and right sides, between individuals) and the residue, which is the variability between replicates, are treated as the measure of error (Klingenberg &amp; McIntyre, 1998).

Procrustes distance analysis was performed using the data obtained from the wing to assess the differences and the asymmetry level between the various localities sampled.

With the Procrustes distances values, we calculated the heritability of wing-shape variation in A. mellifera using the multiple variance analysis model by Goodall (1991) and intraclass correlation coefficient (Sokal &amp; Rohlf, 1995) as reported by Monteiro et al. (2002). The causes for fluctuating asymmetry have been studied, and it was believed that it was completely related to the environment (Palmer, 1994); however, it has a genetic basis that remains unknown with the heritability close to zero (Palmer &amp; Strobeck, 1986; Parsons, 1990; Fuller &amp; Houle, 2002; Leamy &amp; Klingemberg, 2005).

Results

The effect of the ANOVA for the size of the wing, using the centroid size, was not significant for populations from Uniao da Vitoria-PR, Rolim de Moura-RO, (environments with low levels of human disturbance, P = 0.515), Maringa-PR and Piracicaba-SP (highly anthropogenic environments, P = 0.32). The other locations showed significant differences in wing size, in the differences between left and right sides (L-R) and in the effect of interaction between side and individual. In the ANOVA, the centroid size made it possible to verify that some populations, of A. mellifera, showed the presence of directional asymmetry.

On the other hand, the effects of ANOVA, on the wing-shape asymmetry, showed significant results in all populations, regardless of where the samples were obtained and whether they presented a high or low level of disturbance. This result demonstrated the existence of fluctuating asymmetry in all populations of A. mellifera, when wing shape was evaluated. From the interaction analysis performed, to verify the presence of FA and to compare the populations analyzed, we observed that colonies collected close to where insecticide application occurs, such as Maringa-PR, Riachuelo-SE and Ubirata-PR, presented a greater asymmetry value, when compared to localities with low levels of human disturbance (Fig. 2). Furthermore, we found that samples collected from locations, classified with a high degree of human disturbance, such as Aquidauna-MS and Rolim de Moura-RO, also showed a higher degree of asymmetry, compared to the colonies sampled in environments with a low degree of human disturbance. Thus, from all 16 locations sampled, 63 % of the colonies collected in contaminated environments, showed greater asymmetry when compared to others. Significant differences (P < 0.05) were observed when FA was regarded as the environmental response and high and low impacted areas were considered a fixed factor. The only samples that observed a higher value of FA, were those populations located in areas considered with low levels of human disturbance, such as Cassilandia-MS, Itapiuna-CE and Rio Claro-SP. With regards to wing shape heritability, a value of h2 = 0.33 with EPM standard error = 0.0008 was found.

The presence of varying degrees of FA, in the regions studied by evaluating the wing shape of A. mellifera, enables to see that these differences may also be related to the geographic distance between colonies (P = 0.0143 and r = 0.29). This is expected since geographical proximity is associated with greater genetic and climate similarity between colonies.

Discussion

The differences between individuals are often represented using size variation between them. However, according to Schneider, Leamy, Lewis and DeGrandi-Hoffman (2003), size and some structures analysis may have some restrictions for assessing the asymmetry, and may not be good indicators for the analysis of asymmetry in bees. It is consistent with the results found for the size asymmetry in some populations of A. mellifera (Uniao da Vitoria-PR, Rolim de Moura-RO, (P > 0.05), Maringa-PR and Piracicaba-SP). Clarke &amp; Oldroyd (1996), studied males and workers of A. mel-lifera and observed that some characters do not show significant size correlation, usually because the wing size in insects is correlated with the overall size of the body. Therefore, it is possible that identical genotypes exhibit differences associated with the use of nutritional resources and with different environmental conditions. We observed the presence of directional asymmetry for wing size in some populations, a result that corroborates Smith, Crespi and Bookstein (1997) and Schneider et al. (2003), that evaluated hybrid and pure populations of A. mellifera from Procrustes superimposition and symmetry analysis.

Figure 2: Interaction of the presence of fluctuating asymmetry in wing of Apis mellifera in different collection sites, where L are areas that have a low degree of human disturbance, and H are highly anthropogenic environments. 

It is possible to verify the presence of FA in the wing-size in some populations of A. mel-lifera, similar to that found by Clarke, Oldroyd and Hunt (1992), which evaluated asymmetry size between males and females of A. mellifera, and confirmed that inbreeding has no effect on the stability in the development of this species. They also found that in the haplodiploid system, the level of heterozygosity does not seem to be an important factor in the stability of the development.

In the analyses performed for this study, we observed that the shape represents better symmetry differences when compared to size, which are similar to results found by Del Lama, Gruber and Godoy (2002). Those authors also found that the shape of the wing contain better structures to asymmetry study, because it apparently is more affected by the environment when compared with the size. According to Leary and Allendorf (1989), a single feature can provide a good index of stress suffered during the individuals' development. All populations showed significant differences between left and right sides, as well as in the interaction between side and individual, suggesting a directional asymmetry and fluctuating asymmetry.

The results for the shape of the wing of A. mellifera, demonstrate that this feature undergoes more variation during ontogeny compared to the variation in size. These results are consistent with those found by Smith et al. (1997) that noted the presence of fluctuating and directional asymmetry on the shape of the wing of A. mellifera, and affirmed that directional asymmetry is more common than expected, and suggested that this asymmetry is related to the position of the larvae, prepupae and pupae in the brood cell. In addition, studies performed by Abou-Shaara and Al-Ghamdi (2012) showed asymmetric patterns on wings of Carniolan (A. mellifera carnica) and Yemeni honeybees (A. mellifera jementica), using conventional and geometric morphometric analysis, suggesting that temperature, seasons and viruses may have contributed to the presence of asymmetry.

The reason why the locations with low degree of human disturbance (Cassilandia-MS, Itapiuna-CE, and Rio Claro-SP) have shown greater asymmetry may be associated with various environmental factors that generate noise during development. In view of the fact that individuals were collected in forest fragments, could be the factor that is influencing this asymmetry. Studies by Lens, Dongen, Wilder, Brooks and Matthysen (1999) evaluating species of birds found in forest fragments, and comparing them with the birds deposited in museums (collected more than 50 years ago, before the fragmentation of the forest), found an increase in FA. That is, fragmentation can cause stress on these organisms, which may affect their development. A similar cause should be occurring in bee populations collected in environments recognized as preserved.

The greater asymmetry observed in places with insecticide application, corroborates studies by Abaga et al. (2011): a high level of developmental instability, in A. mellifera, was significantly associated with the use of insecticides in cotton farms; it was verified that the FA index, in the tarsus length, was increased, implying that FA can be considered as a biomarker that reflects the stress induced by insecticide treatments. Another factor that may contribute to the existence of FA is that all individuals have been collected in apiaries and the handling could cause changes in behavior, and be a stress generator.

The presence of different degrees in asymmetry, in the regions studied, is similar to that found by Mazeed (2011), evaluating wing patterns in Egyptian (A. mellifera lamarckii) and Carniolan honeybees (A. mellifera carnica). They verified the presence of fluctuating and directional asymmetry, mainly in Carniolan honeybees (as they are not found naturally in Egypt), and therefore, should not be adapted to the region. This suggests that the geographical location and adaptation to different environmental conditions, to which each organism is subjected, affect the degree of asymmetry. Furthermore, this same study draws the attention to the genetic differences between lineages and inbreeding as the possible causes of asymmetry.

The asymmetry exists technically in different scales and all structures will be asymmetrical to some extent (Palmer, 1996). The low heritability observed in the wings of A. mellifera was expected, since FA is a random variation that is related to the environment (Palmer, 1994) to which the genetic basis remains unknown with heritability close to zero (Leamy &amp; Klingenberg, 2005). The heritability of a trait is a parameter that varies in time and space, and the low heritability can be indicative of a strong constraint on development and evolution (Monteiro, Diniz-Filho, dos Reis &amp; Araújo, 2002). These authors also, found significant heritability values for the shape of the forewing in A. mellifera (h2 = 0.2935; P < 0.001), a similar result was found in the present study (h2 = 0.33).

It seems to us that the degree of human disturbance can increase the intensity of fluctuating asymmetry, however, different types of human disturbance may result in different pressures of asymmetry. Environmental factors, location of the apiary and management, can directly influence the results. Studies by Silva, Lomonaco, Augusto and Kerr (2009) found that environmental factors and anthropogenic interference influence the developmental stability in the size of Euglossa pleosticta.

The adjustment to the habitat is the response that causes asymmetry. Populations located in urban areas tend to be more negatively affected and present a greater FA than others, as noted by Weller and Ganzhorn (2004) in populations of carabid beetles, which showed that specimens collected closer to the cities showed greater FA levels. Weller and Ganzhorn (2004) concluded that the urban effect makes the communities move out, causing a predicted loss of species, a reduction in body size, and to the FA, an increase in species that are susceptible to urbanization. Those results are similar to the ones found in the present work, where populations of A. mellifera, collected in urban areas (Aquidauana-MS and Rolim de Moura-RO), showed greater asymmetry effects when compared with populations of the same localities, inferring the contribution of urbanization in the presence of asymmetry in these organisms. Therefore, understanding the influences of human impacts on urban landscapes is important to reduce the negative effects on biodiversity remaining (Pickett et al., 2001).

In this study it was observed that the shape of the wing is more suitable for the study of asymmetry when compared with the size of the wing. We observed that the colonies collected in highly anthropogenic areas exhibit greater asymmetry compared with those collected in environments with low levels of human disturbance. Indeed, morphological characters may also indicate habitat quality. However, since it is a pioneering and still prospective study, the use of genetically closely related colonies and a better description of the environmental variables, can significantly improve the accuracy of the fluctuating asymmetry tests in environmental quality assessment, using the shape of the wings of Apis mellifera.

Acknowledgments

We thank Leandro Rabello Monteiro (UENF) for helpful discussions and suggestions and the Fundacao de Amparo á Pesquisa do Estado de Sao Paulo (FAPESP) and Con-selho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) for the financial support. We also thank the apiarists for supplying the biological material for the study. We thank Sona Jain for the English revision of the manuscript.

REFERENCES

Abaga, N. O. Z., Alibert, P., Dousset, S., Savadogo, P. W., Savadogo, M., &amp; Sedogo, M. (2011). Insecticide residues in cotton soils of Burkina Faso and effects of insecticides on fluctuating asymmetry in honey bees (Apis mellifera Linnaeus). Chemosphere 834585-592. [ Links ]

Abou-Shaara, H. F., &amp; Al-Ghamdi, A. A. (2012). Studies on wings symmetry and honey bee races discrimination by using standard and geometric morphometrics. Biotechnology in Animal Husbandry 28(3), 575-584. doi:10.2298/BAH1203575A [ Links ]

Clarke, G. M., &amp; Oldroyd, B. P. (1996). The genetic basis of developmental stability in Apis mellifera II. Relationships between character size, asymmetry and single-locus heterozygosity. Genetica 97(2), 211224. doi:10.1007/BF00054628 [ Links ]

Clarke, G. M., Oldroyd, B. P., &amp; Hunt, P. (1992). The genetic basis of developmental stability in Apis mellifera: Heterozygosity versus Genic Balance. Evolution 46753-762. [ Links ]

Del Lama, M. A., Gruber, C. V., &amp; Godóy, I. C. de. (2002). Heterozigosidade e assimetria do número de hámulos em operárias adultas de Apis mellifera (Hymenoptera, Apidae). Revista Brasileira de Entomologia 46(4), 591-595. doi:10.1590/S0085-56262002000400014 [ Links ]

Falconer, D. S. (1989). Introduction to Quantitative Genetics (3rd edn, 438 p.). New York: Longmans Halow, Longman Sci e Tech. [ Links ]

Fuller, R. C., &amp; Houle, D. (2002) Inheritance of developmental instability. M. Polak Developmental Instability: Causes and Consequences 157-183). Oxford, United Kingdom Oxford University Press [ Links ]

Goodall, C. R. (1991). Procrustes methods in the statistical analysis of shape. Journal of the Royal Statistical Society B 53285-339. [ Links ]

Gotthard, K., &amp; Nylin, S. (1995). Adaptive Plasticity and Plasticity as an Adaptation: A Selective Review of Plasticity in Animal Morphology and Life History. Oikos 74(1), 3. doi:10.2307/3545669 [ Links ]

Klingenberg, C. P., &amp; Mclntyre, G. S. (1998). Geometric Morphometries of Developmental Instability: Analyzing Patterns of Fluctuating Asymmetry with Procrustes Methods. Evolution 52(5), 1363. doi:10.2307/2411306 [ Links ]

Knierim, U., Van Dongen, S., Forkman, B., Tuyttens, F. A. M., Spinka, M., Campo, J. L., &amp; Weissengruber, G. E. (2007). Fluctuating asymmetry as an animal welfare indicator - A review of methodology and validity. Physiology &amp; Behavior, 92(3), 398-421. doi:10.1016/j.physbeh.2007.02.014 [ Links ]

Leamy, L. J., &amp; Klingenberg, C. P. (2005). The genetics and evolution of Fluctuating Asymmetry. Annual Review of Ecology, Evolution, and Systematics -6(1), 1-21. doi:10.1146/annurev.ecolsys.36.102003.152640 [ Links ]

Leary, R. F., &amp; Allendorf, F. W. (1989). Fluctuating asymmetry as an indicator of stress: Implications for conservation biology. Trends in Ecology &amp;, Evolution 4(7), 214-217. doi:10.1016/0169-5347(89)90077-3 [ Links ]

Lempa, K., Martel, J., Koricheva, J., Haukioja, E., Ossipov, V., Ossipova, S., &amp; Pihlaja, K. (2000). Covariation of fluctuating asymmetry, herbivory and chemistry during birch leaf expansion. Oecologia 1223354360. doi:10.1007/s004420050041 [ Links ]

Lens, L., Dongen, S. V., Wilder, C. M., Brooks, T. M., &amp; Matthysen, E. (1999). Fluctuating asymmetry increases with habitat disturbance in seven bird species of a fragmented afrotropical forest. Proceedings of the Royal Society B: Biological Sciences, 266(1425), 1241-1246. doi:10.1098/rspb.1999.0769 [ Links ]

Mazeed, A. M. M. (2011). Anomalies and asymmetry of wing venation pattern in Carniolan and Egyptian bee populations in Egypt. Egyptian Academic Journal of Biological Sciences 4(1), 149-161. [ Links ]

Monteiro, L. R., Diniz-Filho, J. A. F., dos Reis, S. F., &amp; Araûjo, E. D. (2002). Geometric estimates of heri-tability in biological shape. Evolution 56(3), 563. doi:10.1554/0014-3820(2002)056[0563:GEOHIB]2. 0.CO;2 [ Links ]

Nunes, L. A., Araûjo, E. D. de, Marchini, L. C., &amp; Moreti, A. C. de C. C. (2012). Variation morphogeome-trics of Africanized honey bees (Apis mellifera) in Brazil. Iheringia. Série Zoologia 1023321-326. doi:10.1590/S0073-47212012005000002 [ Links ]

Oliveira, M. L. de, &amp; Cunha, J. A. (2005). Abelhas africanizadas Apis mellifera scutellata Lepeletier, 1836 (Hymenoptera: Apidae: Apinae) exploram recursos na floresta amazônica? Acta Amazonica, -5(3), 389394. doi:10.1590/S0044-59672005000300013 [ Links ]

Palmer, A. R. (1994). Fluctuating Asymmetry Analyses: A Primer. In T. A. Markow (ed.). Developmental Instability: its Origins and Evolution ary Implications (pp. 335-364). Dordrecht, The Netherlands: Kluwer Academic Publishers. [ Links ]

Palmer, A. R. (1996). Waltzing with Asymmetry: Is fluctuating asymmetry a powerful new tool for biologists or just an alluring new dance step? BioScience 46(7), 518-532. doi:10.2307/1312930 [ Links ]

Palmer, A. R., &amp; Strobeck, C. (1986). Fluctuating Asymmetry: Measurement, Analysis, Patterns. Annual Review of Ecology and Systematics 171391-421. doi:10.1146/annurev.es.17.110186.002135 [ Links ]

Palmer, R. A., &amp; Strobeck, C. (2003) Fluctuating asymmetry analyses revisited. In M. Polak (Ed.) Developmental instability (DI): causes and consequences (pp. 279-319). Oxford University Press. [ Links ]

Parsons, P. A. (1990). Fluctuating Asymmetry: an epige-netic measure of stress. Biological Reviews 65(2), 131-145. doi:10.1111/j.1469-185X.1990.tb01186.x [ Links ]

Parsons, P. A. (1992). Fluctuating asymmetry: a biological monitor of environmental and genomic stress. Heredity 68(4), 361-364. doi:10.1038/hdy. 1992.51 [ Links ]

Pickett, S. T. A., Cadenasso, M. L., Grove, J. M., Nilon, C. H., Pouyat, R. V., Zipperer, W. C., &amp; Costanza, R. (2001). URBAN ECOLOGICAL SYSTEMS: Linking Terrestrial Ecological, Physical, and Socioeconomic Components of Metropolitan Areas1. Annual Review of Ecology and Systematics 32(1), 127-157. doi:10.1146/annurev.ecolsys.32.081501.114012 [ Links ]

Polak, M., Opoka, R., &amp; Cartwright, I. L. (2002). Response of fluctuating asymmetry to arsenic toxicity: support for the developmental selection hypothesis. Environmental Pollution 118(1), 19-28. doi:10.1016/ S0269-7491(01)00281-0 [ Links ]

Rohlf, F. J. (1998). On Applications of Geometric Morphometrics to Studies of Ontogeny and Phylogeny. Systematic Biology 47(1), 147-158. doi: [ Links ]

Scheiner, S. M. (1993). Genetics and Evolution of Phe-notypic Plasticity. Annual Review of Ecology and Systematics 24(1), 35-68. doi:10.1146/annurev.es.24.110193.000343 [ Links ]

Schneider, S. S., Leamy, L. J., Lewis, L. A., &amp; DeGrandi-Hoffman, G. (2003). The influence of hybridization between african and european honeybees, Apis melli-fera, on asymmetries in wing size and shape. Evolution 57(10), 2350. doi:10.1554/02-609 [ Links ]

Silva, M. C., Lomonaco, C., Augusto, S. C., &amp; Kerr, W. E. (2009). Climatic and anthropic influence on size and fluctuating asymmetry of Euglossine bees (Hyme-noptera, Apidae) in a semideciduous seasonal forest reserve. Genetics and Molecular Research 8(2), 730737. doi:10.4238/vol8-2kerr037 [ Links ]

Smith, D. R., Crespi, B. J., &amp; Bookstein, F. L. (1997). Fluctuating asymmetry in the honey bee, Apis mellifera: effects of ploidy and hybridization. Journal of Evolutionary Biology 10(4), 551-574. doi:10.1046/j.1420-9101.1997.10040551.x [ Links ]

Sokal, R. R., &amp; Rohlf, F. J. 1995. Biometry: the principles and practice of statistics in biological research (3rd ed., p. 887). New York: Freeman. [ Links ]

Swaddle, J. P., Witter, M. S., &amp; Cuthill, I. C. (1994). The analysis of fluctuating asymmetry. Animal Behaviour 484986-989. doi:10.1006/anbe.1994.1327 [ Links ]

Via, S., Gomulkiewicz, R., De Jong, G., Scheiner, S. M., Schlichting, C. D., &amp; Van Tienderen, P. H. (1995). Adaptive phenotypic plasticity: consensus and controversy. Trends in Ecology &amp;, Evolution 105212-217. doi:10.1016/S0169-5347(00)89061-8 [ Links ]

Via, S., &amp; Lande, R. (1985). Genotype-environment interaction and the evolution of phenotypic plasticity. Evolution 39505-522. [ Links ]

Weller, B., &amp; Ganzhorn, J. U. (2004). Carabid beetle community composition, body size, and fluctuating asymmetry along an urban-rural gradient. Basic and Applied Ecology 5(2), 193-201. doi:10.1078/1439-1791-00220. [ Links ]

Received: September 16, 2014; Revised: February 20, 2015; Accepted: March 23, 2015

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License