Targeting EGLN2/PHD1 protects motor neurons and normalizes the astrocytic interferon response.

Neuroinflammation and dysregulated energy metabolism are linked to motor neuron degeneration in amyotrophic lateral sclerosis (ALS). The egl-9 family hypoxia-inducible factor (EGLN) enzymes, also known as prolyl hydroxylase domain (PHD) enzymes, are metabolic sensors regulating cellular inflammation and metabolism. Using an oligonucleotide-based and a genetic approach, we showed that the downregulation of Egln2 protected motor neurons and mitigated the ALS phenotype in two zebrafish models and a mouse model of ALS. Single-nucleus RNA sequencing of the murine spinal cord revealed that the loss of EGLN2 induced an astrocyte-specific downregulation of interferon-stimulated genes, mediated via the stimulator of interferon genes (STING) protein. In addition, we found that the genetic deletion of EGLN2 restored this interferon response in patient induced pluripotent stem cell (iPSC)-derived astrocytes, confirming the link between EGLN2 and astrocytic interferon signaling. In conclusion, we identified EGLN2 as a motor neuron protective target normalizing the astrocytic interferon-dependent inflammatory axis in vivo, as well as in patient-derived cells.

Nova-ST: Nano-patterned ultra-dense platform for spatial transcriptomics.

Spatial transcriptomics workflows using barcoded capture arrays are commonly used for resolving gene expression in tissues. However, existing techniques are either limited by capture array density or are cost prohibitive for large-scale atlasing. We present Nova-ST, a dense nano-patterned spatial transcriptomics technique derived from randomly barcoded Illumina sequencing flow cells. Nova-ST enables customized, low-cost, flexible, and high-resolution spatial profiling of large tissue sections. Benchmarking on mouse brain sections demonstrates significantly higher sensitivity compared to existing methods at a reduced cost.

PARK7/DJ-1 deficiency impairs microglial activation in response to LPS-induced inflammation.

BACKGROUND: Specific microglia responses are thought to contribute to the development and progression of neurodegenerative diseases, including Parkinson's disease (PD). However, the phenotypic acquisition of microglial cells and their role during the underlying neuroinflammatory processes remain largely elusive. Here, according to the multiple-hit hypothesis, which stipulates that PD etiology is determined by a combination of genetics and various environmental risk factors, we investigate microglial transcriptional programs and morphological adaptations under PARK7/DJ-1 deficiency, a genetic cause of PD, during lipopolysaccharide (LPS)-induced inflammation. METHODS: Using a combination of single-cell RNA-sequencing, bulk RNA-sequencing, multicolor flow cytometry and immunofluorescence analyses, we comprehensively compared microglial cell phenotypic characteristics in PARK7/DJ-1 knock-out (KO) with wildtype littermate mice following 6- or 24-h intraperitoneal injection with LPS. For translational perspectives, we conducted corresponding analyses in human PARK7/DJ-1 mutant induced pluripotent stem cell (iPSC)-derived microglia and murine bone marrow-derived macrophages (BMDMs). RESULTS: By excluding the contribution of other immune brain resident and peripheral cells, we show that microglia acutely isolated from PARK7/DJ-1 KO mice display a distinct phenotype, specially related to type II interferon and DNA damage response signaling, when compared with wildtype microglia, in response to LPS. We also detected discrete signatures in human PARK7/DJ-1 mutant iPSC-derived microglia and BMDMs from PARK7/DJ-1 KO mice. These specific transcriptional signatures were reflected at the morphological level, with microglia in LPS-treated PARK7/DJ-1 KO mice showing a less amoeboid cell shape compared to wildtype mice, both at 6 and 24 h after acute inflammation, as also observed in BMDMs. CONCLUSIONS: Taken together, our results show that, under inflammatory conditions, PARK7/DJ-1 deficiency skews microglia towards a distinct phenotype characterized by downregulation of genes involved in type II interferon signaling and a less prominent amoeboid morphology compared to wildtype microglia. These findings suggest that the underlying oxidative stress associated with the lack of PARK7/DJ-1 affects microglia neuroinflammatory responses, which may play a causative role in PD onset and progression.

Glioblastoma-instructed microglia transition to heterogeneous phenotypic states with phagocytic and dendritic cell-like features in patient tumors and patient-derived orthotopic xenografts.

BACKGROUND: A major contributing factor to glioblastoma (GBM) development and progression is its ability to evade the immune system by creating an immune-suppressive environment, where GBM-associated myeloid cells, including resident microglia and peripheral monocyte-derived macrophages, play critical pro-tumoral roles. However, it is unclear whether recruited myeloid cells are phenotypically and functionally identical in GBM patients and whether this heterogeneity is recapitulated in patient-derived orthotopic xenografts (PDOXs). A thorough understanding of the GBM ecosystem and its recapitulation in preclinical models is currently missing, leading to inaccurate results and failures of clinical trials. METHODS: Here, we report systematic characterization of the tumor microenvironment (TME) in GBM PDOXs and patient tumors at the single-cell and spatial levels. We applied single-cell RNA sequencing, spatial transcriptomics, multicolor flow cytometry, immunohistochemistry, and functional studies to examine the heterogeneous TME instructed by GBM cells. GBM PDOXs representing different tumor phenotypes were compared to glioma mouse GL261 syngeneic model and patient tumors. RESULTS: We show that GBM tumor cells reciprocally interact with host cells to create a GBM patient-specific TME in PDOXs. We detected the most prominent transcriptomic adaptations in myeloid cells, with brain-resident microglia representing the main population in the cellular tumor, while peripheral-derived myeloid cells infiltrated the brain at sites of blood-brain barrier disruption. More specifically, we show that GBM-educated microglia undergo transition to diverse phenotypic states across distinct GBM landscapes and tumor niches. GBM-educated microglia subsets display phagocytic and dendritic cell-like gene expression programs. Additionally, we found novel microglial states expressing cell cycle programs, astrocytic or endothelial markers. Lastly, we show that temozolomide treatment leads to transcriptomic plasticity and altered crosstalk between GBM tumor cells and adjacent TME components. CONCLUSIONS: Our data provide novel insights into the phenotypic adaptation of the heterogeneous TME instructed by GBM tumors. We show the key role of microglial phenotypic states in supporting GBM tumor growth and response to treatment. Our data place PDOXs as relevant models to assess the functionality of the TME and changes in the GBM ecosystem upon treatment.

Xenografted human microglia display diverse transcriptomic states in response to Alzheimer's disease-related amyloid-ß pathology.

Microglia are central players in Alzheimer's disease pathology but analyzing microglial states in human brain samples is challenging due to genetic diversity, postmortem delay and admixture of pathologies. To circumvent these issues, here we generated 138,577 single-cell expression profiles of human stem cell-derived microglia xenotransplanted in the brain of the AppNL-G-F model of amyloid pathology and wild-type controls. Xenografted human microglia adopt a disease-associated profile similar to that seen in mouse microglia, but display a more pronounced human leukocyte antigen or HLA state, likely related to antigen presentation in response to amyloid plaques. The human microglial response also involves a pro-inflammatory cytokine/chemokine cytokine response microglia or CRM response to oligomeric Aß oligomers. Genetic deletion of TREM2 or APOE as well as APOE polymorphisms and TREM2R47H expression in the transplanted microglia modulate these responses differentially. The expression of other Alzheimer's disease risk genes is differentially regulated across the distinct cell states elicited in response to amyloid pathology. Thus, we have identified multiple transcriptomic cell states adopted by human microglia in a multipronged response to Alzheimer's disease-related pathology, which should be taken into account in translational studies.

Unravelling cell type-specific responses to Parkinson's Disease at single cell resolution.

Parkinson's Disease (PD) is the second most common neurodegenerative disorder. The pathological hallmark of PD is loss of dopaminergic neurons and the presence of aggregated α-synuclein, primarily in the substantia nigra pars compacta (SNpc) of the midbrain. However, the molecular mechanisms that underlie the pathology in different cell types is not currently understood. Here, we present a single nucleus transcriptome analysis of human post-mortem SNpc obtained from 15 sporadic Parkinson's Disease (PD) cases and 14 Controls. Our dataset comprises ∼84K nuclei, representing all major cell types of the brain, allowing us to obtain a transcriptome-level characterization of these cell types. Importantly, we identify multiple subpopulations for each cell type and describe specific gene sets that provide insights into the differing roles of these subpopulations. Our findings reveal a significant decrease in neuronal cells in PD samples, accompanied by an increase in glial cells and T cells. Subpopulation analyses demonstrate a significant depletion of tyrosine hydroxylase (TH) enriched astrocyte, microglia and oligodendrocyte populations in PD samples, as well as TH enriched neurons, which are also depleted. Moreover, marker gene analysis of the depleted subpopulations identified 28 overlapping genes, including those associated with dopamine metabolism (e.g., ALDH1A1, SLC6A3 & SLC18A2). Overall, our study provides a valuable resource for understanding the molecular mechanisms involved in dopaminergic neuron degeneration and glial responses in PD, highlighting the existence of novel subpopulations and cell type-specific gene sets.

Single-cell transcriptomics reveals that glial cells integrate homeostatic and circadian processes to drive sleep-wake cycles.

The sleep-wake cycle is determined by circadian and sleep homeostatic processes. However, the molecular impact of these processes and their interaction in different brain cell populations are unknown. To fill this gap, we profiled the single-cell transcriptome of adult Drosophila brains across the sleep-wake cycle and four circadian times. We show cell type-specific transcriptomic changes, with glia displaying the largest variation. Glia are also among the few cell types whose gene expression correlates with both sleep homeostat and circadian clock. The sleep-wake cycle and sleep drive level affect the expression of clock gene regulators in glia, and disrupting clock genes specifically in glia impairs homeostatic sleep rebound after sleep deprivation. These findings provide a comprehensive view of the effects of sleep homeostatic and circadian processes on distinct cell types in an entire animal brain and reveal glia as an interaction site of these two processes to determine sleep-wake dynamics.

Single-cell spatial multi-omics and deep learning dissect enhancer-driven gene regulatory networks in liver zonation.

In the mammalian liver, hepatocytes exhibit diverse metabolic and functional profiles based on their location within the liver lobule. However, it is unclear whether this spatial variation, called zonation, is governed by a well-defined gene regulatory code. Here, using a combination of single-cell multiomics, spatial omics, massively parallel reporter assays and deep learning, we mapped enhancer-gene regulatory networks across mouse liver cell types. We found that zonation affects gene expression and chromatin accessibility in hepatocytes, among other cell types. These states are driven by the repressors TCF7L1 and TBX3, alongside other core hepatocyte transcription factors, such as HNF4A, CEBPA, FOXA1 and ONECUT1. To examine the architecture of the enhancers driving these cell states, we trained a hierarchical deep learning model called DeepLiver. Our study provides a multimodal understanding of the regulatory code underlying hepatocyte identity and their zonation state that can be used to engineer enhancers with specific activity levels and zonation patterns.

The γ-secretase substrate proteome and its role in cell signaling regulation.

γ-Secretases mediate the regulated intramembrane proteolysis (RIP) of more than 150 integral membrane proteins. We developed an unbiased γ-secretase substrate identification (G-SECSI) method to study to what extent these proteins are processed in parallel. We demonstrate here parallel processing of at least 85 membrane proteins in human microglia in steady-state cell culture conditions. Pharmacological inhibition of γ-secretase caused substantial changes of human microglial transcriptomes, including the expression of genes related to the disease-associated microglia (DAM) response described in Alzheimer disease (AD). While the overall effects of γ-secretase deficiency on transcriptomic cell states remained limited in control conditions, exposure of mouse microglia to AD-inducing amyloid plaques strongly blocked their capacity to mount this putatively protective DAM cell state. We conclude that γ-secretase serves as a critical signaling hub integrating the effects of multiple extracellular stimuli into the overall transcriptome of the cell.

Systematic benchmarking of single-cell ATAC-sequencing protocols.

Single-cell assay for transposase-accessible chromatin by sequencing (scATAC-seq) has emerged as a powerful tool for dissecting regulatory landscapes and cellular heterogeneity. However, an exploration of systemic biases among scATAC-seq technologies has remained absent. In this study, we benchmark the performance of eight scATAC-seq methods across 47 experiments using human peripheral blood mononuclear cells (PBMCs) as a reference sample and develop PUMATAC, a universal preprocessing pipeline, to handle the various sequencing data formats. Our analyses reveal significant differences in sequencing library complexity and tagmentation specificity, which impact cell-type annotation, genotype demultiplexing, peak calling, differential region accessibility and transcription factor motif enrichment. Our findings underscore the importance of sample extraction, method selection, data processing and total cost of experiments, offering valuable guidance for future research. Finally, our data and analysis pipeline encompasses 169,000 PBMC scATAC-seq profiles and a best practices code repository for scATAC-seq data analysis, which are freely available to extend this benchmarking effort to future protocols.

SCENIC+: single-cell multiomic inference of enhancers and gene regulatory networks.

Joint profiling of chromatin accessibility and gene expression in individual cells provides an opportunity to decipher enhancer-driven gene regulatory networks (GRNs). Here we present a method for the inference of enhancer-driven GRNs, called SCENIC+. SCENIC+ predicts genomic enhancers along with candidate upstream transcription factors (TFs) and links these enhancers to candidate target genes. To improve both recall and precision of TF identification, we curated and clustered a motif collection with more than 30,000 motifs. We benchmarked SCENIC+ on diverse datasets from different species, including human peripheral blood mononuclear cells, ENCODE cell lines, melanoma cell states and Drosophila retinal development. Next, we exploit SCENIC+ predictions to study conserved TFs, enhancers and GRNs between human and mouse cell types in the cerebral cortex. Finally, we use SCENIC+ to study the dynamics of gene regulation along differentiation trajectories and the effect of TF perturbations on cell state. SCENIC+ is available at scenicplus.readthedocs.io .

Early alterations in the MCH system link aberrant neuronal activity and sleep disturbances in a mouse model of Alzheimer's disease.

Early Alzheimer's disease (AD) is associated with hippocampal hyperactivity and decreased sleep quality. Here we show that homeostatic mechanisms transiently counteract the increased excitatory drive to CA1 neurons in AppNL-G-F mice, but that this mechanism fails in older mice. Spatial transcriptomics analysis identifies Pmch as part of the adaptive response in AppNL-G-F mice. Pmch encodes melanin-concentrating hormone (MCH), which is produced in sleep-active lateral hypothalamic neurons that project to CA1 and modulate memory. We show that MCH downregulates synaptic transmission, modulates firing rate homeostasis in hippocampal neurons and reverses the increased excitatory drive to CA1 neurons in AppNL-G-F mice. AppNL-G-F mice spend less time in rapid eye movement (REM) sleep. AppNL-G-F mice and individuals with AD show progressive changes in morphology of CA1-projecting MCH axons. Our findings identify the MCH system as vulnerable in early AD and suggest that impaired MCH-system function contributes to aberrant excitatory drive and sleep defects, which can compromise hippocampus-dependent functions.

Molecular and cognitive signatures of ageing partially restored through synthetic delivery of IL2 to the brain.

Cognitive decline is a common pathological outcome during aging, with an ill-defined molecular and cellular basis. In recent years, the concept of inflammaging, defined as a low-grade inflammation increasing with age, has emerged. Infiltrating T cells accumulate in the brain with age and may contribute to the amplification of inflammatory cascades and disruptions to the neurogenic niche observed with age. Recently, a small resident population of regulatory T cells has been identified in the brain, and the capacity of IL2-mediated expansion of this population to counter neuroinflammatory disease has been demonstrated. Here, we test a brain-specific IL2 delivery system for the prevention of neurological decline in aging mice. We identify the molecular hallmarks of aging in the brain glial compartments and identify partial restoration of this signature through IL2 treatment. At a behavioral level, brain IL2 delivery prevented the age-induced defect in spatial learning, without improving the general decline in motor skill or arousal. These results identify immune modulation as a potential path to preserving cognitive function for healthy aging.

Glioblastoma-instructed microglia transition to heterogeneous phenotypic states with phagocytic and dendritic cell-like features in patient tumors and patient-derived orthotopic xenografts.

Background: A major contributing factor to glioblastoma (GBM) development and progression is its ability to evade the immune system by creating an immune-suppressive environment, where GBM-associated myeloid cells, including resident microglia and peripheral monocyte-derived macrophages, play critical pro-tumoral roles. However, it is unclear whether recruited myeloid cells are phenotypically and functionally identical in GBM patients and whether this heterogeneity is recapitulated in patient-derived orthotopic xenografts (PDOXs). A thorough understanding of the GBM ecosystem and its recapitulation in preclinical models is currently missing, leading to inaccurate results and failures of clinical trials. Methods: Here, we report systematic characterization of the tumor microenvironment (TME) in GBM PDOXs and patient tumors at the single-cell and spatial levels. We applied single-cell RNA-sequencing, spatial transcriptomics, multicolor flow cytometry, immunohistochemistry and functional studies to examine the heterogeneous TME instructed by GBM cells. GBM PDOXs representing different tumor phenotypes were compared to glioma mouse GL261 syngeneic model and patient tumors. Results: We show that GBM tumor cells reciprocally interact with host cells to create a GBM patient-specific TME in PDOXs. We detected the most prominent transcriptomic adaptations in myeloid cells, with brain-resident microglia representing the main population in the cellular tumor, while peripheral-derived myeloid cells infiltrated the brain at sites of blood-brain barrier disruption. More specifically, we show that GBM-educated microglia undergo transition to diverse phenotypic states across distinct GBM landscapes and tumor niches. GBM-educated microglia subsets display phagocytic and dendritic cell-like gene expression programs. Additionally, we found novel microglial states expressing cell cycle programs, astrocytic or endothelial markers. Lastly, we show that temozolomide treatment leads to transcriptomic plasticity and altered crosstalk between GBM tumor cells and adjacent TME components. Conclusions: Our data provide novel insights into the phenotypic adaptation of the heterogeneous TME instructed by GBM tumors. We show the key role of microglial phenotypic states in supporting GBM tumor growth and response to treatment. Our data place PDOXs as relevant models to assess the functionality of the TME and changes in the GBM ecosystem upon treatment.

Neuronal identity defines α-synuclein and tau toxicity.

Pathogenic α-synuclein and tau are critical drivers of neurodegeneration, and their mutations cause neuronal loss in patients. Whether the underlying preferential neuronal vulnerability is a cell-type-intrinsic property or a consequence of increased expression levels remains elusive. Here, we explore cell-type-specific α-synuclein and tau expression in human brain datasets and use deep phenotyping as well as brain-wide single-cell RNA sequencing of >200 live neuron types in fruit flies to determine which cellular environments react most to α-synuclein or tau toxicity. We detect phenotypic and transcriptomic evidence of differential neuronal vulnerability independent of α-synuclein or tau expression levels. Comparing vulnerable with resilient neurons in Drosophila enabled us to predict numerous human neuron subtypes with increased intrinsic susceptibility to pathogenic α-synuclein or tau. By uncovering synapse- and Ca2+ homeostasis-related genes as tau toxicity modifiers, our work paves the way to leverage neuronal identity to uncover modifiers of neurodegeneration-associated toxic proteins.

Mitochondria metabolism sets the species-specific tempo of neuronal development.

Neuronal development in the human cerebral cortex is considerably prolonged compared with that of other mammals. We explored whether mitochondria influence the species-specific timing of cortical neuron maturation. By comparing human and mouse cortical neuronal maturation at high temporal and cell resolution, we found a slower mitochondria development in human cortical neurons compared with that in the mouse, together with lower mitochondria metabolic activity, particularly that of oxidative phosphorylation. Stimulation of mitochondria metabolism in human neurons resulted in accelerated development in vitro and in vivo, leading to maturation of cells weeks ahead of time, whereas its inhibition in mouse neurons led to decreased rates of maturation. Mitochondria are thus important regulators of the pace of neuronal development underlying human-specific brain neoteny.

Cell type diversity in a developing octopus brain.

Octopuses are mollusks that have evolved intricate neural systems comparable with vertebrates in terms of cell number, complexity and size. The brain cell types that control their sophisticated behavioral repertoire are still unknown. Here, we profile the cell diversity of the paralarval Octopus vulgaris brain to build a cell type atlas that comprises mostly neural cells, but also multiple glial subtypes, endothelial cells and fibroblasts. We spatially map cell types to the vertical, subesophageal and optic lobes. Investigation of cell type conservation reveals a shared gene signature between glial cells of mouse, fly and octopus. Genes related to learning and memory are enriched in vertical lobe cells, which show molecular similarities with Kenyon cells in Drosophila. We construct a cell type taxonomy revealing transcriptionally related cell types, which tend to appear in the same brain region. Together, our data sheds light on cell type diversity and evolution in the octopus brain.

Astrocyte-targeted gene delivery of interleukin 2 specifically increases brain-resident regulatory T cell numbers and protects against pathological neuroinflammation.

The ability of immune-modulating biologics to prevent and reverse pathology has transformed recent clinical practice. Full utility in the neuroinflammation space, however, requires identification of both effective targets for local immune modulation and a delivery system capable of crossing the blood-brain barrier. The recent identification and characterization of a small population of regulatory T (Treg) cells resident in the brain presents one such potential therapeutic target. Here, we identified brain interleukin 2 (IL-2) levels as a limiting factor for brain-resident Treg cells. We developed a gene-delivery approach for astrocytes, with a small-molecule on-switch to allow temporal control, and enhanced production in reactive astrocytes to spatially direct delivery to inflammatory sites. Mice with brain-specific IL-2 delivery were protected in traumatic brain injury, stroke and multiple sclerosis models, without impacting the peripheral immune system. These results validate brain-specific IL-2 gene delivery as effective protection against neuroinflammation, and provide a versatile platform for delivery of diverse biologics to neuroinflammatory patients.

Comparative analysis of antibody- and lipid-based multiplexing methods for single-cell RNA-seq.

BACKGROUND: Multiplexing of samples in single-cell RNA-seq studies allows a significant reduction of the experimental costs, straightforward identification of doublets, increased cell throughput, and reduction of sample-specific batch effects. Recently published multiplexing techniques using oligo-conjugated antibodies or -lipids allow barcoding sample-specific cells, a process called "hashing." RESULTS: Here, we compare the hashing performance of TotalSeq-A and -C antibodies, custom synthesized lipids and MULTI-seq lipid hashes in four cell lines, both for single-cell RNA-seq and single-nucleus RNA-seq. We also compare TotalSeq-B antibodies with CellPlex reagents (10x Genomics) on human PBMCs and TotalSeq-B with different lipids on primary mouse tissues. Hashing efficiency was evaluated using the intrinsic genetic variation of the cell lines and mouse strains. Antibody hashing was further evaluated on clinical samples using PBMCs from healthy and SARS-CoV-2 infected patients, where we demonstrate a more affordable approach for large single-cell sequencing clinical studies, while simultaneously reducing batch effects. CONCLUSIONS: Benchmarking of different hashing strategies and computational pipelines indicates that correct demultiplexing can be achieved with both lipid- and antibody-hashed human cells and nuclei, with MULTISeqDemux as the preferred demultiplexing function and antibody-based hashing as the most efficient protocol on cells. On nuclei datasets, lipid hashing delivers the best results. Lipid hashing also outperforms antibodies on cells isolated from mouse brain. However, antibodies demonstrate better results on tissues like spleen or lung.