Senior Corresponding Authors: Matthew E. Hurles, The Wellcome Trust Sanger Institute,

Disorders of the brain can exhibit considerable epidemiological comorbidity and often share symptoms, provoking debate about their etiologic overlap. We quantified the genetic sharing of 25 brain disorders from genome-wide association studies of 265,218 patients and 784,643 control participants and assessed their relationship to 17 phenotypes from 1,191,588 individuals. Psychiatric disorders share common variant risk, whereas neurological disorders appear more distinct from one another and from the psychiatric disorders. We also identified significant sharing between disorders and a number of brain phenotypes, including cognitive measures. Further, we conducted simulations to explore how statistical power, diagnostic misclassification, and phenotypic heterogeneity affect genetic correlations. These results highlight the importance of common genetic variation as a risk factor for brain disorders and the value of heritability-based methods in understanding their etiology.

Both underweight and obesity have been associated with increased mortality1,2. Underweight, defined as body mass index (BMI) ≤ 18,5 kg/m2 in adults 3 and ≤ −2 standard deviations (SD) in children4,5, is the main sign of a series of heterogeneous clinical conditions such as failure to thrive (FTT) 6–8, feeding and eating disorder and/or anorexia nervosa9,10. In contrast to obesity, few genetic variants underlying these clinical conditions have been reported 11, 12. We previously demonstrated that hemizygosity of a ~600 kb region on the short arm of chromosome 16 (chr16:29.5–30.1Mb), causes a highly-penetrant form of obesity often associated with hyperphagia and intellectual disabilities13. Here we show that the corresponding reciprocal duplication is associated with underweight. We identified 138 (132 novel cases) duplication carriers (108 unrelated carriers) from over 95,000 individuals clinically-referred for developmental or intellectual disabilities (DD/ID), psychiatric disorders or recruited from population-based cohorts. These carriers show significantly reduced postnatal weight (mean Z-score −0.6; p=4.4×10−4) and BMI (mean Z-score −0.5; p=2.0×10−3). In particular, half of the boys younger than 5 years are underweight with a probable diagnosis of FTT, while adult duplication carriers have an 8.7-fold (p=5.9×10−11; CI_95=[4.5–16.6]) increased risk of being clinically underweight. We observe a significant trend towards increased severity in males, as well as a depletion of male carriers among non-medically ascertained cases. These features are associated with an unusually high frequency of selective and restrictive feeding behaviours and a significant reduction in head circumference (mean Z-score −0.9; p=7.8×10−6). Each of the observed phenotypes is the converse of one reported in carriers of deletions at this locus, correlating with changes in transcript levels for genes mapping within the duplication but not within flanking regions. The reciprocal impact of these 16p11.2 copy number variants suggests that severe obesity and being underweight can have mirror etiologies, possibly through contrasting effects on eating behaviour.

Objective To conduct a genome-wide association study (GWAS) of anorexia nervosa and to calculate genetic correlations with a series of psychiatric, educational, and metabolic phenotypes. Method Following uniform quality control and imputation using the 1000 Genomes Project (phase 3) in 12 case-control cohorts comprising 3,495 anorexia nervosa cases and 10,982 controls, we performed standard association analysis followed by a meta-analysis across cohorts. Linkage disequilibrium score regression (LDSC) was used to calculate genome-wide common variant heritability [ hSNP2, partitioned heritability, and genetic correlations (rg)] between anorexia nervosa and other phenotypes. Results Results were obtained for 10,641,224 single nucleotide polymorphisms (SNPs) and insertion-deletion variants with minor allele frequency > 1% and imputation quality scores > 0.6. The hSNP2 of anorexia nervosa was 0.20 (SE=0.02), suggesting that a substantial fraction of the twin-based heritability arises from common genetic variation. We identified one genome-wide significant locus on chromosome 12 (rs4622308, p=4.3×10−9) in a region harboring a previously reported type 1 diabetes and autoimmune disorder locus. Significant positive genetic correlations were observed between anorexia nervosa and schizophrenia, neuroticism, educational attainment, and high density lipoprotein (HDL) cholesterol, and significant negative genetic correlations between anorexia nervosa and body mass index, insulin, glucose, and lipid phenotypes. Conclusions Anorexia nervosa is a complex heritable phenotype for which we have found the first genome-wide significant locus. Anorexia nervosa also has large and significant genetic correlations with both psychiatric phenotypes and metabolic traits. Our results encourage a reconceptualization of this frequently lethal disorder as one with both psychiatric and metabolic etiology.

No abstract

Anorexia nervosa (AN) is a complex and heritable eating disorder characterized by dangerously low body weight. Neither candidate gene studies nor an initial genome wide association study (GWAS) have yielded significant and replicated results. We performed a GWAS in 2,907 cases with AN from 14 countries (15 sites) and 14,860 ancestrally matched controls as part of the Genetic Consortium for AN (GCAN) and the Wellcome Trust Case Control Consortium 3 (WTCCC3). Individual association analyses were conducted in each stratum and meta-analyzed across all 15 discovery datasets. Seventy-six (72 independent) SNPs were taken forward for in silico (two datasets) or de novo (13 datasets) replication genotyping in 2,677 independent AN cases and 8,629 European ancestry controls along with 458 AN cases and 421 controls from Japan. The final global meta-analysis across discovery and replication datasets comprised 5,551 AN cases and 21,080 controls. AN subtype analyses (1,606 AN restricting; 1,445 AN binge-purge) were performed. No findings reached genome-wide significance. Two intronic variants were suggestively associated: rs9839776 (P=3.01×10-7) in SOX2OT and rs17030795 (P=5.84×10-6) in PPP3CA. Two additional signals were specific to Europeans: rs1523921 (P=5.76×10-6) between CUL3 and FAM124B and rs1886797 (P=8.05×10-6) near SPATA13. Comparing discovery to replication results, 76% of the effects were in the same direction, an observation highly unlikely to be due to chance (P=4×10-6), strongly suggesting that true findings exist but that our sample, the largest yet reported, was underpowered for their detection. The accrual of large genotyped AN case-control samples should be an immediate priority for the field.

Disruptions in circadian rhythms have been described in mood disorders (MD), but the involvement of genetic variation in genes pertaining to the molecular circadian machinery in the susceptibility to MD has not been conclusively determined. We examined 209 single-nucleotide polymorphisms (SNPs) covering 19 circadian genes (ADCYAP1, ARNTL, ARNTL2, BHLHB2, BHLHB3, CLOCK, CRY1, CRY2, CSNK1E, DBP, NPAS2, NR1D1, PER1, PER2, PER3, RORA, TIMELESS, VIP, and VIPR2) in a sample of 534 MD patients (335 with unipolar major mood depression (MDD) and 199 with bipolar disorder (BD)) and 440 community-based screened controls. Nominally, statistically significant associations were found in 15 circadian genes. The gene-wide test, corrected for the number of SNPs analyzed in each gene, identified significant associations in CRY1 (rs2287161), NPAS2 (rs11123857), and VIPR2 (rs885861) genes with the combined MD sample. In the MDD subsample, the same SNPs in CRY1 and NPAS2 of the combined sample remained associated, whereas in the BD subsample CLOCK (rs10462028) and VIP (rs17083008) were specifically associated. The association with an SNP located 3 0 near CRY1 gene in MDD remained statistically significant after permutation correction at experiment level (p ¼ 0.007). Significant additive effects were found between the SNPs that were statistically significant at the gene-wide level. We also found evidence of associations between two-marker haplotypes in CRY1 and NPAS2 genes and MD. Our data support the contribution of the circadian system to the genetic susceptibility to MD and suggest that different circadian genes may have specific effects on MD polarity.

Several genes with an essential role in the regulation of eating behavior and body weight are considered candidates involved in the etiology of eating disorders (ED), but no relevant susceptibility genes with a major effect on anorexia nervosa (AN) or bulimia nervosa (BN) have been identified. Brain-derived neurotrophic factor (BDNF) has been implicated in the regulation of food intake and body weight in rodents. We previously reported a strong association of the Met66 allele of the Val66Met BDNF variant with restricting AN (ANR) and low minimum body mass index in Spanish patients. Another single nucleotide polymorphism located in the promoter region of the BDNF gene (-270C>T) showed lack of association with any ED phenotype. In order to replicate these findings in a larger sample, we performed a case-control study in 1142 Caucasian patients with ED consecutively recruited in six different centers from five European countries (France, Germany, Italy, Spain and UK) participating in the 'Factors in Healthy Eating' project. We have found that the Met66 variant is strongly associated to all ED subtypes (AN, ANR, binge-eating/purging AN and BN), and that the -270C BDNF variant has an effect on BN and late age at onset of weight loss. These are the first two variants associated with the pathophysiology of ED in different populations and support a role for BDNF in the susceptibility to aberrant eating behaviors.