BOND: Basic OligoNucleotide Design
© Ilie et al.; licensee BioMed Central Ltd. 2013
Received: 8 August 2012
Accepted: 21 February 2013
Published: 27 February 2013
DNA microarrays have become ubiquitous in biological and medical research. The most difficult problem that needs to be solved is the design of DNA oligonucleotides that (i) are highly specific, that is, bind only to the intended target, (ii) cover the highest possible number of genes, that is, all genes that allow such unique regions, and (iii) are computed fast. None of the existing programs meet all these criteria.
We introduce a new approach with our software program BOND (Basic OligoNucleotide Design). According to Kane’s criteria for oligo design, BOND computes highly specific DNA oligonucleotides, for all the genes that admit unique probes, while running orders of magnitude faster than the existing programs. The same approach enables us to introduce also an evaluation procedure that correctly measures the quality of the oligonucleotides. Extensive comparison is performed to prove our claims. BOND is flexible, easy to use, requires no additional software, and is freely available for non-commercial use from http://www.csd.uwo.ca/~ilie/BOND/.
We provide an improved solution to the important problem of oligonucleotide design, including a thorough evaluation of oligo design programs. We hope BOND will become a useful tool for researchers in biological and medical sciences by making the microarray procedures faster and more accurate.
KeywordsDNA oligonucleotide design Microarray Spaced seeds
DNA microarrays are a very useful tool that has become ubiquitous in biological and medical research. Their ability to simultaneously measure expression levels of thousands of genes depends essentially on the availability of DNA oligonucleotides, or probes, that bind specifically to their targets. A DNA oligonucleotide, henceforth oligo, is a single-stranded piece of DNA that uniquely binds to a given region, called target (e.g., a gene). An oligo can be short (20-25bp) or long (50-70bp); we focus in this work on long oligos as they provide increased performance .
The problem of designing high quality oligos has received considerable attention by the research community and quite a few software programs have been developed for designing oligos [2-21]. The reader is referred to the survey of  for a detailed description of the best existing such programs.
The oligos designed by these programs have to meet several conditions. The first is sensitivity: oligos have to bind strongly to their targets at similar temperature. Sensitivity is achieved by ensuring that their melting temperatures  are not far from each other, with implications on the GC content, and by avoiding secondary structures that would make the oligos bind onto themselves instead of targets.
The second, and much more difficult to fulfill, condition is specificity, that is, the oligos must not bind to non-targets. Avoiding cross-hybridization is the most important, and the most difficult, issue for oligo selection. Precise criteria have been established, and widely adopted by the community, to ensure non-cross-hybridization [23, 24]; they are often referred to in the literature as Kane’s conditions:
(C75) the overall complementarity with non-targets should be less than 75% and
(C15) there should be no contiguous complementary region with non-targets of length 15 or more.
Technically, condition C75 means that the Hamming distance between the oligo and any non-target should be less than 75% of their length. Oligos satisfying these conditions will be called simply “good”, whereas those that do not will be called “bad.” The search for specific oligos implies detection and elimination of regions that are similar. Therefore, many programs employ BLAST , the most widely used tool for similarity search. We show that the specificity that can be achieved this way is quite low. While condition (C15) can be easily fulfilled using text indexes, as done by PICKY , condition (C75) is much more difficult. The only program that achieves good specificity is YODA . However, it employs an exhaustive search that is too slow to be used for most datasets.
Other conditions the oligo design programs have to meet are coverage, that is, designing oligos for as many genes as possible (note here that some genes cannot have their own oligos, due to high similarity with others) and speed; they should not only run in reasonable time but also be able to deal with datasets of any size. It turns out that none of the existing programs meet all these criteria.
We introduce a new approach with our software program BOND (Basic OligoNucleotide Design), that enables us to compute highly specific DNA oligos, for all the genes that admit unique probes, while running orders of magnitude faster than the competing programs. BOND is implemented in C/C++ and OpenMP, does not use any external software, is easy to use, and it allows the user to adjust a variety of parameters.
The same approach enables us to introduce as well an evaluation procedure that correctly measures the quality of the oligos. Extensive comparison is performed to prove our claims.
Oligo design, similarity search, and BLAST
An oligo binds to a non-target when they share considerable complementarity, as described by the conditions (C15) and (C75). That means, the reverse complement of the oligo and the non-target have substantial similarity and hence the search for specific oligos implies detection and elimination of regions that are similar. Therefore, many programs employ BLAST [25, 26], the most widely used tool for similarity search.
The BLAST program has been enormously successful for similarity search (the two papers have received together over 88,500 citations to date) however, it is ill suited for oligo design since in this case exhaustive search is needed and BLAST misses many similarities. Doubts about using BLAST for this purpose have been expressed before [11, 17] but the impact on the quality of the oligos produced has never been properly quantified. We provide such an evaluation in this section.
The above reasoning applies in principle to any hit-and-extend method that uses consecutive matches as hits, assuming the hits are at least as long as the BLAST seed; 16 out of 20 programs listed in  use some form of consecutive matching.
It has been noticed already in [27-29] that spaced matches provide better sensitivity but it was in [30-32] that research on finding the best spaced seeds with such properties was started. The first similarity search software program to use spaced seeds was PatternHunter ; it used the seed 111⋆1⋆⋆1⋆1⋆⋆11⋆111. For comparison, the contiguous seed of BLAST is 11111111111.
A spaced seed s is a sequence of 1’s (matches) and ⋆’s (don’t care positions). The number of 1’s is the weight of a seed and the total number of symbols is its length. Assuming a Bernoulli model , an alignment can be modelled as a sequence R of 1’s, for matches, and 0’s, for mismatches. A hit of s is a region of R of the same length as s such that 1’s in s correspond to 1’s in R. Given the length N of the alignment and the sequence identity p (percentage of expected 1’s in R), the sensitivity of s is the probability that s hits R.
A multiple spaced seed  consists of several spaced seeds. Its sensitivity is computed similarly. We say that a multiple spaced seed consisting of k seeds of weight g is a wg sk seed. For instance, PatternHunter has a w11s1 seed (one seed of weight 11) whereas PatternHunterII has a w11s16 seed (16 seeds of weight 11). Multiple spaced seeds can increase dramatically the sensitivity. As an example, the sensitivity of the BLAST’s, w11s1, and w11s16 seeds mentioned above, for N = 64 and p = 70 %, is 30%, 47%, and 92%, respectively.
Increasing the weight decreases the sensitivity but increases the specificity. To have both sensitivity and specificity high, we need to increase the number of seeds. The crucial observation  is that doubling the number of seeds has an effect on sensitivity that is similar to reducing the weight by one. Since the running time is proportional with the number of seeds and the number of random hits decreases four times when increasing the weight by one, both sensitivity and speed can be increased simultaneously.
For the design of BOND, we have employed multiple spaced seeds, as anticipated in the study of . Properly designed, multiple spaced seeds can achieve near 100% sensitivity and for this reason we used them in BOND. However, the distribution of the matches in seeds is by no means random and finding the optimal arrangement is a hard problem . A number of software programs for designing spaced seeds exist (most notably [36, 37]) but SpEED [38, 39] is the only polynomial-time algorithm that computes good (near optimal) seeds and we have adapted it for our needs. We have computed a w9s8 seed whose sensitivity is plotted in Figure 1(b) (BOND seed); its sensitivity is called “seed sensitivity” to avoid confusion with the sensitivity of the oligos). Its sensitivity at 75% identity is 98.3594% and the probability of detecting bad oligos of length 50 so that they respect condition (C75) is 99.8649%. (That is, the ratio between the area of the green region in Figure 1(b) and the total area of sequence identity 75% or higher; red and green together.) This is nearly perfect. We could have computed better seeds but this one did not make any error in all our tests and decreasing weight or increasing the number of seeds would slow down the program.
There are two reasons why this w9s8 seed is so effective. First, the sensitivity as given by the model is slightly underestimating the actual probability of the seed to find bad oligos since sequence identity below 75% is allowed in the model but not in our application. Second, we screen first the input data using a different w10s8 seed, also computed by SpEED, in order to quickly eliminate long repetitions.
There are several features that make BOND fast. First, in order to achieve the sensitivities of the seeds as computed above, we need to investigate every possible hit in the input data, for each seed. In this context, a hit of a seed is a pair of substrings of the same length as the seed which match on all positions corresponding to 1’s in the seed. This would be time consuming for the w9s8 seed. Therefore, in Phase I, we filter the input data for long repeats, using the w10s8 seed. Moreover, only pairs of consecutive hitting substrings are considered.
Assume L is the required oligo length. A region of length L will be called unique if it does not have 75% sequence identity or more with any other region of the same length. We keep a binary array, initialized with 0’s, where starting positions of non-unique regions are eliminated by marking them with 1’s. Phase I usually eliminates a substantial portion of the input. Phase II uses the more sensitive w9s8 seed. We still do not consider all remaining 0’s as it would still be too time consuming. Instead, we use an effective heuristic. The middle position of the longest contiguous stretch of 0’s is considered first to be tested by the w9s8 seed. All possible hits concerning the potential oligo ending at this position are considered, thus using the entire power of the seed. If this oligo is found to be bad, then the position is marked with a 1 and the algorithm continues trying to find a good oligo until either one is found or the gene is shown not to have one. It is important to notice that the order in which the genes are considered is irrelevant; this helps parallelizing the process.
Encode DNA sequences
Eliminate L-mers not satisfying (C 15)
Eliminate L-mers not satisfying GC content
Eliminate L-mers not satisfying (C 75) - Phase I
Compute melting temperature interval T
Eliminate L-mers with T m not in T
Check selected oligos for (C 75) - Phase II
Output oligos found, with genes and positions
L-mers are eliminated when found not to fulfill one of the required conditions. At any moment in the algorithm, the L-mers that were not eliminated so far are called valid. In step 1, the input data is processed so additional information is eliminated and only DNA sequences are kept. At the same time, all L-mers containing ambiguous or unknown nucleotides are eliminated. Step 2 eliminates L-mers containing 15-mers that occur elsewhere (condition (C15)). Step 3 eliminates L-mers with GC content outside the required range. Step 4 performs Phase I of (C75), using the w10s8 seed as explained above. For step 5, the melting temperatures (T m ) are computed for all currently valid L-mers. The interval T, of size Δ T m , containing the highest number of genes having valid L-mers with T m within T is computed by evaluating, in linear time, all intervals of the given size starting at any hundredth of degree between minimum and maximum T m . All L-mers having their melting temperature outside T are eliminated in step 6. Phase II of (C75) follows at step 7, using the full power of the w9s8 seed as explained above. Finally, at step 8, good oligos are reported.
Evaluating the quality of oligos involves solving problems similar with those faced during oligo design and the lack of tools for the latter implied the nonexistence of adequate procedures for the former. For instance, the evaluation of  uses BLAST with default weight 11, thus having the problems we have mentioned above. We have designed a similar approach for our evaluation procedure, however, in order to accurately measure the performance of oligo design programs, including BOND, we need a much more powerful tool. We have used again SpEED to design the evaluation seed with the estimated probability of error 0.0001%. (This is computed as above, as the ratio between the red area and red plus green together. The areas are not shown for this seed, as the red one would be too small.) The sensitivity of the three seeds are compared in Figure 1(c). This seed is estimated to make one error for each million oligos but, again, it is even more accurate in practice. It is thus safe to assume that it provides perfect evaluation.
All hits of this seed in each oligo are checked against the entire input. Due to very low weight of the seed, the evaluation is time consuming, but the results are quite reliable. It is interesting to notice that our evaluation is still considerably faster than the exhaustive search of YODA. This is due to the performance of the multiple spaced seed and efficient implementation.
The evaluation seed accurately classifies the oligos into “good” or “bad”, depending on whether they satisfy conditions (C75) and (C15) or not, respectively. The ratio between bad oligos and the total number of oligos gives the specificity. In order to evaluate coverage, that is, how many of the maximum possible number of genes can be associated with good oligos, we give an algorithm that approximates very well this maximum number of genes. Due to, again, lack of tools, coverage has not been measured before.
The approximation of the maximum number of genes that can have oligos, needed to evaluate coverage, is done as follows. What we need to find is the melting temperature interval that can produce the highest number of oligos. Steps 1-4 of BOND are performed first. Then, for each interval M of size Δ T m , the number of genes that contain a valid L-mer with melting temperature in M is computed. Those intervals for which this number is less than the number of oligos computed by BOND are eliminated, since they have no chance to perform better than BOND. An interval larger than Δ T m is left and step 6 of BOND is performed with this interval. From the remaining oligos, an extended Phase II is performed, in which all valid L-mers are checked using the w9s8 seed. Finally, the highest number of genes in an interval of size Δ T m is reported.
Both BOND and the evaluation algorithms are given a fast parallel implementation. There are two level of parallelization, one involving bit operations and the other multiple processors.
DNA sequences are stored using two bits per nucleotide. All L-mers including unknown and ambiguous positions are eliminated from the beginning. Seeds are stored as 64-bit integers and hash tables are used to store the positions of all possible spaced w-mers corresponding to each seed; given a seed of weight w, a spaced w-mer is any sequence of w letters in positions that corresponds to the matches in the seed. Each spaced w-mer is computed by AND-ing the mask of the seed with each position in the DNA sequences and stored as a 64-bit integers. A similar strategy is used for 15-substrings (condition (C15)). Since only elimination is performed, there is no need to use a text index, such as the suffix array to sort the 15-mers.
The most time consuming tasks are then performed in parallel. The hash tables for all seeds are computed in parallel for all seeds in each set. The uniqueness checking for all candidate oligos is done in parallel as well.
Results and discussion
Datasets and parameters
Datasets used for comparison
We have used the most common parameters for evaluating the programs: oligo length 50, GC content between 30..70%, size of melting temperature interval (Δ T m ) 10, maximum identity with non-target below 75%, and maximum consecutive matches allowed 14.
Comparison of oligo design programs
All programs were tested on the same machine, Dell Vostro, with Intel Core i7-2600 CPU (3.40GHz, 4 cores, 8 threads), 16 GB of RAM, and 8MB of cache running GNU Linux version 220.127.116.11-desktop-69mib; the compiler used was gcc 4.4.3. The more time consuming evaluation was performed on the SHARCNET high performance computers (http://www.sharcnet.ca).
Out of the six existing programs, ArrayOligoSelector, OligoArray, and OligoWiz showed low performance in all respects. OligoPicker has fair coverage but the specificity and speed are both low. PICKY is the fastest and the specificity is good, however the coverage is low. PICKY uses the suffix array  to correctly fulfill condition (C15) and achieves good speed. YODA has the highest specificity and coverage, which is expected due to the exhaustive search involved. However, its speed is the lowest among all programs. For instance, it took more than five days to complete the rice dataset.
BOND performs significantly better than all the other programs; under the testing conditions, it achieves 100% sensitivity and 100% coverage. All 236,112 oligos produced are good and the estimated maximum possible is 236,123, that is, BOND may have missed 11 oligos. Likely, this is not due to lack of sensitivity but to the slightly suboptimal decision on the melting temperature interval (see Methods). As far as speed is concerned, BOND is six times faster than PICKY, two orders of magnitude faster than OligoPicker, and over three orders of magnitude faster than YODA. BOND completed most tests under one minute.
BOND is freely available for non-commercial use from http://www.csd.uwo.ca/~ilie/BOND/. It requires no additional software and is easy to use. A wide range of input parameters allow the user to customize the oligos been produced in many ways: oligo length, GC content, size of melting temperature interval, oligos at 3’- or 5’-end, maximum number of consecutive matches allowed, maximum sequence identity allowed, secondary structure parameters. Details are given in the user manual provided in the same website.
We provide an improved solution to the important problem of oligonucleotide design. Our BOND software provides highly specific oligos, that cover the highest possible number of genes, while running orders of magnitude faster than existing software. We have provided also a thorough evaluation of oligo design programs.
We have focused here on the basic problem of finding gene-specific oligos. We plan to extend BOND in the near future to compute genome-wide tiling arrays; see  for references to such programs.
We hope BOND will become a useful tool for researchers in biological and medical sciences by making the microarray procedures faster and more accurate. We shall gladly receive their suggestions for further development.
Performance evaluation has been performed using the facilities of the Shared Hierarchical Academic Research Computing Network (SHARCNET: http://www.sharcnet.ca) and Compute/Calcul Canada. L.I., H.M., G.B.G., and W.F.S. have been partially supported by grants from the the Natural Sciences and Engineering Research Council of Canada (NSERC).
- Lemoine S, Combes F, Le Crom S: An evaluation of custom microarray applications: the oligonucleotide design challenge. Nucleic Acids Res 2009,37(6):1726-1739. 10.1093/nar/gkp053PubMed CentralView ArticlePubMedGoogle Scholar
- Li F, Stormo GD: Selection of optimal DNA oligos for gene expression arrays. Bioinformatics 2001,17(11):1067-1076. 10.1093/bioinformatics/17.11.1067View ArticlePubMedGoogle Scholar
- Chen H, Sharp BM: Oliz, a suite of Perl scripts that assist in the design of microarrays using 50mer oligonucleotides from the 3’ untranslated region. BMC Bioinformatics 2002, 3: 27. 10.1186/1471-2105-3-27PubMed CentralView ArticlePubMedGoogle Scholar
- Kaderali L, Schliep A: Selecting signature oligonucleotides to identify organisms using DNA arrays. Bioinformatics 2002,18(10):1340-1349. 10.1093/bioinformatics/18.10.1340View ArticlePubMedGoogle Scholar
- Mrowka R, Schuchhardt J, Gille C: Oligodb-interactive design of oligo DNA for transcription profiling of human genes. Bioinformatics 2002,18(12):1686-1687. 10.1093/bioinformatics/18.12.1686View ArticlePubMedGoogle Scholar
- Bozdech Z, Zhu J, Joachimiak M, Cohen F, Pulliam B, DeRisi J: Expression profiling of the schizont and trophozoite stages of Plasmodium falciparum with a long oligonucleotide microarray. Genome Biol 2003, 4: 1-15. 10.1186/gb-2003-4-2-p1View ArticleGoogle Scholar
- Emrich SJ, Lowe M, Delcher AL: PROBEmer: a web-based software tool for selecting optimal DNA oligos. Nucleic Acids Res 2003,31(13):3746-3750. 10.1093/nar/gkg569PubMed CentralView ArticlePubMedGoogle Scholar
- Rouillard J, Zuker M, Gulari E: OligoArray 2.0: design of oligonucleotide probes for DNA microarrays using a thermodynamic approach. Nucleic Acids Res 2003,31(12):3057-3062. 10.1093/nar/gkg426PubMed CentralView ArticlePubMedGoogle Scholar
- Talla E, Tekaia F, Brino L, Dujon B: A novel design of whole-genome microarray probes for Saccharomyces cerevisiae which minimizes cross-hybridization. BMC Genomics 2003, 4: 38. 10.1186/1471-2164-4-38PubMed CentralView ArticlePubMedGoogle Scholar
- Wang X, Seed B: Selection of oligonucleotide probes for protein coding sequences. Bioinformatics 2003,19(7):796-802. 10.1093/bioinformatics/btg086View ArticlePubMedGoogle Scholar
- Chou HH, Hsia AP, Mooney DL, Schnable PS: Picky: oligo microarray design for large genomes. Bioinformatics 2004,20(17):2893-2902. 10.1093/bioinformatics/bth347View ArticlePubMedGoogle Scholar
- Gordon PMK, Sensen CW: Osprey: a comprehensive tool employing novel methods for the design of oligonucleotides for DNA sequencing and microarrays. Nucleic Acids Res 2004,32(17):e133. 10.1093/nar/gnh127PubMed CentralView ArticlePubMedGoogle Scholar
- Hornshøj H, Stengaard H, Panitz F, Bendixen C: SEPON, a selection and evaluation pipeline for oligonucleotides based on ESTs with a non-target Tm algorithm for reducing cross-hybridization in microarray gene expression experiments. Bioinformatics 2004,20(3):428-429. 10.1093/bioinformatics/btg434View ArticlePubMedGoogle Scholar
- Reymond N, Charles H, Duret L, Calevro F, Beslon G, Fayard JM: ROSO: optimizing oligonucleotide probes for microarrays. Bioinformatics 2004,20(2):271-273. 10.1093/bioinformatics/btg401View ArticlePubMedGoogle Scholar
- Chung WH, Rhee SK, Wan XF, Bae JW, Quan ZX, Park YH: Design of long oligonucleotide probes for functional gene detection in a microbial community. Bioinformatics 2005,21(22):4092-4100. 10.1093/bioinformatics/bti673View ArticlePubMedGoogle Scholar
- Li X, He Z, Zhou J: Selection of optimal oligonucleotide probes for microarrays using multiple criteria, global alignment and parameter estimation. Nucleic Acids Res 2005,33(19):6114-6123. 10.1093/nar/gki914PubMed CentralView ArticlePubMedGoogle Scholar
- Nordberg EK: YODA: selecting signature oligonucleotides. Bioinformatics 2005,21(8):1365-1370. 10.1093/bioinformatics/bti182View ArticlePubMedGoogle Scholar
- Rimour S, Hill DRC, Militon C, Peyret P: GoArrays: highly dynamic and efficient microarray probe design. Bioinformatics 2005,21(7):1094-1103. 10.1093/bioinformatics/bti112View ArticlePubMedGoogle Scholar
- Rouchka EC, Khalyfa A, Cooper NGF: MPrime: efficient large scale multiple primer and oligonucleotide design for customized gene microarrays. BMC Bioinformatics 2005, 6: 175. 10.1186/1471-2105-6-175PubMed CentralView ArticlePubMedGoogle Scholar
- Wernersson R, Nielsen HB: OligoWiz 2.0 - integrating sequence feature annotation into the design of microarray probes. Nucleic Acids Res 2005,33(Web-Server-Issue):611-615.View ArticleGoogle Scholar
- Schretter C, Milinkovitch MC: OligoFaktory: a visual tool for interactive oligonucleotide design. Bioinformatics 2006, 22: 115-116. 10.1093/bioinformatics/bti728View ArticlePubMedGoogle Scholar
- SantaLucia J: A unified view of polymer, dumbbell, and oligonucleotide DNA nearest-neighbor thermodynamics. Proc Natl Acad Sci 1998,95(4):1460-1465. 10.1073/pnas.95.4.1460PubMed CentralView ArticlePubMedGoogle Scholar
- Kane MD, Jatkoe TA, Stumpf CR, Lu J, Thomas JD, Madore SJ: Assessment of the sensitivity and specificity of oligonucleotide (50mer) microarrays. Nucleic Acids Res 2000,28(22):4552-4557. 10.1093/nar/28.22.4552PubMed CentralView ArticlePubMedGoogle Scholar
- Hughes TR, Mao M, Jones AR, Burchard J, Marton MJ, Shannon KW, Lefkowitz SM, Ziman M, Schelter JM, Meyer MR, Kobayashi S, Davis C, Dai H, He YD, Stephaniants SB, Cavet G, Walker WL, West A, Coffey E, Shoemaker DD, Stoughton R, Blanchard AP, Friend SH, Linsley PS: Expression profiling using microarrays fabricated by an ink-jet oligonucleotide synthesizer. Nat Biotechnol 2001,19(4):342-347. 10.1038/86730View ArticlePubMedGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol 1990,215(3):403-410.View ArticlePubMedGoogle Scholar
- Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 1997,25(17):3389-3402. 10.1093/nar/25.17.3389PubMed CentralView ArticlePubMedGoogle Scholar
- Califano A, Rigoutsos I: FLASH: a Fast Look-Up Algorithm for String Homology. ISMB 1993, 56-64.Google Scholar
- Pevzner PA, Waterman MS: Multiple filtration and approximate pattern matching. Algorithmica 1995,13(1/2):135-154.View ArticleGoogle Scholar
- Lehtinen O, Sutinen E, Tarhio J: Experiments on block indexing. Proc. 3rd South American Workshop on String Processing 1996, 183-193.Google Scholar
- Burkhardt S, Kärkkäinen J: Better filtering with gapped q-Grams. CPM 2001, 73-85.Google Scholar
- Ma B, Tromp J, Li M: PatternHunter: faster and more sensitive homology search. Bioinformatics 2002,18(3):440-445. 10.1093/bioinformatics/18.3.440View ArticlePubMedGoogle Scholar
- Burkhardt S, Kärkkäinen J: Better filtering with gapped q-Grams. Fundam. Inform. 2003,56(1-2):51-70.Google Scholar
- Li M, Ma B, Kisman D, Tromp J: Patternhunter II: highly sensitive and fast homology search. J Bioinformatics Comput Biol 2004,2(3):417-440. 10.1142/S0219720004000661View ArticleGoogle Scholar
- Ilie L, Ilie S, Khoshraftar S, Mansouri-Bigvand A: Seeds for effective oligonucleotide design. BMC Genomics 2011, 12: 280. 10.1186/1471-2164-12-280PubMed CentralView ArticlePubMedGoogle Scholar
- Ma B, Yao H: Seed optimization is no easier than optimal Golomb ruler design. In Proceedings of the 6th Asia Pacific Bioinformatics Conference (APBC). London: Imperial College Press; 2008:133-144.Google Scholar
- Buhler J, Keich U, Sun Y: Designing seeds for similarity search in genomic DNA. In Proceedings of the Seventh Annual International Conference on Research in Computational Molecular Biology (RECOMB). New York: ACM; 2003:67-75.View ArticleGoogle Scholar
- Kucherov G, Noé L, Roytberg M: A unifying framework for seed sensitivity and its application to subset seeds. J Bioinformatics Comput Biol 2006,4(2):553-569. 10.1142/S0219720006001977View ArticleGoogle Scholar
- Ilie L, Ilie S: Multiple spaced seeds for homology search. Bioinformatics 2007,23(22):2969-2977. 10.1093/bioinformatics/btm422View ArticlePubMedGoogle Scholar
- Ilie L, Ilie S, Mansouri-Bigvand A: SpEED: fast computation of sensitive spaced seeds. Bioinformatics 2011,27(17):2433-2434. 10.1093/bioinformatics/btr368View ArticlePubMedGoogle Scholar
- Manber U, Myers EW: Suffix arrays: a new method for on-line string searches. SIAM J Comput 1993,22(5):935-948. 10.1137/0222058View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.