The genome of an Encephalitozoon cuniculi type III strain reveals insights into the genetic diversity and mode of reproduction of a ubiquitous vertebrate pathogen

Bioinformatics analysis of the Microsporidia sp. MB genome: a malaria transmission-blocking symbiont of the Anopheles arabiensis mosquito

BACKGROUND

The use of microsporidia as a disease-transmission-blocking tool has garnered significant attention. Microsporidia sp. MB, known for its ability to block malaria development in mosquitoes, is an optimal candidate for supplementing malaria vector control methods. This symbiont, found in Anopheles mosquitoes, can be transmitted both vertically and horizontally with minimal effects on its mosquito host. Its genome, recently sequenced from An. arabiensis, comprises a compact 5.9 Mbp.

RESULTS

Here, we analyze the Microsporidia sp. MB genome, highlighting its major genomic features, gene content, and protein function. The genome contains 2247 genes, predominantly encoding enzymes. Unlike other members of the Enterocytozoonida group, Microsporidia sp. MB has retained most of the genes in the glycolytic pathway. Genes involved in RNA interference (RNAi) were also identified, suggesting a mechanism for host immune suppression. Importantly, meiosis-related genes (MRG) were detected, indicating potential for sexual reproduction in this organism. Comparative analyses revealed similarities with its closest relative, Vittaforma corneae, despite key differences in host interactions.

CONCLUSION

This study provides an in-depth analysis of the newly sequenced Microsporidia sp. MB genome, uncovering its unique adaptations for intracellular parasitism, including retention of essential metabolic pathways and RNAi machinery. The identification of MRGs suggests the possibility of sexual reproduction, offering insights into the symbiont's evolutionary strategies. Establishing a reference genome for Microsporidia sp. MB sets the foundation for future studies on its role in malaria transmission dynamics and host-parasite interactions.

Minimization and complete loss of general transcription factor proteins in the intracellular parasite Encephalitozoon cuniculi

Genome compaction is a common evolutionary feature of parasites. The unicellular, obligate intracellular parasite Encephalitozoon cuniculi has one of smallest known eukaryotic genomes, and is nearly four times smaller than its distant fungi relative, the budding yeast Saccharomyces cerevisiae. Comparison of the proteins encoded by compacted genomes to those encoded by larger genomes can reveal the most highly conserved features of the encoded proteins. In this study, we identified the proteins comprising the RNA polymerases and their corresponding general transcription factors by using several bioinformatic approaches to compare the transcription machinery of E. cuniculi and S. cerevisiae. Surprisingly, our analyses revealed an overall reduction in the size of the proteins comprising transcription machinery of E. cuniculi, which includes the loss of entire regions or functional domains from proteins, as well as the loss of entire proteins and complexes. Unexpectedly, we found that the E. cuniculi ortholog of Rpc37 (a RNA Polymerase III subunit) more closely resembles the H. sapiens ortholog of Rpc37 than the S. cerevisiae ortholog of Rpc37, in both size and structure. Overall, our findings provide new insight into the minimal core eukaryotic transcription machinery and help define the most critical features of Pol components and general transcription factors.

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.

Functional annotation of a divergent genome using sequence and structure-based similarity

BACKGROUND

Microsporidia are a large taxon of intracellular pathogens characterized by extraordinarily streamlined genomes with unusually high sequence divergence and many species-specific adaptations. These unique factors pose challenges for traditional genome annotation methods based on sequence similarity. As a result, many of the microsporidian genomes sequenced to date contain numerous genes of unknown function. Recent innovations in rapid and accurate structure prediction and comparison, together with the growing amount of data in structural databases, provide new opportunities to assist in the functional annotation of newly sequenced genomes.

RESULTS

In this study, we established a workflow that combines sequence and structure-based functional gene annotation approaches employing a ChimeraX plugin named ANNOTEX (Annotation Extension for ChimeraX), allowing for visual inspection and manual curation. We employed this workflow on a high-quality telomere-to-telomere sequenced tetraploid genome of Vairimorpha necatrix. First, the 3080 predicted protein-coding DNA sequences, of which 89% were confirmed with RNA sequencing data, were used as input. Next, ColabFold was used to create protein structure predictions, followed by a Foldseek search for structural matching to the PDB and AlphaFold databases. The subsequent manual curation, using sequence and structure-based hits, increased the accuracy and quality of the functional genome annotation compared to results using only traditional annotation tools. Our workflow resulted in a comprehensive description of the V. necatrix genome, along with a structural summary of the most prevalent protein groups, such as the ricin B lectin family. In addition, and to test our tool, we identified the functions of several previously uncharacterized Encephalitozoon cuniculi genes.

CONCLUSION

We provide a new functional annotation tool for divergent organisms and employ it on a newly sequenced, high-quality microsporidian genome to shed light on this uncharacterized intracellular pathogen of Lepidoptera. The addition of a structure-based annotation approach can serve as a valuable template for studying other microsporidian or similarly divergent species.

Comparative Genomics of Microsporidia

The microsporidia are a phylum of intracellular parasites that represent the eukaryotic cell in a state of extreme reduction, with genomes and metabolic capabilities embodying eukaryotic cells in arguably their most streamlined state. Over the past 20 years, microsporidian genomics has become a rapidly expanding field starting with sequencing of the genome of Encephalitozoon cuniculi, one of the first ever sequenced eukaryotes, to the current situation where we have access to the data from over 30 genomes across 20+ genera. Reaching back further in evolutionary history, to the point where microsporidia diverged from other eukaryotic lineages, we now also have genomic data for some of the closest known relatives of the microsporidia such as Rozella allomycis, Metchnikovella spp. and Amphiamblys sp. Data for these organisms allow us to better understand the genomic processes that shaped the emergence of the microsporidia as a group. These intensive genomic efforts have revealed some of the processes that have shaped microsporidian cells and genomes including patterns of genome expansions and contractions through gene gain and loss, whole genome duplication, differential patterns of invasion and purging of transposable elements. All these processes have been shown to occur across short and longer time scales to give rise to a phylum of parasites with dynamic genomes with a diversity of sizes and organisations.

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

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

Genome-Wide Analyses of Repeat-Induced Point Mutations in the Ascomycota

The Repeat-Induced Point (RIP) mutation pathway is a fungus-specific genome defense mechanism that mitigates the deleterious consequences of repeated genomic regions and transposable elements (TEs). RIP mutates targeted sequences by introducing cytosine to thymine transitions. We investigated the genome-wide occurrence and extent of RIP with a sliding-window approach. Using genome-wide RIP data and two sets of control groups, the association between RIP, TEs, and GC content were contrasted in organisms capable and incapable of RIP. Based on these data, we then set out to determine the extent and occurrence of RIP in 58 representatives of the Ascomycota. The findings were summarized by placing each of the fungi investigated in one of six categories based on the extent of genome-wide RIP. In silico RIP analyses, using a sliding-window approach with stringent RIP parameters, implemented simultaneously within the same genetic context, on high quality genome assemblies, yielded superior results in determining the genome-wide RIP among the Ascomycota. Most Ascomycota had RIP and these mutations were particularly widespread among classes of the Pezizomycotina, including the early diverging Orbiliomycetes and the Pezizomycetes. The most extreme cases of RIP were limited to representatives of the Dothideomycetes and Sordariomycetes. By contrast, the genomes of the Taphrinomycotina and Saccharomycotina contained no detectable evidence of RIP. Also, recent losses in RIP combined with controlled TE proliferation in the Pezizomycotina subphyla may promote substantial genome enlargement as well as the formation of sub-genomic compartments. These findings have broadened our understanding of the taxonomic range and extent of RIP in Ascomycota and how this pathway affects the genomes of fungi harboring it.

Encephalitozoon cuniculi infection in farmed rabbits in Egypt

BACKGROUND

Encephalitozoon cuniculi is an important microsporidian parasite with zoonotic potential. The present study highlights the impact of encephalitozoonosis on rabbit health in Egypt. Three rabbit farms in Giza, with a total of 16,400 rabbits were investigated due to occurrence of rabbits displaying clinical signs consistent with encephalitozoonosis.

RESULTS

Clinical signs observed during a 4 months observation period in 2018 included vestibular disease, paresis, limb paralysis, cataracts, phacoclastic uveitis, frequent urination, marked decrease in body weight and in some pregnant females, also repeated abortions. The total morbidity rates in adult and young rabbits were 76.7% and 81.5%, respectively. The highest mortality rate was recorded in offspring (12.3%), followed by dams (5.6%), and the lowest recorded mortality rate was in males (0.04%). Post-mortem findings included enteritis, pale enlarged kidneys, congested leptomeninges, focal brain necrosis, and endometrial congestion. Histopathological examination revealed nonsuppurative meningoencephalitis and glial nodules with central necrosis in the brain, vacuolation and necrosis of renal tubular epithelium, and corneal ulceration and ruptured lens capsule with fragmentation of lenticular fibres. E. cuniculi were observed in the brain, retinal ganglion cells, kidneys, and liver. Transmission electron microscopy examination revealed the presence of different developmental stages of E. cuniculi in the brain and kidney. Presence of E. cuniculi was confirmed by conventional polymerase chain reaction using a universal 16S gene for Encephalitozoon spp. followed by sequencing and sequence analysis.

CONCLUSIONS

The presence of E. cuniculi in rabbits was confirmed at three farms in Egypt. Nervous signs and ocular lesions were the most predominant findings in these farms.

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.

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.

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.

Evolution of Dicer and Argonaute orthologs in microsporidian parasites

Microsporidia are a group of intracellular parasites which infect animal hosts. The infection can broadly influence the hosts' metabolism, growth as well as immune responses. Recently, a functional RNAi pathway was suggested from the microsporidia parasite Nosema ceranae, whereby the gene Dicer showed strong impact on spore proliferation. Based on sequenced microsporidian species, the RNAi gene orthologs have only been annotated for a few species. In order to study the selection of RNAi gene Dicer and Argonaute orthologs from microsporidian genomes, a phylogenetic analysis was performed based on single copy orthologs of 21 microsporidian parasite species. Of the 21 studied parasite species, 11 parasite species maintained Dicer and Argonaute orthologs, which were further used to build the gene trees. The gene Dicer and Argonaute orthologs were either both maintained or both lost. The topology structures between the phylogenetic Dicer, Argonaute and species trees were consistent. The results suggest that the gene Dicer and Argonaute were selected as a unit, which were selectively maintained/lost during the lineage divergence. The study provides general insights on the selection of RNAi pathway in microsporidian parasites and the evolution of parasitism.

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.

Genome Analysis of Pseudoloma neurophilia: A Microsporidian Parasite of Zebrafish (Danio rerio)

Microsporidia are highly successful parasites that infect virtually all known animal lineages, including the model Danio rerio (zebrafish). The widespread use of this aquatic model for biomedical research has resulted in an unexpected increase in infections from the microsporidium Pseudoloma neurophilia, which can lead to significant physical, behavioral, and immunological modifications, resulting in nonprotocol variation during experimental procedures. Here, we seek to obtain insights into the biology of P. neurophilia by investigating its genome content, which was obtained from only 29 nanograms of DNA using the MiSeq technology and paired-end Illumina sequencing. We found that the genome of P. neurophilia is phylogenetically and genetically related to other fish-microsporidians, but features unique to this intracellular parasite are also found. The small 5.25-Mb genome assembly includes 1,139 unique open-reading frames and an unusually high number of transposable elements for such a small genome. Investigations of intragenomic diversity also provided strong indications that the mononucleate nucleus of this species is diploid. Overall, our study provides insights into the dynamics of microsporidian genomes and a solid sequence reference to be used in future studies of host-parasite interactions using the zebrafish D. rerio and P. neurophilia as a model.

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.