- Open Access
Ultrasensitive quantitation of human papillomavirus type 16 E6 oncogene sequences by nested real time PCR
Infectious Agents and Cancervolume 5, Article number: 9 (2010)
We have developed an ultrasensitive method based on conventional PCR preamplification followed by nested amplification through real time PCR (qPCR) in the presence of the DNA intercalating agent EvaGreen.
Amplification mixtures calibrated with a known number of pHV101 copies carrying a 645 base pair (bp)-long insert of the human papillomavirus type 16 (HPV16) E6 oncogene were used to generate the E6-1 amplicon of 645 bp by conventional PCR and then the E6-2 amplicon of 237 bp by nested qPCR. Direct and nested qPCR mixtures for E6-2 amplification corresponding to 2.5 × 102-2.5 × 106 initial pHV101 copies had threshold cycle (Ct) values in the ranges of 18.7-29.0 and 10.0-25.0, respectively. The Ct of qPCR mixtures prepared with 1/50 volumes of preamplified mixtures containing 50 ng of DNA of the SiHa cell line (derived from an invasive cervical cancer with one HPV16 genome per cell) was 19.9. Thermal fluorescence extinction profiles of E6-2 amplicons generated from pHV101 and SiHa DNA were identical, with a peak at 85.5°C.
Our method based on conventional preamplification for 15 cycles increased 10,750 times the sensitivity of nested qPCR for the quantitation of the E6 viral oncogene and confirmed that the SiHa cell line contains one E6-HPV16 copy per cell.
Invasive cervical cancer develops in women with persistent infection by high-risk human papillomavirus (HR-HPV), among which HPV16 and HPV18 are the most frequent types , HPV16 being the most prevalent and associated with nearly half of the invasive cervical cancer cases in the world [2, 3].
Molecular methods to quantify the number of HPV16 genome copies serve to determine the viral load and the progression of HPV16 neoplastic lesions. Hybrid capture  or "nested" PCR  are used for molecular diagnosis of HPV infection, whereas qPCR is used to assess viral load and integration to the cellular genome [6, 7].
Conventional nested PCR is more sensitive than conventional "direct" PCR to detect HPV16 E6 oncogene sequences [8, 9] whose quantitation is reproducible, specific and more sensitive in qPCR mixtures containing the DNA intercalating agent EvaGreen .
In this work we developed an ultrasensitive method to quantify E6-HPV16 oncogene sequences in two steps: 1) preamplification of the E6-1 (645 bp) sequence through 15 cycles of direct conventional PCR, and 2) amplification of the E6-2 (237 bp) sequence through nested qPCR in the presence of EvaGreen. Preamplification significantly increased the sensitivity and the method confirmed that the SiHa cell line contains only one copy of the HPV16 genome per cell.
E6-2 amplification by direct qPCR
E6-2 amplification from pHV101 in direct qPCR mixtures containing the pU1M/pU2R primer pair and 2.5 × 102 to 2.5 × 106 pHV101 copies in the presence of EvaGreen yielded a consistent family of type curves with Ct values between 18.70 and 29.00, ΔCt values in the 3.1-3.8 range, and lack of fluorescence increase in mixtures devoid of DNA. These experiments indicated that at least 2,500 copies of the E6 oncogene can be quantified by direct qPCR (Fig. 1A).
Preamplification and nested qPCR
To increase qPCR sensitivity E6-1 was first "preamplified" in conventional PCR mixtures and E6-2 was then "amplified" in nested qPCR mixtures containing EvaGreen, the LCRS/E7AS primer pair, and 1/50 volume of the corresponding preamplified mixtures.
The families of type curves required to calibrate the number of E6 copies and to quantify them in problem samples could be generated from nested qPCR mixtures derived from two types of preamplification mixtures: a) "complete" mixtures, prepared with serial logarithmic dilutions of the pHV101 template molecules before preamplification, and b) "minimum" mixtures, prepared with serial logarithmic dilutions from a single preamplified mixture containing the maximum number of template molecules used in complete mixtures.
Since labor, time and the amount of reagents are significantly reduced using minimum mixtures, we compared the results from nested qPCR mixtures prepared with 1/50 volume samples from preamplified mixtures which contained pHV101 molecules in the range of 2.5 × 102-to 2.5 × 106 copies with those prepared using serial logarithmic dilutions from a single mixture containing 2.5 × 106 pHV101 molecules that had been preamplified. We also added 50 ng of "carrier" normal human blood DNA to both complete and minimum positive preamplification mixtures (same amount as SiHa DNA in problem mixtures) and found that all of them generate equivalent nested qPCR amplification results.
Typical Ct values of nested qPCR: 10.80 for 2.5 × 106, 14.35 for 2.5 × 105, 17.71 for 2.5 × 104, 21.62 for 2.5 × 103 and 25.02 for 2.5 × 102 initial pHV101 copies in the preamplification mixtures (Table 1) were consistent (ΔCt range: 3.3-3.9) and reproducible (R2 = 1.000), whereas the negative control mixture devoid of template DNA had a Ct = 27 (Fig. 1B).
The average decrease in Ct values of nested qPCR resulting from preamplification for 15 cycles was 7.78 (Table 1), corresponding to a 215-fold sensitivity increase. Multiplying this factor by the 50-fold dilution factor of preamplified samples in nested reactions implies that under the conditions used the overall sensitivity of nested qPCR increased 10,750 times.
E6 copies per SiHa cell
The sensitivity and reliability of the nested qPCR method was tested by determining the HPV16 genome load through quantitation of E6 sequences present in the SiHa line, known to contain only one copy of the HPV16 genome per cell .
The following PCR mixtures containing the LCRS/E7AS primer pair and the DNA templates indicated were incubated for 15 cycles under conditions to generate ("preamplify") the E6-1 amplicon: a) positive controls with 2.5 × 106-2.5 × 102 copies of pHV101 and 50 ng of "carrier" normal human blood DNA; b) negative controls with 50 ng of "carrier" DNA; c) problem samples with 50 ng of SiHa DNA; and d) blanks, without DNA. Nested qPCR mixtures containing EvaGreen, the pU1M/pU2R primer pair and 1/50 volume from each preamplified mixture were incubated for 35 cycles.
In nested qPCR mixtures for E6-2 containing 1/50 volumes of serial logarithmic dilutions from preamplified mixtures with 2.5 × 106 to 2.5 × 102 pHV101 molecules the Ct values ranged from 10.8 to 25.0 cycles (R2 = 1.000; slope = -3.57), and the ΔCt values were consistent (range: 3.3-3.9). The Ct values of nested qPCR mixtures prepared without preamplified samples were 27.0 and 28.0, whereas the basal fluorescence decreased slightly in preamplification mixtures containing SiHa DNA with the pU1M primer only (Fig. 2A).
The family of type curves obtained from nested qPCR mixtures for E6-2 amplification was consistent with the initial pHV101 copy number. On the other hand, the Ct value (19.9) and E6 copy number obtained from a typical preamplification mixture containing 50 ng of SiHa cell DNA (equivalent to 6.05 × 103 HPV16 genome copies) had the magnitudes expected (i.e., one E6 copy per 7.1 pg of DNA, equivalent to the genomic mass of a human diploid cell) (Fig. 2A).
The thermal denaturation ("fluorescence extinction") profiles of nested qPCR products generated from pHV101 and SiHa DNA were identical, with a major peak around 85.5°C whose magnitude was proportional to the initial number of pHV101 copies, and clearly different from the secondary peak at 78.5°C which appeared in the negative controls containing just the forward primer of the second pair (Fig. 2B).
Our group has already confirmed that detection of HPV infection is more sensitive by nested conventional PCR through the successive synthesis of the E6-1 and E6-2 amplicons used in this work , and that quantitation of E6 oncogene sequences by qPCR in the presence of EvaGreen is reproducible and specific .
The PCR mixtures required to calibrate and quantify the number of E6 copies present in problem samples were prepared using serial logarithmic dilutions of pHV101, a pGEM construct containing a 645 bp-insert whose sequence corresponds to the E6-HPV16 ORF .
To increase amplification sensitivity, E6-1 was "preamplified" for 15 conventional PCR cycles and then E6-2 was "amplified" by 35 "nested" qPCR cycles in the presence of normal human DNA using, as source of template, 1/50 volume of either preamplified mixtures prepared with serial logarithmic dilutions of pHV101 or with serial logarithmic dilutions of a single preamplified mixture prepared with the highest number of pHV101 copies. Direct PCR could quantify as little as 2,500 E6-HPV16 molecules, whereas nested qPCR was around 11,000 times more sensitive, a value close to that expected for 15 preamplification cycles which would increase the initial number of E6 template copies by a factor of 213 = 8,192.
Non-specific fluorescence appears to depend on the interaction of the two E6-2 primers because it increased after cycle 27 in qPCR mixtures prepared with samples from preamplified mixtures containing normal human DNA or HPV16 DNA, but not from those containing SiHa DNA and the E6-2 forward primer only.
Use of the SiHa cell line, originally isolated from an invasive cervical cancer caused by HPV16 and known to contain only one viral genome per cell , confirmed the sensitivity and reliability of our nested qPCR method, since preamplified mixtures containing 50 ng of SiHa DNA generated Ct values corresponding to one copy of the E6-HPV16 oncogene per human diploid genome.
The nested qPCR method may also be used in the future to quantify viral load as well as viral and cellular transcripts with high sensitivity and reliability.
Our method, based on conventional PCR preamplification for 15 cycles increased 10,750 times the sensitivity of nested qPCR for quantitation of the E6 oncogene, and confirmed that the SiHa cell line contains only one E6-HPV16 copy per cell.
E-645, the pHV101 insert sequence spanning the E6-HPV16 open reading frame (ORF), was used as the initial template to generate the E6-1 (645 bp) amplicon through conventional direct PCR and then the E6-2 (237 bp) amplicon through conventional nested PCR (Fig. 3) .
E6-1 and E6-2 amplification was maximized in 50 μL conventional PCR mixtures containing 2.5 × 106 pHV101 copies under the following conditions: 3 mM MgCl2, 0.05 mM of each deoxynucleoside triphosphate (dNTP), 0.15 μM of each of the forward and reverse primers (Table 2), 30 U/ml Taq polymerase (Invitrogen), and annealing temperature of 57°C.
PCR products were analyzed by electrophoresis in sodium borate-1.5% agarose gels  for 2 h at 90 V, followed by ethidium bromide staining; they were visualized by ultraviolet light transillumination and their digital images recorded in a ChemiDoc EQ photodocumenter (BioRad, Hercules, CA) using the Quantity One (BioRad) software.
Preamplification and calibration of nested qPCR
Calculation of the number of E6-HPV16 copies is based on the size of pHV101 (3,645 bp) and the average molecular weight of a deoxynucleotide pair (650 Da). We estimated that 1 ng of pHV101 contains 2.5 × 108 molecules with the formula of Staroscik : number of copies = ((amount in ng) * 6.022 × 1023))/((length in bp) * 109 * 650).
The E6-1 positive control mixtures required to generate the family type curves contained the LCRS/E7AS primer pair, 50 ng of "carrier" normal human DNA and serial logarithmic dilutions of pHV101; they were preamplified for 15 cycles (denaturation at 94°C for 15 sec, annealing at 57°C for 1 min and extension at 72°C for 1 min). Problem preamplification mixtures contained 50 ng of SiHa cell DNA. Negative controls were a blank (i.e., without DNA) and another one with 50 ng of carrier normal human blood DNA.
Positive control preamplification mixtures were prepared either a) as five "complete" PCR mixtures containing serial logarithmic dilutions of pHV101 to attain 2.5 × 106-2.5 × 102 pHV101 copies/mixture (Fig. 4A), or b) as a single "minimum" PCR mixture containing 2.5 × 106 pHV101 copies (Fig. 4B); 1/50 volumes from both types of mixtures were added to prepare nested qPCR amplification mixtures.
qPCR mixtures optimized to amplify E6-2 with the pU1M/pU2R primer pair in the presence of EvaGreen (Biotium, Hayward, CA) were incubated with the same thermocycler program. To control nested amplification and to determine the number of E6 copies in problem samples, additional mixtures containing 1/50 volume of preamplified positive, negative and problem mixtures were prepared; in these experiments an additional negative control mixture was used which contained preamplified SiHa DNA and the forward (pU1M) but not the reverse (pU1R) primer (Fig. 4).
Munoz N, Bosch FX, de Sanjose S, Herrero R, Castellsague X, Shah KV, Snijders PJ, Meijer CJ: Epidemiologic classification of human papillomavirus types associated with cervical cancer. N Engl J Med. 2003, 348 (6): 518-527. 10.1056/NEJMoa021641.
Bosch FX: Epidemiology of human papillomavirus infections: new options for cervical cancer prevention. Salud Publica Mex. 2003, 45 (Suppl 3): S326-339.
Walboomers JM, Jacobs MV, Manos MM, Bosch FX, Kummer JA, Shah KV, Snijders PJ, Peto J, Meijer CJ, Munoz N: Human papillomavirus is a necessary cause of invasive cervical cancer worldwide. J Pathol. 1999, 189 (1): 12-19. 10.1002/(SICI)1096-9896(199909)189:1<12::AID-PATH431>3.0.CO;2-F.
Gravitt PE, Burk RD, Lorincz A, Herrero R, Hildesheim A, Sherman ME, Bratti MC, Rodriguez AC, Helzlsouer KJ, Schiffman M: A comparison between real-time polymerase chain reaction and hybrid capture 2 for human papillomavirus DNA quantitation. Cancer Epidemiol Biomarkers Prev. 2003, 12 (6): 477-484.
Sotlar K, Diemer D, Dethleffs A, Hack Y, Stubner A, Vollmer N, Menton S, Menton M, Dietz K, Wallwiener D, et al: Detection and typing of human papillomavirus by e6 nested multiplex PCR. J Clin Microbiol. 2004, 42 (7): 3176-3184. 10.1128/JCM.42.7.3176-3184.2004.
Hubbard RA: Human papillomavirus testing methods. Arch Pathol Lab Med. 2003, 127 (8): 940-945.
Mayrand MH, Duarte-Franco E, Rodrigues I, Walter SD, Hanley J, Ferenczy A, Ratnam S, Coutlee F, Franco EL: Human papillomavirus DNA versus Papanicolaou screening tests for cervical cancer. N Engl J Med. 2007, 357 (16): 1579-1588. 10.1056/NEJMoa071430.
Noda T, Sasagawa T, Dong Y, Fuse H, Namiki M, Inoue M: Detection of human papillomavirus (HPV) DNA in archival specimens of benign prostatic hyperplasia and prostatic cancer using a highly sensitive nested PCR method. Urol Res. 1998, 26 (3): 165-169. 10.1007/s002400050041.
Lopez-Revilla R, Martinez-Contreras LA, Sanchez-Garza M: Prevalence of high-risk human papillomavirus types in Mexican women with cervical intraepithelial neoplasia and invasive carcinoma. Infect Agent Cancer. 2008, 3: 3-10.1186/1750-9378-3-3.
Hernandez-Arteaga S, Lopez-Revilla R: Quantitation of human papillomavirus type 16 E6 oncogene sequences by real-time or quantitative PCR with EvaGreen. Anal Biochem. 2008, 380 (1): 131-133. 10.1016/j.ab.2008.05.012.
Friedl F, Kimura I, Osato T, Ito Y: Studies on a new human cell line (SiHa) derived from carcinoma of uterus. I. Its establishment and morphology. Proc Soc Exp Biol Med. 1970, 135 (2): 543-545.
Brody JR, Kern SE: Sodium boric acid: a Tris-free, cooler conductive medium for DNA electrophoresis. Biotechniques. 2004, 36 (2): 214-216.
Calculator for determining the number of copies of a template. [http://www.uri.edu/research/gsc/resources/cndna.html]
This work was partially supported by the Fondo Mixto CONACYT-SLP (grant FMSLP-C01-2006-62319) and a scholarship to SHA (No. 176112) from the Consejo Nacional de Ciencia y Tecnología (Mexico).
The authors declare that they have no competing interests.
RLR conceived the study and obtained the funds to carry it out. SHA performed the bioinformatics analyses and the molecular studies. Both authors drafted the manuscript.