Rare coding variants in ten genes confer substantial risk for schizophrenia.

Rare coding variation has historically provided the most direct connections between gene function and disease pathogenesis. By meta-analysing the whole exomes of 24,248 schizophrenia cases and 97,322 controls, we implicate ultra-rare coding variants (URVs) in 10 genes as conferring substantial risk for schizophrenia (odds ratios of 3-50, P < 2.14 × 10-6) and 32 genes at a false discovery rate of <5%. These genes have the greatest expression in central nervous system neurons and have diverse molecular functions that include the formation, structure and function of the synapse. The associations of the NMDA (N-methyl-D-aspartate) receptor subunit GRIN2A and AMPA (α-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid) receptor subunit GRIA3 provide support for dysfunction of the glutamatergic system as a mechanistic hypothesis in the pathogenesis of schizophrenia. We observe an overlap of rare variant risk among schizophrenia, autism spectrum disorders1, epilepsy and severe neurodevelopmental disorders2, although different mutation types are implicated in some shared genes. Most genes described here, however, are not implicated in neurodevelopment. We demonstrate that genes prioritized from common variant analyses of schizophrenia are enriched in rare variant risk3, suggesting that common and rare genetic risk factors converge at least partially on the same underlying pathogenic biological processes. Even after excluding significantly associated genes, schizophrenia cases still carry a substantial excess of URVs, which indicates that more risk genes await discovery using this approach.

Common variants on chromosome 6p22.1 are associated with schizophrenia.

Schizophrenia, a devastating psychiatric disorder, has a prevalence of 0.5-1%, with high heritability (80-85%) and complex transmission. Recent studies implicate rare, large, high-penetrance copy number variants in some cases, but the genes or biological mechanisms that underlie susceptibility are not known. Here we show that schizophrenia is significantly associated with single nucleotide polymorphisms (SNPs) in the extended major histocompatibility complex region on chromosome 6. We carried out a genome-wide association study of common SNPs in the Molecular Genetics of Schizophrenia (MGS) case-control sample, and then a meta-analysis of data from the MGS, International Schizophrenia Consortium and SGENE data sets. No MGS finding achieved genome-wide statistical significance. In the meta-analysis of European-ancestry subjects (8,008 cases, 19,077 controls), significant association with schizophrenia was observed in a region of linkage disequilibrium on chromosome 6p22.1 (P = 9.54 x 10(-9)). This region includes a histone gene cluster and several immunity-related genes--possibly implicating aetiological mechanisms involving chromatin modification, transcriptional regulation, autoimmunity and/or infection. These results demonstrate that common schizophrenia susceptibility alleles can be detected. The characterization of these signals will suggest important directions for research on susceptibility mechanisms.

No significant association of 14 candidate genes with schizophrenia in a large European ancestry sample: implications for psychiatric genetics.

OBJECTIVE: The authors carried out a genetic association study of 14 schizophrenia candidate genes (RGS4, DISC1, DTNBP1, STX7, TAAR6, PPP3CC, NRG1, DRD2, HTR2A, DAOA, AKT1, CHRNA7, COMT, and ARVCF). This study tested the hypothesis of association of schizophrenia with common single nucleotide polymorphisms (SNPs) in these genes using the largest sample to date that has been collected with uniform clinical methods and the most comprehensive set of SNPs in each gene. METHOD: The sample included 1,870 cases (schizophrenia and schizoaffective disorder) and 2,002 screened comparison subjects (i.e. controls), all of European ancestry, with ancestral outliers excluded based on analysis of ancestry-informative markers. The authors genotyped 789 SNPs, including tags for most common SNPs in each gene, SNPs previously reported as associated, and SNPs located in functional domains of genes such as promoters, coding exons (including nonsynonymous SNPs), 3' untranslated regions, and conserved noncoding sequences. After extensive data cleaning, 648 SNPs were analyzed for association of single SNPs and of haplotypes. RESULTS: Neither experiment-wide nor gene-wide statistical significance was observed in the primary single-SNP analyses or in secondary analyses of haplotypes or of imputed genotypes for additional common HapMap SNPs. Results in SNPs previously reported as associated with schizophrenia were consistent with chance expectation, and four functional polymorphisms in COMT, DRD2, and HTR2A did not produce nominally significant evidence to support previous evidence for association. CONCLUSIONS: It is unlikely that common SNPs in these genes account for a substantial proportion of the genetic risk for schizophrenia, although small effects cannot be ruled out.

Genome-wide association study of clinical dimensions of schizophrenia: polygenic effect on disorganized symptoms.

OBJECTIVE: Multiple sources of evidence suggest that genetic factors influence variation in clinical features of schizophrenia. The authors present the first genome-wide association study (GWAS) of dimensional symptom scores among individuals with schizophrenia. METHOD: Based on the Lifetime Dimensions of Psychosis Scale ratings of 2,454 case subjects of European ancestry from the Molecular Genetics of Schizophrenia (MGS) sample, three symptom factors (positive, negative/disorganized, and mood) were identified with exploratory factor analysis. Quantitative scores for each factor from a confirmatory factor analysis were analyzed for association with 696,491 single-nucleotide polymorphisms (SNPs) using linear regression, with correction for age, sex, clinical site, and ancestry. Polygenic score analysis was carried out to determine whether case and comparison subjects in 16 Psychiatric GWAS Consortium (PGC) schizophrenia samples (excluding MGS samples) differed in scores computed by weighting their genotypes by MGS association test results for each symptom factor. RESULTS: No genome-wide significant associations were observed between SNPs and factor scores. Most of the SNPs producing the strongest evidence for association were in or near genes involved in neurodevelopment, neuroprotection, or neurotransmission, including genes playing a role in Mendelian CNS diseases, but no statistically significant effect was observed for any defined gene pathway. Finally, polygenic scores based on MGS GWAS results for the negative/disorganized factor were significantly different between case and comparison subjects in the PGC data set; for MGS subjects, negative/disorganized factor scores were correlated with polygenic scores generated using case-control GWAS results from the other PGC samples. CONCLUSIONS: The polygenic signal that has been observed in cross-sample analyses of schizophrenia GWAS data sets could be in part related to genetic effects on negative and disorganized symptoms (i.e., core features of chronic schizophrenia).

Copy number variants in schizophrenia: confirmation of five previous findings and new evidence for 3q29 microdeletions and VIPR2 duplications.

OBJECTIVE: To evaluate previously reported associations of copy number variants (CNVs) with schizophrenia and to identify additional associations, the authors analyzed CNVs in the Molecular Genetics of Schizophrenia study (MGS) and additional available data. METHOD: After quality control, MGS data for 3,945 subjects with schizophrenia or schizoaffective disorder and 3,611 screened comparison subjects were available for analysis of rare CNVs (<1% frequency). CNV detection thresholds were chosen that maximized concordance in 151 duplicate assays. Pointwise and genewise analyses were carried out, as well as analyses of previously reported regions. Selected regions were visually inspected and confirmed with quantitative polymerase chain reaction. RESULTS: In analyses of MGS data combined with other available data sets, odds ratios of 7.5 or greater were observed for previously reported deletions in chromosomes 1q21.1, 15q13.3, and 22q11.21, duplications in 16p11.2, and exon-disrupting deletions in NRXN1. The most consistently supported candidate associations across data sets included a 1.6-Mb deletion in chromosome 3q29 (21 genes, TFRC to BDH1) that was previously described in a mild-moderate mental retardation syndrome, exonic duplications in the gene for vasoactive intestinal peptide receptor 2 (VIPR2), and exonic duplications in C16orf72. The case subjects had a modestly higher genome-wide number of gene-containing deletions (>100 kb and >1 Mb) but not duplications. CONCLUSIONS: The data strongly confirm the association of schizophrenia with 1q21.1, 15q13.3, and 22q11.21 deletions, 16p11.2 duplications, and exonic NRXN1 deletions. These CNVs, as well as 3q29 deletions, are also associated with mental retardation, autism spectrum disorders, and epilepsy. Additional candidate genes and regions, including VIPR2, were identified. Study of the mechanisms underlying these associations should shed light on the pathophysiology of schizophrenia.

The Internet-based MGS2 control sample: self report of mental illness.

OBJECTIVE: The Molecular Genetics of Schizophrenia (MGS2) project recruited an adult control sample of non-Hispanic European-ancestry (N=3,364) and African American (N=1,301) subjects. METHOD: Subjects gave consent to deposit phenotypic data and blood samples into a repository for general research use, with full anonymization of the sample. The authors compared the control sample with population census data for demographic data and with previous population surveys for anthropometrics and prevalences of psychiatric disorders as estimated by an Internet-administered questionnaire. RESULTS: The full MGS2 control sample includes 4,665 subjects (European-ancestry: N=3,364; African American: N=1,301), of whom 3,626 were included in the MGS2 genome-wide association study (GWAS). The sample is generally demographically representative of the U.S. population, except for being older and more female, educated, and affluent, although all strata are represented. Self-reported ancestry was consistent with genotypic and census data. Lifetime prevalences for depressive, anxiety, and substance use diagnoses were higher than in previous population-based surveys, probably due to use of an abbreviated self-report instrument. However, patterns such as sex ratios, comorbidity, and demographic associations were consistent with previous reports. DNA quality for the Internet collected/evaluated control sample was comparable to that of the face-to-face case sample. CONCLUSIONS: The Internet-based methods facilitated the rapid collection of large and anonymized non-Hispanic European-ancestry and African American control samples that have been validated as being generally representative for many aspects of demography, ancestry, and morbidity. Utilization of clinical screening data shared with the scientific community may permit investigators to select appropriate controls for some studies.

Genome scans and gene expression microarrays converge to identify gene regulatory loci relevant in schizophrenia.

Multiple linkage regions have been reported in schizophrenia, and some appear to harbor susceptibility genes that are differentially expressed in postmortem brain tissue derived from unrelated individuals. We combined traditional genome-wide linkage analysis in a multiplex family with lymphocytic genome-wide expression analysis. A genome scan suggested linkage to a chromosome 4q marker (D4S1530, LOD 2.17, theta = 0) using a dominant model. Haplotype analysis using flanking microsatellite markers delineated a 14 Mb region that cosegregated with all those affected. Subsequent genome-wide scan with SNP genotypes supported the evidence of linkage to 4q33-35.1 (LOD = 2.39) using a dominant model. Genome-wide microarray analysis of five affected and five unaffected family members identified two differentially expressed genes within the haplotype AGA and GALNT7 (aspartylglucosaminidase and UDP-N-acetyl-alpha-D-galactosamine: polypeptide N-acetylgalactosaminyltransferase 7) with nominal significance; however, these genes did not remain significant following analysis of covariance. We carried out genome-wide linkage analyses between the quantitative expression phenotype and genetic markers. AGA expression levels showed suggestive linkage to multiple markers in the haplotype (maximum LOD = 2.37) but to no other genomic region. GALNT7 expression levels showed linkage to regulatory loci at 4q28.1 (maximum LOD = 3.15) and in the haplotype region at 4q33-35.1 (maximum LOD = 2.37). ADH1B (alcohol dehydrogenase IB) was linked to loci at 4q21-q23 (maximum LOD = 3.08) and haplotype region at 4q33-35.1 (maximum LOD = 2.27). Seven differentially expressed genes were validated with RT-PCR. Three genes in the 4q33-35.1 haplotype region were also differentially expressed in schizophrenia in postmortem dorsolateral prefrontal cortex: AGA, HMGB2, and SCRG1. These results indicate that combining differential gene expression with linkage analysis may help in identifying candidate genes and potential regulatory sites. Moreover, they also replicate recent findings of complex trans- and cis- regulation of genes.

Genomewide linkage scan of 409 European-ancestry and African American families with schizophrenia: suggestive evidence of linkage at 8p23.3-p21.2 and 11p13.1-q14.1 in the combined sample.

We report the clinical characteristics of a schizophrenia sample of 409 pedigrees--263 of European ancestry (EA) and 146 of African American ancestry (AA)--together with the results of a genome scan (with a simple tandem repeat polymorphism interval of 9 cM) and follow-up fine mapping. A family was required to have a proband with schizophrenia (SZ) and one or more siblings of the proband with SZ or schizoaffective disorder. Linkage analyses included 403 independent full-sibling affected sibling pairs (ASPs) (279 EA and 124 AA) and 100 all-possible half-sibling ASPs (15 EA and 85 AA). Nonparametric multipoint linkage analysis of all families detected two regions with suggestive evidence of linkage at 8p23.3-q12 and 11p11.2-q22.3 (empirical Z likelihood-ratio score [Z(lr)] threshold >/=2.65) and, in exploratory analyses, two other regions at 4p16.1-p15.32 in AA families and at 5p14.3-q11.2 in EA families. The most significant linkage peak was in chromosome 8p; its signal was mainly driven by the EA families. Z(lr) scores >2.0 in 8p were observed from 30.7 cM to 61.7 cM (Center for Inherited Disease Research map locations). The maximum evidence in the full sample was a multipoint Z(lr) of 3.25 (equivalent Kong-Cox LOD of 2.30) near D8S1771 (at 52 cM); there appeared to be two peaks, both telomeric to neuregulin 1 (NRG1). There is a paracentric inversion common in EA individuals within this region, the effect of which on the linkage evidence remains unknown in this and in other previously analyzed samples. Fine mapping of 8p did not significantly alter the significance or length of the peak. We also performed fine mapping of 4p16.3-p15.2, 5p15.2-q13.3, 10p15.3-p14, 10q25.3-q26.3, and 11p13-q23.3. The highest increase in Z(lr) scores was observed for 5p14.1-q12.1, where the maximum Z(lr) increased from 2.77 initially to 3.80 after fine mapping in the EA families.