Proteomic Analysis of Lymphoblastoid Cells Derived from Monozygotic Twins Discordant for Bipolar Disorder: A Preliminary Study

Bipolar disorder is a severe mental illness characterized by recurrent manic and depressive episodes. In bipolar disorder, family and twin studies suggest contributions from genetic and environmental factors; however, the detailed molecular pathogenesis is yet unknown. Thus, identification of biomarkers may contribute to the clinical diagnosis of bipolar disorder. Monozygotic twins discordant for bipolar disorder are relatively rare but have been reported. Here we performed a comparative proteomic analysis of whole cell lysate derived from lymphoblastoid cells of monozygotic twins discordant for bipolar disorder by using two-dimensional differential in-gel electrophoresis (2D-DIGE). We found approximately 200 protein spots to be significantly differentially expressed between the patient and the co-twin (t test, p<0.05). Some of the proteins were subsequently identified by liquid chromatography tandem mass spectrometry and included proteins involved in cell death and glycolysis. To examine whether these proteins could serve as biomarkers of bipolar disorder, we performed Western blot analysis using case–control samples. Expression of phosphoglycerate mutase 1 (PGAM1), which is involved in glycolysis, was significantly up-regulated in patients with bipolar disorder (t test, p<0.05). Although PGAM1 cannot be regarded as a qualified biomarker of bipolar disorder from this preliminary finding, it could be one of the candidates for further study to identify biomarkers of bipolar disorder.


Introduction
Bipolar disorder is a severe mental illness characterized by recurrent manic and depressive episodes. It affects approximately 1% of the general population. In bipolar disorder, contributions of genetic and environmental factors have been indicated by family and twin studies [1]. Patients with this disease have a high relapse rate, and lifelong treatment with mood stabilizers such as lithium is often needed [2]. Past studies of bipolar disorder have focused on monoamines and intracellular signal transduction pathways as related to the effects of psychotropic drugs [3,4], but recent studies focus more on neuroplasticity or resilience based on reports of neuroprotective effects of mood stabilizers [5]. Recent genomewide association studies identified new candidate genes, such as CACNA1C [6,7], but the effects are relatively small.
Bipolar disorder is often misdiagnosed as depression or schizophrenia, and delayed diagnosis and treatment worsen the course of illness [8]. Thus, early diagnosis is important to prevent deterioration; however, no biomarkers for bipolar disorder are available yet. Identification of biomarkers would be indispensable for early diagnosis.
Much research on biomarkers in psychiatric diseases has been published [9,10]. Though several candidate biomarkers were identified from postmortem brains, these studies are undermined by many confounding factors, including cause of death, postmortem interval, and brain pH [11]. Postmortem interval affects degradation of mRNAs and protein levels [12,13,14,15] and significantly influences phosphorylation of signaling proteins [16]. In particular, the influence of medication cannot be ignored because a number of proteins are affected by mood stabilizers, antipsychotics, or antidepressant medication [17,18,19]. Additionally, large interindividual variations hamper the identification of biomarkers.
Biomarker research in other tissues such as serum, plasma, cerebrospinal fluids, saliva, and urine have also been performed [20,21] but has not led to a diagnostic test.
The concordance rate of bipolar disorder between monozygotic twins is approximately 70% [22,23]. By comparing monozygotic twins discordant for bipolar disorder, biochemical differences associated with bipolar disorder might be detected without interference from interindividual genetic variation [24].
Although several genome, transcriptome, and epigenome analyses in monozygotic twins discordant for bipolar disorder and other diseases have been reported [24][25][26][27][28][29], proteomic analysis has not been applied to identifying the difference between monozygotic twins discordant for bipolar disorder, possibly due to technical difficulties. Because transcript levels do not completely correlate with protein expression levels [30,31] and aberrant posttranslational modifications can cause disease, proteomic analysis is needed to supplement transcriptome and epigenome analyses.
In the present study, we performed proteomic analysis of lymphoblastoid cells derived from monozygotic twins discordant for bipolar disorder.

Subjects
For 2D-DIGE, we used lymphoblastoid cells derived from a pair of 42-year-old male monozygotic twins discordant for bipolar disorder. We initially examined their monozygosity by genotyping microsatellite repeat markers [25] and later confirmed it by single nucleotide polymorphism array [29].
For a case-control study of differentially expressed proteins using Western blot analyses, we used lymphoblastoid cells derived from eight unrelated patients with bipolar I disorder (BPI) (four men and four women, 35.669.0 years old [mean 6 SD], Japanese) and eight unrelated control subjects (six men and two women, 36.9610.0 years old [mean 6 SD], Japanese). Table 1 gives detailed information for each subject used in Western blot analysis. Patients and controls with a history of alcohol or illicit drug abuse were excluded from the study. The patients were treated with various medications. Diagnoses were made by the consensus of two senior psychiatrists using the Diagnostic and Statistical Manual of Mental Disorders, Fourth Edition (DSM-IV) criteria. Controls were selected from among students, nurses, office workers, and doctors in participating institutes, their friends, and other volunteers. A senior psychiatrist interviewed control subjects and found that they did not have major mental disorders. Written informed consent was obtained from all subjects. The ethics committees of RIKEN approved the study.

Cell culture and extraction of proteins
Lymphocytes were separated from the peripheral blood and transformed by Epstein-Barr virus using previously described standard techniques [32]. These cells were cultured and kept frozen until experiments. Lymphoblastoid cells were cultured in RPMI 1640 medium (Sigma-Aldrich, St. Louis, MO, USA) containing 10% fetal bovine serum (Sigma-Aldrich), 50 U/mL penicillin, 50 mg/mL streptomycin (GIBCO, Invitrogen/Life Technologies Corporation, Grand Island, NY), and 60 mg/mL tylosin solution (Sigma-Aldrich). Cells were cultured at approximately 1610 8 cells. Total proteins were extracted from the lymphoblastoid cells using the Q-proteome mammalian protein preparation kit (Qiagen, QIAGEN, Hilden, Germany). After the supernatant was precipitated with acetone, the pellet was dissolved using the lysis buffer (7 M urea, 2 M thiourea, 4% CHAPS, 10 mM Tris-HCl [pH 8.5]) (GE Healthcare Bio-Sciences, San Francisco, CA, USA).

Two-dimensional difference gel electrophoresis (2D-DIGE) and imaging analyses
Protein concentration was determined by the Bradford method using a protein assay kit (Bio-Rad, Hercules, CA, USA) and the Pierce 660 nm protein assay (Pierce, Rockford, IL, USA). Proteins were set to a final concentration of 5 mg/mL with the lysis buffer and labeled separately with 400 pmol of CyDye (Cy3 or Cy5) (GE Healthcare Bio-Sciences), vortexed, and incubated on ice in the dark for 30 min. A mixed sample composed of equal amount of proteins from the patient and the co-twin was labeled with Cy2 and used as an internal standard. After 30 min, the labeling reaction was stopped with 10 mM lysine. To avoid the possible effect of labeling efficiency, the dyes were swapped for each experiment using three gels. Labeled proteins were subjected to SDS-PAGE analysis, and the gels were scanned with the Typhoon 9400 scanner (GE Healthcare Bio-Sciences) at the wavelengths

Protein identification by mass spectrometry and database search
The preparative gels were stained with a SYPROH ruby (Invitrogen/Life Technologies Corporation) and scanned with the Typhoon 9400. Protein spots that showed differences in relative fluorescence were excised from the gel using the automated spot picker (GE Healthcare Bio-Sciences). The picked gel pieces were destained with 50% CH 3 CN in a 50 mM NH 4 HCO 3 solution. After removal of the supernatant, cysteine residues were reduced with dithiothreitol and carbamidomethylated with iodoacetamide. In-gel trypsin digestion was performed at 37uC overnight, using sequencing grade modified trypsin (Promega, Southampton, UK) reconstituted in 100 mM NH 4 HCO 3 . The trypsinized gel was rinsed three times in extraction buffer (5% trifluoroacetic acid in 50% CH 3 CN and 50% H 2 O). The trypsinized peptides solution was dried by speed vacuum, suspended in 2% CH 3 CN with 0.1% trifluoroacetic acid, and analyzed by LTQ (Fisher Scientific, Waltham, MA) liquid chromatography/linear ion trap mass spectrometry (LC-MS/MS) system. Their corresponding proteins were searched using the program Mascot database-searching software (Matrix Science, London, UK), which accesses protein identification by matching mass spectroscopy data with the protein databases NCBI (http://www.ncbi.nlm.nih.gov) and UniProt (http://www.uniprot.org/uniprot). Identification criteria included a Mascot score .45 (selected based on a corrected p-value ,0.05). Functional grouping of altered proteins The IPA software (Ingenuity Systems, Redwood City, CA) was used to identify the key biological relationships and functions of differentially expressed proteins and their interaction networks. For pathway analysis, Swiss-Prot and GenBank accession numbers were used to catalogue the identified protein into known interaction pathways based upon the Ingenuity Pathway Knowledge Base (IPKB). IPA classification data were derived from the published literature in a systematic way using a comprehensive ontology of functional annotations and protein-protein interaction data. The most significant interaction networks, biological functions and pathway associated with the differentially expressed proteins were identified. To confirm biological functions and gene ontology annotation between identified proteins, we used bioinformatics resources, PANTHER (http://www.pantherdb. org/) and DAVID (http://david.abcc.ncifcrf.gov/). DAVID program uses a modified Fisher's exact p value (EASE score) to rank gene clusters by statistical overrepresentation of individual genes, based on the co-occurrence/enrichment of the category within the gene list relative to all genes in the same category on the study.

Western blot analyses for validation
Lymphoblastoid cells in a case-control study were individually cultured. Total proteins were individually extracted using Qproteome Mammalian Protein prep kit (Qiagen), and protein  concentrations were determined by the methods mentioned above. Equal concentration (5 mg per lane) of proteins from control and case samples were separated by 12% or 4-15% SDS-polyacrylamide gel electrophoresis and transferred onto Immun-Blot PVDF membranes (Bio-Rad) using a mini Trans-blot Cell (Bio-Rad). After transfer, the blotted membrane were blocked with 4% w/v ECL Advance Blocking Reagent (GE Healthcare Bio-Sciences) in phosphate-buffered saline containing 0.1% Tween20 (PBST) (MP Biomedicals Inc., Santa Ana, CA) at 4uC overnight and incubated with primary antibody in PBST with 4% w/v skim milk for 1 h at room temperature. PSME1, WARS, and PGAM1 were chosen for quantification by Western blot analysis. As an internal control, NM23A was used because there was no significant difference of protein levels of NM23A between patient and the co-twin in 2D-DIGE. Primary antibodies were as follows: rabbit antibody against human PSME1 (Calbiochem, La Jolla, CA), mouse antibody against human WARS (Abnova, Taipei City, Taiwan), goat antibody against human PGAM1 (Novus Biologicals, Littleton, CO), and rabbit monoclonal and polyclonal against human NM23A (Abcam, Cambridge, MA). The blotted membranes were washed in TBST and incubated at room temperature with each secondary antibodies, Alexa Fluor 488-labeled donkey anti-goat IgG antibody (Molecular Probes/Life Technologies Corporation), Cy5-conjugated affiniPure donkey anti rabbit IgG antibody and Cy3-conjugated affiniPure donkey anti mouse IgG antibody (Jackson Immuno Research, West Grove, PA). The membranes were directly scanned with the Molecular Imager FX (Bio-Rad). Protein bands were analyzed to give a quantitative estimation of intensity change using the Quantity One Software (Bio-Rad) adapted to the Molecular Imager FX. To estimate the relative molecular weight of each protein, molecular markers, Dual Color Precision Plus Protein Standards (Bio-Rad) and ECL Plex Fluorescent Rainbow Markers (GE Healthcare Bio-Sciences) were used. Preliminary experiments indicated that amounts of these proteins in the lysates of lymphoblastoid cells were within the linear range of detection.

Detection of spots differentially expressed between monozygotic twins discordant for bipolar disorder by 2D-DIGE
First, we extracted total protein from cultured lymphoblastoid cells derived from a pair of monozygotic twins discordant for bipolar disorder. The total protein for each twin was separately labeled with different CyDyes (Cy3 or Cy5), and dyes were swapped between gels. A mixed sample composed of an equal amount of proteins from the patient and the co-twin was labeled with Cy2 and used as an internal standard. These processes minimized gel-to-gel variation and improved protein spot statistics at the analysis stage. These labeled proteins were mixed and analyzed by 2D-DIGE. To detect robust differences between the patient and the co-twin, we performed 2D-DIGE. To avoid artifacts, four gels of same condition were simultaneously run for each experiment. Three gels were used as analytical gels to detect differentially expressed spots between the patient and the co-twin, and the remaining gel was used as preparative gel for picking out differentially expressed spots. We performed 2D-DIGE and liquid chromatography tandem mass spectrometry (LC-MS/MS) anal-yses in quadruplicate experiments using protein samples independently extracted from different aliquots of cell culture.
Using 2D-DIGE and DeCyder Ver.6.0 image analysis software, approximately 3200 protein spots were separated (3220651 [mean 6 SD]). The observed spot pattern images in each gel were very similar among quadruplicate experiments. The protein spots were selected if the intensity difference between the affected and nonaffected twin was larger than 1.25-fold (absolute value .1.25 or absolute value ,21.25). Approximately 200 spots (2116115, quadruplicate) were found to be significantly differentially expressed (p,0.05, Student's t test) between the twins per experiment.

Identification of spots by LC-MS/MS
The proteins in preparative gels were stained by SYPRO Ruby after electrophoresis, and the preparative gels were performed matching to analytical gels. The largest 100 spots out of the 200 differentially expressed protein spots were picked from each preparative gel and were considered suitable for subsequent analysis by LC-MS/MS. The spots were chosen sequentially from those with a large absolute value fold change. The 68 spots, averaged for quadruplicate values, were successfully identified as unique proteins through LC-MS/MS and are listed in Table 2. Since Epstein-Barr virus-transformed lymphoblastoid cells were used in this study, all immunoglobulin and B-cell-related proteins were removed from the analysis. Moreover, we also removed keratin and trypsin-related proteins because of the possibility of experimental contamination. Table 2 shows only proteins identified by LC-MS/MS in common with each experiment. Proteins identified twice or more with Mascot search results (ion scores of higher than 45) included PSME1, RPLP0, TPI1, ALDOC, ANXA4, PGAM1, and WARS. Fifty-three proteins had high ion scores and were identified at least twice in four experiments.

Functional grouping of altered proteins
To explore the biological function (protein-interaction network) of the 53 differentially expressed proteins, we performed protein classification using the Ingenuity Pathway Knowledge Base software. A data set containing the gene symbols by The HUGO Gene Nomenclature Committee was uploaded into the application. Each protein identifier was converted to its gene identification and mapped to its corresponding gene object in the Ingenuity Pathway Knowledge Base. These genes were overlaid onto a global molecular network developed from information contained in the Ingenuity pathways analysis (IPA) knowledge base, which is based entirely on findings reported in the literature. Networks of these focus genes were algorithmically generated based on their connectivity. One of the 53 proteins was omitted because it was not included in the database, and the remaining 52 proteins were mapped onto mainly three networks. The largest network, having the highest score, is associated with carbohydrate metabolism, neurological disease, and skeletal and muscular disorders (Fig. 1A). This network includes ALDOA, ALDOC, ANXA4, ANXA5, CAPZB, CORO1A, DNAJB11, ENO1, GAPDH, HIST1H4A (includes others), HNRNPM, HSPB1, LDHA, NPM1, OTUB1, PDIA6, PGAM1, POLR2E, PRDX2, PSME1, RNA polymerase II, RPLP0, STMN1, VDAC1, and WARS. Next, the 53 proteins were classified according to biological function and canonical intensities. The p values were calculated using a t test in all proteins. Number of the subjects is 8 for bipolar disorder and 8 for controls, respectively. The absolute band intensity for the PGAM1 was also significantly higher in patients with bipolar disorder (0.93+/20.23 [mean +/2 standard deviation] [arbitrary unit]) than control subjects (0.39+/20.18, p,0.0005). doi:10.1371/journal.pone.0053855.g002 Figure 3. Verification of the alternation of PGAM1 using the other internal protein, tubulin alpha (TUBA). The methods are similar to the Figure 2, except that tubulin alpha (TUBA) was used for the internal standard. a) PGAM1, b) PSME1, c) WARS. Scatter plots show the ratio of each protein to an internal standard protein, TUBA, measured by densitometric scanning of the band intensities. The PGAM1/TUBA ratio was significantly pathway. The categories pertained to carbohydrate metabolism (8.38E-10-3.27E-02, 14 molecules) and cell death (2.04E-05-4.72E-02, 20 molecules) (Fig. 1B). Pathway analysis and gene ontology classification using PANTHER and DAVID were conducted on the same protein IDs. These analyses also showed pathways and categories associated with glycolysis and antiapoptosis. Taken together, most proteins identified in the present study were related to glycolysis and neurological diseases.

Case-control study by Western blot analysis
We postulated that the differentially expressed proteins might be candidate biomarkers for bipolar disorder. To validate the findings from the proteomic profiling study and to examine the possibility of biomarkers for bipolar disorder, Western blot analyses were performed using a case-control sample set consisting of eight subjects with bipolar disorder and eight healthy control subjects. To compare the protein levels across individuals, protein concentration was measured by the Bradford method, and equal amounts of proteins were loaded onto the gels. Commercially available antibodies for the candidate 7 proteins (PSME1, RPLP0, TPI1, ALDOC, ANXA4, PGAM1, and WARS) were searched for, and among available antibodies, those against PSME1, WARS, and PGAM1 showed good performance, and thus they were chosen for quantification by Western blot analysis. The levels of PGAM1, PSME1, and WARS were quantitatively investigated by Western blot analysis using NM23A as a standard (Fig. 2). Expression of PGAM1 was recognized by the presence of a single band at around 28 kDa and its protein expression was increased by 197% in bipolar disorder compared with controls (p,0.05). However, the levels of the other proteins were similar between bipolar disorder and controls in this case-control sample set (Fig. 2).
The absolute band intensity for the PGAM1 was also significantly higher in patients with bipolar disorder (0.9360.23 [mean 6 standard deviation] [arbitrary unit]) than control subjects (0.3960.18, p,0.0005). In addition, we also performed an independent experiment using the other, more popular housekeeping protein, tubulin alpha (TUBA), as an internal standard. This analysis also showed higher PGAM1 levels in patients with bipolar disorder than controls (p,0.05) (Figure 3).

Discussion
In this study, we identified 53 proteins that were differentially expressed between a pair of monozygotic twins discordant for bipolar disorder; 34 were up-regulated and 19 were downregulated. The differentially expressed proteins included those previously implicated in psychiatric disorders, such as ALDOC, ENO1, and PRDX2 [10,33]. Differences for ALDOC, ANXA4, PGAM1, PSME1, RPLP0, TPI1, and WARS between twins were regarded as robust because they were identified in three of four experiments with high scores.
To evaluate whether identified proteins might be biomarkers for bipolar disorder, we performed a case-control study for several proteins by Western blot analysis using available antibodies. An increased level of PGAM1 was observed in samples from patients with bipolar disorder. PGAM1 is an enzyme of the glycolytic pathway that catalyzes the conversion of 3-phosphoglycerate to 2phosphoglycerate [34]. This enzyme also promotes glycolysis and ATP production via the TCA cycle and the electron transport system. Although previous studies using postmortem brains of patients with bipolar disorder and schizophrenia suggested altered protein expression of glycolysis enzymes, including PGAM1 [10,35], the results were controversial. The differentially expressed proteins between bipolar disorder and healthy control including PGAM1, might be a clue to understand the biological basis of bipolar disorder.
To examine whether the 53 identified proteins were related to each other and constituted a global molecular network, pathway, or category, we applied IPA to our data. The results showed that the networks having a high score belonged to cell death, energy production, and glucose metabolism categories. The cell death category included the following proteins: NPM1, P4HB, LGALS3, CASP3, PDIA3, ATP5A1, GAPDH, ANXA4, HSPA5, RPLP0, UCHL1, STMN1, ENO1, ANXA5, MZB1, PSMB1, ALDOA, VDAC1, LDHA, HSPB1, and PRDX2 (Fig. 1). These results are consistent with previous studies. Benes et al. [36] showed increased expression of pro-apoptotic gene transcripts in postmortem brains of bipolar disorder patients. Furthermore, Herbeth et al. [21] indicated altered cell death and inflammation-related proteins in peripheral blood mononuclear cells and serum from patients with euthymic bipolar disorder. Brain imaging studies demonstrated reductions in the mean gray matter volume of brains from patients with bipolar disorder [37]. Previous studies reported a decreased density of nonpyramidal neurons in layer II of the anterior cingulate and a lower number of glial cells in layer III with bipolar disorder [38]. Meta-analyses of volumetric magnetic resonance imaging studies showed reduced volume of gray matter in the anterior cingulate and bilateral insula [39,40]. Neuropathological studies of bipolar disorder showed decreases of each brain field and neuronal cells. Because mood stabilizers and antidepressants, which are used for treatment of bipolar disorder, have neuroprotective actions [5,41,42], it has been suggested that cells derived from patients with bipolar disorder are more vulnerable to factors related to cell death than those from controls. Patients with unipolar or bipolar depression exhibit decreased brain-derived neurotrophic factor levels [43]. Moreover, mood stabilizers have neuroprotective effects by increasing bcl-2 levels [42,44,45]. These findings suggest cellular vulnerability has a role in the pathology of bipolar disorder. Dysregulation of the apoptotic process found in the monozygotic twins discordant for bipolar disorder might be relevant to this hypothesis.
We examined the relationship of the identified proteins with canonical pathways and found that the proteins were related to the glycolysis pathway. The proteins included PKM2, ALDH2, ENO1, PGAM1, GAPDH, ALDOA, LDHA, and ALDOC. Glycolysis, or anaerobic respiration, is a fundamental metabolic process that produces energy for all cells. In order to maintain its functions, the brain needs an enormous amount of energy compared with other tissues. ALDOC is a brain-specific glycolysis enzyme that catalyzes the reversible aldol cleavage of fructose-1,6biphosphate and fructose-1-phosphate to dihydroxyacetone phosphate and either glyceraldehyde-3-phosphate or glyceraldehyde [46]. In the present study, we found a decrease of the ALDOC protein level in the affected twin. However, previous reports showed that protein expression level of ALDOC was increased in the frontal cortex of patients, including those with mood disorder [35,47]. This discrepancy might reflect differences between tissues. Moreover, we found differential expression of many essential enzymes of glycolysis such as TPI1, ALDOA, and PGAM1. A higher in patients with bipolar disorder compared with controls (p,0.05). Number of the subjects is 8 for bipolar disorder and 8 for controls, respectively. doi:10.1371/journal.pone.0053855.g003 previous report using positron emission tomography showed that familial bipolar depressive patients had decreased blood flow in the cerebrum and a decreased rate of glucose metabolism in the ventral anterior cingulate cortex [37]. As indicated by an alteration in energy metabolism, compromised metabolic function has been reported in bipolar disorder [48,49]. In these studies, alteration of mitochondrial proteins was reported. Mitochondria are involved in processes including the TCA cycle, glycolysis and gluconeogenesis, lipogenesis, and malate-asparate shuttle [50]. Thus, changes in these proteins may lead to major alterations in the energy pathways, thus affecting ATP production. Recently, many reports have suggested that mitochondrial dysfunction is involved in bipolar disorder and other psychiatric disorders [51,52,53]. Mitochondria are also involved in other essential processes such as apoptosis, oxidative stress, and calcium regulation [50]. Thus, a decrease in energy production due to mitochondrial dysfunction in the brains of patients with bipolar disorder may be compensated for by an increase in energy production by glycolysis. It is possible that mitochondrial dysfunction affects neuronal cell death. Further study is needed to know whether these alterations in glycolysis-related proteins are a cause or consequence of the disease process. This is the first study to our knowledge to apply proteomics for the analysis of monozygotic twins discordant for bipolar disorder, and it has major limitations. First of all, we analyzed only a single pair of monozygotic twins. Thus, results cannot be applied to bipolar disorder in general. Another limitation is the tissue examined; that is, lymphoblastoid cells. Although brain samples may be optimal to identify molecules directly related to bipolar disorder, brain samples of twins are difficult to access. In addition, accessible tissues such as body fluid and peripheral cells such as serum, plasma, cerebrospinal fluids, saliva, urine, and peripheral blood cells should be used for biomarkers. In this study, we used lymphoblastoid cells and avoided a possible effect of medication by culturing the cells in drug-free media. However, a possibility that the effect of medication at the collection of blood last even after culturing the cells in drug-free media for a month cannot be totally ruled out. The other major limitation is the small number of casecontrol samples.
In summary, we performed a proteomic analysis of lymphoblastoid cells in a pair of monozygotic twins discordant for bipolar disorder. The identified proteins were mainly categorized as those involved in cell death and glycolysis. In a case-control study, protein expression of PGAM1, which is related to glycolysis, was significantly higher in patients than in healthy controls. The present findings suggest future new targets that may be relevant to the pathology of bipolar disorder. The present results need to be tested in a larger, independent sample set to reach a valid conclusion.