Search

In lactating mothers, the high calcium (Ca2+) demand for milk production triggers significant bone loss1. Although oestrogen normally counteracts excessive bone resorption by promoting bone formation, this sex steroid drops precipitously during this postpartum period. Here we report that brain-derived cellular communication network factor 3 (CCN3) secreted from KISS1 neurons of the arcuate nucleus (ARCKISS1) fills this void and functions as a potent osteoanabolic factor to build bone in lactating females. We began by showing that our previously reported female-specific, dense bone phenotype2 originates from a humoral factor that promotes bone mass and acts on skeletal stem cells to increase their frequency and osteochondrogenic potential. This circulatory factor was then identified as CCN3, a brain-derived hormone from ARCKISS1 neurons that is able to stimulate mouse and human skeletal stem cell activity, increase bone remodelling and accelerate fracture repair in young and old mice of both sexes. The role of CCN3 in normal female physiology was revealed after detecting a burst of CCN3 expression in ARCKISS1 neurons coincident with lactation. After reducing CCN3 in ARCKISS1 neurons, lactating mothers lost bone and failed to sustain their progeny when challenged with a low-calcium diet. Our findings establish CCN3 as a potentially new therapeutic osteoanabolic hormone for both sexes and define a new maternal brain hormone for ensuring species survival in mammals.

Exaggerated airway constriction triggered by repeated exposure to allergen, also called hyperreactivity, is a hallmark of asthma. Whereas vagal sensory neurons are known to function in allergen-induced hyperreactivity1-3, the identity of downstream nodes remains poorly understood. Here we mapped a full allergen circuit from the lung to the brainstem and back to the lung. Repeated exposure of mice to inhaled allergen activated the nuclei of solitary tract (nTS) neurons in a mast cell-, interleukin-4 (IL-4)- and vagal nerve-dependent manner. Single-nucleus RNA sequencing, followed by RNAscope assay at baseline and allergen challenges, showed that a Dbh+ nTS population is preferentially activated. Ablation or chemogenetic inactivation of Dbh+ nTS neurons blunted hyperreactivity whereas chemogenetic activation promoted it. Viral tracing indicated that Dbh+ nTS neurons project to the nucleus ambiguus (NA) and that NA neurons are necessary and sufficient to relay allergen signals to postganglionic neurons that directly drive airway constriction. Delivery of noradrenaline antagonists to the NA blunted hyperreactivity, suggesting noradrenaline as the transmitter between Dbh+ nTS and NA. Together, these findings provide molecular, anatomical and functional definitions of key nodes of a canonical allergen response circuit. This knowledge informs how neural modulation could be used to control allergen-induced airway hyperreactivity.

Fibroblasts are present throughout the body and function to maintain tissue homeostasis. Recent studies have identified diverse fibroblast subsets in healthy and injured tissues1,2, but the origins and functional roles of injury-induced fibroblast lineages remain unclear. Here we show that lung-specialized alveolar fibroblasts take on multiple molecular states with distinct roles in facilitating responses to fibrotic lung injury. We generate a genetic tool that uniquely targets alveolar fibroblasts to demonstrate their role in providing niches for alveolar stem cells in homeostasis and show that loss of this niche leads to exaggerated responses to acute lung injury. Lineage tracing identifies alveolar fibroblasts as the dominant origin for multiple emergent fibroblast subsets sequentially driven by inflammatory and pro-fibrotic signals after injury. We identify similar, but not completely identical, fibroblast lineages in human pulmonary fibrosis. TGFβ negatively regulates an inflammatory fibroblast subset that emerges early after injury and stimulates the differentiation into fibrotic fibroblasts to elicit intra-alveolar fibrosis. Blocking the induction of fibrotic fibroblasts in the alveolar fibroblast lineage abrogates fibrosis but exacerbates lung inflammation. These results demonstrate the multifaceted roles of the alveolar fibroblast lineage in maintaining normal alveolar homeostasis and orchestrating sequential responses to lung injury.

Regular exercise promotes whole-body health and prevents disease, but the underlying molecular mechanisms are incompletely understood1-3. Here, the Molecular Transducers of Physical Activity Consortium4 profiled the temporal transcriptome, proteome, metabolome, lipidome, phosphoproteome, acetylproteome, ubiquitylproteome, epigenome and immunome in whole blood, plasma and 18 solid tissues in male and female Rattus norvegicus over eight weeks of endurance exercise training. The resulting data compendium encompasses 9,466 assays across 19 tissues, 25 molecular platforms and 4 training time points. Thousands of shared and tissue-specific molecular alterations were identified, with sex differences found in multiple tissues. Temporal multi-omic and multi-tissue analyses revealed expansive biological insights into the adaptive responses to endurance training, including widespread regulation of immune, metabolic, stress response and mitochondrial pathways. Many changes were relevant to human health, including non-alcoholic fatty liver disease, inflammatory bowel disease, cardiovascular health and tissue injury and recovery. The data and analyses presented in this study will serve as valuable resources for understanding and exploring the multi-tissue molecular effects of endurance training and are provided in a public repository ( https://motrpac-data.org/ ).

Astrocytes, the most abundant non-neuronal cell type in the mammalian brain, are crucial circuit components that respond to and modulate neuronal activity through calcium (Ca2+) signalling1-7. Astrocyte Ca2+ activity is highly heterogeneous and occurs across multiple spatiotemporal scales-from fast, subcellular activity3,4 to slow, synchronized activity across connected astrocyte networks8-10-to influence many processes5,7,11. However, the inputs that drive astrocyte network dynamics remain unclear. Here we used ex vivo and in vivo two-photon astrocyte imaging while mimicking neuronal neurotransmitter inputs at multiple spatiotemporal scales. We find that brief, subcellular inputs of GABA and glutamate lead to widespread, long-lasting astrocyte Ca2+ responses beyond an individual stimulated cell. Further, we find that a key subset of Ca2+ activity-propagative activity-differentiates astrocyte network responses to these two main neurotransmitters, and may influence responses to future inputs. Together, our results demonstrate that local, transient neurotransmitter inputs are encoded by broad cortical astrocyte networks over a minutes-long time course, contributing to accumulating evidence that substantial astrocyte-neuron communication occurs across slow, network-level spatiotemporal scales12-14. These findings will enable future studies to investigate the link between specific astrocyte Ca2+ activity and specific functional outputs, which could build a consistent framework for astrocytic modulation of neuronal activity.

Animals in the natural world constantly encounter geometrically complex landscapes. Successful navigation requires that they understand geometric features of these landscapes, including boundaries, landmarks, corners and curved areas, all of which collectively define the geometry of the environment1-12. Crucial to the reconstruction of the geometric layout of natural environments are concave and convex features, such as corners and protrusions. However, the neural substrates that could underlie the perception of concavity and convexity in the environment remain elusive. Here we show that the dorsal subiculum contains neurons that encode corners across environmental geometries in an allocentric reference frame. Using longitudinal calcium imaging in freely behaving mice, we find that corner cells tune their activity to reflect the geometric properties of corners, including corner angles, wall height and the degree of wall intersection. A separate population of subicular neurons encode convex corners of both larger environments and discrete objects. Both corner cells are non-overlapping with the population of subicular neurons that encode environmental boundaries. Furthermore, corner cells that encode concave or convex corners generalize their activity such that they respond, respectively, to concave or convex curvatures within an environment. Together, our findings suggest that the subiculum contains the geometric information needed to reconstruct the shape and layout of naturalistic spatial environments.

Determining the structure and phenotypic context of molecules detected in untargeted metabolomics experiments remains challenging. Here we present reverse metabolomics as a discovery strategy, whereby tandem mass spectrometry spectra acquired from newly synthesized compounds are searched for in public metabolomics datasets to uncover phenotypic associations. To demonstrate the concept, we broadly synthesized and explored multiple classes of metabolites in humans, including N-acyl amides, fatty acid esters of hydroxy fatty acids, bile acid esters and conjugated bile acids. Using repository-scale analysis1,2, we discovered that some conjugated bile acids are associated with inflammatory bowel disease (IBD). Validation using four distinct human IBD cohorts showed that cholic acids conjugated to Glu, Ile/Leu, Phe, Thr, Trp or Tyr are increased in Crohn's disease. Several of these compounds and related structures affected pathways associated with IBD, such as interferon-γ production in CD4+ T cells3 and agonism of the pregnane X receptor4. Culture of bacteria belonging to the Bifidobacterium, Clostridium and Enterococcus genera produced these bile amidates. Because searching repositories with tandem mass spectrometry spectra has only recently become possible, this reverse metabolomics approach can now be used as a general strategy to discover other molecules from human and animal ecosystems.

The temporal lobe of the human brain contains the entorhinal cortex (EC). This region of the brain is a highly interconnected integrative hub for sensory and spatial information; it also has a key role in episodic memory formation and is the main source of cortical hippocampal inputs1-4. The human EC continues to develop during childhood5, but neurogenesis and neuronal migration to the EC are widely considered to be complete by birth. Here we show that the human temporal lobe contains many young neurons migrating into the postnatal EC and adjacent regions, with a large tangential stream persisting until the age of around one year and radial dispersal continuing until around two to three years of age. By contrast, we found no equivalent postnatal migration in rhesus macaques (Macaca mulatta). Immunostaining and single-nucleus RNA sequencing of ganglionic eminence germinal zones, the EC stream and the postnatal EC revealed that most migrating cells in the EC stream are derived from the caudal ganglionic eminence and become LAMP5+RELN+ inhibitory interneurons. These late-arriving interneurons could continue to shape the processing of sensory and spatial information well into postnatal life, when children are actively interacting with their environment. The EC is one of the first regions of the brain to be affected in Alzheimers disease, and previous work has linked cognitive decline to the loss of LAMP5+RELN+ cells6,7. Our investigation reveals that many of these cells arrive in the EC through a major postnatal migratory stream in early childhood.

Cerebrospinal fluid (CSF) in the subarachnoid space around the brain has long been known to drain through the lymphatics to cervical lymph nodes1-17, but the connections and regulation have been challenging to identify. Here, using fluorescent CSF tracers in Prox1-GFP lymphatic reporter mice18, we found that the nasopharyngeal lymphatic plexus is a major hub for CSF outflow to deep cervical lymph nodes. This plexus had unusual valves and short lymphangions but no smooth-muscle coverage, whereas downstream deep cervical lymphatics had typical semilunar valves, long lymphangions and smooth muscle coverage that transported CSF to the deep cervical lymph nodes. α-Adrenergic and nitric oxide signalling in the smooth muscle cells regulated CSF drainage through the transport properties of deep cervical lymphatics. During ageing, the nasopharyngeal lymphatic plexus atrophied, but deep cervical lymphatics were not similarly altered, and CSF outflow could still be increased by adrenergic or nitric oxide signalling. Single-cell analysis of gene expression in lymphatic endothelial cells of the nasopharyngeal plexus of aged mice revealed increased type I interferon signalling and other inflammatory cytokines. The importance of evidence for the nasopharyngeal lymphatic plexus functioning as a CSF outflow hub is highlighted by its regression during ageing. Yet, the ageing-resistant pharmacological activation of deep cervical lymphatic transport towards lymph nodes can still increase CSF outflow, offering an approach for augmenting CSF clearance in age-related neurological conditions in which greater efflux would be beneficial.

Medium- and high-entropy alloys (M/HEAs) mix several principal elements with near-equiatomic composition and represent a model-shift strategy for designing previously unknown materials in metallurgy1-8, catalysis9-14 and other fields15-18. One of the core hypotheses of M/HEAs is lattice distortion5,19,20, which has been investigated by different numerical and experimental techniques21-26. However, determining the three-dimensional (3D) lattice distortion in M/HEAs remains a challenge. Moreover, the presumed random elemental mixing in M/HEAs has been questioned by X-ray and neutron studies27, atomistic simulations28-30, energy dispersive spectroscopy31,32 and electron diffraction33,34, which suggest the existence of local chemical order in M/HEAs. However, direct experimental observation of the 3D local chemical order has been difficult because energy dispersive spectroscopy integrates the composition of atomic columns along the zone axes7,32,34 and diffuse electron reflections may originate from planar defects instead of local chemical order35. Here we determine the 3D atomic positions of M/HEA nanoparticles using atomic electron tomography36 and quantitatively characterize the local lattice distortion, strain tensor, twin boundaries, dislocation cores and chemical short-range order (CSRO). We find that the high-entropy alloys have larger local lattice distortion and more heterogeneous strain than the medium-entropy alloys and that strain is correlated to CSRO. We also observe CSRO-mediated twinning in the medium-entropy alloys, that is, twinning occurs in energetically unfavoured CSRO regions but not in energetically favoured CSRO ones, which represents, to our knowledge, the first experimental observation of correlating local chemical order with structural defects in any material. We expect that this work will not only expand our fundamental understanding of this important class of materials but also provide the foundation for tailoring M/HEA properties through engineering lattice distortion and local chemical order.

Recent advances in single-cell technologies have led to the discovery of thousands of brain cell types; however, our understanding of the gene regulatory programs in these cell types is far from complete1-4. Here we report a comprehensive atlas of candidate cis-regulatory DNA elements (cCREs) in the adult mouse brain, generated by analysing chromatin accessibility in 2.3 million individual brain cells from 117 anatomical dissections. The atlas includes approximately 1 million cCREs and their chromatin accessibility across 1,482 distinct brain cell populations, adding over 446,000 cCREs to the most recent such annotation in the mouse genome. The mouse brain cCREs are moderately conserved in the human brain. The mouse-specific cCREs-specifically, those identified from a subset of cortical excitatory neurons-are strongly enriched for transposable elements, suggesting a potential role for transposable elements in the emergence of new regulatory programs and neuronal diversity. Finally, we infer the gene regulatory networks in over 260 subclasses of mouse brain cells and develop deep-learning models to predict the activities of gene regulatory elements in different brain cell types from the DNA sequence alone. Our results provide a resource for the analysis of cell-type-specific gene regulation programs in both mouse and human brains.

The termination of a meal is controlled by dedicated neural circuits in the caudal brainstem. A key challenge is to understand how these circuits transform the sensory signals generated during feeding into dynamic control of behaviour. The caudal nucleus of the solitary tract (cNTS) is the first site in the brain where many meal-related signals are sensed and integrated1-4, but how the cNTS processes ingestive feedback during behaviour is unknown. Here we describe how prolactin-releasing hormone (PRLH) and GCG neurons, two principal cNTS cell types that promote non-aversive satiety, are regulated during ingestion. PRLH neurons showed sustained activation by visceral feedback when nutrients were infused into the stomach, but these sustained responses were substantially reduced during oral consumption. Instead, PRLH neurons shifted to a phasic activity pattern that was time-locked to ingestion and linked to the taste of food. Optogenetic manipulations revealed that PRLH neurons control the duration of seconds-timescale feeding bursts, revealing a mechanism by which orosensory signals feed back to restrain the pace of ingestion. By contrast, GCG neurons were activated by mechanical feedback from the gut, tracked the amount of food consumed and promoted satiety that lasted for tens of minutes. These findings reveal that sequential negative feedback signals from the mouth and gut engage distinct circuits in the caudal brainstem, which in turn control elements of feeding behaviour operating on short and long timescales.

Cytosine DNA methylation is essential in brain development and is implicated in various neurological disorders. Understanding DNA methylation diversity across the entire brain in a spatial context is fundamental for a complete molecular atlas of brain cell types and their gene regulatory landscapes. Here we used single-nucleus methylome sequencing (snmC-seq3) and multi-omic sequencing (snm3C-seq)1 technologies to generate 301,626 methylomes and 176,003 chromatin conformation-methylome joint profiles from 117 dissected regions throughout the adult mouse brain. Using iterative clustering and integrating with companion whole-brain transcriptome and chromatin accessibility datasets, we constructed a methylation-based cell taxonomy with 4,673 cell groups and 274 cross-modality-annotated subclasses. We identified 2.6 million differentially methylated regions across the genome that represent potential gene regulation elements. Notably, we observed spatial cytosine methylation patterns on both genes and regulatory elements in cell types within and across brain regions. Brain-wide spatial transcriptomics data validated the association of spatial epigenetic diversity with transcription and improved the anatomical mapping of our epigenetic datasets. Furthermore, chromatin conformation diversities occurred in important neuronal genes and were highly associated with DNA methylation and transcription changes. Brain-wide cell-type comparisons enabled the construction of regulatory networks that incorporate transcription factors, regulatory elements and their potential downstream gene targets. Finally, intragenic DNA methylation and chromatin conformation patterns predicted alternative gene isoform expression observed in a whole-brain SMART-seq2 dataset. Our study establishes a brain-wide, single-cell DNA methylome and 3D multi-omic atlas and provides a valuable resource for comprehending the cellular-spatial and regulatory genome diversity of the mouse brain.

Single-cell analyses parse the brains billions of neurons into thousands of cell-type clusters residing in different brain structures1. Many cell types mediate their functions through targeted long-distance projections allowing interactions between specific cell types. Here we used epi-retro-seq2 to link single-cell epigenomes and cell types to long-distance projections for 33,034 neurons dissected from 32 different regions projecting to 24 different targets (225 source-to-target combinations) across the whole mouse brain. We highlight uses of these data for interrogating principles relating projection types to transcriptomics and epigenomics, and for addressing hypotheses about cell types and connections related to genetics. We provide an overall synthesis with 926 statistical comparisons of discriminability of neurons projecting to each target for every source. We integrate this dataset into the larger BRAIN Initiative Cell Census Network atlas, composed of millions of neurons, to link projection cell types to consensus clusters. Integration with spatial transcriptomics further assigns projection-enriched clusters to smaller source regions than the original dissections. We exemplify this by presenting in-depth analyses of projection neurons from the hypothalamus, thalamus, hindbrain, amygdala and midbrain to provide insights into properties of those cell types, including differentially expressed genes, their associated cis-regulatory elements and transcription-factor-binding motifs, and neurotransmitter use.

The basic plan of the retina is conserved across vertebrates, yet species differ profoundly in their visual needs1. Retinal cell types may have evolved to accommodate these varied needs, but this has not been systematically studied. Here we generated and integrated single-cell transcriptomic atlases of the retina from 17 species: humans, two non-human primates, four rodents, three ungulates, opossum, ferret, tree shrew, a bird, a reptile, a teleost fish and a lamprey. We found high molecular conservation of the six retinal cell classes (photoreceptors, horizontal cells, bipolar cells, amacrine cells, retinal ganglion cells (RGCs) and Müller glia), with transcriptomic variation across species related to evolutionary distance. Major subclasses were also conserved, whereas variation among cell types within classes or subclasses was more pronounced. However, an integrative analysis revealed that numerous cell types are shared across species, based on conserved gene expression programmes that are likely to trace back to an early ancestral vertebrate. The degree of variation among cell types increased from the outer retina (photoreceptors) to the inner retina (RGCs), suggesting that evolution acts preferentially to shape the retinal output. Finally, we identified rodent orthologues of midget RGCs, which comprise more than 80% of RGCs in the human retina, subserve high-acuity vision, and were previously believed to be restricted to primates2. By contrast, the mouse orthologues have large receptive fields and comprise around 2% of mouse RGCs. Projections of both primate and mouse orthologous types are overrepresented in the thalamus, which supplies the primary visual cortex. We suggest that midget RGCs are not primate innovations, but are descendants of evolutionarily ancient types that decreased in size and increased in number as primates evolved, thereby facilitating high visual acuity and increased cortical processing of visual information.

The transfer of photosynthetically produced organic carbon from surface to mesopelagic waters draws carbon dioxide from the atmosphere1. However, current observation-based estimates disagree on the strength of this biological carbon pump (BCP)2. Earth system models (ESMs) also exhibit a large spread of BCP estimates, indicating limited representations of the known carbon export pathways3. Here we use several decades of hydrographic observations to produce a top-down estimate of the strength of the BCP with an inverse biogeochemical model that implicitly accounts for all known export pathways. Our estimate of total organic carbon (TOC) export at 73.4 m (model euphotic zone depth) is 15.00 ± 1.12 Pg C year-1, with only two-thirds reaching 100 m depth owing to rapid remineralization of organic matter in the upper water column. Partitioned by sequestration time below the euphotic zone, τ, the globally integrated organic carbon production rate with τ > 3 months is 11.09 ± 1.02 Pg C year-1, dropping to 8.25 ± 0.30 Pg C year-1 for τ > 1 year, with 81% contributed by the non-advective-diffusive vertical flux owing to sinking particles and vertically migrating zooplankton. Nevertheless, export of organic carbon by mixing and other fluid transport of dissolved matter and suspended particles remains regionally important for meeting the respiratory carbon demand. Furthermore, the temperature dependence of the sequestration efficiency inferred from our inversion suggests that future global warming may intensify the recycling of organic matter in the upper ocean, potentially weakening the BCP.

Animal studies show aging varies between individuals as well as between organs within an individual1-4, but whether this is true in humans and its effect on age-related diseases is unknown. We utilized levels of human blood plasma proteins originating from specific organs to measure organ-specific aging differences in living individuals. Using machine learning models, we analysed aging in 11 major organs and estimated organ age reproducibly in five independent cohorts encompassing 5,676 adults across the human lifespan. We discovered nearly 20% of the population show strongly accelerated age in one organ and 1.7% are multi-organ agers. Accelerated organ aging confers 20-50% higher mortality risk, and organ-specific diseases relate to faster aging of those organs. We find individuals with accelerated heart aging have a 250% increased heart failure risk and accelerated brain and vascular aging predict Alzheimers disease (AD) progression independently from and as strongly as plasma pTau-181 (ref. 5), the current best blood-based biomarker for AD. Our models link vascular calcification, extracellular matrix alterations and synaptic protein shedding to early cognitive decline. We introduce a simple and interpretable method to study organ aging using plasma proteomics data, predicting diseases and aging effects.

Mouse models are a critical tool for studying human diseases, particularly developmental disorders1. However, conventional approaches for phenotyping may fail to detect subtle defects throughout the developing mouse2. Here we set out to establish single-cell RNA sequencing of the whole embryo as a scalable platform for the systematic phenotyping of mouse genetic models. We applied combinatorial indexing-based single-cell RNA sequencing3 to profile 101 embryos of 22 mutant and 4 wild-type genotypes at embryonic day 13.5, altogether profiling more than 1.6 million nuclei. The 22 mutants represent a range of anticipated phenotypic severities, from established multisystem disorders to deletions of individual regulatory regions4,5. We developed and applied several analytical frameworks for detecting differences in composition and/or gene expression across 52 cell types or trajectories. Some mutants exhibit changes in dozens of trajectories whereas others exhibit changes in only a few cell types. We also identify differences between widely used wild-type strains, compare phenotyping of gain- versus loss-of-function mutants and characterize deletions of topological associating domain boundaries. Notably, some changes are shared among mutants, suggesting that developmental pleiotropy might be 'decomposable' through further scaling of this approach. Overall, our findings show how single-cell profiling of whole embryos can enable the systematic molecular and cellular phenotypic characterization of mouse mutants with unprecedented breadth and resolution.

Decades of previous efforts to develop renal-sparing polyene antifungals were misguided by the classic membrane permeabilization model1. Recently, the clinically vital but also highly renal-toxic small-molecule natural product amphotericin B was instead found to kill fungi primarily by forming extramembraneous sponge-like aggregates that extract ergosterol from lipid bilayers2-6. Here we show that rapid and selective extraction of fungal ergosterol can yield potent and renal-sparing polyene antifungals. Cholesterol extraction was found to drive the toxicity of amphotericin B to human renal cells. Our examination of high-resolution structures of amphotericin B sponges in sterol-free and sterol-bound states guided us to a promising structural derivative that does not bind cholesterol and is thus renal sparing. This derivative was also less potent because it extracts ergosterol more slowly. Selective acceleration of ergosterol extraction with a second structural modification yielded a new polyene, AM-2-19, that is renal sparing in mice and primary human renal cells, potent against hundreds of pathogenic fungal strains, resistance evasive following serial passage in vitro and highly efficacious in animal models of invasive fungal infections. Thus, rational tuning of the dynamics of interactions between small molecules may lead to better treatments for fungal infections that still kill millions of people annually7,8 and potentially other resistance-evasive antimicrobials, including those that have recently been shown to operate through supramolecular structures that target specific lipids9.

To maintain a stable and clear image of the world, our eyes reflexively follow the direction in which a visual scene is moving. Such gaze-stabilization mechanisms reduce image blur as we move in the environment. In non-primate mammals, this behaviour is initiated by retinal output neurons called ON-type direction-selective ganglion cells (ON-DSGCs), which detect the direction of image motion and transmit signals to brainstem nuclei that drive compensatory eye movements1. However, ON-DSGCs have not yet been identified in the retina of primates, raising the possibility that this reflex is mediated by cortical visual areas. Here we mined single-cell RNA transcriptomic data from primate retina to identify a candidate ON-DSGC. We then combined two-photon calcium imaging, molecular identification and morphological analysis to reveal a population of ON-DSGCs in the macaque retina. The morphology, molecular signature and GABA (γ-aminobutyric acid)-dependent mechanisms that underlie direction selectivity in primate ON-DSGCs are highly conserved with those in other mammals. We further identify a candidate ON-DSGC in human retina. The presence of ON-DSGCs in primates highlights the need to examine the contribution of subcortical retinal mechanisms to normal and aberrant gaze stabilization in the developing and mature visual system.

The DNA damage response is essential to safeguard genome integrity. Although the contribution of chromatin in DNA repair has been investigated1,2, the contribution of chromosome folding to these processes remains unclear3. Here we report that, after the production of double-stranded breaks (DSBs) in mammalian cells, ATM drives the formation of a new chromatin compartment (D compartment) through the clustering of damaged topologically associating domains, decorated with γH2AX and 53BP1. This compartment forms by a mechanism that is consistent with polymer-polymer phase separation rather than liquid-liquid phase separation. The D compartment arises mostly in G1 phase, is independent of cohesin and is enhanced after pharmacological inhibition of DNA-dependent protein kinase (DNA-PK) or R-loop accumulation. Importantly, R-loop-enriched DNA-damage-responsive genes physically localize to the D compartment, and this contributes to their optimal activation, providing a function for DSB clustering in the DNA damage response. However, DSB-induced chromosome reorganization comes at the expense of an increased rate of translocations, also observed in cancer genomes. Overall, we characterize how DSB-induced compartmentalization orchestrates the DNA damage response and highlight the critical impact of chromosome architecture in genomic instability.

The molecular mechanisms and evolutionary changes accompanying synapse development are still poorly understood1,2. Here we generate a cross-species proteomic map of synapse development in the human, macaque and mouse neocortex. By tracking the changes of more than 1,000 postsynaptic density (PSD) proteins from midgestation to young adulthood, we find that PSD maturation in humans separates into three major phases that are dominated by distinct pathways. Cross-species comparisons reveal that human PSDs mature about two to three times slower than those of other species and contain higher levels of Rho guanine nucleotide exchange factors (RhoGEFs) in the perinatal period. Enhancement of RhoGEF signalling in human neurons delays morphological maturation of dendritic spines and functional maturation of synapses, potentially contributing to the neotenic traits of human brain development. In addition, PSD proteins can be divided into four modules that exert stage- and cell-type-specific functions, possibly explaining their differential associations with cognitive functions and diseases. Our proteomic map of synapse development provides a blueprint for studying the molecular basis and evolutionary changes of synapse maturation.

Craniosynostosis is a group of disorders of premature calvarial suture fusion. The identity of the calvarial stem cells (CSCs) that produce fusion-driving osteoblasts in craniosynostosis remains poorly understood. Here we show that both physiologic calvarial mineralization and pathologic calvarial fusion in craniosynostosis reflect the interaction of two separate stem cell lineages; a previously identified cathepsin K (CTSK) lineage CSC1 (CTSK+ CSC) and a separate discoidin domain-containing receptor 2 (DDR2) lineage stem cell (DDR2+ CSC) that we identified in this study. Deletion of Twist1, a gene associated with craniosynostosis in humans2,3, solely in CTSK+ CSCs is sufficient to drive craniosynostosis in mice, but the sites that are destined to fuse exhibit an unexpected depletion of CTSK+ CSCs and a corresponding expansion of DDR2+ CSCs, with DDR2+ CSC expansion being a direct maladaptive response to CTSK+ CSC depletion. DDR2+ CSCs display full stemness features, and our results establish the presence of two distinct stem cell lineages in the sutures, with both populations contributing to physiologic calvarial mineralization. DDR2+ CSCs mediate a distinct form of endochondral ossification without the typical haematopoietic marrow formation. Implantation of DDR2+ CSCs into suture sites is sufficient to induce fusion, and this phenotype was prevented by co-transplantation of CTSK+ CSCs. Finally, the human counterparts of DDR2+ CSCs and CTSK+ CSCs display conserved functional properties in xenograft assays. The interaction between these two stem cell populations provides a new biologic interface for the modulation of calvarial mineralization and suture patency.

Social animals live and move through spaces shaped by the presence, motion and sensory cues of multiple other individuals1-6. Neural activity in the hippocampus is known to reflect spatial behaviour7-9 yet its study is lacking in such dynamic group settings, which are ubiquitous in natural environments. Here we studied hippocampal activity in groups of bats engaged in collective spatial behaviour. We find that, under spontaneous conditions, a robust spatial structure emerges at the group level whereby behaviour is anchored to specific locations, movement patterns and individual social preferences. Using wireless electrophysiological recordings from both stationary and flying bats, we find that many hippocampal neurons are tuned to key features of group dynamics. These include the presence or absence of a conspecific, but not typically of an object, at landing sites, shared spatial locations, individual identities and sensory signals that are broadcasted in the group setting. Finally, using wireless calcium imaging, we find that social responses are anatomically distributed and robustly represented at the population level. Combined, our findings reveal that hippocampal activity contains a rich representation of naturally emerging spatial behaviours in animal groups that could in turn support the complex feat of collective behaviour.

Identifying therapeutics to delay, and potentially reverse, age-related cognitive decline is critical in light of the increased incidence of dementia-related disorders forecasted in the growing older population1. Here we show that platelet factors transfer the benefits of young blood to the ageing brain. Systemic exposure of aged male mice to a fraction of blood plasma from young mice containing platelets decreased neuroinflammation in the hippocampus at the transcriptional and cellular level and ameliorated hippocampal-dependent cognitive impairments. Circulating levels of the platelet-derived chemokine platelet factor 4 (PF4) (also known as CXCL4) were elevated in blood plasma preparations of young mice and humans relative to older individuals. Systemic administration of exogenous PF4 attenuated age-related hippocampal neuroinflammation, elicited synaptic-plasticity-related molecular changes and improved cognition in aged mice. We implicate decreased levels of circulating pro-ageing immune factors and restoration of the ageing peripheral immune system in the beneficial effects of systemic PF4 on the aged brain. Mechanistically, we identified CXCR3 as a chemokine receptor that, in part, mediates the cellular, molecular and cognitive benefits of systemic PF4 on the aged brain. Together, our data identify platelet-derived factors as potential therapeutic targets to abate inflammation and rescue cognition in old age.

Sympathetic activation during cold exposure increases adipocyte thermogenesis via the expression of mitochondrial protein uncoupling protein 1 (UCP1)1. The propensity of adipocytes to express UCP1 is under a critical influence of the adipose microenvironment and varies between sexes and among various fat depots2-7. Here we report that mammary gland ductal epithelial cells in the adipose niche regulate cold-induced adipocyte UCP1 expression in female mouse subcutaneous white adipose tissue (scWAT). Single-cell RNA sequencing shows that glandular luminal epithelium subtypes express transcripts that encode secretory factors controlling adipocyte UCP1 expression under cold conditions. We term these luminal epithelium secretory factors 'mammokines'. Using 3D visualization of whole-tissue immunofluorescence, we reveal sympathetic nerve-ductal contact points. We show that mammary ducts activated by sympathetic nerves limit adipocyte UCP1 expression via the mammokine lipocalin 2. In vivo and ex vivo ablation of mammary duct epithelium enhance the cold-induced adipocyte thermogenic gene programme in scWAT. Since the mammary duct network extends throughout most of the scWAT in female mice, females show markedly less scWAT UCP1 expression, fat oxidation, energy expenditure and subcutaneous fat mass loss compared with male mice, implicating sex-specific roles of mammokines in adipose thermogenesis. These results reveal a role of sympathetic nerve-activated glandular epithelium in adipocyte UCP1 expression and suggest that mammary duct luminal epithelium has an important role in controlling glandular adiposity.

Chromosomal instability (CIN) and epigenetic alterations are characteristics of advanced and metastatic cancers1-4, but whether they are mechanistically linked is unknown. Here we show that missegregation of mitotic chromosomes, their sequestration in micronuclei5,6 and subsequent rupture of the micronuclear envelope7 profoundly disrupt normal histone post-translational modifications (PTMs), a phenomenon conserved across humans and mice, as well as in cancer and non-transformed cells. Some of the changes in histone PTMs occur because of the rupture of the micronuclear envelope, whereas others are inherited from mitotic abnormalities before the micronucleus is formed. Using orthogonal approaches, we demonstrate that micronuclei exhibit extensive differences in chromatin accessibility, with a strong positional bias between promoters and distal or intergenic regions, in line with observed redistributions of histone PTMs. Inducing CIN causes widespread epigenetic dysregulation, and chromosomes that transit in micronuclei experience heritable abnormalities in their accessibility long after they have been reincorporated into the primary nucleus. Thus, as well as altering genomic copy number, CIN promotes epigenetic reprogramming and heterogeneity in cancer.

Niche signals maintain stem cells in a prolonged quiescence or transiently activate them for proper regeneration1. Altering balanced niche signalling can lead to regenerative disorders. Melanocytic skin nevi in human often display excessive hair growth, suggesting hair stem cell hyperactivity. Here, using genetic mouse models of nevi2,3, we show that dermal clusters of senescent melanocytes drive epithelial hair stem cells to exit quiescence and change their transcriptome and composition, potently enhancing hair renewal. Nevus melanocytes activate a distinct secretome, enriched for signalling factors. Osteopontin, the leading nevus signalling factor, is both necessary and sufficient to induce hair growth. Injection of osteopontin or its genetic overexpression is sufficient to induce robust hair growth in mice, whereas germline and conditional deletions of either osteopontin or CD44, its cognate receptor on epithelial hair cells, rescue enhanced hair growth induced by dermal nevus melanocytes. Osteopontin is overexpressed in human hairy nevi, and it stimulates new growth of human hair follicles. Although broad accumulation of senescent cells, such as upon ageing or genotoxic stress, is detrimental for the regenerative capacity of tissue4, we show that signalling by senescent cell clusters can potently enhance the activity of adjacent intact stem cells and stimulate tissue renewal. This finding identifies senescent cells and their secretome as an attractive therapeutic target in regenerative disorders.

Pancreatic ductal adenocarcinoma (PDA) is a lethal disease notoriously resistant to therapy1,2. This is mediated in part by a complex tumour microenvironment3, low vascularity4, and metabolic aberrations5,6. Although altered metabolism drives tumour progression, the spectrum of metabolites used as nutrients by PDA remains largely unknown. Here we identified uridine as a fuel for PDA in glucose-deprived conditions by assessing how more than 175 metabolites impacted metabolic activity in 21 pancreatic cell lines under nutrient restriction. Uridine utilization strongly correlated with the expression of uridine phosphorylase 1 (UPP1), which we demonstrate liberates uridine-derived ribose to fuel central carbon metabolism and thereby support redox balance, survival and proliferation in glucose-restricted PDA cells. In PDA, UPP1 is regulated by KRAS-MAPK signalling and is augmented by nutrient restriction. Consistently, tumours expressed high UPP1 compared with non-tumoural tissues, and UPP1 expression correlated with poor survival in cohorts of patients with PDA. Uridine is available in the tumour microenvironment, and we demonstrated that uridine-derived ribose is actively catabolized in tumours. Finally, UPP1 deletion restricted the ability of PDA cells to use uridine and blunted tumour growth in immunocompetent mouse models. Our data identify uridine utilization as an important compensatory metabolic process in nutrient-deprived PDA cells, suggesting a novel metabolic axis for PDA therapy.

Changes in patterns of activity within the medial prefrontal cortex enable rodents, non-human primates and humans to update their behaviour to adapt to changes in the environment-for example, during cognitive tasks1-5. Parvalbumin-expressing inhibitory neurons in the medial prefrontal cortex are important for learning new strategies during a rule-shift task6-8, but the circuit interactions that switch prefrontal network dynamics from maintaining to updating task-related patterns of activity remain unknown. Here we describe a mechanism that links parvalbumin-expressing neurons, a new callosal inhibitory connection, and changes in task representations. Whereas nonspecifically inhibiting all callosal projections does not prevent mice from learning rule shifts or disrupt the evolution of activity patterns, selectively inhibiting only callosal projections of parvalbumin-expressing neurons impairs rule-shift learning, desynchronizes the gamma-frequency activity that is necessary for learning8 and suppresses the reorganization of prefrontal activity patterns that normally accompanies rule-shift learning. This dissociation reveals how callosal parvalbumin-expressing projections switch the operating mode of prefrontal circuits from maintenance to updating by transmitting gamma synchrony and gating the ability of other callosal inputs to maintain previously established neural representations. Thus, callosal projections originating from parvalbumin-expressing neurons represent a key circuit locus for understanding and correcting the deficits in behavioural flexibility and gamma synchrony that have been implicated in schizophrenia and related conditions9,10.

The hippocampus is a mammalian brain structure that expresses spatial representations1 and is crucial for navigation2,3. Navigation, in turn, intricately depends on locomotion; however, current accounts suggest a dissociation between hippocampal spatial representations and the details of locomotor processes. Specifically, the hippocampus is thought to represent mainly higher-order cognitive and locomotor variables such as position, speed and direction of movement4-7, whereas the limb movements that propel the animal can be computed and represented primarily in subcortical circuits, including the spinal cord, brainstem and cerebellum8-11. Whether hippocampal representations are actually decoupled from the detailed structure of locomotor processes remains unknown. To address this question, here we simultaneously monitored hippocampal spatial representations and ongoing limb movements underlying locomotion at fast timescales. We found that the forelimb stepping cycle in freely behaving rats is rhythmic and peaks at around 8 Hz during movement, matching the approximately 8 Hz modulation of hippocampal activity and spatial representations during locomotion12. We also discovered precisely timed coordination between the time at which the forelimbs touch the ground ('plant' times of the stepping cycle) and the hippocampal representation of space. Notably, plant times coincide with hippocampal representations that are closest to the actual position of the nose of the rat, whereas between these plant times, the hippocampal representation progresses towards possible future locations. This synchronization was specifically detectable when rats approached spatial decisions. Together, our results reveal a profound and dynamic coordination on a timescale of tens of milliseconds between central cognitive representations and peripheral motor processes. This coordination engages and disengages rapidly in association with cognitive demands and is well suited to support rapid information exchange between cognitive and sensory-motor circuits.

Motor cortex (M1) has been thought to form a continuous somatotopic homunculus extending down the precentral gyrus from foot to face representations1,2, despite evidence for concentric functional zones3 and maps of complex actions4. Here, using precision functional magnetic resonance imaging (fMRI) methods, we find that the classic homunculus is interrupted by regions with distinct connectivity, structure and function, alternating with effector-specific (foot, hand and mouth) areas. These inter-effector regions exhibit decreased cortical thickness and strong functional connectivity to each other, as well as to the cingulo-opercular network (CON), critical for action5 and physiological control6, arousal7, errors8 and pain9. This interdigitation of action control-linked and motor effector regions was verified in the three largest fMRI datasets. Macaque and pediatric (newborn, infant and child) precision fMRI suggested cross-species homologues and developmental precursors of the inter-effector system. A battery of motor and action fMRI tasks documented concentric effector somatotopies, separated by the CON-linked inter-effector regions. The inter-effectors lacked movement specificity and co-activated during action planning (coordination of hands and feet) and axial body movement (such as of the abdomen or eyebrows). These results, together with previous studies demonstrating stimulation-evoked complex actions4 and connectivity to internal organs10 such as the adrenal medulla, suggest that M1 is punctuated by a system for whole-body action planning, the somato-cognitive action network (SCAN). In M1, two parallel systems intertwine, forming an integrate-isolate pattern: effector-specific regions (foot, hand and mouth) for isolating fine motor control and the SCAN for integrating goals, physiology and body movement.

Skates are cartilaginous fish whose body plan features enlarged wing-like pectoral fins, enabling them to thrive in benthic environments1,2. However, the molecular underpinnings of this unique trait remain unclear. Here we investigate the origin of this phenotypic innovation by developing the little skate Leucoraja erinacea as a genomically enabled model. Analysis of a high-quality chromosome-scale genome sequence for the little skate shows that it preserves many ancestral jawed vertebrate features compared with other sequenced genomes, including numerous ancient microchromosomes. Combining genome comparisons with extensive regulatory datasets in developing fins-including gene expression, chromatin occupancy and three-dimensional conformation-we find skate-specific genomic rearrangements that alter the three-dimensional regulatory landscape of genes that are involved in the planar cell polarity pathway. Functional inhibition of planar cell polarity signalling resulted in a reduction in anterior fin size, confirming that this pathway is a major contributor to batoid fin morphology. We also identified a fin-specific enhancer that interacts with several hoxa genes, consistent with the redeployment of hox gene expression in anterior pectoral fins, and confirmed its potential to activate transcription in the anterior fin using zebrafish reporter assays. Our findings underscore the central role of genome reorganization and regulatory variation in the evolution of phenotypes, shedding light on the molecular origin of an enigmatic trait.

Emerging spatial technologies, including spatial transcriptomics and spatial epigenomics, are becoming powerful tools for profiling of cellular states in the tissue context1-5. However, current methods capture only one layer of omics information at a time, precluding the possibility of examining the mechanistic relationship across the central dogma of molecular biology. Here, we present two technologies for spatially resolved, genome-wide, joint profiling of the epigenome and transcriptome by cosequencing chromatin accessibility and gene expression, or histone modifications (H3K27me3, H3K27ac or H3K4me3) and gene expression on the same tissue section at near-single-cell resolution. These were applied to embryonic and juvenile mouse brain, as well as adult human brain, to map how epigenetic mechanisms control transcriptional phenotype and cell dynamics in tissue. Although highly concordant tissue features were identified by either spatial epigenome or spatial transcriptome we also observed distinct patterns, suggesting their differential roles in defining cell states. Linking epigenome to transcriptome pixel by pixel allows the uncovering of new insights in spatial epigenetic priming, differentiation and gene regulation within the tissue architecture. These technologies are of great interest in life science and biomedical research.

The evolution of new traits enables expansion into new ecological and behavioural niches. Nonetheless, demonstrated connections between divergence in protein structure, function and lineage-specific behaviours remain rare. Here we show that both octopus and squid use cephalopod-specific chemotactile receptors (CRs) to sense their respective marine environments, but structural adaptations in these receptors support the sensation of specific molecules suited to distinct physiological roles. We find that squid express ancient CRs that more closely resemble related nicotinic acetylcholine receptors, whereas octopuses exhibit a more recent expansion in CRs consistent with their elaborated taste by touch sensory system. Using a combination of genetic profiling, physiology and behavioural analyses, we identify the founding member of squid CRs that detects soluble bitter molecules that are relevant in ambush predation. We present the cryo-electron microscopy structure of a squid CR and compare this with octopus CRs1 and nicotinic receptors2. These analyses demonstrate an evolutionary transition from an ancestral aromatic cage that coordinates soluble neurotransmitters or tastants to a more recent octopus CR hydrophobic binding pocket that traps insoluble molecules to mediate contact-dependent chemosensation. Thus, our study provides a foundation for understanding how adaptation of protein structure drives the diversification of organismal traits and behaviour.

Chemotactile receptors (CRs) are a cephalopod-specific innovation that allow octopuses to explore the seafloor via taste by touch1. CRs diverged from nicotinic acetylcholine receptors to mediate contact-dependent chemosensation of insoluble molecules that do not readily diffuse in marine environments. Here we exploit octopus CRs to probe the structural basis of sensory receptor evolution. We present the cryo-electron microscopy structure of an octopus CR and compare it with nicotinic receptors to determine features that enable environmental sensation versus neurotransmission. Evolutionary, structural and biophysical analyses show that the channel architecture involved in cation permeation and signal transduction is conserved. By contrast, the orthosteric ligand-binding site is subject to diversifying selection, thereby mediating the detection of new molecules. Serendipitous findings in the cryo-electron microscopy structure reveal that the octopus CR ligand-binding pocket is exceptionally hydrophobic, enabling sensation of greasy compounds versus the small polar molecules detected by canonical neurotransmitter receptors. These discoveries provide a structural framework for understanding connections between evolutionary adaptations at the atomic level and the emergence of new organismal behaviour.

Mutations in a diverse set of driver genes increase the fitness of haematopoietic stem cells (HSCs), leading to clonal haematopoiesis1. These lesions are precursors for blood cancers2-6, but the basis of their fitness advantage remains largely unknown, partly owing to a paucity of large cohorts in which the clonal expansion rate has been assessed by longitudinal sampling. Here, to circumvent this limitation, we developed a method to infer the expansion rate from data from a single time point. We applied this method to 5,071 people with clonal haematopoiesis. A genome-wide association study revealed that a common inherited polymorphism in the TCL1A promoter was associated with a slower expansion rate in clonal haematopoiesis overall, but the effect varied by driver gene. Those carrying this protective allele exhibited markedly reduced growth rates or prevalence of clones with driver mutations in TET2, ASXL1, SF3B1 and SRSF2, but this effect was not seen in clones with driver mutations in DNMT3A. TCL1A was not expressed in normal or DNMT3A-mutated HSCs, but the introduction of mutations in TET2 or ASXL1 led to the expression of TCL1A protein and the expansion of HSCs in vitro. The protective allele restricted TCL1A expression and expansion of mutant HSCs, as did experimental knockdown of TCL1A expression. Forced expression of TCL1A promoted the expansion of human HSCs in vitro and mouse HSCs in vivo. Our results indicate that the fitness advantage of several commonly mutated driver genes in clonal haematopoiesis may be mediated by TCL1A activation.

Chronic liver disease is a major public health burden worldwide1. Although different aetiologies and mechanisms of liver injury exist, progression of chronic liver disease follows a common pathway of liver inflammation, injury and fibrosis2. Here we examined the association between clonal haematopoiesis of indeterminate potential (CHIP) and chronic liver disease in 214,563 individuals from 4 independent cohorts with whole-exome sequencing data (Framingham Heart Study, Atherosclerosis Risk in Communities Study, UK Biobank and Mass General Brigham Biobank). CHIP was associated with an increased risk of prevalent and incident chronic liver disease (odds ratio = 2.01, 95% confidence interval (95% CI) [1.46, 2.79]; P

Astrocytes and neurons extensively interact in the brain. Identifying astrocyte and neuron proteomes is essential for elucidating the protein networks that dictate their respective contributions to physiology and disease. Here we used cell-specific and subcompartment-specific proximity-dependent biotinylation1 to study the proteomes of striatal astrocytes and neurons in vivo. We evaluated cytosolic and plasma membrane compartments for astrocytes and neurons to discover how these cells differ at the protein level in their signalling machinery. We also assessed subcellular compartments of astrocytes, including end feet and fine processes, to reveal their subproteomes and the molecular basis of essential astrocyte signalling and homeostatic functions. Notably, SAPAP3 (encoded by Dlgap3), which is associated with obsessive-compulsive disorder (OCD) and repetitive behaviours2-8, was detected at high levels in striatal astrocytes and was enriched within specific astrocyte subcompartments where it regulated actin cytoskeleton organization. Furthermore, genetic rescue experiments combined with behavioural analyses and molecular assessments in a mouse model of OCD4 lacking SAPAP3 revealed distinct contributions of astrocytic and neuronal SAPAP3 to repetitive and anxiety-related OCD-like phenotypes. Our data define how astrocytes and neurons differ at the protein level and in their major signalling pathways. Moreover, they reveal how astrocyte subproteomes vary between physiological subcompartments and how both astrocyte and neuronal SAPAP3 mechanisms contribute to OCD phenotypes in mice. Our data indicate that therapeutic strategies that target both astrocytes and neurons may be useful to explore in OCD and potentially other brain disorders.

The APOBEC3 (A3) proteins are host antiviral cellular proteins that hypermutate the viral genome of diverse viral families. In retroviruses, this process requires A3 packaging into viral particles1-4. The lentiviruses encode a protein, Vif, that antagonizes A3 family members by targeting them for degradation. Diversification of A3 allows host escape from Vif whereas adaptations in Vif enable cross-species transmission of primate lentiviruses. How this 'molecular arms race' plays out at the structural level is unknown. Here, we report the cryogenic electron microscopy structure of human APOBEC3G (A3G) bound to HIV-1 Vif, and the hijacked cellular proteins that promote ubiquitin-mediated proteolysis. A small surface explains the molecular arms race, including a cross-species transmission event that led to the birth of HIV-1. Unexpectedly, we find that RNA is a molecular glue for the Vif-A3G interaction, enabling Vif to repress A3G by ubiquitin-dependent and -independent mechanisms. Our results suggest a model in which Vif antagonizes A3G by intercepting it in its most dangerous form for the virus-when bound to RNA and on the pathway to packaging-to prevent viral restriction. By engaging essential surfaces required for restriction, Vif exploits a vulnerability in A3G, suggesting a general mechanism by which RNA binding helps to position key residues necessary for viral antagonism of a host antiviral gene.

Emotional states influence bodily physiology, as exemplified in the top-down process by which anxiety causes faster beating of the heart1-3. However, whether an increased heart rate might itself induce anxiety or fear responses is unclear3-8. Physiological theories of emotion, proposed over a century ago, have considered that in general, there could be an important and even dominant flow of information from the body to the brain9. Here, to formally test this idea, we developed a noninvasive optogenetic pacemaker for precise, cell-type-specific control of cardiac rhythms of up to 900 beats per minute in freely moving mice, enabled by a wearable micro-LED harness and the systemic viral delivery of a potent pump-like channelrhodopsin. We found that optically evoked tachycardia potently enhanced anxiety-like behaviour, but crucially only in risky contexts, indicating that both central (brain) and peripheral (body) processes may be involved in the development of emotional states. To identify potential mechanisms, we used whole-brain activity screening and electrophysiology to find brain regions that were activated by imposed cardiac rhythms. We identified the posterior insular cortex as a potential mediator of bottom-up cardiac interoceptive processing, and found that optogenetic inhibition of this brain region attenuated the anxiety-like behaviour that was induced by optical cardiac pacing. Together, these findings reveal that cells of both the body and the brain must be considered together to understand the origins of emotional or affective states. More broadly, our results define a generalizable approach for noninvasive, temporally precise functional investigations of joint organism-wide interactions among targeted cells during behaviour.

Multiple sclerosis is a chronic inflammatory disease of the central nervous system1. Astrocytes are heterogeneous glial cells that are resident in the central nervous system and participate in the pathogenesis of multiple sclerosis and its model experimental autoimmune encephalomyelitis2,3. However, few unique surface markers are available for the isolation of astrocyte subsets, preventing their analysis and the identification of candidate therapeutic targets; these limitations are further amplified by the rarity of pathogenic astrocytes. Here, to address these challenges, we developed focused interrogation of cells by nucleic acid detection and sequencing (FIND-seq), a high-throughput microfluidic cytometry method that combines encapsulation of cells in droplets, PCR-based detection of target nucleic acids and droplet sorting to enable in-depth transcriptomic analyses of cells of interest at single-cell resolution. We applied FIND-seq to study the regulation of astrocytes characterized by the splicing-driven activation of the transcription factor XBP1, which promotes disease pathology in multiple sclerosis and experimental autoimmune encephalomyelitis4. Using FIND-seq in combination with conditional-knockout mice, in vivo CRISPR-Cas9-driven genetic perturbation studies and bulk and single-cell RNA sequencing analyses of samples from mouse experimental autoimmune encephalomyelitis and humans with multiple sclerosis, we identified a new role for the nuclear receptor NR3C2 and its corepressor NCOR2 in limiting XBP1-driven pathogenic astrocyte responses. In summary, we used FIND-seq to identify a therapeutically targetable mechanism that limits XBP1-driven pathogenic astrocyte responses. FIND-seq enables the investigation of previously inaccessible cells, including rare cell subsets defined by unique gene expression signatures or other nucleic acid markers.

Tissues derive ATP from two pathways-glycolysis and the tricarboxylic acid (TCA) cycle coupled to the electron transport chain. Most energy in mammals is produced via TCA metabolism1. In tumours, however, the absolute rates of these pathways remain unclear. Here we optimize tracer infusion approaches to measure the rates of glycolysis and the TCA cycle in healthy mouse tissues, Kras-mutant solid tumours, metastases and leukaemia. Then, given the rates of these two pathways, we calculate total ATP synthesis rates. We find that TCA cycle flux is suppressed in all five primary solid tumour models examined and is increased in lung metastases of breast cancer relative to primary orthotopic tumours. As expected, glycolysis flux is increased in tumours compared with healthy tissues (the Warburg effect2,3), but this increase is insufficient to compensate for low TCA flux in terms of ATP production. Thus, instead of being hypermetabolic, as commonly assumed, solid tumours generally produce ATP at a slower than normal rate. In mouse pancreatic cancer, this is accommodated by the downregulation of protein synthesis, one of this tissue's major energy costs. We propose that, as solid tumours develop, cancer cells shed energetically expensive tissue-specific functions, enabling uncontrolled growth despite a limited ability to produce ATP.

Diabetes represents a spectrum of disease in which metabolic dysfunction damages multiple organ systems including liver, kidneys and peripheral nerves1,2. Although the onset and progression of these co-morbidities are linked with insulin resistance, hyperglycaemia and dyslipidaemia3-7, aberrant non-essential amino acid (NEAA) metabolism also contributes to the pathogenesis of diabetes8-10. Serine and glycine are closely related NEAAs whose levels are consistently reduced in patients with metabolic syndrome10-14, but the mechanistic drivers and downstream consequences of this metabotype remain unclear. Low systemic serine and glycine are also emerging as a hallmark of macular and peripheral nerve disorders, correlating with impaired visual acuity and peripheral neuropathy15,16. Here we demonstrate that aberrant serine homeostasis drives serine and glycine deficiencies in diabetic mice, which can be diagnosed with a serine tolerance test that quantifies serine uptake and disposal. Mimicking these metabolic alterations in young mice by dietary serine or glycine restriction together with high fat intake markedly accelerates the onset of small fibre neuropathy while reducing adiposity. Normalization of serine by dietary supplementation and mitigation of dyslipidaemia with myriocin both alleviate neuropathy in diabetic mice, linking serine-associated peripheral neuropathy to sphingolipid metabolism. These findings identify systemic serine deficiency and dyslipidaemia as novel risk factors for peripheral neuropathy that may be exploited therapeutically.

To survive, animals must convert sensory information into appropriate behaviours1,2. Vision is a common sense for locating ethologically relevant stimuli and guiding motor responses3-5. How circuitry converts object location in retinal coordinates to movement direction in body coordinates remains largely unknown. Here we show through behaviour, physiology, anatomy and connectomics in Drosophila that visuomotor transformation occurs by conversion of topographic maps formed by the dendrites of feature-detecting visual projection neurons (VPNs)6,7 into synaptic weight gradients of VPN outputs onto central brain neurons. We demonstrate how this gradient motif transforms the anteroposterior location of a visual looming stimulus into the fly's directional escape. Specifically, we discover that two neurons postsynaptic to a looming-responsive VPN type promote opposite takeoff directions. Opposite synaptic weight gradients onto these neurons from looming VPNs in different visual field regions convert localized looming threats into correctly oriented escapes. For a second looming-responsive VPN type, we demonstrate graded responses along the dorsoventral axis. We show that this synaptic gradient motif generalizes across all 20 primary VPN cell types and most often arises without VPN axon topography. Synaptic gradients may thus be a general mechanism for conveying spatial features of sensory information into directed motor outputs.

Systems consolidation-a process for long-term memory stabilization-has been hypothesized to occur in two stages1-4. Whereas new memories require the hippocampus5-9, they become integrated into cortical networks over time10-12, making them independent of the hippocampus. How hippocampal-cortical dialogue precisely evolves during this and how cortical representations change in concert is unknown. Here, we use a skill learning task13,14 to monitor the dynamics of cross-area coupling during non-rapid eye movement sleep along with changes in primary motor cortex (M1) representational stability. Our results indicate that precise cross-area coupling between hippocampus, prefrontal cortex and M1 can demarcate two distinct stages of processing. We specifically find that each animal demonstrates a sharp increase in prefrontal cortex and M1 sleep slow oscillation coupling with stabilization of performance. This sharp increase then predicts a drop in hippocampal sharp-wave ripple (SWR)-M1 slow oscillation coupling-suggesting feedback to inform hippocampal disengagement and transition to a second stage. Notably, the first stage shows significant increases in hippocampal SWR-M1 slow oscillation coupling in the post-training sleep and is closely associated with rapid learning and variability of the M1 low-dimensional manifold. Strikingly, even after consolidation, inducing new manifold exploration by changing task parameters re-engages hippocampal-M1 coupling. We thus find evidence for dynamic hippocampal-cortical dialogue associated with manifold exploration during learning and adaptation.

Multilocular adipocytes are a hallmark of thermogenic adipose tissue1,2, but the factors that enforce this cellular phenotype are largely unknown. Here, we show that an adipocyte-selective product of the Clstn3 locus (CLSTN3β) present in only placental mammals facilitates the efficient use of stored triglyceride by limiting lipid droplet (LD) expansion. CLSTN3β is an integral endoplasmic reticulum (ER) membrane protein that localizes to ER-LD contact sites through a conserved hairpin-like domain. Mice lacking CLSTN3β have abnormal LD morphology and altered substrate use in brown adipose tissue, and are more susceptible to cold-induced hypothermia despite having no defect in adrenergic signalling. Conversely, forced expression of CLSTN3β is sufficient to enforce a multilocular LD phenotype in cultured cells and adipose tissue. CLSTN3β associates with cell death-inducing DFFA-like effector proteins and impairs their ability to transfer lipid between LDs, thereby restricting LD fusion and expansion. Functionally, increased LD surface area in CLSTN3β-expressing adipocytes promotes engagement of the lipolytic machinery and facilitates fatty acid oxidation. In human fat, CLSTN3B is a selective marker of multilocular adipocytes. These findings define a molecular mechanism that regulates LD form and function to facilitate lipid utilization in thermogenic adipocytes.

Immune-related adverse events, particularly severe toxicities such as myocarditis, are major challenges to the utility of immune checkpoint inhibitors (ICIs) in anticancer therapy1. The pathogenesis of ICI-associated myocarditis (ICI-MC) is poorly understood. Pdcd1-/-Ctla4+/- mice recapitulate clinicopathological features of ICI-MC, including myocardial T cell infiltration2. Here, using single-cell RNA and T cell receptor (TCR) sequencing of cardiac immune infiltrates from Pdcd1-/-Ctla4+/- mice, we identify clonal effector CD8+ T cells as the dominant cell population. Treatment with anti-CD8-depleting, but not anti-CD4-depleting, antibodies improved the survival of Pdcd1-/-Ctla4+/- mice. Adoptive transfer of immune cells from mice with myocarditis induced fatal myocarditis in recipients, which required CD8+ T cells. The cardiac-specific protein α-myosin, which is absent from the thymus3,4, was identified as the cognate antigen source for three major histocompatibility complex class I-restricted TCRs derived from mice with fulminant myocarditis. Peripheral blood T cells from three patients with ICI-MC were expanded by α-myosin peptides. Moreover, these α-myosin-expanded T cells shared TCR clonotypes with diseased heart and skeletal muscle, which indicates that α-myosin may be a clinically important autoantigen in ICI-MC. These studies underscore the crucial role for cytotoxic CD8+ T cells, identify a candidate autoantigen in ICI-MC and yield new insights into the pathogenesis of ICI toxicity.

There is considerable interest in screening ultralarge chemical libraries for ligand discovery, both empirically and computationally1-4. Efforts have focused on readily synthesizable molecules, inevitably leaving many chemotypes unexplored. Here we investigate structure-based docking of a bespoke virtual library of tetrahydropyridines-a scaffold that is poorly sampled by a general billion-molecule virtual library but is well suited to many aminergic G-protein-coupled receptors. Using three inputs, each with diverse available derivatives, a one pot C-H alkenylation, electrocyclization and reduction provides the tetrahydropyridine core with up to six sites of derivatization5-7. Docking a virtual library of 75 million tetrahydropyridines against a model of the serotonin 5-HT2A receptor (5-HT2AR) led to the synthesis and testing of 17 initial molecules. Four of these molecules had low-micromolar activities against either the 5-HT2A or the 5-HT2B receptors. Structure-based optimization led to the 5-HT2AR agonists (R)-69 and (R)-70, with half-maximal effective concentration values of 41 nM and 110 nM, respectively, and unusual signalling kinetics that differ from psychedelic 5-HT2AR agonists. Cryo-electron microscopy structural analysis confirmed the predicted binding mode to 5-HT2AR. The favourable physical properties of these new agonists conferred high brain permeability, enabling mouse behavioural assays. Notably, neither had psychedelic activity, in contrast to classic 5-HT2AR agonists, whereas both had potent antidepressant activity in mouse models and had the same efficacy as antidepressants such as fluoxetine at as low as 1/40th of the dose. Prospects for using bespoke virtual libraries to sample pharmacologically relevant chemical space will be considered.

As the United Nations develops a post-2020 global biodiversity framework for the Convention on Biological Diversity, attention is focusing on how new goals and targets for ecosystem conservation might serve its vision of 'living in harmony with nature'1,2. Advancing dual imperatives to conserve biodiversity and sustain ecosystem services requires reliable and resilient generalizations and predictions about ecosystem responses to environmental change and management3. Ecosystems vary in their biota4, service provision5 and relative exposure to risks6, yet there is no globally consistent classification of ecosystems that reflects functional responses to change and management. This hampers progress on developing conservation targets and sustainability goals. Here we present the International Union for Conservation of Nature (IUCN) Global Ecosystem Typology, a conceptually robust, scalable, spatially explicit approach for generalizations and predictions about functions, biota, risks and management remedies across the entire biosphere. The outcome of a major cross-disciplinary collaboration, this novel framework places all of Earth's ecosystems into a unifying theoretical context to guide the transformation of ecosystem policy and management from global to local scales. This new information infrastructure will support knowledge transfer for ecosystem-specific management and restoration, globally standardized ecosystem risk assessments, natural capital accounting and progress on the post-2020 global biodiversity framework.