|
|
||||||||
1 Department of Clinical Molecular Informative Medicine, Nagoya City University Graduate School of Medical Sciences, Kawasumi, Mizuho, Nagoya 467-8601, Japan
2 National Institute of Genetics, Yata 1111, Mishima, Shizuoka, Japan
3 Ogikubo hospital, Tokyo, Japan
Correspondence
Yasuhito Tanaka
ytanaka{at}med.nagoya-cu.ac.jp
| ABSTRACT |
|---|
|
|
|---|
The GenBank/EMBL/DDBJ accession numbers for the sequences reported in this study are AB245555–AB245873.
| INTRODUCTION |
|---|
|
|
|---|
The estimated HCV virion half-life time was, on average, 2.7 h with pre-treatment production and clearance of 1012 virions per day (Neumann et al., 1998
). Such a high rate of HCV replication, combined with lack of an error correction mechanism, results in the development of genetically diverse clones in a patient. The genetic diversity of HCV has provided critical insights into short-term outcomes, including early spontaneous viral clearance (Farci et al., 2000
), interferon-associated viral clearance (Farci et al., 2002
; Pawlotsky et al., 1999
) and HCV emergence following liver transplantation (Lyra et al., 2002
). To infer the genetic diversity of HCV in a patient, we applied two approaches. One method simply assumed that the genetic diversity of HCV is of different divergence of synonymous distance between two time points. The other method applied to coalescent analysis of genetic diversity along the time, assuming that the genetic diversity of HCV represents a heterogeneous viral population in a given carrier. To evaluate the influence on the HCV evolution exerted by the immunosuppression during persistent HIV coinfection, we examined a cohort of HCV carriers, comparing the HIV-positive and -negative groups.
Determination of the antigen-recognition regions associated with HCV-specific immune positive selection is important for understanding selective pressures underlying the evolution of HCV as well as putative therapeutic targets. In this study, we evaluated the genetic diversity of HCV and determined genomic regions associated with positive selection by comparative analysis of selective forces between HIV-positive and -negative groups in a cohort of haemophilia patients followed for more than 6 years.
| METHODS |
|---|
|
|
|---|
Laboratory tests.
Laboratory evaluation included complete blood cell count and serum transaminases [alanine aminotransferase (ALT)]. CD4+ cell counts were examined by fluorescence-activated cell sorting at SRL Inc. (Tokyo, Japan). Serum HCV-RNA levels and HIV-RNA levels were measured by a commercial PCR assay (Amplicor HIV-1 Monitor and Amplicore HCV monitor; Roche Diagnostics). The detection limits of PCR for HCV-RNA and HIV-RNA were 500 IU ml–1 [0.5 kilo international unit (KIU) ml–1] and 50 copies ml–1, respectively.
HCV-RNA isolation and amplification from the core, E1 and E2 regions.
Nucleic acids were extracted from serum samples using a SepaGene RV-R Nucleic Acid Extracting kit (Sanko Junyaku) in accordance with the manufacturer's protocol. Viral RNA was reverse-transcribed to cDNA using SuperScript II RNase H– Reverse Transcriptase (Invitrogen) and random hexamer primer (Takara Shuzo) as described previously (Ohno et al., 1997
).
Partial core, E1 and E2 fragments were amplified by using PCR with primers as described previously (Tanaka et al., 2002
). To reduce the number of artificial substitutions arising during PCR, Platinum Pfx DNA Polymerase (Invitrogen) with a very high fidelity was used.
Cloning and sequencing of cDNA.
The amplified products were ligated into pCR-Blunt II-TOPO Vector and used to transform DH5-
high-efficiency competent cells according to the manufacturer's protocol (Invitrogen). The plasmid DNA was purified using the QIAprep Spin Miniprep kit (Qiagen) and the presence of the inserts confirmed by digestion with EcoRI. Sequencing was performed on more than 10 clones per patient at the baseline (1995–1997) and the end point (2002–2003). All clones were sequenced with Prism Big Dye (Applied Biosystems) in an ABI 3100 DNA automated sequencer.
Construction of phylogenetic trees.
Nucleotide sequences of HCV were aligned by using the program CLUSTAL_X and molecular evolutionary analyses were conducted using Molecular Evolutionary Genetic Analysis software (MEGA version 3.0; Kumar et al., 2001
). The MEGA algorithms were used to calculate the mean Tamura–Nei pairwise distance for all clones as well as a matrix of Tamura–Nei pairwise distances for each patient. To confirm the reliability of the phylogenetic tree, bootstrap resampling tests were performed 1000 times.
Genetic diversity of HCV over a time course.
Two approaches were used to infer the genetic diversity of HCV in each patient. In the first approach, total evolutionary distances among a heterogeneous viral population were compared between the baseline and end point for each patient in the phylogenetic tree. The phylogenetic tree of genetic diversity was constructed by using the maximum-likelihood (ML) method and the ancestral sequence was inferred at every node using the ML method (Yang et al., 1995
). As the evolutionary distance in each branch, the number of synonymous substitutions per synonymous site (synonymous distance) was estimated by the modified Nei–Gojobori method. Total synonymous distances were assumed to represent the genetic diversity of a heterogeneous viral population in each patient.
The other approach is the coalescence theory based on estimation of the genetic diversity. A consensus sequence based on the sequences of all HCV clones isolated from each patient was used as an outgroup to locate the position of the root in each phylogenetic tree. The topology of the phylogenetic tree was estimated by the neighbour-joining method (PHYLIP). Based on the topology, we constructed the phylogenetic tree and inferred the evolutionary rate by the ML method under the premise of the molecular clock (TipDate) (Rambaut, 2000
). Based on the trees and the evolutionary rates estimated by TipDate, the coalescent analysis of genetic diversity was conducted for each patient using the Genie v3.5 software (Pybus et al., 2001
; Pybus & Rambaut, 2002
). In brief, time t was transformed to year using the HCV molecular evolutionary rate, assuming the sample-collection time to be the present. Function N(t) (effective numbers of HCV infections through time) was estimated by the ML method to infer the genetic diversity of HCV (Pybus & Rambaut, 2002
). Although there are several models to infer N(t), the best-fit model was different among patients. Therefore, we chose a simplified model in which the genetic diversity was assumed to be exponentially increased over time (expansion model).
Identification of positively selective regions.
Positively selected regions were identified using the modified method of Suzuki & Gojobori (2001)
. In brief, a phylogenetic tree of sequences from HCV clones was reconstructed in each patient by the ML method. The ancestral sequence was inferred at every node in the phylogenetic tree using the ML method (Yang et al., 1995
). Then, synonymous and non-synonymous substitutions throughout the phylogenetic tree were estimated in each branch for each codon site. Here, to see the differences in selective pressure for HCV between the HIV-positive and -negative groups, we independently summed the total numbers of synonymous (Ns) and non-synonymous (Nn) substitutions occurring at each codon site of the HCV clones from either eight patients infected with HIV or five patients without HIV infection. The mean numbers of synonymous (Cs) and non-synonymous (Cn) sites were calculated for each codon site by the modified Nei–Gojobori method. The genetic distance of synonymous (dS) and non-synonymous (dN) was calculated as Ns/Cs and Nn/Cn, respectively. Although the ratio dN/dS is usually used for estimating selective pressure, we used (dN+0.5)/(dS+0.5) ratio instead in the present study, because no synonymous substitution was found in several codon sites. The ratio was calculated along with the sequence by the sliding-window analysis. Each window size consisted of three codons.
| RESULTS |
|---|
|
|
|---|
|
Long-term intra-host diversity of HCV evaluated on distinct genomic regions
It has been shown previously that the genetic diversity of HCV changes in an oscillatory manner during the natural course of the infection (Devereux et al., 1997
). Taking into account that the genetic diversity of HCV analysed at a single time point might not accurately reflect the dynamic profile of the population over time, we have examined 26 serum samples collected from 13 patients at two distinct time points with intervals of at least 6 years (6–8 years). At least 11 HCV clones were isolated from a single patient at the baseline (1995–1997) and at the end point (2002–2003) of the follow-up. Overall, 325 HCV clones were thus isolated and analysed. All of them belonged to genotype 1a. Phylogenetic relation of the HCV clones isolated from all patients is shown in Fig. 1
. Assuming that HCV is composed of a heterogeneous viral population, which is evolving throughout time in a given host (carrier), we aimed to estimate the size and heterogeneity of the population. Two different methods were used to attain this aim.
|
|
|
|
|
| DISCUSSION |
|---|
|
|
|---|
In the present study, the genetic diversity of HCV was higher in the HIV-positive group compared with the HIV-negative group, which could be associated with either or both higher rate of HCV replication and HIV-associated immunosuppression, leading to less selective pressure on HCV in the HIV-positive group. The increased HCV replication by immune dysregulation in the HIV/HCV-coinfected patients might damage liver cells through apoptosis and other means such as cytokine disruption as reported previously (Puoti et al., 2001
). Although several studies showed that HCV diversity decreases with the degree of HIV-related immunosuppression (Babik & Holodniy, 2003
; Mao et al., 2001
; Martell et al., 1992
; Qin et al., 2005
; Roque-Afonso et al., 2002
; Toyoda et al., 1997
), this might be associated with the different parameters assessed for the genetic diversity; the conclusions in most previous studies were made by analysing total numbers of HCV clones and overall genetic distances at the amino acid level.
To assess the discrepancy with our results by coalescence-based estimation, we further examined the role of HCV-targeted immune pressure. The HCV nucleotide substitution pattern was compared between subjects with intact immune system versus those with HIV-associated immunosupression using a recently designed approach. Previous reports had used the mean of pairwise synonymous and non-synonymous distances within isolates (Ray et al., 2000
; Blackard et al., 2004
). Since most isolates do not have an independent evolutionary process, the mean may not represent overall genetic diversity of the heterogeneous viral population. To examine the non-redundant evolutionary process, recent methods have inferred the evolutionary process throughout the phylogenetic tree (Sheridan et al., 2004
; Suzuki & Gojobori, 2001
; Hanada et al., 2006
). In the present study, we have applied a new approach to examine selection of HCV affected by HIV. The approach supported the theory that the diversified population of HIV-positive patients is due to less selective pressure and allows identification of specific regions indicating the presence of positive selection in HIV-negative patients compared with HIV-positive patients with immunosuppression. Although these positively selective segments, which were observed only in the HIV-negative group, might have relatively weak antigenicity, most of them were located inside potential immunogenic targets and others might be somewhat new antigen-recognition regions associated with HCV-specific immune responses. Interestingly, no influence by HIV coinfection was observed in the HVR1, which contains sequence-specific immunological B-cell epitopes that induce the production of antibodies restricted to the specific viral isolate (Kato et al., 1993
), indicating that the positively selected segments regardless of HIV infection should have strong antigenic epitopes. Taken together, our findings indicate that defenceless HCV clones that are extinct in usual conditions can survive in HIV-positive patients because of less immune pressure leading to HIV infection. Consequently, the genetic diversity of HCV will be greater in HIV-positive patients. In fact, the diversity of both synonymous and non-synonymous substitutions was larger in HIV-positive patients than HIV-negative patients (data not shown).
One limitation of the present study is the general lack of functional immunological data. In this study, HCV-specific CD4 ELISPOT responses were not detected in all subjects and CD8 cell counts were not measured. Therefore, we chose to use HIV load and CD4 cell counts as surrogate markers of immune suppression. Further investigations of HCV diversity in conjunction with HCV-specific cellular responses will be required when more-sensitive immunological assays are available. Another potential limitation of the study is that only two time points were sampled for each individual. However, as we inferred the evolutionary process based on the phylogenetic tree constructed using the number of clones that were isolated through the long period of follow-up (more than 6 years), we believe that the genetic diversity of the heterogeneous viral population may represent an actual evolutionary process.
In conclusion, HIV-positive patients have more diversified HCV populations than HIV-negative patients, possibly because of reduction of selective pressure from the immune system. The positively selective regions determined in this study might be antigen-recognition regions associated with HCV-specific immune responses.
| ACKNOWLEDGEMENTS |
|---|
| REFERENCES |
|---|
|
|
|---|
Benhamou, Y., Bochet, M., Di Martino, V., Charlotte, F., Azria, F., Coutellier, A., Vidaud, M., Bricaire, F., Opolon, P. & other authors (1999). Liver fibrosis progression in human immunodeficiency virus and hepatitis C virus coinfected patients. The Multivirc Group. Hepatology 30, 1054–1058.[CrossRef][Medline]
Bica, I., McGovern, B., Dhar, R., Stone, D., McGowan, K., Scheib, R. & Snydman, D. R. (2001). Increasing mortality due to end-stage liver disease in patients with human immunodeficiency virus infection. Clin Infect Dis 32, 492–497.[CrossRef][Medline]
Blackard, J. T., Yang, Y., Bordoni, P., Sherman, K. E. & Chung, R. T., AIDS Clinical Trials Group 383 Study Team (2004). Hepatitis C virus (HCV) diversity in HIV-HCV-coinfected subjects initiating highly active antiretroviral therapy. J Infect Dis 189, 1472–1481.[CrossRef][Medline]
Braitstein, P., Palepu, A., Dieterich, D., Benhamou, Y. & Montaner, J. S. (2004). Special considerations in the initiation and management of antiretroviral therapy in individuals coinfected with HIV and hepatitis C. AIDS 18, 2221–2234.[CrossRef][Medline]
Devereux, H. L., Brown, D., Dusheiko, G. M., Emery, V. C. & Lee, C. A. (1997). Long-term evolution of the 5'UTR and a region of NS4 containing a CTL epitope of hepatitis C virus in two haemophilic patients. J Gen Virol 78, 583–590.[Abstract]
Eyster, M. E., Fried, M. W., Di Bisceglie, A. M. & Goedert, J. J. (1994). Increasing hepatitis C virus RNA levels in hemophiliacs: relationship to human immunodeficiency virus infection and liver disease. Multicenter Hemophilia Cohort Study. Blood 84, 1020–1023.
Farci, P., Shimoda, A., Coiana, A., Diaz, G., Peddis, G., Melpolder, J. C., Strazzera, A., Chien, D. Y., Munoz, S. J. & other authors (2000). The outcome of acute hepatitis C predicted by the evolution of the viral quasispecies. Science 288, 339–344.
Farci, P., Strazzera, R., Alter, H. J., Farci, S., Degioannis, D., Coiana, A., Peddis, G., Usai, F., Serra, G. & other authors (2002). Early changes in hepatitis C viral quasispecies during interferon therapy predict the therapeutic outcome. Proc Natl Acad Sci U S A 99, 3081–3086.
Goedert, J. J., Eyster, M. E., Lederman, M. M., Mandalaki, T., De Moerloose, P., White, G. C., II, Angiolillo, A. L., Luban, N. L., Sherman, K. E. & other authors (2002). End-stage liver disease in persons with hemophilia and transfusion-associated infections. Blood 100, 1584–1589.
Graham, C. S., Baden, L. R., Yu, E., Mrus, J. M., Carnie, J., Heeren, T. & Koziel, M. J. (2001). Influence of human immunodeficiency virus infection on the course of hepatitis C virus infection: a meta-analysis. Clin Infect Dis 33, 562–569.[CrossRef][Medline]
Hanada, K., Gojobori, T. & Li, W. H. (2006). Radical amino acid change versus positive selection in the evolution of viral envelope proteins. Gene 385, 83–88.[CrossRef][Medline]
Hanada, K., Tanaka, Y., Mizokami, M., Gojobori, T. & Alter, H. J. (2006). A reduction in selective immune pressure during the course of chronic hepatitis C correlates with diminished biochemical evidence of hepatic inflammation. Virology 361, 27–33.[CrossRef][Medline]
Kato, N., Sekiya, H., Ootsuyama, Y., Nakazawa, T., Hijikata, M., Ohkoshi, S. & Shimotohno, K. (1993). Humoral immune response to hypervariable region 1 of the putative envelope glycoprotein (gp70) of hepatitis C virus. J Virol 67, 3923–3930.
Kumar, S., Tamura, K., Jakobsen, I. B. & Nei, M. (2001). MEGA2: Molecular Evolutionary Genetics Analysis software. Bioinformatics 17, 1244–1245.
Lai, C. L., Ratziu, V., Yuen, M. F. & Poynard, T. (2003). Viral hepatitis B. Lancet 362, 2089–2094.[CrossRef][Medline]
Lyra, A. C., Fan, X., Lang, D. M., Yusim, K., Ramrakhiani, S., Brunt, E. M., Korber, B., Perelson, A. S. & Di Bisceglie, A. M. (2002). Evolution of hepatitis C viral quasispecies after liver transplantation. Gastroenterology 123, 1485–1493.[CrossRef][Medline]
Mao, Q., Ray, S. C., Laeyendecker, O., Ticehurst, J. R., Strathdee, S. A., Vlahov, D. & Thomas, D. L. (2001). Human immunodeficiency virus seroconversion and evolution of the hepatitis C virus quasispecies. J Virol 75, 3259–3267.
Martell, M., Esteban, J. I., Quer, J., Genesca, J., Weiner, A., Esteban, R., Guardia, J. & Gomez, J. (1992). Hepatitis C virus (HCV) circulates as a population of different but closely related genomes: quasispecies nature of HCV genome distribution. J Virol 66, 3225–3229.
Neumann, A. U., Lam, N. P., Dahari, H., Gretch, D. R., Wiley, T. E., Layden, T. J. & Perelson, A. S. (1998). Hepatitis C viral dynamics in vivo and the antiviral efficacy of interferon-alpha therapy. Science 282, 103–107.
Ohno, T., Mizokami, M., Wu, R. R., Saleh, M. G., Ohba, K., Orito, E., Mukaide, M., Williams, R. & Lau, J. Y. (1997). New hepatitis C virus (HCV) genotyping system that allows for identification of HCV genotypes 1a, 1b, 2a, 2b, 3a, 3b, 4, 5a, and 6a. J Clin Microbiol 35, 201–207.[Abstract]
Pawlotsky, J. M., Germanidis, G., Frainais, P. O., Bouvier, M., Soulier, A., Pellerin, M. & Dhumeaux, D. (1999). Evolution of the hepatitis C virus second envelope protein hypervariable region in chronically infected patients receiving alpha interferon therapy. J Virol 73, 6490–6499.
Poynard, T., Yuen, M. F., Ratziu, V. & Lai, C. L. (2003). Viral hepatitis C. Lancet 362, 2095–2100.[CrossRef][Medline]
Puoti, M., Bonacini, M., Spinetti, A., Putzolu, V., Govindarajan, S., Zaltron, S., Favret, M., Callea, F., Gargiulo, F. & other authors (2001). Liver fibrosis progression is related to CD4 cell depletion in patients coinfected with hepatitis C virus and human immunodeficiency virus. J Infect Dis 183, 134–137.[CrossRef][Medline]
Pybus, O. G. & Rambaut, A. (2002). GENIE: estimating demographic history from molecular phylogenies. Bioinformatics 18, 1404–1405.
Pybus, O. G., Charleston, M. A., Gupta, S., Rambaut, A., Holmes, E. C. & Harvey, P. H. (2001). The epidemic behavior of the hepatitis C virus. Science 292, 2323–2325.
Qin, H., Shire, N. J., Keenan, E. D., Rouster, S. D., Eyster, M. E., Goedert, J. J., Koziel, M. J. & Sherman, K. E. (2005). HCV quasispecies evolution: association with progression to end-stage liver disease in hemophiliacs infected with HCV or HCV/HIV. Blood 105, 533–541.
Rambaut, A. (2000). Estimating the rate of molecular evolution: incorporating non-contemporaneous sequences into maximum likelihood phylogenies. Bioinformatics 16, 395–399.
Ray, S. C., Mao, Q., Lanford, R. E., Bassett, S., Laeyendecker, O., Wang, Y. M. & Thomas, D. L. (2000). Hypervariable region 1 sequence stability during hepatitis C virus replication in chimpanzees. J Virol 74, 3058–3066.
Roque-Afonso, A. M., Robain, M., Simoneau, D., Rodriguez-Mathieu, P., Gigou, M., Meyer, L. & Dussaix, E. (2002). Influence of CD4 cell counts on the genetic heterogeneity of hepatitis C virus in patients coinfected with human immunodeficiency virus. J Infect Dis 185, 728–733.[CrossRef][Medline]
Sheridan, I., Pybus, O. G., Holmes, E. C. & Klenerman, P. (2004). High-resolution phylogenetic analysis of hepatitis C virus adaptation and its relationship to disease progression. J Virol 78, 3447–3454.
Suzuki, Y. & Gojobori, T. (2001). Positively selected amino acid sites in the entire coding region of hepatitis C virus subtype 1b. Gene 276, 83–87.[CrossRef][Medline]
Tanaka, Y., Hanada, K., Mizokami, M., Yeo, A. E., Shih, J. W., Gojobori, T. & Alter, H. J. (2002). A comparison of the molecular clock of hepatitis C virus in the United States and Japan predicts that hepatocellular carcinoma incidence in the United States will increase over the next two decades. Proc Natl Acad Sci U S A 99, 15584–15589.
Toyoda, H., Fukuda, Y., Koyama, Y., Takamatsu, J., Saito, H. & Hayakawa, T. (1997). Effect of immunosuppression on composition of quasispecies population of hepatitis C virus in patients with chronic hepatitis C coinfected with human immunodeficiency virus. J Hepatol 26, 975–982.[CrossRef][Medline]
Yang, Z., Kumar, S. & Nei, M. (1995). A new method of inference of ancestral nucleotide and amino acid sequences. Genetics 141, 1641–1650.[Abstract]
Received 6 March 2007;
accepted 3 May 2007.
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| INT J SYST EVOL MICROBIOL | MICROBIOLOGY | J GEN VIROL |
| J MED MICROBIOL | ALL SGM JOURNALS | |