Mapping microRNA expression quantitative trait loci in the prenatal human brain implicates miR-1908-5p expression in bipolar disorder and other brain-related traits.

MicroRNA (miRNA) are small non-coding RNA involved in post-transcriptional gene regulation. Given their known involvement in early neurodevelopment processes, we here sought to identify common genetic variants associated with altered miRNA expression in the prenatal human brain. We performed small RNA sequencing on brain tissue from 112 genome-wide genotyped fetuses from the second trimester of gestation, identifying high-confidence (false discovery rate < 0.05) expression quantitative trait loci for 30 mature miRNA. Integrating our findings with genome-wide association study data for brain-related disorders, we implicate increased prenatal expression of miR-1908-5p as a risk mechanism for bipolar disorder and find that predicted mRNA targets of miR-1908-5p that are expressed in the fetal brain are enriched for common variant genetic association with the condition. Extending these analyses to other brain-related traits, we find that common genetic variation associated with increased miR-1908-5p expression in fetal brain is additionally associated with depressive symptoms, irritability, increased right cerebellum exterior volume and increased sleep duration in the general population. Our findings provide support to the view that altered miRNA expression can influence susceptibility to neuropsychiatric illness and suggest an early neurodevelopmental risk mechanism for bipolar disorder.

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.

Cis-effects on gene expression in the human prenatal brain associated with genetic risk for neuropsychiatric disorders.

The majority of common risk alleles identified for neuropsychiatric disorders reside in noncoding regions of the genome and are therefore likely to impact gene regulation. However, the genes that are primarily affected and the nature and developmental timing of these effects remain unclear. Given the hypothesized role for early neurodevelopmental processes in these conditions, we here define genetic predictors of gene expression in the human fetal brain with which we perform transcriptome-wide association studies (TWASs) of attention deficit hyperactivity disorder (ADHD), autism spectrum disorder, bipolar disorder, major depressive disorder, and schizophrenia. We identify prenatal cis-regulatory effects on 63 genes and 166 individual transcripts associated with genetic risk for these conditions. We observe pleiotropic effects of expression predictors for a number of genes and transcripts, including those of decreased DDHD2 expression in association with risk for schizophrenia and bipolar disorder, increased expression of a ST3GAL3 transcript with risk for schizophrenia and ADHD, and increased expression of an XPNPEP3 transcript with risk for schizophrenia, bipolar disorder, and major depression. For the protocadherin alpha cluster genes PCDHA7 and PCDHA8, we find that predictors of low expression are associated with risk for major depressive disorder while those of higher expression are associated with risk for schizophrenia. Our findings support a role for altered gene regulation in the prenatal brain in susceptibility to various neuropsychiatric disorders and prioritize potential risk genes for further neurobiological investigation.

Sites of active gene regulation in the prenatal frontal cortex and their role in neuropsychiatric disorders.

Common genetic variation appears to largely influence risk for neuropsychiatric disorders through effects on gene regulation. It is therefore possible to shed light on the biology of these conditions by testing for enrichment of associated genetic variation within regulatory genomic regions operating in specific tissues or cell types. Here, we have used the assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-Seq) to map open chromatin (an index of active regulatory genomic regions) in bulk tissue, NeuN+ and NeuN- nuclei from the prenatal human frontal cortex, and tested enrichment of single-nucleotide polymorphism (SNP) heritability for five neuropsychiatric disorders (autism spectrum disorder, attention deficit hyperactivity disorder [ADHD], bipolar disorder, major depressive disorder, and schizophrenia) within these regions. We observed significant enrichment of SNP heritability for ADHD, major depressive disorder, and schizophrenia within open chromatin regions (OCRs) mapped in bulk fetal frontal cortex, and for all five tested neuropsychiatric conditions when we restricted these sites to those overlapping histone modifications indicative of enhancers (H3K4me1) or promoters (H3K4me3) in fetal brain. SNP heritability for neuropsychiatric disorders was significantly enriched in OCRs identified in fetal frontal cortex NeuN- as well as NeuN+ nuclei overlapping fetal brain H3K4me1 or H3K4me3 sites. We additionally demonstrate the utility of our mapped OCRs for prioritizing potentially functional SNPs at genome-wide significant risk loci for neuropsychiatric disorders. Our data provide evidence for an early neurodevelopmental component to a range of neuropsychiatric conditions and highlight an important role for regulatory genomic regions active within both NeuN+ and NeuN- cells of the prenatal brain.

Methylomic trajectories across human fetal brain development.

Epigenetic processes play a key role in orchestrating transcriptional regulation during development. The importance of DNA methylation in fetal brain development is highlighted by the dynamic expression of de novo DNA methyltransferases during the perinatal period and neurodevelopmental deficits associated with mutations in the methyl-CpG binding protein 2 (MECP2) gene. However, our knowledge about the temporal changes to the epigenome during fetal brain development has, to date, been limited. We quantified genome-wide patterns of DNA methylation at ∼ 400,000 sites in 179 human fetal brain samples (100 male, 79 female) spanning 23 to 184 d post-conception. We identified highly significant changes in DNA methylation across fetal brain development at >7% of sites, with an enrichment of loci becoming hypomethylated with fetal age. Sites associated with developmental changes in DNA methylation during fetal brain development were significantly underrepresented in promoter regulatory regions but significantly overrepresented in regions flanking CpG islands (shores and shelves) and gene bodies. Highly significant differences in DNA methylation were observed between males and females at a number of autosomal sites, with a small number of regions showing sex-specific DNA methylation trajectories across brain development. Weighted gene comethylation network analysis (WGCNA) revealed discrete modules of comethylated loci associated with fetal age that are significantly enriched for genes involved in neurodevelopmental processes. This is, to our knowledge, the most extensive study of DNA methylation across human fetal brain development to date, confirming the prenatal period as a time of considerable epigenomic plasticity.

5-hydroxymethylcytosine is highly dynamic across human fetal brain development.

BACKGROUND: Epigenetic processes play a key role in orchestrating transcriptional regulation during the development of the human central nervous system. We previously described dynamic changes in DNA methylation (5mC) occurring during human fetal brain development, but other epigenetic processes operating during this period have not been extensively explored. Of particular interest is DNA hydroxymethylation (5hmC), a modification that is enriched in the human brain and hypothesized to play an important role in neuronal function, learning and memory. In this study, we quantify 5hmC across the genome of 71 human fetal brain samples spanning 23 to 184 days post-conception. RESULTS: We identify widespread changes in 5hmC occurring during human brain development, notable sex-differences in 5hmC in the fetal brain, and interactions between 5mC and 5hmC at specific sites. Finally, we identify loci where 5hmC in the fetal brain is associated with genetic variation. CONCLUSIONS: This study represents the first systematic analysis of dynamic changes in 5hmC across human neurodevelopment and highlights the potential importance of this modification in the human brain. A searchable database of our fetal brain 5hmC data is available as a resource to the research community at http://www.epigenomicslab.com/online-data-resources .

Knockdown of the schizophrenia susceptibility gene TCF4 alters gene expression and proliferation of progenitor cells from the developing human neocortex.

BACKGROUND: Common variants in the TCF4 gene are among the most robustly supported genetic risk factors for schizophrenia. Rare TCF4 deletions and loss-of-function point mutations cause Pitt-Hopkins syndrome, a developmental disorder associated with severe intellectual disability. METHODS: To explore molecular and cellular mechanisms by which TCF4 perturbation could interfere with human cortical development, we experimentally reduced the endogenous expression of TCF4 in a neural progenitor cell line derived from the developing human cerebral cortex using RNA interference. Effects on genome-wide gene expression were assessed by microarray, followed by Gene Ontology and pathway analysis of differentially expressed genes. We tested for genetic association between the set of differentially expressed genes and schizophrenia using genome-wide association study data from the Psychiatric Genomics Consortium and competitive gene set analysis (MAGMA). Effects on cell proliferation were assessed using high content imaging. RESULTS: Genes that were differentially expressed following TCF4 knockdown were highly enriched for involvement in the cell cycle. There was a nonsignificant trend for genetic association between the differentially expressed gene set and schizophrenia. Consistent with the gene expression data, TCF4 knockdown was associated with reduced proliferation of cortical progenitor cells in vitro. LIMITATIONS: A detailed mechanistic explanation of how TCF4 knockdown alters human neural progenitor cell proliferation is not provided by this study. CONCLUSION: Our data indicate effects of TCF4 perturbation on human cortical progenitor cell proliferation, a process that could contribute to cognitive deficits in individuals with Pitt-Hopkins syndrome and risk for schizophrenia.

Genome-wide significant schizophrenia risk variation on chromosome 10q24 is associated with altered cis-regulation of BORCS7, AS3MT, and NT5C2 in the human brain.

Chromosome 10q24.32-q24.33 is one of the most robustly supported risk loci to emerge from genome-wide association studies (GWAS) of schizophrenia. However, extensive linkage disequilibrium makes it difficult to distinguish the actual susceptibility gene(s) at the locus, limiting its value for improving biological understanding of the condition. In the absence of coding changes that can account for the association, risk is likely conferred by altered regulation of one or more genes in the region. We, therefore, used highly sensitive measures of allele-specific expression to assess cis-regulatory effects associated with the two best-supported schizophrenia risk variants (SNP rs11191419 and indel ch10_104957618_I/rs202213518) on the primary positional candidates BORCS7, AS3MT, CNNM2, and NT5C2 in the human brain. Heterozygosity at rs11191419 was associated with increased allelic expression of BORCS7 and AS3MT in the fetal and adult brain, and with reduced allelic expression of NT5C2 in the adult brain. Heterozygosity at ch10_104957618_I was associated with reduced allelic expression of NT5C2 in both the fetal and adult brain. Comparisons between cDNA ratios in heterozygotes and homozygotes for the risk alleles indicated that cis-effects on NT5C2 expression in the adult dorsolateral prefrontal cortex could be largely accounted for by genotype at these two risk variants. While not excluding effects on other genes in the region, this study implicates altered neural expression of BORCS7, AS3MT, and NT5C2 in susceptibility to schizophrenia arising from genetic variation at the chromosome 10q24 locus. © 2016 The Authors. American Journal of Medical Genetics Part B: Neuropsychiatric Genetics Published by Wiley Periodicals, Inc.

Knockdown of the psychosis susceptibility gene ZNF804A alters expression of genes involved in cell adhesion.

Genome-wide association studies have convincingly implicated several novel genes in susceptibility to schizophrenia and bipolar disorder. The first genome-wide significant association with the broad phenotype of psychosis was with a polymorphism in the ZNF804A gene. However, the biological function(s) of ZNF804A have, to date, been entirely unknown. In this study, we manipulated the expression of ZNF804A in neural progenitor cells derived from human cortical neuroepithelium and assessed its effects on the cellular transcriptome. Gene ontology analysis of differentially expressed genes indicated a significant effect of ZNF804A knockdown on the expression of genes involved in cell adhesion, suggesting a role for ZNF804A in processes such as neural migration, neurite outgrowth and synapse formation. Several highly significant gene expression changes were confirmed in repeat cell culture experiments. Most consistent gene expression changes were seen for C2ORF80, a gene of as-yet-unknown function, and STMN3, a gene involved in neurite outgrowth and axonal and dendritic branching. These data, generated in a hypothesis-free manner, provide a basis for more targeted investigations of ZNF804A function.

Genetic Implication of Prenatal GABAergic and Cholinergic Neuron Development in Susceptibility to Schizophrenia.

BACKGROUND: The ganglionic eminences (GE) are fetal-specific structures that give rise to gamma-aminobutyric acid (GABA)- and acetylcholine-releasing neurons of the forebrain. Given the evidence for GABAergic, cholinergic, and neurodevelopmental disturbances in schizophrenia, we tested the potential involvement of GE neuron development in mediating genetic risk for the condition. STUDY DESIGN: We combined data from a recent large-scale genome-wide association study of schizophrenia with single-cell RNA sequencing data from the human GE to test the enrichment of schizophrenia risk variation in genes with high expression specificity for developing GE cell populations. We additionally performed the single nuclei Assay for Transposase-Accessible Chromatin with Sequencing (snATAC-Seq) to map potential regulatory genomic regions operating in individual cell populations of the human GE, using these to test for enrichment of schizophrenia common genetic variant liability and to functionally annotate non-coding variants-associated with the disorder. STUDY RESULTS: Schizophrenia common variant liability was enriched in genes with high expression specificity for developing neuron populations that are predicted to form dopamine D1 and D2 receptor-expressing GABAergic medium spiny neurons of the striatum, cortical somatostatin-positive GABAergic interneurons, calretinin-positive GABAergic neurons, and cholinergic neurons. Consistent with these findings, schizophrenia genetic risk was concentrated in predicted regulatory genomic sequence mapped in developing neuronal populations of the GE. CONCLUSIONS: Our study implicates prenatal development of specific populations of GABAergic and cholinergic neurons in later susceptibility to schizophrenia, and provides a map of predicted regulatory genomic elements operating in cells of the GE.

Cross-ancestry atlas of gene, isoform, and splicing regulation in the developing human brain.

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.

Stochastic choice of allelic expression in human neural stem cells.

Monoallelic gene expression, such as genomic imprinting, is well described. Less well-characterized are genes undergoing stochastic monoallelic expression (MA), where specific clones of cells express just one allele at a given locus. We performed genome-wide allelic expression assessment of human clonal neural stem cells derived from cerebral cortex, striatum, and spinal cord, each with differing genotypes. We assayed three separate clonal lines from each donor, distinguishing stochastic MA from genotypic effects. Roughly 2% of genes showed evidence for autosomal MA, and in about half of these, allelic expression was stochastic between different clones. Many of these loci were known neurodevelopmental genes, such as OTX2 and OLIG2. Monoallelic genes also showed increased levels of DNA methylation compared to hypomethylated biallelic loci. Identified monoallelic gene loci showed altered chromatin signatures in fetal brain, suggesting an in vivo correlate of this phenomenon. We conclude that stochastic allelic expression is prevalent in neural stem cells, providing clonal diversity to developing tissues such as the human brain.

Single-Nuclei RNA Sequencing of 5 Regions of the Human Prenatal Brain Implicates Developing Neuron Populations in Genetic Risk for Schizophrenia.

BACKGROUND: While a variety of evidence supports a prenatal component in schizophrenia, there are few data regarding the cell populations involved. We sought to identify cells of the human prenatal brain mediating genetic risk for schizophrenia by integrating cell-specific gene expression measures generated through single-nuclei RNA sequencing with recent large-scale genome-wide association study (GWAS) and exome sequencing data for the condition. METHODS: Single-nuclei RNA sequencing was performed on 5 brain regions (frontal cortex, ganglionic eminence, hippocampus, thalamus, and cerebellum) from 3 fetuses from the second trimester of gestation. Enrichment of schizophrenia common variant genetic liability and rare damaging coding variation was assessed in relation to gene expression specificity within each identified cell population. RESULTS: Common risk variants were prominently enriched within genes with high expression specificity for developing neuron populations within the frontal cortex, ganglionic eminence, and hippocampus. Enrichments were largely independent of genes expressed in neuronal populations of the adult brain that have been implicated in schizophrenia through the same methods. Genes containing an excess of rare damaging variants in schizophrenia had higher expression specificity for developing glutamatergic neurons of the frontal cortex and hippocampus that were also enriched for common variant liability. CONCLUSIONS: We found evidence for a distinct contribution of prenatal neuronal development to genetic risk for schizophrenia, involving specific populations of developing neurons within the second-trimester fetal brain. Our study significantly advances the understanding of the neurodevelopmental origins of schizophrenia and provides a resource with which to investigate the prenatal antecedents of other psychiatric and neurologic disorders.

Genetic implication of prenatal GABAergic and cholinergic neuron development in susceptibility to schizophrenia.

Background: The ganglionic eminences are fetal-specific structures that give rise to gamma-aminobutyric acid (GABA)- and acetylcholine- releasing neurons of the forebrain. Given evidence for GABAergic and cholinergic disturbances in schizophrenia, as well as an early neurodevelopmental component to the disorder, we tested the potential involvement of developing cells of the ganglionic eminences in mediating genetic risk for the condition. Study Design: We combined data from a recent large-scale genome-wide association study of schizophrenia with single cell RNA sequencing data from the human ganglionic eminences to test enrichment of schizophrenia risk variation in genes with high expression specificity for particular developing cell populations within these structures. We additionally performed the single nuclei Assay for Transposase-Accessible Chromatin with Sequencing (snATAC-Seq) to map potential regulatory genomic regions operating in individual cell populations of the human ganglionic eminences, using these to additionally test for enrichment of schizophrenia common genetic variant liability and to functionally annotate non-coding variants associated with the disorder. Study Results: Schizophrenia common variant liability was enriched in genes with high expression specificity for developing neuron populations that are predicted to form dopamine D1 and D2 receptor expressing GABAergic medium spiny neurons of the striatum, cortical somatostatin-positive GABAergic interneurons, calretinin-positive GABAergic neurons and cholinergic neurons. Consistent with these findings, schizophrenia genetic risk was also concentrated in predicted regulatory genomic sequence mapped in developing neuronal populations of the ganglionic eminences. Conclusions: Our study provides evidence for a role of prenatal GABAergic and cholinergic neuron development in later susceptibility to schizophrenia.

Cross-ancestry, cell-type-informed atlas of gene, isoform, and splicing regulation in the developing human brain.

Genomic regulatory elements active in the developing human brain are notably enriched in genetic risk for neuropsychiatric disorders, including autism spectrum disorder (ASD), schizophrenia, and bipolar disorder. However, prioritizing the specific risk genes and candidate molecular mechanisms underlying these genetic enrichments has been hindered by the lack of a single unified large-scale gene regulatory atlas of human brain development. Here, we uniformly process and systematically characterize gene, isoform, and splicing quantitative trait loci (xQTLs) in 672 fetal brain samples from unique subjects across multiple ancestral populations. We identify 15,752 genes harboring a significant xQTL and map 3,739 eQTLs to a specific cellular context. We observe a striking drop in gene expression and splicing heritability as the human brain develops. Isoform-level regulation, particularly in the second trimester, mediates the greatest proportion of heritability across multiple psychiatric GWAS, compared with eQTLs. Via colocalization and TWAS, we prioritize biological mechanisms for ~60% of GWAS loci across five neuropsychiatric disorders, nearly two-fold that observed in the adult brain. Finally, we build a comprehensive set of developmentally regulated gene and isoform co-expression networks capturing unique genetic enrichments across disorders. Together, this work provides a comprehensive view of genetic regulation across human brain development as well as the stage-and cell type-informed mechanistic underpinnings of neuropsychiatric disorders.

Effects of cis-regulatory variation differ across regions of the adult human brain.

Cis-regulatory variation is considered to be an important determinant of human phenotypic variability, including susceptibility to complex disease. Recent studies have shown that the effects of cis-regulatory polymorphism on gene expression can differ widely between tissues. In the present study, we tested whether the effects of cis-regulatory variation can also differ between regions of the adult human brain. We used relative allelic expression to measure cis-effects on the RNA expression of five candidate genes for neuropsychiatric illness (ZNF804A, NOS1, RGS4, AKT1 and TCF4) across multiple discrete brain regions within individual subjects. For all five genes, we observed significant differences in allelic expression between brain regions in several individual subjects, suggesting regional differences in the effects of cis-regulatory polymorphism to be a common phenomenon. As well as highlighting an important caveat for studies of regulatory polymorphism in the brain, our findings indicate that it is possible to delineate brain areas in which cis-regulatory variants are active. This may provide important insights into the fundamental biology of neuropsychiatric phenotypes with which such variants are associated.

Schizophrenia Genomics: Convergence on Synaptic Development, Adult Synaptic Plasticity, or Both?

Large-scale genomic studies of schizophrenia have identified hundreds of genetic loci conferring risk to the disorder. This progress offers an important route toward defining the biological basis of the condition and potentially developing new treatments. In this review, we discuss insights from recent genome-wide association study, copy number variant, and exome sequencing analyses of schizophrenia, together with functional genomics data from the pre- and postnatal brain, in relation to synaptic development and function. These data provide strong support for the view that synaptic dysfunction within glutamatergic and GABAergic (gamma-aminobutyric acidergic) neurons of the cerebral cortex, hippocampus, and other limbic structures is a central component of schizophrenia pathophysiology. Implicated genes and functional genomic data suggest that disturbances in synaptic connectivity associated with susceptibility to schizophrenia begin in utero but continue throughout development, with some alleles conferring risk to the disorder through direct effects on synaptic function in adulthood. This model implies that novel interventions for schizophrenia could include broad preventive approaches aimed at enhancing synaptic health during development as well as more targeted treatments aimed at correcting synaptic function in affected adults.

CSF-resident CD4+ T-cells display a distinct gene expression profile with relevance to immune surveillance and multiple sclerosis.

The CNS has traditionally been considered an immune privileged site, but is now understood to have a system of immune surveillance, predominantly involving CD4+ T-cells. Identifying functional differences between CNS and blood CD4+ T-cells, therefore, have relevance to CNS immune surveillance as well as to neurological conditions, such as multiple sclerosis, in which CD4+ T-cells play a central role. Here, CD4+ T-cells were purified from CSF and blood from 21 patients with newly diagnosed treatment-naïve multiple sclerosis and 20 individuals with non-inflammatory disorders using fluorescence-activated cell sorting, and their transcriptomes were profiled by RNA sequencing. Paired comparisons between CD4+ T-cells from CSF and blood identified 5156 differentially expressed genes in controls and 4263 differentially expressed in multiple sclerosis patients at false discovery rate <5%. Differential expression analysis of CD4+ T-cells collected from the CSF highlighted genes involved in migration, activation, cholesterol biosynthesis and signalling, including those with known relevance to multiple sclerosis pathogenesis and treatment. Expression of markers of CD4+ T-cell subtypes suggested an increased proportion of Th1 and Th17 cells in CSF. Gene ontology terms significant only in multiple sclerosis were predominantly those involved in cellular proliferation. A two-way comparison of CSF versus blood CD4+ T-cells in multiple sclerosis compared with non-inflammatory disorder controls identified four significant genes at false discovery rate <5% (CYP51A1, LRRD1, YES1 and PASK), further implicating cholesterol biosynthesis and migration mechanisms. Analysis of CSF CD4+ T-cells in an extended cohort of multiple sclerosis cases (total N = 41) compared with non-inflammatory disorder controls (total N = 38) identified 140 differentially expressed genes at false discovery rate < 5%, many of which have known relevance to multiple sclerosis, including XBP1, BHLHE40, CD40LG, DPP4 and ITGB1. This study provides the largest transcriptomic analysis of purified cell subpopulations in CSF to date and has relevance for the understanding of CNS immune surveillance, as well as multiple sclerosis pathogenesis and treatment discovery.

Full-length transcript sequencing of human and mouse cerebral cortex identifies widespread isoform diversity and alternative splicing.

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.

Novel epigenetic clock for fetal brain development predicts prenatal age for cellular stem cell models and derived neurons.

Induced pluripotent stem cells (iPSCs) and their differentiated neurons (iPSC-neurons) are a widely used cellular model in the research of the central nervous system. However, it is unknown how well they capture age-associated processes, particularly given that pluripotent cells are only present during the earliest stages of mammalian development. Epigenetic clocks utilize coordinated age-associated changes in DNA methylation to make predictions that correlate strongly with chronological age. It has been shown that the induction of pluripotency rejuvenates predicted epigenetic age. As existing clocks are not optimized for the study of brain development, we developed the fetal brain clock (FBC), a bespoke epigenetic clock trained in human prenatal brain samples in order to investigate more precisely the epigenetic age of iPSCs and iPSC-neurons. The FBC was tested in two independent validation cohorts across a total of 194 samples, confirming that the FBC outperforms other established epigenetic clocks in fetal brain cohorts. We applied the FBC to DNA methylation data from iPSCs and embryonic stem cells and their derived neuronal precursor cells and neurons, finding that these cell types are epigenetically characterized as having an early fetal age. Furthermore, while differentiation from iPSCs to neurons significantly increases epigenetic age, iPSC-neurons are still predicted as being fetal. Together our findings reiterate the need to better understand the limitations of existing epigenetic clocks for answering biological research questions and highlight a limitation of iPSC-neurons as a cellular model of age-related diseases.