Genome-Wide Association with Select Biomarker Traits in the Framingham Heart Study

BACKGROUND: Systemic biomarkers provide insights into disease pathogenesis, diagnosis, and risk stratification. Many systemic biomarker concentrations are heritable phenotypes. Genome-wide association studies (GWAS) provide mechanisms to investigate the genetic contributions to biomarker variability unconstrained by current knowledge of physiological relations. METHODS: We examined the association of Affymetrix 100K GeneChip single nucleotide polymorphisms (SNPs) to 22 systemic biomarker concentrations in 4 biological domains: inflammation/oxidative stress; natriuretic peptides; liver function; and vitamins. Related members of the Framingham Offspring cohort (n = 1012; mean age 59 ± 10 years, 51% women) had both phenotype and genotype data (minimum-maximum per phenotype n = 507–1008). We used Generalized Estimating Equations (GEE), Family Based Association Tests (FBAT) and variance components linkage to relate SNPs to multivariable-adjusted biomarker residuals. Autosomal SNPs (n = 70,987) meeting the following criteria were studied: minor allele frequency ≥ 10%, call rate ≥ 80% and Hardy-Weinberg equilibrium p ≥ 0.001. RESULTS: With GEE, 58 SNPs had p < 10-6: the top SNPs were rs2494250 (p = 1.00*10-14) and rs4128725 (p = 3.68*10-12) for monocyte chemoattractant protein-1 (MCP1), and rs2794520 (p = 2.83*10-8) and rs2808629 (p = 3.19*10-8) for C-reactive protein (CRP) averaged from 3 examinations (over about 20 years). With FBAT, 11 SNPs had p < 10-6: the top SNPs were the same for MCP1 (rs4128725, p = 3.28*10-8, and rs2494250, p = 3.55*10-8), and also included B-type natriuretic peptide (rs437021, p = 1.01*10-6) and Vitamin K percent undercarboxylated osteocalcin (rs2052028, p = 1.07*10-6). The peak LOD (logarithm of the odds) scores were for MCP1 (4.38, chromosome 1) and CRP (3.28, chromosome 1; previously described) concentrations; of note the 1.5 support interval included the MCP1 and CRP SNPs reported above (GEE model). Previous candidate SNP associations with circulating CRP concentrations were replicated at p < 0.05; the SNPs rs2794520 and rs2808629 are in linkage disequilibrium with previously reported SNPs. GEE, FBAT and linkage results are posted at . CONCLUSION: The Framingham GWAS represents a resource to describe potentially novel genetic influences on systemic biomarker variability. The newly described associations will need to be replicated in other studies.National Heart, Lung, and Blood Institute's Framingham Heart Study (N01-HC25195); National Institutes of Health National Center for Research Resources Shared Instrumentation grant (1S10RR163736-01A1); National Institutes of Health (HL064753, HL076784, AG028321, HL71039, 2 K24HL04334, 1K23 HL083102); Doris Duke Charitable Foundation; American Diabetes Association Career Developement Award; National Center for Research Resources (GCRC M01-RR01066); US Department of Agriculture Agricultural Research Service (58-1950-001, 58-1950-401); National Institute of Aging (AG14759

Meta-analysis of genome-wide association studies from the CHARGE consortium identifies common variants associated with carotid intima media thickness and plaque

Carotid intima media thickness (cIMT) and plaque determined by ultrasonography are established measures of subclinical atherosclerosis that each predicts future cardiovascular disease events. We conducted a meta-analysis of genome-wide association data in 31,211 participants of European ancestry from nine large studies in the setting of the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium. We then sought additional evidence to support our findings among 11,273 individuals using data from seven additional studies. In the combined meta-analysis, we identified three genomic regions associated with common carotid intima media thickness and two different regions associated with the presence of carotid plaque (P < 5 × 10 -8). The associated SNPs mapped in or near genes related to cellular signaling, lipid metabolism and blood pressure homeostasis, and two of the regions were associated with coronary artery disease (P < 0.006) in the Coronary Artery Disease Genome-Wide Replication and Meta-Analysis (CARDIoGRAM) consortium. Our findings may provide new insight into pathways leading to subclinical atherosclerosis and subsequent cardiovascular events

Cross-Sectional Associations between Homoarginine, Intermediate Phenotypes, and Atrial Fibrillation in the CommunityThe Gutenberg Health Study

Homoarginine has come into the focus of interest as a biomarker for cardiovascular disease. Atrial fibrillation (AF) causes a substantial increase in morbidity and mortality. Whether circulating homoarginine is associated with occurrence or persistence of AF and may serve as a new predictive biomarker remains unknown. We measured plasma levels of homoarginine in the population-based Gutenberg health study (3761 patients included, of them 51.7% males), mean age 55.6 +/- 10.9 years-old. Associations between homoarginine and intermediate electrocardiographic and echocardiographic phenotypes and manifest AF were examined. Patients with AF (124 patients, of them 73.4% males) had a mean age 64.8 +/- 8.6 years-old compared to a mean age of 55.3 +/- 10.9 in the population without AF (p-value < 0.001) and showed a less beneficial risk factor profile. The median homoarginine levels in individuals with and without AF were 1.9 mol/L (interquartile range (IQR) 1.5-2.5) and 2.0 mol/L (IQR 1.5-2.5), respectively, p = 0.56. In multivariable-adjusted regression analyses homoarginine was not statistically significantly related to electrocardiographic variables. Among echocardiographic variables beta per standard deviation increase was -0.12 (95% confidence interval (CI) -0.23-(-0.02);p = 0.024) for left atrial area and -0.01 (95% CI -0.02-(-0.003);p = 0.013) for E/A ratio. The odds ratio between homoarginine and AF was 0.91 (95% CI 0.70-1.16;p = 0.45). In our large, population-based cross-sectional study, we did not find statistically significant correlations between lower homoarginine levels and occurrence or persistence of AF or most standard electrocardiographic phenotypes, but some moderate inverse associations with echocardiographic left atrial size and E/A. Homoarginine may not represent a strong biomarker to identify individuals at increased risk for AF. Further investigations will be needed to elucidate the role of homoarginine and cardiac function

Anticoagulation with edoxaban in patients with long Atrial High-Rate Episodes ≥24 hours

BACKGROUND AND AIMS: Patients with long atrial high-rate episodes (AHRE) ≥ 24 hours and stroke risk factors are often treated with anticoagulation for stroke prevention. Anticoagulation has never been compared to no anticoagulation in these patients.METHODS: This secondary prespecified analysis of NOAH-AFNET 6 examined interactions between AHRE duration at baseline and anticoagulation with edoxaban compared to placebo in patients with AHRE and stroke risk factors. The primary efficacy outcome was a composite of stroke, systemic embolism, or cardiovascular death. The safety outcome was a composite of major bleeding and death. Key secondary outcomes were components of these outcomes and ECG-diagnosed atrial fibrillation.RESULTS: AHRE ≥24 hours were present at baseline in 259/2389 patients enrolled in NOAH-AFNET 6 (11%, 78 ± 7 years old, 28% women, CHA2DS2-VASc score 4). Clinical characteristics were not different from patients with shorter AHRE. During a median follow-up of 1.8 years, the primary outcome occurred in 9/132 patients with AHRE ≥24 hours (4.3%/patient-year, 2 strokes) treated with anticoagulation and in 14/127 patients treated with placebo (6.9%/patient-year, 2 strokes). AHRE duration did not interact with the efficacy (p-interaction = 0.65) or safety (p-interaction = 0.98) of anticoagulation. Analyses including AHRE as a continuous parameter confirmed this. Patients with AHRE ≥24 hours developed more ECG-diagnosed atrial fibrillation (17.0%/patient-year) than patients with shorter AHRE (8.2%/patient-year; p &lt; 0.001).CONCLUSIONS: This hypothesis-generating analysis does not find an interaction between AHRE duration and anticoagulation therapy in patients with device-detected AHRE and stroke risk factors. Further research is needed to identify patients with long AHRE at high stroke risk.</p

An arrhythmogenic metabolite in atrial fibrillation

Abstract Background Long-chain acyl-carnitines (ACs) are potential arrhythmogenic metabolites. Their role in atrial fibrillation (AF) remains incompletely understood. Using a systems medicine approach, we assessed the contribution of C18:1AC to AF by analysing its in vitro effects on cardiac electrophysiology and metabolism, and translated our findings into the human setting. Methods and results Human iPSC-derived engineered heart tissue was exposed to C18:1AC. A biphasic effect on contractile force was observed: short exposure enhanced contractile force, but elicited spontaneous contractions and impaired Ca2+ handling. Continuous exposure provoked an impairment of contractile force. In human atrial mitochondria from AF individuals, C18:1AC inhibited respiration. In a population-based cohort as well as a cohort of patients, high C18:1AC serum concentrations were associated with the incidence and prevalence of AF. Conclusion Our data provide evidence for an arrhythmogenic potential of the metabolite C18:1AC. The metabolite interferes with mitochondrial metabolism, thereby contributing to contractile dysfunction and shows predictive potential as novel circulating biomarker for risk of AF

RANTES/CCL5 and risk for coronary events:Results from the MONICA/KORA Augsburg case-cohort, Athero-express and CARDIoGRAM studies

Background: The chemokine RANTES (regulated on activation, normal T-cell expressed and secreted)/CCL5 is involved in the pathogenesis of cardiovascular disease in mice, whereas less is known in humans. We hypothesised that its relevance for atherosclerosis should be reflected by associations between CCL5 gene variants, RANTES serum concentrations and protein levels in atherosclerotic plaques and risk for coronary events. Methods and Findings: We conducted a case-cohort study within the population-based MONICA/KORA Augsburg studies. Baseline RANTES serum levels were measured in 363 individuals with incident coronary events and 1,908 non-cases (mean follow-up: 10.2±4.8 years). Cox proportional hazard models adjusting for age, sex, body mass index, metabolic factors and lifestyle factors revealed no significant association between RANTES and incident coronary events (HR [95% CI] for increasing RANTES tertiles 1.0, 1.03 [0.75-1.42] and 1.11 [0.81-1.54]). None of six CCL5 single nucleotide polymorphisms and no common haplotype showed significant associations with coronary events. Also in the CARDIoGRAM study (&gt;22,000 cases, &gt;60,000 controls), none of these CCL5 SNPs was significantly associated with coronary artery disease. In the prospective Athero-Express biobank study, RANTES plaque levels were measured in 606 atherosclerotic lesions from patients who underwent carotid endarterectomy. RANTES content in atherosclerotic plaques was positively associated with macrophage infiltration and inversely associated with plaque calcification. However, there was no significant association between RANTES content in plaques and risk for coronary events (mean follow-up 2.8±0.8 years). Conclusions: High RANTES plaque levels were associated with an unstable plaque phenotype. However, the absence of associations between (i) RANTES serum levels, (ii) CCL5 genotypes and (iii) RANTES content in carotid plaques and either coronary artery disease or incident coronary events in our cohorts suggests that RANTES may not be a novel coronary risk biomarker. However, the potential relevance of RANTES levels in platelet-poor plasma needs to be investigated in further studies.</p

Genome-Wide Association Studies of the PR Interval in African Americans

The PR interval on the electrocardiogram reflects atrial and atrioventricular nodal conduction time. The PR interval is heritable, provides important information about arrhythmia risk, and has been suggested to differ among human races. Genome-wide association (GWA) studies have identified common genetic determinants of the PR interval in individuals of European and Asian ancestry, but there is a general paucity of GWA studies in individuals of African ancestry. We performed GWA studies in African American individuals from four cohorts (n = 6,247) to identify genetic variants associated with PR interval duration. Genotyping was performed using the Affymetrix 6.0 microarray. Imputation was performed for 2.8 million single nucleotide polymorphisms (SNPs) using combined YRI and CEU HapMap phase II panels. We observed a strong signal (rs3922844) within the gene encoding the cardiac sodium channel (SCN5A) with genome-wide significant association (p<2.5×10−8) in two of the four cohorts and in the meta-analysis. The signal explained 2% of PR interval variability in African Americans (beta  = 5.1 msec per minor allele, 95% CI  = 4.1–6.1, p = 3×10−23). This SNP was also associated with PR interval (beta = 2.4 msec per minor allele, 95% CI = 1.8–3.0, p = 3×10−16) in individuals of European ancestry (n = 14,042), but with a smaller effect size (p for heterogeneity <0.001) and variability explained (0.5%). Further meta-analysis of the four cohorts identified genome-wide significant associations with SNPs in SCN10A (rs6798015), MEIS1 (rs10865355), and TBX5 (rs7312625) that were highly correlated with SNPs identified in European and Asian GWA studies. African ancestry was associated with increased PR duration (13.3 msec, p = 0.009) in one but not the other three cohorts. Our findings demonstrate the relevance of common variants to African Americans at four loci previously associated with PR interval in European and Asian samples and identify an association signal at one of these loci that is more strongly associated with PR interval in African Americans than in Europeans

Integrating Genome-Wide Genetic Variations and Monocyte Expression Data Reveals Trans-Regulated Gene Modules in Humans

One major expectation from the transcriptome in humans is to characterize the biological basis of associations identified by genome-wide association studies. So far, few cis expression quantitative trait loci (eQTLs) have been reliably related to disease susceptibility. Trans-regulating mechanisms may play a more prominent role in disease susceptibility. We analyzed 12,808 genes detected in at least 5% of circulating monocyte samples from a population-based sample of 1,490 European unrelated subjects. We applied a method of extraction of expression patterns—independent component analysis—to identify sets of co-regulated genes. These patterns were then related to 675,350 SNPs to identify major trans-acting regulators. We detected three genomic regions significantly associated with co-regulated gene modules. Association of these loci with multiple expression traits was replicated in Cardiogenics, an independent study in which expression profiles of monocytes were available in 758 subjects. The locus 12q13 (lead SNP rs11171739), previously identified as a type 1 diabetes locus, was associated with a pattern including two cis eQTLs, RPS26 and SUOX, and 5 trans eQTLs, one of which (MADCAM1) is a potential candidate for mediating T1D susceptibility. The locus 12q24 (lead SNP rs653178), which has demonstrated extensive disease pleiotropy, including type 1 diabetes, hypertension, and celiac disease, was associated to a pattern strongly correlating to blood pressure level. The strongest trans eQTL in this pattern was CRIP1, a known marker of cellular proliferation in cancer. The locus 12q15 (lead SNP rs11177644) was associated with a pattern driven by two cis eQTLs, LYZ and YEATS4, and including 34 trans eQTLs, several of them tumor-related genes. This study shows that a method exploiting the structure of co-expressions among genes can help identify genomic regions involved in trans regulation of sets of genes and can provide clues for understanding the mechanisms linking genome-wide association loci to disease

Genome-wide and gene-centric analyses of circulating myeloperoxidase levels in the charge and care consortia

Increased systemic levels of myeloperoxidase (MPO) are associated with the risk of coronary artery disease (CAD). To identify the genetic factors that are associated with circulating MPO levels, we carried out a genome-wide association study (GWAS) and a gene-centric analysis in subjects of European ancestry and African Americans (AAs). A locus on chromosome 1q31.1 containing the complement factor H (CFH) gene was strongly associated with serum MPO levels in 9305 subjects of European ancestry (lead SNP rs800292; P = 4.89 × 10−41) and in 1690 AA subjects (rs505102; P = 1.05 × 10−8). Gene-centric analyses in 8335 subjects of European ancestry additionally identified two rare MPO coding sequence variants that were associated with serum MPO levels (rs28730837, P = 5.21 × 10−12; rs35897051, P = 3.32 × 10−8). A GWAS for plasma MPO levels in 9260 European ancestry subjects identified a chromosome 17q22 region near MPO that was significantly associated (lead SNP rs6503905; P = 2.94 × 10−12), but the CFH locus did not exhibit evidence of association with plasma MPO levels. Functional analyses revealed that rs800292 was associated with levels of complement proteins in serum. Variants at chromosome 17q22 also had pleiotropic cis effects on gene expression. In a case–control analysis of ∼80 000 subjects from CARDIoGRAM, none of the identified single-nucleotide polymorphisms (SNPs) were associated with CAD. These results suggest that distinct genetic factors regulate serum and plasma MPO levels, which may have relevance for various acute and chronic inflammatory disorders. The clinical implications for CAD and a better understanding of the functional basis for the association of CFH and MPO variants with circulating MPO levels require further study