Mitochondrial DNA variability in the Canary Islands honeybees (Apis mellifera L.)

The mitochondrial DNA (mtDNA) of individuals from 79 colonies of Apis mellifera from five Canary Islands was studied using the Dra I test based on the restriction of PCR products of the tRNAleu–COII intergenic region. Five haplotypes of the African (A) lineage and one of the west European (C) lineage were found. The haplotypes A14 and A15 are described for the first time. These haplotypes have a new P sequence named P₁. The wide distribution and high frequency of haplotype A15 suggest that it is characteristic of the Canarian Archipelago. Sources of haplotype variability of honeybee mtDNA in the Canary Islands (waves of colonization from Africa, queen importations, habitat diversification) are discussed.     

Genetic structure and distinctness of Apis mellifera L. populations from the Canary Islands

The genetic structure of Apis mellifera populations from the Canary Islands has been assessed by mitochondrial (restriction fragment length polymorphisms of the intergenic transfer RNAleu-COII region) and nuclear (microsatellites) studies. These populations show a low level of genetic variation in terms of average number of alleles and degree of heterozygosity. Significant differences in the distribution of alleles were found in both data sets, confirming the genetic differentiation among some of the islands but not within them. Two mitochondrial haplotypes characteristic of the Canary Islands are found at high frequencies, although populations are introgressed by imported honeybees of eastern European C lineage. This introgression is rather high on Tenerife and El Hierro and low on Gran Canaria and La Gomera, whereas on La Palma it has not been recorded. The finding of microsatellite alleles characteristic of the eastern European lineage corroborates the genetic introgression. Phylogenetic analyses indicate that the Canarian honeybees are differentiated from other lineages and provide genetic evidence of their African origin.     

Turkish honeybees: genetic variation and evidence for a fourth lineage of Apis mellifera mtDNA

The mtDNA of bees from 84 colonies of Turkish honeybees (Apis mellifera) was surveyed for variation at four diagnostic restriction sites and the sequence of a noncoding intergenic region. These colonies came from 16 locations, ranging from European Turkey and the western Mediterranean coast to the Caucasus Mountains along the Georgian border, the eastern Lake Van region, and the extreme south. Combined restriction site and sequence data revealed four haplotypes. Three haplotypes belonged to the eastern Mediterranean mtDNA lineage. The fourth haplotype, which had a novel restriction site pattern and noncoding sequence, was found in samples from the extreme south, near the Syrian border. We found two different noncoding sequences among the eastern Mediterranean haplotypes. The "Caucasian" sequence matches that described from A. m. caucasica, and the "Anatolian" sequence matches that of A. m. carnica. The frequency of the "Caucasian" sequence was highest (98-100%) in sites near the Georgian border and decreased steeply to the south and west. Elsewhere the Anatolian sequence was found. In European Turkey (Thrace) a restriction site polymorphism previously reported from A. m. carnica in Austria and the Balkans was present at high frequency. A novel mtDNA haplotype with a unique restriction site pattern and noncoding sequence was found among bees from Hatay, in the extreme south near the Syrian border. This haplotype differed from the three previously known lineages of honeybee mtDNA--African, western European, and eastern Mediterranean-and may represent a fourth mitochondrial lineage.     

The impact of apiculture on the genetic structure of wild honeybee populations (<Emphasis Type="Italic">Apis mellifera</Emphasis>) in Sudan

Apiculture often relies on the importation of non-native honeybees (Apis mellifera) and large distance migratory beekeeping. These activities can cause biodiversity conflicts with the conservation of wild endemic honeybee subspecies. We studied the impact of large scale honeybee imports on managed and wild honeybee populations in Sudan, a centre of biodiversity of A. mellifera, using as set of linked microsatellite DNA loci and mitochondrial DNA markers. The mitochondrial DNA analyses showed that all wild honey bees exclusively belonged to African haplotypes. However, we revealed non-native haplotypes in managed colonies on apiaries reflecting unambiguous evidence of imports from European stock. Moreover, we found significantly higher linkage disequilibria for microsatellite markers in wild populations in regions with apiculture compared to wild populations which had no contact to beekeeping. Introgression of imported honeybees was only measurable at the population level in close vicinity to apicultural activities but not in wild populations which represent the vast majority of honeybees in Sudan.     

Apis mellifera evolutionary lineages in Northern Africa: Libya, where orient meets occident

Abstract  The distribution of various evolutionary lineages of Apis mellifera subspecies in Africa is still controversial. We sampled honeybees from eight coastal locations and three Saharan oases in Libya and analyzed mtDNA variability with restriction fragment length polymorphisms and the sequence of the tRNAleu-cox2 intergenic region. Haplotypes belonging to the oriental O evolutionary lineage, including four which are newly described, were detected in all investigated locations. Haplotypes belonging to the European M lineage were rarely detected, probably reflecting the effect of sporadic importations. Honeybees belonging to the A lineage were detected in Al Aziziyah and Zlitan close to the Tunisian border. The distribution of the O lineage extends westward up to the border between Libya and Tunisia, a contact area between the O and A lineages. Various Libyan honeybee populations in Saharan oases are characterized by novel and unique haplotypes (O4, O5, O5′ and O5″). These might be natural relic populations that became isolated when the North African Sahara desert was still grassland (0.126–0.168 Myr ago).     

Mitochondrial DNA Diversity of Honey Bees (Apis mellifera) from Unmanaged Colonies and Swarms in the United States

To study the genetic diversity of honey bees (Apis mellifera L.) from unmanaged colonies in the United States, we sequenced a portion of the mitochondrial DNA COI–COII region. From the 530 to 1,230 bp amplicon, we observed 23 haplotypes from 247 samples collected from 12 states, representing three of the four A. mellifera lineages known to have been imported into the United States (C, M, and O). Six of the 13 C lineage haplotypes were not found in previous queen breeder studies in the United States. The O lineage accounted for 9% of unmanaged colonies which have not yet been reported in queen breeder studies. The M lineage accounted for a larger portion of unmanaged samples (7%) than queen breeder samples (3%). Based on our mitochondrial DNA data, the genetic diversity of unmanaged honey bees in the United States differs significantly from that of queen breeder populations (p < 0.00001). The detection of genetically distinct maternal lineages of unmanaged honey bees suggests that these haplotypes may have existed outside the managed honey bee population for a long period.     

Local honeybee (<Emphasis Type="Italic">Apis mellifera mellifera</Emphasis> L.) populations in the Urals

The COI-COII intergenic region of mitochondrial DNA (mtDNA) was studied in local honeybee (Apis mellifera mellifera) L. populations from the Middle and Southern Urals. Analysis of bee colonies in these regions revealed apiaries enriched in families descending from A. m. mellifera in the maternal lineage. These results confirm the suggestion of preservation of A. m. mellifera refuges in the Urals and provide grounds for work on the preservation of the gene pool of this bee variety, valuable for all Russia.     

Geographical patterns of mitochondrial DNA variation in Apis mellifera iberiensis (Hymenoptera: Apidae)

An extensive survey of mitochondrial haplotypes in honeybee colonies from the Iberian Peninsula has corroborated previous hypotheses about the existence of a joint clinal variation of African (A) and west European (M) evolutionary lineages. It has been found that the Iberian Peninsula is the European region with the highest haplotype diversity (12 haplotypes detected of the M lineage and 10 of the A lineage). The frequency of A haplotypes decreases in a SW-NE trend, while that of M haplotypes increases. These results are discussed in relation to hypotheses about the African origin of Apis mellifera and an early colonization of west Europe during intermediate Pleistocene glaciation events, followed by a regional differentiation. The extant pattern of haplotype frequency and distribution seems to be influenced at a regional scale by adaptation to local climatic conditions and the mobile beekeeping that has become a large-scale practice during the last decades. Other previous anthropogenic influences (Greek, Roman and Arab colonizations) are thought to be of minor importance in present day populations.     

Genetic reincarnation of workers as queens in the Eastern honeybee Apis cerana

Thelytokous parthenogenesis, or the asexual production of female offspring, is rare in the animal kingdom, but relatively common in social Hymenoptera. However, in honeybees, it is only known to be ubiquitous in one subspecies of Apis mellifera, the Cape honeybee, A. mellifera capensis. Here we report the appearance of queen cells in two colonies of the Eastern honeybee Apis cerana that no longer contained a queen or queen-produced brood to rear queens from. A combination of microsatellite genotyping and the timing of the appearance of these individuals excluded the possibility that they had been laid by the original queen. Based on the genotypes of these individuals, thelytokous production by natal workers is the most parsimonious explanation for their existence. Thus, we present the first example of thelytoky in a honeybee outside A. mellifera. We discuss the evolutionary and ecological consequences of thelytoky in A. cerana, in particular the role thelytoky may play in the recent invasions by populations of this species.     

Discrimination of Iranian honeybee populations (Apis mellifera meda) from commercial subspecies of Apis mellifera L. using morphometric and genetic methods

The morphological characters of honeybees have an important role for discriminating honeybee subspecies. In the present research, Iranian populations of honeybee (Apis mellifera) were collected from 19 areas in Iran. The samples were collected from stationary beekeeping sites. Moreover, pictures of honeybee forewings held in the Bee Data Bank in Oberursel were compared with Iranian honeybee populations. 19 morphological characters were measured for each forewing of worker honeybee to evaluate differentiation of Iranian honeybee populations from the commercial honeybee subspecies A. m. mellifera, A. m. carnica, A. m. caucasica and A. m. ligustica. Additionally, part of the tRNA^leu gene, an intergenic region and part of COII was used to confirm differentiation of the commercial subspecies from Iranian honeybee populations. Results of the cluster analyses showed that 19 morphological characters of forewings differentiated Iranian populations from the commercial subspecies. Moreover, the phylogenetic tree of part of the tRNA^leu gene, an intergenic region and part of COII differentiated the commercial subspecies from Iranian honeybee populations. Results of the discriminant function analyses (DFA) indicated that the references samples of A. m. meda overlapped with Iranian populations.     

Varying degrees of Apis mellifera ligustica introgression in protected populations of the black honeybee, Apis mellifera mellifera, in northwest Europe

The natural distribution of honeybee subspecies in Europe has been significantly affected by human activities during the last century. Non-native subspecies of honeybees have been introduced and propagated, so that native black honeybee (Apis mellifera mellifera) populations lost their identity by gene-flow or went extinct. After previous studies investigated the remaining gene-pools of native honeybees in France and Spain, we here assess the genetic composition of eight northwest European populations of the black honeybee, using both mitochondrial (restriction fragment length polymorphisms of the intergenic transfer RNAleu-COII region) and nuclear (11 microsatellite loci) markers. Both data sets show that A. m. mellifera populations still exist in Norway, Sweden, Denmark, England, Scotland and Ireland, but that they are threatened by gene flow from commercial honeybees. Both Bayesian admixture analysis of the microsatellite data and DraI-RFLP (restriction fragment length polymorphism) analysis of the intergenic region indicated that gene-flow had hardly occurred in some populations, whereas almost 10% introgression was observed in other populations. The most introgressed population was found on the Danish Island of Laeso, which is the last remaining native Danish population of A. m. mellifera and the only one of the eight investigated populations that is protected by law. We discuss how individual admixture analysis can be used to monitor the restoration of honeybee populations that suffer from unwanted hybridization with non-native subspecies.     

Determining colony densities in wild honeybee populations (<Emphasis Type="Italic">Apis mellifera</Emphasis>) with linked microsatellite DNA markers

Estimating the population size of social bee colonies in the wild is often difficult because nests are highly cryptic. Because of the honeybee (Apis mellifera) mating behaviour, which is characterized by multiple mating of queens at drone congregation areas (DCA), it is possible to use genotypes of drones caught at these areas to infer the number of colonies in a given region. However, DCAs are difficult to locate and we assess the effectiveness of an alternative sampling technique to determine colony density based on inferring male genotypes from queen offspring. We compare these methods in the same population of wild honeybees, Apis mellifera scutellata. A set of linked microsatellite loci is used to decrease the frequency of recombination among marker loci and therefore increase the precision of the estimates. Estimates of population size obtained through sampling of queen offspring is significantly larger than that obtained by sampling drones at DCAs. This difference may be due to the more extensive flying range of queens compared with drones on mating flights. We estimate that the population size sampled through queen offspring is about double that sampled through drones.     

A SNP test to identify Africanized honeybees via proportion of ‘African’ ancestry

The honeybee, Apis mellifera, is the world's most important pollinator and is ubiquitous in most agricultural ecosystems. Four major evolutionary lineages and at least 24 subspecies are recognized. Commercial populations are mainly derived from subspecies originating in Europe (75–95%). The Africanized honeybee is a New World hybrid of A. m. scutellata from Africa and European subspecies, with the African component making up 50–90% of the genome. Africanized honeybees are considered undesirable for bee‐keeping in most countries, due to their extreme defensiveness and poor honey production. The international trade in honeybees is restricted, due in part to bans on the importation of queens (and semen) from countries where Africanized honeybees are extant. Some desirable strains from the United States of America that have been bred for traits such as resistance to the mite Varroa destructor are unfortunately excluded from export to countries such as Australia due to the presence of Africanized honeybees in the USA. This study shows that a panel of 95 single nucleotide polymorphisms, chosen to differentiate between the African, Eastern European and Western European lineages, can detect Africanized honeybees with a high degree of confidence via ancestry assignment. Our panel therefore offers a valuable tool to mitigate the risks of spreading Africanized honeybees across the globe and may enable the resumption of queen and bee semen imports from the Americas.     

A fifth major genetic group among honeybees revealed in Syria

Background

Apiculture has been practiced in North Africa and the Middle-East from antiquity. Several thousand years of selective breeding have left a mosaic of Apis mellifera subspecies in the Middle-East, many uniquely adapted and survived to local environmental conditions. In this study we explore the genetic diversity of A. mellifera from Syria (n?=?1258), Lebanon (n?=?169) and Iraq (n?=?35) based on 14 short tandem repeat (STR) loci in the context of reference populations from throughout the Old World (n?=?732).

Results

Our data suggest that the Syrian honeybee Apis mellifera syriaca occurs in both Syrian and Lebanese territories, with no significant genetic variability between respective populations from Syria and Lebanon. All studied populations clustered within a new fifth independent nuclear cluster, congruent with an mtDNA Z haplotype identified in a previous study. Syrian honeybee populations are not associated with Oriental lineage O, except for sporadic introgression into some populations close to the Turkish and Iraqi borders. Southern Syrian and Lebanese populations demonstrated high levels of genetic diversity compared to the northern populations.

Conclusion

This study revealed the effects of foreign queen importations on Syrian bee populations, especially for the region of Tartus, where extensive introgression of A. m. anatolica and/or A. m. caucasica alleles were identified. The policy of creating genetic conservation centers for the Syrian subspecies should take into consideration the influence of the oriental lineage O from the northern Syrian border and the large population of genetically divergent indigenous honeybees located in southern Syria.

     

Prevalence of African DNA RELP alleles in neotropical African honeybees

A few queens of the honeybee, Apis mellifera scutellata, were imported from Africa and released in Brazil in 1957. Progeny of these bees have now largely colonized the American tropics. Their imminent arrival in the United States poses a serious threat to the beekeeping industry and to agriculture dependent on honeybee pollination. African and European bees are morphologically very similar. DNA restriction fragment length polymorphisms are proving successful in distinguishing between the two. Several DNA markers specific to European honeybees have been described previously. Reported here are three cloned honeybee DNA probes that reveal polymorphisms that appear to be either African or European specific. Of fourteen alleles or haplotypes identified, five were present only in African and neotropical (Venezuelan and Mexican) African bees but absent in European-derived bees, two were present only in European-derived bees but absent in samples from South Africa. Another allele showed apparent frequency differences among populations. Such markers are useful in studying the genetics of neotropical African bee populations. Venezuelan and Mexican honeybee colonies show a preponderance of the African alleles with low levels of the European alleles. These observations of nuclear DNA, revealing limited paternal European introgression, together with previous mitochondrial DNA findings showing negligible European maternal gene flow into feral African populations, indicate that neotropical African bees are primarily African.     

Hybrid origins of honeybees from italy (Apis mellifera ligustica) and sicily (A. m. sicula)

The genetic variability of honeybee populations Apis mellifera ligustica, in continental Italy, and of A. m. sicula, in Sicily, was investigated using nuclear (microsatellite) and mitochondrial markers. Six populations (236 individual bees) and 17 populations (664 colonies) were, respectively, analysed using eight microsatellite loci and DraI restriction fragment length polymorphism (RFLP) of the cytochrome oxidase I (COI)-cytochrome oxidase II (COII) region. Microsatellite loci globally confirmed the southeastern European heritage of both subspecies (evolutionary branch C). However, A. m. ligustica mitochondrial DNA (mtDNA) appeared to be a composite of the two European (M and C) lineages over most of the Italian peninsula, and only mitotypes from the African (A) lineage were found in A. m. sicula samples. This demonstrates a hybrid origin for both subspecies. For A. m. ligustica, the most widely exported subspecies, this hybrid origin has long been obscured by the fact that in the main area of queen production (from which most of the previous ligustica bee samples originated) the M mitochondrial lineage is absent, whereas it is present almost everywhere else in Italy. This presents a new view of the evolutionary history of European honeybees. For instance, the Iberian peninsula was considered as the unique refuge for the M branch during the quaternary ice periods. Our results show that the Apennine peninsula played a similar role. The differential distribution of nuclear and mitochondrial markers observed in Italy seems to be a general feature of introgressed honeybee populations. Presumably, it stems from the social nature of the species in which both genome compartments are differentially affected by the two (individual and colonial) reproduction levels.     

Diversity of Apis mellifera Subspecies from Turkey Revealed by Sequence Analysis of Mitochondrial 16s rDNA Region

Mitochondrial DNA sequence variation can be used to infer honeybee evolutionary relationships. In this study, DNA sequence diversity of the mitochondrial 16s rDNA region was investigated in 112 honeybees from 15 populations in Turkey, which is mainly populated with Apis mellifera anatoliaca, A. m. caucasica, and A. m. meda. The study revealed 11 haplotypes for this segment, with 13 variable sites and nine parsimony informative sites. The haplotypes were not discriminated according to their geographical locations in a neighbor-joining dendrogram based on 16s rDNA sequences available in Genbank, but all the haplotypes obtained in this study are clustered with published haplotypes such as A. mellifera TAS (AF214666) and A. m. ligustica (EF116868) and with some unpublished Genbank records (HQ318928, HQ318934, and HQ318938). This study expands the knowledge of the mitochondrial 16s rDNA region, and it presents the first comprehensive sequence analysis of this region in Turkish honeybees.     

Radiation of Atlantic goldcrests Regulus regulus spp.: evidence of a new taxon from the Canary Islands

Phylogenetic relationships between goldcrest populations from the Atlantic Islands (Azores and Canary Islands) were investigated by two molecular markers (mitochondrial control region and cytochrome b sequences), and partly by morphology and territorial song. The Azorean goldcrest populations are closely related to European nominate R. r. regulus. Most probably, the Azores were colonized by goldcrests in a single late‐pleistocene invasion, while colonization of the Canary Islands presumably occurred in two steps: An early invasion to Tenerife and La Gomera 1.9–2.3 million years (my) ago and a more recent one to El Hierro and La Palma 1.3–1.8 my ago. Distribution of haplotypes on the Azores suggests a division of R. r. azoricus on São Miguel into an eastern population with close affinities to R. r. sanctaemariae and a western population belonging to the lineage of R. r. inermis on the central and western island group. The Canarian populations are genetically substructured into a northeastern group embracing Tenerife and La Gomera and a second, southwestern group including El Hierro and La Palma. Genetic distances between members of the two Canarian clades range at 3.1–3.4% (TrN distance, control region and cytochrome b). Differentiation between the two groups is also supported by morphology and by territorial song. Substitution rate estimates for the both genes range at approximately the same values of 0.0031 and 0.0044 substitutions per site and lineage per my which roughly corresponds 0.61–0.83% divergence between Regulus lineages per my. Highest local rates occur in island clades of the Azorean and the Canarian population and in R. r. japonensis from the Russian Far East and Japan. However, a general acceleration of a molecular clock in island populations is not evident from the Regulus data set due to extremely low local rate estimates in the Canarian clade of Tenerife and La Gomera. As a taxonomic consequence of the marked differentiation of the two Canarian goldcrest clades the populations from El Hierro and La Palma are described as a taxon new to science and are named Regulus regulus ellenthalerae n. ssp.     

Introgression in native populations of <Emphasis Type="Italic">Apis mellifera mellifera</Emphasis> L: implications for conservation

Hybridisation and introgression can have negative impacts on regional biodiversity through the potential erosion of locally adapted lineages. The honey bee (Apis mellifera L.) occurs in twenty-seven subspecies across Europe, is an extremely economically important insect, yet threatened by multifarious impacts. Transhumance of the most commercially appealing varieties threatens native honey bee diversity by introgression and subsequent loss of locally adapted traits, or even by complete removal of some subspecies from parts of the range. Here levels of admixture and introgression are examined in UK honey bees suspected to be from hives of the dark European honey bee (Apis mellifera mellifera). Microsatellite DNA and STRUCTURE analyses reveal that the studied populations are generally admixed, and discriminant analysis of principal components shows them to be intermediate between A. m. mellifera and Apis mellifera carnica populations. However, examining mitochondrial haplotype data (COI-COII intergenic spacer region) and nuclear DNA reveal that some hives are relatively pure (from 4 to 15 hives, depending on the Q-value threshold). Genetic diversity is relatively high in comparison with other European populations. Implications for conservation and management are discussed.     

Molecular genetic analysis of five extant reserves of black honeybee <Emphasis Type="Italic">Apis melifera melifera</Emphasis> in the Urals and the Volga region

Local populations of the black honeybee Apis mellifera mellifera from the Urals and the Volga region were examined in comparison with local populations of southern honeybee subspecies A. m. caucasica and A. m. carpatica from the Caucasus and the Carpathians. Genetic analysis was performed on the basis of the polymorphism of nine microsatellite loci of nuclear DNA and the mtDNA COI–COII locus. On the territory of the Urals and the Volga region, five extant populations (reserves) of the black honeybee A. m. mellifera were identified, including the Burzyanskaya, Tatyshlinskaya, Yuzhno-Prikamskaya, Visherskaya, and Kambarskaya populations. These five populations are the basis of the modern gene pool of the black honeybee A. m. mellifera from the Urals and the Volga region. The greatest proportion of the remaining indigenous gene pool of A. m. mellifera (the core of the gene pool of the population of A. m. mellifera) is distributed over the entire territory of Perm krai and the north of the Republic of Bashkortostan. For the population of A. m. mellifera from the Urals and the Volga region, the genetic standards were calculated, which will be useful for future population studies of honeybees.