Globally defining the effects of mutations in a picornavirus capsid

Mapping mutational fitness effects across the coxsackievirus B3 proteome reveals distinct profiles of mutation tolerability

RNA viruses have notoriously high mutation rates due to error-prone replication by their RNA polymerase. However, natural selection concentrates variability in a few key viral proteins. To test whether this stems from different mutation tolerance profiles among viral proteins, we measured the effect of >40,000 non-synonymous mutations across the full proteome of coxsackievirus B3 as well as >97% of all possible codon deletions in the nonstructural proteins. We find significant variation in mutational tolerance within and between individual viral proteins, which correlated with both general and protein-specific structural and functional attributes. Furthermore, mutational fitness effects remained stable across cell lines, suggesting selection pressures are mostly conserved across environments. In addition to providing a rich dataset for understanding virus biology and evolution, our results illustrate that incorporation of mutational tolerance data into druggable pocket discovery can aid in selecting targets with high barriers to drug resistance.

Within-Host Rhinovirus Evolution in Upper and Lower Respiratory Tract Highlights Capsid Variability and Mutation-Independent Compartmentalization

BACKGROUND

Rhinovirus (RV) infections can progress from the upper (URT) to lower (LRT) respiratory tract in immunocompromised individuals, causing high rates of fatal pneumonia. Little is known about how RV evolves within hosts during infection.

METHODS

We sequenced RV complete genomes from 12 hematopoietic cell transplant patients with infection for up to 190 days from both URT (nasal wash, NW) and LRT (bronchoalveolar lavage, BAL). Metagenomic and amplicon next-generation sequencing were used to track the emergence and evolution of intrahost single nucleotide variants (iSNVs).

RESULTS

Identical RV intrahost populations in matched NW and BAL specimens indicated no genetic adaptation is required for RV to progress from URT to LRT. Coding iSNVs were 2.3-fold more prevalent in capsid over nonstructural genes. iSNVs modeled were significantly more likely to be found in capsid surface residues, but were not preferentially located in known RV-neutralizing antibody epitopes. Newly emergent, genotype-matched iSNV haplotypes from immunocompromised individuals in 2008-2010 could be detected in Seattle-area community RV sequences in 2020-2021.

CONCLUSIONS

RV infections in immunocompromised hosts can progress from URT to LRT with no specific evolutionary requirement. Capsid proteins carry the highest variability and emergent mutations can be detected in other, including future, RV sequences.

Uncovering Structural Plasticity of Enterovirus A through Deep Insertional and Deletional Scanning

Insertions and deletions (InDels) are essential sources of novelty in protein evolution. In RNA viruses, InDels cause dramatic phenotypic changes contributing to the emergence of viruses with altered immune profiles and host engagement. This work aimed to expand our current understanding of viral evolution and explore the mutational tolerance of RNA viruses to InDels, focusing on Enterovirus A71 (EV-A71) as a prototype for Enterovirus A species (EV-A). Using newly described deep InDel scanning approaches, we engineered approximately 45,000 insertions and 6,000 deletions at every site across the viral proteome, quantifying their effects on viral fitness. As a general trend, most InDels were lethal to the virus. However, our screen reproducibly identified a set of InDel-tolerant regions, demonstrating our ability to comprehensively map tolerance to these mutations. Tolerant sites highlighted structurally flexible and mutationally plastic regions of viral proteins that avoid core structural and functional elements. Phylogenetic analysis on EV-A species infecting diverse mammalian hosts revealed that the experimentally-identified hotspots overlapped with sites of InDels across the EV-A species, suggesting structural plasticity at these sites is an important function for InDels in EV speciation. Our work reveals the fitness effects of InDels across EV-A71, identifying regions of evolutionary capacity that require further monitoring, which could guide the development of Enterovirus vaccines.

Comprehensive profiling of neutralizing polyclonal sera targeting coxsackievirus B3

Despite their fundamental role in resolving viral infections, our understanding of how polyclonal neutralizing antibody responses target non-enveloped viruses remains limited. To define these responses, we obtained the full antigenic profile of multiple human and mouse polyclonal sera targeting the capsid of a prototypical picornavirus, coxsackievirus B3. Our results uncover significant variation in the breadth and strength of neutralization sites targeted by individual human polyclonal responses, which contrasted with homogenous responses observed in experimentally infected mice. We further use these comprehensive antigenic profiles to define key structural and evolutionary parameters that are predictive of escape, assess epitope dominance at the population level, and reveal a need for at least two mutations to achieve significant escape from multiple sera. Overall, our data provide a comprehensive analysis of how polyclonal sera target a non-enveloped viral capsid and help define both immune dominance and escape at the population level.

Within-host rhinovirus evolution in upper and lower respiratory tract highlights capsid variability and mutation-independent compartmentalization

Background

Human rhinovirus (HRV) infections can progress from the upper (URT) to lower (LRT) respiratory tract in immunocompromised individuals, causing high rates of fatal pneumonia. Little is known about how HRV evolves within hosts during infection.

Methods

We sequenced HRV complete genomes from 12 hematopoietic cell transplant patients with prolonged infection for up to 190 days from both URT (nasal wash, NW) and LRT (bronchoalveolar lavage, BAL) specimens. Metagenomic (mNGS) and amplicon-based NGS were used to study the emergence and evolution of intra-host single nucleotide variants (iSNVs).

Results

Identical HRV intra-host populations in matched NW and BAL specimens indicated no genetic adaptation is required for HRV to progress from URT to LRT. Microbial composition between matched NW and BAL confirmed no cross-contamination during sampling procedure. Coding iSNVs were 2.3-fold more prevalent in capsid over non-structural genes, adjusted for length. iSNVs modeled onto HRV capsid structures were significantly more likely to be found in surface residues, but were not preferentially located in known HRV neutralizing antibody epitopes. Newly emergent, serotype-matched iSNV haplotypes from immunocompromised individuals from 2008-2010 could be detected in Seattle-area community HRV sequences from 2020-2021.

Conclusion

HRV infections in immunocompromised hosts can progress from URT to LRT with no specific evolutionary requirement. Capsid proteins carry the highest variability and emergent mutations can be detected in other, including future, HRV sequences.

Polymorphisms predicting phylogeny in hepatitis B virus

Hepatitis B viruses (HBVs) are compact viruses with circular genomes of ∼3.2 kb in length. Four genes (HBx, Core, Surface, and Polymerase) generating seven products are encoded on overlapping reading frames. Ten HBV genotypes have been characterised (A-J), which may account for differences in transmission, outcomes of infection, and treatment response. However, HBV genotyping is rarely undertaken, and sequencing remains inaccessible in many settings. We set out to assess which amino acid (aa) sites in the HBV genome are most informative for determining genotype, using a machine learning approach based on random forest algorithms (RFA). We downloaded 5,496 genome-length HBV sequences from a public database, excluding recombinant sequences, regions with conserved indels, and genotypes I and J. Each gene was separately translated into aa, and the proteins concatenated into a single sequence (length 1,614 aa). Using RFA, we searched for aa sites predictive of genotype and assessed covariation among the sites with a mutual information-based method. We were able to discriminate confidently between genotypes A-H using ten aa sites. Half of these sites (5/10) sites were identified in Polymerase (Pol), of which 4/5 were in the spacer domain and one in reverse transcriptase. A further 4/10 sites were located in Surface protein and a single site in HBx. There were no informative sites in Core. Properties of the aa were generally not conserved between genotypes at informative sites. Among the highest co-varying pairs of sites, there were fifty-five pairs that included one of these 'top ten' sites. Overall, we have shown that RFA analysis is a powerful tool for identifying aa sites that predict the HBV lineage, with an unexpectedly high number of such sites in the spacer domain, which has conventionally been viewed as unimportant for structure or function. Our results improve ease of genotype prediction from limited regions of HBV sequences and may have future applications in understanding HBV evolution.

Identification of Cytoplasmic Chaperone Networks Relevant for Respiratory Syncytial Virus Replication

RNA viruses have limited coding capacity and must therefore successfully subvert cellular processes to facilitate their replication. A fundamental challenge faced by both viruses and their hosts is the ability to achieve the correct folding and assembly of their proteome while avoiding misfolding and aggregation. In cells, this process is facilitated by numerous chaperone systems together with a large number of co-chaperones. In this work, we set out to define the chaperones and co-chaperones involved in the replication of respiratory syncytial virus (RSV). Using an RNAi screen, we identify multiple members of cellular protein folding networks whose knockdown alters RSV replication. The reduced number of chaperones and co-chaperones identified in this work can facilitate the unmasking of specific chaperone subnetworks required for distinct steps of the RSV life cycle and identifies new potential targets for antiviral therapy. Indeed, we show that the pharmacological inhibition of one of the genes identified in the RNAi screen, valosin-containing protein (VCP/p97), can impede the replication of RSV by interfering with the infection cycle at multiple steps.

Top Three Strategies of ss(+)RNA Plant Viruses: Great Opportunists and Ecosystem Tuners with a Small Genome

ss(+)RNA viruses represent the dominant group of plant viruses. They owe their evolutionary superiority to the large number of mutations that occur during replication, courtesy of RNA-dependent RNA polymerase. Natural selection rewards successful viral subtypes, whose effective tuning of the ecosystem regulates the interactions between its participants. Thus, ss(+)RNA viruses act as shuttles for the functionally important genes of the participants in symbiotic relationships within the ecosystem, of which the most common ecological triad is “plant–virus–insect”. Due to their short life cycle and large number of offspring, RNA viruses act as skillful tuners of the ecosystem, which benefits both viruses and the system as a whole. A fundamental understanding of this aspect of the role played by viruses in the ecosystem makes it possible to apply this knowledge to the creation of DNA insecticides. In fact, since the genes that viruses are involved in transferring are functionally important for both insects and plants, silencing these genes (for example, in insects) can be used to regulate the pest population. RNA viruses are increasingly treated not as micropathogens but as necessary regulators of ecosystem balance.

Increased RNA virus population diversity improves adaptability

The replication machinery of most RNA viruses lacks proofreading mechanisms. As a result, RNA virus populations harbor a large amount of genetic diversity that confers them the ability to rapidly adapt to changes in their environment. In this work, we investigate whether further increasing the initial population diversity of a model RNA virus can improve adaptation to a single selection pressure, thermal inactivation. For this, we experimentally increased the diversity of coxsackievirus B3 (CVB3) populations across the capsid region. We then compared the ability of these high diversity CVB3 populations to achieve resistance to thermal inactivation relative to standard CVB3 populations in an experimental evolution setting. We find that viral populations with high diversity are better able to achieve resistance to thermal inactivation at both the temperature employed during experimental evolution as well as at a more extreme temperature. Moreover, we identify mutations in the CVB3 capsid that confer resistance to thermal inactivation, finding significant mutational epistasis. Our results indicate that even naturally diverse RNA virus populations can benefit from experimental augmentation of population diversity for optimal adaptation and support the use of such viral populations in directed evolution efforts that aim to select viruses with desired characteristics.