Next Article in Journal
Machine Learning Aided Photonic Diagnostic System for Minimally Invasive Optically Guided Surgery in the Hepatoduodenal Area
Previous Article in Journal
Personalization of Aspirin Therapy Ex Vivo in Patients with Atherosclerosis Using Light Transmission Aggregometry
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Allele-Specific PCR for KRAS Mutation Detection Using Phosphoryl Guanidine Modified Primers

by
Alexey S. Chubarov
*,
Igor P. Oscorbin
,
Maxim L. Filipenko
,
Alexander A. Lomzov
and
Dmitrii V. Pyshnyi
*
Institute of Chemical Biology and Fundamental Medicine, SB RAS, 8 Lavrentiev Avenue, 630090 Novosibirsk, Russia
*
Authors to whom correspondence should be addressed.
Diagnostics 2020, 10(11), 872; https://doi.org/10.3390/diagnostics10110872
Submission received: 29 September 2020 / Revised: 21 October 2020 / Accepted: 22 October 2020 / Published: 26 October 2020
(This article belongs to the Section Pathology and Molecular Diagnostics)

Abstract

:
Establishing the Kirsten rat sarcoma (KRAS) mutational status is essential in terms of managing patients with various types of cancer. Allele-specific real-time polymerase chain reaction (AS-PCR) is a widely used method for somatic mutations detection. To improve the limited sensitivity and specificity, several blocking methods have been introduced in AS-PCR to block the amplification of wild-type templates. Herein, we used a novel modified oligonucleotide with internucleotide phosphates reshaped 1,3-dimethyl-2-imino-imidazolidine moieties (phosphoryl guanidine (PG) groups) as primers and blockers in the AS-PCR method. Four common KRAS mutations were chosen as a model to demonstrate the advantages of the PG primers and blockers utilizing a customized PCR protocol. The methods were evaluated on plasmid model systems providing a KRAS mutation detection limit of 20 copies of mutant DNA in a proportion as low as 0.1% of the total DNA, with excellent specificity. PG-modification can serve as the universal additional mismatch-like disturbance to increase the discrimination between wild-type and mutated DNA. Moreover, PG can serve to increase primer specificity by a synergetic effect with additional mismatch and would greatly facilitate medical research.

Graphical Abstract

1. Introduction

Somatic mutations play a significant role in oncogenesis. They are routinely analyzed for diagnosis and treatment at clinical laboratories worldwide. The Kirsten rat sarcoma (KRAS) is a crucial gene in the development of a variety of cancers [1]. The mutation status of this gene is an important characteristic of many cancer types such as pancreatic, colon, lung and ovarian tumors [1,2]. The most common KRAS-activating mutations arise in codons 12, 13, 61, 117 and 146 [3]. According to the COSMIC database (the Catalogue Of Somatic Mutations In Cancer), a majority (approximately 90–95%) of KRAS-activating mutations occur in codons 12 (p.G12A, p.G12C, p.G12D, p.G12R, p.G12S, and p.G12V) or 13 (p.G13D) of exon 2, whereas less common mutations occur in exon 3 and exon 4. Among the G12 codon mutations, the most frequent ones are G12D and G12V in colon and rectal cancers [3].
Establishing the KRAS mutational status of tumor samples is essential in terms of managing patients with various types of cancer. Predominantly, there are three methods including nucleic acid sequencing, real-time PCR with melt-curve analysis and allele-specific (AS) real-time polymerase chain reaction (AS-PCR), that have been utilized to detect mutations in the KRAS gene [1,2]. At present, sequencing technologies range from classical di-deoxy-Sanger sequencing to next-generation sequencing (NGS) for mutation detection [4,5,6,7]. NGS technology provides high-throughput and base pair-resolution data, good sensitivity and permits the analysis of multiple mutations simultaneously [7,8,9]. However, NGS has a much longer analysis time and relatively high per-sample costs. Thus, it is not cost-effective for routine determination of KRAS in a tumor sample [7,8,9]. PCR with melt-curve analysis is a rapid turnaround, low cost and less labor-intensive method [10]. However, it can detect mutations in samples containing 1% to 10% cells with mutated DNA [1,10].
In AS-PCR, wild-type or mutated DNA molecules are selectively amplified using an AS primer, where the nucleotide substitution site is located at the 3′-end of the primer, i.e., the amplification efficiency of the mismatch allele is lower than that of the perfect match allele. The advantages of this method are low cost, high sensitivity (~1% of allele frequency is detectable) and single closed tube performance (preventing contamination). However, the use of AS-PCR in KRAS mutation detection is restricted by low discrimination efficiency between mutant and wild-type DNA for several types of mismatches [2,11]. A further limitations is that only targeted mutations can be determined, with lack of specificity driven by combining mismatch bases and formation of primer dimers. Therefore, laborious optimization of the AS-PCR conditions and primer structure is desirable for each new mutation [12]. To increase detection sensitivity and specificity for the mutation, various blocking oligonucleotides have been introduced in AS-PCR to decrease the efficiency of wild-type DNA amplification [2,12,13,14].
Over the previous decade, nucleic acid analogues have been extensively described to improve diagnostic analysis by PCR. Several of them, including peptide nucleic acids (PNA), locked nucleic acids (LNA), phosphorothioate, zip nucleic acids (ZNA), methyl phosphotriester oligonucleotide analogues (MPTE), nonnatural bases isoguanine (iG) and 5′-methylisocytosine (iC), are currently used for real-time PCR applications [2,15,16,17,18,19,20,21]. For instance, LNA is widely used for the detection of mutations with miRNA and DNA methylation [2]. Several of these modifications (MPTE, PNA, LNA) due to their physicochemical properties (a neutral charge, conformation, etc.) can bind to the complementary single-stranded nucleic acids. Therefore, they have higher thermal stability of duplexes [2,16,22]. Recently, a new type of uncharged phosphate-modified nucleic acid analog, namely phosphoryl guanidine oligonucleotides (PGOs), has been described [23,24]. A phosphoryl guanidine group may be easily introduced into any desired position of an oligonucleotide by standard automated DNA synthesis using phosphoramidite chemistry. In this class of compounds, replacement of the negatively charged oxygen atom with an electroneutral 1,3-dimethylimidazolidine-2-imine group (Figure 1) increases the hydrophobicity of the oligonucleotide while maintaining good water solubility, forms thermodynamically stable duplexes with DNA, but does not disturb the DNA double helix structure [25]. Previously, PGOs were shown to exhibit resistance to some nucleases [26] and can be used as promising tools for different biological applications, including acting as primers for DNA amplification in PCR [27,28,29,30,31,32,33].
In the current work we investigated the possibility of using PGO analogs as primers for AS-PCR and for the allele-specific competitive blocker PCR (ASB-PCR) to increase discrimination efficiency between mutant and wild-type DNA. Several most common KRAS mutations in codons 12 and 13 (p.G12A, p.G12D, p.G12V, and p.G13D) were used as models to validate the concept. The protocol described herein is standardized with AS-PCR using TaqMan probes for amplification detection. The sensitivity and selectivity of PCR assays have been evaluated on plasmid templates, harboring the aforementioned KRAS mutations or the wild-type KRAS fragment. It has been discovered that primers with one and three PG-modifications are suitable for AS-PCR and ASB-PCR, respectively. Methods are straightforward in assay optimization and can achieve 0.1% sensitivity with excellent specificity in the detection of KRAS mutations. The cumulative effect of the PG-modified primers and PG-blockers has shown great potential for applications in ASB-PCR.

2. Materials and Methods

2.1. Synthesis and Isolation of Oligonucleotides

Oligonucleotides were synthesized in an ASM-800 automated synthesizer (Biosset, Russia) according to the standard protocol of the 2-cyanoethyl phosphoramidite method using commercially available deoxyribonucleoside monomers and appropriate porous glass (Glen Research, Sterling, VA, USA). Oligonucleotides containing phosphoryl guanidine units were synthesized using the protocol described previously by LLC NooGen [23,24]. Oligonucleotides were purified by RP-HPLC performed in the Agilent 1200 series chromatograph (Agilent, Santa Clara, CA, USA) on a column (4.6 × 150 mm) containing the Eclipse XDB-C18 sorbent (5 μm) (Agilent, USA) with a 0–90% linear gradient of acetonitrile concentration in 0.02 M triethylammonium acetate solution for 30 min at a flow rate of 1.5 mL/min. Fractions containing the target product were evaporated in vacuo. The bulk of triethylammonium acetate was removed by coevaporations with ethanol. To remove the protecting dimethoxytrityl group, the isolated oligonucleotides were treated with 80% acetic acid (25 °C, 7 min). Purified oligonucleotides were concentrated following by precipitation with 2% LiClO4 in acetone, washing with pure acetone, and desiccation under vacuum. After desiccation, the oligonucleotides were dissolved in 0.1 mL of deionized water and stored at –20 °C. The concentration of oligonucleotides in stock solutions were determined in accordance with the working process described above [26].

2.2. Plasmid Standards

The control plasmids contained a partial sequence of the wild-type KRAS gene or KRAS gene with mutations in codons 12 and 13 (p.G12A, p.G12D, p.G12V and p.G13D) serving as positive controls, were used to assess method sensitivity and were constructed by Shanghai RealGene Bio-tech, Inc (Shanghai, China). Before use, all control plasmids were purified, linearized by digestion with BamHI restriction endonuclease and quantified using NanoDrop Lite A4 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA).

2.3. Real-Time PCR

Real-time PCR assays were performed in 20 µL containing 1× PCR-buffer (65 mM Tris–HCl, pH 8.9, 24 mM (NH4)2SO4, 0.05 % Tween-20, 3 mM MgSO4), 0.2 mM dNTP, 450 nM primers, 0–2250 nM blockers, 100 nM fluorescent hydrolysis probe, DNA template (exact amount indicated below), and 1 U of Taq-polymerase (Biosan, Novosibirsk, Russia). A control plasmid was used as a DNA template at the concentration indicated below. Amplification was carried out in a CFX96 Real-Time PCR Detection System (Bio-Rad, Hercules, CA, USA) according to the following program: 95 °C for 3 min followed by 45 cycles of 95 °C for 10 s, and 60 °C for 40 s with a collection of fluorescent signals at the FAM channel. Reactions were carried out at least in triplicate and performed several times on separate occasions.
Average Cq ± standard deviation (SD) values are given in the tables. PCR analysis was performed using a reverse primer (k-rev) 5′-CATATTCGTCCACAAAATGATTCTG-3′, probe 5′-FAM-CTGTATCGTCAAGGCACTCTTGC-BHQ1-3′ and a series of forward primers 5′-AAACTTGTGGTAGTTGGAGXXXX-3′ (codon 12), 5’-GTGGTAGTTGGAGCTGXXXX-3′ (codon 13). XXXX—four nucleotides of the primer 3′-terminus, which are presented in the text as an abbreviation of the whole primer. Boldly marked nucleotides represent mismatched nucleotides in the relation to the wild-type DNA sequence. The symbol “*” in primers sequence means phosphate group modified with 1,3-dimethylimidazolidine-2-imine residue (phosphoryl guanidine (PG) modification). PG modification presented in Figure 1.
In each PCR plate forward primer 5′-GACTGAATATAAACTTGTGGTAGTTG-3′ was used as a reference primer (k-ref) to compare the data between various PCR experiments. Corresponding ΔCq values were calculated and used for the further analysis of the primer’s efficacy. NTC—no template control.

3. Results and Discussion

3.1. Effect of the Number of PG Groups and Their Position from 3′ End of the Primer on the AS-PCR Results

Allele-specific amplification, combined with fluorescence probe real-time polymerase chain reaction (real-time AS-PCR), has been widely used for detecting genetic variants, mutations, or single nucleotide polymorphisms (SNPs). The region of the DNA template that interacts with a polymerase via the minor groove of the duplex structure is about 6-8 nucleotides long and is located at the 3′-terminus of the primer strand [34]. In principle, a mutation can be detected using AS-PCR primers based on the 3′-terminal nucleotide of a primer that corresponds to a specific mutation site. To obtain reliable discrimination between wild-type (WT) DNA and DNA with a mutation, AS primers usually have additional nucleotide substitutions within a 2–4 nucleotides region closest to the 3′-end [11,35,36]. To examine the usefulness of phosphoryl guanidine oligonucleotides (PGOs) as primers for AS-PCR, we utilized previously established sequences for KRAS mutation detection [35,37]. PCR assays were evaluated on plasmids harboring KRAS gene fragments of wild-type or with four of the most common mutations in 12-13 codons (p.G12A, p.G12D, p.G12V, and p.G13D).
The influence of the number of PG groups, and their position at 3′-terminus, were investigated using AS-primers CTTC which showed in preliminary experiments the best level of discrimination when used without any modifications (Table 1) and with WT and KRAS G12A mutation plasmids as templates. Bold symbols marked nucleotides represent mismatched nucleotides in relation to the WT DNA sequence, here T/G and C/G substitutions (WT DNA sequence CTGG).
It was discovered that four (*C*T*T*C) and three (C*T*T*C) PG groups at the 3′-end of the primer fully inhibited PCR. The latter fact is in a good correlation with previously obtained results of the authors [32,33]. The Cq value for two PG-modified primers depended on the modification position from the primer’s 3′-end and decreased in the following order: CT*T*C ≥ *CTT*C > C*T*TC > *C*TTC > C*TT*C (Table 1). However, for most of the AS-primers Cq was higher than 35 cycles which indicated an unacceptably low PCR efficiency and is not recommended for use in AS-PCR. The Cq value for single PG-modified primers was much lower. The distancing of the modification from the 3′-end of the primer led to a decrease in Cq value. No amplification with wild-type KRAS plasmid as a template was observed for all PG-modified primers, indicating an excellent specificity. The lowest Cq values with appropriate selectivity were discovered for the primers with PG modification at the third and/or fourth internucleoside phosphate (Table 1). For comparison, the same results were obtained using the primer sequence CTGC without additional mismatch at the second nucleotide from the 3′-end (Table S1).

3.2. Allele-Specific PG-Modified Primer Improves Discrimination of Mutant and WT DNA

As mentioned above, the major limitation of AS-PCR is a lack of primer specificity to discriminate between WT and mutant DNA for some types of mismatches. In the case of somatic mutation detection, the situation is significantly complicated by the heterogeneity of the sample, which may contain both normal and tumor cells leading to low-frequency mutant DNA in the abundant WT DNA background. To investigate the possibility of increasing specificity to discriminate between WT and mutant DNA, primers with PG modification at the third and/or fourth internucleotide phosphate were used as they demonstrated the best performance (lowest Cq values) in the previous experiment. As a control for AS-PCR, a reference primer (k-ref) was used and also primers that had previously been used for KRAS mutation detection (CGGA, GAGA, CTAT, CTTC) [35,37]. During the experiment, the portion of mutant plasmid DNA was reduced to obtain decreased ratios of mutant to WT DNA, based on a total DNA amount of 2 × 104 copies per reaction (Table 2). Primers with PG modification at fourth phosphate from 3′-end usually have significantly higher Cq values (cf. C*TAT and *CTAT, Table 2) which are more than 40 cycles. This could be caused by impaired interaction of Taq-polymerase thumb amino acids (Arg536 and Ser515) with PG-modified phosphate residue in the fourth position [34]. This may lead to the significantly lower coordination of the modified nucleic acid complex by protein. Such high values of Cq do not apply to AS-PCR. For the primers with PG modification at the third internucleotide phosphate for each mutation detection, ΔCq(WT−1%) values (ΔCq = Cq(WT) − Cq(1%) for the same primer, Cq(WT)—cycle quantification value obtained using WT DNA, Cq(1%)—cycle quantification value obtained using 1% mutant DNA on the background of WT DNA) were much higher than the native primers, which indicates better discrimination efficiency. However, in the case of the primer G*AGA the Cq value was three cycles higher than for the primer without the PG group (Table 2). Due to this case, the primer’s sequence design is desirable.

3.3. Simultaneous Efficiency of PG Modification and Mismatches in the Primers Structure

Choosing additional mismatches to increase primer specificity has been a challenge for AS-PCR, and criteria for designing primers are still unclear. To identify the influence of the simultaneous effect of the mismatch and PG modification, mismatched nucleotides at second and/or third position from the 3′-end of the primers were taken randomly (Table 3). It had been described before [34] that an additional mismatch at the third position from the 3′-end proved to be the most effective in increasing the discrimination efficiency without significant loss of PCR efficiency (unpublished data). Nevertheless, the second nucleotide position can be effective in a particular nucleotide context [36]. To evaluate sensitivity and selectivity, DNA from mutant plasmids was diluted into WT plasmids (total DNA amount of 2 × 104 copies per reaction). The proportion of mutant plasmid DNA was gradually reduced to obtain decreasing ratios of mutant to WT DNA.
For G12A mutation detection, an increase occurred in discriminating ability with the introduction of an additional mismatch: 1.5, 5.7, 4.7 of ΔCq(WT−1%) for CTGC, CTTC, CAGC primers, respectively. For G12V, G12D and G13D mutations detection, the same results were registered. This indicates the mutation detection additional mismatch within the 2–3 nucleotides closest to 3′-end is preferable. However, the second mismatch gave the improvement only in several cases (Table 3). Remarkably, that for the primer with PG modification C*TGC, the ΔCq(WT−1%) value was 5.8 versus 1.5 for the CTGC for G12A detection. In this case, the PG group worked better for specific mutant DNA detection than primers with two mismatched nucleotides, such as CTTC or CAGC. However, in terms of the primers G*TGA (G13D), C*TGT (G12V) and C*TGA (G12D), there was no significant effect of the PG modification. A possible explanation is that the obtained results principally depended on the 3′-end mismatch. Polymerases reflect decreasing discrimination of 3′-end nucleotide mismatches in the following order: Pur/Pur > Pyr/Pyr > Pur/Pyr = Pyr/Pur [11,34]. However, there have been exceptions to this general rule reported [34], which could be explained by different thermodynamic destabilization effects that strongly influence primer elongation efficiency [36]. For instance, [34] Taq DNA-polymerase doesn’t efficiently elongate a primer with a Pyr/Pyr like С/С mismatch at the 3′-end of the primer (cf. Ct values of CTGC and CTGA, CTGT, GTGA primers, Table 3).
The simultaneous effect of the PG group and mismatches leads to the extraordinary discrimination effect for C*TTC, and only little benefit for the C*AGC. ΔCq(WT−1%) values for C*TTC and C*AGC were 12.3 and 6.2, respectively (cf. ΔCq for the control primer CTTC was 5.7). Better results for the C*TTC can be explained by the SNP type. It is known [36] that for SNP type G/C (G12A mutation) the mismatches in the second base pair closest to 3′-end of primer shows the highest discrimination efficiency, followed by third. For the SNP type G/A (G12D, G13D), the best location of the mismatch was in the third site and showed the highest discrimination efficiency [36], which was confirmed by the data in Table 3. Obtained results of the WT and the mutant allele discriminating effect for G12A, G12D and G13D mutations were in a good correlation of additional mismatch and PG modification effects in the appropriate site.
For the SNP type G/T (G12V mutation) the mismatches in the third site should have the highest discrimination efficiency followed by the second and fourth sites. However, in the case of G12V mutation, additional mismatch or PG modification in the third internucleotide phosphate didn’t lead to any significant effect of the mismatch discrimination capability. These results are likely to be related to the strong destabilization strength of mismatch types at the 3′-end of primer for G12D, G13D and G12A, while for the G12V mutation the effect was weak [36]. Therefore, primers generated from weak mismatch type C/T with mutated DNA (G12V mutation) had the lowest discrimination efficiency (Table 3). Of all the KRAS mutations, G12V is considered to be one of the most frequent KRAS codon 12 mutations in colorectal cancer patients [3]. Moreover, G12V has been related to a more aggressive phenotype, a worse progression of the disease and a short time of survival [38]. Increasing the specificity of its detection is very important for improving clinical outcome. For the KRAS G12V mutation detection primer CTAT, ΔCq(WT−1%) values were only 1.5,which indicates great difficulty to achieve the appropriate results. By the simultaneous effect of the two unpaired nucleotides and PG modification (C*TTT, C*TAT), much higher ΔCq(WT−1%) values (3.9 and 5.5, respectively) were obtained. However, in comparison with other ΔCq(WT−1%) values for the best PG-primers for other mutations obtained for G12V, the data is not so impressive, which can be explained by the weak destabilization strength of the 3′-end mismatch [36].
In qPCR, amplification efficiency is a very important parameter which can be calculated from a simple linear regression of the Cq values plotted against the log of the initial copy number. As a consequence, it is recommended that the following descriptors of the standard curve are reported for qPCR amplification: amplification efficiency (E), the linear regression coefficient (R2) and especially the y-intercept value, which uniquely describes the standard curve and indicates the sensitivity of the reaction. Standard curves show the PCR efficiency values for primers for G12A mutation detection generated from the dilution series of mutant DNA, e.g., for the C*TTC primer E = 98.3% and R2 = 0.999, y-intercept = −3.3641 (Figure S1). The same results were obtained for G12V, G12D and G13D mutations detection (Figures S2 and S3).
This discovery, while preliminary, suggests that PG-modification can be used as a universal additional mismatch-like disturbance of a primer/DNA complex with good discrimination properties of DNA polymerase action. Primers used in the present experiments contained only one PG-modified nucleotide at their 3′-end. Moreover, the PG group may have served to increase primer specificity by the synergetic effect with additional mismatch and would greatly facilitate and increase the reliability of clinical diagnostics. The PG-modified primers may be synthesized using standard equipment and could be potentially obtained from commercial oligonucleotide suppliers in a low-cost way.

3.4. Detection Limit and Specificity Investigation of the Synthesized PG-Modified Primers

Using plasmid DNA, mutation detection was achieved as low as 1% for each mutation. The best results in discriminating the wild-type and the mutant allele and low Cq values were found for primers C*TAT, C*TTT (G12V mutation), С*TTC (G12A), C*GGA (G12D), G*AGA, G*CGA (G13D) (Table 3). Typical amplification curves obtained from allele-specific PCR of the native and PG modified primers are presented in Figure 2.
To determine assay reproducibility, several experiments were performed in five replicates with various amounts of the mutant template in the samples (SD value~±0.5, Table S2). The next section of the survey was concerned with extended research of the primers with reliable WT discrimination efficiency. The primers were tested to achieve lower than 1% detection sensitivity with good specificity in the abundant wild-type DNA background. A standard real-time PCR protocol and PG modified primers were used to obtain ΔCq(WT−0.1%) values (Table 4). The best result were observed for both KRAS G12A and G13D followed by G12D mutation assays, which showed 0.1% detection sensitivity (~20 copies per reaction) without nonspecific amplification on WT DNA. However, the absolute Cq values for all modified primers were ≥35 cycles. It should be noted that for the primers with and without PG modification сlose Cq values were obtained in most cases. These data suggest that 0.1% of mutant DNA detection can be achieved through assay optimization.

3.5. Improving AS-PCR Specificity by Wild-Type DNA Blocking PCR

Previous studies [32,33] have shown an ability of PGO to terminate enzymatic DNA synthesis. In the present work, a simple approach was described using PG-modified primers to facilitate the mutations detection by blocking WT DNA amplification. The concept of this technique is similar to classical wild-type blocking PCR assays (ASB-PCR) (Figure 3). This approach was tested using several KRAS mutations. As AS-primers, the best in discriminating the WT and the mutant allele primers were chosen (Table 4). The blocking primers had the sequence C*T*G*G for 12 codon and G*T*G*G for 13 codon, which are very similar to AS-primers but lack any nucleotide substitutions at the 3′-end. Therefore, blocker oligonucleotides are perfectly matched to the WT DNA.
The competition-based ASB-PCR approach inherently suffers from lower PCR efficiency with the stronger WT DNA suppression. To obtain good discrimination between wild-type and mutated DNA for a certain type of blocker, very high 10 [39] or 20-fold [40] excess over AS-primers is required. For other systems such as T-blockers [41] or PNA [39], five and 1.2 blocker excesses are sufficient. In the current work, several combinations of AS-primers and blockers were tested to illustrate the largest Cq separation between the mutant and WT reactions. In the reaction mixture, the concentration of AS-primer was constant. The blocking primer performance was sensitive to the amount of blocking primer added. More than equimolar excess of the PG-blocker to an AS-primer inhibited PCR (Table S3). Therefore, fixed primer to blocker ratios 1:0.5, 1:0.25, and 1:0.1 were studied (Table S3). As the AS-primer:blocker ratio increased, the difference between 1% mutant and WT DNA became more pronounced and could clearly be distinguished when the ratio reached 1:0.25 as seen in Table 5 and Table S3. The ratio 1:0.5 partly inhibited PCR for most primers. Keeping blocker concentration at the minimum required amount is a good rule. The best results were obtained using primer/blocker ratios 1:0.25 and/or 1:0.1 (Table 5 and Table S3). In most cases, the system had higher discrimination efficiency and similar Cq values as the reaction without a blocker. For instance, for the primer CTTC, Cq(1%) values were 31.5 and 33.5 (with blocker). ΔCq(WT−1%) values were 5.2, and much higher with the blocker (11.5) which shows increasing mutation selectivity detection. Using 2 × 104 copies of DNA per reaction, no amplification was observed. The (Cq(WT) value was N/A, which made the calculation of ΔCq(WT−1%) impossible. Taking a higher number of recombinant DNA molecules (2 × 105 copies) the Cq value could be calculated (Table 5). Typical amplification curves obtained from allele-specific blocking PCR assays are presented in Figure 4.
Using primers with PG modification, it is possible to reach better discrimination ability of a PCR. For С*TTC (G12A mutation) and C*TTC/C*T*G*G, ΔCq(WT−1%) values were 12.4 and up to 14.2, respectively. At the same time, Cq(1%) values for CTTC and C*TTC/C*T*G*G were almost the same: 31.5 and 32.5, respectively. The same results were obtained for the GAGA primer/blocker system. For the worst discriminated C12V mutation, the primer/blocker assay gave an excellent result in discriminating the WT and the mutant allele. ΔCq(WT−1%) values for CTAT, C*TAT, C*TAT/C*T*G*G were 1.0, 5.5 and 9.0, respectively. These data present clear supporting evidence that PG-blockers play a strong WT DNA suppression role. The cumulative effect of the PG-modified primers and PG-blockers has great potential for applications in ASB-PCR.
Blocking methods have been applied successfully in several studies [2,12,13,14,41] to increase detection sensitivity and specificity for mutation detection. High sensitivity has been achieved in the assay, reaching one copy of mutant plasmid DNA to 1000 wild-type DNA dilution, corresponding to a fractional abundance of 0.1% (Table 5). ΔCq(WT−0.1%) values for the primers/blockers system are ranked in the 5–9 cycles range, which indicates good WT suppression. Recent studies [13,41] have shown sensitivity improvement down to 0.1% of mutant DNA by adopting various primer/blocker ratios, too. Compared to the assays presented in other works [13,41], PG-blockers have some advantages. It is not necessary to design the blocker structure since it can be used by the full complementary to WT DNA primer with three 3′-terminal PG-modifications. Cq values of 0.1% detection have been 37–39 [41] and 36–37 [13] which are comparable to the values in the present work.
Several studies have detected 0.01% of mutant DNA for the detection of KRAS using a highly optimized WT blocking PCR approach [42,43] or a combination of PCR and direct sequencing methods [44]. Comparing the results for the PG-primers and PG-primers/PG-blockers, it can be observed from Table 4 and Table 5 that it is possible to obtain 0.1%, and a bit lower detection sensitivity, with excellent discrimination between mutant/WT DNA by ASB-PCR assay. However, we failed to reach 0.01% sensitivity with sufficient mutant/WT DNA discrimination due to the amount of DNA being close to the theoretical limit of qPCR detection. The differences of the CqWT and Cq0.01% show the potential to obtain 0.01% sensitivity in the presence of more DNA copies per reaction, which represents a remarkably artificial case during analyzing natural samples. As colorectal cancer (CRC) is the most studied with regard to KRAS mutation, due to the implementation of anti-EGFR (epidermal growth factor receptor) antibodies into clinical practice, it has been considered that around 2% of patients with metastatic CRC truly harbor KRAS mutations with a low mutant allele frequency. While the percent of KRAS-positive tumors increases with the increase of the assay sensitivity, it can be explained by artifacts from sample fixation, e.g., formalin treatment for FFPE (formalin-fixed, paraffin-embedded) specimens. Therefore, the high sensitivity of the assay (<1%) itself has its pitfalls and may lead to false-positive results, and such samples with a low mutant allele frequency should be a subject of a more scrupulous analysis with additional procedures, namely, UDG (uracil-DNA glycosylase) treatment [7]. Nevertheless, by adapting the primer/blocker ratio, it is possible to improve the assay in each case to yield optimal Cq and ΔCq values. Using a higher DNA amount, as in previous studies [42,43], it is possible to obtain over 0.1% sensitivity.

4. Conclusions

PG-modification is a noncharged organic residue at the internucleotide phosphodiester fragment which does not significantly alter the secondary structure of oligonucleotide complexes with DNA [26]. However, the presence of the large organic residue (18 N,C,H-atoms instead of a single oxygen atom) in the tight DNA/enzyme interaction site can significantly influence the efficiency and selectivity of the enzymatic reaction. Increased selectivity of DNA-polymerases has been demonstrated during the use of primers modified at the internucleoside phosphodiester residue [16,17,34]. In the current work, PG-modified primers for AS-PCR and an allele-specific competitive blocker PCR (ASB-PCR) for KRAS mutations (G12V, G12A, G12D and G13D) detection, was presented.
Our results suggested that the PG-modified primers and blockers are advanced tools for AS-PCR or ASB-PCR methods. The incorporation of PG modification into PCR primers leads to the identification of single mismatches by AS-PCR and ASB-PCR assays with a good specificity without significant losses in the sensitivity in the abundant wild-type DNA background. Moreover, PG can serve as the universal additional mismatch-like disturbance to increase primer specificity and would greatly facilitate medical research. The method is straightforward for assay optimization and has been evaluated on plasmid model systems, providing a KRAS mutation detection limit of 20 copies of mutant DNA in a proportion as low as 0.1% of the total DNA, without nonspecific WT DNA amplification. A limitation of the proposed KRAS assay is that PG-modification, in some cases, leads to an increase in the Cq value, and partly inhibits PCR that leads to a decrease in PCR efficiency. In the future, it is hoped to use the properties of PGOs to establish optimal reproducible PCR technology in each laboratory, which may stimulate further research and assist clinicians.

Supplementary Materials

The following are available online at https://www.mdpi.com/2075-4418/10/11/872/s1, Figure S1: Cq values of KRAS G12A mutation detection using C*TGC, *CTGC, *C*TGC, CTTC, C*TTC, k-ref primers and various amount of mutant DNA per reaction. Standard curves showing the PCR efficiency values for the several primers generated from a dilution series of mutant DNA; Figure S2: Cq values of KRAS G12V, G12D mutations detection using CTAT, C*TAT, *CTAT, CGGA, C*GGA, *CGGA primers and various amount of mutant DNA per reaction. Standard curves showing the PCR efficiency values for the several primers generated from a dilution series of mutant DNA; Figure S3: Cq values of KRAS G13D mutations detection using GAGA, G*AGA, *GAGA primers and various amount of mutant DNA per reaction. Standard curves showing the PCR efficiency values for the several primers generated from a dilution series of mutant DNA; Table S1: AS-PCR KRAS G12A mutation detection using CTGC and CTTC primers with PG-modification in the different phosphates from the 3′-end and mutant DNA mixed with wild-type DNA (total 2 × 105 copies per reaction) at various ratios (% of mutant DNA); Table S2: AS-PCR reproducibility G12A KRAS mutation detection using wild-type DNA (total 2 × 104 copies per reaction) and 1% mutant DNA in the background of wild-type DNA; Table S3: ASB-PCR KRAS mutations detection of various mutant/wild-type DNA ratio (total 2 × 104 or 2 × 105 copies per reaction).

Author Contributions

Conceptualization, A.S.C., M.L.F. and D.V.P.; methodology, A.S.C. and I.P.O.; investigation, A.S.C. and I.P.O.; resources, I.P.O.; data curation, A.S.C. and I.P.O.; writing—original draft preparation, A.S.C.; writing—review and editing, I.P.O., M.L.F., A.A.L. and D.V.P.; supervision, D.V.P.; project administration, A.A.L.; funding acquisition, D.V.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Russian Science Foundation, grant number 18-14-00357 and the Russian State Federal budget project of ICBFM SB RAS in part of original theoretical conceptualization of the research, grant number АААА-А17-117020210021-7.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Anderson, S.M. Laboratory methods for KRAS mutation analysis. Expert Rev. Mol. Diagn. 2011, 11, 635–642. [Google Scholar] [CrossRef] [PubMed]
  2. Ishige, T.; Itoga, S.; Matsushita, K. Locked Nucleic Acid Technology for Highly Sensitive Detection of Somatic Mutations in Cancer; Elsevier BV: New York, NY, USA, 2018; Volume 83, pp. 53–72. Available online: https://www.sciencedirect.com/science/article/pii/S0065242317300616?via%3Dihub (accessed on 24 October 2020).
  3. Serebriiskii, I.G.; Connelly, C.; Frampton, G.; Newberg, J.; Cooke, M.; Miller, V.; Ali, S.; Ross, J.S.; Handorf, E.; Arora, S.; et al. Comprehensive characterization of RAS mutations in colon and rectal cancers in old and young patients. Nat. Commun. 2019, 10, 1–12. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Nagakubo, Y.; Hirotsu, Y.; Amemiya, K.; Oyama, T.; Mochizuki, H.; Omata, M. Accurate detection of KRAS, NRAS and BRAF mutations in metastatic colorectal cancers by bridged nucleic acid-clamp real-time PCR. BMC Med. Genom. 2019, 12, 162. [Google Scholar] [CrossRef] [Green Version]
  5. Jackson, M.P.; Truong, D.; Brownlow, M.L.; Wagner, J.A.; McKinley, R.A.; Bikson, M.; Jankord, R. Safety parameter considerations of anodal transcranial Direct Current Stimulation in rats. Brain Behav. Immun. 2017, 64, 152–161. [Google Scholar] [CrossRef] [PubMed]
  6. Lin, M.-T.; Mosier, S.L.; Thiess, M.; Beierl, K.F.; Debeljak, M.; Tseng, L.-H.; Chen, G.; Yegnasubramanian, S.; Ho, H.; Cope, L.; et al. Clinical Validation of KRAS, BRAF, and EGFR Mutation Detection Using Next-Generation Sequencing. Am. J. Clin. Pathol. 2014, 141, 856–866. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Tímár, J.; Kashofer, K. Molecular epidemiology and diagnostics of KRAS mutations in human cancer. Cancer Metastasis Rev. 2020, 1–10. [Google Scholar] [CrossRef]
  8. OPATHY Consortium; Arastehfar, A.; Boekhout, T.; Butler, G.; De Cesare, G.B.; Dolk, E.; Gabaldóln, T.; Hafez, A.I.; Hube, B.; Hagen, F.; et al. Recent trends in molecular diagnostics of yeast infections: From PCR to NGS. FEMS Microbiol. Rev. 2019, 43, 517–547. [Google Scholar] [CrossRef] [Green Version]
  9. Van Dijk, E.L.; Jaszczyszyn, Y.; Naquin, D.; Thermes, C. The Third Revolution in Sequencing Technology. Trends Genet. 2018, 34, 666–681. [Google Scholar] [CrossRef]
  10. Liu, Y.-P.; Wu, H.-Y.; Yang, X.; Xu, H.-Q.; Chen, D.; Huang, Q.; Fu, W.-L. Diagnostic accuracy of high resolution melting analysis for detection of KRAS mutations: A systematic review and meta-analysis. Sci. Rep. 2014, 4, 7521. [Google Scholar] [CrossRef] [Green Version]
  11. Rejali, N.A.; Moric, E.; Wittwer, C.T. The Effect of Single Mismatches on Primer Extension. Clin. Chem. 2018, 64, 801–809. [Google Scholar] [CrossRef]
  12. Morlan, J.; Baker, J.; Sinicropi, D. Mutation Detection by Real-Time PCR: A Simple, Robust and Highly Selective Method. PLoS ONE 2009, 4, e4584. [Google Scholar] [CrossRef]
  13. Wang, H.; Jiang, J.; Mostert, B.; Sieuwerts, A.M.; Martens, J.W.; Sleijfer, S.; Foekens, J.; Wang, Y. Allele-Specific, Non-Extendable Primer Blocker PCR (AS-NEPB-PCR) for DNA Mutation Detection in Cancer. J. Mol. Diagn. 2013, 15, 62–69. [Google Scholar] [CrossRef] [PubMed]
  14. McKinzie, P.B.; Parsons, B.L. Detection of rare K-ras codon 12 mutations using allele-specific competitive blocker PCR. Mutat. Res. Mol. Mech. Mutagen. 2002, 517, 209–220. [Google Scholar] [CrossRef]
  15. Navarro, E.; Serrano-Heras, G.; Castaño, M.; Solera, J. Real-time PCR detection chemistry. Clin. Chim. Acta 2015, 439, 231–250. [Google Scholar] [CrossRef] [PubMed]
  16. Li, T.-L.; Wu, M.-W.; Lin, W.-C.; Lai, C.-H.; Chang, Y.-H.; Su, L.-J.; Chen, W.-Y. Designed phosphate-methylated oligonucleotides as PCR primers for SNP discrimination. Anal. Bioanal. Chem. 2019, 411, 3871–3880. [Google Scholar] [CrossRef]
  17. Yang, H.-L.; Jiang, H.-J.; Fang, W.-Y.; Xu, Y.-Y.; Li, K.; Zhang, J.; Liao, D.-F.; He, F.-C. High fidelity PCR with an off/on switch mediated by proofreading polymerases combining with phosphorothioate-modified primer. Biochem. Biophys. Res. Commun. 2005, 328, 265–272. [Google Scholar] [CrossRef]
  18. D’Agata, R.; Giuffrida, M.C.; Spoto, G. Peptide Nucleic Acid-Based Biosensors for Cancer Diagnosis. Molecules 2017, 22, 1951. [Google Scholar] [CrossRef] [Green Version]
  19. Gupta, A.; Mishra, A.; Puri, N. Peptide nucleic acids: Advanced tools for biomedical applications. J. Biotechnol. 2017, 259, 148–159. [Google Scholar] [CrossRef]
  20. Saarbach, J.; Sabale, P.M.; Winssinger, N. Peptide nucleic acid (PNA) and its applications in chemical biology, diagnostics, and therapeutics. Curr. Opin. Chem. Biol. 2019, 52, 112–124. [Google Scholar] [CrossRef]
  21. Ballantyne, K.N.; Van Oorschot, R.; Mitchell, R. Locked nucleic acids in PCR primers increase sensitivity and performance. Genomics 2008, 91, 301–305. [Google Scholar] [CrossRef] [Green Version]
  22. Kuo, T.-C.; Wu, M.-W.; Lin, W.-C.; Matulis, D.; Yang, Y.-S.; Li, S.-Y.; Chen, W.-Y. Reduction of interstrand charge repulsion of DNA duplexes by salts and by neutral phosphotriesters—Contrary effects for harnessing duplex formation. J. Taiwan Inst. Chem. Eng. 2020, 110, 1–7. [Google Scholar] [CrossRef]
  23. Kupryushkin, M.S.; Pyshnyi, D.V.; Stetsenko, D.A. Phosphoryl Guanidines: A New Type of Nucleic Acid Analogues. Acta Nat. 2014, 6, 116–118. [Google Scholar] [CrossRef]
  24. Stetsenko, D.A.; Kupryushkin, M.S.; Pyshnyi, D. Modified Oligonucleotides and Methods for Their Synthesis. WO2016028187A1, 25 February 2016. [Google Scholar]
  25. Dyudeeva, E.S.; Kupryushkin, M.S.; Lomzov, A.A.; Pyshnaya, I.A.; Pyshnyi, D.V. Physicochemical Properties of the Phosphoryl Guanidine Oligodeoxyribonucleotide Analogs. Russ. J. Bioorg. Chem. 2019, 45, 709–718. [Google Scholar] [CrossRef]
  26. Lomzov, A.A.; Kupryushkin, M.S.; Shernyukov, A.V.; Nekrasov, M.D.; Dovydenko, I.S.; Stetsenko, D.A.; Pyshnyi, D.V. Diastereomers of a mono-substituted phosphoryl guanidine trideoxyribonucleotide: Isolation and properties. Biochem. Biophys. Res. Commun. 2019, 513, 807–811. [Google Scholar] [CrossRef] [PubMed]
  27. Dmitrienko, E.V.; Naumova, O.V.; Fomin, B.; Kupryushkin, M.; Volkova, A.; Amirkhanov, N.; Semenov, D.; Pyshnaya, I.; Pyshnyi, D.V. Surface modification of SOI-FET sensors for label-free and specific detection of short RNA analyte. Nanomedicine 2016, 11, 2073–2082. [Google Scholar] [CrossRef] [PubMed]
  28. Garafutdinov, R.R.; Sakhabutdinova, A.R.; Kupryushkin, M.S.; Pyshnyi, D.V. Prevention of DNA multimerization using phosphoryl guanidine primers during isothermal amplification with Bst exo-DNA polymerase. Biochimie 2020, 168, 259–267. [Google Scholar] [CrossRef] [PubMed]
  29. Epanchintseva, A.; Dolodoev, A.; Grigor’Eva, A.; Chelobanov, B.; Pyshnyi, D.V.; Ryabchikova, E.I.; Pyshnaya, I. Non-covalent binding of nucleic acids with gold nanoparticles provides their stability and effective desorption in environment mimicking biological media. Nanotechnology 2018, 29, 355601. [Google Scholar] [CrossRef]
  30. Kuznetsov, N.A.; Kupryushkin, M.S.; Abramova, T.; Kuznetsova, A.A.; Miroshnikova, A.D.; Stetsenko, D.A.; Pyshnyi, D.V.; Fedorova, O.S. New oligonucleotide derivatives as unreactive substrate analogues and potential inhibitors of human apurinic/apyrimidinic endonuclease APE1. Mol. BioSyst. 2016, 12, 67–75. [Google Scholar] [CrossRef] [Green Version]
  31. Lebedeva, N.A.; Anarbaev, R.O.; Kupryushkin, M.S.; Rechkunova, N.I.; Pyshnyi, D.V.; Stetsenko, D.A.; Lavrik, O.I. Design of a New Fluorescent Oligonucleotide-Based Assay for a Highly Specific Real-Time Detection of Apurinic/Apyrimidinic Site Cleavage by Tyrosyl-DNA Phosphodiesterase 1. Bioconjug. Chem. 2015, 26, 2046–2053. [Google Scholar] [CrossRef]
  32. Kupryushkin, M.S.; Pyshnaya, I.A.; Dmitrienko, E.V.; Stetsenko, D.A.; Filipenko, M.L.; Oscorbin, I.P.; Stepanov, G.A.; Richter, V.A.; Ivanov, M.K.; Pyshnyi, D.V. Template-Directed Enzymatic DNA Synthesis Using Phosphoryl Guanidine Oligonucleotides. 2019. Available online: https://patentscope.wipo.int/search/en/detail.jsf?docId=WO2019112485 (accessed on 24 October 2020).
  33. Kupryushkin, M.S.; Pyshnaya, I.A.; Dmitrienko, E.V.; Stetsenko, D.A.; Filipenko, M.L.; Oscorbin, I.P.; Stepanov, G.A.; Richter, V.A.; Ivanov, M.K.; Pyshnyi, D.V. The Method of Amplification of Nucleic Acids Using Phosphoryl Guanidine Oligonucleotides. 2017. Available online: https://patentscope.wipo.int/search/en/detail.jsf?docId=WO2019112485 (accessed on 24 October 2020).
  34. Vinogradova, O.A.; Pyshnyi, D.V. Selectivity of Enzymatic Conversion of Oligonucleotide Probes during Nucleotide Polymorphism Analysis of DN. Acta Nat. 2010, 2, 36–52. [Google Scholar] [CrossRef]
  35. Lang, A.H.; Drexel, H.; Geller-Rhomberg, S.; Stark, N.; Winder, T.; Geiger, K.; Muendlein, A. Optimized Allele-Specific Real-Time PCR Assays for the Detection of Common Mutations in KRAS and BRAF. J. Mol. Diagn. 2011, 13, 23–28. [Google Scholar] [CrossRef] [PubMed]
  36. Liu, J.; Huang, S.; Sun, M.; Liu, S.; Liu, Y.; Wang, W.; Zhang, X.; Wang, H.; Hua, W. An improved allele-specific PCR primer design method for SNP marker analysis and its application. Plant Methods 2012, 8, 34. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Takeda, S.; Ichii, S.; Nakamura, Y. Detection of K-ras mutation in sputum by mutant-allele-specific amplification (MASA). Hum. Mutat. 1993, 2, 112–117. [Google Scholar] [CrossRef] [PubMed]
  38. Olmedillas-López, S.; García-Olmo, D.C.; Garcia-Arranz, M.; Guadalajara, H.; Pastor, C.; García-Olmo, D. KRAS G12V Mutation Detection by Droplet Digital PCR in Circulating Cell-Free DNA of Colorectal Cancer Patients. Int. J. Mol. Sci. 2016, 17, 484. [Google Scholar] [CrossRef] [Green Version]
  39. Vestheim, H.; Deagle, B.E.; Jarman, S.N. Application of Blocking Oligonucleotides to Improve Signal-to-Noise Ratio in a PCR. Methods Mol. Biol. 2011, 687, 265–274. [Google Scholar]
  40. Tolnai, Z.; Harkai, Á.; Szeitner, Z.; Scholz, É.N.; Percze, K.; Gyurkovics, A.; Mészáros, T. A simple modification increases specificity and efficiency of asymmetric PCR. Anal. Chim. Acta 2019, 1047, 225–230. [Google Scholar] [CrossRef]
  41. Kim, H.E.; Ruby, A.; Shandilya, H.G.; Virmani, A.K.; Rahman, N.; Strange, C.M.; Huuskonen, J. T-blocker: A simple and robust probe-free quantitative PCR assay to detect somatic mutations down to 0.1% frequency. Biotechniques 2018, 65, 205–210. [Google Scholar] [CrossRef] [Green Version]
  42. Huang, J.-F.; Zeng, D.-Z.; Duan, G.-J.; Shi, Y.; Deng, G.-H.; Xia, H.; Xu, H.-Q.; Zhao, N.; Fu, W.-L.; Huang, Q. Single-Tubed Wild-Type Blocking Quantitative PCR Detection Assay for the Sensitive Detection of Codon 12 and 13 KRAS Mutations. PLoS ONE 2015, 10, e0145698. [Google Scholar] [CrossRef]
  43. Huang, Q.; Wang, G.-Y.; Huang, J.; Zhang, B.; Fu, W.-L. High sensitive mutation analysis on KRAS gene using LNA/DNA chimeras as PCR amplification blockers of wild-type alleles. Mol. Cell. Probes 2010, 24, 376–380. [Google Scholar] [CrossRef]
  44. Albitar, A.Z.; Ma, W.; Albitar, M. Wild-type Blocking PCR Combined with Direct Sequencing as a Highly Sensitive Method for Detection of Low-Frequency Somatic Mutations. J. Vis. Exp. 2017, e55130. [Google Scholar] [CrossRef]
Figure 1. Taq-polymerase and template/primer complex schematic representation with the primer 3′-end nucleotide numeration. Sphere indicates phosphoryl guanidine (PG) modification of internucleoside phosphate moiety.
Figure 1. Taq-polymerase and template/primer complex schematic representation with the primer 3′-end nucleotide numeration. Sphere indicates phosphoryl guanidine (PG) modification of internucleoside phosphate moiety.
Diagnostics 10 00872 g001
Figure 2. Typical amplification curves obtained in allele-specific PCR assays. KRAS G12A mutation detection was performed using CTTC and C*TTC primers and WT DNA (2 × 104 copies per reaction), 1% and 0.1% of mutant DNA in the background of WT DNA.
Figure 2. Typical amplification curves obtained in allele-specific PCR assays. KRAS G12A mutation detection was performed using CTTC and C*TTC primers and WT DNA (2 × 104 copies per reaction), 1% and 0.1% of mutant DNA in the background of WT DNA.
Diagnostics 10 00872 g002
Figure 3. A schematic representation of wild-type blocking PCR. Wild-type variant is marked by blue, mutant—by red. AS-forward, reverse primer and one WT-blocking primer constitute the PCR reaction mixture. On the 3′-terminus of the blocking primer, inhibiting three PG groups have been attached to prevent extension by DNA polymerization. The AS-primer overlapped with the blocking primer at the target mutation site and was in competition with the blocking primer.
Figure 3. A schematic representation of wild-type blocking PCR. Wild-type variant is marked by blue, mutant—by red. AS-forward, reverse primer and one WT-blocking primer constitute the PCR reaction mixture. On the 3′-terminus of the blocking primer, inhibiting three PG groups have been attached to prevent extension by DNA polymerization. The AS-primer overlapped with the blocking primer at the target mutation site and was in competition with the blocking primer.
Diagnostics 10 00872 g003
Figure 4. Typical amplification curves obtained in allele-specific blocking PCR assays. KRAS G12A mutation detection was performed using CTTC, C*TTC and C*TGC primers and C*T*G*G blocker and WT DNA (2 × 104 copies per reaction), 1%, 0.1% and 0.01% mutant DNA in the background of WT DNA.
Figure 4. Typical amplification curves obtained in allele-specific blocking PCR assays. KRAS G12A mutation detection was performed using CTTC, C*TTC and C*TGC primers and C*T*G*G blocker and WT DNA (2 × 104 copies per reaction), 1%, 0.1% and 0.01% mutant DNA in the background of WT DNA.
Diagnostics 10 00872 g004
Table 1. Allele-specific real-time polymerase chain reaction (AS-PCR) KRAS G12A mutation detection using mutant DNA mixed with wild type (WT) DNA (total 2 × 105 copies per reaction) at various ratios (% of mutant DNA).
Table 1. Allele-specific real-time polymerase chain reaction (AS-PCR) KRAS G12A mutation detection using mutant DNA mixed with wild type (WT) DNA (total 2 × 105 copies per reaction) at various ratios (% of mutant DNA).
PrimersCq
WT1%100%
CTTC34.25 ± 0.0528.7 ± 0.121.72 ± 0.02
C*TTCN/A29.65 ± 0.0723.8 ± 0.2
CT*TCN/A33.24 ± 0.0126.8 ± 0.2
CTT*CN/A33.5 ± 0.126.67 ± 0.01
C*TT*CN/A36.4 ± 0.229.2 ± 0.1
*C*TTCN/A37.9 ± 0.132.0 ± 0.1
C*T*TCN/A40.9 ± 0.234.8 ± 0.1
*CTT*CN/AN/A40.1 ± 0.3
CT*T*CN/AN/AN/A
No template control (NTC) was undetermined in all the reactions, N/A indicates that no Cq was obtained for a typical 45-cycle reaction. Symbol “*” means phosphoryl guanidine (PG) modification location. Boldly marked nucleotides mean mismatched nucleotides in relation to the WT DNA sequence.
Table 2. AS-PCR KRAS mutation detection using WT DNA (total 2 × 104 copies per reaction) and 1% mutant DNA in the background of WT DNA.
Table 2. AS-PCR KRAS mutation detection using WT DNA (total 2 × 104 copies per reaction) and 1% mutant DNA in the background of WT DNA.
PrimersCq
WT1%
G12VCTAT33.4 ± 0.131.9 ± 0.1
C*TAT38.7 ± 0.333.2 ± 0.1
*CTATN/A40.2 ± 1.3
G12DCGGA35.8 ± 0.132.50 ± 0.05
C*GGA41.7 ± 1.832.9 ± 0.1
*CGGAN/A42.9 ± 1.2
G12ACTTC37.7 ± 0.332.0 ± 0.2
С*TTCN/A32.7 ± 0.2
G13DGAGA36.3 ± 0.132.3 ± 0.1
G*AGAN/A35.0 ± 0.1
*GAGAN/A39.2 ± 0.3
No template control (NTC) was undetermined in all the reactions, N/A indicates that no Cq was obtained for a typical 45-cycle reaction. Symbol “*” means PG modification location. Boldly marked nucleotides represent mismatched nucleotides in relation to the WT DNA sequence.
Table 3. AS-PCR KRAS mutation detection using WT DNA (total 2 × 104 copies per reaction) and 1% mutant DNA in the background of WT DNA.
Table 3. AS-PCR KRAS mutation detection using WT DNA (total 2 × 104 copies per reaction) and 1% mutant DNA in the background of WT DNA.
PrimersCq ΔCqPrimersCq ΔCq
WT1%CqWT − Cq1%WT1%CqWT − Cq1%
G12ACTGC33.0 ± 0.531.5 ± 0.11.5G12VCTGT26.2 ± 0.126.0 ± 0.10.2
С*TGC38.4 ± 0.432.6 ± 0.15.8C*TGT30.7 ± 0.129.7 ± 0.11.0
CTTC37.7 ± 0.332.0 ± 0.25.7CTAT33.4 ± 0.131.9 ± 0.11.5
С*TTCN/A32.7 ± 0.212.3 aC*TAT38.7 ± 0.333.2 ± 0.15.5
CAGC37.7 ± 0.233.0 ± 0.14.7CTTT35.3 ± 0.134.2 ± 0.11.1
C*AGC38.7 ± 0.332.5 ± 0.16.2C*TTT41.8 ± 1.137.9 ± 0.33.9
G13DGTGA29.0 ± 0.128.1 ± 0.10.9CGGT28.0 ± 0.127.4 ± 0.10.6
G*TGA33.4 ± 0.132.3 ± 0.11.1C*GGT28.5 ± 0.127.6 ± 0.10.9
GAGA36.3 ± 0.232.3 ± 0.14.0CCGT31.7 ± 0.130.0 ± 0.11.7
G*AGAN/A35.0 ± 0.110.0 aC*CGT33.9 ± 0.230.9 ± 0.13.0
GCGA35.0 ± 0.131.0 ± 0.14.0G12DCTGA27.7 ± 0.127.0 ± 0.10.7
G*CGAN/A33.1 ± 0.111.9 aC*TGA29.2 ± 0.128.7 ± 0.10.5
GGGA38.0 ± 0.432.8 ± 0.15.2CGGA35.8 ± 0.132.50 ± 0.053.3
G*GGAN/A36.4 ± 0.28.6 aC*GGA41.7 ± 1.832.9 ± 0.18.8
No template control (NTC) was undetermined in all the reactions, N/A indicates that no Cq was obtained for a typical 45-cycle reaction. Symbol “*” means PG modification location. Boldly marked nucleotides represent mismatched nucleotides in relation to the WT DNA sequence. SNP G/C (G12A), G/T (G12V), G/A (G12D, G13D); 3′-end mismatch site with mutated DNA C/C, Pyr/Pyr (G12A); C/T, Pyr/Pyr (G12V); C/A, Pyr/Pur (G12D, G13D). a If Cq(WT) = N/A to calculate ΔCq = Cq(WT) − Cq(1%) value Cq(WT) = 45 has been used.
Table 4. KRAS mutations detection by AS-PCR using various mutant/wild-type DNA ratios (total 2 × 104 copies per reaction) and the several best-obtained primers for each mutation detection.
Table 4. KRAS mutations detection by AS-PCR using various mutant/wild-type DNA ratios (total 2 × 104 copies per reaction) and the several best-obtained primers for each mutation detection.
PrimersCqΔCq
WT0.1%1%CqWT − Cq0.1%CqWT − Cq1%
G12AСTTC37.7 ± 0.335.1 ± 0.132.0 ± 0.22.65.7
С*TTCN/A36.9 ± 0.232.7 ± 0.28.1 a12.3 a
G12VCTAT33.4 ± 0.133.2 ± 0.131.9 ± 0.10.21.5
C*TAT38.7 ± 0.437.1 ± 0.233.2 ± 0.11.65.5
G12DCGGA35.8 ± 0.135.1 ± 0.132.5 ± 0.10.73.3
C*GGA41.7 ± 1.836.2 ± 0.232.9 ± 0.15.59.2
G13DGAGA36.3 ± 0.134.6 ± 0.132.3 ± 0.11.74.0
G*AGAN/A38.5 ± 0.335.0 ± 0.16.5 a10.0 a
GCGA35.0 ± 0.134.3 ± 0.131.0 ± 0.10.74.0
G*CGAN/A37.0 ± 0.333.1 ± 0.18.0 a11.9 a
No template control (NTC) was undetermined in all the reactions, N/A indicates that no Cq was obtained for a typical 45-cycle reaction. Symbol “*” means PG modification location. Boldly marked nucleotides represent mismatched nucleotides in relation to the WT DNA sequence. a If Cq(WT) = N/A to calculate ΔCq = Cq(WT) − Cq(1% or 0.1%) value Cq(WT) = 45 was used.
Table 5. ASB-PCR KRAS mutations detection of various mutant/wild-type DNA ratios (total 2 × 104 or 2 × 105 copies per reaction).
Table 5. ASB-PCR KRAS mutations detection of various mutant/wild-type DNA ratios (total 2 × 104 or 2 × 105 copies per reaction).
PrimersCqΔCq
CqWT − Cq1%
2 × 104 Copies2 × 104 2 × 105
1%0.1%0.01%WTCopies
G12ANo blocker CTGG24.6 ± 0.224.8 ± 0.225.0 ± 0.325.0 ± 0.20.4-
Blocker C*T*G*GN/AN/AN/AN/A--
No blocker CTTC31.5 ± 0.133.9 ± 0.134.2 ± 0.236.7 ± 0.25.25.7
CTTC/C*T*G*G = 1/0.133.5 ± 0.138.0 ± 0.3N/AN/A11.5 a10.7
C*TTC/C*T*G*G = 1/0.2532.5 ± 0.136.3 ± 0.338.2 ± 0.4N/A12.5 a14.2
C*TGC/C*T*G*G = 1/0.2533.1 ± 0.137.5 ± 0.240.1 ± 0.542.4 ± 1.59.39.2
G12VNo blocker CTAT30.6 ± 0.130.8 ± 0.231.4 ± 0.231.6 ± 0.21.0-
C*TAT /C*T*G*G = 1/0.136.0 ± 0.237.5 ± 0.440.3 ± 1.0N/A9.0 a8.8
G12DNo blocker CGGA33.2 ± 0.136.0 ± 0.337.0 ± 0.537.13.8-
C*GGA/C*T*G*G = 1/0.133.6 ± 0.137.1 ± 0.639.2 ± 0.9N/A11.4 a9.8
G13DNo blocker GTGG23.5 ± 0.124.0 ± 0.224.1 ± 0.124.2 ± 0.30.7-
Blocker G*T*G*GN/AN/AN/AN/A--
No blocker GAGA31.1 ± 0.133.1 ± 0.233.6 ± 0.335.2 ± 0.24.1-
GAGA/G*T*G*G = 1/0.132.7 ± 0.134.8 ± 0.2N/AN/A12.3 a8.8
G*AGA/G*T*G*G = 1/0.135.3 ± 0.439.2 ± 0.841.0 ± 1.3N/A9.7 a10.9
No template control (NTC) was undetermined in all the reactions, N/A indicates that no Cq was obtained for a typical 45-cycle reaction. Symbol “*” means PG modification location. Boldly marked nucleotides represent mismatched nucleotides in relation to the WT DNA sequence. Real-time PCR assay was done using a constant 450 nM AS-primers concentration and several blocker primers ratios. a If Cq(WT) = N/A to calculate ΔCq = Cq(WT) − Cq(1%) value Cq(WT) = 45 was used.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chubarov, A.S.; Oscorbin, I.P.; Filipenko, M.L.; Lomzov, A.A.; Pyshnyi, D.V. Allele-Specific PCR for KRAS Mutation Detection Using Phosphoryl Guanidine Modified Primers. Diagnostics 2020, 10, 872. https://doi.org/10.3390/diagnostics10110872

AMA Style

Chubarov AS, Oscorbin IP, Filipenko ML, Lomzov AA, Pyshnyi DV. Allele-Specific PCR for KRAS Mutation Detection Using Phosphoryl Guanidine Modified Primers. Diagnostics. 2020; 10(11):872. https://doi.org/10.3390/diagnostics10110872

Chicago/Turabian Style

Chubarov, Alexey S., Igor P. Oscorbin, Maxim L. Filipenko, Alexander A. Lomzov, and Dmitrii V. Pyshnyi. 2020. "Allele-Specific PCR for KRAS Mutation Detection Using Phosphoryl Guanidine Modified Primers" Diagnostics 10, no. 11: 872. https://doi.org/10.3390/diagnostics10110872

APA Style

Chubarov, A. S., Oscorbin, I. P., Filipenko, M. L., Lomzov, A. A., & Pyshnyi, D. V. (2020). Allele-Specific PCR for KRAS Mutation Detection Using Phosphoryl Guanidine Modified Primers. Diagnostics, 10(11), 872. https://doi.org/10.3390/diagnostics10110872

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop