Gene expression signatures of coronary heart disease

Enhancing selection of alcohol consumption-associated genes by random forest

Machine learning methods have been used in identifying omics markers for a variety of phenotypes. We aimed to examine whether a supervised machine learning algorithm can improve identification of alcohol-associated transcriptomic markers. In this study, we analysed array-based, whole-blood derived expression data for 17 873 gene transcripts in 5508 Framingham Heart Study participants. By using the Boruta algorithm, a supervised random forest (RF)-based feature selection method, we selected twenty-five alcohol-associated transcripts. In a testing set (30 % of entire study participants), AUC (area under the receiver operating characteristics curve) of these twenty-five transcripts were 0·73, 0·69 and 0·66 for non-drinkers v. moderate drinkers, non-drinkers v. heavy drinkers and moderate drinkers v. heavy drinkers, respectively. The AUC of the selected transcripts by the Boruta method were comparable to those identified using conventional linear regression models, for example, AUC of 1958 transcripts identified by conventional linear regression models (false discovery rate < 0·2) were 0·74, 0·66 and 0·65, respectively. With Bonferroni correction for the twenty-five Boruta method-selected transcripts and three CVD risk factors (i.e. at P < 6·7e-4), we observed thirteen transcripts were associated with obesity, three transcripts with type 2 diabetes and one transcript with hypertension. For example, we observed that alcohol consumption was inversely associated with the expression of DOCK4, IL4R, and SORT1, and DOCK4 and SORT1 were positively associated with obesity, and IL4R was inversely associated with hypertension. In conclusion, using a supervised machine learning method, the RF-based Boruta algorithm, we identified novel alcohol-associated gene transcripts.

Leveraging the transcriptome to further our understanding of GWAS findings: eQTLs associated with genes related to LDL and LDL subclasses, in a cohort of African Americans

Background: GWAS discoveries often pose a significant challenge in terms of understanding their underlying mechanisms. Further research, such as an integration with expression quantitative trait locus (eQTL) analyses, are required to decipher the mechanisms connecting GWAS variants to phenotypes. An eQTL analysis was conducted on genes associated with low-density lipoprotein (LDL) cholesterol and its subclasses, with the aim of pinpointing genetic variants previously implicated in GWAS studies focused on lipid-related traits. Notably, the study cohort consisted of African Americans, a population characterized by a heightened prevalence of hypercholesterolemia. Methods: A comprehensive differential expression (DE) analysis was undertaken, with a dataset of 17,948 protein-coding mRNA transcripts extracted from the whole-blood transcriptomes of 416 samples to identify mRNA transcripts associated with LDL, with further granularity delineated between small LDL and large LDL subclasses. Subsequently, eQTL analysis was conducted with a subset of 242 samples for which whole-genome sequencing data were available to identify single-nucleotide polymorphisms (SNPs) associated with the LDL-related mRNA transcripts. Lastly, plausible functional connections were established between the identified eQTLs and genetic variants reported in the GWAS catalogue. Results: DE analysis revealed 1,048, 284, and 94 mRNA transcripts that exhibited differential expression in response to LDL, small LDL, and large LDL, respectively. The eQTL analysis identified a total of 9,950 significant SNP-mRNA associations involving 6,955 SNPs including a subset 101 SNPs previously documented in GWAS of LDL and LDL-related traits. Conclusion: Through comprehensive differential expression analysis, we identified numerous mRNA transcripts responsive to LDL, small LDL, and large LDL. Subsequent eQTL analysis revealed a rich landscape of eQTL-mRNA associations, including a subset of eQTL reported in GWAS studies of LDL and related traits. The study serves as a testament to the important role of integrative genomics in unraveling the enigmatic GWAS relationships between genetic variants and the complex fabric of human traits and diseases.

JASPER: fast, powerful, multitrait association testing in structured samples gives insight on pleiotropy in gene expression

Joint association analysis of multiple traits with multiple genetic variants can provide insight into genetic architecture and pleiotropy, improve trait prediction and increase power for detecting association. Furthermore, some traits are naturally high-dimensional, e.g., images, networks or longitudinally measured traits. Assessing significance for multitrait genetic association can be challenging, especially when the sample has population sub-structure and/or related individuals. Failure to adequately adjust for sample structure can lead to power loss and inflated type 1 error, and commonly used methods for assessing significance can work poorly with a large number of traits or be computationally slow. We developed JASPER, a fast, powerful, robust method for assessing significance of multitrait association with a set of genetic variants, in samples that have population sub-structure, admixture and/or relatedness. In simulations, JASPER has higher power, better type 1 error control, and faster computation than existing methods, with the power and speed advantage of JASPER increasing with the number of traits. JASPER is potentially applicable to a wide range of association testing applications, including for multiple disease traits, expression traits, image-derived traits and microbiome abundances. It allows for covariates, ascertainment and rare variants and is robust to phenotype model misspecification. We apply JASPER to analyze gene expression in the Framingham Heart Study, where, compared to alternative approaches, JASPER finds more significant associations, including several that indicate pleiotropic effects, some of which replicate previous results, while others have not previously been reported. Our results demonstrate the promise of JASPER for powerful multitrait analysis in structured samples.

Integrated omics analysis of coronary artery calcifications and myocardial infarction: the Framingham Heart Study

Gene function can be described using various measures. We integrated association studies of three types of omics data to provide insights into the pathophysiology of subclinical coronary disease and myocardial infarction (MI). Using multivariable regression models, we associated: (1) single nucleotide polymorphism, (2) DNA methylation, and (3) gene expression with coronary artery calcification (CAC) scores and MI. Among 3106 participants of the Framingham Heart Study, 65 (2.1%) had prevalent MI and 60 (1.9%) had incident MI, median CAC value was 67.8 [IQR 10.8, 274.9], and 1403 (45.2%) had CAC scores > 0 (prevalent CAC). Prevalent CAC was associated with AHRR (linked to smoking) and EXOC3 (affecting platelet function and promoting hemostasis). CAC score was associated with VWA1 (extracellular matrix protein associated with cartilage structure in endomysium). For prevalent MI we identified FYTTD1 (down-regulated in familial hypercholesterolemia) and PINK1 (linked to cardiac tissue homeostasis and ischemia-reperfusion injury). Incident MI was associated with IRX3 (enhancing browning of white adipose tissue) and STXBP3 (controlling trafficking of glucose transporter type 4 to plasma). Using an integrative trans-omics approach, we identified both putatively novel and known candidate genes associated with CAC and MI. Replication of findings is warranted.

Expression quantitative trait methylation analysis elucidates gene regulatory effects of DNA methylation: the Framingham Heart Study

Expression quantitative trait methylation (eQTM) analysis identifies DNA CpG sites at which methylation is associated with gene expression. The present study describes an eQTM resource of CpG-transcript pairs derived from whole blood DNA methylation and RNA sequencing gene expression data in 2115 Framingham Heart Study participants. We identified 70,047 significant cis CpG-transcript pairs at p < 1E-7 where the top most significant eGenes (i.e., gene transcripts associated with a CpG) were enriched in biological pathways related to cell signaling, and for 1208 clinical traits (enrichment false discovery rate [FDR] ≤ 0.05). We also identified 246,667 significant trans CpG-transcript pairs at p < 1E-14 where the top most significant eGenes were enriched in biological pathways related to activation of the immune response, and for 1191 clinical traits (enrichment FDR ≤ 0.05). Independent and external replication of the top 1000 significant cis and trans CpG-transcript pairs was completed in the Women's Health Initiative and Jackson Heart Study cohorts. Using significant cis CpG-transcript pairs, we identified significant mediation of the association between CpG sites and cardiometabolic traits through gene expression and identified shared genetic regulation between CpGs and transcripts associated with cardiometabolic traits. In conclusion, we developed a robust and powerful resource of whole blood eQTM CpG-transcript pairs that can help inform future functional studies that seek to understand the molecular basis of disease.

Blood transcriptomic biomarkers of alcohol consumption and cardiovascular disease risk factors: the Framingham Heart Study

BACKGROUND

The relations of alcohol consumption and gene expression remain to be elucidated.

MATERIALS AND METHODS

We examined cross-sectional associations between alcohol consumption and whole blood derived gene expression levels and between alcohol-associated genes and obesity, hypertension, and diabetes in 5531 Framingham Heart Study (FHS) participants.

RESULTS

We identified 25 alcohol-associated genes. We further showed cross-sectional associations of 16 alcohol-associated genes with obesity, nine genes with hypertension, and eight genes with diabetes at P < 0.002. For example, we observed decreased expression of PROK2 (β = -0.0018; 95%CI: -0.0021, -0.0007; P = 6.5e - 5) and PAX5 (β = -0.0014; 95%CI: -0.0021, -0.0007; P = 6.5e - 5) per 1 g/day increase in alcohol consumption. Consistent with our previous observation on the inverse association of alcohol consumption with obesity and positive association of alcohol consumption with hypertension, we found that PROK2 was positively associated with obesity (OR = 1.42; 95%CI: 1.17, 1.72; P = 4.5e - 4) and PAX5 was negatively associated with hypertension (OR = 0.73; 95%CI: 0.59, 0.89; P = 1.6e - 3). We also observed that alcohol consumption was positively associated with expression of ABCA13 (β = 0.0012; 95%CI: 0.0007, 0.0017; P = 1.3e - 6) and ABCA13 was positively associated with diabetes (OR = 2.57; 95%CI: 1.73, 3.84; P = 3.5e - 06); this finding, however, was inconsistent with our observation of an inverse association between alcohol consumption and diabetes.

CONCLUSIONS

We showed strong cross-sectional associations between alcohol consumption and expression levels of 25 genes in FHS participants. Nonetheless, complex relationships exist between alcohol-associated genes and CVD risk factors.

Vascular smooth muscle cells in low SYNTAX scores coronary artery disease exhibit proinflammatory transcripts and proteins correlated with IL1B activation

BACKGROUND AND AIMS

The SYNTAX score is clinically validated to stratify number of lesions and pattern of CAD. A better understanding of the underlying molecular mechanisms influencing the pattern and complexity of coronary arteries lesions among CAD patients is needed.

METHODS

Human arterial biopsies from 49 patients (16 low-SYNTAX-score (LSS, <23), 16 intermediate-SYNTAX-score (ISS, 23 to 32) and 17 high-SYNTAX-score (HSS, >32)) were evaluated using Affymetrix GeneChip® Human Genome U133 Plus 2.0 microarray. The data were validated by Next-Generation Sequencing (NGS). Primary VSMC from patients with low and high SYNTAX scores were isolated and compared using immunohistochemistry, qPCR and immunoblotting to confirm mRNA and proteomic results.

RESULTS

The IL1B was verified as the top upstream regulator of 47 inflammatory DEGs in LSS patients and validated by another sets of patient samples using NGS analysis. The upregulated expression of IL1B was translated to increased level of IL1β protein in the LSS tissue based on immunohistochemical quantitative analysis. Plausibility of idea that IL1B in the arterial wall could be originated from VSMC was checked by exposing culture to proinflammatory conditions where IL1B came out as the top DEG (logFC = 7.083, FDR = 1.38 × 10-114). The LSS patient-derived primary VSMCs confirmed higher levels of IL1B mRNA and protein.

CONCLUSIONS

LSS patients could represent a group of patients where IL1B could play a substantial role in disease pathogenesis. The LSS group could represent a plausible cohort of patients for whom anti-inflammatory therapy could be considered.

The impact of alcoholic drinks and dietary factors on epigenetic markers associated with triglyceride levels

Background: Many epigenetic loci have been associated with plasma triglyceride (TG) levels, but epigenetic connections between those loci and dietary exposures are largely unknown. This study aimed to characterize the epigenetic links between diet, lifestyle, and TG. Methods: We first conducted an epigenome-wide association study (EWAS) for TG in the Framingham Heart Study Offspring population (FHS, n = 2,264). We then examined relationships between dietary and lifestyle-related variables, collected four times in 13 years, and differential DNA methylation sites (DMSs) associated with the last TG measures. Third, we conducted a mediation analysis to evaluate the causal relationships between diet-related variables and TG. Finally, we replicated three steps to validate identified DMSs associated with alcohol and carbohydrate intake in the Genetics of Lipid-Lowering Drugs and Diet Network (GOLDN) study (n = 993). Results: In the FHS, the EWAS revealed 28 TG-associated DMSs at 19 gene regions. We identified 102 unique associations between these DMSs and one or more dietary and lifestyle-related variables. Alcohol and carbohydrate intake showed the most significant and consistent associations with 11 TG-associated DMSs. Mediation analyses demonstrated that alcohol and carbohydrate intake independently affect TG via DMSs as mediators. Higher alcohol intake was associated with lower methylation at seven DMSs and higher TG. In contrast, increased carbohydrate intake was associated with higher DNA methylation at two DMSs (CPT1A and SLC7A11) and lower TG. Validation in the GOLDN further supports the findings. Conclusion: Our findings imply that TG-associated DMSs reflect dietary intakes, particularly alcoholic drinks, which could affect the current cardiometabolic risk via epigenetic changes. This study illustrates a new method to map epigenetic signatures of environmental factors for disease risk. Identification of epigenetic markers of dietary intake can provide insight into an individual's risk of cardiovascular disease and support the application of precision nutrition. Clinical Trial Registration: www.ClinicalTrials.gov, the Framingham Heart Study (FHS), NCT00005121; the Genetics of Lipid Lowering Drugs and Diet Network (GOLDN), NCT01023750.

A phenotype driven integrative framework uncovers molecular mechanisms of a rare hereditary thrombophilia

Antithrombin resistance is a rare subtype of hereditary thrombophilia caused by prothrombin gene variants, leading to thrombotic disorders. Recently, the Prothrombin Belgrade variant has been reported as a specific variant that leads to antithrombin resistance in two Serbian families with thrombosis. However, due to clinical data scarcity and the inapplicability of traditional genome-wide association studies (GWAS), a broader perspective on molecular and phenotypic mechanisms associated with the Prothrombin Belgrade variant is yet to be uncovered. Here, we propose an integrative framework to address the lack of genomic samples and support the genomic signal from the full genome sequences of five heterozygous subjects by integrating it with subjects' phenotypes and the genes' molecular interactions. Our goal is to identify candidate thrombophilia-related genes for which our subjects possess germline variants by focusing on the resulting gene clusters of our integrative framework. We applied a Non-negative Matrix Tri-Factorization-based method to simultaneously integrate different data sources, taking into account the observed phenotypes. In other words, our data-integration framework reveals gene clusters involved with this rare disease by fusing different datasets. Our results are in concordance with the current literature about antithrombin resistance. We also found candidate disease-related genes that need to be further investigated. CD320, RTEL1, UCP2, APOA5 and PROZ participate in healthy-specific or disease-specific subnetworks involving thrombophilia-annotated genes and are related to general thrombophilia mechanisms according to the literature. Moreover, the ADRA2A and TBXA2R subnetworks analysis suggested that their variants may have a protective effect due to their connection with decreased platelet activation. The results show that our method can give insights into antithrombin resistance even if a small amount of genetic data is available. Our framework is also customizable, meaning that it applies to any other rare disease.

Assisted clustering of gene expression data using regulatory data from partially overlapping sets of individuals

BACKGROUND

As omics measurements profiled on different molecular layers are interconnected, integrative approaches that incorporate the regulatory effect from multi-level omics data are needed. When the multi-level omics data are from the same individuals, gene expression (GE) clusters can be identified using information from regulators like genetic variants and DNA methylation. When the multi-level omics data are from different individuals, the choice of integration approaches is limited.

METHODS

We developed an approach to improve GE clustering from microarray data by integrating regulatory data from different but partially overlapping sets of individuals. We achieve this through (1) decomposing gene expression into the regulated component and the other component that is not regulated by measured factors, (2) optimizing the clustering goodness-of-fit objective function. We do not require the availability of different omics measurements on all individuals. A certain amount of individual overlap between GE data and the regulatory data is adequate for modeling the regulation, thus improving GE clustering.

RESULTS

A simulation study shows that the performance of the proposed approach depends on the strength of the GE-regulator relationship, degree of missingness, data dimensionality, sample size, and the number of clusters. Across the various simulation settings, the proposed method shows competitive performance in terms of accuracy compared to the alternative K-means clustering method, especially when the clustering structure is due mostly to the regulated component, rather than the unregulated component. We further validate the approach with an application to 8,902 Framingham Heart Study participants with data on up to 17,873 genes and regulation information of DNA methylation and genotype from different but partially overlapping sets of participants. We identify clustering structures of genes associated with pulmonary function while incorporating the predicted regulation effect from the measured regulators. We further investigate the over-representation of these GE clusters in pathways of other diseases that may be related to lung function and respiratory health.

CONCLUSION

We propose a novel approach for clustering GE with the assistance of regulatory data that allowed for different but partially overlapping sets of individuals to be included in different omics data.

In silico discovery of blood cell macromolecular associations

Background

Physical molecular interactions are the basis of intracellular signalling and gene regulatory networks, and comprehensive, accessible databases are needed for their discovery. Highly correlated transcripts may reflect important functional associations, but identification of such associations from primary data are cumbersome. We have constructed and adapted a user-friendly web application to discover and identify putative macromolecular associations in human peripheral blood based on significant correlations at the transcriptional level.

Methods

The blood transcriptome was characterized by quantification of 17,328 RNA species, including 341 mature microRNAs in 105 clinically well-characterized postmenopausal women. Intercorrelation of detected transcripts signal levels generated a matrix with > 150 million correlations recognizing the human blood RNA interactome. The correlations with calculated adjusted p-values were made easily accessible by a novel web application.

Results

We found that significant transcript correlations within the giant matrix reflect experimentally documented interactions involving select ubiquitous blood relevant transcription factors (CREB1, GATA1, and the glucocorticoid receptor (GR, NR3C1)). Their responsive genes recapitulated up to 91% of these as significant correlations, and were replicated in an independent cohort of 1204 individual blood samples from the Framingham Heart Study. Furthermore, experimentally documented mRNAs/miRNA associations were also reproduced in the matrix, and their predicted functional co-expression described. The blood transcript web application is available at http://app.uio.no/med/klinmed/correlation-browser/blood/index.php and works on all commonly used internet browsers.

Conclusions

Using in silico analyses and a novel web application, we found that correlated blood transcripts across 105 postmenopausal women reflected experimentally proven molecular associations. Furthermore, the associations were reproduced in a much larger and more heterogeneous cohort and should therefore be generally representative. The web application lends itself to be a useful hypothesis generating tool for identification of regulatory mechanisms in complex biological data sets.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-022-01077-3.

Identification of Cardiovascular Disease-Related Genes Based on the Co-Expression Network Analysis of Genome-Wide Blood Transcriptome

Inference of co-expression network and identification of disease-related modules and gene sets can help us understand disease-related molecular pathophysiology. We aimed to identify a cardiovascular disease (CVD)-related transcriptomic signature, specifically, in peripheral blood tissue, based on differential expression (DE) and differential co-expression (DcoE) analyses. Publicly available blood sample datasets for coronary artery disease (CAD) and acute coronary syndrome (ACS) statuses were integrated to establish a co-expression network. A weighted gene co-expression network analysis was used to construct modules that include genes with highly correlated expression values. The DE criterion is a linear regression with module eigengenes for module-specific genes calculated from principal component analysis and disease status as the dependent and independent variables, respectively. The DcoE criterion is a paired t-test for intramodular connectivity between disease and matched control statuses. A total of 21 and 23 modules were established from CAD status- and ACS-related datasets, respectively, of which six modules per disease status (i.e., obstructive CAD and ACS) were selected based on the DE and DcoE criteria. For each module, gene–gene interactions with extremely high correlation coefficients were individually selected under the two conditions. Genes displaying a significant change in the number of edges (gene–gene interaction) were selected. A total of 6, 10, and 7 genes in each of the three modules were identified as potential CAD status-related genes, and 14 and 8 genes in each of the two modules were selected as ACS-related genes. Our study identified gene sets and genes that were dysregulated in CVD blood samples. These findings may contribute to the understanding of CVD pathophysiology.

ReDisX, a machine learning approach, rationalizes rheumatoid arthritis and coronary artery disease patients uniquely upon identifying subpopulation differentiation markers from their genomic data

Diseases originate at the molecular-genetic layer, manifest through altered biochemical homeostasis, and develop symptoms later. Hence, symptomatic diagnosis is inadequate to explain the underlying molecular-genetic abnormality and individual genomic disparities. The current trends include molecular-genetic information relying on algorithms to recognize the disease subtypes through gene expressions. Despite their disposition toward disease-specific heterogeneity and cross-disease homogeneity, a gap still exists in describing the extent of homogeneity within the heterogeneous subpopulation of different diseases. They are limited to obtaining the holistic sense of the whole genome-based diagnosis resulting in inaccurate diagnosis and subsequent management. Addressing those ambiguities, our proposed framework, ReDisX, introduces a unique classification system for the patients based on their genomic signatures. In this study, it is a scalable machine learning algorithm deployed to re-categorize the patients with rheumatoid arthritis and coronary artery disease. It reveals heterogeneous subpopulations within a disease and homogenous subpopulations across different diseases. Besides, it identifies granzyme B (GZMB) as a subpopulation-differentiation marker that plausibly serves as a prominent indicator for GZMB-targeted drug repurposing. The ReDisX framework offers a novel strategy to redefine disease diagnosis through characterizing personalized genomic signatures. It may rejuvenate the landscape of precision and personalized diagnosis and a clue to drug repurposing.

Associations of carotid intima media thickness with gene expression in whole blood and genetically predicted gene expression across 48 tissues

Carotid intima media thickness (cIMT) is a biomarker of subclinical atherosclerosis and a predictor of future cardiovascular events. Identifying associations between gene expression levels and cIMT may provide insight to atherosclerosis etiology. Here, we use two approaches to identify associations between mRNA levels and cIMT: differential gene expression analysis in whole blood and S-PrediXcan. We used microarrays to measure genome-wide whole blood mRNA levels of 5647 European individuals from four studies. We examined the association of mRNA levels with cIMT adjusted for various potential confounders. Significant associations were tested for replication in three studies totaling 3943 participants. Next, we applied S-PrediXcan to summary statistics from a cIMT genome-wide association study (GWAS) of 71 128 individuals to estimate the association between genetically determined mRNA levels and cIMT and replicated these analyses using S-PrediXcan on an independent GWAS on cIMT that included 22 179 individuals from the UK Biobank. mRNA levels of TNFAIP3, CEBPD and METRNL were inversely associated with cIMT, but these associations were not significant in the replication analysis. S-PrediXcan identified associations between cIMT and genetically determined mRNA levels for 36 genes, of which six were significant in the replication analysis, including TLN2, which had not been previously reported for cIMT. There was weak correlation between our results using differential gene expression analysis and S-PrediXcan. Differential expression analysis and S-PrediXcan represent complementary approaches for the discovery of associations between phenotypes and gene expression. Using these approaches, we prioritize TNFAIP3, CEBPD, METRNL and TLN2 as new candidate genes whose differential expression might modulate cIMT.

Single-Cell Analysis Uncovers Osteoblast Factor Growth Differentiation Factor 10 as Mediator of Vascular Smooth Muscle Cell Phenotypic Modulation Associated with Plaque Rupture in Human Carotid Artery Disease

(1) Background: Vascular smooth muscle cells (VSMCs) undergo a complex phenotypic switch in response to atherosclerosis environmental triggers, contributing to atherosclerosis disease progression. However, the complex heterogeneity of VSMCs and how VSMC dedifferentiation affects human carotid artery disease (CAD) risk has not been clearly established. (2) Method: A single-cell RNA sequencing analysis of CD45⁻ cells derived from the atherosclerotic aorta of Apolipoprotein E-deficient (Apoe^−/−) mice on a normal cholesterol diet (NCD) or a high cholesterol diet (HCD), respecting the site-specific predisposition to atherosclerosis was performed. Growth Differentiation Factor 10 (GDF10) role in VSMCs phenotypic switch was investigated via flow cytometry, immunofluorescence in human atherosclerotic plaques. (3) Results: scRNAseq analysis revealed the transcriptomic profile of seven clusters, five of which showed disease-relevant gene signature of VSMC macrophagic calcific phenotype, VSMC mesenchymal chondrogenic phenotype, VSMC inflammatory and fibro-phenotype and VSMC inflammatory phenotype. Osteoblast factor GDF10 involved in ossification and osteoblast differentiation emerged as a hallmark of VSMCs undergoing phenotypic switch. Under hypercholesteremia, GDF10 triggered VSMC osteogenic switch in vitro. The abundance of GDF10 expressing osteogenic-like VSMCs cells was linked to the occurrence of carotid artery disease (CAD) events. (4) Conclusions: Taken together, these results provide evidence about GDF10-mediated VSMC osteogenic switch, with a likely detrimental role in atherosclerotic plaque stability.

Comprehensive Statistical and Bioinformatics Analysis in the Deciphering of Putative Mechanisms by Which Lipid-Associated GWAS Loci Contribute to Coronary Artery Disease

The study was designed to evaluate putative mechanisms by which lipid-associated loci identified by genome-wide association studies (GWAS) are involved in the molecular pathogenesis of coronary artery disease (CAD) using a comprehensive statistical and bioinformatics analysis. A total of 1700 unrelated individuals of Slavic origin from the Central Russia, including 991 CAD patients and 709 healthy controls were examined. Sixteen lipid-associated GWAS loci were selected from European studies and genotyped using the MassArray-4 system. The polymorphisms were associated with plasma lipids such as total cholesterol (rs12328675, rs4846914, rs55730499, and rs838880), LDL-cholesterol (rs3764261, rs55730499, rs1689800, and rs838880), HDL-cholesterol (rs3764261) as well as carotid intima-media thickness/CIMT (rs12328675, rs11220463, and rs1689800). Polymorphisms such as rs4420638 of APOC1 (p = 0.009), rs55730499 of LPA (p = 0.0007), rs3136441 of F2 (p < 0.0001), and rs6065906 of PLTP (p = 0.002) showed significant associations with the risk of CAD, regardless of sex, age, and body mass index. A majority of the observed associations were successfully replicated in large independent cohorts. Bioinformatics analysis allowed establishing (1) phenotype-specific and shared epistatic gene–gene and gene–smoking interactions contributing to all studied cardiovascular phenotypes; (2) lipid-associated GWAS loci might be allele-specific binding sites for transcription factors from gene regulatory networks controlling multifaceted molecular mechanisms of atherosclerosis.

Epigenome-wide association study of mitochondrial genome copy number

We conducted cohort- and race-specific epigenome-wide association analyses of mitochondrial deoxyribonucleic acid (mtDNA) copy number (mtDNA CN) measured in whole blood from participants of African and European origins in five cohorts (n = 6182, mean age = 57-67 years, 65% women). In the meta-analysis of all the participants, we discovered 21 mtDNA CN-associated DNA methylation sites (CpG) (P < 1 × 10-7), with a 0.7-3.0 standard deviation increase (3 CpGs) or decrease (18 CpGs) in mtDNA CN corresponding to a 1% increase in DNA methylation. Several significant CpGs have been reported to be associated with at least two risk factors (e.g. chronological age or smoking) for cardiovascular disease (CVD). Five genes [PR/SET domain 16, nuclear receptor subfamily 1 group H member 3 (NR1H3), DNA repair protein, DNA polymerase kappa and decaprenyl-diphosphate synthase subunit 2], which harbor nine significant CpGs, are known to be involved in mitochondrial biosynthesis and functions. For example, NR1H3 encodes a transcription factor that is differentially expressed during an adipose tissue transition. The methylation level of cg09548275 in NR1H3 was negatively associated with mtDNA CN (effect size = -1.71, P = 4 × 10-8) and was positively associated with the NR1H3 expression level (effect size = 0.43, P = 0.0003), which indicates that the methylation level in NR1H3 may underlie the relationship between mtDNA CN, the NR1H3 transcription factor and energy expenditure. In summary, the study results suggest that mtDNA CN variation in whole blood is associated with DNA methylation levels in genes that are involved in a wide range of mitochondrial activities. These findings will help reveal molecular mechanisms between mtDNA CN and CVD.

Large-scale serum analysis identifies unique systemic biomarkers in psoriasis and hidradenitis suppurativa

BACKGROUND

Hidradenitis Suppurativa (HS) is now recognized as a systemic inflammatory disease, sharing molecular similarities with psoriasis. Direct comparison of the systemic inflammation in HS with psoriasis is lacking.

OBJECTIVES

To evaluate the serum proteome of HS and psoriasis, and to identify biomarkers associated with disease severity.

METHODS

In this cross-sectional study, 1,536 serum proteins were assessed using the Olink Explore (Proximity Extension Assay/PEA) high-throughput panel in moderate-to-severe HS (n=11), psoriasis (n=10) and age- and BMI-matched healthy controls (n=10).

RESULTS

HS displayed an overall greater dysregulation of circulating proteins, with 434 differentially expressed proteins (|FCH|≥1.2, p-value≤0.05) in HS versus controls, 138 in psoriasis versus controls, and 503 between HS and psoriasis. IL-17A levels and Th1/Th17 pathway enrichment were comparable between diseases, while HS presented greater TNF and IL-1β-related signaling. Th17-associated markers, PI3 and LCN2, were able to accurately differentiate psoriasis from HS. Both diseases presented increases of atherosclerosis-related proteins. Robust correlations between clinical severity scores and immune and atherosclerosis-related proteins were observed across both diseases.

CONCLUSIONS

HS and psoriasis share significant Th1/Th17 enrichment and upregulation of atherosclerosis-related proteins. Nevertheless, despite the greater body surface area involved in psoriasis, HS presents a greater serum inflammatory burden.

DNA methylation and gene expression integration in cardiovascular disease

BACKGROUND

The integration of different layers of omics information is an opportunity to tackle the complexity of cardiovascular diseases (CVD) and to identify new predictive biomarkers and potential therapeutic targets. Our aim was to integrate DNA methylation and gene expression data in an effort to identify biomarkers related to cardiovascular disease risk in a community-based population. We accessed data from the Framingham Offspring Study, a cohort study with data on DNA methylation (Infinium HumanMethylation450 BeadChip; Illumina) and gene expression (Human Exon 1.0 ST Array; Affymetrix). Using the MOFA2 R package, we integrated these data to identify biomarkers related to the risk of presenting a cardiovascular event.

RESULTS

Four independent latent factors (9, 19, 21-only in women-and 27), driven by DNA methylation, were associated with cardiovascular disease independently of classical risk factors and cell-type counts. In a sensitivity analysis, we also identified factor 21 as associated with CVD in women. Factors 9, 21 and 27 were also associated with coronary heart disease risk. Moreover, in a replication effort in an independent study three of the genes included in factor 27 were also present in a factor identified to be associated with myocardial infarction (CDC42BPB, MAN2A2 and RPTOR). Factor 9 was related to age and cell-type proportions; factor 19 was related to age and B cells count; factor 21 pointed to human immunodeficiency virus infection-related pathways and inflammation; and factor 27 was related to lifestyle factors such as alcohol consumption, smoking and body mass index. Inclusion of factor 21 (only in women) improved the discriminative and reclassification capacity of the Framingham classical risk function and factor 27 improved its discrimination.

CONCLUSIONS

Unsupervised multi-omics data integration methods have the potential to provide insights into the pathogenesis of cardiovascular diseases. We identified four independent factors (one only in women) pointing to inflammation, endothelium homeostasis, visceral fat, cardiac remodeling and lifestyles as key players in the determination of cardiovascular risk. Moreover, two of these factors improved the predictive capacity of a classical risk function.

Subtypes identification on heart failure with preserved ejection fraction via network enhancement fusion using multi-omics data

Heart failure with preserved ejection fraction (HFpEF) is associated with multiple etiologic and pathophysiologic factors. HFpEF leads to significant cardiovascular morbidity and mortality. There are various reasons that fail to identify effective therapeutic interventions for HFpEF, primarily due to its clinical heterogeneity causing significant difficulties in determining physiologic and prognostic implications for this syndrome. Thus, identifying clinical subtypes using multi-omics data has great implications for efficient treatment and prognosis of HFpEF patients. Here we proposed to integrate mRNA, DNA methylation and microRNA (miRNA) expression data of HFpEF with a similarity network fusion (SNF) method following a network enhancement (ne-SNF) denoising technique to form a fused network. A spectral clustering method was then used to obtain clusters of patient subtypes. Experiments on HFpEF datasets demonstrated that ne-SNF significantly outperforms single data subtype analysis and other integrated methods. The identified subgroups were shown to have statistically significant differences in survival. Two HFpEF subtypes were defined: a high-risk group (16.8%) and a low-risk group (83.2%). The 5-year mortality rates were 63.3% and 33.0% for the high- and low-risk group, respectively. After adjusting for the effects of clinical covariates, HFpEF patients in the high-risk group were 2.43 times more likely to die than the low-risk group. A total of 157 differentially expressed (DE) mRNAs, 2199 abnormal methylations and 121 DE miRNAs were identified between two subtypes. They were also enriched in many HFpEF-related biological processes or pathways. The ne-SNF method provides a novel pipeline for subtype identification in integrated analysis of multi-omics data.

Global identification of S-palmitoylated proteins and detection of palmitoylating (DHHC) enzymes in heart

High-throughput experiments suggest that almost 20% of human proteins may be S-palmitoylatable, a post-translational modification (PTM) whereby fatty acyl chains, most commonly palmitoyl chain, are linked to cysteine thiol groups that impact on protein trafficking, distribution and function. In human, protein S-palmitoylation is mediated by a group of 23 palmitoylating 'Asp-His-His-Cys' domain-containing (DHHC) enzymes. There is no information on the scope of protein S-palmitoylation, or the pattern of DHHC enzyme expression, in the heart. We used resin-assisted capture to pull down S-palmitoylated proteins from human, dog, and rat hearts, followed by proteomic search to identify proteins in the pulldowns. We identified 454 proteins present in at least 2 species-specific pulldowns. These proteins are operationally called 'cardiac palmitoylome'. Enrichment analysis based on Gene Ontology terms 'cellular component' indicated that cardiac palmitoylome is involved in cell-cell and cell-substrate junctions, plasma membrane microdomain organization, vesicular trafficking, and mitochondrial enzyme organization. Importantly, cardiac palmitoylome is uniquely enriched in proteins participating in the organization and function of t-tubules, costameres and intercalated discs, three microdomains critical for excitation-contraction coupling and intercellular communication of cardiomyocytes. We validated antibodies targeting DHHC enzymes, and detected eleven of them expressed in hearts across species. In conclusion, we provide resources useful for investigators interested in studying protein S-palmitoylation and its regulation by DHHC enzymes in the heart. We also discuss challenges in these efforts, and suggest methods and tools that should be developed to overcome these challenges.

Genome-wide transcriptome study using deep RNA sequencing for myocardial infarction and coronary artery calcification

Background

Coronary artery calcification (CAC) is a noninvasive measure of coronary atherosclerosis, the proximal pathophysiology underlying most cases of myocardial infarction (MI). We sought to identify expression signatures of early MI and subclinical atherosclerosis in the Framingham Heart Study (FHS). In this study, we conducted paired-end RNA sequencing on whole blood collected from 198 FHS participants (55 with a history of early MI, 72 with high CAC without prior MI, and 71 controls free of elevated CAC levels or history of MI). We applied DESeq2 to identify coding-genes and long intergenic noncoding RNAs (lincRNAs) differentially expressed in MI and high CAC, respectively, compared with the control.

Results

On average, 150 million paired-end reads were obtained for each sample. At the false discovery rate (FDR) < 0.1, we found 68 coding genes and 2 lincRNAs that were differentially expressed in early MI versus controls. Among them, 60 coding genes were detectable and thus tested in an independent RNA-Seq data of 807 individuals from the Rotterdam Study, and 8 genes were supported by p value and direction of the effect. Immune response, lipid metabolic process, and interferon regulatory factor were enriched in these 68 genes. By contrast, only 3 coding genes and 1 lincRNA were differentially expressed in high CAC versus controls. APOD, encoding a component of high-density lipoprotein, was significantly downregulated in both early MI (FDR = 0.007) and high CAC (FDR = 0.01) compared with controls.

Conclusions

We identified transcriptomic signatures of early MI that include differentially expressed protein-coding genes and lincRNAs, suggesting important roles for protein-coding genes and lincRNAs in the pathogenesis of MI.

Mir-30d Regulates Cardiac Remodeling by Intracellular and Paracrine Signaling

RATIONALE

Previous translational studies implicate plasma extracellular microRNA-30d (miR-30d) as a biomarker in left ventricular remodeling and clinical outcome in heart failure (HF) patients, although precise mechanisms remain obscure.

OBJECTIVE

To investigate the mechanism of miR-30d-mediated cardioprotection in HF.

METHODS AND RESULTS

In rat and mouse models of ischemic HF, we show that miR-30d gain of function (genetic, lentivirus, or agomiR-mediated) improves cardiac function, decreases myocardial fibrosis, and attenuates cardiomyocyte (CM) apoptosis. Genetic or locked nucleic acid-based knock-down of miR-30d expression potentiates pathological left ventricular remodeling, with increased dysfunction, fibrosis, and cardiomyocyte death. RNA sequencing of in vitro miR-30d gain and loss of function, together with bioinformatic prediction and experimental validation in cardiac myocytes and fibroblasts, were used to identify and validate direct targets of miR-30d. miR-30d expression is selectively enriched in cardiomyocytes, induced by hypoxic stress and is acutely protective, targeting MAP4K4 (mitogen-associate protein kinase 4) to ameliorate apoptosis. Moreover, miR-30d is secreted primarily in extracellular vesicles by cardiomyocytes and inhibits fibroblast proliferation and activation by directly targeting integrin α5 in the acute phase via paracrine signaling to cardiac fibroblasts. In the chronic phase of ischemic remodeling, lower expression of miR-30d in the heart and plasma extracellular vesicles is associated with adverse remodeling in rodent models and human subjects and is linked to whole-blood expression of genes implicated in fibrosis and inflammation, consistent with observations in model systems.

CONCLUSIONS

These findings provide the mechanistic underpinning for the cardioprotective association of miR-30d in human HF. More broadly, our findings support an emerging paradigm involving intercellular communication of extracellular vesicle-contained miRNAs (microRNAs) to transregulate distinct signaling pathways across cell types. Functionally validated RNA biomarkers and their signaling networks may warrant further investigation as novel therapeutic targets in HF.

miRNA Mechanisms Underlying the Association of Beta Blocker Use and Bone Mineral Density

Osteoporosis is a debilitating and costly disease that causes fractures in 33% of women and 20% of men over the age of 50 years. Recent studies have shown that beta blocker (BB) users have higher bone mineral density (BMD) and decreased risk of fracture compared with non-users. The mechanism underlying this association is thought to be due to suppression of adrenergic signaling in osteoblasts, which leads to increased BMD in rodent models; however, the mechanism in humans is unknown. Also, several miRNAs are associated with adrenergic signaling and BMD in separate studies. To investigate potential miRNA mechanisms, we performed a cross-sectional analysis using clinical data, dual-energy X-ray absorptiometry (DXA) scans, and miRNA and mRNA profiling of whole blood from the Framingham Study's Offspring Cohort. We found nine miRNAs associated with BB use and increased BMD. In parallel network analyses, we discovered a subnetwork associated with BMD and BB use containing two of these nine miRNAs, miR-19a-3p and miR-186-5p. To strengthen this finding, we showed that these two miRNAs had significantly higher expression in individuals without incident fracture compared with those with fracture in an external data set. We also noted a similar trend in association between these miRNA and Z-score as calculated from heel ultrasound measures in two external cohorts (SOS-Hip and SHIP-TREND). Because miR-19a directly targets the ADRB1 mRNA transcript, we propose BB use may downregulate ADRB1 expression in osteoblasts through increased miR-19a-3p expression. We used enrichment analysis of miRNA targets to find potential indirect effects through insulin and parathyroid hormone signaling. This analysis provides a starting point for delineating the role of miRNA on the association between BB use and BMD. © 2020 American Society for Bone and Mineral Research (ASBMR).

Carbohydrate and fat intake associated with risk of metabolic diseases through epigenetics of CPT1A

BACKGROUND

Epigenome-wide association studies identified the cg00574958 DNA methylation site at the carnitine palmitoyltransferase-1A (CPT1A) gene to be associated with reduced risk of metabolic diseases (hypertriglyceridemia, obesity, type 2 diabetes, hypertension, metabolic syndrome), but the mechanism underlying these associations is unknown.

OBJECTIVES

We aimed to elucidate whether carbohydrate and fat intakes modulate cg00574958 methylation and the risk of metabolic diseases.

METHODS

We examined associations between carbohydrate (CHO) and fat (FAT) intake, as percentages of total diet energy, and the CHO/FAT ratio with CPT1A-cg00574958, and the risk of metabolic diseases in 3 populations (Genetics of Lipid Lowering Drugs and Diet Network, n = 978; Framingham Heart Study, n = 2331; and REgistre GIroní del COR study, n = 645) while adjusting for confounding factors. To understand possible causal effects of dietary intake on the risk of metabolic diseases, we performed meta-analysis, CPT1A transcription analysis, and mediation analysis with CHO and FAT intakes as exposures and cg00574958 methylation as the mediator.

RESULTS

We confirmed strong associations of cg00574958 methylation with metabolic phenotypes (BMI, triglyceride, glucose) and diseases in all 3 populations. Our results showed that CHO intake and CHO/FAT ratio were positively associated with cg00574958 methylation, whereas FAT intake was negatively correlated with cg00574958 methylation. Meta-analysis further confirmed this strong correlation, with β = 58.4 ± 7.27, P = 8.98 x 10-16 for CHO intake; β = -36.4 ± 5.95, P = 9.96 x 10-10 for FAT intake; and β = 3.30 ± 0.49, P = 1.48 x 10-11 for the CHO/FAT ratio. Furthermore, CPT1A mRNA expression was negatively associated with CHO intake, and positively associated with FAT intake, and metabolic phenotypes. Mediation analysis supports the hypothesis that CHO intake induces CPT1A methylation, hence reducing the risk of metabolic diseases, whereas FAT intake inhibits CPT1A methylation, thereby increasing the risk of metabolic diseases.

CONCLUSIONS

Our results suggest that the proportion of total energy supplied by CHO and FAT can have a causal effect on the risk of metabolic diseases via the epigenetic status of CPT1A.Study registration at https://www.clinicaltrials.gov/: the Genetics of Lipid Lowering Drugs and Diet Network (GOLDN)-NCT01023750; and the Framingham Heart Study (FHS)-NCT00005121.

Smoking-related changes in DNA methylation and gene expression are associated with cardio-metabolic traits

BACKGROUND

Tobacco smoking is a well-known modifiable risk factor for many chronic diseases, including cardiovascular disease (CVD). One of the proposed underlying mechanism linking smoking to disease is via epigenetic modifications, which could affect the expression of disease-associated genes. Here, we conducted a three-way association study to identify the relationship between smoking-related changes in DNA methylation and gene expression and their associations with cardio-metabolic traits.

RESULTS

We selected 2549 CpG sites and 443 gene expression probes associated with current versus never smokers, from the largest epigenome-wide association study and transcriptome-wide association study to date. We examined three-way associations, including CpG versus gene expression, cardio-metabolic trait versus CpG, and cardio-metabolic trait versus gene expression, in the Rotterdam study. Subsequently, we replicated our findings in The Cooperative Health Research in the Region of Augsburg (KORA) study. After correction for multiple testing, we identified both cis- and trans-expression quantitative trait methylation (eQTM) associations in blood. Specifically, we found 1224 smoking-related CpGs associated with at least one of the 443 gene expression probes, and 200 smoking-related gene expression probes to be associated with at least one of the 2549 CpGs. Out of these, 109 CpGs and 27 genes were associated with at least one cardio-metabolic trait in the Rotterdam Study. We were able to replicate the associations with cardio-metabolic traits of 26 CpGs and 19 genes in the KORA study. Furthermore, we identified a three-way association of triglycerides with two CpGs and two genes (GZMA; CLDND1), and BMI with six CpGs and two genes (PID1; LRRN3). Finally, our results revealed the mediation effect of cg03636183 (F2RL3), cg06096336 (PSMD1), cg13708645 (KDM2B), and cg17287155 (AHRR) within the association between smoking and LRRN3 expression.

CONCLUSIONS

Our study indicates that smoking-related changes in DNA methylation and gene expression are associated with cardio-metabolic risk factors. These findings may provide additional insights into the molecular mechanisms linking smoking to the development of CVD.

Identifying the RNA signatures of coronary artery disease from combined lncRNA and mRNA expression profiles

Coronary artery disease (CAD) is the most common cardiovascular disease. CAD research has greatly progressed during the past decade. mRNA is a traditional and popular pipeline to investigate various disease, including CAD. Compared with mRNA, lncRNA has better stability and thus may serve as a better disease indicator in blood. Investigating potential CAD-related lncRNAs and mRNAs will greatly contribute to the diagnosis and treatment of CAD. In this study, a computational analysis was conducted on patients with CAD by using a comprehensive transcription dataset with combined mRNA and lncRNA expression data. Several machine learning algorithms, including feature selection methods and classification algorithms, were applied to screen for the most CAD-related RNA molecules. Decision rules were also reported to provide a quantitative description about the effect of these RNA molecules on CAD progression. These new findings (CAD-related RNA molecules and rules) can help understand mRNA and lncRNA expression levels in CAD.

An integrative multi-omics analysis to identify candidate DNA methylation biomarkers related to prostate cancer risk

It remains elusive whether some of the associations identified in genome-wide association studies of prostate cancer (PrCa) may be due to regulatory effects of genetic variants on CpG sites, which may further influence expression of PrCa target genes. To search for CpG sites associated with PrCa risk, here we establish genetic models to predict methylation (N = 1,595) and conduct association analyses with PrCa risk (79,194 cases and 61,112 controls). We identify 759 CpG sites showing an association, including 15 located at novel loci. Among those 759 CpG sites, methylation of 42 is associated with expression of 28 adjacent genes. Among 22 genes, 18 show an association with PrCa risk. Overall, 25 CpG sites show consistent association directions for the methylation-gene expression-PrCa pathway. We identify DNA methylation biomarkers associated with PrCa, and our findings suggest that specific CpG sites may influence PrCa via regulating expression of candidate PrCa target genes.

Low oxygen saturation during sleep reduces CD1D and RAB20 expressions that are reversed by CPAP therapy

Background

Sleep Disordered Breathing (SDB) is associated with a wide range of pathophysiological changes due, in part, to hypoxemia during sleep. We sought to identify gene transcription associations with measures of SDB and hypoxemia during sleep, and study their response to treatment.

Methods

In two discovery cohorts, Framingham Offspring Study (FOS; N = 571) and the Multi-Ethnic Study of Atherosclerosis (MESA; N = 580), we studied gene expression in peripheral blood mononuclear cells in association with three measures of SDB: Apnea Hypopnea Index (AHI); average oxyhemoglobin saturation (avgO2) during sleep; and minimum oxyhemoglobin saturation (minO2) during sleep. Associated genes were used for analysis of gene expression in the blood of 15 participants with moderate or severe obstructive sleep apnea (OSA) from the Heart Biomarkers In Apnea Treatment (HeartBEAT) trial. These genes were studied pre- and post-treatment (three months) with continuous positive airway pressure (CPAP). We also performed Gene Set Enrichment Analysis (GSEA) on all traits and cohort analyses.

Findings

Twenty-two genes were associated with SDB traits in both MESA and FOS. Of these, lower expression of CD1D and RAB20 was associated with lower avgO2 in MESA and FOS. CPAP treatment increased the expression of these genes in HeartBEAT participants. Immunity and inflammation pathways were up-regulated in subjects with lower avgO2; i.e., in those with a more severe SDB phenotype (MESA), whereas immuno-inflammatory processes were down-regulated following CPAP treatment (HeartBEAT).

Interpretation

Low oxygen saturation during sleep is associated with alterations in gene expression and transcriptional programs that are partially reversed by CPAP treatment.

Deary I, Liu Q, Liu S, Elosua R, DeMeo DL, Ordovás JM. Epigenomic Assessment of Cardiovascular Disease Risk and Interactions With Traditional Risk Metrics

Background Epigenome-wide association studies for cardiometabolic risk factors have discovered multiple loci associated with incident cardiovascular disease (CVD). However, few studies have sought to directly optimize a predictor of CVD risk. Furthermore, it is challenging to train multivariate models across multiple studies in the presence of study- or batch effects. Methods and Results Here, we analyzed existing DNA methylation data collected using the Illumina HumanMethylation450 microarray to create a predictor of CVD risk across 3 cohorts: Women's Health Initiative, Framingham Heart Study Offspring Cohort, and Lothian Birth Cohorts. We trained Cox proportional hazards-based elastic net regressions for incident CVD separately in each cohort and used a recently introduced cross-study learning approach to integrate these individual scores into an ensemble predictor. The methylation-based risk score was associated with CVD time-to-event in a held-out fraction of the Framingham data set (hazard ratio per SD=1.28, 95% CI, 1.10-1.50) and predicted myocardial infarction status in the independent REGICOR (Girona Heart Registry) data set (odds ratio per SD=2.14, 95% CI, 1.58-2.89). These associations remained after adjustment for traditional cardiovascular risk factors and were similar to those from elastic net models trained on a directly merged data set. Additionally, we investigated interactions between the methylation-based risk score and both genetic and biochemical CVD risk, showing preliminary evidence of an enhanced performance in those with less traditional risk factor elevation. Conclusions This investigation provides proof-of-concept for a genome-wide, CVD-specific epigenomic risk score and suggests that DNA methylation data may enable the discovery of high-risk individuals who would be missed by alternative risk metrics.

Transcriptome analysis reveals the difference between "healthy" and "common" aging and their connection with age-related diseases

A key goal of aging research was to understand mechanisms underlying healthy aging and develop methods to promote the human healthspan. One approach is to identify gene regulations unique to healthy aging compared with aging in the general population (i.e., "common" aging). Here, we leveraged Genotype-Tissue Expression (GTEx) project data to investigate "healthy" and "common" aging gene expression regulations at a tissue level in humans and their interconnection with diseases. Using GTEx donors' disease annotations, we defined a "healthy" aging cohort for each tissue. We then compared the age-associated genes derived from this cohort with age-associated genes from the "common" aging cohort which included all GTEx donors; we also compared the "healthy" and "common" aging gene expressions with various disease-associated gene expressions to elucidate the relationships among "healthy," "common" aging and disease. Our analyses showed that 1. GTEx "healthy" and "common" aging shared a large number of gene regulations; 2. Despite the substantial commonality, "healthy" and "common" aging genes also showed distinct function enrichment, and "common" aging genes had a higher enrichment for disease genes; 3. Disease-associated gene regulations were overall different from aging gene regulations. However, for genes regulated by both, their regulation directions were largely consistent, implying some aging processes could increase the susceptibility to disease development; and 4. Possible protective mechanisms were associated with some "healthy" aging gene regulations. In summary, our work highlights several unique features of GTEx "healthy" aging program. This new knowledge could potentially be used to develop interventions to promote the human healthspan.

Endogenous oxidized phospholipids reprogram cellular metabolism and boost hyperinflammation

Pathogen-associated molecular patterns (PAMPs) have the capacity to couple inflammatory gene expression to changes in macrophage metabolism, both of which influence subsequent inflammatory activities. Similar to their microbial counterparts, several self-encoded damage-associated molecular patterns (DAMPs) induce inflammatory gene expression. However, whether this symmetry in host responses between PAMPs and DAMPs extends to metabolic shifts is unclear. Here we report that the self-encoded oxidized phospholipid oxPAPC alters the metabolism of macrophages exposed to lipopolysaccharide (LPS). While cells activated by LPS rely exclusively on glycolysis, macrophages exposed to oxPAPC also use mitochondrial respiration, feed the Krebs cycle with glutamine and favor the accumulation of oxaloacetate in the cytoplasm: this metabolite potentiates IL-1β production, resulting in hyperinflammation. Similar metabolic adaptions occur in vivo in hypercholesterolemic mice and human subjects. Drugs that interfere with oxPAPC-driven metabolic changes reduce atherosclerotic plaque formation in mice, thereby underscoring the importance of DAMP-mediated activities in pathophysiological conditions.

Taking Systems Medicine to Heart

Systems medicine is a holistic approach to deciphering the complexity of human physiology in health and disease. In essence, a living body is constituted of networks of dynamically interacting units (molecules, cells, organs, etc) that underlie its collective functions. Declining resilience because of aging and other chronic environmental exposures drives the system to transition from a health state to a disease state; these transitions, triggered by acute perturbations or chronic disturbance, manifest as qualitative shifts in the interactions and dynamics of the disease-perturbed networks. Understanding health-to-disease transitions poses a high-dimensional nonlinear reconstruction problem that requires deep understanding of biology and innovation in study design, technology, and data analysis. With a focus on the principles of systems medicine, this Review discusses approaches for deciphering this biological complexity from a novel perspective, namely, understanding how disease-perturbed networks function; their study provides insights into fundamental disease mechanisms. The immediate goals for systems medicine are to identify early transitions to cardiovascular (and other chronic) diseases and to accelerate the translation of new preventive, diagnostic, or therapeutic targets into clinical practice, a critical step in the development of personalized, predictive, preventive, and participatory (P4) medicine.

Epigenome-wide association study of DNA methylation and microRNA expression highlights novel pathways for human complex traits

DNA methylation (DNAm) and microRNAs (miRNAs) have been implicated in a wide-range of human diseases. While often studied in isolation, DNAm and miRNAs are not independent. We analyzed associations of expression of 283 miRNAs with DNAm at >400K CpG sites in whole blood obtained from 3565 individuals and identified 227 CpGs at which differential methylation was associated with the expression of 40 nearby miRNAs (cis-miR-eQTMs) at FDR<0.01, including 91 independent CpG sites at r² < 0.2. cis-miR-eQTMs were enriched for CpGs in promoter and polycomb-repressed state regions, and 60% were inversely associated with miRNA expression. Bidirectional Mendelian randomization (MR) analysis further identified 58 cis-miR-eQTMCpG-miRNA pairs where DNAm changes appeared to drive miRNA expression changes and opposite directional effects were unlikely. Integration of genetic variants in joint analyses revealed an average partial between cis-miR-eQTM CpGs and miRNAs of 2% after conditioning on site-specific genetic variation, suggesting that DNAm is an important epigenetic regulator of miRNA expression. Finally, two-step MR analysis was performed to identify putatively causal CpGs driving miRNA expression in relation to human complex traits. We found that an imprinted region on 14q32 that was previously identified in relation to age at menarche is enriched with cis-miR-eQTMs. Nine CpGs and three miRNAs at this locus tested causal for age at menarche, reflecting novel epigenetic-driven molecular pathways underlying this complex trait. Our study sheds light on the joint genetic and epigenetic regulation of miRNA expression and provides insights into the relations of miRNAs to their targets and to complex phenotypes.

The Cortical Neuroimmune Regulator TANK Affects Emotional Processing and Enhances Alcohol Drinking: A Translational Study

Alcohol abuse is a major public health problem worldwide. Understanding the molecular mechanisms that control regular drinking may help to reduce hazards of alcohol consumption. While immunological mechanisms have been related to alcohol drinking, most studies reported changes in immune function that are secondary to alcohol use. In this report, we analyse how the gene "TRAF family member-associated NF-κB activator" (TANK) affects alcohol drinking behavior. Based on our recent discovery in a large GWAS dataset that suggested an association of TANK, SNP rs197273, with alcohol drinking, we report that SNP rs197273 in TANK is associated both with gene expression (P = 1.16 × 10-19) and regional methylation (P = 5.90 × 10-25). A tank knock out mouse model suggests a role of TANK in alcohol drinking, anxiety-related behavior, as well as alcohol exposure induced activation of insular cortex NF-κB. Functional and structural neuroimaging studies among up to 1896 adolescents reveal that TANK is involved in the control of brain activity in areas of aversive interoceptive processing, including the insular cortex, but not in areas related to reinforcement, reward processing or impulsiveness. Our findings suggest that the cortical neuroimmune regulator TANK is associated with enhanced aversive emotional processing that better protects from the establishment of alcohol drinking behavior.

Integrating genetic, transcriptional, and biological information provides insights into obesity

Objective:

Indices of body fat distribution are heritable, but few genetic signals have been reported from genome-wide association studies (GWAS) of computed tomography (CT) imaging measurements of body fat distribution. We aimed to identify genes associated with adiposity traits and the key drivers that are central to adipose regulatory networks.

Subjects:

We analyzed gene transcript expression data in blood from participants in the Framingham Heart Study, a large community-based cohort (n up to 4,303), as well as implemented an integrative analysis of these data and existing biological information.

Results:

Our association analyses identified unique and common gene expression signatures across several adiposity traits, including body mass index, waist-hip ratio, waist circumference, and CT-measured indices, including volume and quality of visceral and subcutaneous adipose tissues. We identified six enriched KEGG pathways and two co-expression modules for further exploration of adipose regulatory networks. The integrative analysis revealed four gene sets (Apoptosis, p53 signaling pathway, Proteasome, Ubiquitin mediated proteolysis) and two co-expression modules with significant genetic variants and 94 key drivers/genes whose local networks were enriched with adiposity-associated genes, suggesting that these enriched pathways or modules have genetic effects on adiposity. Most identified key driver genes are involved in essential biological processes such as controlling cell cycle, DNA repair and degradation of regulatory proteins and are cancer related.

Conclusion:

Our integrative analysis of genetic, transcriptional and biological information provides a list of compelling candidates for further follow-up functional studies to uncover the biological mechanisms underlying obesity. These candidates highlight the value of examining CT-derived and central adiposity traits.

Genome-wide mapping of plasma protein QTLs identifies putatively causal genes and pathways for cardiovascular disease

Identifying genetic variants associated with circulating protein concentrations (protein quantitative trait loci; pQTLs) and integrating them with variants from genome-wide association studies (GWAS) may illuminate the proteome's causal role in disease and bridge a knowledge gap regarding SNP-disease associations. We provide the results of GWAS of 71 high-value cardiovascular disease proteins in 6861 Framingham Heart Study participants and independent external replication. We report the mapping of over 16,000 pQTL variants and their functional relevance. We provide an integrated plasma protein-QTL database. Thirteen proteins harbor pQTL variants that match coronary disease-risk variants from GWAS or test causal for coronary disease by Mendelian randomization. Eight of these proteins predict new-onset cardiovascular disease events in Framingham participants. We demonstrate that identifying pQTLs, integrating them with GWAS results, employing Mendelian randomization, and prospectively testing protein-trait associations holds potential for elucidating causal genes, proteins, and pathways for cardiovascular disease and may identify targets for its prevention and treatment.

Epigenomics and metabolomics reveal the mechanism of the APOA2-saturated fat intake interaction affecting obesity

Background

The putative functional variant -265T>C (rs5082) within the APOA2 promoter has shown consistent interactions with saturated fatty acid (SFA) intake to influence the risk of obesity.

Objective

The aim of this study was to implement an integrative approach to characterize the molecular basis of this interaction.

Design

We conducted an epigenome-wide scan on 80 participants carrying either the rs5082 CC or TT genotypes and consuming either a low-SFA (<22 g/d) or high-SFA diet (≥22 g/d), matched for age, sex, BMI, and diabetes status in the Boston Puerto Rican Health Study (BPRHS). We then validated the findings in selected participants in the Genetics of Lipid Lowering Drugs and Diet Network (GOLDN) Study (n = 379) and the Framingham Heart Study (FHS) (n = 243). Transcription and metabolomics analyses were conducted to determine the relation between epigenetic status, APOA2 mRNA expression, and blood metabolites.

Results

In the BPRHS, we identified methylation site cg04436964 as exhibiting significant differences between CC and TT participants consuming a high-SFA diet, but not among those consuming low-SFA. Similar results were observed in the GOLDN Study and the FHS. Additionally, in the FHS, cg04436964 methylation was negatively correlated with APOA2 expression in the blood of participants consuming a high-SFA diet. Furthermore, when consuming a high-SFA diet, CC carriers had lower APOA2 expression than those with the TT genotype. Lastly, metabolomic analysis identified 4 pathways as overrepresented by metabolite differences between CC and TT genotypes with high-SFA intake, including tryptophan and branched-chain amino acid (BCAA) pathways. Interestingly, these pathways were linked to rs5082-specific cg04436964 methylation differences in high-SFA consumers.

Conclusions

The epigenetic status of the APOA2 regulatory region is associated with SFA intake and APOA2 -265T>C genotype, promoting an APOA2 expression difference between APOA2 genotypes on a high-SFA diet, and modulating BCAA and tryptophan metabolic pathways. These findings identify potential mechanisms by which this highly reproducible gene-diet interaction influences obesity risk, and contribute new insights to ongoing investigations of the relation between SFA and human health. This study was registered at clinicaltrials.gov as NCT03452787.

EPS-LASSO: test for high-dimensional regression under extreme phenotype sampling of continuous traits

Motivation

Extreme phenotype sampling (EPS) is a broadly-used design to identify candidate genetic factors contributing to the variation of quantitative traits. By enriching the signals in extreme phenotypic samples, EPS can boost the association power compared to random sampling. Most existing statistical methods for EPS examine the genetic factors individually, despite many quantitative traits have multiple genetic factors underlying their variation. It is desirable to model the joint effects of genetic factors, which may increase the power and identify novel quantitative trait loci under EPS. The joint analysis of genetic data in high-dimensional situations requires specialized techniques, e.g. the least absolute shrinkage and selection operator (LASSO). Although there are extensive research and application related to LASSO, the statistical inference and testing for the sparse model under EPS remain unknown.

Results

We propose a novel sparse model (EPS-LASSO) with hypothesis test for high-dimensional regression under EPS based on a decorrelated score function. The comprehensive simulation shows EPS-LASSO outperforms existing methods with stable type I error and FDR control. EPS-LASSO can provide a consistent power for both low- and high-dimensional situations compared with the other methods dealing with high-dimensional situations. The power of EPS-LASSO is close to other low-dimensional methods when the causal effect sizes are small and is superior when the effects are large. Applying EPS-LASSO to a transcriptome-wide gene expression study for obesity reveals 10 significant body mass index associated genes. Our results indicate that EPS-LASSO is an effective method for EPS data analysis, which can account for correlated predictors.

Availability and implementation

The source code is available at https://github.com/xu1912/EPSLASSO.

Supplementary information

Supplementary data are available at Bioinformatics online.

Epigenome-Wide Association Study of Soluble Tumor Necrosis Factor Receptor 2 Levels in the Framingham Heart Study

Background: Transmembrane tumor necrosis factor (TNF) receptors are involved in inflammatory, apoptotic, and proliferative processes. In the bloodstream, soluble TNF receptor II (sTNFR2) can modify the inflammatory response of immune cells and is predictive of cardiovascular disease risk. We hypothesize that sTNFR2 is associated with epigenetic modifications of circulating leukocytes, which may relate to the pathophysiology underlying atherogenic risk. Methods: We conducted an epigenome-wide association study of sTNFR2 levels in the Framingham Heart Study Offspring cohort (examination 8; 2005-2008). sTNFR2 was quantitated by enzyme immunoassay and DNA methylation by microarray. The concentration of sTNFR2 was log_e-transformed and outliers were excluded. We conducted linear mixed effects models to test the association between sTNFR2 level and methylation at over 400,000 CpGs, adjusting for age, sex, BMI, smoking, imputed cell count, technical covariates, and accounting for familial relatedness. Results: The study sample included 2468 participants (mean age: 67 ± 9 years, 52% women, mean sTNFR2 level 2661 ± 1078 pg/ml). After accounting for multiple testing, we identified 168 CpGs (P < 1.2 × 10^-7) that were differentially methylated in relation to sTNFR2. A substantial proportion (27 CpGs; 16%) are in the major histocompatibility complex region and in loci overrepresented for antigen binding molecular functions (P = 1.7 × 10^-4) and antigen processing and presentation biological processes (P = 1.3 × 10^-8). Identified CpGs are enriched in active regulatory regions and associated with expression of 48 cis-genes (±500 kb) in whole blood (P < 1.1 × 10^-5) that coincide with genes identified in GWAS of diseases of immune dysregulation (inflammatory bowel disease, type 1 diabetes, IgA nephropathy). Conclusion: Differentially methylated loci in leukocytes associated with sTNF2 levels reside in active regulatory regions, are overrepresented in antigen processes, and are linked to inflammatory diseases.

Micro RNAs from DNA Viruses are Found Widely in Plasma in a Large Observational Human Population

Viral infections associate with disease risk and select families of viruses encode miRNAs that control an efficient viral cycle. The association of viral miRNA expression with disease in a large human population has not been previously explored. We sequenced plasma RNA from 40 participants of the Framingham Heart Study (FHS, Offspring Cohort, Visit 8) and identified 3 viral miRNAs from 3 different human Herpesviridae. These miRNAs were mostly related to viral latency and have not been previously detected in human plasma. Viral miRNA expression was then screened in the plasma of 2763 participants of the remaining cohort utilizing high-throughput RT-qPCR. All 3 viral miRNAs associated with combinations of inflammatory or prothrombotic circulating biomarkers (sTNFRII, IL-6, sICAM1, OPG, P-selectin) but did not associate with hypertension, coronary heart disease or cancer. Using a large observational population, we demonstrate that the presence of select viral miRNAs in the human circulation associate with inflammatory biomarkers and possibly immune response, but fail to associate with overt disease. This study greatly extends smaller singular observations of viral miRNAs in the human circulation and suggests that select viral miRNAs, such as those for latency, may not impact disease manifestation.

Age-associated microRNA expression in human peripheral blood is associated with all-cause mortality and age-related traits

Recent studies provide evidence of correlations of DNA methylation and expression of protein‐coding genes with human aging. The relations of microRNA expression with age and age‐related clinical outcomes have not been characterized thoroughly. We explored associations of age with whole‐blood microRNA expression in 5221 adults and identified 127 microRNAs that were differentially expressed by age at P < 3.3 × 10⁻⁴ (Bonferroni‐corrected). Most microRNAs were underexpressed in older individuals. Integrative analysis of microRNA and mRNA expression revealed changes in age‐associated mRNA expression possibly driven by age‐associated microRNAs in pathways that involve RNA processing, translation, and immune function. We fitted a linear model to predict ‘microRNA age’ that incorporated expression levels of 80 microRNAs. MicroRNA age correlated modestly with predicted age from DNA methylation (r = 0.3) and mRNA expression (r = 0.2), suggesting that microRNA age may complement mRNA and epigenetic age prediction models. We used the difference between microRNA age and chronological age as a biomarker of accelerated aging (Δage) and found that Δage was associated with all‐cause mortality (hazards ratio 1.1 per year difference, P = 4.2 × 10⁻⁵ adjusted for sex and chronological age). Additionally, Δage was associated with coronary heart disease, hypertension, blood pressure, and glucose levels. In conclusion, we constructed a microRNA age prediction model based on whole‐blood microRNA expression profiling. Age‐associated microRNAs and their targets have potential utility to detect accelerated aging and to predict risks for age‐related diseases.

A whole-blood transcriptome meta-analysis identifies gene expression signatures of cigarette smoking

Cigarette smoking is a leading modifiable cause of death worldwide. We hypothesized that cigarette smoking induces extensive transcriptomic changes that lead to target-organ damage and smoking-related diseases. We performed a meta-analysis of transcriptome-wide gene expression using whole blood-derived RNA from 10,233 participants of European ancestry in six cohorts (including 1421 current and 3955 former smokers) to identify associations between smoking and altered gene expression levels. At a false discovery rate (FDR) <0.1, we identified 1270 differentially expressed genes in current vs. never smokers, and 39 genes in former vs. never smokers. Expression levels of 12 genes remained elevated up to 30 years after smoking cessation, suggesting that the molecular consequence of smoking may persist for decades. Gene ontology analysis revealed enrichment of smoking-related genes for activation of platelets and lymphocytes, immune response, and apoptosis. Many of the top smoking-related differentially expressed genes, including LRRN3 and GPR15, have DNA methylation loci in promoter regions that were recently reported to be hypomethylated among smokers. By linking differential gene expression with smoking-related disease phenotypes, we demonstrated that stroke and pulmonary function show enrichment for smoking-related gene expression signatures. Mediation analysis revealed the expression of several genes (e.g. ALAS2) to be putative mediators of the associations between smoking and inflammatory biomarkers (IL6 and C-reactive protein levels). Our transcriptomic study provides potential insights into the effects of cigarette smoking on gene expression in whole blood and their relations to smoking-related diseases. The results of such analyses may highlight attractive targets for treating or preventing smoking-related health effects.

The Expressed Genome in Cardiovascular Diseases and Stroke: Refinement, Diagnosis, and Prediction: A Scientific Statement From the American Heart Association

There have been major advances in our knowledge of the contribution of DNA sequence variations to cardiovascular disease and stroke. However, the inner workings of the body reflect the complex interplay of factors beyond the DNA sequence, including epigenetic modifications, RNA transcripts, proteins, and metabolites, which together can be considered the "expressed genome." The emergence of high-throughput technologies, including epigenomics, transcriptomics, proteomics, and metabolomics, is now making it possible to address the contributions of the expressed genome to cardiovascular disorders. This statement describes how the expressed genome can currently and, in the future, potentially be used to diagnose diseases and to predict who will develop diseases such as coronary artery disease, stroke, heart failure, and arrhythmias.

Network Analysis of Gene Expression

Studies have pointed out that the expression of genes are highly regulated, which result in a cascade of distinct patterns of coexpression forming a network. Identifying and understanding such patterns is crucial in deciphering molecular mechanisms that underlie the pathophysiology of diseases. With the advance of high throughput assay of messenger RNA (mRNA) and high performance computing, reconstructing such network from molecular data such as gene expression is now possible. This chapter discusses an overview of methods of constructing such networks, practical considerations, and an example.

Transcriptome-wide association study of inflammatory biologic age

Chronic low grade inflammation is a fundamental mechanism of aging. We estimated biologic age using nine biomarkers from diverse inflammatory pathways and we hypothesized that genes associated with inflammatory biological age would provide insights into human aging. In Framingham Offspring Study participants at examination 8 (2005 to 2008), we used the Klemera-Doubal method to estimate inflammatory biologic age and we computed the difference (∆Age) between biologic age and chronologic age. Gene expression in whole blood was measured using the Affymetrix Human Exon 1.0 ST Array. We used linear mixed effect models to test associations between inflammatory ∆Age and gene expression (dependent variable) adjusting for age, sex, imputed cell counts, and technical covariates. Our study sample included 2386 participants (mean age 67A±9 years, 55% women). There were 448 genes significantly were associated with inflammatory ∆Age (P<2.8x10^-6), 302 genes were positively associated and 146 genes were negatively associated. Pathway analysis among the identified genes highlighted the NOD-like receptor signaling and ubiquitin mediated proteolysis pathways. In summary, we identified 448 genes that were significantly associated with inflammatory biologic age. Future functional characterization may identify molecular interventions to delay aging and prolong healthspan in older adults.

High expression of long chain acyl-coenzyme A synthetase 1 in peripheral blood may be a molecular marker for assessing the risk of acute myocardial infarction

The current study aimed to investigate whether the increased expression of long chain acyl-coenzymeA synthetase 1 (ACSL1) in peripheral blood leukocytes (PBL) may be a molecular marker for the genetic evaluation of acute myocardial infarction (AMI). The mechanism of action of ACSL1 in the pathogenesis of AMI was also investigated. A total of 75 patients with AMI and 70 individuals without coronary heart disease were selected to participate in the present study. The demographic and clinical information of the enrolled subjects was recorded. Reverse transcription quantitative polymerase chain reaction and western blot analysis were applied to measure the expression of ACSL1 at the mRNA and protein levels. It was demonstrated that the expression of ACSL1 mRNA and protein in PBL was increased in patients with AMI compared with controls. Logistic regression analysis indicated that ACSL1 expression in PBL was an independent risk factor of AMI. There was a significant positive association between the level of ACSL1 expression and the degree of atherosclerosis in the coronary artery. Furthermore, patients with AMI exhibited an increased risk of atherosclerosis due to increased fasting blood glucose, total cholesterol, triglyceride and lipoprotein levels and decreased high-density lipoprotein levels, compared with controls. Therefore, the current study demonstrated that ACSL1 expression was increased in the PBLs of patients with AMI. The elevated expression of ACSL1 acts an independent risk factor of AMI and may act as a potential biomarker when determining the risk of AMI.

MicroRNA Signature of Cigarette Smoking and Evidence for a Putative Causal Role of MicroRNAs in Smoking-Related Inflammation and Target Organ Damage

BACKGROUND

Cigarette smoking increases risk for multiple diseases. MicroRNAs regulate gene expression and may play a role in smoking-induced target organ damage. We sought to describe a microRNA signature of cigarette smoking and relate it to smoking-associated clinical phenotypes, gene expression, and lung inflammatory signaling.

METHODS AND RESULTS

Expression profiling of 283 microRNAs was conducted on whole blood-derived RNA from 5023 Framingham Heart Study participants (54.0% women; mean age, 55±13 years) using TaqMan assays and high-throughput reverse transcription quantitative polymerase chain reaction. Associations of microRNA expression with smoking status and associations of smoking-related microRNAs with inflammatory biomarkers and pulmonary function were tested with linear mixed effects models. We identified a 6-microRNA signature of smoking. Five of the 6 smoking-related microRNAs were associated with serum levels of C-reactive protein or interleukin-6; miR-1180 was associated with pulmonary function measures at a marginally significant level. Bioinformatic evaluation of smoking-associated genes coexpressed with the microRNA signature of cigarette smoking revealed enrichment for immune-related pathways. Smoking-associated microRNAs altered expression of selected inflammatory mediators in cell culture gain-of-function assays.

CONCLUSIONS

We characterized a novel microRNA signature of cigarette smoking. The top microRNAs were associated with systemic inflammatory markers and reduced pulmonary function, correlated with expression of genes involved in immune function, and were sufficient to modulate inflammatory signaling. Our results highlight smoking-associated microRNAs and are consistent with the hypothesis that smoking-associated microRNAs serve as mediators of smoking-induced inflammation and target organ damage. These findings call for further mechanistic studies to explore the diagnostic and therapeutic use of smoking-related microRNAs.