Genome sequencing of Giardia lamblia genotypes A2 and B isolates (DH and GS) and comparative analysis with the genomes of genotypes A1 and E (WB and Pig)

Comparative genomics of Giardia duodenalis sub-assemblage AI beaver (Be-2) and human (WB-C6) strains show remarkable homozygosity, sequence similarity, and conservation of VSP genes

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.

Harnessing the power of new genetic tools to illuminate Giardia biology and pathogenesis

Giardia is a prevalent single-celled microaerophilic intestinal parasite causing diarrheal disease and significantly impacting global health. Double diploid (essentially tetraploid) Giardia trophozoites have presented a formidable challenge to the development of molecular genetic tools to interrogate gene function. High sequence divergence and the high percentage of hypothetical proteins lacking homology to proteins in other eukaryotes have limited our understanding of Giardia protein function, slowing drug target validation and development. For more than 25 years, Giardia A and B assemblages have been readily amenable to transfection with plasmids or linear DNA templates. Here, we highlight the utility and power of genetic approaches developed to assess protein function in Giardia, with particular emphasis on the more recent clustered regularly interspaced palindromic repeats/Cas9-based methods for knockdowns and knockouts. Robust and reliable molecular genetic approaches are fundamental toward the interrogation of Giardia protein function and evaluation of druggable targets. New genetic approaches tailored for the double diploid Giardia are imperative for understanding Giardia's unique biology and pathogenesis.

Genomic survey maps differences in the molecular complement of vesicle formation machinery between Giardia intestinalis assemblages

Giardia intestinalis is a globally important microbial pathogen with considerable public health, agricultural, and economic burden. Genome sequencing and comparative analyses have elucidated G. intestinalis to be a taxonomically diverse species consisting of at least eight different sub-types (assemblages A-H) that can infect a great variety of animal hosts, including humans. The best studied of these are assemblages A and B which have a broad host range and have zoonotic transmissibility towards humans where clinical Giardiasis can range from asymptomatic to diarrheal disease. Epidemiological surveys as well as previous molecular investigations have pointed towards critical genomic level differences within numerous molecular pathways and families of parasite virulence factors within assemblage A and B isolates. In this study, we explored the necessary machinery for the formation of vesicles and cargo transport in 89 Canadian isolates of assemblage A and B G. intestinalis. Considerable variability within the molecular complement of the endolysosomal ESCRT protein machinery, adaptor coat protein complexes, and ARF regulatory system have previously been reported. Here, we confirm inter-assemblage, but find no intra-assemblage variation within the trafficking systems examined. This variation includes losses of subunits belonging to the ESCRTIII as well as novel lineage specific duplications in components of the COPII machinery, ARF1, and ARFGEF families (BIG and CYTH). Since differences in disease manifestation between assemblages A and B have been controversially reported, our findings may well have clinical implications and even taxonomic, as the membrane trafficking system underpin parasite survival, pathogenesis, and propagation.

Spliceosomal introns in the diplomonad parasite Giardia duodenalis revisited

Complete reference genomes, including correct feature annotations, are a fundamental aspect of genomic biology. In the case of protozoan species such as Giardia duodenalis, a major human and animal parasite worldwide, accurate genome annotation can deepen our understanding of the evolution of parasitism and pathogenicity by identifying genes underlying key traits and clinically relevant cellular mechanisms, and by extension, the development of improved prevention strategies and treatments. This study used bioinformatics analyses of Giardia mRNA libraries to characterize known introns and identify new intron candidates, working towards completion of the G. duodenalis assemblage A strain 'WB' genome and further elucidating Giardia's gene expression. By using a set of experimentally validated positive control loci to calibrate our intron detection pipeline, we were able to detect evidence of previously missed candidate splice junctions directly from expressed transcript data. These intron candidates were further studied in silico using NMDS (non-metric multidimensional scaling) clustering to determine shared characteristics and their relative importance such as secondary structure, splicing efficiency and motif conservation, and thus to refine intron models. Results from this study identified 34 new intron candidates, with several potential introns showing evidence that secondary structure of the mRNA molecule might play a more significant role in splicing than previously reported eukaryotic splicing activity mediated by a reduced spliceosome present in G. duodenalis.

Genetic diversity and molecular diagnosis of Giardia

Giardia is a genus of flagellated protozoan parasites that infect the small intestine of humans and animals, causing the diarrheal illness known as giardiasis. Giardia exhibits significant genetic diversity among its isolates, which can have important implications for disease transmission and clinical presentation. This diversity is influenced by the coevolution of Giardia with its host, resulting in the development of unique genetic assemblages with distinct phenotypic characteristics. Although panmixia has not been observed, some assemblages appear to have a broader host range and exhibit higher transmission rates. Molecular diagnostic methods enable researchers to examine the genetic diversity of Giardia populations, enhancing our understanding of the genetic diversity, population structure, and transmission patterns of this pathogen and providing insights into clinical presentations of giardiasis.

Proposed taxonomic revision of Giardia duodenalis

Giardia duodenalis, Giardia enterica, Giardia intestinalis and Giardia lamblia are the synonyms for a species complex of 8-11 phylogenetically distinct species of Giardia infecting a broad range of animals including humans. Retrospective alignment of 8409 gene sequences from 3 loci confirmed host associations of Assemblages and sub-Assemblages within this species complex and molecular species delimitation testing confirmed that the Assemblages and sub-Assemblages AI and AII should be recognised as distinct species. It is recommended to synonymise the Assemblages with historic species descriptions based on host associations and consider descriptions for new species where no corresponding description exists. Synonyms, Giardia duodenalis, Giardia intestinalis and Giardia enterica, to be removed from synonymy: synonymise "Giardia duodenalis-Assemblage AI" syn. n. to Giardia duodenalis (Davaine, 1875), Kofoid and Christansen, 1915, synonymise "Giardia duodenalis-Assemblage AII" syn. n. to Giardia intestinalis (Lambl, 1859; Blanchard, 1885), Alexeieff, 1914 and synonymise "Giardia duodenalis-Assemblage B" syn. n. to Giardia enterica (Grassi, 1881a), Kofoid, 1920. Host specific Assemblages synonymised: synonymise canid-associated "Giardia duodenalis-Assemblage C" syn. n. to Giardia canisHegner, 1922; synonymise artiodactyl-associated "Giardia duodenalis-Assemblage E" syn. n. to Giardia bovisFantham, 1921; synonymise feline-associated "Giardia duodenalis-Assemblage F" syn. n. to Giardia catiDeschiens, 1925; and synonymise rodent-associated "Giardia duodenalis-Assemblage G" syn. n. to Giardia simoniLavier, 1924. New description for parasite type infecting specific host: canid-associated "Giardia duodenalis-Assemblage D" named Giardia lupus, sp. n. (LSID: urn:lsid:zoobank.org:act:1651A8CB-CBA8-40D9-AB59-D4AB11AC18A3). New proposed names and descriptions for consideration for parasite types infecting specific hosts: cervid-associated "Giardia duodenalis-sub-Assemblage AIII" for consideration "cervus" and Pinnipedia-associated "Giardia duodenalis-Assemblage H" for consideration "pinnipedis".

Genomic comparisons confirm Giardia duodenalis sub-assemblage AII as a unique species

Giardia duodenalis is a parasitic flagellated protozoan which infects a wide range of mammalian hosts, including humans, and is subdivided into at least eight genetic assemblages commonly thought to represent cryptic species. Molecular studies have shown that G. duodenalis assemblage A, which parasitizes humans and animals, contains several phylogenetically distinct groupings known as sub-assemblages. Molecular studies employing poor phylogenetic-resolution markers routinely recover these sub-assemblages, implying that they represent evolutionarily distinct clades and possibly cryptic species, a hypothesis which is supported by epidemiologic trends. Here, we further tested this hypothesis by using available data from 41 whole genomes to characterize sub-assemblages and coalescent techniques for statistical estimation of species boundaries coupled to functional gene content analysis, thereby assessing the stability and distinctiveness of clades. Our analysis revealed two new sub-assemblage clades as well as novel signatures of gene content geared toward differential host adaptation and population structuring via vertical inheritance rather than recombination or panmixia. We formally propose sub-assemblage AII as a new species, Giardia hominis, while preserving the name Giardia duodenalis for sub-assemblage AI. Additionally, our bioinformatic methods broadly address the challenges of identifying cryptic microbial species to advance our understanding of emerging disease epidemiology, which should be broadly applicable to other lower eukaryotic taxa of interest. Giardia hominis n. sp. Zoobank LSID: urn:lsid: zoobank.org:pub:4298F3E1-E3EF-4977-B9DD-5CC59378C80E.

Giardia and growth impairment in children in high-prevalence settings: consequence or co-incidence?

PURPOSE OF REVIEW

Giardia is a common intestinal parasite worldwide, and infection can be associated with clear, and sometimes persistent symptomatology. However, in children in high-prevalence settings, it is most often not associated with or is perhaps even protective against acute diarrhea. Nonetheless, recent longitudinal studies in high-prevalence settings increasingly identify an association with long-term outcomes that has been difficult to discern.

RECENT FINDINGS

Recent studies have made progress in disentangling this apparent paradox. First, prospective, well characterized cohort studies have repeatedly identified associations between Giardia infection, gut function, and child growth. Second, experimental animal and in-vitro models have further characterized the biological plausibility that Giardia could impair intestinal function and subsequently child development through different pathways, depending upon biological and environmental factors. Finally, new work has shed light on the potential for Giardia conspiring with specific other gut microbes, which may explain discrepant findings in the literature, help guide future higher resolution analyses of this pathogen, and inform new opportunities for intervention.

SUMMARY

Recent prospective studies have confirmed a high, if not universal, prevalence of persistent Giardia infections in low-and-middle income countries associated with child-growth shortfalls and altered gut permeability. However, the predominance of subclinical infections limits understanding of the true clinical impact of endemic pediatric giardiasis, and global disease burdens remain uncalculated. Integrating the role of Giardia in multipathogen enteropathies and how nutritional, microbial, metabolic, and pathogen-strain variables influence Giardia infection outcomes could sharpen delineations between pathogenic and potentially beneficial attributes of this enigmatic parasite.

Variation in haplotypes in single cysts of assemblages C and D, but not of assemblage E of Giardia duodenalis

Background

Giardia duodenalis, a single-celled intestinal parasite, is divided into eight assemblages (A-H), with differences in host specificity. Giardia duodenalis reproduces asexually and cycles between the binucleated trophozoite (4 N) and the infectious cyst with four nuclei (16 N). Interaction between the nuclei is limited. Therefore, genetic drift causes differences in genetic make-up between the non-daughter nuclei; the allelic sequence heterozygosity (ASH). The ASH is low (0.01%—0.0023%) for the related assemblages A and E, higher (0.43–0.53) for assemblage B and much higher (0.74% -0.89%) for the assemblage C and D at the root of the phylogenetic tree. The heterozygosity in assemblage F, in the same clade as assemblage A and E, was unknown. The heterozygosity in the sequences of the gdh and dis3 genes was used as proxy for the ASH and whole genome amplification of single cysts followed by cloning and Sanger sequencing of dis3 fragment could reveal the genetic variation within the cyst. The aim of the study was to determine the level of heterozygosity within pooled and single cysts of different assemblages.

Results

The heterozygosity in gdh and dis3 was determined in pooled cysts of the assemblages A to F. Heterozygosity in the isolates of the assemblages C (n = 2) and D (n = 1) ranged from 0.41% to 0.82% for gdh and dis3 and no heterozygosity was found in the isolates of the assemblages A (n = 4), E (n = 3) and F (n = 3). The heterozygosity in assemblage B (n = 7) was intermediate (0% to 0.62%). Next, the number of haplotypes of dis3 was determined for single cysts of assemblages C, D and E. In the assemblages C and D, two to four haplotypes were found per cyst, while in assemblage E only one haplotype was identified.

Conclusions

Having high heterozygosity is characteristic for the assemblages C and D, while having a low heterozygosity is characteristic for the clade with the assemblages A, E and F.

Presence of more than 1 haplotype per cyst in assemblage C and D suggests differences between the non-daughter nuclei, in contrast to the one haplotype in assemblage E.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12866-022-02581-3.

The controversies surrounding Giardia intestinalis assemblages A and B

Giardia intestinalis continues to be one of the most encountered parasitic diseases around the world. Although more frequently detected in developing countries, Giardia infections nonetheless pose significant public health problems in developed countries as well. Molecular characterisation of Giardia isolates from humans and animals reveals that there are two genetically different assemblages (known as assemblage A and B) that cause human infections. However, the current molecular assays used to genotype G. intestinalis isolates are quite controversial. This is in part due to a complex phenomenon where assemblages are incorrectly typed and underreported depending on which targeted locus is sequenced. In this review, we outline current knowledge based on molecular epidemiological studies and raise questions as to the reliability of current genotyping assays and a lack of a globally accepted method. Additionally, we discuss the clinical symptoms caused by G. intestinalis infection and how these symptoms vary depending on the assemblage infecting an individual. We also introduce the host-parasite factors that play a role in the subsequent clinical presentation of an infected person, and explore which assemblages are most seen globally.

A hybrid sequencing and assembly strategy for generating culture free Giardia genomes

Giardia duodenalis is a pathogenic intestinal protozoan parasite of humans and many other animals. Giardia duodenalis is found throughout the world, and infection is known to have adverse health consequences for human and other mammalian hosts. Yet, many aspects of the biology of this ubiquitous parasite remain unresolved. Whole genome sequencing and comparative genomics can provide insight into the biology of G. duodenalis by helping to reveal traits that are shared by all G. duodenalis assemblages or unique to an individual assemblage or strain. However, these types of analyses are currently hindered by the lack of available G. duodenalis genomes, due, in part, to the difficulty in obtaining the genetic material needed to perform whole genome sequencing. In this study, a novel approach using a multistep cleaning procedure coupled with a hybrid sequencing and assembly strategy was assessed for use in producing high quality G. duodenalis genomes directly from cysts obtained from feces of two naturally infected hosts, a cat and dog infected with assemblage A and D, respectively. Cysts were cleaned and concentrated using cesium chloride gradient centrifugation followed by immunomagnetic separation. Whole genome sequencing was performed using both Illumina MiSeq and Oxford Nanopore MinION platforms. A hybrid assembly strategy was found to produce higher quality genomes than assemblies from either platform alone. The hybrid G. duodenalis genomes obtained from fecal isolates (cysts) in this study compare favorably for quality and completeness against reference genomes of G. duodenalis from cultured isolates. The whole genome assembly for assemblage D is the most contiguous genome available for this assemblage and is an important reference genome for future comparative studies. The data presented here support a hybrid sequencing and assembly strategy as a suitable method to produce whole genome sequences from DNA obtained from G. duodenalis cysts which can be used to produce novel reference genomes necessary to perform comparative genomics studies of this parasite.

•

Assemblage A and D genomes were generated directly from cysts isolated from feces.

•

Genomes were sequenced using Illumina and Oxford Nanopore sequencing platforms.

•

A hybrid sequencing/assembly strategy was used to generate G. duodenalis genomes.

•

A hybrid strategy yields reference quality genomes from fecal isolates.

•

These methods have generated the most contiguous Assemblage D genome to date.

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

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

Genomic analysis finds no evidence of canonical eukaryotic DNA processing complexes in a free-living protist

Cells replicate and segregate their DNA with precision. Previous studies showed that these regulated cell-cycle processes were present in the last eukaryotic common ancestor and that their core molecular parts are conserved across eukaryotes. However, some metamonad parasites have secondarily lost components of the DNA processing and segregation apparatuses. To clarify the evolutionary history of these systems in these unusual eukaryotes, we generated a genome assembly for the free-living metamonad Carpediemonas membranifera and carried out a comparative genomics analysis. Here, we show that parasitic and free-living metamonads harbor an incomplete set of proteins for processing and segregating DNA. Unexpectedly, Carpediemonas species are further streamlined, lacking the origin recognition complex, Cdc6 and most structural kinetochore subunits. Carpediemonas species are thus the first known eukaryotes that appear to lack this suite of conserved complexes, suggesting that they likely rely on yet-to-be-discovered or alternative mechanisms to carry out these fundamental processes.

Unexpected organellar locations of ESCRT machinery in Giardia intestinalis and complex evolutionary dynamics spanning the transition to parasitism in the lineage Fornicata

BACKGROUND

Comparing a parasitic lineage to its free-living relatives is a powerful way to understand how that evolutionary transition to parasitism occurred. Giardia intestinalis (Fornicata) is a leading cause of gastrointestinal disease world-wide and is famous for its unusual complement of cellular compartments, such as having peripheral vacuoles instead of typical endosomal compartments. Endocytosis plays an important role in Giardia's pathogenesis. Endosomal sorting complexes required for transport (ESCRT) are membrane-deforming proteins associated with the late endosome/multivesicular body (MVB). MVBs are ill-defined in G. intestinalis, and roles for identified ESCRT-related proteins are not fully understood in the context of its unique endocytic system. Furthermore, components thought to be required for full ESCRT functionality have not yet been documented in this species.

RESULTS

We used genomic and transcriptomic data from several Fornicata species to clarify the evolutionary genome streamlining observed in Giardia, as well as to detect any divergent orthologs of the Fornicata ESCRT subunits. We observed differences in the ESCRT machinery complement between Giardia strains. Microscopy-based investigations of key components of ESCRT machinery such as GiVPS36 and GiVPS25 link them to peripheral vacuoles, highlighting these organelles as simplified MVB equivalents. Unexpectedly, we show ESCRT components associated with the endoplasmic reticulum and, for the first time, mitosomes. Finally, we identified the rare ESCRT component CHMP7 in several fornicate representatives, including Giardia and show that contrary to current understanding, CHMP7 evolved from a gene fusion of VPS25 and SNF7 domains, prior to the last eukaryotic common ancestor, over 1.5 billion years ago.

CONCLUSIONS

Our findings show that ESCRT machinery in G. intestinalis is far more varied and complete than previously thought, associates to multiple cellular locations, and presents changes in ESCRT complement which pre-date adoption of a parasitic lifestyle.

Giardia duodenalis: Biology and Pathogenesis

Giardia duodenalis captured the attention of Leeuwenhoek in 1681 while he was examining his own diarrheal stool, but, ironically, it did not really gain attention as a human pathogen until the 1960s, when outbreaks were reported. Key technological advances, including in vitro cultivation, genomic and proteomic databases, and advances in microscopic and molecular approaches, have led to an understanding that this is a eukaryotic organism with a reduced genome rather than a truly premitochondriate eukaryote. This has included the discovery of mitosomes (vestiges of mitochondria), a transport system with many of the features of the Golgi apparatus, and even evidence for a sexual or parasexual cycle. Cell biology approaches have led to a better understanding of how Giardia survives with two nuclei and how it goes through its life cycle as a noninvasive organism in the hostile environment of the lumen of the host intestine. Studies of its immunology and pathogenesis have moved past the general understanding of the importance of the antibody response in controlling infection to determining the key role of the Th17 response. This work has led to understanding of the requirement for a balanced host immune response that avoids the extremes of an excessive response with collateral damage or one that is unable to clear the organism. This understanding is especially important in view of the remarkable ranges of early manifestations, which range from asymptomatic to persistent diarrhea and weight loss, and longer-term sequelae that include growth stunting in children who had no obvious symptoms and a high frequency of postinfectious irritable bowel syndrome (IBS).

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.

Suitability of current typing procedures to identify epidemiologically linked human Giardia duodenalis isolates

Background

Giardia duodenalis is a leading cause of gastroenteritis worldwide. Humans are mainly infected by two different subtypes, i.e., assemblage A and B. Genotyping is hampered by allelic sequence heterozygosity (ASH) mainly in assemblage B, and by occurrence of mixed infections. Here we assessed the suitability of current genotyping protocols of G. duodenalis for epidemiological applications such as molecular tracing of transmission chains.

Methodology/Principal findings

Two G. duodenalis isolate collections, from an outpatient tropical medicine clinic and from several primary care laboratories, were characterized by assemblage-specific qPCR (TIF, CATH gene loci) and a common multi locus sequence typing (MLST; TPI, BG, GDH gene loci). Assemblage A isolates were further typed at additional loci (HCMP22547, CID1, RHP26, HCMP6372, DIS3, NEK15411).

Of 175/202 (86.6%) patients the G. duodenalis assemblage could be identified: Assemblages A 25/175 (14.3%), B 115/175 (65.7%) and A+B mixed 35/175 (20.0%). By incorporating allelic sequence heterozygosity in the analysis, the three marker MLST correctly identified 6/9 (66,7%) and 4/5 (80.0%) consecutive samples from chronic assemblage B infections in the two collections, respectively, and identified a cluster of five independent patients carrying assemblage B parasites of identical MLST type. Extended MLST for assemblage A altogether identified 5/6 (83,3%) consecutive samples from chronic assemblage A infections and 15 novel genotypes. Based on the observed A+B mixed infections it is estimated that only 75% and 50% of assemblage A or B only cases represent single strain infections, respectively. We demonstrate that typing results are consistent with this prediction.

Conclusions/Significance

Typing of assemblage A and B isolates with resolution for epidemiological applications is possible but requires separate genotyping protocols. The high frequency of multiple infections and their impact on typing results are findings with immediate consequences for result interpretation in this field.

Giardia duodenalis is a leading cause of gastroenteritis worldwide. Humans are mainly infected by the two different genetic subtypes, assemblage A and B. Molecular typing tools for epidemiological applications such as tracking transmission, attribution to a source and outbreak investigations have been developed and are highly desirable. However, to what degree the tetraploid genome with allelic sequence heterogeneity (ASH), and the frequent occurrence of mixed, assemblage A and B infections hamper performance is unclear. Here, we assessed the suitability of current genotyping protocols for deciphering the molecular epidemiology of G. duodenalis. Against a common reporting bias, we incorporated ASH in the analysis and we show that typing with resolution for epidemiological applications is possible for both, assemblage A and B isolates, but requires separate protocols. We also demonstrate how the high frequency of multiple infections overall impacts on typing results, which has immediate consequences for result interpretation in this field.

Molecular diversity of Giardia duodenalis in children under 5 years from the Manhiça district, Southern Mozambique enrolled in a matched case-control study on the aetiology of diarrhoea

Giardia duodenalis is an enteric parasite commonly detected in children. Exposure to this organism may lead to asymptomatic or symptomatic infection. Additionally, early-life infections by this protozoan have been associated with impaired growth and cognitive function in poor resource settings. The Global Enteric Multicenter Study (GEMS) in Mozambique demonstrated that G. duodenalis was more frequent among controls than in diarrhoeal cases (≥3 loosing stools in the previous 24 hours). However, no molecular investigation was conducted to ascertain the molecular variability of the parasite. Therefore, we describe here the frequency and genetic diversity of G. duodenalis infections in children younger than five years of age with and without diarrhoea from the Manhiça district in southern Mozambique enrolled in the context of GEMS. Genomic DNA from 757 G. duodenalis-positive stool samples by immunoassay collected between 2007-2012, were reanalysed by multiplex PCR targeting the E1-HP and C1-P21 genes for the differentiation of assemblages A and B. Overall, 47% (353) of the samples were successfully amplified in at least one locus. Assemblage B accounted for 90% (319/353) of all positives, followed by assemblage A (8%, 29/353) and mixed A+B infections (1%, 5/353). No association between the presence of a given assemblage and the occurrence of diarrhoea could be demonstrated. A total of 351 samples were further analysed by a multi-locus sequence genotyping (MLSG) approach at the glutamate dehydrogenase (gdh), ß-giardin (bg) and triose phosphate isomerase (tpi) genes. Overall, 63% (222/351) of samples were genotyped and/or sub-genotyped in at least one of the three markers. Sequence analysis revealed the presence of assemblages A (10%; 23/222) and B (90%; 199/222) with high molecular diversity at the nucleotide level within the latter; no mixed infections were identified under the MLSG scheme. Assemblage A sequences were assigned to sub-assemblages AI (0.5%, 1/222), AII (7%, 15/222) or ambiguous AII/AIII (3%, 7/222). Within assemblage B, sequences were assigned to sub-assemblages BIII (13%, 28/222), BIV (14%, 31/222) and ambiguous BIII/BIV (59%, 132/222). BIII/BIV sequences accumulated the majority of the single nucleotide polymorphisms detected, particularly in the form of double peaks at chromatogram inspection. This study demonstrated that the occurrence of gastrointestinal illness (diarrhoea) was not associated to a given genotype of G. duodenalis in Mozambican children younger than five years of age. The assemblage B of the parasite was responsible for nine out of ten infections detected in this paediatric population. The extremely high genetic diversity observed within assemblage B isolates was compatible with an hyperendemic epidemiological scenario where infections and reinfections were common. The obtained molecular data may be indicative of high coinfection rates by different G. duodenalis assemblages/sub-assemblages and/or genetic recombination events, although the exact contribution of both mechanisms to the genetic diversity of the parasite remains unknown.

Genetic Diversity of the Flavohemoprotein Gene of Giardia lamblia: Evidence for High Allelic Heterozygosity and Copy Number Variation

Purpose

The flavohemoprotein (gFlHb) in Giardia plays an important role in managing nitrosative and oxidative stress, and potentially also in virulence and nitroimidazole drug tolerance. The aim of this study was to analyze the genetic diversity of gFlHb in Giardia assemblages A and B clinical isolates.

Methods

gFlHb genes from 20 cultured clinical Giardia isolates were subjected to PCR amplification and cloning, followed by Sanger sequencing. Sequences of all cloned PCR fragments from each isolate were analyzed for single nucleotide variants (SNVs) and compared to genomic Illumina sequence data. Identical clone sequences were sorted into alleles, and diversity was further analyzed. The number of gFlHb gene copies was assessed by mining PacBio de novo assembled genomes in eight isolates. Homology models for assessment of SNV's potential impact on protein function were created using Phyre2.

Results

A variable copy number of the gFlHb gene, between two and six copies, depending on isolate, was found. A total of 37 distinct sequences, representing different alleles of the gFlHb gene, were identified in AII isolates, and 41 were identified in B isolates. In some isolates, up to 12 different alleles were found. The total allelic diversity was high for both assemblages (>0.9) and was coupled with a nucleotide diversity of <0.01. The genetic variation (SNVs per CDS length) was 4.8% in sub-assemblage AII and 5.4% in assemblage B. The number of non-synonymous (ns) SNVs was high in gFIHb of both assemblages, 1.6% in A and 3.0% in B, respectively. Some of the identified nsSNV are predicted to alter protein structure and possibly function.

Conclusion

In this study, we present evidence that gFlHb, a putative protective enzyme against oxidative and nitrosative stress in Giardia, is a variable copy number gene with high allelic diversity. The genetic variability of gFlHb may contribute metabolic adaptability against metronidazole toxicity.

Whole-genome sequencing of dog-specific assemblages C and D of Giardia duodenalis from single and pooled cysts indicates host-associated genes

Giardia duodenalis (syn. Giardia intestinalis or Giardia lamblia) infSAects over 280 million people each year and numerous animals. G. duodenalis can be subdivided into eight assemblages with different host specificity. Unculturable assemblages have so far resisted genome sequencing efforts. In this study, we isolated single and pooled cysts of assemblages C and D from dog faeces by FACS, and sequenced them using multiple displacement amplification and Illumina paired-end sequencing. The genomes of assemblages C and D were compared with genomes of assemblages A and B from humans and assemblage E from ruminants and pigs. The genomes obtained from the pooled cysts and from the single cysts were considered complete (>99 % marker genes observed) and the allelic sequence heterozygosity (ASH) values of assemblages C and D were 0.89 and 0.74 %, respectively. These ASH values were slightly higher than for assemblage B (>0.43 %) and much higher than for assemblages A and E, which ranged from 0.002 to 0.037 %. The flavohaemoglobin and 4Fe-4S binding domain family encoding genes involved in O₂ and NO detoxification were only present in assemblages A, B and E. Cathepsin B orthologs were found in all genomes. Six clades of cathepsin B orthologs contained one gene of each genome, while in three clades not all assemblages were represented. We conclude that whole-genome sequencing from a single Giardia cyst results in complete draft genomes, making the genomes of unculturable Giardia assemblages accessible. Observed differences between the genomes of assemblages C and D on one hand and the assemblages A, B and E on the other hand are possibly associated with host specificity.

Benchmarking hybrid assemblies of Giardia and prediction of widespread intra-isolate structural variation

Background

Currently available short read genome assemblies of the tetraploid protozoan parasite Giardia intestinalis are highly fragmented, highlighting the need for improved genome assemblies at a reasonable cost. Long nanopore reads are well suited to resolve repetitive genomic regions resulting in better quality assemblies of eukaryotic genomes. Subsequent addition of highly accurate short reads to long-read assemblies further improves assembly quality. Using this hybrid approach, we assembled genomes for three Giardia isolates, two with published assemblies and one novel, to evaluate the improvement in genome quality gained from long reads. We then used the long reads to predict structural variants to examine this previously unexplored source of genetic variation in Giardia.

Methods

With MinION reads for each isolate, we assembled genomes using several assemblers specializing in long reads. Assembly metrics, gene finding, and whole genome alignments to the reference genomes enabled direct comparison to evaluate the performance of the nanopore reads. Further improvements from adding Illumina reads to the long-read assemblies were evaluated using gene finding. Structural variants were predicted from alignments of the long reads to the best hybrid genome for each isolate and enrichment of key genes was analyzed using random genome sampling and calculation of percentiles to find thresholds of significance.

Results

Our hybrid assembly method generated reference quality genomes for each isolate. Consistent with previous findings based on SNPs, examination of heterozygosity using the structural variants found that Giardia BGS was considerably more heterozygous than the other isolates that are from Assemblage A. Further, each isolate was shown to contain structural variant regions enriched for variant-specific surface proteins, a key class of virulence factor in Giardia.

Conclusions

The ability to generate reference quality genomes from a single MinION run and a multiplexed MiSeq run enables future large-scale comparative genomic studies within the genus Giardia. Further, prediction of structural variants from long reads allows for more in-depth analyses of major sources of genetic variation within and between Giardia isolates that could have effects on both pathogenicity and host range.

Multilocus genotyping of Giardia duodenalis in Southwestern Iran. A community survey

Giardia duodenalis is one of the main enteric pathogens associated with diarrheal disease. In developing countries, giardiasis is a major public health concern, particularly in children under five years of age. This study aimed to evaluate the occurrence and genetic diversity of G. duodenalis causing human infections in Shushtar County, Southwestern Iran. Individual faecal specimens were collected from 1,163 individuals (male/female ratio: 0.9; age range 2-75 years) with (n = 258) and without (n = 905) gastrointestinal symptoms living in rural and urban settings during the period 2017-2018. Conventional (sucrose flotation and microscopy) methods were used for the initial detection of G. duodenalis cysts in faecal specimens. Microscopy-positive samples were confirmed by PCR amplification and sequencing of the small subunit rRNA (ssu rRNA) gene of the parasite. A multilocus genotyping (MLG) scheme targeting the triose phosphate isomerase (tpi), the glutamate dehydrogenase (gdh), and the beta-giardin (bg) genes was used for genotyping purposes. Giardia duodenalis cysts were detected in 7.7% (90/1,163) of samples by microscopy, of which 82 were confirmed by ssu-PCR. Successful amplification and sequencing results were obtained for 23.2% (19/82), 9.8% (8/82), and 8.5% (7/82) of the confirmed samples at the tpi, gdh, and bg loci, respectively. MLG data for the three loci were available for two samples only. Out of the 24 samples genotyped at any loci, 50% (12/24) were identified as assemblage A and the remaining half as assemblage B. Overall, AII was the most prevalent sub-assemblage detected (41.7%, 10/24), followed by BIII (25.0%, 6/24), discordant BIII/BIV (5/24) or AII/AIII (2/24) sequences, and BIV (1/24). No significant correlation was demonstrated between a given assemblage/sub-assemblage and the occurrence of clinical symptoms. No genotypes adapted to animal hosts other than humans (e.g. assemblages C-F) were found circulating in the investigated human population, suggesting that transmission of human giardiasis in this Iranian region is primarily of anthroponotic nature. Further molecular-based studies are needed to confirm and expand these results, and to ascertain the presence and public health relevance of the parasite in environmental (e.g. drinking water) samples.

Recent advances in functional research in Giardia intestinalis

This review considers current advances in tools to investigate the functional biology of Giardia, it's coding and non-coding genes, features and cellular and molecular biology. We consider major gaps in current knowledge of the parasite and discuss the present state-of-the-art in its in vivo and in vitro cultivation. Advances in in silico tools, including for the modelling non-coding RNAs and genomic elements, as well as detailed exploration of coding genes through inferred homology to model organisms, have provided significant, primary level insight. Improved methods to model the three-dimensional structure of proteins offer new insights into their function, and binding interactions with ligands, other proteins or precursor drugs, and offer substantial opportunities to prioritise proteins for further study and experimentation. These approaches can be supplemented by the growing and highly accessible arsenal of systems-based methods now being applied to Giardia, led by genomic, transcriptomic and proteomic methods, but rapidly incorporating advanced tools for detection of real-time transcription, evaluation of chromatin states and direct measurement of macromolecular complexes. Methods to directly interrogate and perturb gene function have made major leaps in recent years, with CRISPr-interference now available. These approaches, coupled with protein over-expression, fluorescent labelling and in vitro and in vivo imaging, are set to revolutionize the field and herald an exciting time during which the field may finally realise Giardia's long proposed potential as a model parasite and eukaryote.

A chromosome-scale reference genome for Giardia intestinalis WB

Giardia intestinalis is a protist causing diarrhea in humans. The first G. intestinalis genome, from the WB isolate, was published more than ten years ago, and has been widely used as the reference genome for Giardia research. However, the genome is fragmented, thus hindering research at the chromosomal level. We re-sequenced the Giardia genome with Pacbio long-read sequencing technology and obtained a new reference genome, which was assembled into near-complete chromosomes with only four internal gaps at long repeats. This new genome is not only more complete but also better annotated at both structural and functional levels, providing more details about gene families, gene organizations and chromosomal structure. This near-complete reference genome will be a valuable resource for the Giardia community and protist research. It also showcases how a fragmented genome can be improved with long-read sequencing technology completed with optical maps.

Measurement(s)	DNA • sequence_assembly • sequence feature annotation
Technology Type(s)	DNA sequencing • sequence assembly process • sequence annotation
Sample Characteristic - Organism	Giardia intestinalis

Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.11695659

Transcriptomic and proteomic analyses of Giardia intestinalis: Intestinal epithelial cell interactions

Giardia intestinalis is a unicellular protozoan parasite that infects the small intestines of humans and animals. Giardiasis, the disease caused by the parasite, occurs globally across socioeconomic boundaries but is mainly endemic in developing countries and particularly within young children, where pronounced effects manifests in a failure to thrive condition. The molecular pathogenesis of Giardia has been studied using in vitro models of human and rat intestinal epithelial cells (IECs) and parasites from the two major human genotypes or assemblages (A and B). High-quality, genome sequencing of representative isolates from assemblages A (WB) and B (GS) has enabled exploration of these host-parasite models using 'omics' technologies, allowing deep and quantitative analyses of global gene expression changes in IECs and parasites during their interactions, cross-talk and competition. These include a major up-regulation of immune-related genes in the IECs early after the start of interactions, as well as competition between host cells and parasites for nutrients like sugars, amino acids and lipids, which is also reflected in their secretome interactions. Unique parasite proteins dominate these interactions, with many major up-regulated genes being either hypothetical proteins or members of Giardia-specific gene families like the high-cysteine-rich membrane proteins (HCMPs), variable surface proteins (VSPs), alpha-giardins and cysteine proteases. Furthermore, these proteins also dominate in the secretomes, suggesting that they are important virulence factors in Giardia and crucial molecular effectors at the host-parasite interface.

Molecular epidemiology of giardiasis from a veterinary perspective

A total of eight Giardia species are accepted. These include: Giardia duodenalis (syn. Giardia intestinalis and Giardia lamblia), which infects humans and animals, Giardia agilis, Giardia ardeae, Giardia psittaci, Giardia muris, Giardia microti, Giardia peramelis and G. cricetidarum, which infect non-human hosts including amphibians, birds, rodents and marsupials. Giardia duodenalis is a species complex consisting of eight assemblages (A-H), with assemblages A and B the dominant assemblages in humans. Molecular studies to date on the zoonotic potential of Giardia in animals are problematic and are hampered by lack of concordance between loci. Livestock (cattle, sheep, goats and pigs) are predominantly infected with G. duodenalis assemblage E, which has recently been shown to be zoonotic, followed by assemblage A. In cats and dogs, assemblages A, B, C, D and F are commonly reported but relatively few studies have conducted molecular typing of humans and their pets and the results are contradictory with some studies support zoonotic transmission but the majority of studies suggesting separate transmission cycles. Giardia also infects a broad range of wildlife hosts and although much less well studied, host-adapted species as well as G. duodenalis assemblages (A-H) have been identified. Fish and other aquatic wildlife represent a source of infection for humans with Giardia via water contamination and/or consumption of undercooked fish and interestingly, assemblage B and A predominated in the two molecular studies conducted to date. Our current knowledge of the transmission dynamics of Giardia is still poor and the development of more discriminatory typing tools such as whole genome sequencing (WGS) of Giardia isolates is therefore essential.

Identification of Conserved Candidate Vaccine Antigens in the Surface Proteome of Giardia lamblia

Giardia lamblia, one of the most common protozoal infections of the human intestine, is an important worldwide cause of diarrheal disease, malabsorption, malnutrition, delayed cognitive development in children, and protracted postinfectious syndromes. Despite its medical importance, no human vaccine is available against giardiasis. A crude veterinary vaccine has been developed, and experimental vaccines based on expression of multiple variant-specific surface proteins have been reported, but poorly defined vaccine components and excessive antigen variability are problematic for pharmaceutical vaccine production. To expand the repertoire of antigen candidates for vaccines, we reasoned that surface proteins may provide an enriched source of such antigens since key host effectors, such as secretory IgA, can directly bind to such antigens in the intestinal lumen and interfere with epithelial attachment. Here, we have applied a proteomics approach to identify 23 novel surface antigens of G. lamblia that show >90% amino acid sequence identity between the two human-pathogenic genetic assemblages (A and B) of the parasite. Surface localization of a representative subset of these proteins was confirmed by immunostaining. Four selected proteins, uridine phosphorylase-like protein-1, protein 21.1 (GL50803_27925), α1-giardin, and α11-giardin, were subsequently produced in recombinant form and shown to be immunogenic in mice and G. lamblia-infected humans and confer protection against G. lamblia infection upon intranasal immunization in rodent models of giardiasis. These results demonstrate that identification of conserved surface antigens provides a powerful approach for overcoming a key rate-limiting step in the design and construction of an effective vaccine against giardiasis.

Genetic variation in metronidazole metabolism and oxidative stress pathways in clinical Giardia lamblia assemblage A and B isolates

Purpose: Treatment-refractory Giardia cases have increased rapidly within the last decade. No markers of resistance nor a standardized susceptibility test have been established yet, but several enzymes and their pathways have been associated with metronidazole (MTZ) resistant Giardia. Very limited data are available regarding genetic variation in these pathways. We aimed to investigate genetic variation in metabolic pathway genes proposed to be involved in MTZ resistance in recently acquired, cultured clinical isolates.

Methods: Whole genome sequencing of 12 assemblage A2 and 8 assemblage B isolates was done, to decipher genomic variation in Giardia. Twenty-nine genes were identified in a literature search and investigated for their single nucleotide variants (SNVs) in the coding/non-coding regions of the genes, either as amino acid changing (non-synonymous SNVs) or non-changing SNVs (synonymous).

Results: In Giardia assemblage B, several genes involved in MTZ activation or oxidative stress management were found to have higher numbers of non-synonymous SNVs (thioredoxin peroxidase, nitroreductase 1, ferredoxin 2, NADH oxidase, nitroreductase 2, alcohol dehydrogenase, ferredoxin 4 and ferredoxin 1) than the average variation. For Giardia assemblage A2, the highest genetic variability was found in the ferredoxin 2, ferredoxin 6 and in nicotinamide adenine dinucleotide phosphate (NADPH) oxidoreductase putative genes. SNVs found in the ferredoxins and nitroreductases were analyzed further by alignment and homology modeling. SNVs close to the iron-sulfur cluster binding sites in nitroreductase-1 and 2 and ferredoxin 2 and 4 could potentially affect protein function. Flavohemoprotein seems to be a variable-copy gene, due to higher, but variable coverage compared to other genes investigated.

Conclusion: In clinical Giardia isolates, genetic variability is common in important genes in the MTZ metabolizing pathway and in the management of oxidative and nitrosative stress and includes high numbers of non-synonymous SNVs. Some of the identified amino acid changes could potentially affect the respective proteins important in the MTZ metabolism.

Molecular genotyping, diversity studies and high-resolution molecular markers unveiled by microsatellites in Giardia duodenalis

BACKGROUND

Giardia duodenalis (synonyms G. lamblia and G. intestinalis) is an enteric protozoan parasite of a wide range of mammalian hosts, including humans and various domestic and wild animals. There is considerable genetic variability in G. duodenalis and isolates of this parasite have been divided into eight genetic assemblages. Microsatellites markers can be used to discriminate isolates with a high level of sensitivity. This study was conducted to identify and characterize genomic microsatellites (simple sequence repeats-SSRs), sequences of one- to six-nucleotide motifs repeated in tandem, present in the available genomes of G. duodenalis and to develop new markers that can serve as a tool for detection and for characterizing the genetic diversity of this parasite.

METHODOLOGY/ PRINCIPAL FINDINGS

For each genetic assemblage, polymorphism levels for the microsatellite markers were evaluated. After performing the analysis using the MISA and SciRoKo software, 1,853 simple sequence repeats (SSRs) were identified. In all the genomes, trinucleotide repeats were the most common class followed by tetranucleotide. Many of the SSR loci are assemblage-specific, and 36 SSR loci shared among all the genomes were identified. Together with hypothetical proteins, variant-specific surface proteins represented nearly half of the annotated SSR loci. The results regarding the most common repeat among the SSRs led us to infer that positive selection occurred to avoid frameshift mutations. Additionally, based on inter- and intra-genetic assemblages polymorphism analyses, we unveiled previously undetected genetic variation, indicating that the microsatellite markers we developed are useful molecular tools for epidemiological inferences based on population genetics patterns and processes.

CONCLUSIONS

There is increasing demand for the development of new molecular markers and for the characterization of pathogens at a higher resolution level. In this study, we present 60 G. duodenalis microsatellites markers that exhibited high polymerase chain reaction (PCR) amplification efficiency among the different genetic assemblages. Twenty of these markers presented nucleotide sequence polymorphisms and may be used as a genotyping tool. The monomorphic markers can be used for the detection of the parasite at the species and genetic assemblage level. These polymorphic markers revealed a genetic diversity that was previously undetectable, thus they can be considered valuable molecular tools for high resolution markers in future studies investigating Giardia and may also be used for epidemiological inferences based on populations genetics patterns and processes.

Robust and stable transcriptional repression in Giardia using CRISPRi

Giardia lamblia is a binucleate protistan parasite causing significant diarrheal disease worldwide. An inability to target Cas9 to both nuclei, combined with the lack of nonhomologous end joining and markers for positive selection, has stalled the adaptation of CRISPR/Cas9-mediated genetic tools for this widespread parasite. CRISPR interference (CRISPRi) is a modification of the CRISPR/Cas9 system that directs catalytically inactive Cas9 (dCas9) to target loci for stable transcriptional repression. Using a Giardia nuclear localization signal to target dCas9 to both nuclei, we developed efficient and stable CRISPRi-mediated transcriptional repression of exogenous and endogenous genes in Giardia. Specifically, CRISPRi knockdown of kinesin-2a and kinesin-13 causes severe flagellar length defects that mirror defects with morpholino knockdown. Knockdown of the ventral disk MBP protein also causes severe structural defects that are highly prevalent and persist in the population more than 5 d longer than defects associated with transient morpholino-based knockdown. By expressing two guide RNAs in tandem to simultaneously knock down kinesin-13 and MBP, we created a stable dual knockdown strain with both flagellar length and disk defects. The efficiency and simplicity of CRISPRi in polyploid Giardia allows rapid evaluation of knockdown phenotypes and highlights the utility of CRISPRi for emerging model systems.

Recent advances in the genomic and molecular biology of Giardia

Giardia duodenalis is the most common gastrointestinal protozoan parasite of humans and a significant contributor to the global burden of both diarrheal disease and post-infectious chronic disorders. Robust tools for analyzing gene function in this parasite have been developed and a range of genetic tools are now available. These together with public databases have provided insights on the function of different genes in Giardia. In this review we provide a current perspective on different molecular aspects of Giardia related to genomics, regulation of encystation, trophozoite transcriptional responses to physiological and xenobiotic (drug-induced) stress, and mechanisms of drug resistance. We also examine recent insights that have contributed to gain knowledge in the study of VSPs, antigenic variation, epigenetics, DNA repair and in the direct manipulation of gene function in Giardia, with a particular focus on the inducible Cre/loxP system.

Beaver Fever: Whole-Genome Characterization of Waterborne Outbreak and Sporadic Isolates To Study the Zoonotic Transmission of Giardiasis

Giardia duodenalis causes large numbers of gastrointestinal illness in humans. Its transmission through the contaminated surface water/wildlife intersect is significant, and the water-dwelling rodents beavers have been implicated as one important reservoir. To trace human infections to their source, we used genome techniques to characterize genetic relationships among 89 Giardia isolates from surface water, humans, and animals. Our study showed the presence of two previously described genetic assemblages, A and B, with mixed infections detected from isolates collected during outbreaks. Study findings also showed that while assemblage A could be divided into A1 and A2, A1 showed little genetic variation among animal and human hosts in isolates collected from across the globe. Assemblage B, the most common type found in the study surface water samples, was shown to be highly variable. Our study demonstrates that the beaver is a possible source of human infections from contaminated surface water, while acknowledging that theirs is only one role in the complex cycle of zoonotic spread. Mixes of parasite groups have been detected in waterborne outbreaks. More information on Giardia diversity and its evolution using genomics will further the understanding of the epidemiology of spread of this disease-causing protozoan.

Giardia causes the diarrheal disease known as giardiasis; transmission through contaminated surface water is common. The protozoan parasite’s genetic diversity has major implications for human health and epidemiology. To determine the extent of transmission from wildlife through surface water, we performed whole-genome sequencing (WGS) to characterize 89 Giardia duodenalis isolates from both outbreak and sporadic infections: 29 isolates from raw surface water, 38 from humans, and 22 from veterinary sources. Using single nucleotide variants (SNVs), combined with epidemiological data, relationships contributing to zoonotic transmission were described. Two assemblages, A and B, were identified in surface water, human, and veterinary isolates. Mixes of zoonotic assemblages A and B were seen in all the community waterborne outbreaks in British Columbia (BC), Canada, studied. Assemblage A was further subdivided into assemblages A1 and A2 based on the genetic variation observed. The A1 assemblage was highly clonal; isolates of surface water, human, and veterinary origins from Canada, United States, and New Zealand clustered together with minor variation, consistent with this being a panglobal zoonotic lineage. In contrast, assemblage B isolates were variable and consisted of several clonal lineages relating to waterborne outbreaks and geographic locations. Most human infection isolates in waterborne outbreaks clustered with isolates from surface water and beavers implicated to be outbreak sources by public health. In-depth outbreak analysis demonstrated that beavers can act as amplification hosts for human infections and can act as sources of surface water contamination. It is also known that other wild and domesticated animals, as well as humans, can be sources of waterborne giardiasis. This study demonstrates the utility of WGS in furthering our understanding of Giardia transmission dynamics at the water-human-animal interface.

IMPORTANCEGiardia duodenalis causes large numbers of gastrointestinal illness in humans. Its transmission through the contaminated surface water/wildlife intersect is significant, and the water-dwelling rodents beavers have been implicated as one important reservoir. To trace human infections to their source, we used genome techniques to characterize genetic relationships among 89 Giardia isolates from surface water, humans, and animals. Our study showed the presence of two previously described genetic assemblages, A and B, with mixed infections detected from isolates collected during outbreaks. Study findings also showed that while assemblage A could be divided into A1 and A2, A1 showed little genetic variation among animal and human hosts in isolates collected from across the globe. Assemblage B, the most common type found in the study surface water samples, was shown to be highly variable. Our study demonstrates that the beaver is a possible source of human infections from contaminated surface water, while acknowledging that theirs is only one role in the complex cycle of zoonotic spread. Mixes of parasite groups have been detected in waterborne outbreaks. More information on Giardia diversity and its evolution using genomics will further the understanding of the epidemiology of spread of this disease-causing protozoan.

The first multilocus genotype analysis of Giardia intestinalis in humans in the Czech Republic

To date, genotyping data on giardiasis have not been available in the Czech Republic. In this study, we characterized 47 human isolates of Giardia intestinalis from symptomatic as well as asymptomatic giardiasis cases. Genomic DNA from trophozoites was tested by PCR-sequence analysis at three loci (β-giardin, glutamate dehydrogenase and triose phosphate isomerase). Sequence analysis showed assemblages A and B in 41 (87.2%) and six (12.8%) isolates, respectively. Two of the 41 assemblage A samples were genotyped as sub-assemblage AI, and 39 were genotyped as sub-assemblage AII. Four previously identified multilocus genotypes (MLGs: AI-1, AII-1, AII-4 and AII-9) and six likely novel variations of MLGs were found. In agreement with previous studies, sequences from assemblage B isolates were characterized by a large genetic variability and by the presence of heterogeneous positions, which prevent the definition of MLGs. This study also investigated whether there was a relationship between the assemblage and clinical data (including drug resistance). However, due to the large number of genotypes and the relatively small number of samples, no significant associations with the clinical data were found.

Genetic variability and transcontinental sharing of Giardia duodenalis infrapopulations determined by glutamate dehydrogenase gene

Microevolutionary data of Giardia duodenalis sub-assemblages is a prerequisite for determining the invasion zoonotic patterns of the parasite. To infer transmission patterns that could not be differentiated by the phenotypic features, a population genetic investigation is crucial for the elucidation of the genetic structure of G. duodenalis among the continents. Forty G. duodenalis positive fecal samples were collected from different foci of Northwest Iran. The specimens were subjected to Trichrome staining and sucrose gradient flotation. DNA samples were extracted, amplified, and sequenced by targeting glutamate dehydrogenase (gdh) gene. The global gdh sequences of sub-assemblages AII and BIV retrieved from NCBI GenBank were analyzed to estimate diversity indices, neutrality indices, and gene migration tests. Sequencing analyses indicated various levels of genetic variability of sub-assemblages AII and BIV among the five continents. Sub-assemblage BIV had greater genetic variability (haplotype diversity: 0.975; nucleotide diversity: 0.04246) than sub-assemblage AII. The statistical Fst value demonstrated that the genetic structure of sub-assemblages AII and BIV are moderately differentiated between European-American populations (Fst: 0.05352-0.15182), whereas a significant differentiation was not seen among other geographical population pairs. We conclude that a high gene flow of G. duodenalis sub-assemblages AII and BIV is unequivocally sharing among the continents. The current findings strengthen our knowledge to assess the evolutionary patterns of G. duodenalis in endemic foci of the world and it will become the basis of public health policy to control human giardiasis.

Host specificity in the Giardia duodenalis species complex

Giardia duodenalis is a unicellular flagellated parasite that infects the gastrointestinal tract of a wide range of mammalian species, including humans. Investigations of protein and DNA polymorphisms revealed that G. duodenalis should be considered as a species complex, whose members, despite being morphologically indistinguishable, can be classified into eight groups, or Assemblages, separated by large genetic distances. Assemblages display various degree of host specificity, with Assemblages A and B occurring in humans and many other hosts, Assemblage C and D in canids, Assemblage E in hoofed animals, Assemblage F in cats, Assemblage G in rodents, and Assemblage H in pinnipeds. The factors determining host specificity are only partially understood, and clearly involve both the host and the parasite. Here, we review the results of in vitro and in vivo experiments, and clinical observations to highlight relevant biological and genetic differences between Assemblages, with a focus on human infection.

Giardia duodenalis Infections in Humans and Other Animals in China

Giardia duodenalis is an important zoonotic pathogen in both public and veterinary health, and has been genotyped into at least eight assemblages (A-H), each with a distinct host range. In recent years, this intestinal protozoan parasite has been identified widely in humans and various other animals, and has even been recorded in environmental contaminants. Along with whole genome sequencing of G. duodenalis, multilocus sequence typing is increasingly being used to characterize G. duodenalis isolates. Here, we review the epidemiology, genotyping, and subtyping of G. duodenalis from humans and a wide range of other animals, as well as from wastewater, in China.

Giardia: a pathogen or commensal for children in high-prevalence settings?

PURPOSE OF REVIEW

Giardia is a common intestinal parasite worldwide, and infection can be associated with clear and sometimes persistent symptomatology. However, in children in high-prevalence settings, it is not associated with or is perhaps even protective against acute diarrhea, and the association with long-term outcomes has been difficult to discern.

RECENT FINDINGS

Recent studies have made progress in helping us disentangle this apparent paradox. First, prospective, well-characterized cohort studies have added to the data on the association between Giardia and diarrhea in these settings and have further characterized associations between Giardia infection and nutrition, gut function, and growth. Second, animal models have further characterized the host response to Giardia and helped elucidate mechanisms by which Giardia could impair child development. Finally, new work has shed light on the heterogeneity of human Giardia strains, which may both explain discrepant findings in the literature and help guide higher-resolution analyses of this pathogen in the future.

SUMMARY

The true clinical impact of endemic pediatric giardiasis remains unclear, but recent prospective studies have confirmed a high prevalence of persistent, subclinical Giardia infections and associated growth shortfalls. Integrating how nutritional, microbial, metabolic, and pathogen-strain variables influence these outcomes could sharpen delineations between pathogenic and potentially beneficial attributes of this enigmatic parasite.

Molecular epidemiologic tools for waterborne pathogens Cryptosporidium spp. and Giardia duodenalis

Molecular diagnostic tools have played an important role in improving our understanding of the transmission of Cryptosporidium spp. and Giardia duodenalis, which are two of the most important waterborne parasites in industrialized nations. Genotyping tools are frequently used in the identification of host-adapted Cryptosporidium species and G. duodenalis assemblages, allowing the assessment of infection sources in humans and public health potential of parasites found in animals and the environment. In contrast, subtyping tools are more often used in case linkages, advanced tracking of infections sources, and assessment of disease burdens attributable to anthroponotic and zoonotic transmission. More recently, multilocus typing tools have been developed for population genetic characterizations of transmission dynamics and delineation of mechanisms for the emergence of virulent subtypes. With the recent development in next generation sequencing techniques, whole genome sequencing and comparative genomic analysis are increasingly used in characterizing Cryptosporidium spp. and G. duodenalis. The use of these tools in epidemiologic studies has identified significant differences in the transmission of Cryptosporidium spp. in humans between developing countries and industrialized nations, especially the role of zoonotic transmission in human infection. Geographic differences are also present in the distribution of G. duodenalis assemblages A and B in humans. In contrast, there is little evidence for widespread zoonotic transmission of giardiasis in both developing and industrialized countries. Differences in virulence have been identified among Cryptosporidium species and subtypes, and possibly between G. duodenalis assemblages A and B, and genetic recombination has been identified as one mechanism for the emergence of virulent C. hominis subtypes. These recent advances are providing insight into the epidemiology of waterborne protozoan parasites in both developing and developed countries.

Cryptosporidium and Giardia Infections in Children: A Review

Diarrheal disease remains the second leading cause of mortality in children in developing countries. Cryptosporidium is a leading cause and its importance stands to increase as rotavirus vaccine becomes used around the world. Cryptosporidium is particularly problematic in children younger than 2 years old and in the immunocompromised. Giardia lamblia is a common intestinal protozoan that is associated with diarrhea and, perhaps, growth faltering in impoverished settings. This review establishes the current prevalence of these infections in global settings and reviews current diagnosis and management approaches.

Curcumin alters the cytoskeleton and microtubule organization on trophozoites of Giardia lamblia

Giardia lamblia is a worldwide protozoan responsible for a significant number of intestinal infections. There are several drugs for the treatment of giardiasis, but they often cause side effects. Curcumin, a component of turmeric, has antigiardial activity; however, the molecular target and mechanism of antiproliferative activity are not clear. The effects of curcumin on cellular microtubules have been widely investigated. Since tubulin is the most abundant protein in the cytoskeleton of Giardia, to elucidate whether curcumin has activity against the microtubules of this parasite, we treated trophozoites with curcumin and the cells were analyzed by scanning electron microscopy and confocal microscopy. Curcumin inhibited Giardia proliferation and adhesion in a time-concentration-dependent mode. The higher inhibitory concentrations of curcumin (3 and 15μM) disrupted the cytoskeletal structures of trophozoites; the damage was evident on the ventral disk, flagella and in the caudal region, also the membrane was affected. The immunofluorescence images showed altered distribution of tubulin staining on ventral disk and flagella. Additionally, we found that curcumin caused a clear reduction of tubulin expression. By docking analysis and molecular dynamics we showed that curcumin has a high probability to bind at the interface of the tubulin dimer close to the vinblastine binding site. All the data presented indicate that curcumin may inhibit Giardia proliferation by perturbing microtubules.

Associations of Giardia lamblia assemblages with HIV infections and symptomatology: HIV virus and assemblage B were they born to each other?

Giardia lamblia is an intestinal parasite that has an extensive genetic variation among isolates. This species is divided into eight different assemblages (A-H), but only assemblages A and B have been associated with human infections. Studies on the associations of G. lamblia assemblages and symptoms have been done but were inconclusive. The aim of this study was to correlate G. lamblia assemblages with symptoms in patients with and without HIV/AIDS and its association with the CD4T cell count. The cross-sectional survey was conducted among patients attending the Evandro Chagas National Institute of Infectious Diseases (INI/FIOCRUZ) in Rio de Janeiro from January 2011 to February 2015. Thirty-eight of 65 microscopically positive stool samples for G. lamblia were from HIV positive patients and 27 were from HIV negative patients. Of the HIV infected patients, 19 (55.9%) were genotyped as assemblage B of which 9 (47.4%) had a CD4Tcell count below 200cells/mm³. In addition, we found a greater number of samples belonging to assemblage B in symptomatic cases (11 of 19; 57.9%). Our data suggest that assemblage B is very likely to be found in HIV infected patients and probably the lower CD4T count gives advantages for assemblage B replication. Furthermore, assemblage B seems to be associated with symptomatology, particularly abdominal pain, asthenia, diarrhea, fever, headache and myalgia. This study provides information on G. lamblia assemblages and symptoms in patients with and without HIV/AIDS virus and their association with CD4Tcell counts.

Giardia Colonizes and Encysts in High-Density Foci in the Murine Small Intestine

Giardia is a single-celled parasite causing significant diarrheal disease in several hundred million people worldwide. Due to limited access to the site of infection in the gastrointestinal tract, our understanding of the dynamics of Giardia infections in the host has remained limited and largely inferred from laboratory culture. To better understand Giardia physiology and colonization in the host, we developed imaging methods to quantify Giardia expressing bioluminescent physiological reporters in two relevant animal models. We discovered that parasites primarily colonize and encyst in the proximal small intestine in discrete, high-density foci. We also show that high parasite density contributes to encystation initiation.

Giardia lamblia is a highly prevalent yet understudied protistan parasite causing significant diarrheal disease worldwide. Hosts ingest Giardia cysts from contaminated sources. In the gastrointestinal tract, cysts excyst to become motile trophozoites, colonizing and attaching to the gut epithelium. Trophozoites later differentiate into infectious cysts that are excreted and contaminate the environment. Due to the limited accessibility of the gut, the temporospatial dynamics of giardiasis in the host are largely inferred from laboratory culture and thus may not mirror Giardia physiology in the host. Here, we have developed bioluminescent imaging (BLI) to directly interrogate and quantify the in vivo temporospatial dynamics of Giardia infection, thereby providing an improved murine model to evaluate anti-Giardia drugs. Using BLI, we determined that parasites primarily colonize the proximal small intestine nonuniformly in high-density foci. By imaging encystation-specific bioreporters, we show that encystation initiates shortly after inoculation and continues throughout the duration of infection. Encystation also initiates in high-density foci in the proximal small intestine, and high density contributes to the initiation of encystation in laboratory culture. We suggest that these high-density in vivo foci of colonizing and encysting Giardia likely result in localized disruption to the epithelium. This more accurate visualization of giardiasis redefines the dynamics of the in vivo Giardia life cycle, paving the way for future mechanistic studies of density-dependent parasitic processes in the host.

IMPORTANCEGiardia is a single-celled parasite causing significant diarrheal disease in several hundred million people worldwide. Due to limited access to the site of infection in the gastrointestinal tract, our understanding of the dynamics of Giardia infections in the host has remained limited and largely inferred from laboratory culture. To better understand Giardia physiology and colonization in the host, we developed imaging methods to quantify Giardia expressing bioluminescent physiological reporters in two relevant animal models. We discovered that parasites primarily colonize and encyst in the proximal small intestine in discrete, high-density foci. We also show that high parasite density contributes to encystation initiation.

The protein subunit of telomerase displays patterns of dynamic evolution and conservation across different metazoan taxa

BACKGROUND

Most animals employ telomerase, which consists of a catalytic subunit known as the telomerase reverse transcriptase (TERT) and an RNA template, to maintain telomere ends. Given the importance of TERT and telomere biology in core metazoan life history traits, like ageing and the control of somatic cell proliferation, we hypothesised that TERT would have patterns of sequence and regulatory evolution reflecting the diverse life histories across the Animal Kingdom.

RESULTS

We performed a complete investigation of the evolutionary history of TERT across animals. We show that although TERT is almost ubiquitous across Metazoa, it has undergone substantial sequence evolution within canonical motifs. Beyond the known canonical motifs, we also identify and compare regions that are highly variable between lineages, but show conservation within phyla. Recent data have highlighted the importance of alternative splice forms of TERT in non-canonical functions and although animals may share some conserved introns, we find that the selection of exons for alternative splicing appears to be highly variable, and regulation by alternative splicing appears to be a very dynamic feature of TERT evolution. We show that even within a closely related group of triclad flatworms, where alternative splicing of TERT was previously correlated with reproductive strategy, we observe highly diverse splicing patterns.

CONCLUSIONS

Our work establishes that the evolutionary history and structural evolution of TERT involves previously unappreciated levels of change and the emergence of lineage specific motifs. The sequence conservation we describe within phyla suggests that these new motifs likely serve essential biological functions of TERT, which along with changes in splicing, underpin diverse functions of TERT important for animal life histories.

Strong genetic structure revealed by multilocus patterns of variation in Giardia duodenalis isolates of patients from Galicia (NW-Iberian Peninsula)

We report a survey of genetic variation at three coding loci in Giardia duodenalis of assemblages A and B obtained from stool samples of patients from Santiago de Compostela (Galicia, NW-Iberian Peninsula). The mean pooled synonymous diversity for assemblage A was nearly five times lower than for assemblage B (0.77%±0.30% and 4.14%±1.65%, respectively). Synonymous variation in both assemblages was in mutation-drift equilibrium and an excess of low-frequency nonsynonymous variants suggested the action of purifying selection at the three loci. Differences between isolates contributed to 40% and 60% of total genetic variance in assemblages A and B, respectively, which revealed a significant genetic structure. These results, together with the lack of evidence for recombination, support that (i) Giardia assemblages A and B are in demographic equilibrium and behave as two genetically isolated populations, (ii) infections are initiated by a reduced number of individuals, which may be genetically diverse and even belong to different assemblages, and (iii) parasites reproduce clonally within the host. However, the observation of invariant loci in some isolates means that mechanisms for the homogenization of the genetic content of the two diploid nuclei in each individual must exist.