The recent Zika virus (ZIKV) outbreak in the Western hemisphere is associated with severe pathology in newborns, including microcephaly and brain damage. The mechanisms underlying these outcomes are under intense investigation. Here, we show that a 2015 ZIKV isolate replicates in multiple cell types, including primary human fetal neural progenitors (hNPs). In immortalized cells, ZIKV is cytopathic and grossly rearranges endoplasmic reticulum membranes similar to other flaviviruses. In hNPs, ZIKV infection has a partial cytopathic phase characterized by cell rounding, pyknosis, and activation of caspase 3. Despite notable cell death, ZIKV did not activate a cytokine response in hNPs. This lack of cell intrinsic immunity to ZIKV is consistent with our observation that virus replication persists in hNPs for at least 28 days. These findings, supported by published fetal neuropathology, establish a proof-of-concept that neural progenitors in the developing human fetus can be direct targets of detrimental ZIKV-induced pathology.

The postsynaptic density (PSD) is a protein dense complex on the postsynaptic membrane of excitatory synapses that is implicated in normal nervous system functions such as synaptic plasticity, and also contains an enrichment of proteins involved in neuropsychiatric disorders. It has recently been reported that the genes encoding PSD proteins evolved more slowly than other genes in the human brain, but the underlying evolutionary advantage for this is not clear. Here, we show that cortical gene expression levels could explain most of this effect, indicating that expression level is a primary contributor to the evolution of these genes in the brain. Furthermore, we identify a positive correlation between the expression of PSD genes and cortical layers, with PSD genes being more highly expressed in deep layers, likely as a result of layer-enriched transcription factors. As the cortical layers of the mammalian brain have distinct functions and anatomical projections, our results indicate that the emergence of the unique six-layered mammalian cortex may have provided differential functional constraints on the evolution of PSD genes. More superficial cortical layers contain PSD genes with less constraint and these layers are primarily involved in intracortical projections, connections that may be particularly important for evolved cognitive functions. Therefore, the differential expression and evolutionary constraint of PSD genes in neocortical laminae may be critical not only for neocortical architecture but the cognitive functions that are dependent on this structure.

The signalling pathways controlling both the evolution and development of language in the human brain remain unknown. So far, the transcription factor FOXP2 (forkhead box P2) is the only gene implicated in Mendelian forms of human speech and language dysfunction. It has been proposed that the amino acid composition in the human variant of FOXP2 has undergone accelerated evolution, and this two-amino-acid change occurred around the time of language emergence in humans. However, this remains controversial, and whether the acquisition of these amino acids in human FOXP2 has any functional consequence in human neurons remains untested. Here we demonstrate that these two human-specific amino acids alter FOXP2 function by conferring differential transcriptional regulation in vitro. We extend these observations in vivo to human and chimpanzee brain, and use network analysis to identify novel relationships among the differentially expressed genes. These data provide experimental support for the functional relevance of changes in FOXP2 that occur on the human lineage, highlighting specific pathways with direct consequences for human brain development and disease in the central nervous system (CNS). Because FOXP2 has an important role in speech and language in humans, the identified targets may have a critical function in the development and evolution of language circuitry in humans.

Methylglyoxal (MG) is a reactive and cytotoxic α-dicarbonyl byproduct of glycolysis. Our bodies have several bio-defense systems to detoxify MG, including an enzymatic system by glyoxalase (GLO) 1 and GLO2. We identified a subtype of schizophrenia patients with novel mutations in the GLO1 gene that results in reductions of enzymatic activity. Moreover, we found that vitamin B6 (VB6) levels in peripheral blood of the schizophrenia patients with GLO1 dysfunction are significantly lower than that of healthy controls. However, the effects of GLO1 dysfunction and VB6 deficiency on the pathophysiology of schizophrenia remains poorly understood. Here, we generated a novel mouse model for this subgroup of schizophrenia patients by feeding Glo1 knockout mice VB6-deficent diets (KO/VB6(-)) and evaluated the combined effects of GLO1 dysfunction and VB6 deficiency on brain function. KO/VB6(-) mice accumulated homocysteine in plasma and MG in the prefrontal cortex (PFC), hippocampus, and striatum, and displayed behavioral deficits, such as impairments of social interaction and cognitive memory and a sensorimotor deficit in the prepulse inhibition test. Furthermore, we found aberrant gene expression related to mitochondria function in the PFC of the KO/VB6(-) mice by RNA-sequencing and weighted gene co-expression network analysis (WGCNA). Finally, we demonstrated abnormal mitochondrial respiratory function and subsequently enhanced oxidative stress in the PFC of KO/VB6(-) mice in the PFC. These findings suggest that the combination of GLO1 dysfunction and VB6 deficiency may cause the observed behavioral deficits via mitochondrial dysfunction and oxidative stress in the PFC.

The hippocampus supports many facets of cognition, including learning, memory, and emotional processing. Anatomically, the hippocampus runs along a longitudinal axis, posterior to anterior in primates. The structure, function, and connectivity of the hippocampus vary along this axis. In human hippocampus, longitudinal functional heterogeneity remains an active area of investigation, and structural heterogeneity has not been described. To understand the cellular and molecular diversity along the hippocampal long axis in human brain and define molecular signatures corresponding to functional domains, we performed single-nuclei RNA sequencing on surgically resected human anterior and posterior hippocampus from epilepsy patients, identifying differentially expressed genes at cellular resolution. We further identify axis- and cell-type-specific gene expression signatures that differentially intersect with human genetic signals, identifying cell-type-specific genes in the posterior hippocampus for cognitive function and the anterior hippocampus for mood and affect. These data are accessible as a public resource through an interactive website.

Forkhead box P2 (FOXP2) is a transcription factor expressed in the human brain that peaks during fetal development, and disruption in its ability to regulate downstream target genes leads to vulnerability to neurodevelopmental disorders. However, the mechanisms by which FOXP2 exerts regulatory control over targets during neuronal maturation have not been fully elucidated. Here, we use genome-wide chromatin accessibility assays and transcriptome-wide expression analyses in differentiating human neurons to show that FOXP2 represses proliferation-promoting genes in a DNA-binding-dependent manner. In contrast, FOXP2 and its cofactors, NFIA and NFIB, activate neuronal maturation genes in a manner that does not require FOXP2 to interact with DNA directly. Moreover, comparisons with expression data from the developing human brain suggest that FOXP2 and NFIA- or NFIB-dependent chromatin alterations drive maturation of excitatory cortical neurons. Thus, FOXP2 and its NFI cofactors may be specifically important for the development of cortical circuits underlying neurodevelopmental disorders.

DNA methylation is a pervasive epigenetic DNA modification that strongly affects chromatin regulation and gene expression. To date, it remains largely unknown how patterns of DNA methylation differ between closely related species and whether such differences contribute to species-specific phenotypes. To investigate these questions, we generated nucleotide-resolution whole-genome methylation maps of the prefrontal cortex of multiple humans and chimpanzees. Levels and patterns of DNA methylation vary across individuals within species according to the age and the sex of the individuals. We also found extensive species-level divergence in patterns of DNA methylation and that hundreds of genes exhibit significantly lower levels of promoter methylation in the human brain than in the chimpanzee brain. Furthermore, we investigated the functional consequences of methylation differences in humans and chimpanzees by integrating data on gene expression generated with next-generation sequencing methods, and we found a strong relationship between differential methylation and gene expression. Finally, we found that differentially methylated genes are strikingly enriched with loci associated with neurological disorders, psychological disorders, and cancers. Our results demonstrate that differential DNA methylation might be an important molecular mechanism driving gene-expression divergence between human and chimpanzee brains and might potentially contribute to the evolution of disease vulnerabilities. Thus, comparative studies of humans and chimpanzees stand to identify key epigenomic modifications underlying the evolution of human-specific traits.

The molecular mechanisms driving brain development at risk in autism spectrum disorders (ASDs) remain mostly unknown. Previous studies have implicated the transcription factor FOXP1 in both brain development and ASD pathophysiology. However, the specific molecular pathways both upstream of and downstream from FOXP1 are not fully understood. To elucidate the contribution of FOXP1-mediated signaling to brain development and, in particular, neocortical development, we generated forebrain-specific Foxp1 conditional knockout mice. We show that deletion of Foxp1 in the developing forebrain leads to impairments in neonatal vocalizations as well as neocortical cytoarchitectonic alterations via neuronal positioning and migration. Using a genomics approach, we identified the transcriptional networks regulated by Foxp1 in the developing neocortex and found that such networks are enriched for downstream targets involved in neurogenesis and neuronal migration. We also uncovered mechanistic insight into Foxp1 function by demonstrating that sumoylation of Foxp1 during embryonic brain development is necessary for mediating proper interactions between Foxp1 and the NuRD complex. Furthermore, we demonstrated that sumoylation of Foxp1 affects neuronal differentiation and migration in the developing neocortex. Together, these data provide critical mechanistic insights into the function of FOXP1 in the developing neocortex and may reveal molecular pathways at risk in ASD.

Malignant gliomas, including glioblastoma multiforme, constitute the most common and aggressive primary brain tumors in adults. The transcription factor signal transducer and activator of transcription 3 (STAT3) plays an essential role in glioblastoma pathogenesis downstream of the major oncogenic protein epidermal growth factor receptor variant III (EGFRvIII). However, the critical gene targets of STAT3 that mediate EGFRvIII-induced glial transformation have remained unknown. Here, we identify inducible nitric oxide synthase (iNOS) as a novel target gene of STAT3 in EGFRvIII-expressing mouse astrocytes. Endogenous STAT3 occupies the endogenous iNOS promoter and stimulates iNOS transcription in EGFRvIII-expressing astrocytes. STAT3 does not appear to control iNOS transcription in astrocytes deficient in the major glioblastoma tumor suppressor protein phosphatase and tensin homolog (PTEN), suggesting that STAT3 regulates iNOS transcription specifically in EGFRvIII-expressing astrocytes. Importantly, inhibition of iNOS by distinct approaches, including knockdown by RNA interference, reduces cell population growth and invasiveness of EGFRvIII-expressing astrocytes. In addition, upon iNOS knockdown or administration of a small-molecule inhibitor of iNOS, EGFRvIII-expressing astrocytes form smaller tumors in vivo. These findings suggest that inhibition of iNOS may have potential therapeutic value for EGFRvIII-activated brain tumors.

Circadian oscillations are generated via transcriptional-translational negative feedback loops. However, individual cells from fibroblast cell lines have heterogeneous rhythms, oscillating independently and with different period lengths. Here we showed that heterogeneity in circadian period is heritable and used a multi-omics approach to investigate underlying mechanisms. By examining large-scale phenotype-associated gene expression profiles in hundreds of mouse clonal cell lines, we identified and validated multiple novel candidate genes involved in circadian period determination in the absence of significant genomic variants. We also discovered differentially co-expressed gene networks that were functionally associated with period length. We further demonstrated that global differential DNA methylation bidirectionally regulated these same gene networks. Interestingly, we found that depletion of DNMT1 and DNMT3A had opposite effects on circadian period, suggesting non-redundant roles in circadian gene regulation. Together, our findings identify novel gene candidates involved in periodicity, and reveal DNA methylation as an important regulator of circadian periodicity.

Concussion is associated with a myriad of deleterious immediate and long-term consequences. Yet the molecular mechanisms and genetic targets promoting the selective vulnerability of different neural subtypes to dysfunction and degeneration remain unclear. Translating experimental models of blunt force trauma in C. elegans to concussion in mice, we identify a conserved neuroprotective mechanism in which reduction of mitochondrial electron flux through complex IV suppresses trauma-induced degeneration of the highly vulnerable dopaminergic neurons. Reducing cytochrome C oxidase function elevates mitochondrial-derived reactive oxygen species, which signal through the cytosolic hypoxia inducing transcription factor, Hif1a, to promote hyperphosphorylation and inactivation of the pyruvate dehydrogenase, PDHE1α. This critical enzyme initiates the Warburg shunt, which drives energetic reallocation from mitochondrial respiration to astrocyte-mediated glycolysis in a neuroprotective manner. These studies demonstrate a conserved process in which glycolytic preconditioning suppresses Parkinson-like hypersensitivity of dopaminergic neurons to trauma-induced degeneration via redox signaling and the Warburg effect.

The suprachiasmatic nucleus (SCN) is the master circadian pacemaker in mammals and is entrained by environmental light. However, the molecular basis of the response of the SCN to light is not fully understood. We used RNA/chromatin immunoprecipitation/single-nucleus sequencing with circadian behavioral assays to identify mouse SCN cell types and explore their responses to light. We identified three peptidergic cell types that responded to light in the SCN: arginine vasopressin (AVP), vasoactive intestinal peptide (VIP), and cholecystokinin (CCK). In each cell type, light-responsive subgroups were enriched for expression of neuronal Per-Arnt-Sim (PAS) domain protein 4 (NPAS4) target genes. Further, mice lacking Npas4 had a longer circadian period under constant conditions, a damped phase response curve to light, and reduced light-induced gene expression in the SCN. Our data indicate that NPAS4 is necessary for normal transcriptional responses to light in the SCN and critical for photic phase-shifting of circadian behavior.

During cortical development, human basal radial glial cells (bRGCs) are highly capable of sustained self-renewal and neurogenesis. Selective pressures on this cell type may have contributed to the evolution of the human neocortex, leading to an increase in cortical size. bRGCs have enriched expression for Forkhead Box P1 (FOXP1), a transcription factor implicated in neurodevelopmental disorders (NDDs) such as autism spectrum disorder. However, the cell type-specific roles of FOXP1 in bRGCs during cortical development remain unexplored. Here, we examine the requirement for FOXP1 gene expression regulation underlying the production of bRGCs using human brain organoids. We examine a developmental time point when FOXP1 expression is highest in the cortical progenitors, and the bRGCs, in particular, begin to actively produce neurons. With the loss of FOXP1, we show a reduction in the number of bRGCs, as well as reduced proliferation and differentiation of the remaining bRGCs, all of which lead to reduced numbers of excitatory cortical neurons over time. Using single-nuclei RNA sequencing and cell trajectory analysis, we uncover a role for FOXP1 in directing cortical progenitor proliferation and differentiation by regulating key signaling pathways related to neurogenesis and NDDs. Together, these results demonstrate that FOXP1 regulates human-specific features in early cortical development.

Long-range glutamatergic inputs from the cortex and thalamus are critical for motor and cognitive processing in the striatum. Transcription factors that orchestrate the development of these inputs are largely unknown. We investigated the role of a transcription factor and high-risk autism-associated gene, FOXP1, in the development of glutamatergic inputs onto spiny projection neurons (SPNs) in the striatum. We find that FOXP1 robustly drives the strengthening and maturation of glutamatergic input onto dopamine receptor 2-expressing SPNs (D2 SPNs) but has a comparatively milder effect on D1 SPNs. This process is cell-autonomous and is likely mediated through postnatal FOXP1 function at the postsynapse. We identified postsynaptic FOXP1-regulated transcripts as potential candidates for mediating these effects. Postnatal reinstatement of FOXP1 rescues electrophysiological deficits, reverses gene expression alterations resulting from embryonic deletion, and mitigates behavioral phenotypes. These results provide support for a possible therapeutic approach for individuals with FOXP1 syndrome.

DNA-based evolutionary comparisons of regulatory genomic elements enable insight into functional changes, overcoming tissue inaccessibility. Here, we harnessed adult and fetal cortex single-cell ATAC-seq datasets to uncover DNA substitutions specific to the human and human-ancestral lineages within apes. We found that fetal microglia identity is evolutionarily divergent in all lineages, whereas other cell types are conserved. Using multiomic datasets, we further identified genes linked to multiple lineage-divergent gene regulatory elements and implicated biological pathways associated with these divergent features. We also uncovered patterns of transcription factor binding site evolution across lineages and identified expansion of bHLH-PAS factor targets in human-hominin lineages, and MEF2 factor targets in the ape lineage. Finally, conserved features were more enriched in brain disease variants, whereas there was no distinct enrichment on the human lineage compared to its ancestral lineages. Our study identifies major evolutionary patterns in the human brain epigenome at cellular resolution.

How do epigenetic modifications change across species and how do these modifications affect evolution? These are fundamental questions at the forefront of our evolutionary epigenomic understanding. Our previous work investigated human and chimpanzee brain methylomes, but it was limited by the lack of outgroup data which is critical for comparative (epi)genomic studies. Here, we compared whole genome DNA methylation maps from brains of humans, chimpanzees and also rhesus macaques (outgroup) to elucidate DNA methylation changes during human brain evolution. Moreover, we validated that our approach is highly robust by further examining 38 human-specific DMRs using targeted deep genomic and bisulfite sequencing in an independent panel of 37 individuals from five primate species. Our unbiased genome-scan identified human brain differentially methylated regions (DMRs), irrespective of their associations with annotated genes. Remarkably, over half of the newly identified DMRs locate in intergenic regions or gene bodies. Nevertheless, their regulatory potential is on par with those of promoter DMRs. An intriguing observation is that DMRs are enriched in active chromatin loops, suggesting human-specific evolutionary remodeling at a higher-order chromatin structure. These findings indicate that there is substantial reprogramming of epigenomic landscapes during human brain evolution involving noncoding regions.

Mutations in the transcription factor Forkhead box p1 (FOXP1) are causative for neurodevelopmental disorders such as autism. However, the function of FOXP1 within the brain remains largely uncharacterized. Here, we identify the gene expression program regulated by FoxP1 in both human neural cells and patient-relevant heterozygous Foxp1 mouse brains. We demonstrate a role for FoxP1 in the transcriptional regulation of autism-related pathways as well as genes involved in neuronal activity. We show that Foxp1 regulates the excitability of striatal medium spiny neurons and that reduction of Foxp1 correlates with defects in ultrasonic vocalizations. Finally, we demonstrate that FoxP1 has an evolutionarily conserved role in regulating pathways involved in striatal neuron identity through gene expression studies in human neural progenitors with altered FOXP1 levels. These data support an integral role for FoxP1 in regulating signaling pathways vulnerable in autism and the specific regulation of striatal pathways important for vocal communication.

Copy-number variants of chromosome 16 region 16p11.2 are linked to neuropsychiatric disorders and are among the most prevalent in autism spectrum disorders. Of many 16p11.2 genes, Kctd13 has been implicated as a major driver of neurodevelopmental phenotypes. The function of KCTD13 in the mammalian brain, however, remains unknown. Here we delete the Kctd13 gene in mice and demonstrate reduced synaptic transmission. Reduced synaptic transmission correlates with increased levels of Ras homolog gene family, member A (RhoA), a KCTD13/CUL3 ubiquitin ligase substrate, and is reversed by RhoA inhibition, suggesting increased RhoA as an important mechanism. In contrast to a previous knockdown study, deletion of Kctd13 or kctd13 does not increase brain size or neurogenesis in mice or zebrafish, respectively. These findings implicate Kctd13 in the regulation of neuronal function relevant to neuropsychiatric disorders and clarify the role of Kctd13 in neurogenesis and brain size. Our data also reveal a potential role for RhoA as a therapeutic target in disorders associated with KCTD13 deletion.

Neuronal activity and gene expression in response to the loss of sleep can provide a window into the enigma of sleep function. Sleep loss is associated with brain differential gene expression, an increase in pyramidal cell mEPSC frequency and amplitude, and a characteristic rebound and resolution of slow wave sleep-slow wave activity (SWS-SWA). However, the molecular mechanism(s) mediating the sleep-loss response are not well understood. We show that sleep-loss regulates MEF2C phosphorylation, a key mechanism regulating MEF2C transcriptional activity, and that MEF2C function in postnatal excitatory forebrain neurons is required for the biological events in response to sleep loss in C57BL/6J mice. These include altered gene expression, the increase and recovery of synaptic strength, and the rebound and resolution of SWS-SWA, which implicate MEF2C as an essential regulator of sleep function.

Microdeletions of the MEF2C gene are linked to a syndromic form of autism termed MEF2C haploinsufficiency syndrome (MCHS). MEF2C hypofunction in neurons is presumed to underlie most of the symptoms of MCHS. However, it is unclear in which cell populations MEF2C functions to regulate neurotypical development.