Filling gaps in the microsporidian tree: rDNA phylogeny of Chytridiopsis typographi (Microsporidia: Chytridiopsida)

Inside the Host: Understanding the Evolutionary Trajectories of Intracellular Parasitism

This review explores the origins of intracellular parasitism, an intriguing facet of symbiosis, where one organism harms its host, potentially becoming deadly. We focus on three distantly related groups of single-celled eukaryotes, namely Kinetoplastea, Holomycota, and Apicomplexa, which contain multiple species-rich lineages of intracellular parasites. Using comparative analysis of morphological, physiological, and molecular features of kinetoplastids, microsporidians, and sporozoans, as well as their closest free-living relatives, we reveal the evolutionary trajectories and adaptations that enabled the transition to intracellular parasitism. Intracellular parasites have evolved various efficient mechanisms for host acquisition and exploitation, allowing them to thrive in a variety of hosts. Each group has developed unique features related to the parasitic lifestyle, involving dedicated protein families associated with host cell invasion, survival, and exit. Indeed, parallel evolution has led to distinct lineages of intracellular parasites employing diverse traits and approaches to achieve similar outcomes.

Classes and phyla of the kingdom Fungi

Fungi are one of the most diverse groups of organisms with an estimated number of species in the range of 2–3 million. The higher-level ranking of fungi has been discussed in the framework of molecular phylogenetics since Hibbett et al., and the definition and the higher ranks (e.g., phyla) of the ‘true fungi’ have been revised in several subsequent publications. Rapid accumulation of novel genomic data and the advancements in phylogenetics now facilitate a robust and precise foundation for the higher-level classification within the kingdom. This study provides an updated classification of the kingdom Fungi, drawing upon a comprehensive phylogenomic analysis of Holomycota, with which we outline well-supported nodes of the fungal tree and explore more contentious groupings. We accept 19 phyla of Fungi, viz. Aphelidiomycota, Ascomycota, Basidiobolomycota, Basidiomycota, Blastocladiomycota, Calcarisporiellomycota, Chytridiomycota, Entomophthoromycota, Entorrhizomycota, Glomeromycota, Kickxellomycota, Monoblepharomycota, Mortierellomycota, Mucoromycota, Neocallimastigomycota, Olpidiomycota, Rozellomycota, Sanchytriomycota, and Zoopagomycota. In the phylogenies, Caulochytriomycota resides in Chytridiomycota; thus, the former is regarded as a synonym of the latter, while Caulochytriomycetes is viewed as a class in Chytridiomycota. We provide a description of each phylum followed by its classes. A new subphylum, Sanchytriomycotina Karpov is introduced as the only subphylum in Sanchytriomycota. The subclass Pneumocystomycetidae Kirk et al. in Pneumocystomycetes, Ascomycota is invalid and thus validated. Placements of fossil fungi in phyla and classes are also discussed, providing examples.

O short-branch Microsporidia, where art thou? Identifying diversity hotspots for future sampling

Short-branch Microsporidia were previously shown to form a basal grade within the expanded Microsporidia clade and to branch near the classical, long-branch Microsporidia. Although they share simpler versions of some morphological characteristics, they do not show accelerated evolutionary rates, making them ideal candidates to study the evolutionary trajectories that have led to long-branch microsporidian unique characteristics. However, most sequences assigned to the short-branch Microsporidia are undescribed, novel environmental lineages for which the identification requires knowledge of where they can be found. To direct future isolation, we used the EukBank database of the global UniEuk initiative that contains the majority of the publicly available environmental V4 SSU rRNA gene sequences of protists. The curated OTU table and corresponding metadata were used to evaluate the occurrence of short-branch Microsporidia across freshwater, hypersaline, marine benthic, marine pelagic, and terrestrial environments. Presence-absence analyses infer that short-branch Microsporidia are most abundant in freshwater and terrestrial environments, and alpha- and beta-diversity measures indicate that focusing our sampling effort on these two environments would cover a large part of their overall diversity. These results can be used to coordinate future isolation and sampling campaigns to better understand the enigmatic evolution of microsporidians' unique characteristics.

Quantification of the dark fungal taxon Cryptomycota using qPCR

Fungi are present in a wide variety of natural environments, and in the last years, various studies have shown that they are quite abundant in aquatic ecosystems. In addition, a whole new highly diverse phylum, the Cryptomycota, was discovered. Nevertheless, research on aquatic fungi and a detailed evaluation of their functions and distribution are still sparse. One of the main reasons is a limitation in reliable identification and quantification methods. To bridge part of the research gap, this study aims to implement a quantitative PCR method to detect and quantify the newly discovered phylum. We developed and validated a Cryptomycota-specific qPCR primer pair targeting the 5.8S region that detects the majority of Cryptomycota, but Microsporidia. The resulting amplicon is 102 bp long. We used different environmental samples to evaluate the primer pair, various fungal sequences as negative control and positive control sequences. Obtained amplicons were sequenced using Illumina, and the obtained ASVs were all classified as Cryptomycota. The qPCR method works reliably and specifically for the quantification of Cryptomycota in environmental samples.

Microsporidians (Microsporidia) parasitic on mosquitoes (Culicidae) in central Europe are often multi-host species

Microsporidians (Microsporidia) are a diverse group of obligate and intracellular parasites of eukaryotes. There is evidence that the real species diversity in the phylum could be greatly underestimated, especially for microsporidians parasitic on invertebrates. Mosquitoes (Culicidae) are among very important microsporidian host groups. However, to date, no extensive survey on the prevalence of microsporidians in European mosquitoes has been performed. Here, we used mosquitoes collected in west-central Poland and a metabarcoding approach to examine the prevalence and diversity of microsporidian species among European mosquitoes. We found that up to one-third of mosquitoes in Europe may be infected with at least 13 microsporidian species belonging to the genera Amblyospora, Hazardia, Encephalitozoon, Enterocytospora, and Nosema and the holding genus Microsporidium. The lack of a difference in microsporidian prevalence between mosquito sexes implies that other factors, e.g., temperature or humidity, affect microsporidian occurrence in adult mosquitoes. Each microsporidian species was found in at least three mosquito species, which suggests that these microsporidians are polyxenic rather than monoxenic parasites. The co-occurrence of at least two different microsporidian species was found in 3.6% of host individuals. The abundance of microsporidian DNA sequences suggests interactions between co-occurring parasites; however, these results should be confirmed by microscopic and quantitative methods. In addition, further histological research is required to describe Microsporidium sp. PL01 or match its DNA to that of an already described species.

Microsporidia: a new taxonomic, evolutionary, and ecological synthesis

Microsporidian diversity is vast. There is a renewed drive to understand how microsporidian pathological, genomic, and ecological traits relate to their phylogeny. We comprehensively sample and phylogenetically analyse 125 microsporidian genera for which sequence data are available. Comparing these results with existing phylogenomic analyses, we suggest an updated taxonomic framework to replace the inconsistent clade numbering system, using informal taxonomic names: Glugeida (previously clades 5/3), Nosematida (4a), Enterocytozoonida (4b), Amblyosporida (3/5), Neopereziida (1), and Ovavesiculida (2). Cellular, parasitological, and ecological traits for 281 well-defined species are compared with identify clade-specific patterns across long-branch Microsporidia. We suggest that future taxonomic circumscriptions of Microsporidia should involve additional markers (SSU/ITS/LSU), and that a comprehensive suite of phenotypic and ecological traits help to predict broad microsporidian functional and lineage diversity.

Proteomic Analysis of Spore Surface Proteins and Characteristics of a Novel Spore Wall Protein and Biomarker, EhSWP3, from the Shrimp Microsporidium Enterocytozoon hepatopenaei (EHP)

Enterocytozoon hepatopenaei, a spore-forming and obligate intracellular microsporidium, mainly infects shrimp and results in growth retardation and body length variation, causing huge economic losses to the Asian shrimp aquaculture industry. However, the lack of a full understanding of the surface proteins of spores associated with host infection has hindered the development of technologies for the detection of EHP. In this study, the surface proteins of EHP spores were extracted using the improved SDS method, and 130 proteins were identified via LC-MS/MS analysis. Bioinformatic analysis revealed that these proteins were enriched in biological processes (67), cellular components (62), and molecular functions (71) based on GO terms. KEGG pathway analysis showed that 20 pathways, including the proteasome (eight proteins) and the fatty acid metabolism (15 proteins), were enriched. Among 15 high-abundance surface proteins (HASPs), EhSWP3 was identified as a novel spore wall protein (SWP), and was localized on the endospore of the EHP spores with an indirect immunofluorescence and immunoelectron microscopy assay. Polyclonal antibodies against EhSWP3 showed strong species specificity and high sensitivity to the hepatopancreas of EHP-infected shrimp. As a specific high-abundance protein, EhSWP3 is therefore a promising target for the development of immunoassay tools for EHP detection, and may play a crucial role in the invasion of EHP into the host.

Metabarcoding reveals low prevalence of microsporidian infections in castor bean tick (Ixodes ricinus)

BACKGROUND

Microsporidia is a large group of eukaryotic obligate intracellular spore-forming parasites, of which 17 species can cause microsporidiosis in humans. Most human-infecting microsporidians belong to the genera Enterocytozoon and Encephalitozoon. To date, only five microsporidian species, including Encephalitozoon-like, have been found in hard ticks (Ixodidae) using microscopic methods, but no sequence data are available for them. Furthermore, no widespread screening for microsporidian-infected ticks based on DNA analysis has been carried out to date. Thus, in this study, we applied a recently developed DNA metabarcoding method for efficient microsporidian DNA identification to assess the role of ticks as potential vectors of microsporidian species causing diseases in humans.

METHODS

In total, 1070 (493 juvenile and 577 adult) unfed host-seeking Ixodes ricinus ticks collected at urban parks in the city of Poznan, Poland, and 94 engorged tick females fed on dogs and cats were screened for microsporidian DNA. Microsporidians were detected by PCR amplification and sequencing of the hypervariable V5 region of 18S rRNA gene (18S profiling) using the microsporidian-specific primer set. Tick species were identified morphologically and confirmed by amplification and sequencing of the shortened fragment of cytochrome c oxidase subunit I gene (mini-COI).

RESULTS

All collected ticks were unambiguously assigned to I. ricinus. Potentially zoonotic Encephalitozoon intestinalis was identified in three fed ticks (3.2%) collected from three different dogs. In eight unfed host-seeking ticks (0.8%), including three males (1.1%), two females (0.7%) and three nymphs (0.7%), the new microsporidian sequence representing a species belonging to the genus Endoreticulatus was identified.

CONCLUSIONS

The lack of zoonotic microsporidians in host-seeking ticks suggests that I. ricinus is not involved in transmission of human-infecting microsporidians. Moreover, a very low occurrence of the other microsporidian species in both fed and host-seeking ticks implies that mechanisms exist to defend ticks against infection with these parasites.

Impact of Genome Reduction in Microsporidia

Microsporidia represent an evolutionary outlier in the tree of life and occupy the extreme edge of the eukaryotic domain with some of their biological features. Many of these unicellular fungi-like organisms have reduced their genomic content to potentially the lowest limit. With some of the most compacted eukaryotic genomes, microsporidia are excellent model organisms to study reductive evolution and its functional consequences. While the growing number of sequenced microsporidian genomes have elucidated genome composition and organization, a recent increase in complementary post-genomic studies has started to shed light on the impacts of genome reduction in these unique pathogens. This chapter will discuss the biological framework enabling genome minimization and will use one of the most ancient and essential macromolecular complexes, the ribosome, to illustrate the effects of extreme genome reduction on a structural, molecular, and cellular level. We outline how reductive evolution in microsporidia has shaped DNA organization, the composition and function of the ribosome, and the complexity of the ribosome biogenesis process. Studying compacted mechanisms, processes, or macromolecular machines in microsporidia illuminates their unique lifestyle and provides valuable insights for comparative eukaryotic structural biology.

Insights into Microsporidia Evolution from Early Diverging Microsporidia

Microsporidia have drastically modified genomes and cytology resulting from their high level of adaptation to intracytoplasmic parasitism. Their origins, which had long remained enigmatic, were placed within the line of Rozella, a primitive endoparasitic chytrid. These origins became more and more refined with the discovery of various parasites morphologically similar to the primitive lines of microsporidia (Metchnikovellids and Chytridiopsids) but which possess fungal-like genomes and functional mitochondria. These various parasites turn out to be distinct missing links between a large assemblage of chytrid-like rozellids and the true microsporidians, which are actually a very evolved branch of the rozellids themselves. The question of how to consider the historically known Microsporidia and the various microsporidia-like organisms within paraphyletic rozellids is discussed.

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.

A phylogenetic framework to investigate the microsporidian communities through metabarcoding and its application to lake ecosystems

Microsporidia are obligate intracellular eukaryotic parasites known to parasitize many species of the animal kingdom as well as some protists. However, their diversity is underestimated, in part as a consequence of the failure of 'universal' primers to detect them in metabarcoding studies. Besides, due to the inconsistency between taxonomy and phylogenetic data, available databases may assign incorrectly sequences obtained with high-throughput sequencing. In this work, we developed a comprehensive reference database which positions microsporidian SSU rRNA gene sequences within a coherent ranked phylogenetic framework. We used this phylogenetic framework to study the microsporidian diversity in lacustrine ecosystems, focusing on < 150 μm planktonic size fractions. Our analysis shows a high diversity of Microsporidia, with the identification of 1531 OTUs distributed within seven clades, of which 76% were affiliated to clade IV2 and 20% to clade I (nomenclature presented hereby). About a quarter of the obtained sequences shared less than 85% identity to the closest known species, which might represent undescribed genera or families infecting small hosts. Variations in the abundance of Microsporidia were recorded between the two lakes sampled and across the sampling period, which might be explained by spatio-temporal variations of their potential hosts such as microeukaryotes and metazooplankton.

Molecular phylogeny and new light microscopic data of Metchnikovella spiralis (Microsporidia: Metchnikovellidae), a hyperparasite of eugregarine Polyrhabdina sp. from the polychaete Pygospio elegans

Metchnikovellids are a deep-branching group of microsporidia, parasites of gregarines inhabiting the alimentary tract of polychaetes and some other invertebrates. The diversity and phylogeny of these hyperparasites remain poorly studied. Modern descriptions and molecular data are still lacking for many species. The results of a light microscopy study and molecular data for Metchnikovella spiralis Sokolova et al., 2014, a hyperparasite of the eugregarine Polyrhabdina sp., isolated from the polychaete Pygospio elegans, were obtained. The original description of M. spiralis was based primarily on the analysis of stained preparations and transmission electron microscopy images. Here, the species description was complemented with the results of in vivo observations and phylogenetic analysis based on the SSU rRNA gene. It was shown that in this species, free sporogony precedes sac-bound sporogony, as it occurs in the life cycle of most other metchnikovellids. Spore sacs are entwined with spirally wound cords, and possess only one polar plug. Phylogenetic analyses did not group M. spiralis with M. incurvata, another metchnikovellid from the same gregarine species, but placed it as a sister branch to Amphiacantha. The paraphyletic nature of the genus Metchnikovella was discussed. The taxonomic summary for M. spiralis was emended.

Infection effects of the new microsporidian species Tubulinosema suzukii on its host Drosophila suzukii

Microsporidian infections of insects are important natural constraints of population growth, often reducing lifespan, fecundity and fertility of the infected host. The recently discovered Tubulinosema suzukii infects Drosophila suzukii (spotted wing drosophila, SWD), an invasive pest of many fruit crops in North America and Europe. In laboratory tests, fitness effects on larval and adult stages were explored. High level infection after larval treatment caused up to 70% pupal mortality, a decreased lifespan and a 70% reduced oviposition of emerging adults in biparental infection clusters. A shift to higher proportion of female offspring compared to controls suggested a potential parthenogenetic effect after microsporidian infection. A clear sex-linkage of effects was noted; females were specifically impaired, as concluded from fecundity tests with only infected female parents. Additive effects were noted when both parental sexes were infected, whereas least effects were found with only infected male parents, though survival of males was most negatively affected if they were fed with T. suzukii spores in the adult stage. Although most negative effects on fitness parameters were revealed after larval treatment, infection of offspring was never higher than 4%, suggesting limited vertical transmission. For that reason, a self-reliant spread in natural SWD populations would probably only occur by spore release from cadavers or frass.

Revisiting the phylogeny of microsporidia

Canonical microsporidians are a group of obligate intracellular parasites of a wide range of hosts comprising ~1,300 species of >220 genera. Microsporidians are related to fungi, and many characterised and uncharacterized groups closely related to them have been discovered recently, filling the knowledge gaps between them. These groups assigned to the superphylum Opisthosporidia have provided several important insights into the evolution of diverse intracellular parasitic lineages within the tree of eukaryotes. The most studied among opisthosporidians, canonical microsporidians, were known to science more than 160 years ago, however, the classification of canonical Microsporidia has been challenging due to common morphological homoplasy, and accelerated evolutionary rates. Instead of morphological characters, ssrRNA sequences have been used as the primary data for the classification of canonical microsporidians. Previous studies have produced a useful backbone of the microsporidian phylogeny, but provided only some nodal support, causing some confusion. Here, we reconstructed phylogenetic trees of canonical microsporidians using Bayesian and Maximum Likelihood inferences. We included rRNA sequences of 126 described/named genera, by far the broadest taxon coverage to date. Overall, our trees show similar topology and recovered four of the five main clades demonstrated in previous studies (Clades 1, 3, 4 and 5). Family level clades were well resolved within each major clade, but many were discordant with the recently revised classification. Therefore, revision and some reshuffling, especially within and between Clades 1 and 3 are required. We also reconstructed phylogenetic trees of Opisthosporidia to better integrate the evolutionary history of canonical microsporidians in a broader context. We discuss several traits shared only by canonical microsporidians that may have contributed to their striking ecological success in diverse metazoans. More targeted studies on the neglected host groups will be of value for a better understanding of the evolutionary history of these interesting intracellular parasites.

Evolutionary relationships of Metchnikovella dogieli Paskerova et al., 2016 (Microsporidia: Metchnikovellidae) revealed by multigene phylogenetic analysis

The species Metchnikovella dogieli (Paskerova et al. Protistology 10:148-157, 2016) belongs to one of the early diverging microsporidian groups, the metchnikovellids (Microsporidia: Metchnikovellidae). In relation to typical ('core') microsporidia, this group is considered primitive. The spores of metchnikovellids have no classical polar sac-anchoring disk complex, no coiled polar tube, no posterior vacuole, and no polaroplast. Instead, they possess a short thick manubrium that expands into a manubrial cistern. These organisms are hyperparasites; they infect gregarines that parasitise marine invertebrates. M. dogieli is a parasite of the archigregarine Selenidium pygospionis (Paskerova et al. Protist 169:826-852, 2018), which parasitises the polychaete Pygospio elegans. This species was discovered in samples collected in the silt littoral zone at the coast of the White Sea, North-West Russia, and was described based on light microscopy. No molecular data are available for this species, and the publicly accessible genomic data for metchnikovellids are limited to two species: M. incurvata Caullery & Mesnil, 1914 and Amphiamblys sp. WSBS2006. In the present study, we applied single-cell genomics methods with whole-genome amplification to perform next-generation sequencing of M. dogieli genomic DNA. We performed a phylogenetic analysis based on the SSU rRNA gene and reconstructed a multigene phylogeny using a concatenated alignment that included 46 conserved single-copy protein domains. The analyses recovered a fully supported clade of metchnikovellids as a basal group to the core microsporidia. Two members of the genus Metchnikovella did not form a clade in our tree. This may indicate that this genus is paraphyletic and requires revision.

Microsporidia Can Acquire Lamin-like Intermediate Filaments and Cell Adhesion Catenin-cadherin Complexes from the Host (?)

On their spore surfaces, Microsporidia often develop a canopy of filaments with characteristics of intermediate filaments (IF), as we demonstrated in previous studies on Thelohania sp., Ameson michaelis, and Spraguea lophii. Genomic studies indicate that among invertebrates, lamins that may localize in the cytoplasm or nucleus, are the only known IF type. These IFs can bind to the substrate containing cell adhesion molecules (CAMs) cadherins, associated with β and γ catenins. The objects of this study were to determine whether microsporidia have CAMs with the attached IFs on their envelopes and to find out if these proteins are provided by the host. An examination was made for localization of lamins and CAMs on the spores of the mentioned above species and Anncaliia algerae, plus in the host animals. Then, we determined whether the spores of A. michaelis and A. algerae could bind vertebrate nuclear lamin onto the spore surface. We also tested transgenic Drosophila melanogaster stocks bearing cadherin-GFP to see whether developing A. algerae parasites in these hosts could acquire host CAMs. The tests were positive for all these experiments. We hypothesize that microsporidia are able to acquire host lamin IFs and cell adhesion catenin-cadherin complexes from the host.

A new method of metabarcoding Microsporidia and their hosts reveals high levels of microsporidian infections in mosquitoes (Culicidae)

DNA metabarcoding offers new perspectives, especially with regard to the high‐throughput identification and diagnostics of pathogens. Microsporidia are an example of widely distributed, opportunistic and pathogenic microorganisms in which molecular identification is important for both environmental research and clinical diagnostics. We have developed a method for parallel detection of both microsporidian infection and the host species. We designed new primer sets: one specific for the classical Microsporidia (targeting the hypervariable V5 region of small subunit [ssu] rDNA), and a second one targeting a shortened fragment of the COI gene (standard metazoan DNA‐barcode); both markers are well suited for next generation sequencing. Analysis of the ssu rDNA data set representing 607 microsporidian species (120 genera) indicated that the V5 region enables identification of >98% species in the data set (596/607). To test the method, we used microsporidians that infect mosquitoes in natural populations. Using mini‐COI data, all field‐collected mosquitoes were unambiguously assigned to seven species; among them almost 60% of specimens were positive for at least 11 different microsporidian species, including a new microsporidian ssu rDNA sequence (Microsporidium sp. PL01). Phylogenetic analysis showed that this species belongs to one of the two main clades in the Terresporidia. We found a high rate of microsporidian co‐infections (9.4%). The numbers of sequence reads for the operational taxonomic units suggest that the occurrence of Nosema spp. in co‐infections could benefit them; however, this observation should be retested using a more intensive host sampling. Our results show that DNA barcoding is a rapid and cost‐effective method for deciphering sample diversity in greater resolution, including the hidden biodiversity that may be overlooked using classical methodology.

Solving an old enigma: Morellospora saccamoebae gen. nov., sp. nov. (Rozellomycota), a Sphaerita-like parasite of free-living amoebae

The Rozellomycota form a lineage basal or sister to the Fungi, ancestor of Microsporidia. Their biodiversity is very rich but remains poorly characterized. The few known species are all parasites, whether of water molds and algae (Rozella), crustaceans (Mitosporidium), or as endonuclear parasites of amoebae (Nucleophaga, Paramicrosporidium). Since the nineteenth century, intracytoplasmic parasites of various protozoa have been described as species of the same genus Sphaerita. However, it was later thought possible to separate these parasites into at least two distinct groups, those forming flagellated zoospores, prevalent in Euglena and other flagellates, and those forming immobile spores, found mainly in free-living and endozoic amoebae. Herein, we report the recovery of a strain of the free-living amoeba species Saccamoeba lacustris, naturally infected by an intracytoplasmic parasite, which under light microscope has a morphology consistent with that of Sphaerita. Biomolecular analyses were thus performed. Our results show that the intracytoplasmic parasite of Saccamoeba belongs to the same subgroup of Mitosporidium and that it forms a new genus within Rozellomycota, Morellospora, that corresponds to the former spore-forming Sphaerita-like parasites of amoebae.

Putative group I introns in the eukaryote nuclear internal transcribed spacers

Group I introns are mobile genetic elements that interrupt genes encoding proteins and RNAs. In the rRNA operon, introns can insert in the small subunit (SSU) and large subunit (LSU) of a wide variety of protists and various prokaryotes, but they were never found in the ITS region. In this study, unusually long ITS regions of fungi and closely related unicellular organisms (Polychytrium aggregatum, Mitosporidium daphniae, Amoeboaphelidium occidentale and Nuclearia simplex) were analysed. While the insertion of repeats is responsible for long ITS in other eukaryotes, the increased size of the sequences analysed herein seems rather due to the presence of introns in ITS-1 or ITS-2. The identified insertions can be folded in secondary structures according to group I intron models, and they cluster within introns in conserved core-based phylogeny. In addition, for Mitosporidium, Amoeboaphelidium and Nuclearia, more conventional ITS-2 structures can be deduced once spacer introns are removed. Sequences of five shark species were also analysed for their structure and included in phylogeny because of unpublished work reporting introns in their ITS, obtaining congruent results. Overall, the data presented herein indicate that spacer regions may contain introns.