Inference with viral quasispecies diversity indices: clonal and NGS approaches

Experimental evolution of an RNA virus in Caenorhabditis elegans

The discovery of Orsay virus (OrV), the first virus infecting wild populations of Caenorhabditis elegans, has boosted studies of viral immunity pathways in this nematode. Considering the many advantages that C. elegans offers for fundamental research in host-pathogen interactions, this pathosystem has high potential to become a model system for experimental virus evolution studies. However, the evolutionary constraints - i.e, the balance between genetic variation, selection, drift and historical contingency- operating in this pathosystem have barely been explored. Here we describe for the first time an evolution experiment of two different OrV strains in C. elegans. Comparison of the two ancestral strains showed differences in infectivity and sequence, and highlighted the importance of consistently normalize viral inocula for meaningful comparisons among strains. After 10 serial passages of evolution, we report slight changes in infectivity and non-synonymous mutations fixed in the evolved viral populations. In addition, we observed numerous minor variants emerging in the viral population. These minor variants were not randomly distributed along the genome but concentrated in polymorphic genomic regions. Overall, our work established the grounds for future experimental virus evolution studies using Caenorhabditis nematodes.

Virus Quasispecies Rarefaction: Subsampling with or without Replacement?

In quasispecies diversity studies, the comparison of two samples of varying sizes is a common necessity. However, the sensitivity of certain diversity indices to sample size variations poses a challenge. To address this issue, rarefaction emerges as a crucial tool, serving to normalize and create fairly comparable samples. This study emphasizes the imperative nature of sample size normalization in quasispecies diversity studies using next-generation sequencing (NGS) data. We present a thorough examination of resampling schemes using various simple hypothetical cases of quasispecies showing different quasispecies structures in the sense of haplotype genomic composition, offering a comprehensive understanding of their implications in general cases. Despite the big numbers implied in this sort of study, often involving coverages exceeding 100,000 reads per sample and amplicon, the rarefaction process for normalization should be performed with repeated resampling without replacement, especially when rare haplotypes constitute a significant fraction of interest. However, it is noteworthy that different diversity indices exhibit distinct sensitivities to sample size. Consequently, some diversity indicators may be compared directly without normalization, or instead may be resampled safely with replacement.

Extreme infectious titer variability in individual Aedes aegypti mosquitoes infected with Sindbis virus is associated with both differences in virus population structure and dramatic disparities in specific infectivity

Variability in how individuals respond to pathogens is a hallmark of infectious disease, yet the basis for individual variation in host response is often poorly understood. The titer of infectious virus among individual mosquitoes infected with arboviruses is frequently observed to vary by several orders of magnitude in a single experiment, even when the mosquitoes are highly inbred. To better understand the basis for this titer variation, we sequenced populations of Sindbis virus (SINV) obtained from individual infected Aedes aegypti mosquitoes that, despite being from a highly inbred laboratory colony, differed in their titers of infectious virus by approximately 10,000-fold. We observed genetic differences between these virus populations that indicated the virus present in the midguts of low titer mosquitoes was less fit than that of high titer mosquitoes, possibly due to founder effects that occurred during midgut infection. Furthermore, we found dramatic differences in the specific infectivity or SI (the ratio of infectious units/viral genome equivalents) between these virus populations, with the SI of low titer mosquitoes being up to 10,000-fold lower than that of high titer mosquitoes. Despite having similar amounts of viral genomes, low titer mosquitoes appeared to contain less viral particles, suggesting that viral genomes were packaged into virions less efficiently than in high titer mosquitoes. Finally, antibiotic treatment, which has been shown to suppress mosquito antiviral immunity, caused an increase in SI. Our results indicate that the extreme variation that is observed in SINV infectious titer between individual Ae. aegypti mosquitoes is due to both genetic differences between virus populations and to differences in the proportion of genomes that are packaged into infectious particles.

Puzzles, challenges, and information reservoir of SARS-CoV-2 quasispecies

Upon the emergence of SARS-CoV-2 in the human population, it was conjectured that for this coronavirus the dynamic intra-host heterogeneity typical of RNA viruses would be toned down. Nothing of this sort is observed. Here we review the main observations on the complexity and diverse composition of SARS-CoV-2 mutant spectra sampled from infected patients, within the framework of quasispecies dynamics. The analyses suggest that the information provided by myriads of genomic sequences within infected individuals may have a predictive value of the genomic sequences that acquire epidemiological relevance. Possibilities to reconcile the presence of broad mutant spectra in the large RNA coronavirus genome with its encoding a 3' to 5' exonuclease proofreading-repair activity are considered. Indeterminations in the behavior of individual viral genomes provide a benefit for the survival of the ensemble. We propose that this concept falls in the domain of "stochastic thinking," a notion that applies also to cellular processes, as a means for biological systems to face unexpected needs.

Fitness-Dependent, Mild Mutagenic Activity of Sofosbuvir for Hepatitis C Virus

The concept of a mild mutagen was coined to describe a minor mutagenic activity exhibited by some nucleoside analogues that potentiated their efficacy as antiretroviral agents. In the present study, we report the mild mutagen activity of sofosbuvir (SOF) for hepatitis C virus (HCV). Serial passages of HCV in human hepatoma cells, in the presence of SOF at a concentration well below its cytotoxic concentration 50 (CC50) led to pre-extinction populations whose mutant spectra exhibited a significant increase of C→U transitions, relative to populations passaged in the absence of SOF. This was reflected in an increase in several diversity indices that were used to characterize viral quasispecies. The mild mutagenic activity of SOF was largely absent when it was tested with isogenic HCV populations that displayed high replicative fitness. Thus, SOF can act as a mild mutagen for HCV, depending on HCV fitness. Possible mechanisms by which the SOF mutagenic activity may contribute to its antiviral efficacy are discussed.

The genetic variability and evolution of red-spotted grouper nervous necrosis virus quasispecies can be associated with its virulence

Nervous necrosis virus, NNV, is a neurotropic virus that causes viral nervous necrosis disease in a wide range of fish species, including European sea bass (Dicentrarchus labrax). NNV has a bisegmented (+) ssRNA genome consisting of RNA1, which encodes the RNA polymerase, and RNA2, encoding the capsid protein. The most prevalent NNV species in sea bass is red-spotted grouper nervous necrosis virus (RGNNV), causing high mortality in larvae and juveniles. Reverse genetics studies have associated amino acid 270 of the RGNNV capsid protein with RGNNV virulence in sea bass. NNV infection generates quasispecies and reassortants able to adapt to various selective pressures, such as host immune response or switching between host species. To better understand the variability of RGNNV populations and their association with RGNNV virulence, sea bass specimens were infected with two RGNNV recombinant viruses, a wild-type, rDl956, highly virulent to sea bass, and a single-mutant virus, Mut270Dl965, less virulent to this host. Both viral genome segments were quantified in brain by RT-qPCR, and genetic variability of whole-genome quasispecies was studied by Next Generation Sequencing (NGS). Copies of RNA1 and RNA2 in brains of fish infected with the low virulent virus were 1,000-fold lower than those in brains of fish infected with the virulent virus. In addition, differences between the two experimental groups in the Ts/Tv ratio, recombination frequency and genetic heterogeneity of the mutant spectra in the RNA2 segment were found. These results show that the entire quasispecies of a bisegmented RNA virus changes as a consequence of a single point mutation in the consensus sequence of one of its segments. Sea bream (Sparus aurata) is an asymptomatic carrier for RGNNV, thus rDl965 is considered a low-virulence isolate in this species. To assess whether the quasispecies characteristics of rDl965 were conserved in another host showing different susceptibility, juvenile sea bream were infected with rDl965 and analyzed as above described. Interestingly, both viral load and genetic variability of rDl965 in seabream were similar to those of Mut270Dl965 in sea bass. This result suggests that the genetic variability and evolution of RGNNV mutant spectra may be associated with its virulence.

Quasispecies Fitness Partition to Characterize the Molecular Status of a Viral Population. Negative Effect of Early Ribavirin Discontinuation in a Chronically Infected HEV Patient

The changes occurring in viral quasispecies populations during infection have been monitored using diversity indices, nucleotide diversity, and several other indices to summarize the quasispecies structure in a single value. In this study, we present a method to partition quasispecies haplotypes into four fractions according to their fitness: the master haplotype, rare haplotypes at two levels (those present at <0.1%, and those at 0.1−1%), and a fourth fraction that we term emerging haplotypes, present at frequencies >1%, but less than that of the master haplotype. We propose that by determining the changes occurring in the volume of the four quasispecies fitness fractions together with those of the Hill number profile we will be able to visualize and analyze the molecular changes in the composition of a quasispecies with time. To develop this concept, we used three data sets: a technical clone of the complete SARS-CoV-2 spike gene, a subset of data previously used in a study of rare haplotypes, and data from a clinical follow-up study of a patient chronically infected with HEV and treated with ribavirin. The viral response to ribavirin mutagenic treatment was selection of a rich set of synonymous haplotypes. The mutation spectrum was very complex at the nucleotide level, but at the protein (phenotypic/functional) level the pattern differed, showing a highly prevalent master phenotype. We discuss the putative implications of this observation in relation to mutagenic antiviral treatment.

SARS-CoV-2 Point Mutation and Deletion Spectra and Their Association with Different Disease Outcomes

Mutant spectra of RNA viruses are important to understand viral pathogenesis and response to selective pressures. There is a need to characterize the complexity of mutant spectra in coronaviruses sampled from infected patients. In particular, the possible relationship between SARS-CoV-2 mutant spectrum complexity and disease associations has not been established. In the present study, we report an ultradeep sequencing (UDS) analysis of the mutant spectrum of amplicons from the nsp12 (polymerase)- and spike (S)-coding regions of 30 nasopharyngeal isolates (diagnostic samples) of SARS-CoV-2 of the first COVID-19 pandemic wave (Madrid, Spain, April 2020) classified according to the severity of ensuing COVID-19. Low-frequency mutations and deletions, counted relative to the consensus sequence of the corresponding isolate, were overwhelmingly abundant. We show that the average number of different point mutations, mutations per haplotype, and several diversity indices was significantly higher in SARS-CoV-2 isolated from patients who developed mild disease than in those associated with moderate or severe disease (exitus). No such bias was observed with RNA deletions. Location of amino acid substitutions in the three-dimensional structures of nsp12 (polymerase) and S suggest significant structural or functional effects. Thus, patients who develop mild symptoms may be a richer source of genetic variants of SARS-CoV-2 than patients with moderate or severe COVID-19. IMPORTANCE The study shows that mutant spectra of SARS-CoV-2 from diagnostic samples differ in point mutation abundance and complexity and that significantly larger values were observed in virus from patients who developed mild COVID-19 symptoms. Mutant spectrum complexity is not a uniform trait among isolates. The nature and location of low-frequency amino acid substitutions present in mutant spectra anticipate great potential for phenotypic diversification of SARS-CoV-2.

Next-Generation Sequencing for Confronting Virus Pandemics

Virus pandemics have happened, are happening and will happen again. In recent decades, the rate of zoonotic viral spillover into humans has accelerated, mirroring the expansion of our global footprint and travel network, including the expansion of viral vectors and the destruction of natural spaces, bringing humans closer to wild animals. Once viral cross-species transmission to humans occurs, transmission cannot be stopped by cement walls but by developing barriers based on knowledge that can prevent or reduce the effects of any pandemic. Controlling a local transmission affecting few individuals is more efficient that confronting a community outbreak in which infections cannot be traced. Genetic detection, identification, and characterization of infectious agents using next-generation sequencing (NGS) has been proven to be a powerful tool allowing for the development of fast PCR-based molecular assays, the rapid development of vaccines based on mRNA and DNA, the identification of outbreaks, transmission dynamics and spill-over events, the detection of new variants and treatment of vaccine resistance mutations, the development of direct-acting antiviral drugs, the discovery of relevant minority variants to improve knowledge of the viral life cycle, strengths and weaknesses, the potential for becoming dominant to take appropriate preventive measures, and the discovery of new routes of viral transmission.

Analysis of HIV quasispecies and virological outcome of an HIV D+/R+ kidney-liver transplantation

INTRODUCTION

Transplantation among HIV positive patients may be a valuable therapeutic intervention. This study involves an HIV D+/R+ kidney-liver transplantation, where PBMC-associated HIV quasispecies were analyzed in donor and transplant recipients (TR) prior to transplantation and thereafter, together with standard viral monitoring.

METHODS

The donor was a 54 year of age HIV infected woman: kidney and liver recipients were two HIV infected men, aged 49 and 61. HIV quasispecies in PBMC was analyzed by ultra-deep sequencing of V3 env region. During TR follow-up, plasma HIV-1 RNA, HIV-1 DNA in PBMC, analysis of proviral integration sites and drug-resistance genotyping were performed. Other virological and immunological monitoring included CMV and EBV DNA quantification in blood and CD4 T cell counts.

RESULTS

Donor and TR were all ART-HIV suppressed at transplantation. Thereafter, TR maintained a nearly suppressed HIV-1 viremia, but HIV-1 RNA blips and the increase of proviral integration sites in PBMC attested some residual HIV replication. A transient peak in HIV-1 DNA occurred in the liver recipient. No major changes of drug-resistance genotype were detected after transplantation. CMV and EBV transient reactivations were observed only in the kidney recipient, but did not require specific treatment. CD4 counts remained stable. No intermixed quasispecies between donor and TR was observed at transplantation or thereafter. Despite signs of viral evolution in TR, HIV genetic heterogeneity did not increase over the course of the months of follow up.

CONCLUSIONS

No evidence of HIV superinfection was observed in the donor nor in the recipients. The immunosuppressive treatment administrated to TR did not result in clinical relevant viral reactivations.

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.

Host-dependent editing of SARS-CoV-2 in COVID-19 patients

A common trait among RNA viruses is their high capability to acquire genetic variability due to viral and host mechanisms. Next-generation sequencing (NGS) analysis enables the deep study of the viral quasispecies in samples from infected individuals. In this study, the viral quasispecies complexity and single nucleotide polymorphisms of the SARS-CoV-2 spike gene of coronavirus disease 2019 (COVID-19) patients with mild or severe disease were investigated using next-generation sequencing (Illumina platform). SARS-CoV-2 spike variability was higher in patients with long-lasting infection. Most substitutions found were present at frequencies lower than 1%, and had an A → G or T → C pattern, consistent with variants caused by adenosine deaminase acting on RNA-1 (ADAR1). ADAR1 affected a small fraction of replicating genomes, but produced multiple, mainly non-synonymous mutations. ADAR1 editing during replication rather than the RNA-dependent RNA polymerase (nsp12) was the predominant mechanism generating SARS-CoV-2 genetic variability. However, the mutations produced are not fixed in the infected human population, suggesting that ADAR1 may have an antiviral role, whereas nsp12-induced mutations occurring in patients with high viremia and persistent infection are the main source of new SARS-CoV-2 variants.

Metagenomic analysis of primary colorectal carcinomas and their metastases identifies potential microbial risk factors

The paucity of microbiome studies at intestinal tissues has contributed to a yet limited understanding of potential viral and bacterial cofactors of colorectal cancer (CRC) carcinogenesis or progression. We analysed whole-genome sequences of CRC primary tumours, their corresponding metastases and matched normal tissue for sequences of viral, phage and bacterial species. Bacteriome analysis showed Fusobacterium nucleatum, Streptococcus sanguinis, F. Hwasookii, Anaerococcus mediterraneensis and further species enriched in primary CRCs. The primary CRC of one patient was enriched for F. alocis, S. anginosus, Parvimonas micra and Gemella sp. 948. Enrichment of Escherichia coli strains IAI1, SE11, K-12 and M8 was observed in metastases together with coliphages enterobacteria phage φ80 and Escherichia phage VT2φ_272. Virome analysis showed that phages were the most preponderant viral species (46%), the main families being Myoviridae, Siphoviridae and Podoviridae. Primary CRCs were enriched for bacteriophages, showing five phages (Enterobacteria, Bacillus, Proteus, Streptococcus phages) together with their pathogenic hosts in contrast to normal tissues. The most frequently detected, and Blast-confirmed, viruses included human endogenous retrovirus K113, human herpesviruses 7 and 6B, Megavirus chilensis, cytomegalovirus (CMV) and Epstein-Barr virus (EBV), with one patient showing EBV enrichment in primary tumour and metastases. EBV was PCR-validated in 80 pairs of CRC primary tumour and their corresponding normal tissues; in 21 of these pairs (26.3%), it was detectable in primary tumours only. The number of viral species was increased and bacterial species decreased in CRCs compared with normal tissues, and we could discriminate primary CRCs from metastases and normal tissues by applying the Hutcheson t-test on the Shannon indices based on viral and bacterial species. Taken together, our results descriptively support hypotheses on microorganisms as potential (co)risk factors of CRC and extend putative suggestions on critical microbiome species in CRC metastasis.

Compartmentalized evolution of hepatitis B virus contributes differently to the prognosis of hepatocellular carcinoma

Serum hepatitis B virus (HBV) mutations can predict hepatocellular carcinoma (HCC) occurrence. We aimed to clarify if HBV evolves synchronously in the sera, adjacent liver and tumors and predict HCC prognosis equally. A total of 203 HBV-positive HCC patients with radical hepatectomy in Shanghai, China, during 2011-15 were enrolled in this prospective study. Quasispecies complexity (QC) in HBV core promoter region was assessed using clone-based sequencing. We performed RNA sequencing on tumors and paired adjacent tissues of another 15 HCC patients and analyzed it with three public data sets containing 127 samples. HBV QC was positively correlated to APOBEC3s' expression level (r = 0.28, P < 0.001), higher in the adjacent tissues than in the tumors (P = 6.50e-3), and higher in early tumors than in advanced tumors (P = 0.039). The evolutionary distance between the sera-derived HBV strains and the tumor-derived ones was significantly longer than that between the sera-derived ones and the adjacent tissue-derived ones (P < 0.001). Multivariate Cox regression analyses indicated that high HBV QC in the sera predicted an unfavorable overall survival (P = 0.002) and recurrence-free survival (RFS; P = 0.004) in HCC, whereas, in the tumors, it predicted a favorable RFS (P < 0.001). APOBECs-related HBV mutations, including G1764A, were more frequent in the sera than in the adjacent tissues. High-frequent A1762T/G1764A in the sera predicted an unfavorable RFS (P < 0.001), whereas, in the tumors, it predicted a favorable RFS (P = 0.035). In conclusion, HBV evolves more advanced in the sera than in the tumors. HBV QC and A1762T/G1764A in the sera and tumors have contrary prognostic effects in HCC.

Quantitative measures of within-host viral genetic diversity

The genetic diversity of virus populations within their hosts is known to influence disease progression, treatment outcome, drug resistance, cell tropism, and transmission risk, and the study of dynamic changes of genetic heterogeneity can provide insights into the evolution of viruses. Several measures to quantify within-host genetic diversity capturing different aspects of diversity patterns in a sample or population are used, based on incidence, relative frequencies, pairwise distances, or phylogenetic trees. Here, we review and compare several of these measures.

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.

Amino Acid Substitutions Associated with Treatment Failure for Hepatitis C Virus Infection

Despite the high virological response rates achieved with current directly acting antiviral agents (DAAs) against hepatitis C virus (HCV), around 2% to 5% of treated patients do not achieve a sustained viral response. The identification of amino acid substitutions associated with treatment failure requires analytical designs, such as subtype-specific ultradeep sequencing (UDS) methods, for HCV characterization and patient management. Using this procedure, we have identified six highly represented amino acid substitutions (HRSs) in NS5A and NS5B of HCV, which are not bona fide resistance-associated substitutions (RAS), from 220 patients who failed therapy. They were present frequently in basal and posttreatment virus of patients who failed different DAA-based therapies. Contrary to several RAS, HRSs belong to the acceptable subset of substitutions according to the PAM250 replacement matrix. Their mutant frequency, measured by the number of deep sequencing reads within the HCV quasispecies that encode the relevant substitutions, ranged between 90% and 100% in most cases. They also have limited predicted disruptive effects on the three-dimensional structures of the proteins harboring them. Possible mechanisms of HRS origin and dominance, as well as their potential predictive value for treatment response, are discussed.

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.

A beginner's guide for FMDV quasispecies analysis: sub-consensus variant detection and haplotype reconstruction using next-generation sequencing

Deep sequencing of viral genomes is a powerful tool to study RNA virus complexity. However, the analysis of next-generation sequencing data might be challenging for researchers who have never approached the study of viral quasispecies by this methodology. In this work we present a suitable and affordable guide to explore the sub-consensus variability and to reconstruct viral quasispecies from Illumina sequencing data. The guide includes a complete analysis pipeline along with user-friendly descriptions of software and file formats. In addition, we assessed the feasibility of the workflow proposed by analyzing a set of foot-and-mouth disease viruses (FMDV) with different degrees of variability. This guide introduces the analysis of quasispecies of FMDV and other viruses through this kind of approach.

Hidden genomic diversity of SARS-CoV-2: implications for qRT-PCR diagnostics and transmission

The COVID-19 pandemic has sparked an urgent need to uncover the underlying biology of this devastating disease. Though RNA viruses mutate more rapidly than DNA viruses, there are a relatively small number of single nucleotide polymorphisms (SNPs) that differentiate the main SARS-CoV-2 clades that have spread throughout the world. In this study, we investigated over 7,000 SARS-CoV-2 datasets to unveil both intrahost and interhost diversity. Our intrahost and interhost diversity analyses yielded three major observations. First, the mutational profile of SARS-CoV-2 highlights iSNV and SNP similarity, albeit with high variability in C>T changes. Second, iSNV and SNP patterns in SARS-CoV-2 are more similar to MERS-CoV than SARS-CoV-1. Third, a significant fraction of small indels fuel the genetic diversity of SARS-CoV-2. Altogether, our findings provide insight into SARS-CoV-2 genomic diversity, inform the design of detection tests, and highlight the potential of iSNVs for tracking the transmission of SARS-CoV-2.

Within-Host Viral Diversity: A Window into Viral Evolution

The evolutionary dynamics of a virus can differ within hosts and across populations. Studies of within-host evolution provide an important link between experimental studies of virus evolution and large-scale phylodynamic analyses. They can determine the extent to which global processes are recapitulated on local scales and how accurately experimental infections model natural ones. They may also inform epidemiologic models of disease spread and reveal how host-level dynamics contribute to a virus's evolution at a larger scale. Over the last decade, advances in viral sequencing have enabled detailed studies of viral genetic diversity within hosts. I review how within-host diversity is sampled, measured, and expressed, and how comparative studies of viral diversity can be leveraged to elucidate a virus's evolutionary dynamics. These concepts are illustrated with detailed reviews of recent research on the within-host evolution of influenza virus, dengue virus, and cytomegalovirus.

The Critical Role of Codon Composition on the Translation Efficiency Robustness of the Hepatitis A Virus Capsid

Hepatoviruses show an intriguing deviated codon usage, suggesting an evolutionary signature. Abundant and rare codons in the cellular genome are scarce in the human hepatitis A virus (HAV) genome, while intermediately abundant host codons are abundant in the virus. Genotype–phenotype maps, or fitness landscapes, are a means of representing a genotype position in sequence space and uncovering how genotype relates to phenotype and fitness. Using genotype–phenotype maps of the translation efficiency, we have shown the critical role of the HAV capsid codon composition in regulating translation and determining its robustness. Adaptation to an environmental perturbation such as the artificial induction of cellular shutoff—not naturally occurring in HAV infection—involved movements in the sequence space and dramatic changes of the translation efficiency. Capsid rare codons, including abundant and rare codons of the cellular genome, slowed down the translation efficiency in conditions of no cellular shutoff. In contrast, rare capsid codons that are abundant in the cellular genome were efficiently translated in conditions of shutoff. Capsid regions very rich in slowly translated codons adapt to shutoff through sequence space movements from positions with highly robust translation to others with diminished translation robustness. These movements paralleled decreases of the capsid physical and biological robustness, and resulted in the diversification of capsid phenotypes. The deviated codon usage of extant hepatoviruses compared with that of their hosts may suggest the occurrence of a virus ancestor with an optimized codon usage with respect to an unknown ancient host.

Broad and Dynamic Diversification of Infectious Hepatitis C Virus in a Cell Culture Environment

Previous studies documented that long-term hepatitis C virus (HCV) replication in human hepatoma Huh-7.5 cells resulted in viral fitness gain, expansion of the mutant spectrum, and several phenotypic alterations. In the present work, we show that mutational waves (changes in frequency of individual mutations) occurred continuously and became more prominent as the virus gained fitness. They were accompanied by an increasing proportion of heterogeneous genomic sites that affected 1 position in the initial HCV population and 19 and 69 positions at passages 100 and 200, respectively. Analysis of biological clones of HCV showed that these dynamic events affected infectious genomes, since part of the fluctuating mutations became incorporated into viable genomes. While 17 mutations were scored in 3 biological clones isolated from the initial population, the number reached 72 in 3 biological clones from the population at passage 200. Biological clones differed in their responses to antiviral inhibitors, indicating a phenotypic impact of viral dynamics. Thus, HCV adaptation to a specific constant environment (cell culture without external influences) broadens the mutant repertoire and does not focus the population toward a limited number of dominant genomes. A retrospective examination of mutant spectra of foot-and-mouth disease virus passaged in cell cultures suggests a parallel behavior here described for HCV. We propose that virus diversification in a constant environment has its basis in the availability of multiple alternative mutational pathways for fitness gain. This mechanism of broad diversification should also apply to other replicative systems characterized by high mutation rates and large population sizes.IMPORTANCE The study shows that extensive replication of an RNA virus in a constant biological environment does not limit exploration of sequence space and adaptive options. There was no convergence toward a restricted set of adapted genomes. Mutational waves and mutant spectrum broadening affected infectious genomes. Therefore, profound modifications of mutant spectrum composition and consensus sequence diversification are not exclusively dependent on environmental alterations or the intervention of population bottlenecks.

Extinction and emergence of genomic haplotypes during the evolution of Avian coronavirus in chicken embryos

Avian coronavirus (AvCoV) is ubiquitously present on poultry as a multitude of virus lineages. Studies on AvCoV phenotypic traits are dependent on the isolation of field strains in chicken embryonated eggs, but the mutant spectrum on each isolate is not considered. This manuscript reports the previously unknown HTS (high throughput sequencing)-based complete genome haplotyping of AvCoV isolates after passages of two field strains in chicken embryonated eggs. For the first and third passages of strain 23/2013, virus loads were 6.699 log copies/ μL and 6 log copies/ μL and, for 38/2013, 5.699 log copies/μL and 2.699 log copies/μL of reaction, respectively. The first passage of strain 23/2013 contained no variant haplotype, while, for the third passage, five putative variant haplotypes were found, with > 99.9% full genome identity with each other and with the dominant genome. Regarding strain 38/2013, five variant haplotypes were found for the first passage, with > 99.9% full genome identity with each other and with the dominant genome, and a single variant haplotype was found. Extinction and emergence of haplotypes with polymorphisms in genes involved in receptor binding and regulation of RNA synthesis were observed, suggesting that phenotypic traits of AvCoV isolates are a result of their mutant spectrum.

Characterization of hepatitis B virus X gene quasispecies complexity in mono-infection and hepatitis delta virus superinfection

BACKGROUND

Hepatitis delta virus (HDV) seems to strongly suppress hepatitis B virus (HBV) replication, although little is known about the mechanism of this interaction. Both these viruses show a dynamic distribution of mutants, resulting in viral quasispecies. Next-generation sequencing is a viable approach for analyzing the composition of these mutant spectra. As the regulatory hepatitis B X protein (HBx) is essential for HBV replication, determination of HBV X gene (HBX) quasispecies complexity in HBV/HDV infection compared to HBV mono-infection may provide information on the interactions between these two viruses.

AIM

To compare HBV quasispecies complexity in the HBX 5’ region between chronic hepatitis delta (CHD) and chronic HBV mono-infected patients.

METHODS

Twenty-four untreated patients were included: 7/24 (29.2%) with HBeAg-negative chronic HBV infection (CI, previously termed inactive carriers), 8/24 (33.3%) with HBeAg-negative chronic hepatitis B (CHB) and 9/24 (37.5%) with CHD. A serum sample from each patient was first tested for HBV DNA levels. The HBX 5’ region [nucleotides (nt) 1255-1611] was then PCR-amplified for subsequent next-generation sequencing (MiSeq, Illumina, United States). HBV quasispecies complexity in the region analyzed was evaluated using incidence-based indices (number of haplotypes and number of mutations), abundance-based indices (Hill numbers of order 1 and 2), and functional indices (mutation frequency and nucleotide diversity). We also evaluated the pattern of nucleotide changes to investigate which of them could be the cause of the quasispecies complexity.

RESULTS

CHB patients showed higher median HBV-DNA levels [5.4 logIU/mL, interquartile range (IQR) 3.5-7.9] than CHD (3.4 logIU/mL, IQR 3-7.6) (P = n.s.) or CI (3.2 logIU/mL, IQR 2.3-3.5) (P < 0.01) patients. The incidence and abundance indices indicated that HBV quasispecies complexity was significantly greater in CI than CHB. A similar trend was observed in CHD patients, although only Hill numbers of order 2 showed statistically significant differences (CHB 2.81, IQR 1.11-4.57 vs CHD 8.87, 6.56-11.18, P = 0.038). There were no significant differences in the functional indices, but CI and CHD patients also showed a trend towards greater complexity than CHB. No differences were found for any HBV quasispecies complexity indices between CHD and CI patients. G-to-A and C-to-T nucleotide changes, characteristic of APOBEC3G, were higher in CHD and CI than in CHB in genotype A haplotypes, but not in genotype D. The proportion of nt G-to-A vs A-to-G changes and C-to-T vs T-to-C changes in genotype A and D haplotypes in CHD patients showed no significant differences. In CHB and CI the results of these comparisons were dependent on HBV genotype.

CONCLUSION

The lower-replication CHD and CI groups show a trend to higher quasispecies complexity than the higher-replication CHB group. The mechanisms associated with this greater complexity require elucidation.

Hepatitis C virus intrinsic molecular determinants may contribute to the development of cholestatic hepatitis after liver transplantation

Cholestatic hepatitis C (CHC) is a severe form of hepatitis C virus (HCV) infection recurrence that leads to high graft loss rates early after liver transplantation (LT). To investigate the pathogenic mechanisms of CHC, we analysed HCV quasispecies in CHC patients compared to a control group (mild hepatitis C recurrence) by deep pyrosequencing. At the time of LT, NS5B quasispecies complexity was similar between the two groups but, after LT, it decreased more sharply in CHC patients than in the control group. Interestingly, the major variant before LT propagated efficiently and remained as the dominant sequence after LT in 62 % of CHC patients versus 11 % of controls (P=0.031). Sequence analysis of the complete non-structural region in a limited number of patients revealed a potential 12 aa signature specific to the CHC group. These data suggest that intrinsic molecular determinants in the circulating HCV quasispecies may provide a fitness advantage, contributing to the development of CHC.

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.

Resistance of high fitness hepatitis C virus to lethal mutagenesis

Viral fitness quantifies the degree of virus adaptation to a given environment. How viral fitness can influence the mutant spectrum complexity of a viral quasispecies subjected to lethal mutagenesis has not been investigated. Here we document that two high fitness hepatitis C virus populations display higher resistance to the mutagenic nucleoside analogues favipiravir and ribavirin than their parental, low fitness HCV. All populations, however, exhibited a mutation transition bias indicative of active mutagenesis. Resistance to the analogues was associated with a limited expansion of mutant spectrum complexity, as evidenced by several diversity indices used to characterize mutant spectra. The results are consistent with a replicative site-drug competition mechanism that was previously proposed for HCV fitness-associated resistance to non-mutagenic inhibitors. Other alternative, non-mutually exclusive mechanisms are considered. The results introduce viral fitness as a relevant parameter to evaluate the response of viruses to lethal mutagenesis, with implications for antiviral designs.

Characterization of intra- and inter-host norovirus P2 genetic variability in linked individuals by amplicon sequencing

Noroviruses are the main cause of epidemics of acute gastroenteritis at a global scale. Although chronically infected immunocompromised individuals are regarded as potential reservoirs for the emergence of new viral variants, viral quasispecies distribution and evolution patterns in acute symptomatic and asymptomatic infections has not been extensively studied. Amplicons of 450 nts from the P2 coding capsid domain were studied using next-generation sequencing (454/GS-Junior) platform. Inter-host diversity between symptomatic and asymptomatic acutely infected individuals linked to the same outbreak as well as their viral intra-host diversity over time were characterized. With an average of 2848 reads per sample and a cutoff frequency of 0.1%, minor variant haplotypes were detected in 5 out of 8 specimens. Transmitted variants could not be confirmed in all infected individuals in one outbreak. The observed initial intra-host viral diversity in asymptomatically infected subjects was higher than in symptomatic ones. Viral quasispecies evolution over time within individuals was host-specific, with an average of 2.8 nt changes per day (0.0062 changes per nucleotide per day) in a given symptomatic case. Nucleotide polymorphisms were detected in 28 out of 450 analyzed nucleotide positions, 32.14% of which were synonymous and 67.86% were non-synonymous. Most observed amino acid changes emerged at or near blockade epitopes A, B, D and E. Our results suggest that acutely infected individuals, even in the absence of symptoms, which go underreported and may enhance transmission, may contribute to norovirus genetic variability and evolution.

Intrapatient viral diversity and treatment outcome in patients with genotype 3a hepatitis C virus infection on sofosbuvir-containing regimens

Treatment with the direct-acting antiviral agent (DAA) sofosbuvir (SOF), an NS5B inhibitor, and velpatasvir (VEL), an NS5A inhibitor, demonstrates viral cure rates of ≥95% in hepatitis C virus (HCV) genotypes (GT) 1-6. Here, we investigated intrapatient HCV diversity in NS5A and NS5B using Shannon entropy to examine the relationship between viral diversity and treatment outcome. At baseline, HCV diversity was lowest in patients infected with HCV GT3 as compared to the other GTs, and viral diversity was greater in NS5A than NS5B (P < .0001). Treatment outcome with SOF/VEL or the comparator regimen of SOF with ribavirin (RBV) was not correlated with baseline diversity. However, among persons treated with SOF/VEL, a decrease in diversity from baseline was observed at relapse in the majority virologic failures, consistent with a viral bottleneck event at relapse. In contrast, an increase in diversity was observed in 27% of SOF+RBV virologic failures. We investigated whether the increase in diversity was due to an increase in the transition rate, one mode of potential RBV-mediated mutagenesis; however, we found no evidence of this mechanism. Overall, we did not observe that viral diversity at baseline influenced treatment outcome, but the diversity changes observed at relapse can improve our understanding of RBV viral suppression in vivo.

Deep sequencing in the management of hepatitis virus infections

The hepatitis viruses represent a major public health problem worldwide. Procedures for characterization of the genomic composition of their populations, accurate diagnosis, identification of multiple infections, and information on inhibitor-escape mutants for treatment decisions are needed. Deep sequencing methodologies are extremely useful for these viruses since they replicate as complex and dynamic quasispecies swarms whose complexity and mutant composition are biologically relevant traits. Population complexity is a major challenge for disease prevention and control, but also an opportunity to distinguish among related but phenotypically distinct variants that might anticipate disease progression and treatment outcome. Detailed characterization of mutant spectra should permit choosing better treatment options, given the increasing number of new antiviral inhibitors available. In the present review we briefly summarize our experience on the use of deep sequencing for the management of hepatitis virus infections, particularly for hepatitis B and C viruses, and outline some possible new applications of deep sequencing for these important human pathogens.

Improving virus production through quasispecies genomic selection and molecular breeding

Virus production still is a challenging issue in antigen manufacture, particularly with slow-growing viruses. Deep-sequencing of genomic regions indicative of efficient replication may be used to identify high-fitness minority individuals suppressed by the ensemble of mutants in a virus quasispecies. Molecular breeding of quasispecies containing colonizer individuals, under regimes allowing more than one replicative cycle, is a strategy to select the fittest competitors among the colonizers. A slow-growing cell culture-adapted hepatitis A virus strain was employed as a model for this strategy. Using genomic selection in two regions predictive of efficient translation, the internal ribosome entry site and the VP1-coding region, high-fitness minority colonizer individuals were identified in a population adapted to conditions of artificially-induced cellular transcription shut-off. Molecular breeding of this population with a second one, also adapted to transcription shut-off and showing an overall colonizer phenotype, allowed the selection of a fast-growing population of great biotechnological potential.

Lethal Mutagenesis of Hepatitis C Virus Induced by Favipiravir

Lethal mutagenesis is an antiviral approach that consists in extinguishing a virus by an excess of mutations acquired during replication in the presence of a mutagen. Here we show that favipiravir (T-705) is a potent mutagenic agent for hepatitis C virus (HCV) during its replication in human hepatoma cells. T-705 leads to an excess of G → A and C → U transitions in the mutant spectrum of preextinction HCV populations. Infectivity decreased significantly in the presence of concentrations of T-705 which are 2- to 8-fold lower than its cytotoxic concentration 50 (CC50). Passaging the virus five times in the presence of 400 μM T-705 resulted in virus extinction. Since T-705 has undergone advanced clinical trials for approval for human use, the results open a new approach based on lethal mutagenesis to treat hepatitis C virus infections. If proven effective for HCV in vivo, this new anti-HCV agent may be useful in patient groups that fail current therapeutic regimens.

Genetic Drift, Purifying Selection and Vector Genotype Shape Dengue Virus Intra-host Genetic Diversity in Mosquitoes

Due to their error-prone replication, RNA viruses typically exist as a diverse population of closely related genomes, which is considered critical for their fitness and adaptive potential. Intra-host demographic fluctuations that stochastically reduce the effective size of viral populations are a challenge to maintaining genetic diversity during systemic host infection. Arthropod-borne viruses (arboviruses) traverse several anatomical barriers during infection of their arthropod vectors that are believed to impose population bottlenecks. These anatomical barriers have been associated with both maintenance of arboviral genetic diversity and alteration of the variant repertoire. Whether these patterns result from stochastic sampling (genetic drift) rather than natural selection, and/or from the influence of vector genetic heterogeneity has not been elucidated. Here, we used deep sequencing of full-length viral genomes to monitor the intra-host evolution of a wild-type dengue virus isolate during infection of several mosquito genetic backgrounds. We estimated a bottleneck size ranging from 5 to 42 founding viral genomes at initial midgut infection, irrespective of mosquito genotype, resulting in stochastic reshuffling of the variant repertoire. The observed level of genetic diversity increased following initial midgut infection but significantly differed between mosquito genetic backgrounds despite a similar initial bottleneck size. Natural selection was predominantly negative (purifying) during viral population expansion. Taken together, our results indicate that dengue virus intra-host genetic diversity in the mosquito vector is shaped by genetic drift and purifying selection, and point to a novel role for vector genetic factors in the genetic breadth of virus populations during infection. Identifying the evolutionary forces acting on arboviral populations within their arthropod vector provides novel insights into arbovirus evolution.

During infection of their arthropod vectors, arthropod-borne viruses (arboviruses) such as dengue viruses traverse several anatomical barriers that are believed to cause dramatic reductions in population size. Such population bottlenecks challenge the maintenance of viral genetic diversity, which is considered critical for fitness and adaptability of arboviruses. Anatomical barriers in the vector were previously associated with both maintenance of arboviral genetic diversity and alteration of the variant repertoire. However, the relative role of random processes and natural selection, and the influence of vector genetic heterogeneity have not been elucidated. In this study, we used high-throughput sequencing to monitor dengue virus genetic diversity during infection of several genetic backgrounds of their mosquito vector. Our results show that initial infection of the vector is randomly founded by only a few tens of individual virus genomes. The overall level of viral genetic diversity generated during infection was predominantly under purifying selection but differed significantly between mosquito genetic backgrounds. Thus, in addition to random evolutionary forces and the purging of deleterious mutations that shape dengue virus genetic diversity during vector infection, our results also point to a novel role for vector genetic factors in the genetic breadth of virus populations.

Viral quasispecies complexity measures

Mutant spectrum dynamics (changes in the related mutants that compose viral populations) has a decisive impact on virus behavior. The several platforms of next generation sequencing (NGS) to study viral quasispecies offer a magnifying glass to study viral quasispecies complexity. Several parameters are available to quantify the complexity of mutant spectra, but they have limitations. Here we critically evaluate the information provided by several population diversity indices, and we propose the introduction of some new ones used in ecology. In particular we make a distinction between incidence, abundance and function measures of viral quasispecies composition. We suggest a multidimensional approach (complementary information contributed by adequately chosen indices), propose some guidelines, and illustrate the use of indices with a simple example. We apply the indices to three clinical samples of hepatitis C virus that display different population heterogeneity. Areas of virus biology in which population complexity plays a role are discussed.

Evidence of an Exponential Decay Pattern of the Hepatitis Delta Virus Evolution Rate and Fluctuations in Quasispecies Complexity in Long-Term Studies of Chronic Delta Infection

Chronic HDV infection can cause a severe form of viral hepatitis for which there is no specific treatment. Characterization of the hepatitis B or C viral quasispecies has provided insight into treatment failure and disease recurrence following liver transplantation, has proven useful to understand hepatitis B e antigen seroconversion, and has helped to predict whether hepatitis C infection will resolve or become chronic. It is likely that characterization of the hepatitis delta virus (HDV) quasispecies will ultimately have similar value for the management of this infection. This study sought to determine the RNA evolution rates in serum of chronic hepatitis delta (CHD) treatment-naïve patients, using next-generation sequencing methods. The region selected for study encompassed nucleotide positions 910 to 1270 of the genome and included the amber/W codon. Amber/W is a substrate of the editing process by the ADAR1 host enzyme and is essential for encoding the 2 delta antigens (HDAg). The amber codon encodes the small (unedited) HDAg form and the W codon the large (edited) HDAg form. The evolution rate was analyzed taking into account the time elapsed between samples, the percentage of unedited and edited genomes, and the complexity of the viral population. The longitudinal studies included 29 sequential samples from CHD patients followed up for a mean of 11.5 years. In total, 121,116 sequences were analyzed. The HDV evolution rate ranged from 9.5x10-3 to 1.2x10-3 substitutions/site/year and showed a negative correlation with the time elapsed between samples (p<0.05). An accumulation of transition-type changes was found to be responsible for higher evolution rates. The percentages of unedited and edited genomes and the quasispecies complexity showed no relationships with the evolution rate, but the fluctuations in the percentages of genomes and in complexity suggest continuous adaptation of HDV to the host conditions.

Hepatitis B virus quasispecies evolution after liver transplantation in patients under long-term lamivudine prophylaxis with or without hepatitis B immune globulin

AIMS

To investigate an optimal long-term prophylactic strategy for prevention of hepatitis B virus (HBV) recurrence after liver transplantation, we conducted a randomized study of 29 transplant recipients receiving a short course of hepatitis B immune globulin (HBIg) + lamivudine (LAM), followed by randomization to long-term prophylaxis with LAM with or without HBIg.

METHODS

The efficacy and safety, and impact on survival and HBV recurrence of these 2 prophylactic regimens were compared over a mean period of 10 years. In patients with viral recurrence, the HBV quasispecies in the surface/polymerase region were studied by ultra-deep pyrosequencing (UDPS).

RESULTS

The 10-year survival rate was 76% and was not affected by the type of prophylaxis. Four patients had hepatitis B surface antigen (HBsAg) recurrence within the first 48 months after orthotopic liver transplantation (OLT). HBsAg-positive and -negative patients showed similar mean survival times, with no differences between the 2 regimens. Low HBV DNA levels were transiently detected in 32% of HBsAg-negative patients. UDPS showed major changes after OLT in the HBV quasispecies of patients with viral recurrence, which may be explained by a "bottleneck" effect of OLT together with prophylactic therapy.

CONCLUSION

Long-term survival after OLT in end-stage chronic hepatitis B patients was good with both prophylactic strategies. However, low, transient HBV DNA levels were detected even in the absence of HBsAg, showing the importance of continuing HBV prophylaxis.

Quasispecies dynamics in hepatitis C liver transplant recipients receiving grafts from hepatitis C virus infected donors

The allocation of liver grafts from hepatitis C virus (HCV)-positive donors in HCV-infected liver transplant (LT) recipients leads to infection with two different viral populations. In a previous study, we examined quasispecies dynamics during reinfection by clonal sequencing, which did not allow an accurate characterization of coexistence and competition events. To overcome this limitation, here we used deep-sequencing analysis of a fragment of the HCV NS5B gene in six HCV-infected LT recipients who received HCV-infected grafts. Successive expansions and contractions of quasispecies complexity were observed, evolving in all cases towards a more homogeneous population. The population that became dominant was the one displaying the highest mutant spectrum complexity. In four patients, coexistence of minority mutants, derived from the donor or the recipient, were detected. In conclusion, our study shows that, during reinfection with a different HCV strain in LT recipients, the viral population with the highest diversity always becomes dominant.

Differentiation of acute from chronic hepatitis C virus infection by nonstructural 5B deep sequencing: a population-level tool for incidence estimation

The ability to classify acute versus chronic hepatitis C virus (HCV) infections at the time of diagnosis is desirable to improve the quality of surveillance information. The aim of this study was to differentiate acute from chronic HCV infections utilizing deep sequencing. HCV nonstructural 5B (NS5B) amplicons (n = 94) were generated from 77 individuals (13 acute and 64 chronic HCV infections) in British Columbia, Canada, with documented seroconversion time frames. Amplicons were deep sequenced and HCV genomic diversity was measured by Shannon entropy (SE) and a single nucleotide variant (SNV) analysis. The relationship between each diversity measure and the estimated days since infection was assessed using linear mixed models, and the ability of each diversity measure to differentiate acute from chronic infections was assessed using generalized estimating equations. Both SE and the SNV diversity measures were significantly different for acute versus chronic infections (P < 0.009). NS5B nucleotide diversity continued to increase for at least 3 years postinfection. Among individuals with the least uncertainty with regard to duration of infection (n = 39), the area under the receiver operating characteristic curve (AUROC) was high (0.96 for SE; 0.98 for SNV). Although the AUROCs were lower (0.86 for SE; 0.80 for SNV) when data for all individuals were included, they remain sufficiently high for epidemiological purposes. Synonymous mutations were the primary discriminatory variable accounting for over 78% of the measured genetic diversity.

CONCLUSIONS

NS5B sequence diversity assessed by deep sequencing can differentiate acute from chronic HCV infections and, with further validation, could become a powerful population-level surveillance tool for incidence estimation.

Experimental evolution of an RNA virus in wild birds: evidence for host-dependent impacts on population structure and competitive fitness

Within hosts, RNA viruses form populations that are genetically and phenotypically complex. Heterogeneity in RNA virus genomes arises due to error-prone replication and is reduced by stochastic and selective mechanisms that are incompletely understood. Defining how natural selection shapes RNA virus populations is critical because it can inform treatment paradigms and enhance control efforts. We allowed West Nile virus (WNV) to replicate in wild-caught American crows, house sparrows and American robins to assess how natural selection shapes RNA virus populations in ecologically relevant hosts that differ in susceptibility to virus-induced mortality. After five sequential passages in each bird species, we examined the phenotype and population diversity of WNV through fitness competition assays and next generation sequencing. We demonstrate that fitness gains occur in a species-specific manner, with the greatest replicative fitness gains in robin-passaged WNV and the least in WNV passaged in crows. Sequencing data revealed that intrahost WNV populations were strongly influenced by purifying selection and the overall complexity of the viral populations was similar among passaged hosts. However, the selective pressures that control WNV populations seem to be bird species-dependent. Specifically, crow-passaged WNV populations contained the most unique mutations (~1.7× more than sparrows, ~3.4× more than robins) and defective genomes (~1.4× greater than sparrows, ~2.7× greater than robins), but the lowest average mutation frequency (about equal to sparrows, ~2.6× lower than robins). Therefore, our data suggest that WNV replication in the most disease-susceptible bird species is positively associated with virus mutational tolerance, likely via complementation, and negatively associated with the strength of selection. These differences in genetic composition most likely have distinct phenotypic consequences for the virus populations. Taken together, these results reveal important insights into how different hosts may contribute to the emergence of RNA viruses.

Deep sequencing and phylogenetic analysis of variants resistant to interferon-based protease inhibitor therapy in chronic hepatitis induced by genotype 1b hepatitis C virus

Because of recent advances in deep sequencing technology, detailed analysis of hepatitis C virus (HCV) quasispecies and their dynamic changes in response to direct antiviral agents (DAAs) became possible, although the role of quasispecies is not fully understood. In this study, to clarify the evolution of viral quasispecies and the origin of drug-resistant mutations induced by interferon (IFN)-based protease inhibitor therapy, the nonstructural-3 (NS3) region of genotype 1b HCV in 34 chronic hepatitis patients treated with telaprevir (TVR)/pegylated interferon (PEG-IFN)/ribavirin (RBV) was subjected to a deep sequencing study coupled with phylogenetic analysis. Twenty-six patients (76.5%) achieved a sustained viral response (SVR), while 8 patients did not (non-SVR; 23.5%). When the complexity of the quasispecies was expressed as the mutation frequency or Shannon entropy value, a significant decrease in the IFNL3 (rs8099917) TT group and a marginal decrease in the SVR group were found soon (12 h) after the introduction of treatment, whereas there was no decrease in the non-SVR group and no significant decrease in mutation frequency in the IFNL3 TG/GG group. In the analysis of viral quasispecies composition in non-SVR patients, major populations greatly changed, accompanied by the appearance of resistance, and the compositions were unlikely to return to the pretreatment composition even after the end of therapy. Clinically TVR-resistant variants were observed in 5 non-SVR patients (5/8, 62.5%), all of which were suspected to have acquired resistance by mutations through phylogenetic analysis. In conclusion, results of the study have important implications for treatment response and outcome in interferon-based protease inhibitor therapy.

IMPORTANCE

In the host, hepatitis C virus (HCV) consists of a variety of populations (quasispecies), and it is supposed that dynamic changes in quasispecies are closely related to pathogenesis, although this is poorly understood. In this study, recently developed deep sequencing technology was introduced, and changes in quasispecies associated with telaprevir (TVR)/pegylated interferon (PEG-IFN)/ribavirin (RBV) triple therapy and their clinical significance were investigated extensively by phylogenetic tree analysis. Through this study, the associations among treatment response, changes in viral quasispecies complexity in the early stage of treatment, changes in the quasispecies composition, and origin of TVR-resistant variant HCV were elucidated.

Advanced molecular surveillance of hepatitis C virus

Hepatitis C virus (HCV) infection is an important public health problem worldwide. HCV exploits complex molecular mechanisms, which result in a high degree of intrahost genetic heterogeneity. This high degree of variability represents a challenge for the accurate establishment of genetic relatedness between cases and complicates the identification of sources of infection. Tracking HCV infections is crucial for the elucidation of routes of transmission in a variety of settings. Therefore, implementation of HCV advanced molecular surveillance (AMS) is essential for disease control. Accounting for virulence is also important for HCV AMS and both viral and host factors contribute to the disease outcome. Therefore, HCV AMS requires the incorporation of host factors as an integral component of the algorithms used to monitor disease occurrence. Importantly, implementation of comprehensive global databases and data mining are also needed for the proper study of the mechanisms responsible for HCV transmission. Here, we review molecular aspects associated with HCV transmission, as well as the most recent technological advances used for virus and host characterization. Additionally, the cornerstone discoveries that have defined the pathway for viral characterization are presented and the importance of implementing advanced HCV molecular surveillance is highlighted.

High-resolution hepatitis C virus subtyping using NS5B deep sequencing and phylogeny, an alternative to current methods

Hepatitis C virus (HCV) is classified into seven major genotypes and 67 subtypes. Recent studies have shown that in HCV genotype 1-infected patients, response rates to regimens containing direct-acting antivirals (DAAs) are subtype dependent. Currently available genotyping methods have limited subtyping accuracy. We have evaluated the performance of a deep-sequencing-based HCV subtyping assay, developed for the 454/GS-Junior platform, in comparison with those of two commercial assays (Versant HCV genotype 2.0 and Abbott Real-time HCV Genotype II) and using direct NS5B sequencing as a gold standard (direct sequencing), in 114 clinical specimens previously tested by first-generation hybridization assay (82 genotype 1 and 32 with uninterpretable results). Phylogenetic analysis of deep-sequencing reads matched subtype 1 calling by population Sanger sequencing (69% 1b, 31% 1a) in 81 specimens and identified a mixed-subtype infection (1b/3a/1a) in one sample. Similarly, among the 32 previously indeterminate specimens, identical genotype and subtype results were obtained by direct and deep sequencing in all but four samples with dual infection. In contrast, both Versant HCV Genotype 2.0 and Abbott Real-time HCV Genotype II failed subtype 1 calling in 13 (16%) samples each and were unable to identify the HCV genotype and/or subtype in more than half of the non-genotype 1 samples. We concluded that deep sequencing is more efficient for HCV subtyping than currently available methods and allows qualitative identification of mixed infections and may be more helpful with respect to informing treatment strategies with new DAA-containing regimens across all HCV subtypes.

Error correction and diversity analysis of population mixtures determined by NGS

The impetus for this work was the need to analyse nucleotide diversity in a viral mix taken from honeybees. The paper has two findings. First, a method for correction of next generation sequencing error in the distribution of nucleotides at a site is developed. Second, a package of methods for assessment of nucleotide diversity is assembled. The error correction method is statistically based and works at the level of the nucleotide distribution rather than the level of individual nucleotides. The method relies on an error model and a sample of known viral genotypes that is used for model calibration. A compendium of existing and new diversity analysis tools is also presented, allowing hypotheses about diversity and mean diversity to be tested and associated confidence intervals to be calculated. The methods are illustrated using honeybee viral samples. Software in both Excel and Matlab and a guide are available at http://www2.warwick.ac.uk/fac/sci/systemsbiology/research/software/, the Warwick University Systems Biology Centre software download site.

Deep sequencing: becoming a critical tool in clinical virology

Population (Sanger) sequencing has been the standard method in basic and clinical DNA sequencing for almost 40 years; however, next-generation (deep) sequencing methodologies are now revolutionizing the field of genomics, and clinical virology is no exception. Deep sequencing is highly efficient, producing an enormous amount of information at low cost in a relatively short period of time. High-throughput sequencing techniques have enabled significant contributions to multiples areas in virology, including virus discovery and metagenomics (viromes), molecular epidemiology, pathogenesis, and studies of how viruses to escape the host immune system and antiviral pressures. In addition, new and more affordable deep sequencing-based assays are now being implemented in clinical laboratories. Here, we review the use of the current deep sequencing platforms in virology, focusing on three of the most studied viruses: human immunodeficiency virus (HIV), hepatitis C virus (HCV), and influenza virus.

Clinical application of estimating hepatitis B virus quasispecies complexity by massive sequencing: correlation between natural evolution and on-treatment evolution

AIM

To evaluate HBV quasispecies (QA) complexity in the preCore/Core regions in relation to HBeAg status, and explore QA changes under natural evolution and nucleoside analogue (NUC) treatment.

METHODS

Ultra-deep pyrosequencing of HBV preCore/Core regions in 30 sequential samples (baseline [diagnosis], treatment-free, and treatment-nonresponse) from 10 retrospectively selected patients grouped according to HBeAg status over time: HBeAg+ (N = 4), HBeAg- (N = 2), and fluctuating HBeAg (transient seroreversion/seroconversion pattern) (N = 4). QA complexity was defined by Shannon entropy, mutation frequency, nucleotide diversity, and mutation frequency of amino acids (MfAA) in preCore and Core.

RESULTS

The QA was less complex in HBeAg+ than in HBeAg- or fluctuating HBeAg. High complexity in preCore was associated with decreased viral replication (preCore MfAA negatively correlated with HBV-DNA, p = 0.005). QA complexity in the treatment-free period negatively correlated with values seen during treatment. Specific variants were mainly selected in the Core region in HBeAg- and fluctuating HBeAg patients, suggesting higher immune pressure than in HBeAg+.

CONCLUSIONS

The negative correlation between QA natural evolution and on-treatment evolution indicates the importance of pre-treatment QA study to predict QA changes in NUC nonresponders. Study of QA complexity could be useful for managing HBV infection.