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.

Robustness of Phylogenetic Inference to Model Misspecification Caused by Pairwise Epistasis.

Likelihood-based phylogenetic inference posits a probabilistic model of character state change along branches of a phylogenetic tree. These models typically assume statistical independence of sites in the sequence alignment. This is a restrictive assumption that facilitates computational tractability, but ignores how epistasis, the effect of genetic background on mutational effects, influences the evolution of functional sequences. We consider the effect of using a misspecified site-independent model on the accuracy of Bayesian phylogenetic inference in the setting of pairwise-site epistasis. Previous work has shown that as alignment length increases, tree reconstruction accuracy also increases. Here, we present a simulation study demonstrating that accuracy increases with alignment size even if the additional sites are epistatically coupled. We introduce an alignment-based test statistic that is a diagnostic for pairwise epistasis and can be used in posterior predictive checks.

Single-molecule long-read 16S sequencing to characterize the lung microbiome from mechanically ventilated patients with suspected pneumonia.

In critically ill patients, the development of pneumonia results in significant morbidity and mortality and additional health care costs. The accurate and rapid identification of the microbial pathogens in patients with pulmonary infections might lead to targeted antimicrobial therapy with potentially fewer adverse effects and lower costs. Major advances in next-generation sequencing (NGS) allow culture-independent identification of pathogens. The present study used NGS of essentially full-length PCR-amplified 16S ribosomal DNA from the bronchial aspirates of intubated patients with suspected pneumonia. The results from 61 patients demonstrated that sufficient DNA was obtained from 72% of samples, 44% of which (27 samples) yielded PCR amplimers suitable for NGS. Out of the 27 sequenced samples, only 20 had bacterial culture growth, while the microbiological and NGS identification of bacteria coincided in 17 (85%) of these samples. Despite the lack of bacterial growth in 7 samples that yielded amplimers and were sequenced, the NGS identified a number of bacterial species in these samples. Overall, a significant diversity of bacterial species was identified from the same genus as the predominant cultured pathogens. The numbers of NGS-identifiable bacterial genera were consistently higher than identified by standard microbiological methods. As technical advances reduce the processing and sequencing times, NGS-based methods will ultimately be able to provide clinicians with rapid, precise, culture-independent identification of bacterial, fungal, and viral pathogens and their antimicrobial sensitivity profiles.

Drug-regulated CD33-targeted CAR T cells control AML using clinically optimized rapamycin dosing.

Chimeric antigen receptor (CAR) designs that incorporate pharmacologic control are desirable; however, designs suitable for clinical translation are needed. We designed a fully human, rapamycin-regulated drug product for targeting CD33+ tumors called dimerizaing agent-regulated immunoreceptor complex (DARIC33). T cell products demonstrated target-specific and rapamycin-dependent cytokine release, transcriptional responses, cytotoxicity, and in vivo antileukemic activity in the presence of as little as 1 nM rapamycin. Rapamycin withdrawal paused DARIC33-stimulated T cell effector functions, which were restored following reexposure to rapamycin, demonstrating reversible effector function control. While rapamycin-regulated DARIC33 T cells were highly sensitive to target antigen, CD34+ stem cell colony-forming capacity was not impacted. We benchmarked DARIC33 potency relative to CD19 CAR T cells to estimate a T cell dose for clinical testing. In addition, we integrated in vitro and preclinical in vivo drug concentration thresholds for off-on state transitions, as well as murine and human rapamycin pharmacokinetics, to estimate a clinically applicable rapamycin dosing schedule. A phase I DARIC33 trial has been initiated (PLAT-08, NCT05105152), with initial evidence of rapamycin-regulated T cell activation and antitumor impact. Our findings provide evidence that the DARIC platform exhibits sensitive regulation and potency needed for clinical application to other important immunotherapy targets.

High-resolution mapping of the neutralizing and binding specificities of polyclonal sera post-HIV Env trimer vaccination.

Mapping polyclonal serum responses is critical to rational vaccine design. However, most high-resolution mapping approaches involve isolating and characterizing individual antibodies, which incompletely defines the polyclonal response. Here we use two complementary approaches to directly map the specificities of the neutralizing and binding antibodies of polyclonal anti-HIV-1 sera from rabbits immunized with BG505 Env SOSIP trimers. We used mutational antigenic profiling to determine how all mutations in Env affected viral neutralization and electron microscopy polyclonal epitope mapping (EMPEM) to directly visualize serum Fabs bound to Env trimers. The dominant neutralizing specificities were generally only a subset of the more diverse binding specificities. Additional differences between binding and neutralization reflected antigenicity differences between virus and soluble Env trimer. Furthermore, we refined residue-level epitope specificity directly from sera, revealing subtle differences across sera. Together, mutational antigenic profiling and EMPEM yield a holistic view of the binding and neutralizing specificity of polyclonal sera.

Vaccines work by stimulating the immune system to produce proteins called antibodies. These antibodies bind to the virus targeted by the vaccine and block the virus from infecting cells. It has been difficult to develop a vaccine for HIV because frequent mutations allow it to evade antibodies. Understanding exactly how these proteins bind to HIV and how various mutations enable the virus to escape them is crucial to designing a successful HIV vaccine. Over the last decade, scientists have developed new techniques for studying individual antibodies and how they bind to viruses. Now, they are using these insights to design vaccines. Most vaccines result in the production of many antibodies that bind to different parts of the virus, making it harder for a virus to escape. But studying many antibodies with different targets on the virus simultaneously remains challenging. By combining two-cutting edge approaches, Dingens et al. catalogued the many antibodies that rabbits produce in response to an experimental vaccine for HIV. In the experiments, they mapped how two types of rabbit antibodies target the virus: those that could bind to the virus, and those that could both bind and neutralize the virus (i.e., block it from infecting cells). The experiments showed that small differences between the HIV virus and the vaccine explained why some rabbit antibodies created in response to the vaccine could bind but not neutralize the virus. Moreover, the ability to stop HIV from infecting the cells appeared to be reserved to antibodies that could bind to several different locations at the virus. Dingens et al. further documented all the virus mutations that would allow it to evade neutralizing antibodies. The techniques used in the experiments may help scientists identify the best sites on the HIV virus to target with vaccines and to better understand the binding and neutralizing activity of antibodies. The results of the experiments may also help to redesign the experimental HIV vaccine ­ which is currently being tested in humans ­ to be even more effective.

Complete Mapping of Mutations to the SARS-CoV-2 Spike Receptor-Binding Domain that Escape Antibody Recognition.

Antibodies targeting the SARS-CoV-2 spike receptor-binding domain (RBD) are being developed as therapeutics and are a major contributor to neutralizing antibody responses elicited by infection. Here, we describe a deep mutational scanning method to map how all amino-acid mutations in the RBD affect antibody binding and apply this method to 10 human monoclonal antibodies. The escape mutations cluster on several surfaces of the RBD that broadly correspond to structurally defined antibody epitopes. However, even antibodies targeting the same surface often have distinct escape mutations. The complete escape maps predict which mutations are selected during viral growth in the presence of single antibodies. They further enable the design of escape-resistant antibody cocktails-including cocktails of antibodies that compete for binding to the same RBD surface but have different escape mutations. Therefore, complete escape-mutation maps enable rational design of antibody therapeutics and assessment of the antigenic consequences of viral evolution.

Hemiptera phylogenomic resources: Tree-based orthology prediction and conserved exon identification.

High-throughput sequencing of transcriptomes and targeted genomic regions are advancing our knowledge of The Tree of Life. Building phylogenies with regions of the genome requires 1-to-1 orthologue resources of genes and noncoding loci. One organismal group that has received little attention in this area is the Hemiptera, the fifth largest insect order represented by ~103,590 named species. Here, we present a set of 3,872 Hemiptera 1-to-1 orthogroups based on tree-based orthology inference of eight Hemiptera species with publicly available genome sequences. We also estimate a set of 406 orthologous exons with similar mRNA splice sites that can be used for Sanger sequencing and develop enrichment probes for targeted genome sequencing for phylogenomic inference. We show this novel set of orthologues is informative at the protein, coding sequence and exon molecular levels and provides robust branch support in both gene tree-species tree methods and concatenated sequence phylogenies. In addition, we demonstrate the utility of these loci to resolve relationships in whiteflies, Bemisia tabaci, a large species complex with few phylogenomic resources. Last, we compare our Hemiptera phylogeny with previously published phylogenies and other orthologue databases, while providing suggestions on further improvement to this phylogenomic resource.

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.

Complete mapping of mutations to the SARS-CoV-2 spike receptor-binding domain that escape antibody recognition.

Antibodies targeting the SARS-CoV-2 spike receptor-binding domain (RBD) are being developed as therapeutics and make a major contribution to the neutralizing antibody response elicited by infection. Here, we describe a deep mutational scanning method to map how all amino-acid mutations in the RBD affect antibody binding, and apply this method to 10 human monoclonal antibodies. The escape mutations cluster on several surfaces of the RBD that broadly correspond to structurally defined antibody epitopes. However, even antibodies targeting the same RBD surface often have distinct escape mutations. The complete escape maps predict which mutations are selected during viral growth in the presence of single antibodies, and enable us to design escape-resistant antibody cocktails-including cocktails of antibodies that compete for binding to the same surface of the RBD but have different escape mutations. Therefore, complete escape-mutation maps enable rational design of antibody therapeutics and assessment of the antigenic consequences of viral evolution.

Modeling site-specific amino-acid preferences deepens phylogenetic estimates of viral sequence divergence.

Molecular phylogenetics is often used to estimate the time since the divergence of modern gene sequences. For highly diverged sequences, such phylogenetic techniques sometimes estimate surprisingly recent divergence times. In the case of viruses, independent evidence indicates that the estimates of deep divergence times from molecular phylogenetics are sometimes too recent. This discrepancy is caused in part by inadequate models of purifying selection leading to branch-length underestimation. Here we examine the effect on branch-length estimation of using models that incorporate experimental measurements of purifying selection. We find that models informed by experimentally measured site-specific amino-acid preferences estimate longer deep branches on phylogenies of influenza virus hemagglutinin. This lengthening of branches is due to more realistic stationary states of the models, and is mostly independent of the branch-length extension from modeling site-to-site variation in amino-acid substitution rate. The branch-length extension from experimentally informed site-specific models is similar to that achieved by other approaches that allow the stationary state to vary across sites. However, the improvements from all of these site-specific but time homogeneous and site independent models are limited by the fact that a protein's amino-acid preferences gradually shift as it evolves. Overall, our work underscores the importance of modeling site-specific amino-acid preferences when estimating deep divergence times-but also shows the inherent limitations of approaches that fail to account for how these preferences shift over time.

Mapping mutational effects along the evolutionary landscape of HIV envelope.

The immediate evolutionary space accessible to HIV is largely determined by how single amino acid mutations affect fitness. These mutational effects can shift as the virus evolves. However, the prevalence of such shifts in mutational effects remains unclear. Here, we quantify the effects on viral growth of all amino acid mutations to two HIV envelope (Env) proteins that differ at [Formula: see text]100 residues. Most mutations similarly affect both Envs, but the amino acid preferences of a minority of sites have clearly shifted. These shifted sites usually prefer a specific amino acid in one Env, but tolerate many amino acids in the other. Surprisingly, shifts are only slightly enriched at sites that have substituted between the Envs-and many occur at residues that do not even contact substitutions. Therefore, long-range epistasis can unpredictably shift Env's mutational tolerance during HIV evolution, although the amino acid preferences of most sites are conserved between moderately diverged viral strains.

phydms: software for phylogenetic analyses informed by deep mutational scanning.

It has recently become possible to experimentally measure the effects of all amino-acid point mutations to proteins using deep mutational scanning. These experimental measurements can inform site-specific phylogenetic substitution models of gene evolution in nature. Here we describe software that efficiently performs analyses with such substitution models. This software, phydms, can be used to compare the results of deep mutational scanning experiments to the selection on genes in nature. Given a phylogenetic tree topology inferred with another program, phydms enables rigorous comparison of how well different experiments on the same gene capture actual natural selection. It also enables re-scaling of deep mutational scanning data to account for differences in the stringency of selection in the lab and nature. Finally, phydms can identify sites that are evolving differently in nature than expected from experiments in the lab. As data from deep mutational scanning experiments become increasingly widespread, phydms will facilitate quantitative comparison of the experimental results to the actual selection pressures shaping evolution in nature.

Metataxonomic and Metagenomic Approaches vs. Culture-Based Techniques for Clinical Pathology.

Diagnoses that are both timely and accurate are critically important for patients with life-threatening or drug resistant infections. Technological improvements in High-Throughput Sequencing (HTS) have led to its use in pathogen detection and its application in clinical diagnoses of infectious diseases. The present study compares two HTS methods, 16S rRNA marker gene sequencing (metataxonomics) and whole metagenomic shotgun sequencing (metagenomics), in their respective abilities to match the same diagnosis as traditional culture methods (culture inference) for patients with ventilator associated pneumonia (VAP). The metagenomic analysis was able to produce the same diagnosis as culture methods at the species-level for five of the six samples, while the metataxonomic analysis was only able to produce results with the same species-level identification as culture for two of the six samples. These results indicate that metagenomic analyses have the accuracy needed for a clinical diagnostic tool, but full integration in diagnostic protocols is contingent on technological improvements to decrease turnaround time and lower costs.