Email updates

Keep up to date with the latest news and content from BMC Genomics and BioMed Central.

Open Access Highly Accessed Research article

Suicide candidate genes associated with bipolar disorder and schizophrenia: An exploratory gene expression profiling analysis of post-mortem prefrontal cortex

Sanghyeon Kim1, Kwang-Ho Choi2, Ali Fuat Baykiz3 and Howard K Gershenfeld14*

Author Affiliations

1 Department of Psychiatry, Univ. of Texas Southwestern Medical Center, Dallas, Texas 75390-9070, USA

2 Stanley Foundation Laboratory of Brain Research, Department of Psychiatry, Uniformed Services University of Health Sciences, Bethesda, MD 20892, USA

3 Elazig Asker Hastanesi Psikiyatri Klinigi 23300 Elazig, Turkey

4 *Department Integrative Biology, Univ. of Texas Southwestern Medical Center, Dallas, Texas 75390-9070, USA

For all author emails, please log on.

BMC Genomics 2007, 8:413  doi:10.1186/1471-2164-8-413


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2164/8/413


Received:11 January 2007
Accepted:12 November 2007
Published:12 November 2007

© 2007 Kim et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Suicide is an important and potentially preventable consequence of serious mental disorders of unknown etiology. Gene expression profiling technology provides an unbiased approach to identifying candidate genes for mental disorders. Microarray studies with post-mortem prefrontal cortex (Brodmann's Area 46/10) tissue require larger sample sizes. This study poses the question: to what extent are differentially expressed genes for suicide a diagnostic specific set of genes (bipolar disorder vs. schizophrenia) vs. a shared common pathway?

Results

In a reanalysis of a large set of Affymetrix Human Genome U133A microarray data, gene expression levels were compared between suicide completers vs. non-suicide groups within a diagnostic group, namely Bipolar disorder (N = 45; 22 suicide completers; 23 non-suicide) or Schizophrenia (N = 45; 10 suicide completers ; 35 non-suicide). Among bipolar samples, 13 genes were found and among schizophrenia samples, 70 genes were found as differentially expressed. Two genes, PLSCR4 (phospholipid scramblase 4) and EMX2 (empty spiracles homolog 2 (Drosophila)) were differentially expressed in suicide groups of both diagnostic groups by microarray analysis. By qRT-PCR, PLSCR4 and EMX2 were significantly down-regulated in the schizophrenia suicide completers, but could not be confirmed in bipolar disorder.

Conclusion

This molecular level analysis suggests that diagnostic specific genes predominate to shared genes in common among suicide vs. non-suicide groups. These differentially expressed, candidate genes are neural correlates of suicide, not necessarily causal. While suicide is a complex endpoint with many pathways, these candidate genes provide entry points for future studies of molecular mechanisms and genetic association studies to test causality.

Background

Suicide is the eleventh leading cause of death for all Americans with an age-adjusted annual rate of 10.5 per 100,000 in 2003 [1]. More than 90% of suicide completers have a psychiatric disorder and mood related disorders are the most common disease associated with suicide [2,3]. Patients suffering with bipolar disorder and schizophrenia have greatly increased rates of suicide with approximately 10% of patients dying of suicide [4,5]. Bipolar disorder and schizophrenia share common risk factors for suicide completion such as depression, previous suicide attempts, hopelessness, substance abuse, agitation, and poor adherence to treatment [4,5]. Suicide is a complex endpoint with many factors and pathways leading to death [6]. The hypothesis of a shared causation for suicide suggests common pathways and genes may function as susceptibility factors in both disorders. Alternatively, there could be specific distinct pathways within a diagnostic group.

Microarray technology provides an unbiased approach to the molecular causes of psychiatric disorders by examining the gene expression profile of cases vs. controls [7,8]. Recent microarray studies identified differentially expressed genes between suicide and depression patients vs. normal controls [9-11]. However, due to the small magnitude of the differential gene expression, the genetic heterogeneity of these mental disorders, and the mixed cellular nature of the brain tissue available [12], microarray studies with small sample sizes are prone to generate many false positive results [13]. Analysis of larger data sets pooled from independent studies increase the statistical power to find differentially expressed genes with small effect sizes in microarray studies. Recently, a large microarray data set generated by the Stanley Medical Research Institute has become available online (SMRIDB) [14]. This database contains clinical information and microarray data from 12 independent studies with post-mortem brain tissues of depression, bipolar disorder, schizophrenia, and unaffected control cohorts [15]. In this study, we reanalyzed this large microarray data set of bipolar disorder and schizophrenia patients. The question posed is, "to what extent are differentially expressed genes for suicide specific to diagnosis (bipolar disorder vs. schizophrenia) vs. a shared common pathway?"

Results

Suicide candidate genes in bipolar disorder and schizophrenia

Chi square tests of association indicated no difference in demographic variables between suicide vs. non-suicide subgroups within bipolar disorder. In contrast, age and smoking showed significant differences with suicide vs. non-suicide in the schizophrenia subgroup. Also, brain pH and sex significantly affected the expression levels of the differentially expressed genes between suicide vs. non-suicide groups within schizophrenia (Table 1). As shown in Table 1, none of these variables met the two criteria for inclusion as covariates in the two disorders (see methods). Therefore, we used no covariates in the omnibus model to provide a generalizable, single, same model for both schizophrenia and bipolar disorder. Between the suicide vs. non-suicide groups within bipolar disorder, a total of 13 genes were differentially expressed (Table 2). Among these genes, 10 genes were down-regulated and 3 genes including gamma-amino butyric acid A receptor, α5 subunit (GABRA5) were up-regulated.

Table 1. Demographic factors of the suicide group and the non-suicide group in bipolar disorder and schizophrenia

Table 2. Summary of differentially expressed genes between the suicide completer group vs. non-suicide group in bipolar disorder with fold change relative to non-suicide group.

Between the suicide vs. non-suicide group within schizophrenia, 70 genes were differentially expressed (Table 3). Most of these genes were down-regulated. From the above lists of differentially expressed genes, within diagnostic groups, two genes overlapped (Fig. 1A). Specifically, the phospholipid scramblase 4 (PLSCR4) and empty spiracles homolog 2, Drosophila (EMX2) genes were down-regulated in both suicide groups compared to the non-suicide groups. As negative controls, the normalization control probe set of 100 genes were tested by the same ANOVA model, and no genes met our statistical criteria (Fold Change ≥ |1.3| and FDR < 0.1) between the suicide group vs. non-suicide group within bipolar or schizophrenia diagnostic categories.

Table 3. Summary of differentially expressed genes between the suicide completer group vs. non-suicide group in schizophrenia with fold change relative to non-suicide group.

thumbnailFigure 1. Differentially expressed genes between suicide completer vs. non-suicide group in bipolar disorder and schizophrenia cohorts. (A) Venn diagram of differentially expressed genes between suicide completer vs. non-suicide groups within bipolar and schizophrenia. (B) Mean expression levels of PLSCR4 (phospholipid scramblase 4) and EMX2 (empty spiracles homolog 2 (Drosophila)) mRNA transcripts were determined for suicide completers vs. non-suicide groups within both bipolar disorder and schizophrenia in prefrontal cortex by qRT-PCR. The bars represent mean ± SEM. In schizophrenia, the mean expression levels of both genes were significantly down-regulated in the suicide group (n = 5) relative to the non-suicide cases (n = 25) by one-tailed, t-tests for unequal variances (EMX2 t(9) = 2.42, p = 0.02; PLSCR4 t(18) = 3.77, p = 0.0005). The estimated fold changes in the suicide group were -1.51 for EMX2 and -2.16 for PLSCR4 relative to the non-suicide group, consistent with our microarray data. In the bipolar disorder samples, no significant differences in mean expression levels for either gene were found between the suicide (n = 14 for EMX2; n = 12 for PLSCR4) vs. non-suicide cases (n = 15 for EMX2 ; n = 11 for PLSCR4) by unequal variance t-tests. * p < .05 ; *** p < .001.

Real time-PCR (RT-PCR) tested the validity of these two shared, differentially expressed genes in the available subset of the microarrayed samples. In the schizophrenia suicide group, the EMX2 and PLSCR4 expression levels were significantly down-regulated by comparison of mean expression levels (EMX2 t(9) = 2.42, p = 0.02; PLSCR4 t(18) = 3.77, p = 0.0005) in the suicide group compared to the non-suicide group (Fig 1B). The estimated fold changes in the suicide group were -1.51 for EMX2 and -2.16 for PLSCR4 relative to the non-suicide group. These differences were consistent with our microarray data. However, in bipolar disorder, these two genes could not be validated in a subset of tissues available from the original microarray study patients (Fig. 1B). These results highlight the overall findings that few common differentially genes for suicide vs. non-suicide exist between diagnostic groups.

Biological process in the differentially expressed genes

Functional annotation of the differentially expressed genes by Gene Ontology indicated that 9 biological processes were significantly overrepresented (at level 4; P < 0.05) among the suicide candidate genes in the schizophrenia cohort (Fig. 2). The transport genes were the most commonly over-represented biological process in our suicide candidate gene list for schizophrenia, including the glial, high affinity, glutamate transporter (SLC1A3). In contrast, no significantly over-represented biological process group emerged with suicide candidate genes in the bipolar disorder group.

thumbnailFigure 2. Distribution of gene ontology groups for biological process at level 4 of differentially expressed genes between suicide vs. non-suicide in schizophrenia cohorts.

Additionally, we analyzed the list of differentially expressed genes for each diagnosis by the Ingenuity Pathways Analysis (IPA) software to identify biological pathways and networks. We identified distinct signaling networks from suicide candidate genes that included the EMX2 gene in both disorders (Fig. 3). In bipolar disorder, the pathway perspective suggested a signaling network related to both cellular movement and cell to cell signaling, with interactions encompassing 10 differentially expressed, suicide candidate genes (Fig. 3A). By contrast, in schizophrenia patients, the differentially expressed genes were related in a cell death signaling network (Fig. 3B).

thumbnailFigure 3. Ingenuity pathway analysis of EMX2 and differentially expressed suicide candidate genes in bipolar disorder (A) and schizophrenia (B).

Discussion

In this re-analysis study, we explored suicide candidate genes associated with bipolar disorder and schizophrenia using an unbiased genome-wide expression profiling strategy. To identify suicide specific effects on the expression level of genes, we compared a suicide completers group to a non-suicide group within the same diagnostic category. The most important finding of this study is the identification of suicide candidate gene lists for bipolar disorder and schizophrenia with only two differentially expressed (suicide vs. non-suicide) genes in both bipolar and schizophrenia cohorts namely PLSCR4 and EMX2 by microarray analysis. The overlap of the two gene lists is small, suggesting few common, shared genes. For schizophrenia, but not bipolar disorder, the differential expression of PLSCR4 and EMX2 was confirmed by RT-PCR. The Phospholipid scramblase (PLSCR) is an integral membrane protein that catalyzes Ca2+-induced bidirectional movement of phospholipids [16]. Four isoforms have been cloned, and PLSCR4 was the major isoform expressed in the brain [17]. However, the biological role of the PLSCR4 remains unknown. While speculative, the changes in phospholipids membrane composition might have pleiotropic effects as evidence suggests that membrane composition can change G protein coupled receptors' functioning and downstream effector signalling [18,19] as well as voltage-dependent K+ channels [20]. EMX2 is a homeodomain containing transcription factor, which plays a crucial role in forebrain patterning and development in mouse models [21]. This finding suggests a possible neurodevelopmental process whereby variation in forebrain development may be a risk factor for suicide completion behaviour associated with schizophrenia. Of note, these differentially expressed genes are neural correlates of suicide and not necessarily causal. They could be epiphenomenon. The questions remain of 1) how these genes function to influence suicide and 2) what intermediate phenotype would be appropriate to demonstrate their possibly causal role.

Microarray studies with small sample sizes result in low statistical power and its attendant "noise discovery". This analysis and post-mortem sample size with replicates is noteworthy for an adequately powered sample to detect 1.3 fold changes, improving sensitivity, reliability, and signal to noise issues.

Previous smaller microarray studies have suggested that GABAA receptor subunits and glutamate-related genes were differentially expressed in bipolar disorder and schizophrenia as well as in suicide completers associated with these disorders [8,10]. We identified the up-regulation of gamma-amino butyric acid A receptor, α5 (GABRA5) subunit gene in suicide associated with bipolar disorder, confirming a previous report [10]. The expression levels of two glutamate-related genes, Glutamate-ammonia ligase (GLUL) and glial high affinity glutamate transporter member 3 (SLC1A3) were decreased in suicide completers with schizophrenia. The serotonergic and noradrenergic systems have been suggested to be associated with suicide [22]. However, no genes related to these two neurotransmitter systems were identified, consistent with a previous report [23]. This negative finding may suggest indirect effects on these neurotransmitter systems.

Genetic linkage studies have identified several loci associated with suicidal behaviors in bipolar disorders. Significant and suggestive linkages for suicide were mapped on chromosome 2, 5, and 10 from 162 bipolar pedigrees [24]. Among this study's suicide candidate genes associated with bipolar disorder, the tripartite motif-containing 23 (TRIM23) gene is located close to the significantly linked D5S1725 marker on chromosome 5. Another large scale genetic linkage study for bipolar disorder using 1060 individuals identified linkage on chromosome 10q25.3 for suicide attempts [25]. The microarray differentially expressed candidate gene EMX2 is included in this region of interest. Therefore, these two genes may be apt for future genetic association studies for suicide associated with bipolar disorder, proving causation.

While this reanalysis study has the strengths of a larger sample size, independent replicates, and well characterized patient samples from specific areas of cortex, the finding should be interpreted cautiously as this study has some limitations. First, the mixed cellular nature of the brain samples might lower sensitivity due to dilutional effects as opposed to pure neuronal cells of a specific cortical layer [12]. In general, most microarray studies with post-mortem brain tissues find fold changes of less than 2 fold, including this study. Second, although smoking, alcohol, and drug abuse were measured as confounding factors, all possible unmeasured, confounding variables for suicide cannot be formally excluded, such as severity of illness (especially the last episode for suicide completers), personality traits, hopelessness, agitation, depressive symptoms, and stress. Third, these findings are correlational and not causal. Fourth, these gene lists should be considered as provisional until confirmed by replications in independent sets of biological samples.

Conclusion

By reanalyzing a large microarray dataset, a list of differentially expressed candidate genes for suicide within bipolar disorder or within schizophrenia have been identified. The overlap of genes in common among these two gene lists is small, with a larger number of disorder specific genes being found. This finding suggests that disorder-specific pathways predominate over common pathways at the molecular level. Two novel candidate genes, PLSCR4 and EMX2, were confirmed as differentially expressed in schizophrenia between suicide completers vs. non-suicide groups.

Methods

Microarray data and Patient Samples

The brain tissues were meticulously collected in a standardized manner via pathologists in the offices of the Medical Examiner in several states with the families' permission under the aegis of the Stanley Foundation Brain Collection (Array Collection plus Consortium Collection) [26]. The selection of specimens, clinical information, diagnoses of patients, and processing of tissues were conducted by Stanley Foundation Consortium as described previously [26]. Gene expression profiling utilized post-mortem prefrontal cortex (Brodmann's Area 46/10) mRNA and Affymetrix Human Genome U133 Set A (HGU133A) using standardized techniques as described [26,27]. The prefrontal cortex was selected as the region of interest due to its role in executive functioning, impulsivity ("lack of premeditation"), and decision making. Disadvantageous decision making and impulsivity have been found to increase the risk of suicide [28,29].

The Stanley Foundation's microarray database is an anonymous, de-identified dataset without any protected health information. Patients' demographic variables used in this study are listed in Table 1.

The robust multi-array averages (RMA)-normalized microarray data from four independent studies were downloaded from the SMRIDB. Microarray data from the same platform, Affymetrix Human Genome U133 Set A (HGU133A), were used to avoid platform-to-platform variation. The platform contains 22,215 probe sets. Quality control analyses for each chip were described previously [15]

For the bipolar disorder cohort, the total dataset consisted of 49 suicide completers' gene chips and 58 non-suicide gene chips, while for schizophrenia cohort, the total dataset consisted of 22 suicide completer gene chips and 89 non-suicide chips. Among 45 bipolar samples, there were 22 suicide cases, and 23 non-suicide cases. Among 45 schizophrenia patients, there were 10 suicide cases, and 35 non-suicide cases. Two to three microarray chip datasets were generated from the each patient's sample. These repeated microarray data from each patient were treated as technical replicates.

Statistical analysis of microarray data

Microarray data was analyzed by a statistical method described previously with slight modifications [26]. Briefly the following steps were followed within each diagnostic group (see Table 1). First, the differentially expressed genes between suicide completers vs. non-suicide groups were filtered by average fold change (FC ≥ |1.3|) using the BRB-array tool[30] without covariates. Second, the influence of continuous demographic variables (such as age, post-mortem interval (PMI) and brain pH) with the nominal variable suicide was tested using ANOVA. Then, categorical variables such as sex, smoking, alcohol and drug abuse were tested using chi square tests of association (Statview software SAS, Cary, NC). In addition, correlation analyses of the demographic factors with expression levels of the differentially expressed probe sets from step 1 were performed. Continuous variables were analyzed by Spearman's rank correlation and categorical variables were tested by ANOVA. P-values were adjusted by False Discovery Rate (FDR) in both tests [31]. Third, significant confounding factors were tested as possible covariates for ANCOVA model inclusion with the following criteria: The variable was required to show both 1) significant association with suicide as well as 2) significant correlation with expression levels of the differentially expressed genes. However, no variables met the criteria in both disorder groups. Therefore, no covariates were used in the omnibus ANOVA, using the factor as suicide vs. non-suicide. As an exploratory analysis, a more liberal FDR P-value (< 0.1) of significance was selected for expression level differences as previously described [32,33], using the BRB array software tool's FDR default setting. For negative controls, we performed statistical analysis with the HGU133A normalization control probe sets using the same ANCOVA models, ensuring the adjustment did not produce "noise" or aberrant false positives. The Microarray Suite, version 5.0 (MAS5) software was used to filter genes with low expression levels as either present or absent, applying the detection call statistical algorithm. This algorithm suggests whether a gene is present or absent.

A power analysis estimated the sample sizes for detection of a 1.3 fold change in a gene with a significance criterion of P-value = 0.001 and a power of 0.90 using a previously described method [13]. This analysis estimated a minimum sample size of 27 cases per group for comparing suicide completers vs. non-suicide groups within bipolar disorder and 21 samples per group for comparing suicide completers vs. non-suicide completers within schizophrenia.

Real-time quantitative PCR

Total RNA from the dorsolateral prefrontal cortex (Brodmann area 46) of the Array Collection was used for this experiment. Complementary DNA was synthesized from DNA-free RNA with a random hexamer primer and Superscript III First-Strand Synthesis System according to the manufacturer's protocol (Invitrogen). Using a 384-well format with the Prism7900HT real-time detector (ABI), 2 μl aliquots of (10×) QuantiTect Primer Assay (validated primers to the specific gene of interest; Qiagen), 10 μl (2×) QuantiTect SYBR PCR Master mix (Qiagen), and 8 μl cDNA were mixed together for 20 μl total reaction volume and pipetted into single wells of the 384 PCR plate. Amplification conditions were: (1) 1 cycle for 2 min at 50°C, (2) 1 cycle for 15 min at 95°C, and (3) 45 cycles for 15 s at 94°C, 30 s at 60°C and 30 s at 72°C and fluorescence was measured during the 72°C step for each cycle as recommended by the manufacturer. The β-2 microglobulin (B2M) was chosen as an endogenous control for the normalization of target genes as it was consistently expressed in microarray samples. A total of 32 samples per each disorder were used for this experiment and run in duplicate. In the bipolar disorder cohort, there were 14 suicides and 18 non-suicide cases. In schizophrenia cohort, there were 5 suicides and 27 non-suicide cases. These samples were matched by age, race, gender, PMI, brain pH, side of the brain and quality of RNA. Reactions were quantified by the comparative Ct method using SDS2.2 software (ABI). This RT-PCR data was also statistically analyzed by the amplification plot method using the Data Analysis for Real Time PCR (DART-PCR) approach [34]. This method identified outliers in amplification efficiency by ANOVA and calculated mean expression levels. Statistical differences in expression levels between groups, namely suicide completers vs. non-suicides within a diagnostic group, were tested by one-tail, t-test with unequal variances (Microsoft Excel) as described [35,36]. To estimate average fold changes between groups, the mean expression values from the DART-PCR approach were used. The alternative 2(-Delta Delta Ct) method [37] for estimating fold change (using all the data without exclusions) verified this fold change estimate. As both methods gave similar estimates, only the DART-PCR approach estimates from mean expression levels were reported.

Functional annotation

The differentially expressed genes were functionally annotated using the DAVID integrated database query tool[38] and by the over-representational analysis method [39]. Functional annotations were based on biological process of Gene Ontology (GO) Consortium[40] at level 4. P-values less than 0.05 were considered significant.

Pathway Analysis

Biologically relevant networks were drawn from the lists of genes that were differentially expressed in bipolar disorder and schizophrenia. This data was generated through the use of Ingenuity Pathways Analysis (IPA) [41], a web-delivered application that enables the visualization and analysis of biologically relevant networks to discover, visualize, and explore relevant networks. Expression data sets containing gene identifiers (Affymetrix probe set ID) and their corresponding expression values as fold changes were uploaded as a tab-delimited text file. Each gene identifier was mapped to its corresponding gene object in the Ingenuity Pathways Knowledge Base. These genes, called Focus Genes, were then used as the starting point for generating biological networks. To start building networks, the application program queries the Ingenuity Pathways Knowledge Base for interactions between Focus Genes and all other gene objects stored in the knowledge base, and generates a set of networks. The program then computes a score for each network according to the fit of the network to the set of focus genes. The score indicates the likelihood of the Focus Genes in a given network being found together due to random chance. A score of greater than 2 indicates that there is a less than 1 in 100 chance that the Focus Genes were assembled randomly into a network due to random chance. The scores of the networks generated from the lists of differentially expressed genes were 24 for Bipolar disorder (Fig 3A) and 40 for Schizophrenia (Fig 3B).

Authors' contributions

SK performed the analysis of the microarray data. HG supervised the study. SK and HG wrote the manuscript. KC performed qRT-PCR validation experiment and SK performed analysis of the data. AFB provided important discussions. All authors read and approved of the final version of manuscript.

Acknowledgements

This work is generously supported by NIH (RO1 MH067211) and NARSAD (HG). The post-mortem brain collection was developed over many years through the generous permission of the families and the long-term support of the Stanley Foundation Brain Collection by the Stanley Medical Research Institute. We thank the Stanley Foundation Brain Collection and its Neuropathology Consortium for making its microarray database publicly available and their foresight in generating this dataset through the collaborative spirit and diligence of many investigators. In particular, we acknowledge The Stanley Medical Research Institute's Investigator's, Drs. E. Fuller Torrey, Robert H. Yolken, Michael B. Knable, Michael Elashoff, and Maree J. Webster along with their many collaborators, who made this work possible. We acknowledge Drs. C. Claasen and A. Kumar for helpful advice.

References

  1. McKeown RE, Cuffe SP, Schulz RM: US suicide rates by age group, 1970-2002: an examination of recent trends.

    Am J Public Health 2006, 96(10):1744-1751. PubMed Abstract | Publisher Full Text OpenURL

  2. Chen YW, Dilsaver SC: Lifetime rates of suicide attempts among subjects with bipolar and unipolar disorders relative to subjects with other Axis I disorders.

    Biol Psychiatry 1996, 39(10):896-899. PubMed Abstract | Publisher Full Text OpenURL

  3. Oquendo MA, Currier D, Mann JJ: Prospective studies of suicidal behavior in major depressive and bipolar disorders: what is the evidence for predictive risk factors?

    Acta Psychiatr Scand 2006, 114(3):151-158. PubMed Abstract | Publisher Full Text OpenURL

  4. Hawton K, Sutton L, Haw C, Sinclair J, Deeks JJ: Schizophrenia and suicide: systematic review of risk factors.

    Br J Psychiatry 2005, 187:9-20. PubMed Abstract | Publisher Full Text OpenURL

  5. Hawton K, Sutton L, Haw C, Sinclair J, Harriss L: Suicide and attempted suicide in bipolar disorder: a systematic review of risk factors.

    J Clin Psychiatry 2005, 66(6):693-704. PubMed Abstract | Publisher Full Text OpenURL

  6. Maris RW, Berman AL, Maltsberger JT: Summary and conclusions: What have we learned about suicide assessment and prediction? In Assessment and prediction of suicide. Edited by Maris RW, Berman AL, Maltsberger JT, Yufit RI. New York , Guilford Press; 1992:640–672. OpenURL

  7. Bunney WE, Bunney BG, Vawter MP, Tomita H, Li J, Evans SJ, Choudary PV, Myers RM, Jones EG, Watson SJ, Akil H: Microarray technology: a review of new strategies to discover candidate vulnerability genes in psychiatric disorders.

    Am J Psychiatry 2003, 160(4):657-666. PubMed Abstract | Publisher Full Text OpenURL

  8. Mirnics K, Levitt P, Lewis DA: Critical appraisal of DNA microarrays in psychiatric genomics.

    Biol Psychiatry 2006, 60(2):163-176. PubMed Abstract | Publisher Full Text OpenURL

  9. Gwadry FG, Sequeira A, Hoke G, Ffrench-Mullen JM, Turecki G: Molecular characterization of suicide by microarray analysis.

    Am J Med Genet C Semin Med Genet 2005, 133(1):48-56. PubMed Abstract | Publisher Full Text OpenURL

  10. Choudary PV, Molnar M, Evans SJ, Tomita H, Li JZ, Vawter MP, Myers RM, Bunney WE Jr., Akil H, Watson SJ, Jones EG: Altered cortical glutamatergic and GABAergic signal transmission with glial involvement in depression.

    Proc Natl Acad Sci U S A 2005, 102(43):15653-15658. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  11. Sequeira A, Gwadry FG, Ffrench-Mullen JM, Canetti L, Gingras Y, Casero RA Jr., Rouleau G, Benkelfat C, Turecki G: Implication of SSAT by gene expression and genetic variation in suicide and major depression.

    Arch Gen Psychiatry 2006, 63(1):35-48. PubMed Abstract | Publisher Full Text OpenURL

  12. Rossner MJ, Hirrlinger J, Wichert SP, Boehm C, Newrzella D, Hiemisch H, Eisenhardt G, Stuenkel C, von Ahsen O, Nave KA: Global transcriptome analysis of genetically identified neurons in the adult cortex.

    J Neurosci 2006, 26(39):9956-9966. PubMed Abstract | Publisher Full Text OpenURL

  13. Wei C, Li J, Bumgarner RE: Sample size for detecting differentially expressed genes in microarray experiments.

    BMC Genomics 2004, 5(1):87. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  14. SMRIDB[https://www.stanleygenomics.org] webcite

  15. Higgs BW, Elashoff M, Richman S, Barci B: An online database for brain disease research.

    BMC Genomics 2006, 7:70. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  16. Ikeda M, Kihara A, Igarashi Y: Lipid asymmetry of the eukaryotic plasma membrane: functions and related enzymes.

    Biol Pharm Bull 2006, 29(8):1542-1546. PubMed Abstract | Publisher Full Text OpenURL

  17. Wiedmer T, Zhou Q, Kwoh DY, Sims PJ: Identification of three new members of the phospholipid scramblase gene family.

    Biochim Biophys Acta 2000, 1467(1):244-253. PubMed Abstract | Publisher Full Text OpenURL

  18. Alves ID, Salgado GF, Salamon Z, Brown MF, Tollin G, Hruby VJ: Phosphatidylethanolamine enhances rhodopsin photoactivation and transducin binding in a solid supported lipid bilayer as determined using plasmon-waveguide resonance spectroscopy.

    Biophys J 2005, 88(1):198-210. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  19. Niu SL, Mitchell DC, Litman BJ: Optimization of receptor-G protein coupling by bilayer lipid composition II: formation of metarhodopsin II-transducin complex.

    J Biol Chem 2001, 276(46):42807-42811. PubMed Abstract | Publisher Full Text OpenURL

  20. Schmidt D, Jiang QX, MacKinnon R: Phospholipids and the origin of cationic gating charges in voltage sensors.

    Nature 2006, 444(7120):775-779. PubMed Abstract | Publisher Full Text OpenURL

  21. Hamasaki T, Leingartner A, Ringstedt T, O'Leary DD: EMX2 regulates sizes and positioning of the primary sensory and motor areas in neocortex by direct specification of cortical progenitors.

    Neuron 2004, 43(3):359-372. PubMed Abstract | Publisher Full Text OpenURL

  22. Mann JJ: Neurobiology of suicidal behaviour.

    Nat Rev Neurosci 2003, 4(10):819-828. PubMed Abstract | Publisher Full Text OpenURL

  23. Sibille E, Arango V, Galfalvy HC, Pavlidis P, Erraji-Benchekroun L, Ellis SP, John Mann J: Gene expression profiling of depression and suicide in human prefrontal cortex.

    Neuropsychopharmacology 2004, 29(2):351-361. PubMed Abstract | Publisher Full Text OpenURL

  24. Willour VL, Zandi PP, Badner JA, Steele J, Miao K, Lopez V, Mackinnon DF, Mondimore FM, Schweizer B, McInnis MG, Miller EB, Depaulo JR Jr., Gershon ES, McMahon FJ, Potash JB: Attempted Suicide in Bipolar Disorder Pedigrees: Evidence for Linkage to 2p12.

    Biol Psychiatry 2007, 61(5):725-727. PubMed Abstract | Publisher Full Text OpenURL

  25. Cheng R, Juo SH, Loth JE, Nee J, Iossifov I, Blumenthal R, Sharpe L, Kanyas K, Lerer B, Lilliston B, Smith M, Trautman K, Gilliam TC, Endicott J, Baron M: Genome-wide linkage scan in a large bipolar disorder sample from the National Institute of Mental Health genetics initiative suggests putative loci for bipolar disorder, psychosis, suicide, and panic disorder.

    Mol Psychiatry 2006, 11(3):252-260. PubMed Abstract | Publisher Full Text OpenURL

  26. Prabakaran S, Swatton JE, Ryan MM, Huffaker SJ, Huang JT, Griffin JL, Wayland M, Freeman T, Dudbridge F, Lilley KS, Karp NA, Hester S, Tkachev D, Mimmack ML, Yolken RH, Webster MJ, Torrey EF, Bahn S: Mitochondrial dysfunction in schizophrenia: evidence for compromised brain metabolism and oxidative stress.

    Mol Psychiatry 2004, 9(7):684-97, 643. PubMed Abstract | Publisher Full Text OpenURL

  27. Iwamoto K, Kakiuchi C, Bundo M, Ikeda K, Kato T: Molecular characterization of bipolar disorder by comparing gene expression profiles of postmortem brains of major mental disorders.

    Mol Psychiatry 2004, 9(4):406-416. PubMed Abstract | Publisher Full Text OpenURL

  28. Raust A, Slama F, Mathieu F, Roy I, Chenu A, Koncke D, Fouques D, Jollant F, Jouvent E, Courtet P, Leboyer M, Bellivier F: Prefrontal cortex dysfunction in patients with suicidal behavior.

    Psychol Med 2007, 37(3):411-419. PubMed Abstract | Publisher Full Text OpenURL

  29. Zermatten A, Van der Linden M, d'Acremont M, Jermann F, Bechara A: Impulsivity and decision making.

    J Nerv Ment Dis 2005, 193(10):647-650. PubMed Abstract | Publisher Full Text OpenURL

  30. BRB array tool[http://linus.nci.nih.gov/BRB-ArrayTools.html] webcite

  31. Benjamini Y, Hochberg Y: Controlling the False Discovery Rate: a practical and powerful approach to multiple testing.

    J Royal Stat Soc Ser B 1995, 57(1):289-300. OpenURL

  32. Hamatani T, Daikoku T, Wang H, Matsumoto H, Carter MG, Ko MS, Dey SK: Global gene expression analysis identifies molecular pathways distinguishing blastocyst dormancy and activation.

    Proc Natl Acad Sci U S A 2004, 101(28):10326-10331. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  33. Raponi M, Belly RT, Karp JE, Lancet JE, Atkins D, Wang Y: Microarray analysis reveals genetic pathways modulated by tipifarnib in acute myeloid leukemia.

    BMC Cancer 2004, 4:56. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  34. Peirson SN, Butler JN, Foster RG: Experimental validation of novel and conventional approaches to quantitative real-time PCR data analysis.

    Nucleic Acids Res 2003, 31(14):e73. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  35. Yuan JS, Reed A, Chen F, Stewart CN Jr.: Statistical analysis of real-time PCR data.

    BMC Bioinformatics 2006, 7:85. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  36. Ruxton GD: The unequal variance t-test is an underused alternative to Student's t-test and the Mann-Whitney U test.

    Behav Ecol 2006, 17(4):688-690. Publisher Full Text OpenURL

  37. Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method.

    Methods 2001, 25(4):402-408. PubMed Abstract | Publisher Full Text OpenURL

  38. DAVID integrated database query tool[http://david.abcc.ncifcrf.gov/] webcite

  39. Dennis G Jr., Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA: DAVID: Database for Annotation, Visualization, and Integrated Discovery.

    Genome Biol 2003, 4(5):P3. PubMed Abstract | BioMed Central Full Text OpenURL

  40. Gene Ontology Consortium[http://www.geneontology.org] webcite

  41. Ingenuity Pathways Analysis[http://www.Ingenuity.com] webcite