Paul Courchesne

Galectin-3, a marker of cardiac fibrosis, predicts incident heart failure in the community.

OBJECTIVES The aim of this study was to examine the relation of galectin-3 (Gal-3), a marker of cardiac fibrosis, with incident heart failure (HF) in the community. BACKGROUND Gal-3 is an emerging prognostic biomarker in HF, and experimental studies suggest that Gal-3 is an important mediator of cardiac fibrosis. Whether elevated Gal-3 concentrations precede the development of HF is unknown. METHODS Gal-3 concentrations were measured in 3,353 participants in the Framingham Offspring Cohort (mean age 59 years; 53% women). The relation of Gal-3 to incident HF was assessed using proportional hazards regression. RESULTS Gal-3 was associated with increased left ventricular mass in age-adjusted and sex-adjusted analyses (p = 0.001); this association was attenuated in multivariate analyses (p = 0.06). A total of 166 participants developed incident HF and 468 died during a mean follow-up period of 11.2 years. Gal-3 was associated with risk for incident HF (hazard ratio [HR]: 1.28 per 1 SD increase in log Gal-3; 95% confidence interval [CI]: 1.14 to 1.43; p < 0.0001) and remained significant after adjustment for clinical variables and B-type natriuretic peptide (HR: 1.23; 95% CI: 1.04 to 1.47; p = 0.02). Gal-3 was also associated with risk for all-cause mortality (multivariable-adjusted HR: 1.15; 95% CI: 1.04 to 1.28; p = 0.01). The addition of Gal-3 to clinical factors resulted in negligible changes to the C-statistic and minor improvements in net reclassification improvement. CONCLUSIONS Higher concentration of Gal-3, a marker of cardiac fibrosis, is associated with increased risk for incident HF and mortality. Future studies evaluating the role of Gal-3 in cardiac remodeling may provide further insights into the role of Gal-3 in the pathophysiology of HF.

A Systems Biology Framework Identifies Molecular Underpinnings of Coronary Heart Disease

Objective—Genetic approaches have identified numerous loci associated with coronary heart disease (CHD). The molecular mechanisms underlying CHD gene–disease associations, however, remain unclear. We hypothesized that genetic variants with both strong and subtle effects drive gene subnetworks that in turn affect CHD. Approach and Results—We surveyed CHD-associated molecular interactions by constructing coexpression networks using whole blood gene expression profiles from 188 CHD cases and 188 age- and sex-matched controls. Twenty-four coexpression modules were identified, including 1 case-specific and 1 control-specific differential module (DM). The DMs were enriched for genes involved in B-cell activation, immune response, and ion transport. By integrating the DMs with gene expression–associated single-nucleotide polymorphisms and with results of genome-wide association studies of CHD and its risk factors, the control-specific DM was implicated as CHD causal based on its significant enrichment for both CHD and lipid expression–associated single-nucleotide polymorphisms. This causal DM was further integrated with tissue-specific Bayesian networks and protein–protein interaction networks to identify regulatory key driver genes. Multitissue key drivers (SPIB and TNFRSF13C) and tissue-specific key drivers (eg, EBF1) were identified. Conclusions—Our network-driven integrative analysis not only identified CHD-related genes, but also defined network structure that sheds light on the molecular interactions of genes associated with CHD risk.

Gene Expression Signatures of Coronary Heart Disease

Objective—To identify transcriptomic biomarkers of coronary heart disease (CHD) in 188 cases with CHD and 188 age- and sex-matched controls who were participants in the Framingham Heart Study. Approach and Results—A total of 35 genes were differentially expressed in cases with CHD versus controls at false discovery rate<0.5, including GZMB, TMEM56, and GUK1. Cluster analysis revealed 3 gene clusters associated with CHD, 2 linked to increased erythrocyte production and a third to reduced natural killer and T cell activity in cases with CHD. Exon-level results corroborated and extended the gene-level results. Alternative splicing analysis suggested that GUK1 and 38 other genes were differentially spliced in cases with CHD versus controls. Gene Ontology analysis linked ubiquitination and T-cell–related pathways with CHD. Conclusions—Two bioinformatically defined groups of genes show consistent associations with CHD. Our findings are consistent with the hypotheses that hematopoesis is upregulated in CHD, possibly reflecting a compensatory mechanism, and that innate immune activity is disrupted in CHD or altered by its treatment. Transcriptomic signatures may be useful in identifying pathways associated with CHD and point toward novel therapeutic targets for its treatment and prevention.

The Liver Toxicity Biomarker Study: Phase I Design and Preliminary Results

Drug-induced liver injury (DILI) is the primary adverse event that results in withdrawal of drugs from the market and a frequent reason for the failure of drug candidates in development. The Liver Toxicity Biomarker Study (LTBS) is an innovative approach to investigate DILI because it compares molecular events produced in vivo by compound pairs that (a) are similar in structure and mechanism of action, (b) are associated with few or no signs of liver toxicity in preclinical studies, and (c) show marked differences in hepatotoxic potential. The LTBS is a collaborative preclinical research effort in molecular systems toxicology between the National Center for Toxicological Research and BG Medicine, Inc., and is supported by seven pharmaceutical companies and three technology providers. In phase I of the LTBS, entacapone and tolcapone were studied in rats to provide results and information that will form the foundation for the design and implementation of phase II. Molecular analysis of the rat liver and plasma samples combined with statistical analyses of the resulting datasets yielded marker analytes, illustrating the value of the broad-spectrum, molecular systems analysis approach to studying pharmacological or toxicological effects.

A Meta-analysis of Gene Expression Signatures of Blood Pressure and Hypertension

Genome-wide association studies (GWAS) have uncovered numerous genetic variants (SNPs) that are associated with blood pressure (BP). Genetic variants may lead to BP changes by acting on intermediate molecular phenotypes such as coded protein sequence or gene expression, which in turn affect BP variability. Therefore, characterizing genes whose expression is associated with BP may reveal cellular processes involved in BP regulation and uncover how transcripts mediate genetic and environmental effects on BP variability. A meta-analysis of results from six studies of global gene expression profiles of BP and hypertension in whole blood was performed in 7017 individuals who were not receiving antihypertensive drug treatment. We identified 34 genes that were differentially expressed in relation to BP (Bonferroni-corrected p<0.05). Among these genes, FOS and PTGS2 have been previously reported to be involved in BP-related processes; the others are novel. The top BP signature genes in aggregate explain 5%–9% of inter-individual variance in BP. Of note, rs3184504 in SH2B3, which was also reported in GWAS to be associated with BP, was found to be a trans regulator of the expression of 6 of the transcripts we found to be associated with BP (FOS, MYADM, PP1R15A, TAGAP, S100A10, and FGBP2). Gene set enrichment analysis suggested that the BP-related global gene expression changes include genes involved in inflammatory response and apoptosis pathways. Our study provides new insights into molecular mechanisms underlying BP regulation, and suggests novel transcriptomic markers for the treatment and prevention of hypertension.

Integrative network analysis reveals molecular mechanisms of blood pressure regulation.

Genome‐wide association studies (GWAS) have identified numerous loci associated with blood pressure (BP). The molecular mechanisms underlying BP regulation, however, remain unclear. We investigated BP‐associated molecular mechanisms by integrating BP GWAS with whole blood mRNA expression profiles in 3,679 individuals, using network approaches. BP transcriptomic signatures at the single‐gene and the coexpression network module levels were identified. Four coexpression modules were identified as potentially causal based on genetic inference because expression‐related SNPs for their corresponding genes demonstrated enrichment for BP GWAS signals. Genes from the four modules were further projected onto predefined molecular interaction networks, revealing key drivers. Gene subnetworks entailing molecular interactions between key drivers and BP‐related genes were uncovered. As proof‐of‐concept, we validated SH2B3, one of the top key drivers, using Sh2b3−/− mice. We found that a significant number of genes predicted to be regulated by SH2B3 in gene networks are perturbed in Sh2b3−/− mice, which demonstrate an exaggerated pressor response to angiotensin II infusion. Our findings may help to identify novel targets for the prevention or treatment of hypertension.

Genome-wide Identification of microRNA Expression Quantitative Trait Loci

Identification of microRNA expression quantitative trait loci (miR-eQTL) can yield insights into regulatory mechanisms of microRNA transcription, and can help elucidate the role of microRNA as mediators of complex traits. Here we present a miR-eQTL mapping study of whole blood from 5239 individuals, and identify 5269 cis-miR-eQTLs for 76 mature microRNAs. Forty-nine percent of cis-miR-eQTLs are located 300–500kb upstream of their associated intergenic microRNAs, suggesting that distal regulatory elements may affect the interindividual variability in microRNA expression levels. We find that cis-miR-eQTLs are highly enriched for cis-mRNA-eQTLs and regulatory SNPs. Among 243 cis-miR-eQTLs that were reported to be associated with complex traits in prior genome-wide association studies, many cis-miR-eQTLs miRNAs display differential expression in relation to the corresponding trait (e.g., rs7115089, miR-125b-5p, and HDL cholesterol). Our study provides a roadmap for understanding the genetic basis of miRNA expression, and sheds light on miRNA involvement in a variety of complex traits.

Association of body mass index with DNA methylation and gene expression in blood cells and relations to cardiometabolic disease: a Mendelian randomization approach

Background The link between DNA methylation, obesity, and adiposity-related diseases in the general population remains uncertain. Methods and Findings We conducted an association study of body mass index (BMI) and differential methylation for over 400,000 CpGs assayed by microarray in whole-blood-derived DNA from 3,743 participants in the Framingham Heart Study and the Lothian Birth Cohorts, with independent replication in three external cohorts of 4,055 participants. We examined variations in whole blood gene expression and conducted Mendelian randomization analyses to investigate the functional and clinical relevance of the findings. We identified novel and previously reported BMI-related differential methylation at 83 CpGs that replicated across cohorts; BMI-related differential methylation was associated with concurrent changes in the expression of genes in lipid metabolism pathways. Genetic instrumental variable analysis of alterations in methylation at one of the 83 replicated CpGs, cg11024682 (intronic to sterol regulatory element binding transcription factor 1 [SREBF1]), demonstrated links to BMI, adiposity-related traits, and coronary artery disease. Independent genetic instruments for expression of SREBF1 supported the findings linking methylation to adiposity and cardiometabolic disease. Methylation at a substantial proportion (16 of 83) of the identified loci was found to be secondary to differences in BMI. However, the cross-sectional nature of the data limits definitive causal determination. Conclusions We present robust associations of BMI with differential DNA methylation at numerous loci in blood cells. BMI-related DNA methylation and gene expression provide mechanistic insights into the relationship between DNA methylation, obesity, and adiposity-related diseases.

Integromic Analysis of Genetic Variation and Gene Expression Identifies Networks for Cardiovascular Disease Phenotypes

Background— Cardiovascular disease (CVD) reflects a highly coordinated complex of traits. Although genome-wide association studies have reported numerous single nucleotide polymorphisms (SNPs) to be associated with CVD, the role of most of these variants in disease processes remains unknown. Methods and Results— We built a CVD network using 1512 SNPs associated with 21 CVD traits in genome-wide association studies (at P⩽5×10−8) and cross-linked different traits by virtue of their shared SNP associations. We then explored whole blood gene expression in relation to these SNPs in 5257 participants in the Framingham Heart Study. At a false discovery rate <0.05, we identified 370 cis–expression quantitative trait loci (eQTLs; SNPs associated with altered expression of nearby genes) and 44 trans-eQTLs (SNPs associated with altered expression of remote genes). The eQTL network revealed 13 CVD-related modules. Searching for association of eQTL genes with CVD risk factors (lipids, blood pressure, fasting blood glucose, and body mass index) in the same individuals, we found examples in which the expression of eQTL genes was significantly associated with these CVD phenotypes. In addition, mediation tests suggested that a subset of SNPs previously associated with CVD phenotypes in genome-wide association studies may exert their function by altering expression of eQTL genes (eg, LDLR and PCSK7), which in turn may promote interindividual variation in phenotypes. Conclusions— Using a network approach to analyze CVD traits, we identified complex networks of SNP-phenotype and SNP-transcript connections. Integrating the CVD network with phenotypic data, we identified biological pathways that may provide insights into potential drug targets for treatment or prevention of CVD.

Protein biomarkers of new-onset cardiovascular disease: prospective study from the systems approach to biomarker research in cardiovascular disease initiative.

Objective— Incorporation of novel plasma protein biomarkers may improve current models for prediction of atherosclerotic cardiovascular disease (ASCVD) risk. Approach and Results— We used discovery mass spectrometry (MS) to determine plasma concentrations of 861 proteins in 135 myocardial infarction (MI) cases and 135 matched controls. Then, we measured 59 markers by targeted MS in 336 ASCVD case–control pairs. Associations with MI or ASCVD were tested in single-marker and multiple-marker analyses adjusted for established ASCVD risk factors. Twelve single markers from discovery MS were associated with MI incidence (at P <0.01), adjusting for clinical risk factors. Seven proteins in aggregate (cyclophilin A, cluster of differentiation 5 molecule [CD5] antigen-like, cell-surface glycoprotein mucin cell surface associated protein 18 [MUC-18], collagen-α 1 [XVIII] chain, salivary α-amylase 1, C-reactive protein, and multimerin-2) were highly associated with MI ( P <0.0001) and significantly improved its prediction compared with a model with clinical risk factors alone (C-statistic of 0.71 versus 0.84). Through targeted MS, 12 single proteins were predictors of ASCVD (at P <0.05) after adjusting for established risk factors. In multiple-marker analyses, 4 proteins in combination (α-1–acid glycoprotein 1, paraoxonase 1, tetranectin, and CD5 antigen-like) predicted incident ASCVD ( P <0.0001) and moderately improved the C-statistic from the model with clinical covariates alone (C-statistic of 0.69 versus 0.73). Conclusions— Proteomics profiling identified single- and multiple-marker protein panels that are associated with new-onset ASCVD and may lead to a better understanding of underlying disease mechanisms. Our findings include many novel protein biomarkers that, if externally validated, may improve risk assessment for MI and ASCVD. # Significance {#article-title-29}Objective—Incorporation of novel plasma protein biomarkers may improve current models for prediction of atherosclerotic cardiovascular disease (ASCVD) risk. Approach and Results—We used discovery mass spectrometry (MS) to determine plasma concentrations of 861 proteins in 135 myocardial infarction (MI) cases and 135 matched controls. Then, we measured 59 markers by targeted MS in 336 ASCVD case–control pairs. Associations with MI or ASCVD were tested in single-marker and multiple-marker analyses adjusted for established ASCVD risk factors. Twelve single markers from discovery MS were associated with MI incidence (at P<0.01), adjusting for clinical risk factors. Seven proteins in aggregate (cyclophilin A, cluster of differentiation 5 molecule [CD5] antigen-like, cell-surface glycoprotein mucin cell surface associated protein 18 [MUC-18], collagen-&agr; 1 [XVIII] chain, salivary &agr;-amylase 1, C-reactive protein, and multimerin-2) were highly associated with MI (P<0.0001) and significantly improved its prediction compared with a model with clinical risk factors alone (C-statistic of 0.71 versus 0.84). Through targeted MS, 12 single proteins were predictors of ASCVD (at P<0.05) after adjusting for established risk factors. In multiple-marker analyses, 4 proteins in combination (&agr;-1–acid glycoprotein 1, paraoxonase 1, tetranectin, and CD5 antigen-like) predicted incident ASCVD (P<0.0001) and moderately improved the C-statistic from the model with clinical covariates alone (C-statistic of 0.69 versus 0.73). Conclusions—Proteomics profiling identified single- and multiple-marker protein panels that are associated with new-onset ASCVD and may lead to a better understanding of underlying disease mechanisms. Our findings include many novel protein biomarkers that, if externally validated, may improve risk assessment for MI and ASCVD.