Search

The function of some genetic variants associated with brain-relevant traits has been explained through colocalization with expression quantitative trait loci (eQTL) conducted in bulk postmortem adult brain tissue. However, many brain-trait associated loci have unknown cellular or molecular function. These genetic variants may exert context-specific function on different molecular phenotypes including post-transcriptional changes. Here, we identified genetic regulation of RNA editing and alternative polyadenylation (APA) within a cell-type-specific population of human neural progenitors and neurons. More RNA editing and isoforms utilizing longer polyadenylation sequences were observed in neurons, likely due to higher expression of genes encoding the proteins mediating these post-transcriptional events. We also detected hundreds of cell-type-specific editing quantitative trait loci (edQTLs) and alternative polyadenylation QTLs (apaQTLs). We found colocalizations of a neuron edQTL in CCDC88A with educational attainment and a progenitor apaQTL in EP300 with schizophrenia, suggesting that genetically mediated post-transcriptional regulation during brain development leads to differences in brain function., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2024 American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

RNA splicing is highly prevalent in the brain and has strong links to neuropsychiatric disorders; yet, the role of cell type-specific splicing and transcript-isoform diversity during human brain development has not been systematically investigated. In this work, we leveraged single-molecule long-read sequencing to deeply profile the full-length transcriptome of the germinal zone and cortical plate regions of the developing human neocortex at tissue and single-cell resolution. We identified 214,516 distinct isoforms, of which 72.6% were novel (not previously annotated in Gencode version 33), and uncovered a substantial contribution of transcript-isoform diversity-regulated by RNA binding proteins-in defining cellular identity in the developing neocortex. We leveraged this comprehensive isoform-centric gene annotation to reprioritize thousands of rare de novo risk variants and elucidate genetic risk mechanisms for neuropsychiatric disorders.

Most genetic variants associated with psychiatric disorders are located in noncoding regions of the genome. To investigate their functional implications, we integrate epigenetic data from the PsychENCODE Consortium and other published sources to construct a comprehensive atlas of candidate brain cis-regulatory elements. Using deep learning, we model these elements' sequence syntax and predict how binding sites for lineage-specific transcription factors contribute to cell type-specific gene regulation in various types of glia and neurons. The elements' evolutionary history suggests that new regulatory information in the brain emerges primarily via smaller sequence mutations within conserved mammalian elements rather than entirely new human- or primate-specific sequences. However, primate-specific candidate elements, particularly those active during fetal brain development and in excitatory neurons and astrocytes, are implicated in the heritability of brain-related human traits. Additionally, we introduce PsychSCREEN, a web-based platform offering interactive visualization of PsychENCODE-generated genetic and epigenetic data from diverse brain cell types in individuals with psychiatric disorders and healthy controls.

Neuropsychiatric genome-wide association studies (GWASs), including those for autism spectrum disorder and schizophrenia, show strong enrichment for regulatory elements in the developing brain. However, prioritizing risk genes and mechanisms is challenging without a unified regulatory atlas. Across 672 diverse developing human brains, we identified 15,752 genes harboring gene, isoform, and/or splicing quantitative trait loci, mapping 3739 to cellular contexts. Gene expression heritability drops during development, likely reflecting both increasing cellular heterogeneity and the intrinsic properties of neuronal maturation. Isoform-level regulation, particularly in the second trimester, mediated the largest proportion of GWAS heritability. Through colocalization, we prioritized mechanisms for about 60% of GWAS loci across five disorders, exceeding adult brain findings. Finally, we contextualized results within gene and isoform coexpression networks, revealing the comprehensive landscape of transcriptome regulation in development and disease.

Few neuropsychiatric disorders have replicable biomarkers, prompting high-resolution and large-scale molecular studies. However, we still lack consensus on a more foundational question: whether quantitative shifts in cell types-the functional unit of life-contribute to neuropsychiatric disorders. Leveraging advances in human brain single-cell methylomics, we deconvolve seven major cell types using bulk DNA methylation profiling across 1270 postmortem brains, including from individuals diagnosed with Alzheimer's disease, schizophrenia, and autism. We observe and replicate cell-type compositional shifts for Alzheimer's disease (endothelial cell loss), autism (increased microglia), and schizophrenia (decreased oligodendrocytes), and find age- and sex-related changes. Multiple layers of evidence indicate that endothelial cell loss contributes to Alzheimer's disease, with comparable effect size to APOE genotype among older people. Genome-wide association identified five genetic loci related to cell-type composition, involving plausible genes for the neurovascular unit ( P2RX5 and TRPV3 ) and excitatory neurons ( DPY30 and MEMO1 ). These results implicate specific cell-type shifts in the pathophysiology of neuropsychiatric disorders.

Single-cell genomics is a powerful tool for studying heterogeneous tissues such as the brain. Yet little is understood about how genetic variants influence cell-level gene expression. Addressing this, we uniformly processed single-nuclei, multiomics datasets into a resource comprising >2.8 million nuclei from the prefrontal cortex across 388 individuals. For 28 cell types, we assessed population-level variation in expression and chromatin across gene families and drug targets. We identified >550,000 cell type-specific regulatory elements and >1.4 million single-cell expression quantitative trait loci, which we used to build cell-type regulatory and cell-to-cell communication networks. These networks manifest cellular changes in aging and neuropsychiatric disorders. We further constructed an integrative model accurately imputing single-cell expression and simulating perturbations; the model prioritized ~250 disease-risk genes and drug targets with associated cell types.

Methods integrating genetics with transcriptomic reference panels prioritize risk genes and mechanisms at only a fraction of trait-associated genetic loci, due in part to an overreliance on total gene expression as a molecular outcome measure. This challenge is particularly relevant for the brain, in which extensive splicing generates multiple distinct transcript-isoforms per gene. Due to complex correlation structures, isoform-level modeling from cis-window variants requires methodological innovation. Here we introduce isoTWAS, a multivariate, stepwise framework integrating genetics, isoform-level expression and phenotypic associations. Compared to gene-level methods, isoTWAS improves both isoform and gene expression prediction, yielding more testable genes, and increased power for discovery of trait associations within genome-wide association study loci across 15 neuropsychiatric traits. We illustrate multiple isoTWAS associations undetectable at the gene-level, prioritizing isoforms of AKT3, CUL3 and HSPD1 in schizophrenia and PCLO with multiple disorders. Results highlight the importance of incorporating isoform-level resolution within integrative approaches to increase discovery of trait associations, especially for brain-relevant traits., (© 2023. The Author(s).)

As recreational use of cannabis is being decriminalized in many places and medical use widely sanctioned, there are growing concerns about increases in cannabis use disorder (CanUD), which is associated with numerous medical comorbidities. Here we performed a genome-wide association study of CanUD in the Million Veteran Program (MVP), followed by meta-analysis in 1,054,365 individuals (n cases  = 64,314) from four broad ancestries designated by the reference panel used for assignment (European n = 886,025, African n = 123,208, admixed American n = 38,289 and East Asian n = 6,843). Population-specific methods were applied to calculate single nucleotide polymorphism-based heritability within each ancestry. Statistically significant single nucleotide polymorphism-based heritability for CanUD was observed in all but the smallest population (East Asian). We discovered genome-wide significant loci unique to each ancestry: 22 in European, 2 each in African and East Asian, and 1 in admixed American ancestries. A genetically informed causal relationship analysis indicated a possible effect of genetic liability for CanUD on lung cancer risk, suggesting potential unanticipated future medical and psychiatric public health consequences that require further study to disentangle from other known risk factors such as cigarette smoking., (© 2023. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)

Human brain size changes dynamically through early development, peaks in adolescence, and varies up to 2-fold among adults. However, the molecular genetic underpinnings of interindividual variation in brain size remain unknown. Here, we leveraged postmortem brain RNA sequencing and measurements of brain weight (BW) in 2,531 individuals across three independent datasets to identify 928 genome-wide significant associations with BW. Genes associated with higher or lower BW showed distinct neurodevelopmental trajectories and spatial patterns that mapped onto functional and cellular axes of brain organization. Expression of BW genes was predictive of interspecies differences in brain size, and bioinformatic annotation revealed enrichment for neurogenesis and cell-cell communication. Genome-wide, transcriptome-wide, and phenome-wide association analyses linked BW gene sets to neuroimaging measurements of brain size and brain-related clinical traits. Cumulatively, these results represent a major step toward delineating the molecular pathways underlying human brain size variation in health and disease., Competing Interests: Declaration of interests J.S., R.A.I.B., J.D.B., and A.F.A.-B. are directors and hold equity in Centile Bioscience. J.D.B. holds positions/equity in UpFront Diagnostics, Treovir, and NeuroX1. A.F.A.-B. receives consulting income from Octave Bioscience., (Copyright © 2023 The Authors. Published by Elsevier Inc. All rights reserved.)

Changes in the function of inhibitory interneurons (INs) during cortical development could contribute to the pathophysiology of neurodevelopmental disorders. Using all-optical in vivo approaches, we find that parvalbumin (PV) INs and their immature precursors are hypoactive and transiently decoupled from excitatory neurons in postnatal mouse somatosensory cortex (S1) of Fmr1 KO mice, a model of fragile X syndrome (FXS). This leads to a loss of parvalbumin INs (PV-INs) in both mice and humans with FXS. Increasing the activity of future PV-INs in neonatal Fmr1 KO mice restores PV-IN density and ameliorates transcriptional dysregulation in S1, but not circuit dysfunction. Critically, administering an allosteric modulator of Kv3.1 channels after the S1 critical period does rescue circuit dynamics and tactile defensiveness. Symptoms in FXS and related disorders could be mitigated by targeting PV-INs., Competing Interests: Declaration of interests A.G.S., C.M., and B.J.H. are employees of H. Lundbeck A/S, who developed the AG00563 compound. A patent was filled by H. Lundbeck A/S (WO 2020/089262)., (Copyright © 2023 The Authors. Published by Elsevier Inc. All rights reserved.)

Our understanding of the genetics of the human cerebral cortex is limited both in terms of the diversity and the anatomical granularity of brain structural phenotypes. Here we conducted a genome-wide association meta-analysis of 13 structural and diffusion magnetic resonance imaging-derived cortical phenotypes, measured globally and at 180 bilaterally averaged regions in 36,663 individuals and identified 4,349 experiment-wide significant loci. These phenotypes include cortical thickness, surface area, gray matter volume, measures of folding, neurite density and water diffusion. We identified four genetic latent structures and causal relationships between surface area and some measures of cortical folding. These latent structures partly relate to different underlying gene expression trajectories during development and are enriched for different cell types. We also identified differential enrichment for neurodevelopmental and constrained genes and demonstrate that common genetic variants associated with cortical expansion are associated with cephalic disorders. Finally, we identified complex interphenotype and inter-regional genetic relationships among the 13 phenotypes, reflecting the developmental differences among them. Together, these analyses identify distinct genetic organizational principles of the cortex and their correlates with neurodevelopment., (© 2023. The Author(s), under exclusive licence to Springer Nature America, Inc.)

We conducted a genome-wide association study (GWAS) analysis of diverticular disease (DivD) of intestine within 724,372 individuals and identified 150 independent genome-wide significant DNA variants. Integration of the GWAS results with human gut single-cell RNA sequencing data implicated gut myocyte, mesothelial and stromal cells, and enteric neurons and glia in DivD development. Ninety-five genes were prioritized based on multiple lines of evidence, including SLC9A3 , a drug target gene of tenapanor used for the treatment of the constipation subtype of irritable bowel syndrome. A DivD polygenic score (PGS) enables effective risk prediction (area under the curve [AUC], 0.688; 95% confidence interval [CI], 0.645-0.732) and the top 20% PGS was associated with ∼3.6-fold increased DivD risk relative to the remaining population. Our statistical and bioinformatic analyses suggest that the mechanism of DivD is through colon structure, gut motility, gastrointestinal mucus, and ionic homeostasis. Our analyses reinforce the link between gastrointestinal disorders and the enteric nervous system through genetics., Competing Interests: The authors declare no competing interests., (© 2023 The Authors.)

Background and Hypothesis: Viral infections are increasingly recognized in the etiology of psychiatric disorders based on epidemiological and serological studies. Few studies have analyzed viruses directly within the brain and no comprehensive investigation of viral infection within diseased brains has been completed. This study aims to determine whether viral infection in brain tissues is a risk factor for 3 major psychiatric disorders, including schizophrenia, bipolar disorder, and autism spectrum disorder., Study Design: This study directly evaluated the presence of viral DNA or RNA in 1569 brains of patients and controls using whole-genome sequencing and RNA sequencing data with 4 independent cohorts. The PathSeq tool was used to identify known human viruses in the genome and transcriptome of patients and controls., Study Results: A variety of DNA and RNA viruses related to the central nervous system were detected in the brains of patients with major psychiatric disorders, including viruses belonging to Herpesviridae, Polyomaviridae, Retroviridae, Flaviviridae, Parvoviridae, and Adenoviridae. However, no consistent significant differences were found between patients and controls in terms of types and amount of virus detected at both DNA and RNA levels., Conclusions: The findings of this study do not suggest an association between viral infection in postmortem brains and major psychiatric disorders., (© The Author(s) 2023. Published by Oxford University Press on behalf of the Maryland Psychiatric Research Center. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

Autism omics research has historically been reductionist and diagnosis centric, with little attention paid to common co-occurring conditions (for example, sleep and feeding disorders) and the complex interplay between molecular profiles and neurodevelopment, genetics, environmental factors and health. Here we explored the plasma lipidome (783 lipid species) in 765 children (485 diagnosed with autism spectrum disorder (ASD)) within the Australian Autism Biobank. We identified lipids associated with ASD diagnosis (n = 8), sleep disturbances (n = 20) and cognitive function (n = 8) and found that long-chain polyunsaturated fatty acids may causally contribute to sleep disturbances mediated by the FADS gene cluster. We explored the interplay of environmental factors with neurodevelopment and the lipidome, finding that sleep disturbances and unhealthy diet have a convergent lipidome profile (with potential mediation by the microbiome) that is also independently associated with poorer adaptive function. In contrast, ASD lipidome differences were accounted for by dietary differences and sleep disturbances. We identified a large chr19p13.2 copy number variant genetic deletion spanning the LDLR gene and two high-confidence ASD genes (ELAVL3 and SMARCA4) in one child with an ASD diagnosis and widespread low-density lipoprotein-related lipidome derangements. Lipidomics captures the complexity of neurodevelopment, as well as the biological effects of conditions that commonly affect quality of life among autistic people., (© 2023. The Author(s).)

Autism spectrum disorder (ASD) is a highly heterogeneous disorder, yet transcriptomic profiling of bulk brain tissue has identified substantial convergence among dysregulated genes and pathways in ASD. However, this approach lacks cell-specific resolution. We performed comprehensive transcriptomic analyses on bulk tissue and laser-capture microdissected (LCM) neurons from 59 postmortem human brains (27 ASD and 32 controls) in the superior temporal gyrus (STG) of individuals ranging from 2 to 73 years of age. In bulk tissue, synaptic signaling, heat shock protein-related pathways, and RNA splicing were significantly altered in ASD. There was age-dependent dysregulation of genes involved in gamma aminobutyric acid (GABA) ( GAD1 and GAD2 ) and glutamate ( SLC38A1 ) signaling pathways. In LCM neurons, AP-1-mediated neuroinflammation and insulin/IGF-1 signaling pathways were upregulated in ASD, while mitochondrial function, ribosome, and spliceosome components were downregulated. GABA synthesizing enzymes GAD1 and GAD2 were both downregulated in ASD neurons. Mechanistic modeling suggested a direct link between inflammation and ASD in neurons, and prioritized inflammation-associated genes for future study. Alterations in small nucleolar RNAs (snoRNAs) associated with splicing events suggested interplay between snoRNA dysregulation and splicing disruption in neurons of individuals with ASD. Our findings supported the fundamental hypothesis of altered neuronal communication in ASD, demonstrated that inflammation was elevated at least in part in ASD neurons, and may reveal windows of opportunity for biotherapeutics to target the trajectory of gene expression and clinical manifestation of ASD throughout the human lifespan.

Background: Increased expression of the complement component 4A (C4A) gene is associated with a greater lifetime risk of schizophrenia. In the brain, C4A is involved in synaptic pruning; yet, it remains unclear the extent to which upregulation of C4A alters brain development or is associated with the risk for psychotic symptoms in childhood. Here, we perform a multi-ancestry phenome-wide association study in 7789 children aged 9-12 years to examine the relationship between genetically regulated expression (GREx) of C4A, childhood brain structure, cognition, and psychiatric symptoms., Results: While C4A GREx is not related to childhood psychotic experiences, cognition, or global measures of brain structure, it is associated with a localized reduction in regional surface area (SA) of the entorhinal cortex. Furthermore, we show that reduced entorhinal cortex SA at 9-10 years predicts a greater number and severity of psychosis-like events at 1-year and 2-year follow-up time points. We also demonstrate that the effects of C4A on the entorhinal cortex are independent of genome-wide polygenic risk for schizophrenia., Conclusions: Our results suggest neurodevelopmental effects of C4A on childhood medial temporal lobe structure, which may serve as a biomarker for schizophrenia risk prior to symptom onset., (© 2023. The Author(s).)

Sex differences in white matter microstructure have been robustly demonstrated in the adult brain using both conventional and advanced diffusion-weighted magnetic resonance imaging approaches. However, sex differences in white matter microstructure prior to adulthood remain poorly understood; previous developmental work focused on conventional microstructure metrics and yielded mixed results. Here, we rigorously characterized sex differences in white matter microstructure among over 6000 children from the Adolescent Brain Cognitive Development study who were between 9 and 10 years old. Microstructure was quantified using both the conventional model-diffusion tensor imaging (DTI)-and an advanced model, restriction spectrum imaging (RSI). DTI metrics included fractional anisotropy (FA) and mean, axial, and radial diffusivity (MD, AD, RD). RSI metrics included normalized isotropic, directional, and total intracellular diffusion (N0, ND, NT). We found significant and replicable sex differences in DTI or RSI microstructure metrics in every white matter region examined across the brain. Sex differences in FA were regionally specific. Across white matter regions, boys exhibited greater MD, AD, and RD than girls, on average. Girls displayed increased N0, ND, and NT compared to boys, on average, suggesting greater cell and neurite density in girls. Together, these robust and replicable findings provide an important foundation for understanding sex differences in health and disease., (© 2022 The Authors. Human Brain Mapping published by Wiley Periodicals LLC.)

Summary: With the continued deluge of results from genome-wide association and functional genomic studies, it has become increasingly imperative to quickly combine and visualize different layers of genetic and genomic data within a given locus to facilitate exploratory and integrative data analyses. While several tools have been developed to visualize locus-level genetic results, the limited speed, scalability and flexibility of current approaches remain a significant bottleneck. Here, we present a Julia package for high-performance genetics and genomics-related data visualization that enables fast, simultaneous plotting of hundreds of association results along with multiple relevant genomic annotations. Leveraging the powerful plotting and layout utilities from Makie.jl facilitates the customization and extensibility of every component of a plot, enabling generation of publication-ready figures., Availability and Implementation: The GeneticsMakie.jl package is open source and distributed under the MIT license via GitHub (https://github.com/mmkim1210/GeneticsMakie.jl). The GitHub repository contains installation instructions as well as examples and documentation for built-in functions., Supplementary Information: Supplementary data are available at Bioinformatics online., (© The Author(s) 2022. Published by Oxford University Press.)

Neuropsychiatric disorders classically lack defining brain pathologies, but recent work has demonstrated dysregulation at the molecular level, characterized by transcriptomic and epigenetic alterations 1-3 . In autism spectrum disorder (ASD), this molecular pathology involves the upregulation of microglial, astrocyte and neural-immune genes, the downregulation of synaptic genes, and attenuation of gene-expression gradients in cortex 1,2,4-6 . However, whether these changes are limited to cortical association regions or are more widespread remains unknown. To address this issue, we performed RNA-sequencing analysis of 725 brain samples spanning 11 cortical areas from 112 post-mortem samples from individuals with ASD and neurotypical controls. We find widespread transcriptomic changes across the cortex in ASD, exhibiting an anterior-to-posterior gradient, with the greatest differences in primary visual cortex, coincident with an attenuation of the typical transcriptomic differences between cortical regions. Single-nucleus RNA-sequencing and methylation profiling demonstrate that this robust molecular signature reflects changes in cell-type-specific gene expression, particularly affecting excitatory neurons and glia. Both rare and common ASD-associated genetic variation converge within a downregulated co-expression module involving synaptic signalling, and common variation alone is enriched within a module of upregulated protein chaperone genes. These results highlight widespread molecular changes across the cerebral cortex in ASD, extending beyond association cortex to broadly involve primary sensory regions., (© 2022. The Author(s).)

Telencephalic organoids generated from human pluripotent stem cells (hPSCs) are a promising system for studying the distinct features of the developing human brain and the underlying causes of many neurological disorders. While organoid technology is steadily advancing, many challenges remain, including potential batch-to-batch and cell-line-to-cell-line variability, and structural inconsistency. Here, we demonstrate that a major contributor to cortical organoid quality is the way hPSCs are maintained prior to differentiation. Optimal results were achieved using particular fibroblast-feeder-supported hPSCs rather than feeder-independent cells, differences that were reflected in their transcriptomic states at the outset. Feeder-supported hPSCs displayed activation of diverse transforming growth factor β (TGFβ) superfamily signaling pathways and increased expression of genes connected to naive pluripotency. We further identified combinations of TGFβ-related growth factors that are necessary and together sufficient to impart broad telencephalic organoid competency to feeder-free hPSCs and enhance the formation of well-structured brain tissues suitable for disease modeling., Competing Interests: Conflicts of interest The authors declare no competing interests., (Copyright © 2022 The Author(s). Published by Elsevier Inc. All rights reserved.)

Gene expression covaries with brain activity as measured by resting state functional magnetic resonance imaging (MRI). However, it is unclear how genomic differences driven by disease state can affect this relationship. Here, we integrate from the ABIDE I and II imaging cohorts with datasets of gene expression in brains of neurotypical individuals and individuals with autism spectrum disorder (ASD) with regionally matched brain activity measurements from fMRI datasets. We identify genes linked with brain activity whose association is disrupted in ASD. We identified a subset of genes that showed a differential developmental trajectory in individuals with ASD compared with controls. These genes are enriched in voltage-gated ion channels and inhibitory neurons, pointing to excitation-inhibition imbalance in ASD. We further assessed differences at the regional level showing that the primary visual cortex is the most affected region in ASD. Our results link disrupted brain expression patterns of individuals with ASD to brain activity and show developmental, cell type, and regional enrichment of activity linked genes., (© 2022. The Author(s).)

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., (© 2022. The Author(s), under exclusive licence to Springer Nature Limited.)

Astrocytes are critical for the development and function of synapses. There are notable species differences between human astrocytes and commonly used animal models. Yet, it is unclear whether astrocytic genes involved in synaptic function are stable or exhibit dynamic changes associated with disease states and age in humans, which is a barrier in understanding human astrocyte biology and its potential involvement in neurologic diseases. To better understand the properties of human astrocytes, we acutely purified astrocytes from the cerebral cortices of over 40 humans across various ages, sexes, and disease states. We performed RNA sequencing to generate transcriptomic profiles of these astrocytes and identified genes associated with these biological variables. We found that human astrocytes in tumor-surrounding regions downregulate genes involved in synaptic function and sensing of signals in the microenvironment, suggesting involvement of peritumor astrocytes in tumor-associated neural circuit dysfunction. In aging, we also found downregulation of synaptic regulators and upregulation of markers of cytokine signaling, while in maturation we identified changes in ionic transport with implications for calcium signaling. In addition, we identified subtle sexual dimorphism in human cortical astrocytes, which has implications for observed sex differences across many neurologic disorders. Overall, genes involved in synaptic function exhibit dynamic changes in the peritumor microenvironment and aging. These data provide powerful new insights into human astrocyte biology in several biologically relevant states that will aid in generating novel testable hypotheses about homeostatic and reactive astrocytes in humans. SIGNIFICANCE STATEMENT Astrocytes are an abundant class of cells playing integral roles at synapses. Astrocyte dysfunction is implicated in a variety of human neurologic diseases. Yet our knowledge of astrocytes is largely based on mouse studies. Direct knowledge of human astrocyte biology remains limited. Here, we present transcriptomic profiles of human cortical astrocytes, and we identified molecular differences associated with age, sex, and disease state. We found that peritumor and aging astrocytes downregulate genes involved in astrocyte-synapse interactions. These data provide necessary insight into human astrocyte biology that will improve our understanding of human disease., (Copyright © 2022 the authors.)

Background: Childhood sleep problems are common and among the most frequent and impairing comorbidities of childhood psychiatric disorders. In adults, sleep disturbances are heritable and show strong genetic associations with brain morphology; however, little is known about the genetic architecture of childhood sleep and potential etiological links between sleep, brain development, and pediatric-onset psychiatric symptoms., Methods: Using data from the Adolescent Brain Cognitive Development Study ( n Phenotype = 4428 for discovery/replication, n Genetics  = 4728; age 9-10 years), we assessed phenotypic relationships, heritability, and genetic correlations between childhood sleep disturbances (insomnia, arousal, breathing, somnolence, hyperhidrosis, sleep-wake transitions), brain size (surface area, cortical thickness, volume), and dimensional psychopathology., Results: Sleep disturbances showed widespread positive associations with multiple domains of childhood psychopathology; however, only insomnia showed replicable associations with smaller brain surface area. Among the sleep disturbances assessed, only insomnia showed significant heritability ( h 2 SNP  = 0.15, p < .05) and showed substantial genetic correlations with externalizing and attention-deficit/hyperactivity disorder symptomatology ( r G s > 0.80, p s < .05). We found no evidence of genetic correlation between childhood insomnia and brain size. Furthermore, polygenic risk scores calculated from genome-wide association studies of adult insomnia and adult brain size did not predict childhood insomnia; instead, polygenic risk scores trained using attention-deficit/hyperactivity disorder genome-wide association studies predicted decreased surface area at baseline as well as insomnia and externalizing symptoms longitudinally., Conclusions: Findings demonstrate a distinct genetic architecture underlying childhood insomnia and brain size and suggest genetic overlap between childhood insomnia and attention-deficit/hyperactivity disorder symptomatology. Additional research is needed to examine how genetic risk manifests in altered developmental trajectories and comorbid sleep/psychiatric symptoms across adolescence., (© 2022 The Authors.)

Human epidemiological studies implicate exposure to infection during gestation in the etiology of neurodevelopmental disorders. Animal models of maternal immune activation (MIA) have identified the maternal immune response as the critical link between maternal infection and aberrant offspring brain and behavior development. Here we evaluate neurodevelopment of male rhesus monkeys ( Macaca mulatta ) born to MIA-treated dams ( n = 14) injected with a modified form of the viral mimic polyinosinic:polycytidylic acid at the end of the first trimester. Control dams received saline injections at the same gestational time points ( n = 10) or were untreated ( n = 4). MIA-treated dams exhibited a strong immune response as indexed by transient increases in sickness behavior, temperature, and inflammatory cytokines. Although offspring born to control or MIA-treated dams did not differ on measures of physical growth and early developmental milestones, the MIA-treated animals exhibited subtle changes in cognitive development and deviated from species-typical brain growth trajectories. Longitudinal MRI revealed significant gray matter volume reductions in the prefrontal and frontal cortices of MIA-treated offspring at 6 months that persisted through the final time point at 45 months along with smaller frontal white matter volumes in MIA-treated animals at 36 and 45 months. These findings provide the first evidence of early postnatal changes in brain development in MIA-exposed nonhuman primates and establish a translationally relevant model system to explore the neurodevelopmental trajectory of risk associated with prenatal immune challenge from birth through late adolescence. SIGNIFICANCE STATEMENT Women exposed to infection during pregnancy have an increased risk of giving birth to a child who will later be diagnosed with a neurodevelopmental disorder. Preclinical maternal immune activation (MIA) models have demonstrated that the effects of maternal infection on fetal brain development are mediated by maternal immune response. Since the majority of MIA models are conducted in rodents, the nonhuman primate provides a unique system to evaluate the MIA hypothesis in a species closely related to humans. Here we report the first longitudinal study conducted in a nonhuman primate MIA model. MIA-exposed offspring demonstrate subtle changes in cognitive development paired with marked reductions in frontal gray and white matter, further supporting the association between prenatal immune challenge and alterations in offspring neurodevelopment., (Copyright © 2021 the authors.)

Alternative splicing is a post-transcriptional regulatory mechanism producing distinct mRNA molecules from a single pre-mRNA with a prominent role in the development and function of the central nervous system. We used long-read isoform sequencing to generate full-length transcript sequences in the human and mouse cortex. We identify novel transcripts not present in existing genome annotations, including transcripts mapping to putative novel (unannotated) genes and fusion transcripts incorporating exons from multiple genes. Global patterns of transcript diversity are similar between human and mouse cortex, although certain genes are characterized by striking differences between species. We also identify developmental changes in alternative splicing, with differential transcript usage between human fetal and adult cortex. Our data confirm the importance of alternative splicing in the cortex, dramatically increasing transcriptional diversity and representing an important mechanism underpinning gene regulation in the brain. We provide transcript-level data for human and mouse cortex as a resource to the scientific community., Competing Interests: Declaration of interests Z.A. and D.A.C. were full-time employees of Eli Lilly & Company, Ltd., and E.T. was a full-time employee of PacBio at the time this work was performed. All other authors declare no competing interests., (Copyright © 2021 The Author(s). Published by Elsevier Inc. All rights reserved.)

Schizophrenia is a heterogeneous disorder, exhibiting variability in presentation and outcomes that complicate treatment and recovery. To explore this heterogeneity, we leverage the comprehensive Danish health registries to conduct a prospective, longitudinal study from birth of 5432 individuals who would ultimately be diagnosed with schizophrenia, building individual trajectories that represent sequences of comorbid diagnoses, and describing patterns in the individual-level variability. We show that psychiatric comorbidity is prevalent among individuals with schizophrenia (82%) and multi-morbidity occur more frequently in specific, time-ordered pairs. Three latent factors capture 79% of variation in longitudinal comorbidity and broadly relate to the number of co-occurring diagnoses, the presence of child versus adult comorbidities and substance abuse. Clustering of the factor scores revealed five stable clusters of individuals, associated with specific risk factors and outcomes. The presentation and course of schizophrenia may be associated with heterogeneity in etiological factors including family history of mental disorders., (© 2021. The Author(s).)

Brain organoids represent a powerful tool for studying human neurological diseases, particularly those that affect brain growth and structure. However, many diseases manifest with clear evidence of physiological and network abnormality in the absence of anatomical changes, raising the question of whether organoids possess sufficient neural network complexity to model these conditions. Here, we explore the network-level functions of brain organoids using calcium sensor imaging and extracellular recording approaches that together reveal the existence of complex network dynamics reminiscent of intact brain preparations. We demonstrate highly abnormal and epileptiform-like activity in organoids derived from induced pluripotent stem cells from individuals with Rett syndrome, accompanied by transcriptomic differences revealed by single-cell analyses. We also rescue key physiological activities with an unconventional neuroregulatory drug, pifithrin-α. Together, these findings provide an essential foundation for the utilization of brain organoids to study intact and disordered human brain network formation and illustrate their utility in therapeutic discovery., (© 2021. The Author(s), under exclusive licence to Springer Nature America, Inc.)

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.

Bipolar disorder is a heritable mental illness with complex etiology. We performed a genome-wide association study of 41,917 bipolar disorder cases and 371,549 controls of European ancestry, which identified 64 associated genomic loci. Bipolar disorder risk alleles were enriched in genes in synaptic signaling pathways and brain-expressed genes, particularly those with high specificity of expression in neurons of the prefrontal cortex and hippocampus. Significant signal enrichment was found in genes encoding targets of antipsychotics, calcium channel blockers, antiepileptics and anesthetics. Integrating expression quantitative trait locus data implicated 15 genes robustly linked to bipolar disorder via gene expression, encoding druggable targets such as HTR6, MCHR1, DCLK3 and FURIN. Analyses of bipolar disorder subtypes indicated high but imperfect genetic correlation between bipolar disorder type I and II and identified additional associated loci. Together, these results advance our understanding of the biological etiology of bipolar disorder, identify novel therapeutic leads and prioritize genes for functional follow-up studies.

Background: Maternal immune activation (MIA) is a proposed risk factor for multiple neuropsychiatric disorders, including schizophrenia. However, the molecular mechanisms through which MIA imparts risk remain poorly understood. A recently developed nonhuman primate model of exposure to the viral mimic poly:ICLC during pregnancy shows abnormal social and repetitive behaviors and elevated striatal dopamine, a molecular hallmark of human psychosis, providing an unprecedented opportunity for studying underlying molecular correlates., Methods: We performed RNA sequencing across psychiatrically relevant brain regions (prefrontal cortex, anterior cingulate, hippocampus) and primary visual cortex for comparison from 3.5- to 4-year-old male MIA-exposed and control offspring-an age comparable to mid adolescence in humans., Results: We identify 266 unique genes differentially expressed in at least one brain region, with the greatest number observed in hippocampus. Co-expression networks identified region-specific alterations in synaptic signaling and oligodendrocytes. Although we observed temporal and regional differences, transcriptomic changes were shared across first- and second-trimester exposures, including for the top differentially expressed genes-PIWIL2 and MGARP. In addition to PIWIL2, several other regulators of retrotransposition and endogenous transposable elements were dysregulated following MIA, potentially connecting MIA to retrotransposition., Conclusions: Together, these results begin to elucidate the brain-level molecular processes through which MIA may impart risk for psychiatric disease., (Copyright © 2020 Society of Biological Psychiatry. Published by Elsevier Inc. All rights reserved.)

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., (Published by Elsevier Inc.)

Autism spectrum disorder (ASD) is a phenotypically and genetically heterogeneous neurodevelopmental disorder. Despite this heterogeneity, previous studies have shown patterns of molecular convergence in post-mortem brain tissue from autistic subjects. Here, we integrate genome-wide measures of mRNA expression, miRNA expression, DNA methylation, and histone acetylation from ASD and control brains to identify a convergent molecular subtype of ASD with shared dysregulation across both the epigenome and transcriptome. Focusing on this convergent subtype, we substantially expand the repertoire of differentially expressed genes in ASD and identify a component of upregulated immune processes that are associated with hypomethylation. We utilize eQTL and chromosome conformation datasets to link differentially acetylated regions with their cognate genes and identify an enrichment of ASD genetic risk variants in hyperacetylated noncoding regulatory regions linked to neuronal genes. These findings help elucidate how diverse genetic risk factors converge onto specific molecular processes in ASD.

Background: 22q11.2 copy number variants are among the most highly penetrant genetic risk variants for developmental neuropsychiatric disorders such as schizophrenia (SCZ) and autism spectrum disorder (ASD). However, the specific mechanisms through which they confer risk remain unclear., Methods: Using a functional genomics approach, we integrated transcriptomic data from the developing human brain, genome-wide association findings for SCZ and ASD, protein interaction data, and gene expression signatures from SCZ and ASD postmortem cortex to 1) organize genes into the developmental cellular and molecular systems within which they operate, 2) identify neurodevelopmental processes associated with polygenic risk for SCZ and ASD across the allelic frequency spectrum, and 3) elucidate pathways and individual genes through which 22q11.2 copy number variants may confer risk for each disorder., Results: Polygenic risk for SCZ and ASD converged on partially overlapping neurodevelopmental modules involved in synaptic function and transcriptional regulation, with ASD risk variants additionally enriched for modules involved in neuronal differentiation during fetal development. The 22q11.2 locus formed a large protein network during development that disproportionately affected SCZ-associated and ASD-associated neurodevelopmental modules, including loading highly onto synaptic and gene regulatory pathways. SEPT5, PI4KA, and SNAP29 genes are candidate drivers of 22q11.2 synaptic pathology relevant to SCZ and ASD, and DGCR8 and HIRA are candidate drivers of disease-relevant alterations in gene regulation., Conclusions: This approach offers a powerful framework to identify neurodevelopmental processes affected by diverse risk variants for SCZ and ASD and elucidate mechanisms through which highly penetrant, multigene copy number variants contribute to disease risk., (Copyright © 2019 Society of Biological Psychiatry. Published by Elsevier Inc. All rights reserved.)

Gene-environment interactions (GxE) can be fundamental in applications ranging from functional genomics to precision medicine and is a conjectured source of substantial heritability. However, unbiased methods to profile GxE genome-wide are nascent and, as we show, cannot accommodate general environment variables, modest sample sizes, heterogeneous noise, and binary traits. To address this gap, we propose a simple, unifying mixed model for gene-environment interaction (GxEMM). In simulations and theory, we show that GxEMM can dramatically improve estimates and eliminate false positives when the assumptions of existing methods fail. We apply GxEMM to a range of human and model organism datasets and find broad evidence of context-specific genetic effects, including GxSex, GxAdversity, and GxDisease interactions across thousands of clinical and molecular phenotypes. Overall, GxEMM is broadly applicable for testing and quantifying polygenic interactions, which can be useful for explaining heritability and invaluable for determining biologically relevant environments., (Copyright © 2019 American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

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., (Copyright © 2019 Elsevier Inc. All rights reserved.)

Many cellular and physiological processes are coordinated by regulatory networks that produce a remarkable complexity of transcript isoforms. In the mammalian nervous system, alternative pre-mRNA splicing generates functionally distinct isoforms that play key roles in normal physiology, supporting development, plasticity, complex behaviors, and cognition. Neuronal splicing programs controlled by RNA-binding proteins, are influenced by chromatin modifications and can exhibit neuronal subtype specificity. As highlighted in recent publications, aberrant alternative splicing is a major contributor to disease phenotypes. Therefore, understanding the underlying mechanisms of alternative splicing regulation and identifying functional splicing isoforms with critical phenotypic roles are expected to provide a comprehensive resource for therapeutic development, as illuminated by recent successful interventions of spinal muscular atrophy. Here, we discuss the latest progress in the study of the emerging complexity of alternative splicing mechanisms in neurons, and how these findings inform new therapies to correct and control splicing defects., (Copyright © 2019 the authors.)

One fundamental but understudied mechanism of gene regulation in disease is allele-specific expression (ASE), the preferential expression of one allele. We leveraged RNA-sequencing data from human brain to assess ASE in autism spectrum disorder (ASD). When ASE is observed in ASD, the allele with lower population frequency (minor allele) is preferentially more highly expressed than the major allele, opposite to the canonical pattern. Importantly, genes showing ASE in ASD are enriched in those downregulated in ASD postmortem brains and in genes harboring de novo mutations in ASD. Two regions, 14q32 and 15q11, containing all known orphan C/D box small nucleolar RNAs (snoRNAs), are particularly enriched in shifts to higher minor allele expression. We demonstrate that this allele shifting enhances snoRNA-targeted splicing changes in ASD-related target genes in idiopathic ASD and 15q11-q13 duplication syndrome. Together, these results implicate allelic imbalance and dysregulation of orphan C/D box snoRNAs in ASD pathogenesis.

We performed a comprehensive assessment of rare inherited variation in autism spectrum disorder (ASD) by analyzing whole-genome sequences of 2,308 individuals from families with multiple affected children. We implicate 69 genes in ASD risk, including 24 passing genome-wide Bonferroni correction and 16 new ASD risk genes, most supported by rare inherited variants, a substantial extension of previous findings. Biological pathways enriched for genes harboring inherited variants represent cytoskeletal organization and ion transport, which are distinct from pathways implicated in previous studies. Nevertheless, the de novo and inherited genes contribute to a common protein-protein interaction network. We also identified structural variants (SVs) affecting non-coding regions, implicating recurrent deletions in the promoters of DLG2 and NR3C2. Loss of nr3c2 function in zebrafish disrupts sleep and social function, overlapping with human ASD-related phenotypes. These data support the utility of studying multiplex families in ASD and are available through the Hartwell Autism Research and Technology portal., (Copyright © 2019 Elsevier Inc. All rights reserved.)

Autism spectrum disorder (ASD) manifests as alterations in complex human behaviors including social communication and stereotypies. In addition to genetic risks, the gut microbiome differs between typically developing (TD) and ASD individuals, though it remains unclear whether the microbiome contributes to symptoms. We transplanted gut microbiota from human donors with ASD or TD controls into germ-free mice and reveal that colonization with ASD microbiota is sufficient to induce hallmark autistic behaviors. The brains of mice colonized with ASD microbiota display alternative splicing of ASD-relevant genes. Microbiome and metabolome profiles of mice harboring human microbiota predict that specific bacterial taxa and their metabolites modulate ASD behaviors. Indeed, treatment of an ASD mouse model with candidate microbial metabolites improves behavioral abnormalities and modulates neuronal excitability in the brain. We propose that the gut microbiota regulates behaviors in mice via production of neuroactive metabolites, suggesting that gut-brain connections contribute to the pathophysiology of ASD., (Copyright © 2019 Elsevier Inc. All rights reserved.)

Background: Williams syndrome (WS) is a neurodevelopmental disorder that has been attributed to heterozygous deletions in chromosome 7q11.23 and exhibits a variety of physical, cognitive, and behavioral features. However, the genetic basis of this phenotypic variability is unclear. In this study, we identified genetic clues underlying these complex phenotypes., Methods: Neurobehavioral function was assessed in WS patients and healthy controls. Total RNA was extracted from peripheral blood and subjected to microarray analysis, RNA-sequencing, and qRT-PCR. Weighted gene co-expression network analysis was performed to identify specific alterations related to intermediate disease phenotypes. To functionally interpret each WS-related module, gene ontology and disease-related gene enrichment were examined. We also investigated the micro (mi)RNA expression profiles and miRNA co-expression networks to better explain the regulation of the transcriptome in WS., Results: Our analysis identified four significant co-expression modules related to intermediate WS phenotypes. Notably, the three upregulated WS-related modules were composed exclusively of genes located outside the 7q11.23 region. They were significantly enriched in genes related to B-cell activation, RNA processing, and RNA transport. BCL11A, which is known for its association with speech disorders and intellectual disabilities, was identified as one of the hub genes in the top WS-related module. Finally, these key upregulated mRNA co-expression modules appear to be inversely correlated with a specific downregulated WS-related miRNA co-expression module., Conclusions: Dysregulation of the mRNA/miRNA network involving genes outside of the 7q11.23 region is likely related to the complex phenotypes observed in WS patients., (© 2018 The Authors. Journal of Child Psychology and Psychiatry published by John Wiley & Sons Ltd on behalf of Association for Child and Adolescent Mental Health.)