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.

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.

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.

TBK1 Suppresses RIPK1-Driven Apoptosis and Inflammation during Development and in Aging.

Aging is a major risk factor for both genetic and sporadic neurodegenerative disorders. However, it is unclear how aging interacts with genetic predispositions to promote neurodegeneration. Here, we investigate how partial loss of function of TBK1, a major genetic cause for amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD) comorbidity, leads to age-dependent neurodegeneration. We show that TBK1 is an endogenous inhibitor of RIPK1 and the embryonic lethality of Tbk1-/- mice is dependent on RIPK1 kinase activity. In aging human brains, another endogenous RIPK1 inhibitor, TAK1, exhibits a marked decrease in expression. We show that in Tbk1+/- mice, the reduced myeloid TAK1 expression promotes all the key hallmarks of ALS/FTD, including neuroinflammation, TDP-43 aggregation, axonal degeneration, neuronal loss, and behavior deficits, which are blocked upon inhibition of RIPK1. Thus, aging facilitates RIPK1 activation by reducing TAK1 expression, which cooperates with genetic risk factors to promote the onset of ALS/FTD.

Comprehensive Classification of Retinal Bipolar Neurons by Single-Cell Transcriptomics.

Patterns of gene expression can be used to characterize and classify neuronal types. It is challenging, however, to generate taxonomies that fulfill the essential criteria of being comprehensive, harmonizing with conventional classification schemes, and lacking superfluous subdivisions of genuine types. To address these challenges, we used massively parallel single-cell RNA profiling and optimized computational methods on a heterogeneous class of neurons, mouse retinal bipolar cells (BCs). From a population of ∼25,000 BCs, we derived a molecular classification that identified 15 types, including all types observed previously and two novel types, one of which has a non-canonical morphology and position. We validated the classification scheme and identified dozens of novel markers using methods that match molecular expression to cell morphology. This work provides a systematic methodology for achieving comprehensive molecular classification of neurons, identifies novel neuronal types, and uncovers transcriptional differences that distinguish types within a class.

Brain Chimeroids reveal individual susceptibility to neurotoxic triggers.

Interindividual genetic variation affects the susceptibility to and progression of many diseases1,2. However, efforts to study how individual human brains differ in normal development and disease phenotypes are limited by the paucity of faithful cellular human models, and the difficulty of scaling current systems to represent multiple people. Here we present human brain Chimeroids, a highly reproducible, multidonor human brain cortical organoid model generated by the co-development of cells from a panel of individual donors in a single organoid. By reaggregating cells from multiple single-donor organoids at the neural stem cell or neural progenitor cell stage, we generate Chimeroids in which each donor produces all cell lineages of the cerebral cortex, even when using pluripotent stem cell lines with notable growth biases. We used Chimeroids to investigate interindividual variation in the susceptibility to neurotoxic triggers that exhibit high clinical phenotypic variability: ethanol and the antiepileptic drug valproic acid. Individual donors varied in both the penetrance of the effect on target cell types, and the molecular phenotype within each affected cell type. Our results suggest that human genetic background may be an important mediator of neurotoxin susceptibility and introduce Chimeroids as a scalable system for high-throughput investigation of interindividual variation in processes of brain development and disease.

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.

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.

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.

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.

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.

Individual brain organoids reproducibly form cell diversity of the human cerebral cortex.

Experimental models of the human brain are needed for basic understanding of its development and disease1. Human brain organoids hold unprecedented promise for this purpose; however, they are plagued by high organoid-to-organoid variability2,3. This has raised doubts as to whether developmental processes of the human brain can occur outside the context of embryogenesis with a degree of reproducibility that is comparable to the endogenous tissue. Here we show that an organoid model of the dorsal forebrain can reliably generate a rich diversity of cell types appropriate for the human cerebral cortex. We performed single-cell RNA-sequencing analysis of 166,242 cells isolated from 21 individual organoids, finding that 95% of the organoids generate a virtually indistinguishable compendium of cell types, following similar developmental trajectories and with a degree of organoid-to-organoid variability comparable to that of individual endogenous brains. Furthermore, organoids derived from different stem cell lines show consistent reproducibility in the cell types produced. The data demonstrate that reproducible development of the complex cellular diversity of the central nervous system does not require the context of the embryo, and that establishment of terminal cell identity is a highly constrained process that can emerge from diverse stem cell origins and growth environments.

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.

Author Correction: Comprehensive comparative analysis of 5'-end RNA-sequencing methods.

The original version of this paper contained an incorrect primer sequence. In the Methods subsection "Rampage libraries," the text for modification 3 stated that the reverse primer used for library indexing was 5'-CAAGCAGAAGACGGCATACGAGATXXXXXXXXGTGACTGGAGT-3'. The correct sequence of the oligonucleotide used is 5'-CAAGCAGAAGACGGCATACGAGATXXXXXXXXGTGACTGGAGTTCAGACGTGTGCTCTTCCGATCT-3'. This error has been corrected in the PDF and HTML versions of the paper.

Comprehensive comparative analysis of 5'-end RNA-sequencing methods.

Specialized RNA-seq methods are required to identify the 5' ends of transcripts, which are critical for studies of gene regulation, but these methods have not been systematically benchmarked. We directly compared six such methods, including the performance of five methods on a single human cellular RNA sample and a new spike-in RNA assay that helps circumvent challenges resulting from uncertainties in annotation and RNA processing. We found that the 'cap analysis of gene expression' (CAGE) method performed best for mRNA and that most of its unannotated peaks were supported by evidence from other genomic methods. We applied CAGE to eight brain-related samples and determined sample-specific transcription start site (TSS) usage, as well as a transcriptome-wide shift in TSS usage between fetal and adult brain.

RIPK1 mediates a disease-associated microglial response in Alzheimer's disease.

Dysfunction of microglia is known to play an important role in Alzheimer's disease (AD). Here, we investigated the role of RIPK1 in microglia mediating the pathogenesis of AD. RIPK1 is highly expressed by microglial cells in human AD brains. Using the amyloid precursor protein (APP)/presenilin 1 (PS1) transgenic mouse model, we found that inhibition of RIPK1, using both pharmacological and genetic means, reduced amyloid burden, the levels of inflammatory cytokines, and memory deficits. Furthermore, inhibition of RIPK1 promoted microglial degradation of Aß in vitro. We characterized the transcriptional profiles of adult microglia from APP/PS1 mice and identified a role for RIPK1 in regulating the microglial expression of CH25H and Cst7, a marker for disease-associated microglia (DAM), which encodes an endosomal/lysosomal cathepsin inhibitor named Cystatin F. We present evidence that RIPK1-mediated induction of Cst7 leads to an impairment in the lysosomal pathway. These data suggest that RIPK1 may mediate a critical checkpoint in the transition to the DAM state. Together, our study highlights a non-cell death mechanism by which the activation of RIPK1 mediates the induction of a DAM phenotype, including an inflammatory response and a reduction in phagocytic activity, and connects RIPK1-mediated transcription in microglia to the etiology of AD. Our results support that RIPK1 is an important therapeutic target for the treatment of AD.

Single-cell transcriptomics reveals bimodality in expression and splicing in immune cells.

Recent molecular studies have shown that, even when derived from a seemingly homogenous population, individual cells can exhibit substantial differences in gene expression, protein levels and phenotypic output, with important functional consequences. Existing studies of cellular heterogeneity, however, have typically measured only a few pre-selected RNAs or proteins simultaneously, because genomic profiling methods could not be applied to single cells until very recently. Here we use single-cell RNA sequencing to investigate heterogeneity in the response of mouse bone-marrow-derived dendritic cells (BMDCs) to lipopolysaccharide. We find extensive, and previously unobserved, bimodal variation in messenger RNA abundance and splicing patterns, which we validate by RNA-fluorescence in situ hybridization for select transcripts. In particular, hundreds of key immune genes are bimodally expressed across cells, surprisingly even for genes that are very highly expressed at the population average. Moreover, splicing patterns demonstrate previously unobserved levels of heterogeneity between cells. Some of the observed bimodality can be attributed to closely related, yet distinct, known maturity states of BMDCs; other portions reflect differences in the usage of key regulatory circuits. For example, we identify a module of 137 highly variable, yet co-regulated, antiviral response genes. Using cells from knockout mice, we show that variability in this module may be propagated through an interferon feedback circuit, involving the transcriptional regulators Stat2 and Irf7. Our study demonstrates the power and promise of single-cell genomics in uncovering functional diversity between cells and in deciphering cell states and circuits.

The African coelacanth genome provides insights into tetrapod evolution.

The discovery of a living coelacanth specimen in 1938 was remarkable, as this lineage of lobe-finned fish was thought to have become extinct 70 million years ago. The modern coelacanth looks remarkably similar to many of its ancient relatives, and its evolutionary proximity to our own fish ancestors provides a glimpse of the fish that first walked on land. Here we report the genome sequence of the African coelacanth, Latimeria chalumnae. Through a phylogenomic analysis, we conclude that the lungfish, and not the coelacanth, is the closest living relative of tetrapods. Coelacanth protein-coding genes are significantly more slowly evolving than those of tetrapods, unlike other genomic features. Analyses of changes in genes and regulatory elements during the vertebrate adaptation to land highlight genes involved in immunity, nitrogen excretion and the development of fins, tail, ear, eye, brain and olfaction. Functional assays of enhancers involved in the fin-to-limb transition and in the emergence of extra-embryonic tissues show the importance of the coelacanth genome as a blueprint for understanding tetrapod evolution.

Use of the MS2 aptamer and coat protein for RNA localization in yeast: A response to "MS2 coat proteins bound to yeast mRNAs block 5' to 3' degradation and trap mRNA decay products: implications for the localization of mRNAs by MS2-MCP system".

The MS2 system has been extensively used to visualize single mRNA molecules in live cells and follow their localization and behavior. In their Letter to the Editor recently published, Garcia and Parker suggest that use of the MS2 system may yield erroneous mRNA localization results due to the accumulation of 3' decay products. Here we cite published works and provide new data which demonstrate that this is not a phenomenon general to endogenously expressed MS2-tagged transcripts, and that some of the results obtained in their study could have arisen from artifacts of gene expression.