Immune pathways and defence mechanisms in honey bees Apis mellifera

The genomes of two key bumblebee species with primitive eusocial organization

Background

The shift from solitary to social behavior is one of the major evolutionary transitions. Primitively eusocial bumblebees are uniquely placed to illuminate the evolution of highly eusocial insect societies. Bumblebees are also invaluable natural and agricultural pollinators, and there is widespread concern over recent population declines in some species. High-quality genomic data will inform key aspects of bumblebee biology, including susceptibility to implicated population viability threats.

Results

We report the high quality draft genome sequences of Bombus terrestris and Bombus impatiens, two ecologically dominant bumblebees and widely utilized study species. Comparing these new genomes to those of the highly eusocial honeybee Apis mellifera and other Hymenoptera, we identify deeply conserved similarities, as well as novelties key to the biology of these organisms. Some honeybee genome features thought to underpin advanced eusociality are also present in bumblebees, indicating an earlier evolution in the bee lineage. Xenobiotic detoxification and immune genes are similarly depauperate in bumblebees and honeybees, and multiple categories of genes linked to social organization, including development and behavior, show high conservation. Key differences identified include a bias in bumblebee chemoreception towards gustation from olfaction, and striking differences in microRNAs, potentially responsible for gene regulation underlying social and other traits.

Conclusions

These two bumblebee genomes provide a foundation for post-genomic research on these key pollinators and insect societies. Overall, gene repertoires suggest that the route to advanced eusociality in bees was mediated by many small changes in many genes and processes, and not by notable expansion or depauperation.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0623-3) contains supplementary material, which is available to authorized users.

A depauperate immune repertoire precedes evolution of sociality in bees

Background

Sociality has many rewards, but can also be dangerous, as high population density and low genetic diversity, common in social insects, is ideal for parasite transmission. Despite this risk, honeybees and other sequenced social insects have far fewer canonical immune genes relative to solitary insects. Social protection from infection, including behavioral responses, may explain this depauperate immune repertoire. Here, based on full genome sequences, we describe the immune repertoire of two ecologically and commercially important bumblebee species that diverged approximately 18 million years ago, the North American Bombus impatiens and European Bombus terrestris.

Results

We find that the immune systems of these bumblebees, two species of honeybee, and a solitary leafcutting bee, are strikingly similar. Transcriptional assays confirm the expression of many of these genes in an immunological context and more strongly in young queens than males, affirming Bateman’s principle of greater investment in female immunity. We find evidence of positive selection in genes encoding antiviral responses, components of the Toll and JAK/STAT pathways, and serine protease inhibitors in both social and solitary bees. Finally, we detect many genes across pathways that differ in selection between bumblebees and honeybees, or between the social and solitary clades.

Conclusions

The similarity in immune complement across a gradient of sociality suggests that a reduced immune repertoire predates the evolution of sociality in bees. The differences in selection on immune genes likely reflect divergent pressures exerted by parasites across social contexts.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0628-y) contains supplementary material, which is available to authorized users.

High throughput profiling of the cotton bollworm Helicoverpa armigera immunotranscriptome during the fungal and bacterial infections

Background

Innate immunity is essential in defending against invading pathogens in invertebrates. The cotton bollworm, Helicoverpa armigera (Hübner) is one of the most destructive lepidopteran pests, which causes enormous economic losses in agricultural production worldwide. The components of the immune system are largely unknown in this insect. The application of entomopathogens is considered as an alternative to the chemical insecticides for its control. However, few studies have focused on the molecular mechanisms of host-pathogen interactions between pest insects and their pathogens. Here, we investigated the immunotranscriptome of H. armigera larvae and examined gene expression changes after pathogen infections. This study provided insights into the potential immunity-related genes and pathways in H. armigera larvae.

Results

Here, we adopted a high throughput RNA-seq approach to determine the immunotranscriptome of H. armigera larvae injected with buffer, fungal pathogen Beauveria bassiana, or Gram-negative bacterium Enterobacter cloacae. Based on sequence similarity to those homologs known to participate in immune responses in other insects, we identified immunity-related genes encoding pattern recognition receptors, signal modulators, immune effectors, and nearly all members of the Toll, IMD and JAK/STAT pathways. The RNA-seq data indicated that some immunity-related genes were activated in fungus- and bacterium-challenged fat body while others were suppressed in B. bassiana challenged hemocytes, including the putative IMD and JAK-STAT pathway members. Bacterial infection elevated the expression of recognition and modulator genes in the fat body and signal pathway genes in hemocytes. Although fat body and hemocytes both are important organs involved in the immune response, our transcriptome analysis revealed that more immunity-related genes were induced in the fat body than that hemocytes. Furthermore, quantitative real-time PCR analysis confirmed that, consistent with the RNA-seq data, the transcript abundances of putative PGRP-SA1, Serpin1, Toll-14, and Spz2 genes were elevated in fat body upon B. bassiana infection, while the mRNA levels of defensin, moricin1, and gloverin1 were up-regulated in hemocytes.

Conclusions

In this study, a global survey of the host defense against fungal and bacterial infection was performed on the non-model lepidopteran pest species. The comprehensive sequence resource and expression profiles of the immunity-related genes in H. armigera are acquired. This study provided valuable information for future functional investigations as well as development of specific and effective agents to control this pest.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1509-1) contains supplementary material, which is available to authorized users.

Honey glycoproteins containing antimicrobial peptides, Jelleins of the Major Royal Jelly Protein 1, are responsible for the cell wall lytic and bactericidal activities of honey

We have recently identified the bacterial cell wall as the cellular target for honey antibacterial compounds; however, the chemical nature of these compounds remained to be elucidated. Using Concavalin A-affinity chromatography, we found that isolated glycoprotein fractions (glps), but not flow-through fractions, exhibited strong growth inhibitory and bactericidal properties. The glps possessed two distinct functionalities: (a) specific binding and agglutination of bacterial cells, but not rat erythrocytes and (b) non-specific membrane permeabilization of both bacterial cells and erythrocytes. The isolated glps induced concentration- and time-dependent changes in the cell shape of both E. coli and B. subtilis as visualized by light and SEM microscopy. The appearance of filaments and spheroplasts correlated with growth inhibition and bactericidal effects, respectively. The time-kill kinetics showed a rapid, >5-log10 reduction of viable cells within 15 min incubation at 1xMBC, indicating that the glps-induced damage of the cell wall was lethal. Unexpectedly, MALDI-TOF and electrospray quadrupole time of flight mass spectrometry, (ESI-Q-TOF-MS/MS) analysis of glps showed sequence identity with the Major Royal Jelly Protein 1 (MRJP1) precursor that harbors three antimicrobial peptides: Jelleins 1, 2, and 4. The presence of high-mannose structures explained the lectin-like activity of MRJP1, while the presence of Jelleins in MRJP1 may explain cell wall disruptions. Thus, the observed damages induced by the MRJP1 to the bacterial cell wall constitute the mechanism by which the antibacterial effects were produced. Antibacterial activity of MRJP1 glps directly correlated with the overall antibacterial activity of honey, suggesting that it is honey's active principle responsible for this activity.

Infection dynamic of symbiotic bacteria in the pea aphid Acyrthosiphon pisum gut and host immune response at the early steps in the infection process

In addition to its obligatory symbiont Buchnera aphidicola, the pea aphid Acyrthosiphon pisum can harbor several facultative bacterial symbionts which can be mutualistic in the context of various ecological interactions. Belonging to a genus where many members have been described as pathogen in invertebrates, Serratia symbiotica is one of the most common facultative partners found in aphids. The recent discovery of strains able to grow outside their host allowed us to simulate environmental acquisition of symbiotic bacteria by aphids. Here, we performed an experiment to characterize the A. pisum response to the ingestion of the free-living S. symbiotica CWBI-2.3T in comparison to the ingestion of the pathogenic Serratia marcescens Db11 at the early steps in the infection process. We found that, while S. marcescens Db11 killed the aphids within a few days, S. symbiotica CWBI-2.3T did not affect host survival and colonized the whole digestive tract within a few days. Gene expression analysis of immune genes suggests that S. symbiotica CWBI-2.3T did not trigger an immune reaction, while S. marcescens Db11 did, and supports the hypothesis of a fine-tuning of the host immune response set-up for fighting pathogens while maintaining mutualistic partners. Our results also suggest that the lysosomal system and the JNK pathway are possibly involved in the regulation of invasive bacteria in aphids and that the activation of the JNK pathway is IMD-independent in the pea aphid.

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.

Evolutionary genetics of insect innate immunity

Patterns of evolution in immune defense genes help to understand the evolutionary dynamics between hosts and pathogens. Multiple insect genomes have been sequenced, with many of them having annotated immune genes, which paves the way for a comparative genomic analysis of insect immunity. In this review, I summarize the current state of comparative and evolutionary genomics of insect innate immune defense. The focus is on the conserved and divergent components of immunity with an emphasis on gene family evolution and evolution at the sequence level; both population genetics and molecular evolution frameworks are considered.

A genomic comparison of two termites with different social complexity

The termites evolved eusociality and complex societies before the ants, but have been studied much less. The recent publication of the first two termite genomes provides a unique comparative opportunity, particularly because the sequenced termites represent opposite ends of the social complexity spectrum. Zootermopsis nevadensis has simple colonies with totipotent workers that can develop into all castes (dispersing reproductives, nest-inheriting replacement reproductives, and soldiers). In contrast, the fungus-growing termite Macrotermes natalensis belongs to the higher termites and has very large and complex societies with morphologically distinct castes that are life-time sterile. Here we compare key characteristics of genomic architecture, focusing on genes involved in communication, immune defenses, mating biology and symbiosis that were likely important in termite social evolution. We discuss these in relation to what is known about these genes in the ants and outline hypothesis for further testing.

Parallel epigenomic and transcriptomic responses to viral infection in honey bees (Apis mellifera)

Populations of honey bees are declining throughout the world, with US beekeepers losing 30% of their colonies each winter. Though multiple factors are driving these colony losses, it is increasingly clear that viruses play a major role. However, information about the molecular mechanisms mediating antiviral immunity in honey bees is surprisingly limited. Here, we examined the transcriptional and epigenetic (DNA methylation) responses to viral infection in honey bee workers. One-day old worker honey bees were fed solutions containing Israeli Acute Paralysis Virus (IAPV), a virus which causes muscle paralysis and death and has previously been associated with colony loss. Uninfected control and infected, symptomatic bees were collected within 20-24 hours after infection. Worker fat bodies, the primary tissue involved in metabolism, detoxification and immune responses, were collected for analysis. We performed transcriptome- and bisulfite-sequencing of the worker fat bodies to identify genome-wide gene expression and DNA methylation patterns associated with viral infection. There were 753 differentially expressed genes (FDR<0.05) in infected versus control bees, including several genes involved in epigenetic and antiviral pathways. DNA methylation status of 156 genes (FDR<0.1) changed significantly as a result of the infection, including those involved in antiviral responses in humans. There was no significant overlap between the significantly differentially expressed and significantly differentially methylated genes, and indeed, the genomic characteristics of these sets of genes were quite distinct. Our results indicate that honey bees have two distinct molecular pathways, mediated by transcription and methylation, that modulate protein levels and/or function in response to viral infections.

Immunogene and viral transcript dynamics during parasitic Varroa destructor mite infection of developing honey bee (Apis mellifera) pupae

The ectoparasitic Varroa destructor mite is a major contributor to the ongoing honey bee health crisis. Varroa interacts with honey bee viruses, exacerbating their pathogenicity. In addition to vectoring viruses, immunosuppression of the developing honey bee hosts by Varroa has been proposed to explain the synergy between viruses and mites. However, the evidence for honey bee immune suppression by V. destructor is contentious. We systematically studied the quantitative effects of experimentally introduced V. destructor mites on immune gene expression at five specific time points during the development of the honey bee hosts. Mites reproduced normally and were associated with increased titers of deformed wing virus in the developing bees. Our data on different immune genes show little evidence for immunosuppression of honey bees by V. destructor. Experimental wounding of developing bees increases relative immune gene expression and deformed wing virus titers. Combined, these results suggest that mite feeding activity itself and not immunosuppression may contribute to the synergy between viruses and mites. However, our results also suggest that increased expression of honey bee immune genes decreases mite reproductive success, which may be explored to enhance mite control strategies. Finally, our expression data for multiple immune genes across developmental time and different experimental treatments indicates co-regulation of several of these genes and thus improves our understanding of the understudied honey bee immune system.

Causes and consequences of crossing-over evidenced via a high-resolution recombinational landscape of the honey bee

Background

Social hymenoptera, the honey bee (Apis mellifera) in particular, have ultra-high crossover rates and a large degree of intra-genomic variation in crossover rates. Aligned with haploid genomics of males, this makes them a potential model for examining the causes and consequences of crossing over. To address why social insects have such high crossing-over rates and the consequences of this, we constructed a high-resolution recombination atlas by sequencing 55 individuals from three colonies with an average marker density of 314 bp/marker.

Results

We find crossing over to be especially high in proximity to genes upregulated in worker brains, but see no evidence for a coupling with immune-related functioning. We detect only a low rate of non-crossover gene conversion, contrary to current evidence. This is in striking contrast to the ultrahigh crossing-over rate, almost double that previously estimated from lower resolution data. We robustly recover the predicted intragenomic correlations between crossing over and both population level diversity and GC content, which could be best explained as indirect and direct consequences of crossing over, respectively.

Conclusions

Our data are consistent with the view that diversification of worker behavior, but not immune function, is a driver of the high crossing-over rate in bees. While we see both high diversity and high GC content associated with high crossing-over rates, our estimate of the low non-crossover rate demonstrates that high non-crossover rates are not a necessary consequence of high recombination rates.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0566-0) 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.

Identification of immunity-related genes in the larvae of Protaetia brevitarsis seulensis (Coleoptera: Cetoniidae) by a next-generation sequencing-based transcriptome analysis

To identify immune-related genes in the larvae of white-spotted flower chafers, next-generation sequencing was conducted with an Illumina HiSeq2000, resulting in 100 million cDNA reads with sequence information from over 10 billion base pairs (bp) and >50× transcriptome coverage. A subset of 77,336 contigs was created, and ∼35,532 sequences matched entries against the NCBI nonredundant database (cutoff, e < 10(-5)). Statistical analysis was performed on the 35,532 contigs. For profiling of the immune response, samples were analyzed by aligning 42 base sequence tags to the de novo reference assembly, comparing levels in immunized larvae to control levels of expression. Of the differentially expressed genes, 3,440 transcripts were upregulated and 3,590 transcripts were downregulated. Many of these genes were confirmed as immune-related genes such as pattern recognition proteins, immune-related signal transduction proteins, antimicrobial peptides, and cellular response proteins, by comparison to published data.

Differential viral levels and immune gene expression in three stocks of Apis mellifera induced by different numbers of Varroa destructor

The viral levels and immune responses of Italian honey bees (IHB), Russian honey bees (RHB) and an outcross of Varroa Sensitive Hygienic bees (POL) deliberately infested with one or two foundress Varroa were compared. We found that the Deformed wing virus (DWV) level in IHB inoculated with one or two foundress Varroa increased to about 10(3) or 10(5) fold the levels of their uninfested brood. In contrast, POL (10(2) or 10(4) fold) and RHB (10(2) or l0(4) fold) supported a lower increase in DWV levels. The feeding of different stages of Varroa nymphs did not increase DWV levels of their pupal hosts. Analyses of their corresponding Varroa mites showed the same trends: two foundress Varroa yielded higher DWV levels than one foundress, and the addition of nymphs did not increase viral levels. Using the same pupae examined for the presence of viruses, 16 out of 24 genes evaluated showed significant differential mRNA expression levels among the three honey bee stocks. However, only four genes (Defensin, Dscam, PPOact and spaetzle), which were expressed at similar levels in uninfested pupae, were altered by the number of feeding foundress Varroa and levels of DWV regardless of stocks. This research provides the first evidence that immune response profiles of different honey bee stocks are induced by Varroa parasitism.

At the brink of eusociality: transcriptomic correlates of worker behaviour in a small carpenter bee

Background

There is great interest in understanding the genomic underpinnings of social evolution, in particular, the evolution of eusociality (caste-containing societies with non-reproductives that care for siblings). Subsociality is a key precursor for the evolution of eusociality and characterized by prolonged parental care and parent-offspring interaction. Here, we provide the first transcriptomic data for the small carpenter bee, Ceratina calcarata. This species is of special interest because it is subsocial and in the same family as the highly eusocial honey bee, Apis mellifera. In addition, some C. calcarata females demonstrate alloparental care without reproduction, which provides a unique opportunity to study worker behaviour in a non-eusocial species.

Results

We uncovered similar gene expression patterns related to maternal care and sibling care in different groups of females. This agrees with the maternal heterochrony hypothesis, specifically, that changes in timing of offspring care gene expression are related to worker behaviour in incipient insect societies. In addition, we also detected some similarity to caste-related gene expression patterns in highly eusocial honey bees, and uncovered large lifetime changes in gene expression that accompany shifts in reproductive and maternal care behaviour.

Conclusions

For Ceratina calcarata, we found that transcript expression profiles were most similar between sibling care and maternal care females. The maternal care behaviour exhibited post-reproductively by Ceratina mothers is concordant in terms of transcript expression with the alloparental care exhibited by workers. In line with theoretical predictions, our data are consistent with the maternal heterochrony hypothesis for the evolutionary development of worker behaviour in subsocial bees.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-014-0260-6) contains supplementary material, which is available to authorized users.

Genomic analysis of the interaction between pesticide exposure and nutrition in honey bees (Apis mellifera)

Populations of pollinators are in decline worldwide. These declines are best documented in honey bees and are due to a combination of stressors. In particular, pesticides have been linked to decreased longevity and performance in honey bees; however, the molecular and physiological pathways mediating sensitivity and resistance to pesticides are not well characterized. We explored the impact of coumaphos and fluvalinate, the two most abundant and frequently detected pesticides in the hive, on genome-wide gene expression patterns of honey bee workers. We found significant changes in 1118 transcripts, including genes involved in detoxification, behavioral maturation, immunity, and nutrition. Since behavioral maturation is regulated by juvenile hormone III (JH), we examined effects of these miticides on hormone titers; while JH titers were unaffected, titers of methyl farnesoate (MF), the precursor to JH, were decreased. We further explored the association between nutrition- and pesticide-regulated gene expression patterns and demonstrated that bees fed a pollen-based diet exhibit reduced sensitivity to a third pesticide, chlorpyrifos. Finally, we demonstrated that expression levels of several of the putative pesticide detoxification genes identified in our study and previous studies are also upregulated in response to pollen feeding, suggesting that these pesticides and components in pollen modulate similar molecular response pathways. Our results demonstrate that pesticide exposure can substantially impact expression of genes involved in several core physiological pathways in honey bee workers. Additionally, there is substantial overlap in responses to pesticides and pollen-containing diets at the transcriptional level, and subsequent analyses demonstrated that pollen-based diets reduce workers' pesticide sensitivity. Thus, providing honey bees and other pollinators with high quality nutrition may improve resistance to pesticides.

Large-scale coding sequence change underlies the evolution of postdevelopmental novelty in honey bees

Whether coding or regulatory sequence change is more important to the evolution of phenotypic novelty is one of biology's major unresolved questions. The field of evo-devo has shown that in early development changes to regulatory regions are the dominant mode of genetic change, but whether this extends to the evolution of novel phenotypes in the adult organism is unclear. Here, we conduct ten RNA-Seq experiments across both novel and conserved tissues in the honey bee to determine to what extent postdevelopmental novelty is based on changes to the coding regions of genes. We make several discoveries. First, we show that with respect to novel physiological functions in the adult animal, positively selected tissue-specific genes of high expression underlie novelty by conferring specialized cellular functions. Such genes are often, but not always taxonomically restricted genes (TRGs). We further show that positively selected genes, whether TRGs or conserved genes, are the least connected genes within gene expression networks. Overall, this work suggests that the evo-devo paradigm is limited, and that the evolution of novelty, postdevelopment, follows additional rules. Specifically, evo-devo stresses that high network connectedness (repeated use of the same gene in many contexts) constrains coding sequence change as it would lead to negative pleiotropic effects. Here, we show that in the adult animal, the converse is true: Genes with low network connectedness (TRGs and tissue-specific conserved genes) underlie novel phenotypes by rapidly changing coding sequence to perform new-specialized functions.

The i5k Workspace@NAL--enabling genomic data access, visualization and curation of arthropod genomes

The 5000 arthropod genomes initiative (i5k) has tasked itself with coordinating the sequencing of 5000 insect or related arthropod genomes. The resulting influx of data, mostly from small research groups or communities with little bioinformatics experience, will require visualization, dissemination and curation, preferably from a centralized platform. The National Agricultural Library (NAL) has implemented the i5k Workspace@NAL (http://i5k.nal.usda.gov/) to help meet the i5k initiative's genome hosting needs. Any i5k member is encouraged to contact the i5k Workspace with their genome project details. Once submitted, new content will be accessible via organism pages, genome browsers and BLAST search engines, which are implemented via the open-source Tripal framework, a web interface for the underlying Chado database schema. We also implement the Web Apollo software for groups that choose to curate gene models. New content will add to the existing body of 35 arthropod species, which include species relevant for many aspects of arthropod genomic research, including agriculture, invasion biology, systematics, ecology and evolution, and developmental research.

Immune responses of honeybees and their fitness costs as compared to bumblebees

Immune responses of invertebrates imply more than developing a merely unspecific response to an infection. Great interest has been raised to unveil whether this investment into immunity also involves fitness costs associated to the individual or the group. Focusing on the immune responses of honeybees, we use the well-studied insect bumblebee for comparison. Bumblebees are capable of producing specific immune responses to infections whereas this has not been assessed for honeybees so far. We investigated whether a prior bacterial encounter provides protection against a later exposure to the same or a different bacterium in honeybees. Additionally, we studied whether the foraging activities of honeybees and bumblebees are affected upon immune stimulation by assessing the flight performance. Finally, the acceptance behavior of nestmates toward immune-challenged honeybees was determined. Results show that despite stimulating the immune system of honeybees, no protective effects to infections were found. Further, honeybees were not affected by an immune challenge in their flight performance whereas bumblebees showed significant flight impairment. Immune-challenged honeybees showed lower survival rates than naive individuals when introduced into a regular colony. Here, we reveal different immune response-cost scenarios in honeybees and bumblebees for the first time.

Honey bee hemocyte profiling by flow cytometry

Multiple stress factors in honey bees are causing loss of bee colonies worldwide. Several infectious agents of bees are believed to contribute to this problem. The mechanisms of honey bee immunity are not completely understood, in part due to limited information about the types and abundances of hemocytes that help bees resist disease. Our study utilized flow cytometry and microscopy to examine populations of hemolymph particulates in honey bees. We found bee hemolymph includes permeabilized cells, plasmatocytes, and acellular objects that resemble microparticles, listed in order of increasing abundance. The permeabilized cells and plasmatocytes showed unexpected differences with respect to properties of the plasma membrane and labeling with annexin V. Both permeabilized cells and plasmatocytes failed to show measurable mitochondrial membrane potential by flow cytometry using the JC-1 probe. Our results suggest hemolymph particulate populations are dynamic, revealing significant differences when comparing individual hive members, and when comparing colonies exposed to diverse conditions. Shifts in hemocyte populations in bees likely represent changing conditions or metabolic differences of colony members. A better understanding of hemocyte profiles may provide insight into physiological responses of honey bees to stress factors, some of which may be related to colony failure.

Diversity of honey stores and their impact on pathogenic bacteria of the honeybee, Apis mellifera

Honeybee colonies offer an excellent environment for microbial pathogen development. The highest virulent, colony killing, bacterial agents are Paenibacillus larvae causing American foulbrood (AFB), and European foulbrood (EFB) associated bacteria. Besides the innate immune defense, honeybees evolved behavioral defenses to combat infections. Foraging of antimicrobial plant compounds plays a key role for this "social immunity" behavior. Secondary plant metabolites in floral nectar are known for their antimicrobial effects. Yet, these compounds are highly plant specific, and the effects on bee health will depend on the floral origin of the honey produced. As worker bees not only feed themselves, but also the larvae and other colony members, honey is a prime candidate acting as self-medication agent in honeybee colonies to prevent or decrease infections. Here, we test eight AFB and EFB bacterial strains and the growth inhibitory activity of three honey types. Using a high-throughput cell growth assay, we show that all honeys have high growth inhibitory activity and the two monofloral honeys appeared to be strain specific. The specificity of the monofloral honeys and the strong antimicrobial potential of the polyfloral honey suggest that the diversity of honeys in the honey stores of a colony may be highly adaptive for its "social immunity" against the highly diverse suite of pathogens encountered in nature. This ecological diversity may therefore operate similar to the well-known effects of host genetic variance in the arms race between host and parasite.

No genetic tradeoffs between hygienic behaviour and individual innate immunity in the honey bee, Apis mellifera

Many animals have individual and social mechanisms for combating pathogens. Animals may exhibit short-term physiological tradeoffs between social and individual immunity because the latter is often energetically costly. Genetic tradeoffs between these two traits can also occur if mutations that enhance social immunity diminish individual immunity, or vice versa. Physiological tradeoffs between individual and social immunity have been previously documented in insects, but there has been no study of genetic tradeoffs involving these traits. There is strong evidence that some genes influence both innate immunity and behaviour in social insects – a prerequisite for genetic tradeoffs. Quantifying genetic tradeoffs is critical for understanding the evolution of immunity in social insects and for devising effective strategies for breeding disease-resistant pollinator populations. We conducted two experiments to test the hypothesis of a genetic tradeoff between social and individual immunity in the honey bee, Apis mellifera. First, we estimated the relative contribution of genetics to individual variation in innate immunity of honey bee workers, as only heritable traits can experience genetic tradeoffs. Second, we examined if worker bees with hygienic sisters have reduced individual innate immune response. We genotyped several hundred workers from two colonies and found that patriline genotype does not significantly influence the antimicrobial activity of a worker’s hemolymph. Further, we did not find a negative correlation between hygienic behaviour and the average antimicrobial activity of a worker’s hemolymph across 30 honey bee colonies. Taken together, our work indicates no genetic tradeoffs between hygienic behaviour and innate immunity in honey bees. Our work suggests that using artificial selection to increase hygienic behaviour of honey bee colonies is not expected to concurrently compromise individual innate immunity of worker bees.

Microbial communities of three sympatric Australian stingless bee species

Bacterial symbionts of insects have received increasing attention due to their prominent role in nutrient acquisition and defense. In social bees, symbiotic bacteria can maintain colony homeostasis and fitness, and the loss or alteration of the bacterial community may be associated with the ongoing bee decline observed worldwide. However, analyses of microbiota associated with bees have been largely confined to the social honeybees (Apis mellifera) and bumblebees (Bombus spec.), revealing – among other taxa – host-specific lactic acid bacteria (LAB, genus Lactobacillus) that are not found in solitary bees. Here, we characterized the microbiota of three Australian stingless bee species (Apidae: Meliponini) of two phylogenetically distant genera (Tetragonula and Austroplebeia). Besides common plant bacteria, we find LAB in all three species, showing that LAB are shared by honeybees, bumblebees and stingless bees across geographical regions. However, while LAB of the honeybee-associated Firm4–5 clusters were present in Tetragonula, they were lacking in Austroplebeia. Instead, we found a novel clade of likely host-specific LAB in all three Australian stingless bee species which forms a sister clade to a large cluster of Halictidae-associated lactobacilli. Our findings indicate both a phylogenetic and geographical signal of host-specific LAB in stingless bees and highlight stingless bees as an interesting group to investigate the evolutionary history of the bee-LAB association.

Transcriptional response of Musca domestica larvae to bacterial infection

The house fly Musca domestica, a cosmopolitan dipteran insect, is a significant vector for human and animal bacterial pathogens, but little is known about its immune response to these pathogens. To address this issue, we inoculated the larvae with a mixture of Escherichia coli and Staphylococcus aureus and profiled the transcriptome 6, 24, and 48 h thereafter. Many genes known to controlling innate immunity in insects were induced following infection, including genes encoding pattern recognition proteins (PGRPs), various components of the Toll and IMD signaling pathways and of the proPO-activating and redox systems, and multiple antimicrobial peptides. Interestingly, we also uncovered a large set of novel immune response genes including two broad-spectrum antimicrobial peptides (muscin and domesticin), which might have evolved to adapt to house-fly's unique ecological environments. Finally, genes mediating oxidative phosphorylation were repressed at 48 h post-infection, suggesting disruption of energy homeostasis and mitochondrial function at the late stages of infection. Collectively, our data reveal dynamic changes in gene expression following bacterial infection in the house fly, paving the way for future in-depth analysis of M. domestica's immune system.

Parasitic and immune modulation of flight activity in honey bees tracked with optical counters

Host-parasite interactions are often characterized by changes in the host behaviour, which are beneficial to either the parasite or the host, or are a non-adaptive byproduct of parasitism. These interactions are further complicated in animal society because individual fitness is associated with group performance. However, a better understanding of host-parasite interaction in animal society first requires the identification of individual host behavioural modification. Therefore, we challenged honey bee (Apis mellifera) workers with the parasite Nosema ceranae or an immune stimulation and tracked their flight activity over their lifetime with an optic counter. We found that bees responded differently to each stress: both Nosema-infected and immune-challenged bees performed a lower number of daily flights compared with control bees, but the duration of their flights increased and decreased over time, respectively. Overall, parasitized bees spent more time in the field each day than control bees, and the inverse was true for immune-challenged bees. Despite the stress of immune challenge, bees had a survival similar to that of control bees likely because of their restricted activity. We discuss how those different behavioural modifications could be adaptive phenotypes. This study provides new insights into how biological stress can affect the behaviour of individuals living in society and how host responses have evolved.

A virulent strain of deformed wing virus (DWV) of honeybees (Apis mellifera) prevails after Varroa destructor-mediated, or in vitro, transmission

The globally distributed ectoparasite Varroa destructor is a vector for viral pathogens of the Western honeybee (Apis mellifera), in particular the Iflavirus Deformed Wing Virus (DWV). In the absence of Varroa low levels DWV occur, generally causing asymptomatic infections. Conversely, Varroa-infested colonies show markedly elevated virus levels, increased overwintering colony losses, with impairment of pupal development and symptomatic workers. To determine whether changes in the virus population were due Varroa amplifying and introducing virulent virus strains and/or suppressing the host immune responses, we exposed Varroa-naïve larvae to oral and Varroa-transmitted DWV. We monitored virus levels and diversity in developing pupae and associated Varroa, the resulting RNAi response and transcriptome changes in the host. Exposed pupae were stratified by Varroa association (presence/absence) and virus levels (low/high) into three groups. Varroa-free pupae all exhibited low levels of a highly diverse DWV population, with those exposed per os (group NV) exhibiting changes in the population composition. Varroa-associated pupae exhibited either low levels of a diverse DWV population (group VL) or high levels of a near-clonal virulent variant of DWV (group VH). These groups and unexposed controls (C) could be also discriminated by principal component analysis of the transcriptome changes observed, which included several genes involved in development and the immune response. All Varroa tested contained a diverse replicating DWV population implying the virulent variant present in group VH, and predominating in RNA-seq analysis of temporally and geographically separate Varroa-infested colonies, was selected upon transmission from Varroa, a conclusion supported by direct injection of pupae in vitro with mixed virus populations. Identification of a virulent variant of DWV, the role of Varroa in its transmission and the resulting host transcriptome changes furthers our understanding of this important viral pathogen of honeybees.

Honeybees are the most important managed pollinating insect, contributing billions of dollars to annual global agricultural production. Over the last century a parasitic mite, Varroa, has spread worldwide, with significant impacts on honeybee colony health as a consequence of its transmission of a cocktail of viruses while feeding on honeybee ‘blood’. The most important virus for colony health is deformed wing virus (DWV), high levels of which cause developmental deformities and premature ageing resulting in high overwintering colony losses. In experiments on individual Varroa-exposed pupae we demonstrate that a single type of virulent DWV is amplified 1,000–10,000 times in the recipient pupae, despite the mite containing a high diversity of replicating DWV strains. We could recapitulate this by direct injection of pupae with mixed virus populations, showing the virulent strain is advantaged by the route of transmission. In parallel, we detected changes in the immune response and developmental gene expression of the honeybee and propose that these contribute to the characteristic pathogenesis of DWV. Identification of a virulent strain of DWV has implications for therapeutic or prophylactic interventions to improve honeybee colony health, as well as contributing to our understanding of the biology of this important honeybee viral pathogen.

Molecular cloning and characterization of a lysozyme cDNA from the mole cricket Gryllotalpa orientalis (Orthoptera: Gryllotalpidae)

A full-length lysozyme cDNA from Gryllotalpa orientalis was cloned and sequenced. The deduced amino acid sequence of the lysozyme protein was 143 amino acids in length, with a calculated molecular mass of 15.84 kDa and an isoelectric point of 4.74. Sequence motifs, together with alignment and phylogenetic results, confirmed that G. orientalis lysozyme belongs to the C (chicken)-type lysozyme family of proteins. The protein sequence of lysozyme from G. orientalis showed high identity to that of Drosophila melanogaster (51.7 %); however, in contrast to D. melanogaster lysozyme, G. orientalis lysozyme was immune inducible and expressed in a wide range of tissues. Expression of G. orientalis lysozyme mRNA was highest at 8 h post-infection and subsequently decreased with time after bacterial infection. We also expressed G. orientalis lysozyme protein in vitro using the pET expression system. Compared with the negative control, over-expressed G. orientalis lysozyme showed antimicrobial activity against Gram-negative bacteria Escherichia coli and Gram-positive bacteria Bacillus subtilis by radial diffusion assay, with minimal inhibitory concentration values of 30.3 and 7.55 µM, respectively. These results indicate that G. orientalis lysozyme may have stronger antimicrobial activity than other lysozymes against a broad range of microorganisms.

Trans-generational immune priming in honeybees

Maternal immune experience acquired during pathogen exposure and passed on to progeny to enhance resistance to infection is called trans-generational immune priming (TgIP). In eusocial insects like honeybees, TgIP would result in a significant improvement of health at individual and colony level. Demonstrated in invertebrates other than honeybees, TgIP has not yet been fully elucidated in terms of intensity and molecular mechanisms underlying this response. Here, we immune-stimulated honeybee queens with Paenibacillus larvae (Pl), a spore-forming bacterium causing American Foulbrood, the most deadly bee brood disease worldwide. Subsequently, offspring of stimulated queens were exposed to spores of Pl and mortality rates were measured to evaluate maternal transfer of immunity. Our data substantiate the existence of TgIP effects in honeybees by direct evaluation of offspring resistance to bacterial infection. A further aspect of this study was to investigate a potential correlation between immune priming responses and prohaemocytes-haemocyte differentiation processes in larvae. The results point out that a priming effect triggers differentiation of prohaemocytes to haemocytes. However, the mechanisms underlying TgIP responses are still elusive and require future investigation.

Genomic analysis of the interactions between social environment and social communication systems in honey bees (Apis mellifera)

Social context is often a primary regulator of social behavior, but genes that affect or are affected by social context have rarely been investigated. In social insects, caste specific pheromones are key modulators of social behavior, e.g., in honey bees the queen mandibular gland (MG) pheromone mediates reproductive dominance, its absence prompting ovary activation and queen pheromone production in workers. Here, we investigate the effect of social environment on genome-wide expression patterns in the MG, to determine how social context modulates expression of genes that, in turn alter social environment. We used microarrays to examine the MGs of virgin and mated queens, and queenright (QR) and queenless (QL) workers with or without activated ovaries. Approximately 2554 transcripts were significantly differentially expressed among these groups, with caste and social context being the main regulators of gene expression patterns, while physiological state (ovary activation) only minimally affecting gene expression. Thus, social context strongly regulates expression of genes, which, in turn, shape social environment. Among these, 25 genes that are putatively involved in caste selective production of the fatty-acid derived MG pheromone were differentially expressed in queens and workers. These genes whose functions correspond with enzymatic or transport processes emphasize the occurrence of disparate pheromone biosynthetic pathways for queens and workers, adding another dimension regarding the regulation of these important pheromones. Gene ontology analysis also revealed genes of different functional categories whose expression was impacted by caste or by the social environment, suggesting that the MG has broader functions than pheromone biosynthesis.

Differential gene expression in laboratory strains of human head and body lice when challenged with Bartonella quintana, a pathogenic bacterium

Human head and body lice are obligatory hematophagous ectoparasites that belong to a single species, Pediculus humanus. Only body lice, however, are vectors of the infectious Gram-negative bacterium Bartonella quintana. Because of their near identical genomes, yet differential vector competence, head and body lice provide a unique model system to study the gain or loss of vector competence. Using our in vitro louse-rearing system, we infected head and body lice with blood containing B. quintana in order to detect both differences in the proliferation of B. quintana and transcriptional differences of immune-related genes in the lice. B. quintana proliferated rapidly in body lice at 6 days post-infection, but plateaued in head lice at 4 days post-infection. RNAseq and quantitative real-time PCR validation analyses determined gene expression differences. Eight immunoresponse genes were observed to be significantly different with many associated with the Toll pathway: Fibrinogen-like protein, Spaetzle, Defensin 1, Serpin, Scavenger receptor A and Apolipoporhrin 2. Our findings support the hypothesis that body lice, unlike head lice, fight infection from B. quintana only at the later stages of its proliferation.

Gene co-citation networks associated with worker sterility in honey bees

BACKGROUND

The evolution of reproductive self-sacrifice is well understood from kin theory, yet our understanding of how actual genes influence the expression of reproductive altruism is only beginning to take shape. As a model in the molecular study of social behaviour, the honey bee Apis mellifera has yielded hundreds of genes associated in their expression with differences in reproductive status of females, including genes directly associated with sterility, yet there has not been an attempt to link these candidates into functional networks that explain how workers regulate sterility in the presence of queen pheromone. In this study we use available microarray data and a co-citation analysis to describe what gene interactions might regulate a worker's response to ovary suppressing queen pheromone.

RESULTS

We reconstructed a total of nine gene networks that vary in size and gene composition, but that are significantly enriched for genes of reproductive function. The networks identify, for the first time, which candidate microarray genes are of functional importance, as evidenced by their degree of connectivity to other genes within each of the inferred networks. Our study identifies single genes of interest related to oogenesis, including eggless, and further implicates pathways related to insulin, ecdysteroid, and dopamine signaling as potentially important to reproductive decision making in honey bees.

CONCLUSIONS

The networks derived here appear to be variable in gene composition, hub gene identity, and the overall interactions they describe. One interpretation is that workers use different networks to control personal reproduction via ovary activation, perhaps as a function of age or environmental circumstance. Alternatively, the multiple networks inferred here may represent segments of the larger, single network that remains unknown in its entirety. The networks generated here are provisional but do offer a new multi-gene framework for understanding how honey bees regulate personal reproduction within their highly social breeding system.

How to kill the honey bee larva: genomic potential and virulence mechanisms of Paenibacillus larvae

Paenibacillus larvae, a Gram positive bacterial pathogen, causes American Foulbrood (AFB), which is the most serious infectious disease of honey bees. In order to investigate the genomic potential of P. larvae, two strains belonging to two different genotypes were sequenced and used for comparative genome analysis. The complete genome sequence of P. larvae strain DSM 25430 (genotype ERIC II) consisted of 4,056,006 bp and harbored 3,928 predicted protein-encoding genes. The draft genome sequence of P. larvae strain DSM 25719 (genotype ERIC I) comprised 4,579,589 bp and contained 4,868 protein-encoding genes. Both strains harbored a 9.7 kb plasmid and encoded a large number of virulence-associated proteins such as toxins and collagenases. In addition, genes encoding large multimodular enzymes producing nonribosomally peptides or polyketides were identified. In the genome of strain DSM 25719 seven toxin associated loci were identified and analyzed. Five of them encoded putatively functional toxins. The genome of strain DSM 25430 harbored several toxin loci that showed similarity to corresponding loci in the genome of strain DSM 25719, but were non-functional due to point mutations or disruption by transposases. Although both strains cause AFB, significant differences between the genomes were observed including genome size, number and composition of transposases, insertion elements, predicted phage regions, and strain-specific island-like regions. Transposases, integrases and recombinases are important drivers for genome plasticity. A total of 390 and 273 mobile elements were found in strain DSM 25430 and strain DSM 25719, respectively. Comparative genomics of both strains revealed acquisition of virulence factors by horizontal gene transfer and provided insights into evolution and pathogenicity.

Establishment of a bacterial infection model using the European honeybee, Apis mellifera L

Injection of human pathogenic bacteria (Pseudomonas aeruginosa, Serratia marcescens, Salmonella enterica, Staphylococcus aureus, and Listeria monocytogenes) into the hemocoel of honeybee (Apis mellifera L.) workers kills the infected bees. The bee-killing effects of the pathogens were affected by temperature, and the LD₅₀ values at 37°C were more than 100-fold lower than those at 15°C. Gene-disrupted S. aureus mutants of virulence genes such as agrA, saeS, arlR, srtA, hla, and hlb had attenuated bee-killing ability. Nurse bees were less susceptible than foragers and drones to S. aureus infection. Injection of antibiotics clinically used for humans had therapeutic effects against S. aureus infections of bees, and the ED₅₀ values of these antibiotics were comparable with those determined in mammalian models. Moreover, the effectiveness of orally administered antibiotics was consistent between honeybees and mammals. These findings suggest that the honeybee could be a useful model for assessing the pathogenesis of human-infecting bacteria and the effectiveness of antibiotics.

Gene expression differences underlying genotype-by-genotype specificity in a host-parasite system

In many systems, host-parasite evolutionary dynamics have led to the emergence and maintenance of diverse parasite and host genotypes within the same population. Genotypes vary in key attributes: Parasite genotypes vary in ability to infect, host genotypes vary in susceptibility, and infection outcome is frequently the result of both parties' genotypic identities. These host-parasite genotype-by-genotype (GH × GP) interactions influence evolutionary and ecological dynamics in important ways. Interactions can be produced through genetic variation; however, here, we assess the role of variable gene expression as an additional source of GH × GP interactions. The bumblebee Bombus terrestris and its trypanosome gut parasite Crithidia bombi are a model system for host-parasite matching. Full-transcriptome sequencing of the bumblebee host revealed that different parasite genotypes indeed induce fundamentally different host expression responses and host genotypes vary in their responses to the infecting parasite genotype. It appears that broadly and successfully infecting parasite genotypes lead to reduced host immune gene expression relative to unexposed bees but induce the expression of genes responsible for controlling gene expression. Contrastingly, a poorly infecting parasite genotype induced the expression of immunologically important genes, including antimicrobial peptides. A targeted expression assay confirmed the transcriptome results and also revealed strong host genotype effects. In all, the expression of a number of genes depends on the host genotype and the parasite genotype and the interaction between both host and parasite genotypes. These results suggest that alongside sequence variation in coding immunological genes, variation that controls immune gene expression can also produce patterns of host-parasite specificity.

Exposure to natural pathogens reveals costly aphid response to fungi but not bacteria

Immune responses are costly, causing trade-offs between defense and other host life history traits. Aphids present a special system to explore the costs associated with immune activation since they are missing several humoral and cellular mechanisms thought important for microbial resistance, and it is unknown whether they have alternative, novel immune responses to deal with microbial threat. Here we expose pea aphids to an array of heat-killed natural pathogens, which should stimulate immune responses without pathogen virulence, and measure changes in life-history traits. We find significant reduction in lifetime fecundity upon exposure to two fungal pathogens, but not to two bacterial pathogens. This finding complements recent genomic and immunological studies indicating that pea aphids are missing mechanisms important for bacterial resistance, which may have important implications for how aphids interact with their beneficial bacterial symbionts. In general, recent exploration of the immune systems of non-model invertebrates has called into question the generality of our current picture of insect immunity. Our data highlight that taking an ecological approach and measuring life-history traits to a broad array of pathogens provides valuable information that can complement traditional approaches.

Honey bees avoid nectar colonized by three bacterial species, but not by a yeast species, isolated from the bee gut

The gut microflora of the honey bee, Apis mellifera, is receiving increasing attention as a potential determinant of the bees' health and their efficacy as pollinators. Studies have focused primarily on the microbial taxa that appear numerically dominant in the bee gut, with the assumption that the dominant status suggests their potential importance to the bees' health. However, numerically minor taxa might also influence the bees' efficacy as pollinators, particularly if they are not only present in the gut, but also capable of growing in floral nectar and altering its chemical properties. Nonetheless, it is not well understood whether honey bees have any feeding preference for or against nectar colonized by specific microbial species. To test whether bees exhibit a preference, we conducted a series of field experiments at an apiary using synthetic nectar inoculated with specific species of bacteria or yeast that had been isolated from the bee gut, but are considered minor components of the gut microflora. These species had also been found in floral nectar. Our results indicated that honey bees avoided nectar colonized by the bacteria Asaia astilbes, Erwinia tasmaniensis, and Lactobacillus kunkeei, whereas the yeast Metschnikowia reukaufii did not affect the feeding preference of the insects. Our results also indicated that avoidance of bacteria-colonized nectar was caused not by the presence of the bacteria per se, but by the chemical changes to nectar made by the bacteria. These findings suggest that gut microbes may not only affect the bees' health as symbionts, but that some of the microbes may possibly affect the efficacy of A. mellifera as pollinators by altering nectar chemistry and influencing their foraging behavior.

Identification of immunity-related genes in Ostrinia furnacalis against entomopathogenic fungi by RNA-seq analysis

BACKGROUND

The Asian corn borer (Ostrinia furnacalis (Guenée)) is one of the most serious corn pests in Asia. Control of this pest with entomopathogenic fungus Beauveria bassiana has been proposed. However, the molecular mechanisms involved in the interactions between O. furnacalis and B. bassiana are unclear, especially under the conditions that the genomic information of O. furnacalis is currently unavailable. So we sequenced and characterized the transcriptome of O. furnacalis larvae infected by B. bassiana with special emphasis on immunity-related genes.

METHODOLOGY/PRINCIPAL FINDINGS

Illumina Hiseq2000 was used to sequence 4.64 and 4.72 Gb of the transcriptome from water-injected and B. bassiana-injected O. furnacalis larvae, respectively. De novo assembly generated 62,382 unigenes with mean length of 729 nt. All unigenes were searched against Nt, Nr, Swiss-Prot, COG, and KEGG databases for annotations using BLASTN or BLASTX algorithm with an E-value cut-off of 10(-5). A total of 35,700 (57.2%) unigenes were annotated to at least one database. Pairwise comparisons resulted in 13,890 differentially expressed genes, with 5,843 up-regulated and 8,047 down-regulated. Based on sequence similarity to homologs known to participate in immune responses, we totally identified 190 potential immunity-related unigenes. They encode 45 pattern recognition proteins, 33 modulation proteins involved in the prophenoloxidase activation cascade, 46 signal transduction molecules, and 66 immune responsive effectors, respectively. The obtained transcriptome contains putative orthologs for nearly all components of the Toll, Imd, and JAK/STAT pathways. We randomly selected 24 immunity-related unigenes and investigated their expression profiles using quantitative RT-PCR assay. The results revealed variant expression patterns in response to the infection of B. bassiana.

CONCLUSIONS/SIGNIFICANCE

This study provides the comprehensive sequence resource and expression profiles of the immunity-related genes of O. furnacalis. The obtained data gives an insight into better understanding the molecular mechanisms of innate immune processes in O. furnacalis larvae against B. bassiana.

Characterizing the infection-induced transcriptome of Nasonia vitripennis reveals a preponderance of taxonomically-restricted immune genes

The innate immune system in insects consists of a conserved core signaling network and rapidly diversifying effector and recognition components, often containing a high proportion of taxonomically-restricted genes. In the absence of functional annotation, genes encoding immune system proteins can thus be difficult to identify, as homology-based approaches generally cannot detect lineage-specific genes. Here, we use RNA-seq to compare the uninfected and infection-induced transcriptome in the parasitoid wasp Nasonia vitripennis to identify genes regulated by infection. We identify 183 genes significantly up-regulated by infection and 61 genes significantly down-regulated by infection. We also produce a new homology-based immune catalog in N. vitripennis, and show that most infection-induced genes cannot be assigned an immune function from homology alone, suggesting the potential for substantial novel immune components in less well-studied systems. Finally, we show that a high proportion of these novel induced genes are taxonomically restricted, highlighting the rapid evolution of immune gene content. The combination of functional annotation using RNA-seq and homology-based annotation provides a robust method to characterize the innate immune response across a wide variety of insects, and reveals significant novel features of the Nasonia immune response.

Honey bee hygienic behaviour does not incur a cost via removal of healthy brood

In the honey bee, hygienic behaviour, the removal of dead or diseased brood from capped cells by workers, is a heritable trait that confers colony-level resistance against brood diseases. This behaviour is quite rare. Only c. 10% of unselected colonies show high levels of hygiene. Previous studies suggested that hygiene might be rare because it also results in the removal of healthy brood, thereby imposing an ongoing cost even when brood diseases are absent. We tested this hypothesis by quantifying hygienic behaviour in 10 colonies using a standard technique, the freeze-killed brood (FKB) bioassay. At the same time, we also quantified the removal of untreated brood. The study colonies showed a wide range in hygienic behaviour, removing 19.7-100% of the FKB. The removal of untreated brood ranged from 2% to 44.4%. However, there was no correlation between the two removal rates for any of the four age groups of untreated brood studied (eggs, young larvae, older larvae from uncapped cells and larvae/pupae from capped cells). These results do not support the cost-to-healthy-brood hypothesis for the rarity of hygienic behaviour.

Effect of copper-induced oxidative stress on sclerotial differentiation and antioxidant properties of Penicillium thomii PT95 strain

Penicillium thomii PT95 strain was able to form abundant orange, sand-shaped sclerotia in which carotenoids were accumulated. The aim of this work was to determine the effects of copper-induced oxidative stress on the sclerotial differentiation and antioxidant properties of PT95 strain. The results showed that the time of exudates initiation, sclerotial initiation and sclerotial maturation of PT95 strain were advanced in 1-2 days under the copper-induced oxidative stress growth conditions. The analytical results of sclerotial biomass, carotenoids content in sclerotia showed that copper-induced oxidative stress favored the sclerotial differentiation and biosynthesis of carotenoids. Under the copper-induced oxidative stress growth conditions, the total phenolics content and DPPH free radical scavenging activity of sclerotia of this fungus were decreased as compared with the control. However, the oxidative stress induced by a lower amount of CuSO₄ in media could enhance significantly the reducing power of sclerotia.

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.

Neonicotinoid clothianidin adversely affects insect immunity and promotes replication of a viral pathogen in honey bees

Large-scale losses of honey bee colonies represent a poorly understood problem of global importance. Both biotic and abiotic factors are involved in this phenomenon that is often associated with high loads of parasites and pathogens. A stronger impact of pathogens in honey bees exposed to neonicotinoid insecticides has been reported, but the causal link between insecticide exposure and the possible immune alteration of honey bees remains elusive. Here, we demonstrate that the neonicotinoid insecticide clothianidin negatively modulates NF-κB immune signaling in insects and adversely affects honey bee antiviral defenses controlled by this transcription factor. We have identified in insects a negative modulator of NF-κB activation, which is a leucine-rich repeat protein. Exposure to clothianidin, by enhancing the transcription of the gene encoding this inhibitor, reduces immune defenses and promotes the replication of the deformed wing virus in honey bees bearing covert infections. This honey bee immunosuppression is similarly induced by a different neonicotinoid, imidacloprid, but not by the organophosphate chlorpyriphos, which does not affect NF-κB signaling. The occurrence at sublethal doses of this insecticide-induced viral proliferation suggests that the studied neonicotinoids might have a negative effect at the field level. Our experiments uncover a further level of regulation of the immune response in insects and set the stage for studies on neural modulation of immunity in animals. Furthermore, this study has implications for the conservation of bees, as it will contribute to the definition of more appropriate guidelines for testing chronic or sublethal effects of pesticides used in agriculture.