- Research article
- Open Access
Lifestyle modifications result in alterations in the gut microbiota in obese children
BMC Microbiology volume 21, Article number: 10 (2021)
The association between the gut microbiota and pediatric obesity was analyzed in a cross-sectional study. A prospective study of obese children was conducted to assess the gut microbial alterations after a weight change. We collected fecal samples from obese children before and after a 2-month weight reduction program that consisted of individual counseling for nutritional education and physical activity, and we performed 16S rRNA gene amplicon sequencing using an Illumina MiSeq platform.
Thirty-six participants, aged 7 to 18 years, were classified into the fat loss (n = 17) and the fat gain (n = 19) groups according to the change in total body fat (%) after the intervention. The baseline analysis of the gut microbiota in the preintervention stages showed dysbiotic features of both groups compared with those of normal-weight children. In the fat loss group, significantly decreased proportions of Bacteroidetes phylum, Bacteroidia class, Bacteroidales order, Bacteroidaceae family, and Bacteroides genus, along with increased proportions of Firmicutes phylum, Clostridia class, and Clostridiales order, were observed after intervention. The microbial richness was significantly reduced, without a change in beta diversity in the fat loss group. The fat gain group showed significantly deceased proportions of Firmicutes phylum, Clostridia class, Clostridiales order, Lachnospiraceae family, and Eubacterium hallii group genus, without a change in diversity after the intervention. According to the functional metabolic analysis by the Phylogenetic Investigation of Communities by Reconstruction of Unobserved States 2, the “Nitrate Reduction VI” and “Aspartate Superpathway” pathways were predicted to increase significantly in the fat loss group. The cooccurring networks of genera were constructed and showed the different microbes that drove the changes between the pre- and postintervention stages in the fat loss and fat gain groups.
This study demonstrated that lifestyle modifications can impact the composition, richness, and predicted functional profiles of the gut microbiota in obese children after weight changes.
ClinicalTrials.govNCT03812497, registration date January 23, 2019, retrospectively registered.
The prevalence of pediatric obesity has increased over the past decade, leading to an increase in concomitant childhood health conditions, including type 2 diabetes mellitus, hypertension, dyslipidemia, fatty liver disease, and psychological problems . The imbalance between energy intake and expenditure is considered the most important cause of obesity . The cornerstone of pediatric obesity management is lifestyle intervention, including dietary modification and increased physical activity, and not weight loss medicines, calorie-restricted diets, or bariatric surgery, which are current treatments for obese adults [3, 4].
Several studies have shown that obesity is associated with gut microbial dysbiosis [5, 6]. The contribution of the microbiome to obesity has been considered using multifactorial approaches, such as supplying additional calories to the host, affecting satiety, favoring fat storage, and disrupting the integrity of the epithelial barrier . Recently, many studies, regarding weight loss interventions in adults, have investigated the relationship between the gut microbiota and obesity . Some reports have demonstrated that weight reduction by bariatric surgery partially reversed obesity-associated microbial alterations in obese adults [9, 10]. A randomized controlled trial involving adults showed that modification of the gut microbiota composition by probiotics could reduce body weight [11, 12]. Regarding childhood obesity, most previous studies of the gut microbiota have been cross-sectional in design [13, 14]. A recent study in obese children revealed that secreted proteins of the gut microbiota affect the microbial composition . One prospective study of children that investigated the microbial changes with weight gain over a four-year period showed that the microbiome-host-diet configuration could be a possible predictor of obesity . However, research on microbial changes over time after weight reduction interventions has remained scarce [13, 17].
In this study, we aimed to investigate changes in the composition, diversity, predicted functional metabolic profiles, and correlation networks of the gut microbiota in obese children after lifestyle modifications.
After screening, 42 obese children participated in the first intervention, and six obese children declined further participation. In this study, adiposity was defined as the measured total body fat percentage in body composition analysis. A total of 36 participants in all of the interventions were classified into two groups: the fat loss group (n = 17, 47.2%), including those who experienced a decrease in total body fat (%) after the intervention, and the fat gain group (n = 19, 52.8%), including those who experienced an increase in total body fat (%) after the intervention.
The mean ages of the children were 10.0 (SD: 2.4) years in the fat loss group and 10.3 (SD: 2.7) years in the fat gain group (t-test, P = 0.733), and 58% of participants were male in each group (chi-square test, P = 0.542, Table 1). Birthweight, delivery type, duration of the intervention, and numbers of exercises and nutritional counseling sessions were not significantly different between the fat loss and fat gain groups (Table 1). The questionnaire on general lifestyle and eating habits was developed for this study, and its contents are shown in Additional file 1. No significant differences were detected in the results of the questionnaires between the fat loss and fat gain groups (Table 1). As expected, the fat loss group showed significantly decreased total body fat (%), total body fat mass (kg), visceral fat area (cm2), and abdominal fat (%) in the body composition analysis, whereas the fat gain group showed significantly increased values after the interventions (paired t-test and Wilcoxon’s signed-rank test, P < 0.05, Table 2). Body mass index (BMI) was significantly decreased in the fat loss group and significantly increased in the fat gain group after lifestyle modifications (paired t-test, P < 0.05, Table 2). In the fat loss group, the level of alanine aminotransferase (ALT) was significantly decreased after the intervention (Wilcoxon’s signed-rank test, P = 0.022, Table 2). The insulin level was significantly decreased in the fat loss group; however, the homeostasis model assessment method-insulin resistance (HOMA-IR) levels were significantly increased in both groups (Wilcoxon’s signed-rank test, P < 0.05, Table 2). Comparison of the anthropometric measurements and blood biochemical profiles before the intervention showed no significant differences between the fat loss and fat gain groups (Table S1).
For the baseline analysis, 16S rRNA gene sequencing data from the feces of the obese groups at the preintervention stages were compared with those of 24 normal-weight children (18 boys and 6 girls, aged 8.1 ± 1.5 years old) from our previous cross-sectional study of pediatric obesity as controls using Quantitative Insights into Microbial Ecology 2 (QIIME2) . Among the results of anthropometric measurements, BMI in the obese group before intervention was significantly higher than that in the control group (Wilcoxon’s rank-sum test, P < 0.05, Table S2). A total of 8,221,270 sequences (mean of 139,343 sequences) were generated from 60 samples. After quality control, the dataset was reduced to a total of 6,074,850 sequences, with a mean of 102,963 sequences per sample, for 13,431 features. At the phylum level, Firmicutes and Bacteroidetes were dominant components of the gut microbiota, followed by Proteobacteria, Actinobacteria, and Verrucomicrobia in the fat loss, fat gain, and control groups (Fig. 1a). We compared the relative abundances of taxa between the control and the preintervention gut microbiota in the obese group using Statistical Analysis of Metagenomic Profiles (STAMP) software . The relative abundance of Bacteroidetes was significantly lower in the preintervention stage in the fat gain group than in the control group (Welch’s t-test, FDR = 0.000846, Fig. 1b). The compositional differences in the gut microbiota of the controls and the preintervention stages in both groups at each phylogenetic level are shown in Fig. 1b. At the genus level, in the preintervention stage, both groups showed increased relative abundances of Blautia, Dorea, Eubacterium hallii group, and Fusicatenibacter compared with the control group, and the fat gain group also showed decreased proportions of Bacteroides, Oscillibacter, and Parabacteroides (Welch’s t-test, FDR < 0.05, Fig. 1b, c). No significantly different taxa at any phylogenetic level were identified between the fat loss and fat gain groups at the preintervention stage. The alpha rarefaction plot of the observed_Operational Taxonomic Unit (OTU) s indices presents the richness of the samples in the control and the preintervention stages in both groups (Fig. S1a). The Shannon diversity index showed no significant differences between the controls and the preintervention stages in both groups (Fig. 1d). In the preintervention stage, both groups showed significantly lower observed_OTUs indices than the control group, and the fat gain group showed significantly lower observed_OTUs indices than the fat loss group (Wilcoxon’s rank-sum test, controls vs fat loss pre q = 0.009, controls vs fat gain pre q = 0.00006, fat loss pre vs fat gain pre q = 0.022, Fig. 1e). Principal coordinate analysis (PCoA) based on the weighted Unifrac distances between the genus-level microbial profiles showed a significant separation between the controls and obese individuals in both groups by permutational multivariate analysis of variance (PERMANOVA) (controls vs fat loss pre: q = 0.002, controls vs fat gain pre: q = 0.001, fat loss pre vs fat gain pre: q = 0.279, Fig. 1f).
The 16S rRNA sequence datasets were collected from feces of the participants before and after the lifestyle modifications. A total of 10,895,632 sequences (mean of 151,328 sequences) were generated from 72 samples. After quality control, the dataset was reduced to a total of 7,296,519 sequences, with a mean of 99,854 sequences per sample, for 9617 features. At the phylum level, Firmicutes and Bacteroidetes were major components of the gut microbiota, followed by Actinobacteria, Proteobacteria, and Verrucomicrobia, in both the pre- and postintervention stages in the fat loss and fat gain groups (Fig. 2a, b). To investigate the changes in the microbial community after intervention, we performed pairwise differential abundance comparison with Analysis of Differential Abundance Taking Sample Variation Into Account (ALDEx2) (Bioconductor v.3.11) . In the fat loss group, a total of 8 taxa (2 phyla, 2 classes, 2 orders, 1 family, and 1 genus) showed significantly different relative abundances after intervention (Wilcoxon signed-rank test, P < 0.05, Fig. 2c). The relative abundance of Firmicutes was significantly increased in the fat loss group (Wilcoxon’s signed-rank test, P = 0.009); conversely, the relative abundance of Bacteroidetes was significantly decreased in the fat loss group after intervention (Wilcoxon’s signed-rank test, P = 0.014, Fig. 2c). Bacteroidia class, Bacteroidales order, Bacteroidaceae family, and Bacteroides genus were significantly decreased, and Clostridiales order and Clostridia class were significantly increased in the fat loss group (Wilcoxon’s signed-rank test, P < 0.05, Fig. 2c). In the fat gain group, 6 taxa (2 phyla, 1 class, 1 order, 1 family, and 1 genus) showed significantly different relative abundances after intervention (Wilcoxon’s signed-rank test, P < 0.05, Fig. 2d). The relative abundance of Firmicutes was significantly decreased (Wilcoxon’s signed-rank test, P = 0.028) and the relative abundance of Actinobacteria was significantly increased in the fat gain group (Wilcoxon’s signed-rank test, P = 0.047, Fig. 2d). The proportions of Clostridia class, Clostridiales order, Lachnospiraceae family, and Eubacterium hallii group genus were significantly decreased in the fat gain group after intervention (Wilcoxon’s signed-rank test, P < 0.05, Fig. 2d).
The alpha rarefaction plot of the observed_OTUs indices clearly illustrates the richness of the samples in the pre- and postintervention stages in both groups (Fig. S1b). In the q2 longitudinal command in QIIME2, the alpha diversity indices, including the Shannon (Fig. 2e), Faith’s PD (Fig. 2f), and observed_OTUs (Fig. 2g) indices, were significantly changed in the fat loss group (Wilcoxon’s signed-rank test, P < 0.05) but not in the fat gain group. The degree of difference of the observed_OTUs indices after intervention was significantly higher in the fat loss group than in the fat gain group (Wilcoxon’s rank-sum test, FDR = 0.045, Fig. 2g). In the linear mixed-effect models, the observed_OTUs index was significantly decreased after lifestyle intervention in the fat loss group (coefficient: -7.5 [IQR: − 14.789, − 0.102], SD; 3.747, P = 0.047, Fig. 2h). The weighted Unifrac distance showed no significant differences between the pre- and postintervention stages in both groups (Fig. S2).
A total of 346 and 363 metabolic pathways in the MetaCyc database were predicted by Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt)2 in the fat loss and fat gain groups, respectively [21, 22]. The weighted nearest sequenced taxon index (NSTI) was used in the fat loss (mean ± SD: 0.127 ± 0.015) and fat gain (mean ± SD: 0.109 ± 0.010) groups . To identify significantly changed predicted pathways in the microbial community after intervention, we performed pairwise differential abundance comparisons with ALDEx2. At the FDR 0.5 level, the “Nitrate Reduction VI” and “Aspartate Superpathway” metabolic pathways were predicted to significantly increase after intervention (Wilcoxon’s signed-rank test, FDR = 0.013, and 0.030, respectively, Fig. 3a, b). In the fat gain group, no significantly changed predicted pathways were identified.
Analyses of the cooccurrence networks of genera in the fat loss and fat gain groups were conducted (Fig. 4a, b). Nodes represent the Amplicon Sequence Variant (ASV) s at the genus level, with size reflecting the relative abundance in the community and color indicating the phylum. The edge represents the correlation between connecting nodes, with edge thickness indicating the correlation value and green and red colors indicating positive and negative correlations, respectively. Only significant correlations (two-sided pseudo P < 0.05 based on permutations of 100 iterations) with correlation thresholds ≥ 0.3 are presented. The microbial dysbiosis index (MD index) was 0.6009 in the pre/postintervention fat loss group and 0.9251 in the pre/postintervention fat gain group . Regarding the network plot properties, the average path length, total nodes, and total edges decreased after the lifestyle intervention (Fig. 4c). On the basis of NetShift analysis, the genera Romboutsia, Ruminococcaeceae_UCG_013, Eubacterium coprostanollgenes-group, and Parabacteroides were identified as driver genera, with key roles in the changes in microbial interactions during the intervention in the fat loss group (Fig. 5a, c). We were also able to identify the genera Romboutsia, Eubacterium_halli_group, and Clostridium_sensu_stricto_1 as driving the microbial community changes from the pre- to the postintervention stages in the fat gain group (Fig. 5b, c).
The main treatment for obesity in children includes lifestyle modifications, such as dietary modification, increased physical activity, and behavioral modification, rather than medications or bariatric surgery . In recent decades, growing evidence has identified an association between the gut microbiota and obesity, and microbiota-targeted strategies have attracted attention in the context of obesity treatment . This study demonstrated that lifestyle modifications could exert a significant influence on the composition, diversity, and predicted functional profiles of the gut microbiota in childhood obesity.
In the baseline analysis, the obese children in the fat gain group showed decreased proportions of Bacteroidetes compared with the controls, as described by previous studies [15, 18, 27, 28]. The gut microbiota in the obese children showed a different microbial composition from that in the controls, including increased proportions of the genera Blautia, Dorea, Eubacterium hallii group, and Fusicatenibacter, which were reported to be associated with obesity, as well as decreased proportions of the genera Bacteroides, Oscillibacter, and Parabacteroides, which were positively associated with leanness [18, 29,30,31]. In both the preintervention groups, reduced microbial richness was noted, and the beta diversity showed significant separation between the control and obese groups, generally implicating dysbiotic features .
Recent studies have reported that the microbiome adapts quickly to lifestyle changes . In the fat loss group, the relative abundances of Firmicutes and Bacteroidetes were changed significantly after lifestyle modifications. Previous clinical weight reduction trials reported no consistent results regarding these major phyla. Some trials showed an increase in the relative abundance of Bacteroidetes along with a decrease in Firmicutes; others reported the opposite effects on these phyla, similar to the results of our study, while a few trials observed no impact [13, 33,34,35]. A population-based study of school-aged children showed that these inconsistent results in microbial communities are related to environmental factors, including socioeconomic status, age, and weight . The members of Bacteroidetes were reported to be more largely influenced by environmental factors rather than by host genetics, and some clinical trials demonstrated that these members are associated positively with a diet rich in protein and animal fat [16, 37, 38]. Our study demonstrated that Bacteroidia class, Bacteroidales order, Bacteroidaceae family, and Bacteroides genus, belonging to the Bacteroidetes phylum, were significantly decreased in the fat loss group. This decrease could be related to diet changes, including a decrease in highly processed foods rich in animal fat for 2 months. Otherwise, among the members of Firmicutes, increased proportions of Clostridiales order and Clostridia class in the fat loss group as well as decreased proportions of Clostridia class, Clostridiales order, Lachnospiraceae family, and Eubacterium hallii group genus in the fat gain were demonstrated after intervention. One possible explanation for these changes in the relative abundance of Firmicutes is that our lifestyle modification program did not include a carbohydrate-restrictive diet. A carbohydrate-restrictive diet can result in a reduction in the abundance of Firmicutes due to low intake of complex carbohydrates, which act as prebiotics . Further studies are required to determine whether the changes in microbiota composition were the result of lifestyle modifications.
In the fat loss group, the bacterial taxa richness was decreased after the application of lifestyle modifications. This observation is contrary to previous studies, which showed that a diet intervention, including calorie and carbohydrate restriction, increases in microbial gene richness with a decrease of adiposity in participants [40, 41]. However, some studies in adults have revealed that microbial richness decreased after a short-term dietary restriction, and other studies in pediatric populations showed no differences [42,43,44]. A recent study in Mexican children revealed significantly greater richness and diversity in the obese group than in the normal-weight group . In one study of 61 adults who underwent bariatric surgery for obesity treatment, low microbial gene richness was correlated with truncal fat mass and remained low at 1 year after surgery, suggesting a complex interplay between the gut microbiome and host obesity . Shoaie et al. analyzed the correlations of reduced bacterial gene counts and the production of several amino acids with the occurrence of metabolic diseases using a computational tool, and they showed that dietary interventions might reduce these products and improve insulin sensitivity . The microbial richness in obesity with weight changes shows an inconsistent pattern, suggesting that more studies are required.
The mechanism by which the microbiota affects energy balance in the human body are not clear. Our results from the PICRUSt2 analysis demonstrated two metabolic pathways in the MetaCyc database that were predicted to significantly increase after lifestyle modifications in the fat loss group. In the “Nitrate Reduction VI” pathway, the final product is L-glutamine, supplementation with which has recently been reported to lead to weight loss [46, 47]. The “Aspartate Superpathway” pathway includes both L-aspartate biosynthesis and degradation to oxaloacetate due to reversible transamination. Oxaloacetate was reported to activate brain mitochondrial biogenesis, leading to enhancement of the insulin pathway . However, these predictions based on the 16S rRNA gene amplicon sequencing data could be limited because the actual functions of the full metagenome likely differ, although amplicon-based predictions might be highly correlated with functional profiles based on shotgun metagenomics sequencing data . To increase the accuracy of predictions, we used the weighted NSTI value, which is a measure of how closely related the ASVs in each sample are to the reference genomes in the database . Analysis was conducted after excluding ASVs presenting weighted NSTI values of 0.15 or greater, indicating generally unreliable prediction . Much more research is required to determine the specific mechanisms associated with functional analysis.
We analyzed the microbial cooccurrence network with Sparse Correlation for Compositional data (SparCC) analysis, which is capable of avoiding the microbial compositional bias introduced when correlating relatively abundant data in Spearman’s and Pearson’s analyses . The MD index in SparCC is the log ratio of the total abundance of genera increased in the preintervention group to the total abundance of genera decreased in the postintervention group . In the previous study, this index showed a strong positive correlation with clinical disease severity and a negative correlation with species richness, resulting in an empirical estimation of the degree of dysbiosis within the microbiome . The fat gain group had a higher MD index value than the fat loss group, representing dysbiosis or an imbalance in the microbial community after intervention. In the network plot properties, the degrees of decrease in the average path length, total nodes, and total edges in the postintervention stage compared to those in the preintervention stage in the fat gain group were greater than those in the fat loss group. This outcome suggests that the relationships of the gut microbiota in the fat gain group decreased to be fewer than those in the fat loss group after intervention.
In the NetShift analysis, the driver genera that played key roles in microbial interactions were different between the fat loss and fat gain groups. The role of Parabacteroides as a driver of community change after lifestyle intervention was greater in the fat loss group than in the fat gain group. This finding suggested that Parabacteroides might play a key role in the process of weight loss during lifestyle modifications. Recent studies have revealed that Parabacteroides exerts an anti-obesogenic effect and could suppress the systemic inflammatory response by regulating IL-10 and Treg cells [31, 43]. Another driver taxon in the fat loss group was Ruminococcaceae_UCG_013, which was demonstrated to be associated with a decreased risk of weight gain, suggesting functional linkage to a lean phenotype in a large-scale, longitudinal adult study . Another driver genus in the fat loss group was Eubacterium_coprostanoligenes_group, which is known to reduce cholesterol by cholesterol-to-coprostanol conversion . The only common driver taxon between the fat loss and fat gain groups was Romboutsia, which was recently described . Romboutsia produces end products such as acetic acid, ethanol, iso-butanoic acid and iso-valeric acid, which are substrates involved in gluconeogenesis and lipogenesis . Other driver taxa in the fat gain group were Eubacterium hallii and Clostridium sensu stricto groups, which produce short-chain fatty acids [56, 57]. These results suggest that lifestyle modifications might exert a different effect on the interactions between microbial communities according to the direction of weight change.
In adult studies, changes in the gut microbiota were associated with weight reduction interventions, including restrictive diets, bariatric surgery, and medications including pre, pro-, and synbiotics and metformin; these interventions are difficult to apply in children due to the risk of nutritional imbalance and surgery . Our study showed BMI reduction in approximately 50% of participants; this result was relatively higher than those in adult studies involving weight loss interventions, which have generally resulted weight loss in 5–20% of participants [59, 60]. Our weight reduction program was a multidisciplinary individualized approach with frequent contact, which is the principal treatment of pediatric obesity, rather than restricted diet or surgery. Frequent management seemed to be effective in reducing the weight of the participants in our study. Most previous cross-sectional studies of childhood obesity have reported altered gut microbiota; however, few clinical trials have investigated gut microbial changes with weight reduction interventions [13, 14]. One study in preschool children participating in a behavioral intervention program reported no significant changes in microbial composition or functional profiles associated with weight loss . However, our study showed altered gut microbial composition, richness, and predicted functional profiles with weight loss. These inconsistent results could be related to the report that the gut microbiota varies with age, ethnicity and diet . Our study of Korean obese children could provide additional valuable information on common traits characterizing pediatric obesity.
Among the markers of the insulin resistance of the fat loss group, the fasting insulin level was improved after the intervention, although HOMA_IR was increased. One possible explanation is that the increase of HOMA-IR levels could be associated with unchanged fasting glucose and HbA1C levels, because our intervention included no carbohydrate restriction. Some studies have reported limitations of HOMA-IR in subjects with high fasting glucose levels and indicated that the fasting insulin value, rather than HOMA-IR, could be a surrogate measure of insulin resistance [61, 62].
Although the duration of the intervention, number of exercises, and nutritional counseling were not significantly different between the fat loss and fat gain groups, the responses to the intervention were different. In this regard, some considerations were noted in our study. In the baseline study, higher microbial richness was identified in the fat loss group than in the fat gain group. This observation agrees with other studies revealing that higher bacterial richness was associated with greater decreases in adiposity for obese adults with bariatric surgery and dietary interventions [10, 42]. The microbial richness at baseline might be a predictive potential factor for the efficacy of interventions. Further, a recent study revealed that baseline microbiota composition is not predictive of weight loss for the intervention . In line with this finding, our baseline study showed no significant differences in the microbial composition between the fat loss and fat gain groups. One limitation of our study was the small number of participants. Another limitation is that we cannot identify whether the gut microbiota will recover and how long it will take. To clarify the function and pathophysiology of the gut microbiota in relation to childhood obesity, further trials that include a larger number of obese children are necessary.
In conclusion, we observed significant alterations in the composition, richness, and expected functional profiles of the gut microbiota with weight loss after lifestyle modifications. Lifestyle modifications could impact microbiota dynamics, although little is known about the effect in obese children.
Participants, questionnaires, and anthropometric measurements
This longitudinal cohort study was an analysis of fecal samples collected from obese children before and after a 2-month weight reduction program. We recruited 50 obese children aged 7 to 18 years old at Hallym University Kangnam Sacred Heart Hospital from August 2018 to August 2019. Obesity was defined as a BMI ≥ 95th percentile based on the 2017 Korean growth chart . Those who had congenital heart disease, chronic inflammatory bowel disease, chronic liver disease, or chronic renal disease were excluded. The participants were required to take no antibiotics, probiotics, or steroids for 1 month prior to the intervention. If participants had acute inflammatory diseases, such as influenza, pneumonia, or acute gastroenteritis, we delayed the lifestyle modification program by 1 month. The participants completed questionnaires, which provided multiple choice questions on general lifestyle (the time spent studying and using electronic devices, the duration and frequency of regular exercise, the presence of easily accessible locations to exercise, and their mode of transportation to school) and eating habits (meal duration, the consumption of late-night snacks, the consumption of breakfast, and the intake of sugar-sweetened beverages) and submitted them at the first hospital visit. Anthropometric measurements, including height, weight, waist circumference, and blood pressure, were performed by professionally trained research nurses at the first and third hospital visits [64, 65]. Body composition analysis of total body fat mass, skeletal muscle mass, total body fat percent, visceral fat area, and abdominal fat percent was performed with an InBody 770 analyzer (Biospace Co. Ltd., Seoul, South Korea) at the first and third hospital visits.
Weight reduction programs
Over 8 weeks, the participants received weight reduction counseling from dietitians, exercise professionals, research nurses, and a pediatric clinician three times. At the first visit, the pediatric clinician, who specialized in childhood obesity, designed individualized feasible lifestyle modification programs in conjunction with dietitians, exercise professionals, and research nurses based on interviews with the participants and their guardians. Dietitians provided practical nutritional counseling on eating habits rather than a restrictive meal menu, and they suggested one or two individualized “must-follow” recommendations, such as eating breakfast, avoiding sugar-sweetened beverages, decreasing processed foods rich in animal fat or lengthening meal duration, to the participants at every counseling session. Exercise professionals found ways for obese children to exercise without a disruption to their general lives after analyzing the participants’ lifestyles, and they suggested one or two individualized “must-follow” recommendations, such as performing daily stretching exercises at home, cycling after school or on the weekends, using stairs instead of elevators if safe, and walking to school or academy, to the participants at every counseling session. The research nurse monitored the children to confirm that they were following the recommendations provided every 2 weeks. At the third visit, the pediatric clinician interviewed the participants and their guardians and suggested practical ways for participants to maintain their “must-follow” recommendations in general life.
Blood sampling and biochemical analysis
At the first and third visits, blood samples were obtained from the participants after an 8-h overnight fast. The levels of glucose, aspartate aminotransferase (AST), ALT, total cholesterol, triglycerides, high-density lipoprotein (HDL) cholesterol, low-density lipoprotein (LDL) cholesterol, high-sensitivity C-reactive protein (hs-CRP) and uric acid were measured using a Hitachi 7600 autoanalyzer (Hitachi, Tokyo, Japan). Concentrations of ferritin, insulin, and 25-OH vitamin D were determined using an ADVIA Centaur XP instrument (Siemens Diagnostics, Deerfield, IL, USA). The hemoglobin A1c (HbA1c) level was determined using a D-100 system (Bio-Rad Laboratories, Hercules, CA, USA). HOMA-IR was calculated as [insulin (μIU / mL) × glucose (mg / dL)]/ 405.
16S rRNA gene amplicon sequencing using an Illumina MiSeq platform and bioinformatics analysis
Fecal samples were collected before and after the interventions. Fecal samples were subsequently stored at − 80 °C within 1 h of sampling until DNA extraction. The genomic DNA in fecal samples was extracted using a FastDNA™ SPIN Kit for Soil (MP Biomedicals, Santa Ana, CA, USA) according to the manufacturer’s instructions. The DNA concentration was measured using a Quant-iT™ PicoGreen™ dsDNA Assay Kit (Invitrogen, Waltham, MA, USA). 16S rRNA gene amplicon sequencing of the V3–4 regions was performed by a commercial company (Chunlab Inc., Seoul, South Korea) using an Illumina MiSeq platform (Illumina, San Diego, CA, USA) . The sequencing outputs were generated as demultiplexed fastq sequences for downstream analysis using QIIME2, which identifies ASVs rather than OTUs. AVS approaches are generally considered to provide more precise identifications of microbes than OTU approaches . The dada2 denoise-paired command was used to filter the low-quality and chimeric sequences in the fastq reads. The taxonomy of these features was assigned via the Silva-genes reference database classifier (version 136) considering 99% similarity. 16S rRNA gene sequencing data of the preintervention stages in both groups were compared with those of 24 normal-weight children in our previous cross-sectional study of pediatric obesity as controls . The differential compositional analyses of the gut microbiota between the controls and both groups in the preintervention stages were performed by STAMP software . Alpha diversity (observed_OTUs index and Shannon index) and beta diversity (weighted Unifrac distance) between the control and the preintervention stages of the obese groups were analyzed at the sequence depth of 5041 in QIIME2 and were visualized by the ggplot2 package in R software. Pairwise significant differences in the relative abundances at phylogenetic levels between the pre- and postintervention stages of both groups were assessed by the ALDEx2 package and were visualized by the ggplot2 package in R software. The q2 longitudinal command in QIIME2 was used to analyze pairwise diversity changes after the intervention at the sequence depth of 5100. Linear mixed-effect models were constructed to analyze the relationships of the statistically significant alpha diversity indices. The PICRUSt2 (v2.3.0 beta) tool was used to infer the functional potentials of the gut microbiota on the basis of Enzyme Classification (EC) numbers in the MetaCyc database from 16S rRNA gene amplicon sequencing [21, 22]. The ASV table generated from QIIME2 was rarefied at the sequence depth of 5100 and then applied to PICRUSt2. After transforming the relative abundance, the ALDEx2 package was used to identify the significantly changed predicted functional pathways in our data after the intervention.
SparCC analysis at the genus level was conducted to detect coabundance and coexclusion correlations with two-sided pseudo P-values (P-values < 0.05 considered significant) based on 100 iterations in the MicrobiomeAnalyst web application . NetShift analyses were performed to detect the driver microbes in the web application based on the networks generated by SparCC .
Paired normally distributed data were analyzed using the paired t-test and are presented as the means and standard deviations; paired nonnormally distributed data were analyzed using the Wilcoxon’s signed-rank test and are presented as medians and interquartile ranges. The t-test and Wilcoxon’s rank-sum test were used to analyze the independent normally distributed data and the independent skewed continuous data, respectively. Categorical variables were analyzed by the chi-square test and are presented as frequencies and percentages. Statistical significance was declared at a P-value < 0.05. Resulting P-values were adjusted for multiple testing with the false discovery rate (FDR) method.
Availability of data and materials
All raw 16S rRNA gene sequencing data have been deposited in the NCBI Sequence Read Archive (SRA) under accession number SUB7459302 (BioProject PRJNA633584).
- F B ratio:
Firmicutes to bacteroidetes ratio
Operational taxonomic unit
Amplicon sequence variant
Quantitative insights into microbial ecology 2
Statistical analysis of metagenomic profiles
Analysis of differential abundance taking sample variation into account 2
Permutational multivariate analysis of variance
Phylogenetic investigation of groups by reconstruction of unobserved states 2
Sparse correlation for compositional data
Brown CL, Halvorson EE, Cohen GM, Lazorick S, Skelton JA. Addressing childhood obesity: opportunities for prevention. Pediatr Clin N Am. 2015;62:1241–61.
Romieu I, Dossus L, Barquera S, Blottière HM, Franks PW, Gunter M, et al. Energy balance and obesity: what are the main drivers? Cancer Causes Control. 2017;28:247–58.
Styne DM, Arslanian SA, Connor EL, Farooqi IS, Murad MH, Silverstein JH, et al. Pediatric obesity-assessment, treatment, and prevention: an endocrine society clinical practice guideline. J Clin Endocrinol Metab. 2017;102:709–57.
Kahan S, Manson JE. Obesity treatment, beyond the guidelines: practical suggestions for clinical practice. JAMA. 2019;321:1349–50.
Crovesy L, Masterson D, Rosado EL. Profile of the gut microbiota of adults with obesity: a systematic review. Eur J Clin Nutr. 2020;74:1251–62.
Le Chatelier E, Nielsen T, Qin J, Prifti E, Hildebrand F, Falony G, et al. Richness of human gut microbiome correlates with metabolic markers. Nature. 2013;500:541–6.
Aoun A, Darwish F, Hamod N. The influence of the gut microbiome on obesity in adults and the role of probiotics, prebiotics, and synbiotics for weight loss. Prev Nutr Food Sci. 2020;25:113–23.
Bryrup T, Thomsen CW, Kern T, Allin KH, Brandslund I, Jørgensen NR, et al. Metformin-induced changes of the gut microbiota in healthy young men: results of a non-blinded, one-armed intervention study. Diabetologia. 2019;62:1024–35.
Liu R, Hong J, Xu X, Feng Q, Zhang D, Gu Y, et al. Gut microbiome and serum metabolome alterations in obesity and after weight-loss intervention. Nat Med. 2017;23:859–68.
Aron-Wisnewsky J, Prifti E, Belda E, Ichou F, Kayser BD, Dao MC, et al. Major microbiota dysbiosis in severe obesity: fate after bariatric surgery. Gut. 2019;68:70–82.
Lin YC, Chen YT, Li KY, Chen MJ. Investigating the mechanistic differences of obesity-inducing Lactobacillus kefiranofaciens M1 and anti-obesity Lactobaillus mali APS1 by microbolomics and metabolomics. Front Microbiol. 2020;11:1454.
Cancello R, Turroni S, Rampelli S, Cattaldo S, Candela M, Cattani L, et al. Effect of short-term dietary intervention and probiotic mix supplementation on the gut microbiota of elderly obese women. Nutrients. 2019;11:3011.
Hollister EB, Foster BA, Dahdouli M, Ramirez J, Lai Z. Characterization of the stool microbiome in hispanic preschool children by weight status and time. Child Obes. 2018;14:122–30.
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:7585989.
Gallardo-Becerra L, Cornejo-Granados F, García-López R, Valdez-Lara A, Bikel S, Canizales-Quinteros S, et al. Metatranscriptomic analysis to define the Secrebiome, and 16S rRNA profiling of the gut microbiome in obesity and metabolic syndrome of Mexican children. Microb Cell Factories. 2020;19:61.
Rampelli S, Guenther K, Turroni S, Wolters M, Veidebaum T, Kourides Y, et al. Pre-obese children's dysbiotic gut microbiome and unhealthy diets may predict the development of obesity. Commun Biol. 2018;1:222.
López-Contreras B, Morán-Ramos S, Villarruel-Vázquez R, Macías-Kauffer L, Villamil-Ramírez H, León-Mimila P, et al. Composition of gut microbiota in obese and normal-weight Mexican school-age children and its association with metabolic traits. Pediatr Obes. 2018;13:381–8.
Shin S, Cho KY. Altered gut microbiota and shift in Bacteroidetes between young obese and normal-weight Korean children: a cross-sectional observational study. Biomed Res Int. 2020;2020:6587136.
Parks DH, Tyson GW, Hugenholtz P, Beiko RG. STAMP: statistical analysis of taxonomic and functional profiles. Bioinformatics. 2014;30:3123–4.
Fernandes AD, Reid JN, Macklaim JM, McMurrough TA, Edgell DR, Gloor GB. Unifying the analysis of high-throughput sequencing datasets: characterizing RNA-seq, 16S rRNA gene sequencing and selective growth experiments by compositional data analysis. Microbiome. 2014;2:15.
Douglas GM, Maffei VJ, Zaneveld J, Yurgel SN, Brown JR, Taylor CM, et al. PICRUSt2: an improved and extensible approach for metagenome inference. BioRxiv. 2019;20:672295. https://doi.org/10.1101/672295.
Caspi R, Billington R, Fulcher CA, Keseler IM, Kothari A, Krummenacker M, et al. The MetaCyc database of metabolic pathways and enzymes. Nucleic Acids Res. 2018;46:D633–9.
Langille MG, Zaneveld J, Caporaso JG, McDonald D, Knights D, Reyes JA, et al. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat Biotechnol. 2013;31:814–21.
Gevers D, Kugathasan S, Denson LA, Vázquez-Baeza Y, Van Treuren W, Ren B, et al. The treatment-naive microbiome in new-onset Crohn’s disease. Cell Host Microbe. 2014;15:382–92.
Paccosi S, Cresci B, Pala L, Rotella CM, Parenti A. Obesity therapy: how and why? Curr Med Chem. 2020;27:174–86.
Muscogiuri G, Cantone E, Cassarano S, Tuccinardi D, Barrea L, Savastano S, et al. Gut microbiota: a new path to treat obesity. Int J Obes Suppl. 2019;9:10–9.
Ley RE, Turnbaugh PJ, Klein S, Gordon JI. Microbial ecology: human gut microbes associated with obesity. Nature. 2006;444:1022–3.
Indiani CMDSP, Rizzardi KF, Castelo PM, Ferraz LFC, Darrieux M, Parisotto TM. Childhood obesity and Firmicutes/Bacteroidetes ratio in the gut microbiota: a systematic review. Child Obes. 2018;14:501–9.
de la Cuesta-Zuluaga J, Corrales-Agudelo V, Velásquez-Mejía EP, Carmona JA, Abad JM, Escobar JS. Gut microbiota is associated with obesity and cardiometabolic disease in a population in the midst of westernization. Sci Rep. 2018;8:11356.
Tseng CH, Wu CY. The gut microbiome in obesity. J Formos Med Assoc. 2019;118:S3–9.
Wang K, Liao M, Zhou N, Bao L, Ma K, Zheng Z, et al. Parabacteroides distasonis alleviates obesity and metabolic dysfunctions via production of succinate and secondary bile acids. Cell Rep. 2019;26:222–35.
Fragiadakis GK, Wastyk HC, Robinson JL, Sonnenburg ED, Sonnenburg JL, Gardner CD. Long-term dietary intervention reveals resilience of the gut microbiota despite changes in diet and weight. Am J Clin Nutr. 2020;111:1127–36.
Remely M, Tesar I, Hippe B, Gnauer S, Rust P, Haslberger A. Gut microbiota composition correlates with changes in body fat content due to weight loss. Benef Microbes. 2015;6:431–9.
Remely M, Hippe B, Zanner J, Aumueller E, Brath H, Haslberger AG. Gut microbiota of obese, type 2 diabetic individuals is enriched in Faecalibacterium prausnitzii, Akkermansia muciniphila and Peptostreptococcus anaerobius after weight loss. Endocr Metab Immune Disord Drug Targets. 2016;16:99–106.
Haro C, Montes-Borrego M, Rangel-Zúñiga OA, Alcalá-Díaz JF, Gómez-Delgado F, Pérez-Martínez P, et al. Two healthy diets modulate gut microbial community improving insulin sensitivity in a human obese population. J Clin Endocrinol Metab. 2016;101:233–42.
Moran-Ramos S, Lopez-Contreras BE, Villarruel-Vazquez R, Ocampo-Medina E, Macias-Kauffer L, Martinez-Medina JN, et al. Environmental and intrinsic factors shaping gut microbiota composition and diversity and its relation to metabolic health in children and early adolescents: a population-based study. Gut Microbes. 2020;11:900–17.
Johnson EL, Heaver SL, Walters WA, Ley RE. Microbiome and metabolic disease: revisiting the bacterial phylum Bacteroidetes. J Mol Med. 2017;95:1–8.
David LA, Maurice CF, Carmody RN, Gootenberg DB, Button JE, Wolfe BE, et al. Diet rapidly and reproducibly alters the human gut microbiome. Nature. 2014;505:559–63.
Conlon MA, Bird AR. The impact of diet and lifestyle on gut microbiota and human health. Nutrients. 2015;7:17–44.
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.
Zhernakova A, Kurilshikov A, Bonder MJ, Tigchelaar EF, Schirmer M, Vatanen T, et al. Population-based metagenomics analysis reveals markers for gut microbiome composition and diversity. Science. 2016;352:565–9.
Cotillard A, Kennedy SP, Kong LC, Prifti E, Pons N, Le Chatelier E, et al. Dietary intervention impact on gut microbial gene richness. Nature. 2013;500:585–8.
Zeng Q, Li D, He Y, Li Y, Yang Z, Zhao X, et al. Discrepant gut microbiota markers for the classification of obesity-related metabolic abnormalities. Sci Rep. 2019;9:13424.
Maya-Lucas O, Murugesan S, Nirmalkar K, Alcaraz LD, Hoyo-Vadillo C, Pizano-Zárate ML, et al. The gut microbiome of Mexican children affected by obesity. Anaerobe. 2019;55:11–23.
Shoaie S, Ghaffari P, Kovatcheva-Datchary P, Mardinoglu A, Sen P, Pujos-Guillot E, et al. Quantifying diet-induced metabolic changes of the human gut microbiome. Cell Metab. 2015;22:320–31.
Abboud KY, Reis SK, Martelli ME, Zordão OP, Tannihão F, de Souza AZZ, et al. Oral glutamine supplementation reduces obesity, pro-inflammatory markers, and improves insulin sensitivity in DIO wistar rats and reduces waist circumference in overweight and obese humans. Nutrients. 2019;11:536.
Palomo-Buitrago ME, Sabater-Masdeu M, Moreno-Navarrete JM, Caballano-Infantes E, Arnoriaga-Rodríguez M, Coll C, et al. Glutamate interactions with obesity, insulin resistance, cognition and gut microbiota composition. Acta Diabetol. 2019;56:569–79.
Wilkins HM, Harris JL, Carl SM, E L, Lu J, Eva Selfridge J, et al. Oxaloacetate activates brain mitochondrial biogenesis, enhances the insulin pathway, reduces inflammation and stimulates neurogenesis. Hum Mol Genet. 2014;23:6528–41. https://doi.org/10.1093/hmg/ddu371.
Douglas GM, Maffei VJ, Zaneveld JR, Yurgel SN, Brown JR, Taylor CM, et al. PICRUSt2 for prediction of metagenome functions. Nat Biotechnol. 2020;38:685–8.
Wang Y, Wang K, Huang L, Dong P, Wang S, Chen H, et al. Fine-scale succession patterns and assembly mechanisms of bacterial community of Litopenaeus vannamei larvae across the developmental cycle. Microbiome. 2020;8:106.
Friedman J, Alm EJ. Inferring correlation networks from genomic survey data. PLoS Comput Biol. 2012;8:e1002687.
Menni C, Jackson MA, Pallister T, Steves CJ, Spector TD, Valdes AM. Gut microbiome diversity and high-fibre intake are related to lower long-term weight gain. Int J Obes. 2017;41:1099–105.
Kenny DJ, Plichta DR, Shungin D, Koppel N, Hall AB, Fu B, et al. Cholesterol metabolism by uncultured human gut bacteria influences host cholesterol level. Cell Host Microbe. 2020;28:245–57.
Ricaboni D, Mailhe M, Khelaifia S, Raoult D, Million M. Romboutsia timonensis, a new species isolated from human gut. New Microbes New Infect. 2016;12:6–7.
Gerritsen J, Hornung B, Ritari J, Paulin L, Rijkers GT, Schaap PJ, et al. A comparative and functional genomics analysis of the genus Romboutsia provides insight into adaptation to an intestinal lifestyle. BioRxiv. 2019;845511.
Engels C, Ruscheweyh H-J, Beerenwinkel N, Lacroix C, Schwab C. The common gut microbe Eubacterium hallii also contributes to intestinal propionate formation. Front Microbiol. 2016;7:713.
Atarashi K, Tanoue T, Oshima K, Suda W, Nagano Y, Nishikawa H, et al. Treg induction by a rationally selected mixture of clostridia strains from the human microbiota. Nature. 2013;500:232–6.
Seganfredo F, Blume C, Moehlecke M, Giongo A, Casagrande D, Spolidoro J, et al. Weight-loss interventions and gut microbiota changes in overweight and obese patients: a systematic review. Obes Rev. 2017;18:832–51.
Louis S, Tappu R-M, Damms-Machado A, Huson DH, Bischoff SC. Characterization of the gut microbial community of obese patients following a weight-loss intervention using whole metagenome shotgun sequencing. PLoS One. 2016;11:e0149564.
Candela M, Biagi E, Soverini M, Consolandi C, Quercia S, Severgnini M, et al. Modulation of gut microbiota dysbioses in type 2 diabetic patients by macrobiotic Ma-pi 2 diet. Br J Nutr. 2016;116:80–93.
Lorenzo C, Haffner SM, Stancáková A, Laakso M. Relation of direct and surrogate measures of insulin resistance to cardiovascular risk factors in nondiabetic finnish offspring of type 2 diabetic individuals. J Clin Endocrinol Metab. 2010;95:5082–90.
Kang ES, Yun YS, Park SW, Kim HJ, Ahn CW, Song YD, et al. Limitation of the validity of the homeostasis model assessment as an index of insulin resistance in Korea. Metabolism. 2005;54:206–11.
Kim JH, Yun S, Hwang S-S, Shim JO, Chae HW, Lee YJ, et al. The 2017 Korean National Growth Charts for children and adolescents: development, improvement, and prospects. Korean J Pediatr. 2018;61:135–49.
Nichols EK, Nichols JS, Selwyn BJ, Coello-Gomez C, Parkerson GR, Brown EL, et al. Implications of the WHO child growth standards in rural Honduras. Public Health Nutr. 2012;15:1015–22.
Kollias A, Ntineri A, Kyriakoulis KG, Stambolliu E, Lagou S, Boubouchairopoulou N, et al. Validation of the professional device for blood pressure measurement microlife WatchBP office in adults and children according to the American National Standards Institute/Association for the Advancement of medical instrumentation/International Organization for Standardization standard. Blood Press Monit. 2018;23:112–4.
Yang J, Park J, Park S, Baek I, Chun J. Introducing murine microbiome database (MMDB): a curated database with taxonomic profiling of the healthy mouse gastrointestinal microbiome. Microorganisms. 2019;7:480.
Callahan BH, McMurdie PJ, Rosen MJ, Han AW, Johnson AJ, Holmes SP. DADA2: high resolution sample inference from Illumina amplicon data. Nat Methods. 2016;13:581–3.
Chong J, Liu P, Zhou G, Xia J. Using MicrobiomeAnalyst for comprehensive statistical, functional, and meta-analysis of microbiome data. Nat Protoc. 2020;15:799–821.
Kuntal BK, Chandrakar P, Sadhu S, Mande SS. ‘NetShift’: a methodology for understanding ‘driver microbes’ from healthy and disease microbiome datasets. ISME J. 2019;13:442–54.
The author has the greatest gratitude for all of the participants in the study.
The research presented in this publication was supported by the Basic Science Research Program through the National Research Foundation of Korea funded by the Ministry of Education (NRF-2018R1D1A1B07046799) and by a Hallym University Research Fund (HURF-2017-81), Anyang si, Republic of Korea. The funding bodies played no role in the study design; the collection, analysis or interpretation of the data; or the writing of the manuscript.
Ethics approval and consent to participate
This study was carried out in accordance with the Declaration of Helsinki and approved by the Hallym University Kangnam Sacred Heart Hospital Institutional Review Board under IRB No. 2018–06–026-003. Written informed consent was obtained from all of the participants and their parents on behalf of any participants under the age of 16 years old prior to the start of the study.
Consent for publication
The author declares no financial or any other competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1.
Questionnaire on general lifestyle and eating habits.
Additional file 2: Table S1
. Comparison of the Anthropometric Measurements and Blood Biochemical Findings between the Fat Loss and Fat Gain Groups in the Pre- and Postintervention Stages. Table S2. The Characteristics of the Normal Weight and the Preitervention Stages in the Fat Loss and Fat Gain Groups.
Additional file 3: Fig. S1.
The alpha rarefaction plots for the observed_OTUs index in the baseline study, including the controls and preintervention stages of the fat loss and fat gain groups (a), and the intervention study, including the pre- and postintervention stages in the fat loss and fat gain groups (b). Fig. S2. Principal coordinate analysis based on the weighted Unifrac distance between the pre- and postintervention stages in the fat loss (a) and fat gain groups (b).
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Cho, K.Y. Lifestyle modifications result in alterations in the gut microbiota in obese children. BMC Microbiol 21, 10 (2021). https://doi.org/10.1186/s12866-020-02002-3
- 16S rRNA gene
- Gut microbiota
- Weight reduction programs