Skip to main content

Enteric parasitic infection disturbs bacterial structure in Mexican children with autoantibodies for type 1 diabetes and/or celiac disease



Intestinal bacterial dysbiosis and increased gut permeability are associated with higher risk of developing type 1 diabetes (T1D) or celiac disease (CD). There is a lack of information on parasitism involved in gut disturbance of predisposed children. We evaluated the effect of enteropathogenic parasites (Cryptosporidium spp., Cyclospora spp. G. lamblia, and Blastocystis spp.) on the bacterial structure of feces from children with autoantibodies for T1D or CD. Participants included 37 children under 18 years of age, from whom stools were analyzed for enteric parasites by qPCR and 22/37 for bacterial profile by sequencing the V3–V4 region of the 16s rRNA gene. Dietary, clinical, and socioeconomic data was recorded.


Pathogens parasitized 28/37 participants, Cryptosporidium spp. was the most prevalent (62.2%), followed by both Cyclospora cayetanensis and Blastocystis spp (37.8%). There were no dietary differences (p > 0.05) attributable to parasitism. Co-infected participants with Cryptosporidium and Cyclospora did not differ (p = 0.064) from non-infected participants in bacterial alpha phylogenetic diversity. The same parasites’ co-infection was associated with a decreased abundance of the Ruminococaceae (p = 0.04) and Verrucomicrobioceae families, of the Akkermansia genus (p = 0.009). There was a lower Firmicutes/Bacteroidetes ratio (p = 0.02) in infected than in uninfected participants.


Cryptosporidium and Cyclospora affected the bacterial structure at family and genus levels, decreasing the ratio between Firmicutes and Bacteroidetes in children with auto-antibodies for T1D or CD, which could increase the risk of illness onset.


Twenty-five years ago, studies about the pathogenesis of celiac disease (CD), an autoimmune enteropathy triggered by wheat consumption, proposed that gut health disturbance influences its onset [30]. Nowadays, the disruption of tight junctions by bacterial-derived products and/or gliadin peptides, which affects the intestinal permeability, is related to the development of CD [11, 23]. Gut health has also been associated to the onset of type 1 diabetes (T1D), through the identification of three gut dynamic components involved in its pathogenesis: the mucosal immunity, the microbiota and the individual’s diet [33, 39]. T1D is another autoimmune disease similar to CD in predisposition factors, characterized by the destruction of the insulin-producing β-cells in pancreas. It has been related to a deficiency of short-chain fatty acids-producing bacteria in the gut, as well as to an increased intestinal permeability [20, 36].

CD and T1D are the two most frequent autoimmune diseases in childhood. The International Diabetes Federation estimates that, globally, more than one million children and adolescents under the age of 19 have T1D, with an annual incidence of almost one hundred thirty thousand new cases [25]. Both T1D and CD incidence have increased in recent decades, mainly attributed to common environmental or pathophysiological mechanisms. Thus, presenting T1D increases the risk of developing CD, in a range between 2.5 and 16.4% in different populations, several times its global prevalence of 2–3% [13].

The intestinal bacterial imbalance (dysbiosis) and the enteropathogenic parasites that affect the gut barrier, have been independently associated with increased risk of developing T1D or CD in genetically predisposed children [9, 23, 35]. However, there is no information regarding the effect of enteropathogenic parasites’ infections on bacterial dysbiosis or vice versa, in children developing autoimmunity. Therefore, our study had two objectives. The first was to evaluate the prevalence of enteropathogenic parasites in children under 18 years with autoantibodies associated with T1D or CD, pre-symptomatic or already diagnosed. The second one was to look for the effect of the enteropathogenic parasites’ infection in the fecal bacterial structure in children developing autoimmunity associated to T1D and/or CD.


Characteristics of participants

Participants were 37 children, 25 of them were girls. Their mean age was 10.93 ± 3.02 years, weight 39.52 ± 17.2 kg and height 141.44 ± 16.05 cm (values are means ± standard deviations). Socioeconomic status was medium for 26 of the participants, with sufficient resources to meet all the family requirements; but 11 belonged to families with poor household circumstances and had access to only some public services.

Autoantibodies and diagnosis status for T1D and CD

All participants presented antibodies associated with T1D or CD because it was an inclusion criterion. Sixteen of 37 participants had a diagnosis for T1D and 2 for CD; 4 of them had autoimmunity (at least two different autoantibodies) and 15 presented only one type of antibodies related to T1D. Ten of the 18 diagnosed participants were at onset of the illness, while one of them had 12 months of evolution and the remaining 7 were diagnosed with T1D at least for 24 months. All of participants had a high genetic risk for T1D or CD, presenting HLA DQ2, DQ8 haplotypes, or any combination of related alleles in their genotype.

Detection and identification of enteropathogenic parasites

Twenty-eight of 37 participants were infected by enteropathogenic parasites. Cryptosporidium spp. showed the highest prevalence (62.2%), followed by Cyclospora spp. and Blastocystis spp. (37.8%). Table 1 shows prevalence of all detected intestinal parasites.

Table 1 Enteric parasites detection in participants with T1D or CD associated antibodies (n = 37)

Dietary evaluation

There were no significant differences (p > 0.05) for energy, macronutrient nor dietary fiber consumption between infected and uninfected participants with enteropathogenic parasites. The energy intake (mean ± standard deviation) of the group with enteric pathogen parasites (n = 28) was 1837.6 ± 919 kcal/day, with 33.8% lipids, 14.6% proteins and 53.2% carbohydrates. In the uninfected group (n = 9), energy intake was 1667.7 ± 489 kcal/day, with 37.2% lipids, 15.3% protein and 48.7% carbohydrates. Almost half of the total children (49%) consumed a healthy dietary pattern, while 35% of them had an ultra-processed pattern, and 16% received a poor dietary pattern.

Fecal bacterial microbiota analysis

Fecal bacterial microbiota sequencing was done in samples of 22 participants, 18 infected with enteropathogenic parasites and 4 uninfected, whose characteristics are shown in Table 2. Infected participants had at least one parasite; 15 with Cryptosporidium spp. 9 with Cyclospora spp. 6 with G. lamblia and 12 with Blastocystis spp. Only 3 of the 22 analyzed samples were from 2 children at the onset of CD, and the third one coursed with more than 2 years of T1D evolution. Seven children presented more than one autoantibody’s type, and the remaining 15 children were only positive for insulin autoantibodies (IAA).

Table 2 General characteristics of the participants in the fecal microbiota analysis (n = 22)

After using quality filters to remove sequences containing ambiguous bases, barcode mismatches, or low quality reads (Phred quality scores < 25), a total 636,876 16S rRNA reads were obtained, with an average of 28,949 sequences per sample. After removing the singletons, there were still 25,328 sequences per sample. The minimum number of sequences reached for any of the samples simultaneously was 16,000; therefore, this value was used as the number of sequences to be generated in the rarefaction process. In this way, a standardized BIOM file was obtained for all samples. The depth of coverage in terms of the number of bacterial taxa per participant was adequate, since the tendency to asymptote was achieved in each curve. No significant differences were found in the rarefaction curves of phylogenetic diversity (p = 0.88) nor in the principal coordinate analysis (PCoA) of fecal communities (p = 0.139) when unweighted Unifrac distances were estimated, between the infected and uninfected groups.

Otherwise, when plotting the PCoA of fecal bacterial communities in participants with diagnosis of T1D or CD (n = 3) against those with only one auto-antibody type (n = 15), the already diagnosed participants had a lower bacterial diversity within group (p = 0.01) than between groups. In the same analysis, autoimmune participants (with more than two positive autoantibodies for T1D), had a similar distribution of their UniFrac distances (p > 0.05) to those who were diagnosed with T1D or CD, as shown in Fig. 1a. However, no differences were found in alpha diversity (p > 0.05), evaluated as rarefaction curves of phylogenetic diversity among the three groups.

Fig. 1

Fecal microbiota structure in Mexican children with positive T1D- or CD-associated antibodies (n = 22). a Weighted UniFrac distances between groups. b Bacterial relative abundance at genus level per participant. c Bacterial relative abundance according to their infection status by enteropathogenic parasites at phylum level and d at genus level (all phyla and genus with ≥ 1% of relative abundance were considered). e Linear Discriminant Analysis Effect Size (LEfSe) analysis is showing those OTUs that were significantly differentially abundant between autoimmune children infected or non-infected with Cyptosporidium and/or Cyclospora, ranked by effect size (all LDA scores > 2). Only IAA: Positives only for T1D-associated anti-insulin autoantibodies, Two Ab: Positives for 2 T1D and CD-associated autoantibodies, T1D: Type 1 diabetes, CD: Celiac disease, n: non-infected, i: infected with Cyptosporidium and/or Cyclospora, NS: non-significant (Student’s t test), p values ≤ 0.05 were considered significant

Figure 1b shows the relative bacterial genus abundance in the fecal microbiota of the 22 evaluated participants. Prevotella was the dominant enterotype (relative abundance of 30–45%) in samples from 12 participants; 9 with only AAI autoantibodies, 1 with autoimmunity but not symptomatic yet and 2 with anti-TG. Seven samples presented more than 20% of Bacteroides abundance, 2 of them from participants at CD onset where one was infected with three enteropathogenic parasites, and the other was uninfected. The rest of the sequenced samples corresponded to 5 participants with only IAA antibodies. Independently of enteropathogenic parasites, the group of participants with T1D or CD (samples 20i, 21i and 22n), presented an increased (p = 0.0001) relative abundance of Ruminococcus (5.99%) as compared to the abundance of the only IAA (1.41%) or 2 Ab (0.65%) groups. In addition, its Bacteroides relative abundance is also significantly higher (p = 0.001) than it of the 2 Ab group (Fig. 1b).

When assessing the individual effect of enteropathogenic parasites infection on bacterial diversity, no differences were detected in the rarefaction curves. However, participants with Cryptosporidium spp. tended (p = 0.057) to have lower bacterial diversity than the uninfected ones, while the p value was higher (p = 0.07) in those participants with Cyclospora spp., as well as for co-infected participants with both parasites (p = 0.064). Additionally, it highlights that parasitized participants with Cryptosporidium spp. and/or Cyclospora spp. had a lower (p > 0.0001) intragroup variability than non-infected children, being its most representative genera Roseburia, Lactobacillus and Haemophilus, as shown in Linear Discriminant Effect Size Analysis (LEfSe), schematized in Fig. 1e.

At the Phyla level, Cryptosporidium spp. and/or Cyclospora spp. infected participants had a lower abundance of Verrucomicrobia (p = 0.013) in feces than the uninfected ones (0.15 vs. 1.29%), as shown in Fig. 1c. This was explained by an increased abundance of the genus Akkermansia (p = 0.009) in uninfected participants when compared to the infected ones (1.21% ± 1.52 vs. 0.1% ± 0.21), as shown in Fig. 1d. In the same way, the Ruminococcaceae bacterial family from phylum Firmicutes, was more abundant (p = 0.04) in participants without Cryptosporidium spp. and/or Cyclospora spp. infections in comparison to the infected ones (16.36% ± 6.38 vs 11.36% ± 4.3). It was the operational taxonomic unit (OTU) with a higher linear discriminant analysis (LDA) score in the LEfSe analysis for that group (Fig. 1c, d). The Firmicutes/Bacteroidetes ratio was higher (p = 0.02) in uninfected participants in contrast with those with Cryptosporidium spp. and/or Cyclospora spp. (0.961 ± 0.142 vs. 0.653 ± 0.048). Figure 2 shows the visualization of the most abundant bacterial genera on a heatmap.

Fig. 2

Heatmap analysis of gut microbiota. a Bacterial relative abundance means, in infected or non-infected children and b per participant. All OTUs with relative abundance greater than 0.01% were considered. i infected, n non infected


Damage of the gut mucosa epithelial barrier by parasites or other factors that affect permeability has been associated with the onset of autoimmune diseases, such as CD and T1D. Also, changes in the gut bacterial composition related to short-chain fatty acid production are involved in the same problem [11, 23, 36]. Therefore, we hypothesize that enteropathogenic parasites adhered or internalized in epithelial cells in the lining of the small intestine of high-risk children, who are developing autoimmunity, can induce bacterial imbalance and trigger the disease onset. This fact could be significant in children populations with a high prevalence of enteropathogenic parasites, such as the investigated Mexican population.

Firstly, our participants with high risk or already suffering from T1D or CD were highly parasitized. In our sampled children, infections with Cryptosporidium spp were present more than twice than the ones previously found in children from an open population in the same area [27]. Although cryptosporidiosis is common in our area, C. parvum and C. hominis infect children principally under five years old [34], and most of our participants were older than ten years. Blastocystis spp and C. cayetanenesis were also higher, and G. lamblia was similar to the reported data. The increased parasitism may be due to the detection method used (qPCR), which is more sensible than the conventional microscopy used in previous studies. However, parasites can opportunistically infect children with altered immunity [6, 38].

Thus, in a subsample of participants, the first part of our hypothesis was contrasted, although it is not possible to separate the effect of autoimmunity itself from the effect of parasitism. We found that the fecal bacterial microbiota structure of studied participants was significantly affected by the autoimmunity status (p = 0.012), in agreement with our previous study [21], as well as, because of the infection by Cryptosporidium spp. alone (p = 0.027) or in co-infection with Cyclospora spp. (p = 0.035). There were no confusion factors, as first year of life and current diet, genetic predisposition, or antibiotic use, involved for such effects to explain the microbiota differences between our groups (p > 0.05).

The effect of infection by Cryptosporidium spp. and Cyclospora spp. on the abundance and structure of some intestinal bacteria is poorly studied. Their previous life forms before going inside of the epithelial cells may interact with the intestinal lumen components, affecting the bacterial microbiota composition [28]. Alternatively, β-defensins 1 and 2, produced by the immune system against the parasites’ infection, may also induce changes in the bacterial composition [16].

The prevalence of Blastocystis spp. infection in our population appears high, although it is unclear if it is a pathogen, a commensal or an opportunistic parasite [6], because it can induce intestinal illness, but it can also produce an asymptomatic or self-limited course in the host. Possibly that was the reason why Audebert et al. [5] found effects of Blastocystis spp. infection on α and β bacterial diversity in human beings, whereas Nagel et al. [22] associated the bacterial changes to the irritable bowel syndrome of their patients. In agreement with Andersen et al. [4], we found that Blastocystis spp. in some fecal samples of infected participants were concomitant with a high relative abundance of Ruminococcus and Prevotella.

Giardia lamblia adheres to the epithelial cells in the small intestine, damaging the intestinal mucosa in a similar spectrum to the one produced by CD [26, 32]. Therefore, giardiasis induces transitory intolerance to gluten proteins, which can evolve to CD in genetically predisposed individuals [10, 24]. There is evidence of the onset of CD after giardiasis [10], and eventually, this infection could increase even more the risk of CD in children with T1D.

According to Maffeis et al. [20], both, the altered intestinal permeability and dysbiosis, are contemporaneously associated with the pre-pathogenic period of T1D condition. From our results, we could infer that infections of intracellular Cryptosporidium spp. and/or Cyclospora spp. could induce both increased gut permeability and dysbiosis, at least in our participants. We knew from our previous study [21] that dysbiosis was very strong in children at onset of T1D, with an increase in Bacteroides abundance. However, we had no information of the previous state (presence of autoantibodies or autoimmunity) in relation to dysbiosis, or about the effect of protozoan parasitic infection of our participants.

The microbiota structure of participants with just AAI antibodies, was similar to that of healthy participants, as well as to the ones of more than 2 years with T1D in our previous study, respecting Prevotella and Bacteroides proportions [21]. An interesting detail is that asymptomatic participants with two antibody types, presented about 6% Elusimicrobia and > 10% Cianobacteria, genera that was not found before in previous analyzed samples from the same population. Besides an increased abundance of Bacteroides, we also found an increased abundance of Ruminococcus (5.99%) in symptomatic participants with T1D or CD. Although data of bacterial species were not accessible in this study,, it could be due to Bacteroides dorei and Ruminococcus gravus found out by Abdellatif et al. [1], prior to the autoimmunity development.

It is remarkable that participants with no Cryptosporidium spp. and/or Cyclospora spp. infections had higher relative abundance of Akkermansia than infected participants. Hu et al. [16] demonstrated that Akkermansia promotes mucin production to protect the intestinal epithelium of the enteropathogenic parasites’ infection, and increases an antimicrobial peptide, Reg3γ [14]. Akkermansia is related also to the delay of T1D onset in a murine model [15]. Possibly, our uninfected children were protected by Akkermansia because they were mainly participants with only AAI positive.

The higher Firmicutes/Bacteroidetes ratio of the uninfected participants from our study, when compared to those infected with Cryptosporidium spp. and/or Cyclospora spp. could be explained, in part, by a higher amount of the Firmicutes’ Ruminococcaceae family in the uninfected participants. However, at least in a murine model, cryptosporidiosis was associated with an increased relative abundance of Bacteroidetes such as the Prevotellaceae and Porhyromondaceae families [28].

Only three of the infected participants with Cryptosporidium spp. and/or Cyclospora spp., referred gastrointestinal symptoms even though two of them were positive for CD-associated autoantibodies. Perhaps it was due to the fact that Lactobacillaceae family and Lactobacillus genera, distinctive of the infected participants (Fig. 1e), protected them from a severe affection. According to Sanad et al. [29], in a murine model, probiotics such as Lactobacillus and Bifidobacterium decreased the number of oocysts of Cryptosporidium spp. in feces, resulting in a reduced damage. However, induced cryptosporidiosis in other murine model [28] did not affect the same Phyla as in our Cryptosporidium-infected children, as expected due to the physiological differences.


In conclusion, the prevalence of enteropathogenic parasites in northwestern Mexican children developing autoimmunity associated to T1D or CD is high; and particularly, Cryptosporidium spp. and/or Cyclospora spp. infection affected their fecal bacterial structure. The Firmicutes/Bacteroidetes ratio was modified by enteropathogenic parasites’ infection since the Verrucomicrobiaceae and Ruminococcaceae families were increased in uninfected children.


Characteristics of the Participants

The convenience non-probabilistic sampling study was carried out collecting data from 2016 to 2018 in northwestern Mexico (Sonora State), a desertic area with dry climate, previously described [2]. Children under 18 years old with a diagnosis of T1D or CD were recruited at the external consultation of the Children Sonora State Hospital and the Mexican Social Security Institute (HIES and IMSS). Additionally, children (7–12 years old) with autoantibodies related to T1D or CD were enrolled from public schools in the same geographical areas. Autoantibodies detection was done after an invitation to the schools’ principals, teachers and parents was made. A screening protocol was applied beginning with genotyping of HLA-DQ2 and DQ8 haplotypes, followed by familiar clinical history of risk, and ending with the autoantibodies’ detection.

After the parents’ consent and confirmation of no antibiotic or anti-parasites use in the last month, all participants were asked for one blood sample and feces samples of 3 consecutive days. The study was approved by the Ethical Committee of the Centro de Investigacion en Alimentacion y Desarrollo, A.C. (CE/016/2014) and an informed consent was signed by the parents. Positive and negative infected groups of participants with enteropathogenic parasites were formed for analysis.

Dried blood spots were used for gDNA extraction and HLA-DQ2 and DQ8 haplotypes analyses according to our technique [3]. After blood serum separation, antibodies analyses were performed. T1D associated autoantibodies, such as anti-insulin (AAI), anti-glutamate decarboxylase (GAD) and anti-tyrosin pancreatic phosphatase (IA-2) were analyzed by ELISA through commercial kits (MyBiosource, Cat. No. MBS772002 and Kronus, Cat. No. KR7770) or related to CD, such as anti-transglutaminase (TG) and anti-gliadins, by the ELISA method of Cabrera-Chávez et al. [7].

Clinical history, diet, anthropometry, and socioeconomic status

Clinical history and socioeconomic information were obtained from hospital data and mother interviews. Three 24-h non consecutive recalls collected dietary intake data. Weight and height were measured.

Enteropathogenic parasites identification

The three feces samples were homogenized together and aliquots were stored at − 70 °C for posterior analyses. Genomic DNA (gDNA) was extracted from fresh feces samples using the QIAmp Fast DNA stool mini kit (Qiagen, Cat. No. 51604), according to the manufacturer’s instructions. DNA was eluted in a final volume of 200 μL and stored at 20 °C. Concentration and quality were measured in a Nanodrop 2000 (Thermo Scientific, Pittsburgh, USA).

Feces’ gDNA samples were analyzed for enteropathogenic parasites Cryptosporidium spp., Blastocystis spp., Cyclospora spp., Entamoeba histolytica (E. histolytica) and Giardia lamblia (G. lamblia) by real-time PCR using specific primers for fragments in the 16S-like rRNA gene. Reactions were done with 5 μM of each primer, 10 uL of SYBR Green Supermix (Bio-Rad, Cat. No. 1708882), 100 ng of template DNA, and milli-Q water up to 20 μL final volume. Reactions were run in a StepOnePlus real-time PCR system (Applied Biosystems, Foster City, USA), to obtain the melting temperature (Tm) of each amplicon. Cycling conditions were different for each pathogen. For Cyclospora spp., conditions began with an initial hold at 95 °C for 10 min, followed by 40 cycles at 95 °C for 15 s and 60 °C for 1 min [37]. For Cryptosporidium spp., conditions were an initial 95 °C for 3 min, followed by 50 cycles at 94 °C for 15 s, 54 °C for 30 s and 72 °C for 30 s [18]. The Blastocystis spp. identification was done according to Grabensteiner and Hess [12], starting with an initial hold at 95 °C for 15 min, followed by 40 cycles at 94 °C for 30 s, 60 °C for 30 s and 72 °C for 60 s. For the E. histolytica, the conditions were the same previously described by ourselves [2]. For G. lamblia, the conditions were an initial 94 °C for 5 min, followed by 35 cycles at 94 °C for 30 s, 58 °C for 90 s, 72 °C for 90 s and 72 °C for 10 min [40].

Bacterial microbiota analysis of feces

The extracted fecal gDNA of a subsample of 22 participants was used for the bacterial microbiota analysis. For this evaluation, the V3-V4 region of the bacterial 16S rRNA was amplified by PCR using the primer pair (S-D-Bact-0341-b-S-17/S-D-Bact-0785-a-A-21) [17], using the Illumina Miseq system (San Diego, CA., USA).

Bioinformatics and statistical comparisons

Descriptive statistics were used for the characteristics of the participants. Student’s t test evaluated the diet composition between the infected and uninfected groups. Analyses were run in NCSS Statistical Software version 2007 considering significant p values ≤ 0.05.

The sequence analysis of the diversity of the microbial communities was made with the QIIME software package (Quantitative Insights Into Microbial Ecology), version 1.8.0 [8]. Quality filters were used to remove sequences containing ambiguous bases, barcode mismatches, or low quality reads (Phred quality scores < 25). After trimming barcodes, the sequences were clustered using UCLUST into operational taxonomic units (OTUs), based on 97% sequence similarity against the Green Genes reference sequence collection (version 13_5). The option of reverse strand matching was enabled. We selected the closed-reference OTU picking method that was a reference-based approach. Thus, chimera removal was not necessary. All analyses were restricted to OTUs comprising ≥ 0.1% of the reads. Weighted UniFrac distances and Principal Coordinate Analysis were used for the analysis of bacterial communities and comparisons [19]. Relative abundances for phylum, class, order and family were obtained and plotted in Excel. Genus whose relative abundance was greater than 1% were represented on a heat map; the visualization was carried out with the Morpheus software ( A linear discriminant analysis was performed to measure the effect (LEfSe) of the statistical significance of the obtained results with the biological consistency and to estimate the size of the effect [31]. Differences between groups were evaluated by ANOVA and Tukey–Kramer multiple comparison tests.

Availability of data and materials

The datasets analyzed during the current study are available from the corresponding author on reasonable request.


  1. 1.

    Abdellatif AM, Jensen Smith H, Harmas RZ, Sarvetnick NE. Human islet response to selected type 1 diabetes-associated bacteria: a transcriptome-based study. Front Immunol. 2019.

    Article  PubMed  PubMed Central  Google Scholar 

  2. 2.

    Aguayo-Patrón S, Castillo-Fimbres R, Quihui-Cota L, Calderón de la Barca AM. Use of real-time polymerase chain reaction to identify Entamoeba histolytica in schoolchildren from northwest Mexico. J Infect Dev Ctries. 2017.

    Article  PubMed  Google Scholar 

  3. 3.

    Aguayo-Patrón S, Beltrán-Sauceda L, Calderón de la Barca AM. A population-wide applicable HLA-DQ2 and DQ8 genotyping using DNA from dried blood spots and duplex allele-specific qPCR amplification. Scand J Clin Lab Invest. 2016;76:581–7.

    CAS  Article  PubMed  Google Scholar 

  4. 4.

    Andersen LO, Bonde I, Nielsen HB, Stensvold CR. A retrospective metagenomics approach to studying Blastocystis. FEMS Microbiol Ecol. 2015.

    Article  PubMed  Google Scholar 

  5. 5.

    Audebert C, Even G, Cian A, The Blastocystis Investigation Group, Loywick A, Merlin S, et al. Colonization with the enteric protozoa Blastocystis is associated with increased diversity of human gut bacterial microbiota. Sci Rep. 2016;1:1.

    CAS  Article  Google Scholar 

  6. 6.

    Bednarska M, Jankowska I, Pawelas A, Piwczynska K, Bajer A, Wolska-Kusnierz B, et al. Prevalence of Cryptosporidium, Blastocystis, and other opportunistic infections in patients with primary and acquired immunodeficiency. Parasitol Res. 2018.

    Article  PubMed  PubMed Central  Google Scholar 

  7. 7.

    Cabrera-Chávez F, Rouzaud-Sández O, Sotelo-Cruz N, Calderón de la Barca AM. Bovine milk caseins and transglutaminase-treated cereal prolamins are differentially recognized by IgA of celiac disease patients according to their age. J Agric Food Chem. 2009.

    Article  PubMed  Google Scholar 

  8. 8.

    Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010.

    Article  PubMed  PubMed Central  Google Scholar 

  9. 9.

    Cinek O, Kramna L, Lin J, Oikarinen S, Kolarova K, Ilonen J, et al. Imbalance of bacteriome profiles within the Finnish Diabetes Prediction and Prevention study: parallel use of 16S profiling and virome sequencing in stool samples from children with islet autoimmunity and matched controls. Pediatr Diabetes. 2017.

    Article  PubMed  Google Scholar 

  10. 10.

    Edling L, Rathsman S, Eriksson S, Bohr J. Celiac disease and giardiasis: a case report. Eur J Gastroenterol Hepatol. 2012.

    Article  PubMed  Google Scholar 

  11. 11.

    Freire R, Ingano L, Serena G, Cetinbas M, Anselmo A, Sapon A, et al. Human gut derived-organoids provide model to study gluten response and effects of microbiota derived molecules in celiac disease. Sci Rep. 2019.

    Article  PubMed  PubMed Central  Google Scholar 

  12. 12.

    Grabensteiner E, Hess M. PCR for the identification and differentiation of Histomonas meleagridis, Tetratrichomonas gallinarum and Blastocystis spp. Vet Parasitol. 2006.

    Article  PubMed  Google Scholar 

  13. 13.

    Hagopian W, Lee HS, Liu E, Rewers M, She J-X, Ziegler A-G, et al. Co-occurrence of Type 1 Diabetes and Celiac Disease Autoimmunity. Pediatrics. 2017.

    Article  PubMed  PubMed Central  Google Scholar 

  14. 14.

    Hänninen A, Toivonen R, Pöysti S, Belzer C, Plovier H, Ouwerkerk JP, et al. Akkermansia muciniphila induces gut microbiota remodeling and controls islet autoimmunity in NOD mice. Gut. 2018.

    Article  PubMed  Google Scholar 

  15. 15.

    Hansen CH, Krych L, Nielsen DS, Vogensen FK, Hansen LH, Sorensen SJ, et al. Early life treatment with vancomycin propagates Akkermansia muciniphila and reduces diabetes incidence in the NOD mouse. Diabetologia. 2012.

    Article  PubMed  PubMed Central  Google Scholar 

  16. 16.

    Hu G, Feng Y, O’Hara SP, Chen XM. Immunology of cryptosporidiosis. In: Cacciò SM, Widmer G, editors. Cryptosporidium: Parasite and Disease. Switzerland: Springer; 2014. p. 423–54.

    Google Scholar 

  17. 17.

    Klindworth A, Pruesse E, Schweer T, Peplies J, Quast C, Horn M, et al. Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies. Nucleic Acids Res. 2013.

    Article  PubMed  Google Scholar 

  18. 18.

    Limor JR, La AA, Xiao L. Detection and differentiation of Cryptosporidium parasites that are pathogenic for humans by real-time PCR. J Clin Microbiol. 2002.

    Article  PubMed  PubMed Central  Google Scholar 

  19. 19.

    Lozupone C, Knight R. UniFrac: a new phylogenetic method for comparing microbial communities. Appl Environ Microbiol. 2005.

    Article  PubMed  PubMed Central  Google Scholar 

  20. 20.

    Maffeis C, Martina A, Corradi M, Quarella S, Nori N, Torriani S, et al. Association between intestinal permeability and faecal microbiota composition in Italian children with beta cell autoimmunity at risk for type 1 diabetes. Diabetes Metab. Res Rev. 2016.

    Article  PubMed  Google Scholar 

  21. 21.

    Mejía-León ME, Petrosino JF, Ajami NJ, Domínguez-Bello MG, Calderón de la Barca AM. Fecal microbiota imbalance in Mexican children with type 1 diabetes. Sci Rep. 2014.

    Article  PubMed  PubMed Central  Google Scholar 

  22. 22.

    Nagel R, Traub RJ, Allcock RJ, Kwan MM, Bielefeldt-Ohmann H. Comparison of faecal microbiota in Blastocystis-positive and Blastocystis-negative irritable bowel syndrome patients. Microbiome. 2016.

    Article  PubMed  PubMed Central  Google Scholar 

  23. 23.

    Olivares M, Benítez-Páez A, de Palma G, Capilla A, Nova E, Castillejo G, et al. Increased prevalence of pathogenic bacteria in the gut microbiota of infants at risk of developing celiac disease: the PROFICEL study. Gut Microbes. 2018.

    Article  PubMed  PubMed Central  Google Scholar 

  24. 24.

    Pal N, Sharma R, Sharma B, Suman R. A case of Cryptosporidium infection in a child of celiac disease. J Postgrad Med. 2012.

    Article  PubMed  Google Scholar 

  25. 25.

    Patterson CC, Karuranga S, Salpea P, Saeedi P, Dahlquist G, Soltesz G, et al. Worldwide estimates of incidence, prevalence and mortality of type 1 diabetes in children and adolescents: results from the International Diabetes Federation Diabetes Atlas. Diabetes Res Clin Pract. 2019;157:107842.

    Article  Google Scholar 

  26. 26.

    Pehlivanoğlu B, Doğanavşargil B, Sezak M, Nalbantoğlu İ, Korkmaz M. Gastrointestinal parasitosis: histopathological insights to rare but intriguing lesions of the gastrointestinal tract. Turk Patoloji Derg. 2016.

    Article  PubMed  Google Scholar 

  27. 27.

    Quihui-Cota L, Lugo-Flores CM, Ponce-Martínez JA, Morales-Figueroa GG. Cryptosporidiosis: a neglected infection and its association with nutritional status in schoolchildren in northwestern Mexico. J Infect Dev Ctries. 2015.

    Article  PubMed  Google Scholar 

  28. 28.

    Ras R, Huynh K, Desoky E, Badawy A, Widmer G. Perturbation of the intestinal microbiota of mice infected with Cryptosporidium parvum. Int J Parasitol. 2015.

    Article  PubMed  Google Scholar 

  29. 29.

    Sanad MM, Al-Malki JS, Al-Ghabban AG. Control of Cryptosporidiosis by probiotic bacteria. AEMS. 2015.

    Article  Google Scholar 

  30. 30.

    Schulzke JD, Schulzke I, Fromm M, Riecken EO. Epithelial barrier and ion transport in coeliac sprue: electrical measurements on intestinal aspiration biopsy specimens. Gut. 1995;37(6):777–82.

    CAS  Article  Google Scholar 

  31. 31.

    Segata N, Izard J, Waldron L, Gevers D, Miropolsky L, Garret WS, et al. Metagenomic biomarker discovery and explanation. Genome Biol. 2011.

    Article  PubMed  PubMed Central  Google Scholar 

  32. 32.

    Tchidjou HK, De Matteis A, Di Iorio L, Finocchi A. Celiac disease in an adoptive child with recurrent Giardia infection. Int J Health Sci. 2015;9:193–7.

    Google Scholar 

  33. 33.

    Tillett BJ, Hamilton-Williams EE. Microbiota derived factors as drivers of type 1 diabetes. Prog Mol Biol Transl Sci. 2020;171:215–35.

    Article  Google Scholar 

  34. 34.

    Urrea-Quezada A, González-Díaz M, Villegas-Gómez I, Durazo M, Hernandez J, Xiao L, et al. Clinical manifestations of cryptosporidiosis and identification of a new Cryptosporidium subtype in patients from Sonora. Mexico. Pediatr Infect Dis J. 2018.

    Article  PubMed  Google Scholar 

  35. 35.

    Vaarala O. Is the origin of type 1 diabetes in the gut? Immunol Cell Biol. 2012.

    Article  PubMed  Google Scholar 

  36. 36.

    Vatanen T, Franzosa EA, Schwager R, Tripathi S, Arthur TD, Vehik K, et al. The human gut microbiome in early-onset type 1 diabetes from the TEDDY study. Nature. 2018.

    Article  PubMed  PubMed Central  Google Scholar 

  37. 37.

    Verweij JJ, Laeijendecker D, Brienen EA, van Lieshout L, Polderman AM. Detection of Cyclospora cayetanensis in travelers returning from the tropics and subtropics using microscopy and real-time PCR. Int J Med Microbiol. 2003.

    Article  PubMed  Google Scholar 

  38. 38.

    Villa-Jiménez OM, Escobedo AE, Hano OM, Wood LW, Pérez FP, González LG. Blastocystis hominis en pacientes celiacos sintomáticos. Acta Gastroenterol Latinoam. 2012;42:175–81.

    Google Scholar 

  39. 39.

    Yap YA, Mariño E. An Insight Into the Intestinal Web of Mucosal Immunity, Microbiota, and Diet in Inflammation. Front Immunol. 2018;9:2617.

    Article  Google Scholar 

  40. 40.

    Zebardasta N, Yeganehb F, Gharavic MJ, Abadid A, Seyyed Tabaeie SJ, Haghighie AK. Simultaneous detection and differentiation of Entamoeba histolytica, E. dispar, E. moshkovskii, Giardia lamblia and Cryptosporidium spp. in human fecal samples using multiplex PCR and qPCR-MCA. Acta Trop. 2016;1:1.

    CAS  Article  Google Scholar 

Download references


Authors are grateful to pediatric endocrinologists Dr. G. Garcia-Galván and Dr. C. Bueno-Castro from IMSS and HIES hospitals and to the participant children and their parents for their good collaboration. We acknowledge also to Ph.D. F. Cornejo, Biol. F. Sánchez Lopez, and MSci. J. Peñuñuri for technical assistance, as well as to MSci. A.L. Cota Ojeda for English edition.


This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors. However, two of the coauthors received scholarship from CONACYT (Mexican Council for Science and Technology); the rest of coauthors’ salary is from the same Council, Universidad Autonoma de Baja California Mexicali or from Universidad Nacional Autonoma de Mexico (UNAM).

Author information




AMC designed, followed and wrote the first draft and edited the final manuscript; RCF was doing all the laboratory analysis and helped write the first draft; EML made the bioinformatics analysis and its description in methods, results and discussion of the manuscript; LQC did the previous analyses of enteric pathogenic parasites, by traditional techniques and edited the last version of the manuscript; AOL was in charge of the sequencing work and SAP made the first tests by real time PCR of the enteric parasites and followed the experiments and description at the manuscript about them. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Ana M. Calderón de la Barca.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the Ethical Committee of the Centro de Investigacion en Alimentacion y Desarrollo, A.C. (CE/016/2014) and an informed consent was signed by the parents of the participant children.

Consent for publication

Not applicable.

Competing interests

Authors declare no conflict of interests.

Additional information

Publisher's Note

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

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Calderón de la Barca, A.M., Castillo-Fimbres, R.S., Mejía-León, M.E. et al. Enteric parasitic infection disturbs bacterial structure in Mexican children with autoantibodies for type 1 diabetes and/or celiac disease. Gut Pathog 12, 37 (2020).

Download citation


  • Enteric pathogens
  • Cryptosporidium
  • Cyclospora
  • Akkermansia
  • Celiac disease
  • Type 1 diabetes