Skip to main content

Main menu

  • Home
  • About
    • Clinical Chemistry
    • Editorial Board
    • Most Read
    • Most Cited
    • Alerts
    • CE Credits
  • Articles
    • Current Issue
    • Early Release
    • Future Table of Contents
    • Archive
    • Browse by Subject
  • Info for
    • Authors
    • Reviewers
    • Subscribers
    • Advertisers
    • Permissions & Reprints
  • Resources
    • AACC Learning Lab
    • Clinical Chemistry Trainee Council
    • Clinical Case Studies
    • Clinical Chemistry Guide to Scientific Writing
    • Clinical Chemistry Guide to Manuscript Review
    • Journal Club
    • Podcasts
    • Q&A
    • Translated Content
  • Abstracts
  • Submit
  • Contact
  • Other Publications
    • The Journal of Applied Laboratory Medicine

User menu

  • Subscribe
  • My alerts
  • Log in

Search

  • Advanced search
Clinical Chemistry
  • Other Publications
    • The Journal of Applied Laboratory Medicine
  • Subscribe
  • My alerts
  • Log in
Clinical Chemistry

Advanced Search

  • Home
  • About
    • Clinical Chemistry
    • Editorial Board
    • Most Read
    • Most Cited
    • Alerts
    • CE Credits
  • Articles
    • Current Issue
    • Early Release
    • Future Table of Contents
    • Archive
    • Browse by Subject
  • Info for
    • Authors
    • Reviewers
    • Subscribers
    • Advertisers
    • Permissions & Reprints
  • Resources
    • AACC Learning Lab
    • Clinical Chemistry Trainee Council
    • Clinical Case Studies
    • Clinical Chemistry Guide to Scientific Writing
    • Clinical Chemistry Guide to Manuscript Review
    • Journal Club
    • Podcasts
    • Q&A
    • Translated Content
  • Abstracts
  • Submit
  • Contact
Research ArticleArticle

Clinical Next Generation Sequencing Outperforms Standard Microbiological Culture for Characterizing Polymicrobial Samples

Lisa A. Cummings, Kyoko Kurosawa, Daniel R. Hoogestraat, Dhruba J. SenGupta, Febrian Candra, Michael Doyle, Stephanie Thielges, Tyler A. Land, Christopher A. Rosenthal, Noah G. Hoffman, Stephen J. Salipante, Brad T. Cookson
DOI: 10.1373/clinchem.2016.258806 Published October 2016
Lisa A. Cummings
Departments of Laboratory Medicine and
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Kyoko Kurosawa
Departments of Laboratory Medicine and
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Daniel R. Hoogestraat
Departments of Laboratory Medicine and
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Dhruba J. SenGupta
Departments of Laboratory Medicine and
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Febrian Candra
Departments of Laboratory Medicine and
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Michael Doyle
Departments of Laboratory Medicine and
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Stephanie Thielges
Departments of Laboratory Medicine and
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Tyler A. Land
Departments of Laboratory Medicine and
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Christopher A. Rosenthal
Departments of Laboratory Medicine and
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Noah G. Hoffman
Departments of Laboratory Medicine and
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Stephen J. Salipante
Departments of Laboratory Medicine and
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Brad T. Cookson
Departments of Laboratory Medicine and Microbiology, University of Washington, Seattle, WA.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: cookson@uw.edu
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

Abstract

BACKGROUND: Humans suffer from infections caused by single species or more complex polymicrobial communities. Identification of infectious bacteria commonly employs microbiological culture, which depends upon the in vitro propagation and isolation of viable organisms. In contrast, detection of bacterial DNA using next generation sequencing (NGS) allows culture-independent microbial profiling, potentially providing important new insights into the microbiota in clinical specimens.

METHODS: NGS 16S rRNA gene sequencing (NGS16S) was compared with culture using (a) synthetic polymicrobial samples for which the identity and abundance of organisms present were precisely defined and (b) primary clinical specimens.

RESULTS: Complex mixtures of at least 20 organisms were well resolved by NGS16S with excellent reproducibility. In mixed bacterial suspensions (107 total genomes), we observed linear detection of a target organism over a 4-log concentration range (500–3 × 106 genomes). NGS16S analysis more accurately recapitulated the known composition of synthetic samples than standard microbiological culture using nonselective media, which distorted the relative abundance of organisms and frequently failed to identify low-abundance pathogens. However, extended quantitative culture using selective media for each of the component species recovered the expected organisms at the proper abundance, validating NGS16S results. In an analysis of sputa from cystic fibrosis patients, NGS16S identified more clinically relevant pathogens than standard culture.

CONCLUSIONS: Biases in standard, nonselective microbiological culture lead to a distorted characterization of polymicrobial mixtures. NGS16S demonstrates enhanced reproducibility, quantification, and classification accuracy compared with standard culture, providing a more comprehensive, accurate, and culture-free analysis of clinical specimens.

Clinical management of infection relies on early detection and identification of the pathogen(s) involved to facilitate appropriate, targeted therapy. The recovery of individual organisms in culture allows for biochemical and phenotypic identification as well as antimicrobial susceptibility testing; culture therefore remains the mainstay of the current clinical microbiology laboratory. However, clinically relevant organisms may be difficult to culture, including fastidious and/or slow-growing bacteria that are vulnerable to improper sample storage, transport, or culture conditions (1–4). Additionally, culture-based identification may take multiple days or weeks, and there are practical limits to the number of different species within a complex mixture that can be isolated and identified using culture-based techniques (5).

Standard molecular methods, in which target genes are PCR-amplified from patient specimens for conventional Sanger sequencing, do not depend on in vitro growth for bacterial detection and have become routine in many clinical microbiology laboratories (6). However, this approach is not readily applied to polymicrobial specimens, where the presence of multiple templates results in superimposed Sanger electropherograms that are generally uninterpretable (5). More recently, next generation sequencing (NGS)3 has been used to overcome this limitation: because millions of DNA molecules are sequenced independently and in parallel, each can be classified individually. With the advent of rapid and affordable bench-top DNA sequencers, NGS is becoming an attractive extension of Sanger sequencing methods for culture-independent microbial profiling in clinical microbiology laboratories.

Although NGS analysis offers theoretical advantages over standard microbiological culture, rigorous comparisons of the 2 methods have not been performed. Defined mixtures of genomic DNA are commercially available, but cannot be used to evaluate culture performance. Comparisons using patient samples are complicated by the fact that the true identity and relative abundance of bacteria present are unknown. In this study, we evaluated NGS and microbiological culture using both clinical specimens and synthetic polymicrobial test samples created by combining bacteria in precisely defined ratios.

Materials and Methods

SAMPLES

Use of clinical microbiological specimens was approved by the University of Washington Human Subjects Review Board (approval number 42541) and was conducted in accordance with the Declaration of Helsinki. As a minimal risk study utilizing surplus microbiological isolates, consent was not required. Control amplicon mixtures (see Supplemental Fig. 1A in the Data Supplement that accompanies the online version of this article at http://www.clinchem.org/content/vol62/issue11) were generated by pooling equal amounts of individual amplicon sequencing libraries. Commercial control mixtures (see online Supplemental Fig. 1B) containing genomic DNA from 20 bacterial species at defined 16S rRNA operon counts (105 copies/μL each, or 103–106 copies/μL), were obtained from BEI Resources (Mock Community B). Synthetic polymicrobial test samples (see online Supplemental Fig. 1C) were created by growing bacteria separately in BHI broth with vitamin K and hemin (Remel R07019), precisely measuring the bacterial concentration for each culture using a Multisizer 4 cell counter (Beckman Coulter), and combining in defined ratios. DNA was extracted from synthetic test samples using the UltraClean microbial DNA isolation kit (MOBIO) and from human-derived specimens using a High Pure PCR template preparation kit (Roche) with mechanical disruption of samples with 1.4-mm ceramic beads followed by enzymatic lysis via proteinase K.

LIBRARY PREPARATION, SEQUENCING, AND DATA ANALYSIS

Sequencing libraries were prepared as previously described (7); a detailed description is provided in the online Supplemental Methods file. Briefly, the 16S rRNA V1–V2 region was amplified using custom primers incorporating Illumina-compatible sequencing adaptors and a sample-specific 8-bp barcode sequence. Paired-end sequencing was performed on an Illumina MiSeq using a 500-cycle sequencing kit (version 2) with custom primers (7) and in the presence of bacterial whole-genome shotgun sequencing libraries or 7% PhiX control (Illumina). Raw data processing and run demultiplexing were performed using on-instrument software. Sequences were processed and classified as previously described (5, 8). Where appropriate, the relative abundance of each classification was corrected for 16S ribosomal RNA operon copy numbers (9) using the NCBI taxonomy database (10) and the Ribosomal RNA Database project (RDP; release 10, update 32), and was expressed as percent corrected reads. For each taxonomic node where 16S rRNA copy numbers were not available, copy numbers were inferred in 1 of 2 ways. First, for taxonomic nodes where 16S copy numbers were available in at least 1 descendent taxonomic classification with lower-order rank (i.e., if genus “Escherichia” had an unknown copy number but species “Escherichia coli” had a known copy number), that node was designated the median copy number of all immediate descendent classifications. This process was repeated recursively at each taxonomic level to the root of the taxonomy. Secondly, and conversely, for taxonomic classifications where 16S rRNA copy numbers of a lower-order rank were not available (i.e., if species “Escherichia coli” had an unknown copy number but genus “Escherichia” had a known copy number) each descendent node was assigned the 16S rRNA copy number of the immediate parent node. This process was repeated recursively, working out from the root of the taxonomy, until every taxonomic node was assigned a copy number.

MICROBIOLOGICAL CULTURE

Microbiological culture was performed by the University of Washington Clinical Microbiology Laboratory, according to standard clinical procedures. Sputum samples were plated on sheep blood, MacConkey, chocolate, mannitol salt, and cepacia agar plates as previously described (5). Synthetic polymicrobial test samples were plated on sheep blood, MacConkey, and chocolate agar plates. For extended quantitative culture of synthetic test samples, 100 μL of various dilutions in PBS were plated onto Hektoen agar (Remel 453572, to isolate Salmonella), mannitol salt agar (Remel 453902, to isolate Staphylococcus), and LB agar containing either 50 μg/mL chloramphenicol (to isolate Pseudomonas) or 40 μg/mL gentamycin (to isolate Achromobacter). Colonies on dilution plates containing 30–300 colonies were counted and used to calculate colony-forming units per milliliter.

Results

ACCURACY OF NGS16S FOR INTERROGATING POLYMICROBIAL SAMPLES

To assess the degree of bias inherent to 16S rRNA NGS analysis (NGS16S) of complex mixtures, we analyzed control mixtures in which the identity and relative abundance of all organism templates was known. We began by sequencing a sample containing equimolar amounts of individual amplicon libraries (see online Supplemental Fig. 1A). In this sample, DNA extraction or PCR biases were eliminated by quantifying and pooling separate products downstream from these steps. Reads corresponding to all organisms were recovered near the expected abundance of 5.5% per organism (Fig. 1A, mean 5.5%, range 4.4%–6.4%), indicating minimal sequencing bias. We next sequenced 2 control mixtures comprised of genomic DNA individually extracted from 20 different organisms and combined in defined 16S rRNA operon counts for each organism (see online Supplemental Fig. 1B). For these samples, extraction biases were eliminated by the pooling strategy, and any differences from the expected values were largely the result of PCR biases introduced during library preparation. When16S templates for each organism were present in equal concentration, all organisms were identified by NGS16S analysis (Fig. 1B). However, the relative abundances of many organisms differed up to 2-fold from the expected value of 5%, indicating organism-specific PCR biases. Despite this observation, all organisms were detected in the sample where the highest and lowest template concentration differed by 4 orders of magnitude (Fig. 1C). This demonstrated that the presence of abundant templates did not obscure the detection of rare templates in the same sample. Importantly, the relative abundance of organisms inferred by NGS16S mirrored the known abundance of templates in this sample (Fig. 1C), indicating its utility as a semiquantitative assay.

Fig. 1.
  • Download figure
  • Open in new tab
Fig. 1. Deconvolution of polymicrobial samples using NGS16S.

A control mixture of pooled amplicons (A) or commercial control genomic DNA mixtures containing defined 16S rRNA operon counts [105copies/μL each (B), or 103–106 copies/μL (C)] were sequenced. Mean read counts are provided for organisms of abundance too low to be visualized. Refer to online Supplemental Fig. 1, A and B for illustrations of each sample type.

REPRODUCIBILITY OF NGS16S ANALYSIS

To test NGS16S reproducibility, we sequenced DNA extracted from a single clinical brain abscess specimen for which multiple templates were detected by conventional Sanger sequencing. Serial dilutions of this DNA were used as template for library preparation and the greatest dilution of DNA that yielded a visible product (1:1000, data not shown) was identified as the detection limit for this specimen. Separate dilutions at above (1:500) and below (1:5000) this limit were subsequently used as templates for duplicate library preparation on separate days by different laboratory personnel. Excellent inter- and intraassay reproducibility was observed (Fig. 2); statistically indistinguishable relative abundance assignments (1-way ANOVA, P <0.05) were generated for each classification reported. We concluded that NGS16S analysis was highly reproducible over a wide range of sample concentrations.

Fig. 2.
  • Download figure
  • Open in new tab
Fig. 2. NGS16S analysis is highly reproducible.

NGS16S libraries were generated from dilutions of the same clinical sample on 3 separate days by different laboratory personnel. Replicates from each day are presented side by side. Classifications at genus level or higher are listed as “Other.” Refer to online Supplemental Table 1 for summary statistics.

DETECTION OF A TARGET ORGANISM IN A DEFINED POLYMICROBIAL SAMPLE

To better define the limits of NGS16S analysis, we created and analyzed defined polymicrobial test mixtures of intact cells. We precisely quantified bacteria in separate cultures of 4 different organisms using a Coulter counter (which accurately counts individual particles in solution), combined 3 of those organisms in defined proportions, and then serially diluted the final organism, Salmonella typhimurium, into the mixture (see online Supplemental Fig. 1C). Thus, the relative abundance of 3 organisms was held constant, while a fourth varied over a wide range. DNA was extracted from each dilution and subjected to NGS16S analysis using 107 total genomes as template. As few as 500 Salmonella (0.005% of the total bacteria) could be detected by NGS16S (Fig. 3A). The recovery of Salmonella sequences was linear over 4 orders of magnitude (Fig. 3B, r2 = 0.86), although the relative abundance of Salmonella by NGS was slightly lower than the known input value. This reflects cumulative bias in DNA extraction, PCR amplification and NGS sequencing (11–18).

Fig. 3.
  • Download figure
  • Open in new tab
Fig. 3. Recovery of target sequences in a defined polymicrobial sample.

(A), S. typhimurium (SALM) was serially diluted into a mixture containing 75% A. xylosoxidans (AXYL), 24% S. aureus (SAUR), and 1% P. aeruginosa (PSAR). (B), The observed abundance of Salmonella detected by NGS (% corrected reads) compared to the known input (% input bacteria). Refer to online Supplemental Fig. 1C for illustration of sample type.

COMPARISON OF NGS16S AND CULTURE FOR ANALYSIS OF POLYMICROBIAL SPECIMENS

To directly compare NGS16S and standard culture, we again used defined polymicrobial test mixtures in which constituent organisms were precisely quantified (see online Supplemental Fig. 1C). Mixtures were divided into multiple aliquots and subjected to (a) DNA extraction and NGS16S analysis, (b) standard clinical culture and scored according to conventional 1+ through 4+ quantification, and (c) extended quantitative culture in which multiple sample dilutions were plated onto various selective media to allow isolation and colony-forming unit enumeration of each of the 4 constituent organisms independently. Studies were performed in quadruplicate and benchmarked against known organism abundance as established using Coulter counting (Fig. 4, A and B). For each replicate, NGS16S properly recapitulated the known composition of the mixtures, including organisms present at 1% relative abundance (Fig. 4, C and D). In contrast, standard culture, where organisms are plated at high density on rich media and must directly compete with one another, was less accurate (Fig. 4, E and F). Staphylococcus aureus and S. typhimurium were universally reported as 4+ (in some cases a higher quantification than the most abundant organism), although they represented no more than 20% of the overall mixture. Organisms at 1% relative abundance were either not detected (in 5 of 8 replicates), or were reported at amounts equal to or higher than the most abundant organism (3 of 8 replicates).

Fig. 4.
  • Download figure
  • Open in new tab
Fig. 4. Comparison of NGS16S and culture for accurate analysis of defined polymicrobial samples.

Bacteria were combined in various proportions to generate polymicrobial test samples of known input composition (A, B), which were subjected to NGS16S analysis (C, D), standard clinical culture (E, F) and quantitative culture on selective medium (G, H). Relative abundance of recovered component organisms (output) is expressed as true percent of input cells by Coulter counting (A, B), percent corrected reads (C, D), 1+, 2+, 3+, or 4+ (E, F), and percent of total recovered colonies by quantitative culture on selective medium (G, H). Mean output (% of total) is provided for organisms of abundance too low to be visualized. PSAR, P. aeruginosa; SALM, S. typhimurium; SAUR, S. aureus; AXYL, A. xylosoxidans; NR, not recovered. Refer to online Supplemental Fig. 1C for illustration of sample type.

Extended quantitative culture, using various selective media to support the growth of a single constituent and to eliminate direct competition from the other constituents, more accurately reflected actual organism abundance than did standard culture (Fig. 4, G and H). The most and least abundant organisms were correctly and reliably identified, and the relative abundance of all organisms mirrored known proportions in the input mixture. These results validated the accuracy of NGS16S analysis and highlight previously underappreciated biases of standard culture practices.

CHARACTERIZATION OF PRIMARY CLINICAL SPECIMENS

To compare the methods in clinical context, we evaluated the ability of NGS16S and culture to identify select bacterial pathogens of clinical significance from cystic fibrosis (CF) sputa samples (see online Supplemental Table 2). A cohort of 15 patient sputa were processed by the clinical laboratory using standard workflows, and DNA was extracted from residual material and subjected to NGS16S analysis (online Supplemental Table 3 provides a comprehensive list of all organisms identified in each sample).

Six of the pathogens of interest were detected in 1 or more cases (Table 1), and all CF pathogens identified by culture were also identified by NGS16S. However, in 3 samples, CF pathogens were identified only by NGS16S (Tables 1 and 2): Streptococcus agalactiae (specimen CF2, 9.55% of corrected reads), Achromobacter xylosoxidans (specimen CF4, 10.26% of corrected reads), and Burkholderia cepacia complex (specimen CF15, 4.48% of corrected reads). The percentage of reads for each of these organisms was well above the detection limit identified in previous experiments with defined samples (0.005%, Fig. 3). It is unlikely that these organisms were the result of reagent or sample contamination since corresponding reads were not observed in any paired negative control sample.

View this table:
  • View inline
  • View popup
Table 1.

Select CF pathogens identified by microbiological culture and NGS16S.a

View this table:
  • View inline
  • View popup
Table 2.

Comparison of culture and NGS16S results for select samples.a

We also compared the quantification of the 2 assays in this context (Table 2). Pseudomonas aeruginosa was the most commonly identified pathogen in the cohort (12 of 15 samples), and its relative abundance by culture was consistent with the quantification by NGS16S in most cases (Table 2). However in 3 samples where P. aeruginosa was reported as 2+ or less by culture, it ranged from 0.47% (specimen CF3) to 90% (specimens CF12, CF13) relative abundance by NGS16S read count. S. aureus, another frequently identified CF pathogen (6 of 16 samples), was consistently overestimated by culture compared to NGS16S; although never more than 3.53% relative abundance by NGS, it was quantified in culture over a range of 1 colony to 4+ abundance (Table 2; also see online Supplemental Fig. 2).

In addition to the identification of known pathogens, NGS16S also provided more information than culture about other organisms in each sample. There were 33 species-level identifications by NGS16S compared to only 10 for standard microbiological culture (see online Supplemental Table 3). Whereas clinical reporting described “contaminating oral-pharyngeal microbiota,” NGS16S was able to resolve the identities and relative abundance of these organisms (see online Supplemental Table 3). Importantly, anaerobes were found in a majority (13 of 15) of these samples by NGS16S, but not identified using standard clinical culture. The most common organisms were Veillonella (8 of 15 samples) and Prevotella (5 of 15 samples).

Finally, culture identified some organisms that were missed by NGS16S (Table 2). 16S rRNA gene sequence analysis is restricted to detection and identification of bacteria, and the majority of discrepancies between the 2 methods (7 of 8) involved yeast or fungi. In only a single case was a bacterium identified by culture but not NGS16S. In CF14, Ewingella americana, was reported in low abundance. The clinical significance of this organism is unclear; it is not known to be part of the normal respiratory microbiota, and can survive in water (19), suggesting that it could be an environmental contaminant.

Discussion

The ease, speed, and decreasing cost of NGS sequencing make the application of this technology increasingly attractive for routine clinical use. Here, we demonstrate that NGS sequencing generates a reproducible (Fig. 2), analytically sensitive (Fig. 3), and accurate assessment of the identity and relative abundance of organisms present in polymicrobial samples, outperforming standard culture (Fig. 4). Our work indicates that the interpretation of standard culture results, where multiple organisms with differing growth requirements compete for resources on rich, nonselective media, is substantially more complicated and biased than commonly presumed.

Culture-based methods are the mainstay of specimen analysis in the clinical laboratory, and the isolation of individual organisms is required for antimicrobial susceptibility testing, strain typing and virulence studies. For an organism to be identified in culture, it must be both viable in situ and culturable in vitro. Individual colonies must also be physically isolated and differentiated from others on a culture plate. Yet, viable organisms may fail to grow in vitro due to prior broad-spectrum antibiotic treatment, improper sample storage, transport, or culture conditions. During culture, there is underappreciated amplification bias: a single organism must replicate to approximately 109 organisms to form a visible colony on a plate. Because growth of a particular organism may also be affected by the presence of other organisms in the sample (1, 20–26), known pathogens in a polymicrobial specimen can be obscured or over-estimated (Fig. 4, E and F). Therefore, standard culture is unlikely to reproducibly provide a complete understanding of the microbial composition of a sample containing a complex bacterial community (27).

NGS technology eliminates the requirements of viability, cultivability and physical isolation for bacterial detection and identification from patient specimens. Despite the potential biases of NGS-based methods, introduced by differential DNA extraction and/or PCR amplification as a result of template concentration and complexity (11–18), our results from defined specimens (Fig. 4), and complex DNA mixtures (Fig. 1), demonstrate that NGS16S is superior to standard culture for accurate detection and relative abundance assignments of specimen constituents in complex mixtures. In applications where bacterial load may be of importance NGS16S could be supplemented with qPCR (28). Extended quantitative culture on selective media, each of which supports the growth of only 1 of the constituent organisms present in a polymicrobial test sample (Fig. 4, G and H), both validated our NGS16S results and confirmed the existence of bias introduced in standard culture practices where multiple organisms may be in direct competition. In our analysis of CF sputa, standard culture consistently overestimated the abundance of S. aureus, suggesting a recovery bias in favor of this fast-growing organism with a clearly distinguishable colony phenotype (Fig. 4, E and F, online Supplemental Fig. 2). Although standard culture may be effective for identification of rapidly growing organisms, it may fail to identify other clinically important community members without extensive effort, and we have found that it is unreliable for determining relative abundance.

This study also highlights other aspects of NGS16S that improve upon existing practice. First, although commonly thought of as single agents, pathogens often exist in a polymicrobial context in vivo. Clinically relevant organisms may not directly cause disease, but nevertheless contribute to pathogenesis by altering the microenvironment to facilitate colonization or virulence gene expression of other bacteria (1, 8, 24, 25, 29–31). Therefore, knowledge of the identity and relative abundance of organisms detected in a clinical sample may have important implications for patient care. For example, accurate and complete identification of all organisms present in a sample (particularly unexpected organisms) will facilitate targeted antimicrobial therapy, reducing broad-spectrum antibiotic use (32, 33). Similarly, the metabiome profiles associated with various disease states may provide diagnostic or prognostic information (1, 4, 5, 34–37). In our evaluation of CF sputa, NGS16S analysis provided a more comprehensive assessment of sample microbiota than did standard culture. As a case in point, 3 pathogens of clinical interest, S. agalactiae, A. xylosoxidans, and B. cepacia complex, were identified at significant levels from NGS16S analysis but not recovered in culture (Table 2). Secondly, NGS16S analysis resolved the identities of organisms commonly reported by clinical laboratories as “contaminating oral-pharyngeal microbiota” (see online Supplemental Table 3). Interestingly, anaerobes were found in a majority of samples, the most common being Veillonella and Prevotella, in agreement with published data (4, 35, 37–40). Although not traditionally part of a standard CF culture workup, anaerobes may play a more important role in CF pathology than has been previously appreciated. Anaerobes can be recovered from CF sputa at equivalent numbers to aerobes (37), and recent studies indicate that production of inflammatory short-chain fatty acids (30) or extended-spectrum β-lactamases (31) by anaerobes may contribute to CF disease.

These data demonstrate that NGS-based analysis is redefining the catalog of potential pathogens in human disease states and expanding our view of polymicrobial contributions to the pathogenesis of infection. Moreover, these studies provide a context for the appropriate evaluation of the accuracy of NGS analysis of patient samples when compared to culture based methods. By performing carefully controlled experiments such as those presented here, we can better appreciate the limitations of both NGS16S and culture-based approaches and understand how these complementary methods can best be used in clinical practice.

Footnotes

  • Disclaimer: The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.

  • ↵3 Nonstandard abbreviations:

    NGS,
    next generation sequencing;
    NGS16S,
    NGS of the 16S rRNA gene;
    CF,
    cystic fibrosis.

  • Author Contributions: All authors confirmed they have contributed to the intellectual content of this paper and have met the following 3 requirements: (a) significant contributions to the conception and design, acquisition of data, or analysis and interpretation of data; (b) drafting or revising the article for intellectual content; and (c) final approval of the published article.

  • Authors' Disclosures or Potential Conflicts of Interest: Upon manuscript submission, all authors completed the author disclosure form. Disclosures and/or potential conflicts of interest:

  • Employment or Leadership: None declared.

  • Consultant or Advisory Role: None declared.

  • Stock Ownership: None declared.

  • Honoraria: None declared.

  • Research Funding: National Center for Advancing Translational Sciences of the National Institutes of Health (grant number UL1TR000423).

  • Expert Testimony: None declared.

  • Patents: None declared.

  • Role of Sponsor: The funding organizations played no role in the design of study, choice of enrolled patients, review and interpretation of data, and final approval of manuscript. See feedback in article file for further information.

  • Received for publication April 7, 2016.
  • Accepted for publication August 17, 2016.
  • © 2016 American Association for Clinical Chemistry

References

  1. 1.↵
    1. Kommedal Ø,
    2. Wilhelmsen MT,
    3. Skrede S,
    4. Meisal R,
    5. Jakovljev A,
    6. Gaustad P,
    7. et al
    . Massive parallel sequencing provides new perspectives on bacterial brain abscesses. J Clin Microbiol 2014;52:1990–7.
    OpenUrlAbstract/FREE Full Text
  2. 2.
    1. Willner D,
    2. Daly J,
    3. Whiley D,
    4. Grimwood K,
    5. Wainwright CE,
    6. Hugenholtz P
    . Comparison of DNA extraction methods for microbial community profiling with an application to pediatric bronchoalveolar lavage samples. PLoS One 2012;7:e34605.
    OpenUrlCrossRefPubMed Order article via Infotrieve
  3. 3.
    1. Tanaka T,
    2. Kawasaki K,
    3. Daimon S,
    4. Kitagawa W,
    5. Yamamoto K,
    6. Tamaki H,
    7. et al
    . A hidden pitfall in the preparation of agar media undermines microorganism cultivability. Appl Environ Microbiol 2014;80:7659–66.
    OpenUrlAbstract/FREE Full Text
  4. 4.↵
    1. Lambiase A,
    2. Catania MR,
    3. Rossano F
    . Anaerobic bacteria infection in cystic fibrosis airway disease. New Microbiol 2010;33:185–94.
    OpenUrlPubMed Order article via Infotrieve
  5. 5.↵
    1. Salipante SJ,
    2. Sengupta DJ,
    3. Rosenthal C,
    4. Costa G,
    5. Spangler J,
    6. Sims EH,
    7. et al
    . Rapid 16s rRNA next-generation sequencing of polymicrobial clinical samples for diagnosis of complex bacterial infections. PLoS One 2013;8:e65226.
    OpenUrlCrossRefPubMed Order article via Infotrieve
  6. 6.↵
    1. Clarridge JE
    . Impact of 16s rRNA gene sequence analysis for identification of bacteria on clinical microbiology and infectious diseases. Clin Microbiol Rev 2004;17:840–62.
    OpenUrlAbstract/FREE Full Text
  7. 7.↵
    1. Salipante SJ,
    2. Kawashima T,
    3. Rosenthal C,
    4. Hoogestraat DR,
    5. Cummings LA,
    6. Sengupta DJ,
    7. et al
    . Performance comparison of Illumina and Ion Torrent next-generation sequencing platforms for 16s rRNA-based bacterial community profiling. Appl Environ Microbiol 2014;80:7583–91.
    OpenUrlAbstract/FREE Full Text
  8. 8.↵
    1. Salipante SJ,
    2. Hoogestraat DR,
    3. Abbott AN,
    4. Sengupta DJ,
    5. Cummings LA,
    6. Butler-Wu SM,
    7. et al
    . Coinfection of Fusobacterium nucleatum and Actinomyces israelii in mastoiditis diagnosed by next-generation DNA sequencing. J Clin Microbiol 2014;52:1789–92.
    OpenUrlAbstract/FREE Full Text
  9. 9.↵
    1. Kembel SW,
    2. Wu M,
    3. Eisen JA,
    4. Green JL
    . Incorporating 16s gene copy number information improves estimates of microbial diversity and abundance. PLoS Comput Biol 2012;8:e1002743.
    OpenUrlCrossRefPubMed Order article via Infotrieve
  10. 10.↵
    National Center for Biotechnology Information. NCBI taxonomy database. http://www.ncbi.nlm.nih.gov/taxonomy (accessed August 2016).
  11. 11.↵
    1. van Dijk EL,
    2. Jaszczyszyn Y,
    3. Thermes C
    . Library preparation methods for next-generation sequencing: tone down the bias. Exp Cell Res 2014;322:12–20.
    OpenUrlCrossRefPubMed Order article via Infotrieve
  12. 12.
    1. Zhao J,
    2. Carmody LA,
    3. Kalikin LM,
    4. Li J,
    5. Petrosino JF,
    6. Schloss PD,
    7. et al
    . Impact of enhanced Staphylococcus DNA extraction on microbial community measures in cystic fibrosis sputum. PLoS One 2012;7:e33127.
    OpenUrlCrossRefPubMed Order article via Infotrieve
  13. 13.
    1. Abusleme L,
    2. Hong BY,
    3. Dupuy AK,
    4. Strausbaugh LD,
    5. Diaz PI
    . Influence of DNA extraction on oral microbial profiles obtained via 16s rRNA gene sequencing. J Oral Microbiol 2014;6.
  14. 14.
    1. Wu JY,
    2. Jiang XT,
    3. Jiang YX,
    4. Lu SY,
    5. Zou F,
    6. Zhou HW
    . Effects of polymerase, template dilution and cycle number on PCR based 16 s rRNA diversity analysis using the deep sequencing method. BMC Microbiol 2010;10:255.
    OpenUrlCrossRefPubMed Order article via Infotrieve
  15. 15.
    1. Engelbrektson A,
    2. Kunin V,
    3. Wrighton KC,
    4. Zvenigorodsky N,
    5. Chen F,
    6. Ochman H,
    7. Hugenholtz P
    . Experimental factors affecting PCR-based estimates of microbial species richness and evenness. ISME J 2010;4:642–7.
    OpenUrlCrossRefPubMed Order article via Infotrieve
  16. 16.
    1. Pinto AJ,
    2. Raskin L
    . PCR biases distort bacterial and archaeal community structure in pyrosequencing datasets. PLoS One 2012;7:e43093.
    OpenUrlCrossRefPubMed Order article via Infotrieve
  17. 17.
    1. Biesbroek G,
    2. Sanders EA,
    3. Roeselers G,
    4. Wang X,
    5. Caspers MP,
    6. Trzciński K,
    7. et al
    . Deep sequencing analyses of low density microbial communities: working at the boundary of accurate microbiota detection. PLoS One 2012;7:e32942.
    OpenUrlCrossRefPubMed Order article via Infotrieve
  18. 18.↵
    1. Kennedy K,
    2. Hall MW,
    3. Lynch MD,
    4. Moreno-Hagelsieb G,
    5. Neufeld JD
    . Evaluating bias of Illumina-based bacterial 16s rRNA gene profiles. Appl Environ Microbiol 2014;80:5717–22.
    OpenUrlAbstract/FREE Full Text
  19. 19.↵
    1. McNeil MM,
    2. Davis BJ,
    3. Solomon SL,
    4. Anderson RL,
    5. Shulman ST,
    6. Gardner S,
    7. et al
    . Ewingella americana: recurrent pseudobacteremia from a persistent environmental reservoir. J Clin Microbiol 1987;25:498–500.
    OpenUrlAbstract/FREE Full Text
  20. 20.↵
    1. Tashiro Y,
    2. Yawata Y,
    3. Toyofuku M,
    4. Uchiyama H,
    5. Nomura N
    . Interspecies interaction between Pseudomonas aeruginosa and other microorganisms. Microbes Environ 2013;28:13–24.
    OpenUrlCrossRefPubMed Order article via Infotrieve
  21. 21.
    1. Costello A,
    2. Reen FJ,
    3. O'Gara F,
    4. Callaghan M,
    5. McClean S
    . Inhibition of co-colonizing cystic fibrosis-associated pathogens by Pseudomonas aeruginosa and Burkholderia multivorans. Microbiology 2014;160:1474–87.
    OpenUrlCrossRefPubMed Order article via Infotrieve
  22. 22.
    1. Patankar CV,
    2. Joshi LM
    . Bacteriocin production in Salmonella. J Postgrad Med 1985;31:46–51.
    OpenUrlPubMed Order article via Infotrieve
  23. 23.
    1. Sizova MV,
    2. Hohmann T,
    3. Hazen A,
    4. Paster BJ,
    5. Halem SR,
    6. Murphy CM,
    7. et al
    . New approaches for isolation of previously uncultivated oral bacteria. Appl Environ Microbiol 2012;78:194–203.
    OpenUrlAbstract/FREE Full Text
  24. 24.↵
    1. Pustelny C,
    2. Komor U,
    3. Pawar V,
    4. Lorenz A,
    5. Bielecka A,
    6. Moter A,
    7. et al
    . Contribution of Veillonella parvula to Pseudomonas aeruginosa-mediated pathogenicity in a murine tumor model system. Infect Immun 2015;83:417–29.
    OpenUrlAbstract/FREE Full Text
  25. 25.↵
    1. McNally L,
    2. Brown SP
    . Building the microbiome in health and disease: niche construction and social conflict in bacteria. Philos Trans R Soc Lond B Biol Sci 2015;370.
  26. 26.↵
    1. Kaeberlein T,
    2. Lewis K,
    3. Epstein SS
    . Isolating “uncultivable” microorganisms in pure culture in a simulated natural environment. Science 2002;296:1127–9.
    OpenUrlAbstract/FREE Full Text
  27. 27.↵
    1. Rogers GB,
    2. Daniels TW,
    3. Tuck A,
    4. Carroll MP,
    5. Connett GJ,
    6. David GJ,
    7. Bruce KD
    . Studying bacteria in respiratory specimens by using conventional and molecular microbiological approaches. BMC Pulm Med 2009;9:14.
    OpenUrlCrossRefPubMed Order article via Infotrieve
  28. 28.↵
    1. Brukner I,
    2. Longtin Y,
    3. Oughton M,
    4. Forgetta V,
    5. Dascal A
    . Assay for estimating total bacterial load: relative qPCR normalisation of bacterial load with associated clinical implications. Diagn Microbiol Infect Dis 2015;83:1–6.
    OpenUrlCrossRefPubMed Order article via Infotrieve
  29. 29.↵
    1. Hajishengallis G,
    2. Darveau RP,
    3. Curtis MA
    . The keystone-pathogen hypothesis. Nat Rev Microbiol 2012;10:717–25.
    OpenUrlCrossRefPubMed Order article via Infotrieve
  30. 30.↵
    1. Mirković B,
    2. Murray MA,
    3. Lavelle GM,
    4. Molloy K,
    5. Azim AA,
    6. Gunaratnam C,
    7. et al
    . The role of short-chain fatty acids, produced by anaerobic bacteria, in the cystic fibrosis airway. Am J Respir Crit Care Med 2015;192:1314–24.
    OpenUrlCrossRefPubMed Order article via Infotrieve
  31. 31.↵
    1. Sherrard LJ,
    2. McGrath SJ,
    3. McIlreavey L,
    4. Hatch J,
    5. Wolfgang MC,
    6. Muhlebach MS,
    7. et al
    . Production of extended-spectrum β-lactamases and the potential indirect pathogenic role of Prevotella isolates from the cystic fibrosis respiratory microbiota. Int J Antimicrob Agents 2016;47:140–5.
    OpenUrlCrossRefPubMed Order article via Infotrieve
  32. 32.↵
    1. Rhoads DD,
    2. Cox SB,
    3. Rees EJ,
    4. Sun Y,
    5. Wolcott RD
    . Clinical identification of bacteria in human chronic wound infections: culturing vs. 16s ribosomal DNA sequencing. BMC Infect Dis 2012;12:321.
    OpenUrlCrossRefPubMed Order article via Infotrieve
  33. 33.↵
    1. Livermore DM,
    2. Wain J
    . Revolutionising bacteriology to improve treatment outcomes and antibiotic stewardship. Infect Chemother 2013;45:1–10.
    OpenUrlCrossRefPubMed Order article via Infotrieve
  34. 34.↵
    1. Zemanick ET,
    2. Sagel SD,
    3. Harris JK
    . The airway microbiome in cystic fibrosis and implications for treatment. Curr Opin Pediatr 2011;23:319–24.
    OpenUrlCrossRefPubMed Order article via Infotrieve
  35. 35.↵
    1. Coburn B,
    2. Wang PW,
    3. Diaz Caballero J,
    4. Clark ST,
    5. Brahma V,
    6. Donaldson S,
    7. et al
    . Lung microbiota across age and disease stage in cystic fibrosis. Sci Rep 2015;5:10241.
    OpenUrlCrossRefPubMed Order article via Infotrieve
  36. 36.
    1. Sprockett DD,
    2. Ammons CG,
    3. Tuttle MS
    . Use of 16s rRNA sequencing and quantitative PCR to correlate venous leg ulcer bacterial bioburden dynamics with wound expansion, antibiotic therapy, and healing. Wound Repair Regen 2015;23:765–71.
    OpenUrlCrossRefPubMed Order article via Infotrieve
  37. 37.↵
    1. Tunney MM,
    2. Field TR,
    3. Moriarty TF,
    4. Patrick S,
    5. Doering G,
    6. Muhlebach MS,
    7. et al
    . Detection of anaerobic bacteria in high numbers in sputum from patients with cystic fibrosis. Am J Respir Crit Care Med 2008;177:995–1001.
    OpenUrlCrossRefPubMed Order article via Infotrieve
  38. 38.
    1. Delhaes L,
    2. Monchy S,
    3. Fréalle E,
    4. Hubans C,
    5. Salleron J,
    6. Leroy S,
    7. et al
    . The airway microbiota in cystic fibrosis: a complex fungal and bacterial community–implications for therapeutic management. PLoS One 2012;7:e36313.
    OpenUrlCrossRefPubMed Order article via Infotrieve
  39. 39.
    1. Field TR,
    2. Sibley CD,
    3. Parkins MD,
    4. Rabin HR,
    5. Surette MG
    . The genus prevotella in cystic fibrosis airways. Anaerobe 2010;16:337–44.
    OpenUrlCrossRefPubMed Order article via Infotrieve
  40. 40.↵
    1. Sibley CD,
    2. Grinwis ME,
    3. Field TR,
    4. Eshaghurshan CS,
    5. Faria MM,
    6. Dowd SE,
    7. et al
    . Culture enriched molecular profiling of the cystic fibrosis airway microbiome. PLoS One 2011;6:e22702.
    OpenUrlCrossRefPubMed Order article via Infotrieve
View Abstract
PreviousNext
Back to top

In this issue

Clinical Chemistry: 62 (11)
Vol. 62, Issue 11
November 2016
  • Table of Contents
  • About the Cover
  • Index by author
  • Table of Contents (PDF)
  • Cover (PDF)
  • Advertising (PDF)
  • Ed Board (PDF)
  • Front Matter (PDF)
  • Audio summary of this issue
Print
Share
Clinical Next Generation Sequencing Outperforms Standard Microbiological Culture for Characterizing Polymicrobial Samples
Lisa A. Cummings, Kyoko Kurosawa, Daniel R. Hoogestraat, Dhruba J. SenGupta, Febrian Candra, Michael Doyle, Stephanie Thielges, Tyler A. Land, Christopher A. Rosenthal, Noah G. Hoffman, Stephen J. Salipante, Brad T. Cookson
Clinical Chemistry Nov 2016, 62 (11) 1465-1473; DOI: 10.1373/clinchem.2016.258806
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
Article Alerts
Sign In to Email Alerts with your Email Address
Citation Tools
Clinical Next Generation Sequencing Outperforms Standard Microbiological Culture for Characterizing Polymicrobial Samples
Lisa A. Cummings, Kyoko Kurosawa, Daniel R. Hoogestraat, Dhruba J. SenGupta, Febrian Candra, Michael Doyle, Stephanie Thielges, Tyler A. Land, Christopher A. Rosenthal, Noah G. Hoffman, Stephen J. Salipante, Brad T. Cookson
Clinical Chemistry Nov 2016, 62 (11) 1465-1473; DOI: 10.1373/clinchem.2016.258806

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero

  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Materials and Methods
    • Results
    • Discussion
    • Footnotes
    • References
  • Figures & Data
  • Info & Metrics
  • PDF

Related Articles

Cited By...

More in this TOC Section

  • Integrated Extreme Real-Time PCR and High-Speed Melting Analysis in 52 to 87 Seconds
  • Species Typing of Nontuberculous Mycobacteria by Use of Deoxyribozyme Sensors
  • Inflammatory Biomarkers and Clinical Judgment in the Emergency Diagnosis of Urgent Abdominal Pain
Show more Articles

Similar Articles

Subjects

  • SUBJECT AREAS
    • Molecular Diagnostics and Genetics

Options

  • Home
  • About
  • Articles
  • Information for Authors
  • Resources
  • Abstracts
  • Submit
  • Contact
  • RSS

Other Publications

  • The Journal of Applied Laboratory Medicine
Footer logo

© 2019 American Association for Clinical Chemistry

Powered by HighWire