JOURNAL OF CLINICAL MICROBIOLOGY, Jan. 1998, p. 227–233 0095-1137/98/$04.0010 Copyright © 1998, American Society for Microbiology
Vol. 36, No. 1
Hepatitis C Virus Heteroduplex Tracking Assay for Genotype Determination Reveals Diverging Genotype 2 Isolates in Italian Hemodialysis Patients PIER LUIGI CALVO,1 JOE KANSOPON,2 KULDIP SRA,2 STELLA QUAN,2 ROBERT DINELLO,2 ROBERTO GUASCHINO,3 GIOVANNI CALABRESE,4 FRANCA DANIELLE,1 MAUIZIA ROSSANA BRUNETTO,5 FERRUCCIO BONINO,5 ANNA LUCIA MASSARO,1 ALAN POLITO,2 MICHAEL HOUGHTON,2 AND AMY J. WEINER2* Chiron Corporation, Emeryville, California,2 Blood Bank3 and Haemodialysis Unit,4 Casale Monferrato Hospital, Casale Monferrato, and Blood Bank1 and Department of Gastroenterology,5 Molinette Hospital, Turin, Italy Received 25 July 1997/Returned for modification 9 September 1997/Accepted 27 October 1997
A heteroduplex tracking assay (HTA) was developed for genetic analyses of the hepatitis C virus (HCV) using single-stranded probes from the core (C)/E1 region. Nucleotide sequencing of reverse transcriptase (RT)-PCR products from 15 Italian dialysis patients confirmed the specificity and accuracy of the HTA genotyping method, which identified 5 of 15 (33.3%) 1b, 7 of 15 (46.7%) 3a, and 3 of 15 (20%) type 2 infections. The genotypes of an additional 12 HCV antibody-positive blood donors from different geographical locations were also in agreement with the genotypes determined by the Inno-LiPA HCV II kit (Innogenetics) and/or restriction fragment length polymorphism (RFLP). Isolates which had between 35 to 40% nucleotide divergence from control subtype 1a, 1b, 2a, 2b, or 3a standards could be typed. Surprisingly, HTA detected one 1b-2 coinfection which was missed by DNA sequencing. Three samples that were designated non-2a or 2b type 2 by HTA were found to be type 2a by both RFLP and direct nucleotide sequencing of the 5* untranslated region. The genetic distance between patient type 2 and control 2a, 2b, and 2c isolates indicated that a new subtype was present in the population being studied. Serotyping (RIBA serotyping strip immunoblot assay kit) of 23 dialysis patients showed that the genotype could be determined in 6 of 8 (75%) C/E1 RT-PCR-negative and 15 of 23 (65.2%) RT-PCR-positive samples, indicating that the two tests complement each other. Hepatitis C virus (HCV) is a single-stranded RNA virus (9) which is known to cause hepatitis and infects an estimated 100 million people worldwide (3). The HCV genome consists of 59 and 39 untranslated regions (UTR) that flank a single long open reading frame which encodes the structural (core [C], envelope 1 [E1], E2, and p7) and nonstructural (NS2, NS3, NS4a, NS4b, NS5a, and NS5b) proteins required for the manufacture of infectious virus (reviewed in references 21, 23, 38, and 43). Typical for RNA viruses, the mutation rate of HCV genomes is reported to be 1.44 3 1023 per genome per site per year (28), and variants circulate in a quasispecies distribution within infected individuals (26, 41). Classification of viral genomes into genotypes, subtypes, and isolates based on the degree of nucleotide heterogeneity of different segments of the viral genome has been proposed (reviewed in reference 4; 33), although the relationship between genetic polymorphism and biological features of the different genotypes has not been clearly established for HCV (reviewed in reference 2). Since genotype may be a factor in designing and administering vaccines and therapeutic agents, reliable assays for hepatitis C are needed. Since sequencing complete viral genomes is not practical, and many clinical laboratories do not have DNA-sequencing facilities, several methods have been described for determining HCV genotypes based on restriction enzyme-, PCR-, or hybridization-based analysis of the 59 UTR (11, 36), core (30, 34), or NS5 (16, 33) region. Current methods for subtyping have limited accuracy either because they assess short, relatively
conserved regions of the genome, such as the 59 UTR, or because they can yield false results due to nonspecificity of PCR primers (25, 30). Although hybridization-based assays may be reasonably accurate (36), information about nonhybridizing samples can be obtained only by DNA sequencing. We describe the application of the hybridization-based heteroduplex tracking assay (HTA) (12–14) to HCV for genotyping and identifying molecular variants which might be missed by other genotyping methods. HTA involves hybridizing probes from known HCV subtypes to reverse transcriptase PCR (RT-PCR) products from the homologous sera or from unknown samples and electrophoresing the hybridization products on gels. The formation of a heteroduplex band on a gel indicates genotype and subtype. In the absence of genetic rearrangements, insertions, and/or deletions, the genetic relationship between isolates can be determined by the relative migration of the heteroduplexes on gels (12). To develop and evaluate an HTA genotyping system for HCV, three methods of typing isolates (restriction fragment length polymorphism [RFLP], DNA sequencing, and RIBA serotyping strip immunoblot assay [SIA]) were performed on 23 well-preserved and carefully characterized serum samples from a cohort of Italian dialysis patients (5, 18, 39). Samples from an additional 12 blood donors from different countries (the United States, Egypt, and Holland) were also analyzed by HTA, and the results were compared to those obtained by the InnoLiPA HCV II test (Innogenetics) and/or RFLP to establish the accuracy of HTA as a method for determining genotype.
* Corresponding author. Mailing address: Chiron Corporation, 4560 Horton St., Emeryville, CA 94608. Phone: (510) 420-4785. Fax: (510) 658-0329. E-mail:
[email protected].
Patient population. Twenty-six of 35 Italian patients who had been undergoing regular hemodialysis in the same unit since 1991 (5, 18, 39) tested positive for
MATERIALS AND METHODS
227
228
CALVO ET AL.
J. CLIN. MICROBIOL.
TABLE 1. HTA RT-PCR primers Primer
Sequence
Nucleotide coordinates (59339)a
PCR round
C170Sb E338A1 E338A2a E338A2b E338A3a E320Ac C179S1 C179S2 C179S3
CCTGGTTGCTCTTTCTCTATCT GATGGCTTGTGGGATCCGGAG GATGACCTCGGGGACGCGCAT GACCAGTTCTGGAACACGAGC CAAGGTCTGGGGTAAACGCAG CCAGTTCATCATCATATCCCA TGGCCCTGCTCTCTTGCTTGAC TTGCTCTTCTGTCGTGCGTCAC TTGCTCTGTTCTCTTGCTTAAT
508–529 1032–1012 1032–1012 1032–1012 1032–1012 978–958 536–557 536–557 536–557
I I I I I II II II II
a b c
According to Choo et al. (10). Universal sense primer for genotype determination. Universal antisense primer for genotype determination.
anti-HCV antibodies in August 1995 by the second- and third-generation Ortho HCV enzyme-linked immunosorbent assay (ELISA). Serum samples were collected, divided into aliquots, and stored at 280°C. Twenty-five patients were both Ortho ELISA 3.0 and RIBA 3.0 HCV SIA reactive, while one patient was indeterminate (C-22-P reactive). Twenty-three ELISA 3.0 and RIBA 3.0 SIAreactive serum samples were available for this study. RT-PCR was performed on 15 of 23 samples with primers from both the 59 UTR (modified from Han et al., [19] sense primer 93B, 59ACCATGAATCACTCCCCT39 and antisense primer JH51, 59CCCAACACTACTGGGCTA39) and C/E1 (described in Table 1). One, nine, and two RT-PCR (59 UTR and C/E1)-positive blood donor serum samples from Egypt (a gift from J. Lau, University of Gainesville), the United States (Sacramento Blood Bank), and Holland (a gift from H. Reesink, C. van der Poel, and N. Lelie; CLB, Amsterdam, Holland), respectively, were analyzed by HTA. RNA extraction, cDNA synthesis, and PCR. RNA was extracted from 100 ml of sera or plasma on at least two different occasions with the Stratagene RNA isolation kit (no. 200345) according to the manufacturer’s instructions (8). A 25-ml cDNA reaction mixture contained 20% plasma RNA extract, 10 mM Tris-HCl (pH 8.3), 10 mM MgCl2, 10 mM dithiothreitol, 75 mM KCl, 1 mM (each) deoxynucleoside triphosphates, 25 U of RNasin (Promega), 5 U of avian myeloblastosis virus RT (no. 80855B; Gibco BRL), and 100 pmol of antisense primer (Table 1). The cDNA mixture was diluted with 50 ml of water, boiled for 5 min, quick-cooled on ice, and added to the PCR reagents, with final concentrations according to the Perkin PCR kit (no. N801-0055) instructions. cDNA primer E338A1, -2a, -2b, or -3a and a universal primer (C170S) were used in the first round of PCR (PCR I) to amplify subtypes 1 (a and b), 2a, 2b, and 3a HCV, respectively. One to ten microliters of the PCR I product was added to a second PCR (PCR II) mixture containing 100 pmol of universal antisense primer (320A) and type- or subtype-specific sense primers (C179S1, 2a, 2b, and 3a). All PCRs consisted of 40 cycles (94°C for 10 s, 55°C for 30 s, and 72°C for 30 s) performed in a Perkin-Elmer 9600 thermocycler, with the exception of type 1b PCR II, in which the annealing temperature was lowered to 45°C. Single-stranded probes. Initially, subtype-specific probes were obtained from RT-PCR products generated from RNA extracted from individuals infected with HCV of known genotype (1a, 1b, 2a, 2b, and 3a) as described above, with the exception of primer 320A, which was biotinylated. Subsequently, RT-PCR products used to synthesize the probes used in Fig. 1 were subcloned in the PCR 2.1 vector (Invitrogen) and used as templates for the biotinylated PCR products used to generate the single-stranded probes shown in Fig. 3 (lanes 1a, 1b, 2a, 2b, and 3a). Approximately 500 ng of biotinylated PCR II product was bound to 20 ml (200 mg) of streptavidin-coated magnetic beads (no. 112.05; Dynal M280) which were first washed in 100 ml of phosphate-buffered saline, pH 7.4–0.1% bovine serum albumin, followed by 100 ml of binding and 13 washing (B&W) buffer containing 5 mM Tris-HCl (pH 7.5), 0.5 mM EDTA, and 1.0 M NaCl and finally resuspended in 50 ml of 23 B&W buffer. After being gently mixed at room temperature for approximately 15 min, the beads were collected in a magnetic field, the supernatant was removed, and the beads were washed three times with 200 ml of 13 B&W buffer to eliminate unbound PCR products, primers, and deoxynucleoside triphosphates. Eight microliters of freshly prepared 0.1 M NaOH was added to the beads for 10 min at room temperature to denature the bound, double-stranded PCR products. Beads were again collected in a magnetic field, and the supernatant containing single-stranded DNA was removed and neutralized with 40 ml of water, 4 ml of 0.2 M HCl, and 1.0 ml of 1.0 M Tris-HCl, pH 8.0 (11a). Twenty nanograms of single-stranded DNA was labeled with 100 mCi of [g-32P]ATP (3,000 Ci/mmol) with T4 kinase and passed over a G-50 column (no. 27.5335.01; Pharmacia) to remove excess [g-32P]ATP. Biotinylated primers and single-stranded probes were stored at 220°C in aliquots for optimal results. HTA. One to two microliters (0.5 ng) of single-stranded DNA probe (specific activity, 1 3 109 cpm/mg) was hybridized to a 100-fold molar excess of colinear,
double-stranded PCR product from either the homologous standard (JK1a, -1b, -2a, -2b, or -3a) or patient sera in the following 10-ml reaction mixture: 4 ml of double-stranded PCR product (;50 ng), 4 ml of 23 SSC (13 SSC is 0.15 M NaCl plus 0.015 mM sodium citrate, 10 mM EDTA [pH 8.0]) and 2 ml of water. The reaction was denatured for 3 min at 95°C and annealed at 55°C for 2 h in a heat block. The entire sample was loaded onto a 1-mm-thick, 6% Hydrolink MDE gel (no. 4739-00; Baker) and electrophoresed for 16 h at 500 V. Gels were vacuum dried at 80°C on filter paper and exposed to X-ray film for 15 min. The relative migration of a heteroduplex was calculated as the distance (in millimeters) the heteroduplex migrated from the top of the gel divided by the distance (in millimeters) the homoduplex migrated from the top of the gel. DNA sequencing. PCR products and cloned DNA were purified by either Wizard PCR or Maxi prep (no. A7170 or A7270; Promega) (no. 12191; Qiagen) according to the manufacturer’s protocol and sequenced by standard methods. PCR products were cloned with the TA cloning vector (no. K2000-01; Invitrogen). A minimum of three clones were sequenced for each patient sample. Phylogenetic trees. The phylogenetic trees in Fig. 2 were constructed by pairwise, progressive alignment of the nucleotide sequences to one another with the computer software program Gene-Works Unweighted Pair Group Method with Arithmetic Mean (41). RFLP analysis of 5* UTR PCR products. cDNA amplified with 59 UTR primers (32) was digested as described by Davidson et al. (11) with ScrF1 to distinguish type 2 HCV subtypes. The MvaI/HinF1 fragment was used to determine genotype for the 12 blood donor samples. Digested DNA was electrophoresed on 6% polyacrylamide gels, stained with ethidium bromide, and visualized under UV light. There were no differences between the restriction fragment patterns of PCR products generated with primers described by Shimizu et al. (32) and those described by Davidson et al. (11), since the PCR primers are almost colinear. Inno-LiPA HCV II test. The Inno-LiPA HCV II test (Innogenetics), which is based on sequence heterogenetics in the 59 UTR, was performed with samples from blood donors from the United States according to the manufacturer’s protocol. RIBA HCV SIA. Eight serotype-specific HCV peptide antigens, five of which are from the NS4 region of HCV subtypes 1a and 1b, 2a and 2b, and 3 and three of which are from the core region of types 1 and 2, were immobilized on a nitrocellulose solid support. An algorithm was defined such that serotype was determined by which type-specific NS4 peptide was reactive in the tested sera. In the absence of NS4 peptide reactivity, core peptide reactivity was considered for the interpretation of the results (15).
RESULTS Specificity of HTA genotyping. Single-stranded probes (nucleotides 536 to 978), derived from sera of known subtype, were hybridized to colinear RT-PCR products from 15 Italian dialysis patients and electrophoresed on MDE gels. As shown in Fig. 1A through E, each probe (lanes 1a, 1b, 2a, 2b, and 3a, respectively) hybridized to RT-PCR products from the homologous sera, forming a homoduplex (lane h). Heteroduplexes were observed for five samples hybridized to the JK1b probe (Fig. 1B, samples 1, 2, 4, 24, and 28), and seven samples hybridized to the JK3a probe (Fig. 1E, samples 3, 7, 18, 20, 22, 26, and 35), while there was no subtype 1a HCV (Fig. 1A). Probe-only lanes (1a, 1b, 2a, 2b, and 3a in Fig. 1 and p in Fig. 3) show the level of background, which can be observed in all lanes, when the probe is degrading. The specificity of HTA probes for their respective subtypes was demonstrated by the lack of cross-hybridization of the subtype-specific probes (JK1a, JK1b, JK2a, JK2b, and JK3a) to each other or to PCR products from more than one patient. Two exceptions were patient 4, who was most likely coinfected with 1b and type 2 HCV (Fig. 1B and D), and blood donor k (Fig. 3A and B), who had both 1a and 1b HCV. The mean relative migration of subtype 1b isolates (0.79; range, 0.72 to 0.83) was less than and did not overlap with that of the 3a isolates (0.89; range, 0.86 to 0.93), indicating that the 1b isolates were less closely related to the 1b probe than the 3a isolates were to the 3a probe. Surprisingly, we found three samples which formed slowly migrating heteroduplexes, two of which hybridized to the JK2a and JK2b probes (samples 30 and 33 in Fig. 1C and D) and one which hybridized to the JK2a probe only (Fig. 1C, sample 23). Samples 23, 30, and 33 (HTA 23, 30, and 33) had relative migrations of 0.45, 0.1, and 0.55 to
VOL. 36, 1998
GENOTYPING HCV BY A HETERODUPLEX TRACKING ASSAY
229
FIG. 1. Single-stranded probes (lanes 1a, 1b, 2a, 2b, and 3a in panels A to E, respectively) were hybridized to PCR products from either the same control serum samples from which they were derived, forming homoduplexes (h), or from dialysis patient sera, forming heteroduplexes (numbered lanes), and were electrophoresed on MDE gels.
the JK2a probe, while HTA 30 and 33 had relative migrations of 0.12 and 0.1 with respect to the JK2b probe. The data indicated that HTA 23, 30, and 33 were substantially divergent from the JK2a and/or JK2b probes and suggested that they were most likely a non-2a, non-2b type 2 subtype. DNA sequencing verifies HTA genotype and subtype assignments and confirms complexity within type 2. The genotypes of all 15 PCR-positive Italian patients and subtype probes were verified by nucleotide sequencing (data not shown) of the partial C/E1 PCR products by direct sequencing and/or by sequencing cloned PCR products. Phylogenetic trees shown in Fig. 2A through C, which were derived from the nucleotide sequences, confirmed that isolates which were designated 1b (HTA 1, 2, 4, 24, and 28) or 3a (HTA 3, 7, 18, 20, 22, 26, and 35) segregated with the appropriate published 1b (HCV-J [22]) and 3a (NZ1-1 [31]) prototype sequences. The mean percentage of nucleotide divergence of 1b and 3a patient isolates from probes JK1a and JK3a of 10.5% (range, 8.7 to 11.6%) and 4.3% (range, 1.4 to 6.5%) was in agreement with the HTA prediction that the 3a isolates were related more closely to the JK3a probe than the 1b isolates were to the JK1b probe, based on the relative migration of the 1b and 3a heteroduplexes on gels. The C/E1 nucleotide sequence data indicated that (i) isolates with $1.4% nucleotide substitutions can be separated by HTA analysis, (ii) isolates within the subtypes 1b and 3a
were ,15% divergent, and (iii) HTA subtype probes do not cross-hybridize with isolates with greater than approximately 35 to 40% nucleotide divergence. Since HTA 23, 30, and 33 were 22.8 to 27% divergent from JK2a, and HTA 30 and 33 were 33 to 34.5% divergent from JK2b, we compared the nucleotide sequences of HTA 23, 30, and 33 to a published 2c isolate (S83 [4]). HTA 23, 30, and 33 were 16.5, 14.7, and 18.6% divergent from S83 and segregated with S83 in Fig. 2B as expected. Interestingly, HTA 33, which was most divergent from S83, was more closely related to HTA 23 and 30 (11.7 and 15.6% nucleotide heterogeneity) than to either JK2a or JK2b, suggesting a greater complexity among type 2 isolates than among type 1 or type 3 isolates. Although most patients had a single virus isolate or a tight cluster of virus isolates, which were not resolved on gels when hybridized to their respective probes, HTA 18 and 22 each had two bands which were not observed in the homoduplex or single-stranded probe control lane (Fig. 1E, lanes h and 3a). The data suggested that these patients each had at least two isolates with greater than 1.4% nucleotide heterogeneity. It is unlikely that the double bands seen in all lanes except one in Fig. 1B resulted from double-stranded probe molecules contaminating the single-stranded probe, since the result was reproducible with different preparations of probe and was not
230
CALVO ET AL.
J. CLIN. MICROBIOL.
FIG. 2. Phylogenetic analysis of partial C/E1, E1, or 59 UTR nucleotide sequences from dialysis patients (HTA 1 to 4, 7, 18, 20 to 26, 28, 30, 33, and 35) and published subtypes 1a (HCV-1 [10]), 1b (HCV-J [22]), 2a (HC-J6 [29]), 2b (HC-J8 [29]), 2c (S83 [4]), and 3a (NZL-1 [31]). Nucleotide coordinates are according to reference 12.
observed for the other probes generated under identical conditions. Comparison of partial C/E1 HTA and RFLP. We compared the genotype results obtained by partial C/E1 HTA with RFLP (59 UTR) and direct nucleotide sequencing of 59 UTR PCR products from all 15 patients. Subtype assignments were all in agreement except for HTA 23, 30, and 33, which were incorrectly identified as subtype 2a by RFLP with ScrFI (Table 2). A dendrogram of the 59 UTR nucleotide sequences of HTA 23, 30, and 33 and published sequences (Fig. 2D) showed HTA 23, 30, and 33 as 2a, in agreement with RFLP results, but was inconsistent with the same analysis performed on the E1 sequence (Fig. 2B). The 59 UTR analysis also inaccurately located 2c distant from subtypes 2a and 2b, a result which is most likely due to the small number of nucleotide substitutions in the 59 UTR. Serotyping analysis of patient sera. The serotype of 15 of 23 (65.2%) patient serum samples could be determined by the RIBA HCV Serotyping SIA (Table 2). Four 59 UTR-, C/E1-, RT-PCR-, and ELISA 3.0-positive samples, which were subtyped by HTA, were nonreactive. Seventy-five percent (6 of 8) of PCR-negative samples were assigned a serotype (Table 2), while 1 of 8 was nonreactive and 1 of 8 was untypeable. Two subtype 3a samples (HTA 20 and 35) were untypeable (designated either type 1 or 3). One sample (HTA 18) was a subtype 3a by HTA, RFLP, and DNA sequencing but was designated type 2 by serotyping. The discrepant sample may be attributed to a past infection or to a coinfection at a time when viremia was undetectable by PCR. Eight of 15 RT-PCR (53.3%) sam-
ples were concordant among all three assays (HTA, RFLP, and serotyping) with respect to genotype. Comparison of genotype determined by HTA to those determined by Inno-LiPA HCV II and/or RFLP in blood donors from different geographical locations. HCV genotypes determined by HTA were in agreement with those determined by RFLP and the Inno-LiPA HCV II test for all 12 samples (Table 3), irrespective of the country from which the serum or probe was obtained. For example, the 1b probe, which came from a Japanese serum sample, hybridized to 1b sera from the United States, Holland, and Italy (Fig. 3B, lanes e and k and Fig. 1B, lanes 1, 2, 4, 24, and 28), and the type 3a probe (of U.S. origin) hybridized to 3a samples from the United States, Egypt, and Italy (Fig. 3E, lanes a, b, and g and Fig. 1E, lanes 3, 7, 18, 20, 22, 26, and 35). The data shown in Fig. 3 indicated that the 1a and 1b probes failed to cross-hybridize to each other (panels A and B, lanes 1a and 1b), and therefore sample k, which hybridized to both probes, must have had a 1a-1b coinfection. HTA also showed that samples c (panel A), e (panel B), g and h (panel C), and d (panel D) had at least one viral variant. The cloned 2b probe in Fig. 3D was derived from blood donor serum sample f, and therefore a homoduplex can be seen at the bottom of the gel. DISCUSSION Based on the work of Delwart et al., who previously established the validity of HTA for determining genotype and molecular evolution studies on human immunodeficiency virus
VOL. 36, 1998
GENOTYPING HCV BY A HETERODUPLEX TRACKING ASSAY
231
TABLE 2. Comparison of serotyping, RFLP and HTA genotype analysis of sera from dialysis patients HTA C/E1 sample PCR no.
1 2 3 4 5 6 7 18 19 20 22 23 24 26 27 28 29 30 31 32 33 34 35
1 1 1 1 2 2 1 1 2 1 1 1 1 1 2 1 2 1 2 2 1 2 1
Results according to RIBA HCV SIAe NS4 subtype 1
2
3
1111 6 2 2 111 1111 2 2 11 6 2 2 1111 2 6 1111 11 2 111 2 2 2 2
2 2 2 2 2 6 2 2 2 2 2 6 2 6 2 2 2 6 2 2 2 11 2
6 2 2 2 6 2 6 2 2 2 11 2 11 1111 2 11 6 2 2 2 2 11 6
Core 1 2(I1J)
111 2 6 2 6 1111 2 11 6 1111 2 2 1111 11 6 6 2 6 2 2 6 2 11
2 2 2 2 2 6 2 111 11 6 2 1 6 11 6 2 2 1111 11 111 11 2 6
Genotype
1 NRb NR NR 1 1 NR 2 1 1–3c 3 2 1 3 NR 1 1 2 1 2 2 UNT(2)c 1–3c
Subtype bya: RFLP HTA
1b 1b 3a 1b ND ND 3a 3a ND 3a 3a 2a 1b 3a ND 1b ND 2a ND ND 2a ND 3a
1b 1b 3a 1b ND ND 3a 3a ND 3a 3a 2d 1b 3a ND 1b ND 2d ND ND 2d ND 3a
a
ND, not done. NR, nonreactive. UNT, untypeable. d Non-2a, non-2b type 2. e 2, negative; 6, positive/negative; 1 to 1111, positive. b c
type 1 (12–14), Gretch et al. applied HTA to the study of HCV quasispecies (42). Unfortunately, the primers used in quasispecies analysis were not appropriate for genotype determination due to a high degree of nucleotide variation in the amino terminus of E2 (20, 40). In order to develop an HTA genotyping method for HCV which would also be suitable for other purposes, RT-PCR primers were chosen from the carboxy terminus of core and the middle two-thirds of E1. A set of universal PCR primers were selected from highly conserved regions of C/E1, along with a set of primers which were derived
TABLE 3. Comparison of Inno LiPA HCV II (IL-II), RFLP, and HTA genotype analyses of blood donor samples Sample
a b c d e f g h i j k l a b
Country of origin
Egypt United States United States United States United States United States United States United States United States United States Holland Holland
Subtype not determined. NT, not tested.
Genotype determined by: IL/II
NT 3a 1a 2b 1b 2b 3a 2a 2a 1a NT NT
b
RFLPa
HTA
3 3 1 2 1 2 3 2 2 1 1 1
3a 3a 1a 2b 1b 2b 3a 2a 2a 1a 1a, 1b 1a
FIG. 3. Single-stranded probes (1a, 1b, 2a, 2b, and 3a) were hybridized to either PCR products from which each single-stranded probe was derived (lanes 1a, 1b, 2a, 2b, 3a in panels A to E) or to HCV antibody-positive sera from the United States (lanes b to j), Egypt (lane a), and Holland (lanes k and l). The single-stranded probe alone is designated by the letter p.
from regions of substantial sequence heterogeneity to ensure subtype specificity. DNA sequencing of the partial C/E1 PCR products from a cohort of Italian dialysis patients confirmed the accuracy of HTA genotype assignment in all 15 individuals but failed to identify type 2 sequences in a 1b-2 coinfected patient (HTA 4). Since direct DNA sequencing or sequencing a relatively small number of cloned PCR products, i.e., 3 to 30, identifies only predominant viral species, sensitive hybridization-based techniques are more capable of detecting low-frequency variants in a population. RFLP and Inno-LiPA HCV II verified C/E1 HTA genotyping of 12 type 1, 2, and 3 blood donors from different geographical locations, indicating that the source of the isolates or probes did not influence the accuracy of HTA. In addition, HTA detected a 1a-1b coinfection in sample k. The inability of certain assays to detect coinfections could be attributed to either the relative abundance and/or efficiency of amplifying different genomes in the sample. We found that single-stranded probes derived from sera containing a virus population of less than 1.4% nucleotide divergence were effective for genotyping; however, the use of single-stranded probes from cloned RT-PCR products guarantees that the pattern of heteroduplex molecules on MDE gels reflects the actual population of variants and avoids potential artifacts
232
CALVO ET AL.
which may be more difficult to discern when double-stranded HTA probes (12–14, 42) or single-strand conformation polymorphism (17) is used. The molecular definitions of genotype, subtype, and isolate were originally based on the nucleotide region of NS5 for types 1 to 3 HCV (33) (the most frequent genotypes found in North America, Europe, and parts of Asia) and was later found to be consistent for other segments of the HCV genome, including the E1 gene (4, 34, 37). For full-length E1, subtypes and genotypes were separated by approximately 24 to 31% and 33 to 46% nucleotide heterogeneity, respectively, while isolates within the same subtype were found to be #15% divergent (4). Recently, sequences of new isolates in Southeast Asia showed that these isolates were more broadly distributed with respect to percent nucleotide divergence than types 1 to 3 HCV and that the boundaries between type, subtype, and isolate might need to be reevaluated in these populations (1, 27, 35). The interpretation of HTA data for Southeast Asian isolates would be expected to be complicated by the same factors that have affected the genotype assignment of these isolates by direct nucleotide sequencing. Our partial C/E1 nucleotide sequencing data for 1b and 3a HCV showed that isolates within the same subtype were #15% divergent and had approximately 35 to 42% nucleotide heterogeneity between isolates of different genotypes, which is in good agreement with the results of Bukh et al. (4). Consistent with previous studies showing that the relative migrations of the heteroduplexes are proportional to genetic distance (12), we found that the mean and range of the relative migration of 1b and 3a isolates were nonoverlapping and reflected the percentage of nucleotide divergence from their respective probes. Although the relative migration of individual isolates was not necessarily directly proportional to the percentage of nucleotide divergence from the probe, the slight differences in the observed relative migration of individual isolates may be due to the distribution of nucleotide substitutions in PCR fragments. In agreement with Delwart et al. (12), a minimum of 1.4% nucleotide divergence was sufficient to separate viral isolates in the HTA described in this study. The substantially retarded migration of HTA 23, 30, and 33 in Fig. 1C and D suggested that they were non-2a, non-2b type 2 subtypes. These samples would have either been scored negative (6) or misassigned as either 2a or 2a-2b coinfections if standard PCR-based genotyping methods had been used, depending on the stringency of the assay. DNA sequence comparisons showed that HTA 23, 30, and 33 were 14.7 to 18.6% divergent from S83 (a 2c subtype known to be prevalent in Italy [6]). These results indicate that HTA 23, 30, and 33 (i) approached or exceeded the expected ;15% divergence between isolates within the same subtype and (ii) were more closely related to each other than to other known type 2 isolates which have been sequenced in the C/E1 region. Interestingly, based on a calculation of the percentage of nucleotide divergence alone, HTA 23, 30, and 33 would be considered borderline type 2 compared to 2b isolates (33 to 34%), while they would clearly be classified as type 2 compared to type 2a or 2c (22.8 to 27.9% versus to 16.5 to 18.6%) HCV. Since percentage of nucleotide divergence is closely related to genetic distance as determined by phylogenetic analysis for types 1, 2, and 3 HCV (33, 34), we demonstrate that, as for HIV-1, HTA is useful for identifying evolutionary variants which can be sequenced and subjected to more sophisticated phylogenetic analyses. These findings also illustrate that viral evolution is a dynamic process and highlight the ambiguities inherent in the classification of HCV strictly on the basis of molecular sequence criteria. A comparison of several genotyping methods demonstrated
J. CLIN. MICROBIOL.
that the partial C/E1 HTA was the most comprehensive for the following reasons: (i) not only was subtype for patients and blood donors accurately assigned, but coinfections, which were missed by DNA sequencing and RFLP, were also identified; and (ii) in contrast to HTA, RFLP and direct sequencing of the 59 UTR incorrectly identified all three type 2 samples (HTA 23, 30, and 33) as either type 2a or as 2a-2b due to inadequate nucleotide heterogeneity in the 59 UTR. HTA has a further advantage in that this method visualizes intrasubtype variants, which tests such as RFLP and Inno-LiPa HCV II are not designed to detect. Such variants could be important in understanding and characterizing therapeutic drug escape mutants. We also conclude that since 65.2% (15 of 23) of total and 75% (6 of 8) of PCR-negative dialysis patients were typeable by the RIBA HCV serotyping SIA, serotyping appears to have independent value in typing HCV and complements nucleic acid-based genotyping assays. ACKNOWLEDGMENTS We thank David Gretch, Jeff Wilson, Steve Polyak, and Kevin Crawford for technical assistance and helpful discussions. DNA sequencing was expertly performed by Chun Ting Lee. We also thank George Lamson for computer-generated dendrograms and Peter Anderson for preparation of the manuscript. This work was supported by Chiron Vaccines. REFERENCES 1. Adams, N. J., R. W. Chamberlain, L. A. Taylor, F. Davidson, C. K. Lin, R. M. Elliott, and P. Simmonds. 1997. Complete coding sequence of hepatitis C virus genotype 6a. Biochem. Biophys. Biophys. Res. Commun. 234:393–396. 2. Alter, H. J. 1995. To C or not to C, these are the questions. Blood 85:1681– 1695. 3. Alter, M. J. 1995. Epidemiology of hepatitis C in the West. Semin. Liver Dis. 15:5–14. 4. Bukh, J., R. H. Miller, and R. H. Purcell. 1995. Genetic heterogeneity of hepatitis C virus: quasispecies and genotypes. Semin. Liv. Dis. 15:41–63. 5. Calabrese, G., G. Vagelli, R. Guaschino, and M. Gomella. 1991. Transmission of anti-HCV within the household of haemodialysis patients. Lancet 337:1466. 6. Cammarota, G., F. Maggi, M. L. Vatteroni, L. Da Prato, L. Barsanti, M. Bendinelli, and M. Pistello. 1995. Partial nucleotide sequencing of six subtype 2c hepatitis C viruses detected in Italy. J. Clin. Microbiol. 33:2781–2784. 7. Chan, S.-W., F. McOmish, E. C. Holmes, B. Dow, J. F. Peutherer, E. Follett, P. L. Yap, and P. Simmonds. 1992. Analysis of a new hepatitis C virus type and its phylogenetic relationship to existing variants. J. Gen. Virol. 73:1131– 1141. 8. Chomczynski, P., and N. Sacchi. 1986. Single-step method of RNA isolation by acid guanidinium thiocyanate-phenol-chloroform extraction. Anal. Biochem. 162:156–159. 9. Choo, Q.-L., G. Kuo, A. J. Weiner, L. R. Overby, D. W. Bradley, and M. Houghton. 1989. Isolation of a cDNA clone derived from a blood-borne non-A, non-B viral hepatitis genome. Science 244:359–362. 10. Choo, Q.-L., K. H. Richman, J. H. Han, K. Berger, C. Lee, C. Dong, C. Gallegos, D. Coit, A. Medina-Selby, P. J. Barr, A. J. Weiner, D. W. Bradley, G. Kuo, and M. Houghton. 1991. Genetic organization and diversity of the hepatitis C virus. Proc. Natl. Acad. Sci. USA 88:2451–2455. 11. Davidson, E., P. Simmonds, J. C. Ferguson, L. M. Jarvis, B. C. Dow, E. A. C. Follett, C. R. Seed, T. Krusius, C. Lin, G. A. Medgyesi, H. Kiyokawa, G. Olim, G. Duraisamy, T. Cuypers, A. A. Saeed, D. Teo, J. Conradie, M. C. Kew, M. Lin, C. Nuchaprayoon, O. K. Ndimbie, and P. L. Yap. 1995. Survey of major genotypes and subtypes of hepatitis C virus using RFLP of sequences amplified from the 59 non-coding region. J. Gen. Virol. 76:1197– 1204. 11a.Delwart, E. Personal communication. 12. Delwart, E. L., E. G. Shpaer, J. Louwagie, F. E. McCutchan, M. Grez, H. Ru ¨bsamen-Waigmann, and J. I. Mullins. 1993. Genetic relationships determined by a DNA heteroduplex mobility assay, Analysis of HIV-1 env genes. Science 262:1257–1261. 13. Delwart, E. L., H. W. Sheppard, B. D. Walker, J. Goudsmit, and J. I. Mullins. 1994. Human immunodeficiency virus type 1 evolution in vivo tracked by DNA heteroduplex mobility assays. J. Virol. 68:6672–6683. 14. Delwart, E. L., M. P. Busch, M. L. Kalish, J. W. Mosley, and J. I. Mullins. 1995. Rapid molecular epidemiology of HIV-1 transmission. AIDS Res. 11:1181–1193. 15. Dixit, V., S. Quan, P. Martin, D. Larson, M. Brezina, R. DiNello, K. Sra, J. Y. N. Lau, D. Chien, J. Kolberg, A. Tagger, G. Davis, A. Polito, and G.
VOL. 36, 1998
16. 17.
18. 19.
20. 21. 22.
23. 24.
25.
26.
27.
28.
29.
GENOTYPING HCV BY A HETERODUPLEX TRACKING ASSAY
Gitnick. 1995. Evaluation of a novel serotyping system for hepatitis C virus: strong correlation with standard genotyping methodologies. J. Clin. Microbiol. 33:2978–2983. Enomoto, N., A. Takada, T. Nakao, and T. Date. 1990. There are two major types of hepatitis C virus in Japan. Biochem. Biophys. Res. Comm. 170:1021– 1025. Enomoto, N., M. Kurosaki, Y. Tanaka, F. Marumo, and C. Sato. 1994. Fluctuation of hepatitis C virus quasispecies in persistent infection and interferon treatment revealed by single-strand conformation polymorphism analysis. J. Gen. Virol. 75:1361–1369. Guaschino, R., G. Calabrese, F. Bonelli, L. Calvo, A. Lazzaro, G. Vagelli, G., and S. Griva. 1994. Should HCV RNA detection affect the isolation policy in dialysis units? Nephrol. Dial. Transplant. 7:1070. Han, J. H., V. Shyamala, K. H. Richman, M. J. Brauer, B. Irvine, M. S. Urdea, P. Tekamp-Olson, G. Kuo, Q.-L. Choo, and M. Houghton. 1991. Characterization of the terminal regions of hepatitis C viral RNA: identification of conserved sequences in the 59 untranslated region and poly(A) tails at the 39 end. Proc. Natl. Acad. Sci. USA 88:1711–1715. Hijikata, M., N. Kato, Y. Ootsuyama, M. Nakagawa, S. Ohkoshi, and K. Shimotohno. 1991. Hypervariable regions in the putative glycoprotein of hepatitis C virus. Biochem. Biophys. Res. Commun. 175:220–228. Houghton, M. 1996. Hepatitis C viruses, p. 1035–1057. In B. N. Fields, D. M. Knipe, P. M. Howley, et al. (ed.), Virology. Lippincott-Raven, Philadelphia, Pa. Kato, N., M. Hijikata, Y. Ootsuyama, M. Nakagawa, S. Ohkoshi, T. Sugimura, and K. Shimotohno. 1990. Molecular cloning of the human hepatitis C virus genome from Japanese patients with non-A, non-B hepatitis. Proc. Natl. Acad. Sci. USA 87:9524–9528. Kolykhalov, A. A., S. M. Feinstone, and C. M. Rice. 1996. Identification of a highly conserved sequence element at the 39 terminus of hepatitis C virus genome RNA. J. Virol. 70:3363–3371. Kostrikis, L. G., E. Bagdades, Y. Cao, L. Zhang, D. Dimitriou, and D. D. Ho. 1995. Genetic analysis of human immunodeficiency virus type 1 strains from patients in Cyprus: identification of a new subtype designated subtype I. J. Virol. 69:6122–6130. Lau, J. Y. N., M. Mizokami, J. A. Kolberg, G. L. Davis, L. E. Prescott, T. Ohno, R. P. Perrillo, K. L. Lindsay, R. G. Gish, K. P. Qian, M. Kohara, P. Simmonds, P. and M. S. Urdea. 1995. Application of six hepatitis C virus genotyping systems to sera from chronic hepatitis C patients in the United States. J. Infect. Dis. 181:281–289. Martell, M., J. I. Esteban, J. Quer, J. Genesca `, A. Weiner, R. Esteban, J. Guardia, and J. Go ´mez. 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. Mellor, J., E. A. Walsh, L. E. Prescott, L. M. Jarvis, F. Davidson, P. L. Yap, P. Simmonds, and the International HCV Collaborative Study Group. 1996. Survey of type 6 group variants of hepatitis C virus in Southeast Asia by using a core-based genotyping assay. J. Clin. Microbiol. 34:417–423. Okamoto, H., M. Kojima, S. Okada, H. Yoshizawa, H. Iizuka, T. Tanaka, E. E. Muchmore, D. A. Peterson, Y. Ity, and S. Mishiro. 1992. Genetic drift of hepatitis C virus during an 8.2-year infection in a chimpanzee: variability and stability. Virology 190:894–899. Okamoto, H., K. Kurai, S. I. Okada, K. Yamamoto, H. Lizuka, T. Tanaka, S.
30. 31. 32. 33. 34.
35.
36.
37.
38. 39. 40.
41.
42.
43.
233
Fukada, F. Tsuda, and S. Mishiro. 1992. Full-length sequence of a hepatitis C virus genome having poor homology to reported isolates: comparative study of four distinct genotypes. Virology 188:331–341. Okamoto, T., Y. Sugiyama, and S. Okada. 1992. Typing hepatitis C virus by polymerase chain reaction with type-specific primers: application to clinical surveys and tracing infectious sources. J. Gen. Virol. 73:673–679. Sakamoto, M., Y. Akahane, F. Tsuda, T. Tanaka, D. G. Woodfield, and H. Okamoto. 1994. Entire nucleotide sequence and characterization of a hepatitis C virus of genotype V/3a. J. Gen. Virol. 75:1761–1768. Shimizu, Y. K., A. Iwamoto, M. Hijikata, R. H. Purcell, and H. Yoshikura. 1992. Evidence for in vitro replication of hepatitis C virus genome in a human T cell line. Proc. Natl. Acad. Sci. USA 89:5477–5481. Simmonds, P., E. C. Holmes, and T.-A. Cha. 1993. Classification of hepatitis C virus into six major genotypes and a series of subtypes by phylogenetic analysis of the NS-5 region. J. Gen. Virol. 74:2391–2399. Simmonds, P., D. B. Smith, F. McOmish, P. L. Yap, J. Kolberg, M. S. Urdea, and E. C. Holmes. 1994. Identification of genotypes of hepatitis C virus by sequence comparisons in the core, E1 and NS-5 regions. J. Gen. Virol. 75:1053–1061. Simmonds, P., J. Mellor, T. Sakuldamrongpanich, C. Nuchaprayoon, S. Tanprasert, E. C. Holmes, and D. B. Smith. 1996. Evolutionary analysis of variants of hepatitis C virus found in Southeast Asia: comparison with classifications based upon sequence similarity. J. Gen. Virol. 77:3013–3024. Stuyver, L., R. Rossau, A. Wyseur, M. Duhamel, B. Vanderborght, H. Van Heuverswyn, and G. Maertens. 1993. Typing of hepatitis C virus isolates and characterization of new subtypes using a line probe assay. J. Gen. Virol. 74:1093–1102. Stuyver, L., W. Van Arnhem, A. Wyseur, F. Hernandez, E. Delaporte, and G. Maertens. 1994. Classification of hepatitis C viruses based on phylogenetic analysis of the envelope 1 and nonstructural 5B regions and identification of five additional subtypes. Proc. Natl. Acad. Sci. USA 91:10134–10138. Tanaka, T., N. Kato, M.-J. Cho, K. Sugiyama, and K. Shimotohno. 1996. Structure of the 39 terminus of the hepatitis C virus genome. J. Virol. 70:3307–3312. Vagelli, G., G. Calabrese, R. Guaschino, and M. Gonella. 1992. Effect of HCV positive patients isolation on HCV infection incidence in a dialysis unit. Nephrol. Dial. Transplant. 7:1070. Weiner, A. J., M. J. Brauer, J. Rosenblatt, K. H. Richman, J. Tung, K. Crawford, F. Bonino, G. Saracco, Q.-L. Choo, M. Houghton, and J. H. Han. 1991. Variable and hypervariable domains are found in the regions of HCV corresponding to the flavivirus envelope and NS1 proteins and the pestivirus envelope glycoproteins. Virology 180:842–848. Weiner, A. J., M. M. Thaler, K. Crawford, K. Ching, J. Kansopon, D. Y. Chien, J. E. Hall, F. Hu, and M. Houghton. 1993. A unique, predominant hepatitis C virus variant found in an infant born to a mother with multiple variants. J. Virol. 67:4365–4368. Wilson, J. J., S. J. Polyak, T. D. Day, and D. R. Gretch. 1995. Characterization of simple and complex hepatitis C virus quasispecies by heteroduplex gel shift analysis: correlation with nucleotide sequencing. J. Gen. Virol. 76:1763– 1771. Yamada, N., K. Tanihara, A. Takada, T. Yorihuzi, M. Tsutsumi, H. Shimomura, T. Tsuji, and T. Date. 1996. Genetic organization and diversity of the 39 noncoding region of the hepatitis C virus genome. Virology 223:255–261.