Molecularly defined and spatially resolved cell atlas of the whole mouse brain.

In mammalian brains, millions to billions of cells form complex interaction networks to enable a wide range of functions. The enormous diversity and intricate organization of cells have impeded our understanding of the molecular and cellular basis of brain function. Recent advances in spatially resolved single-cell transcriptomics have enabled systematic mapping of the spatial organization of molecularly defined cell types in complex tissues1-3, including several brain regions (for example, refs. 1-11). However, a comprehensive cell atlas of the whole brain is still missing. Here we imaged a panel of more than 1,100 genes in approximately 10 million cells across the entire adult mouse brains using multiplexed error-robust fluorescence in situ hybridization12 and performed spatially resolved, single-cell expression profiling at the whole-transcriptome scale by integrating multiplexed error-robust fluorescence in situ hybridization and single-cell RNA sequencing data. Using this approach, we generated a comprehensive cell atlas of more than 5,000 transcriptionally distinct cell clusters, belonging to more than 300 major cell types, in the whole mouse brain with high molecular and spatial resolution. Registration of this atlas to the mouse brain common coordinate framework allowed systematic quantifications of the cell-type composition and organization in individual brain regions. We further identified spatial modules characterized by distinct cell-type compositions and spatial gradients featuring gradual changes of cells. Finally, this high-resolution spatial map of cells, each with a transcriptome-wide expression profile, allowed us to infer cell-type-specific interactions between hundreds of cell-type pairs and predict molecular (ligand-receptor) basis and functional implications of these cell-cell interactions. These results provide rich insights into the molecular and cellular architecture of the brain and a foundation for functional investigations of neural circuits and their dysfunction in health and disease.

A molecularly defined and spatially resolved cell atlas of the whole mouse brain.

In mammalian brains, tens of millions to billions of cells form complex interaction networks to enable a wide range of functions. The enormous diversity and intricate organization of cells in the brain have so far hindered our understanding of the molecular and cellular basis of its functions. Recent advances in spatially resolved single-cell transcriptomics have allowed systematic mapping of the spatial organization of molecularly defined cell types in complex tissues1-3. However, these approaches have only been applied to a few brain regions1-11 and a comprehensive cell atlas of the whole brain is still missing. Here, we imaged a panel of >1,100 genes in ~8 million cells across the entire adult mouse brain using multiplexed error-robust fluorescence in situ hybridization (MERFISH)12 and performed spatially resolved, single-cell expression profiling at the whole-transcriptome scale by integrating MERFISH and single-cell RNA-sequencing (scRNA-seq) data. Using this approach, we generated a comprehensive cell atlas of >5,000 transcriptionally distinct cell clusters, belonging to ~300 major cell types, in the whole mouse brain with high molecular and spatial resolution. Registration of the MERFISH images to the common coordinate framework (CCF) of the mouse brain further allowed systematic quantifications of the cell composition and organization in individual brain regions defined in the CCF. We further identified spatial modules characterized by distinct cell-type compositions and spatial gradients featuring gradual changes in the gene-expression profiles of cells. Finally, this high-resolution spatial map of cells, with a transcriptome-wide expression profile associated with each cell, allowed us to infer cell-type-specific interactions between several hundred pairs of molecularly defined cell types and predict potential molecular (ligand-receptor) basis and functional implications of these cell-cell interactions. These results provide rich insights into the molecular and cellular architecture of the brain and a valuable resource for future functional investigations of neural circuits and their dysfunction in diseases.

A preoptic neuronal population controls fever and appetite during sickness.

During infection, animals exhibit adaptive changes in physiology and behaviour aimed at increasing survival. Although many causes of infection exist, they trigger similar stereotyped symptoms such as fever, warmth-seeking, loss of appetite and fatigue1,2. Yet exactly how the nervous system alters body temperature and triggers sickness behaviours to coordinate responses to infection remains unknown. Here we identify a previously uncharacterized population of neurons in the ventral medial preoptic area (VMPO) of the hypothalamus that are activated after sickness induced by lipopolysaccharide (LPS) or polyinosinic:polycytidylic acid. These neurons are crucial for generating a fever response and other sickness symptoms such as warmth-seeking and loss of appetite. Single-nucleus RNA-sequencing and multiplexed error-robust fluorescence in situ hybridization uncovered the identity and distribution of LPS-activated VMPO (VMPOLPS) neurons and non-neuronal cells. Gene expression and electrophysiological measurements implicate a paracrine mechanism in which the release of immune signals by non-neuronal cells during infection activates nearby VMPOLPS neurons. Finally, we show that VMPOLPS neurons exert a broad influence on the activity of brain areas associated with behavioural and homeostatic functions and are synaptically and functionally connected to circuit nodes controlling body temperature and appetite. Together, these results uncover VMPOLPS neurons as a control hub that integrates immune signals to orchestrate multiple sickness symptoms in response to infection.

Spatially resolved cell atlas of the mouse primary motor cortex by MERFISH.

A mammalian brain is composed of numerous cell types organized in an intricate manner to form functional neural circuits. Single-cell RNA sequencing allows systematic identification of cell types based on their gene expression profiles and has revealed many distinct cell populations in the brain1,2. Single-cell epigenomic profiling3,4 further provides information on gene-regulatory signatures of different cell types. Understanding how different cell types contribute to brain function, however, requires knowledge of their spatial organization and connectivity, which is not preserved in sequencing-based methods that involve cell dissociation. Here we used a single-cell transcriptome-imaging method, multiplexed error-robust fluorescence in situ hybridization (MERFISH)5, to generate a molecularly defined and spatially resolved cell atlas of the mouse primary motor cortex. We profiled approximately 300,000 cells in the mouse primary motor cortex and its adjacent areas, identified 95 neuronal and non-neuronal cell clusters, and revealed a complex spatial map in which not only excitatory but also most inhibitory neuronal clusters adopted laminar organizations. Intratelencephalic neurons formed a largely continuous gradient along the cortical depth axis, in which the gene expression of individual cells correlated with their cortical depths. Furthermore, we integrated MERFISH with retrograde labelling to probe projection targets of neurons of the mouse primary motor cortex and found that their cortical projections formed a complex network in which individual neuronal clusters project to multiple target regions and individual target regions receive inputs from multiple neuronal clusters.

Interactions between cancer cells and immune cells drive transitions to mesenchymal-like states in glioblastoma.

The mesenchymal subtype of glioblastoma is thought to be determined by both cancer cell-intrinsic alterations and extrinsic cellular interactions, but remains poorly understood. Here, we dissect glioblastoma-to-microenvironment interactions by single-cell RNA sequencing analysis of human tumors and model systems, combined with functional experiments. We demonstrate that macrophages induce a transition of glioblastoma cells into mesenchymal-like (MES-like) states. This effect is mediated, both in vitro and in vivo, by macrophage-derived oncostatin M (OSM) that interacts with its receptors (OSMR or LIFR) in complex with GP130 on glioblastoma cells and activates STAT3. We show that MES-like glioblastoma states are also associated with increased expression of a mesenchymal program in macrophages and with increased cytotoxicity of T cells, highlighting extensive alterations of the immune microenvironment with potential therapeutic implications.

MicroRNAs Cause Accelerated Decay of Short-Tailed Target mRNAs.

MicroRNAs (miRNAs) specify the recruitment of deadenylases to mRNA targets. Despite this recruitment, we find that miRNAs have almost no effect on steady-state poly(A)-tail lengths of their targets in mouse fibroblasts, which motivates the acquisition of pre-steady-state measurements of the effects of miRNAs on tail lengths, mRNA levels, and translational efficiencies. Effects on translational efficiency are minimal compared to effects on mRNA levels, even for newly transcribed target mRNAs. Effects on target mRNA levels accumulate as the mRNA population approaches steady state, whereas effects on tail lengths peak for recently transcribed target mRNAs and then subside. Computational modeling of this phenomenon reveals that miRNAs cause not only accelerated deadenylation of their targets but also accelerated decay of short-tailed target molecules. This unanticipated effect of miRNAs largely prevents short-tailed target mRNAs from accumulating despite accelerated target deadenylation. The net result is a nearly imperceptible change to the steady-state tail-length distribution of targeted mRNAs.

The Dynamics of Cytoplasmic mRNA Metabolism.

For all but a few mRNAs, the dynamics of metabolism are unknown. Here, we developed an experimental and analytical framework for examining these dynamics for mRNAs from thousands of genes. mRNAs of mouse fibroblasts exit the nucleus with diverse intragenic and intergenic poly(A)-tail lengths. Once in the cytoplasm, they have a broad (1000-fold) range of deadenylation rate constants, which correspond to cytoplasmic lifetimes. Indeed, with few exceptions, degradation appears to occur primarily through deadenylation-linked mechanisms, with little contribution from either endonucleolytic cleavage or deadenylation-independent decapping. Most mRNA molecules degrade only after their tail lengths fall below 25 nt. Decay rate constants of short-tailed mRNAs vary broadly (1000-fold) and are larger for short-tailed mRNAs that have previously undergone more rapid deadenylation. This coupling helps clear rapidly deadenylated mRNAs, enabling the large range in deadenylation rate constants to impart a similarly large range in stabilities.

Molecular, spatial, and functional single-cell profiling of the hypothalamic preoptic region.

The hypothalamus controls essential social behaviors and homeostatic functions. However, the cellular architecture of hypothalamic nuclei-including the molecular identity, spatial organization, and function of distinct cell types-is poorly understood. Here, we developed an imaging-based in situ cell-type identification and mapping method and combined it with single-cell RNA-sequencing to create a molecularly annotated and spatially resolved cell atlas of the mouse hypothalamic preoptic region. We profiled ~1 million cells, identified ~70 neuronal populations characterized by distinct neuromodulatory signatures and spatial organizations, and defined specific neuronal populations activated during social behaviors in male and female mice, providing a high-resolution framework for mechanistic investigation of behavior circuits. The approach described opens a new avenue for the construction of cell atlases in diverse tissues and organisms.

Genome-Scale Networks Link Neurodegenerative Disease Genes to α-Synuclein through Specific Molecular Pathways.

Numerous genes and molecular pathways are implicated in neurodegenerative proteinopathies, but their inter-relationships are poorly understood. We systematically mapped molecular pathways underlying the toxicity of alpha-synuclein (α-syn), a protein central to Parkinson's disease. Genome-wide screens in yeast identified 332 genes that impact α-syn toxicity. To "humanize" this molecular network, we developed a computational method, TransposeNet. This integrates a Steiner prize-collecting approach with homology assignment through sequence, structure, and interaction topology. TransposeNet linked α-syn to multiple parkinsonism genes and druggable targets through perturbed protein trafficking and ER quality control as well as mRNA metabolism and translation. A calcium signaling hub linked these processes to perturbed mitochondrial quality control and function, metal ion transport, transcriptional regulation, and signal transduction. Parkinsonism gene interaction profiles spatially opposed in the network (ATP13A2/PARK9 and VPS35/PARK17) were highly distinct, and network relationships for specific genes (LRRK2/PARK8, ATXN2, and EIF4G1/PARK18) were confirmed in patient induced pluripotent stem cell (iPSC)-derived neurons. This cross-species platform connected diverse neurodegenerative genes to proteinopathy through specific mechanisms and may facilitate patient stratification for targeted therapy.

mRNA poly(A)-tail changes specified by deadenylation broadly reshape translation in Drosophila oocytes and early embryos.

Because maturing oocytes and early embryos lack appreciable transcription, posttranscriptional regulatory processes control their development. To better understand this control, we profiled translational efficiencies and poly(A)-tail lengths throughout Drosophila oocyte maturation and early embryonic development. The correspondence between translational-efficiency changes and tail-length changes indicated that tail-length changes broadly regulate translation until gastrulation, when this coupling disappears. During egg activation, relative changes in poly(A)-tail length, and thus translational efficiency, were largely retained in the absence of cytoplasmic polyadenylation, which indicated that selective poly(A)-tail shortening primarily specifies these changes. Many translational changes depended on PAN GU and Smaug, and these changes were largely attributable to tail-length changes. Our results also revealed the presence of tail-length-independent mechanisms that maintained translation despite tail-length shortening during oocyte maturation, and prevented essentially all translation of bicoid and several other mRNAs before egg activation. In addition to these fundamental insights, our results provide valuable resources for future studies.

Improved Ribosome-Footprint and mRNA Measurements Provide Insights into Dynamics and Regulation of Yeast Translation.

Ribosome-footprint profiling provides genome-wide snapshots of translation, but technical challenges can confound its analysis. Here, we use improved methods to obtain ribosome-footprint profiles and mRNA abundances that more faithfully reflect gene expression in Saccharomyces cerevisiae. Our results support proposals that both the beginning of coding regions and codons matching rare tRNAs are more slowly translated. They also indicate that emergent polypeptides with as few as three basic residues within a ten-residue window tend to slow translation. With the improved mRNA measurements, the variation attributable to translational control in exponentially growing yeast was less than previously reported, and most of this variation could be predicted with a simple model that considered mRNA abundance, upstream open reading frames, cap-proximal structure and nucleotide composition, and lengths of the coding and 5' UTRs. Collectively, our results provide a framework for executing and interpreting ribosome-profiling studies and reveal key features of translational control in yeast.

mRNA destabilization is the dominant effect of mammalian microRNAs by the time substantial repression ensues.

MicroRNAs (miRNAs) regulate target mRNAs through a combination of translational repression and mRNA destabilization, with mRNA destabilization dominating at steady state in the few contexts examined globally. Here, we extend the global steady-state measurements to additional mammalian contexts and find that regardless of the miRNA, cell type, growth condition, or translational state, mRNA destabilization explains most (66%->90%) miRNA-mediated repression. We also determine the relative dynamics of translational repression and mRNA destabilization for endogenous mRNAs as a miRNA is induced. Although translational repression occurs rapidly, its effect is relatively weak, such that by the time consequential repression ensues, the effect of mRNA destabilization dominates. These results imply that consequential miRNA-mediated repression is largely irreversible and provide other insights into the nature of miRNA-mediated regulation. They also simplify future studies, dramatically extending the known contexts and time points for which monitoring mRNA changes captures most of the direct miRNA effects.

Widespread changes in the posttranscriptional landscape at the Drosophila oocyte-to-embryo transition.

The oocyte-to-embryo transition marks the onset of development. The initial phase of this profound change from the differentiated oocyte to the totipotent embryo occurs in the absence of both transcription and mRNA degradation. Here we combine global polysome profiling, ribosome-footprint profiling, and quantitative mass spectrometry in a comprehensive approach to delineate the translational and proteomic changes that occur during this important transition in Drosophila. Our results show that PNG kinase is a critical regulator of the extensive changes in the translatome, acting uniquely at this developmental window. Analysis of the proteome in png mutants provided insights into the contributions of translation to changes in protein levels, revealing a compensatory dynamic between translation and protein turnover during proteome remodeling at the return to totipotency. The proteome changes additionally suggested regulators of meiosis and early embryogenesis, including the conserved H3K4 demethylase LID, which we demonstrated is required during this period despite transcriptional inactivity.

Altered translation of GATA1 in Diamond-Blackfan anemia.

Ribosomal protein haploinsufficiency occurs in diverse human diseases including Diamond-Blackfan anemia (DBA), congenital asplenia and T cell leukemia. Yet, how mutations in genes encoding ubiquitously expressed proteins such as these result in cell-type- and tissue-specific defects remains unknown. Here, we identify mutations in GATA1, encoding the critical hematopoietic transcription factor GATA-binding protein-1, that reduce levels of full-length GATA1 protein and cause DBA in rare instances. We show that ribosomal protein haploinsufficiency, the more common cause of DBA, can lead to decreased GATA1 mRNA translation, possibly resulting from a higher threshold for initiation of translation of this mRNA in comparison with other mRNAs. In primary hematopoietic cells from patients with mutations in RPS19, encoding ribosomal protein S19, the amplitude of a transcriptional signature of GATA1 target genes was globally and specifically reduced, indicating that the activity, but not the mRNA level, of GATA1 is decreased in patients with DBA associated with mutations affecting ribosomal proteins. Moreover, the defective hematopoiesis observed in patients with DBA associated with ribosomal protein haploinsufficiency could be partially overcome by increasing GATA1 protein levels. Our results provide a paradigm by which selective defects in translation due to mutations affecting ubiquitous ribosomal proteins can result in human disease.

Poly(A)-tail profiling reveals an embryonic switch in translational control.

Poly(A) tails enhance the stability and translation of most eukaryotic messenger RNAs, but difficulties in globally measuring poly(A)-tail lengths have impeded greater understanding of poly(A)-tail function. Here we describe poly(A)-tail length profiling by sequencing (PAL-seq) and apply it to measure tail lengths of millions of individual RNAs isolated from yeasts, cell lines, Arabidopsis thaliana leaves, mouse liver, and zebrafish and frog embryos. Poly(A)-tail lengths were conserved between orthologous mRNAs, with mRNAs encoding ribosomal proteins and other 'housekeeping' proteins tending to have shorter tails. As expected, tail lengths were coupled to translational efficiencies in early zebrafish and frog embryos. However, this strong coupling diminished at gastrulation and was absent in non-embryonic samples, indicating a rapid developmental switch in the nature of translational control. This switch complements an earlier switch to zygotic transcriptional control and explains why the predominant effect of microRNA-mediated deadenylation concurrently shifts from translational repression to mRNA destabilization.

MicroRNA-134 activity in somatostatin interneurons regulates H-Ras localization by repressing the palmitoylation enzyme, DHHC9.

MicroRNA-134 (miR-134) serves as a widely accepted model for microRNA function in synaptic plasticity. In this model, synaptic activity stimulates miR-134 expression, which then regulates dendrite growth and spine formation. By using a ratiometric microRNA sensor, we found, unexpectedly, that miR-134 activity in cortical neurons was restricted to interneurons. Using an assay designed to trap microRNA-mRNA complexes, we determined that miR-134 interacted directly with the mRNA encoding the palmitoylation enzyme, DHHC9. This enzyme is known to palmitoylate H-Ras, a modification required for proper membrane trafficking. Treatment with bicuculline, a GABAA receptor antagonist, decreased DHHC9 expression in somatostatin-positive interneurons and membrane localization of an H-Ras reporter in a manner that depended on miR-134. Thus, although miR-134 has been proposed to affect all types of neurons, we showed that functionally active miR-134 is produced in only a selected population of neurons where it influences the expression of targets, such as DHHC9, that regulate membrane targeting of critical signaling molecules.