Nosema Tolerant Honeybees (Apis mellifera) Escape Parasitic Manipulation of Apoptosis

The untapped potential of cell culture in disentangling insect-microbial relationships

Cell culture is a powerful technique for the investigation of molecular mechanisms fundamental to health and disease in a diverse array of organisms. Cell lines offer several advantages, namely their simplistic approach and high degree of reproducibility. One field where cell culture has proven particularly useful is the study of the microbiome, where cell culture has led to the illumination of microbial influences on host immunity, nutrition, and physiology. Thus far, researchers have focused cell culture work predominantly on humans, but the growing field of insect microbiome research stands to benefit greatly from its application. Insects constitute one of Earth's most diverse and ancient life forms and, just as with humans, possess microbiomes with great significance to their health. Insects, which play critical roles in supporting food security and ecological stability, are facing increasing threats from agricultural intensification, climate change, and pesticide use. As the microbiome is closely tied to host health, gaining a more robust understanding is of increasing importance. In this review, we assert that the cultivation and utilization of insect gut cell lines in microbiome research will bridge critical knowledge gaps essential for informing insect management practices in a world under pressure.

A systematic review of honey bee (Apis mellifera, Linnaeus, 1758) infections and available treatment options

BACKGROUND

Honey bees and honeycomb bees are very valuable for wild flowering plants and economically important crops due to their role as pollinators. However, these insects confront many disease threats (viruses, parasites, bacteria and fungi) and large pesticide concentrations in the environment. Varroa destructor is the most prevalent disease that has had the most negative effects on the fitness and survival of different honey bees (Apis mellifera and A. cerana). Moreover, honey bees are social insects and this ectoparasite can be easily transmitted within and across bee colonies.

OBJECTIVE

This review aims to provide a survey of the diversity and distribution of important bee infections and possible management and treatment options, so that honey bee colony health can be maintained.

METHODS

We used PRISMA guidelines throughout article selection, published between January 1960 and December 2020. PubMed, Google Scholar, Scopus, Cochrane Library, Web of Science and Ovid databases were searched.

RESULTS

We have collected 132 articles and retained 106 articles for this study. The data obtained revealed that V. destructor and Nosema spp. were found to be the major pathogens of honey bees worldwide. The impact of these infections can result in the incapacity of forager bees to fly, disorientation, paralysis, and death of many individuals in the colony. We find that both hygienic and chemical pest management strategies must be implemented to prevent, reduce the parasite loads and transmission of pathogens. The use of an effective miticide (fluvalinate-tau, coumaphos and amitraz) now seems to be an essential and common practice required to minimise the impact of Varroa mites and other pathogens on bee colonies. New, alternative biofriendly control methods, are on the rise, and could be critical for maintaining honey bee hive health and improving honey productivity.

CONCLUSIONS

We suggest that critical health control methods be adopted globally and that an international monitoring system be implemented to determine honey bee colony safety, regularly identify parasite prevalence, as well as potential risk factors, so that the impact of pathogens on bee health can be recognised and quantified on a global scale.

Group Size Buffers against Energetic Stress in Honeybee Workers (Apis mellifera)

Individuals of the Western honeybee species, Apis mellifera, live in large groups of thousands of worker bees, a queen, and a few drones. Workers interact frequently with related individuals while performing various tasks. Although it is well understood why and under which conditions sociality can evolve, the consequences for individuals living in permanent groups are less well understood. As individuals of solitary species become stressed when kept at high density, it might be the opposite in obligate social species. Here, I use an experimental laboratory set-up to study the effect of varying group sizes on the magnitude and within-group variance of stress responses towards energetic and heat stress. While only a weak difference was found in the magnitude of an energetic stress response as a function of group size, the within-group variance showed a statistically significant positive relationship with group size for the glucose/trehalose ratio, a marker for energetic stress. The heat stress marker, hsp70AB gene expression, did not show any relationship to group size. Individuals of obligate social species seem to benefit from adaptations to permanent group living, e.g., buffering against stress, especially at a higher density of individuals. The consequences of infections and immune system activation in isolated individuals are discussed.

Transgenerational genomic analyses reveal allelic oscillation and purifying selection in a gut parasite Nosema ceranae

Standing genetic variation is the predominant source acted on by selection. Organisms with high genetic diversity generally show faster responses toward environmental change. Nosema ceranae is a microsporidian parasite of honey bees, infecting midgut epithelial cells. High genetic diversity has been found in this parasite, but the mechanism for the parasite to maintain this diversity remains unclear. This study involved continuous inoculation of N. ceranae to honey bees. We found that the parasites slowly increased genetic diversity over three continuous inoculations. The number of lost single nucleotide variants (SNVs) was balanced with novel SNVs, which were mainly embedded in coding regions. Classic allele frequency oscillation was found at the regional level along the genome, and the associated genes were enriched in apoptosis regulation and ATP binding. The ratio of synonymous and non-synonymous substitution suggests a purifying selection, and our results provide novel insights into the evolutionary dynamics in microsporidian parasites.

Spillover and genome selection of the gut parasite Nosema ceranae between honey bee species

Nosema ceranae is a honey bee gut parasite that has recently spilled to another honey bee host through trading. The impact of infection on the native host is minor, which is substantial in the novel host. In this study, artificial inoculation simulated the parasite transmission from the native to the novel host. We found that the parasite initiated proliferation earlier in the novel host than in the native host. Additionally, parasite gene expression was significantly higher when infecting the novel host compared with the native host, leading to a significantly higher number of spores. Allele frequencies were similar for spores of parasites infecting both native and novel hosts. This suggests that the high number of spores found in the novel host was not caused by a subset of more fit spores from native hosts. Native hosts also showed a higher number of up-regulated genes in response to infection when compared with novel hosts. Our data further showed that native hosts suppressed parasite gene expression and arguably sacrificed cells to limit the parasite. The results provide novel insights into host defenses and gene selection during a parasite spillover event.

Diet Supplementation Helps Honey Bee Colonies in Combat Infections by Enhancing their Hygienic Behaviour

The hygienic behavior in honey bees is a complex polygenic trait that serves as a natural defense mechanism against bacterial and fungal brood diseases and Varroa destructor mites infesting brood cells. The aim of this study was to evaluate the effect of a dietary amino acids and vitamins supplement “BEEWELL AminoPlus” on hygienic behavior of Apis mellifera colonies combating microsporidial and viral infections. The experiment was performed during a one-year period on 40 colonies alloted to five groups: one supplemented and infected with Nosema ceranae and four viruses (Deformed wing virus - DWV, Acute bee paralysis virus - ABPV, Chronic bee paralysis virus - CBPV and Sacbrood virus – SBV), three not supplemented, but infected with N. ceranae and/ or viruses, and one negative control group. Beside the l isted pathogens, honey bee trypanosomatids were also monitored in all groups.

The supplement “BEEWELL AminoPlus” induced a significant and consistent increase of the hygienic behavior in spite of the negative effects of N. ceranae and viral infections. N. ceranae and viruses significantly and consistently decreased hygienic behavior, but also threatened the survival of bee colonies. The tested supplement showed anti-Nosema effect, since the N. ceranae infection level significantly and consistently declined only in the supplemented group. Among infected groups, only the supplemented one remained Lotmaria passim-free throughout the study. In conclusion, diet supplementation enhances hygienic behavior of honey bee colonies and helps them fight the most common infections of honey bees.

Recent Advances in the Biocontrol of Nosemosis in Honey Bees (Apis mellifera L.)

Nosemosis is a disease triggered by the single-celled spore-forming fungi Nosema apis and Nosema ceranae, which can cause extensive colony losses in honey bees (Apis mellifera L.). Fumagillin is an effective antibiotic treatment to control nosemosis, but due to its toxicity, it is currently banned in many countries. Accordingly, in the beekeeping sector, there is a strong demand for alternative ecological methods that can be used for the prevention and therapeutic control of nosemosis in honey bee colonies. Numerous studies have shown that plant extracts, RNA interference (RNAi) and beneficial microbes could provide viable non-antibiotic alternatives. In this article, recent scientific advances in the biocontrol of nosemosis are summarized.

Comparative Transcriptome Investigation of Nosema ceranae Infecting Eastern Honey Bee Workers

Simple Summary

At present, interaction between Nosema ceranae and Apis cerana is poorly understood, though A. cerana is the original host for N. ceranae. Here, comparative investigation was conducted using transcriptome data from N. ceranae infecting Apis cerana cerana workers at seven days post inoculation (dpi) and 10 dpi (NcT1 and NcT2 groups) as well as N. ceranae spores (NcCK group). There were 1411, 604, and 38 DEGs identified in NcCK vs. NcT1, NcCK vs. NcT2, and NcT1 vs. NcT2 comparison groups. Additionally, 10 upregulated genes and nine downregulated ones were shared by above-mentioned comparison groups. GO classification and KEGG pathway analysis suggested that these DEGs were engaged in a number of key functional terms and pathways such as cell part and glycolysis. Further analysis indicated that most of virulence factor-encoding genes were upregulated, while a few were downregulated during the fungal infection. Findings in this current work provide a basis for clarifying the molecular mechanism udnerlying N. ceranae infection and host-microsporidian interaction during bee nosemosis.

Abstract

Apis cerana is the original host for Nosema ceranae, a widespread fungal parasite resulting in honey bee nosemosis, which leads to severe losses to the apiculture industry throughout the world. However, knowledge of N. ceranae infecting eastern honey bees is extremely limited. Currently, the mechanism underlying N. ceranae infection is still largely unknown. Based on our previously gained high-quality transcriptome datasets derived from N. ceranae spores (NcCK group), N. ceranae infecting Apis cerana cerana workers at seven days post inoculation (dpi) and 10 dpi (NcT1 and NcT2 groups), comparative transcriptomic investigation was conducted in this work, with a focus on virulence factor-associated differentially expressed genes (DEGs). Microscopic observation showed that the midguts of A. c. cerana workers were effectively infected after inoculation with clean spores of N. ceranae. In total, 1411, 604, and 38 DEGs were identified from NcCK vs. NcT1, NcCK vs. NcT2, and NcT1 vs. NcT2 comparison groups. Venn analysis showed that 10 upregulated genes and nine downregulated ones were shared by the aforementioned comparison groups. The GO category indicated that these DEGs were involved in a series of functional terms relevant to biological process, cellular component, and molecular function such as metabolic process, cell part, and catalytic activity. Additionally, KEGG pathway analysis suggested that the DEGs were engaged in an array of pathways of great importance such as metabolic pathway, glycolysis, and the biosynthesis of secondary metabolites. Furthermore, expression clustering analysis demonstrated that the majority of genes encoding virulence factors such as ricin B lectins and polar tube proteins displayed apparent upregulation, whereas a few virulence factor-associated genes such as hexokinase gene and 6-phosphofructokinase gene presented downregulation during the fungal infection. Finally, the expression trend of 14 DEGs was confirmed by RT-qPCR, validating the reliability of our transcriptome datasets. These results together demonstrated that an overall alteration of the transcriptome of N. ceranae occurred during the infection of A. c. cerana workers, and most of the virulence factor-related genes were induced to activation to promote the fungal invasion. Our findings not only lay a foundation for clarifying the molecular mechanism underlying N. ceranae infection of eastern honey bee workers and microsporidian–host interaction.

Nosemosis Prevention and Control

Nosemosis is a serious microsporidian disease of adult European honey bees caused by the spore-forming unicellular fungi Nosema apis and Nosema ceranae. In this paper we describe the currently known techniques for nosemosis prevention and control including Good Beekeeping Practices (GBPs) and biosecurity measures (BMBs). Topics such as queen renewal, nosema-resistant bees and hygienic and control methods are described. Strong efforts are currently provided to find more a sustainable solution than the use of antibiotics. So far, it seems that the best way to approach nosemosis is given by an “integrated pest management strategy”, which foresees the contemporary application of different, specific GBPs and BMBs.

Nosema apis and N. ceranae Infection in Honey bees: A Model for Host-Pathogen Interactions in Insects

There has been increased focus on the role of microbial attack as a potential cause of recent declines in the health of the western honey bee, Apis mellifera. The Nosema species, N. apis and N. ceranae, are microsporidian parasites that are pathogenic to honey bees, and infection by these species has been implicated as a key factor in honey bee losses. Honey bees infected with both Nosema spp. display significant changes in their biology at the cellular, tissue, and organismal levels impacting host metabolism, immune function, physiology, and behavior. Infected individuals lead to colony dysfunction and can contribute to colony disease in some circumstances. The means through which parasite growth and tissue pathology in the midgut lead to the dramatic physiological and behavioral changes at the organismal level are only partially understood. In addition, we possess only a limited appreciation of the elements of the host environment that impact pathogen growth and development. Critical for answering these questions is a mechanistic understanding of the host and pathogen machinery responsible for host-pathogen interactions. A number of approaches are already being used to elucidate these mechanisms, and promising new tools may allow for gain- and loss-of-function experiments to accelerate future progress.

Reproductive Potential Impacts Body Maintenance Parameters and Global DNA Methylation in Honeybee Workers (Apis mellifera L.)

Simple Summary

The queens and sterile workers arise from genetically identical eggs but as imagoes, they differ in their life span, DNA methylation, and their functions. In the absence of the queen, the larvae develop into rebels, i.e., workers with increased reproductive potential. We assumed that since rebels are similar to the queen in many anatomical and behavioral features, they live longer and have lower levels of global DNA methylation, even when infected, e.g., by Nosema spp. Rebels always lived longer in comparison in normal workers and unexpectedly extended longevity of normal workers when they were together, similarly as the presence of a queen did. Rebels became infected more easily but tolerated the infection better. They also had lower level of global DNA methylation than normal workers. These features expand possibilities of the use of honeybees as a model for studies on senescence, nosemosis, eusocial evolution, and epigenetics.

Abstract

The widely accepted hypothesis in life history evolution about the trade-off between fecundity and longevity is not confirmed by long-living and highly fecund queens in eusocial insects. The fact that the queens and facultatively sterile workers usually arise from genetically identical eggs but differ in DNA methylation makes them a good model for studies on senescence, eusocial evolution, and epigenetics. Therefore, honeybees seem to be especially useful here because of long living rebel-workers (RW) with high reproductive potential recently described. Longevity, ovariole number, nosema tolerance, and global DNA methylation have been assayed in normal workers (NW) versus RW in hives and cages. RW always lived longer than NW and unexpectedly extended longevity of NW when they were together, similarly as the presence of a queen did. RW lived longer despite the fact that they had higher Nosema spore load; surprisingly they became infected more easily but tolerated the infection better. Global DNA methylation increased with age, being lower in RW than in NW. Therefore, RW are queen-like considering global DNA methylation and the link between fecundity, longevity, and body maintenance. Presented features of RW expands possibilities of the use of honeybees as a model for studies on senescence, nosemosis, eusocial evolution, and epigenetics.

Fungi of entomopathogenic potential in Chytridiomycota and Blastocladiomycota, and in fungal allies of the Oomycota and Microsporidia

The relationship between entomopathogenic fungi and their insect hosts is a classic example of the co-evolutionary arms race between pathogen and target host. The present review describes the entomopathogenic potential of Chytridiomycota and Blastocladiomycota fungi, and two groups of fungal allies: Oomycota and Microsporidia. The Oomycota (water moulds) are considered as a model biological control agent of mosquito larvae. Due to their shared ecological and morphological similarities, they had long been considered a part of the fungal kingdom; however, phylogenetic studies have since placed this group within the Straminipila. The Microsporidia are parasites of economically-important insects, including grasshoppers, lady beetles, bumblebees, colorado potato beetles and honeybees. They have been found to display some fungal characteristics, and phylogenetic studies suggest that they are related to fungi, either as a basal branch or sister group. The Blastocladiomycota and Chytridiomycota, named the lower fungi, historically were described together; however, molecular phylogenetic and ultrastructural research has classified them in their own phylum. They are considered parasites of ants, and of the larval stages of black flies, mosquitoes and scale insects.

Mitigating Nosema ceranae infection in western honey bee (Apis mellifera) workers using propolis collected from honey bee and stingless bee (Tetrigona apicalis) hives

Beekeepers need sustainable control options to treat Nosema ceranae infection in colonies of western honey bees (Apis mellifera L.) they manage. Propolis is a natural product derived from plant resins and contains chemical compounds with potential antimicrobial activity against N. ceranae. Here, we determined the efficacy of propolis from A. mellifera (USA) and Tetrigona apicalis (stingless bees, Thailand) colonies as treatments for N. ceranae infection in honey bee workers. Newly emerged bees were individually fed 2 μL of 50% (w/v) sucrose solution containing 1 × 10⁵N. ceranae spores. Following this, the infected bees were treated with 50% propolis extracted from A. mellifera or T. apicalis hives and fed in 50% sucrose solution (v/v). All bees were maintained at 34 ± 2 °C and 55 ± 5% RH. Dead bees were counted daily for 30 d to calculate survival. We also determined infection rate (# infected bees/100 bees), infectivity (number of spores per bee) and protein content in the hypopharyngeal glands and hemolymph on 7, 14, and 21 d post infection as measures of bee health. Propolis from both bee species significantly reduced bee mortality, infection rate and infectivity compared with those of untreated bees and led to significantly greater protein contents in hypopharyngeal glands and hemolymph in treated bees than in untreated ones (p < 0.0001). In conclusion, propolis from A. mellifera and T. apicalis colonies shows promise as a control against N. ceranae infection in honey bees.

Immune Response of Eastern Honeybee Worker to Nosema ceranae Infection Revealed by Transcriptomic Investigation

Simple Summary

Currently, knowledge regarding Apis cerana–Nosema ceranae interaction is very limited, though A. cerana is the original host of N. ceranae. Apis cerana cerana is a subspecies of A. cerana and a major bee species used in the beekeeping industry in China and other countries. Here, the effective infection of A. c. cerana workers by N. ceranae was verified, followed by transcriptomic investigation of host responses. Furthermore, immune responses between A. c. cerana and Apis mellifera ligustica were deeply compared and discussed. In total, 1127 and 957 N. ceranae-responsive genes were identified in the infected midguts at 7 d post-inoculation (dpi) and 10 dpi, respectively. Additionally, DEGs in workers’ midguts at both 7 dpi and 10 dpi were associated with six cellular immune pathways and three humoral immune pathways. Noticeably, one up-regulated gene was enriched in the NF-κB signaling pathway in the midgut at 10 dpi. Further analysis indicated that different cellular and humoral immune responses were employed by A. c. cerana and A. m. ligustica workers to combat N. ceranae. Our findings provide a foundation for clarifying the mechanisms regulating the immune response of A. c. cerana workers to N. ceranae invasion and developing new approaches to control bee microsporidiosis.

Abstract

Here, a comparative transcriptome investigation was conducted based on high-quality deep sequencing data from the midguts of Apis cerana cerana workers at 7 d post-inoculation (dpi) and 10 dpi with Nosema ceranae and corresponding un-inoculated midguts. PCR identification and microscopic observation of paraffin sections confirmed the effective infection of A. c. cerana worker by N. ceranae. In total, 1127 and 957 N. ceranae-responsive genes were identified in the infected midguts at 7 dpi and 10 dpi, respectively. RT-qPCR results validated the reliability of our transcriptome data. GO categorization indicated the differentially expressed genes (DEGs) were respectively engaged in 34 and 33 functional terms associated with biological processes, cellular components, and molecular functions. Additionally, KEGG pathway enrichment analysis showed that DEGs at 7 dpi and 10 dpi could be enriched in 231 and 226 pathways, respectively. Moreover, DEGs in workers’ midguts at both 7 dpi and 10 dpi were involved in six cellular immune pathways such as autophagy and phagosome and three humoral immune pathways such as the Toll/Imd signaling pathway and Jak-STAT signaling pathway. In addition, one up-regulated gene (XM_017055397.1) was enriched in the NF-κB signaling pathway in the workers’ midgut at 10 dpi. Further investigation suggested the majority of these DEGs were engaged in only one immune pathway, while a small number of DEGs were simultaneously involved in two immune pathways. These results together demonstrated that the overall gene expression profile in host midgut was altered by N. ceranae infection and some of the host immune pathways were induced to activation during fungal infection, whereas some others were suppressed via host–pathogen interaction. Our findings offer a basis for clarification of the mechanism underlying the immune response of A. c. cerana workers to N. ceranae infection, but also provide novel insights into eastern honeybee-microsporodian interaction.

A SNP assay for assessing diversity in immune genes in the honey bee (Apis mellifera L.)

With a growing number of parasites and pathogens experiencing large-scale range expansions, monitoring diversity in immune genes of host populations has never been so important because it can inform on the adaptive potential to resist the invaders. Population surveys of immune genes are becoming common in many organisms, yet they are missing in the honey bee (Apis mellifera L.), a key managed pollinator species that has been severely affected by biological invasions. To fill the gap, here we identified single nucleotide polymorphisms (SNPs) in a wide range of honey bee immune genes and developed a medium-density assay targeting a subset of these genes. Using a discovery panel of 123 whole-genomes, representing seven A. mellifera subspecies and three evolutionary lineages, 180 immune genes were scanned for SNPs in exons, introns (< 4 bp from exons), 3' and 5´UTR, and < 1 kb upstream of the transcription start site. After application of multiple filtering criteria and validation, the final medium-density assay combines 91 quality-proved functional SNPs marking 89 innate immune genes and these can be readily typed using the high-sample-throughput iPLEX MassARRAY system. This medium-density-SNP assay was applied to 156 samples from four countries and the admixture analysis clustered the samples according to their lineage and subspecies, suggesting that honey bee ancestry can be delineated from functional variation. In addition to allowing analysis of immunogenetic variation, this newly-developed SNP assay can be used for inferring genetic structure and admixture in the honey bee.

Potential of Fumagillin and Agaricus blazei Mushroom Extract to Reduce Nosema ceranae in Honey Bees

Depending on the infection level and colony strength, Nosema ceranae, a microsporidian endoparasite of the honey bee may have significant consequences on the health, reproduction and productivity of bee colonies. Despite exerting some side effects, fumagillin is most often used for Nosema control. In this study, in a cage experiment, N. ceranae infected bees were treated with fumagillin or the extract of Agaricus blazei mushroom, a possible alternative for Nosema control. Bee survival, Nosema spore loads, the expression levels of immune-related genes and parameters of oxidative stress were observed. Fumagillin treatment showed a negative effect on monitored parameters when applied preventively to non-infected bees, while a noticeable anti-Nosema effect and protection from Nosema-induced immunosuppression and oxidative stress were proven in Nosema-infected bees. However, a protective effect of the natural A. blazei extract was detected, without any side effects but with immunostimulatory activity in the preventive application. The results of this research suggest the potential of A. blazei extract for Nosema control, which needs to be further investigated.

A combined LC-MS and NMR approach to reveal metabolic changes in the hemolymph of honeybees infected by the gut parasite Nosema ceranae

Nosema ceranae is an emerging and invasive gut pathogen in Apis mellifera and is considered as a factor contributing to the decline of honeybee populations. Here, we used a combined LC-MS and NMR approach to reveal the metabolomics changes in the hemolymph of honeybees infected by this obligate intracellular parasite. For metabolic profiling, hemolymph samples were collected from both uninfected and N. ceranae-infected bees at two time points, 2 days and 10 days after the experimental infection of emergent bees. Hemolymph samples were individually analyzed by LC-MS, whereas each NMR spectrum was obtained from a pool of three hemolymphs. Multivariate statistical PLS-DA models clearly showed that the age of bees was the parameter with the strongest effect on the metabolite profiles. Interestingly, a total of 15 biomarkers were accurately identified and were assigned as candidate biomarkers representative of infection alone or combined effect of age and infection. These biomarkers included carbohydrates (α/β glucose, α/β fructose and hexosamine), amino acids (histidine and proline), dipeptides (Glu-Thr, Cys-Cys and γ-Glu-Leu/Ile), metabolites involved in lipid metabolism (choline, glycerophosphocholine and O-phosphorylethanolamine) and a polyamine compound (spermidine). Our study demonstrated that this untargeted metabolomics-based approach may be useful for a better understanding of pathophysiological mechanisms of the honeybee infection by N. ceranae.

Innate and Adaptive Immune Responses Against Microsporidia Infection in Mammals

Microsporidia are obligate intracellular and eukaryotic pathogens that can infect immunocompromised and immunocompetent mammals, including humans. Both innate and adaptive immune systems play important roles against microsporidian infection. The innate immune system can partially eliminate the infection by immune cells, such as gamma delta T cell, natural killer cells (NKs), macrophages and dendritic cells (DCs), and present the pathogens to lymphocytes. The innate immune cells can also prime and enhance the adaptive immune response via surface molecules and secreted cytokines. The adaptive immune system is critical to eliminate microsporidian infection by activating cytotoxic T lymphocyte (CTL) and humoral immune responses, and feedback regulation of the innate immune mechanism. In this review, we will discuss the cellular and molecular responses and functions of innate and adaptive immune systems against microsporidian infection.

Brood parasitism in eusocial insects (Hymenoptera): role of host geographical range size and phylogeny

Interspecific brood parasitism is common in many animal systems. Brood parasites enter the nests of other species and divert host resources for producing their own offspring, which can lead to strong antagonistic parasite-host coevolution. Here, we look at commonalities among social insect species that are victims of brood parasites, and use phylogenetic data and information on geographical range size to predict which species are most probably to fall victims to brood parasites in the future. In our analyses, we focus on three eusocial hymenopteran groups and their brood parasites: (i) bumblebees, (ii) Myrmica ants, and (iii) vespine and polistine wasps. In these groups, some, but not all, species are parasitized by obligate workerless inquilines that only produce reproductive-caste descendants. We find phylogenetic signals for geographical range size and the presence of parasites in bumblebees, but not in ants and wasps. Phylogenetic logistic regressions indicate that the probability of being attacked by one or more brood parasite species increases with the size of the geographical range in bumblebees, but the effect is statistically only marginally significant in ants. However, non-phylogenetic logistic regressions suggest that bumblebee species with the largest geographical range sizes may have a lower likelihood of harbouring social parasites than do hosts with medium-sized ranges. Our results provide new insights into the ecology and evolution of host-social parasite systems, and indicate that host phylogeny and geographical range size can be used to predict threats posed by social parasites, as well to design efficient conservation measures for both hosts and their parasites. This article is part of the theme issue 'The coevolutionary biology of brood parasitism: from mechanism to pattern'.

A Pediococcus strain to rescue honeybees by decreasing Nosema ceranae- and pesticide-induced adverse effects

Honeybees ensure a key ecosystemic service by pollinating many agricultural crops and wild plants. However, since few decades, managed bee colonies have declined worldwide. This phenomenon is considered to be multifactorial, with a strong emphasis on both parasites and pesticides. Infection by the parasite Nosema ceranae and exposure to pesticides can contribute to adverse effects, resulting in a perturbation of the honeybee physiology. We thus hypothesized that probiotic treatment could be promising to treat or prevent these disturbances. The aim of this study was to evaluate the effects of probiotics on N. ceranae-infected and intoxicated honeybees (by the insecticide thiamethoxam and the fungicide boscalid). For this purpose, experiments were conducted with five probiotics. Among them, Pediococcus acidilactici (PA) showed the best protective effect against the parasite and pesticides. PA significantly improved the infected honeybee lifespan as prophylactic and curative treatments (respectively 2.3 fold and 1.7 fold). Furthermore, the exposure to pesticides induced an increase of honeybee mortality compared with the control group (p < .001) that was restored by the PA treatment. Despite its beneficial effect on honeybee lifespan, the PA administration did not induce changes in the gut bacterial communities (neither in abundance or diversity). N. ceranae and the pesticides were shown to deregulate genes involved in honeybee development (vitellogenin), immunity (serine protease 40, defensin) and detoxification system (glutathione peroxidase-like 2, catalase), and these effects were corrected by the PA treatment. This study highlights the promising use of PA to protect honeybees from both pathogens and pesticides.

Resistance of native honey bees from Rhodope Mountains and lowland regions of Bulgaria to Nosema ceranae and viral pathogens

The Western honey bee (Apis mellifera L., Hymenoptera: Apidae) is a species of fundamental economic, agricultural and environmental importance. The aim of this study was to compare the prevalence of some parasitic and viral pathogens in local honey bees from the Rodope Mountains and plain regions. To achieve this goal, molecular screening for two of the most distributed Nosema spp. and molecular identification of six honey bee viruses – Deformed wing virus (DWV), Acute bee paralysis virus (ABPV), Chronic bee paralysis virus (CBPV), Sacbrood virus (SBV), Kashmir bee virus (KBV), and Black queen cell virus (BQCV) was performed. Molecular analysis was carried out on 168 honey bee samples from apiaries situated in three different parts of the country where a mix of different honey bee subspecies were reared. In South Bulgaria (the Rhodope Mountains), a local honey bee called Apis mellifera rodopica (a local ecotype of A. m. macedonica) was bred, while in the other two regions (plains) different introduced subspecies existed. The results showed that the samples from the lowland regions in the country were outlined with the highest prevalence (70.5%) of N. ceranae, while those from the mountainous parts had the lowest rate (5.2%). Four of the honey bee viruses were identified – DWV (10/5.9%), followed by SBV (6/3.6%) and ABPV (2/1.2%), and one case of BQCV. In conclusion, the local honey bee A. m. rodopica (despite the higher number of samples) has shown lower prevalence of both nosemosis and viral infections. Therefore, this honey bee has to be preserved as a part of the national biodiversity.

Age and Method of Inoculation Influence the Infection of Worker Honey Bees (Apis mellifera) by Nosema ceranae

The microsporidian parasite Nosema ceranae is a highly prevalent, global honey bee pathogen. Apis mellifera is considered to be a relatively recent host for this microsporidia, which raises questions as to how it affects its host's physiology, behavior and longevity, both at the individual and colony level. As such, honey bees were inoculated with fresh purified spores of this pathogen, both individually (Group A) or collectively (Group B) and they were studied from 0 to 15 days post-emergence (p.e.) to evaluate the effect of bee age and the method of inoculation at 7 days post-infection. The level of infection was analyzed individually by qPCR by measuring the relative amount of the N. ceranae polar tubule protein 3 (PTP3) gene. The results show that the bee's age and the method of infection directly influence parasite load, and thus, early disease development. Significant differences were found regarding bee age at the time of infection, whereby the youngest bees (new-born and 1 day p.e.) developed the highest parasite load, with this load decreasing dramatically in bees infected at 2 days p.e. before increasing again in bees infected at 3-4 days p.e. The parasite load in bees infected when older than 4 days p.e. diminished as they aged. When the age cohort data was pooled and grouped according to the method of infection, a significantly higher mean concentration and lower variation in N. ceranae infection was evident in Group A, indicating greater variation in experimental infection when spores were administered collectively to bees through their food. In summary, these data indicate that both biological and experimental factors should be taken into consideration when comparing data published in the literature.

Comparative Identification of MicroRNAs in Apis cerana cerana Workers' Midguts in Responseto Nosema ceranae Invasion

Here, the expression profiles and differentially expressed miRNAs (DEmiRNAs) in the midguts of Apis cerana cerana workers at 7 d and 10 d post-inoculation (dpi) with N. ceranae were investigated via small RNA sequencing and bioinformatics. Five hundred and twenty nine (529) known miRNAs and 25 novel miRNAs were identified in this study, and the expression of 16 predicted miRNAs was confirmed by Stem-loop RT-PCR. A total of 14 DEmiRNAs were detected in the midgut at 7 dpi, including eight up-regulated and six down-regulated miRNAs, while 12 DEmiRNAs were observed in the midgut at 10 dpi, including nine up-regulated and three down-regulated ones. Additionally, five DEmiRNAs were shared, while nine and seven DEmiRNAs were specifically expressed in midguts at 7 dpi and 10 dpi. Gene ontology analysis suggested some DEmiRNAs and corresponding target mRNAs were involved in various functions including immune system processes and response to stimulus. KEGG pathway analysis shed light on the potential functions of some DEmiRNAs in regulating target mRNAs engaged in material and energy metabolisms, cellular immunity and the humoral immune system. Further investigation demonstrated a complex regulation network between DEmiRNAs and their target mRNAs, with miR-598-y, miR-252-y, miR-92-x and miR-3654-y at the center. Our results can facilitate future exploration of the regulatory roles of miRNAs in host responses to N. ceranae, and provide potential candidates for further investigation of the molecular mechanisms underlying eastern honeybee-microsporidian interactions.

Genome-Wide Identification of Long Non-Coding RNAs and Their Regulatory Networks Involved in Apis mellifera ligustica Response to Nosema ceranae Infection

Long non-coding RNAs (lncRNAs) are a diverse class of transcripts that structurally resemble mRNAs but do not encode proteins, and lncRNAs have been proven to play pivotal roles in a wide range of biological processes in animals and plants. However, knowledge of expression patterns and potential roles of honeybee lncRNA response to Nosema ceranae infection is completely unknown. Here, we performed whole transcriptome strand-specific RNA sequencing of normal midguts of Apis mellifera ligustica workers (Am7CK, Am10CK) and N. ceranae-inoculated midguts (Am7T, Am10T), followed by comprehensive analyses using bioinformatic and molecular approaches. A total of 6353 A. m. ligustica lncRNAs were identified, including 4749 conserved lncRNAs and 1604 novel lncRNAs. These lncRNAs had minimal sequence similarities with other known lncRNAs in other species; however, their structural features were similar to counterparts in mammals and plants, including shorter exon and intron length, lower exon number, and lower expression level, compared with protein-coding transcripts. Further, 111 and 146 N. ceranae-responsive lncRNAs were identified from midguts at 7-days post-inoculation (dpi) and 10 dpi compared with control midguts. Twelve differentially expressed lncRNAs (DElncRNAs) were shared by Am7CK vs. Am7T and Am10CK vs. Am10T comparison groups, while the numbers of unique DElncRNAs were 99 and 134, respectively. Functional annotation and pathway analysis showed that the DElncRNAs may regulate the expression of neighboring genes by acting in cis and trans fashion. Moreover, we discovered 27 lncRNAs harboring eight known miRNA precursors and 513 lncRNAs harboring 2257 novel miRNA precursors. Additionally, hundreds of DElncRNAs and their target miRNAs were found to form complex competitive endogenous RNA (ceRNA) networks, suggesting that these DElncRNAs may act as miRNA sponges. Furthermore, DElncRNA-miRNA-mRNA networks were constructed and investigated, the results demonstrated that a portion of the DElncRNAs were likely to participate in regulating the host material and energy metabolism as well as cellular and humoral immune host responses to N. ceranae invasion. Our findings revealed here offer not only a rich genetic resource for further investigation of the functional roles of lncRNAs involved in the A. m. ligustica response to N. ceranae infection, but also a novel insight into understanding the host-pathogen interaction during honeybee microsporidiosis.

Dicer regulates Nosema ceranae proliferation in honeybees

Nosema ceranae is a microsporidian parasite that infects the honeybee midgut epithelium. The protein-coding gene Dicer is lost in most microsporidian genomes but is present in N. ceranae. By feeding infected honeybees with small interfering RNA targeting the N. ceranae gene coding Dicer (siRNA-Dicer), we found that N. ceranae spore loads were significantly reduced. In addition, over 10% of total parasite protein-coding genes showed significantly divergent expression profiles after siRNA-Dicer treatment. Parasite genes for cell proliferation, ABC transporters and hexokinase were downregulated at 3 days postinfection, a key point in the middle of parasite replication cycles. In addition, genes involved in metabolic pathways of honeybees and N. ceranae showed significant co-expression. Furthermore, the siRNA-Dicer treatment partly reversed the expression patterns of honeybee genes. The honeybee gene mucin-2-like showed significantly upregulation in the siRNA-Dicer group compared with the infection group continually at 4, 5 and 6 days postinfection, suggesting that the siRNA-Dicer feeding promoted the strength of the mucus barrier resulted from interrupted parasite proliferation. As the gene Dicer broadly regulates N. ceranae proliferation and honeybee metabolism, our data suggest the RNA interference pathway is an important infection strategy for N. ceranae.

Encephalitozoon cuniculi and Vittaforma corneae (Phylum Microsporidia) inhibit staurosporine-induced apoptosis in human THP-1 macrophages in vitro

Obligately intracellular microsporidia regulate their host cell life cycles, including apoptosis, but this has not been evaluated in phagocytic host cells such as macrophages that can facilitate infection but also can be activated to kill microsporidia. We examined two biologically dissimilar human-infecting microsporidia species, Encephalitozoon cuniculi and Vittaforma corneae, for their effects on staurosporine-induced apoptosis in the human macrophage-differentiated cell line, THP1. Apoptosis was measured after exposure of THP-1 cells to live and dead mature organisms via direct fluorometric measurement of Caspase 3, colorimetric and fluorometric TUNEL assays, and mRNA gene expression profiles using Apoptosis RT2 Profiler PCR Array. Both species of microsporidia modulated the intrinsic apoptosis pathway. In particular, live E. cuniculi spores inhibited staurosporine-induced apoptosis as well as suppressed pro-apoptosis genes and upregulated anti-apoptosis genes more broadly than V. corneae. Exposure to dead spores induced an opposite effect. Vittaforma corneae, however, also induced inflammasome activation via Caspases 1 and 4. Of the 84 apoptosis-related genes assayed, 42 (i.e. 23 pro-apoptosis, nine anti-apoptosis, and 10 regulatory) genes were more affected including those encoding members of the Bcl2 family, caspases and their regulators, and members of the tumour necrosis factor (TNF)/TNF receptor R superfamily.

Nosemosis control in European honey bees, Apis mellifera, by silencing the gene encoding Nosema ceranae polar tube protein 3

RNA interference (RNAi) is a post-transcriptional gene silencing mechanism triggered by double-stranded RNA (dsRNA) that is homologous in sequence to the silenced gene and is conserved in a wide range of eukaryotic organisms. The RNAi mechanism has provided unique opportunities for combating honey bee diseases caused by various parasites and pathogens. Nosema ceranae is a microsporidian parasite of European honey bees, Apis mellifera, and has been associated with honey bee colony losses in some regions of the world. Here we explored the possibility of silencing the expression of a N. ceranae putative virulence factor encoding polar tube protein 3 (ptp3) which is involved in host cell invasion as a therapeutic strategy for controlling Nosema parasites in honey bees. Our studies showed that the oral ingestion of a dsRNA corresponding to the sequences of N. ceranae ptp3 could effectively suppress the expression of the ptp3 gene in N. ceranae-infected bees and reduce Nosema load. In addition to the knockdown of ptp3 gene expression, ingestion of ptp3-dsRNA also led to improved innate immunity in bees infected with N. ceranae along with an improvement in physiological performance and lifespan compared with untreated control bees. These results strongly suggest that RNAi-based therapeutics hold real promise for the effective treatment of honey bee diseases in the future, and warrant further investigation.

Impact of the microsporidian Nosema ceranae on the gut epithelium renewal of the honeybee, Apis mellifera

The invasive microsporidian species, Nosema ceranae, causes nosemosis in honeybees and is suspected to be involved in Western honeybee (Apis mellifera) declines worldwide. The midgut of honeybees is the site of infection; the microsporidium can disturb the functioning of this organ and, thus, the bee physiology. Host defense against pathogens is not limited to resistance (i.e. the immune response) but also involves resilience. This process implies that the host can tolerate and repair damage inflicted by the infection- by the pathogen itself or by an excessive host immune response. Enterocyte damage caused by N. ceranae can be compensated by proliferation of intestinal stem cells (ISCs) that are under the control of multiple pathways. In the present study, we investigated the impact of N. ceranae on honeybee epithelium renewal by following the mitotic index of midgut stem cells during a 22-day N. ceranae infection. Fluorescence in situ hybridization (FISH) and immunostaining experiments were performed to follow the parasite proliferation/progression in the intestinal tissue, especially in the ISCs as they are key cells for the midgut homeostasis. We also monitored the transcriptomic profile of 7 genes coding for key proteins involved in pathways implicated in the gut epithelium renewal and homeostasis. We have shown for the first time that N. ceranae can negatively alter the gut epithelium renewal rate and disrupt some signaling pathways involved in the gut homeostasis. This alteration is correlated to a reduced longevity of N. ceranae-infected honeybees and we can assume that honeybee susceptibility to N. ceranae could be due to an impaired ability to repair gut damage.

Interactions Among Host-Parasite MicroRNAs During Nosema ceranae Proliferation in Apis mellifera

We previously identified microRNA (miRNA) from Nosema ceranae and found that knockdowns of transcripts for the parasite protein Dicer greatly reduce parasite reproduction. In order to study parasitic miRNA functions and identify candidate target genes, we fed honey bees infected with N. ceranae with small interfering RNA (siRNA) targeting the N. ceranae gene Dicer. We then deep-sequenced honey bee and N. ceranae miRNAs daily across a full 6-day proliferation cycle. We found seven honey bee and five N. ceranae miRNAs that were significantly differently expressed between the infection and siRNA-Dicer groups. N. ceranae miRNA showed potentially strong impacts on the N. ceranae transcriptome, where over 79% of the total protein coding genes were significantly correlated with one or more miRNAs. N. ceranae miRNAs also can regulate honey bee metabolism and immune response, given parasitic miRNAs were secreted into the cytoplasm. Our results suggest that N. ceranae miRNAs regulate both parasite and host gene expression, providing new insights for microsporidia parasitism evolution.

Effects of the gut parasite Nosema ceranae on honey bee physiology and behavior

The common and widespread parasite Nosema ceranae is considered a major threat to the Western honey bee at both the individual and colony levels. Several studies demonstrated that infection by this parasite may affect physiology, behavior, and survival of honey bees. N. ceranae infection impairs midgut integrity and alters the energy demand in honey bees. The infection can also significantly suppress the bee immune response and modify pheromone production in worker and queen honey bees leading to precocious foraging. However, the presence of N. ceranae is not systematically associated with colony weakening and honey bee mortality. This variability depends upon parasite or host genetics, nutrition, climate or interactions with other stressors such as environmental contaminants or other parasites.

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

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

Impact of Nosema ceranae and Nosema apis on individual worker bees of the two host species (Apis cerana and Apis mellifera) and regulation of host immune response

Nosema apis and Nosema ceranae are obligate intracellular microsporidian parasites infecting midgut epithelial cells of host adult honey bees, originally Apis mellifera and Apis cerana respectively. Each microsporidia cross-infects the other host and both microsporidia nowadays have a worldwide distribution. In this study, cross-infection experiments using both N. apis and N. ceranae in both A. mellifera and A. cerana were carried out to compare pathogen proliferation and impact on hosts, including host immune response. Infection by N. ceranae led to higher spore loads than by N. apis in both host species, and there was greater proliferation of microsporidia in A. mellifera compared to A. cerana. Both N. apis and N. ceranae were pathogenic in both host Apis species. N. ceranae induced subtly, though not significantly, higher mortality than N. apis in both host species, yet survival of A. cerana was no different to that of A. mellifera in response to N. apis or N. ceranae. Infections of both host species with N. apis and N. ceranae caused significant up-regulation of AMP genes and cellular mediated immune genes but did not greatly alter apoptosis-related gene expression. In this study, A. cerana enlisted a higher immune response and displayed lower loads of N. apis and N. ceranae spores than A. mellifera, suggesting it may be better able to defend itself against microsporidia infection. We caution against over-interpretation of our results, though, because differences between host and parasite species in survival were insignificant and because size differences between microsporidia species and between host Apis species may alternatively explain the differential proliferation of N. ceranae in A. mellifera.

Comparative analyses of co-evolving host-parasite associations reveal unique gene expression patterns underlying slavemaker raiding and host defensive phenotypes

The transition to parasitism is a drastic shift in lifestyle, involving rapid changes in gene structure, function, and expression. After the establishment of antagonistic relationships, parasites and hosts co-evolve through reciprocal adaptations, often resulting in evolutionary arms-races. Repeated evolution of social parasitism and slavery among Temnothorax ants allows us to examine those gene expression patterns that characterize slavemaker raiding and reciprocal host defensive phenotypes. Previous behavioural studies have established that raiding strategies between Temnothorax slavemakers diverge, while host defense portfolios shift similarly under parasite pressure. We are the first to confirm this at the molecular level, revealing that slavemaking species exhibit a wider variety of genes with species-specific patterns of expression within their raiding phenotypes, whereas expression similarity is commonly found during the non-raiding phenotype. Host species response to slavemaker aggression, however, is indicated by strong changes in the expression of a relatively few number genes. Additionally, the expression of individual genes such as Acyl-CoA-Delta(11) desaturase and Trypsin-7 is strongly associated with the raiding phenotype of all three slavemaking species. Here, we provide novel insight into the gene expression patterns associated with raiding and nest defense behavior in Temnothorax ants, suggesting lineage-specific evolutionary patterns among both slavemakers and hosts.

Disruption of oxidative balance in the gut of the western honeybee Apis mellifera exposed to the intracellular parasite Nosema ceranae and to the insecticide fipronil

The causes underlying the increased mortality of honeybee colonies remain unclear and may involve multiple stressors acting together, including both pathogens and pesticides. Previous studies suggested that infection by the gut parasite Nosema ceranae combined with chronic exposure to sublethal doses of the insecticide fipronil generated an increase in oxidative stress in the midgut of honeybees. To explore the impact of these two stressors on oxidative balance, we experimentally infected bees with N. ceranae and/or chronically exposed to fipronil at low doses for 22 days, and we measured soluble reactive oxygen species (ROS) and ROS damage by quantifying both protein and lipid oxidation in the midgut. Our results revealed a disruption of the oxidative balance, with a decrease in both the amount of ROS and ROS damage in the presence of the parasite alone. However, protein oxidation was significantly increased in the N. ceranae/fipronil combination, revealing an increase in oxidative damage and suggesting higher fipronil toxicity in infected bees. Furthermore, our results highlighted a temporal order in the appearance of oxidation events in the intestinal cells and revealed that all samples tended to undergo protein oxidation during ageing, regardless of treatment.

Characterization of Nosema ceranae Genetic Variants from Different Geographic Origins

In recent years, large-scale colony losses of honey bees (Apis mellifera) have been reported and the infection with the microsporidia Nosema ceranae has been involved. However, the effect of N. ceranae at the colony level and its role in colony losses vary in different geographic areas. This difference may be related to the presence of multiple N. ceranae genetic variants resulting in different biological consequences. In this study, we analyzed the genetic diversity of 75 N. ceranae samples obtained from 13 countries and Hawaii through inter-sequence single repetition (ISSR) and evaluated if two of these genetic variants triggered different immune responses when infecting Apis mellifera iberiensis. The genetic diversity analysis showed that 41% of the samples had the same DNA amplification pattern, including samples from most European countries except Spain, while the remaining samples showed high variability. Infection assays were performed to analyze the infection levels and the immune response of bees infected with N. ceranae from Spain and Uruguay. The infected bees presented similar infection levels, and both isolates downregulated the expression of abaecin, confirming the ability of the microsporidia to depress the immune response. Only N. ceranae from Uruguay downregulated the expression level of imd compared to control bees. On the other hand, both genetic variants triggered different expression levels of lysozyme. As imd and lysozyme play important roles in the response to pathogens, these results could reflect differences in the biological consequences of N. ceranae variants in A. mellifera infection.

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

BACKGROUND

Organisms typically face infection by diverse pathogens, and hosts are thought to have developed specific responses to each type of pathogen they encounter. The advent of transcriptomics now makes it possible to test this hypothesis and compare host gene expression responses to multiple pathogens at a genome-wide scale. Here, we performed a meta-analysis of multiple published and new transcriptomes using a newly developed bioinformatics approach that filters genes based on their expression profile across datasets. Thereby, we identified common and unique molecular responses of a model host species, the honey bee (Apis mellifera), to its major pathogens and parasites: the Microsporidia Nosema apis and Nosema ceranae, RNA viruses, and the ectoparasitic mite Varroa destructor, which transmits viruses.

RESULTS

We identified a common suite of genes and conserved molecular pathways that respond to all investigated pathogens, a result that suggests a commonality in response mechanisms to diverse pathogens. We found that genes differentially expressed after infection exhibit a higher evolutionary rate than non-differentially expressed genes. Using our new bioinformatics approach, we unveiled additional pathogen-specific responses of honey bees; we found that apoptosis appeared to be an important response following microsporidian infection, while genes from the immune signalling pathways, Toll and Imd, were differentially expressed after Varroa/virus infection. Finally, we applied our bioinformatics approach and generated a gene co-expression network to identify highly connected (hub) genes that may represent important mediators and regulators of anti-pathogen responses.

CONCLUSIONS

Our meta-analysis generated a comprehensive overview of the host metabolic and other biological processes that mediate interactions between insects and their pathogens. We identified key host genes and pathways that respond to phylogenetically diverse pathogens, representing an important source for future functional studies as well as offering new routes to identify or generate pathogen resilient honey bee stocks. The statistical and bioinformatics approaches that were developed for this study are broadly applicable to synthesize information across transcriptomic datasets. These approaches will likely have utility in addressing a variety of biological questions.

Microsporidia infection impacts the host cell's cycle and reduces host cell apoptosis

Intracellular parasites can alter the cellular machinery of host cells to create a safe haven for their survival. In this regard, microsporidia are obligate intracellular fungal parasites with extremely reduced genomes and hence, they are strongly dependent on their host for energy and resources. To date, there are few studies into host cell manipulation by microsporidia, most of which have focused on morphological aspects. The microsporidia Nosema apis and Nosema ceranae are worldwide parasites of honey bees, infecting their ventricular epithelial cells. In this work, quantitative gene expression and histology were studied to investigate how these two parasites manipulate their host’s cells at the molecular level. Both these microsporidia provoke infection-induced regulation of genes involved in apoptosis and the cell cycle. The up-regulation of buffy (which encodes a pro-survival protein) and BIRC5 (belonging to the Inhibitor Apoptosis protein family) was observed after infection, shedding light on the pathways that these pathogens use to inhibit host cell apoptosis. Curiously, different routes related to cell cycle were modified after infection by each microsporidia. In the case of N. apis, cyclin B1, dacapo and E2F2 were up-regulated, whereas only cyclin E was up-regulated by N. ceranae, in both cases promoting the G1/S phase transition. This is the first report describing molecular pathways related to parasite-host interactions that are probably intended to ensure the parasite’s survival within the cell.

Differential proteomics reveals novel insights into Nosema-honey bee interactions

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

Approaches and Challenges to Managing Nosema (Microspora: Nosematidae) Parasites in Honey Bee (Hymenoptera: Apidae) Colonies

The microsporidia Nosema apis (Zander) and Nosema ceranae (Fries) are common intestinal parasites in honey bee (Apis mellifera L.) colonies. Though globally prevalent, there are mixed reports of Nosema infection costs, with some regions reporting high parasite virulence and colony losses, while others

REPORT

high Nosema prevalence but few costs. Basic and applied studies are urgently needed to help beekeepers effectively manage Nosema spp., ideally through an integrated pest management approach that allows beekeepers to deploy multiple strategies to control Nosema when Nosema is likely to cause damage to the colonies, rather than using prophylactic treatments. Beekeepers need practical and affordable technologies that facilitate disease diagnosis and science-backed guidelines that recommend when, if at all, to treat infections. In addition, new treatment methods are needed, as there are several problems associated with the chemical use of fumagillin (the only currently extensively studied, but not globally available treatment) to control Nosema parasites. Though selective breeding of Nosema-resistant or tolerant bees may offer a long-term, sustainable solution to Nosema management, other treatments are needed in the interim. Furthermore, the validation of alternative treatment efficacy in field settings is needed along with toxicology assays to ensure that treatments do not have unintended, adverse effects on honey bees or humans. Finally, given variation in Nosema virulence, development of regional management guidelines, rather than universal guidelines, may provide optimal and cost-effective Nosema management, though more research is needed before regional plans can be developed.

Host-Parasite Interactions and Purifying Selection in a Microsporidian Parasite of Honey Bees

To clarify the mechanisms of Nosema ceranae parasitism, we deep-sequenced both honey bee host and parasite mRNAs throughout a complete 6-day infection cycle. By time-series analysis, 1122 parasite genes were significantly differently expressed during the reproduction cycle, clustering into 4 expression patterns. We found reactive mitochondrial oxygen species modulator 1 of the host to be significantly down regulated during the entire infection period. Our data support the hypothesis that apoptosis of honey bee cells was suppressed during infection. We further analyzed genome-wide genetic diversity of this parasite by comparing samples collected from the same site in 2007 and 2013. The number of SNP positions per gene and the proportion of non-synonymous substitutions per gene were significantly reduced over this time period, suggesting purifying selection on the parasite genome and supporting the hypothesis that a subset of N. ceranae strains might be dominating infection.