Gut microbiota dysbiosis contributes to the development of hypertension
- Jing Li†1, 2, 3,
- Fangqing Zhao†4,
- Yidan Wang†1,
- Junru Chen†5,
- Jie Tao†6,
- Gang Tian7,
- Shouling Wu8,
- Wenbin Liu5,
- Qinghua Cui9,
- Bin Geng1,
- Weili Zhang1,
- Ryan Weldon10,
- Kelda Auguste10,
- Lei Yang11,
- Xiaoyan Liu11,
- Li Chen10, 12, 13,
- Xinchun Yang2, 3Email author,
- Baoli Zhu14, 15Email author and
- Jun Cai1Email author
© The Author(s). 2017
Received: 25 July 2016
Accepted: 13 December 2016
Published: 1 February 2017
Recently, the potential role of gut microbiome in metabolic diseases has been revealed, especially in cardiovascular diseases. Hypertension is one of the most prevalent cardiovascular diseases worldwide, yet whether gut microbiota dysbiosis participates in the development of hypertension remains largely unknown. To investigate this issue, we carried out comprehensive metagenomic and metabolomic analyses in a cohort of 41 healthy controls, 56 subjects with pre-hypertension, 99 individuals with primary hypertension, and performed fecal microbiota transplantation from patients to germ-free mice.
Compared to the healthy controls, we found dramatically decreased microbial richness and diversity, Prevotella-dominated gut enterotype, distinct metagenomic composition with reduced bacteria associated with healthy status and overgrowth of bacteria such as Prevotella and Klebsiella, and disease-linked microbial function in both pre-hypertensive and hypertensive populations. Unexpectedly, the microbiome characteristic in pre-hypertension group was quite similar to that in hypertension. The metabolism changes of host with pre-hypertension or hypertension were identified to be closely linked to gut microbiome dysbiosis. And a disease classifier based on microbiota and metabolites was constructed to discriminate pre-hypertensive and hypertensive individuals from controls accurately. Furthermore, by fecal transplantation from hypertensive human donors to germ-free mice, elevated blood pressure was observed to be transferrable through microbiota, and the direct influence of gut microbiota on blood pressure of the host was demonstrated.
Overall, our results describe a novel causal role of aberrant gut microbiota in contributing to the pathogenesis of hypertension. And the significance of early intervention for pre-hypertension was emphasized.
KeywordsHypertension Pre-hypertension Gut microbiota Metabolism Fecal transplant
In recent decades, the potential role of the gut microbiome in altering health status of the hosts has drawn considerable attention. Emerging evidence suggests a link between gut microbiome and various diseases, including colorectal cancer, liver cirrhosis, arthritis, type 2 diabetes, and atherosclerosis [1–5]. A number of microbial biomarkers specific to these diseases have been discovered, and fecal microbiome-targeted strategies are recommended to be a powerful tool for early diagnosis and treatment of different diseases.
More importantly, by fecal transfer experiment and gut microbiota (GM) remodeling, intestinal microbiome has been further indicated to conduce to the pathogenesis of multiple diseases such as obesity, depressive disorder, chronic ileal inflammation, liver diseases, and atherosclerosis [6–12]. Specific mechanisms underlying the causal function of GM have been revealed. For example, the metabolism by intestinal microbiota of dietary L-carnitine, a nutrient in red meat, was demonstrated to promote atherosclerosis and lead to cardiovascular disease risk via producing trimethylamine and trimethylamine-N-oxide . Targeting gut microbial production of trimethylamine specifically and non-lethal microbial inhibitors were confirmed to relieve diet-induced atherosclerotic lesion development . Thus GM may serve as a potential therapeutic approach for the treatment of cardiovascular and metabolic diseases.
Hypertension (HTN) has become a global public health concern and a major risk factor for cardiovascular, cerebrovascular, and kidney diseases [14, 15]. It is believed that the etiology of HTN depends on the complex interplay of both genetic and environmental factors [16, 17], and the precise cause of this morbidity has not been elucidated to date. It has been suggested that the germ-free (GF) mice, in which the intestinal bacteria is completely absent, present relatively lower blood pressure (BP) when compared to conventional mice . And therefore we suspected that GM might have the potential to regulate BP.
Most recently, many lines of seminal evidence, which for the first time demonstrate that aberrant gut microbial community are linked to BP changes of the host, support this hypothesis. For example, disordered GM as a result of decreased microbial richness, diversity, evenness, and increased Firmicutes/Bacteroidetes ratio was reported in hypertensive animals and seven HTN patients, as sequenced by 16S ribosomal RNA . In Dahl rats, distinct metagenomic composition have been revealed between salt-sensitive and salt-resistant strains, and the GM of salt-sensitive rats was suggested to be in a symbiotic relationship with the host . In addition, by rat models of HTN and meta-analyses in randomized human clinical trials, investigators have revealed that administration of probiotics can reduce BP [21, 22]. This drove us to speculate that the alteration in GM by probiotic use may lead to BP changes. Furthermore, it has been proved that transplantation of cecal contents from hypertensive obstructive sleep apnea rats on high-fat diet into recipient rats on normal chow diet lead to higher BP levels, and a major contributor to the gut dysbiosis of obstructive sleep apnea-induced HTN is high-fat diet . These studies have emphasized a strong correlation between gut dysbiosis and HTN, and further implied the significance of GM in BP regulation, yet animal models could not perfectly substitute human disease, and the sample size of human participants for microbial analysis was quite limited.
In consideration of the BP levels being classified into optimal, pre-hypertension (pHTN), and HTN according to the most recent clinical guidelines , it remains obscure how exactly the composition of gut microbes and the products of microbial fermentation change in human patients with HTN, especially in pHTN populations. In addition, decisive evidence is still needed to determine whether gut dysbiosis is a consequence or an important causal factor for the pathogenesis of HTN. Fecal transplantation from human samples into GF mice is required to uncover the involvement of GM dysbiosis in pathophysiology of HTN. Collectively, these key issues are the major goal of the present study.
To address the questions above, we performed deep metagenomic sequencing of stool samples from 196 participants of healthy control, pHTN, and HTN; took metabolomic analyses of their metabolic profiles, further constructed a disease classifier for pHTN and HTN based on GM and metabolites; and demonstrated the crucial role of disordered GM in triggering thigh BP by human fecal microbiota transplantation into GF mice.
GM diversity and enterotype in pHTN and HTN
To identify whether gut microbial changes are associated with HTN, we performed shotgun metagenomic sequencing of fecal samples from a cohort of 196 Chinese individuals. The cohort consisted of 41 healthy controls, 56 subjects with pHTN, and 99 patients with primary HTN. All the participants were from a cohort study among employees of the Kailuan Group Corporation. The Kailuan study is a prospective cohort study focusing on the Kailuan community in Tangshan, a large modern city in northern China. All the subjects in the hypertension group were newly diagnosed hypertensive patients prior to antihypertensive treatment. Patients suffering from cancer, heart failure, renal failure, smoking, stroke, peripheral artery disease, and chronic inflammatory disease were all excluded. Drugs including statins, aspirin, insulin, metformin, nifedipine, and metoprolol were not used on the patients, and other drug consumption was not compared because the sample size was quite small. Individuals were also excluded if they had received antibiotics or probiotics within the last 8 weeks. Other than SBP and DBP, there was no significant difference in other clinical parameters among groups, except for fasting blood glucose level (FBG) (P = 0.026, C vs H; Kruskal-Wallis test, Additional file 1: Table S1). Bacterial DNA was extracted from stool samples, sequenced on the Illumina platform, and a total of 1211 Gb 125-bp paired-end reads were generated, with an average of 6.18 ± 1.43 (s.d.) million reads per sample (Additional file 2: Table S2). For each sample, a majority of high-quality sequencing reads (83.74–97.24%) were de novo assembled into long contigs or scaffolds, which were used for gene prediction, taxonomic classification, and functional annotation.
To explore the difference between the microbial communities at different stages of HTN, enterotypes were identified based on the abundance of genera using Partitioning Around Medoid (PAM) clustering method. The optimal number of enterotypes was two as indicated by Calinski-Harabasz (CH) index (Additional file 3: Figure S1). Then Principal Coordinate Analysis (PCoA) using Jensen-Shannon distance was performed to cluster the 196 samples into two distinct enterotypes (Fig. 1d). Prevotella was the most enriched genus in enterotype 1; Bacteroides was the most enriched genus in enterotype 2 (P = 6.31e−31 and P = 2.09e−15, respectively; Wilcoxon rank sum test, Fig. 1e). Both contributors in the two enterotypes have been reported in European and Chinese populations before [2, 3]. There was a higher percentage of pre-hypertensive and hypertensive patients distributed in enterotype 1 (48.21% for pHTN, and 45.45% for HTN), while more healthy controls (73.17%) were found in enterotype 2 (P = 0.02, C vs P; P = 0.03, C vs H; Fisher’s exact test; Fig. 1f). These findings suggest that enterotype 2 may represent a GM community structure associated with healthy control, while enterotype 1 may be associated with pHTN and HTN.
Considering the higher percentage of HTN patients in enterotype 1, we clustered the genera in this enterotype and further explored the microbial co-occurrence network by Spearman’s correlation. There was a positively interacted network constituted by 12 genera, which were negatively correlated with Prevotella, the core genus in this enterotype (Additional file 4: Figure S2a). All these genera were decreased in enterotype 1 as compared with enterotype 2 (Additional file 4: Figure S2b). Eight out of them were directly linked to Prevotella, while the other four, including Oscillibacter, Faecalibacterium, Butyrivibrio, and Roseburia, were indirectly linked to Prevotella. These findings highlighted the possibility of Prevotella as a key genus associated with pHTN and HTN. The difference in gut enterotype distribution revealed profound changes of the intestinal microbiome structure in both pHTN and HTN, implying the significance of gut microbes in the development of HTN.
pHTN and HTN-associated genera in GM
By contrast, Faecalibacterium, Oscillibacter, Roseburia, Bifidobacterium, Coprococcus, and Butyrivibrio, which were enriched in healthy controls, declined in pHTN and HTN patients (Fig. 2b). Our observations were consistent with the genera negatively correlated with Prevotella in the network of enterotype 1 (Additional file 4: Figure S2), and these bacteria are known to be essential for healthy status. For example, reduced levels of Faecalibacterium and Roseburia in the intestines are associated with Crohn’s disease and ulcerative colitis [30, 31]. Both bacteria are crucial for butyric acid production [30, 32]. Moreover, Bifidobacterium is an important probiotic necessary to intestinal microbial homeostasis, gut barrier, and lipopolysaccharide (LPS) reduction .
The divergence of GM composition in each sample was assessed to explore the correlation of microbial abundance with body mass index (BMI), age, and gender (Additional file 5: Figure S3). Although the gender ratio is discrepant among groups (Additional file 1: Table S1), we found no remarkable regularity of bacterial abundance based on BMI, age or gender.
To further validate the bacterial alterations in HTN, an independent metagenomic analysis was performed using the sequencing data generated from a previous study of type 2 diabetes . From a total of 174 non-diabetic controls in the study, normotensive controls with SBP ≤125 mmHg or DBP ≤80 mmHg were enrolled, and HTN were elected with the inclusion criteria of SBP ≥150 mmHg or DBP ≥100 mmHg. The FBG levels between normotensive controls and HTN were similar. Finally, six subjects (HTNs, n = 3; normotensive controls, n = 3) were included in our analysis (Additional file 2: Table S6). As expected, the microbial diversity was decreased in HTN (Additional file 6: Figure S4a), and there were at least 20 genera showing consistent trends with our findings, including decreased Butyrivibrio, Clostridium, Faecalibacterium, Enterococcus, Roseburia, Blautia, Oscillbacter, and elevated Klebsiella, Prevotella, and Desulfovibrio (Additional file 6: Figure S4b, Additional file 2: Table S7).
Collectively, these results supported our hypothesis that bacteria associated with healthy status were reduced in patients with HTN. This phenomenon together with the overgrowth of bacteria such as Prevotella and Klebsiella may play important role in the pathology of HTN.
Co-abundance groups enriched in pHTN and HTN
Firstly, for each gene, an OR score was calculated according to the abundance of each gene. Then, for the comparative analysis between control and HTN samples, the HTN-associated genes were classified as HTN-enriched (OR >2) or HTN-depleted (OR <0.5) as previously described . When calculating HTN-associated ORs, samples of pHTN were excluded, and samples labeled as HTN were excluded as well when calculating pHTN-associated ORs. A total of 1,120,526 genes significantly different in relative abundance across groups were identified (Additional file 7: Table S8). Secondly, we clustered genes by a rather high threshold (Spearman’s correlation coefficient ≥0.7) according to previous methods [4, 35]. Spearman’s correlation coefficient was analyzed by R. The cluster groups with at least 50 genes were defined as co-abundance groups (CAGs) , and used for further analysis . One thousand ninety-nine distinct CAGs were obtained (Additional file 2: Table S9–S11 and Additional file 8: Figure S5a). Seven hundred fourteen CAGs were assigned to known bacterial genera based on the tracer genes, with at least 80% of the genes mapped to the reference genome at an identity higher than 85% (Additional file 8: Figure S5b).
To further examine the relationship between clinical indices and CAGs of GM, physiological parameters of SBP, DBP, BMI, FBG, total cholesterol (TC), triglyceride (TG), and low-density lipoprotein (LDL) were included in a Spearman’s correlation analysis. We observed that SBP and DBP could negatively influence the CAGs enriched in the control group, such as Firmicutes and Roseburia, and positively interacted with Prevotella and Desulfovibrio, which were abundant in pHTN and HTN (Additional file 9: Figure S6). Whereas, both TC and TG were negatively correlated with Enterobacter, that was enriched in pHTN and HTN groups. Altogether, these results indicated that the bacterial communities in individuals with pHTN and HTN are similar, and the collective effect of these bacteria may account for intestinal dysbiosis in HTN.
Functional alteration in GM of pHTN and HTN
Metabolic profiling of GM in pHTN and HTN
The compositional changes in patients involved 167 analytes that were significantly different between pHTN and control, and 215 analytes altered in HTN (Fig. 5c). There were 26 metabolites which were obviously different in both pHTN and HTN groups as compared to the control (Additional file 2: Table S16). Notably, these metabolites exhibited statistically analogous profiles of alterations in pHTN and HTN, which was consistent with our findings based on gut microbiome (Fig. 5d). Endogenous compounds whose levels significantly decreased in pHTN and HTN include phosphatidylserine (PS), 3,4,5-trimethoxycinnamic acid, lysophosphatidylcholine (LysoPC), S-carboxymethyl-l-cysteine, and lysophosphatidylethanolamine (LysoPE). 3,4,5-Trimethoxycinnamic acid is capable to protect against inflammatory diseases through suppressing cell adhesion molecules in vascular endothelial cells . Also S-Carboxymethyl-l-cysteine exerts anti-inflammatory properties . These observed downregulations could promote the inflammatory environment associated with HTN. On the other hand, endogenous compounds whose levels significantly increased in pHTN and HTN include metabolites such as Nα-acetyl-l-arginine, stearic acid, phosphatidic acid (PA), and glucoside. Elevated levels of Nα-acetyl-l-arginine and stearic acid have been previously observed in uremia and spontaneously hypertensive rats [45, 46]. These compounds may represent possible markers for the development of HTN and might be derived from gut microflora or their fermented products. To explore this idea, the relationship between 26 representative metabolites and the 44 most different genera was examined by correlation analysis (Fig. 5e). Control-enriched trichloroethanol glucuronide was positively correlated with Bifidobacterium and Akkermansia, but negatively linked to Prevotella. Conversely, there was a positive association between 9,10-dichloro-octadecanoic acid (stearic acid) and microflora including Klebsiella, Prevotella, and Enterbacter, which were all overrepresented in HTN. It was accordant that both Bifidobacterium and Roseburia negatively interacted with 9,10-dichloro-octadecanoic acid, which was hence considered as an important GM-influenced metabolic product in HTN. Thus the distinguished metabolic profiling in HTN was closely connected to intestinal microflora variation, although whether these metabolic products were directly metabolized by the intestinal microorganisms remained to be explored.
Identification of pHTN and HTN basing on gut microbiome
To illustrate the microbial and metabolic signature of pHTN and HTN, and further exploit the potential of gut microbiome and metabolites in HTN identification, random forest disease classifier using explanatory variables of CAGs, metabolites, and species abundances were performed. Tenfold cross-validation was repeated for five times and the receiver operating characteristic (ROC) curves for classifying pHTN and HTN patients from controls were developed.
We also investigated the utility of the classifier based on microbial CAGs + species. Consistently, the AUC for identifying pHTN and HTN from the controls was 0.67 (95% CI, 0.39–0.95) and 0.81 (95% CI, 0.53–1), respectively, and the performance on pHTN and HTN individuals was not as satisfactory (AUC, 0.47; 95% CI, 0.19–0.75; Additional file 10: Figure S7a). For HTN classification, CAGs and species taxonomic annotated to Prevotella, including Prevotella sp. CAG:5226.CAG-377, Prevotella bivia, and CAG-184 were typically important (Additional file 10: Figure S7b). Overall, the pHTN- and HTN-associated microbial and metabolic features captured by the classifier offered further evidence for dysbiotic gut microbiome and highlighted great potential ability for detection of pHTN and HTN populations by GM and metabolites.
High BP is transferrable by fecal transplant
Previous studies have revealed that antibiotics and probiotics are potential treatment modalities for BP in both animal models and clinical trials [19, 21, 22, 47]. We speculated that the alterations in GM under pro/antibiotic use may be associated with BP changes. There is evidence that Dahl salt-sensitive rats transplanted with salt-resistant rat microbiota have further exacerbated BP, which indicate that the microbiota resident within the cecum of the Dahl salt-sensitive rat, but not the salt-resistant rat, are in a symbiotic relationship with the host . Thus the differences between Dahl salt-sensitive rats and the salt-resistant rats are highlighted. Investigators have also proved that transplantation of cecal contents from hypertensive obstructive sleep apnea rats on high-fat diet into the same obstructive sleep apnea recipient rats on normal chow diet lead to higher BP similar to the donors . In this study, it seems that a major contributor to the gut dysbiosis of HTN is a high-fat diet. Therefore, direct studies testing if microbial transplantation can transmit changes in BP from hypertensive donors to recipients are still lacking. To further demonstrate whether alterations of GM are a causal factor for the progression of HTN in vivo, fecal bacteria from hypertensive patients were transplanted to GF mice in the present work.
At 10 weeks post-transplantation, BP of recipient mice in HTN and control group was measured by the tail-cuff method. Notably, the HTN mice exhibited significantly higher SBP, DBP, and mean blood pressure (MBP) as compared to controls (P < 0.05), as well as elevated heart rate (P = 0.11) (Fig. 7f). Early studies have shown that when compared to conventional controls, GF rats possess significantly lower cardiac output, relatively diminished regional blood flow, lower level of systemic BP response after blood loss, and hypotonic microvasculature , which might lead to a low systolic BP in the recipient mice. These findings provided novel and direct evidence that GM could influence the BP of host directly. Therefore, changes in the GM might be the mechanism underlying the effect of antibiotics and probiotics on BP control. As the number of donors for transplantation is limited, larger number of fecal transplants from hypertensive, pre-hypertensive, and normotensive control participants should be carried out in the future to further establish the magnitude of BP changes.
To date, there are limited studies indicating a direct association between GM and HTN, especially in human disease. Several important gaps in knowledge of gut and BP remain unexplored, and critical issues should be addressed, such as the microbial profiles of HTN populations in clinical trials, the metabolites signature profiles, the microbial biomarkers for early detection of HTN, and fecal transplantation to make clear the causal relationship between gut dysbiosis and HTN. To make up for these blanks, we applied a strategy based on metagenomic and metabolomic analyses, coupled with GM transplantation. We sequenced the total bacteria DNA of stool samples from a cohort of 196 Chinese individuals and supplemented this analysis with an additional validation cohort. All the individuals in the present study are from a cohort study among employees of the Kailuan Group Corporation. The Kailuan study is a prospective cohort study focusing on the Kailuan community in Tangshan, a large modern city in northern China. As the subjects were from a relatively concentrated environment, the differences in the diets were relatively small. All the subjects in the hypertension group were newly diagnosed hypertensive patients prior to antihypertensive treatment. Patients suffering from cancer, heart failure, renal failure, smoking, stroke, peripheral artery disease, and chronic inflammatory disease were all excluded. Drugs including statins, aspirin, insulin, metformin, nifedipine, and metoprolol were not used on the patients, and other drug consumption was not compared because the sample size was quite small. Hence, it is not likely that the medication use directly influenced the gut metagenome and metabolites, as there was no significant difference in the drugs consumed by these subjects. Our results demonstrate that decreased diversity, altered enterotype distribution, and variation in bacteria populations were associated with both pHTN and HTN. The bacterial metabolic functions and GM-related metabolites in pre-hypertensive and hypertensive adults were closely linked to inflammatory state. Particularly, both pHTN and HTN individuals could be accurately distinguished from the controls by variables of CAGs and metabolites. And most importantly, the direct impact of GM composition on regulating BP was evaluated using an in vivo model of GF mice colonized with human intestinal microbiota. Bacteria such as Prevotella, Klebsiella, Enterobacter, and Fusobacterium are potential candidates for further bacteria transfer experiments to explore the precise mechanisms underlying the effect of GM in BP regulation. Our work provides the first direct evidence that highlights the pivotal role of dysbiotic gut microbiome as an important pathogenic factor for the high BP of the host. Thus GM modulation should be considered during antihypertensive treatment.
Researchers previously suggested that the intestinal bacterium Prevotella copri thrives in a pro-inflammatory environment of rheumatoid arthritis [3, 49]. The superoxide reductase and phosphoadenosine phosphosulphate reductase encoded by Prevotella copri may favor the development of inflammation. In their further demonstration, colonization with Prevotella copri enhances body weight loss and exacerbates epithelial inflammation in colitis mouse model . Interestingly, as shown by our data, the enterotype dominated by Prevotella was enriched with pHTN and HTN populations. Moreover, Prevotella was overrepresented in individuals with pHTN and HTN. And stearic acid, an important metabolite in HTN, was positively linked to Prevotella. Furthermore, CAGs and species taxonomic annotated to Prevotella were the common microbial characteristics of pHTN and HTN, and contributed a lot to classification of HTN. Thus Prevotella may play an essential role in HTN, probably by triggering the inflammatory response. Our findings have consolidated the potential of Prevotella in the pathogenesis of diseases, and call for further exploring whether Prevotella is a causal conducer to inflammation and HTN.
Concomitant with the alteration of gut microbial composition, we observed a dysbiosis in bacterial gene functions. The metagenome of HTN patients were depleted in genes associated with biosynthesis and transport of amino acid, such as lysine, histidine, leucine, and serine, which are essential for human health. Functional annotation also indicated a decline of modules for fatty acid utilization and saccharide transport, suggesting an impaired capacity of energy production. In agreement with previous studies showing a dearth of microbial functions for purine metabolism in arthritis , a significant decrease in purine-metabolizing enzymes was identified to be related to HTN. Indeed, these metabolic functions are quite necessary for healthy populations [3–5, 39, 40]. In contrast, the enrichment of the modules for LPS biosynthesis and export in patients hints at a potential role of GM in causing low-grade inflammation. Inflammation due to immune response triggered by LPS is the cardinal feature of the pathogenesis of gram-negative bacteria, such as Prevotella and Klebsiella [50, 51], and has been identified as an important contributor to the pathogenesis of HTN. Our findings raise the possibility that the low-grade inflammation and increase of gram-negative bacteria, especially Prevotella and Klebsiella, are likely responsible for HTN pathology. Thus, our analysis of bacterial gene functions indicates that functional dysbiosis may contribute to the susceptibility to HTN, and overproduction of LPS by gut bacteria seems to be directly linked to HTN development, whereas amino acid biosynthesis, fatty acid utilization, and purine metabolism by bacteria might have a role in HTN prevention.
Actually, in GF mice, the tail-cuff method has been used for assessment of BP in a recent report, suggesting the methodology is acceptable . In our study, the tail-cuff measurement was performed as the others did previously, and indicated a tendency for higher BP in recipient mice inoculated with stool samples from hypertensive donors as compared to controls. As such results were not obtained by fecal microbiota transplantation in conventionally raised mice, we speculate that the immune inflammatory system might play a crucial role in the pathogenesis of HTN. Further mechanism research to make clear whether gut bacterial metabolites show a contribution to the immune inflammatory system during the development of HTN is being performed.
In HTN studies, most work focused on patients with a clinical definition of HTN, who display a SBP higher than 140 mmHg or DBP ≥90 mmHg. However, population studies suggest that there is an intermediate stage of BP between control and HTN defined as pHTN, which should not be ignored. In our study, we considered subjects with pHTN as an independent group. Surprisingly, the bacterial diversity, enterotype, composition, and metabolic functions, as well as classified characteristics in pHTN highly coincided with those in HTN. As shown in Figs. 1, 2, and 3, there was a little difference in the structure of gut microbiome between pHTN and HTN, indicating that pHTN is not simply a transition stage between normotensive and hypertensive status upon BP levels but rather a state in which gut dysbiosis has already occurred. Moreover, our findings revealed indiscriminate metabolic profilings between pHTN and HTN, consistent with a previous report that the serum spectral profiles of the hosts were similar at a stage of SBP ≥130 mmHg and at SBP ≥150 mmHg . The close correlation of metabolic products and GM further strengthened and highlighted the importance of pHTN. Therefore, early treatment of pHTN has strong clinical value. In agreement with our notion, high BP has become one of the three leading risk factors for death according to the Global Burden of Disease Study . Moreover, our findings fully support the updated conclusion by the Systolic Blood Pressure Intervention Trial (SPRINT) research group, that controlling one’s SBP to an optimal level lower than 120 mmHg rather than a pHTN level below 140 mmHg leads to significantly decreased occurrence of cardiovascular events and death . Thus, more attention should be given to the previously neglected populations in pHTN, and early intervention for pHTN is strongly appealed.
Taken together, we have described clearly the disordered profiles of GM and microbial products in human patients with pHTN and HTN, established the relationship between gut dysbiosis and HTN, and provided important evidence for the novel role of GM dysbiosis as a key factor for BP changes. Our findings point towards a new strategy aimed at preventing the development of HTN and reducing cardiovascular risks through restoring the homeostasis of GM, by improving diet and lifestyle or early intervening with drugs or probiotics.
Study cohort and patient characteristics
All the individuals in the present study were from a cohort study among employees of the Kailuan Group Corporation. The Kailuan study is a prospective cohort study focusing on the Kailuan community in Tangshan, a large modern city in northern China, where 11 hospitals are responsible for the health care of the community, all of which participated in conducting physical examinations. All the subjects in the current work were strictly enrolled and none of them was under antihypertensive treatment. The participants were classified based on the Internal Guidelines for HTN as described in Additional file 1: Table S1. It was composed of 41 healthy controls (SBP ≤125 mmHg, or DBP ≤80 mmHg), 56 pHTNs (125 mmHg < SBP ≤ 139 mmHg, or 80 mmHg < DBP ≤ 89 mmHg), 99 patients of HTN (140 mmHg ≤ SBP, or 90 mmHg ≤ DBP). BP was measured in a sitting position by nurses or physicians. Three readings were recorded at 5-min intervals with a random-zero mercury column sphygmomanometer, and the average was taken as the final measurement.
All clinical information was collected according to standard procedures. Patients suffering from cancer, heart failure, renal failure, smoking, stroke, and peripheral artery disease were excluded, and none of the patients was under antihypertensive treatment. Healthy volunteers had no history of diabetes mellitus or hypercholesterolemia. Individuals were also excluded if they had received antibiotics or probiotics within the last 8 weeks. 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.
Stool sample collection and DNA extraction
Stool samples freshly collected from each participant were immediately frozen at −20 °C and transported to the laboratory with ice pack. Bacterial DNA was extracted at Novogene Bioinformatics Technology Co., Ltd using TIANGEN kit according to the manufacturer’s recommendations.
Metagenomic sequencing and gene catalogue construction
All samples were paired-end sequenced on the Illumina platform (insert size 300 bp, read length 125 bp) at the Novogene Bioinformatics Technology Co., Ltd. After quality control, the reads aligned to the human genome (alignment with SOAP2 , Version 2.21, parameters: -s 135, -l 30, -v 7,-m 200,-x 400) were removed. The remaining high-quality reads were used for further analysis.
The assembly of reads was executed using SOAP denovo (Version 2.04, parameters: -d 1 -M 3 -R -u -F) . For each sample, we used a series of k-mer values (from 49 to 87) and chose the optimal one with the longest N50 value for the remaining scaffolds . We mapped the clean data against scaffolds using SOAP2 (Version 2.21, parameters: -m 200 -x 400 -s 119). Unused reads from each sample were assembled using the same parameters. Genes (minimum length of 100 nucleotides) were predicted on scaftigs (i.e., continuous sequences within scaffolds) longer than 500 bp using MetaGeneMark (prokaryotic GeneMark.hmm version 2.10). Then, a non-redundant gene catalogue was constructed with CD-HIT (version 4.5.8, parameters: -G 0 -aS 0.9 -g 1 -d 0 -c 0.95)  using a sequence identity cut-off of 0.95, with a minimum coverage cut-off of 0.9 for the shorter sequences.
To determine the abundance of genes, reads were realigned to the gene catalogue with SOAP2 using parameters: -m 200 -x 400 -s 119. Only genes with ≥2 mapped reads were deemed to be present in a sample . The abundance of genes was calculated by counting the number of reads and normalizing by gene length.
α diversity and rarefaction curve
To estimate the genera richness of the sample, we calculated the within-sample (α) diversity using Shannon index based on the genera profiles. A high α diversity indicates a high richness of genera within the sample.
Rarefaction analysis was performed to assess the gene richness in the controls, pHTN, and HTN. For a given number of samples, we performed random sampling 100 times in the cohort with replacement and estimated the total number of genes that could be identified from these samples by R (Version 2.15.3, vegan package).
Microbial community types (enterotypes)
The community types of each sample were analyzed by the PAM method using relative abundance of genera. The optimal number of clusters was estimated using the CH index, as previously described . Only genera with an average relative abundance ≥10−4 and existed in at least six samples were considered in the analysis. According to Spearman’s correlation between genera abundances, the genera in enterotype 1 were clustered, and the co-occurrence network of them was visualized by Cytoscape (Version 3.2.1).
Taxonomic annotation and abundance profiling
To assess the taxonomic assignment, genes were aligned to the integrated NR database using DIAMOND (Version 0.7.9.58, default parameter except that -k 50 -sensitive -e 0.00001) . As previously described , for each gene, the significant matches, which were defined by e-values ≤10 × e-value of the top hit, and these retained matches were used to distinguish taxonomic groups. The taxonomical level of each gene was determined by the lowest common ancestor-based algorithm and implemented in MEGAN . The abundance of a taxonomic group was calculated by summing the abundance of genes annotated to a feature.
Metagenomic analysis in the verification phase
All phenotype information of participants were listed in the supplementary tables of Qin J et. al . Subjects with diabetes were excluded. Three HTN patients with SBP ≥150 mmHg or DBP ≥100 mmHg were enrolled, and three normotensive controls with SBP ≤125 mmHg and DBP ≤80 mmHg were included for the analysis.
Co-abundance gene groups
To identify the marker genes associated with pHTN and HTN, the abundance of each gene across groups was compared according to Greenblum S et al. . As previously described , these marker genes were clustered into groups based on their abundance variation across groups. Clusters with more than 50 genes were defined as co-abundance gene groups (CAGs), and used for further analysis. CAG abundance profiles were calculated by the average gene depth signal and weighted by gene length.
Taxonomic assignment of the CAGs was performed according to the taxonomy of tracer genes, as previously described . Briefly, assignment to species requires 90% of the genes in a CAG to align with the species’ genome with 95% identity and 70% overlap of query. Assigning CAG to a genus requires 80% of its genes to align to the genome with 85% identity in both DNA and protein sequences.
Co-occurrence network of CAGs
The enriched CAGs were identified according to Greenblum S et al.. These CAGs were clustered according to Spearman’s correlation. The co-occurrence network was visualized using Cytoscape (Version 3.2.1). The enriched CAGs/genes were identified according to Greenblum S et al. . Briefly, for each CAG, an OR score was calculated according to the abundance in the set of compared samples. Then, for the comparative analysis between control and HTN samples, the HTN-associated CAGs were classified as HTN-enriched (OR >2) or HTN-depleted (OR <0.5). When calculating HTN-associated ORs, samples of pHTN were excluded from the analysis, and when calculating pHTN-associated ORs, samples labeled as HTN were excluded.
Association between CAGs and clinical indices
Based on the clinical indices and enriched CAG profiles, we calculated Spearman’s correlation in all samples. The P values were corrected for multiple testing with Holm method by R (Version 2.15.3, psychpackage). Only 162 samples were considered in the analysis because of the clinical data missing in 34 samples.
All genes in our catalogue were aligned to the KEGG database (Release 73.1, with animal and plant genes removed) and CAZy database (http://www.cazy.org/) using DIAMOND (Version 0.7.9.58, default parameter except that -k 50 --sensitive -e 0.00001). Each protein was assigned to the KEGG orthology and CAZy families by the highest scoring annotated hits containing at least one HSP scoring over 60 bits . The abundance of KEGG orthology/module was calculated by summing the abundance of genes annotated to the same feature.
Metabonomics analysis based on LC/MS
One-hundred twenty-four (n = 124) individuals from our study cohort were subjected to metabonomics analysis based on the LC/MS method. This cohort was composed of 30 healthy controls, 31 pHTNs, and 63 patients of HTN. The whole blood samples were collected and separated into serum by centrifugation. Each serum samples at 100 μL were mixed with 400 μL methanol, and the mixtures were centrifuged at 12,000 rpm for 15 min at 4 °C. For LC/MS analysis, 200 μL of the supernatant was harvested.
The serum metabolic profiles were performed on a Thermo Fisher Ultimate 3000 LC system. For chromatographic separation, C18 (2.1 mm × 100 mm × 1.9 μm) reversed-phase column (Thermo Scientific, USA) preheated at 40 °C was used. A prepared serum sample of 4 μL was injected and maintained at 4 °C for analysis. The gradient conditions for metabolite elution were at 5% B for 0–1 min, 5–40% B for 1–2 min of linear gradient, 40–80% B for 2–11 min of linear gradient, and 95% B for 11–15 min. The mobile phase for positive ion mode (ES+) and negative ion mode (ES−) was composed of water with 0.1% formic acid as solvent A, and acetonitrile with 0.1% formic acid as solvent B, and the flow rate was at 300 μL/min.
For mass spectrometric assay, Orbitrap Elite mass spectrometer (Thermo Scientific, USA) equipped with ESI source was used to analyze the metabolite ions. The spray voltage was set to 3.8 kV in ES+ and 3.2 kV in ES−, the flow rate of sheath gas, aux gas, and sweep gas was 45, 15, and 1 arb, respectively. The ion source temperature was 300 °C, and the capillary temperature was 350 °C. Masses ranging from 50 to 1000 ion mass (m/z) were acquired, and the resolving power was set to 60,000.
The raw ESI data of LC/MS was converted into m/z format and analyzed for non-linear retention time (RT) alignment, peak detection, and filtration. Maximal spectrum of continuous wavelet transform was used to correct baseline and detect peak positions. Impurity peaks and duplicate identifications were eliminated. Compounds significantly different between groups were obtained at a variable influence on projection (VIP) >1.5, and P value of t test statistics <0.05 based on the peak intensities. The m/z values of these compounds were used to identify the metabolites corresponding to the featured peak in the Metlin database.
From the metabolite profile and 44 top differential genera abundance profile, Spearman’s correlation was performed to eliminate multi-collinearity and only one factor will be randomly selected from high correlated clusters (|r| ≥ 0.9) for further analysis. A stepwise regression of linear models was used for modeling the relationship between metabolites and related genera, from the fitting value of individual metabolites, and Spearman’s correlation between metabolic and associated genera was calculated and scaled by coefficients of each respective linear model.
GF C57BL/6L mice were obtained from Shanghai Institutes for Biological Sciences (SLAC Inc., Shanghai, China) and housed under a 12-h light–dark cycle in the gnotobiotic facilities. All mice were fed with sterile food and water ad libitum, and bacterial contamination was monitored by periodic examination of stools. For microbiota transplantation, the fresh fecal samples were collected from donors (Additional file 11: Table S17), resuspended with sterile saline, and centrifuged for supernatant. Male GF mice aged 8–10 weeks were randomly distributed into two groups and orally inoculated (200 μL for each mouse) twice at 1-day interval with prepared fecal contents from control or patients. Recipient mice transferred with microbiota were kept in different Trexler-type flexible film isolators, fed with sterile food and water, and bacterial contamination was strictly controlled. The gut microbial profiles of recipient mice were analyzed by 16S sequencing after 7 days. We chose a time point of 10 weeks post-transplantation for BP measurement. An assessment of BP was performed within 60 min after exporting the mice out of their gnotobiotic facilities, and we could not ensure prevention from bacterial contamination after the measurement; the BP at other time points during 10 weeks was not further examined. The BP was measured by the tail-cuff method and the BP-98A system (Softron, Tokyo, Japan), which was noninvasive and did not require surgery, since using direct invasive methods such as radiotelemetry techniques will immediately expose the mice to a non-sterile condition, which might impact the results. To acclimatize the mice undergoing the measurement procedures and improve measurement reliability, a heat-sterilized dark cover was transported into the germ-free mice isolator, where it was sterilized by spraying with a chlorine dioxide-based disinfectant in the isolator port. Before BP measurement, we have trained the mice by placing them in the dark cover in their sterile flexible film isolators without exporting them out at the same time for 14 days. To minimize contamination, the measurement was performed with UV-sterilized instruments under a sterile hood within 60 min after exporting the mice out of their sterile environment. All animal care and experiments were performed in accordance with the guidelines of Institutional Animal Care and Use Committee of SLAC Inc.
16S ribosomal RNA sequencing
16S rRNA community profiles were characterized using Illumina HiSeq sequencing of the V4 region (insert size 300 bp, read length 250 bp). Sequences were de novo clustered at 97% sequence identity and chimeras were removed using UPARSE . For each representative sequence, the GreenGene Database was used to annotate taxonomic information .
The Shannon index at the genera level was calculated with QIIME (Version 1.7.0). PCA was analyzed using the FactoMineR package in R software (Version 2.15.3). PCoA was performed and displayed by ade4 package, cluster packages, fpc packages, and clusterSim package in R software (Version 2.15.3). PLS-DA was performed using SIMCA-P software to cluster the sample plots across groups.
Differential abundance of genes, genera, and KO modules was tested by Wilcoxon rank sum test, and P values were corrected for multiple testing with the Benjamin & Hochberg method. Only genera with an average relative abundance ≥10−4 and existed in at least six subjects were considered in the analyses. Correlations between enriched CAGs and clinical indices were tested with Spearman’s correlation and visualized by Cytoscape (Version 3.2.1).
Using the profiles of species, CAGs, and metabolites, the samples were randomly divided into training set and test set. A random forest classifier was trained on 80% of the data and tested on the remaining 20% of our data using the random forest package in R. In order to evaluate the performance of the predictive model and get more precise curves, we used a 10-fold cross-validation within the training set. The cross-validational error curves (average of 10 test sets each) from five trials of the 10-fold cross-validation were averaged. Variable importance by mean decrease in accuracy was calculated for the random forest models using the full set of features. The number of variables was 1000 at the lowest cross-validational error. Thus, the predictive model was constructed using the 1000 most important variables, which were further applied for ROC analysis. The performance of the smaller models were measured as AUC when applied to the test set, and the confidence intervals for ROC curves were calculated using the pROC R package.
This work was supported by National Basic Research Program of China (973 Program,2014CB542302), CAMS Innovation Fund for Medical Science (CIFMS, 2016-I2M-1-006), National Natural Science Foundation of China (81470541, 81630014, 91531306, 81500383), National Basic Research Program of China (973 Program,2015CB554200), Beijing Municipal Science and Technology Commission (Z151100002115050, Z151100004015176), Beijing Municipal Commission of Education (KZ201610025028), the International S&T Cooperation Program of China (2013DFB30310), and the AHA Scientist Development Grant (L.C., 12SDG11680011).
Availability of data and materials
The data set supporting the results of this article has been deposited in the EMBL European Nucleotide Archive (ENA) under BioProject accession code PRJEB13870 [http://www.ebi.ac.uk/ena/data/view/PRJEB13870].
JC, XCY, BLZ, JL, and FQZ conceived the study, directed the project, designed the experiments, interpreted the results, and wrote the manuscript; JT, WLZ, TG, SLW, LY, and XYL obtained the samples and clinical details; YDW assisted with microbiome sample processing; JRC, WBL, and QHC performed the computational and metagenomic microbiota analysis; JL and JRC performed the metabonomics analysis; JL, JRC, YDW, and BG supervised the fecal microbiota transplantation and analyzed the data of animal experiments; WR, AK, and CL revised the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
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. All animal care and experiments were performed in accordance with the guidelines of Institutional Animal Care and Use Committee of SLAC Inc.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.
- Karlsson FH, Fak F, Nookaew I, Tremaroli V, Fagerberg B, Petranovic D, et al. Symptomatic atherosclerosis is associated with an altered gut metagenome. Nat Commun. 2012;3:1245.View ArticlePubMed CentralGoogle Scholar
- Qin J, Li Y, Cai Z, Li S, Zhu J, Zhang F, et al. A metagenome-wide association study of gut microbiota in type 2 diabetes. Nature. 2012;490:55–60.View ArticleGoogle Scholar
- Scher JU, Sczesnak A, Longman RS, Segata N, Ubeda C, Bielski C, et al. Expansion of intestinal Prevotella copri correlates with enhanced susceptibility to arthritis. Elife. 2013;2:e01202.View ArticlePubMed CentralGoogle Scholar
- Qin N, Yang F, Li A, Prifti E, Chen Y, Shao L, et al. Alterations of the human gut microbiome in liver cirrhosis. Nature. 2014;513:59–64.View ArticleGoogle Scholar
- 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:766.View ArticlePubMed CentralGoogle Scholar
- Fei N, Zhao L. An opportunistic pathogen isolated from the gut of an obese human causes obesity in germfree mice. ISME J. 2013;7:880–4.View ArticleGoogle Scholar
- Chang CJ, Lin CS, Lu CC, Martel J, Ko YF, Ojcius DM, et al. Ganoderma lucidum reduces obesity in mice by modulating the composition of the gut microbiota. Nat Commun. 2015;6:7489.View ArticlePubMed CentralGoogle Scholar
- Zheng P, Zeng B, Zhou C, Liu M, Fang Z, Xu X, et al. Gut microbiome remodeling induces depressive-like behaviors through a pathway mediated by the host's metabolism. Mol Psychiatry. 2016;21:786–96.View ArticleGoogle Scholar
- Schaubeck M, Clavel T, Calasan J, Lagkouvardos I, Haange SB, Jehmlich N, et al. Dysbiotic gut microbiota causes transmissible Crohn's disease-like ileitis independent of failure in antimicrobial defence. Gut. 2016;65:225–37.View ArticleGoogle Scholar
- Le Roy T, Llopis M, Lepage P, Bruneau A, Rabot S, Bevilacqua C, et al. Intestinal microbiota determines development of non-alcoholic fatty liver disease in mice. Gut. 2013;62:1787–94.View ArticleGoogle Scholar
- Llopis M, Cassard AM, Wrzosek L, Boschat L, Bruneau A, Ferrere G, et al. Intestinal microbiota contributes to individual susceptibility to alcoholic liver disease. Gut. 2016;65:830–9.View ArticleGoogle Scholar
- Koeth RA, Wang Z, Levison BS, Buffa JA, Org E, Sheehy BT, et al. Intestinal microbiota metabolism of L-carnitine, a nutrient in red meat, promotes atherosclerosis. Nat Med. 2013;19:576–85.View ArticlePubMed CentralGoogle Scholar
- Wang Z, Roberts AB, Buffa JA, Levison BS, Zhu W, Org E, et al. Non-lethal inhibition of gut microbial trimethylamine production for the treatment of atherosclerosis. Cell. 2015;163:1585–95.View ArticlePubMed CentralGoogle Scholar
- Hillege HL, Girbes AR, de Kam PJ, Boomsma F, de Zeeuw D, Charlesworth A, et al. Renal function, neurohormonal activation, and survival in patients with chronic heart failure. Circulation. 2000;102:203–10.View ArticleGoogle Scholar
- Pierdomenico SD, Di Nicola M, Esposito AL, Di Mascio R, Ballone E, Lapenna D, et al. Prognostic value of different indices of blood pressure variability in hypertensive patients. Am J Hypertens. 2009;22:842–7.View ArticleGoogle Scholar
- He FJ, Li J, Macgregor GA. Effect of longer-term modest salt reduction on blood pressure. Cochrane Database Syst Rev. 2013, DOI: 10.1002/14651858.CD004937.pub2: CD004937.
- 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:531–8.View ArticlePubMed CentralGoogle Scholar
- Moghadamrad S, McCoy KD, Geuking MB, Sagesser H, Kirundi J, Macpherson AJ, et al. Attenuated portal hypertension in germ-free mice: function of bacterial flora on the development of mesenteric lymphatic and blood vessels. Hepatology. 2015;61:1685–95.View ArticleGoogle Scholar
- Yang T, Santisteban MM, Rodriguez V, Li E, Ahmari N, Carvajal JM, et al. Gut dysbiosis is linked to hypertension. Hypertension. 2015;65:1331–40.View ArticlePubMed CentralGoogle Scholar
- 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:187–97.View ArticlePubMed CentralGoogle Scholar
- Khalesi S, Sun J, Buys N, Jayasinghe R. Effect of probiotics on blood pressure: a systematic review and meta-analysis of randomized, controlled trials. Hypertension. 2014;64:897–903.View ArticleGoogle Scholar
- Gomez-Guzman M, Toral M, Romero M, Jimenez R, Galindo P, Sanchez M, et al. Antihypertensive effects of probiotics Lactobacillus strains in spontaneously hypertensive rats. Mol Nutr Food Res. 2015;59:2326–36.View ArticleGoogle Scholar
- Durgan DJ, Ganesh BP, Cope JL, Ajami NJ, Phillips SC, Petrosino JF, et al. Role of the gut microbiome in obstructive sleep apnea-induced hypertension. Hypertension. 2016;67:469–74.Google Scholar
- Weber MA, Schiffrin EL, White WB, Mann S, Lindholm LH, Kenerson JG, et al. Clinical practice guidelines for the management of hypertension in the community: a statement by the American Society of Hypertension and the International Society of Hypertension. J Clin Hypertens (Greenwich). 2014;16:14–26.View ArticleGoogle Scholar
- Forslund K, Hildebrand F, Nielsen T, Falony G, Le Chatelier E, Sunagawa S, et al. Disentangling type 2 diabetes and metformin treatment signatures in the human gut microbiota. Nature. 2015;528:262–6.View ArticlePubMed CentralGoogle Scholar
- Tanaka S, Yoshida M, Murakami Y, Ogiwara T, Shoji M, Kobayashi S, et al. The relationship of Prevotella intermedia, Prevotella nigrescens and Prevotella melaninogenica in the supragingival plaque of children, caries and oral malodor. J Clin Pediatr Dent. 2008;32:195–200.View ArticleGoogle Scholar
- Podschun R, Ullmann U. Klebsiella spp. as nosocomial pathogens: epidemiology, taxonomy, typing methods, and pathogenicity factors. Clin Microbiol Rev. 1998;11:589–603.PubMedPubMed CentralGoogle Scholar
- Restuccia PA, Cunha BA. Klebsiella. Infect Control. 1984;5:343–7.View ArticlePubMedGoogle Scholar
- Wang X, Yue TL, Barone FC, White RF, Clark RK, Willette RN, et al. Discovery of adrenomedullin in rat ischemic cortex and evidence for its role in exacerbating focal brain ischemic damage. Proc Natl Acad Sci U S A. 1995;92:11480–4.View ArticlePubMedPubMed CentralGoogle Scholar
- Machiels K, Joossens M, Sabino J, De Preter V, Arijs I, Eeckhaut V, et al. A decrease of the butyrate-producing species Roseburia hominis and Fecalibacterium prausnitzii defines dysbiosis in patients with ulcerative colitis. Gut. 2014;63:1275–83.View ArticleGoogle Scholar
- Sokol H, Pigneur B, Watterlot L, Lakhdari O, Bermudez-Humaran LG, Gratadoux JJ, et al. Fecalibacterium prausnitzii is an anti-inflammatory commensal bacterium identified by gut microbiota analysis of Crohn disease patients. Proc Natl Acad Sci U S A. 2008;105:16731–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Bryant MP, Small N. The anaerobic monotrichous butyric acid-producing curved rod-shaped bacteria of the rumen. J Bacteriol. 1956;72:16–21.PubMedPubMed CentralGoogle Scholar
- Pinzone MR, Celesia BM, Di Rosa M, Cacopardo B, Nunnari G. Microbial translocation in chronic liver diseases. Int J Microbiol. 2012;2012:694629.View ArticlePubMedPubMed CentralGoogle Scholar
- Greenblum S, Turnbaugh PJ, Borenstein E. Metagenomic systems biology of the human gut microbiome reveals topological shifts associated with obesity and inflammatory bowel disease. Proc Natl Acad Sci U S A. 2012;109:594–9.View ArticlePubMedGoogle Scholar
- Nielsen HB, Almeida M, Juncker AS, Rasmussen S, Li J, Sunagawa S, et al. Identification and assembly of genomes and genetic elements in complex metagenomic samples without using reference genomes. Nat Biotechnol. 2014;32:822–8.View ArticlePubMedGoogle Scholar
- Aliyu SH, Marriott RK, Curran MD, Parmar S, Bentley N, Brown NM, et al. Real-time PCR investigation into the importance of Fusobacterium necrophorum as a cause of acute pharyngitis in general practice. J Med Microbiol. 2004;53:1029–35.View ArticlePubMedGoogle Scholar
- Gevers D, Kugathasan S, Denson LA, Vazquez-Baeza Y, Van Treuren W, Ren B, et al. The treatment-naive microbiome in new-onset Crohn’s disease. Cell Host Microbe. 2014;15:382–92.View ArticlePubMedPubMed CentralGoogle Scholar
- 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:D480–484.View ArticlePubMedGoogle Scholar
- Feng Q, Liang S, Jia H, Stadlmayr A, Tang L, Lan Z, et al. Gut microbiome development along the colorectal adenoma-carcinoma sequence. Nat Commun. 2015;6:6528.View ArticlePubMedGoogle Scholar
- Zhang X, Zhang D, Jia H, Feng Q, Wang D, Liang D, et al. The oral and gut microbiomes are perturbed in rheumatoid arthritis and partly normalized after treatment. Nat Med. 2015;21:895–905.View ArticlePubMedGoogle Scholar
- Tang XY, Zhu YQ, Wei B, Wang H. Expression and functional research of TLR4 in human colon carcinoma. Am J Med Sci. 2010;339:319–26.View ArticlePubMedGoogle Scholar
- Cario E, Rosenberg IM, Brandwein SL, Beck PL, Reinecker HC, Podolsky DK. Lipopolysaccharide activates distinct signaling pathways in intestinal epithelial cell lines expressing Toll-like receptors. J Immunol. 2000;164:966–72.View ArticlePubMedGoogle Scholar
- Kumar S, Arya P, Mukherjee C, Singh BK, Singh N, Parmar VS, et al. Novel aromatic ester from Piper longum and its analogues inhibit expression of cell adhesion molecules on endothelial cells. Biochemistry. 2005;44:15944–52.View ArticlePubMedGoogle Scholar
- Wang W, Zheng JP, Zhu SX, Guan WJ, Chen M, Zhong NS. Carbocisteine attenuates hydrogen peroxide-induced inflammatory injury in A549 cells via NF-kappaB and ERK1/2 MAPK pathways. Int Immunopharmacol. 2015;24:306–13.View ArticlePubMedGoogle Scholar
- Lu Y, AJ, Wang G, Hao H, Huang Q, Yan B, et al. Gas chromatography/time-of-flight mass spectrometry based metabonomic approach to differentiating hypertension- and age-related metabolic variation in spontaneously hypertensive rats. Rapid Commun Mass Spectrom. 2008;22:2882–8.View ArticlePubMedGoogle Scholar
- De Deyn PP, Robitaille P, Vanasse M, Qureshi IA, Marescau B. Serum guanidino compound levels in uremic pediatric patients treated with hemodialysis or continuous cycle peritoneal dialysis. Correlations between nerve conduction velocities and altered guanidino compound concentrations. Nephron. 1995;69:411–7.View ArticlePubMedGoogle Scholar
- Qi Y, Aranda JM, Rodriguez V, Raizada MK, Pepine CJ. Impact of antibiotics on arterial blood pressure in a patient with resistant hypertension—a case report. Int J Cardiol. 2015;201:157–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Baez S, Gordon HA. Tone and reactivity of vascular smooth muscle in germfree rat mesentery. J Exp Med. 1971;134:846–56.View ArticlePubMedPubMed CentralGoogle Scholar
- Hofer U. Microbiome: anelloviridae go viral. Nat Rev Microbiol. 2014;12:4–5.View ArticlePubMedGoogle Scholar
- Kohn FR, Kung AH. Role of endotoxin in acute inflammation induced by gram-negative bacteria: specific inhibition of lipopolysaccharide-mediated responses with an amino-terminal fragment of bactericidal/permeability-increasing protein. Infect Immun. 1995;63:333–9.PubMedPubMed CentralGoogle Scholar
- Cybulsky MI, Chan MK, Movat HZ. Acute inflammation and microthrombosis induced by endotoxin, interleukin-1, and tumor necrosis factor and their implication in gram-negative infection. Lab Invest. 1988;58:365–78.PubMedGoogle Scholar
- Karbach SH, Schönfelder T, Brandão I, Wilms E, Hörmann N, Jäckel S, et al. Gut microbiota promote angiotensin II-induced arterial hypertension and vascular dysfunction. J Am Heart Assoc. 2016;5(9). doi:10.1161/JAHA.116.003698.
- Brindle JT, Nicholson JK, Schofield PM, Grainger DJ, Holmes E. Application of chemometrics to 1H NMR spectroscopic data to investigate a relationship between human serum metabolic profiles and hypertension. Analyst. 2003;128:32–6.View ArticlePubMedGoogle Scholar
- Lim SS, Vos T, Flaxman AD, Danaei G, Shibuya K, Adair-Rohani H, et al. A comparative risk assessment of burden of disease and injury attributable to 67 risk factors and risk factor clusters in 21 regions, 1990–2010: a systematic analysis for the Global Burden of Disease Study 2010. Lancet. 2012;380:2224–60.View ArticlePubMedPubMed CentralGoogle Scholar
- Group SR, Wright Jr JT, Williamson JD, Whelton PK, Snyder JK, Sink KM, et al. A randomized trial of intensive versus standard blood-pressure control. N Engl J Med. 2015;373:2103–16.View ArticleGoogle Scholar
- Li R, Yu C, Li Y, Lam TW, Yiu SM, Kristiansen K, et al. SOAP2: an improved ultrafast tool for short read alignment. Bioinformatics. 2009;25:1966–7.View ArticlePubMedGoogle Scholar
- Luo R, Liu B, Xie Y, Li Z, Huang W, Yuan J, et al. SOAPdenovo2: an empirically improved memory-efficient short-read de novo assembler. Gigascience. 2012;1:18.View ArticlePubMedPubMed CentralGoogle Scholar
- Li W, Godzik A. Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics. 2006;22:1658–9.View ArticlePubMedGoogle Scholar
- Qin J, Li R, Raes J, Arumugam M, Burgdorf KS, Manichanh C, et al. A human gut microbial gene catalogue established by metagenomic sequencing. Nature. 2010;464:59–65.View ArticlePubMedPubMed CentralGoogle Scholar
- Arumugam M, Raes J, Pelletier E, Le Paslier D, Yamada T, Mende DR, et al. Enterotypes of the human gut microbiome. Nature. 2011;473:174–80.View ArticlePubMedPubMed CentralGoogle Scholar
- Buchfink B, Xie C, Huson DH. Fast and sensitive protein alignment using DIAMOND. Nat Methods. 2015;12:59–60.View ArticlePubMedGoogle Scholar
- Huson DH, Auch AF, Qi J, Schuster SC. MEGAN analysis of metagenomic data. Genome Res. 2007;17:377–86.View ArticlePubMedPubMed CentralGoogle Scholar
- Backhed F, Roswall J, Peng Y, Feng Q, Jia H, Kovatcheva-Datchary P, et al. Dynamics and stabilization of the human gut microbiome during the first year of life. Cell Host Microbe. 2015;17:690–703.View ArticlePubMedGoogle Scholar
- Edgar RC. UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat Methods. 2013;10:996–8.View ArticlePubMedGoogle Scholar
- DeSantis TZ, Hugenholtz P, Larsen N, Rojas M, Brodie EL, Keller K, et al. Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl Environ Microbiol. 2006;72:5069–72.View ArticlePubMedPubMed CentralGoogle Scholar