Characterization of Fecal Microbiota, Short-Chain Fatty Acids and Lactic Acid Concentrations in 5 - 8-Year-Old Children with Cow Milk Protein Allergy

authors:

avatar Yanru Dong 1 , avatar Peng Fei 2 , avatar Ying Han 3 , avatar Ling Guo 1 , *

Key Laboratory of Dairy Science, Ministry of Education, College of Food Science, Northeast Agricultural University, Harbin, China
College of Food and Biological Engineering, Henan University of Science and Technology, Luoyang, China
Department of Hematopathology, Harbin Children’s Hospital, Harbin, China

how to cite: Dong Y, Fei P, Han Y, Guo L. Characterization of Fecal Microbiota, Short-Chain Fatty Acids and Lactic Acid Concentrations in 5 - 8-Year-Old Children with Cow Milk Protein Allergy. Iran J Pediatr. 2018;28(4):e64638. https://doi.org/10.5812/ijp.64638.

Abstract

Background:

Cow milk protein allergy (CMPA) has been reported to be associated with the alterations of gut microbiota. However, changes in the fecal microbial community for children aged 5 - 8 years with CMPA have not been clearly elucidated.

Objectives:

The aim of the study is to assess the differences and relationship between fecal microbiota and levels of fecal short-chain fatty acids (SCFAs) and lactic acid in 5 - 8-year-old children with CMPA, compared with the normal children.

Methods:

We enrolled 6 CMPA children and 8 healthy children from Harbin children’s hospital (Heilongjiang, China). Total bacterial DNA of all fecal samples was extracted and prepared for bioinformatics analysis by 16S rDNA amplicon sequencing. Lactic acids and SCFAs were measured by the high-performance liquid chromatography (HPLC) technique.

Results:

Compared with healthy individuals, the fecal microbial community of CMPA children showed a significant increment of Firmicutes and depletion of Proteobacteria (P < 0.05). At the class level, the number of Clostridia was observed to be significantly increased in the CMPA children’s fecal microbiota. Meanwhile, the diversity of Ruminococcaceae and Subdoligranulum was higher in the fecal microbiome of CMPA children than in healthy children (P < 0.05). Furthermore, linear discriminant analysis of effect size (LEfSe) revealed some significantly enriched taxa of CMPA children and control subjects. There are 40 and 5 enriched taxa in the fecal microbiota of children with CMPA and healthy children respectively. The concentration of total SCFAs was lower in feces of CMPA children than in healthy children. Lactic acid had a significantly negative correlation with the abundance of Bacteroides.

Conclusions:

Alterations in the composition of fecal microbiota and SCFAs in the 5 - 8-year-old children with CMPA may reduce immune function and potentially promote this allergy<b>.</b>

1. Background

Food allergy is one of the most common allergies in infants and children, the prevalence rate of which is increased with time for foods including cow milk protein (1). The prevalence of cow milk protein allergy (CMPA) has been reported to be increased from 1.6% to 3.5% over a decade (2). Although the majority of CMPA children can recover spontaneously by the age of 5 years, a recent report has shown that about 10% of CMPA children continued allergic symptoms until more than 5 years old (3). In addition, CMPA children over 5 years old are often accompanied by other allergic reactions, including asthma and rhinoconjunctivitis (4, 5), and so on.

In recent years, the etiologies of atopic diseases have been hypothesized to be involved in gut microbiota, because it has been observed that some disruptions occur in the gut microbiota between healthy and allergic children (6-8). Compared with healthy children, more Clostridia (9), fewer Bifidobacteria (10, 11) and lower Lactobacilli (12) were detected in the gut microbiota of allergic children. Especially, Bifidobacterium in the feces of infants and children were reported to be associated with the development of food allergy (13). Similarly, our previous study analyzed fecal microbiota in 5 - 8-year-old children with and without CMPA by using polymerase chain reaction-denaturing gradient gel electrophoresis (PCR-DGGE), and revealed that the diversity of Bifidobacterium group was lower in the fecal microbiota of CMPA children than in healthy controls (14). Given an increased prevalence and a long-term persistence of CMPA, it is necessary to conduct a further study on the fecal microbial composition of CMPA children.

High-throughput sequencing technology is commonly used to analyze comprehensively the richness and diversity of the gut microbiota. It can identify approximately 80% more gut microbiota which was hardly identified by culture-based microbiological methods (15). The fecal microbiota characteristics had been known and clarified using the high-throughput sequencing technology, particularly in the allergic diseases, including food allergy (16), atopic eczema (17), and allergic atopic dermatitis (18). The application of high-throughput sequencing technique to the study of gut microbiota in CMPA children has been less fruitful.

Short-chain fatty acids (SCFAs) are major end products metabolized by mainly anaerobic bacteria in the intestine, they play the key role for adjusting the function of the host metabolism and regulating of immune responses (19). A higher level of acetic acid has been observed in allergic 4-year-old children compared with healthy children (20). Thompson-Chagoyan et al. found that fecal butyric acid and branched-chain short fatty acids were higher in CMPA infants than in healthy infants (21). However, a further thorough study on the changes of SCFAs and lactic acid in feces of CMPA children was limited.

In this study, we focus on the alteration of intestinal microflora structure and changes of SCFAs of CMPA children aged 5 to 8 years.

2. Objectives

The study aims to comprehensively compare the composition of fecal microbiota and levels of fecal SCFAs and lactic acid between CMPA children at the age of 5 - 8 years and healthy controls, and analyze the interaction between fecal microflora and its metabolites, namely SCFAs and lactic acid.

3. Methods

3.1. Population Enrollment and Sample Collection

During a 7-month period (March 2016 to October 2016) at the Harbin children’s hospital (Heilongjiang, China), we enrolled 6 CMPA children aged 5 - 8 (mean 6.4) years with about equal number of males and females. The CMPA children in this study showed some symptoms of allergic disease such as angioedema, pruritus, diarrhea or eczema, the details of which are given in supplementary file (Appendix 1). Moreover, CMPA was diagnosed by three tests: a skin prick test, analysis of serum samples and a double-blind placebo-controlled standard cow’s milk formula challenge (5, 22, 23). The patient was considered CMPA when the following criteria were met (21): (i) an obvious immediate hypersensitivity to CMP; (ii) positive skin prick test, specific IgE to CMP; (iii) positive CMP challenge test. Meanwhile, 8 age- and sex-matched controls were recruited from healthy children in Harbin city with no history of CMPA or any other allergies. According to their sensitization, the children were divided into two groups: CMPA group (C) and healthy group (H). In addition, all participants did not ingest antibiotics within 4 weeks before the sampling.

All fecal samples were obtained from all participating children immediately after defecation, and transferred into sterile plastic tubes by physicians. Then the samples were transported quickly to the lab and stored at -80°C until extraction of genomic DNA.

3.2. DNA Extraction, 16S rRNA Gene Amplification and Sequencing

DNA was extracted using a TIANamp Stool DNA Kit (Tiangen, China, cat#DP328-02) according to the manufacturer’s instructions. In order to improve the efficiency of cell wall lysis of gram-positive bacteria in the feces, the samples were incubated at 95°C for 5 minutes. The extracted DNA was stored at -20°C until processed (24).

DNA was amplified with the primers 515 F (5’-GTG CCA GCM GCC GCGG-3’) and 806 R (5’- GGA CTA CNN GGG TAT CTA AT -3’) for the V3-V4 hypervariable regions of 16S rRNA gene (25). All PCR reaction mixture comprised with 15 μL of Phusion® High-Fidelity PCR Master Mix (New England Biolabs), 0.2 μM of forward and reverse primers, and about 10 ng DNA template. The PCR conditions were as follows: thermal cycling consisted of initial denaturation at 98°C for 1 minute, followed by 30 cycles of denaturation at 98°C for 10 seconds, annealing at 50°C for 30 seconds, and elongation at 72°C for 30 seconds. Finally 72°C for 5 minutes. PCR products were detected by electrophoresis on 2% agarose gel. Then, the mixture of PCR products was purified using GeneJET Gel Extraction Kit (Thermo Scientific). After purification, samples with bright main strip were chosen to analyze using the Illumina high-throughput sequencing method (26).

3.3. Bioinformatics Analysis and Statistical Analysis

The paired-end reads were assigned to each sample according to the unique barcodes. Sequences were analyzed using QIIME (Quantitative Insights Into Microbial Ecology) software package. First, reads were filtered by QIIME quality filters (26). Then sequences with ≥ 97% similarity were assigned to the same OTUs by utilizing UPARSE software V7.0.1001. The representative sequences were selected for each OTU and used the Ribosomal Database Project (RDP) classifier to annotate taxonomic information (27). Linear discriminant analysis (LDA) effect size (LEfSe) was used to find significantly different microbial community among the groups, and LDA was performed to estimate the effect size of each feature. The correlation between the abundance of fecal microbiota and SCFAs and lactic acid was evaluated by Spearman correlation analysis.

3.4. Determination of the Fecal Lactic Acid and SCFAs

The concentrations of the fecal lactic acid and SCFAs were analyzed by high-performance liquid chromatography (HPLC) according to Dostal et al. (28). 500 mg faces were mixed with 2 mL of 0.5 mmol/L H2SO4 and centrifuged for 20 minutes (11000 r/min, 4°C). The supernatants were filtered via 0.22 μm organic filter to 1.5 mL vials for HPLC.

HPLC was performed using a COSMOSIL® 5C1818-AR-II (4.6 × 250 mm) column, with the flow rate of 0.6 mL/min at 35°C and the flowing phase 50% methanol: 0.1% acetic acid (1:1). The ultraviolet detector operated at λ = 210, and the total running time was set at 20 minutes. The concentrations of samples were analyzed according to the standard curves and the data was expressed as mmol/kg feces.

3.5. Ethical Considerations

This study was approved by the northeast agricultural university ethics committee of Heilongjiang province, China and informed consents were obtained from the children’s parents. In addition, the study protocol is consistent with the ethical guidelines of the 1975 declaration of Helsinki.

3.6. Statistical Analysis

Independent-samples t-test was used to analyze statistical significance in fecal microbial abundance and alpha diversity, and P < 0.05 was considered statistically different. Analyses were performed using the software SPSS 16.0 for windows.

4. Results

4.1. Sequencing Metrics and Taxonomy

16S rDNA amplicon sequencing was performed to characterize the bacterial lineages presented in the fecal microbiota of the 14 children. A total of 583,198 high-quality sequences reads remained for analysis, after removing low quality and short length sequence reads, chimeric sequences and sequences only appeared once in the dataset. The average length of validated sequences reads was 250 bp, which were clustered into OTUs with a similarity cutoff of 97%. A total of 5, 922 OTUs were obtained and each individual averaged 423 OTUs.

4.2. Microbial Richness and Biodiversity

The alpha diversity of fecal microbiota was analyzed by using Chao1, ACE, Simpson, Shannon, and coverage indices (Table 1). The OTU richness was estimated by the Chao1 estimator and ACE, and the Shannon and Simpson index reflected the diversity of microbial community. Good’s coverage, which presented a measure of sampling completeness, was above 99.4% in all samples. These indexes had no significant difference between the two groups by t-test (P > 0.05).

Table 1.

Analysis of Alpha Diversity at 97% Similaritya

GroupGood’s Coverage, %Richness EstimatorDiversity Index
Chao 1ACESimpsonShannon
C99.585695820.8944.712
H99.733343430.6863.264

4.3. Abundance of Fecal Bacterial Composition

Bioinformatic analysis of 16S rDNA sequence data revealed that there were some differences in the fecal microbiota of children with CMPA and healthy children. More than 96.6% of the sequences in all samples were found to belong to the four most populated bacterial phyla, Actinobacteria, Bacteroidetes, Firmicutes, and Proteobacteria (Figure 1). Bacteroidetes and Firmicutes were more presented in the group C than in the control group (18.84% vs. 10.68%, 59.59% vs. 35.88% respectively), whereas Actinobacteria and Proteobacteria were more abundant in the group H than in the group C (11.74% vs. 9.46%, 40.21% vs. 8.77% respectively). Statistical analysis using independent-samples t-test indicated that Proteobacteria (P = 0.047) and Firmicutes (P = 0.046) significantly differentiate the CMPA children from the healthy children.

The abundance of Actinobacteria, Proteobacteria, Bacteroidetes and Firmicutes in group C and group H (*P < 0.05 )
The abundance of Actinobacteria, Proteobacteria, Bacteroidetes and Firmicutes in group C and group H (*P < 0.05 )

Among the total classes identified in the fecal microbiota, six predominant bacterial classes were detected in two groups. Predominant classes were defined as those accounting for > 1% of the total DNA sequences (Figure 2). The class Clostridia was more dominant in the fecal microbiota of CMPA children (53.50%), Gammaproteobacteria was more dominant in the fecal microbiota of healthy children (38.28%). Compared with group H, the Clostridia class had a significant increment in group C (P < 0.05).

The differences of bacterial distribution at class levels from group C and group H (A indicates group C, B indicates group H)
The differences of bacterial distribution at class levels from group C and group H (A indicates group C, B indicates group H)

The study also found some differences in abundance at the family level between the two groups (Figure 3) Ruminococcaceae (38.01%) were the dominant bacteria in the group C, while Enterobacteriaceae and Streptococcaceae were the dominant bacteria in the group H (28.97% and 25.56%). The percentages of Veillonellaceae and Rikenellaceae were very low (3.38% vs 1.37%, 1.37% vs. 1.01%) in the two groups. Group C had significantly higher numbers of Ruminococcaceae compared with the group H (38.01% vs 10.33%; P = 0.015). Furthermore, Burkholderiaceae, Nannocystaceae, Shewanellaceae, Thermomonosporaceae and Flavobacteriaceae were present exclusively in the group C, but vadinBB60, Methylophilaceae and Dietziaceae were only found in the group H (Figure 4).

The differences in microbial composition at family level in two groups (A indicates group C, B indicates group H)
The differences in microbial composition at family level in two groups (A indicates group C, B indicates group H)
The number of family belonging to the major phyla differing between group C and group H
The number of family belonging to the major phyla differing between group C and group H

There were 11 genera which accounted for more than 2% in each group (Table 2). In the group C, Bacteroides (15.84%, belong to Bacteroidetes) was predominant followed by Subdoligranulum (12.39%, belong to Firmicutes), Faecalibacterium (11.15%, belong to Firmicutes), Bifidobacterium (8.27%, belong to Actinobacteria) and others. In the feces of the healthy children, Escherichia-Shigella (28.78%, belong to Proteobacteria) was predominant followed by Bifidobacterium (10.47%, belong to Actinobacteria), Bacteroides (8.27%, belong to Bacteroidetes), Klebsiella (5.4 %, belong to Proteobacteria) and others. Subdoligranulum (Firmicutes) found in the feces of the group C was more abundant compared with group H (12.39% vs. 2.71%, P < 0.05).

Table 2.

The Distribution of Bacterial Genera Present in Fecal Samples of CMPA and Healthya

GenusGroup CGroup H
Escherichia-Shigella2.8628.78
Bifidobacterium8.2710.47
Bacteroides15.848.27
Klebsiella2.605.40
Lachnoclostridium5.424.03
Faecalibacterium11.152.78
Ruminococcus_25.593.16
Subdoligranulum12.392.71
Coprostanoligenes_group3.462.12
Pseudobutyrivibrio3.76-
Ruminococcus_12.36-
Enterobacter-2.76
Streptococcus-3.32

4.4. Comparison of Relative Abundance Through LDA Score and Cladogram Analysis

To identify the effect size of different abundant taxa between group C and group H, a metagenomic biomarker discovery approach (LEfSe) was utilized to analyze (29). The cladogram (Figure 5) showed that many co-exist taxa (shown in yellow) in two groups, and specific bacterial taxa were highlighted by colored circles and shaded areas (CMPA samples were shown in red and healthy samples were shown in green). The predominant taxa (LDA values > 2) of the two groups were displayed in Figure 6. The distribution histogram of LDA values showed 45 significantly enriched taxa, 40 and 5 taxa enriched in the feces microbiota of group C and group H respectively. Firmicute, Ruminococcaceae, and Subdoligranulum were mainly enriched taxa in group C. Family XI and Gemella were mainly enriched taxa in group H.

A cladogram representation of data is shown in panel (Red color, specific bacterial taxa in group C. Green color, specific bacterial taxa in group H. Yellow color, co-exist taxa in two groups. The brightness of each dot is proportional to its effect size).
A cladogram representation of data is shown in panel (Red color, specific bacterial taxa in group C. Green color, specific bacterial taxa in group H. Yellow color, co-exist taxa in two groups. The brightness of each dot is proportional to its effect size).
Linear discriminant analysis of effect size (LEfSe) identifies the most differentially abundant taxons between group C and group H (The enriched taxa in group C are indicated with a negative LDA score (red), and the enriched taxa in group H are indicated with a positive linear discriminant analysis (LDA) score (green). Only taxa meeting an LDA significant threshold of > 2 are shown).
Linear discriminant analysis of effect size (LEfSe) identifies the most differentially abundant taxons between group C and group H (The enriched taxa in group C are indicated with a negative LDA score (red), and the enriched taxa in group H are indicated with a positive linear discriminant analysis (LDA) score (green). Only taxa meeting an LDA significant threshold of > 2 are shown).

4.5. Analysis of Fecal Lactic Acid and SCFAs

Total SCFAs in feces contents in group C was significantly decreased than in group H (221.33 mmol/kg vs. 255.95 mmol/kg), shown in Figure 7. Acetic acid was the dominant SCFA, followed by isobutyric and propionic acids. However, the difference between each of the SCFA and lactic acid in the two groups was not significant.

Concentrations of short-chain fatty acids (SCFAs) and lactic acid in two groups
Concentrations of short-chain fatty acids (SCFAs) and lactic acid in two groups

The relationship between the composition and abundance of fecal microbiota and SCFAs and lactic acid was analyzed. At family level (Appendix 2 in supplementary file), Lachnospiraceae 1 and Christensenellaceae were positively correlated with acetic acid, butyric acids and total SCFAs. Thermotogaceae were positively correlated with butyric acid, isobutyric acid and total SCFAs. Ruminococcaceae were positively correlated with acetic acid and Pasteurellaceae were significantly and positively correlated with negatively correlated with lactic acid (P < 0.01). Bacteroidaceae were negatively correlated with lactic acid. At genus level (Figure 8), Ruseburia, Christensenellaceae R-7 group, Lachnospiraceae NK4A136 group, and Ruminococcaceae UCG-002 were positively correlated with acetic acid, butyric acid and total SCFAs. Haemophilus was positively correlated with propionic acid, but Escherichia-Shigella was negatively correlated with propionic acids. Moreover, it was noteworthy that there was a negative relationship between Bacteroides and lactic acid.

Spearman rank correlation between bacterial genus and short-chain fatty acids (SCFAs) and lactic acid
Spearman rank correlation between bacterial genus and short-chain fatty acids (SCFAs) and lactic acid

5. Discussion

The majority of bacteria in human gut microbiota come from Proteobacteria, Firmicutes, Actinobacteria and Bacteroidetes (30), which was consistent with our study. The two groups of the research showed that most of sequences belong to the four most populated bacterial phyla mentioned above. Importantly, we found a significant increment of Firmicutes and depletion of Proteobacteria in fecal microbiota of CMPA children (P < 0.05). The facts are consistent with the results of a previous study showing that the high levels of Firmicutes in fecal microbiota of children with other food allergies such as egg allergy, wheat allergy, and peanut allergy (13). Moreover, Bunyavanich et al. found that gut microbiome composition in infancy (3 - 6 months) was associated with CMPA resolution, and the high abundance of Firmicutes was linked to immune tolerance in allergy (31).

The Clostridia class was found to be more enriched in group C than in group H in the present study. The high abundance of Clostridia was previously detected in a culture-based study, in which Clostridia had a connection with the development of allergy (32). Compared with germ-free mice, mice colonized with Clostridia had higher concentrations of fecal IgA. The colonization of Clostridia in gut microbiota reduced IL-22 production, which reduced the concentration of serum allergen (33). These conclusions may suggest that Clostridia took a vital part in the regulation of food allergy.

In a case-control PCR-DGGE-based study, Guo et al. found a lower degree of diversity of Bifidobacterium group in the feces of children with CMPA. Meanwhile, Bacteroides, Clostridium, and Escherichia coli were detected to be more abundant in the fecal microbiota of CMPA children than in healthy controls (14). In an elegant FISH-FC-based study, significantly higher proportions of the Clostridium coccoides group and Atopobium cluster were found in CMPA infant feces (21). However, these bacterial genera to be mentioned in the above studies had no significant difference between group C and group H in the study. The genus, Subdoligranulum, which belongs to the class Clostridia, was more abundant in the CMPA children than in the healthy children. The result implied that Clostridia played an important role in allergy maybe due to the presence of Subdoligranulum.

There was a lower concentration of total SCFAs in CMPA children than in control children in the study. SCFAs were mainly produced by bacteria in the intestine (19). The pattern of SCFAs will be changed with the altered gut microbiota (34). Changes in the fecal microbiota of CMPA children may be one of the reasons leading to a lower concentration of total SCFAs. Therefore, increasing the concentrations of SCFAs may be carried out by changing the composition of gut microbiota.

To our knowledge, some researches analyzed the correlation between gut microbiota and SCFAs, especially in allergy diseases (21, 35). The relationship between some special strains in the gut and SCFAs were clarified. Bacteroides are relatively sensitive to acidic pH, which grow poorly at pH 5.5 but grow well at pH 6.7 (36). An increase in lactic acid will lower the pH and may thus affect the growth of the Bacteroides. So it is reasonable to find that Bacteroides was negatively related to lactic acid in the study. In addition, the results also showed the correlation between SCFAs and bacteria at each taxonomic level of the intestine. Obviously, SCFAs and fecal microbiota were closely related and had a role to play in the remission of allergic diseases. Thus, the treatment that combined increased SCFAs and altered gut microbiota may be more helpful to the resolution of CMPA, and the further research on this may be needed.

5.1. Conclusion

The results of the study revealed the alterations of fecal microbiota in the CMPA of 5 - 8 years old children according to the culture-independent method. The fact that the concentration of total SCFAs was significantly decreased in the feces of CMPA children enriched our knowledge of the relationship between gut microorganisms and their metabolic productions. Furthermore, lactic acid was found to correlate negatively with Bacteroides. The results provide a possible basis for the further study of the relationship between CMPA and gut microbiota in 5 - 8-year-old children.

Acknowledgements

References

  • 1.

    Rachid R, Chatila TA. The role of the gut microbiota in food allergy. Curr Opin Pediatr. 2016;28(6):748-53. [PubMed ID: 27749359]. https://doi.org/10.1097/MOP.0000000000000427.

  • 2.

    Fiocchi A, Schunemann HJ, Brozek J, Restani P, Beyer K, Troncone R, et al. Diagnosis and Rationale for Action Against Cow's Milk Allergy (DRACMA): a summary report. J Allergy Clin Immunol. 2010;126(6):1119-28. [PubMed ID: 21134569]. https://doi.org/10.1016/j.jaci.2010.10.011.

  • 3.

    Levy Y, Segal N, Garty B, Danon YL. Lessons from the clinical course of IgE-mediated cow milk allergy in Israel. Pediatric Allergy and Immunology. 2007;18(7):589-93. https://doi.org/10.1111/j.1399-3038.2007.00562.x.

  • 4.

    Host A, Halken S, Jacobsen HP, Christensen AE, Herskind AM, Plesner K. Clinical course of cow's milk protein allergy/intolerance and atopic diseases in childhood. Pediatr Allergy Immunol. 2002;13 Suppl 15:23-8. [PubMed ID: 12688620].

  • 5.

    Meglio P, Bartone E, Plantamura M, Arabito E, Giampietro PG. A protocol for oral desensitization in children with IgE-mediated cow's milk allergy. Allergy. 2004;59(9):980-7. [PubMed ID: 15291907]. https://doi.org/10.1111/j.1398-9995.2004.00542.x.

  • 6.

    Nakayama J, Kobayashi T, Tanaka S, Korenori Y, Tateyama A, Sakamoto N, et al. Aberrant structures of fecal bacterial community in allergic infants profiled by 16S rRNA gene pyrosequencing. FEMS Immunol Med Microbiol. 2011;63(3):397-406. [PubMed ID: 22029688]. https://doi.org/10.1111/j.1574-695X.2011.00872.x.

  • 7.

    Azad MB, Konya T, Guttman DS, Field CJ, Sears MR, HayGlass KT, et al. Infant gut microbiota and food sensitization: associations in the first year of life. Clin Exp Allergy. 2015;45(3):632-43. [PubMed ID: 25599982]. https://doi.org/10.1111/cea.12487.

  • 8.

    Hua X, Goedert JJ, Pu A, Yu G, Shi J. Allergy associations with the adult fecal microbiota: Analysis of the American Gut Project. EBioMedicine. 2016;3:172-9. [PubMed ID: 26870828]. [PubMed Central ID: PMC4739432]. https://doi.org/10.1016/j.ebiom.2015.11.038.

  • 9.

    Penders J, Thijs C, van den Brandt PA, Kummeling I, Snijders B, Stelma F, et al. Gut microbiota composition and development of atopic manifestations in infancy: the KOALA Birth Cohort Study. Gut. 2007;56(5):661-7. [PubMed ID: 17047098]. [PubMed Central ID: PMC1942165]. https://doi.org/10.1136/gut.2006.100164.

  • 10.

    Kalliomaki M, Kirjavainen P, Eerola E, Kero P, Salminen S, Isolauri E. Distinct patterns of neonatal gut microflora in infants in whom atopy was and was not developing. J Allergy Clin Immunol. 2001;107(1):129-34. [PubMed ID: 11150002]. https://doi.org/10.1067/mai.2001.111237.

  • 11.

    Bjorksten B, Sepp E, Julge K, Voor T, Mikelsaar M. Allergy development and the intestinal microflora during the first year of life. J Allergy Clin Immunol. 2001;108(4):516-20. [PubMed ID: 11590374]. https://doi.org/10.1067/mai.2001.118130.

  • 12.

    Toh ZQ, Anzela A, Tang ML, Licciardi PV. Probiotic therapy as a novel approach for allergic disease. Front Pharmacol. 2012;3:171.

  • 13.

    Chen CC, Chen KJ, Kong MS, Chang HJ, Huang JL. Alterations in the gut microbiotas of children with food sensitization in early life. Pediatr Allergy Immunol. 2016;27(3):254-62. [PubMed ID: 26663491]. https://doi.org/10.1111/pai.12522.

  • 14.

    Guo L, Bai H, Dong Y, Huang DX, Zhang X, Gong S, et al. Comparative Analysis of Fecal Microbiota in 5 - 8-Year-old Children with and without Cow Milk Protein Allergy. Iran J Pediatr. 2016;26(6). https://doi.org/10.5812/ijp.6397.

  • 15.

    Aitoro R, Paparo L, Amoroso A, Di Costanzo M, Cosenza L, Granata V, et al. Gut Microbiota as a Target for Preventive and Therapeutic Intervention against Food Allergy. Nutrients. 2017;9(7). [PubMed ID: 28657607]. [PubMed Central ID: PMC5537787]. https://doi.org/10.3390/nu9070672.

  • 16.

    Ling Z, Li Z, Liu X, Cheng Y, Luo Y, Tong X, et al. Altered fecal microbiota composition associated with food allergy in infants. Appl Environ Microbiol. 2014;80(8):2546-54. [PubMed ID: 24532064]. [PubMed Central ID: PMC3993190]. https://doi.org/10.1128/AEM.00003-14.

  • 17.

    Abrahamsson TR, Jakobsson HE, Andersson AF, Bjorksten B, Engstrand L, Jenmalm MC. Low diversity of the gut microbiota in infants with atopic eczema. J Allergy Clin Immunol. 2012;129(2):434-40. 440 e1-2. [PubMed ID: 22153774]. https://doi.org/10.1016/j.jaci.2011.10.025.

  • 18.

    Drell T, Larionova A, Voor T, Simm J, Julge K, Heilman K, et al. Differences in Gut Microbiota Between Atopic and Healthy Children. Curr Microbiol. 2015;71(2):177-83. [PubMed ID: 25869237]. https://doi.org/10.1007/s00284-015-0815-9.

  • 19.

    Vonk RJ, Priebe M, Meijer K, Venema K, Roelofsen H. The interaction of short-chain fatty acids (SCFA) with adipose tissue; relevance for systemic inflammation. Gastroenterology. 2011;140(5):S-860.

  • 20.

    Sandin A, Braback L, Norin E, Bjorksten B. Faecal short chain fatty acid pattern and allergy in early childhood. Acta Paediatr. 2009;98(5):823-7. [PubMed ID: 19173682]. https://doi.org/10.1111/j.1651-2227.2008.01215.x.

  • 21.

    Thompson-Chagoyan OC, Fallani M, Maldonado J, Vieites JM, Khanna S, Edwards C, et al. Faecal microbiota and short-chain fatty acid levels in faeces from infants with cow's milk protein allergy. Int Arch Allergy Immunol. 2011;156(3):325-32. [PubMed ID: 21720179]. https://doi.org/10.1159/000323893.

  • 22.

    Huang F, Kim JS. IgE-mediated cow's milk allergy in children. Curr Allergy Asthma Rep. 2012;12(6):630-40. [PubMed ID: 22847747]. https://doi.org/10.1007/s11882-012-0289-9.

  • 23.

    El-Agamy EI. The challenge of cow milk protein allergy. Small Ruminant Research. 2007;68(1):64-72.

  • 24.

    Li WH, Zhang H, Guo Q, Wu XD, Xu ZS, Dang CX, et al. Detection of SNCA and FBN1 methylation in the stool as a biomarker for colorectal cancer. Dis Markers. 2015;2015. https://doi.org/10.1155/2015/657570.

  • 25.

    Caporaso JG, Lauber CL, Walters WA, Berg-Lyons D, Lozupone CA, Turnbaugh PJ, et al. Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. Proc Natl Acad Sci U S A. 2011;108 Suppl 1:4516-22. [PubMed ID: 20534432]. [PubMed Central ID: PMC3063599]. https://doi.org/10.1073/pnas.1000080107.

  • 26.

    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;7(5):335-6. [PubMed ID: 20383131]. [PubMed Central ID: PMC3156573]. https://doi.org/10.1038/nmeth.f.303.

  • 27.

    Wang Q, Garrity GM, Tiedje JM, Cole JR. Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol. 2007;73(16):5261-7. [PubMed ID: 17586664]. [PubMed Central ID: PMC1950982]. https://doi.org/10.1128/AEM.00062-07.

  • 28.

    Dostal A, Baumgartner J, Riesen N, Chassard C, Smuts CM, Zimmermann MB, et al. Effects of iron supplementation on dominant bacterial groups in the gut, faecal SCFA and gut inflammation: a randomised, placebo-controlled intervention trial in South African children. Br J Nutr. 2014;112(4):547-56. [PubMed ID: 24916165]. https://doi.org/10.1017/S0007114514001160.

  • 29.

    Segata N, Izard J, Waldron L, Gevers D, Miropolsky L, Garrett WS, et al. Metagenomic biomarker discovery and explanation. Genome Biol. 2011;12(6):R60. [PubMed ID: 21702898]. [PubMed Central ID: PMC3218848]. https://doi.org/10.1186/gb-2011-12-6-r60.

  • 30.

    Backhed F, Ley RE, Sonnenburg JL, Peterson DA, Gordon JI. Host-bacterial mutualism in the human intestine. Science. 2005;307(5717):1915-20. [PubMed ID: 15790844]. https://doi.org/10.1126/science.1104816.

  • 31.

    Bunyavanich S, Shen N, Grishin A, Wood R, Burks W, Dawson P, et al. Early-life gut microbiome composition and milk allergy resolution. J Allergy Clin Immunol. 2016;138(4):1122-30. [PubMed ID: 27292825]. [PubMed Central ID: PMC5056801]. https://doi.org/10.1016/j.jaci.2016.03.041.

  • 32.

    Sepp E, Julge K, Mikelsaar M, Bjorksten B. Intestinal microbiota and immunoglobulin E responses in 5-year-old Estonian children. Clin Exp Allergy. 2005;35(9):1141-6. [PubMed ID: 16164439]. https://doi.org/10.1111/j.1365-2222.2005.02315.x.

  • 33.

    Stefka AT, Feehley T, Tripathi P, Qiu J, McCoy K, Mazmanian SK, et al. Commensal bacteria protect against food allergen sensitization. Proc Natl Acad Sci USA. 2014;111(36):13145-50. [PubMed ID: 25157157]. [PubMed Central ID: PMC4246970]. https://doi.org/10.1073/pnas.1412008111.

  • 34.

    Canfora EE, Jocken JW, Blaak EE. Short-chain fatty acids in control of body weight and insulin sensitivity. Nat Rev Endocrinol. 2015;11(10):577-91. [PubMed ID: 26260141]. https://doi.org/10.1038/nrendo.2015.128.

  • 35.

    Francavilla R, Calasso M, Calace L, Siragusa S, Ndagijimana M, Vernocchi P, et al. Effect of lactose on gut microbiota and metabolome of infants with cow's milk allergy. Pediatr Allergy Immunol. 2012;23(5):420-7. [PubMed ID: 22435727]. https://doi.org/10.1111/j.1399-3038.2012.01286.x.

  • 36.

    Duncan SH, Louis P, Thomson JM, Flint HJ. The role of pH in determining the species composition of the human colonic microbiota. Environ Microbiol. 2009;11(8):2112-22. [PubMed ID: 19397676]. https://doi.org/10.1111/j.1462-2920.2009.01931.x.