Altered Gut Microbiota and Shift in Bacteroidetes between Young Obese and Normal-Weight Korean Children: A Cross-Sectional Observational Study

Emerging data suggest that the gut microbiome is related to the pathophysiology of obesity. This study is aimed at characterizing the gut microbiota composition between obese and normal-weight Korean children aged 5-13. We collected fecal samples from 22 obese and 24 normal-weight children and performed 16S rRNA gene sequencing using the Illumina MiSeq platform. The relative abundance of the phylum Bacteroidetes was lower in the obese group than in the normal-weight group and showed a significant negative correlation with BMI z-score. Linear discriminative analysis (LDA) coupled with effect size measurement (LEfSe) analysis also revealed that the Bacteroidetes population drove the divergence between the groups. There was no difference in alpha diversity, but beta diversity was significantly different between the normal-weight and obese groups. The gut microbial community was linked to BMI z-score; blood biomarkers associated with inflammation and metabolic syndrome; and dietary intakes of niacin, sodium, vitamin B6, and fat. The gut microbiota of the obese group showed more clustering of genera than that of the normal-weight group. Phylogenetic investigation of communities by reconstruction of unobserved states (PICRUSt) analysis revealed that the functions related to carbohydrate and lipid metabolism in the microbiota were more enriched in the normal-weight group than in the obese group. Our data may contribute to the understanding of the gut microbial structure of young Korean children in relation to obesity. These findings suggest that Bacteroidetes may be a potential therapeutic target in pediatric obesity.


Introduction
Childhood obesity is a major public health concern worldwide [1]. Obese children have a high risk of developing adult obesity and obesity-related comorbidities, including type 2 diabetes mellitus, cardiovascular disease, and psychological problems [2,3]. Recently, accumulating evidence suggests that the human gut microbiota is associated with many chronic diseases, including obesity. Previous studies in adults have demonstrated that obese individuals have different microbial compositions than lean individuals [4,5]. The gut microbiota is closely related to energy harvest and metabolism in humans [6,7]. Therefore, microbiota-targeted strat-egies have attracted much attention in the context of obesity treatment [8,9].
Studies on the obesity-related microbiota in children are still scarce. Gut microbial composition is known to vary with age, ethnicity, and diet [6,10]. Although previous studies have suggested that the gut microbiota of infants is converted into adult-like composition in the first 1-3 years [11], some evidence indicates that the microbiota continues to change until adolescence [12]. Therefore, valuable information on the obesity-related microbiota could be obtained from different ethnic and age groups.
The objective of this study was to characterize the composition of the gut microbiota among obese and normal-weight Korean children. Amplicons of the 16S rRNA gene were sequenced using Illumina MiSeq to analyze the composition of the gut microbiota. The dietary and lifestyle patterns of the participants and the levels of blood biochemical markers related to inflammation and metabolic disease were measured and examined in relation to the gut microbiota composition. The correlation structure of the microbiota was shown using network analysis, and predictive functional differences between groups were identified by phylogenetic investigation of communities by reconstruction of unobserved states (PICRUSt) analysis.

Subjects, Questionnaire, and Anthropometric
Measurements. Forty-six children who were 5-13 years of age were enrolled in the study at Hallym University Kangnam Sacred Heart Hospital from December 2017 to March 2018 (NCT03388411). This study protocol was approved by the Institutional Review Board of the Hallym University Kangnam Sacred Heart Hospital (2017-09-015). Written informed consent was obtained from all subjects of the study and their parents, in accordance with the Declaration of Helsinki. Based on the 2017 Korean growth chart [13], subjects with a body mass index ðBMIÞ ≥ 95th percentile were classified as the obese group, and subjects with 5th percentile < BMI < 85th percentile were classified as the normal-weight group. For a month prior to the stool sampling, none of the subjects had taken antibiotics/probiotics/steroids or had diarrhea. None of the subjects had acute infections or chronic disease. Subjects completed questionnaires on lifestyle, bowel habits, and dietary intake and submitted them at the hospital visit. Anthropometric measurements, including height, weight, waist circumference, midarm circumference, hip circumference, thigh circumference, and blood pressure, were performed by professionally trained personnel [14,15]. Body composition analysis including skeletal muscle mass and total body fat content was measured using the inBody 770 analyzer (Biospace Co. Ltd., Seoul, Korea).

Dietary
Assessments. The children and their parents received specific training from a dietitian to describe in a proper way all the foods and the quantities consumed, including the name/brand of the consumed food, recipes of dishes, method of preparation or cooking, and portion sizes. After training with the dietitian, the participants filled out everything the subjects ate and drank for 3 days: 2 weekdays and 1 weekend day. Using data from the dietary records, nutrient intakes were calculated by a dietitian using the Computer-Aided Analysis Program 4.0 for professionals (CAN-pro 4.0, Korean Nutrition Society, Seoul, Korea).

Stool Sampling, Bacterial DNA Extraction, Illumina
MiSeq Sequencing, and Bioinformatics. The stool samples were collected in sterile containers and immediately frozen at -80°C until DNA extraction. DNA was extracted using a QIAamp DNA Stool Mini Kit (Qiagen, Valencia, CA, USA) according to the manufacturer's instructions. Using 2 μL of the extracted DNA, polymerase chain reaction (PCR) amplifications were performed with the primers targeting the V3 to V4 regions of the 16S rRNA gene [16]. The products were then amplified by the second PCR with index primers. Equal concentrations of amplicons were pooled together and purified using an AMPure bead kit (Agencourt Bioscience, Beverly, MA, USA). The product size and quality were assessed on an Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA). Sequencing was carried out at ChunLab, Inc. (Seoul, Korea) using the Illumina MiSeq platform (Illumina, San Diego, CA, USA). The raw reads were checked for quality, and low-quality reads (average score < 25) were filtered by a Trimmomatic tool (version 0.32). Then, the paired-end sequences were merged using a PANDASeq Paired-end Assembler [17]. Chimeric sequences were removed using the UCHIME algorithm [18]. The taxonomic classification of each read was performed based on the EzBioCloud database (http://eztaxon-e.ezbiocloud.net) [19]. Sequences that corresponded to the reference sequence with greater than 97% similarity in EzBioCloud were considered to be identified at the species level. To compare the operational taxonomic units (OTUs) between samples, we determined shared OTUs through the EzBioCloud program (Chunlab Inc.). The functional potential of the microbiota was inferred using PICRUSt metagenomics prediction [20] and was categorized into levels 1-3 based on Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways [21].
2.5. Statistical Analysis. All statistical analyses were carried out using R software (version 3.5.2, http://www.r-project .org/). Based on the Shapiro-Wilk normality test, data are presented as the means and standard deviations (for continuous variables with normal distribution) or medians and interquartile ranges (for continuous variables with skewed distribution). Categorical variables are expressed as frequencies and percentages. In analyzing the characteristics, nutrient intake, and blood biochemical profiles of the participants, the t-test or the Kruskal-Wallis Rank Sum test was used according to the results of the normality test. For comparison of the gut microbiota between groups, the Mann-Whitney U test was used for continuous variables, and chi-square was used for categorical variables. Correlations between continuous variables were calculated using the Pearson correlation test. All analyses were performed after normalizing for the copy number of the bacterial 16S rRNA gene. The ACE index, the number of observed OTUs, the Chao1 richness estimate, and the Jackknife estimate were used to compare gut microbiota richness between samples. The within sample (alpha) diversity was compared using the Simpson diversity index, the Shannon index, and phylogenetic diversity. Cluster analysis was performed with nonmetric multidimensional scaling (NMDS), after computing the Bray-Curtis dissimilarity between each pair of individuals. The Fast UniFrac analysis was used to calculate the (beta) diversity between groups and was visualized with a principal coordinate analysis (PCoA). Differences in beta diversity between the normal-weight and obese groups were tested with nonparametric analysis of variance based on 999 permutations (permutational multivariate analysis of variance (PERMANOVA)). The differently abundant bacterial taxa between the normal-weight and obese groups were identified using the linear discriminant analysis (LDA) coupled with effect size measurement (LEfSe) method [22]. To analyze the ability of specific taxa to predict obesity, we calculated the area under the receiver operating characteristic curve (AUC of ROC). Multivariate analysis for relationships between gut microbial community composition, BMI z -score, and blood biochemical markers was performed using canonical correspondence analysis (CCA). Using the R package qgraph [23], we performed a network analysis to identify the difference in the correlation network of the gut microbiota between the normal-weight and obese groups. To obtain the relative robustness, the sample coverage threshold for the identified genera was set at ≥0.5. The resulting p values were adjusted for multiple testing with the false discovery rate (FDR) method [24]. p values ≤ 0.05 were considered statistically significant.

Results
3.1. Participant Characteristics, Lifestyle Questionnaire, and Blood Biochemical Marker. A total of forty-six children were enrolled in this study (obese group: n = 22; normal-weight group: n = 24). The characteristics of the participants are summarized in Table 1. The sex distribution showed no difference between the normal-weight and obese groups (p = 0:853; Table 1). The individuals in the normal-weight group were slightly younger than those in the obese group (p = 0:048; Table 1). All anthropometric measurements, including height (p = 0:002), weight (p < 0:001), BMI (p < 0:001), BMI z-score (p < 0:001), waist circumference (p < 0:001), midarm circumference (p < 0:001), hip circum-ference (p < 0:001), and thigh circumference (p < 0:001), showed significantly higher values in the obese group than in the normal-weight group (Table 1). Systolic/diastolic blood pressure (p = 0:002/p < 0:001), skeletal muscle mass (p < 0:001), and total body fat content (p < 0:001) also showed significantly higher values in the obese group than in the normal-weight group ( Table 1). The percentage of children born by cesarean section was significantly higher in the obese group (p = 0:012; Table 1). Questionnaires on lifestyle patterns differed between the normal-weight and obese groups. The percentage of children who did not exercise at all was 36.4% in the obese group and 0% in the normal-weight group (p = 0:004; Table 1). The percentage of children who exercised for 30 minutes or more per day was significantly higher in the normal-weight group than in the obese group (p = 0:025; Table 1). The percentage of children who watched television or used electronic devices (video games, smart phones, or computers) more than two hours a day was significantly higher in the obese group than in the normal-weight group (p = 0:017; Table 1). The percentage of children who used electronic devices near bedtime tended to be higher in the obese group than in the normal-weight group (p = 0:072; Table 1).

Gut Microbial Composition and Diversity.
After filtering low-quality, nontarget, and chimeric amplicons, 16S rRNA gene sequencing resulted in a total of 3.9 million highquality reads from 46 fecal samples. The median sequencing read was 78,155 (quartiles: 69,929; 85,731). We obtained a median value of 412.5 OTUs per sample (316; 529.3) after excluding low-abundance OTUs (<1% of total).
We analyzed the correlation between BMI z-score and the relative abundance of bacterial taxa ( Table 5). The relative abundance of the phylum Bacteroidetes was negatively correlated with BMI z-score (Pearson's correlation coefficient r = −0:329, p = 0:026; Figure 1(d)). Furthermore, BMI z-score was negatively correlated with the abundance of Bacteroidia at the class level; Bacteroidales at the order level; and Bacteroidaceae, Devosia_f, Leptotrichiaceae, Odoribacteraceae, Porphyromonadaceae, Rikenellaceae, and Staphylococcaceae at the family level ( Table 5). The relative abundance of the family Lachnospiraceae was positively correlated with BMI z-score (p = 0:05; Table 5). The relative abundance of the phylum Firmicutes was not significantly correlated with BMI z-score (r = 0:24, p = 0:1).
Gut microbiota richness measures showed no obvious difference between samples from normal-weight and obese children (ACE, p = 0:725; the number of observed OTUs, p = 0:692; the Chao1 richness estimate, p = 0:676; and the Jackknife estimate, p = 0:775). In addition, several alpha diversity estimates, including the Simpson diversity index

BioMed Research International
(p = 0:262), the Shannon index (p = 0:272), and phylogenetic diversity (p = 0:982), were not significantly different between normal-weight and obese children. NMDS analysis revealed separation and clustering of the obese group from the normal-weight group along the NMDS1 axis, while naïve tended to cluster along NMDS2 (Figure 1(g)). We used the Fast UniFrac analysis to measure beta diversity. The PCoA plot of the microbiota from all individuals in the normalweight and obese groups is shown in Figure 1(h). The beta diversity showed a statistically significant difference between normal-weight and obese children at the genus level (p = 0:009, PERMANOVA on Fast UniFrac distances; Figure 1(h)).

Taxonomic Differences in the Microbiota between Obese
and Normal-Weight Children. To identify the specific microbial profile distinguishing obese and normal-weight children, a metagenomics biomarker discovery approach using the LEfSe method was applied to assess the effect size of each differently abundant taxon ( Figure 2). The LDA effect size values are shown in Figure 2(f). Using the LEfSe method, we found that Bacteroidetes at the phylum level; Bacteroidia at the class level; Bacteroidales at the order level; Bacteroidaceae, Porphyromonadaceae, and Rikenellaceae at the family level; and Bacteroides, EF404788_g, Desulfovibrio_g3, Anaerofilum, Alistipes, Bacteroidaceae_uc, Hydrogenoanaerobacterium, EF402988_g, Oscillibacter, and Citrobacter at the     Figure 2). This population is dominated by bacteria belonging to the Bacteroidetes phyla in the normalweight group. In addition, Actinomyces, Romboutsia, Weissella, and GL872355_g at the genus level were significantly enriched in the obese group (p < 0:05; Figure 2).

Relationship between Gut Microbial Community
Composition, BMI z-Score, Blood Biochemical Markers, and Dietary Intake. We generated a correlogram to visualize the degree of association between BMI z-score, major phyla, blood biochemical markers, and dietary intake of energy and nutrients (Figure 3(a)). Variables with highly significant differences between the two groups (p value ≤ 0.1) were selected from the univariate analysis. Variables with a correlation coefficient (r) of 0.8 or greater were regarded as the same variables, and then one representative variable was selected from the same variables. The BMI z-score showed a significant positive correlation with inflammatory markers, including HOMA-IR; neutrophil count; and serum levels of triglycerides, hs-CRP, and ferritin, and with increased dietary intake of calories, fat, niacin, vitamin B6, P, Na, and zinc ( Figure 3(a)). On the other hand, serum vitamin D levels and the proportion of Bacteroidetes were negatively correlated with BMI z-score (Figure 3(a)). The Actinobacteria population showed negative correlations with the Bacteroidetes and Proteobacteria populations and positive correlations with blood markers including hs-CRP and neutrophil count and with dietary intake of calories and fat (Figure 3(a)). The correlogram showed that Bacteroidetes and Firmicutes exhibited different orientations in regard to correlation with most of the variables, including blood biomarkers and dietary intakes (Figure 3(a)). The Bacteroidetes population showed a negative correlation with the Firmicutes population and a positive correlation with the Proteobacteria population (Figure 3(a)). Moreover, the proportion of Bacteroidetes was negatively correlated with inflammatory markers, including hs-CRP, ferritin, HOMA-IR, and neutrophil count, and with dietary intake of calories, fat, niacin, vitamin B6, P, Na, and zinc (Figure 3(a)). On the other hand, the Firmicutes population showed a negative correlation with the proportion of Proteobacteria and positive correlations with the inflammatory markers hs-CRP and neutrophil count and with dietary intake of calories, fat, niacin, vitamin B6, and Na (Figure 3(a)).
We performed CCA to visualize the relationship between gut microbial community composition, BMI z-scores, and blood biochemical markers (Figure 3(b)) or dietary intake (Figure 3(c)). Variables were selected from the same standards with the correlogram. The distance between two points shows the significance of the correlation. The distance between Bacteroidetes and the microbial community in the normal-weight group is shorter than that between Firmicutes and the microbial community in the normal-weight group, suggesting that Bacteroidetes has a strong correlation with the normal-weight group (Figures 3(b) and 3(c)). The length of the blue line is proportional to the degree of importance. As shown in Figure 3(b), BMI z-score, HOMA-IR, hs-CRP, ferritin, and neutrophil count are more important in the microbial community of the obese group, whereas vitamin D is an important factor in the microbial community of the normal-weight group. Figure 3(c) shows that fat, Na, and niacin among dietary components, in addition to BMI z-score, are more important in the microbial community of the obese group, whereas vitamin B6 intake is an important factor in the microbial community of the normal-weight group. Firmicutes was negatively correlated with serum vitamin D levels and dietary intake of vitamin B6 (Figures 3(b) and 3(c)). Bacteroidetes was negatively correlated with BMI z -score, serum ferritin level, and fat intake (Figures 3(b) and  3(c)). Actinobacteria showed a negative correlation with HOMA-IR and dietary intake of Zn, P, niacin, and Na (Figures 3(b) and 3(c)). Proteobacteria was negatively correlated with neutrophil count and hs-CRP (Figures 3(b) and 3(c)).

Correlation Network.
We performed a correlation network analysis to investigate whether obesity was associated with alterations in the overall correlation structure of the gut microbiota (Figure 4 and Supplementary Table S1). Constructed networks revealed that samples from the normal-weight group had fewer edges, a lower mean degree, and a longer mean distance than those from the obese group, which indicates that there were fewer significant correlations and less clustering of genera (Supplementary Table S2). The betweenness centrality was higher in the obese group, which indicates that only a few genes were highly correlated in a network (Supplementary  Table S2). Bacteroidetes showed higher positive intraphylum correlations in the normal-weight group, and Firmicutes showed higher positive intraphylum correlations in the obese group (Supplementary Table S2).

Discussion
The present study showed differences in gut microbial composition between young normal-weight and obese Korean children aged 5-13 years. Obese children showed a significant reduction in Bacteroidetes, an elevated F : B ratio, and significantly different beta diversity compared with the same parameters among normal-weight children, as described by previous studies [4,25,26]. The Bacteroidetes population was also detected by LEfSe with a high LDA score, suggesting that it is the key phylotype responsible for the differences between the normal-weight and obese groups. The relative abundance of Firmicutes, however, revealed no significant difference between the groups. The results of our study suggest the importance of Bacteroidetes in pediatric obesity. Recent evidence has indicated that Bacteroidetes is a potentially modifiable therapeutic target because it is more largely influenced by environmental factors rather than host genetics [4,27,28]. In the future, prospective intervention studies will be needed to explore the impact of the specific species or strains belonging to the Bacteroidetes phylum on pediatric obesity modulation. Family Lachnospiraceae was significantly correlated with BMI z-score. This result is consistent with a previous experimental result showing that the colonization of bacteria belonging to Lachnospiraceae induces the development of diabetes in germ-free ob/ob mice [29]. From these findings, it could be assumed that Lachnospiraceae is involved in the development of metabolic dysfunction in children. Akkermansia is a mucin-degrading bacterium, and its abundance has been reported to be negatively correlated with obesity in previous studies with adults [30]. In one study, Akkermansia was reduced in obese children aged 4-5 years living in Sweden, which was analyzed by quantitative PCR [31]. However, the current reports including our study, using 16S rRNA next generation sequencing analysis, revealed no significant difference in Akkermansia levels between normal-weight and obese children [32,33]. This disparity can be explained with the differences in the methodology, ethnicity, and extent of its colonization which starts from early childhood and reaches a similar level to adults [34].
To investigate the relationships between gut microbial community composition, BMI, and selected variables from the biochemical markers and diet intake, CCA analysis was performed. The gut microbial community in the obesity group revealed a strong correlation with BMI z-score, which was in line with previous reports [35,36]. Inflammatory markers, including hs-CRP, neutrophil count, and ferritin, were related to microbial composition in the obese group, suggesting that obesity is closely linked to inflammation [37]. Evidence for a relationship between inflammation and the microbiota continues to be revealed. Bacterial products, such as lipopolysaccharide and short-chain fatty acids (SCFAs), can induce inflammation through immune cell activation and fat accumulation in adipocytes [38,39]. These findings suggest the role of the gut microbiota in the development of inflammation in the pathogenesis of pediatric obesity. Among dietary intakes, niacin, Na, and fat seemed to affect gut microbial composition in the obese group. Higher fat and Na intake is associated with obesity and metabolic syndrome [40,41]. Moreover, a recent study indicated the possible association of chronic niacin overload on pediatric obesity [42]. In our study, dietary intake of vitamin B6 seemed to be important in the microbial community of normal-weight children. Prior work has suggested that the gut microbiota of lean adolescents seems to be more involved in vitamin B6 synthesis [25]. Further research is needed to corroborate the present results. The correlation network of the gut microbiota in this study showed that the normalweight group had less clustering of genera than the obese group. This finding is consistent with prior studies by Riva et al. [32], showing that the gut microbiota in the obese group has a different correlation network structure than the gut microbiota in the normal-weight group.   the PICRUSt analysis showed that gut microbial function in obese children involves energy metabolism, such as photosynthesis and nitrogen metabolism, which can stimulate lipogenesis or gluconeogenesis [43]. Recent research has revealed that Bacteroides ovatus can release monosaccharides from cellulose and hemic cellulose for further metabolism by a wide variety of gut commensals via glycolytic pathways [44]. Our results showed that carbohydrate metabolism was more predicted in normal-weight children than in obese children, which can be speculated to be related to the positive association of Bacteriodes ovatus in normal-weight children.
Our functional analysis also showed a significantly greater presentation of genes involved in amino acid metabolism, such as cysteine and methionine metabolism and tyrosine and phenylalanine biosynthesis, in obese children than in normal-weight children. The fermentation pathways of cysteine and methionine are included in sulfur metabolism, which leads to the production of hydrogen sulfide, which has been known to have detrimental effects on colonic epithelial energy metabolism [45,46]. Tyrosine and phenylalanine biosynthesis are known to be associated with obesity, diabetes, and metabolic syndrome by reducing the activation of alkaline phosphatase [47,48]. This result in our study may have important long-term implications for bowel health in the context of the consumption of excessive protein diets. In addition, higher abundances of microbial communities related to lipid metabolism were observed in normal-weight children than in obese children [49]. These changes in predicted metabolic pathways caused by intestinal microbiota can induce an imbalance between energy production and absorption. The specific mechanism associated with this functional analysis will need to be studied further. In our study, hs-CRP, NLR, and MPV were higher in the obese group than in the normal-weight group. hs-CRP and NLR are well-known inflammatory markers that are associated with obesity because adipose tissue can be the major source of proinflammatory cytokines [50]. Increased MPV, a biomarker of platelet activity, is known to be associated with acute myocardial infarction, stroke, and thrombosis in individuals with morbid obesity [50]. The percentage of children born by cesarean section was significantly higher in the obese group than in the normal-weight group. This finding is the same as those of previous studies, which are explained by the disruption of mother-to-child transmission of gut microbiota associated with cesarean section [52,53]. In regard to lifestyle, the use of electronic devices is becoming a major problem in the context of pediatric obesity rather than a lack of physical activity [54]. The present study also indicated that  Figure 5: Distribution of tier 1 KEGG functional categories from the normal-weight (green) and obese (red) groups (a). Significant differences in the relative abundances of second (b) and third (c) tier KEGG functional categories between the normal-weight (orange) and obese (green) groups. Significance was determined using the Kruskal-Wallis test with FDR correction. Only functional categories with corrected p < 0:05 are shown. 16 BioMed Research International not only the lack of exercise but also the use of electronic devices for more than two hours a day was significantly higher in the obesity group. A recent study reported that electronic device usage close to bedtime can disrupt sleep patterns, which could lead to obesity [55]. Obesity is not only microbiota-driven; thus, a careful evaluation of all factors, including delivery mode, diet, and lifestyle, should be taken into account [56].
One of the strengths of this study is that it included only young Korean children, which makes the study population less influenced by ethnicity and environmental factors such as smoking, drugs, and alcohol use. Instead of simply comparing the microbiota between groups, we analyzed the relationship of multivariate factors in the microbiota, biochemical markers, and diet intake through a correlogram and CCA. Additionally, we showed the results of a network analysis and functional analysis of the gut microbiota in obese and normal-weight young children. The limitations of this study had to do with the small number of participants and the cross-sectional design, which prevents the determination of causality. Therefore, a prospective large-scale study is required to clarify the relationship between the microbiota and childhood obesity.
In conclusion, the microbial communities of obese children exhibited significant differences in beta diversity and a significantly elevated F : B ratio compared to those in normal-weight children. The phylum Bacteroidetes was significantly reduced in the obese group and was negatively correlated with BMI z-score. The LEfSe biomarker discovery analysis also suggested that the Bacteroidetes population was the key phylotype differentiating the two groups. These findings suggest the importance of Bacteroidetes in pediatric obesity. The gut microbial community in the obese group was linked to BMI z-score; blood biomarkers associated with inflammation and metabolic syndrome; and intakes of niacin, Na, and fat. In the microbial network analysis, the gut microbiota in the obese group showed more clustering of genera than the gut microbiota in the normal-weight group. PICRUSt analysis revealed that the functions related to carbohydrate and lipid metabolism were more enriched in the microbiota of the normal-weight group than in that of the obese group. Our data may contribute to the understanding of the gut microbial structure of young Korean children in relation to obesity. Further studies are required to target Bacteroidetes as a new therapeutic intervention for pediatric obesity.

Data Availability
All raw 16S rRNA gene sequencing data were deposited in the NCBI Sequence Read Archive (SRA) under accession number SUB7560699 (BioProject PRJNA637782). Data related to the current study are available from the corresponding author on reasonable request.

Disclosure
The manuscript was presented at the 69th Fall Congress of Korean Pediatrics and Youth Science Society. The funding bodies had no role in the study design; in the collection, analysis, or interpretation of the data; or in the writing of the manuscript.

Conflicts of Interest
The authors declare no conflicts of interest.