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.

Extremely High Mutation Rate of HIV-1 In Vivo.

Rates of spontaneous mutation critically determine the genetic diversity and evolution of RNA viruses. Although these rates have been characterized in vitro and in cell culture models, they have seldom been determined in vivo for human viruses. Here, we use the intrapatient frequency of premature stop codons to quantify the HIV-1 genome-wide rate of spontaneous mutation in DNA sequences from peripheral blood mononuclear cells. This reveals an extremely high mutation rate of (4.1 ± 1.7) × 10-3 per base per cell, the highest reported for any biological entity. Sequencing of plasma-derived sequences yielded a mutation frequency 44 times lower, indicating that a large fraction of viral genomes are lethally mutated and fail to reach plasma. We show that the HIV-1 reverse transcriptase contributes only 2% of mutations, whereas 98% result from editing by host cytidine deaminases of the A3 family. Hypermutated viral sequences are less abundant in patients showing rapid disease progression compared to normal progressors, highlighting the antiviral role of A3 proteins. However, the amount of A3-mediated editing varies broadly, and we find that low-edited sequences are more abundant among rapid progressors, suggesting that suboptimal A3 activity might enhance HIV-1 genetic diversity and pathogenesis.

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.

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.

Effect of host species on the distribution of mutational fitness effects for an RNA virus.

Knowledge about the distribution of mutational fitness effects (DMFE) is essential for many evolutionary models. In recent years, the properties of the DMFE have been carefully described for some microorganisms. In most cases, however, this information has been obtained only for a single environment, and very few studies have explored the effect that environmental variation may have on the DMFE. Environmental effects are particularly relevant for the evolution of multi-host parasites and thus for the emergence of new pathogens. Here we characterize the DMFE for a collection of twenty single-nucleotide substitution mutants of Tobacco etch potyvirus (TEV) across a set of eight host environments. Five of these host species were naturally infected by TEV, all belonging to family Solanaceae, whereas the other three were partially susceptible hosts belonging to three other plant families. First, we found a significant virus genotype-by-host species interaction, which was sustained by differences in genetic variance for fitness and the pleiotropic effect of mutations among hosts. Second, we found that the DMFEs were markedly different between Solanaceae and non-Solanaceae hosts. Exposure of TEV genotypes to non-Solanaceae hosts led to a large reduction of mean viral fitness, while the variance remained constant and skewness increased towards the right tail. Within the Solanaceae hosts, the distribution contained an excess of deleterious mutations, whereas for the non-Solanaceae the fraction of beneficial mutations was significantly larger. All together, this result suggests that TEV may easily broaden its host range and improve fitness in new hosts, and that knowledge about the DMFE in the natural host does not allow for making predictions about its properties in an alternative host.

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.

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.

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.

The fitness effects of synonymous mutations in DNA and RNA viruses.

Despite being silent with respect to protein sequence, synonymous nucleotide substitutions can be targeted by natural selection directly at the DNA or RNA level. However, there has been no systematic assessment of how frequent this type of selection is. Here, we have constructed 53 single random synonymous substitution mutants of the bacteriophages Qß and ΦX174 by site-directed mutagenesis and assayed their fitness. Analysis of this mutant collection and of previous studies undertaken with a variety of single-stranded (ss) viruses demonstrates that selection at synonymous sites is stronger in RNA viruses than in DNA viruses. We estimate that this type of selection contributes approximately 18% of the overall mutational fitness effects in ssRNA viruses under our assay conditions and that random synonymous substitutions have a 5% chance of being lethal to the virus, whereas in ssDNA viruses, these figures drop to 1.4% and 0%, respectively. In contrast, the effects of nonsynonymous substitutions appear to be similar in ssRNA and ssDNA viruses.

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.

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.

Molecular evolution and phylogeography of potato virus Y based on the CP gene.

Potato virus Y (PVY) is an important plant pathogen with a wide host range that includes, among others, potato, tobacco, tomato and pepper. The coat protein (CP) of PVY has been commonly used in phylogenetic studies for strain classification. In this study, we used a pool of 292 CP sequences from isolates collected worldwide. After detecting and removing recombinant sequences, we applied Bayesian techniques to study the influence of geography and host species in CP population structure and dynamics. Finally, we performed selection and covariation analyses to identify specific amino acids involved in adaptation. Our results show that PVY CP diversification is significantly accounted for by both geographical and host-driven adaptations. Amino acid positions detected as positively selected concentrate in the N-terminal region of the protein. Some of these selected positions may discriminate among strains, and to a much lesser extent, between potato and non-potato isolates.

The fitness effects of random mutations in single-stranded DNA and RNA bacteriophages.

Mutational fitness effects can be measured with relatively high accuracy in viruses due to their small genome size, which facilitates full-length sequencing and genetic manipulation. Previous work has shown that animal and plant RNA viruses are very sensitive to mutation. Here, we characterize mutational fitness effects in single-stranded (ss) DNA and ssRNA bacterial viruses. First, we performed a mutation-accumulation experiment in which we subjected three ssDNA (PhiX174, G4, F1) and three ssRNA phages (Qbeta, MS2, and SP) to plaque-to-plaque transfers and chemical mutagenesis. Genome sequencing and growth assays indicated that the average fitness effect of the accumulated mutations was similar in the two groups. Second, we used site-directed mutagenesis to obtain 45 clones of PhiX174 and 42 clones of Qbeta carrying random single-nucleotide substitutions and assayed them for fitness. In PhiX174, 20% of such mutations were lethal, whereas viable ones reduced fitness by 13% on average. In Qbeta, these figures were 29% and 10%, respectively. It seems therefore that high mutational sensitivity is a general property of viruses with small genomes, including those infecting animals, plants, and bacteria. Mutational fitness effects are important for understanding processes of fitness decline, but also of neutral evolution and adaptation. As such, these findings can contribute to explain the evolution of ssDNA and ssRNA viruses.

The evolutionary genetics of emerging plant RNA viruses.

Over the years, agriculture across the world has been compromised by a succession of devastating epidemics caused by new viruses that spilled over from reservoir species or by new variants of classic viruses that acquired new virulence factors or changed their epidemiological patterns. Viral emergence is usually associated with ecological change or with agronomical practices bringing together reservoirs and crop species. The complete picture is, however, much more complex, and results from an evolutionary process in which the main players are ecological factors, viruses' genetic plasticity, and host factors required for virus replication, all mixed with a good measure of stochasticity. The present review puts emergence of plant RNA viruses into the framework of evolutionary genetics, stressing that viral emergence begins with a stochastic process that involves the transmission of a preexisting viral strain into a new host species, followed by adaptation to the new host.

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.

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.

Effect of ribavirin on the mutation rate and spectrum of hepatitis C virus in vivo.

Their extremely error-prone replication makes RNA viruses targets for lethal mutagenesis. In the case of hepatitis C virus (HCV), the standard treatment includes ribavirin, a base analog with an in vitro mutagenic effect, but the in vivo mode of action of ribavirin remains poorly understood. Here, we test the mutagenic effects of ribavirin plus interferon treatment in vivo using a new method to estimate mutation rates based on the analysis of nonsense mutations. We apply this methodology to a large HCV sequence database containing over 15,000 reverse transcription-PCR molecular clone sequences from 74 patients infected with HCV. We obtained an estimate of the spontaneous mutation rate of ca. 10(-4) substitutions per site or lower, a value within the typically accepted range for RNA viruses. A roughly threefold increase in mutation rate and a significant shift in mutation spectrum were observed in samples from patients undergoing 6 months of interferon plus ribavirin treatment. This result is consistent with the known in vitro mutagenic effect of ribavirin and suggests that the antiviral effect of ribavirin plus interferon treatment is at least partly exerted through lethal mutagenesis.

Selection for robustness in mutagenized RNA viruses.

Mutational robustness is defined as the constancy of a phenotype in the face of deleterious mutations. Whether robustness can be directly favored by natural selection remains controversial. Theory and in silico experiments predict that, at high mutation rates, slow-replicating genotypes can potentially outcompete faster counterparts if they benefit from a higher robustness. Here, we experimentally validate this hypothesis, dubbed the "survival of the flattest," using two populations of the vesicular stomatitis RNA virus. Characterization of fitness distributions and genetic variability indicated that one population showed a higher replication rate, whereas the other was more robust to mutation. The faster replicator outgrew its robust counterpart in standard competition assays, but the outcome was reversed in the presence of chemical mutagens. These results show that selection can directly favor mutational robustness and reveal a novel viral resistance mechanism against treatment by lethal mutagenesis.