Carotid–femoral pulse-wave velocity (PWV), a measure of large artery stiffness, is a well known independent predictor of cardiovascular morbidity and mortality [1–4], and it is considered an integrative measure of the impact of cardiovascular risk factors . Though arterial stiffening has been associated with ageing , hypertension , diabetes mellitus , hypercholesterolaemia [9,10], and chronic kidney diseases , its physiopathology is still not fully understood.
Recent advances in metabolomics have allowed for high-throughput assay of an extensive set of small molecules in a number of biological fluids. Low-molecular weight metabolites represent the intermediates and end-products of metabolic pathways that reflect physiological functions and, thus, may mirror the early stages of a pathological state .
A small study followed longitudinally 174 individuals and found that plasma fatty acid composition (assessing levels of 10 fatty acids) correlates with both PWV and mortality . Full metabolomic profiling regarding PWV has not yet been attempted.
In this study, we performed metabolomic screening in a large cohort of women from TwinsUK to identify the novel metabolites that associate with PWV. We further explored the relationship between metabolites associated with PWV and gene-expression data to further understand the molecular mechanisms underlying arterial stiffening. We also explored the association of the identified metabolites with the Framingham 10-year cardiovascular risk assessment .
MATERIALS AND METHODS
The study participants were twins enrolled in the TwinsUK Registry, a national register of adult twins recruited as volunteers without selecting for any particular disease or trait . All recruited twins were of the same sex. In this study, we analysed data from 1797 female twins who had complete data for body composition and metabolomics profiling.
The study was approved by the St. Thomas’ Hospital Research Ethics Committee, and all twins provided informed written consent.
Pulse-wave velocity measurement
Carotid–femoral PWV was calculated from the sequential recordings of carotid and femoral artery pressure waveforms using the same SphygmoCor device and applanation tonometry. Difference in the time of pulse arrival from the R-wave of the electrocardiogram between the two sites was taken as the transit time, and the difference in path length was estimated using surface measurements with a flexible tape measure between the sternal notch and the point of applanation at the femoral artery as previously described. PWV is determined by dividing the distance by transit time. Coefficient of variation between operators was less than 10% . Measurements were made in triplicate, and mean values were used for analysis. PWV measurements were available in 1797 women.
Nontargeted metabolite detection and quantification was conducted by the metabolomics provider Metabolon, Inc. (Durham, North Carolina, USA) on fasting blood samples, as described previously . In this study, we analysed 280 structurally named biochemicals (known metabolites) categorized into the following broad categories – amino acids, acylcarnitines, lysophospholipids, carbohydrates, vitamins, lipids, nucleotides, peptides, xenobiotics, and steroids.
Muther expression data
The Muther gene-expression dataset consists of 825 abdominal fat samples. Gene expression was analysed with the Illumina Human HT-12 V3 (Illumina Inc., San Diego, California, USA) , 586 individuals entered the metabolite association analysis.
Statistical analysis was carried out using Stata version 11 (Stata Corp., College Station, Texas, USA). We inverse normalized the metabolite data, as the metabolite concentrations were not normally distributed. To avoid spurious false-positive associations because of small sample size, we excluded metabolic traits with more than 20% missing values. We imputed the missing values using the minimum run day measures.
We looked for the metabolites associated with PWV by running random intercept linear regression adjusting for age, BMI, metabolite batch, and family relatedness. We corrected for multiple comparisons using Bonferroni correction, thus giving a significant threshold of P = 1.8 × 10−4 (0.05/280 metabolites). We then used a stepwise backward regression model to identify a set of metabolites that were significantly associated with each phenotype using P less than 0.01 as cut-off threshold. As metabolites in their nature can be affected by many factors, in particular dietary factors , we run random intercept linear regressions to test the effect of fruit and vegetable intake and alcohol intake on the association between the metabolites and PWV.
Associations of metabolites with gene-expression levels in fat were tested using random intercept linear regression after adjusting for age, BMI, metabolite batch, expression batch, and family relatedness.
Finally, we explored the association of selected metabolites with Framingham 10-year cardiovascular risk  by running random intercept logistic regression adjusting for age, BMI, metabolite batch, and family relatedness.
The descriptive characteristics of the study participants are shown in Table 1. After adjustment for covariates, we identified 12 metabolites, whose levels significantly correlated with PWV (Table 2). We then proceeded to analyse, in a multivariate model, which of these metabolites contributed independently and identified only three metabolites: phenylacetylglutamine, serine, and uridine (Table 3). The proportion of the variance explained by the circulating levels of these three compounds is R2 = 30%. Of the three metabolites, only phenylacetylglutamine was associated with both SBP and DBP [SBP: −1.6 (−2.36; −0.84), P
= 4.1 × 10−5; DBP: −0.88 (−1.33; −0.44), P = 8.8 × 10−5]; however, the metabolite–PWV association remained significant even after adjusting for mean arterial pressure (MAP) in the linear model and after adjusting for dietary factors (fruit and vegetable intake and alcohol intake). Published studies have shown that PWV can predict cardiovascular risk that is not accounted by the traditional factors included in the Framingham risk score . We therefore proceeded to assess whether these three metabolites were associated with the Framingham risk score. We find that phenylacetylglutamine is associated with the Framingham cardiovascular risk score [Beta = −0.04, standard error (SE) = 0.01, P = 0.004]. However, there was no association with both uridine and serine and the Framingham cardiovascular risk score (uridine: Beta = −0.003, SE = 0.01, P = 0.8; serine: Beta = −0.002, SE = 0.01, P = 0.85). Also, none of the three metabolites were associated with either total or HDL cholesterol. This suggests that some of the molecular pathways contributing to PWV are independent of the traditional cardiovascular disease risk factors measured by the Framingham score (Fig. 1).
Uridine is associated with lower PWV in our data (Table 2). When we tested for correlations between uridine and gene expression, we found that none of the expression probes passed Bonferroni correction for uridine. However, circulating uridine levels were nominally associated with the gene-expression levels in fat of P2RY2 (Beta = −0.010, SE = 0.003, P = 0.007), suggesting that it may play a role in regulating endothelial nitric oxide synthase (eNOS) phosphorylation which may in turn influence arterial stiffness. This association may be mediated through an effect on endothelial function independent of the classical risk factors. Uridine is known to be an agonist of P2 receptors, particularly the P2Y subclass which consists of eight known human P2Y receptors (1, 2, 4, 6, and 11–14). This is relevant as P2Y1, P2Y2, and possibly P2Y4 are the purinergic receptors involved in eNOS phosphorylation during endothelial activation .
Uridine is a nucleotide base used as a dietary supplement for increasing the synthesis of cellular membranes and for other neurological properties. Uridine is able to exert an acute cardioprotective effect against myocardial ischaemia when preloaded, which is abolished by blocking potassium channels on the mitochondria (with 5-hydroxydecanoate); it appears that uridine preloading preserves the levels of energy metabolites (ATP, creatine phosphate, and uridine) and subsequently reduced lipid peroxidation . The nominal association with expression levels of a purinergic receptor suggests that uridine may exert its role on PWV via endothelial dysfunction .
Serine is one of the naturally occurring amino acids and it is synthesized in the body from other metabolites. It participates in the biosynthesis of purines and pyrimidines, and is also the precursor to numerous other metabolites, including sphingolipids and folate, the principal donor of one-carbon fragments in biosynthesis. In our data, we find that circulating levels of serine are significantly associated with lower PWV and with expression levels of a probe on the cystatin-like 1 (CSTL1) gene on chromosome 20 (Beta = 0.09, SE = 0.02, P = 4.26 × 10−8). The cystatin locus on chromosome 20 contains the majority of the type 2 cystatin genes and pseudogenes, and has been associated with cerebral haemorrhage and cerebritis .
Phenylacetylglutamine is a major nitrogenous metabolite that accumulates in uraemia . It is the glutamine conjugate of phenylacetic acid produced in humans and is also a well known gut microbial cometabolite whose levels are significantly different between Asian and North American individuals . Although Holmes et al. did not measure the correlation between phenylacetylglutamine and blood pressure (BP), they did test that between BP and hippurate (another microbial cometabolite significantly different between Asians and Caucasians, positively correlated with phenylacetylglutamine) and reported a negative correlation with DBP. Such published results are consistent with our findings, that is, a negative correlation between phenylacetylglutamine and SBP, cardiovascular risk and PWV. In our data, we also find a weak negative correlation between hippurate and SBP (Beta = −0.97, SE = 0.39, P < 0.012).
Phenylacetylglutamine levels in our data are correlated with adipocyte gene-expression levels of the cell death activator CIDE (CIDEC: Beta = 0.05, SE = 0.01, P = 6.97 × 10−8). This gene is regulated by insulin and its expression is positively correlated with insulin sensitivity . Mutations in this gene may contribute to insulin-resistant diabetes . CIDEC plays an important role in energy metabolism and lipid droplet formation , and its hepatic expression is increased in obese humans and is downregulated by marked weight loss .
Using metabolomic profiling, we searched for the molecular markers and the mechanisms involved in differences in PWV in women in order to investigate the regulation of arterial stiffening. We identified 12 blood metabolites, mainly amino acids, with high statistical significance associated with PWV. We also report three metabolites amongst those identified to be independently associated with PWV: uridine, serine, and phenylacetylglutamine achieving an R2 of 30%. Of the three metabolites identified, the one showing the strongest association is uridine. Uridine triphosphate (UTP, which unfortunately is not measured by the current metabolomic panel) stimulates vasodilatation, automaticity in ventricular myocytes, and release of tissue-plasminogen activator, indicating that UTP may be important in cardiac regulation . Uridine levels may be reflecting lower UTP levels or they may be cardioprotective via some other mechanism. Interestingly, however, uridine levels are not correlated with Framingham risk. We also report that circulating levels of serine are associated with PWV but not with Framingham risk. Our data indicate that the mechanisms underlying the association of PWV with uridine and serine are likely to be independent of the traditional CVD risk factors. Recent studies have shown that PWV improves cardiovascular event prediction [31,32]. Our data are, therefore, consistent and suggest that there are molecular mechanisms related to arterial stiffening and cardiovascular mortality that are not fully encompassed by the traditional cardiovascular risk factors.
We also report a novel association between phenylacetylglutamine and both Framingham risk scores and PWV. The negative correlation with phenylacetylglutamine is consistent with the previous reports on gut-microbiome-derived metabolites  and BP. We find that this metabolite is strongly associated with the gene-expression levels of CIDEC, a gene related to insulin resistance, suggesting that this metabolite may be related to this pathway. Phenylacetylglutamine is related to the gut microbiome composition and a number of reports have linked the function of gut bacteria to insulin resistance. Therefore, it is possible that CIDEC may be linking these two pathways and suggests a new mechanism linking insulin resistance and gut microbiome in BP regulation.
Our study, therefore, on one hand is consistent with the current knowledge of insulin resistance and endothelial activation mechanisms in determining arterial stiffening. On the other hand, our data suggest that some of these mechanisms may be related to the mechanisms that deserve further exploration, for example, those that link with the gut microbiome. These data also suggest that serine and uridine levels are linked to arterial stiffening, possibly via endothelial dysfunction but in a way that is not reflected directly on the traditional CVD risk factors.
The current study has several strengths. It used a nontargeted metabolomic approach that identifies a wide range of biochemicals in addition to lipids. TwinsUK is a very large and accurately phenotyped population, and this allowed us to explore the potential confounders (e.g. diet and MAP). The availability of expression and genetic data enabled us to explore some of the biological implications of the three metabolites identified.
We note some study limitations. Our study sample consisted of women only, and some metabolites could be influenced by sex-specific hormones. In addition, previous studies suggested that traditional risk factors  are less reliable in predicting risk in women than in men. We have only tested individuals of European descent, and the levels of one of the metabolites identified (phenylacetylglutamine) are known to vary between Asians and Caucasians. More importantly, because of the novelty of the phenotypes, we could not validate our results in an independent cohort. The cross-sectional nature of our data does not allow us to draw conclusions as to whether the metabolites identified are causative of arterial stiffness or merely correlated with it. However, our results highlight the relevance of investigating the molecular pathways related to PWV as this may lead to the identification of molecular mechanisms involved in cardiovascular diseases, in particular linked to endothelial activation, that act through other pathways. The identification of key metabolites related to PWV should encourage further research into this field.
The authors wish to express their appreciation to all the study participants of the TwinsUK study.
Sources of funding This work was supported by the EU Framework Programme 7 small-scale focused research collaborative project EurHEALTHAging 277849; Metabolomic analysis was funded by Pfizer; TwinsUK was funded by the Wellcome Trust; European Community's Seventh Framework Programme (FP7/2007–2013). This study also receives support from the National Institute for Health Research (NIHR) Clinical Research Facility at Guy's and St Thomas’ NHS Foundation Trust and NIHR Biomedical Research Centre based at Guy's and St Thomas’ NHS Foundation Trust and King's College London. BHF Special Project Grant SP/08/005/25115 was given to S.P.
Conflicts of interest
R.P.M. is an employee of Metabolon, Inc.
M.J.B. is a full-time employee of Pfizer; J.T. and M.J.B. are share holders of Pfizer.
Reviewer's Summary Evaluation Reviewer 1
This original study proposed by Menni et al. points to molecular pathways and biological processes involved in arterial stiffness in women. The identification of biological mechanisms affecting arterial viscoelastic properties paves the way for eventual specific treatment.
1. Vlachopoulos C, Aznaouridis K, Stefanadis C. Prediction of cardiovascular events and all-cause mortality with arterial stiffness: a systematic review and meta-analysis. J Am Coll Cardiol
2. Willum-Hansen T, Staessen JA, Torp-Pedersen C, Rasmussen S, Thijs L, Ibsen H, et al. Prognostic value of aortic pulse wave velocity as index of arterial stiffness in the general population. Circulation
3. Payne RA, Wilkinson IB, Webb DJ. Arterial stiffness and hypertension: emerging concepts. Hypertension
4. Marti CN, Gheorghiade M, Kalogeropoulos AP, Georgiopoulou VV, Quyyumi AA, Butler J. Endothelial dysfunction, arterial stiffness, and heart failure. J Am Coll Cardiol
5. Cecelja M, Jiang B, McNeill K, Kato B, Ritter J, Spector T, et al. Increased wave reflection rather than central arterial stiffness is the main determinant of raised pulse pressure in women and relates to mismatch in arterial dimensions: a twin study. J Am Coll Cardiol
6. McEniery CM, Yasmin, Hall IR, Qasem A, Wilkinson IB, Cockcroft JR. Normal vascular aging: differential effects on wave reflection and aortic pulse wave velocity: the Anglo-Cardiff Collaborative Trial (ACCT). J Am Coll Cardiol
7. Laurent S, Cockcroft J, Van Bortel L, Boutouyrie P, Giannattasio C, Hayoz D, et al. Expert consensus document on arterial stiffness: methodological issues and clinical applications. Eur Heart J
8. Giannattasio C, Failla M, Capra A, Scanziani E, Amigoni M, Boffi L, et al. Increased arterial stiffness in normoglycemic normotensive offspring of type 2 diabetic parents. Hypertension
9. Riggio S, Mandraffino G, Sardo MA, Iudicello R, Camarda N, Imbalzano E, et al. Pulse wave velocity and augmentation index, but not intima-media thickness, are early indicators of vascular damage in hypercholesterolemic children. Eur J Clin Invest
10. Giannattasio C, Mangoni AA, Failla M, Carugo S, Stella ML, Stefanoni P, et al. Impaired radial artery compliance in normotensive subjects with familial hypercholesterolemia. Atherosclerosis
11. Karras A, Haymann JP, Bozec E, Metzger M, Jacquot C, Maruani G, et al. Large artery stiffening and remodeling are independently associated with all-cause mortality and cardiovascular events in chronic kidney disease. Hypertension
12. Hollywood K, Brison DR, Goodacre R. Metabolomics: current technologies and future trends. Proteomics
13. Anderson SG, Sanders TA, Cruickshank JK. Plasma fatty acid composition as a predictor of arterial stiffness and mortality. Hypertension
14. D’Agostino RB Sr, Vasan RS, Pencina MJ, Wolf PA, Cobain M, Massaro JM, et al. General cardiovascular risk profile for use in primary care: the Framingham Heart Study. Circulation
15. Moayyeri A, Hammond CJ, Valdes AM, Spector TD. Cohort profile: TwinsUK and Healthy Ageing Twin Study. Int J Epidemiol
2013; 42:76–85.doi: 10.1093/ije/dyr207.
16. Menni C, Kastenmuller G, Petersen AK, Bell JT, Psatha M, Tsai PC, et al. Metabolomic markers reveal novel pathways of ageing and early development in human populations. Int J Epidemiol
17. Grundberg E, Small KS, Hedman AK, Nica AC, Buil A, Keildson S, et al. Mapping cis- and trans-regulatory effects across multiple tissues in twins. Nat Genet
18. Menni C, Zhai G, Macgregor A, Prehn C, Romisch-Margl W, Suhre K, et al. Targeted metabolomics profiles are strongly correlated with nutritional patterns in women. Metabolomics
19. Lunder M, Janic M, Kejzar N, Sabovic M. Associations among different functional and structural arterial wall properties and their relations to traditional cardiovascular risk factors in healthy subjects: a cross-sectional study. BMC Cardiovasc Disord
20. Da Silva CG, Specht A, Wegiel B, Ferran C, Kaczmarek E. Mechanism of purinergic activation of endothelial nitric oxide synthase in endothelial cells. Circulation
21. Krylova IB, Bulion VV, Selina EN, Mironova GD, Sapronov NS. Effect of uridine on energy metabolism, LPO, and antioxidant system in the myocardium under conditions of acute coronary insufficiency. Bull Exp Biol Med
22. McEniery CM, Wallace S, Mackenzie IS, McDonnell B, Yasmin, Newby DE, et al. Endothelial function is associated with pulse pressure, pulse wave velocity, and augmentation index in healthy humans. Hypertension
23. Levy E, Jaskolski M, Grubb A. The role of cystatin C in cerebral amyloid angiopathy and stroke: cell biology and animal models. Brain Pathol
24. Zimmerman L, Egestad B, Jornvall H, Bergstrom J. Identification and determination of phenylacetylglutamine, a major nitrogenous metabolite in plasma of uremic patients. Clin Nephrol
25. 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
26. Ito M, Nagasawa M, Omae N, Ide T, Akasaka Y, Murakami K. Differential regulation of CIDEA and CIDEC expression by insulin via Akt1/2- and JNK2-dependent pathways in human adipocytes. J Lipid Res
27. Rubio-Cabezas O, Puri V, Murano I, Saudek V, Semple RK, Dash S, et al. Partial lipodystrophy and insulin resistant diabetes in a patient with a homozygous nonsense mutation in CIDEC. EMBO Mol Med
28. Yonezawa T, Kurata R, Kimura M, Inoko H. Which CIDE are you on? Apoptosis and energy metabolism. Mol Biosyst
29. Hall AM, Brunt EM, Klein S, Finck BN. Hepatic expression of cell death-inducing DFFA-like effector C in obese subjects is reduced by marked weight loss. Obesity (Silver Spring)
30. Erlinge D, Harnek J, van Heusden C, Olivecrona G, Jern S, Lazarowski E. Uridine triphosphate (UTP) is released during cardiac ischemia. Int J Cardiol
31. Ben-Shlomo Y, Spears M, Boustred C, May M, Anderson SG, Benjamin EJ, et al. Aortic pulse wave velocity improves cardiovascular event prediction: an individual participant meta-analysis of prospective observational data from 17,635 subjects. J Am Coll Cardiol
32. Mitchell GF, Hwang SJ, Vasan RS, Larson MG, Pencina MJ, Hamburg NM, et al. Arterial stiffness and cardiovascular events: the Framingham Heart Study. Circulation
33. Kones R. Primary prevention of coronary heart disease: integration of new data, evolving views, revised goals, and role of rosuvastatin in management. A comprehensive survey. Drug Des Dev Ther