MIR137 polygenic risk for schizophrenia and ephrin-regulated pathway: Role in lateral ventricles and corpus callosum volume.

Background/Objective. Enlarged lateral ventricle (LV) volume and decreased volume in the corpus callosum (CC) are hallmarks of schizophrenia (SZ). We previously showed an inverse correlation between LV and CC volumes in SZ, with global functioning decreasing with increased LV volume. This study investigates the relationship between LV volume, CC abnormalities, and the microRNA MIR137 and its regulated genes in SZ, because of MIR137's essential role in neurodevelopment. Methods. Participants were 1224 SZ probands and 1466 unaffected controls from the GENUS Consortium. Brain MRI scans, genotype, and clinical data were harmonized across cohorts and employed in the analyses. Results. Increased LV volumes and decreased CC central, mid-anterior, and mid-posterior volumes were observed in SZ probands. The MIR137-regulated ephrin pathway was significantly associated with CC:LV ratio, explaining a significant proportion (3.42 %) of CC:LV variance, and more than for LV and CC separately. Other pathways explained variance in either CC or LV, but not both. CC:LV ratio was also positively correlated with Global Assessment of Functioning, supporting previous subsample findings. SNP-based heritability estimates were higher for CC central:LV ratio (0.79) compared to CC or LV separately. Discussion. Our results indicate that the CC:LV ratio is highly heritable, influenced in part by variation in the MIR137-regulated ephrin pathway. Findings suggest that the CC:LV ratio may be a risk indicator in SZ that correlates with global functioning.

Evidence supporting the use of a brief cognitive assessment in routine clinical assessment for psychosis.

Cognitive impairment is a core feature of psychosis. Full cognitive assessments are not often conducted in routine clinical practice as administration is time-consuming. Here, we investigated whether brief tests of cognition could be used to predict broader neurocognitive performance in a manner practical for screening use in mental health services. We carried out a principal component analysis (PCA) to obtain an estimate of general cognitive function (N = 415). We investigated whether brief tests of memory accounted for a significant percentage of variation in the PCA scores. We used discriminant function analysis to determine if measures could predict classification as lower, intermediate or higher level of cognitive function and to what extent these groups overlapped with groups based on normative data. Memory tests correctly classified 65% of cases in the highest scoring group, 35% of cases in the intermediate group, and 77% of cases in the lowest scoring group. These PCA-derived groups and groups based on normative scores for the two tests were significantly associated (χ2 = 164.00, p < 0.001). These measures accurately identified three quarters of the low performing group, the group of greatest interest from the perspective of identifying those likely to need greater supports as part of clinical care. In so doing they suggest a potentially useful approach to screening for cognitive impairment in clinical services, upon which further assessment can be built if required.

Esophageal IgG4 levels correlate with histopathologic and transcriptomic features in eosinophilic esophagitis.

BACKGROUND: Recent data associate eosinophilic esophagitis (EoE) with IgG4 rather than IgE, but its significance and function have not been determined. Our aims were to measure esophageal IgG4 levels and to determine functional correlations as assessed by histologic and transcriptome analyses. METHODS: This case-control study included pediatric subjects with EoE (≥15 eosinophils/HPF) and non-EoE controls. Protein lysates were analyzed for IgA, IgM, and IgG1-IgG4 using the Luminex 100 system; IgE was quantified by ELISA. Esophageal biopsies were scored using the EoE histology scoring system. Transcripts were probed by the EoE diagnostic panel, designed to examine the expression of 96 esophageal transcripts. RESULTS: Esophageal IgG subclasses, IgA, and IgM, but not IgE, were increased in subjects with EoE relative to controls. The greatest change between groups was seen in IgG4 (4.2 mg/g protein [interquartile range: 1.0-13.1 mg/g protein] vs 0.2 mg/g protein [0.1-0.9]; P < .0001). Tissue IgG4 levels correlated with esophageal eosinophil counts (P = .0006); histologic grade (P = .0011) and stage (P = .0112) scores; and IL4, IL10, IL13, but not TGFB1, expression and had strong associations with a subset of the EoE transcriptome. Esophageal IgG4 transcript expression was increased and correlated with IgG4 protein levels and IL10 expression. CONCLUSION: These findings extend prior studies on IgG4 in adult EoE to the pediatric population and provide deeper understanding of the potential significance and regulation of IgG4, demonstrating that IgG4 is a relevant feature of the disease; is closely related to esophageal eosinophil levels, type 2 immunity and T regulatory cytokines; and is likely produced locally.

Genetically predicted complement component 4A expression: effects on memory function and middle temporal lobe activation.

BACKGROUND: The longstanding association between the major histocompatibility complex (MHC) locus and schizophrenia (SZ) risk has recently been accounted for, partially, by structural variation at the complement component 4 (C4) gene. This structural variation generates varying levels of C4 RNA expression, and genetic information from the MHC region can now be used to predict C4 RNA expression in the brain. Increased predicted C4A RNA expression is associated with the risk of SZ, and C4 is reported to influence synaptic pruning in animal models. METHODS: Based on our previous studies associating MHC SZ risk variants with poorer memory performance, we tested whether increased predicted C4A RNA expression was associated with reduced memory function in a large (n = 1238) dataset of psychosis cases and healthy participants, and with altered task-dependent cortical activation in a subset of these samples. RESULTS: We observed that increased predicted C4A RNA expression predicted poorer performance on measures of memory recall (p = 0.016, corrected). Furthermore, in healthy participants, we found that increased predicted C4A RNA expression was associated with a pattern of reduced cortical activity in middle temporal cortex during a measure of visual processing (p < 0.05, corrected). CONCLUSIONS: These data suggest that the effects of C4 on cognition were observable at both a cortical and behavioural level, and may represent one mechanism by which illness risk is mediated. As such, deficits in learning and memory may represent a therapeutic target for new molecular developments aimed at altering C4's developmental role.

MiR-137-derived polygenic risk: effects on cognitive performance in patients with schizophrenia and controls.

Variants at microRNA-137 (MIR137), one of the most strongly associated schizophrenia risk loci identified to date, have been associated with poorer cognitive performance. As microRNA-137 is known to regulate the expression of ~1900 other genes, including several that are independently associated with schizophrenia, we tested whether this gene set was also associated with variation in cognitive performance. Our analysis was based on an empirically derived list of genes whose expression was altered by manipulation of MIR137 expression. This list was cross-referenced with genome-wide schizophrenia association data to construct individual polygenic scores. We then tested, in a sample of 808 patients and 192 controls, whether these risk scores were associated with altered performance on cognitive functions known to be affected in schizophrenia. A subgroup of healthy participants also underwent functional imaging during memory (n=108) and face processing tasks (n=83). Increased polygenic risk within the empirically derived miR-137 regulated gene score was associated with significantly lower performance on intelligence quotient, working memory and episodic memory. These effects were observed most clearly at a polygenic threshold of P=0.05, although significant results were observed at all three thresholds analyzed. This association was found independently for the gene set as a whole, excluding the schizophrenia-associated MIR137 SNP itself. Analysis of the spatial working memory fMRI task further suggested that increased risk score (thresholded at P=10-5) was significantly associated with increased activation of the right inferior occipital gyrus. In conclusion, these data are consistent with emerging evidence that MIR137 associated risk for schizophrenia may relate to its broader downstream genetic effects.

Subcortical brain volume abnormalities in 2028 individuals with schizophrenia and 2540 healthy controls via the ENIGMA consortium.

The profile of brain structural abnormalities in schizophrenia is still not fully understood, despite decades of research using brain scans. To validate a prospective meta-analysis approach to analyzing multicenter neuroimaging data, we analyzed brain MRI scans from 2028 schizophrenia patients and 2540 healthy controls, assessed with standardized methods at 15 centers worldwide. We identified subcortical brain volumes that differentiated patients from controls, and ranked them according to their effect sizes. Compared with healthy controls, patients with schizophrenia had smaller hippocampus (Cohen's d=-0.46), amygdala (d=-0.31), thalamus (d=-0.31), accumbens (d=-0.25) and intracranial volumes (d=-0.12), as well as larger pallidum (d=0.21) and lateral ventricle volumes (d=0.37). Putamen and pallidum volume augmentations were positively associated with duration of illness and hippocampal deficits scaled with the proportion of unmedicated patients. Worldwide cooperative analyses of brain imaging data support a profile of subcortical abnormalities in schizophrenia, which is consistent with that based on traditional meta-analytic approaches. This first ENIGMA Schizophrenia Working Group study validates that collaborative data analyses can readily be used across brain phenotypes and disorders and encourages analysis and data sharing efforts to further our understanding of severe mental illness.

De novo mutations in schizophrenia implicate chromatin remodeling and support a genetic overlap with autism and intellectual disability.

Schizophrenia is a serious psychiatric disorder with a broadly undiscovered genetic etiology. Recent studies of de novo mutations (DNMs) in schizophrenia and autism have reinforced the hypothesis that rare genetic variation contributes to risk. We carried out exome sequencing on 57 trios with sporadic or familial schizophrenia. In sporadic trios, we observed a ~3.5-fold increase in the proportion of nonsense DNMs (0.101 vs 0.031, empirical P=0.01, Benjamini-Hochberg-corrected P=0.044). These mutations were significantly more likely to occur in genes with highly ranked probabilities of haploinsufficiency (P=0.0029, corrected P=0.006). DNMs of potential functional consequence were also found to occur in genes predicted to be less tolerant to rare variation (P=2.01 × 10(-)(5), corrected P=2.1 × 10(-)(3)). Genes with DNMs overlapped with genes implicated in autism (for example, AUTS2, CHD8 and MECP2) and intellectual disability (for example, HUWE1 and TRAPPC9), supporting a shared genetic etiology between these disorders. Functionally CHD8, MECP2 and HUWE1 converge on epigenetic regulation of transcription suggesting that this may be an important risk mechanism. Our results were consistent in an analysis of additional exome-based sequencing studies of other neurodevelopmental disorders. These findings suggest that perturbations in genes, which function in the epigenetic regulation of brain development and cognition, could have a central role in the susceptibility to, pathogenesis and treatment of mental disorders.

Effects of ZNF804A on auditory P300 response in schizophrenia.

The common variant rs1344706 within the zinc-finger protein gene ZNF804A has been strongly implicated in schizophrenia (SZ) susceptibility by a series of recent genetic association studies. Although associated with a pattern of altered neural connectivity, evidence that increased risk is mediated by an effect on cognitive deficits associated with the disorder has been equivocal. This study investigated whether the same ZNF804A risk allele was associated with variation in the P300 auditory-evoked response, a cognitively relevant putative endophenotype for SZ. We compared P300 responses in carriers and noncarriers of the ZNF804A risk allele genotype groups in Irish patients and controls (n=97). P300 response was observed to vary according to genotype in this sample, such that risk allele carriers showed relatively higher P300 response compared with noncarriers. This finding accords with behavioural data reported by our group and others. It is also consistent with the idea that ZNF804A may have an impact on cortical efficiency, reflected in the higher levels of activations required to achieve comparable behavioural accuracy on the task used.

Allelic differences between Europeans and Chinese for CREB1 SNPs and their implications in gene expression regulation, hippocampal structure and function, and bipolar disorder susceptibility.

Bipolar disorder (BD) is a polygenic disorder that shares substantial genetic risk factors with major depressive disorder (MDD). Genetic analyses have reported numerous BD susceptibility genes, while some variants, such as single-nucleotide polymorphisms (SNPs) in CACNA1C have been successfully replicated, many others have not and subsequently their effects on the intermediate phenotypes cannot be verified. Here, we studied the MDD-related gene CREB1 in a set of independent BD sample groups of European ancestry (a total of 64,888 subjects) and identified multiple SNPs significantly associated with BD (the most significant being SNP rs6785[A], P=6.32 × 10(-5), odds ratio (OR)=1.090). Risk SNPs were then subjected to further analyses in healthy Europeans for intermediate phenotypes of BD, including hippocampal volume, hippocampal function and cognitive performance. Our results showed that the risk SNPs were significantly associated with hippocampal volume and hippocampal function, with the risk alleles showing a decreased hippocampal volume and diminished activation of the left hippocampus, adding further evidence for their involvement in BD susceptibility. We also found the risk SNPs were strongly associated with CREB1 expression in lymphoblastoid cells (P<0.005) and the prefrontal cortex (P<1.0 × 10(-6)). Remarkably, population genetic analysis indicated that CREB1 displayed striking differences in allele frequencies between continental populations, and the risk alleles were completely absent in East Asian populations. We demonstrated that the regional prevalence of the CREB1 risk alleles in Europeans is likely caused by genetic hitchhiking due to natural selection acting on a nearby gene. Our results suggest that differential population histories due to natural selection on regional populations may lead to genetic heterogeneity of susceptibility to complex diseases, such as BD, and explain inconsistencies in detecting the genetic markers of these diseases among different ethnic populations.

Convergent lines of evidence support CAMKK2 as a schizophrenia susceptibility gene.

Genes that are differentially expressed between schizophrenia patients and healthy controls may have key roles in the pathogenesis of schizophrenia. We analyzed two large-scale genome-wide expression studies, which examined changes in gene expression in schizophrenia patients and their matched controls. We found calcium/calmodulin (CAM)-dependent protein kinase kinase 2 (CAMKK2) is significantly downregulated in individuals with schizophrenia in both studies. To seek the potential genetic variants that may regulate the expression of CAMKK2, we investigated the association between single-nucleotide polymorphisms (SNPs) within CAMKK2 and the expression level of CAMKK2. We found one SNP, rs1063843, which is located in intron 17 of CAMKK2, is strongly associated with the expression level of CAMKK2 in human brains (P=1.1 × 10(-6)) and lymphoblastoid cell lines (the lowest P=8.4 × 10(-6)). We further investigated the association between rs1063843 and schizophrenia in multiple independent populations (a total of 130 623 subjects) and found rs1063843 is significantly associated with schizophrenia (P=5.17 × 10(-5)). Interestingly, we found the T allele of rs1063843, which is associated with lower expression level of CAMKK2, has a higher frequency in individuals with schizophrenia in all of the tested samples, suggesting rs1063843 may be a causal variant. We also found that rs1063843 is associated with cognitive function and personality in humans. In addition, protein-protein interaction (PPI) analysis revealed that CAMKK2 participates in a highly interconnected PPI network formed by top schizophrenia genes, which further supports the potential role of CAMKK2 in the pathogenesis of schizophrenia. Taken together, these converging lines of evidence strongly suggest that CAMKK2 may have pivotal roles in schizophrenia susceptibility.

Excess of rare novel loss-of-function variants in synaptic genes in schizophrenia and autism spectrum disorders.

Schizophrenia (SZ) and autism spectrum disorders (ASDs) are complex neurodevelopmental disorders that may share an underlying pathology suggested by shared genetic risk variants. We sequenced the exonic regions of 215 genes in 147 ASD cases, 273 SZ cases and 287 controls, to identify rare risk mutations. Genes were primarily selected for their function in the synapse and were categorized as: (1) Neurexin and Neuroligin Interacting Proteins, (2) Post-synaptic Glutamate Receptor Complexes, (3) Neural Cell Adhesion Molecules, (4) DISC1 and Interactors and (5) Functional and Positional Candidates. Thirty-one novel loss-of-function (LoF) variants that are predicted to severely disrupt protein-coding sequence were detected among 2 861 rare variants. We found an excess of LoF variants in the combined cases compared with controls (P=0.02). This effect was stronger when analysis was limited to singleton LoF variants (P=0.0007) and the excess was present in both SZ (P=0.002) and ASD (P=0.001). As an individual gene category, Neurexin and Neuroligin Interacting Proteins carried an excess of LoF variants in cases compared with controls (P=0.05). A de novo nonsense variant in GRIN2B was identified in an ASD case adding to the growing evidence that this is an important risk gene for the disorder. These data support synapse formation and maintenance as key molecular mechanisms for SZ and ASD.

Neuropsychological effects of the CSMD1 genome-wide associated schizophrenia risk variant rs10503253.

The single-nucleotide polymorphism (SNP) rs10503253, located within the CUB and Sushi multiple domains-1 (CSMD1) gene on 8p23.2, was recently identified as genome-wide significant for schizophrenia (SZ), but is of unknown function. We investigated the neurocognitive effects of this CSMD1 variant in vivo in patients and healthy participants using behavioral and imaging measures of brain structure and function. We compared carriers and non-carriers of the risk 'A' allele on measures of neuropsychological performance typically impaired in SZ (general cognitive ability, episodic and working memory and attentional control) in independent samples of Irish patients (n = 387) and controls (n = 171) and German patients (205) and controls (n = 533). Across these groups, the risk 'A' allele at CSMD1 was associated with deleterious effects across a number of neurocognitive phenotypes. Specifically, the risk allele was associated with poorer performance on neuropsychological measures of general cognitive ability and memory function but not attentional control. These effects, while significant, were subtle, and varied between samples. Consistent with previous evidence suggesting that CSMD1 may be involved in brain mechanisms related to memory and learning, these data appear to reflect the deleterious effects of the identified 'A' risk allele on neurocognitive function, possibly as part of the mechanism by which CSMD1 is associated with SZ risk.

Mood congruent psychotic symptoms and specific cognitive deficits in carriers of the novel schizophrenia risk variant at MIR-137.

OBJECTIVE: The Schizophrenia Psychiatric Genome-wide Association (GWAS) Consortium recently reported on five novel schizophrenia susceptibility loci. The most significant finding mapped to a micro-RNA, MIR-137, which may be involved in regulating the function of other schizophrenia and bipolar disorder susceptibility genes. METHOD: We genotyped 821 patients with confirmed DSM-IV diagnoses of schizophrenia, bipolar affective disorder I and schizoaffective disorder for the risk SNP (rs1625579) and investigated the clinical profiles of risk allele carriers using a within-case design. We also assessed neurocognitive performance in a subset of cases (n=399) and controls (n=171). RESULTS: Carriers of the risk allele had lower scores for an OPCRIT-derived positive symptom factor (p=0.04) and lower scores on a lifetime measure of psychosis incongruity (p=0.017). Risk allele carriers also had more cognitive deficits involving episodic memory and attentional control. CONCLUSION: This is the first evidence that the MIR-137 risk variant may be associated with a specific subgroup of psychosis patients. Although the effect of this single SNP was not clinically relevant, investigation of the impact of carrying multiple risk SNPs in the MIR-137 regulatory network on diagnosis and illness profile may be warranted.

Preserved cognitive function is associated with suicidal ideation and single suicide attempts in schizophrenia.

OBJECTIVE: Suicide is the leading cause of death in schizophrenia. An association between suicidal behavior and both higher and lower cognitive ability in schizophrenia has been reported. To clarify this relationship, we investigated whether the relationship between suicidality and neurocognition varied according to differences in suicidal ideation and behavior. METHODS: Three hundred and ten patients with DSM-IV diagnoses of schizophrenia or schizoaffective disorder were categorized based on patient and staff interviews as either non-suicide attempters, non-attempters expressing suicidal ideation, single suicide attempters, or multiple suicide attempters. These groups were compared on a neuropsychological battery examining current general cognitive ability, episodic and working memory, and attentional control. RESULTS: Neuropsychological performance in those with a history of suicidal ideation (n=63), and those who had made one suicide attempt (n=48) was comparable. Together, these groups outperformed patients with no history of either suicidal behavior or ideation (n=172) on measures of IQ, episodic memory and working memory. Only differences in global cognition remained significant after controlling for between-group differences in depressive symptoms. Those who had either expressed suicidal ideation and/or made a single suicide attempt demonstrated trend level advantages in neuropsychological tests over those that had made multiple suicide attempts. DISCUSSION: These findings support earlier evidence of an association between suicidality and neurocognitive ability in schizophrenia. Specifically, these data suggest that patients who have contemplated suicide or made a single suicide attempt have better cognitive functioning than those who have not. Suicidality in multiple attempters, who do not perform better in neurocognitive tests than those who have neither contemplated nor attempted suicide, is likely to be influenced by factors other than neurocognitive ability.

Promoter polymorphisms in two overlapping 6p25 genes implicate mitochondrial proteins in cognitive deficit in schizophrenia.

In a previous study, we detected a 6p25-p24 region linked to schizophrenia in families with high composite cognitive deficit (CD) scores, a quantitative trait integrating multiple cognitive measures. Association mapping of a 10 Mb interval identified a 260 kb region with a cluster of single-nucleotide polymorphisms (SNPs) significantly associated with CD scores and memory performance. The region contains two colocalising genes, LYRM4 and FARS2, both encoding mitochondrial proteins. The two tagging SNPs with strongest evidence of association were located around the overlapping putative promoters, with rs2224391 predicted to alter a transcription factor binding site (TFBS). Sequencing the promoter region identified 22 SNPs, many predicted to affect TFBSs, in a tight linkage disequilibrium block. Luciferase reporter assays confirmed promoter activity in the predicted promoter region, and demonstrated marked downregulation of expression in the LYRM4 direction under the haplotype comprising the minor alleles of promoter SNPs, which however is not driven by rs2224391. Experimental evidence from LYRM4 expression in lymphoblasts, gel-shift assays and modelling of DNA breathing dynamics pointed to two adjacent promoter SNPs, rs7752203-rs4141761, as the functional variants affecting expression. Their C-G alleles were associated with higher transcriptional activity and preferential binding of nuclear proteins, whereas the G-A combination had opposite effects and was associated with poor memory and high CD scores. LYRM4 is a eukaryote-specific component of the mitochondrial biogenesis of Fe-S clusters, essential cofactors in multiple processes, including oxidative phosphorylation. LYRM4 downregulation may be one of the mechanisms involved in inefficient oxidative phosphorylation and oxidative stress, increasingly recognised as contributors to schizophrenia pathogenesis.

Fine mapping of ZNF804A and genome-wide significant evidence for its involvement in schizophrenia and bipolar disorder.

A recent genome-wide association study (GWAS) reported evidence for association between rs1344706 within ZNF804A (encoding zinc-finger protein 804A) and schizophrenia (P=1.61 × 10(-7)), and stronger evidence when the phenotype was broadened to include bipolar disorder (P=9.96 × 10(-9)). In this study we provide additional evidence for association through meta-analysis of a larger data set (schizophrenia/schizoaffective disorder N=18 945, schizophrenia plus bipolar disorder N=21 274 and controls N=38 675). We also sought to better localize the association signal using a combination of de novo polymorphism discovery in exons, pooled de novo polymorphism discovery spanning the genomic sequence of the locus and high-density linkage disequilibrium (LD) mapping. The meta-analysis provided evidence for association between rs1344706 that surpasses widely accepted benchmarks of significance by several orders of magnitude for both schizophrenia (P=2.5 × 10(-11), odds ratio (OR) 1.10, 95% confidence interval 1.07-1.14) and schizophrenia and bipolar disorder combined (P=4.1 × 10(-13), OR 1.11, 95% confidence interval 1.07-1.14). After de novo polymorphism discovery and detailed association analysis, rs1344706 remained the most strongly associated marker in the gene. The allelic association at the ZNF804A locus is now one of the most compelling in schizophrenia to date, and supports the accumulating data suggesting overlapping genetic risk between schizophrenia and bipolar disorder.

GWA study data mining and independent replication identify cardiomyopathy-associated 5 (CMYA5) as a risk gene for schizophrenia.

We conducted data-mining analyses using the Clinical Antipsychotic Trials of Intervention Effectiveness (CATIE) and molecular genetics of schizophrenia genome-wide association study supported by the genetic association information network (MGS-GAIN) schizophrenia data sets and performed bioinformatic prioritization for all the markers with P-values ≤0.05 in both data sets. In this process, we found that in the CMYA5 gene, there were two non-synonymous markers, rs3828611 and rs10043986, showing nominal significance in both the CATIE and MGS-GAIN samples. In a combined analysis of both the CATIE and MGS-GAIN samples, rs4704591 was identified as the most significant marker in the gene. Linkage disequilibrium analyses indicated that these markers were in low LD (3 828 611-rs10043986, r(2)=0.008; rs10043986-rs4704591, r(2)=0.204). In addition, CMYA5 was reported to be physically interacting with the DTNBP1 gene, a promising candidate for schizophrenia, suggesting that CMYA5 may be involved in the same biological pathway and process. On the basis of this information, we performed replication studies for these three single-nucleotide polymorphisms. The rs3828611 was found to have conflicting results in our Irish samples and was dropped out without further investigation. The other two markers were verified in 23 other independent data sets. In a meta-analysis of all 23 replication samples (family samples, 912 families with 4160 subjects; case-control samples, 11 380 cases and 15 021 controls), we found that both markers are significantly associated with schizophrenia (rs10043986, odds ratio (OR)=1.11, 95% confidence interval (CI)=1.04-1.18, P=8.2 × 10(-4) and rs4704591, OR=1.07, 95% CI=1.03-1.11, P=3.0 × 10(-4)). The results were also significant for the 22 Caucasian replication samples (rs10043986, OR=1.11, 95% CI=1.03-1.17, P=0.0026 and rs4704591, OR=1.07, 95% CI=1.02-1.11, P=0.0015). Furthermore, haplotype conditioned analyses indicated that the association signals observed at these two markers are independent. On the basis of these results, we concluded that CMYA5 is associated with schizophrenia and further investigation of the gene is warranted.