Identification of putative second genetic hits in schizophrenia carriers of high-risk copy number variants and resequencing in additional samples.

Copy number variants (CNVs) conferring risk of schizophrenia present incomplete penetrance, suggesting the existence of second genetic hits. Identification of second hits may help to find genes with rare variants of susceptibility to schizophrenia. The aim of this work was to search for second hits of moderate/high risk in schizophrenia carriers of risk CNVs and resequencing of the relevant genes in additional samples. To this end, ten patients with risk CNVs at cytobands 15q11.2, 15q11.2-13.1, 16p11.2, or 16p13.11, were subjected to whole-exome sequencing. Rare single nucleotide variants, defined as those absent from main public databases, were classified according to bioinformatic prediction of pathogenicity by CADD scores. The average number of rare predicted pathogenic variants per sample was 13.6 (SD 2.01). Two genes, BFAR and SYNJ1, presented rare predicted pathogenic variants in more than one sample. Follow-up resequencing of these genes in 432 additional cases and 432 controls identified a significant excess of rare predicted pathogenic variants in case samples at SYNJ1. Taking into account its function in clathrin-mediated synaptic vesicle endocytosis at presynaptic terminals, our results suggest an impairment of this process in schizophrenia.

Targeted resequencing of regulatory regions at schizophrenia risk loci: Role of rare functional variants at chromatin repressive states.

There is mounting evidence that regulatory variation plays an important role in genetic risk for schizophrenia. Here, we specifically search for regulatory variants at risk by sequencing promoter regions of twenty-three genes implied in schizophrenia by copy number variant or genome-wide association studies. After strict quality control, a total of 55,206bp per sample were analyzed in 526 schizophrenia cases and 516 controls from Galicia, NW Spain, using the Applied Biosystems SOLiD System. Variants were filtered based on frequency from public databases, chromatin states from the RoadMap Epigenomics Consortium at tissues relevant for schizophrenia, such as fetal brain, mid-frontal lobe, and angular gyrus, and prediction of functionality from RegulomeDB. The proportion of rare variants at polycomb repressive chromatin state at relevant tissues was higher in cases than in controls. The proportion of rare variants with predicted regulatory role was significantly higher in cases than in controls (P=0.0028, OR=1.93, 95% C.I.=1.23-3.04). Combination of information from both sources led to the identification of an excess of carriers of rare variants with predicted regulatory role located at polycomb repressive chromatin state at relevant tissues in cases versus controls (P=0.0016, OR=19.34, 95% C.I.=2.45-2495.26). The variants are located at two genes affected by the 17q12 copy number variant, LHX1 and HNF1B. These data strongly suggest that a specific epigenetic mechanism, chromatin remodeling by histone modification during early development, may be impaired in a subset of schizophrenia patients, in agreement with previous data.

Resequencing and association analysis of coding regions at twenty candidate genes suggest a role for rare risk variation at AKAP9 and protective variation at NRXN1 in schizophrenia susceptibility.

A fraction of genetic risk to develop schizophrenia may be due to low-frequency variants. This multistep study attempted to find low-frequency variants of high effect at coding regions of eleven schizophrenia susceptibility genes supported by genome-wide association studies (GWAS) and nine genes for the DISC1 interactome, a susceptibility gene-set. During the discovery step, a total of 125 kb per sample were resequenced in 153 schizophrenia patients and 153 controls from Galicia (NW Spain), and the cumulative role of low-frequency variants at a gene or at the DISC1 gene-set were analyzed by burden and variance-based tests. Relevant results were meta-analyzed when appropriate data were available. In addition, case-only putative damaging variants were genotyped in a further 419 cases and 398 controls. The discovery step revealed a protective effect of rare missense variants at NRXN1, a result supported by meta-analysis (OR = 0.67, 95% CI: 0.47-0.94, P = 0.021, based on 3848 patients and 3896 controls from six studies). The follow-up step based on case-only putative damaging variants revealed a promising risk variant at AKAP9. This variant, K873R, reached nominal significance after inclusion of 240 additional Spanish controls from databases. The variant, located in an ADCY2 binding region, is absent from large public databases. Interestingly, GWAS revealed an association between common ADCY2 variants and bipolar disorder, a disorder with considerable genetic overlap with schizophrenia. These data suggest a role of rare missense variants at NRXN1 and AKAP9 in schizophrenia susceptibility, probably related to alteration of the excitatory/inhibitory synaptic balance, deserving further investigation.

An efficient screening method for simultaneous detection of recurrent copy number variants associated with psychiatric disorders.

Several recurrent copy number variants (CNVs) increasing risk to neuropsychiatric diseases have been identified in recent years. They show variable clinical expressivity, being associated with different disorders, and incomplete penetrance. However, due to its very low frequency, the full variety of clinical outcomes associated with each one of these CNVs is unknown. Current methods for detection of CNVs are labor intensive, expensive or not suitable for high throughput analysis. Quantitative interspecies competitive PCR linked to variant minisequencing and detection by mass-spectrometry may overcome these limitations. Here, we present two multiplex assays based on this method to screen for eleven psychiatric risk CNVs, such as 1q21, 16p11.2, 3q29, or 16p13.11 regions, among others. The assays were tested in our collection of 514 schizophrenia patients. Results were compared with MLPA at two CNVs. Additional positive results were confirmed by exome sequencing. A total of fourteen patients were CNV carriers. The method presents high sensitivity and specificity, showing its utility as a cheap, accurate, high throughput screening tool for recurrent CNVs. The method may be very useful for management of psychiatric patients as well as screening of different collections of samples to better identify the full spectrum of clinical variability.

Role of DISC1 interacting proteins in schizophrenia risk from genome-wide analysis of missense SNPs.

A balanced translocation affecting DISC1 cosegregates with several psychiatric disorders, including schizophrenia, in a Scottish family. DISC1 is a hub protein of a network of protein-protein interactions involved in multiple developmental pathways within the brain. Gene set-based analysis has been proposed as an alternative to individual analysis of single nucleotide polymorphisms (SNPs) to get information from genome-wide association studies. In this work, we tested for an overrepresentation of the DISC1 interacting proteins within the top results of our ranked list of genes based on our previous genome-wide association study of missense SNPs in schizophrenia. Our data set consisted of 5100 common missense SNPs genotyped in 476 schizophrenic patients and 447 control subjects from Galicia, NW Spain. We used a modification of the Gene Set Enrichment Analysis adapted for SNPs, as implemented in the GenGen software. The analysis detected an overrepresentation of the DISC1 interacting proteins (permuted P-value=0.0158), indicative of the role of this gene set in schizophrenia risk. We identified seven leading-edge genes, MACF1, UTRN, DST, DISC1, KIF3A, SYNE1, and AKAP9, responsible for the overrepresentation. These genes are involved in neuronal cytoskeleton organization and intracellular transport through the microtubule cytoskeleton, suggesting that these processes may be impaired in schizophrenia.

No evidence that major mtDNA European haplogroups confer risk to schizophrenia.

Previous studies suggest that genetic factors could be involved in mitochondrial dysfunction observed in schizophrenia (SZ), some of them claiming a role of mtDNA common variants (mtSNPs) and/or haplogroups (hgs) in developing this disorder. These studies, however, have mainly been undertaken on relatively small cohorts of patients and control individuals and most have not yet been replicated. To further analyze the role of mtSNPs in SZ risk, we have carried out the largest genotyping effort to date using two Spanish case-control samples comprising a total of 942 schizophrenic patients and 1,231 unrelated controls: 454 patients and 616 controls from Santiago de Compostela (Galicia) and 488 patients and 615 controls from Reus (Catalonia). A set of 25 mtSNPs representing main branches of the European mtDNA phylogeny were genotyped in the Galician cohort and a subset of 16 out of these 25 mtSNPs was genotyped in the Catalan cohort. These 16 common variants characterize the most common European branches of the mtDNA phylogeny. We did not observe any positive association of mtSNPs and hgs with SZ. We discuss several deficiencies of previous studies that might explain the false positive nature of previous findings, including the confounding effect of population sub-structure and deficient statistical methodologies. It is unlikely that mtSNPs defining the most common European mtDNA haplogroups are related to SZ.

Association study of nonsynonymous single nucleotide polymorphisms in schizophrenia.

BACKGROUND: Genome-wide association studies using several hundred thousand anonymous markers present limited statistical power. Alternatively, association studies restricted to common nonsynonymous single nucleotide polymorphisms (nsSNPs) have the advantage of strongly reducing the multiple testing problem, while increasing the probability of testing functional single nucleotide polymorphisms (SNPs). METHODS: We performed a case-control association study of common nsSNPs in Galician (northwest Spain) samples using the Affymetrix GeneChip Human 20k cSNP Kit, followed by a replication study of the more promising results. After quality control procedures, the discovery sample consisted of 5100 nsSNPs at minor allele frequency >5% analyzed in 476 schizophrenia patients and 447 control subjects. The replication sample consisted of 4069 cases and 15,128 control subjects of European origin. We also performed multilocus analysis, using aggregated scores of nsSNPs at liberal significance thresholds and cross-validation procedures. RESULTS: The 5 independent nsSNPs with false discovery rate q ≤ .25, as well as 13 additional nsSNPs at p < .01 and located in functional candidate genes, were genotyped in the replication samples. One SNP, rs13107325, located at the metal ions transporter gene SLC39A8, reached significance in the combined sample after Bonferroni correction (trend test, p = 2.7 × 10(-6), allelic odds ratio = 1.32). This SNP presents minor allele frequency of 5% to 10% in many European populations but is rare outside Europe. We also confirmed the polygenic component of susceptibility. CONCLUSIONS: Taking into account that another metal ions transporter gene, SLC39A3, is associated to bipolar disorder, our findings reveal a role for brain metal homeostasis in psychosis.

Interaction between COMT haplotypes and cannabis in schizophrenia: a case-only study in two samples from Spain.

Cannabis use is one of the environmental factors with more solid evidence contributing to schizophrenia risk, especially in genetically susceptible individuals. One of the genes that may interact with cannabis is COMT, although available data are scarce. Here, we present a case-only study of the putative COMT-cannabis interaction in schizophrenia. Two Spanish samples from Santiago de Compostela and Valencia were screened for cannabis use. One hundred and fifty five individuals from a total of 748 patients were identified as cannabis users. Five SNPs in COMT, defining three common functional haplotypes with different enzymatic activities, were genotyped and analyzed for association at the SNP, haplotype and genotype levels. An association was detected between cannabis use and low activity variants (P<0.01) in the joint analysis and results were consistent between the two samples. Schizophrenic subjects homozygous for the Met allele at rs4680 doubled the probability of lifetime prevalence of cannabis use in comparison to Val homozygous (Mantel-Haenszel OR=2.07, 95% CI: 1.27-3.26, P=0.0031, in the combined sample). These data are in contrast to those from Caspi et al. (Biol. Psychiatry 57 (2005)1117-1127) who found association between schizophrenia/schizophreniform disorder and homozygosity at the high activity Val variant of rs4680. The results of our study are discussed in the context of previous findings, suggesting the involvement of COMT polymorphisms in the association between cannabis use and schizophrenia as well as the existence of additional factors mediating this association. However, further research is needed to confirm the COMT-cannabis interaction.

Heterozygosity at catechol-O-methyltransferase Val158Met and schizophrenia: new data and meta-analysis.

Catechol-O-methyltransferase (COMT) has been largely studied in relation to schizophrenia susceptibility. Most studies focused on the functional single nucleotide polymorphism (SNP) rs4680 that causes a substitution of Val by Met at codon 158 of the COMT protein. Recent meta-analyses do not support an association between allelic variants at rs4680 and schizophrenia. However, the putative role of overdominance has not been tested in meta-analyses, despite its biological plausibility. In this work, we tested the overdominant model in two Spanish samples (from Valencia and Santiago de Compostela), representing a total of 762 schizophrenic patients and 1042 controls, and performed a meta-analysis of the available studies under this model. A total of 51 studies comprising 13,894 schizophrenic patients and 16,087 controls were included in the meta-analysis, that revealed a small but significant protective effect for heterozygosity at rs4680 (pooled OR=0.947, P=0.023). Post-hoc analysis on southwestern European samples suggested a stronger effect in these populations (pooled OR=0.813, P=0.0009). Thus, the COMT functional polymorphism rs4680 contributes to schizophrenia genetic susceptibility under an overdominant model, indicating that both too high and too low levels of dopamine (DA) signalling may be risk factors. This effect can be modulated by genetic background.

Testing the antagonistic pleiotropy model of schizophrenia susceptibility by analysis of DAOA, PPP1R1B, and APOL1 genes.

Schizophrenia is a common disease associated with reduced fertility. Therefore, the existence of common susceptibility alleles not removed by natural selection may be considered an evolutionary paradox. The antagonistic pleiotropy model, proposed to explain this paradox, states that an allele may be common because of its overall selective advantage, in spite of deleterious effects on specific traits. Recent work on DAOA, PPP1R1B, and APOL1 suggests that these genes present common alleles associated to increase risk of schizophrenia but conferring an overall selective advantage, related to better cognitive performance (DAOA and PPP1R1B) or protection against pathogens (APOL1). To test if these genes fit the antagonistic pleiotropy model, we searched for recent natural selection at these loci applying the long-range haplotype test on data from the HapMap Project; and performed case-control association analysis in a well-powered sample, including 301 schizophrenic patients and 604 controls from Spain. For DAOA and PPP1R1B, we genotyped the Single-nucleotide polymorphisms (SNPs) needed to replicate previous associations, while for APOL1, we genotyped 15 tagSNPs, and seven putative functional SNPs. We did not detect evidence of recent natural selection. Furthermore, we did not find significant associations. Thus, these genes do not fit the antagonistic pleiotropy model.