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.     

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.     

Genome-wide characterization of long intergenic non-coding RNAs (lincRNAs) provides new insight into viral diseases in honey bees Apis cerana and Apis mellifera

Background

Long non-coding RNAs (lncRNAs) are a class of RNAs that do not encode proteins. Recently, lncRNAs have gained special attention for their roles in various biological process and diseases.

Results

In an attempt to identify long intergenic non-coding RNAs (lincRNAs) and their possible involvement in honey bee development and diseases, we analyzed RNA-seq datasets generated from Asian honey bee (Apis cerana) and western honey bee (Apis mellifera). We identified 2470 lincRNAs with an average length of 1011 bp from A. cerana and 1514 lincRNAs with an average length of 790 bp in A. mellifera. Comparative analysis revealed that 5 % of the total lincRNAs derived from both species are unique in each species. Our comparative digital gene expression analysis revealed a high degree of tissue-specific expression among the seven major tissues of honey bee, different from mRNA expression patterns. A total of 863 (57 %) and 464 (18 %) lincRNAs showed tissue-dependent expression in A. mellifera and A. cerana, respectively, most preferentially in ovary and fat body tissues. Importantly, we identified 11 lincRNAs that are specifically regulated upon viral infection in honey bees, and 10 of them appear to play roles during infection with various viruses.

Conclusions

This study provides the first comprehensive set of lincRNAs for honey bees and opens the door to discover lincRNAs associated with biological and hormone signaling pathways as well as various diseases of honey bee.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1868-7) 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.     

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.     

Differential gene expression of the honey bees Apis mellifera and A. cerana induced by Varroa destructor infection

Varroa destructor mite is currently the most serious threat to the world bee industry. Differences in mite tolerance are reported between two honey bee species Apis mellifera and Apis cerana. Differential gene expression of two honey bee species induced by V. destructor infection was investigated by constructing two suppression subtractive hybridization (SSH) libraries, as first steps toward elucidating molecular mechanisms of Varroa tolerance. From the SSH libraries, we obtained 289 high quality sequences which clustered into 132 unique sequences grouped in 26 contigs and 106 singlets where 49 consisted in A. cerana subtracted library and 83 in A. mellifera. Using BLAST, we found that 85% sequences had counterpart known genes whereas 15% were undescribed. A Gene Ontology analysis classified 51 unique sequences into different functional categories. Eight of these differentially expressed genes, representative of different regulation patterns, were confirmed by qRT-PCR. Upon the mite induction, the differentially expressed genes from both bee species were different, except hex 110 gene, which was up-regulated in A. cerana but down-regulated in A. mellifera, and Npy-r gene, which was down-regulated in both species. In general, most of the differential expression genes were involved in metabolic processes and nerve signaling. The results provide information on the molecular response of these two bee species to Varroa infection.     

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.     

Comparative susceptibility and immune responses of Asian and European honey bees to the American foulbrood pathogen,Paenibacillus larvae

American foulbrood (AFB) disease is caused by Paenibacillus larvae. Currently, this pathogen is widespread in the European honey bee— Apis mellifera. However, little is known about infectivity and pathogenicity of P. lan'ae in the Asiatic cavity-nesting honey bees, Apis cerana. Moreover, comparative knowledge of P. larvae infectivity and pathogenicity between both honey bee species is scarce. In this study, we examined susceptibility, larval mortality, survival rate and expression of genes encoding antimicrobial peptides (AMPs) including defensin, apidaecin, abaecin, and hymenoptaecin in A. mellifera and A. cerana when infected with P. larvae. Our results showed similar effects of P. larvae on the survival rate and patterns of AMP gene expression in both honey bee species when bee larvae are infected with spores at the median lethal concentration (LC5 0 ) for A. mellifera. All AMPs of infected bee larvae showed significant upregulation compared with noninfected bee larvae in both honey bee species. However, larvae of A. cerana were more susceptible than A. mellifera when the same larval ages and spore concentration of P. larvae were used. It also appears that A. cerana showed higher levels of AMP expression than A. mellifera. This research provides the first evidence of survival rate, LC50 and immune response profiles of Asian honey bees, A. cerana, when infected by P. larvae in comparison with the European honey bee, A. mellifera.     

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.     

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.     

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.     

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.     

Flower visitation patterns of the coexisting honey bees Apis cerana japonica and Apis mellifera (Hymenoptera: Apidae)

We compared flower visitation patterns of two coexisting honey bees, Apis mellifera Linnaeus and Apis cerana japonica Radoszkowski, on 20 plant species, including three exotics, under natural conditions in Nara, Japan, from April to August 2012. We also measured flower color based on bee color vision (15 flower species), nectar volume (nine species) and nectar concentration (eight species). Flowers colored white, pink, red, purple and cream were classified as bee‐blue‐green, and yellow was classified as bee‐green. Apis cerana visited 14 plant species and A. mellifera visited 11. Although the two Apis species are similar in morphology, they visited different plants: in particular, A. cerana visited native plant species more often than did A. mellifera. Both A. mellifera and A. cerana visited not only nectariferous flowers but also those with no nectar. We also found different visitation patterns between A. cerana and A. mellifera: Apis cerana more often visited flowers with smaller color angle (bee‐blue‐green), lower chroma and higher brightness, and flowers secreting nectars of higher concentration and smaller volume than did A. mellifera.     

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.     

Genetic diversity and phylogenetic relationship among the western and the Asian honey bees based on two mitochondrial gene segments (COI and ND5)

The Asian honey bee species i.e., Apis cerana (the eastern honey bee), A. dorsata (the giant honey bee), and the western or European honey bee (A. mellifera) collected from Pakistan were studied using partial sequences from two mitochondrial genes (i) the Cytochrome c oxidase I (COI) and (ii) the mitochondrially encoded NADH dehydrogenase 5 (ND5) and then compared with other honey bees sequences (already submitted from different countries around the globe) obtained after the national center for biotechnology information (NCBI). DNA sequences were analyzed employing molecular evolutionary genetics analysis and Kimura 2-parameter model, neighbor-joining method was applied to investigate phylogenetic relationships, and DNA sequence polymorphism was applied to measure the genetic diversity within the genus Apis. The phylogenetic analyses yielded consistent results. Based on COI gene fragment in two Asian and European honey bee species from Pakistan and from other countries showed considerable genetic diversity levels and deviation among the species. While in contrast the phylogenetic analyses based on ND5 gene fragment in Asian and European honey bee species from Pakistan and other countries showed comparatively higher genetic diversity indices and variations than the COI gene. So, in the genus Apis, the mitochondrial ND5 region has shown the possibility to answer the interactions among species. A further detailed work (by linking the analysis of other genomic and mitochondrial genes) is required for good quality solution to establish the concise genetic diversity and interaction among the Apis species. The objective of this study was to explore the extent of genetic differences and phylogenetic links among the three kinds of honey bee species from Pakistan and comparing them with other bee species around the globe.     

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.     

Recipe for a Busy Bee: MicroRNAs in Honey Bee Caste Determination

Social caste determination in the honey bee is assumed to be determined by the dietary status of the young larvae and translated into physiological and epigenetic changes through nutrient-sensing pathways. We have employed Illumina/Solexa sequencing to examine the small RNA content in the bee larval food, and show that worker jelly is enriched in miRNA complexity and abundance relative to royal jelly. The miRNA levels in worker jelly were 7–215 fold higher than in royal jelly, and both jellies showed dynamic changes in miRNA content during the 4^th to 6^th day of larval development. Adding specific miRNAs to royal jelly elicited significant changes in queen larval mRNA expression and morphological characters of the emerging adult queen bee. We propose that miRNAs in the nurse bee secretions constitute an additional element in the regulatory control of caste determination in the honey bee.     

Effect of royal jelly on longevity and memory-related traits of Apis mellifera workers

Royal jelly (RJ) is a key factor for honey bee caste determination. The queen bee is fed with RJ by worker bees throughout her life, while the worker bees eat bee bread themselves. This study was designed to explore the effect of nutrient-rich RJ on longevity and learning and memory abilities of workers of the western honey bee Apis mellifera. The newly emerged worker bees were randomly divided into three groups and were fed 50% sucrose solution containing 0%, 10%, and 20% RJ. We found that worker bees fed with 10% and 20% RJ showed significantly improved longevity and higher proboscis extension response success rate compared to bees fed with 50% sucrose containing 0% RJ. Additionally, bees fed with 20% RJ showed significantly higher level of expression of memory related genes (GluRA and Nmdar1) compared to the control group. Furthermore, expression of the Nmdar1 gene of worker bees fed with 10% RJ was also significantly higher than in the control group. These results indicate that RJ has potential effects on the longevity and learning and memory abilities of A. mellifera.