PPAR Research
Volume 2020 (2020), Article ID 3480412, 6 pages
https://doi.org/10.1155/2020/3480412

Research Article

PPARD May Play a Protective Role against the Development of Schizophrenia

Xinrong Li,1,2 Sha Liu,1,2 Karan Kapoor,³ and Yong Xu1,2

¹Shanxi Key Laboratory of Artificial Intelligence Assisted Diagnosis and Treatment for Mental Disorder, First Hospital of Shanxi Medical University, Taiyuan, China
²Department of Psychiatry, First Hospital/First Clinical Medical College of Shanxi Medical University, Taiyuan, China
³NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA

Correspondence should be addressed to Yong Xu; xuyongsmu@vip.163.com

Received 28 February 2020; Accepted 29 May 2020; Published 7 August 2020

Guest Editor: Anastasia Nesterova

Copyright © 2020 Xinrong Li 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.

Abstract.

PPARD has been suggested to contribute to the etiology of schizophrenia (SCZ) with the underlying mechanisms largely unknown. Here, we first collected and analyzed the PPARD expression profile from three groups: (1) 18 healthy control (HC) subjects, (2) 14 clinical high-risk (CHR) patients, and (3) 19 early onset of SCZ (EOS) patients. After that, we conducted a systematical pathway analysis to explore the potential mechanisms involved in PPARD exerting influence on the pathological development of SCZ. Compared to the HC group, the expression of PPARD was slightly decreased in the EOS group (; ) and increased in the CHR group (; ). However, there was a significant difference between the EOS group and the CHR group (; ), reflecting the amount of variation in PPARD expression before and after the onset of SCZ. Pathway analysis suggested that overexpression of PPARD may regulate ten proteins or molecules to inhibit the pathological development of SCZ, including the deactivation of eight SCZ promoters and stimulation of two SCZ inhibitors. Our results support the association between PPARD and SCZ. The pathways identified may help in the understanding of the potential mechanisms by which PPARD contributes to the etiology of SCZ.

1. Introduction

Schizophrenia (SCZ) is a common and often disabling mental illness characterized not only by a varied group of clinical symptoms [1], but wide-ranging deficits in neurocognitive and neurophysiological functions [2, 3]. The prodromal period is thought to have a high risk of clinical symptoms and precedes illness onset by 1 to 6 years [4, 5]. Subjects with these characteristics are called clinical high-risk (CHR) patients, with about one-third developing SCZ and two-third recovering to normal [4, 5] (PMID: 8782291; PMID: 1571314). Early-onset SCZ (EOS), defined as SCZ with onset before the 21st birthday, shows worse psychosocial disability and poor prognosis [6]. As a neurodevelopment disease [7], SCZ in childhood tends to have a higher possibility of abnormal neural development [8].

PPARD is a nuclear hormone receptor that governs a variety of biological processes [9]. This gene has been suggested to play roles in the development of several chronic diseases, including diabetes, obesity, atherosclerosis, and cancer [10]. Several studies have suggested that PPARG may contribute to the etiology of SCZ [11–13]. For instance, Sun et al.’s study showed that the PPARD polymorphism rs2076169 had an allelic association with SCZ (, ) in a trio study using a transmission disequilibrium test [11]. Maekawa et al. identified a significantly downregulated expression of PPARD () in individuals with SCZ compared with the control subjects [12]. Dzana et al. also discovered linkages between the genetic variants of multiple genes including PPARD and the increased waist circumference in SCZ patients () [13]. However, the underlying mechanism regarding the PPARD-SCZ association is largely unknown.

To explore the relationship between PPARD and SCZ, we studied the expression changes of PPARD in both CHR and EOS groups and compared that to the healthy control (HC) group. CHR is a special state of SCZ before its onset. We hypothesized that the changes in PPARD expression in the CHR group could lead to the regulation of genes not observed in the SCZ group. After that, we conducted a systematic bioinformatics analysis and identified multiple pathways through which PPARD could exert influence on SCZ. Our study provides novel evidence for the association between PPARD and SCZ and adds new insights into the understanding of the roles of PPARD in the etiology of SCZ.

2. Materials and Methods

2.1. Subject Recruitment for Expression Profile Collection

All participants were unrelated Han Chinese recruited from the north of China and under the age of 18 years. EOS patients were diagnosed by 2 associate doctors according to the Diagnostic and Statistical Manual of Mental Disorders: Fourth Edition (DSM-IV) and the Chinese version of the Modified Structured Clinical Interview for DSM-IV, patient version (SCID-I/P). The total score of PANSS was ≥60, and IQ score was ≥70. Exclusions included patients with organic disease of heart, liver, and kidney; all kinds of immune diseases, brain injury, or brain congenital malformation; a tumor of brain and epilepsy; mental retardation, along with anyone taking antipsychotic drugs, antimanic drugs, antidepressants, or mood stabilizers. In addition, we excluded serious excitement or impulsion patients. CHR patients were assessed with the Structured Interview of Prodromal Syndromes (SIPS) [14]. The exclusion criteria of this group were the same as the EOS group. The healthy controls (HC) were selected by matching the age and sex with never having taken any drugs in the latest one month. Exclusion criteria included (1) meeting the standards of inclusion or exclusion term of patients, (2) having the family history of any spirit or nervous system disease, (3) having head injury or newborn related disease, (4) having hyperpyretic convulsion before, and (5) being an adopted child or living in a single-parent family.

All teenage participants’ informed consent was signed by their parents. The study was conducted under the protocols approved by the First Hospital of Shanxi Medical University (Ethical Code: 2019-Y01).

2.2. RNA Extraction and Quantity Control

Total RNA was extracted from all of the samples which had been snap-frozen using TRIzol reagent (Invitrogen, Carlsbad, CA, U.S.) according to the manufacturer’s previous protocol [15]. Total RNA from each sample was quantified by the NanoDrop ND-1000, and RNA integrity was assessed by standard denaturing agarose gel electrophoresis.

2.3. RNA Labeling and Microarray Hybridization

The RNA labeling and microarray hybridization followed the routine process described as follows [16]. The Arraystar Human LncRNA Array v2.0 is designed for researchers who are interested in profiling both lncRNAs and protein-coding RNAs in the human genome. Sample labeling and array hybridization were performed according to the Agilent One-Color Microarray-Based Gene Expression Analysis protocol (Agilent Technology) with minor modifications. Briefly, mRNA was purified from total RNA after removal of rRNA (mRNA-ONLY™ Eukaryotic mRNA Isolation Kit, Epicentre). Then, each sample was amplified and transcribed into fluorescent cRNA along the entire length of the transcripts without 3 bias utilizing a random priming method. The labeled cRNAs were purified by the RNeasy Mini Kit (Qiagen). The concentration and specific activity of the labeled cRNAs (pmol Cy3/μg cRNA) were measured by NanoDrop ND-1000. 1 μg of each labeled cRNA was fragmented by adding 11 μl 10x Blocking Agent and 2.2 μl of 25x Fragmentation Buffer and then heating the mixture at 60°C for 30 min. Finally, 55 μl 2x GE Hybridization buffer was added to dilute the labeled cRNA. 100 μl of hybridization solution was dispensed into the gasket slide and assembled on the RNA expression microarray slide. The slides were incubated for 17 hours at 65°C in an Agilent Hybridization Oven. The hybridized arrays were washed, fixed, and scanned using the Agilent DNA Microarray Scanner (part number G2505C). The microarray work was performed by KangChen Bio-tech (Shanghai).

2.4. Bioinformatics Analysis

To gain a better understanding of the gene expression resulting from different groups of subjects and explore possible roles of PPARD in the etiology of SCZ, we conducted a literature-based pathway analysis to identify the possible molecular pathways connecting PPARD and SCZ. Specifically, by using the tool Pathway studio (version 12.3; http://www.pathwaystudio.com), we identified genes and small molecules that are downstream targets of PPARD and upstream regulators of SCZ with polarity. Then, we constructed the PPARD-driven functional pathways with polarity and direction.

3. Results

3.1. Demographics

The three groups of subjects were comparable in age and gender. In all, we recruited 19 EOS patients (8 males and 11 females, aged years), 14 CHR patients (9 males and 5 females, aged years), and 18 HC (9 males and 9 females, aged years). The demographic information for all participants is provided in Table 1.

Table 1: Demographics and clinical characteristics for all participants.

EOS patients
CHR patients
HC
F/
value

19
14
18

Age (years)

2.012
0.145

Gender (M/F)
8/11
9/5
9/9
1.598
0.450

PANSS total scores

Note: PANSS; Positive and Negative Syndrome Scale.

3.2. Expression Variation of PPARD in Different Groups

Compared with the HC group, PPARD presented increased expression levels in the CHR group () and decreased expression levels in the EOS group (), as shown in Figure 1(a). The changes were milder in terms of the statistical value ( value = 0.23 and 0.20 for EOS vs HC and CHR vs HC, respectively). However, the difference between EOS and CHR group showed statistical significance ( value = 0.015; ). The downregulation of PPARD in the EOS group was consistent with previous study results (PubMed 28872641); however, the increased PPARD expression levels in the CHR groups were not reported before. Usually, there are about two-third of the CHR subjects who do not develop SCZ (PMID: 8782291; PMID: 1571314). Thus, our results indicated that increased PPARD expression might play a protective role against the development of SCZ, with the underlying mechanism explored using subsequent pathway analysis.

(a)

(b)

Figure 1: PPARD expression comparison among three different groups: healthy control (HC) group, clinical high-risk (CHR) group, and early-onset schizophrenia (EOS) group. (a) The -10  log10 transferred values in different comparisons; (b) the log fold change of PPARD expression for different comparisons.

3.3. Genetic Pathways Driven by PPARD

To understand the possible roles that PPARD could exert on the pathologic development of SCZ, we constructed a literature-based genetic pathway connecting PPARD and SCZ, as shown in Figure 2. Our results showed that PPARD deactivates four promoters of SCZ, including CNR1, AGTR1, ACAN, and IL1B. Moreover, one SCZ inhibitor could also get activated by PPARD. These results may partially explain the mechanism regarding the roles of PPARD in the etiology of SCZ. Each relation within Figure 2 was supported by one or more scientific references (Supplementary Material: Ref 4 Figure 2). The corresponding sentences where a relationship has been identified were reviewed to confirm the confidence of the identified relation.

Figure 2: PPARD driven genetic pathways inhibiting schizophrenia. The genes highlighted in red are schizophrenia promoters; green denoting schizophrenia inhibitors.

3.4. Molecule Pathways Driven by PPARD

Besides the genetic pathway, we also identified multiple molecules influencing the pathologic development of SCZ and regulated by PPARG, as shown in Figure 3. Specifically, PPARG promotes the secretion of glutathione, which is an inhibitor of SCZ. Moreover, PPARD suppresses four molecules that could promote the development of SCZ, including superoxide, ROS, glutamate, and fatty acid. These results may further explain the underlying mechanism in which increased expression of PPARD could protect against the development of SCZ. Each relation within Figure 3 was supported by one or more scientific references (Supplementary Material: Ref 4 Figure 3). The corresponding sentences where a relationship has been identified were reviewed to confirm the confidence of the identified relation.

Figure 3: PPARD driven genetic pathways inhibiting schizophrenia. The molecules highlighted in red are schizophrenia promoters; green denoting schizophrenia inhibitors.

4. Discussion

Previous studies have suggested that PPARD may contribute to the etiology of SCZ. However, the underlying mechanism is largely unknown [11, 12]. To better understand the role of PPARD, we tested the expression variation of PPARD in both CHR and EOS groups and compared it to the HC group. We hypothesize that PPARD expression changes in the CHR group may reflect the activities of genes during the development of SCZ but before its onset. Our results showed that the EOS group presented decreased expression levels compared to healthy controls, which was consistent with a previous study [12]. However, we observed increased PPARD expression in the CHR group compared to both HC and EOS groups. The opposite change between CHR and EOS groups may reflect the PPARD variation before and after the onset of SCZ.

Pathway analysis showed that increased expression of PPARD might drive proteins and small molecules to protect against the pathologic development of SCZ, as shown in Figures 2 and 3. These pathways suggested that increased expression of PPARD could inhibit SCZ promoters and activate SCZ inhibitors, consequently influencing the pathophysiology of SCZ. For instance, PPARD has been shown to inhibit CB1 receptor expression (CNR1), which contributes to the pathophysiology of SCZ [17]. PPARD also reduces the upregulation of angiotensin II type 1 receptor (AGTR1) [18], the antagonists of which have been reported to improve clinical symptoms in SCZ patients [19]. In addition, PPARD activation promotes the degradation of aggrecan (ACAN) and attenuates gene expression of IL1B [20, 21]. Both ACAN and IL1B were suggested to contribute to the increased risk of SCZ [22, 23]. Moreover, SCZ patients have been shown to present decreased expression of tumor suppressor gene TXNIP, a thioredoxin-binding protein that is a member of the alpha arrestin protein family [24]. Ratneswaran et al. showed that PPARD could upregulate the expression of TXNIP [25]. These results support the association between PPARD and SCZ and also suggest the possible mechanisms of the protective roles of PPARD in the pathologic development of SCZ.

Our study also shows that PPARD may inhibit the generation and release of four molecules that play important roles in the pathophysiology of SCZ, including superoxide, oxygen free radicals (ROS), glutamate, and fatty acid (Figure 3). Moreover, it has been shown that agonists of PPARD could promote glutathione synthesis [26], with the deficit of glutathione impairing neurotransmission and cerebral connectivity that lead to clinical symptoms of SCZ [27]. These molecule pathways add more support for the mechanisms involved in the PPARD-SCZ association.

To our knowledge, no previous study has explored the PPARD expression variations between CHR patients and the EOS group, which provide new vision of the activity of PPARD in the case of SCZ. This study also has some limitations. Firstly, larger datasets could be collected to validate the expression levels of PPARD in CHR and EOS patients. Secondly, a follow-up study should be conducted to identify the expression of PPARD throughout the clinical high-risk stage to the recovery or disease stage of SCZ. Due to the limitation of supporting sources, this work has been left for future studies.

5. Conclusion

Our results confirmed decreased expression of PPARD in the case of SCZ and revealed increased expression in the clinical high-risk group. Pathway analysis suggested that the overexpression of PPARD in the high-risk group may contribute to the protection of a subject from developing SCZ.

Data Availability

The data of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

All the authors declare no conflict of interest.

Acknowledgments

This work was supported by the National Key Research and Development Program of China (2016YFC1307004), the National Natural Science Foundation of China (81971601 and 81701326), the Support Program of the Youth Sanjin Scholars, the 136 Medical Rejuvenation Project of Shanxi Province, and the Multidisciplinary Team for Cognitive Impairment of Shanxi Science and Technology Innovation Training Team (201705D131027) and Youth Science, Technology Research Fund of Shanxi (201801D221418), and Scientific Research Project of Shanxi Health Commission (201601034).

Supplementary Materials

| Reference for Figure 2 and reference for Figure 3. (Supplementary Materials)

References

C. A. Tamminga and H. H. Holcomb, “Phenotype of schizophrenia: a review and formulation,” Molecular Psychiatry, vol. 10, no. 1, pp. 27–39, 2005.

J. van Os, G. Kenis, and B. P. F. Rutten, “The environment and schizophrenia,” Nature, vol. 468, no. 7321, pp. 203–212, 2010.

B. I. Turetsky, M. E. Calkins, G. A. Light, A. Olincy, A. D. Radant, and N. R. Swerdlow, “Neurophysiological endophenotypes of schizophrenia: the viability of selected candidate measures,” Schizophrenia Bulletin, vol. 33, no. 1, pp. 69–94, 2007.

A. R. Yung and P. D. McGorry, “The prodromal phase of first-episode psychosis: past and current conceptualizations,” Schizophrenia Bulletin, vol. 22, no. 2, pp. 353–370, 1996.

H. Häfner, A. Riecher-Rössler, M. Hambrecht et al., “IRAOS: an instrument for the assessment of onset and early course of schizophrenia,” Schizophrenia Research, vol. 6, no. 3, pp. 209–223, 1992.

L. Clemmensen, D. L. Vernal, and H. C. Steinhausen, “A systematic review of the long-term outcome of early onset schizophrenia,” BMC Psychiatry, vol. 12, no. 1, 2012.

J. L. Rapoport, J. N. Giedd, and N. Gogtay, “Neurodevelopmental model of schizophrenia: update 2012,” Molecular Psychiatry, vol. 17, no. 12, pp. 1228–1238, 2012.

N. Gogtay, “Cortical brain development in schizophrenia: insights from neuroimaging studies in childhood-onset schizophrenia,” Schizophrenia Bulletin, vol. 34, no. 1, pp. 30–36, 2007.

J. Berger and D. E. Moller, “The mechanisms of action of PPARs,” Annual Review of Medicine, vol. 53, no. 1, pp. 409–435, 2002.

J. N. Feige, L. Gelman, L. Michalik, B. Desvergne, and W. Wahli, “From molecular action to physiological outputs: peroxisome proliferator-activated receptors are nuclear receptors at the crossroads of key cellular functions,” Progress in Lipid Research, vol. 45, no. 2, pp. 120–159, 2006.

S. L. Sun, G. Z. Ju, Y. Liu et al., “281 – association of polymorphism of PPARD gene with schizophrenia in Chinese Han population,” Schizophrenia Research, vol. 98, p. 149, 2008.

M. Maekawa, A. Watanabe, Y. Iwayama et al., “Polyunsaturated fatty acid deficiency during neurodevelopment in mice models the prodromal state of schizophrenia through epigenetic changes in nuclear receptor genes,” Translational Psychiatry, vol. 7, no. 9, p. e1229, 2017.

D. S. Hukic, U. Ösby, E. Olsson et al., “Genetic variants of increased waist circumference in psychosis,” Psychiatric Genetics, vol. 27, no. 6, pp. 210–218, 2017.

T. J. Miller, T. H. McGlashan, S. W. Woods et al., “Symptom assessment in schizophrenic prodromal states,” The Psychiatric Quarterly, vol. 70, no. 4, pp. 273–287, 1999.

Y. Xu, F. Li, B. Zhang et al., “MicroRNAs and target site screening reveals a pre-microRNA-30e variant associated with schizophrenia,” Schizophrenia Research, vol. 119, no. 1-3, pp. 219–227, 2010.

Y. Lou, H. Jiang, Z. Cui, X. Wang, L. Wang, and Y. Han, “Gene microarray analysis of lncRNA and mRNA expression profiles in patients with high-grade ovarian serous cancer,” International Journal of Molecular Medicine, vol. 42, no. 1, pp. 91–104, 2018.

A. Seillier, A. A. Martinez, and A. Giuffrida, “Phencyclidine-induced social withdrawal results from deficient stimulation of cannabinoid CB₁ receptors: implications for schizophrenia,” Neuropsychopharmacology, vol. 38, no. 9, pp. 1816–1824, 2013.

M. Romero, R. Jimenez, M. Toral et al., “Vascular and central activation of peroxisome proliferator-activated receptor-β attenuates angiotensin II-induced hypertension: role of RGS-5,” The Journal of Pharmacology and Experimental Therapeutics, vol. 358, no. 1, pp. 151–163, 2016.

S. Mohite, S. M. de Campos-Carli, N. P. Rocha et al., “Lower circulating levels of angiotensin-converting enzyme (ACE) in patients with schizophrenia,” Schizophrenia Research, vol. 202, pp. 50–54, 2018.

A. Ratneswaran, E. A. LeBlanc, E. Walser et al., “Peroxisome Proliferator–Activated Receptor δ Promotes the Progression of Posttraumatic Osteoarthritis in a Mouse Model,” Arthritis & Rhematology, vol. 67, no. 2, pp. 454–464, 2015.

C. Wang, G. Zhou, and Z. Zeng, “Effects of peroxisome proliferator-activated receptor-β/δ on sepsis induced acute lung injury,” Chinese Medical Journal, vol. 127, no. 11, pp. 2129–2137, 2014.

H. Pantazopoulos, M. Markota, F. Jaquet et al., “Aggrecan and chondroitin-6-sulfate abnormalities in schizophrenia and bipolar disorder: a postmortem study on the amygdala,” Translational Psychiatry, vol. 5, no. 1, p. e496, 2015.

M. Yoshida, K. Shiroiwa, K. Mouri et al., “Haplotypes in the expression quantitative trait locus of interleukin-1β gene are associated with schizophrenia,” Schizophrenia Research, vol. 140, no. 1-3, pp. 185–191, 2012.

Y. Su, W. Ding, M. Xing, D. Qi, Z. Li, and D. Cui, “The interaction of TXNIP and AFq1 genes increases the susceptibility of schizophrenia,” Molecular Neurobiology, vol. 54, no. 6, pp. 4806–4812, 2017.

A. Ratneswaran, M. M.-G. Sun, H. Dupuis, C. Sawyez, N. Borradaile, and F. Beier, “Nuclear receptors regulate lipid metabolism and oxidative stress markers in chondrocytes,” Journal of Molecular Medicine (Berlin, Germany), vol. 95, no. 4, pp. 431–444, 2017.

E. Esposito, I. Paterniti, R. Meli, P. Bramanti, and S. Cuzzocrea, “GW0742, a high-affinity PPAR-δ agonist, mediates protection in an organotypic model of spinal cord damage,” Spine, vol. 37, no. 2, pp. E73–E78, 2012.

V. Castagné, M. Cuénod, and K. Q. Do, “An animal model with relevance to schizophrenia: sex-dependent cognitive deficits in osteogenic disorder-Shionogi rats induced by glutathione synthesis and dopamine uptake inhibition during development,” Neuroscience, vol. 123, no. 4, pp. 821–834, 2004.