Recent worldwide expansion of Nosema ceranae (Microsporidia) in Apis mellifera populations inferred from multilocus patterns of genetic variation

Geographic population structure of the honeybee microsporidian parasite Vairimorpha (Nosema) ceranae in the South West Indian Ocean

The microsporidian Vairimorpha (Nosema) ceranae is one of the most common parasites of the honeybee. A single honeybee carries many parasites and therefore multiple alleles of V. ceranae genes that seem to be ubiquitous. As a consequence, nucleotide diversity analyses have not allowed discriminating genetic structure of parasite populations. We performed deep loci-targeted sequencing to monitor the haplotype frequencies of genome markers in isolates from discontinuous territories, namely the tropical islands of the South West Indian Ocean. The haplotype frequency distribution corroborated the suspected tetraploidy of the parasite. Most major haplotypes were ubiquitous in the area but with variable frequency. While oceanic isolates differed from European and Asian outgroups, parasite populations from distinct archipelagoes also differed in their haplotype distribution. Interestingly an original and very divergent Malagasy isolate was detected. The observed population structure allowed formulating hypotheses upon the natural history of V. ceranae in this oceanic area. We also discussed the usefulness of allelic distribution assessment, using multiple informative loci or genome-wide analyses, when parasite population is not clonal within a single host.

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

BACKGROUND

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

OBJECTIVE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Morphological characterization and genetic diversity of a new microsporidium, Neoflabelliforma dubium n. sp. from the adipose tissue of Diaphanosoma dubium (Crustacea: Sididae)

We reported a new microsporidium Neoflabelliforma dubium n. sp. from the adipose tissue of Diaphanosoma dubium in China. The infected daphnids generally appeared opaque due to the presence of numerous spore aggregates located in the adipose tissue. All developmental stages were in direct contact with the host cell cytoplasm. Multinucleate sporogonial plasmodia developed into uninucleate sporoblasts by rosette-like fashion. Mature spores were pyriform and monokaryotic, measuring 4.02 ± 0.24 (3.63-4.53) µm long and 2.27 ± 0.15 (2.12-2.57) µm wide (N = 40). The polaroplast was bipartite with a tightly packed anterior lamellae and a loosely aligned posterior lamellae. Isofilar polar filament was coiled 9-11 turns and arranged in 2-3 rows. The phylogenetic analysis based on the obtained SSU rDNA sequence indicated that the N. dubium n. sp. clustered with the freshwater oligochaete-infecting N. aurantiae to form an independent monophyletic group, positioned at the base of Clade 4. In addition, we analyzed the genetic diversity in three N. dubium n. sp. isolates based on the rDNA (SSU rDNA, ITS and LSU rDNA) and Rpb1 gene. The genetic variation among the rDNA sequences was not distinct, however, high nucleotide diversity could be observed in Rpb1 gene, and a wide variety of Rpb1 haplotypes were identified within each isolate. Genetic recombination detected in the Rpb1 sequences presumes cryptic sexual process occurring in N. dubium n. sp. Statistical evolutionary analyses further indicated that the purifying selection eliminated mutations in the Rpb1 gene.

Sar1 Interacts with Sec23/Sec24 and Sec13/Sec31 Complexes: Insight into Its Involvement in the Assembly of Coat Protein Complex II in the Microsporidian Nosema bombycis

Microsporidia, as unicellular eukaryotes, also have an endomembrane system for transporting proteins, which is essentially similar to those of other eukaryotes. In eukaryotes, coat protein complex II (COPII) consists of Sar1, Sec23, Sec24, Sec13, and Sec31 and mediates protein transport from the endoplasmic reticulum (ER) to the Golgi apparatus. Sar1 is the central player in the regulation of coat protein complex II vesicle formation in the endoplasmic reticulum. In this study, we successfully cloned the NbSar1, NbSec23-1, NbSec23-2, NbSec24-1, NbSec24-2, NbSec13, NbSec31-1, and NbSec31-2 genes and prepared NbSar1 polyclonal antibody. We found that NbSar1 was localized mainly in the perinuclear cytoplasm of Nosema bombycis by immunofluorescence analysis (IFA). Yeast two-hybrid assays demonstrated that NbSar1 interacts with NbSec23-2, NbSec23-2 interacts with NbSec24-1 or NbSec24-2, NbSec23-1 interacts with NbSec31, and NbSec31 interacts with NbSec13. Moreover, the silencing of NbSar1 by RNA interference resulted in the aberrant expression of NbSar1, NbSec23-1, NbSec24-1, NbSec24-2, NbSec13, NbSec31-1, and NbSec31-2 and significantly inhibited the proliferation of N. bombycis. Altogether, these findings indicated that the subunits of coat protein complex II work together to perform functions in the proliferation of N. bombycis and that NbSar1 may play a crucial role in coat protein complex II vesicle formation. IMPORTANCE As eukaryotes, microsporidia have retained the endomembrane system for transporting and sorting proteins throughout their evolution. Whether the microsporidia form coat protein complex II (COPII) vesicles to transport cargo proteins and whether they play other roles besides cargo transport are not fully explained at present. Our results showed that NbSar1, NbSec23-1/NbSec23-2, NbSec24-1/NbSec24-2, NbSec13, and NbSec31 might be assembled to form COPII in the ER of N. bombycis, and the functions of COPII are also closely related to the proliferation of N. bombycis, this may be a new target for the prevention of pébrine disease of the silkworm.

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

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

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

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

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.

Honey Bee Habitat Sharing Enhances Gene Flow of the Parasite Nosema ceranae

Host-parasite co-evolution is a process of reciprocal, adaptive genetic change. In natural conditions, parasites can shift to other host species, given both host and parasite genotypes allow this. Even though host-parasite co-evolution has been extensively studied both theoretically and empirically, few studies have focused on parasite gene flow between native and novel hosts. Nosema ceranae is a native parasite of the Asian honey bee Apis cerana, which infects epithelial cells of mid-guts. This parasite successfully switched to the European honey bee Apis mellifera, where high virulence has been reported. In this study, we used the parasite N. ceranae and both honey bee species as model organisms to study the impacts of two-host habitat sharing on parasite diversity and virulence. SNVs (Single Nucleotide Variants) were identified from parasites isolated from native and novel hosts from sympatric populations, as well as novel hosts from a parapatric population. Parasites isolated from native hosts showed the highest levels of polymorphism. By comparing the parasites isolated from novel hosts between sympatric and parapatric populations, habitat sharing with the native host significantly enhanced parasite diversity, suggesting there is continuing gene flow of parasites between the two host species in sympatric populations.

The Role of Nosema ceranae (Microsporidia: Nosematidae) in Honey Bee Colony Losses and Current Insights on Treatment

Honeybee populations have locally and temporally declined in the last few years because of both biotic and abiotic factors. Among the latter, one of the most important reasons is infection by the microsporidia Nosema ceranae, which is the etiological agent of type C nosemosis. This species was first described in Asian honeybees (Apis cerana). Nowadays, domestic honeybees (Apis mellifera) worldwide are also becoming infected due to globalization. Type C nosemosis can be asymptomatic or can cause important damage to bees, such as changes in temporal polyethism, energy and oxidative stress, immunity loss, and decreased average life expectancy. It causes drastic reductions in workers, numbers of broods, and honey production, finally leading to colony loss. Common treatment is based on fumagillin, an antibiotic with side effects and relatively poor efficiency, which is banned in the European Union. Natural products, probiotics, food supplements, nutraceuticals, and other veterinary drugs are currently under study and might represent alternative treatments. Prophylaxis and management of affected colonies are essential to control the disease. While N. ceranae is one potential cause of bee losses in a colony, other factors must also be considered, especially synergies between microsporidia and the use of insecticides.

Polymorphism of 16s rRNA Gene: Any Effect on the Biomolecular Quantitation of the Honey Bee (Apis mellifera L., 1758) Pathogen Nosema ceranae?

The microsporidian Nosema ceranae is a severe threat to the western honey bee Apis mellifera, as it is responsible for nosemosis type C, which leads the colonies to dwindle and collapse. Infection quantification is essential to clinical and research aims. Assessment is made often with molecular assays based on rRNA genes, which are present in the N. ceranae genome as multiple and polymorphic copies. This study aims to compare two different methods of Real-Time PCR (qPCR), respectively relying on the 16S rRNA and Hsp70 genes, the first of which is described as a multiple and polymorphic gene. Young worker bees, hatched in the laboratory and artificially inoculated with N. ceranae spores, were incubated at 33 °C and subject to different treatment regimens. Samples were taken post-infection and analyzed with both qPCR methods. Compared to Hsp70, the 16S rRNA method systematically detected higher abundance. Straightforward conversion between the two methods is made impossible by erratic 16s rRNA/Hsp70 ratios. The 16s rRNA polymorphism showed an increase around the inoculated dose, where a higher prevalence of ungerminated spores was expected due to the treatment effects. The possible genetic background of that irregular distribution is discussed in detail. The polymorphic nature of 16S rRNA showed to be a limit in the infection quantification. More reliably, the N. ceranae abundance can be assessed in honey bee samples with methods based on the single-copy gene Hsp70.

Genome and Evolutionary Analysis of Nosema ceranae: A Microsporidian Parasite of Honey Bees

Microsporidia comprise a phylum of single cell, intracellular parasites and represent the earliest diverging branch in the fungal kingdom. The microsporidian parasite Nosema ceranae primarily infects honey bee gut epithelial cells, leading to impaired memory, suppressed host immune responses and colony collapse under certain circumstances. As the genome of N. ceranae is challenging to assembly due to very high genetic diversity and repetitive region, the genome was re-sequenced using long reads. We present a robust 8.8 Mbp genome assembly of 2,280 protein coding genes, including a high number of genes involved in transporting nutrients and energy, as well as drug resistance when compared with sister species Nosema apis. We also describe the loss of the critical protein Dicer in approximately half of the microsporidian species, giving new insights into the availability of RNA interference pathway in this group. Our results provided new insights into the pathogenesis of N. ceranae and a blueprint for treatment strategies that target this parasite without harming honey bees. The unique infectious apparatus polar filament and transportation pathway members can help to identify treatments to control this parasite.

Identification of pathogens in the invasive hornet Vespa velutina and in native Hymenoptera (Apidae, Vespidae) from SW-Europe

Invasive species contribute to deteriorate the health of ecosystems due to their direct effects on native fauna and the local parasite-host dynamics. We studied the potential impact of the invasive hornet Vespa velutina on the European parasite-host system by comparing the patterns of diversity and abundance of pathogens (i.e. Microsporidia: Nosematidae; Euglenozoa: Trypanosomatidae and Apicomplexa: Lipotrophidae) in European V. velutina specimens with those in the native European hornet Vespa crabro, as well as other common Hymenoptera (genera Vespula, Polistes and Bombus). We show that (i) V. velutina harbours most common hymenopteran enteropathogens as well as several new parasitic taxa. (ii) Parasite diversity in V. velutina is most similar to that of V. crabro. (iii) No unambiguous evidence of pathogen release by V. velutina was detected. This evidence together with the extraordinary population densities that V. velutina reaches in Europe (around of 100,000 individuals per km² per year), mean that this invasive species could severely alter the native pathogen-host dynamics either by actively contributing to the dispersal of the parasites and/or by directly interacting with them, which could have unexpected long-term harmful consequences on the native entomofauna.

Bioactivity studies of porphyrinoids against microsporidia isolated from honeybees

Microsporidian infections are dangerous to honeybees due to the absence of an efficient treatment for nosemosis. In the present work, the abilities of several porphyrins to directly inactivate microsporidia derived from Nosema-infected honeybees were studied in vitro. Amide derivatives of protoporphyrin IX (PPIX) conjugated with one and two amino acid moieties were synthesized, and their activities were compared with those of two cationic porphyrins, TMePyP and TTMePP. The most active porphyrins, PP[Lys-Asp]₂, PP[Lys-TFA]₂, PP[Asp(ONa)₂]₂ and PP[Lys-Lys]₂ at concentrations as low as 10–50 µM exerted significant effects on microsporidia, reducing the number of spores by 67–80% compared to the control. Live-cell imaging of the spores treated with porphyrins showed that only 1.6% and 3.0% of spores remained alive after 24 h-incubation with 50 µM PP[Asp(ONa)₂]₂ and PP[Lys-Asp]₂, respectively. The length of the amino acid side chains and their identity in the PPIX molecules affected the bioactivity of the porphyrin. Importantly, the irradiation of the porphyrins did not enhance their potency in destroying Nosema spores. We showed that the porphyrins accumulated inside the living spores but not inside dead spores, thus the destruction of the microsporidia by non-metallated porphyrins is not dependent on photosensitization, but is associated with their active transport into the spore cell. When administered to honeybees in vivo, PPIX[Lys-TFA]₂ and PPIX[Lys-Lys]₂ reduced spore loads by 69–76% in infected individuals. They both had no toxic effect on honeybees, in contrast to zinc-coordinated porphyrin.

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.

Genetic and Genome Analyses Reveal Genetically Distinct Populations of the Bee Pathogen Nosema ceranae from Thailand

The recent global decline in Western honeybee (Apis mellifera) populations is of great concern for pollination and honey production worldwide. Declining honeybee populations are frequently infected by the microsporidian pathogen Nosema ceranae. This species was originally described in the Asiatic honeybee (Apis cerana), and its identification in global A. mellifera hives could result from a recent host transfer. Recent genome studies have found that global populations of this parasite are polyploid and that humans may have fueled their global expansion. To better understand N. ceranae biology, we investigated its genetic diversity within part of their native range (Thailand) and among different hosts (A. mellifera, A. cerana) using both PCR and genome-based methods. We find that Thai N. ceranae populations share many SNPs with other global populations and appear to be clonal. However, in stark contrast with previous studies, we found that these populations also carry many SNPs not found elsewhere, indicating that these populations have evolved in their current geographic location for some time. Our genome analyses also indicate the potential presence of diploidy within Thai populations of N. ceranae.

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.

Genotype diversity in the honey bee parasite Nosema ceranae: multi-strain isolates, cryptic sex or both?

Background

There is great controversy as to whether Microsporidia undergo a sexual cycle. In the paradigmatic case of Nosema ceranae, although there is no morphological evidence of sex, some meiosis-specific genes are present in its reduced genome and there is also high intraspecific variability, with incongruent phylogenies having been systematically obtained. The possibility of sexual recombination is important from an epidemiological standpoint, particularly as N. ceranae is considered to be a major factor in the current disquieting epidemic of widespread bee colony losses. This parasite apparently originated in oriental honey bees, spreading out of Asia and Australia to infect honey bees worldwide. This study had three main objectives: i) to obtain genetic markers that are not part of known multi-copy arrays for strain determination; ii) to shed light on the intraspecific variability and recombination of N. ceranae; and iii) to assess the variability in N. ceranae populations. The answers to these questions are critical to understand the capacity of adaptation of microsporidia.

Results

Biallelic polymorphisms were detected at a number of specific points in the five coding loci analyzed from European and Australian isolates of N. ceranae. Heterozygous genotypes were abundant and cloning experiments demonstrate that they reflect the existence of multiple alternative sequences in each isolate. The comparisons of different clones and genotypes clearly indicate that new haplotypes are generated by homologous recombination.

Conclusions

The N. ceranae isolates from honey bees correspond to genotypically distinct populations, revealing that individual honey bees may not be infected by a particular clone but rather, a pool of different strains. Homologous recombination implies the existence of a cryptic sex cycle yet to be described in N. ceranae. There are no diagnostic alleles associated with Australian or European origins, nor are there differences between the two hosts, A. cerana and A. mellifera, supporting the absence of biological barriers for N. ceranae transmission. Diversity is high among microsporidia of both these origins, and the maintenance of a high heterozygosis in the recently invaded European populations, could hypothetically underlie the stronger virulence of N. ceranae observed in A. mellifera.

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

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

REPORT

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

Genetic diversity of two Daphnia-infecting microsporidian parasites, based on sequence variation in the internal transcribed spacer region

Background

Microsporidia are spore-forming obligate intracellular parasites that include both emerging pathogens and economically important disease agents. However, little is known about the genetic diversity of microsporidia. Here, we investigated patterns of geographic population structure, intraspecific genetic variation, and recombination in two microsporidian taxa that commonly infect cladocerans of the Daphnia longispina complex in central Europe. Taken together, this information helps elucidate the reproductive mode and life-cycles of these parasite species.

Methods

Microsporidia-infected Daphnia were sampled from seven drinking water reservoirs in the Czech Republic. Two microsporidia species (Berwaldia schaefernai and microsporidium lineage MIC1) were sequenced at the internal transcribed spacer (ITS) region, using the 454 pyrosequencing platform. Geographical structure analyses were performed applying Fisher’s exact tests, analyses of molecular variance, and permutational MANOVA. To evaluate the genetic diversity of the ITS region, the number of polymorphic sites and Tajima’s and Watterson’s estimators of theta were calculated. Tajima’s D was also used to determine if the ITS in these taxa evolved neutrally. Finally, neighbour similarity score and pairwise homology index tests were performed to detect recombination events.

Results

While there was little variation among Berwaldia parasite strains infecting different host populations, the among-population genetic variation of MIC1 was significant. Likewise, ITS genetic diversity was lower in Berwaldia than in MIC1. Recombination signals were detected only in Berwaldia.

Conclusion

Genetic tests showed that parasite populations could have expanded recently after a bottleneck or that the ITS could be under negative selection in both microsporidia species. Recombination analyses might indicate cryptic sex in Berwaldia and pure asexuality in MIC1. The differences observed between the two microsporidian species present an exciting opportunity to study the genetic basis of microsporidia-Daphnia coevolution in natural populations, and to better understand reproduction in these parasites.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-016-1584-4) contains supplementary material, which is available to authorized users.

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.

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

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

Population Genetics of Nosema apis and Nosema ceranae: One Host (Apis mellifera) and Two Different Histories

Two microsporidians are known to infect honey bees: Nosema apis and Nosema ceranae. Whereas population genetics data for the latter have been released in the last few years, such information is still missing for N. apis. Here we analyze the patterns of nucleotide polymorphism at three single-copy loci (PTP2, PTP3 and RPB1) in a collection of Apis mellifera isolates from all over the world, naturally infected either with N. apis (N = 22) or N. ceranae (N = 23), to provide new insights into the genetic diversity, demography and evolution of N. apis, as well as to compare them with evidence from N. ceranae. Neutral variation in N. apis and N. ceranae is of the order of 1%. This amount of diversity suggests that there is no substantial differentiation between the genetic content of the two nuclei present in these parasites, and evidence for genetic recombination provides a putative mechanism for the flow of genetic information between chromosomes. The analysis of the frequency spectrum of neutral variants reveals a significant surplus of low frequency variants, particularly in N. ceranae, and suggests that the populations of the two pathogens are not in mutation-drift equilibrium and that they have experienced a population expansion. Most of the variation in both species occurs within honey bee colonies (between 62%-90% of the total genetic variance), although in N. apis there is evidence for differentiation between parasites isolated from distinct A. mellifera lineages (20%-34% of the total variance), specifically between those collected from lineages A and C (or M). This scenario is consistent with a long-term host-parasite relationship and contrasts with the lack of differentiation observed among host-lineages in N. ceranae (< 4% of the variance), which suggests that the spread of this emergent pathogen throughout the A. mellifera worldwide population is a recent event.

The roles of sexual and asexual reproduction in the origin and dissemination of strains causing fungal infectious disease outbreaks

Sexual reproduction commonly refers to the reproductive process in which genomes from two sources are combined into a single cell through mating and then the zygote genomes are partitioned to progeny cells through meiosis. Reproduction in the absence of mating and meiosis is referred to as asexual or clonal reproduction. One major advantage of sexual reproduction is that it generates genetic variation among progeny which may allow for faster adaptation of the population to novel and/or stressful environments. However, adaptation to stressful or new environments can still occur through mutation, in the absence of sex. In this review, we analyzed the relative contributions of sexual and asexual reproduction in the origin and spread of strains causing fungal infectious diseases outbreaks. The necessity of sex and the ability of asexual fungi to initiate outbreaks are discussed. We propose a framework that relates the modes of reproduction to the origin and propagation of fungal disease outbreaks. Our analyses suggest that both sexual and asexual reproduction can play critical roles in the origin of outbreak strains and that the rapid spread of outbreak strains is often accomplished through asexual expansion.

Effects of genotype, environment, and their interactions on honey bee health in Europe

There are several reports of honey bee populations in Europe which survive without treatment for Varroa. However, when evaluated outside their native area, higher survival and resistance traits were not observed in colonies of a survivor population. Varroa infestation is strongly influenced by environmental factors, probably affecting threshold levels on a European scale. In a Europe-wide experiment colonies of local origin survived significantly longer than colonies of non-local origin, clearly indicating the presence of genotype-environment interactions. Transmission by Varroa selects for virulent strains of DWV, but it is currently unknown how these may interact with different genotypes of bees. The distribution of Nosema ceranae is significantly affected by environment, but there is at least one Nosema-resistant population.

Genome analyses suggest the presence of polyploidy and recent human-driven expansions in eight global populations of the honeybee pathogen Nosema ceranae

Nosema ceranae is a microsporidian pathogen whose infections have been associated with recent global declines in the populations of western honeybees (Apis mellifera). Despite the outstanding economic and ecological threat that N. ceranae may represent for honeybees worldwide, many aspects of its biology, including its mode of reproduction, propagation and ploidy, are either very unclear or unknown. In the present study, we set to gain knowledge in these biological aspects by re-sequencing the genome of eight isolates (i.e. a population of spores isolated from one single beehive) of this species harvested from eight geographically distant beehives, and by investigating their level of polymorphism. Consistent with previous analyses performed using single gene sequences, our analyses uncovered the presence of very high genetic diversity within each isolate, but also very little hive-specific polymorphism. Surprisingly, the nature, location and distribution of this genetic variation suggest that beehives around the globe are infected by a population of N. ceranae cells that may be polyploid (4n or more), and possibly clonal. Lastly, phylogenetic analyses based on genome-wide single-nucleotide polymorphism data extracted from these parasites and mitochondrial sequences from their hosts all failed to support the current geographical structure of our isolates.