Next Article in Journal
Genetic Causes of Oculocutaneous Albinism in Pakistani Population
Next Article in Special Issue
Family-Based Genome-Wide Association Study of Autism Spectrum Disorder in Middle Eastern Families
Previous Article in Journal
Genome-Wide Association Studies and Prediction of Tan Spot (Pyrenophora tritici-repentis) Infection in European Winter Wheat via Different Marker Platforms
Previous Article in Special Issue
Sensory Reactivity Symptoms Are a Core Feature of ADNP Syndrome Irrespective of Autism Diagnosis
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 12
Issue 4
10.3390/genes12040491
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dissection of the Genetic Association between Anorexia Nervosa and Obsessive–Compulsive Disorder at the Network and Cellular Levels

by
Weichen Song
,
Weidi Wang
,
Shunying Yu
and
Guan Ning Lin
*
School of Biomedical Engineering, Shanghai Jiao Tong University, Shanghai 200030, China
*
Author to whom correspondence should be addressed.
Genes 2021, 12(4), 491; https://doi.org/10.3390/genes12040491
Submission received: 25 February 2021 / Revised: 23 March 2021 / Accepted: 24 March 2021 / Published: 27 March 2021
(This article belongs to the Special Issue Genomics of Neuropsychiatric Disorders)
Browse Figures
Versions Notes

Abstract

:
Anorexia nervosa (AN) and obsessive–compulsive disorder (OCD) exhibit a high co-morbidity rate, similar symptoms, and a shared genetic basis. However, an understanding of the specific underlying mechanisms of these commonalities is currently limited. Here, we collected Genome-Wide Association Analysis results for AN and OCD, and obtained genes hit by the top SNPs as the risk genes. We then carried out an integrative coexpression network analysis to explore the convergence and divergence of AN and OCD risk genes. At first, we observed that the AN risk genes were enriched in coexpression modules that involved extracellular matrix functions and highly are expressed in the postnatal brain, limbic system, and non-neuronal cell types, while the OCD risk genes were enriched in modules of synapse function, the prenatal brain, cortex layers, and neurons. Next, by comparing the expressions from the eating disorder and OCD postmortem patient brain tissues, we observed both disorders have similar prefrontal cortex expression alterations influencing the synapse transmission, suggesting that the two diseases could have similar functional pathways. We found that the AN and OCD risk genes had distinct functional and spatiotemporal enrichment patterns but carried similar expression alterations as a disease mechanism, which may be one of the key reasons they had similar but not identical clinical phenotypes.

1. Introduction

Anorexia nervosa (AN) and obsessive–compulsive disorder (OCD) are two distinct neuropsychiatric disorders but share a genetic basis [1,2,3]. Epidemiological studies have found that the diagnosis of OCD would increase the risk of AN in later life (risk ratio = 3.6) and vice versa (risk ratio = 9.6) [4]. Similarly, the prevalence of OCD is higher in AN patients than in the population [5]. The prevalence of AN is also higher in OCD patients [6], which indicates significant co-morbidity between the two diseases. In accordance with the epidemiological associations, clinical studies revealed that AN and OCD have highly similar phenotypes, such as perfectionism, rigidity, persistency, obsessive thoughts [7], excessive [8], and ritual behaviors [9]. Clinicians have also illustrated the characteristic phenotypes shared by AN and OCD: temporary relief for anxiety after obsessive–compulsive behaviors, followed by increasing long-term anxiety and fear [10]. These significant overlaps between AN and OCD has led some researchers to view them as part of an integrated spectrum, obsessive–compulsive spectrum disorders [11], indicating that they may be the different phenotypes of an identical pathology.
Genetic studies on AN and OCD have also suggested that they have shared underlying mechanisms. A case–control study has found a significant familial aggregation of AN and OCD [12], and the bivariate twin models from the Swedish National Patient Register have reported a significant genetic overlap between the two diseases [4], with a correlation coefficient of 0.52. Consistently, Linkage Disequilibrium Regression analysis using Genome-Wide Association Study (GWAS) data revealed a positive genetic correlation [13,14]. The most recent meta-analysis from Psychiatric Genomics Consortium [15] has revealed that AN and OCD shared a compulsive component of genetic architecture, which contributed to their significant genetic correlation (the correlation coefficient reached 0.50) and commonalities in obsessive–compulsive symptoms. These results have suggested that AN and OCD share a large proportion of their genetic basis and indicated that they might have a common pathology.
However, AN and OCD are clinically recognized as different disorders due to their profound difference in multiple aspects. By comparing AN and OCD patients, researchers have found a significant divergence with regard to cognition [16,17], emotion [18,19], and personality [20,21]. Furthermore, choices of treatments, such as deep brain stimulation [22], selective serotonin reuptake inhibitor [23], and cognition-behavior therapy [24], are different between AN and OCD patients. The situation is even complicated and worse in co-morbidity cases [4,24]. One possible explanation is that AN and OCD exhibit distinct genetic and transcriptomic characteristics at the regional and cellular levels, which served as the foundation of the differential response to cognitive intervention and regional stimulation. These discrepancies urge researchers to substantially investigate the convergence and divergence of the AN and OCD pathologies, to better understand them and thus gain more thorough insight into their genetic etiology.
To achieve this task, we have conducted a systematic analysis combining genetic data [25,26] and transcriptomic data [27,28] to dissect the shared and unique components of the genetic pathology of the two diseases. We first collected the risk genes implicated by the SNPs for the two diseases and explored the spatiotemporal, cell-type specific, and functional patterns of these genes at the network as well as genome-wide levels. As a result, we have found distinct and semblable coexpression patterns across developmental periods, cell types, and phenotype conditions in both disorders. So, these analyses provided insights into the convergence and divergence of AN and OCD, which could help further functional analysis and clinical personal medicine.

2. Materials and Methods

2.1. Data Collection

The brain transcriptome data of the OCD, eating disorder, and control groups were obtained from GSE60190 [28]. This dataset, provided by Jaffe et al., includes normalized expression array data of the prefrontal cortex from patients with eating disorders (ED; including three AN, seven bulimia nervosa patients, two binge–purge type of AN, and three eating disorder not otherwise specified), obsessive–compulsive disorder/obsessive–compulsive personality disorder or tics (OCD/OCPD/Tic; N = 16), and non-psychiatric controls (N = 102). We did not perform extra data preprocessing but reran the differential expression analysis (see section “Differential Expression Analysis”).
The GWAS summary data for AN was collected from Watson et al. [25]. We considered all Linkage Disequilibrium (LD) blocks with a genome-wide significant (p < 5 × 10−8) tag SNP as the target blocks and extracted all genes within the target blocks as the risk genes of AN. The range of LD blocks was estimated by the original GWAS by local LD structure and a p-value threshold. Similarly, the LD blocks for OCD were downloaded from Arnold et al. [26] (Table 1 of ref [26]; tagged region by the index SNP at r2 > 0.2 with 1KG as the reference panel) and the OCD risk genes were retrieved from the blocks, respectively. Only SNPs passing quality control in at least 500 cases and 500 controls were retained. Since this GWAS had limited statistical power, we set the significance threshold at 1 × 10−5 for the target blocks. To minimize the impact of false-positive risk genes, we excluded the most marginal genes for the large blocks. A total of 92 AN risk genes and 51 OCD risk genes were obtained. We downloaded brain transcriptome data during normal brain development from the BrainSpan portal [27].

2.2. Weighted Gene Co-Expression Network Analysis (WGCNA)

We divided the BrainSpan v1.0 [27] RPKM data, which spans 15 consecutive periods of neurodevelopment and adulthood from 5.7 postconceptual weeks (PCW) to 82 years. We used two subsets of prenatal and postnatal samples. For both subsets, only genes with non-zero RPKM in at least half of the samples and a coefficient of variance >0.3 were retained for WGCNA and were log10-transformed. A soft threshold was set at 12 for the prenatal subset and 16 for the postnatal subset by the pickSoftThreshold function in the WGCNA v1.70 [29] R package. Next, we built up the signed networks by the blockwiseModules function in WGCNA R package and partitioned genes into modules with a tree cut height of 0.3, minimum module size = 20, and maximum module size = 1000. We calculated the first principal component of each module, namely Module Eigengene (ME), to represent the overall expression levels of the module in each sample. To compare the temporal trajectory of each module, we fit a curve for ME across different developmental periods using the loess function in R v4.0.0.
After the detection of coexpression modules for both subsets, we tested if the AN and OCD risk genes were enriched in any of the modules by Fisher’s exact test. We removed the risk genes that were not presented in the BrainSpan data. The modules with the lowest p-value were extracted for further analysis. Since the OCD risk genes showed enrichments in one prenatal module (M6) and one postnatal module (M17), we performed EWCE v1.0 (see Section 2.4. “spatiotemporal expression analysis”) on the intersection of these two modules and found that they were more highly expressed in the prenatal period. Thus, we chose prenatal M6 for further analysis. Similarly, for the AN risk genes, we chose postnatal M2 for the following analysis. As a negative control, we repeated the enrichment analysis by genes related to height and BMI curated at traitDB (https://atlas.ctglab.nl/traitDB) (accessed on 10 March 2021). A lack of significant enrichment after p-value adjustment (nominal p > 0.05) was considered the evidence that our enrichment analysis did not introduce false-positive results.

2.3. Gene Ontology (GO) Analysis

We performed a GO enrichment analysis of the gene lists we were interested in by ClusterProfiler v3.12 [30] in the R package. We tested the genes of interest enriched in any GO-BP pathway by the hypergeometric test. Gene background was defined as all genes with a GO annotation. Only pathways with more than ten genes were included in our analysis. The p-values of the hypergeometric tests were adjusted for multiple testing by the Benjamin–Hochberg method. We only retained the pathways with an adjusted p < 0.05, size > 20 and < 500.

2.4. Spatiotemporal Expression Analysis

To analyze the expression patterns of the AN and OCD risk genes across different developmental periods and brain regions, we applied the EWCE v1.0 R package [31] on BrainSpan data, which is an enrichment tool that performed the bootstrap test on the gene specificity matrix based on BrainSpan data. The developmental periods were defined as previously described [32]: early embryo (8–12 weeks), the mid-early embryo (13–16 weeks), the mid-late embryo (17–21 weeks), and late embryo (24–37 weeks) for prenatal data; and infancy (born 3 years), childhood (4–11 years), adolescence (13–19 years), and adult (19 years) for postnatal data. Specificity was defined in the corresponding paper [31]. We first calculated the specificity matrix of the BrainSpan datasets by the “generate.celltype.data” function. Then, the enrichment of genes within the prenatal module M6 and postnatal module M2 was tested by the “bootstrap.enrichment.test” function with the number of bootstraps = 10,000. We then defined the background gene list as all genes annotated in the BrainSpan dataset.

2.5. Network Analysis

To construct the coexpression network for the prenatal module M6 and postnatal module M2, we first calculated the signed Topology Overlap Measures (TOM) [29] for both modules. Then, we retrieved all the genes that connected to at least one risk gene with TOM ≥ 0.1. These genes, together with the risk genes of AN or OCD, were used to construct the network with all connections between them (only connections with TOM ≥ 0.1). Hub genes were identified by two measurements: degree score within the network, and the module eigengene connectivity (kME, defined as the correlation coefficient between gene expression value and the module eigengene). The GO annotations were retrieved from the clusterProfiler v3.12 [30] R package.

2.6. Differential Expression Analysis

We obtained postmortem cortex transcriptome data from Jaffe et al. [28] (GSE60190) and compared the gene expression of the ED and OCD patients to that of a healthy control. GSE60190 contains normalized Illumina HumanHT-12 v3 (San Diego, CA, USA) microarray data collected on postmortem dorsolateral prefrontal cortexes. The acquisition and process of data from the case and control were performed similarly. We applied limma [33] on this transcriptome data to identify the genes that were significantly dysregulated in the cortex of patients with ED or OCD. We first fit a linear model with the following design, which controlled all the provided covariates and compared the ED and OCD groups to the control:
Expression ~ intercept + group + RIN + Ph + age + PMI + sex
Then, we computed empirical Bayesian statistics for comparison between the ED and control groups as well as between the OCD and control groups. We extracted all the genes with a Benjamini–Hochberg adjusted p-value < 0.05 as the Differentially Expressed Genes (DEGs).

2.7. Gene Set Enrichment Analysis

We applied the Gene Set Enrichment Analysis (GSEA) of the Reactome pathways [34] using the ReactomePA v1.34 [35] R package. The input to ReactomePA was the gene list ordered by fold change from the differential expression analysis. GSEA was repeated for both OCD vs control and ED vs control differential expression analysis. We extracted Normalized Enrichment Scores (NES) for each pathway to compare the biological significance between the expression patterns of OCD and ED. GSEA was run with the number of permutations = 10,000. We retained the pathways with an adjusted p < 0.05, and size > 20 and < 500.
Using the same ranked gene lists for OCD and ED, we analyzed the expression alteration of the previously identified M6 module (BrainSpan prenatal) and M2 module (BrainSpan postnatal) using the GSEA function from the clusterProfiler v3.12 [30] R package.

3. Results

To investigate the genetic association between AN and OCD, we first collected the implicated risk genes of AN and OCD identified from GWAS summary statistics [25,26] (Table S1). Since AN and OCD shared a profound genetic basis [1,2,3], our first question is whether their risk genes had a similar coexpression pattern. We utilized the spatially and temporally rich transcriptome data set from BrainSpan [27] and constructed coexpression networks for both the prenatal and postnatal periods. As a result, we identified 19 coexpression modules in the prenatal epoch and 21 in the postnatal epoch (Figure 1A).

3.1. AN and OCD Risk Genes Coexpression Modules Show Different Biological Functions in Prenatal and Postnatal Periods

We observed that the AN and OCD risk genes were enriched in different modules in both subsets. In the prenatal network modules, significantly more AN genes than expected were observed in module 10 (OR = 2.44, Fisher’s exact p = 0.0066), whereas significantly more OCD genes than expected were observed in module M6 (OR = 9.33, Fisher’s exact p = 0.0014) (Figure 1B). For the postnatal analysis, the AN genes were most significantly enriched in the postnatal network M2 module (OR = 2.65, Fisher’s exact p = 0.0013) and the OCD genes in the postnatal network module M17 (OR = 2.52, Fisher’s exact p = 0.012) (Figure 1B). As a negative control, genes related to height and BMI were not enriched in these modules (p > 0.05). These results suggested that the spatiotemporal coexpression patterns of the AN and OCD risk genes during normal brain development was not identical.
Next, we investigated whether modules enriched with the AN or OCD genes differed in biological functions. As shown in Figure 1C, we observed that both prenatal M6 and postnatal M17 modules co-expressed with OCD genes enriched in regulation of neuron projection development (prenatal M6: p-adjust = 1.45 × 10−4; postnatal M17: p-adjust = 1.85 × 10−9) and glutamatergic synaptic transmission (prenatal M6: p-adjust = 0.003; postnatal M17: p-adjust = 0.006). On the other hand, prenatal M10 and postnatal M2 modules co-expressed with AN genes (Figure 1D) had similar functions in extracellular structure organization (prenatal M10: p-adjust = 6.60 × 10−38; postnatal M2: p-adjust = 2.45 × 10−21) and immune-related pathways (prenatal M10: neutrophil activation, p-adjust = 4.31 × 10−20; postnatal M2: defense response to virus, p-adjust = 2.02 × 10−12).
Since the prenatal and postnatal modules of the same disease had similar biological significance, we asked whether these modules were similarly critical to the disease, or one of them served as the primary module of the disease. As shown in Figure S1C,D, the overlapped genes between prenatal M6 and postnatal M17 were preferentially expressed in the prenatal period (fold change = 1.21, p < 0.0001), whereas the overlapped genes between prenatal M10 and postnatal M2 were preferentially expressed in infancy (fold change = 1.33, p < 0.0001) and the adolescence period (fold change = 1.14, p < 0.0001). These results suggest that the genes co-expressed with the OCD risk genes are highly expressed in the prenatal period, and that with the AN risk genes are highly expressed in the postnatal period. These results confirmed prenatal M6 as the primary module for OCD, and postnatal M2 as the primary module for AN. Next, we focused our analyses on these two modules.

3.2. Coexpression Modules of the AN and OCD Risk Genes Exhibited Distinct Spatiotemporal and Cell-Type-Specific Expression Patterns

Having identified the primary modules as well as their biological functions for AN and OCD, we next compared these coexpression modules from the aspects of spatiotemporal expression patterns. Using the EWCE analysis technique [31], we found that the genes in prenatal M6 were expressed significantly higher than expected in the mid-early prenatal period (fold change = 1.3, p < 0.0001), given the fact that module eigengene (ME) of prenatal M6 peaked at mid-early (Figure 2A). Postnatal M2, however, was highly expressed in infancy (fold change = 1.22, p < 0.0001) and adolescence (fold change = 1.06, p < 0.0001) (Figure 2B). A divergence between the modules was also observed for brain regions, where prenatal M6 was enriched in the cerebral cortex (p < 0.0001 for all frontal, parietal, temporal, and occipital cortexes; Figure 2C) and postnatal M2 was enriched in the limbic system (p < 0.0001 for the amygdala, medial thalamus, striatum, and hippocampus; Figure 2D). Lastly, we explored the cell-type-specific expression patterns for these two modules using EWCE. We observed that prenatal M6 was more highly expressed in neurons in the embryo brain than expected (excitatory neuron: fold change = 1.62, p < 0.0001; gamma-Aminobutyric acid[GABA]-ergic neuron: fold change = 1.30, p < 0.0001; Figure 2D). We also found that postnatal M2 was more enriched in the endothelium than expected (fold change = 2.16, p < 0.0001) and in astrocytes (fold change = 1.37, p < 0.0001, Figure 2D) of the adult brain. These results suggest that the AN and OCD modules had distinct spatiotemporal and cell-type-specific expression patterns.

3.3. OCD and Eating Disorder Led to a Similar Brain Transcriptome Alteration Pattern

We have found the distinct functional and spatiotemporal characteristics of AN and OCD risk genes using transcriptome data of the healthy human brain. However, these characteristics may not be mirrored by transcriptomic alteration of brains of the patients with diseases [36], since the risk SNPs may influence normal functions by means other than transcriptomic dysregulation [37]. Thus, we further analyzed the cortex transcriptomic alteration of OCD and eating disorder (ED) using data from Jaffe et al. [28].
Unlike the significant observation in differential coexpression patterns of risk genes in AN and OCD, the log-transformed fold change (logFC) of expression for ED and OCD patients was highly convergent (correlation coefficient = 0.688, p < 2.2 × 10−16). Among these top differential expression genes, we found upregulated genes involved in synapse organization (CPLX1: logFC_OCD = 0.51, logFC_ED = 0.92; DLG4: logFC_OCD = 0.53, logFC_ED = 0.83) and downregulated genes involved in cell adhesion (NINJ2: logFC_OCD = −0.60, logFC_ED = −0.98; COL4A5: logFC_OCD = −0.38, logFC_ED = −0.97). Setting the significance cutoff of p-adjust as 0.05, we identified 1118 significant differentially expressed genes (DEGs) for ED and 238 DEGs for OCD. The significant overlapped DEGs count between them (50; p by hypergeometric test <2.2 × 10−16, Figure 3B) also supported the fact that transcriptome alteration of ED and OCD were nearly parallel. In summary, ED and OCD brought about similar brain expression alterations.
We further applied GSEA to verify that the functional impacts of the expression alteration were in the same pattern between both diseases (Figure 3C and Table S2). As expected, the Normalized Enrichment Score (NES) was also convergent between two diseases (correlation coefficient = 0.649, p < 2.2 × 10−16). The upregulated genes of both diseases were enriched in Transmission across Chemical Synapses (NES_OCD = 1.67, NES_ED = 2.01, Figure 3C) and Phase 0—rapid depolarization (NES_OCD = 1.58, NES_ED = 2.08, Figure 3D), supporting the transcriptomic disruption of synapse functions in both diseases. Moreover, for the downregulated genes, the most significant enrichment was found for the SRP-dependent co-translational protein targeting membranes (NES_OCD = −1.59, NES_ED = −1.88, Figure 3D). These findings demonstrated that the convergent functional pathway from the distinct disorders is unmasked by the SNP risk DEGs across ED and OCD patient samples.

3.4. LAMB2, HTRA1, and KCTD12 as Hubs for AN and OCD Networks

Having compared the characteristics of AN and OCD genetic etiology from multiple aspects, we next aimed to integrate them into genetic networks so that the whole picture of AN and OCD genetics could be drawn. To achieve this goal, we built up the coexpression network of prenatal M6 for OCD and postnatal M2 for AN by Topology Overlap Measures (TOM, see METHOD).
In the network of AN (postnatal M2, Figure 4A), we found the significant enrichment of AN risk genes (p = 0.0013, Figure 1D), genes downregulated in the cortex of ED patients (p < 2.2 × 10−16), genes involved in extracellular structure organization (p-adjust = 2.45 × 10−21, Figure 1D), and defense response to virus (p-adjust = 2.02 × 10−12, Figure 1D). One of the AN risk genes, LAMB2, appeared at the center of the network (degree = 154). It also had high intramodular connectivity (kME = 0.82) and was involved in extracellular structure organization, both indicating it as a hub role gene of postnatal M2. Another extracellular structure organization gene, HTRA1, had a similar central position (kME = 0.86, degree = 182). Unlike LAMB2, HTRA1 was not impacted by AN risk SNPs but was significantly downregulated in the cortex of ED patients (logFC = −0.89, p-adjust = 0.01). As the highest degree of SLC25A20 reached 194, which encodes a mitochondrial solute carrier, it may also be considered a hub risk gene (kME = 0.81).
In the network of OCD (prenatal M6, Figure 4B), we also found the significant enrichment of OCD risk genes (p = 0.0014, Figure 1C) and synapse organization (p-adjust = 3.55 × 10−5), but not of differential expression genes in the cortex of OCD patients. Among these risk genes of OCD, we found KCTD12, as a potential hub with kME = 0.81. Although none of the synapse organization genes fell under the OCD risk genes, one of them (TMEM108) had a central position in the network (kME = 0.91, degree = 60). The only differential expression gene in the cortex of OCD patients in the network was FEZF2, which is located at the border of the network (kME = 0.65, degree = 4). It is of note that FEZF2 is a gene with a potentially high severity level in psychiatric disorders, as evaluated by our published database [38].

4. Discussion

In our current study, we studied the convergence and divergence of the genetic basis underlying AN and OCD from the gene and network levels. We have found that although AN and OCD patients exhibited similar cerebral transcriptomic alterations, their risk genes affected networks with distinct functions and spatiotemporal expression patterns. These results extend the understanding of the shared pathology [13] between the two diseases and partially explain why their symptoms resemble each other but are not identical [8,10].
In the above results, one of the most significant discoveries of this study is the inconsistency between the genome and transcriptome abnormality of AN and OCD: they are caused by distinct risk gene networks, but exhibit a similar transcriptome alteration. The shared transcriptomic patterns across psychiatric diseases have been demonstrated by the meta-analysis of Gandal et al. [39], who found that the prefrontal cortex of patients with autism, schizophrenia, bipolar disorder, and depression exhibited similar differential expression patterns, and their differential expression genes were enriched in similar coexpression modules. Although the expression data of AN and OCD [28] had a relatively smaller sample size, we found a convergent transcriptomic alteration between the two diseases by using the same approach (correlation analysis between logFC) as Gandal et al. This phenomenon approves of the current effort of cross-disorder studies that aim at common symptoms and pathologies shared by psychiatric disorders [40]. However, the divergent roles of risk genes in the different disorders seem contradictory to the convergent transcriptomic patterns. An integrative analysis by Shohat et al. [41] showed that rare mutations found in intellectual disability (ID), autism, and schizophrenia (SCZ) had distinct functions in chromatin organization and neuron development and exhibited different regions and cell-type enrichment. We also found that the AN and OCD risk genes were enriched in distinct coexpression modules.
This contradiction may arise from various influencing factors. First of all, most of the psychiatric transcriptome analyses, including Jaffe et al. and Gandal et al., focused on the same region (prefrontal cortex), whereas genetic factors underlying different diseases pointed toward different brain regions [36,41,42,43]. Similarly, various imaging [44,45] and functional studies [46] also support the idea that the critical region and neural circuits of different diseases are distinct. It could be inferred that genetic variations first introduced transcriptomic abnormalities in the primary brain regions responsible for different diseases, and then caused additional influences in the prefrontal cortex, which are similar across diseases. Secondly, as pointed out by Gandal et al. [39], the linkage between genetic risk and expression alteration is complex and indirect, often a multi-step cascade of different biological processes. Among this cascade, a lot of influencing factors, both biological and environmental, are shared among psychiatric diseases. For example, chromatin modification has been indicated in the pathology of SCZ [47], autism [48], ID [49], OCD [50], and AN [51], and has the capacity of regulating gene expression levels [52]. Stress events and inflammation, two environmental factors that significantly impact the gene expression of central nervous system [53,54], contribute to the risk of various psychiatric diseases [55,56], especially AN [57] and OCD [58]. Lastly, damaging genetic variations do not necessarily lead to expression abnormality. As revealed by the probability of the loss-of-function intolerant (pLi) score [59] and haploinsufficiency score [60], a gene can be sensitive either to truncating variants or dosage change (i.e., expression alteration), or both, without apparent causative relationships. Thus, not all dosage changes of genes are pathologically meaningful, nor do all damaging genetic variations come into effect by causing dosage changes, which explains the divergence between risk genes and transcriptomic alterations.
Using integrative approaches, we found that AN and OCD had different coexpression patterns concerning developmental periods, brain regions, and cell types. Epidemiological studies have demonstrated that the onset age of AN was usually in late adolescence (16~18 years old) [61,62], consistent with our finding that AN risk genes were enriched in modules that peaked in adolescence. For OCD, its risk genes peaked in mid-prenatal, which might partially explain the existence of the early-onset subtype of OCD [63]. Although previous studies did not observe the differences in brain regions (both affecting cortical-limbic circuits [64,65]) or cell types (both affecting serotonergic neuron [66,67]) involved in AN and OCD, we showed that the risk genes of the two diseases have such a distinction. In fact, the AN risk genes mostly affect the subcortical regions and endothelium, whereas the OCD risk genes mainly affect cortices and neurons. These results suggest that the two diseases may have different “primary” regions and cells where the original dysfunctions emerge.
The more encouraging results came from the followed functional enrichment analysis and network analysis. By integration of developmental coexpression patterns, disease GWAS implication and disease transcriptomic alteration, we found that the AN risk genes enriched in the network of extracellular matrix (ECM) organization centered at LAMB2 and HTR1A, whereas the OCD risk genes were enriched in that of synapse functions. Researchers have devoted themselves to studying the dysfunction of neurotransmitter systems for both AN and OCD [68], with a primary focus on the neuronal and synaptic functions [69]; however, a study on the impact of ECM and non-neuronal cell types is still lacking. These studies have been proven as fruitful and provided profound insights into the mechanism of the two diseases. However, our results suggest that applying this identical strategy on two different diseases may introduce some kinds of bias since the risk genes of AN show significant impacts on ECM organization as well as on non-neuronal cell types, such as the endothelium and astrocytes, which is distinct from that of OCD. Large amounts of evidence [70] show that ECM modulates various functions, such as stem cell differentiation, neuron migration, axonal tract formation, myelination, and synapse maturation, all of which serve as essential processes of neurodevelopment. Similarly, non-neuronal cell types modulate neurotransmitter and synaptic functions via molecule clearance, neurotrophic release, and other processes [71]. Both our study and previous studies [72,73] suggest that such critical processes have an essential role in the pathology of AN and should not be neglected. We also infer that such a divergence with OCD may be one of the essential mechanisms underlying the phenotypic difference of the two diseases.
We also notice that the current study has limitations in the inclusion of risk genes. Various genetic factors other than SNP, including de novo mutation and copy number variation, have a fundamental contribution to psychiatric diseases, yet they are excluded in the current study. This is partly because limited trio-based sequencing studies are available publicly for AN or OCD [50,74]. The GWAS used for the current study did not excluded co-morbidity, which might introduce bias into the identified risk genes. Furthermore, the transcriptome data of ED came from various subtypes of ED, with only three samples coming from classic AN patients. Thus, the result of the ED cortex expression might contain signals shared by different types of ED, instead of being specific to AN.

5. Conclusions

The AN risk genes were enriched in coexpression modules that were involved in extracellular matrix functions and highly expressed in the postnatal brain, limbic system, and non-neuronal cell types, whereas the OCD risk genes were enriched in modules of synapse function, the prenatal brain, cortex layers, and neurons. However, they led to similar prefrontal cortex expression alterations, which influenced synapse transmission. We have found that the functional and spatiotemporal enrichment patterns the AN and OCD risk genes were not identical, but they brought about similar expression alterations, which may be one of the reasons why they had similar but not identical clinical phenotypes.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/genes12040491/s1, Figure S1 Modules of AN and OCD had different preference on developmental periods A&B: GO-BP analysis for prenatal AN module and postnatal OCD module. C&D: Period-specific expression of AN and OCD module genes. Table S1 All AN and OCD risk genes included in the current study. Table S2 Full GSEA result.

Author Contributions

Conceptualization, G.N.L.; methodology, software and formal analysis, W.S. and W.W.; investigation, G.N.L.; resources, S.Y.; data curation, S.Y.; writing—original draft preparation, W.S.; writing—review and editing, G.N.L.; visualization, W.S.; supervision, G.N.L.; project administration, G.N.L.; funding acquisition, G.N.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by grants from the National Natural Science Foundation of China (No. 81971292); Program for Professor of Special Appointment (Eastern Scholar) at Shanghai Institutions of Higher Learning (No. 1610000043); Medical-industrial intersection program of Shanghai Jiao Tong University (No. YG2019QNA59); Shanghai Mental Health Center (No. 2019-YJ01 and 2019-QH-03).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data analyzed in this study are curated from public domain.

Acknowledgments

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Trace, S.E.; Baker, J.H.; Peñas-Lledó, E.; Bulik, C.M. The genetics of eating disorders. Annu. Rev. Clin. Psychol. 2013, 9, 589–620. [Google Scholar] [CrossRef] [Green Version]
  2. Pauls, D.L.; Abramovitch, A.; Rauch, S.L.; Geller, D.A. Obsessive-compulsive disorder: An integrative genetic and neurobiological perspective. Nat. Publ. Group 2014, 15, 410–424. [Google Scholar] [CrossRef]
  3. Pettersson, E.; Lichtenstein, P.; Larsson, H.; Song, J.; Agrawal, A.; Borglum, A.D.; Bulik, C.M.; Daly, M.J.; Davis, L.K.; Demontis, D.; et al. Genetic influences on eight psychiatric disorders based on family data of 4 408 646 full and half-siblings, and genetic data of 333 748 cases and controls. Psychol. Med. 2019, 49, 1166–1173. [Google Scholar] [CrossRef] [Green Version]
  4. Cederlöf, M.; Thornton, L.M.; Baker, J.; Lichtenstein, P.; Larsson, H.; Rück, C.; Bulik, C.M.; Mataix-Cols, D. Etiological overlap between obsessive-compulsive disorder and anorexia nervosa: A longitudinal cohort, multigenerational family and twin study. World Psychiatry 2015, 14, 333–338. [Google Scholar] [CrossRef] [Green Version]
  5. Swinbourne, J.M.; Touyz, S.W. The co-morbidity of eating disorders and anxiety disorders: A review. Eur. Eat. Disord. Rev. 2007, 15, 253–274. [Google Scholar] [CrossRef]
  6. Micali, N.; Hilton, K.; Natatani, E.; Heyman, I.; Turner, C.; Mataix-Cols, D. Is childhood OCD a risk factor for eating disorders later in life? A longitudinal study. Psychol. Med. 2011, 41, 2507–2513. [Google Scholar] [CrossRef]
  7. Levinson, C.A.; Brosof, L.C.; Ram, S.S.; Pruitt, A.; Russell, S.; Lenze, E.J. Obsessions are strongly related to eating disorder symptoms in anorexia nervosa and atypical anorexia nervosa. Eat. Behav. 2019, 34. [Google Scholar] [CrossRef] [PubMed]
  8. Naylor, H.; Mountford, V.; Brown, G. Beliefs about excessive exercise in eating disorders: The role of obsessions and compulsions. Eur. Eat. Disord. Rev. 2011, 19, 226–236. [Google Scholar] [CrossRef] [PubMed]
  9. Sunday, S.R.; Halmi, K.A.; Einhorn, A. The Yale-Brown-Cornell Eating Disorder Scale: A new scale to assess eating disorder symptomatology. Int. J. Eat. Disord. 1995, 18, 237–245. [Google Scholar] [CrossRef]
  10. Levinson, C.A.; Sala, M.; Fewell, L.; Brosof, L.C.; Fournier, L.; Lenze, E.J. Meal and snack-time eating disorder cognitions predict eating disorder behaviors and vice versa in a treatment seeking sample: A mobile technology based ecological momentary assessment study. Behav. Res. Ther. 2018, 105, 36–42. [Google Scholar] [CrossRef]
  11. Altman, S.E.; Shankman, S.A. What is the association between obsessive-compulsive disorder and eating disorders? Clin. Psychol. Rev. 2009, 29, 638–646. [Google Scholar] [CrossRef]
  12. Strober, M.; Freeman, R.; Lampert, C.; Diamond, J. The association of anxiety disorders and obsessive compulsive personality disorder with anorexia nervosa: Evidence from a family study with discussion of nosological and neurodevelopmental implications. Int. J. Eat. Disord. 2007, 40. [Google Scholar] [CrossRef] [PubMed]
  13. Yilmaz, Z.; Halvorsen, M.; Bryois, J.; Yu, D.; Thornton, L.M.; Zerwas, S.; Micali, N.; Moessner, R.; Burton, C.L.; Zai, G.; et al. Examination of the shared genetic basis of anorexia nervosa and obsessive-compulsive disorder. Mol. Psychiatry 2018. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Anttila, V.; Bulik-Sullivan, B.; Finucane, H.K.; Walters, R.K.; Bras, J.; Duncan, L.; Escott-Price, V.; Falcone, G.J.; Gormley, P.; Malik, R.; et al. Analysis of shared heritability in common disorders of the brain. Science 2018, 360. [Google Scholar] [CrossRef] [Green Version]
  15. Cross-Disorder Group of the Psychiatric Genomics Consortium; Won, H.; Feng, Y.-C.A.; Rosenthal, J.; Zhu, Z.; Tucker-Drob, E.M.; Nivard, M.G.; Grotzinger, A.D.; Posthuma, D.; Wang, M.M.-J.; et al. Genomic relationships, novel loci, and pleiotropic mechanisms across eight psychiatric disorders. Cell 2019, 179, 1469–1482. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Boisseau, C.L.; Thompson-Brenner, H.; Caldwell-Harris, C.; Pratt, E.; Farchione, T.; Harrison Barlow, D. Behavioral and cognitive impulsivity in obsessive-compulsive disorder and eating disorders. Psychiatry Res. 2012, 200, 1062–1066. [Google Scholar] [CrossRef] [PubMed]
  17. Van Passel, B.; Danner, U.N.; Dingemans, A.E.; Aarts, E.; Sternheim, L.C.; Becker, E.S.; Van Elburg, A.A.; Van Furth, E.F.; Hendriks, G.J.; Cath, D.C. Cognitive remediation therapy does not enhance treatment effect in obsessive-compulsive disorder and anorexia nervosa: A randomized controlled trial. Psychother. Psychosom. 2020, 89, 228–241. [Google Scholar] [CrossRef]
  18. Fischer, M.S.; Baucom, D.H.; Baucom, B.R.; Abramowitz, J.S.; Kirby, J.S.; Bulik, C.M. Disorder-specific patterns of emotion coregulation in couples: Comparing obsessive compulsive disorder and anorexia nervosa. J. Fam. Psychol. 2017, 31, 304–315. [Google Scholar] [CrossRef]
  19. Schneier, F.R.; Kimeldorf, M.B.; Choo, T.H.; Steinglass, J.E.; Wall, M.M.; Fyer, A.J.; Simpson, H.B. Attention bias in adults with anorexia nervosa, obsessive-compulsive disorder, and social anxiety disorder. J. Psychiatr. Res. 2016, 79, 61–69. [Google Scholar] [CrossRef] [Green Version]
  20. Young, S.; Rhodes, P.; Touyz, S.; Hay, P. The relationship between obsessive-compulsive personality disorder traits, obsessive-compulsive disorder and excessive exercise in patients with anorexia nervosa: A systematic review. J. Eat. Disord. 2013, 1, 16. [Google Scholar] [CrossRef] [Green Version]
  21. Flamarique, I.; Plana, M.T.; Castro-Fornieles, J.; Borràs, R.; Moreno, E.; Lázaro, L. Comparison of perfectionism dimensions in adolescents with anorexia nervosa or obsessive-compulsive disorder. J. Can. Acad. Child Adolesc. Psychiatry 2019, 28, 45–54. [Google Scholar]
  22. Beszłej, J.A.; Wieczorek, T.; Kobyłko, A.; Piotrowski, P.; Siwicki, D.; Weiser, A.; Fila-Witecka, K.; Rymaszewska, J.; Tabakow, P. Deep brain stimulation: New possibilities for the treatment of mental disorders. Psychiatr. Pol. 2019, 53, 789–806. [Google Scholar] [CrossRef]
  23. Mayer, L.E.; Walsh, B.T. The use of selective serotonin reuptake inhibitors in eating disorders. J. Clin. Psychiatry 1998, 59, 28–34. [Google Scholar] [PubMed]
  24. Simpson, H.B.; Wetterneck, C.T.; Cahill, S.P.; Steinglass, J.E.; Franklin, M.E.; Leonard, R.C.; Weltzin, T.E.; Riemann, B.C. Treatment of obsessive-compulsive disorder complicated by comorbid eating disorders. Cogn. Behav. Ther. 2013, 42, 64–76. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Watson, H.J.; Yilmaz, Z.; Thornton, L.M.; Hübel, C.; Coleman, J.R.I.; Gaspar, H.A.; Bryois, J.; Hinney, A.; Leppä, V.M.; Mattheisen, M.; et al. Genome-wide association study identifies eight risk loci and implicates metabo-psychiatric origins for anorexia nervosa. Nat. Genet. 2019, 51, 1207–1214. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Arnold, P.D.; Askland, K.D.; Barlassina, C.; Bellodi, L.; Bienvenu, O.J.; Black, D.; Bloch, M.; Brentani, H.; Burton, C.L.; Camarena, B.; et al. Revealing the complex genetic architecture of obsessive-compulsive disorder using meta-analysis. Mol. Psychiatry 2018, 23, 1181–1188. [Google Scholar] [CrossRef]
  27. Miller, J.A.; Ding, S.-L.; Sunkin, S.M.; Smith, K.A.; Ng, L.; Szafer, A.; Ebbert, A.; Riley, Z.L.; Royall, J.J.; Aiona, K.; et al. Transcriptional landscape of the prenatal human brain. Nature 2014, 508, 199–206. [Google Scholar] [CrossRef]
  28. Jaffe, A.E.; Deep-Soboslay, A.; Tao, R.; Hauptman, D.T.; Kaye, W.H.; Arango, V.; Weinberger, D.R.; Hyde, T.M.; Kleinman, J.E. Genetic neuropathology of obsessive psychiatric syndromes. Transl. Psychiatry 2014, 4. [Google Scholar] [CrossRef] [Green Version]
  29. Langfelder, P.; Horvath, S. WGCNA: An R package for weighted correlation network analysis. BMC Bioinform. 2008, 9, 559. [Google Scholar] [CrossRef] [Green Version]
  30. Yu, G.; Wang, L.-G.; Han, Y.; He, Q.-Y. clusterProfiler: An R package for comparing biological themes among gene clusters. Omics 2012, 16, 284–287. [Google Scholar] [CrossRef]
  31. Skene, N.G.; Grant, S.G.N. Identification of vulnerable cell types in major brain disorders using single cell transcriptomes and expression weighted cell type enrichment. Front. Neurosci. 2016, 10, 16. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Lin, G.N.; Corominas, R.; Lemmens, I.; Yang, X.; Tavernier, J.; Hill, D.E.E.; Vidal, M.; Sebat, J.; Iakoucheva, L.M.M. Spatiotemporal 16p11.2 protein network implicates cortical late mid-fetal brain development and KCTD13-Cul3-RhoA pathway in psychiatric diseases. Neuron 2015, 85, 742–754. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Ritchie, M.E.; Phipson, B.; Wu, D.; Hu, Y.; Law, C.W.; Shi, W.; Smyth, G.K. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015, 43, e47. [Google Scholar] [CrossRef] [PubMed]
  34. Fabregat, A.; Jupe, S.; Matthews, L.; Sidiropoulos, K.; Gillespie, M.; Garapati, P.; Haw, R.; Jassal, B.; Korninger, F.; May, B.; et al. The reactome pathway knowledgebase. Nucleic Acids Res. 2018, 46, D649–D655. [Google Scholar] [CrossRef]
  35. Yu, G.; He, Q.Y. ReactomePA: An R/Bioconductor package for reactome pathway analysis and visualization. Mol. Biosyst. 2016, 12, 477–479. [Google Scholar] [CrossRef]
  36. Gandal, M.J.; Zhang, P.; Hadjimichael, E.; Walker, R.L.; Chen, C.; Liu, S.; Won, H.; van Bakel, H.; Varghese, M.; Wang, Y.; et al. Transcriptome-wide isoform-level dysregulation in ASD, schizophrenia, and bipolar disorder. Science 2018, 362, eaat8127. [Google Scholar] [CrossRef] [Green Version]
  37. Gallagher, M.D.; Chen-Plotkin, A.S. The post-GWAS era: From association to function. Am. J. Hum. Genet. 2018, 102, 717–730. [Google Scholar] [CrossRef] [Green Version]
  38. Lin, G.N.; Guo, S.; Tan, X.; Wang, W.; Qian, W.; Song, W.; Wang, J.; Yu, S.; Wang, Z.; Cui, D.; et al. PsyMuKB: An integrative de novo variant knowledge base for developmental disorders. Genom. Proteom. Bioinform. 2019. [Google Scholar] [CrossRef]
  39. Gandal, M.J.; Haney, J.R.; Parikshak, N.N.; Leppa, V.; Ramaswami, G.; Hartl, C.; Schork, A.J.; Appadurai, V.; Buil, A.; Werge, T.M.; et al. Shared molecular neuropathology across major psychiatric disorders parallels polygenic overlap. Science 2018, 359, 693–697. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Owen, M.J. New approaches to psychiatric diagnostic classification. Neuron 2014, 84, 564–571. [Google Scholar] [CrossRef] [Green Version]
  41. Shohat, S.; Ben-David, E.; Shifman, S. Varying intolerance of gene pathways to mutational classes explain genetic convergence across neuropsychiatric disorders. Cell Rep. 2017, 18, 2217–2227. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Li, M.; Santpere, G.; Kawasawa, Y.I.; Evgrafov, O.V.; Gulden, F.O.; Pochareddy, S.; Sunkin, S.M.; Li, Z.; Shin, Y.; Zhu, Y.; et al. Integrative functional genomic analysis of human brain development and neuropsychiatric risks. Science 2018, 362, eaat7615. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Nishioka, M.; Bundo, M.; Iwamoto, K.; Kato, T. Somatic mutations in the human brain: Implications for psychiatric research. Mol. Psychiatry 2019, 24, 839–856. [Google Scholar] [CrossRef] [Green Version]
  44. Crossley, N.A.; Scott, J.; Ellison-Wright, I.; Mechelli, A. Neuroimaging distinction between neurological and psychiatric disorders. Br. J. Psychiatry 2015, 207. [Google Scholar] [CrossRef] [Green Version]
  45. Zhan, X.; Yu, R. A window into the brain: Advances in psychiatric fMRI. BioMed Res. Int. 2015, 2015, 542467. [Google Scholar] [CrossRef] [Green Version]
  46. Bassett, D.S.; Xia, C.H.; Satterthwaite, T.D. Understanding the emergence of neuropsychiatric disorders with network neuroscience. Biol. Psychiatry Cogn. Neurosci. Neuroimaging 2018, 3, 742–753. [Google Scholar] [CrossRef]
  47. McCarthy, S.E.; Gillis, J.; Kramer, M.; Lihm, J.; Yoon, S.; Berstein, Y.; Mistry, M.; Pavlidis, P.; Solomon, R.; Ghiban, E.; et al. De novo mutations in schizophrenia implicate chromatin remodeling and support a genetic overlap with autism and intellectual disability. Mol. Psychiatry 2014, 19, 652–658. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. De Rubeis, S.; He, X.; Goldberg, A.P.; Poultney, C.S.; Samocha, K.; Cicek, A.E.; Kou, Y.; Liu, L.; Fromer, M.; Walker, S.; et al. Synaptic, transcriptional and chromatin genes disrupted in autism. Nature 2014, 515, 209–215. [Google Scholar] [CrossRef] [PubMed]
  49. Hamdan, F.F.; Srour, M.; Capo-Chichi, J.-M.; Daoud, H.; Nassif, C.; Patry, L.; Massicotte, C.; Ambalavanan, A.; Spiegelman, D.; Diallo, O.; et al. De novo mutations in moderate or severe intellectual disability. PLoS Genet. 2014, 10, e1004772. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Cappi, C.; Oliphant, M.E.; Péter, Z.; Zai, G.; do Rosário, M.C.; Sullivan, C.A.W.; Gupta, A.R.; Hoffman, E.J.; Virdee, M.; Olfson, E.; et al. De novo damaging DNA coding mutations are associated with obsessive-compulsive disorder and overlap with Tourette’s disorder and autism. Biol. Psychiatry 2019. [Google Scholar] [CrossRef]
  51. Huckins, L.M.; Hatzikotoulas, K.; Southam, L.; Thornton, L.M.; Steinberg, J.; Aguilera-Mckay, F.; Treasure, J.; Schmidt, U.; Gunasinghe, C.; Romero, A.; et al. Investigation of common, low-frequency and rare genome-wide variation in anorexia nervosa. Mol. Psychiatry 2018, 23, 1169–1180. [Google Scholar] [CrossRef] [Green Version]
  52. Lawrence, M.; Daujat, S.; Schneider, R. Lateral thinking: How histone modifications regulate gene expression. Trends Genet. 2016, 32, 42–56. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Thomson, C.A.; McColl, A.; Cavanagh, J.; Graham, G.J. Peripheral inflammation is associated with remote global gene expression changes in the brain. J. Neuroinflamm. 2014, 11. [Google Scholar] [CrossRef] [Green Version]
  54. McEwen, B.S.; Bowles, N.P.; Gray, J.D.; Hill, M.N.; Hunter, R.G.; Karatsoreos, I.N.; Nasca, C. Mechanisms of stress in the brain. Nat. Neurosci. 2015, 18, 1353–1363. [Google Scholar] [CrossRef] [PubMed]
  55. Carr, C.P.; Martins, C.M.S.; Stingel, A.M.; Lemgruber, V.B.; Juruena, M.F. The role of early life stress in adult psychiatric disorders: A systematic review according to childhood trauma subtypes. J. Nerv. Ment. Dis. 2013, 201, 1007–1020. [Google Scholar] [CrossRef]
  56. Bauer, M.E.; Teixeira, A.L. Inflammation in psychiatric disorders: What comes first? Ann. N. Y. Acad. Sci. 2019, 1437, 57–67. [Google Scholar] [CrossRef]
  57. Solmi, M.; Veronese, N.; Favaro, A.; Santonastaso, P.; Manzato, E.; Sergi, G.; Correll, C.U. Inflammatory cytokines and anorexia nervosa: A meta-analysis of cross-sectional and longitudinal studies. Psychoneuroendocrinology 2015, 51, 237–252. [Google Scholar] [CrossRef] [PubMed]
  58. Malik, R.; Chauhan, G.; Traylor, M.; Sargurupremraj, M.; Okada, Y.; Mishra, A.; Rutten-Jacobs, L.; Giese, A.K.; Van Der Laan, S.W.; Gretarsdottir, S.; et al. Multiancestry genome-wide association study of 520,000 subjects identifies 32 loci associated with stroke and stroke subtypes. Nat. Genet. 2018, 50, 524–537. [Google Scholar] [CrossRef] [Green Version]
  59. Lek, M.; Karczewski, K.J.; Minikel, E.V.; Samocha, K.E.; Banks, E.; Fennell, T.; O’Donnell-Luria, A.H.; Ware, J.S.; Hill, A.J.; Cummings, B.B.; et al. Analysis of protein-coding genetic variation in 60,706 humans. Nature 2016, 536, 285–291. [Google Scholar] [CrossRef] [Green Version]
  60. Huang, N.; Lee, I.; Marcotte, E.M.; Hurles, M.E. Characterising and predicting haploinsufficiency in the human genome. PLoS Genet. 2010, 6, e1001154. [Google Scholar] [CrossRef] [Green Version]
  61. Favaro, A.; Caregaro, L.; Tenconi, E.; Bosello, R.; Santonastaso, P. Time trends in age at onset of anorexia nervosa and bulimia nervosa. J. Clin. Psychiatry 2009, 70, 1715–1721. [Google Scholar] [CrossRef] [PubMed]
  62. Volpe, U.; Tortorella, A.; Manchia, M.; Monteleone, A.M.; Albert, U.; Monteleone, P. Eating disorders: What age at onset? Psychiatry Res. 2016, 238, 225–227. [Google Scholar] [CrossRef] [PubMed]
  63. Taylor, S. Early versus late onset obsessive–compulsive disorder: Evidence for distinct subtypes. Clin. Psychol. Rev. 2011, 31, 1083–1100. [Google Scholar] [CrossRef] [PubMed]
  64. Scharner, S.; Stengel, A. Alterations of brain structure and functions in anorexia nervosa. Clin. Nutr. Exp. 2019, 28, 22–32. [Google Scholar] [CrossRef] [Green Version]
  65. Li, B.; Mody, M. Cortico-striato-thalamo-cortical circuitry, working memory, and obsessive-compulsive disorder. Front. Psychiatry 2016, 7, 78. [Google Scholar] [CrossRef] [Green Version]
  66. Sinopoli, V.M.; Burton, C.L.; Kronenberg, S.; Arnold, P.D. A review of the role of serotonin system genes in obsessive-compulsive disorder. Neurosci. Biobehav. Rev. 2017, 80, 372–381. [Google Scholar] [CrossRef]
  67. Riva, G. Neurobiology of anorexia nervosa: Serotonin dysfunctions link self-starvation with body image disturbances through an impaired body memory. Front. Hum. Neurosci. 2016, 10. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Barbarich, N.C.; Kaye, W.H.; Jimerson, D. Neurotransmitter and imaging studies in anorexia nervosa: New targets for treatment. Curr. Drug Targets CNS Neurol. Disord. 2003, 2, 61–72. [Google Scholar] [CrossRef] [Green Version]
  69. Kaye, W.H.; Fudge, J.L.; Paulus, M. New insights into symptoms and neurocircuit function of anorexia nervosa. Nat. Rev. Neurosci. 2009, 10, 573–584. [Google Scholar] [CrossRef] [PubMed]
  70. Barros, C.S.; Franco, S.J.; Müller, U. Extracellular Matrix: Functions in the nervous system. Cold Spring Harb. Perspect. Biol. 2011, 3, 1–24. [Google Scholar] [CrossRef] [Green Version]
  71. Auld, D.S.; Robitaille, R. Glial cells and neurotransmission: An inclusive view of synaptic function. Neuron 2003, 40, 389–400. [Google Scholar] [CrossRef] [Green Version]
  72. Kesselmeier, M.; Pütter, C.; Volckmar, A.L.; Baurecht, H.; Grallert, H.; Illig, T.; Ismail, K.; Ollikainen, M.; Silén, Y.; Keski-Rahkonen, A.; et al. High-throughput DNA methylation analysis in anorexia nervosa confirms TNXB hypermethylation. World J. Biol. Psychiatry 2018, 19, 187–199. [Google Scholar] [CrossRef] [PubMed]
  73. Frintrop, L.; Trinh, S.; Liesbrock, J.; Leunissen, C.; Kempermann, J.; Etdöger, S.; Kas, M.J.; Tolba, R.; Heussen, N.; Neulen, J.; et al. The reduction of astrocytes and brain volume loss in anorexia nervosa—The impact of starvation and refeeding in a rodent model. Transl. Psychiatry 2019, 9. [Google Scholar] [CrossRef] [PubMed]
  74. Cappi, C.; Brentani, H.; Lima, L.; Sanders, S.J.; Zai, G.; Diniz, B.J.; Reis, V.N.S.; Hounie, A.G.; Conceição do Rosário, M.; Mariani, D.; et al. Whole-exome sequencing in obsessive-compulsive disorder identifies rare mutations in immunological and neurodevelopmental pathways. Transl. Psychiatry 2016, 6, e764. [Google Scholar] [CrossRef] [Green Version]
Genes 12 00491 g001 550
Figure 1. Coexpression and functional enrichments of anorexia nervosa (AN) and obsessive–compulsive disorder (OCD) genes. (A) WGCNA dendrogram for the prenatal brain (left) and postnatal brain (right) expression profiles from BrainSpan. (B): Enrichment of AN and OCD genes in the prenatal (left) and postnatal (right) coexpression modules. Asterisks denoted modules with the lowest p value for each disorder. (C,D) GO-BP enrichment results for the enriched modules. ECM: extracellular matrix organization.
Figure 1. Coexpression and functional enrichments of anorexia nervosa (AN) and obsessive–compulsive disorder (OCD) genes. (A) WGCNA dendrogram for the prenatal brain (left) and postnatal brain (right) expression profiles from BrainSpan. (B): Enrichment of AN and OCD genes in the prenatal (left) and postnatal (right) coexpression modules. Asterisks denoted modules with the lowest p value for each disorder. (C,D) GO-BP enrichment results for the enriched modules. ECM: extracellular matrix organization.
Genes 12 00491 g001
Genes 12 00491 g002 550
Figure 2. Spatiotemporal and cell-type-specific expression patterns for the enriched modules. (A,B) Spatial expression patterns for the OCD module (A, M6) and AN module (B, M2). The y-axis shows the module eigengene. (C) Brain region-specific expression patterns for M6 (left) and M2 (right). The x-axis shows the fold change calculated by the EWCE R package; the colors correspond to the p-values obtained from the EWCE. (D) Cell-type-specific expression patterns for M6 (left) and M2 (right), similar to C. MD: medium dorsal thalamus. HIP: Hippocampus. AMY: Amygdala. STR: Striatum. CB: Cerebellum. OC: Occipital cortex. PFC: Prefrontal cortex. TC: Temporal cortex. PC: Parietal cortex. OPC: oligodendrocyte precursor. Oligo: oligodendrocyte. Ex_neuron: excitatory neuron. In_neuron: inhibitory neuron.
Figure 2. Spatiotemporal and cell-type-specific expression patterns for the enriched modules. (A,B) Spatial expression patterns for the OCD module (A, M6) and AN module (B, M2). The y-axis shows the module eigengene. (C) Brain region-specific expression patterns for M6 (left) and M2 (right). The x-axis shows the fold change calculated by the EWCE R package; the colors correspond to the p-values obtained from the EWCE. (D) Cell-type-specific expression patterns for M6 (left) and M2 (right), similar to C. MD: medium dorsal thalamus. HIP: Hippocampus. AMY: Amygdala. STR: Striatum. CB: Cerebellum. OC: Occipital cortex. PFC: Prefrontal cortex. TC: Temporal cortex. PC: Parietal cortex. OPC: oligodendrocyte precursor. Oligo: oligodendrocyte. Ex_neuron: excitatory neuron. In_neuron: inhibitory neuron.
Genes 12 00491 g002
Genes 12 00491 g003 550
Figure 3. Expression alterations of the prefrontal cortex of ED and OCD patients. (A) Differential Expression Analysis for OCD vs. control (X-axis) and ED vs control (Y-axis). Each point corresponded to one gene; colors correspond to the minimal of the p-value of OCD vs. control and ED vs. control. (B) Venn diagram showing the relation between the differentially expressed genes of ED vs. control (left) and OCD vs. control (right). (C) Gene Set Enrichment Analysis (GSEA) results using gene lists ranked by fold change of OCD vs. control (X-axis) and AN vs. control (Y-axis). Each point corresponds to one REACTOME term; colors correspond to the minimal of the p-value of OCD vs. control and AN vs. control. (D) GSEA plots for top REACTOME terms. NES: Normalized enrichment score. logFC: log fold change. ED: Eating disorder.
Figure 3. Expression alterations of the prefrontal cortex of ED and OCD patients. (A) Differential Expression Analysis for OCD vs. control (X-axis) and ED vs control (Y-axis). Each point corresponded to one gene; colors correspond to the minimal of the p-value of OCD vs. control and ED vs. control. (B) Venn diagram showing the relation between the differentially expressed genes of ED vs. control (left) and OCD vs. control (right). (C) Gene Set Enrichment Analysis (GSEA) results using gene lists ranked by fold change of OCD vs. control (X-axis) and AN vs. control (Y-axis). Each point corresponds to one REACTOME term; colors correspond to the minimal of the p-value of OCD vs. control and AN vs. control. (D) GSEA plots for top REACTOME terms. NES: Normalized enrichment score. logFC: log fold change. ED: Eating disorder.
Genes 12 00491 g003
Genes 12 00491 g004 550
Figure 4. Network analysis for the AN and OCD modules. (A) Network for the AN module. (B) Network for the OCD module. For each plot, points corresponded to genes within the module, TOM: Topology Overlap Measures, a metric that reflect the coexpression strength between two genes. Only edges with a Topology Overlap Measure (TOM) larger than 0.1 are shown. Risk genes (proxy genes of GWAS top loci) for AN and OCD are shown in large print, with the gene symbols printed. ECM: extracellular matrix organization. ED: Eating disorder.
Figure 4. Network analysis for the AN and OCD modules. (A) Network for the AN module. (B) Network for the OCD module. For each plot, points corresponded to genes within the module, TOM: Topology Overlap Measures, a metric that reflect the coexpression strength between two genes. Only edges with a Topology Overlap Measure (TOM) larger than 0.1 are shown. Risk genes (proxy genes of GWAS top loci) for AN and OCD are shown in large print, with the gene symbols printed. ECM: extracellular matrix organization. ED: Eating disorder.
Genes 12 00491 g004
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Song, W.; Wang, W.; Yu, S.; Lin, G.N. Dissection of the Genetic Association between Anorexia Nervosa and Obsessive–Compulsive Disorder at the Network and Cellular Levels. Genes 2021, 12, 491. https://doi.org/10.3390/genes12040491

AMA Style

Song W, Wang W, Yu S, Lin GN. Dissection of the Genetic Association between Anorexia Nervosa and Obsessive–Compulsive Disorder at the Network and Cellular Levels. Genes. 2021; 12(4):491. https://doi.org/10.3390/genes12040491

Chicago/Turabian Style

Song, Weichen, Weidi Wang, Shunying Yu, and Guan Ning Lin. 2021. "Dissection of the Genetic Association between Anorexia Nervosa and Obsessive–Compulsive Disorder at the Network and Cellular Levels" Genes 12, no. 4: 491. https://doi.org/10.3390/genes12040491

APA Style

Song, W., Wang, W., Yu, S., & Lin, G. N. (2021). Dissection of the Genetic Association between Anorexia Nervosa and Obsessive–Compulsive Disorder at the Network and Cellular Levels. Genes, 12(4), 491. https://doi.org/10.3390/genes12040491

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop