Population-specific haplotype association of the postsynaptic density gene DLG4 with schizophrenia, in family-based association studies

Dysregulated Signaling at Postsynaptic Density: A Systematic Review and Translational Appraisal for the Pathophysiology, Clinics, and Antipsychotics' Treatment of Schizophrenia

Emerging evidence from genomics, post-mortem, and preclinical studies point to a potential dysregulation of molecular signaling at postsynaptic density (PSD) in schizophrenia pathophysiology. The PSD that identifies the archetypal asymmetric synapse is a structure of approximately 300 nm in diameter, localized behind the neuronal membrane in the glutamatergic synapse, and constituted by more than 1000 proteins, including receptors, adaptors, kinases, and scaffold proteins. Furthermore, using FASS (fluorescence-activated synaptosome sorting) techniques, glutamatergic synaptosomes were isolated at around 70 nm, where the receptors anchored to the PSD proteins can diffuse laterally along the PSD and were stabilized by scaffold proteins in nanodomains of 50-80 nm at a distance of 20-40 nm creating "nanocolumns" within the synaptic button. In this context, PSD was envisioned as a multimodal hub integrating multiple signaling-related intracellular functions. Dysfunctions of glutamate signaling have been postulated in schizophrenia, starting from the glutamate receptor's interaction with scaffolding proteins involved in the N-methyl-D-aspartate receptor (NMDAR). Despite the emerging role of PSD proteins in behavioral disorders, there is currently no systematic review that integrates preclinical and clinical findings addressing dysregulated PSD signaling and translational implications for antipsychotic treatment in the aberrant postsynaptic function context. Here we reviewed a critical appraisal of the role of dysregulated PSD proteins signaling in the pathophysiology of schizophrenia, discussing how antipsychotics may affect PSD structures and synaptic plasticity in brain regions relevant to psychosis.

Application of per-Residue Energy Decomposition to Design Peptide Inhibitors of PSD95 GK Domain

Specific interaction between the postsynaptic density protein 95 (PSD95) and synapse-associated protein 90/postsynaptic density 95–associated protein (SAPAP) is crucial for excitatory synaptic development and plasticity. Designing inhibitors that target the guanylate kinase (GK) domain of PSD95, which is responsible for the interaction, is a promising manipulation tool for the investigation of the function of PSD95 GK and the etiology of its related psychiatric disorders. Herein, we designed new peptide inhibitors of PSD95 GK/SAPAP with higher binding affinity by using molecular dynamics simulations. First, the interactions between PSD95 GK and their reported phosphorylated and unphosphorylated peptides were explored by molecular dynamics simulations. Besides the hydrogen bonding interactions mediated by the phospho-serine (p-Ser) or corresponding phosphomimic residue Asp/Glu, the hydrophobic interactions from the other amino acids also contribute to the PSD95 GK/SAPAP interaction. As an unphosphorylated synthetic peptide with moderate binding affinity and relatively lower molecular weight, the QSF inhibitory peptide was selected for further modification. Based on per-residue energy decomposition results of the PSD95 GK/QSF complex, ten peptides were designed to enhance the binding interactions, especially the hydrophobic interactions. The top-ranked five peptides with lower binding energy were eventually synthesized. The binding affinities of the synthesized peptides were determined using fluorescence polarization (FP) assay. As expected, all peptides have higher binding affinity than the QSF peptide (Ki = 5.64 ± 0.51 μM). Among them, F10W was the most potent inhibitor (Ki = 0.75 ± 0.25 μM), suggesting that enhancement of the hydrophobic interactions is an important strategy for the design of new inhibitory peptides targeting PSD95 GK.

Integrative Analysis Identified Key Schizophrenia Risk Factors from an Abnormal Behavior Mouse Gene Set

Schizophrenia (SCZ) is a severe chronic psychiatric illness with heterogeneous symptoms. However, the pathogenesis of SCZ is unclear, and the number of well-defined SCZ risk factors is limited. We hypothesized that an abnormal behavior (AB) gene set verified by mouse model experiments can be used to better understand SCZ risks. In this work, we carried out an integrative bioinformatics analysis to study two types of risk genes that are either differentially expressed (DEGs) in the case-control study data or carry reported SCZ genetic variants (MUTs). Next, we used RNA-Seq expression data from the hippocampus (HIPPO) and dorsolateral prefrontal cortex (DLPFC) to define the key genes affected by different types (DEGs and MUTs) in different brain regions (DLPFC and HIPPO): DLPFC-kDEG, DLPFC-kMUT, HIPPO-kDEG, and HIPPO-kMUT. The four hub genes (SHANK1, SHANK2, DLG4, and NLGN3) of the biological functionally enriched terms were strongly linked to SCZ via gene co-expression network analysis. Then, we observed that specific spatial expressions of DLPFC-kMUT and HIPPO-kMUT were convergent in the early stages and divergent in the later stages of development. In addition, all four types of key genes showed significantly larger average protein-protein interaction degrees than the background. Comparing the different cell types, the expression of four types of key genes showed specificity in different dimensions. Together, our results offer new insights into potential risk factors and help us understand the complexity and regional heterogeneity of SCZ.

Genomic analyses reveal distinct genetic architectures and selective pressures in buffaloes

BACKGROUND

The domestic buffalo (Bubalus bubalis) is an essential farm animal in tropical and subtropical regions, whose genomic diversity is yet to be fully discovered.

RESULTS

In this study, we describe the demographic events and selective pressures of buffalo by analyzing 121 whole genomes (98 newly reported) from 25 swamp and river buffalo breeds. Both uniparental and biparental markers were investigated to provide the final scenario. The ancestors of swamp and river buffalo diverged ∼0.23 million years ago and then experienced independent demographic histories. They were domesticated in different regions, the swamp buffalo at the border between southwest China and southeast Asia, while the river buffalo in south Asia. The domestic stocks migrated to other regions and further differentiated, as testified by (at least) 2 ancestral components identified in each subspecies. Different signals of selective pressures were also detected in these 2 types of buffalo. The swamp buffalo, historically used as a draft animal, shows selection signatures in genes associated with the nervous system, while in river dairy breeds, genes under selection are related to heat stress and immunity.

CONCLUSIONS

Our findings substantially expand the catalogue of genetic variants in buffalo and reveal new insights into the evolutionary history and distinct selective pressures in river and swamp buffalo.

Transgenerational Bisphenol A Causes Deficits in Social Recognition and Alters Postsynaptic Density Genes in Mice

Bisphenol A (BPA) is a ubiquitous endocrine-disrupting chemical. Developmental exposure produces changes in behavior and gene expression in the brain. Here, we examined social recognition behaviors in mice from the third familial generation (F3) after exposure to gestational BPA. Second-generation mice were bred in one of four mating combinations to reveal whether characteristics in F3 were acquired via maternal or paternal exposures. After repeated habituation to the same mouse, offspring of dams from the BPA lineage failed to display increased investigation of a novel mouse. Genes involved in excitatory postsynaptic densities (PSDs) were examined in F3 brains using quantitative PCR. Differential expression of genes important for function and stability of PSDs were assessed at three developmental ages. Several related PSD genes-SH3 and multiple ankyrin repeat domains 1 (Shank1), Homer scaffolding protein 1c (Homer1c), DLG associated protein 1 (Gkap), and discs large MAGUK scaffold protein 4 (PSD95)-were differentially expressed in control- vs BPA-lineage brains. Using a second strain of F3 inbred mice exposed to BPA, we noted the same differences in Shank1 and PSD95 expression in C57BL/6J mice. In sum, transgenerational BPA exposure disrupted social interactions in mice and dysregulated normal expression of PSD genes during neural development. The fact that the same genetic effects were found in two different mouse strains and in several brain regions increased potential for translation. The genetic and functional relationship between PSD and abnormal neurobehavioral disorders is well established, and our data suggest that BPA may contribute in a transgenerational manner to neurodevelopmental diseases.

Genetic variability in scaffolding proteins and risk for schizophrenia and autism-spectrum disorders: a systematic review

Scaffolding proteins represent an evolutionary solution to controlling the specificity of information transfer in intracellular networks. They are highly concentrated in complexes located in specific subcellular locations. One of these complexes is the postsynaptic density of the excitatory synapses. There, scaffolding proteins regulate various processes related to synaptic plasticity, such as glutamate receptor trafficking and signalling, and dendritic structure and function. Most scaffolding proteins can be grouped into 4 main families: discs large (DLG), discs-large-associated protein (DLGAP), Shank and Homer. Owing to the importance of scaffolding proteins in postsynaptic density architecture, it is not surprising that variants in the genes that code for these proteins have been associated with neuropsychiatric diagnoses, including schizophrenia and autism-spectrum disorders. Such evidence, together with the clinical, neurobiological and genetic overlap described between schizophrenia and autism-spectrum disorders, suggest that alteration of scaffolding protein dynamics could be part of the pathophysiology of both. However, despite the potential importance of scaffolding proteins in these psychiatric conditions, no systematic review has integrated the genetic and molecular data from studies conducted in the last decade. This review has the following goals: to systematically analyze the literature in which common and/or rare genetic variants (single nucleotide polymorphisms, single nucleotide variants and copy number variants) in the scaffolding family genes are associated with the risk for either schizophrenia or autism-spectrum disorders; to explore the implications of the reported genetic variants for gene expression and/or protein function; and to discuss the relationship of these genetic variants to the shared genetic, clinical and cognitive traits of schizophrenia and autism-spectrum disorders.

The ubiquitin proteasome pathway in neuropsychiatric disorders

The ubiquitin proteasome system (UPS) is a highly conserved pathway that tightly regulates protein turnover in cells. This process is integral to neuronal development, differentiation, and function. Several members of the UPS are disrupted in neuropsychiatric disorders, highlighting the importance of this pathway in brain development and function. In this review, we discuss some of these pathway members, the molecular processes they regulate, and the potential for targeting the UPS in an effort to develop therapeutic strategies in neuropsychiatric and neurodevelopmental disorders.

Association studies of WD repeat domain 3 and chitobiosyldiphosphodolichol beta-mannosyltransferase genes with schizophrenia in a Japanese population

Schizophrenia and schizophrenia-like symptoms induced by the dopamine agonists and N-methyl-D aspartate type glutamate receptor antagonists occur only after the adolescent period. Similarly, animal models of schizophrenia by these drugs are also induced after the critical period around postnatal week three. Based upon the development-dependent onsets of these psychotomimetic effects, by using a DNA microarray technique, we identified the WD repeat domain 3 (WDR3) and chitobiosyldiphosphodolichol beta-mannosyltransferase (ALG1) genes as novel candidates for schizophrenia-related molecules, whose mRNAs were up-regulated in the adult (postnatal week seven), but not in the infant (postnatal week one) rats by an indirect dopamine agonist, and phencyclidine, an antagonist of the NMDA receptor. WDR3 and other related proteins are the nuclear proteins presumably involved in various cellular activities, such as cell cycle progression, signal transduction, apoptosis, and gene regulation. ALG1 is presumed to be involved in the regulation of the protein N-glycosylation. To further elucidate the molecular pathophysiology of schizophrenia, we have evaluated the genetic association of WDR3 and ALG1 in schizophrenia. We examined 21 single nucleotide polymorphisms [SNPs; W1 (rs1812607)-W16 (rs6656360), A1 (rs8053916)-A10 (rs9673733)] from these genes using the Japanese case-control sample (1,808 schizophrenics and 2,170 matched controls). No significant genetic associations of these SNPs were identified. However, we detected a significant association of W4 (rs319471) in the female schizophrenics (allelic P = 0.003, genotypic P = 0.008). Based on a haplotype analysis, the observed haplotypes consisting of W4 (rs319471)–W5 (rs379058) also displayed a significant association in the female schizophrenics (P = 0.016). Even after correction for multiple testing, these associations remained significant. Our findings suggest that the WDR3 gene may likely be a sensitive factor in female patients with schizophrenia, and that modification of the WDR3 signaling pathway warrants further investigation as to the pathophysiology of schizophrenia.

Integrative genomics of microglia implicates DLG4 (PSD95) in the white matter development of preterm infants

Preterm birth places infants in an adverse environment that leads to abnormal brain development and cerebral injury through a poorly understood mechanism known to involve neuroinflammation. In this study, we integrate human and mouse molecular and neuroimaging data to investigate the role of microglia in preterm white matter damage. Using a mouse model where encephalopathy of prematurity is induced by systemic interleukin-1β administration, we undertake gene network analysis of the microglial transcriptomic response to injury, extend this by analysis of protein-protein interactions, transcription factors and human brain gene expression, and translate findings to living infants using imaging genomics. We show that DLG4 (PSD95) protein is synthesised by microglia in immature mouse and human, developmentally regulated, and modulated by inflammation; DLG4 is a hub protein in the microglial inflammatory response; and genetic variation in DLG4 is associated with structural differences in the preterm infant brain. DLG4 is thus apparently involved in brain development and impacts inter-individual susceptibility to injury after preterm birth.Inflammation mediated by microglia plays a key role in brain injury associated with preterm birth, but little is known about the microglial response in preterm infants. Here, the authors integrate molecular and imaging data from animal models and preterm infants, and find that microglial expression of DLG4 plays a role.

Comprehensive association analysis of 27 genes from the GABAergic system in Japanese individuals affected with schizophrenia

Involvement of the gamma-aminobutyric acid (GABA)-ergic system in schizophrenia pathogenesis through disrupted neurodevelopment has been highlighted in numerous studies. However, the function of common genetic variants of this system in determining schizophrenia risk is unknown. We therefore tested the association of 375 tagged SNPs in genes derived from the GABAergic system, such as GABA_A receptor subunit genes, and GABA related genes (glutamate decarboxylase genes, GABAergic-marker gene, genes involved in GABA receptor trafficking and scaffolding) in Japanese schizophrenia case-control samples (n=2926; 1415 cases and 1511 controls). We observed nominal association of SNPs in nine GABA_A receptor subunit genes and the GPHN gene with schizophrenia, although none survived correction for study-wide multiple testing. Two SNPs located in the GABRA1 gene, rs4263535 (P_allele=0.002; uncorrected) and rs1157122 (P_allele=0.006; uncorrected) showed top hits, followed by rs723432 (P_allele=0.007; uncorrected) in the GPHN gene. All three were significantly associated with schizophrenia and survived gene-wide multiple testing. Haplotypes containing associated variants in GABRA1 but not GPHN were significantly associated with schizophrenia. To conclude, we provided substantiating genetic evidence for the involvement of the GABAergic system in schizophrenia susceptibility. These results warrant further investigations to replicate the association of GABRA1 and GPHN with schizophrenia and to discern the precise mechanisms of disease pathophysiology.

Postsynaptic Density-95 Isoform Abnormalities in Schizophrenia

BACKGROUND

Postsynaptic density-95 (PSD-95) protein expression is dysregulated in schizophrenia in a variety of brain regions. We have designed experiments to examine PSD-95 mRNA splice variant expression in the dorsolateral prefrontal cortex from subjects with schizophrenia.

METHODS

We performed quantitative PCR and western blot analysis to measure PSD-95 expression in schizophrenia vs control subjects, rodent haloperidol treatment studies, rodent postmortem interval studies, and GluN1 knockdown (KD) mice vs controls.

RESULTS

We found decreased mRNA expression of beta (t = 4.506, df = 383, P < .0001) and truncated (t = 3.378, df = 383, P = .0008) isoforms of PSD-95, whereas alpha was unchanged. Additionally, we found decreased PSD-95 protein expression in schizophrenia (t = 2.746, df = 71, P = .0076). We found no correlation between PSD-95 protein and alpha, beta, or truncated mRNA isoforms in schizophrenia. PSD-95 beta transcript was increased (t = 3.346, df = 14, P < .05) in the GluN1 KD mouse model of schizophrenia. There was an increase in PSD-95 alpha mRNA expression (t = 2.905, df = 16, P < .05) in rats following long-term haloperidol administration.

CONCLUSIONS

Our findings describe a unique pathophysiology of specific PSD-95 isoform dysregulation in schizophrenia, chronic neuroleptic treatment, and a genetic lesion mouse model of drastically reduced N-methyl-d-aspartate receptor (NMDAR) complex expression. These data indicate that regulation of PSD-95 is multifaceted, may be isoform specific, and biologically relevant for synaptic signaling function. Specifically, NMDAR-mediated synaptic remodeling, and α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid receptor trafficking and interaction may be impaired in schizophrenia by decreased PSD-95 beta and truncated expression (respectively). Further, increased PSD-95 beta transcript in the GluN1 KD mouse model poses a potential compensatory rescue of NMDAR-mediated function via increased postsynaptic throughput of the severely reduced GluN1 signal. Together, these data propose that disruption of excitatory signaling complexes through genetic (GluN1 KD), pharmacologic (antipsychotics), or disease (schizophrenia) mechanisms specifically dysregulates PSD-95 expression.

Integrative approach for inference of gene regulatory networks using lasso-based random featuring and application to psychiatric disorders

Background

Inferring gene regulatory networks is one of the most interesting research areas in the systems biology. Many inference methods have been developed by using a variety of computational models and approaches. However, there are two issues to solve. First, depending on the structural or computational model of inference method, the results tend to be inconsistent due to innately different advantages and limitations of the methods. Therefore the combination of dissimilar approaches is demanded as an alternative way in order to overcome the limitations of standalone methods through complementary integration. Second, sparse linear regression that is penalized by the regularization parameter (lasso) and bootstrapping-based sparse linear regression methods were suggested in state of the art methods for network inference but they are not effective for a small sample size data and also a true regulator could be missed if the target gene is strongly affected by an indirect regulator with high correlation or another true regulator.

Results

We present two novel network inference methods based on the integration of three different criteria, (i) z-score to measure the variation of gene expression from knockout data, (ii) mutual information for the dependency between two genes, and (iii) linear regression-based feature selection.

Based on these criterion, we propose a lasso-based random feature selection algorithm (LARF) to achieve better performance overcoming the limitations of bootstrapping as mentioned above.

Conclusions

In this work, there are three main contributions. First, our z score-based method to measure gene expression variations from knockout data is more effective than similar criteria of related works. Second, we confirmed that the true regulator selection can be effectively improved by LARF. Lastly, we verified that an integrative approach can clearly outperform a single method when two different methods are effectively jointed. In the experiments, our methods were validated by outperforming the state of the art methods on DREAM challenge data, and then LARF was applied to inferences of gene regulatory network associated with psychiatric disorders.

PSD-95 regulates CRFR1 localization, trafficking and β-arrestin2 recruitment

Corticotropin-releasing factor (CRF) is a neuropeptide commonly associated with the hypothalamic-pituitary adrenal axis stress response. Upon release, CRF activates two G protein-coupled receptors (GPCRs): CRF receptor 1 (CRFR1) and CRF receptor 2 (CRFR2). Although both receptors contribute to mood regulation, CRFR1 antagonists have demonstrated anxiolytic and antidepressant-like properties that may be exploited in the generation of new pharmacological interventions for mental illnesses. Previous studies have demonstrated CRFR1 capable of heterologously sensitizing serotonin 2A receptor (5-HT2AR) signaling: another GPCR implicated in psychiatric disease. Interestingly, this phenomenon was dependent on Postsynaptic density 95 (PSD-95)/Disc Large/Zona Occludens (PDZ) interactions on the distal carboxyl termini of both receptors. In the current study, we demonstrate that endogenous PSD-95 can be co-immunoprecipitated with CRFR1 from cortical brain homogenate, and this interaction appears to be primarily via the PDZ-binding motif. Additionally, PSD-95 colocalizes with CRFR1 within the dendritic projections of cultured mouse neurons in a PDZ-binding motif-dependent manner. In HEK 293 cells, PSD-95 overexpression inhibited CRFR1 endocytosis, whereas PSD-95 shRNA knockdown enhanced CRFR1 endocytosis. Although PSD-95 does not appear to play a significant role in CRF-mediated cAMP or ERK1/2 signaling, PSD-95 was demonstrated to suppress β-arrestin2 recruitment: providing a potential mechanism for PSD-95's inhibition of endocytosis. In revisiting previously documented heterologous sensitization, PSD-95 shRNA knockdown did not prevent CRFR1-mediated enhancement of 5-HT2AR signaling. In conclusion, we have identified and characterized a novel functional relationship between CRFR1 and PSD-95 that may have implications in the design of new treatment strategies for mental illness.

IRSp53/BAIAP2 in dendritic spine development, NMDA receptor regulation, and psychiatric disorders

IRSp53 (also known as BAIAP2) is a multi-domain scaffolding and adaptor protein that has been implicated in the regulation of membrane and actin dynamics at subcellular structures, including filopodia and lamellipodia. Accumulating evidence indicates that IRSp53 is an abundant component of the postsynaptic density at excitatory synapses and an important regulator of actin-rich dendritic spines. In addition, IRSp53 has been implicated in diverse psychiatric disorders, including autism spectrum disorders, schizophrenia, and attention deficit/hyperactivity disorder. Mice lacking IRSp53 display enhanced NMDA (N-methyl-d-aspartate) receptor function accompanied by social and cognitive deficits, which are reversed by pharmacological suppression of NMDA receptor function. These results suggest the hypothesis that defective actin/membrane modulation in IRSp53-deficient dendritic spines may lead to social and cognitive deficits through NMDA receptor dysfunction. This article is part of the Special Issue entitled 'Synaptopathy--from Biology to Therapy'.

Genetic analysis of the glyoxalase system in schizophrenia

Recent reports suggest that carbonyl stress might affect a subset of schizophrenia patients suffering from severe symptoms. Carbonyl stress protection is achieved by the glyoxalase system consisting of two enzymes, glyoxalase 1 and 2, which in humans are encoded by the genes GLO1 and HAGH, respectively. Glyoxalase 1 and 2 catalyze the detoxification of reactive alpha-oxoaldehydes such as glyoxal and methylglyoxal, which are particularly damaging components of carbonyl stress. Here, we investigated the role of the glyoxalase system in schizophrenia by performing association analyses of common genetic variants (n=12) in GLO1 and HAGH in a Japanese sample consisting of 2012 schizophrenia patients and 2170 healthy controls. We detected a nominally significant association with schizophrenia (p=0.020) of rs11859266, a SNP in the intronic region of HAGH. However, rs11859266 did not survive multiple testing (empirical p=0.091). The variants in HAGH, rs11859266 and rs3743852, showed significant associations with schizophrenia in males at allelic and genotype levels, which remained persistent after multiple testing with the exception of rs3743852 for the genotype model. We further measured the mRNA expression of both genes in postmortem brain, but did not detect any changes in transcript expression levels between case and control samples or in sex-specific comparisons. Therefore, our findings suggest that an explanation of elevated carbonyl stress in a substantial part (reported as ~20%) of patients with schizophrenia will require the examination of a much larger cohort to detect risk alleles with weak effect size and/or other risk factors.

Schizophrenia and cortical blindness: protective effects and implications for language

The repeatedly noted absence of case-reports of individuals with schizophrenia and congenital/early developed blindness has led several authors to argue that the latter can confer protective effects against the former. In this work, we present a number of relevant case-reports from different syndromes that show comorbidity of congenital and early blindness with schizophrenia. On the basis of these reports, we argue that a distinction between different types of blindness in terms of the origin of the visual deficit, cortical or peripheral, is crucial for understanding the observed patterns of comorbidity. We discuss the genetic underpinnings and the brain structures involved in schizophrenia and blindness, with insights from language processing, laying emphasis on the three structures that particularly stand out: the occipital cortex, the lateral geniculate nucleus (LGN), and the pulvinar. Last, we build on previous literature on the nature of the protective effects in order to offer novel insights into the nature of the protection mechanism from the perspective of the brain structures involved in each type of blindness.

Association study of H2AFZ with schizophrenia in a Japanese case-control sample

It is widely accepted that malfunction of the N-methyl-D-aspartate (NMDA)-type glutamate receptor may be involved in the pathophysiology of schizophrenia. Several recent microRNA (miRNA) studies have demonstrated that the expression of the glutamate system-related miR-132 and miR-212 is changed in postmortem schizophrenic brains. Here we attempted to obtain further insight into the relationships among schizophrenia, the NMDA receptor, the molecular cascades controlled by these miRNAs and commonly predicted target genes of the two miRNAs. We focused on the H2AFZ (encoding H2A histone family, member Z) gene, whose expression was shown in our screening study to be modified by a schizophrenomimetic NMDA antagonist, phencyclidine. By performing polymerase chain reaction with fluorescent signal detention using the TaqMan system, we examined four tag single nucleotide polymorphisms (SNPs; SNP01-04) located around and within the H2AFZ gene for their genetic association with schizophrenia. The subjects were a Japanese cohort (2,012 patients with schizophrenia and 2,170 control subjects). We did not detect any significant genetic association of these SNPs with schizophrenia in this cohort. However, we observed a significant association of SNP02 (rs2276939) in the male patients with schizophrenia (allelic P = 0.003, genotypic P = 0.008). A haplotype analysis revealed that haplotypes consisting of SNP02-SNP03 (rs10014424)-SNP04 (rs6854536) also showed a significant association in the male patients with schizophrenia (P = 0.018). These associations remained significant even after correction for multiple testing. The present findings suggest that the H2AFZ gene may be a susceptibility factor in male subjects with schizophrenia, and that modification of the H2AFZ signaling pathway warrants further study in terms of the pathophysiology of schizophrenia.