The Invasion of the Dwarf Honeybee, Apis florea, along the River Nile in Sudan.

The spread of the dwarf honeybee, Apis florea, in Sudan along the river Nile in a linear fashion provides a good model for studying the population dynamics and genetic effects of an invasion by a honeybee species. We use microsatellite DNA analyses to assess the population structure of both invasive A. florea and native Apis mellifera along the river Nile. The invasive A. florea had significantly higher population densities than the wild, native A. mellifera. Nevertheless, we found no indication of competitive displacement, suggesting that although A. florea had a high invasive potential, it coexisted with the native A. mellifera along the river Nile. The genetic data indicated that the invasion of A. florea was established by a single colony.

A gene for resistance to the Varroa mite (Acari) in honey bee (Apis mellifera) pupae.

Social insect colonies possess a range of defences which protect them against highly virulent parasites and colony collapse. The host-parasite interaction between honey bees (Apis mellifera) and the mite Varroa destructor is unusual, as honey bee colonies are relatively poorly defended against this parasite. The interaction has existed since the mid-20th Century, when Varroa switched host to parasitize A. mellifera. The combination of a virulent parasite and relatively naïve host means that, without acaricides, honey bee colonies typically die within 3 years of Varroa infestation. A consequence of acaricide use has been a reduced selective pressure for the evolution of Varroa resistance in honey bee colonies. However, in the past 20 years, several natural-selection-based breeding programmes have resulted in the evolution of Varroa-resistant populations. In these populations, the inhibition of Varroa's reproduction is a common trait. Using a high-density genome-wide association analysis in a Varroa-resistant honey bee population, we identify an ecdysone-induced gene significantly linked to resistance. Ecdysone both initiates metamorphosis in insects and reproduction in Varroa. Previously, using a less dense genetic map and a quantitative trait loci analysis, we have identified Ecdysone-related genes at resistance loci in an independently evolved resistant population. Varroa cannot biosynthesize ecdysone but can acquire it from its diet. Using qPCR, we are able to link the expression of ecdysone-linked resistance genes to Varroa's meals and reproduction. If Varroa co-opts pupal compounds to initiate and time its own reproduction, mutations in the host's ecdysone pathway may represent a key selection tool for honey bee resistance and breeding.

A Single SNP Turns a Social Honey Bee (Apis mellifera) Worker into a Selfish Parasite.

The evolution of altruism in complex insect societies is arguably one of the major transitions in evolution and inclusive fitness theory plausibly explains why this is an evolutionary stable strategy. Yet, workers of the South African Cape honey bee (Apis mellifera capensis) can reverse to selfish behavior by becoming social parasites and parthenogenetically producing female offspring (thelytoky). Using a joint mapping and population genomics approach, in combination with a time-course transcript abundance dynamics analysis, we show that a single nucleotide polymorphism at the mapped thelytoky locus (Th) is associated with the iconic thelytokous phenotype. Th forms a linkage group with the ecdysis-triggering hormone receptor (Ethr) within a nonrecombining region under strong selection in the genome. A balanced detrimental allele system plausibly explains why the trait is specific to A. m. capensis and cannot easily establish itself into genomes of other honey bee subspecies.

Association of Varroa destructor females in multiply infested cells of the honeybee Apis mellifera.

The genetic diversity of Varroa destructor (Anderson & Trueman) is limited outside its natural range due to population bottlenecks and its propensity to inbreed. In light of the arms race between V. destructor and its honeybee (Apis mellifera L.) host, any mechanism enhancing population admixture of the mite may be favored. One way that admixture can occur is when two genetically dissimilar mites coinvade a brood cell, with the progeny of the foundresses admixing. We determined the relatedness of 393 pairs of V. destructor foundresses, each pair collected from a single bee brood cell (n = five colonies). We used six microsatellites to identify the genotypes of mites coinvading a cell and calculated the frequency of pairs with different or the same genotypes. We found no deviation from random coinvasion, but the frequency of cells infested by mites with different genotypes was high. This rate of recombination, coupled with a high transmission rate of mites, homogenized the allelic pool of mites within the apiary.

The role of epistatic interactions underpinning resistance to parasitic Varroa mites in haploid honey bee (Apis mellifera) drones.

The Red Queen hypothesis predicts that host-parasite coevolutionary dynamics can select for host resistance through increased genetic diversity, recombination and evolutionary rates. However, in haplodiploid organisms such as the honeybee (Apis mellifera), models suggest the selective pressure is weaker than in diploids. Haplodiploid sex determination, found in A. mellifera, can allow deleterious recessive alleles to persist in the population through the diploid sex with negative effects predominantly expressed in the haploid sex. To overcome these negative effects in haploid genomes, epistatic interactions have been hypothesized to play an important role. Here, we use the interaction between A. mellifera and the parasitic mite Varroa destructor to test epistasis in the expression of resistance, through the inhibition of parasite reproduction, in haploid drones. We find novel loci on three chromosomes which explain over 45% of the resistance phenotype. Two of these loci interact only additively, suggesting their expression is independent of each other, but both loci interact epistatically with the third locus. With drone offspring inheriting only one copy of the queen's chromosomes, the drones will only possess one of two queen alleles throughout the years-long lifetime of the honeybee colony. Varroa, in comparison, completes its highly inbred reproductive cycle in a matter of weeks, allowing it to rapidly evolve resistance. Faced with the rapidly evolving Varroa, a diversity of pathways and epistatic interactions for the inhibition of Varroa reproduction could therefore provide a selective advantage to the high levels of recombination seen in A. mellifera. This allows for the remixing of phenotypes despite a fixed queen genotype.

How Honeybees Defy Gravity with Royal Jelly to Raise Queens.

The female sex in honeybees (Apis spp.) comprises a reproductive queen and a sterile worker caste. Nurse bees feed all larvae progressively with a caste-specific food jelly until the prepupal stage. Only those larvae that are exclusively fed a large amount of royal jelly (RJ) develop into queens [1]. RJ is a composite secretion of two specialized head glands: the mandibular glands, which produce mainly fatty acids [2], and the hypopharyngeal glands, which contribute proteins, primarily belonging to the major royal jelly protein (MRJP) family [3]. Past research on RJ has focused on its nutritional function and overlooked its central role with regard to the orientation of the larva in the royal brood cell. Whereas workers are reared in the regular horizontal cells of the comb, the queen cells are specifically built outside of the normal comb area to accommodate for the larger queen [4, 5]. These cells hang freely along the bottom of the comb and are vertically oriented, opening downward [6]. Queen larvae are attached by their RJ diet to the cell ceiling. Thus, the physical properties of RJ are central to successful retention of larvae in the cell. Here, we show that the main protein of RJ (MRJP1) polymerizes in complex with another protein, apisimin, into long fibrous structures that build the basis for the high viscosity of RJ to hold queen larvae on the RJ surface.

Infection dynamics of Nosema ceranae in honey bee midgut and host cell apoptosis.

Nosema ceranae is an intracellular microsporidian parasite that infects epithelial cells of the honey bee (Apis mellifera) midgut. Previous studies have shown that Nosema may alter cell renewal and apoptosis in honey bees. We found that the amount of apoptotic cells progressively declines from the anterior towards posterior regions of the midgut in Nosema-infected sensitive bees. There was no such pattern in the infected Nosema tolerant honey bees and controls. These data provide additional evidence that N. ceranae appears to alter apoptosis in its host cells for its own advantage.

The transcriptomic changes associated with the development of social parasitism in the honeybee Apis mellifera capensis.

Social insects are characterized by the division of labor. Queens usually dominate reproduction, whereas workers fulfill non-reproductive age-dependent tasks to maintain the colony. Although workers are typically sterile, they can activate their ovaries to produce their own offspring. In the extreme, worker reproduction can turn into social parasitism as in Apis mellifera capensis. These intraspecific parasites occupy a host colony, kill the resident queen, and take over the reproductive monopoly. Because they exhibit a queenlike behavior and are also treated like queens by the fellow workers, they are so-called pseudoqueens. Here, we compare the development of parasitic pseudoqueens and social workers at different time points using fat body transcriptome data. Two complementary analysis methods-a principal component analysis and a time course analysis-led to the identification of a core set of genes involved in the transition from a social worker into a highly fecund parasitic pseudoqueen. Comparing our results on pseudoqueens with gene expression data of honeybee queens revealed many similarities. In addition, there was a set of specific transcriptomic changes in the parasitic pseudoqueens that differed from both, queens and social workers, which may be typical for the development of the social parasitism in A. m. capensis.

Neonicotinoid pesticides can reduce honeybee colony genetic diversity.

Neonicotinoid insecticides can cause a variety of adverse sub-lethal effects in bees. In social species such as the honeybee, Apis mellifera, queens are essential for reproduction and colony functioning. Therefore, any negative effect of these agricultural chemicals on the mating success of queens may have serious consequences for the fitness of the entire colony. Queens were exposed to the common neonicotinoid pesticides thiamethoxam and clothianidin during their developmental stage. After mating, their spermathecae were dissected to count the number of stored spermatozoa. Furthermore, their worker offspring were genotyped with DNA microsatellites to determine the number of matings and the genotypic composition of the colony. Colonies providing the male mating partners were also inferred. Both neonicotinoid and control queens mated with drones originating from the same drone source colonies, and stored similar number of spermatozoa. However, queens reared in colonies exposed to both neonicotinoids experienced fewer matings. This resulted in a reduction of the genetic diversity in their colonies (i.e. higher intracolonial relatedness). As decreased genetic diversity among worker bees is known to negatively affect colony vitality, neonicotinoids may have a cryptic effect on colony health by reducing the mating frequency of queens.

Unity in defence: honeybee workers exhibit conserved molecular responses to diverse pathogens.

BACKGROUND: Organisms typically face infection by diverse pathogens, and hosts are thought to have developed specific responses to each type of pathogen they encounter. The advent of transcriptomics now makes it possible to test this hypothesis and compare host gene expression responses to multiple pathogens at a genome-wide scale. Here, we performed a meta-analysis of multiple published and new transcriptomes using a newly developed bioinformatics approach that filters genes based on their expression profile across datasets. Thereby, we identified common and unique molecular responses of a model host species, the honey bee (Apis mellifera), to its major pathogens and parasites: the Microsporidia Nosema apis and Nosema ceranae, RNA viruses, and the ectoparasitic mite Varroa destructor, which transmits viruses. RESULTS: We identified a common suite of genes and conserved molecular pathways that respond to all investigated pathogens, a result that suggests a commonality in response mechanisms to diverse pathogens. We found that genes differentially expressed after infection exhibit a higher evolutionary rate than non-differentially expressed genes. Using our new bioinformatics approach, we unveiled additional pathogen-specific responses of honey bees; we found that apoptosis appeared to be an important response following microsporidian infection, while genes from the immune signalling pathways, Toll and Imd, were differentially expressed after Varroa/virus infection. Finally, we applied our bioinformatics approach and generated a gene co-expression network to identify highly connected (hub) genes that may represent important mediators and regulators of anti-pathogen responses. CONCLUSIONS: Our meta-analysis generated a comprehensive overview of the host metabolic and other biological processes that mediate interactions between insects and their pathogens. We identified key host genes and pathways that respond to phylogenetically diverse pathogens, representing an important source for future functional studies as well as offering new routes to identify or generate pathogen resilient honey bee stocks. The statistical and bioinformatics approaches that were developed for this study are broadly applicable to synthesize information across transcriptomic datasets. These approaches will likely have utility in addressing a variety of biological questions.

Seasonal cycle of inbreeding and recombination of the parasitic mite Varroa destructor in honeybee colonies and its implications for the selection of acaricide resistance.

Varroa destructor is the most devastating parasite of the Western honeybee, Apis mellifera. In the light of the arm race opposing the host and its parasite, the population dynamics and genetic diversity of these organisms are key parameters. However, the life cycle of V. destructor is characterized by extreme inbreeding due to full sibling mating in the host brood cells. We here present an equation reflecting the evolution of inbreeding in such a clonal system, and compare our predictions with empirical data based on the analysis of seven microsatellite markers. This comparison revealed that the mites perform essentially incestuous mating in the beginning of the brood season. However, this pattern changes with the development of mite infestation. Despite the fact that the overall level of genetic diversity of the mites remained low through the season, multiple inbred lineages were identified in the mites we sampled in June. As a response to the decrease of brood availability and the increase of the parasite population in parallel in the colonies, these lineages recombined towards the end of the season as mites co-infest brood cells. Our results suggest that the ratio of the number of mite per brood cell in the colony determines the genetic structure of the populations of V. destructor. This intracolonial population dynamics has great relevance for the selection of acaricide resistance in V. destructor. If chemical treatments occur before the recombination phase, inbreeding will greatly enhance the fixation of resistance alleles at the colony level.

Comparative analyses of the major royal jelly protein gene cluster in three Apis species with long amplicon sequencing.

The western honeybee, Apis mellifera is a prominent model organism in the field of sociogenomics and a recent upgrade substantially improved annotations of the reference genome. Nevertheless, genome assemblies based on short-sequencing reads suffer from problems in regions comprising e.g. multi-copy genes. We used single-molecule nanopore-based sequencing with extensive read-lengths to reconstruct the organization of the major royal jelly protein (mrjp) region in three species of the genus Apis. Long-amplicon sequencing provides evidence for lineage-specific evolutionary fates of Apis mrjps. Whereas the most basal species, A. florea, seems to encode ten mrjps, different patterns of gene loss and retention were observed for A. mellifera and A. dorsata. Furthermore, we show that a previously reported pseudogene in A. mellifera, mrjp2-like, is an assembly artefact arising from short read sequencing.

Food to some, poison to others - honeybee royal jelly and its growth inhibiting effect on European Foulbrood bacteria.

Honeybee colonies (Apis mellifera) serve as attractive hosts for a variety of pathogens providing optimal temperatures, humidity, and an abundance of food. Thus, honeybees have to deal with pathogens throughout their lives and, even as larvae they are affected by severe brood diseases like the European Foulbrood caused by Melissococcus plutonius. Accordingly, it is highly adaptive that larval food jelly contains antibiotic compounds. However, although food jelly is primarily consumed by bee larvae, studies investigating the antibiotic effects of this jelly have largely concentrated on bacterial human diseases. In this study, we show that royal jelly fed to queen larvae and added to the jelly of drone and worker larvae, inhibits not only the growth of European Foulbrood-associated bacteria but also its causative agent M. plutonius. This effect is shown to be caused by the main protein (major royal jelly protein 1) of royal jelly.

Differential proteomics reveals novel insights into Nosema-honey bee interactions.

Host manipulation is a common strategy by parasites to reduce host defense responses, enhance development, host exploitation, reproduction and, ultimately, transmission success. As these parasitic modifications can reduce host fitness, increased selection pressure may result in reciprocal adaptations of the host. Whereas the majority of studies on host manipulation have explored resistance against parasites (i.e. ability to prevent or limit an infection), data describing tolerance mechanisms (i.e. ability to limit harm of an infection) are scarce. By comparing differential protein abundance, we provide evidence of host-parasite interactions in the midgut proteomes of N. ceranae-infected and uninfected honey bees from both Nosema-tolerant and Nosema-sensitive lineages. We identified 16 proteins out of 661 protein spots that were differentially abundant between experimental groups. In general, infections of Nosema resulted in an up-regulation of the bee's energy metabolism. Additionally, we identified 8 proteins that were differentially abundant between tolerant and sensitive honey bees regardless of the Nosema infection. Those proteins were linked to metabolism, response to oxidative stress and apoptosis. In addition to bee proteins, we also identified 3 Nosema ceranae proteins. Interestingly, abundance of two of these Nosema proteins were significantly higher in infected Nosema-sensitive honeybees relative to the infected Nosema-tolerant lineage. This may provide a novel candidate for studying the molecular interplay between N. ceranae and its honey bee host in more detail.

Brain transcriptomes of honey bees (Apis mellifera) experimentally infected by two pathogens: Black queen cell virus and Nosema ceranae.

Regulation of gene expression in the brain plays an important role in behavioral plasticity and decision making in response to external stimuli. However, both can be severely affected by environmental factors, such as parasites and pathogens. In honey bees, the emergence and re-emergence of pathogens and potential for pathogen co-infection and interaction have been suggested as major components that significantly impaired social behavior and survival. To understand how the honey bee is affected and responds to interacting pathogens, we co-infected workers with two prevalent pathogens of different nature, the positive single strand RNA virus Black queen cell virus (BQCV), and the Microsporidia Nosema ceranae, and explored gene expression changes in brains upon single infections and co-infections. Our data provide an important resource for research on honey bee diseases, and more generally on insect host-pathogen and pathogen-pathogen interactions. Raw and processed data are publicly available in the NCBI/GEO database: (http://www.ncbi.nlm.nih.gov/geo/) under accession number GSE81664.

Viral prevalence increases with regional colony abundance in honey bee drones (Apis mellifera L).

Transmission among colonies is a central feature for the epidemiology of honey bee pathogens. High colony abundance may promote transmission among colonies independently of apiary layout, making colony abundance a potentially important parameter determining pathogen prevalence in populations of honey bees. To test this idea, we sampled male honey bees (drones) from seven distinct drone congregation areas (DCA), and used their genotypes to estimate colony abundance at each site. A multiplex ligation dependent probe amplification assay (MLPA) was used to assess the prevalence of ten viruses, using five common viral targets, in individual drones. There was a significant positive association between colony abundance and number of viral infections. This result highlights the potential importance of high colony abundance for pathogen prevalence, possibly because high population density facilitates pathogen transmission. Pathogen prevalence in drones collected from DCAs may be a useful means of estimating the disease status of a population of honey bees during the mating season, especially for localities with a large number of wild or feral colonies.

Parasite resistance and tolerance in honeybees at the individual and social level.

Organisms living in large groups, such as social insects, are particularly vulnerable to parasite transmission. However, they have evolved diverse defence mechanisms which are not only restricted to the individual's immune response, but also include social defences. Here, we review cases of adaptations at the individual and social level in the honeybee Apis mellifera against the ectoparasitic mite Varroa destructor and the endoparasitic microsporidians Nosema ceranae and Nosema apis. They are considered important threats to honeybee health worldwide. We highlight how individual resistance may result in tolerance at the colony level and vice versa.

Nosema spp. infections cause no energetic stress in tolerant honeybees.

Host-pathogen coevolution leads to reciprocal adaptations, allowing pathogens to increase host exploitation or hosts to minimise costs of infection. As pathogen resistance is often associated with considerable costs, tolerance may be an evolutionary alternative. Here, we examined the effect of two closely related and highly host dependent intracellular gut pathogens, Nosema apis and Nosema ceranae, on the energetic state in Nosema tolerant and sensitive honeybees facing the infection. We quantified the three major haemolymph carbohydrates fructose, glucose, and trehalose using high-performance liquid chromatography (HPLC) as a measure for host energetic state. Trehalose levels in the haemolymph were negatively associated with N. apis infection intensity and with N. ceranae infection regardless of the infection intensity in sensitive honeybees. Nevertheless, there was no such association in Nosema spp. infected tolerant honeybees. These findings suggest that energy availability in tolerant honeybees was not compromised by the infection. This result obtained at the individual level may also have implications at the colony level where workers in spite of a Nosema infection can still perform as well as healthy bees, maintaining colony efficiency and productivity.