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.     

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.     

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.     

Varroa destructor changes its cuticular hydrocarbons to mimic new hosts

Varroa destructor (Vd) is a honeybee ectoparasite. Its original host is the Asian honeybee, Apis cerana, but it has also become a severe, global threat to the European honeybee, Apis mellifera. Previous studies have shown that Varroa can mimic a host''s cuticular hydrocarbons (HC), enabling the parasite to escape the hygienic behaviour of the host honeybees. By transferring mites between the two honeybee species, we further demonstrate that Vd is able to mimic the cuticular HC of a novel host species when artificially transferred to this new host. Mites originally from A. cerana are more efficient than mites from A. mellifera in mimicking HC of both A. cerana and A. mellifera. This remarkable adaptability may explain their relatively recent host-shift from A. cerana to A. mellifera.     

Thermoregulation in mixed-species colonies of honeybees (Apis cerana and Apis mellifera)

Apis cerana and Apis mellifera normally display different strategies in cooling hive temperature, raising the question whether they would coordinate their efforts in to achieve stable thermoregulation in mixed colonies. The results show that the normal temperatures in the brood area in mixed colonies are more similar to those of pure A. cerana colonies than pure A. mellifera colonies. Under heat stress, A. cerana workers are more sensitive, and initiate fanning earlier than A. mellifera workers. In mixed colonies, the former become the main force for thermoregulation. When worker bees of both species were fanning together at the entrance, their own species-specific postures were adopted, but due to a significantly smaller number of A. mellifera workers engaged in fanning, the cooling efficiency of mixed colonies were closest to that of pure A. cerana colonies.     

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.     

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.     

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.     

Phylogenetic analysis of Nosema ceranae isolated from European and Asian honeybees in Northern Thailand

Nosema ceranae was found to infect four different host species including the European honeybee (A. mellifera) and the Asian honeybees (Apis florea, A. cerana and Apis dorsata) collected from apiaries and forests in Northern Thailand. Significant sequence variation in the polar tube protein (PTP1) gene of N. ceranae was observed with N. ceranae isolates from A. mellifera and A. cerana, they clustered into the same phylogenetic lineage. N. ceranae isolates from A. dorsata and A. florea were grouped into two other distinct clades. This study provides the first elucidation of a genetic relationship among N. ceranae strains isolated from different host species and demonstrates that the N. ceranae PTP gene was shown to be a suitable and reliable marker in revealing genetic relationships within species.     

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.     

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.     

Pathological effects of the microsporidium Nosema ceranae on honey bee queen physiology (Apis mellifera)

Nosema ceranae, a microsporidian parasite originally described in the Asian honey bee Apis cerana, has recently been found to be cross-infective and to also parasitize the European honey bee Apis mellifera. Since this discovery, many studies have attempted to characterize the impact of this parasite in A. mellifera honey bees. Nosema species can infect all colony members, workers, drones and queens, but the pathological effects of this microsporidium has been mainly investigated in workers, despite the prime importance of the queen, who monopolizes the reproduction and regulates the cohesion of the society via pheromones. We therefore analyzed the impact of N. ceranae on queen physiology. We found that infection by N. ceranae did not affect the fat body content (an indicator of energy stores) but did alter the vitellogenin titer (an indicator of fertility and longevity), the total antioxidant capacity and the queen mandibular pheromones, which surprisingly were all significantly increased in Nosema-infected queens. Thus, such physiological changes may impact queen health, leading to changes in pheromone production, that could explain Nosema-induced supersedure (queen replacement).     

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.     

Infections of Nosema ceranae in four different honeybee species

The microsporidium Nosema ceranae is detected in honeybees in Thailand for the first time. This endoparasite has recently been reported to infect most Apis mellifera honeybee colonies in Europe, the US, and parts of Asia, and is suspected to have displaced the endemic endoparasite species, Nosema apis, from the western A. mellifera. We collected and identified species of microsporidia from the European honeybee (A. mellifera), the cavity nesting Asian honeybee (Apis cerana), the dwarf Asian honeybee (Apis florea) and the giant Asian honeybee (Apis dorsata) from colonies in Northern Thailand. We used multiplex PCR technique with two pairs of primers to differentiate N. ceranae from N. apis. From 80 A. mellifera samples, 62 (77.5%) were positively identified for the presence of the N. ceranae. Amongst 46 feral colonies of Asian honeybees (A. cerana, A. florea and A. dorsata) examined for Nosema infections, only N. ceranae could be detected. No N. apis was found in our samples. N. ceranae is found to be the only microsporidium infesting honeybees in Thailand. Moreover, we found the frequencies of N. ceranae infection in native bees to be less than that of A. mellifera.     

A comprehensive genetic map of sugarcane that provides enhanced map coverage and integrates high-throughput Diversity Array Technology (DArT) markers

Background

Sugarcane genetic mapping has lagged behind other crops due to its complex autopolyploid genome structure. Modern sugarcane cultivars have from 110-120 chromosomes and are in general interspecific hybrids between two species with different basic chromosome numbers: Saccharum officinarum (2n = 80) with a basic chromosome number of 10 and S. spontaneum (2n = 40-128) with a basic chromosome number of 8. The first maps that were constructed utilised the single dose (SD) markers generated using RFLP, more recent maps generated using AFLP and SSRs provided at most 60% genome coverage. Diversity Array Technology (DArT) markers are high throughput allowing greater numbers of markers to be generated.

Results

Progeny from a cross between a sugarcane variety Q165 and a S. officinarum accession IJ76-514 were used to generate 2467 SD markers. A genetic map of Q165 was generated containing 2267 markers, These markers formed 160 linkage groups (LGs) of which 147 could be placed using allelic information into the eight basic homology groups (HGs) of sugarcane. The HGs contained from 13 to 23 LGs and from 204 to 475 markers with a total map length of 9774.4 cM and an average density of one marker every 4.3 cM. Each homology group contained on average 280 markers of which 43% were DArT markers 31% AFLP, 16% SSRs and 6% SNP markers. The multi-allelic SSR and SNP markers were used to place the LGs into HGs.

Conclusions

The DArT array has allowed us to generate and map a larger number of markers than ever before and consequently to map a larger portion of the sugarcane genome. This larger number of markers has enabled 92% of the LGs to be placed into the 8 HGs that represent the basic chromosome number of the ancestral species, S. spontaneum. There were two HGs (HG2 and 8) that contained larger numbers of LGs verifying the alignment of two sets of S. officinarum chromosomes with one set of S. spontaneum chromosomes and explaining the difference in basic chromosome number between the two ancestral species. There was also evidence of more complex structural differences between the two ancestral species.

Electronic supplementary material

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

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.     

Nosema ceranae, a new parasite in Thai honeybees

Adult workers of Apis cerana, Apis florea and Apis mellifera from colonies heavily infected with Nosema ceranae were selected for molecular analyses of the parasite. PCR-specific 16S rRNA primers were designed, cloned, sequenced and compared to GenBank entries. The sequenced products corresponded to N. ceranae. We then infected A. cerana with N. ceranae spores isolated from A. florea workers. Newly emerged bees from healthy colonies were fed 10,000, 20,000 and 40,000 spores/bee. There were significant dosage dependent differences in bee infection and survival rates. The ratio of infected cells to non-infected cells increased at 6, 10 and 14 d post infection. In addition, hypopharyngeal glands of bees from the control group had significantly higher protein concentrations than infected groups. Bees infected with 40,000 spores/bee had the lowest protein concentrations. Thus, N. ceranae isolated from A. florea is capable of infecting another bee species, impairing hypopharyngeal gland protein production and reducing bee survival in A. cerana.     

Antioxidant capacity of major royal jelly proteins of honeybee (Apis mellifera) royal jelly

Major royal jelly proteins (MRJPs) of honeybee royal jelly (RJ) exhibit antimicrobial and antioxidant activities. Although MRJPs of Apis mellifera RJ (AmMRJPs) responsible for antibacterial activity have been identified, AmMRJPs with antioxidant effects remain to be elucidated. Here we identified and compared the antioxidant activities of purified recombinant AmMRJPs 1–7, which are expressed in baculovirus-infected insect cells. Antioxidant assays of recombinant AmMRJPs 1–7 against H₂O₂ revealed that AmMRJPs reduce caspase-3 activity and oxidative stress-induced cell apoptosis and lead to increased cell viability. Consistent with these results, AmMRJPs 1–7 exhibit 1,1-diphenyl-2-picrylhydrazyl radical-scavenging activity and protect against oxidative DNA damage. These results indicate that AmMRJPs play a role as antioxidants in A. mellifera RJ.