Human prefrontal cortex gene regulatory dynamics from gestation to adulthood at single-cell resolution.

Human brain development is underpinned by cellular and molecular reconfigurations continuing into the third decade of life. To reveal cell dynamics orchestrating neural maturation, we profiled human prefrontal cortex gene expression and chromatin accessibility at single-cell resolution from gestation to adulthood. Integrative analyses define the dynamic trajectories of each cell type, revealing major gene expression reconfiguration at the prenatal-to-postnatal transition in all cell types followed by continuous reconfiguration into adulthood and identifying regulatory networks guiding cellular developmental programs, states, and functions. We uncover links between expression dynamics and developmental milestones, characterize the diverse timing of when cells acquire adult-like states, and identify molecular convergence from distinct developmental origins. We further reveal cellular dynamics and their regulators implicated in neurological disorders. Finally, using this reference, we benchmark cell identities and maturation states in organoid models. Together, this captures the dynamic regulatory landscape of human cortical development.

De novo variants predicting haploinsufficiency for DIP2C are associated with expressive speech delay.

The disconnected (disco)-interacting protein 2 (DIP2) gene was first identified in D. melanogaster and contains a DNA methyltransferase-associated protein 1 (DMAP1) binding domain, Acyl-CoA synthetase domain and AMP-binding sites. DIP2 regulates axonal bifurcation of the mushroom body neurons in D. melanogaster and is required for axonal regeneration in the neurons of C. elegans. The DIP2 homologues in vertebrates, Disco-interacting protein 2 homolog A (DIP2A), Disco-interacting protein 2 homolog B (DIP2B), and Disco-interacting protein 2 homolog C (DIP2C), are highly conserved and expressed widely in the central nervous system. Although there is evidence that DIP2C plays a role in cognition, reports of pathogenic variants in these genes are rare and their significance is uncertain. We present 23 individuals with heterozygous DIP2C variants, all manifesting developmental delays that primarily affect expressive language and speech articulation. Eight patients had de novo variants predicting loss-of-function in the DIP2C gene, two patients had de novo missense variants, three had paternally inherited loss of function variants and six had maternally inherited loss-of-function variants, while inheritance was unknown for four variants. Four patients had cardiac defects (hypertrophic cardiomyopathy, atrial septal defects, and bicuspid aortic valve). Minor facial anomalies were inconsistent but included a high anterior hairline with a long forehead, broad nasal tip, and ear anomalies. Brainspan analysis showed elevated DIP2C expression in the human neocortex at 10-24 weeks after conception. With the cases presented herein, we provide phenotypic and genotypic data supporting the association between loss-of-function variants in DIP2C with a neurocognitive phenotype.

Supervised application of internal validation measures to benchmark dimensionality reduction methods in scRNA-seq data.

A typical single-cell RNA sequencing (scRNA-seq) experiment will measure on the order of 20 000 transcripts and thousands, if not millions, of cells. The high dimensionality of such data presents serious complications for traditional data analysis methods and, as such, methods to reduce dimensionality play an integral role in many analysis pipelines. However, few studies have benchmarked the performance of these methods on scRNA-seq data, with existing comparisons assessing performance via downstream analysis accuracy measures, which may confound the interpretation of their results. Here, we present the most comprehensive benchmark of dimensionality reduction methods in scRNA-seq data to date, utilizing over 300 000 compute hours to assess the performance of over 25 000 low-dimension embeddings across 33 dimensionality reduction methods and 55 scRNA-seq datasets. We employ a simple, yet novel, approach, which does not rely on the results of downstream analyses. Internal validation measures (IVMs), traditionally used as an unsupervised method to assess clustering performance, are repurposed to measure how well-formed biological clusters are after dimensionality reduction. Performance was further evaluated over nearly 200 000 000 iterations of DBSCAN, a density-based clustering algorithm, showing that hyperparameter optimization using IVMs as the objective function leads to near-optimal clustering. Methods were also assessed on the extent to which they preserve the global structure of the data, and on their computational memory and time requirements across a large range of sample sizes. Our comprehensive benchmarking analysis provides a valuable resource for researchers and aims to guide best practice for dimensionality reduction in scRNA-seq analyses, and we highlight Latent Dirichlet Allocation and Potential of Heat-diffusion for Affinity-based Transition Embedding as high-performing algorithms.

NeuroCirc: an integrative resource of circular RNA expression in the human brain.

MOTIVATION: CircRNAs are covalently closed RNA molecules that are particularly abundant in the brain. While circRNA expression data from the human brain is rapidly accumulating, integration of large-scale datasets remains challenging and time-consuming, and consequently an integrative view of circRNA expression in the human brain is currently lacking. RESULTS: NeuroCirc is a web-based resource that allows interactive exploration of multiple types of circRNA data from the human brain, including large-scale expression datasets, circQTL data and circRNA expression across neuronal differentiation and cellular maturation time-courses. NeuroCirc also allows users to upload their own circRNA expression data and explore it in the integrative platform, thereby supporting circRNA prioritization for experimental validation and functional studies. AVAILABILITY AND IMPLEMENTATION: NeuroCirc is freely available at: https://voineagulab.github.io/NeuroCirc/. The source code and user documentation are available at: https://github.com/Voineagulab/NeuroCirc. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

De Novo ZMYND8 variants result in an autosomal dominant neurodevelopmental disorder with cardiac malformations.

PURPOSE: ZMYND8 encodes a multidomain protein that serves as a central interactive hub for coordinating critical roles in transcription regulation, chromatin remodeling, regulation of super-enhancers, DNA damage response and tumor suppression. We delineate a novel neurocognitive disorder caused by variants in the ZMYND8 gene. METHODS: An international collaboration, exome sequencing, molecular modeling, yeast two-hybrid assays, analysis of available transcriptomic data and a knockdown Drosophila model were used to characterize the ZMYND8 variants. RESULTS: ZMYND8 variants were identified in 11 unrelated individuals; 10 occurred de novo and one suspected de novo; 2 were truncating, 9 were missense, of which one was recurrent. The disorder is characterized by intellectual disability with variable cardiovascular, ophthalmologic and minor skeletal anomalies. Missense variants in the PWWP domain of ZMYND8 abolish the interaction with Drebrin and missense variants in the MYND domain disrupt the interaction with GATAD2A. ZMYND8 is broadly expressed across cell types in all brain regions and shows highest expression in the early stages of brain development. Neuronal knockdown of the DrosophilaZMYND8 ortholog results in decreased habituation learning, consistent with a role in cognitive function. CONCLUSION: We present genomic and functional evidence for disruption of ZMYND8 as a novel etiology of syndromic intellectual disability.

Chromosome conformation elucidates regulatory relationships in developing human brain.

Three-dimensional physical interactions within chromosomes dynamically regulate gene expression in a tissue-specific manner. However, the 3D organization of chromosomes during human brain development and its role in regulating gene networks dysregulated in neurodevelopmental disorders, such as autism or schizophrenia, are unknown. Here we generate high-resolution 3D maps of chromatin contacts during human corticogenesis, permitting large-scale annotation of previously uncharacterized regulatory relationships relevant to the evolution of human cognition and disease. Our analyses identify hundreds of genes that physically interact with enhancers gained on the human lineage, many of which are under purifying selection and associated with human cognitive function. We integrate chromatin contacts with non-coding variants identified in schizophrenia genome-wide association studies (GWAS), highlighting multiple candidate schizophrenia risk genes and pathways, including transcription factors involved in neurogenesis, and cholinergic signalling molecules, several of which are supported by independent expression quantitative trait loci and gene expression analyses. Genome editing in human neural progenitors suggests that one of these distal schizophrenia GWAS loci regulates FOXG1 expression, supporting its potential role as a schizophrenia risk gene. This work provides a framework for understanding the effect of non-coding regulatory elements on human brain development and the evolution of cognition, and highlights novel mechanisms underlying neuropsychiatric disorders.

TDAview: an online visualization tool for topological data analysis.

SUMMARY: TDAview is an online tool for topological data analysis (TDA) and visualization. It implements the Mapper algorithm for TDA and provides extensive graph visualization options. TDAview is a user-friendly tool that allows biologists and clinicians without programming knowledge to harness the power of TDA. TDAview supports an analysis and visualization mode in which a Mapper graph is constructed based on user-specified parameters, followed by graph visualization. It can also be used in a visualization only mode in which TDAview is used for visualizing the data properties of a Mapper graph generated using other open-source software. The graph visualization options allow data exploration by graphical display of metadata variable values for nodes and edges, as well as the generation of publishable figures. TDAview can handle large datasets, with tens of thousands of data points, and thus has a wide range of applications for high-dimensional data, including the construction of topology-based gene co-expression networks. AVAILABILITY AND IMPLEMENTATION: TDAview is a free online tool available at https://voineagulab.github.io/TDAview/. The source code, usage documentation and example data are available at TDAview GitHub repository: https://github.com/Voineagulab/TDAview.

Transcriptional response to mitochondrial protease IMMP2L knockdown in human primary astrocytes.

IMMP2L encodes the inner membrane peptidase subunit 2, a mitochondrial protease involved in cleaving the space-sorting signals of mitochondrial membrane proteins. IMMP2L has been implicated in Tourette syndrome, but how its dysfunction contributes to the neurodevelopmental phenotype remains unclear. Here we show that IMMP2L transcription requires Topoisomerase I in human primary astrocytes, and characterize the downstream effects of IMMP2L knockdown on gene expression. We demonstrate that IMMP2L knockdown leads to dysregulation of genes involved in central nervous system development. We also find that the transcriptional response to IMMP2L knockdown partially overlaps the one induced by mitochondrial complex III inhibition. Overall, these data bring further insight into the molecular consequences of IMMP2L dysfunction in the brain.

Transcriptomic analysis of autistic brain reveals convergent molecular pathology.

Autism spectrum disorder (ASD) is a common, highly heritable neurodevelopmental condition characterized by marked genetic heterogeneity. Thus, a fundamental question is whether autism represents an aetiologically heterogeneous disorder in which the myriad genetic or environmental risk factors perturb common underlying molecular pathways in the brain. Here, we demonstrate consistent differences in transcriptome organization between autistic and normal brain by gene co-expression network analysis. Remarkably, regional patterns of gene expression that typically distinguish frontal and temporal cortex are significantly attenuated in the ASD brain, suggesting abnormalities in cortical patterning. We further identify discrete modules of co-expressed genes associated with autism: a neuronal module enriched for known autism susceptibility genes, including the neuronal specific splicing factor A2BP1 (also known as FOX1), and a module enriched for immune genes and glial markers. Using high-throughput RNA sequencing we demonstrate dysregulated splicing of A2BP1-dependent alternative exons in the ASD brain. Moreover, using a published autism genome-wide association study (GWAS) data set, we show that the neuronal module is enriched for genetically associated variants, providing independent support for the causal involvement of these genes in autism. In contrast, the immune-glial module showed no enrichment for autism GWAS signals, indicating a non-genetic aetiology for this process. Collectively, our results provide strong evidence for convergent molecular abnormalities in ASD, and implicate transcriptional and splicing dysregulation as underlying mechanisms of neuronal dysfunction in this disorder.

Large-scale expansions of Friedreich's ataxia GAA repeats in yeast.

Large-scale expansions of DNA repeats are implicated in numerous hereditary disorders in humans. We describe a yeast experimental system to analyze large-scale expansions of triplet GAA repeats responsible for the human disease Friedreich's ataxia. When GAA repeats were placed into an intron of the chimeric URA3 gene, their expansions caused gene inactivation, which was detected on the selective media. We found that the rates of expansions of GAA repeats increased exponentially with their lengths. These rates were only mildly dependent on the repeat's orientation within the replicon, whereas the repeat-mediated replication fork stalling was exquisitely orientation dependent. Expansion rates were significantly elevated upon inactivation of the replication fork stabilizers, Tof1 and Csm3, but decreased in the knockouts of postreplication DNA repair proteins, Rad6 and Rad5, and the DNA helicase Sgs1. We propose a model for large-scale repeat expansions based on template switching during replication fork progression through repetitive DNA.

Searching for convergent pathways in autism spectrum disorders: insights from human brain transcriptome studies.

Autism spectrum disorder (ASD) is one of the most heritable neuropsychiatric conditions. The complex genetic landscape of the disorder includes both common and rare variants at hundreds of genetic loci. This marked heterogeneity has thus far hampered efforts to develop genetic diagnostic panels and targeted pharmacological therapies. Here, we give an overview of the current literature on the genetic basis of ASD, and review recent human brain transcriptome studies and their role in identifying convergent pathways downstream of the heterogeneous genetic variants. We also discuss emerging evidence on the involvement of non-coding genomic regions and non-coding RNAs in ASD.

Transcriptome analysis of human brain tissue identifies reduced expression of complement complex C1Q Genes in Rett syndrome.

BACKGROUND: MECP2, the gene mutated in the majority of Rett syndrome cases, is a transcriptional regulator that can activate or repress transcription. Although the transcription regulatory function of MECP2 has been known for over a decade, it remains unclear how transcriptional dysregulation leads to the neurodevelopmental disorder. Notably, little convergence was previously observed between the genes abnormally expressed in the brain of Rett syndrome mouse models and those identified in human studies. METHODS: Here we carried out a comprehensive transcriptome analysis of human brain tissue from Rett syndrome brain using both RNA-seq and microarrays. RESULTS: We identified over two hundred differentially expressed genes, and identified the complement C1Q complex genes (C1QA, C1QB and C1QC) as a point of convergence between gene expression changes in human and mouse Rett syndrome brain. CONCLUSIONS: The results of our study support a role for alterations in the expression level of C1Q complex genes in RTT pathogenesis.

Genome-wide transcriptome profiling reveals the functional impact of rare de novo and recurrent CNVs in autism spectrum disorders.

Copy-number variants (CNVs) are a major contributor to the pathophysiology of autism spectrum disorders (ASDs), but the functional impact of CNVs remains largely unexplored. Because brain tissue is not available from most samples, we interrogated gene expression in lymphoblasts from 244 families with discordant siblings in the Simons Simplex Collection in order to identify potentially pathogenic variation. Our results reveal that the overall frequency of significantly misexpressed genes (which we refer to here as outliers) identified in probands and unaffected siblings does not differ. However, in probands, but not their unaffected siblings, the group of outlier genes is significantly enriched in neural-related pathways, including neuropeptide signaling, synaptogenesis, and cell adhesion. We demonstrate that outlier genes cluster within the most pathogenic CNVs (rare de novo CNVs) and can be used for the prioritization of rare CNVs of potentially unknown significance. Several nonrecurrent CNVs with significant gene-expression alterations are identified (these include deletions in chromosomal regions 3q27, 3p13, and 3p26 and duplications at 2p15), suggesting that these are potential candidate ASD loci. In addition, we identify distinct expression changes in 16p11.2 microdeletions, 16p11.2 microduplications, and 7q11.23 duplications, and we show that specific genes within the 16p CNV interval correlate with differences in head circumference, an ASD-relevant phenotype. This study provides evidence that pathogenic structural variants have a functional impact via transcriptome alterations in ASDs at a genome-wide level and demonstrates the utility of integrating gene expression with mutation data for the prioritization of genes disrupted by potentially pathogenic mutations.

Gene expression studies in autism: moving from the genome to the transcriptome and beyond.

Autism is a clinically and genetically heterogeneous neurodevelopmental disorder. Although multiple genes, risk alleles and copy number variants (CNVs) have been implicated in ASD, none of the currently established genetic causes of ASD accounts for more than 2% of the cases, and a genetic diagnosis is not yet possible for most autism patients. Thus, advancing our understanding of autism genetics requires the integration of genetic information with information on genome function, as provided by transcriptomic data. We review recent autism transcriptome studies, in the context of current knowledge of autism genetics, and discuss the utility of gene expression data in evaluating the functional relevance of genetic variants and identifying common molecular pathways dysregulated in autism.

Alu-minating the Mechanisms Underlying Primate Cortex Evolution.

The higher-order cognitive functions observed in primates correlate with the evolutionary enhancement of cortical volume and folding, which in turn are driven by the primate-specific expansion of cellular diversity in the developing cortex. Underlying these changes is the diversification of molecular features including the creation of human and/or primate-specific genes, the activation of specific molecular pathways, and the interplay of diverse layers of gene regulation. We review and discuss evidence for connections between Alu elements and primate brain evolution, the evolutionary milestones of which are known to coincide along primate lineages. Alus are repetitive elements that contribute extensively to the acquisition of novel genes and the expansion of diverse gene regulatory layers, including enhancers, alternative splicing, RNA editing, and microRNA pathways. By reviewing the impact of Alus on molecular features linked to cortical expansions or gyrification or implications in cognitive deficits, we suggest that future research focusing on the role of Alu-derived molecular events in the context of brain development may greatly advance our understanding of higher-order cognitive functions and neurologic disorders.

Comprehensive evaluation of deconvolution methods for human brain gene expression.

Transcriptome deconvolution aims to estimate the cellular composition of an RNA sample from its gene expression data, which in turn can be used to correct for composition differences across samples. The human brain is unique in its transcriptomic diversity, and comprises a complex mixture of cell-types, including transcriptionally similar subtypes of neurons. Here, we carry out a comprehensive evaluation of deconvolution methods for human brain transcriptome data, and assess the tissue-specificity of our key observations by comparison with human pancreas and heart. We evaluate eight transcriptome deconvolution approaches and nine cell-type signatures, testing the accuracy of deconvolution using in silico mixtures of single-cell RNA-seq data, RNA mixtures, as well as nearly 2000 human brain samples. Our results identify the main factors that drive deconvolution accuracy for brain data, and highlight the importance of biological factors influencing cell-type signatures, such as brain region and in vitro cell culturing.

Replication stalling at unstable inverted repeats: interplay between DNA hairpins and fork stabilizing proteins.

DNA inverted repeats (IRs) are hotspots of genomic instability in both prokaryotes and eukaryotes. This feature is commonly attributed to their ability to fold into hairpin- or cruciform-like DNA structures interfering with DNA replication and other genetic processes. However, direct evidence that IRs are replication stall sites in vivo is currently lacking. Here, we show by 2D electrophoretic analysis of replication intermediates that replication forks stall at IRs in bacteria, yeast, and mammalian cells. We found that DNA hairpins, rather than DNA cruciforms, are responsible for the replication stalling by comparing the effects of specifically designed imperfect IRs with varying lengths of their central spacer. Finally, we report that yeast fork-stabilizing proteins, Tof1 and Mrc1, are required to counteract repeat-mediated replication stalling. We show that the function of the Tof1 protein at DNA structure-mediated stall sites is different from its previously described effect on protein-mediated replication fork barriers.

Evolutionary dynamics of circular RNAs in primates.

Many primate genes produce circular RNAs (circRNAs). However, the extent of circRNA conservation between closely related species remains unclear. By comparing tissue-specific transcriptomes across over 70 million years of primate evolution, we identify that within 3 million years circRNA expression profiles diverged such that they are more related to species identity than organ type. However, our analysis also revealed a subset of circRNAs with conserved neural expression across tens of millions of years of evolution. By comparing to species-specific circRNAs, we identified that the downstream intron of the conserved circRNAs display a dramatic lengthening during evolution due to the insertion of novel retrotransposons. Our work provides comparative analyses of the mechanisms promoting circRNAs to generate increased transcriptomic complexity in primates.

Transcriptomic signature of early life stress in male rat prefrontal cortex.

Early life stress (ELS) is associated with adverse mental health outcomes including anxiety, depression and addiction-like behaviours. While ELS is known to affect the developing brain, leading to increased stress responsiveness and increased glucocorticoid levels, the molecular mechanisms underlying the detrimental effects of ELS remain incompletely characterised. Rodent models have been instrumental in beginning to uncover the molecular and cellular underpinnings of ELS. Limited nesting (LN), an ELS behavioural paradigm with significant improvements over maternal separation, mimics human maternal neglect. We have previously shown that LN leads to an increase in one of the behavioural measures of anxiety like-behaviours in rats (percent of entries in the EPM open arm). Here we assessed gene expression changes induced by ELS in rat prefrontal cortex by RNA-sequencing. We show that LN leads primarily to transcriptional repression and identify a molecular signature of LN in rat PFC that is observed across ELS protocols and replicable across rodent species (mouse and rat).

The Landscape of Circular RNA Expression in the Human Brain.

BACKGROUND: Circular RNAs (circRNAs) are enriched in the mammalian brain and upregulated in response to neuronal differentiation and depolarization. These RNA molecules, formed by noncanonical back-splicing, have both regulatory and translational potential. METHODS: Here, we carried out an extensive characterization of circRNA expression in the human brain, in nearly 200 human brain samples, from both healthy controls and autism cases. RESULTS: We identified hundreds of novel circRNAs and demonstrated that circRNAs are not expressed stochastically, but rather as major isoforms. We characterized interindividual variability of circRNA expression in the human brain and showed that interindividual variability is less pronounced than variability between the cerebral cortex and cerebellum. Finally, we identified a circRNA coexpression module upregulated in autism samples, thereby adding another layer of complexity to the transcriptome changes observed in the autism brain. CONCLUSIONS: These data provide a comprehensive catalog of circRNAs, as well as a deeper insight into their expression in the human brain, and are available as a free resource in browsable format (http://www.voineagulab.unsw.edu.au/circ_rna).