- Research article
- Open Access
- Open Peer Review
Deep sequencing of hepatitis C virus hypervariable region 1 reveals no correlation between genetic heterogeneity and antiviral treatment outcome
BMC Infectious Diseasesvolume 14, Article number: 389 (2014)
Hypervariable region 1 (HVR1) contained within envelope protein 2 (E2) gene is the most variable part of HCV genome and its translation product is a major target for the host immune response. Variability within HVR1 may facilitate evasion of the immune response and could affect treatment outcome. The aim of the study was to analyze the impact of HVR1 heterogeneity employing sensitive ultra-deep sequencing, on the outcome of PEG-IFN-α (pegylated interferon α) and ribavirin treatment.
HVR1 sequences were amplified from pretreatment serum samples of 25 patients infected with genotype 1b HCV (12 responders and 13 non-responders) and were subjected to pyrosequencing (GS Junior, 454/Roche). Reads were corrected for sequencing error using ShoRAH software, while population reconstruction was done using three different minimal variant frequency cut-offs of 1%, 2% and 5%. Statistical analysis was done using Mann–Whitney and Fisher’s exact tests.
Complexity, Shannon entropy, nucleotide diversity per site, genetic distance and the number of genetic substitutions were not significantly different between responders and non-responders, when analyzing viral populations at any of the three frequencies (≥1%, ≥2% and ≥5%). When clonal sample was used to determine pyrosequencing error, 4% of reads were found to be incorrect and the most abundant variant was present at a frequency of 1.48%. Use of ShoRAH reduced the sequencing error to 1%, with the most abundant erroneous variant present at frequency of 0.5%.
While deep sequencing revealed complex genetic heterogeneity of HVR1 in chronic hepatitis C patients, there was no correlation between treatment outcome and any of the analyzed quasispecies parameters.
Hepatitis C virus (HCV) circulates within the infected host as a pool of related but distinct genetic variants (quasispecies); . The genetic variability is mainly generated by viral RNA-dependent RNA polymerase (RdRp) which lacks a proof-reading activity . Genes encoding envelope E1 and E2 proteins, especially the hypervariable 1 region (HVR1) of E2, display the highest genetic variability within the whole HCV genome . HVR1 contains sequences encoding important immune epitopes; thus genetic variability within this region may facilitate evasion of the immune responses and is largely shaped by the immune pressure of the host [4–8]. Complexity and evolution of HVR1 quasispecies was reported to be predictive factor of the outcome of natural infection [9, 10].
Antiviral treatment protocols using interferon and ribavirin have limited efficacy and are plagued by side effects, which often require premature discontinuation of therapy. Factors known to be associated with treatment outcome include both host (i.e. IL28B gene polymorphisms, race, sex, age) as well as viral factors (genotype, serum load and genetic heterogeneity); [11–13].
Interferon and ribavirin treatment is based largely on direct antiviral effect as well as immunomodulation . Thus, HVR1 heterogeneity could facilitate treatment failure since coexistence of multiple antigenic variants could increase the probability of positive selection of those effectively evading immune pressure induced by treatment [15, 16]. However, despite attempts to correlate HVR1 heterogeneity with antiviral treatment outcome, published studies are inconclusive [17–20].
Recent years brought the advent of ultra-deep sequencing techniques which enable parallel sequencing of multiple sequences present in a sample, thus providing better insight into the quasispecies phenomenon. Pyrosequencing (454/Roche), one of the available deep sequencing platforms, is capable of reading sequences up to 1 kb, and it was used successfully for sequence analysis of human immunodeficiency virus (HIV) and HCV [21–25].
Similarly, our previous analysis of HVR1 in chronic HCV infection confirmed the utility of pyrosequencing for HCV haplotypes inference, including identification of very rare variants constituting as little as 0.1% of the whole population .
The present study employed pyrosequencing to explore HVR1 complexity and variability in pretreatment serum samples of patients treated with pegylated interferon α (PEG-IFN α ) and ribavirin. We demonstrated that complexity, Shannon entropy, nucleotide diversity per site, genetic distance and the number of genetic substitutions were not significantly different between responders and non-responders, when analyzing populations present at ≥1%, ≥2% and ≥5% frequency.
Our prospective study involved 95 chronic hepatitis C patients undergoing treatment at the Outpatient Clinic of the Hospital for Infectious Diseases in Warsaw from June 2010 to December 2012. Out of this cohort, twenty five patients were selected according to the following criteria: chronic infection with genotype 1b HCV, no previous antiviral treatment, no co-infection with HBV or HIV, no history of intravenous drugs use. In addition, patients had to achieve complete early viral response (cEVR), defined as undetectable HCV RNA in the serum after 12 weeks of treatment and, subsequently, sustained viral response (SVR) defined as undetectable HCV RNA in the serum 6 months post-treatment (responders, n = 12), or experience no viral load reduction ≥ 2 log at week 12 of treatment (non-responders, NR, n = 13). No statistically significant differences were found between responders and non-responders in mean alanine aminotransferase activity, pretreatment viral load, age, liver grading and staging or sex (Table 1). Viral load was measured by RealTime HCV assay (Abbott), sensitivity: 12 IU/mL, while qualitative evaluation was performed by COBAS Amplicor HCV Test, v2.0 (Roche Diagnostics) which has sensitivity limit of 50 IU/ml. Treatment consisted of pegylated interferon α (Pegasys, Roche), 180 μg per week, n = 13 or Pegintron (Schering-Plough) at dose 1,5 μg/kg of body weight, n = 12 and Ribavirin (Copegus, Roche), 1000 mg/day (body mass < 75 kg) or 1200 mg/day (body mass > 75 kg), n = 13 or Rebetol (Schering-Plough), 800 mg/day (body mass < 64 kg), 1000 mg/day (body mass 65–85 kg), 1200 mg (body mass 86–105 kg) or 1400 mg (body mass > 105 kg), n = 12. Responders were treated for 48 weeks, whereas in non-responders the therapy was stopped after 12 weeks. The study was approved by the Institutional Bioethical Committee (consent No KB/107/2010) and all patients provided informed consent.
HVR1 amplification was performed from pretreatment serum samples as described previously . In brief, viral RNA was extracted from 250 μl of serum by modified guanidinium thiocyanate-phenol/chlorophorm method, then subjected to reverse transcription at 37°C for 30 minutes using AccuScript High Fidelity Reverse Transcriptase (Agilent Technologies). A fragment of E2 region containing HVR1 was amplified in two-step PCR using FastStart High Fidelity Taq DNA Polymerase (Roche). Primers for the second round PCR contained tags recognized by GS Junior sequencing platform, standard 10-nucleotide multiplex identifiers (MID) and target-specific sequence.
Cloned HVR1 sequence
To determine the inherent sequencing error, amplified HVR1 from one sample was purified by Wizard SV Genomic DNA Purification System (Promega) and cloned into TOPO TA vector using TOPO TA Cloning Kit (Invitrogen). Plasmid DNA was extracted from bacterial culture using Quick Plasmid Miniprep Kit (Life technologies). Subsequently, pyrosequencing-specific tags with multiplex identifier (MID) were introduced by means of PCR using plasmid sequence as a target and sample was subjected to pyrosequencing.
Each amplicon was purified from agarose gel by QIAquick Gel Extraction kit (Qiagen) and then by Agencourt AMPure XP beads (Beckman Coulter) using 1.6:1 ratio of beads to sample. Products were quantified by dsDNA HS Qubit® Assay Kit (Life Technologies), fourteen samples were pooled in equivalent amounts and of 3 × 107 DNA copies were subjected to emulsion PCR using GS Junior Titanium emPCR Kit (Lib-A). After initial denaturation at 94°C for 1 minute, the reaction was run for 50 cycles of 94°C for 30 seconds, 58°C for 4 minutes and 30 seconds, and 68°C for 30 seconds. DNA library beads enrichment was carried out according to the emPCR Amplification Method Manual Lib-A (Roche), with the exception that the number of bead washes was 15. The required input of 500 000 enriched beads was loaded onto the Pico Titer Plate (PTP) and sequencing was carried out for 200 cycles using full processing mode for amplicons (GS Junior Sequencer, 454/Roche). In total, two independent pyrosequencing runs were performed (14 samples with specific MID were pooled in each).
Reads of individual samples were demultiplexed, sequencing errors were corrected and haplotypes inferred using the program diri_sampler from the ShoRAH software . Error correction included mismatches as well as insertions and deletions. Subsequently, haplotypes were aligned to the 1b HCV reference sequence (GenBank:AJ406073) and translated into amino acid sequences by MEGA (Molecular Evolutionary Genetics Analysis), version 5.0 . Phylogenetic trees were constructed according to the Maximum Likelihood method based on the Tamura-Nei model  using MEGA 5.0. Genetic diversity parameters were assessed in HVR1 populations of frequency ≥1%, 2% and 5% by DNA SP version 5 . Such cut-off approach facilitated interpatient comparison of sequence populations of different coverages. HVR1 complexity was represented by the number of haplotypes above each frequency cut-off. Nucleotide diversity per site and the number of substitutions were assessed using DNA SP version 5 with respect to the reference sequence (GenBank:AJ406073). Genetic distances in HCV HVR1 populations were assessed by MEGA. Shannon entropy was calculated according to the following equation:
N – number of observations (haplotypes),
f i - frequency of haplotypes
Differences in age, alanine aminotransferase activity, viral load, HVR1 complexity, diversity, number of substitutions within HVR1, Shannon entropy, genetic distance, number of polymorphic amino acid positions and number of inner nodes in phylogenetic trees were compared using Mann–Whitney test, while proportions were compared by Fisher’s exact test.
Estimation of pyrosequencing and amplification errors based on cloned HVR1 sequence
Sequencing of cloned HVR1 fragment provided 3178 reads. After grouping identical reads together, 12 variants were identified (Table 2). Only 96% of reads were identical to the original template. Among 11 erroneous variants, the most abundant constituted 1.48% of all reads, whereas the least abundant was present at a frequency of 0.06% (Figure 1).
Errors included insertions (83.3%), substitutions (12.5%) and deletions (4.2%). Probability of error occurrence per base was estimated to be 0.04% for insertion, 0.006% for substitution and 0.002% for deletion. Fifty one percent of insertions occurred at homopolymeric regions (four repeats of T). Altogether, the probability of any error per base was 0.05%.
After error correction performed with ShoRAH, four variants were identified: one identical to the template at 99.0% frequency, and three erroneous variants present at frequency of 0.5%, 0.3% and 0.2%, respectively.
Characteristics of deep sequencing
Over 15 million nucleotides were sequenced (Table 3). After demultiplexing, the median (IQR) of assigned reads was 2540 (2488) per patient sample - 2540 (1790) in responders and 1230 (2816) in non-responders. Following ShoRAH reconstruction, the mean number of haplotypes obtained per patient was 30.6 (38.4 in responders and 23.4 in non-responders). Most abundant haplotype constituted 57.09%, whereas the least abundant only 0.1%. The number of reconstructed haplotypes depends on several factors, including coverage, frequency of the haplotypes and their distance. In order to make a reliable comparison in different patients, we introduced a threshold to the haplotype frequency. The frequency thresholds explored were 1%, 2% and 5%.
HVR1 genetic heterogeneity
HVR1 complexity at ≥5% haplotype frequency cut-off was slightly lower in responders (R) than non-responders (NR); (4.4 vs 5.3); (Table 4, Figure 2). Likewise, mean Shannon entropy, mean genetic distance of HVR1 populations and mean number of genetic substitutions and nucleotide diversity per site were also lower in the former group (Table 4, Figure 2). However, these differences did not reach statistical significance. Similarly, when the above analysis was repeated at ≥2% and ≥1% frequency cut-offs, no statistically significant differences were either found.
Amino acid variability of HVR1
Within 27 amino acid stretch of HVR1, responders were found to have similar mean number of polymorphic amino acid positions (59.3% ± 9.5%) as non-responders (60% ± 11%); (Table 4). Additional file 1 shows multiple sequence alignment of amino acid sequences of HVR1 populations in responders (R) and non-responders to treatment (NR).
Viral populations ≥5% were also analyzed phylogenetically (Figure 3). As shown, populations in non-responders formed more complex patterns of relatedness as manifested by the higher mean number of inner nodes (4.0 ± 2.9 vs 2.9 ± 0.7). Nevertheless, this difference was not statistically significant.
A number of previous studies attempted to correlate HVR1 heterogeneity with antiviral treatment outcome, but their results were usually inconclusive and occasionally even contradictory. These discrepancies could be partly due to the use of different techniques: two most commonly used were single strand conformational polymorphism (SSCP) and clonal Sanger sequencing [17–20, 31–34]. The latter requires extensive cloning to achieve high sensitivity for minor variants detection, a process that is costly and time-consuming. Thus, studies using this technique rarely included significant number of clones per sample, typically attaining only 15-20% sensitivity. While SSCP has been shown to detect variants constituting as little as 3% of the viral population , it is not informative of the nucleotide sequence, the nature of genetic changes or genetic distances between variants. Furthermore, in a mixture of heterogeneous sequences, certain bands may overlap, resulting in underestimation of viral complexity. Our current study, which was based on deep sequencing, overcomes the above shortcomings and represents a novel approach to analysis of HCV heterogeneity.
While our analysis did not find any significant differences in HVR1 heterogeneity between responders and non-responders to antiviral treatment, these results are largely compatible with some previous studies employing SSCP and clonal sequencing. In the study of Pawlotsky et al.  based on single strand conformational polymorphism and in the study of Saludes et al.  based on clonal sequencing, no significant differences in pretreatment HVR1 complexity were observed between responders and non-responders. Similar results were reported in a study of re-treated patients with advanced fibrosis , while Abbate et al.  found that low pretreatment HVR1 heterogeneity correlated with early response (EVR), but not with SVR. A number of other studies found no correlation between HVR1 complexity and treatment outcome [17, 18, 36].
In our study, such HVR1 heterogeneity parameters, as nucleotide diversity per site, genetic distance, and number of nucleotide substitutions also did not differ significantly between responders and non-responders. These findings are similar to several earlier studies [20, 34, 37, 38]. In the only published study using deep sequencing approach, there were no differences in pretreatment complexity parameters (e.g. Shannon entropy) between immediate virological responders and non-responders. However, the final treatment outcome was not reported .
Lack of statistically significant differences in analyzed heterogeneity parameters between responders and non-responders suggest that the heterogeneity generated by minor variants detectable by deep sequencing has no effect on treatment outcome. Alternatively, it may be speculated that the analyzed depth of frequency is still insufficient to detect minor variants whose heterogeneity would have clinical significance.
Some recent studies brought attention to the problem of inherent ultra-deep sequencing errors affecting the detection of minor variants of the quasispecies population [26, 39, 40]. In our analysis, the internal control experiment using cloned HVR1 revealed the overall sequencing error to be 0.05% per nucleotide, comprising mostly of insertions and occurring predominantly in homopolymeric regions. This error rate contributed to the high proportion of erroneous sequences (4% of total reads, the most abundant erroneous variant being present at a frequency of 1.48%). To minimize the risk of including erroneous variants into analysis, we implemented ShoRAH error correction method, which allowed for correction of 99% of reads reducing both the absolute number and frequency of erroneous variants. Thus, error correction methods should be used to facilitate analysis of minor quasispecies by pyrosequencing.
There were no significant differences in the pretreatment HVR1 heterogeneity parameters such as complexity, Shannon entropy, nucleotide diversity per site, genetic distance and the number of genetic substitutions between responders and non-responders. Thus, pretreatment HVR1 quasispecies composition and heterogeneity analysis seems to have limited value for the prediction of treatment outcome.
Martell M, Esteban JI, Quer J, Genesca J, Weiner A, Esteban R, Guardia J, Gomez J: Hepatitis C virus (HCV) circulates as a population of different but closely related genomes: quasispecies nature of HCV genome distribution. J Virol. 1992, 66 (5): 3225-3229.
Duarte EA, Novella IS, Weaver SC, Domingo E, Wain-Hobson S, Clarke DK, Moya A, Elena SF, de la Torre JC, Holland JJ: RNA virus quasispecies: significance for viral disease and epidemiology. Infect Agents Dis. 1994, 3 (4): 201-214.
Kato N, Ootsuyama Y, Ohkoshi S, Nakazawa T, Sekiya H, Hijikata M, Shimotohno K: Characterization of hypervariable regions in the putative envelope protein of hepatitis C virus. Biochem Biophys Res Commun. 1992, 189 (1): 119-127. 10.1016/0006-291X(92)91533-V.
Di Lorenzo C, Angus AG, Patel AH: Hepatitis C virus evasion mechanisms from neutralizing antibodies. Viruses. 2011, 3 (11): 2280-2300.
Kato N, Sekiya H, Ootsuyama Y, Nakazawa T, Hijikata M, Ohkoshi S, Shimotohno K: Humoral immune response to hypervariable region 1 of the putative envelope glycoprotein (gp70) of hepatitis C virus. J Virol. 1993, 67 (7): 3923-3930.
Guglietta S, Garbuglia AR, Pacciani V, Scotta C, Perrone MP, Laurenti L, Spada E, Mele A, Capobianchi MR, Taliani G, Folgori A, Vitelli A, Ruggeri L, Nicosia A, Piccolella E, Del Porto P: Positive selection of cytotoxic T lymphocyte escape variants during acute hepatitis C virus infection. Eur J Immunol. 2005, 35 (9): 2627-2637. 10.1002/eji.200526067.
Zoulim F, Chevallier M, Maynard M, Trepo C: Clinical consequences of hepatitis C virus infection. Rev Med Virol. 2003, 13 (1): 57-68. 10.1002/rmv.371.
Hoofnagle JH: Hepatitis C: the clinical spectrum of disease. Hepatology. 1997, 26 (3 Suppl 1): 15S-20S.
Farci P, Shimoda A, Coiana A, Diaz G, Peddis G, Melpolder JC, Strazzera A, Chien DY, Munoz SJ, Balestrieri A, Purcell RH, Alter HJ: The outcome of acute hepatitis C predicted by the evolution of the viral quasispecies. Science. 2000, 288 (5464): 339-344. 10.1126/science.288.5464.339.
Laskus T, Wilkinson J, Gallegos-Orozco JF, Radkowski M, Adair DM, Nowicki M, Operskalski E, Buskell Z, Seeff LB, Vargas H, Rakela J: Analysis of hepatitis C virus quasispecies transmission and evolution in patients infected through blood transfusion. Gastroenterology. 2004, 127 (3): 764-776. 10.1053/j.gastro.2004.06.005.
McHutchison JG, Poynard T, Pianko S, Gordon SC, Reid AE, Dienstag J, Morgan T, Yao R, Albrecht J: The impact of interferon plus ribavirin on response to therapy in black patients with chronic hepatitis C. The International Hepatitis Interventional Therapy Group. Gastroenterology. 2000, 119 (5): 1317-1323. 10.1053/gast.2000.19289.
Lam NP, Pitrak D, Speralakis R, Lau AH, Wiley TE, Layden TJ: Effect of obesity on pharmacokinetics and biologic effect of interferon-alpha in hepatitis C. Dig Dis Sci. 1997, 42 (1): 178-185. 10.1023/A:1018865928308.
Asselah T, Estrabaud E, Bieche I, Lapalus M, De Muynck S, Vidaud M, Saadoun D, Soumelis V, Marcellin P: Hepatitis C: viral and host factors associated with non-response to pegylated interferon plus ribavirin. Liver Int. 2010, 30 (9): 1259-1269. 10.1111/j.1478-3231.2010.02283.x.
Ijichi S, Izumo S, Nagai M, Shinmyozu K, Hall WW, Osame M: Anti-viral and immunomodulatory effects of interferon-alpha on cultured lymphocytes from patients with human T lymphotropic virus type I-associated myelopathy (HAM/TSP). J Neuroimmunol. 1995, 61 (2): 213-221. 10.1016/0165-5728(95)00101-7.
Shindo M, Hamada K, Koya S, Arai K, Sokawa Y, Okuno T: The clinical significance of changes in genetic heterogeneity of the hypervariable region 1 in chronic hepatitis C with interferon therapy. Hepatology. 1996, 24 (5): 1018-1023. 10.1002/hep.510240507.
Polyak SJ, McArdle S, Liu SL, Sullivan DG, Chung M, Hofgartner WT, Carithers RL, McMahon BJ, Mullins JI, Corey L, Gretch DR, Diaz G, Balestrieri A, Purcell RH: Evolution of hepatitis C virus quasispecies in hypervariable region 1 and the putative interferon sensitivity-determining region during interferon therapy and natural infection. J Virol. 1998, 72 (5): 4288-4296.
Sandres K, Dubois M, Pasquier C, Payen JL, Alric L, Duffaut M, Vinel JP, Pascal JP, Puel J, Izopet J: Genetic heterogeneity of hypervariable region 1 of the hepatitis C virus (HCV) genome and sensitivity of HCV to alpha interferon therapy. J Virol. 2000, 74 (2): 661-668. 10.1128/JVI.74.2.661-668.2000.
Farci P, Strazzera R, Alter HJ, Farci S, Degioannis D, Coiana A, Peddis G, Usai F, Serra G, Chessa L, et al: Early changes in hepatitis C viral quasispecies during interferon therapy predict the therapeutic outcome. Proc Natl Acad Sci U S A. 2002, 99 (5): 3081-3086. 10.1073/pnas.052712599.
Pawlotsky JM, Germanidis G, Frainais PO, Bouvier M, Soulier A, Pellerin M, Dhumeaux D: Evolution of the hepatitis C virus second envelope protein hypervariable region in chronically infected patients receiving alpha interferon therapy. J Virol. 1999, 73 (8): 6490-6499.
Cuevas JM, Torres-Puente M, Jimenez-Hernandez N, Bracho MA, Garcia-Robles I, Carnicer F, Olmo JD, Ortega E, Moya A, Gonzalez-Candelas F: Refined analysis of genetic variability parameters in hepatitis C virus and the ability to predict antiviral treatment response. J Viral Hepat. 2008, 15 (8): 578-590. 10.1111/j.1365-2893.2008.00991.x.
Wang C, Mitsuya Y, Gharizadeh B, Ronaghi M, Shafer RW: Characterization of mutation spectra with ultra-deep pyrosequencing: application to HIV-1 drug resistance. Genome Res. 2007, 17 (8): 1195-1201. 10.1101/gr.6468307.
Simen BB, Simons JF, Hullsiek KH, Novak RM, Macarthur RD, Baxter JD, Huang C, Lubeski C, Turenchalk GS, Braverman MS, Desany B, Rothberg JM, Egholm M, Kozal MJ, Terry Beirn Community Programs for Clinical Research on AIDS: Low-abundance drug-resistant viral variants in chronically HIV-infected, antiretroviral treatment-naive patients significantly impact treatment outcomes. J Infect Dis. 2009, 199 (5): 693-701. 10.1086/596736.
Spear GT, Sikaroodi M, Zariffard MR, Landay AL, French AL, Gillevet PM: Comparison of the diversity of the vaginal microbiota in HIV-infected and HIV-uninfected women with or without bacterial vaginosis. J Infect Dis. 2008, 198 (8): 1131-1140. 10.1086/591942.
Bull RA, Luciani F, McElroy K, Gaudieri S, Pham ST, Chopra A, Cameron B, Maher L, Dore GJ, White PA, Lloyd AR: Sequential bottlenecks drive viral evolution in early acute hepatitis C virus infection. PLoS Pathog. 2011, 7 (9): e1002243-10.1371/journal.ppat.1002243.
Nasu A, Marusawa H, Ueda Y, Nishijima N, Takahashi K, Osaki Y, Yamashita Y, Inokuma T, Tamada T, Fujiwara T, Sato F, Shimizu K, Chiba T: Genetic heterogeneity of hepatitis C virus in association with antiviral therapy determined by ultra-deep sequencing. PLoS ONE. 2011, 6 (9): e24907-10.1371/journal.pone.0024907.
Caraballo Cortes K, Zagordi O, Laskus T, Ploski R, Bukowska-Osko I, Pawelczyk A, Berak H, Radkowski M: Ultradeep pyrosequencing of hepatitis C virus hypervariable region 1 in quasispecies analysis. Biomed Res Int. 2013, 2013: 626083-
Zagordi O, Bhattacharya A, Eriksson N, Beerenwinkel N: ShoRAH: estimating the genetic diversity of a mixed sample from next-generation sequencing data. BMC Bioinformatics. 2011, 12: 119-10.1186/1471-2105-12-119.
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011, 28 (10): 2731-2739. 10.1093/molbev/msr121.
Tamura K, Nei M: Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Mol Biol Evol. 1993, 10 (3): 512-526.
Librado P, Rozas J: DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009, 25 (11): 1451-1452. 10.1093/bioinformatics/btp187.
Abbate I, Cappiello G, Lo Iacono O, Longo R, Ferraro D, Antonucci G, Di Marco V, Di Stefano R, Craxi A, Solmone MC, Spanò A, Ippolito G, Capobianchi MR: Heterogeneity of HVR-1 quasispecies is predictive of early but not sustained virological response in genotype 1b-infected patients undergoing combined treatment with PEG- or STD-IFN plus RBV. J Biol Regul Homeost Agents. 2003, 17 (2): 162-165.
Abbate I, Lo Iacono O, Di Stefano R, Cappiello G, Girardi E, Longo R, Ferraro D, Antonucci G, Di Marco V, Solmone M, Craxì A, Ippolito G, Capobianchi MR: HVR-1 quasispecies modifications occur early and are correlated to initial but not sustained response in HCV-infected patients treated with pegylated- or standard-interferon and ribavirin. J Hepatol. 2004, 40 (5): 831-836. 10.1016/j.jhep.2004.01.019.
Morishima C, Polyak SJ, Ray R, Doherty MC, Di Bisceglie AM, Malet PF, Bonkovsky HL, Sullivan DG, Gretch DR, Rothman AL, Koziel MJ, Lindsay KL, Hepatitis C Antiviral Long-Term Treatment Against Cirrhosis Trial Group: Hepatitis C virus-specific immune responses and quasi-species variability at baseline are associated with nonresponse to antiviral therapy during advanced hepatitis C. J Infect Dis. 2006, 193 (7): 931-940. 10.1086/500952.
Saludes V, Bracho MA, Valero O, Ardevol M, Planas R, Gonzalez-Candelas F, Ausina V, Martro E: Baseline prediction of combination therapy outcome in hepatitis C virus 1b infected patients by discriminant analysis using viral and host factors. PLoS ONE. 2010, 5 (11): e14132-10.1371/journal.pone.0014132.
Pawlotsky JM, Pellerin M, Bouvier M, Roudot-Thoraval F, Germanidis G, Bastie A, Darthuy F, Remire J, Soussy CJ, Dhumeaux D: Genetic complexity of the hypervariable region 1 (HVR1) of hepatitis C virus (HCV): influence on the characteristics of the infection and responses to interferon alfa therapy in patients with chronic hepatitis C. J Med Virol. 1998, 54 (4): 256-264. 10.1002/(SICI)1096-9071(199804)54:4<256::AID-JMV4>3.0.CO;2-3.
Lopez-Labrador FX, Ampurdanes S, Gimenez-Barcons M, Guilera M, Costa J, Jimenez de Anta MT, Sanchez-Tapias JM, Rodes J, Saiz JC: Relationship of the genomic complexity of hepatitis C virus with liver disease severity and response to interferon in patients with chronic HCV genotype 1b infection [correction of interferon]. Hepatology. 1999, 29 (3): 897-903. 10.1002/hep.510290306.
Cuevas JM, Torres-Puente M, Jimenez-Hernandez N, Bracho MA, Garcia-Robles I, Wrobel B, Carnicer F, del Olmo J, Ortega E, Moya A, González-Candelas F: Genetic variability of hepatitis C virus before and after combined therapy of interferon plus ribavirin. PLoS ONE. 2008, 3 (8): e3058-10.1371/journal.pone.0003058.
Torres-Puente M, Cuevas JM, Jimenez-Hernandez N, Bracho MA, Garcia-Robles I, Wrobel B, Carnicer F, Del Olmo J, Ortega E, Moya A, González-Candelas F: Genetic variability in hepatitis C virus and its role in antiviral treatment response. J Viral Hepat. 2008, 15 (3): 188-199.
Gilles A, Meglecz E, Pech N, Ferreira S, Malausa T, Martin JF: Accuracy and quality assessment of 454 GS-FLX Titanium pyrosequencing. BMC Genomics. 2011, 12: 245-10.1186/1471-2164-12-245.
Beerenwinkel N, Zagordi O: Ultra-deep sequencing for the analysis of viral populations. Curr Opin Virol. 2011, 1 (5): 413-418. 10.1016/j.coviro.2011.07.008.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2334/14/389/prepub
This work was supported by grants NN401 6467 40 and 1M24/PM12/12 from The Polish National Science Centre.
The authors declare that they have no competing interest.
KCC participated in the design of the study and its coordination, amplified HVR1 sequences, carried out the molecular genetic studies deep sequencing of HVR, calculated genetic heterogeneity parameters, performed the statistical analysis and drafted the manuscript. OZ made reconstruction of HVR1 populations and correction of sequencing error, calculated genetic heterogeneity parameters, prepared figures and helped to draft the manuscript. KP prepared HVR1 amplicons for sequencing, made sequence alignments and calculated genetic heterogeneity parameters. TL participated in the design of the study and helped to draft the manuscript. KM calculated genetic heterogeneity parameters. IBO prepared the database of results and tables. AP made RNA isolation from samples, collected and analyzed clinical and virological data of patients. RP participated in the design of the study, interpretation of results and helped to draft the manuscript. HB participated in the design of the study, treated patients and provided information about the study, obtained informed consents and collected clinical and virological data of patients. AH participated in the design of the study, provided information about the study to the patients, collected clinical and virological data of patients, obtained informed consents. MR participated in the design of the study and helped to draft the manuscript. All authors read and approved the final manuscript.