A Method for Studying Protistan Diversity Using Massively Parallel Sequencing of V9 Hypervariable Regions of Small-Subunit Ribosomal RNA Genes

Background

Massively parallel pyrosequencing of amplicons from the V6 hypervariable regions of small-subunit (SSU) ribosomal RNA (rRNA) genes is commonly used to assess diversity and richness in bacterial and archaeal populations. Recent advances in pyrosequencing technology provide read lengths of up to 240 nucleotides. Amplicon pyrosequencing can now be applied to longer variable regions of the SSU rRNA gene including the V9 region in eukaryotes.

Methodology/Principal Findings

We present a protocol for the amplicon pyrosequencing of V9 regions for eukaryotic environmental samples for biodiversity inventories and species richness estimation. The International Census of Marine Microbes (ICoMM) and the Microbial Inventory Research Across Diverse Aquatic Long Term Ecological Research Sites (MIRADA-LTERs) projects are already employing this protocol for tag sequencing of eukaryotic samples in a wide diversity of both marine and freshwater environments.

Conclusions/Significance

Massively parallel pyrosequencing of eukaryotic V9 hypervariable regions of SSU rRNA genes provides a means of estimating species richness from deeply-sampled populations and for discovering novel species from the environment.     

Supragingival Microbial Profiles of Permanent and Deciduous Teeth in Children with Mixed Dentition

Objectives

The present study was designed to investigate the microbial profiles of teeth in different locations in mixed-dentition-stage children, and to compare the microbiomes of permanent and deciduous teeth in the same healthy oral cavity.

Methods

Supragingival plaque samples of teeth in various locations—the first permanent molars, deciduous molars, deciduous canines and incisors and permanent incisors—were collected from 20 healthy mixed-dentition-stage children with 10–12 permanent teeth erupted. Plaque DNA was extracted, and the V3–V4 hypervariable region of the bacterial 16S rRNA gene was amplified and subjected to sequencing.

Results

On average, 18,051 high-quality sequences per sample were generated. Permanent tooth sites tended to host more diverse bacterial communities than those of deciduous tooth sites. A total of 12 phyla, 21 classes, 38 orders, 66 families, 74 genera were detected ultimately. Five predominant phyla (Proteobacteria, Firmicutes, Bacteroidetes, Fusobacteria and Actinobacteria) were highly variable among sites. Of 26 genera with a mean relative abundance of >0.1%, 16 showed significant differences in relative abundance among the groups. More than 20% of the total operational taxonomical units were detected only in permanent or deciduous teeth. The variation in the microbial community composition was due mainly to permanent teeth being enriched in Actinomyces and deciduous teeth in Treponema. The core microbiome of supragingival plaque in mixed dentition comprised 19 genera with complex correlationships.

Conclusion

Our results suggest differences in microbial diversity and composition between permanent and deciduous teeth sites in mixed dentition. Moreover, the core microbiome of these sites was determined. These findings enhance our understanding of the development of the native oral microbiota with age.     

Oral microbiome profiles: 16S rRNA pyrosequencing and microarray assay comparison

Objectives

The human oral microbiome is potentially related to diverse health conditions and high-throughput technology provides the possibility of surveying microbial community structure at high resolution. We compared two oral microbiome survey methods: broad-based microbiome identification by 16S rRNA gene sequencing and targeted characterization of microbes by custom DNA microarray.

Methods

Oral wash samples were collected from 20 individuals at Memorial Sloan-Kettering Cancer Center. 16S rRNA gene survey was performed by 454 pyrosequencing of the V3–V5 region (450 bp). Targeted identification by DNA microarray was carried out with the Human Oral Microbe Identification Microarray (HOMIM). Correlations and relative abundance were compared at phylum and genus level, between 16S rRNA sequence read ratio and HOMIM hybridization intensity.

Results

The major phyla, Firmicutes, Proteobacteria, Bacteroidetes, Actinobacteria, and Fusobacteria were identified with high correlation by the two methods (r = 0.70∼0.86). 16S rRNA gene pyrosequencing identified 77 genera and HOMIM identified 49, with 37 genera detected by both methods; more than 98% of classified bacteria were assigned in these 37 genera. Concordance by the two assays (presence/absence) and correlations were high for common genera (Streptococcus, Veillonella, Leptotrichia, Prevotella, and Haemophilus; Correlation = 0.70–0.84).

Conclusion

Microbiome community profiles assessed by 16S rRNA pyrosequencing and HOMIM were highly correlated at the phylum level and, when comparing the more commonly detected taxa, also at the genus level. Both methods are currently suitable for high-throughput epidemiologic investigations relating identified and more common oral microbial taxa to disease risk; yet, pyrosequencing may provide a broader spectrum of taxa identification, a distinct sequence-read record, and greater detection sensitivity.     

Community Analysis of Chronic Wound Bacteria Using 16S rRNA Gene-Based Pyrosequencing: Impact of Diabetes and Antibiotics on Chronic Wound Microbiota

Background

Bacterial colonization is hypothesized to play a pathogenic role in the non-healing state of chronic wounds. We characterized wound bacteria from a cohort of chronic wound patients using a 16S rRNA gene-based pyrosequencing approach and assessed the impact of diabetes and antibiotics on chronic wound microbiota.

Methodology/Principal Findings

We prospectively enrolled 24 patients at a referral wound center in Baltimore, MD; sampled patients'' wounds by curette; cultured samples under aerobic and anaerobic conditions; and pyrosequenced the 16S rRNA V3 hypervariable region. The 16S rRNA gene-based analyses revealed an average of 10 different bacterial families in wounds—approximately 4 times more than estimated by culture-based analyses. Fastidious anaerobic bacteria belonging to the Clostridiales family XI were among the most prevalent bacteria identified exclusively by 16S rRNA gene-based analyses. Community-scale analyses showed that wound microbiota from antibiotic treated patients were significantly different from untreated patients (p = 0.007) and were characterized by increased Pseudomonadaceae abundance. These analyses also revealed that antibiotic use was associated with decreased Streptococcaceae among diabetics and that Streptococcaceae was more abundant among diabetics as compared to non-diabetics.

Conclusions/Significance

The 16S rRNA gene-based analyses revealed complex bacterial communities including anaerobic bacteria that may play causative roles in the non-healing state of some chronic wounds. Our data suggest that antimicrobial therapy alters community structure—reducing some bacteria while selecting for others.     

The Characterization of Novel Tissue Microbiota Using an Optimized 16S Metagenomic Sequencing Pipeline

Background

Substantial progress in high-throughput metagenomic sequencing methodologies has enabled the characterisation of bacteria from various origins (for example gut and skin). However, the recently-discovered bacterial microbiota present within animal internal tissues has remained unexplored due to technical difficulties associated with these challenging samples.

Results

We have optimized a specific 16S rDNA-targeted metagenomics sequencing (16S metabarcoding) pipeline based on the Illumina MiSeq technology for the analysis of bacterial DNA in human and animal tissues. This was successfully achieved in various mouse tissues despite the high abundance of eukaryotic DNA and PCR inhibitors in these samples. We extensively tested this pipeline on mock communities, negative controls, positive controls and tissues and demonstrated the presence of novel tissue specific bacterial DNA profiles in a variety of organs (including brain, muscle, adipose tissue, liver and heart).

Conclusion

The high throughput and excellent reproducibility of the method ensured exhaustive and precise coverage of the 16S rDNA bacterial variants present in mouse tissues. This optimized 16S metagenomic sequencing pipeline will allow the scientific community to catalogue the bacterial DNA profiles of different tissues and will provide a database to analyse host/bacterial interactions in relation to homeostasis and disease.     

Intra-Genomic Heterogeneity in 16S rRNA Genes in Strictly Anaerobic Clinical Isolates from Periodontal Abscesses

Background

Members of the genera Prevotella, Veillonella and Fusobacterium are the predominant culturable obligate anaerobic bacteria isolated from periodontal abscesses. When determining the cumulative number of clinical anaerobic isolates from periodontal abscesses, ambiguous or overlapping signals were frequently encountered in 16S rRNA gene sequencing chromatograms, resulting in ambiguous identifications. With the exception of the genus Veillonella, the high intra-chromosomal heterogeneity of rrs genes has not been reported.

Methods

The 16S rRNA genes of 138 clinical, strictly anaerobic isolates and one reference strain were directly sequenced, and the chromatograms were carefully examined. Gene cloning was performed for 22 typical isolates with doublet sequencing signals for the 16S rRNA genes, and four copies of the rrs-ITS genes of 9 Prevotella intermedia isolates were separately amplified by PCR, sequenced and compared. Five conserved housekeeping genes, hsp60, recA, dnaJ, gyrB1 and rpoB from 89 clinical isolates of Prevotella were also amplified by PCR and sequenced for identification and phylogenetic analysis along with 18 Prevotella reference strains.

Results

Heterogeneity of 16S rRNA genes was apparent in clinical, strictly anaerobic oral bacteria, particularly in the genera Prevotella and Veillonella. One hundred out of 138 anaerobic strains (72%) had intragenomic nucleotide polymorphisms (SNPs) in multiple locations, and 13 strains (9.4%) had intragenomic insertions or deletions in the 16S rRNA gene. In the genera Prevotella and Veillonella, 75% (67/89) and 100% (19/19) of the strains had SNPs in the 16S rRNA gene, respectively. Gene cloning and separate amplifications of four copies of the rrs-ITS genes confirmed that 2 to 4 heterogeneous 16S rRNA copies existed.

Conclusion

Sequence alignment of five housekeeping genes revealed that intra-species nucleotide similarities were very high in the genera Prevotella, ranging from 94.3–100%. However, the inter-species similarities were relatively low, ranging from 68.7–97.9%. The housekeeping genes rpoB and gyrB1 were demonstrated to be alternative classification markers to the species level based on intra- and inter-species comparisons, whereas based on phylogenetic tree rpoB proved to be reliable phylogenetic marker for the genus Prevotella.     

Molecular Identification of Staphylococcus aureus in Airway Samples from Children with Cystic Fibrosis

Background

Staphylococcus aureus is a common and significant pathogen in cystic fibrosis. We sought to determine if quantitative PCR (qPCR) and 16S rRNA gene sequencing could provide a rapid, culture-independent approach to the identification of S. aureus airway infections.

Methods

We examined the sensitivity and specificity of two qPCR assays, targeting the femA and 16S rRNA gene, using culture as the gold standard. In addition, 16S rRNA gene sequencing to identify S. aureus directly from airway samples was evaluated. DNA extraction was performed with and without prior enzymatic digestion.

Results

87 samples [42 oropharyngeal (OP) and 45 expectorated sputum (ES)] were analyzed. 59 samples (68%) cultured positive for S. aureus. Using standard extraction techniques, sequencing had the highest sensitivity for S. aureus detection (85%), followed by FemA qPCR (52%) and 16SrRNA qPCR (34%). For all assays, sensitivity was higher from ES samples compared to OP swabs. Specificity of the qPCR assays was 100%, but 21.4% for sequencing due to detection of S. aureus in low relative abundance from culture negative samples. Enzymatic digestion increased the sensitivity of qPCR assays, particularly for OP swabs.

Conclusion

Sequencing had a high sensitivity for S. aureus, but low specificity. While femA qPCR had higher sensitivity than 16S qPCR for detection of S. aureus, neither assay was as sensitive as sequencing. The significance of S. aureus detection with low relative abundance by sequencing in culture-negative specimens is not clear.     

Urease and Dental Plaque Microbial Profiles in Children

Objective

Urease enzymes produced by oral bacteria generate ammonia, which can have a significant impact on the oral ecology and, consequently, on oral health. To evaluate the relationship of urease with dental plaque microbial profiles in children as it relates to dental caries, and to identify the main contributors to this activity.

Methods

82 supragingival plaque samples were collected from 44 children at baseline and one year later, as part of a longitudinal study on urease and caries in children. DNA was extracted; the V3–V5 region of the 16S rRNA gene was amplified and sequenced using 454 pyrosequencing. Urease activity was measured using a spectrophotometric assay. Data were analyzed with Qiime.

Results

Plaque urease activity was significantly associated with the composition of the microbial communities of the dental plaque (Baseline P = 0.027, One Year P = 0.012). The bacterial taxa whose proportion in dental plaque exhibited significant variation by plaque urease levels in both visits were the family Pasteurellaceae (Baseline P<0.001; One Year P = 0.0148), especially Haemophilus parainfluenzae. No association was observed between these bacteria and dental caries. Bacteria in the genus Leptotrichia were negatively associated with urease and positively associated with dental caries (Bonferroni P<0.001).

Conclusions

Alkali production by urease enzymes primarily from species in the family Pasteurellaceae can be an important ecological determinant in children’s dental plaque. Further studies are needed to establish the role of urease-associated bacteria in the acid/base homeostasis of the dental plaque, and in the development and prediction of dental caries in children.     

Comparative Analysis of Pyrosequencing and a Phylogenetic Microarray for Exploring Microbial Community Structures in the Human Distal Intestine

Background

Variations in the composition of the human intestinal microbiota are linked to diverse health conditions. High-throughput molecular technologies have recently elucidated microbial community structure at much higher resolution than was previously possible. Here we compare two such methods, pyrosequencing and a phylogenetic array, and evaluate classifications based on two variable 16S rRNA gene regions.

Methods and Findings

Over 1.75 million amplicon sequences were generated from the V4 and V6 regions of 16S rRNA genes in bacterial DNA extracted from four fecal samples of elderly individuals. The phylotype richness, for individual samples, was 1,400–1,800 for V4 reads and 12,500 for V6 reads, and 5,200 unique phylotypes when combining V4 reads from all samples. The RDP-classifier was more efficient for the V4 than for the far less conserved and shorter V6 region, but differences in community structure also affected efficiency. Even when analyzing only 20% of the reads, the majority of the microbial diversity was captured in two samples tested. DNA from the four samples was hybridized against the Human Intestinal Tract (HIT) Chip, a phylogenetic microarray for community profiling. Comparison of clustering of genus counts from pyrosequencing and HITChip data revealed highly similar profiles. Furthermore, correlations of sequence abundance and hybridization signal intensities were very high for lower-order ranks, but lower at family-level, which was probably due to ambiguous taxonomic groupings.

Conclusions

The RDP-classifier consistently assigned most V4 sequences from human intestinal samples down to genus-level with good accuracy and speed. This is the deepest sequencing of single gastrointestinal samples reported to date, but microbial richness levels have still not leveled out. A majority of these diversities can also be captured with five times lower sampling-depth. HITChip hybridizations and resulting community profiles correlate well with pyrosequencing-based compositions, especially for lower-order ranks, indicating high robustness of both approaches. However, incompatible grouping schemes make exact comparison difficult.     

Variable copy number, intra-genomic heterogeneities and lateral transfers of the 16S rRNA gene in Pseudomonas

Even though the 16S rRNA gene is the most commonly used taxonomic marker in microbial ecology, its poor resolution is still not fully understood at the intra-genus level. In this work, the number of rRNA gene operons, intra-genomic heterogeneities and lateral transfers were investigated at a fine-scale resolution, throughout the Pseudomonas genus. In addition to nineteen sequenced Pseudomonas strains, we determined the 16S rRNA copy number in four other Pseudomonas strains by Southern hybridization and Pulsed-Field Gel Electrophoresis, and studied the intra-genomic heterogeneities by Denaturing Gradient Gel Electrophoresis and sequencing. Although the variable copy number (from four to seven) seems to be correlated with the evolutionary distance, some close strains in the P. fluorescens lineage showed a different number of 16S rRNA genes, whereas all the strains in the P. aeruginosa lineage displayed the same number of genes (four copies). Further study of the intra-genomic heterogeneities revealed that most of the Pseudomonas strains (15 out of 19 strains) had at least two different 16S rRNA alleles. A great difference (5 or 19 nucleotides, essentially grouped near the V1 hypervariable region) was observed only in two sequenced strains. In one of our strains studied (MFY30 strain), we found a difference of 12 nucleotides (grouped in the V3 hypervariable region) between copies of the 16S rRNA gene. Finally, occurrence of partial lateral transfers of the 16S rRNA gene was further investigated in 1803 full-length sequences of Pseudomonas available in the databases. Remarkably, we found that the two most variable regions (the V1 and V3 hypervariable regions) had probably been laterally transferred from another evolutionary distant Pseudomonas strain for at least 48.3 and 41.6% of the 16S rRNA sequences, respectively. In conclusion, we strongly recommend removing these regions of the 16S rRNA gene during the intra-genus diversity studies.     

Helicobacter pylori Eradication Causes Perturbation of the Human Gut Microbiome in Young Adults

Background

Accumulating evidence shows that Helicobacter pylori protects against some metabolic and immunological diseases in which the development of these diseases coincide with temporal or permanent dysbiosis. The aim of this study was to assess the effect of H. pylori eradication on the human gut microbiome.

Methods

As part of the currently on-going ESSAY (Eradication Study in Stable Adults/Youths) study, we collected stool samples from 17 H. pylori-positive young adult (18–30 years-old) volunteers. The same cohort was followed up 6, 12 and 18 months-post H. pylori eradication. The impact of H. pylori on the human gut microbiome pre- and post-eradication was investigated using high throughput 16S rRNA gene (V3-V4 region) sequencing using the Illumina Miseq followed by data analysis using Qiime pipeline.

Results

We compared the composition and diversity of bacterial communities in the fecal microbiome of the H. pylori-positive volunteers, before and after H. pylori eradication therapy. The 16S rRNA gene was sequenced at an average of 150,000–170,000 reads/sample. The microbial diversity were similar pre- and post-H. pylori eradication with no significant differences in richness and evenness of bacterial species. Despite that the general profile of the gut microbiome was similar pre- and post-eradication, some changes in the bacterial communities at the phylum and genus levels were notable, particularly the decrease in relative abundance of Bacterioidetes and corresponding increase in Firmicutes after H. pylori eradication. The significant increase of short-chain fatty acids (SCFA)-producing bacteria genera could also be associated with increased risk of metabolic disorders.

Conclusions

Our preliminary stool metagenomics study shows that eradication of H. pylori caused perturbation of the gut microbiome and may indirectly affect the health of human. Clinicians should be aware of the effect of broad spectrum antibiotics used in H. pylori eradication regimen and be cautious in the clinical management of H. pylori infection, particularly in immunocompromised patients.     

Diagnosis of Bacterial Bloodstream Infections: A 16S Metagenomics Approach

Background

Bacterial bloodstream infection (bBSI) is one of the leading causes of death in critically ill patients and accurate diagnosis is therefore crucial. We here report a 16S metagenomics approach for diagnosing and understanding bBSI.

Methodology/Principal Findings

The proof-of-concept was delivered in 75 children (median age 15 months) with severe febrile illness in Burkina Faso. Standard blood culture and malaria testing were conducted at the time of hospital admission. 16S metagenomics testing was done retrospectively and in duplicate on the blood of all patients. Total DNA was extracted from the blood and the V3–V4 regions of the bacterial 16S rRNA genes were amplified by PCR and deep sequenced on an Illumina MiSeq sequencer. Paired reads were curated, taxonomically labeled, and filtered. Blood culture diagnosed bBSI in 12 patients, but this number increased to 22 patients when combining blood culture and 16S metagenomics results. In addition to superior sensitivity compared to standard blood culture, 16S metagenomics revealed important novel insights into the nature of bBSI. Patients with acute malaria or recovering from malaria had a 7-fold higher risk of presenting polymicrobial bloodstream infections compared to patients with no recent malaria diagnosis (p-value = 0.046). Malaria is known to affect epithelial gut function and may thus facilitate bacterial translocation from the intestinal lumen to the blood. Importantly, patients with such polymicrobial blood infections showed a 9-fold higher risk factor for not surviving their febrile illness (p-value = 0.030).

Conclusions/Significance

Our data demonstrate that 16S metagenomics is a powerful approach for the diagnosis and understanding of bBSI. This proof-of-concept study also showed that appropriate control samples are crucial to detect background signals due to environmental contamination.     

16S rRNA gene sequencing on a benchtop sequencer: accuracy for identification of clinically important bacteria

Aims

Test the choice of 16S rRNA gene amplicon and data analysis method on the accuracy of identification of clinically important bacteria utilizing a benchtop sequencer.

Methods and Results

Nine 16S rRNA amplicons were tested on an Ion Torrent PGM to identify 41 strains of clinical importance. The V1–V2 region identified 40 of 41 isolates to the species level. Three data analysis methods were tested, finding that the Ribosomal Database Project's SequenceMatch outperformed BLAST and the Ion Reporter Metagenomics analysis pipeline. Lastly, 16S rRNA gene sequencing mixtures of four species through a six log range of dilution showed species were identifiable even when present as 0·1% of the mixture.

Conclusions

Sequencing the V1–V2 16S rRNA gene region, made possible by the increased read length Ion Torrent PGM sequencer's 400 base pair chemistry, may be a better choice over other commonly used regions for identifying clinically important bacteria. In addition, the SequenceMatch algorithm, freely available from the Ribosomal Database Project, is a good choice for matching filtered reads to organisms. Lastly, 16S rRNA gene sequencing's sensitivity to the presence of a bacterial species at 0·1% of a mixture suggests it has sufficient sensitivity for samples in which important bacteria may be rare.

Significance and Impact of the Study

We have validated 16S rRNA gene sequencing on a benchtop sequencer including simple mixtures of organisms; however, our results highlight deficits for clinical application in place of current identification methods.     

IM-TORNADO: A Tool for Comparison of 16S Reads from Paired-End Libraries

Motivation

16S rDNA hypervariable tag sequencing has become the de facto method for accessing microbial diversity. Illumina paired-end sequencing, which produces two separate reads for each DNA fragment, has become the platform of choice for this application. However, when the two reads do not overlap, existing computational pipelines analyze data from read separately and underutilize the information contained in the paired-end reads.

Results

We created a workflow known as Illinois Mayo Taxon Organization from RNA Dataset Operations (IM-TORNADO) for processing non-overlapping reads while retaining maximal information content. Using synthetic mock datasets, we show that the use of both reads produced answers with greater correlation to those from full length 16S rDNA when looking at taxonomy, phylogeny, and beta-diversity.

Availability and Implementation

IM-TORNADO is freely available at http://sourceforge.net/projects/imtornado and produces BIOM format output for cross compatibility with other pipelines such as QIIME, mothur, and phyloseq.     

Defining reference sequences for Nocardia species by similarity and clustering analyses of 16S rRNA gene sequence data

Background

The intra- and inter-species genetic diversity of bacteria and the absence of ‘reference’, or the most representative, sequences of individual species present a significant challenge for sequence-based identification. The aims of this study were to determine the utility, and compare the performance of several clustering and classification algorithms to identify the species of 364 sequences of 16S rRNA gene with a defined species in GenBank, and 110 sequences of 16S rRNA gene with no defined species, all within the genus Nocardia.

Methods

A total of 364 16S rRNA gene sequences of Nocardia species were studied. In addition, 110 16S rRNA gene sequences assigned only to the Nocardia genus level at the time of submission to GenBank were used for machine learning classification experiments. Different clustering algorithms were compared with a novel algorithm or the linear mapping (LM) of the distance matrix. Principal Components Analysis was used for the dimensionality reduction and visualization.

Results

The LM algorithm achieved the highest performance and classified the set of 364 16S rRNA sequences into 80 clusters, the majority of which (83.52%) corresponded with the original species. The most representative 16S rRNA sequences for individual Nocardia species have been identified as ‘centroids’ in respective clusters from which the distances to all other sequences were minimized; 110 16S rRNA gene sequences with identifications recorded only at the genus level were classified using machine learning methods. Simple kNN machine learning demonstrated the highest performance and classified Nocardia species sequences with an accuracy of 92.7% and a mean frequency of 0.578.

Conclusion

The identification of centroids of 16S rRNA gene sequence clusters using novel distance matrix clustering enables the identification of the most representative sequences for each individual species of Nocardia and allows the quantitation of inter- and intra-species variability.     

Analysis of the Cystic Fibrosis Lung Microbiota via Serial Illumina Sequencing of Bacterial 16S rRNA Hypervariable Regions

The characterization of bacterial communities using DNA sequencing has revolutionized our ability to study microbes in nature and discover the ways in which microbial communities affect ecosystem functioning and human health. Here we describe Serial Illumina Sequencing (SI-Seq): a method for deep sequencing of the bacterial 16S rRNA gene using next-generation sequencing technology. SI-Seq serially sequences portions of the V5, V6 and V7 hypervariable regions from barcoded 16S rRNA amplicons using an Illumina short-read genome analyzer. SI-Seq obtains taxonomic resolution similar to 454 pyrosequencing for a fraction of the cost, and can produce hundreds of thousands of reads per sample even with very high multiplexing. We validated SI-Seq using single species and mock community controls, and via a comparison to cystic fibrosis lung microbiota sequenced using 454 FLX Titanium. Our control runs show that SI-Seq has a dynamic range of at least five orders of magnitude, can classify >96% of sequences to the genus level, and performs just as well as 454 and paired-end Illumina methods in estimation of standard microbial ecology diversity measurements. We illustrate the utility of SI-Seq in a pilot sample of central airway secretion samples from cystic fibrosis patients.     

Assessment of the Microbial Constituents of the Home Environment of Individuals with Cystic Fibrosis (CF) and Their Association with Lower Airways Infections

Introduction

Cystic fibrosis (CF) airways are colonized by a polymicrobial community of organisms, termed the CF microbiota. We sought to define the microbial constituents of the home environment of individuals with CF and determine if it may serve as a latent reservoir for infection.

Methods

Six patients with newly identified CF pathogens were included. An investigator collected repeat sputum and multiple environmental samples from their homes. Bacteria were cultured under both aerobic and anaerobic conditions. Morphologically distinct colonies were selected, purified and identified to the genus and species level through 16S rRNA gene sequencing. When concordant organisms were identified in sputum and environment, pulsed-field gel electrophoresis (PFGE) was performed to determine relatedness. Culture-independent bacterial profiling of each sample was carried out by Illumina sequencing of the V3 region of the 16s RNA gene.

Results

New respiratory pathogens prompting investigation included: Mycobacterium abscessus(2), Stenotrophomonas maltophilia(3), Pseudomonas aeruginosa(3), Pseudomonas fluorescens(1), Nocardia spp.(1), and Achromobacter xylosoxidans(1). A median 25 organisms/patient were cultured from sputum. A median 125 organisms/home were cultured from environmental sites. Several organisms commonly found in the CF lung microbiome were identified within the home environments of these patients. Concordant species included members of the following genera: Brevibacterium(1), Microbacterium(1), Staphylococcus(3), Stenotrophomonas(2), Streptococcus(2), Sphingomonas(1), and Pseudomonas(4). PFGE confirmed related strains (one episode each of Sphinogomonas and P. aeruginosa) from the environment and airways were identified in two patients. Culture-independent assessment confirmed that many organisms were not identified using culture-dependent techniques.

Conclusions

Members of the CF microbiota can be found as constituents of the home environment in individuals with CF. While the majority of isolates from the home environment were not genetically related to those isolated from the lower airways of individuals with CF suggesting alternate sources of infection were more common, a few genetically related isolates were indeed identified. As such, the home environment may rarely serve as either the source of infection or a persistent reservoir for re-infection after clearance.     

The Oral Microbiome of Denture Wearers Is Influenced by Levels of Natural Dentition

Objectives

The composition of dental plaque has been well defined, whereas currently there is limited understanding of the composition of denture plaque and how it directly influences denture related stomatitis (DS). The aims of this study were to compare the microbiomes of denture wearers, and to understand the implications of these towards inter-kingdom and host-pathogen interactions within the oral cavity.

Methods

Swab samples were obtained from 123 participants wearing either a complete or partial denture; the bacterial composition of each sample was determined using bar-coded illumina MiSeq sequencing of the bacterial hypervariable V4 region of 16S rDNA. Sequencing data processing was undertaken using QIIME, clustered in Operational Taxonomic Units (OTUs) and assigned to taxonomy. The dentures were sonicated to remove the microbial flora residing on the prosthesis, sonicate was then cultured using diagnostic colorex Candida media. Samples of unstimulated saliva were obtained and antimicrobial peptides (AMP) levels were measured by ELISA.

Results

We have shown that dental and denture plaques are significantly distinct both in composition and diversity and that the oral microbiome composition of a denture wearer is variable and is influenced by the location within the mouth. Dentures and mucosa were predominantly made up of Bacilli and Actinobacteria. Moreover, the presence of natural teeth has a significant impact on the overall microbial composition, when compared to the fully edentulous. Furthermore, increasing levels of Candida spp. positively correlate with Lactobacillus spp. AMPs were quantified, though showed no specific correlations.

Conclusions

This is the first study to provide a detailed understanding of the oral microbiome of denture wearers and has provided evidence that DS development is more complex than simply a candidal infection. Both fungal and bacterial kingdoms clearly play a role in defining the progression of DS, though we were unable to show a defined role for AMPs.     

NoDe: a fast error-correction algorithm for pyrosequencing amplicon reads

Background

The popularity of new sequencing technologies has led to an explosion of possible applications, including new approaches in biodiversity studies. However each of these sequencing technologies suffers from sequencing errors originating from different factors. For 16S rRNA metagenomics studies, the 454 pyrosequencing technology is one of the most frequently used platforms, but sequencing errors still lead to important data analysis issues (e.g. in clustering in taxonomic units and biodiversity estimation). Moreover, retaining a higher portion of the sequencing data by preserving as much of the read length as possible while maintaining the error rate within an acceptable range, will have important consequences at the level of taxonomic precision.

Results

The new error correction algorithm proposed in this work - NoDe (Noise Detector) - is trained to identify those positions in 454 sequencing reads that are likely to have an error, and subsequently clusters those error-prone reads with correct reads resulting in error-free representative read. A benchmarking study with other denoising algorithms shows that NoDe can detect up to 75% more errors in a large scale mock community dataset, and this with a low computational cost compared to the second best algorithm considered in this study. The positive effect of NoDe in 16S rRNA studies was confirmed by the beneficial effect on the precision of the clustering of pyrosequencing reads in operational taxonomic units.

Conclusions

NoDe was shown to be a computational efficient denoising algorithm for pyrosequencing reads, producing the lowest error rates in an extensive benchmarking study with other denoising algorithms.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0520-5) contains supplementary material, which is available to authorized users.     

Identification of Bacteriology and Risk Factor Analysis of Asymptomatic Bacterial Colonization in Pacemaker Replacement Patients

Background

Recent researches revealed that asymptomatic bacterial colonization on PMs might be ubiquitous and increase the risk of clinical PM infection. Early diagnosis of patients with asymptomatic bacterial colonization could provide opportunity for targeted preventive measures.

Objective

The present study explores the incidence of bacterial colonization of generator pockets in pacemaker replacement patients without signs of infection, and to analyze risk factors for asymptomatic bacterial colonization.

Methods

From June 2011 to December 2013, 118 patients underwent pacemaker replacement or upgrade. Identification of bacteria was carried out by bacterial culture and 16S rRNA sequencing. Clinical risk characteristics were analyzed.

Results

The total bacterial positive rate was 37.3% (44 cases), and the coagulase-negative Staphylococcus aureus detection rate was the highest. Twenty two (18.6%) patients had positive bacterial culture results, of which 50% had coagulase-negative staphylococcus. The bacterial DNA detection rate was 36.4 % (43 cases). Positive bacterial DNA results from pocket tissues and the surface of the devices were 22.0% and 29.7%, respectively. During follow-up (median, 27.0 months), three patients (6.8%, 3/44) became symptomatic with the same genus of microorganism, S. aureus (n=2) and S. epidermidis (n=1). Multivariable logistic regression analysis showed that history of bacterial infection, use of antibiotics, application of antiplatelet drugs, replacement frequency were independent risk factors for asymptomatic bacterial colonization.

Conclusion

There was a high incidence of asymptomatic bacterial colonization in pacemaker patients with independent risk factors. Bacterial culture combined genetic testing could improve the detection rate.