Ultra-deep pyrosequencing (UDPS) data treatment to study amplicon HCV minor variants

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.

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.

Ultra Deep Sequencing of Circulating Cell-Free DNA as a Potential Tool for Hepatocellular Carcinoma Management

Simple Summary

In this unicentric prospective study, we analyzed the most prevalent mutations in HCC (TERT promoter, TP53, CTNNB1, AXIN1 and ARID1A) in plasma cfDNA by next-generation sequencing, aiming to elucidate their value as prognostic noninvasive biomarkers. Total cfDNA (cut-off value 2 ng/µL), number of mutated genes and number of detectable mutations on cfDNA were significantly related to mortality. Number of mutated genes and number of detected mutations in cfDNA and the ratio between number of mutations and total amount of cfDNA were also significantly associated with recurrence. Detection of more than four mutations in cfDNA correlated with a higher risk of death. Dynamic changes in cfDNA mutations were detected prior to radiological detection of HCC recurrence. We believe that these results support the proof of principle and launching of validation studies to confirm that total cfDNA and detection of prevalent HCC mutations could have prognostic implications in early-stage HCC patients.

Abstract

Background: Cell-free DNA (cfDNA) concentrations have been described to be inversely correlated with prognosis in cancer. Mutations in HCC-associated driver genes in cfDNA have been reported, but their relation with patient’s outcome has not been described. Our aim was to elucidate whether mutations found in cfDNA could be representative from those present in HCC tissue, providing the rationale to use the cfDNA to monitor HCC. Methods: Tumoral tissue, paired nontumor adjacent tissue and blood samples were collected from 30 HCC patients undergoing curative therapies. Deep sequencing targeting HCC driver genes was performed. Results: Patients with more than 2 ng/µL of cfDNA at diagnosis had higher mortality (mean OS 24.6 vs. 31.87 months, p = 0.01) (AUC = 0.782). Subjects who died during follow-up, had a significantly higher number of mutated genes (p = 0.015) and number of mutations (p = 0.015) on cfDNA. Number of mutated genes (p = 0.001), detected mutations (p = 0.001) in cfDNA and ratio (number of mutations/cfDNA) (p = 0.003) were significantly associated with recurrence. However, patients with a ratio (number of mutations/cfDNA) above 6 (long-rank p = 0.0003) presented a higher risk of recurrence than those with a ratio under 6. Detection of more than four mutations in cfDNA correlated with higher risk of death (long-rank p = 0.042). Conclusions: In summary, cfDNA and detection of prevalent HCC mutations could have prognostic implications in early-stage HCC patients

[Impact of NGS results on patient outcome with a multiform glioblastoma]

PURPOSE

Although some genetic alterations in glioblastoma (GBM) have been characterized, the prognostic value of these gene mutations is not yet established in patients treated with standard therapy.

PATIENTS AND METHOD

40 patients with newly diagnosed GBM, treated between July 2017 and December 2019, and who had genomic analysis were analyzed. Next-generation sequencing techniques (NGS) were used with a panel of 26 genes. Patients were grouped according to MGMT status, the presence or absence of at least one mutated gene on the panel, and p53 expression by immunohistochemistry.

RESULTS

the median follow-up was 11.5 months (1.0-37). For all patients, the median duration of progression-free survival was 8 months (95% CI, 5.3-10.7) and the median overall survival (OS) was 17 months (95% CI, 7.5-26.5). Progression-free and overall survival were significantly different according to MGMT status but not according to NGS and p53 status. Three groups of patients according to different combined status could be distinguished due to significant differences in overall survival.

CONCLUSION

we have shown that the presence of MGMT promoter methylation is a good prognostic factor. By grouping the patients according to their MGMT, NGS and p53 status, three groups of patients could be separated according to their overall survival. However, these results must be confirmed on a larger number of patients.

Hepatitis B Virus Variants with Multiple Insertions and/or Deletions in the X Open Reading Frame 3′ End: Common Members of Viral Quasispecies in Chronic Hepatitis B Patients

Deletions in the 3′ end region of the hepatitis B virus (HBV) X open reading frame (HBX) may affect the core promoter (Cp) and have been frequently associated with hepatocellular carcinoma (HCC). The aim of this study was to investigate the presence of variants with deletions and/or insertions (Indels) in this region in the quasispecies of 50 chronic hepatitis B (CHB) patients without HCC. We identified 103 different Indels in 47 (94%) patients, in a median of 3.4% of their reads (IQR, 1.3–8.4%), and 25% (IQR, 13.1–40.7%) of unique sequences identified in each quasispecies (haplotypes). Of those Indels, 101 (98.1%) caused 44 different altered stop codons, the most commonly observed were at positions 128, 129, 135, and 362 (putative position). Moreover, 39 (37.9%) Indels altered the TATA-like box (TA) sequences of Cp; the most commonly observed caused TA2 + TA3 fusion, creating a new putative canonical TATA box. Four (8%) patients developed negative clinical outcomes after a median follow-up of 9.4 (8.7–12) years. In conclusion, we observed variants with Indels in the HBX 3′ end in the vast majority of our CHB patients, some of them encoding alternative versions of HBx with potential functional roles, and/or alterations in the regulation of transcription.

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.

Study of Quasispecies Complexity and Liver Damage Progression after Liver Transplantation in Hepatitis C Virus Infected Patients

Cirrhosis derived from chronic hepatitis C virus (HCV) infection is still a common indication for liver transplantation (LT). Reinfection of the engrafted liver is universal in patients with detectable viral RNA at the time of transplant and causes fast progression of cirrhosis (within 5 years) in around one-third of these patients. To prevent damage to the liver graft, effective direct-acting antiviral (DAA) therapy is required as soon as possible. However, because of post-LT clinical instability, it is difficult to determine the optimal time to start DAAs with a low risk of complications. Evaluate changes in quasispecies complexity following LT and seek a predictive index of fast liver damage progression to determine the timing of DAA initiation. HCV genomes isolated from pre-LT and 15-day post-LT serum samples of ten patients, who underwent orthotopic LT, were quantified and sequenced using a next-generation sequencing platform. Sequence alignments, phylogenetic trees, quasispecies complexity measures, biostatistics analyses, adjusted R2 values, and analysis of variance (ANOVA) were carried out. Three different patterns of reinfection were observed (viral bottlenecking, conserved pre-LT population, and mixed populations), suggesting that bottlenecking or homogenization of the viral population is not a generalized effect after liver graft reinfection. None of the quasispecies complexity measures predicted the future degree of liver damage. Higher and more uniform viral load (VL) values were observed in all pre-LT samples, but values were more dispersed in post-LT samples. However, VL increased significantly from the pre-LT to 15-day post-LT samples in patients with advanced fibrosis at 1-year post-LT, suggesting that a VL increase on day 15 may be a predictor of fast liver fibrosis progression. HCV kinetics after LT differ between patients and are not fibrosis-dependent. Higher VL at day 15 post-LT versus pre-LT samples may predict fast liver fibrosis progression.

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.

Sophisticated viral quasispecies with a genotype-related pattern of mutations in the hepatitis B X gene of HBeAg-ve chronically infected patients

Patients with HBeAg-negative chronic infection (CI) have not been extensively studied because of low viremia. The HBx protein, encoded by HBX, has a key role in viral replication. Here, we analyzed the viral quasispecies at the 5' end of HBX in CI patients and compared it with that of patients in other clinical stages. Fifty-eight HBeAg-negative patients were included: 16 CI, 19 chronic hepatitis B, 16 hepatocellular carcinoma and 6 liver cirrhosis. Quasispecies complexity and conservation were determined in the region between nucleotides 1255 and 1611. Amino acid changes detected were tested in vitro. CI patients showed higher complexity in terms of mutation frequency and nucleotide diversity and higher quasispecies conservation (p < 0.05). A genotype D-specific pattern of mutations (A12S/P33S/P46S/T36D-G) was identified in CI (median frequency, 81.7%), which determined a reduction in HBV DNA release of up to 1.5 log in vitro. CI patients showed a more complex and conserved viral quasispecies than the other groups. The genotype-specific pattern of mutations could partially explain the low viremia observed in these patients.

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.

Deep-sequencing reveals broad subtype-specific HCV resistance mutations associated with treatment failure

A percentage of hepatitis C virus (HCV)-infected patients fail direct acting antiviral (DAA)-based treatment regimens, often because of drug resistance-associated substitutions (RAS). The aim of this study was to characterize the resistance profile of a large cohort of patients failing DAA-based treatments, and investigate the relationship between HCV subtype and failure, as an aid to optimizing management of these patients. A new, standardized HCV-RAS testing protocol based on deep sequencing was designed and applied to 220 previously subtyped samples from patients failing DAA treatment, collected in 39 Spanish hospitals. The majority had received DAA-based interferon (IFN) α-free regimens; 79% had failed sofosbuvir-containing therapy. Genomic regions encoding the nonstructural protein (NS) 3, NS5A, and NS5B (DAA target regions) were analyzed using subtype-specific primers. Viral subtype distribution was as follows: genotype (G) 1, 62.7%; G3a, 21.4%; G4d, 12.3%; G2, 1.8%; and mixed infections 1.8%. Overall, 88.6% of patients carried at least 1 RAS, and 19% carried RAS at frequencies below 20% in the mutant spectrum. There were no differences in RAS selection between treatments with and without ribavirin. Regardless of the treatment received, each HCV subtype showed specific types of RAS. Of note, no RAS were detected in the target proteins of 18.6% of patients failing treatment, and 30.4% of patients had RAS in proteins that were not targets of the inhibitors they received. HCV patients failing DAA therapy showed a high diversity of RAS. Ribavirin use did not influence the type or number of RAS at failure. The subtype-specific pattern of RAS emergence underscores the importance of accurate HCV subtyping. The frequency of "extra-target" RAS suggests the need for RAS screening in all three DAA target regions.

Synergistic lethal mutagenesis of hepatitis C virus

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 mutagenic agent, often a nucleotide analogue. One of its advantages is its broad spectrum nature that renders the strategy potentially effective against emergent RNA viral infections. Here we describe synergistic lethal mutagenesis of hepatitis C virus (HCV) by a combination of favipiravir (T-705) and ribavirin. Synergy has been documented over a broad range of analogue concentrations using the Chou-Talalay method as implemented in the CompuSyn graphics, with average dose reduction index (DRI) above 1 (68.02±101.6 for favipiravir, and 5.83±6.07 for ribavirin), and average combination indices (CI) below 1 (0.52±0.28). Furthermore, analogue concentrations that individually did not extinguish high fitness HCV in ten serial infections, when used in combination they extinguished high fitness HCV in one to two passages. Although both analogues display a preference for G→A and C→U transitions, deep sequencing analysis of mutant spectra indicated a different preference of the two analogues for the mutation sites, thus unveiling a new possible synergy mechanism in lethal mutagenesis. Prospects of synergy among mutagenic nucleotides as a strategy to confront emerging viral infections are discussed.

Hepatitis C virus early kinetics and resistance-associated substitution dynamics during antiviral therapy with direct-acting antivirals

The emergence of resistance-associated substitutions (RASs) can compromise the high efficacy of direct-acting antivirals (DAAs). Little is known about RASs selection at very early time points during DAA treatment. Therefore, we analyzed the potential emergence of RASs immediately after therapy initiation. Samples of 71 patients treated with different DAAs were collected at baseline, during therapy (hours 4 and 8; days 1-7; weeks 2-4) or until target not detected. HCV-RNA levels were determined by qPCR, and RASs were detected by deep sequencing. Sixty-three (89%) patients achieved a sustained virological response (SVR), 7 (10%) relapsed, and 1 (1%) experienced a breakthrough. Almost all non-SVR (7/8, 88%) showed RASs either at baseline or relapse. High-frequency RASs detected at baseline (Y93H and L159F+C316N) remained detectable at early time points during therapy and reappeared as most prevalent substitutions at relapse. Conversely, emergent RASs at relapse (Q80R, D168E/V, R155K and L31V) were not observed during the first hours-days, before HCV-RNA became undetectable. HCV-RNA decay and genetic evolution of the quasispecies followed a similar pattern during the first hours of therapy in SVR and non-SVR patients. In conclusion, the absence of early RASs selection and the similar dynamics of HCV kinetics and quasispecies in SVR and non-SVR patients after therapy initiation suggest that RASs selection may occur at later stages in the remaining reservoir, where viral populations persist hidden at very low replication levels. Nevertheless, we cannot completely exclude very early selection, when RASs are present below the sensitivity limit of deep sequencing.

Evidence for positive selection of hepatitis A virus antigenic variants in vaccinated men-having-sex-with men patients: Implications for immunization policies

Background

A huge outbreak in the men-having-sex-with-men (MSM) has hit Europe during the years 2016–2018. Outbreak control has been hampered by vaccine shortages in many countries, and to minimize their impact, reduction of antigen doses has been implemented. However, these measures may have consequences on the evolution of hepatitis A virus (HAV), leading to the emergence of antigenic variants. Cases in vaccinated MSM patients have been detected in Barcelona, opening the possibility to study HAV evolution under immune pressure.

Methods

We performed deep-sequencing analysis of ten overlapping fragments covering the complete capsid coding region of HAV. A total of 14578255 reads were obtained and used for the analysis of virus evolution in vaccinated versus non-vaccinated patients. We estimated maximum and minimum mutation frequencies, and Shannon entropy in the quasispecies of each patient. Non-synonymous (NSyn) mutations affecting residues exposed in the capsid surface were located, with respect to epitopes, using the recently described crystal structure of HAV, as an indication of its potential role in escaping to the effect of vaccines.

Findings

HAV evolution at the quasispecies level, in non-vaccinated and vaccinated patients, revealed higher diversity in epitope-coding regions of the vaccinated group. Although amino acid replacements occurring in and around the epitopes were observed in both groups, their abundance was significantly higher in the quasispecies of vaccinated patients, indicating ongoing processes of fixation.

Interpretation

Our data suggest positive selection of antigenic variants in some vaccinated patients, raising concerns for new vaccination polices directed to the MSM group.

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.

Using high-throughput sequencing for investigating intra-host hepatitis C evolution over long retrospective periods

Collections of biological samples held by hospitals represent invaluable resources for conducting retrospective evolutionary studies of chronic infections. Using high-throughput sequencing, those collections permit analysis of within-host genetic diversity over long follow-up periods, and allow a better understanding of resistance to treatment regimes during disease evolution. Here, we studied the evolution of hepatitis C virus (HCV) populations in two patients with an absence of response to dual therapies. We implemented amplicon sequencing to survey genomic variation at the Core and NS5B regions of HCV over a period of 13 years from blood samples obtained at multiple time points. We observed mixed infection by multiple HCV genotypes in both patients. Genetic heterogeneity and sample composition analysis provided information about the changes in viral population over the course of clinical treatment, with NS5B experiencing an increase in diversity after treatment initiation. Secondary infections were estimated to predate treatment year, and our results pointed towards diversifying selection occurring post-treatment, acting on standing genomic variation and maintaining high genetic heterogeneity during infection. For these two patients infected with multiple HCV genotypes, the maintenance of viral diversity was explained with the hypothesis of soft selective sweep started at the same time as antiviral treatment was initiated.

Quasispecies dynamics and clinical significance of hepatitis C virus (HCV) antiviral resistance

Hepatitis C virus (HCV) follows quasispecies dynamics in infected hosts and this influences its biology, how the virus diversifies into several genotypes and many subtypes, and how viral populations respond to antiviral therapies. Despite current antiviral combinations being able to cure a great percentage of HCV-infected patients, the presence of resistance-associated substitutions (RASs) diminishes the success of antiviral therapies, which is a main concern in the re-treatment of patients treated with direct-acting antiviral agents. Current methodologies such as ultra deep sequencing are ideal tools to obtain a detailed representation of the mutant spectrum composition circulating in infected patients. Such knowledge should allow optimisation of rescue treatments. A new mechanism of antiviral resistance not based on the selection of RASs but on high viral fitness is discussed.

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.

Pipeline for specific subtype amplification and drug resistance detection in hepatitis C virus

BACKGROUND

Despite the high sustained virological response rates achieved with current directly-acting antiviral agents (DAAs) against hepatitis C virus (HCV), around 5-10% of treated patients do not respond to current antiviral therapies, and basal resistance to DAAs is increasingly detected among treatment-naïve infected individuals. Identification of amino acid substitutions (including those in minority variants) associated with treatment failure requires analytical designs that take into account the high diversification of HCV in more than 86 subtypes according to the ICTV website (June 2017).

METHODS

The methodology has involved five sequential steps: (i) to design 280 oligonucleotide primers (some including a maximum of three degenerate positions), and of which 120 were tested to amplify NS3, NS5A-, and NS5B-coding regions in a subtype-specific manner, (ii) to define a reference sequence for each subtype, (iii) to perform experimental controls to define a cut-off value for detection of minority amino acids, (iv) to establish bioinformatics' tools to quantify amino acid replacements, and (v) to validate the procedure with patient samples.

RESULTS

A robust ultra-deep sequencing procedure to analyze HCV circulating in serum samples from patients infected with virus that belongs to the ten most prevalent subtypes worldwide: 1a, 1b, 2a, 2b, 2c, 2j, 3a, 4d, 4e, 4f has been developed. Oligonucleotide primers are subtype-specific. A cut-off value of 1% mutant frequency has been established for individual mutations and haplotypes.

CONCLUSION

The methodological pipeline described here is adequate to characterize in-depth mutant spectra of HCV populations, and it provides a tool to understand HCV diversification and treatment failures. The pipeline can be periodically extended in the event of HCV diversification into new genotypes or subtypes, and provides a framework applicable to other RNA viral pathogens, with potential to couple detection of drug-resistant mutations with treatment planning.

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.

HIV-1 Protease Evolvability Is Affected by Synonymous Nucleotide Recoding

One unexplored aspect of HIV-1 genetic architecture is how codon choice influences population diversity and evolvability. Here we compared the levels of development of HIV-1 resistance to protease inhibitors (PIs) between wild-type (WT) virus and a synthetic virus (MAX) carrying a codon-pair-reengineered protease sequence including 38 (13%) synonymous mutations. The WT and MAX viruses showed indistinguishable replication in MT-4 cells or peripheral blood mononuclear cells (PBMCs). Both viruses were subjected to serial passages in MT-4 cells, with selective pressure from the PIs atazanavir (ATV) and darunavir (DRV). After 32 successive passages, both the WT and MAX viruses developed phenotypic resistance to PIs (50% inhibitory concentrations [IC₅₀s] of 14.6 ± 5.3 and 21.2 ± 9 nM, respectively, for ATV and 5.9 ± 1.0 and 9.3 ± 1.9, respectively, for DRV). Ultradeep sequence clonal analysis revealed that both viruses harbored previously described mutations conferring resistance to ATV and DRV. However, the WT and MAX virus proteases showed different resistance variant repertoires, with the G16E and V77I substitutions observed only in the WT and the L33F, S37P, G48L, Q58E/K, and L89I substitutions detected only in the MAX virus. Remarkably, the G48L and L89I substitutions are rarely found in vivo in PI-treated patients. The MAX virus showed significantly higher nucleotide and amino acid diversity of the propagated viruses with and without PIs (P < 0.0001), suggesting a higher selective pressure for change in this recoded virus. Our results indicate that the HIV-1 protease position in sequence space delineates the evolution of its mutant spectrum. Nevertheless, the investigated synonymously recoded variant showed mutational robustness and evolvability similar to those of the WT virus.IMPORTANCE Large-scale synonymous recoding of virus genomes is a new tool for exploring various aspects of virus biology. Synonymous virus genome recoding can be used to investigate how a virus's position in sequence space defines its mutant spectrum, evolutionary trajectory, and pathogenesis. In this study, we evaluated how synonymous recoding of the human immunodeficiency virus type 1 (HIV-1) protease affects the development of protease inhibitor (PI) resistance. HIV-1 protease is a main target of current antiretroviral therapies. Our present results demonstrate that the wild-type (WT) virus and a virus with recoded protease exhibited different patterns of resistance mutations after PI treatment. Nevertheless, the developed PI resistance phenotypes were indistinguishable between the recoded virus and the WT virus, suggesting that the HIV-1 strain with synonymously recoded protease and the WT virus are equally robust and evolvable.

Quantitative characterization of hepatitis delta virus genome edition by next-generation sequencing

AIM

To determine the capacity of next-generation sequencing (NGS) for quantifying edited and unedited HDV populations, and to confirm if edition is a general phenomenon taking place along the entire HDV region analyzed, as we previously reported (Homs M et al. PLoS One 2016, 11, e0158557).

METHODS

Four serum samples from 4 patients with chronic HDV/HBV infection were included in the study. The region selected for analysis covered 360 nucleotides (nt), positions 910-1270 of the HDV genome, which included the HDAg ORF editing site (nt 1014 within codon 196). Quantification of edited and unedited genomes was performed by molecular cloning and Sanger sequencing and by NGS. To evaluate the reliability of the NGS values obtained, we combined a clone with an edited codon and one with an unedited codon in known percentages in a series of artificial mixtures, which were then analyzed by NGS. In addition, we determined the nt changes occurring over the complete amplified region after excluding the editing codon (196) to evaluate edition along it.

RESULTS

In total, 11,208 quality-filtered sequences were obtained in the 4 samples. The 95% confidence intervals for the proportions of unedited populations by molecular cloning and NGS were overlapping, and those of cloning were wider, indicating that they are comparable and that NGS is more precise than cloning. Unedited genomes predominated over edited ones in all 4 samples analyzed by NGS and in 3 of the 4 samples analyzed by molecular cloning. In total, 83,276 quality-filtered sequences were obtained from the artificial mixtures. Percentages of the two viral populations detected by NGS in these mixtures were comparable to the expected percentages. Evaluation of edition along the HDV coding region showed that transitions were more frequent than transversions, accounting for 63.09% and 36.91%, respectively. Interestingly, among the 4 possible transition-type changes, G:A and A:G accounted for 73.86% of the total.

CONCLUSION

Next-generation sequencing proved useful to quantify edited and unedited HDV genomes, and provided relevant information on the HDV quasispecies.

High-frequency, low-coverage "false positives" mutations may be true in GS Junior sequencing studies

The GS Junior sequencer provides simplified procedures for library preparation and data processing. Errors in pyrosequencing generate some biases during library construction and emulsion PCR amplification. False-positive mutations are identified by related characteristics described in the manufacturer’s manual, and some detected mutations may have ‘borderline’ characteristics when they are detected in few reads or at low frequency. Among these mutations, however, some may be true positives. This study aimed to improve the accuracy of identifying true positives among mutations with borderline false-positive characteristics detected with GS Junior sequencing. Mutations with the borderline features were tested for validity with Sanger sequencing. We examined 10 mutations detected in coverages <20-fold at frequencies >30% (group A) and 16 mutations detected in coverages >20-fold at frequencies < 30% (group B). In group A, two mutations were not confirmed, and two mutations with 100% frequency were confirmed as heterozygous alleles. No mutation in group B was confirmed. The two groups had significantly different false-positive prevalences (p = 0.001). These results suggest that mutations detected at frequencies less than 30% can be confidently identified as false-positives but that mutations detected at frequencies over 30%, despite coverages less than 20-fold, should be verified with Sanger sequencing.

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.

Resistance to hepatitis C virus. Implications and therapeutic options

We are currently living in an unprecedented era of hepatitis C treatment with the availability of highly effective drugs yielding minimal side effects. The problem we currently face is the retreatment of patients refractory to these drugs. Although several factors can influence treatment failure, this review focuses on antiviral resistance. Resistance-associated substitutions may be identified at baseline or be treatment-emergent. The latter seem to be more clinically relevant and must be studied in the event of treatment failure (no virological response). In this article, we present the latest data from clinical trials and studies in real-life clinical practice. Finally, based on this current evidence, we propose some recommendations for the management and retreatment of these patients.

Phylogenetic analysis of an epidemic outbreak of acute hepatitis C in HIV-infected patients by ultra-deep pyrosequencing

BACKGROUND

The incidence of acute hepatitis C (AHC) among HIV-infected men who have sex with men (MSM) has increased significantly in the last 10 years. Several studies point to a social and sexual network of HIV-positive MSM that extends internationally.

OBJECTIVES

The aim of our study was to investigate the dynamics of HCV transmission in an outbreak of AHC in HIV-infected MSM in Barcelona by ultra-deep pyrosequencing.

STUDY DESIGN

Between 2008 and 2013, 113 cases of AHC in HIV-infected MSM were diagnosed in the Infectious Diseases Unit, Hospital Clínic, Barcelona. Massive sequencing was performed using the Roche 454 GS Junior platform. To define possible transmission networks, maximum likelihood phylogenetic trees were constructed, and levels of genetic diversity within and among patients were compared.

RESULTS

Among the 70 cases analyzed, we have identified 16 potential clusters of transmission: 8 for genotype 1a (23 cases involved), 1 for genotype 1b (3 cases) and 7 for genotype 4d (27 cases). Although the initial phylogenetic reconstruction suggested a local transmission cluster of HCV gt4d, our approach based on low genetic differentiation did not corroborate it. Indeed, gt4d strains formed 4 independent groups related to patients from other countries.

CONCLUSIONS

Frequent clustering of HIV-positive MSM shows that HCV infection has spread through a local network in Barcelona. This outbreak is related to a large international HCV transmission network among MSM. Public health efforts are needed to reduce HCV transmission among this high-risk group.

High resolution sequencing of hepatitis C virus reveals limited intra-hepatic compartmentalization in end-stage liver disease

BACKGROUND & AIMS

The high replication and mutation rate of hepatitis C virus (HCV) results in a heterogeneous population of viral sequences in vivo. HCV replicates in the liver and infected hepatocytes occur as foci surrounded by uninfected cells that may promote compartmentalization of viral variants. Given recent reports showing interferon stimulated gene (ISG) expression in chronic hepatitis C, we hypothesized that local interferon responses may limit HCV replication and evolution.

METHODS

To investigate the spatial influence of liver architecture on viral replication we measured HCV RNA and ISG mRNA from each of the 8 Couinaud segments of the liver from 21 patients undergoing liver transplant.

RESULTS

HCV RNA and ISG mRNA levels were comparable across all sites from an individual liver but showed up to 500-fold difference between patients. Importantly, there was no association between ISG and HCV RNA expression across all sites in the liver or plasma. Deep sequencing of HCV RNA isolated from the 8 hepatic sites from two subjects showed a similar distribution of viral quasispecies across the liver and uniform sequence diversity. Single genome amplification of HCV E1E2-envelope clones from 6 selected patients at 2 hepatic sites supported these data and showed no evidence for HCV compartmentalization.

CONCLUSIONS

We found no differences between the hepatic and plasma viral quasispecies in all patients sampled. We conclude that in end-stage liver disease HCV RNA levels and the genetic pool of HCV envelope sequences are indistinguishable between distant sites in the liver and plasma, arguing against viral compartmentalization.

LAY SUMMARY

HCV is an RNA virus that exists as a quasispecies of closely related genomes that are under continuous selection by host innate and adaptive immune responses and antiviral drug therapy. The primary site of HCV replication is the liver and yet our understanding of the spatial distribution of viral variants within the liver is limited. High resolution sequencing of HCV and monitoring of innate immune responses at multiple sites across the liver identified a uniform pattern of diversity and argues against viral compartmentalization.

Management of direct antiviral agent failures

The current standard of care for patients with chronic hepatitis C virus (HCV) infection is a combination of direct-acting antiviral agents (DAAs). Most HCV patients treated with these drugs achieve viral elimination, but 1% to 15% fail to attain this objective. Treatment failures are usually related to relapse, and less often to on-treatment viral breakthrough. HCV drug resistant associated substitutions are detected in most patients who do not eliminate the virus. The risk of developing these variants depends on host- and virus-related factors, the properties of the drugs used, and the treatment strategies applied. Patients who carry Resistant Associated Substitutions (RASs) may not obtain benefits from treatment, and are at a risk of disease progression. Whether HCV RASs persist depends on their type: NS3-4A variants often disappear gradually after DAA therapy is stopped, whereas NS5A variants tend to persist for more than 2 years. The best way to prevent emergence of resistant variants is to eliminate the virus at the first treatment using highly potent DAAs with genetic barriers to resistance. For those who fail an NS5A inhibitor, deferral of treatment is recommended pending the availability of additional data if they do not have cirrhosis or reasons for urgent re-treatment. If re-treatment is needed, the most commonly used strategy is sofosbuvir as backbone therapy plus a drug from a class other than that previously used, for 24 weeks. Unless it is contraindicated, weight-based ribavirin should also be added. If available, nucleotide-based (eg, sofosbuvir) triple or quadruple DAA regimens may be considered. The optimal treatment for patients who fail an NS5A inhibitor and those with multidrug-resistant variants remains to be defined, and research efforts should continue to focus on treatment for these patients.

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.

Hepatitis C virus: life cycle in cells, infection and host response, and analysis of molecular markers influencing the outcome of infection and response to therapy

Hepatitis C virus (HCV) is a major global health burden accounting for around 170 million chronic infections worldwide. Since its discovery, which dates back to about 30 years ago, many details of the viral genome organization and the astonishing genetic diversity have been unveiled but, owing to the difficulty of culturing HCV in vitro and obtaining fully susceptible yet immunocompetent in vivo models, we are still a long way from the full comprehension of viral life cycle, host cell pathways facilitating or counteracting infection, pathogenetic mechanisms in vivo, and host defences. Here, we illustrate the viral life cycle into cells, describe the interplay between immune and genetic host factors shaping the course of infection, and provide details of the molecular approaches currently used to genotype, monitor replication in vivo, and study the emergence of drug-resistant viral variants.

Translating RNA sequencing into clinical diagnostics: opportunities and challenges

RNA-based measurements have the potential for application across diverse areas of human health, including disease diagnosis, prognosis and therapeutic selection. Current clinical applications include infectious diseases, cancer, transplant medicine and fetal monitoring.

RNA sequencing (RNA-seq) allows for the detection of a wide variety of RNA species, including mRNA, non-coding RNA, pathogen RNA, chimeric gene fusions, transcript isoforms and splice variants, and provides the capability to quantify known, pre-defined RNA species and rare RNA transcript variants within a sample. In addition to differential expression and detection of novel transcripts, RNA-seq also supports the detection of mutations and germline variation for hundreds to thousands of expressed genetic variants, facilitating assessment of allele-specific expression of these variants.

Circulating RNAs and small regulatory RNAs, such as microRNAs, are very stable. These RNA species are vigorously being tested for their potential as biomarkers. However, there are currently few agreed upon methods for isolation or quantitative measurements and a current lack of quality controls that can be used to test platform accuracy and sample preparation quality.

Analytical, bioinformatic and regulatory challenges exist, and ongoing efforts toward the establishment of benchmark standards, assay optimization for clinical conditions and demonstration of assay reproducibility are required to expand the clinical utility of RNA-seq.

RNA sequencing (RNA-seq) is a powerful approach for comprehensive analyses of transcriptomes. This Review describes the widespread potential applications of RNA-seq in clinical medicine, such as detecting disease-associated mutations and gene expression disruptions, as well as characteristic non-coding RNAs, circulating extracellular RNAs or pathogen RNAs. The authors also highlight the challenges in adopting RNA-seq routinely into clinical practice.

With the emergence of RNA sequencing (RNA-seq) technologies, RNA-based biomolecules hold expanded promise for their diagnostic, prognostic and therapeutic applicability in various diseases, including cancers and infectious diseases. Detection of gene fusions and differential expression of known disease-causing transcripts by RNA-seq represent some of the most immediate opportunities. However, it is the diversity of RNA species detected through RNA-seq that holds new promise for the multi-faceted clinical applicability of RNA-based measures, including the potential of extracellular RNAs as non-invasive diagnostic indicators of disease. Ongoing efforts towards the establishment of benchmark standards, assay optimization for clinical conditions and demonstration of assay reproducibility are required to expand the clinical utility of RNA-seq.

Mechanisms of Hepatitis C Viral Resistance to Direct Acting Antivirals

There has been a remarkable transformation in the treatment of chronic hepatitis C in recent years with the development of direct acting antiviral agents targeting virus encoded proteins important for viral replication including NS3/4A, NS5A and NS5B. These agents have shown high sustained viral response (SVR) rates of more than 90% in phase 2 and phase 3 clinical trials; however, this is slightly lower in real-life cohorts. Hepatitis C virus resistant variants are seen in most patients who do not achieve SVR due to selection and outgrowth of resistant hepatitis C virus variants within a given host. These resistance associated mutations depend on the class of direct-acting antiviral drugs used and also vary between hepatitis C virus genotypes and subtypes. The understanding of these mutations has a clear clinical implication in terms of choice and combination of drugs used. In this review, we describe mechanism of action of currently available drugs and summarize clinically relevant resistance data.

Evolution of Newcastle Disease Virus Quasispecies Diversity and Enhanced Virulence after Passage through Chicken Air Sacs

It has been reported that lentogenic Newcastle disease virus (NDV) isolates have the potential to become velogenic after their transmission and circulation in chickens, but the underlying mechanism is unclear. In this study, a highly velogenic NDV variant, JS10-A10, was generated from the duck-origin lentogenic isolate JS10 through 10 consecutive passages in chicken air sacs. The velogenic properties of this selected variant were determined using mean death time (MDT) assays, intracerebral pathogenicity index (ICPI), the intravenous pathogenicity index (IVPI), histopathology, and the analysis of host tissue tropism. In contrast, JS10 remained lentogenic after 20 serial passages in chicken eggs (JS10-E20). The JS10, JS10-A10, and JS10-E20 genomes were sequenced and found to be nearly identical, suggesting that both JS10-A10 and JS10-E20 were directly generated from JS10. To investigate the mechanism for virulence enhancement, the partial genome covering the F0 cleavage site of JS10 and its variants were analyzed using ultradeep pyrosequencing (UDPS) and the proportions of virulence-related genomes in the quasispecies were calculated. Velogenic NDV genomes accumulated as a function of JS10 passaging through chicken air sacs. Our data suggest that lentogenic NDV strains circulating among poultry might be a risk factor to future potential velogenic NDV outbreaks in chickens.

IMPORTANCE

An avirulent isolate, JS10, was passaged through chicken air sacs and embryos, and the pathogenicity of the variants was assessed. A virulent variant, JS10-A10, was generated from consecutive passage in air sacs. We developed a deep-sequencing approach to detect low-frequency viral variants across the NDV genome. We observed that virulence enhancement of JS10 was due to the selective accumulation of velogenic quasispecies and the concomitant disappearance of lentogenic quasispecies. Our results suggest that because it is difficult to avoid contact between natural waterfowl reservoirs and sensitive poultry operations, circulating lentogenic NDV strains may represent a potential reservoir for emergent velogenic NDV strains that could cause outbreaks in chickens.

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.

Management of direct-acting antiviral agent failures

Failure to respond to the approved combinations of multiple direct-acting antiviral agents is relatively low in hepatitis C virus treatment registration studies, with rates of 1% to 7%, depending on the patients' baseline characteristics. In real life, failure is slightly higher, likely because of lower compliance. Treatment failures are usually related to relapse and less often to on-treatment viral breakthrough. Hepatitis C drug-resistant variants are detected in most patients who do not achieve viral eradication. The risk of developing these variants depends on host- and virus-related factors, the properties of the drugs used, and the treatment strategies applied. Patients who carry resistance-associated variants may not obtain benefits from treatment and are at risk of disease progression and transmission of the variants. Whether hepatitis C resistance-associated variants persist depends on their type: NS3-4A variants often disappear gradually after therapy is stopped, whereas NS5A variants tend to persist for more than 2 years. The best way to prevent emergence of resistant variants is to eliminate the virus at the first treatment using highly potent antivirals with genetic barriers to resistance. In patients failing first-generation protease inhibitors, combination therapies with sofosbuvir and NS5 inhibitors have proven effective. Some salvage regimens can be shortened to 12 weeks by addition of ribavirin. The optimal treatment for patients who fail an NS5A inhibitor and those with multidrug-resistant variants remains to be defined, and research efforts should continue to focus on treatment for these patients.

Resistance of Hepatitis C Virus to Inhibitors: Complexity and Clinical Implications

Selection of inhibitor-resistant viral mutants is universal for viruses that display quasi-species dynamics, and hepatitis C virus (HCV) is no exception. Here we review recent results on drug resistance in HCV, with emphasis on resistance to the newly-developed, directly-acting antiviral agents, as they are increasingly employed in the clinic. We put the experimental observations in the context of quasi-species dynamics, in particular what the genetic and phenotypic barriers to resistance mean in terms of exploration of sequence space while HCV replicates in the liver of infected patients or in cell culture. Strategies to diminish the probability of viral breakthrough during treatment are briefly outlined.

Dynamic of Mixed HCV Infection in Plasma and PBMC of HIV/HCV Patients Under Treatment With Peg-IFN/Ribavirin

The extent of mixed hepatitis C virus (HCV) genotype in different compartments (plasma and peripheral blood mononuclear cell, PBMC) and possible association with treatment efficacy in HIV/HCV coinfected patients remains to be unknown.The objective of this study was to elucidate the frequency of mixed genotype infection (MG), its profile in different compartments during anti-HCV treatment, and the possible influence of different genotypes on the response rate.The compartmentalization of HCV population was investigated by next-generation sequencing in 19 HIV/HCV coinfected patients under anti-HCV treatment with peginterferon/ribavirin (P-R). Ten individuals were nonresponder (NR) or relapser (RE) to P-R treatment and 9 had a sustained virological response (SVR).Eleven/nineteen (58%) patients had MG in plasma compartment. Ten or 12 patients infected by a difficult to treat genotype (DTG) 1 or 4 as dominant strain, had an MG, whereas only 1/7 individuals infected by easy to treat genotype (ETG) harbored a mixed genotype, P = 0.006. HCV-RNA was more frequently detected in PBMC of NR (10/10) than in those of SVR (5/9), P = 0.032. Mixed genotype infection was detected in 6/15 (40%) PBMC-positive cases and was not associated with P-R treatment response. By multivariate analysis, MG in plasma samples was the most important viral factor affecting the treatment response (P = 0.0237).Detection of MG in plasma of HIV/HCV coinfected patients seems to represent the major determinant of response to P-R treatment. This finding may have important clinical implication in light of the new therapeutic approach in HIV/HCV coinfected individuals suggesting that combination treatment with direct acting antivirals could be less effective in MG.

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.

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.

Ultradeep pyrosequencing of NS3 to predict response to triple therapy with protease inhibitors in previously treated chronic hepatitis C patients

Despite the gain in sustained virological responses (SVR) provided by protease inhibitors (PIs), failures still occur. The aim of this study was to determine if a baseline analysis of the NS3 region using ultradeep pyrosequencing (UDPS) can help to predict an SVR. Serum samples from 40 patients with previously nonresponding genotype 1 chronic hepatitis C who were retreated with triple therapy, including a PI, were analyzed. Baseline UDPS of the NS3 gene was performed on plasma and peripheral blood mononuclear cells (PBMC). Mutations conferring resistance to PIs were sought. The overall diversity of the quasispecies was evaluated by calculating the Shannon entropy (SE). Resistance mutations were found in plasma and PBMC but were not discriminating enough to predict an SVR. NS3 quasispecies heterogeneity was significantly lower at baseline in patients achieving an SVR than in those not achieving an SVR (SE of 26.98 ± 16.64 × 10(-3) versus 44.93 ± 19.58 × 10(-3), P = 0.0047). With multivariate analysis, the independent predictors of an SVR were fibrosis of stage F ≤2 (odds ratio [OR], 13.3; 95% confidence interval [CI], 1.25 to 141.096; P < 0.03) and SE below the median (OR, 5.4; 95% CI, 1.22 to 23.87; P < 0.03). More than the presence of minor mutations at the baseline in plasma or in PBMC, the NS3 viral heterogeneity determined by UDPS is an independent factor for an SVR in previously treated patients receiving triple therapy that includes a PI.

Complex Genotype Mixtures Analyzed by Deep Sequencing in Two Different Regions of Hepatitis B Virus

This study assesses the presence and outcome of genotype mixtures in the polymerase/surface and X/preCore regions of the HBV genome in patients with chronic hepatitis B virus (HBV) infection. Thirty samples from ten chronic hepatitis B patients were included. The polymerase/surface and X/preCore regions were analyzed by deep sequencing (UDPS) in the first available sample at diagnosis, a pre-treatment sample, and a sample while under treatment. HBV genotype was determined by phylogenesis. Quasispecies complexity was evaluated by mutation frequency and nucleotide diversity. The polymerase/surface and X/preCore regions were validated for genotyping from 113 GenBank reference sequences. UDPS yielded a median of 10,960 sequences per sample (IQR 16,645) in the polymerase/surface region and 11,595 sequences per sample (IQR 14,682) in X/preCore. Genotype mixtures were more common in X/preCore (90%) than in polymerase/surface (30%) (p<0.001). On X/preCore genotyping, all samples were genotype A, whereas polymerase/surface yielded genotypes A (80%), D (16.7%), and F (3.3%) (p = 0.036). Genotype changes in polymerase/surface were observed in four patients during natural quasispecies dynamics and in two patients during treatment. There were no genotype changes in X/preCore. Quasispecies complexity was higher in X/preCore than in polymerase/surface (p = 0.004). The results provide evidence of genotype mixtures and differential genotype proportions in the polymerase/surface and X/preCore regions. The genotype dynamics in HBV infection and the different patterns of quasispecies complexity in the HBV genome suggest a new paradigm for HBV genotype classification.

Computational framework for next-generation sequencing of heterogeneous viral populations using combinatorial pooling

MOTIVATION

Next-generation sequencing (NGS) allows for analyzing a large number of viral sequences from infected patients, providing an opportunity to implement large-scale molecular surveillance of viral diseases. However, despite improvements in technology, traditional protocols for NGS of large numbers of samples are still highly cost and labor intensive. One of the possible cost-effective alternatives is combinatorial pooling. Although a number of pooling strategies for consensus sequencing of DNA samples and detection of SNPs have been proposed, these strategies cannot be applied to sequencing of highly heterogeneous viral populations.

RESULTS

We developed a cost-effective and reliable protocol for sequencing of viral samples, that combines NGS using barcoding and combinatorial pooling and a computational framework including algorithms for optimal virus-specific pools design and deconvolution of individual samples from sequenced pools. Evaluation of the framework on experimental and simulated data for hepatitis C virus showed that it substantially reduces the sequencing costs and allows deconvolution of viral populations with a high accuracy.

AVAILABILITY AND IMPLEMENTATION

The source code and experimental data sets are available at http://alan.cs.gsu.edu/NGS/?q=content/pooling.

Analysis of the evolution and structure of a complex intrahost viral population in chronic hepatitis C virus mapped by ultradeep pyrosequencing

Hepatitis C virus (HCV) causes chronic infection in up to 50% to 80% of infected individuals. Hypervariable region 1 (HVR1) variability is frequently studied to gain an insight into the mechanisms of HCV adaptation during chronic infection, but the changes to and persistence of HCV subpopulations during intrahost evolution are poorly understood. In this study, we used ultradeep pyrosequencing (UDPS) to map the viral heterogeneity of a single patient over 9.6 years of chronic HCV genotype 4a infection. Informed error correction of the raw UDPS data was performed using a temporally matched clonal data set. The resultant data set reported the detection of low-frequency recombinants throughout the study period, implying that recombination is an active mechanism through which HCV can explore novel sequence space. The data indicate that polyvirus infection of hepatocytes has occurred but that the fitness quotients of recombinant daughter virions are too low for the daughter virions to compete against the parental genomes. The subpopulations of parental genomes contributing to the recombination events highlighted a dynamic virome where subpopulations of variants are in competition. In addition, we provide direct evidence that demonstrates the growth of subdominant populations to dominance in the absence of a detectable humoral response.

IMPORTANCE

Analysis of ultradeep pyrosequencing data sets derived from virus amplicons frequently relies on software tools that are not optimized for amplicon analysis, assume random incorporation of sequencing errors, and are focused on achieving higher specificity at the expense of sensitivity. Such analysis is further complicated by the presence of hypervariable regions. In this study, we made use of a temporally matched reference sequence data set to inform error correction algorithms. Using this methodology, we were able to (i) detect multiple instances of hepatitis C virus intrasubtype recombination at the E1/E2 junction (a phenomenon rarely reported in the literature) and (ii) interrogate the longitudinal quasispecies complexity of the virome. Parallel to the UDPS, isolation of IgG-bound virions was found to coincide with the collapse of specific viral subpopulations.

Inference with viral quasispecies diversity indices: clonal and NGS approaches

Given the inherent dynamics of a viral quasispecies, we are often interested in the comparison of diversity indices of sequential samples of a patient, or in the comparison of diversity indices of virus in groups of patients in a treated versus control design. It is then important to make sure that the diversity measures from each sample may be compared with no bias and within a consistent statistical framework. In the present report, we review some indices often used as measures for viral quasispecies complexity and provide means for statistical inference, applying procedures taken from the ecology field. In particular, we examine the Shannon entropy and the mutation frequency, and we discuss the appropriateness of different normalization methods of the Shannon entropy found in the literature. By taking amplicons ultra-deep pyrosequencing (UDPS) raw data as a surrogate of a real hepatitis C virus viral population, we study through in-silico sampling the statistical properties of these indices under two methods of viral quasispecies sampling, classical cloning followed by Sanger sequencing (CCSS) and next-generation sequencing (NGS) such as UDPS. We propose solutions specific to each of the two sampling methods-CCSS and NGS-to guarantee statistically conforming conclusions as free of bias as possible.

CONTACT

josep.gregori@gmail.com Supplementary information: Supplementary data are available at Bioinformatics online.

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.