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

Current status of Nosema spp. infection cases in apis mellifera in eurasian countries and Ptp3 gene haplotypes in the Republic of Tatarstan, Russia

The current status of Nosema spp. infections in A. mellifera throughout Eurasia was characterized using electronic databases. Although N. ceranae was predominantly detected in southwestern and south-central regions and N. apis in northwestern and north-central areas, most studies reported the occurrence of both species in Eurasia. In addition, the occurrence of Nosema spp. and Ptp3 gene haplotypes was investigated in the Republic of Tatarstan, Russia. Most of the examined honey bees were infected with both N. apis and N. ceranae. N. apis and N. ceranae isolates were either heterozygous or belonged to different strains and showed infection with more than one strain. New haplotypes were found for N. apis and N. ceranae in the Republic of Tatarstan, Russia. This study expands the data regarding existing haplotypes of Nosema species: there are currently 9 shared and 56 unique Ptp3 nucleotide sequence haplotypes of N. ceranae, and 2 shared and 7 unique haplotypes of N. apis, respectively.

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.

First report of Ovipleistophora ovariae and O. diplostomuri in China provides new insights into the intraspecific genetic variation and extends their distribution

Microsporidia of the genus Ovipleistophora are generally parasites of fishes and aquatic crustaceans. In the current study, Ovipleistophora diplostomuri and O. ovariae were firstly reported from Culter alburnus and Xenocypris argentea and Parabramis pekinensis, respectively. Both of them exclusively infected fish ovary and were morphologically, ultrastructurally and genetically characterized. Sporogony occurred in direct contact with the host cell cytoplasm and sporophorous vesicles were not observed for the new isolates of these two Ovipleistophora species. Spores of O. ovariae were for the first time observed to be dimorphic. Genetic analysis indicated that the genetic variation in the ITS and LSU sequences was distinct among between-host O. diplostomuri isolates. High sequence variation in ITS sequence suggests that it can be a reliable molecular marker to explore the population genetics of O. diplostomuri. This is the first report of these two Ovipleistophora species in China which extends their host and geographical range.

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.

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.

Targeting the honey bee gut parasite Nosema ceranae with siRNA positively affects gut bacteria

BACKGROUND

Gut microbial communities can contribute positively and negatively to host health. So far, eight core bacterial taxonomic clusters have been reported in honey bees. These bacteria are involved in host metabolism and defenses. Nosema ceranae is a gut intracellular parasite of honey bees which destroys epithelial cells and gut tissue integrity. Studies have shown protective impacts of honey bee gut microbiota towards N. ceranae infection. However, the impacts of N. ceranae on the relative abundance of honey bee gut microbiota remains unclear, and has been confounded during prior infection assays which resulted in the co-inoculation of bacteria during Nosema challenges. We used a novel method, the suppression of N. ceranae with specific siRNAs, to measure the impacts of Nosema on the gut microbiome.

RESULTS

Suppressing N. ceranae led to significant positive effects on microbial abundance. Nevertheless, 15 bacterial taxa, including three core taxa, were negatively correlated with N. ceranae levels. In particular, one co-regulated group of 7 bacteria was significantly negatively correlated with N. ceranae levels.

CONCLUSIONS

N. ceranae are negatively correlated with the abundance of 15 identified bacteria. Our results provide insights into interactions between gut microbes and N. ceranae during infection.

Brazilian fungal diversity represented by DNA markers generated over 20 years

Molecular techniques using fungal DNA barcoding (ITS) and other markers have been key to identifying the biodiversity of different geographic areas, mainly in megadiverse countries. Here, we provide an overview of the fungal diversity in Brazil based on DNA markers of phylogenetic importance generated since 1996. We retrieved fungal sequences of ITS, LSU, SSU, tef1-α, β-tubulin, rpb1, rpb2, actin, chitin synthase, and ATP6 from GenBank using different field keywords that indicated their origin in Brazil. A total of 19,440 sequences were recovered. ITS is the most representative marker (11,209 sequences), with 70.1% belonging to Ascomycota, 18.6% Basidiomycota, 10.2% unidentified, 1.1% Mucoromycota, two sequences of Olpidium bornovanus (Fungi incertae sedis), one sequence of Blastocladiomycota (Allomyces arbusculus), and one sequence of Chytridiomycota (Batrachochytrium dendrobatidis). Considering the sequences of all selected markers, only the phyla Cryptomycota and Entorrhizomycota were not represented. Based on ITS, using a cutoff of 98%, all sequences comprise 3047 OTUs, with the majority being Ascomycota (2088 OTUs) and Basidiomycota (681 OTUs). Previous numbers based mainly on morphological and bibliographical data revealed 5264 fungal species from Brazil, with a predominance of Basidiomycota (2741 spp.) and Ascomycota (1881 spp.). The unidentified ITS sequences not assigned to a higher taxonomic level represent 1.61% of all ITS sequences sampled and correspond to 38 unknown class-level lineages (75% cutoff). A maximum likelihood phylogeny based on LSU illustrates the fungal classes occurring in Brazil.

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.

Distribution and prevalence of Nosema apis and N. ceranae in temperate and subtropical eco-regions of Argentina

A total of 361 colonies from 59 apiaries located in two temperate and three subtropical eco-regions were examined during the post-harvest period to determine distribution and prevalence of Nosema spp. Apiaries from subtropical eco-regions showed a lower spore count than those from temperate eco-regions. Pure N. ceranae and co-infection were detected in apiaries from all regions. In contrast, pure N. apis infection was exclusively observed in the subtropical study region. The predominant detection of N. apis in a subtropical region joining a southern temperate region where mainly co-infected apiaries were identified is in contrast to previous reports.

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.

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.