Estimating Fitness of Viral Quasispecies from Next-Generation Sequencing Data

Genomic profile of eGFP-expressing canine distemper virus that undergoes serial plaque-to-plaque transfers

Canine distemper virus (CDV) is classified into the genus Morbillivirus in the family Paramyxoviridae. This virus has a single-stranded genomic RNA with negative polarity. The wild-type CDV genome is generally composed of 15 690 nucleotides. We previously rescued an enhanced green fluorescence protein (eGFP)-tagged recombinant CDV (rCDV-eGFP) using reverse genetics. In this study, the rCDV-eGFP at passage-7 was subjected to 38 serial plaque-to-plaque transfers (or bottleneck passages) and two extra common passages in cells. In theory, the effect of Muller’s ratchet may fix deleterious mutations in a single viral population after consecutive plaque-to-plaque transfers. In order to uncover a mutated landscape of the rCDV-eGFP under the circumstances of bottleneck passages, the passage-47 progeny was collected for the in-depth analysis via next-generation sequencing. The result revealed a total of nine single-nucleotide mutations (SNMs) in the viral antigenome. Out of them, SNMs at nt 1832, 5022, 5536, 5580, 5746, 6913 and 8803 were identified as total single-nucleotide substitution, i.e., 100% of mutation frequency. The result suggested no notable formation of viral quasispecies in the rCDV-eGFP population after consecutive plaque-to-plaque transfers.

A Two-Level, Intramutant Spectrum Haplotype Profile of Hepatitis C Virus Revealed by Self-Organized Maps

RNA viruses replicate as complex mutant spectra termed viral quasispecies. The frequency of each individual genome in a mutant spectrum depends on its rate of generation and its relative fitness in the replicating population ensemble. The advent of deep sequencing methodologies allows for the first-time quantification of haplotype abundances within mutant spectra. There is no information on the haplotype profile of the resident genomes and how the landscape evolves when a virus replicates in a controlled cell culture environment. Here, we report the construction of intramutant spectrum haplotype landscapes of three amplicons of the NS5A-NS5B coding region of hepatitis C virus (HCV). Two-dimensional (2D) neural networks were constructed for 44 related HCV populations derived from a common clonal ancestor that was passaged up to 210 times in human hepatoma Huh-7.5 cells in the absence of external selective pressures. The haplotype profiles consisted of an extended dense basal platform, from which a lower number of protruding higher peaks emerged. As HCV increased its adaptation to the cells, the number of haplotype peaks within each mutant spectrum expanded, and their distribution shifted in the 2D network. The results show that extensive HCV replication in a monotonous cell culture environment does not limit HCV exploration of sequence space through haplotype peak movements. The landscapes reflect dynamic variation in the intramutant spectrum haplotype profile and may serve as a reference to interpret the modifications produced by external selective pressures or to compare with the landscapes of mutant spectra in complex in vivo environments. IMPORTANCE The study provides for the first time the haplotype profile and its variation in the course of virus adaptation to a cell culture environment in the absence of external selective constraints. The deep sequencing-based self-organized maps document a two-layer haplotype distribution with an ample basal platform and a lower number of protruding peaks. The results suggest an inferred intramutant spectrum fitness landscape structure that offers potential benefits for virus resilience to mutational inputs.

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.

Dynamic Variation and Reversion in the Signature Amino Acids of H7N9 Virus During Human Infection

Background

Signature amino acids of H7N9 influenza A virus play critical roles in human adaption and pathogenesis, but their dynamic variation is unknown during disease development.

Methods

We sequentially collected respiratory samples from H7N9 patients at different timepoints and applied next-generation sequencing (NGS) to the whole genome of the H7N9 virus to investigate the variation at signature sites.

Results

A total of 11 patients were involved, from whom 29 samples were successfully sequenced, including samples from multiple timepoints in 9 patients. Neuraminidase (NA) R292K, basic polymerase 2 (PB2) E627K, and D701N were the 3 most dynamic mutations. The oseltamivir resistance-related NA R292K mutation was present in 9 samples from 5 patients, including 1 sample obtained before antiviral therapy. In all patients with the NA 292K mutation, the oseltamivir-sensitive 292R genotype persisted and was not eliminated by antiviral treatment. The PB2 E627K substitution was present in 18 samples from 8 patients, among which 12 samples demonstrated a mixture of E/K and the 627K frequency exhibited dynamic variation. Dual D701N and E627K mutations emerged but failed to achieve predominance in any of the samples.

Conclusions

Signature amino acids in PB2 and NA demonstrated high polymorphism and dynamic variation within individual patients during H7N9 virus infection.

Computational design of hepatitis C virus immunogens from host-pathogen dynamics over empirical viral fitness landscapes

Hepatitis C virus (HCV) afflicts 170 million people and kills 700 000 annually. Vaccination offers the most realistic and cost effective hope of controlling this epidemic, but despite 25 years of research, no vaccine is available. A major obstacle is HCV's extreme genetic variability and rapid mutational escape from immune pressure. Coupling maximum entropy inference with population dynamics simulations, we have employed a computational approach to translate HCV sequence databases into empirical landscapes of viral fitness and simulate the intrahost evolution of the viral quasispecies over these landscapes. We explicitly model the coupled host-pathogen dynamics by combining agent-based models of viral mutation with stochastically-integrated coupled ordinary differential equations for the host immune response. We validate our model in predicting the mutational evolution of the HCV RNA-dependent RNA polymerase (protein NS5B) within seven individuals for whom longitudinal sequencing data is available. We then use our approach to perform exhaustive in silico evaluation of putative immunogen candidates to rationally design tailored vaccines to simultaneously cripple viral fitness and block mutational escape within two selected individuals. By systematically identifying a small number of promising vaccine candidates, our empirical fitness landscapes and host-pathogen dynamics simulator can guide and accelerate experimental vaccine design efforts.

Rare haplotype load as marker for lethal mutagenesis

RNA viruses replicate with a template-copying fidelity, which lies close to an extinction threshold. Increases of mutation rate by nucleotide analogues can drive viruses towards extinction. This transition is the basis of an antiviral strategy termed lethal mutagenesis. We have introduced a new diversity index, the rare haplotype load (RHL), to describe NS5B (polymerase) mutant spectra of hepatitis C virus (HCV) populations passaged in absence or presence of the mutagenic agents favipiravir or ribavirin. The increase in RHL is more prominent in mutant spectra whose expansions were due to nucleotide analogues than to multiple passages in absence of mutagens. Statistical tests for paired mutagenized versus non-mutagenized samples with 14 diversity indices show that RHL provides consistently the highest standardized effect of mutagenic treatment difference for ribavirin and favipiravir. The results indicate that the enrichment of viral quasispecies in very low frequency minority genomes can serve as a robust marker for lethal mutagenesis. The diagnostic value of RHL from deep sequencing data is relevant to experimental studies on enhanced mutagenesis of viruses, and to pharmacological evaluations of inhibitors suspected to have a mutagenic activity.

A mutation-resistant deoxyribozyme OR gate for highly selective detection of viral nucleic acids

Highly selective probes hybridize only to fully complementary DNA or RNA sequences and, therefore, often fail to recognize mutated viral genomes. Here we designed a probe that possesses two seemingly incompatible properties: it tolerates some point mutations in genome, while it remains selective towards others. An OR deoxyribozyme logic gate was designed to fluorescently report the sequences of enterovirus 71 (EV71) covering ∼90% of all known EV71 strains. Importantly, sequences of closely related coxsackieviruses that differed by single nucleotides were reliably differentiated in 7 out of 8 cases.

Comparison of error correction algorithms for Ion Torrent PGM data: application to hepatitis B virus

Ion Torrent Personal Genome Machine (PGM) technology is a mid-length read, low-cost and high-speed next-generation sequencing platform with a relatively high insertion and deletion (indel) error rate. A full systematic assessment of the effectiveness of various error correction algorithms in PGM viral datasets (e.g., hepatitis B virus (HBV)) has not been performed. We examined 19 quality-trimmed PGM datasets for the HBV reverse transcriptase (RT) region and found a total error rate of 0.48% ± 0.12%. Deletion errors were clearly present at the ends of homopolymer runs. Tests using both real and simulated data showed that the algorithms differed in their abilities to detect and correct errors and that the error rate and sequencing depth significantly affected the performance. Of the algorithms tested, Pollux showed a better overall performance but tended to over-correct 'genuine' substitution variants, whereas Fiona proved to be better at distinguishing these variants from sequencing errors. We found that the combined use of Pollux and Fiona gave the best results when error-correcting Ion Torrent PGM viral data.

Internal Disequilibria and Phenotypic Diversification during Replication of Hepatitis C Virus in a Noncoevolving Cellular Environment

Viral quasispecies evolution upon long-term virus replication in a noncoevolving cellular environment raises relevant general issues, such as the attainment of population equilibrium, compliance with the molecular-clock hypothesis, or stability of the phenotypic profile. Here, we evaluate the adaptation, mutant spectrum dynamics, and phenotypic diversification of hepatitis C virus (HCV) in the course of 200 passages in human hepatoma cells in an experimental design that precluded coevolution of the cells with the virus. Adaptation to the cells was evidenced by increase in progeny production. The rate of accumulation of mutations in the genomic consensus sequence deviated slightly from linearity, and mutant spectrum analyses revealed a complex dynamic of mutational waves, which was sustained beyond passage 100. The virus underwent several phenotypic changes, some of which impacted the virus-host relationship, such as enhanced cell killing, a shift toward higher virion density, and increased shutoff of host cell protein synthesis. Fluctuations in progeny production and failure to reach population equilibrium at the genomic level suggest internal instabilities that anticipate an unpredictable HCV evolution in the complex liver environment.IMPORTANCE Long-term virus evolution in an unperturbed cellular environment can reveal features of virus evolution that cannot be explained by comparing natural viral isolates. In the present study, we investigate genetic and phenotypic changes that occur upon prolonged passage of hepatitis C virus (HCV) in human hepatoma cells in an experimental design in which host cell evolutionary change is prevented. Despite replication in a noncoevolving cellular environment, the virus exhibited internal population disequilibria that did not decline with increased adaptation to the host cells. The diversification of phenotypic traits suggests that disequilibria inherent to viral populations may provide a selective advantage to viruses that can be fully exploited in changing environments.

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.

On Eigen's Quasispecies Model, Two-Valued Fitness Landscapes, and Isometry Groups Acting on Finite Metric Spaces

A two-valued fitness landscape is introduced for the classical Eigen's quasispecies model. This fitness landscape can be considered as a direct generalization of the so-called single- or sharply peaked landscape. A general, non-permutation invariant quasispecies model is studied, and therefore the dimension of the problem is [Formula: see text], where N is the sequence length. It is shown that if the fitness function is equal to [Formula: see text] on a G-orbit A and is equal to w elsewhere, then the mean population fitness can be found as the largest root of an algebraic equation of degree at most [Formula: see text]. Here G is an arbitrary isometry group acting on the metric space of sequences of zeroes and ones of the length N with the Hamming distance. An explicit form of this exact algebraic equation is given in terms of the spherical growth function of the G-orbit A. Motivated by the analysis of the two-valued fitness landscapes, an abstract generalization of Eigen's model is introduced such that the sequences are identified with the points of a finite metric space X together with a group of isometries acting transitively on X. In particular, a simplicial analog of the original quasispecies model is discussed, which can be considered as a mathematical model of the switching of the antigenic variants for some bacteria.

Barrier-Independent, Fitness-Associated Differences in Sofosbuvir Efficacy against Hepatitis C Virus

Sofosbuvir displays a high phenotypic barrier to resistance, and it is a component of several combination therapies for hepatitis C virus (HCV) infections. HCV fitness can be a determinant of decreased sensitivity to direct-acting antiviral agents such as telaprevir or daclatasvir, but fitness-dependent decreased drug sensitivity has not been established for drugs with a high phenotypic barrier to resistance. Low- and high-fitness HCV populations and biological clones derived from them were used to infect Huh-7.5 hepatoma cells. Sofosbuvir efficacy was analyzed by measuring virus progeny production during several passages and by selection of possible sofosbuvir resistance mutations determined by sequencing the NS5B-coding region of the resulting populations. Sofosbuvir exhibited reduced efficacy against high-fitness HCV populations, without the acquisition of sofosbuvir-specific resistance mutations. A reduced sofosbuvir efficacy, similar to that observed with the parental populations, was seen for high-fitness individual biological clones. In independently derived high-fitness HCV populations or clones passaged in the presence of sofosbuvir, M289L was selected as the only substitution in the viral polymerase NS5B. In no case was the sofosbuvir-specific resistance substitution S282T observed. High HCV fitness can lead to decreased sensitivity to sofosbuvir, without the acquisition of specific sofosbuvir resistance mutations. Thus, fitness-dependent drug sensitivity can operate with HCV inhibitors that display a high barrier to resistance. This mechanism may underlie treatment failures not associated with selection of sofosbuvir-specific resistance mutations, linked to in vivo fitness of pretreatment viral populations.