Immune gene expression in Bombus terrestris: signatures of infection despite strong variation among populations, colonies, and sister workers

Infection by Crithidia bombi increases relative abundance of Lactobacillus spp. in the gut of Bombus terrestris

Gut microbial communities confer protection against natural pathogens in important pollinators from the genera Bombus and Apis. In commercial species B. terrestris and B. impatiens, the microbiota increases their resistance to the common and virulent trypanosomatid parasite Crithidia bombi. However, the mechanisms by which gut microorganisms protect the host are still unknown. Here, we test two hypotheses: microbiota protect the host (1) through stimulation of its immune response or protection of the gut epithelium and (2) by competing for resources with the parasite inside the gut. To test them, we reduced the microbiota of workers and then rescued the microbial community by feeding them with microbiota supplements. We then exposed them to an infectious dose of C. bombi and characterised gene expression and gut microbiota composition. We examined the expression of three antimicrobial peptide genes and Mucin-5AC, a gene with a putative role in gut epithelium protection, using qPCR. Although a protective effect against C. bombi was observed in bumblebees with supplemented microbiota, we did not observe an effect of the microbiota on gene expression that could explain alone the protective effect observed. On the other hand, we found an increased relative abundance of Lactobacillus bacteria within the gut of infected workers and a negative correlation of this genus with Gilliamella and Snodgrassella genera. Therefore, our results point to a displacement of bumblebee endosymbionts by C. bombi that might be caused by competition for space and nutrients between the parasite and the microbiota within the gut.

Methods matter: the influence of method on infection estimates of the bumblebee parasite Crithidia bombi

The bumblebee gut parasite, Crithidia bombi, is widespread and prevalent in the field. Its interaction with Bombus spp. is a well-established epidemiological model. It is spread faecal-orally between colonies via the shared use of flowers when foraging. Accurately measuring the level of infection in bumblebees is important for assessing its distribution in the field, and also when conducting epidemiological experiments. Studies generally use 1 of 2 methods for measuring infection. One approach measures infection in faeces whereas the other method measures infection in guts. We tested whether the method of measuring infection affected the estimation of infection. Bumblebees were inoculated with a standardized inoculum and infection was measured 1 week later using either the faecal or gut method. We found that when the gut method was used to measure infection intensity estimates were significantly different to and approximately double those from the faecal method. These results have implications for the interpretation of previous study results and for the planning of future studies. Given the importance of bumblebees as pollinators, the impact of C. bombi on bumblebee health, and its use as an epidemiological model, we call on researchers to move towards consistent quantification of infections to enable future comparisons and meta-analyses of studies.

Genetic variation and microbiota in bumble bees cross-infected by different strains of C. bombi

The bumblebee Bombus terrestris is commonly infected by a trypanosomatid gut parasite Crithidia bombi. This system shows a striking degree of genetic specificity where host genotypes are susceptible to different genotypes of parasite. To a degree, variation in host gene expression underlies these differences, however, the effects of standing genetic variation has not yet been explored. Here we report on an extensive experiment where workers of twenty colonies of B. terrestris were each infected by one of twenty strains of C. bombi. To elucidate the host's genetic bases of susceptibility to infection (measured as infection intensity), we used a low-coverage (~2 x) genome-wide association study (GWAS), based on angsd, and a standard high-coverage (~15x) GWAS (with a reduced set from a 8 x 8 interaction matrix, selected from the full set of twenty). The results from the low-coverage approach remained ambiguous. The high-coverage approach suggested potentially relevant genetic variation in cell surface and adhesion processes. In particular, mucin, a surface mucoglycoprotein, potentially affecting parasite binding to the host gut epithelia, emerged as a candidate. Sequencing the gut microbial community of the same bees showed that the abundance of bacterial taxa, such as Gilliamella, Snodgrassella, or Lactobacillus, differed between 'susceptible' and 'resistant' microbiota, in line with earlier studies. Our study suggests that the constitutive microbiota and binding processes at the cell surface are candidates to affect infection intensity after the first response (captured by gene expression) has run its course. We also note that a low-coverage approach may not be powerful enough to analyse such complex traits. Furthermore, testing large interactions matrices (as with the full 20 x 20 combinations) for the effect of interaction terms on infection intensity seems to blur the specific host x parasite interaction effects, likely because the outcome of an infection is a highly non-linear process dominated by variation in individually different pathways of host defence (immune) responses.

Embryonic manipulations modulate differential expressions of heat shock protein, fatty acid metabolism, and antioxidant-related genes in the liver of heat-stressed broilers

In this study, the effects of in ovo feeding of γ-aminobutyric acid (GABA) and embryonic thermal manipulation (TM) on plasma biochemical parameters, organ weights, and hepatic gene expression in broilers exposed to cyclic heat stress (32 ± 1°C for 8 days) (HS) were investigated. A total of 175 chicks were assigned to five treatments: chicks hatched from control eggs (CON); chicks hatched from control eggs but exposed to HS (CON+HS); chicks hatched from eggs injected at 17.5 days of incubation with 0.6mL of 10% GABA and exposed to HS (G10+HS); chicks hatched from thermally manipulated eggs (39.6°C, 6h/d from embryonic days 10 to 18) and exposed to HS (TM+HS); chicks hatched from eggs that received both previous treatments during incubation and exposed to HS (G10+TM+HS). Results revealed that on day 36 post-hatch, hepatic NADPH oxidase 1 (P = 0.034) and 4 (P = 0.021) genes were downregulated in the TM+HS and G10+TM+HS compared to the CON+HS group. In addition, while acetyl-CoA carboxylase gene expression was reduced (P = 0.002) in the G10+TM group, gene expression of extracellular fatty acid-binding protein and peroxisome proliferator-activated receptor-γ was lower (P = 0.045) in the TM+HS group than in the CON+HS group. HS led to higher gene expression of heat shock protein 70 (HSP70) and 90 (HSP90) (P = 0.005, and P = 0.022). On the other hand, the TM+HS group exhibited lower expression of both HSP70 (P = 0.031) and HSP90 (P = 0.043) whereas the G10+TM+HS group had a reduced (P = 0.016) HSP90 expression compared to the CON+HS. MANOVA on different gene sets highlighted an overall lower (P = 0.034) oxidative stress and lower (P = 0.035) heat shock protein expression in the G10+TM+HS group compared to the CON+HS group. Taken together, the current results suggest that the combination of in ovo feeding of GABA with TM can modulate HSPs and antioxidant-related gene expression in heat-stressed broilers.

Poison or Potion: Effects of Sunflower Phenolamides on Bumble Bees and Their Gut Parasite

Simple Summary

Bee declines have been reported worldwide, partly due to parasite spread induced by human activities. However, bees may forage on specific floral resources to face parasite infection. Such natural resources are comparable to ‘natural pharmacies’ and may be favoured in bee conservation strategies. Consumption of sunflower pollen, despite being detrimental for larval development, has been recently shown to reduce the load of a widespread bumble bee gut parasite in the common eastern bumble bee. Although the underlying mechanisms remain unknown, it has been suggested that sunflower phenolamides—a family of molecules found in most flowering plants—may be responsible for such a reduction in parasite load. Here, we tested the impacts of sunflower phenolamides on healthy and infected buff-tailed bumble bees. Expectedly, sunflower pollen had harmful consequences on bumble bee development but surprisingly, it did not alter parasite load. By contrast, sunflower phenolamides had milder effects on bumble bee development but unexpectedly increased parasite load. Phenolamide effects may stem from the physiological stress they induced or the gut microbial community alteration they may have triggered. Since biological models and experimental framework differ greatly in related studies tackling plant–bee–parasite interplays, we challenged the definition of medicinal effects and questioned the way to assess them in controlled conditions.

Abstract

Specific floral resources may help bees to face environmental challenges such as parasite infection, as recently shown for sunflower pollen. Whereas this pollen diet is known to be unsuitable for the larval development of bumble bees, it has been shown to reduce the load of a trypanosomatid parasite (Crithidia bombi) in the bumble bee gut. Recent studies suggested it could be due to phenolamides, a group of compounds commonly found in flowering plants. We, therefore, decided to assess separately the impacts of sunflower pollen and its phenolamides on a bumble bee and its gut parasite. We fed Crithidia-infected and -uninfected microcolonies of Bombus terrestris either with a diet of willow pollen (control), a diet of sunflower pollen (natural diet) or a diet of willow pollen supplemented with sunflower phenolamides (supplemented diet). We measured several parameters at both microcolony (i.e., food collection, parasite load, brood development and stress responses) and individual (i.e., fat body content and phenotypic variation) levels. As expected, the natural diet had detrimental effects on bumble bees but surprisingly, we did not observe any reduction in parasite load, probably because of bee species-specific outcomes. The supplemented diet also induced detrimental effects but by contrast to our a priori hypothesis, it led to an increase in parasite load in infected microcolonies. We hypothesised that it could be due to physiological distress or gut microbiota alteration induced by phenolamide bioactivities. We further challenged the definition of medicinal effects and questioned the way to assess them in controlled conditions, underlining the necessity to clearly define the experimental framework in this research field.

Between-population differences in constitutive and infection-induced gene expression in threespine stickleback

Vertebrate immunity is a complex system consisting of a mix of constitutive and inducible defences. Furthermore, host immunity is subject to selective pressure from a range of parasites and pathogens which can produce variation in these defences across populations. As populations evolve immune responses to parasites, they may adapt via a combination of (1) constitutive differences, (2) shared inducible responses, or (3) divergent inducible responses. Here, we leverage a powerful natural host‐parasite model system (Gasterosteus aculeatus and Schistochephalus solidus) to tease apart the relative contributions of these three types of adaptations to among‐population divergence in response to parasites. Gene expression analyses revealed limited evidence of significant divergence in constitutive expression of immune defence, and strong signatures of conserved inducible responses to the parasite. Furthermore, our results highlight a handful of immune‐related genes which show divergent inducible responses which may contribute disproportionately to functional differences in infection success or failure. In addition to investigating variation in evolutionary adaptation to parasite selection, we also leverage this unique data set to improve understanding of cellular mechanisms underlying a putative resistance phenotype (fibrosis). Combined, our results provide a case study in evolutionary immunology showing that a very small number of genes may contribute to genotype differences in infection response.

Transcriptome profiling of Lymnaea stagnalis (Gastropoda) for ecoimmunological research

BACKGROUND

Host immune function can contribute to numerous ecological/evolutionary processes. Ecoimmunological studies, however, typically use one/few phenotypic immune assays and thus do not consider the complexity of the immune system. Therefore, "omics" resources that allow quantifying immune activity across multiple pathways are needed for ecoimmunological models. We applied short-read based RNAseq (Illumina NextSeq 500, PE-81) to characterise transcriptome profiles of Lymnaea stagnalis (Gastropoda), a multipurpose model snail species. We used a genetically diverse snail stock and exposed individuals to immune elicitors (injury, bacterial/trematode pathogens) and changes in environmental conditions that can alter immune activity (temperature, food availability).

RESULTS

Immune defence factors identified in the de novo assembly covered elements broadly described in other gastropods. For instance, pathogen-recognition receptors (PRR) and lectins activate Toll-like receptor (TLR) pathway and cytokines that regulate cellular and humoral defences. Surprisingly, only modest diversity of antimicrobial peptides and fibrinogen related proteins were detected when compared with other taxa. Additionally, multiple defence factors that may contribute to the phenotypic immune assays used to quantify antibacterial activity and phenoloxidase (PO)/melanisation-type reaction in this species were found. Experimental treatments revealed factors from non-self recognition (lectins) and signalling (TLR pathway, cytokines) to effectors (e.g., antibacterial proteins, PO enzymes) whose transcription depended on immune stimuli and environmental conditions, as well as components of snail physiology/metabolism that may drive these effects. Interestingly, the transcription of many factors (e.g., PRR, lectins, cytokines, PO enzymes, antibacterial proteins) showed high among-individual variation.

CONCLUSIONS

Our results indicate several uniform aspects of gastropod immunity, but also apparent differences between L. stagnalis and some previously examined taxa. Interestingly, in addition to immune defence factors that responded to immune elicitors and changes in environmental conditions, many factors showed high among-individual variation across experimental snails. We propose that such factors are highly important to be included in future ecoimmunological studies because they may be the key determinants of differences in parasite resistance among individuals both within and between natural snail populations.

Workflow of Lotmaria passim isolation: Experimental infection with a low-passage strain causes higher honeybee mortality rates than the PRA-403 reference strain

The impact of trypanosomatid parasites on honeybee health may represent a major threat to bee colonies worldwide. However, few axenic isolates have been generated to date and with no details on cell culture passages, a parameter that could influence parasite virulence. To address this question, a trypanosomatid isolation protocol was developed and a new strain was obtained, named L. passim C1. Using experimental infection of worker honeybees, we compared the virulence and mortality rates of the ATCC PRA-403 reference strain and C1 strain, the latter showing higher virulence from 10 days post-infection onward. This study highlights the impact of cell culture passages on the pathogenicity of L. passim in honeybees, providing new evidence of its negative effects on honeybee health.

Multiple stressors interact to impair the performance of bumblebee Bombus terrestris colonies

Bumblebees are constantly exposed to a wide range of biotic and abiotic stresses which they must defend themselves against to survive. Pathogens and pesticides represent important stressors that influence bumblebee health, both when acting alone or in combination. To better understand bumblebee health, we need to investigate how these factors interact, yet experimental studies to date generally focus on only one or two stressors. The aim of this study is to evaluate how combined effects of four important stressors (the gut parasite Nosema ceranae, the neonicotinoid insecticide thiamethoxam, the pyrethroid insecticide cypermethrin and the EBI fungicide tebuconazole) interact to affect bumblebees at the individual and colony levels. We established seven treatment groups of colonies that we pulse exposed to different combinations of these stressors for 2 weeks under laboratory conditions. Colonies were subsequently placed in the field for 7 weeks to evaluate the effect of treatments on the prevalence of N. ceranae in inoculated bumblebees, expression levels of immunity and detoxification-related genes, food collection, weight gain, worker and male numbers, and production of worker brood and reproductives. Exposure to pesticide mixtures reduced food collection by bumblebees. All immunity-related genes were upregulated in the bumblebees inoculated with N. ceranae when they had not been exposed to pesticide mixtures, and bumblebees exposed to the fungicide and the pyrethroid were less likely to have N. ceranae. Combined exposure to the three-pesticide mixture and N. ceranae reduced bumblebee colony growth, and all treatments had detrimental effects on brood production. The groups exposed to the neonicotinoid insecticide produced 40%-76% fewer queens than control colonies. Our findings show that exposure to combinations of stressors that bumblebees frequently come into contact with have detrimental effects on colony health and performance and could therefore have an impact at the population level. These results also have significant implications for current practices and policies for pesticide risk assessment and use as the combinations tested here are frequently applied simultaneously in the field. Understanding the interactions between different stressors will be crucial for improving our ability to manage bee populations and for ensuring pollination services into the future.

A rank-based normalization method with the fully adjusted full-stage procedure in genetic association studies

In the area of genetic epidemiology, studies of the genotype-phenotype associations have made significant contributions to human complicated trait genetics. These studies depend on specialized statistical methods for uncover the association between traits and genetic variants, both common and rare variants. Often, in analyzing such studies, potentially confounding factors, such as social and environmental conditions, are required to be involved. Multiple linear regression is the most widely used type of regression analysis when the outcome of interest is quantitative traits. Many statistical tests for identifying genotype-phenotype associations using linear regression rely on the assumption that the traits (or the residuals) of the regression follow a normal distribution. In genomic research, the rank-based inverse normal transformation (INT) is one of the most popular approaches to reach normally distributed traits (or normally distributed residuals). Many researchers believe that applying the INT to the non-normality of the traits (or the non-normality of the residuals) is required for valid inference, because the phenotypic (or residual) outliers and non-normality have the significant influence on both the type I error rate control and statistical power, especially under the situation in rare-variant association testing procedures. Here we propose a test for exploring the association of the rare variant with the quantitative trait by using a fully adjusted full-stage INT. Using simulations we show that the fully adjusted full-stage INT is more appropriate than the existing INT methods, such as the fully adjusted two-stage INT and the INT-based omnibus test, in testing genotype-phenotype associations with rare variants, especially when genotypes are uncorrelated with covariates. The fully adjusted full-stage INT retains the advantages of the fully adjusted two-stage INT and ameliorates the problems of the fully adjusted two-stage INT for analysis of rare variants under non-normality of the trait. We also present theoretical results on these desirable properties. In addition, the two available methods with non-normal traits, the quantile/median regression method and the Yeo-Johnson power transformation, are also included in simulations for comparison with these desirable properties.

The Effect of Foraging on Bumble Bees, Bombus terrestris, Reared under Laboratory Conditions

Bumble bees are important pollinators broadly used by farmers in greenhouses and under conditions in which honeybee pollination is limited. As such, bumble bees are increasingly being reared for commercial purposes, which brings into question whether individuals reared under laboratory conditions are fully capable of physiological adaptation to field conditions. To understand the changes in bumble bee organism caused by foraging, we compared the fundamental physiological and immunological parameters of Bombus terrestris workers reared under constant optimal laboratory conditions with workers from sister colonies that were allowed to forage for two weeks in the field. Nutritional status and immune response were further determined in wild foragers of B.terrestris that lived under the constant influence of natural stressors. Both wild and laboratory-reared workers subjected to the field conditions had a lower protein concentration in the hemolymph and increased antimicrobial activity, the detection of which was limited in the non-foragers. However, in most of the tested parameters, specifically the level of carbohydrates, antioxidants, total hemocyte concentration in the hemolymph and melanization response, we did not observe any significant differences between bumble bee workers produced in the laboratory and wild animals, nor between foragers and non-foragers. Our results show that bumble bees reared under laboratory conditions can mount a sufficient immune response to potential pathogens and cope with differential food availability in the field, similarly to the wild bumble bee workers.

Female immunity in response to sexually transmitted opportunistic bacteria in the common bedbug Cimex lectularius

Besides typical sexually transmitted microbes, even environmental, opportunistic microbes have been found in copulatory organs of insects and even humans. To date, only one study has experimentally investigated the sexual transmission of opportunistic microbes from male to female insects, whereas nothing is known about the transmission from females to males. Even if opportunistic microbes do not cause infection upon transmission, they might eventually become harmful if they multiply inside the female. While the immune system of females is often assumed to target sexually transmitted microbes, most studies ignore the role of mating-associated opportunistic microbes. Variation in immunity between populations has been linked to parasite or bacteria prevalence but no study has ever addressed between-population differences in immune responses to sexually transmitted opportunistic microbes. We here show that bacteria applied to the copulatory organs of common bedbugs, Cimex lectularius, are sexually transmitted to the opposite sex at a high rate, including the transmission from female to male. Bacterial growth in the female sperm-receiving organ was inhibited over the first hours after introduction, but after this initial inhibition bacterial numbers increased, suggesting a shift of investment from immune defence towards reproduction. However, 24 h after the injection of bacteria, male components, or saline as a control, the sperm-receiving organ showed lysozyme-like activity and inhibited the growth of Gram-negative and Gram-positive bacteria in vitro, potentially to mop up the remaining bacteria. Contrasting our prediction, neither bacterial growth nor immune responses differed between populations. Future studies should link transmission dynamics, immune responses and fitness effects in both sexes. Experimental manipulation of environmental bacteria could be used to investigate how transmission frequency and toxicity of sexually transmitted opportunistic microbes shapes bacteria clearance and immune responses across populations.

Temperature dependence of parasitic infection and gut bacterial communities in bumble bees

High temperatures (e.g., fever) and gut microbiota can both influence host resistance to infection. However, effects of temperature-driven changes in gut microbiota on resistance to parasites remain unexplored. We examined the temperature dependence of infection and gut bacterial communities in bumble bees infected with the trypanosomatid parasite Crithidia bombi. Infection intensity decreased by over 80% between 21 and 37°C. Temperatures of peak infection were lower than predicted based on parasite growth in vitro, consistent with mismatches in thermal performance curves of hosts, parasites and gut symbionts. Gut bacterial community size and composition exhibited slight but significant, non-linear, and taxon-specific responses to temperature. Abundance of total gut bacteria and of Orbaceae, both negatively correlated with infection in previous studies, were positively correlated with infection here. Prevalence of the bee pathogen-containing family Enterobacteriaceae declined with temperature, suggesting that high temperature may confer protection against diverse gut pathogens. Our results indicate that resistance to infection reflects not only the temperature dependence of host and parasite performance, but also temperature-dependent activity of gut bacteria. The thermal ecology of gut parasite-symbiont interactions may be broadly relevant to infectious disease, both in ectothermic organisms that inhabit changing climates, and in endotherms that exhibit fever-based immunity.

Genomic Variation among Strains of Crithidia bombi and C. expoeki

In this study, we sequenced and analyzed the genomes of 40 strains, in addition to the already-reported two type strains, of two Crithidia species infecting bumblebees in Alaska and Central Europe and demonstrated that different strains of Crithidia bombi and C. expoeki vary considerably in terms of single nucleotide polymorphisms and gene copy number. Based on the genomic structure, phylogenetic analyses, and the pattern of copy number variation, we confirmed the status of C. expoeki as a separate species. The Alaskan populations appear to be clearly separated from those of Central Europe. This pattern fits a scenario of rapid host-parasite coevolution, where the selective advantage of a given parasite strain is only temporary. This study provides helpful insights into possible scenarios of selection and diversification of trypanosomatid parasites.IMPORTANCE A group of trypanosomatid flagellates includes several well-studied medically and economically important parasites of vertebrates and plants. Nevertheless, the vast majority of trypanosomatids infect only insects (mostly flies and true bugs) and, because of that, has attracted little research attention in the past. Of several hundred trypanosomatid species, only four can infect bees (honeybees and bumblebees). Because of such scarcity, these parasites are severely understudied. We analyzed whole-genome information for a total of 42 representatives of bee-infecting trypanosomatids collected in Central Europe and Alaska from a population genetics point of view. Our data shed light on the evolution, selection, and diversification in this important group of trypanosomatid parasites.

Bumblebee olfactory learning affected by task allocation but not by a trypanosome parasite

Parasites can induce behavioural changes in their host organisms. Several parasite species are known to infect bumblebees, an important group of pollinators. Task allocation within bumblebee colonies can also cause differences in behaviour. Thus, task allocation may lead to context-dependent impacts of parasites on host behaviour. This study uses Bombus terrestris and its gut trypanosome Crithidia bombi, to investigate the effects of parasitism, task allocation (foraging or nest-work) and their interactions, on olfactory learning. Prior to undergoing the olfactory learning task, bees were orally infected with a field-realistic dose of C. bombi, and observed to determine task allocation. Parasitism did not significantly affect olfactory learning, but task allocation did, with foragers being significantly more likely to learn than nest bees. There was no significant interaction between parasitism and task. These results suggest that C. bombi is unlikely to affect pollination services via changes in olfactory learning of its host if bees are under no environmental or nutritional stress. However, wild and commercial colonies are likely to face such stressors. Future studies in the field are needed to extrapolate our results to real world effects.

The genomes of Crithidia bombi and C. expoeki, common parasites of bumblebees

Trypanosomatids (Trypanosomatidae, Kinetoplastida) are flagellated protozoa containing many parasites of medical or agricultural importance. Among those, Crithidia bombi and C. expoeki, are common parasites in bumble bees around the world, and phylogenetically close to Leishmania and Leptomonas. They have a simple and direct life cycle with one host, and partially castrate the founding queens greatly reducing their fitness. Here, we report the nuclear genome sequences of one clone of each species, extracted from a field-collected infection. Using a combination of Roche 454 FLX Titanium, Pacific Biosciences PacBio RS, and Illumina GA2 instruments for C. bombi, and PacBio for C. expoeki, we could produce high-quality and well resolved sequences. We find that these genomes are around 32 and 34 MB, with 7,808 and 7,851 annotated genes for C. bombi and C. expoeki, respectively-which is somewhat less than reported from other trypanosomatids, with few introns, and organized in polycistronic units. A large fraction of genes received plausible functional support in comparison primarily with Leishmania and Trypanosoma. Comparing the annotated genes of the two species with those of six other trypanosomatids (C. fasciculata, L. pyrrhocoris, L. seymouri, B. ayalai, L. major, and T. brucei) shows similar gene repertoires and many orthologs. Similar to other trypanosomatids, we also find signs of concerted evolution in genes putatively involved in the interaction with the host, a high degree of synteny between C. bombi and C. expoeki, and considerable overlap with several other species in the set. A total of 86 orthologous gene groups show signatures of positive selection in the branch leading to the two Crithidia under study, mostly of unknown function. As an example, we examined the initiating glycosylation pathway of surface components in C. bombi, finding it deviates from most other eukaryotes and also from other kinetoplastids, which may indicate rapid evolution in the extracellular matrix that is involved in interactions with the host. Bumble bees are important pollinators and Crithidia-infections are suspected to cause substantial selection pressure on their host populations. These newly sequenced genomes provide tools that should help better understand host-parasite interactions in these pollinator pathogens.

Perspectives on the evolutionary ecology of arthropod antimicrobial peptides

Antimicrobial peptides (AMPs) are important elements of the innate immune defence in multicellular organisms that target and kill microbes. Here, we reflect on the various points that are raised by the authors of the 11 contributions to a special issue of Philosophical Transactions on the 'evolutionary ecology of arthropod antimicrobial peptides'. We see five interesting topics emerging. (i) AMP genes in insects, and perhaps in arthropods more generally, evolve much slower than most other immune genes. One explanation refers to the constraints set by AMPs being part of a finely tuned defence system. A new view argues that AMPs are under strong stabilizing selection. Regardless, this striking observation still invites many more questions than have been answered so far. (ii) AMPs almost always are expressed in combinations and sometimes show expression patterns that are dependent on the infectious agent. While it is often assumed that this can be explained by synergistic interactions, such interactions have rarely been demonstrated and need to be studied further. Moreover, how to define synergy in the first place remains difficult and needs to be addressed. (iii) AMPs play a very important role in mediating the interaction between a host and its mutualistic or commensal microbes. This has only been studied in a very small number of (insect) species. It has become clear that the very same AMPs play different roles in different situations and hence are under concurrent selection. (iv) Different environments shape the physiology of organisms; especially the host-associated microbial communities should impact on the evolution host AMPs. Studies in social insects and some organisms from extreme environments seem to support this notion, but, overall, the evidence for adaptation of AMPs to a given environment is scant. (v) AMPs are considered or already developed as new drugs in medicine. However, bacteria can evolve resistance to AMPs. Therefore, in the light of our limited understanding of AMP evolution in the natural context, and also the very limited understanding of the evolution of resistance against AMPs in bacteria in particular, caution is recommended. What is clear though is that study of the ecology and evolution of AMPs in natural systems could inform many of these outstanding questions, including those related to medical applications and pathogen control.This article is part of the themed issue 'Evolutionary ecology of arthropod antimicrobial peptides'.

Larvae act as a transient transmission hub for the prevalent bumblebee parasite Crithidia bombi

Disease transmission networks are key for understanding parasite epidemiology. Within the social insects, structured contact networks have been suggested to limit the spread of diseases to vulnerable members of their society, such as the queen or brood. However, even these complex social structures do not provide complete protection, as some diseases, which are transmitted by workers during brood care, can still infect the brood. Given the high rate of feeding interactions that occur in a social insect colony, larvae may act as disease transmission hubs. Here we use the bumblebee Bombus terrestris and its parasite Crithidia bombi to determine the role of brood in bumblebee disease transmission networks. Larvae that were artificially inoculated with C. bombi showed no signs of infection seven days after inoculation. However, larvae that received either an artificial inoculation or a contaminated feed from brood-caring workers were able to transmit the parasite to naive workers. These results suggest that the developing brood is a potential route of intracolonial disease transmission and should be included when considering social insect disease transmission networks.

Immune response and gut microbial community structure in bumblebees after microbiota transplants

Microbial communities are a key component of host health. As the microbiota is initially 'foreign' to a host, the host's immune system should respond to its acquisition. Such variation in the response should relate not only to host genetic background, but also to differences in the beneficial properties of the microbiota. However, little is known about such interactions. Here, we investigate the gut microbiota of the bumblebee, Bombus terrestris, which has a protective function against the bee's natural trypanosome gut parasite, Crithidia bombi We transplanted 'resistant' and 'susceptible' microbiota into 'resistant' and 'susceptible' host backgrounds, and studied the activity of the host immune system. We found that bees from different resistance backgrounds receiving a microbiota differed in aspects of their immune response. At the same time, the elicited immune response also depended on the received microbiota's resistance phenotype. Furthermore, the microbial community composition differed between microbiota resistance phenotypes (resistant versus susceptible). Our results underline the complex feedback between the host's ability to potentially exert selection on the establishment of a microbial community and the influence of the microbial community on the host immune response in turn.

Sources of variance in immunological traits: evidence of congruent latitudinal trends across species

Among-population differences in immunological traits allow assessment of both evolutionary and plastic changes in organisms' resistance to pathogens. Such knowledge also provides information necessary to predict responses of such traits to environmental changes. Studies on latitudinal trends in insect immunity have so far yielded contradictory results, suggesting that multispecies approaches with highly standardised experimental conditions are needed. Here, we studied among-population differences of two parameters reflecting constitutive immunity-phenoloxidase (PO) and lytic activity, using common-garden design on three distantly related moth species represented by populations ranging from northern Finland to Georgia (Caucasus). The larvae were reared at different temperatures and on different host plants under a crossed factors experimental design. Haemolymph samples for measurement of immune status were taken from the larvae strictly synchronously. Clear among-population differences could be shown only for PO activity in one species (elevated activity in the northern populations). There was some indication that the cases of total absence of lytic activity were more common in southern populations. The effects of temperature, host and sex on the immunological traits studied remained highly species specific. Some evidence was found that lytic activity may be involved in mediating trade-offs between immunity and larval growth performance. In contrast, PO activity rarely covaried with fitness-related traits, and neither were the values of PO and lytic activity correlated with each other. The relatively inconsistent nature of the detected patterns suggests that studies on geographic differences in immunological traits should involve multiple species, and rely on several immunological indices if general trends are a point of interest.

Genomic evidence for population-specific responses to co-evolving parasites in a New Zealand freshwater snail

Reciprocal co-evolving interactions between hosts and parasites are a primary source of strong selection that can promote rapid and often population- or genotype-specific evolutionary change. These host-parasite interactions are also a major source of disease. Despite their importance, very little is known about the genomic basis of co-evolving host-parasite interactions in natural populations, especially in animals. Here, we use gene expression and sequence evolution approaches to take critical steps towards characterizing the genomic basis of interactions between the freshwater snail Potamopyrgus antipodarum and its co-evolving sterilizing trematode parasite, Microphallus sp., a textbook example of natural coevolution. We found that Microphallus-infected P. antipodarum exhibit systematic downregulation of genes relative to uninfected P. antipodarum. The specific genes involved in parasite response differ markedly across lakes, consistent with a scenario where population-level co-evolution is leading to population-specific host-parasite interactions and evolutionary trajectories. We also used an F_ST -based approach to identify a set of loci that represent promising candidates for targets of parasite-mediated selection across lakes as well as within each lake population. These results constitute the first genomic evidence for population-specific responses to co-evolving infection in the P. antipodarum-Microphallus interaction and provide new insights into the genomic basis of co-evolutionary interactions in nature.

No evidence of local adaptation of immune responses to Gyrodactylus in three-spined stickleback (Gasterosteus aculeatus)

Parasitism represents one of the most widespread lifestyles in the animal kingdom, with the potential to drive coevolutionary dynamics with their host population. Where hosts and parasites evolve together, we may find local adaptation. As one of the main host defences against infection, there is the potential for the immune response to be adapted to local parasites. In this study, we used the three-spined stickleback and its Gyrodactylus parasites to examine the extent of local adaptation of parasite infection dynamics and the immune response to infection. We took two geographically isolated host populations infected with two distinct Gyrodactylus species and performed a reciprocal cross-infection experiment in controlled laboratory conditions. Parasite burdens were monitored over the course of the infection, and individuals were sampled at multiple time points for immune gene expression analysis. We found large differences in virulence between parasite species, irrespective of host, and maladaptation of parasites to their sympatric host. The immune system responded to infection, with a decrease in expression of innate and Th1-type adaptive response genes in fish infected with the less virulent parasite, representing a marker of a possible resistance mechanism. There was no evidence of local adaptation in immune gene expression levels. Our results add to the growing understanding of the extent of host-parasite local adaptation, and demonstrate a systemic immune response during infection with a common ectoparasite. Further immunological studies using the stickleback-Gyrodactylus system can continue to contribute to our understanding of the function of the immune response in natural populations.

Chronic neonicotinoid pesticide exposure and parasite stress differentially affects learning in honeybees and bumblebees

Learning and memory are crucial functions which enable insect pollinators to efficiently locate and extract floral rewards. Exposure to pesticides or infection by parasites may cause subtle but ecologically important changes in cognitive functions of pollinators. The potential interactive effects of these stressors on learning and memory have not yet been explored. Furthermore, sensitivity to stressors may differ between species, but few studies have compared responses in different species. Here, we show that chronic exposure to field-realistic levels of the neonicotinoid clothianidin impaired olfactory learning acquisition in honeybees, leading to potential impacts on colony fitness, but not in bumblebees. Infection by the microsporidian parasite Nosema ceranae slightly impaired learning in honeybees, but no interactive effects were observed. Nosema did not infect bumblebees (3% infection success). Nevertheless, Nosema-treated bumblebees had a slightly lower rate of learning than controls, but faster learning in combination with neonicotinoid exposure. This highlights the potential for complex interactive effects of stressors on learning. Our results underline that one cannot readily extrapolate findings from one bee species to others. This has important implications for regulatory risk assessments which generally use honeybees as a model for all bees.

Sexual healing: mating induces a protective immune response in bumblebees

The prevalence of sexual, as opposed to clonal, reproduction given the many costs associated with sexual recombination has been an enduring question in evolutionary biology. In addition to these often discussed costs, there are further costs associated with mating, including the induction of a costly immune response, which leaves individuals prone to infection. Here, we test whether mating results in immune activation and susceptibility to a common, ecologically important, parasite of bumblebees. We find that mating does result in immune activation as measured by gene expression of known immune genes, but that this activation improves resistance to this parasite. We conclude that although mating can corrupt immunity in some systems, it can also enhance immunity in others.

Insect antimicrobial peptides act synergistically to inhibit a trypanosome parasite

The innate immune system provides protection from infection by producing essential effector molecules, such as antimicrobial peptides (AMPs) that possess broad-spectrum activity. This is also the case for bumblebees, Bombus terrestris, when infected by the trypanosome, Crithidia bombi Furthermore, the expressed mixture of AMPs varies with host genetic background and infecting parasite strain (genotype). Here, we used the fact that clones of C. bombi can be cultivated and kept as strains in medium to test the effect of various combinations of AMPs on the growth rate of the parasite. In particular, we used pairwise combinations and a range of physiological concentrations of three AMPs, namely Abaecin, Defensin and Hymenoptaecin, synthetized from the respective genomic sequences. We found that these AMPs indeed suppress the growth of eight different strains of C. bombi, and that combinations of AMPs were typically more effective than the use of a single AMP alone. Furthermore, the most effective combinations were rarely those consisting of maximum concentrations. In addition, the AMP combination treatments revealed parasite strain specificity, such that strains varied in their sensitivity towards the same mixtures. Hence, variable expression of AMPs could be an alternative strategy to combat highly variable infections.This article is part of the themed issue 'Evolutionary ecology of arthropod antimicrobial peptides'.

No effect of low-level chronic neonicotinoid exposure on bumblebee learning and fecundity

In recent years, many pollinators have declined in abundance and diversity worldwide, presenting a potential threat to agricultural productivity, biodiversity and the functioning of natural ecosystems. One of the most debated factors proposed to be contributing to pollinator declines is exposure to pesticides, particularly neonicotinoids, a widely used class of systemic insecticide. Also, newly emerging parasites and diseases, thought to be spread via contact with managed honeybees, may pose threats to other pollinators such as bumblebees. Compared to honeybees, bumblebees could be particularly vulnerable to the effects of stressors due to their smaller and more short-lived colonies. Here, we studied the effect of field-realistic, chronic clothianidin exposure and inoculation with the parasite Nosema ceranae on survival, fecundity, sugar water collection and learning using queenless Bombus terrestris audax microcolonies in the laboratory. Chronic exposure to 1 ppb clothianidin had no significant effects on the traits studied. Interestingly, pesticide exposure in combination with additional stress caused by harnessing bees for Proboscis Extension Response (PER) learning assays, led to an increase in mortality. In contrast to previous findings, the bees did not become infected by N. ceranae after experimental inoculation with the parasite spores, suggesting variability in host resistance or parasite virulence. However, this treatment induced a slight, short-term reduction in sugar water collection, potentially through stimulation of the immune system of the bees. Our results suggest that chronic exposure to 1 ppb clothianidin does not have adverse effects on bumblebee fecundity or learning ability.

Among-lake reciprocal transplants induce convergent expression of immune genes in threespine stickleback

Geographic variation in parasite communities can drive evolutionary divergence in host immune genes. However, biotic and abiotic environmental variation can also induce plastic differences in immune function among populations. At present, there is little information concerning the relative magnitudes of heritable vs. induced immune divergence in natural populations. We examined immune gene expression profiles of threespine stickleback (Gasterosteus aculeatus) from six lakes on Vancouver Island, British Columbia. Parasite community composition differs between lake types (large or small, containing limnetic- or benthic-like stickleback) and between watersheds. We observed corresponding differences in immune gene expression profiles among wild-caught stickleback, using a set of seven immune genes representing distinct branches of the immune system. To evaluate the role of environmental effects on this differentiation, we experimentally transplanted wild-caught fish into cages in their native lake, or into a nearby foreign lake. Transplanted individuals' immune gene expression converged on patterns typical of their destination lake, deviating from their native expression profile. Transplant individuals' source population had a much smaller effect, suggesting relatively weak genetic underpinning of population differences in immunity, as viewed through gene expression. This strong environmental regulation of immune gene expression provides a counterpoint to the large emerging literature documenting microevolution and genetic diversification of immune function. Our findings illustrate the value of studying immunity in natural environmental settings where the immune system has evolved and actively functions.

Bee declines driven by combined stress from parasites, pesticides, and lack of flowers

Bees are subject to numerous pressures in the modern world. The abundance and diversity of flowers has declined; bees are chronically exposed to cocktails of agrochemicals, and they are simultaneously exposed to novel parasites accidentally spread by humans. Climate change is likely to exacerbate these problems in the future. Stressors do not act in isolation; for example, pesticide exposure can impair both detoxification mechanisms and immune responses, rendering bees more susceptible to parasites. It seems certain that chronic exposure to multiple interacting stressors is driving honey bee colony losses and declines of wild pollinators, but such interactions are not addressed by current regulatory procedures, and studying these interactions experimentally poses a major challenge. In the meantime, taking steps to reduce stress on bees would seem prudent; incorporating flower-rich habitat into farmland, reducing pesticide use through adopting more sustainable farming methods, and enforcing effective quarantine measures on bee movements are all practical measures that should be adopted. Effective monitoring of wild pollinator populations is urgently needed to inform management strategies into the future.

Protein-poor diet reduces host-specific immune gene expression in Bombus terrestris

Parasites infect hosts non-randomly as genotypes of hosts vary in susceptibility to the same genotypes of parasites, but this specificity may be modulated by environmental factors such as nutrition. Nutrition plays an important role for any physiological investment. As immune responses are costly, resource limitation should negatively affect immunity through trade-offs with other physiological requirements. Consequently, nutritional limitation should diminish immune capacity in general, but does it also dampen differences among hosts? We investigated the effect of short-term pollen deprivation on the immune responses of our model host Bombus terrestris when infected with the highly prevalent natural parasite Crithidia bombi. Bumblebees deprived of pollen, their protein source, show reduced immune responses to infection. They failed to upregulate a number of genes, including antimicrobial peptides, in response to infection. In particular, they also showed less specific immune expression patterns across individuals and colonies. These findings provide evidence for how immune responses on the individual-level vary with important elements of the environment and illustrate how nutrition can functionally alter not only general resistance, but also alter the pattern of specific host-parasite interactions.

Differential gene expression and alternative splicing in insect immune specificity

Background

Ecological studies routinely show genotype-genotype interactions between insects and their parasites. The mechanisms behind these interactions are not clearly understood. Using the bumblebee Bombus terrestris/trypanosome Crithidia bombi model system (two bumblebee colonies by two Crithidia strains), we have carried out a transcriptome-wide analysis of gene expression and alternative splicing in bees during C. bombi infection. We have performed four analyses, 1) comparing gene expression in infected and non-infected bees 24 hours after infection by Crithidia bombi, 2) comparing expression at 24 and 48 hours after C. bombi infection, 3) determining the differential gene expression associated with the bumblebee-Crithidia genotype-genotype interaction at 24 hours after infection and 4) determining the alternative splicing associated with the bumblebee-Crithidia genotype-genotype interaction at 24 hours post infection.

Results

We found a large number of genes differentially regulated related to numerous canonical immune pathways. These genes include receptors, signaling pathways and effectors. We discovered a possible interaction between the peritrophic membrane and the insect immune system in defense against Crithidia. Most interestingly, we found differential expression and alternative splicing of immunoglobulin related genes (Dscam and Twitchin) are associated with the genotype-genotype interactions of the given bumblebee colony and Crithidia strain.

Conclusions

In this paper we have shown that the expression and alternative splicing of immune genes is associated with specific interactions between different host and parasite genotypes in this bumblebee/trypanosome model.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1031) contains supplementary material, which is available to authorized users.

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.

Antimicrobial peptides play a functional role in bumblebee anti-trypanosome defense

Bumblebees, amongst the most important of pollinators, are under enormous population pressures. One of these is disease. The bumblebee and its gut trypanosome Crithidia bombi are one of the fundamental models of ecological immunology. Although there is previous evidence of increased immune gene expression upon Crithidia infection, recent work has focussed on the bumblebee's gut microbiota. Here, by knocking down gene expression using RNAi, we show for the first time that antimicrobial peptides (AMPs) have a functional role in anti-Crithidia defense.

Qualitatively different immune response of the bumblebee host, Bombus terrestris, to infection by different genotypes of the trypanosome gut parasite, Crithidia bombi

Insects have a complex and highly successful immune system that responds specifically to different types of parasites. Different genotypes of a parasite species can differ in infectivity and virulence; which is important for host-parasite co-evolutionary processes, such as antagonistic, fluctuating selection. Such coevolution obviously requires a genetic basis, but little is known about how hosts immunologically respond to different genotypes. The common European bumblebee Bombus terrestris is infected by the highly prevalent trypanosome gut parasite, Crithidia bombi. Here we examined expression of 26 immunological and metabolic genes in response to infection by two clones of C. bombi and compared that with exposure to injection with a bacterial challenge. Exposure to the two clones of C. bombi elicits qualitatively different immune expression responses. Interestingly, infection with one clone results in up regulation of AMP's similar to bees given the bacterial challenge, while genes related to metabolism, signalling, and other effectors were similar between the two Crithidia exposures. Bees given different challenges were distinct enough to discern using linear discriminant analyses. We also found strong correlations, both positive and negative, among genes, which may shed light on how suites of genes are regulated and trade-offs in expression within this gene set.