Genetic structure of the gentle Africanized honey bee population (gAHB) in Puerto Rico
© Galindo-Cardona et al.; licensee BioMed Central Ltd. 2013
Received: 12 January 2013
Accepted: 1 August 2013
Published: 6 August 2013
The Africanized honey bee is one of the most spectacular invasions in the Americas. African bees escaped from apiaries in Brazil in 1956, spread over Americas and by 1994 they were reported in Puerto Rico. In contrast to other places, the oceanic island conditions in Puerto Rico may mean a single introduction and different dynamics of the resident European and new-coming Africanized bees.
To examine the genetic variation of honey bee feral populations and colonies from different locations in Puerto Rico, we used eight known polymorphic microsatellite loci.
In Puerto Rico, gAHB population does not show any genetic structure (Fst = 0.0783), and is best described as one honey bee population, product of hybridization of AHB and EHB. The genetic variability in this Africanized population was similar to that reported in studies from Texas. We observed that European private allele frequencies are high in all but one locus. This contrasts with mainland Africanized populations, where European allele frequencies are diminished. Two loci with European private alleles, one on Linkage Group 7, known to carry two known defensiveness Quantitative Trait Loci (QTLs), and the other on Linkage Group 1, known to carry three functionally studied genes and 11 candidate genes associated with Varroa resistance mechanisms were respectively, significantly greater or lower in European allele frequency than the other loci with European private alleles.
Genetic structure of Puerto Rico gAHB differs from mainland AHB populations, probably representing evolutionary processes on the island.
Humans and honey bees interacted over history with consequences for distribution of the honey bee, the provider of the first sweetener or honey for humans. Honey bees are naturally distributed throughout Africa, Europe, and the Middle East. In this broad range, they have evolved into 24 subspecies, which have been grouped into three , four  or five  distinct evolutionary branches. Whitfield and colleagues  showed using SNP data that honey bees left Africa three times, first to Asia, second to Europe, and third when the African bee Apis mellifera scutellata was introduced to America in 1956. Before the introduction of the African honey bees, a subset of European subspecies and other African genotypes were introduced into the New World, first by settlers and later by beekeepers [4–6]. Since their introduction to Brazil in 1956, African bees hybridized with different European honey bees in Americas [5–13]. This was an accidental experiment, where natural hybridization between races of honey bees occurred through time. The result was the Africanized honey bees (AHB), shown to be distinct from the original African race or the European Honey Bees (EHB) in its genetics and behavioral characteristics [14–16]. On the oceanic island of Puerto Rico, AHB have been shown to have a mosaic of traits such as European-like low defensiveness behavior and African-like resistance to ectoparasites e.g. to Varroa spp. [14–16]. We dubbed these Puerto Rican bees as gAHB because of their unusual gentle behavior. AHB most likely arrived in Puerto Rico in the first half of 1990s they were reported first in 1994; , and now all feral colonies sampled are Africanized (see [14–16], 112 samples, , and this study-72 samples).
Hybridization between AHB and EHB has been previously reported in Argentina, Brazil, Perú, and Mexico [5, 6, 19, 20]. The current thinking is that when AHB arrives at a location with a resident European population, there is a large amount of hybridization, followed by elimination of European genes over time, perhaps due to selective advantage of African genes under tropical conditions . The mosaic nature of observed characteristics of gAHB may be better understood with data on population genetics [14–16], potentially providing support for different processes discussed by , extending from hybridization to founder effect to selection.
Many studies on population genetics in the Americas have been conducted [5, 6, 19–22], but in Puerto Rico, the population genetic structure of A. mellifera is unclear. Natural hybridization between EHB and AHB has been occurring over the last 17 years, and currently, based on mtDNA and morphology, Puerto Rico has a uniform geographic distribution of gAHB, both feral and domestic [14, 15]. Previous studies on the genetic structure of honey bees in Puerto Rico relied on the use of morphology, mtDNA markers, and single locus nuclear RFLP markers [14, 15]. Other approaches used worldwide to explore the population structure and hybridization zones of honey bees include morphology [13, 23] and different molecular genetic markers to study genetic variation i.e. Restriction Fragment Length Polymorphisms and Short Tandem Repeat-microsatellites, SNPs; e.g. [3, 24, 25]. To further examine the genetic structure of the Puerto Rican honey bee population, we used eight commonly selected microsatellite markers for which historical data on other Africanized and European populations exist [3, 21]. These molecular markers could help detect potential source populations for the gAHB and determine any population substructure on the island. One advantage of microsatellites is their high variability that allowed us to use individual-based measures of relatedness to study gene flow and population boundaries. In addition, the microsatellite loci could be spread across the genome, and occur at identifiable specific locations . Here we chose only eight microsatellite loci shown to be highly informative [21, 26]. More loci could be added since there are many microsatellite loci developed and published for A. mellifera (e.g. 250 reported in ).
The study of genetic structure and variation of gAHB will allow us to test two alternate hypotheses to explain the mix of AHB and EHB traits in this population: Structure and Hybridization hypotheses. Testing of these hypotheses provides a first approach to detect isolation local and global selection:
The “Structure hypothesis” is presence of different composition of European or African alleles in different subpopulations across the island of Puerto Rico. This could be due to local selection or isolation due to geographic barriers (Puerto Rico is composed of one main island and two adjacent islands and exhibit six life zones with different climate and floral phenology). In the local selection scenario, alternate European or African traits that best respond to environmental pressure under different conditions would be selected in different regions across Puerto Rico. According to this hypothesis different levels of European genes could be maintained in populations of feral bees at geographically distinct parts of the island, due to the different environmental conditions. Puerto Rico has a topography that has a strong effect on the climate, with the wetter regions on the windward, northern side of the mountains, and drier climate in the leeward rain shadow . Moreover, Puerto Rico’s six different life zones [28, 29], could provide opportunity for ecologically adapted, isolated population subdivisions. We could therefore expect European traits, for instance, on higher cooler regions of the island. In the isolation scenario, populations with high European genes in isolated regions such as surrounding islands may be observed. Structure hypothesis then predicts genetic differences among delimited geographic areas, in other words, a genetic structure for honey bee subpopulations in Puerto Rico. This will also predict a departure from Hardy-Weinberg Equilibrium (HWE) for different marker loci, indicating population substructure due to genetic drift or due to local selection on associated loci.
Alternately, the “Hybridization hypothesis” based on previous studies, states that there is only one gAHB population with large introgression of European genes, as demonstrated in one nuclear locus [14–16]. According to the hybridization hypothesis, there is no clear differentiation of subpopulations of gAHB in Puerto Rico. This could be due to a brief and effective hybridization of AHB and EHB, with a restricted flow of African genes from outside, followed by a spread of the hybrid bee across the small island. Any selective advantage for alternate AHB or EHB alleles would be present early, and spread to the whole population under a global selection scenario. In absence of a selective advantage for AHB or EHB alleles, hybridization hypothesis predicts that there will be a single, unstructured bee population with similar introgression of European genes across the sampled molecular markers. In the absence of additional Africanized bees breeding and spreading into the newly occupied territory, the sampled African alleles in Puerto Rico should not be increasing over time, but remain similar to the initial hybrid population. Alternately, if certain AHB or EHB alleles were to confer an early and fast selective advantage, hybridization hypothesis predicts that there will be different levels of introgression of EHB genes, commensurate with selective advantage at different loci, but with no current deviation from HWE. Although few known markers would be adequate to test the immediate prediction of a uniform genetic structure as per the hybridization hypothesis, early and fast hybridization would make it very unlikely to detect global selection at a resolution of 8 microsatellite markers. Future greater density marker studies, such as 1000 s of SNPs identified in SNP arrays [3, 30] or RAD-TAG sequencing with potential to generate multiples of 10000 SNPs  is more likely to provide resolution needed to determine regions of genome responding to alternate global selection pressures.
Sample collection and DNA extraction
A total of 72 colony samples were collected from several apiaries and feral colonies from different places in Puerto Rico: Main island (55), Vieques (17). This is a large sample for population genetic studies using microsatellite markers with high levels of polymorphism . We also included 11samples from the typical European bee commercially available to Puerto Rico, the Italian bee or Apis mellifera ligustica, purchased from Hawaii (11). We include more samples from the Vieques island because it is separated from Puerto Rico by a mass of water and may provide an important test for the Structure hypothesis. In addition, in 1988 a program of producing European queen bees was established in Vieques. This discontinued program was established because since that date, Puerto Rico was quarantined against import of European bees from elsewhere (except Hawaii) to prevent introduction of diseases and parasites of the honey bee . However, in our studies we have not encountered EHB morphology or mtDNA on repeated sampling from Vieques [14–16] and this study.
Briefly, we used the protocol from  to identify mtDNA from Africanized (AHB) and Eastern European (A.m. ligustica or Italian or EHB) bees, as described by . Previous studies had demonstrated for Puerto Rico only maternal descent of these two types . We amplified a region of the rDNA of 964 bp. The products from PCR were later subjected to digestion with the restriction enzyme EcoR1. This enzyme digests only the EHB product to two fragments of nearly equal size (484 and 480 bp), resulting in a band that moves further on the gel in comparison to AHB band. The bands were visualized after running in a 1% agarose gel and staining with ethidium bromide.
Bees were stored in 95% ethanol and kept at 4°C until processed. To prevent oversampling of maternal alleles, only one bee per colony was subjected to genetic analysis. Genomic DNA from an individual honey bee thorax was extracted using DNeasy extraction kit from QIAGEN® with the animal tissue protocol.
List of the primers sequences of microsatellite loci and the fluorescent dyes
Loci, sample size (N), no. alleles (Na), no. effective alleles (Ne), no. alleles private to European (NaPE) bees in previous studies, observed heterozygosity (Ho), expected (He) and unbiased expected heterozygosity (UHe), average departure of genotype frequencies within populations ( F is ), correlation between gametes within an individual relative to the entire population ( F it ), and fixation index ( F st ) by population for codominant data Puerto Rico (PR) and Vieques (VIE)
The number of alleles, allele frequencies, and observed (Ho) and expected (He) heterozygosity of each microsatellite locus were estimated using GENEPOP 4.0 . GENEPOP 4.0 was also used to test for linkage disequilibrium and HWE. A regression analysis was used to relate the allelic frequency of European population from Puerto Rico with the reference European population. The reference European population was taken from [21, 22]. Our results were compared with Texas results from  that were divided in two groups: first (1991-1996) and second (1997-2001) group, representing early and later stages of Africanization. Estimates of population structure and gene flow were assessed using F-statistics, also with GENEPOP 4.0. Unbiased gene diversity (Hd), and allelic richness (Rs) per locus and population sample were computed using FSTAT version 2.9.3 package . In FSTAT, the sample size for allelic richness estimation is fixed as the smallest number of individuals typed for a locus in a sample, which in this study was 72 (Locus ED1). Population structure was also analyzed using a Bayesian model-based clustering method provided in STRUCTURE 2.3.3 [41, 42]. This program infers structure of populations using allele frequencies of unlinked markers (microsatellites). The parameter set was programmed for independent alleles and individuals to have a mixed ancestry. The program was asked to place individuals from all populations into one and two groups, and in three runs the analyses consisted of 10000, 100000, and 1000000 burn-in replicates, and run lengths of 10000, 100000, and 1000000 replicates. The best estimate of K or the number of populations was determined by looking at the values of log Pr(X/K) and the value of α [41, 42]. In structure analysis we did not use Culebra samples, because the low sample number dissuaded us from doing the comparison. Single and multiloci Fst values were estimated using the weighed analysis of variance method by  to measure genetic differentiation with the GENETIX 4.02 program . We analyzed distant locations that appeared to be genealogically close on an unrooted neighbor joining tree based on genetic distance, constructed using the Phylip 3.57 program .
Loci, number of alleles per locus (Na), unbiased gene diversity (Hd), allelic richness (Rs) to Texas (TX1991/1996, TX1997/2001, *Pinto 2003) and new alleles to Puerto Rico (PR2007)
Population genetic structure
In Puerto Rico, we found that 18 percent of European alleles are shared with the Texas population, in eight microsatellite loci. Exact tests for linkage disequilibrium at 5% level for each pair of loci did not provide any significant P-value after Bonferroni corrections.
For these eight markers, we identified an average of 12.75 alleles by loci with heterozygosities in concordance with the expectation based on the HWE. We identified 12 (11.8 %) European private alleles in six of the eight microsatellite loci sampled to be present in the gAHB population (Table 3). The presence of all previously reported European private alleles in the gAHB population suggests that individuals have a mixed genetic origin. The locus A107 has six of these private European alleles, and loci A113, A79, A14 and A35 have one private European allele each. The locus A88 has two alleles exclusively European , and that were not present in the population studied by Pinto and colleagues . In Puerto Rico, these two private European alleles are present at a similar frequency to that reported in Hawaii EHB European population .
We found a few alleles only in Vieques that were not found in Puerto Rico, but their frequency was low. Large numbers of markers studied in other AHB populations are needed to determine the African ancestry of Puerto Rican gAHB population. Nevertheless, a genetic distance analysis revealed gAHB to be closer to Texas population than to the Brazil population (results not shown), indicating a more northern source for AHB in Puerto Rico than Brazil, and consistent with speculated Texas origin .
The population genetics analysis indicates the presence of a uniform genetic structure or a single gAHB population in Puerto Rico, consistent with the hybridization hypothesis. Interestingly, two marker loci show introgression levels for EHB alleles significantly different than for other loci tested. This is in agreement with a global selection scenario for AHB or EHB alleles for different traits.
One issue with highly polymorphic marker loci has been the perceived need to obtain a large sample . This probably is important especially if the “population” is actually composed of multiple subdivisions. Because we found a single population in Puerto Rico (including Vieques), the sample size of 72 is well above the recommended 20 to 25 for a single population, in relevant reference for the recent analysis. In addition, we ran multiple tests, based on allele frequencies, and also based on individual genotypes, giving the same result of a single population without subdivisions (see results and Figure 3). In a further attempt to test the accuracy of these analyses, we also performed the individual assignment to the Vieques and main island “populations” using GeneClass2  and also a spatial PCA using Adegenet on R platform. All four analyses gave similar results (analyses not shown), confirming the inference of a single continuous population on Puerto Rico.
The continuous uniform population could be the result of the large range of movement of the bees in comparison to the size of the island. For instance, reproductive drone congregation areas, DCAs are found at about 1.5 miles radius around bee colonies, and less than 10 such DCAs can span the island from North to South . Two other factors that may contribute to the homogenization of bee population in Puerto Rico are the formation of swarms and adoption of feral swarms by beekeepers. A swarm is a colony division, which involves migration of a part of the individuals. This movement of feral colonies occurs throughout the year, but their frequency is high at the beginning of the mating season (May to October, ). Second, an activity that is common for beekeepers around the island is to adopt feral bee colonies to populate their commercial hives. Regularly replacing queens of commercial colonies is considered a good beekeeping practice. However, in Puerto Rico only few beekeepers introduce new queens to their colonies. These queens are mainly from Hawaii (EHB), and rarely from other places in the United States. These factors would contribute to observed genotypic homogenization throughout Puerto Rico.
The variability of the microsatellite loci in the honey bee population in Puerto Rico is similar to that found in the EHB populations for a detected number of alleles and the values of gene diversity , (Table 3). Pinto and colleagues , found that their results were divided in two temporal groups, first (1991-1996) and second (1997-2001). These groups illustrate a pattern of genotypic differentiation, which shows that the first group was more similar to EHB and the second group more similar to AHB. When these two groups were compared with Puerto Rico, we found that Puerto Rican honey bee populations are more similar to the first group (European-like), lacking the temporal component indicating change towards more Africanized features as reported in Texas . In Puerto Rico, we used a paired t-test for eight loci used both in this study, and reported previously  for heterozygosity comparison. We found that the stable Africanized population in Texas (years 1997-2001) had higher heterozygosity, than the Puerto Rico population even when the allelic richness was comparable (positive yet non-significant t). These results may mean that after the original introduction of AHB, there has been no or little influx of new AHB to Puerto Rico, probably because of its island nature. The increased Africanized genes and heterozygosity in AHB population in the mainland (i.e. Texas) could be due to the many spatial and temporal points of entry available to the AHB and selection for Africanized traits in the mainland. In Puerto Rico, it is safe to assume that these arrivals of new AHB are very low. This restricted flow of AHB to the island could explain the genetic difference observed with the late Texas population (Table 3). The early impoverished EHB genes due to loss of feral bees to Varroa mites both in the mainland and in Puerto Rico, could contribute to the similarity . The Puerto Rican Africanized bees would be more similar to early AHB population in Texas in absence of constant AHB gene flow. However, under hybridization without selection, all markers across the genome should show similar introgression levels. Only hybridization followed by alternate alleles selected for different loci could explain the pattern observed in this study.
Introgression is a genetic invasion of a local genome by a foreign genome [48–52]. Although one idea would be that reduced defensiveness is the result of interbreeding throughout the island [14–16], it is intriguing to find that the population maintains Africanized features such as small size, high parasite defenses, and Africanized mtDNA . Therefore selection on particular traits, in particular regions of the genome globally across the island population is the favored hypothesis. In fact, there are two exceptional loci that may indicate presence of selection for European or African genes near the loci. First, there may have been selection for European genes near the locus A107. The locus A107 showed a higher frequency of European alleles in Puerto Rico, when compared with pure EHB populations (Figure 4). Moreover, this locus has 50 percent of all the European private alleles we found in this study. This A107 locus is on the same linkage group (LG7) that also contains the two Quantitative Trait Loci for aggressive behavior (Sting 1 and 2 ).
The deviations in allele frequencies for specific loci may be tentative evidence for global selection for a mosaic of African-like and European-like alleles. This strong global selection may have followed soon after hybridization process since the markers and candidate genes are relatively distant considering the high rate of recombination observed in honey bees . If this is true, for gAHB we predict that LG7 will be mostly European, associated with lower aggressive behavior, and LG1 would be mostly African, associated with greater Varroa mite defenses. This hypothesis of fast evolution or selection was suggested previously . An alternate hypothesis that could address the linkage on the face of high recombination rate is presence of other linked genes associated with the trait, especially for Varroa resistance on LG1, near the A79 marker. Interestingly, in the time this paper was under review, a new publication identified 11 more candidate genes on LG1 in a linkage study for Varroa resistance , supporting the prediction of the more linked genes hypothesis (Figure 5). The two hypotheses, fast evolution and more linked genes, are not mutually exclusive, in fact they could be synergistic in explaining the markers associated with trait genes even with high rate of recombination. Both hypotheses can be tested in future studies by use of high density molecular markers to be produced by next-generation sequencing methods.
The long-term stability of distribution of EHB-AHB allele frequencies suggests a stable equilibrium with reduced effects of genetic drift . The high genetic similarity of honey bee populations in Puerto Rico to initial phase of AHB Texas population, with large EHB contribution is interesting considering the AHB has now been in Puerto Rico for over 18 years. This again suggests an environment different than the mainland subtropical conditions [15, 16]. Because behavioral and now population genetics data indicate formation of a distinct gAHB population in Puerto Rico, the limited entry of European queens to the island, and potential rapid evolution and mosaic selection scenarios , make it important to examine origin and evolution of gAHB using high density molecular markers.
We observed that European private allele frequencies are high in all but one locus. This contrasts with mainland Africanized populations, where European allele frequencies are greatly diminished. Two loci with European private alleles, one linked to two known defensiveness QTLs, and the other to 11 known genes associated with Varroa resistance mechanisms were respectively, significantly greater or lower in European allele frequency than the other loci with European private alleles. Taken together, this evidence supports retention of a mosaic of African and European alleles for certain genes, mirroring the mosaic of African and European traits observed in bees from this gAHB population.
We consider this work is important for population geneticists and evolutionary biologists, because Africanization of honey bee in Puerto Rico is relatively recent, and combination of alleles from European and African bees presented in this study has implications for evolutionary processes on islands. Sociogenomics researchers may be interested in this population with increasing use of functional genomics and next-gen sequencing to study mechanisms underlying described variation in behavior of honey bees.
We thank three anonymous reviewers and Manuel Ramirez, María Alejandra Palacio and Gene E. Robinson for their helpful comments and suggestions that improved the manuscript. We also thank Gene E. Robinson for coining the term gAHB. Special thanks to Ximena Velez, Elsie Rivera, Dania Rodriguez, Karla Maldonado, Silvia Planas, Ivania Cerón, Shakira Quiñones, Nancy Torres, Jay Lou Torres, Mariela, Ailleen for their help in the laboratory bench, analyses and valuable comments. Biol 6996 Behavioral Plasticity seminar students also improved the presentation of this manuscript with their critical reading and comments. Thanks to Caro, Dani and Nico for their patience and impatience. Thanks to beekeepers of Puerto Rico for providing the samples, and the director Manuel Diaz and personnel of the Estación Experimental Agrícola de Gurabo for help with the bee research facility or “Casa Amarilla”. Thanks to the Sequencing and Genotyping Facility in generating and analyzing the data (NCRR-AABRE Grant #P20 RR16470, NIH-SCORE Grant #S06GM08102, University of Puerto Rico Biology Department, NSF-CREST Grant #0206200). Thanks to Humberto Ortiz and High Performance Computing Facility at UPR for additional Structure analyses of data (Puerto Rico INBRE Grant P20 RR-016470 from the National Center for Research Resources - NIH, and the Institute for Functional Nanomaterials award 0701525 from the NSF-EPSCoR). This work has been funded by the USDA NIFA award 2009-05291 to TG and support from the AMP program.
- Ruttner F, Tassencourt L, Louveaux J: Biometrical-statistical analysis of the geographic variability of Apis mellifera L. Apidologie. 1978, 9: 363-381. 10.1051/apido:19780408.View ArticleGoogle Scholar
- Ruttner F: Biogeography and taxonomy of honey bees. 1988, Berlin: SpringerView ArticleGoogle Scholar
- Whitfield CW, Behura SK, Berlocher SH, Clark AG, Johnston JS, Sheppard WS, Smith DR, Suarez AV, Weaver D, Tsutsui ND: Thrice out of Africa: ancient and recent expansions of the honey bee, Apis mellifera. Science. 2006, 314: 642-646. 10.1126/science.1132772.View ArticlePubMedGoogle Scholar
- Morse RA, Burgett DM, Ambrose JP, Conner WE, Fell RD: Early introductions of African bees into Europe and the New World. Bee World. 1973, 54: 57-60.View ArticleGoogle Scholar
- Sheppard WS, Rinderer TE, Mazzolli JA, Stelzer JA, Shimanuki H: Gene flow between African- and European-derived honey bee populations in Argentina. Nature. 1991, 349: 782-784. 10.1038/349782a0.View ArticleGoogle Scholar
- Sheppard WS, Soares AEE, DeJong D, Shimanuki H: Hybrid status of honey bee populations near the historic origin of Africanization in Brazil. Apidologie. 1991, 22: 643-652. 10.1051/apido:19910607.View ArticleGoogle Scholar
- Michener CD: The Brazilian bee problem. Ann Rev Entomol. 1975, 20: 399-416. 10.1146/annurev.en.20.010175.002151.View ArticleGoogle Scholar
- Taylor OR: African Bees: potential impact in the United States. Bull Entomol Soc Am. 1985, 31: 14-24.Google Scholar
- Rinderer TE: Africanized bees: the Africanization process and potential range in the United States. Bull Entomol Soc Am. 1986, 32: 222-227.Google Scholar
- Hall HG, Muralidharan K: Evidence from mitochondrial DNA that African honey bees spread as continuous maternal lineages. Nature. 1989, 339: 211-213. 10.1038/339211a0.View ArticlePubMedGoogle Scholar
- Smith DR, Taylor OR, Brown WM: Neotropical Africanized honey bees have African mitochondrial DNA. Nature. 1989, 339: 213-215. 10.1038/339213a0.View ArticlePubMedGoogle Scholar
- Hall HG: Parental analysis of introgressive hybridization between African and European honeybees using nuclear DNA RFLPs. Genetics. 1990, 125: 611-621.PubMed CentralPubMedGoogle Scholar
- Rinderer TE, Stelzer JA, Oldroyd BP, Buco SM, Rubink WL: Hybridization between European and Africanized Honey bees in the neotropical Yucatan peninsula. Science. 1991, 253: 309-311. 10.1126/science.253.5017.309.View ArticlePubMedGoogle Scholar
- Rivera-Marchand B: PhD thesis. Africanized Honey Bees in Puerto Rico. 2007, Biology Department, Rio Piedras Station: University of Puerto RicoGoogle Scholar
- Rivera-Marchand B, Giray T, Guzman-Novoa E: The cost of defense in social insects: insights from the honey bee. Entomol Expt Appl. 2008, 129: 1-10. 10.1111/j.1570-7458.2008.00747.x.View ArticleGoogle Scholar
- Rivera-Marchand B, Oskay D, Giray T: Gentle Africanized bees on an oceanic island. Evol Appl. 2012, 10.1111/j.1752-4571.2012.00252.x.Google Scholar
- Cox B: Ahb In Puerto-Rico. Am Bee J. 1994, 134: 668-668.Google Scholar
- Galindo-Cardona A: PhD thesis. Male behavior and hybridization of Africanized and European bees. 2010, Biology Department, Rio Piedras Station: University of Puerto RicoGoogle Scholar
- Quezada-Euan JJG, May-Itza WD: Partial seasonal isolation of African and European-derived Apis mellifera (Hymenoptera: Apidae) drones at congregation areas from subtropical Mexico. Ann Entomol Soc Am. 2001, 94: 540-544. 10.1603/0013-8746(2001)094[0540:PSIOAA]2.0.CO;2.View ArticleGoogle Scholar
- Perez-Castro EE, May-Itza WD, Quezada-Euan JJG: Thirty years after: a survey on the distribution and expansion of Africanized honey bees (Apis mellifera) in Peru. J Apic Res. 2002, 41: 69-73.Google Scholar
- Pinto MA, Rubink WL, Patton JC, Coulson RN, Johnston JS: Africanization in the United States: Replacement of Feral European Honey Bees (Apis mellifera L.) by an African Hybrid Swarm. Genetics. 2005, 170: 1653-1665. 10.1534/genetics.104.035030.PubMed CentralView ArticlePubMedGoogle Scholar
- Kraus FB, Franck P, Vandame R: Asymmetric introgression of African genes in honeybee populations (Apis mellifera L.) in Central Mexico. Heredity. 2007, 99: 233-240. 10.1038/sj.hdy.6800988.View ArticlePubMedGoogle Scholar
- Rinderer TE, Rubink BS, Daly WL, Stelzer HV, Riggio JA, Baptista RM: Morphometric identification of Africanized and European honey bees using large reference populations. Apidologie. 1993, 24: 569-585. 10.1051/apido:19930605.View ArticleGoogle Scholar
- Nielsen DI, Ebert PR, Page RE, Hunt GJ, Guzman-Novoa E: Improved polymerase chain reaction-based mitochondrial genotype assay for identification of the Africanized honey bee (Hymenoptera: Apidae). Ann Entomol Soc Am. 2000, 93: 1-6. 10.1603/0013-8746(2000)093[0001:IPCRBM]2.0.CO;2.View ArticleGoogle Scholar
- Suazo A, Hall HG: Nuclear DNA PCR-RFLPs that distinguish African and European honey bee groups of subspecies. I: Comparison of long PCR and standard PCR to screen for polymorphisms. Biochem Genet. 2002, 40: 225-239. 10.1023/A:1019830800282.View ArticlePubMedGoogle Scholar
- Solignac M, Vautrin D, Baudry E, Mougel F, Loiseau A, Cornuet J-M: A microsatellite-based linkage map of the honeybee, Apis mellifera L. Genetics. 2004, 167: 253-262. 10.1534/genetics.167.1.253.PubMed CentralView ArticlePubMedGoogle Scholar
- Gould WA, Alarcon C, Fevold B, Jimenez ME, Martinuzzi S, Potts G, Quiñones M, Solorzano M, Ventosa E: The Puerto Rico Gap Analysis Project. Volume 1: Land cover, vertebrate species distributions, and land stewardship. 2008, U.S. Department of Agriculture, Forest Service, International Institute of Tropical Forestry. U.S.A.: Gen, Tech. Rep. IITF-GTR-39. Rio Piedras, PRGoogle Scholar
- Holdridge LR: Life zone ecology. 1967, San Jose, Costa Rica Tropical: Science CenterGoogle Scholar
- Ewel JJ, Whitmore JL: The ecological life zones of Puerto Rico and the U.S. Virgin Islands. 1973, Institute of Tropical Forestry, Rio Piedras, PR: USDA Forest Service Research Paper No. ITF-18Google Scholar
- Zayed A, Whitfield CW: A genome-wide signature of positive selection in ancient and recent invasive expansions of the honey bee Apis mellifera. PNAS. 2008, 105: 3421-3426. 10.1073/pnas.0800107105.PubMed CentralView ArticlePubMedGoogle Scholar
- Hohenlohe PA, Bassham S, Etter PD, Stiffler N, Johnson EA, Cresko WA: Population Genomics of Parallel Adaptation in Threespine Stickleback using Sequenced RAD Tags. PLoS Genet. 2010, 6 (2): e1000862-10.1371/journal.pgen.1000862.PubMed CentralView ArticlePubMedGoogle Scholar
- Daly C, Helmer EH, Quinones M: Mapping the climate of Puerto Rico, Vieques and Culebra. Int J Climatol. 2003, 23: 1359-1381. 10.1002/joc.937.View ArticleGoogle Scholar
- Hale ML, Burg TM, Steeves TE: Sampling for microsatellite-based population genetic studies: 25 to 30 individuals per population is enough to accurately estimate allele frequencies. PLoS ONE. 2012, 7 (9): e45170-10.1371/journal.pone.0045170.PubMed CentralView ArticlePubMedGoogle Scholar
- Bauzá J: Para evitar la introducción a la isla de Vieques de las enfermedades y plagas peligrosas de las abejas (Apis mellifera). 1988, Departamento de Agricultura de Puerto Rico Cuarentena doméstica No.5Google Scholar
- Estoup A, Solignac M, Cornuet JM: Precise assessment of the number of patrilines and of genetic relatedness in honeybee colonies. Proc R Soc Lond B Biol Sci. 1994, 258: 1-7. 10.1098/rspb.1994.0133.View ArticleGoogle Scholar
- Estoup A, Garnery L, Solignac M, Cornuet JM: Microsatellite variation in Honey Bee (Apis mellifera L.) populations: hierarchical genetic structure and test of the infinite allele and stepwise mutation models. Genetics. 1995, 140: 679-695.PubMed CentralPubMedGoogle Scholar
- Rowe DJ, Rinderer TE, Stelzer JA, Oldroyd BP, Crozier RH: Seven polymorphic microsatellite loci in honeybees (Apis mellifera). Insectes Soc. 1997, 44: 85-93. 10.1007/s000400050032.View ArticleGoogle Scholar
- Schuelke M: An economic method for the fluorescent labeling of PCR fragments. Nat Biotech. 2000, 18: 233-234. 10.1038/72708.View ArticleGoogle Scholar
- Raymond M, Rousset F: Population genetics software for exact test and ecumenism. J Hered. 1995, 86: 248-250.Google Scholar
- Goudet J: FSTAT, a program to estimate and test gene diversities and fixation indices. 2001, version 2.9.3) [http://www2.unil.ch/popgen/softwares/fstat.htmGoogle Scholar
- Pritchard JK, Stephens M, Donnelly P: Inference of population structure using multilocus genotype data. Genetics. 2000, 155: 945-959.PubMed CentralPubMedGoogle Scholar
- Falush D, Stephens M, Pritchard JK: Inference of population structure using multilocus genotype data: Linked loci and correlated allele frequencies. Genetics. 2003, 164: 1567-1587.PubMed CentralPubMedGoogle Scholar
- Weir BS, Cockerham CC: Estimating F-statistics for the analysis of population structure. Evolution. 1984, 38: 1358-1370. 10.2307/2408641.View ArticleGoogle Scholar
- Belkhir KP, Borsa P, Chikhi L, Raufaste N, Bonhomme F: Genetix 4.04, logiciel sous windows TM pour la genetic des populations. 2002, Laboratoire Génome, Populations, Interactions, Montpellier: Université de Montpellier, FranceGoogle Scholar
- Felsenstein J: PHYLIP, phylogeny inference package, Version 3.57. 1993, Seattle: University of WashingtonGoogle Scholar
- Piry S, Alapetite A, Cornuet JM, Paetkau D, Baudouin L, Estoup A: GENECLASS2: A Software for Genetic Assignment and First-Generation Migrant Detection. J Hered. 2004, 95 (6): 536-539. 10.1093/jhered/esh074.View ArticlePubMedGoogle Scholar
- Galindo-Cardona A, Monmany AC, Moreno-Jackson R, Rivera-Rivera C, Huertas-Dones C, Caicedo-Quiroga L, Giray T: Landscape analyses of drone congregation areas of the honey bee Apis mellifera. J Insect Sci. 2012, 12: 104-View ArticleGoogle Scholar
- Barton NH, Hewitt GM: Adaptation, speciation and hybrid zones. Nature. 1989, 341: 497-503. 10.1038/341497a0.View ArticlePubMedGoogle Scholar
- Arnold ML: Natural hybridization as an evolutionary process. Annu Rev Ecol Syst. 1992, 23: 237-261. 10.1146/annurev.es.23.110192.001321.View ArticleGoogle Scholar
- Arnold ML: Natural hybridization and evolution. 1997, New York: Oxford University PressGoogle Scholar
- Rhymer JM, Simberloff D: Extinction by Hybridization and Introgression. Annu Rev Ecol Syst. 1996, 27: 83-10.1146/annurev.ecolsys.27.1.83.View ArticleGoogle Scholar
- Mallet J: Hybridization as an invasion of the genome. Trends Ecol Evol. 2005, 20: 229-10.1016/j.tree.2005.02.010.View ArticlePubMedGoogle Scholar
- Hunt GJ, Amdam GV, Schlipalius D, Emore C, Sardesai N, Williams CE, Rueppell O, Guzmán-Novoa E, Arechavaleta-Velasco M, Chandra S, Fondrk MK, Beye M, Page RE: Behavioral genomics of honeybee foraging and nest defense. Naturwissenschaften. 2007, 94: 247-267. 10.1007/s00114-006-0183-1.PubMed CentralView ArticlePubMedGoogle Scholar
- Navajas M, Migeon A, Alaux C, Martin-Magniette ML, Robinson GE, Evans JD, Cros-Arteil C, Crauser D, Le Cont Y: Differential gene expression of the honey bee Apis mellifera associated with Varroa destructor infection. BMC Genomics. 2008, 9: 301-10.1186/1471-2164-9-301.PubMed CentralView ArticlePubMedGoogle Scholar
- Anderson DL, Trueman JWH: Varroa jacobsoni (Acari: Varroidae) is more than one species. Exp Appl Acarol. 2000, 24: 165-189. 10.1023/A:1006456720416.View ArticlePubMedGoogle Scholar
- Beye M, Gattermeier I, Hasselmann M, Gempe T, Schioett M, Baines JF, Schlipalius D, Mougel F, Emore C, Rueppell O, Sirviö A, Guzmán-Novoa E, Hunt G, Solignac M, Page RE: Exceptionally high levels of recombination across the honey bee genome. Genome Res. 2006, 16: 1339-1344. 10.1101/gr.5680406.PubMed CentralView ArticlePubMedGoogle Scholar
- Tsuruda JM, Harris JW, Bourgeois L, Danka RG, Hunt GJ: High-Resolution Linkage Analyses to Identify Genes That Influence Varroa Sensitive Hygiene Behavior in Honey Bees. PLoS ONE. 2012, 7 (11): e48276-10.1371/journal.pone.0048276.PubMed CentralView ArticlePubMedGoogle Scholar
- Endler JA: Natural selection in the wild. 1986, New Jersey: Princeton University PressGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.