Long-read assembly and comparative evidence-based reanalysis of Cryptosporidium genome sequences reveal expanded transporter repertoire and duplication of entire chromosome ends including subtelomeric regions.

Cryptosporidiosis is a leading cause of waterborne diarrheal disease globally and an important contributor to mortality in infants and the immunosuppressed. Despite its importance, the Cryptosporidium community has only had access to a good, but incomplete, Cryptosporidium parvum IOWA reference genome sequence. Incomplete reference sequences hamper annotation, experimental design, and interpretation. We have generated a new C. parvum IOWA genome assembly supported by Pacific Biosciences (PacBio) and Oxford Nanopore long-read technologies and a new comparative and consistent genome annotation for three closely related species: C. parvum, Cryptosporidium hominis, and Cryptosporidium tyzzeri We made 1926 C. parvum annotation updates based on experimental evidence. They include new transporters, ncRNAs, introns, and altered gene structures. The new assembly and annotation revealed a complete Dnmt2 methylase ortholog. Comparative annotation between C. parvum, C. hominis, and C. tyzzeri revealed that most "missing" orthologs are found, suggesting that the biological differences between the species must result from gene copy number variation, differences in gene regulation, and single-nucleotide variants (SNVs). Using the new assembly and annotation as reference, 190 genes are identified as evolving under positive selection, including many not detected previously. The new C. parvum IOWA reference genome assembly is larger, gap free, and lacks ambiguous bases. This chromosomal assembly recovers all 16 chromosome ends, 13 of which are contiguously assembled. The three remaining chromosome ends are provisionally placed. These ends represent duplication of entire chromosome ends including subtelomeric regions revealing a new level of genome plasticity that will both inform and impact future research.

Hookworms dynamically respond to loss of Type 2 immune pressure.

The impact of the host immune environment on parasite transcription and fitness is currently unknown. It is widely held that hookworm infections have an immunomodulatory impact on the host, but whether the converse is true remains unclear. Immunity against adult-stage hookworms is largely mediated by Type 2 immune responses driven by the transcription factor Signal Transducer and Activator of Transcription 6 (STAT6). This study investigated whether serial passage of the rodent hookworm Nippostrongylus brasiliensis in STAT6-deficient mice (STAT6 KO) caused changes in parasites over time. After adaptation to STAT6 KO hosts, N. brasiliensis increased their reproductive output, feeding capacity, energy content, and body size. Using an improved N. brasiliensis genome, we found that these physiological changes corresponded with a dramatic shift in the transcriptional landscape, including increased expression of gene pathways associated with egg production, but a decrease in genes encoding neuropeptides, proteases, SCP/TAPS proteins, and transthyretin-like proteins; the latter three categories have been repeatedly observed in hookworm excreted/secreted proteins (ESPs) implicated in immunosuppression. Although transcriptional changes started to appear in the first generation of passage in STAT6 KO hosts for both immature and mature adult stages, downregulation of the genes putatively involved in immunosuppression was only observed after multiple generations in this immunodeficient environment. When STAT6 KO-adapted N. brasiliensis were reintroduced to a naive WT host after up to 26 generations, this progressive change in host-adaptation corresponded to increased production of inflammatory cytokines by the WT host. Surprisingly, however, this single exposure of STAT6 KO-adapted N. brasiliensis to WT hosts resulted in worms that were morphologically and transcriptionally indistinguishable from WT-adapted parasites. This work uncovers remarkable plasticity in the ability of hookworms to adapt to their hosts, which may present a general feature of parasitic nematodes.

3, 2, 1, go! Cryptosporidium counts down to sex.

Cryptosporidium is a leading cause of death from childhood diarrhea, but its biology is poorly understood. A recent study in PLOS Biology reveals hitherto unknown aspects of the parasite's life cycle that may lead to improvements in ex vivo culture.

Conservation, abundance, glycosylation profile, and localization of the TSP protein family in Cryptosporidium parvum.

Cryptosporidium parvum is a zoonotic apicomplexan parasite and a common cause of diarrheal disease worldwide. The development of vaccines to prevent or limit infection remains an important goal for tackling cryptosporidiosis. At present, the only approved vaccine against any apicomplexan parasite targets a conserved adhesin possessing a thrombospondin repeat domain. C. parvum possesses 12 orthologous thrombospondin repeat domain-containing proteins known as CpTSP1-12, though little is known about these potentially important antigens. Here, we explore the architecture and conservation of the CpTSP protein family, as well as their abundance at the protein level within the sporozoite stage of the life cycle. We examine the glycosylation states of these proteins using a combination of glycopeptide enrichment techniques to demonstrate that these proteins are modified with C-, O-, and N-linked glycans. Using expansion microscopy, and an antibody against the C-linked mannose that is unique to the CpTSP protein family within C. parvum, we show that these proteins are found both on the cell surface and in structures that resemble the secretory pathway of C. parvum sporozoites. Finally, we generated a polyclonal antibody against CpTSP1 to show that it is found at the cell surface and within micronemes, in a pattern reminiscent of other apicomplexan motility-associated adhesins, and is present both in sporozoites and meronts. This work sheds new light on an understudied family of C. parvum proteins that are likely to be important to both parasite biology and the development of vaccines against cryptosporidiosis.

Precise gene models using long-read sequencing reveal a unique poly(A) signal in Giardia lamblia.

During pre-mRNA processing, the poly(A) signal is recognized by a protein complex that ensures precise cleavage and polyadenylation of the nascent transcript. The location of this cleavage event establishes the length and sequence of the 3' UTR of an mRNA, thus determining much of its post-transcriptional fate. Using long-read sequencing, we characterize the polyadenylation signal and related sequences surrounding Giardia lamblia cleavage sites for over 2600 genes. We find that G. lamblia uses an AGURAA poly(A) signal, which differs from the mammalian AAUAAA. We also describe how G. lamblia lacks common auxiliary elements found in other eukaryotes, along with the proteins that recognize them. Further, we identify 133 genes with evidence of alternative polyadenylation. These results suggest that despite pared-down cleavage and polyadenylation machinery, 3' end formation still appears to be an important regulatory step for gene expression in G. lamblia.

Global Population Genomics of Two Subspecies of Cryptosporidium hominis during 500 Years of Evolution.

Cryptosporidiosis 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.

VIVID: A Web Application for Variant Interpretation and Visualization in Multi-dimensional Analyses.

Large-scale comparative genomics- and population genetic studies generate enormous amounts of polymorphism data in the form of DNA variants. Ultimately, the goal of many of these studies is to associate genetic variants to phenotypes or fitness. We introduce VIVID, an interactive, user-friendly web application that integrates a wide range of approaches for encoding genotypic to phenotypic information in any organism or disease, from an individual or population, in three-dimensional (3D) space. It allows mutation mapping and annotation, calculation of interactions and conservation scores, prediction of harmful effects, analysis of diversity and selection, and 3D visualization of genotypic information encoded in Variant Call Format on AlphaFold2 protein models. VIVID enables the rapid assessment of genes of interest in the study of adaptive evolution and the genetic load, and it helps prioritizing targets for experimental validation. We demonstrate the utility of VIVID by exploring the evolutionary genetics of the parasitic protist Plasmodium falciparum, revealing geographic variation in the signature of balancing selection in potential targets of functional antibodies.

Recent genetic exchanges and admixture shape the genome and population structure of the zoonotic pathogen Cryptosporidium parvum.

Cryptosporidium parvum is a globally distributed zoonotic pathogen and a major cause of diarrhoeal disease in humans and ruminants. The parasite's life cycle comprises an obligatory sexual phase, during which genetic exchanges can occur between previously isolated lineages. Here, we compare 32 whole genome sequences from human- and ruminant-derived parasite isolates collected across Europe, Egypt and China. We identify three strongly supported clusters that comprise a mix of isolates from different host species, geographic origins, and subtypes. We show that: (1) recombination occurs between ruminant isolates into human isolates; (2) these recombinant regions can be passed on to other human subtypes through gene flow and population admixture; (3) there have been multiple genetic exchanges, and most are probably recent; (4) putative virulence genes are significantly enriched within these genetic exchanges, and (5) this results in an increase in their nucleotide diversity. We carefully dissect the phylogenetic sequence of two genetic exchanges, illustrating the long-term evolutionary consequences of these events. Our results suggest that increased globalization and close human-animal contacts increase the opportunity for genetic exchanges between previously isolated parasite lineages, resulting in spillover and spillback events. We discuss how this can provide a novel substrate for natural selection at genes involved in host-parasite interactions, thereby potentially altering the dynamic coevolutionary equilibrium in the Red Queens arms race.

Eukaryote-Conserved Methylarginine Is Absent in Diplomonads and Functionally Compensated in Giardia.

Methylation is a common posttranslational modification of arginine and lysine in eukaryotic proteins. Methylproteomes are best characterized for higher eukaryotes, where they are functionally expanded and evolved complex regulation. However, this is not the case for protist species evolved from the earliest eukaryotic lineages. Here, we integrated bioinformatic, proteomic, and drug-screening data sets to comprehensively explore the methylproteome of Giardia duodenalis-a deeply branching parasitic protist. We demonstrate that Giardia and related diplomonads lack arginine-methyltransferases and have remodeled conserved RGG/RG motifs targeted by these enzymes. We also provide experimental evidence for methylarginine absence in proteomes of Giardia but readily detect methyllysine. We bioinformatically infer 11 lysine-methyltransferases in Giardia, including highly diverged Su(var)3-9, Enhancer-of-zeste and Trithorax proteins with reduced domain architectures, and novel annotations demonstrating conserved methyllysine regulation of eukaryotic elongation factor 1 alpha. Using mass spectrometry, we identify more than 200 methyllysine sites in Giardia, including in species-specific gene families involved in cytoskeletal regulation, enriched in coiled-coil features. Finally, we use known methylation inhibitors to show that methylation plays key roles in replication and cyst formation in this parasite. This study highlights reduced methylation enzymes, sites, and functions early in eukaryote evolution, including absent methylarginine networks in the Diplomonadida. These results challenge the view that arginine methylation is eukaryote conserved and demonstrate that functional compensation of methylarginine was possible preceding expansion and diversification of these key networks in higher eukaryotes.

Tunable Harmonic Flow Patterns in Microfluidic Systems through Simple Tube Oscillation.

Generation of tunable harmonic flows at low cost in microfluidic systems is a persistent and significant obstacle to this field, substantially limiting its potential to address major scientific questions and applications. This work introduces a simple and elegant way to overcome this obstacle. Harmonic flow patterns can be generated in microfluidic structures by simply oscillating the inlet tubes. Complex rib and vortex patterns can be dynamically modulated by changing the frequency and magnitude of tube oscillation and the viscosity of liquid. Highly complex rib patterns and synchronous vortices can be generated in serially connected microfluidic chambers. Similar dynamic patterns can be generated using whole or diluted blood samples without damaging the sample. This method offers unique opportunities for studying complex fluids and soft materials, chemical synthesis of various compounds, and mimicking harmonic flows in biological systems using compact, tunable, and low-cost devices.

Extracellular vesicles from early stage Plasmodium falciparum-infected red blood cells contain PfEMP1 and induce transcriptional changes in human monocytes.

Pathogens can release extracellular vesicles (EVs) for cell-cell communication and host modulation. EVs from Plasmodium falciparum, the deadliest malaria parasite species, can transfer drug resistance genes between parasites. EVs from late-stage parasite-infected RBC (iRBC-EVs) are immunostimulatory and affect endothelial cell permeability, but little is known about EVs from early stage iRBC. We detected the parasite virulence factor PfEMP1, which is responsible for iRBC adherence and a major contributor to disease severity, in EVs, only up to 12-hr post-RBC invasion. Furthermore, using PfEMP1 transport knockout parasites, we determined that EVs originated from inside the iRBC rather than the iRBC surface. Proteomic analysis detected 101 parasite and 178 human proteins in iRBC-EVs. Primary human monocytes stimulated with iRBC-EVs released low levels of inflammatory cytokines and showed transcriptomic changes. Stimulation with iRBC-EVs from PfEMP1 knockout parasites induced more gene expression changes and affected pathways involved in defence response, stress response, and response to cytokines, suggesting a novel function of PfEMP1 when present in EVs. We show for the first time the presence of PfEMP1 in early stage P. falciparum iRBC-EVs and the effects of these EVs on primary human monocytes, uncovering a new mechanism of potential parasite pathogenesis and host interaction.

Ascaris suum draft genome.

Parasitic diseases have a devastating, long-term impact on human health, welfare and food production worldwide. More than two billion people are infected with geohelminths, including the roundworms Ascaris (common roundworm), Necator and Ancylostoma (hookworms), and Trichuris (whipworm), mainly in developing or impoverished nations of Asia, Africa and Latin America. In humans, the diseases caused by these parasites result in about 135,000 deaths annually, with a global burden comparable with that of malaria or tuberculosis in disability-adjusted life years. Ascaris alone infects around 1.2 billion people and, in children, causes nutritional deficiency, impaired physical and cognitive development and, in severe cases, death. Ascaris also causes major production losses in pigs owing to reduced growth, failure to thrive and mortality. The Ascaris-swine model makes it possible to study the parasite, its relationship with the host, and ascariasis at the molecular level. To enable such molecular studies, we report the 273 megabase draft genome of Ascaris suum and compare it with other nematode genomes. This genome has low repeat content (4.4%) and encodes about 18,500 protein-coding genes. Notably, the A. suum secretome (about 750 molecules) is rich in peptidases linked to the penetration and degradation of host tissues, and an assemblage of molecules likely to modulate or evade host immune responses. This genome provides a comprehensive resource to the scientific community and underpins the development of new and urgently needed interventions (drugs, vaccines and diagnostic tests) against ascariasis and other nematodiases.

Divergent Transcriptional Responses to Physiological and Xenobiotic Stress in Giardia duodenalis.

Understanding how parasites respond to stress can help to identify essential biological processes. Giardia duodenalis is a parasitic protist that infects the human gastrointestinal tract and causes 200 to 300 million cases of diarrhea annually. Metronidazole, a major antigiardial drug, is thought to cause oxidative damage within the infective trophozoite form. However, treatment efficacy is suboptimal, due partly to metronidazole-resistant infections. To elucidate conserved and stress-specific responses, we calibrated sublethal metronidazole, hydrogen peroxide, and thermal stresses to exert approximately equal pressure on trophozoite growth and compared transcriptional responses after 24 h of exposure. We identified 252 genes that were differentially transcribed in response to all three stressors, including glycolytic and DNA repair enzymes, a mitogen-activated protein (MAP) kinase, high-cysteine membrane proteins, flavin adenine dinucleotide (FAD) synthetase, and histone modification enzymes. Transcriptional responses appeared to diverge according to physiological or xenobiotic stress. Downregulation of the antioxidant system and α-giardins was observed only under metronidazole-induced stress, whereas upregulation of GARP-like transcription factors and their subordinate genes was observed in response to hydrogen peroxide and thermal stressors. Limited evidence was found in support of stress-specific response elements upstream of differentially transcribed genes; however, antisense derepression and differential regulation of RNA interference machinery suggest multiple epigenetic mechanisms of transcriptional control.

First genetic analysis of Cryptosporidium from humans from Tasmania, and identification of a new genotype from a traveller to Bali.

Little is known about the molecular composition of Cryptosporidium species from humans living in the insular state of Tasmania, Australia. In the present study, we genetically characterized 82 samples of Cryptosporidium from humans following conventional coproscopic testing in a routine, diagnostic laboratory. Using a PCR-coupled single-strand conformation polymorphism (SSCP) technique, targeting portions of the small subunit rRNA (SSU), and 60 kDa glycoprotein (gp60) loci, we identified two species of Cryptosporidium, including C. hominis (subgenotypes IbA10G2, IdA16, IeA12G3T3, and IfA19G1) and C. parvum (IIaA16G1R1 and IIaA18G3), and a new operational taxonomic unit (OTU) that genetically closely resembled C. wrairi. This OTU was further characterized using markers in the actin, Cryptosporidium oocyst wall protein (COWP), and 70 kDa heat shock protein (hsp70) genes. This study provides the first characterization of species and genotypes of Cryptosporidium from Tasmania, and presents clear genetic evidence, using five independent genetic loci, for a new genotype or species of Cryptosporidium in a Tasmanian person with a recent history of travelling to Bali, Indonesia. It would be interesting to undertake detailed molecular-based studies of Cryptosporidium in Indonesia and neighbouring countries, in conjunction with morphological and experimental investigations of new genotypes.

Analyses of mitochondrial amino acid sequence datasets support the proposal that specimens of Hypodontus macropi from three species of macropodid hosts represent distinct species.

BACKGROUND: Hypodontus macropi is a common intestinal nematode of a range of kangaroos and wallabies (macropodid marsupials). Based on previous multilocus enzyme electrophoresis (MEE) and nuclear ribosomal DNA sequence data sets, H. macropi has been proposed to be complex of species. To test this proposal using independent molecular data, we sequenced the whole mitochondrial (mt) genomes of individuals of H. macropi from three different species of hosts (Macropus robustus robustus, Thylogale billardierii and Macropus [Wallabia] bicolor) as well as that of Macropicola ocydromi (a related nematode), and undertook a comparative analysis of the amino acid sequence datasets derived from these genomes. RESULTS: The mt genomes sequenced by next-generation (454) technology from H. macropi from the three host species varied from 13,634 bp to 13,699 bp in size. Pairwise comparisons of the amino acid sequences predicted from these three mt genomes revealed differences of 5.8% to 18%. Phylogenetic analysis of the amino acid sequence data sets using Bayesian Inference (BI) showed that H. macropi from the three different host species formed distinct, well-supported clades. In addition, sliding window analysis of the mt genomes defined variable regions for future population genetic studies of H. macropi in different macropodid hosts and geographical regions around Australia. CONCLUSIONS: The present analyses of inferred mt protein sequence datasets clearly supported the hypothesis that H. macropi from M. robustus robustus, M. bicolor and T. billardierii represent distinct species.

Genetic characterization of selected parasites from people with histories of gastrointestinal disorders using a mutation scanning-coupled approach.

A SSCP analysis and targeted sequencing approach was used for the genetic characterization of some major pathogens from a cohort of 227 people with histories of gastrointestinal disorders. Genomic DNAs from fecal samples were subjected to PCR-amplification of regions in the glycoprotein (gp60) or triose phosphate isomerase (tpi) gene, or the second internal transcribed spacer of nuclear ribosomal DNA (ITS-2). Cryptosporidium, Giardia, and strongylid nematodes were detected in 94, 132 and 12 samples. Cryptosporidium hominis subgenotypes IbA10G2, IdA15G1, IgA17, IgA18, and IfA13G1 were identified in 74.6, 16.9, 5.6, 1.4, and 1.4% of 71 samples, respectively. For Cryptosporidium parvum, subgenotypes IIaA17G2R1 (47.6%) and IIaA18G3R1 (23.8%) were identified in 23 samples. Giardia duodenalis assemblage B (78%) was more common than assemblage A (22%). In addition, DNA of the nematodes Ancylostoma ceylanicum (n = 2), Ancylostoma duodenale (4), Necator americanus (5), and Haemonchus contortus (1) was specifically detected. This is the first report of A. ceylanicum in two persons in Australia and, we provide molecular evidence of H. contortus in a child. This SSCP-based approach should provide a useful diagnostic and analytical tool for a wide range of pathogens.

Assessing calves as carriers of Cryptosporidium and Giardia with zoonotic potential on dairy and beef farms within a water catchment area by mutation scanning.

In the present study, we undertook a molecular epidemiological survey of Cryptosporidium and Giardia in calves on three dairy and two beef farms within an open drinking water catchment area (Melbourne, Australia). Faecal samples (n = 474) were collected from calves at two time points (5 months apart) and tested using a PCR-based mutation scanning-targeted sequencing phylogenetic approach, employing regions within the genes of small subunit (SSU) of ribosomal RNA (designated partial SSU), 60 kDa glycoprotein (pgp60) and triose phosphate isomerase (ptpi) as genetic markers. Using partial SSU, the C. bovis, C. parvum, C. ryanae and a new genotype of Cryptosporidium were characterised from totals of 74 (15.6%), 35 (7.3%), 37 (7.8%) and 9 (1.9%) samples, respectively. Using pgp60, C. parvum genotype IIa subgenotype A18G3R1 was detected in 29 samples. Using ptpi, G. duodenalis assemblages A and E were detected in totals of 10 (2.1%) and 130 (27.4%) samples, respectively. The present study showed that a considerable proportion of dairy and beef calves in this open water catchment region excreted Cryptosporidium (i.e. subgenotype IIaA18G3R1) and Giardia (e.g. assemblage A) that are consistent with those infecting humans, inferring that they are of zoonotic importance. Future work should focus on exploring, in a temporal and spatial way, whether these parasites occur in the environment and water of the catchment reservoir.

Comparative evaluation of two DNA isolation techniques for PCR-based diagnosis of gastrointestinal nematode infections in sheep.

The specific diagnosis of gastrointestinal parasite infections in livestock is central to their control. PCR assays have been developed for routine diagnosis and to overcome limitations of classical methods. Central to the performance of such assays is the effective isolation of the nucleic acids from samples and the elimination of components that are inhibitory to PCR. Here, we directly compared two techniques for the isolation of DNA from strongylid nematode eggs from faecal samples from sheep, and assessed their performance in relation to the sensitivity and specificity of PCR, time required for DNA isolation and ease of use. The results showed differences in the performance of the two isolation techniques, subsequently effecting the PCR results. The main differences related to the time required for DNA isolation, and the elimination of inhibitory substances from the DNA isolated by one technique but not the other.

Rapid, multiplex-tandem PCR assay for automated detection and differentiation of toxigenic cyanobacterial blooms.

Cyanobacterial blooms are a major water quality issue and potential public health risk in freshwater, marine and estuarine ecosystems globally, because of their potential to produce cyanotoxins. To date, a significant challenge in the effective management of cyanobacterial has been an inability of classical microscopy-based approaches to consistently and reliably detect and differentiate toxic from non-toxic blooms. The potential of cyanobacteria to produce toxins has been linked to the presence of specific biosynthetic gene clusters. Here, we describe the application of a robotic PCR-based assay for the semi-automated and simultaneous detection of toxin biosynthesis genes of each of the toxin classes characterized to date for cyanobacteria [i.e., microcystins (MCYs), nodularins (NODs), cylindrospermopsins (CYNs) and paralytic shellfish toxins (PSTs)/saxitoxins (SXTs)]. We demonstrated high sensitivity and specificity for each assay using well-characterized, cultured isolates, and establish its utility as a quantitative PCR using DNA, clone and cell-based dilution series. In addition, we used 206 field-collected samples and 100 known negative controls to compare the performance of each assay with conventional PCR and direct toxin detection. We report a diagnostic specificity of 100% and a sensitivity of ≥97.7% for each assay.

First molecular characterization of Giardia duodenalis from goats in Malaysia.

In the present study, 310 faecal samples from goats from eight different farms in Malaysia were tested for the presence of Giardia using a PCR-coupled approach. The nested PCR for SSU amplified products of the expected size (∼200 bp) from 21 of 310 (6.8%) samples. Sixteen of these 21 products could be sequenced successfully and represented six distinct sequence types. Phylogenetic analysis of the SSU sequence data using Bayesian Inference (BI) identified Giardia assemblages A, B and E. The identification of the 'zoonotic' assemblages A and B suggests that Giardia-infected goats represent a possible reservoir for human giardiasis in Malaysia.