Highlights d Dorsolateral prefrontal cortex and hippocampus gene expression across development d Novel region-specific schizophrenia genetic risk features d Decreased regional functional coherence in schizophrenia d Public brain gene expression and eQTL resource at http:// eqtl.brainseq.org/phase2

Transcriptome analysis of human brain provides fundamental insight about development and disease, but largely relies on existing annotation. We sequenced transcriptomes of 72 prefrontal cortex samples across six life stages, and identified 50,650 differentially expression regions (DERs) associated with developmental and aging, agnostic of annotation. While many DERs annotated to non-exonic sequence (41.1%), most were similarly regulated in cytosolic mRNA extracted from independent samples. The DERs were developmentally conserved across 16 brain regions and within the developing mouse cortex, and were expressed in diverse cell and tissue types. The DERs were further enriched for active chromatin marks and clinical risk for neurodevelopmental disorders like schizophrenia. Lastly, we demonstrate quantitatively that these DERs associate with a changing neuronal phenotype related to differentiation and maturation. These data highlight conserved molecular signatures of transcriptional dynamics across brain development, some potential clinical relevance and the incomplete annotation of the human brain transcriptome.

We characterize the landscape of somatic mutations—mutations occurring after fertilization—in the human brain using ultra-deep (~250X) whole-genome sequencing of prefrontal cortex from 59 autism spectrum disorder (ASD) cases and 15 controls. We observe a mean of 26 somatic single nucleotide variants (sSNVs) per brain present in ≥4% of cells, with enrichment of mutations in coding and putative regulatory regions. Our analysis reveals that the first cell division after fertilization produces ~3.4 mutations, followed by 2–3 mutations in subsequent generations. This suggests that a typical individual possesses ~80 sSNVs present in ≥2% of cells—comparable to the number of de novo germline mutations per generation—with about half of individuals having at least one potentially function-altering somatic mutation somewhere in the cortex. ASD brains show an excess of somatic mutations in neural enhancer sequences compared to controls, suggesting that mosaic enhancer mutations may contribute to ASD risk.

Genetic variations and adverse environmental events in utero or shortly after birth can lead to abnormal brain development and increased risk of schizophrenia. γ-Aminobutyric acid (GABA), the major inhibitory neurotransmitter in the mammalian brain, plays a vital role in normal brain development. GABA synthesis is controlled by enzymes derived from two glutamic acid decarboxylase (GAD) genes, GAD1 and GAD2, both of which produce transcript isoforms. While the full-length GAD1 transcript (GAD67) has been implicated in the neuropathology of schizophrenia, the transcript structure of GAD1 in the human brain has not been fully characterized. In this study, with the use of RNA sequencing and PCR technologies, we report the discovery of 10 novel transcripts of GAD1 in the human brain. Expression levels of four novel GAD1 transcripts (8A, 8B, I80 and I86) showed a lifespan trajectory expression pattern that is anticorrelated with the expression of the full-length GAD1 transcript. In addition, methylation levels of two CpG loci within the putative GAD1 promoter were significantly associated with the schizophrenia-risk SNP rs3749034 and with the expression of GAD25 in dorsolateral prefrontal cortex (DLPFC). Moreover, schizophrenia patients who had completed suicide and/or were positive for nicotine exposure had significantly higher full-length GAD1 expression in the DLPFC. Alternative splicing of GAD1 and epigenetic state appear to play roles in the developmental profile of GAD1 expression and may contribute to GABA dysfunction in the PFC and hippocampus of patients with schizophrenia.