- 1Research Center for High Altitude Medicine, Qinghai University School of Medicine, Xining, China
- 2Department of Public Health, Qinghai University School of Medicine, Xining, China
- 3Key Laboratory for Application of High Altitude Medicine in Qinghai Province, Qinghai University, Xining, China
- 4General Surgery Department, Qinghai Provincial People's Hospital, Xining, China
- 5Qinghai-Utah Joint Research Key Lab for High Altitude Medicine, Qinghai University School of Medicine, Xining, China
Objective: To explore the relationship between intestinal flora and obesity in Tibetan children at different altitudes.
Methods: Using16S rRNA gene sequencing results and blood lipid metabolism indexes to study the characteristics of the intestinal flora present in faeces and changes in blood lipid metabolism in obese children in Tibet who reside at different altitudes and to study correlations between blood lipid metabolism indicators and the intestinal flora composition.
Results: The results showed the following. (a) The triglyceride (TG) and low-density lipoprotein cholesterol (LDL-C) levels in the obesity groups were higher than those in the normal-weight groups, and those in the high-altitude obesity groups were lower than those in the low-altitude obesity groups. (b) The 16S rRNA gene sequencing results showed that altitude affected the composition and relative abundance of the gut microbiota. These parameters were basically the same among the low-altitude groups, while they were significantly lower in the high-altitude groups than in the low-altitude groups. (c) Groups that lived at different altitudes and had different body weights had different dominant bacterial genera. Megamonas was closely related to obesity, and its relative abundance in the low-altitude groups was higher than that in the high-altitude groups. Prevotella was associated with altitude, and its relative abundance in the high-altitude groups was higher than that in the low-altitude groups. In addition, Prevotella elicited changes in the abundance of Escherichia-Shigella. The lower prevalence of obesity and incidence of intestinal inflammation in those living at high altitudes were related to the abundance of Prevotella. (d) There were correlations between the gut microbiota composition and lipid metabolism indicators. The abundance of Romboutsia was positively correlated with TG and LDL-C levels but negatively correlated with high-density lipoprotein cholesterol (HDL-C) levels. The abundance of Akkermansia was negatively correlated with LDL-C levels, and the abundance of Blautia was negatively correlated with body mass index (BMI) and LDL-C levels.
Conclusions: The intestinal flora diversity varied by body weight and altitude, with lower diversity in those at higher altitudes and with lower body weights. Prevotella likely plays a role in suppressing obesity at high altitudes.
Introduction
Obesity is a global health problem. Obesity is a complex, multifactorial and largely preventable disease that affects approximately one-third of people worldwide (1). Obesity is linked to type 2 diabetes (2), respiratory disease (3), cardiovascular disease (4, 5), female infertility (6, 7), and nonalcoholic fatty liver disease (7), and it can also be a risk factor for cancer (8). With improvements in living standards and changes in diets, childhood obesity has become increasingly prevalent worldwide, and children are becoming affected by obesity at a younger age (9, 10). Jia et al. (11) found that the prevalence of overweight and obesity in 6- to 17-year-old children increased from 11.7 to 25.2% and from 2.8 to 10.1% from 1991 to 2011, respectively. Additionally, the prevalence of obesity in children and adolescents in China is rapidly increasing; it increased at a rate of 0.36% annually from 2000 to 2011 (12). Many studies have shown that obesity in childhood is associated with an increased risk of obesity in adulthood (13, 14). This correlation suggests that attention should be given to the problem of childhood obesity.
Obesity is related to not only the imbalance of energy intake but also genetics, the environment and lifestyle; however, the underlying causes of obesity are still unclear (15). Among the nongenetic factors associated with obesity, intestinal flora composition has been recognized as a regulator of obesity, and a correlation between changes in the intestinal flora composition and body weight has been observed in a number of studies conducted in animal models of obesity (16–18). However, the status of the intestinal microbiome in children with obesity has not been well studied. Several studies (19–21) have revealed that the gut microbiome composition might be a factor that affects obesity. Changes in the biological diversity of the intestinal flora can affect the health and growth of the host to a large extent, but the specific factor that regulates the gut microbiome remains poorly understood.
The Tibetan Plateau is the highest plateau in the world, with a mean elevation of 4,500 m (22). It is characterized by unique environmental conditions, such as low pressure, low temperatures, low humidity and high radiation (23). A high-altitude hypoxic environment causes disorder of the intestinal flora; hypoxia causes inflammation, and inflammation aggravates hypoxia in tissues. In a high-altitude hypoxic environment, the intestinal mucosal barrier is highly prone to damage. This damage causes the translocation of intestinal bacteria and toxins, which causes immune system activation and the release of inflammatory mediators, resulting in disruption of the intestinal flora composition. Some studies have shown that factors related to altitude and anoxic environments have significant effects on the composition of the gut microbiome (24), and high altitude may contribute to shaping the human gut microbiota composition (25). The study by Santos (26) showed that the rate of overweight/obesity in children who reside at high altitudes is lower than that in children who reside at low altitudes. The results of our research group's previous epidemiological survey showed that the prevalence of overweight in 1,561 children aged 7–12 years who resided on the Qinghai-Tibet Plateau was 8.5% (9.9% for boys and 7.1% for girls), and the prevalence of obesity was 6.3% (6.0% for boys and 6.5% for girls); these values were much lower than those observed in the rest of China (11). Therefore, we sought to examine whether altitude affects the gut microbiota composition and thus influences obesity.
Materials and methods
Ethics statement
The experimental protocol was established according to the ethical guidelines of the Declaration of Helsinki and approved by the ethics committee of Qing Hai Provincial People's Hospital, Xining, China.
Study design
A cross-sectional study was designed to explore the association between high altitude residency and obesity with respect to lipid levels and gut microbiota diversity by comparing two indigenous communities who live at different altitudes. The study was approved by the Human Rights Committee of the University of Buenos Aires. Each parent provided written informed consent after the study was explained and before it was initiated.
Participants
We enrolled 70 subjects in our research study. All the participants were Tibetan children aged 7–12 years with the same eating habits. Using simple random sampling, we selected one primary school as a research site in each altitude area. Based on the China Center for Disease Control and Prevention (CDC) growth charts (27), we divided the subjects into four groups: obese children (N = 15, high-altitude obese children, HOB) and normal-weight children (N = 20, high-altitude normal-weight children, HN) who lived at an altitude of 4,500 m in Ma Duo County (34°52′12″N, 98°15′36″E), Guoluo Tibetan Autonomous Prefecture, Qinghai–Tibet Plateau; and obese children (N = 15, low-altitude obese children, LOB) and normal-weight children (N = 20, low-altitude normal-weight children, LN) who lived at an altitude of 1,500 m in Hai Dong County (36°18′N, 102°48′E), Hai Dong city, Qinghai–Tibet Plateau. A detailed questionnaire that included questions about age, sex, parental education level, history of breastfeeding, and family history of obesity was administered.
Inclusion and exclusion criteria
The inclusion criteria were as follows: (a) age 7–12 years old (28) and living in the research area for at least 5 years; (b) attendance at the selected boarding school; and (c) a diet containing grains, beef, and fresh vegetables. The exclusion criteria were as follows: (a) use of antibiotics or probiotics within the past 4 weeks; (b) gastrointestinal diseases (e.g., diarrhea) within the past 4 weeks; (c) frequent constipation; or (d) failure to agree to participate in this study. All the enrolled subjects were healthy, with no history of gastrointestinal disease, liver disease, hypertension, or diabetes, as demonstrated by their medical histories and physical examinations.
Anthropometric measurements
Height, weight, waist circumference and hip circumference were measured according to the protocol of the International Society for the Advancement of Kinanthropometry (29). Height was measured to the nearest 0.1 cm while the child stood upright against a mounted stadiometer with bare feet. In addition to height, waist and hip circumferences were measured. A digital weighing scale (Huawei3Pro, ShengZhen, China) that was calibrated regularly to the nearest 0.1 kg (after every 10 measurements) was used to measure body mass. We asked children to wear only a t-shirt and pants. Based on the height and body mass measurements, body mass index [BMI (weight/height2)] was calculated using the following formula:
Blood sample collection and analysis
Venous blood samples (5 ml) were collected by venepuncture into Vacutainer tubes, and plasma levels of glucose (GLU), triglycerides (TGs), total cholesterol (TC), high-density lipoprotein cholesterol (HDL-C) and low-density lipoprotein cholesterol (LDL-C) were analyzed with an LW-C400 automatic biochemical analyser (Landwind, Shenzhen, China).
Stool samples
Using sterile fecal boxes, fresh fecal samples (5 g) were collected from obese and control children within 2 h of defecation. Fecal samples were quickly placed in an ultralow temperature freezer (−80°C) for further processing and testing.
Determination of SCFAs contents in feces
Fresh fecal samples (2 g) were mixed with 3 ml deionized water and homogenized. The mixture was left for 20 min at room temperature and centrifuged at 14,000 rpm for 15 min. Then, the supernatant was transferred to new EP tubes. Then, 3 ml deionized water was added to the precipitate, and all of the above operations were repeated again. The extraction was repeated twice, and the supernatants were combined. Next, it was filtered through a 0.22 μm filter. A 1 ml sample solution was injected into the GC system. Chromatographic analysis of fecal samples was performed using ISQ Single Quadrupole GC–MSMS (Thermo Fisher Scientific), Trace 1300 Series Gas Chromatograph (Thermo Fisher Scientific) and TG-5MS Separation column (30 m × 0.32 mm, 0.25 μm). The injection port temperature was 160°C. The GC oven was programmed with an increasing starting temperature from 50°C for 1 min to 220°C with 10°C/min for 5 min. A volume of 0.2 μl of the sample was injected. Electron impact ionization (EI) was used as an ionization source for the GC/MS analysis at 70 eV. The injection temperature was 230°C, and the transfer line and ion source were set to 250°C. XcaliburTM software (Thermo Fisher Scientific) was used for the automatisation of the GC–MS system and for data acquisition. Sample quantification was obtained by means of acetic, propionic, and butyric acid standard curves (Sinopharm). Construction of Standard Curve The standard calibration curves were constructed using seven concentrations: 9.325, 18.75, 37.5, 75, 150, 300, and 900 μg/ml. The concentration of each sample was obtained by the standard curve. The standard curves of acetic acid, propionic acid, and butyric acid were Y = 3.796e5X + 9.363e4, Y = 3.936e5X + 4.003e5, and Y = 1.366e6X + 1.151e6, respectively (30, 31).
DNA extraction
Total genomic DNA was extracted from the samples using the acetyltrimethylammonium bromide (CTAB) method (28) DNA concentration and purity were analyzed using 1% agarose gels. According to the concentration, the DNA samples were diluted to 1 ng/μl using sterile water. An equal volume of 1 × loading buffer (containing SYBR green) was mixed with the PCR products, and electrophoresis was performed on a 2% agarose gel to visualize the PCR products. PCR products were mixed in equidensity ratios. Then, the PCR products were purified with a Qiagen Gel Extraction Kit (Qiagen, Germany). The DNA purity and concentration were analyzed by measuring the optical density (OD) at wavelengths of 260 and 280 nm with a NanoPhotometer® spectrophotometer (Implen, Munich, Germany) and then calculating the OD260:OD280 ratio. The DNA concentrations were measured with the Qubit® dsDNA Assay Kit in a Qubit® 2.0 Fluorometer (Life Technologies, Camarillo, CA, United States).
Metagenomic bacterial 16S rRNA gene sequencing assay
Based on many previous studies (32–34), we selected the V3–V4 region to study the microbiome through second-generation sequencing. 16S rRNA gene sequencing was performed by Novogene Bioinformatics Technology Co., Ltd., China. In brief, DNA samples were diluted to a concentration of 1 ng/μl in sterile water and then PCR amplified with the 515F/806R primer set (515F: 5′-GTGCCAGCMGCCGCGGTAA-3′, 806R: 5′-XXXXXXGG ACTACHVGGGTATCTAAT-3′). Sequencing libraries were generated using the TruSeq® DNA PCR-Free Sample Preparation Kit (Illumina, USA) following the manufacturer's recommendations, and index codes were added. Library quality was assessed with a Qubit@ 2.0 Fluorometer (Thermo Scientific) and Agilent Bioanalyzer 2100 system. Finally, the library was sequenced on the Illumina NovaSeq platform, and 250-bp paired-end reads were generated.
Data analysis
Paired-end reads were assigned to samples based on their unique barcodes and were truncated by removing the barcode and primer sequences. Paired-end reads were merged using FLASH (V1.2.7) (35), which is a very fast and accurate analysis tool that was designed to merge paired-end reads when at least some of the reads overlapped with the read generated at the opposite end of the same DNA fragment. The splicing sequences are called raw tags. Quality filtering of raw tags was performed under specific filtering conditions to obtain high-quality clean tags (36) according to the QIIME (V1.9.1) (37) quality control process. The tags were compared with those in the reference database (Silva138 database) using the UCHIME (28) algorithm to detect chimeric sequences, and the chimeric sequences were removed (38). Then, effective tags were finally obtained.
The data were analyzed and compared using SPSS version 27.0 (Chicago, Illinois, USA) and R software (v 2.15.3). Data are summarized as the mean (standard deviation, SD)/median (interquartile range, IQR) for continuous variables depending on normality. Student's t test was used to analyse significant differences in age, sex, height, weight, waist circumference, hip circumference, BMI and lipid levels among the schoolchildren who lived at different altitudes. Analysis of variance (ANOVA) was used to compare the content of short-chain fatty acids (SCFAs) in feces among the four groups. Significant differences in categorical variables are expressed as numbers and percentages. The chi-square test was used to analyse count data. The Kruskal–Wallis test was used to investigate significant differences in operational taxonomic units (OTUs) and the abundance-based coverage estimator (ACE), Chao 1, Simpson's, and Shannon's indexes among the four groups. To correct for type I errors, we applied the Bonferroni method for multiple comparisons between two groups. Differences in microbial community abundances between the obesity group and the control group were analyzed using the Wilcoxon rank sum test, and the significance of these differences was assessed using the false discovery rate (FDR). Principal coordinate analysis (PCoA) was performed using the WGCNA package, stat packages and ggplot2 package in R software (v 2.15.3) to compare the similarity of community structures. Multiresponse permutation procedures (MRPPs) (39) were used to analyse differences in the microbial community structure between groups. Linear discriminant analysis effect size (LEfSe) was performed to identify particular taxa with significantly different abundances between the two groups. An unweighted pair-group method with arithmetic means (UPGMA) clustering was performed as a type of hierarchical clustering method to interpret the distance matrix using average linkage by QIIME (2 v2022.2) software. Correlations between blood lipid levels and intestinal flora diversity were analyzed by Spearman's correlation analysis. p ≤ 0.05 indicated statistical significance.
Results
Basic participant information and serum biochemical indexes
A total of 70 subjects were included in the high-altitude groups (n = 35, normal/obese = 20/15) and low-altitude groups (n = 35, normal/obese = 20/15; Table 1). There was no significant difference in age, sex or height (p > 0.05; Table 1). Body weight, BMI, waist circumference and hip circumference were significantly different between the obesity group and the control group from the same altitude (p < 0.05; Table 1). There was also a significant difference between the obese groups from different altitudes (p < 0.05, Table 1). TG levels were significantly different between the obesity group and the control group from different altitudes (p < 0.05, Table 1), and LDL-C levels were significantly different between the obese groups from different altitudes (p < 0.05, Table 1).
High 16S rRNA gene sequence data preprocessing and quality control
A total of 6,130,974 high-quality and classifiable raw reads were obtained from 70 samples. We retained 5,991,598 clean reads after removing the low-quality sequences. A total of 3,109,972 and 2,881,626 clean reads were obtained for the high-altitude groups and the low-altitude groups, respectively. Clustering analysis assigned the microbial sequences to the same OTU if the samples had at least 97% similarity. We obtained 3,843 OTUs by high-throughput sequencing. These reads were classified into 43 phyla, 101 classes, 226 orders, 331 families, 576 genera, and 340 species. The rarefaction curves showed that our sequencing depth was sufficient (Supplementary Figure S1).
α-Diversity analysis
α-diversity analysis was performed using four indexes, namely, OTUs and the Shannon index, Simpson index, Chao 1 index, and ACE index. There were highly significant differences in all four indexes among the four groups according to the Kruskal–Wallis test (p < 0.05; Table 2). The Shannon and Simpson indexes were significantly different in the different BMI groups from the same altitude according to the Bonferroni multiple comparisons test (p < 0.05; Table 2). The OTUs and Chao 1 and ACE indexes were significantly different between groups with the same BMI who were from different altitudes according to the Bonferroni multiple comparisons test (p < 0.05; Table 2).
β-Diversity analysis
β-diversity analysis was used to evaluate differences in species composition among the samples from different groups. Fecal microbial OTU data were examined by principal coordinate analysis (PCoA) to evaluate similarities among the four groups. Multiresponse permutation procedures (MRPPs) were used to analyse differences in microbial community structures between groups. The results showed that the groups from different altitudes had different microbiota compositions. The LN and LOB groups had similar microbiota compositions, and the sequences could be grouped into the same clusters, whereas the HN and HOB groups did not have similar microbiota compositions, and the sequences could not be grouped into the same clusters. In addition, the results revealed obvious differences in the community compositions and structures in the groups from different altitudes (Table 3 and Figure 1).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/963202/fpubh-10-963202-HTML/image_m/fpubh-10-963202-g001.jpg)
Figure 1. Beta diversity analysis based on UniFrac analysis. Red dot represents the high altitude obese children (HOB) group. Blue dot represents the high altitude normal weight children (HN) group. Green dot represents the low altitude obese children (LOB) group. Orange dot represents the low altitude normal weight children (LN) group. Circles in red, blue, orange, and green represent different periodontal bacterial community clusters, respectively. LN, low-altitude obese children; LOB, low-altitude obese children; HN, high-altitude normal weight children; HOB, high-altitude obese children.
Relative abundance analysis
In this study, we found 484 common OTUs among the four groups, and there were 357, 164, 266, and 313 unique OTUs in the LN, LOB, HN, and HOB groups, respectively. The relative abundance of the microbiota was lowest in the LOB group (Figure 2). Further analysis of the microbiota at the phylum, class, genus and species levels in the four groups showed that five of the top 10 most dominant bacteria at any level were the same among the four groups. The microbiota compositions in all four groups were dominated by the following phyla (order: LN, LOB, HN, HOB): Firmicutes (61.2, 60.1, 60.6, and 53.2%, respectively); Bacteroidota (22.2, 15.4, 28.0, and 30.4%, respectively); Actinobacteriota (11.3, 15.7, 8.5, and 9.9%, respectively); Proteobacteria (3.0, 7.2, 0.9, and 4.2%, respectively); and Verrucomicrobiota (0.3, 0.05, 0.7, and 0.2%, respectively; Figure 3). At the phylum level, the relative abundance of Firmicutes (F) was higher than that of Bacteroides (B) in the obese groups. Bacteroides abundances and F/B ratios were significantly different among the four groups (Supplementary Table S1). At the phylum level, the abundances of Bacteroidota, Proteobacteria, and Verrucomicrobiota were significantly different among the four groups (Table 4). At the class level, the abundances of Clostridia, Negativicutes, Gammaproteobacteria, Bacteroidia, Actinobacteria, and Verrucomicrobiae were significantly different among the four groups (Table 4). At the genus level, the abundances of Megamonas, Prevotella, Akkermansia, and Bifidobacterium were significantly different among the four groups (Table 4). At the species level, the abundances of Romboutsia, Bacteroides caccae, and Escherichia coli were significantly different among the four groups (Table 4).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/963202/fpubh-10-963202-HTML/image_m/fpubh-10-963202-g002.jpg)
Figure 2. (A–E) The comparison between the two groups was done using the Venn chart and multiple groups using petal chart. Different colored circles represent different groups. Overlapping parts of the circles represent common OTUs between two groups, none overlapping parts represent unique OTUs. LN, low-altitude obese children; LOB, low-altitude obese children; HN, high-altitude normal weight children; HOB, high-altitude obese children.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/963202/fpubh-10-963202-HTML/image_m/fpubh-10-963202-g003.jpg)
Figure 3. Community composition of fecal microbiota in four groups. In the bar chart, each bar represents the average relative abundance of each bacterial taxon. The taxa with high relative abundances at the phylum level (top 10, A), class level (top 10, B), genus level (top 10, C) and species level (top 10, D) are shown. LN, low-altitude obese children; LOB, low-altitude obese children; HN, high-altitude normal weight children; HOB, high-altitude obese children.
Dominant gut microbiota composition between children living at different altitudes
We further analyzed differences in the microbiota composition between the two different groups at the genus level using the Metastat analysis method. First, between groups from the same altitude, there were significant differences in the abundances of Akkermansia, Streptococcus, Escherichia-Shigella, Megamonas and Sarcina between the LN and LOB groups (p < 0.05; Figure 4A), and there were significant differences in the abundances of Akkermansia, Ruminococcus torques, Holdemanella, Barnesiella, Fusicatenibacter, Megasphaera, Subdoligranulum, Blautia, Eubacterium hallii, Alistiples, Parabacteroides, Ruminococcus, Escherichia-Shigella and Megamonas between the HN and HOB groups (p < 0.05; Figure 4B). Second, between groups with the same BMI, the abundances of Eubacterium hallii, Clostridium sensu stricto, Akkermansia, Ruminococcus torques, Streptococcus, Escherichia-Shigella, Holdemanella, Klebsiella, Romboutsia, Lachnoclostridium, Parabacteroides, Citrobacter, Megamonas, Sarcina, Catenibacterium, Alloprevotella, Succinlvibrlo and Megasphaera were significantly different between the HN and LN groups (p < 0.05; Figure 4C), and the abundances of Blautia, Eubacterium hallii, Holdemanella, Klebsiella, Romboutsia, Fusicatenibacter, Lachnoclostridium, Ruminococcus, Streptococcus, Escherichia-Shigella, Sarcina, Methanobrevibacter, Catenibacterium, Prevotella, Succinlvibrlo, Megasphaera, Akkermansia, Paraprevotella and Alloprevotella (p < 0.01) were significantly different between the HOB and LOB groups (p < 0.05; Figure 4D). In conclusion, we found that Akkermansia, Holdemanella, Streptococcus, Prevotella, Escherichia-Shigella, and Megamonas were strongly associated with BMI, and these six members of the microbiota plus Sarcina and Catenibacterium were associated with altitude.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/963202/fpubh-10-963202-HTML/image_m/fpubh-10-963202-g004.jpg)
Figure 4. (A–D) MetaStats analysis, in which the abundances of the 35 most significantly different taxa in the two groups are shown in heatmap. Purple dot represents a high significant difference (p < 0.01), red dot represents a significant difference (p < 0.05), and yellow dot represents no significant difference (p > 0.05). LN, low-altitude normal weight children; LOB, low-altitude obese children; HN, high-altitude normal weight children; HOB, high-altitude obese children.
Differences in the gut microbiota composition between the different BMI groups
LEfSe was performed to identify particular taxa with significantly different abundances between groups. The results of the UPGMA analysis showed that there were differences in the abundances of the genera Bifidobacterium, Escherichia-Shigella, Blautia, Romboutsia, Streptococcaceae, Megasphaera, and Prevotella in obese children living at different altitudes. Bifidobacterium, Escherichia-Shigella, Blautia, Romboutsia, and Streptococcaceae were highly abundant in children living at low altitudes, while Megasphaera and Prevotella were abundant in children living at high altitudes. These results indicated that these microbes were closely related to altitude when obesity conditions were the same. We further analyzed intestinal differences in normal-weight children living at different altitudes. Prevotella and Holdemanella were highly abundant at high altitudes. It was interesting to further confirm that Prevotella was positively correlated with high altitude (Figure 5).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/963202/fpubh-10-963202-HTML/image_m/fpubh-10-963202-g005.jpg)
Figure 5. Cladogram based on LEfSe analysis showing the different abundance bacteria between two groups. Comparison between high-altitude obese children (HOB) and low-altitude obese children (LOB) (A). Comparison between high-altitude normal weight children(HN) and low-altitude normal weight children (LN) (C), with a linear discriminant analysis (LDA) threshold of 4 (B,D). The bacteria with significantly different abundance at different taxonomic levels among groups were highlighted by colored circles and shadings.
Correlations between the gut microbiota composition and BMI, TC, TG, HDL-C, and LDL-C levels
A Spearman's correlation matrix was generated to analyse correlations between BMI, TC, TG, HDL-C, and LDL-C levels and the dominant genera in the gut microbiota in all children. Significant relationships were observed between the gut microbiota composition and BMI, TC, TG, and LDL-C levels. Correlation analysis revealed that the abundances of the genera Methanobrevibac, Sarcina, Succinivibrio, Catenibacterium, Alloprevotella, Lactobacillus, Megasphaera, Holdemanella and Prevotella were positively correlated with BMI (p < 0.05; Figure 6), and the abundances of Citrobacter, Klebsiella, Lachnoclostridum, Clostridium-sens, X. Eubacterium.-h, Fusicatenibacter, Blautia, Streptococcus, Escherichia. Shige, Romboutsia, Bifidobacterium and Megamonas were negatively correlated with BMI (p < 0.05; Figure 6). The abundance of Megasphaera was positively correlated with the TC level (p < 0.05; Figure 6), and the abundances of Clostridium and Paraprevotella were negatively correlated with the TC level (p < 0.05; Figure 6). The abundances of Sarcina, Megasphaera and Holdemanella were positively correlated with the TG level (p < 0.05; Figure 6), and the abundances of Fusicatenibacter and Romboutsia were negatively correlated with the TG level (p < 0.05; Figure 6). The abundances of Citrobacter, Klebsiella, Subdoligranulum and Romboutsia were positively correlated with the HDL-C level (p < 0.05; Figure 6), and the abundances of Sarcina, Scuccinivibrio, Catenibacterium, Holdemanella and Prevotella were negatively correlated with the HDL-C level (p < 0.05; Figure 6). The abundances of Succinivibrio and Megasphaera were positively correlated with the LDL-C level (p < 0.05; Figure 6), and the abundances of Clostridium_sens, X. eubacterium._h, Akkermansia, Paraprevotella, Alistipes, Blautia, and Romboutsia were negatively correlated with the LDL-C level (p < 0.05; Figure 6). In addition, the difference in the abundance of Prevotella was positively correlated with the difference in the LDL level between the obese groups from different altitudes (Figure 6).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/963202/fpubh-10-963202-HTML/image_m/fpubh-10-963202-g006.jpg)
Figure 6. Spearman's association analysis of bacterial genera and BMI, TC, TG, HDL, and LDL. r indicates Spearman's correlation coefficient. Cells are colored based on the value of r between significantly altered genera and body mass index (BMI), total cholesterol (TC), tri glyceride (TG), high-density lipoprotein cholesterol (HDL-C), and low-density lipoprotein cholesterol (LDL-C). Red represents a significantly positive correlation, blue represents a significant negative correlation, and yellow represents no significant correlation. *p < 0.05; **p < 0.01.
Functional prediction of gut microbiota in obese Tibetan children living at different altitudes
To further explore the function of gut microbes in children at different altitudes, Tax4Fun was used to analyse the functional prediction of the fecal microbiota in different groups. Figure 7 shows the 35 most abundant pathways at the third level of KEGG pathways. Of these 35 pathways, 10 pathways (amino sugar and nucleotide sugar metabolism; quorum sensing; transcription factors; ABC transporters; transporters; two component system; starch and sucrose metabolism; galactose metabolism; secretion system, bacterial motility proteins.) were enriched in the fecal microbial community of the LOB groups. Eighteen pathways (ribosome biogenesis; prokaryotic defense system; aminoacyl tRNA biosynthesis; transfer RNA biogenesis; DNA replication proteins; purine metabolism; glycine, serine and threonine metabolism; carbon fixation pathways in prokaryotes; pyrimidine metabolism; amino acid-related enzymes; ribosome; chromosome and associated proteins; DNA repair and recombination proteins; mitochondrial biogenesis; exosome; chaperones and folding catalysts; alanine, aspartate and glutamate metabolism; peptidases.) were enriched in the HOB groups, and three pathways (glycolysis/gluconeogenesis, butanoate metabolism and pyruvate metabolism) were enriched in the HN groups. The KEGG pathways for predicting the function of gut microbiota were mainly expressed in four aspects: metabolism, genetic information processing, environmental information processing and cellular processes. The differences in these pathways may be related to the influence of altitude on the abundance of intestinal flora. In addition, the prominence of the butyric acid metabolic pathway in normal weight children at high altitude has attracted our attention.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/963202/fpubh-10-963202-HTML/image_m/fpubh-10-963202-g007.jpg)
Figure 7. Functional profiles of the fecal microbial communities in four groups. The 35 most-abundant pathways at the third level of KEGG pathways are shown in the heatmap. LN, low-altitude normal weight children; LOB, low-altitude obese children; HN, high-altitude normal weight children; HOB, high-altitude obese children.
Correlation analysis between KEGG pathways and blood lipids
Through KEGG pathway annotation and enrichment analysis, the results showed that the butyric metabolism pathway was significantly different in the four groups, with the highest abundance in the HN group, followed by the LN group, and the lowest abundance in the LOB group (Figure 7). Combined with this feature, we measured the short-chain fatty acids in the feces of different groups, and the results showed that acetic, propionic and butyric acid levels were significantly different among the four groups (Table 5). We further carried out a correlation analysis between blood lipids and the content of short-chain fatty acids, and the results showed that acetic acid was correlated with LDL, while butyric acid was correlated with TC, TG, and LDL (Table 6).
Discussion
Obesity is linked to many chronic diseases (2–5). The steadily increasing prevalence of excess body weight among children and adolescents is currently one of the greatest challenges for public health authorities worldwide (40–42). Many studies have shown that obesity in childhood is associated with an increased risk of obesity in adulthood. This situation reminds us to pay attention to the problem of childhood obesity. The potential link between obesity and the composition of the gut microbiota has attracted the attention of many researchers. Many studies have shown that disruption of the gut microbiota composition may be an important cause of obesity (16, 43, 44). However, the gut microbiota is also affected by altitude, and the low prevalence rate of obesity in high-altitude areas is of interest (45, 46). Short or chronic exposure of humans or animals to hypoxic conditions (the most typical characteristic of high altitude) can affect the composition of the intestinal microbiota (47–50). In this study, we focused on the composition, structure and diversity of the gut microbiota in obese Tibetan children aged 7–12 years who lived at different altitudes. By exploring the intestinal bacteria that affect weight change in high-altitude areas, we have provided evidence to guide early obesity prevention measures in the future. In addition, we also analyzed correlations between the gut microbiota composition and lipid metabolism indicators.
Obesity is a health problem among children that contributes to the occurrence of lipid disorders and abnormal blood pressure. It is often accompanied by increases in TC, TG, and LDL-C levels and decreases in HDL-C levels (51). We obtained similar results in our study. Obesity can also cause changes in the gut microbiota composition. In recent years, changes in bacterial strains in the human intestine have been proposed to play a causative role in obesity (52–54). Research by Ley et al. (55) showed that obesity affects the diversity of the gut microbiota, and the relative proportion of Bacteroidetes is decreased in obese people compared with that in lean people (16). These authors also obtained the same results in mice (55), and we drew the same conclusion. We found that at the phylum level, the relative abundance of Firmicutes (F) was higher than that of Bacteroides (B) in the obese groups. The F/B ratios were significantly different between the obese groups. Megamonas is a genus of Firmicutes. In our study, Megamonas was present in only the obese groups, and the relative abundance in the LOB group was higher than that in the HOB group and was positively correlated with BMI and negatively correlated with the LDL-C level. A study conducted in Nanjing (sea level), China, showed that height, weight, age, BMI, TG levels and creatinine levels in children were positively correlated with the relative abundance of Megamonas and that TC, HDL-C and LDL-C levels were negatively correlated with the relative abundance of Megamonas (56). Otoniel Maya-Lucas' research also showed that the relative abundance of Megamonas in obese Mexican people living at an altitude of 2,200 m was correlated with height, weight and BMI (57). Therefore, we concluded that Megamonas is associated with obesity. In addition, it has been reported that Megamonas can ferment glucose to form acetate and propionate, which have been shown to be substrates for lipogenesis and cholesterol formation as energy sources for the host (58). This further explains the high relative abundance of Megamonas in obese subjects. These findings are consistent with our findings. The results of Ma (59) showed that the relative abundance of Megamonas in people in the plains area was higher than that in people in the high-altitude area, consistent with the conclusion of our study. The reason may be related to the impairment of gastrointestinal mucosal barrier function caused by hypoxia, but research on the mechanism is still incomplete.
The gut microbiota is affected by altitude. A high-altitude environment is characterized by low pressure, hypoxia, strong radiation, cold temperatures, etc. Hypoxia can affect the behavior and activities of human beings (60). Hypoxia may also alter the composition of the gut microflora (61), and substantial evidence has implicated both aerobic and facultative anaerobic intestinal bacteria in the dynamic configuration and stability of the anaerobic environment inside the gut (62, 63). Li et al. (25) compared the gut microbiomes of people living at different altitudes and found that the composition of the gut microbiota in individuals living at high altitudes was lower than that in individuals living on the plains. In our study, we found that the abundances of Akkermansia, Holdemanella, Streptococcus, Prevotella, Escherichia-Shigella, Megamonas, Sarcina and Catenibacterium were associated with altitude. Akkermansia, Holdemanella, Prevotella, Sarcina, and Catenibacterium were highly abundant in those living at high altitudes, while Streptococcus, Escherichia-Shigella, and Megamonas were present at significant levels in those living at low altitudes. Of these, Prevotella deserves additional attention. The relative abundance of Prevotella in the high-altitude groups was significantly higher than that in the low-altitude groups. However, when comparing the groups with different BMIs at the same altitudes, the abundance of Prevotella was higher in the HOB group than in the HN group, while the opposite was observed between the low-altitude groups. This study has shown (64) that the Tibetan microbiome is characterized by a relative abundance of Prevotella in individuals living at high altitudes (3,600 m). This difference is likely related to dietary structure. The distribution and prevalence of Prevotella in the human gut is influenced by a variety of factors, including body condition, lifestyle, sex and age (65). Some scholars (66) have stated that the abundance of Prevotella is associated with the consumption of a diet rich in carbohydrates. The Qinghai-Tibet Plateau is rich in highland barley (a high-fiber wheat), which is the main source of carbohydrates among local people. Other studies have shown (67) that the abundance of Prevotella is positively correlated with the consumption of high-fiber foods. However, in our study, to focus on the effect of hypoxia on the gut microbiome, we standardized the diets of the study subjects. High-altitude environments may place greater energy demands on mammals due to thermoregulatory stress than low-altitude environments (68–70). Many anaerobic gut bacteria produce short-chain fatty acids (SCFAs) as end products of polysaccharide fermentation (71). A greater abundance of SCFA-producing obligate anaerobes has been reported in populations who live at high altitudes (25, 64, 72–74) A hypoxic environment leads to the enrichment of facultative anaerobic bacteria (Prevotella) that produce SCFAs and enhances the ability of intestinal bacteria to produce SCFAs to provide more energy, regulate blood pressure, and maintain the stability of the intestinal environment in people who live in high-altitude regions (75). The relative abundance of Prevotella increases with altitude, and this change could be explained by differences in diet, climate, atmospheric oxygen, or any other variable correlated with altitude. Moreno's study (76) provided strong evidence that reduced atmospheric oxygen alone caused an increase in Prevotella abundance. A previous study found that high-altitude exposure altered the gastrointestinal system, causing pathological conditions such as upper gastrointestinal bleeding, ulcers, vomiting, diarrhea, and anorexia. This may be the cause of hypobaric hypoxia-induced gut microbial dysbiosis, which can be ameliorated by prebiotic/probiotic treatment (77). Another study explored the relationship between the hypothalamic-pituitary-thyroid (HPT) axis and gut microbiota composition under hypoxic conditions by simulating the effects of 5,500 m altitude on the HPT axis and gut microbiota in rats. The study found that compared with those in the normoxia group, there were significant differences in the relative abundances of 12 genera in the chronic hypoxia exposure group; thyrotropin-releasing hormone (TRH) and thyroid-stimulating hormone (TSH) concentrations were significantly lower in the hypoxia group than in the control group, and TT4 and TT3 concentrations were significantly higher in the hypoxia group than in the control group. There were significant correlations between the differential bacteria and HPT axis hormones in serum, and Prevotella was significantly negatively correlated with TSH (78). These results provide evidence that the increase in Prevotella abundance at high altitudes may be driven by lower atmospheric oxygen levels. Previous studies have shown that the biodiversity and abundance of the microbiota in obese individuals often differ from those in normal individuals. Our study found that the relative abundance of Prevotella in the obese group was lower than that in the normal group on the plain, but the opposite was true at high altitudes. Duan et al. (79) used 16S rRNA sequencing to compare the intestinal flora compositions of 21 obese individuals from Shandong Province (China) and 21 normal individuals from Beijing to study the characteristics of the intestinal flora in the obese population. The team found that at the species level, there were significant differences in nine species between the control and obese groups. The abundance of Prevotella was significantly increased in the obese population. This was consistent with our results in the groups at high altitudes. Barczyńska R studied the relationship between intestinal bacterial composition and weight in 20 obese children and 20 normal weight children and found that the Prevotella abundance was lower in the fecal microflora of obese children, with an average 30% higher abundance of Prevotella in normal weight children than in obese children. However, the difference was significant in only obese children and not in overweight children (80). Fernández-Navarro et al. (81) analyzed the interrelationships among obesity, diet, oxidative stress, inflammation and the intestinal flora composition in 68 healthy adults and found that the lower abundance of Prevotella in the obese group was associated with elevated proinflammatory and pro-oxidative states. Many other studies have also confirmed the effect of obesity on the intestinal flora composition and the low abundance of Prevotella in obese individuals (82–84). In our study, we found that the abundance of Escherichia-Shigella was also related to altitude. A comparative analysis of the four groups revealed that Escherichia-Shigella exhibited the lowest relative abundance in the HN group and the highest relative abundance in the LOB group, and the difference was statistically significant. Escherichia-Shigella is a member of the gram-negative Enterobacteriaceae family and causes enterobacteriosis, usually resulting in diarrhea and dysentery (85). Escherichia-Shigella is transmitted via the fecal-oral route (86), and poor water supply, lack of basic sanitation and unhygienic behavior have all been associated with Escherichia-Shigella infection. The low economic level and harsh climatic conditions of the Tibetan Plateau result in conditions that are favorable for the spread of Escherichia-Shigella. However, our results showed a low relative abundance of Escherichia-Shigella at high altitudes. Successful invasion of Escherichia-Shigella requires overcoming two gut-specific barriers: the microbiota and the mucus layer. Related studies have suggested that a Prevotella-rich microbiota may have a protective effect against Escherichia-Shigella infection (87). This explains the low Escherichia-Shigella abundance in high-altitude populations. In addition, the higher the inflammation score is, the higher the Escherichia-Shigella abundance (88). In a mouse study, Escherichia-Shigella abundance was positively correlated with blood lipid, glucose and insulin levels (89). This finding indicates that in the HOB group, a greater abundance of Prevotella reduces the inflammatory response caused by obesity. All of these results suggest that the Escherichia-Shigella abundance is higher in obese people than in normal-weight people. Therefore, we believe that the reason for the low prevalence of obesity in high-altitude areas may be related to the abundance of Prevotella.
Finally, we found a correlation between the gut microbe composition and blood lipid levels. Our results showed that the abundances of Sarcina, Megasphaera and Holdemanella were positively correlated with the TG level, while the abundances of Fusicatenibacter and Romboutsia were negatively correlated with the TG level. The abundances of Succinivibrio and Megasphaera were positively correlated with the LDL-C level, while the abundances of Clostridium_sens, X. eubacterium._h, Akkermansia, Paraprevotella, Alistipes, Blautia and Romboutsia were negatively correlated with the LDL-C level. Fu et al. (90) showed that the abundances of Bacteroides, Akkermansia, Desulfovibrio and Parabacteroides were negatively correlated with GLU, TG, TC and HDL-C levels and positively correlated with the LDL-C level. The abundance of Ruminiclostridium was positively correlated with GLU, TG, TC and HDL-C levels but negatively correlated with the LDL-C level. Our results are somewhat similar. In this study, we found that although obese children living at a high altitude also had an abnormal BMI and LDL-C level, their blood lipid levels were lower than those in children living at a low altitude. What is the relationship between blood lipid levels and gut microbe composition in populations who live at high altitudes? In our study, we found that the abundance of Faecalibacterium was higher in the high-altitude groups, especially in the HOB group. Faecalibacterium belongs to Firmicutes, which is fermented and metabolism to produce butyric acid in the intestine. Butyrate provides energy for the body through fatty acid oxidation (91). In this study, we conducted a study on the correlation between butyrate and blood lipids and found that butyrate was correlated with TC, TG, and LDL and negatively correlated with TC, indicating that butyrate has a certain regulatory effect on blood lipids. A study has shown that butyrate is involved in diet-induced obesity and insulin resistance (92) by downregulating the expression and activity of PPAR-γ, promoting a change from lipogenesis to lipid oxidation (93). This is consistent with our results. At the same time, we found that the butyric acid metabolic pathway has a high abundance in normal weight children through the KEGG metabolic pathway and the highest abundance in normal weight children who live in high altitude. Therefore, we inferred that altitude changed the abundance of Faecalibacterium, which further affected the production of butyrate.
Conclusions
In this study, we found that the TG and LDL-C levels in the obesity groups were higher than those in the normal-weight groups, and those in the high-altitude obesity groups were lower than those in the low-altitude obesity groups. Altitude affects the composition and relative abundance of the gut microbiota. Groups living at different altitudes and with different body weights have their own dominant bacterial genera. Megamonas was closely related to obesity, and Prevotella was associated with altitude. Prevotella had an inhibitory effect on the abundance of Escherichia-Shigella. There were correlations between the gut microbiota composition and lipid metabolism indicators, and short-chain fatty acids play an important role in lipid regulation.
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 the Ethics Committee of Qing Hai Provincial People's Hospital, Xining, China. Written informed consent to participate in this study was provided by the participants' legal guardian/next of kin. Written informed consent was obtained from the individual(s), and minor(s)' legal guardian/next of kin, for the publication of any potentially identifiable images or data included in this article.
Author contributions
WD and WZ designed the experiments. YM, RJ, YH, and QZ performed the experiments. WD, ZW, XY, and AA analyzed the data. WD wrote the manuscript. LL provided financial support. All authors contributed to the article and approved the submitted version.
Funding
This research was supported by Project of the Basic Research Program of Qinghai Provincial Science and Technology Department Fund (Gran Number: 2018-ZJ-758), High-end Innovative Talents Thousand Talents Program in Qinghai Province, and Leading Talents and Innovation Team Project Funding.
Conflict of interest
The authors 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.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpubh.2022.963202/full#supplementary-material
References
1. Chooi YC, Ding C, Magkos F. The epidemiology of obesity. Metabolism. (2019) 92:6–10. doi: 10.1016/j.metabol.2018.09.005
2. Rubio-Almanza M, Cámara-Gómez R, Merino-Torres JF. Obesity and type 2 diabetes: also linked in therapeutic options. Endocrinol Diabetes Nutr. (2019) 66:140–9. doi: 10.1016/j.endien.2018.11.006
3. Peters U, Dixon AE, Forno E. Obesity and asthma. J Allergy Clin Immunol. (2018) 141:1169–79. doi: 10.1016/j.jaci.2018.02.004
4. Kachur S, Lavie CJ, de Schutter A, Milani RV, Ventura HO. Obesity and cardiovascular diseases. Minerva Med. (2017) 108:212–28. doi: 10.23736/S0026-4806.17.05022-4
5. Seravalle G, Grassi G. Obesity and hypertension. Pharmacol Res. (2017) 122:1–7. doi: 10.1016/j.phrs.2017.05.013
6. Broughton DE, Moley KH. Obesity and female infertility: potential mediators of obesity's impact. Fertil Steril. (2017) 107:840–7. doi: 10.1016/j.fertnstert.2017.01.017
7. Polyzos SA, Kountouras J, Mantzoros CS. Obesity and nonalcoholic fatty liver disease: from pathophysiology to therapeutics. Metabolism. (2019) 92:82–97. doi: 10.1016/j.metabol.2018.11.014
8. Hopkins BD, Goncalves MD, Cantley LC. Obesity and cancer mechanisms: cancer metabolism. J Clin Oncol. (2016) 34:4277–83. doi: 10.1200/JCO.2016.67.9712
9. Bulbul T, Hoque M. Prevalence of childhood obesity and overweight in Bangladesh: findings from a countrywide epidemiological study. BMC Pediatr. (2014) 14:86. doi: 10.1186/1471-2431-14-86
10. Valdés Pizarro J, Royo-Bordonada MA. Prevalence of childhood obesity in Spain: National Health Survey 2006-2007. Nutr Hosp. (2012) 27:154–60. doi: 10.1590/S0212-16112012000100018
11. Jia P, Xue H, Zhang J, Wang Y. Time trend and demographic and geographic disparities in childhood obesity prevalence in china-evidence from twenty years of longitudinal data. Int J Environ Res Public Health. (2017) 14:E369. doi: 10.3390/ijerph14040369
12. Guo C, Wang H, Feng G, Li J, Su C, Zhang J, et al. Spatiotemporal predictions of obesity prevalence in Chinese children and adolescents: based on analyses of obesogenic environmental variability and Bayesian model. Int J Obes. (2019) 43:1380–90. doi: 10.1038/s41366-018-0301-0
13. Dong S-S, Zhang K, Guo Y, Ding J-M, Rong Y, Feng J-C, et al. Phenome-wide investigation of the causal associations between childhood BMI and adult trait outcomes: a two-sample Mendelian randomization study. Genome Med. (2021) 13:48. doi: 10.1186/s13073-021-00865-3
14. Must A, Jacques PF, Dallal GE, Bajema CJ, Dietz WH. Long-term morbidity and mortality of overweight adolescents. A follow-up of the Harvard Growth Study of 1922 to 1935. N Engl J Med. (1992) 327:1350–1355. doi: 10.1056/NEJM199211053271904
15. Albuquerque D, Nóbrega C, Manco L, Padez C. The contribution of genetics and environment to obesity. Br Med Bull. (2017) 123:159–73. doi: 10.1093/bmb/ldx022
16. Ley RE, Turnbaugh PJ, Klein S, Gordon JI. Microbial ecology: human gut microbes associated with obesity. Nature. (2006) 444:1022–3. doi: 10.1038/4441022a
17. Jumpertz R, Le DS, Turnbaugh PJ, Trinidad C, Bogardus C, Gordon JI, et al. Energy-balance studies reveal associations between gut microbes, caloric load, and nutrient absorption in humans. Am J Clin Nutr. (2011) 94:58–65. doi: 10.3945/ajcn.110.010132
18. Blasco-Baque V, Serino M, Burcelin R. Metabolic therapy at the edge between human hosts and gut microbes. Ann Pharm Fr. (2013) 71:34–41. doi: 10.1016/j.pharma.2012.08.003
19. Chang C-J, Lin C-S, Lu C-C, Martel J, Ko Y-F, Ojcius DM, et al. Ganoderma lucidum reduces obesity in mice by modulating the composition of the gut microbiota. Nat Commun. (2015) 6:7489. doi: 10.1038/ncomms8489
20. 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. doi: 10.1038/nm.4358
21. Fei N, Zhao L. An opportunistic pathogen isolated from the gut of an obese human causes obesity in germfree mice. ISME J. (2013) 7:880–4. doi: 10.1038/ismej.2012.153
22. Zhao Q, Ding Y, Wang J, Gao H, Zhang S, Zhao C, et al. Projecting climate change impacts on hydrological processes on the Tibetan Plateau with model calibration against the glacier inventory data and observed streamflow. J Hydrol. (2019) 573:60–81. doi: 10.1016/j.jhydrol.2019.03.043
23. Zhao P, Xu X, Chen F, Guo X, Zheng X, Liu L, et al. The third atmospheric scientific experiment for understanding the earth–atmosphere coupled system over the Tibetan Plateau and its effects. Bull Am Meteorol Soc. (2018) 99:757—76. doi: 10.1175/BAMS-D-16-0050.1
24. Jia Z, Zhao X, Liu X, Zhao L, Jia Q, Shi J, et al. Impacts of the plateau environment on the gut microbiota and blood clinical indexes in Han and Tibetan Individuals. Am Soc Microbiol. (2020) 5:1–16. doi: 10.1128/mSystems.00660-19
25. Li L, Zhao X. Comparative analyses of fecal microbiota in Tibetan and Chinese Han living at low or high altitude by barcoded 454 pyrosequencing. Sci Rep. (2015) 5:14682. doi: 10.1038/srep14682
26. Santos C, Bustamante A, Hedeker D, Vasconcelos O, Garganta R, Katzmarzyk PT, et al. Correlates of overweight in children and adolescents living at different altitudes: the Peruvian health and optimist growth study. J Obes. (2019) 2019:2631713. doi: 10.1155/2019/2631713
27. CDC. BMI for Children and Teens. Centers for Disease Control and Prevention. (2021). Available online at: https://www.cdc.gov/obesity/basics/childhood-defining.html (accessed May 12, 2022).
28. Attitalla IH. Modified CTAB method for high quality genomic DNA extraction from medicinal plants. Pak J Biol Sci. (2011) 14:998–9. doi: 10.3923/pjbs.2011.998.999
29. Musa DI, Toriola AL, Monyeki MA, Lawal B. Prevalence of childhood and adolescent overweight and obesity in Benue State, Nigeria. Trop Med Int Health. (2012) 17:1369–75. doi: 10.1111/j.1365-3156.2012.03083.x
30. Gu H, Jasbi P, Patterson J, Jin Y. Enhanced detection of short-chain fatty acids using gas chromatography mass spectrometry. Curr Protoc. (2021) 1:e177. doi: 10.1002/cpz1.177
31. Zheng X, Qiu Y, Zhong W, Baxter S, Su M, Li Q, et al. A targeted metabolomic protocol for short-chain fatty acids and branched-chain amino acids. Metabolomics. (2013) 9:818–27. doi: 10.1007/s11306-013-0500-6
32. Claesson MJ, Wang Q, O'Sullivan O, Greene-Diniz R, Cole JR, Ross RP, et al. Comparison of two next-generation sequencing technologies for resolving highly complex microbiota composition using tandem variable 16S rRNA gene regions. Nucleic Acids Res. (2010) 38:e200. doi: 10.1093/nar/gkq873
33. Mizrahi-Man O, Davenport ER, Gilad Y. Taxonomic classification of bacterial 16S rRNA genes using short sequencing reads: evaluation of effective study designs. PLoS ONE. (2013) 8:e53608. doi: 10.1371/journal.pone.0053608
34. Jorth P, Turner KH, Gumus P, Nizam N, Buduneli N, Whiteley M. Metatranscriptomics of the human oral microbiome during health and disease. MBio. (2014) 5:e01012-14. doi: 10.1128/mBio.01012-14
35. Magoč T, Salzberg SL. FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics. (2011) 27:2957–63. doi: 10.1093/bioinformatics/btr507
36. Bokulich NA, Subramanian S, Faith JJ, Gevers D, Gordon JI, Knight R, et al. Quality-filtering vastly improves diversity estimates from Illumina amplicon sequencing. Nat Methods. (2013) 10:57–9. doi: 10.1038/nmeth.2276
37. 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:335–6. doi: 10.1038/nmeth.f.303
38. Haas BJ, Gevers D, Earl AM, Feldgarden M, Ward DV, Giannoukos G, et al. Chimeric 16S rRNA sequence formation and detection in Sanger and 454-pyrosequenced PCR amplicons. Genome Res. (2011) 21:494–504. doi: 10.1101/gr.112730.110
39. O'Reilly FJ, Mielke PW Jr. Asymptotic normality of MRPP statistics from invariance principles of u-statistics. Commun Stat. (1980) 9:629–37. doi: 10.1080/03610928008827907
40. Llewellyn A, Simmonds M, Owen CG, Woolacott N. Childhood obesity as a predictor of morbidity in adulthood: a systematic review and meta-analysis. Obes Rev. (2016) 17:56–67. doi: 10.1111/obr.12316
41. Kumar S, Kelly AS. Review of childhood obesity: from epidemiology, etiology, and comorbidities to clinical assessment and treatment. Mayo Clin Proc. (2017) 92:251–65. doi: 10.1016/j.mayocp.2016.09.017
42. Ward ZJ, Long MW, Resch SC, Giles CM, Cradock AL, Gortmaker SL. Simulation of growth trajectories of childhood obesity into adulthood. N Engl J Med. (2017) 377:2145–53. doi: 10.1056/NEJMoa1703860
43. Reinhardt C, Reigstad CS, Bäckhed F. Intestinal microbiota during infancy and its implications for obesity. J Pediatr Gastroenterol Nutr. (2009) 48:249–56. doi: 10.1097/MPG.0b013e318183187c
44. Angelakis E, Armougom F, Million M, Raoult D. The relationship between gut microbiota and weight gain in humans. Future Microbiol. (2012) 7:91–109. doi: 10.2217/fmb.11.142
45. Voss JD, Masuoka P, Webber BJ, Scher AI, Atkinson RL. Association of elevation, urbanization and ambient temperature with obesity prevalence in the United States. Int J Obes. (2013) 37:1407–12. doi: 10.1038/ijo.2013.5
46. Hirschler V, Martin M, Molinari C, Botta E, Tetzlaff WF, Brites F. Activity of lipoprotein-associated enzymes in indigenous children living at different altitudes. Arch Med Res. (2019) 50:98–104. doi: 10.1016/j.arcmed.2019.07.001
47. Kleessen B, Schroedl W, Stueck M, Richter A, Rieck O, Krueger M. Microbial and immunological responses relative to high-altitude exposure in mountaineers. Med Sci Sports Exerc. (2005) 37:1313–8. doi: 10.1249/01.mss.0000174888.22930.e0
48. Adak A, Maity C, Ghosh K, Pati BR, Mondal KC. Dynamics of predominant microbiota in the human gastrointestinal tract and change in luminal enzymes and immunoglobulin profile during high-altitude adaptation. Folia Microbiol. (2013) 58:523–8. doi: 10.1007/s12223-013-0241-y
49. Lucking EF, O'Connor KM, Strain CR, Fouhy F, Bastiaanssen TFS, Burns DP, et al. Chronic intermittent hypoxia disrupts cardiorespiratory homeostasis and gut microbiota composition in adult male guinea-pigs. EBioMedicine. (2018) 38:191–205. doi: 10.1016/j.ebiom.2018.11.010
50. Zhang W, Jiao L, Liu R, Zhang Y, Ji Q, Zhang H, et al. The effect of exposure to high altitude and low oxygen on intestinal microbial communities in mice. PLoS ONE. (2018) 13:e0203701. doi: 10.1371/journal.pone.0203701
51. Bondyra-Wiśniewska B, Myszkowska-Ryciak J, Harton A. Impact of lifestyle intervention programs for children and adolescents with overweight or obesity on body weight and selected cardiometabolic factors-a systematic review. Int J Environ Res Public Health. (2021) 18:2061. doi: 10.3390/ijerph18042061
52. Festi D, Schiumerini R, Eusebi LH, Marasco G, Taddia M, Colecchia A. Gut microbiota and metabolic syndrome. World J Gastroenterol. (2014) 20:16079–94. doi: 10.3748/wjg.v20.i43.16079
53. Kobyliak N, Virchenko O, Falalyeyeva T. Pathophysiological role of host microbiota in the development of obesity. Nutr J. (2016) 15:43. doi: 10.1186/s12937-016-0166-9
54. Castaner O, Goday A, Park Y-M, Lee S-H, Magkos F, Shiow S-ATE, et al. The gut microbiome profile in obesity: a systematic review. Int J Endocrinol. (2018) 2018:4095789. doi: 10.1155/2018/4095789
55. Ley RE, Bäckhed F, Turnbaugh P, Lozupone CA, Knight RD, Gordon JI. Obesity alters gut microbial ecology. Proc Natl Acad Sci U S A. (2005) 102:11070–5. doi: 10.1073/pnas.0504978102
56. Chen X, Sun H, Jiang F, Shen Y, Li X, Hu X, et al. Alteration of the gut microbiota associated with childhood obesity by 16S rRNA gene sequencing. PeerJ. (2020) 8:e8317. doi: 10.7717/peerj.8317
57. 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. doi: 10.1016/j.anaerobe.2018.10.009
58. Zhang Q, Zou R, Guo M, Duan M, Li Q, Zheng H. Comparison of gut microbiota between adults with autism spectrum disorder and obese adults. PeerJ. (2021) 9:e10946. doi: 10.7717/peerj.10946
59. Ma Y, Ga Q, Ge R-L, Ma S. Correlations between intestinal microbial community and hematological profile in Native Tibetans and Han immigrants. Front Microbiol. (2021) 12:615416. doi: 10.3389/fmicb.2021.615416
60. Basnyat B, Murdoch DR. High-altitude illness. Lancet. (2003) 361:1967–74. doi: 10.1016/S0140-6736(03)13591-X
61. Šket R, Debevec T, Kublik S, Schloter M, Schoeller A, Murovec B, et al. Intestinal metagenomes and metabolomes in healthy young males: inactivity and hypoxia generated negative physiological symptoms precede microbial dysbiosis. Front Physiol. (2018) 9:198. doi: 10.3389/fphys.2018.00198
62. Espey MG. Role of oxygen gradients in shaping redox relationships between the human intestine and its microbiota. Free Radic Biol Med. (2013) 55:130–40. doi: 10.1016/j.freeradbiomed.2012.10.554
63. Albenberg L, Esipova TV, Judge CP, Bittinger K, Chen J, Laughlin A, et al. Correlation between intraluminal oxygen gradient and radial partitioning of intestinal microbiota. Gastroenterology. (2014) 147:1055-63.e8. doi: 10.1053/j.gastro.2014.07.020
64. Li K, Dan Z, Gesang L, Wang H, Zhou Y, Du Y, et al. Comparative analysis of gut microbiota of native Tibetan and Han populations living at different altitudes. PLoS ONE. (2016) 11:e0155863. doi: 10.1371/journal.pone.0155863
65. Tett A, Pasolli E, Masetti G, Ercolini D, Segata N. Prevotella diversity, niches and interactions with the human host. Nat Rev Microbiol. (2021) 19:585–99. doi: 10.1038/s41579-021-00559-y
66. Wu GD, Chen J, Hoffmann C, Bittinger K, Chen Y-Y, Keilbaugh SA, et al. Linking long-term dietary patterns with gut microbial enterotypes. Science. (2011) 334:105–8. doi: 10.1126/science.1208344
67. 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
68. Hochachka PW, Buck LT, Doll CJ, Land SC. Unifying theory of hypoxia tolerance: molecular/metabolic defense and rescue mechanisms for surviving oxygen lack. Proc Natl Acad Sci U S A. (1996) 93:9493–8. doi: 10.1073/pnas.93.18.9493
69. Cheviron ZA, Bachman GC, Connaty AD, McClelland GB, Storz JF. Regulatory changes contribute to the adaptive enhancement of thermogenic capacity in high-altitude deer mice. Proc Natl Acad Sci U S A. (2012) 109:8635–40. doi: 10.1073/pnas.1120523109
70. Schippers M-P, Ramirez O, Arana M, Pinedo-Bernal P, McClelland GB. Increase in carbohydrate utilization in high-altitude Andean mice. Curr Biol. (2012) 22:2350–4. doi: 10.1016/j.cub.2012.10.043
71. Topping DL, Clifton PM. Short-chain fatty acids and human colonic function: roles of resistant starch and nonstarch polysaccharides. Physiol Rev. (2001) 81:1031–64. doi: 10.1152/physrev.2001.81.3.1031
72. Zhang Z, Xu D, Wang L, Hao J, Wang J, Zhou X, et al. Convergent evolution of rumen microbiomes in high-altitude mammals. Curr Biol. (2016) 26:1873–9. doi: 10.1016/j.cub.2016.05.012
73. Lan D, Ji W, Lin B, Chen Y, Huang C, Xiong X, et al. Correlations between gut microbiota community structures of Tibetans and geography. Sci Rep. (2017) 7:16982. doi: 10.1038/s41598-017-17194-4
74. Zhao J, Yao Y, Li D, Xu H, Wu J, Wen A, et al. Characterization of the gut microbiota in six geographical populations of Chinese Rhesus Macaques (Macaca mulatta), implying an adaptation to high-altitude environment. Microb Ecol. (2018) 76:565–77. doi: 10.1007/s00248-018-1146-8
75. Suzuki TA, Martins FM, Nachman MW. Altitudinal variation of the gut microbiota in wild house mice. Mol Ecol. (2019) 28:2378–90. doi: 10.1111/mec.14905
76. Moreno-Indias I, Cardona F, Tinahones FJ, Queipo-Ortuño MI. Impact of the gut microbiota on the development of obesity and type 2 diabetes mellitus. Front Microbiol. (2014) 5:190. doi: 10.3389/fmicb.2014.00190
77. Khanna K, Mishra KP, Chanda S, Ganju L, Singh SB, Kumar B. Effect of synbiotics on amelioration of intestinal inflammation under hypobaric hypoxia. High Alt Med Biol. (2021) 22:32–44. doi: 10.1089/ham.2020.0062
78. Xie Y-L, Mei S, Xiong Y-L, Liu S-Y, Xu C-L. [Effects of simulated 5 500 m hypobaric and hypoxia on HPT axis and intestinal flora in rats]. Zhongguo Ying Yong Sheng Li Xue Za Zhi. (2020) 36:432–7. doi: 10.12047/j.cjap.5898.2020.092
79. Duan M, Wang Y, Zhang Q, Zou R, Guo M, Zheng H. Characteristics of gut microbiota in people with obesity. PLoS ONE. (2021) 16:e0255446. doi: 10.1371/journal.pone.0255446
80. Barczyńska R, Litwin M, Slizewska K, Szalecki M, Berdowska A, Bandurska K, et al. Bacterial microbiota and fatty acids in the faeces of overweight and obese Children. Pol J Microbiol. (2018) 67:339–45. doi: 10.21307/pjm-2018-041
81. Fernández-Navarro T, Salazar N, Gutiérrez-Díaz I, de Los Reyes-Gavilán CG, Gueimonde M, González S. Different intestinal microbial profile in over-weight and obese subjects consuming a diet with low content of fiber and antioxidants. Nutrients. (2017) 9:E551. doi: 10.3390/nu9060551
82. Mbakwa CA, Hermes GDA, Penders J, Savelkoul PHM, Thijs C, Dagnelie PC, et al. Gut microbiota and body weight in school-aged children: the KOALA Birth Cohort Study. Obesity. (2018) 26:1767–76. doi: 10.1002/oby.22320
83. Jonduo ME, Wawae L, Masiria G, Suda W, Hattori M, Takayasu L, et al. Gut microbiota composition in obese and non-obese adult relatives from the highlands of Papua New Guinea. FEMS Microbiol Lett. (2020) 367:fnaa161. doi: 10.1093/femsle/fnaa161
84. Pinart M, Dötsch A, Schlicht K, Laudes M, Bouwman J, Forslund SK, et al. Gut microbiome composition in obese and non-obese persons: a systematic review and meta-analysis. Nutrients. (2021) 14:12. doi: 10.3390/nu14010012
85. Nicolas X, Granier H, Le Guen P. [Shigellosis or bacillary dysentery]. Presse Med. (2007) 36:1606–18. doi: 10.1016/j.lpm.2007.03.001
86. Warren BR, Parish ME, Schneider KR. Shigella as a foodborne pathogen and current methods for detection in food. Crit Rev Food Sci Nutr. (2006) 46:551–67. doi: 10.1080/10408390500295458
87. Baker S, The HC. Recent insights into Shigella. Curr Opin Infect Dis. (2018) 31:449–54. doi: 10.1097/QCO.0000000000000475
88. Aranaz P, Ramos-Lopez O, Cuevas-Sierra A, Martinez JA, Milagro FI, Riezu-Boj JI, et al. predictive regression model of the obesity-related inflammatory status based on gut microbiota composition. Int J Obes. (2021) 45:2261–8. doi: 10.1038/s41366-021-00904-4
89. Song H, Shen X, Chu Q, Zheng X. Vaccinium bracteatum Thunb. fruit extract reduces high-fat diet-induced obesity with modulation of the gut microbiota in obese mice. J Food Biochem. (2021) 45:e13808. doi: 10.1111/jfbc.13808
90. Fu J, Wang Y, Tan S, Wang J. Effects of banana resistant starch on the biochemical indexes and intestinal flora of obese rats induced by a high-fat diet and their correlation analysis. Front Bioeng Biotechnol. (2021) 9:575724. doi: 10.3389/fbioe.2021.575724
91. Miquel S, Martín R, Rossi O, Bermúdez-Humarán LG, Chatel JM, Sokol H, et al. Faecalibacterium prausnitzii and human intestinal health. Curr Opin Microbiol. (2013) 16:255–61. doi: 10.1016/j.mib.2013.06.003
92. Ríos-Covián D, Ruas-Madiedo P, Margolles A, Gueimonde M, de Los Reyes-Gavilán CG, Salazar N. Intestinal short chain fatty acids and their link with diet and human health. Front Microbiol. (2016) 7:185. doi: 10.3389/fmicb.2016.00185
Keywords: obesity, high altitude, schoolchildren, gut microbiome, metabolism
Citation: Du W, Liu L, Ma Y, Zhu Q, Jia R, Han Y, Wu Z, Yan X, Ailizire A and Zhang W (2022) Analysis of the gut microbiome in obese native Tibetan children living at different altitudes: A case–control study. Front. Public Health 10:963202. doi: 10.3389/fpubh.2022.963202
Received: 07 June 2022; Accepted: 10 November 2022;
Published: 24 November 2022.
Edited by:
Nur A. Hasan, EzBiome Inc., United StatesReviewed by:
Farzam Vaziri, University of California, Davis, United StatesLifeng Zhu, Nanjing Normal University, China
Copyright © 2022 Du, Liu, Ma, Zhu, Jia, Han, Wu, Yan, Ailizire and Zhang. 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: Wei Zhang, zw0228@qhu.edu.cn