Single-nucleotide polymorphisms of two closely related microsporidian parasites suggest a clonal population expansion after the last glaciation

Polyploidy is widespread in Microsporidia

Microsporidia are obligate intracellular eukaryotic parasites with an extremely broad host range. They have both economic and public health importance. Ploidy in microsporidia is variable, with a few species formally identified as diploid and one as polyploid. Given the increase in the number of studies sequencing microsporidian genomes, it is now possible to assess ploidy levels across all currently explored microsporidian diversity. We estimate ploidy for all microsporidian data sets available on the Sequence Read Archive using k-mer-based analyses, indicating that polyploidy is widespread in Microsporidia and that ploidy change is dynamic in the group. Using genome-wide heterozygosity estimates, we also show that polyploid microsporidian genomes are relatively homozygous, and we discuss the implications of these findings on the timing of polyploidization events and their origin.IMPORTANCEMicrosporidia are single-celled intracellular parasites, distantly related to fungi, that can infect a broad range of hosts, from humans all the way to protozoans. Exploiting the wealth of microsporidian genomic data available, we use k-mer-based analyses to assess ploidy status across the group. Understanding a genome's ploidy is crucial in order to assemble it effectively and may also be relevant for better understanding a parasite's behavior and life cycle. We show that tetraploidy is present in at least six species in Microsporidia and that these polyploidization events are likely to have occurred independently. We discuss why these findings may be paradoxical, given that Microsporidia, like other intracellular parasites, have extremely small, reduced genomes.

Near chromosome-level genome assembly of the microsporidium Hamiltosporidium tvaerminnensis

Microsporidia are intracellular parasitic fungi whose genomes rank among the smallest of all known eukaryotes. A number of outstanding questions remain concerning the evolution of their large-scale variation in genome architecture, responsible for genome size variation of more than an order of magnitude. This genome report presents the first near-chromosomal assembly of a large-genome microsporidium, Hamiltosporidium tvaerminnensis. Combined Oxford Nanopore, Pacific Biosciences (PacBio), and Illumina sequencing led to a genome assembly of 17 contigs, 11 of which represent complete chromosomes. Our assembly is 21.64 Mb in length, has an N50 of 1.44 Mb, and consists of 39.56% interspersed repeats. We introduce a novel approach in microsporidia, PacBio Iso-Seq, as part of a larger annotation pipeline for obtaining high-quality annotations of 3,573 protein-coding genes. Based on direct evidence from the full-length Iso-Seq transcripts, we present evidence for alternative polyadenylation and variation in splicing efficiency, which are potential regulation mechanisms for gene expression in microsporidia. The generated high-quality genome assembly is a necessary resource for comparative genomics that will help elucidate the evolution of genome architecture in response to intracellular parasitism.

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.

Using average nucleotide identity (ANI) to evaluate microsporidia species boundaries based on their genetic relatedness

Microsporidia are obligatory intracellular parasites related to fungi and since their discovery their classification and origin has been controversial due to their unique morphology. Early taxonomic studies of microsporidia were based on ultrastructural spore features, characteristics of their life cycle and transmission modes. However, taxonomy and phylogeny based solely on these characteristics can be misleading. SSU rRNA is a traditional marker used in taxonomical classifications, but the power of SSU rRNA to resolve phylogenetic relationships between microsporidia is considered weak at the species level, as it may not show enough variation to distinguish closely related species. Overall genome relatedness indices (OGRI), such as average nucleotide identity (ANI), allows fast and easy-to-implement comparative measurements between genomes to assess species boundaries in prokaryotes, with a 95% cutoff value for grouping genomes of the same species. Due to the increasing availability of complete genomes, metrics of genome relatedness have been applied for eukaryotic microbes taxonomy such as microsporidia. However, the distribution of ANI values and cutoff values for species delimitation have not yet been fully tested in microsporidia. In this study we examined the distribution of ANI values for 65 publicly available microsporidian genomes and tested whether the 95% cutoff value is a good estimation for circumscribing species based on their genetic relatedness.

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.

Microsporidia with Vertical Transmission Were Likely Shaped by Nonadaptive Processes

Microsporidia have the leanest genomes among eukaryotes, and their physiological and genomic simplicity has been attributed to their intracellular, obligate parasitic life-style. However, not all microsporidia genomes are small or lean, with the largest dwarfing the smallest ones by at least an order of magnitude. To better understand the evolutionary mechanisms behind this genomic diversification, we explore here two clades of microsporidia with distinct life histories, Ordospora and Hamiltosporidium, parasitizing the same host species, Daphnia magna. Based on seven newly assembled genomes, we show that mixed-mode transmission (the combination of horizontal and vertical transmission), which occurs in Hamiltosporidium, is found to be associated with larger and AT-biased genomes, more genes, and longer intergenic regions, as compared with the exclusively horizontally transmitted Ordospora. Furthermore, the Hamiltosporidium genome assemblies contain a variety of repetitive elements and long segmental duplications. We show that there is an excess of nonsynonymous substitutions in the microsporidia with mixed-mode transmission, which cannot be solely attributed to the lack of recombination, suggesting that bursts of genome size in these microsporidia result primarily from genetic drift. Overall, these findings suggest that the switch from a horizontal-only to a mixed mode of transmission likely produces population bottlenecks in Hamiltosporidium species, therefore reducing the effectiveness of natural selection, and allowing their genomic features to be largely shaped by nonadaptive processes.

Asexual parasites and their extraordinary host ranges

In diverse parasite taxa, from scale insects to root-knot nematodes, asexual lineages have exceptionally large host ranges, larger than those of their sexual relatives. Phylogenetic comparative studies of parasite taxa indicate that increases in host range and geographic range increase the probability of establishment of asexual lineages. At first pass, this convergence of traits appears counter-intuitive: intimate, antagonistic association with an enormous range of host taxa correlates with asexual reproduction, which should limit genetic variation within populations. Why would narrow host ranges favor sexual parasites and large host ranges favor asexual parasites? To take on this problem I link theory on ecological specialization to the two predominant hypotheses for the evolution of sex. I argue that both hypotheses predict a positive association between host range and the probability of invasion of asexual parasites, mediated either by variation in population size or in the strength of antagonistic coevolution. I also review hypotheses on colonization and the evolution of niche breadth in asexual lineages. I emphasize parasite taxa, with their diversity of reproductive modes and ecological strategies, as valuable assets in the hunt for solutions to the classic problems of the evolution of sex and geographic parthenogenesis.

Genetic resistance and specificity in sister taxa of Daphnia: insights from the range of host susceptibilities

BACKGROUND

Host genetic diversity can affect various aspects of host-parasite interactions, including individual-level effects on parasite infectivity, production of transmission stages and virulence, as well as population-level effects that reduce disease spread and prevalence, and buffer against widespread epidemics. However, a key aspect of this diversity, the genetic variation in host susceptibility, has often been neglected in interpreting empirical data and in theoretical studies. Daphnia similis naturally coexists with its competitor Daphnia magna and is more resistant to the endoparasitic microsporidium Hamiltosporidium tvaerminnensis, as suggested by a previous survey of waterbodies, which detected this parasite in D. magna, but not in D. similis. However, under laboratory conditions D. similis was sometimes found to be susceptible. We therefore asked if there is genetic variation for disease trait expression, and if the genetic variation in disease traits in D. similis is different from that of D. magna.

METHODS

We exposed ten clones of D. similis and ten clones of D. magna to three isolates of H. tvaerminnensis, and measured infection rates, parasite-induced host mortality and parasite spore production.

RESULTS

The two Daphnia species differ in the range and variation of their susceptibilities. The parasite produced on average two-fold more spores when growing in D. magna clones than in D. similis clones.

CONCLUSIONS

We confirm that D. similis is indeed much more resistant than D. magna and suggest that this could create a dilution effect in habitats where both species coexist.

Parasite-mediated selection in a natural metapopulation of Daphnia magna

Parasite-mediated selection varying across time and space in metapopulations is expected to result in host local adaptation and the maintenance of genetic diversity in disease-related traits. However, nonadaptive processes like migration and extinction-(re)colonization dynamics might interfere with adaptive evolution. Understanding how adaptive and nonadaptive processes interact to shape genetic variability in life-history and disease-related traits can provide important insights into their evolution in subdivided populations. Here we investigate signatures of spatially fluctuating, parasite-mediated selection in a natural metapopulation of Daphnia magna. Host genotypes from infected and uninfected populations were genotyped at microsatellite markers, and phenotyped for life-history and disease traits in common garden experiments. Combining phenotypic and genotypic data a Q_ST -F_ST -like analysis was conducted to test for signatures of parasite mediated selection. We observed high variation within and among populations for phenotypic traits, but neither an indication of host local adaptation nor a cost of resistance. Infected populations have a higher gene diversity (Hs) than uninfected populations and Hs is strongly positively correlated with fitness. These results suggest a strong parasite effect on reducing population level inbreeding. We discuss how stochastic processes related to frequent extinction-(re)colonization dynamics as well as host and parasite migration impede the evolution of resistance in the infected populations. We suggest that the genetic and phenotypic patterns of variation are a product of dynamic changes in the host gene pool caused by the interaction of colonization bottlenecks, inbreeding, immigration, hybrid vigor, rare host genotype advantage and parasitism. Our study highlights the effect of the parasite in ameliorating the negative fitness consequences caused by the high drift load in this metapopulation.

Multilocus Sequence Typing and Population Genetic Analysis of Enterocytozoon bieneusi: Host Specificity and Its Impacts on Public Health

Microsporidia comprise a large class of unicellular eukaryotic pathogens that are medically and agriculturally important, but poorly understood. There have been nearly 1,500 microsporidian species described thus far, which are variable in biology, genetics, genomics, and host specificity. Among those, Enterocytozoon bieneusi is the well-known species responsible for the most recorded cases of human microsporidian affections. The pathogen can colonize a broad range of mammals and birds and most of the animals surveyed share some genotypes with humans, posing a threat to public health. Based on DNA sequence analysis of the ribosomal internal transcribed spacer (ITS) and phylogenetic analysis, several hundreds of E. bieneusi genotypes have been defined and clustered into different genetic groups with varied levels of host specificity. However, single locus-based typing using ITS might have insufficient resolution to discriminate among E. bieneusi isolates with complex genetic or hereditary characteristics and to assess the elusive reproduction or transmission modes of the organism, highlighting the need for exploration and application of multilocus sequence typing (MLST) and population genetic tools. The present review begins with a primer on microsporidia and major microsporidian species, briefly introduces the recent advances on E. bieneusi ITS genotyping and phylogeny, summarizes recent MLST and population genetic data, analyzes the inter- and intragroup host specificity at the MLST level, and interprets the public health implications of host specificity in zoonotic or cross-species transmission of this ubiquitous fungus.

Novel insights on colonization routes and evolutionary potential of Colletotrichum kahawae, a severe pathogen of Coffea arabica

Pathogenic fungi are emerging at an increasing rate on a wide range of host plants, leading to tremendous threats to the global economy and food safety. Several plant pathogens have been considered to be invasive species, rendering large-scale population genomic analyses crucial to better understand their demographic history and evolutionary potential. Colletotrichum kahawae (Ck) is a highly aggressive and specialized pathogen, causing coffee berry disease in Arabica coffee in Africa. This pathogen leads to severe production losses and its dissemination out of Africa is greatly feared. To address this issue, a population genomic approach using thousands of single nucleotide polymorphisms (SNPs) spaced throughout the genome was used to unveil its demographic history and evolutionary potential. The current study confirms that Ck is a true clonal pathogen, perfectly adapted to green coffee berries, with three completely differentiated populations (Angolan, Cameroonian and East African). Two independent clonal lineages were found within the Angolan population as opposed to the remaining single clonal populations. The most probable colonization scenario suggests that this pathogen emerged in Angola and immediately dispersed to East Africa, where these two populations began to differentiate, followed by the introduction in Cameroon from an Angolan population. However, the differentiation between the two Angolan clonal lineages masks the mechanism for the emergence of the Cameroonian population. Our results suggest that Ck is completely differentiated from the ancestral lineage, has a low evolutionary potential and a low dispersion ability, with human transport the most likely scenario for its potential dispersion, which makes the fulfilment of the quarantine measures and management practices implemented crucial.

The genome of an intranuclear parasite, Paramicrosporidium saccamoebae, reveals alternative adaptations to obligate intracellular parasitism

Intracellular parasitism often results in gene loss, genome reduction, and dependence upon the host for cellular functioning. Rozellomycota is a clade comprising many such parasites and is related to the diverse, highly reduced, animal parasites, Microsporidia. We sequenced the nuclear and mitochondrial genomes of Paramicrosporidium saccamoebae [Rozellomycota], an intranuclear parasite of amoebae. A canonical fungal mitochondrial genome was recovered from P. saccamoebae that encodes genes necessary for the complete oxidative phosphorylation pathway including Complex I, differentiating it from most endoparasites including its sequenced relatives in Rozellomycota and Microsporidia. Comparative analysis revealed that P. saccamoebae shares more gene content with distantly related Fungi than with its closest relatives, suggesting that genome evolution in Rozellomycota and Microsporidia has been affected by repeated and independent gene losses, possibly as a result of variation in parasitic strategies (e.g. host and subcellular localization) or due to multiple transitions to parasitism.

QTL mapping of a natural genetic polymorphism for long-term parasite persistence in Daphnia populations

Knowing the determinants of the geographic ranges of parasites is important for understanding their evolutionary ecology, epidemiology and their potential to expand their range. Here we explore the determinants of geographic range in the peculiar case of a parasite species – the microsporidian Hamiltosporidium tvaerminnensis – that has a limited geographic distribution in a wide-spread host – Daphnia magna. We conducted a quantitative trait loci (QTLs) analysis with monoclonal F2D. magna populations originating from a cross between a susceptible northern European genotype and a resistant central European genotype. Contrary to our expectations, long-term persistence turned out to be a quantitative trait across the F2 offspring. Evidence for two QTLs, one epistatic interaction and for further minor QTL was found. This finding contrasts markedly with the previously described bimodal pattern for long-term parasite persistence in natural host genotypes across Europe and leaves open the question of how a quantitative genetic trait could determine the disjunct geographic distribution of the parasite across Europe.

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.

A Recent Whole-Genome Duplication Divides Populations of a Globally Distributed Microsporidian

The Microsporidia are a major group of intracellular fungi and important parasites of animals including insects, fish, and immunocompromised humans. Microsporidian genomes have undergone extreme reductive evolution but there are major differences in genome size and structure within the group: some are prokaryote-like in size and organisation (<3 Mb of gene-dense sequence) while others have more typically eukaryotic genome architectures. To gain fine-scale, population-level insight into the evolutionary dynamics of these tiny eukaryotic genomes, we performed the broadest microsporidian population genomic study to date, sequencing geographically isolated strains of Spraguea, a marine microsporidian infecting goosefish worldwide. Our analysis revealed that population structure across the Atlantic Ocean is associated with a conserved difference in ploidy, with American and Canadian isolates sharing an ancestral whole genome duplication that was followed by widespread pseudogenisation and sorting-out of paralogue pairs. While past analyses have suggested de novo gene formation of microsporidian-specific genes, we found evidence for the origin of new genes from noncoding sequence since the divergence of these populations. Some of these genes experience selective constraint, suggesting the evolution of new functions and local host adaptation. Combining our data with published microsporidian genomes, we show that nucleotide composition across the phylum is shaped by a mutational bias favoring A and T nucleotides, which is opposed by an evolutionary force favoring an increase in genomic GC content. This study reveals ongoing dramatic reorganization of genome structure and the evolution of new gene functions in modern microsporidians despite extensive genomic streamlining in their common ancestor.

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.

Transcriptomic profiling of host-parasite interactions in the microsporidian Trachipleistophora hominis

Background

Trachipleistophora hominis was isolated from an HIV/AIDS patient and is a member of a highly successful group of obligate intracellular parasites.

Methods

Here we have investigated the evolution of the parasite and the interplay between host and parasite gene expression using transcriptomics of T. hominis-infected rabbit kidney cells.

Results

T. hominis has about 30 % more genes than small-genome microsporidians. Highly expressed genes include those involved in growth, replication, defence against oxidative stress, and a large fraction of uncharacterised genes. Chaperones are also highly expressed and may buffer the deleterious effects of the large number of non-synonymous mutations observed in essential T. hominis genes. Host expression suggests a general cellular shutdown upon infection, but ATP, amino sugar and nucleotide sugar production appear enhanced, potentially providing the parasite with substrates it cannot make itself. Expression divergence of duplicated genes, including transporters used to acquire host metabolites, demonstrates ongoing functional diversification during microsporidian evolution. We identified overlapping transcription at more than 100 loci in the sparse T. hominis genome, demonstrating that this feature is not caused by genome compaction. The detection of additional transposons of insect origin strongly suggests that the natural host for T. hominis is an insect.

Conclusions

Our results reveal that the evolution of contemporary microsporidian genomes is highly dynamic and innovative. Moreover, highly expressed T. hominis genes of unknown function include a cohort that are shared among all microsporidians, indicating that some strongly conserved features of the biology of these enormously successful parasites remain uncharacterised.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1989-z) 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.

Interspecific competition in honeybee intracellular gut parasites is asymmetric and favours the spread of an emerging infectious disease

There is increasing appreciation that hosts in natural populations are subject to infection by multiple parasite species. Yet the epidemiological and ecological processes determining the outcome of mixed infections are poorly understood. Here, we use two intracellular gut parasites (Microsporidia), one exotic and one co-evolved in the western honeybee (Apis mellifera), in an experiment in which either one or both parasites were administered either simultaneously or sequentially. We provide clear evidence of within-host competition; order of infection was an important determinant of the competitive outcome between parasites, with the first parasite significantly inhibiting the growth of the second, regardless of species. However, the strength of this 'priority effect' was highly asymmetric, with the exotic Nosema ceranae exhibiting stronger inhibition of Nosema apis than vice versa. Our results reveal an unusual asymmetry in parasite competition that is dependent on order of infection. When incorporated into a mathematical model of disease prevalence, we find asymmetric competition to be an important predictor of the patterns of parasite prevalence found in nature. Our findings demonstrate the wider significance of complex multi-host-multi-parasite interactions as drivers of host-pathogen community structure.

Population structure of a microparasite infecting Daphnia: spatio-temporal dynamics

Background

Detailed knowledge of spatial and temporal variation in the genetic population structure of hosts and parasites is required for understanding of host − parasite coevolution. As hot-spots of contemporary coevolution in natural systems are difficult to detect and long-term studies are restricted to few systems, additional population genetic data from various host − parasite systems may provide important insights into the topic. This is particularly true for parasites, as these players have been under-investigated so far due to the lower availability of suitable molecular markers. Here, we traced genetic variation (based on sequence variants in the internal transcribed spacer region, ITS) among seven geographically isolated populations of the ichthyosporean Caullerya mesnili, a common microparasite of the cladoceran Daphnia (here, the D. longispina hybrid complex). At three sites, we also studied parasite genetic variation over time (three to four sampling points) and tested for associations between parasite genotypes and host species.

Results

Parasite (and host) populations were significantly structured across space, indicating limited dispersal. Moreover, the frequency of parasite genotypes varied significantly over time, suggesting rapid evolutionary change in Caullerya. However, the distribution of parasite genotypes was similar across different host species, which might in turn have important consequences for parasite epidemiology.

Conclusions

The approach proposed here can be applied to track spatial and temporal changes in the population structure of other microparasite species for which sequence variation in the ITS or other highly variable genome regions has been documented but other types of polymorphic markers are lacking. Screening of parasite sequence variants allows for reliable detection of cross-species infections and, using advanced sequencing techniques in the near future, for detailed studies of parasite evolution in natural host − parasite systems.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-014-0247-3) contains supplementary material, which is available to authorized users.

Genetic architecture of resistance in Daphnia hosts against two species of host-specific parasites

Understanding the genetic architecture of host resistance is key for understanding the evolution of host-parasite interactions. Evolutionary models often assume simple genetics based on few loci and strong epistasis. It is unknown, however, whether these assumptions apply to natural populations. Using a quantitative trait loci (QTL) approach, we explore the genetic architecture of resistance in the crustacean Daphnia magna to two of its natural parasites: the horizontally transmitted bacterium Pasteuria ramosa and the horizontally and vertically transmitted microsporidium Hamiltosporidium tvaerminnensis. These two systems have become models for studies on the evolution of host-parasite interactions. In the QTL panel used here, Daphnia's resistance to P. ramosa is controlled by a single major QTL (which explains 50% of the observed variation). Resistance to H. tvaerminnensis horizontal infections shows a signature of a quantitative trait based in multiple loci with weak epistatic interactions (together explaining 38% variation). Resistance to H. tvaerminnensis vertical infections, however, shows only one QTL (explaining 13.5% variance) that colocalizes with one of the QTLs for horizontal infections. QTLs for resistance to Pasteuria and Hamiltosporidium do not colocalize. We conclude that the genetics of resistance in D. magna are drastically different for these two parasites. Furthermore, we infer that based on these and earlier results, the mechanisms of coevolution differ strongly for the two host-parasite systems. Only the Pasteuria-Daphnia system is expected to follow the negative frequency-dependent selection (Red Queen) model. How coevolution works in the Hamiltosporidium-Daphnia system remains unclear.

The expression of virulence for a mixed-mode transmitted parasite in a diapausing host

Many parasites survive harsh periods together with their hosts. Without the possibility of horizontal transmission during host diapause, parasite persistence depends entirely on host survival. We therefore hypothesize that a parasite should be avirulent during its host's diapausing stage. In contrast, the parasite may express higher virulence, i.e. parasite-induced fitness reduction of the host, during host life stages with good opportunities for horizontal transmission. Here we study the effects of a vertically and horizontally transmitted microsporidium parasite, Hamiltosporidium tvaerminnensis, on the quantity and survival of resting eggs of its host Daphnia magna. We find that the parasite did not affect egg volume, hatching success and time to hatching of the Daphnia's resting eggs, although it did strongly reduce the number of resting eggs produced by infected females, revealing high virulence during the non-diapause phase of the host's life cycle. These results also explain another aspect of this system - namely the strong decline in natural population prevalence across diapause. This decline is not caused by mortality in infected resting stages, as was previously hypothesized, but because infected female hosts produce lower rates of resting eggs. Together, these results help explain the epidemiological dynamics of a microsporidian disease and highlight the adaptive nature of life stage-dependent parasite virulence.

(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.

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.

Microsatellite and single-nucleotide polymorphisms indicate recurrent transitions to asexuality in a microsporidian parasite

Assessing the mode of reproduction of microparasites remains a difficult task because direct evidence for sexual processes is often absent and the biological covariates of sex and asex are poorly known. Species with geographically divergent modes of reproduction offer the possibility to explore some of these covariates, for example, the influence of life-history traits, mode of transmission and life-cycle complexity. Here, we present a phylogeographical study of a microsporidian parasite, which allows us to relate population genetic structure and mode of reproduction to its geographically diverged life histories. We show that in microsporidians from the genus Hamiltosporidium, that use the cladoceran Daphnia as host, an epidemic population structure has evolved, most probably since the last Ice Age. We partially sequenced three housekeeping genes (alpha tubulin, beta tubulin and hsp70) and genotyped seven microsatellite loci in 51 Hamiltosporidium isolates sampled within Europe and the Middle East. We found two phylogenetically related asexual parasite lines, one each from Fennoscandia and Israel, which share the unique ability of being transmitted both vertically and horizontally from Daphnia to Daphnia. The sexual forms cannot transmit horizontally among Daphnia, but presumably have a complex life cycle with a second host species. In spite of the similarities between the two asexual lineages, a clustering analysis based on microsatellite polymorphisms shows that asexual Fennoscandian parasites do not share ancestry with any other Hamiltosporidium that we have sampled. Moreover, allele sequence divergence at the hsp70 locus is twice as large in Fennoscandian than in Israeli parasites. Our results indicate that asexual reproduction evolved twice independently, first in Fennoscandian and more recently in the Israeli parasites. We conclude that the independent origin of asexuality in these two populations is associated with the altered parasite mode of transmission and the underlying dynamics of host populations.