Massively parallel functional dissection of schizophrenia-associated noncoding genetic variants.

Schizophrenia (SCZ) is a highly heritable mental disorder with thousands of associated genetic variants located mostly in the noncoding space of the genome. Translating these associations into insights regarding the underlying pathomechanisms has been challenging because the causal variants, their mechanisms of action, and their target genes remain largely unknown. We implemented a massively parallel variant annotation pipeline (MVAP) to perform SCZ variant-to-function mapping at scale in disease-relevant neural cell types. This approach identified 620 functional variants (1.7%) that operate in a highly developmental context and neuronal-activity-dependent manner. Multimodal integration of epigenomic and CRISPRi screening data enabled us to link these functional variants to target genes, biological processes, and ultimately alterations of neuronal physiology. These results provide a multistage prioritization strategy to map functional single-nucleotide polymorphism (SNP)-to-gene-to-endophenotype relations and offer biological insights into the context-dependent molecular processes modulated by SCZ-associated genetic variation.

Human brain organogenesis: Toward a cellular understanding of development and disease.

The construction of the human nervous system is a distinctly complex although highly regulated process. Human tissue inaccessibility has impeded a molecular understanding of the developmental specializations from which our unique cognitive capacities arise. A confluence of recent technological advances in genomics and stem cell-based tissue modeling is laying the foundation for a new understanding of human neural development and dysfunction in neuropsychiatric disease. Here, we review recent progress on uncovering the cellular and molecular principles of human brain organogenesis in vivo as well as using organoids and assembloids in vitro to model features of human evolution and disease.

Transforming lives through genetics.

Mary-Claire King's approach to genetics has had a major impact on breast and ovarian cancer and, more recently, mental illnesses including schizophrenia. Science writer Kendall Morgan talked with Mary-Claire, recipient of a 2021 Canada Gairdner International Award, about her life, her lengthy quest to discover the genetic basis of susceptibility to breast cancer, the struggles for women in science, and much more. An edited version of this conversation is presented below.

Phenotypic Landscape of Schizophrenia-Associated Genes Defines Candidates and Their Shared Functions.

Genomic studies have identified hundreds of candidate genes near loci associated with risk for schizophrenia. To define candidates and their functions, we mutated zebrafish orthologs of 132 human schizophrenia-associated genes. We created a phenotype atlas consisting of whole-brain activity maps, brain structural differences, and profiles of behavioral abnormalities. Phenotypes were diverse but specific, including altered forebrain development and decreased prepulse inhibition. Exploration of these datasets identified promising candidates in more than 10 gene-rich regions, including the magnesium transporter cnnm2 and the translational repressor gigyf2, and revealed shared anatomical sites of activity differences, including the pallium, hypothalamus, and tectum. Single-cell RNA sequencing uncovered an essential role for the understudied transcription factor znf536 in the development of forebrain neurons implicated in social behavior and stress. This phenotypic landscape of schizophrenia-associated genes prioritizes more than 30 candidates for further study and provides hypotheses to bridge the divide between genetic association and biological mechanism.

Genetic Control of Expression and Splicing in Developing Human Brain Informs Disease Mechanisms.

Tissue-specific regulatory regions harbor substantial genetic risk for disease. Because brain development is a critical epoch for neuropsychiatric disease susceptibility, we characterized the genetic control of the transcriptome in 201 mid-gestational human brains, identifying 7,962 expression quantitative trait loci (eQTL) and 4,635 spliceQTL (sQTL), including several thousand prenatal-specific regulatory regions. We show that significant genetic liability for neuropsychiatric disease lies within prenatal eQTL and sQTL. Integration of eQTL and sQTL with genome-wide association studies (GWAS) via transcriptome-wide association identified dozens of novel candidate risk genes, highlighting shared and stage-specific mechanisms in schizophrenia (SCZ). Gene network analysis revealed that SCZ and autism spectrum disorder (ASD) affect distinct developmental gene co-expression modules. Yet, in each disorder, common and rare genetic variation converges within modules, which in ASD implicates superficial cortical neurons. More broadly, these data, available as a web browser and our analyses, demonstrate the genetic mechanisms by which developmental events have a widespread influence on adult anatomical and behavioral phenotypes.

Genomic Dissection of Bipolar Disorder and Schizophrenia, Including 28 Subphenotypes.

Schizophrenia and bipolar disorder are two distinct diagnoses that share symptomology. Understanding the genetic factors contributing to the shared and disorder-specific symptoms will be crucial for improving diagnosis and treatment. In genetic data consisting of 53,555 cases (20,129 bipolar disorder [BD], 33,426 schizophrenia [SCZ]) and 54,065 controls, we identified 114 genome-wide significant loci implicating synaptic and neuronal pathways shared between disorders. Comparing SCZ to BD (23,585 SCZ, 15,270 BD) identified four genomic regions including one with disorder-independent causal variants and potassium ion response genes as contributing to differences in biology between the disorders. Polygenic risk score (PRS) analyses identified several significant correlations within case-only phenotypes including SCZ PRS with psychotic features and age of onset in BD. For the first time, we discover specific loci that distinguish between BD and SCZ and identify polygenic components underlying multiple symptom dimensions. These results point to the utility of genetics to inform symptomology and potential treatment.

Gene Discovery for Complex Traits: Lessons from Africa.

The genetics of African populations reveals an otherwise "missing layer" of human variation that arose between 100,000 and 5 million years ago. Both the vast number of these ancient variants and the selective pressures they survived yield insights into genes responsible for complex traits in all populations.

A concerted neuron-astrocyte program declines in ageing and schizophrenia.

Human brains vary across people and over time; such variation is not yet understood in cellular terms. Here we describe a relationship between people's cortical neurons and cortical astrocytes. We used single-nucleus RNA sequencing to analyse the prefrontal cortex of 191 human donors aged 22-97 years, including healthy individuals and people with schizophrenia. Latent-factor analysis of these data revealed that, in people whose cortical neurons more strongly expressed genes encoding synaptic components, cortical astrocytes more strongly expressed distinct genes with synaptic functions and genes for synthesizing cholesterol, an astrocyte-supplied component of synaptic membranes. We call this relationship the synaptic neuron and astrocyte program (SNAP). In schizophrenia and ageing-two conditions that involve declines in cognitive flexibility and plasticity1,2-cells divested from SNAP: astrocytes, glutamatergic (excitatory) neurons and GABAergic (inhibitory) neurons all showed reduced SNAP expression to corresponding degrees. The distinct astrocytic and neuronal components of SNAP both involved genes in which genetic risk factors for schizophrenia were strongly concentrated. SNAP, which varies quantitatively even among healthy people of similar age, may underlie many aspects of normal human interindividual differences and may be an important point of convergence for multiple kinds of pathophysiology.

Polygenic architecture of rare coding variation across 394,783 exomes.

Both common and rare genetic variants influence complex traits and common diseases. Genome-wide association studies have identified thousands of common-variant associations, and more recently, large-scale exome sequencing studies have identified rare-variant associations in hundreds of genes1-3. However, rare-variant genetic architecture is not well characterized, and the relationship between common-variant and rare-variant architecture is unclear4. Here we quantify the heritability explained by the gene-wise burden of rare coding variants across 22 common traits and diseases in 394,783 UK Biobank exomes5. Rare coding variants (allele frequency < 1 × 10-3) explain 1.3% (s.e. = 0.03%) of phenotypic variance on average-much less than common variants-and most burden heritability is explained by ultrarare loss-of-function variants (allele frequency < 1 × 10-5). Common and rare variants implicate the same cell types, with similar enrichments, and they have pleiotropic effects on the same pairs of traits, with similar genetic correlations. They partially colocalize at individual genes and loci, but not to the same extent: burden heritability is strongly concentrated in significant genes, while common-variant heritability is more polygenic, and burden heritability is also more strongly concentrated in constrained genes. Finally, we find that burden heritability for schizophrenia and bipolar disorder6,7 is approximately 2%. Our results indicate that rare coding variants will implicate a tractable number of large-effect genes, that common and rare associations are mechanistically convergent, and that rare coding variants will contribute only modestly to missing heritability and population risk stratification.

Regulation of Zbp1 by miR-99b-5p in microglia controls the development of schizophrenia-like symptoms in mice.

Current approaches to the treatment of schizophrenia have mainly focused on the protein-coding part of the genome; in this context, the roles of microRNAs have received less attention. In the present study, we analyze the microRNAome in the blood and postmortem brains of schizophrenia patients, showing that the expression of miR-99b-5p is downregulated in both the prefrontal cortex and blood of patients. Lowering the amount of miR-99b-5p in mice leads to both schizophrenia-like phenotypes and inflammatory processes that are linked to synaptic pruning in microglia. The microglial miR-99b-5p-supressed inflammatory response requires Z-DNA binding protein 1 (Zbp1), which we identify as a novel miR-99b-5p target. Antisense oligonucleotides against Zbp1 ameliorate the pathological effects of miR-99b-5p inhibition. Our findings indicate that a novel miR-99b-5p-Zbp1 pathway in microglia might contribute to the pathogenesis of schizophrenia.

Spatial and temporal mapping of de novo mutations in schizophrenia to a fetal prefrontal cortical network.

Genes disrupted in schizophrenia may be revealed by de novo mutations in affected persons from otherwise healthy families. Furthermore, during normal brain development, genes are expressed in patterns specific to developmental stage and neuroanatomical structure. We identified de novo mutations in persons with schizophrenia and then mapped the responsible genes onto transcriptome profiles of normal human brain tissues from age 13 weeks gestation to adulthood. In the dorsolateral and ventrolateral prefrontal cortex during fetal development, genes harboring damaging de novo mutations in schizophrenia formed a network significantly enriched for transcriptional coexpression and protein interaction. The 50 genes in the network function in neuronal migration, synaptic transmission, signaling, transcriptional regulation, and transport. These results suggest that disruptions of fetal prefrontal cortical neurogenesis are critical to the pathophysiology of schizophrenia. These results also support the feasibility of integrating genomic and transcriptome analyses to map critical neurodevelopmental processes in time and space in the brain.

Mapping genomic loci implicates genes and synaptic biology in schizophrenia.

Schizophrenia has a heritability of 60-80%1, much of which is attributable to common risk alleles. Here, in a two-stage genome-wide association study of up to 76,755 individuals with schizophrenia and 243,649 control individuals, we report common variant associations at 287 distinct genomic loci. Associations were concentrated in genes that are expressed in excitatory and inhibitory neurons of the central nervous system, but not in other tissues or cell types. Using fine-mapping and functional genomic data, we identify 120 genes (106 protein-coding) that are likely to underpin associations at some of these loci, including 16 genes with credible causal non-synonymous or untranslated region variation. We also implicate fundamental processes related to neuronal function, including synaptic organization, differentiation and transmission. Fine-mapped candidates were enriched for genes associated with rare disruptive coding variants in people with schizophrenia, including the glutamate receptor subunit GRIN2A and transcription factor SP4, and were also enriched for genes implicated by such variants in neurodevelopmental disorders. We identify biological processes relevant to schizophrenia pathophysiology; show convergence of common and rare variant associations in schizophrenia and neurodevelopmental disorders; and provide a resource of prioritized genes and variants to advance mechanistic studies.

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

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

CNVs: harbingers of a rare variant revolution in psychiatric genetics.

The genetic bases of neuropsychiatric disorders are beginning to yield to scientific inquiry. Genome-wide studies of copy number variation (CNV) have given rise to a new understanding of disease etiology, bringing rare variants to the forefront. A proportion of risk for schizophrenia, bipolar disorder, and autism can be explained by rare mutations. Such alleles arise by de novo mutation in the individual or in recent ancestry. Alleles can have specific effects on behavioral and neuroanatomical traits; however, expressivity is variable, particularly for neuropsychiatric phenotypes. Knowledge from CNV studies reflects the nature of rare alleles in general and will serve as a guide as we move forward into a new era of whole-genome sequencing.

Interplay between DISC1 and GABA signaling regulates neurogenesis in mice and risk for schizophrenia.

How extrinsic stimuli and intrinsic factors interact to regulate continuous neurogenesis in the postnatal mammalian brain is unknown. Here we show that regulation of dendritic development of newborn neurons by Disrupted-in-Schizophrenia 1 (DISC1) during adult hippocampal neurogenesis requires neurotransmitter GABA-induced, NKCC1-dependent depolarization through a convergence onto the AKT-mTOR pathway. In contrast, DISC1 fails to modulate early-postnatal hippocampal neurogenesis when conversion of GABA-induced depolarization to hyperpolarization is accelerated. Extending the period of GABA-induced depolarization or maternal deprivation stress restores DISC1-dependent dendritic regulation through mTOR pathway during early-postnatal hippocampal neurogenesis. Furthermore, DISC1 and NKCC1 interact epistatically to affect risk for schizophrenia in two independent case control studies. Our study uncovers an interplay between intrinsic DISC1 and extrinsic GABA signaling, two schizophrenia susceptibility pathways, in controlling neurogenesis and suggests critical roles of developmental tempo and experience in manifesting the impact of susceptibility genes on neuronal development and risk for mental disorders.

Sequencing chromosomal abnormalities reveals neurodevelopmental loci that confer risk across diagnostic boundaries.

Balanced chromosomal abnormalities (BCAs) represent a relatively untapped reservoir of single-gene disruptions in neurodevelopmental disorders (NDDs). We sequenced BCAs in patients with autism or related NDDs, revealing disruption of 33 loci in four general categories: (1) genes previously associated with abnormal neurodevelopment (e.g., AUTS2, FOXP1, and CDKL5), (2) single-gene contributors to microdeletion syndromes (MBD5, SATB2, EHMT1, and SNURF-SNRPN), (3) novel risk loci (e.g., CHD8, KIRREL3, and ZNF507), and (4) genes associated with later-onset psychiatric disorders (e.g., TCF4, ZNF804A, PDE10A, GRIN2B, and ANK3). We also discovered among neurodevelopmental cases a profoundly increased burden of copy-number variants from these 33 loci and a significant enrichment of polygenic risk alleles from genome-wide association studies of autism and schizophrenia. Our findings suggest a polygenic risk model of autism and reveal that some neurodevelopmental genes are sensitive to perturbation by multiple mutational mechanisms, leading to variable phenotypic outcomes that manifest at different life stages.

Twelve-hour rhythms in transcript expression within the human dorsolateral prefrontal cortex are altered in schizophrenia.

Twelve-hour (12 h) ultradian rhythms are a well-known phenomenon in coastal marine organisms. While 12 h cycles are observed in human behavior and physiology, no study has measured 12 h rhythms in the human brain. Here, we identify 12 h rhythms in transcripts that either peak at sleep/wake transitions (approximately 9 AM/PM) or static times (approximately 3 PM/AM) in the dorsolateral prefrontal cortex, a region involved in cognition. Subjects with schizophrenia (SZ) lose 12 h rhythms in genes associated with the unfolded protein response and neuronal structural maintenance. Moreover, genes involved in mitochondrial function and protein translation, which normally peak at sleep/wake transitions, peak instead at static times in SZ, suggesting suboptimal timing of these essential processes.

Investigating trait variability of gene co-expression network architecture in brain by controlling for genomic risk of schizophrenia.

The effect of schizophrenia (SCZ) genetic risk on gene expression in brain remains elusive. A popular approach to this problem has been the application of gene co-expression network algorithms (e.g., WGCNA). To improve reliability with this method it is critical to remove unwanted sources of variance while also preserving biological signals of interest. In this WCGNA study of RNA-Seq data from postmortem prefrontal cortex (78 neurotypical donors, EUR ancestry), we tested the effects of SCZ genetic risk on co-expression networks. Specifically, we implemented a novel design in which gene expression was adjusted by linear regression models to preserve or remove variance explained by biological signal of interest (GWAS genomic scores for SCZ risk-(GS-SCZ), and genomic scores- GS of height (GS-Ht) as a negative control), while removing variance explained by covariates of non-interest. We calculated co-expression networks from adjusted expression (GS-SCZ and GS-Ht preserved or removed), and consensus between them (representative of a "background" network free of genomic scores effects). We then tested the overlap between GS-SCZ preserved modules and background networks reasoning that modules with reduced overlap would be most affected by GS-SCZ biology. Additionally, we tested these modules for convergence of SCZ risk (i.e., enrichment in PGC3 SCZ GWAS priority genes, enrichment in SCZ risk heritability and relevant biological ontologies. Our results highlight key aspects of GS-SCZ effects on brain co-expression networks, specifically: 1) preserving/removing SCZ genetic risk alters the co-expression modules; 2) biological pathways enriched in modules affected by GS-SCZ implicate processes of transcription, translation and metabolism that converge to influence synaptic transmission; 3) priority PGC3 SCZ GWAS genes and SCZ risk heritability are enriched in modules associated with GS-SCZ effects. Overall, our results indicate that gene co-expression networks that selectively integrate information about genetic risk can reveal novel combinations of biological pathways involved in schizophrenia.

Changes in patterns of age-related network connectivity are associated with risk for schizophrenia.

Alterations in fMRI-based brain functional network connectivity (FNC) are associated with schizophrenia (SCZ) and the genetic risk or subthreshold clinical symptoms preceding the onset of SCZ, which often occurs in early adulthood. Thus, age-sensitive FNC changes may be relevant to SCZ risk-related FNC. We used independent component analysis to estimate FNC from childhood to adulthood in 9,236 individuals. To capture individual brain features more accurately than single-session fMRI, we studied an average of three fMRI scans per individual. To identify potential familial risk-related FNC changes, we compared age-related FNC in first-degree relatives of SCZ patients mostly including unaffected siblings (SIB) with neurotypical controls (NC) at the same age stage. Then, we examined how polygenic risk scores for SCZ influenced risk-related FNC patterns. Finally, we investigated the same risk-related FNC patterns in adult SCZ patients (oSCZ) and young individuals with subclinical psychotic symptoms (PSY). Age-sensitive risk-related FNC patterns emerge during adolescence and early adulthood, but not before. Young SIB always followed older NC patterns, with decreased FNC in a cerebellar-occipitoparietal circuit and increased FNC in two prefrontal-sensorimotor circuits when compared to young NC. Two of these FNC alterations were also found in oSCZ, with one exhibiting reversed pattern. All were linked to polygenic risk for SCZ in unrelated individuals (R2 varied from 0.02 to 0.05). Young PSY showed FNC alterations in the same direction as SIB when compared to NC. These results suggest that age-related neurotypical FNC correlates with genetic risk for SCZ and is detectable with MRI in young participants.

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.