https://doi.org/10.22037/ijpr.2018.2378

Protein-Protein Interaction Analysis of Common Top Genes in Obsessive-Compulsive Disorder (OCD) and Schizophrenia: Towards New Drug Approach Obsessive-Compulsive disorder (OCD) and Schizophrenia Comorbidity Gene Analysis

authors:

avatar Mohammad Mahboubi a , avatar Mona Zamanian Azodi b , avatar Mostafa Rezaei Tavirani b , avatar Vahid Mansouri b , avatar Nayeb Ali Ahmadi b , avatar Mostafa Hamdieh c , avatar Majid Rezaei Tavirani d , * , avatar Bahar Naghavi Gargari e

Abadan School of Medical Sciences, Abadan, Iran.
Proteomics Research Center, Faculty of Paramedical Sciences, Shahid Beheshti University of Medical Sciences, Tehran, Iran.
Psychosomatic Department, Taleghani Hospital, Faculty of Medicine, Shahid Beheshti University of Medical Sciences, Tehran, Iran.
Surgery department, Faculty of Medicine, Iran University of Medical Sciences, Tehran, Iran.
Basic Science Department, Shahid Beheshti University of Medical Sciences, Tehran, Iran.
Iranian Journal of Pharmaceutical Research : IJPR: Vol.17, issue Suppl2; 173-186
article type: Original Article
received: June , 2018
accepted: December , 2018

How To Cite Mahboubi M, Zamanian Azodi M, Rezaei Tavirani M, Mansouri V, Ali Ahmadi N, et al. Protein-Protein Interaction Analysis of Common Top Genes in Obsessive-Compulsive Disorder (OCD) and Schizophrenia: Towards New Drug Approach Obsessive-Compulsive disorder (OCD) and Schizophrenia Comorbidity Gene Analysis. Iran J Pharm Res. 2018;17(Suppl2):e126092. https://doi.org/10.22037/ijpr.2018.2378.

Abstract

Comorbidty is common among psychiatric disorders including obsessive-compulsive disorder and schizophrenia with a high rate. Many studies suggested that the disorders may have same etiological bases. In this regard, shared pathways of glutamate, dopaminergic, and serotonin are the known ones. Here, the common significant genes are examined to understand the possible molecular origin of the disorders in terms of sequence and functional features. Exploring the underling mechanisms of OCD and schizophrenia is important to achieve a better treatment options. Methods of Cytoscape software following R statistical software were applied for this purpose. Needleman-Wunsch global alignment algorithm was used to determine pair-wise similarities followed by clustering methods, AGNES and PAM in R statistical programming software. The results indicate that SLC1A1, DRD2, DRD4, BDNF, ESR1, CDH2, GRIN2B, TNFa, GABBR1, and OLIG2 are significantly common for the two disorders and PPI network analysis showed the important key genes in the interaction profile. ESR1 (estrogen receptor α) as a key hub-bottleneck gene regulates many underling mechanisms of the brain. Application of global alignments indicates some of the genes with sequence similarities also elucidate similar biological terms.

Introduction

Psychiatric disorders affect people around the world noticeably. Some of them tend to manifest together known as comorbidity conditions (2). OCD and schizophrenia as widespread mental disorders are one of the important co-occurrence types (3). OCD is identified by intrusive thoughts and the urge for reducing the unpleasant experience by different types of compulsions (4). The patients suffer from recurrent invasive thoughts or stereotype behaviors at least one hour per day (5). It is with the prevalence of about 1-3% in human population (6, 7). On the other hand, schizophrenia is a sever condition of consciousness abnormalities including hallucinations, delusions, disorganized communication, poor planning, reduced motivation, and blunted affect with 0.5% and 1% lifetime worldwide prevalence (8, 9). Both disorders have demonstrated heterogeneity in genetic concept and clinical phenotype (10, 11). About 30.3% of patients diagnosed with schizophrenia have the obsessive-compulsive symptoms (OCS) and about 13.6% of them shows co-occurrence with OCD (12).  In fact, OCD can be accounted as schizophrenia risk factor (3). There are some common features between OCD and schizophrenia despite the different diagnostic entities. Some symptoms, age of onset, gender, and structural impairment in brain including areas of prefrontal cortex, anterior cingulate, caudate nucleus and thalamus are shared in OCD and schizophrenia. The patients with two disorders, schizo-OCD, have more difficulties in prognosis, treatment responses, and lower quality of life (13). In addition, some of the applied therapeutic methods are also similar (14). Therefore, the disorders may have same molecular origin as they have a tendency of comorbidity (15). Some pathways of serotonin, dopamine and glutamine are reported common between two disorders. The most applied medications for schizophrenia-OCD comorbidity belong to selective serotonin re-uptake inhibitors (SSRI) family (16, 17). It must be noted that some of patients were not responded to treatment. In addition, some other types of drugs affecting on dopamine receptors (antagonism) also showed effective in patients with OCD (18). Antipsychotic medications, which substantially act as dopamine antagonism, are also effective in patients with schizophrenia (19). Particularly, Clozapine and Olanzapine as dopamine antagonism seemed valuable with or without SSRIs application in some schizophrenia patients (20). Both drugs act as D2 receptors antagonisms in the mesolimbic pathway; in addition, they have high affinity to 5-HT2A receptors in the frontal cortex as well (21). However, Clozapine and Olanzapine can increase the risk of obsessive-compulsive symptoms in schizophrenia patients (22). Proteomics, genomics, and protein-protein interaction (1) network analysis can provide a new insight of interactome profile of the two disorders (23, 24). Central genes in OCD comorbid schizophrenia PPI network are essentials in the network integrity and can be targets for more analysis (25). Since there are crucial genes among query and added genes which play significant role in disorder PPI network analysis is a useful tool to explore the critical genes. It is essential to evaluate the structural and functional properties of the related significant genes to decipher the core underling mechanisms of pathogenesis of the OCD and schizophrenia. Furthermore, the analysis of essential genes in the network can provide a new insight in drug targeting strategies.

Experimental

The associated genes with development of OCD and schizophrenia were searched from Google Scholar and PubMed. There are many reported candidate genes for the two types of disorders; however, only limited numbers were significantly common between two disorders. In fact, based on previous investigation, about 39 genes are susceptible to pathogenicity of OCD which among them only 24 candidates show significant relationship with the disorder. The individually emphasized genes for OCD and schizophrenia in the literature were selected and used for more evaluation. About ten polymorphism genes showed significant associations. The Uniprot accession codes were downloaded from http://www.uniprot.org/. Further analysis including PPI network study and sequence global comparison was handled by the applied codes. PPI network was constructed by the use of PSICQUIC (Proteomics Standard Initiative Common QUery InterfaCe) in Cytoscape v: 3.3.0., from public database including Mentha, Reactome-Fls, and Reactome. This algorithm is a universal web service client for importing interactions from public databases (26). Network Analyzer provides topological features including two used in this study, degree and betweenness centralities. This plug-in is integrated into Cytoscape and can provide different types of centrality properties of the analyzed network (27). In the PPI network, the nodes represent genes, while the edges represent physical interactions (28). The number of connected edges to the nodes indicates the degree value of that element. In addition, another evaluated centrality in the network, betweenness centrality signifies the number of shortest path that pass through each node (27). Genes with high betweenness centrality are known as bottlenecks, which can be as essential as hubs. Bottlenecks can connect many protein complexes and processes (29). For gene ontology examination, ClueGO+CluePedia plug-ins in Cytoscape can be helpful. The application explored the biological terms of the investigated genes (30, 31). GO is the reference of controlled vocabulary of the terms of three biological annotations including cellular component, molecular function, and biological process (1). There are some available related statistical features in ClueGo+CluePedia plug-ins. Kappa score was the statistical method for this analysis. The higher the score, the better the terms are grouped significantly together. Kappa score was set to 0.5. Furthermore, other setting was as follows: number of genes per term 2, percent of related genes in each term, 1 and 2 and the minimum and maximum level of ontology was set to 3 to 8 as default options, respectively. The P-value was also set to ≤ 0.05. The correction method for p-value was Bonferroni step down method and the enrichment/depletion test for terms was set to Two-sided (enrichment/depletion) based on hypergeometric. The assessment of sequence similarity of the examined genes with R programming software was another part of the conducted research. By the use of R, clustering profile of the ten genes can be understood. R software is free statistical software available on (www.r-project.org). It possess various packages, one of them is the package cluster can be used to identify subgraphs with maximal density (32). Protein sequences were extracted from Uniprot in FASTA format and then the global pairwise alignment was performed using Needleman-Wunsch algorithm, available in EMBOSS Needle, EBI Database. The construction of similarity matrix by the alignment similarity scores is provided by the alignment. The scores were then deducted from 1, as the dissimilarity matrix is need for clustering purposes in R studio software (Dissimilarity =1-similarity) (33). Agglomerative hierarchical clustering analysis (bottom up) as hierarchical and partitioning around medoids (PAM) non-hierarchical clustering algorithms were the conducted clustering methods. The distance between each cluster is calculated as D = 1-C, where D = distance and C = correlation between spot clusters. Average linkage (UPGMA) is the used method to define the distance between clusters as the average distance. Moreover, agglomerative coefficient signifies clustering structures strength and quality of the dataset. The calculating range is 0 to 1. The agglomerative coefficient (AC) close to 1, shows on good structured cluster, whereas AC close to 0 implies on not a powerful clustering structure (34). PAM clustering finds K representative objects called “medoids” from the dataset. In other words, data structures are represented by medoids as the representatives of the clusters (35). AGNES and PAM were the functions of clustering algorithms in this study. Dendrogram and silhouette plots are the graphical representations of clustering in this study for AGNES and PAM clustering algorithms, respectively. Peter J. Rousseeuw in 1986 presented a new graphical display called “silhouette”. This type of plot shows that which elements belong to their cluster strongly and which are in between clusters. The average silhouette width is obtained by mean of each cluster silhouette width value, which is from 0 to1 and implies on the clustering structure strength. In addition, it is useful for selecting the K (numbers of clusters) in silhouette (36). Validation of clustering algorithms was handled by different indices.

Results

Genes with significant correlation with OCD and schizophrenia are gathered from literature resources (Google Scholar and PubMed) (see Table 1).

Public databases including Mentha, Reactome-Fls, and Reactome were sources for PPI network construction of ten genes by the use of PSICQUIC (Proteomics Standard Initiative Common Query InterfaCe) algorithm in Cytoscape v: 3.3.0. (See Figure 1). As it is depicted in Figure 1 1555 genes are added to the 10 query genes. It is possible that the genes be screened to find the crucial individuals. The degree distribution of nodes which indicates the scale free network is illustrated in Figure 2.

The centralities evaluation provides noteworthy properties of the network topology. The genes with high value are essential elements for integrity of the network. About 20% of nodes with high degree are valuable elements and indicate hubs; similarly, 20% genes with high BC distribution imply on bottleneck genes. On the other hand, in interaction network, degree is a much better predictor of essentiality (29) (see Table 2).

For more resolution, cellular component, molecular function, biological process (BP), and pathway examination of ten genes were analyzed by ClueGo-Cluepedia Plug-ins (see Figures 3-6).

Estrogen receptor α (as a hub-bottleneck protein characterized with the most values of the two centrality indices that are tabulated in Table 2) is in close relationship with progesterone receptor (PGR) in the studied network. For better understanding the correlation is examined by CluePedia application of Cytoscape. The findings indicate the binding relationship between ESR1 and PGR and also activation of PGR by ESR1. Non-partitioning and partitioning clustering algorithms for sequence alignment require dissimilarity matrix. The dissimilarity matrix is calculated as (1-similarity) in which the scores are between 0 and 1. The lower the scores the shorter the distance would be between elements. The similarity and dissimilarity of the matrices of ten common candidate genes obtained from pair-wise sequence global alignment using Needleman-Wunsch algorithm from EBI databases and related dendrogram are shown in and Figure 7.

PAM clustering algorithm is done by the use of silhouette depiction. In each distinct cluster, the genes with the smallest dissimilarity average value are grouped (see Figure 8).

Table 1

List of 10 common genes in OCD and Schizophrenia with their Uniprot accession codes (AC) and related protein names is presented

RGene NamesProtein Names
1SLC1A1Excitatory amino acid transporter 3
2DRD2D(2) dopamine receptor
3DRD4D(4) dopamine receptor
4BDNFBrain-derived neurotrophic factor
5ESR1(37, 38)Estrogen receptor
6CDH2Cadherin-2
7GRIN2BGlutamate receptor ionotropic, NMDA 2B
8TNFaTumor necrosis factor
9GABBR1Gamma-aminobutyric acid type B receptor subunit 1
10OLIG2Oligodendrocyte transcription factor 2
Table 2

Hub-bottleneck nodes with high centrality properties (degree and betweenness) are the determinant of the essentiality of genes in the interaction network. The nodes are descended based on degree value

Hub-bottleneck nodesDegreeBetweenness Centrality
ESR18450.7
DRD22570.1
GABBR12240.1
TNFα2160.2
CDH21960.2
Protein-protein interaction network of ten common genes in OCD and Schizophrenia is shown. The PPI network is comprised of 1565 nodes and 2133 edges using Cytoscape 3.3.0. The yellow colored nodes are the ten query genes and the added genes are colored in green

Protein-protein interaction network of ten common genes in OCD and Schizophrenia is shown. The PPI network is comprised of 1565 nodes and 2133 edges using Cytoscape 3.3.0. The yellow colored nodes are the ten query genes and the added genes are colored in green

Figure 1

A scale-free network. The degree distribution values are significantly inhomogeneous. Just a few nodes show well linked, whereas others possess a small numbers of connections. This distribution implies on the presents of genes with high centrality values computed by Network Analyzer. The red line indicates the power law. The R-squared value is computed on logarithmized values, which is equal to 0.6 and the correlation = 0.9. Genes with high degree are in the right down region of the plot (their location is out of linear range)

A scale-free network. The degree distribution values are significantly inhomogeneous. Just a few nodes show well linked, whereas others possess a small numbers of connections. This distribution implies on the presents of genes with high centrality values computed by Network Analyzer. The red line indicates the power law. The R-squared value is computed on logarithmized values, which is equal to 0.6 and the correlation = 0.9. Genes with high degree are in the right down region of the plot (their location is out of linear range)

Figure 2

Biological process (BP) analysis of ten genes including their inhibition (red), activation (green), expression (yellow) characteristics by ClueGo-Cluepedia Plug-ins. The algorithm considers all the ten genes. Kappa score ≥0.5, P-value ≤ 0.05, gene per term = 2, gene percentage for each term = 2%

Biological process (BP) analysis of ten genes including their inhibition (red), activation (green), expression (yellow) characteristics by ClueGo-Cluepedia Plug-ins. The algorithm considers all the ten genes. Kappa score ≥0.5, P-value ≤ 0.05, gene per term = 2, gene percentage for each term = 2%

Figure 3

Cellular component (4) analysis of the ten genes by ClueGo-Cluepedia Plug-ins. The evaluated algorithm considered only five out of ten genes. Kappa score ≥ 0.5, P-value ≤ 0.05, gene per term = 2, gene percentage for each term =1%.

Cellular component (4) analysis of the ten genes by ClueGo-Cluepedia Plug-ins. The evaluated algorithm considered only five out of ten genes. Kappa score ≥ 0.5, P-value ≤ 0.05, gene per term = 2, gene percentage for each term =1%.

Figure 4

Molecular function (MF) analysis of ten investigated genes by ClueGo-Cluepedia Plug-ins. However, the evaluated algorithm considered only six out of ten genes. Kappa score≥ 0.5, P-value ≤0.05 , gene per term = 2, gene percentage for each term =1%.

Molecular function (MF) analysis of ten investigated genes by ClueGo-Cluepedia Plug-ins. However, the evaluated algorithm considered only six out of ten genes. Kappa score≥ 0.5, P-value ≤0.05 , gene per term = 2, gene percentage for each term =1%.

Figure 5

Pathway analysis of ten candidate genes by ClueGo-Cluepedia Plug-ins. KEGG Pathway and WIKIPathways were considered for pathway analysis. The evaluated algorithm considered only eight out of the ten genes. Kappa score≥0.5, P-value0.05 ≥ , gene per term = 2, gene percentage for each term =1%.

Pathway analysis of ten candidate genes by ClueGo-Cluepedia Plug-ins. KEGG Pathway and WIKIPathways were considered for pathway analysis. The evaluated algorithm considered only eight out of the ten genes. Kappa score≥0.5, P-value0.05 ≥ , gene per term = 2, gene percentage for each term =1%.

Figure 6

Dendrogram plotting of ten analyzed genes by the use of AGNES function. The height indicates the distance between the nodes. The agglomerative coefficient is 0.24

Dendrogram plotting of ten analyzed genes by the use of AGNES function. The height indicates the distance between the nodes. The agglomerative coefficient is 0.24

Figure 7

PAM clustering algorithms of ten genes via silhouette plotting, K = 4, Average silhouette width for each cluster is shown in the right side of the right column in the Figure. The overall average silhouette width is equal to 0.19

PAM clustering algorithms of ten genes via silhouette plotting, K = 4, Average silhouette width for each cluster is shown in the right side of the right column in the Figure. The overall average silhouette width is equal to 0.19

Figure 8

Plot of the APN, AD, ADM and FOM measures with number of cluster. The APN and AD, measures show homogenous decreased trend as the number of clusters increases for hierarchical clustering in which with six clusters has the best overall scores. For ADM, the trend decreases from two to three clusters, then there is increment to four clusters, and finally the score decreases as the number of clusters increases. In FOM, the value decreases from two to three clusters but increases afterwards to five clusters and subsequently decreased again to six clusters. Overall, hierarchical clustering has the lowest value of all evaluated indices

Plot of the APN, AD, ADM and FOM measures with number of cluster. The APN and AD, measures show homogenous decreased trend as the number of clusters increases for hierarchical clustering in which with six clusters has the best overall scores. For ADM, the trend decreases from two to three clusters, then there is increment to four clusters, and finally the score decreases as the number of clusters increases. In FOM, the value decreases from two to three clusters but increases afterwards to five clusters and subsequently decreased again to six clusters. Overall, hierarchical clustering has the lowest value of all evaluated indices

Figure 9

Dunn and Silhouette indices and also stability measure are used for assessment of validity of clustering method. Validation of clustering method via Dunn and Silhouette indices revealed optimal scores for indices with the number of clusters. Hierarchical clustering with four (score = 0.19) and six (score = 0.97) clusters performs the best in each case of Silhouette and Dunn, respectively.

The stability measures include: The average proportion of non-overlap (APN), the average distance (AD), the average distance between means (ADM), and the Figure of merit (FOM). The values should be minimized for each indices (39). Validation of clustering method via stability measure by considering APN, AD, ADM, and FOM is down. Optimal scores are provided for indices with the number of clusters. Hierarchical clustering provides the best scores for all the evaluating indices (APN, AD, and ADM) with six clusters (scores 0.07, 0.46, and 0.15 respectively) except for FOM index, which is best by three cluster (score = 0.30).

In Figure 9, the graphically visualization provides in detail behavior of each indices by the given cluster numbers.

Parallel and multistep analysis in this research showed that SLC1A1, DRD2, DRD4, BDNF, ESR1, CDH2, GRIN2B, TNFa, GABBR1, and OLIG2 are common genes between schizophrenia and OCD which ESR1, DRD2, GABBR1, TNFa, and CDH2 are the most important individuals. It was confirmed that ESR1 acts as a super critical genes in this regard.

Discussion

Epidemiological studies indicated that OCD has frequent incident in schizophrenia (40). Many factors may correlate with this comorbidity. Still, the relationship between disorders remains elusive and it requires more investigations. In this light, it is important to examine the underlying etiological mechanisms to better understanding of the common molecular bases. Earlier investigations highlighted the role of common pathways of (41). Here, by exploring common significant genes between two disorders from literature, it is tried to provide a new insight of other etiological molecular concepts of the disorders. Moreover, by evaluating molecular origins including sequential and functional properties, it may be possible to achieve a better treatment protocol for the comorbid disorders. As mentioned earlier, patients with this comorbidity are more impaired in their daily life and show treatment resistance features (13). Therefore, the choice of treatment is imperative to solve this dilemma. In this regard, ten genes that showed significant association with OCD and schizophrenia were gathered from literature. At first, about thirty-nine genes for OCD were found. Among them, only twenty-four genes demonstrated significant comorbidity based on highlighting by conducted studies. For schizophrenia, there are many reported genes in schizophrenia database (42). Yet, only common OCD twenty four significant genes were evaluated. At last, ten genes (see Table 1) concluded for this study, as they were significant in schizophrenia as well. As depicted in Figure 1, Cytoscape and its plug-ins, analyzed the common significant genes and a network with 1565 nodes and 2133 edges was constructed. The PPI network shows good integrity and all the genes are in connection with each other. As showed in Table 2, centrality analysis of the network indicates that ESR1, DRD2, GABBR1, TNFα, and CDH2 possess high values in degree and betweenness centralities. In other words, these elements are Hub-bottleneck nodes are those with high value in degree and betweenness centralities. As mentioned earlier, those with high amount of degree are assigned as hubs and those with high betweenness centrality values known as bottlenecks of the network (43). These genes are critical points in the network integrity, especially, ESR1 with highest score. It is also important gene in contamination/cleaning subtype of OCD in women (44). ESR1 in schizophrenia has lower expressions in both related wild type and variants. ESR1 is involved in cognative, emotional behavior, as it is active in limbic related areas of the brain. In addition, estrogen and its receptor are vital neurosteroid for regulating many brain processes (37). One of important task of estrogen is regulating serotonin, dopamine, and norepinephrine pathways (45). ESR1 is also a DNA-binding transcription factor in cell that once bound to estrogen; it translocates to the nuclear and therefore regulates the gene expression. It is also confirmed by the high bottleneck value of ESR1 in this study that is presented in Table 2.

Notable biological processes that are pertinent to the function of ESR1and the other common genes, are shown in Figure3. The processes of regulation of oxidoreductase activity Group and phospholipase C-activating G-protein coupled receptor signaling pathway. TNFα is involved in cross-talk between positive regulation of peptidyl-serine phosphorylation, negative regulation of adenylate cyclase activity, and regulation of oxidoreductase activity. These findings are in correlation with bottleneck analysis. The bottleneck score for ESR1 and TNFα was the highest equal to 0.7 and 0.2, respectively. CDH2 is also good linker between negative regulation of adenylate cyclase activity and regulation of gliogenesis Group. BDNF gene with 0.04 betweeness score and degree of 46 is a nonhub-bottleneck node with high interconnecting of central processes including gamma-aminobutyric acid signaling pathway, negative regulation of adenylate cyclase activity and positive regulation of peptidyl-serine phosphorylation. However, the rest of important biological processes facilitated with these genes, in particular, ESR1 with highest bottleneck score are not mentioned. This fact is due to not including other network related genes in the (BP) analysis. What′s more, the BP analysis suggests that the most significant biological process in this relationship is diencephalon development. As it is shown in Figure 4, 5, and 6, neuron projection terminus, dopamine neurotransmitter receptor activity, coupled via Gi/Go and Cocaine addiction are important CC, MF and pathway in the pathogenesis of the comorbidity, respectively. As it is found in results, ESR1 has a fundamental role in activation of progesterone receptor (PGR). The genes are also in close interaction in the network. Polymorphism in ESR1 may lead to malfunction of PGR. On the other hand, studies showed that menstrual cycle is one of the important phases that OCD may trigger (46). ESR1 representation in our research may suggest the decrease of progesterone and its receptor function in this phase; it may correlate to OCD manifestation. Further analysis was the clustering analysis of sequence of the studied genes. In the clustering, it is appear that the use of AGNES and PAM was to compare their action with each other. As it is depicted in Figures 7 and 8, the cluster analysis led to similar results. However, AGNES demonstrated a more powerful cluster as calculated by cluster validation methods. ESR1 and GABBR1 are clustered together and belong to the same pathway named Estrogen signaling pathway. In addition, DRD2 and DRD4 have similarity in their sequences and they tend to contribute to similar CC, BP, MF, and pathway. However, TNF, SLC1A1, and OLIG2 genes included in a distinct cluster did not conclude to a similar corresponding biological processes. In fact, their similarity in their sequence is not very significant as the agglomerative and average silhouette width demonstrates. As mentioned earlier, general prescription for OCD patients with schizophrenia is the application of (SSRIs). Monotherapy in patients with comorbid OCD- schizophrenia may not be adequate. In fact, there are some cases with resistance responses to this type of treatment. This fact implies of the requirement of other drug affecting on central elements of the co-morbid disorders. Here, was shown that some of the key genes in the network can play fundamental role since their impairment can impose vast dysfunction in the interaction network. Additionally, among the significant associated genes, some shows higher critical role in the reactome system. ESR1 fundamental role in the network integrity suggests many clinical researches. Here, it can be concluded that in women with OCD onset in menopausal phase, the hormone therapy (estrone) or selective estrogen receptor modulators (SERM) with a tissue-specific estrogen receptor agonist/antagonist activity, may reduce the OCD symptoms. Likewise, hormone therapy and SERM proved to be effective for schizophrenia patients for both women and men (47). Indeed, estrogen deficiency and ESR1 dysfunction can increase OCD symptoms, and in menopausal state, the amount of estrogen decreases noticeably. Furthermore, many symptoms of menopause alleviates after the mentioned treatments. Yet, hormone therapy can have many side effects such as increasing the risk of breast and ovarian cancers (48, 49). Therefore, (SERMs) can be a better alternative. Based on this evidence, sources of phytoestrogens such as soy beans and flax seeds also showed promising in treatment of menopause symptoms (50, 51). Therefore, they may be useful for reducing OCD symptoms. However, it requires fundamental clinical evaluations. In addition, progesterone receptor dysfunction in patients with ESR1 polymorphism may contribute to OCD onset in the stage of menstrual activity. Therefore, the essential function of ESR1 brings out another relationship between schizophrenia and OCD in females. In other words, gender may play indispensable role in the co-occurrence of schizophrenia and OCD. Our findings were highlighted new molecular aspects of OCD-schizophrenia comorbidity.

Conclusion

In conclusion, regulation of ESR1 may be critical point in schizophrenia patients with OCD incident. It is suggested that the drugs which target ESR1 can be considered beside SSRIs in treatment of schizophrenia-OCD comorbidity or both disorders individually.

Acknowledgements

References

  • 1.

    Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G. Gene Ontology: tool for the unification of biology. Nat. genet. 2000;25:25-9. [PubMed ID: 10802651].

  • 2.

    Khan AA, Jacobson KC, Gardner CO, Prescott CA, Kendler KS. Personality and comorbidity of common psychiatric disorders. Br. J. Psychiatry. 2005;186:190-6. [PubMed ID: 15738498].

  • 3.

    Meier SM, Petersen L, Pedersen MG, Arendt MC, Nielsen PR, Mattheisen M, Mors O, Mortonsen PB. Obsessive-compulsive disorder as a risk factor for schizophrenia: a nationwide study. JAMA. Psychiatry. 2014;71:1215-21. [PubMed ID: 25188738].

  • 4.

    Crino R, Slade T, Andrews G. The changing prevalence and severity of obsessive-compulsive disorder criteria from DSM-III to DSM-IV. Am. J. Psychiatry. 2005;162:876-82. [PubMed ID: 15863788].

  • 5.

    Pittenger C, Kelmendi B, Bloch M, Krystal JH, Coric V. Clinical treatment of obsessive compulsive disorder. Psychiatry (Edgmont). 2005;2:34.

  • 6.

    Tang R, Noh HJ, Wang D, Sigurdsson S, Swofford R, Perloski M, Duxbury M, Patterson EE, Albright J, Castelhano M, Auton A, Boyko AR, Feng G, Lindblad-Toh K, Karlsson EK. Candidate genes and functional noncoding variants identified in a canine model of obsessive-compulsive disorder. Genome. biol. 2014;15:R25. [PubMed ID: 24995881].

  • 7.

    Zamanian-Azodi M, Rezaei-Tavirani M, Kermani-Ranjbar T, Arefi-Oskouie A, Rezaei tavirani M, RahmatiRad S, SotodehAsl N. Pathophysiology, genetics, types, and treatments in obsessive compulsive disorder. Koomesh. 2015;16:475-87.

  • 8.

    Sass LA, Parnas J. Schizophrenia, Consciousness, and the Self. Schizophr Bull. 2003;29:427-44. [PubMed ID: 14609238].

  • 9.

    Saha S, Chant D, Welham J, McGrath J. A systematic review of the prevalence of schizophrenia. PLoS Med. 2005;2:413-33.

  • 10.

    Pauls DL. The genetics of obsessive-compulsive disorder: a review. Dialogues Clin. Neurosci. 2010;12:149-63. [PubMed ID: 20623920].

  • 11.

    Tsutsumi A, Glatt SJ, Kanazawa T, Kawashige S, Uenishi H, Hokyo A, Kaneko T, Moritani M, Kikuyama H, Koh J, Matsumura H, Yoneda H. The genetic validation of heterogeneity in schizophrenia. Behav. Brain Funct. 2011;7:43. [PubMed ID: 21981786].

  • 12.

    Swets M, Dekker J, van Emmerik-van Oortmerssen K, Smid GE, Smit F, de Haan L, Schoevers RA. The obsessive compulsive spectrum in schizophrenia, a meta-analysis and meta-regression exploring prevalence rates. Schizophr. Res. 2014;152:458-68. [PubMed ID: 24361303].

  • 13.

    Poyurovsky M, Faragian S, Pashinian A, Levi A, Viosburd A, Stryjer R, Weizman R, Fuchs C, Weizman A. Neurological soft signs in schizophrenia patients with obsessive-compulsive disorder. J Neuropsychiatry Clin. Neurosci. 2007;19:145-50. [PubMed ID: 17431060].

  • 14.

    Poyurovsky M, Koran LM. Obsessive–compulsive disorder (OCD) with schizotypy vs schizophrenia with OCD: diagnostic dilemmas and therapeutic implications. J. Psychiatry Res. 2005;39:399-408.

  • 15.

    Devi S, Rao NP, Badamath S, Chandrashekhar CR, Janardhan Reddy YC. Prevalence and clinical correlates of obsessive–compulsive disorder in schizophrenia. Compr. Psychiatry. 2015;56:141-8. [PubMed ID: 25308405].

  • 16.

    Kohl S, Ebbert L, Sartorius A, Kuhn J. A step forward in elucidating the mystery of OCD. Eur. Arch. Psychiatry Clin. Neurosci. 2015;265:735-6. [PubMed ID: 25708456].

  • 17.

    Buoli M, Serati M, Ciappolino V, Altamura AC. May selective serotonin reuptake inhibitors (SSRIs) provide some benefit for the treatment of schizophrenia? Expert. Opin. Pharmacother. 2016;17:1375-85.

  • 18.

    Koo M-S, Kim E-J, Roh D, Kim C-H. Role of dopamine in the pathophysiology and treatment of obsessive–compulsive disorder. Expert Rev. Neurother. 2010;10:275-90. [PubMed ID: 20136383].

  • 19.

    Faizi M, Salimi A, Rasoulzadeh M, Naserzadeh P, Pourahmad J. Schizophrenia induces oxidative stress and cytochrome C release in isolated rat brain mitochondria: a possible pathway for induction of apoptosis and neurodegeneration. Iran. J. Pharm .Res. 2014;13:93-100. [PubMed ID: 24711834].

  • 20.

    Chakos M, Lieberman J, Hoffman E, Bradford D, Sheitman B. Effectiveness of second-generation antipsychotics in patients with treatment-resistant schizophrenia: a review and meta-analysis of randomized trials. Am. J. Psychiatry. 2001;158:518-26. [PubMed ID: 11282684].

  • 21.

    Trichard C, Paillére-Martinot M-L, Attar-Levy D, Recassens C, Monnet F, Martinot J-L. Binding of antipsychotic drugs to cortical 5-HT2A receptors: a PET study of chlorpromazine, clozapine, and amisulpride in schizophrenic patients. Am. J .Psychiatry. 1998;155:505-8. [PubMed ID: 9545996].

  • 22.

    Fonseka T, Richter M, Müller D. Second Generation Antipsychotic-Induced Obsessive-Compulsive Symptoms in Schizophrenia: A Review of the Experimental Literature. Curr. Psychiatry Rep. 2014;16:1-17.

  • 23.

    Zamanian-Azodi M, Rezaei-Tavirani M, Rahmati-Rad S, Hasanzadeh H, Tavirani MR, Seyyedi SS. Protein-Protein Interaction Network could reveal the relationship between the breast and colon cancer. Gastroenterol. Hepatol. Bed Bench. 2015;8:215-24. [PubMed ID: 26328044].

  • 24.

    Zamanian Azodi M, Rezaei Tavirani M, Arefi Oskouie A, Hamdieh M, Derakhshan MK, Ahmadzadeh A, Zayeri F, Nejadi N, Rezaei Tavirani M, Mansouri V, Rostami-Nejad M, Vafaee R. Fluoxetine Regulates Ig Kappa Chain C Region Expression Levels in the Serum of Obsessive-Compulsive Disorder Patients: A proteomic Approach. Iran. J. Pharm. Res. 2017;16:1264-71. [PubMed ID: 29201116].

  • 25.

    Peyvandi H, Peyvandi AA, Safaei A, Azodi MZ, Rezaei-Tavirani M. Introducing Potential Key Proteins and Pathways in Human Laryngeal Cancer: A System Biology Approach. Iran. J. Pharm. Res. 2018;17:415-25. [PubMed ID: 29755572].

  • 26.

    Aranda B, Blankenburg H, Kerrien S, Brinkman FS, Ceol A, Chautard E, Dana JM, De Las Rivas J, Dumousseau M, Galeota E, Gaulton A, Goll J, Hancock RE, Isserlin R, Jimenez RC, Kerssemakers J, Khadake J, Lynn DJ, Michaut M, O′Kelly G, Ono K, Orchard S, Prieto C, Razick S, Rigina O, Salwinski L, Simonovic M, Velankar S, Winter A, Wu G, Bader GD, Cesareni G, Donaldson IM, Eisenberg D, Kleywegt GJ, Overington J, Ricard-Blum S, Tyers M, Albrecht M, Hermjakob H. PSICQUIC and PSISCORE: accessing and scoring molecular interactions. Nat. methods. 2011;8:528-9. [PubMed ID: 21716279].

  • 27.

    Assenov Y, Ramírez F, Schelhorn S-E, Lengauer T, Albrecht M. Computing topological parameters of biological networks. Bioinformatics. 2008;24:282-4. [PubMed ID: 18006545].

  • 28.

    Luo X, Huang L, Jia P, Li M, Su B, Zhao Z, Gan L. Protein-protein interaction and pathway analyses of top schizophrenia genes reveal schizophrenia susceptibility genes converge on common molecular networks and enrichment of nucleosome (chromatin) assembly genes in schizophrenia susceptibility loci. Schizophr. Bull. 2014;40:39-49. [PubMed ID: 23671194].

  • 29.

    Yu H, Kim PM, Sprecher E, Trifonov V, Gerstein M. The importance of bottlenecks in protein networks: correlation with gene essentiality and expression dynamics. PLoS Comput. Biol. 2007;3:e59. [PubMed ID: 17447836].

  • 30.

    Bindea G, Mlecnik B, Hackl H, Charoentong P, Tosolini M, Kirilovsky A, Fridman WH, Pagès F, rajanoski Z, Galon J. ClueGO: a Cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks. Bioinformatics. 2009;25:1091-3. [PubMed ID: 19237447].

  • 31.

    Bindea G, Galon J, Mlecnik B. CluePedia Cytoscape plugin: pathway insights using integrated experimental and in silico data. Bioinformatics. 2013;29:631-3.

  • 32.

    Team RC. R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2018. Available online at https://www.R-project.org/.

  • 33.

    Rezaei-Tavirani M, Zamanian-Azodi M, Rajabi S, Masoudi-Nejad A, Rostami-Nejad M, Rahmatirad S. Evaluation of protein clustering of Parkinson and Alzheimer′s diseases. Archives of Iranian medicine. Arch. Iran. Med. 2016;19:101-9.

  • 34.

    Zamanian-Azodi M, Rezaei Tavirani M, Rahmati-Rad S, Rezaei-Tavirani M. Ethanol and cancer induce similar changes on protein expression pattern of human fibroblast cells. Iranian Journal of Pharmaceutical Research. Iran. J. Pharm. Res. 2016;15:175-84.

  • 35.

    Zarnegarnia Y, Majd H, Tavirani M, Khaier N, Maboodi A. Application of fuzzy clustering in analysis of included proteins in esophagus, stomach and colon cancers based on similarity of Gene Ontology annotation. Koomesh. 2010;12:14-21.

  • 36.

    Rousseeuw PJ. Silhouettes: A graphical aid to the interpretation and validation of cluster analysis. J. Comput Appl. Math. 1987;20:53-65.

  • 37.

    Shannon Weickert C, Miranda-Angulo AL, Wong J, Perlman WR, Ward SE, Radhakrishna V, traub RE, Weinberger DR, Kleinman JE. Variants in the estrogen receptor alpha gene and its mRNA contribute to risk for schizophrenia. Hum. Mol. Genet. 2012;21:5238.

  • 38.

    Weickert CS, Miranda-Angulo AL, Wong J, Perlman WR, Ward SE, Radhakrishna V, traub RE, Weinberger DR, Kleinman JE. Variants in the estrogen receptor alpha gene and its mRNA contribute to risk for schizophrenia. Hum. Mol. Genet. 2008;17:2293-309. [PubMed ID: 18424448].

  • 39.

    Brock G, Pihur V, Datta S, Datta S. clValid, an R package for cluster validation. J. Stat. Softw. 2008;25:1-22.

  • 40.

    Wang Y-m, Zou L-q, Xie W-l, Yang Z-y, Zhu X-z, Cheung EFC, Sørensen TA, Møller A, Chan RCK. Altered Functional Connectivity of the Default Mode Network in Patients With Schizo-obsessive Comorbidity: A Comparison Between Schizophrenia and Obsessive-compulsive Disorder. Schizophr. Bull. 2019;45:199-210. [PubMed ID: 29365198].

  • 41.

    Hoffman KL, Cano-Ramírez H. Lost in translation? A critical look at the role that animal models of obsessive compulsive disorder play in current drug discovery strategies. Expert Opin. Drug Discov. 2018;13:211-20. [PubMed ID: 29275640].

  • 42.

    Kim Y, Giusti-Rodriguez P, Crowley JJ, Bryois J, Nonneman RJ, Ryan AK, Quackenbush CR, Iglesias-Ussel MD, Lee PH, Sun W, de Villena FP, Sullivan PF. Comparative genomic evidence for the involvement of schizophrenia risk genes in antipsychotic effects. Mol. Psychiatry. 2018;23:708-12. [PubMed ID: 28555076].

  • 43.

    Rezaei-Tavirani M, Rezaei-Tavirani S, Ahmadi N, Naderi N, Abdi S. Pancreatic adenocarcinoma protein-protein interaction network analysis. Gastroenterol Hepatol Bed Bench. 2017;10:85-92.

  • 44.

    Zamanian Azodi M, Rezaei Tavirani M, Arefi Oskouie A, Mansouri V, Hamdieh M, Nejati N, Hamid-pour M, Alireza Ahmadzadeh A, Rostami-Nejat M, Rezaei Tavirani M, Mortazavi Tabatabaei SA. Introducing Transthyretin as a Differentially Expressed Protein in Washing Subtype of Obsessive-Compulsive Disorder. Basic Clin. Neurosci. 2018;9:187-94. [PubMed ID: 30034649].

  • 45.

    Shepherd JE. Effects of estrogen on congnition mood, and degenerative brain diseases. J. Am. Pharm. Assoc. (Wash). 2001;41:221-8. [PubMed ID: 11297335].

  • 46.

    Yaryura‐Tobias JA, Neziroglu FA, Kaplan S. Self‐mutilation, anorexia and dysmenorrhea in obsessive compulsive disorder. Int. J. Eat. Disord. 1995;17:33-8. [PubMed ID: 7894450].

  • 47.

    Kindler J, Weickert CS, Skilleter AJ, Catts SV, Lenroot R, Weickert TW. Selective Estrogen Receptor Modulation Increases Hippocampal Activity during Probabilistic Association Learning in Schizophrenia. Neuropsychopharmacology. 2015;40:2388-97. [PubMed ID: 25829142].

  • 48.

    Palacios S. Selective estrogen receptor modulators: the future in menopausal treatment. Minerva. ginecol. 2011;63:275-86. [PubMed ID: 21654612].

  • 49.

    Kargozar R, Azizi H, Salari R. A review of effective herbal medicines in controlling menopausal symptoms. Electron Physician. 2017;9:5826-33. [PubMed ID: 29403626].

  • 50.

    Bedell S, Nachtigall M, Naftolin F. The pros and cons of plant estrogens for menopause. J. Steroid. Biochem. Mol. Biol. 2014;139:225-36. [PubMed ID: 23270754].

  • 51.

    Saghafi N, Ghazanfarpour M, Sadeghi R, Najarkolaei AH, Omid MG, Azad A, Bakhtiyari M, Hosseini Najarkolaei E. Effects of phytoestrogens in alleviating the menopausal symptoms: A systematic review and meta-analysis. Iran. J. Pharm. Res. 2017;16:99-111. [PubMed ID: 29844781].