Suitability of partial 16S ribosomal RNA gene sequence analysis for the identification of dangerous bacterial pathogens

AIMS: In a bioterrorism event a rapid tool is needed to identify relevant dangerous bacteria. The aim of the study was to assess the usefulness of partial 16S rRNA gene sequence analysis and the suitability of diverse databases for identifying dangerous bacterial pathogens. METHODS AND RESULTS: For rapid identification purposes a 500-bp fragment of the 16S rRNA gene of 28 isolates comprising Bacillus anthracis, Brucella melitensis, Burkholderia mallei, Burkholderia pseudomallei, Francisella tularensis, Yersinia pestis, and eight genus-related and unrelated control strains was amplified and sequenced. The obtained sequence data were submitted to three public and two commercial sequence databases for species identification. The most frequent reason for incorrect identification was the lack of the respective 16S rRNA gene sequences in the database. CONCLUSIONS: Sequence analysis of a 500-bp 16S rDNA fragment allows the rapid identification of dangerous bacterial species. However, for discrimination of closely related species sequencing of the entire 16S rRNA gene, additional sequencing of the 23S rRNA gene or sequencing of the 16S-23S rRNA intergenic spacer is essential. SIGNIFICANCE AND IMPACT OF THE STUDY: This work provides comprehensive information on the suitability of partial 16S rDNA analysis and diverse databases for rapid and accurate identification of dangerous bacterial pathogens.     

PhytoREF: a reference database of the plastidial 16S rRNA gene of photosynthetic eukaryotes with curated taxonomy

Photosynthetic eukaryotes have a critical role as the main producers in most ecosystems of the biosphere. The ongoing environmental metabarcoding revolution opens the perspective for holistic ecosystems biological studies of these organisms, in particular the unicellular microalgae that often lack distinctive morphological characters and have complex life cycles. To interpret environmental sequences, metabarcoding necessarily relies on taxonomically curated databases containing reference sequences of the targeted gene (or barcode) from identified organisms. To date, no such reference framework exists for photosynthetic eukaryotes. In this study, we built the PhytoREF database that contains 6490 plastidial 16S rDNA reference sequences that originate from a large diversity of eukaryotes representing all known major photosynthetic lineages. We compiled 3333 amplicon sequences available from public databases and 879 sequences extracted from plastidial genomes, and generated 411 novel sequences from cultured marine microalgal strains belonging to different eukaryotic lineages. A total of 1867 environmental Sanger 16S rDNA sequences were also included in the database. Stringent quality filtering and a phylogeny‐based taxonomic classification were applied for each 16S rDNA sequence. The database mainly focuses on marine microalgae, but sequences from land plants (representing half of the PhytoREF sequences) and freshwater taxa were also included to broaden the applicability of PhytoREF to different aquatic and terrestrial habitats. PhytoREF, accessible via a web interface ( http://phytoref.fr ), is a new resource in molecular ecology to foster the discovery, assessment and monitoring of the diversity of photosynthetic eukaryotes using high‐throughput sequencing.     

Novel bacterial taxa in the human microbiome

The human gut harbors thousands of bacterial taxa. A profusion of metagenomic sequence data has been generated from human stool samples in the last few years, raising the question of whether more taxa remain to be identified. We assessed metagenomic data generated by the Human Microbiome Project Consortium to determine if novel taxa remain to be discovered in stool samples from healthy individuals. To do this, we established a rigorous bioinformatics pipeline that uses sequence data from multiple platforms (Illumina GAIIX and Roche 454 FLX Titanium) and approaches (whole-genome shotgun and 16S rDNA amplicons) to validate novel taxa. We applied this approach to stool samples from 11 healthy subjects collected as part of the Human Microbiome Project. We discovered several low-abundance, novel bacterial taxa, which span three major phyla in the bacterial tree of life. We determined that these taxa are present in a larger set of Human Microbiome Project subjects and are found in two sampling sites (Houston and St. Louis). We show that the number of false-positive novel sequences (primarily chimeric sequences) would have been two orders of magnitude higher than the true number of novel taxa without validation using multiple datasets, highlighting the importance of establishing rigorous standards for the identification of novel taxa in metagenomic data. The majority of novel sequences are related to the recently discovered genus Barnesiella, further encouraging efforts to characterize the members of this genus and to study their roles in the microbial communities of the gut. A better understanding of the effects of less-abundant bacteria is important as we seek to understand the complex gut microbiome in healthy individuals and link changes in the microbiome to disease.     

Identification of bacteria from clinical samples using Bartonella alpha-Proteobacteria growth medium

In an effort to overcome historical problems associated with the isolation of Bartonella species from animal and human blood samples, our laboratory developed a novel, chemically modified, insect-based, liquid culture medium (Bartonella alpha-Proteobacteria growth medium, BAPGM). In this study, we describe the isolation of non-Bartonella bacteria from aseptically obtained human blood and tissue samples that were inoculated into BAPGM pre-enrichment culture medium, and were obtained during attempts to define each individuals Bartonella infection status. After incubation for at least 7 days in liquid BAPGM, pre-enriched inoculums were sub-cultured onto a BAPGM/blood agar plate. Bacterial DNA was extracted from pooled plated colonies and amplified using conventional PCR targeting the 16S rRNA gene. Subsequently, amplicons were cloned, sequenced and compared to GenBank database sequences using the BLAST program. Regardless of the patient's Bartonella status, seventeen samples generated only one 16S rDNA sequence, representing the following genera: Arthrobacter, Bacillus, Bartonella, Dermabacter, Methylobacterium, Propionibacterium, Pseudomonas, Staphylococcus and bacteria listed as "non-cultured" in the GenBank database. Alkalibacterium, Arthrobacter, Erwinia, Kineococcus, Methylobacterium, Propionibacterium, Sphingomonas, and Staphylococcus were isolated from nine Bartonella-infected individuals. Co-isolation of Acinetobacter, Sphingomonas, Staphylococcus spp. and bacteria listed as "non-cultured" in the GenBank database was achieved for four samples in which Bartonella spp. were not detected. Despite the phylogenetic limitations of using partial 16S rRNA gene sequencing for species and strain identification, the investigational methodology described in this study may provide a complementary approach for the isolation and identification of bacteria from patient samples.     

Benchmarking of 16S rRNA gene databases using known strain sequences

16S rRNA gene analysis is the most convenient and robust method for microbiome studies. Inaccurate taxonomic assignment of bacterial strains could have deleterious effects as all downstream analyses rely heavily on the accurate assessment of microbial taxonomy. The use of mock communities to check the reliability of the results has been suggested. However, often the mock communities used in most of the studies represent only a small fraction of taxa and are used mostly as validation of sequencing run to estimate sequencing artifacts. Moreover, a large number of databases and tools available for classification and taxonomic assignment of the 16S rRNA gene make it challenging to select the best-suited method for a particular dataset. In the present study, we used authentic and validly published 16S rRNA gene type strain sequences (full length, V3-V4 region) and analyzed them using a widely used QIIME pipeline along with different parameters of OTU clustering and QIIME compatible databases. Data Analysis Measures (DAM) revealed a high discrepancy in ratifying the taxonomy at different taxonomic hierarchies. Beta diversity analysis showed clear segregation of different DAMs. Limited differences were observed in reference data set analysis using partial (V3-V4) and full-length 16S rRNA gene sequences, which signify the reliability of partial 16S rRNA gene sequences in microbiome studies. Our analysis also highlights common discrepancies observed at various taxonomic levels using various methods and databases.     

基于形态特征和ITS序列对7个鹅膏菌属菌株的分类鉴定

以采自浙江省丽水地区的7个鹅膏菌属菌株作为研究材料,在基于形态特征进行初步鉴定的基础上,对7种鹅膏菌的rDNAITS区段进行克隆测序和序列特征比较分析。进一步对ITS序列进行核酸序列数据库GenBank同源性检索比对,将从GenBank检索获得的9个最相似物种的ITS序列连同7种鹅膏菌的ITS序列一起作系统发育分析。结果表明:基于ITS序列对f6、f9和f493个菌株的分子鉴定支持了基于形态特征的鉴定结果,对f5的分子鉴定不支持形态鉴定的结果,f8为鹅膏菌属内某种,f66为鹅膏菌属内某种,并与Amanitafulva,A.atrofusca,A.orientifulva3种鹅膏菌的亲缘关系较近,f7与另外6种鹅膏菌的亲缘关系相差甚远。研究结果提示基于分子水平上的ITS序列分析不能单方面作为大型真菌分类鉴定的可靠依据,可以作为基于传统形态学分类鉴定的辅助参考依据。     

Impact of training sets on classification of high-throughput bacterial 16s rRNA gene surveys

Taxonomic classification of the thousands–millions of 16S rRNA gene sequences generated in microbiome studies is often achieved using a naïve Bayesian classifier (for example, the Ribosomal Database Project II (RDP) classifier), due to favorable trade-offs among automation, speed and accuracy. The resulting classification depends on the reference sequences and taxonomic hierarchy used to train the model; although the influence of primer sets and classification algorithms have been explored in detail, the influence of training set has not been characterized. We compared classification results obtained using three different publicly available databases as training sets, applied to five different bacterial 16S rRNA gene pyrosequencing data sets generated (from human body, mouse gut, python gut, soil and anaerobic digester samples). We observed numerous advantages to using the largest, most diverse training set available, that we constructed from the Greengenes (GG) bacterial/archaeal 16S rRNA gene sequence database and the latest GG taxonomy. Phylogenetic clusters of previously unclassified experimental sequences were identified with notable improvements (for example, 50% reduction in reads unclassified at the phylum level in mouse gut, soil and anaerobic digester samples), especially for phylotypes belonging to specific phyla (Tenericutes, Chloroflexi, Synergistetes and Candidate phyla TM6, TM7). Trimming the reference sequences to the primer region resulted in systematic improvements in classification depth, and greatest gains at higher confidence thresholds. Phylotypes unclassified at the genus level represented a greater proportion of the total community variation than classified operational taxonomic units in mouse gut and anaerobic digester samples, underscoring the need for greater diversity in existing reference databases.     

细菌性阴道病患者阴道分泌物16SrDNA序列分析

目的分析健康妇女及细菌性阴道病(Bacterial vaginosis,BV)患者阴道分泌物16S rDNA序列。方法提取20例健康妇女及40例BV患者阴道分泌物标本中的总DNA,针对细菌16S rDNA保守区设计通用引物进行PCR扩增、克隆、测序,将获得的16S rDNA序列与美国国立生物技术信息中心(NCBI)数据库中的发表序列进行比对,分析克隆群中细菌种类和比例。结果通过阴道分泌物16S rDNA序列分析,发现健康妇女阴道分泌物中以卷曲乳酸杆菌(Lactobacillus crispatus),惰性乳酸杆菌(Lactobacillus iners),加氏乳酸杆菌(Lactobacillus gasseri)为优势菌种,而BV患者阴道菌群种类繁多,以加德纳菌属(Gardnerella)和奇异菌属(Atopobium vaginae)克隆子占较大比例,仅4例患者可见卷曲乳酸杆菌(Lactobacillus crispatus),其他患者均未见有乳酸杆菌克隆子且奇异菌属阴道病患者甲硝唑治疗疗效较差。结论健康妇女和BV患者阴道分泌物菌群种类有较大区别,BV患者在治疗前进行16S rDNA序列分析检测有较大的临床意义。     

High bacterial diversity of a waste gas-degrading community in an industrial biofilter as shown by a 16S rDNA clone library

The bacterial diversity of an industrial biofilter used for waste gas abatement in an animal-rendering plant was investigated. A 16S rDNA clone library was generated and 444 clones were screened using computer-aided amplified ribosomal DNA restriction analysis (ARDRA). Of the screened clones, 60.8% showed unique ARDRA patterns and the remaining 174 clones were clustered into 65 groups. Almost full-length 16S rDNA sequences of 106 clones were determined and 90.5% of the clones were affiliated with the two phyla Proteobacteria and Bacteroidetes. Alpha-, Beta-, and Gammaproteobacteria accounted for 22.1, 17.6 and 18.6% respectively. Minor portions were affiliated with the Actinobacteria (2.0%), Firmicutes and Verrucomicrobia (both 1.0%), and the Deltaproteobacteria and Thermomicrobia (each 0.5%). Only six out of the 106 16S rDNA sequences exhibited similarities of more than 97% to classified bacterial species indicating that a substantial fraction of the clone sequences were derived from unknown taxa. It was also evaluated whether a database containing 281 computer-simulated bacterial rDNA fragment patterns generated from published reference sequences can be used for identification purposes. The data analysis demonstrated that this was possible only for a small number of clones, which were closely related to described bacterial strains. Rarefaction analysis of ARDRA clusters demonstrated that the 444 clones screened are insufficient to describe the entire diversity of the clone library.     

Application of broad-range 16S rRNA PCR amplification and DGGE fingerprinting for detection of tick-infecting bacteria

Ticks play an important role in the transmission of arthropod-borne diseases of viral, protozoal and bacterial origin. The present article describes a molecular-biological based method, which facilitated the broad-range analyses of bacterial communities in ixodid ticks (Ixodes ricinus). DNA was extracted both from single ticks and pooled adult ticks. Eubacterial 16S rRNA gene fragments (16S rDNA) were amplified by polymerase chain reaction (PCR) with broad-range ribosomal primers. Sequences spanning the hypervariable V3 region of the 16S rDNA and representing individual bacterial taxons were separated by denaturing gradient gel electrophoresis (DGGE). For phylogenetic identification, DGGE bands were exised, cloned and sequenced. In addition, we set up a 16S rDNA clone library which was screened by DGGE. Sequences were compared with sequences of known bacteria listed in the GenBank database. A number of bacteria were affiliated with the genera Rickettsia, Bartonella, and Borrelia, which are known to be pathogenic and transmitted by ticks. Two sequences were related to the yet to be cultivated Haemobartonella. To our knowledge, Haemobartonella has never been directly detected in I. ricinus. In addition, members of the genera Staphylococcus, Rhodococcus, Pseudomonas, and Moraxella were detected, which have not been identified in ticks so far. Two bacteria were most closely related to a rickettsial endosymbiont of an Acanthamoeba sp., and to an endosymbiont (Legionellaceae, Coxiella group) of the microarthropod Folsomia candida. The results prove that 16S rDNA genotyping in combination with DGGE analysis is a promising approach for the detection and identification of bacteria infecting ticks, regardless of whether these bacteria are fastidious, obligate intracellular or noncultivable.     

一株番茄青枯菌拮抗菌的鉴定及抗病效果初探

从形态、生理生化、16S rDNA3个方面确定了番茄青枯菌拮抗菌株3-1-16的分类地位。光学显微镜下观察到菌体为杆状细胞,革兰氏染色均匀,并可见菌体染成蓝紫色。透射电镜进一步观察到细胞内有许多颗粒状物质,无伴胞晶体。Biolog鉴定,3-1-16与巨大芽孢杆菌(Bacillus megaterium)具有最高相似率为98%。16S rRNA分析,3-1-16与巨大芽孢杆菌MO31同源性最高为99.4%。聚类分析显示3-1-16与3株巨大芽孢杆菌聚成一支,支持度为100%。生理生化特征及培养特征测定结果表明,菌株3-1-16鉴定为巨大芽孢杆菌(Bacillus megaterium)。盆栽试验表明该菌株对番茄青枯病防病效果达到81.3%。     

A simple binomial test for estimating sequencing errors in public repository 16S rRNA sequences

Sequences in public databases may contain a number of sequencing errors. A double binomial model describing the distribution of indel-excluded similarity coefficients (S) among repeatedly sequenced 16S rRNA was previously developed and it produced a confidence interval of S useful for testing sequence identity among sequences of 400-bp length. We characterized patterns in sequencing errors found in nearly complete 16S rRNA sequences of Vibrionaceae as highly variable in reported sequence length and containing a small number of indels. To accommodate these characteristics, a simple binomial model for distribution of the similarity coefficient (H) that included indels was derived from the double binomial model for S. The model showed good fit to empirical data. By using either a pre-determined or bootstrapping estimated standard probability of base matching, we were able to use the exact binomial test to determine the relative level of sequencing error for a given pair of duplicated sequences. A limitation of the method is the requirement that duplicated sequences for the same template sequence be paired, but this can be overcome by using only conserved regions of 16S rRNA sequences and pairing a given sequence with its highest scoring BLAST search hit from the nr database of GenBank.     

rpoB sequence analysis as a novel basis for bacterial identification

Comparison of the sequences of conserved genes, most commonly those encoding 16S rRNA, is used for bacterial genotypic identification. Among some taxa, such as the Enterobacteriaceae, variation within this gene does not allow confident species identification. We investigated the usefulness of RNA polymerase beta-subunit encoding gene ( rpoB  ) sequences as an alternative tool for universal bacterial genotypic identification. We generated a database of partial rpoB for 14 Enterobacteriaceae species and then assessed the intra- and interspecies divergence between the rpoB and the 16S rRNA genes by pairwise comparisons. We found that levels of divergence between the rpoB sequences of different strains were markedly higher than those between their 16S rRNA genes. This higher discriminatory power was further confirmed by assigning 20 blindly selected clinical isolates to the correct enteric species on the basis of rpoB sequence comparison. Comparison of rpoB sequences from Enterobacteriaceae was also used as the basis for their phylogenetic analysis and demonstrated the genus Klebsiella to be polyphyletic. The trees obtained with rpoB were more compatible with the currently accepted classification of Enterobacteriaceae than those obtained with 16S rRNA. These data indicate that rpoB is a powerful identification tool, which may be useful for universal bacterial identification.     

A Compositional Look at the Human Gastrointestinal Microbiome and Immune Activation Parameters in HIV Infected Subjects

HIV progression is characterized by immune activation and microbial translocation. One factor that may be contributing to HIV progression could be a dysbiotic microbiome. We therefore hypothesized that the GI mucosal microbiome is altered in HIV patients and this alteration correlates with immune activation in HIV. 121 specimens were collected from 21 HIV positive and 22 control human subjects during colonoscopy. The composition of the lower gastrointestinal tract mucosal and luminal bacterial microbiome was characterized using 16S rDNA pyrosequencing and was correlated to clinical parameters as well as immune activation and circulating bacterial products in HIV patients on ART. The composition of the HIV microbiome was significantly different than that of controls; it was less diverse in the right colon and terminal ileum, and was characterized by loss of bacterial taxa that are typically considered commensals. In HIV samples, there was a gain of some pathogenic bacterial taxa. This is the first report characterizing the terminal ileal and colonic mucosal microbiome in HIV patients with next generation sequencing. Limitations include use of HIV-infected subjects on HAART therapy.     

Identifying the plant‐associated microbiome across aquatic and terrestrial environments: the effects of amplification method on taxa discovery

Plants in terrestrial and aquatic environments contain a diverse microbiome. Yet, the chloroplast and mitochondria organelles of the plant eukaryotic cell originate from free‐living cyanobacteria and Rickettsiales. This represents a challenge for sequencing the plant microbiome with universal primers, as ~99% of 16S rRNA sequences may consist of chloroplast and mitochondrial sequences. Peptide nucleic acid clamps offer a potential solution by blocking amplification of host‐associated sequences. We assessed the efficacy of chloroplast and mitochondria‐blocking clamps against a range of microbial taxa from soil, freshwater and marine environments. While we found that the mitochondrial blocking clamps appear to be a robust method for assessing animal‐associated microbiota, Proteobacterial 16S rRNA binds to the chloroplast‐blocking clamp, resulting in a strong sequencing bias against this group. We attribute this bias to a conserved 14‐bp sequence in the Proteobacteria that matches the 17‐bp chloroplast‐blocking clamp sequence. By scanning the Greengenes database, we provide a reference list of nearly 1500 taxa that contain this 14‐bp sequence, including 48 families such as the Rhodobacteraceae, Phyllobacteriaceae, Rhizobiaceae, Kiloniellaceae and Caulobacteraceae. To determine where these taxa are found in nature, we mapped this taxa reference list against the Earth Microbiome Project database. These taxa are abundant in a variety of environments, particularly aquatic and semiaquatic freshwater and marine habitats. To facilitate informed decisions on effective use of organelle‐blocking clamps, we provide a searchable database of microbial taxa in the Greengenes and Silva databases matching various n‐mer oligonucleotides of each PNA sequence.     

Discovery and characterization of a new bacterial candidate division by an anaerobic sludge digester metagenomic approach

We have constructed a large fosmid library from a mesophilic anaerobic digester and explored its 16S rDNA diversity using a high-density filter DNA–DNA hybridization procedure. We identified a group of 16S rDNA sequences forming a new bacterial lineage named WWE3 (Waste Water of Evry 3). Only one sequence from the public databases shares a sequence identity above 80% with the WWE3 group which hence cannot be affiliated to any known or candidate prokaryotic division. Despite representing a non-negligible fraction (5% of the 16S rDNA sequences) of the bacterial population of this digester, the WWE3 bacteria could not have been retrieved using the conventional 16S rDNA amplification procedure due to their unusual 16S rDNA gene sequence. WWE3 bacteria were detected by polymerase chain reaction (PCR) in various environments (anaerobic digesters, swine lagoon slurries and freshwater biofilms) using newly designed specific PCR primer sets. Fluorescence in situ hybridization (FISH) analysis of sludge samples showed that WWE3 microorganisms are oval-shaped and located deep inside sludge flocs. Detailed phylogenetic analysis showed that WWE3 bacteria form a distinct monophyletic group deeply branching apart from all known bacterial divisions. A new bacterial candidate division status is proposed for this group.     

Genomic BLAST: custom-defined virtual databases for complete and unfinished genomes

BLAST (Basic Local Alignment Search Tool) searches against DNA and protein sequence databases have become an indispensable tool for biomedical research. The proliferation of the genome sequencing projects is steadily increasing the fraction of genome-derived sequences in the public databases and their importance as a public resource. We report here the availability of Genomic BLAST, a novel graphical tool for simplifying BLAST searches against complete and unfinished genome sequences. This tool allows the user to compare the query sequence against a virtual database of DNA and/or protein sequences from a selected group of organisms with finished or unfinished genomes. The organisms for such a database can be selected using either a graphic taxonomy-based tree or an alphabetical list of organism-specific sequences. The first option is designed to help explore the evolutionary relationships among organisms within a certain taxonomy group when performing BLAST searches. The use of an alphabetical list allows the user to perform a more elaborate set of selections, assembling any given number of organism-specific databases from unfinished or complete genomes. This tool, available at the NCBI web site http://www.ncbi.nlm.nih.gov/cgi-bin/Entrez/genom_table_cgi, currently provides access to over 170 bacterial and archaeal genomes and over 40 eukaryotic genomes.     

Phylogenetic Diversity of Archaea and Bacteria in a Deep Subsurface Paleosol

Abstract A low-biomass paleosol 188 m below the ground surface at the Department of Energy's Hanford Site in south-central Washington State was recovered and maintained at the in situ temperature (17°C) as an intact core or homogenized sediment for 0, 1, 3, 10, and 21 weeks post-sampling. Bacterial and archaeal 16S rRNA genes were amplified by PCR and cloned. Of 746 bacterial and 190 archaeal clones that were categorized by restriction fragment length polymorphism (RFLP), 242 bacterial and 16 archaeal clones were partially sequenced and compared against the small subunit ribosomal RNA database (RDP) and GenBank. Six bacterial and 16 archaeal clones sequences, with little similarity to those in public databases, were sequenced in their entirety, and subjected to more detained phylogenetic analysis. The most frequently occurring clones types were related to Pseudomonas, Bacillus, Micrococcus, Clavibacter, Nocardioides, Burkholderia, Comamonas, and Erythromicrobium. Clone sequences whose RDP similarity value was ≥0.6 consistently grouped with their nearest RDP neighbor during phylogenetic analysis. Six truly novel eubacterial sequences were identified; they consistently cluster with or near the Chloroflexaceae and sequences recovered from the Sargasso Sea. Sixteen unique archaeal RFLP groups were identified from 190 randomly-sampled clones. The novel archaeal rDNA clones formed a coherent clade along the major Crenarchaea branch containing all previously described mesophilic crenarchae clones, but remained firmly associated with 16S rDNA clones previously obtained from a thermal Fe/S spring in Yellowstone National Park. The wealth of group-specific genetic information identified during this study will now allow us to address specific hypotheses related to in situ stimulation of these deep subsurface microorganisms and changes in microbial community composition resulting from subsurface contamination or remediation processes at the Hanford Site. Revised: 21 October 1997; Accepted: 20 November 1997