Correlation among traits is a fundamental feature of biological systems that remains difficult to study. To address this problem, we developed a flexible approach that allows us to identify factors associated with inter-individual variation in correlation. We use data from three human cohorts to study the effects of genetic and environmental variation on correlations among mRNA transcripts and among NMR metabolites. We first show that environmental exposures (infection and disease) lead to a systematic loss of correlation, which we define as 'decoherence'. Using longitudinal data, we show that decoherent metabolites are better predictors of whether someone will develop metabolic syndrome than metabolites commonly used as biomarkers of this disease. Finally, we demonstrate that correlation itself is under genetic control by mapping hundreds of 'correlation quantitative trait loci (QTLs)'. Together, this work furthers our understanding of how and why coordinated biological processes break down, and points to a potential role for decoherence in disease.

Editorial note

This article has been through an editorial process in which the authors decide how to respond to the issues raised during peer review. The Reviewing Editor's assessment is that all the issues have been addressed (see decision letter).

Background Recent studies have revealed sexually dimorphic associations between the carbamoyl-phosphate synthase 1 locus, intermediates of the metabolic pathway leading from choline to urea, and risk of coronary artery disease ( CAD ) in women. Based on evidence from the literature, the atheroprotective association with carbamoyl-phosphate synthase 1 could be mediated by the strong genetic effect of this locus on increased circulating glycine levels. Methods and Results We sought to identify additional genetic determinants of circulating glycine levels by carrying out a meta-analysis of genome-wide association study data in up to 30 118 subjects of European ancestry. Mendelian randomization and other analytical approaches were used to determine whether glycine-associated variants were associated with CAD and traditional risk factors. Twelve loci were significantly associated with circulating glycine levels, 7 of which were not previously known to be involved in glycine metabolism ( ACADM , PHGDH , COX 18- ADAMTS 3, PSPH , TRIB 1, PTPRD , and ABO ). Glycine-raising alleles at several loci individually exhibited directionally consistent associations with decreased risk of CAD . However, these effects could not be attributed directly to glycine because of associations with other CAD -related traits. By comparison, genetic models that only included the 2 variants directly involved in glycine degradation and for which there were no other pleiotropic associations were not associated with risk of CAD or blood pressure, lipid levels, and obesity-related traits. Conclusions These results provide additional insight into the genetic architecture of glycine metabolism, but do not yield conclusive evidence for a causal relationship between circulating levels of this amino acid and risk of CAD in humans.

Many genetic variants influence complex traits by modulating gene expression, thus altering the abundance of one or multiple proteins. Here we introduce a powerful strategy that integrates gene expression measurements with summary association statistics from large-scale genome-wide association studies (GWAS) to identify genes whose cis-regulated expression is associated with complex traits. We leverage expression imputation from genetic data to perform a transcriptome-wide association study (TWAS) to identify significant expression-trait associations. We applied our approaches to expression data from blood and adipose tissue measured in ∼ 3,000 individuals overall. We imputed gene expression into GWAS data from over 900,000 phenotype measurements to identify 69 new genes significantly associated with obesity-related traits (BMI, lipids and height). Many of these genes are associated with relevant phenotypes in the Hybrid Mouse Diversity Panel. Our results showcase the power of integrating genotype, gene expression and phenotype to gain insights into the genetic basis of complex traits.

Dyslipidemia and obesity are especially prevalent in populations with Amerindian backgrounds, such as Mexican-Americans, which predispose these populations to cardiovascular disease. Here we design an approach, known as the cross-population allele screen (CPAS), which we conduct prior to a genome-wide association study (GWAS) in 19,273 Europeans and Mexicans, in order to identify Amerindian risk genes in Mexicans. Utilizing CPAS to restrict the GWAS input variants to only those differing in frequency between the two populations, we identify novel Amerindian lipid genes, receptor-related orphan receptor alpha (RORA) and salt-inducible kinase 3 (SIK3), and three loci previously unassociated with dyslipidemia or obesity. We also detect lipoprotein lipase (LPL) and apolipoprotein A5 (APOA5) harbouring specific Amerindian signatures of risk variants and haplotypes. Notably, we observe that SIK3 and one novel lipid locus underwent positive selection in Mexicans. Furthermore, after a high-fat meal, the SIK3 risk variant carriers display high triglyceride levels. These findings suggest that Amerindian-specific genetic architecture leads to a higher incidence of dyslipidemia and obesity in modern Mexicans.

Implicating particular genes in the generation of complex brain and behavior phenotypes requires multiple lines of evidence. The rarity of most high-impact genetic variants typically precludes the possibility of accruing statistical evidence that they are associated with a given trait. We found that the enrichment of a rare chromosome 22q11.22 deletion in a recently expanded Northern Finnish sub-isolate enabled the detection of association between TOP3B and both schizophrenia and cognitive impairment. Biochemical analysis of TOP3β revealed that this topoisomerase was a component of cytosolic messenger ribonucleoproteins (mRNPs) and was catalytically active on RNA. The recruitment of TOP3β to mRNPs was independent of RNA cis-elements and was coupled to the co-recruitment of FMRP, the disease gene product in fragile X mental retardation syndrome. Our results indicate a previously unknown role for TOP3β in mRNA metabolism and suggest that it is involved in neurodevelopmental disorders.

Dyslipidemia and obesity are especially prevalent in populations with Amerindian backgrounds, such as Mexican-Americans, which predispose these populations to cardiovascular disease. Here we design an approach, known as the cross-population allele screen (CPAS), which we conduct prior to a genome-wide association study (GWAS) in 19,273 Europeans and Mexicans, in order to identify Amerindian risk genes in Mexicans. Utilizing CPAS to restrict the GWAS input variants to only those differing in frequency between the two populations, we identify novel Amerindian lipid genes, receptor-related orphan receptor alpha (RORA) and salt-inducible kinase 3 (SIK3), and three loci previously unassociated with dyslipidemia or obesity. We also detect lipoprotein lipase (LPL) and apolipoprotein A5 (APOA5) harbouring specific Amerindian signatures of risk variants and haplotypes. Notably, we observe that SIK3 and one novel lipid locus underwent positive selection in Mexicans. Furthermore, after a high-fat meal, the SIK3 risk variant carriers display high triglyceride levels. These findings suggest that Amerindian-specific genetic architecture leads to a higher incidence of dyslipidemia and obesity in modern Mexicans.

Background

There is limited information on genetic factors associated with sudden cardiac arrest (SCA).

Objective

To assess the association of common variation in genes in fatty acid pathways with SCA risk.

Methods

We selected 85 candidate genes and 1155 single nucleotide polymorphisms (SNPs) tagging common variation in each gene. We investigated the SNP associations with SCA in a population-based case-control study. Cases (n = 2160) were from a repository of SCA in the greater Seattle area. Controls (n = 2615), frequency-matched on age and sex, were from the same area. We used linear logistic regression to examine SNP associations with SCA. We performed permutation-based p-min tests to account for multiple comparisons within each gene. The SNP associations with a corrected P value of <.05 were then examined in a meta-analysis of these SNP associations in 9 replication studies totaling 2129 SCA cases and 23,833 noncases.

Results

Eight SNPs in or near 8 genes were associated with SCA risk in the discovery study, one of which was nominally significant in the replication phase (rs7737692, minor allele frequency 36%, near the LPCAT1 gene). For each copy of the minor allele, rs7737692 was associated with 13% lower SCA risk (95% confidence interval -21% to -5%) in the discovery phase and 9% lower SCA risk (95% confidence interval -16% to -1%) in the replication phase.

Conclusions

While none of the associations reached significance with Bonferroni correction, a common genetic variant near LPCAT1, a gene involved in the remodeling of phospholipids, was nominally associated with incident SCA risk. Further study is needed to validate this observation.

Aims

While most patients with myocardial infarction (MI) have underlying coronary atherosclerosis, not all patients with coronary artery disease (CAD) develop MI. We sought to address the hypothesis that some of the genetic factors which establish atherosclerosis may be distinct from those that predispose to vulnerable plaques and thrombus formation.

Methods and results

We carried out a genome-wide association study for MI in the UK Biobank (n∼472 000), followed by a meta-analysis with summary statistics from the CARDIoGRAMplusC4D Consortium (n∼167 000). Multiple independent replication analyses and functional approaches were used to prioritize loci and evaluate positional candidate genes. Eight novel regions were identified for MI at the genome wide significance level, of which effect sizes at six loci were more robust for MI than for CAD without the presence of MI. Confirmatory evidence for association of a locus on chromosome 1p21.3 harbouring choline-like transporter 3 (SLC44A3) with MI in the context of CAD, but not with coronary atherosclerosis itself, was obtained in Biobank Japan (n∼165 000) and 16 independent angiography-based cohorts (n∼27 000). Follow-up analyses did not reveal association of the SLC44A3 locus with CAD risk factors, biomarkers of coagulation, other thrombotic diseases, or plasma levels of a broad array of metabolites, including choline, trimethylamine N-oxide, and betaine. However, aortic expression of SLC44A3 was increased in carriers of the MI risk allele at chromosome 1p21.3, increased in ischaemic (vs. non-diseased) coronary arteries, up-regulated in human aortic endothelial cells treated with interleukin-1β (vs. vehicle), and associated with smooth muscle cell migration in vitro.

Conclusions

A large-scale analysis comprising ∼831 000 subjects revealed novel genetic determinants of MI and implicated SLC44A3 in the pathophysiology of vulnerable plaques.

Background

Short sleep duration has been associated with greater risks of obesity, hypertension, diabetes, and cardiovascular disease. Also, common genetic variants in the human Circadian Locomotor Output Cycles Kaput (CLOCK) show associations with ghrelin and total energy intake.

Objectives

We examined associations between habitual sleep duration, body mass index (BMI), and macronutrient intake and assessed whether CLOCK variants modify these associations.

Design

We conducted inverse-variance weighted, fixed-effect meta-analyses of results of adjusted associations of sleep duration and BMI and macronutrient intake as percentages of total energy as well as interactions with CLOCK variants from 9 cohort studies including up to 14,906 participants of European descent from the Cohorts for Heart and Aging Research in Genomic Epidemiology Consortium.

Results

We observed a significant association between sleep duration and lower BMI (β ± SE = 0.16 ± 0.04, P < 0.0001) in the overall sample; however, associations between sleep duration and relative macronutrient intake were evident in age- and sex-stratified analyses only. We observed a significant association between sleep duration and lower saturated fatty acid intake in younger (aged 20-64 y) adults (men: 0.11 ± 0.06%, P = 0.03; women: 0.10 ± 0.05%, P = 0.04) and with lower carbohydrate (-0.31 ± 0.12%, P < 0.01), higher total fat (0.18 ± 0.09%, P = 0.05), and higher PUFA (0.05 ± 0.02%, P = 0.02) intakes in older (aged 65-80 y) women. In addition, the following 2 nominally significant interactions were observed: between sleep duration and rs12649507 on PUFA intake and between sleep duration and rs6858749 on protein intake.

Conclusions

Our results indicate that longer habitual sleep duration is associated with lower BMI and age- and sex-specific favorable dietary behaviors. Differences in the relative intake of specific macronutrients associated with short sleep duration could, at least in part, explain previously reported associations between short sleep duration and chronic metabolic abnormalities. In addition, the influence of obesity-associated CLOCK variants on the association between sleep duration and macronutrient intake suggests that longer habitual sleep duration could ameliorate the genetic predisposition to obesity via a favorable dietary profile.

Rheumatoid arthritis (RA) and osteoporosis (OP) are two comorbid complex inflammatory conditions with evidence of shared genetic background and causal relationships. We aimed to clarify the genetic architecture underlying RA and various OP phenotypes while additionally considering an inflammatory component, C-reactive protein (CRP). Genome-wide association study summary statistics were acquired from the GEnetic Factors for OSteoporosis Consortium, Cohorts for Heart and Aging Research Consortium and UK Biobank. Mendelian randomization (MR) was used to detect the presence of causal relationships. Colocalization analysis was performed to determine shared genetic variants between CRP and OP phenotypes. Analysis of pleiotropy between traits owing to shared causal single nucleotide polymorphisms (SNPs) was performed using PL eiotropic A nalysis under CO mposite null hypothesis (PLACO). MR analysis was suggestive of horizontal pleiotropy between RA and OP traits. RA was a significant causal risk factor for CRP (β = 0.027, 95% confidence interval = 0.016-0.038). There was no evidence of CRP→OP causal relationship, but horizontal pleiotropy was apparent. Colocalization established shared genomic regions between CRP and OP, including GCKR and SERPINA1 genes. Pleiotropy arising from shared causal SNPs revealed through the colocalization analysis was all confirmed by PLACO. These genes were found to be involved in the same molecular function 'protein binding' (GO:0005515) associated with RA, OP and CRP. We identified three major components explaining the epidemiological relationship among RA, OP and inflammation: (1) Pleiotropy explains a portion of the shared genetic relationship between RA and OP, albeit polygenically; (2) RA contributes to CRP elevation and (3) CRP, which is influenced by RA, demonstrated pleiotropy with OP.