Increased burden of ultra-rare protein-altering variants among 4,877 individuals with schizophrenia

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.

DECO: a framework for jointly analyzing de novo and rare case/control variants, and biological pathways

MOTIVATION

Rare variant-based analyses are beginning to identify risk genes for neuropsychiatric disorders and other diseases. However, the identified genes only account for a fraction of predicted causal genes. Recent studies have shown that rare damaging variants are significantly enriched in specific gene-sets. Methods which are able to jointly model rare variants and gene-sets to identify enriched gene-sets and use these enriched gene-sets to prioritize additional risk genes could improve understanding of the genetic architecture of diseases.

RESULTS

We propose DECO (Integrated analysis of de novo mutations, rare case/control variants and omics information via gene-sets), an integrated method for rare-variant and gene-set analysis. The method can (i) test the enrichment of gene-sets directly within the statistical model, and (ii) use enriched gene-sets to rank existing genes and prioritize additional risk genes for tested disorders. In simulations, DECO performs better than a homologous method that uses only variant data. To demonstrate the application of the proposed protocol, we have applied this approach to rare-variant datasets of schizophrenia. Compared with a method which only uses variant information, DECO is able to prioritize additional risk genes.

AVAILABILITY

DECO can be used to analyze rare-variants and biological pathways or cell types for any disease. The package is available on Github https://github.com/hoangtn/DECO.

Investigating rare pathogenic/likely pathogenic exonic variation in bipolar disorder

Bipolar disorder (BD) is a serious mental illness with substantial common variant heritability. However, the role of rare coding variation in BD is not well established. We examined the protein-coding (exonic) sequences of 3,987 unrelated individuals with BD and 5,322 controls of predominantly European ancestry across four cohorts from the Bipolar Sequencing Consortium (BSC). We assessed the burden of rare, protein-altering, single nucleotide variants classified as pathogenic or likely pathogenic (P-LP) both exome-wide and within several groups of genes with phenotypic or biologic plausibility in BD. While we observed an increased burden of rare coding P-LP variants within 165 genes identified as BD GWAS regions in 3,987 BD cases (meta-analysis OR = 1.9, 95% CI = 1.3-2.8, one-sided p = 6.0 × 10-4), this enrichment did not replicate in an additional 9,929 BD cases and 14,018 controls (OR = 0.9, one-side p = 0.70). Although BD shares common variant heritability with schizophrenia, in the BSC sample we did not observe a significant enrichment of P-LP variants in SCZ GWAS genes, in two classes of neuronal synaptic genes (RBFOX2 and FMRP) associated with SCZ or in loss-of-function intolerant genes. In this study, the largest analysis of exonic variation in BD, individuals with BD do not carry a replicable enrichment of rare P-LP variants across the exome or in any of several groups of genes with biologic plausibility. Moreover, despite a strong shared susceptibility between BD and SCZ through common genetic variation, we do not observe an association between BD risk and rare P-LP coding variants in genes known to modulate risk for SCZ.

Rare variants in the endocytic pathway are associated with Alzheimer's disease, its related phenotypes, and functional consequences

Late-onset Alzheimer's disease (LOAD) is the most common type of dementia causing irreversible brain damage to the elderly and presents a major public health challenge. Clinical research and genome-wide association studies have suggested a potential contribution of the endocytic pathway to AD, with an emphasis on common loci. However, the contribution of rare variants in this pathway to AD has not been thoroughly investigated. In this study, we focused on the effect of rare variants on AD by first applying a rare-variant gene-set burden analysis using genes in the endocytic pathway on over 3,000 individuals with European ancestry from three large whole-genome sequencing (WGS) studies. We identified significant associations of rare-variant burden within the endocytic pathway with AD, which were successfully replicated in independent datasets. We further demonstrated that this endocytic rare-variant enrichment is associated with neurofibrillary tangles (NFTs) and age-related phenotypes, increasing the risk of obtaining severer brain damage, earlier age-at-onset, and earlier age-of-death. Next, by aggregating rare variants within each gene, we sought to identify single endocytic genes associated with AD and NFTs. Careful examination using NFTs revealed one significantly associated gene, ANKRD13D. To identify functional associations, we integrated bulk RNA-Seq data from over 600 brain tissues and found two endocytic expression genes (eGenes), HLA-A and SLC26A7, that displayed significant influences on their gene expressions. Differential expressions between AD patients and controls of these three identified genes were further examined by incorporating scRNA-Seq data from 48 post-mortem brain samples and demonstrated distinct expression patterns across cell types. Taken together, our results demonstrated strong rare-variant effect in the endocytic pathway on AD risk and progression and functional effect of gene expression alteration in both bulk and single-cell resolution, which may bring more insight and serve as valuable resources for future AD genetic studies, clinical research, and therapeutic targeting.

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.

Meta-Analysis of Brain Gene Expression Data from Mouse Model Studies of Maternal Immune Activation Using Poly(I:C)

Maternal immune activation (MIA) is a known risk factor for schizophrenia (SCZ) and autism spectrum disorder (ASD) and is often modelled in animal studies in order to study the effect of prenatal infection on brain function including behaviour and gene expression. Although the effect of MIA on gene expression are highly heterogeneous, combining data from multiple gene expression studies in a robust method may shed light on the true underlying biological effects caused by MIA and this could inform studies of SCZ and ASD. This study combined four RNA-seq and microarray datasets in an overlap analysis and ranked meta-analysis in order to investigate genes, pathways and cell types dysregulated in the MIA mouse models. Genes linked to SCZ and ASD and crucial in neurodevelopmental processes including neural tube folding, regulation of cellular stress and neuronal/glial cell differentiation were among the most consistently dysregulated in these ranked analyses. Gene ontologies including K+ ion channel function, neuron and glial cell differentiation, synaptic structure, axonal outgrowth, cilia function and lipid metabolism were also strongly implicated. Single-cell analysis identified excitatory and inhibitory cell types in the cortex, hippocampus and striatum that may be affected by MIA and are also enriched for genes associated with SCZ, ASD and cognitive phenotypes. This points to the cellular location of molecular mechanisms that may be consistent between the MIA model and neurodevelopmental disease, improving our understanding of its utility to study prenatal infection as an environmental stressor.

De novo variation in bipolar disorder

Bipolar disorder (BD) is a common, highly heritable disorder that affects 1-2% of the world's population. To date, most genetic studies of BD have focused on common gene variation, and while robustly associated loci have been identified, a substantial proportion of the heritability remains missing and could be partially attributable to rare variation. In this study, we apply a de novo paradigm in BD to identify newly arisen variants that have yet to undergo natural selection and may represent highly pathogenic variants. We performed whole genome sequencing of 97 trios of Ashkenazi Jewish descent, selecting "simplex" families with no family history of BD and an early age of onset. We found a total of 6882 de novo variants (an average of 70.9 ± 12.9 S.D. variants per trio), including 107 variants within protein-coding genes. We combined our exonic variations with the results of 79 previously published BD trios, identifying 20 loss-of-function (LoF) and 77 missense damaging de novo variants in BD. These variants showed significant enrichment for constrained genes and for genes located to the postsynaptic density (PSD) (all Bonferroni corrected p < 0.05). Pathway analyses showed enrichment in several pathways, including "Phosphoinositides (PI) and their downstream targets" (Bonferroni p = 4.2 × 10-6), a pathway prominently featured in lithium's hypothesized mechanism of action. In addition, while we found overall evidence for transmission of common variant polygenic risk of BD in our full sample (pTDT p = 2.21 × 10-4), specific trios with LoF variants showed no evidence of polygenic transmission. In sum, our findings support the de novo paradigm as a contributor to the genetic architecture of BD and provide evidence that constrained genes, as well as genes within the PSD and PI pathway harbor rare variation associated with BD.

Genetic contributors to risk of schizophrenia in the presence of a 22q11.2 deletion

Schizophrenia occurs in about one in four individuals with 22q11.2 deletion syndrome (22q11.2DS). The aim of this International Brain and Behavior 22q11.2DS Consortium (IBBC) study was to identify genetic factors that contribute to schizophrenia, in addition to the ~20-fold increased risk conveyed by the 22q11.2 deletion. Using whole-genome sequencing data from 519 unrelated individuals with 22q11.2DS, we conducted genome-wide comparisons of common and rare variants between those with schizophrenia and those with no psychotic disorder at age ≥25 years. Available microarray data enabled direct comparison of polygenic risk for schizophrenia between 22q11.2DS and independent population samples with no 22q11.2 deletion, with and without schizophrenia (total n = 35,182). Polygenic risk for schizophrenia within 22q11.2DS was significantly greater for those with schizophrenia (padj = 6.73 × 10-6). Novel reciprocal case-control comparisons between the 22q11.2DS and population-based cohorts showed that polygenic risk score was significantly greater in individuals with psychotic illness, regardless of the presence of the 22q11.2 deletion. Within the 22q11.2DS cohort, results of gene-set analyses showed some support for rare variants affecting synaptic genes. No common or rare variants within the 22q11.2 deletion region were significantly associated with schizophrenia. These findings suggest that in addition to the deletion conferring a greatly increased risk to schizophrenia, the risk is higher when the 22q11.2 deletion and common polygenic risk factors that contribute to schizophrenia in the general population are both present.

Increased schizophrenia family history burden and reduced premorbid IQ in treatment-resistant schizophrenia: a Swedish National Register and Genomic Study

A high proportion of those with schizophrenia experience treatment non-response, placing them at higher risk for mortality and suicide attempts, compared to treatment responders. The clinical, social, and economic burden of treatment-resistant schizophrenia (TRS) are substantial. Previous genomic and epidemiological studies of TRS were often limited by sample size or lack of comprehensive genomic data. We aimed to systematically understand the clinical, demographic, and genomic correlates of TRS using epidemiological and genetic epidemiological modelling in a Swedish national population sample (n = 24,706) and then in a subgroup with common variant genetic risk scores, rare copy-number variant burden, and rare exonic burden (n = 4936). Population-based analyses identified increasing schizophrenia family history to be significantly associated with TRS (highest quartile of familial burden vs. lowest: adjusted odds ratio (aOR): 1.31, P = 4.8 × 10^-8). In males, a decrease of premorbid IQ of one standard deviation was significantly associated with greater risk of TRS (minimal aOR: 0.94, P = 0.002). In a subset of cases with extensive genomic data, we found no significant association between the genetic risk scores of four psychiatric disorders and two cognitive traits with TRS (schizophrenia genetic risk score: aOR = 1.07, P = 0.067). The association between copy number variant and rare variant burden measures and TRS did not reach the pre-defined statistical significance threshold (all P ≥ 0.005). In conclusion, direct measures of genomic risk were not associated with TRS; however, premorbid IQ in males and schizophrenia family history were significantly correlated with TRS and points to new insights into the architecture of TRS.

Unveiling the Pathogenesis of Psychiatric Disorders Using Network Models

Psychiatric disorders are complex brain disorders with a high degree of genetic heterogeneity, affecting millions of people worldwide. Despite advances in psychiatric genetics, the underlying pathogenic mechanisms of psychiatric disorders are still largely elusive, which impedes the development of novel rational therapies. There has been accumulating evidence suggesting that the genetics of complex disorders can be viewed through an omnigenic lens, which involves contextualizing genes in highly interconnected networks. Thus, applying network-based multi-omics integration methods could cast new light on the pathophysiology of psychiatric disorders. In this review, we first provide an overview of the recent advances in psychiatric genetics and highlight gaps in translating molecular associations into mechanistic insights. We then present an overview of network methodologies and review previous applications of network methods in the study of psychiatric disorders. Lastly, we describe the potential of such methodologies within a multi-tissue, multi-omics approach, and summarize the future directions in adopting diverse network approaches.

Single molecule in situ hybridization reveals distinct localizations of schizophrenia risk-related transcripts SNX19 and AS3MT in human brain

Genome-wide association studies have identified single nucleotide polymorphisms (SNPs) associated with schizophrenia risk. Integration of RNA-sequencing data from postmortem human brains with these risk SNPs identified transcripts associated with increased schizophrenia susceptibility, including a class of exon 9-spliced isoforms of Sorting nexin-19 (SNX19d9) and an isoform of Arsenic methyltransferase (AS3MT) splicing out exons 2 and 3 (AS3MTd2d3). However, the biological function of these transcript variants is unclear. Defining the cell types where these risk transcripts are dominantly expressed is an important step to understand function, in prioritizing specific cell types and/or neural pathways in subsequent studies. To identify the cell type-specific localization of SNX19 and AS3MT in the human dorsolateral prefrontal cortex (DLPFC), we used single-molecule in situ hybridization techniques combined with automated quantification and machine learning approaches to analyze 10 postmortem brains of neurotypical individuals. These analyses revealed that both pan-SNX19 and pan-AS3MT were more highly expressed in neurons than non-neurons in layers II/III and VI of DLPFC. Furthermore, pan-SNX19 was preferentially expressed in glutamatergic neurons, while pan-AS3MT was preferentially expressed in GABAergic neurons. Finally, we utilized duplex BaseScope technology, to delineate the localization of SNX19d9 and AS3MTd2d3 splice variants, revealing consistent trends in spatial gene expression among pan-transcripts and schizophrenia risk-related transcript variants. These findings demonstrate that schizophrenia risk transcripts have distinct localization patterns in the healthy human brains, and suggest that SNX19 transcripts might disrupt the normal function of glutamatergic neurons, while AS3MT may lead to disturbances in the GABAergic system in the pathophysiology of schizophrenia.

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.

Resolving cell state in iPSC-derived human neural samples with multiplexed fluorescence imaging

Human induced pluripotent stem cell-derived (iPSC) neural cultures offer clinically relevant models of human diseases, including Amyotrophic Lateral Sclerosis, Alzheimer’s, and Autism Spectrum Disorder. In situ characterization of the spatial-temporal evolution of cell state in 3D culture and subsequent 2D dissociated culture models based on protein expression levels and localizations is essential to understanding neural cell differentiation, disease state phenotypes, and sample-to-sample variability. Here, we apply PRobe-based Imaging for Sequential Multiplexing (PRISM) to facilitate multiplexed imaging with facile, rapid exchange of imaging probes to analyze iPSC-derived cortical and motor neuron cultures that are relevant to psychiatric and neurodegenerative disease models, using over ten protein targets. Our approach permits analysis of cell differentiation, cell composition, and functional marker expression in complex stem-cell derived neural cultures. Furthermore, our approach is amenable to automation, offering in principle the ability to scale-up to dozens of protein targets and samples.

Tomov et al. utilize DNA-PRISM to allow for multiplexed imaging of cultured cells using antibodies modified with oligonucleotide probes. The differentiation of iPSCs to cortical and motor neurons is characterized in model cultures, relevant for use in disease research and drug screening.

A patient-based model of RNA mis-splicing uncovers treatment targets in Parkinson's disease

Parkinson's disease (PD) is a heterogeneous neurodegenerative disorder with monogenic forms representing prototypes of the underlying molecular pathology and reproducing to variable degrees the sporadic forms of the disease. Using a patient-based in vitro model of PARK7-linked PD, we identified a U1-dependent splicing defect causing a drastic reduction in DJ-1 protein and, consequently, mitochondrial dysfunction. Targeting defective exon skipping with genetically engineered U1-snRNA recovered DJ-1 protein expression in neuronal precursor cells and differentiated neurons. After prioritization of candidate drugs, we identified and validated a combinatorial treatment with the small-molecule compounds rectifier of aberrant splicing (RECTAS) and phenylbutyric acid, which restored DJ-1 protein and mitochondrial dysfunction in patient-derived fibroblasts as well as dopaminergic neuronal cell loss in mutant midbrain organoids. Our analysis of a large number of exomes revealed that U1 splice-site mutations were enriched in sporadic PD patients. Therefore, our study suggests an alternative strategy to restore cellular abnormalities in in vitro models of PD and provides a proof of concept for neuroprotection based on precision medicine strategies in PD.

Systematic analysis of exonic germline and postzygotic de novo mutations in bipolar disorder

Bipolar disorder is a severe mental illness characterized by recurrent manic and depressive episodes. To better understand its genetic architecture, we analyze ultra-rare de novo mutations in 354 trios with bipolar disorder. For germline de novo mutations, we find significant enrichment of loss-of-function mutations in constrained genes (corrected-P = 0.0410) and deleterious mutations in presynaptic active zone genes (FDR = 0.0415). An analysis integrating single-cell RNA-sequencing data identifies a subset of excitatory neurons preferentially expressing the genes hit by deleterious mutations, which are also characterized by high expression of developmental disorder genes. In the analysis of postzygotic mutations, we observe significant enrichment of deleterious ones in developmental disorder genes (P = 0.00135), including the SRCAP gene mutated in two unrelated probands. These data collectively indicate the contributions of both germline and postzygotic mutations to the risk of bipolar disorder, supporting the hypothesis that postzygotic mutations of developmental disorder genes may contribute to bipolar disorder.

The significance of rare and de novo variants in bipolar disorder is not well understood. Here, the authors have analyzed whole exome/genome data from trios to identify deleterious de novo variants associated with bipolar disorder.

A tipping point in neuropsychiatric genetics

Severe neuropsychiatric disorders are so genetically heterogeneous that virtually every unrelated patient harbors different clinically significant alleles. By studying schizophrenia in the Ashkenazi Jewish founder population, Lencz and co-authors identified rare severe alleles each shared by a few patients. Experimental evaluation of an implicated protocadherin allele revealed failure to form homophilic cellular aggregates as a possible mechanism for defective development of neural circuits.

Brain gene co-expression networks link complement signaling with convergent synaptic pathology in schizophrenia

The most significant common variant association for schizophrenia (SCZ) reflects increased expression of the complement component 4A (C4A). Yet, it remains unclear how C4A interacts with other SCZ risk genes or whether the complement system more broadly is implicated in SCZ pathogenesis. Here, we integrate several existing, large-scale genetic and transcriptomic datasets to interrogate the functional role of the complement system and C4A in the human brain. Unexpectedly, we find no significant genetic enrichment among known complement system genes for SCZ. Conversely, brain co-expression network analyses using C4A as a seed gene reveal that genes downregulated when C4A expression increases exhibit strong and specific genetic enrichment for SCZ risk. This convergent genomic signal reflects synaptic processes, is sexually dimorphic and most prominent in frontal cortical brain regions, and is accentuated by smoking. Overall, these results indicate that synaptic pathways-rather than the complement system-are the driving force conferring SCZ risk.

Dysregulation of post-transcriptional modification by copy number variable microRNAs in schizophrenia with enhanced glycation stress

Previously, we identified a subpopulation of schizophrenia (SCZ) showing increased levels of plasma pentosidine, a marker of glycation and oxidative stress. However, its causative genetic factors remain largely unknown. Recently, it has been suggested that dysregulated posttranslational modification by copy number variable microRNAs (CNV-miRNAs) may contribute to the etiology of SCZ. Here, an integrative genome-wide CNV-miRNA analysis was performed to investigate the etiology of SCZ with accumulated plasma pentosidine (PEN-SCZ). The number of CNV-miRNAs and the gene ontology (GO) in the context of miRNAs within CNVs were compared between PEN-SCZ and non-PEN-SCZ groups. Gene set enrichment analysis of miRNA target genes was further performed to evaluate the pathways affected in PEN-SCZ. We show that miRNAs were significantly enriched within CNVs in the PEN-SCZ versus non-PEN-SCZ groups (p = 0.032). Of note, as per GO analysis, the dysregulated neurodevelopmental events in the two groups may have different origins. Additionally, gene set enrichment analysis of miRNA target genes revealed that miRNAs involved in glycation/oxidative stress and synaptic neurotransmission, especially glutamate/GABA receptor signaling, were possibly affected in PEN-SCZ. To the best of our knowledge, this is the first genome-wide CNV-miRNA study suggesting the role of CNV-miRNAs in the etiology of PEN-SCZ, through effects on genes related to glycation/oxidative stress and synaptic function. Our findings provide supportive evidence that glycation/oxidative stress possibly caused by genetic defects related to the posttranscriptional modification may lead to synaptic dysfunction. Therefore, targeting miRNAs may be one of the promising approaches for the treatment of PEN-SCZ.

Downregulation by CNNM2 of ATP5MD expression in the 10q24.32 schizophrenia-associated locus involved in impaired ATP production and neurodevelopment

Genome-wide association studies (GWAS) have accelerated the discovery of numerous genetic variants associated with schizophrenia. However, most risk variants show a small effect size (odds ratio (OR) <1.2), suggesting that more functional risk variants remain to be identified. Here, we employed region-based multi-marker analysis of genomic annotation (MAGMA) to identify additional risk loci containing variants with large OR value from Psychiatry Genomics Consortium (PGC2) schizophrenia GWAS data and then employed summary-data-based mendelian randomization (SMR) to prioritize schizophrenia susceptibility genes. The top-ranked susceptibility gene ATP5MD, encoding an ATP synthase membrane subunit, is observed to be downregulated in schizophrenia by the risk allele of CNNM2-rs1926032 in the schizophrenia-associated 10q24.32 locus. The Atp5md knockout (KO) in mice was associated with abnormal startle reflex and gait, and ATP5MD knockdown (KD) in human induced pluripotent stem cell-derived neurons disrupted the neural development and mitochondrial respiration and ATP production. Moreover, CNNM2-rs1926032 KO could induce downregulation of ATP5MD expression and disruptions of mitochondrial respiration and ATP production. This study constitutes an important mechanistic component that links schizophrenia-associated CNNM2 regions to disruption in energy adenosine system modulation and neuronal function by long-distance chromatin domain downregulation of ATP5MD. This pathogenic mechanism provides therapeutic implications for schizophrenia.

More evidence on TRIO missense mutations in the spectrin repeat domain causing severe developmental delay and recognizable facial dysmorphism with macrocephaly

TRIO is a Dbl family guanine nucleotide exchange factor (GEF) and an important regulator of neuronal development. Most truncating and missense variants affecting the Dbl homology domain of TRIO are associated with a neurodevelopmental disorder with microcephaly (MIM617061). Recently, de novo missense variants affecting the spectrin repeat region of TRIO were associated with a novel phenotype comprising severe developmental delay and macrocephaly (MIM618825). Here, we provide more evidence on this new TRIO-associated phenotype by reporting two severely affected probands with de novo missense variants in TRIO affecting the spectrin repeat region upstream of the typically affected GEF1 domain of the protein.

Refinement of the clinical variant interpretation framework by statistical evidence and machine learning

BACKGROUND

Although the American College of Medical Genetics and Genomics/Association for Molecular Pathology (ACMG/AMP) guidelines for variant interpretation are used widely in clinical genetics, there is room for improvement of these knowledge-based guidelines.

METHODS

Statistical assessment of average deleteriousness of start-lost, stop-lost, and in-frame insertion and deletion (indel) variants and extraction of deleterious subsets was performed, being informed by proportions of rare variants in the general population of the Genome Aggregation Database (gnomAD). A machine learning-based model scoring the pathogenicity of start-lost variants (the PoStaL model) was constructed by predicting possible translation initiation sites on transcripts by deep learning and training a random forest on known pathogenic and likely benign variants.

FINDINGS

The proportion of rare variants was highest in stop-lost variants, followed by in-frame indels and start-lost variants, suggesting that the criteria in the ACMG/AMP guidelines assigning PVS (pathogenic very strong) to start-lost variants and PM (pathogenic moderate) to stop-lost and in-frame indel variants would not be appropriate. Regarding deleterious subsets, stop-lost variants introducing extensions of more than 30 amino acids and in-frame indels computationally predicted to be damaging are enriched for rare and known pathogenic variants. For start-lost variants, we developed the PoStaL model, which outperforms existing tools. We also provide comprehensive lists of the PoStaL scores for start-lost variants and the length of extended amino acids by stop-lost variants.

CONCLUSIONS

Our study could contribute to refinement of the ACMG/AMP guidelines, provides resources for future investigation, and provides an example of how to improve knowledge-based frameworks by data-driven approaches.

FUNDING

The study was supported by grants from the Japan Agency for Medical Research and Development (AMED) and the Japan Society for the Promotion of Science (JSPS).

Revisiting Schizophrenia from an Evolutionary Perspective: An Association Study of Recent Evolutionary Markers and Schizophrenia

The persistence of schizophrenia in human populations at a high prevalence and with a large heritability estimate despite reduced fertility and increased mortality rate is a Darwinian paradox. This may be likely if the genomic components that predispose to schizophrenia are also advantageous for the acquisition of important human traits, such as language and cognition. Accordingly, an emerging group of genomic markers of recent evolution in humans, namely human accelerated regions (HARs), since our divergence from chimpanzees, are gaining importance for neurodevelopmental disorders, such as schizophrenia. We hypothesize that variants within HARs may affect the expression of genes under their control, thus contributing to disease etiology. A total of 49 HAR single nucleotide polymorphisms (SNPs) were prioritized from the complete repertoire of HARs (n = 2737) based on their functional relevance and prevalence in the South Asian population. Test of association using 2 independent schizophrenia case-control cohorts of north Indian ethnicity (discovery: n = 930; replication: n = 1104) revealed 3 SNPs (rs3800926, rs3801844, and rs764453) from chromosome 7 and rs77047799 from chromosome 3 to be significantly associated (combined analysis: Bonferroni corrected P < .002-.000004). Of note, these SNPs were found to alter the expression of neurodevelopmental genes such as SLC25A13, MAD1L1, and ULK4; a few from the HOX gene family; and a few genes that are implicated in mitochondrial function. These SNPs may most likely alter binding sites of transcription factors, including TFCP2, MAFK, SREBF2, E2F1, and/or methylation signatures around these genes. These findings reiterate a neurodevelopmental basis of schizophrenia and also open up a promising avenue to investigate HAR-mediated mitochondrial dysfunction in schizophrenia etiology.

Familial Psychosis Associated With a Missense Mutation at MACF1 Gene Combined With the Rare Duplications DUP3p26.3 and DUP16q23.3, Affecting the CNTN6 and CDH13 Genes

Psychosis is a highly heritable and heterogeneous psychiatric condition. Its genetic architecture is thought to be the result of the joint effect of common and rare variants. Families with high prevalence are an interesting approach to shed light on the rare variant's contribution without the need of collecting large cohorts. To unravel the genomic architecture of a family enriched for psychosis, with four affected individuals, we applied a system genomic approach based on karyotyping, genotyping by whole-exome sequencing to search for rare single nucleotide variants (SNVs) and SNP array to search for copy-number variants (CNVs). We identified a rare non-synonymous variant, g.39914279 C > G, in the MACF1 gene, segregating with psychosis. Rare variants in the MACF1 gene have been previously detected in SCZ patients. Besides, two rare CNVs, DUP3p26.3 and DUP16q23.3, were also identified in the family affecting relevant genes (CNTN6 and CDH13, respectively). We hypothesize that the co-segregation of these duplications with the rare variant g.39914279 C > G of MACF1 gene precipitated with schizophrenia and schizoaffective disorder.

Structural interaction between DISC1 and ATF4 underlying transcriptional and synaptic dysregulation in an iPSC model of mental disorders

Psychiatric disorders are a collection of heterogeneous mental disorders arising from a contribution of genetic and environmental insults, many of which molecularly converge on transcriptional dysregulation, resulting in altered synaptic functions. The underlying mechanisms linking the genetic lesion and functional phenotypes remain largely unknown. Patient iPSC-derived neurons with a rare frameshift DISC1 (Disrupted-in-schizophrenia 1) mutation have previously been shown to exhibit aberrant gene expression and deficits in synaptic functions. How DISC1 regulates gene expression is largely unknown. Here we show that Activating Transcription Factor 4 (ATF4), a DISC1 binding partner, is more abundant in the nucleus of DISC1 mutant human neurons and exhibits enhanced binding to a collection of dysregulated genes. Functionally, overexpressing ATF4 in control neurons recapitulates deficits seen in DISC1 mutant neurons, whereas transcriptional and synaptic deficits are rescued in DISC1 mutant neurons with CRISPR-mediated heterozygous ATF4 knockout. By solving the high-resolution atomic structure of the DISC1-ATF4 complex, we show that mechanistically, the mutation of DISC1 disrupts normal DISC1-ATF4 interaction, and results in excessive ATF4 binding to DNA targets and deregulated gene expression. Together, our study identifies the molecular and structural basis of an DISC1-ATF4 interaction underlying transcriptional and synaptic dysregulation in an iPSC model of mental disorders.

Association of Damaging Variants in Genes With Increased Cancer Risk Among Patients With Congenital Heart Disease

Importance

Patients with congenital heart disease (CHD), the most common birth defect, have increased risks for cancer. Identification of the variables that contribute to cancer risk is essential for recognizing patients with CHD who warrant longitudinal surveillance and early interventions.

Objective

To compare the frequency of damaging variants in cancer risk genes among patients with CHD and control participants and identify associated clinical variables in patients with CHD who have cancer risk variants.

Design, Setting, and Participants

This multicenter case-control study included participants with CHD who had previously been recruited to the Pediatric Cardiac Genomics Consortium based on presence of structural cardiac anomaly without genetic diagnosis at the time of enrollment. Permission to use published sequencing data from unaffected adult participants was obtained from 2 parent studies. Data were collected for this study from December 2010 to April 2019.

Exposures

Presence of rare (allele frequency, <1 × 10-5) loss-of-function (LoF) variants in cancer risk genes.

Main Outcomes and Measures

Frequency of LoF variants in cancer risk genes (defined in the Catalogue of Somatic Mutations in Cancer-Cancer Gene Consensus database), were statistically assessed by binomial tests in patients with CHD and control participants.

Results

A total of 4443 individuals with CHD (mean [range] age, 13.0 [0-84] years; 2225 of 3771 with reported sex [59.0%] male) and 9808 control participants (mean [range] age, 52.1 [1-92] years; 4967 of 9808 [50.6%] male) were included. The frequency of LoF variants in regulatory cancer risk genes was significantly higher in patients with CHD than control participants (143 of 4443 [3.2%] vs 166 of 9808 [1.7%]; odds ratio [OR], 1.93 [95% CI, 1.54-2.42]; P = 1.38 × 10-12), and among CHD genes previously associated with cancer risk (58 of 4443 [1.3%] vs 18 of 9808 [0.18%]; OR, 7.2 [95% CI, 4.2-12.2]; P < 2.2 × 10-16). The LoF variants were also nominally increased in 14 constrained cancer risk genes with high expression in the developing heart. Seven of these genes (ARHGEF12, CTNNB1, LPP, MLLT4, PTEN, TCF12, and TFRC) harbored LoF variants in multiple patients with unexplained CHD. The highest rates for LoF variants in cancer risk genes occurred in patients with CHD and extracardiac anomalies (248 of 1482 individuals [16.7%]; control: 1099 of 9808 individuals [11.2%]; OR, 1.59 [95% CI, 1.37-1.85]; P = 1.3 × 10-10) and/or neurodevelopmental delay (209 of 1393 individuals [15.0%]; control: 1099 of 9808 individuals [11.2%]; OR, 1.40 [95% CI, 1.19-1.64]; P = 9.6 × 10-6).

Conclusions and Relevance

Genotypes of CHD may account for increased cancer risks. In this cohort, damaging variants were prominent in the 216 genes that predominantly encode regulatory proteins. Consistent with their fundamental developmental functions, patients with CHD and damaging variants in these genes often had extracardiac manifestations. These data may also implicate cancer risk genes that are repeatedly varied in patients with unexplained CHD as CHD genes.

Dysregulated Prefrontal Cortical RhoA Signal Transduction in Bipolar Disorder with Psychosis: New Implications for Disease Pathophysiology

While research has identified alterations in dorsolateral prefrontal cortical function as a key factor to the etiology of bipolar disorder, few studies have uncovered robust changes in protein signal transduction pathways in this disorder. Given the direct relevance of protein-based expressional alterations to cellular functions and because many of the key regulatory mechanisms for the disease pathogenesis likely include alterations in protein activity rather than changes in expression alone, the identification of alterations in discrete signal transduction pathways in bipolar disorder would have broad implications for understanding the disease pathophysiology. As prior microarray data point to a previously unrecognized involvement of the RhoA network in bipolar disorder, here we investigate the protein expression and activity of key components of a RhoA signal transduction pathway in dorsolateral prefrontal cortical homogenates from subjects with bipolar disorder. The results of this investigation implicate overactivation of prefrontal cortical RhoA signaling in specific subtypes of bipolar disorder. The specificity of these findings is demonstrated by a lack of comparable changes in schizophrenia; however, our findings do identify convergence between both disorders at the level of activity-mediated actin cytoskeletal regulation. These findings have implications for understanding the altered cortical synaptic connectivity of bipolar disorder.

Novel ultra-rare exonic variants identified in a founder population implicate cadherins in schizophrenia

The identification of rare variants associated with schizophrenia has proven challenging due to genetic heterogeneity, which is reduced in founder populations. In samples from the Ashkenazi Jewish population, we report that schizophrenia cases had a greater frequency of novel missense or loss of function (MisLoF) ultra-rare variants (URVs) compared to controls, and the MisLoF URV burden was inversely correlated with polygenic risk scores in cases. Characterizing 141 "case-only" genes (MisLoF URVs in ≥3 cases with none in controls), the cadherin gene set was associated with schizophrenia. We report a recurrent case mutation in PCDHA3 that results in the formation of cytoplasmic aggregates and failure to engage in homophilic interactions on the plasma membrane in cultured cells. Modeling purifying selection, we demonstrate that deleterious URVs are greatly overrepresented in the Ashkenazi population, yielding enhanced power for association studies. Identification of the cadherin/protocadherin family as risk genes helps specify the synaptic abnormalities central to schizophrenia.

The association between family history and genomic burden with schizophrenia mortality: a Swedish population-based register and genetic sample study

Individuals with schizophrenia (SCZ) have a 2-3-fold higher risk of mortality than the general population. Heritability of mortality in psychiatric disorders has been proposed; however, few have investigated SCZ family history and genetic variation, with all-cause and specific causes of death. We aimed to identify correlates of SCZ mortality using genetic epidemiological and genetic modelling in two samples: a Swedish national population sample and a genotyped subsample. In the Swedish national population sample followed from the first SCZ treatment contact until emigration, death or end of the follow-up, we investigated a standardised measure of SCZ family history. In a subgroup with comprehensive genetic data, we investigated the impact of common and rare genetic variation. Cox proportional hazards regression was used to estimate the association between various factors and mortality (all and specific causes). A total of 13727 SCZ cases fulfilled criteria for the population-based analyses (1268 deaths, 9.2%). The genomic subset contained 4991 cases (1353 deaths, 27.1%). Somatic mutations associated with clonal hematopoiesis with unknown drivers were associated with all-cause mortality (HR 1.77, 95% CI: 1.26-2.49). No other heritable measures were associated with all-cause mortality nor with any specific causes of death. Future studies in larger, comparable cohorts are warranted to further understand the association between hereditary measures and mortality in SCZ.

Tissue-specific enhancer functional networks for associating distal regulatory regions to disease

Systematic study of tissue-specific function of enhancers and their disease associations is a major challenge. We present an integrative machine-learning framework, FENRIR, that integrates thousands of disparate epigenetic and functional genomics datasets to infer tissue-specific functional relationships between enhancers for 140 diverse human tissues and cell types, providing a regulatory-region-centric approach to systematically identify disease-associated enhancers. We demonstrated its power to accurately prioritize enhancers associated with 25 complex diseases. In a case study on autism, FENRIR-prioritized enhancers showed a significant proband-specific de novo mutation enrichment in a large, sibling-controlled cohort, indicating pathogenic signal. We experimentally validated transcriptional regulatory activities of eight enhancers, including enhancers not previously reported with autism, and demonstrated their differential regulatory potential between proband and sibling alleles. Thus, FENRIR is an accurate and effective framework for the study of tissue-specific enhancers and their role in disease. FENRIR can be accessed at fenrir.flatironinstitute.org/.

Potential involvement of DSCAML1 mutations in neurodevelopmental disorders

The molecular mechanisms underlying neurodevelopmental disorders (NDDs) remain unclear. We previously identified Down syndrome cell adhesion molecule like 1 (Dscaml1) as a responsible gene for Ihara epileptic rat (IER), a rat model for human NDDs with epilepsy. However, the relationship between NDDs and DSCAML1 in humans is still elusive. In this study, we screened databases of autism spectrum disorders (ASD), intellectual disability (ID)/developmental disorders (DD) and schizophrenia for genomic mutations in human DSCAML1. We then performed in silico analyses to estimate the potential damage to the mutated DSCAML1 proteins and chose three representative mutations (DSCAML1C729R , DSCAML1R1685* and DSCAML1K2108Nfs*37 ), which lacked a cysteine residue in the seventh Ig domain, the intracellular region and the C-terminal PDZ-binding motif, respectively. In overexpression experiments in a cell line, DSCAML1C729R lost its mature N-glycosylation, whereas DSCAML1K2108Nfs*37 was abnormally degraded via proteasome-dependent protein degradation. Furthermore, in primary hippocampal neurons, the ability of the wild-type DSCAML1 to regulate the number of synapses was lost with all mutant proteins. These results provide insight into understanding the roles of the domains in the DSCAML1 protein and further suggest that these mutations cause functional changes, albeit through different mechanisms, that likely affect the pathophysiology of NDDs.

A Neurodevelopment Approach for a Transitional Model of Early Onset Schizophrenia

In the last decades, the conceptualization of schizophrenia has dramatically changed, moving from a neurodegenerative process occurring in early adult life to a neurodevelopmental disorder starting be-fore birth, showing a variety of premorbid and prodromal symptoms and, in relatively few cases, evolving in the full-blown psychotic syndrome. High rates of co-occurring different neurodevelopmental disorders such as Autism spectrum disorder and ADHD, predating the onset of SCZ, and neurobio-logical underpinning with significant similarities, support the notion of a pan-developmental disturbance consisting of impairments in neuromotor, receptive language, social and cognitive development. Con-sidering that many SCZ risk factors may be similar to symptoms of other neurodevelopmental psychi-atric disorders, transition processes from child & adolescent to adult systems of care should include both high risk people as well as subject with other neurodevelopmental psychiatric disorders with different levels of severity. This descriptive mini-review discuss the need of innovative clinical approaches, re-considering specific diagnostic categories, stimulating a careful analysis of risk factors and promoting the appropriate use of new and safer medications.

A unique missense variant in the E1A-binding protein P400 gene is implicated in schizophrenia by whole-exome sequencing and mutant mouse models

Genetic and epidemiological evidence has suggested that genetic factors are important in schizophrenia, although its pathophysiology is poorly understood. This study used whole-exome sequencing to investigate potential novel schizophrenia-causing genes in a Japanese family containing several members affected by severe or treatment-resistant schizophrenia. A missense variant, chr12:132064747C>T (rs200626129, P2805L), in the E1A-binding protein P400 (EP400) gene completely segregated with schizophrenia in this family. Furthermore, numerous other EP400 mutations were identified in the targeted sequencing of a schizophrenia patient cohort. We also created two lines of Ep400 gene-edited mice, which had anxiety-like behaviours and reduced axon diameters. Our findings suggest that rs200626129 in EP400 is likely to cause schizophrenia in this Japanese family, and may lead to a better understanding and treatment of schizophrenia.

SATB2-LEMD2 interaction links nuclear shape plasticity to regulation of cognition-related genes

SATB2 is a schizophrenia risk gene and is genetically associated with human intelligence. How it affects cognition at molecular level is currently unknown. Here, we show that interactions between SATB2, a chromosomal scaffolding protein, and the inner nuclear membrane protein LEMD2 orchestrate the response of pyramidal neurons to neuronal activation. Exposure to novel environment in vivo causes changes in nuclear shape of CA1 hippocampal neurons via a SATB2-dependent mechanism. The activity-driven plasticity of the nuclear envelope requires not only SATB2, but also its protein interactor LEMD2 and the ESCRT-III/VPS4 membrane-remodeling complex. Furthermore, LEMD2 depletion in cortical neurons, similar to SATB2 ablation, affects neuronal activity-dependent regulation of multiple rapid and delayed primary response genes. In human genetic data, LEMD2-regulated genes are enriched for de novo mutations reported in intellectual disability and schizophrenia and are, like SATB2-regulated genes, enriched for common variants associated with schizophrenia and cognitive function. Hence, interactions between SATB2 and the inner nuclear membrane protein LEMD2 influence gene expression programs in pyramidal neurons that are linked to cognitive ability and psychiatric disorder etiology.

Sequence three million genomes across Africa

Capture the full scope of variation to improve health care, equity and medical research globally.

Genome sequencing broadens the range of contributing variants with clinical implications in schizophrenia

The range of genetic variation with potential clinical implications in schizophrenia, beyond rare copy number variants (CNVs), remains uncertain. We therefore analyzed genome sequencing data for 259 unrelated adults with schizophrenia from a well-characterized community-based cohort previously examined with chromosomal microarray for CNVs (none with 22q11.2 deletions). We analyzed these genomes for rare high-impact variants considered causal for neurodevelopmental disorders, including single-nucleotide variants (SNVs) and small insertions/deletions (indels), for potential clinical relevance based on findings for neurodevelopmental disorders. Also, we investigated a novel variant type, tandem repeat expansions (TREs), in 45 loci known to be associated with monogenic neurological diseases. We found several of these variants in this schizophrenia population suggesting that these variants have a wider clinical spectrum than previously thought. In addition to known pathogenic CNVs, we identified 11 (4.3%) individuals with clinically relevant SNVs/indels in genes converging on schizophrenia-relevant pathways. Clinical yield was significantly enriched in females and in those with broadly defined learning/intellectual disabilities. Genome analyses also identified variants with potential clinical implications, including TREs (one in DMPK; two in ATXN8OS) and ultra-rare loss-of-function SNVs in ZMYM2 (a novel candidate gene for schizophrenia). Of the 233 individuals with no pathogenic CNVs, we identified rare high-impact variants (i.e., clinically relevant or with potential clinical implications) for 14 individuals (6.0%); some had multiple rare high-impact variants. Mean schizophrenia polygenic risk score was similar between individuals with and without clinically relevant rare genetic variation; common variants were not sufficient for clinical application. These findings broaden the individual and global picture of clinically relevant genetic risk in schizophrenia, and suggest the potential translational value of genome sequencing as a single genetic technology for schizophrenia.

Special Article: Translational Science Update. Pharmacological Implications of Emerging Schizophrenia Genetics: Can the Bridge From 'Genomics' to 'Therapeutics' be Defined and Traversed?

Recent schizophrenia genome-wide association studies (GWAS) have identified genomic variants of common and rare frequency, significantly associated with schizophrenia. While numerous functional genomics efforts are ongoing to elucidate the biological effects of schizophrenia risk variants, a consideration of their therapeutic implications is timely and imperative, for patients as well as for an iterative effect on elucidating the underlying biology and pathophysiology of illness. The current article reviews efforts to translate emerging schizophrenia genomics into novel approaches to target discovery and therapeutic intervention. Though the path from 'genetic risk to therapy' is far from straightforward, there are provocative early possibilities that harbor the promise of treatment based on causation rather than phenomenology, as well as 'precision psychiatry,' a basis for stratifying patients to enable more precise and effective, personalized therapy.

The interaction of ASAH1 and NGF gene involving in neurotrophin signaling pathway contributes to schizophrenia susceptibility and psychopathology

The neurodevelopmental hypothesis of schizophrenia has been widely accepted. In light of our previous microarray data, two neurodevelopment-related genes were focused on inclduing the N-acylsphingosine amidohydrolase 1 gene (ASAH1) and the nerve growth factor gene (NGF). The evidence that ASAH1 and NGF are associated with schizophrenia is far from conclusive. Furthermore, their interactions in schizophrenia have not been investigated. Total 413 patients and 578 controls were included. Eleven single-nucleotide polymorphisms (SNPs) in ASAH1 and NGF were selected. A multifactor dimensionality reduction (MDR) was applied to investigate gene-gene interactions in schizophrenia, and the traditional odds ratio methods was applied to validate it. The effects of ASAH1, NGF and their interaction on the severity of the disease were analyzed by 3 × 3 covariance analysis of (ANCOVA). The biological interaction between ASAH1 and NGF was examined. KEGG was used to identify the related signaling pathways. After correction by Bonferroni, there were no differences in the genotypic, allelic, or haplotypic frequencies of 11 SNPs between patients and controls. However, the interaction of certain SNPs had effect on susceptibility to schizophrenia, including two high-risk and one low-risk genotypic combinations (OR = 1.49 [1.11-2.00]; OR = 1.45 [1.09-1.92], and OR = 0.64 [0.41-0.98]). ASAH1-rs7830490 and its interaction with NGF-rs4332358 were associated with the general psychopathological subscale score (F _adjusted = 3.94, p _adjusted = 0.01; F _adjusted = 2.36, p _adjusted = 0.03). We also found that ASAH1 and NGF interacted with CaMK2B involving in the neurotrophin signaling pathway. Our results suggest that the interaction of ASAH1 and NGF with CaMK2B involved in neurotrophin signaling pathway may contribute to schizophrenia susceptibility and psychopathology.

Diagnostic genetic testing for neurodevelopmental psychiatric disorders: closing the gap between recommendation and clinical implementation

Advances in laboratory testing have significantly increased the detection of rare genetic etiologies of neurodevelopmental psychiatric disorders (NPD), particularly developmental delay/ intellectual disability, autism spectrum disorder, and schizophrenia. Establishing a genetic diagnosis has important medical and personal utility for individuals with these conditions. Diagnostic genetic tests for NPD are clinically available but underutilized outside of medical genetics settings. Without clear multidisciplinary consensus recommendations, active involvement of medical specialists working with NPD patients, and practical education and training, the implementation of genetic testing for NPD will continue to lag behind other areas of medicine. In the long-term, collaborative efforts to address educational, logistical, and workforce obstacles will improve patient care and pave the way for targeted, effective NPD treatments.

Multiple Linear Regression Allows Weighted Burden Analysis of Rare Coding Variants in an Ethnically Heterogeneous Population

Weighted burden analysis has been used in exome-sequenced case-control studies to identify genes in which there is an excess of rare and/or functional variants associated with phenotype. Implementation in a ridge regression framework allows simultaneous analysis of all variants along with relevant covariates, such as population principal components. In order to apply the approach to a quantitative phenotype, a weighted burden score is derived for each subject and included in a linear regression analysis. The weighting scheme is adjusted in order to apply differential weights to rare and very rare variants and a score is derived based on both the frequency and predicted effect of each variant. When applied to an ethnically heterogeneous dataset consisting of 49,790 exome-sequenced UK Biobank subjects and using body mass index as the phenotype, the method produces a very inflated test statistic. However, this is almost completely corrected by including 20 population principal components as covariates. When this is done, the top 30 genes include a few which are quite plausibly associated with the phenotype, including LYPLAL1 and NSDHL. This approach offers a way to carry out gene-based analyses of rare variants identified by exome sequencing in heterogeneous datasets without requiring that data from ethnic minority subjects be discarded. This research has been conducted using the UK Biobank Resource.

Convergence and Divergence in the Genetics of Psychiatric Disorders From Pathways to Developmental Stages

In the past decade, the identification of susceptibility genes for psychiatric disorders has become routine, but understanding the biology underlying these discoveries has proven extremely difficult. The large number of potential risk genes and the genetic overlap between disorders are major obstacles for studying the etiology of these conditions. Systems biology approaches relying on gene ontologies, gene coexpression, and protein-protein interactions are used to identify convergence of the genes in relation to biological processes, cell types, brain areas, and developmental stages. Across psychiatric disorders, there is a clear enrichment for genes expressed in the brain and especially in the cortex, but a higher resolution is vastly dependent on sample size and statistical power. There is indication that susceptibility genes tend to be expressed in the brain during periods preceding the typical onset of the disorders. Thus, the role of genes in prenatal brain development is more pronounced for childhood-onset disorders, such as autism spectrum disorder and attention-deficit/hyperactivity disorder, but is much less so for bipolar disorder and depression. One of the most consistent findings across multiple disorders and classes of genetic variants is the role of genes intolerant to mutations in psychiatric disorders, yet this association is more pronounced for disorders with a clear neurodevelopmental component. Notwithstanding, a detailed understanding of the neurobiology of psychiatric disorders is still lacking. It is possible that it will only be revealed by studying the risk genes at the level of the development and function of neuronal networks and circuits.

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.

Distinct genetic patterns of shared and unique genes across four neurodevelopmental disorders

Neurodevelopmental disorders, including autism spectrum disorder (ASD), intellectual disability (ID), developmental disorders (DD) and epileptic encephalopathy (EE), have a strong clinical comorbidity, which indicates a common genetic etiology across various disorders. However, the underlying genetic mechanisms of comorbidity and specificity remain unknown across neurodevelopmental disorders. Based on de novo mutations, we compared systematically the functional characteristics between shared and unique genes under these disorders, as well as the spatiotemporal trajectory of development in brain and common molecular pathways of all shared genes. We observed that shared genes present more constrained against functional rare genetic variation, and harbor more pathogenic rare variants than do unique genes in each disorder. Furthermore, 71 shared genes formed two clusters related to synaptic transmission, transcription regulation and chromatin regulator. Particularly, we also found that two core genes STXBP1 and SCN2A, that were shared by the four neurodevelopmental disorders showed prominent pleiotropy. Our findings shed light on the shared and specific patterns across neurodevelopmental disorders and will enable us to further comprehend the etiology and provide valuable information for the diagnosis of neurodevelopmental disorders.

Transcriptomic Insight Into the Polygenic Mechanisms Underlying Psychiatric Disorders

Over the past decade, large-scale genetic studies have successfully identified hundreds of genetic variants robustly associated with risk for psychiatric disorders. However, mechanistic insight and clinical translation continue to lag the pace of risk variant identification, hindered by the sheer number of targets and their predominant noncoding localization, as well as pervasive pleiotropy and incomplete penetrance. Successful next steps require identification of "causal" genetic variants and their proximal biological consequences; placing variants within biologically defined functional contexts, reflecting specific molecular pathways, cell types, circuits, and developmental windows; and characterizing the downstream, convergent neurobiological impact of polygenicity within an individual. Here, we discuss opportunities and challenges of high-throughput transcriptomic profiling in the human brain, and how transcriptomic approaches can help pinpoint mechanisms underlying genetic risk for psychiatric disorders at a scale necessary to tackle daunting levels of polygenicity. These include transcriptome-wide association studies for risk gene prioritization through integration of genome-wide association studies with expression quantitative trait loci. We outline transcriptomic results that inform our understanding of the brain-level molecular pathology of psychiatric disorders, including autism spectrum disorder, bipolar disorder, major depressive disorder, and schizophrenia. Finally, we discuss systems-level approaches for integration of distinct genetic, genomic, and phenotypic levels, including combining spatially resolved gene expression and human neuroimaging maps. Results highlight the importance of understanding gene expression (dys)regulation across human brain development as a major contributor to psychiatric disease pathogenesis, from common variants acting as expression quantitative trait loci to rare variants enriched for gene expression regulatory pathways.

A transcriptome-wide association study implicates specific pre- and post-synaptic abnormalities in schizophrenia

Schizophrenia is a complex highly heritable disorder. Genome-wide association studies (GWAS) have identified multiple loci that influence the risk of developing schizophrenia, although the causal variants driving these associations and their impacts on specific genes are largely unknown. We identify a significant correlation between schizophrenia risk and expression at 89 genes in the dorsolateral prefrontal cortex (P ≤ 9.43 × 10-6), including 20 novel genes. Genes whose expression correlate with schizophrenia were enriched for those involved in abnormal CNS synaptic transmission (PFDR = 0.02) and antigen processing and presentation of peptide antigen via MHC class I (PFDR = 0.02). Within the CNS synaptic transmission set, we identify individual significant candidate genes to which we assign direction of expression changes in schizophrenia. The findings provide strong candidates for experimentally probing the molecular basis of synaptic pathology in schizophrenia.

Neuronal timescales are functionally dynamic and shaped by cortical microarchitecture

Complex cognitive functions such as working memory and decision-making require information maintenance over seconds to years, from transient sensory stimuli to long-term contextual cues. While theoretical accounts predict the emergence of a corresponding hierarchy of neuronal timescales, direct electrophysiological evidence across the human cortex is lacking. Here, we infer neuronal timescales from invasive intracranial recordings. Timescales increase along the principal sensorimotor-to-association axis across the entire human cortex, and scale with single-unit timescales within macaques. Cortex-wide transcriptomic analysis shows direct alignment between timescales and expression of excitation- and inhibition-related genes, as well as genes specific to voltage-gated transmembrane ion transporters. Finally, neuronal timescales are functionally dynamic: prefrontal cortex timescales expand during working memory maintenance and predict individual performance, while cortex-wide timescales compress with aging. Thus, neuronal timescales follow cytoarchitectonic gradients across the human cortex and are relevant for cognition in both short and long terms, bridging microcircuit physiology with macroscale dynamics and behavior.

Axonal branch patterning and neuronal shape diversity: roles in developmental circuit assembly: Axonal branch patterning and neuronal shape diversity in developmental circuit assembly

Recent progress in human genetics and single cell sequencing rapidly expands the list of molecular factors that offer important new contributions to our understanding of brain wiring. Yet many new molecular factors are being discovered that have never been studied in the context of neuronal circuit development. This is clearly asking for increased efforts to better understand the developmental mechanisms of circuit assembly [1]. Moreover, recent studies characterizing the developmental causes of some psychiatric diseases show impressive progress in reaching cellular resolution in their analysis. They provide concrete support emphasizing the importance of axonal branching and synapse formation as a hotspot for potential defects. Inspired by these new studies we will discuss progress but also challenges in understanding how neurite branching and neuronal shape diversity itself impacts on specificity of neuronal circuit assembly. We discuss the idea that neuronal shape acquisition itself is a key specificity factor in neuronal circuit assembly.

Rare single-nucleotide DAB1 variants and their contribution to Schizophrenia and autism spectrum disorder susceptibility

Disabled 1 (DAB1) is an intracellular adaptor protein in the Reelin signaling pathway and plays an essential role in correct neuronal migration and layer formation in the developing brain. DAB1 has been repeatedly reported to be associated with neurodevelopmental disorders including schizophrenia (SCZ) and autism spectrum disorders (ASD) in genetic, animal, and postmortem studies. Recently, increasing attention has been given to rare single-nucleotide variants (SNVs) found by deep sequencing of candidate genes. In this study, we performed exon-targeted resequencing of DAB1 in 370 SCZ and 192 ASD patients using next-generation sequencing technology to identify rare SNVs with a minor allele frequency <1%. We detected two rare missense mutations (G382C, V129I) and then performed a genetic association study in a sample comprising 1763 SCZ, 380 ASD, and 2190 healthy control subjects. Although no statistically significant association with the detected mutations was observed for either SCZ or ASD, G382C was found only in the case group, and in silico analyses and in vitro functional assays suggested that G382C alters the function of the DAB1 protein. The rare variants of DAB1 found in the present study should be studied further to elucidate their potential functional relevance to the pathophysiology of SCZ and ASD.

A rare mutation could affect nerve development in the brain, potentially increasing susceptibility to schizophrenia and autism spectrum disorder. Hiroki Kimura of Nagoya University Graduate School of Medicine and colleagues in Japan screened the protein-coding parts of the gene DAB1 in people with schizophrenia or autism for rare nucleotide variations. DAB1 is important for neurodevelopment and synaptic function, and has been associated with neurodevelopmental disorders. Their investigations revealed a mutation that can change the structure and stability of the protein coded by DAB1, potentially increasing susceptibility to neurodevelopmental disorders. The mutation was found in people with schizophrenia and autism, but not in healthy controls. Further research with larger patient samples and investigations on mice engineered to include the mutation could clarify its impact on the pathophysiology of schizophrenia and autism.

Contributions of common genetic variants to risk of schizophrenia among individuals of African and Latino ancestry

Schizophrenia is a common, chronic and debilitating neuropsychiatric syndrome affecting tens of millions of individuals worldwide. While rare genetic variants play a role in the etiology of schizophrenia, most of the currently explained liability is within common variation, suggesting that variation predating the human diaspora out of Africa harbors a large fraction of the common variant attributable heritability. However, common variant association studies in schizophrenia have concentrated mainly on cohorts of European descent. We describe genome-wide association studies of 6152 cases and 3918 controls of admixed African ancestry, and of 1234 cases and 3090 controls of Latino ancestry, representing the largest such study in these populations to date. Combining results from the samples with African ancestry with summary statistics from the Psychiatric Genomics Consortium (PGC) study of schizophrenia yielded seven newly genome-wide significant loci, and we identified an additional eight loci by incorporating the results from samples with Latino ancestry. Leveraging population differences in patterns of linkage disequilibrium, we achieve improved fine-mapping resolution at 22 previously reported and 4 newly significant loci. Polygenic risk score profiling revealed improved prediction based on trans-ancestry meta-analysis results for admixed African (Nagelkerke's R2 = 0.032; liability R2 = 0.017; P < 10-52), Latino (Nagelkerke's R2 = 0.089; liability R2 = 0.021; P < 10-58), and European individuals (Nagelkerke's R2 = 0.089; liability R2 = 0.037; P < 10-113), further highlighting the advantages of incorporating data from diverse human populations.

FMRP and CYFIP1 at the Synapse and Their Role in Psychiatric Vulnerability

There is increasing awareness of the role genetic risk variants have in mediating vulnerability to psychiatric disorders such as schizophrenia and autism. Many of these risk variants encode synaptic proteins, influencing biological pathways of the postsynaptic density and, ultimately, synaptic plasticity. Fragile-X mental retardation 1 (FMR1) and cytoplasmic fragile-X mental retardation protein (FMRP)-interacting protein 1 (CYFIP1) contain 2 such examples of highly penetrant risk variants and encode synaptic proteins with shared functional significance. In this review, we discuss the biological actions of FMRP and CYFIP1, including their regulation of (i) protein synthesis and specifically FMRP targets, (ii) dendritic and spine morphology, and (iii) forms of synaptic plasticity such as long-term depression. We draw upon a range of preclinical studies that have used genetic dosage models of FMR1 and CYFIP1 to determine their biological function. In parallel, we discuss how clinical studies of fragile X syndrome or 15q11.2 deletion patients have informed our understanding of FMRP and CYFIP1, and highlight the latest psychiatric genomic findings that continue to implicate FMRP and CYFIP1. Lastly, we assess the current limitations in our understanding of FMRP and CYFIP1 biology and how they must be addressed before mechanism-led therapeutic strategies can be developed for psychiatric disorders.

Genes influenced by MEF2C contribute to neurodevelopmental disease via gene expression changes that affect multiple types of cortical excitatory neurons

Myocyte enhancer factor 2 C (MEF2C) is an important transcription factor during neurodevelopment. Mutation or deletion of MEF2C causes intellectual disability (ID), and common variants within MEF2C are associated with cognitive function and schizophrenia risk. We investigated if genes influenced by MEF2C during neurodevelopment are enriched for genes associated with neurodevelopmental phenotypes and if this can be leveraged to identify biological mechanisms and individual brain cell types affected. We used a set of 1055 genes that were differentially expressed in the adult mouse brain following early embryonic deletion of Mef2c in excitatory cortical neurons. Using genome-wide association studies data, we found these differentially expressed genes (DEGs) to be enriched for genes associated with schizophrenia, intelligence and educational attainment but not autism spectrum disorder (ASD). For this gene set, genes that overlap with target genes of the Fragile X mental retardation protein (FMRP) are a major driver of these enrichments. Using trios data, we found these DEGs to be enriched for genes containing de novo mutations reported in ASD and ID, but not schizophrenia. Using single-cell RNA sequencing data, we identified that a number of different excitatory glutamatergic neurons in the cortex were enriched for these DEGs including deep layer pyramidal cells and cells in the retrosplenial cortex, entorhinal cortex and subiculum, and these cell types are also enriched for FMRP target genes. The involvement of MEF2C and FMRP in synapse elimination suggests that disruption of this process in these cell types during neurodevelopment contributes to cognitive function and risk of neurodevelopmental disorders.