RT-PCR amplification and cloning of large viral sequences

Within-host quantitation of anellovirus genome complexity from clinical samples

Anellovirus (AV) is a ubiquitous and diverse virus in the human population. An individual can be infected with multiple AV genera and species that form a heterogeneous repertoire, called the anellome. Due to its exceptional genetic diversity, efficient evaluation of anellome complexity remains a methodological challenge. In the current study, AV genome was first enriched from patient serum samples through two-phase rolling circle amplification. Following Illumina sequencing, anellome was analyzed with an advanced bioinformatics pipeline, including read extraction at three similarity levels, de novo assembly, species assignment, and determination of relative abundance among AV variants. The method was validated in the mock sample and then applied to 21 hepatitis C virus (HCV) patients with and without hepatocellular carcinoma (HCC). Overall, there was a large variance regarding AV richness, ranging from 2 to 51 AV species. In contrast to HCV patients without HCC, HCC incidence was associated with reduced richness (12.6 ± 14.4 vs. 35.4 ± 13.6, p = 0.001) and Shannon entropy (0.4 ± 0.34 vs. 0.61 ± 0.12, p = 0.095) at the AV species level. Interestingly, AV genus beta and gamma expanded in the anellome in 7 of 10 HCC patients. These observations shed light on the potential association between anellome and HCC incidence in patients with chronic HCV infection. The method presented here represents a valuable tool to investigate the role of anellome in human health and disease.

Molecular basis of genetic variation of viruses

Genetic variation is a necessity of all biological systems. Viruses use all known mechanisms of variation; mutation, several forms of recombination, and segment reassortment in the case of viruses with a segmented genome. These processes are intimately connected with the replicative machineries of viruses, as well as with fundamental physical-chemical properties of nucleotides when acting as template or substrate residues. Recombination has been viewed as a means to rescue viable genomes from unfit parents or to produce large modifications for the exploration of phenotypic novelty. All types of genetic variation can act conjointly as blind processes to provide the raw materials for adaptation to the changing environments in which viruses must replicate. A distinction is made between mechanistically unavoidable and evolutionarily relevant mutation and recombination.

Interpatient mutational spectrum of human coronavirus-OC43 revealed by illumina sequencing

Human coronaviruses (HCoV) are RNA viruses that cause respiratory tract infections with viral replication of limited duration. The host and viral population heterogeneity could influence clinical phenotypes. Employing long RT‐PCR with Illumina sequencing, we quantified the gene mutation load at 0.5% mutation frequency for the 4529 bp‐domain spanning the Spike gene (4086 bp) of HCoV‐OC43 in four upper respiratory clinical specimens obtained during acute illness. There were a total of 121 mutations for all four HCoV samples with the average number of mutations at 30.3 ± 10.2, which is significantly higher than that expected from the Illumina sequencing error rate. There were two mutation peaks, one at the 5′ end and the other near position 1 550 in the S1 subunit. Two coronavirus samples were genotype B and two were genotype D, clustering with HCoV‐OC43 strain AY391777 in neighbor‐joining tree phylogenetic analysis. Nonsynonymous mutations were 76.1 ± 14% of mutation load. Although lower than other RNA viruses such as hepatitis C virus, HCoV‐OC43 did exhibit quasi‐species. The rate of nonsynonymous mutations was higher in the HCoV‐OC43 isolates than in hepatitis C (HCV) virus genotype 1a isolates analyzed for comparison in this study. These characteristics of HCoV‐OC43 may affect viral replication dynamics, receptor binding, antigenicity, evolution, transmission, and clinical illness.

HCV inter-subtype 1a/1b recombinant detected by complete-genome next-generation sequencing

Next-generation sequencing (NGS) provides a practical approach to HCV complete-genome sequencing, detecting low-frequency variants and allowing analysis of viral genetic diversity (quasispecies) in the sample, and so far, it is very useful for identifying preexisting drug-resistant mutants and emerging escape mutations, as well as detecting viral recombinants containing genomic regions from different genotypes and subtypes. The aim of this study was to analyze the complete coding region of hepatitis C virus (HCV) genotype 1 (subtypes 1a and 1b) from patients with chronic infection who were direct-acting antiviral (DAA) naïve. Next-generation sequencing (Ion Torrent™ PGM) was used to determine the sequence of the complete coding region of 100 HCV-monoinfected DAA-naïve patients (51 and 49 subtypes 1a and 1b, respectively). We report the first description of nearly complete HCV genome sequences of subtype 1a and 1b isolates from a large population of Brazilian patients with chronic hepatitis C, and HCV-1a grouped in two different clades. Using this methodology, an inter-subtype 1a/1b recombinant was identified in this study.

Molecular Basis of Genetic Variation of Viruses

Genetic variation is a necessity of all biological systems. Viruses use all known mechanisms of variation: mutation, several forms of recombination, and segment reassortment in the case of viruses with a segmented genome. These processes are intimately connected with the replicative machineries of viruses, as well as with fundamental physico-chemical properties of nucleotides when acting as template or substrate residues. Recombination has been viewed as a means to rescue viable genomes from unfit parents, or to produce large modifications for the exploration of phenotypic novelty. All types of genetic variation can act conjointly as blind processes to provide the raw materials for adaptation to the changing environments in which viruses must replicate.

Evidence for deleterious hepatitis C virus quasispecies mutation loads that differentiate the response patterns in IFN-based antiviral therapy

Viral quasispecies (QS) have long been considered to affect the efficiency of hepatitis C virus (HCV) antiviral therapy, but a correlation between QS diversity and treatment outcomes has not been established conclusively. We previously measured HCV QS diversity by genome-wide quantification of high-resolution mutation load in HCV genotype 1a patients achieving a sustained virological response (1a/SVR) or a null response (1a/null). The current study extended this work into HCV 1a patients experiencing relapse (1a/relapse, n = 19) and genotype 2b patients with SVR (2b/SVR, n = 10). The mean mutation load per patient in 2b/SVR and 1a/relapse was similar, respectively, to 1a/SVR (517.6 ± 174.3 vs 524 ± 278.8 mutations, P = 0.95) and 1a/null (829.2 ± 282.8 vs 805.6 ± 270.7 mutations, P = 0.78). Notably, a deleterious mutation load, as indicated by the percentage of non-synonymous mutations, was highest in 2b/SVR (33.2 ± 8.5%) as compared with 1a/SVR (23.6 ± 7.8%, P = 0.002), 1a/null (18.2 ± 5.1%, P = 1.9 × 10(-7)) or 1a/relapse (17.8 ± 5.3%, P = 1.8) × 10(-6). In the 1a/relapse group, continuous virus evolution was observed with excessive accumulation of a deleterious load (17.8 ± 5.3% vs 35.4 ± 12.9%, P = 3.5 × 10(-6)), supporting the functionality of Muller's ratchet in a treatment-induced population bottleneck. Taken together, the magnitude of HCV mutation load, particularly the deleterious mutation load, provides an evolutionary explanation for the emergence of multiple response patterns as well as an overall high SVR rate in HCV genotype 2 patients. Augmentation of Muller's ratchet represents a potential strategy to reduce or even eliminate viral relapse in HCV antiviral therapy.

High-resolution quantification of hepatitis C virus genome-wide mutation load and its correlation with the outcome of peginterferon-alpha2a and ribavirin combination therapy

Hepatitis C virus (HCV) is a highly mutable RNA virus and circulates as a heterogeneous population in individual patients. The magnitude of such population heterogeneity has long been proposed to be linked with diverse clinical phenotypes, including antiviral therapy. Yet data accumulated thus far are fairly inconclusive. By the integration of long RT-PCR with 454 sequencing, we have built a pipeline optimized for the quantification of HCV genome-wide mutation load at 1% resolution of mutation frequency, followed by a retrospective study to examine the role of HCV mutation load in peginterferon-alpha2a and ribavirin combination antiviral therapy. Genome-wide HCV mutation load varied widely with a range from 92 to 1639 mutations and presented a Poisson distribution among 56 patients (Kolmogorov-Smirnov statistic  = 0.078, p = 0.25). Patients achieving sustained virological response (n = 26) had significantly lower mutation loads than that in null responders (n = 30) (mean and standard derivation: 524±279 vs. 805±271, p = 0.00035). All 36,818 mutations detected in 56 patients displayed a power-law distribution in terms of mutation frequency in viral population. The low-frequency mutation load, but not the high-frequency load, was proportional firmly to the total mutation load. In-depth analyses revealed that intra-patient HCV population structure was shaped by multiple factors, including immune pressure, strain difference and genetic drift. These findings explain previous conflicting reports using low-resolution methods and highlight a dominant role of natural selection in response to therapeutic intervention. By attaining its signatures from complex interaction between host and virus, the high-resolution quantification of HCV mutation load predicts outcomes from interferon-based antiviral therapy and could also be a potential biomarker in other clinical settings.