Genetic diversity and cross-species transmissibility of bat-associated picornaviruses from Spain.

BACKGROUND: Emerging zoonotic diseases arise from cross-species transmission events between wild or domesticated animals and humans, with bats being one of the major reservoirs of zoonotic viruses. Viral metagenomics has led to the discovery of many viruses, but efforts have mainly been focused on some areas of the world and on certain viral families. METHODS: We set out to describe full-length genomes of new picorna-like viruses by collecting feces from hundreds of bats captured in different regions of Spain. Viral sequences were obtained by high-throughput Illumina sequencing and analyzed phylogenetically to classify them in the context of known viruses. Linear discriminant analysis (LDA) was performed to infer likely hosts based on genome composition. RESULTS: We found five complete or nearly complete genomes belonging to the family Picornaviridae, including a new species of the subfamily Ensavirinae. LDA suggested that these were true vertebrate viruses, rather than viruses from the bat diet. Some of these viruses were related to picornaviruses previously found in other bat species from distant geographical regions. We also found a calhevirus genome that most likely belongs to a proposed new family within the order Picornavirales, and for which genome composition analysis suggested a plant host. CONCLUSIONS: Our findings describe new picorna-like viral species and variants circulating in the Iberian Peninsula, illustrate the wide geographical distribution and interspecies transmissibility of picornaviruses, and suggest new hosts for calheviruses.

Full-genome sequencing of dozens of new DNA viruses found in Spanish bat feces.

Bats are natural hosts of multiple viruses, many of which have clear zoonotic potential. The search for emerging viruses has been aided by the implementation of metagenomic tools, which have also enabled the detection of unprecedented viral diversity. Currently, this search is mainly focused on RNA viruses, which are largely over-represented in databases. To compensate for this research bias, we analyzed fecal samples from 189 Spanish bats belonging to 22 different species using viral metagenomics. This allowed us to identify 52 complete or near-complete viral genomes belonging to the families Adenoviridae, Circoviridae, Genomoviridae, Papillomaviridae, Parvoviridae, Polyomaviridae and Smacoviridae. Of these, 30 could constitute new species, doubling the number of viruses currently described in Europe. These findings open the door to a more thorough analysis of bat DNA viruses and their zoonotic potential. IMPORTANCE: Metagenomics has become a fundamental tool to characterize the global virosphere, allowing us not only to understand the existing viral diversity and its ecological implications but also to identify new and emerging viruses. RNA viruses have a higher zoonotic potential, but this risk is also present for some DNA virus families. In our study, we analyzed the DNA fraction of fecal samples from 22 Spanish bat species, identifying 52 complete or near-complete genomes of different viral families with zoonotic potential. This doubles the number of genomes currently described in Europe. Metagenomic data often produce partial genomes that can be difficult to analyze. Our work, however, has characterized a large number of complete genomes, thus facilitating their taxonomic classification and enabling different analyses to be carried out to evaluate their zoonotic potential. For example, recombination studies are relevant since this phenomenon could play a major role in cross-species transmission.

SARS-CoV-2 remodels the landscape of small non-coding RNAs with infection time and symptom severity.

The COVID-19 pandemic caused by the coronavirus SARS-CoV-2 has significantly impacted global health, stressing the necessity of basic understanding of the host response to this viral infection. In this study, we investigated how SARS-CoV-2 remodels the landscape of small non-coding RNAs (sncRNA) from a large collection of nasopharyngeal swab samples taken at various time points from patients with distinct symptom severity. High-throughput RNA sequencing analysis revealed a global alteration of the sncRNA landscape, with abundance peaks related to species of 21-23 and 32-33 nucleotides. Host-derived sncRNAs, including microRNAs (miRNAs), transfer RNA-derived small RNAs (tsRNAs), and small nucleolar RNA-derived small RNAs (sdRNAs) exhibited significant differential expression in infected patients compared to controls. Importantly, miRNA expression was predominantly down-regulated in response to SARS-CoV-2 infection, especially in patients with severe symptoms. Furthermore, we identified specific tsRNAs derived from Glu- and Gly-tRNAs as major altered elements upon infection, with 5' tRNA halves being the most abundant species and suggesting their potential as biomarkers for viral presence and disease severity prediction. Additionally, down-regulation of C/D-box sdRNAs and altered expression of tinyRNAs (tyRNAs) were observed in infected patients. These findings provide valuable insights into the host sncRNA response to SARS-CoV-2 infection and may contribute to the development of further diagnostic and therapeutic strategies in the clinic.

Optimized Recovery of Viral DNA and RNA from Blood Plasma for Viral Metagenomics.

Metagenomics is vastly improving our ability to discover new viruses, as well as their possible associations with disease. However, metagenomics has also changed our understanding of viruses in general. This is because we can find viruses in healthy hosts in the absence of disease, which changes the perspective of viruses as mere pathogens and offers a new perspective in which viruses function as important components of ecosystems. In concrete, human blood metagenomics has revealed the presence of different types of viruses in apparently healthy subjects. These viruses are human anelloviruses and, to a lower extent, human pegiviruses. Viral metagenomics' major challenge is the correct isolation of the viral nucleic acids from a specific sample. For the protocol to be successful, all steps must be carefully chosen, in particular those that optimize the recovery of viral nucleic acids. Here, we present a procedure that allows the recovery of both DNA and RNA viruses from plasma samples.

Human Anelloviruses: Influence of Demographic Factors, Recombination, and Worldwide Diversity.

Anelloviruses represent the major and most diverse component of the healthy human virome, referred to as the anellome. In this study, we determined the anellome of 50 blood donors, forming two sex- and age-matched groups. Anelloviruses were detected in 86% of the donors. The number of detected anelloviruses increased with age and was approximately twice as high in men as in women. A total of 349 complete or nearly complete genomes were classified as belonging to torque teno virus (TTV), torque teno mini virus (TTMV), and torque teno midi virus (TTMDV) anellovirus genera (197, 88, and 64 sequences, respectively). Most donors had intergenus (69.8%) or intragenus (72.1%) coinfections. Despite the limited number of sequences, intradonor recombination analysis showed 6 intragenus recombination events in ORF1. As thousands of anellovirus sequences have been described recently, we finally analyzed the global diversity of human anelloviruses. Species richness and diversity were close to saturation in each anellovirus genus. Recombination was found to be the main factor promoting diversity, although its effect was significantly lower in TTV than in TTMV and TTMDV. Overall, our results suggest that differences in diversity between genera may be caused by variations in the relative contribution of recombination. IMPORTANCE Anelloviruses are the most common human infectious viruses and are considered essentially harmless. Compared to other human viruses, they are characterized by enormous diversity, and recombination is suggested to play an important role in their diversification and evolution. Here, by analyzing the composition of the plasma anellome of 50 blood donors, we find that recombination is also a determinant of viral evolution at the intradonor level. On a larger scale, analysis of anellovirus sequences currently available in databases shows that their diversity is close to saturation and differs among the three human anellovirus genera and that recombination is the main factor explaining this intergenus variability. Global characterization of anellovirus diversity could provide clues about possible associations between certain virus variants and pathologies, as well as facilitate the implementation of unbiased PCR-based detection protocols, which may be relevant for using anelloviruses as endogenous markers of immune status.

Exploring the Diversity of the Human Blood Virome.

Metagenomics is greatly improving our ability to discover new viruses, as well as their possible associations with disease. However, metagenomics has also changed our understanding of viruses in general. The vast expansion of currently known viral diversity has revealed a large fraction of non-pathogenic viruses, and offers a new perspective in which viruses function as important components of many ecosystems. In this vein, studies of the human blood virome are often motivated by the search for new viral diseases, especially those associated with blood transfusions. However, these studies have revealed the common presence of apparently non-pathogenic viruses in blood, particularly human anelloviruses and, to a lower extent, human pegiviruses (HPgV). To shed light on the diversity of the human blood virome, we subjected pooled plasma samples from 587 healthy donors in Spain to a viral enrichment protocol, followed by massive parallel sequencing. This showed that anelloviruses were clearly the major component of the blood virome and showed remarkable diversity. In total, we assembled 332 complete or near-complete anellovirus genomes, 50 of which could be considered new species. HPgV was much less frequent, but we, nevertheless, recovered 17 different isolates that we subsequently used for characterizing the diversity of this virus. In-depth investigation of the human blood virome should help to elucidate the ecology of these viruses, and to unveil potentially associated diseases.

Generation of Photopolymerized Microparticles Based on PEGDA Using Microfluidic Devices. Part 1. Initial Gelation Time and Mechanical Properties of the Material.

Photopolymerized microparticles are made of biocompatible hydrogels like Polyethylene Glycol Diacrylate (PEGDA) by using microfluidic devices are a good option for encapsulation, transport and retention of biological or toxic agents. Due to the different applications of these microparticles, it is important to investigate the formulation and the mechanical properties of the material of which they are made of. Therefore, in the present study, mechanical tests were carried out to determine the swelling, drying, soluble fraction, compression, cross-linking density (Mc) and mesh size (ξ) properties of different hydrogel formulations. Tests provided sufficient data to select the best formulation for the future generation of microparticles using microfluidic devices. The initial gelation times of the hydrogels formulations were estimated for their use in the photopolymerization process inside a microfluidic device. Obtained results showed a close relationship between the amount of PEGDA used in the hydrogel and its mechanical properties as well as its initial gelation time. Consequently, it is of considerable importance to know the mechanical properties of the hydrogels made in this research for their proper manipulation and application. On the other hand, the initial gelation time is crucial in photopolymerizable hydrogels and their use in continuous systems such as microfluidic devices.

Deep viral blood metagenomics reveals extensive anellovirus diversity in healthy humans.

Human blood metagenomics has revealed the presence of different types of viruses in apparently healthy subjects. By far, anelloviruses constitute the viral family that is more frequently found in human blood, although amplification biases and contaminations pose a major challenge in this field. To investigate this further, we subjected pooled plasma samples from 120 healthy donors in Spain to high-speed centrifugation, RNA and DNA extraction, random amplification, and massive parallel sequencing. Our results confirm the extensive presence of anelloviruses in such samples, which represented nearly 97% of the total viral sequence reads obtained. We assembled 114 different viral genomes belonging to this family, revealing remarkable diversity. Phylogenetic analysis of ORF1 suggested 28 potentially novel anellovirus species, 24 of which were validated by Sanger sequencing to discard artifacts. These findings underscore the importance of implementing more efficient purification procedures that enrich the viral fraction as an essential step in virome studies and question the suggested pathological role of anelloviruses.

Genome-scale analysis of evolutionary rate and selection in a fast-expanding Spanish cluster of HIV-1 subtype F1.

This work is aimed at assessing the presence of positive selection and/or shifts of the evolutionary rate in a fast-expanding HIV-1 subtype F1 transmission cluster affecting men who have sex with men in Spain. We applied Bayesian coalescent phylogenetics and selection analyses to 23 full-coding region sequences from patients belonging to that cluster, along with other 19 F1 epidemiologically-unrelated sequences. A shift in the overall evolutionary rate of the virus, explained by positively selected sites in the cluster, was detected. We also found one substitution in Nef (H89F) that was specific to the cluster and experienced positive selection. These results suggest that fast transmission could have been facilitated by some inherent genetic properties of this HIV-1 variant.

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.

[Genetic analysis of a family with Von Hippel-Lindau syndrome]. / Análisis genético en una familia con síndrome de Von Hippel-Lindau.

Von Hippel-Lindau syndrome (VHL) is an autosomal dominant inherited disease associated with mutations in the VHL tumour suppressor gene located on chromosome 3p25. VHL is characterized by the development of multiple malignant and benign tumours in the central nervous system and internal organs, including liver, pancreas and the adrenal gland. More than 823 different mutations of the VHL gene have currently been identified. In the present study we describe the case of a family affected by VHL treated at the University Hospital of La Ribera and the results of the genetic analysis of three relatives, identifying the mutation R167G in exon 3 of VHL gene as the cause of VHL syndrome in this family.

A putative antiviral role of plant cytidine deaminases.

BACKGROUND: A mechanism of innate antiviral immunity operating against viruses infecting mammalian cells has been described during the last decade.  Host cytidine deaminases ( e.g., APOBEC3 proteins) edit viral genomes, giving rise to hypermutated nonfunctional viruses; consequently, viral fitness is reduced through lethal mutagenesis.  By contrast, sub-lethal hypermutagenesis may contribute to virus evolvability by increasing population diversity.  To prevent genome editing, some viruses have evolved proteins that mediate APOBEC3 degradation.  The model plant Arabidopsis thaliana genome encodes nine cytidine deaminases ( AtCDAs), raising the question of whether deamination is an antiviral mechanism in plants as well. METHODS: Here we tested the effects of expression of AtCDAs on the pararetrovirus Cauliflower mosaic virus (CaMV). Two different experiments were carried out. First, we transiently overexpressed each one of the nine A. thalianaAtCDA genes in Nicotianabigelovii plants infected with CaMV, and characterized the resulting mutational spectra, comparing them with those generated under normal conditions.  Secondly, we created A. thaliana transgenic plants expressing an artificial microRNA designed to knock-out the expression of up to six AtCDA genes.  This and control plants were then infected with CaMV.  Virus accumulation and mutational spectra where characterized in both types of plants. RESULTS:  We have shown that the A. thalianaAtCDA1 gene product exerts a mutagenic activity, significantly increasing the number of G to A mutations in vivo, with a concomitant reduction in the amount of CaMV genomes accumulated.  Furthermore, the magnitude of this mutagenic effect on CaMV accumulation is positively correlated with the level of AtCDA1 mRNA expression in the plant. CONCLUSIONS: Our results suggest that deamination of viral genomes may also work as an antiviral mechanism in plants.

Multi-virion infectious units arise from free viral particles in an enveloped virus.

Many animal viruses are enveloped in a lipid bilayer taken up from cellular membranes. Because viral surface proteins bind to these membranes to initiate infection, we hypothesized that free virions may also be capable of interacting with the envelopes of other virions extracellularly. Here, we demonstrate this hypothesis in the vesicular stomatitis virus (VSV), a prototypic negative-strand RNA virus composed of an internal ribonucleocapsid, a matrix protein and an external envelope1. Using microscopy, dynamic light scattering, differential centrifugation and flow cytometry, we show that free viral particles can spontaneously aggregate into multi-virion infectious units. We also show that, following establishment of these contacts, different viral genetic variants are co-transmitted to the same target cell. Furthermore, virion-virion binding can determine key aspects of viral fitness such as antibody escape. In purified virions, this process is driven by protein-lipid interactions probably involving the VSV surface glycoprotein and phosphatidylserine. Whereas we found that multi-virion complexes occurred unfrequently in standard cell cultures, they were abundant in other fluids such as saliva, a natural VSV shedding route2. Our findings contrast with the commonly accepted perception of virions as passive propagules and show the ability of enveloped viruses to establish collective infectious units, which could in turn facilitate the evolution of virus-virus interactions and of social-like traits3.

Análisis genético en una familia con síndrome de Von Hippel-Lindau / Genetic analysis of a family with Von Hippel-Lindau syndrome

El síndrome de Von Hippel-Lindau (VHL) es una enfermedad hereditaria autosómica dominante asociada a mutaciones en el gen supresor de tumores VHL localizado en el cromosoma 3p25. El VHL se caracteriza por el desarrollo de múltiples tumores malignos y benignos en el sistema nervioso central y órganos internos, incluyendo el hígado, el páncreas y la glándula adrenal. Actualmente se conocen más de 823mutaciones distintas del gen VHL. En el presente estudio describimos el caso clínico de una familia afecta por VHL tratada en el Hospital Universitario de la Ribera y los resultados del análisis genético de 3 de sus miembros, permitiendo identificar como causa de la enfermedad a la mutación R167G en el exón 3 del gen VHL (AU)

Von Hippel-Lindau syndrome (VHL) is an autosomal dominant inherited disease associated with mutations in the VHL tumour suppressor gene located on chromosome 3p25. VHL is characterized by the development of multiple malignant and benign tumours in the central nervous system and internal organs, including liver, pancreas and the adrenal gland. More than 823 different mutations of the VHL gene have currently been identified. In the present study we describe the case of a family affected by VHL treated at the University Hospital of La Ribera and the results of the genetic analysis of three relatives, identifying the mutation R167G in exon 3 of VHL gene as the cause of VHL syndrome in this family (AU)

Prevalencia mutacional de los genes BRCA1 y BRCA2 en pacientes con cáncer de mama y ovario atendidos en el Área de Salud de La Ribera (Comunidad Valenciana, España) / Mutational prevalence of BRCA1 and BRCA2 genes in patients with breast and ovarian cancer treated in the Health Area of La Ribera (Comunidad Valenciana, Spain)

El cáncer de mama es la neoplasia más frecuente en la mujer, estando implicadas mutaciones en la línea germinal de los genes BRCA1 y BRCA2 en un 5-10% de los casos. El análisis genético precoz en familias portadoras de mutaciones en BRCA permite la detección de pacientes de riesgo en los que llevar a cabo unas adecuadas medidas de seguimiento o profilácticas. Se han descrito más de 3.000 mutaciones distintas en BRCA relacionadas con la enfermedad, variando su prevalencia alélica según el área geográfica analizada. El objetivo de este trabajo ha sido determinar mediante secuenciación masiva (NGS) qué mutaciones patológicas en los genes BRCA1 y BRCA2 son las prevalentes en los pacientes atendidos en la Unidad de Riesgo de Cáncer de Mama Familiar del Hospital de La Ribera y cuáles son sus características clinicopatológicas. Nuestro estudio muestra que los pacientes con mutaciones en BRCA fueron principalmente mujeres premenopáusicas con algún tipo de cáncer y con antecedentes familiares. La mutación con mayor prevalencia detectada en nuestra población de estudio fue c.9026_9030delATCAT del exón 23 en el gen BRCA2 (AU)

Breast cancer is the most frequent neoplasia in women, with BRCA1 and BRCA2 germline mutations being involved in 5%-10% of cases. Early genetic analyses in families affected by this disease are crucial for the identification of family members at risk and thus candidates for surveillance programmes. More than 3,000 different mutations in BRCA gene have been described with different allelic prevalence depending on the geographic area analyzed. The aim of this study was to determine, by massive sequencing (NGS), which pathological mutations in BRCA1 and BRCA2 genes are prevalent in patients treated in the Health Area of La Ribera (Comunidad Valenciana, Spain) and its relation with clinic-pathological characteristics. Our data show that patients with BRCA mutations are mainly premenopausal women with cancer and a family history. The predominant detected mutation in our region is c.9026_9030delATCAT in exon 23 of BRCA2 gene (AU)

Genome-Wide Estimation of the Spontaneous Mutation Rate of Human Adenovirus 5 by High-Fidelity Deep Sequencing.

Rates of spontaneous mutation determine the ability of viruses to evolve, infect new hosts, evade immunity and undergo drug resistance. Contrarily to RNA viruses, few mutation rate estimates have been obtained for DNA viruses, because their high replication fidelity implies that new mutations typically fall below the detection limits of Sanger and standard next-generation sequencing. Here, we have used a recently developed high-fidelity deep sequencing technique (Duplex Sequencing) to score spontaneous mutations in human adenovirus 5 under conditions of minimal selection. Based on >200 single-base spontaneous mutations detected throughout the entire viral genome, we infer an average mutation rate of 1.3 × 10-7 per base per cell infection cycle. This value is similar to those of other, large double-stranded DNA viruses, but an order of magnitude lower than those of single-stranded DNA viruses, consistent with the possible action of post-replicative repair. Although the mutation rate did not vary strongly along the adenovirus genome, we found several sources of mutation rate heterogeneity. First, two regions mapping to transcription units L3 and E1B-IVa2 were significantly depleted for mutations. Second, several point insertions/deletions located within low-complexity sequence contexts appeared recurrently, suggesting mutational hotspots. Third, mutation probability increased at GpC dinucleotides. Our findings suggest that host factors may influence the distribution of spontaneous mutations in human adenoviruses and potentially other nuclear DNA viruses.

Highly heterogeneous mutation rates in the hepatitis C virus genome.

Spontaneous mutations are the ultimate source of genetic variation and have a prominent role in evolution. RNA viruses such as hepatitis C virus (HCV) have extremely high mutation rates, but these rates have been inferred from a minute fraction of genome sites, limiting our view of how RNA viruses create diversity. Here, by applying high-fidelity ultradeep sequencing to a modified replicon system, we scored >15,000 spontaneous mutations, encompassing more than 90% of the HCV genome. This revealed >1,000-fold differences in mutability across genome sites, with extreme variations even between adjacent nucleotides. We identify base composition, the presence of high- and low-mutation clusters and transition/transversion biases as the main factors driving this heterogeneity. Furthermore, we find that mutability correlates with the ability of HCV to diversify in patients. These data provide a site-wise baseline for interrogating natural selection, genetic load and evolvability in HCV, as well as for evaluating drug resistance and immune evasion risks.

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.

Human norovirus hyper-mutation revealed by ultra-deep sequencing.

Human noroviruses (NoVs) are a major cause of gastroenteritis worldwide. It is thought that, similar to other RNA viruses, high mutation rates allow NoVs to evolve fast and to undergo rapid immune escape at the population level. However, the rate and spectrum of spontaneous mutations of human NoVs have not been quantified previously. Here, we analyzed the intra-patient diversity of the NoV capsid by carrying out RT-PCR and ultra-deep sequencing with 100,000-fold coverage of 16 stool samples from symptomatic patients. This revealed the presence of low-frequency sequences carrying large numbers of U-to-C or A-to-G base transitions, suggesting a role for hyper-mutation in NoV diversity. To more directly test for hyper-mutation, we performed transfection assays in which the production of mutations was restricted to a single cell infection cycle. This confirmed the presence of sequences with multiple U-to-C/A-to-G transitions, and suggested that hyper-mutation contributed a large fraction of the total NoV spontaneous mutation rate. The type of changes produced and their sequence context are compatible with ADAR-mediated editing of the viral RNA.

Constrained evolvability of interferon suppression in an RNA virus.

Innate immunity responses controlled by interferon (IFN) are believed to constitute a major selective pressure shaping viral evolution. Viruses encode a variety of IFN suppressors, but these are often multifunctional proteins that also play essential roles in other steps of the viral infection cycle, possibly limiting their evolvability. Here, we experimentally evolved a vesicular stomatitis virus (VSV) mutant carrying a defect in the matrix protein (M∆51) that abolishes IFN suppression and that has been previously used in the context of oncolytic virotherapy. Serial transfers of this virus in normal, IFN-secreting cells led to a modest recovery of IFN blocking capacity and to weak increases in viral fitness. Full-genome ultra-deep sequencing and phenotypic analysis of population variants revealed that the anti-IFN function of the matrix protein was not restored, and that the Mdelta51 defect was instead compensated by changes in the viral phosphoprotein. We also show that adaptation to IFN-secreting cells can be driven by the selection of fast-growing viruses with no IFN suppression capacity, and that these population variants can be trans-complemented by other, IFN-suppressing variants. Our results thus suggest that virus-virus interactions and alternative strategies of innate immunity evasion can determine the evolution of IFN suppression in a virus.