Effect of alternating passage on adaptation of sindbis virus to vertebrate and invertebrate cells

Mosquito-borne alphaviruses, which replicate alternately and obligately in mosquitoes and vertebrates, appear to experience lower rates of evolution than do many RNA viruses that replicate solely in vertebrates. This genetic stability is hypothesized to result from the alternating host cycle, which constrains evolution by imposing compromise fitness solutions in each host. To test this hypothesis, Sindbis virus was passaged serially, either in one cell type to eliminate host alteration or alternately between vertebrate (BHK) and mosquito (C6/36) cells. Following 20 to 50 serial passages, mutations were identified and changes in fitness were assessed using competition assays against genetically marked, surrogate parent viruses. Specialized viruses passaged in a single cell exhibited more mutations and amino acid changes per passage than those passaged alternately. Single host-adapted viruses exhibited fitness gains in the cells in which they specialized but fitness losses in the bypassed cell type. Most but not all viruses passaged alternately experienced lesser fitness gains than specialized viruses, with fewer mutations per passage. Clonal populations derived from alternately passaged viruses also exhibited adaptation to both cell lines, indicating that polymorphic populations are not required for simultaneous fitness gains in vertebrate and mosquito cells. Nearly all passaged viruses acquired Arg or Lys substitutions in the E2 envelope glycoprotein, but enhanced binding was only detected for BHK cells. These results support the hypothesis that arbovirus evolution may be constrained by alternating host transmission cycles, but they indicate a surprising ability for simultaneous adaptation to highly divergent cell types by combinations of mutations in single genomes.     

Consequences of host adaptation for performance of vesicular stomatitis virus in novel thermal environments

Although it is widely assumed that the selective advantage of niche specialization drives species biodiversity, some theory suggests that generalists are favored over specialists when environments change unexpectedly. But this idea is rarely tested empirically, and its relevance is unknown for microparasites such as RNA viruses. Due to their small genome sizes pleiotropy is not uncommon in RNA viruses. Therefore, the genetic architectures underlying generalist traits may be indirectly molded by selection to better prepare generalist organisms for growth in new environments. Previously, vesicular stomatitis viruses were evolved to specialize on a single host, or to generalize on multiple hosts. Here we test whether virus generalists arising in the context of host adaptation also perform differently than specialists when viruses grow at novel temperatures. We compared thermal reaction norms of performance, within and among groups of viral specialists and generalists. Results showed that host adaptation was consequential for some fitness traits at novel temperatures due to modification of pleiotropic viral genes. Contrary to theoretical predictions, host generalists were selectively disadvantaged at extreme cool and warm environments. Multi-host adaptation may compromise the evolved thermostability of viral proteins, creating a cost of host generalization when viruses replicate at extreme temperatures.     

Differential RNA Sequence Requirement for Dengue Virus Replication in Mosquito and Mammalian Cells

Dengue virus cycles between mosquitoes and humans. Each host provides a different environment for viral replication, imposing different selective pressures. We identified a sequence in the dengue virus genome that is essential for viral replication in mosquito cells but not in mammalian cells. This sequence is located at the viral 3′ untranslated region and folds into a small hairpin structure. A systematic mutational analysis using dengue virus infectious clones and reporter viruses allowed the determination of two putative functions in this cis-acting RNA motif, one linked to the structure and the other linked to the nucleotide sequence. We found that single substitutions that did not alter the hairpin structure did not affect dengue virus replication in mammalian cells but abolished replication in mosquito cells. This is the first sequence identified in a flavivirus genome that is exclusively required for viral replication in insect cells.     

The Interactomes of Influenza Virus NS1 and NS2 Proteins Identify New Host Factors and Provide Insights for ADAR1 Playing a Supportive Role in Virus Replication

Influenza A NS1 and NS2 proteins are encoded by the RNA segment 8 of the viral genome. NS1 is a multifunctional protein and a virulence factor while NS2 is involved in nuclear export of viral ribonucleoprotein complexes. A yeast two-hybrid screening strategy was used to identify host factors supporting NS1 and NS2 functions. More than 560 interactions between 79 cellular proteins and NS1 and NS2 proteins from 9 different influenza virus strains have been identified. These interacting proteins are potentially involved in each step of the infectious process and their contribution to viral replication was tested by RNA interference. Validation of the relevance of these host cell proteins for the viral replication cycle revealed that 7 of the 79 NS1 and/or NS2-interacting proteins positively or negatively controlled virus replication. One of the main factors targeted by NS1 of all virus strains was double-stranded RNA binding domain protein family. In particular, adenosine deaminase acting on RNA 1 (ADAR1) appeared as a pro-viral host factor whose expression is necessary for optimal viral protein synthesis and replication. Surprisingly, ADAR1 also appeared as a pro-viral host factor for dengue virus replication and directly interacted with the viral NS3 protein. ADAR1 editing activity was enhanced by both viruses through dengue virus NS3 and influenza virus NS1 proteins, suggesting a similar virus-host co-evolution.     

Long-range RNA-RNA interactions circularize the dengue virus genome

Secondary and tertiary RNA structures present in viral RNA genomes play essential regulatory roles during translation, RNA replication, and assembly of new viral particles. In the case of flaviviruses, RNA-RNA interactions between the 5' and 3' ends of the genome have been proposed to be required for RNA replication. We found that two RNA elements present at the ends of the dengue virus genome interact in vitro with high affinity. Visualization of individual molecules by atomic force microscopy revealed that physical interaction between these RNA elements results in cyclization of the viral RNA. Using RNA binding assays, we found that the putative cyclization sequences, known as 5' and 3' CS, present in all mosquito-borne flaviviruses, were necessary but not sufficient for RNA-RNA interaction. Additional sequences present at the 5' and 3' untranslated regions of the viral RNA were also required for RNA-RNA complex formation. We named these sequences 5' and 3' UAR (upstream AUG region). In order to investigate the functional role of 5'-3' UAR complementarity, these sequences were mutated either separately, to destroy base pairing, or simultaneously, to restore complementarity in the context of full-length dengue virus RNA. Nonviable viruses were recovered after transfection of dengue virus RNA carrying mutations either at the 5' or 3' UAR, while the RNA containing the compensatory mutations was able to replicate. Since sequence complementarity between the ends of the genome is required for dengue virus viability, we propose that cyclization of the RNA is a required conformation for viral replication.     

The role of evolutionary intermediates in the host adaptation of canine parvovirus

The adaptation of viruses to new hosts is a poorly understood process likely involving a variety of viral structures and functions that allow efficient replication and spread. Canine parvovirus (CPV) emerged in the late 1970s as a host-range variant of a virus related to feline panleukopenia virus (FPV). Within a few years of its emergence in dogs, there was a worldwide replacement of the initial virus strain (CPV type 2) by a variant (CPV type 2a) characterized by four amino acid differences in the capsid protein. However, the evolutionary processes that underlie the acquisition of these four mutations, as well as their effects on viral fitness, both singly and in combination, are still uncertain. Using a comprehensive experimental analysis of multiple intermediate mutational combinations, we show that these four capsid mutations act in concert to alter antigenicity, cell receptor binding, and relative in vitro growth in feline cells. Hence, host adaptation involved complex interactions among both surface-exposed and buried capsid mutations that together altered cell infection and immune escape properties of the viruses. Notably, most intermediate viral genotypes containing different combinations of the four key amino acids possessed markedly lower fitness than the wild-type viruses.     

Structure-altering mutations of the SARS-CoV-2 frameshifting RNA element

With the rapid rate of COVID-19 infections and deaths, treatments and cures besides hand washing, social distancing, masks, isolation, and quarantines are urgently needed. The treatments and vaccines rely on the basic biophysics of the complex viral apparatus. Although proteins are serving as main drug and vaccine targets, therapeutic approaches targeting the 30,000 nucleotide RNA viral genome form important complementary approaches. Indeed, the high conservation of the viral genome, its close evolutionary relationship to other viruses, and the rise of gene editing and RNA-based vaccines all argue for a focus on the RNA agent itself. One of the key steps in the viral replication cycle inside host cells is the ribosomal frameshifting required for translation of overlapping open reading frames. The RNA frameshifting element (FSE), one of three highly conserved regions of coronaviruses, is believed to include a pseudoknot considered essential for this ribosomal switching. In this work, we apply our graph-theory-based framework for representing RNA secondary structures, “RAG (or RNA-As-Graphs),” to alter key structural features of the FSE of the SARS-CoV-2 virus. Specifically, using RAG machinery of genetic algorithms for inverse folding adapted for RNA structures with pseudoknots, we computationally predict minimal mutations that destroy a structurally important stem and/or the pseudoknot of the FSE, potentially dismantling the virus against translation of the polyproteins. Our microsecond molecular dynamics simulations of mutant structures indicate relatively stable secondary structures. These findings not only advance our computational design of RNAs containing pseudoknots, they pinpoint key residues of the SARS-CoV-2 virus as targets for antiviral drugs and gene editing approaches.     

Host-Specific Parvovirus Evolution in Nature Is Recapitulated by In Vitro Adaptation to Different Carnivore Species

Canine parvovirus (CPV) emerged as a new pandemic pathogen of dogs in the 1970s and is closely related to feline panleukopenia virus (FPV), a parvovirus of cats and related carnivores. Although both viruses have wide host ranges, analysis of viral sequences recovered from different wild carnivore species, as shown here, demonstrated that >95% were derived from CPV-like viruses, suggesting that CPV is dominant in sylvatic cycles. Many viral sequences showed host-specific mutations in their capsid proteins, which were often close to sites known to control binding to the transferrin receptor (TfR), the host receptor for these carnivore parvoviruses, and which exhibited frequent parallel evolution. To further examine the process of host adaptation, we passaged parvoviruses with alternative backgrounds in cells from different carnivore hosts. Specific mutations were selected in several viruses and these differed depending on both the background of the virus and the host cells in which they were passaged. Strikingly, these in vitro mutations recapitulated many specific changes seen in viruses from natural populations, strongly suggesting they are host adaptive, and which were shown to result in fitness advantages over their parental virus. Comparison of the sequences of the transferrin receptors of the different carnivore species demonstrated that many mutations occurred in and around the apical domain where the virus binds, indicating that viral variants were likely selected through their fit to receptor structures. Some of the viruses accumulated high levels of variation upon passage in alternative hosts, while others could infect multiple different hosts with no or only a few additional mutations. Overall, these studies demonstrate that the evolutionary history of a virus, including how long it has been circulating and in which hosts, as well as its phylogenetic background, has a profound effect on determining viral host range.     

Generalized selection to overcome innate immunity selects for host breadth in an RNA virus

Virus‐host coevolution has selected for generalized host defense against viruses, exemplified by interferon production/signaling and other innate immune function in eukaryotes such as humans. Although cell‐surface binding primarily limits virus infection success, generalized adaptation to counteract innate immunity across disparate hosts may contribute to RNA virus emergence potential. We examined this idea using vesicular stomatitis virus (VSV) populations previously evolved on strictly immune‐deficient (HeLa) cells, strictly immune competent (MDCK) cells, or on alternating deficient/competent cells. By measuring viral fitness in unselected human cancer cells of differing innate immunity, we confirmed that HeLa‐adapted populations were specialized for innate immune‐deficient hosts, whereas MDCK‐adapted populations were relatively more generalized for fitness on hosts of differing innate immune capacity and of different species origin. We also confirmed that HeLa‐evolved populations maintained fitness in immune‐deficient nonhuman primate cells. These results suggest that innate immunity is more prominent than host species in determining viral fitness at the host‐cell level. Finally, our prediction was inexact that selection on alternating deficient/competent hosts should produce innate viral generalists. Rather, fitness differences among alternating host‐evolved VSV populations indicated variable capacities to evade innate immunity. Our results suggest that the evolutionary history of innate immune selection can affect whether RNA viruses evolve greater host‐breadth.     

Mutagenesis of the dengue virus type 2 NS5 methyltransferase domain

The Flavivirus NS5 protein possesses both (guanine-N7)-methyltransferase and nucleoside-2'-O methyltransferase activities required for sequential methylation of the cap structure present at the 5' end of the Flavivirus RNA genome. Seventeen mutations were introduced into the dengue virus type 2 NS5 methyltransferase domain, targeting amino acids either predicted to be directly involved in S-adenosyl-l-methionine binding or important for NS5 conformation and/or charged interactions. The effects of the mutations on (i) (guanine-N7)-methyltransferase and nucleoside-2'-O methyltransferase activities using biochemical assays based on a bacterially expressed NS5 methyltransferase domain and (ii) viral replication using a dengue virus type 2 infectious cDNA clone were examined. Clustered mutations targeting the S-adenosyl-l-methionine binding pocket or an active site residue abolished both methyltransferase activities and viral replication, demonstrating that both methyltransferase activities utilize a single S-adenosyl-l-methionine binding pocket. Substitutions to single amino acids binding S-adenosyl-l-methionine decreased both methyltransferase activities by varying amounts. However, viruses that replicated at wild type levels could be recovered with mutations that reduced both activities by >75%, suggesting that only a threshold level of methyltransferase activity was required for virus replication in vivo. Mutation of residues outside of regions directly involved in S-adenosyl-l-methionine binding or catalysis also affected methyltransferase activity and virus replication. The recovery of viruses containing compensatory second site mutations in the NS5 and NS3 proteins identified regions of the methyltransferase domain important for overall stability of the protein or likely to play a role in virus replication distinct from that of cap methylation.     

Pleiotropic costs of niche expansion in the RNA bacteriophage phi 6

Natural and experimental systems have failed to universally demonstrate a trade-off between generalism and specialism. When a trade-off does occur it is difficult to attribute its cause to antagonistic pleiotropy without dissecting the genetic basis of adaptation, and few previous experiments provide these genetic data. Here we investigate the evolution of expanded host range (generalism) in the RNA virus phi6, an experimental model system allowing adaptive mutations to be readily identified. We isolated 10 spontaneous host range mutants on each of three novel Pseudomonas hosts and determined whether these mutations imposed fitness costs on the standard laboratory host. Sequencing revealed that each mutant had one of nine nonsynonymous mutations in the phi6 gene P3, important in host attachment. Seven of these nine mutations were costly on the original host, confirming the existence of antagonistic pleiotropy. In addition to this genetically imposed cost, we identified an epigenetic cost of generalism that occurs when phage transition between host types. Our results confirm the existence in phi6 of two costs of generalism, genetic and environmental, but they also indicate that the cost is not always large. The possibility for cost-free niche expansion implies that varied ecological conditions may favor host shifts in RNA viruses.     

Secondary structure of the human T-cell leukemia virus type 1 rex-responsive element is essential for rex regulation of RNA processing and transport of unspliced RNAs.

Rex protein of human T-cell leukemia virus type 1 (HTLV-1) induces cytoplasmic expression of unspliced gag/pol mRNA and singly spliced env mRNA and thus is essential for replication of the virus. This regulation requires a cis-acting rex-responsive element (RXE), located in the 3' region of the viral RNA. By external deletion, we have identified RXE composed of 205 nucleotides. The secondary structure of RXE was confirmed by studies on its susceptibility to nuclease digestions to consist of four stem-loops and a long stretch of stem structure. Substitution and deletion mutations revealed that two regions of the stem-loops and their secondary structures are essential for rex regulation. Similar secondary structures were found in the corresponding regions of HTLV-2, bovine leukemia virus and human immunodeficiency virus. Furthermore, a sequence of 11 nucleotides in the RXE was found to be conserved in the secondary structures of HTLV-1, HTLV-2, and bovine leukemia virus. These observations suggest that the secondary structure as well as the conserved sequence may be important in expression of unspliced RNA even with diverged sequences as observed in these viruses.     

Evolutionary genomics of host-use in bifurcating demes of RNA virus phi-6

ABSTRACT: BACKGROUND: Viruses are exceedingly diverse in their evolved strategies to manipulate hosts for viral replication. However, despite these differences, most virus populations will occasionally experience two commonly-encountered challenges: growth in variable host environments, and growth under fluctuating population sizes. We used the segmented RNA bacteriophage [GREEK PHI SYMBOL]6 as a model for studying the evolutionary genomics of virus adaptation in the face of host switches and parametrically varying population sizes. To do so, we created a bifurcating deme structure that reflected lineage splitting in natural populations, allowing us to test whether phylogenetic algorithms could accurately resolve this 'known phylogeny'. The resulting tree yielded 32 clones at the tips and internal nodes; these strains were fully sequenced and measured for phenotypic changes in selected traits (fitness on original and novel hosts). RESULTS: We observed that RNA segment size was negatively correlated with the extent of molecular change in the imposed treatments; molecular substitutions tended to cluster on the Small and Medium RNA chromosomes of the virus, and not on the Large segment. Our study yielded a very large molecular and phenotypic dataset, fostering possible inferences on genotype-phenotype associations. Using further experimental evolution, we confirmed an inference on the unanticipated role of an allelic switch in a viral assembly protein, which governed viral performance across host environments. CONCLUSIONS: Our study demonstrated that varying complexities can be simultaneously incorporated into experimental evolution, to examine the combined effects of population size, and adaptation in novel environments. The imposed bifurcating structure revealed that some methods for phylogenetic reconstruction failed to resolve the true phylogeny, owing to a paucity of molecular substitutions separating the RNA viruses that evolved in our study.     

Thiouracil Cross-Linking Mass Spectrometry: a Cell-Based Method To Identify Host Factors Involved in Viral Amplification

Eukaryotic RNA viruses are known to utilize host factors; however, the identity of these factors and their role in the virus life cycle remain largely undefined. Here, we report a method to identify proteins bound to the viral RNA during amplification in cell culture: thiouracil cross-linking mass spectrometry (TUX-MS). TUX-MS relies on incorporation of a zero-distance cross-linker into the viral RNA during infection. Proteins bound to viral RNA are cross-linked prior to cell lysis, purified, and identified using mass spectrometry. Using the TUX-MS method, an unbiased screen for poliovirus (PV) host factors was conducted. All host and viral proteins that are known to interact with the poliovirus RNA were identified. In addition, TUX-MS identified an additional 66 host proteins that have not been previously described in poliovirus amplification. From these candidates, eight were selected and validated. Furthermore, we demonstrate that small interfering RNA (siRNA)-mediated knockdown of two of these uncharacterized host factors results in either a decrease in copy number of positive-stranded RNA or a decrease in PV translation. These data demonstrate that TUX-MS is a robust, unbiased method to identify previously unknown host cell factors that influence virus growth. This method is broadly applicable to a range of RNA viruses, such as flaviviruses, alphaviruses, picornaviruses, bunyaviruses, and coronaviruses.     

Group Selection and Contribution of Minority Variants during Virus Adaptation Determines Virus Fitness and Phenotype

Understanding how a pathogen colonizes and adapts to a new host environment is a primary aim in studying emerging infectious diseases. Adaptive mutations arise among the thousands of variants generated during RNA virus infection, and identifying these variants will shed light onto how changes in tropism and species jumps can occur. Here, we adapted Coxsackie virus B3 to a highly permissive and less permissive environment. Using deep sequencing and bioinformatics, we identified a multi-step adaptive process to adaptation involving residues in the receptor footprints that correlated with receptor availability and with increase in virus fitness in an environment-specific manner. We show that adaptation occurs by selection of a dominant mutation followed by group selection of minority variants that together, confer the fitness increase observed in the population, rather than selection of a single dominant genotype.     

Resistance to extinction of low fitness virus subjected to plaque-to-plaque transfers: diversification by mutation clustering.

Plaque-to-plaque transfers of RNA viruses lead to accumulation of mutations and fitness decrease. To test whether continuing plaque-to-plaque transfers would lead to viral extinction, we have subjected several low fitness foot-and-mouth disease virus (FMDV) clones to up to 130 successive plaque transfers, and have analyzed the evolution of plaque titers and genomic nucleotide sequences. No case of viral extinction could be documented. Some low fitness clones that posses an internal poly(A) tract evaded extinction by modifying the length or base composition of the poly(A) tract. The comparison of entire genomic sequences of FMDV clones at increasing plaque transfer number revealed that mutations accumulated at a uniform rate, and that they were distributed unevenly along the genome. Clusters of mutations were identified at different genomic sites in two plaque transfer lineages. Mutation clustering appears to occur stochastically and could not be related to fixation of compensatory mutations. The results document resistance of viral clones to extinction, and suggest that mutation clustering may be a mechanism of genetic diversification of low fitness virus.     

Isolation of Fidelity Variants of RNA Viruses and Characterization of Virus Mutation Frequency

RNA viruses use RNA dependent RNA polymerases to replicate their genomes. The intrinsically high error rate of these enzymes is a large contributor to the generation of extreme population diversity that facilitates virus adaptation and evolution. Increasing evidence shows that the intrinsic error rates, and the resulting mutation frequencies, of RNA viruses can be modulated by subtle amino acid changes to the viral polymerase. Although biochemical assays exist for some viral RNA polymerases that permit quantitative measure of incorporation fidelity, here we describe a simple method of measuring mutation frequencies of RNA viruses that has proven to be as accurate as biochemical approaches in identifying fidelity altering mutations. The approach uses conventional virological and sequencing techniques that can be performed in most biology laboratories. Based on our experience with a number of different viruses, we have identified the key steps that must be optimized to increase the likelihood of isolating fidelity variants and generating data of statistical significance. The isolation and characterization of fidelity altering mutations can provide new insights into polymerase structure and function^1-3. Furthermore, these fidelity variants can be useful tools in characterizing mechanisms of virus adaptation and evolution^4-7.     

Emergence of the Virulence-Associated PB2 E627K Substitution in a Fatal Human Case of Highly Pathogenic Avian Influenza Virus A(H7N7) Infection as Determined by Illumina Ultra-Deep Sequencing

Avian influenza viruses are capable of crossing the species barrier and infecting humans. Although evidence of human-to-human transmission of avian influenza viruses to date is limited, evolution of variants toward more-efficient human-to-human transmission could result in a new influenza virus pandemic. In both the avian influenza A(H5N1) and the recently emerging avian influenza A(H7N9) viruses, the polymerase basic 2 protein (PB2) E627K mutation appears to be of key importance for human adaptation. During a large influenza A(H7N7) virus outbreak in the Netherlands in 2003, the A(H7N7) virus isolated from a fatal human case contained the PB2 E627K mutation as well as a hemagglutinin (HA) K416R mutation. In this study, we aimed to investigate whether these mutations occurred in the avian or the human host by Illumina Ultra-Deep sequencing of three previously uninvestigated clinical samples obtained from the fatal case. In addition, we investigated three chicken samples, two of which were obtained from the source farm. Results showed that the PB2 E627K mutation was not present in any of the chicken samples tested. Surprisingly, the avian samples were characterized by the presence of influenza virus defective RNA segments, suggestive for the synthesis of defective interfering viruses during infection in poultry. In the human samples, the PB2 E627K mutation was identified with increasing frequency during infection. Our results strongly suggest that human adaptation marker PB2 E627K has emerged during virus infection of a single human host, emphasizing the importance of reducing human exposure to avian influenza viruses to reduce the likelihood of viral adaptation to humans.