Previous Article | Next Article ![]()
Journal of Virology, February 2003, p. 1940-1950, Vol. 77, No. 3
0022-538X/03/$08.00+0 DOI: 10.1128/JVI.77.3.1940-1950.2003
Copyright © 2003, American Society for Microbiology. All Rights Reserved.
Division of Infectious Diseases and Geographic Medicine, Stanford University School of Medicine, Stanford, California 94305,1 AIDS Research Center, Veterans Affairs Palo Alto Health Care System, Palo Alto, California 943042
Received 8 May 2002/ Accepted 30 October 2002
|
|
|---|
|
|
|---|
The clinical implications of coinfection have been the focus of intense research. HIV coinfection has been shown to increase the severity of liver disease in patients chronically infected with HCV (10, 34, 47). In addition, many studies have documented that HIV/HCV-coinfected patients have higher HCV loads than do HCV-monoinfected controls (4, 12, 44). In contrast to these well-documented effects of HIV on the course of HCV disease, the effect of HCV on the course of HIV disease is less clear. In some studies, coinfection with HCV has been shown to confer an increased risk for progression to AIDS in HIV-infected individuals (5, 14); however, other studies have failed to demonstrate this increased risk (53).
One of the hallmarks of HCV is its marked genetic diversity. On a population level, genotype 1 infections account for approximately two-thirds of all HCV infections in the United States (25) and may account for up to 83% of infections in the HIV/HCV-coinfected population (45). HCV also exists within an individual as a population of quasispecies (28). The region of the HCV genome with the greatest diversity is hypervariable region 1 (HVR1), located at the N terminus of the E2 envelope gene (19, 24). This region has been implicated to play a role in immune escape by virtue of its high degree of sequence variation (51). The significance of quasispecies diversity, however, is still unclear. Increased diversity within HVR1 has been associated with increased severity of liver disease (16, 21), although some studies have not found such an association (31). Increased quasispecies diversity has also been shown, during acute infection, to predict progression to persistent viremia and chronic HCV infection (9, 36). In addition, many studies have found an association between a higher pretreatment number of quasispecies in HVR1 and a poor response to interferon (IFN) therapy (13, 21).
The investigation into the evolution of HCV quasispecies in the HIV/HCV-coinfected population has yielded conflicting results. When comparing HIV/HCV-coinfected patients with HCV-monoinfected controls, two studies have found that coinfected patients have more quasispecies diversity (6, 43). If HIV infection were a surrogate marker for immunosuppression, this would imply that quasispecies diversity increases with immunosuppression. In support of this theory, Dove et al. have shown in coinfected patients that those with lower CD4+-cell counts have greater quasispecies evolution than those with higher CD4+-cell counts (L. M. Dove, Y. Phung, J. Wrock, M. Kim, and T. L. Wright, Abstract, Hepatology 30:456A, 1999). However, a number of other studies in HIV/HCV-coinfected individuals suggest that the opposite is true, i.e., that quasispecies diversity decreases with immunosuppression, manifested as lower quasispecies variation in those individuals with lower CD4+-cell counts (27, 38, 50).
There has been much focus recently on the administration of HAART in HIV/HCV-coinfected individuals. The effect of HAART on HCV load is still controversial. Most studies have shown no change in HCV RNA titers following HAART (11, 37, 49), although some studies have shown a transient (34, 39) or sustained (35) increase in HCV load; yet others have shown a decrease in HCV RNA levels and in some cases even HCV clearance (55).
To our knowledge, the effect of HAART on the quasispecies profile of HIV/HCV-coinfected individuals has not been described. Given the potential predictive value of quasispecies variation as a marker for HCV-related disease and IFN resistance, determining how HAART affects the genetic diversity of HCV will have important implications in both understanding the course of HCV disease in the setting of HIV infection, as well as in helping to better define the therapeutic management of coinfected individuals. The present data on quasispecies diversity and evolution in coinfected patients are conflicting. However, given that (i) HVR1 responds to immune pressure with increased variability (36, 51) and (ii) in other disease models it has been shown that immunosuppression (and thus decreased immune pressure) is associated with a decrease in HCV genetic diversity (2, 23, 26, 29, 33), we hypothesized that HAART, via immune restoration and increased immune pressure, might cause an increase in HCV quasispecies diversity. As such, we sought to determine the effect of HAART-associated immune restoration on the HCV quasispecies profile in HIV/HCV-coinfected individuals. In addition, we analyzed the effect of HAART on HCV load, as well as the differential effect of HCV genotype on quasispecies evolution during HAART.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Sociodemographic and HCV genotype characteristics of study groupsa
|
Oligonucleotide primers. The E1/E2 region targeted for amplification included the entire 81-bp HVR1 at the 5' end of the E2 gene as well as 170 bp of non-HVR1 sequence at the 3' end of the E1 gene. Because of the diversity of the HCV genome, and in order to avoid selective PCR amplification of different quasispecies, the primers for reverse transcription and amplification of the E1/E2 region were chosen based on the analysis by Toyoda et al. (50), in which primers were designed to be 80% homologous to 128 different reported HCV sequences. Primer sequences were as follows: HVR1 (TGGGACACATGATGATGAACTGGT) was used as the sense primer in both rounds of seminested PCR. HVR4 (CGGTGCTGTTTATGTGCCAACTGCC) was used as the antisense primer for both reverse transcription and the first round of PCR amplification. HVR3 (GATGTGCCAGCTGCCATTGG) was used as the antisense primer in the second round of PCR amplification.
RNA extraction and reverse transcriptase PCR (RT-PCR). Total RNA was extracted from 140 µl of patient sera by using the QIAamp Viral RNA Mini Kit (Qiagen, Valencia, Calif.) and resuspended in 60 µl of RNase-free water containing 0.04% sodium azide. Serial dilutions, when indicated, were made by 1:10, 1:100, and 1:1,000 dilution of this RNA sample. Two and a half microliters of the resulting RNA sample or dilution was then reverse transcribed using the GeneAmp RNA PCR Kit (Applied Biosystems, Foster City, Calif.) in a reaction volume of 20 µl containing 0.75 µM HVR4 primer, 50 mM KCl, 10 mM Tris-HCl (pH 8.3), 25 mM MgCl2, a 1 mM concentration of each deoxynucleoside triphosphate (dNTP), 1 U of RNase inhibitor, and 1 U of murine leukemia virus RT. The E1/E2 region was then amplified from the resulting cDNA by seminested PCR; the first round was performed in a 50-µl reaction volume containing 50 mM KCl, 10 mM Tris-HCl (pH 8.3), 2 mM MgCl2, 0.2 µM HVR1 sense primer, 0.3 µM HVR4 antisense primer, and 1.5 U of AmpliTaq DNA Polymerase (Applied Biosystems). PCR amplification was performed under the following conditions: 95°C for 105 s; 35 cycles of 95°C for 15 s and 60°C for 30 s; and 72°C for 7 min. The second round of PCR was performed in a 100-µl reaction volume containing 1x PCR buffer (Qiagen), 2 mM MgCl2, a 200 µM concentration of each dNTP, 0.3 µM HVR1 sense primer, 0.3 µM HVR3 antisense primer, and 2.5 U of Taq DNA polymerase (Qiagen). Thermal cycling was then performed as follows: 94°C for 105 s; 30 cycles of 94°C for 1 min, 56°C for 1 min, and 72°C for 1 min; and 72°C for 10 min. The resulting 318-bp PCR product was isolated by gel purification using the Qiaex II Gel Extraction Kit (Qiagen).
Cloning and sequencing. The amplified E1/E2 fragments were cloned into the pCR4-TOPO vector, and the resultant plasmids were used to transform chemically competent TOP10 cells using the TOPO TA Cloning Kit for Sequencing (Invitrogen, Carlsbad, Calif.). Colonies were screened for E1/E2-containing plasmids by PCR amplification of colonies in a 40-µl reaction volume containing 1x PCR buffer (Qiagen), 2 mM MgCl2, a 200 µM concentration of each dNTP, 0.07 µM M13 Forward primer (GTAAAACGACGGCCAG), 0.07 µM M13 Reverse primer (CAGGAAACAGCTATGAC), and 1 U of Taq DNA polymerase (Qiagen). Thermal cycling was performed as follows: 94°C for 10 min; 30 cycles of 94°C for 1 min, 55°C for 1 min, and 72°C for 1 min; and 72°C for 10 min. Ten positive clones per time point per patient were then sequenced using AP Biotech DYEnamic ET Terminator cycle sequencing (Amersham Biosciences, Piscataway, N.J.) and ABI PRISM (Applied Biosystems) sequencing technology.
Genetic analysis.
The 10 E1/E2 sequences per time point were aligned at the nucleotide level using EditSeq and MegAlign software (DNASTAR, Inc., Madison, Wis.). All sequence analysis was then performed with the 81-bp (27-amino-acid [aa]) HVR1 as well as with a 162-bp (54-aa) region of non-HVR1 sequence in the E1 gene. The program MEGA version 2.1 (22) was used to analyze the aligned sequences for the following parameters: (i) mean genetic distance, calculated as the mean of all pairwise comparisons of genetic distance using the p distance, defined as the number of amino acid differences divided by the total number of amino acid sites compared; and (ii) numbers of synonymous substitutions per synonymous site (Ks) and nonsynonymous substitutions per nonsynonymous site (Ka) using the Nei-Gojobori method with the Jukes-Cantor correction for multiple substitutions (32). The ratio of Ka/Ks was then calculated. The complexity of the quasispecies population per time point was analyzed at the amino acid level by two different methods: (i) determination of the total number of different clones within the 10-sequence population and (ii) calculation of the normalized Shannon entropy (Sn) to take into account the frequency of each different quasispecies in the population, calculated as follows: Sn = -
(pi ln pi)/ln N, where pi is the frequency of each sequence in the population and N is the total number of sequences analyzed (52). Entropy of the non-HVR1 E1 region was determined by taking the mean of the Sn values calculated for two separate 81-bp sequence blocks. Phylogenetic analysis was performed with the MEGA version 2.1 program (22) by using the neighbor-joining method and the Kimura two-parameter model.
Statistical analysis. Means are expressed as mean plus or minus standard deviation. Comparisons between groups were determined using the Fisher exact probability test, analysis of variance, or the Student t test, where appropriate. Comparison analyses and correlations were performed using StatView 5.0 software (SAS, Cary, N.C.). All reported P values are two-tailed, and a P of less than 0.05 was considered significant.
Nucleotide sequence accession number. Sequences reported herein have been assigned EMBL accession numbers AJ510768 through AJ511257.
|
|
|---|
Quasispecies parameters and comparison between HVR1 and the E1 region. HCV quasispecies were analyzed by assessing the diversity in HVR1, since this area has been shown to have a high degree of variation and has been implicated as playing a role in immune escape (19, 24, 51). Quasispecies parameters analyzed were (i) complexity, measured as both the total number of amino acid clones and Sn; (ii) mean genetic distance at the amino acid level; and (iii) Ka/Ks. Nonsynonymous changes (Ka) in HCV are likely to play a role in immune escape in response to immune pressure, especially in HVR1 (51). On the other hand, synonymous changes (Ks) represent genetic drift, which is more a function of the high replicative capacity and high mutation rate of HCV (36). Thus, the ratio of Ka/Ks is a good measure of immune pressure, since it normalizes for the genetic drift associated with replication and since values of Ka/Ks that are > 1 represent immune pressure (42).
In order to confirm that in our patient population the quasispecies evolution in HVR1 was indeed higher than that of the surrounding nonhypervariable E1 region, we compared the quasispecies parameters between these two regions of the HCV genome (Table 2). When the data for all patients were examined together, HVR1 had significantly greater quasispecies variation than the E1 region for all parameters analyzed: in HVR1, the number of clones was 1.6-fold higher (P < 0.0001), the entropy was 3.4-fold higher (P < 0.0001), the genetic distance was 11.7-fold higher (P < 0.0001), and the Ka/Ks ratio was 5.0-fold higher (P = 0.003). The results were similar when patients were stratified by treatment group (Table 2). Thus, given the significantly higher quasispecies complexity, genetic distance, and Ka/Ks in HVR1 than in an adjacent region of the E1 gene, HVR1 seemed to be an accurate marker of quasispecies variation in our HIV/HCV-coinfected population.
|
View this table: [in a new window] |
TABLE 2. Comparison of quasispecies parameters between HVR1 and the non-HVR1 (E1) regiona
|
![]() View larger version (15K): [in a new window] |
FIG. 1. Phylogenetic analysis of E1/E2 nucleotide sequences from the baseline time point of all patients. A phylogenetic tree was constructed from a total of 160 E1/E2 sequences using the neighbor-joining method and Kimura two-parameter model. The horizontalbranch lengths are drawn to scale, and the scale bar represents 0.1 nucleotide substitutions per site. The clone number is indicated at the end of each horizontal branch tip, and the brackets at the far right indicate groupings based on patient source. Patients in group A are indicated as A1 through A7, patients in group B as B1 through B3, and patients in group C as C1 through C6. Clone numbers are indicated by patient source followed by sequence number (e.g., A1-1, A1-2, etc.).
|
To ensure that analysis of only 10 sequences per time point was an accurate measure of the quasispecies population at that time point, 10 E1/E2 sequences were analyzed from each of two independent RT-PCRs using the same RNA template; five samples were randomly chosen for this analysis. To compare the quasispecies population obtained from each reaction, a phylogenetic tree was constructed for each of the five samples using both sets of 10 clones obtained from each RT-PCR for a total of 20 sequences per tree (Fig. 2). As can be seen in the two representative phylogenetic trees shown in Fig. 2, sequences do not cluster according to RT-PCR, indicating that both sets of sequences were drawn from a similar quasispecies pool. In order to quantitatively compare the HVR1 quasispecies parameters obtained for each RT-PCR, the mean values for the five samples in each group were compared (Table 3), a type of comparison similar to the intergroup analysis used elsewhere in this paper (see below). As such, there were no significant differences between the number of HCV clones, entropy, genetic distance, and Ka/Ks between the independent RT-PCR experiments for the group of five samples examined (P > 0.80). Taken together, the phylogenetic and quantitative quasispecies analyses above indicate that the sampling of 10 clones per time point was sufficient to accurately determine the quasispecies profile of that sample.
![]() View larger version (16K): [in a new window] |
FIG. 2. Phylogenetic analysis of E1/E2 nucleotide sequences from two independent RT-PCRs. Phylogenetic trees were constructed from 10 E1/E2 sequences for each RT-PCR (for a total of 20 sequences per sample) using the neighbor-joining method and Kimura two-parameter model. Representative trees are shown for two samples, A and B. The horizontal branch lengths are drawn to scale, and the scale bar indicates nucleotide substitutions per site. The clone number is indicated at the end of each horizontal branch tip. For each sample, clones from RT-PCR 1 are indicated as 1-1 through 1-10, and clones from RT-PCR 2 are indicated as 2-1 through 2-10.
|
|
View this table: [in a new window] |
TABLE 3. HCV quasispecies parameters determined by two independent RT-PCRsa
|
|
View this table: [in a new window] |
TABLE 4. Immunologic, virologic, and quasispecies characteristics of coinfected patients at baseline and ending time points
|
|
View this table: [in a new window] |
TABLE 5. Mean immunologic, virologic, and quasispecies characteristics of coinfected patients throughout the study perioda
|
Correlation between immunologic, virologic, and quasispecies parameters. In order to further investigate the association between increased CD4+- and CD8+-cell counts with higher HCV load and increased quasispecies diversity, we next assessed the relationship between these parameters by correlation analysis (Table 6). CD4+-cell count was positively correlated with HCV load (r = 0.545, P < 0.0001), genetic distance (r = 0.315, P = 0.02), and Ka (r = 0.316, P = 0.02). The CD4+-cell count was not significantly correlated with the number of HCV clones, entropy, Ks, or Ka/Ks. Similarly, CD8+-cell count was positively correlated with HCV load (r = 0.287, P = 0.04). There was also a trend toward positive correlation between CD8+-cell count and both genetic distance (r = 0.203, P = 0.16) and Ka (r = 0.248, P = 0.08), although neither reached statistical significance. There was no correlation between CD8+-cell count and number of HCV clones, entropy, Ks, or Ka/Ks.
|
View this table: [in a new window] |
TABLE 6. Correlation between immunologic, virologic, and quasispecies parametersa
|
Analysis of immunologic, virologic, and quasispecies parameters by genotype. We then sought to determine whether or not there were any differences in mean immunologic, virologic, or quasispecies parameters between patients infected with HCV genotype 1 and those infected with genotype 2 or 3 (Table 7). Intragroup comparisons were not made for group B because all three patients in this group were infected with genotype 1. Interestingly, patients infected with HCV genotype 2 or 3 had higher mean CD4+-cell counts than those infected with genotype 1. This was true for all patients taken together, as well as when patients were stratified into their respective treatment groups. The difference was most dramatic in group C, where patients infected with genotype 2 or 3 had a mean CD4+-cell count of 864, while those infected with genotype 1 had a mean CD4+-cell count of 264 (P < 0.0001). There was no significant difference in the CD8+-cell count of patients based on HCV genotype. In group C, but not group A, patients infected with genotype 2 or 3 had a higher HCV load (27.7 x 106 copies/ml versus 8.6 x 106 copies/ml; P = 0.005). Overall, patients with genotype 2 or 3 also had an increased genetic distance (0.146 versus 0.070 for all groups taken together; P = 0.01). When patients were analyzed by group, there was only a trend toward an increased genetic distance in patients infected with genotype 2 or 3: 0.12 versus 0.055 for group A (P = 0.14) and 0.174 versus 0.109 for group C (P = 0.20). In group A, but not group C, there was a higher Ka/Ks associated with genotype 2 or 3 infection (0.839 versus 0.399; P = 0.004). There was no significant difference in the number of HCV clones or entropy based on genotype. Taken together, patients infected with genotype 2 or 3 had significantly higher CD4+-cell counts and, in some cases, had a significantly higher HCV load, genetic distance, and Ka/Ks.
|
View this table: [in a new window] |
TABLE 7. Analysis of mean immunologic, virologic, and quasispecies parameters by genotypea
|
|
|
|---|
As a result of the above findings, we propose the following model to describe the relationship of HCV with immune status during HAART: initially, during the virologic containment of HIV, HAART does not have a major impact on quasispecies diversity (group A in our study). However, once virologic containment is stable and immune restoration has been maximized, there is increased immune pressure (increased Ka/Ks), which causes increased HCV quasispecies diversity as a means for immune escape (group C in our study). In addition, there may also be an increase in HCV load during this later stage, which may go hand in hand with increased quasispecies diversity (see below). In summary, this model describes the adaptive evolution of HCV in response to immune pressure: greater immune competence in the host results in increased quasispecies diversity by selecting for escape mutants from the preexisting pool of randomly generated quasispecies. In this manner, a less effective immune response does not exert such pressure for change on the quasispecies pool, resulting in a more homogenous population of viral quasispecies.
The increased immune pressure seen in patients on long-term HAART (group C) is likely a result of both quantitative (i.e., increased CD4+- and CD8+-cell counts) as well as qualitative (e.g., modification of cytokine expression) changes in immunologic status. Patients in group C had already achieved virologic containment of HIV for an average of 9 months prior to the onset of the study period. In comparison to patients who received a shorter duration of HAART (group A), the long-term virologic suppression in these patients may have allowed for better immune recovery (both quantitative and qualitative), consistent with the observation that HIV suppression for >1 year may be necessary to allow for complete and effective recovery of the immune system (30). In addition, it has been shown that HAART normalizes alterations in cytokine patterns induced by HIV infection (17, 30). These qualitative improvements in immune function may explain our observation that CD4+-cell counts, while correlated with genetic distance and Ka, were not correlated with the number of clones, entropy, or Ka/Ks. Absolute CD4+-cell counts may thus be only part of the immune driving force behind increased quasispecies diversity in patients on long-term HAART.
Interestingly, we found that CD4+- and CD8+-cell counts were positively correlated with HCV load. This is in contrast to most other studies, which have shown either an inverse relationship between HCV load and CD4+-cell count (4, 47) or no association at all (44, 45). The reason for these conflicting reports is unclear, although it may be due to differences in the immune function of study subjects or to different methods of studying the relationship between CD4+-cell count and HCV load. The biological significance of any change in HCV load with CD4+-cell count, when assessed for relevance to disease severity or response to IFN therapy, is unclear. In our study, we showed a threefold difference between the HCV load of group C (18.1 x 106 copies/ml) and that of group B (6.6 x 106 copies/ml). Daar et al. have shown that, for every 10-fold increase in HCV load, there is a increased relative risk of 1.66 for progression to AIDS and 1.54 for AIDS-related mortality (5). Therefore, a threefold increase in HCV load is unlikely to cause a substantial increased risk for a more rapid progression to AIDS. In addition, the impact of HCV load on HCV disease is unclear, as many studies have failed to find a correlation between HCV load and severity of liver disease (12, 41). On the other hand, it has been shown that an increase in pretreatment HCV load of only 1.2 x 106 copies/ml may be associated with a poor response to IFN therapy (48). Thus, the increased HCV load associated with long-term HAART in this study may not be predictive of an increased severity of HCV or HIV-related disease but may be important in predicting the response to subsequent IFN therapy.
In this study, we showed a significant positive correlation between HCV load and entropy, genetic distance, and Ka, and a trend toward a positive correlation between HCV load and both the number of clones and Ka/Ks. The data from previous studies on this relationship are not conclusive, with some studies reporting no correlation between quasispecies diversity and HCV load (9, 31, 36, 38) and others showing a positive correlation (13, 18). The reasons for this conflicting data are likely differences in study populations and methods used to evaluate quasispecies variation. The association between increased HCV load and greater quasispecies diversity demonstrated in this study can be explained by the adaptive evolution of HCV in the face of immune pressure: escape mutants generated in this way can replicate more freely, resulting in a higher HCV load. Similarly, a less robust immune response cannot select for escape mutants with such replicative fitness, and so patients with less HCV quasispecies diversity would also have a lower HCV load.
In this study, an association was found between lower CD4+-cell counts and decreased HCV quasispecies diversity in HIV/HCV-coinfected patients. These data are in accordance with results from other studies in coinfected patients showing that HCV quasispecies diversity decreases with the degree of HIV-related immunosuppression (27, 38, 50). Our data are also consistent with what has been shown in other disease models of immunosuppression: studies of HCV quasispecies in patients with agammaglobulinemia/hypogammaglobulinemia (2, 23) and in patients undergoing immunosuppressive therapy for liver transplantation (26, 29) or bone marrow transplantation (33) have all shown a decrease in HCV quasispecies diversity during states of immunosuppression. However, a few studies in HIV/HCV-coinfected patients have found the opposite resultthat HCV quasispecies diversity increases with HIV-induced immunosuppression. One study showed that coinfected patients with lower CD4+-cell counts had a greater percentage of new clones over a 1-year study period than did those with higher CD4+-cell counts (Dove et al., Hepatology 30:456A, 1999); however, this study involved only nine patients, and the authors did not report whether or not there was a change in genetic distance or Ka/Ks. Two other studies compared HIV/HCV-coinfected patients with HCV-monoinfected controls and found that coinfected patients had more quasispecies diversity (6, 43). One study compared only two patients (6), and neither study addressed the effect of CD4+-cell count differences within the coinfected population (6, 43). In addition, none of the aforementioned studies stratified patients based on exposure to HAART (6, 43; Dove et al., abstract).
Our results showing that HCV quasispecies diversity increases with immune pressure are also consistent with what has been described for HIV quasispecies evolution during the course of HIV monoinfection. Wolinsky et al. showed that in HIV-positive individuals over time, rapid CD4+-cell loss was associated with evolutionary stasis of HIV, while lower rates of CD4+-cell loss were associated with a greater accumulation of mutations and, in particular, nonsynonymous substitutions, indicating that selective pressure plays a role in the evolution of HIV quasispecies (52). Thus, HCV and HIV act similarly in response to immune pressure: greater immunologic competence is associated with an adaptive increase in the genetic diversity and evolution of both viruses.
The increased HCV quasispecies diversity observed in our group of patients on long-term HAART has one of two possible effects: (i) the quasispecies diversity is a marker of increased immune pressure by a more effective immune system, which leads to better control of HCV infection and thus less severe liver disease, or (ii) the increased quasispecies diversity results in the potential for more virulent or immunoresistant clones that cannot be adequately contained by the immune system (even though it might be functioning overall at a higher level), thus leading to more severe liver disease. In support of the first hypothesis, data from studies of HCV recurrence after liver transplantation have shown that increased quasispecies diversity and increased Ka/Ks have been associated with less severe HCV recurrence (26, 41). However, the majority of data in HCV-monoinfected patients supports the second hypothesisthat increased diversity within HVR1 will lead to increased severity of liver disease (16, 21). Further investigation into the significance of quasispecies diversity in these patients is warranted.
Many studies have found an association between higher pretreatment HCV quasispecies diversity in HVR1 and a poor response to IFN therapy (13, 21). This has important implications in the debate as to which infectionHIV or HCVto treat first in a coinfected patient. Yokozaki and colleagues (54, 55) have argued that HIV treatment should be undertaken first because HAART would increase CD4+-cell counts and may potentially decrease HCV load, thus providing a better starting point for IFN therapy since it has been shown that a lower HCV load (48) and higher CD4+-cell count (46) both predict a better response to IFN. Others have argued that HCV treatment should be undertaken first, since HAART may in fact increase HCV load, cause hepatotoxicity, and interact with HCV medications in a manner that may limit compliance (3). The data in our study support, in part, both arguments. On the one hand, the increase in CD4+-cell count associated with HAART suggests that treatment of HIV should precede treatment of HCV, as discussed above. However, on the other hand, we have shown that long-term HAART is associated with an increased HCV load and greater quasispecies diversity, both of which predict a poor response to IFNthis suggests that treatment for HCV should precede HAART in order to decrease the probability of a poor response to IFN.
A surprising finding in our study was that patients infected with genotype 2 or 3 had significantly higher CD4+-cell counts and, in some cases, had significantly higher HCV load, genetic distance, and Ka/Ks. Our results are in conflict with other studies in coinfected patients, which have shown an association between genotype 1 infection and increased quasispecies diversity (38) and HCV load (1, 49). However, none of these studies stratified patients based on exposure to HAART. Interestingly, genotype 1 has also been associated with increased severity of liver disease in HCV-monoinfected (56) and -coinfected (10) patients, as well as with a more rapid progression to AIDS (40), although patients in this study were on single-drug antiretroviral therapy and the data may not necessarily extrapolate to patients on HAART. The association of genotype 1 with increased severity of liver disease and more rapid progression to AIDS may be explained by a lower quasispecies diversity (see above discussion) and/or lower CD4+-cell count in these patients, although these parameters were not analyzed in these studies (10, 40, 56). To our knowledge, ours is the first report showing a difference in CD4+-cell count in association with HCV genotype. It may be that HCV, and in particular genotype 1, has an immunomodulatory effect in coinfected patients, possibly by direct interaction of HCV with HIV or via the alteration of cytokine patterns. For example, HCV coinfection has been associated with decreased levels of interleukin 18 (IL-18) and IL-1ß in coinfected patients (15); this in turn could lead to decreased CD4 cell proliferation and, in particular, Th1 cell proliferation and differentiation through the IL-18 pathway. In addition, some studies have found that HCV coinfection is associated with a blunted CD4+-cell response to HAART (14; C. Sabin, B. Dauer, A. N. Phillips, T. Lutz, V. Miller, A. C. Lepri, and S. Staszewski, Abstr. 9th Conf. Retrovir. Opportunistic Infect., abstr. 639-M, 2002), although these studies did not stratify patients by HCV genotype. It is provocative to envision a differential immunomodulatory effect of genotype 1, as is suggested by our results here. However, given the small number of patients in our study, we cannot speculate further on the implications of these results, as they demand confirmation in a larger study.
The potential limitations of our study include the relatively small number of patients in each group and the fact that our study population of veterans is mostly representative of the male HIV/HCV-coinfected population, whose main risk factor for infection is IVDU. In addition, 62% of our patients were infected with genotype 2 or 3; this is a higher percentage than that of the general coinfected population in the United States, in which 16% of HCV infections are with genotype 2 or 3 (45). We also focused our analysis solely on HVR1, and there may be other areas of the HCV genome that respond differently to immune pressure and/or to HAART, such as the IFN sensitivity-determining region. Another potential drawback of our study is that we analyzed only 10 clones per time point. We recognize that, the greater number of clones sequenced per time point, the better the assessment of the quasispecies population. However, given that the HCV quasispecies profile of a given group of samples was nearly identical for two independently obtained sets of 10 sequences (Fig. 2 and Table 3), we felt that the sampling of 10 clones per time point was adequate to achieve an accurate analysis of the HCV quasispecies population. In addition, while single-stranded conformational polymorphism is a better way to analyze the complexity of a larger numbers of clones, it does not permit analysis of genetic distance or Ka/Ks. Finally, the use of PCR amplification can at times lead to the introduction of mutations during cloning. However, in almost all cases, a mutation at a given position in HVR1 was present in more than one distinct clone (data not shown) and so was unlikely to represent a mutation introduced during PCR.
In summary, this study analyzes the effect of HAART and immunologic status on HCV load and quasispecies diversity in HIV/HCV-coinfected patients. To our knowledge, this is the first study analyzing the effect of HAART on HCV quasispecies variation. We have shown that there is no immediate effect of HAART on HCV quasispecies or load but that, after long-term HAART, patients had a higher HCV load and increased quasispecies diversity. We hypothesize that the higher immune pressure associated with maximal immune recovery in a given patient after HAART drives HCV to evolve more extensively in an attempt to create escape mutants. The potential implications for increased HCV quasispecies diversity in the coinfected population have been discussed and demand further investigation, especially given the importance of HCV disease in the long-term management of HIV/HCV-coinfected patients.
We thank Bayer Diagnostics for providing HCV load assay results. We also thank Sharon Lindsay for performing HCV genotype analysis and Suparna Dutt for her assistance in the preparation of HCV RNA.
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»