Using deep learning to quantify neuronal activation from single-cell and spatial transcriptomic data.

Neuronal activity-dependent transcription directs molecular processes that regulate synaptic plasticity, brain circuit development, behavioral adaptation, and long-term memory. Single cell RNA-sequencing technologies (scRNAseq) are rapidly developing and allow for the interrogation of activity-dependent transcription at cellular resolution. Here, we present NEUROeSTIMator, a deep learning model that integrates transcriptomic signals to estimate neuronal activation in a way that we demonstrate is associated with Patch-seq electrophysiological features and that is robust against differences in species, cell type, and brain region. We demonstrate this method's ability to accurately detect neuronal activity in previously published studies of single cell activity-induced gene expression. Further, we applied our model in a spatial transcriptomic study to identify unique patterns of learning-induced activity across different brain regions in male mice. Altogether, our findings establish NEUROeSTIMator as a powerful and broadly applicable tool for measuring neuronal activation, whether as a critical covariate or a primary readout of interest.

Activity-induced gene expression in the human brain.

Activity-induced gene expression underlies synaptic plasticity and brain function. Here, using molecular sequencing techniques, we define activity-dependent transcriptomic and epigenomic changes at the tissue and single-cell level in the human brain following direct electrical stimulation of the anterior temporal lobe in patients undergoing neurosurgery. Genes related to transcriptional regulation and microglia-specific cytokine activity displayed the greatest induction pattern, revealing a precise molecular signature of neuronal activation in the human brain.

Mapping the spatial transcriptomic signature of the hippocampus during memory consolidation.

Memory consolidation involves discrete patterns of transcriptional events in the hippocampus. Despite the emergence of single-cell transcriptomic profiling techniques, mapping the transcriptomic signature across subregions of the hippocampus has remained challenging. Here, we utilized unbiased spatial sequencing to delineate transcriptome-wide gene expression changes across subregions of the dorsal hippocampus of male mice following learning. We find that each subregion of the hippocampus exhibits distinct yet overlapping transcriptomic signatures. The CA1 region exhibited increased expression of genes related to transcriptional regulation, while the DG showed upregulation of genes associated with protein folding. Importantly, our approach enabled us to define the transcriptomic signature of learning within two less-defined hippocampal subregions, CA1 stratum radiatum, and oriens. We demonstrated that CA1 subregion-specific expression of a transcription factor subfamily has a critical functional role in the consolidation of long-term memory. This work demonstrates the power of spatial molecular approaches to reveal simultaneous transcriptional events across the hippocampus during memory consolidation.

Mapping the spatial transcriptomic signature of the hippocampus during memory consolidation.

Memory consolidation involves discrete patterns of transcriptional events in the hippocampus. Despite the emergence of single-cell transcriptomic profiling techniques, defining learning-responsive gene expression across subregions of the hippocampus has remained challenging. Here, we utilized unbiased spatial sequencing to elucidate transcriptome-wide changes in gene expression in the hippocampus following learning, enabling us to define molecular signatures unique to each hippocampal subregion. We find that each subregion of the hippocampus exhibits distinct yet overlapping transcriptomic signatures. Although the CA1 region exhibited increased expression of genes related to transcriptional regulation, the DG showed upregulation of genes associated with protein folding. We demonstrate the functional relevance of subregion-specific gene expression by genetic manipulation of a transcription factor selectively in the CA1 hippocampal subregion, leading to long-term memory deficits. This work demonstrates the power of using spatial molecular approaches to reveal transcriptional events during memory consolidation.

Clinical autism subscales have common genetic liabilities that are heritable, pleiotropic, and generalizable to the general population.

The complexity of autism's phenotypic spectra is well-known, yet most genetic research uses case-control status as the target trait. It is undetermined if autistic symptom domain severity underlying this heterogeneity is heritable and pleiotropic with other psychiatric and behavior traits in the same manner as autism case-control status. In N = 6064 autistic children in the SPARK cohort, we investigated the common genetic properties of twelve subscales from three clinical autism instruments measuring autistic traits: the Social Communication Questionnaire (SCQ), the Repetitive Behavior Scale-Revised (RBS-R), and the Developmental Coordination Disorder Questionnaire (DCDQ). Educational attainment polygenic scores (PGS) were significantly negatively correlated with eleven subscales, while ADHD and major depression PGS were positively correlated with ten and eight of the autism subscales, respectively. Loneliness and neuroticism PGS were also positively correlated with many subscales. Significant PGS by sex interactions were found-surprisingly, the autism case-control PGS was negatively correlated in females and had no strong correlation in males. SNP-heritability of the DCDQ subscales ranged from 0.04 to 0.08, RBS-R subscales ranged from 0.09 to 0.24, and SCQ subscales ranged from 0 to 0.12. GWAS in SPARK followed by estimation of polygenic scores (PGS) in the typically-developing ABCD cohort (N = 5285), revealed significant associations of RBS-R subscale PGS with autism-related behavioral traits, with several subscale PGS more strongly correlated than the autism case-control PGS. Overall, our analyses suggest that the clinical autism subscale traits show variability in SNP-heritability, PGS associations, and significant PGS by sex interactions, underscoring the heterogeneity in autistic traits at a genetic level. Furthermore, of the three instruments investigated, the RBS-R shows the greatest evidence of genetic signal in both (1) autistic samples (greater heritability) and (2) general population samples (strongest PGS associations).

Endoplasmic reticulum chaperone genes encode effectors of long-term memory.

The mechanisms underlying memory loss associated with Alzheimer's disease and related dementias (ADRD) remain unclear, and no effective treatments exist. Fundamental studies have shown that a set of transcriptional regulatory proteins of the nuclear receptor 4a (Nr4a) family serve as molecular switches for long-term memory. Here, we show that Nr4a proteins regulate the transcription of genes encoding chaperones that localize to the endoplasmic reticulum (ER). These chaperones fold and traffic plasticity-related proteins to the cell surface during long-lasting forms of synaptic plasticity and memory. Dysregulation of Nr4a transcription factors and ER chaperones is linked to ADRD, and overexpressing Nr4a1 or the chaperone Hspa5 ameliorates long-term memory deficits in a tau-based mouse model of ADRD, pointing toward innovative therapeutic approaches for treating memory loss. Our findings establish a unique molecular concept underlying long-term memory and provide insights into the mechanistic basis of cognitive deficits in dementia.

p53-mediated neurodegeneration in the absence of the nuclear protein Akirin2.

Proper gene regulation is critical for both neuronal development and maintenance as the brain matures. We previously demonstrated that Akirin2, an essential nuclear protein that interacts with transcription factors and chromatin remodeling complexes, is required for the embryonic formation of the cerebral cortex. Here we show that Akirin2 plays a mechanistically distinct role in maintaining healthy neurons during cortical maturation. Restricting Akirin2 loss to excitatory cortical neurons resulted in progressive neurodegeneration via necroptosis and severe cortical atrophy with age. Comparing transcriptomes from Akirin2-null postnatal neurons and cortical progenitors revealed that targets of the tumor suppressor p53, a regulator of both proliferation and cell death encoded by Trp53, were consistently upregulated. Reduction of Trp53 rescued neurodegeneration in Akirin2-null neurons. These data: (1) implicate Akirin2 as a critical neuronal maintenance protein, (2) identify p53 pathways as mediators of Akirin2 functions, and (3) suggest Akirin2 dysfunction may be relevant to neurodegenerative diseases.

Altered hippocampal transcriptome dynamics following sleep deprivation.

Widespread sleep deprivation is a continuing public health problem in the United States and worldwide affecting adolescents and adults. Acute sleep deprivation results in decrements in spatial memory and cognitive impairments. The hippocampus is vulnerable to acute sleep deprivation with changes in gene expression, cell signaling, and protein synthesis. Sleep deprivation also has long lasting effects on memory and performance that persist after recovery sleep, as seen in behavioral studies from invertebrates to humans. Although previous research has shown that acute sleep deprivation impacts gene expression, the extent to which sleep deprivation affects gene regulation remains unknown. Using an unbiased deep RNA sequencing approach, we investigated the effects of acute sleep deprivation on gene expression in the hippocampus. We identified 1,146 genes that were significantly dysregulated following sleep deprivation with 507 genes upregulated and 639 genes downregulated, including protein coding genes and long non-coding RNAs not previously identified as impacted by sleep deprivation. Notably, genes significantly upregulated after sleep deprivation were associated with RNA splicing and the nucleus. In contrast, downregulated genes were associated with cell adhesion, dendritic localization, the synapse, and postsynaptic membrane. Furthermore, we found through independent experiments analyzing a subset of genes that three hours of recovery sleep following acute sleep deprivation was sufficient to normalize mRNA abundance for most genes, although exceptions occurred for some genes that may affect RNA splicing or transcription. These results clearly demonstrate that sleep deprivation differentially regulates gene expression on multiple transcriptomic levels to impact hippocampal function.

Genetic and morphological estimates of androgen exposure predict social deficits in multiple neurodevelopmental disorder cohorts.

BACKGROUND: Neurodevelopmental disorders (NDDs) such as autism spectrum disorder (ASD) display a strong male bias. Androgen exposure is profoundly increased in typical male development, but it also varies within the sexes, and previous work has sought to connect morphological proxies of androgen exposure, including digit ratio and facial morphology, to neurodevelopmental outcomes. The results of these studies have been mixed, and the relationships between androgen exposure and behavior remain unclear. METHODS: Here, we measured both digit ratio masculinity (DRM) and facial landmark masculinity (FLM) in the same neurodevelopmental cohort (N = 763) and compared these proxies of androgen exposure to clinical and parent-reported features as well as polygenic risk scores. RESULTS: We found that FLM was significantly associated with NDD diagnosis (ASD, ADHD, ID; all [Formula: see text]), while DRM was not. When testing for association with parent-reported problems, we found that both FLM and DRM were positively associated with concerns about social behavior ([Formula: see text], [Formula: see text]; [Formula: see text], [Formula: see text], respectively). Furthermore, we found evidence via polygenic risk scores (PRS) that DRM indexes masculinity via testosterone levels ([Formula: see text], [Formula: see text]), while FLM indexes masculinity through a negative relationship with sex hormone binding globulin (SHBG) levels ([Formula: see text], [Formula: see text]). Finally, using the SPARK cohort (N = 9419) we replicated the observed relationship between polygenic estimates of testosterone, SHBG, and social functioning ([Formula: see text], [Formula: see text], and [Formula: see text], [Formula: see text] for testosterone and SHBG, respectively). Remarkably, when considered over the extremes of each variable, these quantitative sex effects on social functioning were comparable to the effect of binary sex itself (binary male: [Formula: see text]; testosterone: [Formula: see text] from 0.1%-ile to 99.9%-ile; SHBG: [Formula: see text] from 0.1%-ile to 99.9%-ile). LIMITATIONS: In the devGenes and SPARK cohorts, our analyses rely on indirect, rather than direct measurement of androgens and related molecules. CONCLUSIONS: These findings and their replication in the large SPARK cohort lend support to the hypothesis that increasing net androgen exposure diminishes capacity for social functioning in both males and females.

The CBP KIX domain regulates long-term memory and circadian activity.

BACKGROUND: CREB-dependent transcription necessary for long-term memory is driven by interactions with CREB-binding protein (CBP), a multi-domain protein that binds numerous transcription factors potentially affecting expression of thousands of genes. Identifying specific domain functions for multi-domain proteins is essential to understand processes such as cognitive function and circadian clocks. We investigated the function of the CBP KIX domain in hippocampal memory and gene expression using CBPKIX/KIX mice with mutations that prevent phospho-CREB (Ser133) binding. RESULTS: We found that CBPKIX/KIX mice were impaired in long-term memory, but not learning acquisition or short-term memory for the Morris water maze. Using an unbiased analysis of gene expression in the dorsal hippocampus after training in the Morris water maze or contextual fear conditioning, we discovered dysregulation of CREB, CLOCK, and BMAL1 target genes and downregulation of circadian genes in CBPKIX/KIX mice. Given our finding that the CBP KIX domain was important for transcription of circadian genes, we profiled circadian activity and phase resetting in CBPKIX/KIX mice. CBPKIX/KIX mice exhibited delayed activity peaks after light offset and longer free-running periods in constant dark. Interestingly, CBPKIX/KIX mice displayed phase delays and advances in response to photic stimulation comparable to wildtype littermates. Thus, this work delineates site-specific regulation of the circadian clock by a multi-domain protein. CONCLUSIONS: These studies provide insight into the significance of the CBP KIX domain by defining targets of CBP transcriptional co-activation in memory and the role of the CBP KIX domain in vivo on circadian rhythms.

Novel and ultra-rare damaging variants in neuropeptide signaling are associated with disordered eating behaviors.

OBJECTIVE: Eating disorders develop through a combination of genetic vulnerability and environmental stress, however the genetic basis of this risk is unknown. METHODS: To understand the genetic basis of this risk, we performed whole exome sequencing on 93 unrelated individuals with eating disorders (38 restricted-eating and 55 binge-eating) to identify novel damaging variants. Candidate genes with an excessive burden of predicted damaging variants were then prioritized based upon an unbiased, data-driven bioinformatic analysis. One top candidate pathway was empirically tested for therapeutic potential in a mouse model of binge-like eating. RESULTS: An excessive burden of novel damaging variants was identified in 186 genes in the restricted-eating group and 245 genes in the binge-eating group. This list is significantly enriched (OR = 4.6, p<0.0001) for genes involved in neuropeptide/neurotrophic pathways implicated in appetite regulation, including neurotensin-, glucagon-like peptide 1- and BDNF-signaling. Administration of the glucagon-like peptide 1 receptor agonist exendin-4 significantly reduced food intake in a mouse model of 'binge-like' eating. CONCLUSIONS: These findings implicate ultra-rare and novel damaging variants in neuropeptide/neurotropic factor signaling pathways in the development of eating disorder behaviors and identify glucagon-like peptide 1-receptor agonists as a potential treatment for binge eating.

cerebroViz: an R package for anatomical visualization of spatiotemporal brain data.

Summary: Spatiotemporal transcriptomic profiling has provided valuable insight into the patterning of gene expression throughout the human brain from early fetal development to adulthood. When combined with prior knowledge of a disease's age at onset and region-specificity, these expression profiles have provided the necessary context to both strengthen putative gene-disease associations and infer new associations. While a wealth of spatiotemporal expression data exists, there are currently no tools available to visualize this data within the anatomical context of the brain, thus limiting the intuitive interpretation of many such findings. We present cerebroViz, an R package to map spatiotemporal brain data to vector graphic diagrams of the human brain. Our tool allows rapid generation of publication-quality figures that highlight spatiotemporal trends in the input data, while striking a balance between usability and customization. cerebroViz is generalizable to any data quantifiable at a brain-regional resolution and currently supports visualization of up to thirty regions of the brain found in databases such as BrainSpan, GTEx and Roadmap Epigenomics. Availability and Implementation: cerebroViz is freely available through GitHub ( https://github.com/ethanbahl/cerebroViz ). The tutorial is available at http://ethanbahl.github.io/cerebroViz/. Contacts: ethan-bahl@uiowa.edu or jacob-michaelson@uiowa.edu. Supplementary information: Supplementary data are available at Bioinformatics online.