ORIGINAL RESEARCH article
Relationships Between Gastrointestinal Parasite Infections and the Fecal Microbiome in Free-Ranging Western Lowland Gorillas
- 1Department of Pathology and Parasitology, Faculty of Veterinary Medicine, University of Veterinary and Pharmaceutical Sciences Brno, Brno, Czechia
- 2Institute of Vertebrate Biology, Czech Academy of Sciences, Brno, Czechia
- 3Liberec Zoo, Liberec, Czechia
- 4Institute of Parasitology, Biology Centre of the Czech Academy of Sciences, České Budějovice, Czechia
- 5Central European Institute for Technology (CEITEC), University of Veterinary and Pharmaceutical Sciences Brno, Brno, Czechia
- 6WWF, Dzanga Sangha Protected Areas, Bangui, Central African Republic
- 7Department of Animal and Range Sciences, Montana State University, Bozeman, MT, United States
- 8J. Craig Venter Institute, Rockville, MD, United States
- 9Carl R. Woese Institute for Genomic Biology, University of Illinois at Urbana-Champaign, Urbana, IL, United States
- 10Department of Microbiology, University of Illinois at Urbana-Champaign, Urbana, IL, United States
- 11Department of Anthropology, University of Illinois at Urbana-Champaign, Urbana, IL, United States
- 12J. Craig Venter Institute, La Jolla, CA, United States
- 13Department of Anthropology, University of Colorado at Boulder, Boulder, CO, United States
- 14Department of Animal Science, University of Minnesota, St Paul, MN, United States
Relationships between gastrointestinal parasites (GIPs) and the gastrointestinal microbiome (GIM) are widely discussed topics across mammalian species due to their possible impact on the host's health. GIPs may change the environment determining alterations in GIM composition. We evaluated the associations between GIP infections and fecal microbiome composition in two habituated and two unhabituated groups of wild western lowland gorillas (Gorilla g. gorilla) from Dzanga Sangha Protected Areas, Central African Republic. We examined 43 fecal samples for GIPs and quantified strongylid nematodes. We characterized fecal microbiome composition through 454 pyrosequencing of the V1-V3 region of the bacterial 16S rRNA gene. Entamoeba spp. infections were associated with significant differences in abundances of bacterial taxa that likely play important roles in nutrition and metabolism for the host, besides being characteristic members of the gorilla gut microbiome. We did not observe any relationships between relative abundances of several bacterial taxa and strongylid egg counts. Based on our findings, we suggest that there is a significant relationship between fecal microbiome and Entamoeba infection in wild gorillas. This study contributes to the overall knowledge about factors involved in modulating GIM communities in great apes.
The mammalian gastrointestinal microbiome (GIM) is essential in providing access to nutrients, regulating epithelial development, and shaping immunity (Berg, 1996; Isolauri et al., 2004; Eckburg et al., 2005; Round and Mazmanian, 2009; Consortium, 2012). Thus, determining the external and internal factors that impact the composition, metabolic activity and dynamics of the GIM are critical to understand host health (Mackie et al., 1999; Isolauri et al., 2004).
The impact of gastrointestinal parasite (GIP) infections on the composition and activity of the GIM has received increasing attention and several studies have been conducted in humans and animals (Berrilli et al., 2012; Cantacessi et al., 2014; Kreisinger et al., 2015; Morton et al., 2015; Šlapeta et al., 2015; Audebert et al., 2016; Zaiss and Harris, 2016). Great apes, as closest relatives to humans (Kumar et al., 2005), are useful models for microbiome studies, enabling microbiome patterns to be studied in the contexts of ecology, health, and evolution (Gomez et al., 2015, 2016b; Moeller et al., 2016; Vlčková et al., 2016). Previous studies on GIM of western lowland gorillas Gorilla gorilla gorilla showed that feeding behavior, geographical range and increased anthropogenic pressure may be important modulators of the gorilla GIM. Moreover, the GIMs of two gorilla species G. gorilla and G. beringei exhibit significantly different patterns, but the GIMs converge when hosts face similar dietary constraints, associated with low fruit availability in their habitats suggesting that in primates dietary constraints triggered during their adaptive radiation were potential factors behind the species-specific GIM patterns observed today (Ochman et al., 2010; Gomez et al., 2015, 2016b). Moreover, GIM communities of sympatric chimpanzees and lowland gorillas have converged in terms of community composition (Moeller et al., 2015). The overall composition of gorilla GIP corresponds to that of other non-ruminating herbivores with a dominance of entodiniomorphid ciliates (which are in fact commensals), strongylids and anoplocephalid tapeworms, especially in mountain gorillas (Kalousová, 2015). GIP infections of western lowland gorilla have been studied using traditional and advanced techniques at several sites, with most of the studies conducted in Dzanga Sangha Protected Areas (Landsoud-Soukate et al., 1995; Lilly et al., 2002; Masi et al., 2012; Sak et al., 2013; Hasegawa et al., 2014; Pafčo et al., 2017; Vlčková et al., 2018).
To our knowledge, there have been no studies focusing on the relationship between GIP infections and the GIM of free-ranging great apes yet. Human societies have always been challenged by infections by helminth and protozoan parasites and GIP have evolved with humans throughout history (Cox, 2002). It is hypothesized that a reduced exposure to potentially pathogenic organisms, including parasites, in developed countries may result in insufficient stimulation of the human immune system, leading to increased incidence of autoimmune, metabolic and allergic diseases in these populations (Reddy, 2010). Industrialized human populations are hypothesized to have lost a significant number of GIM signatures in comparison with non-human primates and traditional human groups such as hunter-gatherer populations (Moeller et al., 2014; Obregon-Tito et al., 2015; Gomez et al., 2016a). These losses are manifested through decreased abundances of taxa like Treponema and Prevotella and potentially linked to traditional dietary behaviors (i.e., more fibrous food) (Gomez et al., 2016a).
The metabolic products of the resident GIM may strongly interfere with the survival and the physiology of many parasites and consequently, influence the outcome of many parasitic infections. Conversely, GIPs, both protozoans and helminths, constantly excrete and secrete molecules (e.g., carbon dioxide, nitrogenous compounds, enzymes) that may change the gut environment, potentially leading to alterations in GIM composition (Berrilli et al., 2012). For instance, some helminths are able to significantly reduce inflammation during colitis (Elliott et al., 2007), potentially modulating alterations in colon epithelial barrier function (Su et al., 2011). Nematodes can alter the gut habitat, which may lead to alterations in bacterial composition and abundance (Walk et al., 2010). These changes can be either beneficial or harmful to the host. Wu et al. (2012) detected increases in the levels of Campylobacter in pigs infected with Trichuris suis, resulting in campylobacteriosis and decreases in Fibrobacter and Ruminococcus, with subsequent disruptions of fibrolytic activity. Various GIPs like Blastocystis hominis, Giardia spp., Entamoeba histolytica, Dientamoeba fragilis, and Trichinella spp. have been reported to contribute to the progress of irritable bowel syndrome in humans (Mohammadi et al., 2015).
Nematodes such as T. suis and their excretory/secretory products (multiple complex molecules including many (glyco)proteins, which may function as proteases or pore-forming proteins that could affect epithelial integrity) can suppress innate and adaptive pro-inflammatory immune responses (Whelan et al., 2012; Hiemstra et al., 2014), making them potentially suitable for the treatment of inflammatory bowel disease (Moreels and Pelckmans, 2005; Summers et al., 2005). Likewise, hookworm infection (e.g. Necator americanus) has been hypothesized to decrease gluten sensitivity and could serve as a potential treatment for celiac disease (Mohammadi et al., 2015). Moreover, amoebas from the genus Entamoeba significantly impact GIM composition and diversity in humans (Morton et al., 2015) and Giardia intestinalis causes alterations in the bacterial composition of the upper digestive tract (Tandon et al., 1977; Tomkins et al., 1978), with the GIM potentially influencing the clinical outcomes of giardiasis. The impact of GIPs on GIM composition and diversity have been recorded also in the case of Giardia infection in dogs (Šlapeta et al., 2015), helminths in wild rodents (Kreisinger et al., 2015), and Heligmosomoides polygyrus in mice (Walk et al., 2010). Nevertheless, the manner in which GIPs interact with resident gut microbes of primates remains poorly understood (Berrilli et al., 2012). It is assumed that a specific GIM predisposes an individual to parasite colonization, namely Entamoeba (Bracha et al., 1982; Galván-Moroyoqui et al., 2008). Burgess and Petri (2016) observed that GIM is a significant factor influencing the clinical presentation and outcome of E. histolytica infections. The GIM may significantly affect the virulence of the Entamoeba itself, its ability to colonize the gut, and the host's immune responses at baseline, and during amebiasis (Burgess and Petri, 2016; Ngobeni et al., 2017).
As a part of the long-term non-invasive monitoring of GIP infections and fecal microbiome/GIM of western lowland gorillas (Gorilla g. gorilla) in Dzanga Sangha Protected Areas (DSPA), Central African Republic (Sak et al., 2013; Hasegawa et al., 2014; Gomez et al., 2015, 2016b; Vlčková et al., 2016), we evaluate the associations between GIP infections and composition of fecal microbiome of free-ranging western lowland gorillas. We particularly focused on strongylid nematodes (Necator, Oesophagostomum, Mammomonogamus), which are highly prevalent in studied gorillas (Kalousová, 2013; Hasegawa et al., 2014; Kalousová et al., 2016; Červená et al., 2017). We included Mammomonogamus into the study, even though it is not a GIP. As the adults of this strongylid live in the respiratory tract of the host and only eggs are shed in the gastrointestinal tract, we did not expect the impact of Mammomonogamus on fecal microbiome. The following research questions were addressed: (i) Does the fecal microbial diversity differ between individuals positive and negative for particular GIP? (ii) Does the abundance of specific bacterial taxa differ between individuals positive and negative for particular GIP? (iii) Is there a relationship between fecal microbial diversity and surrogate measures of intensity of parasite strongylid infection (egg counts)? (iv) Is there a relationship between abundances of specific bacterial taxa and surrogate measures of intensity of parasite strongylid infection (egg counts)?
Materials and Methods
Study Site and Sample Collection
Fecal samples of western lowland gorillas (G. g. gorilla) were collected at two research sites: Mongambe (2°55′N, 16°23'E) and Bai Hokou (2°50′N, 16°28′E) in the Dzanga sector of the Dzanga-Ndoki National Park (part of the DSPA complex), Central African Republic (Figure S1). The samples were collected in 2011 during June through September, which corresponds to a rainy season with higher fruit availability. Samples from two habituated gorilla groups: H1 (9) and H2 (12) were collected at Bai Hokou and Mongambe, respectively. Fecal samples were obtained non-invasively immediately after defecation during focal follows of individual habituated gorillas. Samples from two unhabituated gorilla groups were collected only once from the gorilla nest sites to avoid duplication: U1 (11) near Bai Hokou and U2 (11) between two research camps. A number of samples and the information about individual gorillas and appropriate gorilla group are listed in Table S1. Data regarding sex and age of individual gorillas was available only for individuals from habituated groups and since we worked also with unhabituated gorillas, this information was excluded. Nest site samples were collected early in the morning, ensuring that only relatively fresh feces were obtained within 3 h from the time the gorillas had left the nests as assessed by expert trackers. The samples were collected from within the core of the feces and placed in two collection tubes, containing RNAlater® (Thermo Fisher Scientific, Waltham, MA) and 4% formaldehyde for microbiome and parasitological analyses, respectively. The samples were kept at room temperature for up to a month before transport to the Department of Pathology and Parasitology, University of Veterinary and Pharmaceutical Sciences Brno, Czech Republic (UVPS), where the samples in RNAlater® were stored at -20°C and samples in 4% formaldehyde were stored at room temperature until coprological examination. Extractions of genomic DNA were performed at the Carl R. Woese Institute for Genomic Biology, University of Illinois at Urbana-Champaign, USA. Resulting genomic DNA was shipped to the J. Craig Venter Institute, Rockville, MD, USA where DNA sequencing analyses were performed. The research adhered to the legal requirements of the Central African Republic and research protocol of DSPA. Importation of the samples to the EU was approved by the State Veterinary Authority of the Czech Republic.
Each sample was homogenized and strained through a sieve into Falcon conical tubes (50 ml). Samples were diluted with water up to 50 ml volume and centrifuged for 5 min at 2,000 rpm. The sediment was weighed and re-suspended to 10 ml with 4% formaldehyde. For coprological examinations, Sheather's flotation with modified sugar solution was used (specific gravity 1.33) (Sheather, 1923).
A modified simple sedimentation method (Kassai, 1999; Zajac and Conboy, 2012) was used for thorough examination and quantification of strongylid eggs. Protozoans are difficult to quantify and most of those occurring in studied gorillas are not pathogenic or even mutualistic. Strongylid nematodes with pathogenic potential are usually present in most of, if not all, examined fecal samples, so further quantification was needed in order to evaluate their impact on the GIM. Each fecal suspension (2 ml) was placed into an Eppendorf tube and centrifuged for 2 min at 1,500 rpm. The supernatant was discarded, and the remaining sediment transferred onto a microscopic slide using a micropipette, covered with a cover-glass and examined using a light microscope. The strongylid eggs present on the slide were counted. This procedure was repeated until all the sediment (from 2 ml of the solution) was examined. The number of eggs per gram of sediment (EPG) was calculated using following equation: EPG = N/(m/5) where N = number of eggs and m = weight of the whole sediment.
We reported parasite presence/absence and for strongylid nematodes, also the intensity of infection expressed as EPG. Although the number of parasite eggs shed in feces may not be linearly correlated with the intensity of infection (i.e., the number of individuals of a particular parasite species in a single infected host Bush et al., 1997) due to high variability in individual and temporal egg output (Anderson and Schad, 1985; Warnick, 1992), several studies did find a linear relationship between egg counts and the number of adult worms in the hosts (Cabaret et al., 1998; Seivwright et al., 2004), and egg output is routinely used to characterize parasite infections in domestic animals (Piekarska et al., 2013; Nielsen et al., 2014), humans (Hodges et al., 2011; Mekonnen et al., 2013), and also in wild animals, including non-human primates (Müller-Graf et al., 1996; Chapman et al., 2006; González-Hernández et al., 2014; Burgunder et al., 2017). Thus egg counts can provide a quantitative description of infection and be used as a proxy for intensity of parasite infection in non-invasive studies of endangered non-human primate populations from which individuals cannot be removed, but the results should be interpreted cautiously owing to variable egg excretion. Negative samples may possibly include samples with very low parasite intensities as given by limitations of coprological methods. With exception of Mammomonogamus found in our samples, the exact determination of strongylids to the genus or species level on the basis of egg morphology is unreliable. However, previous morphological and molecular analyses of L3 larvae developed in coprocultures (Kalousová and Hasegawa, unpublished data), Hasegawa et al. (2014) showed that the thin-walled strongylid eggs belonged to either the genus Oesophagostomum or Necator. Therefore we divided the strongylid eggs found in gorilla samples into two categories: (i) Necator/Oesophagostomum and (ii) Mammomonogamus.
Microbial Community Analyses
Genomic DNA from fecal samples was extracted using the Power Soil DNA Isolation Kit (MoBio Laboratories Inc., Carlsbad, CA, USA) according to the manufacturer's protocol. The V1-V3 region of the 16S rRNA gene was PCR amplified (35 cycles: at 95°C for 30 s, at 55°C for 30 s and at 72°C for 30 s) using primers 27f (5′-AGAGTTTGATYMTGGCTCAG-3′, corresponding to nucleotides 27-47 of the Escherichia coli 16S rRNA gene) and 534r (5′-ATTACCGCGGCTGCTGGCA-3′, corresponding to nucleotides 534-515 of the E. coli 16S rRNA gene). Amplicon purification was done using the QIAquick PCR Purification Kit (QIAGEN, Germany). Amplicons were multiplexed and pyrosequenced using 454 FLX-Titanium technology as described previously (Gomez et al., 2015). Briefly, before following the mothur 454-SOP (Schloss et al., 2011), reads were quality filtered with an in-house perl script according to the following criteria: minimum and maximum length of 200 and 535 nt respectively; maximum number of homopolymers of 6; keeping all reverse and forward reads; removing all primer sequences and reads below Q = 35 over 50 bp windows. Unique sequences were aligned against the silva reference alignment database and chimeras were detected using uchime (Edgar et al., 2011) and removed. Sequences were then classified taxonomically using a Bayesian classifier approach implemented by mothur and reference files from the Ribosomal Database Project (RDP) (Wang et al., 2007) with a minimum cut-off of 80%. Then sequence reads with hits corresponding to unknown, mitochondria, chloroplasts, eukaryotes and archaea were eliminated. The remaining reads were clustered denovo using ModalClust. Reads sharing ≥ 97% 16S rRNA sequence complete-linkage similarity with the most abundant sequence were binned into an operational taxonomic unit (OTU). Taxonomic profiles were determined using the RDP Classifier (Wang et al., 2007) and the phylotype function within mothur. The QC script described above is designed in a way all spurious sequences are discarded and hence the possibility of overestimating diversity is reduced. Besides the astringent processing of the mothur 454 SOP (which emphasizes removal of spurious sequence reads), we removed OTUs present in fewer than 3 individuals and detected fewer than 5 times across the entire dataset to avoid possible 454 sequencing artifacts.
All statistical analyses were performed using the R statistical computing language v. 3.1.2 (R Development Core Team, 2008). UniFrac distance matrices, which take into account patterns of GIM variation based on phylogenetic relationships between community taxa, were calculated in QIIME (Caporaso et al., 2010). To evaluate differences in fecal microbiome of individual samples (individuals) positive and negative for particular parasite, relative abundance of each OTU based on Bray-Curtis dissimilarity matrices were used to perform non-metric multidimensional scaling (NMDS), permutational analysis of variance (PERMANOVA), and analysis of similarities (ANOSIM) in the vegan package (Oksanen et al., 2013). Vegan package was also used to calculate Shannon diversity indices, rarefied richness analysis, and analysis of fecal microbiome dispersion in multivariate space to compare the fecal microbiome diversity between individual samples positive and negative for particular parasite. Shannon diversity indices were used to evaluate the relationships between fecal microbiome diversity and egg counts of Necator/Oesophagostomum and Mammomonogamus. A heatmap showing the differences in abundance of the most abundant bacterial genera was drawn using the gplots package (Warnes et al., 2016). Bacterial indicator taxa in Entamoeba positive and Entamoeba negative individuals were identified using indicator species analysis (Dufrêne and Legendre, 1997) in the labdsv R package (Roberts, 2013) and confirmed using Random forest classifier in the package randomForest (Liaw and Wiener, 2002). These analyses were performed at the microbial genus, family, and phylum levels for all collected samples. To detect differences in relative abundances of indicator bacterial taxa according to Entamoeba status, the Generalized linear model (GLM, negative binomial distribution, log link) with Entamoeba status as a predictor variable and gorilla social group as a covariate was fitted using MASS package (Venables and Ripley, 2002) to evaluate the differences in relative abundances of indicator taxa between Entamoeba positive and Entamoeba negative individuals. Boxplots, ANOVAs, and linear regression models with false discovery rate corrections were calculated using the stats package (R Development Core Team, 2008). Spearman's rank correlations with false discovery rate corrections were used to describe relationships between strongylid infection (egg counts) and levels of bacterial taxa using psych package (Revelle, 2017). The relationships between GIPs and fecal microbiome were determined based on parasite taxa presence/absence data and relative abundances respectively. All analyses of variance, including PERMANOVA, were performed on all individuals (=samples) pooled together with gorilla social group as a control variable and using presence/absence data for each detected parasite as independent variables.
Amplicon sequences have been deposited in the MG-RAST database under project ID 6321. The datasets generated or analyzed during this study are included in Supplementary Information files or are available from the corresponding author on reasonable request.
We detected: trophozoites of Troglodytella/Gorillophilus, Prototapirella gorillae, and other unidentified entodiniomorphid ciliates; cysts of Blastocystis, Entamoeba (including cysts with eight and four nuclei); eggs of Strongyloides, Necator/Oesophagostomum, Mammomonogamus, Bertiella, and Spirurida fam. gen (Table 1, Table S2). The median of the EPG for Necator/Oesophagostomum and Mammomonogamus was 22.96 EPG (Min: 0; Max: 22.96) and 3.13 EPG (Min: 0; Max: 93.75), respectively (see Table S3).
Parasitic Infections and Gastrointestinal Microbiome
We obtained a median sampling depth of 6,458 reads per gorilla fecal sample (Min: 2,114; Max: 116,361) after sequence quality control (for further detail see Gomez et al., 2015, 2016b). The median of the raw reads was 6,559 (Min: 2,204; Max: 117,457) per sample. The total number of reads obtained in total was 421,553 and 416,624 after quality filtering. We found weak differences in the overall GIM composition only between Entamoeba positive (Ent+) and Entamoeba negative (Ent-) individuals (PERMANOVA for OTU relative abundances, based on Bray-Curtis dissimilarity matrices: Pseudo-F = 1.985, R2 = 0.043, p = 0.009, Table 2, Figure S2). Both weighted and unweighted UniFrac distance matrices also indicated differences in GIM composition between Ent+ and Ent- individuals (PERMANOVA: Pseudo-F = 1.668, R2 = 0.037, p = 0.053; Pseudo-F = 1.371, R2 = 0.031, p = 0.006, for weighted and unweighted UniFrac distances respectively; Table 2, Figure 1). We found no differences in overall GIM composition between Ent+ and Ent- gorillas according to ANOSIM analyses: Bray-Curtis dissimilarity matrices (R = 0.091, p = 0.017), weighted UniFrac (R = 0.040, p = 0.107), and unweighted UniFrac (R = 0.010, p = 0.296).
Table 2. Differences in fecal microbiome at the OTU level between individuals positive and negative for particular parasite based on PERMANOVA with gorilla social group as controlled variable.
Figure 1. NMDS plot of unweighted UniFrac distances of fecal microbiome composition of individual gorillas colored by Entamoeba status.
We detected significant differences in relative abundances of indicator taxa, detected by Species indicator analysis and Random forest classifier, in Ent- and Ent+ individuals based on GLM (negative binomial distribution, log link) with Entamoeba status as a predictor variable and gorilla social group as a covariate (Table 3, Figure 2). Specifically, at the phylum level, unknown Tenericutes (GLM: F = 2.992, p = 0.003) reached higher levels in Ent+ individuals. At the family level, Peptostreptococcaceae reached significantly higher levels in Ent− individuals (GLM: F = −2.559, p = 0.010), while unknown Selenomonadales (GLM: F = 3.935, p < 0.001), Erysipelotrichaceae (GLM: F = 3.505, p < 0.001), unknown Mollicutes (GLM: F = 3.282, p = 0.001), and Anaeroplasmataceae (GLM: F = 2.473, p = 0.013) reached higher levels in Ent+ individuals. At the genus level, unknown Selenomonadales (GLM: F = 3.697, p < 0.001), unknown Erysipelotrichaceae (GLM: F = 2.815, p = 0.005), unknown Molicutes (GLM: F = 3.260, p = 0.001), unknown Deltaproteobacteria (GLM: F = 4.217, p < 0.001), Olsenella (GLM: F = 2.852, p = 0.004), and Dorea (GLM: F = 3.167, p = 0.002) reached significantly higher levels in Ent+ individuals. Three out of six indicator taxa on genus level belonged among the 26 most abundant bacterial genera (Figure S3). Although Entamoeba status significantly predicted the abundance of these taxa, social group membership was also a significant factor in some of these taxa (GLM: p < 0.05; data not shown). We did not find significant differences in relative abundances of any other bacterial taxa between individuals positive or negative for other detected parasites (data not shown).
Figure 2. Relative abundances of indicator taxa based on Entamoeba status. Significant differences in relative abundances (%) of indicator taxa in Ent- (A) and Ent+ (B) individuals based on GLM with negative binomial distribution: ***0.001, **0.01, *0.5. Taxonomic rank: P, phyla; F, family; G, genus.
GIM diversity (Shannon diversity index) did not significantly differ between Ent+ and Ent- individuals, even though they reached higher values in Ent+ individuals (ANOVA: F = 2.823, p = 0.101; Figure 3). Moreover, we found no significant differences between Ent+ and Ent- individuals based on rarefied richness analysis (ANOVA: F = 1.129, p = 0.295). Analyses of multivariate dispersion of communities or inter-individual variation within groups also did not show differences between Ent+ and Ent- individuals (ANOVA: F = 0.047, p = 0.829; Figure S4).
Using linear regression models we found no relationships between egg counts of Necator/Oesophagostomum and Mammomonogamus (EPG) and overall fecal microbiome composition, as described by the first NMDS axis, NMDS1 (R2 = 0.009, t = 0.618, p = 0.540; R2 < 0.001, t = 0.025, p = 0.974, for Necator/Oesophagostomum and Mammomonogamus respectively; Figures 4A,B). The results were also confirmed by Spearman's rank correlations (Necator/Oesophagostomum, ρ = −2.795, p = 0.069; Mammomonogamus, ρ = −0.180, p = 0.247). We found weak relationships only between relative abundances of members of phylum Verrucomicrobia, family Verrucomicrobia subdivision 5 and Necator/Oesophagostomum egg counts based on Spearman's rank correlations and linear regression models (Figures S5A,B), but when we applied false discovery rate corrections, we did not find the significant results anymore, even if we defined the fdr-q value threshold as 0.1 (data not shown). The relationships between relative abundances of other detected bacterial taxa and Necator/Oesophagostomum egg counts were not statistically significant (data not shown). We did not find any significant relationship between relative abundances of any bacterial taxa and Mammomonogamus egg counts using Spearmans rank correlation analyses and linear regression model (data not shown).
Figure 4. Relationships between the egg counts of (A) Necator/Oesophagostomum spp., (B) Mammomonogamus sp. (natural log-transformed data) and overall GIM community (as described with the first NMDS axis, NMDS1) in individual gorillas.
We found no relationship between fecal microbiome diversity (Shannon diversity indices) and egg count of Necator/Oesophagostomum (Spearman's rank correlation: ρ = −0.218, p = 0.161; linear regression model: R2 = 0.005, t = −1.434, p = 0.159, Figure S6) or Mammomonogamus (Spearman's rank correlation: ρ = −0.010, p = 0.949; linear regression model: R2 < 0.001, t = 0.095, p = 0.093, Figure S6).
Given the evolutionary relationships between non-human primates and humans, studying associations between parasites and gut bacteria in great apes may help to understand current configurations of the human gut microbiome in the context of health and disease. Although the effect of GIP infection on the bacterial community composition of western lowland gorillas appears to be weak, our data show evidence that presence of Entamoeba spp. is associated with the abundance of specific bacterial taxa in the gorilla gastrointestinal tract. The mechanism behind these observations is unknown, however, this may be a reflection of either direct or indirect interaction between entamoebas, the GIM, and the host immune system. The genus Entamoeba contains both pathogenic and commensal species (Verweij et al., 2003; Stensvold et al., 2010). Although much is known about parasitic E. histolytica, the causative agent of amoebic dysentery in humans, and its impact on GIM (Rani et al., 2006; Galván-Moroyoqui et al., 2008; Verma et al., 2012; Burgess and Petri, 2016), there is no information about the non-pathogenic Entamoeba species and their relationship with host's GIM. In severe cases, E. histolytica trophozoites colonize the large intestine, especially the cecal and sigmoidorectal regions, where they feed on cellular debris and certain species of bacteria, allowing others to proliferate or inducing a host immune response that differentially affects the survival of different microbes (Morton et al., 2015; Burgess and Petri, 2016). It has been shown that direct competition by intestinal commensals can provide protection from invading pathogens, thus a disturbance of the normal GIM may result in increased susceptibility to pathogens (Hayes et al., 2010; Buffie and Pamer, 2013). Moreover, it is assumed that a specific gut microbiota predisposes an individual to Enatmoeba colonization (Bracha et al., 1982; Galván-Moroyoqui et al., 2008).
Our data indicate that the presence of Entamoeba is not associated with heterogeneity and dispersion of the fecal microbiome in lowland gorillas. Thus, either larger sample sizes may be needed to detect such differences or the effects of Entamoeba infections in the fecal microbiome differ between humans and non-human primates since a higher alpha diversity has been previously detected upon infection with Entamoeba spp. in hunter-gatherer and agricultural human populations (Morton et al., 2015). Nevertheless, we observed that in most cases, with the exception of family Peptostreptococcaceae, Ent+ individuals have a higher abundance of taxa that may play an important role in nutrition, and metabolism for the host (e.g., fiber degraders, carbohydrate and glucose fermenters) and taxa that seem to be a signature of the fecal microbiome of gorillas even if they belong among minor taxa (e.g., Olsenella, Dorea, and Verrucomicrobiacea) (Gomez et al., 2015, 2016b). We detected significant changes in members of dominant mammalian GIM phyla: Firmicutes, and Actinobacteria (Lay et al., 2005) associated with Entamoeba infection. We also found differences in members of Bacteroidetes, specifically unknown Bacteroidetes, nevertheless those were not significant. However, lower levels of members of the genus Bacteroides, which are important contributors to carbohydrate metabolism, nutrition and health of human and animals, have been previously observed in Ent+ people (Rani et al., 2006; Verma et al., 2012; Morton et al., 2015). Several genera belonging to Firmicutes reached significantly higher levels in Ent+ gorillas. For instance, Dorea, a common member of the human GIM, and important short-chain fatty acid producing bacteria (Guinane and Cotter, 2013; Schnorr et al., 2014), was more abundant in Ent+ individuals. A higher abundance of the Anaeroplasmataceae family was found in Ent+ individuals. These bacteria are known contributors to fermentative mechanisms in environmental micro-ecosystems (Koike et al., 2003); thus, its abundance may be beneficial for the host in the context of foraging efficiency. Members of family Erysipelotrichaceae, previously associated with several diseases in humans (Kaakoush, 2015), also reached higher levels in Ent+ individuals, however studies examining the direct impact that changes in the abundance of Erysipelotrichaceae have on the host are still required. Olsenella, a common member of GIM of great apes (Gomez et al., 2015; Moeller et al., 2016) reached higher relative abundances in Ent+ individuals. This genus was found to be more abundant in humans suffering from soil-transmitted helminths (Rosa et al., 2018). Olsenella has been previously associated with a reduction in gut inflammation (Wang et al., 2015), suggesting that its association with Entamoeba infection may potentially have positive side effects on the gorilla gut health.
Our observations raise questions about the natural occurrence of Entamoeba infections in the gorilla gastrointestinal ecosystem or the potential advantage of harboring these amoebas in terms of nutrition. Future studies should explore specific molecular interactions between hosts and their resident parasites and bacteria, seeking to determine if parasites help to retain nutritionally useful bacteria and displace those that are not essential or are harmful. Moreover, the prevalence and apparently harmless nature of some parasites (e.g., most of Entamoeba spp.) in gorillas must be tested molecularly, by specifically targeting the presence and expression of known parasite-associated virulence factors. Since large diversity of Entamoeba spp. has been observed in wild great apes (Vlčková et al., 2018), the molecular identification of particular Entamoeba species/genotypes may help to discriminate among particular parasite-bacteria interactions involving selected species, and the way they impact the GIM. Moreover, since many of the detected bacterial taxa, which significantly differed in studied individuals based on Entamoeba status belong to unknown species, it is difficult to evaluate the effect of these differences on the host. Due to the high complexity of the system, the metabolomic analyses should be included in the further studies to evaluate possible changes in the overall functionality of the GIM.
We did not find any associations between bacterial relative abundances as well as fecal microbiome diversity and intensity of Necator/Oesophagostomum and Mammomonogamus infection. However, our results could be affected by the inability to distinguish eggs of strongylid nematodes using the coprological methods as the life cycles and pathogeneses of these strongylids are different. Recently, a hidden diversity in the strongylid nematode communities in this gorilla population has been revealed by implementation ITS-2 metabarcoding at the Illumina MiSeq platform (Pafčo et al., 2018). Such approach allowed the detection of rare, otherwise overlooked taxa such as Ternidens diminutus and Libyostrongylus sp. present in a small number of samples and represented by a low numbers of reads. In future research, advanced molecular tools need to be implemented also to study GIP diversity to assess the associations between strongylids and GIM. Ensuing inflammatory responses from the host to strongylid infections may induce a gut environment with microbe-selective properties, which is further modulated by diet. As we did not record any clinical signs of strongylid infections in the studied gorillas, we may assume that the levels of Necator/Oesophagostomum infection may be too low and their abundance in the gorilla GIM may be facilitated by immune adaptations to pathogen exposure. It is also possible that strongylids may induce modifications of the GIM at the site of infection (e.g., duodenal mucosa for hookworms) that are not reflected in the fecal microbiome composition (Cooper et al., 2013). Along these lines, and in light of the potential therapeutic application of hookworms and suggestions of using chronic hookworm infection as immunomodulatory therapy to counteract detrimental immune responses (Daveson et al., 2011), the contention of parasites as a ubiquitous fraction of the GIM, interacting in synergy with specific GIM signatures in non-human primates and traditional populations should be explored.
Here, we examined the associations between GIP infections and the GIM composition of wild lowland gorillas. Although we did not observe notable changes in global bacterial community profiles in connection with GIP infections, we found significant relationships between Entamoeba infections and abundances of certain gut bacterial taxa. These changes could be mediated by complex interactions of the GIPs with the host ecology (diet), the host immune machinery, and the other microbiota. To gain further insight into the specific nature of these interactions, it will be necessary to apply additional molecular approaches to assess GIP diversity and function at the mucosal level to understand how the host immune system reacts to the GIPs and their virulence factors, with subsequent consequences for gut bacterial communities. Likewise, the role of diet in this interaction should be clarified from a molecular perspective in terms of metabolic impact on the GIPs and the other microbiota. For more complex analyses adding other factors affecting fecal microbiome of great apes, e.g., sex and age would be also beneficial, nevertheless, it would require studying captive or larger number of habituated animals. Moreover, future work may be directed into analyzing the evolutionary correlates of GIP infections for health and disease in humans and great apes. Our results can have implications for understanding host-microbe interactions and gorilla welfare in both wild and captive conditions. Absence of GIP in captive apes may be responsible for particular health or digestive problems (together with other factors like diet), on the other hand, it is well-known that particular GIP in captive apes can reach higher intensities than in the wild and infected animals suffer from clinical outcomes which are rarely observed in the wild.
KV and AG conceived the initial project design, with inputs from KP; BP KP and AT designed the field study and collected the samples; KV, AG, MT performed the microbial community analyses; BP performed parasitological analyses; KV analyzed the data with significant contributions from AG; KP, DM, KN, BreW, BryW, CY, RS, SL provided the financial support; KV wrote the initial manuscript with significant contributions from AG and KP. All authors reviewed the manuscript.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The reviewer AK and handling Editor declared their shared affiliation.
We thank the Government of the Central African Republic and to the World Wildlife Fund for granting permission to conduct our research in the Central African Republic; the Primate Habituation Programme (Dzanga-Ndoki National Park, Dzanga-Sangha Protected Areas) for logistical support in the field; K. A. Shutt and all local trackers and assistants for their help with sample collection. We thank Ran Blekhman for hosting KV for the data analyses in his laboratory at the Department of Genetics, Cell Biology, and Development, Department of Ecology, Evolution, and Behavior, University of Minnesota. This publication derives from the HPI-lab, Laboratory for Infectious Diseases Common to Humans and (non-Human) Primates, Czech Republic. This work was financially supported by the Internal Grant Agency of University of Veterinary and Pharmaceutical Sciences Brno 138/2015/FVL, by the institutional support of the Institute of Vertebrate Biology, Academy of Sciences of the Czech Republic (RVO: 68081766), by the project CEITEC- Central European Institute of Technology (CZ.1.05/1.100/02.0068) from European Regional Development Fund, and co-financed from European Social Fund and State Budget of the Czech Republic (project OPVK CZ.1.07/2.3.00/20.0300). Sequencing was supported by the US National Science Foundation (NSF) 0935347 and the University of Illinois at Urbana-Champaign.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2018.01202/full#supplementary-material
Anderson, R. M., and Schad, G. A. (1985). Hookworm burdens and faecal egg counts: an analysis of the biological basis of variation. Trans. R. Soc. Trop. Med. Hyg. 79, 812–825. doi: 10.1016/0035-9203(85)90128-2
Audebert, C., Even, G., Cian, A., Loywick, A., Merlin, S., et al. (2016). Colonization with the enteric protozoa Blastocystis is associated with increased diversity of human gut bacterial microbiota. Sci. Rep. 6:25255. doi: 10.1038/srep25255
Berrilli, F., Di Cave, D., Cavallero, S., and D'Amelio, S. (2012). Interactions between parasites and microbial communities in the human gut. Front. Cell. Infect. Microbiol. 2:141. doi: 10.3389/fcimb.2012.00141
Burgunder, J., Pafčo, B., Petrželková, K. J., Modrý, D., Hashimoto, C., and MacIntosh, A. J. J. (2017). Complexity in behavioural organization and strongylid infection among wild chimpanzees. Anim. Behav. 129, 257–268. doi: 10.1016/j.anbehav.2017.06.002
Cabaret, J., Gasnier, N., and Jacquiet, P. (1998). Faecal egg counts are representative of digestive-tract strongyle worm burdens in sheep and goats. Parasite 5, 137–142. doi: 10.1051/parasite/1998052137
Cantacessi, C., Giacomin, P., Croese, J., Zakrzewski, M., Sotillo, J., McCann, L., et al. (2014). Impact of experimental hookworm infection on the human gut microbiota. J. Infect. Dis. 210, 1431–1434. doi: 10.1093/infdis/jiu256
Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. doi: 10.1038/nmeth.f.303
Červená, B., Vallo, P., Pafčo, B., Jirku, K., Jirku, M., Petrželková, K. J., et al. (2017). Host specificity and basic ecology of Mammomonogamus (Nematoda, Syngamidae) from lowland gorillas and forest elephants in Central African Republic. Parasitology 144, 1016–1025. doi: 10.1017/S0031182017000221
Chapman, C. A., Speirs, M. L., Gillespie, T. R., Holland, T., and Austad, K. M. (2006). Life on the edge: gastrointestinal parasites from the forest edge and interior primate groups. Am. J. Primatol. 68, 397–409. doi: 10.1002/ajp.20233
Cooper, P., Walker, A. W., Reyes, J., Chico, M., Salter, S. J., Vaca, M., et al. (2013). Patent human infections with the whipworm, Trichuris trichiura, are not associated with alterations in the faecal microbiota. PLoS ONE 8:e76573. doi: 10.1371/journal.pone.0076573
Daveson, A. J., Jones, D. M., Gaze, S., McSorley, H., Clouston, A., Pascoe, A., et al. (2011). Effect of hookworm infection on wheat challenge in celiac disease - a randomised double-blinded placebo controlled trial. PLoS ONE 6:e17366. doi: 10.1371/journal.pone.0017366
Eckburg, P. B., Bik, E. M., Bernstein, C. N., Purdom, E., Dethlefsen, L., Sargent, M., et al. (2005). Diversity of the human intestinal microbial flora. Science 308, 1635 –1638. doi: 10.1126/science.1110591
Edgar, R. C., Haas, B. J., Clemente, J. C., Quince, C., and Knight, R. (2011). UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 27, 2194–2200. doi: 10.1093/bioinformatics/btr381
Gálván-Moroyoqui, J. M., Del Carmen Domínguez-Robles, M., Franco, E., and Meza, I. (2008). The Interplay between Entamoeba and enteropathogenic bacteria modulates epithelial cell damage. PLoS Negl. Trop. Dis. 2:e266. doi: 10.1371/journal.pntd.0000266
Gomez, A., Petrzelkova, K., Yeoman, C. J., Vlckova, K., Mrázek, J., Koppova, I., et al. (2015). Gut microbiome composition and metabolomic profiles of wild western lowland gorillas (Gorilla gorilla gorilla) reflect host ecology. Mol. Ecol. 24, 2551–2565. doi: 10.1111/mec.13181
Gomez, A., Petrzelkova, K. J., Burns, M. B., Yeoman, C. J., Amato, K. R., Vlckova, K., et al. (2016a). Gut microbiome of coexisting BaAka pygmies and Bantu reflects gradients of traditional subsistence patterns. Cell Rep. 14, 2142–2153. doi: 10.1016/j.celrep.2016.02.013
Gomez, A., Rothman, J. M., Petrzelkova, K., Yeoman, C. J., Vlckova, K., Umaña, J. D., et al. (2016b). Temporal variation selects for diet-microbe co-metabolic traits in the gut of Gorilla spp. ISME J. 10, 514–526. doi: 10.1038/ismej.2015.146
González-Hernández, M., Rangel-Negrín, A., Schoof, V. A. M., Chapman, C. A., Canales-Espinosa, D., and Dias, P. A. D. (2014). Transmission patterns of pinworms in two sympatric congeneric primate species. Int. J. Primatol. 35, 445–462. doi: 10.1007/s10764-014-9751-y
Guinane, C. M., and Cotter, P. D. (2013). Role of the gut microbiota in health and chronic gastrointestinal disease: understanding a hidden metabolic organ. Ther. Adv. Gastroenterol. 6, 295–308. doi: 10.1177/1756283X13482996
Hasegawa, H., Modrý, D., Kitagawa, M., Shutt, K. A., Todd, A., Kalousová, B., et al. (2014). Humans and great apes cohabiting the forest ecosystem in central African Republic Harbour the same hookworms. PLoS Negl. Trop. Dis. 8:e2715. doi: 10.1371/journal.pntd.0002715
Hayes, K. S., Bancroft, A. J., Goldrick, M., Portsmouth, C., Roberts, I. S., and Grencis, R. K. (2010). Exploitation of the intestinal microflora by the parasitic nematode Trichuris muris. Science 328, 1391–1394. doi: 10.1126/science.1187703
Hiemstra, H., Klaver, E., Vrijland, K., Kringel, H., Andreasen, A., Bouma, G., et al. (2014). Excreted/secreted Trichuris suis products reduce barrier function andsuppress inflammatory cytokine production of intestinalepithelial cells. Mol. Immunol. 60, 1–7. doi: 10.1016/j.molimm.2014.03.003
Hodges, M., Koroma, M. M., Baldé, M. S., Turay, H., Fofanah, I., Divall, M. J., et al. (2011). Current status of schistosomiasis and soil-transmitted helminthiasis in Beyla and Macenta Prefectures, Forest Guinea. Trans. R. Soc. Trop. Med. Hyg. 105, 672–674. doi: 10.1016/j.trstmh.2011.07.008
Kalousová, B. (2013). Gastrointestinal parasites of Wild Western Lowland Gorillas (Gorilla gorilla gorilla) in Dzanga Sector of Dzanga-Ndoki National Park, Central African Republic. Master thesis. Masaryk University. Czechia.
Kalousová, B. (2015). “Parasites of gorillas,” in Parasites of Great Apes: Atlas of Coproscopic Diagnostics, eds D. Modrý, K. Petrželková, B. Kalousová, and H. Hasegawa (Brno: HPI-lab, University of Veterinary and Pharm. Sciences Brno, Czech Republic), 103–106.
Kalousová, B., Hasegawa, H., Petrželková, K. J., Sakamaki, T., Kooriyma, T., and Modrý, D. (2016). Adult hookworms (Necator spp.) collected from researchers working with wild western lowland gorillas. Parasit. Vectors 9:75. doi: 10.1186/s13071-016-1357-0
Koike, S., Yoshitani, S., Kobayashi, Y., and Tanaka, K. (2003). Phylogenetic analysis of fiber-associated rumen bacterial community and PCR detection of uncultured bacteria. FEMS Microbiol. Lett. 229, 23–30. doi: 10.1016/S0378-1097(03)00760-2
Kreisinger, J., Bastien, G., Hauffe, H. C., Marchesi, J., and Perkins, S. E. (2015). Interactions between multiple helminths and the gut microbiota in wild rodents. Philos. Trans. R. Soc. B 370:20140295. doi: 10.1098/rstb.2014.0295
Kumar, S., Filipski, A., Swarna, V., Walker, A., and Hedges, S. B. (2005). Placing confidence limits on the molecular age of the human-chimpanzee divergence. Proc. Natl. Acad. Sci. U.S.A. 102, 18842–18847. doi: 10.1073/pnas.0509585102
Landsoud-Soukate, J., Tutin, C. E., and Fernandez, M. (1995). Intestinal parasites of sympatric gorillas and chimpanzees in the Lopé Reserve, Gabon. Ann. Trop. Med. Parasit. 89 73–79. doi: 10.1080/00034983.1995.11812931
Lay, C., Sutren, M., Rochet, V., Saunier, K., Doré, J., and Rigottier-Gois, L. (2005). Design and validation of 16S rRNA probes to enumerate members of the Clostridium leptum subgroup in human faecal microbiota. Environ. Microbiol. 7, 933–946. doi: 10.1111/j.1462-2920.2005.00763.x
Lilly, A. A., Mehlman, P. T., and Doran, D. (2002). Intestinal parasites in gorillas, chimpanzees, and humans at Mondika Research Site, Dzanga-Ndoki National Park, Central African Republic. Int. J. Primatol. 23, 555–573. doi: 10.1023/A:1014969617036
Masi, S., Chauffour, S., Bain, O., Todd, A., Guillot, J., and Krief, S. (2012). Seasonal effects on great ape health: a case study of wild chimpanzees and western gorillas. PLoS ONE 7:e49805. doi: 10.1371/journal.pone.0049805
Mekonnen, Z., Levecke, B., Boulet, G., Bogers, J.-P., and Vercruysse, J. (2013). Efficacy of different albendazole and mebendazole regimens against heavy-intensity Trichuris trichiura infections in school children, Jimma Town, Ethiopia. Pathog. Glob. Health 107, 207–209. doi: 10.1179/2047773213Y.0000000092
Moeller, A. H., Li, Y., Mpoudi Ngole, E., Ahuka-Mundeke, S., Lonsdorf, E. V., Pusey, A. E., et al. (2014). Rapid changes in the gut microbiome during human evolution. Proc. Natl. Acad. Sci. U.S.A. 111, 16431–16435. doi: 10.1073/pnas.1419136111
Moeller, A. H., Peeters, M., Ayouba, A., Ngole, E. M., Esteban, A., Hahn, B. H., et al. (2015). Stability of the gorilla microbiome despite simian immunodeficiency virus infection. Mol. Ecol. 24, 690–697. doi: 10.1111/mec.13057
Mohammadi, R., Hosseini-Safa, A., Ehsani Ardakani, M. J., and Rostami-Nejad, M. (2015). The relationship between intestinal parasites and some immune-mediated intestinal conditions. Gastroenterol. Hepatol. Bed Bench 8, 123–131.
Moreels, T. G., and Pelckmans, P. A. (2005). Gastrointestinal parasites. Potential therapy for refractory inflammatory bowel diseases. Inflamm. Bowel Dis. 11, 178–184. doi: 10.1097/00054725-200502000-00012
Morton, E. R., Lynch, J., Froment, A., Lafosse, S., Heyer, E., Przeworski, M., et al. (2015). Variation in rural African gut microbiomes is strongly shaped by parasitism and diet. PLoS Genet. 11:e1005658. doi: 10.1371/journal.pgen.1005658
Müller-Graf, C., Collins, D., and Woolhouse, M. E. (1996). Intestinal parasite burden in five troops of olive baboons (Papio cynocephalus anubis) in Gombe Stream National Park, Tanzania. Parasitology 112, 489–497. doi: 10.1017/S0031182000076952
Ngobeni, R., Samie, A., Moonah, S., Watanabe, K., Petri, Jr, W. A., and Gilchrist, C. (2017). Entamoeba species in South Africa: correlations with the host microbiome, parasite burdens, and first description of Entamoeba bangladeshi outside of asia. J. Infect. Dis. 216, 1592–1600. doi: 10.1093/infdis/jix535
Nielsen, M., Pfisterb, K., and von Samson-Himmelstjerna, G. (2014). Selective therapy in equine parasite control-application and limitations. Vet. Parasitol. 202, 95–103. doi: 10.1016/j.vetpar.2014.03.020
Obregon-Tito, A. J., Tito, R. Y., Metcalf, J., Sankaranarayanan, K., Clemente, J. C., Ursell, L. K., et al. (2015). Subsistence strategies in traditional societies distinguish gut microbiomes. Nat. Commun. 6:6505. doi: 10.1038/ncomms7505
Ochman, H., Worobey, M., Kuo, C.-H., Ndjango, J.-B. N., Peeters, M., Hahn, B. H., et al. (2010). Evolutionary relationships of wild hominids recapitulated by gut microbial communities. PLoS Biol. 8:e1000546. doi: 10.1371/journal.pbio.1000546
Pafčo, B., Benavides, J. A., Pšenková-Profousová, I., Modrý, D., Červená, B., Shutt, K. A., et al. (2017). Do habituation, host traits and seasonality have an impact on protist and helminth infections of wild western lowland gorillas? Parasitol. Res. 116, 3401–3410. doi: 10.1007/s00436-017-5667-8
Pafčo, B., Čížkovaá, D., Kreisinger, J., Hasegawa, H., Vallo, P., Shutt, K., et al. (2018). Metabarcoding analysis of strongylid nematode diversity in two sympatric primate species. Sci. Rep. 8:5933. doi: 10.1038/s41598-018-24126-3
Piekarska, J., Płoneczka-Janeczko, K., Kantyka, M., Kuczaj, M., Gorczykowski, M., and Janeczko, K. (2013). Gastrointestinal nematodes in grazing dairy cattle from small and medium-sized farms in southern Poland. Vet. Parasitol. 198, 250–253. doi: 10.1016/j.vetpar.2013.07.039
Rani, R., Murthy, R. S., Bhattacharya, S., Ahuja, V., Rizvi, M. A., and Paul, J. (2006). Changes in bacterial profile during amebiasis: demonstration of anaerobic bacteria in ALA pus samples. Am. J. Trop. Med. Hyg. 75, 880–885.
Rosa, B. A., Supali, T., Gankpala, L., Djuardi, Y., Sartono, E., Zhou, Y., et al. (2018). Differential human gut microbiome assemblages during soil-transmitted helminth infections in Indonesia and Liberia. Microbiome 6:33. doi: 10.1186/s40168-018-0416-5
Sak, B., Petrzelkova, K. J., Kvetonova, D., Mynarova, A., Shutt, K. A., Pomajbikova, K., et al. (2013). Long-term monitoring of microsporidia, cryptosporidium and giardia infections in Western Lowland Gorillas (Gorilla gorilla gorilla) at different stages of habituation in Dzanga Sangha Protected Areas, Central African Republic. PLoS ONE 8:e71840. doi: 10.1371/journal.pone.0071840
Schloss, P. D., Gevers, D., and Westcott, S. L. (2011). Reducing the effects of PCR amplification and sequencing artifacts on 16S rRNA-based studies. PLoS ONE 6:e27310. doi: 10.1371/journal.pone.0027310
Seivwright, L. J., Redpath, S. M., Mougeot, F., Watt, L., and Hudson, P. J. (2004). Faecal egg counts provide a reliable measure of Trichostrongylus tenuis intensities in free-living red grouse Lagopus lagopus scoticus. J. Helminthol. 78, 69–76. doi: 10.1079/JOH2003220
Šlapeta, J., Dowd, S. E., Alanazi, A. D., Westman, M. E., and Brown, G. K. (2015). Differences in the faecal microbiome of non-diarrhoeic clinically healthy dogs and cats associated with Giardia duodenalis infection: impact of hookworms and coccidia. Int. J. Parasitol. 45, 585–594. doi: 10.1016/j.ijpara.2015.04.001
Stensvold, C. R., Lebbad, M., Verweij, J. J., Jespersgaard, C., von Samson-Himmelstjerna, G., Nielsen, S. S., et al. (2010). Identification and delineation of members of the Entamoeba complex by pyrosequencing. Mol. Cell. Probes 24, 403–406. doi: 10.1016/j.mcp.2010.07.008
Su, C., Cao, Y., Kaplan, J., Zhang, M., Li, W., Conroy, M., et al. (2011). Duodenal helminth infection alters barrier function of the colonic epithelium via adaptive immune activation. Infect. Immun. 79, 2285–2294. doi: 10.1128/IAI.01123-10
Summers, R. W., Elliott, D. E. Jr, J. F. U., Thompson, R. A., and Weinstock, J. V. (2005). Trichuris suis therapy for active ulcerative colitis: a randomized controlled trial. Gastroenterology 128, 825–832. doi: 10.1053/j.gastro.2005.01.005
Tandon, B. N., Tandon, R. K., Satpathy, B. K., and Shriniwas (1977). Mechanism of malabsorption in giardiasis: a study of bacterial flora and bile salt deconjugation in upper jejunum. Gut 18, 176–181. doi: 10.1136/gut.18.3.176
Tomkins, A. M., Wright, S. G., Drasar, B. S., and James, W. P. (1978). Bacterial colonization of jejunal mucosa in giardiasis. Trans. R. Soc. Trop. Med. Hyg. 72, 33–36. doi: 10.1016/0035-9203(78)90294-8
Verma, A. K., Verma, R., Ahuja, V., and Paul, J. (2012). Real-time analysis of gut flora in Entamoeba histolytica infected patients of Northern India. BMC Microbiol. 12:183. doi: 10.1186/1471-2180-12-183
Verweij, J. J., Laeijendecker, D., Brienen, E. A., van Lieshout, L., and Polderman, A. M. (2003). Detection and identification of Entamoeba species in stool samples by a reverse line hybridization assay. J. Clin. Microbiol. 41, 5041–5045. doi: 10.1128/JCM.41.11.5041-5045.2003
Vlčková, K., Gomez, A., Petrželková, K. J., Whittier, C. A., Todd, A. F., Yeoman, C. J., et al. (2016). Effect of antibiotic treatment on the gastrointestinal microbiome of free-ranging western lowland gorillas (Gorilla g. gorilla). Microb. Ecol. 72, 943–954. doi: 10.1007/s00248-016-0745-5
Vlčková, K., Kreisinger, J., Pafčo, B., Čížková, D., Tagg, N., Hehl, A. B., et al. (2018). Diversity of Entamoeba spp. in African great apes and humans: An insight from Illumina MiSeq high-throughput sequencing. Int. J. Parasitol. 48, 519–530. doi: 10.1016/j.ijpara.2017.11.008
Walk, S. T., Blum, A. M., Ewing, S. A.-S., Weinstock, J. V., and Young, V. B. (2010). Alteration of the murine gut microbiota during infection with the parasitic helminth, Heligmosomoides polygyrus. Inflamm. Bowel Dis. 16, 1841–1849. doi: 10.1002/ibd.21299
Wang, J., Tang, H., Zhang, C., Zhao, Y., Derrien, M., Rocher, E., et al. (2015). Modulation of gut microbiota during probiotic-mediated attenuation of metabolic syndrome in high fat diet-fed mice. ISME J. 9, 1–15. doi: 10.1038/ismej.2014.99
Wang, Q., Garrity, G. M., Tiedje, J. M., and Cole, J. R. (2007). Naïve Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl. Environ. Microbiol. 73, 5261–5267. doi: 10.1128/AEM.00062-07
Wu, S., Li, R. W., Li, W., Beshah, E., Dawson, H. D., and Urban, J. F. Jr. (2012). Worm burden-dependent disruption of the porcine colon microbiota by Trichuris suis infection. PLoS ONE 7:e35470. doi: 10.1371/journal.pone.0035470
Keywords: lowland gorilla, fecal microbiome, bacteria, parasite infection, Entamoeba, strongylid nematodes
Citation: Vlčková K, Pafčo B, Petrželková KJ, Modrý D, Todd A, Yeoman CJ, Torralba M, Wilson BA, Stumpf RM, White BA, Nelson KE, Leigh SR and Gomez A (2018) Relationships Between Gastrointestinal Parasite Infections and the Fecal Microbiome in Free-Ranging Western Lowland Gorillas. Front. Microbiol. 9:1202. doi: 10.3389/fmicb.2018.01202
Received: 08 December 2017; Accepted: 16 May 2018;
Published: 15 June 2018.
Edited by:Thomas Dandekar, Universität Würzburg, Germany
Reviewed by:Alexander Keller, Universität Würzburg, Germany
Vikram Saini, University of Alabama at Birmingham, United States
Copyright © 2018 Vlčková, Pafčo, Petrželková, Modrý, Todd, Yeoman, Torralba, Wilson, Stumpf, White, Nelson, Leigh and Gomez. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Klára Vlčková, firstname.lastname@example.org