The 100 most-cited articles in Nature, 2022

SUMMARY Schizophrenia is a psychiatric disorder whose pathophysiology is largely unknown. It has a heritability of 60-80%, much of which is attributable to common risk alleles, suggesting genome-wide association studies can inform our understanding of aetiology 1 . Here, in 69,369 people with schizophrenia and 236,642 controls, we report common variant associations at 270 distinct loci. Using fine-mapping and functional genomic data, we prioritise 19 genes based on protein-coding or UTR variation, and 130 genes in total as likely to explain these associations. 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 autism and developmental disorder. Associations were concentrated in genes expressed in CNS neurons, both excitatory and inhibitory, but not other tissues or cell types, and implicated fundamental processes related to neuronal function, particularly synaptic…

Abstract Magnetic resonance imaging (MRI) has transformed our understanding of the human brain through well-replicated mapping of abilities to specific structures (for example, lesion studies) and functions 1–3 (for example, task functional MRI (fMRI)). Mental health research and care have yet to realize similar advances from MRI. A primary challenge has been replicating associations between inter-individual differences in brain structure or function and complex cognitive or mental health phenotypes (brain-wide association studies (BWAS)). Such BWAS have typically relied on sample sizes appropriate for classical brain mapping 4 (the median neuroimaging study sample size is about 25), but potentially too small for capturing reproducible brain–behavioural phenotype associations 5,6 . Here we used three of the largest neuroimaging datasets currently available—with a total sample size of around 50,000 individuals—to quantify BWAS effect sizes and reproducibility as a function of sample size. BWAS associations were smaller than previously thought, resulting in statistically underpowered…

Abstract 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 weight 1 . 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 trajectories 2 of brain structural changes, and rates of change, over the lifespan. Brain charts identified previously unreported neurodevelopmental…

Abstract The SARS-CoV-2 epidemic in southern Africa has been characterized by three distinct waves. The first was associated with a mix of SARS-CoV-2 lineages, while the second and third waves were driven by the Beta (B.1.351) and Delta (B.1.617.2) variants, respectively 1–3 . In November 2021, genomic surveillance teams in South Africa and Botswana detected a new SARS-CoV-2 variant associated with a rapid resurgence of infections in Gauteng province, South Africa. Within three days of the first genome being uploaded, it was designated a variant of concern (Omicron, B.1.1.529) by the World Health Organization and, within three weeks, had been identified in 87 countries. The Omicron variant is exceptional for carrying over 30 mutations in the spike glycoprotein, which are predicted to influence antibody neutralization and spike function 4 . Here we describe the genomic profile and early transmission dynamics of Omicron, highlighting the rapid spread in regions with high levels of population immunity.

. However, it remains unknown whether the impact of SARS-CoV-2 infection can be detected in milder cases, and whether this can reveal possible mechanisms contributing to brain pathology. Here we investigated brain changes in 785 participants of UK Biobank (aged 51-81 years) who were imaged twice using magnetic resonance imaging, including 401 cases who tested positive for infection with SARS-CoV-2 between their two scans-with 141 days on average separating their diagnosis and the second scan-as well as 384 controls. The availability of pre-infection imaging data reduces the likelihood of pre-existing risk factors being misinterpreted as disease effects. We identified significant longitudinal effects when comparing the two groups, including (1) a greater reduction in grey matter thickness and tissue contrast in the orbitofrontal cortex and parahippocampal gyrus; (2) greater changes in markers of tissue damage in regions that are functionally connected to the primary olfactory cortex; and (3) a greater reduction in global brain size in the SARS-CoV-2 cases. The participants who were…

Abstract Recent emergence of SARS-CoV-2 Omicron sublineages BA.2.12.1, BA.2.13, BA.4 and BA.5 all contain L452 mutations and show potential higher transmissibility over BA.2. The new variants’ receptor binding and immune evasion capability require immediate investigation, especially on the role of L452 substitutions. Herein, coupled with structural comparisons, we showed that BA.2 sublineages, including BA.2.12.1 and BA.2.13, exhibit increased ACE2-binding affinities compared to BA.1; while BA.4/BA.5 shows the weakest receptor-binding activity due to F486V and R493Q reversion. Importantly, compared to BA.2, BA.2.12.1 and BA.4/BA.5 exhibit stronger neutralization escape from the plasma of 3-dose vaccinees and, most strikingly, from vaccinated BA.1 convalescents. To delineate the underlying evasion mechanism, we determined the escaping mutation profiles, epitope distribution and Omicron sub-lineage neutralization efficacy of 1640 RBD-directed neutralizing antibodies (NAbs), including 614 isolated from BA.1 convalescents. Interestingly, post-vaccination BA.1 infection mainly recalls…

The SARS-CoV-2 Omicron BA.1 variant emerged in late 2021 and is characterised by multiple spike mutations across all spike domains. Here we show that Omicron BA.1 has higher affinity for ACE2 compared to Delta, and confers very significant evasion of therapeutic monoclonal and vaccine-elicited polyclonal neutralising antibodies after two doses. mRNA vaccination as a third vaccine dose rescues and broadens neutralisation. Importantly, antiviral drugs remdesevir and molnupiravir retain efficacy against Omicron BA.1. We found that in human nasal epithelial 3D cultures replication was similar for both Omicron and Delta. However, in lower airway organoids, Calu-3 lung cells and gut adenocarcinoma cell lines live Omicron virus demonstrated significantly lower replication in comparison to Delta. We noted that despite presence of mutations predicted to favour spike S1/S2 cleavage, the spike protein is less efficiently cleaved in live Omicron virions compared to Delta virions. We mapped the replication differences between the variants to entry efficiency using spike pseudotyped virus (PV)…

Abstract At least 10,000 species of mammal virus are estimated to have the potential to spread in human populations, but the vast majority are currently circulating in wildlife, largely undescribed and undetected by disease outbreak surveillance 1–3 . In addition, changing climate and land use are already driving geographic range shifts in wildlife, producing novel species assemblages and opportunities for viral sharing between previously isolated species 4, 5 . In some cases, this will inevitably facilitate spillover into humans 6, 7 —a possible mechanistic link between global environmental change and emerging zoonotic disease 8 . Here, we map potential hotspots of viral sharing, using a phylogeographic model of the mammal-virus network, and projections of potential geographic range shifts for 3,139 mammal species under climate change and land use scenarios for the year 2070. Range-shifting mammal species are predicted to aggregate at high elevations, in biodiversity hotspots, and in areas of high human population density in Asia and Africa, driving the novel cross-species…

Abstract The social cost of carbon dioxide (SC-CO 2 ) measures the monetized value of the damages to society caused by an incremental metric tonne of CO 2 emissions and is a key metric informing climate policy. Used by governments and other decision-makers in benefit–cost analysis for over a decade, SC-CO 2 estimates draw on climate science, economics, demography and other disciplines. However, a 2017 report by the US National Academies of Sciences, Engineering, and Medicine 1 (NASEM) highlighted that current SC-CO 2 estimates no longer reflect the latest research. The report provided a series of recommendations for improving the scientific basis, transparency and uncertainty characterization of SC-CO 2 estimates. Here we show that improved probabilistic socioeconomic projections, climate models, damage functions, and discounting methods that collectively reflect theoretically consistent valuation of risk, substantially increase estimates of the SC-CO 2 . Our preferred mean SC-CO 2 estimate is $185 per tonne of CO 2 ($44–$413 per tCO 2 : 5%–95% range, 2020 US dollars) at a near-term…

Abstract As a part of its mandate to compile and disseminate statistics on mortality, the World Health Organization (WHO) has been tracking the progression of the COVID-19 pandemic since the beginning of 2020. However, reported statistics on COVID-19 are problematic for a number of countries due to variations in testing access, differential diagnostic capacity and inconsistencies in the applications of standards to correctly certify COVID-19 as cause-of-death. In addition, the pandemic has caused extensive collateral damage beyond what is directly attributable to it. Consequently, the WHO has estimated excess deaths for each country for the years 2020 and 2021 to quantify the pandemic’s impact more comprehensively and consistently. Defined as the number of deaths in a particular period, relative to that expected during “normal times”, excess deaths capture both the direct and indirect impacts of a crisis. The data required to estimate the excess mortality associated with the COVID-19 pandemic i.e., time-series of known deaths during the pandemic period and historical time-series of…

Abstract A quantum computer attains computational advantage when outperforming the best classical computers running the best-known algorithms on well-defined tasks. No photonic machine offering programmability over all its quantum gates has demonstrated quantum computational advantage: previous machines 1,2 were largely restricted to static gate sequences. Earlier photonic demonstrations were also vulnerable to spoofing 3 , in which classical heuristics produce samples, without direct simulation, lying closer to the ideal distribution than do samples from the quantum hardware. Here we report quantum computational advantage using Borealis, a photonic processor offering dynamic programmability on all gates implemented. We carry out Gaussian boson sampling 4 (GBS) on 216 squeezed modes entangled with three-dimensional connectivity 5 , using a time-multiplexed and photon-number-resolving architecture. On average, it would take more than 9,000 years for the best available algorithms and supercomputers to produce, using exact methods, a single sample from the programmed distribution,…

Abstract By meta-analyzing the whole-exomes of 24,248 cases and 97,322 controls, we implicate ultra-rare coding variants (URVs) in ten genes as conferring substantial risk for schizophrenia (odds ratios 3 - 50, P < 2.14 × 10 -6 ), and 32 genes at a FDR < 5%. These genes have the greatest expression in central nervous system neurons and have diverse molecular functions that include the formation, structure, and function of the synapse. The associations of NMDA receptor subunit GRIN2A and AMPA receptor subunit GRIA3 provide support for the dysfunction of the glutamatergic system as a mechanistic hypothesis in the pathogenesis of schizophrenia. We find significant evidence for an overlap of rare variant risk between schizophrenia, autism spectrum disorders (ASD), and severe neurodevelopmental disorders (DD/ID), supporting a neurodevelopmental etiology for schizophrenia. We show that protein-truncating variants in GRIN2A, TRIO , and CACNA1G confer risk for schizophrenia whereas specific missense mutations in these genes confer risk for DD/ID. Nevertheless, few of the strongly associated…

Abstract Realizing increasingly complex artificial intelligence (AI) functionalities directly on edge devices calls for unprecedented energy efficiency of edge hardware. Compute-in-memory (CIM) based on resistive random-access memory (RRAM) 1 promises to meet such demand by storing AI model weights in dense, analogue and non-volatile RRAM devices, and by performing AI computation directly within RRAM, thus eliminating power-hungry data movement between separate compute and memory 2–5 . Although recent studies have demonstrated in-memory matrix-vector multiplication on fully integrated RRAM-CIM hardware 6–17 , it remains a goal for a RRAM-CIM chip to simultaneously deliver high energy efficiency, versatility to support diverse models and software-comparable accuracy. Although efficiency, versatility and accuracy are all indispensable for broad adoption of the technology, the inter-related trade-offs among them cannot be addressed by isolated improvements on any single abstraction level of the design. Here, by co-optimizing across all hierarchies of the design from algorithms and…

Abstract The tumour-associated microbiota is an intrinsic component of the tumour microenvironment across human cancer types 1,2 . Intratumoral host–microbiota studies have so far largely relied on bulk tissue analysis 1–3 , which obscures the spatial distribution and localized effect of the microbiota within tumours. Here, by applying in situ spatial-profiling technologies 4 and single-cell RNA sequencing 5 to oral squamous cell carcinoma and colorectal cancer, we reveal spatial, cellular and molecular host–microbe interactions. We adapted 10x Visium spatial transcriptomics to determine the identity and in situ location of intratumoral microbial communities within patient tissues. Using GeoMx digital spatial profiling 6 , we show that bacterial communities populate microniches that are less vascularized, highly immuno‑suppressive and associated with malignant cells with lower levels of Ki-67 as compared to bacteria-negative tumour regions. We developed a single-cell RNA-sequencing method that we name INVADEseq (invasion–adhesion-directed expression sequencing) and, by applying this…

Abstract Continuous evolution of Omicron has led to a rapid and simultaneous emergence of numerous variants that display growth advantages over BA. 5. Despite their divergent evolutionary courses, mutations on their receptor-binding domain (RBD) converge on several hotspots. The driving force and destination of such convergent evolution and its impact on humoral immunity remain unclear. Here, we demonstrate that these convergent mutations can cause striking evasion of neutralizing antibody (NAb) drugs and convalescent plasma, including those from BA.5 breakthrough infection, while maintaining sufficient ACE2 binding capability. BQ.1.1.10, BA.4.6.3, XBB, and CH. 1.1 are the most antibody-evasive strain tested, even exceeding SARS-CoV-1 level. To delineate the origin of the convergent evolution, we determined the escape mutation profiles and neutralization activity of monoclonal antibodies (mAbs) isolated from BA.2 and BA.5 breakthrough-infection convalescents. Importantly, due to humoral immune imprinting, BA.2 and especially BA.5 breakthrough infection caused significant reductions…

Abstract Diffuse intrinsic pontine glioma (DIPG) and other H3K27M-mutated diffuse midline gliomas (DMGs) are universally lethal paediatric tumours of the central nervous system 1 . We have previously shown that the disialoganglioside GD2 is highly expressed on H3K27M-mutated glioma cells and have demonstrated promising preclinical efficacy of GD2-directed chimeric antigen receptor (CAR) T cells 2 , providing the rationale for a first-in-human phase I clinical trial (NCT04196413). Because CAR T cell-induced brainstem inflammation can result in obstructive hydrocephalus, increased intracranial pressure and dangerous tissue shifts, neurocritical care precautions were incorporated. Here we present the clinical experience from the first four patients with H3K27M-mutated DIPG or spinal cord DMG treated with GD2-CAR T cells at dose level 1 (1 × 10 6 GD2-CAR T cells per kg administered intravenously). Patients who exhibited clinical benefit were eligible for subsequent GD2-CAR T cell infusions administered intracerebroventricularly 3 . Toxicity was largely related to the location of the…

Abstract 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 sizes 1 . 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 panel 2 ) 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…

Abstract The identification of the Omicron (B.1.1.529.1 or BA.1) variant of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in Botswana in November 2021 1 immediately caused concern owing to the number of alterations in the spike glycoprotein that could lead to antibody evasion. We 2 and others 3–6 recently reported results confirming such a concern. Continuing surveillance of the evolution of Omicron has since revealed the rise in prevalence of two sublineages, BA.1 with an R346K alteration (BA.1+R346K, also known as BA.1.1) and B.1.1.529.2 (BA.2), with the latter containing 8 unique spike alterations and lacking 13 spike alterations found in BA.1. Here we extended our studies to include antigenic characterization of these new sublineages. Polyclonal sera from patients infected by wild-type SARS-CoV-2 or recipients of current mRNA vaccines showed a substantial loss in neutralizing activity against both BA.1+R346K and BA.2, with drops comparable to that already reported for BA.1 (refs. 2,3,5,6 ). These findings indicate that these three sublineages of Omicron are…

Abstract SARS-CoV-2 Omicron subvariants BA.2.12.1 and BA.4/5 have surged dramatically to become dominant in the United States and South Africa, respectively 1,2 . These novel subvariants carrying additional mutations in their spike proteins raise concerns that they may further evade neutralizing antibodies, thereby further compromising the efficacy of COVID-19 vaccines and therapeutic monoclonals. We now report findings from a systematic antigenic analysis of these surging Omicron subvariants. BA.2.12.1 is only modestly (1.8-fold) more resistant to sera from vaccinated and boosted individuals than BA.2. However, BA.4/5 is substantially (4.2-fold) more resistant and thus more likely to lead to vaccine breakthrough infections. Mutation at spike residue L452 found in both BA.2.12.1 and BA.4/5 facilitates escape from some antibodies directed to the so-called class 2 and 3 regions of the receptor-binding domain 3 . The F486V mutation found in BA.4/5 facilitates escape from certain class 1 and 2 antibodies but compromises the spike affinity for the viral receptor. The R493Q reversion…

Abstract Forest ecosystems depend on their capacity to withstand and recover from natural and anthropogenic perturbations (that is, their resilience) 1 . Experimental evidence of sudden increases in tree mortality is raising concerns about variation in forest resilience 2 , yet little is known about how it is evolving in response to climate change. Here we integrate satellite-based vegetation indices with machine learning to show how forest resilience, quantified in terms of critical slowing down indicators 3–5 , has changed during the period 2000–2020. We show that tropical, arid and temperate forests are experiencing a significant decline in resilience, probably related to increased water limitations and climate variability. By contrast, boreal forests show divergent local patterns with an average increasing trend in resilience, probably benefiting from warming and CO 2 fertilization, which may outweigh the adverse effects of climate change. These patterns emerge consistently in both managed and intact forests, corroborating the existence of common large-scale climate drivers.…

Abstract Myocardial infarction is a leading cause of mortality. While advances in the acute treatment have been made, the late-stage mortality is still high, driven by an incomplete understanding of cardiac remodeling processes 1,2 . Here we used single-cell gene expression, chromatin accessibility and spatial transcriptomic profiling of different physiological zones and timepoints of human myocardial infarction and human control myocardium to generate an integrative high-resolution map of cardiac remodeling. This approach allowed us to increase spatial resolution of cell-type composition and provide spatially resolved insights into the cardiac transcriptome and epigenome with identification of distinct cellular zones of injury, repair and remodeling. We here identified and validated mechanisms of fibroblast to myofibroblast differentiation that drive cardiac fibrosis. Our study provides an integrative molecular map of human myocardial infarction and represents a reference to advance mechanistic and therapeutic studies of cardiac disease.

Abstract Nuclear fusion using magnetic confinement, in particular in the tokamak configuration, is a promising path towards sustainable energy. A core challenge is to shape and maintain a high-temperature plasma within the tokamak vessel. This requires high-dimensional, high-frequency, closed-loop control using magnetic actuator coils, further complicated by the diverse requirements across a wide range of plasma configurations. In this work, we introduce a previously undescribed architecture for tokamak magnetic controller design that autonomously learns to command the full set of control coils. This architecture meets control objectives specified at a high level, at the same time satisfying physical and operational constraints. This approach has unprecedented flexibility and generality in problem specification and yields a notable reduction in design effort to produce new plasma configurations. We successfully produce and control a diverse set of plasma configurations on the Tokamak à Configuration Variable 1,2 , including elongated, conventional shapes, as well as advanced…

Abstract Deep-learning models have become pervasive tools in science and engineering. However, their energy requirements now increasingly limit their scalability 1 . Deep-learning accelerators 2–9 aim to perform deep learning energy-efficiently, usually targeting the inference phase and often by exploiting physical substrates beyond conventional electronics. Approaches so far 10–22 have been unable to apply the backpropagation algorithm to train unconventional novel hardware in situ. The advantages of backpropagation have made it the de facto training method for large-scale neural networks, so this deficiency constitutes a major impediment. Here we introduce a hybrid in situ–in silico algorithm, called physics-aware training, that applies backpropagation to train controllable physical systems. Just as deep learning realizes computations with deep neural networks made from layers of mathematical functions, our approach allows us to train deep physical neural networks made from layers of controllable physical systems, even when the physical layers lack any mathematical isomorphism to…

Abstract Ferroptosis, a non-apoptotic form of cell death marked by iron-dependent lipid peroxidation 1 , has a key role in organ injury, degenerative disease and vulnerability of therapy-resistant cancers 2 . Although substantial progress has been made in understanding the molecular processes relevant to ferroptosis, additional cell-extrinsic and cell-intrinsic processes that determine cell sensitivity toward ferroptosis remain unknown. Here we show that the fully reduced forms of vitamin K—a group of naphthoquinones that includes menaquinone and phylloquinone 3 —confer a strong anti-ferroptotic function, in addition to the conventional function linked to blood clotting by acting as a cofactor for γ-glutamyl carboxylase. Ferroptosis suppressor protein 1 (FSP1), a NAD(P)H-ubiquinone reductase and the second mainstay of ferroptosis control after glutathione peroxidase-4 4,5 , was found to efficiently reduce vitamin K to its hydroquinone, a potent radical-trapping antioxidant and inhibitor of (phospho)lipid peroxidation. The FSP1-mediated reduction of vitamin K was also responsible for…

Despite the development and deployment of antibody and vaccine countermeasures, rapidly-spreading SARS-CoV-2 variants with mutations at key antigenic sites in the spike protein jeopardize their efficacy. The recent emergence of B.1.1.529, the Omicron variant1,2, which has more than 30 mutations in the spike protein, has raised concerns for escape from protection by vaccines and therapeutic antibodies. A key test for potential countermeasures against B.1.1.529 is their activity in pre-clinical rodent models of respiratory tract disease. Here, using the collaborative network of the SARS-CoV-2 Assessment of Viral Evolution (SAVE) program of the National Institute of Allergy and Infectious Diseases (NIAID), we evaluated the ability of multiple B.1.1.529 Omicron isolates to cause infection and disease in immunocompetent and human ACE2 (hACE2) expressing mice and hamsters. Despite modeling and binding data suggesting that B.1.1.529 spike can bind more avidly to murine ACE2, we observed attenuation of infection in 129, C57BL/6, and BALB/c mice as compared with previous SARS-CoV-2 variants,…

Abstract The emergence of a new severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variant, Omicron, is the most urgent concern in the global health in December 2021. Our statistical modelling estimates that Omicron is >3.0-fold and >5.6-fold more transmissible than Delta in South Africa and the UK, respectively. Intriguingly, cell culture experiments show that Omicron is less fusogenic than Delta and ancestral SARS-CoV-2. Although the spike (S) protein of Delta is efficiently cleaved into the two subunits, which facilitates cell-cell fusion, Omicron S is faintly cleaved. Further, in hamster model, Omicron shows decreased lung infectivity and is less pathogenic compared to Delta and ancestral SARS-CoV-2. Our data suggest that the efficacy of SARS-CoV-2 S cleavage and viral fusogenicity are closely associated with viral pathogenicity, and Omicron evolved to exhibit increased transmissibility and attenuated pathogenicity.

Abstract Social capital—the strength of an individual’s social network and community—has been identified as a potential determinant of outcomes ranging from education to health1–8. However, efforts to understand what types of social capital matter for these outcomes have been hindered by a lack of social network data. Here, in the first of a pair of papers9, we use data on 21 billion friendships from Facebook to study social capital. We measure and analyse three types of social capital by ZIP (postal) code in the United States: (1) connectedness between different types of people, such as those with low versus high socioeconomic status (SES); (2) social cohesion, such as the extent of cliques in friendship networks; and (3) civic engagement, such as rates of volunteering. These measures vary substantially across areas, but are not highly correlated with each other. We demonstrate the importance of distinguishing these forms of social capital by analysing their associations with economic mobility across areas. The share of high-SES friends among individuals with low SES—which we term…

Abstract Metformin, the most prescribed antidiabetic medicine, has shown other benefits such as anti-ageing and anticancer effects 1–4 . For clinical doses of metformin, AMP-activated protein kinase (AMPK) has a major role in its mechanism of action 4,5 ; however, the direct molecular target of metformin remains unknown. Here we show that clinically relevant concentrations of metformin inhibit the lysosomal proton pump v-ATPase, which is a central node for AMPK activation following glucose starvation 6 . We synthesize a photoactive metformin probe and identify PEN2, a subunit of γ-secretase 7 , as a binding partner of metformin with a dissociation constant at micromolar levels. Metformin-bound PEN2 forms a complex with ATP6AP1, a subunit of the v-ATPase 8 , which leads to the inhibition of v-ATPase and the activation of AMPK without effects on cellular AMP levels. Knockout of PEN2 or re-introduction of a PEN2 mutant that does not bind ATP6AP1 blunts AMPK activation. In vivo, liver-specific knockout of Pen2 abolishes metformin-mediated reduction of hepatic fat content, whereas…

Abstract The SARS-CoV-2 Omicron variant (B.1.1.529) has multiple spike protein mutations 1,2 that contribute to viral escape from antibody neutralization 3–6 and reduce vaccine protection from infection 7,8 . The extent to which other components of the adaptive response such as T cells may still target Omicron and contribute to protection from severe outcomes is unknown. Here we assessed the ability of T cells to react to Omicron spike protein in participants who were vaccinated with Ad26.CoV2.S or BNT162b2, or unvaccinated convalescent COVID-19 patients ( n = 70). Between 70% and 80% of the CD4 + and CD8 + T cell response to spike was maintained across study groups. Moreover, the magnitude of Omicron cross-reactive T cells was similar for Beta (B.1.351) and Delta (B.1.617.2) variants, despite Omicron harbouring considerably more mutations. In patients who were hospitalized with Omicron infections ( n = 19), there were comparable T cell responses to ancestral spike, nucleocapsid and membrane proteins to those in patients hospitalized in previous waves dominated by the ancestral,…

Abstract Covalent organic frameworks (COFs) are distinguished from other organic polymers by their crystallinity 1–3 , but it remains challenging to obtain robust, highly crystalline COFs because the framework-forming reactions are poorly reversible 4,5 . More reversible chemistry can improve crystallinity 6–9 , but this typically yields COFs with poor physicochemical stability and limited application scope 5 . Here we report a general and scalable protocol to prepare robust, highly crystalline imine COFs, based on an unexpected framework reconstruction. In contrast to standard approaches in which monomers are initially randomly aligned, our method involves the pre-organization of monomers using a reversible and removable covalent tether, followed by confined polymerization. This reconstruction route produces reconstructed COFs with greatly enhanced crystallinity and much higher porosity by means of a simple vacuum-free synthetic procedure. The increased crystallinity in the reconstructed COFs improves charge carrier transport, leading to sacrificial photocatalytic hydrogen…

Abstract Previous genome-wide association studies (GWASs) of stroke — the second leading cause of death worldwide — were conducted predominantly in populations of European ancestry 1,2 . Here, in cross-ancestry GWAS meta-analyses of 110,182 patients who have had a stroke (five ancestries, 33% non-European) and 1,503,898 control individuals, we identify association signals for stroke and its subtypes at 89 (61 new) independent loci: 60 in primary inverse-variance-weighted analyses and 29 in secondary meta-regression and multitrait analyses. On the basis of internal cross-ancestry validation and an independent follow-up in 89,084 additional cases of stroke (30% non-European) and 1,013,843 control individuals, 87% of the primary stroke risk loci and 60% of the secondary stroke risk loci were replicated ( P < 0.05). Effect sizes were highly correlated across ancestries. Cross-ancestry fine-mapping, in silico mutagenesis analysis 3 , and transcriptome-wide and proteome-wide association analyses revealed putative causal genes (such as SH3PXD2A and FURIN ) and variants (such as at GRK5 and…

Abstract Although circumstantial evidence supports enhanced Toll-like receptor 7 (TLR7) signalling as a mechanism of human systemic autoimmune disease 1–7 , evidence of lupus-causing TLR7 gene variants is lacking. Here we describe human systemic lupus erythematosus caused by a TLR7 gain-of-function variant. TLR7 is a sensor of viral RNA 8 , 9 and binds to guanosine 10 – 12 . We identified a de novo, previously undescribed missense TLR7 Y264H variant in a child with severe lupus and additional variants in other patients with lupus. The TLR7 Y264H variant selectively increased sensing of guanosine and 2',3'-cGMP 10–12 , and was sufficient to cause lupus when introduced into mice. We show that enhanced TLR7 signalling drives aberrant survival of B cell receptor (BCR)-activated B cells, and in a cell-intrinsic manner, accumulation of CD11c + age-associated B cells and germinal centre B cells. Follicular and extrafollicular helper T cells were also increased but these phenotypes were cell-extrinsic. Deficiency of MyD88 (an adaptor protein downstream of TLR7) rescued autoimmunity,…

. Here we describe a general solution to this problem that starts with a broad exploration of the vast space of possible binding modes to a selected region of a protein surface, and then intensifies the search in the vicinity of the most promising binding modes. We demonstrate the broad applicability of this approach through the de novo design of binding proteins to 12 diverse protein targets with different shapes and surface properties. Biophysical characterization shows that the binders, which are all smaller than 65 amino acids, are hyperstable and, following experimental optimization, bind their targets with nanomolar to picomolar affinities. We succeeded in solving crystal structures of five of the binder-target complexes, and all five closely match the corresponding computational design models. Experimental data on nearly half a million computational designs and hundreds of thousands of point mutants provide detailed feedback on the strengths and limitations of the method and of our current understanding of protein-protein interactions, and should guide improvements of both.…

A hallmark pathological feature of neurodegenerative diseases amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD) is the depletion of RNA-binding protein TDP-43 from the nucleus of neurons in the brain and spinal cord. A major function of TDP-43 is as a repressor of cryptic exon inclusion during RNA splicing. Single nucleotide polymorphisms (SNPs) in UNC13A are among the strongest genome-wide association study (GWAS) hits associated with FTD/ALS in humans, but how those variants increase risk for disease is unknown. Here we show that TDP-43 represses a cryptic exon splicing event in UNC13A . Loss of TDP-43 from the nucleus in human brain, neuronal cell lines, and iPSC-derived motor neurons resulted in the inclusion of a cryptic exon in UNC13A mRNA and reduced UNC13A protein expression. Remarkably, the top variants associated with FTD/ALS risk in humans are located in the cryptic exon harboring intron itself and we show that they increase UNC13A cryptic exon splicing in the face of TDP-43 dysfunction. Together, our data provide a direct functional link between one…

Abstract Comprehensive genome annotation is essential to understand the impact of clinically relevant variants. However, the absence of a standard for clinical reporting and browser display complicates the process of consistent interpretation and reporting. To address these challenges, Ensembl/GENCODE 1 and RefSeq 2 launched a joint initiative, the Matched Annotation from NCBI and EMBL-EBI (MANE) collaboration, to converge on human gene and transcript annotation and to jointly define a high-value set of transcripts and corresponding proteins. Here, we describe the MANE transcript sets for use as universal standards for variant reporting and browser display. The MANE Select set identifies a representative transcript for each human protein-coding gene, whereas the MANE Plus Clinical set provides additional transcripts at loci where the Select transcripts alone are not sufficient to report all currently known clinical variants. Each MANE transcript represents an exact match between the exonic sequences of an Ensembl/GENCODE transcript and its counterpart in RefSeq such that the…

Abstract Obtaining a burning plasma is a critical step towards self-sustaining fusion energy 1 . A burning plasma is one in which the fusion reactions themselves are the primary source of heating in the plasma, which is necessary to sustain and propagate the burn, enabling high energy gain. After decades of fusion research, here we achieve a burning-plasma state in the laboratory. These experiments were conducted at the US National Ignition Facility, a laser facility delivering up to 1.9 megajoules of energy in pulses with peak powers up to 500 terawatts. We use the lasers to generate X-rays in a radiation cavity to indirectly drive a fuel-containing capsule via the X-ray ablation pressure, which results in the implosion process compressing and heating the fuel via mechanical work. The burning-plasma state was created using a strategy to increase the spatial scale of the capsule 2,3 through two different implosion concepts 4–7 . These experiments show fusion self-heating in excess of the mechanical work injected into the implosions, satisfying several burning-plasma metrics 3,8 .…

Abstract Risk management has reduced vulnerability to floods and droughts globally 1,2 , yet their impacts are still increasing 3 . An improved understanding of the causes of changing impacts is therefore needed, but has been hampered by a lack of empirical data 4,5 . On the basis of a global dataset of 45 pairs of events that occurred within the same area, we show that risk management generally reduces the impacts of floods and droughts but faces difficulties in reducing the impacts of unprecedented events of a magnitude not previously experienced. If the second event was much more hazardous than the first, its impact was almost always higher. This is because management was not designed to deal with such extreme events: for example, they exceeded the design levels of levees and reservoirs. In two success stories, the impact of the second, more hazardous, event was lower, as a result of improved risk management governance and high investment in integrated management. The observed difficulty of managing unprecedented events is alarming, given that more extreme hydrological events are…

Abstract COVID-19, which is caused by infection with SARS-CoV-2, is characterized by lung pathology and extrapulmonary complications 1,2 . Type I interferons (IFNs) have an essential role in the pathogenesis of COVID-19 (refs 3–5 ). Although rapid induction of type I IFNs limits virus propagation, a sustained increase in the levels of type I IFNs in the late phase of the infection is associated with aberrant inflammation and poor clinical outcome 5–17 . Here we show that the cyclic GMP-AMP synthase (cGAS)–stimulator of interferon genes (STING) pathway, which controls immunity to cytosolic DNA, is a critical driver of aberrant type I IFN responses in COVID-19 (ref. 18 ). Profiling COVID-19 skin manifestations, we uncover a STING-dependent type I IFN signature that is primarily mediated by macrophages adjacent to areas of endothelial cell damage. Moreover, cGAS–STING activity was detected in lung samples from patients with COVID-19 with prominent tissue destruction, and was associated with type I IFN responses. A lung-on-chip model revealed that, in addition to macrophages, infection…

Abstract Variants of UNC13A , a critical gene for synapse function, increase the risk of amyotrophic lateral sclerosis and frontotemporal dementia 1–3 , two related neurodegenerative diseases defined by mislocalization of the RNA-binding protein TDP-43 4,5 . Here we show that TDP-43 depletion induces robust inclusion of a cryptic exon in UNC13A , resulting in nonsense-mediated decay and loss of UNC13A protein. Two common intronic UNC13A polymorphisms strongly associated with amyotrophic lateral sclerosis and frontotemporal dementia risk overlap with TDP-43 binding sites. These polymorphisms potentiate cryptic exon inclusion, both in cultured cells and in brains and spinal cords from patients with these conditions. Our findings, which demonstrate a genetic link between loss of nuclear TDP-43 function and disease, reveal the mechanism by which UNC13A variants exacerbate the effects of decreased TDP-43 function. They further provide a promising therapeutic target for TDP-43 proteinopathies.

Abstract Thermal insulation under extreme conditions requires materials that can withstand complex thermomechanical stress and retain excellent thermal insulation properties at temperatures exceeding 1,000 degrees Celsius 1–3 . Ceramic aerogels are attractive thermal insulating materials; however, at very high temperatures, they often show considerably increased thermal conductivity and limited thermomechanical stability that can lead to catastrophic failure 4–6 . Here we report a multiscale design of hypocrystalline zircon nanofibrous aerogels with a zig-zag architecture that leads to exceptional thermomechanical stability and ultralow thermal conductivity at high temperatures. The aerogels show a near-zero Poisson’s ratio (3.3 × 10 −4 ) and a near-zero thermal expansion coefficient (1.2 × 10 −7 per degree Celsius), which ensures excellent structural flexibility and thermomechanical properties. They show high thermal stability with ultralow strength degradation (less than 1 per cent) after sharp thermal shocks, and a high working temperature (up to 1,300 degrees Celsius). By…

Abstract High-fidelity control of quantum bits is paramount for the reliable execution of quantum algorithms and for achieving fault tolerance—the ability to correct errors faster than they occur 1 . The central requirement for fault tolerance is expressed in terms of an error threshold. Whereas the actual threshold depends on many details, a common target is the approximately 1% error threshold of the well-known surface code 2,3 . Reaching two-qubit gate fidelities above 99% has been a long-standing major goal for semiconductor spin qubits. These qubits are promising for scaling, as they can leverage advanced semiconductor technology 4 . Here we report a spin-based quantum processor in silicon with single-qubit and two-qubit gate fidelities, all of which are above 99.5%, extracted from gate-set tomography. The average single-qubit gate fidelities remain above 99% when including crosstalk and idling errors on the neighbouring qubit. Using this high-fidelity gate set, we execute the demanding task of calculating molecular ground-state energies using a variational quantum eigensolver…

Abstract Missing heritability in genome-wide association studies defines a major problem in genetic analyses of complex biological traits 1,2 . The solution to this problem is to identify all causal genetic variants and to measure their individual contributions 3,4 . Here we report a graph pangenome of tomato constructed by precisely cataloguing more than 19 million variants from 838 genomes, including 32 new reference-level genome assemblies. This graph pangenome was used for genome-wide association study analyses and heritability estimation of 20,323 gene-expression and metabolite traits. The average estimated trait heritability is 0.41 compared with 0.33 when using the single linear reference genome. This 24% increase in estimated heritability is largely due to resolving incomplete linkage disequilibrium through the inclusion of additional causal structural variants identified using the graph pangenome. Moreover, by resolving allelic and locus heterogeneity, structural variants improve the power to identify genetic factors underlying agronomically important traits leading to, for…

Abstract Myelin is required for the function of neuronal axons in the central nervous system, but the mechanisms that support myelin health are unclear. Although macrophages in the central nervous system have been implicated in myelin health 1 , it is unknown which macrophage populations are involved and which aspects they influence. Here we show that resident microglia are crucial for the maintenance of myelin health in adulthood in both mice and humans. We demonstrate that microglia are dispensable for developmental myelin ensheathment. However, they are required for subsequent regulation of myelin growth and associated cognitive function, and for preservation of myelin integrity by preventing its degeneration. We show that loss of myelin health due to the absence of microglia is associated with the appearance of a myelinating oligodendrocyte state with altered lipid metabolism. Moreover, this mechanism is regulated through disruption of the TGFβ1–TGFβR1 axis. Our findings highlight microglia as promising therapeutic targets for conditions in which myelin growth and integrity are…

ABSTRACT The gains and losses of DNA that emerge as a consequence of mitotic errors and chromosomal instability are prevalent in cancer. These copy number alterations contribute to cancer initiaition, progression and therapeutic resistance. Here, we present a conceptual framework for examining the patterns of copy number alterations in human cancer using whole-genome sequencing, whole-exome sequencing, and SNP6 microarray data making it widely applicable to diverse datasets. Deploying this framework to 9,873 cancers representing 33 human cancer types from the TCGA project revealed a set of 19 copy number signatures that explain the copy number patterns of 93% of TCGA samples. 15 copy number signatures were attributed to biological processes of whole-genome doubling, aneuploidy, loss of heterozygosity, homologous recombination deficiency, and chromothripsis. The aetiology of four copy number signatures are unexplained and some cancer types have unique patterns of amplicon signatures associated with extrachromosomal DNA, disease-specific survival, and gains of proto-oncogenes such as…

, but the reason for such abrupt functional decline after 70 years of age remains unclear. Here we sequenced 3,579 genomes from single cell-derived colonies of haematopoietic cells across 10 human subjects from 0 to 81 years of age. Haematopoietic stem cells or multipotent progenitors (HSC/MPPs) accumulated a mean of 17 mutations per year after birth and lost 30 base pairs per year of telomere length. Haematopoiesis in adults less than 65 years of age was massively polyclonal, with high clonal diversity and a stable population of 20,000-200,000 HSC/MPPs contributing evenly to blood production. By contrast, haematopoiesis in individuals aged over 75 showed profoundly decreased clonal diversity. In each of the older subjects, 30-60% of haematopoiesis was accounted for by 12-18 independent clones, each contributing 1-34% of blood production. Most clones had begun their expansion before the subject was 40 years old, but only 22% had known driver mutations. Genome-wide selection analysis estimated that between 1 in 34 and 1 in 12 non-synonymous mutations were drivers, accruing at…

Abstract Nirmatrelvir, an oral antiviral targeting the 3CL protease of SARS-CoV-2, has been demonstrated to be clinically useful against COVID-19 (refs. 1,2 ). However, because SARS-CoV-2 has evolved to become resistant to other therapeutic modalities 3–9 , there is a concern that the same could occur for nirmatrelvir. Here we examined this possibility by in vitro passaging of SARS-CoV-2 in nirmatrelvir using two independent approaches, including one on a large scale. Indeed, highly resistant viruses emerged from both and their sequences showed a multitude of 3CL protease mutations. In the experiment peformed with many replicates, 53 independent viral lineages were selected with mutations observed at 23 different residues of the enzyme. Nevertheless, several common mutational pathways to nirmatrelvir resistance were preferred, with a majority of the viruses descending from T21I, P252L or T304I as precursor mutations. Construction and analysis of 13 recombinant SARS-CoV-2 clones showed that these mutations mediated only low-level resistance, whereas greater resistance required…

. However, organoids lack the connectivity that exists in vivo, which limits maturation and makes integration with other circuits that control behaviour impossible. Here we show that human stem cell-derived cortical organoids transplanted into the somatosensory cortex of newborn athymic rats develop mature cell types that integrate into sensory and motivation-related circuits. MRI reveals post-transplantation organoid growth across multiple stem cell lines and animals, whereas single-nucleus profiling shows progression of corticogenesis and the emergence of activity-dependent transcriptional programs. Indeed, transplanted cortical neurons display more complex morphological, synaptic and intrinsic membrane properties than their in vitro counterparts, which enables the discovery of defects in neurons derived from individuals with Timothy syndrome. Anatomical and functional tracings show that transplanted organoids receive thalamocortical and corticocortical inputs, and in vivo recordings of neural activity demonstrate that these inputs can produce sensory responses in human cells.…

As SARS-CoV-2 continues to spread and evolve, detecting emerging variants early is critical for public health interventions. Inferring lineage prevalence by clinical testing is infeasible at scale, especially in areas with limited resources, participation, or testing/sequencing capacity, which can also introduce biases. SARS-CoV-2 RNA concentration in wastewater successfully tracks regional infection dynamics and provides less biased abundance estimates than clinical testing. Tracking virus genomic sequences in wastewater would improve community prevalence estimates and detect emerging variants. However, two factors limit wastewater-based genomic surveillance: low-quality sequence data and inability to estimate relative lineage abundance in mixed samples. Here, we resolve these critical issues to perform a high-resolution, 295-day wastewater and clinical sequencing effort, in the controlled environment of a large university campus and the broader context of the surrounding county. We develop and deploy improved virus concentration protocols and deconvolution software that fully…

Abstract Public databases contain a planetary collection of nucleic acid sequences, but their systematic exploration has been inhibited by a lack of efficient methods for searching this corpus, now exceeding multiple petabases and growing exponentially [1, 2]. We developed a cloud computing infrastructure, Serratus , to enable ultra-high throughput sequence alignment at the petabase scale. We searched 5.7 million biologically diverse samples (10.2 petabases) for the hallmark gene RNA dependent RNA polymerase, identifying well over 10 5 novel RNA viruses and thereby expanding the number of known species by roughly an order of magnitude. We characterised novel viruses related to coronaviruses and to hepatitis δ virus, respectively and explored their environmental reservoirs. To catalyse a new era of viral discovery, we established a free and comprehensive database of these data and tools. Expanding the known sequence diversity of viruses can reveal the evolutionary origins of emerging pathogens and improve pathogen surveillance for the anticipation and mitigation of future pandemics.