PRIMEGENS: robust and efficient design of gene-specific probes for microarray analysis

MOTIVATION: DNA microarray is a powerful high-throughput tool for studying gene function and regulatory networks. Due to the problem of potential cross hybridization, using full-length genes for microarray construction is not appropriate in some situations. A bioinformatic tool, PRIMEGENS, has recently been developed for the automatic design of PCR primers using DNA fragments that are specific to individual open reading frames (ORFs). RESULTS: PRIMEGENS first carries out a BLAST search for each target ORF against all other ORFs of the genome to quickly identify possible homologous sequences. Then it performs optimal sequence alignment between the target ORF and each of its homologous ORFs using dynamic programming. PRIMEGENS uses the sequence alignments to select gene- specific fragments, and then feeds the fragments to the Primer3 program to design primer pairs for PCR amplification. PRIMEGENS can be run from the command line on Unix/Linux platforms as a stand-alone package or it can be used from a Web interface. The program runs efficiently, and it takes a few seconds per sequence on a typical workstation. PCR primers specific to individual ORFs from Shewanella oneidensis MR-1 and Deinococcus radiodurans R1 have been designed. The PCR amplification results indicate that this method is very efficient and reliable for designing specific probes for microarray analysis.     

The degenerate primer design problem

A PCR primer sequence is called degenerate if some of its positions have several possible bases. The degeneracy of the primer is the number of unique sequence combinations it contains. We study the problem of designing a pair of primers with prescribed degeneracy that match a maximum number of given input sequences. Such problems occur when studying a family of genes that is known only in part, or is known in a related species. We prove that various simplified versions of the problem are hard, show the polynomiality of some restricted cases, and develop approximation algorithms for one variant. Based on these algorithms, we implemented a program called HYDEN for designing highly-degenerate primers for a set of genomic sequences. We report on the success of the program in an experimental scheme for identifying all human olfactory receptor (OR) genes. In that project, HYDEN was used to design primers with degeneracies up to 10(10) that amplified with high specificity many novel genes of that family, tripling the number of OR genes known at the time.     

G-PRIMER: greedy algorithm for selecting minimal primer set

G-PRIMER, a web-based primer design program, has been developed to compute a minimal primer set specifically annealed to all the open reading frames in a given microbial genome. This program has been successfully used in the microarray experiment for analyzing the expression of genes in the Xanthomonas campestris genome. AVAILABILITY: It is available at http://mammoth.bii.a-star.edu.sg/gprimer/. Its source code is available upon request.     

The degenerate primer design problem: theory and applications.

A PCR primer sequence is called degenerate if some of its positions have several possible bases. The degeneracy of the primer is the number of unique sequence combinations it contains. We study the problem of designing a pair of primers with prescribed degeneracy that match a maximum number of given input sequences. Such problems occur when studying a family of genes that is known only in part, or is known in a related species. We prove that various simplified versions of the problem are hard, show the polynomiality of some restricted cases, and develop approximation algorithms for one variant. Based on these algorithms, we implemented a program called HYDEN for designing highly degenerate primers for a set of genomic sequences. We report on the success of the program in several applications, one of which is an experimental scheme for identifying all human olfactory receptor (OR) genes. In that project, HYDEN was used to design primers with degeneracies up to 10(10) that amplified with high specificity many novel genes of that family, tripling the number of OR genes known at the time.     

Characterizing homologues of crop domestication genes in poorly described wild relatives by high‐throughput sequencing of whole genomes

Wild crop relatives represent a source of novel alleles for crop genetic improvement. Screening biodiversity for useful or diverse gene homologues has often been based upon the amplification of targeted genes using available sequence information to design primers that amplify the target gene region across species. The crucial requirement of this approach is the presence of sequences with sufficient conservation across species to allow for the design of universal primers. This approach is often not successful with diverse organisms or highly variable genes. Massively parallel sequencing (MPS) can quickly produce large amounts of sequence data and provides a viable option for characterizing homologues of known genes in poorly described genomes. MPS of genomic DNA was used to obtain species‐specific sequence information for 18 rice genes related to domestication characteristics in a wild relative of rice, Microlaena stipoides. Species‐specific primers were available for 16 genes compared with 12 genes using the universal primer method. The use of species‐specific primers had the potential to cover 92% of the sequence of these genes, while traditional universal primers could only be designed to cover 80%. A total of 24 species‐specific primer pairs were used to amplify gene homologues, and 11 primer pairs were successful in capturing six gene homologues. The 23 million, 36‐base pair (bp) paired end reads, equated to an average of 2X genome coverage, facilitated the successful amplification and sequencing of six target gene homologues, illustrating an important approach to the discovery of useful genes in wild crop relatives.     

Optimization of primer design for the detection of variable genomic lesions in cancer

Primer approximation multiplex PCR (PAMP) is a new experimental protocol for efficiently assaying structural variation in genomes. PAMP is particularly suited to cancer genomes where the precise breakpoints of alterations such as deletions or translocations vary between patients. The design of PCR primer sets for PAMP is challenging because a large number of primer pairs are required to detect alterations in the hundreds of kilobases range that can occur in cancer. These sets of primers must achieve high coverage of the region of interest, while avoiding primer dimers and satisfying the physico-chemical constraints of good PCR primers. We describe a natural formulation of these constraints as a combinatorial optimization problem. We show that the PAMP primer design problem is NP-hard, and design algorithms based on simulated annealing and integer programming, that provide good solutions to this problem in practice. The algorithms are applied to a test region around the known CDKN2A deletion, which show excellent results even in a 1:49 mixture of mutated:wild-type cells. We use these test results to help set design parameters for larger problems. We can achieve near-optimal designs for regions close to 1 Mb.     

UniPrime: a workflow-based platform for improved universal primer design

UniPrime is an open-source software (http://uniprime.batlab.eu), which automatically designs large sets of universal primers by simply inputting a gene ID reference. UniPrime automatically retrieves and aligns homologous sequences from GenBank, identifies regions of conservation within the alignment and generates suitable primers that can amplify variable genomic regions. UniPrime differs from previous automatic primer design programs in that all steps of primer design are automated, saved and are phylogenetically limited. We have experimentally verified the efficiency and success of this program by amplifying and sequencing four diverse genes (AOF2, EFEMP1, LRP6 and OAZ1) across multiple Orders of mammals. UniPrime is an experimentally validated, fully automated program that generates successful cross-species primers that take into account the biological aspects of the PCR.     

PriSM: a primer selection and matching tool for amplification and sequencing of viral genomes

SUMMARY: PriSM is a set of algorithms designed to select and match degenerate primer pairs for the amplification of viral genomes. The design of panels of hundreds of primer pairs takes just hours using this program, compared with days using a manual approach. PriSM allows for rapid in silico optimization of primers for downstream applications such as sequencing. As a validation, PriSM was used to create an amplification primer panel for human immunodeficiency virus (HIV) Clade B. AVAILABILITY: The program is freely available for use at: www.broadinstitute.org/perl/seq/specialprojects/primerDesign.cgi.     

Cross-Kingdom Amplification Using Bacteria-Specific Primers: Complications for Studies of Coral Microbial Ecology

PCR amplification of pure bacterial DNA is vital to the study of bacterial interactions with corals. Commonly used Bacteria-specific primers 8F and 27F paired with the universal primer 1492R amplify both eukaryotic and prokaryotic rRNA genes. An alternative primer set, 63F/1542R, is suggested to resolve this problem.     

Identification of Lactobacillus crispatus by polymerase chain reaction targeting S-layer protein gene

AIMS: This study aimed to develop a polymerase chain reaction (PCR) method to identify Lactobacillus crispatus. METHODS AND RESULTS: A primer set (CbsA2F-CbsA2R) for amplifying conserved regions of S-layer genes was designed to identify Lact. crispatus and the specificity of this set was compared with that of another primer set (Cri 16SI-Cri 16SII) which has been reported as a species-specific primer set targeting the 16S rRNA gene. Among species in the Lact. acidophilus A1-A4 groups, when KOD polymerase was used for amplification, the primer set CbsA2F-CbsA2R gave PCR products with Lact. crispatus strains only. However, when Taq polymerase was used, this primer set gave products with one Lact. amylovorus strain as well as with Lact. crispatus strains. The primer set Cri 16SI-Cri 16SII gave PCR products with Lact. crispatus strains and two Lact. acidophilus strains, regardless of whether the polymerase used was KOD or Taq. CONCLUSIONS: A PCR targeting the S-layer gene and amplified with KOD polymerase can identify Lact. crispatus accurately and rapidly. SIGNIFICANCE AND IMPACT OF THE STUDY: To the authors' knowledge, this is the first paper to provide a PCR method for the specific identification of Lact. crispatus.     

Characterization of monooxygenase gene diversity in benzene‐amended soils

Aim: To understand soil benzene monooxygenase gene diversity by clone library construction and microarray profiling. Methods and Results: A primer set was designed, and benzene monooxygenase gene diversity was characterized in two benzene‐amended soils. The dominant sequence types in the clone libraries were distinct between the two soils, and both sequences were assigned to novel clusters. Monooxygenase gene richness and diversity increased after benzene degradation. Oligonucleotide probes for microarray analysis were designed to detect a number of sequenced clones and reported monooxygenase genes. The microarray detected several genes that were not detected in the clone libraries of the same samples. Six probes were detected in more than one soil. Conclusions: The primer set designed in this study successfully detected diverse benzene monooxygenase genes. The level of diversity may have increased because the degradation of benzene differed from soil to soil. Microarrays have great potential in the comprehensive detection of gene richness as well as the elucidation of key genes for degradation. Significance and Impact of the Study: This study introduces a new primer set that may be used to identify diverse benzene monooxygenase genes in the environment; moreover, it demonstrates the potential of microarray technology in the profiling of environmental samples.     

GST-PRIME: a genome-wide primer design software for the generation of gene sequence tags

The availability of sequenced genomes has generated a need for experimental approaches that allow the simultaneous analysis of large, or even complete, sets of genes. To facilitate such analyses, we have developed GST-PRIME, a software package for retrieving and assembling gene sequences, even from complex genomes, using the NCBI public database, and then designing sets of primer pairs for use in gene amplification. Primers were designed by the program for the direct amplification of gene sequence tags (GSTs) from either genomic DNA or cDNA. Test runs of GST-PRIME on 2000 randomly selected Arabidopsis and Drosophila genes demonstrate that 93 and 88% of resulting GSTs, respectively, fulfilled imposed length criteria. GST-PRIME primer pairs were tested on a set of 1900 Arabidopsis genes coding for chloroplast-targeted proteins: 95% of the primer pairs used in PCRs with genomic DNA generated the correct amplicons. GST-PRIME can thus be reliably used for large-scale or specific amplification of intron-containing genes of multicellular eukaryotes.     

MultiPLX: automatic grouping and evaluation of PCR primers

SUMMARY: MultiPLX is a new program for automatic grouping of PCR primers. It can use many different parameters to estimate the compatibility of primers, such as primer-primer interactions, primer-product interactions, difference in melting temperatures, difference in product length and the risk of generating alternative products from the template. A unique feature of the MultiPLX is the ability to perform automatic grouping of large number (thousands) of primer pairs. AVAILABILITY: Binaries for Windows, Linux and Solaris are available from http://bioinfo.ebc.ee/download/. A graphical version with limited capabilities can be used through a web interface at http://bioinfo.ebc.ee/multiplx/. The source code of the program is available on request for academic users. CONTACT: maido.remm@ut.ee.     

The metabolic pathway of 2,4-dichlorophenoxyacetic acid degradation involves different families of tfdA and tfdB genes according to PCR-RFLP analysis

Abstract: Twenty-five 2,4-dichlorophenoxyacetic acid (2,4-D) degrading bacteria from geographically diverse locations and presenting various degrees of similarity or no similarity to the tfdA and tfdB genes from Alcaligenes eutrophus JMP134 were analysed by PCR-RFLP (restriction length fragment polymorphism). Primers for the 2,4-D etherase gene were derived by sequence alignment of the tfdA genes from A. eutrophus JMP134 and Burkholderia sp. RASC. Primers for the 2,4-dichlorophenolhydroxylase gene were based on the tfdB gene sequence from A. eutrophus JMP134 by taking codon degeneration and variations in amino acid residue sequences into consideration. PCR amplification using the tfdA primer set produced fragments of 0.3 kb from 17 strains which showed varying degrees of similarity to the tfdA gene probe from A. eutrophus JMP134. Significant variations in the gene sequences were confirmed by PCR-RFLP analysis. DNA amplification using the tfdB primer set produced a 1.1 kb fragment from 19 strains. Amongst them, two did not show any similarity to the tfdB gene probe. The size and restriction pattern of the products obtained from A. eutrophus JMP134 were in accordance with the expected size calculated from the A. eutrophus tfdA and tfdB gene sequence and their theoretical PCR-RFLP patterns. Some strains which did not amplify using the tfdA primer set did however amplify with the tfdB primer set. These results suggest the independent evolution of these two genes in the construction of the 2,4-D metabolic pathway. Our tfdA and tfdB primer sets could be used for the detection of similar sequences in bacteria and soils. Moreover, PCR-RFLP patterns could also be used to select subsets of strains for sequencing to study the phylogeny of the tfdA and tfdB genes.     

用于聚合酶链式反应(PCR)引物设计的计算机程序

 本文报道了两个用于PCR引物设计的计算机程序PCRDESN和PCRDESNA。PCRDESN程序主要从以下4个方面评价用户自己设计的一对引物的质量:(1)引物内的碱基反向重复或发夹结构,(2)两个引物之间的碱基互补配对,(3)两个引物之间的同源性,(4)引物的碱基组成及特点和T_m值计算。通过用多例文献发表的及本院有关实验室提供的引物对序列的验证,确定了程序的运算参数,证明该程序能较好地检验引物对的质量和解释某些PCR实验失败的原因。PCRDESNA程序采用逐级优化的方法和比PCRDESN所选用的更严紧的引物选择参数对用户提供的核酸序列进行快速检索,以确定所有可能的和合适的引物对。     

An improved microcomputer program for finding gene- or gene family-specific oligonucleotides suitable as primers for polymerase chain reactions or as probes

We present here an easy-to-use computer program which findsoligonucleotides suitable as primers in polymerase chain reactions(PCR) or as probes for hybridization. In contrast to other programsused for this purpose, the additional advantage of this oneis the possibility of directly detecting gene- as well as genefamily-specific oligonucleotides. For this purpose, up to 200different DNA sequences, of maximally 65 000 nucleotides each,can be scanned in a single search to ensure either single ormultiple gene binding of the PCR primers or probes. Specificoligonucleotides for genes carrying internal repetitions andfor single genes belonging to a set of highly conserved genescan also be detected. Many parameters such as exclusion of simplesequences, which are known to be highly repeated throughoutvarious genomes or regions of stable secondary structures inboth primer — primer and primer — template, canbe taken into consideration and avoided. Furthermore, the G+ C content and the length of the oligonucleotides can be changedin a broad range by the user. Received on January 14, 1991; accepted on March 14, 1991     

Oligonucleotide Primers for PCR Amplification of Coelomate Introns

Abstract Seven novel oligonucleotide primer pairs for polymerase chain reaction amplification of introns from nuclear genes in coelomates were designed and tested. Each pair bound to adjacent exons that are separated by a single intron in most coelomate species. The primer sets amplified introns in species as widely separated by the course of evolution as oysters (Mollusca: Protostoma) and salmon (Chordata: Deuterostoma). Each primer set was tested on a further 6 coelomate species and found to amplify introns in most cases. These primer sets may therefore be useful tools for developing nuclear DNA markers in diverse coelomate species for studies of population genetics, phylogenetics, or genome mapping.     

Design and Experimental Application of a Novel Non-Degenerate Universal Primer Set that Amplifies Prokaryotic 16S rRNA Genes with a Low Possibility to Amplify Eukaryotic rRNA Genes

The deep sequencing of 16S rRNA genes amplified by universal primers has revolutionized our understanding of microbial communities by allowing the characterization of the diversity of the uncultured majority. However, some universal primers also amplify eukaryotic rRNA genes, leading to a decrease in the efficiency of sequencing of prokaryotic 16S rRNA genes with possible mischaracterization of the diversity in the microbial community. In this study, we compared 16S rRNA gene sequences from genome-sequenced strains and identified candidates for non-degenerate universal primers that could be used for the amplification of prokaryotic 16S rRNA genes. The 50 identified candidates were investigated to calculate their coverage for prokaryotic and eukaryotic rRNA genes, including those from uncultured taxa and eukaryotic organelles, and a novel universal primer set, 342F-806R, covering many prokaryotic, but not eukaryotic, rRNA genes was identified. This primer set was validated by the amplification of 16S rRNA genes from a soil metagenomic sample and subsequent pyrosequencing using the Roche 454 platform. The same sample was also used for pyrosequencing of the amplicons by employing a commonly used primer set, 338F-533R, and for shotgun metagenomic sequencing using the Illumina platform. Our comparison of the taxonomic compositions inferred by the three sequencing experiments indicated that the non-degenerate 342F-806R primer set can characterize the taxonomic composition of the microbial community without substantial bias, and is highly expected to be applicable to the analysis of a wide variety of microbial communities.     

Aromatic hydrocarbon degradation genes from chronically polluted Subantarctic marine sediments

Aim:  The goal of this study was to identify functional targets to detect polycyclic aromatic hydrocarbon (PAH)-degrading bacterial populations in cold marine ecosystems.
Methods and Results:  We designed a degenerate primer set targeting genes encoding the α subunit of PAH-dioxygenases from Gram-positive bacteria. This primer set was used to amplify gene fragments from metagenomic DNA isolated from Subantarctic marine sediments (Ushuaia Bay, Argentina). These gene fragments were cloned and sequenced. We identified 14 distinct groups of genes, most of them showing significant relatedness with dioxygenases from Gram-positive bacteria of the genera Rhodococcus , Mycobacterium , Nocardioides , Terrabacter and Bacillus . The level of identity with these genes, however, was low to moderate (33–62% at the amino acid level).
Conclusion:  These results indicate the presence of a high diversity of hitherto unidentified dioxygenase genes in this cold polluted environment.
Significance and Impact of the Study:  Subantarctic marine ecosystems are particularly vulnerable to hydrocarbon pollution, and the development of environmental restoration strategies for these environments is pressing. The information obtained in this work will be the starting point for the design of quantitative molecular tools to analyse the abundance and dynamics of these aromatic hydrocarbon-degrading bacterial populations in the marine environment.