Two-stage clustering (TSC): a pipeline for selecting operational taxonomic units for the high-throughput sequencing of PCR amplicons

Dynamic changes in microbial communities and flavor during different fermentation stages of proso millet Baijiu, a new product from Shanxi light-flavored Baijiu

Introduction

Proso millet, a high-quality fermentation material used for Chinese yellow wine production, can produce special flavored substances; however, its role in improving the flavor and altering microbial communities of light-flavored Baijiu during fermentation remain unknown. Thus, we aimed to investigate the effect of proso millet on improving the flavor of light-flavored Baijiu and altering microbial communities during different fermentation stages.

Methods

The dynamic changes in the microbial communities and flavor of proso millet (50%) + sorghum (50%) mixed fermentation samples were analyzed through intermittent sampling on days 7, 14, 21, and 28 of the fermentation process. Microbial high-throughput sequencing and the analysis of flavor characteristics were conducted through 16S DNA/ ITS amplicon sequencing and gas chromatography (multi-capillary column)-ion mobility spectrometry, respectively.

Results

Proso millet significantly changed the core flavor compound composition of traditional light-flavored Baijiu from ethyl acetate, ethyl hexanoate, ethyl hexanoate dimer, ethyl butanoate, ethyl lactate, and butyl acetate to oct-2-ene, 2-butanol, propyl propanoate, 2-pentenal, and 4-methylpentanal. The amplicon sequencing analysis revealed that the alpha diversity parameters of bacterial and fungal communities, including the Chao1, Pielou_e, Shannon, and Simpson indices, for proso millet-sorghum mixed fermentation samples were significantly higher than those for sorghum fermentation samples (p < 0.05). Of the 40 most significant microbial genera in two treatments, proso millet significantly increased the abundance of 12 bacterial and 18 fungal genera. Among the 40 most significant bacterial and fungal species, 23 bacterial species belonged to the Lactobacillus genus, whereas the 30 primary fungal species belonged to 28 different genera. The analysis of the relationship between microbial changes and the main flavor compounds of light-flavored Baijiu showed that bacteria from the Weissella, Acinetobacter, Bacteroides, Psychrobacter, Pseudarthrobacter, Lactococcus, Chloroplast, Saccharopolyspora, Psychrobacter, Saccharopolyspora, Pseudonocardiaceae, Bacteroides genera and fungi from the Thermoascus, Aspergillus, Pichia, Rhizomucor, Papiliotrema, Hyphopichia, and Mucor genera significantly inhibited the synthesis of ethyl hexanoate, ethyl butanoate, ethyl lactate ethyl lactate, and butyl acetate but increased the synthesis of ethyl acetate (p < 0.05). Moreover, these microbes exhibited a significantly greater abundance in proso millet-sorghum mixed fermentation samples than in sorghum samples. The synthesis of special flavored compounds in proso millet Baijiu was significantly positively correlated with the presence of fungi from the Rhizopus, Papiliotrema, Wickerhamomyces, Aspergillus, and Thermoascus genera but negative correlated with the presence of bacteria from the Weissella, Acinetobacter, Psychrobacter, Pseudarthrobacter, Bacteroides, and Saccharopolyspora genera. Regarding ethanol content, the low alcohol content of Fenjiu may be due to the significantly high abundance of fungi from the Psathyrella genus and bacteria from the Staphylococcus, Kroppenstedtia, Brevibacterium, and Acetobacter genera during fermentation. In summary, proso millet significantly altered the flavor of light-flavored Baijiu by inducing the formation of a special microbial community; however, it did not increase alcohol concentration.

Discussion

This study lays the foundation for future research on Baijiu fermentation. Additionally, the study findings may help improve the production efficiency and elevate the quality and flavor of the final product.

Collateral Damage in the Human Gut Microbiome - Blastocystis Is Significantly Less Prevalent in an Antibiotic-Treated Adult Population Compared to Non-Antibiotic Treated Controls

Antibiotics can drive the rapid loss of non-target, phylogenetically diverse microorganisms that inhabit the human gut. This so-called “collateral damage” has myriad consequences for host health and antibiotic mediated changes to the gut microbiota have been implicated in the aetiology of many chronic diseases. To date, studies have largely focused on how antibiotics affect the bacterial fraction of the gut microbiome and their impact on non-bacterial members, including prevalent eukaryal species, such as Blastocystis, remains largely unknown. Here we assessed the prevalence and diversity of Blastocystis in an elderly adult group that were in receipt of antibiotics (n = 86) and an equivalent non-antibiotic treated group (n = 88) using a PCR-based approach. This analysis revealed that although similar subtypes were present in both groups, Blastocystis was significantly less prevalent in the antibiotic-treated group (16%) compared to non-antibiotic treated controls (55%); Fisher’s Exact test, p < 0.0001). Given that antibiotics target structures and molecules of prokaryotic cells to kill or inhibit bacterial populations, the most likely explanation for differences in prevalence between both groups is due to secondary extinctions owing to the potential dependence of Blastocystis on bacteria present in the gut microbiome that were negatively affected by antibiotic treatment. Although further work is required to explore this hypothesis in greater detail, these data clearly show that Blastocystis prevalence in human populations is negatively associated with antibiotic treatment. This finding may be relevant to explaining patterns of variation for this microorganism in different human populations and cohorts of interest.

Tetrabromobisphenol A alters soil microbial community via selective antibacterial activity

Tetrabromobisphenol A (TBBPA) is the most widely used brominated flame retardant. Most studies regarding TBBPA have concentrated on its occurrence, distribution, toxicity and degradation in the environment. However, little is known about its ecological effects on soil microbial communities. In this study, we investigated the effect of TBBPA on soil microbial community. Overall, the data suggested that the growth and composition of soil microorganisms were correlated to the TBBPA concentration and exposure time. Phospholipid-derived fatty acid analysis (PLFAs) showed that significant microbial growth inhibitions were 46.1% and 46.9% in 40 mg/kg TBBPA-treated soils after 45-day incubation under aerobic and anaerobic conditions, respectively. Results of PLFAs and llumina sequencing indicated that TBBPA mainly inhibited Gram-positive bacteria, but not Gram-negative bacteria. The selective antibacterial activity of TBBPA toward Gram-positive bacteria was further confirmed in pure bacteria cultures. These data suggested that, in addition to their effect on microbial growth and composition, TBBPA may affect the microbial ecology. Additional research should be carried out to identify the ecological risk of TBBPA in soil.

Specificity in Arabidopsis thaliana recruitment of root fungal communities from soil and rhizosphere

Biotic and abiotic conditions in soil pose major constraints on growth and reproductive success of plants. Fungi are important agents in plant soil interactions but the belowground mycobiota associated with plants remains poorly understood. We grew one genotype each from Sweden and Italy of the widely-studied plant model Arabidopsis thaliana. Plants were grown under controlled conditions in organic topsoil local to the Swedish genotype, and harvested after ten weeks. Total DNA was extracted from three belowground compartments: endosphere (sonicated roots), rhizosphere and bulk soil, and fungal communities were characterized from each by amplification and sequencing of the fungal barcode region ITS2. Fungal species diversity was found to decrease from bulk soil to rhizosphere to endosphere. A significant effect of plant genotype on fungal community composition was detected only in the endosphere compartment. Despite A. thaliana being a non-mycorrhizal plant, it hosts a number of known mycorrhiza fungi in its endosphere compartment, which is also colonized by endophytic, pathogenic and saprotrophic fungi. Species in the Archaeorhizomycetes were most abundant in rhizosphere samples suggesting an adaptation to environments with high nutrient turnover for some of these species. We conclude that A. thaliana endosphere fungal communities represent a selected subset of fungi recruited from soil and that plant genotype has small but significant quantitative and qualitative effects on these communities.

[Longitudinal analysis of vaginal microbiota in women with recurrent vulvovaginal candidiasis]

OBJECTIVE

To investigate the vaginal flora in patients with recurrent vulvovaginal candidiasis (RVVC).

METHODS

Vaginal swabs were collected at different time points from 6 RVVC patients and 5 healthy women of child-bearing age. The dynamic changes, microbiota composition, alpha diversity and beta diversity in the two groups were assessed by analyzing the 16S rRNA V4 hypervariable region amplified from the total genomic DNA from the swabs.

RESULTS

Lactobacillus was the predominant species in healthy women with similar proportions of L.iners and L.crispatus; small proportions of Gardnerella, Prevotella and other genus were also detected. In some healthy women, the vaginal flora showed a high relative abundance of anaerobic bacteria such as Gardnerella, Prevotella, Atopobium, Sneathia. Compared with the healthy women, patients with RVVC showed a significantly reduced diversity of vaginal flora, where L.iners was the predominant species and the content of L.crispatus decreased significantly. In healthy women, the vaginal flora fluctuated with the menstrual cycle, and the fluctuation was the most prominent during menstruation; the dominant species either alternated regularly or maintain an absolute superiority in the menstrual cycle. The vaginal flora showed attenuated fluctuation in women with RVVC, were highly conserved within the menstrual cycle, and maintained a similar composition in the episodes and intermittent periods.

CONCLUSION

The vaginal flora of RVVC patients do not undergo regular variations with the menstrual cycle and shows a similar composition between the episodes and intermittent periods. Promoting the production of L.iners or inhibiting the colonization of L.crispatus to restore the composition of the vaginal flora may help in the treatment of RVVC.

Erratum to: Stability of operational taxonomic units: an important but neglected property for analyzing microbial diversity

[This corrects the article DOI: 10.1186/s40168-015-0081-x.].

Composition and temporal stability of the gut microbiota in older persons

The composition and function of the human gut microbiota has been linked to health and disease. We previously identified correlations between habitual diet, microbiota composition gradients and health gradients in an unstratified cohort of 178 elderly subjects. To refine our understanding of diet–microbiota associations and differential taxon abundance, we adapted an iterative bi-clustering algorithm (iterative binary bclustering of gene sets (iBBiG)) and applied it to microbiota composition data from 732 faecal samples from 371 ELDERMET cohort subjects, including longitudinal samples. We thus identified distinctive microbiota configurations associated with ageing in both community and long-stay residential care elderly subjects. Mixed-taxa populations were identified that had clinically distinct associations. Microbiota temporal instability was observed in both community-dwelling and long-term care subjects, particularly in those with low initial microbiota diversity. However, the stability of the microbiota of subjects had little impact on the directional change of the microbiota as observed for long-stay subjects who display a gradual shift away from their initial microbiota. This was not observed in community-dwelling subjects. This directional change was associated with duration in long-stay. Changes in these bacterial populations represent the loss of the health-associated and youth-associated microbiota components and gain of an elderly associated microbiota. Interestingly, community-associated microbiota configurations were impacted more by the use of antibiotics than the microbiota of individuals in long-term care, as the community-associated microbiota showed more loss but also more recovery following antibiotic treatment. This improved definition of gut microbiota composition patterns in the elderly will better inform the design of dietary or antibiotic interventions targeting the gut microbiota.

Toward accurate molecular identification of species in complex environmental samples: testing the performance of sequence filtering and clustering methods

Metabarcoding has the potential to become a rapid, sensitive, and effective approach for identifying species in complex environmental samples. Accurate molecular identification of species depends on the ability to generate operational taxonomic units (OTUs) that correspond to biological species. Due to the sometimes enormous estimates of biodiversity using this method, there is a great need to test the efficacy of data analysis methods used to derive OTUs. Here, we evaluate the performance of various methods for clustering length variable 18S amplicons from complex samples into OTUs using a mock community and a natural community of zooplankton species. We compare analytic procedures consisting of a combination of (1) stringent and relaxed data filtering, (2) singleton sequences included and removed, (3) three commonly used clustering algorithms (mothur, UCLUST, and UPARSE), and (4) three methods of treating alignment gaps when calculating sequence divergence. Depending on the combination of methods used, the number of OTUs varied by nearly two orders of magnitude for the mock community (60–5068 OTUs) and three orders of magnitude for the natural community (22–22191 OTUs). The use of relaxed filtering and the inclusion of singletons greatly inflated OTU numbers without increasing the ability to recover species. Our results also suggest that the method used to treat gaps when calculating sequence divergence can have a great impact on the number of OTUs. Our findings are particularly relevant to studies that cover taxonomically diverse species and employ markers such as rRNA genes in which length variation is extensive.

Different Dynamic Patterns of β-Lactams, Quinolones, Glycopeptides and Macrolides on Mouse Gut Microbial Diversity

The adverse impact of antibiotics on the gut microbiota has attracted extensive interest, particularly due to the development of microbiome research techniques in recent years. However, a direct comparison of the dynamic effects of various types of antibiotics using the same animal model has not been available. In the present study, we selected six antibiotics from four categories with the broadest clinical usage, namely, β-lactams (Ceftriaxone Sodium, Cefoperazone/Sulbactam and meropenem), quinolones (ofloxacin), glycopeptides (vancomycin), and macrolides (azithromycin), to treat BALB/c mice. Stool samples were collected during and after the administration of antibiotics, and microbial diversity was analyzed through Illumina sequencing and bioinformatics analyses using QIIME. Both α and β diversity analyses showed that ceftriaxone sodium, cefoperazone/sulbactam, meropenem and vancomycin changed the gut microbiota dramatically by the second day of antibiotic administration whereas the influence of ofloxacin was trivial. Azithromycin clearly changed the gut microbiota but much less than vancomycin and the β-lactams. In general, the community changes induced by the three β-lactam antibiotics showed consistency in inhibiting Papillibacter, Prevotella and Alistipes while inducing massive growth of Clostridium. The low diversity and high Clostridium level might be an important cause of Clostridium difficile infection after usage of β-lactams. Vancomycin was unique in that it inhibited Firmicutes, mainly the genus Clostridium. On the other hand, it induced the growth of Escherichia and effect lasted for months afterward. Azithromycin and meropenem induced the growth of Enterococcus. These findings will be useful for understanding the potential adverse effects of antibiotics on the gut microbiome and ensuring their better usage.

First insights into the microbial diversity in the omasum and reticulum of bovine using Illumina sequencing

The digestive systems of mammals harbor a complex gut microbiome, comprising bacteria and other microorganisms that confer metabolic and immunological benefits to the host. Ruminants that digest plant-based foods have a four-compartment stomach consisting of the rumen, reticulum, omasum, and abomasum. The microorganisms in the stomach are essential for providing the host with critical nutrients. However, the majority of these microorganisms are unknown species. The microbiome of the stomach is diverse, and the majority of these organisms cannot be cultured. Next-generation sequencing (NGS) combined with bioinformatic analysis tools have allowed the dissection of the composition of the microbiome in samples collected from a specific environment. In this study, for the first time, the bacterial composition in two compartments, the reticulum and the omasum, of bovine were analyzed using a metagenomic approach and compared to the bacterial composition of the rumen. These data will assist in understanding the biology of ruminants and benefit the agricultural industry. The diversity and composition of the bacterial community in samples collected from the rumen, reticulum, and omasum of bovines in the Changchun Region of Northeast China were analyzed by sequencing the V3 region of the 16S rRNA gene using a barcoded Illumina paired-end sequencing technique, and the primary composition of the microbiome in the rumen, reticulum, and omasum of the bovines was determined. These microbiomes contained 17 phyla and 107 genera in all three samples. Five phyla, Bacteroidetes, Firmicutes, Proteobacteria, Spirochaetes, and Lentisphaerae, were the most abundant taxonomic groups. Additionally, the different stomach compartments harbored different compositions of the microorganisms.

Stability of operational taxonomic units: an important but neglected property for analyzing microbial diversity

BACKGROUND

The operational taxonomic unit (OTU) is widely used in microbial ecology. Reproducibility in microbial ecology research depends on the reliability of OTU-based 16S ribosomal subunit RNA (rRNA) analyses.

RESULTS

Here, we report that many hierarchical and greedy clustering methods produce unstable OTUs, with membership that depends on the number of sequences clustered. If OTUs are regenerated with additional sequences or samples, sequences originally assigned to a given OTU can be split into different OTUs. Alternatively, sequences assigned to different OTUs can be merged into a single OTU. This OTU instability affects alpha-diversity analyses such as rarefaction curves, beta-diversity analyses such as distance-based ordination (for example, Principal Coordinate Analysis (PCoA)), and the identification of differentially represented OTUs. Our results show that the proportion of unstable OTUs varies for different clustering methods. We found that the closed-reference method is the only one that produces completely stable OTUs, with the caveat that sequences that do not match a pre-existing reference sequence collection are discarded.

CONCLUSIONS

As a compromise to the factors listed above, we propose using an open-reference method to enhance OTU stability. This type of method clusters sequences against a database and includes unmatched sequences by clustering them via a relatively stable de novo clustering method. OTU stability is an important consideration when analyzing microbial diversity and is a feature that should be taken into account during the development of novel OTU clustering methods.

Diverse vaginal microbiomes in reproductive-age women with vulvovaginal candidiasis

Vulvovaginal candidiasis (VVC) is one of the most prevalent vaginal infectious diseases, and there are controversial reports regarding the diversity of the associated vaginal microbiota. We determined the vaginal microbial community in patients with VVC, bacterial vaginosis (BV), and mixed infection of VVC and BV using Illumina sequencing of 16S rRNA tags. Our results revealed for the first time the highly variable patterns of the vaginal microbiome from VVC patients. In general, the alpha-diversity results of species richness and evenness showed the following order: normal control < VVC only < mixed BV and VVC infection < BV only. The beta-diversity comparison of community structures also showed an intermediate composition of VVC between the control and BV samples. A detailed comparison showed that, although the control and BV communities had typical patterns, the vaginal microbiota of VVC is complex. The mixed BV and VVC infection group showed a unique pattern, with a relatively higher abundance of Lactobacillus than the BV group and higher abundance of Prevotella, Gardnerella, and Atopobium than the normal control. In contrast, the VVC-only group could not be described by any single profile, ranging from a community structure similar to the normal control (predominated with Lactobacillus) to BV-like community structures (abundant with Gardnerella and Atopobium). Treatment of VVC resulted in inconsistent changes of the vaginal microbiota, with four BV/VVC samples recovering to a higher Lactobacillus level, whereas many VVC-only patients did not. These results will be useful for future studies on the role of vaginal microbiota in VVC and related infectious diseases.

Comparison of microbial diversity determined with the same variable tag sequence extracted from two different PCR amplicons

BACKGROUND

Deep sequencing of the variable region of 16S rRNA genes has become the predominant tool for studying microbial ecology. As sequencing datasets have accumulated, meta-analysis of sequences obtained with different variable 16S rRNA gene targets and by different sequencing methods has become an intriguing prospect that remains to be evaluated experimentally.

RESULTS

We amplified a group of fecal samples using both V4F-V6R and V6F-V6R primer sets, excised the same V6 fragment from the two sets of Illumina sequencing data, and compared the resulting data in terms of the α-diversity, β-diversity, and community structure. Principal component analysis (PCA) comparing the microbial community structures of different datasets, including those with simulated sequencing errors, was very reliable. Procrustes analysis showed a high degree of concordance between the different datasets for both abundance-weighted and binary Jaccard distances (P < 0.05), and a meta-analysis of individual datasets resulted in similar conclusions. The Shannon's diversity index was consistent as well, with comparable values obtained for the different datasets and for the meta-analysis of different datasets. In contrast, richness estimators (OTU and Chao) varied significantly, and the meta-analysis of richness estimators was also biased. The community structures of the two datasets were obviously different and led to significant changes in the biomarkers identified by the LEfSe statistical tool.

CONCLUSIONS

Our results suggest that beta-diversity analysis and Shannon's diversity are relatively reliable for meta-analysis, while community structures and biomarkers are less consistent. These results should be useful for future meta-analyses of microbiomes from different data sources.

Illumina sequencing of 16S rRNA tag revealed spatial variations of bacterial communities in a mangrove wetland

The microbial community plays an essential role in the high productivity in mangrove wetlands. A proper understanding of the spatial variations of microbial communities will provide clues about the underline mechanisms that structure microbial groups and the isolation of bacterial strains of interest. In the present study, the diversity and composition of the bacterial community in sediments collected from four locations, namely mudflat, edge, bulk, and rhizosphere, within the Mai Po Ramsar Wetland in Hong Kong, SAR, China were compared using the barcoded Illumina paired-end sequencing technique. Rarefaction results showed that the bulk sediment inside the mature mangrove forest had the highest bacterial α-diversity, while the mudflat sediment without vegetation had the lowest. The comparison of β-diversity using principal component analysis and principal coordinate analysis with UniFrac metrics both showed that the spatial effects on bacterial communities were significant. All sediment samples could be clustered into two major groups, inner (bulk and rhizosphere sediments collected inside the mangrove forest) and outer mangrove sediments (the sediments collected at the mudflat and the edge of the mangrove forest). With the linear discriminate analysis scores larger than 3, four phyla, namely Actinobacteria, Acidobacteria, Nitrospirae, and Verrucomicrobia, were enriched in the nutrient-rich inner mangrove sediments, while abundances of Proteobacteria and Deferribacterias were higher in outer mangrove sediments. The rhizosphere effect of mangrove plants was also significant, which had a lower α-diversity, a higher amount of Nitrospirae, and a lower abundance of Proteobacteria than the bulk sediment nearby.

Comparison of the levels of bacterial diversity in freshwater, intertidal wetland, and marine sediments by using millions of illumina tags

Sediment, a special realm in aquatic environments, has high microbial diversity. While there are numerous reports about the microbial community in marine sediment, freshwater and intertidal sediment communities have been overlooked. The present study determined millions of Illumina reads for a comparison of bacterial communities in freshwater, intertidal wetland, and marine sediments along Pearl River, China, using a technically consistent approach. Our results show that both taxon richness and evenness were the highest in freshwater sediment, medium in intertidal sediment, and lowest in marine sediment. The high number of sequences allowed the determination of a wide variety of bacterial lineages in all sediments for reliable statistical analyses. Principal component analysis showed that the three types of communities could be well separated from phylum to operational taxonomy unit (OTU) levels, and the OTUs from abundant to rare showed satisfactory resolutions. Statistical analysis (LEfSe) demonstrated that the freshwater sediment was enriched with Acidobacteria, Nitrospira, Verrucomicrobia, Alphaproteobacteria, and Betaproteobacteria. The intertidal sediment had a unique community with diverse primary producers (such as Chloroflexi, Bacillariophyta, Gammaproteobacteria, and Epsilonproteobacteria) as well as saprophytic microbes (such as Actinomycetales, Bacteroidetes, and Firmicutes). The marine sediment had a higher abundance of Gammaproteobacteria and Deltaproteobacteria, which were mainly involved in sulfate reduction in anaerobic conditions. These results are helpful for a systematic understanding of bacterial communities in natural sediment environments.

Comparing clustering and pre-processing in taxonomy analysis

MOTIVATION

Massively parallel sequencing allows for rapid sequencing of large numbers of sequences in just a single run. Thus, 16S ribosomal RNA (rRNA) amplicon sequencing of complex microbial communities has become possible. The sequenced 16S rRNA fragments (reads) are clustered into operational taxonomic units and taxonomic categories are assigned. Recent reports suggest that data pre-processing should be performed before clustering. We assessed combinations of data pre-processing steps and clustering algorithms on cluster accuracy for oral microbial sequence data.

RESULTS

The number of clusters varied up to two orders of magnitude depending on pre-processing. Pre-processing using both denoising and chimera checking resulted in a number of clusters that was closest to the number of species in the mock dataset (25 versus 15). Based on run time, purity and normalized mutual information, we could not identify a single best clustering algorithm. The differences in clustering accuracy among the algorithms after the same pre-processing were minor compared with the differences in accuracy among different pre-processing steps.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

CONTACT

bonder.m.j@gmail.com or b.brandt@acta.nl