Previous Article | Next Article ![]()
Applied and Environmental Microbiology, March 2003, p. 1397-1407, Vol. 69, No. 3
0099-2240/03/$08.00+0 DOI: 10.1128/AEM.69.3.1397-1407.2003
Copyright © 2003, American Society for Microbiology. All Rights Reserved.
Department of Molecular Ecology, Max Planck Institute for Marine Microbiology, 28359 Bremen, Germany
Received 13 September 2002/ Accepted 10 December 2002
|
|
|---|
|
|
|---|
A drawback of the FISH method is the limited number of probes which can be applied in one hybridization experiment. This limitation becomes a distinct bottleneck when this technique is used for community analysis on a high level of phylogenetic resolution. Additionally, there is a need for multiple probes to check for false-positive and false-negative results caused by individual probes used for identification of selected target organisms (2). A highly parallel application of multiple probe sets is facilitated by DNA chips or DNA microarrays. In this reverse hybridization format, matrix-immobilized oligo- or polynucleotides (probes) are used for the specific capture of labeled target molecules. This, in principle, allows the simultaneous application of a nearly unlimited number of probes in a single hybridization experiment (23, 26).
The concept of DNA microarray hybridization was introduced more than a decade ago (5, 14, 29, 31, 45), and the technique has been applied for many years in the pharmaceutical industry (12), clinical diagnostics (25), and many fields of research, e.g., functional genomics (28) and genetic analysis (22). In contrast, DNA microarrays are still not common in microbial ecology. Only a limited number of reports have been published, which mainly showed the "proof of principle" of the method in this field of research (24, 41, 43). What are the reasons? In contrast to "standard" applications in the medical field, in microbial ecology defined nucleic acids have to be identified against an often large and partly unknown genetic background. For this kind of analysis, the discrimination of single mismatches is crucial (27), but the specific hybridization of target molecules to immobilized capture oligonucleotides is one of the major challenges of the DNA microarray approach because large sets of probes with different characteristics are applied under identical hybridization conditions.
Several approaches exist to overcome this problem. An example is the addition of tetramethylammonium chloride (32) or betaine (39) to the hybridization buffer. Both compounds equalize the melting points of different oligonucleotides by stabilization of AT base pairs. Another strategy is based on the acquisition of complete melting curves for every single probe spotted on the chip surface (27), which is a promising but methodically and technically demanding approach. Further methodological investigations and optimizations are needed before DNA microarrays can be used for the fast, simple, and efficient screening of microbial communities. Currently, the usefulness of standard microarray formats is often limited by hard to interpret signal patterns caused by an accumulated number of false-positive signals.
In 1994, Williams et al. (49) pointed out the complexity of a microarray hybridization event and the shortcoming of the existing rules to describe it. In recent years, several systematic studies were conducted on different aspects, such as the secondary structures of target molecules and steric hindrance mediated by the solid support (8, 22, 33, 40, 42, 49), but the complex interaction of all these mechanisms and their influence on hybridization specificity are still mostly unknown. In this study, we systematically evaluated the suitability of the microarray format for 16S rRNA-based detection of bacteria by analyzing a model system including single-mismatch controls. Our objective was to get new insights into the mechanisms underlying a microarray hybridization reaction for the deduction of novel optimization strategies leading to signal patterns with a reduced number of false-positive and false-negative results.
|
|
|---|
Preparation of fluorescently labeled target single-stranded DNA.
Labeled target molecules were prepared by PCR amplification of the nearly complete 16S rRNA genes of the reference strains with the general bacterial primers 8F and 1492R (11). Living cells were picked from culture plates, resuspended in 65 µl of PCR-water (Sigma, Deisenhofen, Germany) and treated three times by freeze-thawing in liquid nitrogen for cell lysis. For PCR amplification, 10x PCR buffer (Eppendorf, Hamburg, Germany), 30 µg of bovine serum albumin, deoxyribonucleoside triphosphates (final concentration, 250 µM each), primers (final concentration, 150 nM each), and 2.5 U of MasterTaq polymerase (Eppendorf, Hamburg, Germany) were added to a final volume of 100 µl.
For subsequent preparation of fluorescently labeled, single-stranded DNA targets, 5'-indocarbocyanine-labeled forward primer 8f and 5'-biotin-labeled reverse primer 1492R were applied, both purchased from Thermo Hybaid (Interactiva Division, Ulm, Germany). Amplification was done according to the protocol of Buchholz-Cleven et al. (11) except that a temperature of 48°C was used for annealing. Successful amplification was confirmed by analyzing PCR products on 1.5% agarose gels. PCR products were purified with the QIAquick PCR purification kit (Qiagen, Hilden, Germany). Subsequently, the biotin-labeled DNA strands were extracted with streptavidin-coated paramagnetic beads (Boehringer, Mannheim, Germany) as described by Niemeyer et al. (36). The remaining, fluorescently labeled single-stranded DNAs are equivalent to the 16S rRNAs of the corresponding strains. Finally, the concentration of the amplicons was determined by UV spectrometry. The labeled single-stranded DNA was stored at -18°C.
Oligonucleotide probe set.
Twenty oligonucleotide probes of 15 to 20 nucleotides in length, targeting the bacterial 16S rRNA and originally designed for FISH, were selected from the literature for reliable differentiation of the six model strains based on the multiple-probe approach (4). The probes and their characteristics are listed in Table 1. One of the probes (KT13) was designed within this study with the software package ARB (http://www.arb-home.de). Current specificity of all probes was evaluated by ARB with the rRNA database of the Technical University of Munich (http://www.arb-home.de, release 08/01). For all 20 probes, a control containing one single central mismatch based on a transversion was included.
|
View this table: [in a new window] |
TABLE 1. Oligonucleotide probes used in this study and their specificitiesa
|
Probes were spotted onto the activated slide surface with the piezo-driven spotting device Robodrop (BIAS, Bremen, Germany). The concentration of the amino-modified oligonucleotides in PCR-water was 10 µM with 1% glycerol. Volumes of deposited probe solutions were about 250 pl, resulting in spots with a diameter of approximately 200 µm.
To complete covalent binding, after being spotted with probe solutions, slides were incubated overnight at room temperature in a wet chamber to restrict evaporation of the spots. Blocking of the slides was performed in 6-amino-1-hexanol (50 mM) and diisopropylethylamine (150 mM) in dimethylformamide as described by Beier et al. (6). Finally, the slides were washed with deionized, particle-free water, air dried, and stored under nitrogen at 4°C.
DNA microarray hybridization.
For hybridization and washing of the microarrays, a standard FISH protocol according to Pernthaler et al. (38) was used. Hybridization time was extended from 1.5 to 3 h. Incubation was performed either at room temperature, at 46°C without formamide, or at 46°C with 20% formamide. In all cases, 3.5 pmol of labeled target single-stranded DNA in a total volume of 200 µl of hybridization buffer was applied to the microarrays. To guarantee a uniform moistening of the slide surface, the sample was covered with a coverslip. Unless otherwise specified, probes were applied for hybridization experiments in the standard way, i.e., all probes immobilized via their 5' end, applied without spacer, and at the same concentration (10 µM).
To investigate the impact of secondary structures, unlabeled "helper" oligonucleotides (18), purchased from Thermo Hybaid (Ulm, Germany), were applied to the hybridization solution at a final concentration of 5 nM. For all probe binding sites analyzed by helper oligonucleotides, a pair of two perfectly matching 20-mer oligonucleotides was applied, one binding adjacent to the 5' end of the corresponding probe and one to the 3' end.
For reducing steric hindrance, capture probes carrying polyadenosine triphosphate spacers of 6, 12, 18, or 24 nucleotides, located at either the 5' end or the 3' end of the capture oligonucleotide, depending on selected probe orientation, were used.
Signal detection and data analysis.
Air-dried slides were imaged at a resolution of 10 µm with a GenePix4000 microarray scanner (Axon, Union City, Calif.) at the same laser power and sensitivity level of the photomultiplier for each slide. Therefore, absolute signal intensities (arbitrary units) presented for independent experiments should be directly comparable.
For automatic spot detection and signal quantification, the image analysis software MetaMorph Offline version 4.6 (Universal Imaging Corporation, West Chester, Pa.) was used. Filters for the automatic removal of spots with poor circularity or low uniformity of pixel intensities were established. Signals were considered positive if values were above zero after local background correction. This strategy for assessing the reliability of positive spots was possible because nonspecific binding of labeled target molecules to the slide surface was not observed for regions where oligonucleotides were deposited. This leads to "black holes" for spots where no detectable hybridization occurred, as shown in Fig. 6, and ultimately results in negative measurements for these spots, which were considered zero. We foresee no biases on individual spots by this method because in general we had homogeneous background noise. Otherwise, experiments were not considered.
![]() View larger version (59K): [in a new window] |
FIG. 6. Image of a scanned microarray taken from the experiment shown in Fig. 3A. Only four of the eight replicates spotted for each probe are shown (arranged in blocks). Columns 1 and 3 represent the 20 perfect-match probes, columns 2 and 4 show the corresponding mismatch controls. Highlighted spots are true-positive (solid line), false-positive (dashed line), or false-negative (dotted line) results.
|
|
|
|---|
![]() View larger version (21K): [in a new window] |
FIG. 1. Signal patterns after microarray hybridization of amplified 16S rDNA of strain KT0202a to the complete set of probes shown in Table 1 at room temperature without formamide (A), at 46°C without formamide (B), and at 46°C with 20% formamide (C). Only probes targeting strain KT0202a are shown. Black bars indicate false-positive signals (controls are only shown if signals were detected).
|
![]() View larger version (32K): [in a new window] |
FIG. 2. Comprehensive overview of signal patterns after microarray hybridization of all six strains analyzed in this study under standard conditions (at 46°C without formamide, all probes 5' immobilized and applied without spacer). False-positive signals are indicated by black bars; false-negative signals are indicated by open rectangles. For clarity, the single-mismatch controls spotted for all 20 probes are not shown (here, no signals could be observed).
|
Analysis of parameters influencing a microarray hybridization reaction.
Strain KT0202a was used for all further systematic evaluations because it was targeted by the most probes of the probe set applied (compare Table 1).
(i) Secondary structures.
With pairs of helper oligonucleotides, the signals for probes ALF968, ROS537, and RC1239 could be specifically increased (Fig. 3A; see also Fig 6). Compared to the corresponding data generated without helper oligonucleotides (Fig. 1B), the probe ALF968 signal was raised from 0 to 1,126 a.u. The signal for probe ROS537 was increased 16-fold, from 219 to 3,511 a.u., whereas the signal of probe RC1239 only doubled from 775 to 1,670 a.u. Interestingly, application of helper pairs also resulted in signals within the mismatch controls ALF968M1c and ROS537M1c. Also probe KT13, which has the same binding site as probe RC1239 and shows only one central mismatch to strain KT0202a, could no longer be discriminated under these conditions. Signals of probes MALF-1, GRb, and RC1031 decreased from 819 to 160 a.u., from 1,982 to 394 a.u., and from 973 to 197 a.u., respectively, compared to the corresponding experiments without helpers (Fig. 1B).
![]() View larger version (25K): [in a new window] |
FIG. 3. Impact of secondary structures of target molecules on hybridization efficiency. Hybridizations were performed with amplified 16S rDNA of strain KT0202a under standard conditions. (A) For probes ALF968, ROS537, and RC1239, helper oligonucleotides were added. Controls are only shown if signals were detected. (B) For all 11 probes targeting strain KT0202a, helpers were added. All corresponding controls plus two additional probes also featuring false-positive signals are shown.
|
(ii) General steric hindrance.
Poly(A) spacers of different lengths (none, 6-mer, and 12-mer) were tested on probes EUB338 and RRP1088. Both probes showed a linear increase of signal intensities with increasing spacer length (Fig. 4A). For probe EUB338, we observed signals within the corresponding mismatch controls when applying a 6- or 12-mer spacer, whereas for probe RRP1088, the hybridization efficiency could be improved by spacer molecules without getting nonspecific signals.
![]() View larger version (35K): [in a new window] |
FIG. 4. Impact of steric hindrance. (A) Hybridization of the amplified 16S rDNA of strain KT0202a at 46°C without formamide to probes EUB338 and RRP1088, both provided with poly(A) spacers of multiple lengths (none, 6-mer, and 12-mer spacers). The corresponding controls are also shown. (B) Hybridization of strain KT0202a at 46°C without formamide to probes ROS537 and ALF968 in the presence of helper pairs for both probes. Probes were provided with spacers of multiple lengths (none, 6-mer, 12-mer, 18-mer, and 24-mer spacers). Probe variants without a spacer and with a 12-mer spacer were also provided by immobilization via the 3' end of the probe instead of the 5' end.
|
(iii) Orientation of immobilized probe.
Probe orientation experiments with strain KT0202a (Fig. 4B) showed an increase in signal intensities after inversion (3' instead of 5' immobilized) from 3,716 to 6,763 a.u. (1.8-fold) and from 2,105 to 2,877 a.u. (1.4-fold) for probes ROS537 and ALF968, respectively. For the 12-mer spacer variants of the same two probes, increases of only about 1.1-fold were observed.
Optimization of signal patterns by directed application of capture oligonucleotides.
By a directed modification of spacer length based on the signal intensities shown in Fig. 1B for strain KT0202a, probe-specific differences in signal intensity and number of false-negative results should be minimized. For example, the two probes featuring former false-negative signals, UNIV1392 and ALF968, were both tested with 12- and 18-mer spacers, whereas probe GRb was not modified. The results for the hybridization of the 16S rDNA of strain KT0202a with this approach are given in Fig. 5. By elongation of probe UNIV1392 with an 18-mer spacer, the earlier observed false-negative signal could be restored without getting a signal for the corresponding mismatch control. In addition, signal intensities could be increased for probes MALF-1 (with a 6-mer spacer from 41 to 116%, normalized to the GRb signal), RRP1088 (12-mer, from 6 to 50%), PAR1457 (6-mer, from 9 to 25%), and RC1239 (6-mer, from 13 to 68%) without losing the specificity of the hybridization reaction. The sensitivity of probes EUB338, ALF968, ROS537, RSB67, and RC1031 could not be improved without the simultaneous appearance of signals within the corresponding controls.
![]() View larger version (14K): [in a new window] |
FIG. 5. Signal pattern optimization by directed application of capture oligonucleotides. Hybridization of the amplified 16S rDNA of strain KT0202a at 46°C without formamide to the 11 probes of the complete set targeting this strain. Probes were provided with spacers of different lengths, depending on former observed signal intensities (Fig. 1B, reference experiment). Only probe variants where no signals could be observed for the corresponding controls are shown.
|
|
|
|---|
The results shown in Fig. 2 indicate that highly specific hybridization of a large number of oligonucleotide probes is feasible in a standard microarray format without any effort to adjust the melting points of the probes. Only one false-positive signal was observed among 199 potential nonspecific hybridization events (0.5%). This particular mismatch was one of the two weakest analyzed in this experimental series, a G:T mismatch three nucleotides away from the 3' end of the probe (see Table 1). Discrimination of a mismatch near the terminus of a short duplex is hard to achieve (46) even under optimized conditions. The fact that we did not find a single nonspecific hybridization event within the central mismatch controls spotted for all 20 probes seems to be in complete contradiction with experiences from other hybridization formats such as FISH, where every probe has to be applied at appropriate stringency to ensure a specific hybridization reaction (4). Considering the results of the systematic investigations, the underlying mechanisms become more evident.
Parameter analysis.
The addition of helper oligonucleotides had a clear impact on the signal intensities of particular probes (Fig. 3). Helpers are unlabeled oligonucleotides designed to bind adjacent to a probe binding site, resulting in an increase in the accessibility of the corresponding target site (18, 35). As shown for probes ALF968, ROS537, and RC1239, they can resolve the secondary structures of the target rDNAs in a directed and selective way for particular probe binding sites (Fig. 3A and 6). Interestingly, the former strong signals of probes MALF-1, GRb, and RC1031 decreased dramatically in this experiment, suggesting that the opening of selected probe binding sites leads to a reorganization of secondary structures in other target regions. The influence of secondary structures on hybridization efficiency becomes much clearer when we applied helper pairs for all 11 probes targeting strain KT0202a (Fig. 3B): former false-negative signals could be restored, overall signal intensities were dramatically increased, and most of the single mismatch controls now showed strong signals. We do not assume that this is based on a change of Td values by base stacking of probe and helper oligonucleotides because no evidences for such an effect was found by Fuchs et al. (18).
Another outstanding feature of the data set shown in Fig. 2 is the distinct signal heterogeneity of different probes targeting the same preparation of 16S rDNA during parallel hybridization. This is a known phenomenon (44) and may even lead to the situation that signals of mismatch controls for particular probes were significantly higher than perfect match signals of other probes, as shown in Fig. 1A, 2, and 3. Considering the principles of DNA duplex stability prediction (9), measured signal intensities should be primarily connected to the nucleotide sequence of a probe. When applying helper oligonucleotides for all probes targeting an analyzed 16S rDNA, the distribution of measured signal intensities changes significantly (compare Fig. 1B and 3B). This observation is in agreement with the common view that secondary structures of target molecules affect the duplex yield of particular probes selectively (44). In conclusion, the measured signal intensity for a particular probe does not directly correlate with the amount of the corresponding target molecule in solution.
Another parameter that has to be considered in the context of signal limitation is steric hindrance. When hybridizing on a solid support, the binding efficiency of target molecules may be reduced by unfavorable steric interactions mediated by the solid matrix (42). The size of this effect should be correlated with the size of the target nucleic acid and the distance between the support and the capture probe. To investigate the influence of steric hindrance in our system, we applied poly(A) spacers of different lengths (6-, 12-, 18-, and 24-mers) between the slide surface and the probe sequence. The disadvantage of this kind of spacer is of course the potential Watson-Crick interaction between the spacer and the target molecule, but in our system, we did not observe any hybridization events indicating such an effect. Without the addition of helper oligonucleotides, in most cases we found a linear correlation between spacer length and measured signal intensity (Fig. 4A; supplementary data are not shown). When adding helper oligonucleotides, signal saturation was reached even with the shorter spacer variants (Fig. 4B, probe ROS537 [6-mer] and probe ALF968 [18-mer]), indicating a complex interaction of the signal-determining parameters. Interestingly, signal enhancement by reduction of steric hindrance may also lead to signals within the corresponding controls (Fig. 4A, probe EUB338).
Furthermore, it must be assumed that the position of the probe binding site also has an effect on the extent of steric hindrance because every nucleic acid binds in a defined orientation to its complement. To address the significance of this parameter, we tested two variants of probes immobilized via the 5' and 3' end. Probe ROS537, whose binding site is located for target KT0202a at nucleotide position 474 of 1428, starting from the 5' end, shows a higher signal for the 3'-immobilized probe. This is in good agreement with our initial assumption that a 5'-immobilized probe binding to the 5' end of the target molecule is leading to an unfavorable steric situation because the longer tail of the target nucleic acid is oriented towards the slide surface and vice versa. Nevertheless, probe ALF968 also shows a slightly increased signal for the 3'-immobilized probe, although its binding site for target KT0202a is located at nucleotide position 907 of 1428. Apparently, the secondary and tertiary structures of the target molecules interfere with simple spatial extrapolations. As for application of spacer molecules, we found little influence on signal intensities for reverse probe immobilization compared to the application of helper oligonucleotides. In conclusion, the secondary structure of a target molecule is the main limiting factor for duplex yield considering the parameters investigated in this study.
Mechanisms enabling observed hybridization specificity.
As shown in Fig. 3 and 4A, signals for the single-mismatch controls can be "generated" under identical hybridization conditions if the influence of signal-limiting parameters such as secondary structures and steric hindrance is reduced. Therefore, we conclude that for the data set shown in Fig. 2, mismatch discrimination was mainly achieved by global signal suppression. This is possible because perfect match and mismatch signals react different on changes of signal limiting parameters as shown for secondary structures: without helpers (Fig. 1B), one mismatch signals were selectively shifted below the detection limit compared to the perfect-match signals.
A full-length 16S rDNA molecule and the resulting secondary structure appear to represent a well-suited "calibrator" for global signal suppression under optimized hybridization conditions. When shorter fragments were hybridized under identical hybridization conditions (data not shown), we also observed an overall increase of signal intensities in combination with an accumulated number of false-positive signals. Here, the reduced influence of two factors becomes evident: (i) shorter fragments lead to reduced steric hindrance and (ii) secondary structures of shorter targets are known to be less stable than those of longer fragments (44). Therefore, our results suggest that fragmentation of target molecules can be disadvantageous in terms of preventing false-positive signals.
In summary, we report on a new concept for improving the specificity of microarray hybridizations. It is based on a converse approach compared to current signal enhancement attempts because we "utilize" signal-limiting parameters to selectively shift nonspecific signals below the detection limit. We would like to point out that the overall sensitivity of the detection system is not affected by this strategy because the detection limit is not shifted but hybridization signals are specifically suppressed.
Prevention of false-negative results by a new optimization approach.
As shown in Fig. 2, we found a large number of false-negative results in our data set. Subsequent application of helper pairs for all 11 probes targeting strain KT0202a under identical hybridization conditions led to signals for all 11 probes and nearly all corresponding controls (Fig. 3B). This indicates that false-negative results and discriminated single-mismatch controls were not based on melting effects. We suppose that nonoccurring signals are mainly caused by signal-suppressing parameters such as secondary structure of the target molecules and steric hindrance, selectively acting for different probe binding sites. Considering the fact that the number of 16S rRNA sequence positions is limited, a false-negative results for a particular probe can be a profound restriction for phylogenetic affiliation because additional probes featuring a redundant target specificity are often not available. For example, strain JP7.1 could be clearly assigned to the marine alpha group by the phylogenetically redundant probes MALF-1 and ROS537 (Fig. 2), but the single strain-specific probe available (Ros7-1029) featured a false-negative signal and therefore inhibited a more exact phylogenetic affiliation of the analyzed 16S rDNA.
Based on our systematic investigations we have deduced an approach for further signal pattern optimization called directed application of capture oligonucleotides. Our study shows that (i) single-mismatch discrimination can also be achieved by "proper calibration" of signal intensities for particular probes and (ii) different parameters are available for a probe-specific signal adjustment. As shown in Fig. 3A and 4A, signals can be selectively enhanced by reducing the influence of signal suppressing factors such as secondary structures and steric hindrance. Our intention was to restore former false-negative signals and to adjust overall signal intensities with these options without losing the status of complete discrimination for the single mismatch controls, as shown in Fig. 4A for probe RRP1088. The usability of this approach, based on a directed variation of spacer length, was successfully tested for strain KT0202a, as shown in Fig. 5. This parameter was chosen because it allows a graduated signal adjustment. Another potential mediator would be the amount of immobilized probe, because a direct correlation to duplex yield was also observed for this parameter (data not shown), in agreement with the results of Guo et al. (22). In the context of directed application, the utilization of helper oligonucleotides is not practical for two reasons: (i) the strong signal increase caused by helpers often also results in signals for the corresponding mismatch controls (Fig. 3) and (ii) design of helpers showing the same specificity as the corresponding probe is often not possible (18).
Data normalization.
An important issue when analyzing microarray experiments is data normalization. In the field of microbial ecology published data were either not normalized (24, 43) or domain specific probes were used for normalization (27). For this study we decided to work with absolute signal intensities (arbitrary units) when comparing data generated under different hybridization conditions or analyzing different target molecules. We assume that by using domain-specific probes as internal standards, an additional error will be introduced, because probe binding site accessibility should not be completely comparable for distantly related targets.
Why use a model system based on 16S rDNA targets?
This study was performed with 16S rDNA targets generated by PCR amplification. We are well aware of the possible biases introduced by PCR during analysis of complex communities (3). For the analysis of complex samples, direct extraction and labeling of the rRNA are essential to circumvent possible artifacts. Furthermore, for environmental studies, aspects such as sensitivity and quantification have to be considered more comprehensively. However, we decided to conduct this initial systematic investigation with exactly defined 16S rDNA targets instead of extracted RNA because only single targets were analyzed in our model system. Our results should help to consider new kinds of optimization strategies for the specific detection of 16S rRNA/rDNA molecules by microarray hybridization and for environmental studies. We foresee PCR-based microarray systems in the (i) fast and low-cost screening of isolates for their phylogenetic affiliation and (ii) prescreening of PCR products for diversity or the presence of particular subpopulations (e.g., for generating gene libraries).
This work was done in the framework of the Centre of Applied Gensensorik (CAG) at the University of Bremen and supported by the German Ministry of Education and Research (BMBF, contract 0311833A) and by the Max Planck Society.
|
|
|---|
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»