Uncovering the novel characteristics of Asian honey bee,Apis cerana,by whole genome sequencing

Background

The honey bee is an important model system for increasing understanding of molecular and neural mechanisms underlying social behaviors relevant to the agricultural industry and basic science. The western honey bee, Apis mellifera, has served as a model species, and its genome sequence has been published. In contrast, the genome of the Asian honey bee, Apis cerana, has not yet been sequenced. A. cerana has been raised in Asian countries for thousands of years and has brought considerable economic benefits to the apicultural industry. A cerana has divergent biological traits compared to A. mellifera and it has played a key role in maintaining biodiversity in eastern and southern Asia. Here we report the first whole genome sequence of A. cerana.

Results

Using de novo assembly methods, we produced a 238 Mbp draft of the A. cerana genome and generated 10,651 genes. A.cerana-specific genes were analyzed to better understand the novel characteristics of this honey bee species. Seventy-two percent of the A. cerana-specific genes had more than one GO term, and 1,696 enzymes were categorized into 125 pathways. Genes involved in chemoreception and immunity were carefully identified and compared to those from other sequenced insect models. These included 10 gustatory receptors, 119 odorant receptors, 10 ionotropic receptors, and 160 immune-related genes.

Conclusions

This first report of the whole genome sequence of A. cerana provides resources for comparative sociogenomics, especially in the field of social insect communication. These important tools will contribute to a better understanding of the complex behaviors and natural biology of the Asian honey bee and to anticipate its future evolutionary trajectory.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-16-1) contains supplementary material, which is available to authorized users.     

Hemolymph proteome changes during worker brood development match the biological divergences between western honey bees (Apis mellifera) and eastern honey bees (Apis cerana)

Background

Hemolymph plays key roles in honey bee molecule transport, immune defense, and in monitoring the physiological condition. There is a lack of knowledge regarding how the proteome achieves these biological missions for both the western and eastern honey bees (Apis mellifera and Apis cerana). A time-resolved proteome was compared using two-dimensional electrophoresis-based proteomics to reveal the mechanistic differences by analysis of hemolymph proteome changes between the worker bees of two bee species during the larval to pupal stages.

Results

The brood body weight of Apis mellifera was significantly heavier than that of Apis cerana at each developmental stage. Significantly, different protein expression patterns and metabolic pathways were observed in 74 proteins (166 spots) that were differentially abundant between the two bee species. The function of hemolymph in energy storage, odor communication, and antioxidation is of equal importance for the western and eastern bees, indicated by the enhanced expression of different protein species. However, stronger expression of protein folding, cytoskeletal and developmental proteins, and more highly activated energy producing pathways in western bees suggests that the different bee species have developed unique strategies to match their specific physiology using hemolymph to deliver nutrients and in immune defense.

Conclusions

Our disparate findings constitute a proof-of-concept of molecular details that the ecologically shaped different physiological conditions of different bee species match with the hemolymph proteome during the brood stage. This also provides a starting point for future research on the specific hemolymph proteins or pathways related to the differential phenotypes or physiology.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-563) contains supplementary material, which is available to authorized users.     

A SNP Based High-Density Linkage Map of Apis cerana Reveals a High Recombination Rate Similar to Apis mellifera

Background

The Eastern honey bee, Apis cerana Fabricius, is distributed in southern and eastern Asia, from India and China to Korea and Japan and southeast to the Moluccas. This species is also widely kept for honey production besides Apis mellifera. Apis cerana is also a model organism for studying social behavior, caste determination, mating biology, sexual selection, and host-parasite interactions. Few resources are available for molecular research in this species, and a linkage map was never constructed. A linkage map is a prerequisite for quantitative trait loci mapping and for analyzing genome structure. We used the Chinese honey bee, Apis cerana cerana to construct the first linkage map in the Eastern honey bee.

Results

F2 workers (N = 103) were genotyped for 126,990 single nucleotide polymorphisms (SNPs). After filtering low quality and those not passing the Mendel test, we obtained 3,000 SNPs, 1,535 of these were informative and used to construct a linkage map. The preliminary map contains 19 linkage groups, we then mapped the 19 linkage groups to 16 chromosomes by comparing the markers to the genome of A. mellfiera. The final map contains 16 linkage groups with a total of 1,535 markers. The total genetic distance is 3,942.7 centimorgans (cM) with the largest linkage group (180 loci) measuring 574.5 cM. Average marker interval for all markers across the 16 linkage groups is 2.6 cM.

Conclusion

We constructed a high density linkage map for A. c. cerana with 1,535 markers. Because the map is based on SNP markers, it will enable easier and faster genotyping assays than randomly amplified polymorphic DNA or microsatellite based maps used in A. mellifera.     

Transcriptome differences in the hypopharyngeal gland between Western Honeybees (Apis mellifera) and Eastern Honeybees (Apis cerana)

Background

Apis mellifera and Apis cerana are two sibling species of Apidae. Apis cerana is adept at collecting sporadic nectar in mountain and forest region and exhibits stiffer hardiness and acarid resistance as a result of natural selection, whereas Apis mellifera has the advantage of producing royal jelly. To identify differentially expressed genes (DEGs) that affect the development of hypopharyngeal gland (HG) and/or the secretion of royal jelly between these two honeybee species, we performed a digital gene expression (DGE) analysis of the HGs of these two species at three developmental stages (newly emerged worker, nurse and forager).

Results

Twelve DGE-tag libraries were constructed and sequenced using the total RNA extracted from the HGs of newly emerged workers, nurses, and foragers of Apis mellifera and Apis cerana. Finally, a total of 1482 genes in Apis mellifera and 1313 in Apis cerana were found to exhibit an expression difference among the three developmental stages. A total of 1417 DEGs were identified between these two species. Of these, 623, 1072, and 462 genes showed an expression difference at the newly emerged worker, nurse, and forager stages, respectively. The nurse stage exhibited the highest number of DEGs between these two species and most of these were found to be up-regulated in Apis mellifera. These results suggest that the higher yield of royal jelly in Apis mellifera may be due to the higher expression level of these DEGs.

Conclusions

In this study, we investigated the DEGs between the HGs of two sibling honeybee species (Apis mellifera and Apis cerana). Our results indicated that the gene expression difference was associated with the difference in the royal jelly yield between these two species. These results provide an important clue for clarifying the mechanisms underlying hypopharyngeal gland development and the production of royal jelly.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-744) contains supplementary material, which is available to authorized users.     

Chronic parasitization by Nosema microsporidia causes global expression changes in core nutritional,metabolic and behavioral pathways in honey bee workers (Apis mellifera)

Background

Chronic infections can profoundly affect the physiology, behavior, fitness and longevity of individuals, and may alter the organization and demography of social groups. Nosema apis and Nosema ceranae are two microsporidian parasites which chronically infect the digestive tract of honey bees (Apis mellifera). These parasites, in addition to other stressors, have been linked to increased mortality of individual workers and colony losses in this key pollinator species. Physiologically, Nosema infection damages midgut tissue, is energetically expensive and alters expression of immune genes in worker honey bees. Infection also accelerates worker transition from nursing to foraging behavior (termed behavioral maturation). Here, using microarrays, we characterized global gene expression patterns in adult worker honey bee midgut and fat body tissue in response to Nosema infection.

Results

Our results indicate that N. apis infection in young workers (1 and 2 days old) disrupts midgut development. At 2 and 7 days post-infection in the fat body tissue, N. apis drives metabolic changes consistent with energetic costs of infection. A final experiment characterizing gene expression in the fat bodies of 14 day old workers parasitized with N. apis and N. ceranae demonstrated that Nosema co-infection specifically alters conserved nutritional, metabolic and hormonal pathways, including the insulin signaling pathway, which is also linked to behavioral maturation in workers. Interestingly, in all experiments, Nosema infection did not appear to significantly regulate overall expression of canonical immune response genes, but infection did alter expression of acute immune response genes identified in a previous study. Comparative analyses suggest that changes in nutritional/metabolic processes precede changes in behavioral maturation and immune processes.

Conclusions

These genome-wide studies of expression patterns can help us disentangle the direct and indirect effects of chronic infection, and understand the molecular pathways that regulate disease symptoms.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-14-799) contains supplementary material, which is available to authorized users.     

Drone and Worker Brood Microclimates Are Regulated Differentially in Honey Bees,Apis mellifera

Background

Honey bee (Apis mellifera) drones and workers show differences in morphology, physiology, and behavior. Because the functions of drones are more related to colony reproduction, and those of workers relate to both survival and reproduction, we hypothesize that the microclimate for worker brood is more precisely regulated than that of drone brood.

Methodology/Principal Findings

We assessed temperature and relative humidity (RH) inside honey bee colonies for both drone and worker brood throughout the three-stage development period, using digital HOBO^® Data Loggers. The major findings of this study are that 1) both drone and worker castes show the highest temperature for eggs, followed by larvae and then pupae; 2) temperature in drones are maintained at higher precision (smaller variance) in drone eggs and larvae, but at a lower precision in pupae than the corresponding stages of workers; 3) RH regulation showed higher variance in drone than workers across all brood stages; and 4) RH regulation seems largely due to regulation by workers, as the contribution from empty honey combs are much smaller compared to that from adult workers.

Conclusions/Significance

We conclude that honey bee colonies maintain both temperature and humidity actively; that the microclimate for sealed drone brood is less precisely regulated than worker brood; and that combs with honey contribute very little to the increase of RH in honey bee colonies. These findings increase our understanding of microclimate regulation in honey bees and may have implications for beekeeping practices.     

Imidacloprid Alters Foraging and Decreases Bee Avoidance of Predators

Concern is growing over the effects of neonicotinoid pesticides, which can impair honey bee cognition. We provide the first demonstration that sublethal concentrations of imidacloprid can harm honey bee decision-making about danger by significantly increasing the probability of a bee visiting a dangerous food source. Apis cerana is a native bee that is an important pollinator of agricultural crops and native plants in Asia. When foraging on nectar containing 40 µg/L (34 ppb) imidacloprid, honey bees (Apis cerana) showed no aversion to a feeder with a hornet predator, and 1.8 fold more bees chose the dangerous feeder as compared to control bees. Control bees exhibited significant predator avoidance. We also give the first evidence that foraging by A. cerana workers can be inhibited by sublethal concentrations of the pesticide, imidacloprid, which is widely used in Asia. Compared to bees collecting uncontaminated nectar, 23% fewer foragers returned to collect the nectar with 40 µg/L imidacloprid. Bees that did return respectively collected 46% and 63% less nectar containing 20 µg/L and 40 µg/L imidacloprid. These results suggest that the effects of neonicotinoids on honey bee decision-making and other advanced cognitive functions should be explored. Moreover, research should extend beyond the classic model, the European honey bee (A. mellifera), to other important bee species.     

Diet and cell size both affect queen-worker differentiation through DNA methylation in honey bees (Apis mellifera, Apidae)

Background

Young larvae of the honey bee (Apis mellifera) are totipotent; they can become either queens (reproductives) or workers (largely sterile helpers). DNA methylation has been shown to play an important role in this differentiation. In this study, we examine the contributions of diet and cell size to caste differentiation.

Methodology/Principal Findings

We measured the activity and gene expression of one key enzyme involved in methylation, Dnmt3; the rates of methylation in the gene dynactin p62; as well as morphological characteristics of adult bees developed either from larvae fed with worker jelly or royal jelly; and larvae raised in either queen or worker cells. We show that both diet type and cell size contributed to the queen-worker differentiation, and that the two factors affected different methylation sites inside the same gene dynactin p62.

Conclusions/Significance

We confirm previous findings that Dnmt3 plays a critical role in honey bee caste differentiation. Further, we show for the first time that cell size also plays a role in influencing larval development when diet is kept the same.     

Genome-wide analysis of signatures of selection in populations of African honey bees (Apis mellifera) using new web-based tools

Background

With the development of inexpensive, high-throughput sequencing technologies, it has become feasible to examine questions related to population genetics and molecular evolution of non-model species in their ecological contexts on a genome-wide scale. Here, we employed a newly developed suite of integrated, web-based programs to examine population dynamics and signatures of selection across the genome using several well-established tests, including F_ST, pN/pS, and McDonald-Kreitman. We applied these techniques to study populations of honey bees (Apis mellifera) in East Africa. In Kenya, there are several described A. mellifera subspecies, which are thought to be localized to distinct ecological regions.

Results

We performed whole genome sequencing of 11 worker honey bees from apiaries distributed throughout Kenya and identified 3.6 million putative single-nucleotide polymorphisms. The dense coverage allowed us to apply several computational procedures to study population structure and the evolutionary relationships among the populations, and to detect signs of adaptive evolution across the genome. While there is considerable gene flow among the sampled populations, there are clear distinctions between populations from the northern desert region and those from the temperate, savannah region. We identified several genes showing population genetic patterns consistent with positive selection within African bee populations, and between these populations and European A. mellifera or Asian Apis florea.

Conclusions

These results lay the groundwork for future studies of adaptive ecological evolution in honey bees, and demonstrate the use of new, freely available web-based tools and workflows (http://usegalaxy.org/r/kenyanbee) that can be applied to any model system with genomic information.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1712-0) contains supplementary material, which is available to authorized users.     

A meta-analysis of effects of Bt crops on honey bees (Hymenoptera: Apidae)

Background

Honey bees (Apis mellifera L.) are the most important pollinators of many agricultural crops worldwide and are a key test species used in the tiered safety assessment of genetically engineered insect-resistant crops. There is concern that widespread planting of these transgenic crops could harm honey bee populations.

Methodology/Principal Findings

We conducted a meta-analysis of 25 studies that independently assessed potential effects of Bt Cry proteins on honey bee survival (or mortality). Our results show that Bt Cry proteins used in genetically modified crops commercialized for control of lepidopteran and coleopteran pests do not negatively affect the survival of either honey bee larvae or adults in laboratory settings.

Conclusions/Significance

Although the additional stresses that honey bees face in the field could, in principle, modify their susceptibility to Cry proteins or lead to indirect effects, our findings support safety assessments that have not detected any direct negative effects of Bt crops for this vital insect pollinator.     

A Comparative Study of Relational Learning Capacity in Honeybees (Apis mellifera) and Stingless Bees (Melipona rufiventris)

Background

Learning of arbitrary relations is the capacity to acquire knowledge about associations between events or stimuli that do not share any similarities, and use this knowledge to make behavioural choices. This capacity is well documented in humans and vertebrates, and there is some evidence it exists in the honeybee (Apis mellifera). However, little is known about whether the ability for relational learning extends to other invertebrates, although many insects have been shown to possess excellent learning capacities in spite of their small brains.

Methodology/Principal Findings

Using a symbolic matching-to-sample procedure, we show that the honeybee Apis mellifera rapidly learns arbitrary relations between colours and patterns, reaching 68.2% correct choice for pattern-colour relations and 73.3% for colour-pattern relations. However, Apis mellifera does not transfer this knowledge to the symmetrical relations when the stimulus order is reversed. A second bee species, the stingless bee Melipona rufiventris from Brazil, seems unable to learn the same arbitrary relations between colours and patterns, although it exhibits excellent discrimination learning.

Conclusions/Significance

Our results confirm that the capacity for learning arbitrary relations is not limited to vertebrates, but even insects with small brains can perform this learning task. Interestingly, it seems to be a species-specific ability. The disparity in relational learning performance between the two bee species we tested may be linked to their specific foraging and recruitment strategies, which evolved in adaptation to different environments.     

Structural and functional differences in the antennal olfactory system of worker honey bees of Apis mellifera and Apis cerana

Olfactory cues are important sensory modalities on individual discrimination, perception, and efficient orientation to food sources in most insects. In honey bees, which are well known as eusocial insects, olfactory cues are mainly used to maintain a colony. Although much research has been reported on olfactory systems in honey bee olfaction, little is known about the differences between two major honey bee species, the European honey bee Apis mellifera and the Asian honey bee Apis cerana. In order to understand the differences of olfactory characteristics in the two species, we compared the distribution of sensory hairs on the antennae and antennal olfactory responses, using electron microscopy, electrophysiological recording and molecular expression level of odorant receptors. Our present study demonstrated that the antennae of A. cerana have more olfactory sensilla than A. mellifera, responding more strongly to various floral volatile compounds. At the molecular level, olfactory co-receptor (Orco), which makes heterodimers with other conventional olfactory receptors, is more abundantly expressed in the antenna of A. cerana than of A. mellifera. These findings extend our understanding of the olfactory systems and behavioral responses to various ecological and biological signals in two closely related honey bee species.     

Pyrosequencing analysis of the bacterial communities in the guts of honey bees Apis cerana and Apis mellifera in Korea

The bacterial communities in the guts of the adults and larvae of the Asian honey bee Apis cerana and the European honey bee Apis mellifera were surveyed by pyrosequencing the 16S rRNA genes. Most of the gut bacterial 16S rRNA gene sequences were highly similar to the known honey bee-specific ones and affiliated with Pasteurellaceae or lactic acid bacteria (LAB). The numbers of operational taxonomic units (OTUs, defined at 97% similarity) were lower in the larval guts (6 or 9) than in the adult guts (18 or 20), and the frequencies of Pasteurellaceae-related OTUs were higher in the larval guts while those of LAB-related OTUs in the adult guts. The frequencies of Lactococcus, Bartonella, Spiroplasma, Enterobacteriaceae, and Flavobacteriaceae-related OTUs were much higher in A. cerana guts while Bifidobacterium and Lachnospiraceae-related OTUs were more abundant in A. mellfera guts. The bacterial community structures in the midguts and hindguts of the adult honey bees were not different for A. cerana, but significantly different for A. mellifera. The above results substantiated the previous observation that honey bee guts are dominated by several specific bacterial groups, and also showed that the relative abundances of OTUs could be markedly changed depending on the developmental stage, the location within the gut, and the honey bee species. The possibility of using the gut bacterial community as an indicator of honey bee health was discussed.     

Comparative Sucrose Responsiveness in Apis mellifera and A. cerana Foragers

In the European honey bee, Apis mellifera, pollen foragers have a higher sucrose responsiveness than nectar foragers when tested using a proboscis extension response (PER) assay. In addition, Africanized honey bees have a higher sucrose responsiveness than European honey bees. Based on the biology of the Eastern honey bee, A. cerana, we hypothesized that A. cerana should also have a higher responsiveness to sucrose than A. mellifera. To test this hypothesis, we compared the sucrose thresholds of pollen foragers and nectar foragers in both A. cerana and A. mellifera in Fujian Province, China. Pollen foragers were more responsive to sucrose than nectar foragers in both species, consistent with previous studies. However, contrary to our hypothesis, A. mellifera was more responsive than A. cerana. We also demonstrated that this higher sucrose responsiveness in A. mellifera was not due to differences in the colony environment by co-fostering two species of bees in the same mixed-species colonies. Because A. mellifera foragers were more responsive to sucrose, we predicted that their nectar foragers should bring in less concentrated nectar compared to that of A. cerana. However, we found no differences between the two species. We conclude that A. cerana shows a different pattern in sucrose responsiveness from that of Africanized bees. There may be other mechanisms that enable A. cerana to perform well in areas with sparse nectar resources.     

Inhibiting DNA methylation alters olfactory extinction but not acquisition learning in Apis cerana and Apis mellifera

DNA methylation plays a key role in invertebrate acquisition and extinction memory. Honey bees have excellent olfactory learning, but the role of DNA methylation in memory formation has, to date, only been studied in Apis mellifera. We inhibited DNA methylation by inhibiting DNA methyltransferase (DNMT) with zebularine (zeb) and studied the resulting effects upon olfactory acquisition and extinction memory in two honey bee species, Apis cerana and A. mellifera. We used the proboscis extension reflex (PER) assay to measure memory. We provide the first demonstration that DNA methylation is also important in the olfactory extinction learning of A. cerana. DNMT did not reduce acquisition learning in either species. However, zeb bidirectionally and differentially altered extinction learning in both species. In particular, zeb provided 1 h before acquisition learning improved extinction memory retention in A. mellifera, but reduced extinction memory retention in A. cerana. The reasons for these differences are unclear, but provide a basis for future studies to explore species-specific differences in the effects of methylation on memory formation.     

Widespread Occurrence of Chemical Residues in Beehive Matrices from Apiaries Located in Different Landscapes of Western France

Background

The honey bee, Apis mellifera, is frequently used as a sentinel to monitor environmental pollution. In parallel, general weakening and unprecedented colony losses have been reported in Europe and the USA, and many factors are suspected to play a central role in these problems, including infection by pathogens, nutritional stress and pesticide poisoning. Honey bee, honey and pollen samples collected from eighteen apiaries of western France from four different landscape contexts during four different periods in 2008 and in 2009 were analyzed to evaluate the presence of pesticides and veterinary drug residues.

Methodology/Findings

A multi-residue analysis of 80 compounds was performed using a modified QuEChERS method, followed by GC-ToF and LC−MS/MS. The analysis revealed that 95.7%, 72.3% and 58.6% of the honey, honey bee and pollen samples, respectively, were contaminated by at least one compound. The frequency of detection was higher in the honey samples (n = 28) than in the pollen (n = 23) or honey bee (n = 20) samples, but the highest concentrations were found in pollen. Although most compounds were rarely found, some of the contaminants reached high concentrations that might lead to adverse effects on bee health. The three most frequent residues were the widely used fungicide carbendazim and two acaricides, amitraz and coumaphos, that are used by beekeepers to control Varroa destructor. Apiaries in rural-cultivated landscapes were more contaminated than those in other landscape contexts, but the differences were not significant. The contamination of the different matrices was shown to be higher in early spring than in all other periods.

Conclusions/Significance

Honey bees, honeys and pollens are appropriate sentinels for monitoring pesticide and veterinary drug environmental pollution. This study revealed the widespread occurrence of multiple residues in beehive matrices and suggests a potential issue with the effects of these residues alone or in combination on honey bee health.