Single-nucleus expression characterization of non-enhancing region of recurrent high-grade glioma.

Background: Non-enhancing (NE) infiltrating tumor cells beyond the contrast-enhancing (CE) bulk of tumor are potential propagators of recurrence after gross total resection of high-grade glioma. Methods: We leveraged single-nucleus RNA sequencing on 15 specimens from recurrent high-grade gliomas (n = 5) to compare prospectively identified biopsy specimens acquired from CE and NE regions. Additionally, 24 CE and 22 NE biopsies had immunohistochemical staining to validate RNA findings. Results: Tumor cells in NE regions are enriched in neural progenitor cell-like cellular states, while CE regions are enriched in mesenchymal-like states. NE glioma cells have similar proportions of proliferative and putative glioma stem cells relative to CE regions, without significant differences in % Ki-67 staining. Tumor cells in NE regions exhibit upregulation of genes previously associated with lower grade gliomas. Our findings in recurrent GBM paralleled some of the findings in a re-analysis of a dataset from primary GBM. Cell-, gene-, and pathway-level analyses of the tumor microenvironment in the NE region reveal relative downregulation of tumor-mediated neovascularization and cell-mediated immune response, but increased glioma-to-nonpathological cell interactions. Conclusions: This comprehensive analysis illustrates differing tumor and nontumor landscapes of CE and NE regions in high-grade gliomas, highlighting the NE region as an area harboring likely initiators of recurrence in a pro-tumor microenvironment and identifying possible targets for future design of NE-specific adjuvant therapy. These findings also support the aggressive approach to resection of tumor-bearing NE regions.

Gene networks and systems biology in Alzheimer's disease: Insights from multi-omics approaches.

Despite numerous studies in the field of dementia and Alzheimer's disease (AD), a comprehensive understanding of this devastating disease remains elusive. Bulk transcriptomics have provided insights into the underlying genetic factors at a high level. Subsequent technological advancements have focused on single-cell omics, encompassing techniques such as single-cell RNA sequencing and epigenomics, enabling the capture of RNA transcripts and chromatin states at a single cell or nucleus resolution. Furthermore, the emergence of spatial omics has allowed the study of gene responses in the vicinity of amyloid beta plaques or across various brain regions. With the vast amount of data generated, utilizing gene regulatory networks to comprehensively study this disease has become essential. This review delves into some techniques employed in the field of AD, explores the discoveries made using these techniques, and provides insights into the future of the field.

Genetic diversity promotes resilience in a mouse model of Alzheimer's disease.

INTRODUCTION: Alzheimer's disease (AD) is a neurodegenerative disorder with multifactorial etiology, including genetic factors that play a significant role in disease risk and resilience. However, the role of genetic diversity in preclinical AD studies has received limited attention. METHODS: We crossed five Collaborative Cross strains with 5xFAD C57BL/6J female mice to generate F1 mice with and without the 5xFAD transgene. Amyloid plaque pathology, microglial and astrocytic responses, neurofilament light chain levels, and gene expression were assessed at various ages. RESULTS: Genetic diversity significantly impacts AD-related pathology. Hybrid strains showed resistance to amyloid plaque formation and neuronal damage. Transcriptome diversity was maintained across ages and sexes, with observable strain-specific variations in AD-related phenotypes. Comparative gene expression analysis indicated correlations between mouse strains and human AD. DISCUSSION: Increasing genetic diversity promotes resilience to AD-related pathogenesis, relative to an inbred C57BL/6J background, reinforcing the importance of genetic diversity in uncovering resilience in the development of AD. HIGHLIGHTS: Genetic diversity's impact on AD in mice was explored. Diverse F1 mouse strains were used for AD study, via the Collaborative Cross. Strain-specific variations in AD pathology, glia, and transcription were found. Strains resilient to plaque formation and plasma neurofilament light chain (NfL) increases were identified. Correlations with human AD transcriptomics were observed.

BIN1K358R suppresses glial response to plaques in mouse model of Alzheimer's disease.

INTRODUCTION: The BIN1 coding variant rs138047593 (K358R) is linked to Late-Onset Alzheimer's Disease (LOAD) via targeted exome sequencing. METHODS: To elucidate the functional consequences of this rare coding variant on brain amyloidosis and neuroinflammation, we generated BIN1K358R knock-in mice using CRISPR/Cas9 technology. These mice were subsequently bred with 5xFAD transgenic mice, which serve as a model for Alzheimer's pathology. RESULTS: The presence of the BIN1K358R variant leads to increased cerebral amyloid deposition, with a dampened response of astrocytes and oligodendrocytes, but not microglia, at both the cellular and transcriptional levels. This correlates with decreased neurofilament light chain in both plasma and brain tissue. Synaptic densities are significantly increased in both wild-type and 5xFAD backgrounds homozygous for the BIN1K358R variant. DISCUSSION: The BIN1 K358R variant modulates amyloid pathology in 5xFAD mice, attenuates the astrocytic and oligodendrocytic responses to amyloid plaques, decreases damage markers, and elevates synaptic densities. HIGHLIGHTS: BIN1 rs138047593 (K358R) coding variant is associated with increased risk of LOAD. BIN1 K358R variant increases amyloid plaque load in 12-month-old 5xFAD mice. BIN1 K358R variant dampens astrocytic and oligodendrocytic response to plaques. BIN1 K358R variant decreases neuronal damage in 5xFAD mice. BIN1 K358R upregulates synaptic densities and modulates synaptic transmission.

Relapse to cocaine seeking is regulated by medial habenula NR4A2/NURR1 in mice.

Drugs of abuse can persistently change the reward circuit in ways that contribute to relapse behavior, partly via mechanisms that regulate chromatin structure and function. Nuclear orphan receptor subfamily4 groupA member2 (NR4A2, also known as NURR1) is an important effector of histone deacetylase 3 (HDAC3)-dependent mechanisms in persistent memory processes and is highly expressed in the medial habenula (MHb), a region that regulates nicotine-associated behaviors. Here, expressing the Nr4a2 dominant negative (Nurr2c) in the MHb blocks reinstatement of cocaine seeking in mice. We use single-nucleus transcriptomics to characterize the molecular cascade following Nr4a2 manipulation, revealing changes in transcriptional networks related to addiction, neuroplasticity, and GABAergic and glutamatergic signaling. The network controlled by NR4A2 is characterized using a transcription factor regulatory network inference algorithm. These results identify the MHb as a pivotal regulator of relapse behavior and demonstrate the importance of NR4A2 as a key mechanism driving the MHb component of relapse.

Spatial and single-nucleus transcriptomic analysis of genetic and sporadic forms of Alzheimer's Disease.

The pathogenesis of Alzheimer's disease (AD) depends on environmental and heritable factors, with remarkable differences evident between individuals at the molecular level. Here we present a transcriptomic survey of AD using spatial transcriptomics (ST) and single-nucleus RNA-seq in cortical samples from early-stage AD, late-stage AD, and AD in Down Syndrome (AD in DS) donors. Studying AD in DS provides an opportunity to enhance our understanding of the AD transcriptome, potentially bridging the gap between genetic mouse models and sporadic AD. Our analysis revealed spatial and cell-type specific changes in disease, with broad similarities in these changes between sAD and AD in DS. We performed additional ST experiments in a disease timecourse of 5xFAD and wildtype mice to facilitate cross-species comparisons. Finally, amyloid plaque and fibril imaging in the same tissue samples used for ST enabled us to directly link changes in gene expression with accumulation and spread of pathology.

hdWGCNA identifies co-expression networks in high-dimensional transcriptomics data.

Biological systems are immensely complex, organized into a multi-scale hierarchy of functional units based on tightly regulated interactions between distinct molecules, cells, organs, and organisms. While experimental methods enable transcriptome-wide measurements across millions of cells, popular bioinformatic tools do not support systems-level analysis. Here we present hdWGCNA, a comprehensive framework for analyzing co-expression networks in high-dimensional transcriptomics data such as single-cell and spatial RNA sequencing (RNA-seq). hdWGCNA provides functions for network inference, gene module identification, gene enrichment analysis, statistical tests, and data visualization. Beyond conventional single-cell RNA-seq, hdWGCNA is capable of performing isoform-level network analysis using long-read single-cell data. We showcase hdWGCNA using data from autism spectrum disorder and Alzheimer's disease brain samples, identifying disease-relevant co-expression network modules. hdWGCNA is directly compatible with Seurat, a widely used R package for single-cell and spatial transcriptomics analysis, and we demonstrate the scalability of hdWGCNA by analyzing a dataset containing nearly 1 million cells.

Integrin-Driven Axon Regeneration in the Spinal Cord Activates a Distinctive CNS Regeneration Program.

The peripheral branch of sensory dorsal root ganglion (DRG) neurons regenerates readily after injury unlike their central branch in the spinal cord. However, extensive regeneration and reconnection of sensory axons in the spinal cord can be driven by the expression of α9 integrin and its activator kindlin-1 (α9k1), which enable axons to interact with tenascin-C. To elucidate the mechanisms and downstream pathways affected by activated integrin expression and central regeneration, we conducted transcriptomic analyses of adult male rat DRG sensory neurons transduced with α9k1, and controls, with and without axotomy of the central branch. Expression of α9k1 without the central axotomy led to upregulation of a known PNS regeneration program, including many genes associated with peripheral nerve regeneration. Coupling α9k1 treatment with dorsal root axotomy led to extensive central axonal regeneration. In addition to the program upregulated by α9k1 expression, regeneration in the spinal cord led to expression of a distinctive CNS regeneration program, including genes associated with ubiquitination, autophagy, endoplasmic reticulum (ER), trafficking, and signaling. Pharmacological inhibition of these processes blocked the regeneration of axons from DRGs and human iPSC-derived sensory neurons, validating their causal contributions to sensory regeneration. This CNS regeneration-associated program showed little correlation with either embryonic development or PNS regeneration programs. Potential transcriptional drivers of this CNS program coupled to regeneration include Mef2a, Runx3, E2f4, and Yy1. Signaling from integrins primes sensory neurons for regeneration, but their axon growth in the CNS is associated with an additional distinctive program that differs from that involved in PNS regeneration.SIGNIFICANCE STATEMENT Restoration of neurologic function after spinal cord injury has yet to be achieved in human patients. To accomplish this, severed nerve fibers must be made to regenerate. Reconstruction of nerve pathways has not been possible, but recently, a method for stimulating long-distance axon regeneration of sensory fibers in rodents has been developed. This research uses profiling of messenger RNAs in the regenerating sensory neurons to discover which mechanisms are activated. This study shows that the regenerating neurons initiate a novel CNS regeneration program which includes molecular transport, autophagy, ubiquitination, and modulation of the endoplasmic reticulum (ER). The study identifies mechanisms that neurons need to activate to regenerate their nerve fibers.

A Trem2R47H mouse model without cryptic splicing drives age- and disease-dependent tissue damage and synaptic loss in response to plaques.

BACKGROUND: The TREM2 R47H variant is one of the strongest genetic risk factors for late-onset Alzheimer's Disease (AD). Unfortunately, many current Trem2 R47H mouse models are associated with cryptic mRNA splicing of the mutant allele that produces a confounding reduction in protein product. To overcome this issue, we developed the Trem2R47H NSS (Normal Splice Site) mouse model in which the Trem2 allele is expressed at a similar level to the wild-type Trem2 allele without evidence of cryptic splicing products. METHODS: Trem2R47H NSS mice were treated with the demyelinating agent cuprizone, or crossed with the 5xFAD mouse model of amyloidosis, to explore the impact of the TREM2 R47H variant on inflammatory responses to demyelination, plaque development, and the brain's response to plaques. RESULTS: Trem2R47H NSS mice display an appropriate inflammatory response to cuprizone challenge, and do not recapitulate the null allele in terms of impeded inflammatory responses to demyelination. Utilizing the 5xFAD mouse model, we report age- and disease-dependent changes in Trem2R47H NSS mice in response to development of AD-like pathology. At an early (4-month-old) disease stage, hemizygous 5xFAD/homozygous Trem2R47H NSS (5xFAD/Trem2R47H NSS) mice have reduced size and number of microglia that display impaired interaction with plaques compared to microglia in age-matched 5xFAD hemizygous controls. This is associated with a suppressed inflammatory response but increased dystrophic neurites and axonal damage as measured by plasma neurofilament light chain (NfL) level. Homozygosity for Trem2R47H NSS suppressed LTP deficits and loss of presynaptic puncta caused by the 5xFAD transgene array in 4-month-old mice. At a more advanced (12-month-old) disease stage 5xFAD/Trem2R47H NSS mice no longer display impaired plaque-microglia interaction or suppressed inflammatory gene expression, although NfL levels remain elevated, and a unique interferon-related gene expression signature is seen. Twelve-month old Trem2R47H NSS mice also display LTP deficits and postsynaptic loss. CONCLUSIONS: The Trem2R47H NSS mouse is a valuable model that can be used to investigate age-dependent effects of the AD-risk R47H mutation on TREM2 and microglial function including its effects on plaque development, microglial-plaque interaction, production of a unique interferon signature and associated tissue damage.

Single-nuclei transcriptome analysis of Huntington disease iPSC and mouse astrocytes implicates maturation and functional deficits.

Huntington disease (HD) is a neurodegenerative disorder caused by expanded CAG repeats in the huntingtin gene that alters cellular homeostasis, particularly in the striatum and cortex. Astrocyte signaling that establishes and maintains neuronal functions are often altered under pathological conditions. We performed single-nuclei RNA-sequencing on human HD patient-induced pluripotent stem cell (iPSC)-derived astrocytes and on striatal and cortical tissue from R6/2 HD mice to investigate high-resolution HD astrocyte cell state transitions. We observed altered maturation and glutamate signaling in HD human and mouse astrocytes. Human HD astrocytes also showed upregulated actin-mediated signaling, suggesting that some states may be cell-autonomous and human specific. In both species, astrogliogenesis transcription factors may drive HD astrocyte maturation deficits, which are supported by rescued climbing deficits in HD drosophila with NFIA knockdown. Thus, dysregulated HD astrocyte states may induce dysfunctional astrocytic properties, in part due to maturation deficits influenced by astrogliogenesis transcription factor dysregulation.

Huntington disease oligodendrocyte maturation deficits revealed by single-nucleus RNAseq are rescued by thiamine-biotin supplementation.

The complexity of affected brain regions and cell types is a challenge for Huntington's disease (HD) treatment. Here we use single nucleus RNA sequencing to investigate molecular pathology in the cortex and striatum from R6/2 mice and human HD post-mortem tissue. We identify cell type-specific and -agnostic signatures suggesting oligodendrocytes (OLs) and oligodendrocyte precursors (OPCs) are arrested in intermediate maturation states. OL-lineage regulators OLIG1 and OLIG2 are negatively correlated with CAG length in human OPCs, and ATACseq analysis of HD mouse NeuN-negative cells shows decreased accessibility regulated by OL maturation genes. The data implicates glucose and lipid metabolism in abnormal cell maturation and identify PRKCE and Thiamine Pyrophosphokinase 1 (TPK1) as central genes. Thiamine/biotin treatment of R6/1 HD mice to compensate for TPK1 dysregulation restores OL maturation and rescues neuronal pathology. Our insights into HD OL pathology spans multiple brain regions and link OL maturation deficits to abnormal thiamine metabolism.

Broad transcriptomic dysregulation occurs across the cerebral cortex in ASD.

Neuropsychiatric disorders classically lack defining brain pathologies, but recent work has demonstrated dysregulation at the molecular level, characterized by transcriptomic and epigenetic alterations1-3. In autism spectrum disorder (ASD), this molecular pathology involves the upregulation of microglial, astrocyte and neural-immune genes, the downregulation of synaptic genes, and attenuation of gene-expression gradients in cortex1,2,4-6. However, whether these changes are limited to cortical association regions or are more widespread remains unknown. To address this issue, we performed RNA-sequencing analysis of 725 brain samples spanning 11 cortical areas from 112 post-mortem samples from individuals with ASD and neurotypical controls. We find widespread transcriptomic changes across the cortex in ASD, exhibiting an anterior-to-posterior gradient, with the greatest differences in primary visual cortex, coincident with an attenuation of the typical transcriptomic differences between cortical regions. Single-nucleus RNA-sequencing and methylation profiling demonstrate that this robust molecular signature reflects changes in cell-type-specific gene expression, particularly affecting excitatory neurons and glia. Both rare and common ASD-associated genetic variation converge within a downregulated co-expression module involving synaptic signalling, and common variation alone is enriched within a module of upregulated protein chaperone genes. These results highlight widespread molecular changes across the cerebral cortex in ASD, extending beyond association cortex to broadly involve primary sensory regions.

Pathway-based approach reveals differential sensitivity to E2F1 inhibition in glioblastoma.

Analysis of tumor gene expression is an important approach for the classification and identification of therapeutic vulnerabilities. However, targeting glioblastoma (GBM) based on molecular subtyping has not yet translated into successful therapies. Here, we present an integrative approach based on molecular pathways to expose new potentially actionable targets. We used gene set enrichment analysis (GSEA) to conduct an unsupervised clustering analysis to condense the gene expression data from bulk patient samples and patient-derived gliomasphere lines into new gene signatures. We identified key targets that are predicted to be differentially activated between tumors and were functionally validated in a library of gliomasphere cultures. Resultant cluster-specific gene signatures associated not only with hallmarks of cell cycle and stemness gene expression, but also with cell-type specific markers and different cellular states of GBM. Several upstream regulators, such as PIK3R1 and EBF1 were differentially enriched in cells bearing stem cell like signatures and bear further investigation. We identified the transcription factor E2F1 as a key regulator of tumor cell proliferation and self-renewal in only a subset of gliomasphere cultures predicted to be E2F1 signaling dependent. Our in vivo work also validated the functional significance of E2F1 in tumor formation capacity in the predicted samples. E2F1 inhibition also differentially sensitized E2F1-dependent gliomasphere cultures to radiation treatment. Our findings indicate that this novel approach exploring cancer pathways highlights key therapeutic vulnerabilities for targeting GBM.

Molecular signatures underlying neurofibrillary tangle susceptibility in Alzheimer's disease.

Tau aggregation in neurofibrillary tangles (NFTs) is closely associated with neurodegeneration and cognitive decline in Alzheimer's disease (AD). However, the molecular signatures that distinguish between aggregation-prone and aggregation-resistant cell states are unknown. We developed methods for the high-throughput isolation and transcriptome profiling of single somas with NFTs from the human AD brain, quantified the susceptibility of 20 neocortical subtypes for NFT formation and death, and identified both shared and cell-type-specific signatures. NFT-bearing neurons shared a marked upregulation of synaptic transmission-related genes, including a core set of 63 genes enriched for synaptic vesicle cycling. Oxidative phosphorylation and mitochondrial dysfunction were highly cell-type dependent. Apoptosis was only modestly enriched, and the susceptibilities of NFT-bearing and NFT-free neurons for death were highly similar. Our analysis suggests that NFTs represent cell-type-specific responses to stress and synaptic dysfunction. We provide a resource for biomarker discovery and the investigation of tau-dependent and tau-independent mechanisms of neurodegeneration.

Protocol for single-nucleus ATAC sequencing and bioinformatic analysis in frozen human brain tissue.

Single-nucleus ATAC sequencing (snATAC-seq) employs a hyperactive Tn5 transposase to gain precise information about the cis-regulatory elements in specific cell types. However, the standard protocol of snATAC-seq is not optimized for all tissues, including the brain. Here, we present a modified protocol for single-nuclei isolation from postmortem frozen human brain tissue, followed by snATAC-seq library preparation and sequencing. We also describe an integrated bioinformatics analysis pipeline using an R package (ArchRtoSignac) to robustly analyze snATAC-seq data. For complete details on the use and execution of this protocol, please refer to Morabito et al. (2021).

Absence of microglia promotes diverse pathologies and early lethality in Alzheimer's disease mice.

Microglia are strongly implicated in the development and progression of Alzheimer's disease (AD), yet their impact on pathology and lifespan remains unclear. Here we utilize a CSF1R hypomorphic mouse to generate a model of AD that genetically lacks microglia. The resulting microglial-deficient mice exhibit a profound shift from parenchymal amyloid plaques to cerebral amyloid angiopathy (CAA), which is accompanied by numerous transcriptional changes, greatly increased brain calcification and hemorrhages, and premature lethality. Remarkably, a single injection of wild-type microglia into adult mice repopulates the microglial niche and prevents each of these pathological changes. Taken together, these results indicate the protective functions of microglia in reducing CAA, blood-brain barrier dysfunction, and brain calcification. To further understand the clinical implications of these findings, human AD tissue and iPSC-microglia were examined, providing evidence that microglia phagocytose calcium crystals, and this process is impaired by loss of the AD risk gene, TREM2.

Cortical diurnal rhythms remain intact with microglial depletion.

Microglia are subject to change in tandem with the endogenously generated biological oscillations known as our circadian rhythm. Studies have shown microglia harbor an intrinsic molecular clock which regulates diurnal changes in morphology and influences inflammatory responses. In the adult brain, microglia play an important role in the regulation of condensed extracellular matrix structures called perineuronal nets (PNNs), and it has been suggested that PNNs are also regulated in a circadian and diurnal manner. We sought to determine whether microglia mediate the diurnal regulation of PNNs via CSF1R inhibitor dependent microglial depletion in C57BL/6J mice, and how the absence of microglia might affect cortical diurnal gene expression rhythms. While we observe diurnal differences in microglial morphology, where microglia are most ramified at the onset of the dark phase, we do not find diurnal differences in PNN intensity. However, PNN intensity increases across many brain regions in the absence of microglia, supporting a role for microglia in the regulation of PNNs. Here, we also show that cortical diurnal gene expression rhythms are intact, with no cycling gene changes without microglia. These findings demonstrate a role for microglia in the maintenance of PNNs, but not in the maintenance of diurnal rhythms.

Systems biology approaches to unravel the molecular and genetic architecture of Alzheimer's disease and related tauopathies.

Over the years, genetic studies have identified multiple genetic risk variants associated with neurodegenerative disorders and helped reveal new biological pathways and genes of interest. However, genetic risk variants commonly reside in non-coding regions and may regulate distant genes rather than the nearest gene, as well as a gene's interaction partners in biological networks. Systems biology and functional genomics approaches provide the framework to unravel the functional significance of genetic risk variants in disease. In this review, we summarize the genetic and transcriptomic studies of Alzheimer's disease and related tauopathies and focus on the advantages of performing systems-level analyses to interrogate the biological pathways underlying neurodegeneration. Finally, we highlight new avenues of multi-omics analysis with single-cell approaches, which provide unparalleled opportunities to systematically explore cellular heterogeneity, and present an example of how to integrate publicly available single-cell datasets. Systems-level analysis has illuminated the function of many disease risk genes, but much work remains to study tauopathies and to understand spatiotemporal gene expression changes of specific cell types.

Microglial dyshomeostasis drives perineuronal net and synaptic loss in a CSF1R+/- mouse model of ALSP, which can be rescued via CSF1R inhibitors.

Adult-onset leukoencephalopathy with axonal spheroids and pigmented glia is an autosomal dominant neurodegenerative disease caused by mutations in colony-stimulating factor 1 receptor (CSF1R). We sought to identify the role of microglial CSF1R haploinsufficiency in mediating pathogenesis. Using an inducible Cx3cr1 CreERT2/+-Csf1r +/fl system, we found that postdevelopmental, microglia-specific Csf1r haploinsufficiency resulted in reduced expression of homeostatic microglial markers. This was associated with loss of presynaptic surrogates and the extracellular matrix (ECM) structure perineuronal nets. Similar phenotypes were observed in constitutive global Csf1r haploinsufficient mice and could be reversed/prevented by microglia elimination in adulthood. As microglial elimination is unlikely to be clinically feasible for extended durations, we treated adult CSF1R+/- mice at different disease stages with a microglia-modulating dose of the CSF1R inhibitor PLX5622, which prevented microglial dyshomeostasis along with synaptic- and ECM-related deficits. These data highlight microglial dyshomeostasis as a driver of pathogenesis and show that CSF1R inhibition can mitigate these phenotypes.