Skip to main content

Characteristics and variation of fecal bacterial communities and functions in isolated systolic and diastolic hypertensive patients

Abstract

Background

Hypertension (HTN) is one of the major cardiovascular risk factors, which contributes to increasing target organ damages and cardiovascular morbidity and mortality worldwide. Isolated systolic HTN (ISH) and isolated diastolic HTN (IDH) are two important subtypes of HTN. Previous researches have demonstrated the alteration of fecal bacteria in HTN, but not down to these two sub-types. In order to identify whether the composition of bacterial taxa and functional modules shift in ISH and IDH, we performed a metagenomic sequencing analysis of fecal samples from 15 controls, 14 ISH, and 11 IDH.

Results

Compared with control and ISH, IDH patients showed decreased gene number, bacterial richness, and evenness, although the bacterial alterations did not reach statistical significance in the Shannon index. Also, at the genus level, the β-diversity for intestinal flora in IDH was distinguishable from those with ISH. Furthermore, the taxonomic composition of ISH or IDH was different from that of healthy control at genus and species levels. Patients with IDH or ISH were confirmed to be enriched with Rothia mucilaginosa, along with reduced Clostridium sp. ASBs410. Lastly, the altered KEGG modules were significantly decreased in IDH compared with the control group, such as sodium transport system; while for ISH, functions relevant to biotin biosynthesis were decreased.

Conclusions

Overall, our results showed the disordered fecal bacteria profiles in subjects with ISH and especially IDH, emphasizing the significance of early intervention for IDH.

Peer Review reports

Introduction

Hypertension (HTN) is one of the major cardiovascular risk factors, which substantially contributed to increasing target organ damage, cardiovascular morbidity and mortality worldwide [1]. Accompanied by increasing urbanization and related changes in lifestyle, rising income, and population aging, the burden of HTN has been on the rise, especially in China [2,3,4]. A large representative survey recently showed that the prevalence of HTN among adults aged ≥18 years reached up to 23.2% in China, affecting approximately 244.5 million individuals [5]. Accumulating evidence has been indicating that both genetic susceptibility and environmental factors are implicated in the development of HTN [6, 7]. In particular, considering the acknowledgment of the role of fecal bacteria in non-communicable diseases, such as metabolic diseases deepen [8,9,10], the relationship between fecal bacteria and HTN has also been advanced in recent years.

Changes in the bacterial composition do not necessary indicate a dysbiosis condition, and bacterial community might find new stable configuration that not result in dysbiotic conditions according to the concepts of microbial resilience proposed by investigators in recent years [11]. Resilience, together with resistance and functional redundancy, are the ability for the intestinal microbiota to restore into the original function, and is known as key properties of a robust microbiota [12]. And dysbiosis as a new state, refers to an imbalance in the taxonomic composition of gut microbiota, which occurs when resilience of the community fails, and would lead to various disorders such as chronic inflammatory diseases and metabolic syndrome [11, 12]. Emerging investigations have identified the aberrant fecal bacteria in various animal models [13, 14], and hypertensive patients [15], such as a decreased bacterial diversity, disordered microbial structure and functions. For instance, significant decrease in microbial richness and evenness was detected in hypertensive animals and individuals [13,14,15]. In addition, compared with the healthy controls, Prevotella-dominated gut enterotype, distinct metagenomic composition with reduced bacteria related to healthy status and overgrowth of bacteria such as Prevotella and Klebsiella in hypertensive populations was reported [15]. Furthermore, elevated blood pressure (BP) was demonstrated to be elicited by transferring intestinal microbiota from hypertensive human donors to germ-free mice. The causal role of fecal bacteria dysbiosis in contributing to HTN was confirmed [15].

Populations suffering HTN display a wide phenotypic variability, and different types of HTN might be attributable to distinct underlying mechanisms. To the best of our knowledge, there are isolated systolic HTN (ISH) and isolated diastolic HTN (IDH) where patients exhibit abnormal elevation in solely systolic BP (SBP) or diastolic BP (DBP). Individuals with these subtypes may be characterized in hemodynamic or structural abnormalities, respectively, during the development of HTN [16]. In detail, ISH has been frequently associated with enhanced stiffness in large artery [17], while IDH is more relevant to increased peripheral vascular resistance [18]. However, the changes in composition and function of gut microbes in human patients suffering ISH or IDH remain unexplored. Therefore, we proposed a hypothesis that the composition of bacterial taxa and functional modules might be altered in ISH or IDH patients. Thus, we herein conducted a shotgun metagenomic analysis of fecal bacteria in ISH and IDH along with healthy controls.

Materials and methods

Study cohort

The subjects were enrolled from a previous study [15]. All the individuals in the present study were from a cohort study among employees of the Kailuan Group Corporation. None of the them was under antihypertensive treatment. Patients suffering from cancer, heart failure, renal failure, stroke, and peripheral artery disease were excluded. Furthermore, individuals were also excluded if they had received antibiotics or probiotics within the last 2 months. In the current study, the participants were divided into three groups: 15 healthy controls (SBP ≤125 mmHg, or DBP ≤80 mmHg), 11 patients with IDH (DBP ≥90 mmHg and SBP < 140 mmHg) [19] and 14 patients with ISH (SBP ≥140 mmHg and DBP < 90 mmHg) [19]. Ethics approval was obtained from Kailuan General Hospital, Beijing Chaoyang Hospital. All the participants signed informed consent forms prior to data collection.

DNA extraction and library preparation

Stool samples were freshly collected in a sterile container from each participant and immediately froze at − 20 °C. Samples were shipped to the laboratory with ice pack and stored at − 80 °C until analysis. On the basis of manufacturer’s recommendations, DNA was extracted using TIANGEN kit from Novogene Bioinformatics Technology Co., Ltd. Qualified DNA was fragmented by ultrasonic processor, and library of approximately 300 bp clone insert sizes was constructed per sample. Paired-end metageomic sequencing was conducted by Illumina platform with read length 150 bp, insert size 300. The reads aligned to the human genome (alignment with Short Oligonucleotide Analysis Package 2 [SOAP2], Version 2.21, parameters: -s 135, −l 30, −v 7, −m 200, −× 400) were removed after quality control, and remaining high-quality reads were used for analysis.

Bacterial diversity analysis

Rarefaction curves were conducted to evaluate whether the sample size is enough by randomly selecting a certain number of individuals from the samples and counting the number of species represented by these individuals. In addition, to identify whether the sequencing quantity is sufficient, rarefaction curves by gradually expanding the sequencing depth of random sampling was also performed. The α-diversity parameters of Shannon index, Chao richness, and Pielou evenness were computed using the vegan R package (version 3.3.3) at the genus and species levels, differences of α-diversity indices among groups were tested by Kruskal-Wallis test. Also, β-diversity metrics, including nonmetric dimensional scaling (NMDS), principal-component analysis (PCA), principal coordinate analysis (PCoA) and Permutational Multivariate Analysis of Variance (PERMANOVA) were carried out based on abundances of the microbes at genus and species levels. NMDS was calculated by the vegan package, PCA by the FactoMineR package, PCoA and PERMANOVA by the vegan and ape packages in R software (version 3.3.3). Furthermore, as previously described [20, 21], we evaluated differences in the distribution of β-diversity on axes by wilcoxon rank sum test.

Taxonomic annotation and abundance profiling

DIAMOND (version 0.7.9.58, default parameters except that −k 50 − sensitive −e 0.00001) was used to align the clean reads to the integrated NR database [22]. Only genes with e-values ≤10 × e-value of the top hit mapped reads were retained to distinguish taxonomic groups [15]. The taxonomic level of each gene was generated following the procedure described by Huson et al. [23]. The taxonomic relative abundance profile was calculated by summing the abundance of genes annotated to the same feature.

Functional analysis

Functional annotation was carried out by DIAMOND (Version 0.7.9.58, default parameters except that −k 50 − sensitive −e 0.00001) against the Kyoto Encyclopedia of Genes and Genomes (KEGG) (Release 73.1, with animal and plant genes removed) database [24]. Each protein was assigned to KEGG modules by the highest scoring annotated hit (s) containing at least one high-scoring segment pair scoring > 60 bits [25].

Statistical analysis

Quantitative variables were expressed as median and interquartile, and Wilcoxon rank-sum test was used for between-group comparisons. Categorical variables were presented as numbers and compared using Chi-square test. Differential abundance of genera, species, and KEGG modules were determined using the Kruskal-Wallis test, and P values were corrected for multiple comparisons using the Benjamini-Hochberg method, shown as q values. Correlations between shared differential signatures in ISH and IDH with clinical factors were estimated by Spearmen’s correlation analysis. All statistical analyses were performed by R software (version 3.3.3), and P or q < 0.05 was regarded as statistical significance.

Results

Baseline characteristics of the study cohort

In the current study, 40 participants, including 15 healthy controls, 14 ISH, and 11 IDH, were recruited from the previous study [15]. The clinical characteristics of all participants are shown in Table 1. Briefly, compared to controls, ISH and IDH patients presented with extremely higher SBP and DBP, respectively (P values < 0.01). Also, both SBP and DBP levels between ISH and IDH were dramatically different as expected. In addition, there was no significant difference in other clinical characteristics, such as body mass index, uric acid, creatinine, fasting blood glucose, total cholesterol, triglyceride, low-density lipoprotein, and high-density lipoprotein between groups.

Table 1 General characteristics of study participants

Gut microbial diversity of ISH and IDH

A total of 223.03 Gb 125 bp paired-end reads were generated from the raw data (244.67 Gb), with an average of 5.5759 ± 0.6260 (s.d.) million reads per sample (Table S1). The assembled long contigs or scaffolds by high-quality sequencing readings were used for gene prediction, taxonomic classification, and functional annotation. Rarefaction curves were conducted to evaluate whether the sample size is enough by randomly selecting a certain number of individuals from the samples and counting the number of species represented by these individuals. We found the curves tend to be flat, indicating that the sample size is sufficient and reasonable, and only a small number of new species characteristics would be produced by more sampling (Fig. S1A). In addition, to identify whether the sequencing quantity is sufficient, rarefaction curves by gradually expanding the sequencing depth of random sampling was performed. Fig. S1B demonstrated that the curves approached saturation as the sample sequencing depth increases, and thus the amount of sequencing data is sufficient and stable. For the α-diversity, although we did not find significant difference among three groups, participants with IDH had deficient genes and a trend toward lower Shannon index, Chao richness and Pielou evenness compared with healthy controls and patients with ISH at both genus (Fig. 1b-d) and species (Fig. 1e-g) levels.

Fig. 1
figure 1

Shifts of intestinal bacteria in α- and β-diversity in patients with IDH and ISH. a-d Box plots show gene number, Shannon index, Chao richness, and Pielou evenness at each group’s genus level. e-g Box plots show Shannon index, Chao richness and Pieloue evenness at each group’s species level. Green, con, control, n = 15; blue, ISH, n = 14; orange, IDH, n = 11. h-j NMDS, PCA, and PCoA plots based on the genera level in groups. Significant differences across groups were established at NMDS2, PCA1, and PCoA2. Circles in green indicate samples from control, squares in blue indicate samples from ISH, rhombus in orange represent individuals from IDH. k-m Scatter diagram showing NMDS, PCA, and PCoA plots at the species level. Boxes represent the interquartile ranges; the inside lines represent the median; circles represent outliers. For gene number and α-diversity indices, P values were obtained by the Kruskal-Wallis test. For β-diversity metrics, P values were obtained by PERMANOVA tests; while the differences in the distribution of β-diversity on axes were obtained by wilcoxon rank sum test

Furthermore, β-diversity observed in NMDS, PCA, and PCoA plots at genus and species levels suggested that participant grouping had a certain effect on bacteria structure (Fig. 1h-m). No obvious separation was seen in bacterial community by patient categories either at the genus (Fig. 1h-j) or species level (Fig. 1k-m). Yet, the distribution of samples in axis (second NMDS, first PCA and second PCoA at genus level; and first PCA at species level) differed when comparing IDH with control and ISH, suggesting a difference of β-diversity to some extent (β-diversity at genus level: Fig. 1h-j; PNMDS2 (Control vs. IDH) = 0.0414, PNMDS2 (ISH vs. IDH) = 0.0442; PPCA1 (Control vs. IDH) = 0.0204, PPCA1 (ISH vs. IDH) = 0.0211; PPCoA2 (Control vs. IDH) = 0.0237, PPCoA2 (ISH vs. IDH) = 0.0179; β-diversity at species level: Fig. 1l; PPCA1 (Control vs. IDH) = 0.0316, PPCA1 (ISH vs. IDH) = 0.0333).

In addition, to identify how heterogeneity among samples within the same category as other investigators previously proposed [26, 27], the dispersion of replicates within each group was evaluated for both bacterial diversity and patient phenotypic variability. A similar within β-diversity for bacterial community was detected (Fig. S2A). Instead, individuals in groups vary more in phenotypic variability (Fig. S2B).

Bacteria profiling and comparison of bacteria composition between groups

Next, taxonomic annotation and abundance profiles of bacteria composition were assessed among groups (Fig. 2). Our results revealed that regardless of the genus (Fig. 2a) or species (Fig. 2b) level, patients with IDH had a lower number of taxa annotated than those with ISH and controls, which is consistent with lower α-diversity obtained in Fig. 1. In addition, a relative abundance of the top ten most abundant genera, including Bacteroides, Prevotella, and Faecalibacterium, and the top 10 predominant species, including Faecalibacterium prausnitzii, Prevotella copri, Prevotella copri CAG:164 were assessed in each group (Fig. 2c-d) and in each sample (Fig. 2e-f). The global taxa composition differed across groups. Compared with ISH patients, those with IDH were characterized by a lower proportion of Alistipes and Faecalibacterium at the genus level and Faecalibacterium prausnitzii at the species level.

Fig. 2
figure 2

Genera and species annotated in the stool of control, ISH, and IDH. a, b Venn diagrams showing the number of genera and species annotated in groups. A total of 1167 genera and 3786 species were shared in control (con), ISH, and IDH groups. c, d Bar plots showed relative abundance and proportion of the top 10 genera and species in each group. Genera and species are differentiated by color. e, f Bar plots showed relative abundance and proportion of the top 10 genera and species in each sample. Genera and species are differentiated by color. Others in Fig. S2C and S2D indicates the sum of all the other genera/species expect the top 10 genera/species

Next, direct comparisons and differential analyses in bacteria abundance among these groups were performed at genus (Fig. 3) and species (Fig. 4) levels so as to identify a number of discriminatory taxa that were differentially abundant. Compared with healthy controls, intestinal bacteria of 100 genera (Fig. 3a) and 389 species (Fig. 4a) significantly differed in abundance from those with IDH. Furthermore, 43 genera (Fig. 3a) and 159 species (Fig. 4a) significantly differed between ISH and healthy controls. Of note, 26 shared taxa with 5 genera and 21 species had significantly different abundance in IDH and ISH, which strongly suggested that their modulation was related to both IDH and ISH (Fig. 3a and Fig. 4a). As compared to controls, the top 30 altered genera in IDH and in ISH are shown in Fig. 3b-c, respectively. Likewise, Fig. 4b-c show the top 30 out of the 389 differential species in IDH and from the 159 in ISH, respectively.

Fig. 3
figure 3

Significantly different Genera between control and ISH, control and IDH, respectively. a Venn diagrams demonstrating the genera number were significantly altered when comparing control (con) and IDH, control (con), and ISH. The overlap identified concurrently altered 5 genera; q<0.05 (P values were corrected using the Benjamini-Hochberg method, shown as q values), Wilcoxon rank-sum test. b Heat map showing the top 30 of the 100 differential genera between control and IDH. The abundance profiles were expressed by Z scores. The Z score was negative (shown in blue) when the row abundance was lower than the mean and is shown in red when the row abundance was higher than the mean. c Heat map for the top 30 of the 43 genera significantly shifted in all individuals with ISH compared to healthy controls. d Scatter plots indicating relative abundance of the 5 differential genera shared between control vs. IDH and control vs. ISH. In order to clearly show it, the relative abundance of genera was transformed into log10 values. Green, control; blue, ISH; and orange, IDH. The dots indicate individual values of the participants, and the horizontal lines from bottom to top represent 25th percentiles, medians, and 75th percentiles, respectively

Fig. 4
figure 4

Significantly distinct Species in IDH and ISH as compared with controls. a Venn diagrams demonstrating the number of differential species on comparing control (con) and IDH, control, and ISH, respectively. The overlap identified 21 concurrently varied species; q<0.05 (P values were corrected using the Benjamini-Hochberg method, shown as q values), Wilcoxon rank-sum test. b Heat map showing the top 30 of the 389 altered species in IDH. The abundance profiles were expressed by Z scores. The Z score was negative (shown in blue) when the row abundance was lower than the mean and is shown in red when the row abundance was higher than the mean. c Heat map showing the top 30 of the 159 different species in control and ISH. d Scatter plots show the relative abundance of the top 10 species shared between control vs. IDH, and control vs. ISH. In order to show it clearly, the relative abundance of species-level was transformed into log10 values. Green, control; blue, ISH; and orange, IDH. The dots indicate individual values of the participants, and the horizontal lines from bottom to top represent 25th percentiles, medians, and 75th percentiles, respectively

Intriguingly, it was worth noting that almost all these bacterial genera and species exhibited a deficiency in IDH, indicating a shift in taxonomic composition of fecal bacteria in patients with IDH. In addition, the relative abundance for the 5 shared genera (Fig. 3)d and top 10 of the 21 species (Fig. 4d) that synchronously shifted in ISH and IDH was further shown, particularly, the bacteria that were more abundant in ISH than controls and that showed higher enrichment in IDH, such as Rothia at the genus level and Rothia mucilaginosa at a species level. On the other hand, in the IDH group, we also found a much lower abundance for the taxa decreased in ISH, including Faecalicoccus and Acetohalobium at the genus level and Clostriduim sp. ASBs410, Bacteroides uniformis, and Bacteroides fluxus at the species level. Thus, more disturbed shifts might be present in populations suffering from IDH.

In addition, to further identify which microbes are correlated with which diseases cues or host metabolic change, the phenotypic variability of patients including sex, age, SBP, DBP and other demographic characteristics was correlated with the bacterial taxonomy by Spearmen’s correlation analysis. We found no matter at genus (Fig. S3) or species (Fig. S4) levels, the abundance of many bacteria were correlated with SBP and DBP. Furthermore, the possible correlation between these shared bacteria and BP levels was also estimated, revealing that among these 26 discriminatory taxa, the abundance of Rothia mucilaginosa was positively correlated with host SBP, while a negative correlation between Thermoanaerobacter indiensis, Clostriduim sp. ASBs410 (marginally significant) and DBP were detected (Fig. S5), suggesting they were likely to be explained by differences in SBP and DBP between groups.

Fecal bacterial function alteration in patients with ISH and IDH

Using the KEGG database, we assessed the potential gut bacterial functions across groups in the present study (Fig. 5). There was no clear separation in β-diversity between IDH and controls, and alteration in the intestinal bacterial functions in axes was not significant enough (Fig. 5a-c). Again, for the IDH patients, all the significantly altered KEGG modules were decreased in comparison with the control group, such as Sodium transport system, Acetyl-CoA pathway, Triacylglycerol biosynthesis, and Methanogenesis (Fig. 5d). The fecal bacteria function that varied in ISH patients was mainly characterized by nine increased functions, such as multidrug resistance and three reduced functions relevant to Biotin biosynthesis (Fig. 5e).

Fig. 5
figure 5

Fecal microbial gene functions were differently enriched between control and IDH, control and ISH. a-c NMDS, PCA, and PCoA plots based on the relative abundance of KEGG modules in groups. Marginal significant (P < 0.1) differences between control (con) and IDH were established at NMDS1, PCA2, and PCoA1; Kruskal-Wallis test. d Heat map showing the top 30 differential KEGG modules in control and IDH. The abundance profiles were expressed by Z scores. The Z score was negative (shown in blue) when the row abundance was lower than the mean and is shown in red when the row abundance was higher than the mean. e Heat map showing the 12 differential KEGG modules between control and ISH

In addition, the correlation between phenotypic variability of patients including sex, age, SBP, DBP and other demographic characteristics and functionality was assessed by Spearmen’s correlation analysis. We found that several functionality were correlated with SBP and DBP (Fig. S6). Meanwhile, in order to examine the functional redundancy, and the extent to which different species exhibit similar functions, we examined the relationship between species and functional diversity. Spearmen’s rank correlation analysis was performed to examine the relationship between number of species and functions is significant or not, and to identify the number of function increases, decreases or remains stable as species diversity increased [28]. We found strong positive relationships between number of species and number of KEGG modules that functional modules increases linearly with increasing species (F = 17.72, d.f. = 61, P < 0.001, linear R2 = 0.23; Fig. S7).

Discussion

The extent to which altered fecal bacteria is involved in ISH and IDH remains unclear to the present date. The results of the study have validated our previous hypothesis and provided some novel insights. By performing analysis based on shotgun metagenomic data, the present study discerned the gut bacterial signals in patients with ISH and IDH. Meanwhile, it compared taxonomic composition as well as potential functional modules with healthy controls. We found that compared with controls and ISH, patients with IDH showed a potential decrease in gene number, bacterial richness, and evenness. However, the bacterial alterations did not reach statistical significance in the Shannon index. Differences in the β-diversity were confirmed, and participants with IDH were clearly distinguished from those with ISH at the genus level. More importantly, a few bacterial taxa vanished in IDH, and apparent deficiency of various bacterial genera and species were detected in the taxonomic composition. Patients with IDH were observed to be enriched with Rothia mucilaginosa, along with reduced Clostridium spp.. Profound bacterial shifts in KEGG functions was noted in patients with IDH, with the sodium transport system and methanogenesis decreased.

The fecal bacteria is known to modify approximately 10% of the host’s transcriptome, especially those related to immunity and metabolism [29, 30], representing the potential signature of HTN. Indeed, a significant decrease in bacterial diversity and an increase of the Firmicutes/Bacteroidets ratio were described in spontaneously hypertensive rats [13]. In both pre-hypertensive and hypertensive populations, a dramatic decrease in microbial diversity, Prevotella-dominated gut enterotype, and disease-linked functions were also demonstrated [15]. To date, existing studies have compared fecal bacteria features between HTN and healthy controls; however, investigation of fecal bacteria in IDH, ISH, and control populations is still lacking. Even though a previous study [31] reported differences in the abundance of a few genera in the intestinal flora between patients with ISH or IDH with controls, from the independent cohort, as sequenced by 16S ribosomal RNA [31], the potential existence of profound bacterial variance of community structure in ISH and IDH patients has also caused increasing interest. Moreover, a comparison of the similarity and difference of ISH and IDH in bacterial composition and predicted functions of fecal bacteria still need to be further investigated. In the present study, we described the alterations in the intestinal bacteria of ISH and IDH by shotgun metagenomic sequencing data, which further allowed exploration in species-level and function-level.

Bacterial diversity has been regarded as a characteristic related to health status and multiple diseases. Previous studies have frequently reported that when the host is in pathological conditions, fecal bacteria’s diversity is often reduced or significantly altered [32]. Decrease of bacterial richness and diversity have been identified in several chronic diseases, such as HTN [13], Behcet’s disease [33], multiple sclerosis [34], and nonalcoholic fatty liver disease [35] by metagenomic shotgun sequencing, while in the current study, patients with IDH exhibit much deficient genes and a trend toward lower Shannon index, Chao richness and Pielou evenness compared with healthy controls and patients with ISH. Still, some indexes did not reach a statistical significance in the Shannon index. In addition, the β-diversity of participants with IDH differed from controls, as well as from patients with ISH. Based on the aforementioned findings, we concluded that the fecal bacteria structure in IDH patients was more severely perturbed.

Moreover, the variation in taxonomic composition was focused, and specific bacteria were distinct in ISH and IDH. It is noteworthy that the number of genus and species annotated in IDH patients was significantly lower compared to ISH patients, with 50 genera and 200 species completely missing, which is further consistent with reduced α-diversity in IDH. Additionally, compared to healthy controls, we identified that genus Rothia and Rothia mucilaginosa were enriched in both ISH and IDH. At the same time, Clostridium spp. was depleted. Rothia mucilaginosa, is a member of the oropharyngeal microbiota and is often detected in the upper respiratory tract [36]. Compelling evidence indicated that certain bacteria can be disseminated from one site to the others and cause systemic diseases [37]. In this regard, numerous studies showed that oral microbes could spread through the body, which have been found in various diseases, such as cardiovascular diseases, adverse pregnancy outcomes, and rheumatoid arthritis [38, 39]. In present study, we found that Rothia mucilaginosa, a member of oropharyngeal microbiota often detected in the upper respiratory tract, was enriched in fecal samples from ISH or IDH, which might be due to a possible cavity translocation.

Meanwhile, along with altering the gut bacterial composition, we observed a shifts in bacterial gene functions. Compared with healthy controls, all the KEGG pathways were significantly altered and deficient in IDH patients, such as methanogenesis and sodium transport system. Biosynthesis of methane is believed to have a crucial role in suppressing the inflammatory response and ameliorating oxidative stress injury in various tissues and organs [40,41,42]. Also, an excessive increase in reactive oxygen production and inflammatory response is implicated in the pathogenesis of HTN [43, 44]. In addition, there is increasing confirmatory evidence that the impaired capacity of sodium transport has an essential role in the development of HTN [45]. Thus, the gastrointestinal dysfunction further verified fecal bacteria’s shifts in IDH and might possibly mediate pathological processes during IDH.

Xie et al. [31] have reported the differences of intestinal flora between ISH patients and controls, and IDH versus controls, respectively by 16S amplicon sequencing. Firstly, the different microbial OTUs in six cases of IDH and six matched cases of normal BP were analyzed. They found five different OTUs between groups, with Acidaminococcus and Megasphaera higher in IDH, while Christensenella, Lactobacillus, and Olsenellathose higher in individuals with normal BP. Furthermore, the different OTUs in 35 cases of ISH and 35 matched controls were also analyzed. There were 38 different OTUs between groups, such as Acetobacteroides, Aestuariispira, and Akkermansia etc. enriched in ISH patients, Aeriscardovia, Alistipes, and Bilophila etc. abundant in the normal BP group. While in present study, we identified the microbial composition and functional modules among healthy controls, ISH and IDH by metagenomic analysis. We found that compared with controls and ISH, IDH patients exhibited decreased gene number, microbial richness and evenness, although the alterations did not reach statistical significance in Shannon index. The taxonomic composition of ISH or IDH was different from that of healthy controls at both genus and species levels. Patients with IDH or ISH were confirmed to be enriched with Rothia mucilaginosa, along with reduced Clostridium sp. ASBs410. Lastly, the altered KEGG modules were significantly decreased in IDH compared with the controls, such as sodium transport system; while for ISH, functions relevant to biotin biosynthesis were decreased.

Multiple clinical and observational studies have demonstrated elevated SBP as a more potent predictor of adverse cardiovascular outcomes than DBP in the elderly [46] and emphasized the importance of treating ISH in the elderly so as to reduce the risk of cardiovascular disease events [47]. However, younger subjects should not be ignored as DBP is increasingly predominant and contributes to many adverse outcomes [47, 48]. Khattar et al [48] have compared the prognostic significance of 24-h intra-arterial ambulatory DBP vs. SBP in untreated subjects between aged < 60 and ≥ 60. It was quite interesting that when SBP and DBP were jointly included in the prediction model, DBP (P = 0.04) instead of SBP (P = 0.67) was related to composite cardiovascular events in the younger subjects. Our findings support a certain extent to previous notions from the aspect of intestinal flora that fecal bacteria of IDH is more disordered in the relatively young subjects.

Our study was the first study to assess the fecal bacteria profiles in participants with ISH or IDH by shoutgun metagenomic sequencing, which allowed to analyze in species-level and functional modules. Meanwhile, there are some limitations. On the one hand, the sample size in our study is relative small. Thus rarefaction curves were performed to evaluate whether the sample size is enough. We found that the curves tend to be flat, indicating that the sample size is sufficient and reasonable, and only a small number of new species characteristics would be produced by more sampling. Moreover, some other studies investigating the gut microbial compositions in diseases, such as type 2 diabetes and obesity, relatively low number of participants with 3 to 6 in each group was also applied [49, 50] Nevertheless, future studies with larger sample size and more patients to validate or check the quantify of these identified species/genus and the findings in the current study by qPCR are still needed. On the other hand, although the reconstructed genomes obtained from the metagenome have the advantage to have longer sequences and cover the entire bacterial 16S rRNA gene for a more detailed taxonomical analysis [12], DNA of dead bacteria could not be excluded from the fecal samples in the whole-genome shotgun sequencing. There is possibility that the DNA in fecal samples detected by DNA-sequencing methods might not be active bacteria. This might hide the real differences among active bacterial members of the different community. As the active bacterial gene expression is affected by gene transcription, metatranscriptomes based methods could be a better indicator of the actual bacteria in future studies [12]. Despite these limitations, the data from present study clearly demonstrate that the fecal bacterial profiles of participants with ISH or IDH.

In conclusion, we described G’s disordered profiles in subjects with IDH and ISH and identified further shifts in IDH. A set of altered bacteria shared between ISH and IDH may comprise features of HTN, which warrants further research. The suppression of specific fecal bacteria, such as translocation of oropharyngeal microbiota, might contribute to pathogenesis. Importantly, our findings also point towards the significance of early intervention for IDH.

Availability of data and materials

The data set supporting the results of this article is from the EMBL European Nucleotide Archive (ENA) under BioProject accession code PRJEB13870.

Abbreviations

HTN:

Hypertension

ISH:

Isolated systolic hypertension

IDH:

Isolated diastolic hypertension

SBP:

Systolic blood pressure

DBP:

Diastolic blood pressure

NMDS:

Nonmetric dimensional scaling

PCA:

Principal-component analysis

PCoA:

Principal coordinate analysis

KEGG:

Kyoto Encyclopedia of Genes and Genomes

Nor:

Healthy controls

References

  1. Coats A, Jain S. Protective effects of nebivolol from oxidative stress to prevent hypertension-related target organ damage. J Hum Hypertens. 2017;31(6):376–81. https://0-doi-org.brum.beds.ac.uk/10.1038/jhh.2017.8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. GBD Causes of Death Collaborators. Global, regional, and national age-sex specific mortality for 264 causes of death, 1980-2016: a systematic analysis for the global burden of disease study 2016. Lancet. 2017;390(10100):1151–210. https://0-doi-org.brum.beds.ac.uk/10.1016/S0140-6736(17)32152-9.

    Article  Google Scholar 

  3. Yang G, Wang Y, Zeng Y, Gao GF, Liang X, Zhou M, et al. Rapid health transition in China, 1990-2010: findings from the global burden of disease study 2010. Lancet. 2013;381(9882):1987–2015. https://0-doi-org.brum.beds.ac.uk/10.1016/S0140-6736(13)61097-1.

    Article  PubMed  PubMed Central  Google Scholar 

  4. GBD Causes of Death. Global, regional, and national under-5 mortality, adult mortality, age-specific mortality, and life expectancy, 1970-2016: a systematic analysis for the global burden of disease study 2016. Lancet. 2017;390(10100):1084–150. https://0-doi-org.brum.beds.ac.uk/10.1016/S0140-6736(17)31833-0.

    Article  Google Scholar 

  5. Wang Z, Chen Z, Zhang L, Wang X, Hao G, Zhang Z, et al. Status of hypertension in China: results from the China hypertension survey, 2012-2015. Circulation. 2018;137(22):2344–56. https://0-doi-org.brum.beds.ac.uk/10.1161/CIRCULATIONAHA.117.032380.

    Article  PubMed  Google Scholar 

  6. He FJ, Li J, Macgregor GA. Effect of longer-term modest salt reduction on blood pressure. Cochrane Database Syst Rev. 2013;(4):Cd004937. https://0-doi-org.brum.beds.ac.uk/10.1002/14651858.CD004937.pub2.

  7. Kato N, Takeuchi F, Tabara Y, Kelly TN, Go MJ, Sim X, et al. Meta-analysis of genome-wide association studies identifies common variants associated with blood pressure variation in east Asians. Nat Genet. 2011;43(6):531–8. https://0-doi-org.brum.beds.ac.uk/10.1038/ng.834.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Karlsson FH, FÃ¥k F, Nookaew I, Tremaroli V, Fagerberg B, Petranovic D, et al. Symptomatic atherosclerosis is associated with an altered gut metagenome. Nat Commun. 2012;3(1):1245. https://0-doi-org.brum.beds.ac.uk/10.1038/ncomms2266.

    Article  CAS  PubMed  Google Scholar 

  9. Tremaroli V, Bäckhed F. Functional interactions between the gut microbiota and host metabolism. Nature. 2012;489(7415):242–9. https://0-doi-org.brum.beds.ac.uk/10.1038/nature11552.

    Article  CAS  PubMed  Google Scholar 

  10. Yamashiro K, Tanaka R, Urabe T, Ueno Y, Yamashiro Y, Nomoto K, et al. Gut dysbiosis is associated with metabolism and systemic inflammation in patients with ischemic stroke. PLoS One. 2017;12(2):e0171521. https://0-doi-org.brum.beds.ac.uk/10.1371/journal.pone.0171521.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Sommer F, Anderson JM, Bharti R, Raes J, Rosenstiel P. The resilience of the intestinal microbiota influences health and disease. Nat Rev Microbiol. 2017;15(10):630–8. https://0-doi-org.brum.beds.ac.uk/10.1038/nrmicro.2017.58.

    Article  CAS  PubMed  Google Scholar 

  12. Moya A, Ferrer M. Functional redundancy-induced stability of gut microbiota subjected to disturbance. Trends Microbiol. 2016;24(5):402–13. https://0-doi-org.brum.beds.ac.uk/10.1016/j.tim.2016.02.002.

    Article  CAS  PubMed  Google Scholar 

  13. 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. https://0-doi-org.brum.beds.ac.uk/10.1161/HYPERTENSIONAHA.115.05315.

    Article  CAS  PubMed  Google Scholar 

  14. Mell B, Jala VR, Mathew AV, Byun J, Waghulde H, Zhang Y, et al. Evidence for a link between gut microbiota and hypertension in the Dahl rat. Physiol Genomics. 2015;47(6):187–97. https://0-doi-org.brum.beds.ac.uk/10.1152/physiolgenomics.00136.2014.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. 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. https://0-doi-org.brum.beds.ac.uk/10.1186/s40168-016-0222-x.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Franklin SS. Arterial stiffness and hypertension: a two-way street? Hypertension. 2005;45(3):349–51. https://0-doi-org.brum.beds.ac.uk/10.1161/01.HYP.0000157819.31611.87.

    Article  CAS  PubMed  Google Scholar 

  17. Chirinos JA, Segers P, Hughes T, Townsend R. Large-artery stiffness in health and disease: JACC state-of-the-art review. J Am Coll Cardiol. 2019;74(9):1237–63. https://0-doi-org.brum.beds.ac.uk/10.1016/j.jacc.2019.07.012.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Romero CA, Alfie J, Galarza C, Waisman G, Peixoto AJ, Tabares AH, et al. Hemodynamic circulatory patterns in young patients with predominantly diastolic hypertension. J Am Soc Hypertens. 2013;7(2):157–62. https://0-doi-org.brum.beds.ac.uk/10.1016/j.jash.2013.01.001.

    Article  PubMed  Google Scholar 

  19. Johnson HM, Bartels CM, Thorpe CT, Schumacher JR, Pandhi N, Smith MA. Differential diagnosis and treatment rates between systolic and diastolic hypertension in young adults: a multidisciplinary observational study. J Clin Hypertens (Greenwich). 2015;17(11):885–94. https://0-doi-org.brum.beds.ac.uk/10.1111/jch.12596.

    Article  Google Scholar 

  20. Zeller G, Tap J, Voigt AY, Sunagawa S, Kultima JR, Costea PI, et al. Potential of fecal microbiota for early-stage detection of colorectal cancer. Mol Syst Biol. 2014;10(11):766. https://0-doi-org.brum.beds.ac.uk/10.15252/msb.20145645.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Yuzefpolskaya M, Bohn B, Nasiri M, Zuver AM, Onat DD, Royzman EA, et al. Gut microbiota, endotoxemia, inflammation, and oxidative stress in patients with heart failure, left ventricular assist device, and transplant. J Heart Lung Transplant. 2020;39(9):880–90. https://0-doi-org.brum.beds.ac.uk/10.1016/j.healun.2020.02.004.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Buchfink B, Xie C, Huson DH. Fast and sensitive protein alignment using DIAMOND. Nat Methods. 2015;12(1):59–60. https://0-doi-org.brum.beds.ac.uk/10.1038/nmeth.3176.

    Article  CAS  PubMed  Google Scholar 

  23. Huson DH, Auch AF, Qi J, Schuster SC. MEGAN analysis of metagenomic data. Genome Res. 2007;17(3):377–86. https://0-doi-org.brum.beds.ac.uk/10.1101/gr.5969107.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, et al. KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008;36(Database issue):D480–4. https://0-doi-org.brum.beds.ac.uk/10.1093/nar/gkm882.

    Article  CAS  PubMed  Google Scholar 

  25. Bäckhed F, Ley RE, Sonnenburg JL, Peterson DA, Gordon JI. Host-bacterial mutualism in the human intestine. Science. 2005;307(5717):1915–20. https://0-doi-org.brum.beds.ac.uk/10.1126/science.1104816.

    Article  CAS  PubMed  Google Scholar 

  26. Zaneveld JR, McMinds R, Vega TR. Stress and stability: applying the Anna Karenina principle to animal microbiomes. Nat Microbiol. 2017;2(9):17121. https://0-doi-org.brum.beds.ac.uk/10.1038/nmicrobiol.2017.121.

    Article  CAS  PubMed  Google Scholar 

  27. Ma ZS. Testing the Anna Karenina principle in human microbiome-associated diseases. iScience. 2020;23(4):101007. https://0-doi-org.brum.beds.ac.uk/10.1016/j.isci.2020.101007.

    Article  CAS  PubMed  Google Scholar 

  28. Guillemot N, Kulbicki M, Chabanet P, Vigliola L. Functional redundancy patterns reveal non-random assembly rules in a species-rich marine assemblage. PLoS One. 2011;6(10):e26735. https://0-doi-org.brum.beds.ac.uk/10.1371/journal.pone.0026735.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Min YW, Rhee PL. The role of microbiota on the gut immunology. Clin Ther. 2015;37(5):968–75. https://0-doi-org.brum.beds.ac.uk/10.1016/j.clinthera.2015.03.009.

    Article  CAS  PubMed  Google Scholar 

  30. Sommer F, Nookaew I, Sommer N, Fogelstrand P, Backhed F. Site-specific programming of the host epithelial transcriptome by the gut microbiota. Genome Biol. 2015;16(1):62. https://0-doi-org.brum.beds.ac.uk/10.1186/s13059-015-0614-4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Dan X, Mushi Z, Baili W, Han L, Enqi W, Huanhu Z, et al. Differential analysis of hypertension-associated intestinal microbiota. Int J Med Sci. 2019;16(6):872–81. https://0-doi-org.brum.beds.ac.uk/10.7150/ijms.29322.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Dingeo G, Brito A, Samouda H, Iddir M, La Frano MR, Bohn T. Phytochemicals as modifiers of gut microbial communities. Food Funct. 2020;11(10):8444–71. https://0-doi-org.brum.beds.ac.uk/10.1039/d0fo01483d.

    Article  CAS  PubMed  Google Scholar 

  33. Ye Z, Zhang N, Wu C, Zhang X, Wang Q, Huang X, et al. A metagenomic study of the gut microbiome in Behcet's disease. Microbiome. 2018;6(1):135. https://0-doi-org.brum.beds.ac.uk/10.1186/s40168-018-0520-6.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Kishikawa T, Ogawa K, Motooka D, Hosokawa A, Kinoshita M, Suzuki K, et al. A Metagenome-wide association study of gut microbiome in patients with multiple sclerosis revealed novel disease pathology. Front Cell Infect Microbiol. 2020;10:585973. https://0-doi-org.brum.beds.ac.uk/10.3389/fcimb.2020.585973.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Schwimmer JB, Johnson JS, Angeles JE, Behling C, Belt PH, Borecki I, et al. Microbiome signatures associated with Steatohepatitis and moderate to severe fibrosis in children with nonalcoholic fatty liver disease. Gastroenterology. 2019;157(4):1109–22. https://0-doi-org.brum.beds.ac.uk/10.1053/j.gastro.2019.06.028.

    Article  CAS  PubMed  Google Scholar 

  36. Luque-Pérez S, Cobos-Carrascosa E, Guarino-Narváez J, Fernández-Puentes V, Eiros JM, Sánchez-Porto A. Rothia mucilaginosa bacteraemia in an immunocompetent paediatric patient: a new pathogen to take into account. A case report. Infez Med. 2017;25(4):371–3.

    PubMed  Google Scholar 

  37. Kitamoto S, Nagao-Kitamoto H, Hein R, Schmidt TM, Kamada N. The bacterial connection between the Oral cavity and the gut diseases. J Dent Res. 2020;99(9):1021–9. https://0-doi-org.brum.beds.ac.uk/10.1177/0022034520924633.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Hajishengallis G. Periodontitis: from microbial immune subversion to systemic inflammation. Nat Rev Immunol. 2015;15(1):30–44. https://0-doi-org.brum.beds.ac.uk/10.1038/nri3785.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Graves DT, Corrêa JD, Silva TA. The Oral microbiota is modified by systemic diseases. J Dent Res. 2019;98(2):148–56. https://0-doi-org.brum.beds.ac.uk/10.1177/0022034518805739.

    Article  CAS  PubMed  Google Scholar 

  40. Shen M, Fan D, Zang Y, Chen Y, Zhu K, Cai Z, et al. Neuroprotective effects of methane-rich saline on experimental acute carbon monoxide toxicity. J Neurol Sci. 2016;369:361–7. https://0-doi-org.brum.beds.ac.uk/10.1016/j.jns.2016.08.055.

    Article  CAS  PubMed  Google Scholar 

  41. Wu J, Wang R, Ye Z, Sun X, Chen Z, Xia F, et al. Protective effects of methane-rich saline on diabetic retinopathy via anti-inflammation in a streptozotocin-induced diabetic rat model. Biochem Biophys Res Commun. 2015;466(2):155–61. https://0-doi-org.brum.beds.ac.uk/10.1016/j.bbrc.2015.08.121.

    Article  CAS  PubMed  Google Scholar 

  42. Zhang X, Li N, Shao H, Meng Y, Wang L, Wu Q, et al. Methane limit LPS-induced NF-κB/MAPKs signal in macrophages and suppress immune response in mice by enhancing PI3K/AKT/GSK-3β-mediated IL-10 expression. Sci Rep. 2016;6(1):29359. https://0-doi-org.brum.beds.ac.uk/10.1038/srep29359.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Djordjević VB, Zvezdanović L, Cosić V. Oxidative stress in human diseases. Srp Arh Celok Lek. 2008;136 Suppl 2:158–65. https://0-doi-org.brum.beds.ac.uk/10.2298/SARH08S2158D.

    Article  PubMed  Google Scholar 

  44. Suzuki H, Suematsu M, Schmid-Schönbein GW. Microvascular oxidative stress, immune reaction and apoptosis in hypertensives. Clin Hemorheol Microcirc. 1999;21(3–4):161–8. https://0-doi-org.brum.beds.ac.uk/10.4049/jimmunol.178.3.1931.

    Article  CAS  PubMed  Google Scholar 

  45. Ayuzawa N, Fujita T. Activation of mineralocorticoid receptor in salt-sensitive hypertension. Curr Hypertens Rep. 2015;17(6):552. https://0-doi-org.brum.beds.ac.uk/10.1007/s11906-015-0552-2.

    Article  CAS  PubMed  Google Scholar 

  46. Bavishi C, Goel S, Messerli FH. Isolated systolic hypertension: an update after SPRINT. Am J Med. 2016;129(12):1251–8. https://0-doi-org.brum.beds.ac.uk/10.1016/j.amjmed.2016.08.032.

    Article  PubMed  Google Scholar 

  47. Staessen JA, Fagard R, Thijs L, Celis H, Arabidze GG, Birkenhäger WH, et al. Randomised double-blind comparison of placebo and active treatment for older patients with isolated systolic hypertension. The systolic hypertension in Europe (Syst-Eur) trial investigators. Lancet. 1997;350(9080):757–64. https://0-doi-org.brum.beds.ac.uk/10.1016/s0140-6736(97)05381-6.

    Article  CAS  PubMed  Google Scholar 

  48. Khattar RS, Swales JD, Dore C, Senior R, Lahiri A. Effect of aging on the prognostic significance of ambulatory systolic, diastolic, and pulse pressure in essential hypertension. Circulation. 2001;104(7):783–9. https://0-doi-org.brum.beds.ac.uk/10.1161/hc3201.094227.

    Article  CAS  PubMed  Google Scholar 

  49. Wang TY, Zhang XQ, Chen AL, Zhang J, Lv BH, Ma MH, et al. A comparative study of microbial community and functions of type 2 diabetes mellitus patients with obesity and healthy people. Appl Microbiol Biotechnol. 2020;104(16):7143–53. https://0-doi-org.brum.beds.ac.uk/10.1007/s00253-020-10689-7.

    Article  CAS  PubMed  Google Scholar 

  50. Li R, Huang X, Liang X, Su M, Lai KP, Chen J. Integrated omics analysis reveals the alteration of gut microbe-metabolites in obese adults. Brief Bioinform. 2020. https://0-doi-org.brum.beds.ac.uk/10.1093/bib/bbaa165.

Download references

Acknowledgments

Not applicable.

Funding

This work was supported by the National Natural Science Foundation of China (81670214, 81870308, and 81970271), the Beijing Natural Science Foundation (7172080), the National Nature Cultivation Fund of Beijing Chaoyang Hospital (CHPY202050).

Author information

Authors and Affiliations

Authors

Contributions

J. Li, X. Yang, P. Wang conceived the study, directed the project, designed the experiments. J. Li, P. Wang and Y. Dong interpreted the results, and wrote the manuscript. J. Jiao, K. Zuo, C. Han recruited and collected the clinical details from the subjects. P. Wang, Y. Dong, and J. Li analyzed the data. X. Yang, P. Wang, Y. Dong, and J. Li revised the manuscript. P. Wang and Y. Dong contributed equally to this work and are joint-first authors. The authors read and approved the final manuscript.

Corresponding authors

Correspondence to Xinchun Yang or Jing Li.

Ethics declarations

Ethics approval and consent to participate

This article does not contain any studies with animals performed by any of the authors. All procedures performed in studies involving human participants were in accordance with the ethical standards of the institutional and/or national research committee and with the 1964 Helsinki declaration and its later amendments or comparable ethical standards. The study was approved by local ethics committees (Kailuan General Hospital, Beijing Chaoyang Hospital, and Beijing Fuwai Hospital) and informed consent was obtained from all subjects.

Consent for publication

Not applicable.

Competing interests

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.

Additional information

Publisher’s Note

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

Supplementary Information

Additional file 1: Table S1.

Data production of fecal samples in control, ISH and IDH.

Additional file 2: Figure S1.

Rarefaction curves. (A) Rarefaction curves by randomly selecting a certain number of individuals from the samples and counting the number of species represented by these individuals in control, ISH and IDH. The curves tend to be flat, indicating the sample size is sufficient and reasonable. (B) Rarefaction curves by gradually expanding the sequencing depth of random sampling. The curves approached saturation as the sample sequencing depth increases, and thus the amount of sequencing data is sufficient and stable.

Additional file 3: Figure S2.

The dispersion of replicates with each group for both bacterial diversity and patient phenotypic variability. (A) PCA plots based on the species level in groups. (B) PCA plots based on the patient phenotypic variability within each group. A similar with beta-diversity for bacterial community was detected. Instead, individuals in groups vary more in phenotypic variability. Circles in green indicate samples from control, squares in blue indicate samples from ISH, rhombus in orange represent individuals from IDH.

Additional file 4: Figure S3.

The correlation between phenotypic variability and bacterial diversity by genus level. (A) Heat map showing the correlation between phenotypic variability and bacterial diversity at genus level. Indices with |correlation| ≥ 0.2 were shown; Blue, negative correlation; red, positive correlation; *, P < 0.05; **, P < 0.01.

Additional file 5: Figure S4.

The correlation between phenotypic variability and bacterial diversity by species level. (A) Heat map showing the correlation between phenotypic variability and bacterial diversity at species level. Indices with |correlation| ≥ 0.2 were shown; Blue, negative correlation; red, positive correlation; *, P < 0.05; **, P < 0.01.

Additional file 6: Figure S5.

The correlation of shared differential bacteria in ISH and IDH with BP. (A) Correlation plots between shared differential bacteria on genus and species level (5 at genus level and 21 at species level) and BP (SBP and DBP). The correlation coefficient is expressed in different colors and sizes. R-value < 0 (negative correlation) are marked with blue circle; r-value > 0 (positive correlation) are in red circle. *P < 0.05, +P < 0.1. (B) A significant negative correlation between Thermoanaerobacter indiensis and DBP (r = − 0.253, P = 0.04568) is shown. (C-F) Significant correlation between zeta proteobacterium SCGC AB-604-B04 (r = − 0.255, P = 0.04339), Morganella (r = − 0.266, P = 0.0349), Rothia and SBP (r = 0.276, P = 0.0286), Rothia mucilaginosa (r = 0.283, P = 0.0246) and SBP is shown.

Additional file 7: Figure S6.

The correlation between phenotypic variability and functional modules. (A) Heat map showing the correlation between phenotypic variability and functional modules. Indices with |correlation| ≥ 0.2 were shown; Blue, negative correlation; red, positive correlation; *, P < 0.05; **, P < 0.01.

Additional file 8: Figure S7.

Functional redundancy was assessed by Spearman’s rank correlation. (A) The relationship between number of species and functions by Spearmen’s rank correlation. To examine the functional redundancy, and the extent to which different species exhibit similar functions, we examined the relationship between species and functional diversity by Spearman’s rank correlation analysis. A strong positive relationships between number of species and number of KEGG modules that functional modules increases linearly with increasing species (F = 17.72, d.f. = 61, P < 0.001, linear R2 = 0.23).

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wang, P., Dong, Y., Zuo, K. et al. Characteristics and variation of fecal bacterial communities and functions in isolated systolic and diastolic hypertensive patients. BMC Microbiol 21, 128 (2021). https://0-doi-org.brum.beds.ac.uk/10.1186/s12866-021-02195-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12866-021-02195-1

Keywords