TEMPORAL PATTERN OF AFRICANIZATION IN A FERAL HONEYBEE POPULATION FROM TEXAS INFERRED FROM MITOCHONDRIAL DNA

Historical Changes in Honey Bee Wing Venation in Romania

Simple Summary

Honey bees, in addition to producing honey, are important pollinators of wild and cultivated plants. Unfortunately, in some places, the population of honey bees is declining. One of the factors that affect their survival is adaptation to the local environment. Bees native to a particular area survive better than those imported. Despite this fact, some beekeepers import non-native bees and use them in their apiaries. Imported bees produce hybrids with bees from surrounding colonies because beekeepers do not control their mating. In consequence, the whole population can change. In this study, we verified how the population of Romanian bees has changed over the last four decades. We found significant temporal changes in wing venation. Despite these changes, the two major subpopulations of bees separated by mountains remain distinct. We provide a tool for the easy identification of native bees from Romania, which can help to protect them.

Abstract

The honey bee (Apis mellifera) is an ecologically and economically important species that provides pollination services to natural and agricultural systems. The biodiversity of the honey bee is being endangered by the mass import of non-native queens. In many locations, it is not clear how the local populations have been affected by hybridisation between native and non-native bees. There is especially little information about temporal changes in hybridisation. In Romania, A. m. carpatica naturally occurs, and earlier studies show that there are two subpopulations separated by the Carpathian Mountains. In this study, we investigated how the arrangement of veins in bees’ wings (venation) has changed in Romanian honey bees in the last four decades. We found that in the contemporary population of Romanian bees, there are still clear differences between the intra- and extra-Carpathian subpopulations, which indicates that natural variation among honey bees is still being preserved. We also found significant differences between bees collected before and after 2000. The observed temporal changes in wing venation are most likely caused by hybridisation between native bees and non-native bees sporadically introduced by beekeepers. In order to facilitate conservation and the monitoring of native Romanian bees, we developed a method facilitating their identification.

Defense Response in Brazilian Honey Bees (Apis mellifera scutellata × spp.) Is Underpinned by Complex Patterns of Admixture

In 1957, an invasive and highly defensive honey bee began to spread across Brazil. In the previous year, Brazilian researchers hoped to produce a subtropical-adapted honey bee by crossing local commercial honey bees (of European origin) with a South African honey bee subspecies (Apis mellifera scutellata; an A-lineage honey bee subspecies). The resulting cross-African hybrid honey bees (AHBs)-escaped from their enclosure and spread through the Americas. Today, AHB is the most common honey bee from Northern Argentina to the Southern United States. AHBs are much more likely to sting nest intruders than managed European-derived honey bee colonies. Previous studies have explored how genetic variation contributes to differences in defense response between European-derived honey bee and AHB. Although this work demonstrated very strong genetic effects on defense response, they have yet to pinpoint which genes influence variation in defense response within AHBs, specifically. We quantified defense response for 116 colonies in Brazil and performed pooled sequencing on the most phenotypically divergent samples. We identified 65 loci containing 322 genes that were significantly associated with defense response. Loci were strongly associated with metabolic function, consistent with previous functional genomic analyses of this phenotype. Additionally, defense-associated loci had nonrandom and unexpected patterns of admixture. Defense response was not simply the product of more A-lineage honey bee ancestry as previously assumed, but rather an interaction between A-lineage and European alleles. Our results suggest that a combination of A-lineage and European alleles play roles in defensive behavior in AHBs.

The Economics of Honey Bee (Hymenoptera: Apidae) Management and Overwintering Strategies for Colonies Used to Pollinate Almonds

Commercial honey bee (Apis mellifera L.) colonies significantly contribute to agricultural productivity through crop pollination. Almond production requires the most colonies because there are more than a million acres of orchards that require cross-pollination for nut set. With the rising costs of managing and transporting colonies to almond orchards combined with the high colony losses beekeepers routinely experience, we asked if renting colonies for almond pollination was profitable. We conducted a longitudinal study on 190 colonies from their establishment in April until they were placed in almond orchards 10 mo later. In the fall, equal numbers of colonies were placed either in cold storage (CS) facilities or in outdoor apiaries for the winter. We found that the cost of overwintering colonies in CS was lower than in apiaries, but CS did not reduce overwintering losses. A key finding from our study is that there is little or no profit in renting colonies for almond pollination once summer management and overwintering costs are considered. Our only profitable venture was honey production in the summer. We propose alternative management strategies to lower costs and make almond pollination profitable. We also developed a decision tool for selecting colonies to overwinter in CS and reduce expenditures on those that will not reach sufficient size for almond pollination. Our study exposes the unsustainable financial burden experienced by migratory beekeepers that is not included in estimates of yearly colony losses, and underscores the urgent need for forage plantings to generate revenue from honey and improve overwinter survival.

Genome Sequencing of Museum Specimens Reveals Rapid Changes in the Genetic Composition of Honey Bees in California

The western honey bee, Apis mellifera, is an enormously influential pollinator in both natural and managed ecosystems. In North America, this species has been introduced numerous times from a variety of different source populations in Europe and Africa. Since then, feral populations have expanded into many different environments across their broad introduced range. Here, we used whole genome sequencing of historical museum specimens and newly collected modern populations from California (USA) to analyze the impact of demography and selection on introduced populations during the past 105 years. We find that populations from both northern and southern California exhibit pronounced genetic changes, but have changed in different ways. In northern populations, honey bees underwent a substantial shift from western European to eastern European ancestry since the 1960s, whereas southern populations are dominated by the introgression of Africanized genomes during the past two decades. Additionally, we identify an isolated island population that has experienced comparatively little change over a large time span. Fine-scale comparison of different populations and time points also revealed SNPs that differ in frequency, highlighting a number of genes that may be important for recent adaptations in these introduced populations.

Africanization of a feral honey bee (Apis mellifera) population in South Texas: does a decade make a difference?

The arrival to the United States of the Africanized honey bee, a hybrid between European subspecies and the African subspecies Apis mellifera scutellata, is a remarkable model for the study of biological invasions. This immigration has created an opportunity to study the dynamics of secondary contact of honey bee subspecies from African and European lineages in a feral population in South Texas. An 11‐year survey of this population (1991–2001) showed that mitochondrial haplotype frequencies changed drastically over time from a resident population of eastern and western European maternal ancestry, to a population dominated by the African haplotype. A subsequent study of the nuclear genome showed that the Africanization process included bidirectional gene flow between European and Africanized honey bees, giving rise to a new panmictic mixture of A. m. scutellata‐ and European‐derived genes. In this study, we examined gene flow patterns in the same population 23 years after the first hybridization event occurred. We found 28 active colonies inhabiting 92 tree cavities surveyed in a 5.14 km² area, resulting in a colony density of 5.4 colonies/km². Of these 28 colonies, 25 were of A. m. scutellata maternal ancestry, and three were of western European maternal ancestry. No colonies of eastern European maternal ancestry were detected, although they were present in the earlier samples. Nuclear DNA revealed little change in the introgression of A. m. scutellata‐derived genes into the population compared to previous surveys. Our results suggest this feral population remains an admixed swarm with continued low levels of European ancestry and a greater presence of African‐derived mitochondrial genetic composition.

Range and Frequency of Africanized Honey Bees in California (USA)

Africanized honey bees entered California in 1994 but few accounts of their northward expansion or their frequency relative to European honey bees have been published. We used mitochondrial markers and morphometric analyses to determine the prevalence of Africanized honeybees in San Diego County and their current northward progress in California west of the Sierra Nevada crest. The northernmost African mitotypes detected were approximately 40 km south of Sacramento in California’s central valley. In San Diego County, 65% of foraging honey bee workers carry African mitochondria and the estimated percentage of Africanized workers using morphological measurements is similar (61%). There was no correlation between mitotype and morphology in San Diego County suggesting Africanized bees result from bidirectional hybridization. Seventy percent of feral hives, but only 13% of managed hives, sampled in San Diego County carried the African mitotype indicating that a large fraction of foraging workers in both urban and rural San Diego County are feral. We also found a single nucleotide polymorphism at the DNA barcode locus COI that distinguishes European and African mitotypes. The utility of this marker was confirmed using 401 georeferenced honey bee sequences from the worldwide Barcode of Life Database. Future censuses can determine whether the current range of the Africanized form is stable, patterns of introgression at nuclear loci, and the environmental factors that may limit the northern range of the Africanized honey bee.

Museum samples reveal rapid evolution by wild honey bees exposed to a novel parasite

Understanding genetic changes caused by novel pathogens and parasites can reveal mechanisms of adaptation and genetic robustness. Using whole-genome sequencing of museum and modern specimens, we describe the genomic changes in a wild population of honey bees in North America following the introduction of the ectoparasitic mite, Varroa destructor. Even though colony density in the study population is the same today as in the past, a major loss of haplotypic diversity occurred, indicative of a drastic mitochondrial bottleneck, caused by massive colony mortality. In contrast, nuclear genetic diversity did not change, though hundreds of genes show signs of selection. The genetic diversity within each bee colony, particularly as a consequence of polyandry by queens, may enable preservation of genetic diversity even during population bottlenecks. These findings suggest that genetically diverse honey bee populations can recover from introduced diseases by evolving rapid tolerance, while maintaining much of the standing genetic variation.

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.

Estimating the density of honeybee colonies across their natural range to fill the gap in pollinator decline censuses

Although pollinator declines are a global biodiversity threat, the demography of the western honeybee (Apis mellifera) has not been considered by conservationists because it is biased by the activity of beekeepers. To fill this gap in pollinator decline censuses and to provide a broad picture of the current status of honeybees across their natural range, we used microsatellite genetic markers to estimate colony densities and genetic diversity at different locations in Europe, Africa, and central Asia that had different patterns of land use. Genetic diversity and colony densities were highest in South Africa and lowest in Northern Europe and were correlated with mean annual temperature. Confounding factors not related to climate, however, are also likely to influence genetic diversity and colony densities in honeybee populations. Land use showed a significantly negative influence over genetic diversity and the density of honeybee colonies over all sampling locations. In Europe honeybees sampled in nature reserves had genetic diversity and colony densities similar to those sampled in agricultural landscapes, which suggests that the former are not wild but may have come from managed hives. Other results also support this idea: putative wild bees were rare in our European samples, and the mean estimated density of honeybee colonies on the continent closely resembled the reported mean number of managed hives. Current densities of European honeybee populations are in the same range as those found in the adverse climatic conditions of the Kalahari and Saharan deserts, which suggests that beekeeping activities do not compensate for the loss of wild colonies. Our findings highlight the importance of reconsidering the conservation status of honeybees in Europe and of regarding beekeeping not only as a profitable business for producing honey, but also as an essential component of biodiversity conservation.

Africanization in the United States: replacement of feral European honeybees (Apis mellifera L.) by an African hybrid swarm

The expansion of Africanized honeybees from South America to the southwestern United States in <50 years is considered one of the most spectacular biological invasions yet documented. In the American tropics, it has been shown that during their expansion Africanized honeybees have low levels of introgressed alleles from resident European populations. In the United States, it has been speculated, but not shown, that Africanized honeybees would hybridize extensively with European honeybees. Here we report a continuous 11-year study investigating temporal changes in the genetic structure of a feral population from the southern United States undergoing Africanization. Our microsatellite data showed that (1) the process of Africanization involved both maternal and paternal bidirectional gene flow between European and Africanized honeybees and (2) the panmitic European population was replaced by panmitic mixtures of A. m. scutellata and European genes within 5 years after Africanization. The post-Africanization gene pool (1998-2001) was composed of a diverse array of recombinant classes with a substantial European genetic contribution (mean 25-37%). Therefore, the resulting feral honeybee population of south Texas was best viewed as a hybrid swarm.

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.