| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Previous Article | Next Article ![]()
Journal of Virology, July 2007, p. 7732-7741, Vol. 81, No. 14
0022-538X/07/$08.00+0 doi:10.1128/JVI.00382-07
Copyright © 2007, American Society for Microbiology. All Rights Reserved.

French National Reference Center for Viral Hepatitis B, C and delta, Department of Virology, Hôpital Henri Mondor, Université Paris 12, Créteil, France,1 INSERM U841, Créteil, France,2 Centro de Astrobiología (INTA-CSIC), Torrejón de Ardoz, Madrid, Spain,3 Department of Hepatology and Gastroenterology, Hôpital Henri Mondor, Université Paris 12, Créteil, France4
Received 22 February 2007/ Accepted 2 May 2007
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
) and ribavirin is the current standard treatment for chronic HCV infection. This treatment yields a sustained virological response (SVR), defined by the lack of detectable HCV RNA in serum 24 weeks after the end of therapy, in 40 to 50% of patients infected by HCV genotype 1 and in 70% to 80% of patients infected by genotype 2 or 3 (14, 17, 27). The vast majority of patients who have an SVR are cured of the infection.
Ribavirin (1-ß-D-ribofuranosyl-1,2,4-triazole-3-carboxamide) is a guanosine analogue with a broad spectrum of activity against DNA and RNA viruses (45). Its use in the treatment of chronic HCV infection has been essentially empirical (36). Indeed, the addition of ribavirin to standard or pegylated IFN-
therapy significantly increases the SVR rate over that with IFN monotherapy (14, 17, 27, 30, 39) and ribavirin can exert this antiviral effect through several possible mechanisms. HCV clearance is biphasic in responders to IFN-
therapy (32, 33). The first rapid phase, which occurs on days 1 and 2, reflects direct inhibition of viral replication by IFN-induced antiviral effectors. Mathematical modeling suggests that the second, slower phase, which starts during the first week of administration, reflects gradual clearance of infected cells in a context of potent inhibition of viral replication (32, 33). Two recent studies by our group suggest that ribavirin impacts both phases (4, 37). Indeed, ribavirin monotherapy exerts significant, moderate, early (days 2 and 3), and transient inhibition of viral replication for approximately half of the patients, and this effect appears to be additive with that of IFN during the first days of therapy (37). This inhibition correlates with the concentration and elimination half-life of ribavirin in peripheral blood (37). We have also recently shown in a prospective randomized trial involving patients infected by HCV genotype 1 that the principal effect of ribavirin is to prevent breakthroughs during and relapses after therapy for individuals who initially respond to the pegylated IFN-
-ribavirin combination (4). This suggests that ribavirin shortens the half-life of infected cells during IFN administration, allowing their complete elimination during the standard 48-week treatment period in the majority of patients who initially respond to combination therapy.
The mechanisms by which ribavirin modestly inhibits HCV replication and prevents relapses during IFN-ribavirin combination therapy are totally unknown. In vivo, ribavirin is converted into ribavirin 5'-monophosphate (RMP), ribavirin 5'-diphosphate, and ribavirin 5'-triphosphate (RTP) by adenosine kinase and possibly cytosolic 5'-nucleotidase II (51). RTP is the principal intracellular metabolite of ribavirin in mammals, and its intracellular concentration appears to vary according to the cell type (15, 35). A number of putative mechanisms have been proposed to explain the efficacy of ribavirin against HCV infection (13, 22).
(i) Direct inhibition of HCV RNA-dependent RNA polymerase: RTP can be incorporated as a substrate by viral polymerases and can thereby inhibit RNA elongation. However, RTP has been shown to be a weak inhibitor of a variety of viral polymerases, including HCV RNA-dependent RNA polymerase (26, 50). Given the modest antiviral effect of ribavirin monotherapy in vivo (37), ribavirin's direct antiviral properties are unlikely to play a major role in HCV elimination during combination therapy.
(ii) Inhibition of IMPDH: ribavirin 5'-monophosphate, which is a minor intracellular metabolite, competitively inhibits both isoforms of IMP dehydrogenase (IMPDH). IMPDH inhibition leads to the depletion of intracellular pools of GTP and dGTP, which are necessary for viral replication (28). IMPDH inhibition has been reported to be the main mechanism whereby ribavirin acts against yellow fever virus and human parainfluenza virus 3 (23, 24). In contrast, recent results with potent specific IMPDH inhibitors, used alone or in combination with ribavirin or IFN-
for patients with chronic HCV infection, suggest that IMPDH inhibition is not the mechanism of ribavirin action in chronic hepatitis C (18, 31).
(iii) Immunomodulatory properties: several reports suggest that at clinically relevant concentrations, ribavirin acts as an immunomodulator, enhancing the production of Th1 cytokines (interleukin-2 and gamma interferon) while suppressing Th2 cytokine production both in murine models and in humans; it is also suggested that these properties might enhance the activity of IFN-
(6, 19, 29, 34, 46). However, the L-eniantomer of ribavirin, levovirin, which has the same immunomodulatory properties as ribavirin but no activity against DNA and RNA viruses in vitro (because it cannot be converted into 5'-phosphate derivatives (41, 47), does not share ribavirin's clinical properties when administered to HCV-infected patients in combination with pegylated IFN-
(38).
(iv) Mutagenic properties leading to "error catastrophe": the most popular hypothesis forwarded to explain the role of ribavirin in hepatitis C therapy is its putative mutagenic potential, which would lead to "error catastrophe," i.e., disorganization of the mutant distribution of the quasipecies and the generation of nonviable viral populations. This concept was recently introduced by Crotty et al., based on the in vitro poliovirus model. These authors reported that ribavirin induced a high frequency of genomic mutations that generated viral progeny with reduced infectivity (7, 8). A similar effect of high concentrations of ribavirin has been found in cell cultures infected by RNA viruses belonging to several families, including foot-and-mouth disease virus, Hantaan virus, West Nile virus, and GB virus B, a member of the Flaviviridae family that is closely related to HCV (1, 9, 16, 20, 44). In addition, ribavirin has been reported to exert an antiviral effect related to error-prone replication in HCV subgenomic replicons, both in Huh7 cell lines and in a full-length binary HCV replication system (5, 20, 21, 48, 53). However, error-prone replication is not the mechanism by which ribavirin acts against arenavirus lymphocytic choriomeningitidis virus or the yellow fever virus, another member of the Flaviviridae family that is related to HCV (24, 42). Conflicting results have been obtained with HCV-infected patients receiving ribavirin (3, 25, 43, 52).
A recent study based on mathematical modeling of viral decay in 17 patients receiving the combination of standard IFN-
and ribavirin suggested that ribavirin acts mainly by reducing the infectivity of HCV virions produced during therapy in a dose-dependent manner and that this effect may appear only when viral replication is profoundly reduced by IFN-
, thereby explaining the poor antiviral effect of ribavirin monotherapy (10). Although mathematical modeling can only raise hypotheses, the authors suggested that such a reduction in virus infectivity would be in keeping with the "error catastrophe" hypothesis (10).
The aims of this study were to determine, by means of thorough quasispecies sequence analysis targeting two independent genomic regions, the following: (i) whether ribavirin acts as an RNA mutagenic agent during treatment of patients with chronic HCV genotype 1 infection and (ii) whether ribavirin's mutagenic effect (if any) occurs only in patients with low-level virus production due to potent inhibition by coadministered IFN.
| MATERIALS AND METHODS |
|---|
|
|
|---|
and/or ribavirin (37). The patients were eight men and three women, and their mean age was 38 ± 5 years (range, 31 to 52 years). All of the patients had anti-HCV antibodies in serum, detectable HCV RNA (>50 IU/ml), elevated alanine aminotransferase levels, and compensated chronic liver disease. Four patients (group A) received ribavirin (Rebetol; Schering-Plough Oncology, Kenilworth, NJ) monotherapy at a dose of 1.0 or 1.2 g/day, depending on body weight (below or above 75 kg, respectively). Four patients (group B) received the combination of standard IFN-
2b (Intron-A; Schering-Plough Oncology) at a dose of 3 million units (MU) three times per week, plus ribavirin, 1.0 or 1.2 g/day depending on body weight. The remaining three patients (group C) were treated with the combination of standard IFN-
2b 3 MU daily plus ribavirin 1.0 or 1.2 g/day, depending on body weight. All patients were treated with the relevant regimen for 24 weeks, as previously reported (37). Written informed consent was obtained from all patients, and the study was approved by the local ethics committee (Comité Consultatif de Protection des Personnes dans la Recherche Biomédicale) in accordance with the Helsinki Declaration.
Materials.
Blood samples were taken at four time points for sequence analyses: 14 days before the initiation of therapy (d14), at the initiation of therapy (d0), on day 14 of therapy (d14) and on day 28 of therapy (d28). At each time point and for each patient, extensive quasispecies sequence characterization (20 clones) of both the full-length nonstructural protein 3 (NS3) protease domain and the full-length NS5A protein was performed. The rate of variation of the consensus sequence, the frequency of mutation, the error generation rate, and the between-sample genetic distance were calculated during the three successive periods for each patient, i.e., between d14 and d0 (natural genetic drift in the absence of therapy, used as the control period for each patient), between d0 and d14 of therapy, and between d14 and d28 of therapy (when viral replication was profoundly inhibited by IFN-
in groups B and C). HCV RNA was quantified at d14, at baseline (d0), and on days 1, 2, 3, 14, and 28 of therapy. Blood samples for serum ribavirin assay were obtained at d0, hours 4, 8, 12, 18, 24, 30, 36, 42, 48, 62, 70, 86, and 94, and on days 14 and 28.
HCV RNA quantification. HCV RNA was quantified by using the SuperQuant assay (National Genetics Institute, Los Angeles, CA), as previously described (37). The detection limit of this assay is 30 IU/ml (1.5 log10 IU/ml).
Serum ribavirin assay. Serum ribavirin concentrations were determined by means of a high-performance liquid chromatography tandem mass spectrometric method at MDS Pharma Services (Saint-Laurent, Québec, Canada), as previously reported (40).
Generation of NS3 and NS5A quasispecies sequences. (i) RNA extraction. RNA was extracted from 200 µl of serum by using the High Pure viral RNA kit (Roche Molecular Biochemicals, Mannheim, Germany), according to the manufacturer's instructions. The RNA pellet was eluted with 50 µl of diethyl pyrocarbonate-treated water and stored at 70°C until analysis.
(ii) cDNA synthesis. For cDNA synthesis, 5 µl of viral RNA was used with 10 pmol of antisense primers 3EAS and R3EAS for NS3 and NS5A, respectively (Table 1). The reaction mixture (20 µl) contained 50 mM Tris-HCl, pH 8.3, 75 mM KCl, 3 mM MgCl2, 10 Mm deoxynucleoside triphosphates, 100 mM dithiothreitol, 40 U of RNase inhibitor (RNaseOUT; Invitrogen, Carlsbad, CA), and 200 U of Superscript III reverse transcriptase (Invitrogen). The reaction mixture was incubated for 1 h at 55°C, and the enzyme was then inactivated for 15 min at 70°C.
|
(iv) Purification of PCR products. The PCR products were purified by using Microcon-100 centrifugal filters (Millipore, Bedford, MA) according to the manufacturer's instructions. Amplicons were eluted in 20 µl of sterile water and stored at 20°C until analysis.
(v) Cloning and sequencing. Purified products were cloned into the pCR4-TOPO vector, using the TOPO TA cloning kit (Invitrogen) according to the manufacturer's protocol. Recombinant plasmids were transformed in Escherichia coli TOP10 competent cells, and transformants were grown on ampicillin plates. Cloned DNAs were PCR amplified with the M13 plasmid universal primers and Pfu DNA polymerase (Qbiogene, Irvine, CA). The amplified clones were purified by means of the Wizard SV96 PCR Clean-Up system (Promega, Madison, WI) according to the manufacturer's protocol. Twenty clones per time point per patient were sequenced with the Big-Dye Terminator Cycle v3.1 sequencing kit on the ABI 3100 sequencer (Applied Biosystems, Foster City, CA), according to the manufacturer's protocol. Both DNA strands were sequenced with the M13 universal primers.
NS3 and NS5A sequence analysis and statistical analysis. The generated sequences were aligned by using CLUSTAL W software (49). For each sample from each patient, the consensus nucleotide sequence of the quasispecies was deduced and its rate of variation was determined by counting the number of nucleotide substitutions at each 14-day interval, dividing it by the length of the sequenced fragment, and multiplying it by 26 to give the variation per site per year.
The error generation rate was calculated by measuring the number of nucleotide substitutions in all of the sequenced clones relative to the consensus sequence at the previous time point divided by the total number of sequenced nucleotides over each 14-day interval (d14 to d0, d0 to d14, and d14 to d28) for each patient. This number was multiplied by 26 to obtain the error generation rate, expressed as the number of substitutions per site per year.
The mutation frequency was determined as the total number of nucleotide substitutions (in the 20 clones analyzed at each time point) relative to the consensus sequence at the same time point, divided by the total number of sequenced nucleotides. The normalized Shannon entropy was calculated as
i(pi ln pi)/ln N, where pi represents the frequency of each sequence in the quasispecies and N is the total number of sequences analyzed. The between-sample genetic distances were calculated on the basis of distances between pairs of sequences generated at 14-day intervals.
Statistical calculations were carried out by using the descriptive statistics, linear regression, hypothesis test, and analysis-of-variance (ANOVA) packages for Mathematica 5.0 (Wolfram Research).
Nucleotide sequence accession numbers. Sequence data from this article have been deposited in GenBank under accession numbers AM400248 to AM400838 and AM401639 to AM402384.
| RESULTS |
|---|
|
|
|---|
three times per week and ribavirin (group B) (Fig. 1B), and those receiving a combination of IFN-
daily and ribavirin (group C) (Fig. 1C). High ribavirin concentrations in serum were achieved within a few hours after the first oral intake, but the levels fluctuated during the first days of administration. A plateau was reached by day 14 in all cases.
|
|
three times per week plus ribavirin) (Fig. 2B) and in group C (IFN-
daily plus ribavirin) (Fig. 2C), HCV RNA showed a typical biphasic decline. The initial, rapid decline was observed on day 1 after the first IFN-
dose. It ranged from 0.22 to 2.13 log10 HCV RNA IU/ml. It was followed by a reincrease in viral replication on day 2 in group B before the second IFN-
injection (Fig. 2B) but not for group C patients, who received daily IFN-
injections (Fig. 2C). A significant second slope of viral decline was observed for all patients in these two groups and was more pronounced in group C than in group B (Fig. 2B and C). In all cases and in both groups, the sample at week 28 was the last in which HCV RNA was detectable (i.e., HCV RNA was undetectable in the next sample, 4 weeks later).
Clonal analysis of HCV quasispecies.
Overall, 1,600 quasispecies variant sequences were generated, including 20 clones per patient of the full-length NS3 protease- and NS5A-encoding regions at 4 time points for the 8 patients in groups A and B and at 3 time points for the 3 patients in group C. On day 28, for the group C patients, who received daily IFN-
and ribavirin and had a more rapid HCV RNA decline, NS3 could not be PCR amplified in one case, NS5A in another case, and both in the remaining case, owing to low HCV RNA levels. In order to determine whether our sampling was biased by a limiting-dilution effect, especially in the samples with a very low HCV RNA content, which would have hidden a possible mutagenic effect of ribavirin, all templates were diluted 1/10 and the same PCR amplification procedures were applied. All samples gave visible PCR bands in both regions, except for two samples with a very low HCV RNA content at D28 for patients receiving the combination of IFN-
and ribavirin. This experiment ruled out a sampling bias.
To determine whether the analysis of 20 clones per time point was sufficient to evaluate total viral populations, we determined the mean, standard deviation (SD), standard error of the mean (SEM), and 95% confidence intervals (CI) for the mutation frequency in both genomic regions analyzed (Table 2). As shown in Table 2, the SEM values were on average 1 order of magnitude lower than the means, suggesting that the sample size was sufficient to ensure excellent accuracy.
|
|
|
.
Additional regression analyses focused on the average difference in mean mutation frequencies for each 14-day period, pooling the results for the 11 patients, and for each treatment group. In the NS3 region, the P values for the slopes of the straight lines best fitting the data were always higher than 0.05 (data not shown), confirming the lack of acceleration of mutation accumulation during ribavirin administration with or without IFN-
. The P values were also nonsignificant for the NS5A region.
In a global analysis based on ANOVA, three or four groups of values were compared. The difference in the mean mutation frequencies for the d14-to-d0 interval for samples from all patients was used as the control. It was compared with the differences in mean mutation frequencies in each treatment group for the d0-to-d14 interval on the one hand and for the d14-to-d28 interval on the other hand (group C could not be included in this last comparison because PCR amplification on d28 was possible with only one patient). This analysis showed that there was no significant difference in the mutation frequencies between the treatment groups or in either group relative to the control during either interval.
Error generation rates and between-sample genetic distances. In order to complete the study of the diversity of HCV quasispecies during treatment with ribavirin, two additional parameters, the error generation rate and the between-sample genetic distance, were analyzed.
The error generation rate was determined as the mutation frequency relative to the consensus sequence at the previous time point, with a correction factor to obtain the rate per year (Tables 5 and 6). The error generation rate did not increase with ribavirin administration for any of the patients (regression analysis). Likewise, ANOVA showed no significant difference in the error generation rates between the controls (d14 to d0) and the treated groups at either time interval (d0 to d14 or d14 to d28).
|
|
A similar global ANOVA analysis of mutation frequencies and error generation rates was performed for the between-sample genetic distances and showed no difference between the three treatment groups or with the controls, either between d0 and d14 or between d14 and d28 (data not shown).
Rate of variation of the consensus sequence. The rate of variation of the consensus sequence was also calculated in order to determine if treatment with ribavirin affects the fixation of mutations in the HCV genome, as would happen if the drug induced the selection of specific nucleotide substitutions. The results are shown in Tables 5 and 6 for the 11 patients. Regression analyses were performed with the differences in the rate of variation of the consensus sequence at 14-day intervals for samples from the 11 patients in the NS3 and NS5A coding regions. The P values for the slopes of the straight lines best fitted to the data indicated no sustained linear trend towards an increase in the fixation of mutations on ribavirin therapy (data not shown). The differences in the mean values at consecutive time points for patients belonging to a given treatment group were also not significant. A similar global ANOVA analysis of mutation frequencies was performed with this parameter and showed no significant difference between the three treatment groups or with the controls, either between d0 and d14 or between d14 and d28.
| DISCUSSION |
|---|
|
|
|---|
. The spontaneous error generation rates observed for our three groups of patients prior to therapy, reflecting natural genetic drift, were higher than the known error rate of the HCV RNA-dependent RNA polymerase (11, 12). This was probably due to the fact that they were measured over a short period of time (14 days) and finally expressed per year.
Ribavirin monotherapy did not increase the mutation frequency, the error generation rate, or the between-sample genetic distances during either the first or the second 14-day period with treatment. Likewise, the accumulation of nucleotide substitutions did not accelerate relative to that in the pretreatment period during combination therapy with ribavirin and IFN-
at a dose of 3 MU given either daily or three times per week. Finally, these parameters did not increase between days 14 and 28 of therapy, when HCV replication was profoundly inhibited by IFN administration (day 28 corresponded to the last sample containing detectable HCV RNA for patients receiving combination therapy). Our painstaking statistical analysis showed no common trend towards an acceleration of the accumulation of mutations during ribavirin monotherapy or during IFN-ribavirin combination therapy even when viral replication was profoundly inhibited by IFN-
. These results show that ribavirin does not act as an RNA mutagen in vivo, at least at the weight-based doses used for therapy and despite relatively high concentrations in plasma (the steady-state ribavirin concentration in peripheral blood was reached on day 14).
Although we examined only viruses circulating in peripheral blood, these viral populations closely reflect virions produced in the liver, since the average half-life of free HCV virions in the general circulation is of the order of 2.7 h and about 1012 virions are produced and cleared daily (32). In addition, there is no clear reason why viruses undergoing random nucleotide substitutions under the putative mutagenic effect of ribavirin should be preferentially sequestered in the liver while nonmutated virions would be normally excreted and released into peripheral blood.
Our findings therefore show that ribavirin has no noteworthy mutagenic effect in vivo, even when HCV replication is strongly inhibited by IFN administration. This is important in light of the recent suggestion by Dixit et al., based on mathematical modeling of viral decay in patients receiving the IFN-ribavirin combination, that the main action of ribavirin is to render HCV virions less infectious and that this property is revealed only when HCV replication is profoundly inhibited by the antiviral properties of coadministered IFN (10). If this is so (and it has not yet been experimentally confirmed), our results do not suggest that the underlying mechanism is error catastrophe due to accelerated mutagenesis.
The very high concentrations of ribavirin used in in vitro studies in which the mutagenic effect of ribavirin was observed could explain the discrepancy with our results obtained in the clinical setting. Indeed, the doses used in vitro were generally near toxic (50 to 1,000 µM, doses that cannot be achieved intracellularly during human therapy), and mutagenesis appeared to be concentration dependent. Interestingly, the ribavirin concentration required to achieve a significant antiviral effect differed between HCV-unrelated viruses (100 to 1,000 µM), GB virus B (100 to 200 µM), and subgenomic HCV replicons (50 to 400 µM) (1, 5, 8, 16, 20, 21, 48). Different cell types might have different capacities to convert ribavirin into RTP, and different rates of division (primary hepatocytes versus cultured cells) could influence the size of the intracellular GTP pool and the affinity of the polymerase for ribavirin. In any case, the ribavirin concentrations achieved with these models were substantially higher than those achieved intracellularly during human therapeutics (the latter cannot easily be measured, owing to the lack of an appropriate assay).
Our results differ from those of Asahina et al. (3), who reported a significantly higher rate of HCV mutations during 4 weeks of ribavirin monotherapy than prior to therapy (3). That study was based on direct sequence analysis; however, the approach is subject to misinterpretation of mixed or ambiguous profiles. In contrast, we used extensive quasispecies analysis to evaluate any possible change in the diversity of the population that could result from an increase of the virus polymerase error rate. In addition, Asahina's study suffered from a major methodological bias, since the pretreatment control period was variable (6 months on average; range, 2 to 48 months), whereas ribavirin monotherapy always lasted 4 weeks (3). We have found that error generation rates bear a negative relation to the interval between the two time points used to calculate them (Chevaliez and Pawlotsky, unpublished observation). This is due to the increase in the numerator (duration of the interval) in a context of slow spontaneous and random accumulation of mutations and to the use of short target sequences. As a result, the error generation rates during the pretreatment period in Asahina's study were probably substantially underestimated. In order to avoid this bias, we chose to use the same duration (14 days) for each of the three studied intervals, comprising the pretreatment period (control period), the first 14 days of treatment, and the second 14-day period, during which patients receiving combination therapy had variable levels of inhibition of HCV replication by IFN-
. With this approach we were able to show that mutation accumulation in two independent HCV genomic regions did not accelerate during treatment with ribavirin either alone or combined with IFN-
. These results are in keeping with preliminary studies which found no difference in error generation rates between patients receiving ribavirin monotherapy, IFN-
monotherapy, or IFN-ribavirin combinations (2, 43, 52).
In conclusion, our study, based on a robust methodology, strongly undermines the hypothesis whereby ribavirin acts as an HCV mutagen in vivo, leading to error catastrophe and viral eradication when combined with IFN. This is not altogether surprising, since the hypothesis was based on the use of in vitro models and near-toxic ribavirin concentrations. Unfortunately, the mechanism of ribavirin action remains as elusive as ever. A direct antiviral effect does not appear to play a major role, and neither does IMPDH inhibition (18, 28, 31, 37). A breakthrough in our understanding of how ribavirin works is urgently needed in order to develop drugs that mimic the antiviral properties of ribavirin without inducing severe hemolytic anemia, an adverse effect that often leads to dose reductions or withdrawal and subsequent treatment failure.
| ACKNOWLEDGMENTS |
|---|
We are grateful to Esteban Domingo for very helpful discussions and for his comments on the manuscript.
| FOOTNOTES |
|---|
Published ahead of print on 9 May 2007. ![]()
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| J. Bacteriol. | Mol. Cell. Biol. | Microbiol. Mol. Biol. Rev. |
|---|
| Clin. Vaccine Immunol. | ALL ASM JOURNALS |
|---|