A shotgun metagenomic analysis of the fecal microbiome in humans infected with Giardia duodenalis

Abstract

Abstract Background The mechanisms underlying the clinical outcome disparity during human infection with Giardia duodenalis are still unclear. In recent years, evidence has pointed to the roles of host factors as well as parasite’s genetic heterogeneity as major contributing factors in the development of symptomatic human giardiasis. However, it remains contested as to how only a small fraction of individuals infected with G. duodenalis develop clinical gastrointestinal manifestations, whereas the majority of infected individuals remain asymptomatic. Here, we demonstrate that diversity in the fecal microbiome correlates with the clinical outcome of human giardiasis. Methods The genetic heterogeneity of G. duodenalis clinical isolates from human subjects with asymptomatic and symptomatic giardiasis was determined using a multilocus analysis approach. We also assessed the genetic proximity of G. duodenalis isolates by constructing phylogenetic trees using the maximum likelihood. Total genomic DNA (gDNA) from fecal specimens was utilized to construct DNA libraries, followed by performing paired-end sequencing using the HiSeq X platform. The Kraken2-generated, filtered FASTQ files were assigned to microbial metabolic pathways and functions using HUMAnN 3.04 and the UniRef90 diamond annotated full reference database (version 201901b). Results from HUMAnN for each sample were evaluated for differences among the biological groups using the Kruskal–Wallis non-parametric test with a post hoc Dunn test. Results We found that a total of 8/11 (72.73%) human subjects were infected with assemblage A (sub-assemblage AII) of G. duodenalis, whereas 3/11 (27.27%) human subjects in the current study were infected with assemblage B of the parasite. We also found that the parasite’s genetic diversity was not associated with the clinical outcome of the infection. Further phylogenetic analysis based on the tpi and gdh loci indicated that those clinical isolates belonging to assemblage A of G. duodenalis subjects clustered compactly together in a monophyletic clade despite being isolated from human subjects with asymptomatic and symptomatic human giardiasis. Using a metagenomic shotgun sequencing approach, we observed that infected individuals with asymptomatic and symptomatic giardiasis represented distinctive microbial diversity profiles, and that both were distinguishable from the profiles of healthy volunteers. Conclusions These findings identify a potential association between host microbiome disparity with the development of clinical disease during human giardiasis, and may provide insights into the mechanisms by which the parasite induces pathological changes in the gut. These observations may also lead to the development of novel selective therapeutic targets for preventing human enteric microbial infections. Graphical Abstract