Genome-wide characterization of circulating metabolic biomarkers.

Genome-wide association analyses using high-throughput metabolomics platforms have led to novel insights into the biology of human metabolism1-7. This detailed knowledge of the genetic determinants of systemic metabolism has been pivotal for uncovering how genetic pathways influence biological mechanisms and complex diseases8-11. Here we present a genome-wide association study for 233 circulating metabolic traits quantified by nuclear magnetic resonance spectroscopy in up to 136,016 participants from 33 cohorts. We identify more than 400 independent loci and assign probable causal genes at two-thirds of these using manual curation of plausible biological candidates. We highlight the importance of sample and participant characteristics that can have significant effects on genetic associations. We use detailed metabolic profiling of lipoprotein- and lipid-associated variants to better characterize how known lipid loci and novel loci affect lipoprotein metabolism at a granular level. We demonstrate the translational utility of comprehensively phenotyped molecular data, characterizing the metabolic associations of intrahepatic cholestasis of pregnancy. Finally, we observe substantial genetic pleiotropy for multiple metabolic pathways and illustrate the importance of careful instrument selection in Mendelian randomization analysis, revealing a putative causal relationship between acetone and hypertension. Our publicly available results provide a foundational resource for the community to examine the role of metabolism across diverse diseases.

Dissecting the polygenic basis of atherosclerosis via disease-associated cell state signatures.

Coronary artery disease (CAD) is a pandemic disease where up to half of the risk is explained by genetic factors. Advanced insights into the genetic basis of CAD require deeper understanding of the contributions of different cell types, molecular pathways, and genes to disease heritability. Here, we investigate the biological diversity of atherosclerosis-associated cell states and interrogate their contribution to the genetic risk of CAD by using single-cell and bulk RNA sequencing (RNA-seq) of mouse and human lesions. We identified 12 disease-associated cell states that we characterized further by gene set functional profiling, ligand-receptor prediction, and transcription factor inference. Importantly, Vcam1+ smooth muscle cell state genes contributed most to SNP-based heritability of CAD. In line with this, genetic variants near smooth muscle cell state genes and regulatory elements explained the largest fraction of CAD-risk variance between individuals. Using this information for variant prioritization, we derived a hybrid polygenic risk score (PRS) that demonstrated improved performance over a classical PRS. Our results provide insights into the biological mechanisms associated with CAD risk, which could make a promising contribution to precision medicine and tailored therapeutic interventions in the future.

Gene expression networks regulated by human personality.

Genome-wide association studies of human personality have been carried out, but transcription of the whole genome has not been studied in relation to personality in humans. We collected genome-wide expression profiles of adults to characterize the regulation of expression and function in genes related to human personality. We devised an innovative multi-omic approach to network analysis to identify the key control elements and interactions in multi-modular networks. We identified sets of transcribed genes that were co-expressed in specific brain regions with genes known to be associated with personality. Then we identified the minimum networks for the co-localized genes using bioinformatic resources. Subjects were 459 adults from the Young Finns Study who completed the Temperament and Character Inventory and provided peripheral blood for genomic and transcriptomic analysis. We identified an extrinsic network of 45 regulatory genes from seed genes in brain regions involved in self-regulation of emotional reactivity to extracellular stimuli (e.g., self-regulation of anxiety) and an intrinsic network of 43 regulatory genes from seed genes in brain regions involved in self-regulation of interpretations of meaning (e.g., production of concepts and language). We discovered that interactions between the two networks were coordinated by a control hub of 3 miRNAs and 3 protein-coding genes shared by both. Interactions of the control hub with proteins and ncRNAs identified more than 100 genes that overlap directly with known personality-related genes and more than another 4000 genes that interact indirectly. We conclude that the six-gene hub is the crux of an integrative network that orchestrates information-transfer throughout a multi-modular system of over 4000 genes enriched in liquid-liquid-phase-separation (LLPS)-related RNAs, diverse transcription factors, and hominid-specific miRNAs and lncRNAs. Gene expression networks associated with human personality regulate neuronal plasticity, epigenesis, and adaptive functioning by the interactions of salience and meaning in self-awareness.

Mitochondrial genome-wide analysis of nuclear DNA methylation quantitative trait loci.

Mitochondria have a complex communication network with the surrounding cell and can alter nuclear DNA methylation (DNAm). Variation in the mitochondrial DNA (mtDNA) has also been linked to differential DNAm. Genome-wide association studies have identified numerous DNAm quantitative trait loci, but these studies have not examined the mitochondrial genome. Herein, we quantified nuclear DNAm from blood and conducted a mitochondrial genome-wide association study of DNAm, with an additional emphasis on sex- and prediabetes-specific heterogeneity. We used the Young Finns Study (n = 926) with sequenced mtDNA genotypes as a discovery sample and sought replication in the Ludwigshafen Risk and Cardiovascular Health study (n = 2317). We identified numerous significant associations in the discovery phase (P < 10-9), but they were not replicated when accounting for multiple testing. In total, 27 associations were nominally replicated with a P < 0.05. The replication analysis presented no evidence of sex- or prediabetes-specific heterogeneity. The 27 associations were included in a joint meta-analysis of the two cohorts, and 19 DNAm sites associated with mtDNA variants, while four other sites showed haplogroup associations. An expression quantitative trait methylation analysis was performed for the identified DNAm sites, pinpointing two statistically significant associations. This study provides evidence of a mitochondrial genetic control of nuclear DNAm with little evidence found for sex- and prediabetes-specific effects. The lack of a comparable mtDNA data set for replication is a limitation in our study and further studies are needed to validate our results.

Bias in O-Information Estimation.

Higher-order relationships are a central concept in the science of complex systems. A popular method of attempting to estimate the higher-order relationships of synergy and redundancy from data is through the O-information. It is an information-theoretic measure composed of Shannon entropy terms that quantifies the balance between redundancy and synergy in a system. However, bias is not yet taken into account in the estimation of the O-information of discrete variables. In this paper, we explain where this bias comes from and explore it for fully synergistic, fully redundant, and fully independent simulated systems of n=3 variables. Specifically, we explore how the sample size and number of bins affect the bias in the O-information estimation. The main finding is that the O-information of independent systems is severely biased towards synergy if the sample size is smaller than the number of jointly possible observations. This could mean that triplets identified as highly synergistic may in fact be close to independent. A bias approximation based on the Miller-Maddow method is derived for the O-information. We find that for systems of n=3 variables the bias approximation can partially correct for the bias. However, simulations of fully independent systems are still required as null models to provide a benchmark of the bias of the O-information.

Meta-analysis of epigenome-wide association studies of carotid intima-media thickness.

Common carotid intima-media thickness (cIMT) is an index of subclinical atherosclerosis that is associated with ischemic stroke and coronary artery disease (CAD). We undertook a cross-sectional epigenome-wide association study (EWAS) of measures of cIMT in 6400 individuals. Mendelian randomization analysis was applied to investigate the potential causal role of DNA methylation in the link between atherosclerotic cardiovascular risk factors and cIMT or clinical cardiovascular disease. The CpG site cg05575921 was associated with cIMT (beta = -0.0264, p value = 3.5 × 10-8) in the discovery panel and was replicated in replication panel (beta = -0.07, p value = 0.005). This CpG is located at chr5:81649347 in the intron 3 of the aryl hydrocarbon receptor repressor gene (AHRR). Our results indicate that DNA methylation at cg05575921 might be in the pathway between smoking, cIMT and stroke. Moreover, in a region-based analysis, 34 differentially methylated regions (DMRs) were identified of which a DMR upstream of ALOX12 showed the strongest association with cIMT (p value = 1.4 × 10-13). In conclusion, our study suggests that DNA methylation may play a role in the link between cardiovascular risk factors, cIMT and clinical cardiovascular disease.

Robust multi-group gene set analysis with few replicates.

BACKGROUND: Competitive gene set analysis is a standard exploratory tool for gene expression data. Permutation-based competitive gene set analysis methods are preferable to parametric ones because the latter make strong statistical assumptions which are not always met. For permutation-based methods, we permute samples, as opposed to genes, as doing so preserves the inter-gene correlation structure. Unfortunately, up until now, sample permutation-based methods have required a minimum of six replicates per sample group. RESULTS: We propose a new permutation-based competitive gene set analysis method for multi-group gene expression data with as few as three replicates per group. The method is based on advanced sample permutation technique that utilizes all groups within a data set for pairwise comparisons. We present a comprehensive evaluation of different permutation techniques, using multiple data sets and contrast the performance of our method, mGSZm, with other state of the art methods. We show that mGSZm is robust, and that, despite only using less than six replicates, we are able to consistently identify a high proportion of the top ranked gene sets from the analysis of a substantially larger data set. Further, we highlight other methods where performance is highly variable and appears dependent on the underlying data set being analyzed. CONCLUSIONS: Our results demonstrate that robust gene set analysis of multi-group gene expression data is permissible with as few as three replicates. In doing so, we have extended the applicability of such approaches to resource constrained experiments where additional data generation is prohibitively difficult or expensive. An R package implementing the proposed method and supplementary materials are available from the website http://ekhidna.biocenter.helsinki.fi/downloads/pashupati/mGSZm.html .

Polygenic risk for schizophrenia, social dispositions, and pace of epigenetic aging: Results from the Young Finns Study.

Schizophrenia is often regarded as a disorder of premature aging. We investigated (a) whether polygenic risk for schizophrenia (PRSsch ) relates to pace of epigenetic aging and (b) whether personal dispositions toward active and emotionally close relationships protect against accelerated epigenetic aging in individuals with high PRSsch . The sample came from the population-based Young Finns Study (n = 1348). Epigenetic aging was measured with DNA methylation aging algorithms such as AgeAccelHannum , EEAAHannum , IEAAHannum , IEAAHorvath , AgeAccelHorvath , AgeAccelPheno , AgeAccelGrim , and DunedinPACE. A PRSsch was calculated using summary statistics from the most comprehensive genome-wide association study of schizophrenia to date. Social dispositions were assessed in terms of extraversion, sociability, reward dependence, cooperativeness, and attachment security. We found that PRSsch did not have a statistically significant effect on any studied indicator of epigenetic aging. Instead, PRSsch had a significant interaction with reward dependence (p = 0.001-0.004), cooperation (p = 0.009-0.020), extraversion (p = 0.019-0.041), sociability (p = 0.003-0.016), and attachment security (p = 0.007-0.014) in predicting AgeAccelHannum , EEAAHannum , or IEAAHannum . Specifically, participants with high PRSsch appeared to display accelerated epigenetic aging at higher (vs. lower) levels of extraversion, sociability, attachment security, reward dependence, and cooperativeness. A rather opposite pattern was evident for those with low PRSsch . No such interactions were evident when predicting the other indicators of epigenetic aging. In conclusion, against our hypothesis, frequent social interactions may relate to accelerated epigenetic aging in individuals at risk for psychosis. We speculate that this may be explained by social-cognitive impairments (perceiving social situations as overwhelming or excessively arousing) or ending up in less supportive or deviant social groups.

Identification of gene networks jointly associated with depressive symptoms and cardiovascular health metrics using whole blood transcriptome in the Young Finns Study.

Background: Studies have shown that cardiovascular health (CVH) is related to depression. We aimed to identify gene networks jointly associated with depressive symptoms and cardiovascular health metrics using the whole blood transcriptome. Materials and methods: We analyzed human blood transcriptomic data to identify gene co-expression networks, termed gene modules, shared by Beck's depression inventory (BDI-II) scores and cardiovascular health (CVH) metrics as markers of depression and cardiovascular health, respectively. The BDI-II scores were derived from Beck's Depression Inventory, a 21-item self-report inventory that measures the characteristics and symptoms of depression. CVH metrics were defined according to the American Heart Association criteria using seven indices: smoking, diet, physical activity, body mass index (BMI), blood pressure, total cholesterol, and fasting glucose. Joint association of the modules, identified with weighted co-expression analysis, as well as the member genes of the modules with the markers of depression and CVH were tested with multivariate analysis of variance (MANOVA). Results: We identified a gene module with 256 genes that were significantly correlated with both the BDI-II score and CVH metrics. Based on the MANOVA test results adjusted for age and sex, the module was associated with both depression and CVH markers. The three most significant member genes in the module were YOD1, RBX1, and LEPR. Genes in the module were enriched with biological pathways involved in brain diseases such as Alzheimer's, Parkinson's, and Huntington's. Conclusions: The identified gene module and its members can provide new joint biomarkers for depression and CVH.

Predicting early-stage coronary artery disease using machine learning and routine clinical biomarkers improved by augmented virtual data.

Aims: Coronary artery disease (CAD) is a highly prevalent disease with modifiable risk factors. In patients with suspected obstructive CAD, evaluating the pre-test probability model is crucial for diagnosis, although its accuracy remains controversial. Machine learning (ML) predictive models can help clinicians detect CAD early and improve outcomes. This study aimed to identify early-stage CAD using ML in conjunction with a panel of clinical and laboratory tests. Methods and results: The study sample included 3316 patients enrolled in the Ludwigshafen Risk and Cardiovascular Health (LURIC) study. A comprehensive array of attributes was considered, and an ML pipeline was developed. Subsequently, we utilized five approaches to generating high-quality virtual patient data to improve the performance of the artificial intelligence models. An extension study was carried out using data from the Young Finns Study (YFS) to assess the results' generalizability. Upon applying virtual augmented data, accuracy increased by approximately 5%, from 0.75 to -0.79 for random forests (RFs), and from 0.76 to -0.80 for Gradient Boosting (GB). Sensitivity showed a significant boost for RFs, rising by about 9.4% (0.81-0.89), while GB exhibited a 4.8% increase (0.83-0.87). Specificity showed a significant boost for RFs, rising by ∼24% (from 0.55 to 0.70), while GB exhibited a 37% increase (from 0.51 to 0.74). The extension analysis aligned with the initial study. Conclusion: Accurate predictions of angiographic CAD can be obtained using a set of routine laboratory markers, age, sex, and smoking status, holding the potential to limit the need for invasive diagnostic techniques. The extension analysis in the YFS demonstrated the potential of these findings in a younger population, and it confirmed applicability to atherosclerotic vascular disease.

Identification and functional characterisation of DNA methylation differences between East- and West-originating Finns.

Eastern and Western Finns show a striking difference in coronary heart disease-related mortality; genetics is a known contributor for this discrepancy. Here, we discuss the potential role of DNA methylation in mediating the discrepancy in cardiometabolic disease-risk phenotypes between the sub-populations. We used data from the Young Finns Study (n = 969) to compare the genome-wide DNA methylation levels of East- and West-originating Finns. We identified 21 differentially methylated loci (FDR < 0.05; Δß >2.5%) and 7 regions (smoothed FDR < 0.05; CpGs ≥ 5). Methylation at all loci and regions associates with genetic variants (p < 5 × 10-8). Independently of genetics, methylation at 11 loci and 4 regions associates with transcript expression, including genes encoding zinc finger proteins. Similarly, methylation at 5 loci and 4 regions associates with cardiometabolic disease-risk phenotypes including triglycerides, glucose, cholesterol, as well as insulin treatment. This analysis was also performed in LURIC (n = 2371), a German cardiovascular patient cohort, and results replicated for the association of methylation at cg26740318 and DMR_11p15 with diabetes-related phenotypes and methylation at DMR_22q13 with triglyceride levels. Our results indicate that DNA methylation differences between East and West Finns may have a functional role in mediating the cardiometabolic disease discrepancy between the sub-populations.

Genetic risk score for coronary artery calcification and its predictive ability for coronary artery disease.

Aim: The modest added predictive value of the existing genetic risk scores (GRSs) for coronary artery disease (CAD) could be partly due to missing genetic components, hidden in the genetic architecture of intermediate phenotypes such as coronary artery calcification (CAC). In this study, we investigated the predictive ability of CAC GRS for CAD. Materials and methods: We investigated the association of CAC GRSs with CAD and coronary calcification among the participants in the Ludwigshafen Risk and Cardiovascular Health study (LURIC) (n = 2742), the Tampere Vascular Study (TVS) (n = 133), and the Tampere Sudden Death Study (TSDS) (n = 660) using summary data from the largest multi-ancestry GWAS meta-analysis of CAC to date. Added predictive value of the CAC GRS over the traditional CVD risk factors as well as metaGRS, a GRS for CAD constructed with 1.7 million genetic variants, was tested with standard train-test machine learning approach using the LURIC data, which had the largest sample size. Results: CAC GRS was significantly associated with CAD in LURIC (OR=1.41, 95 % CI [1.28-1.55]), TVS (OR=1.79, 95 % CI [1.05-3.21]) as well as in TSDS (OR=4.20, 95 % CI [1.74-10.52]). CAC GRS showed strong association with calcification areas in left (OR=1.78, 95 % CI [1.16-2.74]) and right (OR=1.71, 95 % CI [1.98-2.67]) coronary arteries. There was statistically significant added predictive value of the CAC GRS for CAD over the used traditional CVD risk factors (AUC 0.734 vs 0.717, p-value = 0.02). Furthermore, CAC GRS improved the prediction accuracy for CAD when combined with metaGRS. Conclusions: This study showed that CAC GRS is a new risk marker for CAD in three European cohorts, with added predictive value over the traditional CVD risk factors.

Early resilience and epigenetic ageing: Results from the prospective Young Finns Study with a 31-year follow-up.

Evidence is accumulating on the connection of early adversities and harsh family environment with epigenetic ageing. We investigated whether early psychosocial resilience is associated with epigenetic ageing in adulthood. We used the population-based Young Finns data (n = 1593). Early psychosocial resilience was assessed in 1980-1989 across five broad domains: (1) index of psychological strength (self-esteem at home/in general/at school, perceived possibilities to influence at home, internal life control), (2) index of social satisfaction (perceived support from family/friends and life satisfaction), (3) index of leisure time activities (hobbies and physical fitness), (4) index of responsible health behaviors (infrequent smoking or alcohol consumption), and (5) index of school career (school grades and adaptation). Epigenetic ages were calculated for blood samples from 2011, and the analyses were performed with variables describing age deviation (AgeDevHannum, AgeDevHorvath, AgeDevPheno, AgeDevGrim) and DunedinPACE. Covariates included early family environment, polygenic risk scores for schizophrenia and major depression, adulthood education, and adulthood health behaviors. All of the early resilience indexes were associated with lower levels of epigenetic ageing in adulthood, most consistently with AgeDevGrim and DunedinPACE. The associations of psychological strength and social satisfaction, in particular, seemed to be non-linear. In a smaller subsample (n = 289), high early resilience was related to lower AgeDevGrim over a 25-year follow-up in those who had high "baseline" levels of AgeDevGrim. In conclusion, early resilience seems to associate with lower level of epigenetic ageing in adulthood. Our results tentatively suggest that early resilience may increase "equality in epigenetic ageing" in a general population.

Gene Set Based Integrated Methylome and Transcriptome Analysis Reveals Potential Molecular Mechanisms Linking Cigarette Smoking and Related Diseases.

Advanced integrative analysis of DNA methylation and transcriptomics data may provide deeper insights into smoke-induced epigenetic alterations, their effects on gene expression and related biological processes, linking cigarette smoking and related diseases. We hypothesize that accumulation of DNA methylation changes in CpG sites across genomic locations of different genes might have biological significance. We tested the hypothesis by performing gene set based integrative analysis of blood DNA methylation and transcriptomics data to identify potential transcriptomic consequences of smoking via changes in DNA methylation in the Young Finns Study (YFS) participants (n = 1114, aged 34-49 years, women: 54%, men: 46%). First, we performed epigenome-wide association study (EWAS) of smoking. We then defined sets of genes based on DNA methylation status within their genomic regions, for example, sets of genes containing hyper- or hypomethylated CpG sites in their body or promoter regions. Gene set analysis was performed using transcriptomics data from the same participants. Two sets of genes, one containing 49 genes with hypomethylated CpG sites in their body region and the other containing 33 genes with hypomethylated CpG sites in their promoter region, were differentially expressed among the smokers. Genes in the two gene sets are involved in bone formation, metal ion transport, cell death, peptidyl-serine phosphorylation, and cerebral cortex development process, revealing epigenetic-transcriptomic pathways to smoking-related diseases such as osteoporosis, atherosclerosis, and cognitive impairment. These findings contribute to a deeper understanding of the pathophysiology of smoking-related diseases and may provide potential therapeutic targets.

High genetic risk for depression as an independent risk factor for mortality in patients referred for coronary angiography.

Background: Different observations have suggested that patients with depression have a higher risk for a number of comorbidities and mortality. The underlying causes have not been fully understood yet. Aims: The aim of our study was to investigate the association of a genetic depression risk score (GDRS) with mortality [all-cause and cardiovascular (CV)] and markers of depression (including intake of antidepressants and a history of depression) in the Ludwigshafen Risk and Cardiovascular Health (LURIC) study involving 3,316 patients who had been referred for coronary angiography. Methods and results: The GDRS was calculated in 3,061 LURIC participants according to a previously published method and was found to be associated with all-cause (p = 0.016) and CV mortality (p = 0.0023). In Cox regression models adjusted for age, sex, body mass index, LDL-cholesterol, HDL-cholesterol, triglycerides, hypertension, smoking, and diabetes mellitus, the GDRS remained significantly associated with all-cause [1.18 (1.04-1.34, p = 0.013)] and CV [1.31 (1.11-1.55, p = 0.001)] mortality. The GDRS was not associated with the intake of antidepressants or a history of depression. However, this cohort of CV patients had not specifically been assessed for depression, leading to marked underreporting. We were unable to identify any specific biomarkers correlated with the GDRS in LURIC participants. Conclusion: A genetic predisposition for depression estimated by a GDRS was independently associated with all-cause and CV mortality in our cohort of patients who had been referred for coronary angiography. No biomarker correlating with the GDRS could be identified.

Uncovering the complex genetic architecture of human plasma lipidome using machine learning methods.

Genetic architecture of plasma lipidome provides insights into regulation of lipid metabolism and related diseases. We applied an unsupervised machine learning method, PGMRA, to discover phenotype-genotype many-to-many relations between genotype and plasma lipidome (phenotype) in order to identify the genetic architecture of plasma lipidome profiled from 1,426 Finnish individuals aged 30-45 years. PGMRA involves biclustering genotype and lipidome data independently followed by their inter-domain integration based on hypergeometric tests of the number of shared individuals. Pathway enrichment analysis was performed on the SNP sets to identify their associated biological processes. We identified 93 statistically significant (hypergeometric p-value < 0.01) lipidome-genotype relations. Genotype biclusters in these 93 relations contained 5977 SNPs across 3164 genes. Twenty nine of the 93 relations contained genotype biclusters with more than 50% unique SNPs and participants, thus representing most distinct subgroups. We identified 30 significantly enriched biological processes among the SNPs involved in 21 of these 29 most distinct genotype-lipidome subgroups through which the identified genetic variants can influence and regulate plasma lipid related metabolism and profiles. This study identified 29 distinct genotype-lipidome subgroups in the studied Finnish population that may have distinct disease trajectories and therefore could be useful in precision medicine research.

The relationship of trait-like compassion with epigenetic aging: The population-based prospective Young Finns Study.

Introduction: Helping others within and beyond the family has been related to living a healthy and long life. Compassion is a prosocial personality trait characterized by concern for another person who is suffering and the motivation to help. The current study examines whether epigenetic aging is a potential biological mechanism that explains the link between prosociality and longevity. Methods: We used data from the Young Finns Study that follows six birth-cohorts from age 3-18 to 19-49. Trait-like compassion for others was measured with the Temperament and Character Inventory in the years 1997 and 2001. Epigenetic age acceleration and telomere length were measured with five DNA methylation (DNAm) indicators (DNAmAgeHorvath, IEAA_Hannum, EEAA_Hannum, DNAmPhenoAge, and DNAmTL) based on blood drawn in 2011. We controlled for sex, socioeconomic status in childhood and adulthood, and body-mass index. Results and discussion: An association between higher compassion in 1997 and a less accelerated DNAmPhenoAge, which builds on previous work on phenotypic aging, approached statistical significance in a sex-adjusted model (n = 1,030; b = -0.34; p = 0.050). Compassion in 1997 predicted less accelerated epigenetic aging over and above the control variables (n = 843; b = -0.47; p = 0.016). There was no relationship between compassion in 2001 (n = 1108/910) and any of the other four studied epigenetic aging indicators. High compassion for others might indeed influence whether an individual's biological age is lower than their chronological age. The conducted robustness checks partially support this conclusion, yet cannot rule out that there might be a broader prosocial trait behind the findings. The observed associations are interesting but should be interpreted as weak requiring replication.

Dietary and supplemental intake of vitamins C and E is associated with altered DNA methylation in an epigenome-wide association study meta-analysis.

BACKGROUND: Dietary intake of antioxidants such as vitamins C and E protect against oxidative stress, and may also be associated with altered DNA methylation patterns. METHODS: We meta-analysed epigenome-wide association study (EWAS) results from 11,866 participants across eight population-based cohorts to evaluate the association between self-reported dietary and supplemental intake of vitamins C and E with DNA methylation. EWAS were adjusted for age, sex, BMI, caloric intake, blood cell type proportion, smoking status, alcohol consumption, and technical covariates. Significant results of the meta-analysis were subsequently evaluated in gene set enrichment analysis (GSEA) and expression quantitative trait methylation (eQTM) analysis. RESULTS: In meta-analysis, methylation at 4,656 CpG sites was significantly associated with vitamin C intake at FDR ≤ 0.05. The most significant CpG sites associated with vitamin C (at FDR ≤ 0.01) were enriched for pathways associated with systems development and cell signalling in GSEA, and were associated with downstream expression of genes enriched in the immune response in eQTM analysis. Furthermore, methylation at 160 CpG sites was significantly associated with vitamin E intake at FDR ≤ 0.05, but GSEA and eQTM analysis of the top most significant CpG sites associated with vitamin E did not identify significant enrichment of any biological pathways investigated. CONCLUSIONS: We identified significant associations of many CpG sites with vitamin C and E intake, and our results suggest that vitamin C intake may be associated with systems development and the immune response.

Methylation pattern of polymorphically imprinted nc886 is not conserved across mammalia.

BACKGROUND: In humans, the nc886 locus is a polymorphically imprinted metastable epiallele. Periconceptional conditions have an effect on the methylation status of nc886, and further, this methylation status is associated with health outcomes in later life, in line with the Developmental Origins of Health and Disease (DOHaD) hypothesis. Animal models would offer opportunities to study the associations between periconceptional conditions, nc886 methylation status and metabolic phenotypes further. Thus, we set out to investigate the methylation pattern of the nc886 locus in non-human mammals. DATA: We obtained DNA methylation data from the data repository GEO for mammals, whose nc886 gene included all three major parts of nc886 and had sequency similarity of over 80% with the human nc886. Our final sample set consisted of DNA methylation data from humans, chimpanzees, bonobos, gorillas, orangutangs, baboons, macaques, vervets, marmosets and guinea pigs. RESULTS: In human data sets the methylation pattern of nc886 locus followed the expected bimodal distribution, indicative of polymorphic imprinting. In great apes, we identified a unimodal DNA methylation pattern with 50% methylation level in all individuals and in all subspecies. In Old World monkeys, the between individual variation was greater and methylation on average was close to 60%. In guinea pigs the region around the nc886 homologue was non-methylated. Results obtained from the sequence comparison of the CTCF binding sites flanking the nc886 gene support the results on the DNA methylation data. CONCLUSIONS: Our results indicate that unlike in humans, nc886 is not a polymorphically imprinted metastable epiallele in non-human primates or in guinea pigs, thus implying that animal models are not applicable for nc886 research. The obtained data suggests that the nc886 region may be classically imprinted in great apes, and potentially also in Old World monkeys, but not in guinea pigs.

Multi-Omics Integration in a Twin Cohort and Predictive Modeling of Blood Pressure Values.

Abnormal blood pressure is strongly associated with risk of high-prevalence diseases, making the study of blood pressure a major public health challenge. Although biological mechanisms underlying hypertension at the single omic level have been discovered, multi-omics integrative analyses using continuous variations in blood pressure values remain limited. We used a multi-omics regression-based method, called sparse multi-block partial least square, for integrative, explanatory, and predictive interests in study of systolic and diastolic blood pressure values. Various datasets were obtained from the Finnish Twin Cohort for up to 444 twins. Blocks of omics-including transcriptomic, methylation, metabolomic-data as well as polygenic risk scores and clinical data were integrated into the modeling and supported by cross-validation. The predictive contribution of each omics block when predicting blood pressure values was investigated using external participants from the Young Finns Study. In addition to revealing interesting inter-omics associations, we found that each block of omics heterogeneously improved the predictions of blood pressure values once the multi-omics data were integrated. The modeling revealed a plurality of clinical, transcriptomic, and metabolomic factors consistent with the literature and that play a leading role in explaining unit variations in blood pressure. These findings demonstrate (1) the robustness of our integrative method to harness results obtained by single omics discriminant analyses, and (2) the added value of predictive and exploratory gains of a multi-omics approach in studies of complex phenotypes such as blood pressure.