Instruction of microbiome taxonomic profiling based on 16S rRNA sequencing

Journal of Microbiology - Recent studies on microbiome highlighted their importance in various environments including human, where they are involved in multiple biological contexts such as immune...     

MSClust: A Multi-Seeds based Clustering algorithm for microbiome profiling using 16S rRNA sequence

Recent developments of next generation sequencing technologies have led to rapid accumulation of 16S rRNA sequences for microbiome profiling. One key step in data processing is to cluster short sequences into operational taxonomic units (OTUs). Although many methods have been proposed for OTU inferences, a major challenge is the balance between inference accuracy and computational efficiency, where inference accuracy is often sacrificed to accommodate the need to analyze large numbers of sequences. Inspired by the hierarchical clustering method and a modified greedy network clustering algorithm, we propose a novel multi-seeds based heuristic clustering method, named MSClust, for OTU inference. MSClust first adaptively selects multi-seeds instead of one seed for each candidate cluster, and the reads are then processed using a greedy clustering strategy. Through many numerical examples, we demonstrate that MSClust enjoys less memory usage, and better biological accuracy compared to existing heuristic clustering methods while preserving efficiency and scalability.     

Evaluation of the MicroSeq 500 16S rDNA-based gene sequencing for the diagnosis of culture-negative bacterial meningitis

Gene amplification using 16S rDNA primers has been proposed as a strategy for the diagnosis of bacterial meningitis. The aim of this study was to evaluate the performance of the MicroSeq 500 16S ribosomal DNA test (Applied Biosystems) from patients with suspected bacterial meningitis and CSF negative-culture in comparison to traditional methods. Twelve purulent culture-negative CSF samples were collected between January 2005 and January 2007. For DNA extraction, 500 microl of CSF samples were treated using the QIAamp mini kit (QIAGEN). The extracted DNA was examined amplifying 500 bp at the 5' end of 16S rRNA gene using MicroSeq500 16S rDNA Bacterial Identification PCR kit and the sequencing reactions were performed with the MicroSeq500 16S rDNA Bacterial Identification Sequencing kit (Applied Biosystems). The sequences were compared with those available in GenBank. For the culture-negative CSF samples the MicroSeq 500 16S rDNA yielded a positive result in 9 cases (75.0%): three samples were identified as Streptococcus. pneumoniae, three as Neisseria meningitidis, and the remaining 3 as Haemophilus influenzae, Abiotrophia defectiva and Porphyromonas gingivalis. The MicroSeq 500 16S ribosomal DNA test may improve the microbiological diagnosis of bacterial meningitis, especially when spinal fluid samples are obtained after the administration of antimicrobial therapy.     

Variations in 16S rRNA-based microbiome profiling between pyrosequencing runs and between pyrosequencing facilities

Pyrosequencing of 16S rRNA gene amplicons on the 454 FLX Titanium platform has been widely used to analyze microbiomes in various environments. However, different results may stem from variations among sequencing runs or among sequencing facilities. This study aimed to evaluate these variations between different pyrosequencing runs by sequencing 16S rRNA gene amplicon libraries generated from three sets of rumen samples twice each on the 454 FLX Titanium system at two independent sequencing facilities. Similar relative abundances were found for predominant taxa represented by large numbers of sequence reads but not for minor taxa represented by small numbers of sequence reads. The two sequencing facilities revealed different bacterial profiles with respect to both predominant taxa and minor taxa, including the most predominant genus Prevotella, the family Lachnospiraceae, and the phylum Proteobacteria. Differences in primers used to generate amplicon libraries may be a major source of variations in microbiome profiling. Because different primers and regions of 16S rRNA genes are often used by different researchers, significant variations likely exist among studies. Quantitative interpretation for relative abundance of taxa, especially minor taxa, from prevalence of sequence reads and comparisons of results from different studies should be done with caution.     

Influence of 16S rRNA reference databases in amplicon-based environmental microbiome research

Biotechnology Letters - The reference databases play a pivotal role in amplicon microbiome research, however these databases differ in the sequence content and taxonomic information available....     

16S rRNA gene profiling of rhizospheric microbial community of Eichhornia crassipes

Molecular Biology Reports - The rhizosphere of a plant is an important interface for the plant-microbe interaction that plays a significant role in the uptake and removal of heavy metal from...     

Characterization of alkalotolerant bacterial community by 16S rDNA-based denaturing gradient gel electrophoresis method for degradation of dibenzofuran in soil microcosm

An alkalotolerant bacterial community was developed by continuous enrichment in the chemostat in presence of dibenzofuran (DF) as sole carbon source. Six different types of bacterial isolates were cultured on nutrient broth agar plates together with six operational taxonomic units (OTUs) at pH 7.0 and pH 8.0 by 16S rDNA-DGGE method. However, isolates of microbial community was declined from three OTUs (pH 9.0) to two at pH 10.0 after enrichment in alkaline condition. Among the six isolates tested for degradation of DF, Pseudomonas sp. and Bacillus sp. the members of alkalotolerant bacterial community had better potency to degrade dibenzofuran. Alkalotolerant bacterial community introduced in soil microcosm for evaluation of survival of most suitable isolates and degradation of dioxin-like compound indicated more than 90% degradation of dibenzofuran after 45 days by the bacterial community enriched for 180 days in the chemostat at pH 10, however, microbial community was not competent to utilize even 50% DF after day 30, not enriched in the chemostat. The survival of competent bacteria monitored by DGGE method in soil microcosm indicated presence of two major alkalotolerant isolates for utilization of dibenzofuran, substantiated the results and significance of alkalotolerant bacteria for in situ bioremediation of dioxin-like compounds in the environment.     

Development of 16S rDNA-based PCR assay for detecting centipeda periodontii and selenomonas sputigena

To detect oral motile bacteria directly from dental plaque, specific PCR primers for Centipeda periodontii and Selenomonas sputigena were designed based on the sequence analysis of their 16S rDNA. The primers were specific and sensitive enough to amplify DNA fragments from the available oral bacteria. The detection limit was fewer than 10 bacterial cells per sample. It was also possible to detect these bacteria in dental plaque. The prevalence of these bacteria varied in each sample. The specific primers designed in this study may clarify the epidemiology of periodontal disease.     

Assessment of bacterial community structure in the deep sub-seafloor biosphere by 16S rDNA-based techniques: a cautionary tale

Investigations into the deep marine environment have demonstrated the presence of a significant microbial biomass buried deep within sediments on a global scale. It is now believed that this deep biosphere plays a major role in the global cycling of elements and contains a large reservoir of organic carbon. This paper reports the development of a DNA extraction protocol that addresses the particular problems faced in applying molecular ecological techniques to samples containing very low biomass. Sediment samples were collected from different geographical locations within the Pacific Ocean and include the Ocean Drilling Program (ODP) Leg 190, Nankai Trough Accretionary Prism. Seven DNA extraction protocols were tested and a commercially available DNA extraction kit with modifications was shown to produce higher yields of polymerase chain reaction (PCR)-amplifiable DNA than standard laboratory methods. Denaturing gradient gel electrophoresis (DGGE) analysis of 16S rRNA gene diversity revealed that template DNA from these extremely low biomass sediment samples was susceptible to PCR bias and random amplification. We propose that it is essential to screen 16S rRNA gene products for bacterial diversity by DGGE or other rapid fingerprinting methods, prior to their use in establishing a representative clone library of deep sub-seafloor bacteria. This represents a cautionary approach to analysis of microbial diversity in such sub-seafloor ecosystems.     

Characterization of depth-related population variation in microbial communities of a coastal marine sediment using 16S rDNA-based approaches and quinone profiling

Depth-related changes in whole-community structure were evaluated in a coastal marine sediment using a molecular fingerprinting method, terminal restriction fragment length polymorphism (T-RFLP) analysis, and a chemotaxonomic technique (quinone profiling). Dendrograms derived from both T-RFLP analysis and quinone profiling indicated a significant variation in microbial community structure between the 0-2 cm layer and deeper layers. This corresponded to the dramatic change in the redox potential, acid-volatile sulphide-sulphur and bacterial numbers observed at 0-2 cm and 2-4 cm depths. A significant change in the number of terminal restriction fragments (T-RFs) was also detected at this transition depth. However, the change in major T-RFs with depth was not seen in electropherograms. The population changes were primarily variations in minor ribotypes. Most quinone homologues were detected at all depths, although the quinone composition changed with depth. Therefore, quinone profiling also suggested that the depth-related variation was primarily attributable to minor bacterial groups rather than change in the major population structure. 16S rDNA clone library analysis revealed that clones belonging to the genera Vibrio and Serratia predominated as major bacterial groups at all depths. Our data suggested that the sediment community might result from sedimentation effects of sinking particles. Overall, our results demonstrated that the combined methods of T-RFLP analysis and quinone profiling were effective for assessing depth-related microbial populations.     

A core human microbiome as viewed through 16S rRNA sequence clusters

We explore the microbiota of 18 body sites in over 200 individuals using sequences amplified V1-V3 and the V3-V5 small subunit ribosomal RNA (16S) hypervariable regions as part of the NIH Common Fund Human Microbiome Project. The body sites with the greatest number of core OTUs, defined as OTUs shared amongst 95% or more of the individuals, were the oral sites (saliva, tongue, cheek, gums, and throat) followed by the nose, stool, and skin, while the vaginal sites had the fewest number of OTUs shared across subjects. We found that commonalities between samples based on taxonomy could sometimes belie variability at the sub-genus OTU level. This was particularly apparent in the mouth where a given genus can be present in many different oral sites, but the sub-genus OTUs show very distinct site selection, and in the vaginal sites, which are consistently dominated by the Lactobacillus genus but have distinctly different sub-genus V1-V3 OTU populations across subjects. Different body sites show approximately a ten-fold difference in estimated microbial richness, with stool samples having the highest estimated richness, followed by the mouth, throat and gums, then by the skin, nasal and vaginal sites. Richness as measured by the V1-V3 primers was consistently higher than richness measured by V3-V5. We also show that when such a large cohort is analyzed at the genus level, most subjects fit the stool "enterotype" profile, but other subjects are intermediate, blurring the distinction between the enterotypes. When analyzed at the finer-scale, OTU level, there was little or no segregation into stool enterotypes, but in the vagina distinct biotypes were apparent. Finally, we note that even OTUs present in nearly every subject, or that dominate in some samples, showed orders of magnitude variation in relative abundance emphasizing the highly variable nature across individuals.     

Improving the specificity of 16S rDNA-based polymerase chain reaction for detecting Borrelia burgdorferi sensu lato-causative agents of human Lyme disease

AIMS: 16S rDNA sequences of Borrelia burgdorferi sensu lato were aligned with the 16S rDNA sequences of Borrelia hermsii, Borrelia turicatae, and Borrelia lonestari in order to identify primers that might be used to more specifically identify agents of human Lyme disease in ticks in human skin samples. METHODS AND RESULTS: Standard polymerase chain reaction (PCR), using an oligonucleotide sequence, designated TEC1, was shown, in combination with a previously developed primer (LD2) to amplify strains of B. burgdorferi sensu stricto, Borrelia afzelii, and Borrelia garinii, but not the non-Lyme causing B. hermsii or B. turicatae. This primer pair, designated Bbsl, was successfully used to amplify B. burgdorferi sensu lato from skin biopsies of patients with Lyme disease symptoms as well as from Ixodes scapularis, Amblyomma americanum and Dermacentor variabilis ticks. CONCLUSIONS: The primer set Bbsl allows for the rapid detection and differentiation of B. burgdorferi sensu lato from non-Lyme disease-causing Borrelia species in ticks and human tissues. SIGNIFICANCE AND IMPACT OF THE STUDY: The PCR primer set, Bbsl, will greatly facilitate detection of the causative agents of Lyme disease in infected ticks and human skin samples assisting in epidemiological studies, and potentially allowing for a more rapid diagnosis of the disease in patients.     

Multiple displacement amplification of DNA from human colon and rectum biopsies: bacterial profiling and identification of Helicobacter pylori-DNA by means of 16S rDNA-based TTGE and pyrosequencing analysis

Amplifying bacterial DNA by PCR from human biopsy specimens has sometimes proved to be difficult, mainly due to the low amount of bacterial DNA present. Therefore, nested or semi-nested 16S rDNA PCR amplification has been the method of choice. In this study, we evaluate the potential use of whole genome amplification of total DNA isolated from human colon and rectum biopsy specimens, followed by 16S rDNA PCR amplification of multiple displacement amplified (MDA)-DNA. Subsequently, a H. pylori-specific 16S rDNA variable V3 region PCR assay was applied directly on MDA-DNA and, combined with pyrosequencing analysis; the presence of H. pylori in some biopsies from colon in patients with microscopic colitis was confirmed. Furthermore, temporal temperature gradient gel electrophoresis (TTGE) of 16S rDNA amplicons using primers flanking variable regions V3, V4, and V9, was used to establish bacterial profiles from individual biopsies. A variation of the bacterial profiles in the colonic mucosa in microscopic colitis and in normal rectal mucosa was observed. In conclusion we find the MDA technique to be a useful method to overcome the problem of insufficient bacterial DNA in human biopsy specimens.     

16S rDNA-based phylogeny of non-symbiotic bacteria of Entorno-pathogenic nematodes from infected insect cadavers

Using 16S rDNA gene sequencing technique, three different species of non-symbiotic bacteria of entomopatho-genic nematodes (EPNs) (Steinernema sp.and Heterorhabditis sp.) were isolated and identified from infected insect cadavers(Galleria mellonella larvae) after 48-hour post infections.Sequence similarity analysis revealed that the strains SRK3, SRK4 and SRK5 belong to Ochrobactrum cytisi,Schineria larvae and Ochrobactrum anthropi,respectively.The isolates O.anthropi and S.larvae were found to be associated with Heterorhabditis indica strains BDU-17 and Yer-136,respectively,whereas O.cytisi was associated with Steinernema siamkayai strain BDU-87. Phenotypically, temporal EPN bacteria were fairly related to symbiotic EPN bacteria (Photorhabdus and Xenorhabdus genera). The strains SRK3 and SRK5 were phylogeographically similar to several non-symbionts and contaminated EPN bacteria isolated in Germany(LMG3311T) and China (X-14),while the strain SRK4 was identical to the isolates of S.larvae (L1/57,L1/58, L1/68 and L2/11) from Wohlfahrtia magnifica in Hungary.The result was further confirmed by RNA secondary structure and minimum energy calculations of aligned sequences.This study suggested that the non-symbionts of these nematodes are phylogeographically diverged in some extent due to phase variation.Therefore,these strains are not host-dependent, but environment-specific isolates.     

16S rDNA-based characterization of BTX-catabolizing microbial associations isolated from a South African sandy soil

In the presence of different selection pressures, particularly pHand electron donor concentration, indigenous microbial associations which catabolize selected petroleum hydrocarbon components (benzene, toluene and o-, m- and p-xylene (BTX)) were enriched and isolated from a petroleum hydrocarbon-contaminated KwaZulu-Natal sandy soil. Electron microscopy revealed that, numerically, rods constituted the majority of the populations responsible for BTX catabolism. Molecular techniques (polymerase chain reaction (PCR) and 16S rDNA fingerprinting by denaturing-gradient gel electrophoresis (DGGE)) were employed to explore the diversities and analyze the structures of the isolated microbial associations. Pearson product-moment correlation indicated that the different, but chemically similar, petroleum hydrocarbon molecules, effectedthe isolation of different associations. However, some similar numerically-dominant bands characterized the associations. A 30% similarity was evident between the m- and o-xylene-catabolizing associations regardless of the molecule concentration and the enrichment pH. PCR-DGGE was also used to complement conventional culture-based microbiological procedures for environmental parameter optimization. Band pattern differencesindicated profile variations of the isolated associations which possibly accounted for the growth rate changes recorded in response to pH and temperature perturbations.