G protein-linked signaling pathways in bipolar and major depressive disorders

The G-protein linked signaling system (GPLS) comprises a large number of G-proteins, G protein-coupled receptors (GPCRs), GPCR ligands, and downstream effector molecules. G-proteins interact with both GPCRs and downstream effectors such as cyclic adenosine monophosphate (cAMP), phosphatidylinositols, and ion channels. The GPLS is implicated in the pathophysiology and pharmacology of both major depressive disorder (MDD) and bipolar disorder (BPD). This study evaluated whether GPLS is altered at the transcript level. The gene expression in the dorsolateral prefrontal (DLPFC) and anterior cingulate (ACC) were compared from MDD, BPD, and control subjects using Affymetrix Gene Chips and real time quantitative PCR. High quality brain tissue was used in the study to control for confounding effects of agonal events, tissue pH, RNA integrity, gender, and age. GPLS signaling transcripts were altered especially in the ACC of BPD and MDD subjects. Transcript levels of molecules which repress cAMP activity were increased in BPD and decreased in MDD. Two orphan GPCRs, GPRC5B and GPR37, showed significantly decreased expression levels in MDD, and significantly increased expression levels in BPD. Our results suggest opposite changes in BPD and MDD in the GPLS, “activated” cAMP signaling activity in BPD and “blunted” cAMP signaling activity in MDD. GPRC5B and GPR37 both appear to have behavioral effects, and are also candidate genes for neurodegenerative disorders. In the context of the opposite changes observed in BPD and MDD, these GPCRs warrant further study of their brain effects.


INTRODUCTION
The G-protein linked signaling system (GPLS) comprises Gproteins, G protein-coupled receptors (GPCRs), GPCR ligands, and downstream effector molecules. G-proteins interact with both GPCRs and downstream effectors such as cyclic adenosine monophosphate (cAMP), phosphatidylinositols, and ion channels. The GPLS is implicated in the pathophysiology and pharmacology of both major depressive disorder (MDD) and bipolar disorder (BPD). GPCRs are integral membrane proteins with key roles in numerous physiological cellular processes. GPCRs can signal ligand binding events by conformational changes that activate intracellular G proteins. Active G proteins in turn trigger intracellular downstream signaling pathways. Evidence suggests that altered GPCRs and their two major downstream signaling pathways, mediated by effectors such as cAMP and phosphatidylinositol, may be involved in the pathophysiology of BPD and MDD as well as in the mechanism of drug treatment for these disorders. Altered immunoreactivities and functional activities of proteins related to cAMP signaling pathways have been reported in the brain and peripheral blood cells of BPD patients (Jope et al., 1996;Karege et al., 1996;Gould and Manji, 2002;Chang et al., 2003). Also, mood stabilizers and antidepressants have been shown to affect cAMP and phosphatidylinositol signaling (Donati and Rasenick, 2003;Harwood, 2005;Xu et al., 2007;Racagni and Popoli, 2008). There are at least 210 GPCRs for which a natural ligand has been identified. There are 150 orphan GPCRs, with no known ligand or function (Wise et al., 2004), that might be relevant to downstream signaling dysregulation observed in BPD and MDD. Among the signal transmission systems associated with GPCRs, monoaminergic, and neuropeptidergic systems are believed to be dysregulated in BPD and MDD (Lieb et al., 2002;Elhwuegi, 2004;Carvajal et al., 2006;Ruhe et al., 2007).
Messenger RNA expression levels of a number of GPCRs, their ligands, and effector molecules involved in the GPLS pathways have been measured in post-mortem brains of BPD and MDD patients Spleiss et al., 1998;Caberlotto and Hurd, 2001;Kuromitsu et al., 2001;Zhou et al., 2001;Gurevich et al., 2002;Pandey et al., 2002;Agam et al., 2003;Chang et al., 2003;Dwivedi et al., 2004;Escriba et al., 2004;Lopez-Figueroa et al., 2004;Sherrin et al., 2004). Brain imaging studies have implicated anterior cingulate cortex (ACC) and the dorsolateral prefrontal cortex (DLPFC) in the pathophysiology of mood disorders (Harrison, 2002;Rogers et al., 2004;Frazier et al., 2005). In this study, the DLPFC and ACC were selected for microarray gene expression profiling of BPD and MDD patients focusing on the involvement of genes involved in GPLS pathways in these disorders and the effect of drug treatments. These cortical regions were selected because of the integral role of the ACC and DLPFC in the control of impulses, memory, learning, and hedonic evaluation, and depression. This study compares an a priori set of genes, in both mood disorders, using two cortical regions. There have been a number of previous studies comparing MDD and BPD by microarray analysis to controls, e.g., studies involving samples from the Stanley Foundation Neuropathology Brain collection (Kim and Webster, 2010). However, few studies have analyzed gene expression by microarray analysis of the Anterior Cingulate in both mood disorder groups, aside from studies involving the Pritzker Neuropsychiatric Disorders papers which involved selected pathways of glutamate, mitochondria, and growth factors Choudary et al., 2005;Vawter et al., 2006). We hypothesized that there would be dysregulation of the GPLS in both mood disorders compared to controls, and that some overlap between the signature of each mood disorder would be found, given that most of our subjects were in a depressive state near the time of death.

SUBJECTS
The primary cohort consisted of 22 subjects, including seven healthy controls, six patients with bipolar disorder type I (BPD), and nine patients with MDD (Cohort A). The MDD results were validated with an additional independent cohort, including seven controls and five MDD patients (Cohort B). Toxicology data and/or recent clinical data were examined for the BPD and MDD subjects and detailed clinical information is summarized in Table 1. To minimize potential confounding effects of agonal events, we used three criteria in subject selection: brain tissue pH of >6.6, ribosomal RNA 28S/18S ratio of >1.5, and sudden death with the agonal period lasting less than 1 h. We have previously found large differences in gene expression values when using longer agonal periods, >1 h and as long as 1-2 days, and would have found large effects due only to mixing samples with long and short agonal duration Tomita et al., 2004;Vawter et al., 2006;Atz et al., 2007).
Brain tissue and clinical information were obtained by the UCI Brain Bank staff with the informed consent of decedents' next-ofkin, and were processed through a standard protocol approved by the local UCI Institutional Review Board . None of the subjects included in this study had specific agonal conditions including cancer, hypoxia, coma, pyrexia, seizure, dehydration, hypoglycemia, multi-organ failure, skull fracture, ingestion of neurotoxic substances or prolonged agonal duration, which are known to affect tissue pH, RNA integrity and gene expression profiles in post-mortem brain Tomita et al., 2004). The variables for gender, age, post-mortem interval (PMI), tissue pH, and 28S/18S ratio were balanced between diagnostic and control groups. Patients who were prescribed lithium and selective serotonin reuptake inhibitor (SSRI) at the time of death are also indicated ( Table 1).

MICROARRAY GENE EXPRESSION
The experimental procedures are described elsewhere Tomita et al., 2004). Total RNA was extracted from ACC and DLPFC of Cohort A (Table 1), and purified with silicabased mini-spin columns (Qiagen, Valencia, California). The oligonucleotide microarray experiments were conducted utilizing Affymetrix U133A Gene Chips following the manufacturer's protocol (Affymetrix, Santa Clara, CA). To detect gene expression differences between diagnostic groups reliably, we replicated the experiments as follows: Cohort A used 22 subjects' total RNA from ACC which were processed at two laboratories (University of California, Irvine and University of California, Davis) and DLPFC that were processed independently (University of California, Davis and University of Michigan, Ann Arbor). For Cohort B 12 total RNA samples from ACC and DLPFC (additional 12 subjects) were processed at the same independent laboratories as cohort A. In summary, for each BPD subject there were duplicate arrays run, and for MDD subjects there were duplicate arrays run, plus an additional cohort of MDD subjects were run also with duplicate arrays.

DATA ANALYSES
Since a considerable number of probe sequences in the Affymetrix original chip definition file (CDF) were found to incorrectly BLAST to other transcripts, a refined CDF based on an UniGene cluster was developed Dai et al., 2005), and was used for data analyses presented in this paper. The refined CDF file is available on http://brainarray.mhri.med.umich.edu/brainarray/.
The focus of this study was 445 GPLS genes found on the Affymetrix arrays used in this study (Table S1). Signal intensities were extracted with Robust Multi-array Average (RMA) for each probe set and each subject. Gene-wise Pearson's correlation coefficients between experimental duplicates were calculated. Only the genes that were significantly correlated (p < 0.05) between experimental duplicates were considered to be reliably measured genes, and were subjected to downstream analyses. These genes were analyzed in a mixed-model ANOVA utilizing Partek Pro 6.0 (Partek, MO) for the main effect of the diagnostic classification (BPD, MDD, control) and for sex . We calculated the false discovery rate (FDR) in the microarray study by Benjamini-Hochberg step-down procedure using Partek Genomics setting the FDR to 0.05. Since the three diagnostic groups were well matched on age and pH (Table 1), the gene expression values were not tested by ANCOVA. However, in posthoc analysis, we determined whether the tissue pH was correlated with any differentially expressed GPLS genes. The post-hoc analysis showed that there were no significant correlations found with tissue pH and differentially expressed GPLS genes. A post-hoc group comparison of subjects with SSRI treatment and without treatment was conducted in MDD subjects, and similar analysis for lithium was conducted for BPD subjects. The hypergeometric distribution was used for calculating probabilities that the observed number of genes in each gene category for cAMP or PI signaling pathways (listed in Table S2) were detected as differentially expressed genes in the comparisons between patients and control. False discovery rate were evaluated by QVALUE (Storey and Tibshirani, 2003). Calculations for hypergeometric probabilities were performed using SISA online statistics package (http://www.quantitativeskills.com/sisa/ distributions/hypergeo.htm).

QUANTITATIVE PCR
For further technical evaluation of the microarray data, we evaluated mRNA expression levels by quantitative real-time reverse transcriptase PCR (qPCR) for the seven genes listed in Table 3. We followed up on results primarily that met FDR in cohort A for BPD, and both Cohorts A and B for MDD, except for RGS20, which we followed up based upon Cohort A only. These genes selected for qPCR analyses met the criteria of FDR at the level of accepting 5% false positives and percentage change greater than 20% in the ACC microarray results. The threshold cycle (C t ) values for the target transcripts were normalized to averaged C t values for two reference transcripts, Jagged 1 (JAG1) and solute carrier family 9 isoform 1 (SLC9A1), which showed equivalent expression levels among the three diagnostic groups (BPD, MDD, and control) throughout the 3 brain regions on our normalized microarray data. After DNase digestion and purification, total RNA (1 µg) samples from each of the 34 subjects in Cohorts A and B were used as a template for first-strand cDNA synthesis using poly-dT primer. The mRNA for each transcript was measured using the SybrGreen system with Prism model 7000-sequence detection instrument (Applied Biosystems, Foster City, California) and primer sequences listed in Table S1. The C t , which correlatedinversely with the target mRNA concentration, was measured as the cycle number at which the SybrGreen fluorescent signal increased abovea pre-set threshold level. The qPCR experiments were performed three times for each transcript, and the C t triplicate values were averaged.

In situ HYBRIDIZATION
The gene expression fold changes for GPRC5B and GPR37 showed an increase in BPD and decrease in MDD. We selected those genes for in situ hybridization (ISH) using the ACC regions from six representative subjects to evaluate the cellular localization of GPRC5B and GPR37 mRNA expression. A formal statistical analysis of the ISH data was not applied due to low subject number. ISH was performed following a protocol as described elsewhere (Neal et al., 2001;Lopez-Figueroa et al., 2004). Briefly, frozen ACC tissue blocks were cut into 12-µm sections on a cryostat. After fixation with 4% paraformaldehyde for 1 h, sections were hybridized overnight with a 35 S-labeled antisense cRNA probes for GPRC5B and GPR37. Following hybridization, wash, RNase A treatment, dehydration and air-dry, sections were exposed to Kodak XAR-5 X-ray film for 21 days. Technical control studies were performed by (1) using 35 S-labeled sense RNA probes for GPRC5B and GPR37, and (2) pre-treating tissue with RNase A before hybridizing with the antisense probes.

RESULTS
Significantly dysregulated GPLS genes that passed FDR (p < 0.05) in either Cohort A or B and either mood disorder, BPD or MDD, are listed in Table 2. For completeness, the uncorrected p-values are also shown, as some genes only met FDR for a single comparison. In general, equal numbers of GPLS genes passed FDR in BPD and MDD in the ACC (10 genes for BPD and 10 genes for MDD that passed FDR), compared with DLPFC (1 gene for BPD and 5 genes for MDD). Thus, about three fold more gene expression differences in mood disorder were found in the ACC showed compared to the DLPFC. The genes in Table 2 also showed non-significant correlations between gene expression and tissue pH (all p-values > 0.1) suggesting that residual pH-agonal effects on the expression levels of GPCR genes did not remain in the data set, even though complete elimination of these factors may not be possible with matching of subject's age and pH, and inclusion of subjects that had rapid deaths.

CYCLIC AMP SIGNALING PATHWAY
The genes involved in cAMP signaling pathways, which were differentially expressed in ACC of BPD patients compared to controls, are listed in Table 2. The cAMP signaling pathway is shown in Figure 1A. The data show that six significant transcripts out of 39 that repress cAMP signaling activity (Table 1) were increased in ACC from BPD patients at the transcriptional level. The group of molecules that act negatively on cAMP signaling activity was significantly over-represented in BPD based on a hypergeometric distribution (p = 0.037, q = 0.36). Thus, there was a nominal enrichment of genes that act negatively on cAMP signaling activity (G i proteins, phosphodiesterase, protein kinase A inhibitors, and protein phosphatase) in BPD in the ACC.   In contrast, of those genes that repressed cAMP signaling activity in MDD ( Figure 1B and Table 2) there was only one gene (PDE8A) differentially expressed, hence this signaling pathway was not over-represented (p > 0.05).

G s /G i -coupled GPCRs
Significant increases of NPY and NPY1R mRNA expression were found in BPD in ACC. GRM3 expression was also increased in ACC of BPD patients. G i /G 12 -linked EDG2 was increased in ACC of BPD, while G i -linked EDG1 was decreased in ACC of MDD. EDG1 and EDG2 regulate neuronal and glial functions, including neurite retraction, neurogenesis, and axonal guidance (Takuwa et al., 2002).

G proteins
Our data showed an increase of GNAI1 in the ACC of BPD, but no change in GNAI1 expression in the DLPFC of BPD. Consistent with our findings in DLPFC, immunoreactivity and functional activity of G i protein is not altered in the frontal cortex of BPD (Young et al., 1993;Friedman and Wang, 1996). RGS20 was significantly decreased in ACC of MDD cases.

Adenylate cyclase and protein kinase A
Expression levels of transcripts coding adenylate cyclases and cAMP-dependent protein kinases (PKAs) were not altered in BPD and MDD compared with controls in our study, which is consistent with a previous study (Chang et al., 2003). The reported

Frontiers in Genetics | Behavioral and Psychiatric Genetics
December 2013 | Volume 4 | Article 297 | 6 dysregulation in functional activities of adenylate cyclase and PKA (Gould and Manji, 2002) may not be due to alterations at the transcript level in the brain cortices.

Phosphatidylinositol signaling pathway
The phosphatidylinositol (PI) signaling pathway related genes which were differentially expressed in ACC of BPD patients are shown in Figure 1C and are summarized Table 2. The group of phosphatidylinositol 3-kinases was significantly over-represented in BPD based on a hypergeometric distribution (p = 0.042; q = 0.36). The PI signaling pathway related genes which were differentially expressed in ACC of MDD patients compared with controls are summarized in Figure 1D and Table 2. Since antidepressants have been reported to affect PI signaling Quintero et al., 2005), the potential effect of drug treatment on the gene expression was considered. In cohort A, five MDD subjects were prescribed SSRI at the time of death, whereas remaining four MDD subjects were not prescribed and presumably not taking any SSRIs (Table 1). There were no significant differences in mRNA expressions of PKCB1, ITPR1 and INPP5F between SSRItreated MDD and non-SSRI-treated MDD groups. However, ITPKB mRNA expression was significantly lower in the non-SSRI-treated MDD group than the SSRI-treated MDD group, which suggests that SSRIs may attenuate the altered expression in MDD toward levels observed in controls.

G q -coupled GPCRs
NTSR2 was decreased in ACC of BPD. NTSR2 mRNA and NTSR2 binding was reported to be down regulated in transgenic mice over expressing corticotrophin releasing factor, which is linked to anxiety, stress and depression (Peeters et al., 2004). EDNRB was also decreased in ACC of MDD. Endothelin is known to be a vasoconstrictor peptide as well as a neuromodulator regulating neuronal excitability. The data did not show any significant differences in expression levels of G q -linked monoamine GPCRs (serotonin-2 receptors and adrenergic alpha-1 receptors), which was consistent with previous microarray studies (Iwamoto et al., 2004;Konradi et al., 2004;Sibille et al., 2004;Aston et al., 2005).

Lithium and second messengers
The microarray data highlighted several target molecules of lithium treatment (Figure 1C). Lithium is known to block INPP1 (inositol polyphosphate-1-phosphatase) enzyme activity, and INPP1 mRNA expression was increased in ACC of BPD. Among the 6 BPD subjects, 3 BPD subjects treated with lithium at the time of death showed lower signal intensities for INPP1 and PIK3C2B compared with remaining nonlithium treated BPD subjects, although the difference did not reach significant changes. There is a modest genetic association between INPP1 and sensitivity to lithium treatment (Steen et al., 1998;Serretti, 2002). A second activity of lithium is activation of phosphatidylinositol 3-kinases and regulation of glycogen synthesis kinase 3 beta (GSK3B) activity (Sinha et al., 2005).
For MDD subjects, the mRNA expression of PKCB1, ITPR1, and INPP5F were increased in ACC, while ITPKB mRNA was decreased ( Figure 1D). PKCB1 and ITPR1 are receptors for diacylglycerol (DG) and inositol 1, 4, 5-triphosphate (ITP), respectively, both of which are intracellular second messengers produced by PLC-beta through a G q protein-dependent mechanism, while INPP5F and ITPKB metabolize ITP in a signal-terminating reaction. Previous reports suggest decreased PLC-beta activity in depressed patients (Pandey et al., 1997;Frey et al., 1998;Moore et al., 1999), and over expressed PKCB1 may down regulate PLC-beta activity (Pachter et al., 1992). The alterations in expression levels of the molecules interacting with DG and ITP may reflect functional impairments of G q and PLC-beta in MDD.

ORPHAN G PROTEIN-COUPLED RECEPTORS
The two GPCRs that showed the most consistent differential expression patterns through both mood disorder analyses were GPCR family C, group 5, member B (GPRC5B) and G proteincoupled receptor 37 (GPR37)as shown in Table 2 and validated by QPCR (Table 3). GPRC5B was significantly increased in ACC of BPD patients and in DLPFC of BPD. GPR37 was also significantly increased in ACC of BPD. GPRC5B was significantly decreased in ACC and DLPFC of MDD patients.

REAL TIME QUANTITATIVE PCR
We followed up on results primarily that met FDR in cohort A for BPD, and both Cohorts A and B for MDD, except for RGS20, which we followed up based upon Cohort A meeting FDR while Cohort B met uncorrected p-value significance only. Nine GPLS transcripts were selected for validation, five for BPD and four for MDD. The qPCR results showed 100% concordance between qPCR and microarray fold change for BPD in ACC (Table 3) and 4/5 transcripts were statistically significant. Four GPLS related transcripts showed 100% concordance between qPCR and microarray fold change for MDD in ACC ( Table 3) and 4/4 transcripts were statistically significant. NPY, GPRC5B, GPR37, and INPP1 were significantly increased in ACC of BPD, whereas GPRC5B, GPR37, RGS20, and PPP1R3C were significantly decreased in the ACC of MDD.

In situ HYBRIDIZATION HISTOCHEMISTRY
GPR37 mRNA was preferentially expressed in subcortical white matter. The sense probe hybridization signal showed background signal. GPR37 mRNA expression in subcortical white matter was higher in ACC of both of the BPD subjects compared to both of the control subjects, whereas GPR37 mRNA expression in subcortical white matter was lower in ACC of both of the MDD subjects compared to both of the control subjects (Figure 2). This preliminary data suggests that the significant differential microarray and QPCR expressions of GPR37 between BPD/MDD and controls may reflect altered expression levels of GPR37 mRNA in subcortical white matter. mRNA expression of GPRC5B was not specific in ACC samples tested, most likely due to strong background hybridization issues.

DISCUSSION
The results of the study of GPLS transcriptional profiling showed differences in BPD and MDD in both the cAMP-and www.frontiersin.org December 2013 | Volume 4 | Article 297 | 7  phosphatidylinositol signaling pathways. The following pattern of alterations in GPLS were observed: (1) Transcriptional levels of molecules acting as negative regulators in cAMP signaling pathways are increased in BPD and decreased in MDD, which may compensate for the activated cAMP signaling activity in BPD and the blunted cAMP signaling activity in MDD, (3) Expression of INPP1 and phosphatidylinositol 3-kinases were altered in BPD suggesting that these molecules may be involved in both the pathophysiology and the mechanism(s) of action of lithium treatment of BPD, and (4) Two orphan GPCR genes, GPRC5B and GPR37, showed significant increases in ACC of BPD, and decreases in ACC and DLPFC of MDD.

ORPHAN G PROTEIN-COUPLED RECEPTORS
Among all GPCRs, the most consistent differential expression patterns were observed for GPRC5B and GPR37. GPRC5B was significantly increased in ACC and DLPFC of BPD, while significantly decreased in ACC and DLPFC of MDD patients. GPR37 was significantly increased in ACC of BPD, and significantly decreased in ACC and DLPFC of MDD patients. The consistent cortical alterations of GPRC5B and GPR37 expression observed in MDD patients suggest a role of these genes in the pathophysiology of MDD. So far, these two orphan GPCRs have not been highlighted as candidate susceptible genes for BPD and MDD, however previous microarray data supports decrease of both GPRC5B and GPR37 expressions in MDD patients (Aston et al., 2005). Our microarray data strongly suggests these genes might be involved in the pathophysiology of the disorders. GPR37 and GPRC5B are both implicated in Parkinson's disease (Marazziti et al., 2009(Marazziti et al., , 2011, which has depression as one premorbid indicator, possibly concomitant with brain neurodegenerative processes. GPR37 appears to facilitate dopamine neurotransmission, and reductions in brain GPR37 leads to behavioral alterations in conditioned place preferences to reinforcements such as cocaine or amphetamine (Marazziti et al., 2011). GPR37 (Class A receptor) contains a large extracellular domain which led to speculations that the ligand of this receptor may be a protein rather than a peptide (Leng et al., 1999). GPR37 is predominantly expressed in the brain (Zeng et al., 1997;Marazziti et al., 1998;Leng et al., 1999), especially in glial cells of the fiber tracts, Purkinje cells in cerebellum, neuronal cells in hippocampus CA3 region and substantianigra (Zeng et al., 1997;Imai et al., 2001). GPR37 has been shown to be involved in the dopaminergic system (Imai et al., 2001(Imai et al., , 2002Marazziti et al., 2004). GPR37 is also called Pael receptor (PAELR). Unfolded PAELR is a substrate of E3 ubiquitin ligase Parkin, and accumulation of Pael-R in the endoplasmic reticulum (ER) of dopaminergic neurons is considered to induce ER stress leading to neurodegeneration in Parkinson disease (Imai et al., 2002). Gpr37 −/− mice shows reduced striatal dopamine content and enhanced amphetamine sensitivity (Marazziti et al., 2004). Our ISH data suggest that the predominant expression of GPR37 mRNA in ACC occurred mainly in subcortical white matter. GPR37 is a substrate of parkin (PARK2), and its insoluble aggregates accumulate in brain tissue samples of Parkinson's disease patients, AKA PaelR (Marazziti et al., 2009).
GPRC5B is phylogenetically classified into a type C GPCR, however, in contrast to other type C genes, GPRC5B has a short N terminal domain (Brauner-Osborne and Krogsgaard-Larsen, 2000). QPCR based expression analysis shows GPRC5B is predominantly expressed in brain, with most abundant expression in neuronal cells of caudate, putamen, substantia nigra, thalamus and hippocampus and glial cells of the fiber tracts (Brauner-Osborne and Krogsgaard-Larsen, 2000;Robbins et al., 2002). Behavioral abnormalities exist in GPRC5B knock-outs, decreased spontaneous locomotor activity, and decreased responses to a new environment (Sano et al., 2011).
The expression patterns of GPRC5B and GPR37 suggest that both of these orphan GPCRs may be involved in various cellular functional activities in the brain including neuronal and glial cells. Functional characterization of GPRC5B and GPR37 may provide new insights into the pathophysiology of mood disorder, and potential novel strategies for developing methods for diagnoses and treatment of the disorders.
There are limitations to this post-mortem study of GPLS pathway expression. We conducted seven qPCR validations of genes that passed FDR, and although a small number, those genes did represent large enough percent fold changes to validate by qPCR. Other technologies could be used that are less costly, such as Nanostring or Fluidigm for additional validation studies, which are now being conducted on larger sample sizes. In addition, we do not have confirmed toxicological analysis of all drugs of abuse and therapeutic levels to rule out these effects in explaining case-control expression differences. Finally, we are aware that small residual agonal effects might exist that are not corrected in this analysis, however we did confirm that the present data presented in Table 2, was not due to the effects of age or pH on gene expression.

SUMMARY
This study highlights transcript profile differences in GPLS pathways in both BPD and MDD. Since sample sizes were low, we also used a second independent cohort in the case of MDD gene expression to check concordance. We relied upon duplicate microarray analysis at independent laboratories, and used gene expression values that were significantly correlated between sites, and employed FDR to select genes to advance for qPCR experimental validations.
The disturbances observed in cAMP-and phosphatidylinositol-pathway regulation may be due to differential activation of individual GPCRs and individual differences in ligand-receptor sensitivity, suggesting that functional and genetic characterization of these candidate genes may provide further understanding of the pathophysiology of mood disorders. Further studies will also be required to investigate additional GPCR signaling pathway features such as receptor phosphorylation, beta-arrestin recruitment, receptor desensitization, and internalization, and beta-gamma mediated signaling. The results of the present study suggest that future functional activity studies of the GPLS pathways would elucidate the relationship of the observed transcript changes to protein levels and functional outcomes in mood disorders.

ACKNOWLEDGMENTS
Authors of this manuscript were members of the Silvio O. Conte Center for "Genomic Studies in Bipolar and Major Depression" which was supported by the NIMH grant P50 MH60398, and the Pritzker Neuropsychiatric Disorders Research Consortium currently supported by the Pritzker Neuropsychiatric Disorders www.frontiersin.org December 2013 | Volume 4 | Article 297 | 9