Taxonomic precision of different hypervariable regions of 16S rRNA gene and annotation methods for functional bacterial groups in biological wastewater treatment

Fine-scale characterization of the soybean rhizosphere microbiome via synthetic long reads and avidity sequencing

BACKGROUND

The rhizosphere microbiome displays structural and functional dynamism driven by plant, microbial, and environmental factors. While such plasticity is a well-evidenced determinant of host health, individual and community-level microbial activity within the rhizosphere remain poorly understood, due in part to the insufficient taxonomic resolution achieved through traditional marker gene amplicon sequencing. This limitation necessitates more advanced approaches (e.g., long-read sequencing) to derive ecological inferences with practical application. To this end, the present study coupled synthetic long-read technology with avidity sequencing to investigate eukaryotic and prokaryotic microbiome dynamics within the soybean (Glycine max) rhizosphere under field conditions.

RESULTS

Synthetic long-read sequencing permitted de novo reconstruction of the entire 18S-ITS1-ITS2 region of the eukaryotic rRNA operon as well as all nine hypervariable regions of the 16S rRNA gene. All full-length, mapped eukaryotic amplicon sequence variants displayed genus-level classification, and 44.77% achieved species-level classification. The resultant eukaryotic microbiome encompassed five kingdoms (19 genera) of protists in addition to fungi - a depth unattainable with conventional short-read methods. In the prokaryotic fraction, every full-length, mapped amplicon sequence variant was resolved at the species level, and 23.13% at the strain level. Thirteen species of Bradyrhizobium were thereby distinguished in the prokaryotic microbiome, with strain-level identification of the two Bradyrhizobium species most reported to nodulate soybean. Moreover, the applied methodology delineated structural and compositional dynamism in response to experimental parameters (i.e., growth stage, cultivar, and biostimulant application), unveiled a saprotroph-rich core microbiome, provided empirical evidence for host selection of mutualistic taxa, and identified key microbial co-occurrence network members likely associated with edaphic and agronomic properties.

CONCLUSIONS

This study is the first to combine synthetic long-read technology and avidity sequencing to profile both eukaryotic and prokaryotic fractions of a plant-associated microbiome. Findings herein provide an unparalleled taxonomic resolution of the soybean rhizosphere microbiota and represent significant biological and technological advancements in crop microbiome research.

Short-read full-length 16S rRNA amplicon sequencing for characterisation of the respiratory bacteriome of captive and free-ranging African elephants (Loxodonta africana)

Hypervariable region sequencing of the 16S ribosomal RNA (rRNA) gene plays a critical role in microbial ecology by offering insights into bacterial communities within specific niches. While providing valuable genus-level information, its reliance on data from targeted genetic regions limits its overall utility. Recent advances in sequencing technologies have enabled characterisation of the full-length 16S rRNA gene, enhancing species-level classification. Although current short-read platforms are cost-effective and precise, they lack full-length 16S rRNA amplicon sequencing capability. This study aimed to evaluate the feasibility of a modified 150 bp paired-end full-length 16S rRNA amplicon short-read sequencing technique on the Illumina iSeq 100 and 16S rRNA amplicon assembly workflow by utilising a standard mock microbial community and subsequently performing exploratory characterisation of captive (zoo) and free-ranging African elephant (Loxodonta africana) respiratory microbiota. Our findings demonstrate that, despite generating assembled amplicons averaging 869 bp in length, this sequencing technique provides taxonomic assignments consistent with the theoretical composition of the mock community and respiratory microbiota of other mammals. Tentative bacterial signatures, potentially representing distinct respiratory tract compartments (trunk and lower respiratory tract) were visually identified, necessitating further investigation to gain deeper insights into their implication for elephant physiology and health.

Using short-read 16S rRNA sequencing of multiple variable regions to generate high-quality results to a species level

Introduction

Short-read amplicon sequencing studies have typically focused on 1-2 variable regions of the 16S rRNA gene. Species-level resolution is limited in these studies, as each variable region enables the characterisation of a different subsection of the microbiome. Although long-read sequencing techniques take advantage of all 9 variable regions by sequencing the entire 16S rRNA gene, they are substantially more expensive. This work assessed the feasibility of accurate species-level resolution and reproducibility using a relatively new sequencing kit and bioinformatics pipeline developed for short-read sequencing of multiple variable regions of the 16S rRNA gene. In addition, we evaluated the potential impact of different sample collection methods on our outcomes.

Methods

Using xGen™ 16S Amplicon Panel v2 kits, sequencing of all 9 variable regions of the 16S rRNA gene was carried out on an Illumina MiSeq platform. Mock cells and mock DNA for 8 bacterial species were included as extraction and sequencing controls respectively. Within-run and between-run replicate samples, and pairs of stool and rectal swabs collected at 0-5 weeks from the same participants, were incorporated. Observed relative abundances of each species were compared to theoretical abundances provided by ZymoBIOMICS. Paired Wilcoxon rank sum tests and distance-based intraclass correlation coefficients were used to statistically compare alpha and beta diversity measures, respectively, for pairs of replicates and stool/rectal swab sample pairs.

Results

Using multiple variable regions of the 16S ribosomal Ribonucleic Acid (rRNA) gene, we found that we could accurately identify taxa to a species level and obtain highly reproducible results at a species level. Yet, the microbial profiles of stool and rectal swab sample pairs differed substantially despite being collected concurrently from the same infants.

Conclusion

This protocol provides an effective means for studying infant gut microbial samples at a species level. However, sample collection approaches need to be accounted for in any downstream analysis.

Co-occurrence pattern of ARGs and N-functional genes in the aerobic composting system with initial elevated temperature

Animal manure is known to harbor antibiotic resistance genes (ARGs). Aerobic composting is a prevalent cost-effective and sustainable method to treat animal waste. However, the effect of initially elevated temperature on antibiotic resistome during the composting process is unclear. In this study composting was subjected to initial external heating (EHC) for a period of 5 days compared to conventional composting (CC). After composting ARGs abundance was significantly reduced by 2.43 log in EHC and 1.95 log in CC. Mobile genetic elements (MGEs) also exhibited a reduction of 1.95 log in EHC and 1.49 log in CC. However, during the cooling phase, the genes resisting macrolide lincosamide and streptogramin B (MLSB) rebounded by 0.04 log in CC. The potential human pathogenic bacteria Pseudomonas (41.5-61.5%) and Actinobacteria (98.4-98.8%) were significantly reduced in both treatments and the bulk of targeted antibiotics were eliminated by 80.74% in EHC and 68.98% in CC. ARGs and N-functional genes (NFGs), mainly denitrification genes, were carried by the same microbial species, such as Corynebacterium sp. and Bacillus sp., of the dominant phylum. Redundancy analysis (RDA) revealed that CC microbial communities played a key role in the enrichment of ARGs while in EHC the variation of ARGs was attributed to the composting temperature. The number of high-risk ARGs was also lower in EHC (4) compared with CC (6) on day 30. These results provide insight into the effects of an initially enhanced temperature on ARGs removal and the relationship between ARGs and NFGs during the composting process.

Matrix-associated microbial communities in a nitrogen-removing on-site wastewater treatment system are largely structured by niche processes

On-site wastewater treatment systems (OWTSs) can be designed to promote microbial communities with naturally occurring metabolic functions desirable to wastewater treatment. Among such OWTSs are nitrogen-removing biofilters (NRBs), comprising a sand layer overlying a sand-lignocellulose (sand-lc) layer and intended to promote sequential nitrification and denitrification. The design of NRBs is based on the hypothesis that niche processes like environmental selection strongly structure the microbial communities, which predicts that immigrating wastewater communities and matrix-associated communities will be distinct and that the matrix communities in the two layers will be distinct. We characterized NRB microbial communities by 16S ribosomal RNA amplicon sequencing. Selection of the matrix-associated communities was indicated by clear differences from the immigrating community. For matrix-associated communities, alpha and beta diversity differed between the matrix layers, as did the relative abundances of many functional groups and genera. Functional groups with strict metabolisms were nearly exclusively detected in either the sand (ammonia and nitrite oxidizers) or sand-lc layer (methanogens), consistent with the niche hypothesis. Contrary to expectations, denitrifiers as a functional group were not present at greater relative abundance in the sand-lc than sand matrix because of a portfolio effect: some denitrifying genera were more abundant in the sand layer, whereas others were more abundant in the sand-lc layer. This study reveals niche processes acting at different levels of community organization for different biogeochemical functions, a crucial consideration in designing effective and reliable OWTSs to mitigate nitrogen pollution.

Fine-scale evaluation of two standard 16S rRNA gene amplicon primer pairs for analysis of total prokaryotes and archaeal nitrifiers in differently managed soils

The advance of high-throughput molecular biology tools allows in-depth profiling of microbial communities in soils, which possess a high diversity of prokaryotic microorganisms. Amplicon-based sequencing of 16S rRNA genes is the most common approach to studying the richness and composition of soil prokaryotes. To reliably detect different taxonomic lineages of microorganisms in a single soil sample, an adequate pipeline including DNA isolation, primer selection, PCR amplification, library preparation, DNA sequencing, and bioinformatic post-processing is required. Besides DNA sequencing quality and depth, the selection of PCR primers and PCR amplification reactions arguably have the largest influence on the results. This study tested the performance and potential bias of two primer pairs, i.e., 515F (Parada)-806R (Apprill) and 515F (Parada)-926R (Quince) in the standard pipelines of 16S rRNA gene Illumina amplicon sequencing protocol developed by the Earth Microbiome Project (EMP), against shotgun metagenome-based 16S rRNA gene reads. The evaluation was conducted using five differently managed soils. We observed a higher richness of soil total prokaryotes by using reverse primer 806R compared to 926R, contradicting to in silico evaluation results. Both primer pairs revealed various degrees of taxon-specific bias compared to metagenome-derived 16S rRNA gene reads. Nonetheless, we found consistent patterns of microbial community variation associated with different land uses, irrespective of primers used. Total microbial communities, as well as ammonia oxidizing archaea (AOA), the predominant ammonia oxidizers in these soils, shifted along with increased soil pH due to agricultural management. In the unmanaged low pH plot abundance of AOA was dominated by the acid-tolerant NS-Gamma clade, whereas limed agricultural plots were dominated by neutral-alkaliphilic NS-Delta/NS-Alpha clades. This study stresses how primer selection influences community composition and highlights the importance of primer selection for comparative and integrative studies, and that conclusions must be drawn with caution if data from different sequencing pipelines are to be compared.

Multi-locus evaluation of gastrointestinal bacterial communities from Zalophus californianus pups in the Gulf of California, México

Background

The gastrointestinal (GI) bacterial communities of sea lions described to date have occasionally revealed large intraspecific variability, which may originate from several factors including different methodological approaches. Indeed, GI bacterial community surveys commonly rely on the use of a single hypervariable region (HR) of 16S rRNA, which may result in misleading structural interpretations and limit comparisons among studies. Here, we considered a multi-locus analysis by targeting six HRs of 16S rRNA with the aims of (i) comprehensively assessing the GI bacterial consortium in rectal samples from Zalophus californianus pups and (ii) elucidating structural variations among the tested HRs. In addition, we evaluated which HRs may be most suitable for identifying intrinsic, structurally related microbiome characteristics, such as geographic variations or functional capabilities.

Methods

We employed a Short MUltiple Regions Framework (SMURF) approach using the Ion 16S™ Metagenomic Kit. This kit provides different proprietary primers designed to target six HRs of the 16S rRNA gene. To date, the only analytical pipeline available for this kit is the Ion Reporter™ Software of Thermo Fisher Scientific. Therefore, we propose an in-house pipeline to use with open-access tools, such as QIIME2 and PICRUSt 2, in downstream bioinformatic analyses.

Results

As hypothesized, distinctive bacterial community profiles were observed for each analyzed HR. A higher number of bacterial taxa were detected with the V3 and V6-V7 regions. Conversely, the V8 and V9 regions were less informative, as we detected a lower number of taxa. The synergistic information of these HRs suggests that the GI microbiota of Zalophus californianus pups is predominated by five bacterial phyla: Proteobacteria (~50%), Bacteroidetes (~20%), Firmicutes (~18%), Fusobacteria (~7%), and Epsilonbacteraeota (~4%). Notably, our results differ at times from previously reported abundance profiles, which may promote re-evaluations of the GI bacterial compositions in sea lions and other pinniped species that have been reported to date. Moreover, consistent geographic differences were observed only with the V3, V4, and V6-V7 regions. In addition, these HRs also presented higher numbers of predicted molecular pathways, although no significant functional changes were apparent. Together, our results suggests that multi-locus analysis should be encouraged in GI microbial surveys, as single-locus approaches may result in misleading structural results that hamper the identification of structurally related microbiome features.

Seminal Microbiota of Idiopathic Infertile Patients and Its Relationship With Sperm DNA Integrity

The development of new biomarkers for human male infertility is crucial to improve the diagnosis and the prognosis of this disease. Recently, seminal microbiota was shown to be related to sperm quality parameters, suggesting an effect in human fertility and postulating it as a biomarker candidate. However, its relationship to sperm DNA integrity has not been studied yet. The aim of the present study is to characterize the seminal microbiota of a western Mediterranean population and to evaluate its relationship to sperm chromatin integrity parameters, and oxidative stress. For that purpose, 14 samples from sperm donors and 42 samples from infertile idiopathic patients were obtained and were analyzed to assess the composition of the microbiota through full-length 16S rRNA gene sequencing (Illumina MiSeq platform). Microbial diversity and relative abundances were compared to classic sperm quality parameters (macroscopic semen parameters, motility, morphology and concentration), chromatin integrity (global DNA damage, double-stranded DNA breaks and DNA protamination status) and oxidative stress levels (oxidation-reduction potential). The seminal microbiota observed of these samples belonged to the phyla Firmicutes, Proteobacteria, Actinobacteria and Bacteroidetes. The most abundant genera were Finegoldia, Peptoniphilus, Anaerococcus, Campylobacter, Streptococcus, Staphylococcus, Moraxella, Prevotella, Ezakiella, Corynebacterium and Lactobacillus. To our knowledge, this is the first detection of Ezakiella genus in seminal samples. Two clusters of microbial profiles were built based on a clustering analysis, and specific genera were found with different frequencies in relation to seminal quality defects. The abundances of several bacteria negatively correlate with the sperm global DNA fragmentation, most notably Moraxella, Brevundimonas and Flavobacterium. The latter two were also associated with higher sperm motility and Brevundimonas additionally with lower oxidative-reduction potential. Actinomycetaceae, Ralstonia and Paenibacillus correlated with reduced chromatin protamination status and increased double-stranded DNA fragmentation. These effects on DNA integrity coincide in many cases with the metabolism or enzymatic activities of these genera. Significant differences between fertile and infertile men were found in the relative presence of the Propionibacteriaceae family and the Cutibacterium, Rhodopseudomonas and Oligotropha genera, which supports its possible involvement in male fertility. Our findings sustain the hypothesis that the seminal microbiome has an effect on male fertility.

Effects of Waterlogging on Soybean Rhizosphere Bacterial Community Using V4, LoopSeq, and PacBio 16S rRNA Sequence

Waterlogging causes a significant reduction in soil oxygen levels, which in turn negatively affects soil nutrient use efficiency and crop yields. Rhizosphere microbes can help plants to better use nutrients and thus better adapt to this stress, while it is not clear how the plant-associated microbes respond to waterlogging stress. There are also few reports on whether this response is influenced by different sequencing methods and by different soils. In this study, using partial 16S rRNA sequencing targeting the V4 region and two full-length 16S rRNA sequencing approaches targeting the V1 to V9 regions, the effects of waterlogging on soybean rhizosphere bacterial structure in two types of soil were examined. Our results showed that, compared with the partial 16S sequencing, full-length sequencing, both LoopSeq and Pacific Bioscience (PacBio) 16S sequencing, had a higher resolution. On both types of soil, all the sequencing methods showed that waterlogging significantly affected the bacterial community structure of the soybean rhizosphere and increased the relative abundance of Geobacter. Furthermore, modular analysis of the cooccurrence network showed that waterlogging increased the relative abundance of some microorganisms related to nitrogen cycling when using V4 sequencing and increased the microorganisms related to phosphorus cycling when using LoopSeq and PacBio 16S sequencing methods. Core microorganism analysis further revealed that the enriched members of different species might play a central role in maintaining the stability of bacterial community structure and ecological functions. Together, our study explored the role of microorganisms enriched at the rhizosphere under waterlogging in assisting soybeans to resist stress. Furthermore, compared to partial and PacBio 16S sequencing, LoopSeq offers improved accuracy and reduced sequencing prices, respectively, and enables accurate species-level and strain identification from complex environmental microbiome samples. IMPORTANCE Soybeans are important oil-bearing crops, and waterlogging has caused substantial decreases in soybean production all over the world. The microbes associated with the host have shown the ability to promote plant growth, nutrient absorption, and abiotic resistance. High-throughput sequencing of partial 16S rRNA is the most commonly used method to analyze the microbial community. However, partial sequencing cannot provide correct classification information below the genus level, which greatly limits our research on microbial ecology. In this study, the effects of waterlogging on soybean rhizosphere microbial structure in two soil types were explored using partial 16S rRNA and full-length 16S gene sequencing by LoopSeq and Pacific Bioscience (PacBio). The results showed that full-length sequencing had higher classification resolution than partial sequencing. Three sequencing methods all indicated that rhizosphere bacterial community structure was significantly impacted by waterlogging, and the relative abundance of Geobacter was increased in the rhizosphere in both soil types after suffering waterlogging. Moreover, the core microorganisms obtained by different sequencing methods all contain species related to nitrogen cycling. Together, our study not only explored the role of microorganisms enriched at the rhizosphere level under waterlogging in assisting soybean to resist stress but also showed that LoopSeq sequencing is a less expensive and more convenient method for full-length sequencing by comparing different sequencing methods.

Human- and infrastructure-associated bacteria in greywater

Greywater, the wastewater from sinks, showers and laundry, is an understudied environment for bacterial communities. Most greywater studies focus on quantifying pathogens, often via proxies used in other wastewater, like faecal indicator bacteria; there is a need to identify more greywater-appropriate surrogates, like Staphylococcus sp. Sequencing-based studies have revealed distinct communities in different types of greywater as well as in different parts of greywater infrastructure, including biofilms on pipes, holding tanks and filtration systems. The use of metagenomic sequencing provides high resolution on both the taxa and genes present, which may be of interest in cases like identifying pathogens and surrogates relevant to different matrices, monitoring antibiotic resistance genes and understanding metabolic processes occurring in the system. Here, we review what is known about bacterial communities in different types of greywater and its infrastructure. We suggest that wider adoption of environmental sequencing in greywater research is important because it can describe the entire bacterial community along with its metabolic capabilities, including pathways for removal of nutrients and organic materials. We briefly describe a metagenomic dataset comparing different types of greywater samples in a college dormitory building to highlight the type of questions these methods can address. Metagenomic sequencing can help further the understanding of greywater treatment for reuse because it allows for identification of new pathogens or genes of concern.

Species-Level Resolution of Female Bladder Microbiota from 16S rRNA Amplicon Sequencing

The human bladder contains bacteria, even in the absence of infection. Interest in studying these bacteria and their association with bladder conditions is increasing. However, the chosen experimental method can limit the resolution of the taxonomy that can be assigned to the bacteria found in the bladder. 16S rRNA amplicon sequencing is commonly used to identify bacteria in urinary specimens, but it is typically restricted to genus-level identification. Our primary aim here was to determine if accurate species-level identification of bladder bacteria is possible using 16S rRNA amplicon sequencing. We evaluated the ability of different classification schemes, each consisting of combinations of a reference database, a 16S rRNA gene variable region, and a taxonomic classification algorithm to correctly classify bladder bacteria. We show that species-level identification is possible and that the reference database chosen is the most important component, followed by the 16S variable region sequenced. IMPORTANCE Accurate species-level identification from culture-independent techniques is of importance for microbial niches that are less well characterized, such as that of the bladder. 16S rRNA amplicon sequencing, a common culture-independent way to identify bacteria, is often critiqued for lacking species-level resolution. Here, we extensively evaluate classification schemes for species-level bacterial annotation of 16S amplicon data from bladder bacteria. Our results show that the proper choice of taxonomic database and variable region of the 16S rRNA gene sequence makes species level identification possible. We also show that this improvement can be achieved through the more careful application of existing methods and resources. Species-level information may deepen our understanding of associations between bacteria in the bladder and bladder conditions such as lower urinary tract symptoms and urinary tract infections.

Genome-Based Targeted Sequencing as a Reproducible Microbial Community Profiling Assay

Current sequencing-based methods for profiling microbial communities rely on marker gene (e.g., 16S rRNA) or metagenome shotgun sequencing (mWGS) analysis. We present an approach based on a single-primer extension reaction using a highly multiplexed oligonucleotide probe pool. This approach, termed MA-GenTA (microbial abundances from genome tagged analysis), enables quantitative, straightforward, cost-effective microbiome profiling that combines desirable features of both 16S rRNA and mWGS strategies. The use of multiple probes per target genome and rigorous probe design criteria enabled robust determination of relative abundance. To test the utility of the MA-GenTA assay, probes were designed for 830 genome sequences representing bacteria present in mouse stool specimens. Comparison of the MA-GenTA data with mWGS data demonstrated excellent correlation down to 0.01% relative abundance and a similar number of organisms detected per sample. Despite the incompleteness of the reference database, nonmetric multidimensional scaling (NMDS) clustering based on the Bray-Curtis dissimilarity metric of sample groups was consistent between MA-GenTA, mWGS, and 16S rRNA data sets. MA-GenTA represents a potentially useful new method for microbiome community profiling based on reference genomes.

IMPORTANCE New methods for profiling the microbial communities can create new approaches to understanding the composition and function of those communities. In this study, we combined bacterial genome-specific probe design with a highly multiplexed single primer extension reaction as a new method to profile microbial communities, using stool from various mouse strains as a test case. This method, termed MA-GenTA, was benchmarked against 16S rRNA gene sequencing and metagenome sequencing methods and delivered similar relative abundance and clustering data. Since the probes were generated from reference genomes, MA-GenTA was also able to provide functional pathway data for the stool microbiome in the assayed samples. The method is more informative than 16S rRNA analysis while being less costly than metagenome shotgun sequencing.

LightCUD: a program for diagnosing IBD based on human gut microbiome data

BACKGROUND

The diagnosis of inflammatory bowel disease (IBD) and discrimination between the types of IBD are clinically important. IBD is associated with marked changes in the intestinal microbiota. Advances in next-generation sequencing (NGS) technology and the improved hospital bioinformatics analysis ability motivated us to develop a diagnostic method based on the gut microbiome.

RESULTS

Using a set of whole-genome sequencing (WGS) data from 349 human gut microbiota samples with two types of IBD and healthy controls, we assembled and aligned WGS short reads to obtain feature profiles of strains and genera. The genus and strain profiles were used for the 16S-based and WGS-based diagnostic modules construction respectively. We designed a novel feature selection procedure to select those case-specific features. With these features, we built discrimination models using different machine learning algorithms. The machine learning algorithm LightGBM outperformed other algorithms in this study and thus was chosen as the core algorithm. Specially, we identified two small sets of biomarkers (strains) separately for the WGS-based health vs IBD module and ulcerative colitis vs Crohn's disease module, which contributed to the optimization of model performance during pre-training. We released LightCUD as an IBD diagnostic program built with LightGBM. The high performance has been validated through five-fold cross-validation and using an independent test data set. LightCUD was implemented in Python and packaged free for installation with customized databases. With WGS data or 16S rRNA sequencing data of gut microbiome samples as the input, LightCUD can discriminate IBD from healthy controls with high accuracy and further identify the specific type of IBD. The executable program LightCUD was released in open source with instructions at the webpage http://cqb.pku.edu.cn/ZhuLab/LightCUD/ . The identified strain biomarkers could be used to study the critical factors for disease development and recommend treatments regarding changes in the gut microbial community.

CONCLUSIONS

As the first released human gut microbiome-based IBD diagnostic tool, LightCUD demonstrates a high-performance for both WGS and 16S sequencing data. The strains that either identify healthy controls from IBD patients or distinguish the specific type of IBD are expected to be clinically important to serve as biomarkers.

An high-throughput sequencing approach to the preliminary analysis of bacterial communities associated with changes in amino acid nitrogen, organic acid and reducing sugar contents during soy sauce fermentation

Soy sauce is a traditional fermented soy food for enhancing the umami taste in Asian cuisines. In this study, 16S rRNA gene throughput sequencing analysis showed the bacterial communities and the changes in soy sauce during fermentation. Weissella, Bacillus and Lactococcus were the most abundant at genus level. The uncultured bacterium Weissella and Lactococcus had relatively high abundance at species level. Alpha diversity analysis indicated the bacterial community diversity increased at fermentation initiation, while decreased as fermentation progressed. Based on beta-diversity analysis, four clusters including cluster I (time point A-F), cluster II (G,H), cluster III (I,J) and cluster IV(K) were distinctly separated, indicating the fermentation time significantly affected bacterial community diversity. Also, close associations were found between the bacterial communities in soy sauce and its amino acid nitrogen, organic acid and reducing sugar contents during fermentation. Therefore, it will provide important information for optimization of the soy sauce production process.

Analysis of molecular diversity within single cyanobacterial colonies from environmental samples

Attached or floating macroscopic cyanobacteria can be found in shallow waters and can be easily hand-collected, but their identification is often challenging due to their high morphological variability. In addition, many members of environmental samples lose their morphological adaptations under controlled conditions, making the integration of analyses of field populations and derived isolated cultures necessary in order to evaluate phenotypic plasticity for identification purposes. Therefore, in this study, twenty-nine macroscopic field samples were analyzed by Illumina sequencing and parallel optical microscopy. Some colonies showed the typical morphological characteristics of Rivularia biasolettiana, and others showed those of Rivularia haematites. However, other Rivularia-like colonies showed ambiguous morphologies, and some of them showed the phenotypic features of the new genus Cyanomargarita, which is virtually indistinguishable from Rivularia in the field. In all of the colonies, phylotype composition was highly heterogeneous, with abundances varying depending on the analyzed sample. Some colonies were dominated (97-99%) by a single phylotype, while in others, the percentage of the dominant phylotype decreased to approximately 50-60%. Surprisingly, the same dominant phylotype was found in R. biasolettiana and R. haematites colonies. The relationships between environmental and/or biological factors and morphological variability in these colonies are discussed.

Responses of microbial structures, functions, metabolic pathways and community interactions to different C/N ratios in aerobic nitrification

The responses of microbial structures, functional profiles and metabolic pathways during nitrification to four C/N ratios (0, 5, 10 and 15) were investigated in four parallel SBRs denoted as S₀, S₅, S₁₀, S₁₅. Results indicated that microbial diversities were affected by C/N ratios, while the same dominant taxa were observed, mainly including Proteobacteria, Betaproteobacteria, Rhodocyclales, Rhodocyclaceae, Zoogloea. The unique biomarkers were identified in each sludge sample through LEfSe analysis. Functional genera/enzymes responsible for removing organics and nitrogen coexisted in four SBRs at different abundances, except for that ammonia oxidizing bacteria (AOB) Nitrosomonas (0.33%-0.66%) and ammonia monooxygenase (amo) (9.4 × 10^-7-2.8 × 10^-6) were only detected in S₀. Moreover, PICRUSt analysis indicated similar overall patterns of metabolic pathways in four sludge samples. The network analysis revealed that total nitrogen removal positively correlated with hcp (Spearman's ρ of 0.853), and ammonia oxidizing rate was associated with amo (Spearman's ρ of 0.096).

Spider phylosymbiosis: divergence of widow spider species and their tissues' microbiomes

BACKGROUND

Microbiomes can have profound impacts on host biology and evolution, but to date, remain vastly understudied in spiders despite their unique and diverse predatory adaptations. This study evaluates closely related species of spiders and their host-microbe relationships in the context of phylosymbiosis, an eco-evolutionary pattern where the microbial community profile parallels the phylogeny of closely related host species. Using 16S rRNA gene amplicon sequencing, we characterized the microbiomes of five species with known phylogenetic relationships from the family Theridiidae, including multiple closely related widow spiders (L. hesperus, L. mactans, L. geometricus, S. grossa, and P. tepidariorum).

RESULTS

We compared whole animal and tissue-specific microbiomes (cephalothorax, fat bodies, venom glands, silk glands, and ovary) in the five species to better understand the relationship between spiders and their microbial symbionts. This showed a strong congruence of the microbiome beta-diversity of the whole spiders, cephalothorax, venom glands, and silk glands when compared to their host phylogeny. Our results support phylosymbiosis in these species and across their specialized tissues. The ovary tissue microbial dendrograms also parallel the widow phylogeny, suggesting vertical transfer of species-specific bacterial symbionts. By cross-validating with RNA sequencing data obtained from the venom glands, silk glands and ovaries of L. hesperus, L. geometricus, S. grossa, and P. tepidariorum we confirmed that several microbial symbionts of interest are viably active in the host.

CONCLUSION

Together these results provide evidence that supports the importance of host-microbe interactions and the significant role microbial communities may play in the evolution and adaptation of their hosts.

Evaluation of PCR conditions for characterizing bacterial communities with full-length 16S rRNA genes using a portable nanopore sequencer

MinION (Oxford Nanopore Technologies), a portable nanopore sequencer, was introduced in 2014 as a new DNA sequencing technology. MinION is now widely used because of its low initial start-up costs relative to existing DNA sequencers, good portability, easy-handling, real-time analysis and long-read output. However, differences in the experimental conditions used for 16S rRNA-based PCR can bias bacterial community assessments in samples. Therefore, basic knowledge about reliable experimental conditions is needed to ensure the appropriate use of this technology. Our study concerns the reliability of techniques for obtaining accurate and quantitative full-length 16S rRNA amplicon sequencing data for bacterial community structure assessment using MinION. We compared five PCR conditions using three independent mock microbial community standard DNAs and established appropriate, standardized, better PCR conditions among the trials. We then sequenced two mock communities and six environmental samples using Illumina MiSeq for comparison. Modifying the PCR conditions improved the sequencing quality; the optimized conditions were 35 cycles of 95 °C for 1 min, 60 °C for 1 min and 68 °C for 3 min. Our results provide important information for researchers to determine bacterial community using MinION accurately.

Distinct Skin Microbiota Imbalance and Responses to Clinical Treatment in Children With Atopic Dermatitis

Background: Atopic dermatitis (AD) is a common cutaneous disease, associated with imbalances in the skin microbiota. Objective: To explore the characteristics of the cutaneous microbiota and its dynamic changes during clinical treatment. Methods: Cutaneous swab samples were collected from 51 AD patients before treatment, and 40 AD patients remained after a 2-week treatment with mometasone and mupirocin. Results: AD patients exhibited significant enrichments of Prevotella and Desulfovibrio as well as obvious reductions of Corynebacterium, Streptococcus and Parabacteroides. Based on the proportion of Staphylococcus aureus, the AD patients were further classified into S. aureus-predominant group (AD.S) and S. aureus-non-dominant (AD.ND) group. The AD.S group exhibited lower skin microbial diversity and higher atopic dermatitis (SCORAD) index. In the AD.S group, the cutaneous microbial diversity significantly increased from 2.9 ± 0.8 to 3.7 ± 1.0, while the relative abundance of S. aureus decreased from 42.5 ± 20.7 to 10.3 ± 28.4 after treatment. In contrast, no significant skin microbiota changes were detected in the AD.ND group. Conclusions: AD patients with predominant S. aureus had higher disease severity and lower microbiota diversity compared to patients in the AD.ND group. Mometasone and mupirocin therapy had significant effects on skin microbiota in AD.S patients, but had a paradoxical response in the AD.ND patients.

Simultaneous Removal of Dissolved Methane and Nitrogen from Synthetic Mainstream Anaerobic Effluent

Anaerobic technologies have been proposed as a promising solution to enhance bioenergy recovery and to transform a wastewater treatment plant (WWTP) from an energy consumer to an energy exporter. However, 20-60% of the methane produced remains dissolved in the anaerobically treated effluent, which is a potent greenhouse gas and is easily stripped out in the aeration tank. This study aims to develop a solution using dissolved methane to support denitrification, thus simultaneously enhancing nitrogen removal and achieving beneficial use of dissolved methane. By coupling anaerobic ammonium oxidation (anammox) with nitrite/nitrate-dependent anaerobic methane oxidation (n-DAMO), up to 85% of dissolved methane and more than 99% of nitrogen were removed in parallel in a biofilm system. Mass balance was conducted during both long-term operation and short-term batch tests, which indicated that n-DAMO bacteria and n-DAMO archaea indeed contributed jointly to the methane removal. The 16S rRNA gene amplicon sequencing further showed the co-presence of n-DAMO bacteria and n-DAMO archaea, while anammox bacteria were detected with a low relative abundance. This proposed technology can potentially be applied to reduce the carbon footprint and to save the organic carbon consumption in WWTPs.

Influence of 16S rRNA variable region on perceived diversity of marine microbial communities of the Northern North Atlantic

Marine microbes play essential roles in global energy and nutrient cycles. A primary method of determining their diversity and distribution is through sequencing of 16S ribosomal RNA genes from environmental samples. However, the perceived community composition may vary significantly based on differences in methodology, including choice of 16S variable region(s). This study investigated the influence of 16S variable region selection (V4-V5 or V6-V8) on perceived community composition and diversity for bacteria, Archaea and chloroplasts by tag-Illumina sequencing. We used 24 samples from the photic zone of the Scotian Shelf, northwest Atlantic, collected during a spring phytoplankton bloom. Taxonomic assignment and community composition varied greatly depending on the choice of variable regions while observed patterns of beta diversity were reproducible between variable regions. V4-V5 was considered the preferred variable region for future studies based on its superior recognition of Archaea, which has received little attention in bloom dynamics. The V6-V8 region captured more of the bacterial diversity, including the abundant SAR11 clades and, to a lesser extent, that of chloroplasts. However, the magnitude of difference between variable regions for bacteria and chloroplast was less than for Archaea.

Comparative analysis of gut bacterial communities in housefly larvae fed different diets using a high-throughput sequencing approach

Housefly larvae are a synanthropic host for various bacteria, including pathogens and commensals and an important protein source for monogastric animal feed. Many factors, such as diets, life stages, host habitats can influence microbial community structure. In this study, the diversity of bacterial communities in the gut of housefly larvae fed on different artificial diets was comprehensively characterized using high-throughput sequencing with the aim shedding light on an optimal larval diet. The results showed that the dominant bacteria belonging to Proteobacteria, Firmicutes and Bacteroidetes phyla were related to polysaccharide degradation. The comparative analysis indicated that the dominant intestinal bacteria of larvae fed on high-protein were similar to those on high-fat diet. The same was the case in larvae fed high-starch diet and wheat bran alone. In addition, the diversity of intestinal bacteria at genus level in larvae fed high-protein and high-fat diet was higher than in larvae fed the other two diets. Further analysis indicated that the increase of potential commensals and decrease of pathogens in larvae fed on high-fat diet contributed to the increase of housefly larvae immunity. It established a foundation for further research on improvement of nutrition of housefly larvae used for poultry and fish feed.

High-throughput DNA sequencing technologies for water and wastewater analysis

Conventional microbiological water monitoring uses culture-dependent techniques to screen indicator microbial species such as Escherichia coli and fecal coliforms. With high-throughput, second-generation sequencing technologies becoming less expensive, water quality monitoring programs can now leverage the massively parallel nature of second-generation sequencing technologies for batch sample processing to simultaneously obtain compositional and functional information of culturable and as yet uncultured microbial organisms. This review provides an introduction to the technical capabilities and considerations necessary for the use of second-generation sequencing technologies, specifically 16S rDNA amplicon and whole-metagenome sequencing, to investigate the composition and functional potential of microbiomes found in water and wastewater systems.

Re-assessing microbiomes in the low-biomass reproductive niche

The female reproductive tract represents a continuum between the vagina and the upper genital tract. New evidence from cultivation-independent studies suggests that the female upper genital tract is not sterile; however, the significance of this for reproductive health and disease remains to be elucidated fully. Further, diagnosis and treatment of infectious reproductive tract pathologies using cultivation-independent technologies represents a largely unchartered area of modern medical science. The challenge now is to design well-controlled experiments to account for the ease of contamination known to confound molecular-based studies of low-biomass niches, including the uterus and placenta. This will support robust assessment of the potential function of microorganisms, microbial metabolites, and cell-free bacterial DNA on reproductive function in health and disease. TWEETABLE ABSTRACT: Molecular microbial studies of low-biomass niches require stringent experimental controls to reveal causal relations in reproductive health and disease.

Applications of sequencing technology in clinical microbial infection

Infectious diseases are a type of disease caused by pathogenic microorganisms. Although the discovery of antibiotics changed the treatment of infectious diseases and reduced the mortality of bacterial infections, resistant bacterial strains have emerged. Anti‐infective therapy based on aetiological evidence is the gold standard for clinical treatment, but the time lag and low positive culture rate of traditional methods of pathogen diagnosis leads to relative difficulty in obtaining the evidence of pathogens. Compared with traditional methods of pathogenic diagnosis, next‐generation and third‐generation sequencing technologies have many advantages in the detection of pathogenic microorganisms. In this review, we mainly introduce recent progress in research on pathogenic diagnostic technology and the applications of sequencing technology in the diagnosis of pathogenic microorganisms. This review provides new insights into the application of sequencing technology in the clinical diagnosis of microorganisms.

On resolving ambiguities in microbial community analysis of partial nitritation anammox reactors

PCR-based methods have caused a surge for integration of eco-physiological approaches into research on partial nitritation anammox (PNA). However, a lack of rigorous standards for molecular analyses resulted in widespread data misinterpretation and consequently lack of consensus. Data consistency and accuracy strongly depend on the primer selection and data interpretation. An in-silico evaluation of 16S rRNA gene eubacterial primers used in PNA studies from the last ten years unraveled the difficulty of comparing ecological data from different studies due to a variation in the coverage of these primers. Our 16S amplicon sequencing approach, which includes parallel sequencing of six 16S rRNA hypervariable regions, showed that there is no perfect hypervariable region for PNA microbial communities. Using qPCR analysis, we emphasize the significance of primer choice for quantification and caution with data interpretation. We also provide a framework for PCR based analyses that will improve and assist to objectively interpret and compare such results.

Microbial residence time is a controlling parameter of the taxonomic composition and functional profile of microbial communities

A remaining challenge within microbial ecology is to understand the determinants of richness and diversity observed in environmental microbial communities. In a range of systems, including activated sludge bioreactors, the microbial residence time (MRT) has been previously shown to shape the microbial community composition. However, the physiological and ecological mechanisms driving this influence have remained unclear. Here, this relationship is explored by analyzing an activated sludge system fed with municipal wastewater. Using a model designed in this study based on Monod-growth kinetics, longer MRTs were shown to increase the range of growth parameters that enable persistence, resulting in increased richness and diversity in the modeled community. In laboratory experiments, six sequencing batch reactors treating domestic wastewater were operated in parallel at MRTs between 1 and 15 days. The communities were characterized using both 16S ribosomal RNA and non-target messenger RNA sequencing (metatranscriptomic analysis), and model-predicted monotonic increases in richness were confirmed in both profiles. Accordingly, taxonomic Shannon diversity also increased with MRT. In contrast, the diversity in enzyme class annotations resulting from the metatranscriptomic analysis displayed a non-monotonic trend over the MRT gradient. Disproportionately high abundances of transcripts encoding for rarer enzymes occur at longer MRTs and lead to the disconnect between taxonomic and functional diversity profiles.

Contributions and Challenges of High Throughput qPCR for Determining Antimicrobial Resistance in the Environment: A Critical Review

Expansion in whole genome sequencing and subsequent increase in antibiotic resistance targets have paved the way of high throughput qPCR (HT-qPCR) for analyzing hundreds of antimicrobial resistance genes (ARGs) in a single run. A meta-analysis of 51 selected studies is performed to evaluate ARGs abundance trends over the last 7 years. WaferGen^TM SmartChip is found to be the most widely used HT-qPCR platform among others for evaluating ARGs. Up till now around 1000 environmental samples (excluding biological replicates) from different parts of the world have been analyzed on HT-qPCR. Calculated detection frequency and normalized ARGs abundance (ARGs/16S rRNA gene) reported in gut microbiome studies have shown a trend of low ARGs as compared to other environmental matrices. Disparities in the HT-qPCR data analysis which are causing difficulties to researchers in precise interpretation of results have been highlighted and a possible way forward for resolving them is also suggested. The potential of other amplification technologies and point of care or field deployable devices for analyzing ARGs have also been discussed in the review. Our review has focused on updated information regarding the role, current status and future perspectives of HT-qPCR in the field of antimicrobial resistance.

Characterization of microbes and denitrifiers attached to two species of floating plants in the wetlands of Lake Taihu

Biofilms are often observed at the solid-water interface. The leaves of many floating macrophytes have characteristics of both terrestrial plants and submerged macrophytes, because, in general, their upper and lower surfaces are exposed to air and water, respectively. However, little is known about the biofilms attached to floating plants. We investigated biofilms attached to the leaves, stems and roots of the floating plants Nymphoides peltata (in summer and winter) and Trapa natans (in summer) in the Gonghu Bay of Lake Taihu. Bacteria and algae were major components of the biofilm on the leaves of the two species of plants. In addition, 454 pyrosequencing analysis of bacterial 16S rRNA genes revealed that Proteobacteria was the dominant phylum, followed by Bacteroidetes, Firmicutes, Chloroflexi, Acidobacteria, and Verrucomicrobia. Cluster analysis showed that bacterial communities from the same plant source were clustered into the same group. A total of 677 genera were detected, and 47 genera were shared by all samples. Nitrifiers, including Nitrosomonas, Nitrosococcus and Nitrospira were detected in this study. Seven denitrifying genes (napA, napG, nirS, nirK, cnorB, qnorB and nosZ) were used to detect the abundance of denitrifiers. Genes nirK, nirS cnorB and nosZ were the four most abundant genes in all samples. Our results demonstrated that cultivation of floating plants in water column could enlarge the area for biofilm growth, and biofilms might play an important role in denitrification in eutrophic water.

Microbiome differences between river-dwelling and cave-adapted populations of the fish Astyanax mexicanus (De Filippi, 1853)

Symbiotic relationships between host and microbiome can play a major role in local adaptation. Previous studies with freshwater organisms have shown that microbiome performs numerous important biochemical functions for the host, playing a key role in metabolism, physiology or health. Experimental studies in fish groups have found an effect of enzymatic activity of gut microbiota on a variety of metabolic processes. The goal of this study was to compare stomach microbiome from cave and surface Astyanax mexicanus, in order to evaluate the potential response of microbiota to contrasting environmental conditions and physiological adaptations of the host. Stomach microbiota was obtained from three different populations: Pachón cave, and two surface rivers (Rascón and Micos rivers). The stomach microbiome was analyzed using the Ion 16S Metagenomic kit considering seven variable regions: V2, V3, V4, V6-7, V8 and V9. A high diversity was observed across samples, including 16 phyla, 120 families and 178 genera. Gammaproteobacteria, Firmicutes, Bacteroidetes and Betaproteobacteria were the most abundant phyla across the samples. Although the relative abundance of the core OTUs at genus level were highly contrasting among populations, we did not recover differences in stomach microbiome between contrasting habitats (cave vs. surface rivers). Rather, we observed a consistent association between β-diversity and dissolved oxygen concentration in water. Therefore, and unexpectedly, the microbiota of A. mexicanus is not linked with the contrasting conditions of the habitat considered here but is related to water parameters.

16S rRNA Gene Primer Validation for Bacterial Diversity Analysis of Vegetable Products

High-throughput sequencing of the 16S rRNA gene enhances understanding of microbial diversity from complex environmental samples. The 16S rRNA gene is currently the most important target in bacterial evolution and ecology studies, particularly for determination of phylogenetic relationships among taxa, exploration of bacterial diversity in a given environment, and quantification of the relative abundance of taxa at various levels. However, some parts of the conserved region of the bacterial 16S rRNA gene are similar to the conserved regions of plant chloroplasts and eukaryotic mitochondria. Therefore, if DNA contains a large amount of nontarget DNA, this nontarget DNA can be coamplified and consequently produce useless sequence reads. We experimentally assessed the primer pair 335f/769r and the widely used bacterial primer pair SD (S-D-Bact-0341-b-S-17/S-D-Bact-0785-a-A-21). The primer pair 335f/769r was examined for its ability to amplify bacterial DNA in plant and animal feed samples by using the single-strand confirmation polymorphism method. In our present study, these primer pairs were validated for microbial community structure analysis with complex food matrices by using next-generation sequencing. The sequencing results revealed that the primer pair 335f/769r successfully resulted in fewer chloroplast and mitochondrial sequence reads than generated by the universal primer pair SD and therefore is comparatively suitable for metagenomic analyses of complex food matrices, particularly those that are rich in plant DNA. Additionally, some taxonomic groups were missed entirely when only the SD primer pair was used.

Evaluation of different 16S rRNA gene V regions for exploring bacterial diversity in a eutrophic freshwater lake

Massive partial sequencing of 16S rRNA genes has become the predominant tool used for studying microbial ecology. However, determining which hypervariable regions and primer sets should be used for screening microbial communities requires extensive investigation if controversial results are to be avoided. Here, the performances of different variable regions of the 16S rRNA gene on bacterial diversity studies were evaluated in silico with respect to the SILVA non-redundant reference database (SILVA SSU Ref 123NR), and subsequently verified using samples from Lake Taihu in China, a eutrophic lake. We found that the bacterial community composition results were strongly impacted by the different V regions. The results show that V1-V2 and V1-V3 regions were the most reliable regions in the full-length 16S rRNA sequences, while most V3 to V6 regions (including V3, V4, V3-V4, V5, V4-V5, V6, V3-V6, V4-V6, and V5-V6) were more closely aligned with the SILVA SSU Ref 123NR database. Overall, V4 was the most prominent V region for achieving good domain specificity, higher coverage and a broader spectrum in the Bacteria domain, as confirmed by the validation experiments. S-D-Bact-0564-a-S-15/S-D-Bact-0785-b-A-18 is, therefore, a promising primer set for surveying bacterial diversity in eutrophic lakes.

Taxonomy of anaerobic digestion microbiome reveals biases associated with the applied high throughput sequencing strategies

In the past few years, many studies investigated the anaerobic digestion microbiome by means of 16S rRNA amplicon sequencing. Results obtained from these studies were compared to each other without taking into consideration the followed procedure for amplicons preparation and data analysis. This negligence was mainly due to the lack of knowledge regarding the biases influencing specific steps of the microbiome investigation process. In the present study, the main technical aspects of the 16S rRNA analysis were checked giving special attention to the approach used for high throughput sequencing. More specifically, the microbial compositions of three laboratory scale biogas reactors were analyzed before and after addition of sodium oleate by sequencing the microbiome with three different approaches: 16S rRNA amplicon sequencing, shotgun DNA and shotgun RNA. This comparative analysis revealed that, in amplicon sequencing, abundance of some taxa (Euryarchaeota and Spirochaetes) was biased by the inefficiency of universal primers to hybridize all the templates. Reliability of the results obtained was also influenced by the number of hypervariable regions under investigation. Finally, amplicon sequencing and shotgun DNA underestimated the Methanoculleus genus, probably due to the low 16S rRNA gene copy number encoded in this taxon.

Myriophyllum aquaticum Constructed Wetland Effectively Removes Nitrogen in Swine Wastewater

Removal of nitrogen (N) is a critical aspect in the functioning of constructed wetlands (CWs), and the N treatment in CWs depends largely on the presence and activity of macrophytes and microorganisms. However, the effects of plants on microorganisms responsible for N removal are poorly understood. In this study, a three-stage surface flow CW was constructed in a pilot-scale within monospecies stands of Myriophyllum aquaticum to treat swine wastewater. Steady-state conditions were achieved throughout the 600-day operating period, and a high (98.3%) average ammonia removal efficiency under a N loading rate of 9 kg ha^-1 d^-1 was observed. To determine whether this high efficiency was associated with the performance of active microbes, the abundance, structure, and interactions of microbial community were compared in the unvegetated and vegetated samples. Real-time quantitative polymerase chain reactions showed the abundances of nitrifying genes (archaeal and bacterial amoA) and denitrifying genes (nirS, nirK, and nosZ) were increased significantly by M. aquaticum in the sediments, and the strongest effects were observed for the archaeal amoA (218-fold) and nirS genes (4620-fold). High-throughput sequencing of microbial 16S rRNA gene amplicons showed that M. aquaticum greatly changed the microbial community, and ammonium oxidizers (Nitrosospira and Nitrososphaera), nitrite-oxidizing bacteria (Nitrospira), and abundant denitrifiers including Rhodoplanes, Bradyrhizobium, and Hyphomicrobium, were enriched significantly in the sediments. The results of a canonical correspondence analysis and Mantle tests indicated that M. aquaticum may shift the sediment microbial community by changing the sediment chemical properties. The enriched nitrifiers and denitrifiers were distributed widely in the vegetated sediments, showing positive ecological associations among themselves and other bacteria based on phylogenetic molecular ecological networks.

Effect of granular activated carbon on the aerobic granulation of sludge and its mechanism

The granulation of activated sludge and effect of granular activated carbon (GAC) was investigated under the alternative anaerobic and aerobic conditions. The results showed that GAC accelerated the granulation, but had no obvious effect on the bacterial community structure of granules. The whole granulation process could be categorized into three phases, i.e. lag, granulation and granule maturation phase. During lag period GAC provided nuclei for sludge to attach, and thus enhanced the morphological regularization of sludge. During granulation period the granule size increased significantly due to the growth of bacteria in granules. GAC reduced the compression caused by the inter-particle collisions and thus accelerate the granulation. GAC has no negative effect on the performance of SBR, and thus efficient simultaneous removal of COD, nitrogen and phosphorus were obtained during most of the operating time.

Differences in Gut Metabolites and Microbial Composition and Functions between Egyptian and U.S. Children Are Consistent with Their Diets

Previous studies indicated that populations consuming a Mediterranean diet rich in fiber, vegetables, and fruits have a significantly lower risk of cardiovascular and metabolic diseases than populations of industrialized societies consuming diets enriched in processed carbohydrates, animal proteins, and fats. To explore the potential contributions of gut microbiota to the observed diet-related metabolic effects, we conducted an integrative analysis of distal gut microbiota composition and functions and intestinal metabolites in Egyptian and U.S. teenagers. All Egyptian gut microbial communities belonged to the Prevotella enterotype, whereas all but one of the U.S. samples were of the Bacteroides enterotype. The intestinal environment of Egyptians was characterized by higher levels of short-chain fatty acids, a higher prevalence of microbial polysaccharide degradation-encoding genes, and a higher proportion of several polysaccharide-degrading genera. Egyptian gut microbiota also appeared to be under heavier bacteriophage pressure. In contrast, the gut environment of U.S. children was rich in amino acids and lipid metabolism-associated compounds; contained more microbial genes encoding protein degradation, vitamin biosynthesis, and iron acquisition pathways; and was enriched in several protein- and starch-degrading genera. Levels of 1-methylhistamine, a biomarker of allergic response, were elevated in U.S. guts, as were the abundances of members of Faecalibacterium and Akkermansia, two genera with recognized anti-inflammatory effects. The revealed corroborating differences in fecal microbiota structure and functions and metabolite profiles between Egyptian and U.S. teenagers are consistent with the nutrient variation between Mediterranean and Western diets. IMPORTANCE The human gastrointestinal microbiota functions as an important mediator of diet for host metabolism. To evaluate how consumed diets influence the gut environment, we carried out simultaneous interrogations of distal gut microbiota and metabolites in samples from healthy children in Egypt and the United States. While Egyptian children consumed a Mediterranean diet rich in plant foods, U.S. children consumed a Western diet high in animal protein, fats, and highly processed carbohydrates. Consistent with the consumed diets, Egyptian gut samples were enriched in polysaccharide-degrading microbes and end products of polysaccharide fermentation and U.S. gut samples were enriched in proteolytic microbes and end products of protein and fat metabolism. Thus, the intestinal microbiota might be selected on the basis of the diets that we consume, which can open opportunities to affect gut health through modulation of gut microbiota with dietary supplementations.

Evaluation of PacBio sequencing for full-length bacterial 16S rRNA gene classification

Background

Currently, bacterial 16S rRNA gene analyses are based on sequencing of individual variable regions of the 16S rRNA gene (Kozich, et al Appl Environ Microbiol 79:5112–5120, 2013).This short read approach can introduce biases. Thus, full-length bacterial 16S rRNA gene sequencing is needed to reduced biases. A new alternative for full-length bacterial 16S rRNA gene sequencing is offered by PacBio single molecule, real-time (SMRT) technology. The aim of our study was to validate PacBio P6 sequencing chemistry using three approaches: 1) sequencing the full-length bacterial 16S rRNA gene from a single bacterial species Staphylococcus aureus to analyze error modes and to optimize the bioinformatics pipeline; 2) sequencing the full-length bacterial 16S rRNA gene from a pool of 50 different bacterial colonies from human stool samples to compare with full-length bacterial 16S rRNA capillary sequence; and 3) sequencing the full-length bacterial 16S rRNA genes from 11 vaginal microbiome samples and compare with in silico selected bacterial 16S rRNA V1V2 gene region and with bacterial 16S rRNA V1V2 gene regions sequenced using the Illumina MiSeq.

Results

Our optimized bioinformatics pipeline for PacBio sequence analysis was able to achieve an error rate of 0.007% on the Staphylococcus aureus full-length 16S rRNA gene. Capillary sequencing of the full-length bacterial 16S rRNA gene from the pool of 50 colonies from stool identified 40 bacterial species of which up to 80% could be identified by PacBio full-length bacterial 16S rRNA gene sequencing. Analysis of the human vaginal microbiome using the bacterial 16S rRNA V1V2 gene region on MiSeq generated 129 operational taxonomic units (OTUs) from which 70 species could be identified. For the PacBio, 36,000 sequences from over 58,000 raw reads could be assigned to a barcode, and the in silico selected bacterial 16S rRNA V1V2 gene region generated 154 OTUs grouped into 63 species, of which 62% were shared with the MiSeq dataset. The PacBio full-length bacterial 16S rRNA gene datasets generated 261 OTUs, which were grouped into 52 species, of which 54% were shared with the MiSeq dataset. Alpha diversity index reported a higher diversity in the MiSeq dataset.

Conclusion

The PacBio sequencing error rate is now in the same range of the previously widely used Roche 454 sequencing platform and current MiSeq platform. Species-level microbiome analysis revealed some inconsistencies between the full-length bacterial 16S rRNA gene capillary sequencing and PacBio sequencing.

Electronic supplementary material

The online version of this article (doi:10.1186/s12866-016-0891-4) contains supplementary material, which is available to authorized users.

Responses of bacterial community structure and denitrifying bacteria in biofilm to submerged macrophytes and nitrate

Submerged macrophytes play important roles in constructed wetlands and natural water bodies, as these organisms remove nutrients and provide large surfaces for biofilms, which are beneficial for nitrogen removal, particularly from submerged macrophyte-dominated water columns. However, information on the responses of biofilms to submerged macrophytes and nitrogen molecules is limited. In the present study, bacterial community structure and denitrifiers were investigated in biofilms on the leaves of four submerged macrophytes and artificial plants exposed to two nitrate concentrations. The biofilm cells were evenly distributed on artificial plants but appeared in microcolonies on the surfaces of submerged macrophytes. Proteobacteria was the most abundant phylum in all samples, accounting for 27.3–64.8% of the high-quality bacterial reads, followed by Chloroflexi (3.7–25.4%), Firmicutes (3.0–20.1%), Acidobacteria (2.7–15.7%), Actinobacteria (2.2–8.7%), Bacteroidetes (0.5–9.7%), and Verrucomicrobia (2.4–5.2%). Cluster analysis showed that bacterial community structure can be significantly different on macrophytes versus from those on artificial plants. Redundancy analysis showed that electrical conductivity and nitrate concentration were positively correlated with Shannon index and operational taxonomic unit (OTU) richness (log₁₀ transformed) but somewhat negatively correlated with microbial density. The relative abundances of five denitrifying genes were positively correlated with nitrate concentration and electrical conductivity but negatively correlated with dissolved oxygen.

Identification of Comamonas testosteroni as an androgen degrader in sewage

Numerous studies have reported the masculinization of freshwater wildlife exposed to androgens in polluted rivers. Microbial degradation is a crucial mechanism for eliminating steroid hormones from contaminated ecosystems. The aerobic degradation of testosterone was observed in various bacterial isolates. However, the ecophysiological relevance of androgen-degrading microorganisms in the environment is unclear. Here, we investigated the biochemical mechanisms and corresponding microorganisms of androgen degradation in aerobic sewage. Sewage samples collected from the Dihua Sewage Treatment Plant (Taipei, Taiwan) were aerobically incubated with testosterone (1 mM). Androgen metabolite analysis revealed that bacteria adopt the 9, 10-seco pathway to degrade testosterone. A metagenomic analysis indicated the apparent enrichment of Comamonas spp. (mainly C. testosteroni) and Pseudomonas spp. in sewage incubated with testosterone. We used the degenerate primers derived from the meta-cleavage dioxygenase gene (tesB) of various proteobacteria to track this essential catabolic gene in the sewage. The amplified sequences showed the highest similarity (87–96%) to tesB of C. testosteroni. Using quantitative PCR, we detected a remarkable increase of the 16S rRNA and catabolic genes of C. testosteroni in the testosterone-treated sewage. Together, our data suggest that C. testosteroni, the model microorganism for aerobic testosterone degradation, plays a role in androgen biodegradation in aerobic sewage.

Significant loss of sensitivity and specificity in the taxonomic classification occurs when short 16S rRNA gene sequences are used

The classification performance of Kraken was evaluated in terms of sensitivity and specificity when using short and long 16S rRNA sequences. A total of 440,738 sequences from bacteria with complete taxonomic classifications were downloaded from the high quality ribosomal RNA database SILVA. Amplicons produced (86,371 sequences; 1450 bp) by virtual PCR with primers covering the V1–V9 region of the 16S-rRNA gene were used as reference. Virtual PCŔs of internal fragments V3–V4, V4–V5 and V3–V5 were performed. A total of 81,523, 82,334 and 82,998 amplicons were obtained for regions V3–V4, V4–V5 and V3–V5 respectively. Differences in depth of taxonomic classification were detected among the internal fragments. For instance, sensitivity and specificity of sequences classified up to subspecies level were higher when the largest internal fraction (V3–V5) was used (54.0 and 74.6% respectively), compared to V3–V4 (45.1 and 66.7%) and V4–V5 (41.8 and 64.6%) fragments. Similar pattern was detected for sequences classified up to more superficial taxonomic categories (i.e. family, order, class…). Results also demonstrate that internal fragments lost specificity and some could be misclassified at the deepest taxonomic levels (i.e. species or subspecies). It is concluded that the larger V3–V5 fragment could be considered for massive high throughput sequencing reducing the loss of sensitivity and sensibility.

Summer holidays as break-point in shaping a tannery sludge microbial community around a stable core microbiota

Recently, several investigations focused on the discovery of a bacterial consortium shared among different wastewater treatment plants (WWTPs). Nevertheless, the definition of a core microbiota over time represents the necessary counterpart in order to unravel the dynamics of bacterial communities in these environments. Here we performed a monthly survey on the bacterial community of a consortial industrial plant. Objectives of this study were: (1) to identify a core microbiota constant over time; (2) to evaluate the temporal dynamics of the community during one year. A conspicuous and diversified core microbiota is constituted by operational taxonomic units which are present throughout the year in the plant. Community composition data confirm that the presence and abundance of bacteria in WWTPs is highly consistent at high taxonomic level. Our results indicate however a difference in microbial community structure between two groups of samples, identifying the summer holiday period as the break-point. Changes in the structure of the microbial community occur otherwise gradually, one month after another. Further studies will clarify how the size and diversity of the core microbiota could affect the observed dynamics.

Recent developments in detection and enumeration of waterborne bacteria: a retrospective minireview

Waterborne diseases have emerged as global health problems and their rapid and sensitive detection in environmental water samples is of great importance. Bacterial identification and enumeration in water samples is significant as it helps to maintain safe drinking water for public consumption. Culture‐based methods are laborious, time‐consuming, and yield false‐positive results, whereas viable but nonculturable (VBNCs) microorganisms cannot be recovered. Hence, numerous methods have been developed for rapid detection and quantification of waterborne pathogenic bacteria in water. These rapid methods can be classified into nucleic acid‐based, immunology‐based, and biosensor‐based detection methods. This review summarizes the principle and current state of rapid methods for the monitoring and detection of waterborne bacterial pathogens. Rapid methods outlined are polymerase chain reaction (PCR), digital droplet PCR, real‐time PCR, multiplex PCR, DNA microarray, Next‐generation sequencing (pyrosequencing, Illumina technology and genomics), and fluorescence in situ hybridization that are categorized as nucleic acid‐based methods. Enzyme‐linked immunosorbent assay (ELISA) and immunofluorescence are classified into immunology‐based methods. Optical, electrochemical, and mass‐based biosensors are grouped into biosensor‐based methods. Overall, these methods are sensitive, specific, time‐effective, and important in prevention and diagnosis of waterborne bacterial diseases.

Spatial and temporal investigation of the microbiome of the Caribbean octocoral Erythropodium caribaeorum

The octocoral Erythropodium caribaeorum is an important species in the Caribbean coral reef community and a source of the cytotoxic natural product desmethyleleutherobin. We utilized 16S small subunit rRNA gene amplicon pyrosequencing to characterize the microbiome of E. caribaeorum collected from Florida, USA and San Salvador, The Bahamas at multiple time points. This coral was found to have a very high microbial richness with an average Chao1 estimated richness of 1464 ± 707 operational taxonomic units and average Shannon diversity index of 4.26 ± 1.65. The taxonomic class Gammaproteobacteria was a dominant member in all samples and the genus Endozoicomonas accounted for an average of 37.7% ± 30.0% of the total sequence reads. One Endozoicomonas sp. was found to be a stable member of all E. caribaeorum sequence libraries regardless of location or time of collection and accounted for 30.1% of all sequence reads. This is the first report characterizing the microbiome associated with the encrusting octocoral E. caribaeorum.

Development of an Analysis Pipeline Characterizing Multiple Hypervariable Regions of 16S rRNA Using Mock Samples

Objectives

There is much speculation on which hypervariable region provides the highest bacterial specificity in 16S rRNA sequencing. The optimum solution to prevent bias and to obtain a comprehensive view of complex bacterial communities would be to sequence the entire 16S rRNA gene; however, this is not possible with second generation standard library design and short-read next-generation sequencing technology.

Methods

This paper examines a new process using seven hypervariable or V regions of the 16S rRNA (six amplicons: V2, V3, V4, V6-7, V8, and V9) processed simultaneously on the Ion Torrent Personal Genome Machine (Life Technologies, Grand Island, NY). Four mock samples were amplified using the 16S Ion Metagenomics Kit^™ (Life Technologies) and their sequencing data is subjected to a novel analytical pipeline.

Results

Results are presented at family and genus level. The Kullback-Leibler divergence (D_KL), a measure of the departure of the computed from the nominal bacterial distribution in the mock samples, was used to infer which region performed best at the family and genus levels. Three different hypervariable regions, V2, V4, and V6-7, produced the lowest divergence compared to the known mock sample. The V9 region gave the highest (worst) average D_KL while the V4 gave the lowest (best) average D_KL. In addition to having a high D_KL, the V9 region in both the forward and reverse directions performed the worst finding only 17% and 53% of the known family level and 12% and 47% of the genus level bacteria, while results from the forward and reverse V4 region identified all 17 family level bacteria.

Conclusions

The results of our analysis have shown that our sequencing methods using 6 hypervariable regions of the 16S rRNA and subsequent analysis is valid. This method also allowed for the assessment of how well each of the variable regions might perform simultaneously. Our findings will provide the basis for future work intended to assess microbial abundance at different time points throughout a clinical protocol.

Skin Microbiome Surveys Are Strongly Influenced by Experimental Design

Culture-independent studies to characterize skin microbiota are increasingly common, due in part to affordable and accessible sequencing and analysis platforms. Compared to culture-based techniques, DNA sequencing of the bacterial 16S ribosomal RNA (rRNA) gene or whole metagenome shotgun (WMS) sequencing provides more precise microbial community characterizations. Most widely used protocols were developed to characterize microbiota of other habitats (i.e., gastrointestinal) and have not been systematically compared for their utility in skin microbiome surveys. Here we establish a resource for the cutaneous research community to guide experimental design in characterizing skin microbiota. We compare two widely sequenced regions of the 16S rRNA gene to WMS sequencing for recapitulating skin microbiome community composition, diversity, and genetic functional enrichment. We show that WMS sequencing most accurately recapitulates microbial communities, but sequencing of hypervariable regions 1-3 of the 16S rRNA gene provides highly similar results. Sequencing of hypervariable region 4 poorly captures skin commensal microbiota, especially Propionibacterium. WMS sequencing, which is resource and cost intensive, provides evidence of a community's functional potential; however, metagenome predictions based on 16S rRNA sequence tags closely approximate WMS genetic functional profiles. This study highlights the importance of experimental design for downstream results in skin microbiome surveys.