mPSQed: a software for the design of multiplex pyrosequencing assays

Molecular-based diagnostic assays are the gold standard for infectious diseases today, since they allow a rapid and sensitive identification and typing of various pathogens. While PCR can be designed to be specific for a certain pathogen, a subsequent sequence analysis is frequently required for confirmation or typing. The design of appropriate PCR-based assays is a complex task, especially when conserved discriminating polymorphisms are rare or if the number of types which need to be differentiated is high. One extremely useful but underused method for this purpose is the multiplex pyrosequencing technique. Unfortunately there is no software available to aid researchers in designing multiplex pyrosequencing assays. Here, we present mPSQed (Multiplex PyroSeQuencing EDitor), a program targeted at closing this gap. We also present the design of an exemplarily theoretical assay for the differentiation of human adenovirus types A-F using two pyrosequencing primers on two distinct PCR products, designed quickly and easily using our software.     

Predicting the sensitivity and specificity of published real-time PCR assays

Background

In recent years real-time PCR has become a leading technique for nucleic acid detection and quantification. These assays have the potential to greatly enhance efficiency in the clinical laboratory. Choice of primer and probe sequences is critical for accurate diagnosis in the clinic, yet current primer/probe signature design strategies are limited, and signature evaluation methods are lacking.

Methods

We assessed the quality of a signature by predicting the number of true positive, false positive and false negative hits against all available public sequence data. We found real-time PCR signatures described in recent literature and used a BLAST search based approach to collect all hits to the primer-probe combinations that should be amplified by real-time PCR chemistry. We then compared our hits with the sequences in the NCBI taxonomy tree that the signature was designed to detect.

Results

We found that many published signatures have high specificity (almost no false positives) but low sensitivity (high false negative rate). Where high sensitivity is needed, we offer a revised methodology for signature design which may designate that multiple signatures are required to detect all sequenced strains. We use this methodology to produce new signatures that are predicted to have higher sensitivity and specificity.

Conclusion

We show that current methods for real-time PCR assay design have unacceptably low sensitivities for most clinical applications. Additionally, as new sequence data becomes available, old assays must be reassessed and redesigned. A standard protocol for both generating and assessing the quality of these assays is therefore of great value. Real-time PCR has the capacity to greatly improve clinical diagnostics. The improved assay design and evaluation methods presented herein will expedite adoption of this technique in the clinical lab.     

Not all sequence tags are created equal: designing and validating sequence identification tags robust to indels

Ligating adapters with unique synthetic oligonucleotide sequences (sequence tags) onto individual DNA samples before massively parallel sequencing is a popular and efficient way to obtain sequence data from many individual samples. Tag sequences should be numerous and sufficiently different to ensure sequencing, replication, and oligonucleotide synthesis errors do not cause tags to be unrecoverable or confused. However, many design approaches only protect against substitution errors during sequencing and extant tag sets contain too few tag sequences. We developed an open-source software package to validate sequence tags for conformance to two distance metrics and design sequence tags robust to indel and substitution errors. We use this software package to evaluate several commercial and non-commercial sequence tag sets, design several large sets (max_count = 7,198) of edit metric sequence tags having different lengths and degrees of error correction, and integrate a subset of these edit metric tags to polymerase chain reaction (PCR) primers and sequencing adapters. We validate a subset of these edit metric tagged PCR primers and sequencing adapters by sequencing on several platforms and subsequent comparison to commercially available alternatives. We find that several commonly used sets of sequence tags or design methodologies used to produce sequence tags do not meet the minimum expectations of their underlying distance metric, and we find that PCR primers and sequencing adapters incorporating edit metric sequence tags designed by our software package perform as well as their commercial counterparts. We suggest that researchers evaluate sequence tags prior to use or evaluate tags that they have been using. The sequence tag sets we design improve on extant sets because they are large, valid across the set, and robust to the suite of substitution, insertion, and deletion errors affecting massively parallel sequencing workflows on all currently used platforms.     

Set of novel tools for PCR primer design.

We have developed a new package of computer programs and algorithms for different PCR applications, including allele-specific PCR, multiplex PCR, and long PCR. The package is included in the upcoming VectorNTI suite software and attempts to incorporate most of the current knowledge about PCR primer design. A wide range of primer characteristics is available for user manipulation to provide improved efficiency and increased flexibility of primer design. To accelerate the primer calculations, we have optimized algorithms using recent advances in computer science such as dynamic trees and lazy evaluation. Proper structural organization of input parameters provides further program acceleration. New Vector NTI primer design software allows calculations of primer pairs for long PCR amplification of 120-kb genomic DNA in 5 min under most stringent input parameters and clustering 435 primer pairs for multiplex PCR within 30 min on a standard Pentium III PC. Our program allows the user to take advantage of molecule annotation by applying different kinds of filtering features during PCR primer design.     

EasyExonPrimer

EasyExonPrimer is a web-based software that automates the design of PCR primers to amplify exon sequences from genomic DNA. EasyExonPrimer is written in Perl and uses Primer3 to design PCR primers based on the genome builds and annotation databases available at the University of California, Santa Cruz (UCSC) Genome Browser database (http://genome.ucsc.edu/). It masks repeats and known single nucleotide polymorphism (SNP) sites in the genome and designs standardised primers using optimised conditions. Users can input genes by RefSeq mRNA ID, gene name or keyword. The primer design is optimised for large-scale resequencing of exons. For exons larger than 1 kb, the user has the option of breaking the exon sequence down into overlapping smaller fragments. All primer pairs are then verified using the In-Silico PCR software to test for uniqueness in the genome. We have designed >1000 pairs of primers for 90 genes; 95% of the primer pairs successfully amplified exon sequences under standard PCR conditions without requiring further optimisation. AVAILABILITY: EasyExonPrimer is available from http://129.43.22.27/~primer/. The source code is also available upon request. CONTACT: Xiaolin Wu (forestwu@mail.nih.gov).     

A robust PCR primer design platform applied to the detection of Acidobacteria Group 1 in soil

Environmental biosurveillance and microbial ecology studies use PCR-based assays to detect and quantify microbial taxa and gene sequences within a complex background of microorganisms. However, the fragmentary nature and growing quantity of DNA-sequence data make group-specific assay design challenging. We solved this problem by developing a software platform that enables PCR-assay design at an unprecedented scale. As a demonstration, we developed quantitative PCR assays for a globally widespread, ecologically important bacterial group in soil, Acidobacteria Group 1. A total of 33,684 Acidobacteria 16S rRNA gene sequences were used for assay design. Following 1 week of computation on a 376-core cluster, 83 assays were obtained. We validated the specificity of the top three assays, collectively predicted to detect 42% of the Acidobacteria Group 1 sequences, by PCR amplification and sequencing of DNA from soil. Based on previous analyses of 16S rRNA gene sequencing, Acidobacteria Group 1 species were expected to decrease in response to elevated atmospheric CO(2). Quantitative PCR results, using the Acidobacteria Group 1-specific PCR assays, confirmed the expected decrease and provided higher statistical confidence than the 16S rRNA gene-sequencing data. These results demonstrate a powerful capacity to address previously intractable assay design challenges.     

Greene SCPrimer: a rapid comprehensive tool for designing degenerate primers from multiple sequence alignments

Polymerase chain reaction (PCR) is widely applied in clinical and environmental microbiology. Primer design is key to the development of successful assays and is often performed manually by using multiple nucleic acid alignments. Few public software tools exist that allow comprehensive design of degenerate primers for large groups of related targets based on complex multiple sequence alignments. Here we present a method for designing such primers based on tree building followed by application of a set covering algorithm, and demonstrate its utility in compiling Multiplex PCR primer panels for detection and differentiation of viral pathogens.     

URPD: A Specific Product Primer Design Tool

ABSTRACT: BACKGROUND: Polymerase chain reaction (PCR) plays an important role in molecular biology. Primer design fundamentally determines its results. Here, we present a currently available software for a rather straight-forward way of visualizing the primer design process for infrequent users. RESULTS: URPD (yoUR Primer Design), a web-based specific product primer design tool, combines the NCBI Reference Sequences (RefSeq), UCSC In-Silico PCR, memetic algorithm (MA) and genetic algorithm (GA) primer design methods to obtain specific primer sets. A friendly user interface is accomplished by built-in parameter settings. The incorporated smooth pipeline operations effectively guide both occasional and advanced users. URPD contains an automated process, which produces feasible primer pairs that satisfy the specific needs of the experimental design with practical PCR amplifications. Visual virtual gel electrophoresis and in silico PCR provide a simulated PCR environment. The comparison of Practical gel electrophoresis comparison to virtual gel electrophoresis facilitates and verifies the PCR experiment. Wet-laboratory validation proved that the system provides feasible primers. CONCLUSIONS: URPD is a user-friendly tool that provides specific primer design results. The pipeline design path makes it easy to operate for beginners. URPD also provides a high throughput primer design function. Moreover, the advanced parameter settings assist sophisticated researchers in performing experiential PCR. Several novel functions, such as a nucleotide accession number template sequence input, local and global specificity estimation, primer pair redesign, user-interactive sequence scale selection, and virtual and practical PCR gel electrophoresis discrepancies have been developed and integrated into URPD. The URPD program is implemented in JAVA and freely available at http://bio.kuas.edu.tw/urpd/.     

Identification of Hedysarum varieties using amplified fragment length polymorphism on a capillary electrophoresis system.

For non-model organisms that do not have sequence information readily available, amplified fragment length polymorphism (AFLP) is a well-established technique that can be used for genomic mapping applications such as genetic diversity studies or phylogenetic studies. While AFLP can be performed on a variety of systems, including gel-based systems that require multiple labor-intensive steps, the availability of a more automated system that integrates the assay, electrophoresis platform, and analysis software could enable researchers to greatly increase their throughput and facilitate routine AFLP analysis. We demonstrate the use of such a system for AFLP analysis on Hedysarum species. AFLP assays performed on samples belonging to two different species isolated from Utah identified different varieties that clustered as expected from their actual locations.     

Java web tools for PCR, in silico PCR, and oligonucleotide assembly and analysis

The polymerase chain reaction is fundamental to molecular biology and is the most important practical molecular technique for the research laboratory. We have developed and tested efficient tools for PCR primer and probe design, which also predict oligonucleotide properties based on experimental studies of PCR efficiency. The tools provide comprehensive facilities for designing primers for most PCR applications and their combinations, including standard, multiplex, long-distance, inverse, real-time, unique, group-specific, bisulphite modification assays, Overlap-Extension PCR Multi-Fragment Assembly, as well as a programme to design oligonucleotide sets for long sequence assembly by ligase chain reaction. The in silico PCR primer or probe search includes comprehensive analyses of individual primers and primer pairs. It calculates the melting temperature for standard and degenerate oligonucleotides including LNA and other modifications, provides analyses for a set of primers with prediction of oligonucleotide properties, dimer and G-quadruplex detection, linguistic complexity, and provides a dilution and resuspension calculator.     

SNPkit: an efficient approach to systematic evaluation of candidate single nucleotide polymorphisms in public databases

There is widespread interest in devising genotyping methods for SNPs that are robust, inexpensive, and simple to perform. Although several high-throughput SNP genotyping technologies have been developed, including the oligonucleotide ligation assay, real-time PCR, and mass spectrometry, the issues of simplicity and cost-effectiveness have not been adequately addressed. Here we describe the application of a novel computer software package, SNPkit, which designs SNP genotyping assays based on a classical approach for discriminating alleles, restriction enzyme digestion. SNPkit can be used in genotyping assays for almost any SNPs including those that do not alter "natural" restriction sites. Using this method, 164 SNPs have been evaluated in DNA samples from 48 immortalized cell lines of randomly selected Chinese subjects. Sixty-two (37.8%) of the SNPs appeared to be common (frequencies of the minor alleles are > or = 5%) and were subsequently applied to a larger population-based sample. Overall, by using SNPkit, we have been able to validate and genotype accurately a large fraction of publicly available SNPs without sophisticated instrumentation.     

Computation-assisted SiteFinding- PCR for isolating flanking sequence tags in rice

SiteFinding-PCR is a method for isolating flanking sequence tags (FSTs) of T-DNA insertion lines, but the efficiency needs to be improved. Here we report a computation-assisted design for the random primers used in SiteFinding- PCR. A short sequence, GCATG, was screened from the rice genome and used as the 3' end of the random primer. When applying the optimized primer for isolating FSTs from 168 transgenic rice lines, we obtained 107 specific products, including 64 FSTs. The efficiency of obtaining FSTs using the modified version of SiteFinding-PCR increased by 73.0% compared with the method previously reported (P < 0.01, μ test). We also provide computational results for several other plant species such as maize, sorghum, Arabidopsis, foxtail millet, and Brachypodium based on the available genome data, so that the modified method could be easily adapted to other species.     

DNA probes for the identification of microorganisms

Summary The detection and identification of microorganisms is being carried out increasingly using DNA. Each organism has a unique DNA sequence which can be used to distinguish closely related organisms. Using PCR amplification and sequencing of ribosomal RNA genes we have developed DNA probes for a number of pathogenic bacteria and fungi. The development of DNA assays based on PCR has resulted in new questions which must be addressed including process carry-over contamination and inhibition of the PCR amplification reaction once the problems associated with the implementation of DNA assays are ironed out.     

PRaTo: a web-tool to select optimal primer pairs for qPCR

An essential pre-requisite to perform sound quantitative real-time polymerase chain reaction (qPCR) assays is to design outstanding primer pairs. This means they must have a good efficiency and be not prone to produce multiple amplicons or primer dimer products. To circumvent these issues, several softwares are available to help primer design. Although satisfactory computer-aided primer design tools are available for standard PCR, less efforts were done to provide specific methods for selection of optimal primer pairs for qPCR. We have developed PRaTo a web-based tool that enables checking and ranking of primers pairs for their attitude to perform optimally and reliably when used in qPCR experiments. PRaTo is available at http://prato.daapv.unipd.it.     

New softwares for automated microsatellite marker development

Microsatellites are repeated small sequence motifs that are highly polymorphic and abundant in the genomes of eukaryotes. Often they are the molecular markers of choice. To aid the development of microsatellite markers we have developed a module that integrates a program for the detection of microsatellites (TROLL), with the sequence assembly and analysis software, the Staden Package. The module has easily adjustable parameters for microsatellite lengths and base pair quality control. Starting with large datasets of unassembled sequence data in the form of chromatograms and/or text data, it enables the creation of a compact database consisting of the processed and assembled microsatellite containing sequences. For the final phase of primer design, we developed a program that accepts the multi-sequence ‘experiment file’ format as input and produces a list of primer pairs for amplification of microsatellite markers. The program can take into account the quality values of consensus bases, improving success rate of primer pairs in PCR. The software is freely available and simple to install in both Windows and Unix-based operating systems. Here we demonstrate the software by developing primer pairs for 427 new candidate markers for peanut.     

Fluorescent detection techniques for real-time multiplex strand specific detection of Bacillus anthracis using rapid PCR

Speed is a key area in our development of PCR assays for Bacillus anthracis. We believe that the strand specific detection of amplicons within 10 min is a realistic goal and that this will be achieved through fluorescent in-tube assays. We have used the Idaho LightCycler to study and develop candidate assays for B. anthracis. New strand specific fluorescent methods have been developed and a number of formats have been studied for speed and sensitivity. Internal controls have been developed as a method of improving our assay confidence. In this communication we will introduce the field of rapid PCR whilst discussing previous work in the areas described above, the development of our own rapid assay and a novel internal control system for B. anthracis. This work used PCR assays and hardware that are either commercially available, or have been previously described in open literature publications.     

Widespread use of real-time PCR for rickettsial diagnosis

We report 2?years of experience with rickettsial molecular diagnosis using real-time PCR at the French National Reference Center. All Rickettsia genomes available were compared to discover specific sequences to design new sets of primers and probes. The specificity was verified in silico and against a panel of 30 rickettsial species. Sensitivity was determined using 10-fold serial dilutions. Finally, primers and probes that were both specific and sensitive were routinely used for the diagnosis of rickettsial infections from clinical specimens. We retained sets of primers and probes to detect spotted fever group Rickettsia, typhus group Rickettsia,Rickettsia conorii,Rickettsia slovaca,Rickettsia africae and Rickettsia australis; 643 clinical samples were screened for the presence of Rickettsia DNA. Overall, 45 positive samples were detected, including 15 Rickettsia africae, nine R.?conorii, five Rickettsia sibirica mongolitimonae, four R.?slovaca, two R.?australis, four Rickettsia massiliae, one Rickettsia honei, one Rickettsia typhi and eight Rickettsia sp. Positive samples were detected mainly from cutaneous biopsies and swabs (31/45). Widespread use of real-time PCR is inexpensive and reduces delay in the diagnosis of rickettsial infections. These real-time PCR assays could be implemented easily in laboratories that have molecular facilities and may be added to existing molecular tools as a point-of-care strategy.     

Replacing PCR with COLD-PCR enriches variant DNA sequences and redefines the sensitivity of genetic testing

PCR is widely employed as the initial DNA amplification step for genetic testing. However, a key limitation of PCR-based methods is the inability to selectively amplify low levels of mutations in a wild-type background. As a result, downstream assays are limited in their ability to identify subtle genetic changes that can have a profound impact in clinical decision-making and outcome. Here we describe co-amplification at lower denaturation temperature PCR (COLD-PCR), a novel form of PCR that amplifies minority alleles selectively from mixtures of wild-type and mutation-containing sequences irrespective of the mutation type or position on the sequence. We replaced regular PCR with COLD-PCR before sequencing or genotyping assays to improve mutation detection sensitivity by up to 100-fold and identified new mutations in the genes encoding p53, KRAS and epidermal growth factor in heterogeneous cancer samples that had been missed by the currently used methods. For clinically relevant microdeletions, COLD-PCR enabled exclusive amplification and isolation of the mutants. COLD-PCR will transform the capabilities of PCR-based genetic testing, including applications in cancer, infectious diseases and prenatal identification of fetal alleles in maternal blood.