Genetic variation in CADM2 as a link between psychological traits and obesity
0000 0001 2193 314Xgrid.8756.cInstitute of Health and Wellbeing, University of Glasgow, Glasgow, UK
0000 0001 2193 314Xgrid.8756.cSchool of Life Sciences, College of Medical, Veterinary and Life Sciences, University of Glasgow, Glasgow, UK
0000 0004 1757 2822grid.4708.bDepartment of Medical Biotechnology and Translational Medicine, Università degli Studi di Milano, Milan, Italy
0000 0004 1760 1750grid.418230.cCentro Cardiologico Monzino, IRCCS, Milan, Italy
0000 0004 1937 0626grid.4714.6Unit of Cardiovascular Epidemiology, Institute of Environmental Medicine, Karolinska Institutet, Stockholm, Sweden
0000 0004 0636 5158grid.412154.7Division of Cardiovascular Medicine, Department of Clinical Sciences, Danderyd University Hospital, Stockholm, Sweden
Assistance Publique-Hopitaux de Paris, Service Endocrinologie-Metabolisme, Groupe Hôpitalier Pitie-Salpetriere, Unités de Prévention Cardiovasculaire, Paris, France
0000 0004 1936 8948grid.4991.5Division of Cardiovascular Medicine, Radcliffe Department of Medicine, University of Oxford, Oxford, UK
0000 0004 1936 8948grid.4991.5Wellcome Centre for Human Genetics, University of Oxford, Oxford, UK
0000 0004 1937 0626grid.4714.6Cardiovascular Medicine Unit, Department of Medicine Solna, Karolinska Institutet, Stockholm, Sweden
0000000121901201grid.83440.3bCentre for Cardiovascular Genetics, Institute Cardiovascular Science, University College London, London, UK
0000 0004 1936 7988grid.4305.2Division of Psychiatry, College of Medicine, University of Edinburgh, Edinburgh, UK
0000 0004 1936 9457grid.8993.bDepartment of Cell and Molecular Biology, National Bioinformatics Infrastructure Sweden, Science for Life Laboratory, Uppsala University, Uppsala, Sweden
0000 0000 9558 4598grid.4494.dDepartment of Medicine, University Medical Center Groningen and University of Groningen, Groningen, The Netherlands
0000 0004 1757 2822grid.4708.bDipartimento di Scienze Farmacologiche e Biomolecolari, Università di Milano, Milan, Italy
Abstract
CADM2 has been associated with a range of behavioural and metabolic traits, including physical activity, risk-taking, educational attainment, alcohol and cannabis use and obesity. Here, we set out to determine whether CADM2 contributes to mechanisms shared between mental and physical health disorders. We assessed genetic variants in the CADM2 locus for association with phenotypes in the UK Biobank, IMPROVE, PROCARDIS and SCARFSHEEP studies, before performing meta-analyses. A wide range of metabolic phenotypes were meta-analysed. Psychological phenotypes analysed in UK Biobank only were major depressive disorder, generalised anxiety disorder, bipolar disorder, neuroticism, mood instability and risk-taking behaviour. In UK Biobank, four, 88 and 172 genetic variants were significantly (p < 1 × 10−5) associated with neuroticism, mood instability and risk-taking respectively. In meta-analyses of 4 cohorts, we identified 362, 63 and 11 genetic variants significantly (p < 1 × 10−5) associated with BMI, SBP and CRP respectively. Genetic effects on BMI, CRP and risk-taking were all positively correlated, and were consistently inversely correlated with genetic effects on SBP, mood instability and neuroticism. Conditional analyses suggested an overlap in the signals for physical and psychological traits. Many significant variants had genotype-specific effects on CADM2 expression levels in adult brain and adipose tissues. CADM2 variants influence a wide range of both psychological and metabolic traits, suggesting common biological mechanisms across phenotypes via regulation of CADM2 expression levels in adipose tissue. Functional studies of CADM2 are required to fully understand mechanisms connecting mental and physical health conditions.
Introduction
The link between mental illness and poor physical health is well-established but reasons for this are poorly understood. Patients with severe psychiatric illness have 2–3 times higher rate of metabolic syndrome than the general population worldwide, likely due to a combination of factors including the effects of psychotropic drugs, negative health behaviours, hormone dysregulation and shared genetic risk factors1,2. A number of potential shared pathways between mood disorders and cardiometabolic disease have been suggested, including abnormal circadian rhythms, hypothalamic-pituitary-adrenal (HPA) axis dysfunction and inflammation. However, the molecular mechanisms of these pathways are poorly understood.
Single nucleotide polymorphisms (SNPs) in the locus encoding the synaptic cell adhesion molecule 2 (CADM2) on chromosome 3 have been associated with a number of psychological traits, including educational attainment3, alcohol consumption4, cannabis use5, physical activity habits6, risk-taking behaviour7,8, attention-deficit/hyperactivity disorder9 and obesity10. Several lines of evidence point to CADM2 being the gene through which SNPs are having their effects, including genotype-specific effects on CADM2 mRNA expression levels7,8, CADM2 being predominantly expressed in the brain, and cadm2 knockout models demonstrating relevant phenotypes. Specifically, cadm2-knockout mice have reduced adiposity, reduced systemic glucose levels, improved insulin sensitivity, increased locomotor activity, increased energy expenditure rate and raised core body temperature, suggesting an important role in systemic energy homeostasis11.
We set out to systematically evaluate the relationship between CADM2 SNPs and psychological and physical traits, and assess whether there is evidence for distinct signals influencing metabolic versus psychological traits.
Materials and Methods
CADM2 locus
We defined the CADM2 locus as the CADM2 gene plus 250 kb up and downstream (Chromosome 3:84758000-86374000, UCSC genome browser, https://genome-euro.ucsc.edu/).
Study cohorts
High CVD risk population: IMPROVE is a cohort of individuals with no symptoms or history of cardiovascular disease, but with a least three classic risk factors (namely any combination of the following: family history of CVD, type 2 diabetes, smoking, hypertension, dyslipidaemia, male sex or women at least 5 years post-menopause)12. In brief, 3,711 participants were recruited from 7 centres across 5 European countries (Finland, Sweden, the Netherlands, France and Italy) between January 2004 and June 2005. Participants completed a structured medical history and lifestyle questionnaire at baseline, as well as standard biochemical tests and genotyping. Ethics committee approval was granted by the Institutional review board (IRB) at each recruitment centre: Karolinska Institutet, Stockholm, Sweden; University of Milan, Milan, Italy; University of Kuopio and Kuopio Research Institute of Exercise Medicine, Kupio, Finland; University Hospital Groningen, Groningen, The Netherlands; University of Perugia, Perugia, Italy; Groupe Hôpital Pitie-Salpetriere, Paris, France. Informed consent was provided by all participants. The study was conducted in accordance with the Helsinki Declaration.
Young CVD case-control cohort: SCARFSHEEP is a case-control cohort of Swedish participants (N = 2,513) recruited in Stockholm13,14. Cases were those with a first myocardial infarction before 60 years of age. Controls were age and sex-matched from the general population of the same county. Standard biochemical phenotyping was available for all participants. Approval was granted by the Karolinska Hospital and Karolinska Institutet Ethics Committees (for SCARF and SHEEP respectively). Informed consent was provided by all participants. The study conducted in accordance with the Helsinki Declaration.
CVD case-control cohort: PROCARDIS is a case-control cohort15, where cases (n = 5,688) were diagnosed with coronary artery disease before 66 years and controls (n = 2,310) are unrelated participants without coronary artery disease at 66 years. Participants were recruited from 4 centres across 4 European countries (Sweden, the UK, Germany and Italy). Participants completed a questionnaire at baseline. Standard biochemical phenotyping was available for all participants. Ethics Committee approval was granted by the IRB at each recruitment centre: the Regional Ethics Review Board at Karolinska Institutet, Stockholm in Sweden, the IRB at the University of Munster, Munster, in Germany, the IRB at the Mario Negri Institute, Milano in Italy and the IRB at the University of Oxford, Oxford, United Kingdom. Informed consent was provided by all participants. The study conducted in accordance with the Helsinki Declaration.
General population cohort: UK Biobank is a cohort of over 500,000 participants aged 40–69 at baseline16. Participants were recruited from 22 centres across the UK between 2006 and 2010. Participants completed a wide variety of baseline questionnaires and assessments. For consistency with the other cohorts, only white British participants were included. All participants provided informed consent. This study was carried out under the generic approval from the NHS National Research Ethics Service (approval letter dated 13 May 2016, Ref 16/NW/0274) and under UK Biobank applications #6553 (PI Daniel Smith) and #17689 (PI Donald Lyall). The study was conducted in accordance with the Helsinki Declaration.
Genetic data
IMPROVE participants were genotyped using both Illumina Cardio-Metabo17 and Immuno18 arrays at the SNP&SEQ Technology Platform in Uppsala. SCARFSHEEP were genotyped using the Illumina Cardio-Metabo chip17 at the SNP&SEQ Technology Platform in Uppsala. For both cohorts, imputation to the 1000 Genomes reference panel was conducted according to standard protocols, as described previously19.
PROCARDIS participants were genotyped at the Centre National du Genotypage, Paris and the SNP&SEQ Technology Platform in Uppsala, using the Illumina 1 M and 610 K arrays. Imputation to the 1000 Genomes panel was conducted according to standard protocols, as described20.
UK Biobank participants were genotyped using either the Affymetrix UK Biobank Axion or the Affymetrix BiLEVE Axion array16. A modified version of SHAPEIT2 was used for phasing and IMPUTE2 for imputation. The data from UK Biobank was released in two phases. The UK Biobank was imputed to the 1000 Genomes, UK10K haplotype (first release) and Haplotype Reference Consortium (merged with the first release for the second release) reference panels21.
We applied standard quality control procedures to all cohorts, including SNP exclusion for low call rate (<95%), minor allele frequency (MAF < 1%), deviation from Hardy-Weinberg equilibrium (p < 5 × 10−6) or imputation quality score <0·4 and subject exclusion for sex mismatch, cryptic relatedness, low call rate (<95%) and non-Caucasian ancestry (self-reported or based on principle component analysis). For UK Biobank exclusions, further exclusions based on relatedness were applied (one of each pair of individuals with a KING-estimate kindship coefficient >0·0442 was randomly removed). After quality control, 5,684, 2,786, 5,452 and 5,361 SNPs were available for IMPROVE, SCARFSHEEP, PROCARDIS and UK Biobank respectively. In total, 2,123 SNPs were available in all four cohorts, with 2,133 overlapping between the three CAD case-control cohorts and 2,434 overlapping in the three cohorts with biomarker data.
Phenotypes
Psychiatric and psychological phenotypes were only available in the UK Biobank. The baseline questionnaire included questions to assess mood instability (“does your mood often go up and down?” variable #1920) and risk-taking behaviour (“Would you describe yourself as someone who takes risks” variable #2040). Single item questions are imperfect ways to measure psychological traits, however validity of the question used here has been demonstrated relative to more detailed phenotyping (at least for risk-taking)22 and in terms of the expected associations with psychiatric disorders (for mood instability23 and risk-taking7,8,24,25). Neuroticism was assessed using the Eysenck Personality Questionnaire (Revised Short Form), where 12 yes/no questions were asked. These were summed, resulting in a score between one and 12 for each individual26. Phenotyping in relation to psychiatric disorders was based upon the online “Thoughts and Feelings” questionnaire27, which requested information on lifetime symptoms of mental disorders. This enabled classification of likely major depressive disorder (MDD), bipolar disorder (BD), generalised anxiety disorder (GAD) and addiction).
Anthropometric and blood pressure phenotypes (BMI, waist and hip circumferences, SBP and DBP) were assessed in a standardised and comparable manner. Waist to hip circumference ratio adjusted for BMI (WHRadjBMI) was calculated as per Shungin et al.28. For those on anti-hypertensive medication, values of SBP and DBP were adjusted, with 15 and 10 mmHG, respectively, being added prior to analysis29. Current smoking was assessed by questionnaire in all cohorts.
Metabolic parameters were available in IMPROVE, SCARFSHEEP and PROCARDIS, where fasting glucose, lipid (HDL, LDL and TG) and CRP levels were measured using standard methodology at the Department of Clinical Chemistry, Karolinska University Hospital. Fasting insulin levels were measured by radio-immunoassay30,31. HOMA indices were calculated from fasting glucose and insulin levels as described32. Type 2 diabetes was defined as diagnosis, medication and/or fasting glucose levels ≥7 mmol/L for IMPROVE, SCARFSHEEP and PROCARDIS. The definition of T2D in UK Biobank has been described33, and is generally comparable with the assessment used for the other cohorts.
Coronary vascular disease (CVD) was defined as clinically diagnosed myocardial infarction, symptomatic acute coronary syndrome, angina or coronary artery revascularisation before the age of 66 years for PROCARDIS15. Criteria for inclusion as a CVD case in the SCARFSHEEP study was clinical diagnosis of myocardial infarction diagnosis34. For UK Biobank, CVD was defined as clinical diagnosis of heart attack/myocardial infarction or angina (variable # 6150).
Statistical analyses
All continuous phenotypes were assessed for normality and, where necessary, were natural log transformed prior to analysis. For each cohort, phenotypes were analysed in PLINK 1.07 using linear or logistic regression (for continuous vs. binary traits respectively), assuming additive allelic effects. With the exception of WHRadjBMI, all models included age, sex and population structure (3 principal components for PROCARDIS, SCARSHEEP and IMPROVE, 8 for UK Biobank), with further adjustment for genotyping chip being applied for UK Biobank and PROCARDIS analyses. For analysis of lipid traits, lipid-lowering medication and CVD case-control status were included as a covariates. For glucometabolic traits, individuals with type 2 diabetes were excluded and CVD case-control status was included as a covariate.
Results from the individual studies were combined in inverse variance-weighted meta-analyses using METAL35 (with binary effect sizes being analysed as Beta coefficients). Inverse variance-weighted meta-analysis was chosen, as the phenotype measurements and data transformation were comparable (including consistent units) between studies. Averages and standard errors of allele frequencies were computed. No additional filters were applied. Supplementary Table 1 summarises the phenotypes analysed in each cohort, covariates used and total sample number in the meta-analyses. Only SNPs present in 3 (of 3 or 4) cohorts were considered. Despite the prior knowledge implicating this locus in mental and physical health traits, we used a conservative approach, with genome-wide significance being set at p < 5 × 10−8 and suggestive evidence of association being set at p < 1 × 10−5. Locuszoom was used to visualise the results36.
Genetic architecture
In order to determine whether the different traits have distinct signals, or whether there is one signal influencing all traits, two approaches were used:
Firstly, SNPs meeting suggestive or genome-wide significance thresholds for at least one phenotype (candidate SNPs) were identified and linkage disequilibrium (LD) assessed. For analysis of genetic architecture, a random subset of 1000 unrelated white British participants from the UK Biobank were selected. In this subset of UK Biobank, candidate SNPs were filtered to leave only independent SNPs, using PLINK (independent pairwise selection with default settings, including LD r2 threshold 0·5). LD between the independent SNPs and lead/index SNPs was calculated and visualised using Haploview37.
Secondly, conditional analyses were performed to further examine the possibility of multiple signals in the CADM2 locus. Here, the risk-taking and BMI analyses were repeated, with the index SNP (coded as an additive genetic effect, namely 0 for common homozygote, 1 for heterozygotes and 2 for rare homozygotes) from each other phenotype in turn included as a covariate.
Data-mining
The GWAS catalogue (https://www.ebi.ac.uk/gwas/, accessed 2018–09–04, 3:84,758,000–86,374,001) was used to identify CADM2 locus SNPs previously associated with relevant phenotypes (specifically cardio-metabolic and psychiatric disorder-related traits). All SNPs in the CADM2 locus with suggestive or genome-wide evidence for association with at least one phenotype were assessed for predicted functional effects using the Variant Effect Predictor38. For lead and index SNPS, the GTEx portal39 was queried to identify genotype-specific gene expression patterns (or expression quantitative traits loci (eQTLs)).
Results
The cohort characteristics are presented in Table 1 and the phenotypes assessed are presented in Table 2.UK Biobank IMPROVE PROCARDIS SCARFSHEEP N 408961 3390 7998 3417 Male 221052 (54.0) 1634 (48.2) 6009 (75.1) 2459 (72.0) CAD cases 8319 (2.0) 0 (0) 5688 (71.1) 1525 (44.6) T2D cases 17766 (4.3) 908 (26.8) 1157 (14.5) 341 (10.0) Age (years) 56.9 (8.0) 64.2 (5.4) 60.0 (8.3) 58.1 (7.27) BMI (kg/m2) 27.42 (4.76) 27.2 (4.24) 27.9 (4.4) 26.3 (3.9) Waist:hip ratio 0.87 (0.09) 0.91 (0.09) 0.97 (0.08) 0.94 (0.09) SBP (mmHg) 137 (19) 142 (18) 135 (20) 135 (21) DBP (mmHG) 82 (10) 82 (10) 81 (11) 82 (10) Current smoking 144005 (35.3) 507 (15.0) 1137 (15.4) 1133 (33.8) Lipid-lowering medication 43424 (26.0) 1676 (49.5) 4175 (56.6) 235 (6.8) Anti-hypertensive medication 85463 (21.0) 1950 (57.5) 4592 (62.3) 1644 (48.1) LDL cholesterol (mmol/L) na 3.54 (1.00) 2.98 (0.87) 4.02 (0.99) HDL cholesterol (mmol/L) na 1.26 (0.36) 1.22 (0.36) 1.17 (0.36) Triglycerides (mmol/L) na 5.93 (1.66) 1.83 (1.26) 1.70 (1.16) CRP (mmol/L) na 2.97 (5.76) 3.44 (7.20) 2.94 (5.69) Fasting glucose (mmol/L)* na 5.28 (0.66) 5.38 (0.60) 5.10 (0.67) Fasting insulin (pmol/L)* na 44.6 (61.8) 57.9 (54.3) 58.6 (45.2) HOMA B* na 69.1 (54.1) 79.8 (32.6) 25.7 (11.0) HOMA IR* na 0.83 (1.10) 1.02 (0.72) 0.16 (0.11) BD 1899 (1.4) na na na GAD 9251 (10.2) na na na MDD 31338 (28.2) na na na Addiction 7575 (5.8) na na na Mood instability 180743 (44.2) na na na Risk-takers 32735 (25.5) na na na Neuroticism score 4.11 (3.26) na na na Phenotype UK Biobank Procardis IMPROVE SCARFSHEEP Maximum N BMI* Y Y Y Y 416136 WHRadjBMI* Y Y Y Y 416136 SBP* Y Y Y Y 397163 DBP* Y Y Y Y 397163 Diabetes* Y Y Y Y 417615 Current smoking* Y Y Y Y 416136 CVD Y Y Y 414116 TGs* Y Y Y 14799 CRP* Y Y Y 14799 HDL* Y Y Y 14799 LDL* Y Y Y 14799 Glucose*^ Y Y Y 10128 Insulin*^ Y Y Y 10128 HOMA-IR*^ Y Y Y 10128 HOMA-B*^ Y Y Y 10128 BD Y 129366 GAD Y 90536 MDD Y 109436 Addiction Y 129858 Mood instability Y 393367 Risk taking Y 328339 Neuroticism Y 328087
Meta-analysis of cardiovascular and metabolic phenotypes
Only SNPs present in at least 3 of the 4 cohorts were considered. It is worth noting that the heterogeneity I2 value was high for many SNPs in the meta-analysis of UK Biobank, IMPROVE, PROCARDIS and SCARFHSEEP. This could be due to selection of UK vs European or population vs case-control participants. Therefore we present results for both a lead SNP (defined as the SNP with the lowest P-value) and an index SNP (defined as the SNP with the lowest p-value with heterogeneity I2 = 0), but for robustness we focus on the Index SNPs.
Evidence of association was observed for BMI, with 908 SNPs reaching genome-wide significance (Index SNP rs11708632-A, Beta 0·085, p = 2·18 × 10−14) (Fig. 1A and Table 3), of which 166 had heterogeneity I2 = 0. A further 362 SNPs met the threshold for suggestive evidence of association with BMI (p < 1 × 10−5), 140 of which had heterogeneity I2 = 0 (Supplementary Table 1). For SBP, only one SNP reached genome-wide significance but heterogeneity was high (heterogeneity I2 = 69%, rs146071762-A, Beta 0·248, p = 1·14 × 10−8). An additional 63 SNPs demonstrated suggestive associations with SBP, with 10 showing heterogeneity I2 = 0 (Index SNP, rs6803322-A, Beta −0·024, p = 1·41 × 10−7) (Fig. 1B, Table 3, Supplementary Table 1).Phenotype Chr SNP Position A1 A1F Fse Beta Se P Effect I2 Phet Description N GWAS SNPs N suggestive SNPs BMI 3 rs11915747 85699040 C 0.64 0.01 0.090 0.010
1.58E-19
++++ 73 0.011 lead SNP 908 362 BMI 3 rs11708632 85766667 A 0.24 0.00 0.085 0.011
2.18E-14
++++ 0 0.417 Index SNP with I2 = 0 SBP 3 rs146071762 85396778 A 0.51 0.02 0.248 0.044
1.14E-08
+−+? 69 0.039 lead SNP 1 63 SBP 3 rs6803322 84986088 A 0.32 0.00 −0.243 0.046 1.41E-07 −+−? 0 0.619 Index SNP with I2 = 0 CRP 3 3:85906663** 85906663 T 0.31 0.01 −0.072 0.015 2.69E-06 −− 16 0.303 lead SNP 11 CRP 3 rs11708024 85865269 A 0.30 0.01 −0.069 0.015 6.05E-06 −− 0 0.674 Index SNP with I2 = 0 risk-taking* 3 rs4856591 85612550 T 0.38 0.056 0.006
1.02E-19
Index SNP 809 172 mood instability* 3 rs818225 85382140 T 0.46 −0.023 0.005 8.12E-07 Index SNP 4 neuroticism* 3 rs818219 85374589 C 0.46 −0.038 0.008 1.54E-06 Index SNP 88
No associations were observed for WHRadjBMI, DBP, T2D or current smoking (Supplementary Fig. 1A–D). For CVD analysis, SNPs were only considered if they were present in all three cohorts. None met the threshold for suggestive significance (Supplementary Fig. 1E).
CADM2 vs psychological and psychiatric phenotypes in the UK Biobank
As previously reported7,8, CADM2 variants were associated with risk-taking behaviour, with 809 SNPs demonstrating genome-wide evidence of association (p < 5 × 10−8 for risk-taking behaviour (Index SNP rs4856591-T, Beta = 0.056, p = 1.02 × 10−19); Fig. 1D and Table 3). In addition, suggestive evidence of association was observed for neuroticism (88 SNPs, index SNP rs818219-C, beta = −0·038, p = 1.54 × 10−6) and mood instability (4 SNPs, index SNP rs818225-T, beta = −0.02283, p = 8.12 × 10−7; Fig. 1E,F and Table 3). No SNPs demonstrated suggestive evidence of association with MDD, GAD, BD or addiction (Supplementary Fig. 3).
Cross-trait observations
A total of 49 SNPs demonstrated at least suggestive associations with multiple phenotypes (Supplementary Table 2). This observation demonstrates, firstly that the same SNPs influence both metabolic and psychological traits, and secondly that effects on risk-taking, BMI and CRP were positively correlated and these were inversely correlated with effects on neuroticism, mood instability and SBP.
Genetic architecture of CADM2
In order to determine whether the associations with psychological and physical traits reflect the same or distinct signals, two approaches were used. Firstly, filtering the 1,533 candidate SNPs (SNPs meeting suggestive significance of association with any phenotype) in the 1.62 Mb CADM2 locus by LD gave 75 independent loci (Fig. 2A). The index SNPs for neuroticism (rs818219) and mood instability (rs818225) are in perfect LD, therefore represent the same signal, whereas LD between index SNPs for other traits (BMI, rs11708632; SBP, rs6803322; CRP, rs11708024; risk-taking, rs485659) is low (maximum r2 = 0·37, Fig. 2B), which could indicate independent signals for each other phenotype.
Secondly, conditional analysis using the index SNPs (BMI, rs11708632; SBP, rs6803322; CRP, rs11708024; risk-taking, rs4856591; neuroticism, rs818219; mood instability, rs818225) was performed, using the UK Biobank data. If there is only one signal in the locus, then adjusting for the index SNP would remove the effects/significance of other SNPs in the locus. Alternatively, if there are additional signals in the locus which are independent of the index SNP, then adjusting for the index SNP would have little or no impact on the independent signal in terms of effect size or significance. The risk-taking results demonstrated that inclusion of index SNPs had some effect on the p-value of the risk-taking signal (Fig. 3), but the effect size was stable (primary risks analysis Beta = 0.056, conditional Betas = 0.054–0.057), which concurs with the LD analysis suggestion that the signals are independent. The BMI results were similar, but the effect size was less stable (primary BMI analysis Beta = −0.099, conditional Betas = −0.087–0.101) and the plots support the possibility of more than one signal in this region (Supplementary Fig. 4).
eQTL analysis
Firstly, SNPs with genotype-specific effects on mRNA CADM2 levels were identified using the GTEx portal (Fig. 4). The location of the eQTLs appears to be tissue specific (Fig. 4): In subcutaneous adipose tissue, eQTLs are located upstream and centrally relative to the gene location, whereas those in visceral adipose tissue are restricted to the central part of the gene (Fig. 4A,B). In heart (left ventricle) tissue and skeletal muscle, eQTLs are preferentially (but not exclusively) located centrally, whereas eQTLs in lung tissue are mainly upstream (Fig. 4D,F). In contrast, there are no eQTLs for CADM2 in the brain (cerebellum), instead the eQTLs in this brain region are for CADM2-AS1 (Fig. 4C). These findings suggest the potential for differential regulation of CADM2 levels across a range of tissues, thereby explaining how different SNPs within the locus can influence a variety of traits.
Secondly, 1747 SNPs with eQTL effects on CADM2 and CADM2-AS1 were identified and are presented in Supplementary Table 3. No eQTLs for CADM2-AS2 were identified. For CADM2 there were 3702 eQTLS (each SNP can have eQTL effects in more than one tissue), of which 41% were in adipose tissue (subcutaneous or visceral) and 5% in brain tissues. For CADM2-AS1, 1628 eQTLs were identified, of which 90% were in brain tissues and none were in adipose tissues.
Finally, all of the index SNPs available in GTEx (rs1170802 was not available) demonstrated eQTLs for CADM2 and none were eQTLs for CADM2-AS1 (Table 4). Expression levels of CADM2 are highest in the brain (Fig. 5A). However, it was interesting to note that the trait-increasing alleles of rs11708632 (BMI) and rs4856591 (risk-taking) were associated with increased CADM2 expression levels in adipose tissue (Fig. 5B,C). Consistent with their inverse correlations with BMI and risk-taking, the trait-decreasing alleles of rs818225 (mood instability) and rs818219 (neuroticism) were associated with increased CADM2 expression in adipose tissue (Fig. 5D,E).SNP Gene Symbol SNP Id P-Value EA NES Tissue Trait-increasing allele BMI lead
CADM2
rs11915747 8.60E-07 G −0.23 Lung C BMI lead
CADM2-AS1
rs11915747 1.30E-04 G 0.25 Nerve - Tibial C BMI index
CADM2
rs11708632 6.40E-06 A 0.27 Lung A BMI index
CADM2
rs11708632 5.10E-05 A 0.23 Adipose - Visceral (Omentum) A SBP index
CADM2
rs6803322 1.10E-08 A 0.71 Brain - Spinal cord (cervical c-1) C SBP index
CADM2
rs6803322 1.50E-06 A 0.19 Muscle - Skeletal C SBP index
CADM2
rs6803322 7.20E-06 A 0.2 Adipose - Subcutaneous C CRP lead
CADM2-AS1
3:85906663 2.80E-08 T 0.52 Brain - Cerebellum TGTTGCTCAG CRP lead
CADM2-AS1
3:85906663 1.30E-06 T 0.56 Brain - Caudate (basal ganglia) TGTTGCTCAG CRP lead
CADM2-AS1
3:85906663 1.60E-05 T 0.29 Nerve - Tibial TGTTGCTCAG CRP lead
CADM2-AS1
3:85906663 3.60E-05 T 0.49 Brain - Cerebellar Hemisphere TGTTGCTCAG risk-taking
CADM2
rs4856591 6.30E-16 G −0.39 Lung T risk-taking
CADM2
rs4856591 6.30E-07 G −0.22 Adipose - Subcutaneous T risk-taking
CADM2
rs4856591 1.70E-06 G −0.39 Heart - Left Ventricle T risk-taking
CADM2
rs4856591 2.10E-06 G −0.22 Adipose - Visceral (Omentum) T neuroticism
CADM2
rs818219 3.80E-06 C 0.19 Adipose - Subcutaneous T neuroticism
CADM2
rs818219 1.20E-04 C 0.14 Muscle - Skeletal T mood instability
CADM2
rs818225 4.20E-05 T 0.17 Adipose - Subcutaneous C
Data-mining
None of the index SNPs have previously been associated with any traits in the GWAS catalogue. Additionally, none of the candidate SNPs were predicted by Variant Effect Predictor to have more than low or modifier impact.
SNPs in the CADM2 locus which have previously been associated with psychological or metabolic traits are presented in Supplementary Table 4. Where comparison was possible, SNPs with reported effects on measures of obesity demonstrated consistent effect directions in our study compared to those previously reported10,28,40,41, with one exception: rs12495178 has the opposite effect directions in the Japanese population42 compared to our study. It is noteworthy that SNPs associated with educational attainment or intelligence3,43, were associated with increased BMI, which is somewhat surprising, but consistent with effects on risk-taking behaviour7,8. Maybe unsurprisingly, rs62253088-T (associated with strenuous exercise6) had positive associations with risk-taking and negative associations with neuroticism. Other associations with physical activity habits6 and alcohol consumption4 demonstrate less consistent associations with risk-taking.
Discussion
We identified novel associations between CADM2 genetic variants and SBP, CRP levels, neuroticism and mood instability, and have highlighted a possible link between SNPs associated with psychological traits and adiposity via CADM2 expression levels in adipose tissue.
Associations between CADM2 SNPs and obesity have previously been reported10,28,40–42,44–47 and were observed here. It is possible that the associations of CADM2 SNPs with CRP and SBP are secondary to the effects on obesity, as increased fat accumulation is associated with systemic inflammation2 and reduced cardiovascular fitness48. Associations between CADM2 SNPs and risk-taking behaviour have also been reported7,8,25. We previously suggested that the association between risk-taking and obesity might be behavioural, with risk-takers choosing to disregard health-related advice and/or are prone to aberrant reward circuitry predisposing them to poor dietary choices and excessive intake7.
Pleiotropy, where genetic variants influence more than one trait, is a concern in genetic studies. Pleiotropy can be classified as biological (where a genetic variant has true effects on more than one trait), mediated (where a genetic variant has a true effect on one trait, but because there is a causal relationship between that trait and a second trait, effects of the genetic variant are seen on the second trait as well) or spurious (where biases in the study result in genetic effects on multiple traits)49. When considering the effects of CADM2 variants on both psychological and obesity traits, a possible explanation is mediated pleiotropy, for example through physical exercise. If the CADM2 variant effects on behavioural traits such as risk-taking, neuroticism, mood instability (observed here) and physical exercise6 are true, then the effects on obesity might be knock-on effects of physical exercise. Whilst possible, the CADM2 variants associated with increased physical activity were associated with increased BMI6, therefore this logic of this argument is flawed. Spurious pleiotropy is another possibility; however there are consistent effects of CADM2 on psychological and obesity traits in a number of cohorts with different recruitment and study designs (population-based, CVD case-control, high CVD risk) and populations (European, UK, north American, Pakistan), which would be expected to differ in their biases. Whether spurious pleiotropy can result from such a variety of biases is doubtful.
In contrast, biological pleiotropy is supported by the body of evidence indicating a role for CADM2 in psychological and obesity traits from other types of studies. Mouse models demonstrate a clear effect of Cadm2 on obesity and gluco-metabolic parameters: A global Cadm2-knockout mouse demonstrated reduced body weight, improved insulin sensitivity and improved glucose tolerance50. Furthermore, this effect was maintained when Cadm2-knockout mice were crossed with the traditional obesity model, the Leptin-knockout mouse11. Rat models also demonstrate that increased Cadm2 (via a knockdown of a cadm2 regulator) reduced neurite outgrowth in response to ischemic damage51 whilst errors in axon pathfinding and neurite outgrowth were observed in Cadm2-deficient chick embryos52. In vivo and in vitro studies of tumour models have demonstrated that increased expression of CADM2 mRNA or protein is associated with reduced cell viability, proliferation, migration and invasion in glioma53, retinoblastoma54, renal cell55, hepatocellular56,57, endometrial58, prostate59 and oesophageal squamous cell carcinomas60. The wide range of cell types suggests that CADM2 regulation of cell turnover could be ubiquitous. With this in mind, neuronal remodelling is important for health and pathology and requires cell turnover, so levels of CADM2 likely influence plasticity of the brain. Similarly, to increase the fat that can be stored, adipocytes either increase in volume (metabolically detrimental) or in number (metabolically benign)61. CADM2 levels could be a part of the volume vs number fate determination.
Further support for biological pleiotropy comes from the observation that CADM2 SNPs associated with risk-taking (rs4856591), neuroticism (rs818219) and mood instability (rs818225) had eQTLs for adipose tissue (and therefore potentially direct effects on adiposity). This is especially interesting in light of the established risk of obesity in psychiatric disorders and the relevance of these traits to a wider range of psychiatric disorders (MDD, GAD, SCZ, BD) than risk-taking (SCZ and BD). A common biological mechanism, such as CADM2, might also be consistent with the recent observation of a bi-directional link between depression and obesity62.
This study did not find evidence for effects of CADM2 on psychiatric diagnoses, although it should be noted that the number of cases for these analyses were low. We also cannot exclude psychiatric medication as a confounder in our analyses of cardiometabolic variables, however the very low percentage of individuals on psychiatric medication means that this is very unlikely. No associations were identified for CVD, which is unlikely to be due to the number of cases present. We also note that, in comparison to the Cadm2-knockout mice, no effects on glucose-related traits were observed. This may be due to a smaller sample size (N = 10,128) and thus reduced power for these biomarkers, or selection bias due to studying cardiovascular cohorts. It is also possible that effects of CADM2 on insulin sensitivity and glucose levels are secondary to effects on BMI, making it harder to discern.
A surprising finding of this study was the long-range linkage disequilibrium within the CADM2 locus, despite low LD, with effects being evident over a region of nearly 1 Mb. This means that the SNPs identified for psychological and metabolic traits are not independent. Haplotype analysis would be of value here, however the standard approaches rely on higher LD and smaller regions, so this has yet to be attempted. This consideration, combined with a plethora of eQTLs in a variety of tissues, results in complexity regarding the regulation of CADM2 levels. Further functional investigation of CADM2 is required for complete mechanistic understanding of this locus.
There are some limitations to this study, notably incomplete genetic coverage of the locus in the IMPROVE and SCARFSHEEP cohorts and reduced sample sizes for the biomarker analyses. In addition, only the UK Biobank had both psychological and cardio-metabolic phenotyping. As is typical for the majority of cardio-metabolic studies, history of psychiatric illness was an exclusion criterion for IMPROVE, SCARFSHEEP and PROCARDIS, therefore it is possible that these cohorts have lower levels of variants associated with psychiatric disorders than the general population. Conversely, these cohorts have higher rates of cardiovascular risk factors and disease than the general population (average/general population, to moderate/early CVD, to high risk/at least 3 CVD risk factors). Whilst this is a strength when looking at cardiovascular phenotypes, it is likely to contribute to the high I2 demonstrated for some SNPs in the meta-analyses. Despite this, there are a number of variants that show significant associations, with I2 = 0, and these were used for the follow-up analyses. The consistency in effect sizes and directions for these associations with BMI and CRP are striking, especially being irrespective of CVD risk burden. These findings also demonstrate that the effects of CADM2 variants on cardio-metabolic parameters are generalizable to a wider European ancestry population. The same cannot be assumed for the psychological phenotypes, which were analysed exclusively in white British UK Biobank participants. Strengths of the study include large sample sizes for most analyses. The meta-analyses of several cohorts provides robust results, whereas consistent phenotyping is a clear advantage of the UK Biobank study.
In conclusion, we have conducted a systematic, large-scale analysis of multiple datasets providing evidence that CADM2 represents a putative shared biological link between metabolic and psychological disorders. Future work, including animal models which investigate both metabolic and behavioural traits in the same animals, is now needed to understand the functional biological mechanisms that might explain this link.
Supplementary information
Untitled section
Supplementary information
Supplementary information accompanies this paper at 10.1038/s41598-019-43861-9.
Acknowledgements
We thank all participants of the UK Biobank, IMPROVE, SCARFSHEEP and PROCARDIS studies. We acknowledge the contributions of the individuals who initiated and recruited the IMPROVE study (Supplemental Data). The UK Biobank was established by the Wellcome Trust, Medical Research Council, Department of Health, Scottish Government and Northwest Regional Development Agency. UK Biobank has also had funding from the Welsh Assembly Government and the British Heart Foundation. Data collection was funded by UK Biobank. PROCARDIS was supported by the European Community Sixth Framework Program (LSHM-CT- 2007-037273), AstraZeneca, the Swedish Research Council (8691), the Knut and Alice Wallenberg Foundation, the Swedish Heart-Lung Foundation, the Torsten and Ragnar Söderberg Foundation, the Strategic Cardiovascular Program of Karolinska Institutet and Stockholm County Council, the Foundation for Strategic Research and the Stockholm County Council (560283). IMPROVE was funded by the European Commission (LSHM-CT- 2007-037273), the Swedish Heart-Lung Foundation, the Swedish Research Council (8691), the Knut and Alice Wallenberg Foundation, the Foundation for Strategic Research, the Torsten and Ragnar Söderberg Foundation, the Strategic Cardiovascular Programme of Karolinska Institutet and the Stockholm County Council and the Stockholm County Council (560183). S.E.H. is funded by the British Heart Foundation (PG008/08) and is supported by the National Institute for Health Research, University College London Hospitals Biomedical Research Centre. SCARFSHEEP was supported by the Swedish Heart-Lung Foundation, the Swedish Research Council, the Strategic Cardiovascular Programme of Karolinska Institutet and the Stockholm County Council, the Strategic support for epidemiological research at Karolinska Institutet and the Stockholm County Council. JM was supported by a Lister Prize Summer Studentship. JW is supported by the JMAS Sim Fellowship for depression research from the Royal College of Physicians of Edinburgh (173558). AF is supported by an MRC Doctoral Training Programme Studentship at the University of Glasgow (MR/K501335/1). KJAJ is supported by an MRC Doctoral Training Programme Studentship at the Universities of Glasgow and Edinburgh. DJS acknowledges the support of the Brain and Behaviour Research Foundation (Independent Investigator Award 1930) and a Lister Prize Fellowship (173096). RJS is supported by a UKRI Innovation- HDR-UK Fellowship (MR/S003061/1). The funders had no involvement in the study design, analysis, interpretation or manuscript preparation.
Data Availability
Data is available on request, contact either UK Biobank (UK Biobank-only analyses) or the corresponding author (all other analyses).
Notes
Competing Interests
The authors declare no competing interests.