Accumulation Dynamics of Defective Genomes during Experimental Evolution of Two Betacoronaviruses.

Virus-encoded replicases often generate aberrant RNA genomes, known as defective viral genomes (DVGs). When co-infected with a helper virus providing necessary proteins, DVGs can multiply and spread. While DVGs depend on the helper virus for propagation, they can in some cases disrupt infectious virus replication, impact immune responses, and affect viral persistence or evolution. Understanding the dynamics of DVGs alongside standard viral genomes during infection remains unclear. To address this, we conducted a long-term experimental evolution of two betacoronaviruses, the human coronavirus OC43 (HCoV-OC43) and the murine hepatitis virus (MHV), in cell culture at both high and low multiplicities of infection (MOI). We then performed RNA-seq at regular time intervals, reconstructed DVGs, and analyzed their accumulation dynamics. Our findings indicate that DVGs evolved to exhibit greater diversity and abundance, with deletions and insertions being the most common types. Notably, some high MOI deletions showed very limited temporary existence, while others became prevalent over time. We observed differences in DVG abundance between high and low MOI conditions in HCoV-OC43 samples. The size distribution of HCoV-OC43 genomes with deletions differed between high and low MOI passages. In low MOI lineages, short and long DVGs were the most common, with an additional cluster in high MOI lineages which became more prevalent along evolutionary time. MHV also showed variations in DVG size distribution at different MOI conditions, though they were less pronounced compared to HCoV-OC43, suggesting a more random distribution of DVG sizes. We identified hotspot regions for deletions that evolved at a high MOI, primarily within cistrons encoding structural and accessory proteins. In conclusion, our study illustrates the widespread formation of DVGs during betacoronavirus evolution, influenced by MOI and cell- and virus-specific factors.

Decay of HCoV-OC43 infectivity is lower in cell debris-containing media than in fresh culture media.

In the quantitative description of viral dynamics within cell cultures and, more broadly, in modeling within-host viral infections, a question that commonly arises is whether the degradation of a fraction of the virus could be disregarded in comparison with the massive synthesis of new viral particles. Surprisingly, quantitative data on the synthesis and degradation rates of RNA viruses in cell cultures are scarce. In this study, we investigated the decay of the human betacoronavirus OC43 (HCoV-OC43) infectivity in cell culture lysates and in fresh media. Our findings revealed a significantly slower viral decay rate in the medium containing lysate cells compared to the fresh medium. This observation suggests that the presence of cellular debris from lysed cells may offer protection or stabilize virions, slowing down their degradation. Moreover, the growth rate of HCoV-OC43 infectivity is significantly higher than degradation as long as there are productive cells in the medium, suggesting that, as a first approximation, degradation can be neglected during early infection.

A binary interaction map between turnip mosaic virus and Arabidopsis thaliana proteomes.

Viruses are obligate intracellular parasites that have co-evolved with their hosts to establish an intricate network of protein-protein interactions. Here, we followed a high-throughput yeast two-hybrid screening to identify 378 novel protein-protein interactions between turnip mosaic virus (TuMV) and its natural host Arabidopsis thaliana. We identified the RNA-dependent RNA polymerase NIb as the viral protein with the largest number of contacts, including key salicylic acid-dependent transcription regulators. We verified a subset of 25 interactions in planta by bimolecular fluorescence complementation assays. We then constructed and analyzed a network comprising 399 TuMV-A. thaliana interactions together with intravirus and intrahost connections. In particular, we found that the host proteins targeted by TuMV are enriched in different aspects of plant responses to infections, are more connected and have an increased capacity to spread information throughout the cell proteome, display higher expression levels, and have been subject to stronger purifying selection than expected by chance. The proviral or antiviral role of ten host proteins was validated by characterizing the infection dynamics in the corresponding mutant plants, supporting a proviral role for the transcriptional regulator TGA1. Comparison with similar studies with animal viruses, highlights shared fundamental features in their mode of action.

High-throughput sequencing (HTS) for the analysis of viral populations.

The development of High-Throughput Sequencing (HTS) technologies is having a major impact on the genomic analysis of viral populations. Current HTS platforms can capture nucleic acid variation across millions of genes for both selected amplicons and full viral genomes. HTS has already facilitated the discovery of new viruses, hinted new taxonomic classifications and provided a deeper and broader understanding of their diversity, population and genetic structure. Hence, HTS has already replaced standard Sanger sequencing in basic and applied research fields, but the next step is its implementation as a routine technology for the analysis of viruses in clinical settings. The most likely application of this implementation will be the analysis of viral genomics, because the huge population sizes, high mutation rates and very fast replacement of viral populations have demonstrated the limited information obtained with Sanger technology. In this review, we describe new technologies and provide guidelines for the high-throughput sequencing and genetic and evolutionary analyses of viral populations and metaviromes, including software applications. With the development of new HTS technologies, new and refurbished molecular and bioinformatic tools are also constantly being developed to process and integrate HTS data. These allow assembling viral genomes and inferring viral population diversity and dynamics. Finally, we also present several applications of these approaches to the analysis of viral clinical samples including transmission clusters and outbreak characterization.

Role of APOBEC3H in the Viral Control of HIV Elite Controller Patients.

Background APOBEC3H (A3H) gene presents variation at 2 positions (rs139297 and rs79323350) leading to a non-functional protein. So far, there is no information on the role played by A3H in spontaneous control of HIV. The aim of this study was to evaluate the A3H polymorphisms distribution in a well-characterized group of Elite Controller (EC) subjects. Methods We analyzed the genotype distribution of two different SNPs (rs139297 and rs79323350) of A3H in 30 EC patients and compared with 11 non-controller (NC) HIV patients. Genotyping was performed by PCR, cloning and Sanger sequencing. Both polymorphisms were analyzed jointly in order to adequately attribute the active or inactive status of A3H protein. Results EC subjects included in this study were able to maintain a long-term sustained spontaneous HIV-viral control and optimal CD4-T-cell counts; however, haplotypes leading to an active protein were very poorly represented in these patients. We found that the majority of EC subjects (23/30; 77%) presented allelic combinations leading to an inactive A3H protein, a frequency slightly lower than that observed for NC studied patients (10/11; 91%). Conclusions The high prevalence of non-functional protein coding-genotypes in EC subjects seems to indicate that other innate restriction factors different from APOBEC3H could be implicated in the replication control exhibited by these subjects.

The transcriptomics of an experimentally evolved plant-virus interaction.

Models of plant-virus interaction assume that the ability of a virus to infect a host genotype depends on the matching between virulence and resistance genes. Recently, we evolved tobacco etch potyvirus (TEV) lineages on different ecotypes of Arabidopsis thaliana, and found that some ecotypes selected for specialist viruses whereas others selected for generalists. Here we sought to evaluate the transcriptomic basis of such relationships. We have characterized the transcriptomic responses of five ecotypes infected with the ancestral and evolved viruses. Genes and functional categories differentially expressed by plants infected with local TEV isolates were identified, showing heterogeneous responses among ecotypes, although significant parallelism existed among lineages evolved in the same ecotype. Although genes involved in immune responses were altered upon infection, other functional groups were also pervasively over-represented, suggesting that plant resistance genes were not the only drivers of viral adaptation. Finally, the transcriptomic consequences of infection with the generalist and specialist lineages were compared. Whilst the generalist induced very similar perturbations in the transcriptomes of the different ecotypes, the perturbations induced by the specialist were divergent. Plant defense mechanisms were activated when the infecting virus was specialist but they were down-regulated when infecting with generalist.

Evaluating the within-host fitness effects of mutations fixed during virus adaptation to different ecotypes of a new host.

The existence of genetic variation for resistance in host populations is assumed to be essential to the spread of an emerging virus. Models predict that the rate of spread slows down with the increasing frequency and higher diversity of resistance alleles in the host population. We have been using the experimental pathosystem Arabidopsis thaliana-tobacco etch potyvirus (TEV) to explore the interplay between genetic variation in host's susceptibility and virus diversity. We have recently shown that TEV populations evolving in A. thaliana ecotypes that differ in susceptibility to infection gained within-host fitness, virulence and infectivity in a manner compatible with a gene-for-gene model of host-parasite interactions: hard-to-infect ecotypes were infected by generalist viruses, whereas easy-to-infect ecotypes were infected by every virus. We characterized the genomes of the evolved viruses and found cases of host-driven convergent mutations. To gain further insights in the mechanistic basis of this gene-for-gene model, we have generated all viral mutations individually as well as in specific combinations and tested their within-host fitness effects across ecotypes. Most of these mutations were deleterious or neutral in their local ecotype and only a very reduced number had a host-specific beneficial effect. We conclude that most of the mutations fixed during the evolution experiment were so by drift or by selective sweeps along with the selected driver mutation. In addition, we evaluated the ruggedness of the underlying adaptive fitness landscape and found that mutational effects were mostly multiplicative, with few cases of significant epistasis.

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.

Emerging viruses: why they are not jacks of all trades?

In order to limit the impact of the recent pandemics ignited by viral host jumps, it is necessary to better understand the ecological and evolutionary factors influencing the early steps of emergence and the interactions between them. Antagonistic pleiotropy, that is, the negative fitness effect in the primary host of mutations allowing the infection of and the multiplication in a new host, has long been thought to be the main limitation to the evolution of generalist viruses and thus to emergence. However, the accumulation of experimental examples contradicting the hypothesis of antagonistic pleiotropy has highlighted the importance of other factors such as the epistasis between mutations increasing the adaptation to a new host. Epistasis is pervasive in viruses, affects the shape of the adaptive landscape and consequently the accessibility of evolutionary pathways. Finally, recent studies have gone steps further in the complexity of viral fitness determinism and stressed the potential importance of the epistatic pleiotropy and of the impact of host living conditions.

Experimental evolution of an emerging plant virus in host genotypes that differ in their susceptibility to infection.

This study evaluates the extent to which genetic differences among host individuals from the same species condition the evolution of a plant RNA virus. We performed a threefold replicated evolution experiment in which Tobacco etch potyvirus isolate At17b (TEV-At17b), adapted to Arabidopsis thaliana ecotype Ler-0, was serially passaged in five genetically heterogeneous ecotypes of A. thaliana. After 15 passages we found that evolved viruses improved their fitness, showed higher infectivity and stronger virulence in their local host ecotypes. The genome of evolved lineages was sequenced and putative adaptive mutations identified. Host-driven convergent mutations have been identified. Evidences supported selection for increased translational efficiency. Next, we sought for the specificity of virus adaptation by infecting all five ecotypes with all 15 evolved virus populations. We found that some ecotypes were more permissive to infection than others, and that some evolved virus isolates were more specialist/generalist than others. The bipartite network linking ecotypes with evolved viruses was significantly nested but not modular, suggesting that hard-to-infect ecotypes were infected by generalist viruses whereas easy-to-infect ecotypes were infected by all viruses, as predicted by a gene-for-gene model of infection.

Intra-specific variability and biological relevance of P3N-PIPO protein length in potyviruses.

BACKGROUND: Pipo was recently described as a new ORF encoded within the genome of the Potyviridae family members (PNAS 105:5897-5902, 2008). It is embedded within the P3 cistron and is translated in the +2 reading frame relative to the potyviral long ORF as the P3N-PIPO fusion protein. In this work, we first collected pipo nucleotide sequences available for different isolates of 48 Potyvirus species. Second, to determine the biological implications of variation in pipo length, we measured infectivity, viral accumulation, cell-to-cell and systemic movements for two Turnip mosaic virus (TuMV) variants with pipo alleles of different length in three different susceptible host species, and tested for differences between the two variants. RESULTS: In addition to inter-specific variation, there was high variation in the length of the PIPO protein among isolates within species (ranging from 1 to 89 amino acids). Furthermore, selection analyses on the P3 cistron did not account for the existence of stop codons in the pipo ORF, but showed that positive selection was significant in the overlapping region for Potato virus Y (PVY) and TuMV. In some cases, variability in length was associated with host species, geographic provenance and/or other strain features. We found significant empirical differences among the phenotypes associated with TuMV pipo alleles, though the magnitude and sign of the effects were host-dependent. CONCLUSIONS: The combination of computational molecular evolution analyses and experiments stemming from these analyses provide clues about the selective pressures acting upon the different-length pipo alleles and show that variation in length may be maintained by host-driven selection.

Characterization of hepatitis C virus intra- and intergenotypic chimeras reveals a role of the glycoproteins in virus envelopment.

Hepatitis C virus (HCV) is highly variable and associated with chronic liver disease. Viral isolates are grouped into seven genotypes (GTs). Accumulating evidence indicates that viral determinants in the core to NS2 proteins modulate the efficiency of virus production. However, the role of the glycoproteins E1 and E2 in this process is currently poorly defined. Therefore, we constructed chimeric viral genomes to explore the role of E1 and E2 in HCV assembly. Comparison of the kinetics and efficiency of particle production by intragenotypic chimeras highlighted core and p7 as crucial determinants for efficient virion release. Glycoprotein sequences, however, had only a minimal impact on this process. In contrast, in the context of intergenotypic HCV chimeras, HCV assembly was profoundly influenced by glycoprotein genes. On the one hand, insertion of GT1a-derived (H77) E1-E2 sequences into a chimeric GT2a virus (Jc1) strongly suppressed virus production. On the other hand, replacement of H77 glycoproteins within the GT1a-GT2a chimeric genome H77/C3 by GT2a-derived (Jc1) E1-E2 increased infectious particle production. Thus, within intergenotypic chimeras, glycoprotein features strongly modulate virus production. Replacement of Jc1 glycoprotein genes by H77-derived E1-E2 did not grossly affect subcellular localization of core, E2, and NS2. However, it caused an accumulation of nonenveloped core protein and increased abundance of nonenveloped core protein structures with slow sedimentation. These findings reveal an important role for the HCV glycoproteins E1 and E2 in membrane envelopment, which likely depends on a genotype-specific interplay with additional viral factors.

Transcript Profiling of Different Arabidopsis thaliana Ecotypes in Response to Tobacco etch potyvirus Infection.

The use of high-throughput transcript profiling techniques has opened the possibility of identifying, in a single experiment, multiple host mRNAs whose levels of accumulation are altered in response to virus infection. Several studies have used this approach to analyze the response of Arabidopsis thaliana to the infection by different RNA and DNA viruses. However, the possible differences in response of genetically heterogeneous ecotypes of the plant to the same virus have never been addressed before. Here we have used a strain of Tobacco etch potyvirus (TEV) experimentally adapted to A. thaliana ecotype Ler-0 and a set of seven plant ecotypes to tackle this question. Each ecotype was inoculated with the same amount of the virus and the outcome of infection characterized phenotypically (i.e., virus infectivity, accumulation, and symptoms development). Using commercial microarrays containing probes for more than 43,000 A. thaliana transcripts, we explored the effect of viral infection on the plant transcriptome. In general, we found that ecotypes differ in the way they perceive and respond to the virus. Some ecotypes developed strong symptoms and accumulated large amounts of viral genomes, while others only developed mild symptoms and accumulated less virus. At the transcriptomic level, ecotypes could be classified into two groups according to the particular genes whose expression was altered upon infection. Moreover, a functional enrichment analyses showed that the two groups differed in the nature of the altered biological processes. For the group constituted by ecotypes developing milder symptoms and allowing for lower virus accumulation, genes involved in abiotic stresses and in the construction of new tissues tend to be up-regulated. For those ecotypes in which infection was more severe and productive, defense genes tend to be up-regulated, deviating the necessary resources from building new tissues.

Luria-delbruck estimation of turnip mosaic virus mutation rate in vivo.

A potential drawback of recent antiviral therapies based on the transgenic expression of artificial microRNAs is the ease with which viruses may generate escape mutations. Using a variation of the classic Luria-Delbrück fluctuation assay, we estimated that the spontaneous mutation rate in the artificial microRNA (amiR) target of a plant virus was ca. 6 × 10(-5) per replication event.

Characterization of the interaction between hepatitis C virus NS5B and the human oestrogen receptor alpha.

The RNA-dependent RNA polymerase (NS5B) of hepatitis C virus (HCV) is part of the viral replicative complex and plays a crucial role in HCV replication. It has been described that NS5B interacts with cellular proteins, and that interactions between NS5B and host proteins are crucial for viral replication. Some of the host factors involved in the HCV replication cycle include the oestrogen receptor alpha (ESR1), protein kinases (c-Src) and chaperones (Hsp70). In this report, we determine the requirements for the interplay between NS5B and the domain C of ESR1 (ESR1C) by using Förster Resonance Energy Transfer. NS5B-ESR1C and ESR1C-ESR1C interactions are dependent on ionic strength, indicating that contacts are mainly electrostatic. Additionally, NS5B residues involved in NS5B oligomerization were also essential for NS5B-ESR1C interaction. The study of the interactions among viral and host factors will provide data to establish innovative therapeutic strategies and the development of new antiviral drugs.