Origins and Proliferative States of Human Oligodendrocyte Precursor Cells.

Human cerebral cortex size and complexity has increased greatly during evolution. While increased progenitor diversity and enhanced proliferative potential play important roles in human neurogenesis and gray matter expansion, the mechanisms of human oligodendrogenesis and white matter expansion remain largely unknown. Here, we identify EGFR-expressing "Pre-OPCs" that originate from outer radial glial cells (oRGs) and undergo mitotic somal translocation (MST) during division. oRG-derived Pre-OPCs provide an additional source of human cortical oligodendrocyte precursor cells (OPCs) and define a lineage trajectory. We further show that human OPCs undergo consecutive symmetric divisions to exponentially increase the progenitor pool size. Additionally, we find that the OPC-enriched gene, PCDH15, mediates daughter cell repulsion and facilitates proliferation. These findings indicate properties of OPC derivation, proliferation, and dispersion important for human white matter expansion and myelination.

Neurotoxic microglia promote TDP-43 proteinopathy in progranulin deficiency.

Aberrant aggregation of the RNA-binding protein TDP-43 in neurons is a hallmark of frontotemporal lobar degeneration caused by haploinsufficiency in the gene encoding progranulin1,2. However, the mechanism leading to TDP-43 proteinopathy remains unclear. Here we use single-nucleus RNA sequencing to show that progranulin deficiency promotes microglial transition from a homeostatic to a disease-specific state that causes endolysosomal dysfunction and neurodegeneration in mice. These defects persist even when Grn-/- microglia are cultured ex vivo. In addition, single-nucleus RNA sequencing reveals selective loss of excitatory neurons at disease end-stage, which is characterized by prominent nuclear and cytoplasmic TDP-43 granules and nuclear pore defects. Remarkably, conditioned media from Grn-/- microglia are sufficient to promote TDP-43 granule formation, nuclear pore defects and cell death in excitatory neurons via the complement activation pathway. Consistent with these results, deletion of the genes encoding C1qa and C3 mitigates microglial toxicity and rescues TDP-43 proteinopathy and neurodegeneration. These results uncover previously unappreciated contributions of chronic microglial toxicity to TDP-43 proteinopathy during neurodegeneration.

Neuronal vulnerability and multilineage diversity in multiple sclerosis.

Multiple sclerosis (MS) is a neuroinflammatory disease with a relapsing-remitting disease course at early stages, distinct lesion characteristics in cortical grey versus subcortical white matter and neurodegeneration at chronic stages. Here we used single-nucleus RNA sequencing to assess changes in expression in multiple cell lineages in MS lesions and validated the results using multiplex in situ hybridization. We found selective vulnerability and loss of excitatory CUX2-expressing projection neurons in upper-cortical layers underlying meningeal inflammation; such MS neuron populations exhibited upregulation of stress pathway genes and long non-coding RNAs. Signatures of stressed oligodendrocytes, reactive astrocytes and activated microglia mapped most strongly to the rim of MS plaques. Notably, single-nucleus RNA sequencing identified phagocytosing microglia and/or macrophages by their ingestion and perinuclear import of myelin transcripts, confirmed by functional mouse and human culture assays. Our findings indicate lineage- and region-specific transcriptomic changes associated with selective cortical neuron damage and glial activation contributing to progression of MS lesions.

Tropism of SARS-CoV-2 for human cortical astrocytes.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) readily infects a variety of cell types impacting the function of vital organ systems, with particularly severe impact on respiratory function. Neurological symptoms, which range in severity, accompany as many as one-third of COVID-19 cases, indicating a potential vulnerability of neural cell types. To assess whether human cortical cells can be directly infected by SARS-CoV-2, we utilized stem-cell-derived cortical organoids as well as primary human cortical tissue, both from developmental and adult stages. We find significant and predominant infection in cortical astrocytes in both primary tissue and organoid cultures, with minimal infection of other cortical populations. Infected and bystander astrocytes have a corresponding increase in inflammatory gene expression, reactivity characteristics, increased cytokine and growth factor signaling, and cellular stress. Although human cortical cells, particularly astrocytes, have no observable ACE2 expression, we find high levels of coronavirus coreceptors in infected astrocytes, including CD147 and DPP4. Decreasing coreceptor abundance and activity reduces overall infection rate, and increasing expression is sufficient to promote infection. Thus, we find tropism of SARS-CoV-2 for human astrocytes resulting in inflammatory gliosis-type injury that is dependent on coronavirus coreceptors.

Cross-regional homeostatic and reactive glial signatures in multiple sclerosis.

Multiple sclerosis (MS) is a multifocal and progressive inflammatory disease of the central nervous system (CNS). However, the compartmentalized pathology of the disease affecting various anatomical regions including gray and white matter and lack of appropriate disease models impede understanding of the disease. Utilizing single-nucleus RNA-sequencing and multiplex spatial RNA mapping, we generated an integrated transcriptomic map comprising leukocortical, cerebellar and spinal cord areas in normal and MS tissues that captures regional subtype diversity of various cell types with an emphasis on astrocytes and oligodendrocytes. While we found strong cross-regional diversity among glial subtypes in control tissue, regional signatures become more obscure in MS. This suggests that patterns of transcriptomic changes in MS are shared across regions and converge on specific pathways, especially those regulating cellular stress and immune activation. In addition, we found evidence that a subtype of white matter oligodendrocytes appearing across all three CNS regions adopt pro-remyelinating gene signatures in MS. In summary, our data suggest that cross-regional transcriptomic glial signatures overlap in MS, with different reactive glial cell types capable of either exacerbating or ameliorating pathology.

Cocaine alters Homer1 natural antisense transcript in the nucleus accumbens.

Natural antisense transcripts (NATs) are an abundant class of long noncoding RNAs that have recently been shown to be key regulators of chromatin dynamics and gene expression in nervous system development and neurological disorders. However, it is currently unclear if NAT-based mechanisms also play a role in drug-induced neuroadaptations. Aberrant regulation of gene expression is one critical factor underlying the long-lasting behavioral abnormalities that characterize substance use disorder, and it is possible that some drug-induced transcriptional responses are mediated, in part, by perturbations in NAT activity. To test this hypothesis, we used an automated algorithm that mines the NCBI AceView transcriptomics database to identify NAT overlapping genes linked to addiction. We found that 22% of the genes examined contain NATs and that expression of Homer1 natural antisense transcript (Homer1-AS) was altered in the nucleus accumbens (NAc) of mice 2h and 10days following repeated cocaine administration. In in vitro studies, depletion of Homer1-AS lead to an increase in the corresponding sense gene expression, indicating a potential regulatory mechanisms of Homer1 expression by its corresponding antisense transcript. Future in vivo studies are needed to definitely determine a role for Homer1-AS in cocaine-induced behavioral and molecular adaptations.

The long non-coding RNA FMR4 promotes proliferation of human neural precursor cells and epigenetic regulation of gene expression in trans.

Triplet repeat expansions in the Fragile X mental retardation 1 (FMR1) gene cause either intellectual disability and autism, or adult-onset neurodegeneration, with poorly understood variability in presentation. Previous studies have identified several long noncoding RNAs (lncRNAs) at the FMR1 locus, including FMR4. Similarly to FMR1, FMR4 is silenced by large-repeat expansions that result in enrichment of DNA and histone methylation within the shared promoter and repeat sequence, suggesting a possible role for this noncoding RNA in the pathophysiology of Fragile X. We therefore assessed the functional role of FMR4 to gain further insight into the molecular processes in Fragile X-associated disorders. Previous work showed that FMR4 does not exhibit cis-regulation of FMR1. Here, we found that FMR4 is a chromatin-associated transcript and, using genome-wide chromatin immunoprecipitation experiments, showed that FMR4 alters the chromatin state and the expression of several hundred genes in trans. Among the genes regulated by FMR4, we found enrichment for those involved in neural development and cellular proliferation. S-phase marker assays further demonstrated that FMR4 may promote cellular proliferation, rather than differentiation, of human neural precursor cells (hNPCs). By establishing this novel function for FMR4 in hNPCs, we lend support to existing evidence of the epigenetic involvement of lncRNA in nervous system development, and increase our understanding of the complex pathogenesis underlying neurological disorders associated with FMR1 repeat expansions.

The N-terminal Set-ß Protein Isoform Induces Neuronal Death.

Set-ß protein plays different roles in neurons, but the diversity of Set-ß neuronal isoforms and their functions have not been characterized. The expression and subcellular localization of Set-ß are altered in Alzheimer disease, cleavage of Set-ß leads to neuronal death after stroke, and the full-length Set-ß regulates retinal ganglion cell (RGC) and hippocampal neuron axon growth and regeneration in a subcellular localization-dependent manner. Here we used various biochemical approaches to investigate Set-ß isoforms and their role in the CNS, using the same type of neurons, RGCs, across studies. We found multiple alternatively spliced isoforms expressed from the Set locus in purified RGCs. Set transcripts containing the Set-ß-specific exon were the most highly expressed isoforms. We also identified a novel, alternatively spliced Set-ß transcript lacking the nuclear localization signal and demonstrated that the full-length (∼39-kDa) Set-ß is localized predominantly in the nucleus, whereas a shorter (∼25-kDa) Set-ß isoform is localized predominantly in the cytoplasm. Finally, we show that an N-terminal Set-ß cleavage product can induce neuronal death.

CANEapp: a user-friendly application for automated next generation transcriptomic data analysis.

BACKGROUND: Next generation sequencing (NGS) technologies are indispensable for molecular biology research, but data analysis represents the bottleneck in their application. Users need to be familiar with computer terminal commands, the Linux environment, and various software tools and scripts. Analysis workflows have to be optimized and experimentally validated to extract biologically meaningful data. Moreover, as larger datasets are being generated, their analysis requires use of high-performance servers. RESULTS: To address these needs, we developed CANEapp (application for Comprehensive automated Analysis of Next-generation sequencing Experiments), a unique suite that combines a Graphical User Interface (GUI) and an automated server-side analysis pipeline that is platform-independent, making it suitable for any server architecture. The GUI runs on a PC or Mac and seamlessly connects to the server to provide full GUI control of RNA-sequencing (RNA-seq) project analysis. The server-side analysis pipeline contains a framework that is implemented on a Linux server through completely automated installation of software components and reference files. Analysis with CANEapp is also fully automated and performs differential gene expression analysis and novel noncoding RNA discovery through alternative workflows (Cuffdiff and R packages edgeR and DESeq2). We compared CANEapp to other similar tools, and it significantly improves on previous developments. We experimentally validated CANEapp's performance by applying it to data derived from different experimental paradigms and confirming the results with quantitative real-time PCR (qRT-PCR). CANEapp adapts to any server architecture by effectively using available resources and thus handles large amounts of data efficiently. CANEapp performance has been experimentally validated on various biological datasets. CANEapp is available free of charge at http://psychiatry.med.miami.edu/research/laboratory-of-translational-rna-genomics/CANE-app . CONCLUSIONS: We believe that CANEapp will serve both biologists with no computational experience and bioinformaticians as a simple, timesaving but accurate and powerful tool to analyze large RNA-seq datasets and will provide foundations for future development of integrated and automated high-throughput genomics data analysis tools. Due to its inherently standardized pipeline and combination of automated analysis and platform-independence, CANEapp is an ideal for large-scale collaborative RNA-seq projects between different institutions and research groups.

A comparative transcriptomic analysis of astrocytes differentiation from human neural progenitor cells.

Astrocytes are a morphologically and functionally heterogeneous population of cells that play critical roles in neurodevelopment and in the regulation of central nervous system homeostasis. Studies of human astrocytes have been hampered by the lack of specific molecular markers and by the difficulties associated with purifying and culturing astrocytes from adult human brains. Human neural progenitor cells (NPCs) with self-renewal and multipotent properties represent an appealing model system to gain insight into the developmental genetics and function of human astrocytes, but a comprehensive molecular characterization that confirms the validity of this cellular system is still missing. Here we used an unbiased transcriptomic analysis to characterize in vitro culture of human NPCs and to define the gene expression programs activated during the differentiation of these cells into astrocytes using FBS or the combination of CNTF and BMP4. Our results demonstrate that in vitro cultures of human NPCs isolated during the gliogenic phase of neurodevelopment mainly consist of radial glial cells (RGCs) and glia-restricted progenitor cells. In these cells the combination of CNTF and BMP4 activates the JAK/STAT and SMAD signaling cascades, leading to the inhibition of oligodendrocytes lineage commitment and activation of astrocytes differentiation. On the other hand, FBS-derived astrocytes have properties of reactive astrocytes. Our work suggests that in vitro culture of human NPCs represents a valuable cellular system to study human disorders characterized by impairment of astrocytes development and function. Our datasets represent an important resource for researchers studying human astrocytes development and might set the basis for the discovery of novel human-specific astrocyte markers.

Proximity analysis of native proteomes reveals phenotypic modifiers in a mouse model of autism and related neurodevelopmental conditions.

One of the main drivers of autism spectrum disorder is risk alleles within hundreds of genes, which may interact within shared but unknown protein complexes. Here we develop a scalable genome-editing-mediated approach to target 14 high-confidence autism risk genes within the mouse brain for proximity-based endogenous proteomics, achieving the identification of high-specificity spatial proteomes. The resulting native proximity proteomes are enriched for human genes dysregulated in the brain of autistic individuals, and reveal proximity interactions between proteins from high-confidence risk genes with those of lower-confidence that may provide new avenues to prioritize genetic risk. Importantly, the datasets are enriched for shared cellular functions and genetic interactions that may underlie the condition. We test this notion by spatial proteomics and CRISPR-based regulation of expression in two autism models, demonstrating functional interactions that modulate mechanisms of their dysregulation. Together, these results reveal native proteome networks in vivo relevant to autism, providing new inroads for understanding and manipulating the cellular drivers underpinning its etiology.

Single molecule array measures of LRRK2 kinase activity in serum link Parkinson's disease severity to peripheral inflammation.

Background: LRRK2-targeting therapeutics that inhibit LRRK2 kinase activity have advanced to clinical trials in idiopathic Parkinson's disease (iPD). LRRK2 phosphorylates Rab10 on endolysosomes in phagocytic cells to promote some types of immunological responses. The identification of factors that regulate LRRK2-mediated Rab10 phosphorylation in iPD, and whether phosphorylated-Rab10 levels change in different disease states, or with disease progression, may provide insights into the role of Rab10 phosphorylation in iPD and help guide therapeutic strategies targeting this pathway. Methods: Capitalizing on past work demonstrating LRRK2 and phosphorylated-Rab10 interact on vesicles that can shed into biofluids, we developed and validated a high-throughput single-molecule array assay to measure extracellular pT73-Rab10. Ratios of pT73-Rab10 to total Rab10 measured in biobanked serum samples were compared between informative groups of transgenic mice, rats, and a deeply phenotyped cohort of iPD cases and controls. Multivariable and weighted correlation network analyses were used to identify genetic, transcriptomic, clinical, and demographic variables that predict the extracellular pT73-Rab10 to total Rab10 ratio. Results: pT73-Rab10 is absent in serum from Lrrk2 knockout mice but elevated by LRRK2 and VPS35 mutations, as well as SNCA expression. Bone-marrow transplantation experiments in mice show that serum pT73-Rab10 levels derive primarily from circulating immune cells. The extracellular ratio of pT73-Rab10 to total Rab10 is dynamic, increasing with inflammation and rapidly decreasing with LRRK2 kinase inhibition. The ratio of pT73-Rab10 to total Rab10 is elevated in iPD patients with greater motor dysfunction, irrespective of disease duration, age, sex, or the usage of PD-related or anti-inflammatory medications. pT73-Rab10 to total Rab10 ratios are associated with neutrophil activation, antigenic responses, and the suppression of platelet activation. Conclusions: The extracellular ratio of pT73-Rab10 to total Rab10 in serum is a novel pharmacodynamic biomarker for LRRK2-linked innate immune activation associated with disease severity in iPD. We propose that those iPD patients with higher serum pT73-Rab10 levels may benefit from LRRK2-targeting therapeutics to mitigate associated deleterious immunological responses.

Single molecule array measures of LRRK2 kinase activity in serum link Parkinson's disease severity to peripheral inflammation.

BACKGROUND: LRRK2-targeting therapeutics that inhibit LRRK2 kinase activity have advanced to clinical trials in idiopathic Parkinson's disease (iPD). LRRK2 phosphorylates Rab10 on endolysosomes in phagocytic cells to promote some types of immunological responses. The identification of factors that regulate LRRK2-mediated Rab10 phosphorylation in iPD, and whether phosphorylated-Rab10 levels change in different disease states, or with disease progression, may provide insights into the role of Rab10 phosphorylation in iPD and help guide therapeutic strategies targeting this pathway. METHODS: Capitalizing on past work demonstrating LRRK2 and phosphorylated-Rab10 interact on vesicles that can shed into biofluids, we developed and validated a high-throughput single-molecule array assay to measure extracellular pT73-Rab10. Ratios of pT73-Rab10 to total Rab10 measured in biobanked serum samples were compared between informative groups of transgenic mice, rats, and a deeply phenotyped cohort of iPD cases and controls. Multivariable and weighted correlation network analyses were used to identify genetic, transcriptomic, clinical, and demographic variables that predict the extracellular pT73-Rab10 to total Rab10 ratio. RESULTS: pT73-Rab10 is absent in serum from Lrrk2 knockout mice but elevated by LRRK2 and VPS35 mutations, as well as SNCA expression. Bone-marrow transplantation experiments in mice show that serum pT73-Rab10 levels derive primarily from circulating immune cells. The extracellular ratio of pT73-Rab10 to total Rab10 is dynamic, increasing with inflammation and rapidly decreasing with LRRK2 kinase inhibition. The ratio of pT73-Rab10 to total Rab10 is elevated in iPD patients with greater motor dysfunction, irrespective of disease duration, age, sex, or the usage of PD-related or anti-inflammatory medications. pT73-Rab10 to total Rab10 ratios are associated with neutrophil degranulation, antigenic responses, and suppressed platelet activation. CONCLUSIONS: The extracellular serum ratio of pT73-Rab10 to total Rab10 is a novel pharmacodynamic biomarker for LRRK2-linked innate immune activation associated with disease severity in iPD. We propose that those iPD patients with higher serum pT73-Rab10 levels may benefit from LRRK2-targeting therapeutics that mitigate associated deleterious immunological responses.

Calorie restriction alleviates the age-related decrease in neural progenitor cell division in the aging brain.

Production of new neurons from stem cells is important for cognitive function, and the reduction of neurogenesis in the aging brain may contribute to the accumulation of age-related cognitive deficits. Restriction of calorie intake and prolonged treatment with rapamycin have been shown to extend the lifespan of animals and delay the onset of the age-related decline in tissue and organ function. Using a reporter line in which neural stem and progenitor cells are marked by the expression of green fluorescent protein (GFP), we examined the effect of prolonged exposure to calorie restriction (CR) or rapamycin on hippocampal neural stem and progenitor cell proliferation in aging mice. We showed that CR increased the number of dividing cells in the dentate gyrus of female mice. The majority of these cells corresponded to nestin-GFP-expressing neural stem or progenitor cells; however, this increased proliferative activity of stem and progenitor cells did not result in a significant increase in the number of doublecortin-positive newborn neurons. Our results suggest that restricted calorie intake may increase the number of divisions that neural stem and progenitor cells undergo in the aging brain of females.

Single cell analysis of dup15q syndrome reveals developmental and postnatal molecular changes in autism.

Duplication 15q (dup15q) syndrome is the most common genetic cause of autism spectrum disorder (ASD). Due to a higher genetic and phenotypic homogeneity compared to idiopathic autism, dup15q syndrome provides a well-defined setting to investigate ASD mechanisms. Previous bulk gene expression studies identified shared molecular changes in ASD. However, how cell type specific changes compare across different autism subtypes and how they change during development is largely unknown. In this study, we used single cell and single nucleus mRNA sequencing of dup15q cortical organoids from patient iPSCs, as well as post-mortem patient brain samples. We find cell-type specific dysregulated programs that underlie dup15q pathogenesis, which we validate by spatial resolved transcriptomics using brain tissue samples. We find degraded identity and vulnerability of deep-layer neurons in fetal stage organoids and highlight increased molecular burden of postmortem upper-layer neurons implicated in synaptic signaling, a finding shared between idiopathic ASD and dup15q syndrome. Gene co-expression network analysis of organoid and postmortem excitatory neurons uncovers modules enriched with autism risk genes. Organoid developmental modules were involved in transcription regulation via chromatin remodeling, while postmortem modules were associated with synaptic transmission and plasticity. The findings reveal a shifting landscape of ASD cellular vulnerability during brain development.

Single-cell analysis of prenatal and postnatal human cortical development.

We analyzed >700,000 single-nucleus RNA sequencing profiles from 106 donors during prenatal and postnatal developmental stages and identified lineage-specific programs that underlie the development of specific subtypes of excitatory cortical neurons, interneurons, glial cell types, and brain vasculature. By leveraging single-nucleus chromatin accessibility data, we delineated enhancer gene regulatory networks and transcription factors that control commitment of specific cortical lineages. By intersecting our results with genetic risk factors for human brain diseases, we identified the cortical cell types and lineages most vulnerable to genetic insults of different brain disorders, especially autism. We find that lineage-specific gene expression programs up-regulated in female cells are especially enriched for the genetic risk factors of autism. Our study captures the molecular progression of cortical lineages across human development.

Astroglial toxicity promotes synaptic degeneration in the thalamocortical circuit in frontotemporal dementia with GRN mutations.

Mutations in the human progranulin (GRN) gene are a leading cause of frontotemporal lobar degeneration (FTLD). While previous studies implicate aberrant microglial activation as a disease-driving factor in neurodegeneration in the thalamocortical circuit in Grn-/- mice, the exact mechanism for neurodegeneration in FTLD-GRN remains unclear. By performing comparative single-cell transcriptomics in the thalamus and frontal cortex of Grn-/- mice and patients with FTLD-GRN, we have uncovered a highly conserved astroglial pathology characterized by upregulation of gap junction protein GJA1, water channel AQP4, and lipid-binding protein APOE, and downregulation of glutamate transporter SLC1A2 that promoted profound synaptic degeneration across the two species. This astroglial toxicity could be recapitulated in mouse astrocyte-neuron cocultures and by transplanting induced pluripotent stem cell-derived astrocytes to cortical organoids, where progranulin-deficient astrocytes promoted synaptic degeneration, neuronal stress, and TDP-43 proteinopathy. Together, these results reveal a previously unappreciated astroglial pathology as a potential key mechanism in neurodegeneration in FTLD-GRN.

Tropism of SARS-CoV-2 for Developing Human Cortical Astrocytes.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) readily infects a variety of cell types impacting the function of vital organ systems, with particularly severe impact on respiratory function. It proves fatal for one percent of those infected. Neurological symptoms, which range in severity, accompany a significant proportion of COVID-19 cases, indicating a potential vulnerability of neural cell types. To assess whether human cortical cells can be directly infected by SARS-CoV-2, we utilized primary human cortical tissue and stem cell-derived cortical organoids. We find significant and predominant infection in cortical astrocytes in both primary and organoid cultures, with minimal infection of other cortical populations. Infected astrocytes had a corresponding increase in reactivity characteristics, growth factor signaling, and cellular stress. Although human cortical cells, including astrocytes, have minimal ACE2 expression, we find high levels of alternative coronavirus receptors in infected astrocytes, including DPP4 and CD147. Inhibition of DPP4 reduced infection and decreased expression of the cell stress marker, ARCN1. We find tropism of SARS-CoV-2 for human astrocytes mediated by DPP4, resulting in reactive gliosis-type injury.

Cell-Type-Specific Analysis of Molecular Pathology in Autism Identifies Common Genes and Pathways Affected Across Neocortical Regions.

Despite its heterogeneity, autism is characterized by a defined behavioral phenotype, suggesting that the molecular pathology affects specific neural substrates to cause behavioral dysfunction. Previous studies identified genes dysregulated in autism cortex but did not address their cell-type specificity. Moreover, it is unknown whether there is a core of genes dysregulated across multiple neocortical regions. We applied RNA sequencing to postmortem brain tissue samples from autism patients and neurologically normal controls and combined our data with previously published datasets. We then identified genes, pathways, and alternative splicing events which are dysregulated in five neocortical regions in autism. To gain information about cell-type specificity of the dysregulated genes, we analyzed single-nuclei RNA sequencing data of adult human cortex and intersected cell-type-specific gene signatures with genes dysregulated in autism in specific cortical regions. We found that autism-associated gene expression changes across 4 frontal and temporal cortex regions converge on 27 genes related to immune response and enriched in human astrocytes, microglia, and brain endothelium. Shared splicing changes, however, are found in genes predominantly associated with synaptic function and adult interneurons and projection neurons. Finally, we demonstrate that regions of DNA differentially methylated in autism overlap genes associated with development and enriched in human cortical oligodendrocytes. Our study identifies signatures of autism molecular pathology shared across neocortical regions, as well as neural cell types enriched for common dysregulated genes, thus paving way for assessing cell-type-specific mechanisms of autism pathology.

Identification of amygdala-expressed genes associated with autism spectrum disorder.

BACKGROUND: Studies of individuals with autism spectrum disorder (ASD) have revealed a strong multigenic basis with the identification of hundreds of ASD susceptibility genes. ASD is characterized by social deficits and a range of other phenotypes, implicating complex genetics and involvement of a variety of brain regions. However, how mutations and mis-expression of select gene sets are associated with the behavioral components of ASD remains unknown. We reasoned that for genes to be associated with ASD core behaviors they must be: (1) expressed in brain regions relevant to ASD social behaviors and (2) expressed during the ASD susceptible window of brain development. METHODS: Focusing on the amygdala, a brain region whose dysfunction has been highly implicated in the social component of ASD, we mined publicly available gene expression databases to identify ASD-susceptibility genes expressed during human and mouse amygdala development. We found that a large cohort of known ASD susceptibility genes is expressed in the developing human and mouse amygdala. We further performed analysis of single-nucleus RNA-seq (snRNA-seq) data from microdissected amygdala tissue from five ASD and five control human postmortem brains ranging in age from 4 to 20 years to elucidate cell type specificity of amygdala-expressed genes and their dysregulation in ASD. RESULTS: Our analyses revealed that of the high-ranking ASD susceptibility genes, 80 are expressed in both human and mouse amygdala during fetal to early postnatal stages of development. Our human snRNA-seq analyses revealed cohorts of genes with altered expression in the ASD amygdala postnatally, especially within excitatory neurons, with dysregulated expression of seven genes predicted from our datamining pipeline. LIMITATIONS: We were limited by the ages for which we were able to obtain human tissue; therefore, the results from our datamining pipeline approach will require validation, to the extent possible, in human tissue from earlier developmental stages. CONCLUSIONS: Our pipeline narrows down the number of amygdala-expressed genes possibly involved in the social pathophysiology of ASD. Our human single-nucleus gene expression analyses revealed that ASD is characterized by changes in gene expression in specific cell types in the early postnatal amygdala.