Search

We performed a genome-wide association study of 6447 bipolar disorder (BD) cases and 12 639 controls from the International Cohort Collection for Bipolar Disorder (ICCBD). Meta-analysis was performed with prior results from the Psychiatric Genomics Consortium Bipolar Disorder Working Group for a combined sample of 13 902 cases and 19 279 controls. We identified eight genome-wide significant, associated regions, including a novel associated region on chromosome 10 (rs10884920; P=3.28 × 10−8) that includes the brain-enriched cytoskeleton protein adducin 3 (ADD3), a non-coding RNA, and a neuropeptide-specific aminopeptidase P (XPNPEP1). Our large sample size allowed us to test the heritability and genetic correlation of BD subtypes and investigate their genetic overlap with schizophrenia and major depressive disorder. We found a significant difference in heritability of the two most common forms of BD (BD I SNP-h2=0.35; BD II SNP-h2=0.25; P=0.02). The genetic correlation between BD I and BD II was 0.78, whereas the genetic correlation was 0.97 when BD cohorts containing both types were compared. In addition, we demonstrated a significantly greater load of polygenic risk alleles for schizophrenia and BD in patients with BD I compared with patients with BD II, and a greater load of schizophrenia risk alleles in patients with the bipolar type of schizoaffective disorder compared with patients with either BD I or BD II. These results point to a partial difference in the genetic architecture of BD subtypes as currently defined.

A number of large, rare copy number variants (CNVs) are deleterious for neurodevelopmental disorders, but large, rare, protective CNVs have not been reported for such phenotypes. Here we show in a CNV analysis of 47 005 individuals, the largest CNV analysis of schizophrenia to date, that large duplications (1.5-3.0 Mb) at 22q11.2--the reciprocal of the well-known, risk-inducing deletion of this locus--are substantially less common in schizophrenia cases than in the general population (0.014% vs 0.085%, OR=0.17, P=0.00086). 22q11.2 duplications represent the first putative protective mutation for schizophrenia.

The gap junction subunit connexin permits direct intercellular exchange of ions and molecules including glutamate, and plays an important role in the central nervous system. The connexin 40 (Cx40) and connexin 50 (Cx50) genes are located on chromosome 1q21.1, a region strongly linked with schizophrenia. These lines of evidence suggest that Cx40 and Cx50 may play a role in schizophrenia.Using an allele-specific PCR assay, four polymorphisms each were genotyped for Cx40 and Cx50 in 190 Caucasian patients with schizophrenia and 190 controls matched for sex, age and ethnicity. Following up, Cx50 rs989192 and rs4950495 were investigated in 99 Canadian and 163 Portuguese trios and nuclear families with schizophrenia probands. Hardy-Weinberg equilibrium and linkage disequilibrium (LD) block identification was carried out with HaploView, and association analysis for alleles and haplotypes with a permutation test of 10 000 simulations was carried out using the UNPHASED software program.Distributions of genotype frequencies of all markers were in Hardy-Weinberg equilibrium in Caucasian patients, controls and families. One rs989192-rs4950495 LD block was found in patients but not in controls. We found a significant association between the Cx50 rs989192-rs4950495 haplotype and schizophreniay (chi(2) = 29.55, p0.01). The A-C haplotype had a higher frequency in patients (chi(2) = 7.153, p0.01). Family studies also showed that the A-C haplotype was transmitted more often to patients with schizophrenia (chi(2) = 8.43, p0.01). No association of Cx40 with schizophrenia was found for allele, genotype or haplotype analyses.Our matched case-control and family study indicate that Cx50, but not Cx40, may play a role in the genetic susceptibility to schizophrenia.

We previously performed a genome-wide linkage scan in Portuguese schizophrenia families that identified a risk locus on chromosome 5q31-q35. This finding was supported by meta-analysis of 20 other schizophrenia genome-wide scans that identified 5q23.2-q34 as the second most compelling susceptibility locus in the genome. In the present report, we took a two-stage candidate gene association approach to investigate a group of gamma-aminobutyric acid (GABA) A receptor subunit genes (GABRA1, GABRA6, GABRB2, GABRG2, and GABRP) within our linkage peak. These genes are plausible candidates based on prior evidence for GABA system involvement in schizophrenia. In the first stage, associations were detected in a Portuguese patient sample with single nucleotide polymorphisms (SNPs) and haplotypes in GABRA1 (P=0.00062-0.048), GABRP (P=0.0024-0.042), and GABRA6 (P=0.0065-0.0088). The GABRA1 and GABRP findings were replicated in the second stage in an independent German family-based sample (P=0.0015-0.043). Supportive evidence for association was also obtained for a previously reported GABRB2 risk haplotype. Exploratory analyses of the effects of associated GABRA1 haplotypes on transcript levels found altered expression of GABRA6 and coexpressed genes of GABRA1 and GABRB2. Comparison of transcript levels in schizophrenia patients and unaffected siblings found lower patient expression of GABRA6 and coexpressed genes of GABRA1. Interestingly, the GABRA1 coexpressed genes include synaptic and vesicle-associated genes previously found altered in schizophrenia prefrontal cortex. Taken together, these results support the involvement of the chromosome 5q GABAA receptor gene cluster in schizophrenia, and suggest that schizophrenia-associated haplotypes may alter expression of GABA-related genes.

Large and rare copy number variants (CNVs) at several loci have been shown to increase risk for schizophrenia. Aiming to discover novel susceptibility CNV loci, we analyzed 6882 cases and 11 255 controls genotyped on Illumina arrays, most of which have not been used for this purpose before. We identified genes enriched for rare exonic CNVs among cases, and then attempted to replicate the findings in additional 14 568 cases and 15 274 controls. In a combined analysis of all samples, 12 distinct loci were enriched among cases with nominal levels of significance (P < 0.05); however, none would survive correction for multiple testing. These loci include recurrent deletions at 16p12.1, a locus previously associated with neurodevelopmental disorders (P = 0.0084 in the discovery sample and P = 0.023 in the replication sample). Other plausible candidates include non-recurrent deletions at the glutamate transporter gene SLC1A1, a CNV locus recently suggested to be involved in schizophrenia through linkage analysis, and duplications at 1p36.33 and CGNL1. A burden analysis of large (>500 kb), rare CNVs showed a 1.2% excess in cases after excluding known schizophrenia-associated loci, suggesting that additional susceptibility loci exist. However, even larger samples are required for their discovery.

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, r2 = 0.008; rs10043986–rs4704591, r2 = 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.

Schizophrenia is a common, multigenic psychiatric disorder. Linkage studies, including a recent meta-analysis of genome scans, have repeatedly implicated chromosome 8p12-p23.1 in schizophrenia susceptibility. More recently, significant association with a candidate gene on 8p12, neuregulin 1 (NRG1), has been reported in several European and Chinese samples. We investigated NRG1 for association in schizophrenia patients of Portuguese descent to determine whether this gene is a risk factor in this population. We tested NRG1 markers and haplotypes for association in 111 parent-proband trios, 321 unrelated cases, and 242 control individuals. Associations were found with a haplotype that overlaps the risk haplotype originally reported in the Icelandic population ("Hap(ICE)"), and two haplotypes located in the 3' end of NRG1 (all P0.05). However, association was not detected with Hap(ICE) itself. Comparison of NRG1 transcript expression in peripheral leukocytes from schizophrenia patients and unaffected siblings identified 3.8-fold higher levels of the SMDF variant in patients (P=0.039). Significant positive correlations (P0.001) were found between SMDF and HRG-beta 2 expression and between HRG-gamma and ndf43 expression, suggesting common transcriptional regulation of NRG1 variants. In summary, our results suggest that haplotypes across NRG1 and multiple NRG1 variants are involved in schizophrenia.

Linkage disequilibrium and haplotype structure of five GABA A receptor subunit genes investigated for association with schizophrenia