New insights on the genetic diversity of the honeybee parasite Nosema ceranae based on multilocus sequence analysis

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.

Intraspecific genetic diversity of the fish-infecting microsporidian parasite Pseudokabatana alburnus (Microsporidia)

Pseudokabatana alburnus is a xenoma-forming fish microsporidium, firstly described from the liver of the Culter alburnus from Poyang Lake in China. In the present study, P. alburnus was firstly reported from the ovary of 6 other East Asian minnows, including Squaliobarbus curriculus, Hemiculter leucisculus, Cultrichthys erythropterus, Pseudolaubuca engraulis, Toxabramis swinhonis, and Elopichthys bambusa. Genetic analysis revealed high sequence diversity in the ribosomal internal transcribed spacer region (ITS) and the largest subunit of RNA polymerase II (Rpb1) loci of P. alburnus isolated from different hosts and locations. The variation of Rpb1 mainly occurred in the 1,477–1737 bp regions. The presence of a wide variety of Rpb1 haplotypes within a single fish host, together with evidence of genetic recombination suggested that P. alburnus may have the intergenomic variation and sexual reproduction might be present in other hosts (possibly freshwater shrimp). Phylogenetic analysis and population genetic analysis showed that there was no geographical population divergence for P. alburnus. Homogeneity and high variability of ITS sequences indicates that ITS may be a suitable molecular marker to distinguish different P. alburnus isolates. Our data confirm the broad geographical distribution and host range of P. alburnus in the middle and lower reaches of the Yangtze River. Additionally, we emendated the genus Pseudokabatana to exclude the infection site, liver as one of the taxonomic criteria, and proposed that fish ovary was be the general infection site of P. alburnus.

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.

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.

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.

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.

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

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

Characterization of Nosema ceranae Genetic Variants from Different Geographic Origins

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

Geographical distribution and molecular detection of Nosema ceranae from indigenous honey bees of Saudi Arabia

The aim of the study was to detect the infection level of honey bees with Nosema apis and/or Nosema ceranae using microscopic and molecular analysis from indigenous honeybee race of eight Saudi Arabian geographical regions. A detailed survey was conducted and fifty apiaries were chosen at random from these locations. Infection level was determined both by microscope and Multiplex-PCR and data were analyzed using bioinformatics tools and phylogenetic analysis. Result showed that N. ceranae was the only species infecting indigenous honeybee colonies in Saudi Arabia. As determined by microscope, Nosema spores were found to be in 20.59% of total samples colonies, while 58% of the samples evaluated by PCR were found to be positive for N. ceranae, with the highest prevalence in Al-Bahah, a tropical wet and dry climatic region, whereas low prevalence was found in the regions with hot arid climate. Honeybees from all eight locations surveyed were positive for N. ceranae. This is the first report about the N. ceranae detection, contamination level and distribution pattern in Saudi Arabia.

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.

Clonal Evolution of Enterocytozoon bieneusi Populations in Swine and Genetic Differentiation in Subpopulations between Isolates from Swine and Humans

Enterocytozoon bieneusi is a widespread parasite with high genetic diversity among hosts. Its natural reservoir remains elusive and data on population structure are available only in isolates from primates. Here we describe a population genetic study of 101 E. bieneusi isolates from pigs using sequence analysis of the ribosomal internal transcribed spacer (ITS) and four mini- and microsatellite markers. The presence of strong linkage disequilibrium (LD) and limited genetic recombination indicated a clonal structure for the population. Bayesian inference of phylogeny, structural analysis, and principal coordinates analysis separated the overall population into three subpopulations (SP3 to SP5) with genetic segregation of the isolates at some geographic level. Comparative analysis showed the differentiation of SP3 to SP5 from the two known E. bieneusi subpopulations (SP1 and SP2) from primates. The placement of a human E. bieneusi isolate in pig subpopulation SP4 supported the zoonotic potential of some E. bieneusi isolates. Network analysis showed directed evolution of SP5 to SP3/SP4 and SP1 to SP2. The high LD and low number of inferred recombination events are consistent with the possibility of host adaptation in SP2, SP3, and SP4. In contrast, the reduced LD and high genetic diversity in SP1 and SP5 might be results of broad host range and adaptation to new host environment. The data provide evidence of the potential occurrence of host adaptation in some of E. bieneusi isolates that belong to the zoonotic ITS Group 1.

Single and multi-gene phylogeny of Hepatospora (Microsporidia) - a generalist pathogen of farmed and wild crustacean hosts

Almost half of all known microsporidian taxa infect aquatic animals. Of these, many cause disease in arthropods. Hepatospora, a recently erected genus, infects epithelial cells of the hepatopancreas of wild and farmed decapod crustaceans. We isolated Hepatospora spp. from three different crustacean hosts, inhabiting different habitats and niches; marine edible crab (Cancer pagurus), estuarine and freshwater Chinese mitten crab (Eriocheir sinensis) and the marine mussel symbiont pea crab (Pinnotheres pisum). Isolates were initially compared using histology and electron microscopy revealing variation in size, polar filament arrangement and nuclear development. However, sequence analysis of the partial SSU rDNA gene could not distinguish between the isolates (~99% similarity). In an attempt to resolve the relationship between Hepatospora isolated from E. sinensis and C. pagurus, six additional gene sequences were mined from on-going unpublished genome projects (RNA polymerase, arginyl tRNA synthetase, prolyl tRNA synthetase, chitin synthase, beta tubulin and heat shock protein 70). Primers were designed based on the above gene sequences to analyse Hepatospora isolated from pea crab. Despite application of gene sequences to concatenated phylogenies, we were unable to discriminate Hepatospora isolates obtained from these hosts and concluded that they likely represent a single species or, at least subspecies thereof. In this instance, concatenated phylogenetic analysis supported the SSU-based phylogeny, and further, demonstrated that microsporidian taxonomies based upon morphology alone are unreliable, even at the level of the species. Our data, together with description of H. eriocheir in Asian crab farms, reveal a preponderance for microvariants of this parasite to infect the gut of a wide array of decapods crustacean hosts and the potential for Hepatospora to exist as a cline across wide geographies and habitats.

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.

Environmental microbiology as a mosaic of explored ecosystems and issues

Microbes are phylogenetically (Archaea, Bacteria, Eukarya, and viruses) and functionally diverse. They colonize highly varied environments and rapidly respond to and evolve as a response to local and global environmental changes, including those induced by pollutants resulting from human activities. This review exemplifies the Microbial Ecology EC2CO consortium's efforts to explore the biology, ecology, diversity, and roles of microbes in aquatic and continental ecosystems.

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.

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

Nosema ceranae has been found infecting Apismellifera colonies with increasing frequency and it now represents a major threat to the health and long-term survival of these honeybees worldwide. However, so far little is known about the population genetics of this parasite. Here, we describe the patterns of genetic variation at three genomic loci in a collection of isolates from all over the world. Our main findings are: (i) the levels of genetic polymorphism (πS≈1%) do not vary significantly across its distribution range, (ii) there is substantial evidence for recombination among haplotypes, (iii) the best part of the observed genetic variance corresponds to differences within bee colonies (up to 88% of the total variance), (iv) parasites collected from Asian honeybees (Apis cerana and Apis florea) display significant differentiation from those obtained from Apismellifera (8-16% of the total variance, p<0.01) and (v) there is a significant excess of low frequency variants over neutral expectations among samples obtained from A. mellifera, but not from Asian honeybees. Overall these results are consistent with a recent colonization and rapid expansion of N. ceranae throughout A. mellifera colonies.

Exploiting the architecture and the features of the microsporidian genomes to investigate diversity and impact of these parasites on ecosystems

Fungal species play extremely important roles in ecosystems. Clustered at the base of the fungal kingdom are Microsporidia, a group of obligate intracellular eukaryotes infecting multiple animal lineages. Because of their large host spectrum and their implications in host population regulation, they influence food webs, and accordingly, ecosystem structure and function. Unfortunately, their ecological role is not well understood. Present also as highly resistant spores in the environment, their characterisation requires special attention. Different techniques based on direct isolation and/or molecular approaches can be considered to elucidate their role in the ecosystems, but integrating environmental and genomic data (for example, genome architecture, core genome, transcriptional and translational signals) is crucial to better understand the diversity and adaptive capacities of Microsporidia. Here, we review the current status of Microsporidia in trophic networks; the various genomics tools that could be used to ensure identification and evaluate diversity and abundance of these organisms; and how these tools could be used to explore the microsporidian life cycle in different environments. Our understanding of the evolution of these widespread parasites is currently impaired by limited sampling, and we have no doubt witnessed but a small subset of their diversity.

Evidence for weak genetic recombination at the PTP2 locus of Nosema ceranae

The microsporidian Nosema ceranae is an emergent pathogen that threatens the health of honeybees and other pollinators all over the world. Its recent rapid spread across a wide variety of host species and environments demonstrated an enhanced ability of adaptation, which seems to contradict the lack of evidence for genetic recombination and the absence of a sexual stage in its life cycle. Here we retrieved fresh data of the patterns of genetic variation at the PTP2 locus in naturally infected Apis mellifera colonies, by means of single genome amplification. This technique, designed to prevent the formation of chimeric haplotypes during polymerase chain reaction (PCR), provides more reliable estimates of the diversity levels and haplotype structure than standard PCR-cloning methods. Our results are consistent with low but significant rates of recombination in the history of the haplotypes detected: estimates of the population recombination rate are of the order of 30 and support recent evidence for unexpectedly high levels of variation of the parasites within honeybee colonies. These observations suggest the existence of a diploid stage at some point in the life cycle of this parasite and are relevant for our understanding of the dynamics of its expanding population.

Virulence and polar tube protein genetic diversity of Nosema ceranae (Microsporidia) field isolates from Northern and Southern Europe in honeybees (Apis mellifera iberiensis)

Infection of honeybees by the microsporidian Nosema ceranae is considered to be one of the factors underlying the increased colony losses and decreased honey production seen in recent years. However, these effects appear to differ in function of the climatic zone, the distinct beekeeping practices and the honeybee species employed. Here, we compared the response of Apis mellifera iberiensis worker bees to experimental infection with field isolates of N. ceranae from an Oceanic climate zone in Northern Europe (Netherlands) and from a Mediterranean region of Southern Europe (Spain). We found a notable but non-significant trend (P = 0.097) towards higher honeybee survival for bees infected with N. ceranae from the Netherlands, although no differences were found between the two isolates in terms of anatomopathological lesions in infected ventricular cells or the morphology of the mature and immature stages of the parasite. In addition, the population genetic survey of the N. ceranae PTP3 locus revealed high levels of genetic diversity within each isolate, evidence for meiotic recombination, and no signs of differentiation between the Dutch and Spanish populations. A cross-infection study is needed to further explore the differences in virulence observed between the two N. ceranae populations in field conditions.

High levels of genetic diversity inNosema ceranaewithinApis melliferacolonies

Nosema ceranaeis a widespread honeybee parasite, considered to be one of the pathogens involved in the colony losses phenomenon. To date, little is known about its intraspecific genetic variability. The few studies onN. ceranaevariation have focused on the subunits of ribosomal DNA, which are not ideal for this purpose and have limited resolution. Here we characterized three single copy loci (Actin, Hsp70andRPB1) in threeN. ceranaeisolates from Hungary and Hawaii. Our results provide evidence of unexpectedly high levels of intraspecific polymorphism, the coexistence of a wide variety of haplotypes within each bee colony, and the occurrence of genetic recombination inRPB1. Most haplotypes are not shared across isolates and derive from a few frequent haplotypes by a reduced number of singletons (mutations that appear usually just once in the sample), which suggest that they have a fairly recent origin. Overall, our data indicate that this pathogen has experienced a recent population expansion. The presence of multiple haplotypes within individual isolates could be explained by the existence of different strains ofN. ceranaeinfecting honeybee colonies in the field which complicates, and must not be overlooked, further analysis of host–parasite interactions.