VITCOMIC: visualization tool for taxonomic compositions of microbial communities based on 16S rRNA gene sequences

Metagenomic Thermometer

Various microorganisms exist in environments, and each of them has its optimal growth temperature (OGT). The relationship between genomic information and OGT of each species has long been studied, and one such study revealed that OGT of prokaryotes can be accurately predicted based on the fraction of seven amino acids (IVYWREL) among all encoded amino-acid sequences in its genome. Extending this discovery, we developed a 'Metagenomic Thermometer' as a means of predicting environmental temperature based on metagenomic sequences. Temperature prediction of diverse environments using publicly available metagenomic data revealed that the Metagenomic Thermometer can predict environmental temperatures with small temperature changes and little influx of microorganisms from other environments. The accuracy of the Metagenomic Thermometer was also confirmed by a demonstration experiment using an artificial hot water canal. The Metagenomic Thermometer was also applied to human gut metagenomic samples, yielding a reasonably accurate value for human body temperature. The result further suggests that deep body temperature determines the dominant lineage of the gut community. Metagenomic Thermometer provides a new insight into temperature-driven community assembly based on amino-acid composition rather than microbial taxa.

CDEMI: characterizing differences in microbial composition and function in microbiome data

Microbial communities influence host phenotypes through microbiota-derived metabolites and interactions between exogenous active substances (EASs) and the microbiota. Owing to the high dynamics of microbial community composition and difficulty in microbial functional analysis, the identification of mechanistic links between individual microbes and host phenotypes is complex. Thus, it is important to characterize variations in microbial composition across various conditions (for example, topographical locations, times, physiological and pathological conditions, and populations of different ethnicities) in microbiome studies. However, no web server is currently available to facilitate such characterization. Moreover, accurately annotating the functions of microbes and investigating the possible factors that shape microbial function are critical for discovering links between microbes and host phenotypes. Herein, an online tool, CDEMI, is introduced to discover microbial composition variations across different conditions, and five types of microbe libraries are provided to comprehensively characterize the functionality of microbes from different perspectives. These collective microbe libraries include (1) microbial functional pathways, (2) disease associations with microbes, (3) EASs associations with microbes, (4) bioactive microbial metabolites, and (5) human body habitats. In summary, CDEMI is unique in that it can reveal microbial patterns in distributions/compositions across different conditions and facilitate biological interpretations based on diverse microbe libraries. CDEMI is accessible at http://rdblab.cn/cdemi/.

High-Throughput 16 S rRNA Gene Sequencing Reveals Bacterial Diversity of Infant Formula Production Line Samples in Spring and Summer

The trend of low breastfeeding rates increases the demand for infant milk formula (IMF) worldwide, but the use of IMF may be one of the causes of bacterial infections in infants. Complete sterility in the whole production line of IMF cannot be guaranteed; therefore, it is necessary to closely monitor the microbial content in the process. In the present study, an IMF powder production line based on the wet mixing process was sampled at 27 suspicious points in spring and summer to analyze the bacterial diversity by high-throughput sequencing. We found that 70 and 69 different bacterial phyla were present in spring and summer samples, respectively, with Proteobacteria and Firmicutes being the dominant phyla (>80% relative abundance). Moreover, 13 dominant genera each were present in spring (e. g., Pseudomonas and Lactococcus) and summer (e. g., Pseudomonas, Bacillus, and Streptococcus). Samples associated with workers showed higher bacterial species diversity (Shannon index) and richness (Chao1 index) in summer than in spring. The bacterial community composition showed high similarity between liquid milk after pasteurization and concentrated milk after evaporation. The potential bacterial pathogens were identified as Pseudomonas aeruginosa in spring and Acinetobacter baumannii in summer. Through retrospective analysis of the two opportunistic pathogens identified, it was found that the workshop environment was the potential contamination point in spring, whereas the auxiliary ingredients were the potential source of contamination in summer. The results highlight the effect of season on bacterial diversity associated with the production process of IMF and are useful in controlling the microbial quality and safety of infant dairy products.

Molecular characterization of bacteria and archaea in a bioaugmented zero-water exchange shrimp pond

In the zero-water exchange shrimp culture pond maintained with the application of indigenous bioaugmentor, low levels of total ammonia–nitrogen were reported, indicating the relevance of indigenous microbial communities. Sediments (0–5 cm layer) were sampled from the pond (85th day) and the bacterial and archaeal communities; specifically, the ammonia oxidizers (ammonia-oxidizing bacteria, ammonia-oxidizing archaea, and anaerobic ammonia-oxidizing bacteria) in the sediment metagenome of the pond were analysed using the 16S rRNA and functional genes. Bacterial and archaeal 16S rRNA genes showed the relative abundance of Delta-Proteobacteria and Bacteroidetes groups performing sulphur respiration and organic matter degradation, archaeal groups of anaerobic sulphur respiring Crenarchaeotae, and chemolithoautotrophic ammonia oxidizers belonging to Thaumarchaeota. The presence of these diverse bacterial and archaeal communities denotes their significant roles in the cycling the carbon, nitrogen, and sulphur thereby bringing out efficient bioremediation in the bioaugmented zero-water exchange shrimp culture pond. Similarly, the functional gene-specific study showed the predominance of Nitrosomonas sp. (ammonia-oxidizing bacteria), Nitrosopumilus maritimus (ammonia-oxidizing archaea), and Candidatus Kuenenia (anaerobic ammonia-oxidizing bacteria) in the system, which points to their importance in the removal of accumulated ammonia. Thus, this study paves the way for understanding the microbial communities, specifically the ammonia oxidizers responsible for maintaining healthy and optimal environmental conditions in the bioaugmented zero-water exchange shrimp culture pond.

Community composition of marine and brackish water ammonia-oxidizing consortia developed for aquaculture application

To mitigate the toxicity of ammonia in aquaculture systems, marine and brackish water ammonia-oxidizing bacterial consortia have been developed and are used for activation of nitrifying bioreactors integrated to recirculating aquaculture systems. To shed more light on to these biological entities, diversity of both the consortia were analyzed based on random cloning of 16S rRNA gene and ammonia-oxidizing bacterial specific amoA gene sequences. The dendrograms of representative clones on the basis of amplified ribosomal DNA restriction analysis generated 22 and 19 clusters for marine and brackish water nitrifying consortia, respectively. Phylogenetic analysis demonstrated the presence of various autotrophic nitrifiers belonging to α-, β- and γ-Proteobacteria, anaerobic ammonia oxidizers, heterotrophic denitrifiers, Bacteroidetes, and Actinobacteria. Distribution patterns of the organisms within the two consortia were determined using the software Geneious and diversity indices were investigated using Mega 5.0, VITCOMIC and Primer 7. The abundance of ammonia oxidizers was found in the order of 2.21 ± 0.25 × 10⁹ copies/g wet weight of marine consortium and 6.20 ± 0.23 × 10⁷ copies/g of brackish water consortium. Besides, marine ammonia-oxidizing consortium exhibited higher mean population diversity and Shannon Wiener diversity than the brackish water counterparts.

Subtractive proteomics to identify novel drug targets and reverse vaccinology for the development of chimeric vaccine against Acinetobacter baumannii

The emergence of drug-resistant Acinetobacter baumannii is the global health problem associated with high mortality and morbidity. Therefore it is high time to find a suitable therapeutics for this pathogen. In the present study, subtractive proteomics along with reverse vaccinology approaches were used to predict suitable therapeutics against A. baumannii. Using subtractive proteomics, we have identified promiscuous antigenic membrane proteins that contain the virulence factors, resistance factors and essentiality factor for this pathogenic bacteria. Selected promiscuous targeted membrane proteins were used for the design of chimeric-subunit vaccine with the help of reverse vaccinology. Available best tools and servers were used for the identification of MHC class I, II and B cell epitopes. All selected epitopes were further shortlisted computationally to know their immunogenicity, antigenicity, allergenicity, conservancy and toxicity potentials. Immunogenic predicted promiscuous peptides used for the development of chimeric subunit vaccine with immune-modulating adjuvants, linkers, and PADRE (Pan HLA-DR epitopes) amino acid sequence. Designed vaccine construct V4 also interact with the MHC, and TLR4/MD2 complex as confirm by docking and molecular dynamics simulation studies. Therefore designed vaccine construct V4 can be developed to control the host-pathogen interaction or infection caused by A. baumannii.

Latent environment allocation of microbial community data

As data for microbial community structures found in various environments has increased, studies have examined the relationship between environmental labels given to retrieved microbial samples and their community structures. However, because environments continuously change over time and space, mixed states of some environments and its effects on community formation should be considered, instead of evaluating effects of discrete environmental categories. Here we applied a hierarchical Bayesian model to paired datasets containing more than 30,000 samples of microbial community structures and sample description documents. From the training results, we extracted latent environmental topics that associate co-occurring microbes with co-occurring word sets among samples. Topics are the core elements of environmental mixtures and the visualization of topic-based samples clarifies the connections of various environments. Based on the model training results, we developed a web application, LEA (Latent Environment Allocation), which provides the way to evaluate typicality and heterogeneity of microbial communities in newly obtained samples without confining environmental categories to be compared. Because topics link words and microbes, LEA also enables to search samples semantically related to the query out of 30,000 microbiome samples.

In the past decade, microbiomes from various natural and human symbiotic environments have been thoroughly studied. However, our knowledge is limited as to what types of environments affect the structure of a microbial community. In the first place, how can we define “environments”, in particular, the environmental entities that are often continuously varying and difficult to discretely categorize? We assumed that environments could be represented from microbiome data because the structure of microbial communities reflect the state of the environment. We applied a probabilistic topic model to a dataset containing taxonomic composition data and natural language sample descriptions of >30,000 microbiome samples and extracted “latent environments” of the microbial communities, which are core elements of environmental mixtures. Integrating the training results of the model, we developed a web application to explore the microbiome universe and to place new metagenomic data on this universe like a global positioning system. Our tool shows what kinds of the environment naturally exist and are similar to each other on the perspective of the structural patterns of microbiome, and provides the way to evaluate the commonality and the heterogeneity of users’ microbiome samples.

VITCOMIC2: visualization tool for the phylogenetic composition of microbial communities based on 16S rRNA gene amplicons and metagenomic shotgun sequencing

Background

The 16S rRNA gene-based amplicon sequencing analysis is widely used to determine the taxonomic composition of microbial communities. Once the taxonomic composition of each community is obtained, evolutionary relationships among taxa are inferred by a phylogenetic tree. Thus, the combined representation of taxonomic composition and phylogenetic relationships among taxa is a powerful method for understanding microbial community structure; however, applying phylogenetic tree-based representation with information on the abundance of thousands or more taxa in each community is a difficult task. For this purpose, we previously developed the tool VITCOMIC (VIsualization tool for Taxonomic COmpositions of MIcrobial Community), which is based on the genome-sequenced microbes’ phylogenetic information. Here, we introduce VITCOMIC2, which incorporates substantive improvements over VITCOMIC that were necessary to address several issues associated with 16S rRNA gene-based analysis of microbial communities.

Results

We developed VITCOMIC2 to provide (i) sequence identity searches against broad reference taxa including uncultured taxa; (ii) normalization of 16S rRNA gene copy number differences among taxa; (iii) rapid sequence identity searches by applying the graphics processing unit-based sequence identity search tool CLAST; (iv) accurate taxonomic composition inference and nearly full-length 16S rRNA gene sequence reconstructions for metagenomic shotgun sequencing; and (v) an interactive user interface for simultaneous representation of the taxonomic composition of microbial communities and phylogenetic relationships among taxa. We validated the accuracy of processes (ii) and (iv) by using metagenomic shotgun sequencing data from a mock microbial community.

Conclusions

The improvements incorporated into VITCOMIC2 enable users to acquire an intuitive understanding of microbial community composition based on the 16S rRNA gene sequence data obtained from both metagenomic shotgun and amplicon sequencing.

Electronic supplementary material

The online version of this article (10.1186/s12918-018-0545-2) contains supplementary material, which is available to authorized users.

Thermophiles in the genomic era: Biodiversity, science, and applications

Thermophiles and hyperthermophiles are present in various regions of the Earth, including volcanic environments, hot springs, mud pots, fumaroles, geysers, coastal thermal springs, and even deep-sea hydrothermal vents. They are also found in man-made environments, such as heated compost facilities, reactors, and spray dryers. Thermophiles, hyperthermophiles, and their bioproducts facilitate various industrial, agricultural, and medicinal applications and offer potential solutions to environmental damages and the demand for biofuels. Intensified efforts to sequence the entire genome of hyperthermophiles and thermophiles are increasing rapidly, as evidenced by the fact that over 120 complete genome sequences of the hyperthermophiles Aquificae, Thermotogae, Crenarchaeota, and Euryarchaeota are now available. In this review, we summarise the major current applications of thermophiles and thermozymes. In addition, emphasis is placed on recent progress in understanding the biodiversity, genomes, transcriptomes, metagenomes, and single-cell sequencing of thermophiles in the genomic era.

Metagenomic analysis of the stool microbiome in patients receiving allogeneic stem cell transplantation: loss of diversity is associated with use of systemic antibiotics and more pronounced in gastrointestinal graft-versus-host disease

Next-generation sequencing of the hypervariable V3 region of the 16s rRNA gene isolated from serial stool specimens collected from 31 patients receiving allogeneic stem cell transplantation (SCT) was performed to elucidate variations in the composition of the intestinal microbiome in the course of allogeneic SCT. Metagenomic analysis was complemented by strain-specific enterococcal PCR and indirect assessment of bacterial load by liquid chromatography-tandem mass spectrometry of urinary indoxyl sulfate. At the time of admission, patients showed a predominance of commensal bacteria. After transplantation, a relative shift toward enterococci was observed, which was more pronounced under antibiotic prophylaxis and treatment of neutropenic infections. The shift was particularly prominent in patients that developed subsequently or suffered from active gastrointestinal (GI) graft-versus-host disease (GVHD). The mean proportion of enterococci in post-transplant stool specimens was 21% in patients who did not develop GI GVHD as compared with 46% in those that subsequently developed GI GVHD and 74% at the time of active GVHD. Enterococcal PCR confirmed predominance of Enterococcus faecium or both E. faecium and Enterococcus faecalis in these specimens. As a consequence of the loss of bacterial diversity, mean urinary indoxyl sulfate levels dropped from 42.5 ± 11 μmol/L to 11.8 ± 2.8 μmol/L in all post-transplant samples and to 3.5 ± 3 μmol/L in samples from patients with active GVHD. Our study reveals major microbiome shifts in the course of allogeneic SCT that occur in the period of antibiotic treatment but are more prominent in association with GI GVHD. Our data indicate early microbiome shifts and a loss of diversity of the intestinal microbiome that may affect intestinal inflammation in the setting of allogeneic SCT.

Design and experimental application of a novel non-degenerate universal primer set that amplifies prokaryotic 16S rRNA genes with a low possibility to amplify eukaryotic rRNA genes

The deep sequencing of 16S rRNA genes amplified by universal primers has revolutionized our understanding of microbial communities by allowing the characterization of the diversity of the uncultured majority. However, some universal primers also amplify eukaryotic rRNA genes, leading to a decrease in the efficiency of sequencing of prokaryotic 16S rRNA genes with possible mischaracterization of the diversity in the microbial community. In this study, we compared 16S rRNA gene sequences from genome-sequenced strains and identified candidates for non-degenerate universal primers that could be used for the amplification of prokaryotic 16S rRNA genes. The 50 identified candidates were investigated to calculate their coverage for prokaryotic and eukaryotic rRNA genes, including those from uncultured taxa and eukaryotic organelles, and a novel universal primer set, 342F-806R, covering many prokaryotic, but not eukaryotic, rRNA genes was identified. This primer set was validated by the amplification of 16S rRNA genes from a soil metagenomic sample and subsequent pyrosequencing using the Roche 454 platform. The same sample was also used for pyrosequencing of the amplicons by employing a commonly used primer set, 338F-533R, and for shotgun metagenomic sequencing using the Illumina platform. Our comparison of the taxonomic compositions inferred by the three sequencing experiments indicated that the non-degenerate 342F-806R primer set can characterize the taxonomic composition of the microbial community without substantial bias, and is highly expected to be applicable to the analysis of a wide variety of microbial communities.

Bacterial diversity of surface sand samples from the Gobi and Taklamaken deserts

Arid regions represent nearly 30 % of the Earth's terrestrial surface, but their microbial biodiversity is not yet well characterized. The surface sands of deserts, a subset of arid regions, are generally subjected to large temperature fluctuations plus high UV light exposure and are low in organic matter. We examined surface sand samples from the Taklamaken (China, three samples) and Gobi (Mongolia, two samples) deserts, using pyrosequencing of PCR-amplified 16S V1/V2 rDNA sequences from total extracted DNA in order to gain an assessment of the bacterial population diversity. In total, 4,088 OTUs (using ≥97 % sequence similarity levels), with Chao1 estimates varying from 1,172 to 2,425 OTUs per sample, were discernable. These could be grouped into 102 families belonging to 15 phyla, with OTUs belonging to the Firmicutes, Proteobacteria, Bacteroidetes, and Actinobacteria phyla being the most abundant. The bacterial population composition was statistically different among the samples, though members from 30 genera were found to be common among the five samples. An increase in phylotype numbers with increasing C/N ratio was noted, suggesting a possible role in the bacterial richness of these desert sand environments. Our results imply an unexpectedly large bacterial diversity residing in the harsh environment of these two Asian deserts, worthy of further investigation.

Characteristics of microbial communities in crustal fluids in a deep-sea hydrothermal field of the suiyo seamount

To directly access the sub-seafloor microbial communities, seafloor drilling has been done in a deep-sea hydrothermal field of the Suiyo Seamount, Izu-Bonin Arc, Western Pacific. In the present study, crustal fluids were collected from the boreholes, and the bacterial and archaeal communities in the fluids were investigated by culture-independent molecular analysis based on 16S rRNA gene sequences. Bottom seawater, sands, rocks, sulfide mound, and chimneys were also collected around the boreholes and analyzed for comparisons. Comprehensive analysis revealed the characteristics of the microbial community composition in the crustal fluids. Phylotypes closely related to cultured species, e.g., Alteromonas, Halomonas, Marinobacter, were relatively abundant in some crustal fluid samples, whereas the phylotypes related to Pelagibacter and the SUP05-group were relatively abundant in the seawater samples. Phylotypes related to other uncultured environmental clones in Alphaproteobacteria and Gammaproteobacteria were relatively abundant in the sand, rock, sulfide mound, and chimney samples. Furthermore, comparative analysis with previous studies of the Suiyo Seamount crustal fluids indicates the change in the microbial community composition for 3 years. Our results provide novel insights into the characteristics of the microbial communities in crustal fluids beneath a deep-sea hydrothermal field.

Phylogenetic affiliation of SSU rRNA genes generated by massively parallel sequencing: new insights into the freshwater protist diversity

Recent advances in next-generation sequencing (NGS) technologies spur progress in determining the microbial diversity in various ecosystems by highlighting, for example, the rare biosphere. Currently, high-throughput pyrotag sequencing of PCR-amplified SSU rRNA gene regions is mainly used to characterize bacterial and archaeal communities, and rarely to characterize protist communities. In addition, although taxonomic assessment through phylogeny is considered as the most robust approach, similarity and probabilistic approaches remain the most commonly used for taxonomic affiliation. In a first part of this work, a tree-based method was compared with different approaches of taxonomic affiliation (BLAST and RDP) of 18S rRNA gene sequences and was shown to be the most accurate for near full-length sequences and for 400 bp amplicons, with the exception of amplicons covering the V5-V6 region. Secondly, the applicability of this method was tested by running a full scale test using an original pyrosequencing dataset of 18S rRNA genes of small lacustrine protists (0.2-5 µm) from eight freshwater ecosystems. Our results revealed that i) fewer than 5% of the operational taxonomic units (OTUs) identified through clustering and phylogenetic affiliation had been previously detected in lakes, based on comparison to sequence in public databases; ii) the sequencing depth provided by the NGS coupled with a phylogenetic approach allowed to shed light on clades of freshwater protists rarely or never detected with classical molecular ecology approaches; and iii) phylogenetic methods are more robust in describing the structuring of under-studied or highly divergent populations. More precisely, new putative clades belonging to Mamiellophyceae, Foraminifera, Dictyochophyceae and Euglenida were detected. Beyond the study of protists, these results illustrate that the tree-based approach for NGS based diversity characterization allows an in-depth description of microbial communities including taxonomic profiling, community structuring and the description of clades of any microorganisms (protists, Bacteria and Archaea).

Metagenomic profiles of free-living archaea, bacteria and small eukaryotes in coastal areas of Sichang island, Thailand

Background

Tha Wang and Tham Phang coasts, though situated at similar oceanographic positions on Sichang island, Chonburi province, Thailand, are different in bay geography and amount of municipal disturbances. These affect the marine ecosystems. The study used metagenomics combined with 16S and 18S rDNA pyrosequencing to identify types and distributions of archaea, bacteria, fungi and small eukaryotes of sizes ranges 0.45 and ~30 μm.

Results

Following the open bay geography and minimal municipal sewages, Tham Phang coast showed the cleaner water properties, described by color, salinity, pH, conductivity and percent dissolved oxygen. The 16S and 18S rDNA metagenomic profiles for Tha Wang and Tham Phang coasts revealed many differences, highlighting by low Lennon and Yue & Clayton theta similarity indices (66.03-73.03% for 16S rDNA profiles, 2.85-25.38% for 18S rDNA profiles). For 16S rDNA, the percent compositions of species belonging to Proteobacteria, Bacteroidetes, Cyanobacteria, Firmicutes, Verrucomicrobia, Gammatimonadetes, Tenericutes, Acidobacteria, Spirochaetes, Chlamydiae, Euryarchaeota, Nitrospirae, Planctomycetes, Thermotogae and Aquificae were higher or distinctly present in Tha Wang. In Tham Phang, except Actinobacteria, the fewer number of prokaryotic species existed. For 18S rDNA, fungi represented 74.745% of the species in Tha Wang, whereas only 6.728% in Tham Phang. Basidiomycota (71.157%) and Ascomycota (3.060%) were the major phyla in Tha Wang. Indeed, Tha Wang-to-Tham Phang percent composition ratios for fungi Basidiomycota and Chytridiomycota were 1264.701 and 25.422, respectively. In Tham Phang, Brachiopoda (lamp shells) and Mollusca (snails) accounted for 80.380% of the 18S rDNA species detected, and their proportions were approximately tenfold greater than those in Tha Wang. Overall, coastal Tham Phang comprised abundant animal species.

Conclusions

Tha Wang contained numerous archaea, bacteria and fungi, many of which could synthesize useful biotechnology gas and enzymes that could also function in high-saline and high-temperature conditions. Tham Phang contained less abundant archaea, bacteria and fungi, and the majority of the extracted metagenomes belonged to animal kingdom. Many microorganisms in Tham Phang were essential for nutrient-recycling and pharmaceuticals, for instances, Streptomyces, Pennicilium and Saccharomyces. Together, the study provided metagenomic profiles of free-living prokaryotes and eukaryotes in coastal areas of Sichang island.

Assessment of Microbial Richness in Pelagic Sediment of Andaman Sea by Bacterial Tag Encoded FLX Titanium Amplicon Pyrosequencing (bTEFAP)

Microbial diversity of 1,000 m deep pelagic sediment from off Coast of Andaman Sea was analyzed by a culture independent technique, bacterial tag encoded FLX titanium amplicon pyrosequencing. The hypervariable region of small subunit ribosomal rRNA gene covering V6-V9, was amplified from the metagenomic DNA and sequenced. We obtained 19,271 reads, of which 18,206 high quality sequences were subjected to diversity analysis. A total of 305 operational taxonomic units (OTUs) were obtained corresponding to the members of firmicutes, proteobacteria, plantomycetes, actinobacteria, chloroflexi, bacteroidetes, and verucomicrobium. Firmicutes was the predominant phylum, which was largely represented with the family bacillaceae. More than 44 % of sequence reads could not be classified up to the species level and more than 14 % of the reads could not be assigned to any genus. Thus, the data indicates the possibility for the presence of uncultivable or unidentified novel bacterial species. In addition, the community structure identified in this study significantly differs with other reports from marine sediments.

Estimation of bacterial diversity using next generation sequencing of 16S rDNA: a comparison of different workflows

Background

Next generation sequencing (NGS) enables a more comprehensive analysis of bacterial diversity from complex environmental samples. NGS data can be analysed using a variety of workflows. We test several simple and complex workflows, including frequently used as well as recently published tools, and report on their respective accuracy and efficiency under various conditions covering different sequence lengths, number of sequences and real world experimental data from rhizobacterial populations of glyphosate-tolerant maize treated or untreated with two different herbicides representative of differential diversity studies.

Results

Alignment and distance calculations affect OTU estimations, and multiple sequence alignment exerts a major impact on the computational time needed. Generally speaking, most of the analyses produced consistent results that may be used to assess differential diversity changes, however, dataset characteristics dictate which workflow should be preferred in each case.

Conclusions

When estimating bacterial diversity, ESPRIT as well as the web-based workflow, RDP pyrosequencing pipeline, produced good results in all circumstances, however, its computational requirements can make method-combination workflows more attractive, depending on sequence variability, number and length.

STITCH: algorithm to splice, trim, identify, track, and capture the uniqueness of 16S rRNAs sequence pairs using public or in-house database

A comparison of variable regions within the 16S rRNA gene is widely used to characterize relationships between bacteria and to identify phylogenetic affiliation of unknown bacteria. In environmental studies, polymerase chain reaction amplification of 16S rRNA followed by cloning and sequencing of numerous individual clones is an extensively used molecular method for elucidating microbial diversity. The sequencing process typically utilizes a forward and reverse primer pair to produce two partial reads (~700 to 800 base pairs each) that overlap and in total cover a large region of the full 16S rRNA sequence (~1.5 k base). In a typical application, this approach rapidly generates very large numbers of 16S rRNA datasets that can overwhelm manual processing efforts leading to both delays and errors. In particular, the approach presents two computational challenges: (1) the assembly of a composite sequence from the two partial reads and (2) the subsequent appropriate identification of the organism represented by the newly sequenced clones. Herein, we describe a software package, search, trim, identify, track, and capture the uniqueness of 16S rRNAs using public and in-house database (STITCH), which offers automated sequence pair splicing and genetic identification, thus simplifying the computationally intensive analysis of large sequencing libraries. The STITCH software is freely accessible over the Internet at: http://prion.bchs.uh.edu/stitch/.