key: cord-010056-zfin4bko
authors: Mejia, Rojelio; Damania, Ashish; Jeun, Rebecca; Bryan, Patricia E.; Vargas, Paola; Juarez, Marisa; Cajal, Pamela S.; Nasser, Julio; Krolewiecki, Alejandro; Lefoulon, Emilie; Long, Courtney; Drake, Evan; Cimino, Rubén O.; Slatko, Barton
title: Impact of intestinal parasites on microbiota and cobalamin gene sequences: a pilot study
date: 2020-04-19
journal: Parasit Vectors
DOI: 10.1186/s13071-020-04073-7
sha: 
doc_id: 10056
cord_uid: zfin4bko

BACKGROUND: Approximately 30% of children worldwide are infected with gastrointestinal parasites. Depending on the species, parasites can disrupt intestinal bacterial microbiota affecting essential vitamin biosynthesis. METHODS: Stool samples were collected from 37 asymptomatic children from a previous cross-sectional Argentinian study. A multi-parallel real-time quantitative PCR was implemented for Ascaris lumbricoides, Ancylostoma duodenale, Necator americanus, Strongyloides stercoralis, Trichuris trichiura, Cryptosporidium spp., Entamoeba histolytica and Giardia duodenalis. In addition, whole-genome sequencing analysis was conducted for bacterial microbiota on all samples and analyzed using Livermore Metagenomic Analysis Toolkit and DIAMOND software. Separate analyses were carried out for uninfected, Giardia-only, Giardia + helminth co-infections, and helminth-only groups. RESULTS: For Giardia-only infected children compared to uninfected children, DNA sequencing data showed a decrease in microbiota biodiversity that correlated with increasing Giardia burden and was statistically significant using Shannonʼs alpha diversity (Giardia-only > 1 fg/µl 2.346; non-infected group 3.253, P = 0.0317). An increase in diversity was observed for helminth-only infections with a decrease in diversity for Giardia + helminth co-infections (P = 0.00178). In Giardia-only infections, microbiome taxonomy changed from Firmicutes towards increasing proportions of Prevotella, with the degree of change related to the intensity of infection compared to uninfected (P = 0.0317). The abundance of Prevotella bacteria was decreased in the helminths-only group but increased for Giardia + helminth co-infections (P = 0.0262). Metagenomic analysis determined cobalamin synthesis was decreased in the Giardia > 1 fg/µl group compared to both the Giardia < 1 fg/µl and the uninfected group (P = 0.0369). Giardia + helminth group also had a decrease in cobalamin CbiM genes from helminth-only infections (P = 0.000754). CONCLUSION: The study results may provide evidence for an effect of parasitic infections enabling the permissive growth of anaerobic bacteria such as Prevotella, suggesting an altered capacity of vitamin B12 (cobalamin) biosynthesis and potential impact on growth and development in children [Image: see text].

chronic diarrhea, severe anemia, and can lead to intestinal obstruction. Economically disadvantaged children have recurrent infections and malnutrition that may lead to growth and cognitive delays [4] . These children have more difficulties in school and, subsequently, in the job market [4] . The cycle continues when they remain in poverty and have their children [4] . The link between intestinal helminths and malnutrition leading to growth stunting and anemia has been found by others [5] [6] [7] [8] [9] [10] [11] [12] , and a Global Burden of Disease Study points to evidence that hookworm is a leading cause of anemia in resource-poor settings [13] . Valuable information from the Global Enteric Multicenter Study (GEMS) and studies related to returning travelers also reveals an unexpected global health impact caused by some protozoans, possibly including giardiasis [14, 15] . The Etiology, Risk Factors, and Interactions of Enteric Infections and Malnutrition and the Consequences for Child Health and Development (MAL-ED) study found an association with subclinical, non-diarrheal giardiasis and decreases in growth of children [16] . There are few studies attributing gut microbiome changes to giardiasis [17] [18] [19] and no published studies showing the impact on the human intestinal microbiome using multi-parallel real-time quantitative (qPCR) to detect the presence of Giardia and quantitating the burden of infection [20] . To date, most studies examining intestinal parasitism have not been able to study intestinal worms and protozoans simultaneously and to successfully dissect the relative contribution of each of the significant intestinal helminth or protozoan pathogens to specific diseases. The current state of diagnosing gastrointestinal parasites in resourcelimited areas uses the subjective method of microscopy. Depending on the parasite, single stool microscopy exam sensitivity ranges from 50-80% [21] . As a result, large numbers of infected children are not being diagnosed correctly and treated.

Gastrointestinal parasites may modulate intestinal inflammation, malabsorption, and microbiome changes [22] [23] [24] [25] [26] . The microbiome is associated with digestion, nutrition and health, but alterations in biodiversity can increase disease states and induce intestinal inflammation [27] . Animal studies further show changes in microbial diversity due to G. duodenalis [18, 28] . There is a lack of literature studying the relationship between GI parasitesʼ impact on the human intestinal microbiome [19] . These few studies have presented discordant results of the impact parasites have on human intestinal microbiota biodiversity [29, 30] .

The relationship between intestinal microbiota biodiversity may depend upon which specific parasite is present in the gut [30] . Since the burden of parasite infection is directly correlated to morbidity and disease, the intensity of infection may also impact the intestinal microbiota [30] . The qPCR quantitates the burden of helminths and protozoans, determining the correlation of burden to changes in intestinal microbiota biodiversity. Alterations in intestinal microbiota alter bacterial metabolites, such as vitamin B12 (cobalamin), reducing their availability for human use.

Interactions between parasites and intestinal microbiota may have a direct impact on child nutrition. Giardia duodenalis is known to cause malabsorption, steatorrhea, and diarrhea [31] , with preliminary studies finding improvements in vitamin B12 serum levels after treatment for giardiasis [32, 33] . Vitamin B12 is a crucial microbiota-derived co-enzyme for humans who cannot produce it [34] [35] [36] . As vitamin B12 production is unique to specific intestinal bacteria, alterations in intestinal microbiota could diminish vitamin B12 availability for human use [37] . Advances in next-generation DNA sequencing allow for precise taxonomic comparisons between intestinal microbiotas and can simultaneously be used to scan the intestinal microbiota meta-genome for the presence of functional genes necessary for the specific functions, like cobalamin synthesis. In this pilot study, parasite qPCR and next-generation DNA sequencing was used to explore whether quantitative burden of specific parasites (Giardia duodenalis and soil-transmitted helminths) influence the composition of intestinal microbial communities. Using vitamin B12 as a representative bacteria-generated nutrient, we analyzed bacterial metagenomes as a surrogate for changes in intestinal bacteria functions associated with intestinal parasitic infections. This was a preliminary study using a population with high prevalence of intestinal parasites. It is meant as an introductory for the future direction of our research.

This descriptive study aimed to determine the effect of G. duodenalis and other intestinal parasites on bacterial microbiota and subsequent cobalamin metagenomics. Samples were randomly selected from a previously published study using qPCR in peri-urban Argentina [38] . No previous antiparasitics or antibiotics were administered 3 months prior to the sample collection. Samples consisted of four groups: (i) a control group with no parasites detected by qPCR (uninfected); (ii) a Giardia-only infected group; (iii) a Giardia and helminth co-infection group; and (iv) a helminth-only infected group. Helminths included in this study were Ascaris lumbricoides, Ancylostoma duodenale, Necator americanus and Strongyloides stercoralis (Table 1) .

All stools were collected and immediately stored on ice and sealed in air-tight containers, frozen within 3 h, and DNA extracted using MP FastDNA Spin Kits for Soil (MP Biomedicals, Irvine, CA) [38] . This collection method was found to yield adequate microbiome data for study analysis [29] . Samples were processed in Salta, Argentina, using a modified bead-beating process described previously [38, 39] . DNA was processed for qPCR in Houston, Texas, as previously described [38, 39] . The intensity of parasite DNA was calculated using reference plasmids to create a standard curve, as previously reported [38] . The qPCR results from the previous Argentinian study were used for these analyses [38] .

DNA sequencing and library construction were performed at New England Biolabs (NEB). All the samples underwent removal of methylated DNA with NEBNext ® Microbiome enrichment kits (New England Biolabs, Ipswich, MA). No DNA size selection was made. DNA Library prep was done as per the manufacturerʼs protocol using 1 µg of sample DNA and eight cycles of PCR enrichment (NEBNext Ultra DNA Library Prep Kit for Illumina, Version 5.1, 9/17). Sequencing was performed using an Illumina NextSeq (Illumina, San Diego, CA) with paired or single ends, and 151 bp reads.

Fastq reads were filtered for quality at a minimum Phred score of 20 (99% accuracy) and a minimum length of 50 using Cutadapt version 1.8.3 [40] . Paired-end reads were interleaved using merge_fastq_reads_with_N_separator. Pl Perl script included with Livermore Metagenomics Analysis Toolkit (LMAT) software version 1.2.6 [41] . Fastq files were converted to fasta files using seqtk software version 1.0 (https ://githu b.com/lh3/seqtk ). Fasta files were processed by LMAT for taxonomic classification using LMAT kFull database. LMAT output text files were filtered for LMAT defined confidence score of 1 and minimum reads of 500 using tolineage.py script. Subsequently, output files were combined using merge_met-aphlan_tables.py script using Metaphlan [42] .

Alpha diversity was calculated using the Phyloseq R package [43] . Abundant different operational taxonomic units (OTU) among the four groups were identified using the LEfSe algorithm [44] . Most abundant bacteria for Lefse were run with a logarithmic LDA score threshold of 4.5 and other parameters set to default. Metagenomics analysis was performed using Diamond v0.8.4 using blastx mode with 90% minimum identity and e-value of 10 −5 against nr database fasta file [45] . Results from Diamond analysis were exported to Megan version 6 using daa-meganizer program [46] . GenInfo identifier to Interpro identifier mapping within Megan program was used to annotate the vitamin B12 synthesis gene [47] . STAMP software was used for statistical analysis pertaining to taxonomic and metagenomic differences [48] . Microbial attributes were derived from (LMAT) taxonomic output using the Megan program.

qPCR results were recorded for each patient as positive or negative, including the concentration of DNA (fg/ µl) for each parasite. All statistics were performed using Prism v. 7.0b (GraphPad, La Jolla, CA). Mann-Whitney and one-way ANOVA tests were used to compare two and multiple groups, respectively. Spearmanʼs rank test was used to correlate the Giardia DNA concentration to Shannonʼs alpha diversity, including the "uninfected" and "Giardia-only" group. Shannonʼs alpha diversity is a commonly reported diversity metric that weights the numbers of species by their relative evenness [49] . Proportion of targeted sequences was used in calculating OTU and metagenomics. Vitamin B12 data was calculated comparing specific genes to all other vitamin B12 synthesis genes. All statistical models used P-values less than 0.05 as significant.

Parasite DNA intensity with the intestinal parasite detection by stool qPCR as described by Cimino et al. [38] ( Table 1) .

Measuring bacterial diversity in relation from children with no Giardia DNA to increasing Giardia DNA intensity of burden (fg/µl) showed an inverse correlation of bacterial diversity to Giardia DNA (fg/µl) (Spearmanʼs r (18) = − 0.4781, P = 0.0447) (Fig. 1a) . From the scatter plot ( Fig. 1a ) a noticeable change in slope occurred at levels > 1 fg/µl. Since there are no guidelines for heavy Giardia burden, greater than 1 fg/µl was selected, and all further analyses validated this threshold. Giardia > 1 fg/µl group had a significantly lower biodiversity compared to the uninfected children using Shannonʼs alpha diversity (Giardia-only > 1 fg/µl: 2.346; uninfected group: 3.253, U (10) = 2, Z = − 2.1934, P = 0.0317 (Fig. 1b) . Giardia > 1 fg/µl had lower diversity than the helminth-only group U (17) = 0, Z = − 3.1096, P = 0.0003 (Fig. 1b) . Giardia > 1 fg/µl also had lower diversity compared to all other groups F (3, 25) = 9.982, P = 0.000166 (Fig. 1b) . Alpha diversity was also lower for the Giardia + helminth co-infected group compared to the helminth-only group U (19) = 7, Z = 2.958, P = 0.00178 (Fig. 1b) .

In Giardia-only infections, microbiome analysis data indicates decreased biodiversity in the infected parasite group compared to the non-infected group, a bias toward increased Prevotella, with the degree of change related to the intensity of infection (Fig. 2a) . Giardia-only children had significantly higher proportions of the genus Prevotella bacteria directly correlating to above 1 fg/µl Giardia DNA versus uninfected group U (10) = 2, Z = 2.1934, P = 0.0317 (Fig. 2a) . The helminth-only group with decreased Prevotella proportions relative to the Giardiaonly group U (25) = 41, Z = 2.0125, P = 0.0457, but similar to the control group (Fig. 2b) . Interestingly, the Giardia + helminth co-infected group had increased Prevotella proportions compared to the helminth-only group U (19) = 8, Z = − 2.8735, P = 0.00262 (Fig. 2b) . Prevotella copri was the major species in each group, Giardia (37%), helminths-only (17%), uninfected (22%), Giardia + helminths (36%) (data not shown).

All four groups had different bacteria genera as their most abundant microbiota. Giardia-infected children, including Giardia + helminth co-infected, had higher Bacteroidales, including Prevotella species (Fig. 3) . Specific cobalamin producing bacteria, such as Lactobacillus and Bifidobacterium, were found in higher abundance for the non-Giardia infected children (Fig. 3) .

Vitamin B12 InterPro identifier IPR002751 biosynthesis CbiM gene was used in the analysis. Children with Giardia > 1 fg/µl infections had fewer quantities of Fig. 1 a There was an inverse correlation of Giardia DNA to decreased Shannonʼs alpha diversity (Spearmanʼs r = − 0.4781, P = 0.0447); > 1 fg/ µl was selected as heavy Giardia DNA burden since a noticeable change in slope occurred after 1 fg/µl. b The mean of the Giardia > 1 fg/µl group (dotted line) has a significant decrease in bacterial diversity compared with the uninfected group (P = 0.0317) and all other groups (P = 0.000166). *P < 0.05 cobalamin DNA sequences than the Giardia < 1 fg/µl and uninfected group combined F (2, 15) = 4.145, P = 0.0369 (Fig. 4a) . Children with Giardia + helminth co-infections had a less proportion of vitamin B12 pathway DNA sequences, compared to helminth-group U (19) = 5, Z = 3.1271, P = 0.000754 (Figs. 4b, 5 ).

The association of G. duodenalis with microbiome diversity was observed in this study using whole-genome sequencing. Giardia duodenalis plays a prominent role, perhaps as its primary site of infection and replication is in the small intestine.

This study determined that children infected with Giardia DNA > 1 fg/µl are associated with decreased microbial diversity and increases of Prevotella. DNA levels of Giardia > 1 fg/µl implies that there are more parasites to alter the intestinal microbiome, and thus have a higher impact on intestinal bacterial species.

It is unclear whether the parasites are impacting the intestinal microbiota or external factors such as age, diet, or sex differences altering the intestinal microbiota and making the subjects more susceptible to enabling a G. duodenalis infection. Certain bacteria can permit G. duodenalis colonization, as evidenced from a mouse model study where the mouse intestinal microbiota (enteroaggregative Escherichia coli), independent of G. duodenalis infection, can promote inflammation, but together synergistically increased signals of intestinal injury [18] .

Vitamin B12 synthesis primarily occurs in anaerobes, including Bifidobacterium and Lactobacillus species [37, [50] [51] [52] [53] . These microorganisms may promote intestinal homeostasis and may protect against inflammatory diseases [54] [55] [56] [57] . Vitamin B12 is absorbed in the small intestines [58] [59] [60] while the majority of microbiota reside in the colon [58] , although, the small intestine is not sterile and does contain a robust microbiota that influences the absorption of vitamins [61] [62] [63] . Specific bacteria produce vitamin B12 [64] , and the children infected with G. duodenalis with DNA levels above 1 fg/µl may be unable to synthesize the required amounts of vitamin B12 for nutritional benefit. The group with less than 1 fg/µl Giardia DNA infections had equal vitamin B12 bacterial genes to the uninfected group, both being higher than the Giardia > 1 fg/µl group, likely showing that lower intensity of Giardia infections has similar effects as the uninfected group on vitamin B12 synthesis.

Evidence of the impact of G. duodenalis on the diversity and available micronutrients was also observed in the Giardia + helminth co-infection group, compared to the alterations of the microbiome seen in the helminthonly group. While the helminth-only group did not have changes in diversity or decreased cobalamin synthesis genes, a possible explanation is that most of these helminths reside in the colon and do not alter the microenvironments as does G. duodenalis.

A limitation of this cross-sectional study is the small sample size. However, the results are consistent, and this potential limitation can be viewed in terms of the sheer number, and fidelity of, using enriched microbe DNA for whole gene sequencing, producing millions of reads for analysis. This study also did not take into account Giardia assemblages that can have differing amount of pathogenicity. Another limitation is that serum vitamin B12 levels were not measures in children and thus the decrease of vitamin B12 genes in Giardia infected children could not be translated into a loss of vitamin B12.

In this study, there is a possible link as to why G. duodenalis and other parasites may cause growth and developmental delays in infected children. Giardia-infected children with > 1 fg/µl DNA concentrations were associated with less microbiome diversity, and a higher abundance of Prevotella associated with the diminished presence of cobalamin synthesis genes. The influence of G. duodenalis appears to be evident regardless of the presence of STH, and was associated with altered microbiome composition or function as measured by cobalamin synthesis. This descriptive study is a preliminary evidence for future studies looking at the cobalamin biosynthesis pathway. We are currently extending these results to additional populations where G. duodenalis and other intestinal parasite infections are endemic, determining Giardia assemblages, and confirming the biochemical analysis of the vitamin B12 pathway products in affected children.

Abbreviations DNA: Deoxyribonucleic acid; GEMS: Global enteric multicenter study; GI: Gastrointestinal; LMAT: Livermore metagenomics analysis toolkit; OTU: Operational taxonomic unit; qPCR: Multi-parallel real-time quantitative polymerase chain reaction; STH: Soil-transmitted helminths.

We wish to thank Don Comb, Rich Roberts, Tom Evans, Andy Gardner and Tilde Carlow (from NEB) for support of this project and the NEB DNA Sequencing Core (Laurie Mazzola, Danielle Fuchs, Joanna Bybee).

Authorsʼ contributions RM and AD were responsible for bioinformatics and statistical analysis. RM, RJ and PEB performed DNA extraction, qPCR, demethylation, and library preparation. BS, EL ED and CL were involved in the DNA sequencing sample reactions and preliminary analysis. PV, MJ, PSC, JN, AK and ROC collected patient samples in Argentina. RM, AD and BS wrote the manuscript. All authors read and approved the final manuscript.

Internal funds from New England Biolabs for DNA sequencing. Research funding support for RM was provided by the US Department of Health and Human Fig. 4 a Children infected with > 1 fg/µl Giardia DNA concentrations had decreased cobalamin biosynthesis genes compared to the combined Giardia < 1 fg/µl and uninfected groups (P = 0.0369) with b the Giardia + helminth group having less cobalamin biosynthesis genes than the helminth-only infections (P = 0.000754). *P < 0.05 

Global epidemiology, ecology and control of soil-transmitted helminth infections

Soil-transmitted helminth infections: ascariasis, trichuriasis, and hookworm

Waterborne transmission of protozoan parasites: review of worldwide outbreaks-an update

Forgotten people, forgotten diseases: the neglected tropical diseases and their impact on global health and development

The effect of deworming on growth in 1-year-old children living in a soil-transmitted helminth-endemic area of Peru: a randomized controlled trial

Impacts of dietary intake and helminth infection on diversity in growth among schoolchildren in rural south China: a 4-year longitudinal study

Soil-transmitted helminth infection, skin infection, anaemia, and growth retardation in schoolchildren of Taveuni Island. Fiji

Treatment for intestinal helminth infection. Review needed to take account of all relevant evidence, not only effects on growth and cognitive performance

Effects of treatment for intestinal helminth infection on growth and cognitive performance in children: systematic review of randomised trials

Improvements of growth, appetite, and physical activity in helminthinfected schoolboys 6 months after single dose of albendazole

Improvements in appetite and growth in helminth-infected schoolboys 3 and 7 weeks after a single dose of pyrantel pamoate

Influence of helminth parasite exposure and strategic application of anthelmintics on the development of immunity and growth of swine

A systematic analysis of global anemia burden from 1990 to

Burden and aetiology of diarrhoeal disease in infants and young children in developing countries (the Global Enteric Multicenter Study, GEMS): a prospective, case-control study

Enteropathogens and chronic illness in returning travelers

Use of quantitative molecular diagnostic methods to investigate the effect of enteropathogen infections on linear growth in children in lowresource settings: longitudinal analysis of results from the MAL-ED cohort study

Intestinal parasitic infection alters bacterial gut microbiota in children

Cross-modulation of pathogen-specific pathways enhances malnutrition during enteric co-infection with Giardia lamblia and enteroaggregative Escherichia coli

Gut microbiota related to Giardia duodenalis, Entamoeba spp. and Blastocystis hominis infections in humans from Côte dʼIvoire

A prospective longitudinal cohort to investigate the effects of early life giardiasis on growth and all cause diarrhea

Nonutility of routine testing of stool for ova and parasites in a tertiary care Canadian centre

The impact of anthelmintic treatment on human gut microbiota based on cross-sectional and pre-and postdeworming comparisons in western Kenya

Co-infection with Ascaris lumbricoides modulates protective immune responses against Giardia duodenalis in school Venezuelan rural children

Comparison of cytokine responses in Ecuadorian children infected with Giardia, Ascaris, or both parasites

Giardia intestinalis and fructose malabsorption: a frequent association

Determinants and impact of Giardia infection in the first 2 years of life in the MAL-ED birth cohort

Iron fortification adversely affects the gut microbiome, increases pathogen abundance and induces intestinal inflammation in Kenyan infants

Giardia alters commensal microbial diversity throughout the murine gut

Patent human infections with the whipworm, Trichuris trichiura, are not associated with alterations in the faecal microbiota

Helminth colonization is associated with increased diversity of the gut microbiota

Malabsorption associated with Giardia lamblia infestation

Vitamin B12 and folic acid in children with intestinal parasitic infection

Impact of protein malnutrition on histological parameters of experimentally infected animals with Giardia lamblia

Microbial production of vitamin B12

Cobalamin (coenzyme B12): synthesis and biological significance

Evolution of coenzyme B12 synthesis among enteric bacteria: evidence for loss and reacquisition of a multigene complex

Bacteria as vitamin suppliers to their host: a gut microbiota perspective

Identification of human intestinal parasites affecting an asymptomatic peri-urban Argentinian population using multi-parallel quantitative real-time polymerase chain reaction

A novel, multi-parallel, real-time polymerase chain reaction approach for eight gastrointestinal parasites provides improved diagnostic capabilities to resource-limited at-risk populations

Trimmomatic: a flexible trimmer for illumina sequence data

Scalable metagenomic taxonomy classification using a reference genome database

Metagenomic microbial community profiling using unique cladespecific marker genes

Phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data

Metagenomic biomarker discovery and explanation

Fast and sensitive protein alignment using DIAMOND

MEGAN analysis of metagenomic data

The InterPro protein families database: the classification resource after 15 years

Identifying biologically relevant differences between metagenomic communities

Drivers of microbiome biodiversity: a review of general rules, feces, and ignorance

Microbiomic analysis of the bifidobacterial population in the human distal gut

Exploring the diversity of the bifidobacterial population in the human intestinal tract

Host-microbe interactions that facilitate gut colonization by commensal bifidobacteria

Diversity of bifidobacteria within the infant gut microbiota

Interaction of probiotics and pathogens-benefits to human health?

Probiotics, immune function, infection and inflammation: a review of the evidence from studies conducted in humans

Systematic review of randomized controlled trials of probiotics, prebiotics, and synbiotics in inflammatory bowel disease

Low-complexity microbiota in the duodenum of children with newly diagnosed ulcerative colitis

Vitamin B12 in meat and dairy products

Intestinal absorption of vitamin B12 in man

Intestinal absorption of vitamin B12 in humans as studies by isotope technic

Bacterial communities in the small intestine respond differently to those in the caecum and colon in mice fed low-and high-fat diets

Small intestinal bacterial overgrowth and irritable bowel syndrome: a bridge between functional organic dichotomy

Small intestinal bacterial overgrowth and other intestinal disorders

Systematic genome assessment of B-vitamin biosynthesis suggests co-operation among gut microbes

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations

• fast, convenient online submission 

Ready to submit your research ? Choose BMC and benefit from: