The causal associations of circulating amino acids with blood pressure: a Mendelian randomization study
BMC Medicine volume 20, Article number: 414 (2022)
Circulating levels of amino acids were associated with blood pressure (BP) in observational studies. However, the causation of such associations has been hypothesized but is difficult to prove in human studies. Here, we aimed to use two-sample Mendelian randomization analyses to evaluate the potential causal associations of circulating levels of amino acids with BP and risk of hypertension.
We generated genetic instruments for circulating levels of nine amino acids by conducting meta-analyses of genome-wide association study (GWAS) in UK Biobank participants with metabolomic data (n = 98,317) and another published metabolomics GWAS (n = 24,925). Data on the associations of the genetic variants with BP and hypertension were obtained in the UK Biobank participants without metabolomic data (n = 286,390). The causal effects were estimated using inverse-variance weighted method.
Significant evidence consistently supported the causal effects of increased branched-chain amino acids (BCAAs, i.e., leucine, isoleucine, and valine) levels on higher BP and risk of hypertension (all P < 0.006 after Bonferroni correction except for Pleucine-on-diastolicBP = 0.008). For example, per standard deviation higher of genetically predicted isoleucine levels were associated with 2.71 ± 0.78 mmHg higher systolic BP and 1.24 ± 0.34 mmHg higher diastolic BP, as well as with 7% higher risk of hypertension (odds ratio: 1.07, [95% CI: 1.04–1.10]). In addition, per standard deviation higher of genetically predicted glycine level was associated with lower systolic BP (− 0.70 ± 0.17 mmHg, P = 4.04 × 10−5) and a lower risk of hypertension (0.99 [0.98–0.99], P = 6.46 × 10−5). In the reverse direction, genetically predicted higher systolic BP was associated with lower circulating levels of glycine (− 0.025±0.008, P = 0.001).
This study provides evidence for causal impacts of genetically predicted circulating BCAAs and glycine levels on BP. Meanwhile, genetically predicted higher BP was associated with lower glycine levels. Further investigations are warranted to clarify the underlying mechanisms.
Cardiometabolic diseases remain the leading cause of morbidity and mortality globally, and high blood pressure (BP) is not only an endpoint but also a significant risk factor [1,2,3]. Although the application of biomarkers has augmented traditional cardiometabolic prediction and helped prevention, the aim of biomarker-guided prevention is still to be fully realized . With the maturation of comprehensive metabolomic profiling technologies in the last decade, the pace of discoveries of disease-related biomarkers has been accelerated .
Identification of novel predictive biomarkers may help facilitate the detection and management of hypertension. To date, abundant evidence indicates that amino acids play a crucial role in cardiovascular physiology and pathology. Observational findings have suggested associations of circulating branched-chain amino acids (BCAAs) with cardiometabolic diseases and related risk factors [6,7,8,9,10], including hypertension [11, 12]. Increased dietary intake of aromatic amino acids, especially phenylalanine obtained from animal protein source, was associated with higher BP . Other amino acids such as glycine, glutamine, histidine, and alanine also have been associated with hypertension [14,15,16,17,18,19]. Taken together, these pieces of observational evidence (Additional file 1: Table S1) have raised the hypotheses that these amino acids may causally contribute to the development of hypertension, and the other way around, higher BP levels may also cause a change in the circulating levels of amino acids. However, population-based intervention studies are lacking due to challenges such as cost and safety issues.
Mendelian randomization (MR) analysis has emerged as a valuable approach to estimate causal associations between exposures and outcomes . In brief, MR uses single nucleotide polymorphisms (SNPs) from genome-wide association studies (GWASs) as proxies for exposure, namely instrumental variables (IVs), to make causal inferences. Based on the truth that the random allocation of parental alleles at conception, MR design is analogous to a randomized controlled trial and could therefore address conventional biases in observational studies, such as confounding and reverse causation. To date, the MR approach has contributed to several major inferences and validations of causal relationships between traditional risk factors and disease .
In the present study, we conducted two-sample MR analyses to infer the causality of circulating levels of amino acids with BP and risk of hypertension using individual-level data of European participants in UK Biobank and additional summary-level data of metabolomics GWAS by Kettunen et al. .
Study population and design
The overview of the study is presented in Fig. 1. We used the individual-level data from UK Biobank, a large prospective cohort study of about 500,000 participants across the United Kingdom . The detailed procedures of genotyping, imputation, and quality control of the genetic information in the UK Biobank have been described elsewhere . Plasma metabolic biomarkers were measured in about 120,000 randomly selected participants using a high-throughput nuclear magnetic resonance (NMR)-based metabolite profiling platform developed by Nightingale Health Ltd (Helsinki, Finland). To date, the first release of metabolic data covers 249 metabolites, including nine amino acids. The unique data identifiers within the UK Biobank repository were shown in Additional file 1: Table S2.
In this study, we focused on the circulating amino acids reported to be associated with BP in observational studies [9,10,11,12, 14,15,16,17,18,19], including nine available in UK Biobank: leucine, isoleucine, valine, alanine, phenylalanine, tyrosine, histidine, glutamine, and glycine (Additional file 1: Table S2). Following central quality control procedures of UK Biobank , we excluded (1) individuals who had withdrawn consent, (2) non-Caucasian participants, (3) those without genotyping data, (4) those with sex chromosome aneuploidy, (5) those with a mismatch between genetically inferred sex and self-reported sex, and (6) outliers for heterozygosity or missing rate (Fig. 1) in the following analyses. We included two non-overlapping populations to conduct the two-sample MR analysis. One population included the individual-level data of 98,317 UK Biobank participants who had available metabolomics measurements together with another published summary-level data. The other population consisted of the remaining 309,827 UK Biobank participants without metabolomics data (Fig. 1). The protocol of the present study was approved by UK Biobank (#54294). This study followed the Strengthening the reporting of observational studies in epidemiology using Mendelian randomization (STROBE-MR, Additional file 2: STROBE-MR checklist) guideline .
Data sources for exposures
We conducted GWAS analysis for circulating levels of amino acids in the UK Biobank samples of European ancestry with available metabolomic data (n = 98,317). We applied a rank-based inverse normal transformation to the circulating levels of amino acids. We used the Bayesian mixed-model association method [25, 26] allowing related individuals to perform GWAS in BOLT-LMM (2.3) software, with the adjustment of age, sex, fasting time, genotyping chips, assessment centers, and the top ten principal components.
For each amino acid, the results of GWAS in the UK Biobank were then meta-analyzed with the summary-level results from a published NMR-based metabolomics GWAS [21, 27] using METAL software . Briefly, this external GWAS included 24,925 European participants from 14 cohorts whose blood metabolites were profiled by using quantitative NMR metabolomics platform. The covariates included age, sex, and fasting time . Meta-analyses of these GWASs were conducted with genomic control correction under a fixed-effect model and the weights were based on standard error. We finally reached a sample size of 123,242 participants and restricted the current analyses to the 8,247,959 SNPs that were captured in both UK Biobank and this GWAS with a minor allele frequency > 1%.
To construct a valid IV, three basic assumptions must be satisfied: (1) the IV is associated with the exposure, (2) the IV is independent of any confounder in the association between the exposure and the outcome, and (3) the IV influences the outcome only through the exposure pathway. To meet assumption 1, we selected the corresponding SNPs for each circulating amino acid at the threshold of genome-wide significance (P < 5 × 10−8). Linkage disequilibrium was estimated between SNPs to select independent genetic variants using clump parameter in PLINK (1.9) software  (window size = 10,000 kb, r2 < 0.01 using 1000G Phase3 v5 EUR [GRCh37/hg19] as a reference panel). For those variants in linkage disequilibrium, we chose the one with the lowest P-value. We calculated the F statistic for each IV to assess the strength and avoid weak instrument bias . As commonly used, the cut-off was set to F statistic > 10 . The detailed information of the IV for each amino acid is shown in Additional file 1: Table S3. We excluded the SNPs that were associated with potential confounders (namely, coronary heart diseases, type 2 diabetes, and body mass index) at genome-wide significance threshold (P < 5 × 10−8) using available large-scale GWAS summary data from FinnGen study R5 release (https://www.finngen.fi/fi)  and GIANT consortium [33, 34] (Additional file 1: Table S4), respectively. The SNPs that were directly associated with outcomes were also excluded from analyses. In order to test the influence of the removal of SNPs on the results, we constructed broad sets of IVs that used all independent SNPs without any exclusion as a sensitivity analysis (Additional file 1: Table S5). All IVs in current analyses were robust instruments with F statistic >10 (Additional file 1: Tables S3 and S5).
Notably, several genetic variants for glycine are in (or near) carbamoyl-phosphate synthase 1 gene (CPS1, Additional file 1: Table S3). The CPS1 gene exhibits broad effects on the circulating levels of multiple metabolites, such as trimethylamine N-oxide and its predecessors, which have been associated with cardiovascular events, including hypertension [35,36,37,38]. Due to the known pleiotropic effects of CPS1 gene on human metabolism, we further generated one additional IV excluding variants in CPS1 gene for glycine.
Data sources for outcomes
Effect estimates for corresponding SNPs on systolic blood pressure (SBP), diastolic blood pressure (DBP), and hypertension were calculated in the UK Biobank participants without metabolomic data. Sitting BPs were measured twice using the Omron digital blood pressure monitor (Hoofddorp, Netherlands) or manual sphygmomanometer in the assessment center. Two automated and two manual BP measurements were averaged to derive SBP and DBP values for most participants, and the measured value was directly used for the participants who had single BP measurement only (0.084% of participants). Among the participants who reported taking BP-lowering medication, 15 and 10 mmHg were added to the measured SBP and DBP values, respectively . We defined hypertension cases according to the usage of BP-lowering medication, the International Classification of Diseases (ICD-9 code 401 and ICD-10 code I10), or the average BP values of two measurements (average SBP>140 mmHg or average DBP > 90 mmHg). As many medications on coronary artery disease would affect the blood pressure levels, the participants who had a history of coronary artery disease or related procedures (Additional file 1: Table S2) were excluded in the association analysis of BP or hypertension, which is consistently with published GWASs of blood pressure and hypertension risk [40, 41].
The genetic estimates for BP traits and hypertension were calculated using linear mixed models with adjustment of age, sex, genotyping chips, assessment centers, and the top ten principal components in the BOLT-LMM (2.3) software .
Data harmonization was conducted to make sure that the effect alleles were the same for both exposures and outcomes before every MR analysis. Because the summary statistics were calculated from the UK Biobank’s imputed genotype data, which were consistently aligned to the forward strand, we did not exclude palindromic SNPs with intermediate allele frequencies. For the primary analysis, we used the inverse variance weighted (IVW) method with multiplicative random effects , which is a powerful MR method under valid IVs and balanced pleiotropy assumptions, to estimate each circulating amino acid's causal effect on BP or risk of hypertension. Several sensitivity analyses were performed to assess the reliability of our findings, including MR-Egger regression , the median-based methods , and the mode-based method . Specifically, under the weaker InSIDE (INstrument Strength Independent of Direct Effect) assumption, MR-Egger regression allows the IV to have horizontal pleiotropy if such pleiotropic effect is independent of the effect size of the genetic variant on the risk factor . The directional pleiotropy can be estimated by conducting a hypothesis test for the intercept . For the two median-based methods, the simple median method assumes that at least half of the variants are valid IVs, while weighted median method assumes that valid IVs provide more than half of the weight . The mode-based method assumes a plurality of genetic variants is valid. Of note, these sensitivity analyses mentioned above are robust to violations of MR assumptions in some extent, but result in a relatively lower power to detect a causal effect.
Additionally, because circulating levels of glycine are correlated with that of BCAAs in human [47, 48], we performed multivariable MR through the IVW method to estimate the independent causal effect of circulating amino acids on outcomes. Because of the strong correlations among the BCAAs (Additional file 3: Fig. S1), we included each BCAA at one time in the multivariable MR model. The conditional F statistics  for each exposure was calculated to describe the instrument strength.
Results were presented as β ± SE mmHg change in BP measurements and odds ratios (ORs) (95% confidence intervals [95% CI]) for hypertension risk per one standard deviation (SD) increase in genetically predicted circulating amino acids level. A robust causal effect was defined by a consistency of directions across the above-mentioned MR methods. We used the MR-Egger method to detect the presence of potential pleiotropy. Heterogeneity was assessed using Cochran’s Q statistics for all circulating amino acids. If heterogeneity was detected, a multiplicative random-effects IVW method would be applied . Meanwhile, a leave-one-out analysis was performed to check if any single SNP strongly drove the causal effects. To address the multiple testing issues, the threshold of statistical significance was defined as P < 0.006 (0.05/9 amino acids) after Bonferroni correction. Findings with 0.006 ≤ P < 0.05 were considered as associations of suggestively significant evidence. All the MR analyses were performed using R software (3.6.1) with the TwoSampleMR package (0.5.6) .
Reverse MR analysis
To estimate the causal effects of genetically predicted BP or risk of hypertension on circulating levels of amino acids, we further performed bi-directional MR analyses using the IVW method. Independent genetic variants associated with the two BP traits and hypertension (P < 5 × 10−8; r2 < 0.01, European panel from the 1000 Genome Project) came from the GWASs in UK Biobank samples without metabolomic data as described above. The SNPs that were directly associated with circulating levels of amino acids were excluded from analyses. We constructed IVs including 360 SNPs for SBP, 358 for DBP, and 238 for hypertension, respectively (Additional file 1: Table S6). Data harmonization was also performed as described in forward MR analysis. The results were presented as β ± SE and the estimates represent SD change in circulating levels of amino acids per 10 mmHg of genetically predicted higher BP or per unit higher in log odds of hypertension risk. Again, MR-Egger regression, the median-based and mode-based methods were applied to test the reliability of the findings in sensitivity analyses. All the MR analyses were also performed using R software (3.6.1) with the TwoSampleMR package (0.5.6) .
Summary of population characteristics
Among all included participants of UK Biobank (n = 384,707), the average age was 56.6 years, 44.4% were men, and 18.8% reported use of antihypertensive drugs (Additional file 1: Table S7). In the 286,390 participants without metabolomic data, the mean (SD) values of SBP and DBP were 138.02 (18.58) and 82.35 (10.06) mmHg, respectively, and 146 345 hypertensive participants were documented (Additional file 1: Table S7 and Additional file 3: Fig. S2). The distributions of circulating levels of amino acids are shown in Additional file 1: Table S7 and Additional file 3: Fig. S3.
Association of genetically predicted circulating levels of amino acids with BP and risk of hypertension
Several potential causal effects of circulating amino acids on BP traits were observed (Fig. 2 and Additional file 1: Table S8). Genetically predicted higher levels of circulating BCAAs, leucine, isoleucine, and valine, were significantly associated with higher SBP (2.22 ± 0.70 mmHg per 1-SD increase in leucine; 2.71 ± 0.78 mmHg per 1-SD increase in isoleucine; 1.96 ± 0.55 mmHg per 1-SD increase in valine, all P < 0.006). Meanwhile, genetically predicted higher levels of isoleucine and valine were significantly associated with higher DBP (1.24 ± 0.34 mmHg per 1-SD increase in isoleucine; 0.94 ± 0.32 mmHg DBP per 1-SD increase in valine, all P < 0.006). We also observed suggestive associations of genetically predicted higher circulating leucine with higher DBP (1.01 ± 0.38 mmHg per 1-SD, P = 0.008). On the contrary, genetically predicted higher circulating glycine was associated with lower SBP (− 0.70 ± 0.17 mmHg per 1-SD, P = 4.04 × 10−5) and was suggestively associated with lower DBP (− 0.23 ± 0.10 mmHg per 1-SD, P = 0.02). Meanwhile, the causal effects of circulating levels of amino acids on the risk of hypertension were in accordance with the results from BP measurements. Genetically predicted higher circulating levels of three BCAAs were associated with a higher risk of hypertension (OR 1.05 [95% CI 1.03–1.08] per 1-SD increase in leucine, 1.07 [1.04-1.10] per 1-SD increase in isoleucine, 1.05 [1.02-1.07] per 1-SD increase in valine, all P < 0.006). Likewise, genetically predicted higher level of circulating glycine was associated with a lower risk of hypertension (0.99 [0.98–0.99], P < 0.006). Analyses using the broad sets of IVs showed consistent results, providing robust evidence on the causal effects of genetically determined circulating levels of BCAAs and glycine on BP (Additional file 1: Table S9 and Additional file 3: Fig. S4).
In addition, the association directions in sensitivity analyses using different methods were consistent with the results from the IVW method, and no evidence of pleiotropy was detected in the MR-Egger regression test (P > 0.05 for intercepts, Additional file 1: Table S8 and Additional file 3: Fig. S5). The leave-one-out analyses suggested that the associations were robust to exclusion of any single genetic variant (Additional file 3: Figs. S6-S8), though the Cochran’s Q test reported the presence of significant heterogeneity (Additional file 1: Table S8).
Because circulating levels of glycine are correlated with that of BCAAs in human [47, 48], we next conducted multivariable MR analyses, including glycine with one BCAA at a time as exposure to estimate their independent causal effects on BP traits and hypertension. After adjusting for isoleucine, the BCAA with the largest effect size, genetically predicted higher glycine was still significantly associated with a lower level of SBP, whereas the association with DBP was attenuated to null. Meanwhile, the positive association of the circulating level of isoleucine with BP traits and hypertension persisted (Additional file 1: Table S10). The results were similar when glycine was additionally adjusted for leucine or valine (Additional file 1: Table S10).
Association of genetically predicted circulating glycine with BP and hypertension with respect to CPS1 gene
In consideration of the remarkable pleiotropy of CPS1 variants  and their strong effects on glycine (Additional file 1: Table S3), we investigated the associations of an additional IV for glycine without CPS1 variants with outcomes. When using IV without the pleiotropic CPS1 locus, only a suggestive association between circulating level of glycine and SBP (− 0.79 ± 0.30 mmHg per 1-SD, P = 0.009) or risk of hypertension (0.98 [0.97–0.99] per 1-SD, P = 0.009) were observed, and there was no causal effect of circulating glycine on DBP (P > 0.05, Fig. 3). The results were similar in sensitivity analyses, and the MR-Egger regression suggested that the results were not influenced by horizontal pleiotropy (all P > 0.05 for intercepts, Additional file 1: Table S8). Additionally, the combination of CPS1 variants for glycine exhibited significant associations with lower BP levels and risk of hypertension (Additional file 1: Table S8).
Reverse MR analysis assessing the causal effect of BP on circulating amino acids
We found that genetically predicted 10 mmHg higher in SBP was associated with a lower circulating level of glycine (− 0.025 ± 0.008, P = 0.001, Table 1). Genetically predicted higher BP or hypertension risk were not associated with other circulating levels of amino acids after Bonferroni correction (all P>0.006, Table 1 and Additional file 1: Table S11), though some suggestive associations were presented. For example, at suggestive significance level, genetically predicted higher SBP was positively associated with tyrosine, genetically predicted higher DBP was positively associated with isoleucine and valine but inversely associated with glycine, and genetically predicted higher risk of hypertension was associated with higher levels of circulating isoleucine and tyrosine, but lower levels of glycine (0.006 ≤ P < 0.05, Table 1 and Additional file 1: Table S11). In addition, although genetically predicted higher SBP was suggestively associated with histidine in the IVW method, the MR-Egger method found potential pleiotropy in the analysis with null causal effect (Additional file 1: Table S11). Besides, the directions were consistent in sensitivity analyses, and no other evidence of pleiotropy was detected in the MR-Egger regression test (Additional file 1: Table S11 and Additional file 3: Fig. S9). The Cochran’s Q test identified significant heterogeneity in the reverse MR analyses (Additional file 1: Table S11).
Using the two-sample MR approach, this study showed that genetically predicted higher circulating levels of BCAAs were associated with higher levels of BP and increased risk of hypertension. In contrast, genetically predicted higher circulating level of glycine was consistently associated with lower SBP and DBP, supporting a potential protective role in hypertension. Although the effect size was attenuated, genetically predicted glycine has a robust BP-lowering effect in the sensitivity analysis when pleiotropic CPS1 variant were excluded. In the reverse MR analysis, we found that genetically predicted higher SBP was associated with a lower circulating level of glycine. Heterogeneity was detected among the instrumental variables. Though random-effects IVW method was applied and the results from the sensitivity analysis indicated the general robustness of our results, the biological interpretations should be made with caution.
The results from the present MR analysis are in line with previous observational studies, which suggested that high levels of circulating BCAAs were associated with higher BP measurements , or with an increased risk of incident hypertension . Intriguingly, circulating levels of BCAAs are highly correlated and may have common genetic background. In the MR analyses, four common SNPs were identified as the shared IVs for leucine, isoleucine, and valine. The most significant variant, rs10014755 in PPM1K gene, encodes the mitochondrial phosphatase that regulate a rate-limiting enzyme for BCAA catabolism [8, 53], and two other variants, rs117643180 and rs2422358 in solute-carrier gene, affect amino acid transport . These shared genetic mechanisms influence the transport and catabolism of BCAAs, thereby impact the circulating levels of BCAAs simultaneously.
One possible mechanism linking circulating BCAAs to BP involves the persistent activation of mTORC1 pathway in vascular system by BCAAs . The activation of mTORC1 leads to changes of sympathetic nerve traffic , proliferation of vascular smooth muscle cell , kidney injury , and eventually the development of hypertension . Another possible explanation is that higher levels of circulating BCAAs promote inflammation and oxidative stress in endothelial cells and elicit inflammatory cell adhesion and endothelial dysfunction, which contribute to increased risk of hypertension [59, 60]. Despite of these existing speculations, the investigation of the detailed mechanisms is still ongoing and further studies are warranted.
Furthermore, glycine is a non-essential amino acid that has been inversely associated with the risk of cardiovascular diseases  and type 2 diabetes [62, 63]. The potential protective role of glycine in hypertension prevention was supported by both cross-sectional studies  and prospective studies . Of note, our results corroborated with previous observational findings and suggested a causal effect of genetically predicted higher circulating glycine on lower levels of BP or risk of hypertension. However, the original IV for glycine in the present study was largely driven by the CPS1 locus, which was not only restricted to the glycine pathway but also exhibited significant effects on plasma trimethylamine N-oxide, betaine, and other intermediate metabolites in urea cycle . A recent MR analysis, which aimed to evaluate the causal association between circulating glycine and coronary artery disease, also highlighted the potential pleiotropy of the CPS1 rs1047891 SNP and did not obtain conclusive evidence for causality . In addition, suggestive effects of glycine on SBP and risk of hypertension were still observed even when all CPS1 variants were removed from the IV, and such findings were in line with another MR study based on six variants for glycine on cardiometabolic diseases when excluding the CPS1 variant . These consistent associations in analysis using IVs with and without CPS1 variants suggested a robust and direct role of circulating glycine in BP management, independent of potential traditional risk factors.
While emerging observational studies implied the effects of circulating amino acids on hypertension risk , few studies explored whether hypertension would affect the circulating amino acids levels. In our reverse MR analyses, genetically predicted BP was associated with circulating levels of glycine and suggestively associated with circulating levels of isoleucine, tyrosine, and histidine. One possible explanation linking BP and glycine is that hypertension may reduce peripheral glucose utilization and promote insulin resistance and gluconeogenesis  and consequently accelerate the conversion of glycine to serine . However, to our knowledge, there is no obvious mechanism well explaining the observed suggestive effects of higher BP on increased levels of BCAAs and aromatic amino acids. Further studies are needed to validate our observations and provide clues to the underlying mechanisms.
Our MR study supports a causal role of BCAAs in the pathophysiology of hypertension, and proposes a potential BP-lowering effect of glycine, implying the potential of clinical measurements and monitoring of circulating amino acids. This MR study has several strengths. First, to the best of our knowledge, this is the first study that has comprehensively scanned the causal effects of circulating amino acids on BP or hypertension using an MR approach. The two-sample MR design can minimize the bias inherent to observational study, such as residual confounding and reverse causality. Second, to calculate the genetic effect estimates for circulating amino acids, we performed meta-analyses using individual-level data in UK Biobank and summary statistics from another external GWAS, resulting in the largest sample size for the metabolomic GWAS to date. By splitting UK Biobank population, the current combination of data sources can avoid bias due to sample overlap and maintain a relatively large sample size with high SNP consistency. Meanwhile, the selection of SNPs was restricted by several procedures to meet basic assumptions for the MR approach and to ensure the validity and strength of IVs. Third, we conducted a thorough evaluation considering the pleiotropic effects of the CPS1 variants, and our results further emphasized the robust roles of glycine in BP regulation. Finally, several sensitivity analyses were conducted to ensure consistency and robustness for our findings.
Of note, there are limitations in our study. First, due to the widespread pleiotropic effects of metabolite-related SNPs, it is hard to rule out the possibility that selected genetic variants may affect BP outcomes through an alternative causal pathway than through amino acid exposure. Nevertheless, we used several MR approaches as sensitivity analyses to ensure the robustness of our results against violation of MR assumptions. Another limitation is that the effect sizes of these amino acids on BP are relatively small. Finally, the present MR study relied on genetic estimates that derived from GWAS in participants of European descent, which though reduced population stratification, might limit the generalization of our findings to non-European populations.
In summary, our findings support the causal association of higher circulating BCAAs levels and lower glycine levels with higher BP and elevated risk of hypertension. Furthermore, we are the first to report the associations of genetically predicted hypertension with circulating amino acids, suggesting an impact of hypertension on systematic metabolism. Our findings may improve the understanding of the connections between circulating levels of amino acids and BP. However, further studies are needed to clarify the potential role of amino acids in the pathophysiology of hypertension and identify plausible biological mechanisms.
Availability of data and materials
The UK Biobank data are available from the UK Biobank on application (www.ukbiobank.ac.uk/). All the supporting data for MR analyses are available within the article and supporting information. Further data can be obtained by a reasonable request to the corresponding author.
Branched-chain amino acid
- CPS1 :
Carbamoyl-phosphate synthase 1
Diastolic blood pressure
Genome-wide association study
International Classification of Diseases
Inverse variance weighted
Nuclear magnetic resonance
Systolic blood pressure
Single nucleotide polymorphism
STrengthening the Reporting of Observational Studies in Epidemiology using Mendelian Randomisation
Fuchs FD, Whelton PK. High blood pressure and cardiovascular disease. Hypertension. 2020;75(2):285–92.
Sattar N, Gill JMR, Alazawi W. Improving prevention strategies for cardiometabolic disease. Nat Med. 2020;26(3):320–5.
Flint AC, Conell C, Ren X, Banki NM, Chan SL, Rao VA, et al. Effect of systolic and diastolic blood pressure on cardiovascular outcomes. N Engl J Med. 2019;381(3):243–51.
Roberts LD, Gerszten RE. Toward new biomarkers of cardiometabolic diseases. Cell Metab. 2013;18(1):43–50.
Johnson CH, Ivanisevic J, Siuzdak G. Metabolomics: beyond biomarkers and towards mechanisms. Nat Rev Mol Cell Biol. 2016;17(7):451–9.
Lynch CJ, Adams SH. Branched-chain amino acids in metabolic signalling and insulin resistance. Nat Rev Endocrinol. 2014;10(12):723–36.
Guasch-Ferre M, Hruby A, Toledo E, Clish CB, Martinez-Gonzalez MA, Salas-Salvado J, et al. Metabolomics in prediabetes and diabetes: a systematic review and meta-analysis. Diabetes Care. 2016;39(5):833–46.
Neinast M, Murashige D, Arany Z. Branched chain amino acids. Annu Rev Physiol. 2019;81:139–64.
Poggiogalle E, Fontana M, Giusti AM, Pinto A, Iannucci G, Lenzi A, et al. Amino acids and hypertension in adults. Nutrients. 2019;11(7):1459.
Chakraborty S, Mandal J, Yang T, Cheng X, Yeo JY, McCarthy CG, et al. Metabolites and hypertension: insights into hypertension as a metabolic disorder: 2019 Harriet Dustan Award. Hypertension. 2020;75(6):1386–96.
Flores-Guerrero JL, Groothof D, Connelly MA, Otvos JD, Bakker SJL, Dullaart RPF. Concentration of branched-chain amino acids is a strong risk marker for incident hypertension. Hypertension. 2019;74(6):1428–35.
Mahbub MH, Yamaguchi N, Hase R, Takahashi H, Ishimaru Y, Watanabe R, et al. Plasma branched-chain and aromatic amino acids in relation to hypertension. Nutrients. 2020;12(12):3791. https://doi.org/10.3390/nu12123791.
Teymoori F, Asghari G, Mirmiran P, Azizi F. High dietary intake of aromatic amino acids increases risk of hypertension. J Am Soc Hypertens. 2018;12(1):25–33.
Dietrich S, Floegel A, Weikert C, Prehn C, Adamski J, Pischon T, et al. Identification of serum metabolites associated with incident hypertension in the European Prospective Investigation into Cancer and Nutrition-Potsdam Study. Hypertension. 2016;68(2):471–7.
Wang L, Hou E, Wang L, Wang Y, Yang L, Zheng X, et al. Reconstruction and analysis of correlation networks based on GC-MS metabolomics data for young hypertensive men. Anal Chim Acta. 2015;854:95–105.
Ntzouvani A, Nomikos T, Panagiotakos D, Fragopoulou E, Pitsavos C, McCann A, et al. Amino acid profile and metabolic syndrome in a male Mediterranean population: a cross-sectional study. Nutr Metab Cardiovasc Dis. 2017;27(11):1021–30.
Holmes E, Loo RL, Stamler J, Bictash M, Yap IK, Chan Q, et al. Human metabolic phenotype diversity and its association with diet and blood pressure. Nature. 2008;453(7193):396–400.
Liu Y, Chen T, Qiu Y, Cheng Y, Cao Y, Zhao A, et al. An ultrasonication-assisted extraction and derivatization protocol for GC/TOFMS-based metabolite profiling. Anal Bioanal Chem. 2011;400(5):1405–17.
Zhong L, Zhang JP, Nuermaimaiti AG, Yunusi KX. Study on plasmatic metabolomics of Uygur patients with essential hypertension based on nuclear magnetic resonance technique. Eur Rev Med Pharmacol Sci. 2014;18(23):3673–80.
Holmes MV, Ala-Korpela M, Smith GD. Mendelian randomization in cardiometabolic disease: challenges in evaluating causality. Nat Rev Cardiol. 2017;14(10):577–90.
Kettunen J, Demirkan A, Wurtz P, Draisma HH, Haller T, Rawal R, et al. Genome-wide study for circulating metabolites identifies 62 loci and reveals novel systemic effects of LPA. Nat Commun. 2016;7:11122.
Sudlow C, Gallacher J, Allen N, Beral V, Burton P, Danesh J, et al. UK Biobank: an open access resource for identifying the causes of a wide range of complex diseases of middle and old age. PLoS Med. 2015;12(3):e1001779.
Bycroft C, Freeman C, Petkova D, Band G, Elliott LT, Sharp K, et al. The UK Biobank resource with deep phenotyping and genomic data. Nature. 2018;562(7726):203–9.
Skrivankova VW, Richmond RC, Woolf BAR, Davies NM, Swanson SA, VanderWeele TJ, et al. Strengthening the reporting of observational studies in epidemiology using mendelian randomisation (STROBE-MR): explanation and elaboration. BMJ. 2021;375:n2233.
Loh PR, Kichaev G, Gazal S, Schoech AP, Price AL. Mixed-model association for biobank-scale datasets. Nat Genet. 2018;50(7):906–8.
Loh PR, Tucker G, Bulik-Sullivan BK, Vilhjalmsson BJ, Finucane HK, Salem RM, et al. Efficient Bayesian mixed-model analysis increases association power in large cohorts. Nat Genet. 2015;47(3):284–90.
Kettunen J, Demirkan A, Wurtz P, Draisma HH, Haller T, Rawal R, Vaarhorst A, Kangas AJ, Lyytikainen LP, Pirinen M et al. MAGNETIC NMR-GWAS summary statistics. http://www.computationalmedicine.fi/data. Accessed 10 August 2020.
Willer CJ, Li Y, Abecasis GR. METAL: fast and efficient meta-analysis of genomewide association scans. Bioinformatics. 2010;26(17):2190–1.
Chang CC, Chow CC, Tellier LC, Vattikuti S, Purcell SM, Lee JJ. Second-generation PLINK: rising to the challenge of larger and richer datasets. Gigascience. 2015;4:7.
Burgess S, Thompson SG, Collaboration CCG. Avoiding bias from weak instruments in Mendelian randomization studies. Int J Epidemiol. 2011;40(3):755–64.
Staiger D, Stock JH. Instrumental variables regression with weak instruments. Econometrica. 1997;65(3):557–86.
FinnGen consortium. FinnGen GWAS summary statistics, R5. https://www.finngen.fi/en/access_results. Accessed 1 October 2021
Yengo L, Sidorenko J, Kemper KE, Zheng Z, Wood AR, Weedon MN, et al. Meta-analysis of genome-wide association studies for height and body mass index in approximately 700000 individuals of European ancestry. Hum Mol Genet. 2018;27(20):3641–9.
Yengo L, Sidorenko J, Kemper KE, Zheng Z, Wood AR, Weedon MN, Frayling TM, Hirschhorn J, Yang J, Visscher PM et al. 2018 GIANT and UK BioBank Meta-analysis. https://portals.broadinstitute.org/collaboration/giant/index.php/GIANT_consortium_data_files. Accessed 1 October 2021.
Meijer AJ, Lamers WH, Chamuleau RA. Nitrogen metabolism and ornithine cycle function. Physiol Rev. 1990;70(3):701–48.
Hartiala JA, Tang WH, Wang Z, Crow AL, Stewart AF, Roberts R, et al. Genome-wide association study and targeted metabolomics identifies sex-specific association of CPS1 with coronary artery disease. Nat Commun. 2016;7:10558.
Senthong V, Li XS, Hudec T, Coughlin J, Wu Y, Levison B, et al. Plasma Trimethylamine N-oxide, a gut microbe-generated phosphatidylcholine metabolite, is associated with atherosclerotic burden. J Am Coll Cardiol. 2016;67(22):2620–8.
Ge X, Zheng L, Zhuang R, Yu P, Xu Z, Liu G, et al. The gut microbial metabolite trimethylamine N-oxide and hypertension risk: a systematic review and dose-response meta-analysis. Adv Nutr. 2020;11(1):66–76.
Evangelou E, Warren HR, Mosen-Ansorena D, Mifsud B, Pazoki R, Gao H, et al. Genetic analysis of over 1 million people identifies 535 new loci associated with blood pressure traits. Nat Genet. 2018;50(10):1412–25.
International Consortium for Blood Pressure Genome-Wide Association Studies, Ehret GB, Munroe PB, Rice KM, Bochud M, Johnson AD, et al. Genetic variants in novel pathways influence blood pressure and cardiovascular disease risk. Nature. 2011;478(7367):103–9.
Newton-Cheh C, Johnson T, Gateva V, Tobin MD, Bochud M, Coin L, et al. Genome-wide association study identifies eight loci associated with blood pressure. Nat Genet. 2009;41(6):666–76.
Bowden J, Del Greco MF, Minelli C, Davey Smith G, Sheehan N, Thompson J. A framework for the investigation of pleiotropy in two-sample summary data Mendelian randomization. Stat Med. 2017;36(11):1783–802.
Bowden J, Davey Smith G, Burgess S. Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int J Epidemiol. 2015;44(2):512–25.
Bowden J, Davey Smith G, Haycock PC, Burgess S. Consistent estimation in Mendelian randomization with some invalid instruments using a weighted median estimator. Genet Epidemiol. 2016;40(4):304–14.
Hartwig FP, Davey Smith G, Bowden J. Robust inference in summary data Mendelian randomization via the zero modal pleiotropy assumption. Int J Epidemiol. 2017;46(6):1985–98.
Burgess S, Thompson SG. Interpreting findings from Mendelian randomization using the MR-Egger method. Eur J Epidemiol. 2017;32(5):377–89.
Glynn EL, Piner LW, Huffman KM, Slentz CA, Elliot-Penry L, AbouAssi H, et al. Impact of combined resistance and aerobic exercise training on branched-chain amino acid turnover, glycine metabolism and insulin sensitivity in overweight humans. Diabetologia. 2015;58(10):2324–35.
White PJ, Lapworth AL, McGarrah RW, Kwee LC, Crown SB, Ilkayeva O, et al. Muscle-liver trafficking of BCAA-derived nitrogen underlies obesity-related glycine depletion. Cell Rep. 2020;33(6):108375.
Sanderson E, Davey Smith G, Windmeijer F, Bowden J. An examination of multivariable Mendelian randomization in the single-sample and two-sample summary data settings. Int J Epidemiol. 2019;48(3):713–27.
Bowden J, Holmes MV. Meta-analysis and Mendelian randomization: a review. Res Synth Methods. 2019;10(4):486–96.
Hemani G, Zheng J, Elsworth B, Wade KH, Haberland V, Baird D, et al. The MR-Base platform supports systematic causal inference across the human phenome. Elife. 2018;7:e34408. https://doi.org/10.7554/eLife.34408.
Siomkajlo M, Rybka J, Mierzchala-Pasierb M, Gamian A, Stankiewicz-Olczyk J, Bolanowski M, et al. Specific plasma amino acid disturbances associated with metabolic syndrome. Endocrine. 2017;58(3):553–62.
Oyarzabal A, Martinez-Pardo M, Merinero B, Navarrete R, Desviat LR, Ugarte M, et al. A novel regulatory defect in the branched-chain alpha-keto acid dehydrogenase complex due to a mutation in the PPM1K gene causes a mild variant phenotype of maple syrup urine disease. Hum Mutat. 2013;34(2):355–62.
Yahyaoui R, Perez-Frias J. Amino acid transport defects in human inherited metabolic disorders. Int J Mol Sci. 2019;21(1):119. https://doi.org/10.3390/ijms21010119.
Harlan SM, Guo DF, Morgan DA, Fernandes-Santos C, Rahmouni K. Hypothalamic mTORC1 signaling controls sympathetic nerve activity and arterial pressure and mediates leptin effects. Cell Metab. 2013;17(4):599–606.
Lightell DJ Jr, Moss SC, Woods TC. Loss of canonical insulin signaling accelerates vascular smooth muscle cell proliferation and migration through changes in p27Kip1 regulation. Endocrinology. 2011;152(2):651–8.
Kumar V, Wollner C, Kurth T, Bukowy JD, Cowley AW Jr. Inhibition of mammalian target of rapamycin complex 1 attenuates salt-induced hypertension and kidney injury in Dahl salt-sensitive rats. Hypertension. 2017;70(4):813–21.
Reho JJ, Guo DF, Morgan DA, Rahmouni K. mTORC1 (mechanistic target of rapamycin complex 1) signaling in endothelial and smooth muscle cells is required for vascular function. Hypertension. 2021;77(2):594–604.
Zhenyukh O, Civantos E, Ruiz-Ortega M, Sanchez MS, Vazquez C, Peiro C, et al. High concentration of branched-chain amino acids promotes oxidative stress, inflammation and migration of human peripheral blood mononuclear cells via mTORC1 activation. Free Radic Biol Med. 2017;104:165–77.
Zhenyukh O, Gonzalez-Amor M, Rodrigues-Diez RR, Esteban V, Ruiz-Ortega M, Salaices M, et al. Branched-chain amino acids promote endothelial dysfunction through increased reactive oxygen species generation and inflammation. J Cell Mol Med. 2018;22(10):4948–62.
Ding Y, Svingen GF, Pedersen ER, Gregory JF, Ueland PM, Tell GS, et al. Plasma glycine and risk of acute myocardial infarction in patients with suspected stable angina pectoris. J Am Heart Assoc. 2015;5(1):e002621. https://doi.org/10.1161/JAHA.115.002621.
Floegel A, Stefan N, Yu Z, Muhlenbruch K, Drogan D, Joost HG, et al. Identification of serum metabolites associated with risk of type 2 diabetes using a targeted metabolomic approach. Diabetes. 2013;62(2):639–48.
Merino J, Leong A, Liu CT, Porneala B, Walford GA, von Grotthuss M, et al. Metabolomics insights into early type 2 diabetes pathogenesis and detection in individuals with normal fasting glucose. Diabetologia. 2018;61(6):1315–24.
Jia Q, Han Y, Huang P, Woodward NC, Gukasyan J, Kettunen J, et al. Genetic determinants of circulating glycine levels and risk of coronary artery disease. J Am Heart Assoc. 2019;8(10):e011922.
Wittemans LBL, Lotta LA, Oliver-Williams C, Stewart ID, Surendran P, Karthikeyan S, et al. Assessing the causal association of glycine with risk of cardio-metabolic diseases. Nat Commun. 2019;10(1):1060.
Polak-Iwaniuk A, Harasim-Symbor E, Golaszewska K, Chabowski A. How hypertension affects heart metabolism. Front Physiol. 2019;10:435.
Dasarathy S, Kasumov T, Edmison JM, Gruca LL, Bennett C, Duenas C, et al. Glycine and urea kinetics in nonalcoholic steatohepatitis in human: effect of intralipid infusion. Am J Physiol Gastrointest Liver Physiol. 2009;297(3):G567–75.
UK Biobank approved this present research under application number 54294. We gratefully acknowledge the contributors to the data used in this work.
This research was supported by the National Key R&D Program of China (2020YFC2005000) and the Program for Professor of Special Appointment (Eastern Scholar) at Shanghai Institutions of Higher Learning. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Ethics approval and consent to participate
The UK Biobank study was approved by the North West Multi-Centre Research Ethics Committee. Written consents were obtained from all participants in the UK Biobank.
Consent for publication
The authors declared that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Overview of observational studies investigating associations of amino acids with blood pressure or hypertension. Table S2. Definitions and the unique data identifiers within UK Biobank data. Table S3. Characteristics of the instrumental variables for circulating levels of amino acids in primary analyses. Table S4. Association of selected SNPs in our primary analyses with potential confounders. Table S5. Characteristics of the broad sets of instrumental variables for circulating levels of amino acids in sensitivity analysis. Table S6. Characteristics of instrumental variables for SBP, DBP and hypertension in reverse MR analyses. Table S7. Characteristics of participants from UK Biobank included in the MR analyses. Table S8. Causal effects of circulating levels of amino acids on BP and hypertension using different methods. Table S9. Causal effects of circulating levels of amino acids on BP and hypertension using the broad sets of instrumental variables. Table S10. MVMR estimates the independent effects of BCAA and glycine on SBP, DBP and risk of hypertension. Table S11. Causal effects of BP and hypertension on circulating levels of amino acid using different methods.
Heatmap of Pearson correlation coefficients between circulating levels of amino acids in UK Biobank. Figure S2. Distribution of systolic blood pressure (SBP) and diastolic blood pressure (DBP) in 286 390 European participants without metabolomic data from UK Biobank. Figure S3. Distribution of circulating levels of amino acids in 98 317 European participants with metabolomic data from UK Biobank. Figure S4. Estimated causal effects of circulating levels of amino acids on blood pressure and risk of hypertension. Figure S5. Scatter plots of SNPs used as instrumental variables for the MR analyses of circulating amino acids with (A) systolic blood pressure, (B) diastolic blood pressure and (C)hypertension. Figure S6. Leave-one-out plots to assess if a single SNP is driving the causal effects of circulating amino acids on systolic blood pressure (SBP). Figure S7. Leave-one-out plots to assess if a single SNP is driving the causal effects of circulating amino acids on diastolic blood pressure (DBP). Figure S8. Leave-one-out plots to assess if a single SNP is driving the causal effects of circulating amino acids on hypertension. Figure S9. Scatter plots of SNPs used as IVs for the reverse MR analyses of (A) systolic blood pressure, (B) diastolic blood pressure and (C) hypertension with circulating levels of amino acids.
About this article
Cite this article
Lin, C., Sun, Z., Mei, Z. et al. The causal associations of circulating amino acids with blood pressure: a Mendelian randomization study. BMC Med 20, 414 (2022). https://doi.org/10.1186/s12916-022-02612-w