Admixture increases diversity in managed honey bees: Reply to De la Rúa et al. (2013)

De la Rúa et al. (2013) express some concerns about the conclusions of our recent study showing that management increases genetic diversity of honey bees (Apis mellifera) by promoting admixture (Harpur et al. 2012). We provide a brief review of the literature on the population genetics of A. mellifera and show that we utilized appropriate sampling methods to estimate genetic diversity in the focal populations. Our finding of higher genetic diversity in two managed A. mellifera populations on two different continents is expected to be the norm given the large number of studies documenting admixture in honey bees. Our study focused on elucidating how management affects genetic diversity in honey bees, not on how to best manage bee colonies. We do not endorse the intentional admixture of honey bee populations, and we agree with De la Rúa et al. (2013) that native honey bee subspecies should be conserved.     

Higher immunocompetence is associated with higher genetic diversity in feral honey bee colonies (<Emphasis Type="Italic">Apis mellifera</Emphasis>)

Honey bees are the most important managed pollinators as they provide key ecosystem services for crop production worldwide. Recent losses of honey bee colonies in North America and Europe have demonstrated a need to develop strategies to improve their health and conserve their populations. Previously, we showed that feral honey bees—colonies that live in the wild without human assistance—exhibit higher levels of immunocompetence than managed colonies in North Carolina (USA). In a first attempt to investigate the underlying mechanisms of this difference in immune response, here we characterize the genetic composition of feral and managed honey bees using microsatellite markers. Our results reveal significant but small genetic differentiation between feral and managed honey bee colonies (?_CT = 0.047, P?=?0.03) indicating admixture between these two groups. Higher genetic diversity was correlated with higher immune response in feral (P _MANOVA = 0.011) but not managed bees, despite the fact that the latter group showed significantly higher average genetic diversity (P _ANCOVA < 0.001). These findings suggest that genetic diversity is positively associated with immunocompetence in feral honey bee colonies, and that the benefits of genetic diversity are obscured in managed bees, perhaps as a result of artificial selection. We hypothesize that high genetic variability provides the raw material upon which natural selection acts and generates adaptive genotypes in unmanaged populations. Feral populations could be useful sources of genetic variation to use in breeding programs that aim to improve honey bee health.     

Conserving genetic diversity in the honeybee: Comments on Harpur et al. (2012)

The article by Harpur et al. (2012) ‘Management increases genetic diversity of honey bees via admixture’ concludes that ‘…honey bees do not suffer from reduced genetic diversity caused by management and, consequently, that reduced genetic diversity is probably not contributing to declines of managed Apis mellifera populations’. In the light of current honeybee and beekeeping declines and their consequences for honeybee conservation and the pollination services they provide, we would like to express our concern about the conclusions drawn from the results of Harpur et al. (2012). While many honeybee management practices do not imply admixture, we are convinced that the large-scale genetic homogenization of admixed populations could drive the loss of valuable local adaptations. We also point out that the authors did not account for the extensive gene flow that occurs between managed and wild/feral honeybee populations and raise concerns about the data set used. Finally, we caution against underestimating the importance of genetic diversity for honeybee colonies and highlight the importance of promoting the use of endemic honeybee subspecies in apiculture.     

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.     

Farming practices influence wild pollinator populations on squash and pumpkin

Recent declines in managed honey bee, Apis mellifera L., colonies have increased interest in the current and potential contribution of wild bee populations to the pollination of agricultural crops. Because wild bees often live in agricultural fields, their population density and contribution to crop pollination may be influenced by farming practices, especially those used to reduce the populations of other insects. We took a census of pollinators of squash and pumpkin at 25 farms in Virginia, West Virginia, and Maryland to see whether pollinator abundance was related to farming practices. The main pollinators were Peponapis pruinosa Say; honey bees, and bumble bees (Bombus spp.). The squash bee was the most abundant pollinator on squash and pumpkin, occurring at 23 of 25 farms in population densities that were commonly several times higher than that of other pollinators. Squash bee density was related to tillage practices: no-tillage farms hosted three times as great a density of squash bees as tilled farms. Pollinator density was not related to pesticide use. Honey bee density on squash and pumpkin was not related to the presence of managed honey bee colonies on farms. Farms with colonies did not have more honey bees per flower than farms that did not keep honey bees, probably reflecting the lack of affinity of honey bees for these crops. Future research should examine the economic impacts of managing farms in ways that promote pollinators, particularly pollinators of crops that are not well served by managed honey bee colonies.     

Admixture in Africanized honey bees (Apis mellifera) from Panamá to San Diego,California (U.S.A.)

The Africanized honey bee (AHB) is a New World amalgamation of several subspecies of the western honey bee (Apis mellifera), a diverse taxon historically grouped into four major biogeographic lineages: A (African), M (Western European), C (Eastern European), and O (Middle Eastern). In 1956, accidental release of experimentally bred “Africanized” hybrids from a research apiary in Sao Paulo, Brazil initiated a hybrid species expansion that now extends from northern Argentina to northern California (U.S.A.). Here, we assess nuclear admixture and mitochondrial ancestry in 60 bees from four countries (Panamá; Costa Rica, Mexico; U.S.A) across this expansive range to assess ancestry of AHB several decades following initial introduction and test the prediction that African ancestry decreases with increasing latitude. We find that AHB nuclear genomes from Central America and Mexico have predominately African genomes (76%–89%) with smaller contributions from Western and Eastern European lineages. Similarly, nearly all honey bees from Central America and Mexico possess mitochondrial ancestry from the African lineage with few individuals having European mitochondria. In contrast, AHB from San Diego (CA) shows markedly lower African ancestry (38%) with substantial genomic contributions from all four major honey bee lineages and mitochondrial ancestry from all four clades as well. Genetic diversity measures from all New World populations equal or exceed those of ancestral populations. Interestingly, the feral honey bee population of San Diego emerges as a reservoir of diverse admixture and high genetic diversity, making it a potentially rich source of genetic material for honey bee breeding.     

EVOLUTION AND POPULATION STRUCTURE OF AFRICANIZED HONEY BEES IN BRAZIL: EVIDENCE FROM SPATIAL ANALYSIS OF MORPHOMETRIC DATA

In recent years, studies based on isoenzymatic patterns of geographic variation have revealed that what is usually called the Africanized honey bee does not constitute a single population. Instead, several local populations exist with various degrees of admixture with European honey bees. In this paper, we evaluated new data on morphometric patterns of Africanized honey bees collected at 42 localities in Brazil, using univariate and multivariate (canonical) trend surface and spatial autocorrelation analyses. The clinal patterns of variation found for genetically independent characters (wing size characters and some wing venation angles) are concordant with previous studies of malate dehydrogenase (MDH) allelic frequencies and support the hypothesis that larger honey bees in southern and southeastern Brazil originated by racial admixture in the initial phases of African honey bee colonization. Geographic variation patterns of Africanized honey bee populations reflect a demic diffusion process in which European genes were gradually lost because of the higher fitness of the African gene pool in Neotropical environmental conditions.     

Evaluation of the Distribution and Impacts of Parasites,Pathogens, and Pesticides on Honey Bee (Apis mellifera) Populations in East Africa

In East Africa, honey bees (Apis mellifera) provide critical pollination services and income for small-holder farmers and rural families. While honey bee populations in North America and Europe are in decline, little is known about the status of honey bee populations in Africa. We initiated a nationwide survey encompassing 24 locations across Kenya in 2010 to evaluate the numbers and sizes of honey bee colonies, assess the presence of parasites (Varroa mites and Nosema microsporidia) and viruses, identify and quantify pesticide contaminants in hives, and assay for levels of hygienic behavior. Varroa mites were present throughout Kenya, except in the remote north. Levels of Varroa were positively correlated with elevation, suggesting that environmental factors may play a role in honey bee host-parasite interactions. Levels of Varroa were negatively correlated with levels of hygienic behavior: however, while Varroa infestation dramatically reduces honey bee colony survival in the US and Europe, in Kenya Varroa presence alone does not appear to impact colony size. Nosema apis was found at three sites along the coast and one interior site. Only a small number of pesticides at low concentrations were found. Of the seven common US/European honey bee viruses, only three were identified but, like Varroa, were absent from northern Kenya. The number of viruses present was positively correlated with Varroa levels, but was not correlated with colony size or hygienic behavior. Our results suggest that Varroa, the three viruses, and Nosema have been relatively recently introduced into Kenya, but these factors do not yet appear to be impacting Kenyan bee populations. Thus chemical control for Varroa and Nosema are not necessary for Kenyan bees at this time. This study provides baseline data for future analyses of the possible mechanisms underlying resistance to and the long-term impacts of these factors on African bee populations.     

Crop Pollination Exposes Honey Bees to Pesticides Which Alters Their Susceptibility to the Gut Pathogen Nosema ceranae

Recent declines in honey bee populations and increasing demand for insect-pollinated crops raise concerns about pollinator shortages. Pesticide exposure and pathogens may interact to have strong negative effects on managed honey bee colonies. Such findings are of great concern given the large numbers and high levels of pesticides found in honey bee colonies. Thus it is crucial to determine how field-relevant combinations and loads of pesticides affect bee health. We collected pollen from bee hives in seven major crops to determine 1) what types of pesticides bees are exposed to when rented for pollination of various crops and 2) how field-relevant pesticide blends affect bees’ susceptibility to the gut parasite Nosema ceranae. Our samples represent pollen collected by foragers for use by the colony, and do not necessarily indicate foragers’ roles as pollinators. In blueberry, cranberry, cucumber, pumpkin and watermelon bees collected pollen almost exclusively from weeds and wildflowers during our sampling. Thus more attention must be paid to how honey bees are exposed to pesticides outside of the field in which they are placed. We detected 35 different pesticides in the sampled pollen, and found high fungicide loads. The insecticides esfenvalerate and phosmet were at a concentration higher than their median lethal dose in at least one pollen sample. While fungicides are typically seen as fairly safe for honey bees, we found an increased probability of Nosema infection in bees that consumed pollen with a higher fungicide load. Our results highlight a need for research on sub-lethal effects of fungicides and other chemicals that bees placed in an agricultural setting are exposed to.     

Thinking beyond Western commercial honeybee hives: towards improved conservation of honey bee diversity

A decline of wild pollinators, along with a decline of bee diversity, has been a cause of concern among academics and governmental organizations. According to IPBES, a lack of wild pollinator data contributes to difficulties in comprehensively analyzing the regional status of wild pollinators in Africa, Latin America, Asia and Oceania. It may have also contributed to the prevailing lack of awareness of the diversity of honey bees, of which the managed Apis mellifera is often considered as “the (only) honey bee,” despite the fact that there are eight other honey bee species extant in Asia. A survey of 100 journal articles published in 2016 shows that 57% of the studies still identified A. mellifera as “the honey bee.” In total, 80% of studies were conducted solely on A. mellifera. This focus on A. mellifera has also caused the honey standard of Codex Alimentarius and the European Union to be based solely on A. mellifera, causing improper evaluation of honeys from other species. We recommend adapting current standards to reflect the diversity of honey bees and in the process correct failures in the honey market and pave the way towards improved protection of honey bee species and their habitats.     

Admixture and diversity in West African cattle populations

We present a population genetic analysis of microsatellite variation in 16 West African cattle populations. West Africa represents a unique juxtaposition of different climatic and ecological zones in a relatively small geographical area. While more humid coastal regions are inhabited by the tsetse fly, a vector which spreads trypanosomiasis among cattle, the disease is not transmitted in the drier areas outside this zone. This is the most thorough study of genetic diversity in cattle within this area, which contains genetically important trypanotolerant Bos taurus breeds. Genetic relationships among the many breeds are examined and levels of diversity are assessed. Admixture levels were determined using a variety of methods. Ancestry informative or population-associated alleles (PAAs) were selected using populations from India, the Near East and Europe. Multivariate analysis, the admix program and model-based Bayesian admixture analysis approaches were also employed. These analyses reveal the direct impact of ecological factors and the profound effect of admixture on the cattle of this region. They also highlight the importance of efforts to prevent further dilution of African taurine breeds by B. indicus cattle.     

Honey bee‐collected pollen in agro‐ecosystems reveals diet diversity,diet quality,and pesticide exposure

European honey bees Apis mellifera are important commercial pollinators that have suffered greater than normal overwintering losses since 2007 in North America and Europe. Contributing factors likely include a combination of parasites, pesticides, and poor nutrition. We examined diet diversity, diet nutritional quality, and pesticides in honey bee‐collected pollen from commercial colonies in the Canadian Maritime Provinces in spring and summer 2011. We sampled pollen collected by honey bees at colonies in four site types: apple orchards, blueberry fields, cranberry bogs, and fallow fields. Proportion of honey bee‐collected pollen from crop versus noncrop flowers was high in apple, very low in blueberry, and low in cranberry sites. Pollen nutritional value tended to be relatively good from apple and cranberry sites and poor from blueberry and fallow sites. Floral surveys ranked, from highest to lowest in diversity, fallow, cranberry, apple, and blueberry sites. Pesticide diversity in honey bee‐collected pollen was high from apple and blueberry sites and low from cranberry and fallow sites. Four different neonicotinoid pesticides were detected, but neither these nor any other pesticides were at or above LD₅₀ levels. Pollen hazard quotients were highest in apple and blueberry sites and lowest in fallow sites. Pollen hazard quotients were also negatively correlated with the number of flower taxa detected in surveys. Results reveal differences among site types in diet diversity, diet quality, and pesticide exposure that are informative for improving honey bee and land agro‐ecosystem management.     

Higher vitellogenin concentrations in honey bee workers may be an adaptation to life in temperate climates

The honey bee originated in tropical Africa and later dispersed to northern Europe. It has been suggested that a higher hemolymph storage capacity for the glycolipoprotein vitellogenin evolved in temperate regions, and that the trait constitutes an adaptation to a strongly seasonal environment. We have investigated whether the relative vitellogenin levels of European and African honey bees are in accordance with this hypothesis. Our data indicate that European workers have a higher set-point concentration for vitellogenin compared to their African origin. Considered together with available life history information and physiological data, the results lend support to the view that “winter bees”, a longlived honey bee worker caste that survives winter in temperate regions, evolved through an increase in the worker bees’ capacity for vitellogenin accumulation. Received 20 September 2004; revised 25 March 2005; accepted 13 April 2005.     

Viral transmission in honey bees and native bees,supported by a global black queen cell virus phylogeny

In recent decades, we have realized that honey bee viruses are not, in fact, exclusive to honey bees. The potential impact of Apis-affiliated viruses on native pollinators is prompting concern. Our research addresses the issue of virus crossover between honey bees and native bees foraging in the same localities. We measured the presence of black queen cell virus (BQCV), deformed wing virus (DWV) and sacbrood virus (SBV) in managed Apis mellifera (honey bees) and native Andrena spp. (subgenus Melandrena) bee populations in five commercial orchards. We identified viral presence across sites and bees and related these data to measures of bee community diversity. All viruses were found in both managed and native bees, and BQCV was the most common virus in each. To establish evidence for viral crossover between taxa, we undertook an additional examination of BQCV where 74 samples were sequenced and placed in a global phylogenic framework of hundreds of BQCV strains. We demonstrate pathogen sharing across managed honey bees and distantly related wild bees. This phylogenetic analysis contributes to growing evidence for host switching and places local incidence patterns in a worldwide context, revealing multispecies viral transmission.     

Genetic network analysis between Apis mellifera subspecies based on mtDNA argues the purity of specimens from North Africa,the Levant and Saudi Arabia

ObjectivesThis study aimed to analyze the genetic relationships between honey bee subspecies using reference specimens and recently collected specimens from different parts of the world. The purity of these specimens was discussed in light of the obtained results.MethodsThe genetic networks were constructed between 21 subspecies of honey bees, Apis mellifera L.: 9 in Africa, 7 in Europe and 5 in Asia. The analysis was performed using the mtDNA of these subspecies and the Population Analysis with Reticulate Trees software. Some subspecies were represented by more than two specimens based on the available online sequences.Results and conclusionsThe subspecies A. m. sahariensis from Africa showed unique characteristics and is genetically isolated than all other studied bee subspecies. Specimens collected from Saudi Arabia showed genetic relatedness to A. m. jemenitica, A. m. lamarckii, and some European subspecies, suggesting high degree of hybridization. The close genetic relationship between the Egyptian bees, A. m. lamarckii, and the Syrian bees, A. m. syriaca, were emphasized. The overall genetic network showed the presence of three distinct branches in relation to geographical locations. The high accurateness of the used analysis was confirmed by previous phylogenetic studies as well as the genetic relationships between hybrid bees of A. m. capensis and A. m. scutellata. The genetic networks showed the presence of bee subspecies from Africa in all branches including Europe and Asia. The study suggests the impurity of some specimens mostly due to the hybridization between subspecies. Specific recommendations for future conservation efforts of bees were presented in light of this study.     

Geographical segregation in Dunlin Calidris alpina populations wintering along the East Atlantic migratory flyway – evidence from mitochondrial DNA analysis

Dunlin Calidris alpina is one of the most abundant shorebirds using coastal habitats in the East Atlantic migratory flyway, that links arctic breeding locations (Greenland to Siberia) with wintering grounds (West Europe to West Africa). Differential migration and winter segregation between populations have been indicated by morphometrics and ringing recoveries. Here, we analyse the potential of genetic markers (mitochondrial DNA – mtDNA) to validate and enhance such findings. We compared mtDNA haplotypes frequencies at different wintering sites (from north-west Europe to West Africa). All birds from West Africa had western (European) haplotypes, while the eastern (Siberian) haplotypes were only present in European winter samples, reaching higher frequencies further north in Europe. Compilation of published results from migrating birds also confirmed these differences, with the sole presence of European haplotypes in Iberia and West Africa and increasingly higher frequencies of Siberian haplotypes from south-west to north-west Europe. Comparison with published haplotype frequencies of breeding populations shows that birds from Greenland, Iceland, and North Europe were predominant in wintering grounds in West Africa, while populations wintering in West Europe originated from more eastern breeding grounds (e.g. North Russia). These results show that genetic markers can be used to enhance the integrative monitoring of wintering and breeding populations, by providing biogeographical evidence that validate the winter segregation of breeding populations.     

Colonization history and population differentiation of the Honey Bees (Apis mellifera L.) in Puerto Rico

Honey bees (Apis mellifera L.) are the primary commercial pollinators across the world. The subspecies A. m. scutellata originated in Africa and was introduced to the Americas in 1956. For the last 60 years, it hybridized successfully with European subspecies, previous residents in the area. The result of this hybridization was called Africanized honey bee (AHB). AHB has spread since then, arriving to Puerto Rico (PR) in 1994. The honey bee population on the island acquired a mosaic of features from AHB or the European honey bee (EHB). AHB in Puerto Rico shows a major distinctive characteristic, docile behavior, and is called gentle Africanized honey bees (gAHB). We used 917 SNPs to examine the population structure, genetic differentiation, origin, and history of range expansion and colonization of gAHB in PR. We compared gAHB to populations that span the current distribution of A. mellifera worldwide. The gAHB population is shown to be a single population that differs genetically from the examined populations of AHB. Texas and PR groups are the closest genetically. Our results support the hypothesis that the Texas AHB population is the source of gAHB in Puerto Rico.     

Reduced SNP Panels for Genetic Identification and Introgression Analysis in the Dark Honey Bee (Apis mellifera mellifera)

Beekeeping activities, especially queen trading, have shaped the distribution of honey bee (Apis mellifera) subspecies in Europe, and have resulted in extensive introductions of two eastern European C-lineage subspecies (A. m. ligustica and A. m. carnica) into the native range of the M-lineage A. m. mellifera subspecies in Western Europe. As a consequence, replacement and gene flow between native and commercial populations have occurred at varying levels across western European populations. Genetic identification and introgression analysis using molecular markers is an important tool for management and conservation of honey bee subspecies. Previous studies have monitored introgression by using microsatellite, PCR-RFLP markers and most recently, high density assays using single nucleotide polymorphism (SNP) markers. While the latter are almost prohibitively expensive, the information gained to date can be exploited to create a reduced panel containing the most ancestry-informative markers (AIMs) for those purposes with very little loss of information. The objective of this study was to design reduced panels of AIMs to verify the origin of A. m. mellifera individuals and to provide accurate estimates of the level of C-lineage introgression into their genome. The discriminant power of the SNPs using a variety of metrics and approaches including the Weir & Cockerham’s F_ST, an F_ST-based outlier test, Delta, informativeness (I_n), and PCA was evaluated. This study shows that reduced AIMs panels assign individuals to the correct origin and calculates the admixture level with a high degree of accuracy. These panels provide an essential tool in Europe for genetic stock identification and estimation of admixture levels which can assist management strategies and monitor honey bee conservation programs.     

Effects of Imidacloprid and Varroa destructor on survival and health of European honey bees,Apis mellifera

There has been growing concern over declines in populations of honey bees and other pollinators which are a vital part to our food security. It is imperative to identify factors responsible for accelerated declines in bee populations and develop solutions for reversing bee losses. While exact causes of colony losses remain elusive, risk factors thought to play key roles are ectoparasitic mites Varroa destructor and neonicotinoid pesticides. The present study aims to investigate effects of a neonicotinoid pesticide Imidacloprid and Varroa mites individually on survivorship, growth, physiology, virus dynamics and immunity of honey bee workers. Our study provides clear evidence that the exposure to sublethal doses of Imidacloprid could exert a significantly negative effect on health and survival of honey bees. We observed a significant reduction in the titer of vitellogenin (Vg), an egg yolk precursor that regulates the honey bees development and behavior and often are linked to energy homeostasis, in bees exposed to Imidacloprid. This result indicates that sublethal exposure to neonicotinoid could lead to increased energy usage in honey bees as detoxification is a energy‐consuming metabolic process and suggests that Vg could be a useful biomarker for measuring levels of energy stress and sublethal effects of pesticides on honey bees. Measurement of the quantitative effects of different levels of Varroa mite infestation on the replication dynamic of Deformed wing virus (DWV), an RNA virus associated with Varroa infestation, and expression level of immune genes yields unique insights into how honey bees respond to stressors under laboratory conditions.