Potential of a 16S rRNA-based taxonomic microarray for analyzing the rhizosphere effects of maize on Agrobacterium spp. and bacterial communities

Bacterial diversity is central to ecosystem sustainability and soil biological function, for which the role of roots is important. The high-throughput analysis potential of taxonomic microarray should match the breadth of bacterial diversity. Here, the power of this technology was evidenced through methodological verifications and analysis of maize rhizosphere effect based on a 16S rRNA-based microarray developed from the prototype of H. Sanguin et al. (Environ. Microbiol. 8:289-307, 2006). The current probe set was composed of 170 probes (41 new probes in this work) that targeted essentially the Proteobacteria. Cloning and sequencing of 16S rRNA amplicons were carried out on maize rhizosphere and bulk soil DNA. All tested clones that had a perfect match with corresponding probes were positive in the hybridization experiment. The hierarchically nested probes were reliable, but the level of taxonomic identification was variable, depending on the probe set specificity. The comparison of experimental and theoretical hybridizations revealed 0.91% false positives and 0.81% false negatives. The microarray detection threshold was estimated at 0.03% of a given DNA type based on DNA spiking experiments. A comparison of the maize rhizosphere and bulk soil hybridization results showed a significant rhizosphere effect, with a higher predominance of Agrobacterium spp. in the rhizosphere, as well as a lower prevalence of Acidobacteria, Bacteroidetes, Verrucomicrobia, and Planctomycetes, a new taxon of interest in soil. In addition, well-known taxonomic groups such as Sphingomonas spp., Rhizobiaceae, and Actinobacteria were identified in both microbial habitats with strong hybridization signals. The taxonomic microarray developed in the present study was able to discriminate and characterize bacterial community composition in related biological samples, offering extensive possibilities for systematic exploration of bacterial diversity in ecosystems.     

Development of a 16S rRNA gene-based prototype microarray for the detection of selected actinomycetes genera

Actinomycetes are known for their secondary metabolites, which have been successfully used as drugs in human and veterinary medicines. However, information on the distribution of this group of Gram-positive bacteria in diverse ecosystems and a comprehension of their activities in ecosystem processes are still scarce. We have developed a 16S rRNA-based taxonomic microarray that targets key actinomycetes at the genus level. In total, 113 actinomycete 16S rRNA probes, corresponding to 55 of the 202 described genera, were designed. The microarray accuracy was evaluated by comparing signal intensities with probe/target-weighted mismatch values and the Gibbs energy of the probe/target duplex formation by hybridizing 17 non-actinomycete and 29 actinomycete strains/clones with the probe set. The validation proved that the probe set was specific, with only 1.3% of false results. The incomplete coverage of actinomycetes by a genus-specific probe was caused by the limited number of 16S rRNA gene sequences in databases or insufficient 16S rRNA gene polymorphism. The microarray enabled discrimination between actinomycete communities from three forest soil samples collected at one site. Cloning and sequencing of 16S rRNA genes from one of the soil samples confirmed the microarray results. We propose that this newly constructed microarray will be a valuable tool for genus-level comparisons of actinomycete communities in various ecological conditions.     

Application and validation of DNA microarrays for the 16S rRNA-based analysis of marine bacterioplankton

An oligonucleotide probe-based DNA microarray was evaluated for its ability to detect 16S rRNA targets in marine bacterioplankton samples without prior amplification by polymerase chain reaction (PCR). The results obtained were compared with those of quantitative fluorescence in situ hybridization (FISH). For extraction and direct labelling of total RNA, a fast and efficient protocol based on commercially available kits was established. A set of redundant and hierarchically structured probes was applied, and specificity of hybridization was assessed by additional control oligonucleotides comprising single central mismatches. The protocol was initially tested by microarray analysis of bacterial pure cultures. Complete discrimination of all control oligonucleotides was achieved, indicating a high degree of hybridization specificity. In a co-culture, abundant members were detected by microarray analysis, but signal ratios of positive probes did not correlate well with quantitative data from FISH experiments. A marine picoplankton sample from the German Bight was analysed. Bacterial populations with relative abundances of at least 5% were detected by hybridizing 0.1 microg of total RNA extracted from a sample of 375 ml equivalent to 4.1 x 10(8) cells. Our results demonstrate that major populations of marine bacterioplankton can be identified by microarray analysis in a fast and reliable way, even in relatively low volumes of sea water.     

Construction and validation of a prototype microarray for efficient and high-throughput genotyping of angiosperms

Until recently, the identification of plants relied on conventional techniques, such as morphological, anatomical and chemical profiling, that are often inefficient or unfeasible in certain situations. Extensive literature exists describing the use of polymerase chain reaction (PCR) DNA-based identification techniques, which offer a reliable platform, but their broad application is often limited by a low throughput. However, hybridization-based microarray technology represents a rapid and high-throughput tool for genotype identification at a molecular level. Using an innovative technique, a 'Subtracted Diversity Array' (SDA) of 376 features was constructed from a pooled genomic DNA library of 49 angiosperm species, from which pooled non-angiosperm genomic DNA was subtracted. Although not the first use of a subtraction technique for genotyping, the SDA method was superior in accuracy, sensitivity and efficiency, and showed high-throughput capacity and broad application. The SDA technique was validated for potential genotyping use, and the results indicated a successful subtraction of non-angiosperm DNA. Statistical analysis of the polymorphic features from the pilot study enabled the establishment of accurate phylogenetic relationships, confirming the potential use of the SDA technique for genotyping. Further, the technique substantially enriched the presence of polymorphic sequences; 68% were polymorphic when using the array to differentiate six angiosperm clades (Asterids, Rosids, Caryophyllids, Ranunculids, Monocots and Eumagnoliids). The 'proof of concept' experiments demonstrate the potential of establishing a highly informative, reliable and high-throughput microarray-based technique for novel application to sequence independent genotyping of major angiosperm clades.     

Designation of Streptomycete 16S and 23S rRNA-based target regions for oligonucleotide probes

The 16S and 23S rRNA of various Streptomyces species were partially sequenced and screened for the presence of stretches that could define all members of the genus, groups of species, or individual species. Nucleotide 929 (Streptomyces ambofaciens nomenclature [J.L. Pernodet, M.T. Alegre, F. Boccard, and M. Guerineau, Gene 79:33-46, 1989]) is a nucleotide highly unique to Streptomyces species which, in combination with flanking regions, allowed the designation of a genus-specific probe. Regions 158 through 203 of the 16S rRNA and 1518 through 1645 of the 23S rRNA (helix 54 [Pernodet et al., Gene 79:33-46, 1989]) have a high potential to define species, whereas the degree of variation in regions 982 through 998 and 1102 through 1122 of the 16S rRNA is less pronounced but characteristic for at least certain species. Alone or in combination with each other, these regions may serve as target sites for synthetic oligonucleotide probes and primers to be used in the determination of pure cultures and in the characterization of community structures. The specificity of several probes is demonstrated by dot blot hybridization.     

Designation of Streptomycete 16S and 23S rRNA-based target regions for oligonucleotide probes.

The 16S and 23S rRNA of various Streptomyces species were partially sequenced and screened for the presence of stretches that could define all members of the genus, groups of species, or individual species. Nucleotide 929 (Streptomyces ambofaciens nomenclature [J.L. Pernodet, M.T. Alegre, F. Boccard, and M. Guerineau, Gene 79:33-46, 1989]) is a nucleotide highly unique to Streptomyces species which, in combination with flanking regions, allowed the designation of a genus-specific probe. Regions 158 through 203 of the 16S rRNA and 1518 through 1645 of the 23S rRNA (helix 54 [Pernodet et al., Gene 79:33-46, 1989]) have a high potential to define species, whereas the degree of variation in regions 982 through 998 and 1102 through 1122 of the 16S rRNA is less pronounced but characteristic for at least certain species. Alone or in combination with each other, these regions may serve as target sites for synthetic oligonucleotide probes and primers to be used in the determination of pure cultures and in the characterization of community structures. The specificity of several probes is demonstrated by dot blot hybridization.     

Development and validation of a diagnostic microbial microarray for methanotrophs

The potential of DNA microarray technology in high-throughput detection of bacteria and quantitative assessment of their community structures is widely acknowledged but has not been fully realised yet. A generally applicable set of techniques, based on readily available technologies and materials, was developed for the design, production and application of diagnostic microbial microarrays. A microarray targeting the particulate methane monooxygenase (pmoA) gene was developed for the detection and quantification of methanotrophs and functionally related bacteria. A microarray consisting of a set of 59 probes that covers the whole known diversity of these bacteria was validated with a representative set of extant strains and environmental clones. The potential of the pmoA microarray was tested with environmental samples. The results were in good agreement with those of clone library sequence analyses. The approach can currently detect less dominant bacteria down to 5% of the total community targeted. Initial tests assessing the quantification potential of this system with artificial PCR mixtures showed very good correlation with the expected results with standard deviations in the range of 0.4-17.2%. Quantification of environmental samples with this method requires the design of a reference mixture consisting of very close relatives of the strains within the sample and is currently limited by biases inherent in environmental DNA extraction and universal PCR amplification.     

Novel primers for 16S rRNA-based archaeal community analyses in environmental samples

Next generation sequencing technologies for in depth analyses of complex microbial communities rely on rational primer design based on up-to-date reference databases. Most of the 16S rRNA-gene based analyses of environmental Archaea community composition use PCR primers developed from small data sets several years ago, making an update long overdue. Here we present a new set of archaeal primers targeting the 16S rRNA gene designed from 8500 aligned archaeal sequences in the SILVA database. The primers 340F-1000R showed a high archaeal specificity (< 1% bacteria amplification) covering 93 and 97% of available sequences for Crenarchaeota and Euryarchaeota respectively. In silico tests of the primers revealed at least 38% higher coverage for Archaea compared to other commonly used primers. Empirical tests with clone libraries confirmed the high specificity of the primer pair to Archaea in three biomes: surface waters in the Arctic Ocean, the pelagic zone of a temperate lake and a methanogenic bioreactor. The clone libraries featured both Euryarchaeota and Crenarchaeota in variable proportions and revealed dramatic differences in the archaeal community composition and minimal phylogenetic overlap between samples.     

Development of a 23S rRNA-based PCR assay for the identification of Pasteurella multocida

AIMS: The aim of this work was to develop a rapid diagnostic test for Pasteurella multocida. METHODS AND RESULTS: A polymerase chain reaction (PCR) assay using primers derived from the 23S rRNA gene sequence of Past. multocida was developed. The PCR assay correctly identified all 144 isolates of Past. multocida tested, including type strains of the three subspecies as well as the reference strains for the Heddleston and Carter typing schemes. Of 20 closely related bacteria from the family Pasteurellaceae tested, only the type strains of Past. canis biovar 2 and Past. avium biovar 2 were positive. These two bacteria, formerly known as Bisgaard Taxon 13, are the closest phylogenetic relatives of Past. multocida based on 16S ribosomal rRNA. All phylogenetically unrelated avian and porcine organisms tested were negative. CONCLUSION: This PCR enables rapid identification of Past. multocida colonies from avian or porcine origin. SIGNIFICANCE AND IMPACT OF THE STUDY: Veterinary diagnostic laboratories can use this PCR to rapidly and accurately diagnose fowl cholera and porcine pasteurellosis.     

Polyphyly in 16S rRNA-based LVTree Versus Monophyly in Whole-genome-based CVTree

We report an important but long-overlooked manifestation of low-resolution power of 16S rRNA sequence analysis at the species level, namely, in 16S rRNA-based phylogenetic trees polyphyletic placements of closely-related species are abundant compared to those in genome-based phylogeny. This phenomenon makes the demarcation of genera within many families ambiguous in the 16S rRNA-based taxonomy. In this study, we reconstructed phylogenetic relationship for more than ten thousand prokaryote genomes using the CVTree method, which is based on whole-genome information. And many such genera, which are polyphyletic in 16S rRNA-based trees, are well resolved as monophyletic clusters by CVTree. We believe that with genome sequencing of prokaryotes becoming a commonplace, genome-based phylogeny is doomed to play a definitive role in the construction of a natural and objective taxonomy.     

Extensive set of 16S rRNA-based probes for detection of bacteria in human feces

For the detection of six groups of anaerobic bacteria in human feces, we designed seven new 16S rRNA-based oligonucleotide probes. This set of probes extends the current set of probes and gives more data on the composition of the human gut flora. Probes were designed for Phascolarctobacterium and relatives (Phasco741), Veillonella (Veil223), Eubacterium hallii and relatives (Ehal1469), Lachnospira and relatives (Lach571), and Eubacterium cylindroides and relatives (Ecyl387), and two probes were designed for Ruminococcus and relatives (Rbro730 and Rfla729). The hybridization conditions for the new probes were optimized for fluorescent in situ hybridization, and the probes were validated against a set of reference organisms. The probes were applied to fecal samples of 11 volunteers to enumerate their target bacterial groups. The Phasco741 and Veil223 probes both detected average numbers below 1% of the total number of bacteria as determined with the bacterial kingdom-specific Bact338 probe. The Ecyl387 probe detected about 1.4%, the Lach571 and Ehal1469 probes detected 3.8 and 3.6%, respectively, and a combination of the Rbro730 and Rfla729 probes detected 10.3%. A set of 15 probes consisting of probes previously described and those presented here were evaluated in hybridization with the fecal samples of the same volunteers. Together, the group-specific probes detected 90% of the total bacterial cells.     

Retrieval of nearly complete 16S rRNA gene sequences from environmental DNA following 16S rRNA-based community fingerprinting

16S rRNA-based fingerprinting techniques allow rapid analyses of overall bacterial community structure but suffer from a lack of phylogenetic information hitherto retrievable from the short 16S rRNA gene sequences obtained from excised bands. An approach is presented that allows nearly complete 16S rRNA gene sequences to be retrieved for abundant components of the bacterial community as obtained by the community fingerprint, i.e. those reflected by major fingerprint bands. This was achieved by designing a pair of highly specific primers derived from the sequence of an excised band. Combined with universal 16S rRNA primers, these specific primers were applied directly to environmental DNA serving as template. This procedure allowed the generation of a nearly complete 16S rRNA gene sequence of the target taxon by specific polymerase chain reaction (PCR) followed by cycle sequencing down to a relative abundance of at least 1.5% of the environmental DNA. The procedure was exemplified for an epsilonproteobacterium related to Thiomicrospira denitrificans occurring in the central Baltic Sea. This approach is based only on PCR without any cloning step involved. It allows focussing on specific target taxa and is thus rather efficient. This approach should be applicable in general to 16S rRNA or 16S rRNA gene-based fingerprinting techniques and their respective environmental DNA.     

Development and evaluation of 16S rDNA microarray for detecting bacterial pathogens in cerebrospinal fluid

The rapid identification of bacteria in cerebrospinal fluid (CSF) is very important for patient management and antimicrobial therapies. We developed a 16S DNA microarray-based method that targets 16S rDNA and can directly detect bacteria from CSF without cultivation. Universal primers and specific probes were designed from the 16S rDNA sequence data retrieved directly from the GenBank database. The specificity of the assay is obtained through a combination of microarray hybridization and enzymatic labeling of the constructed specific probes. Cultivation-dependent assays were used as reference methods in the development and evaluation of the method. With the exception of Mycobacterium tuberculosis and Proteus mirabilis, forty-five positive blood culture media were successfully differentiated. When this procedure was applied directly to 100 CSF specimens, 29 specimens from 16 patients were positive by bacterial culture and 3 culture-positive CSF specimens produced no hybridized signals. The remaining 26 specimens were correctly identified, including one with mixed infection. The accuracy, sensitivity, and specificity of the assay can be increased further by designing more oligonucleotides for the microarray. This method is versatile and makes it possible to detect more bacteria in a single assay and discriminate different bacterial genera.     

C16S - a Hidden Markov Model based algorithm for taxonomic classification of 16S rRNA gene sequences

Recent advances in high throughput sequencing technologies and concurrent refinements in 16S rDNA isolation techniques have facilitated the rapid extraction and sequencing of 16S rDNA content of microbial communities. The taxonomic affiliation of these 16S rDNA fragments is subsequently obtained using either BLAST-based or word frequency based approaches. However, the classification accuracy of such methods is observed to be limited in typical metagenomic scenarios, wherein a majority of organisms are hitherto unknown. In this study, we present a 16S rDNA classification algorithm, called C16S, that uses genus-specific Hidden Markov Models for taxonomic classification of 16S rDNA sequences. Results obtained using C16S have been compared with the widely used RDP classifier. The performance of C16S algorithm was observed to be consistently higher than the RDP classifier. In some scenarios, this increase in accuracy is as high as 34%. A web-server for the C16S algorithm is available at http://metagenomics.atc.tcs.com/C16S/.     

Development and validation of a modified broad-range 16S rDNA PCR for diagnostic purposes in clinical microbiology

Broad-range PCR followed by sequencing identifies bacterial pathogens, even in challenging settings such as patients receiving antibiotics or infected with fastidious or non-cultivable organisms. The major problem with broad-range PCR is the risk of sample contamination. Risk is present at every step of the procedure, starting from sample collection. Contaminating bacterial DNA may be present not only in laboratory reagents but also at the surface of plastic consumables and containers used for specimen drawing and transport to the diagnostic laboratory. Contaminating DNA is amplified efficiently, leading to false-positive results. Thus, high specificity depends on eliminating such spurious targets, an awkward problem given the abundance of such targets and a highly sensitive method that detects very small numbers of molecules. Several investigators have reported strategies for eliminating the amplification of contaminating DNA sequences. So far, none of these methods has been entirely effective and reproducible. Here we describe a method that uses Exonuclease III (ExoIII) to disable contaminating sequences from acting as templates, while maintaining the high sensitivity of PCR for pathogen DNA. We use this assay in 144 clinical specimens from normally sterile sites, identifying pathogens from 24 (17%). Conventional methods identified pathogens in only four of these specimens, all of which were positive for the same pathogen by PCR. Compared with conventional methods, broad-range PCR with ExoIII pre-treatment of reagents substantially improves the diagnostic yield of bacterial pathogen identification from normally sterile sites.     

16S rRNA-based analysis of microbiota from the cecum of broiler chickens.

The microbiota of the intestinal tract of chickens plays an important role in inhibiting the establishment of intestinal pathogens. Earlier culturing and microscopic examinations indicated that only a fraction of the bacteria in the cecum of chickens could be grown in the laboratory. Therefore, a survey of cecal bacteria was done by retrieval of 16S rRNA gene sequences from DNA isolated from the cecal content and the cecal mucosa. The ribosomal gene sequences were amplified with universal primers and cloned or subjected to temporal temperature gradient gel electrophoresis (TTGE). Partial 16S rRNA gene sequences were determined from the clones and from the major bands in TTGE gels. A total of 1,656 partial 16S rRNA gene sequences were obtained and compared to sequences in the GenBank. The comparison indicated that 243 different sequences were present in the samples. Overall, sequences representing 50 phylogenetic groups or subgroups of bacteria were found, but approximately 89% of the sequences represented just four phylogenetic groups (Clostridium leptum, Sporomusa sp., Clostridium coccoides, and enterics). Sequences of members of the Bacteroides group, the Bifidobacterium infantis subgroup, and of Pseudomonas sp. each accounted for less than 2% of the total. Sequences related to those from the Escherichia sp. subgroup and from Lactobacillus, Pseudomonas, and Bifidobacterium spp. were generally between 98 and 100% identical to sequences already deposited in the GenBank. Sequences most closely related to those of the other bacteria were generally 97% or less identical to those in the databases and therefore might be from currently unknown species. TTGE and random cloning indicated that certain phylogenetic subgroups were common to all birds analyzed, but sequence data from random cloning also provided evidence for qualitative and quantitative differences among the cecal microbiota of individual birds reared under very similar conditions.     

Reducing the effects of PCR amplification and sequencing artifacts on 16S rRNA-based studies

The advent of next generation sequencing has coincided with a growth in interest in using these approaches to better understand the role of the structure and function of the microbial communities in human, animal, and environmental health. Yet, use of next generation sequencing to perform 16S rRNA gene sequence surveys has resulted in considerable controversy surrounding the effects of sequencing errors on downstream analyses. We analyzed 2.7×10(6) reads distributed among 90 identical mock community samples, which were collections of genomic DNA from 21 different species with known 16S rRNA gene sequences; we observed an average error rate of 0.0060. To improve this error rate, we evaluated numerous methods of identifying bad sequence reads, identifying regions within reads of poor quality, and correcting base calls and were able to reduce the overall error rate to 0.0002. Implementation of the PyroNoise algorithm provided the best combination of error rate, sequence length, and number of sequences. Perhaps more problematic than sequencing errors was the presence of chimeras generated during PCR. Because we knew the true sequences within the mock community and the chimeras they could form, we identified 8% of the raw sequence reads as chimeric. After quality filtering the raw sequences and using the Uchime chimera detection program, the overall chimera rate decreased to 1%. The chimeras that could not be detected were largely responsible for the identification of spurious operational taxonomic units (OTUs) and genus-level phylotypes. The number of spurious OTUs and phylotypes increased with sequencing effort indicating that comparison of communities should be made using an equal number of sequences. Finally, we applied our improved quality-filtering pipeline to several benchmarking studies and observed that even with our stringent data curation pipeline, biases in the data generation pipeline and batch effects were observed that could potentially confound the interpretation of microbial community data.     

Endophytic Pseudomonasspp. populations of pathogen-infected potato plants analysed by 16S rDNA- and 16S rRNA-based denaturating gradient gel electrophoresis

The diversity of abundant and metabolically active pseudomonads in potato plants was analysed using a culture-independent approach. The effect of two plant varieties, Agria and Bionta, as well as the presence of a plant pathogen Erwinia carotovora ssp. atroseptica on this bacterial group was tested. A combination of Pseudomonas-specific PCR, DGGE analysis, cloning and sequencing of partial 16S rDNA genes was performed using DNA and RNA extracted from potato stem tissue. Sequence analysis revealed a high species diversity, with the most prominent ones being Pseudomonas stutzeri and Pseudomonas gingeri. Some species showed high rRNA contents indicating high metabolic activity. Both, highly abundant and metabolically active Pseudomonaspopulations were more affected by the plant genotype than by the presence of E. carotovora.