Previous Article | Next Article ![]()
Journal of Virology, April 2006, p. 4196-4199, Vol. 80, No. 8
0022-538X/06/$08.00+0 doi:10.1128/JVI.80.8.4196-4199.2006
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
Division of Infectious Diseases and Geographic Medicine, Stanford University, Stanford, California,1 Department of Epidemiology, Boston University School of Public Health, Boston, Massachusetts,2 AIDS Research Center, VA Palo Alto Health Care System, Palo Alto, California3
Received 7 October 2005/ Accepted 1 February 2006
|
|
|---|
|
|
|---|
Several reports have described NS3 diversity using bulk sequencing; no significant difference in diversity between groups with HCV monoinfection or with human immunodeficiency virus (HIV) coinfection has been demonstrated (1, 3, 9, 12). However, the levels of complexity and diversity of HCV protease gene variants and the resultant impact on HIV coinfection are unknown. Understanding this diversity and the impact of naturally occurring polymorphisms, which could result in reduced susceptibility, is crucial for studies of PI development. We performed a clonal analysis of HCV protease sequence diversity across subtypes in 24 chronically infected patients with HIV coinfection.
The Stanford University Institutional Review Board approved this study. HCV genotyping was performed using an INNO-LiPA HCV II assay (Bayer Diagnostics, Tarrytown, NY). Quantification of HIV and of HCV RNA was performed using an Amplicor HIV-1 Monitor v1.5 and HCV v2.0, respectively (Roche Diagnostics, Branchburg, NY).
RNA was extracted from 200 µl of patient sera using TRIzol (Invitrogen, Carlsbad, CA), resuspended in water, and reverse transcribed using Superscript II (Invitrogen) according to the manufacturer's instructions. The resulting cDNA was amplified by PCR using genotype-specific primers (Operon Technologies, Alameda, CA) and PlatinumTaq DNA polymerase (Invitrogen), generating fragments of approximately 1,000 bp. Primers for PCR amplification were as follows: for genotype 1a, G and J (9) or 1a-3258 (CCAGTCGTCTTCTCCCAAATG) and 1a-4150R (TTGGACATGTAAGCACCAAAGC); for genotype 1b, 1b-3145 (AAATGGCTCTCATGAAGCTGG) and 1b-4031R (AGCGTGTAGATGGGCCACTT); for genotype 2b, 2b-3107f (GCTACGAGTGTGTACCCTGGTG) and 2b-4109r (GAGTACTTTATACCCCTGACTGGCATA); and for genotype 3a, 3a-3290f (CCATGGAAATCAAGGTCATCAC) and 3a-4298r (ATAATCACATCATATGCCCCCC).
Amplified NS3 fragments were cloned using a TOPO TA cloning kit (Invitrogen) according to the manufacturer's instructions. This amplification and cloning approach was used previously to study HCV envelope gene diversity (3). Clones were sequenced with M13 primers using ABI PRISM (Applied Biosystems, Foster City, CA) sequencing technology (Macrogen, Seoul, Korea). Nucleotide sequences were assembled and aligned using AutoAssembler (Applied Biosystems). Individual patient consensus sequences were prepared by alignment of individual clones using Sequence Navigator (Applied Biosystems). Genetic distances between consensus or individual clone nucleic acid or amino acid sequences were determined using codon-aligned multiple sequence alignments and distance matrices. Sequences were submitted to GenBank.
Statistical analyses included Kruskal-Wallis rank sum tests to evaluate differences in the distribution of demographics, virologic and immunologic indicators of the stage of HIV infection, and nucleotide and amino acid diversity (genetic distances, calculated using GCG software [Madison, WI]) expressed as median values and the associated 25% to 75% interquartile ranges (IQR) for each group. Fixed- and random-effect mixed models were used to evaluate the contributions of HCV genotype and within-patient clustering of HCV clone nucleotide or amino acid sequences to observed nucleotide or amino acid differences by pairwise comparisons among clones from a given patient. Mixed models (SAS, version 8; Cary, NC) were constructed to include effects of genotype (fixed) and within-patient effects (random) as dependent variables, with nucleotide or amino acid differences as the independent variable. Mixed models (Y = ßX +
Z +
) were a variation of a general linear model, where the dependent variable Y (the percent difference in nucleotide or amino acid sequence among clones from a given patient) is a function of HCV genotype (ß, the fixed effect) or clustering within patients (
, the random effect);
is the error term. Therefore, there are two levels of clustering, with HCV clones being clustered within patients and patients being clustered within HCV genotypes.
Since up to 10 virus clones were sequenced per patient, mixed models permitted adjustments for potential within-patient clustering (correlation) of percent nucleotide and amino acid differences derived from individual HCV clone sequence data. Because dependent variables (nucleotide and amino acid differences) were initially non-normally distributed, they were log transformed prior to inclusion in mixed models. Least square means were used to ascertain nucleotide and amino acid differences by specific genotype by using pairwise comparisons.
Patient demographics, laboratory data, and treatment history are presented in Table 1. The distribution of HCV genotypes presented is representative of the population seen at our institution (5). HCV viral loads were significantly greater in the patients with genotype 1a and 2b than in those with genotypes 1b and 3a (Table 1). Mean CD4 count and percentage were significantly lower for genotype 1a than for other genotypes. Patient consensus sequences created from each individual's clone sequences demonstrated that between-genotype variability ranged from 24% (1a versus 1b) to 51% (2b versus 3a) at the nucleotide level and from 10% (1a versus 1b) to 40% (2b versus 3a) at the amino acid level. These results are similar to those found in another study of HCV NS3 diversity that used bulk sequencing (3, 9), and the inter- and intragenotype diversity is similar to interclade and interpatient differences in HIV strains (15, 17).
|
View this table: [in a new window] |
TABLE 1. Patient characteristics by HCV genotypea
|
Individual clone sequences for each individual (Table 2) showed an overall median percent nucleotide difference of 0.74% (IQR, 0.65 to 1.49%). However, individual clones had nucleotide differences as high as 5.34% compared to other clones for the same patient. When data from all clones across patients were combined, genotype 1a had a higher nucleotide variability between clones than the other genotypes evaluated (P = 0.01). At the amino acid level, there was a median percent amino acid difference of 0.55% (IQR, 0 to 1.11%) between clones within individuals. However, some clones were as much as 5.75% different in their amino acid sequences compared to other clones from the same patient (Table 2). Genotype 1b clones were significantly less variable than those of the other genotypes (P = 0.007).
|
View this table: [in a new window] |
TABLE 2. Nucleotide and amino acid differences between NS3 protease clones
|
The total numbers of positions with at least one amino acid change (compared to other patients' consensus sequences in the same genotype) were as follows: genotype 1a, 14 residues; 1b, 11 residues; 2a, 11 residues; and 3a, 9 residues. Two clones had amino acid changes at active site residues, a genotype 1a with S139L and 3a with H57N. There were also two clones (from the same genotype 1a patient), C97R and H149R, that had mutations at metal binding positions. There were no specific amino acid mutations that were found to confer reduced susceptibility to HCV PIs (A156T, R155Q, D168V/Y/A) (11, 13, 19). However, one 3a clone had a R155T mutation, one 1a clone had a D168E mutation, and all 3a clones had a D168Q substitution.
Other reports have described NS3 protease sequence variability using bulk sequencing (3, 9, 12). Our clonal analysis provides a more detailed evaluation of sequence variability by detecting minority variants and establishing linkage as demonstrated in studies of HIV diversity (7, 14). While there are significant differences in nucleic and amino acid sequences and enzymatic parameters (Ki, Km, etc.) (4) between genotypes, the lack of significant nucleotide and amino acid sequence variability within each genotype in our findings does not provide a basis for a genotype-dependent response to current PIs in development, as has been suggested (18).
While it is possible that bias due to differences in viral load could influence measures of diversity, the likelihood of sample bias or founder effects in this study is low, as the minimum virus load of all patients was 205,000 HCV RNA copies/ml, and only a single round of PCR was employed. Indeed, no significant correlation between HCV viral load and genetic diversity was observed. Overall, changes in important protease positions were rare, and some may have been introduced during Taq-mediated amplification. Thus, further analysis of these variants may be warranted to evaluate their prevalence in other patient groups and their effect on enzyme function.
Resistance to HCV PIs will be an important factor in managing the long-term use of these agents, similar to that seen with HIV protease gene resistance after HIV PI use (8). None of the mutations definitively associated with resistance to HCV PIs was seen in our untreated patients. However, amino acid substitutions at some of those positions were found, suggesting that genotype 1a may develop resistant strains early during PI therapy. Observations with HIV PI resistance suggest that multiple amino acid variants at known drug resistance positions will affect susceptibility (16).
The additional impact of HIV infection or treatment on NS3 diversity is largely unknown, as we did not analyze patients with HCV monoinfection. Previously, we have shown that HIV treatment significantly increases HCV E2 gene diversity and complexity in patients coinfected with HIV and HCV (2). Longitudinal studies will be required to determine whether the evolution of NS3 genetic diversity is related to virologic, immunologic, or clinical progression and HIV and/or HCV treatment responses.
Nucleotide sequence accession numbers. Sequences reported in this paper have been submitted to GenBank under accession numbers DQ355522 to DQ355749.
|
|
|---|
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»