Previous Article | Next Article ![]()
Applied and Environmental Microbiology, July 2004, p. 4216-4221, Vol. 70, No. 7
0099-2240/04/$08.00+0 DOI: 10.1128/AEM.70.7.4216-4221.2004
Copyright © 2004, American Society for Microbiology. All Rights Reserved.
Department of Veterinary Microbiology and Pathology,1 Washington Animal Disease Diagnostic Laboratory,2 Department of Civil and Environmental Engineering,3 WSU & UI Center for Reproductive Biology, Washington State University, Pullman, Washington4
Received 15 December 2003/ Accepted 31 March 2004
|
|
|---|
|
|
|---|
PCR-based assays are designed to amplify specific fragments of DNA that are subsequently identified by size using gel electrophoresis (3, 9, 13). By combining primer sets for multiple species-specific markers, a single PCR can be used to simultaneously detect several microbes. There are, however, practical limits to these assays. It can be difficult to incorporate more than six primer sets because of challenges in optimizing reaction conditions and challenges inherent in size discrimination between PCR products by conventional electrophoresis. Additional methods, such as Southern blotting or sequencing, are often needed to confirm product identity. An alternative approach is to use species-specific 16S ribosomal DNA (rDNA) PCR (17, 26) or universal PCR based on the 16S rDNA gene. Products from a universal reaction can be screened to identify species-specific polymorphisms in the DNA sequence using oligonucleotide probes (15), nested PCR (16), or automated sequencing. Sequencing is practical only if one 16S rDNA sequence is present in a sample. If more than one 16S rDNA sequence is present in a sample, then DNA microarrays can be used to efficiently screen a complex mixture of different sequences (27).
DNA microarrays permit simultaneous product interrogation with a large number of probe sequences (5, 7, 8, 27). In this format, probes for specific targets are typically deposited on a glass substrate to which PCR product or genomic DNA is then hybridized and detected (4). Because product size is no longer a factor in product identification, PCR assays can be designed to amplify similar-sized products and thus reduce PCR template bias. Furthermore, PCR assays can be designed to generate small products and thus maximize PCR efficiency and maximize the probability of template recovery from degraded DNA. Differences as small as 1-nucleotide base can be distinguished with an oligonucleotide-based microarray, although this degree of specificity is dependent on the sequence context (e.g., local melting temperature), hybridization conditions, and detection chemistry. In this study, we demonstrate a microbial detection system that couples PCR amplification of the 16S rDNA gene with an oligonucleotide-based microarray to simultaneously screen for 18 bacteria, 15 of which are fish pathogens.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Oligonucleotides used in this study
|
DNA was prepared by one of three methods. Cells were washed from agar plates using 1x phosphate-buffered saline (137 mM NaCl, 2.7 mM KCl, 10 mM Na2HPO4 [dibasic anhydrous], 2 mM KH2PO4 [monobasic anhydrous] [pH 7.4]), and DNA was extracted using a DNeasy kit (Qiagen, Valencia, Calif.). Alternatively, a single isolate was picked, put into fresh broth, and grown to stationary phase in the appropriate medium, and DNA was extracted using a DNeasy kit. DNA template was also prepared as boiled lysate when a colony was picked, resuspended in 200 µl of TE buffer (10 mM Tris base [pH 8], 1 mM EDTA [pH 8]), and boiled for 5 min. Template was removed after centrifuging the lysed cells for 5 min at 14,000 x g.
P. salmonis is an obligate intracellular pathogen and consequently was cultivated in CHSE-214 cells without antibiotics as described in the Office International des Epizooties (22). The CHSE-214 cells were harvested after a cytopathic effect was observed, approximately 10 to 14 days postinoculation. Harvested cells were extracted using the DNeasy kit, and species identification was confirmed using a PCR assay (19). The same kit was used to extract DNA from Chinook salmon (Oncorhynchus tshawytscha) kidneys obtained from a Pacific Northwest hatchery.
Our choice of such a diverse set of organisms reflects our interest in applying this technique across a phylogenetically disparate group of organisms. This selection also reflects broad regulatory concerns and the interests of clients that regularly use the Aquatic Health Inspection Service at the Washington Animal Disease Diagnostic Laboratory. E. coli, S. aureus, and P. phosphoreum were included as control organisms for this and other ongoing research projects (D. R. Call, A. Warsen, and M. Soule, unpublished data).
16S rDNA probe sequences.
16S rDNA sequences were retrieved from GenBank (www.ncbi.nlm.nih.gov). In some cases, the GenBank data were ambiguous; therefore, we generated new sequences for this project (F. columnare, F. psychrophilum, S. iniae, V. salmoninarum, and P. phosphoreum) using 16S universal primers 16s_008fwd and 16S_517rvs (Table 1). The resulting ca. 513- to 534-bp PCR products were subsequently cloned (TA kit; Invitrogen Corp., Carlsbad, Calif.) and sequenced (GenBank accession numbers AY577821 to AY577825) (Amplicon Express, Pullman, Wash.). Sequence alignments were used to identify polymorphic regions (Vector NTI version 7.1; InforMax, Inc., Bethesda, Md.) in the vicinity of position 453 (E. coli) (GenBank accession number AE000452). Oligonucleotide probes were designed for the polymorphic region using Primer Premier software (Biosoft International, Palo Alto, Calif.) (Table 1). All oligonucleotides were normalized to have a calculated annealing temperature of 63°C ± 3°C. Unmodified, desalted oligonucleotide probes were commercially synthesized (Invitrogen Corp.).
Slide preparation.
Multiple microarrays were printed on glass slides so that independent microarrays were contained within eight individual wells defined by Teflon masking (Erie Scientific, Portsmouth, N.H.); the hydrophobic nature of the masking permitted independent samples to be hybridized within each well. Slides were derivatized with an epoxy-silane monolayer (6). Prior to printing, the slides were sonicated in 2.5% Contrad 70 detergent (Fisher Scientific, Pittsburgh, Pa.) for 2 min, rinsed three times with distilled water, and dried using compressed air. Slides were then soaked for 1 h in 3 N HCl, rinsed three times with deionized water, and dried with compressed air. Slides were derivatized with a 2% solution of 3-glycidoxypropyltrimethoxysilane (Sigma Aldrich, St. Louis, Mo.) in methanol (high-performance liquid chromatography grade) for 15 min, rinsed twice with 100% methanol, and dried with compressed air.
Microarray construction.
Probes were reconstituted in TE buffer, quantified using a biophotometer (Eppendorf Scientific, Westburg, N.Y.), diluted to 60 µM in print buffer (0.1 M Na2HPO4, 0.2 M NaCl, 0.01% sodium dodecyl sulfate) with a pH of
11 and transferred to 96-well plates for printing. Arbitrary biotinylated oligonucleotides (25-mer; 5 µM) were included with every array. These biotin pseudoprobes served as positive controls for the detection chemistry and to orient the array for image processing. All probes were deposited as four replicates at a fixed location within each masked well using a Microgrid II arrayer (BioRobotics, Woburn, Mass.) with humidity held at 45%. Printing parameters included washing the pins in a recirculating bath (four pins washed twice for 4 s each time), followed by 0.5 s of flushing and 6 s of drying. This washing procedure was repeated twice between probes to minimize possible probe carryover. Printed slides were baked under a vacuum (22 Hg/mm) for 1 h (130°C) and stored away from light at room temperature.
Target preparation.
Biotinylated primers 16s_336 (this study) and 16s_517 (21) were used to amplify a ca.199-bp product from the 16S rDNA gene in both gram-positive and gram-negative bacteria. Each 50-µl PCR mixture contained reaction buffer (Fisher Scientific), 0.2 mM concentration of each deoxynucleoside triphosphate, 2.0 mM MgCl2, 2 to 5 U of Taq polymerase, 0.4 µM concentration of each biotinylated primer, and ca. 20 to 50 ng of genomic DNA or 10 µl of boiled lysate. PCR cycling conditions included an initial 5-min incubation at 95°C, followed by 25 to 35 cycles of PCR, with 1 cycle consisting of denaturation (95°C for 30 s), annealing (62°C for 60 s), and extension (72°C for 60 s), with a final extension step at 72°C for 10 min. Agarose gel electrophoresis was used to confirm the presence of PCR products.
Hybridization.
Prior to hybridization, slides were preblocked at room temperature for 30 min in TNB buffer. TNB buffer consists of 100 mM Tris-HCl (pH 7.5), 150 mM NaCl, and 0.5% blocking reagent (Biotin Tyramide Signal Amplication [TSA] kit; PerkinElmer Life Sciences, Boston, Mass.). PCR product (21.4 µl) was combined with hybridization buffer (final concentration, 4x SSC [1x SSC is 15 mM NaCl plus 0.15 mM sodium citrate, pH 7.0]) and 5x Denhardt's solution [0.001% Ficoll, 0.001% polyvinylpyrrolidone, 0.001% bovine serum albumin] in a final volume of 100 µl. This mixture was denatured by heat treatment (95°C for 2 min), cooled, and held on ice until application to two independent wells on the slide. Blocking buffer was aspirated and immediately replaced with 45 µl of denatured target per well. Slides were placed in a humidified chamber (50-ml conical tube with paper moistened with hybridization buffer). Chambers were then submerged in a 55°C water bath for incubation overnight. After removal from hybridization chambers, targets were removed by aspiration. Throughout these procedures, it was critical that the surface of the slide did not dry out. Slides were washed three times for 1 min each time in TNT buffer (0.1 M Tris-HCl [pH 7.5], 0.15 M NaCl, 0.05% Tween 20). A streptavidin-horseradish peroxidase conjugate (TSA kit) diluted 1:100 in TNB buffer was then incubated in each well for 30 min, and then each well was washed with TNT buffer three times for 1 min each time. Fetal equine serum (10%) in 2x SSC was incubated in each well for 30 min to provide a protein surface for tyramide binding, and then each well was washed with TNT buffer three times for 1 min each time. Biotinyl tyramide (1:50 in amplification diluent) was incubated in each well for 10 min. Slides were washed and incubated with 2 µg of streptavidin Alexa Fluor 546 conjugate (Molecular Probes, Eugene, Oreg.) per ml prepared in a mixture consisting of 1x SSC and 5x Denhardt's solution (total volume, 45 µl) for 1 h at room temperature (shielded from ambient light). After incubation, slides were washed and dried with a Microarray high-speed centrifuge (ArrayIt, Sunnyvale, Calif.) and imaged with an arrayWoRxe scanner (Applied Precision, Issaquah, Wash.).
Assay sensitivity.
We determined the baseline sensitivity using serial dilution curves (10-fold) between 6.6 x 108 and 6.6 x 1017 g. We arbitrarily selected genomic DNA from Y. ruckeri for these experiments, and each dilution series was repeated three times. Dilution series were assessed using (i) Y. ruckeri DNA only, (ii) Y. ruckeri DNA plus 5 ng of S. aureus DNA in each dilution, or (iii) Y. ruckeri DNA plus 500 ng of fish DNA (O. tshawytscha) in each dilution. NCSS 2004 (Number Cruncher Statistical Systems, Kaysville, Utah) was used for statistical analysis.
Signal detection and array analysis.
SoftWoRx Tracker software (Applied Precision) was used to quantify median spot intensity, and data were managed with a custom relational database (Access; Microsoft Corp., Redmond, Wash.). All experiments were performed in triplicate. For most purposes, we were able to identify positive hybridization simply by visual inspection. We also developed a cluster algorithm to score microarray data and detect positive hybridization. For each hybridization experiment, probe intensity data were sorted from the lowest signal intensity (ca. 200 U) to the highest signal intensity (65,535 U). The lowest and highest intensity values were identified as cluster seeds. The next lowest intensity value was compared to the cluster seeds and assigned to the nearest cluster. Cluster seeds were then recalculated as the average of all cluster members. This process continued until all probes were assigned to either the low (CL1) or high (CH1) cluster. Because this procedure will produce a higher variance for the CL1 cluster, the procedure was repeated, except that cluster membership was assigned in the opposite order, from the highest intensity to the lowest intensity to calculate CL2 and CH2. Final cluster averages were calculated as CL = (CL1 + CL2)/2 and CH = (CH1 + CH2)/2. All probe intensity values were then reassigned to the low and high clusters on the basis of the proximity to CL and CH. Standard deviations were then calculated for the low and high clusters and if CL and CH were separated by at least 3 standard deviations, then probes belonging to the high cluster were classified as positive. This procedure was programmed using Visual Basic (Microsoft Corp.) and is available upon request.
|
|
|---|
![]() View larger version (67K): [in a new window] |
FIG. 1. Probe specificity. Positive-control hybridization signals for A. hydrophila (AEHY), A. salmonicida subsp. salmonicida (AESA), E. ictaluri (EDIC), E. coli (ESCO), F. branchiophilum (FLBR), F. columnare (FLCO), F. psychrophilum (FLPS), M. chelonae (MYCH), M. fortuitum (MYFO), M. marinum (MYMA), R. salmoninarum (RESA), P. salmonis (PISA), S. aureus (STAU), S. iniae (STIN), T. maritimum (TEMA), V. salmoninarum (VASA), P. phosphoreum (PHPH), and Y. ruckeri type I (YERU) are shown. The lower right row of spots in each panel (not shown for MYMA or MYCH) is a biotinylated oligonucleotide control for detection chemistry.
|
![]() View larger version (33K): [in a new window] |
FIG. 2. Representative examples of microarray hybridization from no-template control PCRs. Little to no background is evident with 30 cycles or less, but considerable background is evident with 35 cycles of PCR.
|
The sensitivity of the assay, when limited to 28 PCR cycles, was equivalent to <150 genomes for purified template DNA (Fig. 3). This conclusion was based on the magnitude of probe signal (>0 U; P < 0.05 by the t test) and our cluster algorithm, which indicated positive detection with this starting concentration of template DNA. Importantly, this level of sensitivity was not compromised by the addition of a competing template DNA (5 ng of S. aureus DNA added; ca. 1.1 x 106 genomes), and the addition of 500 ng of fish DNA did not affect assay sensitivity (Fig. 3). Hybridization signals from nontarget probes (excluding Y. ruckeri type I, S. aureus, and E. coli) were significantly lower for dilution curves mixed with fish DNA (3,733 U) than for dilution curves with purified DNA (9,789 U) and purified DNA mixed with competing bacterial DNA (10,022 U) (analysis of variance with Tukey-Kramer multiple-comparison test; P < 0.05).
![]() View larger version (14K): [in a new window] |
FIG. 3. Assay sensitivity. Data showing assay sensitivity for purified pathogen DNA alone () or for DNA with competing bacterial ( ) or fish ( ) DNA. Competing DNA does not impact assay sensitivity (<150 copies; P < 0.05) using this protocol.
|
|
|
|---|
Our assay clearly detects and discriminates between 18 bacterial species. It is, however, important to note a potential "polymorphism paradox" in detecting these genetic differences. That is, we are purposely selecting a region of the 16S rDNA gene that has a high degree of variation; therefore, this region may be prone to frequent mutations. A strain that should be identified using the microarray might be missed due to a random mutation within this region and thus produce a false negative. We do not know what the frequency of these events might be at a population level, but we have evidence that this might be the case for one strain of F. pyschophilum (Call et al., unpublished data). An additional challenge with this strategy is that two "distinct" species may share the same sequence identity within the probe region (10), although this was not an issue for the 18 species tested here.
Because we had to reduce the number of PCR cycles in this assay, the assay sensitivity (<150 genomes) was not as high as can be achieved with PCR (<100 genomes [5, 12]). We submit, however, that the most effective gains in assay sensitivity can be achieved "upstream" of the microarray detector itself. That is, processing larger sample volumes (water) or mass (tissue) combined with more efficient DNA extraction techniques may be the most effective means to improve overall assay sensitivity (23).
For the assay detection sensitivity to be 100%, all truly positive samples must be detected and differentiated with the array. Nevertheless, depending on the sample matrix and "upstream" activities, false-negative results can be an important problem (18). In the format described here, coextracted sample impurities can interfere with PCR amplification. Failed amplification would not be readily apparent if only the microarray were used with no additional controls. One strategy to test for failed amplification is to spike samples with a known 16S rDNA sequence concurrently with the PCR. The internal control sequence would be amplified by PCR along with templates from the sample, and all products would be differentiated with the microarray. Adding template DNA runs the risk of reduced sensitivity if the spiked template is preferentially amplified during the PCR (Call et al., unpublished data). Alternatively, when complex prokaryotic DNA mixtures are amplified (e.g., from a water sample), we routinely see an amplified product by gel electrophoresis even when no targets of interest hybridize to the microarray. In these cases, nonspecific amplification of the 16S rDNA gene serves as an internal control for false-negative results due to PCR inhibition. Amplification of an independent eukaryotic marker would be very useful if tissues are being tested, although this is not part of the current assay.
We developed a relatively simple assay for screening PCR products for species-specific sequence polymorphisms. Our multiwell format permits higher sample throughput than conventional microarrays where a separate slide is required for each test. We have also found that the multiwell format is relatively simple to use and requires less time to develop necessary technical skills compared with unmasked slides that require coverslips. Ultimately, detection methods such as the one described here will become even more useful when the probes are used in conjunction with microsphere beads (14). Suspension arrays of this nature can be processed more efficiently, and they permit beginning-to-end sample processing in a 96-well format without need for posthybridization image processing.
This work was supported in part by funding from the WSU & UI Aquaculture Initiative, WSU & UI Center for Reproductive Biology, and Agricultural Animal Health Program, Pullman, Wash.
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»