Several hundred ciliate species live in animals\' guts as a part of their microbiome. Among them, Muniziella cunhai (Trichostomatia, Pycnotrichidae), the largest described ciliate, is found exclusively associated with Hydrochoerus hydrochaeris (capybara), the largest known rodent reaching up to 90 kg. Here, we present the sequence, structural and functional annotation of this giant microeukaryote macronuclear genome and discuss its phylogenetic placement. The 85 Mb genome is highly AT rich (GC content 25.71 %) and encodes a total of 11 397 protein-coding genes, of which 2793 could have their functions predicted with automated functional assignments. Functional annotation showed that M. cunhai can digest recalcitrant structural carbohydrates, non-structural carbohydrates, and microbial cell walls, suggesting a role in diet metabolization and in microbial population control in the capybara\'s intestine. Moreover, the phylogenetic placement of M. cunhai provides insights on the origins of gigantism in the subclass Trichostomatia.

The high-throughput proteomics data generated by increasingly more sensible mass spectrometers greatly contribute to our better understanding of molecular and cellular mechanisms operating in live beings. Nevertheless, proteomics analyses are based on accurate genomic and protein annotations, and some information may be lost if these resources are incomplete. Here, we show that most proteomics data may be recovered by interconnecting genomics and proteomics approaches (i.e., following a proteogenomic strategy), resulting, in turn, in an improvement of gene/protein models. In this study, we generated proteomics data from Leishmania donovani (HU3 strain) promastigotes that allowed us to detect 1908 proteins in this developmental stage on the basis of the currently annotated proteins available in public databases. However, when the proteomics data were searched against all possible open reading frames existing in the L. donovani genome, twenty new protein-coding genes could be annotated. Additionally, 43 previously annotated proteins were extended at their N-terminal ends to accommodate peptides detected in the proteomics data. Also, different post-translational modifications (phosphorylation, acetylation, methylation, among others) were found to occur in a large number of Leishmania proteins. Finally, a detailed comparative analysis of the L. donovani and Leishmania major experimental proteomes served to illustrate how inaccurate conclusions can be raised if proteomes are compared solely on the basis of the listed proteins identified in each proteome. Finally, we have created data entries (based on freely available repositories) to provide and maintain updated gene/protein models. Raw data are available via ProteomeXchange with the identifier PXD051920.

Dicytostelium firmibasis is a member of Dictyostelia, a group of social amoebae that upon starvation display aggregative multicellularity where the amoebae transition from uni- to multicellular life. The D. firmibasis genome assembly that is currently available is of limited use due to its low contiguity, large number of undetermined bases, and lack of annotations. Here we used Nanopore long read sequencing, complemented with Illumina sequencing, and developmental transcriptomics as well as small RNA-sequencing, to present a new, fully annotated, chromosome-level D. firmibasis genome assembly. The new assembly contains no undetermined bases, and consists mainly of six large contigs representing the chromosomes, as well as a complete mitochondrial genome. This new genome assembly will be a valuable tool, allowing comprehensive comparison to Dictyostelium discoideum, the dictyostelid genetically tractable model. Further, the new genome will be important for studies of evolutionary processes governing the transition from unicellular to multicellular organisms and will aid in the sequencing and annotation of other dictyostelids genomes, many of which are currently of poor quality.

Malaria-causing protozoa of the genus Plasmodium have exerted one of the strongest selective pressures on the human genome, and resistance alleles provide biomolecular footprints that outline the historical reach of these species1. Nevertheless, debate persists over when and how malaria parasites emerged as human pathogens and spread around the globe1,2. To address these questions, we generated high-coverage ancient mitochondrial and nuclear genome-wide data from P. falciparum, P. vivax and P. malariae from 16 countries spanning around 5,500 years of human history. We identified P. vivax and P. falciparum across geographically disparate regions of Eurasia from as early as the fourth and first millennia BCE, respectively; for P. vivax, this evidence pre-dates textual references by several millennia3. Genomic analysis supports distinct disease histories for P. falciparum and P. vivax in the Americas: similarities between now-eliminated European and peri-contact South American strains indicate that European colonizers were the source of American P. vivax, whereas the trans-Atlantic slave trade probably introduced P. falciparum into the Americas. Our data underscore the role of cross-cultural contacts in the dissemination of malaria, laying the biomolecular foundation for future palaeo-epidemiological research into the impact of Plasmodium parasites on human history. Finally, our unexpected discovery of P. falciparum in the high-altitude Himalayas provides a rare case study in which individual mobility can be inferred from infection status, adding to our knowledge of cross-cultural connectivity in the region nearly three millennia ago.

Giardia duodenalis, a major cause of waterborne infection, infects a wide range of mammalian hosts and is subdivided into eight genetically well-defined assemblages named A through H. However, fragmented genomes and a lack of comparative analysis within and between the assemblages render unclear the molecular mechanisms controlling host specificity and differential disease outcomes. To address this, we generated a near-complete de novo genome of AI assemblage using the Oxford Nanopore platform by sequencing the Be-2 genome. We generated 148,144 long-reads with quality scores of > 7. The final genome assembly consists of only nine contigs with an N50 of 3,045,186 bp. This assembly agrees closely with the assembly of another strain in the AI assemblage (WB-C6). However, a critical difference is that a region previously placed in the five-prime region of Chr5 belongs to Chr4 of Be-2. We find a high degree of conservation in the ploidy, homozygosity, and the presence of cysteine-rich variant-specific surface proteins (VSPs) within the AI assemblage. Our assembly provides a nearly complete genome of a member of the AI assemblage of G. duodenalis, aiding population genomic studies capable of elucidating Giardia transmission, host range, and pathogenicity.

Plasmodiophora brassicae (Woronin, 1877), a biotrophic, obligate parasite, is the causal agent of clubroot disease in brassicas. The clubroot pathogen has been reported in more than 80 countries worldwide, causing economic losses of hundreds of millions every year. Despite its widespread impact, very little is known about the molecular strategies it employs to induce the characteristic clubs in the roots of susceptible hosts during infection, nor about the mechanisms it uses to overcome genetic resistance. Here, we provide the first telomere-to-telomere complete genome of P. brassicae. We generated ∼27 Gb of Illumina, Oxford Nanopore, and PacBio HiFi data from resting spores of strain Pb3A and produced a 25.3 Mb assembly comprising 20 chromosomes, with an N50 of 1.37 Mb. The BUSCO score, the highest reported for any member of the group Rhizaria (Eukaryota: 88.2%), highlights the limitations within the Eukaryota database for members of this lineage. Using available transcriptomic data and protein evidence, we annotated the Pb3A genome, identifying 10,521 protein-coding gene models. This high-quality, complete genome of P. brassicae will serve as a crucial resource for the plant pathology community to advance the much-needed understanding of the evolution of the clubroot pathogen.

Microorganisms play a central role in sustaining soil ecosystems and agriculture, and these functions are usually associated with their complex life history. Yet, the regulation and evolution of life history have remained enigmatic and poorly understood, especially in protozoa, the third most abundant group of organisms in the soil. Here, we explore the life history of a cosmopolitan species-Colpoda steinii. Our analysis has yielded a high-quality macronuclear genome for C. steinii, with size of 155 Mbp and 37,123 protein-coding genes, as well as mean intron length of ~93 bp, longer than most other studied ciliates. Notably, we identify two possible whole-genome duplication events in C. steinii, which may account for its genome being about twice the size of C. inflata\'s, another co-existing species. We further resolve the gene expression profiles in diverse life stages of C. steinii, which are also corroborated in C. inflata. During the resting cyst stage, genes associated with cell death and vacuole formation are upregulated, and translation-related genes are downregulated. While the translation-related genes are upregulated during the excystment of resting cysts. Reproductive cysts exhibit a significant reduction in cell adhesion. We also demonstrate that most genes expressed in specific life stages are under strong purifying selection. This study offers a deeper understanding of the life history evolution that underpins the extraordinary success and ecological functions of microorganisms in soil ecosystems.IMPORTANCEColpoda species, as a prominent group among the most widely distributed and abundant soil microorganisms, play a crucial role in sustaining soil ecosystems and promoting plant growth. This investigation reveals their exceptional macronuclear genomic features, including significantly large genome size, long introns, and numerous gene duplications. The gene expression profiles and the specific biological functions associated with the transitions between various life stages are also elucidated. The vast majority of genes linked to life stage transitions are subject to strong purifying selection, as inferred from multiple natural strains newly isolated and deeply sequenced. This substantiates the enduring and conservative nature of Colpoda\'s life history, which has persisted throughout the extensive evolutionary history of these highly successful protozoa in soil. These findings shed light on the evolutionary dynamics of microbial eukaryotes in the ever-fluctuating soil environments. This integrative research represents a significant advancement in understanding the life histories of these understudied single-celled eukaryotes.

Toxoplasma gondii is a global protozoan pathogen. Clonal lineages predominate in Europe, North America, Africa, and China, whereas highly recombinant parasites are endemic in South/Central America. Far East Asian T. gondii isolates are not included in current global population genetic structure analyses at WGS resolution. Here we report a genome-wide population study that compared eight Japanese and two Chinese isolates against representative worldwide T. gondii genomes using POPSICLE, a novel population structure analyzing software. Also included were 7 genomes resurrected from non-viable isolates by target enrichment sequencing. Visualization of the genome structure by POPSICLE shows a mixture of Chinese haplogroup (HG) 13 haploblocks introgressed within the genomes of Japanese HG2 and North American HG12. Furthermore, two ancestral lineages were identified in the Japanese strains; one lineage shares a common ancestor with HG11 found in both Japanese strains and North American HG12. The other ancestral lineage, found in T. gondii isolates from a small island in Japan, is admixed with genetically diversified South/Central American strains. Taken together, this study suggests multiple ancestral links between Far East Asian and American T. gondii strains and provides insight into the transmission history of this cosmopolitan organism.

BACKGROUND: Predation is a fundamental mechanism for organisms to acquire energy, and various species have evolved diverse tools to enhance their hunting abilities. Among protozoan predators, raptorial Haptorian ciliates are particularly fascinating as they possess offensive extrusomes known as toxicysts, which are rapidly discharged upon prey contact. However, our understanding of the genetic processes and specific toxins involved in toxicyst formation and discharge is still limited.
RESULTS: In this study, we investigated the predation strategies and subcellular structures of seven Haptoria ciliate species and obtained their genome sequences using single-cell sequencing technology. Comparative genomic analysis revealed distinct gene duplications related to membrane transport proteins and hydrolytic enzymes in Haptoria, which play a crucial role in the production and discharge of toxicysts. Transcriptomic analysis further confirmed the abundant expression of genes related to membrane transporters and cellular toxins in Haptoria compared to Trichostomatia. Notably, polyketide synthases (PKS) and L-amino acid oxidases (LAAO) were identified as potentially toxin genes that underwent extensive duplication events in Haptoria.
CONCLUSIONS: Our results shed light on the evolutionary and genomic adaptations of Haptorian ciliates for their predation strategies in evolution and provide insights into their toxic mechanisms.

In dinoflagellates, a unique and extremely divergent genomic and nuclear organization has evolved. The highly unusual features of dinoflagellate nuclei and genomes include permanently condensed liquid crystalline chromosomes, primarily packaged by proteins other than histones, genes organized in very long unidirectional gene arrays, a general absence of transcriptional regulation, high abundance of the otherwise very rare DNA modification 5-hydroxymethyluracil (5-hmU), and many others. While most of these fascinating properties are originally identified in the 1970s and 1980s, they have not yet been investigated using modern genomic tools.
In this work, we address some of the outstanding questions regarding dinoflagellate genome organization by mapping the genome-wide distribution of 5-hmU (using both immunoprecipitation-based and basepair-resolution chemical mapping approaches) and of chromatin accessibility in the genome of the Symbiodiniaceae dinoflagellate Breviolum minutum. We find that the 5-hmU modification is preferentially enriched over certain classes of repetitive elements, often coincides with the boundaries between gene arrays, and is generally correlated with decreased chromatin accessibility, the latter otherwise being largely uniform along the genome. We discuss the potential roles of 5-hmU in the functional organization of dinoflagellate genomes and its relationship to the transcriptional landscape of gene arrays.
Our results provide the first window into the 5-hmU and chromatin accessibility landscapes in dinoflagellates.