Open Access Article
Jian Gong
a,
Yishuai Lib,
Ting Linc,
Xiaoyan Fengc and
Li Chu*ad
aHebei Medical University, Hebei 050017, China. E-mail: gongjian2001@126.com
bDepartment of Thoracic Surgery, Hebei Chest Hospital, Hebei 050000, China. E-mail: liyishuai66@126.com
cApexbio Biotech Co., LTD, Beijing 100176, China
dDepartment of Pharmacology, Hebei University of Chinese Medicine, Hebei 050200, China. E-mail: chuli0614@126.com
First published on 1st August 2018
With the continuous development and application of targeted drugs, it is particularly desirable to find a non-invasive diagnostic approach to screen patients for precision treatment. Specifically, detection of multiple cancer-related mutations is very important for targeted therapy and prediction of drug resistance. Although numerous advanced PCR methods have been developed to discriminate single nucleotide polymorphisms, their drawbacks significantly limit their application, such as low sensitivity and throughput, complicated operations, and expensive costs. In order to overcome these challenges, in this study, we developed a method combining multiplex and sensitive real-time PCR assay with rolling circle amplification. This allows specific and sensitive discrimination of the single nucleotide mutation and provides convenient multiplex detection by real-time PCR assay. The clinical potential of the MPRP assay was further demonstrated by comparing samples from 8 patients with a digital PCR assay. The coincident results between these two methods indicated that the MPRP assay can provide a specific, sensitive, and convenient method for multiplex detection of cancer-related mutations.
Therefore, developing a multiplex, sensitive and convenient method to detect tumor-related mutations is an urgent need for clinical precision treatment. Liquid biopsies of blood or urine were recently established to be used for the diagnosis of cancer and other diseases. The clinical application of cell-free DNA, which contains the genetic information of tumors, took center stage.7 However, the low amount and shorter length of cell-free DNA limited the effective detection of tumor-related mutations.8 Moreover, sensitive detection of the mutation under the highly abundant background of wild-type DNA remains difficult. At present, several real-time PCR methods such as ARMS,9,10 COBAS,11,12 etc., have been developed to detect this mutation. Unfortunately, low sensitivity and false-positive results have limited their clinical application. Due to the increased number of approved targeted medications and the subsequent increase in DNA targets that need to be detected, more input samples need to be prepared. Excitingly, the padlock RCA method has been established to specifically identify single nucleotide polymorphisms (SNP).13 After complete annealing to the target DNA strand, two ends of the padlock probe can be ligated to circularize. A single mismatched base present at the 3′ terminal of the probe will abolish the ligation, thus ensuring high specificity in DNA and RNA detection. Then, rolling circle amplification (RCA), which is based on rolling circle replication of circular DNA pathogens in nature, was carried out to achieve exponential amplification of the circularized DNA.14,15 Combined with RCA technology, padlock RCA can detect SNP from 1 ng (300 copies) of genomic DNA and the accuracy rate of SNP genotyping results can reach 100%.16 However, the effectiveness of ligation and the stability of the reaction conditions can still be problematic. On the other hand, some methods require the removal of the residual probe which cannot hybridize to the target template, thereby increasing the complexity of detection.17 Hence, padlock RCA-based assay needs further improvement.
In this study, we presented the MPRP method to combine padlock RCA assay and real-time PCR technology to detect multiplex tumor-related mutations. We found that this combined method can specifically and sensitively distinguish the genetic mutation and allow convenient multi-detection of tumor-related mutations by real-time PCR. Furthermore, we demonstrated the clinical applicability of this combined method by multi-detection of EGFR mutations L858R and T790M and BRAF mutation V600E in patients' plasma.
| Name | Sequence (5′–3′) | |
|---|---|---|
| Padlock probe | pEG858A | 5′–P-gcccaaaatctgtgatcttgaacataggtctcagtccagccattttagccgttcctcaca tcagactcgcactcgttagcaatcactttttacccagcagtttggcca-OH–3′ |
| pEG858C | 5′–P-gcccaaaatctgtgatcttgaacataggtctcagtccagccattttagccgttcctca catcagactcgcactcgttagcaatcactttttacccagcagtttggccc–OH-3′ | |
| pEG790G | 5′–P-tgatgagctgcacggtggaggacataggtctcagtccagccattttagccgttcctc acatcagactcgcactcgttagcaatcactttttcagccgaagggcatgagctgcg–OH-3′ | |
| pEG790A | 5′–P-tgatgagctgcacggtggaggacataggtctcagtccagccattttagccgttcctca catcagactcgcactcgttagcaatcactttttcagccgaagggcatgagctgca–OH-3′ | |
| pBR600A | 5′–P-ctgtagctagaccaaaatcacctatacataggtctcagtccagccattttagccgttcct cacatcagactcgcactcgttagcaatcactttttggacccactccatcgagatttca–OH-3′ | |
| pBR600T | 5′–P-ctgtagctagaccaaaatcacctatacataggtctcagtccagccattttagccgttcct cacatcagactcgcactcgttagcaatcactttttggacccactccatcgagatttct–OH-3′ | |
| Template DNA for padlock probe | tBR600T | 5′-ataggtgattttggtctagctacagtgaaatctcgatggagtgggtcccat-3′ |
| tBR600A | 5′-ataggtgattttggtctagctacagagaaatctcgatggagtgggtcccat-3′ | |
| tEG790C | 5′-ctcacctccaccgtgcagctcatcacgcagctcatgcccttcggctgcctc-3′ | |
| tEG790T | 5′-ctcacctccaccgtgcagctcatcatgcagctcatgcccttcggctgcctc-3′ | |
| tEG858T | 5′-agcatgtcaagatcacagattttgggctggccaaactgctgggtgcggaagagaa-3′ | |
| tEG858G | 5′-agcatgtcaagatcacagattttgggcgggccaaactgctgggtgcggaagagaa-3′ | |
| Digital PCR primers | 858F | 5′-gcagcatgtcaagatcacagatt-3′ |
| 858R | 5′-cctccttctgcatggtattctttct-3′ | |
| 790F | 5′-gcctgctgggcatctg-3′ | |
| 790R | 5′-tctttgtgttcccggacatagtc-3′ | |
| BF | 5′-acctcagatatatttcttcatg-3′ | |
| BR | 5′-ccagacaactgttcaaac-3′ | |
| PCR probes | P858M | 5′-FAM-agtttggcccgcccaa-MGB-3′ |
| P858W | 5′-VIC-agtttggccagcccaa-MGB-3′ | |
| P790W | 5′-VIC-tgagctgcgtgatga-MGB-3′ | |
| P790M | 5′-CY5-tgagctgcatgatga-MGB-3′ | |
| PBW | 5′-CY5-tcgagatttcactgtagct-MGB-3′ | |
| PBM | 5′-VIC-tcgagatttctctgtagct-MGB-3′ | |
| Universal primers | uniF | 5′-tggctggactgagacctatgt-3′ |
| uniR | 5′-agccgttcctcacatcagac-3′ |
The successful circularization and specific discrimination of padlock probes depends on the activity of DNA ligase and the optimization of the reaction. T4 DNA ligase, which is one of the most popular DNA ligases used for joining the 5′ phosphate and 3′ hydroxyl terminal of duplex DNA or RNA,25 has the capability of SNP discrimination in RNA, but no capacity of SNP detection in DNA.24,26 HiFi Taq DNA ligase, which demonstrates strong capability of SNP discrimination in DNA, is quite suitable for the padlock probe.27
Moreover, a couple of factors affecting the padlock RCA reaction for SNP discrimination have been reported.28–30 Ligase temperature is one of the most important factors for the efficiency of ligation, which can be predicted using the Thermostable Ligase Reaction Temperature Calculator v0.8.4 tool from NEB (http://ligasecalc.neb.com/#!/ligation). In our experiment, 43 °C was selected as the ligase temperature. Moreover, the bases at the 3′ end of the padlock probe were LNA modified to enhance the base stacking of perfectly matched base pairs and decrease the stacking stability of the mismatched pairs.31 The optimal concentration of the padlock probe for padlock RCA was calculated to be 0.1 pM according to the Ct value of template DNA determined by real time PCR, which was chosen in the following experiments (Fig. S1†). In addition, the stability of the padlock RCA reaction was further improved by endonuclease digestion and betaine solution as previously reported.17,32
In order to investigate the feasibility and specificity of the MPRP system, three padlock probes were utilized to detect three different target DNAs with point mutations in a Padlock RCA reaction. For example, probe pEG858A and pEG858C were perfectly matched with target tEG858T and tEG858G, respectively. Similarly, probe pEG790G and pEG790A were perfectly matched with target tEG790C and tEG790T, while probe pBR600A and pBR600T were perfectly matched with target tBR600T and tBR600A, respectively (Table 1). Finally, one pair of universal primers and three different probes labeled with different fluorescent moieties were used to perform multiplex real time PCR. It was shown in Table 2 that in the presence of mutant DNA and corresponding padlock probes, the Ct value of mutant DNA obtained from the MPRP assay was able to reach 9.74 ± 0.565, 7.99 ± 0.143 and 15.27 ± 0.243, respectively (Fig. S2†). No valid Ct value was observed with those padlock probes complementary to the wild type DNA. The results demonstrated the feasibility and specificity of the MPRP assay to detect multiplex DNA point mutations.
| Sample type | Padlock probe type | Average Ct value (mean ± standard error) |
|---|---|---|
| a ND: not detected.b The Ct value beyond 35 means that the data is unbelievable and can be considered as undetected. | ||
| EGFR L858R mutation (tEG858G) | L858R MU | 9.74 ± 0.565 |
| L858R WT | NDa | |
| L858R RCA NTC | ND | |
| L858R NTC | ND | |
| EGFR T790M mutation (tEG790T) | T790M MU | 7.99 ± 0.143 |
| T790M WT | 36.98 ± 0.660b | |
| T790M RCA NTC | ND | |
| T790M NTC | ND | |
| BRAF V600E mutation (tBR600A) | BRAF MU | 15.27 ± 0.243 |
| BRAF WT | ND | |
| BRAF RCA NTC | ND | |
| BRAF NTC | ND | |
Besides the specificity, sensitivity also affects the application potential of the SNP detection approach. To evaluate the sensitivity of the MPRP assay for mutant DNA detection, we prepared 0.05%, 0.1%, 1% and 5% of mutant standard DNA templates including the EGFR mutation L858R, T790M and BRAF mutation V600E, through dilution with wild type DNA according to the copy numbers. Then, the multiplex real-time PCR assay was performed to semi-quantify the level of mutant DNA, which was negatively associated with the Ct value. As was observed in Table 3, the corresponding Ct values of 0.1%, 1% and 5% of the EGFR T790M mutation were determined to be 31.72 ± 0.319, 25.39 ± 0.044 and 18.74 ± 0.458, respectively. Consistently, the Ct value of 0.1%, 1% and 5% of the EGFR L858R mutation was 34.95 ± 0.139, 17.69 ± 0.127 and 18.02 ± 0.234, respectively. The MPRP assay could detect 1% and 5% of the BRAF V600E mutation with Ct values of 20.85 ± 0.057 and 23.94 ± 0.063, respectively, but was not sensitive enough to detect the lower dilutions of the BRAF V600E mutation and the 0.05% of the two EGFR mutations (Fig. S3†).
| Gene type | Mutation ratio | Average Ct value (mean ± standard error) |
|---|---|---|
| a ND: not detected.b The Ct value beyond 35 means that the data is unbelievable and can be considered as undetected. | ||
| EGFR T790M | 5% | 18.74 ± 0.458 |
| 1% | 25.39 ± 0.044 | |
| 0.1% | 31.72 ± 0.319 | |
| 0.05% | 38.36 ± 0.091b | |
| EGFR L858R | 5% | 18.02 ± 0.234 |
| 1% | 17.69 ± 0.127 | |
| 0.1% | 34.95 ± 0.139 | |
| 0.05% | 38.97 ± 0.389b | |
| BRAF V600E | 5% | 20.85 ± 0.057 |
| 1% | 23.94 ± 0.063 | |
| 0.1% | NDa | |
| 0.05% | ND | |
After enrichment of the single-base mutation, multiplex real time PCR is used for the quantification of gene mutations, which further increases the specificity of SNP detection. Dual mutation discrimination mode prevents the false positive results that can be caused by high sensitivity. Furthermore, instead of the traditional fluorescent detection method utilizing a fluorescent reader or microscope that is difficult to calculate and has low sensitivity,17,24,26 the real time PCR method provides automatic detection of signals which can be calculated and presented with visual results using professional analysis software. In our study, we detected the mutation with the mutant ratio as low as 0.1%, which is much better than the sensitivity of traditional real time PCR such as ARMS-PCR.33,34
| Patients ID | Mutation type | Mutation detected | |
|---|---|---|---|
| MPRP assay (Ct value) | Digital PCR (mutation ratio) | ||
| a ND: not detected.b The Ct value beyond 35 means that the data is unbelievable and can be considered as undetected. | |||
| P6 | EGFR L858R | 36.76 ± 0.783b | 2.38% |
| EGFR T790M | 27.09 ± 0.525 | 7.59% | |
| BRAF V600E | NDa | 0 | |
| P3 | EGFR L858R | 24.37 ± 0.120 | 4.21% |
| EGFR T790M | 20.97 ± 0.182 | 8.25% | |
| BRAF V600E | ND | 0 | |
| P2 | EGFR L858R | ND | 0 |
| EGFR T790M | 12.62 ± 0.091 | 25.98% | |
| BRAF V600E | ND | 0 | |
| P9 | EGFR L858R | ND | 0 |
| EGFR T790M | 38.31 ± 1.046b | 0 | |
| BRAF V600E | ND | 0 | |
| P5 | EGFR L858R | 21.00 ± 0.059 | 2.30% |
| EGFR T790M | 37.58 ± 1.358b | 0 | |
| BRAF V600E | ND | 0 | |
| P1 | EGFR L858R | ND | 0 |
| EGFR T790M | 36.80 ± 0.128b | 0.26% | |
| BRAF V600E | ND | 0 | |
| P7 | EGFR L858R | ND | 0 |
| EGFR T790M | 36.80 ± 0.801b | 0 | |
| BRAF V600E | ND | 0 | |
| P8 | EGFR L858R | 37.48 ± 0.399b | 0 |
| EGFR T790M | ND | 0 | |
| BRAF V600E | ND | 0 | |
According to our knowledge, it's an innovative strategy to combine padlock RCA with multiplex real-time PCR for the non-invasive detection of lung cancer-related mutations, which appears to be very promising in clinical application.
000 g for 10 min. The supernatant was then transferred to fresh tubes and stored at −80 °C. Genomic DNA was extracted from the cells using a DNA Extraction Kit (Apexbio, Beijing, China) according to the user manual. Circulating DNA from plasma was extracted with the QIAamp Circulating Nucleic Acid Kit (Qiagen, Valencia, CA) according to the manufacturer's protocol. DNA quantification was performed in a Colibri microvolume spectrophotometer (Titertek-Berthold, Pforzheim, Germany).
For the padlock RCA reaction of specificity test, the following were added into three PCR tubes for EGFR L858R, EGFR T790M and BRAF V600E, respectively: 1 μL of each 0.1 pM corresponding padlock probe (both of mutant type and wild type), 1 μL of 0.1 pM corresponding mutant template which was synthesized, 22 μL of ligation solution containing 0.5 μL HiFi Taq DNA ligase (NEB, Ipswich, MA), and 1× HiFi Taq DNA ligase reaction buffer. For the padlock RCA reaction of sensitivity test, the following were added into the PCR tubes for each dilution series of standard DNA: 1 μL of each 0.1 pM padlock probe of mutant types (pEG858C, pEG790A and pBR600A), 1 μL of 0.1 pM standard DNA, 21 μL of ligation solution containing 0.5 μL HiFi Taq DNA ligase (NEB, Ipswich, MA), and 1× HiFi Taq DNA ligase reaction buffer. For the padlock RCA reaction of clinical test, the following were added into the PCR tubes for each sample: 1 μL of each 0.1 pM padlock probe of mutant types (pEG858C, pEG790A and pBR600A), 5 μL of plasma DNA, 17 μL of ligation solution containing 0.5 μL HiFi Taq DNA ligase (NEB, Ipswich, MA), and 1× HiFi Taq DNA ligase reaction buffer. The mixture was incubated at 55 °C for 1 hour after being heated at 95 °C for 3 min. After the circularization of the padlock probe, 10 U exonuclease I (NEB, Ipswich, MA) and 40 U exonuclease III (NEB, Ipswich, MA) were added to digest the ssDNA and dsDNA, respectively. Then, the RCA reaction was immediately performed to elongate the padlock probe circulated for subsequent real time PCR assay in 25 μL of 1× Thermopol buffer (50 mM Tris–HCL, 10 mM MgCL2, 10 mM (NH4)2SO4, 4 mM DTT, PH 7.5@RT), 4 μL 5 M betaine (ThermoFisher, USA), 14 mM dNTPs mix, 1 μL 100 μM uniF, 1 μL 100 μM uniR and 5 μL ligation products. After heating at 95 °C for 3 min and cooling down in ice, 1× BSA, 0.5 μL betaine, 6 mM MgSO4 and 10 U μL−1 of Bst DNA polymerase (NEB, Ipswich, US) were added. The mixture was incubated at 65 °C for 2 h before heating at 85 °C for 10 min to deactivate the polymerase. Then, quantitative real-time PCR testing was performed in a 20 μL reaction containing 1× PerFecTa Multiplex qPCR ToughMix (Quanta Biosciences, Gaithersburg, USA), 400 nM of each uniF and uniR, 200 nM of each multiplex probes, and 1 μL RCA product. For specificity test, the PCR probes included both wild type and mutant type for EGFR L858R, EGFR T790M and BRAF V600E, respectively. For sensitivity test and clinical test, the PCR probes included three probes of mutant type (P858M, P790M and PBM) simultaneously. This was programmed as: 95 °C for 3 min, followed by 40 cycles of 94 °C for 10 s, 60 °C for 20 s on a CFX96 real-time PCR instrument (Bio-Rad, Hercules, CA). All experiments were replicated to ensure reproducibility.
Footnote |
| † Electronic supplementary information (ESI) available. See DOI: 10.1039/c8ra05259j |
| This journal is © The Royal Society of Chemistry 2018 |