Association study of 182 candidate genes in anorexia nervosa

Anti-hypothalamus autoantibodies in anorexia nervosa: a possible new mechanism in neuro-physiological derangement?

PURPOSE

Anorexia nervosa (AN) is a serious and complex mental disorder affecting mainly young adult women. AN patients are characterized by low body weight in combination with self-induced starvation, intense fear of gaining weight, and distortion of body image. AN is a multifactorial disease, linked by recent evidence to a dysregulation of the immune system.

METHODS

In this pilot study, 22 blood serums from AN patients were tested for the presence of autoantibodies against primate hypothalamic periventricular neurons by immunofluorescence and by a home-made ELISA assay. Cellular fluorescence suggests the presence of autoantibodies which are able to recognize these neurons (both to body cell and fiber levels). By means of ELISA, these autoantibodies are quantitatively evaluated. In addition, orexigenic and anorexigenic molecules were measured by ELISA. As control, 18 blood serums from healthy age matched woman were analysed.

RESULTS

All AN patients showed a reactivity against hypothalamic neurons both by immunofluorescence and ELISA. In addition, ghrelin, pro-opiomelanocortin (POMC), and agouti-related peptide (AGRP) were significantly higher than in control serums (p < 0.0001). In contrast, leptin was significantly lower in AN patients than controls (p < 0.0001).

CONCLUSIONS

Immunoreaction and ELISA assays on AN blood serum suggest the presence of autoantibodies AN related. However, it is not easy to determine the action of these antibodies in vivo: they could interact with specific ligands expressed by hypothalamic cells preventing their physiological role, however, it is also possible that they could induce an aspecific stimulation in the target cells leading to an increased secretion of anorexigenic molecules. Further studies are needed to fully understand the involvement of the immune system in AN pathogenesis.

LEVEL OF EVIDENCE

V, descriptive study.

Childhood trauma determines different clinical and biological manifestations in patients with eating disorders

PURPOSE

There is a significant relationship between childhood trauma and the development of an eating disorder in adolescence or adulthood, possibly influenced by circulating levels of inflammatory parameters. The main objective is to identify and describe a subgroup of patients with eating disorders and a history of trauma in childhood or adolescence with differential clinical features.

METHODS

An observational study on a sample of 55 patients who met the diagnostic criteria for any DSM-5 eating disorder was carried out. Inflammatory parameters in white blood cells were examined. Patients underwent different assessments, including clinical and personality scales.

RESULTS

Patients with a history of trauma had higher scores in the delirious and narcissistic items of the Millon Clinical Multiaxial Inventory (MCMI-II) (p < 0.05) and a higher score in the paranoid item of the SCID-5 Personality Disorders Version (SCID-5-PD) (p < 0.05). Patients with distinguishing personality features were grouped according to the Childhood Trauma Questionnaire sexual subscale. Tumor necrosis alpha (TNF-α) showed a significant association with childhood trauma history.

CONCLUSIONS

There is a profile of patients with eating disorders who have increased activity in the inflammatory pathways that, if identified precociously, can benefit from specifically aimed interventions.

LEVEL OF EVIDENCE

Level V, observational study.

Variants in the Obesity-Linked FTO gene locus modulates psychopathological features of patients with Anorexia Nervosa

Our aim was to determine whether variability in the fat mass obesity (FTO) gene locus, consistently related to obesity, affects the risk of eating disorders (ED) and/or the psychopathology displayed by these patients. We analyzed 26 tag-single nucleotide polymorphisms (SNPs) that capture FTO variability in 352 ED patients [233 with Anorexia Nervosa (AN) and 119 with binge-eating] and 396 controls. Psychopathological symptoms and traits were evaluated by the Eating Disorders Inventory Test 2 (EDI-2) and Symptoms Checklist 90 Revised (SCL-90R) questionnaires. No associations were found for ED risk. The rs7205987 CC genotype correlated with higher scores in all but one of the EDI-2 scales in the AN group. Associations with Bulimia (p = 0.0019) and Interoceptive Awareness (p = 0.00007) retained significance after False Discovery Rate (FDR) correction for multiple testing. A 3-SNP sliding window analysis showed that FTO haplotypes were again highly associated with Interoceptive Awareness (rs9921255/rs6499662/rs7205987 haplotype; FDR-q = 0.04), Bulimia (rs1125338/rs2192872/rs708258; FDR-q = 0.00037), and Maturity Fears (rs708258/rs12599672/rs11076017; FDR-q = 0.041). In addition, a distal region of the gene between rs9924877 (position 53947509) and rs2192872 (54040715) was linked to Anxiety, Depression and Phobic Anxiety in AN patients, with FDR-q values ranging from 0.023 to 0.045. The results suggest that the FTO gene might be an important locus regarding traits and psychopathological symptoms often displayed by AN patients.

Genetic risk factors for eating disorders: an update and insights into pathophysiology

Genome-wide-association studies (GWASs), epigenetic, gene-expression and gene-gene interaction projects, nutritional genomics and investigations of the gut microbiota have increased our knowledge of the pathophysiology of eating disorders (EDs). However, compared with anorexia nervosa, genetic studies in patients with bulimia nervosa and binge-eating disorder are relatively scarce, with the exception of a few formal genetic and small-sized candidate-gene-association studies. In this article, we review important findings derived from formal and molecular genetics in order to outline a genetics-based pathophysiological model of EDs. This model takes into account environmental and nutritional factors, genetic factors related to the microbiome, the metabolic and endocrine system, the immune system, and the brain, in addition to phenotypical traits of EDs. Shortcomings and advantages of genetic research in EDs are discussed against the historical background, but also in light of potential future treatment options for patients with EDs.

Detecting significant genotype-phenotype association rules in bipolar disorder: market research meets complex genetics

Background

Disentangling the etiology of common, complex diseases is a major challenge in genetic research. For bipolar disorder (BD), several genome-wide association studies (GWAS) have been performed. Similar to other complex disorders, major breakthroughs in explaining the high heritability of BD through GWAS have remained elusive. To overcome this dilemma, genetic research into BD, has embraced a variety of strategies such as the formation of large consortia to increase sample size and sequencing approaches. Here we advocate a complementary approach making use of already existing GWAS data: a novel data mining procedure to identify yet undetected genotype–phenotype relationships. We adapted association rule mining, a data mining technique traditionally used in retail market research, to identify frequent and characteristic genotype patterns showing strong associations to phenotype clusters. We applied this strategy to three independent GWAS datasets from 2835 phenotypically characterized patients with BD. In a discovery step, 20,882 candidate association rules were extracted.

Results

Two of these rules—one associated with eating disorder and the other with anxiety—remained significant in an independent dataset after robust correction for multiple testing. Both showed considerable effect sizes (odds ratio ~ 3.4 and 3.0, respectively) and support previously reported molecular biological findings.

Conclusion

Our approach detected novel specific genotype–phenotype relationships in BD that were missed by standard analyses like GWAS. While we developed and applied our method within the context of BD gene discovery, it may facilitate identifying highly specific genotype–phenotype relationships in subsets of genome-wide data sets of other complex phenotype with similar epidemiological properties and challenges to gene discovery efforts.

Electronic supplementary material

The online version of this article (10.1186/s40345-018-0132-x) contains supplementary material, which is available to authorized users.

The Anorexia Nervosa Genetics Initiative (ANGI): Overview and methods

BACKGROUND

Genetic factors contribute to anorexia nervosa (AN); and the first genome-wide significant locus has been identified. We describe methods and procedures for the Anorexia Nervosa Genetics Initiative (ANGI), an international collaboration designed to rapidly recruit 13,000 individuals with AN and ancestrally matched controls. We present sample characteristics and the utility of an online eating disorder diagnostic questionnaire suitable for large-scale genetic and population research.

METHODS

ANGI recruited from the United States (US), Australia/New Zealand (ANZ), Sweden (SE), and Denmark (DK). Recruitment was via national registers (SE, DK); treatment centers (US, ANZ, SE, DK); and social and traditional media (US, ANZ, SE). All cases had a lifetime AN diagnosis based on DSM-IV or ICD-10 criteria (excluding amenorrhea). Recruited controls had no lifetime history of disordered eating behaviors. To assess the positive and negative predictive validity of the online eating disorder questionnaire (ED100K-v1), 109 women also completed the Structured Clinical Interview for DSM-IV (SCID), Module H.

RESULTS

Blood samples and clinical information were collected from 13,363 individuals with lifetime AN and from controls. Online diagnostic phenotyping was effective and efficient; the validity of the questionnaire was acceptable.

CONCLUSIONS

Our multi-pronged recruitment approach was highly effective for rapid recruitment and can be used as a model for efforts by other groups. High online presence of individuals with AN rendered the Internet/social media a remarkably effective recruitment tool in some countries. ANGI has substantially augmented Psychiatric Genomics Consortium AN sample collection. ANGI is a registered clinical trial: clinicaltrials.govNCT01916538; https://clinicaltrials.gov/ct2/show/NCT01916538?cond=Anorexia+Nervosa&draw=1&rank=3.

Novel and ultra-rare damaging variants in neuropeptide signaling are associated with disordered eating behaviors

Objective

Eating disorders develop through a combination of genetic vulnerability and environmental stress, however the genetic basis of this risk is unknown.

Methods

To understand the genetic basis of this risk, we performed whole exome sequencing on 93 unrelated individuals with eating disorders (38 restricted-eating and 55 binge-eating) to identify novel damaging variants. Candidate genes with an excessive burden of predicted damaging variants were then prioritized based upon an unbiased, data-driven bioinformatic analysis. One top candidate pathway was empirically tested for therapeutic potential in a mouse model of binge-like eating.

Results

An excessive burden of novel damaging variants was identified in 186 genes in the restricted-eating group and 245 genes in the binge-eating group. This list is significantly enriched (OR = 4.6, p<0.0001) for genes involved in neuropeptide/neurotrophic pathways implicated in appetite regulation, including neurotensin-, glucagon-like peptide 1- and BDNF-signaling. Administration of the glucagon-like peptide 1 receptor agonist exendin-4 significantly reduced food intake in a mouse model of ‘binge-like’ eating.

Conclusions

These findings implicate ultra-rare and novel damaging variants in neuropeptide/neurotropic factor signaling pathways in the development of eating disorder behaviors and identify glucagon-like peptide 1-receptor agonists as a potential treatment for binge eating.

Convergent and divergent genetic changes in the genome of Chinese and European pigs

Since 10,000 BC, continuous human selection has led to intense genetic and phenotypic changes in pig (Sus scrofa) domestication. Through whole genome analysis of 257 individuals, we demonstrated artificial unidirectional and bidirectional selection as the primary force to shape the convergent and divergent changes between Chinese domestic pigs (CHD) and European domestic pigs (EUD). We identified 31 genes in unidirectional selection regions that might be related to fundamental domestication requirements in pigs. And these genes belong predominantly to categories related to the nervous system, muscle development, and especially to metabolic diseases. In addition, 35 genes, representing different breeding preference, were found under bidirectional selection for the distinct leanness and reproduction traits between CHD and EUD. The convergent genetic changes, contributing physical and morphological adaption, represent the common concerns on pig domestication. And the divergent genetic changes reflect distinct breeding goals between Chinese and European pigs. Using ITPR3, AHR and NMU as examples, we explored and validated how the genetic variations contribute to the phenotype changes.

A novel relationship for schizophrenia, bipolar and major depressive disorder Part 5: a hint from chromosome 5 high density association screen

Familial clustering of schizophrenia (SCZ), bipolar disorder (BPD), and major depressive disorder (MDD) was systematically reported (Aukes, M. F. Genet Med 2012, 14, 338-341) and any two or even three of these disorders could co-exist in some families. In addition, evidence from symptomatology and psychopharmacology also imply that there are intrinsic connections between these three major disorders. A total of 56,569 single nucleotide polymorphism (SNPs) on chromosome 5 were genotyped by Affymetrix Genome-Wide Human SNP array 6.0 on 119 SCZ, 253 BPD (type-I), 177 MDD patients and 1000 controls. Associated SNPs and flanking genes was screen out systematically, and cadherin pathway genes (CDH6, CDH9, CDH10, CDH12, and CDH18) belong to outstanding genes. Unexpectedly, nearly all flanking genes of the associated SNPs distinctive for BPD and MDD were replicated in an enlarged cohort of 986 SCZ patients (P ≤ 9.9E-8). Considering multiple bits of evidence, our chromosome 5 analyses implicated that bipolar and major depressive disorder might be subtypes of schizophrenia rather than two independent disease entities. Also, cadherin pathway genes play important roles in the pathogenesis of the three major mental disorders.

The Eating-Disorder Associated HDAC4A778T Mutation Alters Feeding Behaviors in Female Mice

BACKGROUND

While eating disorders (EDs) are thought to result from a combination of environmental and psychological stressors superimposed on genetic vulnerability, the neurobiological basis of EDs remains incompletely understood. We recently reported that a rare missense mutation in the gene for the transcriptional repressor histone deacetylase 4 (HDAC4) is associated with the risk of developing an ED in humans.

METHODS

To understand the biological consequences of this missense mutation, we created transgenic mice carrying this mutation by introducing the alanine to threonine mutation at position 778 of mouse Hdac4 (corresponding to position 786 of the human protein). Bioinformatic analysis to identify Hdac4-regulated genes was performed using available databases.

RESULTS

Male mice heterozygous for HDAC4^A778T did not show any metabolic or behavioral differences. In contrast, female mice heterozygous for HDAC4^A778T display several ED-related feeding and behavioral deficits depending on housing condition. Individually housed HDAC4^A778T female mice exhibit reduced effortful responding for high-fat diet and compulsive grooming, whereas group-housed female mice display increased weight gain on high-fat diet, reduced behavioral despair, and increased anxiety-like behaviors. Bioinformatic analysis identifies mitochondrial biogenesis including synthesis of glutamate/gamma-aminobutyric acid as a potential transcriptional target of HDAC4^A778T activity relevant to the behavioral deficits identified in this new mouse model of disordered eating.

CONCLUSIONS

The HDAC4^A778T mouse line is a novel model of ED-related behaviors and identifies mitochondrial biogenesis as a potential molecular pathway contributing to behavioral deficits.

Integrating multi-omics biomarkers and postprandial metabolism to develop personalized treatment for anorexia nervosa

BACKGROUND

Anorexia Nervosa (AN) is a serious mental illness characterized by emaciation, an intense fear of gaining weight despite being underweight, and distorted body image. Few treatments reverse the core symptoms in AN such as profound aversion to food and food avoidance. Consequently, AN has a chronic and relapsing course and the highest mortality rate of any psychiatric illness. A more complete understanding of the disease pathogenesis is needed in order to develop better treatments and improve AN outcome. The pathogenesis and psychopathophysiology of AN can be better elucidated by combining longitudinal phenotyping with multiple "omics" techniques, including genomics, proteomics, lipidomics, and metabolomics.

DESIGN

This paper summarizes the key findings of a series of interrelated studies including new experimental data and previously published data, and describes our current initiatives and future directions.

RESULTS

Exon sequencing data was analyzed in 1205 AN and 1948 controls. Targeted metabolomics, lipidomics, and proteomics data were collected in two independent convenience samples consisting of 75 subjects with eating disorders and 61 sex- and age-matched healthy controls. Study participants were female and the mean age was 22.9 (4.9 [SD]) years. Epoxide hydrolase 2 (EPHX2) genetic variations were significantly associated with AN risk, and epoxide hydrolase (sEH) activity was elevated in AN compared to controls. The polyunsaturated fatty acids (PUFAs) and eicosanoids data revealed that cytochrome P450 pathway was implicated in AN, and AN displayed a dysregulated postprandial metabolism of PUFAs and sEH-dependent eicosanoids.

IMPLICATION AND CURRENT INITIATIVES

Collectively, our data suggest that dietary factors may contribute to the burden of EPHX2-associated AN susceptibility and affect disease outcome. We are implementing new investigations using a longitudinal study design in order to validate and develop an EPHX2 multi-omics biomarker system. We will test whether sEH-associated postprandial metabolism increases AN risk and affects treatment outcome through an ω-6 rich breakfast challenge. Participants will include 100 ill AN patients, 100 recovered AN patients, and 100 age- and race-matched healthy women. These data will allow us to investigate 1) how genetic and dietary factors independently and synergistically contribute to AN risk and progression, and 2) if clinical severity and treatment response in AN are affected by sEH activity and eicosanoid dysregulation. Results of our study will 1) identify clinically relevant biomarkers, 2) unravel mechanistic functions of sEH, and 3) delineate contributory roles of dietary PUFAs and cytochrome P450 pathway eicosanoids for the purpose of developing novel AN treatments and improving disease prognosis.

The First Pilot Genome-Wide Gene-Environment Study of Depression in the Japanese Population

Stressful events have been identified as a risk factor for depression. Although gene–environment (G × E) interaction in a limited number of candidate genes has been explored, no genome-wide search has been reported. The aim of the present study is to identify genes that influence the association of stressful events with depression. Therefore, we performed a genome-wide G × E interaction analysis in the Japanese population. A genome-wide screen with 320 subjects was performed using the Affymetrix Genome-Wide Human Array 6.0. Stressful life events were assessed using the Social Readjustment Rating Scale (SRRS) and depression symptoms were assessed with self-rating questionnaires using the Center for Epidemiologic Studies Depression (CES-D) scale. The p values for interactions between single nucleotide polymorphisms (SNPs) and stressful events were calculated using the linear regression model adjusted for sex and age. After quality control of genotype data, a total of 534,848 SNPs on autosomal chromosomes were further analyzed. Although none surpassed the level of the genome-wide significance, a marginal significant association of interaction between SRRS and rs10510057 with depression were found (p = 4.5 × 10⁻⁸). The SNP is located on 10q26 near Regulators of G-protein signaling 10 (RGS10), which encodes a regulatory molecule involved in stress response. When we investigated a similar G × E interaction between depression (K6 scale) and work-related stress in an independent sample (n = 439), a significant G × E effect on depression was observed (p = 0.015). Our findings suggest that rs10510057, interacting with stressors, may be involved in depression risk. Incorporating G × E interaction into GWAS can contribute to find susceptibility locus that are potentially missed by conventional GWAS.

New Insights in Anorexia Nervosa

Anorexia nervosa (AN) is classically defined as a condition in which an abnormally low body weight is associated with an intense fear of gaining weight and distorted cognitions regarding weight, shape, and drive for thinness. This article reviews recent evidences from physiology, genetics, epigenetics, and brain imaging which allow to consider AN as an abnormality of reward pathways or an attempt to preserve mental homeostasis. Special emphasis is put on ghrelino-resistance and the importance of orexigenic peptides of the lateral hypothalamus, the gut microbiota and a dysimmune disorder of neuropeptide signaling. Physiological processes, secondary to underlying, and premorbid vulnerability factors-the "pondero-nutritional-feeding basements"- are also discussed.

Catechol-O-methyltransferase activity in erythrocytes from patients with eating disorders

PURPOSE

Abnormal feeding has been linked to disruptions in brain dopaminergic activity and recent studies have assessed the role of catechol-O-methyltransferase (COMT) in eating disorders. This is the first study to quantify the soluble catechol-O-methyltransferase (S-COMT) activity in erythrocytes from patients with anorexia nervosa (AN), bulimia nervosa (BN) and binge-eating disorder (BED) and the first study at all to evaluate the COMT on patients with BED.

METHODS

Forty blood samples from patients with AN, BN and BED and healthy controls were drawn to evaluate S-COMT activity in erythrocytes by high-performance liquid chromatography and mass spectrometry. Since several patients were being treated with fluoxetine 20 mg, they were included in a different group (BN MED and BED MED). Liver homogenates from rats were used to evaluate baseline S-COMT activity in the presence of fluoxetine by the same in vitro procedures and assays.

RESULTS

Erythrocyte S-COMT activity (pmol/mg prt/h) was significantly increased in patients with BN and BED (41.3 ± 6.8 and 41.4 ± 14, respectively) compared to control group (25.3 ± 9.7). In fluoxetine-treated patients with BN, S-COMT activity (15.9 ± 8.8) was decreased compared to the other BN group; however, in BED group, the difference between BED MED and BED was not observed. In patients with AN, no significant difference was found compared to controls.

CONCLUSION

Patients with BN and BED presented higher S-COMT activity in erythrocytes, which is in agreement with previous studies on the literature addressing the high-activity COMT allele, Val158, as risk factor for eating disorders. Although in fluoxetine-treated patients with BN the activity of S-COMT was similar to the controls, this is not explained by a direct interaction between fluoxetine and S-COMT as verified in in vitro assays.

BDNF-Val66Met variant and adolescent stress interact to promote susceptibility to anorexic behavior in mice

There is an urgent need to identify therapeutic targets for anorexia nervosa (AN) because current medications do not impact eating behaviors that drive AN's high mortality rate. A major obstacle to developing new treatments is the lack of animal models that recapitulate the pattern of disease onset typically observed in human populations. Here we describe a translational mouse model to study interactions between genetic, psychological and biological risk factors that promote anorexic behavior. We combined several factors that are consistently associated with increased risk of AN-adolescent females, genetic predisposition to anxiety imposed by the BDNF-Val66Met gene variant, social isolation stress and caloric restriction (CR). Approximately 40% of the mice with all of these risk factors will exhibit severe self-imposed dietary restriction, sometimes to the point of death. We systematically varied the risk factors outlined above to explore how they interact to influence anorexic behavior. We found that the Val66Met genotype markedly increases the likelihood and severity of abnormal feeding behavior triggered by CR, but only when CR is imposed in the peri-pubertal period. Incidence of anorexic behavior in our model is dependent on juvenile exposure to social stress and can be extinguished by adolescent handling, but is discordant from anxiety-like behavior. Thus, this study characterized gene × environment interactions during adolescence that could be the underlying driver of abnormal eating behavior in certain AN patients, and represents a promising system to identify possible targets for therapeutic intervention.

Dysregulation of soluble epoxide hydrolase and lipidomic profiles in anorexia nervosa

Individuals with anorexia nervosa (AN) restrict eating and become emaciated. They tend to have an aversion to foods rich in fat. Because epoxide hydrolase 2 (EPHX2) was identified as a novel AN susceptibility gene, and because its protein product, soluble epoxide hydrolase (sEH), converts bioactive epoxides of polyunsaturated fatty acid (PUFA) to the corresponding diols, lipidomic and metabolomic targets of EPHX2 were assessed to evaluate the biological functions of EPHX2 and their role in AN. Epoxide substrates of sEH and associated oxylipins were measured in ill AN, recovered AN and gender- and race-matched controls. PUFA and oxylipin markers were tested as potential biomarkers for AN. Oxylipin ratios were calculated as proxy markers of in vivo sEH activity. Several free- and total PUFAs were associated with AN diagnosis and with AN recovery. AN displayed elevated n-3 PUFAs and may differ from controls in PUFA elongation and desaturation processes. Cytochrome P450 pathway oxylipins from arachidonic acid, linoleic acid, alpha-linolenic acid and docosahexaenoic acid PUFAs are associated with AN diagnosis. The diol:epoxide ratios suggest the sEH activity is higher in AN compared with controls. Multivariate analysis illustrates normalization of lipidomic profiles in recovered ANs. EPHX2 influences AN risk through in vivo interaction with dietary PUFAs. PUFA composition and concentrations as well as sEH activity may contribute to the pathogenesis and prognosis of AN. Our data support the involvement of EPHX2-associated lipidomic and oxylipin dysregulations in AN, and reveal their potential as biomarkers to assess responsiveness to future intervention or treatment.

Machine learning derived risk prediction of anorexia nervosa

Background

Anorexia nervosa (AN) is a complex psychiatric disease with a moderate to strong genetic contribution. In addition to conventional genome wide association (GWA) studies, researchers have been using machine learning methods in conjunction with genomic data to predict risk of diseases in which genetics play an important role.

Methods

In this study, we collected whole genome genotyping data on 3940 AN cases and 9266 controls from the Genetic Consortium for Anorexia Nervosa (GCAN), the Wellcome Trust Case Control Consortium 3 (WTCCC3), Price Foundation Collaborative Group and the Children’s Hospital of Philadelphia (CHOP), and applied machine learning methods for predicting AN disease risk. The prediction performance is measured by area under the receiver operating characteristic curve (AUC), indicating how well the model distinguishes cases from unaffected control subjects.

Results

Logistic regression model with the lasso penalty technique generated an AUC of 0.693, while Support Vector Machines and Gradient Boosted Trees reached AUC’s of 0.691 and 0.623, respectively. Using different sample sizes, our results suggest that larger datasets are required to optimize the machine learning models and achieve higher AUC values.

Conclusions

To our knowledge, this is the first attempt to assess AN risk based on genome wide genotype level data. Future integration of genomic, environmental and family-based information is likely to improve the AN risk evaluation process, eventually benefitting AN patients and families in the clinical setting.

Electronic supplementary material

The online version of this article (doi:10.1186/s12920-016-0165-x) contains supplementary material, which is available to authorized users.

Canonical genetic signatures of the adult human brain

The structure and function of the human brain are highly stereotyped, implying a conserved molecular program responsible for its development, cellular structure and function. We applied a correlation-based metric called differential stability to assess reproducibility of gene expression patterning across 132 structures in six individual brains, revealing mesoscale genetic organization. The genes with the highest differential stability are highly biologically relevant, with enrichment for brain-related annotations, disease associations, drug targets and literature citations. Using genes with high differential stability, we identified 32 anatomically diverse and reproducible gene expression signatures, which represent distinct cell types, intracellular components and/or associations with neurodevelopmental and neurodegenerative disorders. Genes in neuron-associated compared to non-neuronal networks showed higher preservation between human and mouse; however, many diversely patterned genes displayed marked shifts in regulation between species. Finally, highly consistent transcriptional architecture in neocortex is correlated with resting state functional connectivity, suggesting a link between conserved gene expression and functionally relevant circuitry.

Genetics and Epigenetics of Eating Disorders

Eating disorders (EDs) are serious psychiatric conditions influenced by biological, psychological, and sociocultural factors. A better understanding of the genetics of these complex traits and the development of more sophisticated molecular biology tools have advanced our understanding of the etiology of EDs. The aim of this review is to critically evaluate the literature on the genetic research conducted on three major EDs: anorexia nervosa (AN), bulimia nervosa (BN), and binge eating disorder (BED). We will first review the diagnostic criteria, clinical features, prevalence, and prognosis of AN, BN, and BED, followed by a review of family, twin, and adoption studies. We then review the history of genetic studies of EDs covering linkage analysis, candidate gene association studies, genome-wide association studies, and the study of rare variants in EDs. Our review also incorporates a translational perspective by covering animal models of ED-related phenotypes. Finally, we review the nascent field of epigenetics of EDs and a look forward to future directions for ED genetic research.

Impact of NEGR1 genetic variability on psychological traits of patients with eating disorders

Genetics variants in the NEGR1 gene, strongly expressed in the brain, have been reported to affect the neuronal control of food intake therefore inducing obesity. With the same rationale, we hypothesized that this genetic variability may be associated with psychological traits commonly displayed by eating disorder (ED) patients and/or with the risk for the disorder. We analyzed 21 tag-single-nucleotide polymorphisms (SNPs) in the coding sequence and adjacent regions of the NEGR1 gene. A total of 169 ED patients (106 with anorexia nervosa (AN) and 63 with bulimia nervosa (BN)) and 312 healthy subjects were genotyped. Personality traits and general psychopathological symptoms were assessed by the Eating Disorders Inventory Test-2 (EDI-2) and Symptom Checklist 90 Revised inventories. None of the SNPs or haplotypes analyzed were associated with a greater risk of ED or correlated with anthropometric parameters. However, in patients with BN, four SNPs (rs12740031, rs10789322, rs6659202 and rs591540) correlated with the scores in Drive for Thinness (DT), Ineffectiveness (I) and Interoceptive Awareness (IA) (Bonferroni-P<0.05 in all instances). The first two SNPs along with rs954299 and rs2422021 formed a haplotype block, which showed a consistent association with the EDI-2 score in BN patients (Bonferroni-P=0.01). A subsequent three-SNP sliding-window approach identified a central area, encompassing both the haplotype block and the individually relevant SNPs that strongly correlated with the scores of BN patients in DT, I, IA and Bulimia. No associations were identified in the AN group. These preliminary results indicate that NEGR1 could be an important locus influencing certain personality dimensions in BN patients.

Genetic neuropathology of obsessive psychiatric syndromes

Anorexia nervosa (AN), bulimia nervosa (BN) and obsessive-compulsive disorder (OCD) are complex psychiatric disorders with shared obsessive features, thought to arise from the interaction of multiple genes of small effect with environmental factors. Potential candidate genes for AN, BN and OCD have been identified through clinical association and neuroimaging studies; however, recent genome-wide association studies of eating disorders (ED) so far have failed to report significant findings. In addition, few, if any, studies have interrogated postmortem brain tissue for evidence of expression quantitative trait loci (eQTLs) associated with candidate genes, which has particular promise as an approach to elucidating molecular mechanisms of association. We therefore selected single-nucleotide polymorphisms (SNPs) based on candidate gene studies for AN, BN and OCD from the literature, and examined the association of these SNPs with gene expression across the lifespan in prefrontal cortex of a nonpsychiatric control cohort (N=268). Several risk-predisposing SNPs were significantly associated with gene expression among control subjects. We then measured gene expression in the prefrontal cortex of cases previously diagnosed with obsessive psychiatric disorders, for example, ED (N=15) and OCD/obsessive-compulsive personality disorder or tics (OCD/OCPD/Tic; N=16), and nonpsychiatric controls (N=102) and identified 6 and 286 genes that were differentially expressed between ED compared with controls and OCD cases compared with controls, respectively (false discovery rate (FDR) <5%). However, none of the clinical risk SNPs were among the eQTLs and none were significantly associated with gene expression within the broad obsessive cohort, suggesting larger sample sizes or other brain regions may be required to identify candidate molecular mechanisms of clinical association in postmortem brain data sets.

BDNF genetic variability modulates psychopathological symptoms in patients with eating disorders

The brain-derived neurotrophic factor (BDNF) gene may influence eating behavior, body weight and cognitive impairments. We aimed to investigate whether BDNF genetic variability may affect anthropometric and psychological parameters in patients with anorexia or bulimia nervosa (AN, BN) and/or modulate the risk for the disorder. A total of 169 unrelated female patients and 312 healthy controls were genotyped for two common BDNF single-nucleotide polymorphisms (SNPs), Val66Met and C-270T, and several selected tag-SNPs. Associated personality characteristics and psychopathological symptoms were assessed by the EDI-2 and SCL-90R inventories, respectively. No single SNP or haplotype played a relevant role in the risk for AN or BN. The rs16917237 TT genotype was significantly associated with increased weight (74.63 ± 16.58 vs. 57.93 ± 13.02) and body mass index (28.94 ± 6.22 vs. 22.23 ± 4.77) in the BN group after correcting for multiple testing. Haplotype analyses using a sliding window approach with three adjacent SNPs produced four loci of interest. Locus 3 (rs10835210/rs16917237/C-270T) showed a broad impact on the measured psychopathological symptoms. Haplotypes CGC and CGT in this locus correlated with scores in all three scales of the SCL-90R inventory, both in AN and BN patients. In contrast, the results of the EDI-2 inventory were largely unaffected. These preliminary results suggest that variability in the BDNF gene locus may contribute to anthropometric characteristics and also psychopathological symptoms that are common but not exclusive of ED patients.

Interpreting genetics in the context of eating disorders: evidence of disease, not diversity

How is genetic involvement interpreted for disorders whose medicalisation is contested? Framing psychiatric and behavioural disorders in terms of genetics is expected to make them seem more medical. Yet a genetic aetiology can also be used to frame behaviour as acceptable human variation, rather than a medical problem (for example, sexual orientation). I analyse responses to the idea that there is a genetic component in anorexia and bulimia nervosa (AN or BN) via semi-structured interviews with a sample of 50 women diagnosed with an eating disorder (25 had recovered). All but three volunteered that genetics would medicalise AN or BN by (i) making eating disorders seem more like 'real diseases'; implying that these disorders need (ii) professional treatment or (iii) a biologically based treatment. The results also indicate there are several counter-logics by which genetic framing could support non-medical definitions of AN or BN. I argue that genetic framing reduces perceived individual responsibility, which can support definitions of behaviour as either a reflection of disease (which entails intervention) or a reflection of normal human diversity (which does not). In the context of public scepticism as to the 'reality' of AN or BN, genetic involvement was taken as evidence of disease in ongoing negotiations about the medical and moral status of people with eating disorders.

Evidence for the role of EPHX2 gene variants in anorexia nervosa

Anorexia nervosa (AN) and related eating disorders are complex, multifactorial neuropsychiatric conditions with likely rare and common genetic and environmental determinants. To identify genetic variants associated with AN, we pursued a series of sequencing and genotyping studies focusing on the coding regions and upstream sequence of 152 candidate genes in a total of 1205 AN cases and 1948 controls. We identified individual variant associations in the Estrogen Receptor-ß (ESR2) gene, as well as a set of rare and common variants in the Epoxide Hydrolase 2 (EPHX2) gene, in an initial sequencing study of 261 early-onset severe AN cases and 73 controls (P=0.0004). The association of EPHX2 variants was further delineated in: (1) a pooling-based replication study involving an additional 500 AN patients and 500 controls (replication set P=0.00000016); (2) single-locus studies in a cohort of 386 previously genotyped broadly defined AN cases and 295 female population controls from the Bogalusa Heart Study (BHS) and a cohort of 58 individuals with self-reported eating disturbances and 851 controls (combined smallest single locus P<0.01). As EPHX2 is known to influence cholesterol metabolism, and AN is often associated with elevated cholesterol levels, we also investigated the association of EPHX2 variants and longitudinal body mass index (BMI) and cholesterol in BHS female and male subjects (N=229) and found evidence for a modifying effect of a subset of variants on the relationship between cholesterol and BMI (P<0.01). These findings suggest a novel association of gene variants within EPHX2 to susceptibility to AN and provide a foundation for future study of this important yet poorly understood condition.

Association of the c.385C>A (p.Pro129Thr) polymorphism of the fatty acid amide hydrolase gene with anorexia nervosa in the Japanese population

The functional c.385C>A single-nucleotide polymorphism (SNP) in the fatty acid amide hydrolase (FAAH) gene, one of the major degrading enzymes of endocannabinoids, is reportedly associated with anorexia nervosa (AN). We genotyped the c.385C>A SNP (rs324420) in 762 lifetime AN and 605 control participants in Japan. There were significant differences in the genotype and allele frequencies of c.385C>A between the AN and control groups. The minor 385A allele was less frequent in the AN participants than in the controls (allele-wise, odds ratio = 0.799, 95% confidence interval [CI] 0.653–0.976, P = 0.028). When the cases were subdivided into lifetime restricting subtype AN and AN with a history of binge eating or purging, only the restricting AN group exhibited a significant association (allele-wise, odds ratio = 0.717, 95% CI 0.557–0.922, P = 0.0094). Our results suggest that having the minor 385A allele of the FAAH gene may be protective against AN, especially restricting AN. This finding supports the possible role of the endocannabinoid system in susceptibility to AN.

Functional genetic variants in the vesicular monoamine transporter 1 modulate emotion processing

Emotional behavior is in part heritable and often disrupted in psychopathology. Identification of specific genetic variants that drive this heritability may provide important new insight into molecular and neurobiological mechanisms involved in emotionality. Our results demonstrate that the presynaptic vesicular monoamine transporter 1 (VMAT1) Thr136Ile (rs1390938) polymorphism is functional in vitro, with the Ile allele leading to increased monoamine transport into presynaptic vesicles. Moreover, we show that the Thr136Ile variant predicts differential responses in emotional brain circuits consistent with its effects in vitro. Lastly, deep sequencing of bipolar disorder (BPD) patients and controls identified several rare novel VMAT1 variants. The variant Phe84Ser was only present in individuals with BPD and leads to marked increase monoamine transport in vitro. Taken together, our data show that VMAT1 polymorphisms influence monoamine signaling, the functional response of emotional brain circuits and risk for psychopathology.

Constitutional thinness and anorexia nervosa: a possible misdiagnosis?

Clinical and biological aspects of restrictive anorexia nervosa (R-AN) are well documented. More than 10,000 articles since 1911 and more than 600 in 2013 have addressed R-AN psychiatric, somatic, and biological aspects. Genetic background, ineffectiveness of appetite regulating hormones on refeeding process, bone loss, and place of amenorrhea in the definition are widely discussed and reviewed. Oppositely, constitutional thinness (CT) is an almost unknown entity. Only 32 articles have been published on this topic since 1953. Similar symptoms associating low body mass index, low fat, and bone mass are reported in both CT and R-AN subjects. Conversely, menses are preserved in CT women and almost the entire hormonal profile is normal, except for leptin and PYY. The aim of the present review is to alert the clinician on the confusing clinical presentation of these two situations, a potential source of misdiagnosis, especially since R-AN definition has changed in DSM5.

Eating disorder predisposition is associated with ESRRA and HDAC4 mutations

Anorexia nervosa and bulimia nervosa are common and severe eating disorders (EDs) of unknown etiology. Although genetic factors have been implicated in the psychopathology of EDs, a clear biological pathway has not been delineated. DNA from two large families affected by EDs was collected, and mutations segregating with illness were identified by whole-genome sequencing following linkage mapping or by whole-exome sequencing. In the first family, analysis of twenty members across three generations identified a rare missense mutation in the estrogen-related receptor α (ESRRA) gene that segregated with illness. In the second family, analysis of eight members across four generations identified a missense mutation in the histone deacetylase 4 (HDAC4) gene that segregated with illness. ESRRA and HDAC4 were determined to interact both in vitro in HeLa cells and in vivo in mouse cortex. Transcriptional analysis revealed that HDAC4 potently represses the expression of known ESRRA-induced target genes. Biochemical analysis of candidate mutations revealed that the identified ESRRA mutation decreased its transcriptional activity, while the HDAC4 mutation increased transcriptional repression of ESRRA. Our findings suggest that mutations that result in decreased ESRRA activity increase the risk of developing EDs.

Genome-wide association analysis of eating disorder-related symptoms, behaviors, and personality traits

Eating disorders (EDs) are common, complex psychiatric disorders thought to be caused by both genetic and environmental factors. They share many symptoms, behaviors, and personality traits, which may have overlapping heritability. The aim of the present study is to perform a genome-wide association scan (GWAS) of six ED phenotypes comprising three symptom traits from the Eating Disorders Inventory 2 [Drive for Thinness (DT), Body Dissatisfaction (BD), and Bulimia], Weight Fluctuation symptom, Breakfast Skipping behavior and Childhood Obsessive-Compulsive Personality Disorder trait (CHIRP). Investigated traits were derived from standardized self-report questionnaires completed by the TwinsUK population-based cohort. We tested 283,744 directly typed SNPs across six phenotypes of interest in the TwinsUK discovery dataset and followed-up signals from various strata using a two-stage replication strategy in two independent cohorts of European ancestry. We meta-analyzed a total of 2,698 individuals for DT, 2,680 for BD, 2,789 (821 cases/1,968 controls) for Bulimia, 1,360 (633 cases/727 controls) for Childhood Obsessive-Compulsive Personality Disorder trait, 2,773 (761 cases/2,012 controls) for Breakfast Skipping, and 2,967 (798 cases/2,169 controls) for Weight Fluctuation symptom. In this GWAS analysis of six ED-related phenotypes, we detected association of eight genetic variants with P < 10(-5) . Genetic variants that showed suggestive evidence of association were previously associated with several psychiatric disorders and ED-related phenotypes. Our study indicates that larger-scale collaborative studies will be needed to achieve the necessary power to detect loci underlying ED-related traits.

"Not all my fault": genetics, stigma, and personal responsibility for women with eating disorders

Medical researchers and clinicians increasingly understand and present eating disorders (anorexia and bulimia nervosa) as biologically-based psychiatric disorders, with genetic risk factors established by high heritability estimates in twin studies. But there has been no research on interpretation of genetic involvement by people with eating disorders, who may hold other views. Their interpretations are particularly important given the frequent presumption that biogenetic framing will reduce stigma, and recent findings that it exacerbates stigma for other mental illnesses. To identify implications of genetic framing in eating disorders, I conducted semi-structured interviews with 50 US women with a history of eating disorders (half recovered, half in treatment; interviewed 2008-9 in the USA). Interviews introduced the topic of genetics, but not stigma per se. Analysis followed the general principles of grounded theory to identify perceived implications of genetic involvement; those relevant to stigma are reported here. Most anticipated that genetic reframing would help reduce stigma from personal responsibility (i.e., blame and guilt for eating disorder as ongoing choice). A third articulated ways it could add stigma, including novel forms of stigma related to genetic-essentialist effacing of social factors. Despite welcoming reductions in blame and guilt, half also worried genetic framing could hamper recovery, by encouraging fatalistic self-fulfilling prophecies and genetic excuses. This study is the first to elicit perceptions of genetic involvement by those with eating disorders, and contributes to an emerging literature on perceptions of psychiatric genetics by people with mental illness.

Examining associations between disordered eating and serotonin transporter gene polymorphisms

OBJECTIVE

The serotonin system has been implicated in mood and appetite regulation, and the serotonin transporter gene (SLC6A4) is a commonly studied candidate gene for eating pathology. However, most studies have focused on a single polymorphism (5-HTTLPR) in SLC6A4; little research has utilized multiple single nucleotide polymorphisms (SNPs) to investigate associations between SLC6A4 and eating pathology more comprehensively.

METHOD

Family-based association tests were conducted for seven polymorphisms in or near SLC6A4, using families from the Colorado Center for Antisocial Drug Dependence. Data were available for 135 families, with phenotypic data available for female twins and female nontwin siblings. Seven items assessed two disordered eating characteristics: weight and shape concerns and behaviors (WSCB) and binge eating (BE).

RESULTS

No significant associations were found between any genetic variant and the two disordered eating characteristics.

DISCUSSION

This study suggests that utilizing polymorphisms in and near SLC6A4, including 5-HTTLPR, may not be useful in identifying genetic risk factors for disordered eating.

Genome-wide association study of comorbid depressive syndrome and alcohol dependence

OBJECTIVE

Depression and alcohol dependence (AD) are common psychiatric disorders that often co-occur. Both disorders are genetically influenced, with heritability estimates in the range of 35-60%. In addition, evidence from twin studies suggests that AD and depression are genetically correlated. Herein, we report results from a genome-wide association study of a comorbid phenotype, in which cases meet the Diagnostic and Statistical Manual of Mental Disorders-IV symptom threshold for major depressive symptomatology and the Diagnostic and Statistical Manual of Mental Disorders-IV criteria for AD.

METHODS

Samples (N=467 cases and N=407 controls) were of European-American descent and were genotyped using the Illumina Human 1M BeadChip array.

RESULTS

Although no single-nucleotide polymorphism (SNP) meets genome-wide significance criteria, we identified 10 markers with P values less than 1 × 10(-5), seven of which are located in known genes, which have not been previously implicated in either disorder. Genes harboring SNPs yielding P values less than 1 × 10(-5) are functionally enriched for a number of gene ontology categories, notably several related to glutamatergic function. Investigation of expression localization using online resources suggests that these genes are expressed across a variety of tissues, including behaviorally relevant brain regions. Genes that have been previously associated with depression, AD, or other addiction-related phenotypes - such as CDH13, CSMD2, GRID1, and HTR1B - were implicated by nominally significant SNPs. Finally, the degree of overlap of significant SNPs between a comorbid phenotype and an AD-only phenotype is modest.

CONCLUSION

These results underscore the complex genomic influences on psychiatric phenotypes and suggest that a comorbid phenotype is partially influenced by genetic variants that do not affect AD alone.

Association analysis of 94 candidate genes and schizophrenia-related endophenotypes

While it is clear that schizophrenia is highly heritable, the genetic basis of this heritability is complex. Human genetic, brain imaging, and model organism studies have met with only modest gains. A complementary research tactic is to evaluate the genetic substrates of quantitative endophenotypes with demonstrated deficits in schizophrenia patients. We used an Illumina custom 1,536-SNP array to interrogate 94 functionally relevant candidate genes for schizophrenia and evaluate association with both the qualitative diagnosis of schizophrenia and quantitative endophenotypes for schizophrenia. Subjects included 219 schizophrenia patients and normal comparison subjects of European ancestry and 76 schizophrenia patients and normal comparison subjects of African ancestry, all ascertained by the UCSD Schizophrenia Research Program. Six neurophysiological and neurocognitive endophenotype test paradigms were assessed: prepulse inhibition (PPI), P50 suppression, the antisaccade oculomotor task, the Letter-Number Span Test, the California Verbal Learning Test-II, and the Wisconsin Card Sorting Test-64 Card Version. These endophenotype test paradigms yielded six primary endophenotypes with prior evidence of heritability and demonstrated schizophrenia-related impairments, as well as eight secondary measures investigated as candidate endophenotypes. Schizophrenia patients showed significant deficits on ten of the endophenotypic measures, replicating prior studies and facilitating genetic analyses of these phenotypes. A total of 38 genes were found to be associated with at least one endophenotypic measure or schizophrenia with an empirical p-value<0.01. Many of these genes have been shown to interact on a molecular level, and eleven genes displayed evidence for pleiotropy, revealing associations with three or more endophenotypic measures. Among these genes were ERBB4 and NRG1, providing further support for a role of these genes in schizophrenia susceptibility. The observation of extensive pleiotropy for some genes and singular associations for others in our data may suggest both converging and independent genetic (and neural) pathways mediating schizophrenia risk and pathogenesis.

Genetic association of recovery from eating disorders: the role of GABA receptor SNPs

Follow-up studies of eating disorders (EDs) suggest outcomes ranging from recovery to chronic illness or death, but predictors of outcome have not been consistently identified. We tested 5151 single-nucleotide polymorphisms (SNPs) in approximately 350 candidate genes for association with recovery from ED in 1878 women. Initial analyses focused on a strictly defined discovery cohort of women who were over age 25 years, carried a lifetime diagnosis of an ED, and for whom data were available regarding the presence (n=361 ongoing symptoms in the past year, ie, 'ill') or absence (n=115 no symptoms in the past year, ie, 'recovered') of ED symptoms. An intronic SNP (rs17536211) in GABRG1 showed the strongest statistical evidence of association (p=4.63 × 10(-6), false discovery rate (FDR)=0.021, odds ratio (OR)=0.46). We replicated these findings in a more liberally defined cohort of women age 25 years or younger (n=464 ill, n=107 recovered; p=0.0336, OR=0.68; combined sample p=4.57 × 10(-6), FDR=0.0049, OR=0.55). Enrichment analyses revealed that GABA (γ-aminobutyric acid) SNPs were over-represented among SNPs associated at p<0.05 in both the discovery (Z=3.64, p=0.0003) and combined cohorts (Z=2.07, p=0.0388). In follow-up phenomic association analyses with a third independent cohort (n=154 ED cases, n=677 controls), rs17536211 was associated with trait anxiety (p=0.049), suggesting a possible mechanism through which this variant may influence ED outcome. These findings could provide new insights into the development of more effective interventions for the most treatment-resistant patients.

Association of candidate genes with phenotypic traits relevant to anorexia nervosa

This analysis is a follow-up to an earlier investigation of 182 genes selected as likely candidate genetic variations conferring susceptibility to anorexia nervosa (AN). As those initial case-control results revealed no statistically significant differences in single nucleotide polymorphisms, herein, we investigate alternative phenotypes associated with AN. In 1762 females, using regression analyses, we examined the following: (i) lowest illness-related attained body mass index; (ii) age at menarche; (iii) drive for thinness; (iv) body dissatisfaction; (v) trait anxiety; (vi) concern over mistakes; and (vii) the anticipatory worry and pessimism versus uninhibited optimism subscale of the harm avoidance scale. After controlling for multiple comparisons, no statistically significant results emerged. Although results must be viewed in the context of limitations of statistical power, the approach illustrates a means of potentially identifying genetic variants conferring susceptibility to AN because less complex phenotypes associated with AN are more proximal to the genotype and may be influenced by fewer genes.

Absence of association between specific common variants of the obesity-related FTO gene and psychological and behavioral eating disorder phenotypes

Extensive population-based genome-wide association studies have identified an association between the FTO gene and BMI; however, the mechanism of action is still unknown. To determine whether FTO may influence weight regulation through psychological and behavioral factors, seven single-nucleotide polymorphisms (SNPs) of the FTO gene were genotyped in 1,085 individuals with anorexia nervosa (AN) and 677 healthy weight controls from the international Price Foundation Genetic Studies of Eating Disorders. Each SNP was tested in association with eating disorder phenotypes and measures that have previously been associated with eating behavior pathology: trait anxiety, harm-avoidance, novelty seeking, impulsivity, obsessionality, compulsivity, and concern over mistakes. After appropriate correction for multiple comparisons, no significant associations between individual FTO gene SNPs and eating disorder phenotypes or related eating behavior pathology were identified in cases or controls. Thus, this study found no evidence that FTO gene variants associated with weight regulation in the general population are associated with eating disorder phenotypes in AN participants or matched controls.