http://orcid.org/0000-0002-1075-3773
Figures
Abstract
Single-nucleotide polymorphisms (SNPs) is associated with efficacy of specific drugs. Although there are several methods for SNP genotyping in clinical settings, alternative methods with lower cost, higher throughout and less complexity are still needed. In this study, we modified Kompetitive Allele Specific PCR to enable multiplex SNP genotyping by introducing additional fluorescent cassettes that specifically help to differentiate more amplification signals in a single reaction. This new format of assay achieved a limit of detection down to 310 copies/ reactions for simultaneous detection of 2 SNPs with only standard end-point PCR workflow for synthetic controls, and genotyped 117 clinical samples with results that were in 100% agreement with hospital reports. This study presented a simplified, cost-effective high-throughput SNP genotyping alternative for pharmacogenetic variants, and enabled easier access to pharmaceutical guidance when needed.
Citation: Suo W, Shi X, Xu S, Li X, Lin Y (2020) Towards low cost, multiplex clinical genotyping: 4-fluorescent Kompetitive Allele-Specific PCR and its application on pharmacogenetics. PLoS ONE 15(3): e0230445. https://doi.org/10.1371/journal.pone.0230445
Editor: Ruslan Kalendar, University of Helsinki, FINLAND
Received: December 2, 2019; Accepted: February 29, 2020; Published: March 16, 2020
Copyright: © 2020 Suo et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All relevant data are within the paper.
Funding: The experiment was supported from Special fund support project for clinical medicine development of Beijing Municipal Hospital Administration(ZYLX201805) and Major new drug creation national major science and technology special project "construction of clinical research demonstration platform for cardiovascular disease drugs"(2017ZX09304017). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Background
Single-nucleotide polymorphisms (SNPs) is associated with efficacy of specific drugs. For pharmacogenomic purpose SNP test is extensively used to guide clinical decisions. For example, Polymorphisms of the CYP2C9 (cytochrome P450, family 2, subfamily C, polypeptide 9) gene (CYP2C9*2, CYP2C9*3) and the VKORC1 (vitamin K epoxide reductase complex, subunit 1) gene (−1639G>A) can pose major impact on the maintenance dose of certain drugs such as Warfarin. Genotyping target SNPs for patients provides a time-saving method to guide proper maintenance dose, while reducing the adverse reaction risks and reoccurrence of thromboembolic episodes. [1] The genotyping of glucose-6-phosphate dehydrogenase (G6PD) is another typical example which demonstrates how population genotyping provides guidance in personalized medicine. G6PD deficiency prevalence is found to be high in malaria-endemic regions, while dominant antimalarial drugs can cause hemolysis to various degrees in G6PD deficient patients. Screening of G6PD deficiency alleles in high-risk populations would provide valuable estimation of hemolysis risk when using drugs such as primaquine, thereby improving public drug safety.[2]
Many methods exist for clinical SNP genotyping. Majority of methods that can be used in clinical settings are PCR-based [3], which include, probe based qPCR assay that detect SNP by allele specific probes labeled with difference fluorescent dye [4]; allele-specific PCR (ASPCR) using SNP-containing primers only direct amplification on their complementary allele [5]; melting curve analysis (high-resolution melting analysis) utilizing saturated intercalating dye that would display the change of melting curve caused by slight difference of Tm between two SNP alleles [6]; and Cold-PCR which amplifies minority alleles selectively from mixtures of wild-type, utilizing the difference of dsDNA-dissociation by precisely controlled denaturing temperature [7]. These chemistries represent the most commercial kit now used for clinical diagnosis of SNP. Other methods also exist, such as Mass Array [2], Capillary electrophoresis [8] and digital PCR [9]. Although being high-throughput or super sensitive, these methods highly rely on the post-PCR analysis, which is time-consuming, vulnerable to cross contamination and equipment-demanding, therefore hindering the application in clinical settings where turnaround time, robustness and simplicity is crucial.
Improved from ASPCR, Kompetitive Allele Specific PCR (KASP, LGC Group) introduced fluorescence resonance energy transfer (FRET) for signal generation, where 2 fluorescent cassettes are used for identification of allele-specific amplification for a single bi-allelic SNP. Compared with other alternatives, KASP is more flexible [10], cost-effective [11] and robust [12], and therefore is now widely applied in large scale genotyping for many areas including agriculture [13] and human health [14, 15]. However, due to KASP’s closed platform, the multiplexity is limited to only single bi-allelic SNP in each reaction, restricting the further improvement of the throughput and efficiency.
In Beijing Anzhen Hospital, we provided 22,593 patients with genotype test and pharmaceutical guidance report service in 2018. Average number of SNPs tested for each patient is between 2~5. The current kit adopted for these tests utilize a ligase chain reaction-based technology that detect single SNP in each reaction, which means at least over 22,000 separate reactions are needed for the above tests. If 2 or even more SNPs can be tested in 1 reaction, more than 10, 000 reactions can be saved, enabling lower cost and reduced turnaround time. To address this problem, we aim to develop an enhanced version of KASP: 4-fluorescent KASP, which simultaneously utilize more fluorescent channels than just FAM/ Hex, and thereby enabling multiplexing genotyping. Here we describe the feasibility of genotyping with the new cassette with Cal Red 610 and Quasar 670 labels, and the compatibility of the new cassette with those contained in standard KASP reagent.
Method
Design of fluorescent cassette
Following FRET principle, one fluorescent dye and one BHQ-2 quencher were attached to 5’ and 3’ ends of two reverse complementary single stranded DNA separately. The sequence was blasted to confirm that it would not amplify any natural sequences. These sequences were sent to LGC, Biosearch Technologies for synthesis. Multiple pairs of cassettes were designed for empirical selection. The final selected sequences are as follows,
Design of KASP primers
The allele specific primer was designed by following the principles described previously[5], and validated with Primer-BLAST tool against NCBI database. The sequence of the used primers is available upon request. For standard KASP, sequences of specific primers were respectively added at 5’ end by FAM tag or Hex tag following previous published papers [16]. For modified KASP which involved fluorescent other than FAM/ Hex, and tails for specific fluorescent cassette were added, which were CAL Fluor Red 610 tag or Quasar 670 tag. Oligos were ordered through LGC, Biosearch Technologies. The sequences are as follows,
Design of positive control
Positive control for each target was designed and ordered through LGC, Biosearch Technologies. The sequence was designed by including primer binding site for allele specific sequences of each primer, with ‘tt’ as spacer. The sequences are as follows,
Clinical samples
This study was approved by the Clinical Research Ethics Committee of Beijing Anzhen Hospital. As the samples were collected for clinical tests following routine protocol and we used only small portion of the retention before it was discarded, the Clinical Research Ethics Committee of Beijing Anzhen Hospital agree to exempt the informed consent. The clinical DNA, extracted from blood samples between Oct. 12th to Jan. 6th 2020, was provided with clinically diagnosed genotype. The information was accompanied with sample ID only, personal information were not revealed for privacy purpose.
Clinical diagnostics.
Clinical diagnostic in our hospital was done using a commercial kit from Sino Era Genotech (http://www.sino-era.com/, Beijing, China). The kit includes reagents for sample preparation, molecular beacon probes and hybridization buffer. Briefly, 200 μL of whole blood was added to 1 mL of red blood cell lysis buffer and incubated for 5 minutes at room temperature. The mixture was centrifuged (3000 rpm) for 5 minutes and the supernatant removed. 50 μL preservation solution was added to resuspend the precipitate. 1.5 μL of sample was then combined with detecting reagent and analyzed with the TL998A real-time PCR system (Tianlong, Xi’an, China). The detecting reagent includes hybridization buffer and two molecular beacon probes labeled with FAM and HEX for wild-type and mutation, respectively.
PCR protocol and data analysis
For standard KASP, the reaction was done with KASP V4.0 2X Master Mix (Low Rox, LGC Biosearch Technologies) following manufacturer’s instructions. For assays with additional fluorescent channels, 0.08 μM of both fluorescent oligos and 0.32 μM of corresponding quencher oligos were added to BHQ probe mastermix (LGC Biosearch Technologies), or KASP V4.0 2X Master Mix. For each SNP, 0.15 μM of each allele specific primer and 0.4 μM of common reverse primer were used. The reaction volume in this study was 10 μl unless other wisely specified. The PCR protocol was the same as standard KASP (36 cycles of amplification in total), except that all fluorescent channel was read at the final step. The reaction was done with Biorad CFX96 and analyzed by Bio-Rad CFX Maestro 1.1 (Version:4.1.2433.1219) software.
Result
Evaluation of cassettes with Cal Red 610 and Quasar 670
Feasibility of doing self-made KASP with fluorescent channels other than FAM/ Hex was evaluated by adding newly designed fluorescent cassette labeled with Cal Red 610 and Quasar 670 to routinize PCR master mix (BHQ probe Mastermix, LGC Biosearch Technologies). Synthetic controls for rs1057910 and rs4986893 were assayed with this self-made KASP system, and all 16 controls were correctly clustered for both assays (Fig 1), indicating that this new pair of cassettes with Cal Red 610 and Quasar 670 works fine both for rs1057910 and rs4986893 in routine qPCR master mix.
Genotyping of rs1057910 (left) and rs4986893 using cassette with Cal Red 610 and Quasar 670 in routine qPCR master mix following KASP protocol. Allele specific primers were added with tail for Cal Red 610 (G for rs1057910, and C for rs4986893) or Quasar 670 (T for rs1057910, and T for rs4986893). Synthetic positive control of each gene was genotyped individually (homozygotes being approximately 300,000 copies/ reactions, and heterozygotes being 150,000 copies of each allele per reaction). One None-Template Control (NTC) was included in each run.
Evaluation of new cassette in standard KASP reagent
To confirm the compatibility of the new cassette with standard KASP chemistry, cassettes with Cal Red 610 and Quasar 670 were added to standard KASP, and genotyped clinical samples as well as synthetic controls for rs1057910 and rs4986893. The controls are all correctly genotyped, and the result for the clinical samples is in 100% agreement with that from hospital report (Table 1), indicating that the newly designed cassettes with Cal Red 610 and Quasar 670 work well within standard KASP reagent.
Limit of detection of 4-fluorescent KASP
The limit of detection (LOD) was determined as the minimal amount of positive control added to each reaction that produced a correct clustering. The LOD was determined to be 310 copies/ reactions, or 5 μl of 62 copies/ μl positive control (Fig 2), indicating that this new 4-fluorescent KASP was highly sensitive in all four channels, and could meet the needs of clinical detection.
Allele specific primers were added with tails for Cal Red 610 (G for rs1057910), Quasar 670 (T for rs1057910), FAM (A for rs9923231) and Hex (G for rs9923231). Positive control of rs1057910 and rs9923231 were diluted by 10-fold from 6.23 x 1010 copies/ μl (100 nM stock solution) to 6.2 copies/ μl. The two synthetic genes were diluted together to make duplex-SNP controls, where allele-G for rs1057910 and allele-A for rs9923231 were put in one solution and the other two alleles were put in another solution. The heterozygote controls were made by mixing homozygote solution of both alleles with 1:1 ratio at each concentration. 5 μl of each concentration was used in each reaction. For each concentration, each genotype was run in duplicates, with a total of 12 NTCs included. The data above showed the concentrations from 62 copies/ μl to 6.23 x 106 copies/ μl (6 concentrations in total). The 6.2 copies/ μl samples were not clustered correctly and therefore not included.
Discussion
Clinical SNP test is critical for prescribing specific drugs [1, 2]. In Beijing Anzhen Hospital, we provided 22,593 patients with genotype test and pharmaceutical guidance report service in 2018. The currently used genotyping assays were approved by the State Drug Administration of China, all of which involved more than just 1 SNP. These assays can only detect 1 SNP in each reaction, therefore multiple reactions are required for each patient. In addition, the cost is up to $ 60 per SNP, limiting the access of the low-income patients and thus exposing them to higher risks. Alternative method with higher throughout and lower cost is imperatively needed.
Although Mass-array and NGS offer super high multiplexibility, ranging from tens to even thousands of SNPs with reasonable price, they do not fit clinical needs at Anzhen hospital. As only 2~5 SNPs are required for each patient, and the result shall be returned within 24 hours, rendering both Mass-array and NGS incompatible in this scenario. To address this problem, we also considbuquered probe-based multiplex qPCR, which efficiently detects 2~3 SNPs in one reaction, but the perplexity of designing both primers and probes for each SNP, and the burden of fine tune PCR conditions for each assay is challenging.
With this new 4-fluorescent KASP assay available, the output for SNP genotyping will be doubled without requiring additional labor or significant cost. The workflow is exactly the same as standard KASP, which is already well known for its simplicity and flexibly [10, 12], and the additional cassette added to the assay is only single-labeled oligos, which is significantly less expensive as compared with dual-labeled oligos used for routine multiplex genotyping. The cost can be further reduced by making this single-labeled oligos in bulk; these single-labeled oligos are not gene-specific, and therefore applicable for genotyping of all SNPs. The universal-nature of these oligos bring in at least two benefit: 1, the unit cost (per nmol) of synthesizing oligos goes down significantly as the scale goes up. The universal oligos that can be made in one large batch for all SNPs will be more cost-effective than those individually customized, small scale probes; 2, the perplexity of designing both primers and probes for each SNP no longer exists, as only gene specific primers are needed, leading to improved convenience and flexibility for designing new assays for multiplexing SNP analysis.
The limitations of current study include limited scale of clinical validation, and small number of SNP sites evaluated for dual-SNP assays. To further support 4-fluorescent KASP’s capacity on clinical genotyping, experiments that involve more clinical samples, and various SNP sites should be conducted.
Current result, to our knowledge, is the first demonstration of performing KASP for multiple SNPs in one reaction. Although we only used 2 pairs of cassettes with 4-fluorescent channels for simultaneously genotyping of 2 SNPs in one reaction, there is the potential that more cassettes can be added and thus more SNPs and even multi-allelic SNPs can be evaluated in a single reaction, considering that 6 or more channels plate reader and qPCR machines are emerging.
References
- 1. Aomori T, Yamamoto K, Oguchi-Katayama A, Kawai Y, Ishidao T, Mitani Y, et al. Rapid single-nucleotide polymorphism detection of cytochrome P450 (CYP2C9) and vitamin K epoxide reductase (VKORC1) genes for the warfarin dose adjustment by the SMart-amplification process version 2. Clin Chem. 2009;55(4):804–12. Epub 2009/02/03. pmid:19181737.
- 2. Tian X, Zhou J, Zhao Y, Cheng Z, Song W, Sun Y, et al. High-Throughput, Multiplex Genotyping Directly from Blood or Dried Blood Spot without DNA Extraction for the Screening of Multiple G6PD Gene Variants at Risk for Drug-Induced Hemolysis. J Mol Diagn. 2017;19(5):638–50. Epub 2017/08/23. pmid:28826608.
- 3. Matsuda K. PCR-Based Detection Methods for Single-Nucleotide Polymorphism or Mutation: Real-Time PCR and Its Substantial Contribution Toward Technological Refinement. Adv Clin Chem. 2017;80:45–72. Epub 2017/04/23. pmid:28431642.
- 4. Shen GQ, Abdullah KG, Wang QK. The TaqMan method for SNP genotyping. Methods Mol Biol. 2009;578:293–306. Epub 2009/09/22. pmid:19768602.
- 5. Wu DY, Ugozzoli L, Pal BK, Wallace RB. Allele-specific enzymatic amplification of beta-globin genomic DNA for diagnosis of sickle cell anemia. Proc Natl Acad Sci U S A. 1989;86(8):2757–60. Epub 1989/04/01. pmid:2704745; PubMed Central PMCID: PMC286997.
- 6. Wittwer CT, Reed GH, Gundry CN, Vandersteen JG, Pryor RJ. High-resolution genotyping by amplicon melting analysis using LCGreen. Clin Chem. 2003;49(6 Pt 1):853–60. Epub 2003/05/27. pmid:12765979.
- 7. Li J, Wang L, Mamon H, Kulke MH, Berbeco R, Makrigiorgos GM. Replacing PCR with COLD-PCR enriches variant DNA sequences and redefines the sensitivity of genetic testing. Nat Med. 2008;14(5):579–84. Epub 2008/04/15. pmid:18408729.
- 8. Elek Z, Kovacs Z, Keszler G, Szabo M, Csanky E, Luo J, et al. High Throughput Multiplex SNP-analysis In Chronic Obstructive Pulmonary Disease and Lung CancerHigh Throughput Multiplex SNP-analysis In Chronic Obstructive Pulmonary Disease and Lung Cancer. Curr Mol Med. 2019. Epub 2019/10/24. pmid:31644407.
- 9. Hindson BJ, Ness KD, Masquelier DA, Belgrader P, Heredia NJ, Makarewicz AJ, et al. High-throughput droplet digital PCR system for absolute quantitation of DNA copy number. Anal Chem. 2011;83(22):8604–10. Epub 2011/11/01. pmid:22035192; PubMed Central PMCID: PMC3216358.
- 10. He C, Holme J, Anthony J. SNP genotyping: the KASP assay. Methods Mol Biol. 2014;1145:75–86. Epub 2014/05/13. pmid:24816661.
- 11. Steele KA, Quinton-Tulloch MJ, Amgai RB, Dhakal R, Khatiwada SP, Vyas D, et al. Accelerating public sector rice breeding with high-density KASP markers derived from whole genome sequencing of indica rice. Mol Breed. 2018;38(4):38. Epub 2018/03/23. pmid:29563850; PubMed Central PMCID: PMC5842261.
- 12. Ayalew H, Tsang PW, Chu C, Wang J, Liu S, Chen C, et al. Comparison of TaqMan, KASP and rhAmp SNP genotyping platforms in hexaploid wheat. PLoS One. 2019;14(5):e0217222. Epub 2019/05/23. pmid:31116793; PubMed Central PMCID: PMC6530864 following competing interests: Pak Wah Tsang, Junzhou Wang and Caifu Chen were employed by Integrated DNA Technologies at the time of study. Integrated DNA Technologies has patent application (20170145486, 20170253925, and 20170260583). This does not alter our adherence to PLOS ONE policies on sharing data and materials.
- 13. Khanal C, Kularathna MT, Ray JD, Stetina SR, McGawley EC, Overstreet C. Single Nucleotide Polymorphism Analysis Using KASP Assay Reveals Genetic Variability in Rotylenchulus reniformis. Plant Dis. 2019;103(8):1835–42. Epub 2019/06/14. pmid:31194618.
- 14. Shitaye G, Gadisa E, Grignard L, Shumie G, Chali W, Menberu T, et al. Low and heterogeneous prevalence of glucose-6-phosphate dehydrogenase deficiency in different settings in Ethiopia using phenotyping and genotyping approaches. Malar J. 2018;17(1):281. Epub 2018/08/04. pmid:30071859; PubMed Central PMCID: PMC6071387.
- 15. Zhang G, Xu Q, Liu J, Lv Z, Lu Y, Yang H, et al. Five P53 SNPs Involved in Low Rectal Cancer Risk and Prognosis in a Chinese Population. J Cancer. 2019;10(7):1772–80. Epub 2019/06/18. pmid:31205533; PubMed Central PMCID: PMC6547998.
- 16. Gascuel Q, Bordat A, Sallet E, Pouilly N, Carrere S, Roux F, et al. Effector Polymorphisms of the Sunflower Downy Mildew Pathogen Plasmopara halstedii and Their Use to Identify Pathotypes from Field Isolates. PLoS One. 2016;11(2):e0148513. Epub 2016/02/06. pmid:26845339; PubMed Central PMCID: PMC4742249.