Bacterial discrimination by means of a universal array approach mediated by LDR (ligase detection reaction)
© Busti et al; licensee BioMed Central Ltd. 2002
Received: 24 June 2002
Accepted: 20 September 2002
Published: 20 September 2002
PCR amplification of bacterial 16S rRNA genes provides the most comprehensive and flexible means of sampling bacterial communities. Sequence analysis of these cloned fragments can provide a qualitative and quantitative insight of the microbial population under scrutiny although this approach is not suited to large-scale screenings. Other methods, such as denaturing gradient gel electrophoresis, heteroduplex or terminal restriction fragment analysis are rapid and therefore amenable to field-scale experiments. A very recent addition to these analytical tools is represented by microarray technology.
Here we present our results using a Universal DNA Microarray approach as an analytical tool for bacterial discrimination. The proposed procedure is based on the properties of the DNA ligation reaction and requires the design of two probes specific for each target sequence. One oligo carries a fluorescent label and the other a unique sequence (cZipCode or complementary ZipCode) which identifies a ligation product. Ligated fragments, obtained in presence of a proper template (a PCR amplified fragment of the 16s rRNA gene) contain either the fluorescent label or the unique sequence and therefore are addressed to the location on the microarray where the ZipCode sequence has been spotted. Such an array is therefore "Universal" being unrelated to a specific molecular analysis. Here we present the design of probes specific for some groups of bacteria and their application to bacterial diagnostics.
The combined use of selective probes, ligation reaction and the Universal Array approach yielded an analytical procedure with a good power of discrimination among bacteria.
The detection, identification, and characterization of bacterial populations is an important goal in analytical microbiology. Culture-independent techniques represent a rapid and flexible mean to study bacterial communities; in fact, the use of 16S rRNAs as molecular marker has became routine for microbial ecologists. The most comprehensive strategy to characterize bacterial populations probably consists in 16S rDNA clones sequencing and phylogenetic reconstruction . However, analysis of individual clones in multiple libraries is expensive and time consuming and therefore not suited to large-scale screenings. Other methods to assess the molecular composition of an environmental DNA sample, such as thermal or denaturing gradient gel electrophoresis (DGGE) , single stranded conformational polymorphism (SSCP)  heteroduplex analysis [4, 5], or terminal restriction fragment (T-RFLP or TRF) analysis [6–9], are more rapid and therefore amenable to large-scale experiments.
Moreover, the employment of group-specific DNA probes complementary to 16S rRNA has provided a framework to study microbial populations in complex systems. The recent development of the DNA microarray technology has added a high throughput experimental format, potentially with great sensitivity [10–12].
In the microarray format, the most commonly used procedure is the differential hybridization of a fluorescently labelled target, often a PCR product, with microarray-immobilized oligonucleotide probes. This method, in order to gain good probe specificity, requires very careful probe design and optimized experimental set up.
This approach presents some advantages. Ligase detection reaction had been shown to be a sensitive assay for detecting Single Nucleotide Polymorphisms , therefore a difference in a single nucleotide along the 16S rRNA can be employed to distinguish between sequences of different microorganisms. Moreover the system maintain the positive characteristics of the microarray format without requiring the optimization of the hybridization conditions for each probe set. Using such an approach we targeted the 16S rRNA genomic region using 223 sequences of cyanobacteria, 987 of actinomycetes, 284 of clostridia, 281 of bacilli, 69 of myxobacteria and 270 of pseudomonads, selected from the Ribosomal Database Project to identify consensus sequences for each group. Group-specific consensi were used to design selective molecular probes. These probes in a LDR on DNA from pure bacterial cultures, gave excellent selectivity for the target group and sensitivity down to 10 fmol of amplified 16S DNA.
Sequence analysis of 16S rDNA and ligation probes design
Secondly, among this set of probes, we selected only those pairs differing from all representatives of the other five groups at least for the 3' terminal position of the discriminating oligos, but invariant in all members of their group. This second criterion significantly reduced the number of actinomycete, clostridium and cyanobacteria-specific probe pairs, as shown in Table 1, column B.
Finally, in order to discard potentially aspecific probe pairs, we analyzed each common probe and discriminating oligo using the Probe Match tool on RDPII database, which permits to verify probes against all the bacteria sequences not considered in our alignments . This analysis significantly reduced the number of pseudomonads and myxobacteria-specific probe pairs. Furthermore, the identification of a clostridium-specific probe was not possible, while more than one was found for some of the other groups (Table 1, column C).
For the subsequent experimental work, we decided to select just one probe pair for each group of interest (Tab.2). When more than one base was present in the same position of the consensus, we included inosine, during oligo synthesis, at these degenerate positions.
In order to have a positive control for the Ligation Detection Reaction, a universal probe pair, matching all the studied groups, was designed according to the process described above, and the corresponding Zip code was included in the Universal Array.
Zip Codes assignment and quality control of the universal microarray
Selected group-specific probes
DISCRIMINATING OLIGO (5'–3')
COMMON PROBE (5'–3')a cZIP code
In order to verify the quality of deposition of the Zip Code oligos to the slides, we performed hybridizations with Cy5 labeled poly(dT) which is complementary to the poly(dA)10 sequence present in each Zip Code.
Ligation reaction detection set up
Using purified, PCR amplified 16S segment from Pseudomonas putida DNA as substrate, the Ligation Detection Reaction was set up changing probe concentration (100–400 fmol/μl), temperature (60–65°C) and number of cycles (20–60) (data not shown). Optimized working conditions are those described in Material and Methods. For each of the six groups, we amplified 16S rDNA according to the following scheme. When the group included a single genus, the template for the PCR amplification consisted of purified DNA from a single strain (Bacillus subtilis, Pseudomonas putida, Clostridium perfringens). When more than one genus was included within a group, PCR amplifications were conducted on a mixture consisting of DNA prepared from two (bacteria belonging to the genus Myxococcus), three (S. cinnabarinum, D. matsuzakiense, A. teichomyceticus) or four (Anabena, Nostoc, Microystis, Synechocystis) strains belonging to selected genera within mixobacteria, actinomycetes or cyanobacteria, respectively. LDRs were conducted in the presence of each group-specific PCR product as template and of all the probes (6 discriminating oligos and 6 common probes).
This result indicated that, in the absence of the perfectly matching PCR product, the probes present in the LDR mix do not generate false signals. Thus, the LDR reaction proceed in a template-specific manner.
Use of complex molecular targets
In order to determine the efficiency of the LDR technique, we carried out different assays varying the complexity of the molecular target. In details, we used artificial mixes of genomic DNA samples or mixes of PCR products.
Similar experiments were performed with balanced (1:1) and unbalanced (1:10) mixes of two out of five groups yielding the expected results (data not shown).
The main goal of this work was the development of a flexible method to detect bacterial groups. Ligation Detection Reaction, combined with a Universal Microarray, appeared an interesting approach suited for this application [13, 14]. It requires PCR amplification of a target region, in this case the bacterial 16S rDNA, which is then subjected to a multiplex cycled LDR. The LDR is achieved using Pfu DNA ligase, a thermostable enzyme which seals the nick between two adjacent oligonucleotides (the common probe and the discriminating oligo), annealed to a complementary target, only if the oligonucleotides are perfectly base-paired, in particular at the junction site (Fig. 1A). Therefore a single mismatch in 3' terminal position of the discriminating oligo is able to prevent ligation, thus conferring total selectivity . This feature confers an high resolution power to hybridisation, decreasing the effort for the search for stringency conditions. As shown in Fig. 1B, the presence of a specific target is determined by hybridizing the content of a LDR to an addressable DNA universal array, on which every single spot contains oligos with a unique Zip Code. A complementary cZip code is affixed to the 3' end of each common probe. During hybridization, the cZip Code drives the LDR product to the corresponding Zip Code on the chip surface. As every discriminating oligo carries a Cy3 molecule on its 5' end, detection of hybridized LDR products can be accomplished by laser scanning.
Probe design can be considered a crucial point: during the definition of subgroup and group-specific consensi, a cut off of 75% allowed preserving as much sequence information as possible, but required the inclusion of some degenerated positions in the probe sequences. Probes containing too many ambiguous residues were discarded, while a limited number of inosine residues was included in the oligos. Furthermore, we adopted a three-step selection process to ensure as much specificity as possible and this involved the rejection of about 80% of the probes identified after the first step. Due to these stringent criteria and ligation assay requirements, a low number of suitable group-specific probes was identified.
In fact, at this level of phylogenetic resolution, it was very difficult finding unique positions, with the ability to discriminate between groups, that fell in conserved region inside the group itself. In our experience, if the target groups are phylogenetically less distant, like members of the same order, probe design can be much more fruitful. In this case, in fact, a relevant part of diversity can be eliminated by the use of more specific PCR primer, instead of the "universal" primers F27-R1492. Moreover, in a more related subgroup we found that the 16S sequences are more similar thus simplifying probe construction (Castiglioni et al. unpublished results). This suggests that this approach seems particularly appealing for "fishing out" certain bacterial species within a complex microbial community
The combined use of selective probes and LDR gave satisfactory results. LDR combines the specificity of the hybridization base pairing with the selectivity introduced by the enzymatic reaction , resulting in good power of discrimination as demonstrated by the presented results. It should be emphasized that perfect pairing in the 3' terminus of the discriminating oligos and the 5' terminus of the common probe is crucial for ligation. On the contrary mismatches placed along the remaining part of the two sequences are easily tolerated by ligase, conferring a certain flexibility in probe design for test on complex samples.
As described we were able to detect the presence of different groups in balanced and unbalanced mixes (1:1 and 1:10 molar ratio respectively).
The optimized LDR method can be performed starting from low amounts of substrate. As little as 1 fmol of PCR-amplified material can be observed in our conditions. Below this limit not enough signal can be observed even increasing the amount of probes and the number of cycles (data not shown). Apart from any consideration regarding overall sensitivity and the feasibility of our procedure for quantitation, these results (sensitivity down to 1 fmol and proper results with unbalanced mixes) suggest the possibility of detecting a low amount of a specific 16S molecular fragment within complex 16S molecular mixtures.
We think this approach is particularly appealing for different reasons. First of all, since the ZipCodes sequences are not related with a specific molecular analysis, they remain constant and their complements can be appended to any set of LDR primers. In this sense, the array can be defined Universal. Moreover, the optimization of hybridization conditions for each probe set is not required, therefore the Universal chip become a versatile tool as new probe pairs can be added to the system without further optimization, thus reducing costs and set up time. Presented results suggest that a combination of careful probe design, PCR and LDR can be a valuable tool for the detection of bacterial groups in the environment although an intensive validation is required in order to ascertain potential interferences in complex natural samples.
All chemicals and solvents were purchased from Sigma-Aldrich (Italy) and used without further purification. Oligonucleotides were purchased from Interactiva Biotechnologie GmbH (Germany).
Genomic DNA from E. coli (ATCC 10536), Bacillus subtilis (ATCC 8185), Pseudomonas putida (ATCC 11250), Streptosporangium "cinnabarinum" (DSM 44094), Dactylosporangium matsuzakiense (ATCC 31570), Actinoplanes teichomyiceticus (ATCC 31121) and two myxobacteria belonging to genus Myxococcus (F. Gaspari, personal communication), was extracted using PUREGENE™ DNA isolation kit (Gentra Systems, Minneapolis, MN) according to manufacturers instructions. Clostridium perfringens DNA was purchased from Sigma-Aldrich (Italy). Genomic DNA from bacteria belonging to the generaAnabena, Nostoc, Microcystis and Synechocystis, and DNA extracted from water samples coming from European lakes were kindly supplied by Dr. Stefano Ventura (CNR-CSMA, Firenze, Italy).
Ligation probe design
The probes for Ligation Detection Reaction were designed to be specific to the rDNA 16S sequences of six different bacterial groups: actinomycetes, bacilli, clostridia, cyanobacteria, myxobacteria, and pseudomonads.
For each of these groups, a substantial number of 16S rRNA sequences (see Table 1), chosen among those available in the Ribosomal Database Project II, release 8.0 http://rdp.cme.msu.edu/html/, were imported in GCG Omiga 2.0 (Oxford Molecular Ltd.). In every group, adopting the RDP taxonomic classification, the sequences were assembled in sub-groups and aligned using the Clustal W algorithm, yielding a consensus sequence with a cut off of 75% (meaning that 3 out of 4 sequences determined the consensus at a given position). Then, sub-group consensi were aligned within each group to extract a "group-specific" consensus, adopting the same cut off of 75%. Group-specific probe design was carried out on these "group-specific" consensi. The specificity of each probe pair (common probe and discriminating oligo) was controlled on the RDP II database, using the Probe Match tool. All oligos were designed to have melting temperature (Tm values between 64 and 70°C. Discriminating oligos were purchased with a Cy3 molecule at their 5' terminal position, while common probes with a phosphate in the same position.
Universal microarray preparation
Each of our universal arrays consists of six rows, each corresponding to a group and containing ten replicas. Microarrays were prepared using Code-Link™ activated slides (Motorola Life Sciences), designed to covalently immobilize NH2-modified oligonucleotides. 5' amino-modified Zip Code oligonucleotides, carrying an additional poly(dA)10 tail at their 5'end, were diluted to 25 μM in Printing Buffer (pH 8.5). Spotting was performed using a non-contact piezo-driven dispensing system (Nanoplotter, GeSim, Germany). Printed slides were left overnight in a saturated NaCl chamber with a relative humidity of 75% (this was obtained adding as much solid NaCl to water as needed to form a 1 cm deep slurry in the bottom of a plastic container with an airtight lid). Slide were subsequently placed 20 minutes in a pre-warmed solution (50°C) containing 50 mM ethanolamine, 0,1 M Tris pH 9, 0,1% SDS. They were rinsed twice with water and washed on a shaker for 40 min in 4X SSC/0.1% SDS at 50°C. Finally they were rinsed twice in distilled water and centrifuged at 800 rpm using microplate carriers.
Quality control of printed surfaces was performed by sampling one slide for each deposition batch. The printed slide was hybridized with 1 μM 5' Cy5 labeled poly(dT)10 in a solution containing 5X SSC and 0.1 mg/ml salmon sperm DNA at RT for 3 h, then washed for 15 min in 1X SSC. The fluorescent signal was controlled by laser scanning following procedures described in "Array hybridization and detection".
The DNA region coding for 16S ribosomal RNA was amplified by using universal primers F27 (5'AGAGTTTGATCMTGGCTCAG 3') and R1492 (5'TACGGYTACCTTGTTACGACTT3') which are targeted to universally conserved regions  and permit the amplification of a 1500 bp fragment. For cyanobacteria-specific amplifications, we employed primer 23S30R (5'cctcgcctctgtgtgcctaggt3')  instead of R1492.
PCRs were performed in a PCR Express thermal cycler (Hybaid, England). The reaction mixtures include 500 nM each primer, 200 μM each dNTP, 10 mM Tris-HCl (pH 8.8), 1.5 mM MgCl2, 50 mM KCl, 0.1% (wt/vol) Triton X-100, 1 U of DynaZyme DNA polymerase (Finnzymes OY, Espoo, Finland) and 5–8 ng of genomic DNA, in a final volume of 50 μl. Prior to amplification, DNA was denatured for 5 min at 95°C. Amplification consisted of 30 cycles of 94°C for 45 s, 61°C for 45 s and 72°C for 2 min. After the cycles, an extension step (10 min at 72°C) was performed.
After thermal cycling was complete, 1 μl of proteinase K (1 mg/ml) was added, and the reaction heated at 70°C for 10 min and then quenched at 94°C for 15 min. After this, PCR products were purified by GFX PCR DNA purification kit (Amersham Pharmacia Biotech Inc, Piscataway-NJ), eluted in 50 μl of autoclaved water and quantified by a spectrophotometer.
Ligation detection reaction
Ligation Reaction was carried out in a final volume of 20 μl containing 20 mM Tris-HCl (pH 7.5), 20 mM KCl, 10 mM MgCl2, 0.1% NP40, 0.01 mM ATP, 1 mM DTT, 2 pmol of each discriminating oligo, 2 pmol of each common probe and 1–500 fmol of purified PCR products. The reaction mixture was preheated for 2 min at 94°C and centrifuged in a microcentrifuge for 1 min; then 1 μl of 4 U/μl Pfu DNA ligase (Stratagene, La Jolla, California) was added. The LDR was cycled for 40 rounds of 94°C for 30 sec and 64°C for 4 min in a PCR Express thermal cycler (Hybaid, England).
Array hybridization and detection
In a 0.5-ml microcentrifuge tube, the LDR mix (20 μ) was diluted to obtain 65 μl of hybridization mixture containing 5X SSC and 0.1 mg/ml salmon sperm DNA. The mix, after heating at 94°C for 2 min and chilling on ice, was applied onto the slide under an EasiSeal encase of 2.4 cm2 (Hybaid, England). Hybridization was carried out in the dark at 65°C for one hour and a half, in a temperature-controlled water bath. After removal of the chamber, the microarray was washed for 15 min in pre-warmed (65°C) 1X SSC, 0.1% SDS. Finally, the slide was spinned at 800 rpm for 3 min.
The fluorescent signal was detected at 5 um resolution using a ScanArray® 4000 laser scanning system (Packard GSI Lumonics, Billerica, MA) with green laser for Cy3 dye (λ ex 543 nm/λ em 570 nm). Both the laser and the photomultiplier (PMT) tube power were set at 70–95%.
To quantify the fluorescent intensity of spots the QuantArray® quantitative microarray analysis software was employed (Packard GSI Lumonics).
The authors gratefully acknowledge CNR Target Project Biotecnologie and Biosearch Italia for financially sustaining this work. Elena Busti and Roberta Bordoni were recipient of a fellowship granted by Italian MURST and Biosearch Italia.
- Weissburg WG, Barns SM, Pelletier DA, Lane DJ: 16S ribosomal DNA amplification for phylogenetic study. J Bacteriol. 1991, 173: 697-703.Google Scholar
- Muyzer GG: DGGE/TGGE a method for identifying genes from natural ecosystems. Curr Opin Microbiol. 1999, 2: 317-322. 10.1016/S1369-5274(99)80055-1.View ArticlePubMedGoogle Scholar
- Schwieger F, Tebbe CC: A new approach to utilize PCR-single-strand-conformation polymorphism for 16S rRNA gene-based microbial community analysis. Appl Environ Microbiol. 1998, 64: 4870-4876.PubMed CentralPubMedGoogle Scholar
- Espejo RT, Feijoo CG, Romero J, Vasquez M: PAGE analysis of the heteroduplexes formed between PCR-amplified 16S rRNA genes: estimation of sequence similarity and rDNA complexity. Microbiology. 1998, 144: 1611-1617.View ArticlePubMedGoogle Scholar
- Delwart EL, Shpaer EG, Louwagie J, McCutcham FE, Grez M, Rubsamen-Waigmann H, Mullins JI: Genetic relationships determined by a DNA heteroduplex mobility assay analysis of HIV-1 env genes. Science. 1993, 262: 1257-1261.View ArticlePubMedGoogle Scholar
- Brunk CF, Avaniss-Aghajani E, Brunk CA: A computer analysis of primer and probe hybridization potential with bacterial small subunit rRNA sequences. Appl Environ Microbiol. 1996, 62: 872-879.PubMed CentralPubMedGoogle Scholar
- Clement BG, Kehl LE, DeBord KL, Kitts CL: Terminal Restriction Fragment Patterns (TRFPs), a rapid, PCR-based method for the comparison of complex bacterial communities. J Microbiol Methods. 1998, 31: 135-142. 10.1016/S0167-7012(97)00105-X.View ArticleGoogle Scholar
- Liu W, Marsh LT, Cheng H, Forney LJ: Characterization of microbial diversity by determining terminal restriction fragment length polymorphisms of genes encoding 16S rRNA. Appl Environ Microbiol. 1997, 63: 4516-4522.PubMed CentralPubMedGoogle Scholar
- Moeseneder MM, Arrieta JM, Muyzer G, Winter C, Herndl GJ: Optimization of terminal-restriction fragment length polymorphism analysis for complex marine bacterioplankton communities and comparison with denaturing gradient gel electrophoresis. Appl Environ Microbiol. 1999, 65: 3518-3525.PubMed CentralPubMedGoogle Scholar
- Guschin DY, Mobarry BK, Proudnikov D, Stahl DA, Rittmann BE, Mirzabekov AD: Oligonucleotide microchips as genosensors for determinative and environmental studies in microbiology. Appl Environ Microbiol. 1997, 63: 2397-2402.PubMed CentralPubMedGoogle Scholar
- Rudi K, Skulberg OM, Skulberg R, Jakobsen KS: Application of Sequence-Specific Labeled 16S rRNA Gene Oligonucleotide Probes for Genetic Profiling of Cyanobacterial Abundance and Diversity by Array Hybridization. Appl Environ Microbiol. 2000, 66: 4004-4011. 10.1128/AEM.66.9.4004-4011.2000.PubMed CentralView ArticlePubMedGoogle Scholar
- Small J, Call DR, Brockman FJ, Straub TM, Chandler DP: Direct Detection of 16S rRNA in Soil Extracts by Using Oligonucleotide Microarrays. Appl Environ Microbiol. 2001, 67: 4708-4716. 10.1128/AEM.67.10.4708-4716.2001.PubMed CentralView ArticlePubMedGoogle Scholar
- Gerry NP, Witowsky NE, Day J, Hammer RP, Barany G, Barany F: Universal DNA microarray method for multiplex detection of low abundance point mutations. J Mol Biol. 1999, 292: 251-262. 10.1006/jmbi.1999.3063.View ArticlePubMedGoogle Scholar
- Favis R, Day J, Gerry NP, Phelan C, Narod SA, Barany F: Universal DNA array detection of small insertion and deletions in BRCA1 and BRCA2. Nature Biotechnology. 2000, 18: 561-564. 10.1038/75452.View ArticlePubMedGoogle Scholar
- Maidak BL, Cole JR, Lilburn TG, Parker CT, Saxman PR, Farris RJ, Garrity GM, Olsen GJ, Schmidt TM, Tiedje JM: The RDP-II (Ribosomal Database Project). Nucleic Acids Res. 2001, 29 (1): 173-174. 10.1093/nar/29.1.173.PubMed CentralView ArticlePubMedGoogle Scholar
- Wilmotte A, Van der Auwera G, De Wachter R: Structure of the 16 S ribosomal RNA of the thermophilic cyanobacterium Chlorogloeopsis HTF ('Mastigocladus laminosus HTF') strain PCC7518, and phylogenetic analysis. FEBS Lett. 1993, 317 (1–2): 96-100. 10.1016/0014-5793(93)81499-P.View ArticlePubMedGoogle Scholar
- Lepere C, Wilmotte A, Meyer B: Molecular diversity of Microcystis strains (Cyanophyceae, Chroococcales) based on 16S rDNA sequences. Syst Geogr Pl. 2000, 70: 275-283.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article: verbatim copying and redistribution of this article are permitted in all media for any purpose, provided this notice is preserved along with the article's original URL.