SARS-CoV-2 variant of concern fitness and adaptation in primary human airway epithelia

The severe acute respiratory syndrome coronavirus 2 pandemic is characterized by the emergence of novel variants of concern (VOCs) that replace ancestral strains. Here, we dissect the complex selective pressures by evaluating variant fitness and adaptation in human respiratory tissues. We evaluate viral properties and host responses to reconstruct forces behind D614G through Omicron (BA.1) emergence. We observe differential replication in airway epithelia, differences in cellular tropism, and virus-induced cytotoxicity. D614G accumulates the most mutations after infection, supporting zoonosis and adaptation to the human airway. We perform head-to-head competitions and observe the highest fitness for Gamma and Delta. Under these conditions, RNA recombination favors variants encoding the B.1.617.1 lineage 3' end. Based on viral growth kinetics, Alpha, Gamma, and Delta exhibit increased fitness compared to D614G. In contrast, the global success of Omicron likely derives from increased transmission and antigenic variation. Our data provide molecular evidence to support epidemiological observations of VOC emergence.

U-CAN-seq: A Universal Competition Assay by Nanopore Sequencing

RNA viruses quickly evolve subtle genotypic changes that can have major impacts on viral fitness and host range, with potential consequences for human health. It is therefore important to understand the evolutionary fitness of novel viral variants relative to well-studied genotypes of epidemic viruses. Competition assays are an effective and rigorous system with which to assess the relative fitness of viral genotypes. However, it is challenging to quickly and cheaply distinguish and quantify fitness differences between very similar viral genotypes. Here, we describe a protocol for using reverse transcription PCR in combination with commercial nanopore sequencing services to perform competition assays on untagged RNA viruses. Our assay, called the Universal Competition Assay by Nanopore Sequencing (U-CAN-seq), is relatively cheap and highly sensitive. We used a well-studied N24A mutation in the chikungunya virus (CHIKV) nsp3 gene to confirm that we could detect a competitive disadvantage using U-CAN-seq. We also used this approach to show that mutations to the CHIKV 5' conserved sequence element that disrupt sequence but not structure did not affect the fitness of CHIKV. However, similar mutations to an adjacent CHIKV stem loop (SL3) did cause a fitness disadvantage compared to wild-type CHIKV, suggesting that structure-independent, primary sequence determinants in this loop play an important role in CHIKV biology. Our novel findings illustrate the utility of the U-CAN-seq competition assay.

Properties and Mechanisms of Deletions, Insertions, and Substitutions in the Evolutionary History of SARS-CoV-2

SARS-CoV-2 has accumulated many mutations since its emergence in late 2019. Nucleotide substitutions leading to amino acid replacements constitute the primary material for natural selection. Insertions, deletions, and substitutions appear to be critical for coronavirus's macro- and microevolution. Understanding the molecular mechanisms of mutations in the mutational hotspots (positions, loci with recurrent mutations, and nucleotide context) is important for disentangling roles of mutagenesis and selection. In the SARS-CoV-2 genome, deletions and insertions are frequently associated with repetitive sequences, whereas C>U substitutions are often surrounded by nucleotides resembling the APOBEC mutable motifs. We describe various approaches to mutation spectra analyses, including the context features of RNAs that are likely to be involved in the generation of recurrent mutations. We also discuss the interplay between mutations and natural selection as a complex evolutionary trend. The substantial variability and complexity of pipelines for the reconstruction of mutations and the huge number of genomic sequences are major problems for the analyses of mutations in the SARS-CoV-2 genome. As a solution, we advocate for the development of a centralized database of predicted mutations, which needs to be updated on a regular basis.

Neutralizing antibodies evolve to exploit vulnerable sites in the HCV envelope glycoprotein E2 and mediate spontaneous clearance of infection

Individuals who clear primary hepatitis C virus (HCV) infections clear subsequent reinfections more than 80% of the time, but the mechanisms are poorly defined. Here, we used HCV variants and plasma from individuals with repeated clearance to characterize longitudinal changes in envelope glycoprotein E2 sequences, function, and neutralizing antibody (NAb) resistance. Clearance of infection was associated with early selection of viruses with NAb resistance substitutions that also reduced E2 binding to CD81, the primary HCV receptor. Later, peri-clearance plasma samples regained neutralizing capacity against these variants. We identified a subset of broadly NAbs (bNAbs) for which these loss-of-fitness substitutions conferred resistance to unmutated bNAb ancestors but increased sensitivity to mature bNAbs. These data demonstrate a mechanism by which neutralizing antibodies contribute to repeated immune-mediated HCV clearance, identifying specific bNAbs that exploit fundamental vulnerabilities in E2. The induction of bNAbs with these specificities should be a goal of HCV vaccine development.

Cross-neutralization and viral fitness of SARS-CoV-2 Omicron sublineages

The rapid evolution of SARS-CoV-2 Omicron sublineages mandates a better understanding of viral replication and cross-neutralization among these sublineages. Here we used K18-hACE2 mice and primary human airway cultures to examine the viral fitness and antigenic relationship among Omicron sublineages. In both K18-hACE2 mice and human airway cultures, Omicron sublineages exhibited a replication order of BA.5 ≥ BA.2 ≥ BA.2.12.1 > BA.1; no difference in body weight loss was observed among different sublineage-infected mice. The BA.1-, BA.2-, BA.2.12.1-, and BA.5-infected mice developed distinguishable cross-neutralizations against Omicron sublineages, but exhibited little neutralization against the index virus (i.e. USA-WA1/2020) or the Delta variant. Surprisingly, the BA.5-infected mice developed higher neutralization activity against heterologous BA.2 and BA.2.12.1 than that against homologous BA.5; serum neutralizing titres did not always correlate with viral replication levels in infected animals. Our results revealed a distinct antigenic cartography of Omicron sublineages and support the bivalent vaccine approach.

Effect of the L74I Polymorphism on Fitness of Cabotegravir-Resistant Variants of Human Immunodeficiency Virus 1 Subtype A6

The presence of human immunodeficiency virus (HIV) 1 subtype A6, characterized by the L74I integrase (IN) polymorphism, is associated with confirmed virologic failure in clinical trials of long-acting cabotegravir and rilpivirine. We investigated the effect of L74I on replication capacity (RC) of recombinant viruses carrying this polymorphism in combination with various IN stand-transfer inhibitor resistance mutations. The presence of L74I conferred greater RC to recombinant viruses expressing HIV-1 A6 IN when present together with G118R, G140R, Q148H, and R263K; no significant difference in RC was observed for the Q148K or R mutants. These findings may explain, in part, the association of HIV-1 subtype A6 with virologic failure.

A mutation in the coronavirus nsp13-helicase impairs enzymatic activity and confers partial remdesivir resistance

Coronaviruses (CoVs) encode nonstructural proteins 1-16 (nsps 1-16) which form replicase complexes that mediate viral RNA synthesis. Remdesivir (RDV) is an adenosine nucleoside analog antiviral that inhibits CoV RNA synthesis. RDV resistance mutations have been reported only in the nonstructural protein 12 RNA-dependent RNA polymerase (nsp12-RdRp). We here show that a substitution mutation in the nsp13-helicase (nsp13-HEL A335V) of the betacoronavirus murine hepatitis virus (MHV) that was selected during passage with the RDV parent compound confers partial RDV resistance independently and additively when expressed with co-selected RDV resistance mutations in the nsp12-RdRp. The MHV A335V substitution did not enhance replication or competitive fitness compared to WT MHV and remained sensitive to the active form of the cytidine nucleoside analog antiviral molnupiravir (MOV). Biochemical analysis of the SARS-CoV-2 helicase encoding the homologous substitution (A336V) demonstrates that the mutant protein retained the ability to associate with the core replication proteins nsps 7, 8, and 12 but had impaired helicase unwinding and ATPase activity. Together, these data identify a novel determinant of nsp13-HEL enzymatic activity, define a new genetic pathway for RDV resistance, and demonstrate the importance of surveillance for and testing of helicase mutations that arise in SARS-CoV-2 genomes. IMPORTANCE Despite the development of effective vaccines against COVID-19, the continued circulation and emergence of new variants support the need for antivirals such as RDV. Understanding pathways of antiviral resistance is essential for surveillance of emerging variants, development of combination therapies, and for identifying potential new targets for viral inhibition. We here show a novel RDV resistance mutation in the CoV helicase also impairs helicase functions, supporting the importance of studying the individual and cooperative functions of the replicase nonstructural proteins 7-16 during CoV RNA synthesis. The homologous nsp13-HEL mutation (A336V) has been reported in the GISAID database of SARS-CoV-2 genomes, highlighting the importance of surveillance of and genetic testing for nucleoside analog resistance in the helicase.

Impact of mAb-induced A475V substitution on viral fitness and antibody neutralization of SARS-CoV-2 omicron variants in the presence of monoclonal antibodies and human convalescent sera

The emergence and rapid evolution of SARS-CoV-2 variants have posed a major challenge to the global efforts to control the COVID -19 pandemic. In this study, we investigated the potential of two SARS-CoV-2 variants, BA.2 and BA.5, to evade neutralization by a human monoclonal antibody targeting the virus's spike RBD (mAb 1D1). By subjecting the viruses to serial propagation in the presence of the antibody, we found that BA.2 exhibited poor growth, whereas BA.5 regained robust growth with significantly higher kinetics than the parental virus. Genetic analysis identified a single mutation, A475V, in the spike protein of BA.5 that substantially reduced the neutralizing activities of monoclonal antibodies and convalescent sera. In addition, the A475V mutation alone in BA.2 moderately reduced the neutralizing activity but completely abolished the neutralizing effect of mAb 1D1 when F486V or L452R were also present. Our results shed light on the possible evolutionary development of SARS-CoV-2 variants under selection pressure by monoclonal antibodies and have implications for the development of effective antibody therapies and vaccines against the virus.

Fitness-Dependent, Mild Mutagenic Activity of Sofosbuvir for Hepatitis C Virus

The concept of a mild mutagen was coined to describe a minor mutagenic activity exhibited by some nucleoside analogues that potentiated their efficacy as antiretroviral agents. In the present study, we report the mild mutagen activity of sofosbuvir (SOF) for hepatitis C virus (HCV). Serial passages of HCV in human hepatoma cells, in the presence of SOF at a concentration well below its cytotoxic concentration 50 (CC50) led to pre-extinction populations whose mutant spectra exhibited a significant increase of C→U transitions, relative to populations passaged in the absence of SOF. This was reflected in an increase in several diversity indices that were used to characterize viral quasispecies. The mild mutagenic activity of SOF was largely absent when it was tested with isogenic HCV populations that displayed high replicative fitness. Thus, SOF can act as a mild mutagen for HCV, depending on HCV fitness. Possible mechanisms by which the SOF mutagenic activity may contribute to its antiviral efficacy are discussed.

Enhanced fitness of SARS-CoV-2 B.1.617.2 Delta variant in ferrets

Competition assays were conducted in vitro and in vivo to examine how the Delta (B.1.617.2) variant displaced the prototype Washington/1/2020 (WA/1) strain. While WA/1 virus exhibited a moderately increased proportion compared to that in the inoculum following co-infection in human respiratory cells, Delta variant possessed a substantial in vivo fitness advantage as this virus becoming predominant in both inoculated and contact animals. This work identifies critical traits of the Delta variant that likely played a role in it becoming a dominant variant and highlights the necessities of employing multiple model systems to assess the fitness of newly emerged SARS-CoV-2 variants.

Human airway and nasal organoids reveal escalating replicative fitness of SARS-CoV-2 emerging variants

The high transmissibility of SARS-CoV-2 Omicron subvariants was generally ascribed to immune escape. It remained unclear whether the emerging variants have gradually acquired replicative fitness in human respiratory epithelial cells. We sought to evaluate the replicative fitness of BA.5 and earlier variants in physiologically active respiratory organoids. BA.5 exhibited a dramatically increased replicative capacity and infectivity than B.1.1.529 and an ancestral strain wildtype (WT) in human nasal and airway organoids. BA.5 spike pseudovirus showed a significantly higher entry efficiency than that carrying WT or B.1.1.529 spike. Notably, we observed prominent syncytium formation in BA.5-infected nasal and airway organoids, albeit elusive in WT- and B.1.1.529-infected organoids. BA.5 spike-triggered syncytium formation was verified by lentiviral overexpression of spike in nasal organoids. Moreover, BA.5 replicated modestly in alveolar organoids, with a significantly lower titer than B.1.1.529 and WT. Collectively, the higher entry efficiency and fusogenic activity of BA.5 spike potentiated viral spread through syncytium formation in the human airway epithelium, leading to enhanced replicative fitness and immune evasion, whereas the attenuated replicative capacity of BA.5 in the alveolar organoids may account for its benign clinical manifestation.

Viral Fitness of Baloxavir-Resistant Recombinant Influenza B/Victoria- and B/Yamagata-like Viruses Harboring the I38T PA Change, In Vitro, Ex Vivo and in Guinea Pigs

Seasonal influenza A and B viruses may cause severe infections requiring therapeutic interventions. Baloxavir, the latest antiviral drug approved against those infections, targets the endonuclease activity encoded by the polymerase acidic (PA) protein. While appearing effective at cessation of viral shedding, baloxavir demonstrated a low barrier of resistance. Herein, we aimed to assess the impact of PA-I38T substitution, a major marker of baloxavir-resistance, on the fitness of contemporary influenza B viruses. Recombinant wild-type (WT) influenza B/Phuket/2073/13 (B/Yamagata/16/88-like) and B/Washington/02/19 (B/Victoria/2/87-like) viruses and their respective PA-I38T mutants were used to evaluate replication kinetics in vitro, using A549 and Calu3 cells, and ex vivo, using nasal human airway epithelium (HAE) cells. Infectivity was also assessed in guinea pigs. In the B/Washington/02/19 background, there were no major differences between the recombinant WT virus and its I38T mutant when viral replication kinetics were evaluated in human lung cell lines and HAE as well as in nasal washes of experimentally infected guinea pigs. By contrast, the I38T mutation moderately impacted the B/Phuket/2073/13 viral fitness. In conclusion, contemporary influenza B viruses that may acquire baloxavir-resistance through the PA-I38T substitution could retain a significant level of fitness, highlighting the importance of monitoring the emergence of such variant.

Evaluation of SARS-CoV-2 ORF7a Deletions from COVID-19-Positive Individuals and Its Impact on Virus Spread in Cell Culture

The spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), causing the COVID-19 outbreak, posed a primary concern of public health worldwide. The most common changes in SARS-CoV-2 are single nucleotide substitutions, also reported insertions and deletions. This work investigates the presence of SARS-CoV-2 ORF7a deletions identified in COVID-19-positive individuals. Sequencing of SARS-CoV-2 complete genomes showed three different ORF7a size deletions (190-nt, 339-nt and 365-nt). Deletions were confirmed through Sanger sequencing. The ORF7a∆190 was detected in a group of five relatives with mild symptoms of COVID-19, and the ORF7a∆339 and ORF7a∆365 in a couple of co-workers. These deletions did not affect subgenomic RNAs (sgRNA) production downstream of ORF7a. Still, fragments associated with sgRNA of genes upstream of ORF7a showed a decrease in size when corresponding to samples with deletions. In silico analysis suggests that the deletions impair protein proper function; however, isolated viruses with partial deletion of ORF7a can replicate in culture cells similarly to wild-type viruses at 24 hpi, but with less infectious particles after 48 hpi. These findings on deleted ORF7a accessory protein gene, contribute to understanding SARS-CoV-2 phenotypes such as replication, immune evasion and evolutionary fitness as well insights into the role of SARS-CoV-2_ORF7a in the mechanism of virus-host interactions.

Efficacy decrease of antiviral agents when administered to ongoing hepatitis C virus infections in cell culture

We report a quantification of the decrease of effectiveness of antiviral agents directed to hepatitis C virus, when the agents are added during an ongoing infection in cell culture vs. when they are added at the beginning of the infection. Major determinants of the decrease of inhibitory activity are the time post-infection of inhibitor administration and viral replicative fitness. The efficacy decrease has been documented with antiviral assays involving the combination of the direct-acting antiviral agents, daclatasvir and sofosbuvir, and with the combination of the lethal mutagens, favipiravir and ribavirin. The results suggest that strict antiviral effectiveness assays in preclinical trials may involve the use of high fitness viral populations and the delayed administration of the agents, relative to infection onset.

Inspecting the Ribozyme Region of Hepatitis Delta Virus Genotype 1: Conservation and Variability

The hepatitis delta virus (HDV) genome has an autocatalytic region called the ribozyme, which is essential for viral replication. The aim of this study was to use next-generation sequencing (NGS) to analyze the ribozyme quasispecies (QS) in order to study its evolution and identify highly conserved regions potentially suitable for a gene-silencing strategy. HDV RNA was extracted from 2 longitudinal samples of chronic HDV patients and the ribozyme (nucleotide, nt 688-771) was analyzed using NGS. QS conservation, variability and genetic distance were analyzed. Mutations were identified by aligning sequences with their specific genotype consensus. The main relevant mutations were tested in vitro. The ribozyme was conserved overall, with a hyper-conserved region between nt 715-745. No difference in QS was observed over time. The most variable region was between nt 739-769. Thirteen mutations were observed, with three showing a higher frequency: T23C, T69C and C64 deletion. This last strongly reduced HDV replication by more than 1 log in vitro. HDV Ribozyme QS was generally highly conserved and was maintained during follow-up. The most conserved portion may be a valuable target for a gene-silencing strategy. The presence of the C64 deletion may strongly impair viral replication, as it is a potential mechanism of viral persistence.

Spontaneous Mutation in the Movement Protein of Citrus Leprosis Virus C2, in a Heterologous Virus Infection Context, Increases Cell-to-Cell Transport and Generates Fitness Advantage

Previous results using a movement defective alfalfa mosaic virus (AMV) vector revealed that citrus leprosis virus C (CiLV-C) movement protein (MP) generates a more efficient local movement, but not more systemic transport, than citrus leprosis virus C2 (CiLV-C2) MP, MPs belonging to two important viruses for the citrus industry. Here, competition experiment assays in transgenic tobacco plants (P12) between transcripts of AMV constructs expressing the cilevirus MPs, followed by several biological passages, showed the prevalence of the AMV construct carrying the CiLV-C2 MP. The analysis of AMV RNA 3 progeny recovered from P12 plant at the second viral passage revealed the presence of a mix of progeny encompassing the CiLV-C2 MP wild type (MP^WT) and two variants carrying serines instead phenylalanines at positions 72 (MP^S72F) or 259 (MP^S259F), respectively. We evaluated the effects of each modified residue in virus replication, and cell-to-cell and long-distance movements. Results indicated that phenylalanine at position 259 favors viral cell-to-cell transport with an improvement in viral fitness, but has no effect on viral replication, whereas mutation at position 72 (MP^S72F) has a penalty in the viral fitness. Our findings indicate that the prevalence of a viral population may be correlated with its greater efficiency in cell-to-cell and systemic movements.

No Evidence of the Vertical Transmission of Non-Virulent Infectious Salmon Anaemia Virus (ISAV-HPR0) in Farmed Atlantic Salmon

The nonvirulent infectious salmon anaemia virus (ISAV-HPR0) is the putative progenitor for virulent-ISAV, and a potential risk factor for the development of infectious salmon anaemia (ISA). Understanding the transmission dynamics of ISAV-HPR0 is fundamental to proper management and mitigation strategies. Here, we demonstrate that ISAV-HPR0 causes prevalent and transient infections in all three production stages of Atlantic salmon in the Faroe Islands. Phylogenetic analysis of the haemagglutinin-esterase gene from 247 salmon showed a clear geographical structuring into two significantly distinct HPR0-subgroups, which were designated G2 and G4. Whereas G2 and G4 co-circulated in marine farms, Faroese broodfish were predominantly infected by G2, and smolt were predominantly infected by G4. This infection pattern was confirmed by our G2- and G4-specific RT-qPCR assays. Moreover, the HPR0 variants detected in Icelandic and Norwegian broodfish were never detected in the Faroe Islands, despite the extensive import of ova from both countries. Accordingly, the vertical transmission of HPR0 from broodfish to progeny is uncommon. Phylogenetic and statistical analysis suggest that HPR0 persists in the smolt farms as “house-strains”, and that new HPR0 variants are occasionally introduced from the marine environment, probably by HPR0-contaminated sea-spray. Thus, high biosecurity—including water and air intake—is required to avoid the introduction of pathogens to the smolt farms.

Revisiting dengue virus-mosquito interactions: molecular insights into viral fitness

Dengue virus (DENV), like other viruses, closely interacts with the host cell machinery to complete its life cycle. Over the course of infection, DENV interacts with several host factors with pro-viral activities to support its infection. Meanwhile, it has to evade or counteract host factors with anti-viral activities which inhibit its infection. These molecular virus-host interactions play a crucial role in determining the success of DENV infection. Deciphering such interactions is thus paramount to understanding viral fitness in its natural hosts. While DENV-mammalian host interactions have been extensively studied, not much has been done to characterize DENV-mosquito host interactions despite its importance in controlling DENV transmission. Here, to provide a snapshot of our current understanding of DENV-mosquito interactions, we review the literature that identified host factors and cellular processes related to DENV infection in its mosquito vectors, Aedes aegypti and Aedes albopictus, with a particular focus on DENV-mosquito omics studies. This knowledge provides fundamental insights into the DENV life cycle, and could contribute to the development of novel antiviral strategies.

Viral evolution sustains a dengue outbreak of enhanced severity

Compared to the previous 2013-2014 outbreak, dengue 2016-2017 outbreak in New Caledonia was characterized by an increased number of severe forms associated with hepatic presentations. In this study, we assessed the virological factors associated with this enhanced severity. Whole-genome sequences were retrieved from dengue virus (DENV)-1 strains collected in 2013-2014 and from severe and non-severe patients in 2016-2017. Fitness, hepatic tropism and cytopathogenicity of DENV 2016-2017 strains were compared to those of 2013-2014 strains using replication kinetics in the human hepatic cell line HuH7. Whole-genome sequencing identified four amino acid substitutions specific to 2016-2017 strains and absent from 2013-2014 strains. Three of these mutations occurred in predicted T cell epitopes, among which one was also a B cell epitope. Strains retrieved from severe forms did not exhibit specific genetic features. DENV strains from 2016-2017 exhibited a trend towards reduced replicative fitness and cytopathogenicity in vitro compared to strains from 2013-2014. Overall, the 2016-2017 dengue outbreak in New Caledonia was associated with a viral genetic evolution which had limited impact on DENV hepatic tropism and cytopathogenicity. These mutations, however, may have modified DENV strains antigenicity, altering the anti-DENV immune response in some patients, in turn favoring the development of severe forms.Trial registration: ClinicalTrials.gov identifier: NCT04615364.

Epistatic interactions promote persistence of NS3-Q80K in HCV infection by compensating for protein folding instability

The Q80K polymorphism in the NS3-4A protease of the hepatitis C virus is associated with treatment failure of direct-acting antiviral agents. This polymorphism is highly prevalent in genotype 1a infections and stably transmitted between hosts. Here, we investigated the underlying molecular mechanisms of evolutionarily conserved coevolving amino acids in NS3-Q80K and revealed potential implications of epistatic interactions in immune escape and variants persistence. Using purified protein, we characterized the impact of epistatic amino acid substitutions on the physicochemical properties and peptide cleavage kinetics of the NS3-Q80K protease. We found that Q80K destabilized the protease protein fold (p < 0.0001). Although NS3-Q80K showed reduced peptide substrate turnover (p < 0.0002), replicative fitness in an H77S.3 cell culture model of infection was not significantly inferior to the WT virus. Epistatic substitutions at residues 91 and 174 in NS3-Q80K stabilized the protein fold (p < 0.0001) and leveraged the WT protease stability. However, changes in protease stability inversely correlated with enzymatic activity. In infectious cell culture, these secondary substitutions were not associated with a gain of replicative fitness in NS3-Q80K variants. Using molecular dynamics, we observed that the total number of residue contacts in NS3-Q80K mutants correlated with protein folding stability. Changes in the number of contacts reflected the compensatory effect on protein folding instability by epistatic substitutions. In summary, epistatic substitutions in NS3-Q80K contribute to viral fitness by mechanisms not directly related to RNA replication. By compensating for protein-folding instability, epistatic interactions likely protect NS3-Q80K variants from immune cell recognition.

Prediction and Evolution of the Molecular Fitness of SARS-CoV-2 Variants: Introducing SpikePro

The understanding of the molecular mechanisms driving the fitness of the SARS-CoV-2 virus and its mutational evolution is still a critical issue. We built a simplified computational model, called SpikePro, to predict the SARS-CoV-2 fitness from the amino acid sequence and structure of the spike protein. It contains three contributions: the inter-human transmissibility of the virus predicted from the stability of the spike protein, the infectivity computed in terms of the affinity of the spike protein for the ACE2 receptor, and the ability of the virus to escape from the human immune response based on the binding affinity of the spike protein for a set of neutralizing antibodies. Our model reproduces well the available experimental, epidemiological and clinical data on the impact of variants on the biophysical characteristics of the virus. For example, it is able to identify circulating viral strains that, by increasing their fitness, recently became dominant at the population level. SpikePro is a useful, freely available instrument which predicts rapidly and with good accuracy the dangerousness of new viral strains. It can be integrated and play a fundamental role in the genomic surveillance programs of the SARS-CoV-2 virus that, despite all the efforts, remain time-consuming and expensive.

Prediction and Evolution of the Molecular Fitness of SARS-CoV-2 Variants: Introducing SpikePro

The understanding of the molecular mechanisms driving the fitness of the SARS-CoV-2 virus and its mutational evolution is still a critical issue. We built a simplified computational model, called SpikePro, to predict the SARS-CoV-2 fitness from the amino acid sequence and structure of the spike protein. It contains three contributions: the inter-human transmissibility of the virus predicted from the stability of the spike protein, the infectivity computed in terms of the affinity of the spike protein for the ACE2 receptor, and the ability of the virus to escape from the human immune response based on the binding affinity of the spike protein for a set of neutralizing antibodies. Our model reproduces well the available experimental, epidemiological and clinical data on the impact of variants on the biophysical characteristics of the virus. For example, it is able to identify circulating viral strains that, by increasing their fitness, recently became dominant at the population level. SpikePro is a useful, freely available instrument which predicts rapidly and with good accuracy the dangerousness of new viral strains. It can be integrated and play a fundamental role in the genomic surveillance programs of the SARS-CoV-2 virus that, despite all the efforts, remain time-consuming and expensive.

The Signatures of Natural Selection and Molecular Evolution in Fusarium graminearum Virus 1

Fusarium graminearum virus 1 (FgV1) is a positive-sense ssRNA virus that confers hypovirulence in its fungal host, Fusarium graminearum. Like most mycoviruses, FgV1 exists in fungal cells, lacks an extracellular life cycle, and is therefore transmitted during sporulation or hyphal anastomosis. To understand FgV1 evolution and/or adaptation, we conducted mutation accumulation (MA) experiments by serial passage of FgV1 alone or with FgV2, 3, or 4 in F. graminearum. We expected that the effects of positive selection would be highly limited because of repeated bottleneck events. To determine whether selection on the virus was positive, negative, or neutral, we assessed both the phenotypic traits of the host fungus and the RNA sequences of FgV1. We inferred that there was positive selection on beneficial mutations in FgV1 based on the ratio of non-synonymous to synonymous substitutions (d_N/d_S), on the ratio of radical to conservation amino acid replacements (p_NR/p_NC), and by changes in the predicted protein structures. In support of this inference, we found evidence of positive selection only in the open reading frame 4 (ORF4) protein of DK21/FgV1 (MA line 1); mutations at amino acids 163A and 289H in the ORF4 of MA line 1 affected the entire structure of the protein predicted to be under positive selection. We also found, however, that deleterious mutations were a major driving force in viral evolution during serial passages. Linear relationships between changes in viral fitness and the number of mutations in each MA line demonstrated that some deleterious mutations resulted in fitness decline. Several mutations in MA line 1 were not shared with any of the other four MA lines (PH-1/FgV1, PH-1/FgV1 + 2, PH-1/FgV1 + 3, and PH-1/FgV1 + 4). This suggests that evolutionary pathways of the virus could differ with respect to hosts and also with respect to co-infecting viruses. The data also suggested that the differences among MA lines might also be explained by mutational robustness and other unidentified factors. Additional research is needed to clarify the effects of virus co-infection on the adaptation or evolution of FgV1 to its environments.

HIV-1 Subtype C Tier 3 Viruses Have Increased Infectivity Compared to Tier 2 Viruses

A primary concern of an antibody-based HIV-1 therapy is the virus' ability to rapidly escape antibody responses. Therefore, we investigated the relationships between antibody neutralization sensitivity, viral phenotype, and infectivity in 13 subtype C viruses using a HeLa transfectant-based assay. We observed that the seven tier 3 viruses exhibited higher infectivity than the tier 2 viruses, suggesting that higher neutralization resistance did not have a substantial entry cost. There was no relationship between neutralization resistance and susceptibility to entry inhibitors Maraviroc, PSC RANTES, or the fusion inhibitor T20, indicating that neutralization resistance may not alter these inhibitor target sites. By analyzing glycosylation patterns in 82 subtype C viruses, we found that the presence of an N-linked glycan motif at position N413 and its absence at N332 were the most important predictors of neutralization resistance. In a set of 200 subtype C viruses, tier 3 strains were more resistant than tier 2 or 1B viruses to several broadly neutralizing monoclonal antibodies targeting three different epitopes. This suggests that it is unlikely that resistance to antibodies targeting a single epitope drives overall resistance. In the context of an antibody-based intervention, highly resistant viruses with increased infectivity, circulating in the population, could hinder HIV-1 control since entry of tier 3 viruses is not always selected against. Therefore, for any long-term antibody-based intervention to be globally relevant, it must elicit responses that limit the occurrence of resistance.

Amino Acid Substitutions Associated with Treatment Failure for Hepatitis C Virus Infection

Despite the high virological response rates achieved with current directly acting antiviral agents (DAAs) against hepatitis C virus (HCV), around 2% to 5% of treated patients do not achieve a sustained viral response. The identification of amino acid substitutions associated with treatment failure requires analytical designs, such as subtype-specific ultradeep sequencing (UDS) methods, for HCV characterization and patient management. Using this procedure, we have identified six highly represented amino acid substitutions (HRSs) in NS5A and NS5B of HCV, which are not bona fide resistance-associated substitutions (RAS), from 220 patients who failed therapy. They were present frequently in basal and posttreatment virus of patients who failed different DAA-based therapies. Contrary to several RAS, HRSs belong to the acceptable subset of substitutions according to the PAM250 replacement matrix. Their mutant frequency, measured by the number of deep sequencing reads within the HCV quasispecies that encode the relevant substitutions, ranged between 90% and 100% in most cases. They also have limited predicted disruptive effects on the three-dimensional structures of the proteins harboring them. Possible mechanisms of HRS origin and dominance, as well as their potential predictive value for treatment response, are discussed.

West Nile Virus fidelity modulates the capacity for host cycling and adaptation

The fidelity of flaviviruses is thought to be tightly regulated for optimal fitness within and between hosts. West Nile virus (WNV) high-fidelity (HiFi) mutations V793I and G806R within the RNA-dependent RNA polymerase, and low-fidelity (LoFi) mutation T248I within the methyltransferase, were previously shown to attenuate infectivity and replicative fitness in Culex mosquitoes and Culex tarsalis (CXT) cells but not in mammalian cells. We hypothesized that fidelity alterations would modify adaptation and maintenance in a host-specific manner. To test this hypothesis, wild-type (WT), HiFi (V793I/G806R) and LoFi (T248I) variants were sequentially passaged eight times in avian (PDE) or mosquito cells, or alternately between the two. Initial characterization confirmed that fidelity mutants are attenuated in mosquito, but not avian, cells. Deep sequencing revealed mutations unique to both cell lines and fidelity mutants, including ENV G1378A, a mutation associated with avian cell adaptation. To characterize maintenance and adaptation, viral outputs were monitored throughout passaging and viral fitness was assessed. The results indicate that fidelity mutants can at times recover fitness during mosquito cell passage, but remain attenuated relative to WT. Despite similar initial fitness, LoFi mutants were impaired during sequential passage in avian cells. Conversely, HiFi mutants passaged in avian cells showed increased adaptation, suggesting that increased fidelity may be advantageous in avian hosts. Although some adaptation occurred with individual mutants, the output titres of fidelity mutants were on average lower and were often lost during host switching. These data confirm that arbovirus fidelity is likely fine-tuned to maximize survival in disparate hosts.

In vitro evolution of herpes simplex virus 1 (HSV-1) reveals selection for syncytia and other minor variants in cell culture

The large dsDNA virus herpes simplex virus 1 (HSV-1) is considered to be genetically stable, yet it can rapidly evolve in response to strong selective pressures such as antiviral treatment. Deep sequencing has revealed that clinical and laboratory isolates of this virus exist as populations that contain a mixture of minor alleles or variants, similar to many RNA viruses. The classic virology approach of plaque purifying virus creates a genetically homogenous population, but it is not clear how closely this represents the mixed virus populations found in nature. We sought to study the evolution of mixed versus highly purified HSV-1 populations in controlled cell culture conditions, to examine the impact of this genetic diversity on evolution. We found that a mixed population of HSV-1 acquired more genetic diversity and underwent a more dramatic phenotypic shift than a plaque-purified population, producing a viral population that was almost entirely syncytial after just ten passages. At the genomic level, adaptation and genetic diversification occurred at the level of minor alleles or variants in the viral population. Certain genetic variants in the mixed viral population appeared to be positively selected in cell culture, and this shift was also observed in clinical samples during their first passages in vitro. In contrast, the plaque-purified viral population did not appear to change substantially in phenotype or overall quantity of minor allele diversity. These data indicate that HSV-1 is capable of evolving rapidly in a given environment, and that this evolution is facilitated by diversity in the viral population.

The PB2 and M genes are critical for the superiority of genotype S H9N2 virus to genotype H in optimizing viral fitness of H5Nx and H7N9 avian influenza viruses in mice

Genotype S H9N2 avian influenza virus, which has been predominant in China since 2010, contributed its entire internal gene cassette to the genesis of novel reassortant influenza viruses, including H5Nx, H7N9 and H10N8 viruses that pose great threat to poultry and humans. A key feature of the genotype S H9N2 virus is the substitution of G1-like M and PB2 genes for the earlier F/98-like M and PB2 of genotype H virus. However, how this gene substitution has influenced viral adaptability of emerging influenza viruses in mammals remains unclear. We report here that reassortant H5Nx and H7N9 viruses with the genotype S internal gene cassette displayed enhanced replication and virulence over those with genotype H internal gene cassette in cell cultures as well as in the mouse models. We showed that the G1-like PB2 gene was associated with increased polymerase activity and improved nuclear accumulation compared with the F/98-like counterpart, while the G1-like M gene facilitated effective translocation of RNP to cytoplasm. Our findings suggest that the genotype S H9N2 internal gene cassette, which possesses G1-like M and PB2 genes, is superior to that of genotype H, in optimizing viral fitness, and thus have implications for assessing the potential risk of these gene introductions to generate emerging influenza viruses.

Emergence and Persistence of Letermovir-Resistant Cytomegalovirus in a Patient With Primary Immunodeficiency

Background

Letermovir is a novel cytomegalovirus antiviral that is approved for prophylaxis in hematopoietic stem cell transplantation recipients

Methods

After obtaining informed consent, letermovir prophylaxis was started in a patient with a presumed late-onset primary, combined T- and B-cell immunodeficiency. Plasma CMV DNAemia was monitored with real-time polymerase chain reaction, and letermovir resistance analyses were performed using Sanger sequencing and Illumina MiSeq next-generation sequencing.

Results

A letermovir-resistant cytomegalovirus variant (C325Y mutation in UL56) emerged 17 weeks after start of prophylaxis. The letermovir-resistant variant was able to reactivate without drug selective pressure as this variant was again detected in plasma 20.6 weeks after stopping of letermovir.

Conclusions

This case indicates that the C325Y mutation in UL56 does not significantly alter fitness of cytomegalovirus in vivo.

Dengue genetic divergence generates within-serotype antigenic variation, but serotypes dominate evolutionary dynamics

Dengue virus (DENV) exists as four genetically distinct serotypes, each of which is historically assumed to be antigenically uniform. Recent analyses suggest that antigenic heterogeneity may exist within each serotype, but its source, extent and impact remain unclear. Here, we construct a sequence-based model to directly map antigenic change to underlying genetic divergence. We identify 49 specific substitutions and four colinear substitution clusters that robustly predict dengue antigenic relationships. We report moderate antigenic diversity within each serotype, resulting in genotype-specific patterns of heterotypic cross-neutralization. We also quantify the impact of antigenic variation on real-world DENV population dynamics, and find that serotype-level antigenic fitness is a dominant driver of dengue clade turnover. These results provide a more nuanced understanding of the relationship between dengue genetic and antigenic evolution, and quantify the effect of antigenic fitness on dengue evolutionary dynamics.

In Vitro and In Vivo Characterization of Novel Neuraminidase Substitutions in Influenza A(H1N1)pdm09 Virus Identified Using Laninamivir-Mediated In Vitro Selection

Neuraminidase (NA) inhibitors (NAIs) are widely used antiviral drugs for the treatment of humans with influenza virus infections. There have been widespread reports of NAI resistance among seasonal A(H1N1) viruses, and most have been identified in oseltamivir-exposed patients or those treated with other NAIs. Thus, monitoring and identifying NA markers conferring resistance to NAIs-particularly newly introduced treatments-are critical to the management of viral infections. Therefore, we screened and identified substitutions conferring resistance to laninamivir by enriching random mutations in the NA gene of the 2009 pandemic influenza [A(H1N1)pdm09] virus followed by deep sequencing of the laninamivir-selected variants. After the generation of single mutants possessing each identified mutation, two A(H1N1)pdm09 recombinants possessing novel NA gene substitutions (i.e., D199E and P458T) were shown to exhibit resistance to more than one NAI. Of note, mutants possessing P458T-which is located outside of the catalytic or framework residue of the NA active site-exhibited highly reduced inhibition by all four approved NAIs. Using MDCK cells, we observed that the in vitro viral replication of the two recombinants was lower than that of the wild type (WT). Additionally, in infected mice, decreased mortality and/or mean lung viral titers were observed in mutants compared with the WT. Reverse mutations to the WT were observed in lung homogenate samples from D199E-infected mice after 3 serial passages. Overall, the novel NA substitutions identified could possibly emerge in influenza A(H1N1)pdm09 viruses during laninamivir therapy and the viruses could have altered NAI susceptibility, but the compromised in vitro/in vivo viral fitness may limit viral spreading.IMPORTANCE With the widespread emergence of NAI-resistant influenza virus strains, continuous monitoring of mutations that confer antiviral resistance is needed. Laninamivir is the most recently approved NAI in several countries; few data exist related to the in vitro selection of viral mutations conferring resistance to laninamivir. Thus, we screened and identified substitutions conferring resistance to laninamivir by random mutagenesis of the NA gene of the 2009 pandemic influenza [A(H1N1)pdm09] virus strain followed by deep sequencing of the laninamivir-selected variants. We found several novel substitutions in NA (D199E and P458T) in an A(H1N1)pdm09 background which conferred resistance to NAIs and which had an impact on viral fitness. Our study highlights the importance of continued surveillance for potential antiviral-resistant variants and the development of alternative therapeutics.

Contemporary Zika Virus Isolates Induce More dsRNA and Produce More Negative-Strand Intermediate in Human Astrocytoma Cells

The recent emergence and rapid geographic expansion of Zika virus (ZIKV) poses a significant challenge for public health. Although historically causing only mild febrile illness, recent ZIKV outbreaks have been associated with more severe neurological complications, such as Guillain-Barré syndrome and fetal microcephaly. Here we demonstrate that two contemporary (2015) ZIKV isolates from Puerto Rico and Brazil may have increased replicative fitness in human astrocytoma cells. Over a single infectious cycle, the Brazilian isolate replicates to higher titers and induces more severe cytopathic effects in human astrocytoma cells than the historical African reference strain or an early Asian lineage isolate. In addition, both contemporary isolates induce significantly more double-stranded RNA in infected astrocytoma cells, despite similar numbers of infected cells across isolates. Moreover, when we quantified positive- and negative-strand viral RNA, we found that the Asian lineage isolates displayed substantially more negative-strand replicative intermediates than the African lineage isolate in human astrocytoma cells. However, over multiple rounds of infection, the contemporary ZIKV isolates appear to be impaired in cell spread, infecting a lower proportion of cells at a low MOI despite replicating to similar or higher titers. Taken together, our data suggests that contemporary ZIKV isolates may have evolved mechanisms that allow them to replicate with increased efficiency in certain cell types, thereby highlighting the importance of cell-intrinsic factors in studies of viral replicative fitness.

Characterization of Influenza B Virus Variants with Reduced Neuraminidase Inhibitor Susceptibility

Treatment options for influenza B virus infections are limited to neuraminidase inhibitors (NAIs), which block the neuraminidase (NA) glycoprotein on the virion surface. The development of NAI resistance would therefore result in a loss of antiviral treatment options for influenza B virus infections. This study characterized two contemporary influenza B viruses with known resistance-conferring NA amino acid substitutions, D197N and H273Y, detected during routine surveillance. The D197N and H273Y variants were characterized in vitro by assessing NA enzyme activity and affinity, as well as replication in cell culture compared to those of NAI-sensitive wild-type viruses. In vivo studies were also performed in ferrets to assess the replication and transmissibility of each variant. Mathematical models were used to analyze within-host and between-host fitness of variants relative to wild-type viruses. The data revealed that the H273Y variant had NA enzyme function similar to that of its wild type but had slightly reduced replication and transmission efficiency in vivo The D197N variant had impaired NA enzyme function, but there was no evidence of reduction in replication or transmission efficiency in ferrets. Our data suggest that the influenza B virus variant with the H273Y NA substitution had a more notable reduction in fitness compared to wild-type viruses than the influenza B variant with the D197N NA substitution. Although a D197N variant is yet to become widespread, it is the most commonly detected NAI-resistant influenza B virus in surveillance studies. Our results highlight the need to carefully monitor circulating viruses for the spread of influenza B viruses with the D197N NA substitution.

Evolution and persistence of resistance-associated substitutions of hepatitis C virus after direct-acting antiviral treatment failures

Daclatasvir plus asunaprevir (DCV+ASV) treatment is an all-oral direct-acting antiviral (DAA) therapy for the genotype 1b HCV-infected patients. In this study, we investigated how resistance-associated substitutions (RASs) evolved after treatment failures and assessed the effect of those substitutions on viral fitness. Sequencing of NS5A and NS3 revealed typical RASs after treatment failures. Interestingly, the RASs of NS3 reverted to the wild-type amino acid within 1 year after treatment failures. However, the RASs of NS5A were stable and did not change. The effect of NS5A and NS3 RASs on viral RNA replication was assessed after mutagenic substitution in the genotype 1b HCV RNA. Among single substitutions, the effect of D168V was more substantial than the others and the effect of the triple mutant combination (D168V+L31V+Y93H) was the most severe. The RAS at NS5A Y93 affected both viral RNA replication and virus production. Finally, the effect of trans-complementation of NS5A was demonstrated in our co-transfection experiments and these results suggest that such a trans-complementation effect of NS5A may help maintain the NS5A RASs for a long time even after cessation of the DAA treatment. In conclusion, the results from this investigation would help understand the emergence and persistence of RASs.

CD8+ T-Cell Response-Associated Evolution of Hepatitis B Virus Core Protein and Disease Progress. J Virol. 2018;92:e02120-17. [PMID: 29950410 DOI: 10.1128/jvi.02120-17]

Under the immune pressure of cytotoxic T cells (CTLs), hepatitis B virus (HBV) evolves to accumulate mutations more likely within epitopes to evade immune detection. However, little is known about the specific patterns of the immune pressure-associated HBV mutation of T-cell epitopes and their link to disease progression. Here, we observed a correlation of the accumulated variants on HBV core protein (HBc) with the disease severity of HBV infection. Further analysis indicated that these substitutions were mostly located within CD8⁺ T-cell epitopes of HBc protein, which were systematically screened and identified in an unbiased manner in our study. From individual peptide level to the human leukocyte antigen I (HLA-I)-restricted population level, we elucidated that the mutations in these well-defined HLA-I-restricted T-cell epitopes significantly decreased antiviral activity-specific CTLs and were positively associated with clinical parameters and disease progression in HBV-infected patients. The molecular pattern for viral epitope variations based on the sequencing of 105 HBV virus genomes indicated that the C-terminal portion (Pc), especially the Pc-1 and Pc-2 positions, have the highest mutation rates. Further structural analysis of HLA-A*02 complexed to diverse CD8⁺ T-cell epitopes revealed that the highly variable C-terminal bulged peak of M-shaped HBc-derived epitopes are solvent exposed, and most of the CDR3βs of the T-cell receptor hover over them. These data shed light on the molecular and immunological mechanisms of T-cell immunity-associated viral evolution in hepatitis B progression, which is beneficial for designing immunotherapies and vaccines.IMPORTANCE The specific patterns of sequence polymorphisms of T-cell epitopes and the immune mechanisms of the HBV epitope mutation-linked disease progression are largely unclear. In this study, we systematically evaluated the contribution of CD8⁺ T cells to the disease progress-associated evolution of HBV. By evaluation of patient T-cell responses based on the peptide repertoire, we comprehensively characterized the association of clinical parameters in chronic hepatitis B with the antiviral T-cell response-associated mutations of the viruses from the single-epitope level to the overall HLA-I-restricted peptide levels. Furthermore, we investigated the molecular basis of the HLA-A2-restricted peptide immune escape and found that the solvent-exposed C-terminal portion of the epitopes is highly variable under CDR3β recognition. Our work may provide a comprehensive evaluation of viral mutations impacted by the host CTL response in HBV disease progression in the context of the full repertoire of HBc-derived epitopes.

The HIV-1 Reverse Transcriptase A62V Mutation Influences Replication Fidelity and Viral Fitness in the Context of Multi-Drug-Resistant Mutations

Emergence of human immunodeficiency virus type 1 (HIV-1) drug resistance arises from mutation fixation in the viral genome during antiretroviral therapy. Primary mutations directly confer antiviral drug resistance, while secondary mutations arise that do not confer drug resistance. The A62V amino acid substitution in HIV-1 reverse transcriptase (RT) was observed to be associated with multi-drug resistance, but is not known to be a resistance-conferring mutation. In particular, A62V was observed in various multi-dideoxynucleoside resistant (MDR) mutation complexes, including the Q151M complex (i.e., A62V, V75I, F77L, F116Y, and Q151M), and the T69SSS insertion complex, which has a serine–serine insertion between amino acid positions 69 and 70 (i.e., M41L, A62V, T69SSS, K70R, and T215Y). However, what selective advantage is conferred to the virus remains unresolved. In this study, we hypothesized that A62V could influence replication fidelity and viral fitness with viruses harboring the Q151M and T69SSS MDR mutation complexes. A single-cycle replication assay and a dual-competition fitness assay were used to assess viral mutant frequency and viral fitness, respectively. A62V was found to increase the observed lower mutant frequency identified with each of the viruses harboring the MDR mutation complexes in the single-cycle assay. Furthermore, A62V was observed to improve viral fitness of replication-competent MDR viruses. Taken together, these observations indicate an adaptive role of A62V in virus replication fidelity and viral fitness, which would likely enhance virus persistence during drug-selective pressure.

Viral Entry Properties Required for Fitness in Humans Are Lost through Rapid Genomic Change during Viral Isolation

Human parainfluenza viruses cause a large burden of human respiratory illness. While much research relies upon viruses grown in cultured immortalized cells, human parainfluenza virus 3 (HPIV-3) evolves in culture. Cultured viruses differ in their properties compared to clinical strains. We present a genome-wide survey of HPIV-3 adaptations to culture using metagenomic next-generation sequencing of matched pairs of clinical samples and primary culture isolates (zero passage virus). Nonsynonymous changes arose during primary viral isolation, almost entirely in the genes encoding the two surface glycoproteins-the receptor binding protein hemagglutinin-neuraminidase (HN) or the fusion protein (F). We recovered genomes from 95 HPIV-3 primary culture isolates and 23 HPIV-3 strains directly from clinical samples. HN mutations arising during primary viral isolation resulted in substitutions at HN's dimerization/F-interaction site, a site critical for activation of viral fusion. Alterations in HN dimer interface residues known to favor infection in culture occurred within 4 days (H552 and N556). A novel cluster of residues at a different face of the HN dimer interface emerged (P241 and R242) and imply a role in HPIV-3-mediated fusion. Functional characterization of these culture-associated HN mutations in a clinical isolate background revealed acquisition of the fusogenic phenotype associated with cultured HPIV-3; the HN-F complex showed enhanced fusion and decreased receptor-cleaving activity. These results utilize a method for identifying genome-wide changes associated with brief adaptation to culture to highlight the notion that even brief exposure to immortalized cells may affect key viral properties and underscore the balance of features of the HN-F complex required for fitness by circulating viruses.IMPORTANCE Human parainfluenza virus 3 is an important cause of morbidity and mortality among infants, the immunocompromised, and the elderly. Using deep genomic sequencing of HPIV-3-positive clinical material and its subsequent viral isolate, we discover a number of known and novel coding mutations in the main HPIV-3 attachment protein HN during brief exposure to immortalized cells. These mutations significantly alter function of the fusion complex, increasing fusion promotion by HN as well as generally decreasing neuraminidase activity and increasing HN-receptor engagement. These results show that viruses may evolve rapidly in culture even during primary isolation of the virus and before the first passage and reveal features of fitness for humans that are obscured by rapid adaptation to laboratory conditions.

Differential Shape of Geminivirus Mutant Spectra Across Cultivated and Wild Hosts With Invariant Viral Consensus Sequences

Geminiviruses (family Geminiviridae) possess single-stranded circular DNA genomes that are replicated by cellular polymerases in plant host cell nuclei. In their hosts, geminivirus populations behave as ensembles of mutant and recombinant genomes, known as viral quasispecies. This favors the emergence of new geminiviruses with altered host range, facilitating new or more severe diseases or overcoming resistance traits. In warm and temperate areas several whitefly-transmitted geminiviruses of the genus Begomovirus cause the tomato yellow leaf curl disease (TYLCD) with significant economic consequences. TYLCD is frequently controlled in commercial tomatoes by using the dominant Ty-1 resistance gene. Over a 45 day period we have studied the diversification of three begomoviruses causing TYLCD: tomato yellow leaf curl virus (TYLCV), tomato yellow leaf curl Sardinia virus (TYLCSV) and tomato yellow leaf curl Malaga virus (TYLCMaV, a natural recombinant between TYLCV and TYLCSV). Viral quasispecies resulting from inoculation of geminivirus infectious clones were examined in plants of susceptible tomato (ty-1/ty-1), heterozygous resistant tomato (Ty-1/ty-1), common bean, and the wild reservoir Solanum nigrum. Differences in virus fitness across hosts were observed while viral consensus sequences remained invariant. However, the complexity and heterogeneity of the quasispecies were high, especially in common bean and the wild host. Interestingly, the presence or absence of the Ty-1 allele in tomato did not lead to differences in begomovirus mutant spectra. However, the fitness decrease of TYLCSV and TYLCV in tomato at 45 dpi might be related to an increase in CP (Coat protein) mutation frequency. In Solanum nigrum the recombinant TYLCMaV, which showed lower fitness than TYLCSV, at 45 dpi actively explored Rep (Replication associated protein) ORF but not the overlapping C4. Our results underline the importance of begomovirus mutant spectra during infections. This is especially relevant in the wild reservoir of the viruses, which has the potential to maintain highly diverse mutant spectra without modifying their consensus sequences.

Viral Fitness Correlates with the Magnitude and Direction of the Perturbation Induced in the Host's Transcriptome: The Tobacco Etch Potyvirus-Tobacco Case Study

Determining the fitness of viral genotypes has become a standard practice in virology as it is essential to evaluate their evolutionary potential. Darwinian fitness, defined as the advantage of a given genotype with respect to a reference one, is a complex property that captures, in a single figure, differences in performance at every stage of viral infection. To what extent does viral fitness result from specific molecular interactions with host factors and regulatory networks during infection? Can we identify host genes in functional classes whose expression depends on viral fitness? Here, we compared the transcriptomes of tobacco plants infected with seven genotypes of tobacco etch potyvirus that differ in fitness. We found that the larger the fitness differences among genotypes, the more dissimilar the transcriptomic profiles are. Consistently, two different mutations, one in the viral RNA polymerase and another in the viral suppressor of RNA silencing, resulted in significantly similar gene expression profiles. Moreover, we identified host genes whose expression showed a significant correlation, positive or negative, with the virus' fitness. Differentially expressed genes which were positively correlated with viral fitness activate hormone- and RNA silencing-mediated pathways of plant defense. In contrast, those that were negatively correlated with fitness affect metabolism, reducing growth, and development. Overall, these results reveal the high information content of viral fitness and suggest its potential use to predict differences in genomic profiles of infected hosts.

High-Resolution Sequencing of Viral Populations during Early Simian Immunodeficiency Virus Infection Reveals Evolutionary Strategies for Rapid Escape from Emerging Env-Specific Antibody Responses

Primate lentiviruses, including the human and simian immunodeficiency viruses (HIV and SIV), produce infections marked by persistent, ongoing viral replication. This occurs despite the presence of virus-specific adaptive immune responses, including antibodies targeting the viral envelope glycoprotein (Env), and evolution of antibody-escape variants is a well-documented feature of lentiviral infection. Here, we examined the evolutionary dynamics of the SIV env gene during early infection (≤29 weeks postinfection) in a cohort of four SIV_mac251-infected rhesus macaques. We tracked env evolution during acute and early infection using frequent sampling and ultradeep sequencing of viral populations, capturing a transmission bottleneck and the subsequent reestablishment of Env diversity. A majority of changes in the gp120 subunit mapped to two short clusters, one in the first variable region (V1) and one in V4, while most changes in the gp41 subunit appeared in the cytoplasmic domain. Variation in V1 was dominated by short duplications and deletions of repetitive sequence, while variation in V4 was marked by short in-frame deletions and closely overlapping substitutions. The most common substitutions in both patches did not alter viral replicative fitness when tested using a highly sensitive, deep-sequencing-based competition assay. Our results, together with the observation that very similar or identical patterns of sequence evolution also occur in different macaque species infected with related but divergent strains of SIV, suggest that resistance to early, strain-specific anti-Env antibodies is the result of temporally and mutationally predictable pathways of escape that occur during the early stages of infection.IMPORTANCE The envelope glycoprotein (Env) of primate lentiviruses mediates entry by binding to host cell receptors followed by fusion of the viral membrane with the cell membrane. The exposure of Env complexes on the surface of the virion results in targeting by antibodies, leading to selection for virus escape mutations. We used the SIV/rhesus macaque model to track in vivo evolution of variation in Env during acute/early infection in animals with and without antibody responses to Env, uncovering remarkable variation in animals with antibody responses within weeks of infection. Using a deep-sequencing-based fitness assay, we found substitutions associated with antibody escape had little to no effect on inherent replicative capacity. The ability to readily propagate advantageous changes that incur little to no replicative fitness costs may be a mechanism to maintain continuous replication under constant immune selection, allowing the virus to persist for months to years in the infected host.

Incomplete inhibition of HIV infection results in more HIV infected lymph node cells by reducing cell death

HIV has been reported to be cytotoxic in vitro and in lymph node infection models. Using a computational approach, we found that partial inhibition of transmissions of multiple virions per cell could lead to increased numbers of live infected cells. If the number of viral DNA copies remains above one after inhibition, then eliminating the surplus viral copies reduces cell death. Using a cell line, we observed increased numbers of live infected cells when infection was partially inhibited with the antiretroviral efavirenz or neutralizing antibody. We then used efavirenz at concentrations reported in lymph nodes to inhibit lymph node infection by partially resistant HIV mutants. We observed more live infected lymph node cells, but with fewer HIV DNA copies per cell, relative to no drug. Hence, counterintuitively, limited attenuation of HIV transmission per cell may increase live infected cell numbers in environments where the force of infection is high.

The HIVvirus infects cells of the immune system. Once inside, it hijacks the cellular molecular machineries to make more copies of itself, which are then transmitted to new host cells. HIV eventually kills most cells it infects, either in the steps leading to the infection of the cell, or after the cell is already producing virus. HIV can spread between cells in two ways, known as cell-free or cell-to-cell. In the first, individual viruses are released from infected cells and move randomly through the body in the hope of finding new cells to infect. In the second, infected cells interact directly with uninfected cells. The second method is often much more successful at infecting new cells since they are exposed to multiple virus particles.

HIV infections can be controlled by using combinations of antiretroviral drugs, such as efavirenz, to prevent the virus from making more of itself. With a high enough dose, the drugs can in theory completely stop HIV infections, unless the virus becomes resistant to treatment. However, some patients continue to use these drugs even after the virus they are infected with develops resistance. It is not clear what effect taking ineffective, or partially effective, drugs has on how HIV progresses.

Using efavirenz, Jackson, Hunter et al. partially limited the spread of HIV between human cells grown in the laboratory. The experiments mirrored the situation where a partially resistant HIV strain spreads through the body. The results show that the success of cell-free infection is reduced as drug dose increases. Yet paradoxically, in cell-to-cell infection, the presence of drug caused more cells to become infected. This can be explained by the fact that, in cell-to-cell spread, each cell is exposed to multiple copies of the virus. The drug dose reduced the number of viral copies per cell without stopping the virus from infecting completely. The reduced number of viral copies per cell made it more likely that infected cells would survive the infection long enough to produce virus particles themselves.

Viruses that can kill cells, such as HIV, must balance the need to make more of themselves against the speed that they kill their host cell to maximize the number of infected cells. If transmission between cells is too effective and too many virus particles are delivered to the new cell, the virus may not manage to infect new hosts before killing the old ones. These findings highlight this delicate balance. They also indicate a potential issue in using drugs to treat partially resistant virus strains. Without care, these treatments could increase the number of infected cells in the body, potentially worsening the effects of living with HIV.

Natural Reassortants of Potentially Zoonotic Avian Influenza Viruses H5N1 and H9N2 from Egypt Display Distinct Pathogenic Phenotypes in Experimentally Infected Chickens and Ferrets

The cocirculation of zoonotic highly pathogenic avian influenza virus (HPAIV) of subtype H5N1 and avian influenza virus (AIV) of subtype H9N2 among poultry in Egypt for at least 6 years should render that country a hypothetical hot spot for the emergence of reassortant, phenotypically altered viruses, yet no reassortants have been detected in Egypt. The present investigations proved that reassortants of the Egyptian H5N1 clade 2.2.1.2 virus and H9N2 virus of the G1-B lineage can be generated by coamplification in embryonated chicken eggs. Reassortants were restricted to the H5N1 subtype and acquired between two and all six of the internal segments of the H9N2 virus. Five selected plaque-purified reassortant clones expressed a broad phenotypic spectrum both in vitro and in vivo Two groups of reassortants were characterized to have retarded growth characteristics in vitro compared to the H5N1 parent virus. One clone provoked reduced mortality in inoculated chickens, although the characteristics of a highly pathogenic phenotype were retained. Enhanced zoonotic properties were not predicted for any of these clones, and this prediction was confirmed by ferret inoculation experiments: neither the H5N1 parent virus nor two selected clones induced severe clinical symptoms or were transmitted to sentinel ferrets by contact. While the emergence of reassortants of Egyptian HPAIV of subtype H5N1 with internal gene segments of cocirculating H9N2 viruses is possible in principle, the spread of such viruses is expected to be governed by their fitness to outcompete the parental viruses in the field. The eventual spread of attenuated phenotypes, however, would negatively impact syndrome surveillance on poultry farms and might foster enzootic virus circulation.IMPORTANCE Despite almost 6 years of the continuous cocirculation of highly pathogenic avian influenza virus H5N1 and avian influenza virus H9N2 in poultry in Egypt, no reassortants of the two subtypes have been reported. Here, the principal compatibility of the two subtypes is shown by forcing the reassortment between copassaged H5N1 und H9N2 viruses in embryonated chicken eggs. The resulting reassortant viruses displayed a wide range of pathogenicity including attenuated phenotypes in chickens, but did not show enhanced zoonotic propensities in the ferret model.

Short Communication: Potential Risk of Replication-Competent Virus in HIV-1 Env-Pseudotyped Virus Preparations

Env-pseudotyped viruses are valuable reagents for studies of HIV-1 neutralizing antibodies. It is often assumed that all pseudovirus particles are capable of only a single round of infection, making them a safe alternative to work with live HIV-1. In this study, we show that some Env-pseudotyped virus preparations give rise to low levels of replication-competent virus. These levels did not compromise results in the TZM-bl neutralization assay; however, their presence highlights a need to adhere to the same level of biosafety when working with Env-pseudotyped viruses that are required for work with replication competent HIV-1.

Role of HLA Adaptation in HIV Evolution

Killing of HIV-infected cells by CD8⁺ T-cells imposes strong selection pressure on the virus toward escape. The HLA class I molecules that are successful in mediating some degree of control over the virus are those that tend to present epitopes in conserved regions of the proteome, such as in p24 Gag, in which escape also comes at a significant cost to viral replicative capacity (VRC). In some instances, compensatory mutations can fully correct for the fitness cost of such an escape variant; in others, correction is only partial. The consequences of these events within the HIV-infected host, and at the population level following transmission of escape variants, are discussed. The accumulation of escape mutants in populations over the course of the epidemic already shows instances of protective HLA molecules losing their impact, and in certain cases, a modest decline in HIV virulence in association with population-level increase in mutants that reduce VRC.

Resistance against Integrase Strand Transfer Inhibitors and Relevance to HIV Persistence

Drug resistance prevents the successful treatment of HIV-positive individuals by decreasing viral sensitivity to a drug or a class of drugs. In addition to transmitted resistant viruses, treatment-naïve individuals can be confronted with the problem of drug resistance through de novo emergence of such variants. Resistant viruses have been reported for every antiretroviral drug tested so far, including the integrase strand transfer inhibitors raltegravir, elvitegravir and dolutegravir. However, de novo resistant variants against dolutegravir have been found in treatment-experienced but not in treatment-naïve individuals, a characteristic that is unique amongst antiretroviral drugs. We review here the issue of drug resistance against integrase strand transfer inhibitors as well as both pre-clinical and clinical studies that have led to the identification of the R263K mutation in integrase as a signature resistance substitution for dolutegravir. We also discuss how the topic of drug resistance against integrase strand transfer inhibitors may have relevance in regard to the nature of the HIV reservoir and possible HIV curative strategies.

HIV-1 resistance to dolutegravir: update and new insights

Integrase strand transfer inhibitors (INSTIs) are the latest class of potent anti-HIV drugs. Currently, three INSTIs have been approved by the US Food and Drug Administration: raltegravir (RAL), elvitegravir (EVG) and dolutegravir (DTG). Resistance mutations to RAL and EVG emerge rapidly, and significant cross-resistance between these compounds has been documented. In addition, limited cross-resistance has been observed among DTG, a newer INSTI, and RAL and EVG even though clinical resistance to DTG, or mutations associated with DTG resistance in treatment-naïve patients, has not yet been observed. This review summarises progress in studies on understanding resistance to DTG, mechanisms of possible resistance to DTG, and reasons for the absence of DTG-associated resistance mutations when the drug has been used in first-line therapy.

A restriction enzyme based cloning method to assess the in vitro replication capacity of HIV-1 subtype C Gag-MJ4 chimeric viruses

The protective effect of many HLA class I alleles on HIV-1 pathogenesis and disease progression is, in part, attributed to their ability to target conserved portions of the HIV-1 genome that escape with difficulty. Sequence changes attributed to cellular immune pressure arise across the genome during infection, and if found within conserved regions of the genome such as Gag, can affect the ability of the virus to replicate in vitro. Transmission of HLA-linked polymorphisms in Gag to HLA-mismatched recipients has been associated with reduced set point viral loads. We hypothesized this may be due to a reduced replication capacity of the virus. Here we present a novel method for assessing the in vitro replication of HIV-1 as influenced by the gag gene isolated from acute time points from subtype C infected Zambians. This method uses restriction enzyme based cloning to insert the gag gene into a common subtype C HIV-1 proviral backbone, MJ4. This makes it more appropriate to the study of subtype C sequences than previous recombination based methods that have assessed the in vitro replication of chronically derived gag-pro sequences. Nevertheless, the protocol could be readily modified for studies of viruses from other subtypes. Moreover, this protocol details a robust and reproducible method for assessing the replication capacity of the Gag-MJ4 chimeric viruses on a CEM-based T cell line. This method was utilized for the study of Gag-MJ4 chimeric viruses derived from 149 subtype C acutely infected Zambians, and has allowed for the identification of residues in Gag that affect replication. More importantly, the implementation of this technique has facilitated a deeper understanding of how viral replication defines parameters of early HIV-1 pathogenesis such as set point viral load and longitudinal CD4+ T cell decline.

What if HIV were unable to develop resistance against a new therapeutic agent?

BACKGROUND

The HIV integrase inhibitor, Dolutegravir (DTG), was recently approved by the Food and Drug Administration in the United States and is the only HIV drug that has not selected for resistance mutations in the clinic when used as part of first-line therapy. This has led to speculation that DTG might have a higher genetic barrier for the development of drug resistance than the other compounds that are used in therapy.

DISCUSSION

In this Opinion article, we speculate that this is due to greatly diminished replication capacity on the part of viruses that might become resistant to DTG when the drug is used in initial therapy and that DTG might be able to be used in HIV prevention and eradication strategies. We also note that no compensatory mutation that might restore viral replication fitness to HIV in the aftermath of the appearance of a single drug resistance mutation has yet to be observed.

SUMMARY

DTG is a valuable addition to the anti-HIV armamentarium of drugs and its long-term utility may potentially exceed its obvious use in treatment of HIV disease.

In vitro phenotypes to elvitegravir and dolutegravir in primary macrophages and lymphocytes of clonal recombinant viral variants selected in patients failing raltegravir

OBJECTIVES

The cross-resistance profiles of elvitegravir and dolutegravir on raltegravir-resistant variants is still controversial or not available in macrophages and lack extensive evaluations on wide panels of clonal variants. Thus, a complete evaluation in parallel with all currently available integrase inhibitors (INIs) was performed.

METHODS

The integrase coding region was RT-PCR-amplified from patient-derived plasma samples and cloned into an HIV-1 molecular clone lacking the integrase region. Twenty recombinant viruses bearing mutations to all primary pathways of resistance to raltegravir were phenotypically evaluated with each integrase inhibitor in freshly purified CD4+ T cells or monocyte-derived macrophages.

RESULTS

Y143R single mutants conferred a higher level of raltegravir resistance in macrophages [fold change (FC) 47.7-60.24] compared with CD4+ T cells (FC 9.55-11.56). All other combinations had similar effects on viral susceptibility to raltegravir in both cell types. Elvitegravir displayed a similar behaviour both in lymphocytes and macrophages with all the tested patterns. When compared with raltegravir, none to modest increases in resistance were observed for the Y143R/C pathways. Dolutegravir maintained its activity and cross-resistance profile in macrophages. Only Q148H/R variants had a reduced level of susceptibility (FC 5.48-18.64). No variations were observed for the Y143R/C (+/-T97A) or N155H variants.

CONCLUSIONS

All INIs showed comparable antiretroviral activity in both cell types even if single mutations were associated with a different level of susceptibility in vitro to raltegravir and elvitegravir in macrophages. In particular, dolutegravir was capable of inhibiting with similar potency infection of raltegravir-resistant variants with Y143 or N155 pathways in both HIV-1 major cell reservoirs.