Fetal brain response to maternal inflammation requires microglia.

In utero infection and maternal inflammation can adversely impact fetal brain development. Maternal systemic illness, even in the absence of direct fetal brain infection, is associated with an increased risk of neuropsychiatric disorders in affected offspring. The cell types mediating the fetal brain response to maternal inflammation are largely unknown, hindering the development of novel treatment strategies. Here, we show that microglia, the resident phagocytes of the brain, highly express receptors for relevant pathogens and cytokines throughout embryonic development. Using a rodent maternal immune activation (MIA) model in which polyinosinic:polycytidylic acid is injected into pregnant mice, we demonstrate long-lasting transcriptional changes in fetal microglia that persist into postnatal life. We find that MIA induces widespread gene expression changes in neuronal and non-neuronal cells; importantly, these responses are abolished by selective genetic deletion of microglia, indicating that microglia are required for the transcriptional response of other cortical cell types to MIA. These findings demonstrate that microglia play a crucial durable role in the fetal response to maternal inflammation, and should be explored as potential therapeutic cell targets.

Massively parallel in vivo Perturb-seq reveals cell-type-specific transcriptional networks in cortical development.

Leveraging AAVs' versatile tropism and labeling capacity, we expanded the scale of in vivo CRISPR screening with single-cell transcriptomic phenotyping across embryonic to adult brains and peripheral nervous systems. Through extensive tests of 86 vectors across AAV serotypes combined with a transposon system, we substantially amplified labeling efficacy and accelerated in vivo gene delivery from weeks to days. Our proof-of-principle in utero screen identified the pleiotropic effects of Foxg1, highlighting its tight regulation of distinct networks essential for cell fate specification of Layer 6 corticothalamic neurons. Notably, our platform can label >6% of cerebral cells, surpassing the current state-of-the-art efficacy at <0.1% by lentivirus, to achieve analysis of over 30,000 cells in one experiment and enable massively parallel in vivo Perturb-seq. Compatible with various phenotypic measurements (single-cell or spatial multi-omics), it presents a flexible approach to interrogate gene function across cell types in vivo, translating gene variants to their causal function.

Lupus autoantibodies initiate neuroinflammation sustained by continuous HMGB1:RAGE signaling and reversed by increased LAIR-1 expression.

Cognitive impairment is a frequent manifestation of neuropsychiatric systemic lupus erythematosus, present in up to 80% of patients and leading to a diminished quality of life. In the present study, we used a model of lupus-like cognitive impairment that is initiated when antibodies that crossreact with excitatory neuronal receptors penetrate the hippocampus, causing immediate, self-limited, excitotoxic death of hippocampal neurons, which is then followed by a significant loss of dendritic complexity in surviving neurons. This injury creates a maladaptive equilibrium that is sustained in mice for at least 1 year. We identified a feedforward loop of microglial activation and microglia-dependent synapse elimination dependent on neuronal secretion of high mobility group box 1 protein (HMGB1) which binds the receptor for advanced glycation end products (RAGE) and leads to microglial secretion of C1q, upregulation of interleukin-10 with consequent downregulation of leukocyte-associated immunoglobulin-like receptor 1 (LAIR-1), an inhibitory receptor for C1q. Treatment with a centrally acting angiotensin-converting enzyme inhibitor or with an angiotensin-receptor blocker restored a healthy equilibrium, microglial quiescence and intact spatial memory.

Massively parallel in vivo Perturb-seq reveals cell type-specific transcriptional networks in cortical development.

Systematic analysis of gene function across diverse cell types in vivo is hindered by two challenges: obtaining sufficient cells from live tissues and accurately identifying each cell's perturbation in high-throughput single-cell assays. Leveraging AAV's versatile cell type tropism and high labeling capacity, we expanded the resolution and scale of in vivo CRISPR screens: allowing phenotypic analysis at single-cell resolution across a multitude of cell types in the embryonic brain, adult brain, and peripheral nervous system. We undertook extensive tests of 86 AAV serotypes, combined with a transposon system, to substantially amplify labeling and accelerate in vivo gene delivery from weeks to days. Using this platform, we performed an in utero genetic screen as proof-of-principle and identified pleiotropic regulatory networks of Foxg1 in cortical development, including Layer 6 corticothalamic neurons where it tightly controls distinct networks essential for cell fate specification. Notably, our platform can label >6% of cerebral cells, surpassing the current state-of-the-art efficacy at <0.1% (mediated by lentivirus), and achieve analysis of over 30,000 cells in one experiment, thus enabling massively parallel in vivo Perturb-seq. Compatible with various perturbation techniques (CRISPRa/i) and phenotypic measurements (single-cell or spatial multi-omics), our platform presents a flexible, modular approach to interrogate gene function across diverse cell types in vivo, connecting gene variants to their causal functions.

Brain-region-specific changes in neurons and glia and dysregulation of dopamine signaling in Grin2a mutant mice.

A genetically valid animal model could transform our understanding of schizophrenia (SCZ) disease mechanisms. Rare heterozygous loss-of-function (LoF) mutations in GRIN2A, encoding a subunit of the NMDA receptor, greatly increase the risk of SCZ. By transcriptomic, proteomic, and behavioral analyses, we report that heterozygous Grin2a mutant mice show (1) large-scale gene expression changes across multiple brain regions and in neuronal (excitatory and inhibitory) and non-neuronal cells (astrocytes and oligodendrocytes), (2) evidence of hypoactivity in the prefrontal cortex (PFC) and hyperactivity in the hippocampus and striatum, (3) an elevated dopamine signaling in the striatum and hypersensitivity to amphetamine-induced hyperlocomotion (AIH), (4) altered cholesterol biosynthesis in astrocytes, (5) a reduction in glutamatergic receptor signaling proteins in the synapse, and (6) an aberrant locomotor pattern opposite of that induced by antipsychotic drugs. These findings reveal potential pathophysiologic mechanisms, provide support for both the "hypo-glutamate" and "hyper-dopamine" hypotheses of SCZ, and underscore the utility of Grin2a-deficient mice as a genetic model of SCZ.

Lupus autoantibodies initiate a maladaptive equilibrium sustained by HMGB1:RAGE signaling and reversed by LAIR-1:C1q signaling.

Cognitive impairment is a frequent manifestation of neuropsychiatric systemic lupus erythematosus (NPSLE), present in up to 80% of patients and leading to a diminished quality of life. We have developed a model of lupus-like cognitive impairment which is initiated when anti-DNA, anti-N-methyl D-aspartate receptor (NMDAR) cross- reactive antibodies, which are present in 30% of SLE patients, penetrate the hippocampus1. This leads to immediate, self-limited excitotoxic death of CA1 pyramidal neurons followed by a significant loss of dendritic arborization in the remaining CA1 neurons and impaired spatial memory. Both microglia and C1q are required for dendritic loss1. Here we show that this pattern of hippocampal injury creates a maladaptive equilibrium that is sustained for at least one year. It requires HMGB1 secretion by neurons to bind RAGE, a receptor for HMGB1 expressed on microglia, and leads to decreased expression of microglial LAIR-1, an inhibitory receptor for C1q. The angiotensin converting enzyme (ACE) inhibitor captopril, which can restore a healthy equilibrium, microglial quiescence, and intact spatial memory, leads to upregulation of LAIR-1. This paradigm highlights HMGB1:RAGE and C1q:LAIR-1 interactions as pivotal pathways in the microglial-neuronal interplay that defines a physiologic versus a maladaptive equilibrium.

The effect of background noise and its removal on the analysis of single-cell expression data.

BACKGROUND: In droplet-based single-cell and single-nucleus RNA-seq experiments, not all reads associated with one cell barcode originate from the encapsulated cell. Such background noise is attributed to spillage from cell-free ambient RNA or barcode swapping events. RESULTS: Here, we characterize this background noise exemplified by three scRNA-seq and two snRNA-seq replicates of mouse kidneys. For each experiment, cells from two mouse subspecies are pooled, allowing to identify cross-genotype contaminating molecules and thus profile background noise. Background noise is highly variable across replicates and cells, making up on average 3-35% of the total counts (UMIs) per cell and we find that noise levels are directly proportional to the specificity and detectability of marker genes. In search of the source of background noise, we find multiple lines of evidence that the majority of background molecules originates from ambient RNA. Finally, we use our genotype-based estimates to evaluate the performance of three methods (CellBender, DecontX, SoupX) that are designed to quantify and remove background noise. We find that CellBender provides the most precise estimates of background noise levels and also yields the highest improvement for marker gene detection. By contrast, clustering and classification of cells are fairly robust towards background noise and only small improvements can be achieved by background removal that may come at the cost of distortions in fine structure. CONCLUSIONS: Our findings help to better understand the extent, sources and impact of background noise in single-cell experiments and provide guidance on how to deal with it.

Cell-type specific defects in PTEN-mutant cortical organoids converge on abnormal circuit activity.

De novo heterozygous loss-of-function mutations in phosphatase and tensin homolog (PTEN) are strongly associated with autism spectrum disorders; however, it is unclear how heterozygous mutations in this gene affect different cell types during human brain development and how these effects vary across individuals. Here, we used human cortical organoids from different donors to identify cell-type specific developmental events that are affected by heterozygous mutations in PTEN. We profiled individual organoids by single-cell RNA-seq, proteomics and spatial transcriptomics and revealed abnormalities in developmental timing in human outer radial glia progenitors and deep-layer cortical projection neurons, which varied with the donor genetic background. Calcium imaging in intact organoids showed that both accelerated and delayed neuronal development phenotypes resulted in similar abnormal activity of local circuits, irrespective of genetic background. The work reveals donor-dependent, cell-type specific developmental phenotypes of PTEN heterozygosity that later converge on disrupted neuronal activity.

Heterochronic parabiosis reprograms the mouse brain transcriptome by shifting aging signatures in multiple cell types.

Aging is a complex process involving transcriptomic changes associated with deterioration across multiple tissues and organs, including the brain. Recent studies using heterochronic parabiosis have shown that various aspects of aging-associated decline are modifiable or even reversible. To better understand how this occurs, we performed single-cell transcriptomic profiling of young and old mouse brains after parabiosis. For each cell type, we cataloged alterations in gene expression, molecular pathways, transcriptional networks, ligand-receptor interactions and senescence status. Our analyses identified gene signatures, demonstrating that heterochronic parabiosis regulates several hallmarks of aging in a cell-type-specific manner. Brain endothelial cells were found to be especially malleable to this intervention, exhibiting dynamic transcriptional changes that affect vascular structure and function. These findings suggest new strategies for slowing deterioration and driving regeneration in the aging brain through approaches that do not rely on disease-specific mechanisms or actions of individual circulating factors.

Mostly natural sequencing-by-synthesis for scRNA-seq using Ultima sequencing.

Here we introduce a mostly natural sequencing-by-synthesis (mnSBS) method for single-cell RNA sequencing (scRNA-seq), adapted to the Ultima genomics platform, and systematically benchmark it against current scRNA-seq technology. mnSBS uses mostly natural, unmodified nucleotides and only a low fraction of fluorescently labeled nucleotides, which allows for high polymerase processivity and lower costs. We demonstrate successful application in four scRNA-seq case studies of different technical and biological types, including 5' and 3' scRNA-seq, human peripheral blood mononuclear cells from a single individual and in multiplex, as well as Perturb-Seq. Benchmarking shows that results from mnSBS-based scRNA-seq are very similar to those using Illumina sequencing, with minor differences in results related to the position of reads relative to annotated gene boundaries, owing to single-end reads of Ultima being closer to gene ends than reads from Illumina. The method is thus compatible with state-of-the-art scRNA-seq libraries independent of the sequencing technology. We expect mnSBS to be of particular utility for cost-effective large-scale scRNA-seq projects.

Proper acquisition of cell class identity in organoids allows definition of fate specification programs of the human cerebral cortex.

Realizing the full utility of brain organoids to study human development requires understanding whether organoids precisely replicate endogenous cellular and molecular events, particularly since acquisition of cell identity in organoids can be impaired by abnormal metabolic states. We present a comprehensive single-cell transcriptomic, epigenetic, and spatial atlas of human cortical organoid development, comprising over 610,000 cells, from generation of neural progenitors through production of differentiated neuronal and glial subtypes. We show that processes of cellular diversification correlate closely to endogenous ones, irrespective of metabolic state, empowering the use of this atlas to study human fate specification. We define longitudinal molecular trajectories of cortical cell types during organoid development, identify genes with predicted human-specific roles in lineage establishment, and uncover early transcriptional diversity of human callosal neurons. The findings validate this comprehensive atlas of human corticogenesis in vitro as a resource to prime investigation into the mechanisms of human cortical development.

Pyramidal neuron subtype diversity governs microglia states in the neocortex.

Microglia are specialized macrophages in the brain parenchyma that exist in multiple transcriptional states and reside within a wide range of neuronal environments1-4. However, how and where these states are generated remains poorly understood. Here, using the mouse somatosensory cortex, we demonstrate that microglia density and molecular state acquisition are determined by the local composition of pyramidal neuron classes. Using single-cell and spatial transcriptomic profiling, we unveil the molecular signatures and spatial distributions of diverse microglia populations and show that certain states are enriched in specific cortical layers, whereas others are broadly distributed throughout the cortex. Notably, conversion of deep-layer pyramidal neurons to an alternate class identity reconfigures the distribution of local, layer-enriched homeostatic microglia to match the new neuronal niche. Leveraging the transcriptional diversity of pyramidal neurons in the neocortex, we construct a ligand-receptor atlas describing interactions between individual pyramidal neuron subtypes and microglia states, revealing rules of neuron-microglia communication. Our findings uncover a fundamental role for neuronal diversity in instructing the acquisition of microglia states as a potential mechanism for fine-tuning neuroimmune interactions within the cortical local circuitry.

Temporally divergent regulatory mechanisms govern neuronal diversification and maturation in the mouse and marmoset neocortex.

Mammalian neocortical neurons span one of the most diverse cell type spectra of any tissue. Cortical neurons are born during embryonic development, and their maturation extends into postnatal life. The regulatory strategies underlying progressive neuronal development and maturation remain unclear. Here we present an integrated single-cell epigenomic and transcriptional analysis of individual mouse and marmoset cortical neuron classes, spanning both early postmitotic stages of identity acquisition and later stages of neuronal plasticity and circuit integration. We found that, in both species, the regulatory strategies controlling early and late stages of pan-neuronal development diverge. Early postmitotic neurons use more widely shared and evolutionarily conserved molecular regulatory programs. In contrast, programs active during later neuronal maturation are more brain- and neuron-specific and more evolutionarily divergent. Our work uncovers a temporal shift in regulatory choices during neuronal diversification and maturation in both mice and marmosets, which likely reflects unique evolutionary constraints on distinct events of neuronal development in the neocortex.

The evolution, evolvability and engineering of gene regulatory DNA.

Mutations in non-coding regulatory DNA sequences can alter gene expression, organismal phenotype and fitness1-3. Constructing complete fitness landscapes, in which DNA sequences are mapped to fitness, is a long-standing goal in biology, but has remained elusive because it is challenging to generalize reliably to vast sequence spaces4-6. Here we build sequence-to-expression models that capture fitness landscapes and use them to decipher principles of regulatory evolution. Using millions of randomly sampled promoter DNA sequences and their measured expression levels in the yeast Saccharomyces cerevisiae, we learn deep neural network models that generalize with excellent prediction performance, and enable sequence design for expression engineering. Using our models, we study expression divergence under genetic drift and strong-selection weak-mutation regimes to find that regulatory evolution is rapid and subject to diminishing returns epistasis; that conflicting expression objectives in different environments constrain expression adaptation; and that stabilizing selection on gene expression leads to the moderation of regulatory complexity. We present an approach for using such models to detect signatures of selection on expression from natural variation in regulatory sequences and use it to discover an instance of convergent regulatory evolution. We assess mutational robustness, finding that regulatory mutation effect sizes follow a power law, characterize regulatory evolvability, visualize promoter fitness landscapes, discover evolvability archetypes and illustrate the mutational robustness of natural regulatory sequence populations. Our work provides a general framework for designing regulatory sequences and addressing fundamental questions in regulatory evolution.

Autism genes converge on asynchronous development of shared neuron classes.

Genetic risk for autism spectrum disorder (ASD) is associated with hundreds of genes spanning a wide range of biological functions1-6. The alterations in the human brain resulting from mutations in these genes remain unclear. Furthermore, their phenotypic manifestation varies across individuals7,8. Here we used organoid models of the human cerebral cortex to identify cell-type-specific developmental abnormalities that result from haploinsufficiency in three ASD risk genes-SUV420H1 (also known as KMT5B), ARID1B and CHD8-in multiple cell lines from different donors, using single-cell RNA-sequencing (scRNA-seq) analysis of more than 745,000 cells and proteomic analysis of individual organoids, to identify phenotypic convergence. Each of the three mutations confers asynchronous development of two main cortical neuronal lineages-γ-aminobutyric-acid-releasing (GABAergic) neurons and deep-layer excitatory projection neurons-but acts through largely distinct molecular pathways. Although these phenotypes are consistent across cell lines, their expressivity is influenced by the individual genomic context, in a manner that is dependent on both the risk gene and the developmental defect. Calcium imaging in intact organoids shows that these early-stage developmental changes are followed by abnormal circuit activity. This research uncovers cell-type-specific neurodevelopmental abnormalities that are shared across ASD risk genes and are finely modulated by human genomic context, finding convergence in the neurobiological basis of how different risk genes contribute to ASD pathology.

A RIPK1-regulated inflammatory microglial state in amyotrophic lateral sclerosis.

Microglial-derived inflammation has been linked to a broad range of neurodegenerative and neuropsychiatric conditions, including amyotrophic lateral sclerosis (ALS). Using single-cell RNA sequencing, a class of Disease-Associated Microglia (DAMs) have been characterized in neurodegeneration. However, the DAM phenotype alone is insufficient to explain the functional complexity of microglia, particularly with regard to regulating inflammation that is a hallmark of many neurodegenerative diseases. Here, we identify a subclass of microglia in mouse models of ALS which we term RIPK1-Regulated Inflammatory Microglia (RRIMs). RRIMs show significant up-regulation of classical proinflammatory pathways, including increased levels of Tnf and Il1b RNA and protein. We find that RRIMs are highly regulated by TNFα signaling and that the prevalence of these microglia can be suppressed by inhibiting receptor-interacting protein kinase 1 (RIPK1) activity downstream of the TNF receptor 1. These findings help to elucidate a mechanism by which RIPK1 kinase inhibition has been shown to provide therapeutic benefit in mouse models of ALS and may provide an additional biomarker for analysis in ongoing phase 2 clinical trials of RIPK1 inhibitors in ALS.

Pluripotent stem cell-derived models of neurological diseases reveal early transcriptional heterogeneity.

BACKGROUND: Many neurodegenerative diseases develop only later in life, when cells in the nervous system lose their structure or function. In many forms of neurodegenerative diseases, this late-onset phenomenon remains largely unexplained. RESULTS: Analyzing single-cell RNA sequencing from Alzheimer's disease (AD) and Huntington's disease (HD) patients, we find increased transcriptional heterogeneity in disease-state neurons. We hypothesize that transcriptional heterogeneity precedes neurodegenerative disease pathologies. To test this idea experimentally, we use juvenile forms (72Q; 180Q) of HD iPSCs, differentiate them into committed neuronal progenitors, and obtain single-cell expression profiles. We show a global increase in gene expression variability in HD. Autophagy genes become more stable, while energy and actin-related genes become more variable in the mutant cells. Knocking down several differentially variable genes results in increased aggregate formation, a pathology associated with HD. We further validate the increased transcriptional heterogeneity in CHD8+/- cells, a model for autism spectrum disorder. CONCLUSIONS: Overall, our results suggest that although neurodegenerative diseases develop over time, transcriptional regulation imbalance is present already at very early developmental stages. Therefore, an intervention aimed at this early phenotype may be of high diagnostic value.

Reduction of mNAT1/hNAT2 Contributes to Cerebral Endothelial Necroptosis and Aß Accumulation in Alzheimer's Disease.

The contribution and mechanism of cerebrovascular pathology in Alzheimer's disease (AD) pathogenesis are still unclear. Here, we show that venular and capillary cerebral endothelial cells (ECs) are selectively vulnerable to necroptosis in AD. We identify reduced cerebromicrovascular expression of murine N-acetyltransferase 1 (mNat1) in two AD mouse models and hNat2, the human ortholog of mNat1 and a genetic risk factor for type-2 diabetes and insulin resistance, in human AD. mNat1 deficiency in Nat1-/- mice and two AD mouse models promotes blood-brain barrier (BBB) damage and endothelial necroptosis. Decreased mNat1 expression induces lysosomal degradation of A20, an important regulator of necroptosis, and LRP1ß, a key component of LRP1 complex that exports Aß in cerebral ECs. Selective restoration of cerebral EC expression of mNAT1 delivered by adeno-associated virus (AAV) rescues cerebromicrovascular levels of A20 and LRP1ß, inhibits endothelial necroptosis and activation, ameliorates mitochondrial fragmentation, reduces Aß deposits, and improves cognitive function in the AD mouse model.

In vivo Perturb-Seq reveals neuronal and glial abnormalities associated with autism risk genes.

The number of disease risk genes and loci identified through human genetic studies far outstrips the capacity to systematically study their functions. We applied a scalable genetic screening approach, in vivo Perturb-Seq, to functionally evaluate 35 autism spectrum disorder/neurodevelopmental delay (ASD/ND) de novo loss-of-function risk genes. Using CRISPR-Cas9, we introduced frameshift mutations in these risk genes in pools, within the developing mouse brain in utero, followed by single-cell RNA-sequencing of perturbed cells in the postnatal brain. We identified cell type-specific and evolutionarily conserved gene modules from both neuronal and glial cell classes. Recurrent gene modules and cell types are affected across this cohort of perturbations, representing key cellular effects across sets of ASD/ND risk genes. In vivo Perturb-Seq allows us to investigate how diverse mutations affect cell types and states in the developing organism.

Distinct subnetworks of the thalamic reticular nucleus.

The thalamic reticular nucleus (TRN), the major source of thalamic inhibition, regulates thalamocortical interactions that are critical for sensory processing, attention and cognition1-5. TRN dysfunction has been linked to sensory abnormality, attention deficit and sleep disturbance across multiple neurodevelopmental disorders6-9. However, little is known about the organizational principles that underlie its divergent functions. Here we performed an integrative study linking single-cell molecular and electrophysiological features of the mouse TRN to connectivity and systems-level function. We found that cellular heterogeneity in the TRN is characterized by a transcriptomic gradient of two negatively correlated gene-expression profiles, each containing hundreds of genes. Neurons in the extremes of this transcriptomic gradient express mutually exclusive markers, exhibit core or shell-like anatomical structure and have distinct electrophysiological properties. The two TRN subpopulations make differential connections with the functionally distinct first-order and higher-order thalamic nuclei to form molecularly defined TRN-thalamus subnetworks. Selective perturbation of the two subnetworks in vivo revealed their differential role in regulating sleep. In sum, our study provides a comprehensive atlas of TRN neurons at single-cell resolution and links molecularly defined subnetworks to the functional organization of thalamocortical circuits.