Viral Fitness Research Articles

Enhanced Diversifying Selection on Polymerase Genes in H5N1 Clade 2.3.4.4b: A Key Driver of Altered Species Tropism and Host Range Expansion.

Highly pathogenic avian influenza H5N1 clade 2.3.4.4b viruses have shown unprecedented host range and pathogenicity, including infections in cattle, previously not susceptible to H5N1. We investigated whether selection pressures on clade 2.3.4.4b viral genes could shed light on their unique epidemiological features. Our analysis revealed that while the gene products of clade 2.3.4.4b H5N1 primarily undergo purifying selection, there are notable instances of episodic diversifying selection. Specifically, the polymerase genes PB2, PB1, and PA exhibit significantly greater selection pressures in clade 2.3.4.4b than all earlier H5N1 virus clades. Polymerases play critical roles in influenza virus adaptation, including viral fitness, interspecies transmission, and virulence. Our findings provide evidence that significant selection pressures have shaped the evolution of the H5N1 clade 2.3.4.4b viruses, facilitating their expanded host tropism and the potential for further adaptation to mammalian hosts. We discuss how exogenous factors, such as altered bird migration patterns and increased host susceptibility, may have contributed to the expanded host range. As H5N1 viruses continue to infect new hosts, there is a greater risk of emergent novel variants with increased pathogenicity in humans and animals. Thus, comprehensive One Health surveillance is critical to monitor transmission among avian and mammalian hosts.

Dynamic clade transitions and the influence of vaccination on the spatiotemporal circulation of SARS-CoV-2 variants

Since 2021, the emergence of variants of concern (VOC) has led Brazil to experience record numbers of in COVID-19 cases and deaths. The expanded spread of the SARS-CoV-2 combined with a low vaccination rate has contributed to the emergence of new mutations that may enhance viral fitness, leading to the persistence of the disease. Due to limitations in the real-time genomic monitoring of new variants in some Brazilian states, we aimed to investigate whether genomic surveillance, coupled with epidemiological data and SARS-CoV-2 variants spatiotemporal spread in a smaller region, can reflect the pandemic progression at a national level. Our findings revealed three SARS-CoV-2 variant replacements from 2021 to early 2022, corresponding to the introduction and increase in the frequency of Gamma, Delta, and Omicron variants, as indicated by peaks of the Effective Reproductive Number (Reff). These distinct clade replacements triggered two waves of COVID-19 cases, influenced by the increasing vaccine uptake over time. Our results indicated that the effectiveness of vaccination in preventing new cases during the Delta and Omicron circulations was six and eleven times higher, respectively, than during the period when Gamma was predominant, and it was highly efficient in reducing the number of deaths. Furthermore, we demonstrated that genomic monitoring at a local level can reflect the national trends in the spread and evolution of SARS-CoV-2.

Phylogenetic Network Analyses Reveal the Influence of Transmission Clustering on the Spread of HIV Drug Resistance in Quebec from 2002 to 2022.

HIV drug resistance (HIV-DR) may jeopardize the benefit of antiretroviral therapy (ART) in treatment and prevention. This study utilized viral phylogenetics to resolve the influence of transmission networks on sustaining the spread of HIV-DR in Quebec spanning 2002 to 2022. Time trends in acquired (ADR) and transmitted drug resistance (TDR) were delineated in treatment-experienced (n = 3500) and ART-naïve persons (n = 6011) with subtype B infections. Similarly, non-B-subtype HIV-DR networks were assessed pre- (n = 1577) and post-ART experience (n = 488). Risks of acquisition of resistance-associated mutations (RAMs) were related to clustering using 1, 2-5, vs. 6+ members per cluster as categorical variables. Despite steady declines in treatment failure and ADR since 2007, rates of TDR among newly infected, ART-naive persons remained at 14% spanning the 2007-2011, 2012-2016, and 2017-2022 periods. Notably, half of new infections among men having sex with men and heterosexual groups were linked in large, clustered networks having a median of 35 (14-73 IQR) and 16 (9-26 IQR) members per cluster, respectively. Cluster membership and size were implicated in forward transmission of non-nucleoside reverse transcriptase inhibitor NNRTI RAMs (9%) and thymidine analogue mutations (TAMs) (5%). In contrast, transmission of M184V, K65R, and integrase inhibitors (1-2%) remained rare. Levels of TDR reflected viral replicative fitness. The median baseline viremia in ART-naïve groups having no RAMs, NNRTI RAMs, TAMs, and M184VI were 46.088, 38,447, 20,330, and 6811 copies/mL, respectively (p < 0.0001). Phylogenetics emphasize the need to prioritize ART and pre-exposure prophylaxis strategies to avert the expansion of transmission cascades of HIV-DR.

Quantifying the strength of viral fitness trade-offs between hosts: a meta-analysis of pleiotropic fitness effects

Abstract The range of hosts a given virus can infect is widely presumed to be limited by fitness trade-offs between alternative hosts. These fitness trade-offs may arise naturally due to antagonistic pleiotropy if mutations that increase fitness in one host tend to decrease fitness in alternate hosts. Yet there is also growing recognition that positive pleiotropy may be more common than previously appreciated. With positive pleiotropy, mutations have concordant fitness effects such that a beneficial mutation can simultaneously increase fitness in different hosts, providing a genetic mechanism by which selection can overcome fitness trade-offs. How readily evolution can overcome fitness trade-offs therefore depends on the overall distribution of mutational fitness effects between hosts, including the relative frequency of antagonistic versus positive pleiotropy. We therefore conducted a systematic meta-analysis of the pleiotropic fitness effects of viral mutations reported in different hosts. Our analysis indicates that while both antagonistic and positive pleiotropy are common, fitness effects are overall positively correlated between hosts and unconditionally beneficial mutations are not uncommon. Moreover, the relative frequency of antagonistic versus positive pleiotropy may simply reflect the underlying frequency of beneficial and deleterious mutations in individual hosts. Given a mutation is beneficial in one host, the probability that it is deleterious in another host is roughly equal to the probability that any mutation is deleterious, suggesting there is no natural tendency toward antagonistic pleiotropy. The widespread prevalence of positive pleiotropy suggests that many fitness trade-offs may be readily overcome by evolution given the right selection pressures.

Molecular basis of RNA recombination in the 3'UTR of chikungunya virus genome.

Chikungunya virus (CHIKV) is a rapidly spreading re-emergent virus transmitted from mosquitoes to humans. The emergence of epidemic variants has been associated with changes in the viral genome, such as the duplication of repeated sequences in the 3' untranslated region (UTR). Indeed, blocks of repeated sequences seemingly favor RNA recombination, providing the virus with a unique ability to continuously change the 3'UTR architecture during host switching. In this work, we provide experimental data on the molecular mechanism of RNA recombination and describe specific sequence and structural elements in the viral 3'UTR that favor template switching of the viral RNA-dependent RNA polymerase on the 3'UTR. Furthermore, we found that a 3'UTR deletion mutant that exhibits markedly delayed replication in mosquito cells and impaired transmission in vivo, recombines in reference laboratory strains of mosquitoes. Altogether, our data provide novel experimental evidence indicating that RNA recombination can act as a nucleic acid repair mechanism to add repeated sequences that are associated to high viral fitness in mosquito during chikungunya virus replication.

The interferon-regulated host factor hnRNPA0 modulates HIV-1 production by interference with LTR activity, mRNA trafficking, and programmed ribosomal frameshifting.

The interplay between host factors and viral components impacts viral replication efficiency profoundly. Members of the cellular heterogeneous nuclear ribonucleoprotein family (hnRNPs) have been extensively studied as HIV-1 host dependency factors, but whether they play a role in innate immunity is currently unknown. This study aimed to identify hnRNPA0 as a type I interferon (IFN)-repressed host factor in HIV-1-infected cells. Knockdown of hnRNPA0, a situation that mirrors conditions under IFN stimulation, increased LTR activity, export of unspliced HIV-1 mRNA, viral particle production, and thus, increased infectivity. Conversely, hnRNPA0 overexpression primarily reduced plasmid-driven and integrated HIV-1 long terminal repeat (LTR) activity, significantly decreasing total viral mRNA and protein levels. In addition, high levels of hnRNPA0 significantly reduced the HIV-1 programmed ribosomal frameshifting efficiency, resulting in a shift in the HIV-1 p55/p15 ratio. The HIV-1 alternative splice site usage remained largely unaffected by altered hnRNPA0 levels suggesting that the synergistic inhibition of the LTR activity and viral mRNA transcription, as well as impaired ribosomal frameshifting efficiency, are critical factors for efficient HIV-1 replication regulated by hnRNPA0. The pleiotropic dose-dependent effects under high or low hnRNPA0 levels were further confirmed in HIV-1-infected Jurkat cells. Finally, our study revealed that hnRNPA0 levels in PBMCs were lower in therapy-naive HIV-1-infected individuals compared to healthy controls. Our findings highlight a significant role for hnRNPA0 in HIV-1 replication and suggest that its IFN-I-regulated expression levels are critical for viral fitness allowing replication in an antiviral environment.IMPORTANCERNA-binding proteins, in particular, heterogeneous nuclear ribonucleoproteins (hnRNPs), have been extensively studied. Some act as host dependency factors for HIV-1 since they are involved in multiple cellular gene expression processes. Our study revealed hnRNPA0 as an IFN-regulated host factor, that is differently expressed after IFN-I treatment in HIV-1 target cells and lower expressed in therapy-naïve HIV-1-infected individuals. Our findings demonstrate the significant pleiotropic role of hnRNPA0 in viral replication: In high concentrations, hnRNPA0 limits viral replication by negatively regulating Tat-LTR transcription, retaining unspliced mRNA in the nucleus, and significantly impairing programmed ribosomal frameshifting. Low hnRNPA0 levels as observed in IFN-treated THP-1 cells, particularly facilitate HIV LTR activity and unspliced mRNA export, suggesting a role in innate immunity in favor of HIV replication. Understanding the mode of action between hnRNPA0 and HIV-1 gene expression might help to identify novel therapeutically strategies against HIV-1 and other viruses.

Computational Analysis of the Accumulation of Mutations in Therapeutically Important RNA Viral Proteins During Pandemics with Special Emphasis on SARS-CoV-2

Single stranded RNA viruses are primary causative agents for pandemics, causing extensive morbidity and mortality worldwide. A pivotal question in pandemic preparedness and therapeutic intervention is what are the specific mutations which are more likely to emerge during such global health crises? This study aims to identify markers for mutations with the highest probability of emergence in these pandemics, focusing on the SARS-CoV-2 spike protein, an essential and therapeutically significant viral protein, starting from sequence information from the onset of the pandemic until July 2022. Quite consistently, we observed that emerged mutations tended to demonstrate a high genetic score, which reflects high similarity of the type of codon required for translation between an amino acid and to the mutated one. Further, this pattern is also observed in therapeutically significant proteins of other ssRNA pandemic viruses, including influenza (HA, NA), spike proteins of Ebola, envelope of Dengue and Chikungunya. We propose that the genetic score serves as an initial indicator, preceding the actual impact of the mutation on viral fitness. Finally, we developed a comprehensive computational pipeline to further explore and predict the subsequent effects of mutations on viral fitness. We believe that our pipeline can narrow down and predict future mutations in therapeutically important viral proteins during a pandemic.

Dengue virus preferentially uses human and mosquito non-optimal codons

Codon optimality refers to the effect that codon composition has on messenger RNA (mRNA) stability and translation level and implies that synonymous codons are not silent from a regulatory point of view. Here, we investigated the adaptation of virus genomes to the host optimality code using mosquito-borne dengue virus (DENV) as a model. We demonstrated that codon optimality exists in mosquito cells and showed that DENV preferentially uses nonoptimal (destabilizing) codons and avoids codons that are defined as optimal (stabilizing) in either human or mosquito cells. Human genes enriched in the codons preferentially and frequently used by DENV are upregulated during infection, and so is the tRNA decoding the nonoptimal and DENV preferentially used codon for arginine. We found that adaptation during single-host passaging in human or mosquito cells results in the selection of synonymous mutations towards DENV’s preferred nonoptimal codons that increase virus fitness. Finally, our analyses revealed that hundreds of viruses preferentially use nonoptimal codons, with those infecting a single host displaying an even stronger bias, suggesting that host–pathogen interaction shapes virus-synonymous codon choice.

Identification of distinct genotypes in circulating RSV A strains based on variants in the virus replication-associated genes.

Respiratory syncytial virus (RSV) is a common cause of respiratory infection that often leads to hospitalization of infected younger children and older adults. RSV is classified into two strains, A and B, each with several subgroups or genotypes. One issue with the definition of these subgroups is the lack of a unified method of identification or genotyping. We propose that genotyping strategies based on the genes coding for replication-associated proteins could provide critical information on the replication capacity of the distinct subgroups, while clearly distinguishing genotypes. Here, we analyzed the virus replication-associated genes N, P, M2, and L from de novo assembled RSV A sequences obtained from 31 newly sequenced samples from hospitalized patients in Philadelphia and 78 additional publicly available sequences from different geographic locations within the United States. In-depth analysis and annotation of variants in the replication-associated proteins identified the polymerase protein L as a robust target for genotyping RSV subgroups. Importantly, our analysis revealed non-synonymous variations in L that were consistently accompanied by conserved changes in its co-factor P or the M2-2 protein, suggesting associations and interactions between specific domains of these proteins. Similar associations were seen among sequences of the related human metapneumovirus. These results highlight L as an alternative to other RSV genotyping targets and demonstrate the value of in-depth analyses and annotations of RSV sequences as it can serve as a foundation for subsequent in vitro and clinical studies on the efficiency of the polymerase and fitness of different virus isolates.IMPORTANCEGiven the historical heterogeneity of respiratory syncytial virus (RSV) and the disease it causes, there is a need to understand the properties of the circulating RSV strains each season. This information would benefit from an informative and consensus method of genotyping the virus. Here, we carried out a variant analysis that shows a pattern of specific variations among the replication-associated genes of RSV A across different seasons. Interestingly, these variation patterns, which were also seen in human metapneumovirus sequences, point to previously defined interactions of domains within these genes, suggesting co-variation in the replication-associated genes. Our results also suggest a genotyping strategy that can prove to be particularly important in understanding the genotype-phenotype correlation in the era of RSV vaccination, where selective pressure on the virus to evolve is anticipated. More importantly, the categorization of pneumoviruses based on these patterns may be of prognostic value.

High throughput profiling identified PA-L106R amino acid substitution in A(H1N1)pdm09 influenza virus that confers reduced susceptibility to baloxavir in vitro

Baloxavir acid (BXA) is a pan-influenza antiviral that targets the cap-dependent endonuclease of the polymerase acidic (PA) protein required for viral mRNA synthesis. To gain a comprehensive understanding on the molecular changes associated with reduced susceptibility to BXA and their fitness profile, we performed a deep mutational scanning at the PA endonuclease domain of an A (H1N1)pdm09 virus. The recombinant virus libraries were serially passaged in vitro under increasing concentrations of BXA followed by next-generation sequencing to monitor PA amino acid substitutions with increased detection frequencies. Enriched PA amino acid changes were each introduced into a recombinant A (H1N1)pdm09 virus to validate their effect on BXA susceptibility and viral replication fitness in vitro. The I38 T/M substitutions known to confer reduced susceptibility to BXA were invariably detected from recombinant virus libraries within 5 serial passages. In addition, we identified a novel L106R substitution that emerged in the third passage and conferred greater than 10-fold reduced susceptibility to BXA. PA-L106 is highly conserved among seasonal influenza A and B viruses. Compared to the wild-type virus, the L106R substitution resulted in reduced polymerase activity and a minor reduction of the peak viral load, suggesting the amino acid change may result in moderate fitness loss. Our results support the use of deep mutational scanning as a practical tool to elucidate genotype-phenotype relationships, including mapping amino acid substitutions with reduced susceptibility to antivirals.

Atomistic Prediction of Structures, Conformational Ensembles and Binding Energetics for the SARS-CoV-2 Spike JN.1, KP.2 and KP.3 Variants Using AlphaFold2 and Molecular Dynamics Simulations: Mutational Profiling and Binding Free Energy Analysis Reveal Epistatic Hotspots of the ACE2 Affinity and Immune Escape.

The most recent wave of SARS-CoV-2 Omicron variants descending from BA.2 and BA.2.86 exhibited improved viral growth and fitness due to convergent evolution of functional hotspots. These hotspots operate in tandem to optimize both receptor binding for effective infection and immune evasion efficiency, thereby maintaining overall viral fitness. The lack of molecular details on structure, dynamics and binding energetics of the latest FLiRT and FLuQE variants with the ACE2 receptor and antibodies provides a considerable challenge that is explored in this study. We combined AlphaFold2-based atomistic predictions of structures and conformational ensembles of the SARS-CoV-2 Spike complexes with the host receptor ACE2 for the most dominant Omicron variants JN.1, KP.1, KP.2 and KP.3 to examine the mechanisms underlying the role of convergent evolution hotspots in balancing ACE2 binding and antibody evasion. Using the ensemble-based mutational scanning of the spike protein residues and computations of binding affinities, we identified binding energy hotspots and characterized molecular basis underlying epistatic couplings between convergent mutational hotspots. The results suggested that the existence of epistatic interactions between convergent mutational sites at L455, F456, Q493 positions that enable to protect and restore ACE2 binding affinity while conferring beneficial immune escape. To examine immune escape mechanisms, we performed structure-based mutational profiling of the spike protein binding with several classes of antibodies that displayed impaired neutralization against BA.2.86, JN.1, KP.2 and KP.3. The results confirmed the experimental data that JN.1, KP.2 and KP.3 harboring the L455S and F456L mutations can significantly impair the neutralizing activity of class-1 monoclonal antibodies, while the epistatic effects mediated by F456L can facilitate the subsequent convergence of Q493E changes to rescue ACE2 binding. Structural and energetic analysis provided a rationale to the experimental results showing that BD55-5840 and BD55-5514 antibodies that bind to different binding epitopes can retain neutralizing efficacy against all examined variants BA.2.86, JN.1, KP.2 and KP.3. The results support the notion that evolution of Omicron variants may favor emergence of lineages with beneficial combinations of mutations involving mediators of epistatic couplings that control balance of high ACE2 affinity and immune evasion.

Evolution at Spike protein position 519 in SARS-CoV-2 facilitated adaptation to humans

As the COVID-19 pandemic enters its fourth year, the pursuit of identifying a progenitor virus to SARS-CoV-2 and understanding the mechanism of its emergence persists, albeit against the backdrop of intensified efforts to monitor the ongoing evolution of the virus and the influx of new mutations. Surprisingly, few residues hypothesized to be essential for SARS-CoV-2 emergence and adaptation to humans have been validated experimentally, despite the importance that these mutations could contribute to the development of effective antivirals. To remedy this, we searched for genomic regions in the SARS-CoV-2 genome that show evidence of past selection around residues unique to SARS-CoV-2 compared with closely related coronaviruses. In doing so, we identified a residue at position 519 in Spike within the receptor binding domain that holds a static histidine in human-derived SARS-CoV-2 sequences but an asparagine in SARS-related coronaviruses from bats and pangolins. In experimental validation, the SARS-CoV-2 Spike protein mutant carrying the putatively ancestral H519N substitution showed reduced replication in human lung cells, suggesting that the histidine residue contributes to viral fitness in the human host. Structural analyses revealed a potential role of Spike residue 519 in mediating conformational transitions necessary for Spike prior to binding with ACE2. Pseudotyped viruses bearing the putatively ancestral N519 also demonstrated significantly reduced infectivity in cells expressing the human ACE2 receptor compared to H519. ELISA data corroborated that H519 enhances Spike binding affinity to the human ACE2 receptor compared to the putatively ancestral N519. Collectively, these findings suggest that the evolutionary transition at position 519 of the Spike protein played a critical role in SARS-CoV-2 emergence and adaptation to the human host. Additionally, this residue presents as a potential drug target for designing small molecule inhibitors tailored to this site.

Virological characteristics of the SARS-CoV-2 Omicron EG.5.1 variant.

In middle to late 2023, a sublineage of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron XBB, EG.5.1 (a progeny of XBB.1.9.2), is spreading rapidly around the world. We performed multiscale investigations, including phylogenetic analysis, epidemic dynamics modeling, infection experiments using pseudoviruses, clinical isolates, and recombinant viruses in cell cultures and experimental animals, and the use of human sera and antiviral compounds, to reveal the virological features of the newly emerging EG.5.1 variant. Our phylogenetic analysis and epidemic dynamics modeling suggested that two hallmark substitutions of EG.5.1, S:F456L and ORF9b:I5T are critical to its increased viral fitness. Experimental investigations on the growth kinetics, sensitivity to clinically available antivirals, fusogenicity, and pathogenicity of EG.5.1 suggested that the virological features of EG.5.1 are comparable to those of XBB.1.5. However, cryo-electron microscopy revealed structural differences between the spike proteins of EG.5.1 and XBB.1.5. We further assessed the impact of ORF9b:I5T on viral features, but it was almost negligible in our experimental setup. Our multiscale investigations provide knowledge for understanding the evolutionary traits of newly emerging pathogenic viruses, including EG.5.1, in the human population.

Spike deep mutational scanning helps predict success of SARS-CoV-2 clades

SARS-CoV-2 variants acquire mutations in the spike protein that promote immune evasion1 and affect other properties that contribute to viral fitness, such as ACE2 receptor binding and cell entry2,3. Knowledge of how mutations affect these spike phenotypes can provide insight into the current and potential future evolution of the virus. Here we use pseudovirus deep mutational scanning4 to measure how more than 9,000 mutations across the full XBB.1.5 and BA.2 spikes affect ACE2 binding, cell entry or escape from human sera. We find that mutations outside the receptor-binding domain (RBD) have meaningfully affected ACE2 binding during SARS-CoV-2 evolution. We also measure how mutations to the XBB.1.5 spike affect neutralization by serum from individuals who recently had SARS-CoV-2 infections. The strongest serum escape mutations are in the RBD at sites 357, 420, 440, 456 and 473; however, the antigenic effects of these mutations vary across individuals. We also identify strong escape mutations outside the RBD; however, many of them decrease ACE2 binding, suggesting they act by modulating RBD conformation. Notably, the growth rates of human SARS-CoV-2 clades can be explained in substantial part by the measured effects of mutations on spike phenotypes, suggesting our data could enable better prediction of viral evolution.

Temperature-dependent Spike-ACE2 interaction of Omicron subvariants is associated with viral transmission.

The continued evolution of severe acute respiratory syndrome 2 (SARS-CoV-2) requires persistent monitoring of its subvariants. Omicron subvariants are responsible for the vast majority of SARS-CoV-2 infections worldwide, with XBB and BA.2.86 sublineages representing more than 90% of circulating strains as of January 2024. To better understand parameters involved in viral transmission, we characterized the functional properties of Spike glycoproteins from BA.2.75, CH.1.1, DV.7.1, BA.4/5, BQ.1.1, XBB, XBB.1, XBB.1.16, XBB.1.5, FD.1.1, EG.5.1, HK.3, BA.2.86 and JN.1. We tested their capacity to evade plasma-mediated recognition and neutralization, binding to angiotensin-converting enzyme 2 (ACE2), their susceptibility to cold inactivation, Spike processing, as well as the impact of temperature on Spike-ACE2 interaction. We found that compared to the early wild-type (D614G) strain, most Omicron subvariants' Spike glycoproteins evolved to escape recognition and neutralization by plasma from individuals who received a fifth dose of bivalent (BA.1 or BA.4/5) mRNA vaccine and improve ACE2 binding, particularly at low temperatures. Moreover, BA.2.86 had the best affinity for ACE2 at all temperatures tested. We found that Omicron subvariants' Spike processing is associated with their susceptibility to cold inactivation. Intriguingly, we found that Spike-ACE2 binding at low temperature was significantly associated with growth rates of Omicron subvariants in humans. Overall, we report that Spikes from newly emerged Omicron subvariants are relatively more stable and resistant to plasma-mediated neutralization, present improved affinity for ACE2 which is associated, particularly at low temperatures, with their growth rates.IMPORTANCEThe persistent evolution of SARS-CoV-2 gave rise to a wide range of variants harboring new mutations in their Spike glycoproteins. Several factors have been associated with viral transmission and fitness such as plasma-neutralization escape and ACE2 interaction. To better understand whether additional factors could be of importance in SARS-CoV-2 variants' transmission, we characterize the functional properties of Spike glycoproteins from several Omicron subvariants. We found that the Spike glycoprotein of Omicron subvariants presents an improved escape from plasma-mediated recognition and neutralization, Spike processing, and ACE2 binding which was further improved at low temperature. Intriguingly, Spike-ACE2 interaction at low temperature is strongly associated with viral growth rate, as such, low temperatures could represent another parameter affecting viral transmission.

R229I substitution from oseltamivir induction in HA1 region significantly increased the fitness of a H7N9 virus bearing NA 292K

ABSTRACT The proportion of human isolates with reduced neuraminidase inhibitors (NAIs) susceptibility in highly pathogenic avian influenza (HPAI) H7N9 virus was high. These drug-resistant strains showed good replication capacity without serious loss of fitness. In the presence of oseltamivir, R229I substitution were found in HA1 region of the HPAI H7N9 virus before NA R292K appeared. HPAI H7N9 or H7N9/PR8 recombinant viruses were developed to study whether HA R229I could increase the fitness of the H7N9 virus bearing NA 292K. Replication efficiency was assessed in MDCK or A549 cells. Neuraminidase enzyme activity and receptor-binding ability were analyzed. Pathogenicity in C57 mice was evaluated. Antigenicity analysis was conducted through a two-way HI test, in which the antiserum was obtained from immunized ferrets. Transcriptomic analysis of MDCK infected with HPAI H7N9 24hpi was done. It turned out that HA R229I substitution from oseltamivir induction in HA1 region increased (1) replication ability in MDCK(P < 0.05) and A549(P < 0.05), (2) neuraminidase enzyme activity, (3) binding ability to both α2,3 and α2,6 receptor, (4) pathogenicity to mice(more weight loss; shorter mean survival day; viral titer in respiratory tract, P < 0.05; Pathological changes in pneumonia), (5) transcriptome response of MDCK, of the H7N9 virus bearing NA 292K. Besides, HA R229I substitution changed the antigenicity of H7N9/PR8 virus (>4-fold difference of HI titre). It indicated that through the fine-tuning of HA-NA balance, R229I increased the fitness and changed the antigenicity of H7N9 virus bearing NA 292K. Public health attention to this mechanism needs to be drawn.

An efficient plasmid-based system for the recovery of recombinant vesicular stomatitis virus encoding foreign glycoproteins

Viral glycoproteins mediate entry into host cells, thereby dictating host range and pathogenesis. In addition, they constitute the principal target of neutralizing antibody responses, making them important antigens in vaccine development. Recombinant vesicular stomatitis virus (VSV) encoding foreign glycoproteins can provide a convenient and safe surrogate system to interrogate the function, evolution, and antigenicity of viral glycoproteins from viruses that are difficult to manipulate or those requiring high biosafety level containment. However, the production of recombinant VSV can be technically challenging. In this work, we present an efficient and robust plasmid-based system for the production of recombinant VSV encoding foreign glycoproteins. We validate the system using glycoproteins from different viral families, including arenaviruses, coronaviruses, and hantaviruses, as well as highlight their utility for studying the effects of mutations on viral fitness. Overall, the methods described herein can facilitate the study of both native and recombinant VSV encoding foreign glycoproteins and can serve as the basis for the production of VSV-based vaccines.

Natural variation in neuraminidase activity influences the evolutionary potential of the seasonal H1N1 lineage hemagglutinin.

The antigenic evolution of the influenza A virus hemagglutinin (HA) gene poses a major challenge for the development of vaccines capable of eliciting long-term protection. Prior efforts to understand the mechanisms that govern viral antigenic evolution mainly focus on HA in isolation, ignoring the fact that HA must act in concert with the viral neuraminidase (NA) during replication and spread. Numerous studies have demonstrated that the degree to which the receptor-binding avidity of HA and receptor-cleaving activity of NA are balanced with each other influences overall viral fitness. We recently showed that changes in NA activity can significantly alter the mutational fitness landscape of HA in the context of a lab-adapted virus strain. Here, we test whether natural variation in relative NA activity can influence the evolutionary potential of HA in the context of the seasonal H1N1 lineage (pdmH1N1) that has circulated in humans since the 2009 pandemic. We observed substantial variation in the relative activities of NA proteins encoded by a panel of H1N1 vaccine strains isolated between 2009 and 2019. We comprehensively assessed the effect of NA background on the HA mutational fitness landscape in the circulating pdmH1N1 lineage using deep mutational scanning and observed pronounced epistasis between NA and residues in or near the receptor-binding site of HA. To determine whether NA variation could influence the antigenic evolution of HA, we performed neutralizing antibody selection experiments using a panel of monoclonal antibodies targeting different HA epitopes. We found that the specific antibody escape profiles of HA were highly contingent upon NA background. Overall, our results indicate that natural variation in NA activity plays a significant role in governing the evolutionary potential of HA in the currently circulating pdmH1N1 lineage.

Spike N354 glycosylation augments SARS-CoV-2 fitness for human adaptation through structural plasticity.

Selective pressures have given rise to a number of SARS-CoV-2 variants during the prolonged course of the COVID-19 pandemic. Recently evolved variants differ from ancestors in additional glycosylation within the spike protein receptor-binding domain (RBD). Details of how the acquisition of glycosylation impacts viral fitness and human adaptation are not clearly understood. Here, we dissected the role of N354-linked glycosylation, acquired by BA.2.86 sub-lineages, as a RBD conformational control element in attenuating viral infectivity. The reduced infectivity is recovered in the presence of heparin sulfate, which targets the 'N354 pocket' to ease restrictions of conformational transition resulting in a 'RBD-up' state, thereby conferring an adjustable infectivity. Furthermore, N354 glycosylation improved spike cleavage and cell-cell fusion, and in particular escaped one subset of ADCC antibodies. Together with reduced immunogenicity in hybrid immunity background, these indicate a single spike amino acid glycosylation event provides selective advantage in humans through multiple mechanisms.

Role of N343 glycosylation on the SARS-CoV-2 S RBD structure and co-receptor binding across variants of concern.

Glycosylation of the SARS-CoV-2 spike (S) protein represents a key target for viral evolution because it affects both viral evasion and fitness. Successful variations in the glycan shield are difficult to achieve though, as protein glycosylation is also critical to folding and structural stability. Within this framework, the identification of glycosylation sites that are structurally dispensable can provide insight into the evolutionary mechanisms of the shield and inform immune surveillance. In this work, we show through over 45 μs of cumulative sampling from conventional and enhanced molecular dynamics (MD) simulations, how the structure of the immunodominant S receptor binding domain (RBD) is regulated by N-glycosylation at N343 and how this glycan's structural role changes from WHu-1, alpha (B.1.1.7), and beta (B.1.351), to the delta (B.1.617.2), and omicron (BA.1 and BA.2.86) variants. More specifically, we find that the amphipathic nature of the N-glycan is instrumental to preserve the structural integrity of the RBD hydrophobic core and that loss of glycosylation at N343 triggers a specific and consistent conformational change. We show how this change allosterically regulates the conformation of the receptor binding motif (RBM) in the WHu-1, alpha, and beta RBDs, but not in the delta and omicron variants, due to mutations that reinforce the RBD architecture. In support of these findings, we show that the binding of the RBD to monosialylated ganglioside co-receptors is highly dependent on N343 glycosylation in the WHu-1, but not in the delta RBD, and that affinity changes significantly across VoCs. Ultimately, the molecular and functional insight we provide in this work reinforces our understanding of the role of glycosylation in protein structure and function and it also allows us to identify the structural constraints within which the glycosylation site at N343 can become a hotspot for mutations in the SARS-CoV-2 S glycan shield.