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.     

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.     

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.     

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.     

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.     

Asymmetrical coexistence of Nosema ceranae and Nosema apis in honey bees

Globalization has provided opportunities for parasites/pathogens to cross geographic boundaries and expand to new hosts. Recent studies showed that Nosema ceranae, originally considered a microsporidian parasite of Eastern honey bees, Apis cerana, is a disease agent of nosemosis in European honey bees, Apis mellifera, along with the resident species, Nosema apis. Further studies indicated that disease caused by N. ceranae in European honey bees is far more prevalent than that caused by N. apis. In order to gain more insight into the epidemiology of Nosema parasitism in honey bees, we conducted studies to investigate infection of Nosema in its original host, Eastern honey bees, using conventional PCR and duplex real time quantitative PCR methods. Our results showed that A. cerana was infected not only with N. ceranae as previously reported [Fries, I., Feng, F., Silva, A.D., Slemenda, S.B., Pieniazek, N.J., 1996. Nosema ceranae n. sp. (Microspora, Nosematidae), morphological and molecular characterization of a microsporidian parasite of the Asian honey bee Apis cerana (Hymenoptera, Apidae). Eur. J. Protistol. 32, 356-365], but also with N. apis. Both microsporidia produced single and mixed infections. Overall and at each location alone, the prevalence of N. ceranae was higher than that of N. apis. In all cases of mixed infections, the number of N. ceranae gene copies (corresponding to the parasite load) significantly out numbered those of N. apis. Phylogenetic analysis based on a variable region of small subunit ribosomal RNA (SSUrRNA) showed four distinct clades of N. apis and five clades of N. ceranae and that geographical distance does not appear to influence the genetic diversity of Nosema populations. The results from this study demonstrated that duplex real-time qPCR assay developed in this study is a valuable tool for quantitative measurement of Nosema and can be used to monitor the progression of microsprodian infections of honey bees in a timely and cost efficient manner.     

The role of genetic diversity in nest cooling in a wild honey bee, Apis florea

Simulation studies of the task threshold model for task allocation in social insect colonies suggest that nest temperature homeostasis is enhanced if workers have slightly different thresholds for engaging in tasks related to nest thermoregulation. Genetic variance in task thresholds is one way a distribution of task thresholds can be generated. Apis mellifera colonies with large genetic diversity are able to maintain more stable brood nest temperatures than colonies that are genetically uniform. If this phenomenon is generalizable to other species, we would predict that patrilines should vary in the threshold in which they engage in thermoregulatory tasks. We exposed A. florea colonies to different temperatures experimentally, and retrieved fanning workers at these different temperatures. In many cases we found statistically significant differences in the proportion of fanning workers of different patrilines at different experimental temperatures. This suggests that genetically different workers have different thresholds for performing the thermoregulatory task of fanning. We suggest, therefore, that genetically based variance in task threshold is a widespread phenomenon in the genus Apis.     

Deformed wing virus in western honey bees (Apis mellifera) from Atlantic Canada and the first description of an overtly-infected emerging queen

Deformed wing virus (DWV) in western honey bees (Apis mellifera) often remains asymptomatic in workers and drones, and symptoms have never been described from queens. However, intense infections linked to parasitism by the mite Varroa destructor can cause worker wing deformity and death within 67 h of emergence. Ten workers (eight with deformed wings and two with normal wings) and three drones (two with deformed wings and one with normal wings) from two colonies infected with V. destructor from Nova Scotia, Canada, and two newly-emerged queens (one with deformed wings and one with normal wings) from two colonies infected with V. destructor from Prince Edward Island, Canada, were genetically analyzed for DWV. We detected DWV in all workers and drones, regardless of wing morphology, but only in the deformed-winged queen. This is the first report of DWV from Atlantic Canada and the first detection of a symptomatic queen with DWV from anywhere.     

Nestmate Recognition Differences between Honeybee Colonies of <Emphasis Type="BoldItalic">Apis cerana</Emphasis> and <Emphasis Type="BoldItalic">Apis mellifera</Emphasis>

Nestmate recognition in Apis cerana and Apis mellifera was studied by introducing sealed queen cells heterospecifically between queenless colonies. No A. cerana queens were accepted by queenless A. mellifera; but A. mellifera queens were accepted in queenless A. cerana colonies. A. mellifera queens oviposited in queenless A. cerana colonies, but A. cerana workers removed most eggs. In time, egg removals declined, and some A. mellifera larvae that hatched from these eggs reached adulthood, and eventually about half of the workers were newly emerged A. mellifera. Eventually, the colonies consisted only of A. mellifera after A. cerana workers died by attrition. A. mellifera workers are more sensitive to nestmate recognition and killed the A. cerana virgin queens. In mixed-species colonies, after newly emerged A. mellifera workers matured, they removed eggs laid by the A. cerana queens until there were no workers to replace the old ones.     

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.     

Differences in venom and cuticular peptides in individuals of Apis mellifera (Hymenoptera: Apidae) determined by MALDI-TOF MS

The fraction between 950 and 4000 Da of the venom of Apis mellifera has been analyzed with MALDI-TOF mass spectrometry and statistical facilities of the ClinProTools™ software. Consistent differences in the composition of this venom fraction were observed between queens and workers while younger and older workers (nurses and guards as well as foragers) differ for the relative percentages of two well known cytolytic peptides, namely Melittin and Apamin. Total in situ body methanol extracts and methanol micro-extractions on the cuticle of various parts of the body of drones and females confirmed that venom peptides are smeared on the body surface of females in a not yet clarified way. The observation that venom peptides have been found also on comb wax rises the hypothesis that the use of venom as antimicrobial agent makes part of the social immunity system of A. mellifera.     

Expression and characterization of the chitinases from Serratia marcescens GEI strain for the control of Varroa destructor, a honey bee parasite

Serratia marcescens GEI strain was isolated from the gut of the workers of Chinese honey bee Apis cerana and evaluated in the laboratory for the control of Varroa destructor, a parasite of western honey bee A. mellifera. The supernatant and the collected proteins by ammonium sulfate from the bacterial cultures showed a strong miticidal effect on the female mites, with 100% mite mortality in 5 days. Heat (100 °C for 10 min) and proteinase K treatment of the collected proteins destroyed the miticidal activity. The improved miticial activity of this bacterial strain on chitin medium indicated the involvement of chitinases. The expressed chitinases ChiA, ChiB and ChiC1 from S. marcescens GEI by recombinant Escherichia coli showed pathogenicity against the mites in the laboratory. These chitinases were active in a broad pH range (5-9) and the optimum temperatures were between 60 and 75 °C. Synergistic effects of ChiA and ChiB on the miticidal activity against V. destructor were observed. The workers of both honey bee species were not sensitive to the spraying and feeding chitinases. These results provided alternative control strategies for Varroa mites, by formulating chitinase agents and by constructing transgenetic honey bees.     

Effects at Nearctic north-temperate latitudes of indoor versus outdoor overwintering on the microsporidium Nosema ceranae and western honey bees (Apis mellifera)

In northern temperate climates, western honey bee (Apis mellifera) colonies can be wintered outdoors exposed to ambient conditions, or indoors in a controlled setting. Because very little is known about how this affects the recently-detected microsporidium Nosema ceranae, we investigated effects of indoor versus outdoor overwintering on spring N. ceranae intensity (spores per bee), and on winter and spring colony mortality. For colonies medicated with Fumagilin-B® to control N. ceranae, overwintering treatment did not affect N. ceranae intensity, despite outdoor-wintered colonies having significantly greater mortality. These findings suggest that N. ceranae may not always pose the most significant threat to western honey bees, and that indoor-wintering may ensure that a greater number of colonies are available for honey production and pollination services during the summer.     

Crossfostered drones ofApis cerana (Fabricius, 1793) andApis koschevnikovi (Buttel-Reepen, 1906) fly at their species specific mating times

Summary Reciprocal transfer of sealed drone brood between colonies ofApis cerana andApis koschevnikovi was successful and resulted in four colonies (two of each species) with a mixed drone population. Flights ofApis cerana drones occurred between 14.00 and 16.15 regardless whether they were in a conspecific or alien colony.Apis koschevnikovi drones also flew at their species specific time from 16.45 to 18.30. A variance estimation revealed that 99.4% of the total variance depended on the species of the drone. In contrast to theApis drone's general biological dependence upon the colony, crossfostered drones ofApis cerana andApis koschevnikovi showed an unexpected autonomy in chosing their mating flight time.     

The pheromones of laying workers in two honeybee sister species: <Emphasis Type="Italic">Apis cerana</Emphasis> and <Emphasis Type="Italic">Apis mellifera</Emphasis>

When a honeybee colony loses its queen, workers activate their ovaries and begin to lay eggs. This is accompanied by a shift in their pheromonal bouquet, which becomes more queen like. Workers of the Asian hive bee Apis cerana show unusually high levels of ovary activation and this can be interpreted as evidence for a recent evolutionary arms race between queens and workers over worker reproduction in this species. To further explore this, we compared the rate of pheromonal bouquet change between two honeybee sister species of Apis cerana and Apis mellifera under queenright and queenless conditions. We show that in both species, the pheromonal components HOB, 9-ODA, HVA, 9-HDA, 10-HDAA and 10-HDA have significantly higher amounts in laying workers than in non-laying workers. In the queenright colonies of A. mellifera and A. cerana, the ratios (9-ODA)/(9-ODA + 9-HDA + 10-HDAA + 10-HDA) are not significantly different between the two species, but in queenless A. cerana colonies the ratio is significant higher than in A. mellifera, suggesting that in A. cerana, the workers’ pheromonal bouquet is dominated by the queen compound, 9-ODA. The amount of 9-ODA in laying A. cerana workers increased by over 585% compared with the non-laying workers, that is 6.75 times higher than in A. mellifera where laying workers only had 86% more 9-ODA compared with non-laying workers.     

Presence of Nosema ceranae in honeybees (Apis mellifera) in Uruguay

The microsporidium Nosema ceranae is an emergent pathogen of European honeybees Apis mellifera. Using a PCR-RFLP diagnosis, 29 samples of infected honeybees obtained in 2007-2008 (N = 26), 2004 (N = 2) and before 1990 (N = 1) were analyzed for the presence of Nosema apis and N. ceranae. Only N. ceranae was found in all samples, indicating that this species dispersed to Uruguay (and likely the region) at some time before 1990. The presence of N. ceranae in Uruguay is not associated with an increase of Nosemosis, and its role in colony loss seems to be irrelevant.     

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.     

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

Nestmate recognition by guards of the Asian hive bee <Emphasis Type="Italic">Apis cerana</Emphasis>

When a honey bee colony becomes queenless and broodless its only reproductive option is for some of its workers to produce sons before the colony perishes. However, for this to be possible the policing of worker-laid eggs must be curtailed and this provides the opportunity for queenless colonies to be reproductively parasitized by workers from other nests. Such reproductive parasitism is known to occur in Apis florea and A. cerana. Microsatellite analyses of worker samples have demonstrated that the proportion of non-natal workers present in an A. cerana colony declines after a colony is made queenless. This observation suggests that queenless A. cerana colonies may be more vigilant in repelling potentially parasitic non-natal workers than queenright colonies. We compared rates of nestmate and non-nestmate acceptance in both queenright and queenless A. cerana colonies using standard assays and showed that there is no statistical difference between the proportion of non-nestmate workers that are rejected in queenless and queenright colonies. We also show that, contrary to earlier reports, A. cerana guards are able to discriminate nestmate workers from non-nestmates, and that they reject significantly more non-nestmate workers than nestmate workers. Received 25 February 2008; revised 21 May 2008; accepted 25 June 2008.     

Differential hygienic behavior of Apis cerana F. and Apis mellifera L. to Sacbrood virus infection

Beekeeping with Apis cerana of Korean apiculture is facing with serious colony collapse caused by invasive Sacbrood virus (SBV) disease. This fatal brood disease was the main reason of more than 90% colony lost in Korea leading almost the extinct crisis. Sacbrood virus can infect either larvae or adult honeybees, with a higher sensibility of larvae to the infection. Since SBV has spread to all over the country, efforts have been made to treat and prevent this devastating disease although no effective results have so far been obtained. Several studies have demonstrated that Apis mellifera bee colonies that express an efficient hygienic behavior exhibit a higher resistance to the brood disease. In this study we demonstrated that the differences of hygienic behavior between A. cerana and A. mellifera. A. cerana more efficiently removed the pin-killed brood than A. mellifera. On the other hand, A. mellifera more efficiently removed SBV-infected larvae and SBV-dead brood than A. cerana. However, it remains unclear whether the advantage of hygienic bee could have efficacy against Sacbrood disease on A. cerana colonies.