Research Article | Host-Microbe Biology

Asymptomatic Intestinal Colonization with Protist Blastocystis Is Strongly Associated with Distinct Microbiome Ecological Patterns

M. E. Nieves-Ramírez, O. Partida-Rodríguez, I. Laforest-Lapointe, L. A. Reynolds, E. M. Brown, A. Valdez-Salazar, P. Morán-Silva, L. Rojas-Velázquez, E. Morien, L. W. Parfrey, M. Jin, J. Walter, J. Torres, M. C. Arrieta, C. Ximénez-García, B. B. Finlay
Catherine Lozupone, Editor
DOI: 10.1128/mSystems.00007-18

ABSTRACT

Blastocystis is the most prevalent protist of the human intestine, colonizing approximately 20% of the North American population and up to 100% in some nonindustrialized settings. Blastocystis is associated with gastrointestinal and systemic disease but can also be an asymptomatic colonizer in large populations. While recent findings in humans have shown bacterial microbiota changes associated with this protist, it is unknown whether these occur due to the presence of Blastocystis or as a result of inflammation. To explore this, we evaluated the fecal bacterial and eukaryotic microbiota in 156 asymptomatic adult subjects from a rural population in Xoxocotla, Mexico. Colonization with Blastocystis was strongly associated with an increase in bacterial alpha diversity and broad changes in beta diversity and with more discrete changes to the microbial eukaryome. More than 230 operational taxonomic units (OTUs), including those of dominant species Prevotella copri and Ruminococcus bromii, were differentially abundant in Blastocystis-colonized individuals. Large functional changes accompanied these observations, with differential abundances of 202 (out of 266) predicted metabolic pathways (PICRUSt), as well as lower fecal concentrations of acetate, butyrate, and propionate in colonized individuals. Fecal calprotectin was markedly decreased in association with Blastocystis colonization, suggesting that this ecological shift induces subclinical immune consequences to the asymptomatic host. This work is the first to show a direct association between the presence of Blastocystis and shifts in the gut bacterial and eukaryotic microbiome in the absence of gastrointestinal disease or inflammation. These results prompt further investigation of the role Blastocystis and other eukaryotes play within the human microbiome.

IMPORTANCE Given the results of our study and other reports of the effects of the most common human gut protist on the diversity and composition of the bacterial microbiome, Blastocystis and, possibly, other gut protists should be studied as ecosystem engineers that drive community diversity and composition.

INTRODUCTION

Blastocystis is the most prevalent protozoan resident of the human intestine, colonizing approximately 0.5 to 30% of humans in industrialized countries and 30 to 100% in nonindustrialized societies (1–7). Blastocystis is transmitted via the oral-fecal route, and colonization has been linked to poor hygiene, animal contact, and contamination of food and water sources (4, 8). And yet, Blastocystis is also common in industrialized populations, with 20 to 50% prevalence among healthy cohorts in Europe (9, 10), and interestingly, is associated with high income in Denmark (11). Similarly, its role in human health and disease is controversial. Blastocystis has been associated with a variety of gastrointestinal symptoms (12–14) and inflammation, both within the gut (15) and systemically (16, 17). Several studies have suggested that Blastocystis may play a role in the development of irritable bowel syndrome (IBS) (11, 18, 19). However, Blastocystis has often been described as an asymptomatic member of the normal intestinal microbiota (7, 9, 10, 20) and has even been inversely associated with body mass index and Crohn’s disease (7). It remains unclear how Blastocystis is capable of acting as an opportunistic pathogen or a commensal or a beneficial microbe and to what extent this depends on the specific circumstances of its host.

The plethora of human-associated microbiome studies documenting the role of the microbiome in health and disease has directed considerable attention to the factors that influence microbiome diversity and composition. To date, diet (21, 22), antibiotics (23), age (24), inflammation (25), and to a lesser extent, host genetics (26–28) have been shown to impact the ecology of the gut bacterial community, yet most studies are limited to profiling bacteria and seldom consider interdomain, multitrophic interactions as one of these factors. The relationship between human-associated gut protists and the resident gut bacterial community has only recently begun to be explored. Pathogenic protists, such as Giardia and Entamoeba histolytica, are associated with an increase in gut bacterial diversity or with compositional differences (29–32), likely due to the mucosal damage and inflammation they cause in the intestinal mucosa. Interestingly, nonpathogenic protists may also cause major shifts in the bacterial microbiota: Morton et al. document an association between the presence of Entamoeba and increased diversity and compositional shifts in the bacterial microbiota in asymptomatic individuals, a portion of whom likely harbor nonpathogenic Entamoeba dispar (33). Recent reports in humans also suggest that Blastocystis is associated with compositional and diversity differences of the bacterial gut microbiota (2, 7, 34). However, given that this protist is known to cause disease in certain individuals and that the reported results are confounded by disease status, it remains unknown whether the changes in bacterial microbiomes are due to the presence of Blastocystis or to ongoing inflammation. Differentiating whether Blastocystis causes direct ecological effects on the bacterial microbiome or if microbiota shifts are mainly mediated through immune and physiologic changes associated with Blastocystis-associated inflammation will help determine the ecological consequences of decreased Blastocystis incidence in the westernized gut human microbiome.

In light of this, we aimed to study the gut microbiome in relation to Blastocystis colonization in a large sample of healthy individuals residing in a semi-industrialized setting in rural Mexico. Through this investigation, we detected marked taxonomic and functional bacterial differences associated with asymptomatic colonization of Blastocystis that support its role as a common eukaryotic gut commensal that drastically influences the bacterial microbiota through currently unknown mechanisms.

RESULTS

Prevalence of Blastocystis.We studied 156 individuals from the town of Xoxocotla, in the state of Morelos, Mexico. This population was deemed free of gastrointestinal symptoms associated with disease, according to physical examination by a clinician and data obtained from the Rome III questionnaire (see Materials and Methods for criteria in this questionnaire). Besides the Rome III questionnaire, a complete medical history was carried out for each individual, emphasizing inflammatory symptoms and antibiotic use in the 6 months preceding the clinical assessment and sampling. None of the selected samples came from individuals that reported antibiotic use or previous inflammatory conditions. The gender distribution was 84 women (53.8%) and 72 men (46.2%), with a median age of 27 years (10 to 53 years). The majority of individuals in the study had some degree of education (84.6%). The most frequent level was secondary school (33.0%), followed by primary school (32%), preschool (26%), technical education (8.4%), and professional studies (0.49%) (see Table S1 in the supplemental material for baseline characteristics of study participants). The status of colonization by Blastocystis and other protozoa (Cryptosporidium parvum, Entamoeba histolytica/E. dispar, and Giardia intestinalis) was determined by microscopy and quantitative PCR (qPCR) assays targeting regions of the small subunit (SSU) rRNA gene to detect and quantify levels directly from human stool specimens. Blastocystis colonization was detected in 65% of individuals (102/156), 52% of whom (53/102) were women. There was no difference in the mean ages of Blastocystis-positive versus -negative individuals (27.8 ± 8.4 years [mean ± standard deviation] versus 28.1 ± 7.9 years, respectively; Student’s t test, P = 0.83). Targeted 18S rRNA gene sequencing results allowed further resolution of Blastocystis taxonomy, showing that Blastocystis subtype 3 (ST3) is the prevalent subtype colonizing this human population. We did not identify any individuals infected by Crysptosporidium parvum or Giardia intestinalis but did observe coinfection with Entamoeba histolytica/E. dispar in five individuals, who were excluded from the bioinformatics analysis in order to assess the effect of Blastocystis only.

TABLE S1 

Characteristics of study population. Download TABLE S1, DOCX file, 0.04 MB.
Copyright © 2018 Nieves-Ramírez et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Intestinal microbiota composition and Blastocystis infection. (i) Prokaryotes.We determined the bacterial community by amplification and sequencing of the 16S rRNA gene (V3 region). Blastocystis colonization was associated with profound changes in bacterial alpha and beta diversity (Fig. 1) (for analysis of alpha diversity, Chao1 and Shannon indices and Mann-Whitney test were used [P < 0.001], and for analysis of beta diversity, principal-component analysis [PCoA] and permutational multivariate analysis of variance [PERMANOVA] were used [P < 0.001]). The marked change in beta diversity can be explained by large differences in the abundances of predominant bacterial taxa, including Prevotella copri, Prevotella stercorea, Ruminoccoccus bromii, Alistipes putredinis, Bacteroides species, Bifidobacterium longum, and Oscillospira species (DESeq2, Wald test, and false discovery rate [FDR], P < 0.05) (Table 1; the full list of differential operational taxonomic units [OTUs] is in Table S2). The presence of Blastocystis was associated with significant differences in 7 of the 10 most abundant OTUs, including Prevotella copri and Ruminococcus bromii, the two most abundant taxa in the data set (Fig. 1D; Table 1). Similarly, at the genus level, 7 of the 10 most abundant genera, comprising ~65% of all bacterial sequences, were significantly different in relation to their abundances in the bacterial community associated with Blastocystis colonization (Fig. S1).

FIG S1 

Relative abundances of the 10 most abundant bacterial genera relative to the presence of Blastocystis. Significant differences are shown, calculated by DESeq2 (Wald test plus FDR). Download FIG S1, TIF file, 2.3 MB.
Copyright © 2018 Nieves-Ramírez et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

TABLE S2 

Differential abundances of microbial OTUs (FDR < 0.05) between colonized and noncolonized individuals, calculated by DESeq2. Download TABLE S2, XLSX file, 0.1 MB.
Copyright © 2018 Nieves-Ramírez et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

TABLE S3 

Differential PICRUSt-predicted metabolic pathways between Blastocystis-positive and -negative individuals. Download TABLE S3, XLSX file, 0.03 MB.
Copyright © 2018 Nieves-Ramírez et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

FIG 1 
FIG 1 

Variations in beta and alpha diversity of gut microbiome bacterial communities in relation to presence of Blastocystis. (A) Principal-component analysis (PCoA) ordination of variation in beta diversity of human gut bacterial communities based on Bray-Curtis dissimilarities among fecal samples. Colors represent the presence of Blastocystis in gut microbial communities (red for negative and blue for positive), and arrows represent the significant (P < 0.001) correlations between PCoA axes versus the relative abundances of bacterial genera in communities. (B and C) Shannon diversity (B) and Chao1 estimated richness (C) display differences in alpha diversity. Significant differences are shown by P values of Mann-Whitney tests for comparison between 2 groups. (D) Relative abundances of the 10 most abundant bacterial OTUs relative to the presence of Blastocystis. Significant differences are shown, calculated by DESeq2 (Wald test plus FDR).

View this table:
TABLE 1 

Differential abundances of taxa in relation to Blastocystis colonization calculated by DESeq2a

Given the metadata collected in this cohort (age, education level, and family relationships), we captured the effects of the Blastocystis colonization on the gut bacterial microbiome while controlling for potential confounding variables using MaAsLin (35). MaAsLin is a multivariate linear modeling tool with boosting that tests for associations between specific microbial taxa and continuous and/or Boolean metadata. This method reduces the total amount of correlations to be tested, therefore allowing for improvements in the speed and the robustness of the additive general linear models. With MaAsLin, we found a significant association between most of the differential taxa (previously obtained through the DESeq2 analysis) and Blastocystis status, whereas no other variable explained the taxonomic differences observed in colonized individuals (Table S4).

TABLE S4 

Bacterial taxa associated with Blastocystis colonization corrected for possible confounding effects using MaAsLin. Direction of the correlation coefficient is in relation to the data for the noncolonized group. Findings are FDR corrected. Download TABLE S4, DOCX file, 0.1 MB.
Copyright © 2018 Nieves-Ramírez et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

We also controlled for the effect of sequencing depth influencing our results, as previously described by Weiss et al. (36), and found that while sequencing depth did impact the ordination of our data, the influence of Blastocystis species described in our study is not driven by the variation in sequencing depth across samples (Fig. S3A) (PERMANOVA, P = 0.5).

(ii) Eukaryotes.We determined the eukaryotic microbiota by amplification and sequencing of the 18S rRNA gene (V4 region). In the 18S rRNA gene data set, 63 samples had fewer than 1,000 sequences per sample after applying filtering steps (see Materials and Methods), which prompted their removal from the data set. Thus, 93 samples (Npos = 60 and Nneg = 33, from colonized and noncolonized individuals, respectively) were retained for ecological analysis of this data set. There were no significant differences with respect to baseline characteristics between the reduced number of samples and the original data set used for 16S rRNA gene analysis (Table S5).

TABLE S5 

Baseline characteristics in reduced cohort used in 18S sequencing analysis. Download TABLE S5, DOCX file, 0.04 MB.
Copyright © 2018 Nieves-Ramírez et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Blastocystis colonization was associated with more discrete differences in the eukaryotic microbiota (Fig. 2). As expected, our analysis detected Blastocystis subtype 3 as the protist most significantly different in abundance between Blastocystis-positive and -negative individuals (Table 1; see full list in Table S2). Blastocystis colonization was also associated with increases in yeast and fungal species (Debaryomyces hansenii, Mucor mucedo, Aspergillus flavus, Mucor racemosus, and Issatchenkia terricola) and a decrease of Hymenolepis nana (Table 1; Table S2). Since H. nana is a known cestode colonizer in the human gut, often with pathological consequences, we evaluated whether the presence of H. nana influenced bacterial community structure and confounded our results. However, we found that H. nana did not drive a significant effect in bacterial community structure (Fig. S2A). The detection of a significant difference in H. nana abundances between Blastocystis-positive and -negative individuals is likely due to a high variation in abundance in both groups and relatively few samples where they are represented. Likewise, we investigated the impact of the presence of Debaryomyces hansenii in the microbiota because this yeast species has been associated with human infections (37). Interestingly, the presence of this yeast explained 3.3% of bacterial community variation (PERMANOVA, P > 0.001). However, the distribution of D. hansenni did not overlap the distribution of Blastocystis (Fig. S2B), and their effects on the bacterial community are independent from each other (PERMANOVA, P = 0.68).

FIG S2 

Influence of Hymenolepis nana and Debaryomyces hansenii on bacterial community structure. Principal-component analysis (PCoA) ordination of variations in beta diversity of human gut bacterial communities based on Bray-Curtis dissimilarities among samples. (A) Colors represent the presence of Blastocystis in gut microbial communities (red for negative and blue for positive). Circles highlighted in yellow represent samples positive for Hymenolepis nana. (B) Ordination of bacterial community structures showing samples positive (yellow) or negative (black) for fungus D. hansenii and also showing the presence (triangles) and absence (circles) of Blastocystis. Download FIG S2, TIF file, 2.5 MB.
Copyright © 2018 Nieves-Ramírez et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

FIG 2 
FIG 2 

Variations in beta and alpha diversity of gut microbiome eukaryotic communities explained by presence of Blastocystis. (A) Principal-component analysis (PCoA) ordination of variation in beta diversity of human gut eukaryotic communities based on Bray-Curtis dissimilarities among samples. Colors represent the presence of Blastocystis in gut microbial communities (red for negative and blue for positive). (B and C) Shannon diversity (B) and Chao1 estimated richness (C) display differences in alpha diversity. Significant differences are shown by P values of Mann-Whitney tests for comparison between 2 groups. (D) Relative abundances of the 10 most abundant eukaryotic OTUs relative to the presence of Blastocystis. Significant differences, calculated by DESeq2 (Wald test plus FDR), are shown.

While controlling for sequencing depth, we observed that it interacted significantly with the effect of Blastocystis species colonization on community structure (PERMANOVA, P = 0.02) (Fig. S3B). This result is biologically intuitive since the presence of Blastocystis species in itself is contributing to the number of sequences detected. However, because of this result, we carried out the DESeq2 analysis using variance-stabilizing-transformed data (38). This analysis yielded the same numbers and identities of differential OTUs between Blastocystis-colonized and -noncolonized individuals as we obtained with nontransformed data.

FIG S3 

Influence of sequencing depth on bacterial and eukaryotic community structures. (A and B) Ordination of the bacterial (A) and eukaryotic (B) community structures showing samples positive (triangles) and negative (circles) for Blastocystis and also showing sequencing depth using a color scale (increasing sequence numbers from blue to yellow). Download FIG S3, TIF file, 2.4 MB.
Copyright © 2018 Nieves-Ramírez et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Altogether, the eukaryomes’ compositional differences in relation to Blastocystis colonization resulted in only a marginal difference in eukaryome beta diversity (PCoA, PERMANOVA, P = 0.01) (Fig. 2A). Blastocystis colonization was associated with a smaller, yet significant increase in eukaryome alpha diversity (Shannon index, Mann-Whitney, P = 0.014) (Fig. 2B) compared to the changes in bacterial diversity. There were no changes in community richness (Chao1, Mann-Whitney, P = 0.49).

(iii) Interdomain associations.Correlative analysis of the abundances of the top 100 most abundant taxa of the 16S versus 18S rRNA gene data sets showed, as expected, that Blastocystis subtype 3 was positively correlated with members of the Ruminoccaceae family and inversely correlated with Prevotella copri (Fig. 3).

FIG 3 
FIG 3 

Heatmap of biweight correlations (bicor method) between top 100 bacterial (x axis) and top 100 taxon (y axis) OTUs in fecal samples from study participants. Colors denote positive (red) and negative (blue) correlation values. Significant correlations are denoted with a plus sign (P < 0.05; FDR).

Bacterial-community-derived functional changes associated with Blastocystis colonization.To determine the functional differences associated with Blastocystis colonization, we inferred metagenomics potential using PICRUSt (39). Metagenomes predicted from the 16S rRNA gene data revealed 202 differential biochemical pathways out of 266 biochemical pathways associated with Blastocystis colonization (Welch’s t test) (Table S3), strongly suggesting that the microbiomes in the two groups were functionally distinct. The functional differences involved diverse metabolic functions, including metabolism of secondary bile acids, lipids, tryptophan and other amino acids, and carbohydrates. To confirm the predicted changes in carbohydrate metabolism, we next directly measured the concentration of fecal short-chain fatty acids (SCFA) in all samples. As predicted by our PICRUSt analysis, colonization with Blastocystis was associated with vastly different bacterial colonic fermentation patterns, as measured by production of SCFA. The three main SCFA, acetate, butyrate, and propionate, were reduced in colonized individuals, whereas caproate was increased (Fig. 4).

FIG 4 
FIG 4 

Short-chain fatty acid production in relation to Blastocystis colonization. Concentrations of fecal acetate, propionate, butyrate, isobutyrate, isovalerate, valerate, and caproate were measured by gas chromatography (Npos = 102, Nneg = 54; *, P < 0.05; ***, P < 0.001 [Mann-Whitney]).

Immunological differences in Blastocystis-colonized individuals. (i) Fecal calprotectin.To determine whether Blastocystis colonization was not associated with intestinal inflammation (as inferred by the Rome III criteria questionnaires), we measured fecal calprotectin in a subset of samples (Npos = 26 and Nneg = 17). Calprotectin is a neutrophil cytosolic protein commonly used as a marker of gut inflammation (40). Strikingly, lower concentrations of fecal calprotectin were observed in Blastocystis*-colonized individuals than in noncolonized individuals (Mann-Whitney, P = 0.0003) (Fig. 5A). All calprotectin values were well under levels consistent with pathological gut inflammation (>200 µg/mg), in accordance with the lack of clinical symptoms noted in our study cohort.

FIG 5 
FIG 5 

Gut mucosal immune changes associated with Blastocystis colonization. Concentrations of calprotectin (A) and total IgA (B) were determined by ELISA in feces in a subset of samples from this study (Nneg = 17, Npos = 26; *, P < 0.05; ***, P < 0.001 [Mann-Whitney]).

(ii) Fecal IgA.Immunoglobulin A (IgA) is the most abundant antibody in mucosal surfaces, and it has a critical role in maintaining homeostasis with the microbiome (41). It acts by binding to and neutralizing invading pathogens, and it can specially target microbes close to the mucus layer (42). IgA secretion into the lumen is elicited by parasitic infections with helminths, and it is important in limiting parasite fecundity and providing immune protection against reinfection (43).

We measured the total IgA concentration in feces to determine if asymptomatic Blastocystis colonization was associated with changes in this important humoral immune factor. In a subset of individuals (Npos = 26 and Nneg = 17), we observed that, similar to the calprotectin finding, colonized individuals exhibited lower levels of fecal IgA (Mann-Whitney, P = 0.03) (Fig. 5B).

(iii) Oxidative stress markers in urine.It has been described that oxidative stress and the presence of free-radical activities in hosts colonized by parasites correlate with the production of specific metabolites, such as advanced oxidative protein product (AOPP) and malondialdehyde (MDA) (44–46). Previous studies have evaluated MDA levels among healthy individuals and people infected by Ascaris lumbricoides (45), Entamoeba histolytica, and Plasmodium vivax (44), and all revealed increased levels of MDA in the infected individuals, suggesting that oxidative stress acts as a mediator of tissue damage concurrent with different parasites. We measured these biomarkers in urine samples of most study participants (Npos = 102 and Nneg = 53), as described previously (13, 47), and found no significant differences between colonized and noncolonized individuals (Table S6), suggesting that Blastocystis colonization is not accompanied by oxidative stress in this population.

TABLE S6 

Oxidative stress biomarkers measured in urine. Download TABLE S6, DOCX file, 0.1 MB.
Copyright © 2018 Nieves-Ramírez et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

(iv) Serum cytokines.Blastocystis ST7 has been shown to induce the expression of proinflammatory cytokines in mouse intestinal explants (15) and in serum of patients with IBS (16). However, in a reduced subset of our population, no differences in the serum concentrations of cytokines interleukin-2 (IL-2), IL-4, IL-6, IL-10, IL-17, tumor necrosis factor alpha (TNF-α), and gamma interferon (IFN-γ) were detected between Blastocystis-infected and noninfected individuals (Npos = 16 and  Nneg = 9, Mann-Whitney) (Table S7), further supporting the finding that Blastocystis colonization in our cohort is not associated with inflammation. This may reflect on the different immune effects exerted by different Blastocystis subtypes; previous work reported differential induction of IL-10 but not of IL-12 or IL-8 by subtype (17). However, Blastocystis pathogenicity does not directly correlate with subtype identity (4, 10), and pathogenicity is likely a better predictor of cytokine induction. These results reflect data from a very reduced number of blood samples available in this study and should be confirmed in a larger sample size.

TABLE S7 

Cytokine concentrations in serum samples for a subset of study participants, measured by cytometric bead array. Download TABLE S7, DOCX file, 0.1 MB.
Copyright © 2018 Nieves-Ramírez et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

DISCUSSION

The intestinal microbiota is highly variable among individual humans, and its diversity is affected by factors like diet, sociogeographic setting, antibiotic use, disease, age, and to a lesser degree, genetics (24, 27, 48). A direct association of intestinal parasites, such as helminths and protozoa, with human intestinal microbiota composition and diversity has been previously reported (33, 49). However, the influence of gastrointestinal inflammation had not been previously accounted for, and thus, it remained unanswered whether the microbiome differences associated with parasite colonization were linked to ecological parameters and/or were a direct result of gut inflammation. Our study focused on residents from a rural Mexican community that were confirmed to be healthy and without clinical symptoms. This population provided an ideal cohort to study how the most common human protist altered the ecology of the bacterial and eukaryotic microbiome in the absence of inflammation. Similar to previous studies (2, 34), our results confirm that colonization with Blastocystis is strongly associated with broad shifts in the gut-resident bacterial community and an increase in bacterial alpha diversity, providing very strong evidence that the presence of this protist may influence human gut microbial ecology through mechanisms that remain unknown.

In this population, Blastocystis colonization was very strongly correlated with a decrease in Prevotella copri and an increase in Ruminococcus bromii, the first and second most abundant bacterial phylotypes in this data set, respectively. These taxa are so commonly observed within the human gut microbiome that individuals are often described as either Prevotella rich or Ruminococcus rich (50). Prevotella species are dominant colonizers of the human gut in nonwesternized populations that consume a plant-rich diet (24, 51), whereas Ruminococcus bromii has been proposed to be a keystone species of the human gut, due to its colonization capacity and exceptional ability to degrade resistant starches (52, 53). In this context, this study reveals a potential role of Blastocystis as a discriminant feature between two fundamental bacterial species of the human gut, as Prevotella copri is positively correlated with Blastocystis subtype 3 (ST3), while Ruminococcus bromii is negatively correlated with Blastocystis ST3. Blastocystis colonization is thus a potentially important factor in structuring the gut microbiota architecture.

While other variables could explain these striking differences in microbial communities, age, family relationships, and health status were controlled for in this study. Although a full dietary assessment was not carried out in this population, dietary patterns in this rural population have historically been observed to be highly homogenous, characterized by a high intake of both plant-based sources (corn, legumes, chilies, tomatoes, and other fruits and vegetables) and meat (mainly chicken and beef). Thus, we do not expect our results to be explained by dietary differences among individuals in our study.

Instead, we postulate that Blastocystis may exert a predatory effect on bacteria of highly abundant taxa, such as Prevotella copri, which accounts for ~50% of the bacterial abundance in Blastocystis-negative individuals. Blastocystis is known to graze on bacteria. Dunn et al. (54) reported that the ameboid form of the protist was capable of bacterial engulfment, a process that has been suggested to serve the nutritional needs of encystation (55). Additional evidence of the grazing of Blastocystis on bacteria includes (i) the low frequency of the ameboid form in axenic cultures (56) and (ii) the contact between Blastocystis pseudopodia and bacteria (57). Bacterial predation is also known to occur in the case of a group of free-living amoebas, which are the only group of organisms that can cause a decrease of a bacterial population from 108 down to 105 per gram of soil (58, 59). Bacterial engulfment has also been described during colonization by the intestinal protist Entamoeba histolytica (60). Interestingly, a previous study showed that the presence of Entamoeba (likely including E. histolytica and E. dispar) was associated with an effect on Prevotellaceae and Ruminococcaceae that was similar to what we found in this study (33).

In the absence of Blastocystis, a strong bacterial competitor like Prevotella copri may dominate the community, which limits species richness and community evenness, as was observed in this study (Fig. 1). When present, Blastocystis predation on the most abundant bacterial taxa (also known as density-dependent predation) may lower the competition for nutrients and space, which may lead to an increase in bacterial richness and community evenness, as has been shown in this and previous studies (2, 61). This is well supported by food web theory and examples from macro- and microecology (62, 63), where the increase in community diversity through grazing or predation occurs through a top-down control on the strongest competitors, which consequently allows for the colonization and persistence of weaker competitors in the community. Our results do not eliminate the possibility of a bottom-up control, where specific bacterial community structures may favor colonization of the intestine by Blastocystis. However, this alternative explanation, where a more diverse community could favor the entry of a predator into the ecosystem, is not well supported by ecology theory (64). Future studies considering both bottom-up and top-down mechanisms are necessary to test various scenarios, including a deterministic influence of competition and predation on the gut microbial communities, which are impossible to assess in an observational cross-sectional study like this one.

Other potential indirect mechanisms for Blastocystis to affect the abundance and diversity of the bacterial microbiome are via direct interaction with the intestinal epithelium and the underlying immune tissue. Indeed, our study revealed differences in calprotectin and IgA. IgA-mediated modulation of the bacterial microbiome has been reported in mice lacking IgA (65, 66) or when comparing low-IgA mice to high-IgA mice (67). However, the result observed for IgA in this study was subtle and unlikely to cause such drastic changes in microbiome composition and function. Likewise, the levels detected for fecal calprotectin (a granulocyte marker protein produced at high levels in the cytosol of neutrophils and a well-known marker of intestinal inflammation [40]) are well below levels consistent with clinical inflammation (68), and Blastocystis was associated with lower levels. Still, indirect mechanisms need to be ruled out experimentally before concluding that they are not at play in the results obtained in this study.

These subclinical immune differences may also be a result of (i) the large differences in bacterial community structures observed in relation to Blastocystis colonization and/or (ii) a direct anti-inflammatory effect of Blastocystis on the intestinal mucosa. Evidence of the former mechanism clearly exists. Changes in IgA production have been observed in cohousing and fecal transfer experiments carried out in mice (67), where high-IgA mice that adopted the microbiota from low-IgA mice experienced a decrease in IgA levels, demonstrating the IgA modulatory capacity of the bacterial microbiome (67). There is also evidence of direct effects of Blastocystis on the immune system. Blastocystis is associated with reduced neutrophil counts in blood (69) and is known to produce serine proteases that degrade secretory IgA (sIgA) (70). These mechanisms must be further studied to understand if Blastocystis creates an anti-inflammatory environment in the intestine and if this may explain the asymptomatic status in this population. Ideally, this could be studied in a longitudinal cohort study that includes participants that transition from a noncolonized to a colonized status.

The large taxonomic differences observed were, unsurprisingly, accompanied by broad functional differences, inferred computationally by PICRUSt and measured directly through the determination of fecal SCFA concentration. Interestingly, the difference in the bacterial microbiome associated with Blastocystis colonization resulted in an overall decrease in carbohydrate fermentative metabolism. We attribute this to (i) the decrease in abundance of Prevotella copri, an efficient fiber degrader and acetate producer (71), and (ii) the reduced cross-feeding by other bacterial taxa due to the decrease in Prevotella copri. Cross-feeding, whereby a species produces a metabolite that can be used by other species, is a common feature in many microbial communities (72). In the mammalian gut, fiber is degraded into hexoses and pentoses, which are subsequently used as substrates in pathways that lead to the production of acetate, propionate, and butyrate (71). Butyrate can also be produced from acetate (via acetyl coenzyme A [acetyl-CoA]) (71), which could explain the larger decrease in butyrate concentrations in Blastocystis-colonized individuals.

Our study found similarities with previous reports. For example, Audebert et al. also reported increases in bacterial alpha diversity and in the abundance of Ruminococcus species in association with Blastocystis colonization (2). However, our study found taxonomic discrepancies with previous reports. Previously, O’Brien Andersen et al. (34) had reported an increase in Prevotella species associated with Blastocystis in a cohort of Danish patients referred to their study by doctors. Likewise, Audebert et al. (2) reported a positive association between Prevotellaceae and Blastocystis colonization in French individuals. These discrepancies may arise from comparing westernized populations, expected to have different microbiomes than rural Mexicans. Members of the Prevotellaceae family account for ~30% of the total gut bacterial taxa in this study, compared to ~5% of the French population studied previously (2). A study that compared the microbiota in relation to Blastocystis using shotgun metagenomics data from 11 studies across four continents found an association between Blastocysis, Prevotella copri, and Ruminococcus gnavus that was the opposite of what we observed in this study (7) and did not detect a difference in bacterial alpha diversity. These differences may be due to comparing our results to studies that have analyzed a combination of healthy and diseased individuals or because of differences in Blastocystis subtypes (2, 34). For example, while colonization was mostly explained by Blastocystis subtype 3 in this study, the samples analyzed by Audebert et al. included colonization with many more subtypes, including ST1, -2, -3, -4, -6, and -7 (2). Notably, an ongoing study in asymptomatic individuals from Cameroon has revealed findings similar to those of our study in relation to Blastocystis colonization (Laure Segurel, personal communication), suggesting that the effect of Blastocystis on the bacterial microbiome may be dependent on the sociogeographical setting studied or on similarities in subtype colonization. The important interstudy discrepancies reflect the need to study Blastocystis in a context-specific manner, accounting for the effects of subtypes, inflammation, geographical setting, and other covariates that may influence microbiome composition and diversity.

Altogether, our work is the first to show the important ecological association between the most common human-associated protist and the fundamental species of the bacterial microbiome in the absence of gastrointestinal disease or inflammation. While our findings are in agreement with the well-supported top-down effect of a predator on ecosystem diversity, this remains to be tested experimentally. If proven, the presence of Blastocystis or similar bacterial predators may be an important mechanism of microbial diversity maintenance that is not currently being considered in the human microbiome field. Including microbial eukaryotes in human gut microbiome surveys will provide a more integrated understanding of the ecological processes that shape the gut microbiome and the mechanisms by which they relate to health and disease.

MATERIALS AND METHODS

Study population and study design.This cross-sectional study was carried out in Xoxocotla, Morelos, Mexico, about 120 km south of Mexico City (longitude 99°19′W, latitude 18°3′N). This semirural area spans 29,917 km2 and hosts a total population of 5,163 people. The region is characterized by tropical climate (warm subhumid), and agriculture is the main source of economic income. Sample collection was carried out between January and September 2014. A total of 156 urine and stool samples were collected from individuals at Xoxocotla Health Center as a part of a routine medical examination of people living in the town or within the vicinity. Urine and stool samples were collected in sterile plastic containers, immediately placed at 4°C for transport, and then stored at −20°C until analysis. Prior to sample and data collection, the nature of the study was explained to the participants and informed verbal and written consents were obtained. Demographic, socioeconomic, and environmental data and the history of gastrointestinal symptoms in the 3 months preceding sampling (abdominal pain or discomfort and defecation patterns) were obtained based on the standardized Rome III questionnaire in face-to-face interviews. The Rome classification system is based on the symptom clusters that remain consistent across clinical and population groups and is currently the gold standard to determine gut health (73). The protocol was approved by the Medical Ethics Committee of the National Autonomous University of Mexico.

Real-time qPCR analysis for detection of Blastocystis and other parasites.The presence of Blastocystis, Cryptosporidium parvum, Entamoeba histolytica/E. dispar, and Giardia intestinalis was assessed by quantitative PCR (qPCR) using QuantiTect SYBR green master mix (Qiagen). The Blastocystis-specific primer BhRDr (5′ GAG CTT TTT AAC TGC AAC AAC G 3′ [74]) and the broad-specificity eukaryote-specific primer RD5 (5′ ATC TGG TTG ATC CTG CCA GT 3′ [75]), the Cryptosporidium-specific primers CrF (forward [F], 5′ CGC TTC TCT AGC CTT TCA TGA 3′) and CrR (reverse [R], 5′ CTT CAC GTG TGT TTG CCA AT 3′), the Entamoeba histolytica/E. dispar-specific primers Ehd-239F (5′ ATT GTC GTG GCA TCC TAA CTC A 3′) and Ehd-88R (5′ GCG GAC GGC TCA TTA TAA CA 3′), and the Giardia-specific primers Giardia-80F (F, 5′ GAC GGC TCA GGA CAA CGG TT 3′) and Giardia-127R (R, 5′ TTG CCA GCG GTG TCC G 3′) (76) were used. Amplification reactions were performed in 10-µl reaction mixture volumes with the Taq PCR master mix kit (Qiagen) and 6.25 pmol each of primers BhRDr-RD5, CrF-CrR, Ehd-239F–Ehd-88R, and Giardia-80F–Giardia-127R. The amplification conditions consisted of 35 cycles of 1 min each at 94, 59, and 72°C, with an additional step of 95°C for 15 s, 60°C for 1 min, 95°C for 30 s, and 60°C for 15 s (77). The qPCR plates included positive controls (samples known to be positive for each parasite), as well as standard curves using DNA from each parasite from ATCC’s enteric protozoon DNA panel.

The 18S rRNA gene was amplified using primers 5′ GTA CAC ACC GCC CGT C 3′ (F) and 5′ TGA TCC TTC TGC AGG TTC ACC TAC 3′ (R). The amplification conditions consisted of 35 cycles of 1 min each at 94, 59, and 72°C, with an additional step consisting of 95°C for 15 s, 60°C for 1 min, 95°C for 30 s, and 60°C for 15 s. All qPCRs were performed on an Applied Biosystems 7500 machine. The parasitic loads in the samples were calculated as the difference between the average cycle threshold (CT) value of each parasite and the average CT value of the 18S rRNA gene of each sample.

Fecal microbial community analysis. (i) Extraction.DNA was extracted from ~50 mg of stool. Samples were mechanically lysed using Mo Bio dry bead tubes (Mo Bio Laboratories, Inc.) and the FastPrep homogenizer (FastPrep instrument; MP Biochemicals) before DNA extraction with the Qiagen DNA stool minikit.

(ii) Amplification.For 16S rRNA gene amplification, samples were amplified by PCR in triplicate using bar-coded primer pairs flanking the V3 region of the 16S rRNA gene as previously described (78). Each 50 ml of PCR mixture contained 22 ml of water, 25 ml of TopTaq master mix, 0.5 ml of each forward and reverse bar-coded primer, and 2 ml of template DNA. The PCR program consisted of an initial DNA denaturation step at 95°C (5 min), 25 cycles of DNA denaturation at 95°C (1 min), an annealing step at 50°C (1 min), an elongation step at 72°C (1 min), and a final elongation step at 72°C (7 min). Controls without template DNA were included to ensure that no contamination occurred. Amplicons were run on a 2% agarose gel to ensure adequate amplification. Amplicons displaying bands at ~160 bp were purified using the illustra GFX PCR DNA purification kit. Purified samples were diluted 1:50 and quantified using PicoGreen (Invitrogen) in the Tecan M200 plate reader (excitation at 480 nm and emission at 520 nm).

(iii) Sequencing.For 16S rRNA gene sequencing, each PCR pool was analyzed on the Agilent Bioanalyzer using the high-sensitivity double-stranded DNA (dsDNA) assay to determine approximate library fragment size and verify library integrity. Pooled-library concentrations were determined using the TruSeq DNA sample preparation kit, version 2 (Illumina). Library pools were diluted to 4 nM and denatured into single strands using fresh 0.2 N NaOH. The final library loading concentration was 8 pM, with an additional PhiX spike-in of 20%. Sequencing was carried out using a Hi-Seq 2000 bidirectional Illumina sequencing and cluster kit, version 4 (Macrogen, Inc.). The 18S rRNA gene was sent to the Integrated Microbiome Resource at Dalhousie University for amplification and sequencing. The 18S rRNA gene was amplified with the primers E572F (5′ CYGCGGTAATTCCAGCTC 3′) and E1009R (5′ AYGGTATCTRATCRTCTTYG 3′), and the reaction mixture included a PNA blocking primer (5′ TCTTAATCATGGCCTCAGTT 3′) to reduce amplification of mammalian sequences. Amplification was carried out in duplicate, with one reaction mixture using undiluted DNA and the other using DNA diluted 1:10 in PCR water. Amplification was conducted according to previously described protocols (79). PCR products were visualized on E-gels, quantified using Invitrogen Qubit with PicoGreen, and pooled at equal concentrations, according to a previous report (79). PhiX was spiked in at 5%, and the resulting library was sequenced at Dalhousie University on the Illumina MiSeq using the MiSeq 500-cycle reagent kit, version 2 (250 × 2).

Bioinformatics.Sequences were preprocessed, denoised, and quality filtered by size using the Mothur MiSeq SOP (standard operation protocol) (16S rRNA gene [80]) or QIIME (18S rRNA gene [81]).

(i) 16S rRNA gene sequences.All sequences were processed using Mothur according to the standard operating procedure as previously described (82). Quality sequences were obtained by removing sequences with ambiguous bases, a low-quality read length, and/or chimeras identified using chimera.uchime. Quality sequences were aligned to the SILVA bacterial reference alignment, and OTUs were generated using a dissimilarity cutoff of 0.03. Sequences were classified using the classify.seqs command.

(ii) 18S rRNA gene sequences.Demultiplexed reads were trimmed to a uniform length of 250 bp using the FastX-Toolkit (http://hannonlab.cshl.edu/fastx_toolkit/) and clustered into operational taxonomic units (OTUs) using the minimum entropy decomposition (MED) method (83) as implemented in the oligotyping microbial analysis software package (84). MED performs de novo taxonomic clustering using Shannon entropy to separate biologically meaningful patterns of nucleotide diversity from sequencing noise; the processed data are partitioned into phylogenetically homogeneous units (MED nodes) for downstream bacterial diversity analyses. This analysis was carried out with the minimum substantive abundance parameter (−M) set at 250 reads. All other parameters were run with default settings; the maximum variation allowed per node (−V) was automatically set at 3 nucleotides.

Representative sequences were classified by clustering against the Greengenes Database at 97% similarity (16S rRNA gene [85]) or SILVA release 123 at 99% similarity (18S rRNA gene [86]). The 16S rRNA gene data set was filtered to remove mitochondrion and chloroplast sequences and OTUs present in fewer than three samples. The 18S rRNA gene data set was filtered to remove mammalian and plant sequences and all OTUs present in fewer than three samples. Both data sets were filtered to exclude singletons and doubletons. Following filtering, a cutoff of 1,000 reads per sample was applied. All 16S rRNA gene samples passed the cutoff, while 63 samples were excluded from the 18S rRNA gene data set.

PICRUSt.We used PICRUSt (39) to generate a profile of putative functions (via metagenome prediction) from the 16S rRNA OTU data. Predicted metagenomes from all the samples were categorized by function at KEGG Orthology level 3. To test for significant differences in functional category abundances between colonized and noncolonized individuals, we used the Welch’s t test implementation of STAMP (87).

SCFA analysis.Fecal samples were combined with 25% phosphoric acid, vortexed, and centrifuged until a clear supernatant was obtained. Supernatants were submitted for GC analysis to the Department of Agricultural, Food and Nutritional Science of the University of Alberta. Samples were analyzed as previously described (88), with modifications. Briefly, samples were combined with 4-methyl-valeric acid as an internal standard, and 0.2 ml was injected into the Bruker Scion 456 gas chromatograph, using a Stabilwax-DA 30-m by 0.53-mm by 0.5-µm column (Restek). A standard solution containing acetic acid, propionic acid, isobutyric acid, butyric acid, isovaleric acid, valeric acid, and caproic acid combined with internal standard was injected in every run.

The PTV (programmable temperature vaporization) injector and FID (flame ionization detector) detector temperatures were held at 250°C for the entire run. The oven was started at 80°C and immediately ramped to 210°C at 45°C/min, where it was held for 5.11 min. Total run time was 8.00 min. Helium was used at a constant flow of 20.00 ml/min. Sample concentrations were normalized to the wet weight of feces.

Calprotectin assay.Fecal calprotectin was determined using a human calprotectin enzyme-linked immunosorbent assay (ELISA) kit (Hycult Biotech, Inc., Uden, Netherlands). Stool extracts were prepared and analyzed according to the manufacturer’s instructions. After normalizing fecal weight to 1-mg/µl solutions, samples were diluted 10× in phosphate-buffered saline (PBS). The standard curve ranged between 1.6 and 50 ng/ml.

IgA assay.Fecal IgA concentrations were measured by enzyme-linked immunosorbent assay (Chondrex) in samples diluted 1:1,000. After normalizing fecal weight to 1-mg/µl solutions, samples were diluted 1,000× in PBS. The standard curve ranged between 1.6 and 50 ng/ml.

Oxidative stress biomarkers in urine.Biochemical assays for AOPP and MDA were carried out on collected urine specimens according to methods previously established (13, 47). AOPP was measured spectrophotometrically according to a previously described method (89), using a microplate reader (Tecan Infinite M200; Tecan, Switzerland) with the following modifications. Briefly, 200-µl amounts of urine diluted 1:5 in PBS or chloramine-T standard solutions (0 to 100 µmol/liter) were mixed with 20 µl of acetic acid and 10 µl of 1.16 M potassium iodide (KI; Sigma). The absorbance was immediately read at 340 nm against a blank containing 200 µl of PBS, 20 µl of acetic acid, and 10 µl of KI. AOPP concentrations were expressed in µmol/liter of chloramine-T equivalents.

The free radical-induced lipid peroxidation level in urine was determined by measuring MDA with an assay modified from the thiobarbituric acid-reactive-substance method (90). Briefly, 600-µl amounts of urine or 1,1,3,3-tetraethoxypropane standard solutions (0.15, 0.30, 0.60, 1.20, 2.40, 5.0, and 10.0 nmol/ml) were mixed with 300 µl of 37% hydrochloric acid. After centrifugation (800 × g for 10 min), 500 µl of 0.65% N-methyl-2-phenylindole (MPI) (acetonitrile-methanol [3:1] diluent; Sigma) and 500 µl of deproteinized serum were added, followed by a 45-min incubation at 45°C. The resulting blue-colored chromophore was measured at an absorbance of 586 nm in a microplate reader (Tecan Infinite M200; Tecan, Switzerland) using 1,1,3,3-tetraethoxypropane as a standard and water as the blank. Concentrations were expressed as nmol/ml.

Serum IL measurement.Serum was extracted from 5 ml of peripheral blood of each individual to measure IL-2, IL-4, IL-6, IL-10, TNF, IFN-γ, and IL-17A protein levels according to the manufacturer’s recommendations for the BD cytometric bead array (CBA) human Th1/Th2/Th17 cytokine kit (BD Biosciences, San Jose CA, USA).

Statistical analysis.Differences in frequencies for categorical and continuous variables between cases and controls were evaluated using the chi-squared and Student’s t test, respectively. We assessed fecal microbial diversity and the relative abundances of bacterial and eukaryotic taxa using phyloseq (91), along with additional R-based computational tools (92–97). PCoAs were conducted using phyloseq (Bray-Curtis dissimilarities as distance metric) on both variance-stabilizing-transformed and rarefied OTU matrices and then statistically confirmed by a permutational multivariate analysis of variance (PERMANOVA) to confirm our results were not a consequence of heteroscedastic dispersion between groups (38). The Shannon and Chao1 alpha diversity indices were calculated using phyloseq and statistically confirmed by the Mann-Whitney test (GraphPad Prism software, version 5c). The R packages DESeq2 (98) and MaAsLin (35) were used to calculate differentially abundant OTUs. Correlation analysis was performed using the bicor method in the R package microbiome to correlate the 100 most abundant OTUs from the 16S and 18S rRNA gene data sets. Features in the analysis were included as OTUs and as OTUs combined into taxonomic families.

To confirm that the effects of Blastocystis on the gut bacterial community structure were not due to the presence and interaction of other protozoans or influenced by variation in sequencing depth, we used PERMANOVA models. We built three models, one each for Blastocystis with Hymenolepis nana, Blastocystis with Debaryomyces hansenii, and Blastocystis with sequencing. Neither protozoan interacted significantly with the effect of Blastocystis on the gut bacterial community structure, and Blastocystis remained a significant driver of community structure in both models. For the 18S rRNA gene data set, a PERMANOVA on gut eukaryotic community structure shows that the sequencing depth appears to interact significantly with the Blastocystis species colonization effect on community structure. This result is biologically intuitive, since the presence of Blastocystis species in itself is contributing to the number of sequences detected. However, to further control for this result, we carried out the DESeq2 analysis using variance-stabilizing-transformed data. This analysis yielded the same number and identity of differential OTUs between Blastocystis-colonized and -noncolonized individuals as we obtained with nontransformed data.

The statistical analyses for calprotectin, IgA, serum cytokines, SCFA, and urine oxidative stress markers were performed with GraphPad Prism software, version 6.00 (GraphPad Software, Inc., San Diego, CA, USA). For comparison between two groups, Student’s t tests and Mann-Whitney tests were used for normally and abnormally distributed data sets, respectively, and no samples were excluded from statistical analysis. Statistically significant differences were set at a P value of ≤0.05. A P value of >0.05 was considered not significant.

Data availability.Sequences are publicly available at https://doi.org/10.6084/m9.figshare.6359306. The R code supporting the results of this article is available in the Xoxocotla Gut Microbiome project repository on figshare at https://figshare.com/s/a99310b6d9be4e9de4fe.

ACKNOWLEDGMENTS

We thank all study participants from Xoxocotla as well as the personnel involved in this cohort.

Funding for this study was obtained from grant numbers 272601, 283522, and 257091 from the National Council of Sciences and Technology in Mexico (CONACyT), IN218214 from PAPIIT UNAM, and FIS/IMSS/PROT/1368 from IMSS to C.X.-G. and J.T., as well as Canadian Health Research Institutes grants to B.B.F. and a Human Frontier Science Program grant to L.W.P. (grant number RGY0078/2015). M.E.N.-R. (proposal number 235618) and O.P.-R. (proposal number 208253) each received a 1-year scholarship from the Estancias Posdoctorales al Extranjero para la Consolidacion de Grupos de Investigación program of CONACyT.

M.E.N.-R., O.P.-R., M.C.A., C.X.-G., and B.B.F. designed the study. A.V.-S., P.M.-S., J.T, L.R.-V., and C.X.-G. coordinated and facilitated the cohort study in Xoxocotla. P.M.-S. and C.X.-G. conducted medical examinations and Rome III questionnaires. A.V.-S. curated the database and metadata. M.C.A. and L.W.P. optimized sequencing strategies. M.E.N.-R. prepared samples for sequencing analysis. M.C.A. and E.M. performed the bioinformatics analysis of sequencing data. M.E.N.-R. performed oxidative stress marker and qPCR assays. O.P.-R. and E.M.B. conducted IgA and calprotectin assays. L.A.R. performed cytokine assays. M.J. and J.W. conducted SCFA assays and analysis. M.-C.A. and I.L.-L. designed and performed statistical analyses and created figures for the paper. M.E.N.-R., I.L.-L. and M.C.A. wrote the paper. M.E.N.-R., I.L.-L., L.W.P, J.W., B.B.F., C.X.-G., and M.C.A. edited the paper. All authors contributed extensively to the work presented in the paper.

The authors declare that this study complies with the current laws of the countries in which the sample collection was performed. The authors declare that they have no competing interest.

FOOTNOTES

    • Received January 23, 2018.
    • Accepted May 2, 2018.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

REFERENCES

  1. 1.↵
    1. Tan KS,
    2. Mirza H,
    3. Teo JD,
    4. Wu B,
    5. Macary PA
    . 2010. Current views on the clinical relevance of Blastocystis spp. Curr Infect Dis Rep 12:28–35. doi:10.1007/s11908-009-0073-8.
    OpenUrlCrossRefPubMed
  2. 2.↵
    1. Audebert C,
    2. Even G,
    3. Cian A, Blastocystis Investigation Group,
    4. Loywick A,
    5. Merlin S,
    6. Viscogliosi E,
    7. Chabe M
    . 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.
    OpenUrlCrossRef
  3. 3.↵
    1. Bart A,
    2. Wentink-Bonnema EM,
    3. Gilis H,
    4. Verhaar N,
    5. Wassenaar CJ,
    6. van Vugt M,
    7. Goorhuis A,
    8. van Gool T
    . 2013. Diagnosis and subtype analysis of Blastocystis sp. in 442 patients in a hospital setting in the Netherlands. BMC Infect Dis 13:389. doi:10.1186/1471-2334-13-389.
    OpenUrlCrossRefPubMed
  4. 4.↵
    1. Clark CG,
    2. van der Giezen M,
    3. Alfellani MA,
    4. Stensvold CR
    . 2013. Recent developments in Blastocystis research. Adv Parasitol 82:1–32. doi:10.1016/B978-0-12-407706-5.00001-0.
    OpenUrlCrossRefPubMed
  5. 5.↵
    1. El Safadi D,
    2. Gaayeb L,
    3. Meloni D,
    4. Cian A,
    5. Poirier P,
    6. Wawrzyniak I,
    7. Delbac F,
    8. Dabboussi F,
    9. Delhaes L,
    10. Seck M,
    11. Hamze M,
    12. Riveau G,
    13. Viscogliosi E
    . 2014. Children of Senegal River Basin show the highest prevalence of Blastocystis sp. ever observed worldwide. BMC Infect Dis 14:164. doi:10.1186/1471-2334-14-164.
    OpenUrlCrossRefPubMed
  6. 6.↵
    1. Wawrzyniak I,
    2. Poirier P,
    3. Viscogliosi E,
    4. Dionigia M,
    5. Texier C,
    6. Delbac F,
    7. Alaoui HE
    . 2013. Blastocystis, an unrecognized parasite: an overview of pathogenesis and diagnosis. Ther Adv Infect Dis 1:167–178. doi:10.1177/2049936113504754.
    OpenUrlCrossRefPubMed
  7. 7.↵
    1. Beghini F,
    2. Pasolli E,
    3. Truong TD,
    4. Putignani L,
    5. Cacciò SM,
    6. Segata N
    . 2017. Large-scale comparative metagenomics of Blastocystis, a common member of the human gut microbiome. ISME J 11:2848–2863. doi:10.1038/ismej.2017.139.
    OpenUrlCrossRef
  8. 8.↵
    1. Tan KS
    . 2008. New insights on classification, identification, and clinical relevance of Blastocystis spp. Clin Microbiol Rev 21:639–665. doi:10.1128/CMR.00022-08.
    OpenUrlAbstract/FREE Full Text
  9. 9.↵
    1. Krogsgaard LR,
    2. Engsbro AL,
    3. Stensvold CR,
    4. Nielsen HV,
    5. Bytzer P
    . 2015. The prevalence of intestinal parasites is not greater among individuals with irritable bowel syndrome: a population-based case-control study. Clin Gastroenterol Hepatol 13:507–513.e2. doi:10.1016/j.cgh.2014.07.065.
    OpenUrlCrossRef
  10. 10.↵
    1. Scanlan PD,
    2. Stensvold CR,
    3. Rajilić-Stojanović M,
    4. Heilig HG,
    5. De Vos WM,
    6. O’Toole PW,
    7. Cotter PD
    . 2014. The microbial eukaryote Blastocystis is a prevalent and diverse member of the healthy human gut microbiota. FEMS Microbiol Ecol 90:326–330. doi:10.1111/1574-6941.12396.
    OpenUrlCrossRefPubMed
  11. 11.↵
    1. Vargas-Sanchez GB,
    2. Romero-Valdovinos M,
    3. Ramirez-Guerrero C,
    4. Vargas-Hernandez I,
    5. Ramirez-Miranda ME,
    6. Martinez-Ocaña J,
    7. Valadez A,
    8. Ximenez C,
    9. Lopez-Escamilla E,
    10. Hernandez-Campos ME,
    11. Villalobos G,
    12. Martinez-Hernandez F,
    13. Maravilla P
    . 2015. Blastocystis isolates from patients with irritable bowel syndrome and from asymptomatic carriers exhibit similar parasitological loads, but significantly different generation times and genetic variability across multiple subtypes. PLoS One 10:e0124006. doi:10.1371/journal.pone.0124006.
    OpenUrlCrossRef
  12. 12.↵
    1. Bálint A,
    2. Dóczi I,
    3. Bereczki L,
    4. Gyulai R,
    5. Szűcs M,
    6. Farkas K,
    7. Urbán E,
    8. Nagy F,
    9. Szepes Z,
    10. Wittmann T,
    11. Molnár T
    . 2014. Do not forget the stool examination!—cutaneous and gastrointestinal manifestations of Blastocystis sp. infection. Parasitol Res 113:1585–1590. doi:10.1007/s00436-014-3805-0.
    OpenUrlCrossRef
  13. 13.↵
    1. Chandramathi S,
    2. Suresh K,
    3. Shuba S,
    4. Mahmood A,
    5. Kuppusamy UR
    . 2010. High levels of oxidative stress in rats infected with Blastocystis hominis. Parasitology 137:605–611. doi:10.1017/S0031182009991351.
    OpenUrlCrossRefPubMed
  14. 14.↵
    1. Chandramathi S,
    2. Suresh K,
    3. Sivanandam S,
    4. Kuppusamy UR
    . 2014. Stress exacerbates infectivity and pathogenicity of Blastocystis hominis: in vitro and in vivo evidences. PLoS One 9:e94567. doi:10.1371/journal.pone.0094567.
    OpenUrlCrossRef
  15. 15.↵
    1. Lim MX,
    2. Png CW,
    3. Tay CY,
    4. Teo JD,
    5. Jiao H,
    6. Lehming N,
    7. Tan KS,
    8. Zhang Y
    . 2014. Differential regulation of proinflammatory cytokine expression by mitogen-activated protein kinases in macrophages in response to intestinal parasite infection. Infect Immun 82:4789–4801. doi:10.1128/IAI.02279-14.
    OpenUrlAbstract/FREE Full Text
  16. 16.↵
    1. Azizian M,
    2. Basati G,
    3. Abangah G,
    4. Mahmoudi MR,
    5. Mirzaei A
    . 2016. Contribution of Blastocystis hominis subtypes and associated inflammatory factors in development of irritable bowel syndrome. Parasitol Res 115:2003–2009. doi:10.1007/s00436-016-4942-4.
    OpenUrlCrossRef
  17. 17.↵
    1. Yakoob J,
    2. Abbas Z,
    3. Usman MW,
    4. Sultana A,
    5. Islam M,
    6. Awan S,
    7. Ahmad Z,
    8. Hamid S,
    9. Jafri W
    . 2014. Cytokine changes in colonic mucosa associated with Blastocystis spp. subtypes 1 and 3 in diarrhoea-predominant irritable bowel syndrome. Parasitology 141:957–969. doi:10.1017/S003118201300173X.
    OpenUrlCrossRef
  18. 18.↵
    1. Nagel R,
    2. Traub RJ,
    3. Kwan MM,
    4. Bielefeldt-Ohmann H
    . 2015. Blastocystis specific serum immunoglobulin in patients with irritable bowel syndrome (IBS) versus healthy controls. Parasit Vectors 8:453. doi:10.1186/s13071-015-1069-x.
    OpenUrlCrossRef
  19. 19.↵
    1. Ragavan ND,
    2. Kumar S,
    3. Chye TT,
    4. Mahadeva S,
    5. Shiaw-Hooi H
    . 2015. Blastocystis sp. in irritable bowel syndrome (IBS)—detection in stool aspirates during colonoscopy. PLoS One 10:e0121173. doi:10.1371/journal.pone.0121173.
    OpenUrlCrossRef
  20. 20.↵
    1. Parfrey LW,
    2. Walters WA,
    3. Lauber CL,
    4. Clemente JC,
    5. Berg-Lyons D,
    6. Teiling C,
    7. Kodira C,
    8. Mohiuddin M,
    9. Brunelle J,
    10. Driscoll M,
    11. Fierer N,
    12. Gilbert JA,
    13. Knight R
    . 2014. Communities of microbial eukaryotes in the mammalian gut within the context of environmental eukaryotic diversity. Front Microbiol 5:298. doi:10.3389/fmicb.2014.00298.
    OpenUrlCrossRefPubMed
  21. 21.↵
    1. David LA,
    2. Maurice CF,
    3. Carmody RN,
    4. Gootenberg DB,
    5. Button JE,
    6. Wolfe BE,
    7. Ling AV,
    8. Devlin AS,
    9. Varma Y,
    10. Fischbach MA,
    11. Biddinger SB,
    12. Dutton RJ,
    13. Turnbaugh PJ
    . 2014. Diet rapidly and reproducibly alters the human gut microbiome. Nature 505:559–563. doi:10.1038/nature12820.
    OpenUrlCrossRefPubMedWeb of Science
  22. 22.↵
    1. Wu GD,
    2. Chen J,
    3. Hoffmann C,
    4. Bittinger K,
    5. Chen YY,
    6. Keilbaugh SA,
    7. Bewtra M,
    8. Knights D,
    9. Walters WA,
    10. Knight R,
    11. Sinha R,
    12. Gilroy E,
    13. Gupta K,
    14. Baldassano R,
    15. Nessel L,
    16. Li H,
    17. Bushman FD,
    18. Lewis JD
    . 2011. Linking long-term dietary patterns with gut microbial enterotypes. Science 334:105–108. doi:10.1126/science.1208344.
    OpenUrlAbstract/FREE Full Text
  23. 23.↵
    1. Keeney KM,
    2. Yurist-Doutsch S,
    3. Arrieta MC,
    4. Finlay BB
    . 2014. Effects of antibiotics on human microbiota and subsequent disease. Annu Rev Microbiol 68:217–235. doi:10.1146/annurev-micro-091313-103456.
    OpenUrlCrossRefPubMedWeb of Science
  24. 24.↵
    1. Yatsunenko T,
    2. Rey FE,
    3. Manary MJ,
    4. Trehan I,
    5. Dominguez-Bello MG,
    6. Contreras M,
    7. Magris M,
    8. Hidalgo G,
    9. Baldassano RN,
    10. Anokhin AP,
    11. Heath AC,
    12. Warner B,
    13. Reeder J,
    14. Kuczynski J,
    15. Caporaso JG,
    16. Lozupone CA,
    17. Lauber C,
    18. Clemente JC,
    19. Knights D,
    20. Knight R,
    21. Gordon JI
    . 2012. Human gut microbiome viewed across age and geography. Nature 486:222–227. doi:10.1038/nature11053.
    OpenUrlCrossRefPubMedWeb of Science
  25. 25.↵
    1. Slingerland AE,
    2. Schwabkey Z,
    3. Wiesnoski DH,
    4. Jenq RR
    . 2017. Clinical evidence for the microbiome in inflammatory diseases. Front Immunol 8:400. doi:10.3389/fimmu.2017.00400.
    OpenUrlCrossRef
  26. 26.↵
    1. Goodrich JK,
    2. Davenport ER,
    3. Beaumont M,
    4. Jackson MA,
    5. Knight R,
    6. Ober C,
    7. Spector TD,
    8. Bell JT,
    9. Clark AG,
    10. Ley RE
    . 2016. Genetic determinants of the gut microbiome in UK twins. Cell Host Microbe 19:731–743. doi:10.1016/j.chom.2016.04.017.
    OpenUrlCrossRefPubMed
  27. 27.↵
    1. Goodrich JK,
    2. Waters JL,
    3. Poole AC,
    4. Sutter JL,
    5. Koren O,
    6. Blekhman R,
    7. Beaumont M,
    8. Van Treuren W,
    9. Knight R,
    10. Bell JT,
    11. Spector TD,
    12. Clark AG,
    13. Ley RE
    . 2014. Human genetics shape the gut microbiome. Cell 159:789–799. doi:10.1016/j.cell.2014.09.053.
    OpenUrlCrossRefPubMedWeb of Science
  28. 28.↵
    1. Zhernakova A,
    2. Kurilshikov A,
    3. Bonder MJ,
    4. Tigchelaar EF,
    5. Schirmer M,
    6. Vatanen T,
    7. Mujagic Z,
    8. Vila AV,
    9. Falony G,
    10. Vieira-Silva S,
    11. Wang J,
    12. Imhann F,
    13. Brandsma E,
    14. Jankipersadsing SA,
    15. Joossens M,
    16. Cenit MC,
    17. Deelen P,
    18. Swertz MA, Lifelines cohort study,
    19. Weersma RK,
    20. Feskens EJ,
    21. Netea MG,
    22. Gevers D,
    23. Jonkers D,
    24. Franke L,
    25. Aulchenko YS,
    26. Huttenhower C,
    27. Raes J,
    28. Hofker MH,
    29. Xavier RJ,
    30. Wijmenga C,
    31. Fu J
    . 2016. Population-based metagenomics analysis reveals markers for gut microbiome composition and diversity. Science 352:565–569. doi:10.1126/science.aad3369.
    OpenUrlAbstract/FREE Full Text
  29. 29.↵
    1. Barash NR,
    2. Maloney JG,
    3. Singer SM,
    4. Dawson SC
    . 2017. Giardia alters commensal microbial diversity throughout the murine gut. Infect Immun 85:e00948-16. doi:10.1128/IAI.00948-16.
    OpenUrlAbstract/FREE Full Text
  30. 30.↵
    1. Beatty JK,
    2. Akierman SV,
    3. Motta JP,
    4. Muise S,
    5. Workentine ML,
    6. Harrison JJ,
    7. Bhargava A,
    8. Beck PL,
    9. Rioux KP,
    10. McKnight GW,
    11. Wallace JL,
    12. Buret AG
    . 2017. Giardia duodenalis induces pathogenic dysbiosis of human intestinal microbiota biofilms. Int J Parasitol 47:311–326. doi:10.1016/j.ijpara.2016.11.010.
    OpenUrlCrossRef
  31. 31.↵
    1. Burgess SL,
    2. Petri WA, Jr.
    2016. The intestinal bacterial microbiome and E. histolytica infection. Curr Trop Med Rep 3:71–74. doi:10.1007/s40475-016-0083-1.
    OpenUrlCrossRef
  32. 32.↵
    1. Verma AK,
    2. Verma R,
    3. Ahuja V,
    4. 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.
    OpenUrlCrossRefPubMed
  33. 33.↵
    1. Morton ER,
    2. Lynch J,
    3. Froment A,
    4. Lafosse S,
    5. Heyer E,
    6. Przeworski M,
    7. Blekhman R,
    8. Ségurel L
    . 2015. Variation in rural African gut microbiota is strongly correlated with colonization by Entamoeba and subsistence. PLoS Genet 11:e1005658. doi:10.1371/journal.pgen.1005658.
    OpenUrlCrossRefPubMed
  34. 34.↵
    1. O’Brien Andersen L,
    2. Karim AB,
    3. Roager HM,
    4. Vigsnæs LK,
    5. Krogfelt KA,
    6. Licht TR,
    7. Stensvold CR
    . 2016. Associations between common intestinal parasites and bacteria in humans as revealed by qPCR. Eur J Clin Microbiol Infect Dis 35:1427–1431. doi:10.1007/s10096-016-2680-2.
    OpenUrlCrossRef
  35. 35.↵
    1. Morgan XC,
    2. Tickle TL,
    3. Sokol H,
    4. Gevers D,
    5. Devaney KL,
    6. Ward DV,
    7. Reyes JA,
    8. Shah SA,
    9. LeLeiko N,
    10. Snapper SB,
    11. Bousvaros A,
    12. Korzenik J,
    13. Sands BE,
    14. Xavier RJ,
    15. Huttenhower C
    . 2012. Dysfunction of the intestinal microbiome in inflammatory bowel disease and treatment. Genome Biol 13:R79. doi:10.1186/gb-2012-13-9-r79.
    OpenUrlCrossRefPubMed
  36. 36.↵
    1. Weiss S,
    2. Xu ZZ,
    3. Peddada S,
    4. Amir A,
    5. Bittinger K,
    6. Gonzalez A,
    7. Lozupone C,
    8. Zaneveld JR,
    9. Vázquez-Baeza Y,
    10. Birmingham A,
    11. Hyde ER,
    12. Knight R
    . 2017. Normalization and microbial differential abundance strategies depend upon data characteristics. Microbiome 5:27. doi:10.1186/s40168-017-0237-y.
    OpenUrlCrossRef
  37. 37.↵
    1. Desnos-Ollivier M,
    2. Ragon M,
    3. Robert V,
    4. Raoux D,
    5. Gantier JC,
    6. Dromer F
    . 2008. Debaryomyces hansenii (Candida famata), a rare human fungal pathogen often misidentified as Pichia guilliermondii (Candida guilliermondii). J Clin Microbiol 46:3237–3242. doi:10.1128/JCM.01451-08.
    OpenUrlAbstract/FREE Full Text
  38. 38.↵
    1. McMurdie PJ,
    2. Holmes S
    . 2014. Waste not, want not: why rarefying microbiome data is inadmissible. PLoS Comput Biol 10:e1003531. doi:10.1371/journal.pcbi.1003531.
    OpenUrlCrossRefPubMed
  39. 39.↵
    1. Langille MG,
    2. Zaneveld J,
    3. Caporaso JG,
    4. McDonald D,
    5. Knights D,
    6. Reyes JA,
    7. Clemente JC,
    8. Burkepile DE,
    9. Vega Thurber RL,
    10. Knight R,
    11. Beiko RG,
    12. Huttenhower C
    . 2013. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat Biotechnol 31:814–821. doi:10.1038/nbt.2676.
    OpenUrlCrossRefPubMed
  40. 40.↵
    1. Walsham NE,
    2. Sherwood RA
    . 2016. Fecal calprotectin in inflammatory bowel disease. Clin Exp Gastroenterol 9:21–29. doi:10.2147/CEG.S51902.
    OpenUrlCrossRefPubMed
  41. 41.↵
    1. Gutzeit C,
    2. Magri G,
    3. Cerutti A
    . 2014. Intestinal IgA production and its role in host-microbe interaction. Immunol Rev 260:76–85. doi:10.1111/imr.12189.
    OpenUrlCrossRef
  42. 42.↵
    1. Johansen FE,
    2. Pekna M,
    3. Norderhaug IN,
    4. Haneberg B,
    5. Hietala MA,
    6. Krajci P,
    7. Betsholtz C,
    8. Brandtzaeg P
    . 1999. Absence of epithelial immunoglobulin A transport, with increased mucosal leakiness, in polymeric immunoglobulin receptor/secretory component-deficient mice. J Exp Med 190:915–922. doi:10.1084/jem.190.7.915.
    OpenUrlAbstract/FREE Full Text
  43. 43.↵
    1. McCoy KD,
    2. Stoel M,
    3. Stettler R,
    4. Merky P,
    5. Fink K,
    6. Senn BM,
    7. Schaer C,
    8. Massacand J,
    9. Odermatt B,
    10. Oettgen HC,
    11. Zinkernagel RM,
    12. Bos NA,
    13. Hengartner H,
    14. Macpherson AJ,
    15. Harris NL
    . 2008. Polyclonal and specific antibodies mediate protective immunity against enteric helminth infection. Cell Host Microbe 4:362–373. doi:10.1016/j.chom.2008.08.014.
    OpenUrlCrossRefPubMedWeb of Science
  44. 44.↵
    1. Bilgin R,
    2. Yalcin MS,
    3. Yucebilgic G,
    4. Koltas IS,
    5. Yazar S
    . 2012. Oxidative stress in vivax malaria. Korean J Parasitol 50:375–377. doi:10.3347/kjp.2012.50.4.375.
    OpenUrlCrossRef
  45. 45.↵
    1. Kilic E,
    2. Yazar S,
    3. Saraymen R,
    4. Ozbilge H
    . 2003. Serum malondialdehyde level in patients infected with Ascaris lumbricoides. World J Gastroenterol 9:2332–2334. doi:10.3748/wjg.v9.i10.2332.
    OpenUrlCrossRefPubMed
  46. 46.↵
    1. Namıduru ES,
    2. Tarakçıoğlu M,
    3. Namıduru M,
    4. Kocabaş R,
    5. Erbağcı B,
    6. Meram I,
    7. Karaoğlan I,
    8. Yılmaz N,
    9. Cekmen M
    . 2011. Increased serum nitric oxide and malondialdehyde levels in patients with acute intestinal amebiasis. Asian Pac J Trop Biomed 1:478–481. doi:10.1016/S2221-1691(11)60104-4.
    OpenUrlCrossRef
  47. 47.↵
    1. Chandramathi S,
    2. Suresh K,
    3. Anita ZB,
    4. Kuppusamy UR
    . 2009. Comparative assessment of urinary oxidative indices in breast and colorectal cancer patients. J Cancer Res Clin Oncol 135:319–323. doi:10.1007/s00432-008-0462-7.
    OpenUrlCrossRefPubMed
  48. 48.↵
    1. Cho I,
    2. Blaser MJ
    . 2012. The human microbiome: at the interface of health and disease. Nat Rev Genet 13:260–270. doi:10.1038/nrg3182.
    OpenUrlCrossRefPubMed
  49. 49.↵
    1. Zaiss MM,
    2. Rapin A,
    3. Lebon L,
    4. Dubey LK,
    5. Mosconi I,
    6. Sarter K,
    7. Piersigilli A,
    8. Menin L,
    9. Walker AW,
    10. Rougemont J,
    11. Paerewijck O,
    12. Geldhof P,
    13. McCoy KD,
    14. Macpherson AJ,
    15. Croese J,
    16. Giacomin PR,
    17. Loukas A,
    18. Junt T,
    19. Marsland BJ,
    20. Harris NL
    . 2015. The intestinal microbiota contributes to the ability of helminths to modulate allergic inflammation. Immunity 43:998–1010. doi:10.1016/j.immuni.2015.09.012.
    OpenUrlCrossRefPubMed
  50. 50.↵
    1. Gorvitovskaia A,
    2. Holmes SP,
    3. Huse SM
    . 2016. Interpreting Prevotella and Bacteroides as biomarkers of diet and lifestyle. Microbiome 4:15. doi:10.1186/s40168-016-0160-7.
    OpenUrlCrossRefPubMed
  51. 51.↵
    1. De Filippo C,
    2. Cavalieri D,
    3. Di Paola M,
    4. Ramazzotti M,
    5. Poullet JB,
    6. Massart S,
    7. Collini S,
    8. Pieraccini G,
    9. Lionetti P
    . 2010. Impact of diet in shaping gut microbiota revealed by a comparative study in children from Europe and rural Africa. Proc Natl Acad Sci U S A 107:14691–14696. doi:10.1073/pnas.1005963107.
    OpenUrlAbstract/FREE Full Text
  52. 52.↵
    1. Ze X,
    2. Ben David Y,
    3. Laverde-Gomez JA,
    4. Dassa B,
    5. Sheridan PO,
    6. Duncan SH,
    7. Louis P,
    8. Henrissat B,
    9. Juge N,
    10. Koropatkin NM,
    11. Bayer EA,
    12. Flint HJ
    . 2015. Unique organization of extracellular amylases into amylosomes in the resistant starch-utilizing human colonic Firmicutes bacterium Ruminococcus bromii. mBio 6:e01058-15. doi:10.1128/mBio.01058-15.
    OpenUrlAbstract/FREE Full Text
  53. 53.↵
    1. Ze X,
    2. Duncan SH,
    3. Louis P,
    4. Flint HJ
    . 2012. Ruminococcus bromii is a keystone species for the degradation of resistant starch in the human colon. ISME J 6:1535–1543. doi:10.1038/ismej.2012.4.
    OpenUrlCrossRefPubMedWeb of Science
  54. 54.↵
    1. Dunn LA,
    2. Boreham PF,
    3. Stenzel DJ
    . 1989. Ultrastructural variation of Blastocystis hominis stocks in culture. Int J Parasitol 19:43–56. doi:10.1016/0020-7519(89)90020-9.
    OpenUrlCrossRefPubMed
  55. 55.↵
    1. Singh M,
    2. Suresh K,
    3. Ho LC,
    4. Ng GC,
    5. Yap EH
    . 1995. Elucidation of the life cycle of the intestinal protozoan Blastocystis hominis. Parasitol Res 81:446–450. doi:10.1007/BF00931510.
    OpenUrlCrossRefPubMedWeb of Science
  56. 56.↵
    1. Lanuza MD,
    2. Carbajal JA,
    3. Villar J,
    4. Borrás R
    . 1997. Description of an improved method for Blastocystis hominis culture and axenization. Parasitol Res 83:60–63. doi:10.1007/s004360050209.
    OpenUrlCrossRefPubMed
  57. 57.↵
    1. Tan TC,
    2. Suresh KG
    . 2006. Predominance of amoeboid forms of Blastocystis hominis in isolates from symptomatic patients. Parasitol Res 98:189–193. doi:10.1007/s00436-005-0033-7.
    OpenUrlCrossRefPubMedWeb of Science
  58. 58.↵
    1. Clarholm M
    . 1981. Protozoan grazing of bacteria in soil-impact and importance. Microb Ecol 7:343–350. doi:10.1007/BF02341429.
    OpenUrlCrossRefPubMedWeb of Science
  59. 59.↵
    1. Habte M,
    2. Alexander M
    . 1977. Further evidence for the regulation of bacterial populations in soil by protozoa. Arch Microbiol 113:181–183. doi:10.1007/BF00492022.
    OpenUrlCrossRefPubMedWeb of Science
  60. 60.↵
    1. Brewer MT,
    2. Agbedanu PN,
    3. Zamanian M,
    4. Day TA,
    5. Carlson SA
    . 2013. Evidence for a bacterial lipopolysaccharide-recognizing G-protein-coupled receptor in the bacterial engulfment by Entamoeba histolytica. Eukaryot Cell 12:1433–1438. doi:10.1128/EC.00150-13.
    OpenUrlAbstract/FREE Full Text
  61. 61.↵
    1. Andersen LO,
    2. Bonde I,
    3. Nielsen HB,
    4. Stensvold CR
    . 2015. A retrospective metagenomics approach to studying Blastocystis. FEMS Microbiol Ecol 91:fiv072. doi:10.1093/femsec/fiv072.
    OpenUrlCrossRefPubMed
  62. 62.↵
    1. Paine CE,
    2. Beck H,
    3. Terborgh J
    . 2016. How mammalian predation contributes to tropical tree community structure. Ecology 97:3326–3336. doi:10.1002/ecy.1586.
    OpenUrlCrossRef
  63. 63.↵
    1. Morán XAG,
    2. Gasol JM,
    3. Pernice MC,
    4. Mangot JF,
    5. Massana R,
    6. Lara E,
    7. Vaqué D,
    8. Duarte CM
    . 2017. Temperature regulation of marine heterotrophic prokaryotes increases latitudinally as a breach between bottom-up and top-down controls. Glob Chang Biol 23:3956–3964. doi:10.1111/gcb.13730.
    OpenUrlCrossRef
  64. 64.↵
    1. Terborgh JW
    . 2015. Toward a trophic theory of species diversity. Proc Natl Acad Sci U S A 112:11415–11422. doi:10.1073/pnas.1501070112.
    OpenUrlAbstract/FREE Full Text
  65. 65.↵
    1. Fagarasan S,
    2. Muramatsu M,
    3. Suzuki K,
    4. Nagaoka H,
    5. Hiai H,
    6. Honjo T
    . 2002. Critical roles of activation-induced cytidine deaminase in the homeostasis of gut flora. Science 298:1424–1427. doi:10.1126/science.1077336.
    OpenUrlAbstract/FREE Full Text
  66. 66.↵
    1. Suzuki K,
    2. Meek B,
    3. Doi Y,
    4. Muramatsu M,
    5. Chiba T,
    6. Honjo T,
    7. Fagarasan S
    . 2004. Aberrant expansion of segmented filamentous bacteria in IgA-deficient gut. Proc Natl Acad Sci U S A 101:1981–1986. doi:10.1073/pnas.0307317101.
    OpenUrlAbstract/FREE Full Text
  67. 67.↵
    1. Moon C,
    2. Baldridge MT,
    3. Wallace MA,
    4. Burnham CAD,
    5. Virgin HW,
    6. Stappenbeck TS
    . 2015. Vertically transmitted faecal IgA levels determine extra-chromosomal phenotypic variation. Nature 521:90–93. doi:10.1038/nature14139.
    OpenUrlCrossRefPubMed
  68. 68.↵
    1. Bjarnason I
    . 2017. The use of fecal calprotectin in inflammatory bowel disease. Gastroenterol Hepatol 13:53–56. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5390326/.
    OpenUrl
  69. 69.↵
    1. Cheng HS,
    2. Guo YL,
    3. Shin JW
    . 2003. Hematological effects of Blastocystis hominis infection in male foreign workers in Taiwan. Parasitol Res 90:48–51. doi:10.1007/s00436-002-0804-3.
    OpenUrlCrossRefPubMed
  70. 70.↵
    1. Puthia MK,
    2. Vaithilingam A,
    3. Lu J,
    4. Tan KS
    . 2005. Degradation of human secretory immunoglobulin A by Blastocystis. Parasitol Res 97:386–389. doi:10.1007/s00436-005-1461-0.
    OpenUrlCrossRefPubMed
  71. 71.↵
    1. Koh A,
    2. De Vadder F,
    3. Kovatcheva-Datchary P,
    4. Bäckhed F
    . 2016. From dietary fiber to host physiology: short-chain fatty acids as key bacterial metabolites. Cell 165:1332–1345. doi:10.1016/j.cell.2016.05.041.
    OpenUrlCrossRefPubMed
  72. 72.↵
    1. Seth EC,
    2. Taga ME
    . 2014. Nutrient cross-feeding in the microbial world. Front Microbiol 5:350. doi:10.3389/fmicb.2014.00350.
    OpenUrlCrossRefPubMed
  73. 73.↵
    1. Drossman DA
    . 2006. The functional gastrointestinal disorders and the Rome III process. Gastroenterology 130:1377–1390. doi:10.1053/j.gastro.2006.03.008.
    OpenUrlCrossRefPubMedWeb of Science
  74. 74.↵
    1. Scicluna SM,
    2. Tawari B,
    3. Clark CG
    . 2006. DNA barcoding of Blastocystis. Protist 157:77–85. doi:10.1016/j.protis.2005.12.001.
    OpenUrlCrossRefPubMedWeb of Science
  75. 75.↵
    1. Clark CG
    . 1997. Extensive genetic diversity in Blastocystis hominis. Mol Biochem Parasitol 87:79–83. doi:10.1016/S0166-6851(97)00046-7.
    OpenUrlCrossRefPubMedWeb of Science
  76. 76.↵
    1. Verweij JJ,
    2. Blangé RA,
    3. Templeton K,
    4. Schinkel J,
    5. Brienen EA,
    6. van Rooyen MA,
    7. van Lieshout L,
    8. Polderman AM
    . 2004. Simultaneous detection of Entamoeba histolytica, Giardia lamblia, and Cryptosporidium parvum in fecal samples by using multiplex real-time PCR. J Clin Microbiol 42:1220–1223. doi:10.1128/JCM.42.3.1220-1223.2004.
    OpenUrlAbstract/FREE Full Text
  77. 77.↵
    1. Nazeer JT,
    2. El Sayed Khalifa K,
    3. von Thien H,
    4. El-Sibaei MM,
    5. Abdel-Hamid MY,
    6. Tawfik RA,
    7. Tannich E
    . 2013. Use of multiplex real-time PCR for detection of common diarrhea causing protozoan parasites in Egypt. Parasitol Res 112:595–601. doi:10.1007/s00436-012-3171-8.
    OpenUrlCrossRefPubMed
  78. 78.↵
    1. Bartram AK,
    2. Lynch MD,
    3. Stearns JC,
    4. Moreno-Hagelsieb G,
    5. Neufeld JD
    . 2011. Generation of multimillion-sequence 16S rRNA gene libraries from complex microbial communities by assembling paired-end Illumina reads. Appl Environ Microbiol 77:3846–3852. doi:10.1128/AEM.02772-10.
    OpenUrlAbstract/FREE Full Text
  79. 79.↵
    1. Comeau AM,
    2. Douglas GM,
    3. Langille MG
    . 2017. Microbiome helper: a custom and streamlined workflow for microbiome research. mSystems 2:e00127-16. doi:10.1128/mSystems.00127-16.
    OpenUrlAbstract/FREE Full Text
  80. 80.↵
    1. Schloss PD,
    2. Westcott SL,
    3. Ryabin T,
    4. Hall JR,
    5. Hartmann M,
    6. Hollister EB,
    7. Lesniewski RA,
    8. Oakley BB,
    9. Parks DH,
    10. Robinson CJ,
    11. Sahl JW,
    12. Stres B,
    13. Thallinger GG,
    14. Van Horn DJ,
    15. Weber CF
    . 2009. Introducing Mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol 75:7537–7541. doi:10.1128/AEM.01541-09.
    OpenUrlAbstract/FREE Full Text
  81. 81.↵
    1. Kuczynski J,
    2. Stombaugh J,
    3. Walters WA,
    4. González A,
    5. Caporaso JG,
    6. Knight R
    . 2011. Using QIIME to analyze 16S rRNA gene sequences from microbial communities. Curr Protoc Bioinformatics Chapter 10:Unit 10.7. doi:10.1002/0471250953.bi1007s36.
    OpenUrlCrossRef
  82. 82.↵
    1. Schloss PD,
    2. Gevers D,
    3. Westcott SL
    . 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.
    OpenUrlCrossRefPubMed
  83. 83.↵
    1. Eren AM,
    2. Morrison HG,
    3. Lescault PJ,
    4. Reveillaud J,
    5. Vineis JH,
    6. Sogin ML
    . 2015. Minimum entropy decomposition: unsupervised oligotyping for sensitive partitioning of high-throughput marker gene sequences. ISME J 9:968–979. doi:10.1038/ismej.2014.195.
    OpenUrlCrossRefPubMed
  84. 84.↵
    1. Eren AM,
    2. Maignien L,
    3. Sul WJ,
    4. Murphy LG,
    5. Grim SL,
    6. Morrison HG,
    7. Sogin ML
    . 2013. Oligotyping: differentiating between closely related microbial taxa using 16S rRNA gene data. Methods Ecol Evol 4:1111–1119. doi:10.1111/2041-210X.12114.
    OpenUrlCrossRef
  85. 85.↵
    1. DeSantis TZ,
    2. Hugenholtz P,
    3. Larsen N,
    4. Rojas M,
    5. Brodie EL,
    6. Keller K,
    7. Huber T,
    8. Dalevi D,
    9. Hu P,
    10. Andersen GL
    . 2006. Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl Environ Microbiol 72:5069–5072. doi:10.1128/AEM.03006-05.
    OpenUrlAbstract/FREE Full Text
  86. 86.↵
    1. Yilmaz P,
    2. Parfrey LW,
    3. Yarza P,
    4. Gerken J,
    5. Pruesse E,
    6. Quast C,
    7. Schweer T,
    8. Peplies J,
    9. Ludwig W,
    10. Glöckner FO
    . 2014. The SILVA and “All-species Living Tree Project (LTP)” taxonomic frameworks. Nucleic Acids Res 42:D643–D648. doi:10.1093/nar/gkt1209.
    OpenUrlCrossRefPubMedWeb of Science
  87. 87.↵
    1. Parks DH,
    2. Tyson GW,
    3. Hugenholtz P,
    4. Beiko RG
    . 2014. STAMP: statistical analysis of taxonomic and functional profiles. Bioinformatics 30:3123–3124. doi:10.1093/bioinformatics/btu494.
    OpenUrlCrossRefPubMed
  88. 88.↵
    1. Campbell JM,
    2. Fahey GC, Jr,
    3. Wolf BW
    . 1997. Selected indigestible oligosaccharides affect large bowel mass, cecal and fecal short-chain fatty acids, pH and microflora in rats. J Nutr 127:130–136. doi:10.1093/jn/127.1.130.
    OpenUrlAbstract/FREE Full Text
  89. 89.↵
    1. Witko-Sarsat V,
    2. Friedlander M,
    3. Nguyen Khoa T,
    4. Capeillère-Blandin C,
    5. Nguyen AT,
    6. Canteloup S,
    7. Dayer JM,
    8. Jungers P,
    9. Drüeke T,
    10. Descamps-Latscha B
    . 1998. Advanced oxidation protein products as novel mediators of inflammation and monocyte activation in chronic renal failure. J Immunol 161:2524–2532.
    OpenUrlAbstract/FREE Full Text
  90. 90.↵
    1. Das NP,
    2. Ratty AK
    . 1987. Studies on the effects of the narcotic alkaloids, cocaine, morphine, and codeine on nonenzymatic lipid peroxidation in rat brain mitochondria. Biochem Med Metab Biol 37:258–264. doi:10.1016/0885-4505(87)90035-1.
    OpenUrlCrossRefPubMedWeb of Science
  91. 91.↵
    1. McMurdie PJ,
    2. Holmes S
    . 2013. Phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS One 8:e61217. doi:10.1371/journal.pone.0061217.
    OpenUrlCrossRefPubMed
  92. 92.↵
    1. Neuwirth E
    . 2011. RColorBrewer: ColorBrewer palettes. https://CRAN.R-project.org/package=RColorBrewer.
  93. 93.↵
    1. Oksanen J,
    2. Blanchet FG,
    3. Kindt R,
    4. Legendre P,
    5. Minchin PR,
    6. O’Hara RB,
    7. Simpson GL,
    8. Solymos P,
    9. Henry M,
    10. Stevens H,
    11. Wagner H
    . 2013. vegan: community ecology package. https://CRAN.R-project.org/package=vegan.
  94. 94.↵
    1. Ploner A
    . 2012. Heatplus: heatmaps with row and/or column covariates and colored clusters. https://bioconductor.org/packages/release/bioc/html/Heatplus.html.
  95. 95.↵
    1. Pollard K,
    2. Gilbert HN,
    3. Ge Y,
    4. Taylor S,
    5. Dudoit S
    . 2005. Multiple testing procedures: R multtst package and applications to genomics, p 249–271. In Gentleman R, Carey V, Huber W, Irizarry R, Dudoit S (ed), Bioinformatics and Computational Biology Solutions Using R and Bioconductor. Springer-Verlag, New York, NY.
  96. 96.↵
    1. Wickham H
    . 2009. ggplot2: elegant graphics for data analysis. Springer-Verlag, New York, NY.
  97. 97.↵
    1. Wickham H
    . 2011. The split-apply-combine strategy for data analysis. J Stat Soft 40:1–29. doi:10.18637/jss.v040.i01.
    OpenUrlCrossRef
  98. 98.↵
    1. Love MI,
    2. Huber W,
    3. Anders S
    . 2014. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol 15:550. doi:10.1186/s13059-014-0550-8.
    OpenUrlCrossRefPubMed
Back to top
Download PDF
Citation Tools
Print
Alerts
Email

KEYWORDS

Blastocystis
gut microbiome
eukaryome
host-microbe interactions
microbial ecology

Related Articles

Cited By...