Sarah E. Bergen

Genome-wide association analysis identifies 13 new risk loci for schizophrenia

Schizophrenia is an idiopathic mental disorder with a heritable component and a substantial public health impact. We conducted a multi-stage genome-wide association study (GWAS) for schizophrenia beginning with a Swedish national sample (5,001 cases and 6,243 controls) followed by meta-analysis with previous schizophrenia GWAS (8,832 cases and 12,067 controls) and finally by replication of SNPs in 168 genomic regions in independent samples (7,413 cases, 19,762 controls and 581 parent-offspring trios). We identified 22 loci associated at genome-wide significance; 13 of these are new, and 1 was previously implicated in bipolar disorder. Examination of candidate genes at these loci suggests the involvement of neuronal calcium signaling. We estimate that 8,300 independent, mostly common SNPs (95% credible interval of 6,300–10,200 SNPs) contribute to risk for schizophrenia and that these collectively account for at least 32% of the variance in liability. Common genetic variation has an important role in the etiology of schizophrenia, and larger studies will allow more detailed understanding of this disorder.

A polygenic burden of rare disruptive mutations in schizophrenia

Schizophrenia is a common disease with a complex aetiology, probably involving multiple and heterogeneous genetic factors. Here, by analysing the exome sequences of 2,536 schizophrenia cases and 2,543 controls, we demonstrate a polygenic burden primarily arising from rare (less than 1 in 10,000), disruptive mutations distributed across many genes. Particularly enriched gene sets include the voltage-gated calcium ion channel and the signalling complex formed by the activity-regulated cytoskeleton-associated scaffold protein (ARC) of the postsynaptic density, sets previously implicated by genome-wide association and copy-number variation studies. Similar to reports in autism, targets of the fragile X mental retardation protein (FMRP, product of FMR1) are enriched for case mutations. No individual gene-based test achieves significance after correction for multiple testing and we do not detect any alleles of moderately low frequency (approximately 0.5 to 1 per cent) and moderately large effect. Taken together, these data suggest that population-based exome sequencing can discover risk alleles and complements established gene-mapping paradigms in neuropsychiatric disease.

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.

Extremely low-coverage sequencing and imputation increases power for genome-wide association studies

Genome-wide association studies (GWAS) have proven to be a powerful method to identify common genetic variants contributing to susceptibility to common diseases. Here, we show that extremely low-coverage sequencing (0.1–0.5×) captures almost as much of the common (>5%) and low-frequency (1–5%) variation across the genome as SNP arrays. As an empirical demonstration, we show that genome-wide SNP genotypes can be inferred at a mean r2 of 0.71 using off-target data (0.24× average coverage) in a whole-exome study of 909 samples. Using both simulated and real exome-sequencing data sets, we show that association statistics obtained using extremely low-coverage sequencing data attain similar P values at known associated variants as data from genotyping arrays, without an excess of false positives. Within the context of reductions in sample preparation and sequencing costs, funds invested in extremely low-coverage sequencing can yield several times the effective sample size of GWAS based on SNP array data and a commensurate increase in statistical power.

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).

Genetic modifiers and subtypes in schizophrenia: Investigations of age at onset, severity, sex and family history

Schizophrenia is a genetically and clinically heterogeneous disorder. Genetic risk factors for the disorder may differ between the sexes or between multiply affected families compared to cases with no family history. Additionally, limited data support a genetic basis for variation in onset and severity, but specific loci have not been identified. We performed genome-wide association studies (GWAS) examining genetic influences on age at onset (AAO) and illness severity as well as specific risk by sex or family history status using up to 2762 cases and 3187 controls from the International Schizophrenia Consortium (ISC). Subjects with a family history of schizophrenia demonstrated a slightly lower average AAO that was not significant following multiple testing correction (p=.048), but no differences in illness severity were observed by family history status (p=.51). Consistent with prior reports, we observed earlier AAO (p=.005) and a more severe course of illness for men (p=.002). Family history positive analyses showed the greatest association with KIF5C (p=1.96×10(-8)), however, genetic risk burden overall does not differ by family history. Separate association analyses for males and females revealed no significant sex-specific associations. The top GWAS hit for AAO was near the olfactory receptor gene OR2K2 (p=1.52×10(-7)). Analyses of illness severity (episodic vs. continuous) implicated variation in ST18 (p=8.24×10(-7)). These results confirm recognized demographic relationships but do not support a simplified genetic architecture for schizophrenia subtypes based on these variables.

Socioeconomic status and social support following illicit drug use: causal pathways or common liability?

The negative social attributes associated with drug use and abuse/dependence may arise as a result of shared genetic or environmental factors rather than through causal pathways. To evaluate this possibility, structured interviews were conducted for 3969 male and female twins from the Mid-Atlantic Twin Registry and evaluations of their socioeconomic status (SES), social interactions, and use of drugs were obtained. Drug involvement was categorized as never used, tried, or met criteria for abuse or dependence. A co-twin control design was implemented using hierarchical linear modeling to assess whether twins who used drugs experienced lower SES and social support than non-using co-twins. Poorer social functioning in the drug-exposed twin is consistent with a causal relationship, while similar functioning in the drug exposed versus naive twins imply shared genetic or common environmental factors. Use of drugs was not significantly related to any SES measures. However, education and job status appear to share genetic influences with drug abuse/dependence. Lower income was not related to abuse/dependence of drugs. Negative interactions with friends and relatives share genetic factors with use of drugs, but the escalation from trying drugs to abusing them appears to generate discord between the abuser and friends and relatives in a causal fashion. These results indicate that presumptive causal influences of drug abuse/dependence on low SES may actually be mediated by shared genes. Drug use and social discord also appear to have shared genetic factors, but increased levels of drug involvement seem to causally influence social interactions.

Genome-wide association study identifies SESTD1 as a novel risk gene for lithium-responsive bipolar disorder

Lithium is the mainstay prophylactic treatment for bipolar disorder (BD), but treatment response varies considerably across individuals. Patients who respond well to lithium treatment might represent a relatively homogeneous subtype of this genetically and phenotypically diverse disorder. Here, we performed genome-wide association studies (GWAS) to identify (i) specific genetic variations influencing lithium response and (ii) genetic variants associated with risk for lithium-responsive BD. Patients with BD and controls were recruited from Sweden and the United Kingdom. GWAS were performed on 2698 patients with subjectively defined (self-reported) lithium response and 1176 patients with objectively defined (clinically documented) lithium response. We next conducted GWAS comparing lithium responders with healthy controls (1639 subjective responders and 8899 controls; 323 objective responders and 6684 controls). Meta-analyses of Swedish and UK results revealed no significant associations with lithium response within the bipolar subjects. However, when comparing lithium-responsive patients with controls, two imputed markers attained genome-wide significant associations, among which one was validated in confirmatory genotyping (rs116323614, P=2.74 × 10−8). It is an intronic single-nucleotide polymorphism (SNP) on chromosome 2q31.2 in the gene SEC14 and spectrin domains 1 (SESTD1), which encodes a protein involved in regulation of phospholipids. Phospholipids have been strongly implicated as lithium treatment targets. Furthermore, we estimated the proportion of variance for lithium-responsive BD explained by common variants (‘SNP heritability’) as 0.25 and 0.29 using two definitions of lithium response. Our results revealed a genetic variant in SESTD1 associated with risk for lithium-responsive BD, suggesting that the understanding of BD etiology could be furthered by focusing on this subtype of BD.

Suicidal Behavior During Lithium and Valproate Treatment: A Within-Individual 8-Year Prospective Study of 50,000 Patients With Bipolar Disorder

OBJECTIVE Conclusions regarding lithiums antisuicidal effect for bipolar disorder have been limited due to nonrepresentative subjects and potential confounding factors, including varying severity of illness. Findings regarding the effect of valproate, the most common alternative to lithium, are inconsistent for suicidal behavior. This study investigated the associations of these two drugs with the risk of suicide-related events, and possible differences between drugs, by using within-individual designs in a register-based longitudinal cohort. METHOD Through linkage of multiple Swedish national registers, 51,535 individuals with bipolar disorder were followed from 2005 to 2013 for treatment with lithium and valproate. Stratified Cox regression was used to estimate the hazard ratios of suicide-related events during treated periods compared with untreated periods. For significant associations between medication and suicide-related events, the population attributable fraction was estimated to assess the public health impact for patients with bipolar disorder. RESULTS During follow-up, 10,648 suicide-related events occurred. The incidence rate was significantly decreased by 14% during lithium treatment (hazard ratio 0.86, 95% confidence interval [CI] 0.78-0.95) but not during valproate treatment (hazard ratio 1.02, 95% CI 0.89-1.15). The difference in hazard ratios of suicide-related events between lithium and valproate was statistically significant. Estimates of the population attributable fraction suggested that 12% (95% CI 4%-20%) of suicide-related events could have been avoided if patients had taken lithium during the entire follow-up. CONCLUSIONS The results suggest that lithium should be considered for patients with bipolar disorder with suspected suicidal intentions, although risk for suicide is only one of the considerations when providing clinical care.

Acute intermittent porphyria: comorbidity and shared familial risks with schizophrenia and bipolar disorder in Sweden

Acute intermittent porphyria (AIP) has been associated with schizophrenia in some studies, but prior research is limited by the absence of comparison populations. Here, we linked Swedish registers to examine the risk of schizophrenia and bipolar disorder in 717 individuals diagnosed with AIP and their first-degree relatives, compared with matched individuals without AIP and their first-degree relatives. Individuals with AIP had a fourfold increased risk of schizophrenia or bipolar disorder. Similarly, relatives of individuals with AIP had double the risk of schizophrenia or bipolar disorder, suggesting that these associations may be as a result of common genetic influences.