Genetic divergence between Melipona quadrifasciata Lepeletier (Hymenoptera, Apidae) populations

RESEARCH ARTICLE

Genetic divergence between Melipona quadrifasciata Lepeletier (Hymenoptera, Apidae) populations

Mara Garcia Tavares; Nathalia Teixeira Pietrani; Maxwell de Castro Durvale; Helder Canto Resende; Lucio Antonio de Oliveira Campos

Departamento de Biologia Geral, Universidade Federal de Viçosa, Viçosa, MG, Brazil

^{Send correspondence to}

ABSTRACT

Melipona quadrifasciata is a stingless bee widely found throughout the Brazilian territory, with two recognized subspecies, M. quadrifasciata anthidioides, that exhibits interrupted metasomal stripes, and M. quadrifasciata quadrifasciata, with continuous metasomal stripes. This study aimed to estimate the genetic variability of these subspecies. For this purpose, 127 colonies from 15 Brazilian localities were analyzed, using nine species-specific microsatellite primers. At these loci, the number of alleles ranged from three to 15 (mean: 7.2), and the observed heterozygosity (H_o) ranged from 0.03-0.21, while the expected heterozygosity (H_e) ranged from 0.23-0.47. The genetic distances among populations ranged from 0.03-0.45. The F_ST multilocus value (0.23) indicated that the populations sampled were structured, and the clustering analysis showed the formation of two subgroups and two more distant populations. The first group contained the subspecies M. quadrifasciata quadrifasciata, and the other, the subspecies M. quadrifasciata anthidioides and the two M. quadrifasciata populations with continuous metasomal stripes from northern Minas Gerais. These results confirmed that the yellow metasomal stripes alone are not a good means for correctly identifying the different subspecies of M. quadrifasciata.

Keywords: genetic differentiation, microsatellites, population genetics, stingless bees.

Introduction

Despite their ecological importance as pollinators of native trees in Brazil (Kerr et al., 2001; Imperatriz-Fonseca et al., 2006), stingless bee populations have decreased in the Neotropics, similarly to other bee populations, due to habitat degradation and fragmentation, agricultural intensification, indiscriminate use of agrochemicals and pesticides, and introduction/spread of exotic species (Freitas et al., 2009).

Melipona quadrifasciata, in particular, is a stingless bee species with a wide geographic distribution in Brazil (Silveira et al., 2002; revised in Camargo and Pedro, 2007). Few molecular studies, however, have been conducted on this species to analyze its genetic variability and structuring of populations (Waldschmidt et al., 2000, 2002; Moreto and Arias, 2005; Souza et al., 2008; Batalha-Filho et al., 2009, 2010; Nascimento et al., 2010).

In general, these studies have provided evidence that it is possible to recognize molecular differences between the two M. quadrifasciata subspecies, confirming the usefulness of the yellow metasomal stripes, which are continuous in M. quadrifasciata quadrifasciata and interrupted in M. quadrifasciata anthidioides, to distinguish them (Schwarz, 1932).

However, according to Moure and Kerr (1950), although M. quadrifasciata quadrifasciata is found in the colder regions located from Rio Grande do Sul to southern São Paulo, and M. quadrifasciata anthidioides occurs from northeastern São Paulo along eastern Brazil up to Paraíba, a hybridization zone exists in the states of São Paulo and Minas Gerais, with specimens presenting an intermediary pattern of yellow metasomal stripes.

More recently, M. quadrifasciata individuals from Januária and Urucuia (northern Minas Gerais) have been observed exhibiting continuous metasomal stripes identical to those of M. quadrifasciata quadrifasciata instead of the interrupted pattern characteristic to M. quadrifasciata anthidioides, as was expected based on their geographical localization (review in Batalha-Filho et al., 2009). The specimens from Januária, however, did not exhibit the RAPD marker present in individuals of the subspecies M. quadrifasciata quadrifasciata, being therefore genetically more similar to M. quadrifasciata anthidioides (Waldschmidt et al., 2000).

Another population bearing a continuous pattern of metasomal stripes was later on described in Sergipe and northeastern Bahia (Batalha-Filho et al., 2009). RFLP analysis and sequencing of the mtDNA COI gene again demonstrated that this population grouped together with populations from northern Minas Gerais and populations of

M. quadrifasciata anthidioides (Batalha-Filho et al., 2010). These data demonstrate that a distinction between M. quadrifasciata anthidioides and M. quadrifasciata quadrifasciata based exclusively on morphological characteristics is not reliable. Thus, the main objective of this study was to investigate the genetic variability of some M. quadrifasciata populations from different locations in Brazil, in order to improve the knowledge about their taxonomic status.

Materials and Methods

Bee samples

Melipona quadrifasciata adult workers from 127 colonies were collected in eight Brazilian states. M. quadrifasciata quadrifasciata samples were collected in the states of Rio Grande do Sul (Porto Alegre), Santa Catarina (Içara) and Paraná (Curitiba), while M. quadrifasciata anthidioides samples were obtained in Rio de Janeiro (Miguel Pereira), Minas Gerais (Cristiano Otoni, Piranga, Caeté, Rio Vermelho and Poté), Espírito Santo (Domingos Martins) and Bahia (Caitité). We also included samples from Urucuia/Januária and São Cristovão in the analyses, representing M. quadrifasciata populations with continuous metasomal stripes from northern Minas Gerais and Sergipe, respectively. It is important to emphasize that these populations were located completely outside the distribution area of M. quadrifasciata quadrifasciata which typically exhibits continuous yellow metasomal stripes.

Samples from Ourolândia (Bahia) presenting a hybrid pattern of metasomal stripes were also analyzed.

Samples from Cristiano Otoni, Caeté, Urucuia and Januária were collected from feral colonies, directly from nature. Samples from Piranga, Rio Vermelho, Poté, São Cristovão, Ourolândia and Caitité were provided by beekeepers that maintain a small number of colonies in the original trunks or in rational boxes. These colonies were collected in their propriety wood or in the neighborhood, approximately 500-1000 m from their present-day localization. Samples from Domingos Martins, Miguel Pereira, Içara, Porto Alegre and Curitiba were obtained from colonies kept in local meliponaries. These colonies, however, were also obtained from the respective regions and did not represent translocated colonies.

The geographic distribution and the number of colonies sampled in each locality are shown in Table 1. Samples were stored in absolute ethanol at -80 ºC until DNA extraction.

DNA extraction and amplification

Genomic DNA from 127 individuals (one per colony) was extracted according to the protocol of Waldschmidt et al. (1997) and amplified using nine microsatellite primers specifically designed for M. quadrifasciata in our laboratory (Table 2).

To perform the amplification, the forward (F) sequence of each primer was modified at the 5' end with the addition of a sequence of the M13 universal primer (5'-TTTTCCCAGTCACGA-3') (Schuelke, 2000). Additionally, the M13 universal primer was labeled with a fluorescent dye (FAM-6-carboxy-fluorescein or HEX hexacloro-6-carboxy-fluorescein).

The amplification reaction solutions consisted of 10 ng DNA, 200 µM of each dNTP, 1.0 µM of each primer (F, R and M13), 1.5 mM MgCl₂, and 0.5 unit of Platinum Taq polymerase (Invitrogen) in 10x buffer. Amplifications were performed by an initial denaturation step of 94 ºC for 3 min, followed by 10 cycles of 94 ºC for 15 s, with a specific annealing temperature for each primer for 20 s and 72 ºC for 30 s. Subsequently, 25 cycles at 89 ºC for 15 s, 53 ºC for 20 s, 72 ºC for 30 s, and a final extension at 72 ºC for 30 min were performed. The polymerase chain reaction (PCR) products were identified using a MegaBACE automated capillary sequencer. Allele sizes were scored against the size standard ET 550R (GE Healthcare Life Sciences, Pittsburgh, PA, USA), and the peaks were analyzed using the MegaBACE Fragment Profiler program.

Statistical analysis

The Popgene version 1.32 program (Yeah et al., 1999) was used to calculate the following estimators of genetic diversity: allele frequency, expected (H_e) and observed (H_o) heterozygosity, genetic distance between populations (Nei, 1978), and F-statistics (Wright, 1978). For these statistical analyses, colonies were grouped according to locality, and each of these groups was considered to be a population.

Grouping analysis based on the UPGMA algorithm and the analysis of molecular variance (AMOVA; Excoffier et al., 1992) were carried out using the Genes program (Cruz, 2012). The cophenetic correlation coefficient that measures to which extent the clustering result corresponds to the original dissimilarity matrix was used to test the efficiency of the clustering analysis for the data (Rohlf and Fisher, 1968), and AMOVA was calculated considering two hierarchic levels and the groups formed in the clustering analysis.

Results

The nine microsatellite loci analyzed presented 65 alleles, ranging from three (Mquad24) to 15 (Mquad7) alleles (mean: 7.2 alleles/locus). All loci showed deviations from the Hardy-Weinberg equilibrium (p < 0.05) at one or more sites, but no single locus showed deviations in all cases. Table S1 shows the allele frequencies observed for these loci, considering all the analyzed populations.

The most frequent alleles of the different loci were present in both populations of M. quadrifasciata quadrifasciata and M. quadrifasciata anthidioides, as well as in the populations with continuous metasomal stripes from northern Minas Gerais and Sergipe and in the hybrid morphotype from Bahia. In general, however, the populations exhibited different allelic constitutions, with some fixed and other private alleles (Table S2).

In the different analyzed populations, the observed heterozygosity (Ho) ranged from 0.03-0.21, while the expected heterozygosity (He) ranged from 0.23-0.47.

The genetic distances among populations ranged from 0.05 (Piranga and Rio Vermelho, Cristiano Otoni and Caeté, Cristiano Otoni and Urucuia) to 0.45 (Januária and São Cristovão) (Table 3).

Clustering analysis demonstrated that populations of both subspecies, although clustering together in a single major group, can be divided into two subgroups, one containing the populations of the subspecies M. quadrifasciata quadrifasciata and the other the populations of the subspecies M. quadrifasciata anthidioides and the two M. quadrifasciata populations with continuous metasomal stripes from northern Minas Gerais (Figure 1). Only populations from São Cristovão and Ourolândia, which showed higher genetic distance values than the other populations (Table 3), were isolated in this analysis. The cophenetic correlation of this clustering was high (0.87), further showing a good adjustment between the genetic distances and the graphically represented distances.

The FST multilocus value obtained when populations of both subspecies were considered together (FST: 0.23) also showed that they were structured, and the AMOVA test demonstrated that the interpopulation variation (52.5%) was higher than the intrapopulation variation (47.5%) of the subgroups obtained in the clustering analysis. Differences in the allele frequencies among populations as well as their private alleles probably contributed to these results.

Discussion

Our results showed a relatively high level of multiallelism in the nine microsatellite loci analyzed, which is characteristic when species-specific microsatellite primers are used (Peters et al., 1998; Paxton et al., 2003; Brito et al., 2009; Kapheim et al., 2009; Lopes et al., 2009a,b; Oliveira et al., 2009; Francisco et al., 2011).

This high number of alleles, however, was evidenced only when all populations were considered together. When each population was considered separately, only few alleles were found in each of the loci analyzed, always with one of them at a higher frequency. Consequently, a low frequency of heterozygotes was detected, which may explain the differences between the expected and observed heterozygosity values found. Additionally, despite having sampled 13 to 18 colonies in some localities, this sample size can be considered relatively small, which might also explain this difference, as well as the observed deviations from the Hardy-Weinberg equilibrium.

Notwithstanding, the lower level of genetic variability detected in the M. quadrifasciata populations (that generally mate once) compared to the values reported for M. bicolor (Ho: 0.40, Peters et al., 1998) and Apis mellifera (Ho: 0.35-0.59, De La Rúa et al., 2003), for example, may be due to biological differences between the species, such as the facultative polygyny of M. bicolor (Cepeda, 2006; Velthuis et al., 2006), the multiple mating of honeybee queens (Estoup et al., 1994) or the limited flight dispersion distance of stingless bees (Araújo et al., 2004).

Although M. quadrifasciata presents a widespread geographic distribution, occupying distinct biomes such as the Atlantic rainforest, gallery forests (seasonal forests) and arborous caatinga, environment degradation may also be contributing to their reduced genetic variability. This hypothesis should be considered, as it was already observed in

M. rufiventris and M. mondury populations (Tavares et al., 2007) that environmental degradation may isolate populations in small forest fragments, leading to inbreeding and local extinction. An indication of this local extinction is provided by Batalha-Filho et al. (2009) who, in contrast to Kerr (1951), found no M. quadrifasciata specimens in Pernambuco and Paraiba, despite intensive searching.

It is worthy of note that many of the private alleles detected in the present study were found in samples from Urucuia/Januária (northern Minas Gerais). Similarly, allele A of locus Mquad24 was found only in M. quadrifasciata population from São Cristovão (Sergipe) and in the hybrid population from Ourolândia (Bahia). As mentioned above, although the samples from Urucuia, Januária and São Cristovão contain M. quadrifasciata populations with continuous metasomal stripes, they were located completely outside the distribution area of M. quadrifasciata quadrifasciata.

Therefore, the presence of these private alleles is evidence of restricted genetic flow between these and the other populations of M. quadrifasciata anthidioides and could explain the genetic differentiation detected for the populations of São Cristovão and Ourolândia. Genetic differentiation among some M. quadrifasciata populations has also been detected using ISSR markers (Nascimento et al., 2010), PCR-RFLP assays of mtDNA (Batalha-Filho et al., 2009), and mtDNA COI gene sequencing (Batalha-Filho et al., 2010). According to these authors, geographic barriers separating close neighbors of this species can explain this differentiation.

The clustering analyses showed that M. quadrifasciata populations from Curitiba, Içara and Porto Alegre, which correspond to the subspecies M. quadrifasciata quadrifasciata, form a separate subgroup from the M. quadrifasciata populations of Urucuia and Januária, which exhibit a continuous metasomal stripe pattern similar to this subspecies, but are genetically more similar to M. quadrifasciata anthidioides populations.

The two other morphologically atypical M. quadrifasciata populations (São Cristovão and Ourolândia), however, were found to be isolated in the clustering analyses performed in this study. These same populations were clustered with M. quadrifasciata anthidioides in the phylogeographic study of Batalha-Filho et al. (2010), although the presence of exclusive haplotypes was detected in populations with continuous metasomal stripes from Sergipe. This result may represent differences between the two molecular markers used (mitochondrial vs. nuclear). In the present case, as microsatellites represent nuclear markers with a high mutation rate, the presence of two or three exclusive alleles was sufficient to differentiate populations. Further studies are necessary to elucidate the taxonomic status of these populations.

Overall, the results obtained in the present study corroborate those of previous studies demonstrating genetic differences between the two M. quadrifasciata subspecies (Waldschmidt et al., 2000; Moreto and Arias, 2005; Souza et al., 2008; Batalha-Filho et al., 2009, 2010). Additionally, as already observed by Batalha-Filho et al. (2009, 2010), our results also indicated that the M. quadrifasciata populations from northeastern Bahia and Sergipe are genetically different from the M. quadrifasciata quadrifasciata samples, even though they show the characteristic continuous metasomal stripes pattern.

Acknowledgments

The authors are grateful to the Fundação de Amparo à Pesquisa do Estado de Minas Gerais (FAPEMIG) for the financial support granted. We also thank Dr. Leonardo Lopes Bhering and Dr. Cosme Damião Cruz from Universidade Federal de Viçosa for statistical assistance with the GENES program. We also wish to thank Dr. Fabrício Rodrigues dos Santos from Universidade Federal de Minas Gerais for the help with sequencing analyses.

Supplementary Material

The following material concerning this article is available online:

Table S1 -Allele frequencies of the microsatellite loci.

Table S2 -Distribution of the alleles present in the nine loci analyzed.

This material is available as part of the online article at http://www.scielo.br/gmb.

Received: May 2, 2012

Accepted: August 13, 2012.

Associate Editor: Klaus Hartfelder

License information: This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

References

• Araújo ED, Chaud-Neto J and Fowler HG (2004) Body size and flight distance in stingless bees (Hymenoptera, Meliponini): Inference of flight range and possible ecological implications. Braz J Biol 64:563-568.
• Batalha-Filho H, Melo GAR, Waldschmidt AM, Campos LAO and Fernandes-Salomão TM (2009) Geographic distribution and spatial differentiation in the color pattern of abdominal stripes of the Neotropical stingless bee Melipona quadrifasciata (Hymenoptera, Apidae). Zoologia 26:213-219.
• Batalha-Filho H, Waldschmidt AM, Campos LAO, Tavares MG and Fernandes-Salomão TM (2010) Phylogeography and historical demography of the Neotropical stingless bee Melipona quadrifasciata (Hymenoptera, Apidae): Incongruence between morphology and mitochondrial DNA. Apidologie 41:534-547.
• Brito RM, Francisco FO, Domingues-Yamada AMT, Gonçalves PHP, Pioker FC, Soares AEE and Arias MC (2009) Characterization of microsatellite loci of Tetragonisca angustula (Hymenoptera, Apidae, Meliponini). Conserv Genet Resour 1:183-187.
• Camargo JMF and Pedro SEM (2007) Meliponini Lepeletier, 1836. In: Moure JS, Urban D and Melo GAR (eds) Catalogue of Bees (Hymenoptera, Apoidea) in the Neotropical Region. Sociedade Brasileira de Entomologia, Curitiba, pp 272-578.
• Cepeda OI (2006) Division of labor during brood production in stingless bees with special reference to individual participation. Apidologie 37:175-190.
• Cruz CD (2012) Programa Genes: Aplicativo Computacional em Genética e Estatística. Versão Windows -2012. UFV, Viçosa.
• De La Rúa P, Galián J, Serrano J and Moritz RFA (2003) Genetic structure of Baleartic honeybee populations based on microsatellite polymorphism. Genet Sel Evol 35:339-350.
• Estoup A, Solignac M and Cornuet JM (1994) Precise assessment of the number of patrilines and of genetic relatedness in honeybee colonies. Proc R Soc Lond B 258:1-7.
• Excoffier L, Smouse PE and Quattro JM (1992) Analyses of molecular variance inferred from metric distances among DNA haplotypes: Application to human mitochondrial DNA restriction data. Genetics 131:479-491.
• Francisco FO, Brito RM, Santiago LR, Gonçalves PHP, Pioker FC, Domingues-Yamada AMT and Arias MC (2011) Isolation and characterization of 15 microsatellite loci in the stingless bee Plebeia remota (Apidae, Meliponini). Conserv Genet Resour 3:417-419.
• Freitas BM, Imperatriz-Fonseca VL, Medina LM, Kleinert AMP, Galetto L, Nates-Parra G and Queza-Eúan JJG (2009) Diversity, threats and conservation of native bees in the Neotropics. Apidologie 40:332-346.
• Imperatriz-Fonseca VL, Saraiva AM and De Jong D (2006) Bees as Pollinators in Brazil: Assessing the Status and Suggesting Best Practices. Holos, Ribeirão Preto, 111 pp.
• Kapheim KM, Pollinger JP, Wcislo WT and Wayne RK (2009) Characterization of 12 polymorphic microsatellite markers for a facultatively eusocial sweat bee (Megalopta genalis). Mol Ecol Res 9:1527-1529.
• Kerr WE (1951) Estudos sobre a genética de populações de Himenópteros em geral e dos Apíneos sociais em particular. Tese para livre docência. An Esc Super Agric Luiz de Queiroz 8:219-354.
• Kerr WE, Carvalho GA, Da Silva AC and Assis MG (2001) Aspectos pouco mencionados da biodiversidade amazônica. Parc Estrat 12:20-41.
• Lopes DM, Silva FO, Fernandes-Salomão TM, Campos LAO and Tavares MG (2009a) Microsatellite loci for the stingless bee Melipona rufiventris (Hymenoptera, Apidae). Mol Ecol Res 9:923-925.
• Lopes DM, Silva FO, Fernandes-Salomão TM, Campos LAO and Tavares MG (2009b) A scientific note on the characterization of microsatellite loci for Melipona mondury (Hymenoptera, Apidae). Apidologie 41:138-140.
• Moreto G and Arias MC (2005) Detection of mitochondrial DNA restriction site differences between the subspecies of Melipona quadrifasciata Lepeletier (Hymenoptera, Apidae, Meliponini). Neotrop Entomol 34:381-385.
• Moure JS and Kerr WE (1950) Sugestões para a modificação da sistemática do gênero Melipona (Hymen.-Apoidea). Dusenia 1:105-129.
• Nascimento MA, Batalha-Filho H, Waldschmidt AM, Tavares MG, Campos LAO and Fernandes-Salomão TM (2010) Variation and genetic structure of Melipona quadrifasciata Lepeletier (Hymenoptera, Apidae) populations based on ISSR pattern. Genet Mol Biol 33:394-397.
• Nei M (1978) Estimation of average heterozygosity and genetic distance from a small number of individuals. Genetics 89:583-590.
• Oliveira EJF, Freitas GS, Fonseca AS, Sousa ACB, Campos T, Assis AF, Souza AP, Contel EPB and Soares AEE (2009) Isolation and characterization of microsatellite markers from the stingless bee Nannotrigona testaceicornis Conserv Genet Resour 1:97-99.
• Paxton RJ, Arévalo E and Field J (2003) Microsatellite loci for the eusocial Lasioglossum malachurum and other sweat bees (Hymenoptera, Halictidae). Mol Ecol Notes 3:82-84.
• Peters JM, Queller DC, Imperatriz-Fonseca VL and Strassmann JE (1998) Microsatellite loci for stingless bees. Mol Ecol 7:784-787.
• Rohlf FJ and Fisher DL (1968) Test for hierarchical structure in random data sets. Syst Zool 17:407-412.
• Schuelke M (2000) An economic method for the fluorescence labeling of PCR fragments. Nat Biotechnol 18:233-234.
• Schwarz H (1932) The genus Melipona IV. The type genus of the Meliponidae or stingless bees. Bull Am Mus Nat Hist 63:231-460.
• Silveira FA, Melo GAR and Almeida EAB (2002) Abelhas Brasileiras: Sistemática e Identificação, Fundação Araucária, Belo Horizonte, 253 pp.
• Souza RO, Moretto G, Arias MC and Del Lama MA (2008) Differentiation of Melipona quadrifasciata L. (Hymenoptera, Apidae, Meliponini) subspecies using cytochrome b PCRRFLP patterns. Genet Mol Biol 30:445-450.
• Tavares MG, Dias LAD, Borges AA, Lopes DM, Busse AHP, Costa RG, Salomão TMF and Campos LAO (2007) Genetic divergence between populations of the stingless bee uruçu amarela (Melipona rufiventris group, Hymenoptera, Meliponini): Is there a new Melipona species in the Brazilian state of Minas Gerais? Genet Mol Biol 30:667-675.
• Velthuis HW, De Vries H and Imperatriz-Fonseca VL (2006) The polygyny of Melipona bicolor: Scramble competition among queens. Apidologie 37:222-239.
• Waldschmidt AM, Salomão TMF, Barros EG and Campos LAO (1997) Extraction of genomic DNA from Melipona quadrifasciata (Hymenoptera, Apidae, Meliponinae). Braz J Genet 20:421-423.
• Waldschmidt M, Barros EG and Campos LAO (2000) A molecular marker distinguishes the subspecies Melipona quadrifasciata quadrifasciata and Melipona quadrifasciata anthidioides (Hymenoptera, Apidae, Meliponini). Genet Mol Biol 23:609-611.
• Waldschmidt AM, Marco-Júnior P, Barros EG and Campos LAO (2002) Genetic analysis of Melipona quadrifasciata Lep. (Hymenoptera, Apidae, Meliponini) with RAPD markers. Braz J Biol 62:923-928.
• Wright S (1978) Evolution and the Genetics of Populations. Variability Within and Among Natural Populations. University of Chicago Press, Chicago, 590 pp.
• Yeah FC, Yang R and Boyle T (1999) Popgene ver. 1.32: Microsoft Windows-based freeware for population genetic analysis. Quick User Guide.

Publication Dates

• Publication in this collection
  28 Feb 2013
• Date of issue
  2013

History

• Received
  02 May 2012
• Accepted
  13 Aug 2012