20090521biotechniques Prim SNPing
20090521biotechniques Prim SNPing
net/publication/26251611
CITATIONS READS
28 2,927
7 authors, including:
All content following this page was uploaded by Yu-Huei Cheng on 20 May 2014.
SQL server database system. Java server the MySQL format and a local database Design, (4) Annealing Check, (5) RFLP
pages and Java applets are used to input copy, respectively. We plan to update Analysis, and (6) Output.
data and process files between the user these databases annually.
and the application, and to parse the data, Input module
respectively. The database structure is Program workflow Four kinds of primer design types (i.e.,
mainly set up by REBASE v. 607 (http:// The schematic program workflow of general PD, natural PD, mutagenic PD
rebase.neb.com) (18) and NCBI dbSNP Prim-SNPing (Figure 1) consists of six and CTPP PD) are provided. Except for
Build 125 (19), and is transformed into modules: (1) Input, (2) Query, (3) Primer the general PD, all other forms of PD in
Prim-SNPing provide the primers for SNP
genotyping in a cost-effective manner.
Sequence input types include SNP ID rs#
feed-in, sequence copy-and-paste, and text
file upload. Primer design constraints are
shown in detail in Figure 1.
Query module
When a user inputs the NCBI SNP ID rs#,
the sequence information for the targeted
SNP is retrieved from the local NCBI
dbSNP database. The contig position
is retrieved from NCBI’s Reference
Sequence (RefSeq) (20) to increase the
length of short flanking sequences, which
are usually provided in NCBI dbSNP.
The Prim-SNPing system can provide a
maximum flanking sequence length of
1000 bp for each SNP. This service for
retrieval of flanking sequence improves
the system’s flexibility for primer design.
Figure 1. Flowchart of web-based Prim-SNPing. Three kinds of functions are incorporated into the
input module of the Prim-SNPing system: the primer design (PD) type (general PD, natural PD, Primer design module
mutagenic PD and CTPP PD), the sequence input type (rs#, sequence and file input), and the The primer design module provides
primer design constraints, as indicated. The query module provides the SNP FASTA sequence and primers based on the settings of the input
flanking sequence nearby the SNP using rs# search via NCBI dbSNP (19). In a next step, the flank-
ing sequence of the SNP is fed into the primer design module for usage with the four PDs. Subse-
module (Figure 1). Default parameters
quently, the annealing check and RFLP analysis modules provide further analysis if selected. In the such as primer length (∼18–26 bp), primer
RFLP module, the SNP-containing sequences are transferred to a local database downloaded from length difference (5 bp), GC proportion
REBASE (18), and then the RFLP availability for the sense and antisense sequences is determined. (∼40–60%), Tm calculation, hairpin,
Finally, the output module displays the RFLP enzyme and primer information as indicated. Tm, melt- and GC clamp of primers are similar to
ing temperature; Tm-diff, Tm difference.
the ones used in the literature (21–25).
In general PD, primers are designed only
Table 1. Comparison of the Features of Primer Software Tools
for regular PCR without SNP genotyping
consideration. For natural PD, mutagenic
Primer3-
Primer Z
SNPbox
SNPing
WASP
Prim-
PDA
RFLP-mutagenic primers •
by (i) GC clamp, (ii) dimer check, (iii)
CTPP primers •
hairpin check, and (iv) BLAST specificity
Output data
By clicking the corresponding box in
Figure 4A, the designed primers can be
sorted by Tm difference and processed by
the BLAST function. The user-defined
primer conditions are shown in Figure
4B. Ten primers are provided for each
primer design (PD) in Prim-SNPing
(only one or two representative primers
are shown in Figure 4). In Figures 4C
and 4D, both the commercial and
Figure 2. Introduction of the PCR-CTPP method (modified from References 9 and 10). In the example, non-commercial restriction enzymes
the PCR-CTPP is designed to confront two-pair primers like F1, F2, R1, and R2, where F2 and R1 for SNPs in the sense (input) strand are
are opposite to the alternative nucleotide of the same SNP in both sense and antisense strands, mined from the local REBASE database
respectively. If the SNP is A/G in the sense strand and the primer F2 only covers SNP-A (sense
strand), the nucleotide in primer R1 is SNP-C (antisense strand) for its corresponding SNP-G in the (18) incorporated in Prim-SNPing. The
sense strand. (B) Primers F1/R2 are used for PCR as a positive control designed for both alternative BLAST function of a primer is acces-
nucleotides of the SNP (i.e., SNP-A and SNP-G). Primers F1/R1 and F2/R2 are preferentially used sible by clicking on it individually. The
for SNP-A and SNP-G, respectively. When electrophoresis is performed after the PCR, the AA type BLAST time for a primer is rather long
should be positive for PCR using F1/R1 and F2/R1 primers and the GG type should be positive for due to the short primer length. In Figures
PCR using F2/R2 and F2/R1 primers. For heterozygous type AG, PCR using F1/R1, F2/R2, and F2/
R1 primers should all be positive. Prim-SNPing provides primers to create different lengths for these
4C, 4D, and 4E, the Tm difference for the
three PCR reactions. Accordingly, the SNP genotype data obtained by PCR-CTPP primers contains primer designed in Prim-SNPing is very
additional information. small, ensuring a highly effective PCR
reaction. The primer sequence, product
for the secondary structure and specificity flanking sequence is adjustable in a range size, GC number, GC percentage, Tm,
of primers (Figure 1). These tests ensure of up to 1000 bp. The SNP flanking Tm difference, and primer length are
quality control of the designed primers. sequences provided by Prim-SNPing provided. In Figure 4D (mutagenic PD),
are pre-computed; they are longer than the position of the changed nucleotide
RFLP analysis module the SNP FASTA sequences provided (blue color) is adjustable. Figure 2 serves
The RFLP analysis module is constructed by NCBI dbSNP. This increased length as a reference to determine the relative
by REBASE as described (17). Natural is beneficial for the primer design. location of the primers designed in CTPP
and mutagenic primers are processed, For example, the length of SNP ID PD (Figure 4E). In brief, the primers F2
and the RFLP analysis module is used rs4930098 (http://www.ncbi.nlm.nih. and R1 as described (9) represent the
to verify the availability of restriction gov/SNP/snp_ref.cgi?rs=4930098) is default SNP-containing primers that
enzymes in the SNP sequence. 601 bp. When updating this SNP ID, confront each other for the alternative
Prim-SNPing provides a maximum of nucleotide of the SNP. The Prim-SNPing
Output module 1000 bp. The output results are reported system narrows the Tm difference down
Finally, information from the primer online and/or via email (Figure 3D). to within 1°C via computation, a value in
design modules, the RFLP analysis Before the primer design output, accordance with requirements listed in
module, and the annealing check module Prim-SNPing provides further analysis the PCR-CTPP literature (9). Finally, the
is integrated in the output module and using the advanced options listed in system provides the BLAST check for the
made available online and via email. Figure 3, E–G and an annealing check sequences of primers rather than BLAST
as shown in Figure 3H. Figure 3E shows check for the sequence of PCR amplicon
the common part for all primer designs (Figure 4F).
Results and discussion provided by Prim-SNPing. In principle,
Input data PCR-RFLP for SNP genotyping depends comparison of some primer design
Input data for Prim-SNPing is line-fed on the restriction enzyme availability software tools with the Prim-SNPing
through the web interface for the human, occurring in only one of the alternative Many primer design tools have been
rat or mouse SNP genotyping assay. The nucleotides of the targeted SNP (17). developed. Primer 3 (25) is one of the
input mode contains sequence and file Therefore, the lengths of the undigested earliest and most commonly used primer
inputs like the ones shown in Figure 3A. and digested PCR products vary design software tools. Primer Design
When clicking on the input button, the depending on the location of the primers. Assistant (PDA) (26) is a web interface
sequence data with the SNP in [dNTP1/ Prim-SNPing provides a novel function primer design service, which incorporates
dNTP2] or IUPAC format is fed into for the user-defined RFLP pattern for thermodynamic theory to evaluate the
the sequence window for primer design natural PD (Figure 3F). The “Ratio for fitness of primers. Recently, Primer3Plus
(Figure 3B). NCBI SNP ID rs# input RE (restriction enzyme)–digested RFLP (27), an enhanced web interface of
is also acceptable (Figure 3C) and the length” is set to 3:1 by default but can be Primer 3, was introduced. WASP (28)
References
1. Kwok, P.Y. 2000. Approaches to allele
frequency determination. Pharmacogenomics
1:231-235.
2. Kwok, P.Y. 2002. SNP genotyping with
fluorescence polarization detection. Hum.
Mutat. 19:315-323.
3. Olivier, M. 2005. The Invader assay for SNP
genotyping. Mutat. Res. 573:103-110.
4. Lamy, P., C.L. Andersen, F.P. Wikman,
and C. Wiuf. 2006. Genotyping and
annotation of Affymetrix SNP arrays. Nucleic
Acids Res. 34:e100.
5. Meaburn, E., L.M. Butcher, L.C. Schalkwyk,
and R. Plomin. 2006. Genotyping pooled DNA
using 100K SNP microarrays: a step towards
genomewide association scans. Nucleic Acids
Res. 34:e28.
6. Ahmadian, A., M. Ehn, and S. Hober. 2006.
Pyrosequencing: history, biochemistry and
future. Clin. Chim. Acta 363:83-94.
7. Chuang, L.Y., C.H. Yang, K.H. Tsui, Y.H.
Cheng, P.L. Chang, C.H. Wen, and H.W.
Chang. 2008. Restriction enzyme mining for
SNPs in genomes. Anticancer Res. 28:2001-
2007.
8. Atsuta, Y. and N. Hamajima. 2003. A new
Figure 4. Output data of Prim-SNPing. (A) List order and BLAST choice. The list order of primer genotyping method, PCR-CTPP. Rinsho Byori
appears randomly. Primers sorted by Tm difference (Tm-diff) are acceptable. BLAST function 51:886-891.
of primer is performed by clicking the respective button. (B) User-defined primer condition. 9. Hamajima, N. 2001. PCR-CTPP: a new
The conditions are the same as those for the user input in the primer design constraint under genotyping technique in the era of genetic
advanced options shown in Figure 3E. (C) Results of natural primer design (PD). The system epidemiology. Expert Rev. Mol. Diagn. 1:119-
shows only commercial and non-commercial restriction enzymes for the sense (input) strand 123.
only by default. Natural PD is similar to general PD, but general PD (not shown) is not capable of 10. Hamajima, N., T. Saito, K. Matsuo,
RFLP. (D) Results of mutagenic PD. The blue letter indicates the artificially changed nucleotide, K. Kozaki, T. Takahashi, and K. Tajima.
and the nucleotide of the IUPAC code indicates the SNP of interest (marked by red color). (E) 2000. Polymerase chain reaction with
Results of CTPP PD. The red-colored letters for nucleotides C and A at the 3′ ends of forward 1 confronting two-pair primers for polymor-
(F1) and reverse 2 (R2) primers indicates the alternative nucleotide of the same SNP in oppos- phism genotyping. Jpn. J. Cancer Res. 91:865-
ing orientation. (F) BLAST check for primer sequences. The system-designed primers can be 868.
selected individually to perform the BLAST function. 11. Hamajima, N., T. Saito, K. Matsuo,
and K. Tajima. 2002. Competitive ampli-
fication and unspecific amplification in
polymerase chain reaction with confronting
(9). We suggest that the melting enough, correct genotyping is not easy. two-pair primers. J. Mol. Diagn. 4:103-107.
temperature (Tm) should be the same Accordingly, the CTPP method elimi- 12. Linhart, C. and R. Shamir. 2002. The
degenerate primer design problem. Bioinfor-
for the four primers. However, PCR nates the need for a separate digestion matics 18(Suppl 1):S172-S181.
is often conducted using two or three step in RFLP, and is therefore faster and 13. Linhart, C. and R. Shamir. 2005. The
different annealing temperatures to cheaper to perform. degenerate primer design problem: theory and
determine ideal conditions and it is Prim-SNPing also provides a applications. J. Comput. Biol. 12:431-456.
difficult to predict which combination BLAST function for primer sequences, 14. Vallone, P.M. and J.M. Butler. 2004.
is superior for genotyping. The process so the specificity of designed primers AutoDimer: a screening tool for primer-dimer
and hairpin structures. BioTechniques 37:226-
is very similar to designing ordinary is controlled and improved. An added 231.
PCR workflows. Accordingly, the PCR benefit of Prim-SNPing is the integrated 15. Fredman, D., M. Jobs, L. Stromqvist,
protocol requires gradient PCR to be primer design function for cost-effective and A.J. Brookes. 2004. DFold: PCR design
performed in order to determine optimal SNP genotyping, including PCR-RFLP that minimizes secondary structure and
conditions. Specificity is also of impor- and PCR-CTPP. optimizes downstream genotyping applica-
tance when setting PCR-CTPP condi- tions. Hum. Mutat. 24:1-8.
16. Zhang, R., Z. Zhu, H. Zhu, T. Nguyen,
tions (9). Often, a weak band for the
nonexistent allele is observed. As long Acknowledgements F. Yao, K. Xia, D. Liang, and C. Liu. 2005.
SNP Cutter: a comprehensive tool for SNP
as the contrast of the bands between the This work was supported in part by the PCR-RFLP assay design. Nucleic Acids Res.
existing and nonexistent alleles is high National Science Council in Taiwan 33:W489-W492.
BioTechniques
2005. SNPbox: a modular software package for large-scale primer
design. Bioinformatics 21:385-387.
30. Tsai, M.F., Y.J. Lin, Y.C. Cheng, K.H. Lee, C.C. Huang,
Y.T. Chen, and A. Yao. 2007. PrimerZ: streamlined primer design
for promoters, exons and human SNPs. Nucleic Acids Res. 35:W63- E-Alerts
W65.
31. Niu, T. and Z. Hu. 2004. SNPicker: a graphical tool for primer
picking in designing mutagenic endonuclease restriction assays. Bioin-
BioTechniques Weekly
formatics 20:3263-3265. Weekly newsletter with industry insight for life
32. Kang, H.J., K.O. Choi, B.D. Kim, S. Kim, and Y.J. Kim. 2005. scientists: news, new products, grants deadlines,
FESD: a Functional Element SNPs Database in human. Nucleic Acids
Res. 33:D518-D522.
employment opportunities, and events.
Received 3 May 2008; accepted 7 January 2009. Advance Online Publication (AOP)
Address correspondence to Li-Yeh Chuang, Department of Chemical
Monthly alert with select research papers published
Engineering, I-Shou University, Taiwan, email: [email protected]; online ahead of the print.
or Cheng-Hong Yang, Department of Electronic Engineering, National
Kaohsiung University of Applied Sciences, Taiwan, email: chyang@
cc.kuas.edu.tw.
Table of Contents (TOC)
Monthly alert listing current papers and their links
allowing you to browse the latest issue immediately.
Register today at
www.BioTechniques.com/newsletters