Gregarine single-cell transcriptomics reveals differential mitochondrial remodeling and adaptation in apicomplexans

Plasmodium LCCL domain-containing modular proteins have their origins in the ancestral alveolate

Plasmodium species encode a unique set of six modular proteins named LCCL lectin domain adhesive-like proteins (LAPs) that operate as a complex and that are essential for malaria parasite transmission from mosquito to vertebrate. LAPs possess complex architectures obtained through unique assemblies of conserved domains associated with lipid, protein and carbohydrate interactions, including the name-defining LCCL domain. Here, we assessed the prevalence of Plasmodium LAP orthologues across eukaryotic life. Our findings show orthologous conservation in all apicomplexans, with lineage-specific repertoires acquired through differential lap gene loss and duplication. Besides Apicomplexa, LAPs are found in their closest relatives: the photosynthetic chromerids, which encode the broadest repertoire including a novel membrane-bound LCCL protein. LAPs are notably absent from other alveolate lineages (dinoflagellates, perkinsids and ciliates), but are encoded by predatory colponemids, a sister group to the alveolates. These results reveal that the LAPs are much older than previously thought and pre-date not only the Apicomplexa but the Alveolata altogether.

Parallel functional reduction in the mitochondria of apicomplexan parasites

Extreme functional reduction of mitochondria has taken place in parallel in many distantly related lineages of eukaryotes, leading to a number of recurring metabolic states with variously lost electron transport chain (ETC) complexes, loss of the tricarboxylic acid (TCA) cycle, and/or loss of the mitochondrial genome. The resulting mitochondria-related organelles (MROs) are generally structurally reduced and in the most extreme cases barely recognizable features of the cell with no role in energy metabolism whatsoever (e.g., mitosomes, which generally only make iron-sulfur clusters). Recently, a wide diversity of MROs were discovered to be hiding in plain sight: in gregarine apicomplexans. This diverse group of invertebrate parasites has been known and observed for centuries, but until recent applications of culture-free genomics, their mitochondria were unremarkable. The genomics, however, showed that mitochondrial function has reduced in parallel in multiple gregarine lineages to several different endpoints, including the most reduced mitosomes. Here we review this remarkable case of parallel evolution of MROs, and some of the interesting questions this work raises.

Phylogenomics reveals Adeleorina are an ancient and distinct subgroup of Apicomplexa

Apicomplexans are a diverse phylum of unicellular eukaryotes that share obligate relationships with terrestrial and aquatic animal hosts. Many well-studied apicomplexans are responsible for several deadly zoonotic and human diseases, most notably malaria caused by Plasmodium. Interest in the evolutionary origin of apicomplexans has also spurred recent work on other more deeply-branching lineages, especially gregarines and sister groups like squirmids and chrompodellids. But a full picture of apicomplexan evolution is still lacking several lineages, and one major, diverse lineage that is notably absent is the adeleorinids. Adeleorina apicomplexans comprises hundreds of described species that infect invertebrate and vertebrate hosts across the globe. Although historically considered coccidians, phylogenetic trees based on limited data have shown conflicting branch positions for this subgroup, leaving this question unresolved. Phylogenomic trees and large-scale analyses comparing cellular functions and metabolism between major subgroups of apicomplexans have not incorporated Adeleorina because only a handful of molecular markers and a couple organellar genomes are available, ultimately excluding this group from contributing to our understanding of apicomplexan evolution and biology. To address this gap, we have generated complete genomes from mitochondria and plastids, as well as multiple deep-coverage single-cell transcriptomes of nuclear genes from two Adeleorina species, Klossia helicina and Legerella nova, and inferred a 206-protein phylogenomic tree of Apicomplexa. We observed distinct structures reported in species descriptions as remnant host structures surrounding adeleorinid oocysts. Klossia helicina and L. nova branched, as expected, with monoxenous adeleorinids within the Adeleorina and their mitochondrial and plastid genomes exhibited similarity to published organellar adeleorinid genomes. We show with a phylogeneomic tree and subsequent phylogenomic analyses that Adeleorina are not closely related to any of the currently sampled apicomplexan subgroups, and instead fall as a sister to a large clade encompassing Coccidia, Protococcidia, Hematozoa, and Nephromycida, collectively. This resolves Adeleorina as a key independently-branching group, separate from coccidians, on the tree of Apicomplexa, which now has all known major lineages sampled.

Encyclopedia of Family A DNA Polymerases Localized in Organelles: Evolutionary Contribution of Bacteria Including the Proto-Mitochondrion

DNA polymerases synthesize DNA from deoxyribonucleotides in a semiconservative manner and serve as the core of DNA replication and repair machinery. In eukaryotic cells, there are 2 genome-containing organelles, mitochondria, and plastids, which were derived from an alphaproteobacterium and a cyanobacterium, respectively. Except for rare cases of genome-lacking mitochondria and plastids, both organelles must be served by nucleus-encoded DNA polymerases that localize and work in them to maintain their genomes. The evolution of organellar DNA polymerases has yet to be fully understood because of 2 unsettled issues. First, the diversity of organellar DNA polymerases has not been elucidated in the full spectrum of eukaryotes. Second, it is unclear when the DNA polymerases that were used originally in the endosymbiotic bacteria giving rise to mitochondria and plastids were discarded, as the organellar DNA polymerases known to date show no phylogenetic affinity to those of the extant alphaproteobacteria or cyanobacteria. In this study, we identified from diverse eukaryotes 134 family A DNA polymerase sequences, which were classified into 10 novel types, and explored their evolutionary origins. The subcellular localizations of selected DNA polymerases were further examined experimentally. The results presented here suggest that the diversity of organellar DNA polymerases has been shaped by multiple transfers of the PolI gene from phylogenetically broad bacteria, and their occurrence in eukaryotes was additionally impacted by secondary plastid endosymbioses. Finally, we propose that the last eukaryotic common ancestor may have possessed 2 mitochondrial DNA polymerases, POP, and a candidate of the direct descendant of the proto-mitochondrial DNA polymerase I, rdxPolA, identified in this study.

Single cell transcriptomics reveals UAR codon reassignment in Palmarella salina (Metopida, Armophorea) and confirms Armophorida belongs to APM clade

Anaerobes have emerged in several major lineages of ciliates, but the number of independent transitions to anaerobiosis among ciliates is unknown. The APM clade (Armophorea, Muranotrichea, Parablepharismea) represents the largest clade of obligate anaerobes among ciliates and contains free-living marine and freshwater representatives as well as gut endobionts of animals. The evolution of APM group has only recently started getting attention, and our knowledge on its phylogeny and genetics is still limited to a fraction of taxa. While ciliates portray a wide array of alternatives to the standard genetic code across numerous classes, the APM ciliates were considered to be the largest group using exclusively standard nuclear genetic code. In this study, we present a pan-ciliate phylogenomic analysis with emphasis on the APM clade, bringing the first phylogenomic analysis of the family Tropidoatractidae (Armophorea) and confirming the position of Armophorida within Armophorea. We include five newly sequenced single cell transcriptomes from marine, freshwater, and endobiotic APM ciliates - Palmarella salina, Anteclevelandella constricta, Nyctotherus sp., Caenomorpha medusula, and Thigmothrix strigosa. We report the first discovery of an alternative nuclear genetic code among APM ciliates, used by Palmarella salina (Tropidoatractidae, Armophorea), but not by its close relative, Tropidoatractus sp., and provide a comparative analysis of stop codon identity and frequency indicating the precedency to the UAG codon loss/reassignment over the UAA codon reassignment in the specific ancestor of Palmarella. Comparative genomic and proteomic studies of this group may help explain the constraints that underlie UAR stop-to-sense reassignment, the most frequent type of alternative nuclear genetic code, not only in ciliates, but eukaryotes in general.

Multiple pathways for glucose phosphate transport and utilization support growth of Cryptosporidium parvum

Cryptosporidium parvum is an obligate intracellular parasite with a highly reduced mitochondrion that lacks the tricarboxylic acid cycle and the ability to generate ATP, making the parasite reliant on glycolysis. Genetic ablation experiments demonstrated that neither of the two putative glucose transporters CpGT1 and CpGT2 were essential for growth. Surprisingly, hexokinase was also dispensable for parasite growth while the downstream enzyme aldolase was required, suggesting the parasite has an alternative way of obtaining phosphorylated hexose. Complementation studies in E. coli support a role for direct transport of glucose-6-phosphate from the host cell by the parasite transporters CpGT1 and CpGT2, thus bypassing a requirement for hexokinase. Additionally, the parasite obtains phosphorylated glucose from amylopectin stores that are released by the action of the essential enzyme glycogen phosphorylase. Collectively, these findings reveal that C. parvum relies on multiple pathways to obtain phosphorylated glucose both for glycolysis and to restore carbohydrate reserves.

Morphological and molecular phylogenetic characterization of Sarcocystis kani sp. nov. and other novel, closely related Sarcocystis spp. infecting small mammals and colubrid snakes in Asia

We investigated the morphology and phylogenetic relationships of novel and previously recognized Sarcocystis spp. infecting small mammals and colubrid snakes in Asia. The nuclear 18S rRNA and mitochondrial cox1 of Sarcocystis sp.1 from mangrove snakes (Boiga dendrophila) in Thailand and Sarcocystis sp.2 from a ricefield rat (Rattus argentiventer) in Sumatra were partially sequenced. Sporocysts of Sarcocystis sp.1 induced development of sarcocysts in experimentally infected rats, which showed a unique ultrastructure that was observed previously by S.P. Kan in rats from Malaysia; therefore, we describe this species as Sarcocystis kani sp. nov. Its integration into the 18S rRNA phylogeny of Sarcocystis spp. cycling between small mammals and colubrid snakes helped clarify relationships among the so-called S. zuoi-complex of molecularly cryptic species: Sarcocystis kani sp. nov., S. sp.2, S. attenuati, S. scandentiborneensis, and S. zuoi were all included in this clade. Tree topology was resolved into dichotomies congruent with the morphological disparities between the taxa. However, cox1 gene sequencing (including newly sequenced S. singaporensis and S. zamani) revealed that Sarcocystis kani, S. attenuati, and S. scandentiborneensis were identical suggesting a recent, common ancestry. To identify other distinctive features, lineage-specific molecular patterns within both genes were examined revealing that all 18S rRNA sequences of the S. zuoi - complex possess a unique, 7-nt long motif in helix 38 of domain V7 that was different in S. clethrionomyelaphis which branched off basally from the complex. Three-dimensional homology modelling of COX1 protein structure identified amino acid substitutions within the barcode area specific for the S. zuoi-complex and substantial divergence in structurally important amino acids between Sarcocystis species of snakes as definitive hosts and other lineages of the Sarcocystidae. We discuss the utility of selected genes for species delimitation of the Sarcocystis spp. under investigation, which probably evolved during recent radiations of their intermediate and definitive hosts.

The multifaceted roles of Myb domain-containing proteins in apicomplexan parasites

Apicomplexan parasites are a large and diverse clade of protists responsible for significant diseases of humans and animals. Central to the ability of these parasites to colonize their host and evade immune responses is an expanded repertoire of gene-expression programs that requires the coordinated action of complex transcriptional networks. DNA-binding proteins and chromatin regulators are essential orchestrators of apicomplexan gene expression that often act in concert. Although apicomplexan genomes encode various families of putative DNA-binding proteins, most remain functionally and mechanistically unexplored. This review highlights the versatile role of myeloblastosis (Myb) domain-containing proteins in apicomplexan parasites as transcription factors and chromatin regulators. We explore the diversity of Myb domain structure and use phylogenetic analysis to identify common features across the phylum. This provides a framework to discuss functional heterogeneity and regulation of Myb domain-containing proteins particularly emphasizing their role in parasite differentiation.

The evolution of endoparasitism and complex life cycles in parasitic platyhelminths

Within flatworms, the vast majority of parasitism is innate to Neodermata, the most derived and diversified group of the phylum Platyhelminthes.1,2 The four major lineages of Neodermata maintain various combinations of life strategies.3 They include both externally (ecto-) and internally feeding (endo-) parasites. Some lineages complete their life cycles directly by infecting a single host, whereas others succeed only through serial infections of multiple hosts of various vertebrate and invertebrate groups. Food sources and modes of digestion add further combinatorial layers to the often incompletely understood mosaic of neodermatan life histories. Their evolutionary trajectories have remained molecularly unresolved because of conflicting evolutionary inferences and a lack of genomic data.4 Here, we generated transcriptomes for nine early branching neodermatan representatives and performed detailed phylogenomic analyses to address these critical gaps. Polyopisthocotylea, mostly hematophagous ectoparasites, form a group with the mostly hematophagous but endoparasitic trematodes (Trematoda), rather than sharing a common ancestor with Monopisthocotylea, ectoparasitic epithelial feeders. Phylogenetic placement of the highly specialized endoparasitic Cestoda alters depending on the model. Regardless of this uncertainty, this study brings an unconventional perspective on the evolution of platyhelminth parasitism, rejecting a common origin for the endoparasitic lifestyle intrinsic to cestodes and trematodes. Instead, our data indicate that complex life cycles and invasion of vertebrates' gut lumen, the hallmark features of these parasites, evolved independently within Neodermata. We propose the demise of the traditionally recognized class Monogenea and the promotion of its two subclasses to the class level as Monopisthocotyla new class and Polyopisthocotyla new class.

Microbiota-produced indole metabolites disrupt mitochondrial function and inhibit Cryptosporidium parvum growth

Cryptosporidiosis is a leading cause of life-threatening diarrhea in young children in resource-poor settings. To explore microbial influences on susceptibility, we screened 85 microbiota-associated metabolites for their effects on Cryptosporidium parvum growth in vitro. We identify eight inhibitory metabolites in three main classes: secondary bile salts/acids, a vitamin B6 precursor, and indoles. Growth restriction of C. parvum by indoles does not depend on the host aryl hydrocarbon receptor (AhR) pathway. Instead, treatment impairs host mitochondrial function and reduces total cellular ATP, as well as directly reducing the membrane potential in the parasite mitosome, a degenerate mitochondria. Oral administration of indoles, or reconstitution of the gut microbiota with indole-producing bacteria, delays life cycle progression of the parasite in vitro and reduces the severity of C. parvum infection in mice. Collectively, these findings indicate that microbiota metabolites impair mitochondrial function and contribute to colonization resistance to Cryptosporidium infection.

A Large-Scale Study into Protist-Animal Interactions Based on Public Genomic Data Using DNA Barcodes

With the birth of next-generation sequencing (NGS) technology, genomic data in public databases have increased exponentially. Unfortunately, exogenous contamination or intracellular parasite sequences in assemblies could confuse genomic analysis. Meanwhile, they can provide a valuable resource for studies of host-microbe interactions. Here, we used a strategy based on DNA barcodes to scan protistan contamination in the GenBank WGS/TSA database. The results showed a total of 13,952 metazoan/animal assemblies in GenBank, where 17,036 contigs were found to be protistan contaminants in 1507 assemblies (10.8%), with even higher contamination rates in taxa of Cnidaria (150/281), Crustacea (237/480), and Mollusca (107/410). Taxonomic analysis of the protists derived from these contigs showed variations in abundance and evenness of protistan contamination across different metazoan taxa, reflecting host preferences of Apicomplexa, Ciliophora, Oomycota and Symbiodiniaceae for mammals and birds, Crustacea, insects, and Cnidaria, respectively. Finally, mitochondrial proteins COX1 and CYTB were predicted from these contigs, and the phylogenetic analysis corroborated the protistan origination and heterogeneous distribution of the contaminated contigs. Overall, in this study, we conducted a large-scale scan of protistan contaminant in genomic resources, and the protistan sequences detected will help uncover the protist diversity and relationships of these picoeukaryotes with Metazoa.

NAD(P) transhydrogenase isoform distribution provides insight into apicomplexan evolution

Membrane-located NAD(P) transhydrogenase (NTH) catalyses reversible hydride ion transfer between NAD(H) and NADP(H), simultaneously translocating a proton across the membrane. The enzyme is structurally conserved across prokaryotes and eukaryotes. In heterotrophic bacteria NTH proteins reside in the cytoplasmic membrane, whereas in animals they localise in the mitochondrial inner membrane. Eukaryotic NTH proteins exists in two distinct configurations (isoforms) and have non-mitochondrial functions in unicellular eukaryotes like Plasmodium, the causative agent of malaria. In this study, we carried out a systematic analysis of nth genes across eukaryotic life to determine its prevalence and distribution of isoforms. The results reveal that NTH is found across all major lineages, but that some organisms, notably plants, lack nth genes altogether. Isoform distribution and phylogenetic analysis reveals different nth gene loss scenarios in apicomplexan lineages, which sheds new light on the evolution of the Piroplasmida and Eimeriidae.

Multiple pathways for glucose phosphate transport and utilization support growth of Cryptosporidium parvum

Cryptosporidium parvum is an obligate intracellular parasite with a highly reduced mitochondrion that lacks the TCA cycle and the ability to generate ATP, making the parasite reliant on glycolysis. Genetic ablation experiments demonstrated that neither of the two putative glucose transporters CpGT1 and CpGT2 were essential for growth. Surprisingly, hexokinase was also dispensable for parasite growth while the downstream enzyme aldolase was required, suggesting the parasite has an alternative way of obtaining phosphorylated hexose. Complementation studies in E. coli support a role for direct transport of glucose-6-phosphate from the host cell by the parasite transporters CpGT1 and CpGT2, thus bypassing a requirement for hexokinase. Additionally, the parasite obtains phosphorylated glucose from amylopectin stores that are released by the action of the essential enzyme glycogen phosphorylase. Collectively, these findings reveal that C. parvum relies on multiple pathways to obtain phosphorylated glucose both for glycolysis and to restore carbohydrate reserves.

Microbiota produced indole metabolites disrupt host cell mitochondrial energy production and inhibit Cryptosporidium parvum growth

Cryptosporidiosis is a leading cause of life-threatening diarrhea in young children in resource-poor settings. Susceptibility rapidly declines with age, associated with changes in the microbiota. To explore microbial influences on susceptibility, we screened 85 microbiota- associated metabolites enriched in the adult gut for their effects on C. parvum growth in vitro. We identified eight inhibitory metabolites in three main classes: secondary bile salts/acids, a vitamin B 6 precursor, and indoles. Growth restriction of C. parvum by indoles did not depend on the host aryl hydrocarbon receptor (AhR) pathway. Instead, treatment impaired host mitochondrial function and reduced total cellular ATP, as well as directly reduced the membrane potential in the parasite mitosome, a degenerate mitochondria. Oral administration of indoles, or reconstitution of the gut microbiota with indole producing bacteria, delayed life cycle progression of the parasite in vitro and reduced severity of C. parvum infection in mice. Collectively, these findings indicate that microbiota metabolites contribute to colonization resistance to Cryptosporidium infection.

The Tubulin Superfamily in Apicomplexan Parasites

Microtubules and specialized microtubule-containing structures are assembled from tubulins, an ancient superfamily of essential eukaryotic proteins. Here, we use bioinformatic approaches to analyze features of tubulins in organisms from the phylum Apicomplexa. Apicomplexans are protozoan parasites that cause a variety of human and animal infectious diseases. Individual species harbor one to four genes each for α- and β-tubulin isotypes. These may specify highly similar proteins, suggesting functional redundancy, or exhibit key differences, consistent with specialized roles. Some, but not all apicomplexans harbor genes for δ- and ε-tubulins, which are found in organisms that construct appendage-containing basal bodies. Critical roles for apicomplexan δ- and ε-tubulin are likely to be limited to microgametes, consistent with a restricted requirement for flagella in a single developmental stage. Sequence divergence or the loss of δ- and ε-tubulin genes in other apicomplexans appears to be associated with diminished requirements for centrioles, basal bodies, and axonemes. Finally, because spindle microtubules and flagellar structures have been proposed as targets for anti-parasitic therapies and transmission-blocking strategies, we discuss these ideas in the context of tubulin-based structures and tubulin superfamily properties.

A Mitosome With Distinct Metabolism in the Uncultured Protist Parasite Paramikrocytos canceri (Rhizaria, Ascetosporea)

Ascetosporea are endoparasites of marine invertebrates that include economically important pathogens of aquaculture species. Owing to their often-minuscule cell sizes, strict intracellular lifestyle, lack of cultured representatives and minimal availability of molecular data, these unicellular parasites remain poorly studied. Here, we sequenced and assembled the genome and transcriptome of Paramikrocytos canceri, an endoparasite isolated from the European edible crab Cancer pagurus. Using bioinformatic predictions, we show that P. canceri likely possesses a mitochondrion-related organelle (MRO) with highly reduced metabolism, resembling the mitosomes of other parasites but with key differences. Like other mitosomes, this MRO is predicted to have reduced metabolic capacity and lack an organellar genome and function in iron-sulfur cluster (ISC) pathway-mediated Fe-S cluster biosynthesis. However, the MRO in P. canceri is uniquely predicted to produce ATP via a partial glycolytic pathway and synthesize phospholipids de novo through the CDP-DAG pathway. Heterologous gene expression confirmed that proteins from the ISC and CDP-DAG pathways retain mitochondrial targeting sequences that are recognized by yeast mitochondria. This represents a unique combination of metabolic pathways in an MRO, including the first reported case of a mitosome-like organelle able to synthesize phospholipids de novo. Some of these phospholipids, such as phosphatidylserine, are vital in other protist endoparasites that invade their host through apoptotic mimicry.

Multilocus Sequence Typing as a Useful Tool for the Study of the Genetic Diversity and Population Structure of Cryptosporidium Spp

One of the most important aquatic parasites in industrialized countries, Cryptosporidium spp., is a major cause of diarrheal disease in humans and animals worldwide. The contingent evolution of cryptosporidia with hosts, host adaptation, and geographic variation contributed to the creation of species subtypes, thereby shaping their population genetic structures. Multilocus typing tools for population genetic characterizations of transmission dynamics and delineation of mechanisms for the emergence of virulent subtypes have played an important role in improving our understanding of the transmission of this parasite. However, to better understand the significance of different subtypes with clinical disease manifestations and transmission risks, a large number of samples and preferably from different geographical areas need to be analyzed. This review provides an analysis of genetic variation through multilocus sequence typing, provides an overview of subtypes, typing gene markers for Cryptosporidium parvum, Cryptosporidium hominis, Cryptosporidium muris and Cryptosporidium andersoni genotypes and an overview of the hosts of these parasites.

Reconstruction of Plastid Proteomes of Apicomplexans and Close Relatives Reveals the Major Evolutionary Outcomes of Cryptic Plastids

Apicomplexans and related lineages comprise many obligate symbionts of animals; some of which cause notorious diseases such as malaria. They evolved from photosynthetic ancestors and transitioned into a symbiotic lifestyle several times, giving rise to species with diverse non-photosynthetic plastids. Here, we sought to reconstruct the evolution of the cryptic plastids in the apicomplexans, chrompodellids, and squirmids (ACS clade) by generating five new single-cell transcriptomes from understudied gregarine lineages, constructing a robust phylogenomic tree incorporating all ACS clade sequencing datasets available, and using these to examine in detail, the evolutionary distribution of all 162 proteins recently shown to be in the apicoplast by spatial proteomics in Toxoplasma. This expanded homology-based reconstruction of plastid proteins found in the ACS clade confirms earlier work showing convergence in the overall metabolic pathways retained once photosynthesis is lost, but also reveals differences in the degrees of plastid reduction in specific lineages. We show that the loss of the plastid genome is common and unexpectedly find many lineage- and species-specific plastid proteins, suggesting the presence of evolutionary innovations and neofunctionalizations that may confer new functional and metabolic capabilities that are yet to be discovered in these enigmatic organelles.

Reduced mitochondria provide an essential function for the cytosolic methionine cycle

The loss of mitochondria in oxymonad protists has been associated with the redirection of the essential Fe-S cluster assembly to the cytosol. Yet as our knowledge of diverse free-living protists broadens, the list of functions of their mitochondrial-related organelles (MROs) expands. We revealed another such function in the closest oxymonad relative, Paratrimastix pyriformis, after we solved the proteome of its MRO with high accuracy, using localization of organelle proteins by isotope tagging (LOPIT). The newly assigned enzymes connect to the glycine cleavage system (GCS) and produce folate derivatives with one-carbon units and formate. These are likely to be used by the cytosolic methionine cycle involved in S-adenosyl methionine recycling. The data provide consistency with the presence of the GCS in MROs of free-living species and its absence in most endobionts, which typically lose the methionine cycle and, in the case of oxymonads, the mitochondria.

Marine gregarine genomes reveal the breadth of apicomplexan diversity with a partially conserved glideosome machinery

Our current view of the evolutionary history, coding and adaptive capacities of Apicomplexa, protozoan parasites of a wide range of metazoan, is currently strongly biased toward species infecting humans, as data on early diverging apicomplexan lineages infecting invertebrates is extremely limited. Here, we characterized the genome of the marine eugregarine Porospora gigantea, intestinal parasite of Lobsters, remarkable for the macroscopic size of its vegetative feeding forms (trophozoites) and its gliding speed, the fastest so far recorded for Apicomplexa. Two highly syntenic genomes named A and B were assembled. Similar in size (~ 9 Mb) and coding capacity (~ 5300 genes), A and B genomes are 10.8% divergent at the nucleotide level, corresponding to 16-38 My in divergent time. Orthogroup analysis across 25 (proto)Apicomplexa species, including Gregarina niphandrodes, showed that A and B are highly divergent from all other known apicomplexan species, revealing an unexpected breadth of diversity. Phylogenetically these two species branch sisters to Cephaloidophoroidea, and thus expand the known crustacean gregarine superfamily. The genomes were mined for genes encoding proteins necessary for gliding, a key feature of apicomplexans parasites, currently studied through the molecular model called glideosome. Sequence analysis shows that actin-related proteins and regulatory factors are strongly conserved within apicomplexans. In contrast, the predicted protein sequences of core glideosome proteins and adhesion proteins are highly variable among apicomplexan lineages, especially in gregarines. These results confirm the importance of studying gregarines to widen our biological and evolutionary view of apicomplexan species diversity, and to deepen our understanding of the molecular bases of key functions such as gliding, well known to allow access to the intracellular parasitic lifestyle in Apicomplexa.

Protists in the Insect Rearing Industry: Benign Passengers or Potential Risk?

As the insects for food and feed industry grows, a new understanding of the industrially reared insect microbiome is needed to better comprehend the role that it plays in both maintaining insect health and generating disease. While many microbiome projects focus on bacteria, fungi or viruses, protists (including microsporidia) can also make up an important part of these assemblages. Past experiences with intensive invertebrate rearing indicate that these parasites, whilst often benign, can rapidly sweep through populations, causing extensive damage. Here, we review the diversity of microsporidia and protist species that are found in reared insect hosts and describe the current understanding of their host spectra, life cycles and the nature of their interactions with hosts. Major entomopathogenic parasite groups with the potential to infect insects currently being reared for food and feed include the Amoebozoa, Apicomplexa, Ciliates, Chlorophyta, Euglenozoa, Ichtyosporea and Microsporidia. However, key gaps exist in the understanding of how many of these entomopathogens affect host biology. In addition, for many of them, there are very limited or even no molecular data, preventing the implementation of molecular detection methods. There is now a pressing need to develop and use novel molecular tools, coupled with standard molecular diagnostic methods, to help unlock their biology and predict the effects of these poorly studied protist parasites in intensive insect rearing systems.

The Elusive Mitochondrial Genomes of Apicomplexa: Where Are We Now?

Mitochondria are vital organelles of eukaryotic cells, participating in key metabolic pathways such as cellular respiration, thermogenesis, maintenance of cellular redox potential, calcium homeostasis, cell signaling, and cell death. The phylum Apicomplexa is entirely composed of obligate intracellular parasites, causing a plethora of severe diseases in humans, wild and domestic animals. These pathogens include the causative agents of malaria, cryptosporidiosis, neosporosis, East Coast fever and toxoplasmosis, among others. The mitochondria in Apicomplexa has been put forward as a promising source of undiscovered drug targets, and it has been validated as the target of atovaquone, a drug currently used in the clinic to counter malaria. Apicomplexans present a single tubular mitochondria that varies widely both in structure and in genomic content across the phylum. The organelle is characterized by massive gene migrations to the nucleus, sequence rearrangements and drastic functional reductions in some species. Recent third generation sequencing studies have reignited an interest for elucidating the extensive diversity displayed by the mitochondrial genomes of apicomplexans and their intriguing genomic features. The underlying mechanisms of gene transcription and translation are also ill-understood. In this review, we present the state of the art on mitochondrial genome structure, composition and organization in the apicomplexan phylum revisiting topological and biochemical information gathered through classical techniques. We contextualize this in light of the genomic insight gained by second and, more recently, third generation sequencing technologies. We discuss the mitochondrial genomic and mechanistic features found in evolutionarily related alveolates, and discuss the common and distinct origins of the apicomplexan mitochondria peculiarities.

Taxonomy and molecular epidemiology of Cryptosporidium and Giardia - a 50 year perspective (1971-2021)

The protozoan parasites Cryptosporidium and Giardia are significant causes of diarrhoea worldwide and are responsible for numerous waterborne and foodborne outbreaks of diseases. Over the last 50 years, the development of improved detection and typing tools has facilitated the expanding range of named species. Currently at least 44 Cryptosporidium spp. and >120 genotypes, and nine Giardia spp., are recognised. Many of these Cryptosporidium genotypes will likely be described as species in the future. The phylogenetic placement of Cryptosporidium at the genus level is still unclear and further research is required to better understand its evolutionary origins. Zoonotic transmission has long been known to play an important role in the epidemiology of cryptosporidiosis and giardiasis, and the development and application of next generation sequencing tools is providing evidence for this. Comparative whole genome sequencing is also providing key information on the genetic mechanisms for host specificity and human infectivity, and will enable One Health management of these zoonotic parasites in the future.

PhyloFisher: A phylogenomic package for resolving eukaryotic relationships

Phylogenomic analyses of hundreds of protein-coding genes aimed at resolving phylogenetic relationships is now a common practice. However, no software currently exists that includes tools for dataset construction and subsequent analysis with diverse validation strategies to assess robustness. Furthermore, there are no publicly available high-quality curated databases designed to assess deep (>100 million years) relationships in the tree of eukaryotes. To address these issues, we developed an easy-to-use software package, PhyloFisher (https://github.com/TheBrownLab/PhyloFisher), written in Python 3. PhyloFisher includes a manually curated database of 240 protein-coding genes from 304 eukaryotic taxa covering known eukaryotic diversity, a novel tool for ortholog selection, and utilities that will perform diverse analyses required by state-of-the-art phylogenomic investigations. Through phylogenetic reconstructions of the tree of eukaryotes and of the Saccharomycetaceae clade of budding yeasts, we demonstrate the utility of the PhyloFisher workflow and the provided starting database to address phylogenetic questions across a large range of evolutionary time points for diverse groups of organisms. We also demonstrate that undetected paralogy can remain in phylogenomic "single-copy orthogroup" datasets constructed using widely accepted methods such as all vs. all BLAST searches followed by Markov Cluster Algorithm (MCL) clustering and application of automated tree pruning algorithms. Finally, we show how the PhyloFisher workflow helps detect inadvertent paralog inclusions, allowing the user to make more informed decisions regarding orthology assignments, leading to a more accurate final dataset.

Nutrient Acquisition and Attachment Strategies in Basal Lineages: A Tough Nut to Crack in the Evolutionary Puzzle of Apicomplexa

Apicomplexa are unicellular eukaryotes that parasitise a wide spectrum of invertebrates and vertebrates, including humans. In their hosts, they occupy a variety of niches, from extracellular cavities (intestine, coelom) to epicellular and intracellular locations, depending on the species and/or developmental stages. During their evolution, Apicomplexa thus developed an exceptionally wide range of unique features to reach these diversified parasitic niches and to survive there, at least long enough to ensure their own transmission or that of their progeny. This review summarises the current state of knowledge on the attachment/invasive and nutrient uptake strategies displayed by apicomplexan parasites, focusing on trophozoite stages of their so far poorly studied basal representatives, which mostly parasitise invertebrate hosts. We describe their most important morphofunctional features, and where applicable, discuss existing major similarities and/or differences in the corresponding mechanisms, incomparably better described at the molecular level in the more advanced Apicomplexa species, of medical and veterinary significance, which mainly occupy intracellular niches in vertebrate hosts.