|
|
||||||||
Molecular Diagnostics and Genetics |
1 Department of Biotechnology and Biosciences, University of Milano Bicocca, Milano, Italy.
2 DiaSorin SpA, Saluggia (VC), Italy.
3 Dana Farber-Brigham and Womens Cancer Center, Harvard Medical School, Boston, MA.
aAddress correspondence to this author at: Diasorin SpA, Piazza della Scienza 4, 20126 Milano, Italy. Fax 0039-02-6448-3565; e-mail daniel.adlerstein{at}diasorin.it.
| Abstract |
|---|
|
|
|---|
Materials and Methods: We analyzed KRAS codon 12 sequence variants in 106 lung cancer surgical specimens. We conducted a PNA-PCR reaction that suppresses wild-type KRAS amplification and genotyped the product with a set of OCEAN reactions carried out in fluorescence microplate format. The isothermal OCEAN assay enabled a 3-way DNA junction to form between the specific target nucleic acid, a fluorescently labeled "amplifier", and an "anchor". The amplifier-anchor contact contains the recognition site for a restriction enzyme. Digestion produces a cleaved amplifier and generation of a fluorescent signal. The cleaved amplifier dissociates from the 3-way DNA junction, allowing a new amplifier to bind and propagate the reaction.
Results: The system detected and genotyped KRAS sequence variants down to
0.3% variant-to-wild-type alleles. PNA-PCR/OCEAN had a concordance rate with PNA-PCR/sequencing of 93% to 98%, depending on the exact implementation. Concordance rate with restriction endonuclease-mediated selective-PCR/sequencing was 89%.
Conclusion: OCEAN is a practical and low-cost novel technology for sequence-specific signal generation. Reliable analysis of KRAS sequence alterations in human specimens circumvents the requirement for sequencing. Application is expected in genotyping KRAS codon 12 sequence variants in surgical specimens or in bodily fluids, as well as single-base variations and sequence alterations in other genes.
| Introduction |
|---|
|
|
|---|
In most current methods, nucleic acid amplification via PCR (17) from genomic DNA is the first step in the identification of single-nucleotide variants and sequence alterations. Detection is performed either during PCR [e.g., TaqMan real-time PCR (1)] or after PCR [e.g., primer-extension, dHPLC, matrix-assisted laser desorption/ionization time-of-flight (MALDI-TOF), or beads-based approaches (18)(19)(20)(21)]. Detection of single-nucleotide differences during PCR is more convenient but usually requires an additional screening step to provide the specific genotype, and some postPCR assays, such as dHPLC, may require a third step. We describe a novel signal generation and amplification technology that enables simultaneous detection and genotyping of single-nucleotide variants and sequence alterations in a single postPCR step and combines most of the requirements for convenient single-nucleotide genotyping.
| Materials and Methods |
|---|
|
|
|---|
|
Double-labeling of the amplifier with a fluorescent probe and a quencher (Fig. 1
) generates a fluorescent signal by eliminating fluorescent quenching each time the amplifier is digested (22). Repeated digestion of the amplifier leads to a continuously increasing fluorescent signal that is monitored homogeneously. The stability of the 3-way DNA junction at the selected temperature is dependent on the correct hybridization of the 3 sequences. Therefore, if the amplifier and the interrogated DNA sequence differ by a single nucleotide in the match-region, the junction is destabilized and the fluorescent signal is diminished. Thus, with appropriate design of the amplifier sequence, any chosen nucleotide change on the target DNA sequence can be genotyped.
To validate this OneCutEventAmplificatioN (OCEAN) signal-generation method (23) and adapt it to the detection of clinically relevant sequence variants, we designed assays appropriate for genotyping codon-12 sequence alterations in the KRAS gene. We tested the application of OCEAN in model-oligonucleotide systems, in cell lines with KRAS sequence variants, and in lung cancer surgical specimens. We compared our results with those obtained with peptide-nucleic acid (PNA)-PCR (24), restriction endonuclease-mediated selective (REMS)-PCR (25), and sequencing. We demonstrate that OCEAN combined with asymmetric PNA-PCR variant enrichment (see Fig. S1 in the online Data Supplement), can detect and simultaneously genotype KRAS codon 12 sequence alterations in human tissue specimens.
source and extraction of synthetic and genomic dna
Synthetic DNA.
To develop and optimize the OCEAN assays, the synthetic oligonucleotides GA61, GA62, GA93, GA94, and GA95 (see Table S1 in the online Data Supplement), representing DNA containing different genotypes of KRAS codon 12 (GGT-wild-type, GAT, GCT, GTT, and TGT), were synthesized by Sigma-Proligo.
Plasmids.
For the wild-type and the 4 most common KRAS codon 12 sequence alterations (GGT, GTT, GAT, GCT, and TGT), synthetic 166-bp fragments from 95 bp upstream to 68 bp downstream of this codon were cloned into pCR-Blunt plasmids (Invitrogen) previously digested with Stu I restriction endonuclease (New England Biolabs). The resulting recombinant plasmids were cultured in Escherichia coli strain DH5a and extracted with a QIAprep Spin Miniprep Kit (Qiagen).
Control cell lines.
To verify the specificity and selectivity of PNA-PCR/OCEAN, genomic DNA from cultured cell lines SW480 (homozygous GTT), PL45 (heterozygous GAT/GGT), and CALU 1 (heterozygous TGT/GGT), obtained from the American Type Culture Collection, (Manassas, VA) was extracted with the NucleoSpinTM Tissue method (Macherey-Nagel).
Clinical samples.
We used surgical specimens from 106 patients who underwent surgery for lung adenocarcinoma and nonsmall cell lung cancer (NSCLC). The specimens were obtained after we received internal review board approval from the Massachussets General Hospital Tumor Bank and from the Cooperative Human Tissue Network. Specimens were removed immediately after surgery and snap-frozen in liquid nitrogen. The specimens were accompanied by a pathology report that characterized them as normal lung tissue (NL) or tumor lung tissue (TL). Genomic DNA from lung samples was extracted at the Dana Farber Cancer Institute with a QIAampTM DNA Mini Kit, (Qiagen Inc.). Typically, 30 mg lung tissue yielded 2030 µg of purified genomic DNA.
ocean assays on synthetic oligonucleotides
To validate the OCEAN assay for KRAS sequence variant detection, we used the synthetic oligonucleotides GA61, GA62, GA93, GA94, and GA95, the anchor, and amplifier-oligonucleotides listed in Table S1 in the online Data Supplement. We designed all OCEAN probes with VisualOmp4 software (DNA Software, Inc.) with a custom made add-on for calculating melting temperatures of 3-way DNA junctions, based on the methods of SantaLucia et al. (26)(27). Each reaction contained 1 nmol/L anchor probe, 200 nmol/L amplifier internally labeled with fluorescein, and with a quencher (BHQ) at the 3' end, 200 nmol/L unlabeled amplifier complementary to the wild-type, and 0.5 units/µL BsoBI enzyme that digests the sequence 5'-CTCGGG-3', in 20 µL total volume. After 1 h incubation at 65 °C, the fluorescence signals from the amplifier cleavage fragments were detected in 15 µL samples and quantified on a laser scanner (Typhoon 9200, Amersham Biosciences).
enrichment of variants using pna-pcr on plasmids and on cell lines with kras sequence variants
PNAs that form a stable duplex with the wild-type sequence of KRAS codon 12 suppress wild-type amplification and prevent annealing and extension of the reverse primer during PCR (24). When a KRAS variant is present, the resulting mismatch reduces the stability of the PNA-DNA duplex; thus, amplification occurs as the primer competes with the mismatched PNA-DNA duplex. We designed and synthesized (Eurogentec) a PNA probe (see Table S1 in the online Data Supplement) to be used in combination with OCEAN-based identification of the genotype, and the preferential amplification of KRAS variant sequences was tested in recombinant plasmids and human cell lines containing KRAS sequence alterations. For plasmids, the PNA-PCR reaction was conducted on 500 pM of plasmid containing 1:1 KRAS codon 12 wild-type (GGT) or variant (GAT). The experiment was performed in the absence and the presence of 1 µmol/L PNA probe complementary to the wild-type sequence. The reaction was performed in 20 µL as follows: 120 s of initial denaturation at 94 °C, 60 s at 94 °C, 50 s at 70 °C, 60 s at 53 °C, and 60 s at 68 °C, for a total of 25 cycles followed by a 10-min final extension at 68°. For the human cell lines, 68 ng of genomic DNA was used in PNA-PCR reactions containing 1 µmol/L PNA, 1 x BD TitaniumTM Taq polymerase, PCR buffer (BD Biosciences), 200 µmol/L dNTPs, 3% formamide, and 2.5 mmol/L MgCl2. The reaction was carried out at 95 °C for 5 min, followed by 35 cycles of 94 °C for 30 s (denaturation), 54 °C for 30 s (primers and PNA annealing), and 72 °C for 30 s (extension), with a final extension step of 72 °C for 5 min. The products were examined via ethidium-stained gel electrophoresis.
asymmetric pna-pcr-ocean
We used PNA-PCR followed by OCEAN as a combined method for detection of the most common sequence variants of KRAS codon 12 in clinical lung cancer samples (see Fig. S1 in the online Data Supplement). An excess (5:1) of the forward primer was used in the PNA-PCR reactions to accumulate the upper DNA strand (single-stranded product), which is the substrate for the OCEAN amplifier and anchor. After PNA-PCR, we examined the sample via gel electrophoresis. For each possible type of KRAS codon 12 sequence variant, we used the doubly labeled amplifiers GA85QF-GA88QF (see Table S1 in the online Data Supplement) to perform a series of specific OCEAN reactions on 1:200 dilutions of the PCR product. The reactions were performed at 65 °C for 1 h in the presence of 200 nmol/L unlabeled wild-type amplifier, 200 nmol/L labeled amplifier specific for the sequence variant, 1 nmol/L anchor probe, and 0.5 units/µL BsoBI enzyme. A 15-µL aliquot from each OCEAN reaction was loaded onto 384-well microplates (OptiPlate-384 Black, Perkin-Elmer) and read on a fluorescence microplate reader (FluoDia T70, PTI). For each measurement we used the instrument software to calculate the mean of 5 fluorescence readings. The fluorescence values obtained were corrected by subtracting the mean value obtained from identical experiments carried out in parallel on wild-type genomic DNA. After subtraction of the wild-type, we empirically selected a threshold of 1000 fluorescent units to distinguish wild-type from variant samples. Thus, negative signals or signals <1000 a.u. were determined to be wild-type. Experiments were repeated in triplicate and the SDs were calculated.
rems-pcr protocol
We adapted the REMS-PCR protocol to a TaqMan real-time PCR format and used the thermostable BstNI enzyme, which digests wild-type KRAS alleles (25). The REMS-PCR primers used were: 5'-ACT GAA TAT AAA CTT GTG GTA GTT GGA CCT-3' (forward) and 5'-GCA TAT TAA AAC AAG ATT TAC-3' (reverse). We used Oligo 6 software (Molecular Biology Insights) to design a Taqman probe sequence close to but not overlapping the KRAS codon 12 sequence variant region 5'-FAM-CAC TCT TGC CTA CGC CAC CAG-TAMRA-3'. We used REMS-PCR to detect sequence variants in codon 12 of the KRAS gene in the same samples analyzed by PNA-PCR/OCEAN and compared the results. In each PCR reaction, we used 30 ng of genomic DNA as a template. Along with each test of clinical samples, we used 30 ng of human male genomic DNA (wild-type negative control), SW480 genomic DNA (variant control), and no-template as positive controls. We performed PCR for 55 cycles (94 °C for 15 s and 55 °C for 1 min) in a reaction containing PCR buffer (10 mmol/L Tris-HCl, pH 8.3, 50 mmol/L KCl, 0.01% gelatin, 1.5 mmol/L MgCl2), 0.2 mmol/L dNTPs, 0.26 µmol/L of probe (from IDT Inc.), 0.2 µmol/L each of KRAS forward and reverse primers, 0.625 units Taq DNA polymerase (Promega Corporation), and 10 units of BstNI (New England Biolabs). Experiments were repeated in triplicate and SDs were calculated.
| Results |
|---|
|
|
|---|
|
We incubated a KRAS sequence corresponding to a GAT codon 12 sequence variant with either fully matched (GAT) or mismatched (TGT) amplifier and used fluorescence gel imaging to quantify the percentage of input amplifier digested by the enzyme as a function of time. The gradual increase of the fluorescent signal from the cut amplifier is depicted in Fig. 2B
. For the fully matched sequence, the fluorescent signal from the digested amplifier increased with time, consistent with the expectation that after enzymatic digestion the cut amplifier dissociates and is replaced by a new, intact amplifier to propagate the reaction. In contrast, for the mismatched amplifier, the signal amplification was substantially weaker. After a 2 h enzymatic reaction, 81.2% or 1.7%, respectively, of the total input of matched or mismatched amplifier was digested.
pna-pcr detection of variant kras sequences with on plasmids and human cell lines
Before OCEAN-genotyping of the sequence variants, we enriched the variant sequences by PNA to detect small concentrations of somatic KRAS codon 12 sequence alterations in the presence of an excess of wild-type sequences. After amplification from plasmids engineered to contain the wild-type (GGT) KRAS codon 12 sequence, the wild-type sequences were suppressed in the presence of PNA (Fig. 3A
). In contrast, despite the presence of PNA plasmids containing a variant (GAT) KRAS codon, 12 sequences displayed efficient amplification of the KRAS sequence. Next, we amplified KRAS codon 12 sequences with genomic DNA from 4 human cell lines containing wild-type or variant KRAS (K562, wild-type GGT; SW480, homozygous GTT; PL45, heterozygous GAT/GGT; and CALU 1, heterozygous TGT/GGT), in the presence or the absence of PNA. The PNA-PCR reaction suppressed the amplification of the wild-type gene, leading to a selective amplification of KRAS variant templates (Fig. 3B
). The selective amplification of variant KRAS in the presence of PNA was more prominent for the homozygous variant cell line SW480 than for the heterozygous cell lines PL45 and Calu 1. The latter 2 cell lines also contain wild-type sequences not expected to amplify efficiently in the presence of PNA.
|
pna-pcr/ocean on genomic dna from human cells
To combine selective PNA-PCR amplification of KRAS codon 12 sequence variants with OCEAN-based genotyping, we carried out the PNA-PCR reaction with an unbalanced primer ratio (forward:reverse, 5:1) to preferentially generate the upper DNA strand, which acts as the template for the OCEAN reaction. After the asymmetric PNA-PCR reaction on the variant cell lines SW480, PL45, and CALU 1, we carried out OCEAN reactions with oligonucleotide amplifiers designed to screen for the most common KRAS codon 12 sequence alterations encountered in human lung adenocarcinoma and NSCLC, (GTT, GAT, GCT, TGT) for 1 h on 1:200 dilutions of the amplification product. We read the fluorescence on a fluorescence microplate reader and corrected the fluorescence values obtained by subtracting the mean value from identical experiments carried out in parallel on wild-type genomic DNA. The fluorescence signals obtained when oligonucleotide amplifiers matching the correct genotype were used were substantially stronger than for nonmatching amplifiers (Fig. 4A
). The PNA-PCR/OCEAN assay detected the presence and genotype of sequence alterations correctly for the 3 variant human cell lines.
|
selectivity of pna-pcr/ocean
To examine the selectivity of the combined PNA-PCR/OCEAN assay, we serially diluted genomic DNA from the homozygous GTT variant cell line SW480 into wild-type KRAS codon 12 DNA from K562 cells to obtain decreasing ratios of variant-to-wild-type alleles, 100%, 30%, 10%, 3%, 1%, 0.3%, and 0%. We performed triplicate independent experiments with an asymmetric PNA-PCR reaction with a primer ratio of 2:1, 50 ng of total genomic DNA, 1 µmol/L probe PNA complementary to the wild-type sequence, followed by a 5-h OCEAN assay containing the labeled GTT amplifier, with fluorescence detected on a microplate reader. The combined PNA-PCR/OCEAN methodology enabled reproducible differentiation of the variant from the wild-type down to a 0.3% variant-to-wild-type ratio (Fig. 4B
).
rems-pcr screening for kras codon 12 sequence alterations
The REMS-PCR method uses a thermostable restriction enzyme, BstN1, that destroys the PCR product from wild-type KRAS codon 12 sequences (25). REMS-PCR was performed essentially as described in (25), with the modification that a fluorescently labeled TaqMan probe was used to monitor the reaction in real time. Control experiments with genomic DNA from the KRAS-variant cell line SW480 and wild-type DNA demonstrate the presence and absence, respectively, of PCR product (see Fig S3A in the online Data Supplement). Dilution of the variant (SW480) into wild-type genomic DNA was then performed and REMS-PCR was applied. A linear change of the PCR threshold vs concentration of the variant DNA was observed down to 0.1% variant-to-wild-type ratio (see Fig. S3B in the online Data Supplement).
detection of kras sequence alterations in clinical samples using pna-pcr/ocean, pna-pcr/sequencing, or rems-pcr/sequencing
To examine the reliability of PNA-PCR/OCEAN in detecting KRAS codon 12 sequence variants in clinical samples, we screened 106 clinical lung surgical samples in triplicate independent experiments. We first performed a PNA-PCR reaction on these samples. When a product visible on gel electrophoresis was obtained, we processed the DNA for direct sequencing at the Dana Farber Cancer Institute sequencing Core facility to identify the sequence alteration, but when there was no product visible, we used regular PCR amplification followed by sequencing to verify the anticipated wild-type sequence. At the same time, we diagnosed the genotype of the clinical samples with PNA-PCR followed by fluorescence microplate-based OCEAN reactions with amplifiers addressing a different KRAS codon 12 sequence variant in each reaction. Representative results for clinical lung samples are depicted in Fig. 5
. We screened 5 samples harboring KRAS codon 12 sequence variants (TL3, TL30, TL42, TL37, and TL39) and 5 wild-type samples (TL55, TL72, TL92, NL15, and NL37) with amplifiers corresponding to the 4 most common sequence alterations encountered in lung samples. In each case, OCEAN generated a genotype-specific signal corresponding to the correct KRAS codon 12 sequence in agreement with the sequencing results. The 106 surgical samples were also screened via real-time REMS-PCR in triplicate experiments. Samples generating a real-time PCR signal (sequence variant-positive) were then processed for sequencing. The detailed results of the 3 independent methods (PNA-PCR/sequencing, PNA-PCR/OCEAN, and REMS-PCR/sequencing) are presented in Table S2 in the online Supplement Data, and a summary of results is presented in Table 1
. There was generally an agreement among the 3 independent methods. Of the 3 methods, PNA-PCR/sequencing, which does not require enzymatic steps for sequence alteration identification and therefore is assumed to be more robust, was taken to be the calibrated method against which the 2 other methods were compared. The concordance rate between PNA-PCR/OCEAN and PNA-PCR/sequencing was 98% but dropped to 93% when OCEAN was conducted blindly after PNA-PCR, i.e., without prior examination of the presence or absence of a visible band on agarose gel electrophoresis after the PNA-PCR step. In 5% of cases, there may have been nonspecific signals that exceeded the selected threshold of 1 000 fluorescence units. The concordance rate of real-time REMS-PCR/sequencing with PNA-PCR/sequencing was 93%. The 7% discrepancy could in part be because of failure of the restriction enzyme to digest efficiently during PCR, owing to the presence of impurities. In a few cases (e.g., sample TL60), both PNA-PCR/OCEAN and REMS-PCR scored KRAS sequence alterations, but the genotypes detected after sequencing were different, possibly because these clinical samples contained a mixture of wild-type and variant genomic DNA, and the selection efficiency of the 2 methods differed for the 2 genotypes, leading to distinct genotype calling. Overall, the 2 methods, PNA-PCR/OCEAN and REMS-PCR/sequencing agreed in 89% of cases.
|
|
| Discussion |
|---|
|
|
|---|
Genotyping technologies include methods with signal generation by the cleavage of branched DNA structures by the 5' nuclease activity of enzymes, such as TaqMan assays (1), in which specific probes are cleaved by Taq DNA polymerase, and Invader assays (34), in which Cleavase (a thermostable flap endonuclease) recognizes and cleaves a 3-dimensional branched complex. In the OCEAN method, the specific reaction leading to enzymatic digestion of the amplifiers leads to a loss of quenching monitorable on a fluorescence microplate reader.
Identification of KRAS sequence variant genotypes after PCR by OCEAN could circumvent the requirement for sequencing. Identification of specific genotypes of KRAS codon 12 sequence alterations in surgical tumor specimens or bodily fluids is necessary for differentiating benign from malignant tissues and determining disease prognosis. For example, in NSCLC patients, conversion of wild-type GGT to GTT or CGT has been associated with a poorer prognosis than conversion to TGT (8)(9). PNA-PCR/OCEAN is a low-cost approach for genotyping that can replace expensive sequencing technologies in small or modest-sized molecular diagnostics laboratories. We envision adaptation of the 3-way DNA digestion principle for detecting other somatic sequence variants [such as specific epidermal growth factor receptor deletions and single-nucleotide substitutions (46)] and genotyping blood-related disorders. In other applications such as single nucleotide variant genotyping, OCEAN does not require any additional selective procedures (unpublished data).
In conclusion, OCEAN is novel and practical technology for rapid and reliable postPCR detection of DNA sequence changes in human specimens. Combined PNA-PCR/OCEAN simultaneously detects the presence and type of KRAS sequence alterations in clinical samples. The method has high selectivity and sensitivity required for analysis of surgical/biopsy tumor samples and can be used for rapid diagnosis of KRAS sequence alterations in human tissue and, potentially, fluid specimens. We anticipate the use of this method for identification and genotyping of variants in other genes.
| Footnotes |
|---|
2 Human gene: KRAS, v-Ki-ras2 Kirsten rat sarcoma viral oncogene homolog. ![]()
| References |
|---|
|
|
|---|
The following articles in journals at HighWire Press have cited this article:
![]() |
A. Di Nicola, E. Ghezzi, F. Gillio, F. Zerilli, E. Shehi, D. Maritano, M. Panizzo, F. Bonelli, and D. Adlerstein Anchor-Based Fluorescent Amplicon Generation Assays (FLAG) for Real-Time Measurement of Human Cytomegalovirus, Epstein-Barr Virus, and Varicella-Zoster Virus Viral Loads Clin. Chem., November 1, 2008; 54(11): 1900 - 1907. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. Amicarelli, E. Shehi, G. M. Makrigiorgos, and D. Adlerstein FLAG assay as a novel method for real-time signal generation during PCR: application to detection and genotyping of KRAS codon 12 mutations Nucleic Acids Res., October 11, 2007; (2007) gkm809v1. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |