HBeeID: a molecular tool that identifies honey bee subspecies from different geographic populations.

BACKGROUND: Honey bees are the principal commercial pollinators. Along with other arthropods, they are increasingly under threat from anthropogenic factors such as the incursion of invasive honey bee subspecies, pathogens and parasites. Better tools are needed to identify bee subspecies. Genomic data for economic and ecologically important organisms is increasing, but in its basic form its practical application to address ecological problems is limited. RESULTS: We introduce HBeeID a means to identify honey bees. The tool utilizes a knowledge-based network and diagnostic SNPs identified by discriminant analysis of principle components and hierarchical agglomerative clustering. Tests of HBeeID showed that it identifies African, Americas-Africanized, Asian, and European honey bees with a high degree of certainty even when samples lack the full 272 SNPs of HBeeID. Its prediction capacity decreases with highly admixed samples. CONCLUSION: HBeeID is a high-resolution genomic, SNP based tool, that can be used to identify honey bees and screen species that are invasive. Its flexible design allows for future improvements via sample data additions from other localities.

Identification of Multiple Loci Associated with Social Parasitism in Honeybees.

In colonies of the honeybee Apis mellifera, the queen is usually the only reproductive female, which produces new females (queens and workers) by laying fertilized eggs. However, in one subspecies of A. mellifera, known as the Cape bee (A. m. capensis), worker bees reproduce asexually by thelytoky, an abnormal form of meiosis where two daughter nucleii fuse to form single diploid eggs, which develop into females without being fertilized. The Cape bee also exhibits a suite of phenotypes that facilitate social parasitism whereby workers lay such eggs in foreign colonies so their offspring can exploit their resources. The genetic basis of this switch to social parasitism in the Cape bee is unknown. To address this, we compared genome variation in a sample of Cape bees with other African populations. We find genetic divergence between these populations to be very low on average but identify several regions of the genome with extreme differentiation. The regions are strongly enriched for signals of selection in Cape bees, indicating that increased levels of positive selection have produced the unique set of derived phenotypic traits in this subspecies. Genetic variation within these regions allows unambiguous genetic identification of Cape bees and likely underlies the genetic basis of social parasitism. The candidate loci include genes involved in ecdysteroid signaling and juvenile hormone and dopamine biosynthesis, which may regulate worker ovary activation and others whose products localize at the centrosome and are implicated in chromosomal segregation during meiosis. Functional analysis of these loci will yield insights into the processes of reproduction and chemical signaling in both parasitic and non-parasitic populations and advance understanding of the process of normal and atypical meiosis.

Differences in Varroa destructor infestation rates of two indigenous subspecies of Apis mellifera in the Republic of South Africa.

Varroa destructor Anderson & Trueman (Varroa) is a damaging pest of the Western honey bee, Apis mellifera, in North America, Europe, and Asia. However, Varroa infestations have not produced equivalent colony losses of African subspecies of honey bee throughout Africa and parts of the Americas. We surveyed the Varroa infestation rates (number of Varroa per 100 adult honey bees) in colonies of A. m. scutellata, A. m. capensis, and hybrids of the two subspecies throughout the Republic of South Africa in the fall of 2014. We found that A. m. scutellata colonies had significantly higher Varroa infestations than did A. m. capensis colonies. Furthermore, hybridized colonies of the two subspecies had Varroa infestations intermediate to those of A. m. scutellata and A. m. capensis. This is the first documentation of a clear difference in Varroa infestation rates of A. m. scutellata, A. m. capensis, and hybridized colonies in South Africa. Furthermore, our data confirm that Varroa populations in A. m. scutellata colonies are within the range of populations that are damaging to European honey bees.

Exploring the impact of fungicide exposure and nutritional stress on the microbiota and immune response of the Cape honey bee (Apis mellifera capensis).

Honey bees (Apis mellifera) harbour a stable core microbial community within their gut, that is suggested to play a role in metabolic functioning, immune regulation, and host homeostasis. This microbiota presents a unique opportunity to observe the effects of stressors on honey bee health. We examined the effects of two common honey bee stressors: indirect fungicide contamination and nutrient limitation. These effects were observed through changes in their hind- and midgut microbiota using Automated Ribosomal Intergenic Spacer Analysis (ARISA), alongside high-throughput amplicon sequencing. Expression of the honey bees' immune response was examined through the expression of three immune-related genes, namely, immune deficiency (imd), proPhenolOxidase (proPO), and spaetzle (spz). Additionally, longevity of the honey bees was monitored through observation of the expression levels of Vitellogenin (Vg). Both treatment groups were compared to a negative control, and a diseased positive control. There was no effect on the hindgut microbiota due to the stressors, while significant changes in the midgut was observed. This was also observed in the expression of the immune-related genes within the treatment groups. The Imd pathway was substantially downregulated, with upregulation in the prophenoloxidase pathway. However, no significant effect was observed in the expression of spz, and only the pollen treatment group showed reduced longevity through a downregulation of Vg. Overall, the effect of these two common stressors indicate a compromise in honey bee immunity, and potential vulnerabilities within the immune defence mechanisms.

Valuation of pollinator forage services provided by Eucalyptus cladocalyx.

We assess the monetary value of forage provisioning services for honeybees as provided by an alien tree species in the Western Cape province of South Africa. Although Eucalyptus cladocalyx is not an officially declared invader, it is cleared on a regular basis along with other invasive Eucalyptus species such as Eucalyptus camaldulensis, and Eucalyptus conferruminata (which have been prioritised for eradication in South Africa). We present some of the trade-offs associated with the clearing of E. cladocalyx by means of a practical example that illustrates a situation where the benefits of the species to certain stakeholders could support the containment of the species in demarcated areas, while allowing clearing outside such areas. Given the absence of market prices for such forage provisioning services, the replacement cost is used to present the value of the loss in forage as provided by E. cladocalyx if the alien tree species is cleared along with invasive alien tree species. Two replacement scenarios formed the basis for our calculations. The first scenario was an artificial diet as replacement for the forage provisioning service, which yielded a direct cost estimate of US$7.5 m per year. The second was based on a Fynbos cultivation/restoration initiative aimed at substituting the forage provisioning service of E. cladocalyx, which yielded a direct cost of US$20.2 m per year. These figures provide estimates of the potential additional cost burden on the beekeeping industry if E. cladocalyx is completely eradicated from the Western Cape. The cost estimates should be balanced against the negative impacts of E. cladocalyx on ecosystem services in order to make an informed decision with regard to appropriate management strategies for this species. The findings therefore serve as useful inputs to balance trade-offs for alien species that are considered as beneficial to some, but harmful to other.

Application of ANOVA-simultaneous component analysis to quantify and characterise effects of age, temperature, syrup adulteration and irradiation on near-infrared (NIR) spectral data of honey.

NIR spectroscopy combined with chemometric analysis has proven to be a rapid and cost-effective screening tool for the detection of syrup-adulterated honey. Processing and storage conditions which alter the chemical and physical state of honey may affect the spectra. The effects of age, storage temperature, syrup adulteration (10 and 20% w/w) and irradiation treatment on the NIR spectra of honey were investigated as a function of time with ANOVA-simultaneous component analysis (ASCA), an experimental design-focused exploratory data analysis method. The factors 'time', 'temperature' and 'adulteration' were found to have significant effects (p < 0.05), but no significant effect was observed for irradiation treatment. A significant interaction effect was found between factors 'time' and 'adulteration', with the greatest disparity between authentic and adulterated class signals found immediately after adulteration and decreasing within three months thereafter.

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.

Substantial Heritable Variation in Recombination Rate on Multiple Scales in Honeybees and Bumblebees.

Meiotic recombination shuffles genetic variation and promotes correct segregation of chromosomes. Rates of recombination vary on several scales, both within genomes and between individuals, and this variation is affected by both genetic and environmental factors. Social insects have extremely high rates of recombination, although the evolutionary causes of this are not known. Here, we estimate rates of crossovers and gene conversions in 22 colonies of the honeybee, Apis mellifera, and 9 colonies of the bumblebee, Bombus terrestris, using direct sequencing of 299 haploid drone offspring. We confirm that both species have extremely elevated crossover rates, with higher rates measured in the highly eusocial honeybee than the primitively social bumblebee. There are also significant differences in recombination rate between subspecies of honeybee. There is substantial variation in genome-wide recombination rate between individuals of both A. mellifera and B. terrestris and the distribution of these rates overlap between species. A large proportion of interindividual variation in recombination rate is heritable, which indicates the presence of variation in trans-acting factors that influence recombination genome-wide. We infer that levels of crossover interference are significantly lower in honeybees compared to bumblebees, which may be one mechanism that contributes to higher recombination rates in honeybees. We also find a significant increase in recombination rate with distance from the centromere, mirrored by methylation differences. We detect a strong transmission bias due to GC-biased gene conversion associated with noncrossover gene conversions. Our results shed light on the mechanistic causes of extreme rates of recombination in social insects and the genetic architecture of recombination rate variation.

Whole-genome scan in thelytokous-laying workers of the Cape honeybee (Apis mellifera capensis): central fusion, reduced recombination rates and centromere mapping using half-tetrad analysis.

While workers of almost all subspecies of honeybee are able to lay only haploid male eggs, Apis mellifera capensis workers are able to produce diploid female eggs by thelytokous parthenogenesis. Cytological analyses have shown that during parthenogenesis, egg diploidy is restored by fusion of the two central meiotic products. This peculiarity of the Cape bee preserves two products of a single meiosis in the daughters and can be used to map centromere positions using half-tetrad analysis. In this study, we use the thelytokous progenies of A. m. capensis workers and a sample of individuals from a naturally occurring A. m. capensis thelytokous clone to map centromere position for most of the linkage groups of the honeybee. We also show that the recombination rate is reduced by >10-fold during the meiosis of A. m. capensis workers. This reduction is restricted to thelytokous parthenogenesis of capensis workers and is not observed in the meiosis of queen within the same subspecies or in arrhenotokous workers of another subspecies. The reduced rate of recombination seems to be associated with negative crossover interference. These results are discussed in relation to evolution of thelytokous parthenogenesis and maintenance of heterozygosity and female sex after thelytoky.

A worldwide survey of genome sequence variation provides insight into the evolutionary history of the honeybee Apis mellifera.

The honeybee Apis mellifera has major ecological and economic importance. We analyze patterns of genetic variation at 8.3 million SNPs, identified by sequencing 140 honeybee genomes from a worldwide sample of 14 populations at a combined total depth of 634×. These data provide insight into the evolutionary history and genetic basis of local adaptation in this species. We find evidence that population sizes have fluctuated greatly, mirroring historical fluctuations in climate, although contemporary populations have high genetic diversity, indicating the absence of domestication bottlenecks. Levels of genetic variation are strongly shaped by natural selection and are highly correlated with patterns of gene expression and DNA methylation. We identify genomic signatures of local adaptation, which are enriched in genes expressed in workers and in immune system- and sperm motility-related genes that might underlie geographic variation in reproduction, dispersal and disease resistance. This study provides a framework for future investigations into responses to pathogens and climate change in honeybees.

Valuing insect pollination services with cost of replacement.

Value estimates of ecosystem goods and services are useful to justify the allocation of resources towards conservation, but inconclusive estimates risk unsustainable resource allocations. Here we present replacement costs as a more accurate value estimate of insect pollination as an ecosystem service, although this method could also be applied to other services. The importance of insect pollination to agriculture is unequivocal. However, whether this service is largely provided by wild pollinators (genuine ecosystem service) or managed pollinators (commercial service), and which of these requires immediate action amidst reports of pollinator decline, remains contested. If crop pollination is used to argue for biodiversity conservation, clear distinction should be made between values of managed- and wild pollination services. Current methods either under-estimate or over-estimate the pollination service value, and make use of criticised general insect and managed pollinator dependence factors. We apply the theoretical concept of ascribing a value to a service by calculating the cost to replace it, as a novel way of valuing wild and managed pollination services. Adjusted insect and managed pollinator dependence factors were used to estimate the cost of replacing insect- and managed pollination services for the Western Cape deciduous fruit industry of South Africa. Using pollen dusting and hand pollination as suitable replacements, we value pollination services significantly higher than current market prices for commercial pollination, although lower than traditional proportional estimates. The complexity associated with inclusive value estimation of pollination services required several defendable assumptions, but made estimates more inclusive than previous attempts. Consequently this study provides the basis for continued improvement in context specific pollination service value estimates.

Development of infectious transcripts and genome manipulation of Black queen-cell virus of honey bees.

The South African isolate of Black queen-cell virus (BQCV), a honey bee virus, was previously found to have an 8550 nucleotide genome excluding the poly(A) tail. Its genome contained two ORFs, a 5'-proximal ORF encoding a putative replicase protein and a 3'-proximal ORF encoding a capsid polyprotein. Long reverse transcription (RT)-PCR was used to produce infectious transcripts for BQCV and to manipulate its genome. Primers were designed for the amplification of the complete genome, the in vitro transcription of infectious RNA and PCR-directed mutagenesis. An 18-mer antisense primer was designed for RT to produce full-length single-stranded cDNA (ss cDNA). Unpurified ss cDNA from the RT reaction mixture was used directly as a template to amplify the full genome by long high-fidelity PCR. The SP6 promoter sequence was introduced into the sense primer to transcribe RNA directly from the amplicon. RNA was transcribed in vitro with and without the presence of a cap analogue and injected directly into bee pupae, which were then incubated for 8 days. In vitro transcripts were infectious but the presence of a cap analogue did not increase the amount of virus recovered. A single base mutation abolishing an EcoRI restriction site was introduced by fusion-PCR, to distinguish viral particles recovered from infectious transcripts from wild-type virus (wtBQCV). Mutant virus (mutBQCV) and wtBQCV were indistinguishable by electron microscopy and Western blot analysis. The EcoRI restriction site was present in wtBQCV and not in mutBQCV.

A non-policing honey bee colony (Apis mellifera capensis).

In the Cape honey bee Apis mellifera capensis, workers lay female eggs without mating by thelytokous parthenogenesis. As a result, workers are as related to worker-laid eggs as they are to queen-laid eggs and therefore worker policing is expected to be lower, or even absent. This was tested by transferring worker- and queen-laid eggs into three queenright A. m. capensis discriminator colonies and monitoring their removal. Our results show that worker policing is variable in A. m. capensis and that in one colony worker-laid eggs were not removed. This is the first report of a non-policing queenright honey bee colony. DNA microsatellite and morphometric analysis suggests that the racial composition of the three discriminator colonies was different. The variation in policing rates could be explained by differences in degrees of hybridisation between A. m. capensis and A. m. scutellata, although a larger survey is needed to confirm this.