Sexual recombination is a hallmark of eukaryotic evolution. Without recombination, asexual eukaryotes should succumb to deleterious mutations and more rapidly evolving pathogens. Giardia duodenalis, a parasitic protist, sits within one of the earliest-branching eukaryotic lineages and has no known sexual stage. Whether Giardia are 'ancient asexuals' has been long explored but is unresolved. Here, we find clear evidence of sex in Giardia and also discover an asexual sublineage that has a broader host range than its sexual ancestor. This asexual lineage is not ancient, and is accumulating deleterious mutations. Unlike its sexual counterparts, its genetic variation lacks the signatures of selection and Red Queen coevolution. We propose a new hypothesis that explains how a mutational meltdown during Muller's Ratchet might enable asexual pathogens to expand their host ranges transiently. Fittingly, our results suggest that Giardia is not the last exception to, but rather further evidence of, the essentiality of eukaryotic sex.

0.0AbstractCryptosporidium, a protozoan parasite of significant public health concern, is responsible for severe diarrheal diseases, particularly in immunocompromised individuals and young children in resource-limited settings. Analysis of whole genome next generation sequencing (NGS) data is a critical next step in improving our understanding ofCryptosporidiumepidemiology, transmission dynamics, and genetic diversity. However, effective analysis of NGS data in a public health context necessitates the development of robust, validated bioinformatics tools. Here, we present Parapipe, a modular ISO accreditable bioinformatics pipeline designed for high-throughput processing and analysis ofCryptosporidiumNGS datasets. Built using Nextflow DSL2 and containerized with Singularity, Parapipe is portable, scalable, and capable of end-to-end analyses, including quality control, variant calling, multiplicity of infection (MOI) investigations and phylogenomic clustering analysis.Using both simulated and real-world datasets, we demonstrate Parapipe’s ability to resolve genetic heterogeneity, identify mixed infections, and generate high-resolution phylogenomic insights. Here, we use it to carry out a comparison between whole genome single nucleotide polymorphism (wgSNP) typing and the conventionally used gp60 molecular typing scheme. Compared to existing pipelines, Parapipe uniquely integrates MOI analysis, enabling the differentiation of mixed infections and supporting epidemiological investigations. Parapipe’s design facilitates integration with geographic, demographic, epidemiological and environmental data, enhancing its utility for tracking transmission pathways and outbreak sources.Parapipe represents a significant advance in utilising genomics for public health surveillance ofCryptosporidium, offering a streamlined and reproducible framework for analysis with potential application to other pathogenic protozoa. By automating complex workflows and enabling detailed genomic characterization, Parapipe provides a valuable tool for public health agencies and researchers, supporting efforts to mitigate the global burden of cryptosporidiosis.

Cryptosporidium, a protozoan parasite of significant public health concern, is responsible for severe diarrhoeal disease, particularly in immunocompromised individuals and young children in resource-limited settings. Analysis of whole-genome next-generation sequencing (NGS) data is critical in improving our understanding of Cryptosporidium epidemiology, transmission and diversity. However, effective analysis of NGS data in a public health context necessitates the development of robust, validated computational tools. We present Parapipe, an ISO-accreditable bioinformatic pipeline for high-throughput analysis of NGS data from Cryptosporidium and related taxa. Built using Nextflow DSL2 and containerized with Singularity, Parapipe is modular, portable, scalable and designed for use by public health laboratories. Using both simulated and real Cryptosporidium datasets, we demonstrate the power of Parapipe's genomic analysis for generating epidemiological insights. We highlight how whole-genome analysis yields substantially greater phylogenetic resolution than conventional gp60 molecular typing in Cryptosporidium parvum. Uniquely, Parapipe facilitates the integration of mixed infection analysis and phylogenomic clustering with epidemiological metadata, representing a powerful tool in the investigation of complex transmission pathways and identification of outbreak sources. Parapipe significantly advances genomic surveillance of Cryptosporidium, offering a streamlined, reproducible analytical framework. By automating a complex workflow and delivering detailed genomic characterization, Parapipe provides a valuable tool for public health agencies and researchers, supporting efforts to mitigate the global burden of cryptosporidiosis.

Many helminths and protozoa are transmitted to humans through the consumption of contaminated food or water, and this underlines the importance of methods for their detection in these matrices. Due to the difficulties in isolating parasites prior to their identification, indirect detection methods are used, mostly relying upon targeted amplification of nucleic acids via PCR and/or qPCR. With the development of high throughput sequencing technologies, an untargeted detection method, shotgun metagenomics, became available. By sequencing the total DNA extracted from a given source, and through bioinformatics analyses of the sequencing reads, shotgun metagenomics allows profiling the entire microbial community therein present, including eukaryotes and, therefore, parasites. In this article, we reviewed the studies that specifically addressed the detection of parasites in food (n = 2) and water matrices (n = 10) by shotgun metagenomics. Most studies focused on wastewater samples and reported the detection of many parasites of human and veterinary importance from various areas of the world, highlighting the potential of shotgun metagenomics to provide important data for parasitic pathogens surveillance. After examining the different analytical workflows employed in these studies, which were not developed for detection of eukaryotes (or parasites), we identified two aspects deserving attention. First, that assignment based on short reads matching ribosomal sequences may generate false positives due to high sequence conservation among eukaryotic organisms. Second, that reassessing the relatively small number of reads of eukaryotic origin by a BLAST search can confirm, or deny, identification of parasitic pathogens.

Limited molecular data exist on zoonotic parasites Cryptosporidium spp., Giardia duodenalis, and Enterocytozoon bieneusi in Angus calves from Guizhou, China. This study constitutes the first molecular epidemiological survey of these pathogens in this region. 817 fecal samples from Angus calves across 7 intensive beef farms (Bijie City). Nested PCR methods targeting SSU rRNA (Cryptosporidium spp.), gp60 (Cryptosporidium bovis subtyping), bg/gdh/tpi (G. duodenalis), and ITS (E. bieneusi) coupled with DNA sequencing were employed. DNA sequences were analyzed against the NCBI. database. Statistical differences were assessed via a generalized linear mixed-effects model. Cryptosporidium spp. prevalence 23.5% (192/817; 95% CI 28.1–34.6%), with C. bovis predominating 89.6% (172/192; 95% CI 84.4–93.5%) and six subtypes (XXVIa-XXVIf). Highest infection in 4–8-week-olds 29.9% (143/479; 95% CI 25.8–34.1%) (p < 0.01). G. duodenalis: 31.3% (256/817; 95% CI 28.1–34.6%) positive, overwhelmingly assemblage E 97.6% (6/256; 95% CI 0.9–5.0%), zoonotic assemblage A was marginal 0.7% (6/817; 95% CI 0.3–1.6%). Farm-level variation exceeded 10-fold (e.g., Gantang: 55.0% (55/100; 95% CI 44.7–65.0%) vs. Tieshi: 4.9% (5/102; 95% CI 1.6–11.1%). E. bieneusi: prevalence 19.7% (161/817; 95% CI 17.0–22.6%), exclusively zoonotic genotypes BEB4: 49.7% (80/161; 95% CI 41.7–57.7%); I: 40.4% (65/161; 95% CI 32.7–48.4%). Strong diarrhea association (p < 0.01) and site-specific patterns (e.g., Guanyindong: 39.2%). While Giardia exhibited the highest prevalence (31.3%) with minimal zoonotic risk, Enterocytozoon—despite lower prevalence (19.7%)—posed the greatest public health threat due to exclusive circulation of human-pathogenic genotypes (BEB4/I) and significant diarrhea association, highlighting divergent control priorities for these enteric parasites in Guizhou calves. Management/Public health impact: Dominant zoonotic E. bieneusi genotypes (BEB4/I) necessitate: 1. Targeted treatment of 4–8-week-old Angus calves. 2. Manure biofermentation (≥55 °C, 3 days), and 3. UV-disinfection (≥1 mJ/cm2) for karst water to disrupt transmission in this high-humidity region.

The zoonotic parasite Cryptosporidium parvum is an important global cause of diarrheal disease in humans and young ruminants. Molecular typing is essential to track transmission routes and identify clusters of cases. Here, we developed a novel multi-locus sequence typing (MLST) scheme based on single nucleotide polymorphisms (SNPs) in unlinked markers. Coding regions with high variability were identified by comparing whole genome sequences (WGS) from 43 human- and 92 ruminant-derived C. parvum samples collected across Europe. We first selected 18 markers and showed that they provide high discrimination among the samples with WGS data, with 88% of the MLSTs being singletons. Next, we defined a MLST based on eight genetically unlinked markers and generated sequence data from 305 C. parvum samples, collected from four different host species and 13 European countries. We consolidated a set of 365 fully genotyped samples, characterized by the presence of 154 different MLSTs, 105 of which were singletons. Network analyses showed no complete clustering of samples by host species or country of origin at the European scale. We further showed that samples with gp60 subtypes that are common in Europe are divided into many MLSTs by the new scheme, highlighting its increased discriminatory ability. However, the applicability of the scheme in public health settings is limited by its cost, turnaround time, and scalability. To achieve discrimination of C. parvum samples based on SNPs, a large number of loci needs to be analysed, and this is feasible using amplicon sequencing technologies.

Abstract Background Giardiasis, caused by the protozoan parasite Giardia intestinalis, often presents a treatment challenge, particularly in terms of resistance to metronidazole. Despite extensive research, markers for metronidazole resistance have not yet been identified. Methods This study analysed 28 clinical samples of G. intestinalis from sub-assemblage AII, characterised by varying responses to metronidazole treatment. We focussed on copy number variation (CNV) of the multi-copy flavohemoprotein gene, analysed using digital polymerase chain reaction (dPCR) and next generation sequencing (NGS). Additionally, chromosomal ploidy was tested in 18 of these samples. Flavohemoprotein CNV was also assessed in 17 samples from other sub-assemblages. Results Analyses revealed variable CNVs of the flavohemoprotein gene among the isolates, with no correlation to clinical metronidazole resistance. Discrepancies in CNVs detected from NGS data were attributed to biases linked to the whole genome amplification. However, dPCR helped to clarify these discrepancies by providing more consistent CNV data. Significant differences in flavohemoprotein CNVs were observed across different G. intestinalis sub-assemblages. Notably, Giardia exhibits a propensity for aneuploidy, contributing to genomic variability within and between sub-assemblages. Conclusions The complexity of the clinical metronidazole resistance in Giardia is influenced by multiple genetic factors, including CNVs and aneuploidy. No significant differences in the CNV of the flavohemoprotein gene between isolates from metronidazole-resistant and metronidazole-sensitive cases of giardiasis were found, underscoring the need for further research to identify reliable genetic markers for resistance. We demonstrate that dPCR and NGS are robust methods for analysing CNVs and provide cross-validating results, highlighting their utility in the genetic analyses of this parasite. Graphical Abstract

The cestode Echinococcus multilocularis is the causative agent of alveolar echinococcosis, a fatal zoonotic parasitic disease of the northern hemisphere. Red foxes are the main reservoir hosts and, likely, the main drivers of the geographic spread of the disease in Europe. Knowledge of genetic relationships among E. multilocularis isolates at a European scale is key to understanding the dispersal characteristics of E. multilocularis. Hence, the present study aimed to describe the genetic diversity of E. multilocularis isolates obtained from different host species in 19 European countries. Based on the analysis of complete nucleotide sequences of the cob, atp6, nad2, nad1 and cox1 mitochondrial genes (4,968 bp), 43 haplotypes were inferred. Four haplotypes represented 62.56 % of the examined isolates (142/227), and one of these four haplotypes was found in each country investigated, except Svalbard, Norway. While the haplotypes from Svalbard were markedly different from all the others, mainland Europe appeared to be dominated by two main clusters, represented by most western, central and eastern European countries, and the Baltic countries and northeastern Poland, respectively. Moreover, one Asian-like haplotype was identified in Latvia and northeastern Poland. To better elucidate the presence of Asian genetic variants of E. multilocularis in Europe, and to obtain a more comprehensive Europe-wide coverage, further studies, including samples from endemic regions not investigated in the present study, especially some eastern European countries, are needed. Further, the present work proposes historical causes that may have contributed to shaping the current genetic variability of E. multilocularis in Europe.

Diet impacts human health, influencing body adiposity and the risk of developing cardiometabolic diseases. The gut microbiome is a key player in the diet-health axis, but while its bacterial fraction is widely studied, the role of micro-eukaryotes, including Blastocystis, is underexplored. We performed a global-scale analysis on 56,989 metagenomes and showed that human Blastocystis exhibits distinct prevalence patterns linked to geography, lifestyle, and dietary habits. Blastocystis presence defined a specific bacterial signature and was positively associated with more favorable cardiometabolic profiles and negatively with obesity (p < 1e-16) and disorders linked to altered gut ecology (p < 1e-8). In a diet intervention study involving 1,124 individuals, improvements in dietary quality were linked to weight loss and increases in Blastocystis prevalence (p = 0.003) and abundance (p < 1e-7). Our findings suggest a potentially beneficial role for Blastocystis, which may help explain personalized host responses to diet and downstream disease etiopathogenesis.

Purpose A new high-resolution next-generation sequencing (NGS)-based method was established to type closely related European type II Toxoplasma gondii strains. Methods T. gondii field isolates were collected from different parts of Europe and assessed by whole genome sequencing (WGS). In comparison to ME49 (a type II reference strain), highly polymorphic regions (HPRs) were identified, showing a considerable number of single nucleotide polymorphisms (SNPs). After confirmation by Sanger sequencing, 18 HPRs were used to design a primer panel for multiplex PCR to establish a multilocus Ion AmpliSeq typing method. Toxoplasma gondii isolates and T. gondii present in clinical samples were typed with the new method. The sensitivity of the method was tested with serially diluted reference DNA samples. Results Among type II specimens, the method could differentiate the same number of haplotypes as the reference standard, microsatellite (MS) typing. Passages of the same isolates and specimens originating from abortion outbreaks were identified as identical. In addition, seven different genotypes, two atypical and two recombinant specimens were clearly distinguished from each other by the method. Furthermore, almost all SNPs detected by the Ion AmpliSeq method corresponded to those expected based on WGS. By testing serially diluted DNA samples, the method exhibited a similar analytical sensitivity as MS typing. Conclusion The new method can distinguish different T. gondii genotypes and detect intra-genotype variability among European type II T. gondii strains. Furthermore, with WGS data additional target regions can be added to the method to potentially increase typing resolution.

While waters might be contaminated by oocysts from >40 Cryptosporidium species, only viable oocysts of C. parvum and C. hominis truly pose the main health risk to the immunocompetent population. Oocyst viability is also an important but often neglected risk factor in monitoring waterborne parasites. However, commonly used methods in water monitoring and surveys cannot distinguish species (microscopic observation) or oocyst viability (PCR), as dead oocysts in water could retain gross structure and DNA content for weeks to months. Here, we report new TaqMan qRT-PCR/qPCR assays for quantitative detection of viable C. parvum and C. hominis oocysts. By targeting a hypothetical protein-encoding gene cgd6_3920 that is highly expressed in oocysts and variable between species, the qRT-PCR/qPCR assays achieve excellent analytical specificity and sensitivity (limit of quantification [LOQ] = 0.25 and 1.0 oocyst/reaction). Using calibration curves, the number and ratio of viable oocysts in specimens could be calculated. Additionally, we also establish a TaqMan-18S qPCR for cost-effective screening of pan-Cryptosporidium-positive specimens (LOQ = 0.1 oocyst/reaction). The assay feasibility is validated using field water (N = 43) and soil (79) specimens from 17 locations in Changchun, China, which detects four Cryptosporidium species from seven locations, including three gp60-subtypes (i.e., IIdA19G1, IIdA17G1 and IIdA24G2) of C. parvum oocysts showing varied viability ratios. These new TaqMan q(RT)-PCR assays supplement current methods in the survey of waters and other samples (e.g., surfaces, foods and beverages), and are applicable to assessing the efficiency of oocyst deactivation protocols.

The zoonotic parasiteCryptosporidium parvumis a global cause of gastrointestinal disease in humans and ruminants. Sequence analysis of the highly polymorphicgp60gene enabled the classification ofC. parvumisolates into multiple groups (e.g., IIa, IIc, Id) and a large number of subtypes. In Europe, subtype IIaA15G2R1 is largely predominant and has been associated with many water- and food-borne outbreaks. In this study, we generated new whole-genome sequence (WGS) data from 123 human- and ruminant-derived isolates collected in 13 European countries and included other available WGS data from Europe, Egypt, China, and the United States (n = 72) in the largest comparative genomics study to date. We applied rigorous filters to exclude mixed infections and analyzed a data set from 141 isolates from the zoonotic groups IIa (n = 119) and IId (n = 22). Based on 28,047 high-quality, biallelic genomic SNPs, we identified three distinct and strongly supported populations: Isolates from China (IId) and Egypt (IIa and IId) formed population 1; a minority of European isolates (IIa and IId) formed population 2; and the majority of European (IIa, including all IIaA15G2R1 isolates) and all isolates from the United States (IIa) clustered in population 3. Based on analyses of the population structure, population genetics, and recombination, we show that population 3 has recently emerged and expanded throughout Europe to then, possibly from the United Kingdom, reach the United States, where it also expanded. The reason(s) for the successful spread of population 3 remain elusive, although genes under selective pressure uniquely in this population were identified.

The Inter European Union Reference Laboratories (EURLs) Working Group on Next Generation Sequencing (NGS) involves eight EURLs for microbiological food and feed hazards and has been working since 2017 to promote the adoption of NGS by the National Reference Laboratories (NRLs) in the European Union. This work illustrates the results of the first 5 years of activity. By working together, the EURLs involved have released guidance documents for assisting NRLs in all the steps of NGS, helping the transition from classical molecular methods towards whole genome sequencing while ensuring harmonization, with the final aim of improving preparedness in the use of NGS to characterize microbial hazards and trace the sources of infection.

Giardia duodenalis (syn. G. intestinalis, G. lamblia ) is a widespread gastrointestinal protozoan parasite with debated taxonomic status. Currently, eight distinct genetic sub-groups, termed assemblages A–H, are defined based on a few genetic markers. Assemblages A and B may represent distinct species and are both of human public health relevance. Genomic studies are scarce and the few reference genomes available, in particular for assemblage B, are insufficient for adequate comparative genomics. Here, by combining long- and short-read sequences generated by PacBio and Illumina sequencing technologies, we provide nine annotated genome sequences for reference from new clinical isolates (four assemblage A and five assemblage B parasite isolates). Isolates chosen represent the currently accepted classification of sub-assemblages AI, AII, BIII and BIV. Synteny over the whole genome was generally high, but we report chromosome-level translocations as a feature that distinguishes assemblage A from B parasites. Orthologue gene group analysis was used to define gene content differences between assemblage A and B and to contribute a gene-set-based operational definition of respective taxonomic units. Giardia is tetraploid, and high allelic sequence heterogeneity (ASH) for assemblage B vs. assemblage A has been observed so far. Noteworthy, here we report an extremely low ASH (0.002%) for one of the assemblage B isolates (a value even lower than the reference assemblage A isolate WB-C6). This challenges the view of low ASH being a notable feature that distinguishes assemblage A from B parasites, and low ASH allowed assembly of the most contiguous assemblage B genome currently available for reference. In conclusion, the description of nine highly contiguous genome assemblies of new isolates of G. duodenalis assemblage A and B adds to our understanding of the genomics and species population structure of this widespread zoonotic parasite.

AbstractCryptosporidiosis is a major global health problem and a primary cause of diarrhea, particularly in young children in low- and middle-income countries (LMICs). The zoonotic Cryptosporidium parvum and anthroponotic Cryptosporidium hominis cause most human infections. Here, we present a comprehensive whole-genome study of C. hominis, comprising 114 isolates from 16 countries within five continents. We detect two lineages with distinct biology and demography, which diverged circa 500 years ago. We consider these lineages two subspecies and propose the names C. hominis hominis and C. hominis aquapotentis (gp60 subtype IbA10G2). In our study, C. h. hominis is almost exclusively represented by isolates from LMICs in Africa and Asia and appears to have undergone recent population contraction. In contrast, C. h. aquapotentis was found in high-income countries, mainly in Europe, North America, and Oceania, and appears to be expanding. Notably, C. h. aquapotentis is associated with high rates of direct human-to-human transmission, which may explain its success in countries with well-developed environmental sanitation infrastructure. Intriguingly, we detected genomic regions of introgression following secondary contact between the subspecies. This resulted in high diversity and divergence in genomic islands of putative virulence genes, including muc5 (CHUDEA2_430) and a hypothetical protein (CHUDEA6_5270). This diversity is maintained by balancing selection, suggesting a co-evolutionary arms race with the host. Finally, we find that recent gene flow from C. h. aquapotentis to C. h. hominis, likely associated with increased human migration, maybe driving the evolution of more virulent C. hominis variants.

Cryptosporidium is a leading global cause of waterborne disease, with many reported outbreaks related to main water supplies. In August 2019, an outbreak of cryptosporidiosis involving 80 cases occurred among 114 vacationers in a small municipality located in the Tuscan-Emilian Apennines, north-eastern Italy. After excluding a potential food-borne outbreak, the epidemiological investigation focussed on the hypothesis of a waterborne outbreak. This was confirmed by the finding of Cryptosporidium oocysts in stools of the cases and in water samples from the municipal water network. Molecular characterisation revealed the zoonotic species Cryptosporidium parvum as the causative agent. A single subtype (IIdA25G1) was found among all cases, and in one of two positive water samples. The municipality's water supply used spring water that only received a disinfection treatment insufficient to inactivate the parasite. Possible entry means into the water mains were found through further environmental investigations. As these types of water supplies are particularly vulnerable to various environmental factors, a control system based on the risk assessment of each phase of the water supply chain is required to guarantee water safety. Effective methods for detection of protozoan pathogens, which are generally excluded from routine water supply analysis, should be applied.

The flagellated parasite Giardia duodenalis is a major and global cause of diarrhoeal disease. Eight genetically very distinct groups, known as assemblages A to H, have been recognized in the G. duodenalis species complex, two of which (assemblages A and B) infect humans and other mammalian hosts. Informative typing schemes are essential to understand transmission pathways, characterize outbreaks and trace zoonotic transmission. In this study, we evaluated a published multi-locus sequence typing (MLST) scheme for G. duodenalis assemblage A, which is based on six polymorphic markers. We genotyped 60 human-derived and 11 animal-derived G. duodenalis isolates collected in Europe and on other continents based on the published protocol. After retrieving previously published genotyping data and excluding isolates whose sequences showed allelic sequence heterozygosity, we analysed a dataset comprising 146 isolates. We identified novel variants at five of the six markers and identified 78 distinct MLST types in the overall dataset. Phylogenetic interpretation of typing data confirmed that sub-assemblage AII only comprises human-derived isolates, whereas sub-assemblage AI comprises all animal-derived isolates and a few human-derived isolates, suggesting limited zoonotic transmission. Within sub-assemblage AII, isolates from two outbreaks, which occurred in Sweden and Italy, respectively, had unique and distinct MLST types. Population genetic analysis showed a lack of clustering by geographical origin of the isolates. The MLST scheme evaluated provides sufficient discriminatory power for epidemiological studies of G. duodenalis assemblage A.