Next Article in Journal
Predicting Risk of Antenatal Depression and Anxiety Using Multi-Layer Perceptrons and Support Vector Machines
Next Article in Special Issue
‘Statistical Irreproducibility’ Does Not Improve with Larger Sample Size: How to Quantify and Address Disease Data Multimodality in Human and Animal Research
Previous Article in Journal
GBP5 Serves as a Potential Marker to Predict a Favorable Response in Triple-Negative Breast Cancer Patients Receiving a Taxane-Based Chemotherapy
Previous Article in Special Issue
Intestinal Dysbiosis in Young Cystic Fibrosis Rabbits
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
JPM
Volume 11
Issue 3
10.3390/jpm11030198
jpm-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Differences in the Microbial Composition of Hemodialysis Patients Treated with and without β-Blockers

by
Yi-Ting Lin
1,2,3,
Ting-Yun Lin
4,5,
Szu-Chun Hung
4,5,
Po-Yu Liu
6,
Wei-Chun Hung
7,
Wei-Chung Tsai
8,
Yi-Chun Tsai
2,9,10,
Rachel Ann Delicano
11,
Yun-Shiuan Chuang
1,
Mei-Chuan Kuo
2,10,12,
Yi-Wen Chiu
2,10,12 and
Ping-Hsun Wu
2,3,12,*
1
Department of Family Medicine, Kaohsiung Medical University Hospital, Kaohsiung 80708, Taiwan
2
Faculty of Medicine, College of Medicine, Kaohsiung Medical University, Kaohsiung 80708, Taiwan
3
Graduate Institute of Clinical Medicine, College of Medicine, Kaohsiung Medical University, Kaohsiung 80708, Taiwan
4
Division of Nephrology, Taipei Tzu Chi Hospital, Buddhist Tzu Chi Medical Foundation, New Taipei City 23142, Taiwan
5
School of Medicine, Tzu Chi University, Hualien 97071, Taiwan
6
Department of Internal Medicine, College of Medicine, National Taiwan University, Taipei 100225, Taiwan
7
Department of Microbiology and Immunology, Kaohsiung Medical University, Kaohsiung 80708, Taiwan
8
Division of Cardiology, Department of Internal Medicine, Kaohsiung Medical University Hospital, Kaohsiung Medical University, Kaohsiung 80708, Taiwan
9
Division of General Medicine, Kaohsiung Medical University Hospital, Kaohsiung Medical University, Kaohsiung 80708, Taiwan
10
Faculty of Renal Care, College of Medicine, Kaohsiung Medical University, Kaohsiung 80708, Taiwan
11
Institute of Surgical Sciences, Uppsala University, 752 36 Uppsala, Sweden
12
Division of Nephrology, Department of Internal Medicine, Kaohsiung Medical University Hospital, Kaohsiung Medical University, Kaohsiung 80708, Taiwan
 
add Show full affiliation list
 
remove Hide full affiliation list
*
Author to whom correspondence should be addressed.
J. Pers. Med. 2021, 11(3), 198; https://0-doi-org.brum.beds.ac.uk/10.3390/jpm11030198
Submission received: 6 February 2021 / Revised: 21 February 2021 / Accepted: 6 March 2021 / Published: 12 March 2021
(This article belongs to the Special Issue Targeting the Microbiome for Disease Diagnosis and Therapy: New Frontiers for Personalized Medicine)
Browse Figures
Versions Notes

Abstract

:
β-blockers are commonly prescribed to treat cardiovascular disease in hemodialysis patients. Beyond the pharmacological effects, β-blockers have potential impacts on gut microbiota, but no study has investigated the effect in hemodialysis patients. Hence, we aim to investigate the gut microbiota composition difference between β-blocker users and nonusers in hemodialysis patients. Fecal samples collected from hemodialysis patients (83 β-blocker users and 110 nonusers) were determined by 16S ribosomal RNA amplification sequencing. Propensity score (PS) matching was performed to control confounders. The microbial composition differences were analyzed by the linear discriminant analysis effect size, random forest, and zero-inflated Gaussian fit model. The α-diversity (Simpson index) was greater in β-blocker users with a distinct β-diversity (Bray–Curtis Index) compared to nonusers in both full and PS-matched cohorts. There was a significant enrichment in the genus Flavonifractor in β-blocker users compared to nonusers in full and PS-matched cohorts. A similar finding was demonstrated in random forest analysis. In conclusion, hemodialysis patients using β-blockers had a different gut microbiota composition compared to nonusers. In particular, the Flavonifractor genus was increased with β-blocker treatment. Our findings highlight the impact of β-blockers on the gut microbiota in hemodialysis patients.

Graphical Abstract

1. Introduction

The gut microbiota has a crucial role in metabolic, nutritional, physiological, defensive, and immunological processes in the human body, with its composition linked to human health and the development of diseases [1,2]. Human-microbiome association can be considered as integration in evolution. The microbiome can modulate and restore human health [3]. Changes in this microbial equilibrium, that is, dysbiosis, promotes and influences the course of many intestinal and extra-intestinal diseases [4,5,6]. In addition to genetic and environmental factors, several common medications (e.g., proton pump inhibitors, nonsteroidal anti-inflammatory drugs, atypical antipsychotics, selective serotonin reuptake inhibitors, antibiotics, statins, and antidiabetic drugs) are associated with the specific gut microbiota composition [7,8,9,10,11,12,13]. Indeed, drug-microbiome-host interactions are complex and multifactorial, impacting host metabolism [14,15]. Hence, they should be part of the core phenotype set for human gut microbiota research [16].
Patients with chronic kidney disease (CKD) have altered gut microbiota, with a bidirectional causal effect relationship [17,18]. Among the cardiovascular preventive drugs for patients with end-stage renal disease (ESRD), β-blockers are commonly prescribed in higher cardiovascular risk patients to prevent sudden cardiac death [19,20]. Beyond the clinical effect of β-blockers in ESRD patients, they also have a potential impact on gut microbiota [7,16]. Besides, the benefit of beta-blockers may be attributed to preventing the activity of the gut microbe-generated metabolite, such as phenylacetylglutamine [21]. However, limited study has investigated the impact on ESRD patients. Herein, we evaluate the gut microbiota composition of β-blocker users and nonusers in Taiwanese hemodialysis patients.

2. Materials and Methods

2.1. Study Participants

The Ethics Committee approved the study protocols of Kaohsiung Medical University Hospital (KMUHIRB-E(I)-20160095 and KMUHIRB-E(I)-20180118) and Taipei Tzu Chi Hospital (07-X01-002). All participants provided written informed consent. Hemodialysis patients were recruited from the dialysis unit of Taipei Tzu Chi Hospital and Kaohsiung Medical University Hospital in Taiwan from August 2017 to February 2018. The inclusion criteria were patients with age more than 18 years old and received regular hemodialysis three times per week, 3.5–4 h with high-flux dialyzers. The exclusion criteria included patients with partial or total colectomy, inflammatory bowel diseases, active malignancies, or patients who were prescribed antibiotics within three months before enrollment. Fecal samples were collected from 193 stable hemodialysis patients and analyzed by high-throughput 16S ribosomal RNA gene sequencing to compared participants with and without β-blocker treatment. All β-blocker users were prescribed for at least one month.

2.2. Comorbidity, Laboratory, and Clinical Variables

All baseline characteristics of sociodemographic data, age, sex, body mass index, dialysis vintage, arteriovenous shunt type, comorbidities, medications, and biochemical data were collected in the built-in electronic health care system. Blood samples were collected after overnight fasting through the arteriovenous fistula or graft before scheduled hemodialysis sessions. The biochemical data included serum values for hemoglobin, albumin, high sensitivity C reactive protein, total cholesterol, low-density lipoprotein, triglycerides, ion calcium, and phosphate from routine blood samples obtained within 30 days before stool sample collection. Diet was evaluated by a licensed dietitian using a modified short-form food frequency questionnaire. No specific antioxidant supplements (i.e., tea, cocoa products, or wine) were recorded because of strict dietary restrictions in hemodialysis patients. Participants have followed the nutrition guideline of the National Kidney Foundation’s Kidney Disease Outcomes Quality Initiative (KDOQI™) [22], which recommends a high-protein intake (1.1–1.4 g/kg/day) and reduced consumption of fruits, vegetables, and dietary fiber to avoid potassium overload. Diabetes was defined as HbA1C 6.5% or higher or use of oral antidiabetic agents or insulin. Hypertension was defined as 140/90 mmHg or higher or taking blood pressure-lowering drugs. A history of myocardial infarction or documented by coronary angiography, class III or IV congestive heart failure, or a cerebrovascular accident were defined as cardiovascular disease.

2.3. Fecal Sample Collection and Bacterial 16S rRNA Amplicon Sequencing and Processing

All stool samples were frozen immediately after collection by each participant, then delivered in cooler bags to the laboratory (Germark Biotechnology, Taichung, Taiwan) within 24 h. A QIAamp DNA Stool Mini Kit (Qiagen, MD, USA) was used to extract DNA from fecal samples. Barcode-indexed PCR primers (341F and 805R) were used to create an amplicon library by amplifying the variable regions 3 and 4 (V3–V4) of the 16S rRNA gene [23]. The amplicons were sequenced (300 bp paired-end) using an Illumina MiSeq sequencer at the same time in the same laboratory to avoid batch effects (Germark Biotechnology, Taichung, Taiwan). The 16S-amplicon pipeline was adapted from 16S Bacteria/Archaea SOP v1 of Microbiome Helper workflows [24]. Paired-End reAd mergeR (PEAR; version 0.9.8) [25] was used to merge paired-end reads to raw reads, then filtered low-quality reads by thresholds of sequence length ≥400 bp and quality score of 90% bases of reads ≥20. Quantitative Insight Into Microbial Ecology (QIIME; version 1.9.1) software was used to select operational taxonomic units (OTU) [26]. The SILVA (version 123) 16S database [27,28] was applied to cluster OTUs and assign taxonomy using the UCLUST algorithm (version v1.2.22q) [29] with a 97% sequence identity threshold. Reads were dereplicated, and singletons were discarded. The final OTU table was rarefied into minimum sequencing depth in the data set.

2.4. Propensity Score Matching

Propensity score (PS) matching [30,31] was performed to balance confounders between the comparisons of interest (i.e., β-blocker users versus nonusers) and minimize the confounding by indication resulting from nonrandom treatment study. Using a logistic regression model, β-blocker use was accessed to estimate the propensity to receive a β-blocker for each participant based on potential confounders, including age, sex, body mass index, dialysis vintage, smoking history, vascular access type, Bristol stool scale, dietary intake, comorbidities (diabetes mellitus, hypertension, dyslipidemia, coronary artery disease, heart failure, cerebrovascular disease, and parathyroidectomy history), concomitant drugs used (including ACEI (angiotensin converting enzyme inhibitors)/ARB (angiotensin-receptor blockers), glucose-lowering drugs (such as sulfonylurea, dipeptidyl peptidase-4 inhibitors, insulin), statin, calcium carbonate, and proton pump inhibitors), and clinical laboratory data (hemoglobin, albumin, total cholesterol, triglyceride, high sensitivity C reactive protein (hsCRP), sodium, potassium, total calcium, phosphate, parathyroid hormone, serum iron, ferritin, normalized protein catabolic rate (nPCR), and single pool Kt/V). In this study, 193 hemodialysis patients were enrolled, including 83 β-blocker users and 110 nonusers (full cohort). PS-matched (1:1) analysis was used to match participants with β-blocker treatment (N = 62) to participants without β-blocker treatment (N = 62) (PS-matched cohort, Figure 1).

2.5. Statistical and Bioinformatics Analyses of Microbiota

The study design is presented in Figure 1. Demographic characteristics are shown as the mean, median, or frequency, with differences between β-blocker users and nonusers determined using an independent T-test or chi-squared test, as appropriate. A rarefaction curve was built to prevent methodological artifacts originating from variations in sequencing depth. The α-diversity indices (Shannon and Simpson’s indices) estimated the evenness of taxa within each sample and were generated using the R “vegan” package and calculated the p-value by the Kruskal–Wallis test. The β-diversity provides a comparison of the taxonomic profiles’ differences between pairs of individual samples. The β-diversity was calculated based on the Bray–Curtis distance matrices and was visualized using principal coordinates analysis (PCoA) and calculated using homogeneity of group dispersions by Permutational Analysis of Multivariate Dispersions (PERMDISP) [32].
Co-occurrence analysis was used to determine the relationships within communities, with core microbiome analysis performed at the genus level using MicrobiomeAnalyst [33], in which sample prevalence and relative abundance cut-off values were set at 20 and 0.2%, respectively. For visualization of the internal interactions and further measurement of the microbial community, Sparse Correlations for Compositional data (SparCC) was used to calculate the Spearman correlation coefficient with the corresponding p-value between every two taxa. Microbiota community structure was assessed by co-occurrence networks built by the SparCC algorithm [34]. The p-values were estimated by 100 random permutations and iterations for each SparCC calculation, and correlation matrices were computed from the resampled data matrices. Only OTUs with correlation scores greater than 0.4 and p-value less than 0.05 were categorized into co-abundance groups (CAGs); these coefficients were also used to assess the length of edges on the network. An undirected network, weighted by SparCC correlation magnitude, was generated using bioinformatics tools in MicrobiomeAnalyst [33].
The bacterial OTU difference between β-blocker users and nonusers was analyzed by the linear discriminant analysis (LDA) of effect size (LEfSe) analysis with samples presenting more than 0.1% relative abundance and found >30% of all samples. The LEfSe analysis employed the nonparametric factorial Kruskal–Wallis test or Wilcoxon rank-sum test and LDA to identify differentially abundant taxa between the groups. Only taxa with LDA score greater than two or less than two at a p-value < 0.05 were considered significantly enriched. All statistical tests are two-tailed, and a p-value < 0.05 was considered statistically significant. The random forest method [35] was performed to determine a ranked list of all bacterial taxa to identify the most predictive bacterial community to classify β-blocker users and nonusers. The random forest is a supervised learning algorithm ranking OTUs based on their ability to discriminate among the groups, while accounting for the complex interrelationships in high dimensional data. The MetagenomeSeq method was also used to evaluate differential abundance in sparse marker-gene survey data using a zero-inflated Gaussian (ZIG) fit model to account for undersampling and sparsity in OTU count data after normalizing the data through cumulative sum scaling (CSS) [36]. Finally, the log-transformed read counts difference of the top selected genera from the ZIG fit model between β-blocker users and nonusers was analyzed in the full and PS-matched cohorts.
Co-occurrence and random forest analyses were performed by MicrobiomeAnalyst [33]. The other statistical analyses were performed using R statistical software (version 3.5.1) and STATA statistical software (version 14; StataCorp LLC, College Station, TX, USA).

2.6. Functional Prediction Analysis

Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt2) [37] was used to predict the metagenome, which was based on Integrated Microbial Genomes (IMG) database [38] to evaluate the functions of gut microbiota among β blocker users and nonusers in the full cohort and PS-matched cohort. An OTU table was used for predicting metagenome based on Kyoto Encyclopedia of Genes and Genomes (KEGG) Orthology (KO) annotations. Metabolic module enrichment analysis was done with functional sets enrichment analysis (FSEA) described by Liu et al. [39]. The ‘FSEA’ function in the MARco R package based on the Liu et al. paper was applied in this study [39]. The ‘FSEA’ was embedded with the ‘gage’ R-package [40]. Enrichment scores were scored based on the GSEA algorithm of the Database for Annotation, Visualization, and Integrated Discovery (DAVID) bioinformatics resources [41,42].

3. Results

3.1. Patient Characteristics

Patient characteristics are shown in Table 1, with those receiving β-blockers having a higher proportion of diabetes, hypertension, dyslipidemia, coronary artery disease, heart failure, cerebrovascular disease, and more commonly used ACEI/ARB, glucose-lowering drugs (such as dipeptidyl peptidase-4 inhibitors or insulin) and statin. PS matching resulted in 62 matched pairs with balanced baseline characteristics (Table 1).

3.2. Gut Microbiota Profile Differs in Hemodialysis Patients with and without β Blocker Treatment

The rarefaction curves were close to asymptotic based on the number of OTUs observed. To represent the microbiome community with enough coverage, the rarefaction curves reached saturation at a cutoff point of 45,000 sequences per sample (Supplementary Figure S1). Compared to the gut microbiota composition and structure between β-blocker users and nonusers, no substantial differences were observed in the relative abundance proportion in the full and PS-matched cohorts (Supplementary Figure S2). Hemodialysis patients taking β-blockers had a higher α-diversity and a distinct β-diversity compared to nonusers in the full and PS-matched cohorts (Figure 2). The core microbiome was Bacteroides in hemodialysis patients (Supplementary Figure S3A), with a similar core microbiome in β-blocker users and nonusers (Supplementary Figure S3B).

3.3. Specific Microbial Taxa Differences between β-Blocker Users and Nonusers

Discriminant analysis using LEfSe identified the significant differentiating taxa between study groups. In the full cohort, the genera Ruminococcus 2, Collinsella, Ruminococcaceae UCG-004, Ruminiclostridium 5, Anaerotruncus, Eisenbergiella, and Flavonifractor were enriched in β-blocker users compared to nonusers (Figure 3A). In the PS-matched cohort, the enriched genera were Faecalibacterium, Subdoligranulum, Tyzzerella, Pantoea, Lachnospiraceae UCG-004, and Flavonifractor were found (Figure 3B). Using random forest models for taxonomy prediction, the top three ranked genera to discriminate between β-blocker users and nonusers were Parabacteroides, Flavonifractor, and Ruminococcaceae UCG-004 in the full cohort (Figure 4A), Prevotella 9, Flavonifractor, and Tyzzerella in the PS-matched cohort (Figure 4B).
To reduce the effect of zero-inflation in the microbiome data, we performed the MetagenomeSeq algorithm integrating the CSS method and a statistical model based on the ZIG distribution. Evaluating the significant difference in genus taxonomy between β-blocker users and nonusers, we found eight genera differences in the full cohort and PS-matched cohort (Supplementary Table S1). There were three different genera (Flavonifractor, Tyzzerella, and Prevotellaceae NK3B31 group) in both the full and PS-matched cohorts (Figure 5A). Focusing on the ZIG fit model to predict specific genera, there was an increased Flavonifractor genus in β-blocker users compared to nonusers using a classical univariate test (Kruskal–Wallis test) in the full (p = 0.023) and PS-matched cohorts (p = 0.01) (Figure 5B). However, no differences were found in Tyzzerella or Prevotellaceae NK3B31 group (Figure 5B).
Using PICRUSt2 as a metagenome predictive exploratory tool, genes were categorized into KEGG Orthology metabolic pathways. All predicted KEGG Orthology (KOs) were mapped to KEGG metabolic pathways. Each pathway was tested with gene-set enrichment by comparing expected gene abundance between β blocker users and nonusers in full and PS-matched cohorts. However, no significant KEGG enriched pathways were observed (Figures S4 and S5).

4. Discussion

In the present study, hemodialysis patients treated with β-blockers had a higher α-diversity and a distinct β-diversity compared to nonusers. The microbial communities contained higher levels of Bacteroidetes and lower levels of Firmicutes in all hemodialysis patients, which is similar to CKD rat microbial communities [43] and in a human CKD microbiota study [44]. Co-occurrence analysis revealed no difference in keystone taxa Bacteroides between β-blocker users and nonusers. Overall, there was an enriched genus Flavonifractor in β-blocker users in the full and PS-matched cohorts. Furthermore, LEfSe analysis, random forest algorithm, ZIG fit model, and univariate test all confirmed this difference between groups. However, we did not determine KEGG metabolic pathways between β blocker users and nonusers using PICRUSt2 functional prediction analysis.
β-blocker use was associated with a higher α-diversity than nonusers in hemodialysis patients, which was linked to a favorable healthy state [45]. Increased α-diversity has been associated with foods generally considered healthy, such as plant consumption or red wine [46,47,48]. Furthermore, commonly used medications such as antibiotics or proton pump inhibitors can decrease gut α-diversity [49]. Regarding the specific taxonomy of the gut microbiome, the genus Flavonifractor was enriched in β-blocker users in both the full and PS-matched cohorts. Flavonifractor is associated with several diseases, such as obesity [50], atrial fibrillation [51], coronary artery disease [52], and medications (antidiabetic drugs, such as Metformin and Glucagon-like peptide 1 Receptor agonist [53]). It can convert quercetin or other flavonoids into acetic acid and butyric acid [54] and is also correlated with oxidative stress and inflammation [55]. The presence of Flavonifractor was found in association with circulating inflammatory markers (i.e., interleukin-6, interleukin-8, interleukin-1β) [56], which were linked to cardiovascular disease. Besides, oral administration of Flavonifractor plautii was involved in the inhibition of tumor necrosis factor-α expression in obese adipose tissue inflammatory environments [57]. Thus, the increased abundance of Flavonifractor by β-blocker treatment may have a potential benefit in cardiovascular disease via gut microbiota regulation.
We also identified a potential link between β-blocker use and the genus Tyzzerella in the PS-matched cohort. Importantly, Tyzzerella was enriched in those with a high cardiovascular risk profile [58]. However, the small sample size limited the potential association between β-blocker use and Tyzzerella in univariate analysis, so more extensive studies are needed to confirm this association. Regarding the link between β-blocker and microbiota changes; a chimera mouse model suggested bone marrow beta1/2 adrenergic receptor signaling can regulate host-microbiota interactions, leading to the generation of novel anti-inflammatory treatments for gut dysbiosis [59]. Therefore, depletion of this sympathetic regulation in bone marrow promotes beneficial shifts in gut microbiota associated with gut immune suppression [59]. It is proposed that beta-blockers may provide a beneficial microbiome in such conditions.
We compared the microbiota differences between β-blocker users and nonusers using PS matching analysis in the present study. Since β-blocker intake is highly correlated with age, cardiovascular risk, comorbidities, and concurrent medication, each factor represents a relevant confounder for microbiome analyses [16,60]. Most observational studies have controlled for possible confounding variables, but even rigorous data adjustment cannot eliminate the risk of bias. PS matching is an alternative to reduce the effect of influencing factors on gut microbiota analysis [30,31]; thus, we selected variables of interest as potential confounders and then performed PS matching to reduce these effects deviations and confounding variables to conduct a reasonable comparison between groups. The intestinal microbiota was affected by various factors, including demographic data, comorbidities, concomitant medications, and clinical laboratory data, and the application of PS matching eliminated confounding factors. Using PS analysis, there was still a higher α-diversity and different β-diversity in β-blocker users compared to nonusers. We also identified six genera associated explicitly with the β-blocker user in the LEfSe analysis, four top-ranked genera in random forest analysis, and eight genera in ZIG fit model analysis. Although there were some differences in bacterial associations with β-blocker use in our full (before PS matching) and PS-matched cohorts, we investigated the taxa represented in both the full and PS-matched cohorts. Importantly, three genera (Flavonifractor, Tyzzerella, and Prevotellaceae NK3B31 group) were both significant differences in ZIG fit model among the full and PS-matched cohorts. The genera abundance differences between β-blocker users and nonusers were changed in the PS matching procedure. The genera abundance significant differences in Ruminiclostridium 9, Ruminococcaceae UCG-004, Anaerotruncus, and Ruminiclostridium 5 were attenuated after PS matching, suggesting that these genera abundances may be more strongly associated with other confounding variables, such as comorbidities or concomitant medications, which was accounted for in the PS models. In specific, genus Anaerotruncus was reported related to hypertension [61,62,63], diabetes mellitus [64], and dipeptidyl peptidase 4 inhibitors used [65], which were unbalance in the pre-matched cohort. The genus Ruminiclostridium and Ruminococcaceae were correlated to hypertension in the previous study [66,67,68]. Thus, the change of gut microbiome difference in β-blockers users and nonusers before and after PS matching reflected confounders’ influence. Since many factors influencing gut microbiota, we performed PS matching as an alternative technique to account for multiple confounders in this study.
In addition, there were more zeros than expected under the assumption of Poisson or negative binomial distributions for microbiome OTU counts, known as zero-inflation. One popular strategy to circumvent the zero-inflation problem is to add a pseudo-count [69]; however, this assumption may not be appropriate due to the large extent of structural zeros due to physical absence. Moreover, the pseudo-count choice is arbitrary, and the clustering results can be highly dependent upon the choice [70]. Thus, CSS was developed for microbiome sequencing data, and a zero-inflated model was used to model read counts that have an excess of zeros [36,71,72]. In CSS, raw counts are divided by the cumulative sum of counts up to a percentile determined using a data-driven approach to capture the relatively invariant count distribution for a dataset. To solve the zero-inflation issue, we applied the ZIG fit model and calculated the CSS. Interestingly, the four genera in the full (Ruminococcaceae UCG-004, Ruminiclostridium 5, Anaerotruncus, and Flavonifractor) and PS-matched cohorts (Flavonifractor, Tyzzerella, Faecalibacterium, Subdoligranulum) overlapped in the LEfSe analysis and ZIG fit model analysis.
Several limitations should be mentioned. First, cross-sectional studies only provide an impression of the relative abundance of bacterial taxa at a single time point, so causal inference cannot be addressed. Besides, the observational study only demonstrates the association rather than the causality. Second, the microbiota was assessed with a fecal sample, which may differ from microbiota from other parts of the intestine. Besides, 16S rRNA sequencing is limited as it cannot differentiate viable from non-viable bacteria. A significant portion of the taxa identified by sequencing may not be metabolically active. Thus, further study is needed to investigate various samples, such as small intestine or colon mucosal bacteria. Third, PS matching might not fully balance the overall effects of medications or disease severity, such as the dose of medications or the status between controlled and uncontrolled DM. Finally, the study was performed in Asia hemodialysis patients whose diet is different from other populations, so dietary effects on the gut microbiome should be interpreted with caution.

5. Conclusions

This study demonstrated that the composition of the gut microbiota was different in hemodialysis patients treated with β-blockers, with a higher level of α-diversity and genus Flavonifractor. These findings support the additional benefits of β-blocker treatment, which may mediate the microbiota in hemodialysis patients. However, the functional relevance of the β-blocker induced microbial differences is unclear. Hence, larger prospective treatment naïve studies are warranted to understand the impact of β-blockers on the gut microbiome of CKD patients and their implications for health and disease.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2075-4426/11/3/198/s1, Figure S1: Rarefaction curves of the number of OTUs versus the sequencing effort per sample in the full cohort. Figure S2: The relative abundance percentage of intestinal microbiota between β-blocker users and nonusers in the full cohort and propensity score matching cohort. (A) Phylum level (B) Class level (C) Order level. Figure S3: Core microbiome analysis in hemodialysis patients with and without β-blocker used. (A) SparCC correlation analysis (genus level using 100 SparCC permutations, 0.35 correlation threshold, and 0.05 p-value threshold) in all hemodialysis patients with and without β-blocker used (B) Relative abundance and sample prevalence of bacterial genus in β-blocker users and nonusers. Figure S4: Enrichment analysis for predictive Kyoto Encyclopedia of Genes and Genomes (KEGG) metabolic modules between β-blocker users and nonusers in full (before propensity score matching) cohort. No significant KEGG enriched pathways were observed (all p-value > 0.05). Figure S5: Enrichment analysis for predictive Kyoto Encyclopedia of Genes and Genomes (KEGG) metabolic modules between β-blocker users and nonusers in propensity score-matched cohort. No significant KEGG enriched pathways were observed (all p-value > 0.05). Table S1: Summary table of significant genus difference in hemodialysis patients with and without β-blocker treatment in zero-inflated Gaussian fit model.

Author Contributions

Conceptualization, P.-H.W., Y.-T.L., Y.-W.C., and M.-C.K.; methodology, P.-H.W., Y.-T.L., P.-Y.L., Y.-W.C., and M.-C.K.; formal analysis, P.-H.W., Y.-T.L., and P.-Y.L.; investigation, P.-H.W., Y.-T.L., P.-Y.L.; data curation, P.-H.W., Y.-T.L., and P.-Y.L.; writing—original draft preparation, P.-H.W., Y.-T.L., and P.-Y.L.; writing—review and editing, W.-C.H., W.-C.T., Y.-C.T., Y.-W.C., M.-C.K., T.-Y.L., S.-C.H., R.A.D., and Y.-S.C.; funding acquisition, P.-H.W. and Y.-T.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by grants from the Ministry of Science and Technology, Taiwan (MOST 107-2314-B-037-104, MOST 107-2314-B-037-098-MY3, and MOST 109-2314-B-037-088), Kaohsiung Medical University Hospital, Taiwan (KMUH108-8M11, KMUH108-8R70, KMUH109-9R16, KMUH109-9R81, KMUH-DK(B)110003-1, and KMUH-DK(B)110003-4), Kaohsiung Medical University, Taiwan (KMU-Q108024, KMU-Q108027, KMU-DK(B)110003, and KMU-DK(B)110003-2), and NSYSU-KMU JOINT RESEARCH PROJECT (NSYSUKMU 105-I005 and NSYSUKMU 106-I005).

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Ethics Committee of Kaohsiung Medical University Hospital (KMUHIRB-E(I)-20160095 and KMUHIRB-E(I)-20180118) and Taipei Tzu Chi Hospital (07-X01-002).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Raw FASTQ data files of samples from hemodialysis patients have been deposited in the NCBI Sequence Read Archive database under BioProject accession number PRJNA648014.

Acknowledgments

We acknowledge Germark Biotechnology Lab for the stool sample analysis of gut microbiota.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lynch, S.V.; Pedersen, O. The Human Intestinal Microbiome in Health and Disease. N. Engl. J. Med. 2016, 375, 2369–2379. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Jandhyala, S.M.; Talukdar, R.; Subramanyam, C.; Vuyyuru, H.; Sasikala, M.; Nageshwar Reddy, D. Role of the normal gut microbiota. World J. Gastroenterol. 2015, 21, 8787–8803. [Google Scholar] [CrossRef]
  3. Salvucci, E. The human-microbiome superorganism and its modulation to restore health. Int. J. Food Sci. Nutr. 2019, 70, 781–795. [Google Scholar] [CrossRef]
  4. Belizario, J.E.; Faintuch, J.; Garay-Malpartida, M. Gut Microbiome Dysbiosis and Immunometabolism: New Frontiers for Treatment of Metabolic Diseases. Mediators Inflamm. 2018, 2018, 2037838. [Google Scholar] [CrossRef] [PubMed]
  5. DeGruttola, A.K.; Low, D.; Mizoguchi, A.; Mizoguchi, E. Current Understanding of Dysbiosis in Disease in Human and Animal Models. Inflamm. Bowel. Dis. 2016, 22, 1137–1150. [Google Scholar] [CrossRef] [Green Version]
  6. Carding, S.; Verbeke, K.; Vipond, D.T.; Corfe, B.M.; Owen, L.J. Dysbiosis of the gut microbiota in disease. Microb. Ecol. Health Dis. 2015, 26, 26191. [Google Scholar] [CrossRef]
  7. Jackson, M.A.; Verdi, S.; Maxan, M.E.; Shin, C.M.; Zierer, J.; Bowyer, R.C.E.; Martin, T.; Williams, F.M.K.; Menni, C.; Bell, J.T.; et al. Gut microbiota associations with common diseases and prescription medications in a population-based cohort. Nat. Commun. 2018, 9, 2655. [Google Scholar] [CrossRef] [Green Version]
  8. Imhann, F.; Bonder, M.J.; Vich Vila, A.; Fu, J.; Mujagic, Z.; Vork, L.; Tigchelaar, E.F.; Jankipersadsing, S.A.; Cenit, M.C.; Harmsen, H.J.; et al. Proton pump inhibitors affect the gut microbiome. Gut 2016, 65, 740–748. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Makivuokko, H.; Tiihonen, K.; Tynkkynen, S.; Paulin, L.; Rautonen, N. The effect of age and non-steroidal anti-inflammatory drugs on human intestinal microbiota composition. Br. J. Nutr. 2010, 103, 227–234. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Bahr, S.M.; Tyler, B.C.; Wooldridge, N.; Butcher, B.D.; Burns, T.L.; Teesch, L.M.; Oltman, C.L.; Azcarate-Peril, M.A.; Kirby, J.R.; Calarge, C.A. Use of the second-generation antipsychotic, risperidone, and secondary weight gain are associated with an altered gut microbiota in children. Transl. Psychiatry 2015, 5, e652. [Google Scholar] [CrossRef]
  11. Forslund, K.; Hildebrand, F.; Nielsen, T.; Falony, G.; Le Chatelier, E.; Sunagawa, S.; Prifti, E.; Vieira-Silva, S.; Gudmundsdottir, V.; Pedersen, H.K.; et al. Disentangling type 2 diabetes and metformin treatment signatures in the human gut microbiota. Nature 2015, 528, 262–266. [Google Scholar] [CrossRef]
  12. Freedberg, D.E.; Toussaint, N.C.; Chen, S.P.; Ratner, A.J.; Whittier, S.; Wang, T.C.; Wang, H.H.; Abrams, J.A. Proton Pump Inhibitors Alter Specific Taxa in the Human Gastrointestinal Microbiome: A Crossover Trial. Gastroenterology 2015, 149, 883–885.e9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Jackson, M.A.; Goodrich, J.K.; Maxan, M.E.; Freedberg, D.E.; Abrams, J.A.; Poole, A.C.; Sutter, J.L.; Welter, D.; Ley, R.E.; Bell, J.T.; et al. Proton pump inhibitors alter the composition of the gut microbiota. Gut 2016, 65, 749–756. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Fu, J.; Bonder, M.J.; Cenit, M.C.; Tigchelaar, E.F.; Maatman, A.; Dekens, J.A.; Brandsma, E.; Marczynska, J.; Imhann, F.; Weersma, R.K.; et al. The Gut Microbiome Contributes to a Substantial Proportion of the Variation in Blood Lipids. Circ. Res. 2015, 117, 817–824. [Google Scholar] [CrossRef] [PubMed]
  15. Buffie, C.G.; Pamer, E.G. Microbiota-mediated colonization resistance against intestinal pathogens. Nat. Rev. Immunol. 2013, 13, 790–801. [Google Scholar] [CrossRef] [Green Version]
  16. Zhernakova, A.; Kurilshikov, A.; Bonder, M.J.; Tigchelaar, E.F.; Schirmer, M.; Vatanen, T.; Mujagic, Z.; Vila, A.V.; Falony, G.; Vieira-Silva, S.; et al. Population-based metagenomics analysis reveals markers for gut microbiome composition and diversity. Science 2016, 352, 565–569. [Google Scholar] [CrossRef] [Green Version]
  17. Evenepoel, P.; Poesen, R.; Meijers, B. The gut-kidney axis. Pediatr. Nephrol. 2017, 32, 2005–2014. [Google Scholar] [CrossRef] [PubMed]
  18. Vaziri, N.D.; Wong, J.; Pahl, M.; Piceno, Y.M.; Yuan, J.; DeSantis, T.Z.; Ni, Z.; Nguyen, T.H.; Andersen, G.L. Chronic kidney disease alters intestinal microbial flora. Kidney Int. 2013, 83, 308–315. [Google Scholar] [CrossRef] [Green Version]
  19. Weir, M.A.; Herzog, C.A. Beta blockers in patients with end-stage renal disease-Evidence-based recommendations. Semin. Dial. 2018, 31, 219–225. [Google Scholar] [CrossRef]
  20. Bakris, G.L.; Hart, P.; Ritz, E. Beta blockers in the management of chronic kidney disease. Kidney Int. 2006, 70, 1905–1913. [Google Scholar] [CrossRef] [Green Version]
  21. Nemet, I.; Saha, P.P.; Gupta, N.; Zhu, W.; Romano, K.A.; Skye, S.M.; Cajka, T.; Mohan, M.L.; Li, L.; Wu, Y.; et al. A Cardiovascular Disease-Linked Gut Microbial Metabolite Acts via Adrenergic Receptors. Cell 2020, 180, 862–877.e22. [Google Scholar] [CrossRef]
  22. Ikizler, T.A.; Burrowes, J.D.; Byham-Gray, L.D.; Campbell, K.L.; Carrero, J.J.; Chan, W.; Fouque, D.; Friedman, A.N.; Ghaddar, S.; Goldstein-Fuchs, D.J.; et al. KDOQI Clinical Practice Guideline for Nutrition in CKD: 2020 Update. Am. J. Kidney Dis. 2020, 76, S1–S107. [Google Scholar] [CrossRef] [PubMed]
  23. Herlemann, D.P.; Labrenz, M.; Jurgens, K.; Bertilsson, S.; Waniek, J.J.; Andersson, A.F. Transitions in bacterial communities along the 2000 km salinity gradient of the Baltic Sea. ISME J. 2011, 5, 1571–1579. [Google Scholar] [CrossRef] [Green Version]
  24. Comeau, A.M.; Douglas, G.M.; Langille, M.G. Microbiome Helper: A Custom and Streamlined Workflow for Microbiome Research. mSystems 2017, 2, e00127-16. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Zhang, J.; Kobert, K.; Flouri, T.; Stamatakis, A. PEAR: A fast and accurate Illumina Paired-End reAd mergeR. Bioinformatics 2014, 30, 614–620. [Google Scholar] [CrossRef] [Green Version]
  26. Caporaso, J.G.; Kuczynski, J.; Stombaugh, J.; Bittinger, K.; Bushman, F.D.; Costello, E.K.; Fierer, N.; Pena, A.G.; Goodrich, J.K.; Gordon, J.I.; et al. QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 2010, 7, 335–336. [Google Scholar] [CrossRef] [Green Version]
  27. Quast, C.; Pruesse, E.; Yilmaz, P.; Gerken, J.; Schweer, T.; Yarza, P.; Peplies, J.; Glockner, F.O. The SILVA ribosomal RNA gene database project: Improved data processing and web-based tools. Nucleic. Acids. Res. 2013, 41, D590–D596. [Google Scholar] [CrossRef]
  28. Yilmaz, P.; Parfrey, L.W.; Yarza, P.; Gerken, J.; Pruesse, E.; Quast, C.; Schweer, T.; Peplies, J.; Ludwig, W.; Glockner, F.O. The SILVA and “All-species Living Tree Project (LTP)” taxonomic frameworks. Nucleic. Acids. Res. 2014, 42, D643–D648. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Edgar, R.C. Search and clustering orders of magnitude faster than BLAST. Bioinformatics 2010, 26, 2460–2461. [Google Scholar] [CrossRef] [Green Version]
  30. Sturmer, T.; Wyss, R.; Glynn, R.J.; Brookhart, M.A. Propensity scores for confounder adjustment when assessing the effects of medical interventions using nonexperimental study designs. J. Intern. Med. 2014, 275, 570–580. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Ali, M.S.; Groenwold, R.H.; Belitser, S.V.; Pestman, W.R.; Hoes, A.W.; Roes, K.C.; Boer, A.; Klungel, O.H. Reporting of covariate selection and balance assessment in propensity score analysis is suboptimal: A systematic review. J. Clin. Epidemiol. 2015, 68, 112–121. [Google Scholar] [CrossRef]
  32. Lozupone, C.A.; Hamady, M.; Kelley, S.T.; Knight, R. Quantitative and qualitative beta diversity measures lead to different insights into factors that structure microbial communities. Appl. Environ. Microbiol. 2007, 73, 1576–1585. [Google Scholar] [CrossRef] [Green Version]
  33. Chong, J.; Liu, P.; Zhou, G.; Xia, J. Using MicrobiomeAnalyst for comprehensive statistical, functional, and meta-analysis of microbiome data. Nat. Protoc. 2020, 15, 799–821. [Google Scholar] [CrossRef] [PubMed]
  34. Friedman, J.; Alm, E.J. Inferring correlation networks from genomic survey data. PLoS Comput. Biol. 2012, 8, e1002687. [Google Scholar] [CrossRef] [Green Version]
  35. Svetnik, V.; Liaw, A.; Tong, C.; Culberson, J.C.; Sheridan, R.P.; Feuston, B.P. Random forest: A classification and regression tool for compound classification and QSAR modeling. J. Chem. Inf. Comput. Sci. 2003, 43, 1947–1958. [Google Scholar] [CrossRef] [PubMed]
  36. Paulson, J.N.; Stine, O.C.; Bravo, H.C.; Pop, M. Differential abundance analysis for microbial marker-gene surveys. Nat. Methods 2013, 10, 1200–1202. [Google Scholar] [CrossRef] [Green Version]
  37. Douglas, G.M.; Maffei, V.J.; Zaneveld, J.R.; Yurgel, S.N.; Brown, J.R.; Taylor, C.M.; Huttenhower, C.; Langille, M.G.I. PICRUSt2 for prediction of metagenome functions. Nat. Biotechnol. 2020, 38, 685–688. [Google Scholar] [CrossRef]
  38. Markowitz, V.M.; Chen, I.M.; Palaniappan, K.; Chu, K.; Szeto, E.; Grechkin, Y.; Ratner, A.; Jacob, B.; Huang, J.; Williams, P.; et al. IMG: The Integrated Microbial Genomes database and comparative analysis system. Nucleic. Acids. Res. 2012, 40, D115–D122. [Google Scholar] [CrossRef] [Green Version]
  39. Liu, P.-Y. poyuliu/MARco: MARco: Microbiome Analysis RcodeDB (Version v1.0). Zenodo. Available online: http://0-doi-org.brum.beds.ac.uk/10.5281/zenodo.4589898 (accessed on 1 February 2021).
  40. Luo, W.; Brouwer, C. Pathview: An R/Bioconductor package for pathway-based data integration and visualization. Bioinformatics 2013, 29, 1830–1831. [Google Scholar] [CrossRef] [Green Version]
  41. Huang da, W.; Sherman, B.T.; Lempicki, R.A. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc. 2009, 4, 44–57. [Google Scholar] [CrossRef]
  42. Huang, D.W.; Sherman, B.T.; Tan, Q.; Collins, J.R.; Alvord, W.G.; Roayaei, J.; Stephens, R.; Baseler, M.W.; Lane, H.C.; Lempicki, R.A. The DAVID Gene Functional Classification Tool: A novel biological module-centric algorithm to functionally analyze large gene lists. Genome Biol. 2007, 8, R183. [Google Scholar] [CrossRef] [Green Version]
  43. Lau, W.L.; Vaziri, N.D.; Nunes, A.C.F.; Comeau, A.M.; Langille, M.G.I.; England, W.; Khazaeli, M.; Suematsu, Y.; Phan, J.; Whiteson, K. The Phosphate Binder Ferric Citrate Alters the Gut Microbiome in Rats with Chronic Kidney Disease. J. Pharmacol. Exp. Ther. 2018, 367, 452–460. [Google Scholar] [CrossRef] [Green Version]
  44. Lun, H.; Yang, W.; Zhao, S.; Jiang, M.; Xu, M.; Liu, F.; Wang, Y. Altered gut microbiota and microbial biomarkers associated with chronic kidney disease. Microbiologyopen 2019, 8, e00678. [Google Scholar] [CrossRef]
  45. Le Chatelier, E.; Nielsen, T.; Qin, J.; Prifti, E.; Hildebrand, F.; Falony, G.; Almeida, M.; Arumugam, M.; Batto, J.M.; Kennedy, S.; et al. Richness of human gut microbiome correlates with metabolic markers. Nature 2013, 500, 541–546. [Google Scholar] [CrossRef]
  46. McDonald, D.; Hyde, E.; Debelius, J.W.; Morton, J.T.; Gonzalez, A.; Ackermann, G.; Aksenov, A.A.; Behsaz, B.; Brennan, C.; Chen, Y.; et al. American Gut: An Open Platform for Citizen Science Microbiome Research. mSystems 2018, 3, e00031-18. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Le Roy, C.I.; Wells, P.M.; Si, J.; Raes, J.; Bell, J.T.; Spector, T.D. Red Wine Consumption Associated With Increased Gut Microbiota alpha-Diversity in 3 Independent Cohorts. Gastroenterology 2020, 158, 270–272. [Google Scholar] [CrossRef]
  48. Leeming, E.R.; Johnson, A.J.; Spector, T.D.; Le Roy, C.I. Effect of Diet on the Gut Microbiota: Rethinking Intervention Duration. Nutrients 2019, 11, 2862. [Google Scholar] [CrossRef] [Green Version]
  49. Le Bastard, Q.; Al-Ghalith, G.A.; Gregoire, M.; Chapelet, G.; Javaudin, F.; Dailly, E.; Batard, E.; Knights, D.; Montassier, E. Systematic review: Human gut dysbiosis induced by non-antibiotic prescription medications. Aliment. Pharmacol. Ther. 2018, 47, 332–345. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Kasai, C.; Sugimoto, K.; Moritani, I.; Tanaka, J.; Oya, Y.; Inoue, H.; Tameda, M.; Shiraki, K.; Ito, M.; Takei, Y.; et al. Comparison of the gut microbiota composition between obese and non-obese individuals in a Japanese population, as analyzed by terminal restriction fragment length polymorphism and next-generation sequencing. BMC Gastroenterol. 2015, 15, 100. [Google Scholar] [CrossRef] [Green Version]
  51. Zuo, K.; Li, J.; Li, K.; Hu, C.; Gao, Y.; Chen, M.; Hu, R.; Liu, Y.; Chi, H.; Wang, H.; et al. Disordered gut microbiota and alterations in metabolic patterns are associated with atrial fibrillation. Gigascience 2019, 8, giz058. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Zhu, Q.; Gao, R.; Zhang, Y.; Pan, D.; Zhu, Y.; Zhang, X.; Yang, R.; Jiang, R.; Xu, Y.; Qin, H. Dysbiosis signatures of gut microbiota in coronary artery disease. Physiol. Genom. 2018, 50, 893–903. [Google Scholar] [CrossRef] [PubMed]
  53. Gurung, M.; Li, Z.; You, H.; Rodrigues, R.; Jump, D.B.; Morgun, A.; Shulzhenko, N. Role of gut microbiota in type 2 diabetes pathophysiology. EBioMedicine 2020, 51, 102590. [Google Scholar] [CrossRef] [Green Version]
  54. Carlier, J.P.; Bedora-Faure, M.; K’Ouas, G.; Alauzet, C.; Mory, F. Proposal to unify Clostridium orbiscindens Winter et al. 1991 and Eubacterium plautii (Seguin 1928) Hofstad and Aasjord 1982, with description of Flavonifractor plautii gen. nov., comb. nov., and reassignment of Bacteroides capillosus to Pseudoflavonifractor capillosus gen. nov., comb. nov. Int. J. Syst. Evol. Microbiol. 2010, 60, 585–590. [Google Scholar] [CrossRef] [Green Version]
  55. Coello, K.; Hansen, T.H.; Sorensen, N.; Munkholm, K.; Kessing, L.V.; Pedersen, O.; Vinberg, M. Gut microbiota composition in patients with newly diagnosed bipolar disorder and their unaffected first-degree relatives. Brain Behav. Immun. 2019, 75, 112–118. [Google Scholar] [CrossRef] [PubMed]
  56. Huang, S.; Mao, J.; Zhou, L.; Xiong, X.; Deng, Y. The imbalance of gut microbiota and its correlation with plasma inflammatory cytokines in pemphigus vulgaris patients. Scand. J. Immunol. 2019, 90, e12799. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Mikami, A.; Ogita, T.; Namai, F.; Shigemori, S.; Sato, T.; Shimosato, T. Oral administration of Flavonifractor plautii attenuates inflammatory responses in obese adipose tissue. Mol. Biol. Rep. 2020, 47, 6717–6725. [Google Scholar] [CrossRef] [PubMed]
  58. Kelly, T.N.; Bazzano, L.A.; Ajami, N.J.; He, H.; Zhao, J.; Petrosino, J.F.; Correa, A.; He, J. Gut Microbiome Associates With Lifetime Cardiovascular Disease Risk Profile Among Bogalusa Heart Study Participants. Circ. Res. 2016, 119, 956–964. [Google Scholar] [CrossRef] [Green Version]
  59. Yang, T.; Ahmari, N.; Schmidt, J.T.; Redler, T.; Arocha, R.; Pacholec, K.; Magee, K.L.; Malphurs, W.; Owen, J.L.; Krane, G.A.; et al. Shifts in the Gut Microbiota Composition Due to Depleted Bone Marrow Beta Adrenergic Signaling Are Associated with Suppressed Inflammatory Transcriptional Networks in the Mouse Colon. Front. Physiol. 2017, 8, 220. [Google Scholar] [CrossRef] [Green Version]
  60. Falony, G.; Joossens, M.; Vieira-Silva, S.; Wang, J.; Darzi, Y.; Faust, K.; Kurilshikov, A.; Bonder, M.J.; Valles-Colomer, M.; Vandeputte, D.; et al. Population-level analysis of gut microbiome variation. Science 2016, 352, 560–564. [Google Scholar] [CrossRef]
  61. Dan, X.; Mushi, Z.; Baili, W.; Han, L.; Enqi, W.; Huanhu, Z.; Shuchun, L. Differential Analysis of Hypertension-Associated Intestinal Microbiota. Int. J. Med. Sci. 2019, 16, 872–881. [Google Scholar] [CrossRef] [Green Version]
  62. Robles-Vera, I.; Toral, M.; Duarte, J. Microbiota and Hypertension: Role of the Sympathetic Nervous System and the Immune System. Am. J. Hypertens. 2020, 33, 890–901. [Google Scholar] [CrossRef] [PubMed]
  63. Palmu, J.; Lahti, L.; Niiranen, T. Targeting Gut Microbiota to Treat Hypertension: A Systematic Review. Int. J. Environ. Res. Public Health 2021, 18, 1248. [Google Scholar] [CrossRef]
  64. Liu, C.; Shao, W.; Gao, M.; Liu, J.; Guo, Q.; Jin, J.; Meng, F. Changes in intestinal flora in patients with type 2 diabetes on a low-fat diet during 6 months of follow-up. Exp. Ther. Med. 2020, 20, 40. [Google Scholar] [CrossRef] [PubMed]
  65. Whang, A.; Nagpal, R.; Yadav, H. Bi-directional drug-microbiome interactions of anti-diabetics. EBioMedicine 2019, 39, 591–602. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Zuo, K.; Li, J.; Xu, Q.; Hu, C.; Gao, Y.; Chen, M.; Hu, R.; Liu, Y.; Chi, H.; Yin, Q.; et al. Dysbiotic gut microbes may contribute to hypertension by limiting vitamin D production. Clin. Cardiol. 2019, 42, 710–719. [Google Scholar] [CrossRef] [Green Version]
  67. Durgan, D.J.; Ganesh, B.P.; Cope, J.L.; Ajami, N.J.; Phillips, S.C.; Petrosino, J.F.; Hollister, E.B.; Bryan, R.M., Jr. Role of the Gut Microbiome in Obstructive Sleep Apnea-Induced Hypertension. Hypertension 2016, 67, 469–474. [Google Scholar] [CrossRef] [Green Version]
  68. Li, J.; Zhao, F.; Wang, Y.; Chen, J.; Tao, J.; Tian, G.; Wu, S.; Liu, W.; Cui, Q.; Geng, B.; et al. Gut microbiota dysbiosis contributes to the development of hypertension. Microbiome 2017, 5, 14. [Google Scholar] [CrossRef] [Green Version]
  69. Mandal, S.; Van Treuren, W.; White, R.A.; Eggesbo, M.; Knight, R.; Peddada, S.D. Analysis of composition of microbiomes: A novel method for studying microbial composition. Microb. Ecol. Health Dis. 2015, 26, 27663. [Google Scholar] [CrossRef] [Green Version]
  70. Costea, P.I.; Zeller, G.; Sunagawa, S.; Bork, P. A fair comparison. Nat. Methods 2014, 11, 359. [Google Scholar] [CrossRef]
  71. Jonsson, V.; Osterlund, T.; Nerman, O.; Kristiansson, E. Modelling of zero-inflation improves inference of metagenomic gene count data. Stat. Methods Med. Res. 2019, 28, 3712–3728. [Google Scholar] [CrossRef]
  72. Xu, L.; Paterson, A.D.; Turpin, W.; Xu, W. Assessment and Selection of Competing Models for Zero-Inflated Microbiome Data. PLoS ONE 2015, 10, e0129606. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Jpm 11 00198 g001 550
Figure 1. Study design.
Figure 1. Study design.
Jpm 11 00198 g001
Jpm 11 00198 g002 550
Figure 2. The α-diversity and β-diversity in hemodialysis patients with and without β blocker used in full cohort (A,B) and propensity score matching cohort (C,D). β blocker users had a higher α-diversity than β blocker nonusers in full cohort (A) and propensity score matching cohort (C) β blocker users had a different β-diversity (Bray–Curtis index) compared to β blocker nonusers in full cohort (B) and propensity score matching cohort (D). The β-diversity p-value was calculated using the homogeneity of group dispersions by the Permutational Analysis of Multivariate Dispersions (PERMDISP) test.
Figure 2. The α-diversity and β-diversity in hemodialysis patients with and without β blocker used in full cohort (A,B) and propensity score matching cohort (C,D). β blocker users had a higher α-diversity than β blocker nonusers in full cohort (A) and propensity score matching cohort (C) β blocker users had a different β-diversity (Bray–Curtis index) compared to β blocker nonusers in full cohort (B) and propensity score matching cohort (D). The β-diversity p-value was calculated using the homogeneity of group dispersions by the Permutational Analysis of Multivariate Dispersions (PERMDISP) test.
Jpm 11 00198 g002
Jpm 11 00198 g003 550
Figure 3. Taxonomic differences were detected between β blocker users and nonusers in the full cohort (A) and propensity score matching cohort (B). Linear discriminative analysis (LDA) effect size (LEfSe) analysis between β blocker users (red) and nonusers (blue) with an LDA score > 2.0 or < −2 with p-value > 0.1 among β blocker users and nonusers.
Figure 3. Taxonomic differences were detected between β blocker users and nonusers in the full cohort (A) and propensity score matching cohort (B). Linear discriminative analysis (LDA) effect size (LEfSe) analysis between β blocker users (red) and nonusers (blue) with an LDA score > 2.0 or < −2 with p-value > 0.1 among β blocker users and nonusers.
Jpm 11 00198 g003
Jpm 11 00198 g004 550
Figure 4. Determination of specific bacteria for discriminatory across hemodialysis patients with and without β blocker treatment in full cohort (A) and propensity score matching cohort (B). The discriminatory taxa were determined by applying Random Forest analysis using the genus-level abundance.
Figure 4. Determination of specific bacteria for discriminatory across hemodialysis patients with and without β blocker treatment in full cohort (A) and propensity score matching cohort (B). The discriminatory taxa were determined by applying Random Forest analysis using the genus-level abundance.
Jpm 11 00198 g004
Jpm 11 00198 g005 550
Figure 5. The genera difference between β blocker users and nonusers in the full cohort and propensity score matching cohort using zero-inflated Gaussian fit model. (A) The Venn diagram showed the different significant genera in the full cohort and propensity score-matched cohort. (B) Univariate test between selected genera from zero-inflated Gaussian fit model. Significance was considered for p < 0.05.
Figure 5. The genera difference between β blocker users and nonusers in the full cohort and propensity score matching cohort using zero-inflated Gaussian fit model. (A) The Venn diagram showed the different significant genera in the full cohort and propensity score-matched cohort. (B) Univariate test between selected genera from zero-inflated Gaussian fit model. Significance was considered for p < 0.05.
Jpm 11 00198 g005
Table
Table 1. Baseline characteristics of hemodialysis patients with and without β blocker treatment.
Table 1. Baseline characteristics of hemodialysis patients with and without β blocker treatment.
Baseline CharacteristicsBefore Propensity Score MatchingAfter Propensity Score Matching
β-Blocker
Users
(N = 83)		β-Blocker
Nonusers
(N = 110)		p-Valueβ-Blocker
Users
(N = 62)		β-Blocker
Nonusers
(N = 62)		p-Value
Age (years)64.3 ± 11.465.4 ± 11.20.51164.7 ± 11.666.3 ± 11.80.446
Male49 (59.0%)57 (51.8%)0.31837 (59.7%)28 (45.2%0.106
Body mass index23.4 ± 3.2523.6 ± 3.910.70823.5 ± 3.3423.5 ± 3.930.988
Dialysis vintage (months)86.24 ± 56.5396.54 ± 63.210.24393.22 ± 57.6185.4 ± 55.670.444
Smoking history15 (18.1%)12 (10.9%)0.1569 (14.5%)6 (9.7%)0.409
Arteriovenous fistula75 (90.4%)99 (90.0%)0.93457 (91.9%)57 (91.9%)>0.999
Comorbidities
Diabetes mellitus45 (54.2%)34 (30.9%)0.00124 (38.7%)30 (48.4%)0.277
Hypertension80 (96.4%)87 (79.1%)<0.00159 (95.2%)59 (95.2%)>0.999
Dyslipidemia31 (37.3%)24 (21.8%)0.01816 (25.8%)15 (24.2%)0.836
Coronary artery disease34 (41.0%)22 (20.0%)0.00121 (33.9%)18 (29.0%)0.562
Heart failure22 (26.5%)15 (13.6%)0.02514 (22.6%)11 (17.7%)0.502
Cerebrovascular disease31 (37.3%)24 (21.8%)0.0185 (8.1%)8 (12.9%)0.379
Parathyroidectomy history7 (8.4%)18 (16.4%)0.1046 (9.7%)6 (9.7%)>0.999
Medications
ACEI/ARB29 (34.9%)24 (21.8%)0.04323 (37.1%)15 (24.2%)0.119
Glucose lowering drugs34 (41.0%)23 (20.9%)0.00320 (32.3%)19 (30.6%)0.847
Sulfonylurea14 (16.9%)13 (11.8%)0.3176 (9.7%)11 (17.7%)0.192
Dipeptidyl peptidase 4 inhibitors28 (33.7%)13 (11.8%)<0.00117 (27.4%)11 (17.7%)0.198
Insulin17 (20.5%)10 (9.1%)0.0249 (14.5%)8 (12.9%)0.794
Statin29 (34.9%)17 (15.5%)0.00217 (27.4%)12 (19.4%)0.289
Calcium carbonate67 (80.7%)94 (85.5%)0.38251 (82.3%)50 (80.6%)0.817
Proton pump inhibitors13 (15.7%)10 (9.1%)0.1639 (14.5%)7 (11.3%)0.592
Clinical laboratory data
Hemoglobin (g/dL)10.62 ± 1.1410.71 ± 1.410.65010.6 ± 1.0510.74 ± 1.490.555
Albumin (g/dL)3.52 ± 0.513.56 ± 0.460.5383.53 ± 0.463.54 ± 0.470.902
Total cholesterol (mg/dL)154.01 ± 33.75161.89 ± 33.620.109151.94 ± 33.57163.51 ± 35.300.064
Triglyceride (mg/dL)140.52 ± 103.77129.61 ± 90.350.437136.21 ± 105.99131.14 ± 95.510.780
High sensitivity CRP (mg/dL)2.15 ± 4.652.5 ± 4.210.5892.45 ± 5.232.21 ± 3.950.779
Sodium (mmol/L)136.92 ± 2.68137.07 ± 2.620.700137.19 ± 2.80136.64 ± 2.440.241
Potassium (mmol/L)4.73 ± 0.684.61 ± 0.620.1954.77 ± 0.664.65 ± 0.650.294
Total calcium (mg/dL)9.15 ± 0.869.29 ± 0.940.2779.19 ± 0.929.25 ± 0.860.683
Phosphate (mg/dL)5.08 ± 1.214.95 ± 1.240.4535.16 ± 1.155.09 ± 1.350.768
Parathyroid hormone (pg/mL)376.53 ± 338.79383.5 ± 278.130.876394.16 ± 370.62357.29 ± 245.840.515
Serum iron (μg/dL)63.57 ± 26.7365.85 ± 21.160.50863.94 ± 26.6167.52 ± 22.930.424
Ferritin (ng/mL)567.53 ± 549.64496.67 ± 377.330.291534.93 ± 330.67538.54 ± 413.540.957
nPCR (g/kg/day)1.12 ± 0.211.16 ± 0.270.3261.12 ± 0.201.18 ± 0.280.180
Single pool Kt/V1.67 ± 0.271.65 ± 0.270.5911.67 ± 0.281.68 ± 0.270.817
Dietary intake (serving/day)
Meat0.86 ± 0.570.82 ± 0.530.6520.86 ± 0.570.74 ± 0.520.241
Vegetable2.01 ± 1.091.86 ± 1.110.2652.05 ± 1.061.91 ± 1.180.499
Fruit0.93 ± 0.720.95 ± 0.720.5830.86 ± 0.630.89 ± 0.750.837
Bristol stool scale3.94 ± 1.863.74 ± 1.760.4484 ± 1.783.71 ± 1.670.352
Abbreviation: ACEI/ARB, angiotensin-converting enzyme inhibitors/angiotensin-receptor blockers; CRP, C reactive protein; nPCR, normalized protein catabolic rate.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lin, Y.-T.; Lin, T.-Y.; Hung, S.-C.; Liu, P.-Y.; Hung, W.-C.; Tsai, W.-C.; Tsai, Y.-C.; Delicano, R.A.; Chuang, Y.-S.; Kuo, M.-C.; et al. Differences in the Microbial Composition of Hemodialysis Patients Treated with and without β-Blockers. J. Pers. Med. 2021, 11, 198. https://0-doi-org.brum.beds.ac.uk/10.3390/jpm11030198

AMA Style

Lin Y-T, Lin T-Y, Hung S-C, Liu P-Y, Hung W-C, Tsai W-C, Tsai Y-C, Delicano RA, Chuang Y-S, Kuo M-C, et al. Differences in the Microbial Composition of Hemodialysis Patients Treated with and without β-Blockers. Journal of Personalized Medicine. 2021; 11(3):198. https://0-doi-org.brum.beds.ac.uk/10.3390/jpm11030198

Chicago/Turabian Style

Lin, Yi-Ting, Ting-Yun Lin, Szu-Chun Hung, Po-Yu Liu, Wei-Chun Hung, Wei-Chung Tsai, Yi-Chun Tsai, Rachel Ann Delicano, Yun-Shiuan Chuang, Mei-Chuan Kuo, and et al. 2021. "Differences in the Microbial Composition of Hemodialysis Patients Treated with and without β-Blockers" Journal of Personalized Medicine 11, no. 3: 198. https://0-doi-org.brum.beds.ac.uk/10.3390/jpm11030198

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop