Recent advances in population and quantitative genomics of honey bees

The molecular determinants of pesticide sensitivity in bee pollinators

Bees carry out vital ecosystem services by pollinating both wild and economically important crop plants. However, while performing this function, bee pollinators may encounter potentially harmful xenobiotics in the environment such as pesticides (fungicides, herbicides and insecticides). Understanding the key factors that influence the toxicological outcomes of bee exposure to these chemicals, in isolation or combination, is essential to safeguard their health and the ecosystem services they provide. In this regard, recent work using toxicogenomic and phylogenetic approaches has begun to identify, at the molecular level, key determinants of pesticide sensitivity in bee pollinators. These include detoxification systems that convert pesticides to less toxic forms and key residues in insecticide target-sites that underlie species-specific insecticide selectivity. Here we review this emerging body of research and summarise the state of knowledge of the molecular determinants of pesticide sensitivity in bee pollinators. We identify gaps in our knowledge for future research and examine how an understanding of the genetic basis of bee sensitivity to pesticides can be leveraged to, a) predict and avoid negative bee-pesticide interactions and facilitate the future development of pest-selective bee-safe insecticides, and b) inform traditional effect assessment approaches in bee pesticide risk assessment and address issues of ecotoxicological concern.

Genetic differentiation and local adaptation of the Japanese honeybee, Apis cerana japonica

We examine the population genetic structure and divergence among the regional populations of the Japanese honeybee, Apis cerana japonica, by re-sequencing the genomes of 105 individuals from the three main Japanese islands with diverse climates. The genetic structure results indicated that these individuals are distinct from the mainland Chinese A. cerana samples. Furthermore, population structure analyses have identified three genetically distinct geographic regions in Japan: Northern (Tohoku-Kanto-Chubu districts), Central (Chugoku district), and Southern (Kyushu district). In some districts, "possible non-native" individuals, likely introduced from other regions in recent years, were discovered. Then, genome-wide scans were conducted to detect candidate genes for adaptation by two different approaches. We performed a population branch statistics (PBS) analysis to identify candidate genes for population-specific divergence. A latent factor mixed model (LFMM) was used to identify genes associated with climatic variables along a geographic gradient. The PBSmax analysis identified 25 candidate genes for population-specific divergence whereas the LFMM analysis identified 73 candidate genes for adaptation to climatic variables along a geographic gradient. However, no common genes were identified by both methods.

Genetics of tolerance in honeybees to the neonicotinoid clothianidin

The effects of neonicotinoid insecticides (NNIs) on honeybee health are intensely debated, with numerous studies showing negative effects of exposure, while others report no such effects. We carried out experiments to study the genetic and molecular basis of NNI tolerance in honeybees, which may underlie the discrepancies observed in the literature. We discovered that worker survival post-exposure to an acute oral dose of clothianidin is heritable (H ² = 37.8%). Tolerance to clothianidin was not associated with differences in the expression of detoxification enzymes in our experiments. Instead, mutations in the primary neonicotinoid detoxification genes CYP9Q1 and CYP9Q3 were strongly associated with worker survival post-clothianidin exposure. In some instances, the strong association between CYP9Q haplotypes and worker survival was associated with the protein's predicted binding affinity for clothianidin. Our findings have implications regarding future toxicological studies utilizing honeybees as a model pollinator.

A genotyping by sequencing approach can disclose Apis mellifera population genomic information contained in honey environmental DNA

Awareness has been raised over the last years on the genetic integrity of autochthonous honey bee subspecies. Genomic tools available in Apis mellifera can make it possible to measure this information by targeting individual honey bee DNA. Honey contains DNA traces from all organisms that contributed or were involved in its production steps, including the honey bees of the colony. In this study, we designed and tested a genotyping by sequencing (GBS) assay to analyse single nucleotide polymorphisms (SNPs) of A. mellifera nuclear genome using environmental DNA extracted from honey. A total of 121 SNPs (97 SNPs informative for honey bee subspecies identification and 24 SNPs associated with relevant traits of the colonies) were used in the assay to genotype honey DNA, which derives from thousands of honey bees. Results were integrated with information derived from previous studies and whole genome resequencing datasets. This GBS method is highly reliable in estimating honey bee SNP allele frequencies of the whole colony from which the honey derived. This assay can be used to identify the honey bee subspecies of the colony that produced the honey and, in turn, to authenticate the entomological origin of the honey.

Thrice out of Asia and the adaptive radiation of the western honey bee

The origin of the western honey bee Apis mellifera has been intensely debated. Addressing this knowledge gap is essential for understanding the evolution and genetics of one of the world’s most important pollinators. By analyzing 251 genomes from 18 native subspecies, we found support for an Asian origin of honey bees with at least three expansions leading to African and European lineages. The adaptive radiation of honey bees involved selection on a few genomic “hotspots.” We found 145 genes with independent signatures of selection across all bee lineages, and these genes were highly associated with worker traits. Our results indicate that a core set of genes associated with worker and colony traits facilitated the adaptive radiation of honey bees across their vast distribution.

Tissue-specific transcriptional patterns underlie seasonal phenotypes in honey bees (Apis mellifera)

Faced with adverse conditions, such as winter in temperate regions or hot and dry conditions in tropical regions, many insect species enter a state of diapause, a period of dormancy associated with a reduction or arrest of physical activity, development and reproduction. Changes in common physiological pathways underlie diapause phenotypes in different insect species. However, most transcriptomic studies of diapause have not simultaneously evaluated and compared expression patterns in different tissues. Honey bees (Apis mellifera) represent a unique model system to study the mechanisms underpinning diapause-related phenotypes. In winter, honey bees exhibit a classic diapause phenotype, with reduced metabolic activity, increased physiological nutritional resources and altered hormonal profiles. However, winter bees actively heat their colony by vibrating their wing muscles; thus, this tissue is not quiescent. Here, we evaluated the transcriptional profiles of flight muscle tissue and fat body tissue (involved in nutrient storage, metabolism and immune function) of winter bees. We also evaluated two behavioural phenotypes of summer bees: nurses, which exhibit high nutritional stores and low flight activity, and foragers, which exhibit low nutritional stores and high flight activity. We found winter bees and nurses have similar fat body transcriptional profiles, whereas winter bees and foragers have similar flight muscle transcriptional profiles. Additionally, differentially expressed genes were enriched in diapause-related gene ontology terms. Thus, honey bees exhibit tissue-specific transcriptional profiles associated with seasonal phenotypes, laying the groundwork for future studies evaluating the mechanisms, evolution and consequences of this tissue-specific regulation.

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

Objectives

This 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.

Methods

The 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 conclusions

The 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.

Disentangling Ethiopian Honey Bee (Apis mellifera) Populations Based on Standard Morphometric and Genetic Analyses

Simple Summary

We conducted this population study of Ethiopian honey bees, using morphometric and genetic methods, to decipher their controversial classification. These honey bees are highly diverse and showed differentiation based on size and genetic information according to prevailing agro-ecological conditions, demonstrating morphological and molecular signatures of local adaptation. The results of both morphometric and genetic analyses suggest that Ethiopian honey bees differ from populations in the neighboring geographic regions and are characterized by extensive gene flow within the country, enhanced by honey bee colony trade. Consequently, future research that includes studying traits of vitality, behavior and colony performance of honey bees in remaining pocket areas of highland agro-ecological zones could contribute to the development of appropriate conservation management.

Abstract

The diversity and local differentiation of honey bees are subjects of broad general interest. In particular, the classification of Ethiopian honey bees has been a subject of debate for decades. Here, we conducted an integrated analysis based on classical morphometrics and a putative nuclear marker (denoted r7-frag) for elevational adaptation to classify and characterize these honey bees. Therefore, 660 worker bees were collected out of 66 colonies from highland, midland and lowland agro-ecological zones (AEZs) and were analyzed in reference to populations from neighboring countries. Multivariate morphometric analyses show that our Ethiopian samples are separate from Apis mellifera scutellata, A. m. jemenitica, A. m. litorea and A. m. monticola, but are closely related to A. m. simensis reference. Linear discriminant analysis showed differentiation according to AEZs in the form of highland, midland and lowland ecotypes. Moreover, size was positively correlated with elevation. Similarly, our Ethiopian samples were differentiated from A. m. monticola and A. m. scutellata based on r7-frag. There was a low tendency towards genetic differentiation between the Ethiopian samples, likely impacted by increased gene flow. However, the differentiation slightly increased with increasing elevational differences, demonstrated by the highland bees that showed higher differentiation from the lowland bees (F_ST = 0.024) compared to the midland bees (F_ST = 0.015). An allelic length polymorphism was detected (denoted as d) within r7-frag, showing a patterned distribution strongly associated with AEZ (X² = 11.84, p < 0.01) and found predominantly in highland and midland bees of some pocket areas. In conclusion, the Ethiopian honey bees represented in this study are characterized by high gene flow that suppresses differentiation.

Authoritative subspecies diagnosis tool for European honey bees based on ancestry informative SNPs

BACKGROUND

With numerous endemic subspecies representing four of its five evolutionary lineages, Europe holds a large fraction of Apis mellifera genetic diversity. This diversity and the natural distribution range have been altered by anthropogenic factors. The conservation of this natural heritage relies on the availability of accurate tools for subspecies diagnosis. Based on pool-sequence data from 2145 worker bees representing 22 populations sampled across Europe, we employed two highly discriminative approaches (PCA and F_ST) to select the most informative SNPs for ancestry inference.

RESULTS

Using a supervised machine learning (ML) approach and a set of 3896 genotyped individuals, we could show that the 4094 selected single nucleotide polymorphisms (SNPs) provide an accurate prediction of ancestry inference in European honey bees. The best ML model was Linear Support Vector Classifier (Linear SVC) which correctly assigned most individuals to one of the 14 subspecies or different genetic origins with a mean accuracy of 96.2% ± 0.8 SD. A total of 3.8% of test individuals were misclassified, most probably due to limited differentiation between the subspecies caused by close geographical proximity, or human interference of genetic integrity of reference subspecies, or a combination thereof.

CONCLUSIONS

The diagnostic tool presented here will contribute to a sustainable conservation and support breeding activities in order to preserve the genetic heritage of European honey bees.

Mitochondrial genomes illuminate the evolutionary history of the Western honey bee (Apis mellifera)

Western honey bees (Apis mellifera) are one of the most important pollinators of agricultural crops and wild plants. Despite the growth in the availability of sequence data for honey bees, the phylogeny of the species remains a subject of controversy. Most notably, the geographic origin of honey bees is uncertain, as are the relationships among its constituent lineages and subspecies. We aim to infer the evolutionary and biogeographical history of the honey bee from mitochondrial genomes. Here we analyse the full mitochondrial genomes of 18 A. mellifera subspecies, belonging to all major lineages, using a range of gene sampling strategies and inference models to identify factors that may have contributed to the recovery of incongruent results in previous studies. Our analyses support a northern African or Middle Eastern origin of A. mellifera. We show that the previously suggested European and Afrotropical cradles of honey bees are the result of phylogenetic error. Monophyly of the M, C, and O lineages is strongly supported, but the A lineage appears paraphyletic. A. mellifera colonised Europe through at least two pathways, across the Strait of Gibraltar and via Asia Minor.

Gene-regulatory context of honey bee worker sterility

The highly organized societies of the Western honey bee Apis mellifera feature a highly reproductive queen at the center of attention and a large cohort of daughters that suppress their own reproduction to help rear more sisters, some of whom become queens themselves. This reproductive altruism is peculiar because in theory it evolves via indirect selection on genes for altruism that are expressed in the sterile workers but not in the reproductive queens. In this study we attempt to situate lists of genes previously implicated in queenright worker sterility into a broader regulatory framework. To do so we use a model bee brain transcriptional regulatory network as a template to infer how sets of genes responsive to ovary-suppressing queen pheromone are functionally interconnected over the model's topology. We predict that genes jointly involved in the regulation of worker sterility should be tightly networked, relative to genes whose functions are unrelated to each other. We find that sets of mapped genes - ranging in size from 17 to 250 - are well dispersed across the network's substructural scaffolds, suggesting that ovary de-activation involves genes that reside within more than one transcriptional regulatory module. For some sets, however, this dispersion is biased into certain areas of the network's substructure. Our analysis identifies the regions enriched for sterility genes and likewise identifies local hub genes that are presumably critical to subnetwork function. Our work offers a glimpse into the gene regulatory context of honey bee worker sterility and uses this context to identify new candidate gene targets for functional analysis. Finally, to the extent that any sterility-related modules identified here have evolved via selection for worker altruism, we can assume that this selection was indirect and of the type specifically invoked by inclusive fitness theory.

The complete mitochondrial genome of the Xinyuan honey bee, Apis mellifera sinisxinyuan (Insecta: Hymenoptera: Apidae)

We analyzed the complete mitochondrial genome of the recently discovered Xinyuan honey bee, Apis mellifera sinisxinyuan using single molecule real-time sequencing. The mitochondrial genome of A. m. sinisxinyuan is a circular molecule of 16,886 bp, comprising 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes and a control region rich in A + T. Phylogenetic analysis using 13 protein-coding genes supports a close relationship to another M-lineage honey bee, A. m. mellifera.