Phylogenomics Identifies a New Major Subgroup of Apicomplexans, Marosporida class nov., with Extreme Apicoplast Genome Reduction

Survival of the fittest: genomic investigations of the bay scallop reveal a shift in population structure through a summer mortality event

BACKGROUND

Understanding the genetic basis of resilience in marine organisms is critical for conservation and management, particularly in the face of escalating environmental stress and disease outbreaks. The bay scallop Argopecten irradians is a commercially and recreationally important shellfish species found in estuarine and coastal environments of the United States from New England to the Gulf of Mexico. In New York, adult bay scallop populations have been decimated every summer since 2019 leading to the collapse of their fishery. These mortality events were associated with annual outbreaks of an undescribed apicomplexan parasite recently named Bay Scallop Marosporida (BSM) that disrupts scallop kidneys.

RESULTS

This study investigates host-pathogen interactions and assesses changes in population structure during BSM-associated mortality events. The research compared wild and aquacultured scallops used for stock enhancement in New York, revealing significant change in population structures throughout the mortality outbreak. The results underscore the selective pressures exerted by BSM infection and environmental stressors, as evidenced by shifts in genetic divergence and allele frequencies particularly in genes associated with kidney function, stress and infection response. Through a detailed genomic and population genetic approach, this research represents a unique case study highlighting the impact of disease on marine biodiversity and advances our understanding of the impact of summer mortality events on the scallop population in NY.

CONCLUSIONS

This study highlights changes in the genomic structure of bay scallops during a BSM-associated mortality event. Identified mutations (such as the one in the nephrocystin-3-like gene) represent prime candidates for specific targeted investigations to link genotypes to phenotypes. By integrating genomic and epidemiological data, the research provides a basis for understanding the impact of disease on scallop biodiversity. These findings may help guide conservation strategies for sustainable fisheries in the face of environmental change and disease outbreaks.

Adaptations and metabolic evolution of myzozoan protists across diverse lifestyles and environments

SUMMARYMyzozoans encompass apicomplexans and dinoflagellates that manifest diverse lifestyles in highly varied environments. They show enormous propensity to employ different metabolic programs and exploit different nutrient resources and niches, and yet, they share much core biology that underlies this evolutionary success and impact. This review discusses apicomplexan parasites of medical significance and the traits and properties they share with non-pathogenic myzozoans. These include the versatility of myzozoan plastids, which scale from fully photosynthetic organelles to the site of very select key metabolic pathways. Pivotal evolutionary innovations, such as the apical complex, have allowed myzozoans to shift from predatory to parasitic and other symbiotic lifestyles multiple times in both apicomplexan and dinoflagellate branches of the myzozoan evolutionary tree. Such traits, along with shared mechanisms for nutrient acquisition, appear to underpin the prosperity of myzozoans in their varied habitats. Understanding the mechanisms of these shared traits has the potential to spawn new strategic interventions against medically and veterinary relevant parasites within this grouping.

The apicoplast biogenesis and metabolism: current progress and questions

Many apicomplexan parasites have a chloroplast-derived apicoplast containing several metabolic pathways. Recent studies have greatly expanded our understanding of apicoplast biogenesis and metabolism while also raising new questions. Here, we review recent progress on the biological roles of individual metabolic pathways, focusing on two medically important parasites, Plasmodium spp. and Toxoplasma gondii. We highlight the similarities and differences in how similar apicoplast metabolic pathways are utilized to adapt to different parasitic lifestyles. The execution of apicoplast metabolic functions requires extensive interactions with other subcellular compartments, but the underlying mechanisms remain largely unknown. Apicoplast metabolic functions have historically been considered attractive drug targets, and a comprehensive understanding of their metabolic capacities and interactions with other organelles is essential to fully realize their potential.

Chromosome-level genome assembly of the bay scallop Argopecten irradians

The bay scallop, Argopecten irradians, is a species of major commercial, cultural, and ecological importance. It is endemic to the eastern coast of the United States, but has also been introduced to China, where it supports a significant aquaculture industry. Here, we provide an annotated chromosome-level reference genome assembly for the bay scallop, assembled using PacBio and Hi-C data. The total genome size is 845.9 Mb, distributed over 1,503 scaffolds with a scaffold N50 of 44.3 Mb. The majority (92.9%) of the assembled genome is contained within the 16 largest scaffolds, corresponding to the 16 chromosomes confirmed by Hi-C analysis. The assembly also includes the complete mitochondrial genome. Approximately 36.2% of the genome consists of repetitive elements. The BUSCO analysis showed a completeness of 96.2%. We identified 33,772 protein-coding genes. This genome assembly will be a valuable resource for future research on evolutionary dynamics, adaptive mechanisms, and will support genome-assisted breeding, contributing to the conservation and management of this iconic species in the face of environmental and pathogenic challenges.

Genome evolution in intracellular parasites: Microsporidia and Apicomplexa

Microsporidia and Apicomplexa are eukaryotic, single-celled, intracellular parasites with huge public health and economic importance. Typically, these parasites are studied separately, emphasizing their uniqueness and diversity. In this review, we explore the huge amount of genomic data that has recently become available for the two groups. We compare and contrast their genome evolution and discuss how their transitions to intracellular life may have shaped it. In particular, we explore genome reduction and compaction, genome expansion and ploidy, gene shuffling and rearrangements, and the evolution of centromeres and telomeres.

Towards high-throughput parallel imaging and single-cell transcriptomics of microbial eukaryotic plankton

Single-cell transcriptomics has the potential to provide novel insights into poorly studied microbial eukaryotes. Although several such technologies are available and benchmarked on mammalian cells, few have been tested on protists. Here, we applied a microarray single-cell sequencing (MASC-seq) technology, that generates microscope images of cells in parallel with capturing their transcriptomes, on three species representing important plankton groups with different cell structures; the ciliate Tetrahymena thermophila, the diatom Phaeodactylum tricornutum, and the dinoflagellate Heterocapsa sp. Both the cell fixation and permeabilization steps were adjusted. For the ciliate and dinoflagellate, the number of transcripts of microarray spots with single cells were significantly higher than for background spots, and the overall expression patterns were correlated with that of bulk RNA, while for the much smaller diatom cells, it was not possible to separate single-cell transcripts from background. The MASC-seq method holds promise for investigating "microbial dark matter", although further optimizations are necessary to increase the signal-to-noise ratio.

Patterns of host-parasite associations between marine meiofaunal flatworms (Platyhelminthes) and rhytidocystids (Apicomplexa)

Microturbellarians are abundant and ubiquitous members of marine meiofaunal communities around the world. Because of their small body size, these microscopic animals are rarely considered as hosts for parasitic organisms. Indeed, many protists, both free-living and parasitic ones, equal or surpass meiofaunal animals in size. Despite several anecdotal records of "gregarines", "sporozoans", and "apicomplexans" parasitizing microturbellarians in the literature-some of them dating back to the nineteenth century-these single-celled parasites have never been identified and characterized. More recently, the sequencing of eukaryotic microbiomes in microscopic invertebrates have revealed a hidden diversity of protist parasites infecting microturbellarians and other meiofaunal animals. Here we show that apicomplexans isolated from twelve taxonomically diverse rhabdocoel taxa and one species of proseriate collected in four geographically distinct areas around the Pacific Ocean (Okinawa, Hokkaido, and British Columbia) and the Caribbean Sea (Curaçao) all belong to the apicomplexan genus Rhytidocystis. Based on comprehensive molecular phylogenies of Rhabdocoela and Proseriata inferred from both 18S and 28S rDNA sequences, as well as a molecular phylogeny of Marosporida inferred from 18S rDNA sequences, we determine the phylogenetic positions of the microturbellarian hosts and their parasites. Multiple lines of evidence, including morphological and molecular data, show that at least nine new species of Rhytidocystis infect the microturbellarian hosts collected in this study, more than doubling the number of previously recognized species of Rhytidocystis, all of which infect polychaete hosts. A cophylogenetic analysis examining patterns of phylosymbiosis between hosts and parasites suggests a complex picture of overall incongruence between host and parasite phylogenies, and varying degrees of geographic signals and taxon specificity.

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.

p25alpha Domain-Containing Proteins of Apicomplexans and Related Taxa

TPPP (tubulin polymerization promoting protein)-like proteins contain one or more p25alpha (Pfam05517) domains. TPPP-like proteins occur in different types as determined by their length (e.g., long-, short-, truncated-, and fungal-type TPPP) and include the protein apicortin, which possesses another domain, doublecortin (DCX, Pfam 03607). These various TPPP-like proteins are found in various phylogenomic groups. In particular, short-type TPPPs and apicortin are well-represented in the Myzozoa, which include apicomplexans and related taxa, chrompodellids, dinoflagellates, and perkinsids. The long-, truncated-, and fungal-type TPPPs are not found in the myzozoans. Apicortins are found in all apicomplexans except one piroplasmid species, present in several other myzozoans, and seem to be correlated with the conoid and apical complex. Short-type TPPPs are predominantly found in myzozoans that have flagella, suggesting a role in flagellum assembly or structure.

Characterisation of the OTU domain deubiquitinase complement of Toxoplasma gondii

The phylum Apicomplexa contains several parasitic species of medical and agricultural importance. The ubiquitination machinery remains, for the most part, uncharacterised in apicomplexan parasites, despite the important roles that it plays in eukaryotic biology. Bioinformatic analysis of the ubiquitination machinery in apicomplexan parasites revealed an expanded ovarian tumour domain-containing (OTU) deubiquitinase (DUB) family in Toxoplasma, potentially reflecting functional importance in apicomplexan parasites. This study presents comprehensive characterisation of Toxoplasma OTU DUBs. AlphaFold-guided structural analysis not only confirmed functional orthologues found across eukaryotes, but also identified apicomplexan-specific enzymes, subsequently enabling discovery of a cryptic OTU DUB in Plasmodium species. Comprehensive biochemical characterisation of 11 Toxoplasma OTU DUBs revealed activity against ubiquitin- and NEDD8-based substrates and revealed ubiquitin linkage preferences for Lys6-, Lys11-, Lys48-, and Lys63-linked chain types. We show that accessory domains in Toxoplasma OTU DUBs impose linkage preferences, and in case of apicomplexan-specific TgOTU9, we discover a cryptic ubiquitin-binding domain that is essential for TgOTU9 activity. Using the auxin-inducible degron (AID) to generate knockdown parasite lines, TgOTUD6B was found to be important for Toxoplasma growth.

Genome sequence of Ophryocystis elektroscirrha, an apicomplexan parasite of monarch butterflies: cryptic diversity and response to host-sequestered plant chemicals

Apicomplexa are ancient and diverse organisms which have been poorly characterized by modern genomics. To better understand the evolution and diversity of these single-celled eukaryotes, we sequenced the genome of Ophryocystis elektroscirrha, a parasite of monarch butterflies, Danaus plexippus. We contextualize our newly generated resources within apicomplexan genomics before answering longstanding questions specific to this host-parasite system. To start, the genome is miniscule, totaling only 9 million bases and containing fewer than 3,000 genes, half the gene content of two other sequenced invertebrate-infecting apicomplexans, Porospora gigantea and Gregarina niphandrodes. We found that O. elektroscirrha shares different orthologs with each sequenced relative, suggesting the true set of universally conserved apicomplexan genes is very small indeed. Next, we show that sequencing data from other potential host butterflies can be used to diagnose infection status as well as to study diversity of parasite sequences. We recovered a similarly sized parasite genome from another butterfly, Danaus chrysippus, that was highly diverged from the O. elektroscirrha reference, possibly representing a distinct species. Using these two new genomes, we investigated potential evolutionary response by parasites to toxic phytochemicals their hosts ingest and sequester. Monarch butterflies are well-known to tolerate toxic cardenolides thanks to changes in the sequence of their Type II ATPase sodium pumps. We show that Ophryocystis completely lacks Type II or Type 4 sodium pumps, and related proteins PMCA calcium pumps show extreme sequence divergence compared to other Apicomplexa, demonstrating new avenues of research opened by genome sequencing of non-model Apicomplexa.

An apicomplexan parasite drives the collapse of the bay scallop population in New York

The bay scallop, Argopecten irradians, represents a commercially, culturally and ecologically important species found along the United States' Atlantic and Gulf coasts. Since 2019, scallop populations in New York have been suffering large-scale summer mortalities resulting in 90-99% reduction in biomass of adult scallops. Preliminary investigations of these mortality events showed 100% prevalence of an apicomplexan parasite infecting kidney tissues. This study was designed to provide histological, ultrastructural and molecular characteristics of a non-described parasite, member of the newly established Marosporida clade (Apicomplexa) and provisionally named BSM (Bay Scallop Marosporida). Molecular diagnostics tools (quantitative PCR, in situ hybridization) were developed and used to monitor disease development. Results showed that BSM disrupts multiple scallop tissues including kidney, adductor muscle, gill, and gonad. Microscopy observations allowed the identification of both intracellular and extracellular stages of the parasite. Field surveys demonstrated a strong seasonal signature in disease prevalence and intensity, as severe cases and mortality increase as summer progresses. These results strongly suggest that BSM infection plays a major role in the collapse of bay scallop populations in New York. In this framework, BSM may synergistically interact with stressful environmental conditions to impair the host and lead to mortality.

Morphological and molecular analysis of Aggregata aspera n. sp. (Apicomplexa: Aggregatidae) in Amphioctopus ovulum and Amphioctopus marginatus (Mollusca: Cephalopoda) from the Western Pacific Ocean

Aggregata Frenzel, 1885 (Apicomplexa) are dangerous protozoan parasites that cause malabsorption syndrome in wild and reared cephalopod species, resulting in significant economic loss to fishery and aquaculture industries. The new parasitic species, Aggregata aspera n. sp., in the digestive tract of Amphioctopus ovulum and Amphioctopus marginatus from an area in the Western Pacific Ocean was identified, it is the second two-host parasite species of Aggregata. Mature oocysts and sporocysts were spherical to ovoid in shape. Sporulated oocysts were 380.6-1,158.4 μm in length and 284.0-1,090.6 μm in width. The mature sporocysts were 16.2-18.3 μm in length and 15.7-17.6 μm in width, with irregular protuberances on the lateral wall of the sporocysts. Sporozoites within mature sporocysts were curled in shape and measured 13.0-17.0 μm in length and 1.6-2.4 μm in width. Each sporocyst contained 12-16 sporozoites. Phylogenetic tree analysis, based on 18S rRNA gene partial sequences, indicated that Ag. aspera forms a monophyletic cluster within the genus Aggregata and has a sister relationship with Ag. sinensis. These findings will provide the theoretical basis for the histopathology and diagnosis of coccidiosis in cephalopods.

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.

Coinfection patterns of two marine apicomplexans are not associated with genetic diversity of their polychaete host

Coinfections of two or more parasites within one host are more of a rule than an exception in nature. Interactions between coinfecting parasites can greatly affect their abundance and prevalence. Characteristics of the host, such as genetic diversity, can also affect the infection dynamics of coinfecting parasites. Here, we investigate for the first time the association of coinfection patterns of two marine apicomplexans, Rhytidocystis sp. and Selenidium pygospionis, with the genetic diversity of their host, the polychaete Pygospio elegans, from natural populations. Host genetic diversity was determined with seven microsatellite loci and summarized as allelic richness, inbreeding coefficient, and individual heterozygosity. We detected nonsignificant correlations between infection loads and both individual host heterozygosity and population genetic diversity. Prevalence and infection load of Rhytidocystis sp. were higher than those of S. pygospionis, and both varied spatially. Coinfections were common, and almost all hosts infected by S. pygospionis were also infected by Rhytidocystis sp. Rhytidocystis sp. infection load was significantly higher in dual infections. Our results suggest that factors other than host genetic diversity might be more important in marine apicomplexan infection patterns and experimental approaches would be needed to further determine how interactions between the apicomplexans and their host influence infection.

Composition and organization of kinetochores show plasticity in apicomplexan chromosome segregation

Kinetochores are multiprotein assemblies directing mitotic spindle attachment and chromosome segregation. In apicomplexan parasites, most known kinetochore components and associated regulators are apparently missing, suggesting a minimal structure with limited control over chromosome segregation. In this study, we use interactomics combined with deep homology searches to identify 13 previously unknown components of kinetochores in Apicomplexa. Apicomplexan kinetochores are highly divergent in sequence and composition from animal and fungal models. The nanoscale organization includes at least four discrete compartments, each displaying different biochemical interactions, subkinetochore localizations and evolutionary rates across the phylum. We reveal alignment of kinetochores at the metaphase plate in both Plasmodium berghei and Toxoplasma gondii, suggestive of a conserved "hold signal" that prevents precocious entry into anaphase. Finally, we show unexpected plasticity in kinetochore composition and segregation between apicomplexan lifecycle stages, suggestive of diverse requirements to maintain fidelity of chromosome segregation across parasite modes of division.

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.

Differential Regulation of an OsIspH1, the Functional 4-Hydroxy-3-Methylbut-2-Enyl Diphosphate Reductase, for Photosynthetic Pigment Biosynthesis in Rice Leaves and Seeds

The methylerythritol 4-phosphate (MEP) pathway is responsible for providing common precursors for the biosynthesis of diverse plastidial terpenoids, including chlorophylls, carotenoids, and phytohormones, in plants. In rice (Oryza sativa), the last-step genes encoding 4-hydroxy-3-methylbut-2-enyl diphosphate reductase [HDR/isoprenoid synthesis H (IspH)] have been annotated in two genes (OsIspH1 and OsIspH2) in the rice genome. The spatial transcript levels indicated that OsIspH1 is highly expressed in all tissues at different developmental stages, whereas OsIspH2 is barely expressed due to an early stop in exon 1 caused by splicing error. OsIspH1 localized into plastids and osisph1, a T-DNA inserted knockout mutant, showed an albino phenotype, indicating that OsIspH1 is the only functional gene. To elucidate the role of OsIspH1 in the MEP pathway, we created two single (H145P and K407R) and double (H145P/K407R) mutations and performed complementation tests in two hdr mutants, including Escherichia coli DLYT1 strains and osisph1 rice plants. The results showed that every single mutation retained HDR function, but a double mutation lost it, proposing that the complementary relations of two residues might be important for enzyme activity but not each residue. When overexpressed in rice plants, the double-mutated gene, OsIspH1^MUT , reduced chlorophyll and carotenoid biosynthesis in the leaves and seeds. It confirmed the crucial role of OsIspH1 in plastidic terpenoid biosynthesis, revealing organ-specific differential regulation of OsIspH1 in rice plants.

Droplet digital PCR as a tool for investigating dynamics of cryptic symbionts

Interactions among symbiotic organisms and their hosts are major drivers of ecological and evolutionary processes. Monitoring the infection patterns among natural populations and identifying factors affecting these interactions are critical for understanding symbiont-host relationships. However, many of these interactions remain understudied since the knowledge about the symbiont species is lacking, which hinders the development of appropriate tools. In this study, we developed a digital droplet PCR (ddPCR) assay based on apicomplexan COX1 gene to detect an undescribed agamococcidian symbiont. We show that the method gives precise and reproducible results and enables detecting cryptic symbionts in low target concentration. We further exemplify the assay's use to survey seasonally sampled natural host (Pygospio elegans) populations for symbiont infection dynamics. We found that symbiont prevalence differs spatially but does not show seasonal changes. Infection load differed between populations and was low in spring and significantly increased towards fall in all populations. We also found that the symbiont prevalence is affected by host length and population density. Larger hosts were more likely to be infected, and high host densities were found to have a lower probability of infection. The observed variations could be due to characteristics of both symbiont and host biology, especially the seasonal variation in encounter rates. Our findings show that the developed ddPCR assay is a robust tool for detecting undescribed symbionts that are otherwise difficult to quantify, enabling further insight into the impact cryptic symbionts have on their hosts.

Aggregata polibraxiona n. sp. (Apicomplexa: Aggregatidae) from Octopus bimaculatus Verrill, 1883 (Mollusca: Cephalopoda) from the Gulf of California, Mexico

The Apicomplexa Aggregata spp. are intracellular parasites of cephalopods that infect the intestinal tract of commercially important species such as Octopus bimaculatus, which sustains the octopus fishery in Baja California (B.C.), Mexico. In this study, Aggregata polibraxiona n. sp. was described from the cecum of O. bimaculatus collected from Bahia de Los Angeles, B. C. Light and electron microscopy revealed that oocysts and sporocysts were spherical to ovoid in shape. Sporulated oocysts (293-835 × 177-688 μm) contained 135-674 sporocysts (12-24 × 11-22 μm). The sporocyst wall was covered by tubular projections (0.55-2.19 μm in length) bifurcated in the top, unevenly distributed, covered by a thin membrane. Each sporocyst contains 11-13 sporozoites (16-26 × 1.20-3 μm). Three partial sequences of the 18S rDNA gene were obtained, and two phylogenetic approaches were performed according to Bayesian inference and Maximum Likelihood. In both phylogenetic reconstructions, the sequences of A. polibraxiona n. sp. were recovered as a monophyletic group within the genus Aggregata and placed as a sister group to Aggregata octopiana Lineage II. Aggregata polibraxiona n. sp. is the first Apicomplexa described from a cephalopod host from Mexico and extends the geographical range of Apicomplexa infecting cephalopods.

Parallel functional reduction in the mitochondria of apicomplexan parasites

Gregarines are an early-diverging lineage of apicomplexan parasites that hold many clues into the origin and evolution of the group, a remarkable transition from free-living phototrophic algae into obligate parasites of animals.¹ Using single-cell transcriptomics targeting understudied lineages to complement available sequencing data, we characterized the mitochondrial metabolic repertoire across the tree of apicomplexans. In contrast to the large suite of proteins involved in aerobic respiration in well-studied parasites like Toxoplasma or Plasmodium,² we find that gregarine trophozoites have significantly reduced energy metabolism: most lack respiratory complexes III and IV, and some lack the electron transport chains (ETCs) and tricarboxylic acid (TCA) cycle entirely. Phylogenomic analyses show that these reductions took place several times in parallel, resulting in a functional range from fully aerobic organelles to extremely reduced "mitosomes" restricted to Fe-S cluster biosynthesis. The mitochondrial genome has also been lost repeatedly: in species with severe functional reduction simply by gene loss but in one species with a complete ETC by relocating cox1 to the nuclear genome. Severe functional reduction of mitochondria is generally associated with structural reduction, resulting in small, nondescript mitochondrial-related organelles (MROs).³ By contrast, gregarines retain distinctive mitochondria with tubular cristae, even the most functionally reduced cases that also lack genes associated with cristae formation. Overall, the parallel, severe reduction of gregarine mitochondria expands the diversity of organisms that contain MROs and further emphasizes the role of parallel transitions in apicomplexan evolution.

The apicoplast of Haemoproteus columbae: A comparative study of this organelle genome in Haemosporida

Apicomplexa is a phylum of parasitic protozoa; among them are the order Haemosporida, vector-borne parasites that include those that cause malaria (genus Plasmodium). Most Apicomplexa species have a non-photosynthetic plastid or apicoplast. Given its unique metabolic pathways, this organelle is considered a target for malaria therapeutics. Regardless of its importance, there is a paucity of complete apicoplast genome data hindering comparative studies. Here, the Haemoproteus (Haemoproteus) columbae apicoplast genome (lineage HAECOL1) was obtained using next-generation sequencing. This genome was included in a comparative analysis with other plastids. This 29.8 kb circular genome shares the same structure found in Plasmodium parasites. It is A + T rich (87.7%), comparable but at the higher end of A + T content observed in Plasmodium species (85.5-87.2%). As expected, considering its high A + T content, the synonymous codon usage (RSCU) and the effective number of codons (ENc) showed a moderate codon bias. Several apicoplast genes have a phylogenetic signal. However, unlike mitochondrial genes, single-gene phylogenies have low support in haemosporidian clades that diverged recently. The H. columbae apicoplast genome suggests that the apicoplast function may be conserved across Haemosporida. This parasite could be a model to study this organelle in a non-mammalian system.

Peptidylarginine Deiminase (PAD) and Post-Translational Protein Deimination-Novel Insights into Alveolata Metabolism, Epigenetic Regulation and Host-Pathogen Interactions

The alveolates (Superphylum Alveolata) comprise a group of primarily single-celled eukaryotes that have adopted extremely diverse modes of nutrition, such as predation, photoautotrophy and parasitism. The alveolates consists of several major phyla including the apicomplexans, a large group of unicellular, spore forming obligate intracellular parasites, and chromerids, which are believed to be the phototrophic ancestors of the parasitic apicomplexans. Molecular pathways involved in Alveolata host-pathogen interactions, epigenetic regulation and metabolism in parasite development remain to be fully understood. Peptidylarginine deiminases (PADs) are a phylogenetically conserved enzyme family which causes post-translational protein deimination, affecting protein function through the conversion of arginine to citrulline in a wide range of target proteins, contributing to protein moonlighting in physiological and pathological processes. The identification of deiminated protein targets in alveolate parasites may therefore provide novel insight into pathogen survival and host-pathogen interactions. The current study assessed PAD homologues and deiminated protein profiles of two alveolate parasites, Piridium sociabile (Chromerida) and Merocystis kathae (Apicomplexa). Histological analysis verified strong cytoplasmic PAD expression in both Alveolates, detected deiminated proteins in nuclear and cytoplasmic compartments of the alveolate parasites and verified the presence of citrullinated histone H3 in Alveolata nucleus, indicating roles in epigenetic regulation. Histone H3 citrullination was also found significantly elevated in the host tissue, indicative of neutrophil extracellular trap formation, a host-defence mechanism against a range of pathogens, particularly those that are too large for phagocytosis. Proteomic analysis of deiminated proteins from both Alveolata identified GO and KEGG pathways strongly relating to metabolic and genetic regulation, with some species-specific differences between the apicomplexan and the chromerid. Our findings provide novel insights into roles for the conserved PAD/ADI enzyme family in the regulation of metabolic and epigenetic pathways in alveolate parasites, possibly also relating to their life cycle and host-pathogen interactions.