RNA virus populations as quasispecies

Investigating Variability in Viral Presence and Abundance across Soybean Seed Development Stages Using Transcriptome Analysis

Plant transcriptomes offer a valuable resource for studying viral communities (viromes). In this study, we explore how plant transcriptome data can be applied to virome research. We analyzed 40 soybean transcriptomes across different growth stages and identified six viruses: broad bean wilt virus 2 (BBWV2), brassica yellow virus (BrYV), beet western yellow virus (BWYV), cucumber mosaic virus (CMV), milk vetch dwarf virus (MDV), and soybean mosaic virus (SMV). SMV was the predominant virus in both Glycine max (GM) and Glycine soja (GS) cultivars. Our analysis confirmed its abundance in both, while BBWV2 and CMV were more prevalent in GS than GM. The viral proportions varied across developmental stages, peaking in open flowers. Comparing viral abundance measured by viral reads and fragments per kilobase of transcript per million (FPKM) values revealed insights. SMV showed similar FPKM values in GM and GS, but BBWV2 and CMV displayed higher FPKM proportions in GS. Notably, the differences in viral abundance between GM and GS were generally insignificant based on the FPKM values across developmental stages, except for the apical bud stage in four GM cultivars. We also detected MDV, a multi-segmented virus, in two GM samples, with variable proportions of its segments. In conclusion, our study demonstrates the potential of plant transcriptomes for virome research, highlighting their strengths and limitations.

Reconstruction of Viral Variants via Monte Carlo Clustering

Identifying viral variants through clustering is essential for understanding the composition and structure of viral populations within and between hosts, which play a crucial role in disease progression and epidemic spread. This article proposes and validates novel Monte Carlo (MC) methods for clustering aligned viral sequences by minimizing either entropy or Hamming distance from consensuses. We validate these methods on four benchmarks: two SARS-CoV-2 interhost data sets and two HIV intrahost data sets. A parallelized version of our tool is scalable to very large data sets. We show that both entropy and Hamming distance-based MC clusterings discern the meaningful information from sequencing data. The proposed clustering methods consistently converge to similar clusterings across different runs. Finally, we show that MC clustering improves reconstruction of intrahost viral population from sequencing data.

Viral infectious diseases severity: co-presence of transcriptionally active microbes (TAMs) can play an integral role for disease severity

Humans have been challenged by infectious diseases for all of their recorded history, and are continually being affected even today. Next-generation sequencing (NGS) has enabled identification of, i) culture independent microbes, ii) emerging disease-causing pathogens, and iii) understanding of the genome architecture. This, in turn, has highlighted that pathogen/s are not a monolith, and thereby allowing for the differentiation of the wide-ranging disease symptoms, albeit infected by a primary pathogen. The conventional 'one disease - one pathogen' paradigm has been positively revisited by considering limited yet important evidence of the co-presence of multiple transcriptionally active microbes (TAMs), potential pathogens, in various infectious diseases, including the COVID-19 pandemic. The ubiquitous microbiota presence inside humans gives reason to hypothesize that the microbiome, especially TAMs, contributes to disease etiology. Herein, we discuss current evidence and inferences on the co-infecting microbes particularly in the diseases caused by the RNA viruses - Influenza, Dengue, and the SARS-CoV-2. We have highlighted that the specific alterations in the microbial taxonomic abundances (dysbiosis) is functionally connected to the exposure of primary infecting pathogen/s. The microbial presence is intertwined with the differential host immune response modulating differential disease trajectories. The microbiota-host interactions have been shown to modulate the host immune responses to Influenza and SARS-CoV-2 infection, wherein the active commensal microbes are involved in the generation of virus-specific CD4 and CD8 T-cells following the influenza virus infection. Furthermore, COVID-19 dysbiosis causes an increase in inflammatory cytokines such as IL-6, TNF-α, and IL-1β, which might be one of the important predisposing factors for severe infection. Through this article, we aim to provide a comprehensive view of functional microbiomes that can have a significant regulatory impact on predicting disease severity (mild, moderate and severe), as well as clinical outcome (survival and mortality). This can offer fresh perspectives on the novel microbial biomarkers for stratifying patients for severe disease symptoms, disease prevention and augmenting treatment regimens.

Strainline: full-length de novo viral haplotype reconstruction from noisy long reads

Haplotype-resolved de novo assembly of highly diverse virus genomes is critical in prevention, control and treatment of viral diseases. Current methods either can handle only relatively accurate short read data, or collapse haplotype-specific variations into consensus sequence. Here, we present Strainline, a novel approach to assemble viral haplotypes from noisy long reads without a reference genome. Strainline is the first approach to provide strain-resolved, full-length de novo assemblies of viral quasispecies from noisy third-generation sequencing data. Benchmarking on simulated and real datasets of varying complexity and diversity confirm this novelty and demonstrate the superiority of Strainline.

TNet: Transmission Network Inference Using Within-Host Strain Diversity and its Application to Geographical Tracking of COVID-19 Spread

The inference of disease transmission networks is an important problem in epidemiology. One popular approach for building transmission networks is to reconstruct a phylogenetic tree using sequences from disease strains sampled from infected hosts and infer transmissions based on this tree. However, most existing phylogenetic approaches for transmission network inference are highly computationally intensive and cannot take within-host strain diversity into account. Here, we introduce a new phylogenetic approach for inferring transmission networks, TNet, that addresses these limitations. TNet uses multiple strain sequences from each sampled host to infer transmissions and is simpler and more accurate than existing approaches. Furthermore, TNet is highly scalable and able to distinguish between ambiguous and unambiguous transmission inferences. We evaluated TNet on a large collection of 560 simulated transmission networks of various sizes and diverse host, sequence, and transmission characteristics, as well as on 10 real transmission datasets with known transmission histories. Our results show that TNet outperforms two other recently developed methods, phyloscanner and SharpTNI, that also consider within-host strain diversity. We also applied TNet to a large collection of SARS-CoV-2 genomes sampled from infected individuals in many countries around the world, demonstrating how our inference framework can be adapted to accurately infer geographical transmission networks. TNet is freely available from https://compbio.engr.uconn.edu/software/TNet/.

Population Disequilibrium as Promoter of Adaptive Explorations in Hepatitis C Virus

Replication of RNA viruses is characterized by exploration of sequence space which facilitates their adaptation to changing environments. It is generally accepted that such exploration takes place mainly in response to positive selection, and that further diversification is boosted by modifications of virus population size, particularly bottleneck events. Our recent results with hepatitis C virus (HCV) have shown that the expansion in sequence space of a viral clone continues despite prolonged replication in a stable cell culture environment. Diagnosis of the expansion was based on the quantification of diversity indices, the occurrence of intra-population mutational waves (variations in mutant frequencies), and greater individual residue variations in mutant spectra than those anticipated from sequence alignments in data banks. In the present report, we review our previous results, and show additionally that mutational waves in amplicons from the NS5A-NS5B-coding region are equally prominent during HCV passage in the absence or presence of the mutagenic nucleotide analogues favipiravir or ribavirin. In addition, by extending our previous analysis to amplicons of the NS3- and NS5A-coding region, we provide further evidence of the incongruence between amino acid conservation scores in mutant spectra from infected patients and in the Los Alamos National Laboratory HCV data banks. We hypothesize that these observations have as a common origin a permanent state of HCV population disequilibrium even upon extensive viral replication in the absence of external selective constraints or changes in population size. Such a persistent disequilibrium—revealed by the changing composition of the mutant spectrum—may facilitate finding alternative mutational pathways for HCV antiviral resistance. The possible significance of our model for other genetically variable viruses is discussed.

Principles of dengue virus evolvability derived from genotype-fitness maps in human and mosquito cells

Dengue virus (DENV) cycles between mosquito and mammalian hosts. To examine how DENV populations adapt to these different host environments, we used serial passage in human and mosquito cell lines and estimated fitness effects for all single-nucleotide variants in these populations using ultra-deep sequencing. This allowed us to determine the contributions of beneficial and deleterious mutations to the collective fitness of the population. Our analysis revealed that the continuous influx of a large burden of deleterious mutations counterbalances the effect of rare, host-specific beneficial mutations to shape the path of adaptation. Beneficial mutations preferentially map to intrinsically disordered domains in the viral proteome and cluster to defined regions in the genome. These phenotypically redundant adaptive alleles may facilitate host-specific DENV adaptation. Importantly, the evolutionary constraints described in our simple system mirror trends observed across DENV and Zika strains, indicating it recapitulates key biophysical and biological constraints shaping long-term viral evolution.

Analysis of heterogeneous genomic samples using image normalization and machine learning

BACKGROUND

Analysis of heterogeneous populations such as viral quasispecies is one of the most challenging bioinformatics problems. Although machine learning models are becoming to be widely employed for analysis of sequence data from such populations, their straightforward application is impeded by multiple challenges associated with technological limitations and biases, difficulty of selection of relevant features and need to compare genomic datasets of different sizes and structures.

RESULTS

We propose a novel preprocessing approach to transform irregular genomic data into normalized image data. Such representation allows to restate the problems of classification and comparison of heterogeneous populations as image classification problems which can be solved using variety of available machine learning tools. We then apply the proposed approach to two important problems in molecular epidemiology: inference of viral infection stage and detection of viral transmission clusters using next-generation sequencing data. The infection staging method has been applied to HCV HVR1 samples collected from 108 recently and 257 chronically infected individuals. The SVM-based image classification approach achieved more than 95% accuracy for both recently and chronically HCV-infected individuals. Clustering has been performed on the data collected from 33 epidemiologically curated outbreaks, yielding more than 97% accuracy.

CONCLUSIONS

Sequence image normalization method allows for a robust conversion of genomic data into numerical data and overcomes several issues associated with employing machine learning methods to viral populations. Image data also help in the visualization of genomic data. Experimental results demonstrate that the proposed method can be successfully applied to different problems in molecular epidemiology and surveillance of viral diseases. Simple binary classifiers and clustering techniques applied to the image data are equally or more accurate than other models.

Amplicon-Based, Next-Generation Sequencing Approaches to Characterize Single Nucleotide Polymorphisms of Orthohantavirus Species

Whole-genome sequencing (WGS) of viruses from patient or environmental samples can provide tremendous insight into the epidemiology, drug resistance or evolution of a virus. However, we face two common hurdles in obtaining robust sequence information; the low copy number of viral genomes in specimens and the error introduced by WGS techniques. To optimize detection and minimize error in WGS of hantaviruses, we tested four amplification approaches and different amplicon pooling methods for library preparation and examined these preparations using two sequencing platforms, Illumina MiSeq and Oxford Nanopore Technologies MinION. First, we tested and optimized primers used for whole segment PCR or one kilobase amplicon amplification for even coverage using RNA isolated from the supernatant of virus-infected cells. Once optimized we assessed two sources of total RNA, virus-infected cells and supernatant from the virus-infected cells, with four variations of primer pooling for amplicons, and six different amplification approaches. We show that 99–100% genome coverage was obtained using a one-step RT-PCR reaction with one forward and reverse primer. Using a two-step RT-PCR with three distinct tiling approaches for the three genomic segments (vRNAs), we optimized primer pooling approaches for PCR amplification to achieve a greater number of aligned reads, average depth of genome, and genome coverage. The single nucleotide polymorphisms identified from MiSeq and MinION sequencing suggested intrinsic mutation frequencies of ~10⁻⁵-10⁻⁷ per genome and 10⁻⁴-10⁻⁵ per genome, respectively. We noted no difference in the coverage or accuracy when comparing WGS results with amplicons amplified from RNA extracted from infected cells or supernatant of these infected cells. Our results show that high-throughput diagnostics requiring the identification of hantavirus species or strains can be performed using MiSeq or MinION using a one-step approach. However, the two-step MiSeq approach outperformed the MinION in coverage depth and accuracy, and hence would be superior for assessment of genomes for epidemiology or evolutionary questions using the methods developed herein.

Quasi-species nature and differential gene expression of severe acute respiratory syndrome coronavirus 2 and phylogenetic analysis of a novel Iranian strain

A novel coronavirus related to severe acute respiratory syndrome virus, (SARS-CoV-2) is the causal agent of the COVID-19 pandemic. Despite the genetic mutations across the SARS-CoV-2 genome being recently investigated, its transcriptomic genetic polymorphisms at inter-host level and the viral gene expression level based on each Open Reading Frame (ORF) remains unclear. Using available High Throughput Sequencing (HTS) data and based on SARS-CoV-2 infected human transcriptomic data, this study presents a high-resolution map of SARS-CoV-2 single nucleotide polymorphism (SNP) hotspots in a viral population at inter-host level. Four throat swab samples from COVID-19 infected patients were pooled, with RNA-Seq read retrieved from SRA NCBI to detect 21 SNPs and a replacement across the SARS-CoV-2 genomic population. Twenty-two RNA modification sites on viral transcripts were identified that may cause inter-host genetic diversity of this virus. In addition, the canonical genomic RNAs of N ORF showed higher expression in transcriptomic data and reverse transcriptase quantitative PCR compared to other SARS-CoV-2 ORFs, indicating the importance of this ORF in virus replication or other major functions in virus cycle. Phylogenetic and ancestral sequence analyses based on the entire genome revealed that SARS-CoV-2 is possibly derived from a recombination event between SARS-CoV and Bat SARS-like CoV. Ancestor analysis of the isolates from different locations including Iran suggest shared Chinese ancestry. These results propose the importance of potential inter-host level genetic variations to the evolution of SARS-COV-2, and the formation of viral quasi-species. The RNA modifications discovered in this study may cause amino acid sequence changes in polyprotein, spike protein, product of ORF8 and nucleocapsid (N) protein, suggesting further insights to understanding the functional impacts of mutations in the life cycle and pathogenicity of SARS-CoV-2.

Dissimilar Conservation Pattern in Hepatitis C Virus Mutant Spectra, Consensus Sequences, and Data Banks

The influence of quasispecies dynamics on long-term virus diversification in nature is a largely unexplored question. Specifically, whether intra-host nucleotide and amino acid variation in quasispecies fit the variation observed in consensus sequences or data bank alignments is unknown. Genome conservation and dynamics simulations are used for the computational design of universal vaccines, therapeutic antibodies and pan-genomic antiviral agents. The expectation is that selection of escape mutants will be limited when mutations at conserved residues are required. This strategy assumes long-term (epidemiologically relevant) conservation but, critically, does not consider short-term (quasispecies-dictated) residue conservation. We calculated mutant frequencies of individual loci from mutant spectra of hepatitis C virus (HCV) populations passaged in cell culture and from infected patients. Nucleotide or amino acid conservation in consensus sequences of the same populations, or in the Los Alamos HCV data bank did not match residue conservation in mutant spectra. The results relativize the concept of sequence conservation in viral genetics and suggest that residue invariance in data banks is an insufficient basis for the design of universal viral ligands for clinical purposes. Our calculations suggest relaxed mutational restrictions during quasispecies dynamics, which may contribute to higher calculated short-term than long-term viral evolutionary rates.

Animal and human RNA viruses: genetic variability and ability to overcome vaccines

RNA viruses, in general, exhibit high mutation rates; this is mainly due to the low fidelity displayed by the RNA-dependent polymerases required for their replication that lack the proofreading machinery to correct misincorporated nucleotides and produce high mutation rates. This lack of replication fidelity, together with the fact that RNA viruses can undergo spontaneous mutations, results in genetic variants displaying different viral morphogenesis, as well as variation on their surface glycoproteins that affect viral antigenicity. This diverse viral population, routinely containing a variety of mutants, is known as a viral ‘quasispecies’. The mutability of their virions allows for fast evolution of RNA viruses that develop antiviral resistance and overcome vaccines much more rapidly than DNA viruses. This also translates into the fact that pathogenic RNA viruses, that cause many diseases and deaths in humans, represent the major viral group involved in zoonotic disease transmission, and are responsible for worldwide pandemics.

Long-term context-dependent genetic adaptation of the viral genetic cloud

MOTIVATION

RNA viruses generate a cloud of genetic variants within each host. This cloud contains high-frequency genotypes, and many rare variants. The dynamics of these variants is crucial to understand viral evolution and their effect on their host.

RESULTS

We use an experimental evolution system to show that the genetic cloud surrounding the Coxsackie virus master sequence slowly, but steadily, evolves over hundreds of generations. This movement is determined by strong context-dependent mutations, where the frequency and type of mutations are affected by neighboring positions, even in silent mutations. This context-dependent mutation pattern serves as a spearhead for the viral population's movement within the adaptive landscape and affects which new dominant variants will emerge. The non-local mutation patterns affect the mutated dinucleotide distribution, and eventually lead to a non-uniform dinucleotide distribution in the main viral sequence. We tested these results on other RNA viruses with similar conclusions.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Limited evolution of the yellow fever virus 17d in a mouse infection model

By infecting mice with the yellow fever virus vaccine strain 17D (YFV-17D; Stamaril®), the dose dependence and evolutionary consequences of neurotropic yellow fever infection was assessed. Highly susceptible AG129 mice were used to allow for a maximal/unlimited expansion of the viral populations. Infected mice uniformly developed neurotropic disease; the virus was isolated from their brains, plaque purified and sequenced. Viral RNA populations were overall rather homogenous [Shannon entropies 0−0.15]. The remaining, yet limited intra-host population diversity (0−11 nucleotide exchanges per genome) appeared to be a consequence of pre-existing clonal heterogeneities (quasispecies) of Stamaril®. In parallel, mice were infected with a molecular clone of YFV-17D which was in vivo launched from a plasmid. Such plasmid-launched YFV-17D had a further reduced and almost clonal evolution. The limited intra-host evolution during unrestricted expansion in a highly susceptible host is relevant for vaccine and drug development against flaviviruses in general. Firstly, a propensity for limited evolution even upon infection with a (very) low inoculum suggests that fractional dosing as implemented in current YF-outbreak control may pose only a limited risk of reversion to pathogenic vaccine-derived virus variants. Secondly, it also largely lowers the chance of antigenic drift and development of resistance to antivirals.

Dynamic evolution in the key honey bee pathogen deformed wing virus: Novel insights into virulence and competition using reverse genetics

The impacts of invertebrate RNA virus population dynamics on virulence and infection outcomes are poorly understood. Deformed wing virus (DWV), the main viral pathogen of honey bees, negatively impacts bee health, which can lead to colony death. Despite previous reports on the reduction of DWV diversity following the arrival of the parasitic mite Varroa destructor, the key DWV vector, we found high genetic diversity of DWV in infested United States honey bee colonies. Phylogenetic analysis showed that divergent US DWV genotypes are of monophyletic origin and were likely generated as a result of diversification after a genetic bottleneck. To investigate the population dynamics of this divergent DWV, we designed a series of novel infectious cDNA clones corresponding to coexisting DWV genotypes, thereby devising a reverse-genetics system for an invertebrate RNA virus quasispecies. Equal replication rates were observed for all clone-derived DWV variants in single infections. Surprisingly, individual clones replicated to the same high levels as their mixtures and even the parental highly diverse natural DWV population, suggesting that complementation between genotypes was not required to replicate to high levels. Mixed clone–derived infections showed a lack of strong competitive exclusion, suggesting that the DWV genotypes were adapted to coexist. Mutational and recombination events were observed across clone progeny, providing new insights into the forces that drive and constrain virus diversification. Accordingly, our results suggest that Varroa influences DWV dynamics by causing an initial selective sweep, which is followed by virus diversification fueled by negative frequency-dependent selection for new genotypes. We suggest that this selection might reflect the ability of rare lineages to evade host defenses, specifically antiviral RNA interference (RNAi). In support of this hypothesis, we show that RNAi induced against one DWV strain is less effective against an alternate strain from the same population.

Deformed wing virus, a key pathogen of honey bees, shows rapid diversification after genetic bottlenecks; a novel reverse-genetic system provides insights into the forces that shape virus diversity, suggesting that virus quasi-species diversification may be driven by selection of genotypes capable of evading host RNAi defences.

Viral quasispecies

Viral quasispecies refers to a population structure that consists of extremely large numbers of variant genomes, termed mutant spectra, mutant swarms or mutant clouds. Fueled by high mutation rates, mutants arise continually, and they change in relative frequency as viral replication proceeds. The term quasispecies was adopted from a theory of the origin of life in which primitive replicons) consisted of mutant distributions, as found experimentally with present day RNA viruses. The theory provided a new definition of wild type, and a conceptual framework for the interpretation of the adaptive potential of RNA viruses that contrasted with classical studies based on consensus sequences. Standard clonal analyses and deep sequencing methodologies have confirmed the presence of myriads of mutant genomes in viral populations, and their participation in adaptive processes. The quasispecies concept applies to any biological entity, but its impact is more evident when the genome size is limited and the mutation rate is high. This is the case of the RNA viruses, ubiquitous in our biosphere, and that comprise many important pathogens. In virology, quasispecies are defined as complex distributions of closely related variant genomes subjected to genetic variation, competition and selection, and that may act as a unit of selection. Despite being an integral part of their replication, high mutation rates have an upper limit compatible with inheritable information. Crossing such a limit leads to RNA virus extinction, a transition that is the basis of an antiviral design termed lethal mutagenesis.

In vivo evolution of viral virulence: switching of deformed wing virus between hosts results in virulence changes and sequence shifts

The health of the Western honey bee is threatened by a global epidemic of deformed wing virus (DWV) infections driven by the ectoparasitic mite Varroa destructor acting as mechanical and biological virus vector. Three different variants of DWV, DWV-A, -B and -C exist. Virulence differences between these variants and their relation to V. destructor are still controversially discussed. We performed laboratory experiments to analyze the virulence of DWV directly isolated from crippled bees (DWV_P0 ) or after one additional passage in bee pupae (DWV_P1 ). We demonstrated that DWV_P0 was more virulent than DWV_P1 for pupae, when pupal mortality was taken as virulence marker, and for adult bees, when neurotropism and cognitive impairment were taken as virulence markers. Phylogenetic analysis supported that DWV exists as quasispecies and showed that DWV_P0 clustered with DWV-B and DWV_P1 with DWV-A when the phylogeny was based on the master sequences of the RNA-dependent RNA polymerase but not so when it was based on the VP3 region master sequences. We propose that switching of DWV between the bee and the mite host is accompanied by changes in viral sequence, tissue tropism and virulence and that the RNA-dependent RNA polymerase is involved in determining host range and virulence.

Evaluating the performance of tools used to call minority variants from whole genome short-read data

Background: High-throughput whole genome sequencing facilitates investigation of minority virus sub-populations from virus positive samples. Minority variants are useful in understanding within and between host diversity, population dynamics and can potentially assist in elucidating person-person transmission pathways. Several minority variant callers have been developed to describe low frequency sub-populations from whole genome sequence data. These callers differ based on bioinformatics and statistical methods used to discriminate sequencing errors from low-frequency variants. Methods: We evaluated the diagnostic performance and concordance between published minority variant callers used in identifying minority variants from whole-genome sequence data from virus samples. We used the ART-Illumina read simulation tool to generate three artificial short-read datasets of varying coverage and error profiles from an RSV reference genome. The datasets were spiked with nucleotide variants at predetermined positions and frequencies. Variants were called using FreeBayes, LoFreq, Vardict, and VarScan2. The variant callers' agreement in identifying known variants was quantified using two measures; concordance accuracy and the inter-caller concordance. Results: The variant callers reported differences in identifying minority variants from the datasets. Concordance accuracy and inter-caller concordance were positively correlated with sample coverage. FreeBayes identified the majority of variants although it was characterised by variable sensitivity and precision in addition to a high false positive rate relative to the other minority variant callers and which varied with sample coverage. LoFreq was the most conservative caller. Conclusions: We conducted a performance and concordance evaluation of four minority variant calling tools used to identify and quantify low frequency variants. Inconsistency in the quality of sequenced samples impacts on sensitivity and accuracy of minority variant callers. Our study suggests that combining at least three tools when identifying minority variants is useful in filtering errors when calling low frequency variants.

Evaluating the performance of tools used to call minority variants from whole genome short-read data

Background: High-throughput whole genome sequencing facilitates investigation of minority sub-populations from virus positive samples. Minority variants are useful in understanding within and between host diversity, population dynamics and can potentially help to elucidate person-person transmission chains. Several minority variant callers have been developed to describe the minority variants sub-populations from whole genome sequence data. However, they differ on bioinformatics and statistical approaches used to discriminate sequencing errors from low-frequency variants. Methods: We evaluated the diagnostic performance and concordance between published minority variant callers used in identifying minority variants from whole-genome sequence data. The ART-Illumina read simulation tool was used to generate three artificial short-read datasets of varying coverage and error profiles from an RSV reference genome. The datasets were spiked with nucleotide variants at predetermined positions and frequencies. Variants were called using FreeBayes, LoFreq, Vardict, and VarScan2. The variant callers' agreement in identifying known variants was quantified using two measures; concordance accuracy and the inter-caller concordance. Results: The variant callers reported differences in identifying minority variants from the datasets. Concordance accuracy and inter-caller concordance were positively correlated with sample coverage. FreeBayes identified majority of the variants although it was characterised by variable sensitivity and precision in addition to a high false positive rate relative to the other minority variant callers and which varied with sample coverage. LoFreq was the most conservative caller. Conclusions: We conducted a performance and concordance evaluation of four minority variant calling tools used to identify and quantify low frequency variants. Inconsistency in the quality of sequenced samples impact on sensitivity and accuracy of minority variant callers. Our study suggests that combining at least three tools when identifying minority variants is useful in filtering errors when calling low frequency variants.

5' UTR and NS5B-based genotyping of hepatitis C virus in patients from Damietta governorate, Egypt

Chronic hepatitis C virus (HCV) infection is a main health problem in Egypt causing high rates of mortalities. Egypt has the highest HCV prevalence in the world, with specific HCV subtypes epidemic and circulating extensively in the country. Different antiviral therapy protocols have been implemented for treating Egyptian HCV patients. Due to the limited data about HCV in Egypt, this study aimed to genotype HCV strains circulating in the Nile Delta Damietta governorate and to investigate the variation in the nonstructural 5B (NS5B) region targeted by the newly approved antiviral drugs. Thirty HCV samples from treatment-naïve patients were genotyped by restriction fragment length polymorphism. Some samples were genotyped by direct sequencing of their 5' untranslated region (UTR) and NS5B regions. Phylogenetic analysis was also performed on the sequences of their NS5B regions. Fourteen new sequences have been deposited in the GenBank database. Results showed that subtype 4a was prevalent in addition to subtype 1g. None of the previously reported NS5B substitutions were detected in the sequenced isolates from treatment-naïve patients, which may be a good predictor for efficient treatment of HCV Egyptian patients with Sofosbuvir. Further studies on Sofosbuvir treated-HCV Egyptian patients are required to investigate whether any NS5B substitutions can confer resistance to treatment.

Quasispecies composition and evolution of a typical Zika virus clinical isolate from Suriname

The arthropod-borne Zika virus (ZIKV) is currently causing a major international public health threat in the Americas. This study describes the isolation of ZIKV from the plasma of a 29-year-old female traveler that developed typical symptoms, like rash, fever and headache upon return from Suriname. The complete genome sequence including the 5′ and 3′ untranslated regions was determined and phylogenetic analysis showed the isolate clustering within the Asian lineage, close to other viruses that have recently been isolated in the Americas. In addition, the viral quasispecies composition was analyzed by single molecule real time sequencing, which suggested a mutation frequency of 1.4 × 10⁻⁴ for this ZIKV isolate. Continued passaging of the virus in cell culture led to the selection of variants with mutations in NS1 and the E protein. The latter might influence virus binding to cell surface heparan sulfate.

Analysis of Quasispecies of Avain Leukosis Virus Subgroup J Using Sanger and High-throughput Sequencing

Background

Avian leukosis viruses subgroup J (ALV-J) exists as a complex mixture of different, but closely related genomes named quasispecies subjected to continuous change according to the Principles of Darwinian evolution.

Method

The present study seeks to compare conventional Sanger sequencing with deep sequencing using MiSeq platform to study quasispecies dynamics of ALV-J.

Results

The accuracy and reproducibility of MiSeq sequencing was determined better than Sanger sequencing by running each experiment in duplicate. According to the mutational rate of single position and the ability to distinguish dominant quasispecies with two sequencing methods, conventional Sanger sequencing technique displayed high randomness due to few sequencing samples, while deep sequencing could reflect the composition of the quasispecies more accurately. In the mean time, the research of quasispecies via Sanger sequencing was simulated and analyzed with the aid of re-sampling strategy with replacement for 1000 times repeat from high-throughput sequencing data, which indicated that the higher antibody titer, the higher sequence entropy, the harder analyzing with the conventional Sanger sequencing, resulted in lower ratios of dominant variants.

Conclusions

In sum, deep sequencing is better suited for detecting rare variants comprehensively. The simulation of Sanger sequencing that we propose here will also help to standardize quasispecies researching under different selection pressure based on next-generation sequencing data.

Constraints on the Genetic and Antigenic Variability of Measles Virus

Antigenic drift and genetic variation are significantly constrained in measles virus (MeV). Genetic stability of MeV is exceptionally high, both in the lab and in the field, and few regions of the genome allow for rapid genetic change. The regions of the genome that are more tolerant of mutations (i.e., the untranslated regions and certain domains within the N, C, V, P, and M proteins) indicate genetic plasticity or structural flexibility in the encoded proteins. Our analysis reveals that strong constraints in the envelope proteins (F and H) allow for a single serotype despite known antigenic differences among its 24 genotypes. This review describes some of the many variables that limit the evolutionary rate of MeV. The high genomic stability of MeV appears to be a shared property of the Paramyxovirinae, suggesting a common mechanism that biologically restricts the rate of mutation.

Analysis of the Enzymatic Activity of an NS3 Helicase Genotype 3a Variant Sequence Obtained from a Relapse Patient

The hepatitis C virus (HCV) is a species of diverse genotypes that infect over 170 million people worldwide, causing chronic inflammation, cirrhosis and hepatocellular carcinoma. HCV genotype 3a is common in Brazil, and it is associated with a relatively poor response to current direct-acting antiviral therapies. The HCV NS3 protein cleaves part of the HCV polyprotein, and cellular antiviral proteins. It is therefore the target of several HCV drugs. In addition to its protease activity, NS3 is also an RNA helicase. Previously, HCV present in a relapse patient was found to harbor a mutation known to be lethal to HCV genotype 1b. The point mutation encodes the amino acid substitution W501R in the helicase RNA binding site. To examine how the W501R substitution affects NS3 helicase activity in a genotype 3a background, wild type and W501R genotype 3a NS3 alleles were sub-cloned, expressed in E. coli, and the recombinant proteins were purified and characterized. The impact of the W501R allele on genotype 2a and 3a subgenomic replicons was also analyzed. Assays monitoring helicase-catalyzed DNA and RNA unwinding revealed that the catalytic efficiency of wild type genotype 3a NS3 helicase was more than 600 times greater than the W501R protein. Other assays revealed that the W501R protein bound DNA less than 2 times weaker than wild type, and both proteins hydrolyzed ATP at similar rates. In Huh7.5 cells, both genotype 2a and 3a subgenomic HCV replicons harboring the W501R allele showed a severe defect in replication. Since the W501R allele is carried as a minor variant, its replication would therefore need to be attributed to the trans-complementation by other wild type quasispecies.

Nature and Extent of Genetic Diversity of Dengue Viruses Determined by 454 Pyrosequencing

Dengue virus (DENV) populations are characteristically highly diverse. Regular lineage extinction and replacement is an important dynamic DENV feature, and most DENV lineage turnover events are associated with increased incidence of disease. The role of genetic diversity in DENV lineage extinctions is not understood. We investigated the nature and extent of genetic diversity in the envelope (E) gene of DENV serotype 1 representing different lineages histories. A region of the DENV genome spanning the E gene was amplified and sequenced by Roche/454 pyrosequencing. The pyrosequencing results identified distinct sub-populations (haplotypes) for each DENV-1 E gene. A phylogenetic tree was constructed with the consensus DENV-1 E gene nucleotide sequences, and the sequences of each constructed haplotype showed that the haplotypes segregated with the Sanger consensus sequence of the population from which they were drawn. Haplotypes determined through pyrosequencing identified a recombinant DENV genome that could not be identified through Sanger sequencing. Nucleotide level sequence diversities of DENV-1 populations determined from SNP analysis were very low, estimated from 0.009–0.01. There were also no stop codon, frameshift or non-frameshift mutations observed in the E genes of any lineage. No significant correlations between the accumulation of deleterious mutations or increasing genetic diversity and lineage extinction were observed (p>0.5). Although our hypothesis that accumulation of deleterious mutations over time led to the extinction and replacement of DENV lineages was ultimately not supported by the data, our data does highlight the significant technical issues that must be resolved in the way in which population diversity is measured for DENV and other viruses. The results provide an insight into the within-population genetic structure and diversity of DENV-1 populations.

Widespread recombination, reassortment, and transmission of unbalanced compound viral genotypes in natural arenavirus infections

Arenaviruses are one of the largest families of human hemorrhagic fever viruses and are known to infect both mammals and snakes. Arenaviruses package a large (L) and small (S) genome segment in their virions. For segmented RNA viruses like these, novel genotypes can be generated through mutation, recombination, and reassortment. Although it is believed that an ancient recombination event led to the emergence of a new lineage of mammalian arenaviruses, neither recombination nor reassortment has been definitively documented in natural arenavirus infections. Here, we used metagenomic sequencing to survey the viral diversity present in captive arenavirus-infected snakes. From 48 infected animals, we determined the complete or near complete sequence of 210 genome segments that grouped into 23 L and 11 S genotypes. The majority of snakes were multiply infected, with up to 4 distinct S and 11 distinct L segment genotypes in individual animals. This S/L imbalance was typical: in all cases intrahost L segment genotypes outnumbered S genotypes, and a particular S segment genotype dominated in individual animals and at a population level. We corroborated sequencing results by qRT-PCR and virus isolation, and isolates replicated as ensembles in culture. Numerous instances of recombination and reassortment were detected, including recombinant segments with unusual organizations featuring 2 intergenic regions and superfluous content, which were capable of stable replication and transmission despite their atypical structures. Overall, this represents intrahost diversity of an extent and form that goes well beyond what has been observed for arenaviruses or for viruses in general. This diversity can be plausibly attributed to the captive intermingling of sub-clinically infected wild-caught snakes. Thus, beyond providing a unique opportunity to study arenavirus evolution and adaptation, these findings allow the investigation of unintended anthropogenic impacts on viral ecology, diversity, and disease potential.

Evaluation of canonical siRNA and Dicer substrate RNA for inhibition of hepatitis C virus genome replication--a comparative study

Hepatitis C virus (HCV) frequently establishes persistent infections in the liver, leading to the development of chronic hepatitis and, potentially, to liver cirrhosis and hepatocellular carcinoma at later stages. The objective of this study was to test the ability of five Dicer substrate siRNAs (DsiRNA) to inhibit HCV replication and to compare these molecules to canonical 21 nt siRNA. DsiRNA molecules were designed to target five distinct regions of the HCV genome - the 5' UTR and the coding regions for NS3, NS4B, NS5A or NS5B. These molecules were transfected into Huh7.5 cells that stably harboured an HCV subgenomic replicon expressing a firefly luciferase/neoR reporter (SGR-Feo-JFH-1) and were also tested on HCVcc-infected cells. All of the DsiRNAs inhibited HCV replication in both the subgenomic system and HCVcc-infected cells. When DsiRNAs were transfected prior to infection with HCVcc, the inhibition levels reached 99.5%. When directly compared, canonical siRNA and DsiRNA exhibited similar potency of virus inhibition. Furthermore, both types of molecules exhibited similar dynamics of inhibition and frequencies of resistant mutants after 21 days of treatment. Thus, DsiRNA molecules are as potent as 21 nt siRNAs for the inhibition of HCV replication and may provide future approaches for HCV therapy if the emergence of resistant mutants can be addressed.

Identification of mutations in the HVR1 and PKR-BD regions in HCV-infected patients resistant to PEG-IFNα/RBV therapy

The identification of mutations in the HVR1 region of hepatitis type C virus (HCV) is time-consuming and expensive, and there is a need for a rapid, inexpensive method of screening for these mutations to predict the ineffectiveness of pegylated interferon alpha combined with ribavirin (PEG-IFNα/RBV) therapy. The project was designed to evaluate the usefulness of the high resolution melting (HRM) technique to screen for mutation in the cDNAs encoding the HVR1 and protein kinase R-binding domain (PKR-BD) regions in a group of 36 patients infected with HCV and resistant to 12 months of combined therapy with PEG-IFNα/RBV. Viral RNA was isolated, reverse transcribed, and the fragments encoding the HVR1 and PKR-BD regions were polymerase chain reaction (PCR)-amplified, cloned, sequenced, and the melting profiles and the melting temperature (T_m) were determined by the HRM technique. After the treatment, the melting profiles of HVR1 cDNAs revealed a dominant peak corresponding to the T_m of about 85 °C (HCVs85) in almost all patients. One or more minor peaks were also observed, indicating the existence of cDNA(s) of different T_m. The HMR analysis suggested four typical forms of response to treatment. These suppositions were supported by sequencing. The HRM analysis revealed no changes in the melting profiles of PKR-BD cDNAs in the same patient before and after the therapy, suggesting that, within 12 months of treatment, new mutations were not introduced in PKR-BD. These findings were substantiated by sequencing. The HRM technique can be applied for the rapid screening for mutations in the cDNAs encoding the HVR and PKR-BD regions of HCV. We suggest that the detection of HCVs85 peak before the IFNα/RBV therapy might predict the ineffectiveness of treatment.

Quasispecies in population of compositional assemblies

BACKGROUND

The quasispecies model refers to information carriers that undergo self-replication with errors. A quasispecies is a steady-state population of biopolymer sequence variants generated by mutations from a master sequence. A quasispecies error threshold is a minimal replication accuracy below which the population structure breaks down. Theory and experimentation of this model often refer to biopolymers, e.g. RNA molecules or viral genomes, while its prebiotic context is often associated with an RNA world scenario. Here, we study the possibility that compositional entities which code for compositional information, intrinsically different from biopolymers coding for sequential information, could show quasispecies dynamics.

RESULTS

We employed a chemistry-based model, graded autocatalysis replication domain (GARD), which simulates the network dynamics within compositional molecular assemblies. In GARD, a compotype represents a population of similar assemblies that constitute a quasi-stationary state in compositional space. A compotype's center-of-mass is found to be analogous to a master sequence for a sequential quasispecies. Using single-cycle GARD dynamics, we measured the quasispecies transition matrix (Q) for the probabilities of transition from one center-of-mass Euclidean distance to another. Similarly, the quasispecies' growth rate vector (A) was obtained. This allowed computing a steady state distribution of distances to the center of mass, as derived from the quasispecies equation. In parallel, a steady state distribution was obtained via the GARD equation kinetics. Rewardingly, a significant correlation was observed between the distributions obtained by these two methods. This was only seen for distances to the compotype center-of-mass, and not to randomly selected compositions. A similar correspondence was found when comparing the quasispecies time dependent dynamics towards steady state. Further, changing the error rate by modifying basal assembly joining rate of GARD kinetics was found to display an error catastrophe, similar to the standard quasispecies model. Additional augmentation of compositional mutations leads to the complete disappearance of the master-like composition.

CONCLUSIONS

Our results show that compositional assemblies, as simulated by the GARD formalism, portray significant attributes of quasispecies dynamics. This expands the applicability of the quasispecies model beyond sequence-based entities, and potentially enhances validity of GARD as a model for prebiotic evolution.

Evaluation of anti-HIV-1 mutagenic nucleoside analogues

Because of their high mutation rates, RNA viruses and retroviruses replicate close to the threshold of viability. Their existence as quasi-species has pioneered the concept of "lethal mutagenesis" that prompted us to synthesize pyrimidine nucleoside analogues with antiviral activity in cell culture consistent with an accumulation of deleterious mutations in the HIV-1 genome. However, testing all potentially mutagenic compounds in cell-based assays is tedious and costly. Here, we describe two simple in vitro biophysical/biochemical assays that allow prediction of the mutagenic potential of deoxyribonucleoside analogues. The first assay compares the thermal stabilities of matched and mismatched base pairs in DNA duplexes containing or not the nucleoside analogues as follows. A promising candidate should display a small destabilization of the matched base pair compared with the natural nucleoside and the smallest gap possible between the stabilities of the matched and mismatched base pairs. From this assay, we predicted that two of our compounds, 5-hydroxymethyl-2'-deoxyuridine and 5-hydroxymethyl-2'-deoxycytidine, should be mutagenic. The second in vitro reverse transcription assay assesses DNA synthesis opposite nucleoside analogues inserted into a template strand and subsequent extension of the newly synthesized base pairs. Once again, only 5-hydroxymethyl-2'-deoxyuridine and 5-hydroxymethyl-2'-deoxycytidine are predicted to be efficient mutagens. The predictive potential of our fast and easy first line screens was confirmed by detailed analysis of the mutation spectrum induced by the compounds in cell culture because only compounds 5-hydroxymethyl-2'-deoxyuridine and 5-hydroxymethyl-2'-deoxycytidine were found to increase the mutation frequency by 3.1- and 3.4-fold, respectively.

Application of a high-resolution melting technique for the rapid detection of partial replacement of HCV-1b by HCV-1a after PEG-IFNα/RBV therapy

A modified method which can be used for the rapid screening of mutations in the protein kinase R-binding domain (PKR-BD) region and the hypervariable region 1 (HVR1) of hepatitis C virus (HCV) is described. This method is based on a high-resolution melting (HRM) technique used for genotyping single nucleotide polymorphisms and allows the detection of single nucleotide substitutions in the DNA sequence by measuring its T_m. The modified method, in addition to precisely measuring the T_m, allows the recording of the melting curve of the investigated cDNA fragment, which can provide provisional information about the number of different quasi-species present in the sample. The HRM analysis of the amplified cDNAs encoding the PKR-BD and HVR1 allowed the detection of partial replacement of HCV-1b by HCV-1a subspecies in one of our patients, as well as evaluation of the effectiveness of pegylated interferon α/ribavirin (PEG-IFNα/RBV) therapy. The HRM technique has never been used for the rapid screening of sequence variations in these regions and may be used for a similar purpose in any viral genome.

Current perspectives on HIV-1 antiretroviral drug resistance

Current advancements in antiretroviral therapy (ART) have turned HIV-1 infection into a chronic and manageable disease. However, treatment is only effective until HIV-1 develops resistance against the administered drugs. The most recent antiretroviral drugs have become superior at delaying the evolution of acquired drug resistance. In this review, the viral fitness and its correlation to HIV-1 mutation rates and drug resistance are discussed while emphasizing the concept of lethal mutagenesis as an alternative therapy. The development of resistance to the different classes of approved drugs and the importance of monitoring antiretroviral drug resistance are also summarized briefly.

One severe acute respiratory syndrome coronavirus protein complex integrates processive RNA polymerase and exonuclease activities

In addition to members causing milder human infections, the Coronaviridae family includes potentially lethal zoonotic agents causing severe acute respiratory syndrome (SARS) and the recently emerged Middle East respiratory syndrome. The ∼30-kb positive-stranded RNA genome of coronaviruses encodes a replication/transcription machinery that is unusually complex and composed of 16 nonstructural proteins (nsps). SARS-CoV nsp12, the canonical RNA-dependent RNA polymerase (RdRp), exhibits poorly processive RNA synthesis in vitro, at odds with the efficient replication of a very large RNA genome in vivo. Here, we report that SARS-CoV nsp7 and nsp8 activate and confer processivity to the RNA-synthesizing activity of nsp12. Using biochemical assays and reverse genetics, the importance of conserved nsp7 and nsp8 residues was probed. Whereas several nsp7 mutations affected virus replication to a limited extent, the replacement of two nsp8 residues (P183 and R190) essential for interaction with nsp12 and a third (K58) critical for the interaction of the polymerase complex with RNA were all lethal to the virus. Without a loss of processivity, the nsp7/nsp8/nsp12 complex can associate with nsp14, a bifunctional enzyme bearing 3'-5' exoribonuclease and RNA cap N7-guanine methyltransferase activities involved in replication fidelity and 5'-RNA capping, respectively. The identification of this tripartite polymerase complex that in turn associates with the nsp14 proofreading enzyme sheds light on how coronaviruses assemble an RNA-synthesizing machinery to replicate the largest known RNA genomes. This protein complex is a fascinating example of the functional integration of RNA polymerase, capping, and proofreading activities.

High rates of human immunodeficiency virus type 1 mutational profiles by single-genome amplification after 48-hour propagation in peripheral blood mononuclear cells at different levels of cell activation

Human immunodeficiency virus type 1 (HIV-1) genetic diversity is one of the most important features of HIV-1 infections and the result of error accumulation during reverse transcription and of high viral turnover. HIV-1 reverse transcription is influenced by factors such as the level of nucleotides and/or the cellular activation state. HIV-1 diversity was investigated after 48 h of viral propagation in peripheral blood mononuclear cells (PBMCs) obtained from healthy donors in three different cell culture conditions: (1) resting PBMCs, (2) simultaneous infection and PBMC activation, and (3) PBMC activation 72 h before infection. Cellular DNA was extracted and proviruses of each culture condition were amplified. Single-genome PCR clones were obtained and the protease and reverse transcriptase of the pol gene were sequenced. An elevated number of nucleotide substitutions in all three culture conditions were observed. In condition 1, the mutational rate observed ranged from 1.0 × 10(-3) to 2.1 × 10(-2), the genetic diversity was 0.6%, and hypermutation was observed in 7.1% of sequenced clones. In condition 2, the mutational rate ranged from 1.0 × 10(-3) to 1.0 × 10(-2), the genetic diversity was 0.8%, and hypermutation affected 6.7% of clones. In condition 3, the mutational rate ranged from 2.8 × 10(-3) to 1.1 × 10(-2), the genetic diversity was 1%, and 5.9% of clones were hypermutated. Substitutions occurred more frequently in some specific nucleotide stretches, and a common pattern for substitutions in all the different conditions was identified. There was a significant accumulation of mutations during the initial periods of in vitro HIV-1 propagation irrespective of culture conditions. The rapid accumulation of virus diversity might represent a viral strategy when colonizing new hosts. Complementary studies are necessary to allow for a better understanding of the initial periods of infection, which represent a crucial event related to disease progression.

Genetic structure and molecular variability of Cucumber mosaic virus isolates in the United States

Cucumber mosaic virus (CMV) has a worldwide distribution and the widest host range of any known plant virus. From 2000 to 2012, epidemics of CMV severely affected the production of snap bean (Phaseulos vulgaris L.) in the Midwest and Northeastern United States. Virus diversity leading to emergence of new strains is often considered a significant factor in virus epidemics. In addition to epidemics, new disease phenotypes arising from genetic exchanges or mutation can compromise effectiveness of plant disease management strategies. Here, we captured a snapshot of genetic variation of 32 CMV isolates collected from different regions of the U.S including new field as well as historic isolates. Nucleotide diversity (π) was low for U.S. CMV isolates. Sequence and phylogenetic analyses revealed that CMV subgroup I is predominant in the US and further showed that the CMV population is a mixture of subgroups IA and IB. Furthermore, phylogenetic analysis suggests likely reassortment between subgroups IA and IB within five CMV isolates. Based on phylogenetic and computational analysis, recombination between subgroups I and II as well as IA and IB in RNA 3 was detected. This is the first report of recombination between CMV subgroups I and II. Neutrality tests illustrated that negative selection was the major force operating upon the CMV genome, although some positively selected sites were detected for all encoded proteins. Together, these data suggest that different regions of the CMV genome are under different evolutionary constraints. These results also delineate composition of the CMV population in the US, and further suggest that recombination and reassortment among strain subgroups does occur but at a low frequency, and point towards CMV genomic regions that differ in types of selection pressure.

Comparison of test methodologies for foot-and-mouth disease virus serotype A vaccine matching

Vaccination has been one of the most important interventions in disease prevention and control. The impact of vaccination largely depends on the quality and suitability of the chosen vaccine. To determine the suitability of a vaccine strain, antigenic matching is usually studied by in vitro analysis. In this study, we performed three in vitro test methods to determine which one gives the lowest variability and the highest discriminatory capacity. Binary ethylenimine inactivated vaccines, prepared from 10 different foot-and-mouth disease (FMD) virus serotype A strains, were used to vaccinate cattle (5 animals for each strain). The antibody titers in blood serum samples 3 weeks postvaccination (w.p.v.) were determined by a virus neutralization test, neutralization index test, and liquid-phase blocking enzyme-linked immunosorbent assay (ELISA). The titers were then used to calculate relationship coefficient (r1) values. These r1 values were compared to the genetic lineage using receiver operating characteristic (ROC) analysis. In the two neutralization test methods, the median titers observed against the test strains differed considerably, and the sera of the vaccinated animals did not always show the highest titers against their respective homologous virus strains. When the titers were corrected for test strain effect (scaling), the variability (standard error of the mean per vaccinated group) increased because the results were on a different scale, but the discriminatory capacity improved. An ROC analysis of the r1 value calculated on both observed and scaled titers showed that only r1 values of the liquid-phase blocking ELISA gave a consistent statistically significant result. Under the conditions of the present study, the liquid-phase blocking ELISA showed less variation and still had a higher discriminatory capacity than the other tests.

Alphavirus mutator variants present host-specific defects and attenuation in mammalian and insect models

Arboviruses cycle through both vertebrates and invertebrates, which requires them to adapt to disparate hosts while maintaining genetic integrity during genome replication. To study the genetic mechanisms and determinants of these processes, we use chikungunya virus (CHIKV), a re-emerging human pathogen transmitted by the Aedes mosquito. We previously isolated a high fidelity (or antimutator) polymerase variant, C483Y, which had decreased fitness in both mammalian and mosquito hosts, suggesting this residue may be a key molecular determinant. To further investigate effects of position 483 on RNA-dependent RNA-polymerase (RdRp) fidelity, we substituted every amino acid at this position. We isolated novel mutators with decreased replication fidelity and higher mutation frequencies, allowing us to examine the fitness of error-prone arbovirus variants. Although CHIKV mutators displayed no major replication defects in mammalian cell culture, they had reduced specific infectivity and were attenuated in vivo. Unexpectedly, mutator phenotypes were suppressed in mosquito cells and the variants exhibited significant defects in RNA synthesis. Consequently, these replication defects resulted in strong selection for reversion during infection of mosquitoes. Since residue 483 is conserved among alphaviruses, we examined the analogous mutations in Sindbis virus (SINV), which also reduced polymerase fidelity and generated replication defects in mosquito cells. However, replication defects were mosquito cell-specific and were not observed in Drosophila S2 cells, allowing us to evaluate the potential attenuation of mutators in insect models where pressure for reversion was absent. Indeed, the SINV mutator variant was attenuated in fruit flies. These findings confirm that residue 483 is a determinant regulating alphavirus polymerase fidelity and demonstrate proof of principle that arboviruses can be attenuated in mammalian and insect hosts by reducing fidelity.

Chikungunya (CHIKV) is a re-emerging mosquito-borne virus that constitutes a major and growing human health burden. Like all RNA viruses, during viral replication CHIKV copies its genome using a polymerase that makes an average of one mistake per replication cycle. Therefore, a single virus generates millions of viral progeny that carry a multitude of distinct mutations in their genomes. In this study, we isolated CHIKV mutators (strains that make more errors than the wildtype virus), to study how higher mutation rates affect fitness in arthropod-borne viruses (arboviruses). CHIKV mutators have reduced virulence in mice and severe replication defects in Aedes mosquito cells. However, these replication defects result in selective pressure for reversion of mutators to a wildtype polymerase in mosquito hosts. To examine how mutators would behave in an insect model in absence of this genetic instability, we isolated mutators of a related virus, Sindbis virus (SINV). SINV mutators had no replication defect in fruit fly (Drosophila) cells, and a SINV mutator strain was stable and attenuated in fruit flies. This work shows proof of principle that arbovirus mutators can exhibit attenuation in both mammalian and insect hosts, and may remain a viable vaccine strategy.

GeneSV - an Approach to Help Characterize Possible Variations in Genomic and Protein Sequences

A computational approach for identification and assessment of genomic sequence variability (GeneSV) is described. For a given nucleotide sequence, GeneSV collects information about the permissible nucleotide variability (changes that potentially preserve function) observed in corresponding regions in genomic sequences, and combines it with conservation/variability results from protein sequence and structure-based analyses of evaluated protein coding regions. GeneSV was used to predict effects (functional vs. non-functional) of 37 amino acid substitutions on the NS5 polymerase (RdRp) of dengue virus type 2 (DENV-2), 36 of which are not observed in any publicly available DENV-2 sequence. 32 novel mutants with single amino acid substitutions in the RdRp were generated using a DENV-2 reverse genetics system. In 81% (26 of 32) of predictions tested, GeneSV correctly predicted viability of introduced mutations. In 4 of 5 (80%) mutants with double amino acid substitutions proximal in structure to one another GeneSV was also correct in its predictions. Predictive capabilities of the developed system were illustrated on dengue RNA virus, but described in the manuscript a general approach to characterize real or theoretically possible variations in genomic and protein sequences can be applied to any organism.

Comparison of the live attenuated yellow fever vaccine 17D-204 strain to its virulent parental strain Asibi by deep sequencing

BACKGROUND

The first comparison of a live RNA viral vaccine strain to its wild-type parental strain by deep sequencing is presented using as a model the yellow fever virus (YFV) live vaccine strain 17D-204 and its wild-type parental strain, Asibi.

METHODS

The YFV 17D-204 vaccine genome was compared to that of the parental strain Asibi by massively parallel methods. Variability was compared on multiple scales of the viral genomes. A modeled exploration of small-frequency variants was performed to reconstruct plausible regions of mutational plasticity.

RESULTS

Overt quasispecies diversity is a feature of the parental strain, whereas the live vaccine strain lacks diversity according to multiple independent measurements. A lack of attenuating mutations in the Asibi population relative to that of 17D-204 was observed, demonstrating that the vaccine strain was derived by discrete mutation of Asibi and not by selection of genomes in the wild-type population.

CONCLUSIONS

Relative quasispecies structure is a plausible correlate of attenuation for live viral vaccines. Analyses such as these of attenuated viruses improve our understanding of the molecular basis of vaccine attenuation and provide critical information on the stability of live vaccines and the risk of reversion to virulence.

The role of mutational robustness in RNA virus evolution

RNA viruses face dynamic environments and are masters at adaptation. During their short 'lifespans', they must surmount multiple physical, anatomical and immunological challenges. Central to their adaptative capacity is the enormous genetic diversity that characterizes RNA virus populations. Although genetic diversity increases the rate of adaptive evolution, low replication fidelity can present a risk because excess mutations can lead to population extinction. In this Review, we discuss the strategies used by RNA viruses to deal with the increased mutational load and consider how this mutational robustness might influence viral evolution and pathogenesis.

Arenavirus variations due to host-specific adaptation

Arenavirus particles are enveloped and contain two single-strand RNA genomic segments with ambisense coding. Genetic plasticity of the arenaviruses comes from transcription errors, segment reassortment, and permissive genomic packaging, and results in their remarkable ability, as a group, to infect a wide variety of hosts. In this review, we discuss some in vitro studies of virus genetic and phenotypic variation after exposure to selective pressures such as high viral dose, mutagens and antivirals. Additionally, we discuss the variation in vivo of selected isolates of Old World arenaviruses, particularly after infection of different animal species. We also discuss the recent emergence of new arenaviruses in the context of our observations of sequence variations that appear to be host-specific.

Longitudinal study of a human drug-induced model of autoantibody to cytoplasmic rods/rings following HCV therapy with ribavirin and interferon-α

Background

A novel pattern in the indirect immunofluorescence antinuclear antibody assay on HEp-2 cells (IIF-HEp-2) characterized by cytoplasmic rods and rings (RR) was reported in HCV patients, but stringent disease specificity studies and longitudinal analysis are lacking. We investigated the clinical significance of anti-RR in an HCV cohort with up to a 12-month treatment follow up.

Methodology/Results

597 patients (342 HCV, 55 HCV/HIV, 200 non-HCV) were screened and titered for anti-RR. Serial samples were available from 78 of 176 treated and 27 of 166 untreated patients. Anti-RR was detected in 14.1% of 342 HCV patients, 9.1% of 55 HCV/HIV, 3.4% of 29 Hepatitis B, and none of 171 non-HCV (p<0.0001; HCV versus non-HCV). Anti-RR was present in 38% of 108 patients receiving interferon-α/ribavirin, but none in 26 receiving either interferon-α or ribavirin, or 166 untreated patients (p<0.0001). Other IIF-HEp-2 patterns were more frequently associated with interferon-α treatment alone (52.2%) as compared to interferon-α/ribavirin (25%), ribavirin alone (33.3%), and no therapy (26.5%). Anti-RR frequency was not associated with sex, age, ethnicity, HCV genotype or viral load. Anti-RR occurred only after initiation of treatment, beginning as early as 1 month (6%), but by the sixth month >47% tested positive for anti-RR. The anti-RR titer generally increased with sustained treatment and remained high in 53% of patients. After treatment, anti-RR titer was negative in 41%. Non-responders to HCV therapy were 77% in anti-RR-positive versus 64% in anti-RR-negative patients. Response to treatment was not associated with anti-RR titer or the dynamics of anti-RR reactivity during and after treatment.

Conclusions

The exquisite association of anti-RR reactivity with combined interferon-α/ribavirin therapy in HCV patients represents a unique model for drug-induced autoantibody generation in humans as demonstrated by the fact that a significant fraction of patients who have anti-RR during therapy becomes anti-RR-negative after completion of therapy.

Emerging complexities of APOBEC3G action on immunity and viral fitness during HIV infection and treatment

The enzyme APOBEC3G (A3G) mutates the human immunodeficiency virus (HIV) genome by converting deoxycytidine (dC) to deoxyuridine (dU) on minus strand viral DNA during reverse transcription. A3G restricts viral propagation by degrading or incapacitating the coding ability of the HIV genome. Thus, this enzyme has been perceived as an innate immune barrier to viral replication whilst adaptive immunity responses escalate to effective levels. The discovery of A3G less than a decade ago led to the promise of new anti-viral therapies based on manipulation of its cellular expression and/or activity. The rationale for therapeutic approaches has been solidified by demonstration of the effectiveness of A3G in diminishing viral replication in cell culture systems of HIV infection, reports of its mutational footprint in virions from patients, and recognition of its unusually robust enzymatic potential in biochemical studies in vitro. Despite its effectiveness in various experimental systems, numerous recent studies have shown that the ability of A3G to combat HIV in the physiological setting is severely limited. In fact, it has become apparent that its mutational activity may actually enhance viral fitness by accelerating HIV evolution towards the evasion of both anti-viral drugs and the immune system. This body of work suggests that the role of A3G in HIV infection is more complex than heretofore appreciated and supports the hypothesis that HIV has evolved to exploit the action of this host factor. Here we present an overview of recent data that bring to light historical overestimation of A3G's standing as a strictly anti-viral agent. We discuss the limitations of experimental systems used to assess its activities as well as caveats in data interpretation.

Genetic relationships and epidemiological links between wild type 1 poliovirus isolates in Pakistan and Afghanistan

Background/Aim

Efforts have been made to eliminate wild poliovirus transmission since 1988 when the World Health Organization began its global eradication campaign. Since then, the incidence of polio has decreased significantly. However, serotype 1 and serotype 3 still circulate endemically in Pakistan and Afghanistan. Both countries constitute a single epidemiologic block representing one of the three remaining major global reservoirs of poliovirus transmission. In this study we used genetic sequence data to investigate transmission links among viruses from diverse locations during 2005-2007.

Methods

In order to find the origins and routes of wild type 1 poliovirus circulation, polioviruses were isolated from faecal samples of Acute Flaccid Paralysis (AFP) patients. We used viral cultures, two intratypic differentiation methods PCR, ELISA to characterize as vaccine or wild type 1 and nucleic acid sequencing of entire VP1 region of poliovirus genome to determine the genetic relatedness.

Results

One hundred eleven wild type 1 poliovirus isolates were subjected to nucleotide sequencing for genetic variation study. Considering the 15% divergence of the sequences from Sabin 1, Phylogenetic analysis by MEGA software revealed that active inter and intra country transmission of many genetically distinct strains of wild poliovirus type 1 belonged to genotype SOAS which is indigenous in this region. By grouping wild type 1 polioviruses according to nucleotide sequence homology, three distinct clusters A, B and C were obtained with multiple chains of transmission together with some silent circulations represented by orphan lineages.

Conclusion

Our results emphasize that there was a persistent transmission of wild type1 polioviruses in Pakistan and Afghanistan during 2005-2007. The epidemiologic information provided by the sequence data can contribute to the formulation of better strategies for poliomyelitis control to those critical areas, associated with high risk population groups which include migrants, internally displaced people, and refugees. The implication of this study is to maintain high quality mass immunization with oral polio vaccine (OPV) in order to interrupt chains of virus transmission in both countries to endorse substantial progress in Eastern-Mediterranean region.

Hepatitis C virus evasion mechanisms from neutralizing antibodies

Hepatitis C virus (HCV) represents a major public health problem, affecting 3% of the world's population. The majority of infected individuals develop chronic hepatitis, which can progress to cirrhosis and hepatocellular carcinoma. To date, a vaccine is not available and current therapy is limited by resistance, adverse effects and high costs. Although it is very well established that cell-mediated immunity is necessary for viral clearance, the importance of host antibodies in clearing HCV infection is being increasingly recognized. Indeed, recent studies indicate that neutralizing antibodies are induced in the early phase of infection by patients who subsequently clear viral infection. Conversely, patients who do not clear the virus develop high titers of neutralizing antibodies during the chronic stage. Surprisingly, these antibodies are not able to control HCV infection. HCV has therefore developed mechanisms to evade immune elimination, allowing it to persist in the majority of infected individuals. A detailed understanding of the mechanisms by which the virus escapes immune surveillance is therefore necessary if novel preventive and therapeutic treatments have to be designed. This review summarizes the current knowledge of the mechanisms used by HCV to evade host neutralizing antibodies.

Inferring viral quasispecies spectra from 454 pyrosequencing reads

BACKGROUND

RNA viruses infecting a host usually exist as a set of closely related sequences, referred to as quasispecies. The genomic diversity of viral quasispecies is a subject of great interest, particularly for chronic infections, since it can lead to resistance to existing therapies. High-throughput sequencing is a promising approach to characterizing viral diversity, but unfortunately standard assembly software was originally designed for single genome assembly and cannot be used to simultaneously assemble and estimate the abundance of multiple closely related quasispecies sequences.

RESULTS

In this paper, we introduce a new Viral Spectrum Assembler (ViSpA) method for quasispecies spectrum reconstruction and compare it with the state-of-the-art ShoRAH tool on both simulated and real 454 pyrosequencing shotgun reads from HCV and HIV quasispecies. Experimental results show that ViSpA outperforms ShoRAH on simulated error-free reads, correctly assembling 10 out of 10 quasispecies and 29 sequences out of 40 quasispecies. While ShoRAH has a significant advantage over ViSpA on reads simulated with sequencing errors due to its advanced error correction algorithm, ViSpA is better at assembling the simulated reads after they have been corrected by ShoRAH. ViSpA also outperforms ShoRAH on real 454 reads. Indeed, 7 most frequent sequences reconstructed by ViSpA from a real HCV dataset are viable (do not contain internal stop codons), and the most frequent sequence was within 1% of the actual open reading frame obtained by cloning and Sanger sequencing. In contrast, only one of the sequences reconstructed by ShoRAH is viable. On a real HIV dataset, ShoRAH correctly inferred only 2 quasispecies sequences with at most 4 mismatches whereas ViSpA correctly reconstructed 5 quasispecies with at most 2 mismatches, and 2 out of 5 sequences were inferred without any mismatches. ViSpA source code is available at http://alla.cs.gsu.edu/~software/VISPA/vispa.html.

CONCLUSIONS

ViSpA enables accurate viral quasispecies spectrum reconstruction from 454 pyrosequencing reads. We are currently exploring extensions applicable to the analysis of high-throughput sequencing data from bacterial metagenomic samples and ecological samples of eukaryote populations.

"Boom" and "Bust" cycles in virus growth suggest multiple selective forces in influenza a evolution

Background

Influenza A virus evolution in humans is driven at least in part by mutations allowing the virus to escape antibody neutralization. Little is known about the evolution of influenza in birds, a major reservoir of influenza A.

Methods

Neutralizing polyclonal antiserum was raised in chicken against reassortant influenza virus, CalX, bearing the hemagglutinin (HA) and neuraminidase (NA) of A/California/7/2004 [H3N2]. CalX was serially passaged in the presence of anti-CalX polyclonal IgY to derive viruses capable of growth in the presence of antibody.

Results

Polyclonal chicken antibody neutralized both HA activity and infection by CalX, but had no effect on a strain bearing an earlier human H3 and an irrelevant neuraminidase (A/Memphis/71-Bellamy/42 [H3N1]). Surprisingly, most of the antibody-resistant viruses were still at least partially sensitive to neutralization of HA activity and viral infection. Although mutant HA genes bearing changes that might affect antibody neutralization were identified, the vast majority of HA sequences obtained were identical to wild type, and no individual mutant sequence was found in more than one passage, suggesting that those mutations that were observed did not confer sufficient selective advantage to come to dominate the population. Different passages yielded infectious foci of varying size and plaques of varying size and morphology. Yields of infectious virus and relative frequency of different morphologies changed markedly from passage to passage. Sequences of bulk, uncloned PCR products from antibody-resistant passages indicated changes in the PB2 and PA proteins with respect to the wild type virus.

Conclusions

Each antibody-selected passage consisted of a variety of different cocirculating populations, rather than pure populations of virus able to escape antibody by changes in antibody epitopes. The ability to escape antibody is apparently due to changes in genes encoding the viral polymerase complex, probably resulting in more robust viral replication, allowing the few virus particles not completely neutralized by antibody to rapidly produce large numbers of progeny. Our data suggest that the relative success of an individual variant may depend on both its own gain and loss of fitness, as well as that of its cocirculating variants.

Identification of the 3' and 5' terminal sequences of the 8 rna genome segments of European and North American genotypes of infectious salmon anemia virus (an orthomyxovirus) and evidence for quasispecies based on the non-coding sequences of transcripts

Background

Infectious salmon anemia (ISA) virus (ISAV) is a pathogen of marine-farmed Atlantic salmon (Salmo salar); a disease first diagnosed in Norway in 1984. This virus, which was first characterized following its isolation in cell culture in 1995, belongs to the family Orthomyxoviridae, genus, Isavirus. The Isavirus genome consists of eight single-stranded RNA segments of negative sense, each with one to three open reading frames flanked by 3' and 5' non-coding regions (NCRs). Although the terminal sequences of other members of the family Orthomyxoviridae such as Influenzavirus A have been extensively analyzed, those of Isavirus remain largely unknown, and the few reported are from different ISAV strains and on different ends of the different RNA segments. This paper describes a comprehensive analysis of the 3' and 5' end sequences of the eight RNA segments of ISAV of both European and North American genotypes, and evidence of quasispecies of ISAV based on sequence variation in the untranslated regions (UTRs) of transcripts.

Results

Two different ISAV strains and two different RNA preparations were used in this study. ISAV strain ADL-PM 3205 ISAV-07 (ADL-ISAV-07) of European genotype was the source of total RNA extracted from ISAV-infected TO cells, which contained both viral mRNA and cRNA. ISAV strain NBISA01 of North American genotype was the source of vRNA extracted from purified virus. The NCRs of each segment were identified by sequencing cDNA prepared by three different methods, 5' RACE (Rapid amplification of cDNA ends), 3' RACE, and RNA ligation mediated PCR. Sequence analysis of five clones each derived from one RT-PCR product from each NCR of ISAV transcripts of segments 1 to 8 revealed significant heterogeneity among the clones of the same segment end, providing unequivocal evidence for presence of intra-segment ISAV quasispecies. Both RNA preparations (mRNA/cRNA and vRNA) yielded complementary sequence information, allowing the simultaneous identification and confirmation of the 3' and 5' NCR sequences of the 8 RNA genome segments of both genotypes of ISAV. The 3' sequences of the mRNA transcripts of ADL-ISAV-07 terminated 13-18 nucleotides from the full 3' terminus of cRNA, continuing as a poly(A) tail, which corresponded with the location of the polyadenylation signal. The lengths of the 3' and 5' NCRs of the vRNA were variable in the different genome segments, but the terminal 7 and 11 nucleotides of the 3' and 5' ends, respectively, were highly conserved among the eight genomic segments of ISAV. The first three nucleotides at the 3' end are GCU-3' (except in segment 5 with ACU-3'), whereas at the 5' end are 5'-AGU with the polyadenylation signal of 3-5 uridines 13-15 nucleotides downstream of the 5' end terminus of the vRNA. Exactly the same features were found in the respective complementary 5' and 3' end NCR sequences of the cRNA transcripts of ADL-ISAV-07, indicating that the terminal sequences of the 8 RNA genome segments are highly conserved among the two ISAV genotypes. The 5' NCR sequences of segments 1, 2, 3, 5, and 7, and the 3' NCR sequences of segments 3 and 4 cRNA were 100% identical in the two genotypes, and the 3' NCR sequences of segment 5 cRNA was the most divergent, with a sequence identity of 77.2%.

Conclusions

We report for the first time, the presence of intra-segment ISAV quasispecies, based on sequence variation in the NCR sequences of transcripts. In addition, this is the first report of a comprehensive unambiguous analysis of the 3' and 5' NCR sequences of all 8 RNA genome segments from two strains of ISAV representing the two genotypes of ISAV. Because most ISAV sequences are of cDNA to mRNA, they do not contain the 3' end sequences, which are removed during polyadenylation of the mRNA transcripts. We report for the first time the ISAV consensus sequence CA^T/_ATTTTTACT-3' (in the message sense 5'-3') in all segments of both ISAV genotypes.

Models of RNA virus evolution and their roles in vaccine design

Viruses are fast evolving pathogens that continuously adapt to the highly variable environments they live and reproduce in. Strategies devoted to inhibit virus replication and to control their spread among hosts need to cope with these extremely heterogeneous populations and with their potential to avoid medical interventions. Computational techniques such as phylogenetic methods have broadened our picture of viral evolution both in time and space, and mathematical modeling has contributed substantially to our progress in unraveling the dynamics of virus replication, fitness, and virulence. Integration of multiple computational and mathematical approaches with experimental data can help to predict the behavior of viral pathogens and to anticipate their escape dynamics. This piece of information plays a critical role in some aspects of vaccine development, such as viral strain selection for vaccinations or rational attenuation of viruses. Here we review several aspects of viral evolution that can be addressed quantitatively, and we discuss computational methods that have the potential to improve vaccine design.