A Novel Schizophrenia Diagnostic Model Based on Statistically Significant Changes in Gene Methylation in Specific Brain Regions

Zou, Donghua; Qiu, Yufen; Li, Rongjie; Meng, Youshi; Wu, Yuan

doi:https://doi.org/10.1155/2020/8047146

BioMed Research International

On this page

Abstract Introduction Materials and Methods Results Discussion Data Availability Conflicts of Interest Authors’ Contributions Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2020 | Article ID 8047146 | https://doi.org/10.1155/2020/8047146

A Novel Schizophrenia Diagnostic Model Based on Statistically Significant Changes in Gene Methylation in Specific Brain Regions

Donghua Zou,¹Yufen Qiu,²Rongjie Li,¹Youshi Meng,¹and Yuan Wu³

Academic Editor: Nikolai E. Lazarov

Received26 Jul 2019

Revised28 Oct 2019

Accepted15 Nov 2019

Published12 Feb 2020

Abstract

Objective. The present study identified methylation patterns of schizophrenia- (SCZ-) related genes in different brain regions and used them to construct a novel DNA methylation-based SCZ diagnostic model. Methods. Four DNA methylation datasets representing different brain regions were downloaded from the Gene Expression Omnibus. The common differentially methylated genes (CDMGs) in all datasets were identified to perform functional enrichment analysis. The differential methylation sites of 10 CDMGs involved in the largest numbers of neurological or psychiatric-related biological processes were used to construct a DNA methylation-based diagnostic model for SCZ in the respective datasets. Results. A total of 849 CDMGs were identified in the four datasets, but the methylation sites as well as degree of methylation differed across the brain regions. Functional enrichment analysis showed CDMGs were significantly involved in biological processes associated with neuronal axon development, intercellular adhesion, and cell morphology changes and, specifically, in PI3K-Akt, AMPK, and MAPK signaling pathways. Four DNA methylation-based classifiers for diagnosing SCZ were constructed in the four datasets, respectively. The sample recognition efficiency of the classifiers showed an area under the receiver operating characteristic curve of 1.00 in three datasets and >0.9 in one dataset. Conclusion. DNA methylation patterns in SCZ vary across different brain regions, which may be a useful epigenetic characteristic for diagnosing SCZ. Our novel model based on SCZ-gene methylation shows promising diagnostic power.

1. Introduction

Schizophrenia (SCZ) is a serious mental illness [1]. The World Health Organization estimates the global lifetime prevalence of SCZ at 3.8–8.4% [2]. SCZ is a severe psychosis induced by multiple factors and it manifests as a clinical syndrome with many symptoms. The course of the disease can include repeated relapses that aggravate disease and reduce quality of life. Some patients suffer from depression or mental disability [3]. Currently, clinical diagnosis of SCZ is based on the diagnostic scale of International Mental Disorders Classification [4]. The heterogeneous nature of SCZ pathogenesis has made it impossible so far to identify a single, reliable diagnostic biomarker or model.

Previous studies have shown that epigenetic changes may be related to the pathology of SCZ [5]. DNA methylation is the most stable epigenetic modification, and it can lead to changes in phenotype although the DNA sequence remains unchanged [6, 7]. DNA methylation may affect neuronal activity, transcriptional output, and synaptic function. Thus, methylation may be important in the pathology of SCZ [8]. Abnormal methylation of some genes, such as DRD2 [9], DLGAP2 [10], or COMT [11], may be associated with the occurrence and development of SCZ. In addition, a few DNA methylation-based classifiers for SCZ diagnosis have been reported [12, 13]. However, these studies were limited because they relied more on statistical associations without in-depth analysis of biological function.

In the present study, we identified common differentially methylated genes (CDMGs) shared across different brain regions in SCZ patients. We used the 10 CDMGs involved in the greatest number of neurological or psychiatric-related biological processes to construct a DNA methylation-based classifier for SCZ diagnosis. This method may be more biologically relevant than previous ones and may provide new insights to guide future research in SCZ.

2. Materials and Methods

2.1. DNA Methylation in SCZ

In the present study, four SCZ methylation datasets (GSE89702, GSE89703, GSE89705, and GSE89706) were downloaded from the Gene Expression Omnibus database (https://www.ncbi.nlm.nih.gov/geo/). DNA was isolated from postmortem brain samples. GSE89702 was derived from cerebellum from the Douglas Bell-Canada Brain Bank, and it included 16 SCZ samples and 17 normal controls. GSE89703 was derived from hippocampus, and it included 14 SCZ samples and 13 normal controls. GSE89705 was derived from striatum, and it included 16 SCZ samples and 17 normal controls. GSE89706 was derived from striatum from the London Brain Bank for Neurodegenerative Disorders, and it included 21 SCZ samples and 28 normal controls. The platform was GPL13534 and included each probe, the position on the chromosome, and the corresponding gene name of each probe. The normalized methylation data matrix was shown as beta values (ranging 0 to 1) of each probe with probe ID in row and patient ID in column. The workflow of the present study is shown in Figure 1.

2.2. Differential Methylation Analysis

Although GSE89705 and GSE89706 were both taken from the region of striatum, they were not combined because we did not know if there is any difference in the processing of samples between the two brain banks and that it is possible that combining datasets may result in some residual inflation according to a published study [14]. Thus, differential methylation analysis was performed separately in the four datasets using the limma package [15] in R software. The size of samples in each dataset was relatively small, and thus, the significance of differential methylation sites may be relatively low. If rigorously value filtering, genes with potential biological function may be filtered out. Therefore, in the present study, differences with a value <0.05 were considered significant.

2.3. Enrichment Analysis

To explore the biological functions of CDMGs that may be related to SCZ, Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) enrichment analyses were performed using the ClusterProfiler package [16] in R. GO terms and KEGG pathways with values <0.05 were considered significant.

2.4. Variable Selection and LASSO Classifier

Differentially methylated genes were involved in larger numbers of neurological or psychiatric-related biological processes and their methylation levels are thought to be more likely associated with SCZ than genes involved in fewer such processes. The differential methylation genes are different in different brain regions, and the differential methylation sites of the same gene in different brain regions may also be different. So we tried to build a brain-specific methylation-based classifier for SCZ. The corresponding different methylation sites of ten CDMGs involved in the most neurological or psychiatric-related biological processes were used to construct a diagnostic model of SCZ. These samples of the four datasets were randomly assigned to the training set (75%) and test set (25%), respectively. The four training sets were, respectively, used to select variables (different methylation sites) for establishing a DNA methylation-based classifier, and the test sets were used to validate the four classifiers. The glmnet package [17] in R used the least absolute shrinkage and selection operator (LASSO) [18] was used to select variables and construct a DNA methylation-based diagnostic classifier.

2.5. Evaluation of Methylation-Based SCZ Diagnostic Model

The diagnostic performance of the DNA methylation-based SCZ diagnostic model was evaluated by accuracy, sensitivity, specificity, positive predictive value, negative predictive value, and area under the receiver operating characteristic curve (AUC) as analyzed in the pROC package [19] in R software.

3. Results

3.1. SCZ-Related CDMGs in Different Brain Regions

We assessed DNA methylation characteristics from the four datasets taken from different regions of the brain. Compared with control groups from the respective datasets, we identified 7887 hypermethylated positions and 7484 hypomethylated positions in GSE89702 (Figure 2(a)), 5087 hypermethylated positions and 4227 hypomethylated positions in GSE89703 (Figure 2(b)), 5116 hypermethylated positions and 4275 hypomethylated positions in GSE89705 (Figure 2(c)), and 5569 hypermethylated positions and 4801 hypomethylated positions in GSE89706 (Figure 2(d)). A total of 849 genes were identified as CDMGs shared by all four datasets (Figure 2(e) and Table S1). However, these CDMGs were methylated at different sites and to different degrees in different brain regions.

(a)

(b)

(c)

(d)

(e)

3.2. CDMGs Involved in Multiple Neurological or Psychiatric-Related Biological Processes and Pathways

GO analysis of the SCZ-related CDMGs revealed these genes were involved in 244 biological processes, 43 cellular components, and 31 molecular functions (Figures 3(a)–3(d)). KEGG pathway analysis showed that the CDMGs were involved in 80 signaling pathways, most significantly in PI3K-Akt, AMPK, and m-activated protein kinases/extracellular regulated protein kinases (MAPK/ERK), as well as several pathways related to neurological or psychiatric-related biological processes (Figure 3(e)). Notably, the CDMGs were involved in multiple neuronal axon-related biological processes (Table S2); this suggests that the methylation level of CDMGs in neuronal axons may be associated with SCZ.

(a)

(b)

(c)

(d)

(e)

3.3. DNA Methylation-Based Diagnostic Model for SCZ

The following CDMGs were involved in the greatest numbers of neurological or psychiatric-related biological processes: SHANK3, WNT5A, NLGN1, GLI3, PTPRS, DISC1, SHH, BAIAP2, GLI2, and PAX6. The beta values of the differential methylation sites corresponding to these ten genes were used to construct a methylation template for diagnosing SCZ. Since our four datasets were taken from different regions of the brain, we found that these genes were differentially methylated based on their location—both in the brain region and on the DNA methylation position (Tables 1–4). Therefore, the beta values of the differential methylation sites corresponding to each gene were collected, and the LASSO method was used for variable selection and construction of the DNA methylation-based classifier. The results suggested that the counts of differential methylation sites selected by LASSO in different brain regions were various (Figure 4). The diagnostic ability of the same CDMG in different brain regions varied, which indicated that the epigenetic dysregulation of SCZ is complicated.

(a)

(b)

(c)

(d)

3.4. Diagnostic Efficiency of the Composite Model in Training and Validation

In order to evaluate the diagnostic efficiency of the DNA methylation-based classifier, receiver operating characteristic curves were analyzed (Figure 5). In GES89702, GES89703, and GES89705, the accuracy, sensitivity, specificity, positive predictive value, negative predictive value, and AUC in the training and test sets were 1 (Tables 5–7). In GES89706, the accuracy was 0.950 in the training set and 0.920 in the test set, while the AUC was 0.994 in the training set and 0.943 in the test set (Table 8). The results suggest the DNA methylation-based classifier is a potential biomarker for diagnosing SCZ.

(a)

(b)

(c)

(d)

4. Discussion

In recent years, the morbidity and mortality rates of SCZ have increased, and the many health- and social-related problems for patients with SCZ are cause for much concern. However, the effective molecular diagnostic methods are unmet. In particular, diagnostic models that take into account both the molecular statistical and biological significance have not received much attention.

In the present study, a total of 849 CDMGs were identified in different brain regions. Functional enrichment analysis indicated the CDMGs were involved in various neurological or psychiatric-related biological processes and pathways, specifically signaling pathways PI3K-Akt, AMPK, and MAPK. The methylation levels of CDMGs may affect these biological processes and pathways. Our study identified biological processes with confirmed roles in mental diseases, including SCZ [20, 21]. MAPK/ERK and PI3K/Akt signaling pathways can regulate protuberant growth and protein synthesis related to neural plasticity, and it can assist in the normal development of nerve cells, which may protect against SCZ [22]. The initiation of nerve axon regeneration is regulated by the MAPK pathway and this initiates a neuronal response [23]. Further studies should explore directly whether the CDMGs in the present study are associated with SCZ.

The beta values of the differential methylation sites corresponding to 10 CDMGs (SHANK3, WNT5A, NLGN1, GLI3, PTPRS, DISC1, SHH, BAIAP2, GLI2, and PAX6) which were most often involved in neurological or psychiatric-related biological processes were used to construct a DNA methylation-based brain region-specific for SCZ. The methylation sites within genes and degree of methylation varied in different brain regions, suggesting that the methylation patterns of SCZ-related genes are extremely complex. The DISC1 protein regulates the development, maturation, and migration of brain neurons and synaptic signal transmission [24–27], and disruption can lead to various mental diseases, including SCZ [28, 29]. The process of nerve development and synaptic transmission regulated by DISC1 can be affected by its degree of methylation [30]. SHANK3 knock out may affect neuronal development and induce SCZ [31]. SHANK3 and NLGN1 are also related to the progression of SCZ [32, 33]. However, few reports exist on the association of genes WNT5A, GLI3, PTPRS, SHH, BAIAP2, GLI2, or PAX6 with SCZ. Our results suggest that the methylation level of these genes may be related to the disease. Indeed, our DNA methylation-based classifier showed strong diagnostic potential based on AUC analysis. It is worth noting that due to the characteristics of LASSO method, the more the inclusion variables, the better the effect of the classifier. However, from the perspective of economics, the more the variables (differential methylation sites) are included, the higher the cost is. So taking into account the effectiveness of the model and the cost of economics, we started to perform feature section and classifier construction from 10 genes instead of fewer or more. From the results of the present results, when we included 10 genes, LASSO method identified differential methylation sites of 7-8 genes for us and obtained the best AUC value (close to 1). So we may foresee that the inclusion of fewer genes may greatly reduce classification efficiency, while the inclusion of more genes may be not necessary because it would increase costs but not increase the effect.

Some critical limitations exist in the present study. Due to the small sample size, our DNA methylation-based SCZ diagnostic model needs to be further validated and improved in larger, independent datasets. The potential roles of CDMGs in SCZ need to be explored experimentally.

Despite these limitations, our findings suggest that gene methylation patterns are significantly associated with SCZ and may be a promising diagnostic method. Methylation levels and sites in CDMGs varied widely across different brain regions, and future studies should explore the potential relevance of this variation for SCZ onset and progression.

Data Availability

The data used to support the findings of this study are included within Supplementary Table S1.

Conflicts of Interest

The authors report no conflicts of interest in this work.

Authors’ Contributions

Donghua Zou and Yufen Qiu contributed equally to this work.

Acknowledgments

This study was supported by the Guangxi Natural Science Foundation (2016GXNSFCA380012), the Project of Nanning Scientific Research and Technology Development Plan (20193093), the High-Level Medical Expert Training Program of Guangxi “139” Plan Funding (G201903049) and sponsored by Nanning Excellent Young Scientist Program (RC20190103).

Supplementary Materials

Table S1: common differential methylated genes in four datasets. Table S2: common differential methylated genes involved in neurological or psychiatric-related biological processes. (Supplementary Materials.)

References

G. Pergola, P. Di Carlo, A. E. Jaffe et al., “Prefrontal coexpression of schizophrenia risk genes is associated with treatment response in patients,” Biological Psychiatry, vol. 86, no. 1, pp. 45–55, 2019.
View at: Publisher Site | Google Scholar
M. Avram, F. Brandl, J. Cabello et al., “Reduced striatal dopamine synthesis capacity in patients with schizophrenia during remission of positive symptoms,” Brain, vol. 142, no. 6, pp. 1813–1826, 2019.
View at: Publisher Site | Google Scholar
Z. Hu, X. Xiao, Z. Zhang, and M. Li, “Genetic insights and neurobiological implications from NRXN1 in neuropsychiatric disorders,” Molecular Psychiatry, vol. 24, no. 10, pp. 1400–1414, 2019.
View at: Publisher Site | Google Scholar
S. Guloksuz, L. K. Pries, P. Delespaul et al., “Examining the independent and joint effects of molecular genetic liability and environmental exposures in schizophrenia: results from the EUGEI study,” World Psychiatry, vol. 18, no. 2, pp. 173–182, 2019.
View at: Publisher Site | Google Scholar
C. Zhuo, D. Wang, C. Zhou et al., “Double-edged sword of tumour suppressor genes in schizophrenia,” Frontiers in Molecular Neuroscience, vol. 12, p. 1, 2019.
View at: Publisher Site | Google Scholar
M. V. Micioni Di Bonaventura, M. Pucci, M. E. Giusepponi et al., “Regulation of adenosine A2A receptor gene expression in a model of binge eating in the amygdaloid complex of female rats,” Journal of Psychopharmacology, vol. 33, no. 12, pp. 1550–1561, 2019.
View at: Publisher Site | Google Scholar
M. Ehrlich, “DNA hypermethylation in disease: mechanisms and clinical relevance,” Epigenetics, vol. 14, no. 12, pp. 1141–1163, 2019.
View at: Publisher Site | Google Scholar
M. Herre and E. Korb, “The chromatin landscape of neuronal plasticity,” Current Opinion in Neurobiology, vol. 59, pp. 79–86, 2019.
View at: Publisher Site | Google Scholar
Y. Funahashi, Y. Yoshino, K. Yamazaki et al., “Analysis of methylation and -141C Ins/Del polymorphisms of the dopamine receptor D2 gene in patients with schizophrenia,” Psychiatry Research, vol. 278, pp. 135–140, 2019.
View at: Publisher Site | Google Scholar
Y. Li, K. Wang, P. Zhang et al., “Quantitative DNA methylation analysis of DLGAP2 gene using pyrosequencing in schizophrenia with tardive dyskinesia: a linear mixed model approach,” Scientific Reports, vol. 8, no. 1, p. 17466, 2018.
View at: Publisher Site | Google Scholar
S. Gao, J. Cheng, G. Li et al., “Catechol-O-methyltransferase gene promoter methylation as a peripheral biomarker in male schizophrenia,” European Psychiatry, vol. 44, pp. 39–46, 2017.
View at: Publisher Site | Google Scholar
A. Yoshimi, S. Yamada, S. Kunimoto et al., “Proteomic analysis of lymphoblastoid cell lines from schizophrenic patients,” Translational Psychiatry, vol. 9, no. 1, p. 126, 2019.
View at: Publisher Site | Google Scholar
V. L. Nimgaonkar, F. Dickerson, J. G. Pouget et al., “Joint evaluation of serum C-reactive protein levels and polygenic risk scores as risk factors for schizophrenia,” Psychiatry Research, vol. 261, pp. 148–153, 2018.
View at: Publisher Site | Google Scholar
J. Viana, E. Hannon, E. Dempster et al., “Schizophrenia-associated methylomic variation: molecular signatures of disease and polygenic risk burden across multiple brain regions,” Human Molecular Genetics, vol. 26, no. 1, pp. 210–225, 2017.
View at: Publisher Site | Google Scholar
M. E. Ritchie, B. Phipson, D. Wu et al., “Limma powers differential expression analyses for RNA-sequencing and microarray studies,” Nucleic Acids Research, vol. 43, no. 7, p. e47, 2015.
View at: Publisher Site | Google Scholar
G. Yu, L.-G. Wang, Y. Han, and Q.-Y. He, “Clusterprofiler: an R package for comparing biological themes among gene clusters,” OMICS: A Journal of Integrative Biology, vol. 16, no. 5, pp. 284–287, 2012.
View at: Publisher Site | Google Scholar
T. T. Wu, Y. F. Chen, T. Hastie, E. Sobel, and K. Lange, “Genome-wide association analysis by lasso penalized logistic regression,” Bioinformatics, vol. 25, no. 6, pp. 714–721, 2009.
View at: Publisher Site | Google Scholar
N. Simon, J. Friedman, T. Hastie, and R. Tibshirani, “Regularization paths for cox’s proportional hazards model via coordinate descent,” Journal of Statistical Software, vol. 39, no. 5, pp. 1–13, 2011.
View at: Publisher Site | Google Scholar
X. Robin, N. Turck, A. Hainard et al., “pROC: an open-source package for R and S+ to analyze and compare ROC curves,” BMC Bioinformatics, vol. 12, no. 1, p. 77, 2011.
View at: Publisher Site | Google Scholar
A. Mishra, M. H. Traut, L. Becker, T. Klopstock, V. Stein, and R. Klein, “Genetic evidence for the adhesion protein IgSF9/Dasm1 to regulate inhibitory synapse development independent of its intracellular domain,” Journal of Neuroscience, vol. 34, no. 12, pp. 4187–4199, 2014.
View at: Publisher Site | Google Scholar
A. R. Filous, A. Tran, C. J. Howell et al., “Entrapment via synaptic-like connections between NG2 proteoglycan+ cells and dystrophic axons in the lesion plays a role in regeneration failure after spinal cord injury,” The Journal of Neuroscience, vol. 34, no. 49, pp. 16369–16384, 2014.
View at: Publisher Site | Google Scholar
R. S. Ferreira, N. A. G. Dos Santos, C. P. Bernardes et al., “Caffeic acid phenethyl ester (CAPE) protects PC₁₂ cells against cisplatin-induced neurotoxicity by activating the AMPK/SIRT₁, MAPK/Erk, and PI3k/Akt signaling pathways,” Neurotoxicity Research, vol. 36, no. 1, pp. 175–192, 2019.
View at: Publisher Site | Google Scholar
N. Hisamoto, T. Shimizu, K. Asai et al., “C. elegans tensin promotes axon regeneration by linking the met-like SVH-2 and integrin signaling pathways,” The Journal of Neuroscience, vol. 39, no. 29, pp. 5662–5672, 2019.
View at: Publisher Site | Google Scholar
K. K. Singh, X. Ge, Y. Mao et al., “Dixdc1 is a critical regulator of DISC1 and embryonic cortical development,” Neuron, vol. 67, no. 1, pp. 33–48, 2010.
View at: Publisher Site | Google Scholar
K. K. Singh, G. De Rienzo, L. Drane et al., “Common DISC1 polymorphisms disrupt Wnt/GSK3β signaling and brain development,” Neuron, vol. 72, no. 4, pp. 545–558, 2011.
View at: Publisher Site | Google Scholar
G. De Rienzo, J. A. Bishop, Y. Mao et al., “Disc1 regulates both β-catenin-mediated and noncanonical Wnt signaling during vertebrate embryogenesis,” The FASEB Journal, vol. 25, no. 12, pp. 4184–4197, 2011.
View at: Publisher Site | Google Scholar
D. Deng, C. Jian, L. Lei et al., “A prenatal interruption of DISC1 function in the brain exhibits a lasting impact on adult behaviors, brain metabolism, and interneuron development,” Oncotarget, vol. 8, no. 49, pp. 84798–84817, 2017.
View at: Publisher Site | Google Scholar
S. D. Greenhill, K. Juczewski, A. M. de Haan, G. Seaton, K. Fox, and N. R. Hardingham, “Adult cortical plasticity depends on an early postnatal critical period,” Science, vol. 349, no. 6246, pp. 424–427, 2015.
View at: Publisher Site | Google Scholar
Y. Mao, X. Ge, C. L. Frank et al., “Disrupted in schizophrenia 1 regulates neuronal progenitor proliferation via modulation of GSK3β/β-catenin signaling,” Cell, vol. 136, no. 6, pp. 1017–1031, 2009.
View at: Publisher Site | Google Scholar
E. Bentea, E. A. K. Depasquale, S. M. O’Donovan et al., “Kinase network dysregulation in a human induced pluripotent stem cell model of DISC1 schizophrenia,” Molecular Omics, vol. 15, no. 3, pp. 173–188, 2019.
View at: Publisher Site | Google Scholar
G. Huang, S. Chen, X. Chen et al., “Uncovering the functional link between SHANK3 deletions and deficiency in neurodevelopment using iPSC-derived human neurons,” Frontiers in Neuroanatomy, vol. 13, p. 23, 2019.
View at: Publisher Site | Google Scholar
X. Chen, F. Long, B. Cai, X. Chen, and G. Chen, “A novel relationship for schizophrenia, bipolar and major depressive disorder part 3: evidence from chromosome 3 high density association screen,” Journal of Comparative Neurology, vol. 526, no. 1, pp. 59–79, 2018.
View at: Publisher Site | Google Scholar
B. McKinney, Y. Ding, D. A. Lewis, and R. A. Sweet, “DNA methylation as a putative mechanism for reduced dendritic spine density in the superior temporal gyrus of subjects with schizophrenia,” Translational Psychiatry, vol. 7, no. 2, p. e1032, 2017.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2020 Donghua Zou et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

761

Downloads

1290

Citations