Diversity and recombination of dispersed ribosomal DNA and protein coding genes in microsporidia

Widespread infection, diversification and old host associations of Nosema Microsporidia in European freshwater gammarids (Amphipoda)

The microsporidian genus Nosema is primarily known to infect insects of economic importance stimulating high research interest, while other hosts remain understudied. Nosema granulosis is one of the formally described Nosema species infecting amphipod crustaceans, being known to infect only two host species. Our first aim was to characterize Nosema spp. infections in different amphipod species from various European localities using the small subunit ribosomal DNA (SSU) marker. Second, we aimed to assess the phylogenetic diversity, host specificity and to explore the evolutionary history that may explain the diversity of gammarid-infecting Nosema lineages by performing a phylogenetic reconstruction based on RNA polymerase II subunit B1 (RPB1) gene sequences. For the host species Gammarus balcanicus, we also analyzed whether parasites were in excess in females to test for sex ratio distortion in relation with Nosema infection. We identified Nosema spp. in 316 individuals from nine amphipod species being widespread in Europe. The RPB1-based phylogenetic reconstruction using newly reported sequences and available data from other invertebrates identified 39 haplogroups being associated with amphipods. These haplogroups clustered into five clades (A-E) that did not form a single amphipod-infecting monophyletic group. Closely related sister clades C and D correspond to Nosema granulosis. Clades A, B and E might represent unknown Nosema species infecting amphipods. Host specificity seemed to be variable with some clades being restricted to single hosts, and some that could be found in several host species. We show that Nosema parasite richness in gammarid hosts is much higher than expected, illustrating the advantage of the use of RPB1 marker over SSU. Finally, we found no hint of sex ratio distortion in Nosema clade A infecting G. balcanicus. This study shows that Nosema spp. are abundant, widespread and diverse in European gammarids. Thus, Nosema is as diverse in aquatic as in terrestrial hosts.

Possible seasonal and diurnal modulation of Gammarus pulex (Crustacea, Amphipoda) drift by microsporidian parasites

In lotic freshwater ecosystems, the drift or downstream movement of animals (e.g., macroinvertebrates) constitutes a key dispersal pathway, thus shaping ecological and evolutionary patterns. There is evidence that macroinvertebrate drift may be modulated by parasites. However, most studies on parasite modulation of host drifting behavior have focused on acanthocephalans, whereas other parasites, such as microsporidians, have been largely neglected. This study provides new insight into possible seasonal and diurnal modulation of amphipod (Crustacea: Gammaridae) drift by microsporidian parasites. Three 72 h drift experiments were deployed in a German lowland stream in October 2021, April, and July 2022. The prevalence and composition of ten microsporidian parasites in Gammarus pulex clade E varied seasonally, diurnally, and between drifting and stationary specimens of G. pulex. Prevalence was generally higher in drifting amphipods than in stationary ones, mainly due to differences in host size. However, for two parasites, the prevalence in drift samples was highest during daytime suggesting changes in host phototaxis likely related to the parasite's mode of transmission and site of infection. Alterations in drifting behavior may have important implications for G. pulex population dynamics and microsporidians' dispersal. The underlying mechanisms are more complex than previously thought.

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.

The emerging pathogen Enterocytozoon hepatopenaei drives a degenerative cyclic pattern in the hepatopancreas microbiome of the shrimp (Penaeus vannamei)

The microsporidian Enterocytozoon hepatopenaei (EHP) is an emerging pathogen that causes substantial economic losses in shrimp (Penaeus spp.) aquaculture worldwide. To prevent diseases in shrimp, the manipulation of the gut microbiota has been suggested. However, prior knowledge of the host-microbiome is necessary. We assessed the modulation of the microbiome (bacteria/fungi) and its predicted functions over the course of disease progression in shrimp experimentally challenged with EHP for 30 days using high throughput 16S rRNA and ITS amplicon sequencing. Infection grade was assessed for the first time by quantitative digital histopathology. According to the infection intensity, three disease-stages (early/developmental/late) were registered. During the early-stage, EHP was not consistently detected, and a high diversity of potentially beneficial microorganisms related to nutrient assimilation were found. In the development-stage, most of the shrimp start to register a high infection intensity related to a decrease in beneficial microorganisms and an increase in opportunistic/pathogenic fungi. During late-stage, animals displayed different infection intensities, showed a displacement of beneficial microorganisms by opportunistic/pathogenic bacteria and fungi related to pathogen infection processes and depletion of energetic reserves. The degenerative cyclic pattern of EHP infection and its effects on beneficial microorganisms and beneficial functions of the shrimp hepatopancreas microbiome are discussed.

Generalist parasites persist in degraded environments: a lesson learned from microsporidian diversity in amphipods

The present study provides new insight into suitable microsporidian–host associations. It relates regional and continental-wide host specialization in microsporidians infecting amphipods to degraded and recovering habitats across 2 German river catchments. It provides a unique opportunity to infer the persistence of parasites following anthropogenic disturbance and their establishment in restored rivers. Amphipods were collected in 31 sampling sites with differing degradation and restoration gradients. Specimens were morphologically (hosts) and molecularly identified (host and parasites). Amphipod diversity and abundance, microsporidian diversity, host phylogenetic specificity and continental-wide β-specificity were investigated and related to each other and/or environmental variables. Fourteen microsporidian molecular operational taxonomic units (MOTUs), mainly generalist parasites, infecting 6 amphipod MOTUs were detected, expanding the current knowledge on the host range by 17 interactions. There was no difference in microsporidian diversity and host specificity among restored and near-natural streams (Boye) or between those located in urban and rural areas (Kinzig). Similarly, microsporidian diversity was generally not influenced by water parameters. In the Boye catchment, host densities did not influence microsporidian MOTU richness across restored and near-natural sites. High host turnover across the geographical range suggests that neither environmental conditions nor host diversity plays a significant role in the establishment into restored areas. Host diversity and environmental parameters do not indicate the persistence and dispersal of phylogenetic host generalist microsporidians in environments that experienced anthropogenic disturbance. Instead, these might depend on more complex mechanisms such as the production of resistant spores, host switching and host dispersal acting individually or conjointly.

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.

Intragenomic sequence variations in the second internal transcribed spacer (ITS2) ribosomal DNA of the malaria vector Anopheles stephensi

Second Internal Transcribed Spacer (ITS2) ribosomal DNA (rDNA) sequence is a widely used molecular marker for species-identification or -delimitation due to observed concerted evolution which is believed to homogenize rDNA copies in an interbreeding population. However, intra-specific differences in ITS2 of Anopheles stephensi have been reported. This study reports the presence of intragenomic sequence variation in the ITS2-rDNA of An. stephensi and hypothesizes that observed intra-specific differences in this species may have resulted due to ambiguous DNA sequence-chromatogram resulting from intragenomic heterogeneity. Anopheles stephensi collected from different parts of India were sequenced for complete ITS2 and the variable region of 28S-rDNA (d1-d3 domains). Intragenomic variations were found in ITS2 region of all An. stephensi sequenced, but no such variation was observed in d1 to d3 domains of 28S-rDNA. Cloning and sequencing of ITS2 through the d3 domain of the 28S region of rDNA from representative samples from northern, central, and southern India confirmed the presence of intragenomic variation in ITS2 due to transitions at three loci and two bp indel in a di-nucleotide microsatellite locus. Multiple haplotypes were observed in ITS2 raised from such variations. Due to the absence of detectable intragenomic sequence variation in the d1 to d3 domain of 28S rDNA of An. stephensi, this region can serve as an ideal reference sequence for taxonomic and phylogenetic studies. The presence of intragenomic variation in rDNA should be carefully examined before using this as a molecular marker for species delimitation or phylogenetic analyses.

BEExact: a Metataxonomic Database Tool for High-Resolution Inference of Bee-Associated Microbial Communities

High-throughput 16S rRNA gene sequencing technologies have robust potential to improve our understanding of bee (Hymenoptera: Apoidea)-associated microbial communities and their impact on hive health and disease. Despite recent computation algorithms now permitting exact inferencing of high-resolution exact amplicon sequence variants (ASVs), the taxonomic classification of these ASVs remains a challenge due to inadequate reference databases. To address this, we assemble a comprehensive data set of all publicly available bee-associated 16S rRNA gene sequences, systematically annotate poorly resolved identities via inclusion of 618 placeholder labels for uncultivated microbial dark matter, and correct for phylogenetic inconsistencies using a complementary set of distance-based and maximum likelihood correction strategies. To benchmark the resultant database (BEExact), we compare performance against all existing reference databases in silico using a variety of classifier algorithms to produce probabilistic confidence scores. We also validate realistic classification rates on an independent set of ∼234 million short-read sequences derived from 32 studies encompassing 50 different bee types (36 eusocial and 14 solitary). Species-level classification rates on short-read ASVs range from 80 to 90% using BEExact (with ∼20% due to "bxid" placeholder names), whereas only ∼30% at best can be resolved with current universal databases. A series of data-driven recommendations are developed for future studies. We conclude that BEExact (https://github.com/bdaisley/BEExact) enables accurate and standardized microbiota profiling across a broad range of bee species-two factors of key importance to reproducibility and meaningful knowledge exchange within the scientific community that together, can enhance the overall utility and ecological relevance of routine 16S rRNA gene-based sequencing endeavors.IMPORTANCE The failure of current universal taxonomic databases to support the rapidly expanding field of bee microbiota research has led to many investigators relying on "in-house" reference sets or manual classification of sequence reads (usually based on BLAST searches), often with vague identity thresholds and subjective taxonomy choices. This time-consuming, error- and bias-prone process lacks standardization, cripples the potential for comparative cross-study analysis, and in many cases is likely to incorrectly sway study conclusions. BEExact is structured on and leverages several complementary bioinformatic techniques to enable refined inference of bee host-associated microbial communities without any other methodological modifications necessary. It also bridges the gap between current practical outcomes (i.e., phylotype-to-genus level constraints with 97% operational taxonomic units [OTUs]) and the theoretical resolution (i.e., species-to-strain level classification with 100% ASVs) attainable in future microbiota investigations. Other niche habitats could also likely benefit from customized database curation via implementation of the novel approaches introduced in this study.

Evolutionary Diversity in the Intracellular Microsporidian Parasite Nosema sp. Infecting Wild Silkworm Revealed by IGS Nucleotide Sequence Diversity

Intracellular microsporidian Nosema mylitta infects Indian wild silkworm Antheraea mylitta causing pebrine disease. Genetic structure and phylogeny of N. mylitta are analysed using nucleotide variability in 5S ribosomal DNA and intergenic spacer (IGS) sequence from 20 isolates collected from Southern, Northern and Central regions of Jharkhand State. Nucleotide diversity (π) and genetic differentiation Gst were highest in the Central isolates whereas lowest in the North. Among the isolates, absence of nucleotides, transitions and transversions were observed. Haplotyping showed nucleotide variability at 83 positions in IGS and 13 positions in 5S rDNA. Haplotype-based genetic differentiation was 0.96 to 0.97 whereas nucleotide sequence-based genetic differentiation was higher (Ks = 22.29) between Southern and Central isolates. Bottleneck analysis showed negative value for Tajima's D and other summary statistics revealing induction of loss of rare alleles and population explosion. From IGS, 17 ancestral sequences were inferred by Network algorithm. Core of nine closely related nodes having ancient nucleotides and peripheral nodes with highly divergent nucleotides were derived. Most diverged peripheral haplotype was Bero (H11) from the Central region whereas Deoghar (H3) of the Northern region diverged early. Phylogeny of N. mylitta grouped Southern and Northern isolates together revealed weak phylogenetic signal for these locations. Phylogeny of N. mylitta with Nosema sp. infecting other lepidopterans clustered N. mylitta isolates with N. antheraea and N. philosamiae of China indicating genetic similarity whereas other species were dissimilar showing diversity irrespective of country of origin.

Widespread prevalence but contrasting patterns of intragenomic rRNA polymorphisms in nematodes: Implications for phylogeny, species delimitation and life history inference

Ribosomal RNA genes have long been a favoured locus in phylogenetic and metabarcoding studies. Within a genome, rRNA loci are organized as tandem repeated arrays and the copies are homogenized through the process of concerted evolution. However, some level of rRNA variation (intragenomic polymorphism) is known to persist and be maintained in the genomes of many species. In nematode worms, the extent of rRNA polymorphism (RP) across species and the evolutionary and life history factors that contribute to the maintenance of intragenomic RP is largely unknown. Here, we present an extensive analysis across 30 terrestrial nematode species representing a range of free-living and parasitic taxa isolated worldwide. Our results indicate that RP is common and widespread, ribosome function appears to be maintained despite mutational changes, and intragenomic variants are stable in the genome and neutrally evolving. However, levels of variation were varied widely across rRNA locus and species, with some taxa observed to lack RP entirely. Higher levels of RP were significantly correlated with shorter generation time and high reproductive rates, and population-level factors may play a role in the geographic and phylogenetic structuring of rRNA variants observed in genera such as Rotylenchulus and Pratylenchus. Although RP did not dramatically impact the clustering and recovery of taxa in mock metabarcoding analyses, the present study has significant implications for global biodiversity estimates of nematode species derived from environmental rRNA amplicon studies, as well as our understanding of the evolutionary and ecological factors shaping genetic diversity across the nematode Tree of Life.

Microsporidian infections in the species complex Gammarus roeselii (Amphipoda) over its geographical range: evidence for both host-parasite co-diversification and recent host shifts

BACKGROUND

Microsporidians are obligate endoparasites infecting taxonomically diverse hosts. Both vertical (from mother to eggs) and horizontal (between conspecifics or between species) transmission routes are known. While the former may promote co-speciation and host-specificity, the latter may promote shifts between host species. Among aquatic arthropods, freshwater amphipod crustaceans are hosts for many microsporidian species. However, despite numerous studies, no general pattern emerged about host specificity and co-diversification. In south-eastern Europe, the gammarid Gammarus roeselii is composed of 13 cryptic lineages of Miocene to Pleistocene age but few genotypes of one lineage have spread postglacially throughout north-western Europe. Based on nearly 100 sampling sites covering its entire range, we aim to: (i) explore the microsporidian diversity present in G. roeselii and their phylogenetic relationships, especially in relation to the parasites infecting other Gammaridae; (ii) test if the host phylogeographical history might have impacted host-parasite association (e.g. co-diversifications or recent host shifts from local fauna).

METHODS

We used part of the small subunit rRNA gene as source of sequences to identify and determine the phylogenetic position of the microsporidian taxa infecting G. roeselii.

RESULTS

Microsporidian diversity was high in G. roeselii with 24 detected haplogroups, clustered into 18 species-level taxa. Ten microsporidian species were rare, infecting a few individual hosts in a few populations, and were mostly phylogenetically related to parasites from other amphipods or various crustaceans. Other microsporidians were represented by widespread genera with high prevalence: Nosema, Cucumispora and Dictyocoela. Two contrasting host association patterns could be observed. First, two vertically transmitted microsporidian species, Nosema granulosis and Dictyocoela roeselum, share the pattern of infecting G. roeselii over most of its range and are specific to this host suggesting the co-diversification scenario. This pattern contrasted with that of Dictyocoela muelleri, the three species of Cucumispora, and the rare parasites, present only in the recently colonised region by the host. These patterns suggest recent acquisitions from local host species, after the recent spread of G. roeselii.

CONCLUSIONS

Microsporidians infecting G. roeselii revealed two scenarios of host-parasite associations: (i) ancient associations with vertically transmitted parasites that probably co-diversified with their hosts, and (ii) host shifts from local host species, after the postglacial spread of G. roeselii.

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.

Development and comparison of real-time and conventional PCR tools targeting β-tubulin gene for detection of Nosema infection in silkworms

Microsporidiosis (Pebrine) caused by the microsporidian parasite is one of the important devastating disease which affect the silk production leading to an unprofitable harvest. Till date ribosomal RNA (rRNA) gene was used as a target for detection of microsporidian species. In this study, we describe conventional and SYBR green based real-time PCR techniques alternatively targeting β-tubulin gene for quantitative detection of microsporidia infecting both the mulberry and non-mulberry silkworms. The modified DNA extraction method followed in our study was found to be easy, economical and could be used for both conventional and real time PCR as template. The real time qPCR revealed the expression of β-tubulin gene in different infected tissues of the silkworm Bombyx mori. The sensitivity of the SYBR green based real time PCR was found to be 100 times more than the conventional PCR and PCR was found more sensitive than the microscopic examination. The developed method did not produce any false positive results with the other silkworm pathogens and healthy silkworm. The data suggest that both the developed PCR methods targeting β-tubulin gene could be used effectively in quarantine process at seed centres for early detection of microsporidian infection in silkworms.

Hexokinase as a versatile molecular genetic marker for Microsporidia

Hexokinase (HK) is a core glycolytic enzyme of Microsporidia which regulates host cell metabolic processes. The goal of the present study was to test for the utility of HK for molecular phylogenetics, species identification and molecular detection of microsporidia in infected insects. HK sequence-based reconstructions were essentially similar to those based upon largest subunit RNA polymerase (RPB1) gene sequences, as well as previously published rRNA gene and genome-based trees. Comparing HK sequences allowed clear differentiation of closely related taxa, such as Nosema bombycis and Nosema pyrausta. In Nosema ceranae, unique SNPs were found for an isolate from wild colonies of the Burzyan dark honey bee as compared with the isolates from domesticated European honey bee. Similarly, in Encephalitozoon cuniculi, HK was as effective as RPB1 for discrimination of isolates belonging to different ITS genotypes. Amplification using species-specific primers flanking short fragments at the 3′-end of HK gene showed the presence of infection in insect tissues infected with N. pyrausta, Nosema ceranae and Paranosema (Antonospora) locustae. For the latter parasite species, HK expression was also demonstrated at early stages of infection using total mRNA extracts of locust larvae. These results indicate the suitability of HK as a novel tool for molecular genetic studies of Microsporidia.

Vairimorpha ephestiae is a synonym of Vairimorpha necatrix (Opisthosporidia: Microsporidia) based on multilocus sequence analysis

An isolate of the microsporidium Vairimorpha ephestiae (originally isolated from Ephestia kühniella) from collection of Prof. J. Weiser was propagated in a laboratory culture of Galleria mellonella. Only disporoblastic sporogony was observed and formation of octospores, characteristic of the genus Vairimorpha, never occurred. A partial nucleotide sequence of the small subunit rRNA gene (1247 bp) for this microsporidium showed 100% identity to the homologous sequences of Vairimorpha (Nosema) necatrix (Genbank accession # U11051 and # DQ996241), a microsporidium with a broad host range within the Lepidoptera. Sequence similarity of protein-coding genes (RPB1, HSP70 and actin) between V. ephestiae and V. necatrix was about 98-100%. The level of genetic polymorphism in the RPB1 locus between these two species was essentially the same as between isolates of V. necatrix. It is therefore concluded that V. ephestiae is in fact an isolate of V. necatrix and the former species should be synonymized with the latter. Though described later, V. necatrix has prevailing usage and its precedence over V. ephestiae is proposed to conserve stability and avoid confusion.

Error-prone protein synthesis in parasites with the smallest eukaryotic genome

Microsporidia are parasitic fungi-like organisms that invade the interior of living cells and cause chronic disorders in a broad range of animals, including humans. These pathogens have the tiniest known genomes among eukaryotic species, for which they serve as a model for exploring the phenomenon of genome reduction in obligate intracellular parasites. Here we report a case study to show an apparent effect of overall genome reduction on the primary structure and activity of aminoacyl-tRNA synthetases, indispensable cellular proteins required for protein synthesis. We find that most microsporidian synthetases lack regulatory and eukaryote-specific appended domains and have a high degree of sequence variability in tRNA-binding and catalytic domains. In one synthetase, LeuRS, an apparent sequence degeneration annihilates the editing domain, a catalytic center responsible for the accurate selection of leucine for protein synthesis. Unlike accurate LeuRS synthetases from other eukaryotic species, microsporidian LeuRS is error-prone: apart from leucine, it occasionally uses its near-cognate substrates, such as norvaline, isoleucine, valine, and methionine. Mass spectrometry analysis of the microsporidium Vavraia culicis proteome reveals that nearly 6% of leucine residues are erroneously replaced by other amino acids. This remarkably high frequency of mistranslation is not limited to leucine codons and appears to be a general property of protein synthesis in microsporidian parasites. Taken together, our findings reveal that the microsporidian protein synthesis machinery is editing-deficient, and that the proteome of microsporidian parasites is more diverse than would be anticipated based on their genome sequences.

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.

Integrative Approach for the Reliable Detection and Specific Identification of the Microsporidium Loma morhua in Atlantic Cod (Gadus morhua)

Microsporidia are fungal parasites that infect diverse invertebrate and vertebrate hosts. Finfish aquaculture supports epizootics due to high host density and the high biotic potential of these parasites. Reliable methods for parasite detection and identification are a necessary precursor to empirical assessment of strategies to mitigate the effects of these pathogens during aquaculture. We developed an integrative approach to detect and identify Loma morhua infecting Atlantic cod. We show that the spleen is more reliable than the commonly presumed gills as best organ for parasite detection in spite of substantial morphological plasticity in xenoma complexes. We developed rDNA primers with 100% sensitivity in detecting L. morhua and with utility in distinguishing some congeneric Loma species. ITS sequencing is necessary to distinguish L. morhua from other congeneric microsporidia due to intraspecific nucleotide variation. 64% of L. morhua ITS variants from Atlantic cod have a 9-nucleotide motif that distinguishes it from Loma spp. infecting non-Gadus hosts. The remaining 36% of ITS variants from Atlantic cod are distinguished from currently represented Loma spp., particularly those infecting Gadus hosts, based on a 14-nucleotide motif. This research approach is amenable to developing templates in support of reliable detection and identification of other microsporidian parasites in fishes.

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.

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.

Invaders, natives and their enemies: distribution patterns of amphipods and their microsporidian parasites in the Ruhr Metropolis, Germany

Background

The amphipod and microsporidian diversity in freshwaters of a heterogeneous urban region in Germany was assessed. Indigenous and non-indigenous host species provide an ideal framework to test general hypotheses on potentially new host-parasite interactions, parasite spillback and spillover in recently invaded urban freshwater communities.

Methods

Amphipods were sampled in 17 smaller and larger streams belonging to catchments of the four major rivers in the Ruhr Metropolis (Emscher, Lippe, Ruhr, Rhine), including sites invaded and not invaded by non-indigenous amphipods. Species were identified morphologically (hosts only) and via DNA barcoding (hosts and parasites). Prevalence was obtained by newly designed parasite-specific PCR assays.

Results

Three indigenous and five non-indigenous amphipod species were detected. Gammarus pulex was further distinguished into three clades (C, D and E) and G. fossarum more precisely identified as type B. Ten microsporidian lineages were detected, including two new isolates (designated as Microsporidium sp. nov. RR1 and RR2). All microsporidians occurred in at least two different host clades or species. Seven genetically distinct microsporidians were present in non-invaded populations, six of those were also found in invaded assemblages. Only Cucumispora dikerogammari and Dictyocoela berillonum can be unambiguously considered as non-indigenous co-introduced parasites. Both were rare and were not observed in indigenous hosts. Overall, microsporidian prevalence ranged from 50 % (in G. roeselii and G. pulex C) to 73 % (G. fossarum) in indigenous and from 10 % (Dikerogammarus villosus) to 100 % (Echinogammarus trichiatus) in non-indigenous amphipods. The most common microsporidians belonged to the Dictyocoela duebenum- /D. muelleri- complex, found in both indigenous and non-indigenous hosts. Some haplotype clades were inclusive for a certain host lineage.

Conclusions

The Ruhr Metropolis harbours a high diversity of indigenous and non-indigenous amphipod and microsporidian species, and we found indications for an exchange of parasites between indigenous and non-indigenous hosts. No introduced microsporidians were found in indigenous hosts and prevalence of indigenous parasites in non-indigenous hosts was generally low. Therefore, no indication for parasite spillover or spillback was found. We conclude that non-indigenous microsporidians constitute only a minimal threat to the native amphipod fauna. However, this might change e.g. if C. dikerogammari adapts to indigenous amphipod species or if other hosts and parasites invade.

Electronic supplementary material

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

Microsporidia: Eukaryotic Intracellular Parasites Shaped by Gene Loss and Horizontal Gene Transfers

Microsporidia are eukaryotic parasites of many animals that appear to have adapted to an obligate intracellular lifestyle by modifying the morphology and content of their cells. Living inside other cells, they have lost many, or all, metabolic functions, resulting in genomes that are always gene poor and often very small. The minute content of microsporidian genomes led many to assume that these parasites are biochemically static and uninteresting. However, recent studies have demonstrated that these organisms can be surprisingly complex and dynamic. In this review I detail the most significant recent advances in microsporidian genomics and discuss how these have affected our understanding of many biological aspects of these peculiar eukaryotic intracellular pathogens.

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.

(Cryptic) sex in the microsporidian Nosema granulosis--evidence from parasite rDNA and host mitochondrial DNA

Microsporidia are single-celled, intracellular eukaryotes that parasitise a wide range of animals. The Nosema/Vairimorpha group includes some putative asexual species, and asexuality is proposed to have originated multiple times from sexual ancestors. Here, we studied the variation in the ribosomal DNA (rDNA) of 14 isolates of the presumed apomictic and vertically transmitted Nosema granulosis to evaluate its sexual status. The analysed DNA fragment contained a part of the small-subunit ribosomal gene (SSU) and the entire intergenic spacer (IGS). The mitochondrial cox1 gene of the host Gammarus duebeni (Crustacea) was analysed to temporally calibrate the system and to test the expectation of cophylogeny of host and parasite genealogies. Genetic variability of the SSU gene was very low within and between the isolates. In contrast, intraisolate (within a single host) variability of the IGS felt in two categories, because 12 isolates possess a very high IGS genetic diversity and two isolates were almost invariable in the IGS. This difference suggests variable models of rDNA evolution involving birth-and-death and unexpectedly concerted evolution. An alternative explanation could be a likewise unattended mixed infection of host individuals by more than one parasite strain. Despite considerable genetic divergence between associated host mitochondrial haplotypes, some N. granulosis 'IGS populations' seem not to belong to different gene pools; the relevant tests failed to show significant differences between populations. A set of recombinant IGS sequences made our data incompatible with the model of a solely maternally inherited, asexual species. In line with recent reports, our study supports the hypothesis that some assumed apomictic Microsporidia did not entirely abstain from the evolutionary advantages of sex. In addition, the presented data indicate that horizontal transmission may occur occasionally. This transmission mode could be a survival strategy of N. granulosis whose host often populates ephemeral habitats.

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.

New insights on the genetic diversity of the honeybee parasiteNosema ceranaebased onmultilocussequence analysis

The microsporidian parasiteNosema ceranaeis a common pathogen of the Western honeybee (Apis mellifera) whose variable virulence could be related to its genetic polymorphism and/or its polyphenism responding to environmental cues. Since the genotyping ofN. ceranaebased on unique marker sequences had been unsuccessful, we tested whether amultilocusapproach, assessing the diversity of ten genetic markers – encoding nine proteins and the small ribosomal RNA subunit – allowed the discrimination betweenN. ceranaevariants isolated from singleA. melliferaindividuals in four distant locations. High nucleotide diversity and allele content were observed for all genes. Most importantly, the diversity was mainly present within parasite populations isolated from single honeybee individuals. In contrast the absence of isolate differentiation precluded anytaxadiscrimination, even through amultilocusapproach, but suggested that similar populations of parasites seem to infect honeybees in distant locations. As statistical evolutionary analyses showed that the allele frequency is under selective pressure, we discuss the origin and consequences ofN. ceranaeheterozygosity in a single host and lack of population divergence in the context of the parasite natural and evolutionary history.

Identification of two new polar tube proteins related to polar tube protein 2 in the microsporidian Antonospora locustae

Microsporidia are obligate intracellular eukaryotic parasites with a broad host spectrum characterized by a unique and highly sophisticated invasion apparatus, the polar tube (PT). In a previous study, two PT proteins, named AlPTP1 (50 kDa) and AlPTP2 (35 kDa), were identified in Antonospora locustae, an orthoptera parasite that is used as a biological control agent against locusts. Antibodies raised against AlPTP2 cross-reacted with a band migrating at ~70 kDa, suggesting that this 70-kDa antigen is closely related to AlPTP2. A blastp search against the A. locustae genome database allowed the identification of two further PTP2-like proteins named AlPTP2b (568 aa) and AlPTP2c (599 aa). Both proteins are characterized by a specific serine- and glycine-rich N-terminal extension with elastomeric structural features and share a common C-terminal end conserved with AlPTP2 (~88% identity for the last 250 aa). MS analysis of the 70-kDa band revealed the presence of AlPTP2b. Specific anti-AlPTP2b antibodies labelled the extruded PTs of the A. locustae spores, confirming that this antigen is a PT component. Finally, we showed that several PTP2-like proteins are also present in other phylogenetically related insect microsporidia, including Anncaliia algerae and Paranosema grylli.

Microsporidia and 'the art of living together'

Parasitism, aptly defined as one of the 'living-together' strategies (Trager, 1986), presents a dynamic system in which the parasite and its host are under evolutionary pressure to evolve new and specific adaptations, thus enabling the coexistence of the two closely interacting partners. Microsporidia are very frequently encountered obligatory intracellular protistan parasites that can infect both animals and some protists and are a consummate example of various aspects of the 'living-together' strategy. Microsporidia, relatives of fungi in the superkingdom Opisthokonta, belong to the relatively small group of parasites for which the host cell cytoplasm is the site of both reproduction and maturation. The structural and physiological reduction of their vegetative stage, together with the manipulation of host cell physiology, enables microsporidia to live in the cytosolic environment for most of their life cycle in a way resembling endocytobionts. The ability to form structurally complex spores and the invention and assembly of a unique injection mechanism enable microsporidia to disperse within host tissues and between host organisms, resulting in long-lasting infections. Microsporidia have adapted their genomes to the intracellular way of life, evolved strategies how to obtain nutrients directly from the host and how to manipulate not only the infected cells, but also the hosts themselves. The enormous variability of host organisms and their tissues provide microsporidian parasites a virtually limitless terrain for diversification and ecological expansion. This review attempts to present a general overview of microsporidia, emphasising some less known and/or more recently discovered facets of their biology.