Rapid 16S rRNA Next-Generation Sequencing of Polymicrobial Clinical Samples for Diagnosis of Complex Bacterial Infections

Classifying individual bacterial species comprising complex, polymicrobial patient specimens remains a challenge for culture-based and molecular microbiology techniques in common clinical use. We therefore adapted practices from metagenomics research to rapidly catalog the bacterial composition of clinical specimens directly from patients, without need for prior culture. We have combined a semiconductor deep sequencing protocol that produces reads spanning 16S ribosomal RNA gene variable regions 1 and 2 (∼360 bp) with a de-noising pipeline that significantly improves the fraction of error-free sequences. The resulting sequences can be used to perform accurate genus- or species-level taxonomic assignment. We explore the microbial composition of challenging, heterogeneous clinical specimens by deep sequencing, culture-based strain typing, and Sanger sequencing of bulk PCR product. We report that deep sequencing can catalog bacterial species in mixed specimens from which usable data cannot be obtained by conventional clinical methods. Deep sequencing a collection of sputum samples from cystic fibrosis (CF) patients reveals well-described CF pathogens in specimens where they were not detected by standard clinical culture methods, especially for low-prevalence or fastidious bacteria. We also found that sputa submitted for CF diagnostic workup can be divided into a limited number of groups based on the phylogenetic composition of the airway microbiota, suggesting that metagenomic profiling may prove useful as a clinical diagnostic strategy in the future. The described method is sufficiently rapid (theoretically compatible with same-day turnaround times) and inexpensive for routine clinical use.

[1]  W R Pearson,et al.  Flexible sequence similarity searching with the FASTA3 program package. , 2000, Methods in molecular biology.

[2]  C. Quince,et al.  Accurate determination of microbial diversity from 454 pyrosequencing data , 2009, Nature Methods.

[3]  H. Swerdlow,et al.  A tale of three next generation sequencing platforms: comparison of Ion Torrent, Pacific Biosciences and Illumina MiSeq sequencers , 2012, BMC Genomics.

[4]  S. Sagel,et al.  The airway microbiome in cystic fibrosis and implications for treatment , 2011, Current opinion in pediatrics.

[5]  Frederick A. Matsen IV,et al.  Edge Principal Components and Squash Clustering: Using the Special Structure of Phylogenetic Placement Data for Sample Comparison , 2011, PloS one.

[6]  R. Knight,et al.  Rapid denoising of pyrosequencing amplicon data: exploiting the rank-abundance distribution , 2010, Nature Methods.

[7]  Robert Schlaberg,et al.  A Systematic Approach for Discovering Novel, Clinically Relevant Bacteria , 2012, Emerging infectious diseases.

[8]  P. Schreckenberger,et al.  The Role of 16S rRNA Gene Sequencing in Identification of Microorganisms Misidentified by Conventional Methods , 2005, Journal of Clinical Microbiology.

[9]  Jens Stoye,et al.  Bacterial Community Shift in Treated Periodontitis Patients Revealed by Ion Torrent 16S rRNA Gene Amplicon Sequencing , 2012, PloS one.

[10]  M. Guibert,et al.  Ralstonia pickettii Traced in Blood Culture Bottles , 2002, Journal of Clinical Microbiology.

[11]  Kun Tang,et al.  Comparative analysis of human saliva microbiome diversity by barcoded pyrosequencing and cloning approaches. , 2009, Analytical biochemistry.

[12]  Frederick Albert Matsen IV,et al.  A Format for Phylogenetic Placements , 2012, PloS one.

[13]  D. Cowan,et al.  Review and re-analysis of domain-specific 16S primers. , 2003, Journal of microbiological methods.

[14]  P. Kämpfer,et al.  Herminiimonas contaminans sp. nov., isolated as a contaminant of biopharmaceuticals. , 2013, International journal of systematic and evolutionary microbiology.

[15]  J. Clarridge,et al.  Impact of 16 S rRNA Gene Sequence Analysis for Identification of Bacteria on Clinical Microbiology and Infectious Diseases , 2004 .

[16]  Marcus J. Claesson,et al.  Comparison of two next-generation sequencing technologies for resolving highly complex microbiota composition using tandem variable 16S rRNA gene regions , 2010, Nucleic acids research.

[17]  A. van Belkum,et al.  Comparison of conventional and molecular methods for the detection of bacterial pathogens in sputum samples from cystic fibrosis patients. , 2000, FEMS immunology and medical microbiology.

[18]  Hanlee P. Ji,et al.  Next-generation DNA sequencing , 2008, Nature Biotechnology.

[19]  E. Virginia Armbrust,et al.  pplacer: linear time maximum-likelihood and Bayesian phylogenetic placement of sequences onto a fixed reference tree , 2010, BMC Bioinformatics.

[20]  M. Surette,et al.  Culture Enriched Molecular Profiling of the Cystic Fibrosis Airway Microbiome , 2011, PloS one.

[21]  P. Vandamme,et al.  Taxonomy of the genus Cupriavidus: a tale of lost and found. , 2004, International journal of systematic and evolutionary microbiology.

[22]  Zaid Abdo,et al.  Temporal Dynamics of the Human Vaginal Microbiota , 2012, Science Translational Medicine.

[23]  James R. Cole,et al.  The Ribosomal Database Project: improved alignments and new tools for rRNA analysis , 2008, Nucleic Acids Res..

[24]  P. Salamon,et al.  Cystic fibrosis therapy: a community ecology perspective. , 2013, American journal of respiratory cell and molecular biology.

[25]  Daniel J. Wilson,et al.  Transforming clinical microbiology with bacterial genome sequencing , 2012, Nature Reviews Genetics.

[26]  J. Kellogg,et al.  Corynebacterium pseudodiphtheriticum: a respiratory tract pathogen in adults. , 1995, Clinical infectious diseases : an official publication of the Infectious Diseases Society of America.

[27]  S. Batzoglou,et al.  Bacterial flora-typing with targeted, chip-based Pyrosequencing , 2007, BMC Microbiology.

[28]  K. Kinzler,et al.  Detection and quantification of rare mutations with massively parallel sequencing , 2011, Proceedings of the National Academy of Sciences.

[29]  D. Persing,et al.  Elimination of contaminating DNA within polymerase chain reaction reagents: implications for a general approach to detection of uncultured pathogens , 1993, Journal of clinical microbiology.

[30]  D. Vullo,et al.  Bacterial swimming, swarming and chemotactic response to heavy metal presence: which could be the influence on wastewater biotreatment efficiency? , 2012, World Journal of Microbiology and Biotechnology.

[31]  A. J. Jones,et al.  At Least 1 in 20 16S rRNA Sequence Records Currently Held in Public Repositories Is Estimated To Contain Substantial Anomalies , 2005, Applied and Environmental Microbiology.

[32]  J. Clarridge,et al.  Impact of 16S rRNA Gene Sequence Analysis for Identification of Bacteria on Clinical Microbiology and Infectious Diseases , 2004, Clinical Microbiology Reviews.

[33]  N. Pace,et al.  Molecular identification of bacteria in bronchoalveolar lavage fluid from children with cystic fibrosis , 2007, Proceedings of the National Academy of Sciences.

[34]  Bernard P. Puc,et al.  An integrated semiconductor device enabling non-optical genome sequencing , 2011, Nature.

[35]  Jeet Sukumaran,et al.  DendroPy: a Python library for phylogenetic computing , 2010, Bioinform..

[36]  Sean R. Eddy,et al.  Infernal 1.0: inference of RNA alignments , 2009, Bioinform..

[37]  D. Raoult,et al.  Prospects for the future using genomics and proteomics in clinical microbiology. , 2011, Annual review of microbiology.

[38]  S. Giovannoni,et al.  Bias caused by template annealing in the amplification of mixtures of 16S rRNA genes by PCR , 1996, Applied and environmental microbiology.

[39]  G. Olsen,et al.  Critical Evaluation of Two Primers Commonly Used for Amplification of Bacterial 16S rRNA Genes , 2008, Applied and Environmental Microbiology.

[40]  M. Surette,et al.  Modulation of Pseudomonas aeruginosa gene expression by host microflora through interspecies communication , 2003, Molecular microbiology.

[41]  Didier Raoult,et al.  16S Ribosomal DNA Sequence Analysis of a Large Collection of Environmental and Clinical Unidentifiable Bacterial Isolates , 2000, Journal of Clinical Microbiology.

[42]  S. Dowd,et al.  Comparison of Culture and Molecular Identification of Bacteria in Chronic Wounds , 2012, International journal of molecular sciences.

[43]  Susan Murray,et al.  Decade-long bacterial community dynamics in cystic fibrosis airways , 2012, Proceedings of the National Academy of Sciences.

[44]  N. Pace,et al.  Microbial ecology and evolution: a ribosomal RNA approach. , 1986, Annual review of microbiology.

[45]  Sean R. Eddy,et al.  Infernal 1.0: inference of RNA alignments , 2009, Bioinform..

[46]  Lauren M. Bragg,et al.  Fast, accurate error-correction of amplicon pyrosequences using Acacia , 2012, Nature Methods.

[47]  F. Bushman,et al.  Sampling and pyrosequencing methods for characterizing bacterial communities in the human gut using 16S sequence tags , 2010, BMC Microbiology.

[48]  Euan A Ashley,et al.  Performance comparison of whole-genome sequencing platforms , 2011, Nature Biotechnology.

[49]  D. Thaler,et al.  Optimizing Taq Polymerase Concentration for Improved Signal-to-Noise in the Broad Range Detection of Low Abundance Bacteria , 2009, PloS one.

[50]  Russell J. Davenport,et al.  Removing Noise From Pyrosequenced Amplicons , 2011, BMC Bioinformatics.

[51]  G. Rogers,et al.  Studying bacterial infections through culture-independent approaches. , 2009, Journal of medical microbiology.

[52]  D. Raoult,et al.  Metagenomic analysis of brain abscesses identifies specific bacterial associations. , 2012, Clinical infectious diseases : an official publication of the Infectious Diseases Society of America.

[53]  S. Vermeire,et al.  Coamplification of Eukaryotic DNA with 16S rRNA Gene-Based PCR Primers: Possible Consequences for Population Fingerprinting of Complex Microbial Communities , 2008, Current Microbiology.

[54]  Robert C. Edgar,et al.  BIOINFORMATICS APPLICATIONS NOTE , 2001 .

[55]  Robert C. Edgar,et al.  MUSCLE: multiple sequence alignment with high accuracy and high throughput. , 2004, Nucleic acids research.

[56]  A. Matlow,et al.  Comparison of Culture and PCR for Detection of Burkholderia cepacia in Sputum Samples of Patients with Cystic Fibrosis , 1998, Journal of Clinical Microbiology.

[57]  T. Dallman,et al.  Performance comparison of benchtop high-throughput sequencing platforms , 2012, Nature Biotechnology.

[58]  Anthony O'Donnell,et al.  Microbial 16S rRNA Ion Tag and community metagenome sequencing using the Ion Torrent (PGM) Platform. , 2012, Journal of microbiological methods.

[59]  M. Wolfgang,et al.  Detection of anaerobic bacteria in high numbers in sputum from patients with cystic fibrosis. , 2008, American journal of respiratory and critical care medicine.

[60]  S. Lory,et al.  Phylogenetic and metabolic diversity of bacteria associated with cystic fibrosis , 2011, The ISME Journal.

[61]  M. Pop,et al.  Metagenomic Analysis of the Human Distal Gut Microbiome , 2006, Science.

[62]  S. Leroy,et al.  The Airway Microbiota in Cystic Fibrosis: A Complex Fungal and Bacterial Community—Implications for Therapeutic Management , 2012, PloS one.

[63]  T. Kudo,et al.  Evaluation of primers and PCR conditions for the analysis of 16S rRNA genes from a natural environment. , 2003, FEMS Microbiology Letters.

[64]  G. Wang,et al.  The frequency of chimeric molecules as a consequence of PCR co-amplification of 16S rRNA genes from different bacterial species. , 1996, Microbiology.

[65]  Fiona Powrie,et al.  Microbiota, Disease, and Back to Health: A Metastable Journey , 2012, Science Translational Medicine.

[66]  B. Haas,et al.  Chimeric 16S rRNA sequence formation and detection in Sanger and 454-pyrosequenced PCR amplicons. , 2011, Genome research.

[67]  Paramvir S. Dehal,et al.  FastTree 2 – Approximately Maximum-Likelihood Trees for Large Alignments , 2010, PloS one.