Development of a sensitive RT-PCR method for amplifying and sequencing near full-length HCV genotype 1 RNA from patient samples
© Zhang et al.; licensee BioMed Central Ltd. 2013
Received: 19 September 2012
Accepted: 7 February 2013
Published: 12 February 2013
Direct-acting antiviral (DAAs) agents for hepatitis C virus (HCV) span a variety of targets, including proteins encoded by the NS3/4A, NS4B, NS5A, and NS5B genes. Treatment with DAAs has been shown to select variants with sequence changes in the HCV genome encoding amino acids that may confer resistance to the treatment. In order to assess these effects in patients, a Reverse Transcription Polymerase Chain Reaction (RT-PCR) method was developed to sequence these regions of HCV from patient plasma.
A method was developed to amplify and sequence genotype 1 HCV RNA from patient plasma. Optimization of HCV RNA isolation, cDNA synthesis, and nested PCR steps were performed. The optimization of HCV RNA isolation, design of RT-PCR primers, optimization of RT-PCR amplification conditions and reagents, and the evaluation of the RT-PCR method performance is described.
The optimized method is able to successfully, accurately, and reproducibly amplify near full-length genotype 1 HCV RNA containing a wide range of concentrations (103 to 108 IU/mL) with a success rate of 97%. The lower limit of detection was determined to be 1000 IU/mL HCV RNA.
This assay allows viral sequencing of all regions targeted by the most common DAAs currently in development, as well as the possibility to determine linkage between variants conferring resistance to multiple DAAs used in combination therapy.
KeywordsHepatitis C Virus (HCV) Reverse Transcription Polymerase Chain Reaction (RT-PCR) Genotype Direct Acting Antiviral Agents (DAA) Lower Limit of Detection (LLOD) Quasispecies Sequencing Resistance NS (Non-structural) NS3/4A NS4B NS5A NS5B
More than 170 million people (3% of the world population) are infected with hepatitis C virus (HCV) . Chronic HCV infection is a major cause of liver cirrhosis and hepatocellular carcinoma (HCC), as well as liver failure. Liver failure is one of the most common indications for liver transplantation in the United States [2, 3]. The recent development of potent, direct-acting antiviral agents (DAAs), used in combination with pegylated-interferon (Peg-IFN) and ribavirin (RBV) has significantly improved the sustained virologic response rates compared with Peg-IFN and RBV alone [4, 5].
HCV has high genetic variability across the six major genotypes as well as in different subtypes within each genotype. Genotype 1, the predominant cause of chronic HCV infections in the United States, Europe, and Japan, has more diversity within subtypes 1a and 1b than the entire diversity of HBV and HIV . In an HCV-infected patient, this pre-existing genetic variability provides a pool of variants that can adapt to new selection pressures, such as antiviral treatment. The variants are a natural part of the viral quasispecies, which arise from the error-prone nature of the HCV RNA-dependent RNA polymerase and the high replication rate of the virus . In patients who fail treatment with a potent DAA regimen, the elimination of the majority of the DAA-sensitive quasispecies can lead to selection and detection of DAA-resistant variants .
The nonstructural (NS) HCV proteins NS3, NS4A, NS4B, NS5A, and NS5B are essential for viral replication and are the major targets of DAAs currently in development . Thus, monitoring viral sequence changes across multiple NS regions is beneficial in understanding the potential association between therapeutic resistance and clinical outcome .
Various methods for amplification and sequence analysis of specific HCV genomes have been reported, and most of them have focused on individual NS segments or smaller amplicons [11–14]. In 2006, a long reverse transcription polymerase chain reaction (RT-PCR) amplification protocol of a near full-length HCV genome was developed by Fan, X, Xu, Y, and Di Bisceglie, AM, that amplified a limited number of patient samples [15–17]. While conducting telaprevir (INCIVEK™, INCIVO®) clinical trials [18, 19], the need arose to develop a method to analyze the HCV genotype 1 sequence of NS3, NS4A,NS5A, and NS5B from a single genome or the catalytic domain of the NS3·4A protease domain from a large number of patients. The challenge was that samples from patients are extremely variable and existing methods described in the literature required patient-by-patient optimization to efficiently amplify near full-length HCV RNA, particularly in samples with low HCV RNA levels, which were common in patients who were responding to telaprevir and/or VX-222. In this study, we describe the development of an efficient RT-PCR method that was used in telaprevir clinical trials to amplify a nearly full-length HCV genotype 1 genome from 7800 unique patient samples containing a wide range of HCV RNA concentrations (103 to 108 IU/mL) with a success rate of 97% [4, 5, 20–23]. This method enables HCV genes, targeted by the most common DAAs in clinical trials today, to be efficiently sequenced using clinical samples with low HCV RNA levels and permits the assessment of linkage between multiple amino acid changes in different HCV genes in samples from patients receiving DAA drug combinations.
Results and discussion
Optimization of HCV RNA extraction
The use of the 96-well format QIAamp Virus BioRobot 9604 Kit (QIAGEN, Valencia, CA) was optimized to produce high-quality viral RNA from HCV RNA-positive plasma. Several modifications were made, with the primary change involving the replacement of the carrier RNA in the RNA isolation. Comparison of the use of poly A RNA (rA) provided in the QIAGEN kit with transfer RNA (tRNA) (Roche Applied Science, Indianapolis, IN) or ribosomal RNA (rRNA) (Roche Applied Science, Indianapolis, IN) using plasma samples with a wide range of HCV RNA copy numbers (10 ~ 18,500 HCV RNA copies/reaction) showed that 60 μg of tRNA per column for 660 μL plasma led to the highest level of sensitivity in the subsequent amplification reaction. In addition, several other minor changes to the standard QIAamp protocol also increased viral RNA purity and quality, which helped the performance of subsequent amplifications including mixing via gentle rocking, instead of pipetting (per the BioRobot 9604 Kit protocol), to minimize shearing of the 9.6 kb HCV RNA. Also, for optimal recovery of viral RNA, the samples were lysed at 56°C and an additional spin step was added following the final wash step to eliminate residual ethanol carry-over.
Synthesis of cDNA is a critical step in the ability to amplify large PCR products, such as the 9.6 kb HCV genome. Optimization of the cDNA synthesis from HCV RNA involved testing different oligonucleotide primers including random hexamers, oligo d(A) (20 nucleotide), and primers that annealed in the NS5B or 3’ untranslated region (UTR). Surprisingly, most of the primers binding in the NS5B or 3’UTR did not allow amplification of full-length viral RNA or lead to very low cDNA yields. In contrast, oligo dA20, which anneals to the highly conserved poly (U) tail at the 3’ end of the HCV RNA, did allow for high sensitivity and reproducible amplification of full-length cDNA. We also tested several conditions and concentrations of Reverse Transcriptase enzyme (Superscript™ III Reverse Transcriptase [Life Technologies Corporation, Carlsbad, CA]) and it was determined that the use of 400 units with a 2.5 hour protocol, containing increasing annealing/elongation temperatures, was preferred.
Amplification of HCV via nested PCR
List of primers used for the amplification of the HCV coding region
Sensitivity of RT-PCR amplification
To assess the sensitivity of the method, four genotype 1 HCV RNA-positive plasma samples (two subtypes 1a and two subtypes 1b) were serially diluted with concentrations of HCV RNA ranging from 11–649,000 IU/mL, and two analysts performed the amplification on different days. The HCV RNA-positive plasma samples were processed through the RT-PCR step and subsequent PCR amplifications were performed in triplicate. Samples with HCV RNA greater than or equal to 1000 IU/mL had an amplification success rate of 70% (100 out of 143 samples) in the first attempt. Samples with HCV RNA below 1000 IU/mL had a much lower amplification success rate of 1.6% (2 out of 121 samples). Therefore, the lower limit of detection (LLOD) of the RT-PCR amplification method was established to be ~1000 IU/mL. No significant differences in success rates were observed either between subtype 1a and 1b samples or between analysts.
Robustness and success rate of amplification and sequencing
Summary of RT-PCR amplification and sequencing success rate of clinical samples from telaprevir and VX-222 clinical trials
HCV RNA (IU/mL)a
Numbers of clinical samples
Success rate of RT-PCR amplification % (n/N)
Success rate of sequencing % (n/N)
Comparison of HCV NS3(4A) region of Los Alamos and vertex sequences
Number of sequences°
Length of alignment (base pair)
Mean Shannon entropy
Average pairwise difference (Standard deviation)
Los Alamos Sequences
7.5% (± 1.9%)
9.1% (± 4.6%)
Los Alamos Sequences
9.3% (± 1.5%)
10.6% (± 3.4%)
Despite the high sequence diversity of HCV, successful amplification of a near full-length (~9 kb) HCV genotype 1 genome was completed on 97% (7566/7800) of patient plasma samples with HCV RNA concentrations greater than or equal to 1000 IU/mL. This method facilitated the examination of sequence changes in the viral population during DAA treatment and aided the development of telaprevir . The ability to amplify and sequence the near full length HCV genome and to look at mutations across multiple DAA targets will be important in the development of multiple DAA-based regimens.
Studies included in these analyses were conducted in full compliance with the guidelines of Good Clinical Practice and of the World Medical Assembly Declaration of Helsinki. Prior to study initiation, protocols and informed consent forms were reviewed and approved by institutional review boards at each study site. All patients provided written informed consent before participating in any study-related activity. ClinicalTrials Identifiers for studies included in the manuscript are as follows: NCT00336479, NCT00372385, NCT00420784, NCT00535847, NCT00627926, NCT00983853, NCT00758043, NCT00916474, NCT00528528, NCT00703118, NCT00911963, NCT01080222.
This method is based on the amplification protocol originally described in Kwong et al. .
HCV RNA extraction
Viral RNA was extracted from plasma in a 96-well format according to the manufacturer’s instruction with several modifications. Plasma (220 ~ 660 μL, the volume was adjusted based on viral load, i.e., 220 μL for HCV RNA greater than or equal to 50,000 IU/mL, while 660 μL for HCV RNA less than 50,000 IU/mL) was mixed with QIAGEN Protease K (40 ~ 120 μL) and QIAamp AL buffer (240 ~ 720 μL) supplemented with 20 ~ 60 μg of carrier RNA (tRNA) per column. The mixture was incubated for 15 min at 56°C, and absolute ethanol (293 ~ 875 μL) was added to each well and mixed well. The mixture was loaded into QIAamp column and passed through the column by suction with a peristaltic micropump (IPS-16; Ismatec, Zurich, Switzerland). Subsequently, the column was washed with AW1 buffer (1000 μL) and AW2 buffer (1000 μL) (provided in QIAGEN kit), respectively, by applying vacuum. A second wash with AW2 buffer (1000 μL) was spun at 6000 × g for 10 min and an additional spin at 6000 × g for 15 min was applied to remove residual ethanol. 40 μL of RNA storage solution (Ambion) was loaded into each column and incubated at room temperature for 5 min. The RNA was eluted by centrifugation at 6000 × g for 10 min. The elution was repeated once to increase the yield. A total of 80 μL viral RNA was isolated from 220 to 660 μL of plasma.
A complementary DNA (cDNA) fragment was synthesized from HCV RNA that was diluted (1:2 or 1:4, the dilution factor was based on viral load, i.e., 1:2 and 1:4 for HCV RNA less than 50,000 IU/mL, while 1:4 for HCV RNA greater than or equal to 50,000 IU/mL) into a total RT reaction volume of 20 μL composed of 2.5 μM of an oligo-dA20 primer (Table 1), 400 units of Superscript™ III Reverse Transcriptase (Life Technologies Corporation, Carlsbad, CA), 40 units of RNAseOUT (Life Technologies Corporation, Carlsbad, CA), PC2 reaction buffer (50 mM Tris–HCl pH 9.1, 16 mM ammonium sulfate, 3.5 mM magnesium chloride, and 150 μg/mL BSA) (AB Peptides, St. Louis, MO), 500 μM dNTPs (Clontech, Mountain View, CA), and 5 mM DTT (Life Technologies Corporation, Carlsbad, CA). The RNA was denatured at 65°C for 5 min, followed by ramping extension temperatures at 25°C for 10 min, 42°C for 60 min, 50°C for 30 min, 55°C for 30 min, and 70°C for 15 min.
Amplification of HCV via nested PCR
The cDNA from the RT reaction was diluted 1:1 into a 40 μL PCR1 reaction mixture, containing PC2 reaction buffer (AB Peptides, St. Louis, MO), 200 μM dNTPs (Clontech, Mountain View, CA), 1.5 M betaine (Sigma Aldrich, St. Louis, MO), 2.56 units Klentaq DNA polymerase (AB Peptides, St. Louis, MO), 1.28 units pfu DNA polymerase (Stratagene, La Jolla, CA), and 400 μM each subtype-specific primer (Table 1). The PCR1 reaction was incubated at 94°C for 2 min, followed by 30 cycles at 94°C for 15 sec, 68°C −0.5°C/cycle (“touchdown” PCR) for 20 sec, followed by an incubation at 68°C for 12 min. The completed PCR1 reaction was diluted 1:10 into a second PCR mixture (PCR2, 50 μL). The reaction composition and PCR cycling parameters were identical to the PCR1 reaction, with the following changes; 3.2 units Klentaq DNA polymerase, 1.6 units of Pfu DNA polymerase, and nested subtype-specific primers were utilized (Table 1).
Preparation of sample for sequencing
The PCR2 product (8991 bp) was analyzed and verified by agarose gel electrophoresis, and purified using the QIAquick 96 PCR Purification kit (QIAGEN, Valencia, CA). The purified PCR2 product was subsequently quantified using a NanoDrop 8000 Spectrophotometer (Thermo Scientific, Hudson, NH) prior to sequence analysis.
Benjamin Shames and Ann Kwong were employees of Vertex Pharmaceuticals Incorporated at time of study.
Lower limit of detection
Hepatitis C virus
International units per milliliter
Reverse-transcriptase polymerase chain reaction
We thank the study coordinators, nurses, physicians, and patients involved in the clinical trials and Ann Tigges for assistance in generating Shannon Entropy data. We thank Elizabeth Dorn, PhD, Kristin Stephan, PhD, and Erika D. Reynoso, PhD, current employees of Vertex Pharmaceuticals Incorporated, who may own stock or stock options in that company, for editorial support and coordination.
- Butt AA: Hepatitis C virus infection: the new global epidemic. Expert Rev Anti Infect Ther 2005, 3: 241-249. 10.1586/1478718.104.22.168PubMedView Article
- Kim WR: The burden of hepatitis C in the United States. Hepatology 2002, 36: S30-S34.PubMedView Article
- Lavanchy D: The global burden of hepatitis C. Liver Int 2009,29(Suppl 1):74-81.PubMedView Article
- Jacobson IM, McHutchison JG, Dusheiko G, Di Bisceglie AM, Reddy KR, Bzowej NH, Marcellin P, Muir AJ, Ferenci P, Flisiak R, George J, Rizzetto M, Shouval D, Sola R, Terg RA, Yoshida EM, Adda N, Bengtsson L, Sankoh AJ, Kieffer TL, George S, Kauffman RS, Zeuzem S: Telaprevir for previously untreated chronic hepatitis C virus infection. N Engl J Med 2011, 364: 2405-2416. 10.1056/NEJMoa1012912PubMedView Article
- Zeuzem S, Andreone P, Pol S, Lawitz E, Diago M, Roberts S, Focaccia R, Younossi Z, Foster GR, Horban A, Ferenci P, Nevens F, Mullhaupt B, Pockros P, Terg R, Shouval D, van Hoek B, Weiland O, Van Heeswijk R, De Meyer S, Luo D, Boogaerts G, Polo R, Picchio G, Beumont M: Telaprevir for retreatment of HCV infection. N Engl J Med 2011, 364: 2417-2428. 10.1056/NEJMoa1013086PubMedView Article
- Ray SC, Thomas DL: Hepatitis C. In Mandell, Douglas, and Bennett's Principles and Practice of Infectious Diseases. 7th edition. Edited by: Mandell G, Bennett JE, Dolin R. Elsevier, Amsterdam; 2009:2157-2187.
- Simmonds P: Genetic diversity and evolution of hepatitis C virus–15 years on. J Gen Virol 2004, 85: 3173-3188. 10.1099/vir.0.80401-0PubMedView Article
- Adiwijaya BS, Hare B, Caron PR, Randle JC, Neumann AU, Reesink HW, Zeuzem S, Herrmann E: Rapid decrease of wild-type hepatitis C virus on telaprevir treatment. Antivir Ther 2009, 14: 591-595.PubMed
- Kwong AD, Najera I, Bechtel J, Bowden S, Fitzgibbon J, Harrington P, Kempf D, Kieffer TL, Koletzki D, Kukolj G, Lim S, Pilot-Matias T, Lin K, Mani N, Mo H, O'Rear J, Otto M, Parkin N, Pawlotsky JM, Petropoulos C, Picchio G, Ralston R, Reeves JD, Schooley RT, Seiwert S, Standring D, Stuyver L, Sullivan J, Miller V: Sequence and phenotypic analysis for resistance monitoring in hepatitis C virus drug development: recommendations from the HCV DRAG. Gastroenterology 2011, 140: 755-760. 10.1053/j.gastro.2011.01.029PubMedView Article
- Welzel TM, Zeuzem S: Mixing and matching drugs: what makes sense? Clin Liver Dis 2011, 15: 657-664. 10.1016/j.cld.2011.05.012PubMedView Article
- Aizaki H, Aoki Y, Harada T, Ishii K, Suzuki T, Nagamori S, Toda G, Matsuura Y, Miyamura T: Full-length complementary DNA of hepatitis C virus genome from an infectious blood sample. Hepatology 1998, 27: 621-627. 10.1002/hep.510270242PubMedView Article
- Liu Z, Netski DM, Mao Q, Laeyendecker O, Ticehurst JR, Wang XH, Thomas DL, Ray SC: Accurate representation of the hepatitis C virus quasispecies in 5.2-kilobase amplicons. J Clin Microbiol 2004, 42: 4223-4229. 10.1128/JCM.42.9.4223-4229.2004PubMedPubMed CentralView Article
- Lu L, Nakano T, Smallwood GA, Heffron TG, Robertson BH, Hagedorn CH: A refined long RT-PCR technique to amplify complete viral RNA genome sequences from clinical samples: application to a novel hepatitis C virus variant of genotype 6. J Virol Methods 2005, 126: 139-148. 10.1016/j.jviromet.2005.01.031PubMedView Article
- Tellier R, Bukh J, Emerson SU, Miller RH, Purcell RH: Long PCR and its application to hepatitis viruses: amplification of hepatitis A, hepatitis B, and hepatitis C virus genomes. J Clin Microbiol 1996, 34: 3085-3091.PubMedPubMed Central
- Fan X, Xu Y, Di Bisceglie AM: Efficient amplification and cloning of near full-length hepatitis C virus genome from clinical samples. Biochem Biophys Res Comm 2006, 346: 1163-1172. 10.1016/j.bbrc.2006.06.039PubMedView Article
- Xu Z, Fan X, Xu Y, Di Bisceglie AM: Comparative analysis of nearly full-length hepatitis C virus quasispecies from patients experiencing viral breakthrough during antiviral therapy: clustered mutations in three functional genes, E2, NS2, and NS5a. J Virol 2008, 82: 9417-9424. 10.1128/JVI.00896-08PubMedPubMed CentralView Article
- Zhou D, Fan X, Tan D, Xu Y, Tavis JE, Di Bisceglie AM: Separation of near full-length hepatitis C virus quasispecies variants from a complex population. J Virol Methods 2007, 141: 220-224. 10.1016/j.jviromet.2006.12.002PubMedView Article
- INCIVEK™: US perscribing information. Vertex Pharmaceuticals Incorporated, Cambridge, MA; 2012.
- INCIVO®: EU Summary of Product Characteristics. Janssen Pharmaceuticals, Beerse, Belgium; 2012.
- Hézode C, Forestier N, Dusheiko G, Ferenci P, Pol S, Goeser T, Bronowicki JP, Bourliere M, Gharakhanian S, Bengtsson L, McNair L, George S, Kieffer T, Kwong A, Kauffman RS, Alam J, Pawlotsky JM, Zeuzem S: Telaprevir and peginterferon with or without ribavirin for chronic HCV infection. N Engl J Med 2009, 360: 1839-1850. 10.1056/NEJMoa0807650PubMedView Article
- McHutchison JG, Everson GT, Gordon SC, Jacobson IM, Sulkowski M, Kauffman R, McNair L, Alam J, Muir AJ: Telaprevir with peginterferon and ribavirin for chronic HCV genotype 1 infection. N Engl J Med 2009, 360: 1827-1838. Erratum: N Engl J Med. 2009, 361:1516 10.1056/NEJMoa0806104PubMedView Article
- McHutchison JG, Manns MP, Muir AJ, Terrault NA, Jacobson IM, Afdhal NH, Heathcote EJ, Zeuzem S, Reesink HW, Garg J, Bsharat M, George S, Kauffman RS, Adda N, Di Bisceglie AM: Telaprevir for previously treated chronic HCV infection. N Engl J Med 2010, 362: 1292-1303. Erratum: N Engl J Med. 2010, 362:1647. 10.1056/NEJMoa0908014PubMedView Article
- Sherman KE, Flamm SL, Afdhal NH, Nelson DR, Sulkowski MS, Everson GT, Fried MW, Adler M, Reesink HW, Martin M, Sankoh AJ, Adda N, Kauffman RS, George S, Wright CI, Poordad F: Response-guided telaprevir combination treatment for hepatitis C virus infection. N Engl J Med 2011, 365: 1014-1024. Erratum: N Engl J Med. 2011, 365:1551 10.1056/NEJMoa1014463PubMedView Article
- Kuiken C, Yusim K, Boykin L, Richardson R: The Los Alamos HCV sequence database. Bioinformatics 2005,21(3):379-384. 10.1093/bioinformatics/bth485PubMedView Article
- Henke W, Herdel K, Jung K, Schnorr D, Loening SA: Betaine improves the PCR amplification of GC-rich DNA sequences. Nucleic Acids Res 1997, 25: 3957-3958. 10.1093/nar/25.19.3957PubMedPubMed CentralView Article
- Rees WA, Yager TD, Korte J, von Hippel PH: Betaine can eliminate the base pair composition dependence of DNA melting. Biochemistry 1993, 32: 137-144. 10.1021/bi00052a019PubMedView Article
- Shannon Claude E: A Mathematical Theory of Communication. Bell Syst Tech J 1948,27(3):379-423.View Article
- Kieffer TL, De Meyer S, Bartels DJ, Sullivan JC, Zhang EZ, Tigges A, Dierynck I, Spanks J, Dorrian J, Jiang M, Adiwijaya B, Ghys A, Beumont M, Kauffman RS, Adda N, Jacobson IM, Sherman KE, Zeuzem S, Kwong AD, Picchio G: Hepatitis C viral evolution in genotype 1 treatment-naive and treatment-experienced patients receiving telaprevir-based therapy in clinical trials. PLoS One 2012, 7: e34372. 10.1371/journal.pone.0034372PubMedPubMed CentralView Article
- Kwong AD, Frantz JD, Bartels DJ, Lin C, Shames B, Seepersaud S, Lippke JA, Kieffer TL, Zhou Y, Zhang EZ, Sullivan JC, inventors; Vertex Pharmaceuticals Incorporated, assignee: Methods for amplifying hepatitis C virus nucleic acids. 2010. US patent WO/20120/090857 2010
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.