Extremely diversified haplotypes observed among assemblage B population of Giardia intestinalis in Kenya

Assemblages and Subassemblages of Giardia duodenalis in Rural Western, Kenya: Association with Sources, Signs, and Symptoms

Background

Giardia duodenalis causes sporadic or epidemic infections in humans. The parasite comprises assemblages A-H with A and B subdivided further into AI-IV and BI-IV subassemblages. Attempts aimed at linking these genotypes with sources and gastrointestinal manifestations of the infection are largely unexplored in rural communities.

Methods

In this cross-sectional study, G. duodenalis infection was genotyped and associated with sources, and gastrointestinal signs and symptoms of the disease among residents of Busia County, a rural setting in western Kenya. Demographic and clinical information were captured using standardized forms. Stool specimens were obtained from the patients and used for genotyping at glutamate dehydrogenase and triose-phosphate isomerase loci using the polymerase chain reaction and restriction fragment length polymorphism.

Results

Assemblage B (63.6%) was the most prevalent G. duodenalis infection, while A (20.5%) and mixed A/B (15.9%) were also detected. Among the subassemblages, AI (5.7%), AII (8.0%), AIII (3.4), BIII (30.7%), and BIV (17.0%) were diagnosed including the mixed AII/BIII (15.9%), BIII/BIV (15.9%), AI/AIII (2.3%), and AI/AII (1.1%) infections. Binary logistic regression indicated associations for assemblage A with stomach upset, history of nitroimidazole treatment, and residing in a homestead with cattle and B with age < 18 years, history of eating outdoors, vomiting, steatorrhea, and residing in a homestead with cattle, goats, and poultry (p < 0.05 for all). Among the subassemblages, associations were found for AI with residing in a homestead having cattle and history of nitroimidazole treatment, BIII with residing in a homestead having cattle and poultry, and BIV with steatorrhea (p < 0.05 for all). Altogether, this study illustrates that G. duodenalis assemblage B and subassemblage BIII are the most predominant and are linked to age < 18 years, gastrointestinal manifestations, and living in a homestead with domestic ruminants and poultry.

Conclusion

Targeted mass prophylactic treatment of domestic animals and utilization of gastrointestinal presentations, age < 18 years, and a history of nitroimidazole use are useful in the diagnosis and prevention of giardiasis among residents of rural communities.

Genetic diversity in the metronidazole metabolism genes nitroreductases and pyruvate ferredoxin oxidoreductases in susceptible and refractory clinical samples of Giardia lamblia

The effectiveness of metronidazole against the tetraploid intestinal parasite Giardia lamblia is dependent on its activation/inactivation within the cytoplasm. There are several activating enzymes, including pyruvate ferredoxin reductase (PFOR) and nitroreductase (NR) 1 which metabolize metronidazole into toxic forms, while NR2 on the other hand inactivates it. Metronidazole treatment failures have been increasing rapidly over the last decade, indicating genetic resistance mechanisms. Analyzing genetic variation in the PFOR and NR genes in susceptible and refractory Giardia isolates may help identify potential markers of resistance. Full length PFOR1, PFOR2, NR1 and NR2 genes from clinical culturable isolates and non-cultured clinical Giardia assemblage B samples were cloned, sequenced and single nucleotide variants (SNVs) were analyzed to assess genetic diversity and alleles. A similar ratio of amino acid changing SNVs per gene length was found for the NRs; 4.2% for NR1 and 6.4% for NR2, while the PFOR1 and PFOR2 genes had less variability with a ratio of 1.1% and 1.6%, respectively. One of the samples from a refractory case had a nonsense mutation which caused a truncated NR1 gene in one out of six alleles. Further, we found three NR2 alleles with frameshift mutations, possibly causing a truncated protein in two susceptible isolates. One of these isolates was homozygous for the affected NR2 allele. Three nsSNVs with potential for affecting protein function were found in the ferredoxin domain of the PFOR2 gene. The considerable variation and discovery of mutations possibly causing dysfunctional NR proteins in clinical Giardia assemblage B isolates, reveal a potential for genetic link to metronidazole susceptibility and resistance.

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.

Giardia duodenalis in Wildlife: Exploring Genotype Diversity in Italy and across Europe

Fragmented data are so far available on genotype diversity of G. duodenalis in wildlife in different countries in Europe, in particular, in Italy. In the present study, G. duodenalis sequences obtained from different Italian wild animals [12 porcupines (Hystrix cristata), 4 wild boars (Sus scrofa), 1 wolf (Canis lupus italicus), 6 Alpine chamois (Rupicapra rupicapra rupicapra)] were compared with those available from wild host species in Europe to add new data on the geographic distribution of Giardia assemblages/sub-assemblages and their transmission patterns among natural hosts. Thirty-eight sequences were obtained by MLG analysis (SSU-rRNA, bg, gdh, and tpi genes) and subsequently compared by phylogenetic and network analyses with those from wild species monitored in the last decades in Europe. The results revealed the presence of potentially zoonotic (A-AI, A-AII from wild boar; B from porcupine) and host-adapted (D from wolf; E, A-AIII from chamois) assemblages and sub-assemblages and represent the first report for Italian wild boar. The analysis did not find any evidence of spatial or host segregation for specific genetic variants, mostly shared between different hosts from different European countries. However, conflicting evidence was found in genotypic assignment, advocating for data improvement and new genomic approaches.

The sexual side of parasitic protists

Much of the vast evolutionary landscape occupied by Eukaryotes is dominated by protists. Though parasitism has arisen in many lineages, there are three main groups of parasitic protists of relevance to human and livestock health: the Apicomplexa, including the malaria parasite Plasmodium and coccidian pathogens of livestock such as Eimeria; the excavate flagellates, encompassing a diverse range of protist pathogens including trypanosomes, Leishmania, Giardia and Trichomonas; and the Amoebozoa, including pathogenic amoebae such as Entamoeba. These three groups represent separate, deep branches of the eukaryote tree, underlining their divergent evolutionary histories. Here, I explore what is known about sex in these three main groups of parasitic protists.

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.

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.