- Research article
- Open Access
A generic approach for the design of whole-genome oligoarrays, validated for genomotyping, deletion mapping and gene expression analysis on Staphylococcus aureus
BMC Genomics volume 6, Article number: 95 (2005)
DNA microarray technology is widely used to determine the expression levels of thousands of genes in a single experiment, for a broad range of organisms. Optimal design of immobilized nucleic acids has a direct impact on the reliability of microarray results. However, despite small genome size and complexity, prokaryotic organisms are not frequently studied to validate selected bioinformatics approaches. Relying on parameters shown to affect the hybridization of nucleic acids, we designed freely available software and validated experimentally its performance on the bacterial pathogen Staphylococcus aureus.
We describe an efficient procedure for selecting 40–60 mer oligonucleotide probes combining optimal thermodynamic properties with high target specificity, suitable for genomic studies of microbial species. The algorithm for filtering probes from extensive oligonucleotides libraries fitting standard thermodynamic criteria includes positional information of predicted target-probe binding regions. This algorithm efficiently selected probes recognizing homologous gene targets across three different sequenced genomes of Staphylococcus aureus. BLAST analysis of the final selection of 5,427 probes yielded >97%, 93%, and 81% of Staphylococcus aureus genome coverage in strains N315, Mu50, and COL, respectively. A manufactured oligoarray including a subset of control Escherichia coli probes was validated for applications in the fields of comparative genomics and molecular epidemiology, mapping of deletion mutations and transcription profiling.
This generic chip-design process merging sequence information from several related genomes improves genome coverage even in conserved regions.
Current hybridization technologies allow assaying thousands of nucleic acid sequences in a single reaction on a solid substrate. Such massively parallel systems offer unprecedented opportunities for basic research and diagnostic applications, including gene sequencing , detection of genetic polymorphisms , genome-composition analysis [3, 4] and measurement of gene expression profiles in prokaryotes [5, 6] or cancer cells . Oligonucleotide probes (up to 70-mer) offer more flexibility than cDNA probes since they can be tailored according to optimal in silico physico-chemical and specificity properties, and applied to any sequence data.
Early available probe design software identified sets of probes sharing homogeneous thermodynamic properties for probe-target hybridization . More elaborated software tools include cross-homology testing of probes against a reference database by BLAST (Basic Local Alignment Search Tool) [9, 10] or prediction of secondary structures into the thermodynamically-based approach [11–14]. A frequent drawback of some of these algorithms is to yield an excessive number of unprocessed BLAST outputs that complicates final selection of the most specific probes. Furthermore, these approaches do not take into consideration probe interaction with microarray surface, in particular the impact of mismatches position between the target and probes, as shown by Hughes et al . Designing reliable oligonucleotide probes with available software is quite difficult for bacterial genomes with low GC content , low complexity in sequence composition, or frequent conserved repeats leading to erroneous target identification by cross-hybridization.
The reported method (OliCheck) implements an algorithm for filtering oligonucleotide probes libraries sharing homogeneous thermodynamic properties by using positional information of predicted target-probe binding regions. An additional characteristic of OliCheck is to annotate probes recognizing highly conserved targets shared by different genomes. Staphylococcus aureus (S. aureus) was selected as a model organism for implementing and experimentally validating this approach. The choice of this clinically important pathogen for fundamental and applied genomic studies is prompted by the availability of several fully or partially sequenced strain genomes [16–18]. A set of feature elements was designed by OliCheck to yield an extensive S. aureus genome coverage. This S. aureus specific probe set together with control probes were used to manufacture an oligoarray that was extensively validated for comparative genomics, molecular epidemiology, mapping of deletion mutations, and transcription profiling applications. The specificity, signal-response linearity, and influence of hybridization temperatures for transcript profiling are also described. Further genomic oligoarrays of several distinct microbial species have been successfully designed using this generic methodological approach.
Design of a S. aureus oligoarray
The major genomic component of S. aureus is a 2.8 million base pairs (bp) circular chromosome with a low GC content (32.8%) which represents 2,595 ORFs in strain N315. The median length of strain N315 ORFs is 768 bp and the 25th–75th percentile ranges from 444-1,152 bp . A probe design software (i.e. ArrayDesigner™) generated a comprehensive list of candidate 40–60 mer probes for each ORF of strain N315 according to the preset thermodynamic hybridization parameters (Fig. 1, step A). A similar process was performed with strain Mu50 and COL. This step yielded 417,776, 607,461 and 321,286 candidate probes for strains N315, Mu50 and COL, respectively.
Further selection of candidate oligonucleotides to sort out cross-reactive probes within the genome of each strain was achieved by OliCheck. This specificity filtering test is a mathematical transposition of Hughes et al observations  on the impact of mismatches position between the target and probe with respect to the microarray surface. The occurrence of mismatch(es) in the distal portion (solution end) of the probe leads to a strong decrease in signal intensity, as opposed to the proximal portion (surface end). The results of the specificity filtering test, when performed separately on each strain, yielded a list of 48,415, 48,510, and 34,303 probes for N315, Mu50 and COL, respectively (Fig. 1, step B). In contrast to OliCheck-filtered probes that are expected to be devoid of cross-hybridization, >85% of the probes selected by ArrayDesigner™ alone displayed significant cross-hybridization with multiple ORFs.
For each strain, all accepted probes were further annotated by OliCheck against heterologous S. aureus genomes to identify probes common to the different genomes (Fig. 1, step C).
To fulfill the manufacturing requirements of a S. aureus genome-wide oligoarray, a further probe selection was performed. This selection used a spreadsheet program to rank probes according to optimal strain coverage and thermodynamic criteria, for providing one or more non-overlapping probes per target ORF (Fig. 1, step D).
In silico properties of the S. aureus oligoarray and manufacturing of StaphChip
The final set of 5,335 S. aureus OliCheck-filtered probes recognized 97.5, 93.0, and 81.0% of N315, Mu50, and COL ORFs, respectively. The low residual percentage of ORFs (2.5% for strain N315 and 7.0% for Mu50) that escaped recognition by our final probe set mostly included (51/65 ORFs for N315) mobile genetic elements, located either on prophage or transposon elements. A likely explanation for exclusion of these targets by OliCheck is the occurrence of highly homologous (>98% nucleotide identity) sequence repeats found in other ORFs. Accordingly, 92 residual probes covering these ORFs had to be selected relying on step A (Fig. 1) only.
To manufacture StaphChip, a total of 5,427 S. aureus probes were synthesized on the array together with a subset of control probes designed by OliCheck against E. coli K12 genome  (see methods).
Comparison of OliCheck with a currently used method
To validate the properties of OliCheck by comparison with an established method, we generated 60-mer oligonucleotides probes with homogenous thermodynamic criteria for N315 genome by the recently published ArrayOligoSelector  tool. The final set of probes generated according to the ArrayOligoSelector procedure  (n = 2,592) was further tested for cross-homology by using OliCheck. A large percentage of these probes (83%) were sorted out by OliCheck because they failed specificity criteria defined by Hughes et al .
Comparison of in silico-predicted with experimentally detected hybridization signals
Among the total number of 5,427 S. aureus probes, 4,812 (89%) recognized targets common to all three strains (Fig. 2), thus affording StaphChip valuable properties for comparative genomic and transcriptomic studies on S. aureus. The finding of a significant number of probes (n = 401) commonly identified in N315 and Mu50 but not COL may be explained by their closely-related genomic contents .
The number of probes predicted to detect genomic elements from each S. aureus strain was compared to those experimentally recorded using StaphChip (Table 1). For each strain, >99.5% of the in silico- predicted hybridization signals were indeed detected by hybridization of genomic DNA on StaphChip probes. Most of the false-positive signals were recorded on probes that were generated using ArrayDesigner™ only (n = 92/104). In silico analysis determined that these false-positive signals did not originate from recognition of intergenic regions (data not shown). Thus, such false-positive frequency is not transferable to the whole array.
Use of StaphChip for deletion mapping and genomotyping applications
To evaluate the accuracy of StaphChip for deletion mapping, the Cy-3 labelled DNA from the SA113ica deletion mutant  was competitively hybridized with the Cy-5 labelled DNA from its isogenic parent. Figure 3 maps the ica-related signals which are missing in the ica mutant, as opposed to the unique positive signal generated by the tetracycline resistance marker used for the construction and selection of strain SA113ica. Fluorescent intensities of all other signals except ica- related genes were equivalent for both strains.
The potential of StaphChip for epidemiological typing was analysed by comparative genomic hybridization (CGH) of S. aureus strains from various epidemiological origins. Genomic DNA from each individual S. aureus strain labelled with Cy3 was co-hybridized with equivalent amounts of Cy5-labelled control genomic DNA (pooled from N315, Mu50 and COL) and analyzed by two-way hierarchical clustering (Fig. 4).
The recently sequenced community-acquired MRSA strain MW2 , not included in StaphChip probe design, revealed important differences with strains N315, Mu50, and COL, as shown by a major yellow region on Figure 4. The set of probes yielding negative signals with MW2 DNA (Fig. 4, cluster a) corresponds to ORFs coding for antibiotic and heavy metal resistance determinants, bacterial adhesins and DNA modification enzymes lacking in this epidemic strain . In contrast, extensively conserved genomic regions (Fig. 4, cluster b) are composed of house-keeping genes contributing to cell-wall, DNA synthesis, as well as essential metabolic enzymes.
Further analysis of strain SA113 and SA113ica compared to the other strains, showed no hybridization signals (Fig. 4, cluster c) for a number of genes (e.g. exotoxins and antibiotic resistance determinants) known to be missing in pathogenicity islands I-II-III of the NCTC8325 family .
Application of StaphChip for expression profiling studies
Control experiments were performed to study the dose-response of labelled cDNA and influence of hybridization temperature on the linearity and intensities of fluorescent signals. Two portions of N315 cDNA were labelled during reverse-transcription, with either Cy-3 or Cy-5 and competitively hybridized on StaphChip. Characteristics of fluorescent signals obtained on N315 specific probes were compared to those recorded on control E. coli oligonucleotide probes, at either 55 or 60°C. The median level of fluorescence intensities were approximately 5–10 fold higher for S. aureus probes compared to E. coli probes.
On S. aureus capture elements, log-transformed signal intensities recorded with equivalent input amounts (5 or 10 μg) of Cy-3 and Cy-5 cDNA were highly correlated (R > 0.95). Linear regression of Cy-3 versus Cy-5 scatter plots yielded slopes from 0.94 to 1.02 at 55 or 60°C (Table 2). When 5 μg of Cy-3 was competitively hybridized with 10 μg of Cy-5 labelled cDNA, slopes from 0.48 to 0.53 were recorded at 55 or 60°C, as expected (Table 2). Since log-transformed signal intensities remained highly correlated (R>0.85), this assessed the robustness of the recorded signals. Furthermore, the median intensity of fluorescent signals was marginally altered by temperature changes (not shown).
In contrast, signal intensities from E. coli control probes were highly disturbed by temperature changes or by altering the ratios of fluorescently-labelled S. aureus cDNA (Table 2). Furthermore, the median intensity of fluorescent signals decreased by >2 fold at 60 compared to 55°C (not shown). Dye-swap experiments yielded equivalent results (not shown).
Reproducibility of StaphChip for expression profiling
To evaluate the reproducibility of fluorescent signals (Fig. 5), eight independent hybridization experiments were performed under identical conditions using Cy-3 labelled N315 cDNA, derived from overnight cultures. Maximal relative errors [(average-min)/average] of fluorescent probe signal intensities were more widely scattered on E. coli (Fig. 5B) compared to S. aureus (Fig. 5A) capture elements, thus reflecting the instability of poorly specific interactions with S. aureus targets and E. coli probes. The 90th percentile of maximal relative errors from S. aureus probes represented <25% of average signal intensities. In contrast, the same percentile of maximal relative errors recorded from E. coli probes represented >100% of average signal intensities. Furthermore, the 50th percentile on E. coli probes was superior to the 90th percentile of maximal relative errors on S. aureus probes.
Finally, we compared signal intensities generated from individual gene transcripts covered by two or more adjacent but non-overlapping S. aureus probes. Figure 6 shows average fluorescence intensities and their maximal relative errors (A) and the cumulative distribution of maximal relative errors (B) For nearly 90% of the gene transcripts (n = 2,269), maximal relative errors of ORF-related signals were <60%. This yields evidence that multiple probes recognizing the same transcripts provide reproducible signals and that StaphChip provides reliable and robust determinations of genome-wide transcripts.
Validation of relative transcript levels with RT-qPCR
Figure 7 shows the fold changes estimated by quantitative RT-PCR for 18 transcripts found to be either up regulated, down regulated or unchanged on the StaphChip. A good fitting was observed between both methods, with a correlation coefficient of 0.91 obtained from linear fitting. Of note, the quantitative changes recorded with RT-qPCR tended to be higher than those quantified with the microarray. This finding is supported by previous observations that the dynamic range of RT qPCR is higher than that of microarray [22, 23], that only reach 2–3 orders of magnitude. These data validate the use of our method for quantitative gene-expression analysis.
Several recent studies have shown the usefulness of microarrays for genomic and global transcriptomic studies of microbial pathogens [23–25]. Each step of microarray experiments needs to be optimized and validated, from array design and manufacture to data collection and analysis. Among critical technical parameters that need to be controlled are microarray surface chemistry, probe sequence, probe deposition process, and hybridization conditions. Accuracy of microarray-generated data can be improved by using multiple replicates, dye swaps and statistical data analysis. Compared to cDNA microarrays, oligoarrays provide a flexible design and are considered more reliable in terms of sensitivity and specificity [26–29]. Reported drawbacks of cDNA or PCR probes include: i) unpredicted secondary structures, ii) uncontrolled cross-hybridization occurring on repeats or partially homologous regions of PCR probes, and iii) varying amounts and purity of spotted products [27, 30]. Recent software development allows genome-wide selection of sub-sequences that uniquely identify genes. Ideally, these approaches should amplify fragments of constant length, thus minimizing the differences in PCR amplification efficiency as well as in hybridization kinetics . However, the extent of cross-hybridization has rarely been evaluated and reported, and thus may lead to severe errors in higher level data analysis, such as clustering , genome composition analysis and genotyping for molecular epidemiology .
Most oligoarray applications dedicated to prokaryotes were developed by companies using proprietary algorithms [33, 34] whose detailed properties are rarely available. Furthermore, the lack of published validation data prevents adequate comparison of those short-probe oligoarrays with investigator-designed oligonucleotide arrays.
To date, several strategies for oligoarray design have been described, but their experimental validation is often limited  or absent . A drawback of these approaches is to apply thermodynamics laws on probe/target interaction as determined in solution , but ignoring effects related to oligonucleotide immobilization on a microarray surface [15, 38].
To address this issue, OliCheck considered the influence of predicted probe/target binding with respect to its position along the immobilized probe, as demonstrated by Hughes et al . This tool allowed selecting probes from large oligonucleotide libraries in order to provide multiple genome coverage, suitable for epidemiological and transcriptomic studies.
As a specific application, OliCheck was used to design StaphChip, an oligoarray dedicated to genomic studies on S. aureus, a clinically important pathogen with a low GC content and numerous sequence repeats. The 5,427 feature elements were selected to ideally cover all ORFs of three S. aureus genomes with two non-overlapping probes, as validated under experimental conditions. A particular achievement of this strategy is to yield cross-annotations between the designed probes and homologous ORF regions conserved across several genomes. Any new genome sequence can be screened by OliCheck to identify probes that can specifically detect homologous ORFs. Cross-annotations on the recently released S. aureus MW2 genome  yielded 78% gene coverage. It should be mentioned that the recently published S. aureus COL genome  composition and annotation have changed significantly since the early release by TIGR in 2003.
The properties of StaphChip design were confirmed by comparative genome hybridization and global gene expression studies. Work in progress assesses the reliability of StaphChip for monitoring S. aureus transcription changes during biofilm formation, endocytic stage, or expression of antibiotic resistance. Another achievement (unpublished data) was the design of oligonucleotide probes for the genomes of Neisseria meningitidis A and B, E. coli K12, Erwinia carotovora and E. chrysanthemi; having GC contents ranging from 32.8 to 52%. Furthermore, OliCheck design expands oligoarray use for the study of host-pathogen interactions by potentially preventing cross-hybridization between bacterial probes and contaminating host nucleic acids.
In summary, this work describes a validated approach to select optimal oligoarray capture elements for S. aureus expression analysis and comparative genome hybridization studies. This generic approach will enable researchers to develop customized oligoarrays for genomic studies of any sequenced microorganism.
Design of specific oligonucleotide probes
An initial set of candidate oligonucleotide probes was generated by ArrayDesigner™ 1.17 (Premier Biosoft Intl) using the following probe parameters: (i) 40–60 bp probe length, (ii) 65 ± 10°C target Tm, (iii) <-5.0 kcal/mol for hairpins, (iv) <-8.0 kcal/mol for self-dimers, and (v) dinucleotide repeats shorter than 5 bp (Supplementary material provides OliCheck input format description, for using other probe design software or probe list). The program tested separately each open reading frame (ORF) of the different S. aureus genomes freely available at NCBI (S. aureus N315 [Genbank# BA000018], Mu50 [Genbank#BA000017]), and TIGR (S. aureus COL [TIGR unfinished microbial genome, released in March 27, 2003]). A comprehensive list of all possible probes ranked according to thermodynamic criteria was provided for each genome (Fig. 1a).
Further selection of specific oligonucleotide probes was performed by the design of an original program called OliCheck. This approach is derived from the experimental findings of Hughes et al. that microarray hybridization signals are mostly influenced by mismatches in the solution-end (distal part) rather than surface-end (proximal part) portion of the oligonucleotide probe . OliCheck queries the locally available S. aureus genome databases by performing a BLAST (BLASTN for Windows, version 2.2.2)  analysis for each probe.
Each BLAST output is analyzed to extract alignment information. The best scored target-probe alignment (perfect match) is considered as the target of the probe, therefore exempted from further testing. To avoid cross-hybridization, all other alignments are tested using the specificity filtering test. For this test, any 60-mer probe that would align to the target with >11 matched nucleotides by its distal third (solution-end), or >31 matched nucleotides by its proximal two-thirds (surface-end) is rejected. An additional condition is the exclusion of irrelevant target-probe pairs with >18 consecutive nucleotide matches.
Further probe selection is achieved by aligning each candidate probe from each genome (e.g. N315) against the genomes of the other S. aureus strains (e.g. Mu50 and COL) (Fig. 1d, and 1e). This process allows annotating probes detecting homologous targets in the other genomes.
Each BLAST output is analyzed to extract alignment information. The best scored target-probe alignment is tested to predict a high hybridization signal. This efficiency test requires the absence of any mismatch in the distal half and <20 mismatches in the proximal half of the probe. If those requirements are not fulfilled, the probe is rejected; otherwise the alignment is considered appropriate for detecting a homologous target and the process is continued. To avoid cross-hybridization with targets from other genomes, further alignments obtained with that probe are checked by the specificity filtering test, as defined in step B. Each probe fulfilling these requirements is annotated as detecting a unique homologous sequence target.
Using a spreadsheet program, the probes were sorted by their ORF target names and their best thermodynamic criteria. Probes showing the best combination of thermodynamic criteria and strain coverage were selected for microarray manufacturing (Fig. 1f). In addition to the selected S. aureus probes, OliCheck was used with the same parameters to design 2,873 probes specific for Escherichia coli K12 (E. coli K12 [Genbank# U00096]).
A compiled version of OliCheck compatible with Windows 2000 or XP and written in the Delphi programming language (Delphi 7, Borland) is freely downloadable for non-profit use at Genomic Research Laboratory website .
In silico comparison of our algorithm with ArrayOligoSelector
Three probes set of 60-mer probes with homogenous thermodynamic criteria (Tm = 60°C) were generated using default parameters for N315 genome by: i) ArrayDesigner ii) ArrayOligoSelector for all candidate probes, iii) ArrayOligoSelector for the best candidate probes (one per ORF). The output lists generated by ArrayOligoSelector were reformatted to match OliCheck input file format. The list of probes generated by either software was further processed by OliCheck for cross-homology filtering. The sets of probes selected by each method were further compared for homology using BLAST. Alignment showing E-value <1E-20 were considered as homologous.
The StaphChip microarray was manufactured by in situ synthesis of 8,455 long oligonucleotide probes (Agilent). It consists of 5,427 S. aureus and 2,873 E. coli specific probes, together with A. thaliana control probes for spiked controls.
Preparation of the labelled nucleic acids
For comparative genome hybridization, each S. aureus strain was grown overnight in 2 ml Mueller-Hinton broth (MHB) and total DNA was extracted and purified using DNeasy columns (QIagen) following manufacturer's instructions. DNA purity and concentration were assayed by spectrophotometer. 2 μg DNA were labelled by the Klenow fragment of DNA polymerase I (BioPrime, Invitrogen) with Cyanine-3 or Cyanine-5 coupled dCTP (NEN) for 2 hours at 37°C, then stopped by the addition of 5 μl 0.5 M EDTA. Labelled DNA was purified on QiaQuick columns (Qiagen).
For gene expression analysis, total RNA was extracted from 2 ml exponential or overnight cultures using the Rneasy kit (Qiagen) as previously described . Batches of 5 or 10 μg of total S. aureus RNA were spiked with increasing amounts of different Arabidopsis thaliana mRNAs (SpotReport, Strategene), used as external calibrators. The RNA mixture was labelled by Cy-3 dCTP or Cy-5 dCTP, using the SuperScript II (Invitrogen) following manufacturer instructions. Labelled cDNA was then purified onto QiaQuick columns.
Hybridization and scanning parameters
Unless specified, equivalent amounts of cDNA (or genomic DNA) labelled with Cyanine-3 or Cyanine-5, were diluted in 250 μl Agilent hybridization buffer, and hybridized at a temperature of 60°C for 17 hours in a dedicated hybridization oven (Robbins Scientific). For comparative genome hybridization, genomic DNA from each individual S. aureus strain was labelled with Cy3 and co-hybridized with equivalent amounts of Cy5-labelled genomic DNA pooled from N315, Mu50 and COL . Slides were washed, dried under Nitrogen flow and scanned (Agilent) using 100% PMT power for both wavelengths. Data were extracted and processed using Feature Extraction™ software (version 5.0, Agilent).
For gene expression analysis, saturated spots were excluded from subsequent analysis. Local background-subtracted signals were corrected for unequal dye incorporation or unequal load of labelled product. The algorithm consisted of a rank consistency filter and a curve fit using the default LOWESS (locally weighted linear regression) method. Spots showing a reference signal lower than background plus two standard deviations were also excluded from subsequent analyses.
For comparative genome hybridization, local background-subtracted data were expressed as Log10 ratios and analyzed by two-way clustering using GeneSpring 6.1 (SiliconGenetics).
Real time, quantitative PCR analysis
The expression of 18 genes involved in metabolic pathways was quantitatively assayed using by 1 step RT-qPCR using a SDS7700 (Applied Biosystems, Framingham, MA). Primers and probes were identified by scanning each gene sequence using the software Primer Express 2.0(Applied Biosystems). All identified sequences were further aligned on the whole genome of sequenced strains to ensure gene specificity and conservation of the target sequence between strains. Optimal concentration of primers and Taqman probes (labelled with FAM in 3' and coupled to TAMRA in 5' as quencher and purchased from Eurogentec, Seraing, Belgium), determined accordingly to manufacturer's instructions were 200 nM and 100 nM per reaction, respectively. Primers and probes were mixed in Platinum qRT-PCR Thermoscript kit (Invitrogen) with 0.4 ng of total purified RNA. Fold changes were calculated after normalization with the expression level of the 16s rRNA gene as previously described .
Yershov G, Barsky V, Belgovskiy A, Kirillov E, Kreindlin E, Ivanov I, Parinov S, Guschin D, Drobishev A, Dubiley S, Mirzabekov A: DNA analysis and diagnostics on oligonucleotide microchips. Proc Natl Acad Sci U S A. 1996, 93: 4913-4918. 10.1073/pnas.93.10.4913.
Dubiley S, Kirillov E, Mirzabekov A: Polymorphism analysis and gene detection by minisequencing on an array of gel-immobilized primers. Nucleic Acids Res. 1999, 27: e19-10.1093/nar/27.18.e19.
Kim CC, Joyce EA, Chan K, Falkow S: Improved analytical methods for microarray-based genome-composition analysis. Genome Biol. 2002, 3: RESEARCH0065-10.1186/gb-2002-3-11-research0065.
Lucchini S, Thompson A, Hinton JC: Microarrays for microbiologists. Microbiology. 2001, 147: 1403-1414.
Eriksson S, Lucchini S, Thompson A, Rhen M, Hinton JC: Unravelling the biology of macrophage infection by gene expression profiling of intracellular Salmonella enterica. Mol Microbiol. 2003, 47: 103-118. 10.1046/j.1365-2958.2003.03313.x.
Staudinger BJ, Oberdoerster MA, Lewis PJ, Rosen H: mRNA expression profiles for Escherichia coli ingested by normal and phagocyte oxidase-deficient human neutrophils. J Clin Invest. 2002, 110: 1151-1163. 10.1172/JCI200215268.
Golub TR, Slonim DK, Tamayo P, Huard C, Gaasenbeek M, Mesirov JP, Coller H, Loh ML, Downing JR, Caligiuri MA, Bloomfield CD, Lander ES: Molecular classification of cancer: class discovery and class prediction by gene expression monitoring. Science. 1999, 286: 531-537. 10.1126/science.286.5439.531.
Array Designer 3.01. http://www premierbiosoft com/dnamicroarray/index html. 2004, [http://www.premierbiosoft.com/dnamicroarray/index.html]
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215: 403-410. 10.1006/jmbi.1990.9999.
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-38.
Rouillard JM, Herbert CJ, Zuker M: OligoArray: genome-scale oligonucleotide design for microarrays. Bioinformatics. 2002, 18: 486-487. 10.1093/bioinformatics/18.3.486.
Wright MA, Church GM: An open-source oligomicroarray standard for human and mouse. Nat Biotechnol. 2002, 20: 1082-1083. 10.1038/nbt1102-1082.
Bozdech Z, Zhu J, Joachimiak MP, Cohen FE, Pulliam B, DeRisi JL: Expression profiling of the schizont and trophozoite stages of Plasmodium falciparum with a long-oligonucleotide microarray. Genome Biol. 2003, 4: R9-10.1186/gb-2003-4-2-r9.
Nielsen HB, Wernersson R, Knudsen S: Design of oligonucleotides for microarrays and perspectives for design of multi-transcriptome arrays. Nucleic Acids Res. 2003, 31: 3491-3496. 10.1093/nar/gkg622.
Hughes TR, Mao M, Jones AR, Burchard J, Marton MJ, Shannon KW, Ziman M, Meyer MR, Kobayashi S, 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: 342-347. 10.1038/86730.
Kuroda M, Ohta T, Uchiyama I, Baba T, Yuzawa H, Kobayashi I, Cui L, Oguchi A, Aoki K, Nagai Y, Lian J, Ito T, Kanamori M, Matsumaru H, Maruyama A, Murakami H, Hosoyama A, Mizutani-Ui Y, Takahashi NK, Sawano T, Inoue R, Kaito C, Sekimizu K, Hirakawa H, Kuhara S, Goto S, Yabuzaki J, Kanehisa M, Yamashita A, Oshima K, Furuya K, Yoshino C, Shiba T, Hattori M, Ogasawara N, Hayashi H, Hiramatsu K: Whole genome sequencing of meticillin-resistant Staphylococcus aureus. Lancet. 2001, 357: 1225-1240. 10.1016/S0140-6736(00)04403-2.
Baba T, Takeuchi F, Kuroda M, Yuzawa H, Aoki K, Oguchi A, Nagai Y, Iwama N, Asano K, Naimi T, Kuroda H, Cui L, Yamamoto K, Hiramatsu K: Genome and virulence determinants of high virulence community-acquired MRSA. Lancet. 2002, 359: 1819-1827. 10.1016/S0140-6736(02)08713-5.
Gill SR: Staphylococcus aureus COL genome. http://www tigr org/tigr-scripts/CMR2/GenomePage3 spl?database=gsa. 2004, [http://www.tigr.org/tigr-scripts/CMR2/GenomePage3.spl?database=gsa]
Blattner FR, Plunkett GIII, Bloch CA, Perna NT, Burland V, Riley M, Collado-Vides J, Glasner JD, Rode CK, Mayhew GF, Gregor J, Davis NW, Kirkpatrick HA, Goeden MA, Rose DJ, Mau B, Shao Y: The complete genome sequence of Escherichia coli K-12. Science. 1997, 277: 1453-1474. 10.1126/science.277.5331.1453.
Cramton SE, Gerke C, Schnell NF, Nichols WW, Gotz F: The intercellular adhesion (ica) locus is present in Staphylococcus aureus and is required for biofilm formation. Infect Immun. 1999, 67: 5427-5433.
Novick RP, Schlievert P, Ruzin A: Pathogenicity and resistance islands of staphylococci. Microbes Infect. 2001, 3: 585-594. 10.1016/S1286-4579(01)01414-9.
Ramakrishnan R, Dorris D, Lublinsky A, Nguyen A, Domanus M, Prokhorova A, Gieser L, Touma E, Lockner R, Tata M, Zhu X, Patterson M, Shippy R, Sendera TJ, Mazumder A: An assessment of Motorola CodeLink microarray performance for gene expression profiling applications. Nucleic Acids Res. 2002, 30: e30-10.1093/nar/30.7.e30.
Nuwaysir EF, Huang W, Albert TJ, Singh J, Nuwaysir K, Pitas A, Richmond T, Gorski T, Berg JP, Ballin J, McCormick M, Norton J, Pollock T, Sumwalt T, Butcher L, Porter D, Molla M, Hall C, Blattner F, Sussman MR, Wallace RL, Cerrina F, Green RD: Gene expression analysis using oligonucleotide arrays produced by maskless photolithography. Genome Res. 2002, 12: 1749-1755. 10.1101/gr.362402.
Wang D, Coscoy L, Zylberberg M, Avila PC, Boushey HA, Ganem D, DeRisi JL: Microarray-based detection and genotyping of viral pathogens. Proc Natl Acad Sci U S A. 2002, 99: 15687-15692. 10.1073/pnas.242579699.
Okamoto T, Suzuki T, Yamamoto N: Microarray fabrication with covalent attachment of DNA using bubble jet technology. Nat Biotechnol. 2000, 18: 438-441. 10.1038/74507.
Blanchard AP, Friend SH: Cheap DNA arrays-it's not all smoke and mirrors. Nat Biotechnol. 1999, 17: 953-10.1038/13644.
Barczak A, Rodriguez MW, Hanspers K, Koth LL, Tai YC, Bolstad BM, Speed TP, Erle DJ: Spotted long oligonucleotide arrays for human gene expression analysis. Genome Res. 2003, 13: 1775-1785. 10.1101/gr.1048803.
Kothapalli R, Yoder SJ, Mane S, Loughran TPJ: Microarray results: how accurate are they?. BMC Bioinformatics. 2002, 3: 22-10.1186/1471-2105-3-22.
Li J, Pankratz M, Johnson JA: Differential gene expression patterns revealed by oligonucleotide versus long cDNA arrays. Toxicol Sci. 2002, 69: 383-390. 10.1093/toxsci/69.2.383.
Yang IV, Chen E, Hasseman JP, Liang W, Frank BC, Wang S, Sharov V, Saeed AI, White J, Li J, Lee NH, Yeatman TJ, Quackenbush J: Within the fold: assessing differential expression measures and reproducibility in microarray assays. Genome Biol. 2002, 3: research0062-
Haas SA, Hild M, Wright AP, Hain T, Talibi D, Vingron M: Genome-scale design of PCR primers and long oligomers for DNA microarrays. Nucleic Acids Res. 2003, 31: 5576-5581. 10.1093/nar/gkg752.
Quackenbush J: Microarray data normalization and transformation. Nat Genet. 2002, 32 Suppl: 496-501. 10.1038/ng1032.
Rosenow C, Saxena RM, Durst M, Gingeras TR: Prokaryotic RNA preparation methods useful for high density array analysis: comparison of two approaches. Nucleic Acids Research. 2001, 29: U35-U42. 10.1093/nar/29.22.e112.
Dunman PM, Murphy E, Hanney S, Palacio D, Tucker-Kellogg G, Wu S, Brown EL, Zagurski RJ, Shlaes D, Projan SJ: Transcription Profiling-Based Identification of Staphylococcus aureus Genes Regulated by the agr and/or sarA Loci. J Bacteriol. 2001, 183: 7341-7353. 10.1128/JB.183.24.7341-7353.2001.
Tolstrup N, Nielsen PS, Kolberg JG, Frankel AM, Vissing H, Kauppinen S: OligoDesign: Optimal design of LNA (locked nucleic acid) oligonucleotide capture probes for gene expression profiling. Nucleic Acids Res. 2003, 31: 3758-3762. 10.1093/nar/gkg580.
Rouillard JM, Zuker M, Gulari E: OligoArray 2.0: design of oligonucleotide probes for DNA microarrays using a thermodynamic approach. Nucleic Acids Res. 2003, 31: 3057-3062. 10.1093/nar/gkg426.
SantaLucia JJ: A unified view of polymer, dumbbell, and oligonucleotide DNA nearest-neighbor thermodynamics. Proc Natl Acad Sci U S A. 1998, 95: 1460-1465. 10.1073/pnas.95.4.1460.
Ketomaki K, Hakala H, Kuronen O, Lonnberg H: Hybridization properties of support-bound oligonucleotides: the effect of the site of immobilization on the stability and selectivity of duplex formation. Bioconjug Chem. 2003, 14: 811-816. 10.1021/bc0340058.
Gill SR, Fouts DE, Archer GL, Mongodin EF, Deboy RT, Ravel J, Paulsen IT, Kolonay JF, Brinkac L, Beanan M, Dodson RJ, Daugherty SC, Madupu R, Angiuoli SV, Durkin AS, Haft DH, Vamathevan J, Khouri H, Utterback T, Lee C, Dimitrov G, Jiang L, Qin H, Weidman J, Tran K, Kang K, Hance IR, Nelson KE, Fraser CM: Insights on evolution of virulence and resistance from the complete genome analysis of an early methicillin-resistant Staphylococcus aureus strain and a biofilm-producing methicillin-resistant Staphylococcus epidermidis strain. J Bacteriol. 2005, 187: 2426-2438. 10.1128/JB.187.7.2426-2438.2005.
Genomic Research Laboratory. http://www genomic ch/software php. 2004, [http://www.genomic.ch/software.php]
Vaudaux P, Francois P, Bisognano C, Kelley WL, Lew DP, Schrenzel J, Proctor RA, McNamara PJ, Peters G, von Eiff C: Increased expression of clumping factor and fibronectin-binding proteins by hemB mutants of Staphylococcus aureus expressing small colony variant phenotypes. Infect Immun. 2002, 70: 5428-5437. 10.1128/IAI.70.10.5428-5437.2002.
Talaat AM, Howard ST, Hale W, Lyons R, Garner H, Johnston SA: Genomic DNA standards for gene expression profiling in Mycobacterium tuberculosis. Nucleic Acids Res. 2002, 30: e104-10.1093/nar/gnf103.
This work was supported by grants from the Swiss National Science Foundation PP00B–103002/1 (JS), 4049-63250 and 3200BO-103951 (PV), National Research Program 49 4049-63250 (PV and JS) and the Fondation pour Recherches Médicales (Gene Expression in Health and Diseases Program). We thank Prof K. Hiramatsu and J. Etienne for sending strains N315, Mu50, and MW2 as well as Prof F. Götz and P. Moreillon for providing strains SA113, SA113ica, and COL. We also thank Prof D. Lew for his constructive discussion throughout the work.
YC performed Olicheck implementation, designed microarray and wrote the manuscript. BG performed Olicheck implementation, designed microarray and helped writing the manuscript. PF participated in the design of the study, designed experiment protocols, performed microarray analysis. MB contributed to design experiment protocols and performed microarray experiments. AR has been involved in the CGH project. PV has been involved in drafting the article and revising it critically for intellectual content. WS has made substantial contributions to conception and design. JS initiated the study and helped writing the manuscript. All authors read and approved the final manuscript.
Yvan Charbonnier, Brian Gettler contributed equally to this work.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Charbonnier, Y., Gettler, B., François, P. et al. A generic approach for the design of whole-genome oligoarrays, validated for genomotyping, deletion mapping and gene expression analysis on Staphylococcus aureus. BMC Genomics 6, 95 (2005) doi:10.1186/1471-2164-6-95
- Comparative Genome Hybridization
- Oligonucleotide Probe
- Maximal Relative Error
- Candidate Probe
- Thermodynamic Criterion