Molecular detection of protozoan parasites infecting Apis mellifera colonies in Japan

Effects of Nosema ceranae and Lotmaria passim infections on honey bee foraging behaviour and physiology

Nosema ceranae and Lotmaria passim are two commonly encountered digestive tract parasites of the adult honey bee (Apis mellifera L.). Although these parasites are associated with colony losses, little is known about how they affect individual bee physiology and behaviour at the colony level. Using locally obtained isolates, we investigated the effects of both single and mixed infections of L. passim and N. ceranae on honey bee vitellogenin (Vg) expression and foraging behaviour. At the first instance of foraging, bees inoculated with either parasite had significantly lower Vg expression than uninoculated bees, with bees from the mixed infection treatment having the lowest Vg expression. Bees from the mixed infection treatment also had significantly higher densities of N. ceranae spores and numerically greater densities of L. passim cells per bee compared with bees inoculated with either parasite alone. In addition, bees from the mixed infection treatment had a significantly younger average foraging age compared with uninoculated bees from the same cohort. Although we did not find any effect of treatment on foraging effort, we discovered that bees inoculated with L. passim alone, or together with N. ceranae, had higher returning rates of foragers than control bees or bees inoculated with N. ceranae alone. Our findings indicate that both parasites can alter individual bee physiology, leading to individual changes in behaviour that could alter colony foraging dynamics. These have the potential to result in smaller, less productive colonies, decreased colony survivorship and reduced income for beekeepers.

Somy evolution in the honey bee infecting trypanosomatid parasite Lotmaria passim

Lotmaria passim is a ubiquitous trypanosomatid parasite of honey bees nestled within the medically important subfamily Leishmaniinae. Although this parasite is associated with honey bee colony losses, the original draft genome-which was completed before its differentiation from the closely related Crithidia mellificae-has remained the reference for this species despite lacking improvements from newer methodologies. Here, we report the updated sequencing, assembly, and annotation of the BRL-type (Bee Research Laboratory) strain (ATCC PRA-422) of Lotmaria passim. The nuclear genome assembly has been resolved into 31 complete chromosomes and is paired with an assembled kinetoplast genome consisting of a maxicircle and 30 minicircle sequences. The assembly spans 33.7 Mb and contains very little repetitive content, from which our annotation of both the nuclear assembly and kinetoplast predicted 10,288 protein-coding genes. Analyses of the assembly revealed evidence of a recent chromosomal duplication event within chromosomes 5 and 6 and provided evidence for a high level of aneuploidy in this species, mirroring the genomic flexibility employed by other trypanosomatids as a means of adaptation to different environments. This high-quality reference can therefore provide insights into adaptations of trypanosomatids to the thermally regulated, acidic, and phytochemically rich honey bee hindgut niche, which offers parallels to the challenges faced by other Leishmaniinae during the challenges they undergo within insect vectors, during infection of mammals, and exposure to antiparasitic drugs throughout their multi-host life cycles. This reference will also facilitate investigations of strain-specific genomic polymorphisms, their role in pathogenicity, and the development of treatments for pollinator infection.

Phylogenetic analysis of the trypanosomatid parasite Lotmaria passim in honey bees (Apis mellifera) in Poland

Introduction

Lotmaria passim (L. passim) is a single-celled flagellate which colonises the bee gastrointestinal tract and is highly prevalent in honey bees. This parasite is associated with colony losses. Honey bee (Apis mellifera) colonies were sampled from five apiaries in the north-eastern part of Poland for the phylogenetic analysis of L. passim.

Material and Methods

Each apiary consisted of approximately 60 bee colonies, of which 20 were randomly selected. Samples of 60 differently aged worker bees were collected from each colony and pooled. A total of 100 bee colonies from five apiaries were examined. Protozoa of the Trypanosomatidae family were identified by PCR. L. passim was detected in 47 (47%) of the samples. The 18S ribosomal (r) RNA amplicons of L. passim were sequenced by a commercial service. Their sequences were analysed with BLASTN and noted to be compatible with the GenBank sequences of this region of the organism's genome. A sequence analysis was performed using the BioEdit Sequence Alignment Editor and Clustal W software.

Results

The amplicon sequences of L. passim were 100% homologous with the sequences deposited in GenBank under accession numbers KM066243.1., KJ684964.1 and KM980181.1.

Conclusion

This is the first study to perform a phylogenetic analysis of L. passim in Polish honey bees. The analysis demonstrated high levels of genetic similarity between isolates of L. passim colonising apiaries in the north-eastern region of Poland.

A tale of two parasites: Responses of honey bees infected with Nosema ceranae and Lotmaria passim

Nosema ceranae and Lotmaria passim are two commonly encountered digestive tract parasites of the honey bee that have been associated with colony losses in Canada, the United States, and Europe. Though honey bees can be co-infected with these parasites, we still lack basic information regarding how they impact bee health at the individual and colony level. Using locally-isolated parasite strains, we investigated the effect of single and co-infections of these parasites on individual honey bee survival, and their responsiveness to sucrose. Results showed that a single N. ceranae infection is more virulent than both single L. passim infections and co-infections. Honey bees singly infected with N. ceranae reached < 50% survival eight days earlier than those inoculated with L. passim alone, and four days earlier than those inoculated with both parasites. Honey bees infected with either one, or both, parasites had increased responsiveness to sucrose compared to uninfected bees, which could correspond to higher levels of hunger and increased energetic stress. Together, these findings suggest that N. ceranae and L. passim pose threats to bee health, and that the beekeeping industry should monitor for both parasites in an effort correlate pathogen status with changes in colony-level productivity and survival.

Frequent Parasitism of Apis mellifera by Trypanosomatids in Geographically Isolated Areas with Restricted Beekeeping Movements

Trypanosomatids form a group of high prevalence protozoa that parasitise honey bees, with Lotmaria passim as the predominant species worldwide. However, the knowledge about the ecology of trypanosomatids in isolated areas is limited. The Portuguese archipelagos of Madeira and Azores provide an interesting setting to investigate these parasites because of their geographic isolation, and because they harbour honey bee populations devoid of two major enemies: Varroa destructor and Nosema ceranae. Hence, a total of 661 honey bee colonies from Madeira and the Azores were analysed using different molecular techniques, through which we found a high prevalence of trypanosomatids despite the isolation of these islands. L. passim was the predominant species and, in most colonies, was the only one found, even on islands free of V. destructor and/or N. ceranae with severe restrictions on colony movements to prevent the spread of them. However, islands with V. destructor had a significantly higher prevalence of L. passim and, conversely, islands with N. ceranae did not shown any significant correlation with the trypanosomatid. Crithidia bombi was detected in Madeira and on three islands of the Azores, almost always coincident with L. passim. By contrast, Crithidia mellificae was not detected in any sample. A high-throughput sequencing analysis distinguished two main haplotypes of L. passim, which accounted for 98% of the total sequence reads. This work suggests that L. passim and C. bombi are parasites that have been associated with honey bees predating the spread of V. destructor and N. ceranae.

Trypanosomatids in honey bee colonies in Spain: A new specific qPCR method for specific quantification of Lotmaria passim, Crithidia mellificae and Crithidia bombi

Bee trypanosomatids have not been widely studied due to the original belief that these organisms were not pathogenic to honey bees. However, trypanosomatids have been linked to increased winter mortality in honey bee colonies in recent years and it has been shown that these pathogens can shorten a honey bee worker's lifespan in laboratory conditions. These studies found that this mortality corresponded to dose-dependent infection. Although Lotmaria passim is the most prevalent species worldwide, the natural load in colonies remains poorly investigated. Here we describe a new highly specific and sensitive qPCR method that allows the differentiation and quantification of the parasitic load of each of the three most common trypanosomatid species described to date in honey bee colonies: L. passim, Crithidia mellificae, and Crithidia bombi. We have used this new method to analyze honey bee colonies in central Spain and confirm that L. passim is the most common species and the one with higher parasitic loads in the colonies, which increased over the years, being higher in spring than in autumn. Crithidia mellificae was present along the study, with the highest prevalence in autumn 2019 and lately it was only found in non-quantifiable loads. Crithidia bombi was not detected in any of the colonies analyzed.

High temperatures augment inhibition of parasites by a honey bee gut symbiont

Temperature affects growth, metabolism, and interspecific interactions in microbial communities. Within animal hosts, gut bacterial symbionts can provide resistance to parasitic infections. Both infection and populations of symbionts can be shaped by the host body temperature. However, the effects of temperature on the antiparasitic activities of gut symbionts have seldom been explored. The Lactobacillus-rich gut microbiota of facultatively endothermic honey bees is subject to seasonal and ontogenetic changes in host temperature that could alter the effects of symbionts against parasites. We used cell cultures of a Lactobacillus symbiont and an important trypanosomatid gut parasite of honey bees to test the potential for temperature to shape parasite-symbiont interactions. We found that symbionts showed greater heat tolerance than parasites and chemically inhibited parasite growth via production of acids. Acceleration of symbiont growth and acid production at high temperatures resulted in progressively stronger antiparasitic effects across a temperature range typical of bee colonies. Consequently, the presence of symbionts reduced both the peak growth rate and heat tolerance of parasites. Substantial changes in parasite-symbiont interactions were evident over a temperature breadth that parallels changes in diverse animals exhibiting infection-related fevers and the amplitude of circadian temperature variation typical of endothermic birds and mammals, implying the frequent potential for temperature to alter symbiont-mediated resistance to parasites in endo- and ectothermic hosts. Results suggest that the endothermic behavior of honey bees could enhance the impacts of gut symbionts on parasites, implicating thermoregulation as a reinforcer of core symbioses and possibly microbiome-mediated antiparasitic defense. IMPORTANCE Two factors that shape the resistance of animals to infection are body temperature and gut microbiota. However, temperature can also alter interactions among microbes, raising the question of whether and how temperature changes the antiparasitic effects of gut microbiota. Honey bees are agriculturally important hosts of diverse parasites and infection-mitigating gut microbes. They can also socially regulate their body temperatures to an extent unusual for an insect. We show that high temperatures found in honey bee colonies augment the ability of a gut bacterial symbiont to inhibit the growth of a common bee parasite, reducing the parasite's ability to grow at high temperatures. This suggests that fluctuations in colony and body temperatures across life stages and seasons could alter the protective value of bees' gut microbiota against parasites, and that temperature-driven changes in gut microbiota could be an underappreciated mechanism by which temperature-including endothermy and fever-alters animal infection.

Microbiota and pathogens in an invasive bee: Megachile sculpturalis from native and invaded regions

The present study aimed to characterise the bacterial, fungal and parasite gut community of the invasive bee Megachile sculpturalis sampled from native (Japan) and invaded (USA and France) regions via 16S rRNA and ITS2 amplicon sequencing and PCR detection of bee microparasites. The bacterial and fungal gut microbiota communities in bees from invaded regions were highly similar and differed strongly from those obtained in Japan. Core amplicon sequence variants (ASVs) within each population represented environmental micro-organisms commonly present in bee-associated niches that likely provide beneficial functions to their host. Although the overall bacterial and fungal communities of the invasive M. sculpturalis in France and the co-foraging native bees Anthidium florentinum and Halictus scabiosae, were significantly different, five out of eight core ASVs were shared suggesting common environmental sources and potential transmission. None of the 46 M. sculpturalis bees analysed harboured known bee pathogens, while microparasite infections were common in A. florentinum, and rare in H. scabiosae. A common shift in the gut microbiota of M. sculpturalis in invaded regions as a response to changed environmental conditions, or a founder effect coupled to population re-establishment in the invaded regions may explain the observed microbial community profiles and the absence of parasites. While the role of pathogen pressure in shaping biological invasions is still debated, the absence of natural enemies may contribute to the invasion success of M. sculpturalis.

Molecular detection of Lotmaria passim in honeybees in Japan

Crithidia mellificae (C. mellificae) and Lotmaria passim (L. passim) are trypanosomatids that infect Apis mellifera. We analyzed the prevalence of C. mellificae and L. passim in six regions of Japan from 2018 to 2019. The detection rate of C. mellificae was 0.0% in all regions, whereas L. passim was detected in 16.7%-66.7% of the honeybees. L. passim was detected at a significantly lower rate in the Cyugoku-Shikoku region than in other regions. Furthermore, phylogenetic analysis of the internal transcribed spacer 1 (ITS1) locus of related species was performed. All the samples in this study could be assigned to the L. passim clade. This study reveals that L. passim infection is predominantly prevalent in Japan. Further epidemiological surveys are needed to clarify the prevalence of C. mellificae infection in honeybees in Japan.

Endosymbionts that threaten commercially raised and wild bumble bees (Bombus spp.)

Bumble bees (Bombus spp.) are important pollinators for both wild and agriculturally managed plants. We give an overview of what is known about the diverse community of internal potentially deleterious bumble bee symbionts, including viruses, bacteria, protozoans, fungi, and nematodes, as well as methods for their detection, quantification, and control. We also provide information on assessment of risk for select bumble bee symbionts and highlight key knowledge gaps. This information is crucial for ongoing efforts to establish parasite- conscious programs for future commerce in bumble bees for crop pollination, and to mitigate the problems with pathogen spillover to wild populations.

Population genetic diversity and dynamics of the honey bee brood pathogen Melissococcus plutonius in a region with high prevalence

European foulbrood (EFB) is a honey bee brood disease caused by the bacterium Melissococcus plutonius. Large-scale EFB outbreaks have been reported in several countries in recent decades, which entail costly sanitation measures of affected apiaries to restrict the spread of this contagious pathogen. To mitigate its impact, a better understanding of the population dynamics of the etiological agent is required. We here used multi-locus sequence typing (MLST) to infer the genetic diversity and geographical distribution of 160 M. plutonius isolates collected from EFB symptomatic honey bee colonies seven years apart. Isolates belonged to three clonal complexes (CCs) known worldwide and to 12 sequence types (STs), of which five were novel. Phylogenetic and clustering analyses showed that some of these novel sequence types have likely evolved locally during a period of outbreak, but most disappeared again. We further screened the isolates for melissotoxin A (mtxA), a putative virulence gene. The prevalence of STs in which mtxA was frequent increased over time, suggesting that this gene promotes spread. Despite the increased frequency of this gene in the population, the total number of cases decreased, which could be due to stricter control measures implemented before the second sampling period. Our results provide a better understanding of M. plutonius population dynamics and help identify knowledge gaps that limit efficient control of this emerging disease.

Molecular Detection and Differentiation of Arthropod, Fungal, Protozoan, Bacterial and Viral Pathogens of Honeybees

The honeybee Apis mellifera is highly appreciated worldwide because of its products, but also as it is a pollinator of crops and wild plants. The beehive is vulnerable to infections due to arthropods, fungi, protozoa, bacteria and/or viruses that manage to by-pass the individual and social immune mechanisms of bees. Due to the close proximity of bees in the beehive and their foraging habits, infections easily spread within and between beehives. Moreover, international trade of bees has caused the global spread of infections, several of which result in significant losses for apiculture. Only in a few cases can infections be diagnosed with the naked eye, by direct observation of the pathogen in the case of some arthropods, or by pathogen-associated distinctive traits. Development of molecular methods based on the amplification and analysis of one or more genes or genomic segments has brought significant progress to the study of bee pathogens, allowing for: (i) the precise and sensitive identification of the infectious agent; (ii) the analysis of co-infections; (iii) the description of novel species; (iv) associations between geno- and pheno-types and (v) population structure studies. Sequencing of bee pathogen genomes has allowed for the identification of new molecular targets and the development of specific genotypification strategies.

Presence of Known and Emerging Honey Bee Pathogens in Apiaries of Veneto Region (Northeast of Italy) during Spring 2020 and 2021

A progressive honey bee population decline has been reported worldwide during the last decades, and it could be attributed to several causes, in particular to the presence of pathogens and parasites that can act individually or in synergy. The health status of nine apiaries located in different areas of the Veneto region (northeast of Italy) was assessed for two consecutive years (2020 and 2021) in spring, during the resumption of honey bee activity, for determining the presence of known (Nosema spp., Varroa mite and viruses) and less known or emerging pathogens (Lotmaria passim and Crithidia mellificae) in honey bees. After honey bees sampling from each of the nine apiaries, Nosema apis, Nosema ceranae, L. passim, C. mellificae, ABPV, CBPV, IAPV, KBV, BQCV, SBV, DWV-A, DWV-B and V. destructor were investigated either by microscopic observation or PCR protocols. The viruses BQCV, SBV, CBPV followed by N. ceranae and L. passim were the most prevalent pathogens, and many of the investigated hives, despite asymptomatic, had different degrees of co-infection. This study aimed to highlight, during the resumption of honey bee activity in spring, the prevalence and spreading in the regional territory of different honey bee pathogens, which could alone or synergistically alter the homeostasis of bees colonies. The information gathered would increase our knowledge about the presence of these microorganisms and parasites in the territory and could contribute to improve beekeepers practice.

Hot and sour: parasite adaptations to honeybee body temperature and pH

Host temperature and gut chemistry can shape resistance to parasite infection. Heat and acidity can limit trypanosomatid infection in warm-blooded hosts and could shape infection resistance in insects as well. The colony-level endothermy and acidic guts of social bees provide unique opportunities to study how temperature and acidity shape insect-parasite associations. We compared temperature and pH tolerance between three trypanosomatid parasites from social bees and a related trypanosomatid from poikilothermic mosquitoes, which have alkaline guts. Relative to the mosquito parasites, all three bee parasites had higher heat tolerance that reflected body temperatures of hosts. Heat tolerance of the honeybee parasite Crithidia mellificae was exceptional for its genus, implicating honeybee endothermy as a plausible filter of parasite establishment. The lesser heat tolerance of the emerging Lotmaria passim suggests possible spillover from a less endothermic host. Whereas both honeybee parasites tolerated the acidic pH found in bee intestines, mosquito parasites tolerated the alkaline conditions found in mosquito midguts, suggesting that both gut pH and temperature could structure host-parasite specificity. Elucidating how host temperature and gut pH affect infection-and corresponding parasite adaptations to these factors-could help explain trypanosomatids' distribution among insects and invasion of mammals.

Crithidia mellificae infection in different mammalian species in Brazil

Crithidia mellificae, a monoxenous trypanosomatid considered restricted to insects, was recently reported to infect a bat. Herein, C. mellificae has been demonstrated to have a wider range of vertebrate hosts and distribution in Brazilian biomes than once thought. Parasites isolated from haemocultures were characterized using V7V8 SSU rDNA and glyceraldehyde 3-phosphate dehydrogenase genes. Coatis (Nasua nasua) in the Cerrado; marmosets (Callithrix sp.) and bats (Carollia perspicillata, Myotis lavali, M. izecksohni, Artibeus lituratus) in the Atlantic Forest; crab-eating foxes (Cerdocyon thous) and ocelot (Leopardus pardalis) in the Pantanal biomes were infected by trypanosomatids that displayed choanomastigote forms in haemoculture in Giemsa-stained slide smears. Molecular characterization and phylogenetic inference confirmed the infection of C. mellificae in these animals. Moreover, slight differences in C. mellificae sequences were observed. Crithidia mellificae growth curves were counted at 27°C, 36°C and 37°C, and the morphotypes were able to grow and survive for up to 16 days. Serological titers for C. mellificae were observed in nonhuman primates, demonstrating that this parasite is able to induce a humoral immune response in an infected mammal. These results showed that host specificity in trypanosomatids is complex and far from understood.

Putative determinants of virulence in Melissococcus plutonius, the bacterial agent causing European foulbrood in honey bees

Melissococcus plutonius

is a bacterial pathogen that causes epidemic outbreaks of European foulbrood (EFB) in honey bee populations. The pathogenicity of a bacterium depends on its virulence, and understanding the mechanisms influencing virulence may allow for improved disease control and containment. Using a standardized in vitro assay, we demonstrate that virulence varies greatly among sixteen M. plutonius isolates from five European countries. Additionally, we explore the causes of this variation. In this study, virulence was independent of the multilocus sequence type of the tested pathogen, and was not affected by experimental co-infection with Paenibacillus alvei, a bacterium often associated with EFB outbreaks. Virulence in vitro was correlated with the growth dynamics of M. plutonius isolates in artificial medium, and with the presence of a plasmid carrying a gene coding for the putative toxin melissotoxin A. Our results suggest that some M. plutonius strains showed an increased virulence due to the acquisition of a toxin-carrying mobile genetic element. We discuss whether strains with increased virulence play a role in recent EFB outbreaks.

Comparison between Apicystis cryptica sp. n. and Apicystis bombi (Arthrogregarida, Apicomplexa): Gregarine parasites that cause fat body hypertrophism in bees

The molecular divergence, morphology and pathology of a cryptic gregarine that is related to the bee parasite Apicystis bombi Lipa and Triggiani, 1996 is described. The 18S ribosomal DNA gene sequence of the new gregarine was equally dissimilar to that of A. bombi and the closest related genus Mattesia Naville, 1930, although phylogenetic analysis supported a closer relation to A. bombi. Pronounced divergence with A. bombi was found in the ITS1 sequence (69.6% similarity) and seven protein-coding genes (nucleotide 78.05% and protein 90.2% similarity). The new gregarine was isolated from a Bombus pascuorum Scopoli, 1763 female and caused heavy hypertrophism of the fat body tissue in its host. In addition, infected cells of the hypopharyngeal gland tissue, an important excretory organ of the host, were observed. Mature oocysts were navicular in shape and contained four sporozoites, similar to A. bombi oocysts. Given these characteristics, we proposed the name Apicystis cryptica sp. n. Detections so far indicated that distribution and host species occupation of Apicystis spp. overlap at least in Europe, and that historical detections could not discriminate between them. Specific molecular assays were developed that can be implemented in future pathogen screens that aim to discriminate Apicystis spp. in bees.

Trypanosomatid parasite dynamically changes the transcriptome during infection and modifies honey bee physiology

It is still not understood how honey bee parasite changes the gene expression to adapt to the host environment and how the host simultaneously responds to the parasite infection by modifying its own gene expression. To address this question, we studied a trypanosomatid, Lotmaria passim, which can be cultured in medium and inhabit the honey bee hindgut. We found that L. passim decreases mRNAs associated with protein translation, glycolysis, detoxification of radical oxygen species, and kinetoplast respiratory chain to adapt to the anaerobic and nutritionally poor honey bee hindgut during the infection. After the long term infection, the host appears to be in poor nutritional status, indicated by the increase and decrease of take-out and vitellogenin mRNAs, respectively. Simultaneous gene expression profiling of L. passim and honey bee during infection by dual RNA-seq provided insight into how both parasite and host modify their gene expressions.

The toxic unit approach as a risk indicator in honey bees surveillance programmes: A case of study in Apis mellifera iberiensis

The influence of genetic diversity and exposure to xenobiotics on the prevalence of pathogens was studied within the context of a voluntary epidemiological study in Spanish apiaries of Apis mellifera iberiensis, carried out during the spring season of years 2014 and 2015. As such, the evolutionary lineages of the honey bee colonies were identified, a multiresidue analysis of xenobiotics was carried out in beebread and worker bee samples, and the Toxic Unit (TUm) was estimated for each sampled apiary. The relationship between lineages and the most prevalent pathogens (Nosema ceranae, Varroa destructor, trypanosomatids, Black Queen Cell Virus; and Deformed Wing Virus) was analysed with contingency tables, and the possible relationships between TUm and the prevalence of these pathogens were studied by using a factor analysis. The statistical analysis supported the associations between V. destructor and Deformed Wing Virus (DWV), and between N. ceranae and Black Queen Cell Virus (BQCV), but the association between these pathogens and trypanosomatids was not observed. TUm values varied between 5.5 × 10^-6 and 3.65 × 10^-1. When TUm < 3.35 × 10^-4, it was mainly determined by coumaphos, tau-fluvalinate and/or chlorfenvinphos. At higher values, other insecticides also contributed to TUm, although a clear predominance was not seen up to TUm ≥ 1.83 × 10^-2, when it was mainly defined by acrinathrin, spinosad and/or imidacloprid. The possible cumulative effect from the joint action of xenobiotics was >10% in the 63% of the cases. The prevalence of pathogens did not appear to be influenced by the distribution of evolutionary lineages and, while the prevalence of V. destructor was not found to be determined by TUm, there was a trend towards an increasing prevalence of N. ceranae when TUm ≥ 23 10^-4. This study is an example of using TUm approach beyond the field of the ecotoxicology.

Gene Disruption of Honey Bee Trypanosomatid Parasite, Lotmaria passim, by CRISPR/Cas9 System

Two trypanosomatid species, Lotmaria passim and Crithidia mellificae, have been shown to parasitize honey bees to date. L. passim appears to be more prevalent than C. mellificae and specifically infects the honey bee hindgut. Although the genomic DNA has been sequenced, the effects of infection on honey bee health and colony are poorly understood. To identify the genes that are important for infecting honey bees and to understand their functions, we applied the CRISPR/Cas9 system to establish a method to manipulate L. passim genes. By electroporation of plasmid DNA and subsequent selection by drug, we first established an L. passim clone expressing tdTomato or Cas9. We also successfully disrupted the endogenous miltefosine transporter and tyrosine aminotransferase genes by replacement with drug (hygromycin) resistant gene using the CRISPR/Cas9-induced homology-directed repair pathway. The L. passim clone expressing fluorescent marker, as well as the simple method for editing specific genes, could become useful approaches to understand the underlying mechanisms of honey bee-trypanosomatid parasite interactions.

Detection of Lotmaria passim in Africanized and European honey bees from Uruguay, Argentina and Chile

Trypanosomatids affecting honey bees, Crithidia mellificae and Lotmaria passim, have been poorly studied in South America. We therefore analyzed their presence in Africanized and European honeybees from Uruguay, Argentina and Chile collected between 1990 and 2011 and assessed their association with other bee parasites and pathogens. Crithidia mellificae was not detected while L. passim was wide-spread. This report shows that L. passim has been present in this region at least since 2007 and it infects both Africanized and European honey bees. L. passim infected colonies showed high V. destructor parasitization levels, suggesting an association between them.

Early gut colonizers shape parasite susceptibility and microbiota composition in honey bee workers

Microbial symbionts living within animal guts are largely composed of resident bacterial species, forming communities that often provide benefits to the host. Gut microbiomes of adult honey bees (Apis mellifera) include core residents such as the betaproteobacterium Snodgrassella alvi, alongside transient parasites such as the protozoan Lotmaria passim To test how these species affect microbiome composition and host physiology, we administered S alvi and/or L passim inocula to newly emerged worker bees from four genetic backgrounds (GH) and reared them in normal (within hives) or stressed (protein-deficient, asocial) conditions. Microbiota acquired by normal bees were abundant but quantitatively differed across treatments, indicating treatment-associated dysbiosis. Pretreatment with S. alvi made normal bees more susceptible to L. passim and altered developmental and detoxification gene expression. Stressed bees were more susceptible to L. passim and were depauperate in core microbiota, yet supplementation with S. alvi did not alter this susceptibility. Microbiomes were generally more variable by GH in stressed bees, which also showed opposing and comparatively reduced modulation of gene expression responses to treatments compared with normal bees. These data provide experimental support for a link between altered gut microbiota and increased parasite and pathogen prevalence, as observed from honey bee colony collapse disorder.

Species-specific diagnostics of Apis mellifera trypanosomatids: A nine-year survey (2007-2015) for trypanosomatids and microsporidians in Serbian honey bees

In this study, honey bees collected in Serbia over 9 consecutive years (2007-2015) were retrospectively surveyed to determine the prevalence of eukaryotic gut parasites by molecular screening of archival DNA samples. We developed species-specific primers for PCR to detect the two known honey bee trypanosomatid species, Crithidia mellificae and the recently described Lotmaria passim. These primers were validated for target specificity under single and mixed-species conditions as well as against the bumblebee trypanosomatid Crithidia bombi. Infections by Nosema apis and Nosema ceranae (Microsporidia) were also determined using PCR. Samples from 162 colonies (18 from each year) originating from 57 different localities were surveyed. L. passim was detected in every year with an overall frequency of 62.3% and annual frequencies ranging from 38.9% to 83.3%. This provides the earliest confirmed record to date for L. passim and the first report of this species in Serbia. N. ceranae was ubiquitous, occurring in every year and at 95.7% overall frequency, ranging annually from 83.3% to 100%. The majority of colonies (60.5%) were co-infected with L. passim and N. ceranae, but colony infections by each species were statistically independent of one another over the nine years. Although C. mellificae and N. apis have both been reported recently at low frequency in Europe, neither of these species was detected in Serbia. These results support the hypothesis that L. passim has predominated over C. mellificae in A. mellifera during the past decade.

The effects of single and mixed infections of Apicystis bombi and deformed wing virus in Bombus terrestris

Many pollinators are currently suffering from declines, diminishing their gene pool and increasing their vulnerability to parasites. Recently, an increasing diversity of parasites has been recorded in bumblebees, yet for many, knowledge of their virulence and hence the risk their presence poses, is lacking. The deformed wing virus (DWV), known to be ubiquitous in honey bees, has now been detected in bumblebees. In addition, the neogregarine Apicystis bombi has been discovered to be more prevalent than previously thought. Here, we assess for the first time the lethal and sublethal effects of these parasites during single and mixed infections of worker bumblebees (Bombus terrestris). Fifteen days after experimental exposure, 22% of bees exposed to A. bombi, 50% of bees exposed to DWV and 86% of bees exposed to both parasites had died. Bumblebees that had ingested A. bombi had increased sucrose sensitivity (SS) and a lower lipid:body size ratio than control bees. While dual infected bumblebees showed no increase in SS. Overall, we find that A. bombi exhibits both lethal and sublethal effects. DWV causes lethal effect and may reduce the sub lethal effects imposed by A. bombi. The results show that both parasites have significant, negative effects on bumblebee health, making them potentially of conservation concern.

Hologenome theory and the honey bee pathosphere

Recent research has provided improved genome-level views of diversity across global honey bee populations, the gut microbiota residing within them, and the expanding pathosphere challenging honey bees. Different combinations of bee/microbiota/pathosphere genome complexes may explain regional variation in apiculture productivity and mortality. To understand this, we must consider management and research approaches in light of a hologenome paradigm: that honey bee fitness is determined by the composite bee and microbiota genomes. Only by considering the hologenome can we truly interpret and address impacts from the pathosphere, pesticides, toxins, nutrition, climate and other stressors affecting bee health.

Differential diagnosis of the honey bee trypanosomatids Crithidia mellificae and Lotmaria passim

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.

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

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

Metatranscriptomic analyses of honey bee colonies

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

Holistic screening of collapsing honey bee colonies in Spain: a case study

Background

Here we present a holistic screening of collapsing colonies from three professional apiaries in Spain. Colonies with typical honey bee depopulation symptoms were selected for multiple possible factors to reveal the causes of collapse.

Results

Omnipresent were Nosema ceranae and Lake Sinai Virus. Moderate prevalences were found for Black Queen Cell Virus and trypanosomatids, whereas Deformed Wing Virus, Aphid Lethal Paralysis Virus strain Brookings and neogregarines were rarely detected. Other viruses, Nosema apis, Acarapis woodi and Varroa destructor were not detected. Palinologic study of pollen demonstrated that all colonies were foraging on wild vegetation. Consequently, the pesticide residue analysis was negative for neonicotinoids. The genetic analysis of trypanosomatids GAPDH gene, showed that there is a large genetic distance between Crithidia mellificae ATCC30254, an authenticated cell strain since 1974, and the rest of the presumed C. mellificae sequences obtained in our study or published. This means that the latter group corresponds to a highly differentiated taxon that should be renamed accordingly.

Conclusion

The results of this study demonstrate that the drivers of colony collapse may differ between geographic regions with different environmental conditions, or with different beekeeping and agricultural practices. The role of other pathogens in colony collapse has to bee studied in future, especially trypanosomatids and neogregarines. Beside their pathological effect on honey bees, classification and taxonomy of these protozoan parasites should also be clarified.

Electronic supplementary material

The online version of this article (doi:10.1186/1756-0500-7-649) contains supplementary material, which is available to authorized users.

Widespread occurrence of honey bee pathogens in solitary bees

Solitary bees and honey bees from a neighbouring apiary were screened for a broad set of putative pathogens including protists, fungi, spiroplasmas and viruses. Most sampled bees appeared to be infected with multiple parasites. Interestingly, viruses exclusively known from honey bees such as Apis mellifera Filamentous Virus and Varroa destructor Macula-like Virus were also discovered in solitary bees. A microsporidium found in Andrena vaga showed most resemblance to Nosema thomsoni. Our results suggest that bee hives represent a putative source of pathogens for other pollinators. Similarly, solitary bees may act as a reservoir of honey bee pathogens.

A draft genome of the honey bee trypanosomatid parasite Crithidia mellificae

Since 2006, honey bee colonies in North America and Europe have experienced increased annual mortality. These losses correlate with increased pathogen incidence and abundance, though no single etiologic agent has been identified. Crithidia mellificae is a unicellular eukaryotic honey bee parasite that has been associated with colony losses in the USA and Belgium. C. mellificae is a member of the family Trypanosomatidae, which primarily includes other insect-infecting species (e.g., the bumble bee pathogen Crithidia bombi), as well as species that infect both invertebrate and vertebrate hosts including human pathogens (e.g.,Trypanosoma cruzi, T. brucei, and Leishmania spp.). To better characterize C. mellificae, we sequenced the genome and transcriptome of strain SF, which was isolated and cultured in 2010. The 32 megabase draft genome, presented herein, shares a high degree of conservation with the related species Leishmania major. We estimate that C. mellificae encodes over 8,300 genes, the majority of which are orthologs of genes encoded by L. major and other Leishmania or Trypanosoma species. Genes unique to C. mellificae, including those of possible bacterial origin, were annotated based on function and include genes putatively involved in carbohydrate metabolism. This draft genome will facilitate additional investigations of the impact of C. mellificae infection on honey bee health and provide insight into the evolution of this unique family.

Genetic variability of the neogregarine Apicystis bombi, an etiological agent of an emergent bumblebee disease

The worldwide spread of diseases is considered a major threat to biodiversity and a possible driver of the decline of pollinator populations, particularly when novel species or strains of parasites emerge. Previous studies have suggested that populations of introduced European honeybee (Apis mellifera) and bumblebee species (Bombus terrestris and Bombus ruderatus) in Argentina share the neogregarine parasite Apicystis bombi with the native bumblebee (Bombus dahlbomii). In this study we investigated whether A. bombi is acting as an emergent parasite in the non-native populations. Specifically, we asked whether A. bombi, recently identified in Argentina, was introduced by European, non-native bees. Using ITS1 and ITS2 to assess the parasite's intraspecific genetic variation in bees from Argentina and Europe, we found a largely unstructured parasite population, with only 15% of the genetic variation being explained by geographic location. The most abundant haplotype in Argentina (found in all 9 specimens of non-native species) was identical to the most abundant haplotype in Europe (found in 6 out of 8 specimens). Similarly, there was no evidence of structuring by host species, with this factor explaining only 17% of the genetic variation. Interestingly, parasites in native Bombus ephippiatus from Mexico were genetically distant from the Argentine and European samples, suggesting that sufficient variability does exist in the ITS region to identify continent-level genetic structure in the parasite. Thus, the data suggest that A. bombi from Argentina and Europe share a common, relatively recent origin. Although our data did not provide information on the direction of transfer, the absence of genetic structure across space and host species suggests that A. bombi may be acting as an emergent infectious disease across bee taxa and continents.

Comprehensive bee pathogen screening in Belgium reveals Crithidia mellificae as a new contributory factor to winter mortality

Since the last decade, unusually high honey bee colony losses have been reported mainly in North-America and Europe. Here, we report on a comprehensive bee pathogen screening in Belgium covering 363 bee colonies that were screened for 18 known disease-causing pathogens and correlate their incidence in summer with subsequent winter mortality. Our analyses demonstrate that, in addition to Varroa destructor, the presence of the trypanosomatid parasite Crithidia mellificae and the microsporidian parasite Nosema ceranae in summer are also predictive markers of winter mortality, with a negative synergy being observed between the two in terms of their effects on colony mortality. Furthermore, we document the first occurrence of a parasitizing phorid fly in Europe, identify a new fourth strain of Lake Sinai Virus (LSV), and confirm the presence of other little reported pathogens such as Apicystis bombi, Aphid Lethal Paralysis Virus (ALPV), Spiroplasma apis, Spiroplasma melliferum and Varroa destructor Macula-like Virus (VdMLV). Finally, we provide evidence that ALPV and VdMLV replicate in honey bees and show that viruses of the LSV complex and Black Queen Cell Virus tend to non-randomly co-occur together. We also noticed a significant correlation between the number of pathogen species and colony losses. Overall, our results contribute significantly to our understanding of honey bee diseases and the likely causes of their current decline in Europe.

Molecular and phylogenetic characterization of honey bee viruses, Nosema microsporidia, protozoan parasites, and parasitic mites in China

China has the largest number of managed honey bee colonies, which produce the highest quantity of honey and royal jelly in the world; however, the presence of honey bee pathogens and parasites has never been rigorously identified in Chinese apiaries. We thus conducted a molecular survey of honey bee RNA viruses, Nosema microsporidia, protozoan parasites, and tracheal mites associated with nonnative Apis mellifera ligustica and native Apis cerana cerana colonies in China. We found the presence of black queen cell virus (BQCV), chronic bee paralysis virus (CBPV), deformed wing virus (DWV), Israeli acute paralysis virus (IAPV), and sacbrood virus (SBV), but not that of acute bee paralysis virus (ABPV) or Kashmir bee virus (KBV). DWV was the most prevalent in the tested samples. Phylogenies of Chinese viral isolates demonstrated that genetically heterogeneous populations of BQCV, CBPV, DWV, and A. cerana-infecting SBV, and relatively homogenous populations of IAPV and A. meliifera-infecting new strain of SBV with single origins, are spread in Chinese apiaries. Similar to previous observations in many countries, Nosema ceranae, but not Nosema apis, was prevalent in the tested samples. Crithidia mellificae, but not Apicystis bombi was found in five samples, including one A. c. cerana colony, demonstrating that C. mellificae is capable of infecting multiple honey bee species. Based on kinetoplast-encoded cytochrome b sequences, the C. mellificae isolate from A. c. cerana represents a novel haplotype with 19 nucleotide differences from the Chinese and Japanese isolates from A. m. ligustica. This suggests that A. c. cerana is the native host for this specific haplotype. The tracheal mite, Acarapis woodi, was detected in one A. m. ligustica colony. Our results demonstrate that honey bee RNA viruses, N. ceranae, C. mellificae, and tracheal mites are present in Chinese apiaries, and some might be originated from native Asian honey bees.