- 1Danone Nutricia Research, Singapore, Singapore
- 2Laboratory of Microbiology, Wageningen University, Wageningen, Netherlands
- 3Danone Nutricia, São Paulo, Brazil
- 4Department of Pediatrics, Universidade Federal Do Triângulo Mineiro, Uberaba, Brazil
- 5Department of Pediatrics, Universidade Federal Do Rio Grande Do Norte, Natal, Brazil
- 6Faculdade De Medicina Do ABC, São Paulo, Brazil
- 7Danone Nutricia Research, Utrecht, Netherlands
- 8Escola Paulista de Medicina, Universidade Federal de São Paulo, São Paulo, Brazil
Recent studies have demonstrated that gut microbiota development is influenced by human biogeographic factors such as race, ethnicity, diet, lifestyle or culture-specific variations, and other environmental influences. However, biogeographic variation in gut microbiota assembly remains largely unexplored in Latin America. In this paper, we compared food recall information and microbiota composition of toddlers living in geographically separated urban populations within four states of Brazil. 16S RNA gene sequencing revealed that alpha diversity was similar between the four different populations. Gut microbiota compositions were dominated by members of the phyla Firmicutes and Bacteroidetes, resembling a more adult-like microbiota as compared with those of Western European toddlers of similar age. These findings suggest that inter-individual and nutrition-induced differences were apparent in the fecal microbiota. We conclude that urban dietary pattern plays a larger role in influencing the gut microbiota composition than do biogeographic factors.
Introduction
As interest in the gut microbiome has rapidly increased in the last decade, attention has turned to the role of early life microbiota development in later life health (1–3). It is now widely recognized that the gut microbiome plays an important role in human metabolic, immunity, and even cognitive development. Nevertheless, the ecological succession and function of the gut microbiome during childhood are not completely understood, as research mainly focuses on infants. Although it has been suggested that the gut microbiome reaches a relatively stable, adult-like state after the first 1–3 years of life (4–6), there are other studies that suggest that its composition continues to develop into the teenage years (7–9).
Early childhood is the most critical phase in various aspects of human development. Factors such as nutrition, environmental changes, and infection can influence the maturation of the immune system. This in turn may severely impact a child's resilience to pathogenic infections. The microbiota plays critical roles in the maturation and priming of the human innate and adaptive immune systems, while the immune system reciprocally orchestrates the stability and function of the host–microbiota symbiosis. However, poor nutrition can directly and indirectly compromise immune function and subsequent susceptibility to infection (10).
Several studies have investigated the nutritional status and feeding patterns of Brazilian children over a broad age range (11–13). The most prevalent form of deficiencies identified were incidences of low vitamin and micronutrient status (14). This observation has urged medical professionals and policymakers to better engage and educate parents and caregivers on the importance of establishing good eating habits that will in the long term benefit their children's health (15).
The gut microbiota also act as a protective physical barrier to pathogenic bacteria, harvest energy from food, and are crucial to synthesis of metabolites required for cellular functions (16–18). Many research groups have focused their investigation on factors such as mode of delivery, type of feeding, and environmental impacts in the 1st year of life (19–22). After weaning, the child is exposed to greater range of food, which slowly evolves the gut microbiota toward adult-like profile.
A report by De Filippo et al. (23, 24) suggests that sanitation and increased amounts of readily accessible acellular nutrients typical for a modern Western-type diet and lifestyle might be key to the selection of microbial dynamics and functionality, which differentiates urban and rural living. Traditionally, the main source of food comes from unprocessed agriculture crops, and progressive urbanization and industrialization have increased accessibility to processed foods. This transition to modern lifestyles has been associated with reduced intestinal diversity and changing microbiota (25). A study that evaluated Brazilian data from the National Demographic and Health Survey (2006–2007) reported that the diet of children in developed urbanized areas consists of higher processed foods and more foods not recommended for their age, like sweets and soft drinks (soda), as compared with children from rural areas. This is particularly prominent in regions such as South, Southeastern, and Midwest Brazil relative to the North and Northeast.
In this cross-sectional cohort study, we conduct a descriptive analysis of microbiota composition and dietary intake by simple food recall information in healthy, young children aged 1–3 years (n = 99) within four urban cities in four states in Western Brazil.
Materials and Methods
Subjects and Study Population
This cross-sectional observational multicenter study included a total of 99 children aged between 1 and 3 years, enrolled full time in day care centers from different regions of Brazil [Santo André (SAN), State of São Paulo; Uberaba (UBE), State of Minas Gerais; Porto Alegre (POA), State of Rio Grande do Sul; and Natal (NAT), State of Rio Grande do Norte]. Subjects with missing data or no fecal samples were excluded. A 2-day recall food frequency questionnaire (FFQ) was used to collect information on health and diet of the children, and non-compulsory fecal samples were collected for intestinal microbiota composition and function analyses. Participants with prior chronic diseases and congenital malformations or who received antibiotics <31 days prior to sample collection were excluded. Participants' samples that were collected without traceable subject information (e.g., geographical origin) and/or no fecal sample collected at all were excluded from the subsequent analyses. Full written informed consent was obtained from all included participants parents prior to collection of biological samples and any personal information. This study was approved by the ethics committee (CEP-FMABC, Comitê de Ética, Santo Andre) and registered at ClinicalTrials.gov (November 1, 2016; NCT02950740).
Sample Collection and DNA Extraction
Fecal samples were collected at a single cross-sectional time point from subjects in day care centers in the four cities and frozen immediately at 4°C. Samples were then transported in insulated storage containers and stored at −80°C at São Paulo in GC-2 company laboratory until further analysis.
Samples were processed using QIAGEN DNA Stool Mini-Kit + RNase step (QIAGEN, Hilden, Germany), modified with an additional bead-beating step. Approximately 200 mg of 0.1-mm glass beads were added to 200-mg stool sample and resuspended. Samples were processed using QIAGEN DNA Stool Mini-Kit + RNase step (QIAGEN, Hilden, Germany), modified with an additional bead-beating step. Approximately 200 mg of 0.1-mm glass beads were added to 200-mg stool sample and resuspended in QIAGEN ASL buffer; bead-beating was done using FastPrep24 (M.P. Biomedicals, USA) for three repetitions of 1-min bead-beating with 5-min incubation on ice. Samples were then heated at 95°C for 15 min before centrifugation at 20,000 × g for 1 min. The supernatant was transferred into a clean tube containing an InhibitEX tablet and vortexed to mix proceeding with manufacturer instructions. Isolated fecal DNA was eluded in 50 μl of AE buffer, quantified using a NanoDrop 2000 (Thermo Scientific, DE, USA), and stored at −20°C before used for further analyses. The remaining fecal samples were resuspended in phosphate-buffered saline (PBS) for short-chain fatty acid (SCFA) analysis.
16S rDNA Sequencing and Bioinformatics Analysis
From the purified fecal DNA extracts, the V3–V4 hypervariable regions of the bacterial 16S rRNA gene were amplified, using universal primers S-D-Bact-0341-b-S-17 (forward 5′-CCTACGGGNGGCWGCAG-3′) and S-D-Bact-0785-a-A-21 (reverse 5′-GACTACHVGGGTATCTAATCC-3′) (44). Sequencing was performed by LifeSequencing S.L. (Valencia, Spain) on an Illumina MiSeq instrument (San Diego, California, USA).
An adaptation of the “Quantitative Insights Into Microbial Ecology” (QIIME) v1.9.0 package (16) was used to analyze the sequence data (26). Briefly, sequences were clustered into Operational Taxonomic Units (OTUs) based on 97% sequence identity as proxy for bacterial species using VSEARCH v2.03 with exclusion of chimeric sequences identified against the RDP gold database. Taxonomic assignment will be performed using the RDP classifier (27) against the SILVA119 database (28). The microbial species diversity (α-diversity) was calculated using Faith's phylogenetic diversity (PD) (29); and the Shannon index for diversity was calculated by applying correction for the differences in sequencing depths by rarefaction. Sequences are deposited in EMBL-EBI European Nucleotide Archive (Accession number: PRJEB39067).
Short-Chain Fatty Acid Analyses
SCFA analysis was performed as described by Bakker-Zierikzee et al. (30). Briefly, SCFAs were extracted from fecal PBS homogenate. The homogenates were mixed in formic acid, ethyl butyric acid (1.25 g/L; Sigma USA), and MilliQ water. Some fecal SCFAs such as acetic, propionic, n-butyric, iso-butyric, n-valeric, and iso-valeric acids were detected through gas chromatography with capillary column (Varian 3800, Varian Inc., Walnut Creek, CA, USA).
Statistics Analysis
All statistical analyses were performed using PRISM GraphPad 6 (version 6.02). Samples analyses were done with pairwise Mann–Whitney non-parametric comparison. The interaction of the intestinal microbiota composition in relation to geographical regions was visualized spatially using distance matrices on principal component analysis (PCA), which is an unconstrained ordination method, using CANOCO 5 (2012) software.
Results
Subject Background Information
From the 116 included participants recruited from public day care centers from four different geographical regions Santo André (n = 35) (Southeast region, SAN), Natal (n = 31) (Northeast region, NAT), Porto Alegre (n = 16) (Southern region, POA), and Uberaba (n = 17) (Southeast region, UBE), 99 met all requirements for subsequent analyses. Children in this cohort are evenly distributed in gender (Supplementary Table 1). Based on data collected from 2-day food recall questionnaires, we split each field into five categories: milk (inclusive of cow's milk/formula), vegetables, drinks (excluding water), animal proteins, and junk food (Figure 1). SAN has shown significant later age introduction of milk consumption (p < 0.001). All four cities show a high level of junk food introduction before a child is 1 year old, and POA was significantly later than the other cities (Figure 2).
Figure 1. Food diary average frequency consumed by child collected using 2-day food frequency questionnaire (FFQ) across the four cities (SAN, Santo Andre; UBE, Uberaba; POA, Porto Alegre; and NAT, Natal) showing that junk food consumption is prevalent in all four cites after the age of 1 year, while milk consumption varied.
Figure 2. Food groups based on 2-day food frequency questionnaire (FFQ) documented. Milk group includes cow's milk, breastmilk (if any), formula, and yogurt; drinks include sweetened juices and beverages; vegetables are inclusive of legumes; junk food includes sweets, biscuits, and chocolate. Differences in frequencies of complementary foods consumption, such as food puree and family diet, were not shown, as they were not statistically significant. Where *p < 0.05, **p < 0.005, ****p < 0.0001.
Microbiota Profiling in Brazilian Children
Initial processing and filtering of raw 16S rRNA gene sequence data resulted in minimum 25,546 high-quality reads per sample, and rarefaction analyses showed slopes reaching a plateau after a rapid initial increase, indicating a sufficient sequencing depth (Table 1). Diversity and richness of the fecal microbiota composition were the highest in participants from Uberaba. There are significant differences between Uberaba and Natal (p = 0.012) and Porto Alegre (p = 0.0056). Similar results were obtained by other diversity indices (Shannon, Simpson, and Chao1).
Table 1. Number of subjects and sequencing statistics by geographical region and Firmicutes/Bacteroidets ratio.
Taxonomic classification of high-quality sequence data showed that children aged between 1 and 3 years have microbiota dominated by taxa affiliated to the bacterial phyla Firmicutes and Bacteroidetes followed by Verrucomicrobia and Actinobacteria (Figure 3A). There were slight significant difference in seen in Actinobacteria & Proteobacteria in the overall community composition which might be due to inter-individual variations. We observed no difference between other subject characteristics such as birth mode (cesarean section or vaginal delivery) or antibiotic usage prior to the collection of samples (data not shown). Table 1 shows the average Firmicutes/Bacteroidetes ratio (F/B ratio) for each region ranging from 1.47 to 2.0 (Santo Andre, Natal, Porto Alegre, and Uberaba). F/B ratio has been shown by multiple studies to correlate with obesity and other metabolic diseases (31).
Figure 3. Relative abundance at the (A) phylum and (B) family level shows dominance of taxa affiliated to the phyla Firmicutes and Bacteroidetes (shown in gray and light blue, respectively).
Comparison of relative abundances at the bacterial family level revealed that the six most abundant bacterial families within the cohort were members of the phyla Firmicutes (Ruminococcaceae and Lachnospiraceae), Bacteroidetes (Bacteroidaceae and Prevotellaceae), Verrucomicrobium (Verrucomicrobiaceae), and Actinobacteria (Bifidobacteriaceae) (in order of relative abundance) (Figures 3, 4).
Figure 4. Comparison of various phyla across the four regions reveals significant differences in Actinobacteria and Proteobacteria. Where *p < 0.05, **p < 0.005.
PCA was used to explore the taxonomic compositions of fecal samples. No clear clustering by the geographic origin was observed (Figure 5), which suggests that microbiota composition was not primarily driven by biogeographic factors.
Figure 5. Unconstrained principal component analysis (PCA) of subjects grouped into four geographical regions: San Andre (SAN), Uberaba (UBE), Porto Alegre (POA), and Natal (NAT). Each dot represents one sample, and each color corresponds to a region. The PCA plot shows scattering of the samples across the quadrants. The arrows indicate the bacterial genera or families that potentially influence the scattering of the samples.
Microbiota Metabolite Functionality Indicators (Short-Chain Fatty Acid)
The major detected metabolites were acetic, butyric, and propionic acids. There were small amounts of other valeric acid and the iso-forms of butyric and valeric acids, which are significantly different. It is also noted that Uberaba has the highest level of SCFA than the other three regions, although the profiles are very similar (see Figure 6 and Table 2).
Figure 6. Comparison of short-chain fatty acids (SCFAs) across the geographical locations indicates that subject variation contributes to the significant differences observed in the different fatty acids. Where *p < 0.05, **p < 0.005, ***p < 0.0005.
Table 2. Average (minimum–maximum) fecal short chain fatty acid concentrations in mmol/kg of fecal material measured in subjects from the four sample sites.
Discussion
In this study, we have explored the microbiota composition of young children from four cities in Brazil located in three geographical regions. Although the locations are not near, the average microbiota composition was similar. Uberaba, based on demographics data (IBGE, Brazil 2019), has the most agriculture activity of the four regions. In our study, children from all centers are residing in urban areas and had access to different imported produce similarly for all day care centers. The Natal region is also characterized by exponential population growth. However, the degree of urbanization is less than that of the Southern parts of the country.
Bacterial diversity was found to be the greatest in children from Uberaba, indicated by highest alpha-diversity indices. The microbiota composition of these Brazilian children was dominated by Firmicutes and Bacteroidetes, which resembled more closely mature adult-like composition (31, 32). The F/B ratio has been extensively examined for human and mouse gut microbiota and has been linked to metabolic health by multiple studies utilizing this ratio to investigate obesity and age (31, 33). De Filippo et al. (23) compared fecal microbiota compositions of young children from a rural African village in Burkina Faso with fecal microbiota compositions of European children and proposed that children from less developed countries are characterized by lower F/B ratios (higher Bacteroidetes abundances), which may be an adaptation to a more polysaccharide-rich diet. The urbanized populations investigated in this study are also characterized by increased consumption of high-calorie food, which is low in nutritional value (junk food). This supports the hypotheses that energy-dense ultra-processed foods, a hallmark of the Western diet, indeed have pervasive effects on the microbiota development.
Other studies have reported positive correlations between F/B ratios and body mass index (BMI) and obesity (34, 35, 45). A study performed in Mexico reported that malnourished children have a higher F/B ratio than normal children of the same age. They have attributed this to the high amount of sugar and low-fiber diet that these children consumed (36). This may also be true for our study where children also consumed a lot of sweet juices and confectioneries (Figure 2). These habits might potentially increase children's predisposition to obesity in later life (25, 34, 37). In a further investigation, we found no significant interactions between BMI, dietary intake, and microbiota in our study. However, there are notable limitations to using only a 2-day recall questionnaire sampling diet in school and at home. These measurements might be biased due to caregiver's recollection and may not truly reflect a child's diet.
Interestingly, the Actinobacteria levels in the young children of this cohort are lower across all regions as compared with the abundances of this phylum observed in other microbiota studies. De Filippo et al. (23) reported high Bifidobacterium in European countries. In developed countries, Bifidobacterium has been documented to be an abundant and keystone member of the microbiota from infancy and all the way into toddlerhood (7, 9). Several Brazilian research groups have investigated Bifidobacterium spp. in various cities such as Porto Alegre (maternal-infant abundance) and Minas Gerais (1st months of life) (20, 38, 39). Another member of the phylum Actinobacteria, Collinsella was most abundant in children from Santo Andre. The genus Collinsella belongs to the family Coriobacteriaceae within the phylum Actinobacteria and has been positively correlated with insulin sensitivity and intestinal cholesterol absorption (35, 40). There is some evidence that Collinsella abundances depend on the dietary fiber intake (40, 41).
In line with the higher taxonomy diversity, the fecal SCFA profile also indicated that children from Uberaba have overall higher levels than children from the other three regions. This suggests that diet consumed by Uberaba children is richer in plant-based fibers, which facilitates microbial SCFA production. Fecal acetate, which is commonly associated with the activity of Bifidobacterium species, has been linked with decreased enteropathogenic infections by acidification of the gut (42). In our cohort, we observed low abundance of bifidobacterial species; this suggests the presence of other acetic acid-producing bacteria from the different genera. This suggests that the microbiota of these Brazilian children are different in this aspect than the children from the Amazonas of Venezuela, rural Malawi, and US metropolitan areas described Yatsunenko et al. (6), who report decreasing levels of Bifidobacterium with age.
A large-scale nutritional study across Brazil concluded that public schools have standardized food guidelines, which were the main factor driving food consumption to be similar (43). In the same study, household surveys showed a large population of children aged between 2 and 6 years who were not full time enrolled in the schools and required more careful assessment of their food intake and nutritional status. This study also highlighted that in several regions, consumption of processed sugary foods is high, even though nutrient levels are adequate.
To our knowledge, our study is one of the first to describe the gut microbiota composition, dietary intake, and function in young children across several cities in Brazil. We hypothesize that the microbiota of these Brazilian young children rapidly transited to a more adult-like microbiota profile unlike what has been reported for children of similar ages from either Western countries or “rural” areas in Africa and South America characterized by lower consumption of processed foods.
In conclusion, our study suggests that urban dietary factors play a larger role in influencing the gut microbiota composition, although they are separated geographically. It also illustrates the importance of considering the microbiota composition in diet, nutrition, and lifestyle recommendations and evaluating the impact of Westernization on health. This study is limited in that it only includes cross sections of these population at a single time point. A longitudinal sampling will be required to understand potential difference in microbiota of children from the same state with different socioeconomic status and to investigate how these environmental and lifestyle factors may contribute to later life health.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics Statement
The studies involving human participants were reviewed and approved by CEP-FMABC, Comitê de Ética, Santo Andre. Written informed consent to participate in this study was provided by the participants' legal guardian/next of kin.
Author Contributions
All authors have participated in the concept and design, analysis, and interpretation of data, drafting or revising of the manuscript. All authors have approved the manuscript and agree with the submission to Pediatrics.
Funding
The authors declare that this research project was funded by Danone Nutricia Brazil.
Conflict of Interest
The study design, data collection, analysis was carried out initially in Brazil with the participation of KB (full-time Danone Nutricia Brazil employee), MM, HM, ML, and VW. The study protocol is not directly related to any commercial product of the sponsor. CC (Danone Nutricia and Wageningen University), JK (Danone Nutricia and Wageningen University), and GR (Danone Nutricia) participated actively in the laboratory analysis of the intestinal microbiota and short chain fatty acids, interpretation of results, statistical analysis, and writing of the paper. MM, HM, ML, and VW declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We thank all the toddlers and their parents and caregivers for participating in this study. We also thank the staff of the early childhood public schools [Santo André (SA), Porto Alegre (POA), Uberaba (MG), and Natal (RN)] for their support.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fped.2020.573815/full#supplementary-material
References
1. Arrieta M-C, Stiemsma LT, Dimitriu PA, Thorson L, Russell S, Yurist-Doutsch S, et al. Early infancy microbial and metabolic alterations affect risk of childhood asthma. Sci Transl Med. (2015) 7:307ra152. doi: 10.1126/scitranslmed.aab2271
2. Hill CJ, Lynch DB, Murphy K, Ulaszewska M, Jeffery IB, O'Shea CA, et al. Evolution of gut microbiota composition from birth to 24 weeks in the INFANTMET Cohort. Microbiome. (2017) 5:4. doi: 10.1186/s40168-016-0213-y
3. Stiemsma LT, Michels KB. The role of the microbiome in the developmental origins of health and disease. Pediatrics. (2018) 141:e20172437. doi: 10.1542/peds.2017-2437
4. Palmer C, Bik E, DiGiulio D, Relman D, Brown P. Development of the human infant intestinal microbiota. PLoS Biol. (2007) 5:e0050177. doi: 10.1371/journal.pbio.0050177
5. Koenig JE, Spor A, Scalfone N, Fricker AD, Stombaugh J, Knight R, et al. Succession of microbial consortia in the developing infant gut microbiome. Proc Natl Acad Sci U S A. (2011) 108(Suppl 1):4578–85. doi: 10.1073/pnas.1000081107
6. Yatsunenko T, Rey FE, Manary MJ, Trehan I, Dominguez-Bello MG, Contreras M, et al. Human gut microbiome viewed across age and geography. Nature. (2012) 486:222–7. doi: 10.1038/nature11053
7. Agans R, Rigsbee L, Kenche H, Michail S, Khamis HJ, Paliy O. Distal gut microbiota of adolescent children is different from that of adults. FEMS Microbiol Ecol. (2011) 77:404–12. doi: 10.1111/j.1574-6941.2011.01120.x
8. Ringel-Kulka T, Cheng J, Ringel Y, Salojärvi J, Carroll I, Palva A, et al. Intestinal microbiota in healthy U.S. young children and adults-a high throughput microarray analysis. PLoS One. (2013) 8:e64315. doi: 10.1371/journal.pone.0064315
9. Hollister EB, Riehle K, Luna RA, Weidler EM, Rubio-Gonzales M, Mistretta T-A, et al. Structure and function of the healthy pre-adolescent pediatric gut microbiome. Microbiome. (2015) 3:1–13. doi: 10.1186/s40168-015-0101-x
10. Maggini S, Pierre A, Calder PC. Immune function and micronutrient requirements change over the life course. Nutrients. (2018) 10:1531. doi: 10.3390/nu10101531
11. Fidelis CMF, Osório MM. Consumo alimentar de macro e micronutrientes de crianças menores de cinco anos no Estado de Pernambuco, Brasil. Rev Bras Saúde Materno Infantil. (2007) 7:63–74. doi: 10.1590/S1519-38292007000100008
12. de Castro MA, Verly E, Fisberg M, Fisberg RM. Children's nutrient intake variability is affected by age and body weight status according to results from a Brazilian multicenter study. Nutr Res. (2014) 34:74–84. doi: 10.1016/j.nutres.2013.09.006
13. de Carvalho CA, de Almeida FPC, Priore SE, Franceschini SCC, de Novaes JF. Food consumption and nutritional adequacy in Brazilian children: a systematic review. Rev Paul Pediatr. (2015) 33:211–21. doi: 10.1016/j.rpped.2015.03.002
14. Hilger J, Goerig T, Hoeft B, Weber P, CostaCarvalho N, Goldberger U, et al. Micronutrient intake status in healthy toddlers: a multinational perspective: Kristina Hoffmann. Eur J Public Health. (2014) 24(Suppl 2):6939–55. doi: 10.1093/eurpub/cku165.085
15. Lopes WC, Marques F, Oliveira CF, Rodrigues JA, Silveira MF, Caldeira AP, et al. Infant feeding in the first two years of life. Rev Paul Pediatr. (2018) 36:164–70. doi: 10.1590/1984-0462/;2018;36;2;00004
16. Samuel BS, Shaito A, Motoike T, Rey FE, Backhed F, Manchester JK, et al. Effects of the gut microbiota on host adiposity are modulated by the short-chain fatty-acid binding G proteincoupled receptor, Gpr41. Proc Natl Acad Sci U S A. (2008) 105:1676716772. doi: 10.1073/pnas.0808567105
17. Canani RB, Di Costanzo M, Leone L. The epigenetic effects of butyrate: potential therapeutic implications for clinical practice. Clin Epigenet. (2012) 4:4. doi: 10.1186/1868-7083-4-4
18. Korpela K, Flint HJ, Johnstone AM, Lappi J, Poutanen K, Dewulf E, et al. Gut microbiota signatures predict host and microbiota responses to dietary interventions in obese individuals. PLoS One. (2014) 9:e0090702. doi: 10.1371/journal.pone.0090702
19. Taddei CR, Oliveira FF, Duarte RTD, Talarico ST, Takagi EH, Carvalho IR, et al. High abundance of Escherichia during the establishment of fecal microbiota in Brazilian children. Microb Ecol. (2014) 67:624–34. doi: 10.1007/s00248-014-0381-x
20. Grzeskowiak Ł, Teixeira TFS, Bigonha SM, Lobo G, Salminen S, Ferreira CL, et al. Gut Bifidobacterium microbiota in one-month-old Brazilian newborns. Anaerobe. (2015) 35(Part B):54–8. doi: 10.1016/j.anaerobe.2015.07.004
21. Mello CS, Barros KV, de Morais MB. Brazilian infant and preschool children feeding: literature review. J Pediatr (Rio J). (2016) 92:451–63. doi: 10.1016/j.jped.2016.02.013
22. Mueller NT, Shin H, Pizoni A, Werlang IC, Matte U, Goldani MZ, et al. Delivery mode and the transition of pioneering gut-microbiota structure, composition and predicted metabolic function. Genes. (2017) 8:364. doi: 10.3390/genes8120364
23. De Filippo C, Cavalieri D, Di Paola M, Ramazzotti M, Poullet JB, Massart S, et al. 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. (2010) 107:14691–6. doi: 10.1073/pnas.1005963107
24. De Filippo C, Paola MD, Ramazzotti M, Albanese D, Pieraccini G, Banci E, et al. Diet, environments, and gut microbiota. A preliminary investigation in children living in rural and urban Burkina Faso and Italy. Front Microbiol. (2017) 8:1979. doi: 10.3389/fmicb.2017.01979
25. Ayeni FA, Biagi E, Rampelli S, Fiori J, Soverini M, Audu HJ, et al. Infant and adult gut microbiome and metabolome in rural Bassa and Urban settlers from Nigeria. Cell Rep. (2018) 23:3056–67. doi: 10.1016/j.celrep.2018.05.018
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:335336. doi: 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:5264–7. doi: 10.1128/AEM.00062-07
28. Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. (2012) 41:D590–D6. doi: 10.1093/nar/gks1219
29. Faith DP, Baker AM. Phylogenetic diversity (PD) and biodiversity conservation: some bioinformatics challenges. Evol Bioinform Online. (2007) 2:121–8. doi: 10.1177/117693430600200007
30. Bakker-Zierikzee A, Alles M, Knol J, Kok F, Tolboom J, Bindels J. Prebiotic oligosaccharides in infant nutrition: addition of only galactooligosaccharides does not induce faecal acidic pH and SCFA-spectrum typical for breast fed infants: PN1-20. J Pediatr Gastroenterol Nutr. (2005) 40:693. doi: 10.1097/00005176-200505000-00235
31. Ley RE, Peterson DA, Gordon JI. Ecological and evolutionary forces shaping microbial diversity in the human intestine. Cell. (2006) 124:837–48. doi: 10.1016/j.cell.2006.02.017
32. Costello EK, Lauber CL, Hamady M, Fierer N, Gordon JI, Knight R. Bacterial community variation in human body habitats across space and time. Science. (2009) 326:16941697. doi: 10.1126/science.1177486
33. Mariat D, Firmesse O, Levenez F, Guimarăes VD, Sokol H, Doré J, et al. The Firmicutes/Bacteroidetes ratio of the human microbiota changes with age. BMC Microbiol. (2009) 9:123. doi: 10.1186/1471-2180-9-123
34. Koliada A, Syzenko G, Moseiko V, Budovska L, Puchkov K, Perederiy V, et al. Association between body mass index and Firmicutes/Bacteroidetes ratio in an adult Ukrainian population. BMC Microbiol. (2017) 17:120. doi: 10.1186/s12866-017-1027-1
35. Hou Y-P, He Q-Q, Ouyang H-M, Peng H-S, Wang Q, Li J, et al. Human gut microbiota associated with obesity in Chinese children and adolescents. BioMed Res Int. (2017) 2017:8. doi: 10.1155/2017/7585989
36. Méndez-Salazar EO, Ortiz-López MG, Granados-Silvestre MLÁ, Palacios-González B, Menjivar M. Altered gut microbiota and compositional changes in firmicutes and proteobacteria in mexican undernourished and obese children. Front Microbiol. (2018) 9:2693. doi: 10.3389/fmicb.2018.02693
37. Bergström A, Skov TH, Bahl MI, Roager HM, Christensen LB, Ejlerskov KT, et al. Establishment of intestinal microbiota during early life: a longitudinal, explorative study of a large cohort of Danish infants. Appl Environ Microbiol. (2014) 80:2889–900. doi: 10.1128/AEM.00342-14
38. Brandt K, Taddei CR, Takagi EH, Oliveira FF, Duarte RTD, Irino I, et al. Establishment of the bacterial fecal community during the first month of life in Brazilian newborns. Clinics. (2012) 67:113–23. doi: 10.6061/clinics/2012(02)05
39. Mueller NT, Shin H, Pizoni A, Werlang IC, Matte U, Goldani MZ, et al. Birth modedependent association between pre-pregnancy maternal weight status and the neonatal intestinal microbiome. Sci Rep. (2016) 6:23133. doi: 10.1038/srep23133
40. Gomez-Arango LF, Barrett HL, Wilkinson SA, Callaway LK, McIntyre HD, Morrison M, et al. Low dietary fiber intake increases Collinsella abundance in the gut microbiota of overweight and obese pregnant women. Gut Microbes. (2018) 9:189–201. doi: 10.1080/19490976.2017.1406584
41. Walker AW, Ince J, Duncan SH, Webster LM, Holtrop G, Ze X, et al. Dominant and dietresponsive groups of bacteria within the human colonic microbiota. ISME J. (2010) 5:220–30. doi: 10.1038/ismej.2010.118
42. Fukuda S, Toh H, Hase K, Oshima K, Nakanishi Y, Yoshimura K, et al. Bifidobacteria can protect from enteropathogenic infection through production of acetate. Nature. (2011) 469:543547. doi: 10.1038/nature09646
43. Bueno MB, Fisberg RM, Maximino P, Rodrigues Gde P, Fisberg M. Nutritional risk among Brazilian children 2 to 6 years old: a multicenter study. Nutrition. (2013) 29:405–10. doi: 10.1016/j.nut.2012.06.012
44. 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) 41:e1. doi: 10.1093/nar/gks808
Keywords: gut microbiome, gut microbiota, young children, toddler, preschool
Citation: Chew C, Barros KV, Weffort VRS, Maranhão HdS, Laranjeira M, Knol J, Roeselers G and de Morais MB (2020) Gut Microbiota of Young Children Living in Four Brazilian Cities. Front. Pediatr. 8:573815. doi: 10.3389/fped.2020.573815
Received: 26 June 2020; Accepted: 26 October 2020;
Published: 07 December 2020.
Edited by:
Pietro Vajro, University of Salerno, ItalyReviewed by:
Tudor Lucian Pop, Iuliu Haţieganu University of Medicine and Pharmacy, RomaniaClaudia Mandato, AORN Santobono-Pausilipon, Italy
Copyright © 2020 Chew, Barros, Weffort, Maranhão, Laranjeira, Knol, Roeselers and de Morais. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Mauro Batista de Morais, maurobmorais@gmail.com