Three Halloween genes from the Varroa mite, Varroa destructor (Anderson & Trueman) and their expression during reproduction

The ecdysteroid biosynthetic pathway involves sequential enzymatic hydroxylations by a group of enzymes collectively known as Halloween gene proteins. Complete sequences for three Halloween genes, spook (Vdspo), disembodied (Vddib) and shade (Vdshd), were identified in varroa mites and sequenced. Phylogenetic analyses of predicted amino acid sequences for Halloween orthologues showed that the acarine orthologues were distantly associated with insect and crustacean clades indicating that acarine genes had more ancestral characters. The lack of orthologues or pseudogenes for remaining genes suggests these pathway elements had not evolved in ancestral arthropods. Vdspo transcript levels were highest in gut tissues, while Vddib transcript levels were highest in ovary-lyrate organs. In contrast, Vdshd transcript levels were lower overall but present in both gut and ovary-lyrate organs. All three transcripts were present in eggs removed from gravid female mites. A brood cell invasion assay was developed for acquiring synchronously staged mites. Mites within 4 h of entering a brood cell had transcript levels of all three that were not significantly different from mites on adult bees. These analyses suggest that varroa mites may be capable of modifying 7-dehydro-cholesterol precursor and hydroxylations of other steroid precursors, but whether the mites directly produce ecdysteroid precursors and products remains undetermined.

Characterization of gut bacteria at different developmental stages of Asian honey bees, Apis cerana

Previous surveys have shown that adult workers of the Asian honey bee Apis cerana harbor four major gut microbes (Bifidobacterium, Snodgrassella alvi, Gilliamella apicola, and Lactobacillus). Using quantitative PCR we characterized gut bacterial communities across the life cycle of A. cerana from larvae to workers. Our results indicate that the presence and quantity of these four bacteria were low on day 1, increased rapidly after day 5, and then peaked during days 10-20. They stabilized from days 20-25 or days 25-30, then dropped to a low level at day 30. In addition, the larvae infected by Sacbrood virus or European foulbrood had significantly lower copies of 16S rRNA genes than healthy individuals.

Characterization of Two Species of Trypanosomatidae from the Honey Bee Apis mellifera: Crithidia mellificae Langridge and McGhee, and Lotmaria passim n. gen., n. sp

Trypanosomatids are increasingly recognized as prevalent in European honey bees (Apis mellifera) and by default are attributed to one recognized species, Crithidia mellificae Langridge and McGhee, 1967. We provide reference genetic and ultrastructural data for type isolates of C. mellificae (ATCC 30254 and 30862) in comparison with two recent isolates from A. mellifera (BRL and SF). Phylogenetics unambiguously identify strains BRL/SF as a novel taxonomic unit distinct from C. mellificae strains 30254/30862 and assign all four strains as lineages of a novel clade within the subfamily Leishmaniinae. In vivo analyses show strains BRL/SF preferably colonize the hindgut, lining the lumen as adherent spheroids in a manner identical to previous descriptions from C. mellificae. Microscopy images show motile forms of C. mellificae are distinct from strains BRL/SF. We propose the binomial Lotmaria passim n. gen., n. sp. for this previously undescribed taxon. Analyses of new and previously accessioned genetic data show C. mellificae is still extant in bee populations, however, L. passim n. gen., n. sp. is currently the predominant trypanosomatid in A. mellifera globally. Our findings require that previous reports of C. mellificae be reconsidered and that subsequent trypanosomatid species designations from Hymenoptera provide genetic support.

Metatranscriptomic analyses of honey bee colonies

Honey bees face numerous biotic threats from viruses to bacteria, fungi, protists, and mites. Here we describe a thorough analysis of microbes harbored by worker honey bees collected from field colonies in geographically distinct regions of Turkey. Turkey is one of the World's most important centers of apiculture, harboring five subspecies of Apis mellifera L., approximately 20% of the honey bee subspecies in the world. We use deep ILLUMINA-based RNA sequencing to capture RNA species for the honey bee and a sampling of all non-endogenous species carried by bees. After trimming and mapping these reads to the honey bee genome, approximately 10% of the sequences (9–10 million reads per library) remained. These were then mapped to a curated set of public sequences containing ca. Sixty megabase-pairs of sequence representing known microbial species associated with honey bees. Levels of key honey bee pathogens were confirmed using quantitative PCR screens. We contrast microbial matches across different sites in Turkey, showing new country recordings of Lake Sinai virus, two Spiroplasma bacterium species, symbionts Candidatus Schmidhempelia bombi, Frischella perrara, Snodgrassella alvi, Gilliamella apicola, Lactobacillus spp.), neogregarines, and a trypanosome species. By using metagenomic analysis, this study also reveals deep molecular evidence for the presence of bacterial pathogens (Melissococcus plutonius, Paenibacillus larvae), Varroa destructor-1 virus, Sacbrood virus, and fungi. Despite this effort we did not detect KBV, SBPV, Tobacco ringspot virus, VdMLV (Varroa Macula like virus), Acarapis spp., Tropilaeleps spp. and Apocephalus (phorid fly). We discuss possible impacts of management practices and honey bee subspecies on microbial retinues. The described workflow and curated microbial database will be generally useful for microbial surveys of healthy and declining honey bees.

Effects of different vaccine combinations against Mycoplasma gallisepticum on the internal egg and eggshell characteristics of commercial layer chickens 1,2,3

Live F-strain Mycoplasma gallisepticum (FMG) vaccines are presently being used to help control field-strain MG outbreaks. However, they may exert some adverse effects on egg production. Live strains of MG of lesser virulence as well as killed vaccines have little or no effect on egg production, but afford lower levels of protection. This has led to research investigating their use in combination with a subsequent overlay vaccination of FMG given later in the production cycle. In the present study, 2 trials were conducted to investigate the effects of prelay vaccinations of live and killed MG vaccines or their combination, in conjunction with an FMG vaccine overlay after peak production, on the egg characteristics of commercial layers. The following vaccination treatments were administered at 10 wk of age (woa): 1) unvaccinated (Control), 2) MG-Bacterin (MGBac) vaccine, 3) ts-11 strain MG (ts11MG) vaccine, and 4) MGBac and ts11MG combination (MGBac + ts11MG). At 45 woa, half of the birds were overlaid with an FMG vaccine. In each trial, internal egg and eggshell parameters including egg weight (EW), Haugh unit score (HU), eggshell breaking strength (EBS), percentage yolk weight (PYW), percentage albumen weight (PAW), percentage eggshell weight (PSW), eggshell weight per unit surface area (SWUSA), percentage yolk moisture (PYM), and percent total lipids (PYL) were determined at various time periods between 21 and 52 woa. At 28 woa, SWUSA was lower in the ts11MG and MGBac + ts11MG groups compared to the Control group. Conversely, at 43 woa, SWUSA was higher in the ts11MG than in the MGBac group. Between 23 and 43 woa, PYL was higher in the MGBac and ts11MG groups in comparison to the Control group. In conclusion, vaccination with MGBac alone or in combination with ts11MG at 10 woa with or without an FMG vaccine overlay at 45 woa does not adversely affect the internal egg or eggshell quality of commercial layers throughout lay.

Israeli acute paralysis virus: epidemiology, pathogenesis and implications for honey bee health

Israeli acute paralysis virus (IAPV) is a widespread RNA virus of honey bees that has been linked with colony losses. Here we describe the transmission, prevalence, and genetic traits of this virus, along with host transcriptional responses to infections. Further, we present RNAi-based strategies for limiting an important mechanism used by IAPV to subvert host defenses. Our study shows that IAPV is established as a persistent infection in honey bee populations, likely enabled by both horizontal and vertical transmission pathways. The phenotypic differences in pathology among different strains of IAPV found globally may be due to high levels of standing genetic variation. Microarray profiles of host responses to IAPV infection revealed that mitochondrial function is the most significantly affected biological process, suggesting that viral infection causes significant disturbance in energy-related host processes. The expression of genes involved in immune pathways in adult bees indicates that IAPV infection triggers active immune responses. The evidence that silencing an IAPV-encoded putative suppressor of RNAi reduces IAPV replication suggests a functional assignment for a particular genomic region of IAPV and closely related viruses from the Family Dicistroviridae, and indicates a novel therapeutic strategy for limiting multiple honey bee viruses simultaneously and reducing colony losses due to viral diseases. We believe that the knowledge and insights gained from this study will provide a new platform for continuing studies of the IAPV–host interactions and have positive implications for disease management that will lead to mitigation of escalating honey bee colony losses worldwide.

The mysterious outbreak of honey bee Colony Collapse Disorder (CCD) in the US in 2006–2007 has attracted massive media attention and created great concerns over the effects of various risk factors on bee health. Understanding the factors that are linked to the honey bee colony declines may provide insights for managing similar incidents in the future. We conducted this study to elucidate traits of a key honey bee virus, Israeli acute paralysis virus. We then developed an innovative strategy to control virus levels. The knowledge and insights gained from this study will have positive implications for bee disease management, helping to mitigate worldwide colony losses.

Effects of live and killed vaccines against Mycoplasma gallisepticum on the performance characteristics of commercial layer chickens

Different vaccine strains of Mycoplasma gallisepticum have been used on multiple-age commercial layer farms in an effort to protect birds against virulent field-strain infections. Use of the F-strain of M. gallisepticum (FMG), as an overlay vaccine during lay, may be necessary because of the lower level of protection afforded by M. gallisepticum vaccines of low virulence given before lay. Two replicate trials were conducted to investigate effects of live and killed M. gallisepticum vaccines administered individually and in combination before lay, in conjunction with an FMG vaccine overlay after peak egg production (EP), on the performance characteristics of commercial layers. The following treatments were utilized at 10 wk of age (woa): 1) control (no vaccinations); 2) ts11 strain M. gallisepticum (ts11MG) vaccine; 3) M. gallisepticum-Bacterin vaccine (MG-Bacterin); and 4) ts11MG and MG-Bacterin vaccines combination. At 45 woa, half of the birds were overlaid with an FMG vaccine. Hen mortality, BW, egg weight, percentage hen-day EP, egg blood spots, and egg meat spots were determined at various time periods between 18 and 52 woa. The data from each trial were pooled. Treatment did not affect performance in interval I (23 to 45 woa). However, during interval II (46 to 52 woa), the EP of control and MG-Bacterin-vaccinated birds that later received an FMG vaccine overlay was lower than that in the other treatment groups. Furthermore, treatment application reduced bird BW during interval II. Despite the effects on BW and EP, no differences were observed for egg blood or meat spots among the various treatments. It is suggested that the vaccination of commercial layers before lay with ts11MG, but not MG-Bacterin, may reduce the negative impacts of an FMG overlay vaccination given during lay. These results establish that the vaccination of pullets with ts11MG in combination with the vaccination of hens with an FMG overlay, for continual protection against field-strain M. gallisepticum infections, may be used without suppressing performance.

Honey bee colonies act as reservoirs for two Spiroplasma facultative symbionts and incur complex, multiyear infection dynamics

Two species of Spiroplasma (Mollicutes) bacteria were isolated from and described as pathogens of the European honey bee, Apis mellifera, ~30 years ago but recent information on them is lacking despite global concern to understand bee population declines. Here we provide a comprehensive survey for the prevalence of these two Spiroplasma species in current populations of honey bees using improved molecular diagnostic techniques to assay multiyear colony samples from North America (U.S.A.) and South America (Brazil). Significant annual and seasonal fluctuations of Spiroplasma apis and Spiroplasma melliferum prevalence in colonies from the U.S.A. (n = 616) and Brazil (n = 139) occurred during surveys from 2011 through 2013. Overall, 33% of U.S.A. colonies and 54% of Brazil colonies were infected by Spiroplasma spp., where S. melliferum predominated over S. apis in both countries (25% vs. 14% and 44% vs. 38% frequency, respectively). Colonies were co-infected by both species more frequently than expected in both countries and at a much higher rate in Brazil (52%) compared to the U.S.A. (16.5%). U.S.A. samples showed that both species were prevalent not only during spring, as expected from prior research, but also during other seasons. These findings demonstrate that the model of honey bee spiroplasmas as springtime-restricted pathogens needs to be broadened and their role as occasional pathogens considered in current contexts.

Examining the role of foraging and malvolio in host-finding behavior in the honey bee parasite, Varroa destructor (Anderson & Trueman)

When a female varroa mite, Varroa destructor (Anderson & Trueman), invades a honey bee brood cell, the physiology rapidly changes from feeding phoretic to reproductive. Changes in foraging and malvolio transcript levels in the brain have been associated with modulated intra-specific food searching behaviors in insects and other invertebrates. Transcription profiles for both genes were examined during and immediately following brood cell invasion to assess their role as potential control elements. Vdfor and Vdmvl transcripts were found in all organs of varroa mites with the highest Vdfor transcript levels in ovary-lyrate organs and the highest Vdmvl in Malpighian tubules. Changes in transcript levels of Vdfor and Vdmvl in synganglia were not associated with the cell invasion process, remaining comparable between early reproductive mites (collected from the pre-capping brood cells) and phoretic mites. However, Vdfor and Vdmvl transcript levels were lowered by 37 and 53%, respectively, in synganglia from reproductive mites compared to early reproductive mites, but not significantly different to levels in synganglia from phoretic mites. On the other hand, in whole body preparations the Vdfor and Vdmvl had significantly higher levels of transcript in reproductive mites compared to phoretic and early reproductive, mainly due to the presence of both transcripts accumulating in the eggs carried by the ovipositing mite. Varroa mites are a critical component for honey bee population decline and finding varroa mite genes associated with brood cell invasion, reproduction, ion balance and other physiological processes will facilitate development of novel control avenues for this honey bee parasite.

Evaluation of cage designs and feeding regimes for honey bee (Hymenoptera: Apidae) laboratory experiments

The aim of this study was to improve cage systems for maintaining adult honey bee (Apis mellifera L.) workers under in vitro laboratory conditions. To achieve this goal, we experimentally evaluated the impact of different cages, developed by scientists of the international research network COLOSS (Prevention of honey bee COlony LOSSes), on the physiology and survival of honey bees. We identified three cages that promoted good survival of honey bees. The bees from cages that exhibited greater survival had relatively lower titers of deformed wing virus, suggesting that deformed wing virus is a significant marker reflecting stress level and health status of the host. We also determined that a leak- and drip-proof feeder was an integral part of a cage system and a feeder modified from a 20-ml plastic syringe displayed the best result in providing steady food supply to bees. Finally, we also demonstrated that the addition of protein to the bees' diet could significantly increase the level ofvitellogenin gene expression and improve bees' survival. This international collaborative study represents a critical step toward improvement of cage designs and feeding regimes for honey bee laboratory experiments.

Expression of insulin/insulin-like signalling and TOR pathway genes in honey bee caste determination

The development of queen and worker castes in honey bees is induced by differential nutrition, with future queens and workers receiving diets that are qualitatively and quantitatively different. We monitored the gene expression of 14 genes for components of the insulin/insulin-like signalling and TOR pathways in honey bee larvae from 40-88 h after hatching. We compared normally fed queen and normally fed worker larvae and found that three genes showed expression differences in 40-h-old larvae. Genes that show such early differences in expression may be part of the mechanism that transduces nutrition level into a hormone signal. We then compared changes in expression after shifts in diet with those in normally developing queens and workers. Following a shift to the worker diet, the expression of 9/14 genes was upregulated in comparison with queens. Following a shift to the queen diet, expression of only one gene changed. The honey bee responses may function together as a homeostatic mechanism buffering larvae from caste-disrupting variation in nutrition. The different responses would be part of the canalization of both the queen and worker developmental pathways, and as such, a signature of advanced sociality.

Finding the missing honey bee genes: lessons learned from a genome upgrade

Background

The first generation of genome sequence assemblies and annotations have had a significant impact upon our understanding of the biology of the sequenced species, the phylogenetic relationships among species, the study of populations within and across species, and have informed the biology of humans. As only a few Metazoan genomes are approaching finished quality (human, mouse, fly and worm), there is room for improvement of most genome assemblies. The honey bee (Apis mellifera) genome, published in 2006, was noted for its bimodal GC content distribution that affected the quality of the assembly in some regions and for fewer genes in the initial gene set (OGSv1.0) compared to what would be expected based on other sequenced insect genomes.

Results

Here, we report an improved honey bee genome assembly (Amel_4.5) with a new gene annotation set (OGSv3.2), and show that the honey bee genome contains a number of genes similar to that of other insect genomes, contrary to what was suggested in OGSv1.0. The new genome assembly is more contiguous and complete and the new gene set includes ~5000 more protein-coding genes, 50% more than previously reported. About 1/6 of the additional genes were due to improvements to the assembly, and the remaining were inferred based on new RNAseq and protein data.

Conclusions

Lessons learned from this genome upgrade have important implications for future genome sequencing projects. Furthermore, the improvements significantly enhance genomic resources for the honey bee, a key model for social behavior and essential to global ecology through pollination.

Genomic organization and reproductive regulation of a large lipid transfer protein in the varroa mite, Varroa destructor (Anderson & Trueman)

The complete genomic region and corresponding transcript of the most abundant protein in phoretic varroa mites, Varroa destructor (Anderson & Trueman), were sequenced and have homology with acarine hemelipoglycoproteins and the large lipid transfer protein (LLTP) super family. The genomic sequence of VdLLTP included 14 introns and the mature transcript coded for a predicted polypeptide of 1575 amino acid residues. VdLLTP shared a minimum of 25% sequence identity with acarine LLTPs. Phylogenetic assessment showed VdLLTP was most closely related to Metaseiulus occidentalis vitellogenin and LLTP proteins of ticks; however, no heme binding by VdLLTP was detected. Analysis of lipids associated with VdLLTP showed that it was a carrier for free and esterified C12 -C22 fatty acids from triglycerides, diacylglycerides and monoacylglycerides. Additionally, cholesterol and β-sitosterol were found as cholesterol esters linked to common fatty acids. Transcript levels of VdLLTP were 42 and 310 times higher in phoretic female mites when compared with males and quiescent deutonymphs, respectively. Coincident with initiation of the reproductive phase, VdLLTP transcript levels declined to a third of those in phoretic female mites. VdLLTP functions as an important lipid transporter and should provide a significant RNA interference target for assessing the control of varroa mites.

In vitro infection of pupae with Israeli acute paralysis virus suggests disturbance of transcriptional homeostasis in honey bees (Apis mellifera)

The ongoing decline of honey bee health worldwide is a serious economic and ecological concern. One major contributor to the decline are pathogens, including several honey bee viruses. However, information is limited on the biology of bee viruses and molecular interactions with their hosts. An experimental protocol to test these systems was developed, using injections of Israeli Acute Paralysis Virus (IAPV) into honey bee pupae reared ex-situ under laboratory conditions. The infected pupae developed pronounced but variable patterns of disease. Symptoms varied from complete cessation of development with no visual evidence of disease to rapid darkening of a part or the entire body. Considerable differences in IAPV titer dynamics were observed, suggesting significant variation in resistance to IAPV among and possibly within honey bee colonies. Thus, selective breeding for virus resistance should be possible. Gene expression analyses of three separate experiments suggest IAPV disruption of transcriptional homeostasis of several fundamental cellular functions, including an up-regulation of the ribosomal biogenesis pathway. These results provide first insights into the mechanisms of IAPV pathogenicity. They mirror a transcriptional survey of honey bees afflicted with Colony Collapse Disorder and thus support the hypothesis that viruses play a critical role in declining honey bee health.

Single and mixed-species trypanosome and microsporidia infections elicit distinct, ephemeral cellular and humoral immune responses in honey bees

Frequently encountered parasite species impart strong selective pressures on host immune system evolution and are more apt to concurrently infect the same host, yet molecular impacts in light of this are often overlooked. We have contrasted immune responses in honey bees to two common eukaryotic endoparasites by establishing single and mixed-species infections using the long-associated parasite Crithidia mellificae and the emergent parasite Nosema ceranae. Quantitative polymerase chain reaction was used to screen host immune gene expression at 9 time points post inoculation. Systemic responses in abdomens during early stages of parasite establishment revealed conserved receptor (Down syndrome cell adhesion molecule, Dscam and nimrod C1, nimC1), signaling (MyD88 and Imd) and antimicrobial peptide (AMP) effector (Defensin 2) responses. Late, established infections were distinct with a refined 2 AMP response to C. mellificae that contrasted starkly with a 5 AMP response to N. ceranae. Mixed species infections induced a moderate 3 AMPs. Transcription in gut tissues highlighted important local roles for Dscam toward both parasites and Imd signaling toward N. ceranae. At both systemic and local levels Dscam, MyD88 and Imd transcription was consistently correlated based on clustering analysis. Significant gene suppression occurred in two cases from midgut to ileum tissue: Dscam was lowered during mixed infections compared to N. ceranae infections and both C. mellificae and mixed infections had reduced nimC1 transcription compared to uninfected controls. We show that honey bees rapidly mount complex immune responses to both Nosema and Crithidia that are dynamic over time and that mixed-species infections significantly alter local and systemic immune gene transcription.

Transcriptional response of honey bee larvae infected with the bacterial pathogen Paenibacillus larvae

American foulbrood disease of honey bees is caused by the bacterium Paenibacillus larvae. Infection occurs per os in larvae and systemic infection requires a breaching of the host peritrophic matrix and midgut epithelium. Genetic variation exists for both bacterial virulence and host resistance, and a general immunity is achieved by larvae as they age, the basis of which has not been identified. To quickly identify a pool of candidate genes responsive to P. larvae infection, we sequenced transcripts from larvae inoculated with P. larvae at 12 hours post-emergence and incubated for 72 hours, and compared expression levels to a control cohort. We identified 75 genes with significantly higher expression and six genes with significantly lower expression. In addition to several antimicrobial peptides, two genes encoding peritrophic-matrix domains were also up-regulated. Extracellular matrix proteins, proteases/protease inhibitors, and members of the Osiris gene family were prevalent among differentially regulated genes. However, analysis of Drosophila homologs of differentially expressed genes revealed spatial and temporal patterns consistent with developmental asynchrony as a likely confounder of our results. We therefore used qPCR to measure the consistency of gene expression changes for a subset of differentially expressed genes. A replicate experiment sampled at both 48 and 72 hours post infection allowed further discrimination of genes likely to be involved in host response. The consistently responsive genes in our test set included a hymenopteran-specific protein tyrosine kinase, a hymenopteran specific serine endopeptidase, a cytochrome P450 (CYP9Q1), and a homolog of trynity, a zona pellucida domain protein. Of the known honey bee antimicrobial peptides, apidaecin was responsive at both time-points studied whereas hymenoptaecin was more consistent in its level of change between biological replicates and had the greatest increase in expression by RNA-seq analysis.

Population-genomic variation within RNA viruses of the Western honey bee, Apis mellifera, inferred from deep sequencing

Background

Deep sequencing of viruses isolated from infected hosts is an efficient way to measure population-genetic variation and can reveal patterns of dispersal and natural selection. In this study, we mined existing Illumina sequence reads to investigate single-nucleotide polymorphisms (SNPs) within two RNA viruses of the Western honey bee (Apis mellifera), deformed wing virus (DWV) and Israel acute paralysis virus (IAPV). All viral RNA was extracted from North American samples of honey bees or, in one case, the ectoparasitic mite Varroa destructor.

Results

Coverage depth was generally lower for IAPV than DWV, and marked gaps in coverage occurred in several narrow regions (< 50 bp) of IAPV. These coverage gaps occurred across sequencing runs and were virtually unchanged when reads were re-mapped with greater permissiveness (up to 8% divergence), suggesting a recurrent sequencing artifact rather than strain divergence. Consensus sequences of DWV for each sample showed little phylogenetic divergence, low nucleotide diversity, and strongly negative values of Fu and Li’s D statistic, suggesting a recent population bottleneck and/or purifying selection. The Kakugo strain of DWV fell outside of all other DWV sequences at 100% bootstrap support. IAPV consensus sequences supported the existence of multiple clades as had been previously reported, and Fu and Li’s D was closer to neutral expectation overall, although a sliding-window analysis identified a significantly positive D within the protease region, suggesting selection maintains diversity in that region. Within-sample mean diversity was comparable between the two viruses on average, although for both viruses there was substantial variation among samples in mean diversity at third codon positions and in the number of high-diversity sites. F_ST values were bimodal for DWV, likely reflecting neutral divergence in two low-diversity populations, whereas IAPV had several sites that were strong outliers with very low F_ST.

Conclusions

This initial survey of genetic variation within honey bee RNA viruses suggests future directions for studies examining the underlying causes of population-genetic structure in these economically important pathogens.

Variable induction of vitellogenin genes in the varroa mite, Varroa destructor (Anderson & Trueman), by the honeybee, Apis mellifera L, host and its environment

Transcript levels of vitellogenins (Vgs) in the varroa mite, Varroa destructor (Anderson & Trueman), were variably induced by interactions between the developing honeybee, Apis mellifera L, as a food source and the capped honeybee cell environment. Transcripts for two Vgs of varroa mites were sequenced and putative Vg protein products characterized. Sequence analysis of VdVg1 and VdVg2 proteins showed that each had greater similarity with Vg1 and Vg2 proteins from ticks, respectively, than between themselves and were grouped separately by phylogenetic analyses. This suggests there was a duplication of the ancestral acarine Vg gene prior to the divergence of the mites and ticks. Low levels of transcript were detected in immature mites, males and phoretic females. Following cell invasion by phoretic females, VdVg1 and VdVg2 transcript levels were up-regulated after cell capping to a maximum at the time of partial cocoon formation by the honeybee. During oviposition the two transcripts were differentially expressed with higher levels of VdVg2 being observed. A bioassay based on assessing the transcript levels was established. Increases in VdVg1 and VdVg2 transcripts were induced experimentally in phoretic females when they were placed inside a cell containing an early metamorphosing last instar bee but not when exposed to the metamorphosing bee alone. The variable response of Vg expression to the food source as well as environmental cues within the capped cell demonstrates that perturbation of host-parasite interactions may provide avenues to disrupt the reproductive cycle of the varroa mites and prevent varroasis.

The prevalence of parasites and pathogens in Asian honeybees Apis cerana in China

Pathogens and parasites represent significant threats to the health and well-being of honeybee species that are key pollinators of agricultural crops and flowers worldwide. We conducted a nationwide survey to determine the occurrence and prevalence of pathogens and parasites in Asian honeybees, Apis cerana, in China. Our study provides evidence of infections of A. cerana by pathogenic Deformed wing virus (DWV), Black queen cell virus (BQCV), Nosema ceranae, and C. bombi species that have been linked to population declines of European honeybees, A. mellifera, and bumble bees. However, the prevalence of DWV, a virus that causes widespread infection in A. mellifera, was low, arguably a result of the greater ability of A. cerana to resist the ectoprasitic mite Varroa destructor, an efficient vector of DWV. Analyses of microbial communities from the A. cerana digestive tract showed that Nosema infection could have detrimental effects on the gut microbiota. Workers infected by N. ceranae tended to have lower bacterial quantities, with these differences being significant for the Bifidobacterium and Pasteurellaceae bacteria groups. The results of this nationwide screen show that parasites and pathogens that have caused serious problems in European honeybees can be found in native honeybee species kept in Asia. Environmental changes due to new agricultural practices and globalization may facilitate the spread of pathogens into new geographic areas. The foraging behavior of pollinators that are in close geographic proximity likely have played an important role in spreading of parasites and pathogens over to new hosts. Phylogenetic analyses provide insights into the movement and population structure of these parasites, suggesting a bidirectional flow of parasites among pollinators. The presence of these parasites and pathogens may have considerable implications for an observed population decline of Asian honeybees.

Differential expression of immune genes of adult honey bee (Apis mellifera) after inoculated by Nosema ceranae

Nosema ceranae is a microsporidium parasite infecting adult honey bees (Apis mellifera) and is known to affects at both the individual and colony level. In this study, the expression levels were measured for four antimicrobial peptide encoding genes that are associated with bee humoral immunity (defensin, abaecin, apidaecin, and hymenoptaecin), eater gene which is a transmembrane protein involved cellular immunity and gene encoding female-specific protein (vitellogenin) in honey bees when inoculated by N. ceranae. The results showed that four of these genes, defensin, abaecin, apidaecin and hymenoptaecin were significantly down-regulated 3 and 6days after inoculations. Additionally, antimicrobial peptide expressions did not significantly differ between control and inoculated bees after 12days post inoculation. Moreover, our results revealed that the mRNA levels of eater and vitellogenin did not differ significantly following N. ceranae inoculation. Therefore, in this study we reaffirmed that N. ceranae infection induces host immunosuppression.

Gene expression in honey bee (Apis mellifera) larvae exposed to pesticides and Varroa mites (Varroa destructor)

Honey bee (Apis mellifera) larvae reared in vitro were exposed to one of nine pesticides and/or were challenged with the parasitic mite, Varroa destructor. Total RNA was extracted from individual larvae and first strand cDNAs were generated. Gene-expression changes in larvae were measured using quantitative PCR (qPCR) targeting transcripts for pathogens and genes involved in physiological processes, bee health, immunity, and/or xenobiotic detoxification. Transcript levels for Peptidoglycan Recognition Protein (PGRPSC), a pathogen recognition gene, increased in larvae exposed to Varroa mites (P<0.001) and were not changed in pesticide treated larvae. As expected, Varroa-parasitized brood had higher transcripts of Deformed Wing Virus than did control larvae (P<0.001). Varroa parasitism, arguably coupled with virus infection, resulted in significantly higher transcript abundances for the antimicrobial peptides abaecin, hymenoptaecin, and defensin1. Transcript levels for Prophenoloxidase-activating enzyme (PPOact), an immune end product, were elevated in larvae treated with myclobutanil and chlorothalonil (both are fungicides) (P<0.001). Transcript levels for Hexameric storage protein (Hsp70) were significantly upregulated in imidacloprid, fluvalinate, coumaphos, myclobutanil, and amitraz treated larvae. Definitive impacts of pesticides and Varroa parasitism on honey bee larval gene expression were demonstrated. Interactions between larval treatments and gene expression for the targeted genes are discussed.

Transcriptome analysis of the honey bee fungal pathogen, Ascosphaera apis: implications for host pathogenesis

Background

We present a comprehensive transcriptome analysis of the fungus Ascosphaera apis, an economically important pathogen of the Western honey bee (Apis mellifera) that causes chalkbrood disease. Our goals were to further annotate the A. apis reference genome and to identify genes that are candidates for being differentially expressed during host infection versus axenic culture.

Results

We compared A. apis transcriptome sequence from mycelia grown on liquid or solid media with that dissected from host-infected tissue. 454 pyrosequencing provided 252 Mb of filtered sequence reads from both culture types that were assembled into 10,087 contigs. Transcript contigs, protein sequences from multiple fungal species, and ab initio gene predictions were included as evidence sources in the Maker gene prediction pipeline, resulting in 6,992 consensus gene models. A phylogeny based on 12 of these protein-coding loci further supported the taxonomic placement of Ascosphaera as sister to the core Onygenales. Several common protein domains were less abundant in A. apis compared with related ascomycete genomes, particularly cytochrome p450 and protein kinase domains. A novel gene family was identified that has expanded in some ascomycete lineages, but not others. We manually annotated genes with homologs in other fungal genomes that have known relevance to fungal virulence and life history. Functional categories of interest included genes involved in mating-type specification, intracellular signal transduction, and stress response. Computational and manual annotations have been made publicly available on the Bee Pests and Pathogens website.

Conclusions

This comprehensive transcriptome analysis substantially enhances our understanding of the A. apis genome and its expression during infection of honey bee larvae. It also provides resources for future molecular studies of chalkbrood disease and ultimately improved disease management.

Direct effect of acaricides on pathogen loads and gene expression levels in honey bees Apis mellifera

The effect of using acaricides to control varroa mites has long been a concern to the beekeeping industry due to unintended negative impacts on honey bee health. Irregular ontogenesis, suppression of immune defenses, and impairment of normal behavior have been linked to pesticide use. External stressors, including parasites and the pathogens they vector, can confound studies on the effects of pesticides on the metabolism of honey bees. This is the case of Varroa destructor, a mite that negatively affects honey bee health on many levels, from direct parasitism, which diminishes honey bee productivity, to vectoring and/or activating other pathogens, including many viruses. Here we present a gene expression profile comprising genes acting on diverse metabolic levels (detoxification, immunity, and development) in a honey bee population that lacks the influence of varroa mites. We present data for hives treated with five different acaricides; Apiguard (thymol), Apistan (tau-fluvalinate), Checkmite (coumaphos), Miteaway (formic acid) and ApiVar (amitraz). The results indicate that thymol, coumaphos and formic acid are able to alter some metabolic responses. These include detoxification gene expression pathways, components of the immune system responsible for cellular response and the c-Jun amino-terminal kinase (JNK) pathway, and developmental genes. These could potentially interfere with the health of individual honey bees and entire colonies.

Predictive markers of honey bee colony collapse

Across the Northern hemisphere, managed honey bee colonies, Apis mellifera, are currently affected by abrupt depopulation during winter and many factors are suspected to be involved, either alone or in combination. Parasites and pathogens are considered as principal actors, in particular the ectoparasitic mite Varroa destructor, associated viruses and the microsporidian Nosema ceranae. Here we used long term monitoring of colonies and screening for eleven disease agents and genes involved in bee immunity and physiology to identify predictive markers of honeybee colony losses during winter. The data show that DWV, Nosema ceranae, Varroa destructor and Vitellogenin can be predictive markers for winter colony losses, but their predictive power strongly depends on the season. In particular, the data support that V. destructor is a key player for losses, arguably in line with its specific impact on the health of individual bees and colonies.