Abstract Background Several theories have been proposed to conceptualize the pathological processes inherent to schizophrenia. The 'prostaglandin deficiency' hypothesis postulates that defective enzyme systems converting essential fatty acids to prostaglandins lead to diminished levels of prostaglandins, which in turn affect synaptic transmission. Methods Here we sought to determine the lipidomic profiles associated with schizophrenia in twin pairs discordant for schizophrenia as well as unaffected twin pairs. The study included serum samples from 19 twin pairs discordant for schizophrenia (mean age 51 ± 10 years; 7 monozygotic pairs; 13 female pairs) and 34 age and gender matched healthy twins as controls. Neurocognitive assessment data and gray matter density measurements taken from high-resolution magnetic resonance images were also obtained. A lipidomics platform using ultra performance liquid chromatography coupled to time-of-flight mass spectrometry was applied for the analysis of serum samples. Results In comparison to their healthy co-twins, the patients had elevated triglycerides and were more insulin resistant. They had diminished lysophosphatidylcholine levels, which associated with decreased cognitive speed. Conclusions Our findings may be of pathophysiological relevance since lysophosphatidylcholines, byproducts of phospholipase A2-catalyzed phospholipid hydrolysis, are preferred carriers of polyunsaturated fatty acids across the blood-brain barrier. Furthermore, diminishment of lysophosphatidylcholines suggests that subjects at risk of schizophrenia may be more susceptible to infections. Their association with cognitive speed supports the view that altered neurotransmission in schizophrenia may be in part mediated by reactive lipids such as prostaglandins.

Large-scale studies spanning diverse project sites, populations, languages, and measurements are increasingly important to relate psychological to biological variables. National and international consortia already are collecting and executing mega-analyses on aggregated data from individuals, with different measures on each person. In this research, we show that Asparouhov and Muthén's alignment method can be adapted to align data from disparate item sets and response formats. We argue that with these adaptations, the alignment method is well suited for combining data across multiple sites even when they use different measurement instruments. The approach is illustrated using data from the Whole Genome Sequencing in Psychiatric Disorders consortium and a real-data-based simulation is used to verify accurate parameter recovery. Factor alignment appears to increase precision of measurement and validity of scores with respect to external criteria. The resulting parameter estimates may further inform development of more effective and efficient methods to assess the same constructs in prospectively designed studies.

Several lines of evidence indicate the involvement of neuroinflammatory processes in the pathophysiology of schizophrenia (SCZ). Microglia are brain resident immune cells responding toward invading pathogens and injury-related products, and additionally, have a critical role in improving neurogenesis and synaptic functions. Aberrant activation of microglia in SCZ is one of the leading hypotheses for disease pathogenesis, but due to the lack of proper human cell models, the role of microglia in SCZ is not well studied. We used monozygotic twins discordant for SCZ and healthy individuals to generate human induced pluripotent stem cell-derived microglia to assess the transcriptional and functional differences in microglia between healthy controls, affected twins and unaffected twins. The microglia from affected twins had increased expression of several common inflammation-related genes compared to healthy individuals. Microglia from affected twins had also reduced response to interleukin 1 beta (IL1β) treatment, but no significant differences in migration or phagocytotic activity. Ingenuity Pathway Analysis (IPA) showed abnormalities related to extracellular matrix signaling. RNA sequencing predicted downregulation of extracellular matrix structure constituent Gene Ontology (GO) terms and hepatic fibrosis pathway activation that were shared by microglia of both affected and unaffected twins, but the upregulation of major histocompatibility complex (MHC) class II receptors was observed only in affected twin microglia. Also, the microglia of affected twins had heterogeneous response to clozapine, minocycline, and sulforaphane treatments. Overall, despite the increased expression of inflammatory genes, we observed no clear functional signs of hyperactivation in microglia from patients with SCZ. We conclude that microglia of the patients with SCZ have gene expression aberrations related to inflammation response and extracellular matrix without contributing to increased microglial activation.

Implicating particular genes in the generation of complex brain and behavior phenotypes requires multiple lines of evidence. The rarity of most high-impact genetic variants typically precludes the possibility of accruing statistical evidence that they are associated with a given trait. We found that the enrichment of a rare chromosome 22q11.22 deletion in a recently expanded Northern Finnish sub-isolate enabled the detection of association between TOP3B and both schizophrenia and cognitive impairment. Biochemical analysis of TOP3β revealed that this topoisomerase was a component of cytosolic messenger ribonucleoproteins (mRNPs) and was catalytically active on RNA. The recruitment of TOP3β to mRNPs was independent of RNA cis-elements and was coupled to the co-recruitment of FMRP, the disease gene product in fragile X mental retardation syndrome. Our results indicate a previously unknown role for TOP3β in mRNA metabolism and suggest that it is involved in neurodevelopmental disorders.

By analyzing the whole-exome sequences of 4,264 schizophrenia cases, 9,343 controls and 1,077 trios, we identified a genome-wide significant association between rare loss-of-function (LoF) variants in SETD1A and risk for schizophrenia (P = 3.3 × 10(-9)). We found only two heterozygous LoF variants in 45,376 exomes from individuals without a neuropsychiatric diagnosis, indicating that SETD1A is substantially depleted of LoF variants in the general population. Seven of the ten individuals with schizophrenia carrying SETD1A LoF variants also had learning difficulties. We further identified four SETD1A LoF carriers among 4,281 children with severe developmental disorders and two more carriers in an independent sample of 5,720 Finnish exomes, both with notable neuropsychiatric phenotypes. Together, our observations indicate that LoF variants in SETD1A cause a range of neurodevelopmental disorders, including schizophrenia. Combining these data with previous common variant evidence, we suggest that epigenetic dysregulation, specifically in the histone H3K4 methylation pathway, is an important mechanism in the pathogenesis of 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 disorders¹, epilepsy and severe neurodevelopmental disorders², 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 risk³, 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.

Copy number variants (CNVs) have been strongly implicated in the genetic etiology of schizophrenia (SCZ). However, genome-wide investigation of the contribution of CNV to risk has been hampered by limited sample sizes. We sought to address this obstacle by applying a centralized analysis pipeline to a SCZ cohort of 21,094 cases and 20,227 controls. A global enrichment of CNV burden was observed in cases (odds ratio (OR) = 1.11, P = 5.7 × 10^-15), which persisted after excluding loci implicated in previous studies (OR = 1.07, P = 1.7 × 10^-6). CNV burden was enriched for genes associated with synaptic function (OR = 1.68, P = 2.8 × 10^-11) and neurobehavioral phenotypes in mouse (OR = 1.18, P = 7.3 × 10^-5). Genome-wide significant evidence was obtained for eight loci, including 1q21.1, 2p16.3 (NRXN1), 3q29, 7q11.2, 15q13.3, distal 16p11.2, proximal 16p11.2 and 22q11.2. Suggestive support was found for eight additional candidate susceptibility and protective loci, which consisted predominantly of CNVs mediated by nonallelic homologous recombination.

Genome-wide association studies (GWASs) are highly effective at identifying common risk variants for schizophrenia. Rare risk variants are also important contributors to schizophrenia etiology but, with the exception of large copy number variants, are difficult to detect with GWAS. Exome and genome sequencing, which have accelerated the study of rare variants, are expensive so alternative methods are needed to aid detection of rare variants. Here we re-analyze an Irish schizophrenia GWAS dataset (n = 3,473) by performing identity-by-descent (IBD) mapping followed by exome sequencing of individuals identified as sharing risk haplotypes to search for rare risk variants in coding regions. We identified 45 rare haplotypes (>1 cM) that were significantly more common in cases than controls. By exome sequencing 105 haplotype carriers, we investigated these haplotypes for functional coding variants that could be tested for association in independent GWAS samples. We identified one rare missense variant in PCNT but did not find statistical support for an association with schizophrenia in a replication analysis. However, IBD mapping can prioritize both individual samples and genomic regions for follow-up analysis but genome rather than exome sequencing may be more effective at detecting risk variants on rare haplotypes.

Regulatory and coding variants are known to be enriched with associations identified by genome-wide association studies (GWASs) of complex disease, but their contributions to trait heritability are currently unknown. We applied variance-component methods to imputed genotype data for 11 common diseases to partition the heritability explained by genotyped SNPs (hg(2)) across functional categories (while accounting for shared variance due to linkage disequilibrium). Extensive simulations showed that in contrast to current estimates from GWAS summary statistics, the variance-component approach partitions heritability accurately under a wide range of complex-disease architectures. Across the 11 diseases DNaseI hypersensitivity sites (DHSs) from 217 cell types spanned 16% of imputed SNPs (and 24% of genotyped SNPs) but explained an average of 79% (SE = 8%) of hg(2) from imputed SNPs (5.1× enrichment; p = 3.7 × 10(-17)) and 38% (SE = 4%) of hg(2) from genotyped SNPs (1.6× enrichment, p = 1.0 × 10(-4)). Further enrichment was observed at enhancer DHSs and cell-type-specific DHSs. In contrast, coding variants, which span 1% of the genome, explained <10% of hg(2) despite having the highest enrichment. We replicated these findings but found no significant contribution from rare coding variants in independent schizophrenia cohorts genotyped on GWAS and exome chips. Our results highlight the value of analyzing components of heritability to unravel the functional architecture of common disease.

Schizophrenia is a highly heritable disorder. Genetic risk is conferred by a large number of alleles, including common alleles of small effect that might be detected by genome-wide association studies. Here we report a multi-stage schizophrenia genome-wide association study of up to 36,989 cases and 113,075 controls. We identify 128 independent associations spanning 108 conservatively defined loci that meet genome-wide significance, 83 of which have not been previously reported. Associations were enriched among genes expressed in brain, providing biological plausibility for the findings. Many findings have the potential to provide entirely new insights into aetiology, but associations at DRD2 and several genes involved in glutamatergic neurotransmission highlight molecules of known and potential therapeutic relevance to schizophrenia, and are consistent with leading pathophysiological hypotheses. Independent of genes expressed in brain, associations were enriched among genes expressed in tissues that have important roles in immunity, providing support for the speculated link between the immune system and schizophrenia.