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.

Polymorphisms in ERBB4 and TACR1 associated with dry mouth in clozapine-treated patients.

BACKGROUND: Sialorrhea is a common and uncomfortable adverse effect of clozapine, and its severity varies between patients. The aim of the study was to select broadly genes related to the regulation of salivation and study associations between sialorrhea and dry mouth and polymorphisms in the selected genes. METHODS: The study population consists of 237 clozapine-treated patients, of which 172 were genotyped. Associations between sialorrhea and dry mouth with age, sex, BMI, smoking, clozapine dose, clozapine and norclozapine serum levels, and other comedication were studied. Genetic associations were analyzed with linear and logistic regression models explaining sialorrhea and dry mouth with each SNP added separately to the model as coefficients. RESULTS: Clozapine dose, clozapine or norclozapine concentration and their ratio were not associated with sialorrhea or dryness of mouth. Valproate use (p = 0.013) and use of other antipsychotics (p = 0.015) combined with clozapine were associated with excessive salivation. No associations were found between studied polymorphisms and sialorrhea. In analyses explaining dry mouth with logistic regression with age and sex as coefficients, two proxy-SNPs were associated with dry mouth: epidermal growth factor receptor 4 (ERBB4) rs3942465 (adjusted p = 0.025) and tachykinin receptor 1 (TACR1) rs58933792 (adjusted p = 0.029). CONCLUSION: Use of valproate or antipsychotic polypharmacy may increase the risk of sialorrhea. Genetic variations in ERBB4 and TACR1 might contribute to experienced dryness of mouth among patients treated with clozapine.

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.

Genetic risk scores associated with temperament clusters in Finnish depression patients.

OBJECTIVE: Cloninger's temperament dimensions have been studied widely in relation to genetics. In this study, we examined Cloninger's temperament dimensions grouped with cluster analyses and their association with single nucleotide polymorphisms (SNPs). This study included 212 genotyped Finnish patients from the Ostrobothnia Depression Study. METHODS: The temperament clusters were analysed at baseline and at six weeks from the beginning of the depression intervention study. We selected depression-related catecholamine and serotonin genes based on a literature search, and 59 SNPs from ten different genes were analysed. The associations of single SNPs with temperament clusters were studied. Using the selected genes, genetic risk score (GRS) analyses were conducted considering appropriate confounding factors. RESULTS: No single SNP had a significant association with the temperament clusters. Associations between GRSs and temperament clusters were observed in multivariate models that were significant after permutation analyses. Two SNPs from the DRD3 gene, two SNPs from the SLC6A2 gene, one SNP from the SLC6A4 gene, and one SNP from the HTR2A gene associated with the HHA/LRD/LP (high harm avoidance, low reward dependence, low persistence) cluster at baseline. Two SNPs from the HTR2A gene were associated with the HHA/LRD/LP cluster at six weeks. Two SNPs from the HTR2A gene and two SNPs from the COMT gene were associated with the HP (high persistence) cluster at six weeks. CONCLUSION: GRSs seem to associate with an individual's temperament profile, which can be observed in the clusters used. Further research needs to be conducted on these types of clusters and their clinical applicability.

Genetic insights into resting heart rate and its role in cardiovascular disease.

Resting heart rate is associated with cardiovascular diseases and mortality in observational and Mendelian randomization studies. The aims of this study are to extend the number of resting heart rate associated genetic variants and to obtain further insights in resting heart rate biology and its clinical consequences. A genome-wide meta-analysis of 100 studies in up to 835,465 individuals reveals 493 independent genetic variants in 352 loci, including 68 genetic variants outside previously identified resting heart rate associated loci. We prioritize 670 genes and in silico annotations point to their enrichment in cardiomyocytes and provide insights in their ECG signature. Two-sample Mendelian randomization analyses indicate that higher genetically predicted resting heart rate increases risk of dilated cardiomyopathy, but decreases risk of developing atrial fibrillation, ischemic stroke, and cardio-embolic stroke. We do not find evidence for a linear or non-linear genetic association between resting heart rate and all-cause mortality in contrast to our previous Mendelian randomization study. Systematic alteration of key differences between the current and previous Mendelian randomization study indicates that the most likely cause of the discrepancy between these studies arises from false positive findings in previous one-sample MR analyses caused by weak-instrument bias at lower P-value thresholds. The results extend our understanding of resting heart rate biology and give additional insights in its role in cardiovascular disease development.

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.

Gene set analysis of transcriptomics data identifies new biological processes associated with early markers of atherosclerosis but not with those of osteoporosis: Atherosclerosis-osteoporosis co/multimorbidity study in the Young Finns Study.

AIM: We aimed at identifying the shared biological processes underlying atherosclerosis-osteoporosis co/multimorbidity. METHODS: We performed gene set analysis (GSA) of whole-blood transcriptomic data to identify biological processes shared by the early markers of these two diseases. Early markers of diseases, carotid intima-media thickness (CIMT) for atherosclerosis and trabecular bone mineral density (BMD) from distal radius and tibia for osteoporosis, were used to categorize the study participants into cases and controls. Participants with high CIMT (>90th percentile) were defined as cases for subclinical atherosclerosis. Study population-based T-scores for BMD were calculated and T-score ≤ -1 was used for the definition of low BMD cases i.e., early indicator of osteoporosis. RESULTS: We did not identify any gene sets jointly associated with early markers of atherosclerosis and osteoporosis. We identified three novel and replicated 234 gene sets significantly associated with high CIMT with false discovery rate (FDR) ≤ 0.01. Only two genes, both related to the immune system, were identified to be associated with high CIMT by traditional differential gene expression analysis. However, none of the studied gene sets or individual genes were significantly associated with tibial or radial BMD. The three novel CIMT associated gene sets contained genes involved in copper homeostasis, neural crest cell migration and nicotinate and nicotinamide metabolism. The 234 replicated gene sets in this study are related to the immune system, hypoxia and apoptosis, consistent with the existing literature on atherosclerosis. CONCLUSIONS: This study identified novel biological processes associated with high CIMT but not with reduced BMD.

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.

Methylation status of nc886 epiallele reflects periconceptional conditions and is associated with glucose metabolism through nc886 RNAs.

BACKGROUND: Non-coding RNA 886 (nc886) is coded from a maternally inherited metastable epiallele. We set out to investigate the determinants and dynamics of the methylation pattern at the nc886 epiallele and how this methylation status associates with nc886 RNA expression. Furthermore, we investigated the associations between the nc886 methylation status or the levels of nc886 RNAs and metabolic traits in the YFS and KORA cohorts. The association between nc886 epiallele methylation and RNA expression was also validated in induced pluripotent stem cell (iPSC) lines. RESULTS: We confirm that the methylation status of the nc886 epiallele is mostly binomial, with individuals displaying either a non- or hemi-methylated status, but we also describe intermediately and close to fully methylated individuals. We show that an individual's methylation status is associated with the mother's age and socioeconomic status, but not with the individual's own genetics. Once established, the methylation status of the nc886 epiallele remains stable for at least 25 years. This methylation status is strongly associated with the levels of nc886 non-coding RNAs in serum, blood, and iPSC lines. In addition, nc886 methylation status associates with glucose and insulin levels during adolescence but not with the indicators of glucose metabolism or the incidence of type 2 diabetes in adulthood. However, the nc886-3p RNA levels also associate with glucose metabolism in adulthood. CONCLUSIONS: These results indicate that nc886 metastable epiallele methylation is tuned by the periconceptional conditions and it associates with glucose metabolism through the expression of the ncRNAs coded in the epiallele region.

Human Prostate Tissue MicroRNAs and Their Predicted Target Pathways Linked to Prostate Cancer Risk Factors.

MicroRNAs are important in prostate cancer development, progression and metastasis. The aim of this study was to test microRNA expression profile in prostate tissue obtained from prostate cancer patients for associations with various prostate cancer related factors and to pinpoint the predicted target pathways for these microRNAs. Prostate tissue samples were obtained at prostatectomy from patients participating in a trial evaluating impact of pre-operative atorvastatin on serum prostate specific antigen (PSA) and Ki-67 expression in prostate tissue. Prostate tissue microRNA expression profiles were analyzed using OpenArray® MicroRNA Panel. Pathway enrichment analyses were conducted for predicted target genes of microRNAs that correlated significantly with studied factors. Eight microRNAs correlated significantly with studied factors of patients after Bonferroni multiple testing correction. MiR-485-3p correlated with serum HDL-cholesterol levels. In atorvastatin-treated subjects, miR-34c-5p correlated with a change in serum PSA and miR-138-3p with a change in total cholesterol. In the placebo arm, both miR-576-3p and miR-550-3p correlated with HDL-cholesterol and miR-627 with PSA. In pathway analysis, these eight microRNAs related significantly to several pathways relevant to prostate cancer. This study brings new evidence from the expression of prostate tissue microRNAs and related pathways that may link risk factors to prostate cancer and pinpoint new therapeutic possibilities.

Uncovering the shared lipidomic markers of subclinical osteoporosis-atherosclerosis comorbidity: The Young Finns Study.

BACKGROUND: Osteoporosis and atherosclerosis are complex multifactorial diseases sharing common risk factors and pathophysiological mechanisms suggesting that these are comorbidities. Omics studies identifying joint molecular markers associated with these diseases are sparse. SUBJECTS AND METHODS: Using liquid chromatography-tandem mass spectrometry, we quantified 437 molecular lipid species from the Young Finns Study cohort (aged 30-45 years and 57% women) and performed lipidome-wide multivariate analysis of variance (MANOVA) with early markers for both diseases. Carotid intima-media thickness for atherosclerosis measured with ultrasound and bone mineral density from distal radius and tibia for osteoporosis measured with peripheral quantitative computed tomography were used as early markers of the diseases. RESULTS: MANOVA adjusted with age, sex and body mass index, identified eight statistically significant (adjusted p-value (padj) < 0.05) and 15 suggestively significant (padj < 0.25) molecular lipid species associated with the studied markers. Similar analysis adjusted additionally for smoking habit, physical activity and alcohol consumption identified four significant and six suggestively significant molecular lipid species. These most significant lipid classes/species jointly associated with the studied markers were glycerolipid/TAG(18:0/18:0/18:1), glycerophospholipid/PC(40:3), sphingolipid/Gb3(d18:1/22:0), and sphingolipid/Gb3(d18:1/24:0). CONCLUSION: Our results support the osteoporosis-atherosclerosis comorbidity hypothesis and present potential new joint lipid biomarkers for these diseases.

Modular genome-wide gene expression architecture shared by early traits of osteoporosis and atherosclerosis in the Young Finns Study.

We analysed whole blood genome-wide expression data to identify gene co-expression modules shared by early traits of osteoporosis and atherosclerosis. Gene expression was profiled for the Young Finns Study participants. Bone mineral density and content were measured as early traits of osteoporosis. Carotid and bulbus intima media thickness were measured as early traits of atherosclerosis. Joint association of the modules, identified with weighted co-expression analysis, with early traits of the diseases was tested with multivariate analysis. Among the six modules significantly correlated with early traits of both the diseases, two had significant (adjusted p-values (p.adj) < 0.05) and another two had suggestively significant (p.adj < 0.25) joint association with the two diseases after adjusting for age, sex, body mass index, smoking habit, alcohol consumption, and physical activity. The three most significant member genes from the significant modules were NOSIP, GXYLT2, and TRIM63 (p.adj ≤ 0.18). Genes in the modules were enriched with biological processes that have separately been found to be involved in either bone metabolism or atherosclerosis. The gene modules and their most significant member genes identified in this study support the osteoporosis-atherosclerosis comorbidity hypothesis and can provide new joint biomarkers for both diseases and their dual prevention.

Examining the effect of mitochondrial DNA variants on blood pressure in two Finnish cohorts.

High blood pressure (BP) is a major risk factor for many noncommunicable diseases. The effect of mitochondrial DNA single-nucleotide polymorphisms (mtSNPs) on BP is less known than that of nuclear SNPs. We investigated the mitochondrial genetic determinants of systolic, diastolic, and mean arterial BP. MtSNPs were determined from peripheral blood by sequencing or with genome-wide association study SNP arrays in two independent Finnish cohorts, the Young Finns Study and the Finnish Cardiovascular Study, respectively. In total, over 4200 individuals were included. The effects of individual common mtSNPs, with an additional focus on sex-specificity, and aggregates of rare mtSNPs grouped by mitochondrial genes were evaluated by meta-analysis of linear regression and a sequence kernel association test, respectively. We accounted for the predicted pathogenicity of the rare variants within protein-encoding and the tRNA regions. In the meta-analysis of 87 common mtSNPs, we did not observe significant associations with any of the BP traits. Sex-specific and rare-variant analyses did not pinpoint any significant associations either. Our results are in agreement with several previous studies suggesting that mtDNA variation does not have a significant role in the regulation of BP. Future studies might need to reconsider the mechanisms thought to link mtDNA with hypertension.

Sex-specific associations of TCF7L2 variants with fasting glucose, type 2 diabetes and coronary heart disease among Turkish adults.

OBJECTIVE: TCF7L2 is a repressor and transactivator of genes, and its variants are strongly associated with diabetes. This study aimed to evaluate the sex-specific relationship between the most common TCF7L2 gene variants (-98368G>T, rs12255372 and -47833C>T, rs7903146) with diabetes and coronary heart disease in Turkish Adult Risk Factor (TARF) Study. METHODS: Single nucleotide variants (SNVs) have been genotyped using the TaqMan allelic discrimination assays in 2,024 (51.3% in women, age: 55±11.8) Turkish adults participating in the TARF study. Statistical analyses were used to investigate the association of genotypes with clinical and biochemical measurements. RESULTS: Among the TARF study participants, 11.7%, 24.3%, 14.1%, and 38.3% had diabetes, hypertension, coronary heart disease (CHD), and obesity, respectively. The frequencies of T allele for -47833C>T and -98368G>T in Turkish adults were determined to be 0.35 and 0.33, respectively. -47833C>T was significantly associated with higher fasting glucose concentrations in all participants, especially in men. Both SNVs were significantly associated with diabetes and CHD in all participants (p<0.05). When study population was stratified according to sex, -98368G>T was associated with diabetes in women (p=0.041) and -47833C>T was associated with diabetes and CHD in men (p=0.018 and p=0.032, respectively). Also, both SNVs and the diplotypes of common haplotype (H1) remained strongly associated with type 2 diabetes after risk factors were adjusted (p<0.05). CONCLUSION: T allele homozygosity of two SNVs as well as the diplotype H1-/H1- reflects risk of diabetes primarily in men. Enhanced CHD risk is determined by the presence of diplotype H1-/H1- among nondiabetic participants.

Identification, Heritability, and Relation With Gene Expression of Novel DNA Methylation Loci for Blood Pressure.

We conducted an epigenome-wide association study meta-analysis on blood pressure (BP) in 4820 individuals of European and African ancestry aged 14 to 69. Genome-wide DNA methylation data from peripheral leukocytes were obtained using the Infinium Human Methylation 450k BeadChip. The epigenome-wide association study meta-analysis identified 39 BP-related CpG sites with P<1×10-5. In silico replication in the CHARGE consortium of 17 010 individuals validated 16 of these CpG sites. Out of the 16 CpG sites, 13 showed novel association with BP. Conversely, out of the 126 CpG sites identified as being associated (P<1×10-7) with BP in the CHARGE consortium, 21 were replicated in the current study. Methylation levels of all the 34 CpG sites that were cross-validated by the current study and the CHARGE consortium were heritable and 6 showed association with gene expression. Furthermore, 9 CpG sites also showed association with BP with P<0.05 and consistent direction of the effect in the meta-analysis of the Finnish Twin Cohort (199 twin pairs and 4 singletons; 61% monozygous) and the Netherlands Twin Register (266 twin pairs and 62 singletons; 84% monozygous). Bivariate quantitative genetic modeling of the twin data showed that a majority of the phenotypic correlations between methylation levels of these CpG sites and BP could be explained by shared unique environmental rather than genetic factors, with 100% of the correlations of systolic BP with cg19693031 (TXNIP) and cg00716257 (JDP2) determined by environmental effects acting on both systolic BP and methylation levels.

Epigenome-450K-wide methylation signatures of active cigarette smoking: The Young Finns Study.

Smoking as a major risk factor for morbidity affects numerous regulatory systems of the human body including DNA methylation. Most of the previous studies with genome-wide methylation data are based on conventional association analysis and earliest threshold-based gene set analysis that lacks sensitivity to be able to reveal all the relevant effects of smoking. The aim of the present study was to investigate the impact of active smoking on DNA methylation at three biological levels: 5'-C-phosphate-G-3' (CpG) sites, genes and functionally related genes (gene sets). Gene set analysis was done with mGSZ, a modern threshold-free method previously developed by us that utilizes all the genes in the experiment and their differential methylation scores. Application of such method in DNA methylation study is novel. Epigenome-wide methylation levels were profiled from Young Finns Study (YFS) participants' whole blood from 2011 follow-up using Illumina Infinium HumanMethylation450 BeadChips. We identified three novel smoking related CpG sites and replicated 57 of the previously identified ones. We found that smoking is associated with hypomethylation in shore (genomic regions 0-2 kilobases from CpG island). We identified smoking related methylation changes in 13 gene sets with false discovery rate (FDR) ≤ 0.05, among which is olfactory receptor activity, the flagship novel finding of the present study. Overall, we extended the current knowledge by identifying: (i) three novel smoking related CpG sites, (ii) similar effects as aging on average methylation in shore, and (iii) a novel finding that olfactory receptor activity pathway responds to tobacco smoke and toxin exposure through epigenetic mechanisms.

Multi-ancestry GWAS of the electrocardiographic PR interval identifies 202 loci underlying cardiac conduction.

The electrocardiographic PR interval reflects atrioventricular conduction, and is associated with conduction abnormalities, pacemaker implantation, atrial fibrillation (AF), and cardiovascular mortality. Here we report a multi-ancestry (N = 293,051) genome-wide association meta-analysis for the PR interval, discovering 202 loci of which 141 have not previously been reported. Variants at identified loci increase the percentage of heritability explained, from 33.5% to 62.6%. We observe enrichment for cardiac muscle developmental/contractile and cytoskeletal genes, highlighting key regulation processes for atrioventricular conduction. Additionally, 8 loci not previously reported harbor genes underlying inherited arrhythmic syndromes and/or cardiomyopathies suggesting a role for these genes in cardiovascular pathology in the general population. We show that polygenic predisposition to PR interval duration is an endophenotype for cardiovascular disease, including distal conduction disease, AF, and atrioventricular pre-excitation. These findings advance our understanding of the polygenic basis of cardiac conduction, and the genetic relationship between PR interval duration and cardiovascular disease.

Lipidomic architecture shared by subclinical markers of osteoporosis and atherosclerosis: The Cardiovascular Risk in Young Finns Study.

BACKGROUND: Studies have shown that osteoporosis and atherosclerosis are comorbid conditions sharing common risk factors and pathophysiological mechanisms. Understanding these is crucial in order to develop shared methods for risk stratification, prevention, diagnosis and treatment. The aim of this study was to apply a system-level bioinformatics approach to lipidome-wide data in order to pinpoint the lipidomic architecture jointly associated with surrogate markers of these complex comorbid diseases. SUBJECTS AND METHODS: The study was based on the Cardiovascular Risk in Young Finns Study cohort from the 2007 follow-up (n = 1494, aged 30-45 years, women: 57%). Liquid chromatography-tandem mass spectrometry (LC-MS/MS) was used to analyse the serum lipidome, involving 437 molecular lipid species. The subclinical osteoporotic markers included indices of bone mineral density and content, measured using peripheral quantitative computer tomography from the distal and shaft sites of both the tibia and the radius. The subclinical atherosclerotic markers included carotid and bulbus intima media thickness measured with high-resolution ultrasound. Weighted co-expression network analysis was performed to identify networks of densely interconnected lipid species (i.e. lipid modules) associated with subclinical markers of both osteoporosis and atherosclerosis. The levels of lipid species (lipid profiles) of each of the lipid modules were summarized by the first principal component termed as module eigenlipid. Then, Pearson's correlation (r) was calculated between the module eigenlipids and the markers. Lipid modules that were significantly and jointly correlated with subclinical markers of both osteoporosis and atherosclerosis were considered to be related to the comorbidities. The hypothesis that the eigenlipids and profiles of the constituent lipid species in the modules have joint effects on the markers was tested with multivariate analysis of variance (MANOVA). RESULTS: Among twelve studied molecular lipid modules, we identified one module with 105 lipid species significantly and jointly associated with both subclinical markers of both osteoporosis (r = 0.24, p-value = 2 × 10-20) and atherosclerosis (r = 0.16, p-value = 2 × 10-10). The majority of the lipid species in this module belonged to the glycerolipid (n = 60), glycerophospholipid (n = 13) and sphingolipid (n = 29) classes. The module was also enriched with ceramides (n = 20), confirming their significance in cardiovascular outcomes and suggesting their joint role in the comorbidities. The top three of the 37 statistically significant (adjusted p-value < 0.05) lipid species jointly associated with subclinical markers of both osteoporosis and atherosclerosis within the module were all triacylglycerols (TAGs) - TAG(18:0/18:0/18:1) with an adjusted p-value of 8.6 × 10-8, TAG(18:0/18:1/18:1) with an adjusted p-value of 3.7 × 10-6, and TAG(16:0/18:0/18:1) with an adjusted p-value of 8.5 × 10-6. CONCLUSION: This study identified a novel lipid module associated with both surrogate markers of both subclinical osteoporosis and subclinical atherosclerosis. Alterations in the metabolism of the identified lipid module and, more specifically, the TAG related molecular lipids within the module may provide potential new biomarkers for testing the comorbidities, opening avenues for the emergence of dual-purpose prevention measures.

New evidence from plasma ceramides links apoE polymorphism to greater risk of coronary artery disease in Finnish adults.

apoE, a key regulator of plasma lipids, mediates altered functionalities in lipoprotein metabolism and thus affects the risk of coronary artery disease (CAD). The significance of different apoE polymorphisms remains unclear; although the ε4 allele is clearly associated with increased cholesterol levels (which inform CAD risk), direct studies about apoE polymorphisms on CAD risk and development have yielded controversial results. Furthermore, certain species of ceramides-complex lipids abundant in plasma LDL-are markers of increased risk of myocardial infarction and cardiovascular death. Using a high-throughput MS approach, we quantified 30 molecular plasma ceramide species from a cohort of 2,160 apoE-genotyped (rs7412, rs429358) young adults enrolled in the population-based Cardiovascular Risk in Young Finns Study. We then searched this lipidome data set to identify new indications of pathways influenced by apoE polymorphisms and possibly related to CAD risk. This approach revealed a previously unreported association between apoE polymorphism and a consistently documented high-risk CAD marker, Cer(d18:1/16:0). Compared with the apoE ε3/3 reference group, plasma levels of apoE ε4 were elevated and those of apoE ε2 were lowered in all subjects without evidence of apoE-by-sex interactions. apoE associated with seven ceramides that are connected to atherogenically potent macrophages and/or lipoprotein particles; these associations could indicate a plausible linkage between apoE polymorphism and ceramide metabolism, leading to adverse plasma LDL metabolism and atherogenesis. In conclusion, new evidence from plasma ceramides links apoE polymorphism with an increased risk of CAD and extends our understanding of the role of apoE in health and disease.

Whole blood microRNA levels associate with glycemic status and correlate with target mRNAs in pathways important to type 2 diabetes.

We analyzed the associations between whole blood microRNA profiles and the indices of glucose metabolism and impaired fasting glucose and examined whether the discovered microRNAs correlate with the expression of their mRNA targets. MicroRNA and gene expression profiling were performed for the Young Finns Study participants (n = 871). Glucose, insulin, and glycated hemoglobin (HbA1c) levels were measured, the insulin resistance index (HOMA2-IR) was calculated, and the glycemic status (normoglycemic [n = 534]/impaired fasting glucose [IFG] [n = 252]/type 2 diabetes [T2D] [n = 24]) determined. Levels of hsa-miR-144-5p, -122-5p, -148a-3p, -589-5p, and hsa-let-7a-5p associated with glycemic status. hsa-miR-144-5p and -148a-3p associated with glucose levels, while hsa-miR-144-5p, -122-5p, -184, and -339-3p associated with insulin levels and HOMA2-IR, and hsa-miR-148a-3p, -15b-3p, -93-3p, -146b-5p, -221-3p, -18a-3p, -642a-5p, and -181-2-3p associated with HbA1c levels. The targets of hsa-miR-146b-5p that correlated with its levels were enriched in inflammatory pathways, and the targets of hsa-miR-221-3p were enriched in insulin signaling and T2D pathways. These pathways showed indications of co-regulation by HbA1c-associated miRNAs. There were significant differences in the microRNA profiles associated with glucose, insulin, or HOMA-IR compared to those associated with HbA1c. The HbA1c-associated miRNAs also correlated with the expression of target mRNAs in pathways important to the development of T2D.