Search

Life-threatening thrombotic events and neurological symptoms are prevalent in COVID-19 and are persistent in patients with long COVID experiencing post-acute sequelae of SARS-CoV-2 infection1-4. Despite the clinical evidence1,5-7, the underlying mechanisms of coagulopathy in COVID-19 and its consequences in inflammation and neuropathology remain poorly understood and treatment options are insufficient. Fibrinogen, the central structural component of blood clots, is abundantly deposited in the lungs and brains of patients with COVID-19, correlates with disease severity and is a predictive biomarker for post-COVID-19 cognitive deficits1,5,8-10. Here we show that fibrin binds to the SARS-CoV-2 spike protein, forming proinflammatory blood clots that drive systemic thromboinflammation and neuropathology in COVID-19. Fibrin, acting through its inflammatory domain, is required for oxidative stress and macrophage activation in the lungs, whereas it suppresses natural killer cells, after SARS-CoV-2 infection. Fibrin promotes neuroinflammation and neuronal loss after infection, as well as innate immune activation in the brain and lungs independently of active infection. A monoclonal antibody targeting the inflammatory fibrin domain provides protection from microglial activation and neuronal injury, as well as from thromboinflammation in the lung after infection. Thus, fibrin drives inflammation and neuropathology in SARS-CoV-2 infection, and fibrin-targeting immunotherapy may represent a therapeutic intervention for patients with acute COVID-19 and long COVID.

The genomes of living lungfishes can inform on the molecular-developmental basis of the Devonian sarcopterygian fish-tetrapod transition. We de novo sequenced the genomes of the African (Protopterus annectens) and South American lungfishes (Lepidosiren paradoxa). The Lepidosiren genome (about 91 Gb, roughly 30 times the human genome) is the largest animal genome sequenced so far and more than twice the size of the Australian (Neoceratodus forsteri)1 and African2 lungfishes owing to enlarged intergenic regions and introns with high repeat content (about 90%). All lungfish genomes continue to expand as some transposable elements (TEs) are still active today. In particular, Lepidosiren's genome grew extremely fast during the past 100 million years (Myr), adding the equivalent of one human genome every 10 Myr. This massive genome expansion seems to be related to a reduction of PIWI-interacting RNAs and C2H2 zinc-finger and Krüppel-associated box (KRAB)-domain protein genes that suppress TE expansions. Although TE abundance facilitates chromosomal rearrangements, lungfish chromosomes still conservatively reflect the ur-tetrapod karyotype. Neoceratodus' limb-like fins still resemble those of their extinct relatives and remained phenotypically static for about 100 Myr. We show that the secondary loss of limb-like appendages in the Lepidosiren-Protopterus ancestor was probably due to loss of sonic hedgehog limb-specific enhancers.

The photoelectric effect is not truly instantaneous but exhibits attosecond delays that can reveal complex molecular dynamics1-7. Sub-femtosecond-duration light pulses provide the requisite tools to resolve the dynamics of photoionization8-12. Accordingly, the past decade has produced a large volume of work on photoionization delays following single-photon absorption of an extreme ultraviolet photon. However, the measurement of time-resolved core-level photoionization remained out of reach. The required X-ray photon energies needed for core-level photoionization were not available with attosecond tabletop sources. Here we report measurements of the X-ray photoemission delay of core-level electrons, with unexpectedly large delays, ranging up to 700 as in NO near the oxygen K-shell threshold. These measurements exploit attosecond soft X-ray pulses from a free-electron laser to scan across the entire region near the K-shell threshold. Furthermore, we find that the delay spectrum is richly modulated, suggesting several contributions, including transient trapping of the photoelectron owing to shape resonances, collisions with the Auger-Meitner electron that is emitted in the rapid non-radiative relaxation of the molecule and multi-electron scattering effects. The results demonstrate how X-ray attosecond experiments, supported by comprehensive theoretical modelling, can unravel the complex correlated dynamics of core-level photoionization.

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.

Major migration events in Holocene Eurasia have been characterized genetically at broad regional scales1-4. However, insights into the population dynamics in the contact zones are hampered by a lack of ancient genomic data sampled at high spatiotemporal resolution5-7. Here, to address this, we analysed shotgun-sequenced genomes from 100 skeletons spanning 7,300 years of the Mesolithic period, Neolithic period and Early Bronze Age in Denmark and integrated these with proxies for diet (13C and 15N content), mobility (87Sr/86Sr ratio) and vegetation cover (pollen). We observe that Danish Mesolithic individuals of the Maglemose, Kongemose and Ertebølle cultures form a distinct genetic cluster related to other Western European hunter-gatherers. Despite shifts in material culture they displayed genetic homogeneity from around 10,500 to 5,900 calibrated years before present, when Neolithic farmers with Anatolian-derived ancestry arrived. Although the Neolithic transition was delayed by more than a millennium relative to Central Europe, it was very abrupt and resulted in a population turnover with limited genetic contribution from local hunter-gatherers. The succeeding Neolithic population, associated with the Funnel Beaker culture, persisted for only about 1,000 years before immigrants with eastern Steppe-derived ancestry arrived. This second and equally rapid population replacement gave rise to the Single Grave culture with an ancestry profile more similar to present-day Danes. In our multiproxy dataset, these major demographic events are manifested as parallel shifts in genotype, phenotype, diet and land use.

Numerous studies have shown reduced performance in plants that are surrounded by neighbours of the same species1,2, a phenomenon known as conspecific negative density dependence (CNDD)3. A long-held ecological hypothesis posits that CNDD is more pronounced in tropical than in temperate forests4,5, which increases community stabilization, species coexistence and the diversity of local tree species6,7. Previous analyses supporting such a latitudinal gradient in CNDD8,9 have suffered from methodological limitations related to the use of static data10-12. Here we present a comprehensive assessment of latitudinal CNDD patterns using dynamic mortality data to estimate species-site-specific CNDD across 23 sites. Averaged across species, we found that stabilizing CNDD was present at all except one site, but that average stabilizing CNDD was not stronger toward the tropics. However, in tropical tree communities, rare and intermediate abundant species experienced stronger stabilizing CNDD than did common species. This pattern was absent in temperate forests, which suggests that CNDD influences species abundances more strongly in tropical forests than it does in temperate ones13. We also found that interspecific variation in CNDD, which might attenuate its stabilizing effect on species diversity14,15, was high but not significantly different across latitudes. Although the consequences of these patterns for latitudinal diversity gradients are difficult to evaluate, we speculate that a more effective regulation of population abundances could translate into greater stabilization of tropical tree communities and thus contribute to the high local diversity of tropical forests.

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 Holocene (beginning around 12,000 years ago) encompassed some of the most significant changes in human evolution, with far-reaching consequences for the dietary, physical and mental health of present-day populations. Using a dataset of more than 1,600 imputed ancient genomes1, we modelled the selection landscape during the transition from hunting and gathering, to farming and pastoralism across West Eurasia. We identify key selection signals related to metabolism, including that selection at the FADS cluster began earlier than previously reported and that selection near the LCT locus predates the emergence of the lactase persistence allele by thousands of years. We also find strong selection in the HLA region, possibly due to increased exposure to pathogens during the Bronze Age. Using ancient individuals to infer local ancestry tracts in over 400,000 samples from the UK Biobank, we identify widespread differences in the distribution of Mesolithic, Neolithic and Bronze Age ancestries across Eurasia. By calculating ancestry-specific polygenic risk scores, we show that height differences between Northern and Southern Europe are associated with differential Steppe ancestry, rather than selection, and that risk alleles for mood-related phenotypes are enriched for Neolithic farmer ancestry, whereas risk alleles for diabetes and Alzheimer's disease are enriched for Western hunter-gatherer ancestry. Our results indicate that ancient selection and migration were large contributors to the distribution of phenotypic diversity in present-day Europeans.

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.

Close-in giant exoplanets with temperatures greater than 2,000 K (ultra-hot Jupiters) have been the subject of extensive efforts to determine their atmospheric properties using thermal emission measurements from the Hubble Space Telescope (HST) and Spitzer Space Telescope1-3. However, previous studies have yielded inconsistent results because the small sizes of the spectral features and the limited information content of the data resulted in high sensitivity to the varying assumptions made in the treatment of instrument systematics and the atmospheric retrieval analysis3-12. Here we present a dayside thermal emission spectrum of the ultra-hot Jupiter WASP-18b obtained with the NIRISS13 instrument on the JWST. The data span 0.85 to 2.85 μm in wavelength at an average resolving power of 400 and exhibit minimal systematics. The spectrum shows three water emission features (at >6σ confidence) and evidence for optical opacity, possibly attributable to H-, TiO and VO (combined significance of 3.8σ). Models that fit the data require a thermal inversion, molecular dissociation as predicted by chemical equilibrium, a solar heavy-element abundance (metallicity, [Formula: see text] times solar) and a carbon-to-oxygen (C/O) ratio less than unity. The data also yield a dayside brightness temperature map, which shows a peak in temperature near the substellar point that decreases steeply and symmetrically with longitude towards the terminators.

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.

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.

Photochemistry is a fundamental process of planetary atmospheres that regulates the atmospheric composition and stability1. However, no unambiguous photochemical products have been detected in exoplanet atmospheres so far. Recent observations from the JWST Transiting Exoplanet Community Early Release Science Program2,3 found a spectral absorption feature at 4.05 μm arising from sulfur dioxide (SO2) in the atmosphere of WASP-39b. WASP-39b is a 1.27-Jupiter-radii, Saturn-mass (0.28 MJ) gas giant exoplanet orbiting a Sun-like star with an equilibrium temperature of around 1,100 K (ref. 4). The most plausible way of generating SO2 in such an atmosphere is through photochemical processes5,6. Here we show that the SO2 distribution computed by a suite of photochemical models robustly explains the 4.05-μm spectral feature identified by JWST transmission observations7 with NIRSpec PRISM (2.7σ)8 and G395H (4.5σ)9. SO2 is produced by successive oxidation of sulfur radicals freed when hydrogen sulfide (H2S) is destroyed. The sensitivity of the SO2 feature to the enrichment of the atmosphere by heavy elements (metallicity) suggests that it can be used as a tracer of atmospheric properties, with WASP-39b exhibiting an inferred metallicity of about 10× solar. We further point out that SO2 also shows observable features at ultraviolet and thermal infrared wavelengths not available from the existing observations.

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.

Measuring the metallicity and carbon-to-oxygen (C/O) ratio in exoplanet atmospheres is a fundamental step towards constraining the dominant chemical processes at work and, if in equilibrium, revealing planet formation histories. Transmission spectroscopy (for example, refs. 1,2) provides the necessary means by constraining the abundances of oxygen- and carbon-bearing species; however, this requires broad wavelength coverage, moderate spectral resolution and high precision, which, together, are not achievable with previous observatories. Now that JWST has commenced science operations, we are able to observe exoplanets at previously uncharted wavelengths and spectral resolutions. Here we report time-series observations of the transiting exoplanet WASP-39b using JWSTs Near InfraRed Camera (NIRCam). The long-wavelength spectroscopic and short-wavelength photometric light curves span 2.0-4.0 micrometres, exhibit minimal systematics and reveal well defined molecular absorption features in the planets spectrum. Specifically, we detect gaseous water in the atmosphere and place an upper limit on the abundance of methane. The otherwise prominent carbon dioxide feature at 2.8 micrometres is largely masked by water. The best-fit chemical equilibrium models favour an atmospheric metallicity of 1-100-times solar (that is, an enrichment of elements heavier than helium relative to the Sun) and a substellar C/O ratio. The inferred high metallicity and low C/O ratio may indicate significant accretion of solid materials during planet formation (for example, refs. 3,4,) or disequilibrium processes in the upper atmosphere (for example, refs. 5,6).

Marine heatwaves have been linked to negative ecological effects in recent decades1,2. If marine heatwaves regularly induce community reorganization and biomass collapses in fishes, the consequences could be catastrophic for ecosystems, fisheries and human communities3,4. However, the extent to which marine heatwaves have negative impacts on fish biomass or community composition, or even whether their effects can be distinguished from natural and sampling variability, remains unclear. We investigated the effects of 248 sea-bottom heatwaves from 1993 to 2019 on marine fishes by analysing 82,322 hauls (samples) from long-term scientific surveys of continental shelf ecosystems in North America and Europe spanning the subtropics to the Arctic. Here we show that the effects of marine heatwaves on fish biomass were often minimal and could not be distinguished from natural and sampling variability. Furthermore, marine heatwaves were not consistently associated with tropicalization (gain of warm-affiliated species) or deborealization (loss of cold-affiliated species) in these ecosystems. Although steep declines in biomass occasionally occurred after marine heatwaves, these were the exception, not the rule. Against the highly variable backdrop of ocean ecosystems, marine heatwaves have not driven biomass change or community turnover in fish communities that support many of the world's largest and most productive fisheries.

Many critics raise concerns about the prevalence of echo chambers on social media and their potential role in increasing political polarization. However, the lack of available data and the challenges of conducting large-scale field experiments have made it difficult to assess the scope of the problem1,2. Here we present data from 2020 for the entire population of active adult Facebook users in the USA showing that content from like-minded sources constitutes the majority of what people see on the platform, although political information and news represent only a small fraction of these exposures. To evaluate a potential response to concerns about the effects of echo chambers, we conducted a multi-wave field experiment on Facebook among 23,377 users for whom we reduced exposure to content from like-minded sources during the 2020 US presidential election by about one-third. We found that the intervention increased their exposure to content from cross-cutting sources and decreased exposure to uncivil language, but had no measurable effects on eight preregistered attitudinal measures such as affective polarization, ideological extremity, candidate evaluations and belief in false claims. These precisely estimated results suggest that although exposure to content from like-minded sources on social media is common, reducing its prevalence during the 2020 US presidential election did not correspondingly reduce polarization in beliefs or attitudes.

The SARS-CoV-2 Delta (Pango lineage B.1.617.2) variant of concern spread globally, causing resurgences of COVID-19 worldwide1,2. The emergence of the Delta variant in the UK occurred on the background of a heterogeneous landscape of immunity and relaxation of non-pharmaceutical interventions. Here we analyse 52,992 SARS-CoV-2 genomes from England together with 93,649 genomes from the rest of the world to reconstruct the emergence of Delta and quantify its introduction to and regional dissemination across England in the context of changing travel and social restrictions. Using analysis of human movement, contact tracing and virus genomic data, we find that the geographic focus of the expansion of Delta shifted from India to a more global pattern in early May 2021. In England, Delta lineages were introduced more than 1,000 times and spread nationally as non-pharmaceutical interventions were relaxed. We find that hotel quarantine for travellers reduced onward transmission from importations; however, the transmission chains that later dominated the Delta wave in England were seeded before travel restrictions were introduced. Increasing inter-regional travel within England drove the nationwide dissemination of Delta, with some cities receiving more than 2,000 observable lineage introductions from elsewhere. Subsequently, increased levels of local population mixing-and not the number of importations-were associated with the faster relative spread of Delta. The invasion dynamics of Delta depended on spatial heterogeneity in contact patterns, and our findings will inform optimal spatial interventions to reduce the transmission of current and future variants of concern, such as Omicron (Pango lineage B.1.1.529).

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

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.

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.

Neurotransmitters play essential roles in regulating neural circuit dynamics both in the central nervous system as well as at the peripheral, including the gastrointestinal tract1-3. Their real-time monitoring will offer critical information for understanding neural function and diagnosing disease1-3. However, bioelectronic tools to monitor the dynamics of neurotransmitters in vivo, especially in the enteric nervous systems, are underdeveloped. This is mainly owing to the limited availability of biosensing tools that are capable of examining soft, complex and actively moving organs. Here we introduce a tissue-mimicking, stretchable, neurochemical biological interface termed NeuroString, which is prepared by laser patterning of a metal-complexed polyimide into an interconnected graphene/nanoparticle network embedded in an elastomer. NeuroString sensors allow chronic in vivo real-time, multichannel and multiplexed monoamine sensing in the brain of behaving mouse, as well as measuring serotonin dynamics in the gut without undesired stimulations and perturbing peristaltic movements. The described elastic and conformable biosensing interface has broad potential for studying the impact of neurotransmitters on gut microbes, brain-gut communication and may ultimately be extended to biomolecular sensing in other soft organs across the body.

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.

Measuring the abundances of carbon and oxygen in exoplanet atmospheres is considered a crucial avenue for unlocking the formation and evolution of exoplanetary systems1,2. Access to the chemical inventory of an exoplanet requires high-precision observations, often inferred from individual molecular detections with low-resolution space-based3-5 and high-resolution ground-based6-8 facilities. Here we report the medium-resolution (R ≈ 600) transmission spectrum of an exoplanet atmosphere between 3 and 5 μm covering several absorption features for the Saturn-mass exoplanet WASP-39b (ref. 9), obtained with the Near Infrared Spectrograph (NIRSpec) G395H grating of JWST. Our observations achieve 1.46 times photon precision, providing an average transit depth uncertainty of 221 ppm per spectroscopic bin, and present minimal impacts from systematic effects. We detect significant absorption from CO2 (28.5σ) and H2O (21.5σ), and identify SO2 as the source of absorption at 4.1 μm (4.8σ). Best-fit atmospheric models range between 3 and 10 times solar metallicity, with sub-solar to solar C/O ratios. These results, including the detection of SO2, underscore the importance of characterizing the chemistry in exoplanet atmospheres and showcase NIRSpec G395H as an excellent mode for time-series observations over this critical wavelength range10.

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.

Phosphoinositide 3-kinase δ (PI3Kδ) has a key role in lymphocytes, and inhibitors that target this PI3K have been approved for treatment of B cell malignancies1-3. Although studies in mouse models of solid tumours have demonstrated that PI3Kδ inhibitors (PI3Kδi) can induce anti-tumour immunity4,5, its effect on solid tumours in humans remains unclear. Here we assessed the effects of the PI3Kδi AMG319 in human patients with head and neck cancer in a neoadjuvant, double-blind, placebo-controlled randomized phase II trial (EudraCT no. 2014-004388-20). PI3Kδ inhibition decreased the number of tumour-infiltrating regulatory T (Treg) cells and enhanced the cytotoxic potential of tumour-infiltrating T cells. At the tested doses of AMG319, immune-related adverse events (irAEs) required treatment to be discontinued in 12 out of 21 of patients treated with AMG319, suggestive of systemic effects on Treg cells. Accordingly, in mouse models, PI3Kδi decreased the number of Treg cells systemically and caused colitis. Single-cell RNA-sequencing analysis revealed a PI3Kδi-driven loss of tissue-resident colonic ST2 Treg cells, accompanied by expansion of pathogenic T helper 17 (TH17) and type 17 CD8+ T (TC17) cells, which probably contributed to toxicity; this points towards a specific mode of action for the emergence of irAEs. A modified treatment regimen with intermittent dosing of PI3Kδi in mouse models led to a significant decrease in tumour growth without inducing pathogenic T cells in colonic tissue, indicating that alternative dosing regimens might limit toxicity.

Small-ring cage hydrocarbons are popular bioisosteres (molecular replacements) for commonly found para-substituted benzene rings in drug design1. The utility of these cage structures derives from their superior pharmacokinetic properties compared with their parent aromatics, including improved solubility and reduced susceptibility to metabolism2,3. A prime example is the bicyclo[1.1.1]pentane motif, which is mainly synthesized by ring-opening of the interbridgehead bond of the strained hydrocarbon [1.1.1]propellane with radicals or anions4. By contrast, scaffolds mimicking meta-substituted arenes are lacking because of the challenge of synthesizing saturated isosteres that accurately reproduce substituent vectors5. Here we show that bicyclo[3.1.1]heptanes (BCHeps), which are hydrocarbons for which the bridgehead substituents map precisely onto the geometry of meta-substituted benzenes, can be conveniently accessed from [3.1.1]propellane. We found that [3.1.1]propellane can be synthesized on a multigram scale, and readily undergoes a range of radical-based transformations to generate medicinally relevant carbon- and heteroatom-substituted BCHeps, including pharmaceutical analogues. Comparison of the absorption, distribution, metabolism and excretion (ADME) properties of these analogues reveals enhanced metabolic stability relative to their parent arene-containing drugs, validating the potential of this meta-arene analogue as an sp3-rich motif in drug design. Collectively, our results show that BCHeps can be prepared on useful scales using a variety of methods, offering a new surrogate for meta-substituted benzene rings for implementation in drug discovery programmes.

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.

The evolution of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus leads to new variants that warrant timely epidemiological characterization. Here we use the dense genomic surveillance data generated by the COVID-19 Genomics UK Consortium to reconstruct the dynamics of 71 different lineages in each of 315 English local authorities between September 2020 and June 2021. This analysis reveals a series of subepidemics that peaked in early autumn 2020, followed by a jump in transmissibility of the B.1.1.7/Alpha lineage. The Alpha variant grew when other lineages declined during the second national lockdown and regionally tiered restrictions between November and December 2020. A third more stringent national lockdown suppressed the Alpha variant and eliminated nearly all other lineages in early 2021. Yet a series of variants (most of which contained the spike E484K mutation) defied these trends and persisted at moderately increasing proportions. However, by accounting for sustained introductions, we found that the transmissibility of these variants is unlikely to have exceeded the transmissibility of the Alpha variant. Finally, B.1.617.2/Delta was repeatedly introduced in England and grew rapidly in early summer 2021, constituting approximately 98% of sampled SARS-CoV-2 genomes on 26 June 2021.

Bacteria encode myriad defences that target the genomes of infecting bacteriophage, including restriction-modification and CRISPR-Cas systems1. In response, one family of large bacteriophages uses a nucleus-like compartment to protect its replicating genomes by excluding host defence factors2-4. However, the principal composition and structure of this compartment remain unknown. Here we find that the bacteriophage nuclear shell assembles primarily from one protein, which we name chimallin (ChmA). Combining cryo-electron tomography of nuclear shells in bacteriophage-infected cells and cryo-electron microscopy of a minimal chimallin compartment in vitro, we show that chimallin self-assembles as a flexible sheet into closed micrometre-scale compartments. The architecture and assembly dynamics of the chimallin shell suggest mechanisms for its nucleation and growth, and its role as a scaffold for phage-encoded factors mediating macromolecular transport, cytoskeletal interactions, and viral maturation.

The grey wolf (Canis lupus) was the first species to give rise to a domestic population, and they remained widespread throughout the last Ice Age when many other large mammal species went extinct. Little is known, however, about the history and possible extinction of past wolf populations or when and where the wolf progenitors of the present-day dog lineage (Canis familiaris) lived1-8. Here we analysed 72 ancient wolf genomes spanning the last 100,000 years from Europe, Siberia and North America. We found that wolf populations were highly connected throughout the Late Pleistocene, with levels of differentiation an order of magnitude lower than they are today. This population connectivity allowed us to detect natural selection across the time series, including rapid fixation of mutations in the gene IFT88 40,000-30,000 years ago. We show that dogs are overall more closely related to ancient wolves from eastern Eurasia than to those from western Eurasia, suggesting a domestication process in the east. However, we also found that dogs in the Near East and Africa derive up to half of their ancestry from a distinct population related to modern southwest Eurasian wolves, reflecting either an independent domestication process or admixture from local wolves. None of the analysed ancient wolf genomes is a direct match for either of these dog ancestries, meaning that the exact progenitor populations remain to be located.

Translation initiation defines the identity and quantity of a synthesized protein. The process is dysregulated in many human diseases1,2. A key commitment step is when the ribosomal subunits join at a translation start site on a messenger RNA to form a functional ribosome. Here, we combined single-molecule spectroscopy and structural methods using an in vitro reconstituted system to examine how the human ribosomal subunits join. Single-molecule fluorescence revealed when the universally conserved eukaryotic initiation factors eIF1A and eIF5B associate with and depart from initiation complexes. Guided by single-molecule dynamics, we visualized initiation complexes that contained both eIF1A and eIF5B using single-particle cryo-electron microscopy. The resulting structure revealed how eukaryote-specific contacts between the two proteins remodel the initiation complex to orient the initiator aminoacyl-tRNA in a conformation compatible with ribosomal subunit joining. Collectively, our findings provide a quantitative and architectural framework for the molecular choreography orchestrated by eIF1A and eIF5B during translation initiation in humans.

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 this Article, the authors 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) were mistakenly omitted from the ENCODE Project Consortium author list. The original Article has been corrected online.

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.

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.

The land ice contribution to global mean sea level rise has not yet been predicted1 using ice sheet and glacier models for the latest set of socio-economic scenarios, nor using coordinated exploration of uncertainties arising from the various computer models involved. Two recent international projects generated a large suite of projections using multiple models2-8, but primarily used previous-generation scenarios9 and climate models10, and could not fully explore known uncertainties. Here we estimate probability distributions for these projections under the new scenarios11,12 using statistical emulation of the ice sheet and glacier models. We find that limiting global warming to 1.5 degrees Celsius would halve the land ice contribution to twenty-first-century sea level rise, relative to current emissions pledges. The median decreases from 25 to 13 centimetres sea level equivalent (SLE) by 2100, with glaciers responsible for half the sea level contribution. The projected Antarctic contribution does not show a clear response to the emissions scenario, owing to uncertainties in the competing processes of increasing ice loss and snowfall accumulation in a warming climate. However, under risk-averse (pessimistic) assumptions, Antarctic ice loss could be five times higher, increasing the median land ice contribution to 42 centimetres SLE under current policies and pledges, with the 95th percentile projection exceeding half a metre even under 1.5 degrees Celsius warming. This would severely limit the possibility of mitigating future coastal flooding. Given this large range (between 13 centimetres SLE using the main projections under 1.5 degrees Celsius warming and 42 centimetres SLE using risk-averse projections under current pledges), adaptation planning for twenty-first-century sea level rise must account for a factor-of-three uncertainty in the land ice contribution until climate policies and the Antarctic response are further constrained.

In Supplementary Table 1 of this Article, 23 samples (B10K-DU-029-32, B10K-DU-029-33, B10K-DU-029-36 to B10K-DU-029-44, B10K-DU- 029-46, B10K-DU-029-47, B10K-DU-029-49 to B10K-DU-029-53, B10K-DU- 029-75 to B10K-DU-029-77, B10K-DU-029-80, and B10K-DU-030-03; styled in boldface in the revised table) were assigned to the incorrect institution. Supplementary Table 1 has been amended to reflect the correct source institution for these samples, and associated data (tissue, museum ID/source specimen ID, site, state/province, latitude, longitude, date collected and sex) have been updated accordingly. The original table is provided as Supplementary Information to this Amendment, and the original Article has been corrected online.

Differentiation proceeds along a continuum of increasingly fate-restricted intermediates, referred to as canalization1,2. Canalization is essential for stabilizing cell fate, but the mechanisms that underlie robust canalization are unclear. Here we show that the BRG1/BRM-associated factor (BAF) chromatin-remodelling complex ATPase gene Brm safeguards cell identity during directed cardiogenesis of mouse embryonic stem cells. Despite the establishment of a well-differentiated precardiac mesoderm, Brm-/- cells predominantly became neural precursors, violating germ layer assignment. Trajectory inference showed a sudden acquisition of a non-mesodermal identity in Brm-/- cells. Mechanistically, the loss of Brm prevented de novo accessibility of primed cardiac enhancers while increasing the expression of neurogenic factor POU3F1, preventing the binding of the neural suppressor REST and shifting the composition of BRG1 complexes. The identity switch caused by the Brm mutation was overcome by increasing BMP4 levels during mesoderm induction. Mathematical modelling supports these observations and demonstrates that Brm deletion affects cell fate trajectory by modifying saddle-node bifurcations2. In the mouse embryo, Brm deletion exacerbated mesoderm-deleted Brg1-mutant phenotypes, severely compromising cardiogenesis, and reveals an in vivo role for Brm. Our results show that Brm is a compensable safeguard of the fidelity of mesoderm chromatin states, and support a model in which developmental canalization is not a rigid irreversible path, but a highly plastic trajectory.

Quantum many-body systems display rich phase structure in their low-temperature equilibrium states1. However, much of nature is not in thermal equilibrium. Remarkably, it was recently predicted that out-of-equilibrium systems can exhibit novel dynamical phases2-8 that may otherwise be forbidden by equilibrium thermodynamics, a paradigmatic example being the discrete time crystal (DTC)7,9-15. Concretely, dynamical phases can be defined in periodically driven many-body-localized (MBL) systems via the concept of eigenstate order7,16,17. In eigenstate-ordered MBL phases, the entire many-body spectrum exhibits quantum correlations and long-range order, with characteristic signatures in late-time dynamics from all initial states. It is, however, challenging to experimentally distinguish such stable phases from transient phenomena, or from regimes in which the dynamics of a few select states can mask typical behaviour. Here we implement tunable controlled-phase (CPHASE) gates on an array of superconducting qubits to experimentally observe an MBL-DTC and demonstrate its characteristic spatiotemporal response for generic initial states7,9,10. Our work employs a time-reversal protocol to quantify the impact of external decoherence, and leverages quantum typicality to circumvent the exponential cost of densely sampling the eigenspectrum. Furthermore, we locate the phase transition out of the DTC with an experimental finite-size analysis. These results establish a scalable approach to studying non-equilibrium phases of matter on quantum processors.

Increased blood lipid levels are heritable risk factors of cardiovascular disease with varied prevalence worldwide owing to different dietary patterns and medication use1. Despite advances in prevention and treatment, in particular through reducing low-density lipoprotein cholesterol levels2, heart disease remains the leading cause of death worldwide3. Genome-wideassociation studies (GWAS) of blood lipid levels have led to important biological and clinical insights, as well as new drug targets, for cardiovascular disease. However, most previous GWAS4-23 have been conducted in European ancestry populations and may have missed genetic variants that contribute to lipid-level variation in other ancestry groups. These include differences in allele frequencies, effect sizes and linkage-disequilibrium patterns24. Here we conduct a multi-ancestry, genome-wide genetic discovery meta-analysis of lipid levels in approximately 1.65 million individuals, including 350,000 of non-European ancestries. We quantify the gain in studying non-European ancestries and provide evidence to support the expansion of recruitment of additional ancestries, even with relatively small sample sizes. We find that increasing diversity rather than studying additional individuals of European ancestry results in substantial improvements in fine-mapping functional variants and portability of polygenic prediction (evaluated in approximately 295,000 individuals from 7 ancestry groupings). Modest gains in the number of discovered loci and ancestry-specific variants were also achieved. As GWAS expand emphasis beyond the identification of genes and fundamental biology towards the use of genetic variants for preventive and precision medicine25, we anticipate that increased diversity of participants will lead to more accurate and equitable26 application of polygenic scores in clinical practice.

Physical exercise is generally beneficial to all aspects of human and animal health, slowing cognitive ageing and neurodegeneration1. The cognitive benefits of physical exercise are tied to an increased plasticity and reduced inflammation within the hippocampus2-4, yet little is known about the factors and mechanisms that mediate these effects. Here we show that 'runner plasma', collected from voluntarily running mice and infused into sedentary mice, reduces baseline neuroinflammatory gene expression and experimentally induced brain inflammation. Plasma proteomic analysis revealed a concerted increase in complement cascade inhibitors including clusterin (CLU). Intravenously injected CLU binds to brain endothelial cells and reduces neuroinflammatory gene expression in a mouse model of acute brain inflammation and a mouse model of Alzheimer's disease. Patients with cognitive impairment who participated in structured exercise for 6 months had higher plasma levels of CLU. These findings demonstrate the existence of anti-inflammatory exercise factors that are transferrable, target the cerebrovasculature and benefit the brain, and are present in humans who engage in exercise.

Long-term climate change and periodic environmental extremes threaten food and fuel security1 and global crop productivity2-4. Although molecular and adaptive breeding strategies can buffer the effects of climatic stress and improve crop resilience5, these approaches require sufficient knowledge of the genes that underlie productivity and adaptation6-knowledge that has been limited to a small number of well-studied model systems. Here we present the assembly and annotation of the large and complex genome of the polyploid bioenergy crop switchgrass (Panicum virgatum). Analysis of biomass and survival among 732 resequenced genotypes, which were grown across 10 common gardens that span 1,800 km of latitude, jointly revealed extensive genomic evidence of climate adaptation. Climate-gene-biomass associations were abundant but varied considerably among deeply diverged gene pools. Furthermore, we found that gene flow accelerated climate adaptation during the postglacial colonization of northern habitats through introgression of alleles from a pre-adapted northern gene pool. The polyploid nature of switchgrass also enhanced adaptive potential through the fractionation of gene function, as there was an increased level of heritable genetic diversity on the nondominant subgenome. In addition to investigating patterns of climate adaptation, the genome resources and gene-trait associations developed here provide breeders with the necessary tools to increase switchgrass yield for the sustainable production of bioenergy.

Whole-genome sequencing projects are increasingly populating the tree of life and characterizing biodiversity1-4. Sparse taxon sampling has previously been proposed to confound phylogenetic inference5, and captures only a fraction of the genomic diversity. Here we report a substantial step towards the dense representation of avian phylogenetic and molecular diversity, by analysing 363 genomes from 92.4% of bird families-including 267 newly sequenced genomes produced for phase II of the Bird 10,000 Genomes (B10K) Project. We use this comparative genome dataset in combination with a pipeline that leverages a reference-free whole-genome alignment to identify orthologous regions in greater numbers than has previously been possible and to recognize genomic novelties in particular bird lineages. The densely sampled alignment provides a single-base-pair map of selection, has more than doubled the fraction of bases that are confidently predicted to be under conservation and reveals extensive patterns of weak selection in predominantly non-coding DNA. Our results demonstrate that increasing the diversity of genomes used in comparative studies can reveal more shared and lineage-specific variation, and improve the investigation of genomic characteristics. We anticipate that this genomic resource will offer new perspectives on evolutionary processes in cross-species comparative analyses and assist in efforts to conserve species.

The maritime expansion of Scandinavian populations during the Viking Age (about AD 750-1050) was a far-flung transformation in world history1,2. Here we sequenced the genomes of 442 humans from archaeological sites across Europe and Greenland (to a median depth of about 1×) to understand the global influence of this expansion. We find the Viking period involved gene flow into Scandinavia from the south and east. We observe genetic structure within Scandinavia, with diversity hotspots in the south and restricted gene flow within Scandinavia. We find evidence for a major influx of Danish ancestry into England; a Swedish influx into the Baltic; and Norwegian influx into Ireland, Iceland and Greenland. Additionally, we see substantial ancestry from elsewhere in Europe entering Scandinavia during the Viking Age. Our ancient DNA analysis also revealed that a Viking expedition included close family members. By comparing with modern populations, we find that pigmentation-associated loci have undergone strong population differentiation during the past millennium, and trace positively selected loci-including the lactase-persistence allele of LCT and alleles of ANKA that are associated with the immune response-in detail. We conclude that the Viking diaspora was characterized by substantial transregional engagement: distinct populations influenced the genomic makeup of different regions of Europe, and Scandinavia experienced increased contact with the rest of the continent.

We are a group of archaeologists, anthropologists, curators and geneticists representing diverse global communities and 31 countries. All of us met in a virtual workshop dedicated to ethics in ancient DNA research held in November 2020. There was widespread agreement that globally applicable ethical guidelines are needed, but that recent recommendations grounded in discussion about research on human remains from North America are not always generalizable worldwide. Here we propose the following globally applicable guidelines, taking into consideration diverse contexts. These hold that: (1) researchers must ensure that all regulations were followed in the places where they work and from which the human remains derived; (2) researchers must prepare a detailed plan prior to beginning any study; (3) researchers must minimize damage to human remains; (4) researchers must ensure that data are made available following publication to allow critical re-examination of scientific findings; and (5) researchers must engage with other stakeholders from the beginning of a study and ensure respect and sensitivity to stakeholder perspectives. We commit to adhering to these guidelines and expect they will promote a high ethical standard in DNA research on human remains going forward.

The primary motor cortex (M1) is essential for voluntary fine-motor control and is functionally conserved across mammals1. Here, using high-throughput transcriptomic and epigenomic profiling of more than 450,000 single nuclei in humans, marmoset monkeys and mice, we demonstrate a broadly conserved cellular makeup of this region, with similarities that mirror evolutionary distance and are consistent between the transcriptome and epigenome. The core conserved molecular identities of neuronal and non-neuronal cell types allow us to generate a cross-species consensus classification of cell types, and to infer conserved properties of cell types across species. Despite the overall conservation, however, many species-dependent specializations are apparent, including differences in cell-type proportions, gene expression, DNA methylation and chromatin state. Few cell-type marker genes are conserved across species, revealing a short list of candidate genes and regulatory mechanisms that are responsible for conserved features of homologous cell types, such as the GABAergic chandelier cells. This consensus transcriptomic classification allows us to use patch-seq (a combination of whole-cell patch-clamp recordings, RNA sequencing and morphological characterization) to identify corticospinal Betz cells from layer 5 in non-human primates and humans, and to characterize their highly specialized physiology and anatomy. These findings highlight the robust molecular underpinnings of cell-type diversity in M1 across mammals, and point to the genes and regulatory pathways responsible for the functional identity of cell types and their species-specific adaptations.

Loss of skeletal integrity during ageing and disease is associated with an imbalance in the opposing actions of osteoblasts and osteoclasts1. Here we show that intrinsic ageing of skeletal stem cells (SSCs)2 in mice alters signalling in the bone marrow niche and skews the differentiation of bone and blood lineages, leading to fragile bones that regenerate poorly. Functionally, aged SSCs have a decreased bone- and cartilage-forming potential but produce more stromal lineages that express high levels of pro-inflammatory and pro-resorptive cytokines. Single-cell RNA-sequencing studies link the functional loss to a diminished transcriptomic diversity of SSCs in aged mice, which thereby contributes to the transformation of the bone marrow niche. Exposure to a youthful circulation through heterochronic parabiosis or systemic reconstitution with young haematopoietic stem cells did not reverse the diminished osteochondrogenic activity of aged SSCs, or improve bone mass or skeletal healing parameters in aged mice. Conversely, the aged SSC lineage promoted osteoclastic activity and myeloid skewing by haematopoietic stem and progenitor cells, suggesting that the ageing of SSCs is a driver of haematopoietic ageing. Deficient bone regeneration in aged mice could only be returned to youthful levels by applying a combinatorial treatment of BMP2 and a CSF1 antagonist locally to fractures, which reactivated aged SSCs and simultaneously ablated the inflammatory, pro-osteoclastic milieu. Our findings provide mechanistic insights into the complex, multifactorial mechanisms that underlie skeletal ageing and offer prospects for rejuvenating the aged skeletal system.