Search

Predicting elemental cycles and maintaining water quality under increasing anthropogenic influence requires knowledge of the spatial drivers of river microbiomes. However, understanding of the core microbial processes governing river biogeochemistry is hindered by a lack of genome-resolved functional insights and sampling across multiple rivers. Here we used a community science effort to accelerate the sampling, sequencing and genome-resolved analyses of river microbiomes to create the Genome Resolved Open Watersheds database (GROWdb). GROWdb profiles the identity, distribution, function and expression of microbial genomes across river surface waters covering 90% of United States watersheds. Specifically, GROWdb encompasses microbial lineages from 27 phyla, including novel members from 10 families and 128 genera, and defines the core river microbiome at the genome level. GROWdb analyses coupled to extensive geospatial information reveals local and regional drivers of microbial community structuring, while also presenting foundational hypotheses about ecosystem function. Building on the previously conceived River Continuum Concept1, we layer on microbial functional trait expression, which suggests that the structure and function of river microbiomes is predictable. We make GROWdb available through various collaborative cyberinfrastructures2,3, so that it can be widely accessed across disciplines for watershed predictive modelling and microbiome-based management practices.

Sterile inflammation after myocardial infarction is classically credited to myeloid cells interacting with dead cell debris in the infarct zone1,2. Here we show that cardiomyocytes are the dominant initiators of a previously undescribed type I interferon response in the infarct borderzone. Using spatial transcriptomics analysis in mice and humans, we find that myocardial infarction induces colonies of interferon-induced cells (IFNICs) expressing interferon-stimulated genes decorating the borderzone, where cardiomyocytes experience mechanical stress, nuclear rupture and escape of chromosomal DNA. Cardiomyocyte-selective deletion of Irf3 abrogated IFNIC colonies, whereas mice lacking Irf3 in fibroblasts, macrophages, neutrophils or endothelial cells, Ccr2-deficient mice or plasmacytoid-dendritic-cell-depleted mice did not. Interferons blunted the protective matricellular programs and contractile function of borderzone fibroblasts, and increased vulnerability to pathological remodelling. In mice that died after myocardial infarction, IFNIC colonies were immediately adjacent to sites of ventricular rupture, while mice lacking IFNICs were protected from rupture and exhibited improved survival3. Together, these results reveal a pathological borderzone niche characterized by a cardiomyocyte-initiated innate immune response. We suggest that selective inhibition of IRF3 activation in non-immune cells could limit ischaemic cardiomyopathy while avoiding broad immunosuppression.

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/ ).

Developmental myelination is a protracted process in the mammalian brain1. One theory for why oligodendrocytes mature so slowly posits that myelination may stabilize neuronal circuits and temper neuronal plasticity as animals age2-4. We tested this theory in the visual cortex, which has a well-defined critical period for experience-dependent neuronal plasticity5. During adolescence, visual experience modulated the rate of oligodendrocyte maturation in visual cortex. To determine whether oligodendrocyte maturation in turn regulates neuronal plasticity, we genetically blocked oligodendrocyte differentiation and myelination in adolescent mice. In adult mice lacking adolescent oligodendrogenesis, a brief period of monocular deprivation led to a significant decrease in visual cortex responses to the deprived eye, reminiscent of the plasticity normally restricted to adolescence. This enhanced functional plasticity was accompanied by a greater turnover of dendritic spines and coordinated reductions in spine size following deprivation. Furthermore, inhibitory synaptic transmission, which gates experience-dependent plasticity at the circuit level, was diminished in the absence of adolescent oligodendrogenesis. These results establish a critical role for oligodendrocytes in shaping the maturation and stabilization of cortical circuits and support the concept of developmental myelination acting as a functional brake on neuronal plasticity.

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.

Systemic lupus erythematosus (SLE) is prototypical autoimmune disease driven by pathological T cell-B cell interactions1,2. Expansion of T follicular helper (TFH) and T peripheral helper (TPH) cells, two T cell populations that provide help to B cells, is a prominent feature of SLE3,4. Human TFH and TPH cells characteristically produce high levels of the B cell chemoattractant CXCL13 (refs. 5,6), yet regulation of T cell CXCL13 production and the relationship between CXCL13+ T cells and other T cell states remains unclear. Here, we identify an imbalance in CD4+ T cell phenotypes in patients with SLE, with expansion of PD-1+/ICOS+ CXCL13+ T cells and reduction of CD96hi IL-22+ T cells. Using CRISPR screens, we identify the aryl hydrocarbon receptor (AHR) as a potent negative regulator of CXCL13 production by human CD4+ T cells. Transcriptomic, epigenetic and functional studies demonstrate that AHR coordinates with AP-1 family member JUN to prevent CXCL13+ TPH/TFH cell differentiation and promote an IL-22+ phenotype. Type I interferon, a pathogenic driver of SLE7, opposes AHR and JUN to promote T cell production of CXCL13. These results place CXCL13+ TPH/TFH cells on a polarization axis opposite from T helper 22 (TH22) cells and reveal AHR, JUN and interferon as key regulators of these divergent T cell states.

By directly altering microscopic interactions, pressure provides a powerful tuning knob for the exploration of condensed phases and geophysical phenomena1. The megabar regime represents an interesting frontier, in which recent discoveries include high-temperature superconductors, as well as structural and valence phase transitions2-6. However, at such high pressures, many conventional measurement techniques fail. Here we demonstrate the ability to perform local magnetometry inside a diamond anvil cell with sub-micron spatial resolution at megabar pressures. Our approach uses a shallow layer of nitrogen-vacancy colour centres implanted directly within the anvil7-9; crucially, we choose a crystal cut compatible with the intrinsic symmetries of the nitrogen-vacancy centre to enable functionality at megabar pressures. We apply our technique to characterize a recently discovered hydride superconductor, CeH9 (ref. 10). By performing simultaneous magnetometry and electrical transport measurements, we observe the dual signatures of superconductivity: diamagnetism characteristic of the Meissner effect and a sharp drop of the resistance to near zero. By locally mapping both the diamagnetic response and flux trapping, we directly image the geometry of superconducting regions, showing marked inhomogeneities at the micron scale. Our work brings quantum sensing to the megabar frontier and enables the closed-loop optimization of superhydride materials synthesis.

Electrified solid-liquid interfaces (ESLIs) play a key role in various electrochemical processes relevant to energy1-5, biology6 and geochemistry7. The electron and mass transport at the electrified interfaces may result in structural modifications that markedly influence the reaction pathways. For example, electrocatalyst surface restructuring during reactions can substantially affect the catalysis mechanisms and reaction products1-3. Despite its importance, direct probing the atomic dynamics of solid-liquid interfaces under electric biasing is challenging owing to the nature of being buried in liquid electrolytes and the limited spatial resolution of current techniques for in situ imaging through liquids. Here, with our development of advanced polymer electrochemical liquid cells for transmission electron microscopy (TEM), we are able to directly monitor the atomic dynamics of ESLIs during copper (Cu)-catalysed CO2 electroreduction reactions (CO2ERs). Our observation reveals a fluctuating liquid-like amorphous interphase. It undergoes reversible crystalline-amorphous structural transformations and flows along the electrified Cu surface, thus mediating the crystalline Cu surface restructuring and mass loss through the interphase layer. The combination of real-time observation and theoretical calculations unveils an amorphization-mediated restructuring mechanism resulting from charge-activated surface reactions with the electrolyte. Our results open many opportunities to explore the atomic dynamics and its impact in broad systems involving ESLIs by taking advantage of the in situ imaging capability.

Forests are a substantial terrestrial carbon sink, but anthropogenic changes in land use and climate have considerably reduced the scale of this system1. Remote-sensing estimates to quantify carbon losses from global forests2-5 are characterized by considerable uncertainty and we lack a comprehensive ground-sourced evaluation to benchmark these estimates. Here we combine several ground-sourced6 and satellite-derived approaches2,7,8 to evaluate the scale of the global forest carbon potential outside agricultural and urban lands. Despite regional variation, the predictions demonstrated remarkable consistency at a global scale, with only a 12% difference between the ground-sourced and satellite-derived estimates. At present, global forest carbon storage is markedly under the natural potential, with a total deficit of 226 Gt (model range = 151-363 Gt) in areas with low human footprint. Most (61%, 139 Gt C) of this potential is in areas with existing forests, in which ecosystem protection can allow forests to recover to maturity. The remaining 39% (87 Gt C) of potential lies in regions in which forests have been removed or fragmented. Although forests cannot be a substitute for emissions reductions, our results support the idea2,3,9 that the conservation, restoration and sustainable management of diverse forests offer valuable contributions to meeting global climate and biodiversity targets.

Correction to: Nature Published online 4 July 2022 In the version of this article initially published, the ATLAS Collaboration author names, affiliations and acknowledgements were omitted and have now been included in the HTML and PDF versions of the article.

The early evolution of a supernova (SN) can reveal information about the environment and the progenitor star. When a star explodes in vacuum, the first photons to escape from its surface appear as a brief, hours-long shock-breakout flare1,2, followed by a cooling phase of emission. However, for stars exploding within a distribution of dense, optically thick circumstellar material (CSM), the first photons escape from the material beyond the stellar edge and the duration of the initial flare can extend to several days, during which the escaping emission indicates photospheric heating3. Early serendipitous observations2,4 that lacked ultraviolet (UV) data were unable to determine whether the early emission is heating or cooling and hence the nature of the early explosion event. Here we report UV spectra of the nearby SN 2023ixf in the galaxy Messier 101 (M101). Using the UV data as well as a comprehensive set of further multiwavelength observations, we temporally resolve the emergence of the explosion shock from a thick medium heated by the SN emission. We derive a reliable bolometric light curve that indicates that the shock breaks out from a dense layer with a radius substantially larger than typical supergiants.

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.

Ecosystems generate a wide range of benefits for humans, including some market goods as well as other benefits that are not directly reflected in market activity1. Climate change will alter the distribution of ecosystems around the world and change the flow of these benefits2,3. However, the specific implications of ecosystem changes for human welfare remain unclear, as they depend on the nature of these changes, the value of the affected benefits and the extent to which communities rely on natural systems for their well-being4. Here we estimate country-level changes in economic production and the value of non-market ecosystem benefits resulting from climate-change-induced shifts in terrestrial vegetation cover, as projected by dynamic global vegetation models (DGVMs) driven by general circulation climate models. Our results show that the annual population-weighted mean global flow of non-market ecosystem benefits valued in the wealth accounts of the World Bank will be reduced by 9.2% in 2100 under the Shared Socioeconomic Pathway SSP2-6.0 with respect to the baseline no climate change scenario and that the global population-weighted average change in gross domestic product (GDP) by 2100 is -1.3% of the baseline GDP. Because lower-income countries are more reliant on natural capital, these GDP effects are regressive. Approximately 90% of these damages are borne by the poorest 50% of countries and regions, whereas the wealthiest 10% experience only 2% of these losses.

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.

In the published version of this paper, the list of members of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium and their affiliations contained minor errors in the affiliations. The original Article has been corrected to include the corrected list.

Patients from historically under-represented racial and ethnic groups are enrolled in cancer clinical trials at disproportionately low rates in the USA1-3. As these patients often have limited English proficiency4-7, we hypothesized that one barrier to their inclusion is the cost to investigators of translating consent documents. To test this hypothesis, we evaluated more than 12,000 consent events at a large cancer centre and assessed whether patients requiring translated consent documents would sign consent documents less frequently in studies lacking industry sponsorship (for which the principal investigator pays the translation costs) than for industry-sponsored studies (for which the translation costs are covered by the sponsor). Here we show that the proportion of consent events for patients with limited English proficiency in studies not sponsored by industry was approximately half of that seen in industry-sponsored studies. We also show that among those signing consent documents, the proportion of consent documents translated into the patient's primary language in studies without industry sponsorship was approximately half of that seen in industry-sponsored studies. The results suggest that the cost of consent document translation in trials not sponsored by industry could be a potentially modifiable barrier to the inclusion of patients with limited English proficiency.

The basic helix-loop-helix (bHLH) family of transcription factors recognizes DNA motifs known as E-boxes (CANNTG) and includes 108 members1. Here we investigate how chromatinized E-boxes are engaged by two structurally diverse bHLH proteins: the proto-oncogene MYC-MAX and the circadian transcription factor CLOCK-BMAL1 (refs. 2,3). Both transcription factors bind to E-boxes preferentially near the nucleosomal entry-exit sites. Structural studies with engineered or native nucleosome sequences show that MYC-MAX or CLOCK-BMAL1 triggers the release of DNA from histones to gain access. Atop the H2A-H2B acidic patch4, the CLOCK-BMAL1 Per-Arnt-Sim (PAS) dimerization domains engage the histone octamer disc. Binding of tandem E-boxes5-7 at endogenous DNA sequences occurs through direct interactions between two CLOCK-BMAL1 protomers and histones and is important for circadian cycling. At internal E-boxes, the MYC-MAX leucine zipper can also interact with histones H2B and H3, and its binding is indirectly enhanced by OCT4 elsewhere on the nucleosome. The nucleosomal E-box position and the type of bHLH dimerization domain jointly determine the histone contact, the affinity and the degree of competition and cooperativity with other nucleosome-bound factors.

Systems of correlated particles appear in many fields of modern science and represent some of the most intractable computational problems in nature. The computational challenge in these systems arises when interactions become comparable to other energy scales, which makes the state of each particle depend on all other particles1. The lack of general solutions for the three-body problem and acceptable theory for strongly correlated electrons shows that our understanding of correlated systems fades when the particle number or the interaction strength increases. One of the hallmarks of interacting systems is the formation of multiparticle bound states2-9. Here we develop a high-fidelity parameterizable fSim gate and implement the periodic quantum circuit of the spin-½ XXZ model in a ring of 24 superconducting qubits. We study the propagation of these excitations and observe their bound nature for up to five photons. We devise a phase-sensitive method for constructing the few-body spectrum of the bound states and extract their pseudo-charge by introducing a synthetic flux. By introducing interactions between the ring and additional qubits, we observe an unexpected resilience of the bound states to integrability breaking. This finding goes against the idea that bound states in non-integrable systems are unstable when their energies overlap with the continuum spectrum. Our work provides experimental evidence for bound states of interacting photons and discovers their stability beyond the integrability limit.

Chromothripsis, the shattering and imperfect reassembly of one (or a few) chromosome(s)1, is an ubiquitous2 mutational process generating localized and complex chromosomal rearrangements that drive genome evolution in cancer. Chromothripsis can be initiated by mis-segregation errors in mitosis3,4 or DNA metabolism5-7 that lead to entrapment of chromosomes within micronuclei and their subsequent fragmentation in the next interphase or following mitotic entry6,8-10. Here we use inducible degrons to demonstrate that chromothriptically produced pieces of a micronucleated chromosome are tethered together in mitosis by a protein complex consisting of mediator of DNA damage checkpoint 1 (MDC1), DNA topoisomerase II-binding protein 1 (TOPBP1) and cellular inhibitor of PP2A (CIP2A), thereby enabling en masse segregation to the same daughter cell. Such tethering is shown to be crucial for the viability of cells undergoing chromosome mis-segregation and shattering after transient inactivation of the spindle assembly checkpoint. Transient, degron-induced reduction in CIP2A following chromosome micronucleation-dependent chromosome shattering is shown to drive acquisition of segmental deletions and inversions. Analyses of pancancer tumour genomes showed that expression of CIP2A and TOPBP1 was increased overall in cancers with genomic rearrangements, including copy number-neutral chromothripsis with minimal deletions, but comparatively reduced in cancers with canonical chromothripsis in which deletions were frequent. Thus, chromatin-bound tethers maintain the proximity of fragments of a shattered chromosome enabling their re-encapsulation into, and religation within, a daughter cell nucleus to form heritable, chromothriptically rearranged chromosomes found in the majority of human cancers.

In the ongoing debates about eukaryogenesis-the series of evolutionary events leading to the emergence of the eukaryotic cell from prokaryotic ancestors-members of the Asgard archaea play a key part as the closest archaeal relatives of eukaryotes1. However, the nature and phylogenetic identity of the last common ancestor of Asgard archaea and eukaryotes remain unresolved2-4. Here we analyse distinct phylogenetic marker datasets of an expanded genomic sampling of Asgard archaea and evaluate competing evolutionary scenarios using state-of-the-art phylogenomic approaches. We find that eukaryotes are placed, with high confidence, as a well-nested clade within Asgard archaea and as a sister lineage to Hodarchaeales, a newly proposed order within Heimdallarchaeia. Using sophisticated gene tree and species tree reconciliation approaches, we show that analogous to the evolution of eukaryotic genomes, genome evolution in Asgard archaea involved significantly more gene duplication and fewer gene loss events compared with other archaea. Finally, we infer that the last common ancestor of Asgard archaea was probably a thermophilic chemolithotroph and that the lineage from which eukaryotes evolved adapted to mesophilic conditions and acquired the genetic potential to support a heterotrophic lifestyle. Our work provides key insights into the prokaryote-to-eukaryote transition and a platform for better understanding the emergence of cellular complexity in eukaryotic cells.

In natural photosynthesis, the light-driven splitting of water into electrons, protons and molecular oxygen forms the first step of the solar-to-chemical energy conversion process. The reaction takes place in photosystem II, where the Mn4CaO5 cluster first stores four oxidizing equivalents, the S0 to S4 intermediate states in the Kok cycle, sequentially generated by photochemical charge separations in the reaction center and then catalyzes the O-O bond formation chemistry1-3. Here, we report room temperature snapshots by serial femtosecond X-ray crystallography to provide structural insights into the final reaction step of Kok's photosynthetic water oxidation cycle, the S3→[S4]→S0 transition where O2 is formed and Kok's water oxidation clock is reset. Our data reveal a complex sequence of events, which occur over micro- to milliseconds, comprising changes at the Mn4CaO5 cluster, its ligands and water pathways as well as controlled proton release through the hydrogen-bonding network of the Cl1 channel. Importantly, the extra O atom Ox, which was introduced as a bridging ligand between Ca and Mn1 during the S2→S3 transition4-6, disappears or relocates in parallel with Yz reduction starting at approximately 700 μs after the third flash. The onset of O2 evolution, as indicated by the shortening of the Mn1-Mn4 distance, occurs at around 1,200 μs, signifying the presence of a reduced intermediate, possibly a bound peroxide.

The spatiotemporal structure of the human microbiome1,2, proteome3 and metabolome4,5 reflects and determines regional intestinal physiology and may have implications for disease6. Yet, little is known about the distribution of microorganisms, their environment and their biochemical activity in the gut because of reliance on stool samples and limited access to only some regions of the gut using endoscopy in fasting or sedated individuals7. To address these deficiencies, we developed an ingestible device that collects samples from multiple regions of the human intestinal tract during normal digestion. Collection of 240 intestinal samples from 15 healthy individuals using the device and subsequent multi-omics analyses identified significant differences between bacteria, phages, host proteins and metabolites in the intestines versus stool. Certain microbial taxa were differentially enriched and prophage induction was more prevalent in the intestines than in stool. The host proteome and bile acid profiles varied along the intestines and were highly distinct from those of stool. Correlations between gradients in bile acid concentrations and microbial abundance predicted species that altered the bile acid pool through deconjugation. Furthermore, microbially conjugated bile acid concentrations exhibited amino acid-dependent trends that were not apparent in stool. Overall, non-invasive, longitudinal profiling of microorganisms, proteins and bile acids along the intestinal tract under physiological conditions can help elucidate the roles of the gut microbiome and metabolome in human physiology and disease.

Gliomas synaptically integrate into neural circuits1,2. Previous research has demonstrated bidirectional interactions between neurons and glioma cells, with neuronal activity driving glioma growth1-4 and gliomas increasing neuronal excitability2,5-8. Here we sought to determine how glioma-induced neuronal changes influence neural circuits underlying cognition and whether these interactions influence patient survival. Using intracranial brain recordings during lexical retrieval language tasks in awake humans together with site-specific tumour tissue biopsies and cell biology experiments, we find that gliomas remodel functional neural circuitry such that task-relevant neural responses activate tumour-infiltrated cortex well beyond the cortical regions that are normally recruited in the healthy brain. Site-directed biopsies from regions within the tumour that exhibit high functional connectivity between the tumour and the rest of the brain are enriched for a glioblastoma subpopulation that exhibits a distinct synaptogenic and neuronotrophic phenotype. Tumour cells from functionally connected regions secrete the synaptogenic factor thrombospondin-1, which contributes to the differential neuron-glioma interactions observed in functionally connected tumour regions compared with tumour regions with less functional connectivity. Pharmacological inhibition of thrombospondin-1 using the FDA-approved drug gabapentin decreases glioblastoma proliferation. The degree of functional connectivity between glioblastoma and the normal brain negatively affects both patient survival and performance in language tasks. These data demonstrate that high-grade gliomas functionally remodel neural circuits in the human brain, which both promotes tumour progression and impairs cognition.

Here the Human Pangenome Reference Consortium presents a first draft of the human pangenome reference. The pangenome contains 47 phased, diploid assemblies from a cohort of genetically diverse individuals1. These assemblies cover more than 99% of the expected sequence in each genome and are more than 99% accurate at the structural and base pair levels. Based on alignments of the assemblies, we generate a draft pangenome that captures known variants and haplotypes and reveals new alleles at structurally complex loci. We also add 119 million base pairs of euchromatic polymorphic sequences and 1,115 gene duplications relative to the existing reference GRCh38. Roughly 90 million of the additional base pairs are derived from structural variation. Using our draft pangenome to analyse short-read data reduced small variant discovery errors by 34% and increased the number of structural variants detected per haplotype by 104% compared with GRCh38-based workflows, which enabled the typing of the vast majority of structural variant alleles per sample.

In quantum materials, degeneracies and frustrated interactions can have a profound impact on the emergence of long-range order, often driving strong fluctuations that suppress functionally relevant electronic or magnetic phases1-7. Engineering the atomic structure in the bulk or at heterointerfaces has been an important research strategy to lift these degeneracies, but these equilibrium methods are limited by thermodynamic, elastic and chemical constraints8. Here we show that all-optical, mode-selective manipulation of the crystal lattice can be used to enhance and stabilize high-temperature ferromagnetism in YTiO3, a material that shows only partial orbital polarization, an unsaturated low-temperature magnetic moment and a suppressed Curie temperature, Tc = 27 K (refs. 9-13). The enhancement is largest when exciting a 9 THz oxygen rotation mode, for which complete magnetic saturation is achieved at low temperatures and transient ferromagnetism is realized up to Tneq > 80 K, nearly three times the thermodynamic transition temperature. We interpret these effects as a consequence of the light-induced dynamical changes to the quasi-degenerate Ti t2g orbitals, which affect the magnetic phase competition and fluctuations found in the equilibrium state14-20. Notably, the light-induced high-temperature ferromagnetism discovered in our work is metastable over many nanoseconds, underscoring the ability to dynamically engineer practically useful non-equilibrium functionalities.

Oncogene amplification on extrachromosomal DNA (ecDNA) drives the evolution of tumours and their resistance to treatment, and is associated with poor outcomes for patients with cancer1-6. At present, it is unclear whether ecDNA is a later manifestation of genomic instability, or whether it can be an early event in the transition from dysplasia to cancer. Here, to better understand the development of ecDNA, we analysed whole-genome sequencing (WGS) data from patients with oesophageal adenocarcinoma (EAC) or Barrett's oesophagus. These data included 206 biopsies in Barrett's oesophagus surveillance and EAC cohorts from Cambridge University. We also analysed WGS and histology data from biopsies that were collected across multiple regions at 2 time points from 80 patients in a case-control study at the Fred Hutchinson Cancer Center. In the Cambridge cohorts, the frequency of ecDNA increased between Barrett's-oesophagus-associated early-stage (24%) and late-stage (43%) EAC, suggesting that ecDNA is formed during cancer progression. In the cohort from the Fred Hutchinson Cancer Center, 33% of patients who developed EAC had at least one oesophageal biopsy with ecDNA before or at the diagnosis of EAC. In biopsies that were collected before cancer diagnosis, higher levels of ecDNA were present in samples from patients who later developed EAC than in samples from those who did not. We found that ecDNAs contained diverse collections of oncogenes and immunomodulatory genes. Furthermore, ecDNAs showed increases in copy number and structural complexity at more advanced stages of disease. Our findings show that ecDNA can develop early in the transition from high-grade dysplasia to cancer, and that ecDNAs progressively form and evolve under positive selection.

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.

Monoterpene indole alkaloids (MIAs) are a diverse family of complex plant secondary metabolites with many medicinal properties, including the essential anti-cancer therapeutics vinblastine and vincristine1. As MIAs are difficult to chemically synthesize, the world's supply chain for vinblastine relies on low-yielding extraction and purification of the precursors vindoline and catharanthine from the plant Catharanthus roseus, which is then followed by simple in vitro chemical coupling and reduction to form vinblastine at an industrial scale2,3. Here, we demonstrate the de novo microbial biosynthesis of vindoline and catharanthine using a highly engineered yeast, and in vitro chemical coupling to vinblastine. The study showcases a very long biosynthetic pathway refactored into a microbial cell factory, including 30 enzymatic steps beyond the yeast native metabolites geranyl pyrophosphate and tryptophan to catharanthine and vindoline. In total, 56 genetic edits were performed, including expression of 34 heterologous genes from plants, as well as deletions, knock-downs and overexpression of ten yeast genes to improve precursor supplies towards de novo production of catharanthine and vindoline, from which semisynthesis to vinblastine occurs. As the vinblastine pathway is one of the longest MIA biosynthetic pathways, this study positions yeast as a scalable platform to produce more than 3,000 natural MIAs and a virtually infinite number of new-to-nature analogues.

When protons and neutrons (nucleons) are bound into atomic nuclei, they are close enough to feel significant attraction, or repulsion, from the strong, short-distance part of the nucleon-nucleon interaction. These strong interactions lead to hard collisions between nucleons, generating pairs of highly energetic nucleons referred to as short-range correlations (SRCs). SRCs are an important but relatively poorly understood part of nuclear structure1-3, and mapping out the strength and the isospin structure (neutron-proton (np) versus proton-proton (pp) pairs) of these virtual excitations is thus critical input for modelling a range of nuclear, particle and astrophysics measurements3-5. Two-nucleon knockout or 'triple coincidence' reactions have been used to measure the relative contribution of np-SRCs and pp-SRCs by knocking out a proton from the SRC and detecting its partner nucleon (proton or neutron). These measurements6-8 have shown that SRCs are almost exclusively np pairs, but they had limited statistics and required large model-dependent final-state interaction corrections. Here we report on measurements using inclusive scattering from the mirror nuclei hydrogen-3 and helium-3 to extract the np/pp ratio of SRCs in systems with a mass number of three. We obtain a measure of the np/pp SRC ratio that is an order of magnitude more precise than previous experiments, and find a marked deviation from the near-total np dominance observed in heavy nuclei. This result implies an unexpected structure in the high-momentum wavefunction for hydrogen-3 and helium-3. Understanding these results will improve our understanding of the short-range part of the nucleon-nucleon interaction.

Calcium imaging with protein-based indicators1,2 is widely used to follow neural activity in intact nervous systems, but current protein sensors report neural activity at timescales much slower than electrical signalling and are limited by trade-offs between sensitivity and kinetics. Here we used large-scale screening and structure-guided mutagenesis to develop and optimize several fast and sensitive GCaMP-type indicators3-8. The resulting 'jGCaMP8' sensors, based on the calcium-binding protein calmodulin and a fragment of endothelial nitric oxide synthase, have ultra-fast kinetics (half-rise times of 2 ms) and the highest sensitivity for neural activity reported for a protein-based calcium sensor. jGCaMP8 sensors will allow tracking of large populations of neurons on timescales relevant to neural computation.

Neoantigens are peptides derived from non-synonymous mutations presented by human leukocyte antigens (HLAs), which are recognized by antitumour T cells1-14. The large HLA allele diversity and limiting clinical samples have restricted the study of the landscape of neoantigen-targeted T cell responses in patients over their treatment course. Here we applied recently developed technologies15-17 to capture neoantigen-specific T cells from blood and tumours from patients with metastatic melanoma with or without response to anti-programmed death receptor 1 (PD-1) immunotherapy. We generated personalized libraries of neoantigen-HLA capture reagents to single-cell isolate the T cells and clone their T cell receptors (neoTCRs). Multiple T cells with different neoTCR sequences (T cell clonotypes) recognized a limited number of mutations in samples from seven patients with long-lasting clinical responses. These neoTCR clonotypes were recurrently detected over time in the blood and tumour. Samples from four patients with no response to anti-PD-1 also demonstrated neoantigen-specific T cell responses in the blood and tumour to a restricted number of mutations with lower TCR polyclonality and were not recurrently detected in sequential samples. Reconstitution of the neoTCRs in donor T cells using non-viral CRISPR-Cas9 gene editing demonstrated specific recognition and cytotoxicity to patient-matched melanoma cell lines. Thus, effective anti-PD-1 immunotherapy is associated with the presence of polyclonal CD8+ T cells in the tumour and blood specific for a limited number of immunodominant mutations, which are recurrently recognized over time.

Mitochondria are critical to the governance of metabolism and bioenergetics in cancer cells1. The mitochondria form highly organized networks, in which their outer and inner membrane structures define their bioenergetic capacity2,3. However, in vivo studies delineating the relationship between the structural organization of mitochondrial networks and their bioenergetic activity have been limited. Here we present an in vivo structural and functional analysis of mitochondrial networks and bioenergetic phenotypes in non-small cell lung cancer (NSCLC) using an integrated platform consisting of positron emission tomography imaging, respirometry and three-dimensional scanning block-face electron microscopy. The diverse bioenergetic phenotypes and metabolic dependencies we identified in NSCLC tumours align with distinct structural organization of mitochondrial networks present. Further, we discovered that mitochondrial networks are organized into distinct compartments within tumour cells. In tumours with high rates of oxidative phosphorylation (OXPHOSHI) and fatty acid oxidation, we identified peri-droplet mitochondrial networks wherein mitochondria contact and surround lipid droplets. By contrast, we discovered that in tumours with low rates of OXPHOS (OXPHOSLO), high glucose flux regulated perinuclear localization of mitochondria, structural remodelling of cristae and mitochondrial respiratory capacity. Our findings suggest that in NSCLC, mitochondrial networks are compartmentalized into distinct subpopulations that govern the bioenergetic capacity of tumours.

Carbon dioxide (CO2) is a key chemical species that is found in a wide range of planetary atmospheres. In the context of exoplanets, CO2 is an indicator of the metal enrichment (that is, elements heavier than helium, also called 'metallicity')1-3, and thus the formation processes of the primary atmospheres of hot gas giants4-6. It is also one of the most promising species to detect in the secondary atmospheres of terrestrial exoplanets7-9. Previous photometric measurements of transiting planets with the Spitzer Space Telescope have given hints of the presence of CO2, but have not yielded definitive detections owing to the lack of unambiguous spectroscopic identification10-12. Here we present the detection of CO2 in the atmosphere of the gas giant exoplanet WASP-39b from transmission spectroscopy observations obtained with JWST as part of the Early Release Science programme13,14. The data used in this study span 3.0-5.5 micrometres in wavelength and show a prominent CO2 absorption feature at 4.3 micrometres (26-sigma significance). The overall spectrum is well matched by one-dimensional, ten-times solar metallicity models that assume radiative-convective-thermochemical equilibrium and have moderate cloud opacity. These models predict that the atmosphere should have water, carbon monoxide and hydrogen sulfide in addition to CO2, but little methane. Furthermore, we also tentatively detect a small absorption feature near 4.0 micrometres that is not reproduced by these models.

Thwaites Glacier represents 15% of the ice discharge from the West Antarctic Ice Sheet and influences a wider catchment1-3. Because it is grounded below sea level4,5, Thwaites Glacier is thought to be susceptible to runaway retreat triggered at the grounding line (GL) at which the glacier reaches the ocean6,7. Recent ice-flow acceleration2,8 and retreat of the ice front8-10 and GL11,12 indicate that ice loss will continue. The relative impacts of mechanisms underlying recent retreat are however uncertain. Here we show sustained GL retreat from at least 2011 to 2020 and resolve mechanisms of ice-shelf melt at the submetre scale. Our conclusions are based on observations of the Thwaites Eastern Ice Shelf (TEIS) from an underwater vehicle, extending from the GL to 3 km oceanward and from the ice-ocean interface to the sea floor. These observations show a rough ice base above a sea floor sloping upward towards the GL and an ocean cavity in which the warmest water exceeds 2 °C above freezing. Data closest to the ice base show that enhanced melting occurs along sloped surfaces that initiate near the GL and evolve into steep-sided terraces. This pronounced melting along steep ice faces, including in crevasses, produces stratification that suppresses melt along flat interfaces. These data imply that slope-dependent melting sculpts the ice base and acts as an important response to ocean warming.

The global emergence of many severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants jeopardizes the protective antiviral immunity induced after infection or vaccination. To address the public health threat caused by the increasing SARS-CoV-2 genomic diversity, the National Institute of Allergy and Infectious Diseases within the National Institutes of Health established the SARS-CoV-2 Assessment of Viral Evolution (SAVE) programme. This effort was designed to provide a real-time risk assessment of SARS-CoV-2 variants that could potentially affect the transmission, virulence, and resistance to infection- and vaccine-induced immunity. The SAVE programme is a critical data-generating component of the US Government SARS-CoV-2 Interagency Group to assess implications of SARS-CoV-2 variants on diagnostics, vaccines and therapeutics, and for communicating public health risk. Here we describe the coordinated approach used to identify and curate data about emerging variants, their impact on immunity and effects on vaccine protection using animal models. We report the development of reagents, methodologies, models and notable findings facilitated by this collaborative approach and identify future challenges. This programme is a template for the response to rapidly evolving pathogens with pandemic potential by monitoring viral evolution in the human population to identify variants that could reduce the effectiveness of countermeasures.

In particle collider experiments, elementary particle interactions with large momentum transfer produce quarks and gluons (known as partons) whose evolution is governed by the strong force, as described by the theory of quantum chromodynamics (QCD)1. These partons subsequently emit further partons in a process that can be described as a parton shower2, which culminates in the formation of detectable hadrons. Studying the pattern of the parton shower is one of the key experimental tools for testing QCD. This pattern is expected to depend on the mass of the initiating parton, through a phenomenon known as the dead-cone effect, which predicts a suppression of the gluon spectrum emitted by a heavy quark of mass mQ and energy E, within a cone of angular size mQ/E around the emitter3. Previously, a direct observation of the dead-cone effect in QCD had not been possible, owing to the challenge of reconstructing the cascading quarks and gluons from the experimentally accessible hadrons. We report the direct observation of the QCD dead cone by using new iterative declustering techniques4,5 to reconstruct the parton shower of charm quarks. This result confirms a fundamental feature of QCD. Furthermore, the measurement of a dead-cone angle constitutes a direct experimental observation of the non-zero mass of the charm quark, which is a fundamental constant in the standard model of particle physics.

In the version of this article initially published, two members of the ENCODE Project Consortium were missing from the author list. Rizi Ai (Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, CA, USA) and Shantao Li (Program in Computational Biology and Bioinformatics, Yale University, New Haven, CT, USA) are now included in the author list. These errors have been corrected in the online version of the article.

Over the past few decades, neuroimaging has become a ubiquitous tool in basic research and clinical studies of the human brain. However, no reference standards currently exist to quantify individual differences in neuroimaging metrics over time, in contrast to growth charts for anthropometric traits such as height and weight1. Here we assemble an interactive open resource to benchmark brain morphology derived from any current or future sample of MRI data ( http://www.brainchart.io/ ). With the goal of basing these reference charts on the largest and most inclusive dataset available, acknowledging limitations due to known biases of MRI studies relative to the diversity of the global population, we aggregated 123,984 MRI scans, across more than 100 primary studies, from 101,457 human participants between 115 days post-conception to 100 years of age. MRI metrics were quantified by centile scores, relative to non-linear trajectories2 of brain structural changes, and rates of change, over the lifespan. Brain charts identified previously unreported neurodevelopmental milestones3, showed high stability of individuals across longitudinal assessments, and demonstrated robustness to technical and methodological differences between primary studies. Centile scores showed increased heritability compared with non-centiled MRI phenotypes, and provided a standardized measure of atypical brain structure that revealed patterns of neuroanatomical variation across neurological and psychiatric disorders. In summary, brain charts are an essential step towards robust quantification of individual variation benchmarked to normative trajectories in multiple, commonly used neuroimaging phenotypes.

Schizophrenia has a heritability of 60-80%1, much of which is attributable to common risk alleles. Here, in a two-stage genome-wide association study of up to 76,755 individuals with schizophrenia and 243,649 control individuals, we report common variant associations at 287 distinct genomic loci. Associations were concentrated in genes that are expressed in excitatory and inhibitory neurons of the central nervous system, but not in other tissues or cell types. Using fine-mapping and functional genomic data, we identify 120 genes (106 protein-coding) that are likely to underpin associations at some of these loci, including 16 genes with credible causal non-synonymous or untranslated region variation. We also implicate fundamental processes related to neuronal function, including synaptic organization, differentiation and transmission. Fine-mapped candidates were enriched for genes associated with rare disruptive coding variants in people with schizophrenia, including the glutamate receptor subunit GRIN2A and transcription factor SP4, and were also enriched for genes implicated by such variants in neurodevelopmental disorders. We identify biological processes relevant to schizophrenia pathophysiology; show convergence of common and rare variant associations in schizophrenia and neurodevelopmental disorders; and provide a resource of prioritized genes and variants to advance mechanistic studies.

Rare coding variation has historically provided the most direct connections between gene function and disease pathogenesis. By meta-analysing the whole exomes of 24,248 schizophrenia cases and 97,322 controls, we implicate ultra-rare coding variants (URVs) in 10 genes as conferring substantial risk for schizophrenia (odds ratios of 3-50, P

Tobacco and alcohol use are heritable behaviours associated with 15% and 5.3% of worldwide deaths, respectively, due largely to broad increased risk for disease and injury1-4. These substances are used across the globe, yet genome-wide association studies have focused largely on individuals of European ancestries5. Here we leveraged global genetic diversity across 3.4 million individuals from four major clines of global ancestry (approximately 21% non-European) to power the discovery and fine-mapping of genomic loci associated with tobacco and alcohol use, to inform function of these loci via ancestry-aware transcriptome-wide association studies, and to evaluate the genetic architecture and predictive power of polygenic risk within and across populations. We found that increases in sample size and genetic diversity improved locus identification and fine-mapping resolution, and that a large majority of the 3,823 associated variants (from 2,143 loci) showed consistent effect sizes across ancestry dimensions. However, polygenic risk scores developed in one ancestry performed poorly in others, highlighting the continued need to increase sample sizes of diverse ancestries to realize any potential benefit of polygenic prediction.

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.

The current human reference genome, GRCh38, represents over 20 years of effort to generate a high-quality assembly, which has benefitted society1,2. However, it still has many gaps and errors, and does not represent a biological genome as it is a blend of multiple individuals3,4. Recently, a high-quality telomere-to-telomere reference, CHM13, was generated with the latest long-read technologies, but it was derived from a hydatidiform mole cell line with a nearly homozygous genome5. To address these limitations, the Human Pangenome Reference Consortium formed with the goal of creating high-quality, cost-effective, diploid genome assemblies for a pangenome reference that represents human genetic diversity6. Here, in our first scientific report, we determined which combination of current genome sequencing and assembly approaches yield the most complete and accurate diploid genome assembly with minimal manual curation. Approaches that used highly accurate long reads and parent-child data with graph-based haplotype phasing during assembly outperformed those that did not. Developing a combination of the top-performing methods, we generated our first high-quality diploid reference assembly, containing only approximately four gaps per chromosome on average, with most chromosomes within ±1% of the length of CHM13. Nearly 48% of protein-coding genes have non-synonymous amino acid changes between haplotypes, and centromeric regions showed the highest diversity. Our findings serve as a foundation for assembling near-complete diploid human genomes at scale for a pangenome reference to capture global genetic variation from single nucleotides to structural rearrangements.

Infectious diseases are among the strongest selective pressures driving human evolution1,2. This includes the single greatest mortality event in recorded history, the first outbreak of the second pandemic of plague, commonly called the Black Death, which was caused by the bacterium Yersinia pestis3. This pandemic devastated Afro-Eurasia, killing up to 30-50% of the population4. To identify loci that may have been under selection during the Black Death, we characterized genetic variation around immune-related genes from 206 ancient DNA extracts, stemming from two different European populations before, during and after the Black Death. Immune loci are strongly enriched for highly differentiated sites relative to a set of non-immune loci, suggesting positive selection. We identify 245 variants that are highly differentiated within the London dataset, four of which were replicated in an independent cohort from Denmark, and represent the strongest candidates for positive selection. The selected allele for one of these variants, rs2549794, is associated with the production of a full-length (versus truncated) ERAP2 transcript, variation in cytokine response to Y. pestis and increased ability to control intracellular Y. pestis in macrophages. Finally, we show that protective variants overlap with alleles that are today associated with increased susceptibility to autoimmune diseases, providing empirical evidence for the role played by past pandemics in shaping present-day susceptibility to disease.

Common single-nucleotide polymorphisms (SNPs) are predicted to collectively explain 40-50% of phenotypic variation in human height, but identifying the specific variants and associated regions requires huge sample sizes1. Here, using data from a genome-wide association study of 5.4 million individuals of diverse ancestries, we show that 12,111 independent SNPs that are significantly associated with height account for nearly all of the common SNP-based heritability. These SNPs are clustered within 7,209 non-overlapping genomic segments with a mean size of around 90 kb, covering about 21% of the genome. The density of independent associations varies across the genome and the regions of increased density are enriched for biologically relevant genes. In out-of-sample estimation and prediction, the 12,111 SNPs (or all SNPs in the HapMap 3 panel2) account for 40% (45%) of phenotypic variance in populations of European ancestry but only around 10-20% (14-24%) in populations of other ancestries. Effect sizes, associated regions and gene prioritization are similar across ancestries, indicating that reduced prediction accuracy is likely to be explained by linkage disequilibrium and differences in allele frequency within associated regions. Finally, we show that the relevant biological pathways are detectable with smaller sample sizes than are needed to implicate causal genes and variants. Overall, this study provides a comprehensive map of specific genomic regions that contain the vast majority of common height-associated variants. Although this map is saturated for populations of European ancestry, further research is needed to achieve equivalent saturation in other ancestries.

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.

Multiple sclerosis (MS) is a heterogenous autoimmune disease in which autoreactive lymphocytes attack the myelin sheath of the central nervous system. B lymphocytes in the cerebrospinal fluid (CSF) of patients with MS contribute to inflammation and secrete oligoclonal immunoglobulins1,2. Epstein-Barr virus (EBV) infection has been epidemiologically linked to MS, but its pathological role remains unclear3. Here we demonstrate high-affinity molecular mimicry between the EBV transcription factor EBV nuclear antigen 1 (EBNA1) and the central nervous system protein glial cell adhesion molecule (GlialCAM) and provide structural and in vivo functional evidence for its relevance. A cross-reactive CSF-derived antibody was initially identified by single-cell sequencing of the paired-chain B cell repertoire of MS blood and CSF, followed by protein microarray-based testing of recombinantly expressed CSF-derived antibodies against MS-associated viruses. Sequence analysis, affinity measurements and the crystal structure of the EBNA1-peptide epitope in complex with the autoreactive Fab fragment enabled tracking of the development of the naive EBNA1-restricted antibody to a mature EBNA1-GlialCAM cross-reactive antibody. Molecular mimicry is facilitated by a post-translational modification of GlialCAM. EBNA1 immunization exacerbates disease in a mouse model of MS, and anti-EBNA1 and anti-GlialCAM antibodies are prevalent in patients with MS. Our results provide a mechanistic link for the association between MS and EBV and could guide the development of new MS therapies.

Medulloblastoma (MB) comprises a group of heterogeneous paediatric embryonal neoplasms of the hindbrain with strong links to early development of the hindbrain1-4. Mutations that activate Sonic hedgehog signalling lead to Sonic hedgehog MB in the upper rhombic lip (RL) granule cell lineage5-8. By contrast, mutations that activate WNT signalling lead to WNT MB in the lower RL9,10. However, little is known about the more commonly occurring group 4 (G4) MB, which is thought to arise in the unipolar brush cell lineage3,4. Here we demonstrate that somatic mutations that cause G4 MB converge on the core binding factor alpha (CBFA) complex and mutually exclusive alterations that affect CBFA2T2, CBFA2T3, PRDM6, UTX and OTX2. CBFA2T2 is expressed early in the progenitor cells of the cerebellar RL subventricular zone in Homo sapiens, and G4 MB transcriptionally resembles these progenitors but are stalled in developmental time. Knockdown of OTX2 in model systems relieves this differentiation blockade, which allows MB cells to spontaneously proceed along normal developmental differentiation trajectories. The specific nature of the split human RL, which is destined to generate most of the neurons in the human brain, and its high level of susceptible EOMES+KI67+ unipolar brush cell progenitor cells probably predisposes our species to the development of G4 MB.