107 Population Divergence of Melipona scutellaris (Hymenoptera: Meliponina) in Two Restricted Areas in Bahia, Brazil by Cristovam Alves de Lima Junior1, Carlos Alfredo Lopes de Carvalho1, Lorena Andrade Nunes2 & Tiago Maurício Francoy3 ABSTRACT Melipona scutellaris Latreille has great economic and ecological importance, especially because it is a pollinator of native plant species. Despite the impor- tance of this species, there is little information about the conservation status of their populations. The objective of this study was to assess the diversity in populations of M. scutellaris coming from a Semideciduous Forest Fragment and an Atlantic Forest Fragment in the Northeast Brazil, through geometric morphometric analysis of wings in worker bees. In each area, worker bees were collected from 10 colonies, 10 workers per colony. To assess the diversity on the right wings of worker bees, 15 landmarks were plotted and the measures were used in analysis of variance and multivariate analysis, principal component analysis, discriminant analysis and clustering analysis. There were significant differences in the shape of the wing venation patterns between colonies of two sites (Wilk’s λ = 0.000006; p < 0.000001), which is probably due to the geographical distance between places of origin which impedes the gene flow between them. It indicates that inter and intrapopulation morphometric variability exists (p <0.000001) in M. scutellaris coming from two different biomes, revealing the existence of diversity in these populations, which is necessary for the conservation of this bee species. KEY WORDS: Geometric morphometrics, diversity, conservation. 1 Insecta Research Group, Centro de Ciências Agrárias, Ambientais e Biológicas, UFRB, C. Postal: 118, 44380-000, Cruz das Almas-BA. Brazil. E-mail: cristovamalves@yahoo.com.br; calfredo@ufrb. edu.br 2 Escola Superior de Agricultura “Luiz de Queiroz”, Universidade de São Paulo, 13418-900, Piracicaba- SP. Brazil. E-mail; lorenunes1@yahoo.com.br 3 Escola de Artes Ciências e Humanidades, Universidade de São Paulo, 03828-000, São Paulo-SP. Brazil. E-mail: tfrancoy@usp.br 108 Sociobiolog y Vol. 59, No. 1, 2012 INTRODUCTION The bees of the genus Melipona Illiger, 1806 occur throughout the Neo- tropical region, and this genus that has the largest number of species in the subtribe Meliponina (Silveira et al. 2002). Among the species of Melipona occurring in Bahia State, in Northeast Brazil, M. scutellaris Latreille 1811, M. quadrifasciata Lepeletier 1836 and M. mandacaia Smith 1863, have an emphasis as honey producers (Carvalho et al. 2003). Despite its economic and ecological importance, there are few studies about the current state of the populations conservation of these species. According to Nunes et al. (2007), this kind of information is important in the preparation of management plans for the preservation of M. scutellaris, since human actions, such as deforestation, use of insecticides and inadequate management of colonies, result in a decrease in populations of stingless bees and the consequential loss of biodiversity. In this context, the use of morphometric techniques have been used as a tool for morphological, taxonomic, phenotypic and geographic studies among and within bee populations, as in the works carried out by Crewe et al. (1994), Diniz-Filho & Malaspina (1996 ), Araujo et al. (2004), Mendes et al. (2007) and Francoy et al. (2011). Among the morphometric methods, geometric morphometrics provides powerful tools for the study of morphological variation among populations and therefore offers promising prospects for many problems such as the study of biodiversity patterns (Albert et al. 2001). In this technique, the first step is the definition of a structure and landmarks that are homologous among organisms under study (Monteiro & Reis 1999). In bees, wings have been used in morphometric analysis because they are quite flat structures and due to the ease of measuring the size and shape (Diniz-Filho & Bini 1994, Nunes et al. 2007). The objective of this work was to study the population divergence through the shape of wing venation of worker bees in M. scutellaris, originating from colonies living in two restricted areas, a semideciduous forest fragment in the Chapada Diamantina and an Atlantic forest fragment area in the coastal zone in the State of Bahia. 109 Lima, C.A. et al. — Population Divergence of M. Scutellaris in Brazil MATERIAL AND METHODS Populations of M. scutellaris were sampled between January and February 2009 in two restricted areas, one in an Atlantic forest fragment coastal zone in the city of Vera Cruz, Bahia, Itaparica Island (12°57’37’’S; 38°36’31’’W; Altitude: 13 m) and the other in a semideciduous forest fragment in the Chapada Diamantina region, in the municipality of Novo Mundo, Bahia (12°02’59”S and 40°29’43”W, altitude 604 m). These areas have a distance of 230 km between them. Traits of the right forewings of worker bees were used, with 10 workers of each colony and 10 colonies per area, totaling 200 subjects from 20 colonies. The colonies were named according to the place of origin, and to the Semi- deciduous Forest Fragment – colonies SFF 01 to 10, and the Atlantic Forest Fragment – colonies AFF 01 to 10. The wings were removed and placed be- tween two plates to capture images with the Motic 2.0 ML program using a digital camera coupled to a stereomicroscope, with 7.5 X magnification. The wings of each subject were identified and stored in Eppendorf tubes of 1.5 ml and deposited at the Núcleo de Estudos dos Insetos (Experimental Center for the Studies of Insects) (INSECTA), at the Centro de Ciências Agrárias, Ambientais e Biológicas (Center of Agricultural, Environmental and Biological Sciences) of the Federal University of Reconcavo of Bahia. Geometric Morphometric Analysis 15 landmarks on the forewing were defined and recorded using the tpsDig program, version 2.12 (Rohlf 2008a) as shown in Fig. 1. Fig. 1. Landmarks used in morphometric analysis of the forewing of Melipona scutellaris. 110 Sociobiolog y Vol. 59, No. 1, 2012 With the tpsRelw program, version 1.46 (Rohlf 2008b), using the archive of recorded images of landmarks, we extracted the x and y coordinates of each landmark of the aligned wings (alx and aly), the centroid size of the wings and the matrix W which consists of the deformed partial scores plus the estimate of the uniform components (Rohlf 1996). This matrix W can be used as a data matrix in order to perform various multivariate analyzes (Rohlf 2008b). Centroid size is the square root of the sum of squared distances from each landmark to the centroid (center of mass of the configuration), and this measure is, in the absence of allometry (correlation of size and shape), the only variable size that is not correlated with the variables of form (Rohlf 2008b). The partial deformation vectors are generated from the projection of the coordinates of each subject in the main deformation, in this way the partial distortions are the characteristics of each subject and are interpreted as the coordinates of each specimen (Ferreira 2009). The uniform components are calculated from all the landmarks in relation to the average shape and express information on a global scale (uniform) of the shape variation (Monteiro & Reis 1999). The visualization of the displacement of landmarks and the deformation grids was performed using the MorphoJ software, version 1.02g (Klingenberg 2011). The programs used for morphometric analysis are available on http://life. bio.sunysb.edu/morph. Statistical Analysis The coordinates of the landmarks of the aligned wings and the centroid size were used in principal component analysis (PCA). The data matrix W were used in multivariate analysis of variance (MANOVA) and cluster analysis by UPGMA (Unweighted Pair-Group Method using Arithmetic Average) performed with the STATISTICA 9.0 (Stat Soft 2010) program and discriminant analysis performed using the XLSTAT software, Version 2011.1.02. The cophenetic correlation coefficient was calculated by the sta- tistical Genes program (Cruz 2006). The centroid size was also assessed using an analysis of variance in the statistical STATISTICA 9.0 (Stat Soft 2010) program. 111 Lima, C.A. et al. — Population Divergence of M. Scutellaris in Brazil RESULTS AND DISCUSSION Multivariate Statistics Principal Component Analysis and Analysis of Multivariate Variance There were 27 principal components generated using the Cartesian coor- dinates of landmarks of the aligned wings and the centroid size of the wings. It took the first nine components to explain over 76.0% of the variation of the data (Table 1). The separation of the bees by area of origin was plotted graphically in a two- dimensional space formed by the scores of the first two principal components that explained 16.95% and 13.57% of the variation, respectively (Fig.2). Fig. 2. Dispersal of populations of Melipona scutellaris of two restricted areas in the state of Bahia with the principal component analysis. Points with the same format in the graph represent the worker bees belonging to the respective areas as shown in the legend. SFF = Semideciduous Forest Fragment; AFF = Atlantic Forest Fragment. 112 Sociobiolog y Vol. 59, No. 1, 2012 Variables aly4, aly9, aly11 and aly15 were the ones that contributed most to the first principal component, while for the second principal component, the variables that most influenced it were alx1, aly1, and alx2 aly13 (Table 2). The Multivariate Analysis of Variance (MANOVA) showed the existence of a highly significant difference in the shape of wings between colonies of different areas (Wilk’s λ = 0.000006; p < 0.000001). We also found differ- ence among colonies from the same location (SFF: Wilk’s λ = 0.000045; Table 1. Eigenvalues and total and cumulative variation of the principal components obtained from the Carte- sian coordinates of the aligned wing landmarks and the centroid size of the wings of Melipona scutellaris of two restricted areas in the state of Bahia. Principal Component Eigenvalue Total Variation (%) Cumulative Eigenvalue Cumulative Variation (%) 1 5.2537 16.9475 5.2537 16.95 2 4.2060 13.5677 9.4597 30.52 3 3.6318 11.7155 13.0915 42.23 4 2.6447 8.5313 15.7362 50.76 5 2.1901 7.0647 17.9263 57.83 6 2.0854 6.7270 20.0116 64.55 7 1.3303 4.2914 21.3420 68.85 8 1.2040 3.8840 22.5460 72.73 9 1.0800 3.4839 23.6260 76.21 10 0.9275 2.9920 24.5535 79.20 11 0.7939 2.5610 25.3474 81.77 12 0.6834 2.2044 26.0308 83.97 13 0.5996 1.9340 26.6304 85.90 14 0.5601 1.8069 27.1905 87.71 15 0.5082 1.6395 27.6987 89.35 16 0.4417 1.4249 28.1405 90.78 17 0.3966 1.2795 28.5371 92.06 18 0.3549 1.1448 28.8920 93.20 19 0.3355 1.0824 29.2275 94.28 20 0.2971 0.9584 29.5246 95.24 21 0.2918 0.9411 29.8164 96.18 22 0.2362 0.7620 30.0526 96.94 23 0.2296 0.7405 30.2822 97.68 24 0.2270 0.7322 30.5091 98.42 25 0.2107 0.6797 30.7198 99.10 26 0.1537 0.4959 30.8736 99.59 27 0.1264 0.4079 31.0000 100.00 113 Lima, C.A. et al. — Population Divergence of M. Scutellaris in Brazil p < 0.000001 and AFF: Wilk’s λ = 0.000168; p < 0.000001), demonstrating the existence of morphometric variability between and within the areas of origin of these colonies. All 26 variables used of the matrix W used in MANOVA contributed significantly (α = 0.05%) to the separation of the colony. In a similar context, Francoy et al. (2009) found that the obtained partial deformation of the wings of both males and the worker bees have contributed significantly to the separation of five species of stingless bees. It indicates that the variables obtained in the geometric morphometric analysis can be used to assess the morphometric variability between and within species of bees. The change of the overall shape shows the greatest variations in the landmarks 04, 11 and 05 (Fig. 3). Since landmarks 04 and 11 are located in the extreme points of the wing, it indicates a variation in the total length of this structure among the sampled bees. The size variation of the wings from the two populations was confirmed whit the significant differences found in the variance analysis of the centroid size (F = 70.3579; P <0.001). The bees originating from the semideciduous forest fragment have the centroid size greater than the bees from the Atlantic forest fragment, indicating that the wings of workers of M. scutellaris from the Semideciduous Forest Table 2. Contribution of variables to the princip a l c omp onents o bta ine d from cartesian coordinates of the aligned wing landmarks and the centroid size of Melipona scutellaris. alx and aly = x and y coordinates of each landmark of the aligned wings. Variable Principal Component CP01 CP02 Centroid Size -0.550453 -0.243003 alx1 -0.042075 0.776222 aly1 -0.306264 0.641001 alx2 -0.188666 -0.606984 aly2 0.202185 -0.121721 alx3 -0.177111 -0.518128 aly3 -0.250889 -0.136957 alx4 0.656952 0.434852 aly4 -0.787486 0.351908 alx5 -0.630051 -0.532158 aly5 -0.199672 -0.377015 alx6 0.343194 -0.241252 aly6 0.165304 0.175745 alx7 0.236879 0.066574 aly7 0.276928 -0.080970 alx8 0.207306 -0.229745 aly8 0.524310 -0.306929 alx9 0.194469 -0.113525 aly9 0.754889 -0.058321 alx10 0.148413 -0.162990 aly10 0.533427 -0.180559 alx11 -0.687989 0.309534 aly11 0.542941 0.173048 alx12 0.381650 0.135767 aly12 0.202381 -0.396282 alx13 -0.425072 -0.473367 aly13 -0.273201 -0.676573 alx14 0.002385 -0.037379 aly14 -0.014791 -0.282355 alx15 -0.073616 0.384956 aly15 -0.685634 0.435803 114 Sociobiolog y Vol. 59, No. 1, 2012 Fig. 3. Deformation grid showing the change in the overall shape of the wings of Melipona scutellaris of two restricted areas in the State of Bahia. Fig. 4. Wing centroid size of Melipona scutellaris of two restricted areas in the State of Bahia. SFF = Semideciduous Forest Fragment; AFF = Atlantic Forest Fragment. 115 Lima, C.A. et al. — Population Divergence of M. Scutellaris in Brazil Fragment are larger (Fig. 4). Araújo et al. (2004) in studies with six species of stingless bees, among which there were three species of Melipona, revealed that the generalized wing size is strongly correlated to the flight distance, which suggests that the stingless bees occupy a maximum effective space proportional to body size, especially with wing dimensions, which might constitute strong constraints on local populations restricted to forest fragments. Discriminant Analysis In the discriminant analysis performed using the matrix W, it can be seen that 91.0% of the bees were correctly classified within their respective colo- nies (Table 3). The Semideciduous Forest Fragment colonies stood out with five colonies that had correct classification rates of 100.0% (colonies SFF 03, 04, 06, 07 and 10). As for the Atlantic Forest Fragment area, there were four colonies with 100.0% of correct classification (colonies AFF 01, 03, 04 and 07). The colonies that had the lowest rates of correct classification (70.0%) were the colonies SFF 01 and AFF 10. Mendes et al. (2007), noted that in stingless bees, geometric morphometrics can be used as a technique for the discriminant analysis even in small populations. By cross-validation (Table 3), the average of correct identification of individuals to its correct colony of origin was 72.0%, with the colonies of the AFF 01 and 03 and SFF 06 standing out, which had 100.0% of subjects correctly identified. The discriminant analysis performed for the classification of subjects be- tween areas revealed that in the Atlantic Forest Fragment a correct classifica- tion rate of 100.0% was obtained, while the Semideciduous Forest Fragment rate was 99.0%, and the overall average was 99.5%. In cross-validation tests the overall average was 97.5%, and for the Semideciduous Forest Fragment the correct identification of its subjects was 98.0%, while for the Atlantic Forest Fragment the rate was 97.0%. Despite the absence of a record of discriminant analysis in populations of M. scutellaris, rates of correct classification and identification obtained in this study are considered high when compared to the works of Tofilski (2008) and Francoy et al. (2008), which had rates of 84.9 and 97.8% of correct identifi- cation by cross-validation for discrimination of subspecies of Apis mellifera with geometric morphometric variables, but also with those of Mendes et al. 116 Sociobiolog y Vol. 59, No. 1, 2012 (2007) and Francoy et al. (2009) in studies of species of stingless bees that reached from 77.6 to 100.0% of ratings and from 69.79 to 100.0% of correct identifications of subjects in respective groups. Comparing the rates of correct classification and identification considering the areas and also the colonies as a separate group, the highest correct rates found in the discrimination per area can be verified, since these higher rates of classifications and identifications are due to the greater variability in shape that are occurring between the areas of origins of the colonies than among the colonies themselves. Table 3. Classification and correct identification of subjects within their respective colonies and areas of origin by discriminant analysis and cross validation, respectively. SFF = Semideciduous Forest Fragment; AFF = Atlantic Forest Fragment. Colony Correct Classification (%) Correct Identification (%) SFF 01 70.00 60.00 SFF 02 90.00 70.00 SFF 03 100.00 80.00 SFF 04 100.00 40.00 SFF 05 90.00 80.00 SFF 06 100.00 100.00 SFF 07 100.00 80.00 SFF 08 80.00 50.00 SFF 09 90.00 90.00 SFF 10 100.00 80.00 AFF 01 100.00 100.00 AFF 02 90.00 70.00 AFF 03 100.00 100.00 AFF 04 100.00 90.00 AFF 05 90.00 80.00 AFF 06 80.00 70.00 AFF 07 80.00 40.00 AFF 08 100.00 80.00 AFF 09 90.00 40.00 AFF 10 70.00 40.00 Average 91.00 72.00 Area Correct Classification (%) Correct Identification (%) Semideciduous Forest Fragment 99.00 98.00 Atlantic Forest Fragment 100.00 97.00 Average 99.50 97.50 117 Lima, C.A. et al. — Population Divergence of M. Scutellaris in Brazil Cluster Analysis The Mahalanobis D2 distances between colonies were calculated for the cluster analysis (Table 4). All distances were statistically significant (α = 0.05) indicating the existence of significant differences among the colonies. Based on the Mahalanobis distance the clustering of the colonies was performed by the method of UPGMA. In the dendrogram obtained (Fig. 5) it can be observed the formation of nine groups consisting of: I – colonies SFF 01 and 02; II – colonies SFF 04, 08 and 09; III – colonies SFF 10; IV – colonies SFF 03 and 05; V – colonies SFF 06 and 07; VI - Colony AFF 01; VII – colony AFF 02; VIII - Colony AFF 03, and IX – colonies AFF 04, 05, 06, 07, 08, 09 and 10. The cophenetic correlation coefficient (CCC) for the dendrogram was 0.80 (p <0.01). The CCC shows a good fit between the graphical represen- tation of the original distances and its original matrix according to Bussab Fig. 5. Dendrogram generated by UPGMA through the Mahalanobis distances between colonies of the Semideciduous Forest Fragment (SFF) and the Atlantic Forest Fragment (AFF). Roman numbers represent the groups formed. 118 Sociobiolog y Vol. 59, No. 1, 2012 Ta bl e 4 . D is si m ila ri ty b et w ee n co lo ni es o f M el ip on a s cu te lla ri s c om in g f ro m a Se m id ec id uo us F or es t F ra gm en t a nd an A tla nt ic F or es t F ra gm en t i n th e S ta te of B ah ia . A t t he to p of th e ta bl e ar e th e M ah al an ob is d is ta nc es D 2 a nd a t t he b ot to m o f t he ta bl e is th e st at is tic al si gn ifi ca nc e (p ). SF F = S em id ec id uo us Fo re st F ra gm en t; A FF = A tla nt ic F or es t F ra gm en t. C ol on y SF F0 1 SF F0 2 SF F0 3 SF F0 4 SF F0 5 SF F0 6 SF F0 7 SF F0 8 SF F0 9 SF F1 0 A FF 01 A FF 02 A FF 03 A FF 04 A FF 05 A FF 06 A FF 07 A FF 08 A FF 09 A FF 10 SF F0 1 9. 46 48 .8 7 24 .6 8 61 .3 6 45 .3 8 51 .9 9 33 .0 3 29 .0 1 28 .2 4 87 .5 4 83 .9 2 10 5. 16 57 .1 2 49 .3 5 83 .5 5 83 .5 6 58 .8 3 56 .5 1 66 .0 1 SF F0 2 0. 05 38 .9 7 23 .6 1 48 .6 0 41 .1 2 45 .8 2 26 .7 4 34 .0 4 30 .0 7 97 .9 3 89 .9 6 10 8. 78 60 .3 5 57 .3 6 85 .4 4 94 .2 3 68 .2 7 64 .6 3 76 .4 0 SF F0 3 0. 00 0. 00 27 .6 3 18 .3 1 38 .3 5 33 .2 2 25 .0 4 27 .6 4 35 .1 2 86 .1 4 63 .0 4 57 .4 0 54 .0 7 62 .9 8 50 .1 6 73 .1 3 55 .3 5 42 .9 4 43 .8 3 SF F0 4 0. 00 0. 00 0. 00 39 .1 3 39 .7 0 45 .3 6 13 .6 5 28 .0 3 32 .9 9 77 .9 5 56 .7 4 78 .6 2 38 .8 0 33 .9 8 58 .0 6 68 .6 0 45 .2 5 38 .5 4 41 .8 8 SF F0 5 0. 00 0. 00 0. 00 0. 00 45 .8 9 23 .2 6 38 .8 0 43 .2 3 51 .4 9 11 1. 66 67 .1 2 67 .3 9 63 .4 5 73 .4 7 68 .6 9 93 .7 9 58 .6 9 53 .5 0 66 .1 6 SF F0 6 0. 00 0. 00 0. 00 0. 00 0. 00 27 .0 2 32 .5 3 30 .1 8 54 .2 1 92 .2 8 71 .0 3 75 .9 3 58 .7 8 56 .1 9 53 .7 8 73 .0 0 67 .2 6 45 .7 0 58 .2 0 SF F0 7 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 30 .9 1 32 .4 2 43 .8 2 10 2. 32 75 .6 8 89 .3 4 62 .0 6 76 .9 5 70 .4 6 85 .4 5 66 .3 5 55 .5 4 76 .4 7 SF F0 8 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 23 .7 1 38 .7 1 87 .0 3 61 .5 0 83 .5 9 42 .2 6 50 .2 1 66 .6 3 73 .5 1 59 .9 7 45 .5 1 55 .5 4 SF F0 9 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 37 .0 4 63 .9 8 59 .8 7 78 .3 3 50 .6 6 43 .7 0 53 .9 5 67 .9 0 63 .5 1 35 .0 2 43 .8 5 SF F1 0 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 99 .1 5 84 .0 3 10 3. 57 54 .0 4 72 .8 5 84 .9 4 87 .9 1 64 .2 0 57 .3 5 75 .7 6 A FF 01 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 38 .6 0 56 .2 0 40 .1 9 32 .3 4 53 .7 0 49 .3 3 51 .0 7 36 .5 6 37 .8 8 A FF 02 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 30 .2 7 28 .1 0 34 .9 2 47 .7 9 51 .4 6 32 .8 0 27 .7 5 43 .5 9 A FF 03 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 39 .7 7 50 .6 6 32 .4 4 42 .7 7 33 .9 3 32 .3 7 32 .2 8 A FF 04 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 19 .8 1 33 .1 8 26 .1 4 21 .1 1 19 .5 6 32 .0 1 A FF 05 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 28 .7 8 26 .9 7 22 .5 3 17 .3 0 20 .5 7 A FF 06 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 12 .4 6 26 .3 1 18 .4 6 15 .8 1 A FF 07 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 20 .6 3 19 .8 9 22 .6 6 A FF 08 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 22 .8 1 21 .1 9 A FF 09 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 17 .0 9 A FF 10 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 0. 00 119 Lima, C.A. et al. — Population Divergence of M. Scutellaris in Brazil et al. 1990, thus enabling the realization of inferences by means of visual assessment of the figure. In this dendrogram the formation of groups composed of a variable number of colonies stands out, with the existence of clusters with up to seven colo- nies (Group IX) and also colonies that are in isolation (Groups III and VI). This group formation reveals the similarity of the shape of wings among the colonies, and colonies grouped into a single set of wings are more similar to each other, while colonies that were isolated have wings more divergent in the populations studied. It was also noted that there is not any clustering of colonies from different origins at any time, which also shows the divergence between these colonies with respect to the location from which they were collected. According to Francoy & Imperatriz-Fonseca (2010), the variables generated in geometric morphometrics are independent of the size of structures, which eliminate part of the environmental influence that affects the development of organisms, especially insects. Thus, the shape variation found among the colonies, mainly coming from different areas may be due to genetic differences between these colonies. The subject’s phenotype (e.g. morpholog y, in this case the shape of the wing ) is the result of their interaction with genetic and environmental factors and when the effect of these environmental factors is eliminated there is a greater effect of genetics on the phenotype of the subject. As the distance of the stingless bee swarming is short (Kerr et al. 1996), the geographical distance between these environments can be considered as a factor that prevents the exchange of genetic material between the colonies, as there is no reference that there has been any handling or transport of the colonies from one area to another. Gonçalves (2010), when studying variability in Frieseomelitta varies by geometric morphometrics and molecular analysis of wings, suggested that the genetic diversity found in the small population under study is the result of gene flow via males, which have greater ability to fly and thus may originate from further distant regions and may have fertilized queens of the isolated population, thus introducing new alleles in that population. This hypothesis can also be adopted as the source of the diversity found among colonies of the same geographical origin, because near the locations where the bees were collected from there is a forest area with colonies in their natural state. 120 Sociobiolog y Vol. 59, No. 1, 2012 The wing shape variations found in this work show that the use of geomet- ric morphometrics of wings is a simple, practical, inexpensive and important tool in assessing diversity within populations of bees. CONCLUSION The geometric morphometric analysis of wings indicates that inter and intrapopulation morphometric variability exists in Melipona scutellaris com- ing from a semideciduous forest fragment and an Atlantic forest fragment in the State of Bahia, revealing the existence of diversity in these populations, which is necessary for the conservation of this bee species. ACKNOWLEDGMENTS The authors thank Conselho Nacional de Desenvolvimento Científico e Tec- nológico (CNPq) (Proc. 303237/2010-4), Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) (AUX-PE-PNPD-1598/2008) for the scholarships and financial support and Dr. Wyratan da Silva Santos for his critiques and suggestions. REFERENCES Addinsoft Sarl 2011. XLSTAT versão 2011.1.02. para Windows. Alibert P., B. Moureau, J.L. Dommergues & B. David 2001. Differentiation at a microgeographical scale within two species of ground beetle, Carabus auronitens and C. nemoralis (Coleoptera: Carabidae): a geometrical morphometric approach. Zoologica Scripta 30: 299-311. Araújo, E.D., M. Costa, J. Chaud-Netto & H.G. Fowler 2004. Body size and flight distance in stingless bees (Hymenoptera: Meliponini): Interference of flight range and possible ecological implications. Brazilian Journal of Biolog y 64: 563-368. Bussab, W.O., E.S. Miazaki & D.F. Andrade 1990. Introdução à análise de agrupamentos. São Paulo: IME/USP. 105p. Carvalho, C.A.L., R.M.O. Alves & B.A. Souza 2003. Criação de abelhas sem ferrão: aspectos práticos. Cruz das Almas: Universidade Federal da Bahia/SEAGRI. 42 p. (Série Meliponicultura – 01). Crewe, R.M., H.R. Hepburn & R.F.A. Moritz 1994. Morphometric analysis of 2 southern African races of honey bee. Apidologie 25: 61-70. Cruz, C.D. 2006. Programa Genes - Biometria. 1. ed. Viçosa, MG: Editora UFV. 382 p. Diniz-Filho, J.A.F. & L.M. Bini 1994. Space-free correlation between morphometric and climatic data: a multivariate analysis of Africanized honey bees (Apis mellifera L.) in Brazil. Global Ecolog y and Biogeography. Letters 4: 195-202. 121 Lima, C.A. et al. — Population Divergence of M. Scutellaris in Brazil Diniz-Filho, J.A.F. & O. Malaspina 1996. Geographic variation of Africanized honey bees (Apis mellifera) in Brazil - multivariate morphometrics and racial admixture. Brazilian Journal of Genetics 19: 217-224. Ferreira, V.S. 2009. Análise morfométrica e genética da variabilidade em populações de Centris (Centris) aenea Lepeletier, 1841 (Hymenoptera, Apidae, Centridini), uma abelha polinizadora de fruteiras. Dissertação de Mestrado, Universidade Estadual de Santa Cruz, Ilhéus–BA. 78 p. Francoy, T.M., D. Wittmann, M. Drauschke, S. Müller, V. Steinhage, M.A.F. Bezerra-Laure, D. De Jong & L.S. Gonçalves 2008. Identification of Africanized honey bees through wing morphometrics: two fast and efficient procedures. Apidologie, 39: 488- 494. Francoy, T.M., R.A.O. Silva, P. Nunes-Silva, C. Menezes & V.L. Imperatriz-Fonseca 2009. Gender Identification of five genera of stingless bees (Apidae, Meliponini) based on wing morpholog y. Genetics and Molecular Research 8: 207-214. Francoy ,T.M. & V.L. Imperatriz-Fonseca 2010. A morfometria geométrica de asas e a identificação automática de espécies de abelhas. Oecologia Australis 14: 317-321. Francoy, T.M., M.L. Grassi, V.L. Imperatriz-Fonseca, W.J. May-Itzá & J.J.G. Quezada-Euán 2011. Geometric morphometrics of the wing as a tool for assigning genetic lineages and geographic origin to Melipona beecheii (Hymenoptera: Meliponini). Apidologie, 42: 499- 507. Gonçalves, P.H.P. 2010. Análise da variabilidade de uma pequena população de Frieseomelitta varia (Hymenoptera, Apidae, Meliponini), por meio de análise do DNA mitocondrial, microddatélites e morfometria geométrica das asas. Dissertação de Mestrado, Universidade de São Paulo, São Paulo–SP. 140 p. Kerr, W.E., G.A. Carvalho & V.A. Nascimento. 1996. Abelha Uruçu: Biologia, Manejo e Conservação. Belo Horizonte- MG, Ancagaú. 143p. Klingenberg , C. P. 2011. MorphoJ: an integrated software package for geometric morphometrics. Molecular Ecolog y Resources 11: 353–357. Mendes, M.F.M., T.M. Francoy, P. Nunes-Silva, C. Menezes & V.L. Imperatriz-Fonseca 2007. Intra-populational variability of Nannotrigona testaceicornis Lepeletier 1836 (Hymenoptera, Meliponini) using relative warp analysis. Bioscience Journal 23: 147- 152. Monteiro, L.R. & S.F. Reis 1999. Princípios de morfometria geométrica. Ribeirão Preto: Holos. 188p. Nunes, L.A., M.F.F.C. Pinto, P. Carneiro, D.G. Pereira & A.M. Waldschmidt 2007. Divergência Genética em Melipona scutellaris Latreille (Hymenoptera: Apidae) com base em Caracteres Morfológicos. Bioscience Journal 23: 1-9. Rohlf, F.J. 1996. Morphometric spaces, shape components and the effect of linera transformations, p. 117-130. In Marcus LF, Corti M, Loy A, Naylor G, Slice D (Ed.) Advances in Morphometrics. New York: Plenum Publishind Corp. 587p. Rohlf, F.J. 2008a. Program tpsDIG2 for Windows version 2.12. Departament of Ecolog y and Evolution, State University of New York, Stony Book. (http://life.bio.sunysb.edu/morph/index.html, acessado em 21/outubro/2008). 122 Sociobiolog y Vol. 59, No. 1, 2012 Rohlf, F.J. 2008b. Program tpsRelw version 1.46. Department of Ecolog y and Evolution, State University of New York, Stony Book. (http://life.bio.sunysb.edu/morph/index. html, acessado em 21/outubro/2008). Stat Soft 2010. Statistica for windows: versão 9.0. Tulsa. Tofilski, A. 2008. Using geometric morphometrics and standard morphometry to discriminate three honeybee subspecies. Apidologie 39:558-563.