Skip to main content

Alteration of gut microbiota associated with hypertension in children

Abstract

Background

The association of disturbance in gut microbiota with hypertension (HTN) defined on three separate occasions among children and adolescents remains unclear. In this study, we aimed to compare the differences in gut microbiota composition and diversity between children with HTN and those with normal blood pressure (BP).

Methods

Data and stool samples were collected from the second follow-up of a childhood cardiovascular health cohort study in 2021. 16 S ribosomal RNA gene sequencing was conducted to determine the relative abundance of microbial taxa in 51 children aged 10–14 years with HTN and 51 children with normal BP.

Results

Compared with children with normal BP, those with HTN had decreased gut microbiome diversity. At the genus level, after adjusting for the false discovery rate (FDR), the proportions of several gut microbiota such as Blautia (PFDR=0.042), Coprococcus (PFDR=0.042), Eubacterium_ventriosum_group (PFDR=0.027), Christensenellaceae_R-7_group (PFDR=0.027), and norank_f__Lachnospiraceae (PFDR=0.015) significantly decreased in children with HTN compared to those with normal BP. Receiver operating characteristic analysis, net reclassification improvement (NRI) and integrated discrimination improvement (IDI) were performed and showed that the genera norank_f__Lachnospiraceae and Dorea significantly enhanced the ability of body mass index to differentiate between children with HTN and those with normal BP (area under the receiver operating characteristic curve: 0.95, 95% confidence interval 0.91–0.99; NRI > 0; IDI = 0.12, P < 0.05). Phylogenetic Investigation of Communities by Reconstruction of Unobserved States showed that the mean proportions of cofactors and vitamins metabolism pathway and the glycan anabolism pathway were higher in children with HTN.

Conclusions

Disturbances in the abundance and diversity of gut microbiota may contribute to the development of HTN in children. Gut microbiota biomarkers may be of significant importance in the early identification and diagnosis of childhood HTN.

Clinical trial number

Not applicable.

Peer Review reports

Introduction

Cardiovascular disease (CVD), the leading cause of mortality worldwide, accounted for nearly one-third (i.e., 17.9 million) of total deaths in 2019 [1]. Hypertension (HTN), a major risk factor for CVD, affects 1.28 billion adults aged 30 to 79 years worldwide, but only 42% are diagnosed and treated [2]. Moreover, an estimated 16 million people die from HTN-related complications each year [3]. Therefore, early prevention of HTN is likely to lead to a significant reduction in the burden of CVD and mortality later in life.

HTN can develop in childhood and track into adulthood [4,5,6]. Moreover, elevated blood pressure (BP) and HTN in childhood are associated with subclinical markers of CVD [7, 8] and fatal and non-fatal CVD events [9]. A diagnosis of HTN in children relies on measurements from at least three separate occasions followed by the application of complex tables that stratify individuals’ systolic and diastolic BP levels by age, sex, and height percentile [10,11,12]. This approach poses numerous challenges to implementation [13] in community- and clinic-based settings, which leads to the underdiagnosis of pediatric HTN [14, 15]. Therefore, exploring alternatives to detect HTN could facilitate diagnosis and lead to the development of new therapeutic strategies to reduce the short- and long-term CVD burden associated with HTN.

With only 5.7% of BP variability accounted for by genetic variants [16], environmental factors, such as diet [17], physical activity [18], and sleep [19], are important modifiable factors in the development of HTN. Recently, it has been reported that disturbance of gut microbial community diversity, taxonomic composition, and related pathways was associated with high BP in animal models [20, 21] and human adults [22,23,24]. Despite these associations observed in adult populations and animals, the association between gut microbiota disturbances and HTN, as determined through three separate occasions among children and adolescents, remains poorly understood.

In this study, we aimed to compare the diversity, composition, and abundance of gut microbiota by sequencing the 16 S ribosomal RNA (16 S rRNA) of stool samples from 51 children with HTN, as determined through three separate occasions, and 51 children with normal BP from the second follow-up of the “Huantai Childhood Cardiovascular Health Cohort Study”. Furthermore, we selected gut microbiota biomarkers to identify HTN in children.

Methods

Participants and sample collection

Children were recruited from the second follow-up of the “Huantai Childhood Cardiovascular Health Cohort Study”, which was conducted in Zibo, Shandong, China, from November to December 2021. A total of 51 children aged 10–14 years with HTN (i.e., HTN group) and 51 randomly selected children with normal BP (i.e., NBP group) were included in this study. All participants had no history of antibiotic use, vomiting, diarrhea, and gastrointestinal diseases within the past three months. Informed consent was obtained from all participants and their guardians. This study was approved by the Ethics Committee of Shandong University (Number: 20160308), China.

Measurement and definition

After a rest of at least 10 mins, seated BP values were measured by trained staff using a verified electronic sphygmomanometer (Omron HBP-1300) [25] in a quiet room using a suitable cuff according to the individual’s right upper arm circumference. Three consecutive BP measurements were recorded and subsequent measurements were performed after at least 1 min. If the difference between any two of the three BP readings exceeded 10 mmHg, an additional BP measurement was performed. The mean value of the last two readings was used for analysis.

Elevated BP was defined using sex-, age-, and height-specific 95th percentile cutoffs of systolic BP (SBP) and/or diastolic BP (DBP) for Chinese children [10]. Participants with elevated BP at the first visit underwent a second visit at least 2 weeks later. A third visit was conducted at least another 2 weeks later if elevated BP persisted at the second visit, according to the same procedures. Children who had elevated BP on all three separate occasions were identified as HTN [11, 12].

Covariates

A self-designed structured questionnaire, as used previously [26], was completed by participants to collect information on demographic characteristics (e.g., sex, age), diet (frequency of breakfast consumption, fruit, vegetables, and carbonated beverages during the past 30 days, intake of western food during the past 7 days, intake of dairy products, yogurt, aquatic products, nuts, spicy foods, and beans/bean products during the past 12 months, numbers of eggs consumed per day, and the salty taste of participants), and lifestyle variables (e.g., the average frequency of exercise per week [more than 10 mins each time], screen time per day, sleep duration per day, and secondhand smoke exposure at home during the past 7 days).

Weight and height were measured twice using an electronic ultrasonic height and weight scale (HGM-300) with participants in light clothing after removing hats and shoes. Weight was measured to the nearest 0.1 kg and height to the nearest 0.1 cm. The mean of the two measurements was used for data analysis. Body mass index (BMI, kg/m2) was calculated as weight divided by the square of height. Waist circumference (WC) was measured twice using a standard non-elastic tape positioned horizontally 1 cm above the umbilicus at the end of expiration. WC was measured to the nearest 0.1 cm, with the mean of replicate measures used for analysis.

Blood samples were collected after fasting for at least 10 h. Blood biochemistry (e.g., fasting plasma glucose [FPG], high-density lipoprotein cholesterol [HDL-C], low-density lipoprotein cholesterol [LDL-C], total cholesterol [TC], and triglyceride [TG]) was determined using an automatic electrochemiluminescence immunoanalyzer (Cobas e 601) according to the kit instructions.

Fecal sample collection and 16 S rRNA gene sequencing

A sampling spoon was used to pry off the surface of the stool, and fresh internal fecal samples were collected from each participant using a 2 ml enzyme-free cryopreservation tube without enzymatic substances that can catalyze specific chemical reactions, thus effectively preventing contamination of samples with enzymes and protecting the integrity and activity of the samples. During the entire operation process, avoid direct contact between the sampling spoon and hands or other objects. Fecal samples were then flash-frozen with liquid nitrogen and stored at -80 ºC until microbiome sequencing and analysis.

Microbial DNA was extracted from each homogenized fecal sample using DNA Stool Kit and then detected using 1% agarose gel electrophoresis. The primers 338 F (5’-ACTCCTACGGGAGGCAGCAG-3’) and 806R (5’-GGACTACHVGGGTWTCTAAT-3’) were used to amplify the V3-V4 region of the bacterial 16 S rRNA gene [27, 28]. The designated hypervariable sequencing region was amplified according to the protocol of TransStart FastPfu DNA Polymerase (TransGen AP221-02). Each sample was carried out with three replicates. After recovering with AxyPrep DNA Gel Recovery Kit (Axygen Corp., USA) and eluting with Tris-HCl, the PCR products of the same sample were mixed and then detected by 2% agarose gel electrophoresis [29].

According to the preliminary electrophoretic quantitative results, the PCR products were detected and quantified utilizing the QuantiFluor™ - St Blue Fluorescence Quantitative System (Promega Corp., USA). The resulting PCR amplicons were prepared using TruSeq™ DNA Sample Prep Kit (Illumina, USA), and the library was then sequenced on an Illumina MiSeq platform.

Bioinformatics analysis

After quality control, adapter trimming, quality filtering, and data filtering of the sequencing data using the Fastp tool (version 0.19.6, https://github.com/OpenGene/fastp) [30], Fast Length Adjustment of Short Reads (FLASH, version 1.2.11, https://ccb.jhu.edu/software/FLASH/index.shtml) was used to paired-end reads merging [31]. The subsequent microbial community and diversity analysis was conducted using the Quantitative Insights Into Microbial Ecology (QIIME, version 1.9.1, http://qiime.org/install/index.html).

First, during preprocessing, we conducted Operational Taxonomic Units (OTU) clustering to assign sequences to different microbial taxa. Second, at a similarity level of 97%, non-repeating sequences (excluding single sequences) were clustered into OTU using Uparse software (version 7.0.1090, http://drive5.com/uparse/). After removing the chimera during clustering, we obtained representative sequences of OTUs. Third, we calculated the relative abundance of each OTU within each sample, which was the number of sequences assigned to that OTU divided by the total number of sequences in the sample. This approach eliminates the influence of different sequencing depths on comparisons across samples, resulting in the normalization of the microbial abundances. Based on the RDP classifier Bayesian algorithm Classifier (version 2.11, http://sourceforge.net/projects/rdp-classifier/) at a confidence threshold of 0.7, species classification annotation of representative OTU sequences was conducted according to the SILVA database (Release138, http://www.arb-silva.de).

Statistical analyses

Continuous variables including age, height, weight, BMI, WC, SBP, DBP, screen time, sleep duration, FPG, HDL-C, LDL-C, TC, and TG were presented as mean ± standard deviation or median (P25, P75), and the differences between HTN and NBP were examined using t-test (for data normally distributed or nearly normally distributed) or rank sum test (for data non-normally distributed). Categorical variables, such as sex, were presented as n (%) and the difference between the two groups was examined using the chi-square test. SAS v9.4 (SAS Institute, Cary, NC, USA) was used for data analyses and a two-sided P < 0.05 was considered statistically significant.

16 S rRNA microbial analysis was conducted by Majorbio Bio-Pharm Technology Co. Ltd. (Shanghai, China). Alpha diversity indexes were calculated to assess the community richness (e.g., Chao and Ace) and diversity (e.g., Shannon) at the OTU level based on Mothur v1.30.2 (https://www.mothur.org/wiki/Download_mothur) [32]. Rarefaction curves were constructed by randomly subsampling a certain number of sequences from each sample and plotting the alpha diversity index (not limited to the number of OTUs) against the number of sampled sequences, providing an overview of the diversity trends across different groups. To minimize the effects of sequencing depth on alpha and beta diversity measures, the numbers of 16 S rRNA gene sequences from each sample were rarefied to 26,837.

Species Venn diagram analysis (R package “VennDiagram”) [33] was used to count the number of shared and unique OTUs between HTN and NBP groups. The composition and relative abundance of gut microbiota at the phylum level and genus level were visualized based on the community pie plots. The Principal Co-ordinate Analysis (PCoA, R package “Vegan”) [34] was performed to visualize the beta diversity of bacterial communities based on weighted-unifrac distance matrices (FastUniFrac, http://UniFrac.colorado.edu/). Adonis analysis (QIIME, http://qiime.org/install/index.html) was used to compare the differences between groups.

After adjusting for false discovery rate (FDR), we used the non-parametric Wilcoxon rank-sum test to compare the relative abundance of gut microbiota at the genus, phylum, and OTU levels, respectively, between NBP and HTN groups. The redundancy analysis (RDA, R package “Vegan”) [35] was used to reflect the association of gut microbiota with multiple environmental factors, including WC, BMI, intake frequency of breakfast, fruit, vegetables, carbonated beverages, western food, dairy products, yogurt, aquatic product, nuts, spicy foods, and beans/bean products, numbers of eggs consumed per day, the salty taste of participants, the frequency of exercise, secondhand smoke exposure, screen time, sleep duration, and levels of FPG, HDL-C, LDL-C, TC, and TG. The Random Forest algorithm was used, via the R package “randomForest” [36], to screen out predominant biomarkers that could distinguish HTN from NBP based on the area under the curve (AUC) validation method. We used different combinations of top-ranked species to construct random forest models and evaluated their ability to discriminate between the two groups. The AUC values indicate the predictive power of these models. By varying the number of species included, we found that certain combinations resulted in higher AUC values, suggesting better discrimination between groups. Among the top n important biomarkers in the barplot of variable importance, logistic regression models were used to evaluate the odds ratio (OR) and its 95% confidence intervals (CI) between the biomarkers and childhood HTN, after adjusting for key environmental factors. Those significant biomarkers were selected to perform receiver operating characteristic (ROC) curve analyses to evaluate their performance in discriminating HTN from NBP children. In addition, net reclassification improvement (NRI) and integrated discrimination improvement (IDI) [37] were calculated to quantitatively evaluate the improvement in discriminating HTN from NBP children using significant intestinal biomarkers, based solely on BMI, which is an important predictor of childhood HTN. The aforementioned analyses were conducted using SAS v9.4 software. To illustrate the metabolic functions of the gut microbiota, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis was performed using Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt) 2 v2.2.0 (https://github.com/picrust/picrust2/) [38].

Results

Basic characteristics

Compared with children with NBP, those with HTN had higher levels of BMI, WC, and TG (P < 0.05), shorter sleep duration, and lower HDL-C (P < 0.05). No significant differences were found in age, sex, FPG, LDL-C, TC, screen time, diet variables, physical activity, and exposure to secondhand smoke at home (Table 1).

Table 1 Characteristics of children with hypertension and children with normal blood pressure

Rarefaction curve analysis

A total of 5,160,048 optimization sequences and 2,125,606,197 bases with an average length of 412 (min length: 404; max length: 426) were generated. The rarefaction curve (Ace, Chao, and Shannon indices) tended to be flat, indicating that the amount of sequencing data in this study was sufficient (Fig. 1A-C).

Gut microbial diversity, richness, and composition

We evaluated alpha diversity based on measurements of the community richness (Chao and Ace) and diversity (Shannon) of the gut microbiota in two groups. The values of Chao (P = 0.006) and Ace (P = 0.022) showed that children with NBP had a higher community richness than those with HTN (Fig. 1D-E). Shannon index also indicated a higher diversity in children with NBP compared with those with HTN (P = 0.005, Fig. 1F). Higher indices of alpha diversity were observed in boys and older children (> 12 years) who had NBP, similar to the main results. However, no significant differences in alpha diversity were observed between HTN and NBP groups in girls and younger children (≤ 12 years) (Table S1). Microbial community analysis identified a total of 874 OTUs from 102 fecal samples. Of these, 620 (70.94%) were shared by all samples, while 153 (17.51%) and 101 (11.56%) OTUs, respectively, were uniquely identified in children with NBP and HTN (Fig. 1G).

Community analysis showed that at the phylum level, the proportion of Bacteroidota (17.18% vs. 6.69%, P = 0.006) increased, but the proportion of Firmicutes (62.88% vs. 71.82%, P = 0.030) decreased in children with HTN compared with those with NBP, and the differences disappeared after adjustment of FDR (Fig. 2A, B, C). The proportions of Bacteroidota and Firmicutes, along with the P-values after adjustment of FDR, are shown in Table S2. The proportion of gut microbiota at the genus level was also different in the HTN group compared to the NBP group (Fig. 2D, E). For example, the abundance of Blautia (PFDR=0.042), Coprococcus (PFDR=0.042), Eubacterium_ventriosum_group (PFDR=0.027), Christensenellaceae_R-7_group (PFDR=0.027), norank_f__Lachnospiraceae (PFDR=0.015) decreased in the HTN group (Fig. 2F). Figure 2G presents the differences in the relative abundance of OTUs between groups. The PCoA analysis also revealed that the gut microbiota based on the OTU profile in the HTN group could be significantly distinguished from the NBP group (R2 = 0.052, P = 0.002; Fig. 2H).

Fig. 1
figure 1

Differences in essential features and α diversity of gut microbiota between children with HTN and those with NBP. The rarefaction curves based on the (A) Ace index, (B) Chao index, and (C) Shannon index were used to demonstrate the amount of sequencing data in this study; Alpha diversity was evaluated using (D) Ace; (E) Chao; and (F) Shannon indexes. (G) Venn diagram illustrating the distribution of unique or common OTUs. * P < 0.05; ** P < 0.01

Fig. 2
figure 2

Differences in the composition of gut microbiota and β diversity between children with HTN and those with NBP. Composition of gut microbiota in children with HTN (A, D) and those with NBP (B, E) at the phylum and genus levels; Wilcoxon rank-sum test for examining the differences in gut microbiota at phylum (C), genus (F), and OTU (G) levels, respectively; PCoA analysis for evaluating β diversity based on the Bray-Curtis distance algorithm on the OTU level (H). Original P-values are shown in (C) since the differences disappeared after the adjustment of FDR.* P < 0.05, ** P < 0.01

Identification of children with hypertension based on gut microbiota

The RDA analysis revealed that the 7 environmental factors, including frequency of breakfast consumption (r2 = 0.106, P = 0.007), sleep duration (r2 = 0.087, P = 0.011), secondhand smoke exposure (r2 = 0.095, P = 0.014), BMI (r2 = 0.078, P = 0.022), LDL-C (r2 = 0.070, P = 0.031), WC (r2 = 0.074, P = 0.035), and fruit intake (r2 = 0.070, P = 0.038), were significantly associated with the community distribution of gut microbiota (Fig. 3A). Forty-six important features at the genus level were selected using random forest analyses (Fig. 3B) and the models for distinguishing hypertensive children achieved a maximum AUC of 0.92 when the top 8 important features were included sequentially (Fig. 3C). After adjusting for these 7 environmental factors, the logistic regression analysis showed that genus norank_f__Lachnospiraceae (OR = 0.984, 95% CI: 0.972–0.997, P  = 0.015) and genus Dorea (OR = 0.997, 95% CI: 0.994–0.999, P  = 0.019) were inversely associated with childhood HTN after the further adjustment of lifestyle factors and BMI (Table S3). Our ROC analysis showed that BMI had high utility in distinguishing children with HTN from those with NBP (AUC = 0.91, 95% CI: 0.85–0.97, Fig. 3D). When genus norank_f__Lachnospiraceae and genus Dorea were combined, the AUC increased to 0.95 (95% CI = 0.91–0.99, Fig. 3E). We further calculated NRI and IDI and found that genera norank_f__Lachnospiraceae and Dorea could significantly enhance the ability of BMI to distinguish between children with HTN and children with NBP (category NRI = 0.49, 95% CI: 0.32–0.66; category-free NRI = 0.75, 95% CI: 0.39–1.10). In addition, the combination of BMI with genera norank_f__Lachnospiraceae and Dorea significantly enhanced the comprehensive discriminative ability by 12% (IDI = 0.12, P < 0.05) (Table 2).

Table 2 NRI and IDI analysis to assess the predictive performance of BMI alone and BMI combined gut microbiota in identifying children with hypertension
Fig. 3
figure 3

Screening for biomarkers of gut microbiota after adjustment for environmental factors, and evaluation of the utility of gut microbiota in distinguishing children with HTN from those with NBP. (A) The RDA analysis examines the association between environmental factors and the distribution of gut microbiota. (B) The random forest algorithm for screening intestinal microbial markers at the genus level. (C) The AUC of models developed using top-ranked important features sequentially. (D) The utility of BMI in distinguishing children with HTN from those with NBP (E) The utility of the combination of BMI with genera norank_f__Lachnospiraceae and Dorea in distinguishing children with HTN from those with NBP

Metabolic pathway analyses by KEGG

PICRUSt2 showed that cofactors and vitamins metabolism (such as riboflavin metabolism, vitamin B6 metabolism, etc.), as well as glycan biosynthesis and metabolism (such as lipopolysaccharide biosynthesis, glycosphingolipid biosynthesis, etc.), were enriched in the HTN group, while terpenoids and polyketides metabolism (such as limonene and pinene degradation) were enriched in the NBP group (Fig. 4A, B).

Fig. 4
figure 4

Difference in metabolic pathway analyses between children with HTN and those with NBP by KEGG. KEGG pathway analysis to illustrate the metabolic function differences of the gut microbiota between HTN and NBP in the (A) primary pathway and the (B) secondary pathway

Discussion

In this study, we found that alterations in gut microbiota composition and diversity were associated with HTN in children, varied by age and sex. In addition, the combination of BMI, genus norank_f__Lachnospiraceae, and genus Dorea had a high ability to distinguish children with HTN from those with NBP, with an AUC of 0.95. Our study suggests potential associations between the cofactors and vitamins metabolism pathway, glycan anabolism pathway, and terpenoids and polyketides metabolism pathway and HTN in children, which require further validation to establish causal associations.

It has been demonstrated that gut microbial diversity has decreased in hypertensive animal models and adult humans [21,22,23,24]. For example, a microbiome study of the Coronary Artery Risk Development in Young Adults Study found that microbial diversity and richness were inversely associated with SBP and HTN [23]. Another cross-sectional observational study including 315 “metabolically healthy” and “metabolically unhealthy” children aged 8–14 years showed a significantly lower alpha diversity in children with elevated SBP than children with normal SBP [39]. However, elevated SBP mentioned above was defined as ≥ 1.5SD above the mean SBP z-score based on BP measurements taken on only one occasion, which might overestimate the true prevalence of HTN [40, 41]. In this study, we first found that, among children with HTN defined by BP measurement on three separate occasions, there was a significantly lower diversity and richness in gut microbiota compared with children with NBP. Our findings indicate an association between lower diversity and richness in children with HTN, which warrants further investigation as a potential target for intervention.

When we conducted an analysis stratified by sex, only boys with HTN exhibited a lower alpha diversity in gut microbiota, while girls did not. Reduction in diversity of gut microbiota was found among older children with HTN at puberty (> 12 years), but not among prepubertal children. This finding may be attributed to the role of gut microbiota in regulating pubertal timing and sex hormone levels [42]. Several studies have reported that sex may influence gut microbiota differences, however, findings have been inconsistent [28, 43, 44]. A large cohort study including 551 Chinese adults reported that females with obesity exhibited increased alpha diversity, whereas this was not observed among males [28]. However, another study involving 75 coronary heart disease patients from the CORDIOPREV study (39 men and 36 post-menopausal women) revealed that when BMI was ≤ 33 kg/m2, women exhibited lower alpha diversity, whereas males did not. The opposite was observed when BMI was > 33 kg/m2 [43]. Given the differences in sample size, disease outcomes, and race/ethnicity, the previous findings related to sex differences should be interpreted with caution. Further studies are warranted to elucidate the role of sex in the association of gut microbiota with HTN among children.

Similar to previous studies conducted in animals [21, 45] and humans [22, 46, 47], we found that Coprococcus, a fiber-fermenting bacterium that produces butyrate and propionate, was decreased in children with HTN compared with those with NBP, as well as another butyric acid-producing bacterium Blautia. Although a study of 92 Caucasian participants in Moscow aged 25–76 years reported that glucose metabolism disturbance and arterial HTN were associated with a higher abundance of Blautia [48], other studies have reported decreased abundance among Chinese patients with preeclampsia [49], increased abundance among Chinese young adults on a lower-fat diet [50], and decreased abundance among adults with higher TC, LDL-C, and non-HDL-C levels [50]. These findings might help explain the association between decreased abundance of Blautia and HTN in children. In addition, two previous Mendelian randomization studies found an association between the decreased abundance of Eubacterium_ventriosum_group and elevated SBP level [51] and eclampsia [52]. A cross-sectional study of 29 hypertensive adults and 32 adults with normal BP from Spain showed that the abundance of Christensenellaceae R-7 group was negatively correlated with levels of SBP and DBP [53]. We similarly found that Chinese children with NBP had a higher abundance of Eubacterium_ventriosum_group and Christensenellaceae_R-7_group, suggesting that these gut microbes might have protective effects on cardiovascular health.

To the best of our knowledge, no prior study has examined the association of norank_f_Lachnospiraceae and Dorea with HTN in both adults and children. Experimental studies conducted on mice have demonstrated that the abundance of norank_f_Lachnospiraceae decreased in immunosuppressive or high-fat-fed mice [54, 55] as well as in mice with elevated levels of inflammatory cytokines (e.g., TNF-α, IL-6, and IL-1β) [56]. Dorea belongs to the class Clostridia within the phylum Firmicutes whose abundance positively associates with energy harvesting and obesity [57]. In this study, we found that after adjusting for environmental factors and BMI, genus norank_f_Lachnospiraceae and Dorea exhibited an independent negative association with childhood HTN. Furthermore, we found that the combination of BMI, genera norank_f_Lachnospiraceae and Dorea, could accurately differentiate children with HTN from those with NBP. Our findings indicate that the combination of BMI with norank_f_Lachnospiraceae and Dorea offers high diagnostic utility in identifying children with HTN. To compensate for the limited sensitivity of AUC to model changes, we further conducted NRI and IDI analysis [58, 59] to quantitatively evaluate the ability of genera norank_f_Lachnospiraceae and Dorea to enhance the ability of BMI in distinguishing between children with HTN and those with NBP. Consistent with the results of ROC analysis, the combination of BMI with genera norank_f_Lachnospiraceae and Dorea increased the discriminatory ability between children with HTN and those with NBP by 12%, suggesting that these genera significantly enhance the screening ability of BMI for children with HTN. This indirectly suggests an association of norank_f_Lachnospiraceae and Dorea with the development of childhood HTN.

We identified potential associations between HTN in children and the metabolic pathways of cofactors and vitamins metabolism pathway, glycan anabolism pathway, and terpenoids and polyketides metabolism pathway. Studies have demonstrated that glycans play multiple roles in the initiation and progression of stroke, coronary artery disease, as well as peripheral vascular disease [60]. Several glycans exhibited differences between hypertensive individuals and healthy controls [61]. Specifically, a positive correlation was observed between triantennary glycans and BP levels in adults from China and Croatia, while a negative correlation was observed in glycans containing core-fucose [62]. Cofactors and vitamins play pivotal roles in modulating enzyme activity. It is increasingly recognized that deficiencies and disruptions in cofactors (e.g., magnesium [63], coenzyme Q10 [64]) and vitamins (e.g., vitamin D [65, 66]) are associated with an elevated risk of HTN. However, most of the abovementioned findings were derived from studies conducted on adults with HTN. To the best of our knowledge, our study is the first to demonstrate that the cofactors and vitamins metabolism pathway is enriched in children with HTN compared to children with NBP. However, these findings require further validation. Our findings simply a potential link between HTN and gut microbiota via the abovementioned pathways.

To the best of our knowledge, we are the first to evaluate the gut microbiota associated with HTN in children defined by BP measurements taken on three separate occasions. In addition, we found sex- and age-related differences in gut microbiota among children with HTN. However, several limitations should be considered. First, this study was conducted in one city in China. Given the variation in the prevalence of HTN in different race/ethnicity [67] and age [68] groups, our results should be cautiously extended to populations from other regions or races/ethnicities. Second, as fecal sample collection occurred simultaneously with anthropometric and laboratory measurements, we could not investigate the temporal association between gut microbiota and HTN. Prospective studies are necessary to confirm if the identified microbiota alterations occur before the onset and progression of HTN in childhood. Third, the metabolic pathway analyses were based on PICRUSt, which provides speculative functional inferences at the species level, and thus these findings should be interpreted with caution. Finally, this study, employing 16 S rRNA gene sequencing technology, was limited to the sequencing resolution at the genus level and prevented us from performing in-depth analyses of the associated genetic, functional, and metabolic pathways.

Conclusion

In conclusion, disturbance in gut microbiota was associated with HTN in children, varied by age and sex. In addition, the combination of BMI and specific gut microbes, including genera norank_f__Lachnospiraceae and Dorea provided high diagnostic utility in distinguishing HTN from NBP children. Further research could explore the underlying mechanisms of this association and its long-term implications for HTN management in children.

Data availability

The datasets generated and/or analyzed during the current study are available in the [National Library of Medicine; National Center for Biotechnology Information; BioProject] repository, [Persistent web link: https://www.ncbi.nlm.nih.gov/bioproject/PRJNA939951 OR Accession number: PRJNA939951].

Abbreviations

16 S rRNA:

16 S ribosomal RNA

AUC:

Area under the curve

BMI:

Body mass index

BP:

Blood pressure

CI:

Confidence intervals

CVD:

Cardiovascular disease

DBP:

Diastolic blood pressure

FDR:

False discovery rate

FPG:

Fasting plasma glucose

HDL-C:

High-density lipoprotein cholesterol

HTN:

Hypertension

IDI:

Integrated discrimination improvement

KEGG:

Kyoto encyclopedia of genes and genomes

LDL-C:

Low-density lipoprotein cholesterol

NRI:

Net reclassification improvement

OR:

Odds ratio

OTU:

Operational taxonomic units

PCoA:

Principal co-ordinate analysis

PICRUSt:

Phylogenetic investigation of communities by reconstruction of unobserved states

RDA:

Redundancy analysis

ROC:

Receiver operating characteristic

SBP:

Systolic blood pressure

TC:

Total cholesterol

TG:

Triglyceride

WC:

Waist circumference

References

  1. World Health Organization. Cardiovascular diseases (CVDs). https://www.who.int/news-room/fact-sheets/detail/cardiovascular-diseases-(cvds). Accessed 7 January 2023.

  2. World Health Organization. Hypertension. https://www.who.int/news-room/fact-sheets/detail/hypertension. Accessed 7 January 2023.

  3. GBD 2017 Causes of Death Collaborators. Global, regional, and National age-sex-specific mortality for 282 causes of death in 195 countries and territories, 1980–2017: a systematic analysis for the global burden of disease study 2017. Lancet. 2018;392(10160):2170–2170.

    Google Scholar 

  4. Gartlehner G, Vander Schaaf EB, Orr C, Kennedy SM, Clark R, Viswanathan M. Screening for hypertension in children and adolescents: updated evidence report and systematic review for the US preventive services task force. JAMA. 2020;324(18):1884–95.

    Article  PubMed  Google Scholar 

  5. Chen X, Wang Y. Tracking of blood pressure from childhood to adulthood: a systematic review and meta-regression analysis. Circulation. 2008;117(25):3171–80.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Juhola J, Oikonen M, Magnussen CG, Mikkila V, Siitonen N, Jokinen E, et al. Childhood physical, environmental, and genetic predictors of adult hypertension: the cardiovascular risk in young Finns study. Circulation. 2012;126(4):402–9.

    Article  PubMed  Google Scholar 

  7. Daniels SR, Loggie JM, Khoury P, Kimball TR. Left ventricular geometry and severe left ventricular hypertrophy in children and adolescents with essential hypertension. Circulation. 1998;97(19):1907–11.

    Article  CAS  PubMed  Google Scholar 

  8. Lande MB, Carson NL, Roy J, Meagher CC. Effects of childhood primary hypertension on carotid intima media thickness: a matched controlled study. Hypertension. 2006;48(1):40–4.

    Article  CAS  PubMed  Google Scholar 

  9. Jacobs DR Jr., Woo JG, Sinaiko AR, Daniels SR, Ikonen J, Juonala M, et al. Childhood cardiovascular risk factors and adult cardiovascular events. N Engl J Med. 2022;386(20):1877–88.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Fan H, Yan Y, Mi J, On behalf of the Chiese Child Blood Pressure References Collaborative Group. Updating blood pressure references for Chinese children aged 3–17 years. Chin J Hypertens. 2017;25(5):428–35.

    Google Scholar 

  11. Flynn JT, Kaelber DC, Baker-Smith CM, Blowey D, Carroll AE, Daniels SR, et al. Clinical practice guideline for screening and management of high blood pressure in children and adolescents. Pediatrics. 2017;140(3):e20171904.

    Article  PubMed  Google Scholar 

  12. Xi B, Zong XN, Kelishadi R, Hong YM, Khadilkar A, Steffen LM, et al. Establishing international blood pressure references among nonoverweight children and adolescents aged 6 to 17 years. Circulation. 2016;133(4):398–408.

    Article  PubMed  Google Scholar 

  13. Bello JK, Mohanty N, Bauer V, Rittner SS, Rao G. Pediatric hypertension: provider perspectives. Glob Pediatr Health. 2017;4:2333794X17712637.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Hansen ML, Gunn PW, Kaelber DC. Underdiagnosis of hypertension in children and adolescents. JAMA. 2007;298(8):874–9.

    Article  CAS  PubMed  Google Scholar 

  15. Moin A, Mohanty N, Tedla YG, Carroll AJ, Padilla R, Langman CB, et al. Under-recognition of pediatric hypertension diagnosis: examination of 1 year of visits to community health centers. J Clin Hypertens (Greenwich). 2021;23(2):257–64.

    Article  PubMed  Google Scholar 

  16. Evangelou E, Warren HR, Mosen-Ansorena D, Mifsud B, Pazoki R, Gao H, et al. Genetic analysis of over 1 million people identifies 535 new loci associated with blood pressure traits. Nat Genet. 2018;50(10):1412–25.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Sacks FM, Svetkey LP, Vollmer WM, Appel LJ, Bray GA, Harsha D, et al. Effects on blood pressure of reduced dietary sodium and the dietary approaches to stop hypertension (DASH) diet. DASH-Sodium collaborative research group. N Engl J Med. 2001;344(1):3–10.

    Article  CAS  PubMed  Google Scholar 

  18. Whelton SP, Chin A, Xin X, He J. Effect of aerobic exercise on blood pressure: a meta-analysis of randomized, controlled trials. Ann Intern Med. 2002;136(7):493–503.

    Article  PubMed  Google Scholar 

  19. St-Onge MP, Grandner MA, Brown D, Conroy MB, Jean-Louis G, Coons M, et al. Sleep duration and quality: impact on lifestyle behaviors and cardiometabolic health: A scientific statement from the American heart association. Circulation. 2016;134(18):e367–86.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Adnan S, Nelson JW, Ajami NJ, Venna VR, Petrosino JF, Bryan RM Jr., et al. Alterations in the gut microbiota can elicit hypertension in rats. Physiol Genomics. 2017;49(2):96–104.

    Article  CAS  PubMed  Google Scholar 

  21. Yang T, Santisteban MM, Rodriguez V, Li E, Ahmari N, Carvajal JM, et al. Gut dysbiosis is linked to hypertension. Hypertension. 2015;65(6):1331–40.

    Article  CAS  PubMed  Google Scholar 

  22. Li J, Zhao F, Wang Y, Chen J, Tao J, Tian G, et al. Gut microbiota dysbiosis contributes to the development of hypertension. Microbiome. 2017;5(1):14.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Sun S, Lulla A, Sioda M, Winglee K, Wu MC, Jacobs DR Jr., et al. Gut microbiota composition and blood pressure. Hypertension. 2019;73(5):998–1006.

    Article  CAS  PubMed  Google Scholar 

  24. Yan Q, Gu Y, Li X, Yang W, Jia L, Chen C, et al. Alterations of the gut Microbiome in hypertension. Front Cell Infect Microbiol. 2017;7:381.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Meng L, Zhao D, Pan Y, Ding W, Wei Q, Li H, et al. Validation of Omron HBP-1300 professional blood pressure monitor based on auscultation in children and adults. BMC Cardiovasc Disord. 2016;16:9.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Wang H, Sun J, Zhang Z, Yang L, Zhao M, Bovet P, et al. Waist circumference change and risk of high carotid intima-media thickness in a cohort of Chinese children. J Hypertens. 2021;39(9):1901–7.

    Article  CAS  PubMed  Google Scholar 

  27. Zeng S, Cao J, Wei C, Chen Y, Liu Q, Li C, et al. Polysaccharides from Artocarpus heterophyllus lam. (jackfruit) pulp alleviate obesity by modulating gut microbiota in high fat diet-induced rats. Food Hydrocolloids. 2023;139:108521.

    Article  CAS  Google Scholar 

  28. Gao X, Zhang M, Xue J, Huang J, Zhuang R, Zhou X, et al. Body mass index differences in the gut microbiota are gender specific. Front Microbiol. 2018;9:1250.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Li JW, Dong T, Keerthisinghe TP, Chen H, Li M, Chu WQ, et al. Long-term Oxytetracycline exposure potentially alters brain thyroid hormone and serotonin homeostasis in zebrafish. J Hazard Mater. 2020;399:123061.

    Article  CAS  PubMed  Google Scholar 

  30. Chen SF, Zhou YQ, Chen YR, Gu J. Fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics. 2018;34(17):884–90.

    Article  Google Scholar 

  31. Magoc T, Salzberg SL, Flash. Fast length adjustment of short reads to improve genome assemblies. Bioinformatics. 2011;27(21):2957–63.

  32. Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, et al. Introducing Mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol. 2009;75(23):7537–41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Chen H, Boutros PC. VennDiagram: a package for the generation of highly-customizable Venn and Euler diagrams in R. BMC Bioinformatics. 2011;12:35.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Jing Z, Zheng W, Jianwen S, Hong S, Xiaojian Y, Qiang W, et al. Gut microbes on the risk of advanced adenomas. BMC Microbiol. 2024;24(1):264.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Kambura AK, Mwirichia RK, Kasili RW, Karanja EN, Makonde HM, Boga HI. Bacteria and Archaea diversity within the hot springs of lake Magadi and little Magadi in Kenya. BMC Microbiol. 2016;16(1):136.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Ren Z, Fan Y, Li A, Shen Q, Wu J, Ren L, et al. Alterations of the human gut Microbiome in chronic kidney disease. Adv Sci (Weinh). 2020;7(20):2001936.

    Article  CAS  PubMed  Google Scholar 

  37. Li J, Gao M, D’Agostino R. Evaluating classification accuracy for modern learning approaches. Stat Med. 2019;38(13):2477–503.

    Article  PubMed  Google Scholar 

  38. Zhao R, Liu J, Feng J, Li X, Li B. Microbial community composition and metabolic functions in landfill leachate from different landfills of China. Sci Total Environ. 2021;767:144861.

    Article  CAS  PubMed  Google Scholar 

  39. Alcazar M, Escribano J, Ferre N, Closa-Monasterolo R, Selma-Royo M, Feliu A, et al. Gut microbiota is associated with metabolic health in children with obesity. Clin Nutr. 2022;41(8):1680–8.

    Article  CAS  PubMed  Google Scholar 

  40. Liu K, Li C, Gong HB, Guo Y, Hou BJ, Chen LY, et al. Prevalence and risk factors for hypertension in adolescents aged 12 to 17 years: A School-Based study in China. Hypertension. 2021;78(5):1577–85.

    Article  CAS  PubMed  Google Scholar 

  41. Dong J, Dong H, Yan Y, Cheng H, Zhao X, Mi J, et al. Prevalence of hypertension and hypertension phenotypes after three visits in Chinese urban children. J Hypertens. 2022;40(7):1270–7.

    CAS  PubMed  Google Scholar 

  42. Korpela K, Kallio S, Salonen A, Hero M, Kukkonen AK, Miettinen PJ, et al. Gut microbiota develop towards an adult profile in a sex-specific manner during puberty. Sci Rep. 2021;11(1):23297.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Haro C, Rangel-Zuniga OA, Alcala-Diaz JF, Gomez-Delgado F, Perez-Martinez P, Delgado-Lista J, et al. Intestinal microbiota is influenced by gender and body mass index. PLoS ONE. 2016;11(5):e0154090.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Markle JG, Frank DN, Mortin-Toth S, Robertson CE, Feazel LM, Rolle-Kampczyk U, et al. Sex differences in the gut Microbiome drive hormone-dependent regulation of autoimmunity. Science. 2013;339(6123):1084–8.

    Article  CAS  PubMed  Google Scholar 

  45. Li HB, Yang T, Richards EM, Pepine CJ, Raizada MK. Maternal treatment with Captopril persistently alters Gut-Brain communication and attenuates hypertension of male offspring. Hypertension. 2020;75(5):1315–24.

    Article  CAS  PubMed  Google Scholar 

  46. Kim S, Rigatto K, Gazzana MB, Knorst MM, Richards EM, Pepine CJ, et al. Altered gut Microbiome profile in patients with pulmonary arterial hypertension. Hypertension. 2020;75(4):1063–71.

    Article  CAS  PubMed  Google Scholar 

  47. Silveira-Nunes G, Durso DF Jr, Cunha L, Maioli EHM, Vieira TU. Hypertension is associated with intestinal microbiota dysbiosis and inflammation in a Brazilian population. Front Pharmacol. 2020;11:258.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Kashtanova DA, Tkacheva ON, Doudinskaya EN, Strazhesko ID, Kotovskaya YV, Popenko AS, et al. Gut microbiota in patients with different metabolic statuses: Moscow study. Microorganisms. 2018;6(4):98.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Chang Y, Chen Y, Zhou Q, Wang C, Chen L, Di W, et al. Short-chain fatty acids accompanying changes in the gut Microbiome contribute to the development of hypertension in patients with preeclampsia. Clin Sci (Lond). 2020;134(2):289–302.

    Article  CAS  PubMed  Google Scholar 

  50. Wan Y, Wang FL, Yuan JH, Li J, Jiang DD, Zhang JJ, et al. Effects of dietary fat on gut microbiota and faecal metabolites, and their relationship with cardiometabolic risk factors: a 6-month randomised controlled-feeding trial. Gut. 2019;68(8):1417–29.

    Article  CAS  PubMed  Google Scholar 

  51. Zhang J, He J, Liao Y, Xia X, Yang F. Genetic association between gut Microbiome and blood pressure and blood cell count as mediator: A two-step Mendelian randomization analysis. Gene. 2024;925:148573.

    Article  CAS  PubMed  Google Scholar 

  52. Wu X, Li Q, Lin D, Cai J, Huang H, Tan H. Gut microbiota and hypertensive disorders in pregnancy: evidence from the Mendelian randomization study. Aging. 2023;15(17):9105–27.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Calderon-Perez L, Gosalbes MJ, Yuste S, Valls RM, Pedret A, Llaurado E, et al. Gut metagenomic and short chain fatty acids signature in hypertension: a cross-sectional study. Sci Rep. 2020;10(1):6436.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Sun Y, Wang F, Liu Y, Liu S, An Y, Xue H, et al. Microbiome-metabolome responses of Fuzhuan brick tea crude polysaccharides with immune-protective benefit in cyclophosphamide-induced immunosuppressive mice. Food Res Int. 2022;157:111370.

    Article  CAS  PubMed  Google Scholar 

  55. Jing N, Liu X, Jin M, Yang X, Hu X, Li C, et al. Fubrick tea attenuates high-fat diet induced fat deposition and metabolic disorder by regulating gut microbiota and caffeine metabolism. Food Funct. 2020;11(8):6971–86.

    Article  CAS  PubMed  Google Scholar 

  56. Wang K, Li W, Wang K, Du B, Hu Z, Zhao L. Litchi thaumatin-like protein induced the liver inflammation and altered the gut microbiota community structure in mice. Food Res Int. 2022;161:111868.

    Article  CAS  PubMed  Google Scholar 

  57. Ley RE, Backhed F, Turnbaugh P, Lozupone CA, Knight RD, Gordon JI. Obesity alters gut microbial ecology. P Natl Acad Sci USA. 2005;102(31):11070–5.

    Article  CAS  Google Scholar 

  58. Pencina MJ, D’Agostino RB, Sr., D’Agostino RB Jr., Vasan RS. Evaluating the added predictive ability of a new marker: from area under the ROC curve to reclassification and beyond. Stat Med. 2008;27(2):157–72. discussion 207– 112.

  59. Kennedy KF, Pencina MJ. A SAS® macro to compute added predictive ability of mew markers predicting a dichotomous outcome. In: SESUG 2010: Savannah, Georgia; Proceeding of SouthEast SAS User Groups: 2010.

  60. Varki AP, Baum LG, Bellis SL, Cummings RD, Esko JD, Hart GW, et al. Working group report: the roles of glycans in hemostasis, inflammation and vascular biology. Glycobiology. 2008;18(10):747–9.

  61. Wang YX, Klaric L, Yu XW, Thaqi K, Dong J, Novokmet M, et al. The association between glycosylation of Immunoglobulin G and hypertension A multiple ethnic Cross-Sectional study. Medicine. 2016;95(17):e3379.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Lu JP, Knezevic A, Wang YX, Rudan I, Campbell H, Zou ZK, et al. Screening novel biomarkers for metabolic syndrome by profiling human plasma N-glycans in Chinese Han and Croatian populations. J Proteome Res. 2011;10(11):4959–69.

    Article  CAS  PubMed  Google Scholar 

  63. Joosten MM, Gansevoort RT, Mukamal KJ, Kootstra-Ros JE, Feskens EJ, Geleijnse JM, et al. Urinary magnesium excretion and risk of hypertension: the prevention of renal and vascular end-stage disease study. Hypertension. 2013;61(6):1161–7.

    Article  CAS  PubMed  Google Scholar 

  64. Gonenc A, Hacisevki A, Tavil Y, Cengel A, Torun M. Oxidative stress in patients with essential hypertension: a comparison of dippers and non-dippers. Eur J Intern Med. 2013;24(2):139–44.

    Article  CAS  PubMed  Google Scholar 

  65. Wang TJ, Pencina MJ, Booth SL, Jacques PF, Ingelsson E, Lanier K, et al. Vitamin D deficiency and risk of cardiovascular disease. Circulation. 2008;117(4):503–11.

    Article  CAS  PubMed  Google Scholar 

  66. Forman JP, Giovannucci E, Holmes MD, Bischoff-Ferrari HA, Tworoger SS, Willett WC, et al. Plasma 25-hydroxyvitamin D levels and risk of incident hypertension. Hypertension. 2007;49(5):1063–9.

    Article  CAS  PubMed  Google Scholar 

  67. Sorof JM, Lai D, Turner J, Poffenbarger T, Portman RJ. Overweight, ethnicity, and the prevalence of hypertension in school-aged children. Pediatrics. 2004;113(3 Pt 1):475–82.

    Article  PubMed  Google Scholar 

  68. Song P, Zhang Y, Yu J, Zha M, Zhu Y, Rahimi K, et al. Global prevalence of hypertension in children: A systematic review and Meta-analysis. JAMA Pediatr. 2019;173(12):1154–63.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Thank you to the Center for Disease Control and Prevention of Zibo City and the Center for Disease Control and Prevention of Huantai County for their support of this study. Thank you to all the participating children and their parents/guardians from Huantai, China.

Funding

This work was supported by the National Natural Science Foundation of China (82204066 and 82173538); and the Natural Science Foundation of Shandong Province (ZR2021QH272).

Author information

Authors and Affiliations

Authors

Contributions

LY drafted the manuscript. JS and BX contributed to the study design. JS analyzed, accessed and verified the data and had full access to all the data in the study. All authors contributed to the interpretation of the data, revised the manuscript, and approved the final version of the manuscript.

Corresponding author

Correspondence to Bo Xi.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Ethics Committee of Shandong University (Number: 20160308), China. Informed consent was obtained from all participants and their guardians. This study was conducted in compliance with the principles of the Declaration of Helsinki.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

: Supplementary material 1.doc: Table S1 shows the comparison of alpha diversity between children with hypertension and children with normal blood pressure by sex and age. Table S2 shows the comparison of community proportion in gut microbiota between children with hypertension and children with normal blood pressure at the phylum level. Table S3 shows the association between the genera identified by the Random Forest analysis and childhood hypertension.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, 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 you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. 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-nc-nd/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Sun, J., Yang, L., Ma, C. et al. Alteration of gut microbiota associated with hypertension in children. BMC Microbiol 25, 282 (2025). https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12866-025-03999-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12866-025-03999-1

Keywords