Meta-analysis of sex differences in gene expression in schizophrenia
© Qin et al. 2016
Published: 11 January 2016
Schizophrenia is a severe psychiatric disorder which influences around 1 % of the worldwide population. Differences between male and female patients with schizophrenia have been noted. There is an earlier age of onset in males compared with females with this diagnosis, and in addition, there are differences in symptom profiles between the sexes. The underlying molecular mechanism of sex difference remains unclear. Here we present a comprehensive analysis to reveal the sex differences in gene expression in schizophrenia with stringent statistics criteria. We compiled a data set consisting of 89 male controls, 90 male schizophrenia patients, 35 female controls and 32 female schizophrenia patients from six independent studies of the prefrontal cortex (PFC) in postmortem brain. When we tested for a sex by diagnosis interaction on gene expression, 23 genes were up-regulated and 23 genes were down-regulated in the male group (q-value < 0.05), several genes are related to energy metabolism, while 4 genes are located on sex chromosome. No genes were statistically significant in the female group when multiple testing correction were conducted (q-value <0.05), most likely due to the small sample size. Our protocol and results from the male group provide a starting point for identifying the underlying different mechanism between male and female schizophrenia patients.
Schizophrenia is a severe psychiatric disorder with a population frequency of approximately 1 % . Schizophrenia is a syndrome characterized by positive symptoms such as delusions, hallucinations, disorganized speech and grossly disorganized or catatonic behavior; negative symptoms such as affective flatterning, alogia, or avolition , . The etiology and pathophysiological mechanisms of the disorder are not well understood. Research to date indicates that schizophrenia is a multi-factorial neurodevelopmental impairment of the brain that could be attributed to both genetic and environmental factors , , , , . Gene expression is readout of both the genetic and the environmental factors that contribute to the pathophysiology of schizophrenia. Analysis of human postmortem brain is a powerful approach for the identification of risk factors for schizophrenia, because unlike studies of living patients, detailed molecular investigations can be performed directly in the critical brain regions of interest.
The pathophysiology of schizophrenia is likely to be different between males and females. Sex differences have been noted in several epidemiological analyses. For example, several studies indicate that men have a slightly higher incidence of schizophrenia compared with women. In addition, males have an earlier age of onset of schizophrenia, between 18–25 years of age, compared with the female age of onset which is 25–35 years . The symptoms exhibited by male and female patients with schizophrenia also differ. Males tend to have a greater vulnerability to negative symptoms and traits of disorganization, while females more frequently exhibit depressive symptoms . These findings suggest that different underlying mechanisms of schizophrenia occur in males and females. Therefore, we have investigated sex differences in schizophrenia to gain a better understanding of the pathophysiological mechanisms underpinning this disorder.
To identify the biological factors involved in the pathogenesis of schizophrenia and how they are differentially influenced in the sexes, we have investigated microarray expression data from the prefrontal cortex (PFC) in postmortem brain. The PFC region has been strongly associated with deficits of executive function and other cognitive symptoms that occur in patients with schizophrenia. Gene expression within the PFC has been studied extensively using the microarray approach -, , , , , , , , . However, statistical analyses of gene expression data from individual small cohorts have lacked sufficient statistical power to avoid conflicting data in these different studies , . Meta-analysis is a strategy by which these problems could be addressed, because the data from multiple studies can be combined, thus increasing the statistical power available. In this study, we have collected gene expression data generated from post-mortem cohorts of schizophrenia cases and psychiatrically healthy comparison groups that are included in publicly available databases. We have tested for sex differences in PFC gene expression in schizophrenia, using the meta-analysis paradigm.
Methods and materials
Public microarray datasets of postmortem gene expression in schizophrenia
Public postmortem microarray datasets used in this study
Stanley Bahn 
Stanley AltarC 
Mclean Harvard Brain Bank 
Prefrontal Cortex (BA9)
Prefrontal Cortex (BA46)
Maycox GSE17612 
Anterior prefrontal cortex (BA10)
Narayan GSE21138 
Differential expression analysis of each probe set
Expression values of each probe set were modeled using a fixed effect linear model approach, where disease status and imbalanced covariates between two groups are treated as fixed effects to be estimated from data. A model selection procedure was also employed for each probe set to address the confounding effect of imbalanced covariates. The details of the procedure will be described in next section. For each probe set, the t-statistic for the disease effect was then extracted from the selected model. P-values were computed using two-sided t-test. The resulting P-values were further converted to q-values using the qvalue package in R  which were defined as the minimum positive False Discovery Rate (q-value(t): Pr(H = 0|T > =t) in Bayesian interpretation over a nested rejection region containing observed statistics t) for multiple testing correction.
The observed covariate imbalance between two groups was analyzed to avoid generating misleading result. However modelling covariates for every gene could unnecessarily diminish statistical power if the covariate does not influence the expression of the gene. In our study, we consider only modelling imbalanced covariates for each probe set. We first obtained a probe set (gene) list where the covariate influence on the gene expression is determined with confidence. This refined gene list was generated using a previous method in a postmortem brain gene expression study, in which a correlation analysis was used to evaluate covariates such as age, post-mortem interval, brain PH etc. that influenced the expression of specific genes across multiple post-mortem normal brain datasets . We extracted genes with meta-Q ≤ 0.01 to indicate that the gene was significantly influenced by a particular covariate, and separate it into two lists: positively-correlated and negatively-correlated. Notice that one probe set could be mapped to multiple genes, we excluded those probe sets appearing in the two lists. For each predefined covariate influenced probe set, we first modelled the expression value with a linear model including that covariate. If the direction of the fitted covariate estimates was inconsistent with the pre-refined list, we exclude this covariate and re-estimate a reduced model for this probe set. For the rest of the probe sets, no covariate adjustment was performed.
Selection of sex-specific differentially expressed genes
To identify genes that had a sex by diagnosis interaction with schizophrenia, we used a strategy similar to . Each individual in our dataset is assigned to one of four subgroups: Schizophrenia Male, Control Male, Schizophrenia Female and Control Female. Individuals could also be combined into two groups based on their diagnosis: Schizophrenia and Control group, or based on their sex: Male and Female group. Differential expression analysis was first performed using the procedure described in section 2.2 within each sex. After the initial probe set list was obtained from each sex, we further eliminated those probe sets that are associated with schizophrenia regardless of sex when all of the following criteria were met: (a) the difference between Schizophrenia and Control groups was statistically significant after multiple test correction (q-value < 0.05); (b) the fold change of Schizophrenia Female vs. Control Females, and Schizophrenia Male vs Control Males should be in the same direction i.e. both higher or both lower; (c) the expression difference was not significant between Male and Female groups (defined as p > 0.05 between Female and Male group). After removal of these probe sets from the initial probe set list, we sorted the remaining probe sets within each sex by q-value and report the top ranked probe sets (q-value < 0.05) as sex-specific differentially expressed genes.
Function enrichment of differentially expressed genes
All differentially expressed genes, along with their Affymetrix ID numbers were imported into EASE (Expression Analysis Systematic Explorer) in DAVID (Database for Annotation, Visualization and Integrated Discovery), and were used to identify functionally significant gene classes (https://david.ncifcrf.gov/) . This webserver uses statistical methods to map and identify functional gene categories (for example, Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) or BioCarta), which are enriched in the significant gene list compared with their presence on the array.
Batch and covariate adjustment
A total of 22277 probe sets is analyzed in this study. Before further analysis, by using hierarchical clustering analysis and Principal Component Analysis (PCA), we found it was necessary to correct for “batch effect” as samples from the same study were clustered together. ComBat  was used to correct this technical bias by treating each study as a batch. After ComBat adjustment, the hierarchical clustering and PCA results show that no significant clustering remained in the dataset and the adjusted expression data was suitable for further analysis (Additional file 1).
Demographic data of postmortem subjects
50.29 ± 17.2
51.09 ± 18.90
47.85 ± 16.44
48.45 ± 17.71
56.49 ± 17.89
58.53 ± 20.38
89 M:35 F
90 M:32 F
6.48 ± 0.29
6.36 ± 0.29
6.50 ± 0.29
6.35 ± 0.26
6.42 ± 0.29
6.40 ± 0.38
25.49 ± 14.8
25.09 ± 15.9
26.12 ± 14.37
25.79 ± 16.45
23.91 ± 15.98
23.14 ± 14.34
Covariate adjustment summary
Two group comparison
Schizophrenia vs control
Schizophrenia-male vs control-male
Male vs Female
Imbalanced covariate between two groups
Proportion of probe sets adjusteda (%)
Genes with altered PFC expression in schizophrenia
We also examined whether the identified 778 probe sets are associated with schizophrenia in case control studies of genetic polymorphisms. We compared our gene list with those deposited in SZGene (www.SZgene.org) database which contains the most comprehensive review of schizophrenia association studies . 80 probe sets representing 68 unique genes are found to be genetically associated with schizophrenia. Genes previously showing strong genetic evidence implicated in schizophrenia identified in this study include: regulator of G-protein signaling 4 (RGS4) ; discoidin domain receptor family, member 1(DDR1) ; and the selenium binding protein 1(SELENBP1) . The full result is included in Additional file 3.
Sex differences in PFC gene expression in schizophrenia
Genes with altered expression in males with schizophrenia: downregulated probe sets
ATP-binding cassette, sub-family G (WHITE), member 2
GABA(A) receptor-associated protein like 1
glutamic-oxaloacetic transaminase 1, soluble (aspartate aminotransferase 1)
kinesin light chain 1
acyl-CoA thioesterase 7
mitochondrial ribosomal protein L23
transforming, acidic coiled-coil containing protein 2
ATP synthase, H+ transporting, mitochondrial F1 complex, alpha subunit 1, cardiac muscle
ring finger protein, transmembrane 2
ATP synthase, H+ transporting, mitochondrial F1 complex, beta polypeptide
tumor suppressor candidate 2
brain expressed, X-linked 1
NHP2 non-histone chromosome protein 2-like 1 (S. cerevisiae)
pleckstrin homology domain containing, family B (evectins) member 2
lin-7 homolog B (C. elegans)
fibroblast growth factor 22
unconventional SNARE in the ER 1 homolog (S. cerevisiae)
glutathione S-transferase alpha 4
transmembrane protein 8B
AFG3 ATPase family gene 3-like 2 (S. cerevisiae)
Genes with altered expression in males with schizophrenia: upregulated probe sets
phosphatidic acid phosphatase type 2B
Rho-related BTB domain containing 3
SRY (sex determining region Y)-box 9
DNA-damage-inducible transcript 4
palladin, cytoskeletal associated protein
palladin, cytoskeletal associated protein
fermitin family member 2
bobby sox homolog (Drosophila)
geminin, DNA replication inhibitor
fermitin family member 2
H3 histone, family 3B (H3.3B)
stimulator of chondrogenesis 1
CDC14 cell division cycle 14 homolog B (S. cerevisiae)
H3 histone, family 3B (H3.3B)
SAM and SH3 domain containing 1
acyl-CoA oxidase 1, palmitoyl
family with sequence similarity 38, member A
midline 1 (Opitz/BBB syndrome)
palladin, cytoskeletal associated protein
Yes-associated protein 1
FYN oncogene related to SRC, FGR, YES
PABPC1 /// RLIM
poly(A) binding protein, cytoplasmic 1 /// ring finger protein, LIM domain interacting
par-3 partitioning defective 3 homolog (C. elegans)
Function annotation of male-specific differentially expressed genes
The functions of 46 genes associated with schizophrenia in male patients were manually inspected. Of the genes with significantly lower expression in males with schizophrenia, several were related to energy metabolism (ATP5B, ATP5A1, MRPL23, AFG3L2, ABCG2). 4 genes (BEX1, UBL4A, CD99 and MID1) located on the sex chromosome are identified. We also detected differential expression of 4 genes in the males with schizophrenia that were previously identified by Mistry et al.  in which sex was not considered. These were Rho-Related BTB Domain-Containing Protein 3 (RHOBTB3), Bobby Sox homolog (BBX), H3 Histone, Family 3B (H3F3B) and pleckstrin homology domain containing, family B (evectins) member 2 (PLEKHB2). Finally, using DAVID webserver, we performed an enrichment analysis to systematically identify over-representation of biological processes or pathways that are altered in the PFC in male schizophrenia patients. After correction for multiple comparisons, we were unable to identify any significant biological process (False Discovery Rate <0.05) in the GO term Biological Process database. We achieved similarly negative results using the KEGG pathway database.
We have reported the gene expression differences that show a sex by diagnosis interaction in the PFC in schizophrenia. To our knowledge, this is the first study using a meta-analytical approach to identify sex differences in this brain region in patients with schizophrenia. There are limited data on sex differences in schizophrenia at a molecular level , , although evidence from epidemiological and animal studies indicates that sex differences exist in this disorder , , . Individual post-mortem gene expression studies have low statistical power to identify gene expression differences in schizophrenia. This is, most often due to the small sample sizes and moderate gene expression differences between the diagnostic groups. Meta-analysis, on the other hand, addresses this problem and increases statistical power by combining samples from different subject cohorts. The results obtained from our meta-analysis are robust at the statistical level . Our findings open a new window to understand the different pathophysiological mechanisms that lead to schizophrenia in males and females.
In the differential analyses between schizophrenia and control group, we identified the most number of genes most of which overlapped with the results of two previously published studies, while the other two published studies had little overlap between each other (Fig. 1). To find the reason for the difference, we run the same procedure as in Mistry’s for our dataset. We find that the main difference comes from different treatment of “batch effect”. In Mistry’s study, they treated each experiment date and study as batch variable. A total of 50 batch date and 6 study was modelled in their linear model framework. Introduction of too many predictor variables will decrease the degrees of freedom of t-test on the estimated disease effect coefficient, thus decreasing the likelihood of rejecting the null hypothesis, leading to higher Type-II error. Based on the observation of hierarchical clustering and PCA result, we consider that modelling each study alone would be sufficient to correct the batch effect. The second reason explaining the difference is covariate adjustment. We only include a covariate in the model when its influence on the gene is confirmed with confidence; while in Mistry’s study, they assumed that the covariate influenced every gene resulting in unnecessary inclusion of unrelated variables in the model. Therefore the approach taken in this study has more detection power and leads to discovery of more genes.
Our analyses identified 46 genes that were differentially expressed specifically in male patients with schizophrenia. This finding of 50 probe sets is much larger than the expected number of false positives according to our selection procedure. The expected number of false positives was calculated to be 6.9 (138*0.05 = 6.9). Five genes were related to energy metabolism (ATP5B, ATP5A1, MRPL23, AFG3L2, ABCG2). Genes from this function category are consistently implicated in studies of schizophrenia. Another gene that had altered expression in the male schizophrenia group encodes γ-aminobutyric acid receptor-associated protein-like 1 (GABARAPL1) is an early estrogen-induced gene that when overexpressed, interacts with GABA-A or κ-opioid receptors, and plays a role in cell proliferation and cellular metabolic processes . Function enrichment analysis generates negative results for these genes, as indicated in other microarray studies of the PFC in schizophrenia . These results suggest that a diverse number of molecular functions are disrupted in males with schizophrenia.
No genes could be identified in female group after multiple testing correction (q-value < 0.05). To determine if this is due to a much smaller sample size than in male group, we randomly picked the same number of control and schizophrenia subjects from the male groups ten times and run differential analysis on these samples. The procedure was repeated 100 times. No significant difference could be detected in the expression level of any gene in the phenotype groups (data not shown). We then gradually increase the number of samples until genes could be detected with q-value < 0.05 in each run. We started to identify differentially expressed genes when there are 60 controls and 60 schizophrenia in the sample pool. This analysis showed that increasing the number of samples would improve the likelihood of identifying schizophrenia-associated genes in the females.
Genes reported by Vawter et al.  in this study
In summary, this is the most comprehensive and up-to-date analysis of sex differences in prefrontal cortex gene expression in schizophrenia. Some of our data consolidate the reports of previously published papers. Our results indicate that further investigation of sex differences in schizophrenia is required . These data bring us closer to understanding the different molecular mechanisms underpinning schizophrenia in males and females, so that novel targets for antipsychotic drug development can be identified. Genes with altered expression in schizophrenia can also serve as biological markers for the disorder, so that biochemical diagnostic tools can facilitate the practice of clinical psychiatry.
This work was partially supported by the Hans W. Vahlteich Fund (MS), and UIC Department of Bioengineering (WQ, CL, HL). Publication fee are paid by Department of Bioengineering at UIC and by Shanghai Jiaotong University.
- Boyer-Guittaut M, Poillet L, Liang Q, Bôle-Richard E, Ouyang X, Benavides GA, et al. The Role of Gabarapl1/Gec1 in Autophagic Flux and Mitochondrial Quality Control in Mda-Mb-436 Breast Cancer Cells. Autophagy. 2014;10(6).Google Scholar
- Brown AS: The Environment and Susceptibility to Schizophrenia. Prog Neurobiol. 2011, 93 (1): 23-58. 10.1016/j.pneurobio.2010.09.003.View ArticlePubMedGoogle Scholar
- Cantuti-Castelvetri I, Keller-McGandy C, Bouzou B, Asteris G, Clark TW, Frosch MP, et al: Effects of Gender on Nigral Gene Expression and Parkinson Disease. Neurobiol Dis. 2007, 26 (3): 606-14. 10.1016/j.nbd.2007.02.009.View ArticlePubMedPubMed CentralGoogle Scholar
- Chu TT, Liu Y, Kemether E: Thalamic Transcriptome Screening in Three Psychiatric States. J Hum Genet. 2009, 54 (11): 665-75. 10.1038/jhg.2009.93.View ArticlePubMedGoogle Scholar
- Dean B, Keriakous D, Scarr E, Thomas EA, et al: Gene Expression Profiling in Brodmann’s Area 46 from Subjects with Schizophrenia. Aust N Z J Psychiatry. 2007, 41 (4): 308-20. 10.1080/00048670701213245.View ArticlePubMedGoogle Scholar
- Garbett KR, Gal-Chis RO, Gaszner GA, Lewis DA, Mirnics KÁ. Transcriptome Alterations in the Prefrontal Cortex of Subjects with Schizophrenia Who Committed Suicide. Neuropsychopharmacol Hung. 2008;10(1):9–14.PubMedGoogle Scholar
- Gejman PV, Sanders AR, Kendler KS, et al: Genetics of Schizophrenia: New Findings and Challenges. Annu Rev Genomics Hum Genet. 2011, 12: 121-144. 10.1146/annurev-genom-082410-101459.View ArticlePubMedGoogle Scholar
- Mirnics K, Middleton FA, Stanwood GD, Lewis DA, Levitt P. Disease-specific changes in regulator of G-protein signaling 4 (RGS4) expression in schizophrenia. Mol Psychiatry 2001; 6: 293–301.View ArticlePubMedGoogle Scholar
- Higgs BW, Elashoff M, Richman S, Barci B: An Online Database for Brain Disease Research. BMC Genomics. 2006, 7: 70-10.1186/1471-2164-7-70.View ArticlePubMedPubMed CentralGoogle Scholar
- Huang DW, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID Bioinformatics Resources. Nature Protoc. 2009, 4 (1): 44-57. 10.1038/nprot.2008.211.View ArticleGoogle Scholar
- Horvath S, Janka Z, Mirnics K: Analyzing Schizophrenia by DNA Microarrays. Biol Psychiatry. 2011, 69 (2): 157-62. 10.1016/j.biopsych.2010.07.017.View ArticlePubMedPubMed CentralGoogle Scholar
- Johnson WE, Li C, Rabinovic A: Adjusting Batch Effects in Microarray Expression Data Using Empirical Bayes Methods. Biostatistics. 2007, 8 (1): 118-27. 10.1093/biostatistics/kxj037.View ArticlePubMedGoogle Scholar
- Katsel P, Davis KL, Gorman JM, Haroutunian V: Variations in Differential Gene Expression Patterns across Multiple Brain Regions in Schizophrenia. Schizophr Res. 2005, 77 (2–3): 241-52. 10.1016/j.schres.2005.03.020.View ArticlePubMedGoogle Scholar
- Kumarasinghe N, Tooney PA, Schall U: Finding the Needle in the Haystack: A Review of Microarray Gene Expression Research into Schizophrenia. Aust N Z J Psychiatry. 2012, 46 (7): 598-610. 10.1177/0004867412442405.View ArticlePubMedGoogle Scholar
- Martins-de-Souza D, Schmitt A, Röder R, Lebar M, Schneider-Axmann T, Falkai P, et al: Sex-Specific Proteome Differences in the Anterior Cingulate Cortex of Schizophrenia. J Psychiatr Res. 2010, 44 (14): 989-91. 10.1016/j.jpsychires.2010.03.003.View ArticlePubMedGoogle Scholar
- Maycox PR, Kelly F, Taylor A, Bates S, Reid J, Logendra R, et al: Analysis of Gene Expression in Two Large Schizophrenia Cohorts Identifies Multiple Changes Associated with Nerve Terminal Function. Mol Psychiatry. 2009, 14 (12): 1083-94. 10.1038/mp.2009.18.View ArticlePubMedGoogle Scholar
- Mistry M, Gillis J, Pavlidis P: Genome-Wide Expression Profiling of Schizophrenia Using a Large Combined Cohort. Mol Psychiatry. 2013, 18 (2): 215-25. 10.1038/mp.2011.172.View ArticlePubMedGoogle Scholar
- Mistry M, Pavlidis P: A Cross-Laboratory Comparison of Expression Profiling Data from Normal Human Postmortem Brain. Neuroscience. 2010, 167 (2): 384-95. 10.1016/j.neuroscience.2010.01.016.View ArticlePubMedPubMed CentralGoogle Scholar
- Narayan S, Tang B, Head SR, Gilmartin TJ, Sutcliffe JG, Dean B, et al: Molecular Profiles of Schizophrenia in the Cns at Different Stages of Illness. Brain Res. 2008, 1239: 235-48. 10.1016/j.brainres.2008.08.023.View ArticlePubMedPubMed CentralGoogle Scholar
- Ochoa S, Usall J, Cobo J, Labad X, Kulkarni J: Gender Differences in Schizophrenia and First-Episode Psychosis: A Comprehensive Literature Review. Schizophr Res Treatment. 2012, 2012: 916198-PubMedPubMed CentralGoogle Scholar
- Owen MJ, Craddock N, O’Donovan MC: Suggestion of Roles for Both Common and Rare Risk Variants in Genome-Wide Studies of Schizophrenia. Arch Gen Psychiatry. 2010, 67 (7): 667-673. 10.1001/archgenpsychiatry.2010.69.View ArticlePubMedGoogle Scholar
- Pérez-Santiago J, Diez-Alarcia R, Callado LF, Zhang JX, Chana G, White CH, et al: A Combined Analysis of Microarray Gene Expression Studies of the Human Prefrontal Cortex Identifies Genes Implicated in Schizophrenia. J Psychiatr Res. 2012, 46 (11): 1464-74. 10.1016/j.jpsychires.2012.08.005.View ArticlePubMedGoogle Scholar
- Roth TL, Lubin FD, Sodhi M, Kleinman JE: Epigenetic Mechanisms in Schizophrenia. Biochimica et Biophysica Acta (BBA)-General Subjects. 2009, 1790 (9): 869-877. 10.1016/j.bbagen.2009.06.009.View ArticleGoogle Scholar
- Sequeira PA, Martin MV, Vawter MP: The First Decade and Beyond of Transcriptional Profiling in Schizophrenia. Neurobiol Dis. 2012, 45 (1): 23-36. 10.1016/j.nbd.2011.03.001.View ArticlePubMedGoogle Scholar
- Storey JD. The Positive False Discovery Rate: A Bayesian Interpretation and the Q-Value.The Annals of Statistics. 2003; 31(6): 2013-2035.Google Scholar
- Tandon R, Gaebel W, Barch DM, Bustillo J, Gur RE, Heckers S, et al: Definition and Description of Schizophrenia in the Dsm-5. Schizophr Res. 2013, 150 (1): 3-10. 10.1016/j.schres.2013.05.028.View ArticlePubMedGoogle Scholar
- Teaching Network, EURON and De Crespigny Park London: Jim Van Os, Shitij Kapur. Lancet. 2009, 374: 635-45. 10.1016/S0140-6736(09)60995-8.View ArticleGoogle Scholar
- Uher R. Gene–Environment Interactions in Severe Mental Illness. Frontiers in Psychiatry. 2014;5(48).Google Scholar
- Vawter MP, Weickert CS, et al: Gene Expression of Metabolic Enzymes and a Protease Inhibitor in the Prefrontal Cortex Are Decreased in Schizophrenia. Neurochem Res. 2004, 29 (6): 1245-1255. 10.1023/B:NERE.0000023611.99452.47.View ArticlePubMedGoogle Scholar
- Wu YC, Hill RA, Gogos A, Van Den Buuse M: Sex Differences and the Role of Estrogen in Animal Models of Schizophrenia: Interaction with Bdnf. Neuroscience. 2013, 239: 67-83. 10.1016/j.neuroscience.2012.10.024.View ArticlePubMedGoogle Scholar
- Allen NC, Bagade S, McQueen MB, Ioannidis JPA, Kavvoura FK, Khoury MJ, et al: Systematic Meta-Analyses and Field Synopsis of Genetic Association Studies in Schizophrenia: The SzGene Database. Nat Genet. 2014, 40 (7): 827-34. 10.1038/ng.171.View ArticleGoogle Scholar
- Prasad KMR, Chowdari KV, Nimgaonkar VL, Talkowski ME, Lewis DA, Keshavan MS: Genetic polymorphisms of the RGS4 and dorsolateral prefrontal cortex morphometry among first episode schizophrenia patients. Molec Psychiat. 2005, 10: 213-219. 10.1038/sj.mp.4001562.View ArticlePubMedGoogle Scholar
- Roig B, Virgos C, Franco N, Martorell L, Valero J, Costas J, et al: The discoidin domain receptor 1 as a novel susceptibility gene for schizophrenia. Mol Psychiatry. 2007, 12 (9): 833-41. 10.1038/sj.mp.4001995.View ArticlePubMedGoogle Scholar
- Kanazawa T, Glatt SJ, Faraone SV, Hwu HG, Yoneda H, Tsuang MT: Family-based association study of SELENBP1 in schizophrenia. Schizophr Res. 2009, 113 (2–3): 268-72. 10.1016/j.schres.2009.06.011.View ArticlePubMedGoogle Scholar
- Stanley Medical Research Institute. https://www.stanleygenomics.org/. Accessed 1 Feb 2013.
- Harvard Brain Tissure Resource Center (Mclean 66 collection). http://www.brainbank.mclean.org/. Accessed 1 Feb 2013.
- Garbett K, Gal-Chis R, Gaszner G, Lewis DA, Mirnics K: Transcriptome alterations in the prefrontal cortex of subjects withschizophrenia who committed suicide. Neuropsychopharmacol Hung. 2008, 10: 9-14.PubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.