Deep Mutational Scanning of SARS-CoV-2 Receptor Binding Domain Reveals Constraints on Folding and ACE2 Binding.

The receptor binding domain (RBD) of the SARS-CoV-2 spike glycoprotein mediates viral attachment to ACE2 receptor and is a major determinant of host range and a dominant target of neutralizing antibodies. Here, we experimentally measure how all amino acid mutations to the RBD affect expression of folded protein and its affinity for ACE2. Most mutations are deleterious for RBD expression and ACE2 binding, and we identify constrained regions on the RBD's surface that may be desirable targets for vaccines and antibody-based therapeutics. But a substantial number of mutations are well tolerated or even enhance ACE2 binding, including at ACE2 interface residues that vary across SARS-related coronaviruses. However, we find no evidence that these ACE2-affinity-enhancing mutations have been selected in current SARS-CoV-2 pandemic isolates. We present an interactive visualization and open analysis pipeline to facilitate use of our dataset for vaccine design and functional annotation of mutations observed during viral surveillance.

A White-Box Machine Learning Approach for Revealing Antibiotic Mechanisms of Action.

Current machine learning techniques enable robust association of biological signals with measured phenotypes, but these approaches are incapable of identifying causal relationships. Here, we develop an integrated "white-box" biochemical screening, network modeling, and machine learning approach for revealing causal mechanisms and apply this approach to understanding antibiotic efficacy. We counter-screen diverse metabolites against bactericidal antibiotics in Escherichia coli and simulate their corresponding metabolic states using a genome-scale metabolic network model. Regression of the measured screening data on model simulations reveals that purine biosynthesis participates in antibiotic lethality, which we validate experimentally. We show that antibiotic-induced adenine limitation increases ATP demand, which elevates central carbon metabolism activity and oxygen consumption, enhancing the killing effects of antibiotics. This work demonstrates how prospective network modeling can couple with machine learning to identify complex causal mechanisms underlying drug efficacy.

Organoid Modeling of the Tumor Immune Microenvironment.

In vitro cancer cultures, including three-dimensional organoids, typically contain exclusively neoplastic epithelium but require artificial reconstitution to recapitulate the tumor microenvironment (TME). The co-culture of primary tumor epithelia with endogenous, syngeneic tumor-infiltrating lymphocytes (TILs) as a cohesive unit has been particularly elusive. Here, an air-liquid interface (ALI) method propagated patient-derived organoids (PDOs) from >100 human biopsies or mouse tumors in syngeneic immunocompetent hosts as tumor epithelia with native embedded immune cells (T, B, NK, macrophages). Robust droplet-based, single-cell simultaneous determination of gene expression and immune repertoire indicated that PDO TILs accurately preserved the original tumor T cell receptor (TCR) spectrum. Crucially, human and murine PDOs successfully modeled immune checkpoint blockade (ICB) with anti-PD-1- and/or anti-PD-L1 expanding and activating tumor antigen-specific TILs and eliciting tumor cytotoxicity. Organoid-based propagation of primary tumor epithelium en bloc with endogenous immune stroma should enable immuno-oncology investigations within the TME and facilitate personalized immunotherapy testing.

Inborn Errors of RNA Lariat Metabolism in Humans with Brainstem Viral Infection.

Viruses that are typically benign sometimes invade the brainstem in otherwise healthy children. We report bi-allelic DBR1 mutations in unrelated patients from different ethnicities, each of whom had brainstem infection due to herpes simplex virus 1 (HSV1), influenza virus, or norovirus. DBR1 encodes the only known RNA lariat debranching enzyme. We show that DBR1 expression is ubiquitous, but strongest in the spinal cord and brainstem. We also show that all DBR1 mutant alleles are severely hypomorphic, in terms of expression and function. The fibroblasts of DBR1-mutated patients contain higher RNA lariat levels than control cells, this difference becoming even more marked during HSV1 infection. Finally, we show that the patients' fibroblasts are highly susceptible to HSV1. RNA lariat accumulation and viral susceptibility are rescued by wild-type DBR1. Autosomal recessive, partial DBR1 deficiency underlies viral infection of the brainstem in humans through the disruption of tissue-specific and cell-intrinsic immunity to viruses.

Clinical practice guideline and expert consensus recommendations for rehabilitation among children with cancer: A systematic review.

Increased attention to the rehabilitation needs of children with cancer is vital to enhance health, quality-of-life, and productivity outcomes. Among adults with cancer, rehabilitation recommendations are frequently incorporated into guidelines, but the extent to which recommendations exist for children is unknown. Reports included in this systematic review are guideline or expert consensus reports containing recommendations related to rehabilitation referral, evaluation, and/or intervention for individuals diagnosed with cancer during childhood (younger than 18 years). Eligible reports were published in English from January 2000 to August 2022. Through database searches, 42,982 records were identified; 62 records were identified through citation and website searching. Twenty-eight reports were included in the review: 18 guidelines and 10 expert consensus reports. Rehabilitation recommendations were identified in disease-specific (e.g., acute lymphoblastic leukemia), impairment-specific (e.g., fatigue, neurocognition, pain), adolescent and young adult, and long-term follow-up reports. Example recommendations included physical activity and energy-conservation techniques to address fatigue, referral to physical therapy for chronic pain management, ongoing psychosocial surveillance, and referral to speech-language pathology for those with hearing loss. High-level evidence supported rehabilitation recommendations for long-term follow-up care, fatigue, and psychosocial/mental health screening. Few intervention recommendations were included in guideline and consensus reports. In this developing field, it is critical to include pediatric oncology rehabilitation providers in guideline and consensus development initiatives. This review enhances the availability and clarity of rehabilitation-relevant guidelines that can help prevent and mitigate cancer-related disability among children by supporting access to rehabilitation services.

Synthetic biology looks good on paper.

Tremendous progress has been made in the design and implementation of synthetic gene circuits, but real-world applications of such circuits have been limited. Cell-free circuits embedded on paper developed by Pardee et al. promise to deliver specific and rapid diagnostics on a low-cost, highly scalable platform.

Molecular fate-mapping of serum antibody responses to repeat immunization.

The protective efficacy of serum antibodies results from the interplay of antigen-specific B cell clones of different affinities and specificities. These cellular dynamics underlie serum-level phenomena such as original antigenic sin (OAS)-a proposed propensity of the immune system to rely repeatedly on the first cohort of B cells engaged by an antigenic stimulus when encountering related antigens, in detriment to the induction of de novo responses1-5. OAS-type suppression of new, variant-specific antibodies may pose a barrier to vaccination against rapidly evolving viruses such as influenza and SARS-CoV-26,7. Precise measurement of OAS-type suppression is challenging because cellular and temporal origins cannot readily be ascribed to antibodies in circulation; its effect on subsequent antibody responses therefore remains unclear5,8. Here we introduce a molecular fate-mapping approach with which serum antibodies derived from specific cohorts of B cells can be differentially detected. We show that serum responses to sequential homologous boosting derive overwhelmingly from primary cohort B cells, while later induction of new antibody responses from naive B cells is strongly suppressed. Such 'primary addiction' decreases sharply as a function of antigenic distance, allowing reimmunization with divergent viral glycoproteins to produce de novo antibody responses targeting epitopes that are absent from the priming variant. Our findings have implications for the understanding of OAS and for the design and testing of vaccines against evolving pathogens.

Quantum critical dynamics in a 5,000-qubit programmable spin glass.

Experiments on disordered alloys1-3 suggest that spin glasses can be brought into low-energy states faster by annealing quantum fluctuations than by conventional thermal annealing. Owing to the importance of spin glasses as a paradigmatic computational testbed, reproducing this phenomenon in a programmable system has remained a central challenge in quantum optimization4-13. Here we achieve this goal by realizing quantum-critical spin-glass dynamics on thousands of qubits with a superconducting quantum annealer. We first demonstrate quantitative agreement between quantum annealing and time evolution of the Schrödinger equation in small spin glasses. We then measure dynamics in three-dimensional spin glasses on thousands of qubits, for which classical simulation of many-body quantum dynamics is intractable. We extract critical exponents that clearly distinguish quantum annealing from the slower stochastic dynamics of analogous Monte Carlo algorithms, providing both theoretical and experimental support for large-scale quantum simulation and a scaling advantage in energy optimization.

ACE2 binding is an ancestral and evolvable trait of sarbecoviruses.

Two different sarbecoviruses have caused major human outbreaks in the past two decades1,2. Both of these sarbecoviruses, SARS-CoV-1 and SARS-CoV-2, engage ACE2 through the spike receptor-binding domain2-6. However, binding to ACE2 orthologues of humans, bats and other species has been observed only sporadically among the broader diversity of bat sarbecoviruses7-11. Here we use high-throughput assays12 to trace the evolutionary history of ACE2 binding across a diverse range of sarbecoviruses and ACE2 orthologues. We find that ACE2 binding is an ancestral trait of sarbecovirus receptor-binding domains that has subsequently been lost in some clades. Furthermore, we reveal that bat sarbecoviruses from outside Asia can bind to ACE2. Moreover, ACE2 binding is highly evolvable-for many sarbecovirus receptor-binding domains, there are single amino-acid mutations that enable binding to new ACE2 orthologues. However, the effects of individual mutations can differ considerably between viruses, as shown by the N501Y mutation, which enhances the human ACE2-binding affinity of several SARS-CoV-2 variants of concern12 but substantially decreases it for SARS-CoV-1. Our results point to the deep ancestral origin and evolutionary plasticity of ACE2 binding, broadening the range of sarbecoviruses that should be considered to have spillover potential.

Living with Two Genomes: Grafting and Its Implications for Plant Genome-to-Genome Interactions, Phenotypic Variation, and Evolution.

Plant genomes interact when genetically distinct individuals join, or are joined, together. Individuals can fuse in three contexts: artificial grafts, natural grafts, and host-parasite interactions. Artificial grafts have been studied for decades and are important platforms for studying the movement of RNA, DNA, and protein. Yet several mysteries about artificial grafts remain, including the factors that contribute to graft incompatibility, the prevalence of genetic and epigenetic modifications caused by exchanges between graft partners, and the long-term effects of these modifications on phenotype. Host-parasite interactions also lead to the exchange of materials, and RNA exchange actively contributes to an ongoing arms race between parasite virulence and host resistance. Little is known about natural grafts except that they can be frequent and may provide opportunities for evolutionary innovation through genome exchange. In this review, we survey our current understanding about these three mechanisms of contact, the genomic interactions that result, and the potential evolutionary implications.

Detection of a SARS-CoV-2 variant of concern in South Africa.

Continued uncontrolled transmission of SARS-CoV-2 in many parts of the world is creating conditions for substantial evolutionary changes to the virus1,2. Here we describe a newly arisen lineage of SARS-CoV-2 (designated 501Y.V2; also known as B.1.351 or 20H) that is defined by eight mutations in the spike protein, including three substitutions (K417N, E484K and N501Y) at residues in its receptor-binding domain that may have functional importance3-5. This lineage was identified in South Africa after the first wave of the epidemic in a severely affected metropolitan area (Nelson Mandela Bay) that is located on the coast of the Eastern Cape province. This lineage spread rapidly, and became dominant in Eastern Cape, Western Cape and KwaZulu-Natal provinces within weeks. Although the full import of the mutations is yet to be determined, the genomic data-which show rapid expansion and displacement of other lineages in several regions-suggest that this lineage is associated with a selection advantage that most plausibly results from increased transmissibility or immune escape6-8.

SARS-CoV-2 RBD antibodies that maximize breadth and resistance to escape.

An ideal therapeutic anti-SARS-CoV-2 antibody would resist viral escape1-3, have activity against diverse sarbecoviruses4-7, and be highly protective through viral neutralization8-11 and effector functions12,13. Understanding how these properties relate to each other and vary across epitopes would aid the development of therapeutic antibodies and guide vaccine design. Here we comprehensively characterize escape, breadth and potency across a panel of SARS-CoV-2 antibodies targeting the receptor-binding domain (RBD). Despite a trade-off between in vitro neutralization potency and breadth of sarbecovirus binding, we identify neutralizing antibodies with exceptional sarbecovirus breadth and a corresponding resistance to SARS-CoV-2 escape. One of these antibodies, S2H97, binds with high affinity across all sarbecovirus clades to a cryptic epitope and prophylactically protects hamsters from viral challenge. Antibodies that target the angiotensin-converting enzyme 2 (ACE2) receptor-binding motif (RBM) typically have poor breadth and are readily escaped by mutations despite high neutralization potency. Nevertheless, we also characterize a potent RBM antibody (S2E128) with breadth across sarbecoviruses related to SARS-CoV-2 and a high barrier to viral escape. These data highlight principles underlying variation in escape, breadth and potency among antibodies that target the RBD, and identify epitopes and features to prioritize for therapeutic development against the current and potential future pandemics.

A Fast, Reproducible, High-throughput Variant Calling Workflow for Population Genomics.

The increasing availability of genomic resequencing data sets and high-quality reference genomes across the tree of life present exciting opportunities for comparative population genomic studies. However, substantial challenges prevent the simple reuse of data across different studies and species, arising from variability in variant calling pipelines, data quality, and the need for computationally intensive reanalysis. Here, we present snpArcher, a flexible and highly efficient workflow designed for the analysis of genomic resequencing data in nonmodel organisms. snpArcher provides a standardized variant calling pipeline and includes modules for variant quality control, data visualization, variant filtering, and other downstream analyses. Implemented in Snakemake, snpArcher is user-friendly, reproducible, and designed to be compatible with high-performance computing clusters and cloud environments. To demonstrate the flexibility of this pipeline, we applied snpArcher to 26 public resequencing data sets from nonmammalian vertebrates. These variant data sets are hosted publicly to enable future comparative population genomic analyses. With its extensibility and the availability of public data sets, snpArcher will contribute to a broader understanding of genetic variation across species by facilitating the rapid use and reuse of large genomic data sets.

Unraveling the features of somatic transposition in the Drosophila intestine.

Transposable elements (TEs) play a significant role in evolution, contributing to genetic variation. However, TE mobilization in somatic cells is not well understood. Here, we address the prevalence of transposition in a somatic tissue, exploiting the Drosophila midgut as a model. Using whole-genome sequencing of in vivo clonally expanded gut tissue, we have mapped hundreds of high-confidence somatic TE integration sites genome-wide. We show that somatic retrotransposon insertions are associated with inactivation of the tumor suppressor Notch, likely contributing to neoplasia formation. Moreover, applying Oxford Nanopore long-read sequencing technology we provide evidence for tissue-specific differences in retrotransposition. Comparing somatic TE insertional activity with transcriptomic and small RNA sequencing data, we demonstrate that transposon mobility cannot be simply predicted by whole tissue TE expression levels or by small RNA pathway activity. Finally, we reveal that somatic TE insertions in the adult fly intestine are enriched in genic regions and in transcriptionally active chromatin. Together, our findings provide clear evidence of ongoing somatic transposition in Drosophila and delineate previously unknown features underlying somatic TE mobility in vivo.

Illuminati: a form of gene expression plasticity in Drosophila neural stem cells.

While testing for genome instability in Drosophila as reported by unscheduled upregulation of UAS-GFP in cells that co-express GAL80 and GAL4, we noticed that, as expected, background levels were low in most developing tissues. However, GFP-positive clones were frequent in the larval brain. Most of these clones originated from central brain neural stem cells. Using imaging-based approaches and genome sequencing, we show that these unscheduled clones do not result from chromosome loss or mutations in GAL80. We have named this phenomenon 'Illuminati'. Illuminati is strongly enhanced in brat tumors and is also sensitive to environmental conditions such as food content and temperature. Illuminati is suppressed by Su(var)2-10, but it is not significantly affected by several modifiers of position effect variegation or Gal4::UAS variegation. We conclude that Illuminati identifies a previously unknown type of functional instability that may have important implications in development and disease.

Validation of Ligand Tetramers for the Detection of Antigen-Specific Lymphocytes.

The study of Ag-specific lymphocytes has been a key advancement in immunology over the past few decades. The development of multimerized probes containing Ags, peptide:MHC complexes, or other ligands was one innovation allowing the direct study of Ag-specific lymphocytes by flow cytometry. Although these types of study are now common and performed by thousands of laboratories, quality control and assessment of probe quality are often minimal. In fact, many of these types of probe are made in-house, and protocols vary between laboratories. Although peptide:MHC multimers can often be obtained from commercial sources or core facilities, few such services exist for Ag multimers. To ensure high quality and consistency with ligand probes, we have developed an easy and robust multiplexed approach using commercially available beads able to bind Abs specific for the ligand of interest. Using this assay, we have sensitively assessed the performance of peptide:MHC and Ag tetramers and have found considerable batch-to-batch variability in performance and stability over time more easily than using murine or human cell-based assays. This bead-based assay can also reveal common production errors such as miscalculation of Ag concentration. This work could set the stage for the development of standardized assays for all commonly used ligand probes to limit laboratory-to-laboratory technical variation and experimental failure caused by probe underperformance.

Metabolism in Pulmonary Hypertension.

Pulmonary arterial hypertension (PAH) is characterized by impaired regulation of pulmonary hemodynamics and vascular growth. Alterations of metabolism and bioenergetics are increasingly recognized as universal hallmarks of PAH, as metabolic abnormalities are identified in lungs and hearts of patients, animal models of the disease, and cells derived from lungs of patients. Mitochondria are the primary organelle critically mediating the complex and integrative metabolic pathways in bioenergetics, biosynthetic pathways, and cell signaling. Here, we review the alterations in metabolic pathways that are linked to the pathologic vascular phenotype of PAH, including abnormalities in glycolysis and glucose oxidation, fatty acid oxidation, glutaminolysis, arginine metabolism, one-carbon metabolism, the reducing and oxidizing cell environment, and the tricarboxylic acid cycle, as well as the effects of PAH-associated nuclear and mitochondrial mutations on metabolism. Understanding of the metabolic mechanisms underlying PAH provides important knowledge for the design of new therapeutics for treatment of patients.

Distinct roles of Bdnf I and Bdnf IV transcript variant expression in hippocampal neurons.

Brain-derived neurotrophic factor (Bdnf) plays a critical role in brain development, dendritic growth, synaptic plasticity, as well as learning and memory. The rodent Bdnf gene contains nine 5' non-coding exons (I-IXa), which are spliced to a common 3' coding exon (IX). Transcription of individual Bdnf variants, which all encode the same BDNF protein, is initiated at unique promoters upstream of each non-coding exon, enabling precise spatiotemporal and activity-dependent regulation of Bdnf expression. Although prior evidence suggests that Bdnf transcripts containing exon I (Bdnf I) or exon IV (Bdnf IV) are uniquely regulated by neuronal activity, the functional significance of different Bdnf transcript variants remains unclear. To investigate functional roles of activity-dependent Bdnf I and IV transcripts, we used a CRISPR activation system in which catalytically dead Cas9 fused to a transcriptional activator (VPR) is targeted to individual Bdnf promoters with single guide RNAs, resulting in transcript-specific Bdnf upregulation. Bdnf I upregulation is associated with gene expression changes linked to dendritic growth, while Bdnf IV upregulation is associated with genes that regulate protein catabolism. Upregulation of Bdnf I, but not Bdnf IV, increased mushroom spine density, volume, length, and head diameter, and also produced more complex dendritic arbors in cultured rat hippocampal neurons. In contrast, upregulation of Bdnf IV, but not Bdnf I, in the rat hippocampus attenuated contextual fear expression. Our data suggest that while Bdnf I and IV are both activity-dependent, BDNF produced from these promoters may serve unique cellular, synaptic, and behavioral functions.

Excess waitlist mortality among candidates for multivisceral liver-intestine transplant in acuity circle allocation.

Candidates for multivisceral transplant (MVT) have experienced decreased access to transplant in recent years. Using Organ Procurement and Transplantation Network data, transplant and waiting list outcomes for MVT (ie, liver-intestine, liver-intestine-pancreas, and liver-intestine-kidney-pancreas) candidates listed between February 4, 2018, and February 3, 2022, were analyzed, including model for end-stage liver disease/pediatric end-stage liver disease and exception scores by era (before and after acuity circle [AC] implementation on February 4, 2020) and age group (pediatric and adult). Of 284 MVT waitlist registrations (45.6% pediatric), fewer had exception points at listing post-AC compared to pre-AC (10.0% vs 19.1%), and they were less likely to receive transplant (19.1% vs 35.9% at 90 days; 35.7% vs 57.2% at 1 year). Of 177 MVT recipients, exception points at transplant were more common post-AC compared to pre-AC (30.8% vs 20.2%). Postpolicy, adult MVT candidates were more likely to be removed due to death/too sick compared with liver-alone candidates (13.5% vs 5.6% at 90 days; 24.2% vs 9.8% at 1 year), whereas no excess waitlist mortality was observed among pediatric MVT candidates. Under current allocation policy, multivisceral candidates experience inferior waitlist outcomes compared with liver-alone candidates. Clarification of guidance around submission and approval of multivisceral exception requests may help improve their access to transplantation and achieve equity between multivisceral and liver-alone candidates on the liver transplant waiting list.

OPTN/SRTR 2022 Annual Data Report: Liver.

In 2022, liver transplant activity continued to increase in the United States, with an all-time high of 9,527 transplants performed, representing a 52% increase over the past decade (2012-2022). Of these transplants, 8,924 (93.7%) were from deceased donors and 603 (6.3%) were from living donors. Liver transplant recipients were 94.5% adult and 5.5% pediatric. The overall size of the liver transplant waiting list contracted, with more patients being removed than added, although 10,548 adult patients still remained on the waiting list at the end of 2022. Alcohol-associated liver disease continued to be the leading diagnosis among both candidates and recipients, followed by metabolic dysfunction-associated steatohepatitis. Simultaneous liver-kidney transplant was the most common multiorgan combination, with 800 liver-kidney transplants performed in 2022; in addition, there were 303 new listings for kidney transplant via the safety net mechanism. Among adults added to the liver waiting list in 2021, 39.9% received a deceased donor liver transplant within 3 months; 45.7%, within 6 months; and 54.5%, within 1 year. Pretransplant mortality decreased to 12.3 deaths per 100 patient-years in 2022, although still 15.6% of removals from the waiting list were for death or being too sick for transplant. Graft and patient survival outcomes after deceased donor liver transplant improved, approximating pre-COVID-19 pandemic levels, with 5.1% mortality observed at 6 months; 6.8%, at 1 year; 12.7%, at 3 years; 19.8%, at 5 years; and 35.7%, at 10 years. Five-year graft and patient survival rates after living donor liver transplant exceeded those of deceased donor liver transplant. Candidates receiving model for end-stage liver disease exception points for hepatocellular carcinoma constituted 15.5% of transplants performed in 2022, with similar transplant rates and posttransplant outcomes compared to cases without hepatocellular carcinoma exception. In 2022, more pediatric liver transplant candidates were added to the waiting list and underwent transplant compared with either of the preceding 2 years, with an uptick in living donor liver transplant volume. Although pretransplant mortality has improved after the recent policy change prioritizing pediatric donors for pediatric recipients, still, in 2022, 50 children died or were removed from the waiting list for being too sick to undergo transplant. Posttransplant mortality among pediatric liver transplant recipients remained notable, with death occurring in 4.0% at 6 months, 6.0% at 1 year, 8.2% at 3 years, 9.8% at 5 years, and 13.9% at 10 years. Similar to adult living donor recipients, pediatric living donor recipients had better 5-year patient survival compared with deceased donor recipients.