Nucleic Acids Research, 2003, Vol. 31, No. 13 3491-3496
© 2003 Oxford University Press
Design of oligonucleotides for microarrays and perspectives for design of multi-transcriptome arrays
Center for Biological Sequence Analysis, BioCentrum-DTU, Technical University of Denmark, Building 208, 2800 Lyngby, Denmark
*To whom correspondence should be addressed. Tel: +45 45252489; Fax: +45 45931585; Email: hbjorn{at}cbs.dtu.dk
Received January 25, 2003; Revised and Accepted April 7, 2003
 |
ABSTRACT |
|---|
 TOP ABSTRACT INTRODUCTIONMETHODDESIGNING OLIGONUCLEOTIDESRESULTSDISCUSSIONREFERENCES |
|---|
Optimal design of oligonucleotides for microarrays involves tedious and laborious work evaluating potential oligonucleotides relative to a series of parameters. The currently available tools for this purpose are limited in their flexibility and do not present the oligonucleotide designer with an overview of these parameters. We present here a flexible tool named OligoWiz for designing oligonucleotides for multiple purposes. OligoWiz presents a set of parameter scores in a graphical interface to facilitate an overview for the user. Additional custom parameter scores can easily be added to the program to extend the default parameters: homology,
Tm, low-complexity, position and GATC-only. Furthermore we present an analysis of the limitations in designing oligonucleotide sets that can detect transcripts from multiple organisms. OligoWiz is available at www.cbs.dtu.dk/services/OligoWiz/. |
INTRODUCTION |
|---|
|
TOPABSTRACTINTRODUCTIONMETHODDESIGNING OLIGONUCLEOTIDESRESULTSDISCUSSIONREFERENCES |
|---|
Oligonucleotides of 50–70 bp, or even shorter in the range of 20–30 bp, are now widely used for microarrays in gene expression studies. An ideal oligonucleotide must discriminate well between its intended target and all other targets in the pool. It must detect concentration differences under the hybridization conditions with minimal variation and detect within the range of its target concentration under the given conditions. However in reality the optimal oligonucleotide will often have to be a compromise between these criteria.
Here we present a flexible program for designing optimal oligonucleotides for microarrays named OligoWiz. The program will allow for the design of highly specialized arrays as well as general use arrays. OligoWiz evaluates and graphically presents all potential oligonucleotides in the input sequences according to a collection of parameters, each of which can be assigned a different weight. In addition to a standard collection of parameters (homology, Tm, GATC-only, position within transcript and complexity), custom parameters may be of relevance to specialized microarrays. Therefore OligoWiz has been designed to be flexible towards additional user defined parameters. We present an analysis of the limitations in designing oligonucleotide sets that can detect transcripts from multiple organisms.
OligoWiz is implemented as a client–server solution. The server is responsible for the calculation of the scores. The client is used for submitting jobs to the server as well as for visualizing the scores and fine-tuning the placement of the oligonucleotides. The OligoWiz client is freely available from the OligoWiz web page: http://www.cbs.dtu.dk/services/OligoWiz/ and the server is commercially available.
|
METHOD |
|---|
|
TOPABSTRACTINTRODUCTIONMETHODDESIGNING OLIGONUCLEOTIDESRESULTSDISCUSSIONREFERENCES |
|---|
The OligoWiz client is written in Java 1.3.1 with cross-platform functionality in mind. It is tested on MacOS X, Linux and Windows, and should be able to run on any platform with Java 1.3.1 or better.
The OligoWiz server was developed on a SGI Unix system and is written in Perl5. The server utilizes the BLAST program for homology search (1) and saco_patterns for generating a pattern information database (2). The server program is parallelized using the Perl module ChildManager. Additional information about the sequence databases available for OligoWiz can be found on the OligoWiz web page.
Multi-transcriptome analysis
Predicted genes from two strains of Campylobacter jejuni (RM1221 from TIGR and NCTC11168 from GenBank) were extracted based on EasyGene annotations with R-values <2 (Larsen and Krogh, manuscript submitted). The rat, mouse, Streptococcaceae pneumoniae and Streptococcaceae pyogenes coding DNA sequences (CDS) were extracted based on the ftp://ftp.ncbi.nih.gov/genbank/genomes/M_musculus/RNA/rna.gbk.gz, ftp://ftp.ncbi.nih.gov/genbank/genomes/R_norvegicus/rn_mrna.gz files, NC_003028
[GenBank]
and NC_002737
[GenBank]
annotations, respectively.
Patterns were extracted and redundancy reduced for stretches of 15 bp using saco_patterns (2) and custom Perl scripts. The oligonucleotides were further redundancy reduced using BLAST (1). The oligonucleotides were mapped back on the transcriptomes using Perl scripts.
|
DESIGNING OLIGONUCLEOTIDES |
|---|
|
TOPABSTRACTINTRODUCTIONMETHODDESIGNING OLIGONUCLEOTIDESRESULTSDISCUSSIONREFERENCES |
|---|
The protocol for oligonucleotide design using OligoWiz is as follows: the program calculates a series of independent parameter scores that describes all potential oligonucleotides. Each parameter score has a value between 0 and 1, where 1 is optimal. For each oligonucleotide, the parameter scores are weighted and summed to form a total score, which is normalized to a value between 0 and 1. The graphical interface presents the user with graphs of all the parameter scores including the total score, along a virtual sequence (Fig. 1). OligoWiz points out the oligonucleotide with the highest total score in each input sequence, but the graphical interface also allows for direct user intervention. In addition to the scores the user interface can also visualize sequence annotations in a color bar. An example could be exon/intron annotation. The user can then select an oligonucleotide based on the information presented. Below we describe the parameter scores that are used for the oligonucleotide design.
|
Detailed descriptions of how to use the program are found at: www.cbs.dtu.dk/services/OligoWiz/instructions.html.
Cross-hybridization
In order for the oligonucleotide to avoid cross-hybridization, the affinity difference between the intended target and all other targets should ideally be maximized. However, such a calculation requires that the concentration of all targets are known, which they rarely are. Even if the concentrations were ignored, an affinity calculation would be very time consuming for large collections of oligonucleotides. In practice, a homology search between all targets and the oligonucleotide is a good estimate of the affinity of the oligonucleotide to alternative targets. Experimental evidence suggests that a significantly false signal can be detected if a 50 bp oligonucleotide has >75–80% of the bases complementary or if continuous stretches of >15 bp are complementary to a false target (3). Similar results have been found for 60mers (4).
To estimate the degree of cross-hybridization, OligoWiz calculates a homology score for each oligonucleotide in a given sequence, based on a BLAST search of the input sequence against a species-specific database. The resulting BLAST hits are then evaluated relative to each oligo along the input sequence.
Let m be the number of BLAST hits considered in position i of the oligonucleotide and h={h1i,..., hmi} be the BLAST hits in position i.
 |
To avoid including BLAST hits against the gene sequence that the oligonucleotide originated from OligoWiz allows for filtering against BLAST hits that cover more than a fraction (default 0.8) of the entire length of the input sequence and have a homology higher than a threshold (default 97%). The filter also enables filtering of homology lower than a threshold (default 70%) or shorter than a given length (default 15 bp).
Tm
To further enhance the ability of the oligonucleotides to discriminate between the targets, the hybridization and washing conditions need to be optimal. Therefore it is crucial that all oligonucleotides perform well under similar hybridization conditions. The melting temperature of the DNA:DNA duplex (Tm) is a good description of an oligonucleotide hybridization property and as such the minimal difference between the Tm of the oligonucleotides is desired.
OligoWiz uses a nearest-neighbor model for Tm estimation:
 |
H is the enthalpy, S is the entropy change of the nucleation reaction; A is a constant correcting for helix initiation (-10.8), R is the universal gas constant (1.99 cal K-1 mol-1), Ct is the total molar concentration of strands (5). Since the total molar concentration of strands is unknown for most microarray experiments, OligoWiz uses a constant of 2.5x10-10 M. Based on the Tm estimation a Tm score is calculated.
 |
To facilitate maximal freedom in selecting oligonucleotides according to other requirements, OligoWiz calculates melting temperatures for all oligonucleotides in a user defined length interval. For each 5' position along the input sequence the oligonucleotide length (extending toward the 3' end) with the best Tm score is chosen. This narrows the Tm difference distribution of the oligonucleotides significantly as can be seen in Figure 2. Thus, the Tm score and oligonucleotide length are connected. The oligonucleotide length influences the remaining score calculations (except the position score). Therefore the Tm score is the first calculation the OligoWiz server performs.
|
Position within transcript
The position within the target transcript can be of importance, especially for long transcripts. The reverse transcriptase will fall off the transcript with a certain probability and the further away from the starting point the less signal will be generated. Therefore oligonucleotides with target positions closer to the starting point of the reverse transcriptase are preferable.
If the labeling commences from the 3' end (poly A tail) the following score is used:
 |
3'end is the oligonucleotide distance to the 3' end of the input sequence.
In cases where the labeling is done with random primers, as would be the case under prokaryote mRNA labeling, the chance of having an oligonucleotide upstream of a given position should be accounted for:
 |
Low-complexity filtering
To avoid oligonucleotides composed of very common sequence fragments in probe design a low-complexity score was implemented. Different sequences are common in different species, therefore to estimate a low-complexity measure for an oligonucleotide a list of sequence subfragments and the information content of these is generated specifically for each species. The information content can be calculated by the following equation (6):
 |
 |
A low-complexity score of 0 indicates an oligonucleotide with very low complexity. By far the majority of oligonucleotides have a low-complexity score between 1 and 0.8, and in fact the distribution is exponential. As a consequence the low-complexity score only significantly influences few oligonucleotides, with very low low-complexity measures.
GATC-only score
To allow for filtering out sequence containing ambiguity annotation OligoWiz has a score called ‘GATC-only’. Oligonucleotides containing R, Y, M, K, X, S, W, H, B, V, D, N or anything else but A, C, T or G will be given a score of 0, while oligonucleotides only containing G, A, T and C will be assigned a score of 1.
Custom parameter scores
In addition to the parameter scores calculated by the OligoWiz server, custom scores can be added to the list and will be handled by the graphically interface. (For detailed description see the program documentation.)
|
RESULTS |
|---|
|
TOPABSTRACTINTRODUCTIONMETHODDESIGNING OLIGONUCLEOTIDESRESULTSDISCUSSIONREFERENCES |
|---|
To evaluate OligoWiz, oligonucleotides were designed to detect the expression from the 6600 genes annotated in the Saccharomyces cerevisiae genome (7). The design was performed using standard settings of OligoWiz, designing oligonucleotides in the length interval 45–55 bp. The homology search and complexity score was based on whole genome databases. The mean Tm of the oligonucleotides was 75.7°C. OligoWiz was able to find oligonucleotides with good parameters for almost all genes as shown in Table 1.
|
Less than 150 oligonucleotides had low complexity. The calculations of the parameters for the 6600 yeast genes were done in just 20 min.
Multi-transcriptome arrays
Microarray experiments using oligonucleotides covering a whole transcriptome are still very expensive to initiate and it is of course best to reduce these costs if possible. One possible way to reduce the setup cost for several organisms is to design oligonucleotides that enable detection of transcripts from more than one organism. Furthermore a microarray that has been designed to be robust toward strain differences is of interest for laboratories that focus on more than one strain of a species. To shed light on this topic the following analysis was engaged.
If it is assumed that the sequences in two transcriptomes are random and independent, the expected number of oligonucleotides of length L that will be perfectly complementary to a sequence in each of two transcriptomes of lengths N1 and N2 would be approximately N1 * N2/4L. If we further restrict the oligonucleotides not to have more than 70% of the sequence in common with any other sequence within the transcriptomes and to be unique in any stretch of length k (k being less than 70% of L). Then the expected number of oligonucleotides that meet these criteria can be approximated by:
 |
A shared 50mer where any of its 36 possible 15mers are unique within each transcriptome and overall does not share >70% homology to any other 50mer (3), is unlikely even in large transcriptomes like the mouse and rat transcriptomes. Furthermore, we find that the expected number of such shared oligonucleotides in random transcriptomes is very dependent on the length of the oligonucleotide.
However, transcriptomes are not composed of random sequences and related transcriptomes are not independent. To investigate the impact of such biases the number of 50mers in the transcriptomes of rat and mouse as well as in S.pneumoniae and S.pyogenes and the two C.jejuni strains (NCTC11168 and RM1221) were counted. Hereof only few were shared within the Streptococcaceae (<0.1%) and the rodents (<2%) species relative to the 55–59% between the two C.jejuni strains, indicating how related the two C.jejuni strains are. Still all the compared groups shared far more oligonucleotides than expected for random transcriptomes with the respective sizes.
The theoretical model predicts that only a very moderate fraction of the oligonucleotides will be rejected by the redundancy criteria. In reality we find that 73, 34 and 11% of the Streptococcaceae, C.jejuni and rodent oligonucleotides pass the redundancy criteria respectively (i.e. the oligonucleotide does not share any stretch with continuous complementary of 15 bp or more, and have <70% overall complementary, to any other transcript sequence from either transcriptome). This is considerably less than expected for a random sequence, indicating some redundancy in the sequences. The differences in redundancy between the groups may to some degree be explained by the AT skew in the C.jejuni transcriptomes (68%) and the large size of the rodent transcriptomes.
The shared and non-redundant 50mers were mapped back on the transcriptome and the transcripts that could be measured by these were counted. A total of 619 out of the 9538 and 9487 transcripts in mouse and rat, respectively, and only 23 transcripts out of the 2094 and 1697 S.pneumoniae and S.pyogenes transcriptomes, respectively. In the two highly related C.jejuni strains on the other hand 88 and 76% or 1437 transcripts out of 1640 and 1880 were found to be measurable by common oligonucleotides.
To investigate the potential for shorter oligonucleotides, all oligonucleotides of 25 bp in the transcriptome pairs were likewise counted. As expected more conserved 25mers than 50mers were found, especially between the more distantly related Streptococcaceae species (Table 2). Using the same redundancy criteria as for the 50mers reduced the potential numbers of 25mers relatively less than was the case for the 50mers. This is expected as a 25mer consists of only 11 stretches of 15 bp, in contrast to the 36 stretches in a 50mer. The shared and non-redundant 25mers were mapped back on the transcriptome and the transcripts that could be measured by these were counted (Table 3). Approximately twice as many rodent and more than seven times more Streptococcaceae transcripts were measurable by common 25mers.
|
|
Finally, we investigated the potential of even shorter oligonucleotides of 18 bp to allow for designing multi-transcriptome arrays for mouse and rat. It was possible to detect 2901 transcripts with common 18mers.
The shorter oligonucleotides enabled more transcripts to be measured by common oligonucleotides, but the impact of the oligonucleotide length does not have the impact expected based on the random transcriptome model. This suggests that the majority of the common oligonucleotides originate from long stretches that are conserved between the transcriptomes.
To demonstrate the flexibility of OligoWiz, we added a custom score indicating oligonucleotides that were suitable for multi-transcriptome arrays. The custom score was assigned a 1 where oligonucleotides were shared and not redundant between two transcriptomes, otherwise a 0 was assigned. The custom scores were added as an extra column in the input file for the OligoWiz client program (ASCII file) and the score was named in the file header. The OligoWiz client automatically included the new parameter and allowed for oligonucleotide selection using the new parameter.
|
DISCUSSION |
|---|
|
TOPABSTRACTINTRODUCTIONMETHODDESIGNING OLIGONUCLEOTIDESRESULTSDISCUSSIONREFERENCES |
|---|
The design of oligonucleotides for microarray analysis is a laborious and tedious work if the analysis is not automated in some way. A couple of tools for this purpose already exist (8–10). Li and Stormo (8) have developed a program that focuses on avoiding cross-hybridization. OligoArray by Rouillard et al. (9) is, as OligoWiz, a client server application that runs locally and sends jobs to a web server. However OligoArray requires a local BLAST function. Mrowka et al. (10) describes a method to design oligonucleotides for the human ENSEMBL project entries. Neither of these tools supports extensive user interaction and they give only limited overview of the design process.
The purposes of microarrays are numerous and, in many cases, this has an impact on the design. OligoWiz is designed to have maximum degrees of freedom and to give the designer an overview of a number of parameters in a graphical display. In particular, OligoWiz gives the designer extended freedom by minimizing the Tm of the oligonucleotides by varying their length within an interval.
Additional custom parameters can easily be added to the parameter collection and displayed in the graphical display. This may be a phred score to indicate sequence quality (11); indications of intron or exon regions or some other parameter the designer finds important for the microarray design at hand.
The amount of shared oligonucleotides between the three transcriptome pairs (mouse/rat, S.pneumoniae/S.pyogenes and two C.jejuni strains NCTC11168/RM1221) indicates how related the transcriptomes are (Table 2). This by far exceeds the amount expected for random transcriptomes. However, the fraction of oligonucleotides shared between the transcriptomes is a major bottleneck for the success of multi-transcriptome array design. Only relatively few potential multi-transcriptome oligonucleotides are rejected because of redundancy: 27–66% in the two bacterial groups and 89% in the larger rodent transcriptomes. For transcriptomes with many paralogs this may be of greater importance. To the extent this is caused by paralogs, grouping paralogs, and thereby accepting not to be able to distinguish these in the resulting array, can diminish loss of potential oligonucleotide fragments.
Our analysis also indicates the impact the length of the oligonucleotides has on the possibility for designing multi-transcriptome microarrays. Furthermore, shorter oligonucleotides may be better at discriminating between perfect complementary sequence and sequences with mismatches. This is indicated by the increasing difference in melting temperature between such sequences, as they get shorter (data not shown).
Alternative hybridization chemistries that will allow for shorter oligonucleotides to be used as microarray probes like LNA (reviewed in 12) may allow for multi-transcriptome arrays of more distantly related organisms. However our analysis indicates the strong limitations this concept has. But, for very related organisms, it is possible to save a number of oligonucleotides and make the array measurements more strain robust for a considerable number of genes.
|
ACKNOWLEDGEMENTS |
|---|
The authors wish to thank L. Jensen, A. Fausbøll, U. de Lichtenberg and T. Jensen for discussion on the theoretical multi-transcriptome model; Peter Hallin for discussions on multi-transcriptome design, and David Ussery for proofreading the manuscript. This work was funded by grants from Novozymes A/S, The Danish National Research Foundation and The Danish Biotechnology Instrument Center.
|
REFERENCES |
|---|
|
TOPABSTRACTINTRODUCTIONMETHODDESIGNING OLIGONUCLEOTIDESRESULTSDISCUSSIONREFERENCES |
|---|
- Altschul,S.F., Gish,W., Miller,W., Myers,E.W. and Lipman,D.J. (1990) Basic local alignment search tool. J. Mol. Biol., 215, 403–410.[CrossRef][ISI][Medline]
- Jensen,L.J. and Knudsen,S. (2000) Automatic discovery of regulatory patterns in promoter regions based on whole cell expression data and functional annotation. Bioinformatics, 16, 326–333.
[Abstract/Free Full Text] - Kane,M.D., Jatkoe,T.A., Stumpf,C.R., Lu,J., Thomas,J.D. and Madore,S.J. (2000) Assessment of the sensitivity and specificity of oligonucleotide (50mer) microarrays. Nucleic Acids Res., 28, 4552–4557.
[Abstract/Free Full Text] - Hughes,T.R., Mao,M., Jones.A.R., Burchard,J., Marton,M.J., Shannon,K.W., Lefkowitz,S.M., Ziman,M., Schelter,J.M., Meyer,M.R. et al. (2001) Expression profiling using microarrays fabricated by an ink-jet oligonucleotide synthesizer. Nat. Biotechnol., 19, 342–347.[CrossRef][ISI][Medline]
- Breslauer,K.J., Frank,R., Blocker,H. and Marky,L.A. (1986) Predicting DNA duplex stability from the base sequence. Proc. Natl Acad. Sci. USA, 83, 3746–3750.
[Abstract/Free Full Text] - Shannon,C.E. (1948) The mathematical theory of communication. Bell System Tech. J., 27, 379–423.
- Goffeau,A., Barrell,B.G., Bussey,H., Davis,R.W., Dujon,B., Feldmann,H., Galibert,F., Hoheisel,J.D., Jacq,C., Johnston,M. et al. (1996) Life with 6000 genes. Science, 274, 563–567.
- Li,F. and Stormo,G.D. (2001) Selection of optimal DNA oligos for gene expression arrays. Bioinformatics, 17, 1067–1076.
[Abstract/Free Full Text] - Rouillard,J.M., Herbert,C.J. and Zuker,M. (2002) OligoArray: genome-scale oligonucleotide design for microarrays. Bioinformatics, 18, 486–487.
[Abstract/Free Full Text] - Mrowka,R., Schuchhardt,J. and Gille,C. (2002) Oligodb-interactive design of oligo DNA for transcription profiling of human genes. Bioinformatics, 18, 1686–1687.
[Abstract/Free Full Text] - Ewing,B. and Green,P. (1998) Base-calling of automated sequencer traces using phred. II. Error probabilities. Genome Res., 8, 186–194.
[Abstract/Free Full Text] - Braasch,D.A. and Corey,D.R. (2001) Locked nucleic acid (LNA): fine-tuning the recognition of DNA and RNA. Chem. Biol., 8, 1–7.[CrossRef][ISI][Medline]
CiteULike 
Connotea 
Del.icio.us What's this?
This article has been cited by other articles:
|
|
 |
|
  T. B. Rasmussen, M. Danielsen, O. Valina, C. Garrigues, E. Johansen, and M. B. Pedersen Streptococcus thermophilus Core Genome: Comparative Genome Hybridization Study of 47 Strains Appl. Envir. Microbiol., August 1, 2008; 74(15): 4703 - 4710. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 M. B. Pedersen, C. Garrigues, K. Tuphile, C. Brun, K. Vido, M. Bennedsen, H. Mollgaard, P. Gaudu, and A. Gruss Impact of Aeration and Heme-Activated Respiration on Lactococcus lactis Gene Expression: Identification of a Heme-Responsive Operon J. Bacteriol., July 15, 2008; 190(14): 4903 - 4911. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 M. R. Andersen, W. Vongsangnak, G. Panagiotou, M. P. Salazar, L. Lehmann, and J. Nielsen A trispecies Aspergillus microarray: Comparative transcriptomics of three Aspergillus species PNAS, March 18, 2008; 105(11): 4387 - 4392. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
C. J. Paredes, R. S. Senger, I. S. Spath, J. R. Borden, R. Sillers, and E. T. Papoutsakis A General Framework for Designing and Validating Oligomer-Based DNA Microarrays and Its Application to Clostridium acetobutylicum Appl. Envir. Microbiol., July 15, 2007; 73(14): 4631 - 4638. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 D. Pyeon, M. A. Newton, P. F. Lambert, J. A. den Boon, S. Sengupta, C. J. Marsit, C. D. Woodworth, J. P. Connor, T. H. Haugen, E. M. Smith, et al. Fundamental Differences in Cell Cycle Deregulation in Human Papillomavirus-Positive and Human Papillomavirus-Negative Head/Neck and Cervical Cancers Cancer Res., May 15, 2007; 67(10): 4605 - 4619. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 S. Feng and E. R.M. Tillier A fast and flexible approach to oligonucleotide probe design for genomes and gene families Bioinformatics, May 15, 2007; 23(10): 1195 - 1202. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 I. Brukner, R. El-Ramahi, I. Gorska-Flipot, M. Krajinovic, and D. Labuda An in vitro selection scheme for oligonucleotide probes to discriminate between closely related DNA sequences Nucleic Acids Res., May 14, 2007; 35(9): e66 - e66. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
G. M. Bruun, R. Wernersson, A. S. Juncker, H. Willenbrock, and H. B. Nielsen Improving comparability between microarray probe signals by thermodynamic intensity correction Nucleic Acids Res., April 1, 2007; 35(7): e48 - e48. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
H. Willenbrock, A. Petersen, C. Sekse, K. Kiil, Y. Wasteson, and D. W. Ussery Design of a Seven-Genome Escherichia coli Microarray for Comparative Genomic Profiling J. Bacteriol., November 15, 2006; 188(22): 7713 - 7721. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
L. A. Marcelino, V. Backman, A. Donaldson, C. Steadman, J. R. Thompson, S. P. Preheim, C. Lien, E. Lim, D. Veneziano, and M. F. Polz Accurately quantifying low-abundant targets amid similar sequences by revealing hidden correlations in oligonucleotide microarray data PNAS, September 12, 2006; 103(37): 13629 - 13634. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 M. Mehlmann, E. D. Dawson, M. B. Townsend, J. A. Smagala, C. L. Moore, C. B. Smith, N. J. Cox, R. D. Kuchta, and K. L. Rowlen Robust Sequence Selection Method Used To Develop the FluChip Diagnostic Microarray for Influenza Virus. J. Clin. Microbiol., August 1, 2006; 44(8): 2857 - 2862. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
T. Yamada, H. Soma, and S. Morishita PrimerStation: a highly specific multiplex genomic PCR primer design server for the human genome. Nucleic Acids Res., July 1, 2006; 34(Web Server issue): W665 - W669. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
W.-H. Chung, S.-K. Rhee, X.-F. Wan, J.-W. Bae, Z.-X. Quan, and Y.-H. Park Design of long oligonucleotide probes for functional gene detection in a microbial community Bioinformatics, November 15, 2005; 21(22): 4092 - 4100. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
X. Li, Z. He, and J. Zhou Selection of optimal oligonucleotide probes for microarrays using multiple criteria, global alignment and parameter estimation Nucleic Acids Res., October 24, 2005; 33(19): 6114 - 6123. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
Z. He, L. Wu, X. Li, M. W. Fields, and J. Zhou Empirical Establishment of Oligonucleotide Probe Design Criteria Appl. Envir. Microbiol., July 1, 2005; 71(7): 3753 - 3760. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
R. Wernersson and H. B. Nielsen OligoWiz 2.0--integrating sequence feature annotation into the design of microarray probes Nucleic Acids Res., July 1, 2005; 33(suppl_2): W611 - W615. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
E. K. Nordberg YODA: selecting signature oligonucleotides Bioinformatics, April 15, 2005; 21(8): 1365 - 1370. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
S. Rimour, D. Hill, C. Militon, and P. Peyret GoArrays: highly dynamic and efficient microarray probe design Bioinformatics, April 1, 2005; 21(7): 1094 - 1103. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 P. Hilson, J. Allemeersch, T. Altmann, S. Aubourg, A. Avon, J. Beynon, R. P. Bhalerao, F. Bitton, M. Caboche, B. Cannoot, et al. Versatile Gene-Specific Sequence Tags for Arabidopsis Functional Genomics: Transcript Profiling and Reverse Genetics Applications Genome Res., October 1, 2004; 14(10b): 2176 - 2189. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
P. M. K. Gordon and C. W. Sensen Osprey: a comprehensive tool employing novel methods for the design of oligonucleotides for DNA sequencing and microarrays Nucleic Acids Res., September 29, 2004; 32(17): e133 - e133. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
P. Pancoska, Z. Moravek, and U. M. Moll Rational design of DNA sequences for nanotechnology, microarrays and molecular computers using Eulerian graphs Nucleic Acids Res., August 27, 2004; 32(15): 4630 - 4645. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 B. C. Meyers, D. W. Galbraith, T. Nelson, and V. Agrawal Methods for Transcriptional Profiling in Plants. Be Fruitful and Replicate Plant Physiology, June 1, 2004; 135(2): 637 - 652. [Full Text] [PDF]
|
|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||