Search

Background/Objectives: This observational study dissects the complex temporal associations between body-mass index (BMI), waist-hip ratio (WHR) and circulating metabolomics using a combination of longitudinal and cross-sectional population-based datasets and new systems epidemiology tools. Subjects/Methods: Firstly, a data-driven subgrouping algorithm was employed to simplify high-dimensional metabolic profiling data into a single categorical variable: a self-organizing map (SOM) was created from 174 metabolic measures from cross-sectional surveys (FINRISK, n = 9708, ages 25–74) and a birth cohort (NFBC1966, n = 3117, age 31 at baseline, age 46 at follow-up) and an expert committee defined four subgroups of individuals based on visual inspection of the SOM. Secondly, the subgroups were compared regarding BMI and WHR trajectories in an independent longitudinal dataset: participants of the Young Finns Study (YFS, n = 1286, ages 24–39 at baseline, 10 years follow-up, three visits) were categorized into the four subgroups and subgroup-specific age-dependent trajectories of BMI, WHR and metabolic measures were modelled by linear regression. Results: The four subgroups were characterised at age 39 by high BMI, WHR and dyslipidemia (designated TG-rich); low BMI, WHR and favourable lipids (TG-poor); low lipids in general (Low lipid) and high low-density-lipoprotein cholesterol (High LDL-C). Trajectory modelling of the YFS dataset revealed a dynamic BMI divergence pattern: despite overlapping starting points at age 24, the subgroups diverged in BMI, fasting insulin (three-fold difference at age 49 between TG-rich and TG-poor) and insulin-associated measures such as triglyceride-cholesterol ratio. Trajectories also revealed a WHR progression pattern: despite different starting points at the age of 24 in WHR, LDL-C and cholesterol-associated measures, all subgroups exhibited similar rates of change in these measures, i.e. WHR progression was uniform regardless of the

Cardiometabolic diseases are a leading cause of morbidity and mortality worldwide. Although common risk factors are well known, there has been no success in changing the growing trend in these diseases. An earlier detection of underlying harmful changes in the body, as well as a deeper understanding of the biological changes related to these diseases, could help in the early prevention and development of novel treatment strategies. Different omics technologies combined with a large population and clinical cohorts allow for a more detailed exploration of the biological changes related to various diseases. In the current thesis, these methodologies were first utilised to explore the metabolic changes associated with incident stroke and its subtypes. The first study identified 10 metabolic measures associated with incident strokes. Second, metabolomic data, together with clinical variables, were used to investigate the metabolic characteristics of mortality risk among coronary angiography patients. Furthermore, the predictive properties of metabolic measures in mortality risk prediction were assessed. The second study demonstrated that metabolic changes associated with all-cause mortality were present also when controlling for important clinical variables. The metabolic measures slightly improved discrimination at the cost of calibration compared to clinical variables when assessing the utility of these measures in the prediction of all-cause mortality. Third, the potential role of circulating inflammatory cytokines as mediators in the pathway from increased adiposity measured by body mass index to relevant disease outcomes was assessed. The third study identified three cytokines driven by body mass index. Only weak associations were seen when further assessing the causal role of the identified cytokines in different inflammatory-related diseases. In conclusion, this thesis illustrates the utility of the application of omics data in three different epidemiolog, Tiivistelmä Kardiometaboliset sairaudet ovat suurimpia sairastavuuden ja kuolleisuuden aiheuttajia maailmanlaajuisesti. Vaikka yleisimmät riskitekijät näille sairauksille tunnetaan hyvin, kasvavaa trendiä näiden sairauksien osalta ei ole kyetty muuttamaan. Kehossa olevien piilevien muutosten aikaisempi tunnistaminen sekä näihin sairauksiin liittyvien biologisten muutosten syvempi ymmärtäminen voisi auttaa varhaisessa ehkäisyssä sekä uusien hoitomuotojen kehittämisessä. Erilaiset omiikkateknologiat yhdistettyinä suuriin kliinisiin ja väestöaineistoihin mahdollistavat sairauksiin liittyvien erilaisten biologisten muutosten entistä tarkemman tarkastelun. Tässä väitöskirjassa näitä menetelmiä on hyödynnetty ensinnäkin aivohalvauksiin ja sen alatyyppeihin liittyvien aineenvaihduntamuutosten selvittämiseen. Ensimmäinen osatyö tunnisti kymmenen aineenvaihduntamuuttujaa, jotka olivat yhteydessä aivohalvauksiin. Toisessa työssä aineenvaihduntamuuttujia käytettiin kuolleisuuteen yhteydessä olevien metabolisten muutoksien selvittämiseen sydämen varjoainekuvauspotilailla. Lisäksi metabolisten muuttujien tuomaa lisähyötyä kliinisten muuttujien lisäksi kuolleisuusriskin ennustamisessa arvioitiin. Toinen osatyö osoitti, että aineenvaihduntamuutoksia on havaittavissa, vaikka analyysissä huomioitaisiin tärkeät kliiniset muuttujat. Aineenvaihduntamuuttujat paransivat hieman ennustemallin erottelukykyä kokonaiskuolleisuuden ennustamisessa, mutta samalla systemaattinen harha hieman lisääntyi verrattuna pelkkien kliinisten muuttujien käyttöön. Kolmanneksi verenkierron tulehdusproteiinien roolia painoindeksillä mitatun liikalihavuuden ja erilaisten relevanttien sairauksien syy-seuraussuhteiden välittäjinä arvioitiin. Kolmas osatyö tunnisti kolme painoindeksiin yhteydessä olevaa tulehdusproteiinia. Arvioitaessa näiden tulehdusproteiinien syy-seuraussuhteita eri tulehduksellisiin sairauksiin, nähtiin vain heikkoja yhteyksiä. Yhteenvetona tämä väitöskirjatyö osoittaa erilaisten omiikka-aine

Circulating cell-free DNA (cf-DNA) has emerged as a promising biomarker of ageing, tissue damage and cellular stress. However, less is known about health behaviours, ageing phenotypes and metabolic processes that lead to elevated cf-DNA levels. We sought to analyse the relationship of circulating cf-DNA level to age, sex, smoking, physical activity, vegetable consumption, ageing phenotypes (physical functioning, the number of diseases, frailty) and an extensive panel of biomarkers including blood and urine metabolites and inflammatory markers in three human cohorts (N = 5385; 17‐82 years). The relationships were assessed using correlation statistics, and linear and penalised regressions (the Lasso), also stratified by sex. cf-DNA levels were significantly higher in men than in women, and especially in middle-aged men and women who smoke, and in older more frail individuals. Correlation statistics of biomarker data showed that cf-DNA level was higher with elevated inflammation (C-reactive protein, interleukin-6), and higher levels of homocysteine, and proportion of red blood cells and lower levels of ascorbic acid. Inflammation (C-reactive protein, glycoprotein acetylation), amino acids (isoleucine, leucine, tyrosine), and ketogenesis (3-hydroxybutyrate) were included in the cf-DNA level-related biomarker profiles in at least two of the cohorts. In conclusion, circulating cf-DNA level is different by sex, and related to health behaviour, health decline and metabolic processes common in health and disease. These results can inform future studies where epidemiological and biological pathways of cf-DNA are to be analysed in details, and for studies evaluating cf-DNA as a potential clinical marker.

Context: Aging varies between individuals, with profound consequences for chronic diseases and longevity. One hypothesis to explain the diversity is a genetically regulated molecular clock that runs differently between individuals. Large human studies with long enough follow-up to test the hypothesis are rare due to practical challenges, but statistical models of aging are built as proxies for the molecular clock by comparing young and old individuals cross-sectionally. These models remain untested against longitudinal data. Objective: We applied novel methodology to test if cross-sectional modeling can distinguish slow vs accelerated aging in a human population. Methods: We trained a machine learning model to predict age from 153 clinical and cardiometabolic traits. The model was tested against longitudinal data from another cohort. The training data came from cross-sectional surveys of the Finnish population (n = 9708; ages 25–74 years). The validation data included 3 time points across 10 years in the Young Finns Study (YFS; n = 1009; ages 24–49 years). Predicted metabolic age in 2007 was compared against observed aging rate from the 2001 visit to the 2011 visit in the YFS dataset and correlation between predicted vs observed metabolic aging was determined. Results: The cross-sectional proxy failed to predict longitudinal observations (R2 = 0.018%, P = 0.67). Conclusion: The finding is unexpected under the clock hypothesis that would produce a positive correlation between predicted and observed aging. Our results are better explained by a stratified model where aging rates per se are similar in adulthood but differences in starting points explain diverging metabolic fates.

A systematic comparison is presented for the effects of seven different normalization schemes in quantitative urinary metabolomics. Morning spot urine samples were analyzed with nuclear magnetic resonance (NMR) spectroscopy from a population-based group of 994 individuals. Forty-four metabolites were quantified and the metabolite–metabolite associations and the associations of metabolite concentrations with two representative clinical measures, body mass index and mean arterial pressure, were analyzed. Distinct differences were observed when comparing the effects of normalization for the intra-urine metabolite associations with those for the clinical associations. The metabolite–metabolite associations show quite complex patterns of similarities and dissimilarities between the different normalization methods, while the epidemiological association patterns are consistent, leading to the same overall biological interpretations. The results indicate that, in general, the normalization method appears to have only minor influences on standard epidemiological regression analyses with clinical/physiological measures. Multimetabolite normalization schemes showed consistent results with the customary creatinine reference. Nevertheless, interpretations of intra-urine metabolite associations and nuanced understanding of the epidemiological associations call for comparisons with different normalizations and accounting for the physiology, metabolism and kidney function related to the normalization schemes.

We assigned 329,908 UK Biobank participants into six subgroups based on a self‑organizing map of 51 biochemical measures (blinded for clinical outcomes). The subgroup with the most favorable metabolic traits was chosen as the reference. Hazard ratios (HR) for incident disease were modeled by Cox regression. Enrichment ratios (ER) of incident multi‑morbidity versus randomly expected co‑occurrence were evaluated by permutation tests; ER is like HR but captures co‑occurrence rather than event frequency. The subgroup with high urinary excretion without kidney stress (HR = 1.24) and the subgroup with the highest apolipoprotein B and blood pressure (HR = 1.52) were associated with ischemic heart disease (IHD). The subgroup with kidney stress, high adiposity and inflammation was associated with IHD (HR = 2.11), cancer (HR = 1.29), dementia (HR = 1.70) and mortality (HR = 2.12). The subgroup with high liver enzymes and triglycerides was at risk of diabetes (HR = 15.6). Multimorbidity was enriched in metabolically favorable subgroups (3.4 ≤ ER ≤ 4.0) despite lower disease burden overall; the relative risk of co‑occurring disease was higher in the absence of obvious metabolic dysfunction. These results provide synergistic insight into metabolic health and its associations with cardiovascular disease in a large population sample.

Background: The prevalence of atrial fibrillation (AF) is increasing with an aging worldwide population, yet a comprehensive understanding of its causes and consequences remains limited. We aim to assess the causes and consequences of AF via a bidirectional Mendelian randomization (MR) analysis. Methods: We used publicly available genome-wide association study (GWAS) summary data, centralized and harmonized by an open GWAS database. We assessed the genetically predicted effects of 5048 exposures on risk of AF, and the genetically predicted effects of genetic liability to AF, on 10 308 outcomes via two-sample MR analysis. Multivariable MR analysis was further conducted to explore the comparative roles of identified risk factors. Results: MR analysis suggested that 55 out of 5048 exposure traits, including four proteins, play a causal role in AF (P <1e-5 allowing for multiple comparisons). Multivariable analysis suggested that higher body mass index, height and systolic blood pressure as well as genetic liability to coronary artery diseases independently cause AF. Three out of the four proteins (DUSP13, TNFSF12 and IL6R) had a drug prioritizing score for atrial fibrillation of 0.26, 0.38 and 0.88, respectively (values closer to 1 indicating stronger evidence of the protein as a potential drug target). Genetic liability to AF was linked to a higher risk of cardio-embolic ischaemic stroke. Conclusions: Our results suggest body mass index, height, systolic blood pressure and genetic liability to coronary artery disease are independent causal risk factors for AF. Several proteins, including DUSP13, IL-6R and TNFSF12, may have therapeutic potential for AF.

Genetic and early environmental differences including early health habits associate with future health. To provide insight on the causal nature of these associations, monozygotic (MZ) twin pairs discordant for health habits provide an interesting natural experiment. Twin pairs discordant for leisure-time physical activity (LTPA) in early adult life is thus a powerful study design to investigate the associations between long-term LTPA and indicators of health and wellbeing. We have identified 17 LTPA discordant twin pairs from two Finnish twin cohorts and summarize key findings of these studies in this paper. The carefully characterized rare long-term LTPA discordant MZ twin pairs have participated in multi-dimensional clinical examinations. Key findings highlight that compared with less active twins in such MZ twin pairs, the twins with higher long-term LTPA have higher physical fitness, reduced body fat, reduced visceral fat, reduced liver fat, increased lumen diameters of conduit arteries to the lower limbs, increased bone mineral density in loaded bone areas, and an increased number of large high-density lipoprotein particles. The findings increase our understanding on the possible site-specific and system-level effects of long-term LTPA.

Background: Observational findings for high-density lipoprotein (HDL)-mediated cholesterol efflux capacity (HDL-CEC) and coronary heart disease (CHD) appear inconsistent, and knowledge of the genetic architecture of HDL-CEC is limited. Objectives: A large-scale observational study on the associations of HDL-CEC and other HDL-related measures with CHD and the largest genome-wide association study (GWAS) of HDL-CEC. Participants/Methods: Six independent cohorts were included with follow-up data for 14,438 participants to investigate the associations of HDL-related measures with incident CHD (1,570 events). The GWAS of HDL-CEC was carried out in 20,372 participants. Results: HDL-CEC did not associate with CHD when adjusted for traditional risk factors and HDL cholesterol (HDL-C). In contradiction, almost all HDL-related concentration measures associated consistently with CHD after corresponding adjustments. There were no genetic loci associated with HDL-CEC independent of HDL-C and triglycerides. Conclusions: HDL-CEC is not unequivocally associated with CHD in contrast to HDL-C, apolipoprotein A-I, and most of the HDL subclass particle concentrations.

Background: Low‐grade inflammation in the young may contribute to the early development of cardiovascular disease. We assessed whether circulating levels of glycoprotein acetyls (GlycA) were better able to predict the development of adverse cardiovascular disease risk profiles compared with the more commonly used biomarker high‐sensitivity CRP (C‐reactive protein). Methods: A total of 3306 adolescents and young adults from the Avon Longitudinal Study of Parents and Children (mean age, 15.4±0.3; n=1750) and Cardiovascular Risk in Young Finns Study (mean age, 32.1±5.0; n=1556) were included. Baseline associations between inflammatory biomarkers, body composition, cardiovascular risk factors, and subclinical measures of vascular dysfunction were assessed cross‐sectionally in both cohorts. Prospective risk of developing hypertension and metabolic syndrome during 9‐to‐10‐year follow‐up were also assessed as surrogate markers for future cardiovascular risk. GlycA showed greater within‐subject correlation over 9‐to‐10‐year follow‐up in both cohorts compared with CRP, particularly in the younger adolescent group (r=0.36 versus 0.07). In multivariable analyses, GlycA was found to associate with multiple lifestyle‐related cardiovascular disease risk factors, cardiometabolic risk factor burden, and vascular dysfunction (eg, mean difference in flow‐mediated dilation=−1.2 [−1.8, −0.7]% per z‐score increase). In contrast, CRP levels appeared predominantly driven by body mass index and showed little relationship to any measured cardiovascular risk factors or phenotypes. In both cohorts, only GlycA predicted future risk of both hypertension (risk ratio [RR], ≈1.1 per z‐score increase for both cohorts) and metabolic syndrome (RR, ≈1.2–1.3 per z‐score increase for both cohorts) in 9‐to‐10‐year follow‐up. Conclusions: Low‐grade inflammation captured by the novel biomarker GlycA is associated with adverse cardiovascular risk profiles from as early as adolescence and predicts f

Objective: Previous studies on the association between metabolic biomarkers and hypertension have been limited by small sample sizes, low number of studied biomarkers, and cross-sectional study design. In the largest study to date, we assess the cross-sectional and longitudinal associations between high-abundance serum biomarkers and blood pressure (BP). Methods: We studied cross-sectional (N = 36 985; age 50.5 ± 14.2; 53.1% women) and longitudinal (N = 4197; age 49.4 ± 11.8, 55.3% women) population samples of Finnish individuals. We included 53 serum biomarkers and other detailed lipoprotein subclass measures in our analyses. We studied the associations between serum biomarkers and BP using both conventional statistical methods and a machine learning algorithm (gradient boosting) while adjusting for clinical risk factors. Results: Fifty-one of 53 serum biomarkers were cross-sectionally related to BP (adjusted P < 0.05 for all). Conventional linear regression modeling demonstrated that LDL cholesterol, remnant cholesterol, apolipoprotein B, and acetate were positively, and HDL particle size was negatively, associated with SBP change over time (adjusted P < 0.05 for all). Adding serum biomarkers (cross-sectional root-mean-square error: 16.27 mmHg; longitudinal: 17.61 mmHg) in the model with clinical measures (cross-sectional: 16.70 mmHg; longitudinal 18.52 mmHg) improved the machine learning model fit. Glucose, albumin, triglycerides in LDL, glycerol, VLDL particle size, and acetoacetate had the highest importance scores in models related to current or future BP. Conclusions: Our results suggest that serum lipids, and particularly LDL-derived and VLDL-derived cholesterol measures, and glucose metabolism abnormalities are associated with hypertension onset. Use of serum metabolite determination could improve identification of individuals at high risk of developing hypertension.

Background: Individuals who are obese in childhood have an elevated risk of disease in adulthood. However, whether childhood adiposity directly impacts intermediate markers of this risk, independently of adult adiposity, is unclear. In this study, we have simultaneously evaluated the effects of childhood and adulthood body size on 123 systemic molecular biomarkers representing multiple metabolic pathways. Methods: Two-sample Mendelian randomization (MR) was conducted to estimate the causal effect of childhood body size on a total of 123 nuclear magnetic resonance-based metabolic markers using summary genome-wide association study (GWAS) data from up to 24 925 adults. Multivariable MR was then applied to evaluate the direct effects of childhood body size on these metabolic markers whilst accounting for adult body size. Further MR analyses were undertaken to estimate the potential mediating effects of these circulating metabolites on the risk of coronary artery disease (CAD) in adulthood using a sample of 60 801 cases and 123 504 controls. Results: Univariable analyses provided evidence that childhood body size has an effect on 42 of the 123 metabolic markers assessed (based on P < 4.07 × 10−4). However, the majority of these effects (35/42) substantially attenuated when accounting for adult body size using multivariable MR. We found little evidence that the biomarkers that were potentially influenced directly by childhood body size (leucine, isoleucine and tyrosine) mediate this effect onto adult disease risk. Very-low-density lipoprotein markers provided the strongest evidence of mediating the long-term effect of adiposity on CAD risk. Conclusions: Our findings suggest that childhood adiposity predominantly exerts its detrimental effect on adult systemic metabolism along a pathway that involves adulthood body size.

Objective: Oxidation of low-density lipoprotein (LDL) may promote atherosclerosis, whereas the reverse transport of oxidized lipids by high-density lipoprotein (HDL) may contribute to atheroprotection. To provide insights into the associations of lipoprotein lipid oxidation markers with lipoprotein subclasses at the population level, we investigated the associations of oxidized HDL lipids (oxHDLlipids) and oxidized LDL lipids (oxLDLlipids) with lipoprotein subclasses in a population-based cross-sectional study of 1395 Finnish adults ages 24–39 years. Methods: The analysis of oxidized lipids was based on the determination of the baseline level of conjugated dienes in lipoprotein lipids. A high-throughput nuclear magnetic resonance (NMR) platform was used to quantify circulating lipoprotein subclass concentrations and analyze their lipid compositions. Results: OxHDLlipids were mainly not associated with lipoprotein subclass lipid concentrations and lipid composition after adjustment for Apolipoprotein-A1 (Apo-A1), waist circumference and age. OxLDLlipids were associated with several markers of lipoprotein subclass lipid concentrations and composition after adjustment for Apolipoprotein-B (Apo-B), age and waist circumference. Several measures of HDL and LDL subclasses, including phospholipid and triglyceride composition, associated directly with oxLDLlipids. Cholesterol ester and free cholesterol composition in HDL and LDL associated inversely with oxLDLlipids. Conclusions: We conclude that these results do not support the idea that HDL’s particle size or composition would reflect its functional capacity in the reverse transport of oxidized lipids. On the contrary, oxLDLlipids were associated with the entire lipoprotein subclass profile, including numerous associations with the compositional descriptors of the particles. This is in line with the suggested role of LDL oxidation in atherogenesis.

Background: Genetic variants can be used to prioritize risk factors as potential therapeutic targets via Mendelian randomization (MR). An agnostic statistical framework using Bayesian model averaging (MR-BMA) can disentangle the causal role of correlated risk factors with shared genetic predictors. Here, our objective is to identify lipoprotein measures as mediators between lipid-associated genetic variants and coronary artery disease (CAD) for the purpose of detecting therapeutic targets for CAD. Methods: As risk factors we consider 30 lipoprotein measures and metabolites derived from a high-throughput metabolomics study including 24 925 participants. We fit multivariable MR models of genetic associations with CAD estimated in 453 595 participants (including 113 937 cases) regressed on genetic associations with the risk factors. MR-BMA assigns to each combination of risk factors a model score quantifying how well the genetic associations with CAD are explained. Risk factors are ranked by their marginal score and selected using false-discovery rate (FDR) criteria. We perform supplementary and sensitivity analyses varying the dataset for genetic associations with CAD. Results: In the main analysis, the top combination of risk factors ranked by the model score contains apolipoprotein B (ApoB) only. ApoB is also the highest ranked risk factor with respect to the marginal score (FDR <0.005). Additionally, ApoB is selected in all sensitivity analyses. No other measure of cholesterol or triglyceride is consistently selected otherwise. Conclusions: Our agnostic genetic investigation prioritizes ApoB across all datasets considered, suggesting that ApoB, representing the total number of hepatic-derived lipoprotein particles, is the primary lipid determinant of CAD.

Activation of the hypoxia-inducible factor (HIF) pathway reprograms energy metabolism. Hemoglobin (Hb) is the main carrier of oxygen. Using its normal variation as a surrogate measure for hypoxia, we explored whether lower Hb levels could lead to healthier metabolic profiles in mice and humans (nn = 7175) and used Mendelian randomization (MR) to evaluate potential causality (n = 173,480). The results showed evidence for lower Hb levels being associated with lower body mass index, better glucose tolerance and other metabolic profiles, lower inflammatory load, and blood pressure. Expression of the key HIF target genes SLC2A4 and Slc2a1 in skeletal muscle and adipose tissue, respectively, associated with systolic blood pressure in MR analyses and body weight, liver weight, and adiposity in mice. Last, manipulation of murine Hb levels mediated changes to key metabolic parameters. In conclusion, low-end normal Hb levels may be favorable for metabolic health involving mild chronic activation of the HIF response.

Background and purpose: Many drugs and environmental contaminants induce hypercholesterolemia and promote the risk of atherosclerotic cardiovascular disease. We tested the hypothesis that pregnane X receptor (PXR), a xenobiotic-sensing nuclear receptor, regulates the level of circulating atherogenic lipids in humans and utilized mouse experiments to identify the mechanisms involved. Experimental approach: We performed serum NMR metabolomics in healthy volunteers administered rifampicin, a prototypical human PXR ligand or placebo in a crossover setting. We used high-fat diet fed wild-type and PXR knockout mice to investigate the mechanisms mediating the PXR-induced alterations in cholesterol homeostasis. Key results: Activation of PXR induced cholesterogenesis both in pre-clinical and clinical settings. In human volunteers, rifampicin increased intermediate-density lipoprotein (IDL), low-density lipoprotein (LDL) and total cholesterol and lathosterol–cholesterol ratio, a marker of cholesterol synthesis, suggesting increased cholesterol synthesis. Experiments in mice indicated that PXR activation causes widespread induction of the cholesterol synthesis genes including the rate-limiting Hmgcr and upregulates the intermediates in the Kandutsch–Russell cholesterol synthesis pathway in the liver. Additionally, PXR activation induced plasma proprotein convertase subtilisin/kexin type 9 (PCSK9), a negative regulator of hepatic LDL uptake, in both mice and humans. We propose that these effects were mediated through increased proteolytic activation of sterol regulatory element-binding protein 2 (SREBP2) in response to PXR activation. Conclusions and implications: PXR activation induces cholesterol synthesis, elevating LDL and total cholesterol in humans. The PXR–SREBP2 pathway is a novel regulator of the cholesterol and PCSK9 synthesis and a molecular mechanism for drug- and chemical-induced hypercholesterolemia.

Our work is motivated by the search for metabolite quantitative trait loci (QTL) in a cohort of more than 5000 people. There are 158 metabolites measured by NMR spectroscopy in the 31-year follow-up of the Northern Finland Birth Cohort 1966 (NFBC66). These metabolites, as with many multivariate phenotypes produced by high-throughput biomarker technology, exhibit strong correlation structures. Existing approaches for combining such data with genetic variants for multivariate QTL analysis generally ignore phenotypic correlations or make restrictive assumptions about the associations between phenotypes and genetic loci. We present a computationally efficient Bayesian seemingly unrelated regressions model for high-dimensional data, with cell-sparse variable selection and sparse graphical structure for covariance selection. Cell sparsity allows different phenotype responses to be associated with different genetic predictors and the graphical structure is used to represent the conditional dependencies between phenotype variables. To achieve feasible computation of the large model space, we exploit a factorisation of the covariance matrix. Applying the model to the NFBC66 data with 9000 directly genotyped single nucleotide polymorphisms, we are able to simultaneously estimate genotype–phenotype associations and the residual dependence structure among the metabolites. The R package BayesSUR with full documentation is available at https://cran.r-project.org/web/packages/BayesSUR

Preterm birth (PTB) is associated with increased risk of type 2 diabetes and neurocognitive impairment later in life. We analyzed for the first time the associations of PTB with blood miRNA levels in adulthood. We also investigated the relationship of PTB associated miRNAs and adulthood phenotypes previously linked with premature birth. Blood MicroRNA profiling, genome-wide gene expression analysis, computer-based cognitive testing battery (CANTAB) and serum NMR metabolomics were performed for Young Finns Study subjects (aged 34–49 years, full-term n = 682, preterm n = 84). Preterm birth (vs. full-term) was associated with adulthood levels of hsa-miR-29b-3p in a fully adjusted regression model (p = 1.90 × 10–4, FDR = 0.046). The levels of hsa-miR-29b-3p were down-regulated in subjects with PTB with appropriate birthweight for gestational age (p = 0.002, fold change [FC] = − 1.20) and specifically in PTB subjects with small birthweight for gestational age (p = 0.095, FC = − 1.39) in comparison to individuals born full term. Hsa-miR-29b-3p levels correlated with the expressions of its target-mRNAs BCL11A and CS and the gene set analysis results indicated a target-mRNA driven association between hsa-miR-29b-3p levels and Alzheimer’s disease, Parkinson’s disease, Insulin signaling and Regulation of Actin Cytoskeleton pathway expression. The level of hsa-miR-29b-3p was directly associated with visual processing and sustained attention in CANTAB test and inversely associated with serum levels of VLDL subclass component and triglyceride levels. In conlcusion, adult blood levels of hsa-miR-29b-3p were lower in subjects born preterm. Hsa-miR-29b-3p associated with cognitive function and may be linked with adulthood morbidities in subjects born preterm, possibly through regulation of gene sets related to neurodegenerative diseases and insulin signaling as well as VLDL and triglyceride metabolism.

Background: The discovery of robust and trans-ethnically replicated DNA methylation markers of metabolic phenotypes, has hinted at a potential role of epigenetic mechanisms in lipid metabolism. However, DNA methylation and the lipid compositions and lipid concentrations of lipoprotein sizes have been scarcely studied. Here, we present an epigenome-wide association study (EWAS) (N = 5414 total) of mostly lipid-related metabolic measures, including a fine profiling of lipoproteins. As lipoproteins are the main players in the different stages of lipid metabolism, examination of epigenetic markers of detailed lipoprotein features might improve the diagnosis, prognosis, and treatment of metabolic disturbances. Results: We conducted an EWAS of leukocyte DNA methylation and 226 metabolic measurements determined by nuclear magnetic resonance spectroscopy in the population-based KORA F4 study (N = 1662) and replicated the results in the LOLIPOP, NFBC1966, and YFS cohorts (N = 3752). Follow-up analyses in the discovery cohort included investigations into gene transcripts, metabolic-measure ratios for pathway analysis, and disease endpoints. We identified 161 associations (p value < 4.7 × 10⁻¹⁰), covering 16 CpG sites at 11 loci and 57 metabolic measures. Identified metabolic measures were primarily medium and small lipoproteins, and fatty acids. For apolipoprotein B-containing lipoproteins, the associations mainly involved triglyceride composition and concentrations of cholesterol esters, triglycerides, free cholesterol, and phospholipids. All associations for HDL lipoproteins involved triglyceride measures only. Associated metabolic measure ratios, proxies of enzymatic activity, highlight amino acid, glucose, and lipid pathways as being potentially epigenetically implicated. Five CpG sites in four genes were associated with differential expression of transcripts in blood or adipose tissue. CpG sites in ABCG1 and PHGDH showed associations with metabolic measures

Background: Low socio-economic position (SEP) is a risk factor for multiple health outcomes, but its molecular imprints in the body remain unclear. Methods: We examined SEP as a determinant of serum nuclear magnetic resonance metabolic profiles in ∼30 000 adults and 4000 children across 10 UK and Finnish cohort studies. Results: In risk-factor-adjusted analysis of 233 metabolic measures, low educational attainment was associated with 37 measures including higher levels of triglycerides in small high-density lipoproteins (HDL) and lower levels of docosahexaenoic acid (DHA), omega-3 fatty acids, apolipoprotein A1, large and very large HDL particles (including levels of their respective lipid constituents) and cholesterol measures across different density lipoproteins. Among adults whose father worked in manual occupations, associations with apolipoprotein A1, large and very large HDL particles and HDL-2 cholesterol remained after adjustment for SEP in later life. Among manual workers, levels of glutamine were higher compared with non-manual workers. All three indicators of low SEP were associated with lower DHA, omega-3 fatty acids and HDL diameter. At all ages, children of manual workers had lower levels of DHA as a proportion of total fatty acids. Conclusions: Our work indicates that social and economic factors have a measurable impact on human physiology. Lower SEP was independently associated with a generally unfavourable metabolic profile, consistent across ages and cohorts. The metabolites we found to be associated with SEP, including DHA, are known to predict cardiovascular disease and cognitive decline in later life and may contribute to health inequalities.

The Special Turku Coronary Risk Factor Intervention Project (STRIP) is a prospective infancy-onset randomized dietary intervention trial targeting dietary fat quality and cholesterol intake, and favoring consumption of vegetables, fruit, and whole-grains. Diet (food records) and circulating metabolites were studied at six time points between the ages of 9–19 years (n = 549–338). Dietary targets for this study were defined as (1) the ratio of saturated fat (SAFA) to monounsaturated and polyunsaturated fatty acids (MUFA + PUFA) < 1:2, (2) intake of SAFA < 10% of total energy intake, (3) fiber intake ≥ 80th age-specific percentile, and (4) sucrose intake ≤ 20th age-specific percentile. Metabolic biomarkers were quantified by high-throughput nuclear magnetic resonance metabolomics. Better adherence to the dietary targets, regardless of study group allocation, was assoiated with higher serum proportion of PUFAs, lower serum proportion of SAFAs, and a higher degree of unsaturation of fatty acids. Achieving ≥ 1 dietary target resulted in higher low-density lipoprotein (LDL) particle size, lower circulating LDL subclass lipid concentrations, and lower circulating lipid concentrations in medium and small high-density lipoprotein subclasses compared to meeting 0 targets. Attaining more dietary targets (≥2) was associated with a tendency to lower lipid concentrations of intermediate-density lipoprotein and very low-density lipoprotein subclasses. Thus, adherence to dietary targets is favorably associated with multiple circulating fatty acids and lipoprotein subclass lipid concentrations, indicative of better cardio-metabolic health.

Accurate differentiation between neurodegenerative diseases is developing quickly and has reached an effective level in disease recognition. However, there has been less focus on effectively distinguishing the prodromal state from later dementia stages due to a lack of suitable biomarkers. We utilized the Disease State Index (DSI) machine learning classifier to see how well quantified metabolomics data compares to clinically used cerebrospinal fluid (CSF) biomarkers of Alzheimer’s disease (AD). The metabolic profiles were quantified for 498 serum and CSF samples using proton nuclear magnetic resonance spectroscopy. The patient cohorts in this study were dementia (with a clinical AD diagnosis) (N = 359), mild cognitive impairment (MCI) (N = 96), and control patients with subjective memory complaints (N = 43). DSI classification was conducted for MCI (N = 51) and dementia (N = 214) patients with low CSF amyloid-β levels indicating AD pathology and controls without such amyloid pathology (N = 36). We saw that the conventional CSF markers of AD were better at classifying controls from both dementia and MCI patients. However, quantified metabolic subclasses were more effective in classifying MCI from dementia. Our results show the consistent effectiveness of traditional CSF biomarkers in AD diagnostics. However, these markers are relatively ineffective in differentiating between MCI and the dementia stage, where the quantified metabolomics data provided significant benefit.

BACKGROUND: Chronic obstructive pulmonary disease (COPD) is a common lung disorder characterized by persistent and progressive airflow limitation as well as systemic changes. Metabolic changes in blood may help detect COPD in an earlier stage and predict prognosis. METHODS: We conducted a comprehensive study of circulating metabolites, measured by proton Nuclear Magnetic Resonance Spectroscopy, in relation with COPD and lung function. The discovery sample consisted of 5557 individuals from two large population-based studies in the Netherlands, the Rotterdam Study and the Erasmus Rucphen Family study. Significant findings were replicated in 12,205 individuals from the Lifelines-DEEP study, FINRISK and the Prospective Investigation of the Vasculature in Uppsala Seniors (PIVUS) studies. For replicated metabolites further investigation of causality was performed, utilizing genetics in the Mendelian randomization approach. RESULTS: There were 602 cases of COPD and 4955 controls used in the discovery meta-analysis. Our logistic regression results showed that higher levels of plasma Glycoprotein acetyls (GlycA) are significantly associated with COPD (OR = 1.16, P = 5.6 × 10- 4 in the discovery and OR = 1.30, P = 1.8 × 10- 6 in the replication sample). A bi-directional two-sample Mendelian randomization analysis suggested that circulating blood GlycA is not causally related to COPD, but that COPD causally i

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].

Circulating proteins are vital in human health and disease and are frequently used as biomarkers for clinical decision-making or as targets for pharmacological intervention. Here, we map and replicate protein quantitative trait loci (pQTL) for 90 cardiovascular proteins in over 30,000 individuals, resulting in 451 pQTLs for 85 proteins. For each protein, we further perform pathway mapping to obtain trans-pQTL gene and regulatory designations. We substantiate these regulatory findings with orthogonal evidence for trans-pQTLs using mouse knockdown experiments (ABCA1 and TRIB1) and clinical trial results (chemokine receptors CCR2 and CCR5), with consistent regulation. Finally, we evaluate known drug targets, and suggest new target candidates or repositioning opportunities using Mendelian randomization. This identifies 11 proteins with causal evidence of involvement in human disease that have not previously been targeted, including EGF, IL-16, PAPPA, SPON1, F3, ADM, CASP-8, CHI3L1, CXCL16, GDF15 and MMP-12. Taken together, these findings demonstrate the utility of large-scale mapping of the genetics of the proteome and provide a resource for future precision studies of circulating proteins in human health.

Markers of biological aging have potential utility in primary care and public health. We developed a model of age based on untargeted metabolic profiling across multiple platforms, including nuclear magnetic resonance spectroscopy and liquid chromatography–mass spectrometry in urine and serum, within a large sample (N = 2,239) from the UK Airwave cohort. We validated a subset of model predictors in a Finnish cohort including repeat measurements from 2,144 individuals. We investigated the determinants of accelerated aging, including lifestyle and psychological risk factors for premature mortality. The metabolomic age model was well correlated with chronological age (mean r = 0.86 across independent test sets). Increased metabolomic age acceleration (mAA) was associated after false discovery rate (FDR) correction with overweight/obesity, diabetes, heavy alcohol use and depression. DNA methylation age acceleration measures were uncorrelated with mAA. Increased DNA methylation phenotypic age acceleration (N = 1,110) was associated after FDR correction with heavy alcohol use, hypertension and low income. In conclusion, metabolomics is a promising approach for the assessment of biological age and appears complementary to established epigenetic clocks.

Background and aims: Apolipoprotein A-I (apoA-I) infusions represent a potential novel therapeutic approach for the prevention of coronary artery disease (CAD). Although circulating apoA-I concentrations inversely associate with risk of CAD, the evidence base of this representing a causal relationship is lacking. The aim was to assess the causal role of apoA-I using human genetics. Methods: We identified a variant (rs12225230) in APOA1 locus that associated with circulating apoA-I concentrations (p < 5 × 10−8) in 20,370 Finnish participants, and meta-analyzed our data with a previous GWAS of apoA-I. We obtained genetic estimates of CAD from UK Biobank and CARDIoGRAMplusC4D (totaling 122,733 CAD cases) and conducted a two-sample Mendelian randomization analysis. We compared our genetic findings to observational associations of apoA-I with risk of CAD in 918 incident CAD cases among 11,535 individuals from population-based prospective cohorts. Results: ApoA-I was associated with a lower risk of CAD in observational analyses (HR 0.81; 95%CI: 0.75, 0.88; per 1-SD higher apoA-I), with the association showing a dose-response relationship. Rs12225230 associated with apoA-I concentrations (per-C allele beta 0.076 SD; SE: 0.013; p = 1.5 × 10−9) but not with confounders. In Mendelian randomization analyses, apoA-I was not related to risk of CAD (OR 1.13; 95%CI: 0.98,1.30 per 1-SD higher apoA-I), which was different from the observational association. Similar findings were observed using an independent ABCA1 variant in sensitivity analysis. Conclusions: Genetic evidence fails to support a cardioprotective role for apoA-I. This is in line with the cumulative evidence showing that HDL-related phenotypes are unlikely to have a protective role in CAD.

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.

Apolipoprotein E (apoE) is the key regulator of plasma lipids, mediating altered functionalities in lipoprotein metabolism – affecting the risk of coronary artery (CAD) and Alzheimer’s diseases, as well as longevity. Searching pathways influenced by apoE prior to adverse manifestations, we utilized a metabolome dataset of 228 nuclear-magnetic-resonance-measured serum parameters with a 10-year follow-up from the population-based Young Finns Study cohort of 2,234 apoE-genotyped (rs7412, rs429358) adults, aged 24–39 at baseline. At the end of our follow-up, by limiting FDR-corrected p < 0.05, regression analyses revealed 180/228 apoE-polymorphism-related associations with the studied metabolites, in all subjects – without indications of apoE x sex interactions. Across all measured apoE- and apoB-containing lipoproteins, ε4 allele had consistently atherogenic and ε2 protective effect on particle concentrations of free/esterified cholesterol, triglycerides, phospholipids and total lipids. As novel findings, ε4 associated with glycoprotein acetyls, LDL-diameter and isoleucine – all reported biomarkers of CAD-risk, inflammation, diabetes and total mortality. ApoE-subgroup differences persisted through our 10-year follow-up, although some variation of individual metabolite levels was noticed. In conclusion, apoE polymorphism associate with a complex metabolic change, including aberrations in multiple novel biomarkers related to elevated cardiometabolic and all-cause mortality risk, extending our understanding about the role of apoE in health and disease.