Relationship between obesity and the risk of clinically significant depression: Mendelian randomisation study.

BACKGROUND: Obesity has been shown to be associated with depression and it has been suggested that higher body mass index (BMI) increases the risk of depression and other common mental disorders. However, the causal relationship remains unclear and Mendelian randomisation, a form of instrumental variable analysis, has recently been employed to attempt to resolve this issue. AIMS: To investigate whether higher BMI increases the risk of major depression. METHOD: Two instrumental variable analyses were conducted to test the causal relationship between obesity and major depression in RADIANT, a large case-control study of major depression. We used a single nucleotide polymorphism (SNP) in FTO and a genetic risk score (GRS) based on 32 SNPs with well-established associations with BMI. RESULTS: Linear regression analysis, as expected, showed that individuals carrying more risk alleles of FTO or having higher score of GRS had a higher BMI. Probit regression suggested that higher BMI is associated with increased risk of major depression. However, our two instrumental variable analyses did not support a causal relationship between higher BMI and major depression (FTO genotype: coefficient -0.03, 95% CI -0.18 to 0.13, P = 0.73; GRS: coefficient -0.02, 95% CI -0.11 to 0.07, P = 0.62). CONCLUSIONS: Our instrumental variable analyses did not support a causal relationship between higher BMI and major depression. The positive associations of higher BMI with major depression in probit regression analyses might be explained by reverse causality and/or residual confounding.

Uric Acid and Cardiovascular Events: A Mendelian Randomization Study

Obesity and diets rich in uric acid-raising components appear to account for the increased prevalence of hyperuricemia in Westernized populations. Prevalence rates of hypertension, diabetes mellitus, CKD, and cardiovascular disease are also increasing. We used Mendelian randomization to examine whether uric acid is an independent and causal cardiovascular risk factor. Serum uric acid was measured in 3315 patients of the Ludwigshafen Risk and Cardiovascular Health Study. We calculated a weighted genetic risk score (GRS) for uric acid concentration based on eight uric acid-regulating single nucleotide polymorphisms. Causal odds ratios and causal hazard ratios (HRs) were calculated using a two-stage regression estimate with the GRS as the instrumental variable to examine associations with cardiometabolic phenotypes (cross-sectional) and mortality (prospectively) by logistic regression and Cox regression, respectively. Our GRS was not consistently associated with any biochemical marker except for uric acid, arguing against pleiotropy. Uric acid was associated with a range of prevalent diseases, including coronary artery disease. Uric acid and the GRS were both associated with cardiovascular death and sudden cardiac death. In a multivariate model adjusted for factors including medication, causal HRs corresponding to each 1-mg/dl increase in genetically predicted uric acid concentration were significant for cardiovascular death (HR, 1.77; 95% confidence interval, 1.12 to 2.81) and sudden cardiac death (HR, 2.41; 95% confidence interval, 1.16 to 5.00). These results suggest that high uric acid is causally related to adverse cardiovascular outcomes, especially sudden cardiac death.

Large-scale metabolomic profiling identifies novel biomarkers for incident coronary heart disease

Analyses of circulating metabolites in large prospective epidemiological studies could lead to improved prediction and better biological understanding of coronary heart disease (CHD). We performed a mass spectrometry-based non-targeted metabolomics study for association with incident CHD events in 1,028 individuals (131 events; 10 y. median follow-up) with validation in 1,670 individuals (282 events; 3.9 y. median follow-up). Four metabolites were replicated and independent of main cardiovascular risk factors [lysophosphatidylcholine 18∶1 (hazard ratio [HR] per standard deviation [SD] increment = 0.77, P-value<0.001), lysophosphatidylcholine 18∶2 (HR = 0.81, P-value<0.001), monoglyceride 18∶2 (MG 18∶2; HR = 1.18, P-value = 0.011) and sphingomyelin 28∶1 (HR = 0.85, P-value = 0.015)]. Together they contributed to moderate improvements in discrimination and re-classification in addition to traditional risk factors (C-statistic: 0.76 vs. 0.75; NRI: 9.2%). MG 18∶2 was associated with CHD independently of triglycerides. Lysophosphatidylcholines were negatively associated with body mass index, C-reactive protein and with less evidence of subclinical cardiovascular disease in additional 970 participants; a reverse pattern was observed for MG 18∶2. MG 18∶2 showed an enrichment (P-value = 0.002) of significant associations with CHD-associated SNPs (P-value = 1.2×10-7 for association with rs964184 in the ZNF259/APOA5 region) and a weak, but positive causal effect (odds ratio = 1.05 per SD increment in MG 18∶2, P-value = 0.05) on CHD, as suggested by Mendelian randomization analysis. In conclusion, we identified four lipid-related metabolites with evidence for clinical utility, as well as a causal role in CHD development.

Genome-wide association analysis identifies novel loci for chronotype in 100,420 individuals from the UK Biobank

Our sleep timing preference, or chronotype, is a manifestation of our internal biological clock. Variation in chronotype has been linked to sleep disorders, cognitive and physical performance, and chronic disease. Here we perform a genome-wide association study of self-reported chronotype within the UK Biobank cohort (n=100,420). We identify 12 new genetic loci that implicate known components of the circadian clock machinery and point to previously unstudied genetic variants and candidate genes that might modulate core circadian rhythms or light-sensing pathways. Pathway analyses highlight central nervous and ocular systems and fear-response-related processes. Genetic correlation analysis suggests chronotype shares underlying genetic pathways with schizophrenia, educational attainment and possibly BMI. Further, Mendelian randomization suggests that evening chronotype relates to higher educational attainment. These results not only expand our knowledge of the circadian system in humans but also expose the influence of circadian characteristics over human health and life-history variables such as educational attainment.

The shared allelic architecture of adiponectin levels and coronary artery disease.

OBJECTIVE: A large body of epidemiologic data strongly suggests an association between excess adiposity and coronary artery disease (CAD). Low adiponectin levels, a hormone secreted only from adipocytes, have been associated with an increased risk of CAD in observational studies. However, these associations cannot clarify whether this relationship is causal or due to a shared set of causal factors or even confounding. Genome-wide association studies have identified common variants that influence adiponectin levels, providing valuable tools to examine the genetic relationship between adiponectin and CAD. METHODS: Using 145 genome wide significant SNPs for adiponectin from the ADIPOGen consortium (n = 49,891), we tested whether adiponectin-decreasing alleles influenced risk of CAD in the CARDIoGRAM consortium (n = 85,274). RESULTS: In single-SNP analysis, 5 variants among 145 SNPs were associated with increased risk of CAD after correcting for multiple testing (P < 4.4 × 10(-4)). Using a multi-SNP genotypic risk score to test whether adiponectin levels and CAD have a shared genetic etiology, we found that adiponectin-decreasing alleles increased risk of CAD (P = 5.4 × 10(-7)). CONCLUSION: These findings demonstrate that adiponectin levels and CAD have a shared allelic architecture and provide rationale to undertake a Mendelian randomization studies to understand if this relationship is causal.

Genetic Loci Associated with C-reactive Protein Levels and Risk of Coronary Heart Disease.

CONTEXT: Plasma levels of C-reactive protein (CRP) are independently associated with risk of coronary heart disease, but whether CRP is causally associated with coronary heart disease or merely a marker of underlying atherosclerosis is uncertain. OBJECTIVE: To investigate association of genetic loci with CRP levels and risk of coronary heart disease. DESIGN, SETTING, AND PARTICIPANTS: We first carried out a genome-wide association (n = 17,967) and replication study (n = 13,615) to identify genetic loci associated with plasma CRP concentrations. Data collection took place between 1989 and 2008 and genotyping between 2003 and 2008. We carried out a mendelian randomization study of the most closely associated single-nucleotide polymorphism (SNP) in the CRP locus and published data on other CRP variants involving a total of 28,112 cases and 100,823 controls, to investigate the association of CRP variants with coronary heart disease. We compared our finding with that predicted from meta-analysis of observational studies of CRP levels and risk of coronary heart disease. For the other loci associated with CRP levels, we selected the most closely associated SNP for testing against coronary heart disease among 14,365 cases and 32,069 controls. MAIN OUTCOME MEASURE: Risk of coronary heart disease. RESULTS: Polymorphisms in 5 genetic loci were strongly associated with CRP levels (% difference per minor allele): SNP rs6700896 in LEPR (-14.8%; 95% confidence interval [CI], -17.6% to -12.0%; P = 6.2 x 10(-22)), rs4537545 in IL6R (-11.5%; 95% CI, -14.4% to -8.5%; P = 1.3 x 10(-12)), rs7553007 in the CRP locus (-20.7%; 95% CI, -23.4% to -17.9%; P = 1.3 x 10(-38)), rs1183910 in HNF1A (-13.8%; 95% CI, -16.6% to -10.9%; P = 1.9 x 10(-18)), and rs4420638 in APOE-CI-CII (-21.8%; 95% CI, -25.3% to -18.1%; P = 8.1 x 10(-26)). Association of SNP rs7553007 in the CRP locus with coronary heart disease gave an odds ratio (OR) of 0.98 (95% CI, 0.94 to 1.01) per 20% lower CRP level. Our mendelian randomization study of variants in the CRP locus showed no association with coronary heart disease: OR, 1.00; 95% CI, 0.97 to 1.02; per 20% lower CRP level, compared with OR, 0.94; 95% CI, 0.94 to 0.95; predicted from meta-analysis of the observational studies of CRP levels and coronary heart disease (z score, -3.45; P < .001). SNPs rs6700896 in LEPR (OR, 1.06; 95% CI, 1.02 to 1.09; per minor allele), rs4537545 in IL6R (OR, 0.94; 95% CI, 0.91 to 0.97), and rs4420638 in the APOE-CI-CII cluster (OR, 1.16; 95% CI, 1.12 to 1.21) were all associated with risk of coronary heart disease. CONCLUSION: The lack of concordance between the effect on coronary heart disease risk of CRP genotypes and CRP levels argues against a causal association of CRP with coronary heart disease.

Corrigendum : Novel loci affecting iron homeostasis and their effects in individuals at risk for hemochromatosis

Absorption, transport and storage of iron are tightly regulated, as expected for an element, which is both essential and potentially toxic. Iron deficiency is the leading cause of anaemia, and it also compromises immune function and cognitive development. Iron overload damages the liver and other organs in hereditary hemochromatosis, and in thalassaemia patients with both transfusion and non-transfusionrelated iron accumulation. Excess iron has harmful effects in chronic liver diseases caused by excessive alcohol, obesity or viruses. There is evidence for involvement of iron in neurodegenerative diseases and in Type 2 diabetes. Variation in transferrin saturation, a biomarker of iron status, has been associated with mortality in patients with diabetes and in the general population13. All these associations between iron and either clinical disease or pathological processes make it important to understand the causes of variation in iron status. Importantly, information on genetic causes of variation can be used in Mendelian randomization studies to test whether variation in iron status is a cause or consequence of disease. We have used biomarkers of iron status (serum iron, transferrin, transferrin saturation and ferritin), which are commonly used clinically and readily measurable in thousands of individuals, and carried out a meta-analysis of human genomewide association study (GWAS) data from 11 discovery and eight replication cohorts. Our aims were to identify additional loci affecting markers of iron status in the general population and to relate the significant loci to information on gene expression to identify relevant genes. We also made an initial assessment of whether any such loci affect iron status in HFE C282Y homozygotes, who are at genetic risk of HFE-related iron overload (hereditary hemochromatosis type 1, OMIM #235200)

Instrumental variable method for the causal relationship between soluble receptor for advanced glycation end-products (sRAGE) and risk of pancreatic cancer

This thesis project is motivated by the potential problem of using observational data to draw inferences about a causal relationship in observational epidemiology research when controlled randomization is not applicable. Instrumental variable (IV) method is one of the statistical tools to overcome this problem. Mendelian randomization study uses genetic variants as IVs in genetic association study. In this thesis, the IV method, as well as standard logistic and linear regression models, is used to investigate the causal association between risk of pancreatic cancer and the circulating levels of soluble receptor for advanced glycation end-products (sRAGE). Higher levels of serum sRAGE were found to be associated with a lower risk of pancreatic cancer in a previous observational study (255 cases and 485 controls). However, such a novel association may be biased by unknown confounding factors. In a case-control study, we aimed to use the IV approach to confirm or refute this observation in a subset of study subjects for whom the genotyping data were available (178 cases and 177 controls). Two-stage IV method using generalized method of moments-structural mean models (GMM-SMM) was conducted and the relative risk (RR) was calculated. In the first stage analysis, we found that the single nucleotide polymorphism (SNP) rs2070600 of the receptor for advanced glycation end-products (AGER) gene meets all three general assumptions for a genetic IV in examining the causal association between sRAGE and risk of pancreatic cancer. The variant allele of SNP rs2070600 of the AGER gene was associated with lower levels of sRAGE, and it was neither associated with risk of pancreatic cancer, nor with the confounding factors. It was a potential strong IV (F statistic = 29.2). However, in the second stage analysis, the GMM-SMM model failed to converge due to non- concaveness probably because of the small sample size. Therefore, the IV analysis could not support the causality of the association between serum sRAGE levels and risk of pancreatic cancer. Nevertheless, these analyses suggest that rs2070600 was a potentially good genetic IV for testing the causality between the risk of pancreatic cancer and sRAGE levels. A larger sample size is required to conduct a credible IV analysis.^

Caffeine intake and CYP1A2 variants associated with high caffeine intake protect non-smokers from hypertension.

The 15q24.1 locus, including CYP1A2, is associated with blood pressure (BP). The CYP1A2 rs762551 C allele is associated with lower CYP1A2 enzyme activity. CYP1A2 metabolizes caffeine and is induced by smoking. The association of caffeine consumption with hypertension remains controversial. We explored the effects of CYP1A2 variants and CYP1A2 enzyme activity on BP, focusing on caffeine as the potential mediator of CYP1A2 effects. Four observational (n = 16 719) and one quasi-experimental studies (n = 106) including European adults were conducted. Outcome measures were BP, caffeine intake, CYP1A2 activity and polymorphisms rs762551, rs1133323 and rs1378942. CYP1A2 variants were associated with hypertension in non-smokers, but not in smokers (CYP1A2-smoking interaction P = 0.01). Odds ratios (95% CIs) for hypertension for rs762551 CC, CA and AA genotypes were 1 (reference), 0.78 (0.59-1.02) and 0.66 (0.50-0.86), respectively, P = 0.004. Results were similar for the other variants. Higher CYP1A2 activity was linearly associated with lower BP after quitting smoking (P = 0.049 and P = 0.02 for systolic and diastolic BP, respectively), but not while smoking. In non-smokers, the CYP1A2 variants were associated with higher reported caffeine intake, which in turn was associated with lower odds of hypertension and lower BP (P = 0.01). In Mendelian randomization analyses using rs1133323 as instrument, each cup of caffeinated beverage was negatively associated with systolic BP [-9.57 (-16.22, -2.91) mmHg]. The associations of CYP1A2 variants with BP were modified by reported caffeine intake. These observational and quasi-experimental results strongly support a causal role of CYP1A2 in BP control via caffeine intake.

Stratification by smoking status reveals an association of CHRNA5-A3-B4 genotype with body mass index in never smokers.

We previously used a single nucleotide polymorphism (SNP) in the CHRNA5-A3-B4 gene cluster associated with heaviness of smoking within smokers to confirm the causal effect of smoking in reducing body mass index (BMI) in a Mendelian randomisation analysis. While seeking to extend these findings in a larger sample we found that this SNP is associated with 0.74% lower body mass index (BMI) per minor allele in current smokers (95% CI -0.97 to -0.51, P = 2.00 × 10(-10)), but also unexpectedly found that it was associated with 0.35% higher BMI in never smokers (95% CI +0.18 to +0.52, P = 6.38 × 10(-5)). An interaction test confirmed that these estimates differed from each other (P = 4.95 × 10(-13)). This difference in effects suggests the variant influences BMI both via pathways unrelated to smoking, and via the weight-reducing effects of smoking. It would therefore be essentially undetectable in an unstratified genome-wide association study of BMI, given the opposite association with BMI in never and current smokers. This demonstrates that novel associations may be obscured by hidden population sub-structure. Stratification on well-characterized environmental factors known to impact on health outcomes may therefore reveal novel genetic associations.

Association between C-reactive protein and adiposity in women.

Context: The link between C-reactive protein (CRP) and adiposity deserves to be further explored considering the controversial diabetogenic role of CRP. Objective: We explored the potential causal role of CRP on measures of adiposity. Design: We used a Mendelian randomization approach with the CRP and LEPR genes as instrumental variables in a cross-sectional Caucasian population-based study comprising 2526 men and 2836 women. Adiposity was measured using body mass index (BMI), fat and lean mass estimated by bioelectrical impedance, and waist circumference. Results: Log-transformed CRP explained by the rs7553007 SNP tagging the CRP gene was significantly associated with BMI (regression coefficient: 1.22 [0.18;2.25], P=0.02) and fat mass (2.67 [0.65;4.68], P=0.01), but not with lean mass in women, whereas no association was found in men. Log-transformed CRP explained by the rs1805096 LEPR SNP was also positively associated, although not significantly, with BMI or fat mass. The combined CRP-LEPR instrument explained 2.24% and 0.77% of CRP variance in women and in men, respectively. Log-transformed CRP explained by this combined instrument was significantly associated with BMI (0.98 [0.32;1.63], P=0.004), fat mass (2.07 [0.79;3.34], P=0.001) and waist (2.09 [0.39;3.78], P=0.01) in women, but not in men. Conclusion: Our data suggest that CRP is causally and positively related to BMI in women, and that this is mainly due to fat mass. Results on the combined CRP-LEPR instrument suggest that leptin may play a role in the causal association between CRP and adiposity in women. Results in men were not significant.

Uric acid and cardio-metabolic risk factors : an epidemiological approach

The role of serum uric acid (SUA) in cardio-metabolic conditions has long been contentious. It is still unclear if SUA is an independent risk factor or marker of cardio-metabolic conditions and most observed associations are not necessarily causal. This study aimed to further understand and explore the causal role of SUA in cardio-metabolic conditions using genetic and non-genetic epidemiological methods in population-based data. In the first part of this study, we found moderate to high heritability estimates for SUA and fractional excretion of urate (FEUA) suggesting the role of genetic factors in the etiology of hyperuricemia. With regards to the role of SUA on inflammatory markers (IMs), a strong positive association of SUA with C-reactive protein (CRP) and a weaker positive association with tumor necrosis factor alpha (TNF-α) and interleukin 6 (IL-6) was observed, which was in part mediated by body mass index (BMI). These findings suggest that SUA may have a role in sterile inflammation. In view of the inconsistency surrounding the causal nature and direction of the relation between SUA and adiposity, we applied a bidirectional Mendelian randomization approach using genetic variants to decipher the association. The finding that elevated SUA is a consequence rather than a cause of adiposity was not totally unexpected and is compatible with the hypothesis that hyperinsulinemia, accompanying obesity, enhances renal proximal tubular reabsorption of uric acid. The fourth part of this study examined the relationship between SUA and blood pressure (BP) in young adults. The association between SUA and BP, significant only in females, was strongly attenuated upon adjustment for BMI. The possibility that BMI lies in the causal pathway may explain the attenuation observed in the associations of SUA with BP and IMs. Finally, a significant hockey-stick shaped association of SUA with social phobia in our data suggests a protective effect of SUA only up to a certain concentration. Although our study findings have shed some light on the uncertainty underlying the pathophysiology of SUA, more compelling evidence using longitudinal designs, randomized controlled trials and the use of robust genetic tools is warranted to increase our understanding of the clinical significance of SUA.

Role of vitamin D deficiency in extraskeletal complications : predictor of health outcome or marker of health status?

The relationship of vitamin D with extraskeletal complications, such as cardiovascular disease, cancer, and autoimmune disease, is of major interest considering its roles in key biological processes and the worldwide high prevalence of vitamin D deficiency. However, the causal relationships between vitamin D and most extraskeletal complications are weak. Currently, a heated debate over vitamin D is being conducted according to two hypotheses. In this review, we first present the different arguments that suggest a major role of vitamin D in a very broad type of extraskeletal complications (hypothesis #1). We then present results from recent meta-analyses of randomized controlled trials indicating a lack of association of vitamin D with major extraskeletal complications (hypothesis #2). We discuss different issues (e.g., causality, confounding, reverse causation, misclassification, and Mendelian randomization) that contribute to the favoring of one hypothesis over the other. While ultimately only one hypothesis is correct, we anticipate that the results from the ongoing randomized controlled trials will be unlikely to reconcile the divided experts.