Alterations of the Gut Microbiome in Hypertension
Qiulong Yan1,2 
Yifang Gu3 Xiangchun Li4 Wei Yang5 Liqiu Jia1 
Changming Chen1 Xiuyan Han1 Yukun Huang1 Lizhe Zhao1 
Peng Li3 Zhiwei Fang3 Junpeng Zhou3 Xiuru Guan5 Yanchun Ding6 Shaopeng Wang7 Muhammad Khan8 
Yi Xin9 
Shenghui Li3* 
Yufang Ma1*- 1Department of Biochemistry and Molecular Biology, Dalian Medical University, Dalian, China
- 2Department of Microbiology, Dalian Medical University, Dalian, China
- 3Shenzhen Puensum Genetech Institute, Shenzhen, China
- 4Beijing Genomics Institute, Shenzhen, China
- 5Department of Laboratory Diagnostics, The First Affiliated Hospital of Harbin Medical University, Harbin, China
- 6Department of Cardiology V, The Second Affiliated Hospital of Dalian Medical University, Dalian, China
- 7Department of Cardiology, The First Affiliated Hospital of Dalian Medical University, Dalian, China
- 8College of Basic Medical Sciences, Dalian Medical University, Dalian, China
- 9Department of Biotechnology, Dalian Medical University, Dalian, China
Introduction: Human gut microbiota is believed to be directly or indirectly involved in cardiovascular diseases and hypertension. However, the identification and functional status of the hypertension-related gut microbe(s) have not yet been surveyed in a comprehensive manner.
Methods: Here we characterized the gut microbiome in hypertension status by comparing fecal samples of 60 patients with primary hypertension and 60 gender-, age-, and body weight-matched healthy controls based on whole-metagenome shotgun sequencing.
Results: Hypertension implicated a remarkable gut dysbiosis with significant reduction in within-sample diversity and shift in microbial composition. Metagenome-wide association study (MGWAS) revealed 53,953 microbial genes that differ in distribution between the patients and healthy controls (false discovery rate, 0.05) and can be grouped into 68 clusters representing bacterial species. Opportunistic pathogenic taxa, such as, Klebsiella spp., Streptococcus spp., and Parabacteroides merdae were frequently distributed in hypertensive gut microbiome, whereas the short-chain fatty acid producer, such as, Roseburia spp. and Faecalibacterium prausnitzii, were higher in controls. The number of hypertension-associated species also showed stronger correlation to the severity of disease. Functionally, the hypertensive gut microbiome exhibited higher membrane transport, lipopolysaccharide biosynthesis and steroid degradation, while in controls the metabolism of amino acid, cofactors and vitamins was found to be higher. We further provided the microbial markers for disease discrimination and achieved an area under the receiver operator characteristic curve (AUC) of 0.78, demonstrating the potential of gut microbiota in prediction of hypertension.
Conclusion: These findings represent specific alterations in microbial diversity, genes, species and functions of the hypertensive gut microbiome. Further studies on the causality relationship between hypertension and gut microbiota will offer new prospects for treating and preventing the hypertension and its associated diseases.
Introduction
Hypertension is a global public health problem. In 2010, about 31% of the world's population has been estimated to suffer from hypertension and over 1 billon of this population is living in low- and middle- income countries(Mittal and Singh, 2010; Mills et al., 2016). Hypertension is one of the major risk factors for cardiovascular diseases, such as, stroke and heart failure (Lim et al., 2012; Faraco and Iadecola, 2013). Moreover, it is believed to be one of the most common comorbidities associated with chronic renal disease (Lash et al., 2009), obesity and type 2 diabetes (Landsberg and Molitch, 2004; Kotchen, 2010). Presently, genome-wide association studies (GWAS) have identified a series of genetic loci and pathways associated with blood pressure (Xu et al., 2015; Liu et al., 2016). The environmental factors, such as, dietary salt intake, alcohol consumption and lack of exercise, are also linked to the occurrence of hypertension (Fuchs et al., 2001; Karppanen and Mervaala, 2006). Recent practice of metabolomics also identified new pathogenic pathways involved in blood pressure regulation (Menni et al., 2015; Galla et al., 2017). Nevertheless, due to the complexity and heterogeneity of hypertension, identification of the causes of this disease is still challenging.
Recent studies have demonstrated that the gut microflora plays an essential role in development of cardiovascular diseases, via metabolizing dietary choline, phosphatidylcholine and L-carnitine to produce trimethylamine (TMA), which is further oxidized into TMA N-oxide (TMAO, a metabolite that enhances atherosclerosis; Wang et al., 2011; Koeth et al., 2013; Tang et al., 2013). Even though the direct link between hypertension and TMAO has not been established currently, TMAO's role to prolong the hypertensive effect of angiotensin II were reported (Ufnal et al., 2014). Inhibition of gut microbiota-mediated TMAO production may serve as a potential therapeutic approach for the treatment of cardiometabolic diseases (Wang Z. et al., 2015). These findings suggest an intricate and predictable correlation between hypertension and gut microbiota. To validate this, a recent study based on metagenomic analyses of the fecal samples of 41 healthy controls, 56 pre-hypertension subjects, and 99 hypertension individuals described a novel causal role of aberrant gut microbiota in contributing to the pathogenesis of hypertension, and emphasized the significance of early intervention for pre-hypertension (Li et al., 2017). Moreover, rat experiments have linked gut microbial dysbiosis with hypertension (Mell et al., 2015; Yang et al., 2015; Adnan et al., 2017; Santisteban et al., 2017). The causal role of gut microbiome in obstructive sleep apnea-induced hypertension have been reported (Durgan et al., 2016). Here, to investigate the alteration of the human gut microbiome underlying hypertension, we compared the microbial composition of fecal samples obtained from 60 patients with primary hypertension and 60 healthy counterparts of Chinese origin. We used quantitative metagenomic analysis to identify genic, microbial, and functional characteristics underlying hypertension.
Methods
Subjects and Sample Collection
Sixty primary hypertensive patients (current blood pressure ≥140/90 mm Hg) and sixty gender-, age-, and body weight-matched healthy controls (current blood pressure ≤ 120/80 mm Hg) were recruited for this study. Other than systolic blood pressure (SBP) and diastolic blood pressure (DBP), the other clinical parameters have no significant differences in the two groups of populations, except for triglyceride (TG). The characteristics of the subjects are summarized in Table 1, and detailed information is given in Table S1. Subjects were excluded if they had symptoms of respiratory infection or digestive tract disease, or if they were treated with antibiotics or anti-inflammatory agents in recent 2 months before sampling. Subjects with hypertension or severe cardiovascular diseases (such as, coronary artery disease or stroke) history in previous 5 years were also excluded from healthy controls. Fresh fecal samples were collected from each subject and were stored at a −80°C freezer immediately.
Table 1. Characteristics of subjects.
Ethics Statement
This study received approval from the ethics committee of The First Affiliated Hospital of Harbin Medical University, and written informed consent was obtained from each participant. The methods were carried out in accordance with the approved guidelines.
DNA Preparation and Sequencing
Genomic DNA was extracted from all samples according to a modified protocol provided in the QIAamp DNA mini kit (Qiagen, Manchester, UK; Yan et al., 2016). Briefly, ASL buffer (1.4 ml) was added to 220 mg of fecal sample and the pellets were homogenized in a 2 ml screw cap tubes (Axygen) by vortex. The suspension was incubated at 95°C for 5 min to lyse bacterial cells. After centrifugation (13,000 × g, 1 min) and incubation with an InhibitEx tablet, the supernatant was treated with 15 μl proteinase K and 200 μl Buffer AL at 70°C for 10 min. The extracted DNA was dissolved in 100 μl sterile water. Paired-end DNA libraries (insert size 350 bp, read length 150 bp) were constructed according to the manufacturer's instructions (Illumina, USA). Whole-metagenome shotgun sequencing was performed on the Illumina HiSeq3000 platform. Further methodological detail is available in the Supplementary Methods.
Bioinformatic Analysis
Quantification of Metagenomic Genes and Species
High-quality reads from each sample were aligned to the integrated non-redundant human gut gene catalog (IGC; Li et al., 2014) using SOAP2 (Li et al., 2009; >90% similarity). The relative abundance of a gene in a sample was estimated by dividing the number of reads uniquely mapped to that gene by the length of gene region and by the total number of reads from the sample. We also aligned the sequencing reads against the available microbial genomes (bacteria, archaea, and virus) from the National Center for Biotechnology Information (NCBI) database and generated the taxonomic compositions (i.e., phylum and genus composition) for all samples.
Alpha Diversity
The gene count (Le Chatelier et al., 2013) of a metagenomic sample were calculated based on their mapped reads number on the gene catalog (to eliminate the influence of sequencing amount fluctuation, 10 million reads were randomly extracted from each sample for mapping). The Shannon index (within-sample diversity) was calculated based on the gene relative abundance profiles, using the method described previously (Qin et al., 2012).
Functional Annotation and Profiling
The Kyoto Encyclopedia of Genes and Genomes (KEGG, downloaded Jan-2016) database (Kanehisa et al., 2014) was used for functional annotation of genes. Amino acid sequences were searched against the databases using USEARCH v8.1 (Edgar, 2010) with a minimum similarity of 30%. Each gene was assigned a KEGG ortholog (KO) based on the best hit protein. The abundance profiles of KO were calculated by summing the relative abundance of its genes. The choline-trimethylamine lyase (cutC, KO: K20038; Craciun and Balskus, 2012) was used to evaluate the gut microbiota-mediated TMA production in subjects, and the short-chain fatty acid (SCFA)-producing enzymes were represented by acetyl-CoA decarbonylase/synthase (K00193, K00194, K00197, K14138, which are key enzymes of the acetate biosynthesis pathways: KEGG modules M00377 and M00422; Koh et al., 2016), propionyl-CoA:succinate-CoA transferase (Reichardt et al., 2014), butyryl-CoA:acetate CoA-transferase (K01034, K01035), and butyrate kinase (K00929) (Pryde et al., 2002; Louis et al., 2010).
Metagenome-Wide Association Study
We used the metagenome-wide association study (MGWAS) method to identify gene markers that showed significant abundance differences between hypertensive patients and control subjects. The MGWAS was performed using the methodology developed by Qin et al. (2012). Co-abundance genes were clustered into metagenomic linkage groups (MLGs) based on the previous methods (Qin et al., 2012). Taxonomic assignment and abundance profiling of the MLGs were performed according to the taxonomy and the relative abundance of their constituent genes (see Supplementary Methods for detail). MLGs were considered to be interacted if absolute value of Spearman's correlation coefficient between them is greater than 0.4, and the co-occurrence network of MLGs was visualized by Cytoscape (Shannon et al., 2003).
Statistical Analyses
Statistical analyses were implemented using the R platform. Distance-based redundancy analysis (dbRDA) was performed on normalized taxa abundance matrices with R vegan package (Dixon, 2003) according to Bray-Curtis distance, then visualized with R ggplot2 package. Random forest models were trained with R randomForest package (10,000 trees) to predict hypertension status according to MLG abundance profiles. The performance of the predictive model was evaluated with cross-validation error. Receiver operator characteristic (ROC) analysis was performed using R pROC package. P-value < 0.05 was considered statistical significance, and the q-value was calculated to evaluate the false discovery rate (FDR) for correction of multiple comparisons.
Data Availability
The raw whole-metagenomic shotgun sequencing data acquired in this study have been deposited to the European Bioinformatics Institute (EBI) database under the accession code ERP023883.
Results
Comparison of the Gut Microbiota between Hypertensive Patients and Controls
To investigate the gut microbial composition of 60 hypertensive patients and 60 healthy controls, we obtained 652.9 Gbp high-quality data (5.4 ± 1.1 Gbp per sample) via whole-metagenome shotgun sequencing on their fecal samples. When we quantified the microbial (alpha) diversity within each subject, the patients showed significantly lower gene count and Shannon index compared with the controls (Figure 1A). Multivariate analysis based on Bray-Curtis distance between microbial genera revealed remarkable differences between patients and controls (Figure 1B). At the phylum level, patients had higher levels of Proteobacteria (p < 0.01), but fewer Actinobacteria (p = 0.02). At the genus level, Klebsiella, Clostridium, Streptococcus, Parabacteroides, Eggerthella, and Salmonella were frequently distributed in hypertensive gut compared to normal controls while Faecalibacterium, Roseburia, and Synergistetes were found to be higher in control group compared to hypertensive patients (Figure 1C). These findings demonstrated considerable gut microbial dysbiosis in hypertensive patients.
Figure 1. Difference of gut microbial community between hypertensive patients and controls. (A), Difference of alpha diversity between hypertensive (HT) patients and controls. (B), dbRDA based on the Bray-Curtis distances between microbial genera, revealing a hypertensive microbial dysbiosis which overlaps only in part with taxonomic composition in patients and controls. The first two principle components and the ratio of variance contributed is shown. Genera (square) as the main contributors are plotted by their loadings in these two components. Lines connect samples in the same group, and circles cover samples near the center of gravity for each group. (C), Boxplot shows the significantly different genera between patients and controls. Genera with q < 0.05 (Mann-Whitney U-test corrected by FDR) are shown. Only the genera with average relative abundances greater than 0.05% of total abundance in all samples are shown for clarity. Gray and white boxes represent the patients and controls, respectively. For A and C, the boxes represent the interquartile range (IQR) between first and third quartiles and the line inside represents the median. The whiskers denote the lowest and highest values within 1.5 times IQR from the first and third quartiles, respectively. The dots represent outliers beyond the whiskers.
Identification of Hypertension-Associated Markers from Gut Microbiome
To explore signatures of the gut microbiome in hypertensive patients and controls, we integrated the sequencing data into an existing gut microbial reference gene catalog to obtain a set of 5.3 million genes, which allowed for saturation mapping of the reads (80.3%). Using the MGWAS methods, we identified 53,953 genes that showed a significant difference between two groups (FDR corrected q < 0.05). Approximately, 69% of these genes were clustered into 68 metagenome linkage groups (MLGs, Table S2), that allowed to species level description for the microbiome differences. Thirty-one MLGs were higher in patients while 37 in controls. Consistent with the genus level observations, MLGs of Klebsiella (mainly consisting of K. pneumoniae and K. variicola), Streptococcus (S. infantarius, S. pasteurianus and S. salivarius), and Parabacteroides merdae were found to be higher in hypertensive samples, whereas MLGs of Roseburia (mainly consisting of R. intestinalis and R. hominis) and Faecalibacterium prausnitzii were higher in controls. Moreover, the MLGs enriched in hypertensive patients also contain several Bacteroides spp. (including B. eggerthii and B. cellulosilyticus), Sutterella wadsworthensis and Pyramidobacter piscolens, and the MLGs enriched in controls include several other Bacteroides spp. (including B. uniformis, B. nordii and B. dorei), Megasphaera spp. (M. micronuciformis), and Aeromicrobium massiliense. A co-occurrence network on these MLGs revealed a large number of interconnections within hypertension-enriched and control-enriched MLGs (Figure 2A), as well as some MLGs derived from two groups negatively correlated. This result suggested that the MLGs did not occur independently and interacted with the taxa in its environment.
Figure 2. Characterization and interconnection of hypertension-associated MLGs. (A) Co-occurrence network shows the interconnection of the hypertension- and control-enriched MLGs. Nodes depict MLGs with their ID or taxonomic assignment (unclassified MLGs under genus or higher taxonomy rank are marked by “*”) displayed in the center. The size of the nodes indicates the number of gene within the MLG. Connecting lines represent Spearman correlation coefficient ρ > 0.40 (represented by blue line) or < −0.40 (represented by red line). (B), Correlation of gross abundance of hypertension- and control-enriched MLGs with hypertension stage. NS, not significant; *, q < 0.05; **, q < 0.01; Mann-Whitney U-test corrected by FDR.
We next found that the gross abundances of hypertension- and control-enriched MLGs are correlated to the severity of hypertension (Figure 2B), suggesting that the bacterial relative abundance of these MLGs could be related to the development and disease progress of hypertension.
Functional Characterization of Gut Microbiota
Based on the KEGG pathway comparison, we revealed that the hypertensive gut microbiomes were more abundant in membrane transport, lipopolysaccharide (LPS) biosynthesis, and steroid degradation (Figure 3A and Table S3), while the controls were enriched in metabolism of “other amino acids,” cofactors and vitamins (including folate biosynthesis and metabolism, riboflavin metabolism, and ubiquinone biosynthesis). In addition, the gut microbial enzymes involved in TMA production were enriched in the hypertensive patients compared to controls, whereas the SCFA-producing enzymes were depleted (Figure 3B).
Figure 3. Functional comparison of the gut microbiomes between hypertensive patients and healthy controls. (A), Distributions of relative abundances of KEGG pathway categories in hypertensive patients and controls. *, q < 0.05; **, q < 0.01; Mann-Whitney U-test corrected by FDR. (B), Difference of the relative abundance of cutC (TMA-producing) and SCFA-producing enzymes between hypertensive (HT) patients and controls.
Gut Microbiota-Based Classification of Hypertension
We evaluated the performance of gut microbiota composition to identify hypertension status in the MLG profiles using the Random Forest model, and obtained the discriminatory power of the area under the ROC curve (AUC) of 0.78 (95% CI 0.73–0.82; Figure 4A). Several control-enriched MLGs (including Clostridiales, Blautia hansenii, Megasphaera) and two hypertension-enriched members of Streptococcus (S. salivarius and S. infantarius) featured the highest score for the discrimination of hypertensive patients and healthy controls (Figure 4B).
Figure 4. Classification of hypertension status by the abundances of MLGs. (A), ROC analysis for classification of hypertensive status by MLGs, assessed by AUC. (B), The 20 most discriminant MLGs in the model classifying hypertensive patients and healthy controls. The bar lengths indicate the importance of the variable, and colors represent enrichment in patients (black) or controls (white).
Discussion
To identify and analyze the differences of the gut microbiota in hypertension, we characterized the genic, microbial, and functional repertoire of the microbiomes of 60 hypertensive patients and 60 gender-, age-, and body weight-matched controls. Our study strengthened previous metagenomic study on gut microbiome of hypertension (Li et al., 2017) by adding more information. Furthermore, we observed significant differences in microbial community dysbiosis, taxonomic shifts, and functional changes between hypertensive- and control-gut microbiome.
Previous studies showed that the gut microbes participate in choline and phosphatidylcholine metabolism to form circulating and urinary TMAO, while high levels of plasma TMAO promote accelerated atherosclerosis and increase the risk of cardiovascular disorders (Tang et al., 2013; Wang Z. et al., 2015). The choline utilization (cutC) gene, a critical gene that coverts the choline to trimethylamine, was identified in a variety of human gut commensals belonging to Firmicutes, Proteobacteria, and Actinobacteria (Craciun and Balskus, 2012). Notably, several genera such as, Klebsiella, Clostridium, and Streptococcus, which are highly distributed in hypertensive patients are choline degraders (Hakenbeck et al., 2009; Craciun and Balskus, 2012; Kalnins et al., 2015). Functional analysis also showed that the abundance of cutC gene was enriched in the gut microbiota of the hypertensive patients. These findings suggested that the dietary choline intake and TMAO production via gut microflora would be a probable pathway for hypertensive pathogenesis.
Klebsiella, is a pathogen routinely found in human gut that causes pneumonia, diarrhea, and urinary tract infection. The distribution of Klebsiella was found to be significantly higher in hypertensive patients compared to healthy controls as evident from Figure 1C. Overgrowth of Klebsiella usually foreshadows gut flora dysbiosis, which leads to a variety of serious chronic disease, such as, colitis (Garrett et al., 2010), Crohn's disease and ankylosing spondylitis (Ebringer et al., 2007). The present study demonstrates that Klebsiella species which are highly distributed in hypertensive patients are K. pneumoniae (the main component of Klebsiella that associated with nosocomial infection and multiple diseases), K. variicola [a human and animal opportunistic pathogen that is associated with bovine mastitis (Brisse and Duijkeren, 2005)], and four unclassified MLGs. Based on these information, however, the potential correlation between Klebsiella and hypertension is still unclear.
Streptococcus, the dominant species of human oral microbiome (Wade, 2013) that causes upper respiratory tract infection, were also found highly distributed in gut microbiota of hypertensive patients as compared to the controls. Gut streptococci is also associated with diseases, such as, inflammatory bowel disease (Conte et al., 2006) and liver cirrhosis (Qin et al., 2014). It has been reported previously that oral cavity and/or gut might be the source of streptococci found in the majority of atherosclerotic plaque microbiota (Koren et al., 2011). These findings suggest that possible correlation of gut streptococci in hypertension.
F. prausnitzii and Roseburia spp., which were abundantly distributed in controls compared to hypertensive patients, were also distributed abundantly in the healthy control microbiomes of many chronic diseases, including type 2 diabetes (Qin et al., 2012), liver cirrhosis (Qin et al., 2014), Crohn's disease (Gevers et al., 2014), and ulcerative colitis (Machiels et al., 2014). F. prausnitzii and Roseburia (both R. intestinalis and R. hominis) are the major SCFA producer in human colon (Shoaie et al., 2015), which might explain the depletion of SCFA-producing enzymes in hypertensive gut microbiome. Functionally, SCFAs modulates the gut inflammation and metabolism via functioning as important colonocytes energy source and signaling molecules (Donohoe et al., 2011), suggesting that low level of SCFA production in gut microbiota may be a considerable risk factor of multiple metabolic syndromes and hypertension.
Several other bacteria also played important function in human gut and showed potential function in hypertension, such as, the patient-enriched Bacteroides (including B. eggerthii, B. cellulosilyticus, and 3 unclassified Bacteroides MLGs) and Parabacteroides (P. merdae) which are generally opportunistic pathogens in infectious diseases and are able to develop antimicrobial drug resistance (Boente et al., 2010), and the control-enriched Megasphaera spp. (M. micronuciformis and two unclassified MLGs) which are producer of SCFAs, vitamins and essential amino acids (Shetty et al., 2013). In addition, co-abundance analysis (Figure 2A) generated a striking number of positive correlations within the patient/control-enriched MLGs and negative correlations between the two groups, revealing that a comprehensive bacterial synergism and antagonism existed in the human gut. In this case, the microbial dysbiosis of hypertensive gut microbiome would not be determined by independent pathogens (e.g., the patient-enriched MLGs), but more likely to be caused by a series of risk factors (e.g., improper diet or lifestyle that inhibit the growth of beneficial bacterium) that change the balance of ecosystem. Intriguingly, the severity of hypertension was positively correlated with the total abundance of patient-enriched MLGs and negatively correlated with those of control-enriched MLGs (Figure 2B), suggesting that the bacterial relative abundance may also be a potential risk factor of hypertension development. Such a “dose response” was also found in the gut microbiome of liver cirrhosis (Qin et al., 2014) and colorectal adenoma-carcinoma patients (Feng et al., 2015).
Our study further provided the microbial markers for hypertension discrimination, and achieved an AUC of 0.78 for identifying disease status based on 68 species-level MLGs. This discriminatory power was higher than that from the prediction models based on genomic markers identified by GWAS (Evans et al., 2009; Fava et al., 2013), and was almost at same level with the phenotype-based models (AUC 0.71-0.81) (Echouffo-Tcheugui et al., 2013; Wang A. et al., 2015). Thus, the fecal microbiota showed a good potential on prediction and early diagnosis of hypertension, however, systematic investigations of key species and gene markers identified here might be helpful in the future.
Drug-induced gut microbiome shifts were observed during the treatment of multiple diseases, such as, the metformin therapy in type 2 diabetes (Forslund et al., 2015) and antirheumatic drugs therapy in rheumatoid arthritis (Zhang et al., 2015). In this study, a part of patients (~35%) had taken antihypertensive drugs or specific nutritious supplementary, however, the relationship between drug treatment and gut microbiota is still unclear. Another significant limitation of this study is that the gut microbial community would be sensitive to environmental factors, such as, host race, geography, life and diet style, and so on. Although our samples were age-, gender-, BMI-matched, some phenotype differences were still unobservable. To avoid this, larger cohort containing multi-types of hypertensive patients are needed for further investigation. Generally, hypertension is a highly complex and heterogeneous disease, it is still infeasible at this moment to draw any conclusions about causal relationships of gut microbiota and hypertension, and direct experimental studies (e.g., the animal model studies) are needed to show causality of proposed microbes or pathways.
In summary, our finding extends previous knowledge of correlation between gut microbiota and hypertension in animal models (Yang et al., 2015; Durgan et al., 2016) and provides a range of signatures in metagenomic diversity, genes, species, and functions of the hypertensive gut microbiome. Further studies on the causality relationship between hypertension and gut microbiota will lead to a better understanding of the mutual interaction.
Author Contributions
YM, SL, and QY designed experiments; QY, LJ, CC, XH, YH, LZ, PL, and MK carried out experiments; QY, WY, YD, SW, and YX analyzed experimental results. SL, YG, XL, ZF, and JZ analyzed sequencing data. WY, XG, YD, and SW collected the samples. YM, SL, and QY wrote the manuscript.
Funding
This study was supported by grants from the National Naturel Science Foundation of China (81573469) and the National Basic Research Program of China (2012CB518803).
Conflict of Interest Statement
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.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fcimb.2017.00381/full#supplementary-material
References
Adnan, S., Nelson, J. W., Ajami, N. J., Venna, V. R., Petrosino, J. F., Bryan, R. M. Jr., et al. (2017). Alterations in the gut microbiota can elicit hypertension in rats. Physiol. Genomics 49, 96–104. doi: 10.1152/physiolgenomics.00081.2016
Boente, R. F., Ferreira, L. Q., Falcao, L. S., Miranda, K. R., Guimaraes, P. L., Domingues, R. M., et al. (2010). Detection of resistance genes and susceptibility patterns in Bacteroides and Parabacteroides strains. Anaerobe 16, 190–194. doi: 10.1016/j.anaerobe.2010.02.003
Brisse, S., and Duijkeren, E. (2005). Identification and antimicrobial susceptibility of 100 Klebsiella animal clinical isolates. Vet. Microbiol. 105, 307–312. doi: 10.1016/j.vetmic.2004.11.010
Conte, M. P., Schippa, S., Zamboni, I., Penta, M., Chiarini, F., Cucchiara, S., et al. (2006). Gut-associated bacterial microbiota in paediatric patients with inflammatory bowel disease. Gut 55, 1760–1767. doi: 10.1136/gut.2005.078824
Craciun, S., and Balskus, E. P. (2012). Microbial conversion of choline to trimethylamine requires a glycyl radical enzyme. Proc. Natl. Acad. Sci. U.S.A. 109, 21307–21312. doi: 10.1073/pnas.1215689109
Dixon, P. (2003). VEGAN, a package of R functions for community ecology. J. Veg. Sci. 14, 927–930. doi: 10.1111/j.1654-1103.2003.tb02228.x
Donohoe, D. R., Garge, N., Zhang, X., Sun, W., O'Connell, T. M., Bultman, S. J., et al. (2011). The microbiome and butyrate regulate energy metabolism and autophagy in the mammalian colon. Cell Metab. 13, 517–526. doi: 10.1016/j.cmet.2011.02.018
Durgan, D. J., Ganesh, B. P., Cope, J. L., Ajami, N. J., Phillips, S. C., Petrosino, J. F., et al. (2016). Role of the gut microbiome in obstructive sleep apnea-induced hypertension. Hypertension 67, 469–474. doi: 10.1161/HYPERTENSIONAHA.115.06672
Ebringer, A., Rashid, T., Tiwana, H., and Wilson, C. (2007). A possible link between Crohn's disease and ankylosing spondylitis via Klebsiella infections. Clin. Rheumatol. 26, 289–297. doi: 10.1007/s10067-006-0391-2
Echouffo-Tcheugui, J. B., Batty, G. D., Kivimaki, M., and Kengne, A. P. (2013). Risk models to predict hypertension: a systematic review. PLoS ONE 8:e67370. doi: 10.1371/journal.pone.0067370
Edgar, R. C. (2010). Search and clustering orders of magnitude faster than BLAST. Bioinformatics 26, 2460–2461. doi: 10.1093/bioinformatics/btq461
Evans, D. M., Visscher, P. M., and Wray, N. R. (2009). Harnessing the information contained within genome-wide association studies to improve individual prediction of complex disease risk. Hum. Mol. Genet. 18, 3525–3531. doi: 10.1093/hmg/ddp295
Faraco, G., and Iadecola, C. (2013). Hypertension: a harbinger of stroke and dementia. Hypertension 62, 810–887. doi: 10.1161/HYPERTENSIONAHA.113.01063
Fava, C., Sjogren, M., Montagnana, M., Danese, E., Almgren, P., Engstrom, G., et al. (2013). Prediction of blood pressure changes over time and incidence of hypertension by a genetic risk score in Swedes. Hypertension 61, 319–326. doi: 10.1161/HYPERTENSIONAHA.112.202655
Feng, Q., Liang, S., Jia, H., Stadlmayr, A., Tang, L., Lan, Z., et al. (2015). Gut microbiome development along the colorectal adenoma-carcinoma sequence. Nat. Commun. 6:6528. doi: 10.1038/ncomms7528
Forslund, K., Hildebrand, F., Nielsen, T., Falony, G., Le Chatelier, E., Sunagawa, S., et al. (2015). Disentangling type 2 diabetes and metformin treatment signatures in the human gut microbiota. Nature 528, 262–266. doi: 10.1038/nature15766
Fuchs, F. D., Chambless, L. E., Whelton, P. K., Nieto, F. J., and Heiss, G. (2001). Alcohol consumption and the incidence of hypertension: the atherosclerosis risk in communities study. Hypertension 37, 1242–1250. doi: 10.1161/01.HYP.37.5.1242
Galla, S., Chakraborty, S., Mell, B., Vijay-Kumar, M., and Joe, B. (2017). Microbiotal-host interactions and hypertension. Physiology 32, 224–233. doi: 10.1152/physiol.00003.2017
Garrett, W. S., Gallini, C. A., Yatsunenko, T., Michaud, M., DuBois, A., Glimcher, L. H., et al. (2010). Enterobacteriaceae act in concert with the gut microbiota to induce spontaneous and maternally transmitted colitis. Cell Host Microbe 8, 292–300. doi: 10.1016/j.chom.2010.08.004
Gevers, D., Kugathasan, S., Denson, L. A., Vazquez-Baeza, Y., Van Treuren, W., Xavier, R. J., et al. (2014). The treatment-naive microbiome in new-onset Crohn's disease. Cell Host Microbe 15, 382–392. doi: 10.1016/j.chom.2014.02.005
Hakenbeck, R., Madhour, A., Denapaite, D., and Bruckner, R. (2009). Versatility of choline metabolism and choline-binding proteins in Streptococcus pneumoniae and commensal streptococci. FEMS Microbiol. Rev. 33, 572–586. doi: 10.1111/j.1574-6976.2009.00172.x
Kalnins, G., Kuka, J., Grinberga, S., Makrecka-Kuka, M., Liepinsh, E., Dambrova, M., et al. (2015). Structure and function of CutC Choline Lyase from human microbiota bacterium Klebsiella pneumoniae. J. Biol. Chem. 290, 21732–21740. doi: 10.1074/jbc.M115.670471
Kanehisa, M., Goto, S., Sato, Y., Kawashima, M., Furumichi, M., and Tanabe, M. (2014). Data, information, knowledge and principle: back to metabolism in KEGG. Nucleic Acids Res. 42, D199–D205. doi: 10.1093/nar/gkt1076
Karppanen, H., and Mervaala, E. (2006). Sodium intake and hypertension. Prog. Cardiovasc. Dis. 49, 59–75. doi: 10.1016/j.pcad.2006.07.001
Koeth, R. A., Wang, Z., Levison, B. S., Buffa, J. A., Org, E., Hazen, S. L., et al. (2013). Intestinal microbiota metabolism of L-carnitine, a nutrient in red meat, promotes atherosclerosis. Nat. Med. 19, 576–585. doi: 10.1038/nm.3145
Koh, A., De Vadder, F., Kovatcheva-Datchary, P., and Backhed, F. (2016). From dietary fiber to host physiology: short-chain fatty acids as key bacterial metabolites. Cell 165, 1332–1345. doi: 10.1016/j.cell.2016.05.041
Koren, O., Spor, A., Felin, J., Fak, F., Stombaugh, J., Tremaroli, V., et al. (2011). Human oral, gut, and plaque microbiota in patients with atherosclerosis. Proc. Natl. Acad. Sci. U.S.A. 108(Suppl. 1), 4592–4598. doi: 10.1073/pnas.1011383107
Kotchen, T. A. (2010). Obesity-related hypertension: epidemiology, pathophysiology, and clinical management. Am. J. Hypertens. 23, 1170–1178. doi: 10.1038/ajh.2010.172
Landsberg, L., and Molitch, M. (2004). Diabetes and hypertension: pathogenesis, prevention and treatment. Clin. Exp. Hypertens. 26, 621–628. doi: 10.1081/CEH-200031945
Lash, J. P., Go, A. S., Appel, L. J., He, J., Ojo, A., Chronic Renal Insufficiency Cohort Study, G., et al. (2009). Chronic Renal Insufficiency Cohort (CRIC) study: baseline characteristics and associations with kidney function. Clin. J. Am. Soc. Nephrol. 4, 1302–1311. doi: 10.2215/CJN.00070109
Le Chatelier, E., Nielsen, T., Qin, J., Prifti, E., Hildebrand, F., Falony, G., et al. (2013). Richness of human gut microbiome correlates with metabolic markers. Nature 500, 541–546. doi: 10.1038/nature12506
Li, J., Jia, H., Cai, X., Zhong, H., Feng, Q., Sunagawa, S., et al. (2014). An integrated catalog of reference genes in the human gut microbiome. Nat. Biotechnol. 32, 834–841. doi: 10.1038/nbt.2942
Li, J., Zhao, F., Wang, Y., Chen, J., Tao, J., Tian, G., et al. (2017). Gut microbiota dysbiosis contributes to the development of hypertension. Microbiome 5:14. doi: 10.1186/s40168-016-0222-x
Li, R., Yu, C., Li, Y., Lam, T. W., Yiu, S. M., Wang, J., et al. (2009). SOAP2: an improved ultrafast tool for short read alignment. Bioinformatics 25, 1966–1967. doi: 10.1093/bioinformatics/btp336
Lim, S. S., Vos, T., Flaxman, A. D., Danaei, G., Shibuya, K., Memish, Z. A., et al. (2012). 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 380, 2224–2260. doi: 10.1016/S0140-6736(12)61766-8
Liu, C., Kraja, A. T., Smith, J. A., Brody, J. A., Franceschini, N., Chasman, D. I., et al. (2016). Meta-analysis identifies common and rare variants influencing blood pressure and overlapping with metabolic trait loci. Nat. Genet. 48, 1162–1170. doi: 10.1038/ng.3660
Louis, P., Young, P., Holtrop, G., and Flint, H. J. (2010). Diversity of human colonic butyrate-producing bacteria revealed by analysis of the butyryl-CoA:acetate CoA-transferase gene. Environ. Microbiol. 12, 304–314. doi: 10.1111/j.1462-2920.2009.02066.x
Machiels, K., Joossens, M., Sabino, J., De Preter, V., Arijs, I., Eeckhaut, V., et al. (2014). A decrease of the butyrate-producing species Roseburia hominis and Faecalibacterium prausnitzii defines dysbiosis in patients with ulcerative colitis. Gut 63, 1275–1283. doi: 10.1136/gutjnl-2013-304833
Mell, B., Jala, V. R., Mathew, A. V., Byun, J., Waghulde, H., Joe, B., et al. (2015). Evidence for a link between gut microbiota and hypertension in the Dahl rat. Physiol. Genomics 47, 187–197. doi: 10.1152/physiolgenomics.00136.2014
Menni, C., Graham, D., Kastenmuller, G., Alharbi, N. H., Alsanosi, S. M., Valdes, A. M., et al. (2015). Metabolomic identification of a novel pathway of blood pressure regulation involving hexadecanedioate. Hypertension 66, 422–429. doi: 10.1161/HYPERTENSIONAHA.115.05544
Mills, K. T., Bundy, J. D., Kelly, T. N., Reed, J. E., Kearney, P. M., He, J., et al. (2016). Global disparities of hypertension prevalence and control: a systematic analysis of population-based studies from 90 countries. Circulation 134, 441–450. doi: 10.1161/CIRCULATIONAHA.115.018912
Mittal, B. V., and Singh, A. K. (2010). Hypertension in the developing world: challenges and opportunities. Am. J. Kidney Dis. 55, 590–598. doi: 10.1053/j.ajkd.2009.06.044
Pryde, S. E., Duncan, S. H., Hold, G. L., Stewart, C. S., and Flint, H. J. (2002). The microbiology of butyrate formation in the human colon. FEMS Microbiol. Lett. 217, 133–139. doi: 10.1111/j.1574-6968.2002.tb11467.x
Qin, J., Li, Y., Cai, Z., Li, S., Zhu, J., Zhang, F., Wang, J., et al. (2012). A metagenome-wide association study of gut microbiota in type 2 diabetes. Nature 490, 55–60. doi: 10.1038/nature11450
Qin, N., Yang, F., Li, A., Prifti, E., Chen, Y., Shao, L., et al. (2014). Alterations of the human gut microbiome in liver cirrhosis. Nature 513, 59–64. doi: 10.1038/nature13568
Reichardt, N., Duncan, S. H., Young, P., Belenguer, A. C., McWilliam Leitch, C., Scott, K., et al. (2014). Phylogenetic distribution of three pathways for propionate production within the human gut microbiota. ISME J. 8, 1323–1335. doi: 10.1038/ismej.2014.14
Santisteban, M. M., Qi, Y., Zubcevic, J., Kim, S., Yang, T., Raizada, M. K., et al. (2017). Hypertension-linked pathophysiological alterations in the gut. Circ. Res. 120, 312–323. doi: 10.1161/CIRCRESAHA.116.309006
Shannon, P., Markiel, A., Ozier, O., Baliga, N. S., Wang, J. T., Ideker, T., et al. (2003). Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 13, 2498–2504. doi: 10.1101/gr.1239303
Shetty, S. A., Marathe, N. P., Lanjekar, V., Ranade, D., and Shouche, Y. S. (2013). Comparative genome analysis of Megasphaera sp. reveals niche specialization and its potential role in the human gut. PLoS ONE 8:e79353. doi: 10.1371/journal.pone.0079353
Shoaie, S., Ghaffari, P., Kovatcheva-Datchary, P., Mardinoglu, A., Sen, P., Pujos-Guillot, E., et al. (2015). Quantifying diet-induced metabolic changes of the human gut microbiome. Cell Metab. 22, 320–331. doi: 10.1016/j.cmet.2015.07.001
Tang, W. H., Wang, Z., Levison, B. S., Koeth, R. A., Britt, E. B., Hazen, S. L., et al. (2013). Intestinal microbial metabolism of phosphatidylcholine and cardiovascular risk. N. Engl. J. Med. 368, 1575–1584. doi: 10.1056/NEJMoa1109400
Ufnal, M., Jazwiec, R., Dadlez, M., Drapala, A., Sikora, M., and Skrzypecki, J. (2014). Trimethylamine-N-oxide: a carnitine-derived metabolite that prolongs the hypertensive effect of angiotensin II in rats. Can. J. Cardiol. 30, 1700–1705. doi: 10.1016/j.cjca.2014.09.010
Wade, W. G. (2013). The oral microbiome in health and disease. Pharmacol. Res. 69, 137–143. doi: 10.1016/j.phrs.2012.11.006
Wang, A., An, N., Chen, G., Li, L., and Alterovitz, G. (2015). Predicting hypertension without measurement: a non-invasive, questionnaire-based approach. Expert Syst. Appl. 42, 7601–7609. doi: 10.1016/j.eswa.2015.06.012
Wang, Z., Klipfell, E., Bennett, B. J., Koeth, R., Levison, B. S., Hazen, S. L., et al. (2011). Gut flora metabolism of phosphatidylcholine promotes cardiovascular disease. Nature 472, 57–63. doi: 10.1038/nature09922
Wang, Z., Roberts, A. B., Buffa, J. A., Levison, B. S., Zhu, W., Hazen, S. L., et al. (2015). Non-lethal inhibition of gut microbial trimethylamine production for the treatment of atherosclerosis. Cell 163, 1585–1595. doi: 10.1016/j.cell.2015.11.055
Xu, K., Ma, L., Li, Y., Wang, F., Zheng, G. Y., Tian, X. L., et al. (2015). Genetic and functional evidence supports LPAR1 as a susceptibility gene for hypertension. Hypertension 66, 641–646. doi: 10.1161/HYPERTENSIONAHA.115.05515
Yan, Q., Cui, S., Chen, C., Li, S., Sha, S., Wan, X., et al. (2016). Metagenomic analysis of sputum microbiome as a tool toward culture-independent pathogen detection of patients with ventilator-associated pneumonia. Am. J. Respir. Crit. Care Med. 194, 636–639. doi: 10.1164/rccm.201601-0034LE
Yang, T., Santisteban, M. M., Rodriguez, V., Li, E., Ahmari, N., Mohamadzadeh, M., et al. (2015). Gut dysbiosis is linked to hypertension. Hypertension 65, 1331–1340. doi: 10.1161/HYPERTENSIONAHA.115.05315
Keywords: hypertension, gut microbiome, microbial dysbiosis, metagenome-wide association study
Citation: Yan Q, Gu Y, Li X, Yang W, Jia L, Chen C, Han X, Huang Y, Zhao L, Li P, Fang Z, Zhou J, Guan X, Ding Y, Wang S, Khan M, Xin Y, Li S and Ma Y (2017) Alterations of the Gut Microbiome in Hypertension. Front. Cell. Infect. Microbiol. 7:381. doi: 10.3389/fcimb.2017.00381
Received: 16 May 2017; Accepted: 09 August 2017;
Published: 24 August 2017.
Edited by:
Venkatakrishna Rao Jala, University of Louisville, United StatesReviewed by:
Bina Joe, University of Toledo, United StatesMorgan Langille, Dalhousie University, Canada
Copyright © 2017 Yan, Gu, Li, Yang, Jia, Chen, Han, Huang, Zhao, Li, Fang, Zhou, Guan, Ding, Wang, Khan, Xin, Li and Ma. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Shenghui Li, lishenghui@puensum.com
Yufang Ma, yufang_ma@hotmail.com