Clinical and biochemical associations of urinary metabolites: quantitative epidemiological approach on renal-cardiometabolic biomarkers.

BACKGROUND: Urinary metabolomics has demonstrated considerable potential to assess kidney function and its metabolic corollaries in health and disease. However, applications in epidemiology remain sparse due to technical challenges. METHODS: We added 17 metabolites to an open-access urinary nuclear magnetic resonance metabolomics platform, extending the panel to 61 metabolites (n = 994). We also introduced automated quantification for 11 metabolites, extending the panel to 12 metabolites (+creatinine). Epidemiological associations between these 12 metabolites and 49 clinical measures were studied in three independent cohorts (up to 5989 participants). Detailed regression analyses with various confounding factors are presented for body mass index (BMI) and smoking. RESULTS: Sex-specific population reference concentrations and distributions are provided for 61 urinary metabolites (419 men and 575 women), together with methodological intra-assay metabolite variations as well as the biological intra-individual and epidemiological population variations. For the 12 metabolites, 362 associations were found. These are mostly novel and reflect potential molecular proxies to estimate kidney function, as the associations cannot be simply explained by estimated glomerular filtration rate. Unspecific renal excretion results in leakage of amino acids (and glucose) to urine in all individuals. Seven urinary metabolites associated with smoking, providing questionnaire-independent proxy measures of smoking status in epidemiological studies. Common confounders did not affect metabolite associations with smoking, but insulin had a clear effect on most associations with BMI, including strong effects on 2-hydroxyisobutyrate, valine, alanine, trigonelline and hippurate. CONCLUSIONS: Urinary metabolomics provides new insight on kidney function and related biomarkers on the renal-cardiometabolic system, supporting large-scale applications in epidemiology.

Apt interpretation of comprehensive lipoprotein data in large-scale epidemiology: disclosure of fundamental structural and metabolic relationships.

BACKGROUND: Quantitative lipoprotein analytics using nuclear magnetic resonance (NMR) spectroscopy is currently commonplace in large-scale studies. One methodology has become widespread and is currently being utilized also in large biobanks. It allows the comprehensive characterization of 14 lipoprotein subclasses, clinical lipids, apolipoprotein A-I and B. The details of these data are conceptualized here in relation to lipoprotein metabolism with particular attention on the fundamental characteristics of subclass particle numbers, lipid concentrations and compositional measures. METHODS AND RESULTS: The NMR methodology was applied to fasting serum samples from Northern Finland Birth Cohorts 1966 and 1986 with 5651 and 5605 participants, respectively. All results were highly consistent between the cohorts. Circulating lipid concentrations in a particular lipoprotein subclass arise predominantly as the result of the circulating number of those subclass particles. The spherical lipoprotein particle shape, with a radially oriented surface monolayer, imposes size-dependent biophysical constraints for the lipid composition of individual subclass particles and inherently restricts the accommodation of metabolic changes via compositional modifications. The new finding that the relationship between lipoprotein subclass particle concentrations and the particle size is log-linear reveals that circulating lipoprotein particles are also under rather strict metabolic constraints for both their absolute and relative concentrations. CONCLUSIONS: The fundamental structural and metabolic relationships between lipoprotein subclasses elucidated in this study empower detailed interpretation of lipoprotein metabolism. Understanding the intricate details of these extensive data is important for the precise interpretation of novel therapeutic opportunities and for fully utilizing the potential of forthcoming analyses of genetic and metabolic data in large biobanks.

Correction: Lipoprotein signatures of cholesteryl ester transfer protein and HMG-CoA reductase inhibition.

[This corrects the article DOI: 10.1371/journal.pbio.3000572.].

EpiMetal: an open-source graphical web browser tool for easy statistical analyses in epidemiology and metabolomics.

MOTIVATION: An intuitive graphical interface that allows statistical analyses and visualizations of extensive data without any knowledge of dedicated statistical software or programming. IMPLEMENTATION: EpiMetal is a single-page web application written in JavaScript, to be used via a modern desktop web browser. GENERAL FEATURES: Standard epidemiological analyses and self-organizing maps for data-driven metabolic profiling are included. Multiple extensive datasets with an arbitrary number of continuous and category variables can be integrated with the software. Any snapshot of the analyses can be saved and shared with others via a www-link. We demonstrate the usage of EpiMetal using pilot data with over 500 quantitative molecular measures for each sample as well as in two large-scale epidemiological cohorts (N >10 000). AVAILABILITY: The software usage exemplar and the pilot data are open access online at [http://EpiMetal.computationalmedicine.fi]. MIT licensed source code is available at the Github repository at [https://github.com/amergin/epimetal].

Data-driven multivariate population subgrouping via lipoprotein phenotypes versus apolipoprotein B in the risk assessment of coronary heart disease.

BACKGROUND AND AIMS: Population subgrouping has been suggested as means to improve coronary heart disease (CHD) risk assessment. We explored here how unsupervised data-driven metabolic subgrouping, based on comprehensive lipoprotein subclass data, would work in large-scale population cohorts. METHODS: We applied a self-organizing map (SOM) artificial intelligence methodology to define subgroups based on detailed lipoprotein profiles in a population-based cohort (n = 5789) and utilised the trained SOM in an independent cohort (n = 7607). We identified four SOM-based subgroups of individuals with distinct lipoprotein profiles and CHD risk and compared those to univariate subgrouping by apolipoprotein B quartiles. RESULTS: The SOM-based subgroup with highest concentrations for non-HDL measures had the highest, and the subgroup with lowest concentrations, the lowest risk for CHD. However, apolipoprotein B quartiles produced better resolution of risk than the SOM-based subgroups and also striking dose-response behaviour. CONCLUSIONS: These results suggest that the majority of lipoprotein-mediated CHD risk is explained by apolipoprotein B-containing lipoprotein particles. Therefore, even advanced multivariate subgrouping, with comprehensive data on lipoprotein metabolism, may not advance CHD risk assessment.

Lipoprotein signatures of cholesteryl ester transfer protein and HMG-CoA reductase inhibition.

Cholesteryl ester transfer protein (CETP) inhibition reduces vascular event risk, but confusion surrounds its effects on low-density lipoprotein (LDL) cholesterol. Here, we clarify associations of genetic inhibition of CETP on detailed lipoprotein measures and compare those to genetic inhibition of 3-hydroxy-3-methylglutaryl-coenzyme A reductase (HMGCR). We used an allele associated with lower CETP expression (rs247617) to mimic CETP inhibition and an allele associated with lower HMGCR expression (rs12916) to mimic the well-known effects of statins for comparison. The study consists of 65,427 participants of European ancestries with detailed lipoprotein subclass profiling from nuclear magnetic resonance spectroscopy. Genetic associations were scaled to 10% reduction in relative risk of coronary heart disease (CHD). We also examined observational associations of the lipoprotein subclass measures with risk of incident CHD in 3 population-based cohorts totalling 616 incident cases and 13,564 controls during 8-year follow-up. Genetic inhibition of CETP and HMGCR resulted in near-identical associations with LDL cholesterol concentration estimated by the Friedewald equation. Inhibition of HMGCR had relatively consistent associations on lower cholesterol concentrations across all apolipoprotein B-containing lipoproteins. In contrast, the associations of the inhibition of CETP were stronger on lower remnant and very-low-density lipoprotein (VLDL) cholesterol, but there were no associations on cholesterol concentrations in LDL defined by particle size (diameter 18-26 nm) (-0.02 SD LDL defined by particle size; 95% CI: -0.10 to 0.05 for CETP versus -0.24 SD, 95% CI -0.30 to -0.18 for HMGCR). Inhibition of CETP was strongly associated with lower proportion of triglycerides in all high-density lipoprotein (HDL) particles. In observational analyses, a higher triglyceride composition within HDL subclasses was associated with higher risk of CHD, independently of total cholesterol and triglycerides (strongest hazard ratio per 1 SD higher triglyceride composition in very large HDL 1.35; 95% CI: 1.18-1.54). In conclusion, CETP inhibition does not appear to affect size-specific LDL cholesterol but is likely to lower CHD risk by lowering concentrations of other atherogenic, apolipoprotein B-containing lipoproteins (such as remnant and VLDLs). Inhibition of CETP also lowers triglyceride composition in HDL particles, a phenomenon reflecting combined effects of circulating HDL, triglycerides, and apolipoprotein B-containing particles and is associated with a lower CHD risk in observational analyses. Our results reveal that conventional composite lipid assays may mask heterogeneous effects of emerging lipid-altering therapies.

Heat shock protein 90 is downregulated in calcific aortic valve disease.

BACKGROUND: Calcific aortic valve disease (CAVD) is an atheroinflammatory process; finally it leads to progressive calcification of the valve. There is no effective pharmacological treatment for CAVD and many of the underlying molecular mechanisms remain unknown. We conducted a proteomic study to reveal novel factors associated with CAVD. METHODS: We compared aortic valves from patients undergoing valvular replacement surgery due to non-calcified aortic insufficiency (control group, n = 5) to a stenotic group (n = 7) using two-dimensional difference gel electrophoresis (2D-DIGE). Protein spots were identified with mass spectrometry. Western blot and immunohistochemistry were used to validate the results in a separate patient cohort and Ingenuity Pathway Analysis (IPA) was exploited to predict the regulatory network of CAVD. RESULTS: We detected an upregulation of complement 9 (C9), serum amyloid P-component (APCS) and transgelin as well as downregulation of heat shock protein (HSP90), protein disulfide isomerase A3 (PDIA3), annexin A2 (ANXA2) and galectin-1 in patients with aortic valve stenosis. The decreased protein expression of HSP90 was confirmed with Western blot. CONCLUSIONS: We describe here a novel data set of proteomic changes associated with CAVD, including downregulation of the pro-inflammatory cytosolic protein, HSP90.

Direct Estimation of HDL-Mediated Cholesterol Efflux Capacity from Serum.

BACKGROUND: HDL-mediated cholesterol efflux capacity (HDL-CEC) is a functional attribute that may have a protective role in atherogenesis. However, the estimation of HDL-CEC is based on in vitro cell assays that are laborious and hamper large-scale phenotyping. METHODS: Here, we present a cost-effective high-throughput nuclear magnetic resonance (NMR) spectroscopy method to estimate HDL-CEC directly from serum. We applied the new method in a population-based study of 7603 individuals including 574 who developed incident coronary heart disease (CHD) during 15 years of follow-up, making this the largest quantitative study for HDL-CEC. RESULTS: As estimated by NMR-spectroscopy, a 1-SD higher HDL-CEC was associated with a lower risk of incident CHD (hazards ratio, 0.86; 95%CI, 0.79-0.93, adjusted for traditional risk factors and HDL-C). These findings are consistent with published associations based on in vitro cell assays. CONCLUSIONS: These corroborative large-scale findings provide further support for a potential protective role of HDL-CEC in CHD and substantiate this new method and its future applications.

Increased mesenchymal podoplanin expression is associated with calcification in aortic valves.

BACKGROUND AND AIM OF THE STUDY: Calcific aortic valve disease (CAVD) is a progressive disease starting from mild valvular sclerosis and progressing to severe aortic stenosis (AS) with calcified valves. The origin of the calcification is proposed to be mesenchymal cells which have differentiated towards an osteoblastic phenotype. Podoplanin is a glycoprotein expressed in the endothelium of lymphatic vessels and in osteoblasts and osteocytes, mesenchymal cells, as well as in many carcinomas and aortic atherosclerotic lesions. In CAVD, its expression has been evaluated only as a marker of the lymphatic vasculature. MATERIALS AND METHODS: We determined podoplanin expression in human aortic valves in four patient groups: control (C, n=7), aortic regurgitation (AR, n=8), aortic regurgitation and fibrosis (AR + f, n=15) and AS (n=49) by immunohistochemistry and quantitative real-time PCR (RT-PCR). RESULTS: Immunohistochemically, podoplanin expression was significantly increased in AR + f and AS groups when compared with the control and AR groups and the level of expression positively correlated with the extent of calcification and vascularity. Podoplanin mRNA levels were 1.7-fold higher in the AS group as compared with the control group (P=.05). Podoplanin-positivity was present not only in lymphatic vessel endothelium but also in osteoblasts, osteocytes, chondrocytes, macrophages and extracellular matrix. The majority of the podoplanin-positivity was in spindle cells with a myofibroblastic phenotype, often associated with calcifications. Tricuspid valves had more calcification-associated podoplanin than bi/unicuspid valves (median 1.52 vs 1.16, P<.001). CONCLUSIONS: CAVD is characterized by an increased expression of podoplanin; this is associated with the differentiation of valvular interstitial cells into calcium-producing, myofibroblast-like cells. In addition, tricuspid valves express relatively more podoplanin than bi/unicuspid valves.

Proof of concept for quantitative urine NMR metabolomics pipeline for large-scale epidemiology and genetics.

BACKGROUND: Quantitative molecular data from urine are rare in epidemiology and genetics. NMR spectroscopy could provide these data in high throughput, and it has already been applied in epidemiological settings to analyse urine samples. However, quantitative protocols for large-scale applications are not available. METHODS: We describe in detail how to prepare urine samples and perform NMR experiments to obtain quantitative metabolic information. Semi-automated quantitative line shape fitting analyses were set up for 43 metabolites and applied to data from various analytical test samples and from 1004 individuals from a population-based epidemiological cohort. Novel analyses on how urine metabolites associate with quantitative serum NMR metabolomics data (61 metabolic measures; n = 995) were performed. In addition, confirmatory genome-wide analyses of urine metabolites were conducted (n = 578). The fully automated quantitative regression-based spectral analysis is demonstrated for creatinine and glucose (n = 4548). RESULTS: Intra-assay metabolite variations were mostly <5%, indicating high robustness and accuracy of urine NMR spectroscopy methodology per se. Intra-individual metabolite variations were large, ranging from 6% to 194%. However, population-based inter-individual metabolite variations were even larger (from 14% to 1655%), providing a sound base for epidemiological applications. Metabolic associations between urine and serum were found to be clearly weaker than those within serum and within urine, indicating that urinary metabolomics data provide independent metabolic information. Two previous genome-wide hits for formate and 2-hydroxyisobutyrate were replicated at genome-wide significance. CONCLUSION: Quantitative urine metabolomics data suggest broad novelty for systems epidemiology. A roadmap for an open access methodology is provided.

Cellular Mechanisms of Valvular Thickening in Early and Intermediate Calcific Aortic Valve Disease.

BACKGROUND: Calcific aortic valve disease is common in an aging population. It is an active atheroinflammatory process that has an initial pathophysiology and similar risk factors as atherosclerosis. However, the ultimate disease phenotypes are markedly different. While coronary heart disease results in rupture-prone plaques, calcific aortic valve disease leads to heavily calcified and ossified valves. Both are initiated by the retention of low-density lipoprotein particles in the subendothelial matrix leading to sterile inflammation. In calcific aortic valve disease, the process towards calcification and ossification is preceded by valvular thickening, which can cause the first clinical symptoms. This is attributable to the accumulation of lipids, inflammatory cells and subsequently disturbances in the valvular extracellular matrix. Fibrosis is also increased but the innermost extracellular matrix layer is simultaneously loosened. Ultimately, the pathological changes in the valve cause massive calcification and bone formation - the main reasons for the loss of valvular function and the subsequent myocardial pathology. CONCLUSION: Calcification may be irreversible, and no drug treatments have been found to be effective, thus it is imperative to emphasize lifestyle prevention of the disease. Here we review the mechanisms underpinning the early stages of the disease.

Cardiac Actions of a Small Molecule Inhibitor Targeting GATA4-NKX2-5 Interaction.

Transcription factors are fundamental regulators of gene transcription, and many diseases, such as heart diseases, are associated with deregulation of transcriptional networks. In the adult heart, zinc-finger transcription factor GATA4 is a critical regulator of cardiac repair and remodelling. Previous studies also suggest that NKX2-5 plays function role as a cofactor of GATA4. We have recently reported the identification of small molecules that either inhibit or enhance the GATA4-NKX2-5 transcriptional synergy. Here, we examined the cardiac actions of a potent inhibitor (3i-1000) of GATA4-NKX2-5 interaction in experimental models of myocardial ischemic injury and pressure overload. In mice after myocardial infarction, 3i-1000 significantly improved left ventricular ejection fraction and fractional shortening, and attenuated myocardial structural changes. The compound also improved cardiac function in an experimental model of angiotensin II -mediated hypertension in rats. Furthermore, the up-regulation of cardiac gene expression induced by myocardial infarction and ischemia reduced with treatment of 3i-1000 or when micro- and nanoparticles loaded with 3i-1000 were injected intramyocardially or intravenously, respectively. The compound inhibited stretch- and phenylephrine-induced hypertrophic response in neonatal rat cardiomyocytes. These results indicate significant potential for small molecules targeting GATA4-NKX2-5 interaction to promote myocardial repair after myocardial infarction and other cardiac injuries.

Targeting vasoactive peptides for managing calcific aortic valve disease.

Calcific aortic valve disease (CAVD) represents a spectrum of disease spanning from milder degrees of calcification of valve leaflets, i.e., aortic sclerosis, to severe calcification i.e., aortic stenosis (AS) with hemodynamic instability. The prevalence of CAVD is increasing rapidly due to the aging of the population, being up to 2.8% among patients over 75 years of age. Even without significant aortic valve stenosis, aortic sclerosis is associated with a 50% increased risk of myocardial infarction and death from cardiovascular causes. To date, there is no pharmacological treatment available to reverse or hinder the progression of CAVD. So far, the cholesterol-lowering therapies (statins) and renin-angiotensin system (RAS) blocking drugs have been the major pharmacological agents investigated for treatment of CAVD. Especially angiotensin receptor blockers (ARB)s and angiotensin convertase enzyme inhibitors (ACEI)s, have been under active investigation in clinical trials, but have proven to be unsuccessful in slowing the progression of CAVD. Several studies have suggested that other vasoactive hormones, including endothelin and apelin systems are also associated with development of AS. In the present review, we discuss the role of vasoactive factors in the pathogenesis of CAVD as novel pharmacological targets for the treatment of aortic valve calcification. Key messages Vasoactive factors are involved in the progression of calcific aortic valve disease. Endothelin and renin-angiotensin systems seem to be most prominent targets for therapeutic interventions in the view of valvular pathogenesis. Circulating vasoactive factors may provide targets for diagnostic tools of calcified aortic valve disease.

Mitogen-activated protein kinase p38 target regenerating islet-derived 3γ expression is upregulated in cardiac inflammatory response in the rat heart.

Regenerating islet-derived 3γ (Reg3γ) is a multifunctional protein, associated with various tissue injuries and inflammatory states. Since chronic inflammation is characteristics also for heart failure, the aim of this study was to characterize Reg3γ expression in cardiac inflammatory conditions. Reg3γ expression was studied in experimental rat models of myocardial infarction (MI) and pressure overload in vivo. For cell culture studies neonatal rat cardiac myocytes (NRCMs) were used. In addition, adenovirus-mediated gene transfer of upstream mitogen-activated protein kinase (MAPK) kinase 3b and p38α MAPK in vivo and in vitro was performed. Reg3γ mRNA (12.8-fold, P < 0.01) and protein (5.8-fold, P < 0.001) levels were upregulated during the postinfarction remodeling at day 1 after MI, and angiotensin II (Ang II) markedly increased Reg3γ mRNA levels from 6 h to 2 weeks. Immunohistochemistry revealed that the Ang II-induced expression of Reg3γ was localized into the cardiac fibroblasts and myofibroblasts of the proliferating connective tissue in the heart. Stretching and treatments with endothelin-1, lipopolysaccharide (LPS), and fibroblast growth factor-1 increased Reg3γ mRNA levels in NRCMs. SB203580, a selective p38 MAPK inhibitor, markedly attenuated LPS and mechanical stretch-induced upregulation of Reg3γ gene expression. Moreover, combined overexpression of MKK3bE and WT p38α increased Reg3γ gene expression in cultured cardiomyocytes in vitro and in the rat heart in vivo. Our study shows that cardiac stress activates Reg3γ expression and p38 MAPK is an upstream regulator of Reg3γ gene expression in heart. Altogether our data suggest Reg3γ is associated with cardiac inflammatory signaling.

MicroRNA-125b and chemokine CCL4 expression are associated with calcific aortic valve disease.

Calcific aortic valve disease (CAVD) is a progressive pathological condition with no effective pharmacological therapy. To identify novel molecular pathways as potential targets for pharmacotherapy, we studied microRNA (miRNA) profiles of heavily stenotic aortic valves (AS). One of the most upregulated miRNAs in AS valves compared to control valves was miR-125b (1.4-fold; P < 0.05). To identify CAVD-related changes in gene expression, DNA microarray analysis was performed, including an intermediate fibro(sclero)tic stage of the disease. This revealed changes especially in genes related to inflammation and immune response, including chemokine (C-C motif) ligand 3 (CCL3) and 4 (CCL4). CCL3 mRNA level was increased 3.9-fold (P < 0.05) when AS valves were compared to control valves, and a 2.5-fold increase (P < 0.05) in CCL4 gene expression was observed when fibro(sclero)tic valves were compared to control valves. Both CCL3 and CCL4 localized to macrophages by immunofluorescence. To identify chemokine-miRNA target pairs, data from miRNA target prediction databases were combined with valvular miRNA and mRNA expression profiles. MiR-125b was computationally predicted to target CCL4, as confirmed experimentally in cultured human THP-1 macrophages. Collectively, miR-125b and CCL4 appear to be involved in the progression of CAVD and may offer novel therapeutic and diagnostic strategies related to this disease.

The Early-Onset Myocardial Infarction Associated PHACTR1 Gene Regulates Skeletal and Cardiac Alpha-Actin Gene Expression.

The phosphatase and actin regulator 1 (PHACTR1) locus is a very commonly identified hit in genome-wide association studies investigating coronary artery disease and myocardial infarction (MI). However, the function of PHACTR1 in the heart is still unknown. We characterized the mechanisms regulating Phactr1 expression in the heart, used adenoviral gene delivery to investigate the effects of Phactr1 on cardiac function, and analyzed the relationship between MI associated PHACTR1 allele and cardiac function in human subjects. Phactr1 mRNA and protein levels were markedly reduced (60%, P<0.01 and 90%, P<0.001, respectively) at 1 day after MI in rats. When the direct myocardial effects of Phactr1 were studied, the skeletal α-actin to cardiac α-actin isoform ratio was significantly higher (1.5-fold, P<0.05) at 3 days but 40% lower (P<0.05) at 2 weeks after adenovirus-mediated Phactr1 gene delivery into the anterior wall of the left ventricle. Similarly, the skeletal α-actin to cardiac α-actin ratio was lower at 2 weeks in infarcted hearts overexpressing Phactr1. In cultured neonatal cardiac myocytes, adenovirus-mediated Phactr1 overexpression for 48 hours markedly increased the skeletal α-actin to cardiac α-actin ratio, this being associated with an enhanced DNA binding activity of serum response factor. Phactr1 overexpression exerted no major effects on the expression of other cardiac genes or LV structure and function in normal and infarcted hearts during 2 weeks' follow-up period. In human subjects, MI associated PHACTR1 allele was not associated significantly with cardiac function (n = 1550). Phactr1 seems to regulate the skeletal to cardiac α-actin isoform ratio.

WDR12, a Member of Nucleolar PeBoW-Complex, Is Up-Regulated in Failing Hearts and Causes Deterioration of Cardiac Function.

AIMS: In a recent genome-wide association study, WD-repeat domain 12 (WDR12) was associated with early-onset myocardial infarction (MI). However, the function of WDR12 in the heart is unknown. METHODS AND RESULTS: We characterized cardiac expression of WDR12, used adenovirus-mediated WDR12 gene delivery to examine effects of WDR12 on left ventricular (LV) remodeling, and analyzed relationship between MI associated WDR12 allele and cardiac function in human subjects. LV WDR12 protein levels were increased in patients with dilated cardiomyopathy and rats post-infarction. In normal adult rat hearts, WDR12 gene delivery into the anterior wall of the LV decreased interventricular septum diastolic and systolic thickness and increased the diastolic and systolic diameters of the LV. Moreover, LV ejection fraction (9.1%, P<0.05) and fractional shortening (12.2%, P<0.05) were declined. The adverse effects of WDR12 gene delivery on cardiac function were associated with decreased cellular proliferation, activation of p38 mitogen-activated protein kinase (MAPK)/heat shock protein (HSP) 27 pathway, and increased protein levels of Block of proliferation 1 (BOP1), essential for ribosome biogenesis. Post-infarction WDR12 gene delivery decreased E/A ratio (32%, P<0.05) suggesting worsening of diastolic function. In human subjects, MI associated WDR12 allele was associated significantly with diastolic dysfunction and left atrial size. CONCLUSIONS: WDR12 triggers distinct deterioration of cardiac function in adult rat heart and the MI associated WDR12 variant is associated with diastolic dysfunction in human subjects.