A giant virus infecting the amoeboflagellate Naegleria.

Giant viruses (Nucleocytoviricota) are significant lethality agents of various eukaryotic hosts. Although metagenomics indicates their ubiquitous distribution, available giant virus isolates are restricted to a very small number of protist and algal hosts. Here we report on the first viral isolate that replicates in the amoeboflagellate Naegleria. This genus comprises the notorious human pathogen Naegleria fowleri, the causative agent of the rare but fatal primary amoebic meningoencephalitis. We have elucidated the structure and infection cycle of this giant virus, Catovirus naegleriensis (a.k.a. Naegleriavirus, NiV), and show its unique adaptations to its Naegleria host using fluorescence in situ hybridization, electron microscopy, genomics, and proteomics. Naegleriavirus is only the fourth isolate of the highly diverse subfamily Klosneuvirinae, and like its relatives the NiV genome contains a large number of translation genes, but lacks transfer RNAs (tRNAs). NiV has acquired genes from its Naegleria host, which code for heat shock proteins and apoptosis inhibiting factors, presumably for host interactions. Notably, NiV infection was lethal to all Naegleria species tested, including the human pathogen N. fowleri. This study expands our experimental framework for investigating giant viruses and may help to better understand the basic biology of the human pathogen N. fowleri.

Transcriptomic, proteomic, and functional consequences of codon usage bias in human cells during heterologous gene expression.

Differences in codon frequency between genomes, genes, or positions along a gene, modulate transcription and translation efficiency, leading to phenotypic and functional differences. Here, we present a multiscale analysis of the effects of synonymous codon recoding during heterologous gene expression in human cells, quantifying the phenotypic consequences of codon usage bias at different molecular and cellular levels, with an emphasis on translation elongation. Six synonymous versions of an antibiotic resistance gene were generated, fused to a fluorescent reporter, and independently expressed in HEK293 cells. Multiscale phenotype was analyzed by means of quantitative transcriptome and proteome assessment, as proxies for gene expression; cellular fluorescence, as a proxy for single-cell level expression; and real-time cell proliferation in absence or presence of antibiotic, as a proxy for the cell fitness. We show that differences in codon usage bias strongly impact the molecular and cellular phenotype: (i) they result in large differences in mRNA levels and protein levels, leading to differences of over 15 times in translation efficiency; (ii) they introduce unpredicted splicing events; (iii) they lead to reproducible phenotypic heterogeneity; and (iv) they lead to a trade-off between the benefit of antibiotic resistance and the burden of heterologous expression. In human cells in culture, codon usage bias modulates gene expression by modifying mRNA availability and suitability for translation, leading to differences in protein levels and eventually eliciting functional phenotypic changes.

Defensive symbiosis against giant viruses in amoebae.

Protists are important regulators of microbial communities and key components in food webs with impact on nutrient cycling and ecosystem functioning. In turn, their activity is shaped by diverse intracellular parasites, including bacterial symbionts and viruses. Yet, bacteria-virus interactions within protists are poorly understood. Here, we studied the role of bacterial symbionts of free-living amoebae in the establishment of infections with nucleocytoplasmic large DNA viruses (Nucleocytoviricota). To investigate these interactions in a system that would also be relevant in nature, we first isolated and characterized a giant virus (Viennavirus, family Marseilleviridae) and a sympatric potential Acanthamoeba host infected with bacterial symbionts. Subsequently, coinfection experiments were carried out, using the fresh environmental isolates as well as additional amoeba laboratory strains. Employing fluorescence in situ hybridization and qPCR, we show that the bacterial symbiont, identified as Parachlamydia acanthamoebae, represses the replication of the sympatric Viennavirus in both recent environmental isolates as well as Acanthamoeba laboratory strains. In the presence of the symbiont, virions are still taken up, but viral factory maturation is inhibited, leading to survival of the amoeba host. The symbiont also suppressed the replication of the more complex Acanthamoeba polyphaga mimivirus and Tupanvirus deep ocean (Mimiviridae). Our work provides an example of an intracellular bacterial symbiont protecting a protist host against virus infections. The impact of virus-symbiont interactions on microbial population dynamics and eventually ecosystem processes requires further attention.

Genomic and phylogenetic characterization of ChPV2, a novel goat PV closely related to the Xi-PV1 species infecting bovines.

BACKGROUND: Papillomaviruses (PVs) infecting artiodactyls are very diverse, and only second in number to PVs infecting primates. PVs associated to lesions in economically important ruminant species have been isolated from cattle and sheep. METHODS: Potential PV DNA from teat lesions of a Damascus goat was isolated, cloned and sequenced. The PV genome was analyzed using bioinformatics approaches to detect open reading frames and to predict potential features of encoded proteins as well as putative regulatory elements. Sequence comparison and phylogenetic analyses using the concatenated E1E2L2L1 nucleotide and amino acid alignments was used to reveal the relationship of the new PV to the known PV diversity and its closest relevants. RESULTS: We isolated and characterized the full-genome of novel Capra hircus papillomavirus. We identified the E6, E7, E1, E2, L2, L1 open reading frames with protein coding potential and putative active elements in the ChPV2 proteins and putative regulatory genome elements. Sequence similarities of L1 and phylogenetic analyses using concatenated E1E2L2L1 nucleotide and amino acid alignments suggest the classification as a new PV type designated ChPV2 with a phylogenetic position within the XiPV genus, basal to the XiPV1 species. ChPV2 is not closely related to ChPV1, the other known goat PV isolated from healthy skin, although both of them belong confidently into a clade composed of PVs infecting cervids and bovids. Interestingly, ChPV2 contains an E6 open reading frame whereas all closely related PVs do not CONCLUSION: ChPV2 is a novel goat PV closely related to the Xi-PV1 species infecting bovines. Phylogenetic relationships and genome architecture of ChPV2 and closely related PV types suggest at least two independent E6 losses within the XiPV clade.

Papillomaviruses infecting cetaceans exhibit signs of genome adaptation following a recombination event.

Papillomaviruses (PVs) have evolved through a complex evolutionary scenario where virus-host co-evolution alone is not enough to explain the phenotypic and genotypic PV diversity observed today. Other evolutionary processes, such as host switch and recombination, also appear to play an important role in PV evolution. In this study, we have examined the genomic impact of a recombination event between distantly related PVs infecting Cetartiodactyla (even-toed ungulates and cetaceans). Our phylogenetic analyses suggest that one single recombination was responsible for the generation of extant 'chimeric' PV genomes infecting cetaceans. By correlating the phylogenetic relationships to the genomic content, we observed important differences between the recombinant and non-recombinant cetartiodactyle PV genomes. Notably, recombinant PVs contain a unique set of conserved motifs in the upstream regulatory region (URR). We interpret these regulatory changes as an adaptive response to drastic changes in the PV genome. In terms of codon usage preferences (CUPrefs), we did not detect any particular differences between orthologous open reading frames in recombinant and non-recombinant PVs. Instead, our results are in line with previous observations suggesting that CUPrefs in PVs are rather linked to gene expression patterns as well as to gene function. We show that the non-coding URR of PVs infecting cetaceans, the central regulatory element in these viruses, exhibits signs of adaptation following a recombination event. Our results suggest that also in PVs, the evolution of gene regulation can play an important role in speciation and adaptation to novel environments.

Detection of Alpha, Beta, Gamma, and Unclassified Human Papillomaviruses in Cervical Cancer Samples From Mexican Women.

Background: Cervical cancer (CC) is associated to high-risk human papillomavirus (HPV) infections, for this reason it is crucial to have sensitive and accurate HPV diagnostic tests. To date, most research is focused on HPVs within the Alphapapillomavirus (α-PVs) genus and little attention has been paid to cervical infections with other HPV genotypes, like those of the Betapapillomavirus (ß-PVs) and Gammapapillomavirus (γ-PVs) genera. The aim of this study was to determine the HPV genotypes from different genera in women with CC using Next-Generation Sequencing (NGS). Methods: The study comprised 48 HPV positive CC samples evaluated with the Linear Array HPV Genotyping test and individually sequenced by 454 NGS using PGMY09/11 and FAP primers. To determine the HPV genotypes present in each sample, the obtained sequences were compared with all HPV L1 gene reference sequences from the Papillomavirus Episteme database (PaVE). Moreover, 50 HPV positive low-grade cervical lesion samples individually genotyped with NGS were also included to determine the genotypes present preferentially in CC patients. Results: Among the 48 CC samples, 68.75% consisted of multiple HPV infections, 51 different genotypes were detected, of which 7 are still unclassified, 28 belong to α-PVs (6, 11, 16, 18, 26, 30, 33, 35, 39, 42, 43, 44, 45, 51, 52, 53, 54, 59, 62, 66, 68, 69, 70, 71, 74, 81, 102, 114), 10 to ß-PVs (5, 12, 21, 37, 38b, 47, 80, 107, 118, 122), and 6 to γ-PVs (101, 103, 123, 135, 147, 214). Among them, HPV16 was the most prevalent genotype (54.2%), followed by HPV18 (16.7%), HPV38b (14.6%), and HPVs 52/62/80 (8.3%). Some genotypes were exclusively found in CC when compared with Cervical Intraepithelial Neoplasia grade 1 (CIN1) samples, such as HPVs 5, 18, 38b, 107, 122, FA39, FA116, mSK_120, and mSK_136. Conclusions: This work demonstrates the great diversity of HPV genotypes detected by combining PGMY and FAP primers with NGS in cervical swabs. The relatively high attribution of ß- and γ- PVs in CC samples suggest their possible role as carcinogenic cofactors, but deeper studies need to be performed to determine if they have transforming properties and the significance of HPV-coinfections.

Erratum: On the stability of sequences inserted into viral genomes.

[This corrects the article DOI: 10.1093/ve/vez045.].

On the stability of sequences inserted into viral genomes.

Viruses are widely used as vectors for heterologous gene expression in cultured cells or natural hosts, and therefore a large number of viruses with exogenous sequences inserted into their genomes have been engineered. Many of these engineered viruses are viable and express heterologous proteins at high levels, but the inserted sequences often prove to be unstable over time and are rapidly lost, limiting heterologous protein expression. Although virologists are aware that inserted sequences can be unstable, processes leading to insert instability are rarely considered from an evolutionary perspective. Here, we review experimental work on the stability of inserted sequences over a broad range of viruses, and we present some theoretical considerations concerning insert stability. Different virus genome organizations strongly impact insert stability, and factors such as the position of insertion can have a strong effect. In addition, we argue that insert stability not only depends on the characteristics of a particular genome, but that it will also depend on the host environment and the demography of a virus population. The interplay between all factors affecting stability is complex, which makes it challenging to develop a general model to predict the stability of genomic insertions. We highlight key questions and future directions, finding that insert stability is a surprisingly complex problem and that there is need for mechanism-based, predictive models. Combining theoretical models with experimental tests for stability under varying conditions can lead to improved engineering of viral modified genomes, which is a valuable tool for understanding genome evolution as well as for biotechnological applications, such as gene therapy.

Genome Plasticity in Papillomaviruses and De Novo Emergence of E5 Oncogenes.

The clinical presentations of papillomavirus (PV) infections come in many different flavors. While most PVs are part of a healthy skin microbiota and are not associated to physical lesions, other PVs cause benign lesions, and only a handful of PVs are associated to malignant transformations linked to the specific activities of the E5, E6, and E7 oncogenes. The functions and origin of E5 remain to be elucidated. These E5 open reading frames (ORFs) are present in the genomes of a few polyphyletic PV lineages, located between the early and the late viral gene cassettes. We have computationally assessed whether these E5 ORFs have a common origin and whether they display the properties of a genuine gene. Our results suggest that during the evolution of Papillomaviridae, at least four events lead to the presence of a long noncoding DNA stretch between the E2 and the L2 genes. In three of these events, the novel regions evolved coding capacity, becoming the extant E5 ORFs. We then focused on the evolution of the E5 genes in AlphaPVs infecting primates. The sharp match between the type of E5 protein encoded in AlphaPVs and the infection phenotype (cutaneous warts, genital warts, or anogenital cancers) supports the role of E5 in the differential oncogenic potential of these PVs. In our analyses, the best-supported scenario is that the five types of extant E5 proteins within the AlphaPV genomes may not have a common ancestor. However, the chemical similarities between E5s regarding amino acid composition prevent us from confidently rejecting the model of a common origin. Our evolutionary interpretation is that an originally noncoding region entered the genome of the ancestral AlphaPVs. This genetic novelty allowed to explore novel transcription potential, triggering an adaptive radiation that yielded three main viral lineages encoding for different E5 proteins, displaying distinct infection phenotypes. Overall, our results provide an evolutionary scenario for the de novo emergence of viral genes and illustrate the impact of such genotypic novelty in the phenotypic diversity of the viral infections.

Origin and evolution of papillomavirus (onco)genes and genomes.

Papillomaviruses (PVs) are ancient viruses infecting vertebrates, from fishes to mammals. Although the genomes of PVs are small and show conserved synteny, PVs display large genotypic diversity and ample variation in the phenotypic presentation of the infection. Most PV genomes contain two small early genes E6 and E7. In a bunch of closely related human papillomaviruses (HPVs), the E6 and E7 proteins provide the viruses with oncogenic potential. The recent discoveries of PVs without E6 and E7 in different fish species place a new root on the PV tree, and suggest that ancestral PVs consisted of the minimal PV backbone E1-E2-L2-L1. Bayesian phylogenetic analyses date the most recent common ancestor of the PV backbone to 424 million years ago (Ma). Common ancestry tests on extant E6 and E7 genes indicate that they share a common ancestor dating back to at least 184 Ma. In AlphaPVs infecting Old World monkeys and apes, the appearance of the E5 oncogene 53-58 Ma concurred with (i) a significant increase in substitution rate, (ii) a basal radiation and (iii) key gain of functions in E6 and E7. This series of events was instrumental to construct the extant phenotype of oncogenic HPVs. Our results assemble the current knowledge on PV diversity and present an ancient evolutionary timeline punctuated by evolutionary innovations in the history of this successful viral family. This article is part of the theme issue 'Silent cancer agents: multi-disciplinary modelling of human DNA oncoviruses'.

Going, going, gone: predicting the fate of genomic insertions in plant RNA viruses.

Horizontal gene transfer is common among viruses, while they also have highly compact genomes and tend to lose artificial genomic insertions rapidly. Understanding the stability of genomic insertions in viral genomes is therefore relevant for explaining and predicting their evolutionary patterns. Here, we revisit a large body of experimental research on a plant RNA virus, tobacco etch potyvirus (TEV), to identify the patterns underlying the stability of a range of homologous and heterologous insertions in the viral genome. We obtained a wide range of estimates for the recombination rate-the rate at which deletions removing the insertion occur-and these appeared to be independent of the type of insertion and its location. Of the factors we considered, recombination rate was the best predictor of insertion stability, although we could not identify the specific sequence characteristics that would help predict insertion instability. We also considered experimentally the possibility that functional insertions lead to higher mutational robustness through increased redundancy. However, our observations suggest that both functional and non-functional increases in genome size decreased the mutational robustness. Our results therefore demonstrate the importance of recombination rates for predicting the long-term stability and evolution of viral RNA genomes and suggest that there are unexpected drawbacks to increases in genome size for mutational robustness.

2b or Not 2b: Experimental Evolution of Functional Exogenous Sequences in a Plant RNA Virus.

Horizontal gene transfer (HGT) is pervasive in viruses and thought to be a key mechanism in their evolution. On the other hand, strong selective constraints against increasing genome size are an impediment for HGT, rapidly purging horizontally transferred sequences and thereby potentially hindering evolutionary innovation. Here, we explore experimentally the evolutionary fate of viruses with simulated HGT events, using the plant RNA virus Tobacco etch virus (TEV), by separately introducing two functional, exogenous sequences to its genome. One of the events simulates the acquisition of a new function though HGT of a conserved AlkB domain, responsible for the repair of alkylation or methylation damage in many organisms. The other event simulates the acquisition of a sequence that duplicates an existing function, through HGT of the 2b RNA silencing suppressor from Cucumber mosaic virus. We then evolved these two viruses, tracked the maintenance of the horizontally transferred sequences over time, and for the final virus populations, sequenced their genome and measured viral fitness. We found that the AlkB domain was rapidly purged from the TEV genome, restoring fitness to wild-type levels. Conversely, the 2b gene was stably maintained and did not have a major impact on viral fitness. Moreover, we found that 2b is functional in TEV, as it provides a replicative advantage when the RNA silencing suppression domain of HC-Pro is mutated. These observations suggest a potentially interesting role for HGT of short functional sequences in ameliorating evolutionary constraints on viruses, through the duplication of functions.

High virulence does not necessarily impede viral adaptation to a new host: a case study using a plant RNA virus.

BACKGROUND: Theory suggests that high virulence could hinder between-host transmission of microparasites, and that virulence therefore will evolve to lower levels. Alternatively, highly virulent microparasites could also curtail host development, thereby limiting both the host resources available to them and their own within-host effective population size. In this case, high virulence might restrain the mutation supply rate and increase the strength with which genetic drift acts on microparasite populations. Thereby, this alternative explanation limits the microparasites' potential to adapt to the host and ultimately the ability to evolve lower virulence. As a first exploration of this hypothesis, we evolved Tobacco etch virus carrying an eGFP fluorescent marker in two semi-permissive host species, Nicotiana benthamiana and Datura stramonium, for which it has a large difference in virulence. We compared the results to those previously obtained in the natural host, Nicotiana tabacum, where we have shown that carriage of eGFP has a high fitness cost and its loss serves as a real-time indicator of adaptation. RESULTS: After over half a year of evolution, we sequenced the genomes of the evolved lineages and measured their fitness. During the evolution experiment, marker loss leading to viable virus variants was only observed in one lineage of the host for which the virus has low virulence, D. stramonium. This result was consistent with the observation that there was a fitness cost of eGFP in this host, while surprisingly no fitness cost was observed in the host for which the virus has high virulence, N. benthamiana. Furthermore, in both hosts we observed increases in viral fitness in few lineages, and host-specific convergent evolution at the genomic level was only found in N. benthamiana. CONCLUSIONS: The results of this study do not lend support to the hypothesis that high virulence impedes microparasites' evolution. Rather, they exemplify that jumps between host species can be game changers for evolutionary dynamics. When considering the evolution of genome architecture, host species jumps might play a very important role, by allowing evolutionary intermediates to be competitive.

Predicting the Stability of Homologous Gene Duplications in a Plant RNA Virus.

One of the striking features of many eukaryotes is the apparent amount of redundancy in coding and non-coding elements of their genomes. Despite the possible evolutionary advantages, there are fewer examples of redundant sequences in viral genomes, particularly those with RNA genomes. The factors constraining the maintenance of redundant sequences in present-day RNA virus genomes are not well known. Here, we use Tobacco etch virus, a plant RNA virus, to investigate the stability of genetically redundant sequences by generating viruses with potentially beneficial gene duplications. Subsequently, we tested the viability of these viruses and performed experimental evolution. We found that all gene duplication events resulted in a loss of viability or in a significant reduction in viral fitness. Moreover, upon analyzing the genomes of the evolved viruses, we always observed the deletion of the duplicated gene copy and maintenance of the ancestral copy. Interestingly, there were clear differences in the deletion dynamics of the duplicated gene associated with the passage duration and the size and position of the duplicated copy. Based on the experimental data, we developed a mathematical model to characterize the stability of genetically redundant sequences, and showed that fitness effects are not enough to predict genomic stability. A context-dependent recombination rate is also required, with the context being the duplicated gene and its position. Our results therefore demonstrate experimentally the deleterious nature of gene duplications in RNA viruses. Beside previously described constraints on genome size, we identified additional factors that reduce the likelihood of the maintenance of duplicated genes.

Multiple Barriers to the Evolution of Alternative Gene Orders in a Positive-Strand RNA Virus.

The order in which genes are organized within a genome is generally not conserved between distantly related species. However, within virus orders and families, strong conservation of gene order is observed. The factors that constrain or promote gene-order diversity are largely unknown, although the regulation of gene expression is one important constraint for viruses. Here we investigate why gene order is conserved for a positive-strand RNA virus encoding a single polyprotein in the context of its authentic multicellular host. Initially, we identified the most plausible trajectory by which alternative gene orders could evolve. Subsequently, we studied the accessibility of key steps along this evolutionary trajectory by constructing two virus intermediates: (1) duplication of a gene followed by (2) loss of the ancestral gene. We identified five barriers to the evolution of alternative gene orders. First, the number of viable positions for reordering is limited. Second, the within-host fitness of viruses with gene duplications is low compared to the wild-type virus. Third, after duplication, the ancestral gene copy is always maintained and never the duplicated one. Fourth, viruses with an alternative gene order have even lower fitness than viruses with gene duplications. Fifth, after more than half a year of evolution in isolation, viruses with an alternative gene order are still vastly inferior to the wild-type virus. Our results show that all steps along plausible evolutionary trajectories to alternative gene orders are highly unlikely. Hence, the inaccessibility of these trajectories probably contributes to the conservation of gene order in present-day viruses.

Genetic variability and evolutionary analysis of parietaria mottle virus: role of selection and genetic exchange.

The genetic variability and evolution of parietaria mottle virus (PMoV) of the genus Ilarvirus was studied by analyzing nucleotide sequences of 2b and CP genes from isolates collected in different countries. Phylogenetic analysis showed that PMoV isolates clustered in different clades: one (clade I) composed of only Italian isolates and three clades (clades II-IV) including the Spanish isolates. The Greek isolate GrT-1 used in this study was in clade IV for the CP phylogenetic tree whereas it formed a separate branch in the 2b phylogenetic tree. The nucleotide sequence diversity of both the 2b and CP genes was low (0.062 ± 0.006 and 0.063 ± 0.006 for 2b and CP, respectively) but higher than those of other ilarviruses. Distribution of synonymous and nonsynonymous substitutions revealed that 2b and CP proteins are under purifying selection, with some positions under diversifying selection. Genetic exchange among Spanish isolates was also detected.

Temporal dynamics of intrahost molecular evolution for a plant RNA virus.

Populations of plant RNA viruses are highly polymorphic in infected plants, which may allow rapid within-host evolution. To understand tobacco etch potyvirus (TEV) evolution, longitudinal samples from experimentally evolved populations in the natural host tobacco and from the alternative host pepper were phenotypically characterized and genetically analyzed. Temporal and compartmental variabilities of TEV populations were quantified using high throughput Illumina sequencing and population genetic approaches. Of the two viral phenotypic traits measured, virulence increased in the novel host but decreased in the original one, and viral load decreased in both hosts, though to a lesser extent in the novel one. Dynamics of population genetic diversity were also markedly different among hosts. Population heterozygosity increased in the ancestral host, with a dominance of synonymous mutations fixed, whereas it did not change or even decreased in the new host, with an excess of nonsynonymous mutations. All together, these observations suggest that directional selection is the dominant evolutionary force in TEV populations evolving in a novel host whereas either diversifying selection or random genetic drift may play a fundamental role in the natural host. To better understand these evolutionary dynamics, we developed a computer simulation model that incorporates the effects of mutation, selection, and drift. Upon parameterization with empirical data from previous studies, model predictions matched the observed patterns, thus reinforcing our idea that the empirical patterns of mutation accumulation represent adaptive evolution.

Genetic variation and possible mechanisms driving the evolution of worldwide fig mosaic virus isolates.

Fig mosaic virus (FMV) is a multipartite negative-sense RNA virus infecting fig trees worldwide. FMV is transmitted by vegetative propagation and grafting of plant materials, and by the eriophyid mite Aceria ficus. In this work, the genetic variation and evolutionary mechanisms shaping FMV populations were characterized. Nucleotide sequences from four genomic regions (each within the genomic RNAs 1, 2, 3, and 4) from FMV isolates from different countries were determined and analyzed. FMV genetic variation was low, as is seen for many other plant viruses. Phylogenetic analysis showed some geographically distant FMV isolates which clustered together, suggesting long-distance migration. The extent of migration was limited, although varied, between countries, such that FMV populations of different countries were genetically differentiated. Analysis using several recombination algorithms suggests that genomes of some FMV isolates originated by reassortment of genomic RNAs from different genetically similar isolates. Comparison between nonsynonymous and synonymous substitutions showed selection acting on some amino acids; however, most evolved neutrally. This and neutrality tests together with the limited gene flow suggest that genetic drift plays an important role in shaping FMV populations.

Relocation of the NIb gene in the tobacco etch potyvirus genome.

Potyviruses express most of their proteins from a long open reading frame that is translated into a large polyprotein processed by three viral proteases. To understand the constraints on potyvirus genome organization, we relocated the viral RNA-dependent RNA polymerase (NIb) cistron to all possible intercistronic positions of the Tobacco etch virus (TEV) polyprotein. Only viruses with NIb at the amino terminus of the polyprotein or in between P1 and HC-Pro were viable in tobacco plants.

Experimental evolution of pseudogenization and gene loss in a plant RNA virus.

Viruses have evolved highly streamlined genomes and a variety of mechanisms to compress them, suggesting that genome size is under strong selection. Horizontal gene transfer has, on the other hand, played an important role in virus evolution. However, evolution cannot integrate initially nonfunctional sequences into the viral genome if they are rapidly purged by selection. Here we report on the experimental evolution of pseudogenization in virus genomes using a plant RNA virus expressing a heterologous gene. When long 9-week passages were performed, the added gene was lost in all lineages, whereas viruses with large genomic deletions were fixed in only two out of ten 3-week lineages and none in 1-week lineages. Illumina next-generation sequencing revealed considerable convergent evolution in the 9- and 3-week lineages with genomic deletions. Genome size was correlated to within-host competitive fitness, although there was no correlation with virus accumulation or virulence. Within-host competitive fitness of the 3-week virus lineages without genomic deletions was higher than for the 1-week lineages. Our results show that the strength of selection for a reduced genome size and the rate of pseudogenization depend on demographic conditions. Moreover, for the 3-week passage condition, we observed increases in within-host fitness, whereas selection was not strong enough to quickly remove the nonfunctional heterologous gene. These results suggest a demographically determined "sweet spot" might exist, where heterologous insertions are not immediately lost while evolution can act to integrate them into the viral genome.