How populations respond to climatic change depends on how climate affects individual growth, survival and reproduction; traits that ultimately shape population dynamics. However, climate can have complex effects on these traits and long-term data linking insect life-histories and climatic variation are rare. Here, we test how climate affects health, survival and reproduction in wild Bombus terrestris queens caught over 15 years after emerging from hibernation. We relate climate variation experienced by queens as they developed in their natal colonies and entered hibernation, to fitness traits assayed under constant lab conditions. We show that wet years consistently reduced queen fitness, while warm temperatures had positive and negative impacts. Behind these annual effects, lies strong seasonality. In particular, climatic conditions experienced by young queens as they forage, mate and enter hibernation are vital in determining whether they reproduce the following spring. Results suggest that these climate effects are partly due to resources: climatic drivers that reduce resource acquisition prior to hibernation, or disrupt diapause and thus accelerate resource loss over winter, are especially detrimental to spring queen fitness. While results show potential for increasingly negative effects of climate change on bumblebees, they also inform mitigation strategies; providing high-quality forage in late summer, before queens enter diapause. ### Competing Interest Statement The authors have declared no competing interest.

Climate change is a global biodiversity threat. To understand how a changing climate affects individual fitness and gain insights into which mechanisms are responsible, we need to establish how climatic variation affects individual life-history traits (e.g. growth, survival and reproduction). Long-term data linking insect life histories and climate parameters are therefore valuable but, unfortunately, rare. Here, we test how climatic variation affects health, survival and reproduction in a key European pollinator, the buff-tailed bumblebee Bombus terrestris. We relate climatic variation experienced by developing queens in the field, caught between 2000 and 2014, to fitness traits assayed in these queens under otherwise constant lab conditions. We show that wet years consistently reduce queen fitness, while warm temperatures have positive and negative impacts. Behind these annual effects lies strong seasonality. In particular, climatic conditions experienced by young queens as they forage, mate and enter hibernation are vital determinants of whether they reproduce the following spring. Results suggest that climatic drivers that reduce queen resource acquisition prior to hibernation, or accelerate resource loss over winter, are especially detrimental to spring queen fitness. This suggests a strategy to mitigate the negative effects of climate change on bumblebees: ensuring high-quality forage late in summer before queens enter diapause.

Leishbuviridae (Bunyavirales) are a diverse monophyletic group of negative-sense single-stranded RNA virus infecting parasitic flagellates of the family Trypanosomatidae. The presence of RNA viruses in trypanosomatids can influence the virulence of the latter. Here, we performed a screening for viruses in Crithidia bombi - a common parasite of important pollinators Bombus spp. (bumblebees) that negatively affects its host in stressful conditions. The majority (8/10) of C. bombi isolates collected in Europe and North America were positive for a virus that we named Crithidia bombi leishbuvirus 1 with high conservation of amino acid sequences between isolates. The results of our comparative phylogenetic analyses of the trypanosomatids and their viruses suggest that the high mobility of bumblebees and frequent coinfections by different strains of C. bombi determine an extensive viral exchange between the latter.

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.

Climate change is predicted to affect host-parasite interactions, and for some hosts, parasite infection is expected to increase with rising temperatures. Global population declines of important pollinators already have been attributed to climate change and parasitism. However, the role of climate in driving parasite infection and the genetic basis for pollinator hosts to respond often remain obscure. Based on decade-long field data, we investigated the association between climate and Nosema bombi (Microsporidia) infection of buffed-tailed bumblebees (Bombus terrestris), and whether host genotypes play a role. For this, we genotyped 876 wild bumblebee queens and screened for N. bombi infection of those queens between 2000 and 2010. We recorded seven climate parameters during those 11 years and tested for correlations between climate and infection prevalence. Here we show that climatic factors drive N. bombi infection and that the impact of climate depends on mitochondrial DNA cytochrome oxidase I (COI) haplotypes of the host. Infection prevalence was correlated with climatic variables during the time when queens emerge from hibernation. Remarkably, COI haplotypes best predict this association between climatic factors and infection. In particular, two host haplotypes ("A" and "B") displayed phenotypic plasticity in response to climatic variation: Temperature was positively correlated with infection of host haplotype B, but not haplotype A. The likelihood of infection of haplotype A was associated with moisture, conferring greater resistance to parasite infection during wetter years. In contrast, infection of haplotype B was unrelated to moisture. To the best of our knowledge, this is the first study that identifies specific host genotypes that confer differential parasite resistance under variable climatic conditions. Our results underscore the importance of mitochondrial haplotypes to ward off parasites in a changing climate. More broadly, this also suggests that COI may play a pertinent role in climate change adaptations of insect pollinators.

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.

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.

The antagonistic relationship between parasites and their hosts is strongly influenced by genotype-by-genotype interactions. Defense against parasitism is commonly studied in the context of immune system-based mechanisms and, thus, the focus in the search for candidate genes in host-parasite interactions is often on immune genes. In this study, we investigated the association between prevalence of parasite infection and host mitochondrial DNA (mtDNA) haplotypes in two natural populations of bumblebees (Bombus terrestris). The two most common haplotypes of the host populations, termed A and B, differ by a single nonsynonymous nucleotide substitution within the coding region of cytochrome oxidase I, an important player in metabolic pathways. We screened infection by Nosema bombi, a common endoparasite of bumblebees, and the corresponding host mtDNA-haplotype frequencies in over 1400 bumblebees between 2000 and 2010. The island population of Gotland showed lower mtDNA diversity compared to the mainland population in Switzerland. Over time, we observed large fluctuations in infection prevalence, as well as variation in host haplotype frequencies in both populations. Our long-term observation revealed that N. bombi infection of specific host genotypes is transient: We found that with increasing infection prevalence, proportionally more individuals with haplotype B, but fewer individuals with haplotype A were infected. This suggests that the presence of N. bombi in specific host genotypes relates to infection prevalence. This may be a result of parasite competition, or differential resilience of host types to ward off infections. The findings highlight the important role of host mtDNA haplotypes in the interaction with parasites.

A total of 1940 isolates from gut samples of 60 bumblebees representing Bombus pascuorum, Bombus terrestris, Bombus lucorum and Bombus lapidarius was collected and identified through state-of the-art taxonomic methods. The bacterial species diversity in these Bombus species exceeded that suggested by phylotype analysis through 16S rRNA amplicon sequencing, and revealed that B. pascuorum and B. terrestris had a unique microbiota composition, each. Representatives of most phylotypes reported earlier and detected in the present study were effectively isolated, and included several novel bacterial taxa and species reported for the first time in the bumblebee gut. Isolates were screened in pectin degradation assays and growth inhibition assays against the honeybee pathogens Paenibacillus larvae, Melissococcus plutonius and Ascosphaera apis and the bumblebee parasite Crithidia bombi. While inhibitory activity against each of these pathogens was observed, only one single culture was able to degrade pectin and polygalacturonic acid in vitro. The availability of accurately identified microbial isolates will facilitate future evaluation of the functional potential of the bumblebee gut microbiota.

Bumblebees are highly valued for their pollination services in natural ecosystems as well as for agricultural crops. These precious pollinators are known to be declining worldwide, and one major factor contributing to this decline are infections by parasites. Knowledge about parasites in wild bumblebee populations is thus of paramount importance for conservation purposes. We here report the geographical distribution of Crithidia and Nosema, two common parasites of bumblebees, in a yet poorly investigated country: Mexico. Based on sequence divergence of the Cytochrome b and Glycosomal glyceraldehyde phosphate deshydrogenase (gGPDAH) genes, we discovered the presence of a new Crithidia species, which is mainly distributed in the southern half of the country. It is placed by Bayesian inference as a sister species to C. bombi. We suggest the name Crithidia mexicana for this newly discovered organism. A population of C. expoeki was encountered concentrated on the flanks of the dormant volcanic mountain, Iztaccihuatl, and microsatellite data showed evidence of a bottleneck in this population. This study is the first to provide a large-scale insight into the health status of endemic bumblebees in Mexico, based on a large sample size (n = 3,285 bees examined) over a variety of host species and habitats.

ABSTRACT As pollinators, bees are cornerstones for terrestrial ecosystem stability and key components in agricultural productivity. All animals, including bees, are associated with a diverse community of microbes, commonly referred to as the microbiome. The bee microbiome is likely to be a crucial factor affecting host health. However, with the exception of a few pathogens, the impacts of most members of the bee microbiome on host health are poorly understood. Further, the evolutionary and ecological forces that shape and change the microbiome are unclear. Here, we discuss recent progress in our understanding of the bee microbiome, and we present challenges associated with its investigation. We conclude that global coordination of research efforts is needed to fully understand the complex and highly dynamic nature of the interplay between the bee microbiome, its host, and the environment. High-throughput sequencing technologies are ideal for exploring complex biological systems, including host-microbe interactions. To maximize their value and to improve assessment of the factors affecting bee health, sequence data should be archived, curated, and analyzed in ways that promote the synthesis of different studies. To this end, the BeeBiome consortium aims to develop an online database which would provide reference sequences, archive metadata, and host analytical resources. The goal would be to support applied and fundamental research on bees and their associated microbes and to provide a collaborative framework for sharing primary data from different research programs, thus furthering our understanding of the bee microbiome and its impact on pollinator health.

BACKGROUND:Sociality has many rewards, but can also be dangerous, as high population density and low genetic diversity, common in social insects, is ideal for parasite transmission. Despite this risk, honeybees and other sequenced social insects have far fewer canonical immune genes relative to solitary insects. Social protection from infection, including behavioral responses, may explain this depauperate immune repertoire. Here, based on full genome sequences, we describe the immune repertoire of two ecologically and commercially important bumblebee species that diverged approximately 18 million years ago, the North American Bombus impatiens and European Bombus terrestris.RESULTS:We find that the immune systems of these bumblebees, two species of honeybee, and a solitary leafcutting bee, are strikingly similar. Transcriptional assays confirm the expression of many of these genes in an immunological context and more strongly in young queens than males, affirming Bateman's principle of greater investment in female immunity. We find evidence of positive selection in genes encoding antiviral responses, components of the Toll and JAK/STAT pathways, and serine protease inhibitors in both social and solitary bees. Finally, we detect many genes across pathways that differ in selection between bumblebees and honeybees, or between the social and solitary clades.CONCLUSIONS:The similarity in immune complement across a gradient of sociality suggests that a reduced immune repertoire predates the evolution of sociality in bees. The differences in selection on immune genes likely reflect divergent pressures exerted by parasites across social contexts.

Background: The shift from solitary to social behavior is one of the major evolutionary transitions. Primitively eusocial bumblebees are uniquely placed to illuminate the evolution of highly eusocial insect societies. Bumblebees are also invaluable natural and agricultural pollinators, and there is widespread concern over recent population declines in some species. High-quality genomic data will inform key aspects of bumblebee biology, including susceptibility to implicated population viability threats.Results: We report the high quality draft genome sequences of Bombus terrestris and Bombus impatiens, two ecologically dominant bumblebees and widely utilized study species. Comparing these new genomes to those of the highly eusocial honeybee Apis mellifera and other Hymenoptera, we identify deeply conserved similarities, as well as novelties key to the biology of these organisms. Some honeybee genome features thought to underpin advanced eusociality are also present in bumblebees, indicating an earlier evolution in the bee lineage. Xenobiotic detoxification and immune genes are similarly depauperate in bumblebees and honeybees, and multiple categories of genes linked to social organization, including development and behavior, show high conservation. Key differences identified include a bias in bumblebee chemoreception towards gustation from olfaction, and striking differences in microRNAs, potentially responsible for gene regulation underlying social and other traits.Conclusions: These two bumblebee genomes provide a foundation for post-genomic research on these key pollinators and insect societies. Overall, gene repertoires suggest that the route to advanced eusociality in bees was mediated by many small changes in many genes and processes, and not by notable expansion or depauperation.

Trypanosomatids infecting honey bees have been poorly studied with molecular methods until recently. After the description of Crithidia mellificae (Langridge and McGhee, 1967) it took about forty years until molecular data for honey bee trypanosomatids became available and were used to identify and describe a new trypanosomatid species from honey bees, Lotmaria passim (Evans and Schwarz, 2014). However, an easy method to distinguish them without sequencing is not yet available. Research on the related bumble bee parasites Crithidia bombi and Crithidia expoeki revealed a fragment length polymorphism in the internal transcribed spacer 1 (ITS1), which enabled species discrimination. In search of fragment length polymorphisms for differential diagnostics in honey bee trypanosomatids, we studied honey bee trypanosomatid cell cultures of C. mellificae and L. passim. This research resulted in the identification of fragment length polymorphisms in ITS1 and ITS1-2 markers, which enabled us to develop a diagnostic method to differentiate both honey bee trypanosomatid species without the need for sequencing. However, the amplification success of the ITS1 marker depends probably on the trypanosomatid infection level. Further investigation confirmed that L. passim is the dominant species in Belgium, Japan and Switzerland. We found C. mellificae only rarely in Belgian honey bee samples, but not in honey bee samples from other countries. C. mellificae was also detected in mason bees (Osmia bicornis and Osmia cornuta) besides in honey bees. Further, the characterization and comparison of additional markers from L. passim strain SF (published as C. mellificae strain SF) and a Belgian honey bee sample revealed very low divergence in the 18S rRNA, ITS1-2, 28S rRNA and cytochrome b sequences. Nevertheless, a variable stretch was observed in the gp63 virulence factor.