The Phenome-wide Consequences of Anorexia Nervosa Genes

Jessica S. Johnson, Alanna C. Cote, Amanda Dobbyn, Laura G. Sloofman, Jiayi Xu, Liam Cotter, Alexander W. Charney, Eating Disorders Working Group of the Psychiatric Genomics Consortium, Jennifer Jordan, Martin Kennedy, Mikael Landén, Sarah L Maguire, Nicholas G Martin, Preben Bo Mortensen, Cynthia M. Bulik, Laura M. Huckins
doi: https://doi.org/10.1101/2021.02.12.21250941

ABSTRACT

Anorexia nervosa (AN) is a psychiatric disorder with complex etiology, with a significant portion of disease risk imparted by genetics. Traditional GWAS studies produce principal evidence for the association of genetic variants with disease, and provide a jumping-off point for downstream functional analyses. Transcriptomic imputation (TI) allows for the translation of SNPs into regulatory mechanisms, which can then be used to assess the functional outcome of genetically regulated gene expression (GReX) in a more broad setting through the use of phenome-wide association studies (PheWAS) in large and diverse clinical biobank populations with electronic health record (EHR) phenotypes. Here, we applied TI using S-PrediXcan to translate the most recent PGC-ED AN GWAS findings into AN-GReX. For significant genes, we imputed AN-GReX in the Mount Sinai BioMe™ Biobank and performed PheWAS on over 2000 clinical outcomes to test the clinical consequences of aberrant expression of these genes. We performed a secondary analysis to assess the impact of BMI on AN-GReX clinical associations.

Our S-PrediXcan analysis identified 47 genes associated with AN, including what is, to our knowledge, the first genetic association of AN with the Major Histocompatibility Complex (MHC). AN-GReX was associated with autoimmune, anthropometric, metabolic, psychiatric and gastrointestinal diagnoses in our biobank cohort, as well as measures of anthropometry, substance use, and pain score. Our analyses reveal that AN-GReX associations with measures of weight and substance use are modified by BMI, and indicate potential avenues of functional mechanism to investigate further.

INTRODUCTION

Anorexia nervosa (AN) is a severe eating disorder characterized by extreme low body weight and fear of gaining weight, and is often accompanied by compensatory behaviors to lose weight such as dietary restriction, purging, and higher than average physical activity(1). With a lifetime prevalence of 0.9-4%, AN has high mortality, with few effective treatments and an increased risk of suicide(2–4). Clinical and epidemiological studies have shown higher risk of other psychiatric disorders such as OCD, anxiety, and major depression in individuals with AN(3,5), and higher risk of substance use disorders in AN subtypes (6), as well as higher risk of metabolic and autoimmune disorders such as type 1 diabetes(7,8),. Twin studies have established the heritability of AN to be between 50 to 60%, indicating a considerable contribution of genetic factors to AN disease liability(9,10).

Both metabolic and psychiatric factors contribute to anorexia nervosa. Epidemiological studies indicate high rates of psychiatric comorbidity, particularly with depression, anxiety disorders, and OCD(5), and many of these disorders share substantial symptomatology and environmental risk factors(11). Evidence is mounting that AN may lie at the extreme of a BMI spectrum, with obesity lying at the opposite end, and that AN may be a disorder of “fundamental metabolic dysregulation”(2). In line with this hypothesis, individuals with AN are able to reach and maintain incredibly low body weight, and often revert to this weight even after re-nourishment treatment, possibly indicating a fundamentally low BMI biological settling point(12,13). Additionally, individuals with AN have been shown to have lower concentrations of leptin, a hormone involved in regulation of body weight and other metabolic processes(13–16). Microbiome studies have revealed abnormal microflora in the gut of AN patients(17), and evidence from many studies have pointed to a possible brain-gut etiology of AN(9,18–20).

Genetic studies of AN have provided further evidence of both psychiatric and metabolic etiology. The largest AN GWAS to date (NCases=16992) uncovered eight loci associated with AN risk(9), and determined SNP-based heritability (h2SNP) of 11-17%, similar to other psychiatric disorders, indicating that common variants contribute to polygenic risk of AN. In addition, genetic correlations demonstrate significant overlap with other psychiatric disorders, including schizophrenia(9,21), MDD(9), OCD(9), alcohol use disorder(22) and anxiety disorders(9), as well as physical activity(9,23), indicating shared genetic variation between these traits. Significant negative genetic correlation between AN and anthropometric and metabolic traits, such as BMI(9,21,23), fat mass(9,21,23), obesity(9), type 2 diabetes(9), insulin resistance(9,21), fasting insulin(9,21) and leptin(9) have also been observed, further indicating metabolic components to AN disease risk. One study of AN polygenic risk scores (PRS) has indicated additional genetic associations of AN risk variants with anthropometric, behavioral, and psychiatric traits(24).

Traditional observational genetic studies provide substantial data for factors contributing to risk of a particular disease, but may not capture the genetic factors that are specifically contributing to subthreshold or prodromal disease states. Electronic Health Records (EHR) contain information about an individual’s health history such as diagnoses, medications, laboratory tests, vital signs, and family medical history, and can be utilized for scalable disease research. The vast amounts of data from the EHR can be queried to provide an understanding of the clinical spectrum of disease and disease progression across the lifetime of the patient. The initiation of biobanks to collect biological specimens from large numbers of patients within a healthcare system allows us to connect massive amounts of phenotype data from the EHR to genetic data. A phenome-wide association study or PheWAS is an effective method of querying the EHR to look for the associations of a trait of interest with the clinical phenome,(25–27). Previous studies have demonstrated the efficacy of PheWAS methods in determining clinical outcomes associated with genomic disease risk variants, and have replicated associations from GWASs of strictly defined disease traits, including psychiatric disorders(28–30),. Exploring the phenotypic associations with AN genetic architecture has the potential to clarify how some of these GWAS variants functionally contribute to AN disease risk, symptomatology, and clinical presentation.

Transcriptomic imputation (TI) approaches leverage expression quantitative trait loci (eQTLs) derived from large, well-curated gene expression datasets (for example, the Genotype-Tissue Expression project (GTEx)(31), CommonMind Consortium (CMC)(32), and Depression Genes and Networks (DGN)(33)) in order to derive gene expression predictor models. These models may be applied to predict genetically regulated gene expression (GReX) in large genotyped cohorts, without the need to collect tissue samples or RNA(19,34–36). Here, we explored the association between AN-associated GReX and the clinical phenome. We performed transcriptomic imputation using S-PrediXcan on the most recent PGC-ED AN GWAS to first find GReX associated with AN, and then tested for clinical associations of these genes with structured EHR-encoded phenotypes using PheWAS. Given the relationship of AN with BMI, we further investigated the effect of BMI on the GReX-phenotype associations by stratifying biobank individuals based on sex- and ancestry-adjusted categories of low, mid, and high BMI. Understanding the clinical consequences of aberrant gene expression across the phenome may clarify the biological mechanisms of clinically relevant AN GWAS risk variants.

METHODS

AN GENELIST and TRANSCRIPTOMIC IMPUTATION

We performed transcriptomic imputation using S-PrediXcan(35) on the largest available summary statistics of Anorexia Nervosa (AN) (16992 AN cases and 55525 controls) downloaded from the Psychiatric Genomics Consortium (PGC) website(9). We tested for association of genetically-regulated gene expression (GReX) using available GTEx, CMC and DGN predictor models(34–36) for a total of 50 different tissues with AN case-control status. We established two thresholds for significance; first, correcting for all genes tested within each tissue (ptissue, values shown in Table S1, and second, correcting for all tissues and genes tested (pbonferroni=3.75 × 10−8).

BioMe™

The BioMe™ Biobank at the Icahn School of Medicine at Mount Sinai includes genotype and EHR data from 31704 individuals recruited from the general patient population of the Mount Sinai healthcare system. Individuals were genotyped on the Illumina Global Screening Array (GSA), an array designed for multi-ethnic populations. The BioMe™ population is very diverse, with the proportion of major ancestral population groups representative of the diversity of New York City and Mount Sinai’s patient population: Hispanic (HA, 35%), European (EA, 34%), African (AA, 25%), East Asian (ESA, 3%), South Asian (SAS, 2%), and Native American (NA, 0.3%) (Full BioMe™ demographics are shown in Table 1). Quality control and imputation of the genotyping data for BioMe™ is described in elsewhere(28).

Ancestry for BioMe™ was initially assigned from self-reported values into six ancestral groupings: African (AA), European (EA), Hispanic (HA), East Asian (ESA), Native American (NA), and O (Other). We plotted the first two principal components from the genotyping to confirm ancestral grouping. For individuals not clearly falling in an ancestral group, and for those designated “Other”, we merged genotype information with phase 3 of 1000Genomes and reran the principal components analysis (PCA) using PLINK(37,38) in order to assign individuals to 1000Genomes super-populations. From this we were able to assign individuals to one of six ancestral groups: African (AA), European (EA), East Asian (ESA), Hispanic (HA), Native American (NA), and South Asian (SAS). A total of 119 samples that could not be grouped were removed, leaving 31585 individuals for analysis.

PheWAS

We calculated GReX for all ptissue significant genes across GTEx(31), DGN(33), and CMC(32,39) tissues (N=50) in BioMe™ individuals using PrediXcan-2(34,35) software, and performed a phenome-wide association (PheWAS) analysis. Logistic regression between the calculated AN-GReX and phecodes was performed using the PheWAS R package(40), adjusting for age, biological sex, and the first five genotype-derived principal components. PheWAS were run per BioMe™ ancestry cohort of (NAncestry=6) and results meta-analyzed using an inverse-variance based approach in METAL(41). Due to the nature of EHR data, sample sizes can vary greatly depending on the phenotype, which can, in turn, affect our ability to detect associations. We believe our study is sufficiently powered given results from a previous study(42) that used simulations of MAF, penetrance and sample number parameters to estimate power for pheWAS studies.

Encounter Diagnosis

Current diagnoses for each individual were recorded at each visit using the International Classification of Disease (ICD) coding system and recorded as “Encounter Diagnosis”. Phecodes were assigned from Encounter Diagnosis information by grouping ICD-9 and ICD-10 diagnostic codes(43). Individuals with at least two counts of a particular diagnostic code were considered “cases”, those with zero counts were marked as “controls”, and those with only one count were set to missing. After QC, there were a total of 2178 unique Encounter Diagnosis codes and 1093 unique phecodes in the dataset.

For both encounter diagnoses and phecodes, we required at least 10 occurrences of each phenotype within each population group for inclusion in the analysis, and an overall effective sample size Neff > 100 (Equation 1). Embedded Image Due to the high correlation between Encounter Diagnosis and phecode files, we combined all results from both PheWAS and performed an FDR correction in R to determine overall significant associations with diagnostic phenotypes. Significant associations were set at an FDR-corrected p threshold of 0.05.

BioMe-R Structured Phenotype Database (BRSPD) PheWAS

In addition to diagnosis code data, we also had access to other types of EHR-derived data that included information on allergy status, vital signs (weight, height, blood pressure, pulse, pulse oximetry, respirations, and temperature), family history, personal history, medication use, obstetrics and gynecology outcomes (OB/GYN), and other social and behavioral traits. This information derived from the EHR was collected into files under the heading BRSPD or BioMe-R Structured Phenotype Database. Summaries and descriptive statistics of all of the BRSPD phenotypes can be found in Supplementary Methods. For all BRSPD traits, we required at least ten cases in each ancestry group, and a total Neff>100 for inclusion in our analysis.

Phenotype QC and Analysis for Continuous Traits

We plotted the distribution of continuous traits across the cohort, and removed implausible or impossible values (likely due to errors in entry; for example, height listed as 1745 inches). Continuous traits were then analyzed as follows:

  1. Height (cm): mean measure per individual.

  2. Weight (kg): (i) highest recorded weight, (ii) lowest recorded weight, (iii) weight change over time. Weight change over time was calculated by taking the difference between the highest and lowest recorded weights, divided by the number of years between those two measures.

  3. Pain scores were self-reported on a 1 to 10 scale and assessed for 25 different body locations. First, we tested pain score regardless of location, using three measures: (i) highest ever pain score, (ii) mean pain score, (iii) total pain score (sum), to encompass both acute and chronic pain conditions. We then repeated this analysis to each reported pain location, with the same three measures.

  4. Measures of alcohol, tobacco, and illicit drug use were categorized as follows: alcohol-ounces (oz) per day; tobacco- (i) packs per day and (ii) years used; illicit drug use-frequency of use.

  5. Vital signs (blood pressure, pulse, pulse oximetry, respiration, temperature (Celsius)) were analyzed using (i) highest ever value, (ii) lowest ever value, (iii) mean, (iv) variance per individual.

Categorical Traits

Categorical traits were converted to binary TRUE/FALSE values (phenotypes included are shown in Supplemental Information). For allergy and medication phenotypes, we collapsed and combined data on equivalent root components and/or active ingredients (e.g., allergies for the antibiotic Vancomycin were coded as “VANCOMYCIN”, “VANCOMYCIN HCL”, “VANCOMYCIN IN D5W”, “VANCOMYCIN IN 0.9% SODIUM CL”, “VANCOMYCIN IN DEXTROSE ISO-OSM”, and “VANCOMYCIN (BULK)”), removed extraneous symbols and corrected typing errors. A full account of QC steps is documented in the Supplementary Methods.

STRATIFICATION BY BMI

Given the nature of AN symptomology and association with weight phenotypes, we tested whether AN-GReX associations varied based on BMI. We stratified BioMe™ individuals into three main BMI categories, defined based on ancestry- and sex-specific BMI distributions within our data (Figures S1 and S2; Tables S2 and S3). BMI values used for these analyses were measured at intake into BioMe™. We categorized individuals with BMI within the 1st quartile as “Low” BMI; those above the 3rd quartile as “High” BMI, and those falling between the 1st and 3rd quartiles as “Mid” BMI (Figure S1). We then repeated our analysis in these three groups for all significant associations from our full analysis.

Testing for hidden case contamination

It is possible that some of the clinical associations within our study stem from undiagnosed ED- and AN-cases, rather than a direct effect of gene expression on phenomic expression. We term this “diagnostic contamination”. To test whether this effect may drive the associations we observe, we simulated the effects of diagnostic contamination occurring at rate p within our biobank sample, with effect β on GReX.

Full derivations and simulations are shown in our Supplementary Material; briefly, we derive the expression among controls, and cases contaminated at rate p; the difference in expression between the two groups; the expected variance among cases, and pooled across all samples, and the expected statistical significance (T-score, and p-value). Next, we simulated gene expression distributions for (i) 1000 cases and 1000 controls; (ii) 1000 cases and 10000 controls; (iii) 1000 cases and 30000 controls; each for a range of p (0.1%, 0.5%, 1%, 2%, 5%, 10%, 20%, 30%, 40%, 50%) and β (1/10, 1/5, ¼, 1/3, ½, 1, 2,3,4,5,10) values. We repeated our simulations 10000 times at each case-control proportion-p-β combination, and demonstrated that our formulae accurately estimate the desired values (Figure S3).

Next, in order to test whether diagnostic contamination may account for the associations observed within our PheWAS, we calculated the expected impact of diagnostic contamination for two genes (NCKIPSD-Aorta; SEMA3F-Spinal Cord) with the largest effect sizes β observed in our S-PrediXcan analysis, under two scenarios:

  1. Diagnostic contamination occurs among our cases only, assuming the normal population rate for AN: 0.9% among women; 0.3% among men.

  2. Contamination occurs at significantly higher levels than might be expected in the population (1.2%, 3%, 6%), and all of these samples fall into our case category.

For each scenario, we repeated our calculation using 500 and 1000 cases, matched with 1000, 5000, and 30000 controls.

RESULTS

S-PrediXcan identifies gene-tissue associations with AN

We applied S-PrediXcan to PGC-AN GWAS summary statistics, and identified 47 genes across 12 loci, two loci on chromosomes 2 and 3 that reached experiment-wide significance (p<3.75 × 10−8; Figure 1, Table 2, Table S4), and another 10 loci that reached tissue-specific significance (Each tissue type has its own p value threshold of significance; Figure 1, Table 2, Table S4). These 12 loci encompass 218 gene-tissue associations (which may include the same gene associated across multiple tissues). The most significant peak on chromosome 3 includes 29 genes with significant associations with AN; this region of the genome overlaps with a GWAS peak for inflammatory bowel disease (IBD)(44), Crohn’s disease, and ulcerative colitis. Many of the AN-associated genes within this peak have been found to be associated with IBD. Genes included in our AN GReX gene set have been previously associated not only with IBD, but with other psychiatric traits such as schizophrenia (ARNTL, CLIC1), depression and/or depressive symptoms (C3orf62, CCDC36, CLIC1, DAG1, GPX1, MST1R), autism (CLIC1), and substance use (RBM5, RBM6), as well as anthropometric measures such as BMI, waist-hip ratio, and lean body mass (Table S5). Of these genes, CLIC1 lies within the major histocompatibility complex (MHC), a large region on chromosome 6 that includes genes that code for components of the adaptive immune system (Tibial nerve-CLIC1, p=1.47 × 10−6). To our knowledge, this is the first genetic association for AN within the MHC region.

Figure

AN-GREX is associated with anthropometric measures

We tested whether GReX of AN-associated genes was associated with anthropometric phenotypes in our BioMe™ biobank (Table 3). In the overall cohort, upregulation of MHC-gene CLIC1 was associated with lowest recorded weight (Spleen-CLIC1, p=6.63 × 10−4; Breast Mammary tissue-CLIC1, p=4.18 × 10−4) and mean height (Subcutaneous Adipose-CLIC1, p=1.39 × 10−5; Table S6). Downregulation of multiple genes, ARIH2, NCKIPSD, NDUFAF3, NICN1, P4HTM, RNF123, SLC26A10, was associated with decreased measures of weight change over time in the overall cohort (p<9.86 × 10−4)(Table 3, Figure S4). Upregulation of WDR6 was associated with increased measures of weight change over time (p<9.81 × 10−4)(Table 3, Figure S4). CTNNB1, GPR75 and LINC00324 GReX was associated with mean height in BioMe™ (p<6.51 × 10−4)(Table S6).

Patient BMI mediates AN-GREX-PheWAS associations

Given the relationship of AN with BMI, we performed a secondary analysis to assess the effect of BMI on the association between AN-GReX and anthropometric phenotypes. Specifically we wanted to look at whether our AN-GReX-phenotype associations in individuals of High and Low BMI were still present when BMI was within a normal range. We split the BioMe™ cohort into three groups: High, Mid, and Low BMI (adjusted for ancestry and sex; see Methods and Supplementary Information for more details), and re-ran our PheWAS within these groups. Indeed, our AN-GReX anthropometry associations (with highest recorded weight, lowest recorded weight, and weight change over time) varied according to the BMI group tested (Figure 2, Table 3). The association between CLIC1-GReX and highest and lowest recorded weight measures was strongest in individuals of High BMI compared with the other BMI groups (Figure 2, Table 3). Likewise, for multiple genes, associations between AN-GReX and weight change over time differed markedly in the High BMI group (Table 3, Figure S5).

Figure

AN-GREX is associated with autoimmune and autoinflammatory diseases

At a false discovery rate (FDR) p<0.05, we found 17 FDR-significant gene-tissue associations with four Phecode- and Encounter Diagnosis phenotypes: type I diabetes, celiac disease, peptic ulcer, and unspecified immunodeficiency (Figure 3). Upregulation of CLIC1 in spleen was associated with two autoimmune/autoinflammatory phenotypes in the overall BioMe™ cohort: celiac disease and type I diabetes (p=1.30 × 10−11, p=5.08 × 10−20, respectively) (Tables S7 and S8). Downregulation of SLC26A6 was associated with the phenotype “peptic ulcer, site unspecified” in the overall cohort (Liver-SLC26A6, p=8.38 × 10−8; Skeletal muscle-SLC26A6, p=1.71 × 10−7; Thyroid-SLC26A6, p=7.16 × 10−8). Downregulation of PFKFB4 was associated with ICD10 code D84.9, “immunodeficiency, unspecified” (Skin not sun exposed-PFKFB4, p=1.34 × 10−7). We stratified our cohort by BMI (see Methods) to assess whether BMI was modifying the association of AN-GReX with Encounter Diagnosis and phecodes. We saw significant associations of AN-GReX phecodes within our Mid BMI group; however, when plotted across the BMI groups, there was no notable difference in the magnitude and direction of effect of AN-GReX on either Encounter Diagnosis or phecodes across the BMI groups (Figures S6 and S7). Within individuals with Mid BMI, upregulation of CLIC1 was similarly associated with type 1 diabetes (Mid-Subcutaneous adipose, Mid-Breast Mammary tissue, Mid-Spleen-CLIC1, p<1.96 × 10−8).

Figure

AN-GREX is associated with tobacco, alcohol, and illicit drug use

Multiple AN genes were associated with both categorical and continuous measures of tobacco use (Table S9, Figure S8). APEH, CCDC36, DALRD3, GPX1, LINC00324, NICN1, and P4HTM were associated with years used/smoked of tobacco in the overall BioMe™ cohort (p<4.97 × 10−4, Table S9) and upregulation of MGMT was associated with being a smoker (Heart Atrial Appendage-MGMT, p=8.87 × 10− 5; Putamen basal ganglia-MGMT, p=1.09 × 10−4). Again, we note mediating effects of patient BMI category on the association between AN-GReX and smoking phenotypes. Among the Mid and High BMI groups, different genes were associated with both continuous and categorical measures of tobacco use (Table S9, Figure S8), including years of tobacco use, tobacco cigarette packs per day, and smoking every day. In addition to tobacco use measures, we see an association between AN genes and tobacco use related ICD codes. Downregulation of LINC00324 in Low BMI individuals was associated with ICD10 diagnosis code F17.210 “Nicotine dependence, cigarettes, uncomplicated” (Low-Cells EBV-tranformed-LINC00324, p=3.03 × 10−5).

Downregulation of MGMT was associated with decreased alcohol use (ounces per week) in individuals in the Mid-BMI group (p<9.75 × 10−4; Figure S8). In individuals of Low BMI, downregulation of GPR75 and GPX1, and upregulation of NICN1 were associated with decreased and increased continuous measures of illicit drug use frequency respectively (Liver-GPR75, p=3.34 × 10−4; Nucleus accumbens-GPX1, p=2.04 × 10−4; Testis-NICN1, p=3.34 × 10−4; Figure S8).

AN-GREX is associated with measures of pain score and pain location

We tested for association between AN-GReX and three sets of pain score measures: (1) pain location, (2) pain score overall, (3) pain score by body location. AN-GReX was associated with multiple measures of pain score measurement and location across all BMI groups (Table S10). Most significantly, upregulation of TNFSF12 was associated with total knee pain score (sum) in individuals with High BMI (High-Brain Substantia nigra-TNFSF12, p=1.50 × 10−3211). Downregulation of TNFSF12 was associated with total abdomen pain score (sum) in individuals with Low BMI (Low-Brain Substantia nigra-TNFSF12, p=3.93 × 10−58) (Table S10).

When looking at pain location irrespective of pain score, upregulation of TUSC2 in the overall cohort is associated with presence of knee pain (Pancreas-TUSC2, p=7.62 × 10−6) and downregulation of CTNNB1 in Mid BMI individuals is associated with the presence of left leg pain (Mid-Stomach-CTNNB1, p=1.93 × 10−5) (Figure 4, Table S10). Repeating this analysis according to BMI category reveals a very large number of associations specific to the Low BMI category (Figure 4, Table S10); in particular, we note 129 gene-tissue associations with presence of foot pain, nine of which pass phenotype-wide significance, which may indicate a propensity to excess exercise, or exercise-related injuries among these individuals.

Figure

Upregulation of CLIC1 is associated with glucagon medication

AN-GReX PheWAS with medication phenotypes was associated with glucagon, a hormone used to treat severe hypoglycemia, in the overall BioMe™ cohort. Upregulation of CLIC1 was associated with glucagon in the overall cohort (Spleen-CLIC1, p=4.20 × 10−9) and also in individuals with Mid BMI (Mid-Spleen-CLIC1, p=2.09 × 10−6) (Figure S9, Table S11).

Identifying potential case contamination effects

Our BioMe™ patients have not all been explicitly assessed for eating disorders (ED), and information regarding ED diagnoses and assessments earlier in life may be omitted from the records due to incomplete clinical history assessments. Therefore, it is possible that diagnostic contamination within some of our sample is responsible for the associations observed within our data. In order to test this, we performed a simple experiment.

For this simple experiment, we assume that diagnostic contamination occurs among true cases only, at prevalence p. In this case, we expect gene expression among those individuals to be up-or down-regulated by effect size beta; consequently, we expect a gene expression difference between cases and controls proportional to the prevalence of the contamination within the sample (p) and the effect size (beta); (Supplementary information). We tested whether such contamination may drive the associations observed in our study, assuming the following different contamination scenarios (Table S12), for two genes with large S-PrediXcan effect sizes: (i) that contamination occurs among cases, assuming the highest common estimate for AN prevalence (0.9% among women, 0.3% among men); (ii) that contamination occurs within our biobank, assuming the highest common estimate for AN prevalence (0.6%, 190 cases), and that all of these samples fall into our ‘case’ category.

Our model indicates that diagnostic contamination among cases is unlikely to result in significant PheWAS associations in our data. Contamination occurring at population ED prevalence (0.6%) did not result in significant associations for any of the case-control scenarios in our model (Table S12); nor did contamination at 1.2% or 3%. Assuming 6% contamination resulted in potentially significant associations (estimated p<2.8×10−8; see Methods and Supplementary Methods for more detail).

DISCUSSION

AN-associated genes are enriched for associations with psychiatric, autoimmune, and anthropometric phenotypes

Our S-PrediXcan gene results illustrate the contribution of genes associated with metabolic, anthropometric, autoimmune, and psychiatric phenotypes to AN, and 53% of our genes (25/47) overlap with those from a previous prediXcan analysis of AN cohorts(9). 24/47 AN genes (51%) had variants that were previously associated with some autoimmune disorder in GWAS studies (Table S5). The experiment-wide significant loci on chromosome 3 overlaps with a known GWAS peak for inflammatory bowel disease(44) (IBD), and many of the genes within that peak that we found associated with AN were also associated with disorders under the IBD umbrella such as Crohn’s and ulcerative colitis (Table S5), demonstrating further shared genetic architecture of AN with autoimmune and autoinflammatory disease(7,8,45,46). Earlier AN GWAS studies have found AN risk loci that overlap with loci implicated in autoimmune disease(47), but the association of CLIC1 gene expression in tibial nerve that we have found is the first known association of AN with the MHC locus, a region that has been associated with many other psychiatric disorders(48,49) as well as autoimmune disorders. In particular, CLIC1 has previously been identified as associated with schizophrenia(50), autism(50), major depressive disorder(51), post-traumatic stress disorder(52), neuroticism(53) and depressive phenotypes(53). Importantly, CLIC1 variants have also been associated with complement component C4 and C3 protein levels in the blood(54). C4 is a gene with two isotypes, C4A and C4B, that, through various structural allelic combinations, give rise to complement component 4, which, along with complement component 3 (C3), is part of the complement cascade of proteins(55). Previous research has shown that the complement cascade is involved in not only immunological functions of pathogen clearance (e.g., tagging pathogens for destruction), but also in synaptic pruning and neuronal connectivity(56). The C4A isotype has been associated with schizophrenia risk(55), and both C4A and C4B genes had subthreshold associations with AN in our dataset.

CLIC1 (Chloride Intracellular 1) encodes a chloride ion channel protein involved in many necessary cellular functions including the regulation of cell membrane potential, and the proliferation and differentiation of cells(57), including a role in axonal outgrowth of neurons51. In addition to neurons, CLIC1 has been found to be expressed in many different tissue types(59), and in particular, CLIC1 protein has been found to accumulate in peripheral blood mononuclear cells (PBMCs) when the central nervous system is in a state of chronic inflammation(58), and increased CLIC1 expression has been shown to occur in microglia stimulated by the presence of β-amyloid in vitro(60,61). Therefore, CLIC1 may potentially act as a biomarker for neurodegenerative disorders such as Alzheimer’s disease, and downregulation of CLIC1 has been shown to have a neuroprotective effect in the presence of β-amyloid proteins(58,60–62). Variants in CLIC1 may contribute to risk of multiple psychiatric disorders and behaviors, such as schizophrenia, MDD and PTSD, as well as many autoimmune and immune-related phenotypes such as asthma and systemic lupus erythematosus(50–53),(63) (Table S5). Many studies have posed the role of inflammation in the development of psychiatric disorders, and in the etiology of AN(7,8,45,46,64).

Multiple genes (8/47 or 17%) in our results were associated with regular gym attendance (Table S5), or measures of physical activity. Excessive and compulsive exercise is a behavior often seen in individuals with AN(65,66), and evidence is suggestive that hyperactivity increases risk of chronic AN(67,68). Likewise, genetic studies have shown a strong genetic correlation between AN and physical activity (rg=0.17)(9,23,69). These genes may reflect a genetic liability to compulsive behaviors of physical activity. Our pain location results point to possible symptoms of excessive activity, with the strong association of AN-GReX with foot pain in Low BMI individuals. Given the known association of AN disease with low bone mineral density (BMD) and propensity for bone fracture(70), our results may reflect the result of genes associated with compulsive exercise and BMD contributing to increased osteoarticular pain. We also see multiple genes that have been previously associated with decreased sleep and insomnia phenotypes (including a core circadian clock gene ARNTL), which have previously been implicated in a range of psychiatric disorders(71–82), including eating disorders(83), as well as in satiety and hunger(84,85).

In a patient population, expression of AN genes is associated with psychiatric, metabolic, and autoimmune phenotypes

Unlike GWAS, which include carefully constructed case-control cohorts, PheWAS encompass all individuals within a given healthcare system, including patients with subthreshold or partial presentations of a disorder, or individuals with commonly comorbid or co-diagnosed conditions. Importantly, because individuals enrolled in the biobank are not ascertained for any specific disorder, they likely represent a comprehensive clinical picture of the comorbidities and symptomatology associated with AN gene expression. Examining the consequences of aberrant predicted gene expression among these patients (i.e., testing for GReX-associations) may reveal clinical and biological consequences of these genes; for example, studying whether AN-associated genes have anthropometric and metabolic consequences among adults with no evidence for previous AN-or ED-diagnoses may disentangle whether certain endophenotypes present as a cause, or consequence of AN. For example, it has been postulated that behaviors of food avoidance and restriction may arise due to gastrointestinal complaints and distress that provoke these behaviors and precede development of AN(86). Likewise, autoimmune disorders of the gastrointestinal tract, such as celiac disease and Crohn’s disease, show a bidirectional relationship with AN, with previous diagnosis of a GI-associated autoimmune disorders increasing the risk of AN and vice versa(8). Our pheWAS results of AN-GReX associations with gastrointestinal symptoms such as abdominal pain, ascites, and peptic ulcer, as well as GI-related autoimmune disorders like celiac disease, suggest AN-GReX may contribute directly to these diseases and phenotypes, and that food aversive behaviors and gastric distress may be genetically regulated in these individuals, rather than occurring as a consequence of AN. While these associations could be the consequence of undiagnosed AN individuals in our biobank, they more likely reflect real biological associations of expression of those particular genes with the phenotype.

Our results further confirm the contribution of metabolic factors to AN etiology as we see a very robust association of AN-GReX with type 1 diabetes, and our top findings for AN-GReX with medication prescriptions include the hyperglycemic hormone glucagon and various forms of insulin. The association of AN with metabolic traits and abnormalities has been fairly well established. The Anorexia Nervosa Genetics Initiative reported significant negative genetic correlation between AN and a number of metabolic traits including insulin resistance (rg=-0.29), fasting insulin (rg=-0.24), leptin (rg=-0.26) and type 2 diabetes (rg=-0.22), and a strong positive genetic correlation with HDL cholesterol (rg=0.21)(9). Additional studies have shown associations of AN with fasting insulin(87) and insulin sensitivity(88). Notably, our results point to a similar role of aberrant glycemic regulation in the etiology of AN. Future analyses including EHR-derived lab results (LabWAS) studies may further elucidate AN genes associated with abnormal metabolic regulation and clinical features.

BMI modification of AN-GReX on clinical outcomes

We stratified our pheWAS analyses by BMI in order to observe whether the effects of predicted AN gene expression on clinical outcomes were modified by BMI status, and whether this modification had a significant impact on clinical phenotypes. Specifically, are the AN-GReX-phenotype associations that we see in individuals of High and Low BMI still present when BMI is within the normal “healthy” range? Understanding how BMI is influencing the association of AN-GReX with the clinical phenome might give us further insight into the biological pathways leading to those outcomes and the conferral of AN disease risk. Additionally, we want to understand how BMI could be interacting with AN genetic risk to influence clinical outcomes. We know from our S-PrediXcan analyses and previous GWAS that many genes found to be associated with AN have variants that have been associated with BMI and other anthropometric traits (Table S5). 20/47 AN-genes (42.6%) had evidence from previous GWAS of association with anthropometric measures such as BMI, body fat distribution, waist-to-hip ratio, fat-free mass, obesity, and hip circumference. BMI polygenic risk scores (PRS) have been shown to be associated with disordered eating, with high BMI PRS negatively associated with thin ideal internalization and external eating, where individuals tend to eat in response to external cues of the sight or smell of food(89,90).

We postulate that the modification of AN-GReX-weight associations by BMI could be due to actual biological differences in the regulation of AN gene expression and downstream processes in the presence of different measures of BMI. CLIC1 is involved in pathways of inflammation, and individuals with a higher BMI have been shown to have increased systemic inflammation(91), especially in adipose tissues, so it is possible that the modification effect by BMI on the effect of CLIC1-GReX on weight measures is happening via a tissue-specific inflammatory process. If we view AN and obesity as metabolic disorders on opposite ends of the spectrum, then perhaps the interaction of BMI with risk for either disorder may lead to differences in metabolic regulation resulting in higher weight phenotypes for those with increased obesity risk, and lower weight phenotypes in those with higher AN risk.

Although the hypothesis-free, phenome-wide design of our study allows for powerful detection of clinical and biological consequences of AN risk genes, the same design also bears some notable limitations and caveats. One key caveat to our results is the lack of diagnostic detail and insights available for the patients in our study. EHR analyses leverage large, existing data sets to rapidly amass cohorts for analysis, and to yield insights into whole phenome consequences of genotype and GReX associations; however, the scale and scope of these studies precludes deep phenotyping or performing clinical interviews. This lack of diagnostic precision may arise from a number of factors.

First, we make use of ICD codes within the medical record in order to infer diagnoses; since these are primarily used by clinicians for billing purposes, they likely provide an imperfect proxy for true disease state. In order to mitigate spurious results stemming from imperfectly assigned codes, we established criteria for inclusion of diagnoses in our analyses: First, diagnosis codes with fewer than 10 counts in the biobank were removed. Second, we only included diagnoses in our pheWAS with at least 10 counts of a code in the ancestral population (i.e., ancestral groups with fewer than 10 case counts were not included in the analysis for that diagnosis), and with a total effective sample size after meta-analysis that was greater than 100. Lastly, within each pheWAS, individuals had to have two or more counts of a diagnosis code to count as a “case”, whereas individuals with only one count of a diagnosis code were set to “missing” and not included in the analysis.

Second, it is possible that our patients regularly receive treatment at multiple different hospital systems; as such, we may be capturing only partial data for each of our patients. In order to mitigate this, we restricted our analysis to patients with multiple data points within our EHR. Future analyses that seek to harmonize or meta-analyze patient data across EHR (e.g., NYC-CDRN, PsycheMERGE) are ongoing, and will be vital to disentangling this effect further.

Importantly, our patients have not all been explicitly assessed for EDs, and information regarding ED diagnoses and assessments earlier in life may be omitted from the records due to incomplete clinical history assessments. Therefore, it is possible that case-contamination within some of our sample is responsible for the associations observed within our data. In order to address this, we performed a simple experiment to simulate the effect sizes expected if undiagnosed AN contamination drives our result. We estimated the expected association statistics observed due to two possible levels of contamination among PheWAS cases; i, that 0.06% of PheWAS-trait cases are missing a true AN diagnosis-i.e., the proportion of individuals we expect to receive an AN diagnosis based on standard prevalence estimates; ii, that 1.2%, 3% and 6% of individuals in BioMe are missing an AN diagnosis (i.e., double, five times, and ten times the expected population rate), and that all of these individuals fall within our PheWAS trait case group, and applied this model for two genes with large S-PrediXcan effect sizes (NCKIPSD-Aorta; SEMA3F-Spinal Cord). Our model indicates that diagnostic contamination at 0.06%, 1.2% or 3% is insufficient to account for the effect sizes observed within our PheWAS analysis. Moreover, since these genes are selected due to their large S-PrediXcan effect sizes, we expect that the contamination effects observed for these two genes will be greater than any others in our study; as such, we do not believe that our findings are attributable to hidden case contamination. Contamination at 6% resulted in significant differences; however, since this estimate represents a >five-fold increase in expected numbers of AN within a population, we do not expect that this level of contamination occurs in our sample.

Our results demonstrate that there are real clinical consequences to differences in AN gene expression. Characterization of the phenotypic consequences of AN gene expression in a clinical setting can give us more insight into the biological mechanisms underlying AN and, consequently, how to diagnose and treat the disorder. By understanding the associations of AN gene expression with symptomatology, prodromal or subthreshold disease states, we may gain insights into the biology of the disease, and perhaps identify therapeutic targets and opportunities for clinical interventions. For example, if gastrointestinal complaints are truly the consequence of aberrant AN gene expression, and contribute to disordered eating due to gastrointestinal distress, treatment of those symptoms may help alleviate other AN symptoms or prevent development of later AN(92,93). An understanding of the clinical consequences of AN gene expression can further augment the definition of AN, and could allow clinicians to more broadly identify individuals at greater risk of AN, or those who present with symptom constellations that do not yet meet the established diagnostic threshold.

Data Availability

Summary data is made available in the supplementary tables of this article.

Competing Interests

CMB reports: Shire (grant recipient, Scientific Advisory Board member); Idorsia (consultant); Pearson (author, royalty recipient). ML declares that over the past 36 months, he has received lecture honoraria from Lundbeck pharmaceutical (No other equity ownership, profit-sharing agreements, royalties, or patent). The remaining authors declare no competing interests.

Acknowledgements

We are deeply grateful for the mentorship of Pamela Sklar, whose guidance and wisdom we miss daily. We strive to continue her legacy of thoughtful, innovative, and collaborative science.

JJ and LMH were supported by funding from the Klarman Family Foundation (2019 Eating Disorders Research Grants Program) and the NIMH (R01MH118278). CMB is supported by NIMH (R01MH120170; R01MH119084; R01MH118278; U01 MH109528); Brain and Behavior Research Foundation Distinguished Investigator Grant; Swedish Research Council (Vetenskapsrådet, award: 538-2013-8864); Lundbeck Foundation (Grant no. R276-2018-4581).

This work was supported in part through the resources and staff expertise provided by the Charles Bronfman Institute for Personalized Medicine and The BioMe™ Biobank Program at the Icahn School of Medicine at Mount Sinai.

Research reported in this paper was supported by the Office of Research Infrastructure of the National Institutes of Health under award numbers S10OD018522 and S10OD026880. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.

Footnotes

  • ↵* Full list of PGC-ED working group authors is available in the supplementary materials.

  • An analysis to test for hidden case contamination has been added to the main manuscript and supplemental text. Supplemental files have been updated to reflect this additional analysis.

REFERENCES

  1. 1.↵
    Zipfel S, Giel KE, Bulik CM, Hay P, Schmidt U. Anorexia nervosa: aetiology, assessment, and treatment. Lancet Psychiatry. 2015 Dec;2(12):1099–111.
    OpenUrl
  2. 2.↵
    Bulik CM, Flatt R, Abbaspour A, Carroll I. Reconceptualizing anorexia nervosa. Psychiatry Clin Neurosci. 2019;73(9):518–25.
    OpenUrlCrossRefPubMed
  3. 3.↵
    Mitchell JE, Peterson CB. Anorexia Nervosa. N Engl J Med. 2020 02;382(14):1343–51.
    OpenUrl
  4. 4.↵
    Chesney E, Goodwin GM, Fazel S. Risks of all-cause and suicide mortality in mental disorders: a meta-review. World Psychiatry Off J World Psychiatr Assoc WPA. 2014 Jun;13(2):153–60.
    OpenUrl
  5. 5.↵
    Ulfvebrand S, Birgegård A, Norring C, Högdahl L, von Hausswolff-Juhlin Y. Psychiatric comorbidity in women and men with eating disorders results from a large clinical database. Psychiatry Res. 2015 Dec 15;230(2):294–9.
    OpenUrl
  6. 6.↵
    Root TL, Pinheiro AP, Thornton L, Strober M, Fernandez-Aranda F, Brandt H, et al. Substance use disorders in women with anorexia nervosa. Int J Eat Disord. 2010 Jan;43(1):14–21.
    OpenUrlPubMed
  7. 7.↵
    Raevuori A, Haukka J, Vaarala O, Suvisaari JM, Gissler M, Grainger M, et al. The increased risk for autoimmune diseases in patients with eating disorders. PloS One. 2014;9(8):e104845.
    OpenUrlCrossRefPubMed
  8. 8.↵
    Hedman A, Breithaupt L, Hübel C, Thornton LM, Tillander A, Norring C, et al. Bidirectional relationship between eating disorders and autoimmune diseases. J Child Psychol Psychiatry. 2019;60(7):803– 12.
    OpenUrl
  9. 9.↵
    Anorexia Nervosa Genetics Initiative, Eating Disorders Working Group of the Psychiatric Genomics Consortium, Watson HJ, Yilmaz Z, Thornton LM, Hübel C, et al. Genome-wide association study identifies eight risk loci and implicates metabo-psychiatric origins for anorexia nervosa. Nat Genet. 2019 Aug;51(8):1207–14.
    OpenUrl
  10. 10.↵
    Yilmaz Z, Hardaway JA, Bulik CM. Genetics and Epigenetics of Eating Disorders. Adv Genomics Genet. 2015;5:131–50.
    OpenUrlCrossRef
  11. 11.↵
    Levinson CA, Zerwas SC, Brosof LC, Thornton LM, Strober M, Pivarunas B, et al. Associations between dimensions of anorexia nervosa and obsessive-compulsive disorder: An examination of personality and psychological factors in patients with anorexia nervosa. Eur Eat Disord Rev J Eat Disord Assoc. 2019;27(2):161–72.
    OpenUrl
  12. 12.↵
    Kaye WH, Wierenga CE, Bailer UF, Simmons AN, Bischoff-Grethe A. Nothing tastes as good as skinny feels: the neurobiology of anorexia nervosa. Trends Neurosci. 2013 Feb 1;36(2):110–20.
    OpenUrlCrossRefPubMedWeb of Science
  13. 13.↵
    Speakman JR, Levitsky DA, Allison DB, Bray MS, Castro JM de, Clegg DJ, et al. Set points, settling points and some alternative models: theoretical options to understand how genes and environments combine to regulate body adiposity. Dis Model Mech. 2011 Nov 1;4(6):733–45.
    OpenUrlAbstract/FREE Full Text
  14. 14.
    Müller MJ, Bosy-Westphal A, Heymsfield SB. Is there evidence for a set point that regulates human body weight? F1000 Med Rep [Internet]. 2010 Aug 9 [cited 2020 Oct 20];2. Available from: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC2990627/
  15. 15.
    Hebebrand J, Muller TD, Holtkamp K, Herpertz-Dahlmann B. The role of leptin in anorexia nervosa: clinical implications. Mol Psychiatry. 2007 Jan;12(1):23–35.
    OpenUrlCrossRefPubMedWeb of Science
  16. 16.↵
    Hübel C, Yilmaz Z, Schaumberg KE, Breithaupt L, Hunjan A, Horne E, et al. Body composition in anorexia nervosa: Meta-analysis and meta-regression of cross-sectional and longitudinal studies. Int J Eat Disord. 2019;52(11):1205–23.
    OpenUrl
  17. 17.↵
    Igudesman D, Sweeney M, Carroll IM, Mayer-Davis EJ, Bulik CM. Gut-Brain Interactions: Implications for a Role of the Gut Microbiota in the Treatment and Prognosis of Anorexia Nervosa and Comparison to Type I Diabetes. Gastroenterol Clin North Am. 2019;48(3):343–56.
    OpenUrl
  18. 18.↵
    Huckins LM, Hatzikotoulas K, Southam L, Thornton LM, Steinberg J, Aguilera-Mckay F, et al. Investigation of common, low-frequency and rare genome-wide variation in anorexia nervosa. Mol Psychiatry. 2018;23(5).
  19. 19.↵
    Huckins LM, Dobbyn A, McFadden W, Ruderfer D, Wang W, Gamazon E, et al. Identifying tissues implicated in Anorexia Nervosa using Transcriptomic Imputation. bioRxiv. 2018;265017.
  20. 20.↵
    Hinney A, Kesselmeier M, Jall S, Volckmar A-L, Föcker M, Antel J, et al. Evidence for three genetic loci involved in both anorexia nervosa risk and variation of body mass index. Mol Psychiatry. 2017;22(2).
  21. 21.↵
    Duncan L, Yilmaz Z, Gaspar H, Walters R, Goldstein J, Anttila V, et al. Significant Locus and Metabolic Genetic Correlations Revealed in Genome-Wide Association Study of Anorexia Nervosa. Am J Psychiatry. 2017 01;174(9):850–8.
    OpenUrl
  22. 22.↵
    Munn-Chernoff MA, Johnson EC, Chou Y-L, Coleman JRI, Thornton LM, Walters RK, et al. Shared genetic risk between eating disorder- and substance-use-related phenotypes: Evidence from genome-wide association studies. Addict Biol. 2020 Feb 16;e12880.
  23. 23.↵
    Hübel C, Gaspar HA, Coleman JRI, Finucane H, Purves KL, Hanscombe KB, et al. Genomics of body fat percentage may contribute to sex bias in anorexia nervosa. Am J Med Genet Part B Neuropsychiatr Genet Off Publ Int Soc Psychiatr Genet. 2019;180(6):428–38.
    OpenUrl
  24. 24.↵
    C H, M A, M H, Rj L, G B, Cm B, et al. Binge-eating disorder, anorexia nervosa, and constitutional thinness differ in their associations with anthropometric and psychiatric polygenic scores. 2020 Mar 26 [cited 2020 Aug 20]; Available from: https://europepmc.org/article/ppr/ppr126993
  25. 25.↵
    Smoller JW. The use of electronic health records for psychiatric phenotyping and genomics. Am J Med Genet Part B Neuropsychiatr Genet Off Publ Int Soc Psychiatr Genet. 2018;177(7):601–12.
    OpenUrl
  26. 26.
    Denny JC, Ritchie MD, Basford MA, Pulley JM, Bastarache L, Brown-Gentry K, et al. PheWAS: demonstrating the feasibility of a phenome-wide scan to discover gene-disease associations. Bioinforma Oxf Engl. 2010 May;26(9):1205–10.
    OpenUrl
  27. 27.↵
    Pendergrass SA, Brown-Gentry K, Dudek SM, Torstenson ES, Ambite JL, Avery CL, et al. The use of phenome-wide association studies (PheWAS) for exploration of novel genotype-phenotype relationships and pleiotropy discovery. Genet Epidemiol. 2011 Jul;35(5):410–22.
    OpenUrlCrossRefPubMed
  28. 28.↵
    Zheutlin AB, Dennis J, Karlsson Linnér R, Moscati A, Restrepo N, Straub P, et al. Penetrance and Pleiotropy of Polygenic Risk Scores for Schizophrenia in 106,160 Patients Across Four Health Care Systems. Am J Psychiatry. 2019 Aug 16;176(10):846–55.
    OpenUrlCrossRefPubMed
  29. 29.
    Leppert B, Millard LAC, Riglin L, Davey Smith G, Thapar A, Tilling K, et al. A cross-disorder PRS-pheWAS of 5 major psychiatric disorders in UK Biobank. PLoS Genet. 2020 May;16(5):e1008185.
    OpenUrl
  30. 30.↵
    Denny JC, Bastarache L, Ritchie MD, Carroll RJ, Zink R, Mosley JD, et al. Systematic comparison of phenome-wide association study of electronic medical record data and genome-wide association study data. Nat Biotechnol. 2013 Dec;31(12):1102–10.
    OpenUrlCrossRefPubMed
  31. 31.↵
    Lonsdale J, Thomas J, Salvatore M, Phillips R, Lo E, Shad S, et al. The Genotype-Tissue Expression (GTEx) project. Nat Genet. 2013 Jun;45(6):580–5.
    OpenUrlCrossRefPubMed
  32. 32.↵
    Hoffman GE, Bendl J, Voloudakis G, Montgomery KS, Sloofman L, Wang Y-C, et al. CommonMind Consortium provides transcriptomic and epigenomic data for Schizophrenia and Bipolar Disorder. Sci Data. 2019 Sep 24;6(1):180.
    OpenUrl
  33. 33.↵
    Battle A, Mostafavi S, Zhu X, Potash JB, Weissman MM, McCormick C, et al. Characterizing the genetic basis of transcriptome diversity through RNA-sequencing of 922 individuals. Genome Res. 2014 Jan;24(1):14–24.
    OpenUrlAbstract/FREE Full Text
  34. 34.↵
    Gamazon ER, Wheeler HE, Shah KP, Mozaffari SV, Aquino-Michaels K, Carroll RJ, et al. A gene-based association method for mapping traits using reference transcriptome data. Nat Genet. 2015 Sep;47(9):1091–8.
    OpenUrlCrossRefPubMed
  35. 35.↵
    Barbeira AN, Dickinson SP, Bonazzola R, Zheng J, Wheeler HE, Torres JM, et al. Exploring the phenotypic consequences of tissue specific gene expression variation inferred from GWAS summary statistics. Nat Commun. 2018 May 8;9(1):1825.
    OpenUrlCrossRefPubMed
  36. 36.↵
    Huckins LM, Dobbyn A, Ruderfer DM, Hoffman G, Wang W, Pardiñas AF, et al. Gene expression imputation across multiple brain regions provides insights into schizophrenia risk. Nat Genet. 2019;51(4):659–74.
    OpenUrlCrossRefPubMed
  37. 37.↵
    Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MAR, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007 Sep;81(3):559–75.
    OpenUrlCrossRefPubMed
  38. 38.↵
    Chang CC, Chow CC, Tellier LC, Vattikuti S, Purcell SM, Lee JJ. Second-generation PLINK: rising to the challenge of larger and richer datasets. GigaScience. 2015;4:7.
    OpenUrlCrossRefPubMed
  39. 39.↵
    Fromer M, Roussos P, Sieberts SK, Johnson JS, Kavanagh DH, Perumal TM, et al. Gene expression elucidates functional impact of polygenic risk for schizophrenia. Nat Neurosci. 2016;19(11):1442–53.
    OpenUrlCrossRefPubMed
  40. 40.↵
    PheWAS/PheWAS [Internet]. PheWAS; 2019 [cited 2019 Nov 26]. Available from: https://github.com/PheWAS/PheWAS
  41. 41.↵
    Willer CJ, Li Y, Abecasis GR. METAL: fast and efficient meta-analysis of genomewide association scans. Bioinformatics. 2010 Sep 1;26(17):2190–1.
    OpenUrlCrossRefPubMedWeb of Science
  42. 42.↵
    Verma A, Bradford Y, Dudek S, Lucas AM, Verma SS, Pendergrass SA, et al. A simulation study investigating power estimates in phenome-wide association studies. BMC Bioinformatics. 2018 Apr 4;19(1):120.
    OpenUrlCrossRefPubMed
  43. 43.↵
    Wu P, Gifford A, Meng X, Li X, Campbell H, Varley T, et al. Developing and Evaluating Mappings of ICD-10 and ICD-10-CM Codes to PheCodes. bioRxiv. 2019 Jul 3;462077.
  44. 44.↵
    de Lange KM, Moutsianas L, Lee JC, Lamb CA, Luo Y, Kennedy NA, et al. Genome-wide association study implicates immune activation of multiple integrin genes in inflammatory bowel disease. Nat Genet. 2017 Feb;49(2):256–61.
    OpenUrlCrossRefPubMed
  45. 45.↵
    Mårild K, Størdal K, Bulik CM, Rewers M, Ekbom A, Liu E, et al. Celiac Disease and Anorexia Nervosa: A Nationwide Study. Pediatrics. 2017 May;139(5).
  46. 46.↵
    S Z, Jt L, L P, Lm T, M Q, Sv K, et al. Eating Disorders, Autoimmune, and Autoinflammatory Disease. Pediatrics [Internet]. 2017 Nov 9 [cited 2020 Nov 16];140(6). Available from: https://europepmc.org/article/med/29122972
  47. 47.↵
    Gibson D, Mehler PS. Anorexia Nervosa and the Immune System—A Narrative Review. J Clin Med. 2019 Nov;8(11):1915.
    OpenUrl
  48. 48.↵
    Stahl EA, Breen G, Forstner AJ, McQuillin A, Ripke S, Trubetskoy V, et al. Genome-wide association study identifies 30 loci associated with bipolar disorder. Nat Genet. 2019 May;51(5):793.
    OpenUrlCrossRefPubMed
  49. 49.↵
    Schizophrenia Working Group of the Psychiatric Genomics Consortium. Biological insights from 108 schizophrenia-associated genetic loci. Nature. 2014 Jul;511(7510):421–7.
    OpenUrlCrossRefPubMedWeb of Science
  50. 50.↵
    The Autism Spectrum Disorders Working Group of The Psychiatric Genomics Consortium. Meta-analysis of GWAS of over 16,000 individuals with autism spectrum disorder highlights a novel locus at 10q24.32 and a significant overlap with schizophrenia. Mol Autism. 2017 Dec;8(1):21.
    OpenUrlCrossRefPubMed
  51. 51.↵
    Zhu Z, Zhu X, Liu C-L, Shi H, Shen S, Yang Y, et al. Shared genetics of asthma and mental health disorders: a large-scale genome-wide cross-trait analysis. Eur Respir J [Internet]. 2019 Dec 1 [cited 2020 May 21];54(6). Available from: https://erj.ersjournals.com/content/54/6/1901507-0
  52. 52.↵
    Marchese S, Hoffman G, Pietrzak RH, Lebovitch D, Daskalakis NP, Yehuda R, et al. Altered gene expression and PTSD symptom dimensions in World Trade Center responders. In Prep. 2020;
  53. 53.↵
    Baselmans BML, Jansen R, Ip HF, van Dongen J, Abdellaoui A, van de Weijer MP, et al. Multivariate Genome-wide and integrated transcriptome and epigenome-wide analyses of the Well-being spectrum [Internet]. Genetics; 2017 Mar [cited 2020 May 21]. Available from: http://biorxiv.org/lookup/doi/10.1101/115915
  54. 54.↵
    Yang X, Sun J, Gao Y, Tan A, Zhang H, Hu Y, et al. Genome-Wide Association Study for Serum Complement C3 and C4 Levels in Healthy Chinese Subjects. PLOS Genet. 2012 Sep 13;8(9):e1002916.
    OpenUrlCrossRefPubMed
  55. 55.↵
    Sekar A, Bialas AR, de Rivera H, Davis A, Hammond TR, Kamitaki N, et al. Schizophrenia risk from complex variation of complement component 4. Nature. 2016;530(7589):177–83.
    OpenUrlCrossRefPubMed
  56. 56.↵
    Stephan AH, Barres BA, Stevens B. The complement system: an unexpected role in synaptic pruning during development and disease. Annu Rev Neurosci. 2012;35:369–89.
    OpenUrlCrossRefPubMedWeb of Science
  57. 57.↵
    Li B, Mao Y, Wang Z, Chen Y, Wang Y, Zhai C, et al. CLIC1 Promotes the Progression of Gastric Cancer by Regulating the MAPK/AKT Pathways. Cell Physiol Biochem. 2018;46(3):907–24.
    OpenUrlCrossRef
  58. 58.↵
    Carlini V, Verduci I, Cianci F, Cannavale G, Fenoglio C, Galimberti D, et al. CLIC1 Protein Accumulates in Circulating Monocyte Membrane during Neurodegeneration. Int J Mol Sci. 2020 Jan;21(4):1484.
    OpenUrl
  59. 59.↵
    Facal F, Costas J. Evidence of association of the DISC1 interactome gene set with schizophrenia from GWAS. Prog Neuropsychopharmacol Biol Psychiatry. 2019 20;95:109729.
    OpenUrlCrossRef
  60. 60.↵
    Milton RH, Abeti R, Averaimo S, DeBiasi S, Vitellaro L, Jiang L, et al. CLIC1 function is required for beta-amyloid-induced generation of reactive oxygen species by microglia. J Neurosci Off J Soc Neurosci. 2008 Nov 5;28(45):11488–99.
    OpenUrlAbstract/FREE Full Text
  61. 61.↵
    Novarino G, Fabrizi C, Tonini R, Denti MA, Malchiodi-Albedi F, Lauro GM, et al. Involvement of the intracellular ion channel CLIC1 in microglia-mediated beta-amyloid-induced neurotoxicity. J Neurosci Off J Soc Neurosci. 2004 Jun 9;24(23):5322–30.
    OpenUrlAbstract/FREE Full Text
  62. 62.↵
    Swarup V, Chang TS, Duong DM, Dammer EB, Dai J, Lah JJ, et al. Identification of Conserved Proteomic Networks in Neurodegenerative Dementia. Cell Rep. 2020 Jun 23;31(12):107807.
    OpenUrl
  63. 63.↵
    Lee YH, Bae S-C, Choi SJ, Ji JD, Song GG. Genome-wide pathway analysis of genome-wide association studies on systemic lupus erythematosus and rheumatoid arthritis. Mol Biol Rep. 2012 Dec 1;39(12):10627–35.
    OpenUrlCrossRefPubMed
  64. 64.
    Nilsson IAK, Millischer V, Göteson A, Hübel C, Thornton LM, Bulik CM, et al. Aberrant inflammatory profile in acute but not recovered anorexia nervosa. Brain Behav Immun. 2020 Aug 1;88:718–24.
    OpenUrl
  65. 65.↵
    1. Thompson JK
    Bulik CM. Genetic and Biological Risk Factors. In: Thompson JK, editor. Handbook of Eating Disorders and Obesity [Internet]. Hoboken, NJ: Wiley; 2004. p. 3. Available from: https://books.google.com/books/about/Handbook%7B_%7Dof%7B_%7DEating%7B_%7DDisorders%7B_%7Dand%7B_%7DObesity.html?id=KZhHAAAAMAAJ%7B&%7Dpgis=1
  66. 66.↵
    Shroff H, Reba L, Thornton LM, Tozzi F, Klump KL, Berrettini WH, et al. Features associated with excessive exercise in women with eating disorders. Int J Eat Disord. 2006 Sep;39(6):454–61.
    OpenUrlCrossRefPubMed
  67. 67.↵
    Achamrah N, Coëffier M, Déchelotte P. Physical activity in patients with anorexia nervosa. Nutr Rev. 2016 May 1;74(5):301–11.
    OpenUrlCrossRefPubMed
  68. 68.↵
    Strober M, Freeman R, Morrell W. The long-term course of severe anorexia nervosa in adolescents: Survival analysis of recovery, relapse, and outcome predictors over 10–15 years in a prospective study. Int J Eat Disord. 1997;22(4):339–60.
    OpenUrlCrossRefPubMedWeb of Science
  69. 69.↵
    Hübel C, Gaspar HA, Coleman JRI, Hanscombe KB, Purves K, Prokopenko I, et al. Genetic correlations of psychiatric traits with body composition and glycemic traits are sex-and age-dependent. Nat Commun. 2019 18;10(1):5765.
    OpenUrl
  70. 70.↵
    Solmi M, Veronese N, Correll CU, Favaro A, Santonastaso P, Caregaro L, et al. Bone mineral density, osteoporosis, and fractures among people with eating disorders: a systematic review and meta-analysis. Acta Psychiatr Scand. 2016 May;133(5):341–51.
    OpenUrlCrossRefPubMed
  71. 71.↵
    Etain B, Jamain S, Milhiet V, Lajnef M, Boudebesse C, Dumaine A, et al. Association between circadian genes, bipolar disorders and chronotypes. Chronobiol Int. 2014 Aug;31(7):807–14.
    OpenUrl
  72. 72.
    Melo MCA, Abreu RLC, Linhares Neto VB, de Bruin Pfc, de Bruin Vms. Chronotype and circadian rhythm in bipolar disorder: A systematic review. Sleep Med Rev. 2017 Aug;34:46–58.
    OpenUrl
  73. 73.
    Huckins L, Dobbyn A, McFadden W, Wang W, Ruderfer D, Hoffman G, et al. Transcriptomic Imputation of Bipolar Disorder and Bipolar subtypes reveals 29 novel associated genes. bioRxiv. 2017 Nov;222786.
  74. 74.
    Meyrer R, Demling J, Kornhuber J, Nowak M. Effects of night shifts in bipolar disorders and extreme morningness. Bipolar Disord. 2009 Dec;11(8):897–9.
    OpenUrl
  75. 75.
    Alloy LB, Ng TH, Titone MK, Boland EM. Circadian Rhythm Dysregulation in Bipolar Spectrum Disorders. Curr Psychiatry Rep. 2017 Apr;19(4):21.
    OpenUrl
  76. 76.
    Boivin DB. Influence of sleep-wake and circadian rhythm disturbances in psychiatric disorders. J Psychiatry Neurosci JPN. 2000 Nov;25(5):446–58.
    OpenUrl
  77. 77.
    Bunney BG, Bunney WE. Mechanisms of Rapid Antidepressant Effects of Sleep Deprivation Therapy: Clock Genes and Circadian Rhythms. Biol Psychiatry. 2013 Jun;73(12):1164–71.
    OpenUrlCrossRefPubMedWeb of Science
  78. 78.
    Bunney BG, Li JZ, Walsh DM, Stein R, Vawter MP, Cartagena P, et al. Circadian dysregulation of clock genes: clues to rapid treatments in major depressive disorder. Mol Psychiatry. 2015 Feb;20(1):48–55.
    OpenUrlCrossRefPubMed
  79. 79.
    Karatsoreos IN. Links between Circadian Rhythms and Psychiatric Disease. Front Behav Neurosci. 2014;8:162.
    OpenUrlCrossRefPubMed
  80. 80.
    Kripke DF, Mullaney DJ, Atkinson M, Wolf S. Circadian rhythm disorders in manic-depressives. Biol Psychiatry. 1978 Jun;13(3):335–51.
    OpenUrlPubMedWeb of Science
  81. 81.
    Mansour HA, Talkowski ME, Wood J, Chowdari KV, McClain L, Prasad K, et al. Association study of 21 circadian genes with bipolar I disorder, schizoaffective disorder, and schizophrenia. Bipolar Disord. 2009 Nov;11(7):701–10.
    OpenUrl
  82. 82.↵
    Murray G, Allen NB, Trinder J, Burgess H. Is weakened circadian rhythmicity a characteristic of neuroticism? J Affect Disord. 2002 Dec;72(3):281–9.
    OpenUrlCrossRefPubMed
  83. 83.↵
    Allison KC, Spaeth A, Hopkins CM. Sleep and Eating Disorders. Curr Psychiatry Rep. 2016;18(10):92.
    OpenUrl
  84. 84.↵
    Arble DM, Bass J, Laposky AD, Vitaterna MH, Turek FW. Circadian timing of food intake contributes to weight gain. Obesity. 2009;17(11):2100–2.
    OpenUrl
  85. 85.↵
    Herpertz S, Albers N, Wagner R, Pelz B, Köpp W, Mann K, et al. Longitudinal changes of circadian leptin, insulin and cortisol plasma levels and their correlation during refeeding in patients with anorexia nervosa. Eur J Endocrinol. 2000 Apr;142(4):373–9.
    OpenUrlAbstract
  86. 86.↵
    Zucker NL, Bulik CM. On bells, saliva, and abdominal pain or discomfort: Early aversive visceral conditioning and vulnerability for anorexia nervosa. Int J Eat Disord. 2020;53(4):508–12.
    OpenUrl
  87. 87.↵
    Adams DM, Reay WR, Geaghan MP, Cairns MJ. Investigation of glycaemic traits in psychiatric disorders using Mendelian randomisation revealed a causal relationship with anorexia nervosa. Neuropsychopharmacology. 2020 Sep 13;1–10.
  88. 88.↵
    Ilyas A, Hübel C, Stahl D, Stadler M, Ismail K, Breen G, et al. The metabolic underpinning of eating disorders: A systematic review and meta-analysis of insulin sensitivity. Mol Cell Endocrinol. 2019 Nov 1;497:110307.
    OpenUrl
  89. 89.↵
    Abdulkadir M, Herle M, De Stavola BL, Hübel C, Santos Ferreira DL, Loos RJF, et al. Polygenic Score for Body Mass Index Is Associated with Disordered Eating in a General Population Cohort. J Clin Med [Internet]. 2020 Apr 21 [cited 2020 Nov 3];9(4). Available from: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC7231239/
  90. 90.↵
    Nagata JM, Braudt DB, Domingue BW, Bibbins-Domingo K, Garber AK, Griffiths S, et al. Genetic risk, body mass index, and weight control behaviors: Unlocking the triad. Int J Eat Disord. 2019;52(7):825–33.
    OpenUrl
  91. 91.↵
    Kolb R, Sutterwala FS, Zhang W. Obesity and cancer: inflammation bridges the two. Curr Opin Pharmacol. 2016 Aug 1;29:77–89.
    OpenUrlCrossRef
  92. 92.↵
    Wiklund CA, Kuja-Halkola R, Thornton LM, Hübel C, Leppä V, Bulik CM. Prolonged constipation and diarrhea in childhood and disordered eating in adolescence. J Psychosom Res. 2019;126:109797.
    OpenUrl
  93. 93.↵
    Riedlinger C, Schmidt G, Weiland A, Stengel A, Giel KE, Zipfel S, et al. Which Symptoms, Complaints and Complications of the Gastrointestinal Tract Occur in Patients With Eating Disorders? A Systematic Review and Quantitative Analysis. Front Psychiatry. 2020;11:195.
    OpenUrl
Back to top
Posted March 13, 2021.
Download PDF

Supplementary Material

Data/Code
Email
The Phenome-wide Consequences of Anorexia Nervosa Genes
Jessica S. Johnson, Alanna C. Cote, Amanda Dobbyn, Laura G. Sloofman, Jiayi Xu, Liam Cotter, Alexander W. Charney, Eating Disorders Working Group of the Psychiatric Genomics Consortium, Jennifer Jordan, Martin Kennedy, Mikael Landén, Sarah L Maguire, Nicholas G Martin, Preben Bo Mortensen, Cynthia M. Bulik, Laura M. Huckins
medRxiv 2021.02.12.21250941; doi: https://doi.org/10.1101/2021.02.12.21250941
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools

Subject Area