Increased exonic de novo mutation rate in individuals with schizophrenia

De Novo Variants Found in Three Distinct Schizophrenia Populations Hit a Common Core Gene Network Related to Microtubule and Actin Cytoskeleton Gene Ontology Classes

Schizophrenia (SZ) is a heterogeneous and debilitating psychiatric disorder with a strong genetic component. To elucidate functional networks perturbed in schizophrenia, we analysed a large dataset of whole-genome studies that identified SNVs, CNVs, and a multi-stage schizophrenia genome-wide association study. Our analysis identified three subclusters that are interrelated and with small overlaps: GO:0007017~Microtubule-Based Process, GO:00015629~Actin Cytoskeleton, and GO:0007268~SynapticTransmission. We next analysed three distinct trio cohorts of 75 SZ Algerian, 45 SZ French, and 61 SZ Japanese patients. We performed Illumina HiSeq whole-exome sequencing and identified de novo mutations using a Bayesian approach. We validated 88 de novo mutations by Sanger sequencing: 35 in French, 21 in Algerian, and 32 in Japanese SZ patients. These 88 de novo mutations exhibited an enrichment in genes encoding proteins related to GO:0051015~actin filament binding (p = 0.0011) using David, and enrichments in GO: 0003774~transport (p = 0.019) and GO:0003729~mRNA binding (p = 0.010) using Amigo. One of these de novo variant was found in CORO1C coding sequence. We studied Coro1c haploinsufficiency in a Coro1c+/- mouse and found defects in the corpus callosum. These results could motivate future studies of the mechanisms surrounding genes encoding proteins involved in transport and the cytoskeleton, with the goal of developing therapeutic intervention strategies for a subset of SZ cases.

A Kpna1-deficient psychotropic drug-induced schizophrenia model mouse for studying gene-environment interactions

KPNA1 is a mediator of nucleocytoplasmic transport that is abundantly expressed in the mammalian brain and regulates neuronal differentiation and synaptic function. De novo mutations in Kpna1 have been identified using genome-wide association studies in humans with schizophrenia; however, it remains unclear how KPNA1 contributes to schizophrenia pathogenesis. Recent studies have suggested a complex combination of genetic and environmental factors that are closely related to psychiatric disorders. Here, we found that subchronic administration of phencyclidine, a psychotropic drug, induced vulnerability and behavioral abnormalities consistent with the symptoms of schizophrenia in Kpna1-deficient mice. Microarray assessment revealed that the expression levels of dopamine d1/d2 receptors, an RNA editing enzyme, and a cytoplasmic dynein component were significantly altered in the nucleus accumbens brain region in a gene-environment (G × E) interaction-dependent manner. Our findings demonstrate that Kpna1-deficient mice may be useful as a G × E interaction mouse model for psychiatric disorders and for further investigation into the pathogenesis of such diseases and disorders.

The Association of Hippocampal Long-Term Potentiation-Induced Gene Expression with Genetic Risk for Psychosis

Genomic studies focusing on the contribution of common and rare genetic variants of schizophrenia and bipolar disorder support the view that substantial risk is conferred through molecular pathways involved in synaptic plasticity in the neurons of cortical and subcortical brain regions, including the hippocampus. Synaptic long-term potentiation (LTP) is central to associative learning and memory and depends on a pattern of gene expression in response to neuronal stimulation. Genes related to the induction of LTP have been associated with psychiatric genetic risk, but the specific cell types and timepoints responsible for the association are unknown. Using published genomic and transcriptomic datasets, we studied the relationship between temporally defined gene expression in hippocampal pyramidal neurons following LTP and enrichment for common genetic risk for schizophrenia and bipolar disorder, and for copy number variants (CNVs) and de novo coding variants associated with schizophrenia. We observed that upregulated genes in hippocampal pyramidal neurons at 60 and 120 min following LTP induction were enriched for common variant association with schizophrenia and bipolar disorder subtype I. At 60 min, LTP-induced genes were enriched in duplications from patients with schizophrenia, but this association was not specific to pyramidal neurons, perhaps reflecting the combined effects of CNVs in excitatory and inhibitory neuron subtypes. Gene expression following LTP was not related to enrichment for de novo coding variants from schizophrenia cases. Our findings refine our understanding of the role LTP-related gene sets play in conferring risk to conditions causing psychosis and provide a focus for future studies looking to dissect the molecular mechanisms associated with this risk.

Genes positively regulated by Mef2c in cortical neurons are enriched for common genetic variation associated with IQ and educational attainment

The myocyte enhancer factor 2 C (MEF2C) gene encodes a transcription factor important for neurogenesis and synapse development and contains common variants associated with intelligence (IQ) and educational attainment (EA). Here, we took gene expression data from the mouse cortex of a Mef2c mouse model with a heterozygous DNA binding-deficient mutation of Mef2c (Mef2c-het) and combined these data with MEF2C ChIP-seq data from cortical neurons and single-cell data from the mouse brain. This enabled us to create a set of genes that were differentially regulated in Mef2c-het mice, represented direct target genes of MEF2C and had elevated in expression in cortical neurons. We found this gene-set to be enriched for genes containing common genetic variation associated with IQ and EA. Genes within this gene-set that were down-regulated, i.e. have reduced expression in Mef2c-het mice versus controls, were specifically significantly enriched for both EA and IQ associated genes. These down-regulated genes were enriched for functionality in the adenylyl cyclase signalling system, which is known to positively regulate synaptic transmission and has been linked to learning and memory. Within the adenylyl cyclase signalling system, three genes regulated by MEF2C, CRHR1, RGS6, and GABRG3, are associated at genome-wide significant levels with IQ and/or EA. Our results indicate that genetic variation in MEF2C and its direct target genes within cortical neurons contribute to variance in cognition within the general population, and the molecular mechanisms involved include the adenylyl cyclase signalling system's role in synaptic function.

Brain microvascular endothelial cells and blood-brain barrier dysfunction in psychotic disorders

Although there is convergent evidence for blood-brain barrier (BBB) dysfunction and peripheral inflammation in schizophrenia (SZ) and bipolar disorder (BD), it is unknown whether BBB deficits are intrinsic to brain microvascular endothelial cells (BMECs) or arise via effects of peripheral inflammatory cytokines. We examined BMEC function using stem cell-based models to identify cellular and molecular deficits associated with BBB dysfunction in SZ and BD. Induced pluripotent stem cells (iPSCs) from 4 SZ, 4 psychotic BD and 4 healthy control (HC) subjects were differentiated into BMEC-"like" cells. Gene expression and protein levels of tight junction proteins were assessed. Transendothelial electrical resistance (TEER) and permeability were assayed to evaluate BBB function. Cytokine levels were measured from conditioned media. BMECs derived from human iPSCs in SZ and BD did not show differences in BBB integrity or permeability compared to HC BMECs. Outlier analysis using TEER revealed a BBB-deficit (n = 3) and non-deficit (n = 5) group in SZ and BD lines. Stratification based on BBB function in SZ and BD patients identified a BBB-deficit subtype with reduced barrier function, tendency for increased permeability to smaller molecules, and decreased claudin-5 (CLDN5) levels. BMECs from the BBB-deficit group show increased matrix metallopeptidase 1 (MMP1) activity, which correlated with reduced CLDN5 and worse BBB function, and was improved by tumor necrosis factor α (TNFα) and MMP1 inhibition. These results show potential deficits in BMEC-like cells in psychotic disorders that result in BBB disruption and further identify TNFα and MMP1 as promising targets for ameliorating BBB deficits.

Experimental estimates of germline mutation rate in eukaryotes: a phylogenetic meta-analysis

Mutation is the ultimate source of all genetic variation, and over the last 10 years the ready availability of whole-genome sequencing has permitted direct estimation of mutation rate for many non-model species across the tree of life. In this meta-analysis, we make a comprehensive search of the literature for mutation rate estimates in eukaryotes, identifying 140 mutation accumulation (MA) and parent-offspring (PO) sequencing studies covering 134 species. Based on these data, we revisit differences in the single-nucleotide mutation (SNM) rate between different phylogenetic lineages and update the known relationships between mutation rate and generation time, genome size, and nucleotide diversity-while accounting for phylogenetic nonindependence. We do not find a significant difference between MA and PO in estimated mutation rates, but we confirm that mammal and plant lineages have higher mutation rates than arthropods and that unicellular eukaryotes have the lowest mutation rates. We find that mutation rates are higher in species with longer generation times and larger genome sizes, even when accounting for phylogenetic relationships. Moreover, although nucleotide diversity is positively correlated with mutation rate, the gradient of the relationship is significantly less than one (on a logarithmic scale), consistent with higher mutation rates in populations with smaller effective size. For the 29 species for which data are available, we find that indel mutation rates are positively correlated with nucleotide mutation rates and that short deletions are generally more common than short insertions. Nevertheless, despite recent progress, no estimates of either SNM or indel mutation rates are available for the majority of deeply branching eukaryotic lineages-or even for most animal phyla. Even among charismatic megafauna, experimental mutation rate estimates remain unknown for amphibia and scarce for reptiles and fish.

The complete chloroplast genome of Ulmus mianzhuensis with insights into structural variations, adaptive evolution, and phylogenetic relationships of Ulmus (Ulmaceae)

BACKGROUND

Ulmus mianzhuensis is an endemic tree species in China with high ornamental and economic value. Currently, little is known regarding its genomic architecture, phylogenetic position, or adaptive evolution. Here, we sequenced the complete chloroplast genome (cp genome) of U. mianzhuensis and further compared the variations in gene organization and structure within Ulmus species to define their genomic evolution, then reconstructed the phylogenomic relationship of 31 related Ulmus species to explore the systematic position of U. mianzhuensis and the utility of cp genome for resolving phylogenetics among Ulmus species.

RESULTS

Our results revealed that all the Ulmus species exhibited a typical quadripartite structure, with a large single copy (LSC) region of 87,170 - 88,408 bp, a small single copy (SSC) region of 18,650 - 19,038 bp and an inverted repeat (IR) region of 26,288 - 26,546 bp. Within Ulmus species, gene structure and content of cp genomes were highly conserved, although slight variations were found in the boundary of SC/IR regions. Moreover, genome-wide sliding window analysis uncovered the variability of ndhC-trnV-UAC, ndhF-rpl32, and psbI-trnS-GCU were higher among 31 Ulmus that may be useful for the population genetics and potential DNA barcodes. Two genes (rps15 and atpF) were further detected under a positive selection of Ulmus species. Comparative phylogenetic analysis based on the cp genome and protein-coding genes revealed consistent topology that U. mianzhuensis is a sister group to U. parvifolia (sect. Microptelea) with a relatively low-level nucleotide variation of the cp genome. Additionally, our analyses also found that the traditional taxonomic system of five sections in Ulmus is not supported by the current phylogenomic topology with a nested evolutionary relationship between sections.

CONCLUSIONS

Features of the cp genome length, GC content, organization, and gene order were highly conserved within Ulmus. Furthermore, molecular evidence from the low variation of the cp genome suggested that U. mianzhuensis should be merged into U. parvifolia and regarded as a subspecies of U. parvifolia. Overall, we demonstrated that the cp genome provides valuable information for understanding the genetic variation and phylogenetic relationship in Ulmus.

The genetic basis of onset age in schizophrenia: evidence and models

Schizophrenia is a heritable neurocognitive disorder affecting about 1% of the population, and usually has an onset age at around 21-25 in males and 25-30 in females. Recent advances in genetics have helped to identify many common and rare variants for the liability to schizophrenia. Earlier evidence appeared to suggest that younger onset age is associated with higher genetic liability to schizophrenia. Clinical longitudinal research also found that early and very-early onset schizophrenia are associated with poor clinical, neurocognitive, and functional profiles. A recent study reported a heritability of 0.33 for schizophrenia onset age, but the genetic basis of this trait in schizophrenia remains elusive. In the pre-Genome-Wide Association Study (GWAS) era, genetic loci found to be associated with onset age were seldom replicated. In the post-Genome-Wide Association Study era, new conceptual frameworks are needed to clarify the role of onset age in genetic research in schizophrenia, and to identify its genetic basis. In this review, we first discussed the potential of onset age as a characterizing/subtyping feature for psychosis, and as an important phenotypic dimension of schizophrenia. Second, we reviewed the methods, samples, findings and limitations of previous genetic research on onset age in schizophrenia. Third, we discussed a potential conceptual framework for studying the genetic basis of onset age, as well as the concepts of susceptibility, modifier, and "mixed" genes. Fourth, we discussed the limitations of this review. Lastly, we discussed the potential clinical implications for genetic research of onset age of schizophrenia, and how future research can unveil the potential mechanisms for this trait.

De novo mutations disturb early brain development more frequently than common variants in schizophrenia

Investigating functional, temporal, and cell-type expression features of mutations is important for understanding a complex disease. Here, we collected and analyzed common variants and de novo mutations (DNMs) in schizophrenia (SCZ). We collected 2,636 missense and loss-of-function (LoF) DNMs in 2,263 genes across 3,477 SCZ patients (SCZ-DNMs). We curated three gene lists: (a) SCZ-neuroGenes (159 genes), which are intolerant to LoF and missense DNMs and are neurologically important, (b) SCZ-moduleGenes (52 genes), which were derived from network analyses of SCZ-DNMs, and (c) SCZ-commonGenes (120 genes) from a recent GWAS as reference. To compare temporal gene expression, we used the BrainSpan dataset. We defined a fetal effect score (FES) to quantify the involvement of each gene in prenatal brain development. We further employed the specificity indexes (SIs) to evaluate cell-type expression specificity from single-cell expression data in cerebral cortices of humans and mice. Compared with SCZ-commonGenes, SCZ-neuroGenes and SCZ-moduleGenes were highly expressed in the prenatal stage, had higher FESs, and had higher SIs in fetal replicating cells and undifferentiated cell types. Our results suggested that gene expression patterns in specific cell types in early fetal stages might have impacts on the risk of SCZ during adulthood.

The genetic architecture of schizophrenia: review of large-scale genetic studies

Schizophrenia is a complex and often chronic psychiatric disorder with high heritability. Diagnosis of schizophrenia is still made clinically based on psychiatric symptoms; no diagnostic tests or biomarkers are available. Pathophysiology-based diagnostic scheme and treatments are also not available. Elucidation of the pathogenesis is needed for development of pathology-based diagnostics and treatments. In the past few decades, genetic research has made substantial advances in our understanding of the genetic architecture of schizophrenia. Rare copy number variations (CNVs) and rare single-nucleotide variants (SNVs) detected by whole-genome CNV analysis and whole-genome/-exome sequencing analysis have provided the great advances. Common single-nucleotide polymorphisms (SNPs) detected by large-scale genome-wide association studies have also provided important information. Large-scale genetic studies have been revealed that both rare and common genetic variants play crucial roles in this disorder. In this review, we focused on CNVs, SNVs, and SNPs, and discuss the latest research findings on the pathogenesis of schizophrenia based on these genetic variants. Rare variants with large effect sizes can provide mechanistic hypotheses. CRISPR-based genetics approaches and induced pluripotent stem cell technology can facilitate the functional analysis of these variants detected in patients with schizophrenia. Recent advances in long-read sequence technology are expected to detect variants that cannot be detected by short-read sequence technology. Various studies that bring together data from common variant and transcriptomic datasets provide biological insight. These new approaches will provide additional insight into the pathophysiology of schizophrenia and facilitate the development of pathology-based therapeutics.

Delineation of a KDM2B-related neurodevelopmental disorder and its associated DNA methylation signature

PURPOSE

Pathogenic variants in genes involved in the epigenetic machinery are an emerging cause of neurodevelopment disorders (NDDs). Lysine-demethylase 2B (KDM2B) encodes an epigenetic regulator and mouse models suggest an important role during development. We set out to determine whether KDM2B variants are associated with NDD.

METHODS

Through international collaborations, we collected data on individuals with heterozygous KDM2B variants. We applied methylation arrays on peripheral blood DNA samples to determine a KDM2B associated epigenetic signature.

RESULTS

We recruited a total of 27 individuals with heterozygous variants in KDM2B. We present evidence, including a shared epigenetic signature, to support a pathogenic classification of 15 KDM2B variants and identify the CxxC domain as a mutational hotspot. Both loss-of-function and CxxC-domain missense variants present with a specific subepisignature. Moreover, the KDM2B episignature was identified in the context of a dual molecular diagnosis in multiple individuals. Our efforts resulted in a cohort of 21 individuals with heterozygous (likely) pathogenic variants. Individuals in this cohort present with developmental delay and/or intellectual disability; autism; attention deficit disorder/attention deficit hyperactivity disorder; congenital organ anomalies mainly of the heart, eyes, and urogenital system; and subtle facial dysmorphism.

CONCLUSION

Pathogenic heterozygous variants in KDM2B are associated with NDD and a specific epigenetic signature detectable in peripheral blood.

A rigorous in silico genomic interrogation at 1p13.3 reveals 16 autosomal dominant candidate genes in syndromic neurodevelopmental disorders

Genome-wide chromosomal microarray is extensively used to detect copy number variations (CNVs), which can diagnose microdeletion and microduplication syndromes. These small unbalanced chromosomal structural rearrangements ranging from 1 kb to 10 Mb comprise up to 15% of human mutations leading to monogenic or contiguous genomic disorders. Albeit rare, CNVs at 1p13.3 cause a variety of neurodevelopmental disorders (NDDs) including development delay (DD), intellectual disability (ID), autism, epilepsy, and craniofacial anomalies (CFA). Most of the 1p13.3 CNV cases reported in the pre-microarray era encompassed a large number of genes and lacked the demarcating genomic coordinates, hampering the discovery of positional candidate genes within the boundaries. In this study, we present four subjects with 1p13.3 microdeletions displaying DD, ID, autism, epilepsy, and CFA. In silico comparative genomic mapping with three previously reported subjects with CNVs and 22 unreported DECIPHER CNV cases has resulted in the identification of four different sub-genomic loci harboring five positional candidate genes for DD, ID, and CFA at 1p13.3. Most of these genes have pathogenic variants reported, and their interacting genes are involved in NDDs. RT-qPCR in various human tissues revealed a high expression pattern in the brain and fetal brain, supporting their functional roles in NDDs. Interrogation of variant databases and interacting protein partners led to the identification of another set of 11 potential candidate genes, which might have been dysregulated by the position effect of these CNVs at 1p13.3. Our studies define 1p13.3 as a genomic region harboring 16 NDD candidate genes and underscore the critical roles of small CNVs in in silico comparative genomic mapping for disease gene discovery. Our candidate genes will help accelerate the isolation of pathogenic heterozygous variants from exome/genome sequencing (ES/GS) databases.

Evolution and dysfunction of human cognitive and social traits: A transcriptional regulation perspective

Evolutionary changes in brain and craniofacial development have endowed humans with unique cognitive and social skills, but also predisposed us to debilitating disorders in which these traits are disrupted. What are the developmental genetic underpinnings that connect the adaptive evolution of our cognition and sociality with the persistence of mental disorders with severe negative fitness effects? We argue that loss of function of genes involved in transcriptional regulation represents a crucial link between the evolution and dysfunction of human cognitive and social traits. The argument is based on the haploinsufficiency of many transcriptional regulator genes, which makes them particularly sensitive to loss-of-function mutations. We discuss how human brain and craniofacial traits evolved through partial loss of function (i.e. reduced expression) of these genes, a perspective compatible with the idea of human self-domestication. Moreover, we explain why selection against loss-of-function variants supports the view that mutation-selection-drift, rather than balancing selection, underlies the persistence of psychiatric disorders. Finally, we discuss testable predictions.

De novo mutations, genetic mosaicism and human disease

Mosaicism—the existence of genetically distinct populations of cells in a particular organism—is an important cause of genetic disease. Mosaicism can appear as de novo DNA mutations, epigenetic alterations of DNA, and chromosomal abnormalities. Neurodevelopmental or neuropsychiatric diseases, including autism—often arise by de novo mutations that usually not present in either of the parents. De novo mutations might occur as early as in the parental germline, during embryonic, fetal development, and/or post-natally, through ageing and life. Mutation timing could lead to mutation burden of less than heterozygosity to approaching homozygosity. Developmental timing of somatic mutation attainment will affect the mutation load and distribution throughout the body. In this review, we discuss the timing of de novo mutations, spanning from mutations in the germ lineage (all ages), to post-zygotic, embryonic, fetal, and post-natal events, through aging to death. These factors can determine the tissue specific distribution and load of de novo mutations, which can affect disease. The disease threshold burden of somatic de novo mutations of a particular gene in any tissue will be important to define.

Developmental disruption to the cortical transcriptome and synaptosome in a model of SETD1A loss-of-function

Large-scale genomic studies of schizophrenia implicate genes involved in the epigenetic regulation of transcription by histone methylation and genes encoding components of the synapse. However, the interactions between these pathways in conferring risk to psychiatric illness are unknown. Loss-of-function (LoF) mutations in the gene encoding histone methyltransferase, SETD1A, confer substantial risk to schizophrenia. Among several roles, SETD1A is thought to be involved in the development and function of neuronal circuits. Here, we employed a multi-omics approach to study the effects of heterozygous Setd1a LoF on gene expression and synaptic composition in mouse cortex across five developmental timepoints from embryonic day 14 to postnatal day 70. Using RNA sequencing, we observed that Setd1a LoF resulted in the consistent downregulation of genes enriched for mitochondrial pathways. This effect extended to the synaptosome, in which we found age-specific disruption to both mitochondrial and synaptic proteins. Using large-scale patient genomics data, we observed no enrichment for genetic association with schizophrenia within differentially expressed transcripts or proteins, suggesting they derive from a distinct mechanism of risk from that implicated by genomic studies. This study highlights biological pathways through which SETD1A LOF may confer risk to schizophrenia. Further work is required to determine whether the effects observed in this model reflect human pathology.

Exploration of Mutated Genes and Prediction of Potential Biomarkers for Childhood-Onset Schizophrenia Using an Integrated Bioinformatic Analysis

Childhood-onset schizophrenia (COS) is an unusual severe neurodevelopmental disorder of unknown etiology. In this study, we aimed to survey the missense variants in new cases of COS and also identify possible pathology biomarkers for COS. We found one list of mutated genes such as TTN, MUC12, and MUC2, which are the candidates to be involved in the etiology of COS. Next, we used WGSNA to predict COS disease-related genes and identified differential DNA methylation among COS disease groups, COS dangerous groups, and normal groups and found eight methylation sites that can be used as the diagnostic biomarkers. A total of six key genes are obtained through the intersection analysis between weighted correlation network analysis (WGCNA) mode, methylation-related genes, and differentially expressed genes (DGenes). These genes may play important roles in the progression of COS and serve as the potential biomarkers for future diagnosis. Our results might help to design the molecule or gene-targeted drugs for COS.

In vitro fluorescence assay to measure GDP/GTP exchange of guanine nucleotide exchange factors of Rho family GTPases

Guanine nucleotide exchange factors (GEFs) are enzymes that promote the activation of GTPases through GTP loading. Whole exome sequencing has identified rare variants in GEFs that are associated with disease, demonstrating that GEFs play critical roles in human development. However, the consequences of these rare variants can only be understood through measuring their effects on cellular activity. Here, we provide a detailed, user-friendly protocol for purification and fluorescence-based analysis of the two GEF domains within the protein, Trio. This analysis offers a straight-forward, quantitative tool to test the activity of GEF domains on their respective GTPases, as well as utilize high-throughput screening to identify regulators and inhibitors. This protocol can be adapted for characterization of other Rho family GEFs. Such analyses are crucial for the complete understanding of the roles of GEF genetic variants in human development and disease.

Effects of Importin α1/KPNA1 deletion and adolescent social isolation stress on psychiatric disorder-associated behaviors in mice

Importin α1/KPNA1 is a member of the Importin α family widely present in the mammalian brain and has been characterized as a regulator of neuronal differentiation, synaptic functionality, and anxiety-like behavior. In humans, a de novo mutation of the KPNA1 (human Importin α5) gene has been linked with schizophrenia; however, the precise roles of KPNA1 in disorder-related behaviors are still unknown. Moreover, as recent studies have highlighted the importance of gene-environment interactions in the development of psychiatric disorders, we investigated the effects of Kpna1 deletion and social isolation stress, a paradigm that models social stress factors found in human patients, on psychiatric disorder-related behaviors in mice. Through assessment in a behavioral battery, we found that Kpna1 knockout resulted in the following behavioral phenotype: (1) decreased anxiety-like behavior in an elevated plus maze test, (2) short term memory deficits in novel object recognition test (3) impaired sensorimotor gating in a prepulse inhibition test. Importantly, exposure to social isolation stress resulted in additional behavioral abnormalities where isolated Kpna1 knockout mice exhibited: (1) impaired aversive learning and/or memory in the inhibitory avoidance test, as well as (2) increased depression-like behavior in the forced swim test. Furthermore, we investigated whether mice showed alterations in plasma levels of stress-associated signal molecules (corticosterone, cytokines, hormones, receptors), and found that Kpna1 knockout significantly altered levels of corticosterone and LIX (CXCL5). Moreover, significant decreases in the level of prolactin were found in all groups except for group-housed wild type mice. Our findings demonstrate that Kpna1 deletion can trigger widespread behavioral abnormalities associated with psychiatric disorders, some of which were further exacerbated by exposure to adolescent social isolation. The use of Kpna1 knockout mice as a model for psychiatric disorders may show promise for further investigation of gene-environment interactions involved in the pathogenesis of psychiatric disorders.

Schizophrenia, autism spectrum disorders and developmental disorders share specific disruptive coding mutations

People with schizophrenia are enriched for rare coding variants in genes associated with neurodevelopmental disorders, particularly autism spectrum disorders and intellectual disability. However, it is unclear if the same changes to gene function that increase risk to neurodevelopmental disorders also do so for schizophrenia. Using data from 3444 schizophrenia trios and 37,488 neurodevelopmental disorder trios, we show that within shared risk genes, de novo variants in schizophrenia and neurodevelopmental disorders are generally of the same functional category, and that specific de novo variants observed in neurodevelopmental disorders are enriched in schizophrenia (P = 5.0 × 10-6). The latter includes variants known to be pathogenic for syndromic disorders, suggesting that schizophrenia be included as a characteristic of those syndromes. Our findings imply that, in part, neurodevelopmental disorders and schizophrenia have shared molecular aetiology, and therefore likely overlapping pathophysiology, and support the hypothesis that at least some forms of schizophrenia lie on a continuum of neurodevelopmental disorders.

Rare germline variants in individuals diagnosed with schizophrenia within multiplex families

An extensive catalog of common and rare genetic variants contributes to overall risk for schizophrenia and related disorders. As a complement to population genetics efforts, here we present whole genome sequences of multiple affected probands within individual families to search for possible high penetrance driver variants. From a total of 15 families diagnostically evaluated by a single research psychiatrist, we performed whole genome sequencing of a total of 61 affected individuals, called SNPs, indels, and copy number variants, and compared to reference genomes. In fourteen out of fifteen families, the schizophrenia polygenic risk score for each proband was within the control range defined by the Thousand Genomes cohort. In six families, each affected member carried a very rare or private, predicted-damaging, variant in at least one gene. Among these genes, variants in LRP1 and TENM2 suggest these are candidate disease-related genes when taken into context with existing population genetic studies and biological information. Results add to the number of pedigree sequences reported, suggest pathways for the investigation of biological mechanisms, and are consistent with the overall accumulating evidence that very rare damaging variants contribute to the heritability of schizophrenia.

Transcriptome Analysis of Post-Mortem Brain Tissue Reveals Up-Regulation of the Complement Cascade in a Subgroup of Schizophrenia Patients

Schizophrenia is a genetically complex neuropsychiatric disorder with largely unresolved mechanisms of pathology. Identification of genes and pathways associated with schizophrenia is important for understanding the development, progression and treatment of schizophrenia. In this study, pathways associated with schizophrenia were explored at the level of gene expression. The study included post-mortem brain tissue samples from 68 schizophrenia patients and 44 age and sex-matched control subjects. Whole transcriptome poly-A selected paired-end RNA sequencing was performed on tissue from the prefrontal cortex and orbitofrontal cortex. RNA expression differences were detected between case and control individuals, focusing both on single genes and pathways. The results were validated with RT-qPCR. Significant differential expression between patient and controls groups was found for 71 genes. Gene ontology analysis of differentially expressed genes revealed an up-regulation of multiple genes in immune response among the patients (corrected p-value = 0.004). Several genes in the category belong to the complement system, including C1R, C1S, C7, FCN3, SERPING1, C4A and CFI. The increased complement expression is primarily driven by a subgroup of patients with increased expression of immune/inflammatory response genes, pointing to important differences in disease etiology within the patient group. Weighted gene co-expression network analysis highlighted networks associated with both synaptic transmission and activation of the immune response. Our results demonstrate the importance of immune-related pathways in schizophrenia and provide evidence for elevated expression of the complement cascade as an important pathway in schizophrenia pathology.

DNA methylation differences in cortical grey and white matter in schizophrenia

Aim: Identify grey- and white-matter-specific DNA-methylation differences between schizophrenia (SCZ) patients and controls in postmortem brain cortical tissue. Materials & methods: Grey and white matter were separated from postmortem brain tissue of the superior temporal and medial frontal gyrus from SCZ (n = 10) and control (n = 11) cases. Genome-wide DNA-methylation analysis was performed using the Infinium EPIC Methylation Array (Illumina, CA, USA). Results: Four differentially methylated regions associated with SCZ status and tissue type (grey vs white matter) were identified within or near KLF9, SFXN1, SPRED2 and ALS2CL genes. Gene-expression analysis showed differential expression of KLF9 and SFXN1 in SCZ. Conclusion: Our data show distinct differences in DNA methylation between grey and white matter that are unique to SCZ, providing new leads to unravel the pathogenesis of SCZ.

Low-Level Brain Somatic Mutations Are Implicated in Schizophrenia

BACKGROUND

Somatic mutations arising from the brain have recently emerged as significant contributors to neurodevelopmental disorders, including childhood intractable epilepsy and cortical malformations. However, whether brain somatic mutations are implicated in schizophrenia (SCZ) is not well established.

METHODS

We performed deep whole exome sequencing (average read depth > 550×) of matched dorsolateral prefrontal cortex and peripheral tissues from 27 patients with SCZ and 31 age-matched control individuals, followed by comprehensive and strict analysis of somatic mutations, including mutagenesis signature, substitution patterns, and involved pathways. In particular, we explored the impact of deleterious mutations in GRIN2B through primary neural culture.

RESULTS

We identified an average of 4.9 and 5.6 somatic mutations per exome per brain in patients with SCZ and control individuals, respectively. These mutations presented with average variant allele frequencies of 8.0% in patients with SCZ and 7.6% in control individuals. Although mutational profiles, such as the number and type of mutations, showed no significant difference between patients with SCZ and control individuals, somatic mutations in SCZ brains were significantly enriched for SCZ-related pathways, including dopamine receptor, glutamate receptor, and long-term potentiation pathways. Furthermore, we showed that brain somatic mutations in GRIN2B (encoding glutamate ionotropic NMDA receptor subunit 2B), which were found in two patients with SCZ, disrupted the location of GRIN2B across the surface of dendrites among primary cultured neurons.

CONCLUSIONS

Taken together, this study shows that brain somatic mutations are associated with the pathogenesis of SCZ.

Genetic association of FMRP targets with psychiatric disorders

Genes encoding the mRNA targets of fragile X mental retardation protein (FMRP) are enriched for genetic association with psychiatric disorders. However, many FMRP targets possess functions that are themselves genetically associated with psychiatric disorders, including synaptic transmission and plasticity, making it unclear whether the genetic risk is truly related to binding by FMRP or is alternatively mediated by the sampling of genes better characterised by another trait or functional annotation. Using published common variant, rare coding variant and copy number variant data, we examined the relationship between FMRP binding and genetic association with schizophrenia, major depressive disorder and bipolar disorder. High-confidence targets of FMRP, derived from studies of multiple tissue types, were enriched for common schizophrenia risk alleles, as well as rare loss-of-function and de novo nonsynonymous variants in schizophrenia cases. Similarly, through common variation, FMRP targets were associated with major depressive disorder, and we present novel evidence of association with bipolar disorder. These relationships could not be explained by other functional annotations known to be associated with psychiatric disorders, including those related to synaptic structure and function. This study reinforces the evidence that targeting by FMRP captures a subpopulation of genes enriched for genetic association with a range of psychiatric disorders.

CHD4 variants are associated with childhood idiopathic epilepsy with sinus arrhythmia

Aims

CHD4 gene, encoding chromodomain helicase DNA‐binding protein 4, is a vital gene for fetal development. In this study, we aimed to explore the association between CHD4 variants and idiopathic epilepsy.

Methods

Trios‐based whole‐exome sequencing was performed in a cohort of 482 patients with childhood idiopathic epilepsy. The Clinical Validity Framework of ClinGen and an evaluating method from five clinical‐genetic aspects were used to determine the association between CHD4 variants and epilepsy.

Results

Four novel heterozygous missense mutations in CHD4, including two de novo mutations (c.1597A>G/p.K533E and c.4936G>A/p.E1646K) and two inherited mutations with co‐segregation (c.856C>G/p.P286A and c.4977C>G/p.D1659E), were identified in four unrelated families with eight individuals affected. Seven affected individuals had sinus arrhythmia. From the molecular sub‐regional point of view, the missense mutations located in the central regions from SNF2‐like region to DUF1087 domain were associated with multisystem developmental disorders, while idiopathic epilepsy‐related mutations were outside this region. Strong evidence from ClinGen Clinical Validity Framework and evidences from four of the five clinical‐genetic aspects suggested an association between CHD4 variants and epilepsy.

Conclusions

CHD4 was potentially a candidate pathogenic gene of childhood idiopathic epilepsy with arrhythmia. The molecular sub‐regional effect of CHD4 mutations helped explaining the mechanisms underlying phenotypic variations.

Genetic investigations of 100 inherited cardiac disease-related genes in deceased individuals with schizophrenia

Cardiac diseases and sudden cardiac death (SCD) are more prevalent in individuals diagnosed with schizophrenia compared to the general population, with especially coronary artery disease (CAD) as the major cardiovascular cause of death. Antipsychotic medications, genetics, and lifestyle factors may contribute to the increased SCD in individuals with schizophrenia. The role of antipsychotic medications and lifestyle factors have been widely investigated, while the genetic predisposition to inherited cardiac diseases in schizophrenia is poorly understood. In this study, we examined 100 genes associated with inherited cardiomyopathies and cardiac channelopathies in 97 deceased individuals diagnosed with schizophrenia for the prevalence of genetic variants associated with SCD. The deceased individuals had various causes of death and were included in the SURVIVE project, a prospective, autopsy-based study of mentally ill individuals in Denmark. This is the first study of multiple inherited cardiac disease-related genes in deceased individuals with diagnosed schizophrenia to shed light on the genetic predisposition to SCD in individuals with schizophrenia. We found no evidence for an overrepresentation of rare variants with high penetrance in inherited cardiac diseases, following the American College of Medical Genetics and Genomics and the Association for Molecular Pathology (ACMG) consensus guidelines. However, we found that the deceased individuals had a statistically significantly increased polygenic burden caused by variants in the investigated heart genes compared to the general population. This indicates that common variants with smaller effects in heart genes may play a role in schizophrenia.

Integrative Analysis Identified Key Schizophrenia Risk Factors from an Abnormal Behavior Mouse Gene Set

Schizophrenia (SCZ) is a severe chronic psychiatric illness with heterogeneous symptoms. However, the pathogenesis of SCZ is unclear, and the number of well-defined SCZ risk factors is limited. We hypothesized that an abnormal behavior (AB) gene set verified by mouse model experiments can be used to better understand SCZ risks. In this work, we carried out an integrative bioinformatics analysis to study two types of risk genes that are either differentially expressed (DEGs) in the case-control study data or carry reported SCZ genetic variants (MUTs). Next, we used RNA-Seq expression data from the hippocampus (HIPPO) and dorsolateral prefrontal cortex (DLPFC) to define the key genes affected by different types (DEGs and MUTs) in different brain regions (DLPFC and HIPPO): DLPFC-kDEG, DLPFC-kMUT, HIPPO-kDEG, and HIPPO-kMUT. The four hub genes (SHANK1, SHANK2, DLG4, and NLGN3) of the biological functionally enriched terms were strongly linked to SCZ via gene co-expression network analysis. Then, we observed that specific spatial expressions of DLPFC-kMUT and HIPPO-kMUT were convergent in the early stages and divergent in the later stages of development. In addition, all four types of key genes showed significantly larger average protein-protein interaction degrees than the background. Comparing the different cell types, the expression of four types of key genes showed specificity in different dimensions. Together, our results offer new insights into potential risk factors and help us understand the complexity and regional heterogeneity of SCZ.

Positional cloning and comprehensive mutation analysis identified a novel KDM2B mutation in a Japanese family with minor malformations, intellectual disability, and schizophrenia

The importance of epigenetic control in the development of the central nervous system has recently been attracting attention. Methylation patterns of lysine 4 and lysine 36 in histone H3 (H3K4 and H3K36) in the central nervous system are highly conserved among species. Numerous complications of body malformations and neuropsychiatric disorders are due to abnormal histone H3 methylation modifiers. In this study, we analyzed a Japanese family with a dominant inheritance of symptoms including Marfan syndrome-like minor physical anomalies (MPAs), intellectual disability, and schizophrenia (SCZ). We performed positional cloning for this family using a single nucleotide polymorphism (SNP) array and whole-exome sequencing, which revealed a missense coding strand mutation (rs1555289644, NM_032590.4: c.2173G>A, p.A725T) in exon 15 on the plant homeodomain of the KDM2B gene as a possible cause of the disease in the family. The exome sequencing revealed that within the coding region, only a point mutation in KDM2B was present in the region with the highest logarithm of odds score of 2.41 resulting from whole genome linkage analysis. Haplotype analysis revealed co-segregation with four affected family members (IV-9, III-4, IV-5, and IV-8). Lymphoblastoid cell lines from the proband with this mutation showed approximately halved KDM2B expression in comparison with healthy controls. KDM2B acts as an H3K4 and H3K36 histone demethylase. Our findings suggest that haploinsufficiency of KDM2B in the process of development, like other H3K4 and H3K36 methylation modifiers, may have caused MPAs, intellectual disability, and SCZ in this Japanese family.

Association analysis of potentially functional variants within 8p12 with schizophrenia in the Han Chinese population

OBJECTIVES

Chromosome 8p12 was first identified as a schizophrenia (SCZ) risk locus in Chinese populations and replicated in European populations. However, the underlying functional variants still need to be further explored. In this study, we sought to identify plausible causal variants within this locus.

METHODS

A total of 386 potentially functional variants from 29 genes within the 8p12 locus were analysed in 2403 SCZ cases and 2594 control subjects in the Han Chinese population using Affymetrix customised genotyping assays. SHEsisplus was used for association analysis. A multiple testing corrected p value (false discovery rate (FDR)) < .05 was considered significant, and an unadjusted p value < .05 was considered nominal evidence of an association.

RESULTS

We did not find significant associations between the tested variants and SCZ. However, nominal associations were found for rs201292574 (unadjusted p = .033, FDR p = .571; 95% confidence interval (CI): 0.265-0.945; TACC1, NP_006274.2:p.Ala211Thr) and rs45563241 (unadjusted p = .039, FDR p = .571; 95% CI: 1.023-1.866; a synonymous mutation in ADRB3).

CONCLUSIONS

Our results provide limited evidence for the associations between variants from protein coding regions in 8p12 and SCZ in the Chinese population. Analyses of both coding and regulatory variants in larger sample sizes are required to further clarify the causal variants for SCZ with this risk locus.

Inferring the molecular and phenotypic impact of amino acid variants with MutPred2

Identifying pathogenic variants and underlying functional alterations is challenging. To this end, we introduce MutPred2, a tool that improves the prioritization of pathogenic amino acid substitutions over existing methods, generates molecular mechanisms potentially causative of disease, and returns interpretable pathogenicity score distributions on individual genomes. Whilst its prioritization performance is state-of-the-art, a distinguishing feature of MutPred2 is the probabilistic modeling of variant impact on specific aspects of protein structure and function that can serve to guide experimental studies of phenotype-altering variants. We demonstrate the utility of MutPred2 in the identification of the structural and functional mutational signatures relevant to Mendelian disorders and the prioritization of de novo mutations associated with complex neurodevelopmental disorders. We then experimentally validate the functional impact of several variants identified in patients with such disorders. We argue that mechanism-driven studies of human inherited disease have the potential to significantly accelerate the discovery of clinically actionable variants.

The Role of Brain Microvascular Endothelial Cell and Blood-Brain Barrier Dysfunction in Schizophrenia

BACKGROUND

Despite decades of research, little clarity exists regarding pathogenic mechanisms related to schizophrenia. Investigations on the disease biology of schizophrenia have primarily focused on neuronal alterations. However, there is substantial evidence pointing to a significant role for the brain's microvasculature in mediating neuroinflammation in schizophrenia.

SUMMARY

Brain microvascular endothelial cells (BMEC) are a central element of the microvasculature that forms the blood-brain barrier (BBB) and shields the brain against toxins and immune cells via paracellular, transcellular, transporter, and extracellular matrix proteins. While evidence for BBB dysfunction exists in brain disorders, including schizophrenia, it is not known if BMEC themselves are functionally compromised and lead to BBB dysfunction.

KEY MESSAGES

Genome-wide association studies, postmortem investigations, and gene expression analyses have provided some insights into the role of the BBB in schizophrenia pathophysiology. However, there is a significant gap in our understanding of the role that BMEC play in BBB dysfunction. Recent advances differentiating human BMEC from induced pluripotent stem cells (iPSC) provide new avenues to examine the role of BMEC in BBB dysfunction in schizophrenia.

Genes regulated by BCL11B during T-cell development are enriched for de novo mutations found in schizophrenia patients

While abnormal neurodevelopment contributes to schizophrenia (SCZ) risk, there is also evidence to support a role for immune dysfunction in SCZ. BCL11B, associated with SCZ in genome-wide association study (GWAS), is a transcription factor that regulates the differentiation and development of cells in the central nervous and immune systems. Here, we use functional genomics data from studies of BCL11B to investigate the contribution of neuronal and immune processes to SCZ pathophysiology. We identified the gene targets of BCL11B in brain striatal cells (n = 223 genes), double negative 4 (DN4) developing T cells (n = 114 genes) and double positive (DP) developing T cells (n = 518 genes) using an integrated analysis of RNA-seq and ChIP-seq data. No gene-set was enriched for genes containing common variants associated with SCZ but the DP gene-set was enriched for genes containing missense de novo mutations (DNMs; p = .001) using data from 3,447 SCZ trios. Post hoc analysis revealed the enrichment to be stronger for DP genes negatively regulated by BCL11B. Biological processes enriched for genes negatively regulated by BCL11B in DP gene-set included immune system development and cytokine signaling. These analyses, leveraging a GWAS-identified SCZ risk gene and data on gene expression and transcription factor binding, indicate that DNMs in immune pathways contribute to SCZ risk.

Identification of de novo mutations in prenatal neurodevelopment-associated genes in schizophrenia in two Han Chinese patient-sibling family-based cohorts

Schizophrenia (SCZ) is a severe psychiatric disorder with a strong genetic component. High heritability of SCZ suggests a major role for transmitted genetic variants. Furthermore, SCZ is also associated with a marked reduction in fecundity, leading to the hypothesis that alleles with large effects on risk might often occur de novo. In this study, we conducted whole-genome sequencing for 23 families from two cohorts with unaffected siblings and parents. Two nonsense de novo mutations (DNMs) in GJC1 and HIST1H2AD were identified in SCZ patients. Ten genes (DPYSL2, NBPF1, SDK1, ZNF595, ZNF718, GCNT2, SNX9, AACS, KCNQ1, and MSI2) were found to carry more DNMs in SCZ patients than their unaffected siblings by burden test. Expression analyses indicated that these DNM implicated genes showed significantly higher expression in prefrontal cortex in prenatal stage. The DNM in the GJC1 gene is highly likely a loss function mutation (pLI = 0.94), leading to the dysregulation of ion channel in the glutamatergic excitatory neurons. Analysis of rare variants in independent exome sequencing dataset indicates that GJC1 has significantly more rare variants in SCZ patients than in unaffected controls. Data from genome-wide association studies suggested that common variants in the GJC1 gene may be associated with SCZ and SCZ-related traits. Genes co-expressed with GJC1 are involved in SCZ, SCZ-associated pathways, and drug targets. These evidences suggest that GJC1 may be a risk gene for SCZ and its function may be involved in prenatal and early neurodevelopment, a vulnerable period for developmental disorders such as SCZ.

Investigation of Schizophrenia with Human Induced Pluripotent Stem Cells

Schizophrenia is a chronic and severe neuropsychiatric condition manifested by cognitive, emotional, affective, perceptual, and behavioral abnormalities. Despite decades of research, the biological substrates driving the signs and symptoms of the disorder remain elusive, thus hampering progress in the development of treatments aimed at disease etiologies. The recent emergence of human induced pluripotent stem cell (hiPSC)-based models has provided the field with a highly innovative approach to generate, study, and manipulate living neural tissue derived from patients, making possible the exploration of fundamental roles of genes and early-life stressors in disease-relevant cell types. Here, we begin with a brief overview of the clinical, epidemiological, and genetic aspects of the condition, with a focus on schizophrenia as a neurodevelopmental disorder. We then highlight relevant technical advancements in hiPSC models and assess novel findings attained using hiPSC-based approaches and their implications for disease biology and treatment innovation. We close with a critical appraisal of the developments necessary for both further expanding knowledge of schizophrenia and the translation of new insights into therapeutic innovations.

Gene4Denovo: an integrated database and analytic platform for de novo mutations in humans

De novo mutations (DNMs) significantly contribute to sporadic diseases, particularly in neuropsychiatric disorders. Whole-exome sequencing (WES) and whole-genome sequencing (WGS) provide effective methods for detecting DNMs and prioritizing candidate genes. However, it remains a challenge for scientists, clinicians, and biologists to conveniently access and analyse data regarding DNMs and candidate genes from scattered publications. To fill the unmet need, we integrated 580 799 DNMs, including 30 060 coding DNMs detected by WES/WGS from 23 951 individuals across 24 phenotypes and prioritized a list of candidate genes with different degrees of statistical evidence, including 346 genes with false discovery rates <0.05. We then developed a database called Gene4Denovo (http://www.genemed.tech/gene4denovo/), which allowed these genetic data to be conveniently catalogued, searched, browsed, and analysed. In addition, Gene4Denovo integrated data from >60 genomic sources to provide comprehensive variant-level and gene-level annotation and information regarding the DNMs and candidate genes. Furthermore, Gene4Denovo provides end-users with limited bioinformatics skills to analyse their own genetic data, perform comprehensive annotation, and prioritize candidate genes using custom parameters. In conclusion, Gene4Denovo conveniently allows for the accelerated interpretation of DNM pathogenicity and the clinical implication of DNMs in humans.

Whole-exome sequencing of 81 individuals from 27 multiply affected bipolar disorder families

Bipolar disorder (BD) is a highly heritable neuropsychiatric disease characterized by recurrent episodes of depression and mania. Research suggests that the cumulative impact of common alleles explains 25-38% of phenotypic variance, and that rare variants may contribute to BD susceptibility. To identify rare, high-penetrance susceptibility variants for BD, whole-exome sequencing (WES) was performed in three affected individuals from each of 27 multiply affected families from Spain and Germany. WES identified 378 rare, non-synonymous, and potentially functional variants. These spanned 368 genes, and were carried by all three affected members in at least one family. Eight of the 368 genes harbored rare variants that were implicated in at least two independent families. In an extended segregation analysis involving additional family members, five of these eight genes harbored variants showing full or nearly full cosegregation with BD. These included the brain-expressed genes RGS12 and NCKAP5, which were considered the most promising BD candidates on the basis of independent evidence. Gene enrichment analysis for all 368 genes revealed significant enrichment for four pathways, including genes reported in de novo studies of autism (p_adj < 0.006) and schizophrenia (p_adj = 0.015). These results suggest a possible genetic overlap with BD for autism and schizophrenia at the rare-sequence-variant level. The present study implicates novel candidate genes for BD development, and may contribute to an improved understanding of the biological basis of this common and often devastating disease.

Exome sequencing in schizophrenia-affected parent-offspring trios reveals risk conferred by protein-coding de novo mutations

Protein-coding de novo mutations (DNMs) are significant risk factors in many neurodevelopmental disorders, whereas schizophrenia (SCZ) risk associated with DNMs has thus far been shown to be modest. We analyzed DNMs from 1,695 SCZ-affected trios and 1,077 published SCZ-affected trios to better understand the contribution to SCZ risk. Among 2,772 SCZ probands, exome-wide DNM burden remained modest. Gene set analyses revealed that SCZ DNMs were significantly concentrated in genes that were highly expressed in the brain, that were under strong evolutionary constraint and/or overlapped with genes identified in other neurodevelopmental disorders. No single gene surpassed exome-wide significance; however, 16 genes were recurrently hit by protein-truncating DNMs, corresponding to a 3.15-fold higher rate than the mutation model expectation (permuted 95% confidence interval: 1-10 genes; permuted P = 3 × 10-5). Overall, DNMs explain a small fraction of SCZ risk, and larger samples are needed to identify individual risk genes, as coding variation across many genes confers risk for SCZ in the population.

Transcriptome analysis of fibroblasts from schizophrenia patients reveals differential expression of schizophrenia-related genes

Schizophrenia is a complex neurodevelopmental disorder with high rate of morbidity and mortality. While the heritability rate is high, the precise etiology is still unknown. Although schizophrenia is a central nervous system disorder, studies using peripheral tissues have also been established to search for patient specific biomarkers and to increase understanding of schizophrenia etiology. Among all peripheral tissues, fibroblasts stand out as they are easy to obtain and culture. Furthermore, they keep genetic stability for long period and exhibit molecular similarities to cells from nervous system. Using a unique set of fibroblast samples from a genetically isolated population in northern Sweden, we performed whole transcriptome sequencing to compare differentially expressed genes in seven controls and nine patients. We found differential fibroblast expression between cases and controls for 48 genes, including eight genes previously implicated in schizophrenia or schizophrenia related pathways; HGF, PRRT2, EGR1, EGR3, C11orf87, TLR3, PLEKHH2 and PIK3CD. Weighted gene correlation network analysis identified three differentially co-expressed networks of genes significantly-associated with schizophrenia. All three modules were significantly suppressed in patients compared to control, with one module highly enriched in genes involved in synaptic plasticity, behavior and synaptic transmission. In conclusion, our results support the use of fibroblasts for identification of differentially expressed genes in schizophrenia and highlight dysregulation of synaptic networks as an important mechanism in schizophrenia.

Novel genetic susceptibility loci identified by family based whole exome sequencing in Han Chinese schizophrenia patients

Schizophrenia (SCZ) is a highly heritable psychiatric disorder that affects approximately 1% of population around the world. However, early relevant studies did not reach clear conclusions of the genetic mechanisms of SCZ, suggesting that additional susceptibility loci that exert significant influence on SCZ are yet to be revealed. So, in order to identify novel susceptibility genes that account for the genetic risk of SCZ, we performed a systematic family-based study using whole exome sequencing (WES) in 65 Han Chinese families. The analysis of 51 SCZ trios with both unaffected parents identified 22 exonic and 1 splice-site de novo mutations (DNMs) on a total of 23 genes, and showed that 12 genes carried rare protein-altering compound heterozygous mutations in more than one trio. In addition, we identified 26 exonic or splice-site single nucleotide polymorphisms (SNPs) on 18 genes with nominal significance (P < 5 × 10-4) using a transmission disequilibrium test (TDT) in all the families. Moreover, TDT result confirmed a SCZ susceptibility locus on 3p21.1, encompassing the multigenetic region NEK4-ITIH1-ITIH3-ITIH4. Through several different strategies to predict the potential pathogenic genes in silico, we revealed 4 previous discovered susceptibility genes (TSNARE1, PBRM1, STAB1 and OLIG2) and 4 novel susceptibility loci (PSEN1, TLR5, MGAT5B and SSPO) in Han Chinese SCZ patients. In summary, we identified a list of putative candidate genes for SCZ using a family-based WES approach, thus improving our understanding of the pathology of SCZ and providing critical clues to future functional validation.

Synaptic and Gene Regulatory Mechanisms in Schizophrenia, Autism, and 22q11.2 Copy Number Variant-Mediated Risk for Neuropsychiatric Disorders

BACKGROUND

22q11.2 copy number variants are among the most highly penetrant genetic risk variants for developmental neuropsychiatric disorders such as schizophrenia (SCZ) and autism spectrum disorder (ASD). However, the specific mechanisms through which they confer risk remain unclear.

METHODS

Using a functional genomics approach, we integrated transcriptomic data from the developing human brain, genome-wide association findings for SCZ and ASD, protein interaction data, and gene expression signatures from SCZ and ASD postmortem cortex to 1) organize genes into the developmental cellular and molecular systems within which they operate, 2) identify neurodevelopmental processes associated with polygenic risk for SCZ and ASD across the allelic frequency spectrum, and 3) elucidate pathways and individual genes through which 22q11.2 copy number variants may confer risk for each disorder.

RESULTS

Polygenic risk for SCZ and ASD converged on partially overlapping neurodevelopmental modules involved in synaptic function and transcriptional regulation, with ASD risk variants additionally enriched for modules involved in neuronal differentiation during fetal development. The 22q11.2 locus formed a large protein network during development that disproportionately affected SCZ-associated and ASD-associated neurodevelopmental modules, including loading highly onto synaptic and gene regulatory pathways. SEPT5, PI4KA, and SNAP29 genes are candidate drivers of 22q11.2 synaptic pathology relevant to SCZ and ASD, and DGCR8 and HIRA are candidate drivers of disease-relevant alterations in gene regulation.

CONCLUSIONS

This approach offers a powerful framework to identify neurodevelopmental processes affected by diverse risk variants for SCZ and ASD and elucidate mechanisms through which highly penetrant, multigene copy number variants contribute to disease risk.

Integrating Clinical Data and Imputed Transcriptome from GWAS to Uncover Complex Disease Subtypes: Applications in Psychiatry and Cardiology

Classifying subjects into clinically and biologically homogeneous subgroups will facilitate the understanding of disease pathophysiology and development of targeted prevention and intervention strategies. Traditionally, disease subtyping is based on clinical characteristics alone, but subtypes identified by such an approach may not conform exactly to the underlying biological mechanisms. Very few studies have integrated genomic profiles (e.g., those from GWASs) with clinical symptoms for disease subtyping. Here we proposed an analytic framework capable of finding complex diseases subgroups by leveraging both GWAS-predicted gene expression levels and clinical data by a multi-view bicluster analysis. This approach connects SNPs to genes via their effects on expression, so the analysis is more biologically relevant and interpretable than a pure SNP-based analysis. Transcriptome of different tissues can also be readily modeled. We also proposed various evaluation metrics for assessing clustering performance. Our framework was able to subtype schizophrenia subjects into diverse subgroups with different prognosis and treatment response. We also applied the framework to the Northern Finland Birth Cohort (NFBC) 1966 dataset and identified high and low cardiometabolic risk subgroups in a gender-stratified analysis. The prediction strength by cross-validation was generally greater than 80%, suggesting good stability of the clustering model. Our results suggest a more data-driven and biologically informed approach to defining metabolic syndrome and subtyping psychiatric disorders. Moreover, we found that the genes "blindly" selected by the algorithm are significantly enriched for known susceptibility genes discovered in GWASs of schizophrenia or cardiovascular diseases. The proposed framework opens up an approach to subject stratification.

Double hits in schizophrenia

The co-occurrence of a copy number variant (CNV) and a functional variant on the other allele may be a relevant genetic mechanism in schizophrenia. We hypothesized that the cumulative burden of such double hits-in particular those composed of a deletion and a coding single-nucleotide variation (SNV)-is increased in patients with schizophrenia. We combined CNV data with coding variants data in 795 patients with schizophrenia and 474 controls. To limit false CNV-detection, only CNVs called by two algorithms were included. CNV-affected genes were subsequently examined for coding SNVs, which we termed "CNV-SNVs." Correcting for total queried sequence, we assessed the CNV-SNV-burden and the combined predicted deleterious effect. We estimated P-values by permutation of the phenotype. We detected 105 CNV-SNVs; 67 in duplicated and 38 in deleted genic sequence. Although the difference in CNV-SNVs rates was not significant, the combined deleteriousness inferred by CNV-SNVs in deleted sequence was almost 4-fold higher in cases compared with controls (nominal P = 0.009). This effect may be driven by a higher number of CNV-SNVs and/or by a higher degree of predicted deleteriousness of CNV-SNVs. No such effect was observed for duplications. We provide early evidence that deletions co-occurring with a functional variant may be relevant, albeit of modest impact, for the genetic etiology of schizophrenia. Large-scale consortium studies are required to validate our findings. Sequence-based analyses would provide the best resolution for detection of CNVs as well as coding variants genome-wide.

Whole-Exome Sequencing Identifies Damaging de novo Variants in Anencephalic Cases

BACKGROUND

Anencephaly is a lethal neural tube defect (NTD). Although variants in several genes have been implicated in the development of anencephaly, a more complete picture of variants in the genome, especially de novo variants (DNVs), remains unresolved. We aim to identify DNVs that play an important role in the development of anencephaly by performing whole-exome DNA sequencing (WES) of proband-parent trios.

RESULTS

A total of 13 DNVs were identified in 8 anencephaly trios by WES, including two loss of function (LoF) variants detected in pLI > 0.9 genes (SPHKAP, c.2629_2633del, and NCOR1, p.Y1907X). Damaging DNVs were identified in 61.5% (8/13) of the anencephalic cases. Independent validation was conducted in an additional 502 NTD cases. Gene inactivation using targeted morpholino antisense oligomers and rescue assays were conducted in zebrafish, and transfection expression in HEK293T cells. Four DNVs in four cases were identified and predicted to alter protein function, including p.R328Q in WD repeat domain phosphoinositide-interacting 1 (WIPI1). Three variants, p.G313R, p.T418M, and p.L406P, in the WIPI1 gene were identified from the independent replication cohort consisting of 502 cases. Functional analysis suggested that the wipi1 p.L406P and p.R328Q variants most likely displayed loss-of-function effects during embryonic development.

CONCLUSION

De novo damaging variants are the main culprit for majority of anencephalic cases. Missense variants in WIPI1 may play a role in the genetic etiology of anencephaly, and LoF variants in SPHKAP and NCOR1 may also contribute to anencephaly. These findings add to our existing understanding of the genetic mechanisms of NTD formation.

Targeted exon sequencing in deceased schizophrenia patients in Denmark

Schizophrenia patients have higher mortality rates and lower life expectancy than the general population. However, forensic investigations of their deaths often fail to determine the cause of death, hindering prevention. As schizophrenia is a highly heritable condition and given recent advances in our understanding of the genetics of schizophrenia, it is now possible to investigate how genetic factors may contribute to mortality. We made use of findings from genome-wide association studies (GWAS) to design a targeted panel (PsychPlex) for sequencing of exons of 451 genes near index single nucleotide polymorphisms (SNPs) identified with GWAS. We sequenced the DNA of 95 deceased schizophrenia patients included in SURVIVE, a prospective, autopsy-based study of mentally ill persons in Denmark. We compared the allele frequencies of 1039 SNPs in these cases with the frequencies of 2000 Danes without psychiatric diseases and calculated their deleteriousness (CADD) scores. For 81 SNPs highly associated with schizophrenia and CADD scores above 15, expression profiles in the Genotype-Tissue Expression (GTEx) Project indicated that these variants were in exons, whose expressions are increased in several types of brain tissues, particularly in the cerebellum. Molecular pathway analysis indicated the involvement of 163 different pathways. As for rare SNP variants, most variants were scored as either benign or likely benign with an average of 17 variants of unknown significance per individual and no pathogenic variant. Our results highlight the potential of DNA sequencing of an exon panel to discover genetic factors that may be involved in the development of schizophrenia.

Plastome sequencing of Myripnois dioica and comparison within Asteraceae

Myripnois is a monotypic shrub genus in the daisy family constricted to northern China. Although wild populations of Myripnois dioica are relatively rare, this plant may potentially be cultured as a fine ornamental. In the present study, we sequenced the complete plastome of M. dioica, generating the first plastome sequences of the subfamily Pertyoideae. The plastome of M. dioica has a typical quadripartite circular structure. A large ∼20-kb and a small ∼3-kb inversion were detected in the large single copy (LSC) region and shared by other Asteraceae species. Plastome phylogenomic analyses based on 78 Asteraceae species and three outgroups revealed four groups, corresponding to four Asteraceae subfamilies: Asteroideae, Cichorioideae, Pertyoideae and Carduoideae. Among these four subfamilies, Pertyoideae is sister to Asteroideae + Cichorioideae; Carduoideae is the most basal clade. In addition, we characterized 13 simple sequence repeats (SSRs) that may be useful in future studies on population genetics.

Dendritic structural plasticity and neuropsychiatric disease

The structure of neuronal circuits that subserve cognitive functions in the brain is shaped and refined throughout development and into adulthood. Evidence from human and animal studies suggests that the cellular and synaptic substrates of these circuits are atypical in neuropsychiatric disorders, indicating that altered structural plasticity may be an important part of the disease biology. Advances in genetics have redefined our understanding of neuropsychiatric disorders and have revealed a spectrum of risk factors that impact pathways known to influence structural plasticity. In this Review, we discuss the importance of recent genetic findings on the different mechanisms of structural plasticity and propose that these converge on shared pathways that can be targeted with novel therapeutics.

Temporal dynamics of miRNAs in human DLPFC and its association with miRNA dysregulation in schizophrenia

Brain development is dependent on programmed gene expression, which is both genetically and epigenetically regulated. Post-transcriptional regulation of gene expression by microRNAs (miRNAs) is essential for brain development. As abnormal brain development is hypothesized to be associated with schizophrenia, miRNAs are an intriguing target for this disorder. The aims of this study were to determine the temporal dynamics of miRNA expression in the human dorsolateral prefrontal cortex (DLPFC), and the relationship between miRNA's temporal expression pattern and dysregulation in schizophrenia. This study used next-generation sequencing to characterize the temporal dynamics of miRNA expression in the DLPFC of 109 normal subjects (second trimester-74 years of age) and miRNA expression changes in 34 schizophrenia patients. Unlike mRNAs, the majority of which exhibits a wave of change in fetuses, most miRNAs are preferentially expressed during a certain period before puberty. It is noted that in schizophrenia patients, miRNAs normally enriched in infants tend to be upregulated, while those normally enriched in prepuberty tend to be downregulated, and the targets of these miRNAs are enriched for genes encoding synaptic proteins and those associated with schizophrenia. In addition, miR-936 and miR-3162 were found to be increased in the DLPFC of patients with schizophrenia. These findings reveal the temporal dynamics of miRNAs in the human DLPFC, implicate the importance of miRNAs in DLPFC development, and suggest a possible link between schizophrenia and dysregulation of miRNAs enriched in infancy and prepuberty.

Biomarker Identification Through Integrating fMRI and Epigenetics

OBJECTIVE

Integration of multiple datasets is a hot topic in many fields. When studying complex mental disorders, great effort has been dedicated to fusing genetic and brain imaging data. However, an increasing number of studies have pointed out the importance of epigenetic factors in the cause of psychiatric diseases. In this study, we endeavor to fill the gap by combining epigenetics (e.g., DNA methylation) with imaging data (e.g., fMRI) to identify biomarkers for schizophrenia (SZ).

METHODS

We propose to combine linear regression with canonical correlation analysis (CCA) in a relaxed yet coupled manner to extract discriminative features for SZ that are co-expressed in the fMRI and DNA methylation data.

RESULT

After validation through simulations, we applied our method to real imaging epigenetics data of 184 subjects from the Mental Illness and Neuroscience Discovery Clinical Imaging Consortium. After significance test, we identified 14 brain regions and 44 cytosine-phosphate-guanine(CpG) sites. Average classification accuracy is [Formula: see text]. By linking the CpG sites to genes, we identified pathways Guanosine ribonucleotides de novo biosynthesis and Guanosine nucleotides de novo biosynthesis, and a GO term Perikaryon.

CONCLUSION

This imaging epigenetics study has identified both brain regions and genes that are associated with neuron development and memory processing. These biomarkers contribute to a good understanding of the mechanism underlying SZ but are overlooked by previous imaging genetics studies.

SIGNIFICANCE

Our study sheds light on the understanding and diagnosis of SZ with a imaging epigenetics approach, which is demonstrated to be effective in extracting novel biomarkers associated with SZ.