Skip to main content

Evidence of recombination in quasispecies populations of a Hepatitis C Virus patient undergoing anti-viral therapy

Abstract

Background/Aim

Hepatitis C virus (HCV) has been the subject of intense research and clinical investigation as its major role in human disease has emerged. HCV circulates in vivo as a complex population of different but closely related viral variants, commonly referred to as a quasispecies. The extent to which recombination plays a role in the evolution of HCV quasispecies when patients are undergoing anti-viral therapy is currently unknown. In order to gain insight into these matters, we have performed a phylogenetic analysis of HCV quasispecies populations from six patients undergoing anti-viral therapy.

Methods

Putative recombinant sequences were identified with the use of SimPlot program. Recombination events were confirmed by bootscaning, using putative recombinant sequence as a query. Statistical support for the presence of a recombination event was done by the use of LARD program.

Results

A crossing-over event in the NS5A gene in a HCV strain recovered after four weeks of treatment was identified in quasispecies from a patient with sustained response. Putative parental-like strains were identified as strains circulating in previous weeks on the same patient.

Conclusion

Only one recombinant strain was detected in all patient quasispecies populations studied. The recombination break-point is situated on the PKR-binding region of NS5A. Although recombination may not appeared to be extensive in NS5A genes of HCV quasispecies populations of patients undergoing antiviral therapy, this possibility should be taken into account as a mechanism of genetic variation for HCV.

Background

Hepatitis C virus (HCV) is estimated to infect 170 million people worldwide and creates a huge disease burden from chronic, progressive liver disease [1] HCV has become a major cause of liver cancer and one of the commonest indications of liver transplantation [2, 3].

HCV has been classified in the family Flaviviridae, although it differs from other members of the family in many details of its genome organization [1]. Like most RNA viruses, HCV circulates in vivo as a complex population of different but closely related viral variants, commonly referred to as a quasispecies [47].

HCV is an enveloped virus with an RNA genome of approximately 9400 bp in length [8, 9]. Comparison of nucleotide sequences of variants recovered from different individuals and geographical regions has revealed the existence of at least six major genetic groups [1, 1012]. Each of the six major genetic groups of HCV contains a series of more closely related sub-types [1].

Interferon monotherapy provided the first hope for patients with chronic hepatitis C that the virus could be permanently eradicated. An important development in treating this disease was the recognition that the effects of interferon could be greatly enhanced by combining it with ribavirin, a nucleoside analogue. This combination regimen essentially doubled the sustained virological response rates seen with interferon alone. Recently, pegylated forms of interferon have been developed, and when that forms of interferon are used in combination with ribavirin, it demonstrates even better efficacy. For that reason, peginterferon alfa-2a and peginterferon alfa-2b are the latest innovations for the treatment of chronic hepatitis C [3, 13, 14].

Recombination plays a significant role in the evolution of RNA viruses by creating genetic variation. For example, the frequent recovery of poliovirus that results from recombination has the potential to produce "escape mutants" in nature as well as in experiments [15]. Recombination has also been detected in other RNA viruses for which multivalent vaccines are in use or in trials [16, 17].

Recently, a natural intergenotypic recombinant (2 k/1b) of HCV has been identified in Saint Petersburg (Russia) [18, 19]. Phylogenetic analyses of HCV strains circulating in Peru, demonstrated the existence of natural intra-genotypic HCV recombinant strains (1a/1b) circulating in the Peruvian population [20]. In these cases, the recombination events have taken place in the non-structural region of the HCV genome. Recombination break-points in HCV structural capsid genomic region has been recently identified [21].

Given the implications that recombination has for RNA virus evolution [16], it is clearly important to determine the extent to which recombination plays a role in the evolution of HCV quasispecies populations in vivo, when patients are undergoing anti-viral therapy.

Results

To gain insight into possible recombination events, a phylogenetic profile analysis was carried out using HCV NS5A sequences from HCV quasispecies populations obtained by Puig-Basgoiti et al. [22] from patients undergoing anti-viral therapy. Sequences were obtained by the use of the LANL database [23] (for patients, strains accession numbers, quasispecies obtained at different time points during therapy and therapy outcome, see Table 1). Phylogenetic profile analysis was done by the use of the SimPlot program [24]. Interesting, when the analysis was carried out for strain AY378694 (obtained on week 4 in patient No. 7, see Table 1), a recombination point was detected at position 286 of the NS5A sequence alignment and two putative parental-like strains (AY378615 and AY378641, obtained on weeks 0 and 2 from the same patient, respectively) were identified (see Fig. 1 and Table 1).

Table 1 Quasispecies population in HCV patients undergoing anti-viral therapya
Figure 1
figure 1

Phylogenetic profiles of HCV sequences. Results from SimPlot analysis are shown. The query sequence (AY378694) is indicated on the upper part of the figure. Sequences to be compared with the query sequence are indicated on the right side of the figure. When comparisons were done, SimPlot generates a similarity plot using the Kimura-two parameter distance model in a sliding window of 200 nucleotides, moving 20 nucleotides between plots. The y-axis gives the percentage of identity found. Comparison of HCV strain AY378694 with strains AY378615, AY378641 and AY378635 is shown. The red vertical line shows the recombination point at position 286. Red numbers on the bottom part of the figure denote the number of informative sites that support clustering of the query sequence with the respective strains indicated in red on the bottom left side of the figure.

In order to confirm these results, the same sequences were used for a bootscanning study [25]. The basic principle of bootscanning is that mosaicism is suggested when one observes high levels of phylogenetic relatedness between a query sequence and more than one reference sequence [25].When the putative recombinant strain identified in the previous analysis (AY378694) is used as a query, this is observed for this strain and the two putative parental-like strains previously detected (see Fig. 2). The same recombination break-point position is observed in the bootscanning analysis (see Figs. 1 and 2), confirming a recombination break-point at position 286 of the NS5A of strain AY378694.

Figure 2
figure 2

Bootscanning of HCV sequences. Query sequence (AY378694) is shown on the upper part of the figure. Sequences to be compared with the query sequence are indicated on the right side of the figure. When comparisons were done, SimPlot generates a graph of percentage of permutated trees obtained using a sliding window of 200 nucleotides, moving 20 nucleotides at a time. The y-axis gives the percentage of permutated trees. This approach permits to observe levels of phylogenetic relatedness between a query sequence and a reference sequence in different genomic regions. The rest same as Fig. 1A.

To assess whether the recombination model we obtained gave a significantly better fit to the data than the null hypothesis of no recombination, we used LARD [26]. The results of these studies are shown in Fig. 3.

Figure 3
figure 3

Distribution of the likelihood ratios expected by chance. The distribution of likelihood ratios for the null hypothesis (i.e. no recombination) is shown. The y-axis shows the number of simulations. Likelihood ratios are shown at the bottom of the figure. The arrow shows the likelihood ratio obtained for the real dataset for the putative recombinant strain.

As it can be seen in the figure, simulations of sequence evolution under the null hypothesis (i.e., no recombination) gave strong statistical support for the alternative hypothesis of recombination (P < 0.004, Fig. 3).

Discussion

The extent to which recombination plays a role in the evolution of HCV quasispecies when patients are undergoing anti-viral therapy is currently unknown. In order to get insight into this issue, we performed phylogenetic studies on the NS5A gene of HCV quasispecies populations from patients undergoing antiviral therapy. We have selected NS5A since previous work by Enomoto et al. [27] suggested that the genetic heterogeneity of a specific domain of HCV NS5A, termed IFN sensitivity-determining region (ISDR), was related closely to response in Japanese patients. Although this issue continues to be controversial [22], analysis of the published information supports the hypothesis that a relationship exists between NS5A and response to therapy [28, 29].

The results of these studies reveal that recombination can not be denied as an evolutionary mechanism for generating diversity in HCV in vivo, in patients undergoing anti-viral therapy (see Figs. 1 and 2). Recombination does not seems to play an extensive roll in the evolution of HCV quasispecies populations, at least by the study of NS5A genes, since only one recombinant isolate was observed among all HCV quasispecies populations studied. On the other hand, the true frequency of recombination may be underestimated, since a statistical significant number of differences among recombinant and putative parental-like strains are needed in order to achieve detection by current methods applied to detect recombination events. Recombination may serve two opposite purposes: exploration of a new combination of genomic region from different origins or rescuing of viable genomes from debilitated parental genomes [30]. Interestingly, the recombination break-point identified in recombinant strain AY378694 is situated on the PKR-binding region of NS5A ISDR [27, 31] (see Fig. 4). Recent evidence suggests that HCV NS5A protein can repress PKR function in vivo, possibly allowing HCV to escape the antiviral effects of interferon [1, 3, 32, 33]. An analysis of NS5A translated sequences of recombinant and putative parental-like virus suggests the possibility that recombinant isolate AY378694 may have acquired amino acids known to be present in HCV strains resistant to interferon treatment (see Fig. 4), although more work is needed in order to test this hypothesis, since the results found for the ISDR still remains a controversial issue [22]. For that reason, although recombination may not appeared to be extensive in NS5A genes of HCV quasispecies populations of patients undergoing antiviral therapy, this possibility should be taken into account as a mechanism of genetic variation for HCV.

Figure 4
figure 4

Alignment of amino acid sequences of the PKR-binding region and ISDR of recombinant and putative parental-like strains. Strains are indicated by accession number in the left side of the figure. Identity to strain recombinant strain AY378694 in putative parental-like strains AY378615 and AY 378641 is shown by a dash. PKR-binding region and ISDR sequences are shown in bold. Recombination break-point is indicated by an arrow. Amino acids known to be present in interferon resistant viruses are indicated in red in recombinant strain AY378694 (see refs. [31, 35]).

Conclusion

Only one recombinant strain was detected in all patient quasispecies populations studied. The recombination break-point identified in strain AY378694 is situated on the PKR-binding region of NS5A. The results of these studies reveal that recombination events can be observed in patients undergoing anti-viral therapy. Recombination can not be denied as a mechanism of genetic variation for HCV.

Methods

Patients and strains

Patients and strains refereed in these studies belong to a recent study done by Puig-Basagoiti et al. [22]. All patients had genotype 1b infection. Treatment consisted of the administration of IFN-α-2b and ribavirin (see ref. [22]).

Sequences

NS5A sequences from the study done by Puig-Basagoity et al. [22] were obtained from the HCV LANL database [23] (for patient identification, quasispecies accession numbers obtained at each week of therapy and therapy outcome, see Table 1). Sequences were aligned using the CLUSTAL W program [34].

Recombination analysis

Putative recombinant sequences were identified with the SimPlot program [24]. This program is based on a sliding window method and constitutes a way of graphically displaying the coherence of the sequence relationship over the entire length of a set of aligned homologous sequences. The window width and the step size were set to 200 bp and 20 bp, respectively.

Bootscanning analysis

The results obtained in the recombination analysis were confirmed using a bootscanning analysis [25]. The window width and the step size were set to 200 bp and 20 pb, respectively.

LARD analysis

To assess whether the recombination model we obtained gave a significantly better fit to the data than the null hypothesis of no recombination, we used LARD [26]. Briefly, for every possible breakpoint, the sequence alignment was divided into two independent regions for which the branch lengths of a tree of the putative recombinant and its two parent sequences were optimised. The two results (likelihoods) obtained by using the separate regions were then combined to give a likelihood score for that breakpoint position and the breakpoint position that yielded the highest likelihood then was compared, by using a likelihood ratio test, to the likelihood obtained from the same data under a model that permitted no recombination. The likelihood ratio obtained by using the real data were evaluated for significance against a null distribution of likelihood ratios produced by using Monte Carlo simulation of sequences generated without recombination. Sequences were simulated 1,000 times by using the maximum likelihood model parameters and sequence lengths from the real data using Seq-Gen [35].

References

  1. Simmonds P: Genetic diversity and evolution of hepatitis C virus 15 years on. J Gen Virol 2004, 85: 3173-3188. 10.1099/vir.0.80401-0

    Article  CAS  PubMed  Google Scholar 

  2. Hoofnagle JH: Course and outcome of hepatitis C. Hepatology 2003, 36: S21-S29. 10.1002/hep.1840360704

    Article  Google Scholar 

  3. Pawlotsky JM: The nature of interferon-alfa resistance in hepatitis C virus infection. Curr Opin Infect Dis 2003, 16: 587-592. 10.1097/00001432-200312000-00012

    Article  CAS  PubMed  Google Scholar 

  4. Chambers TJ, Fan X, Droll DA, Hembrador E, Slater T, Nickells MW, Dustin LB, Dibisceglie AM: Quasispecies heterogeneity within the E1/E2 region as a pretreatment variable during pegylated interferon therapy of chronic hepatitis C virus infection. J Virol 2005, 79: 3071-3083. 10.1128/JVI.79.5.3071-3083.2005

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  5. Laskus T, Wilkinson J, Gallegos-Orozco JF, Radkowski M, Adair DM, Nowicki M, Operskalsi 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: 764-776. 10.1053/j.gastro.2004.06.005

    Article  CAS  PubMed  Google Scholar 

  6. Feliu A, Gay E, Garcia-Retortillo M, Saiz JC, Foms X: Evolution of hepatitis C virus quasispecies immediately following liver transplantation. Liver Transpl 2004, 10: 1131-1139. 10.1002/lt.20206

    Article  PubMed  Google Scholar 

  7. Martell M, Esteban JL, 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: 3225-3229.

    PubMed Central  CAS  PubMed  Google Scholar 

  8. You S, Stump DD, Branch AD, Rice CM: A cis -actingreplication element in the sequence encoding the NS5B RNA-dependent RNA polymerase is required for hepatitis C virus RNA replication. J Virol 2004, 78: 1352-1366. 10.1128/JVI.78.3.1352-1366.2004

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  9. Pestova TV, Shatsky IN, Fletcher SP, Jackson RJ, Hellen CUT: A prokaryotic-like mode of cytoplasmic eukaryotic ribosome binging to the initiation codon during internal translation initiation of hepatitis C and classical swine fever virus RNAs. Genes Dev 1998, 12: 6783.

    Article  Google Scholar 

  10. Simmonds P, Bukh J, Combet C, Deleage G, Enomoto N, Feinstone S, Halfon P, Inchauspe G, Kuiken C, Maertens G, Mizokami M, Murphy DG, Okamoto H, Pawlotsky JM, Penin F, Sablon E, Shin IT, Stuyver LJ, Thiel HJ, Viazov S, Weiner AJ, Widell A: Consensus proposals for a unified system of nomenclature of hepatitis C virus genotypes. Hematology 2005, 42: 962-973.

    CAS  Google Scholar 

  11. Simmonds P: The origin and evolution of hepatitis viruses in humans. J Gen Virol 2001, 82: 693-712.

    Article  CAS  PubMed  Google Scholar 

  12. Simmonds P, Holmes EC, Cha TA, Chan SW, McOmish F, Irvine B, Beall E, Yap PL, Kolberg J, Urdea MS: Classification of hepatitis C virus into six major genotypes and a series of subtypes by phylogenetic analysis of the NS-5 region. J Gen Virol 1993, 74: 2391-2399.

    Article  CAS  PubMed  Google Scholar 

  13. Zeuzem S: Heterogenous virologic response rates to interferon-based therapy in patients with chronic hepatitis C: who responds less well? Ann Intern Med 2004, 140: 370-381.

    CAS  PubMed  Google Scholar 

  14. Fried MW, Hadziyannis SJ: Treatment of chronic hepatitis C infection with peginterferons plus ribavirin. Semin Liver Dis 2004, 24: 47-54. 10.1055/s-2004-832928

    Article  CAS  PubMed  Google Scholar 

  15. Kew OM, Nottay BK: Evolution of the oral polio vaccine strains in humans occurs by both mutation and intra-molecular recombination. In Modern Approaches to Vaccines. Edited by: Chanock RM, Lerner RA. Cold Spring Harbor, NY: Cold Spring Harbor Laboratory; 1984:357-363.

    Google Scholar 

  16. Worobey M, Holmes EC: Evolutionary aspects of recombination in RNA viruses. J Gen Virol 1999, 80: 2535-2543.

    Article  CAS  PubMed  Google Scholar 

  17. Suzuki Y, Gojobori T, Nakagomi O: Intragenic recombinations in rotaviruses. FEBS Lett 1998, 427: 183-187. 10.1016/S0014-5793(98)00415-3

    Article  CAS  PubMed  Google Scholar 

  18. Kalinina O, Norder H, Magnius O: Full-length open reading frame of a recombinant hepatitis C virus strain from St. Petersburg: proposed mechanism of its formation. J Gen Virol 2004, 85: 1853-1857. 10.1099/vir.0.79984-0

    Article  CAS  PubMed  Google Scholar 

  19. Kalinina O, Norder H, Mukomolov S, Magnius LO: A natural intergenotypic recombinant of hepatitis C virus identified in St. Petersburg. J Virol 2002, 76: 4034-4043. 10.1128/JVI.76.8.4034-4043.2002

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  20. Colina R, Casane D, Vasquez S, Garcia-Aguirre L, Chunga A, Romero H, Khan B, Cristina J: Evidence of intratypic recombination in natural populations of hepatitis C virus. J Gen Virol 2004, 85: 31-37. 10.1099/vir.0.19472-0

    Article  CAS  PubMed  Google Scholar 

  21. Cristina J, Colina R: Evidence of structural genomic region recombination in Hepatitis C virus. Virology J 2006, 3: 53-56. 10.1186/1743-422X-3-53

    Article  Google Scholar 

  22. Puig-Basagoiti F, Forns X, Furcic I, Ampurdanes S, Gimenez-Barcons M, Franco S, Sanchez-Tapias JM, Saiz JC: Dynamics of hepatitis C virus NS5A quasispecies during interferon and ribavirin therapy in responder and non-responder patients with genotype 1b chonic hepatitis C. J Gen Virol 2005, 86: 1067-1075. 10.1099/vir.0.80526-0

    Article  CAS  PubMed  Google Scholar 

  23. Kuiken C, Yusim K, Boykin L, Richardon R: The HCV Sequence Database. Bioinformatics 2005, 21: 379-384. 10.1093/bioinformatics/bth485

    Article  CAS  PubMed  Google Scholar 

  24. Lole KS, Bollinger RC, Parnjape RS, Gadkari D, Kulkarni SS: Full-length human immunodeficiency virus type 1 genomes from subtype C-infected seroconverters in India, with evidence of intersubtype recombination. J Virol 1999, 73: 152-160.

    PubMed Central  CAS  PubMed  Google Scholar 

  25. Salminen MO, Carr JK, Burke DS, McCutchan FE: Identification of breakpoints in intergenotypic recombinants of HIV type 1 by bootscanning. AIDS Res Hum Retroviruses 1995, 11: 1423-1425.

    Article  CAS  PubMed  Google Scholar 

  26. Holmes EC, Worobey M, Rambaut A: Phylogenetic evidence for recombination in dengue virus. Mol Biol Evol 1999, 16: 405-406.

    Article  CAS  PubMed  Google Scholar 

  27. Enomoto N, Sakuma I, Asahina Y, Kurosaki M, Murakami T, Yamamoto C, Izumi N, Marumo F, Sato C: Comparison of full-length sequences of interferon-sensitive and resistant hepatitis C virus 1b. Sensitivity to interferon is conferred by amino acid substitutions in the NS5A gene. J Clin Investig 1995, 96: 224-230.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  28. Gimenez-Barcons M, Franco S, Suarez Y, Forns X, Ampurdanes S, Puig-Basagoiti F, Sanchez-Fueyo A, Barrera JM, Llovet JM, Bruix J, Sanchez-Tapias JM, Rodes J, Saiz JC: High amino acid variability within the NS5A of hepatitis C virus (HCV) is associated with hepatocellular carcinoma in patients with HCV-1b-related cirrhosis. Hepatology 2001, 34: 158-167. 10.1053/jhep.2001.25512

    Article  CAS  PubMed  Google Scholar 

  29. Witherell GW, Beineke P: Statistical analysis of combined substitutions in non-structural 5A region of hepatitis C virus and interferon response. J Med Virol 2001, 63: 8-16. 10.1002/1096-9071(200101)63:1<8::AID-JMV1001>3.0.CO;2-K

    Article  CAS  PubMed  Google Scholar 

  30. Domingo E, Holland JJ: RNA virus mutations and fitness for survival. Annu Rev Microbiol 1997, 51: 151-178. 10.1146/annurev.micro.51.1.151

    Article  CAS  PubMed  Google Scholar 

  31. Gale MJ, Korth MJ, Tang NM, Tan SL, Hopkins DA, Dever TE, Polyak SJ, Gretch DR, Katze MG: Evidence that hepatitis C virus resistance to interferon is mediated through repression of the PKR protein kinase by the nonstructural 5A protein. Virology 1997, 230: 217-227. 10.1006/viro.1997.8493

    Article  CAS  PubMed  Google Scholar 

  32. Gale MJ, Blakely CM, Kwieciszewski B, Tan SL, Dossett M, Tang NM, Korth MJ, Polyak SJ, Gretch DR, Katze MG: Control of PKR protein kinase by Hepatitis C Virus nonstructural 5A protein: molecular mechanisms of kinase regulation. Mol Cel Biol 1998, 18: 5208-5218.

    Article  CAS  Google Scholar 

  33. Pawlotsky JM, Germanidis G, Neumann AU, Pellerin M, Frainais PO, Dhumeaux D: Interferon Resistance of Hepatitis C Virus genotype 1b: relationship to nonstructural 5A gene quasispecies mutations. J Virol 1998, 72: 2795-2805.

    PubMed Central  CAS  PubMed  Google Scholar 

  34. Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acid Res 1994, 22: 4673-4680.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  35. Rambaut A, Grassly NC: Seq-Gen: an application for the Monte Carlo simulation of DNA sequence evolution along phylogenetic trees. Comp Appl Biosci 1997, 13: 235-238.

    CAS  PubMed  Google Scholar 

  36. Enomoto N, Sakuma I, Asahina Y, Kurosaki M, Murakami T, Yamamoto C, Ogura Y, Izumi N, Maruno F, Sato C: Mutations in the non-structural protein 5A gene and response to interferon in patients with chronic hepatitis C virus 1b infection. N Engl J Med 1996, 344: 77-81. 10.1056/NEJM199601113340203

    Article  Google Scholar 

  37. Pawlotsky JM, Germanidis G, Neumann AU, Pellerin M, Frainais PO, Dhumeaux D: Interferon Resistance of Hepatitis C Virus genotype 1b: relationship to nonstructural 5A gene quasispecies mutations. J Virol 1998, 72: 2795-2805.

    PubMed Central  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

This work was supported by ICGEB, PAHO, and RELAB through Project CRP.LA/URU03-032. International Atomic Energy Agency through project ARCAL 6050 and from DINACYT (Uruguay) through Project Fondo Clemente Estable 05/039 and Project No. 8006.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Juan Cristina.

Additional information

Competing interests

The author(s) declare that they have no competing interests.

Authors' contributions

JC and PM conceived the study. JC, PM and LL designed the analysis. JC and PM performed the SimPlot and Bootscanning analysis. DC performed the LARD analysis and contributed to the discussion of all results found in this work. JC wrote the paper.

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Moreno, M.P., Casane, D., López, L. et al. Evidence of recombination in quasispecies populations of a Hepatitis C Virus patient undergoing anti-viral therapy. Virol J 3, 87 (2006). https://doi.org/10.1186/1743-422X-3-87

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1743-422X-3-87

Keywords