DADA2: High-resolution sample inference from Illumina amplicon data

Assessment of urinary 3-indoxyl sulfate as a marker for gut microbiota diversity and abundance of Clostridiales

OBJECTIVES

After allogeneic hematopoietic stem cell transplantation (allo-HCT), urinary levels of 3-indoxyl sulfate (3-IS) correlate with the relative abundance of bacteria from the class Clostridia (RAC), and antibiotic treatment is considered the major determinant of this outcome. A high RAC has been associated with favorable outcome after allo-HCT and protection from Clostridium difficile infection (CDI). We assessed correlations between alpha diversity, RAC and urinary 3-IS levels in a non-allo-HCT clinical cohort of antibiotic treated patients to further explore 3-IS as a biomarker of reduced diversity and predisposition to CDI.

METHODS

Fecal and urinary specimens were analyzed from 40 non-allo-HCT hospitalized patients before and 9 ± 2 days after initiation of intravenous antibiotic treatment. Fecal microbiota were analyzed by 16s RNA sequencing and urinary 3-IS was analyzed by liquid chromatography-tandem mass spectrometry. Receiver operating characteristic (ROC) analysis was performed to assess the predictive value of 3-IS.

RESULTS

At a RAC cutoff of <30%, the binary logarithm of 3-IS (medium 3-IS: ≤2.5; high 3-IS: >2.5) was predictive with an accuracy of 82% (negative predictive value: 87%, positive predictive value 67%). Accuracy was improved by combing antibiotic history with 3-IS levels (accuracy 89%, npv 88%, ppv 92%).

CONCLUSION

In conjunction with patient antibiotic history, 3-IS is a candidate marker to predict RAC.

Links between environment, diet, and the hunter-gatherer microbiome

The study of traditional populations provides a view of human-associated microbes unperturbed by industrialization, as well as a window into the microbiota that co-evolved with humans. Here we discuss our recent work characterizing the microbiota from the Hadza hunter-gatherers of Tanzania. We found seasonal shifts in bacterial taxa, diversity, and carbohydrate utilization by the microbiota. When compared to the microbiota composition from other populations around the world, the Hadza microbiota shares bacterial families with other traditional societies that are rare or absent from microbiotas of industrialized nations. We present additional observations from the Hadza microbiota and their lifestyle and environment, including microbes detected on hands, water, and animal sources, how the microbiota varies with sex and age, and the short-term effects of introducing agricultural products into the diet. In the context of our previously published findings and of these additional observations, we discuss a path forward for future work.

Biodiversity estimates and ecological interpretations of meiofaunal communities are biased by the taxonomic approach

Accurate assessments of biodiversity are crucial to advising ecosystem-monitoring programs and understanding ecosystem function. Nevertheless, a standard operating procedure to assess biodiversity accurately and consistently has not been established. This is especially true for meiofauna, a diverse community (>20 phyla) of small benthic invertebrates that have fundamental ecological roles. Recent studies show that metabarcoding is a cost-effective and time-effective method to estimate meiofauna biodiversity, in contrast to morphological-based taxonomy. Here, we compare biodiversity assessments of a diverse meiofaunal community derived by applying multiple taxonomic methods based on comparative morphology, molecular phylogenetic analysis, DNA barcoding of individual specimens, and metabarcoding of environmental DNA. We show that biodiversity estimates are strongly biased across taxonomic methods and phyla. Such biases affect understanding of community structures and ecological interpretations. This study supports the urgency of improving aspects of environmental high-throughput sequencing and the value of taxonomists in correctly understanding biodiversity estimates.

Identification of Key Bacteria Involved in the Induction of Incident Bacterial Vaginosis: A Prospective Study

Background

The sequence of events preceding incident bacterial vaginosis (iBV) is unclear.

Methods

African American women who have sex with women, who had no Amsel criteria and Nugent scores of 0-3, were followed for 90 days to detect iBV (defined as a Nugent score of 7-10 on at least 2-3 consecutive days), using self-collected vaginal swab specimens. For women with iBV (cases) and women maintaining normal vaginal flora (healthy women), 16S ribosomal RNA gene sequencing targeting V4 was performed. Longitudinal vaginal microbiome data were analyzed.

Results

Of 204 women screened, 42 enrolled; of these, 45% developed iBV. Sequencing was performed on 448 specimens from 14 cases and 8 healthy women. Among healthy women, Lactobacillus crispatus dominated the vaginal microbiota in 75%. In contrast, prior to iBV, the vaginal microbiota in 79% of cases was dominated by Lactobacillus iners and/or Lactobacillus jensenii/Lactobacillus gasseri. The mean relative abundance of Prevotella bivia, Gardnerella vaginalis, Atopobium vaginae, and Megasphaera type I became significantly higher in cases 4 days before (P. bivia), 3 days before (G. vaginalis), and on the day of (A. vaginae and Megasphaera type I) iBV onset. The mean relative abundance of Sneathia sanguinegens, Finegoldia magna, BV-associated bacteria 1-3, and L. iners was not significantly different between groups before onset of iBV.

Conclusion

G. vaginalis, P. bivia, A. vaginae, and Megasphaera type I may play significant roles in iBV.

Fecal microbiota transplantation in a kidney transplant recipient with recurrent urinary tract infection

PURPOSE

We report on a kidney transplant recipient treated with fecal microbiota transplantation (FMT) for recurrent urinary tract infections.

METHODS

FMT was administered via frozen capsulized microbiota. Before and after FMT, urinary, fecal and vaginal microbiota compositions were analyzed.

RESULTS

The patient remained without symptoms after FMT.

CONCLUSIONS

Underlying mechanisms of action need to be addressed in depth by future research.

Breast Cancer and Its Relationship with the Microbiota

The microorganisms that live symbiotically in human beings are increasingly recognized as important players in health and disease. The largest collection of these microorganisms is found in the gastrointestinal tract. Microbial composition reflects both genetic and lifestyle variables of the host. This microbiota is in a dynamic balance with the host, exerting local and distant effects. Microbial perturbation (dysbiosis) could contribute to the risk of developing health problems. Various bacterial genes capable of producing estrogen-metabolizing enzymes have been identified. Accordingly, gut microbiota is capable of modulating estrogen serum levels. Conversely, estrogen-like compounds may promote the proliferation of certain species of bacteria. Therefore, a crosstalk between microbiota and both endogenous hormones and estrogen-like compounds might synergize to provide protection from disease but also to increase the risk of developing hormone-related diseases. Recent research suggests that the microbiota of women with breast cancer differs from that of healthy women, indicating that certain bacteria may be associated with cancer development and with different responses to therapy. In this review, we discuss recent knowledge about the microbiome and breast cancer, identifying specific characteristics of the human microbiome that may serve to develop novel approaches for risk assessment, prevention and treatment for this disease.

The Microbial Landscape of Sea Stars and the Anatomical and Interspecies Variability of Their Microbiome

Sea stars are among the most important predators in benthic ecosystems worldwide which is partly attributed to their unique gastrointestinal features and feeding behaviors. Despite their ecological importance, the microbiome of these animals and its influence on adult host health and development largely remains unknown. To begin to understand such interactions we sought to understand what bacteria are associated with these animals, how the microbiome is partitioned across regions of the body and how seawater influences their microbiome. We analyzed the microbiome composition of a geographically and taxonomically diverse set of sea star taxa by using 16S rRNA gene amplicon sequencing and compared microorganisms associated with different regions of their body and to their local environment. In addition, we estimated the bacterial and coelomocyte abundance in the sea star coelomic fluid and bacterioplankton abundance in the surrounding seawater via epifluorescence microscopy. The average bacterial cell abundance observed in the coelomic fluid was one to two orders of magnitude lower than the bacterioplankton abundance in the surrounding seawater suggesting a selection against the presence of microorganisms in the coelomic fluid. The sea star microbiome was also significantly different from seawater with relatively few shared microbial taxa. Microbial communities were found to be significantly different between the pyloric caeca, gonads, coelomic fluid, and body wall of the animals. The most noticeable difference between anatomical sites was the greater relative abundance of Spirochaetae and Tenericutes found in hard tissues (gonads, pyloric caeca, and body wall) than in the coelomic fluid. The microbiome of sea stars thus appears to be anatomically partitioned, distinct from the microbial community of seawater and contains a relatively low abundance of bacteria within the coelomic cavity.

Dynamics and Functional Potential of Stormwater Microorganisms Colonizing Sand Filters

Stormwater management is increasingly relying on engineered infiltration systems (EIS) to reduce the volume and improve the quality of managed stormwater. Yet, EIS in the field will be colonized by a diverse array of environmental microorganisms that change the physiochemical properties of the EIS and provide a habitat for microorganisms with harmful or beneficial qualities. Understanding factors influencing the composition and stability of microbial communities could open up strategies for more efficient management of stormwater. Here, we analyzed the potential pathogenic and metabolic capabilities of stormwater microorganisms colonizing idealized EIS (i.e., sand columns) under laboratory conditions over time. The diversity of microbial communities was analyzed using 16S rRNA gene sequencing, and potential pathogens and denitrifying microbes were identified from taxonomic match to known species. Denitrification potential as determined by nosZ abundance was also assessed with quantitative polymerase chain reaction PCR. Our findings demonstrate that replicate microbial communities colonizing sand columns change in a similar way over time, distinct from control columns and the source community. Potential pathogens were initially more abundant on the columns than in the stormwater but returned to background levels by 24 days after inoculation. The conditions within sand columns select for potential denitrifying microorganisms, some of which were also potential pathogens. These results demonstrate that a diverse suite of stormwater microorganisms colonize sand filters, including a transient population of potential pathogens and denitrifiers. Manipulating the inoculating microbial community of EIS could prove an effective mechanism for changing both potential pathogens and denitrifying bacteria.

Initial description of the core ocular surface microbiome in dogs: Bacterial community diversity and composition in a defined canine population

OBJECTIVE

To characterize the bacterial community residing on the conjunctiva of clinically healthy dogs.

METHODS

Bacterial DNA from conjunctival swabs of 10 dogs with normal ocular examinations (both OD and OS, n = 20) was extracted, and 16S rRNA amplicons were sequenced using Illumina MiSeq 600. Resulting data were subjected to quality control steps, and analyzed for bacterial community richness and diversity, within- and between-group dissimilarity, and relative taxonomic composition.

RESULTS

High-quality reads (2.22 million bp) resulted in a mean of 159 068 sequences per sample. Bacterial community evenness and diversity was high when compared to other species, and did not significantly differ when samples were grouped by dogs or eyes. As expected, within-dog samples were more similar than between-dog samples. Taxonomic classification revealed that >95% of the community consisted of Firmicutes (34.9 ± 8.8%), Actinobacteria (26.3 ± 7.1%), Proteobacteria (26.2 ± 6.6%), and Bacteroidetes (9.4 ± 2.4%). Key members of the dog ocular surface microbiome, found in all dogs and corresponding to >25% of all identified OTUs (operational taxonomic units), were part of the Bifidobacteriaceae, Lachnospiraceae, Moraxellaceae, Corynebacteriaceae families. Genera previously thought to account for the majority of the core ocular surface microbiome in the dog (Staphylococcus sp., Streptococcus sp., and Bacillus sp.) were associated with only 2.63% of overall reads.

CONCLUSIONS

This study shows the feasibility of conjunctival swabs and high-throughput sequencing to profile the bacterial community structure of the canine ocular surface. A core ocular surface microbiome was identified for this canine population.

Impact of lincosamides antibiotics on the composition of the rat gut microbiota and the metabolite profile of plasma and feces

The importance of the gut microorganisms and their wide range of interactions with the host are well-acknowledged. In this study, lincomycin and clindamycin were used to modulate microbial communities of Wistar rats to gain a comprehensive understanding of the implications of microbiome alterations. A metabolomics approach and taxonomic profiling were applied to characterize the effects of these antibiotics on the functionality of the microbiome and to identify microbiome-related metabolites. After treatment, the diversity of the microbial community was drastically reduced. Bacteroidetes and Verrucomicrobia were drastically reduced, Tenericutes and Deferribacteres completely disappeared, while abundance of Firmicutes and Proteobacteria were highly increased. Changes in plasma and feces metabolites were observed for metabolites belonging mainly to the class of complex lipids, fatty acids and related metabolites as well as amino acids and related compounds. Bile acid metabolism was markedly affected: taurocholic acid, glycochenodeoxycholic acid and cholic acid presented abrupt changes showing a specific metabolite pattern indicating disruption of the microbial community. In both plasma and feces taurocholic acid was highly upregulated upon treatment whereas glycochenodeoxycholic acid was downregulated. Cholic acid was upregulated in feces but downregulated in plasma. These results show that changes in the gut microbial community lead to alterations of the metabolic profile in blood and feces of the host and can be used to identify potentially microbiome-related metabolites. This implies that metabolomics could be a suitable tool to estimate the extent of changes induced in the intestinal microbiome with respect to consequences for the host.

Effects of transportation to and co-mingling at an auction market on nasopharyngeal and tracheal bacterial communities of recently weaned beef cattle

The objective was to study effects of transportation to and co-mingling at an auction market on nasopharyngeal and tracheal bacterial communities of feedlot cattle. Two groups of 30 Angus-cross heifers were studied from weaning to 28 d after arrival at a feedlot. For each group, half the heifers were either transported directly to a feedlot after weaning (RANC) or transported to and co-mingled at an auction market for 24 h before being placed in a feedlot (AUCT). Deep nasal swabs (DNS) and trans-tracheal aspirates (TTA) were collected at weaning (d0) and at on-arrival processing at the feedlot (d2). At 7 (d9) and 28 d (d30) after arrival, DNS were repeated. The DNA was extracted from DNS and TTA and the V4 region of the 16S rRNA gene sequenced (MiSeq). Alpha diversity analysis did not reveal differences between AUCT and RANC. However, bacterial diversity decreased over time in the nasopharynx, especially at d9. Although beta-diversity was not different between AUCT and RANC, interval after arrival and feedlot where heifers were placed affected composition of the nasopharyngeal bacterial communities. In both groups, a large increase in Mycoplasma was observed after arrival; in one group, Mycoplasma bovis was dominant at d9 and remained dominant until d30. However, in the other group, Mycoplasma dispar dominated at d9 and was replaced by Moraxella at d30. We concluded that transportation to and co-mingling at an auction market for 24 h did not significantly influence diversity or composition of nasopharyngeal or tracheal bacterial communities.

Stationary and portable sequencing-based approaches for tracing wastewater contamination in urban stormwater systems

Urban sewer systems consist of wastewater and stormwater sewers, of which only wastewater is processed before being discharged. Occasionally, misconnections or damages in the network occur, resulting in untreated wastewater entering natural water bodies via the stormwater system. Cultivation of faecal indicator bacteria (e.g. Escherichia coli; E. coli) is the current standard for tracing wastewater contamination. This method is cheap but has limited specificity and mobility. Here, we compared the E. coli culturing approach with two sequencing-based methodologies (Illumina MiSeq 16S rRNA gene amplicon sequencing and Oxford Nanopore MinION shotgun metagenomic sequencing), analysing 73 stormwater samples collected in Stockholm. High correlations were obtained between E. coli culturing counts and frequencies of human gut microbiome amplicon sequences, indicating E. coli is indeed a good indicator of faecal contamination. However, the amplicon data further holds information on contamination source or alternatively how much time has elapsed since the faecal matter has entered the system. Shotgun metagenomic sequencing on a subset of the samples using a portable real-time sequencer, MinION, correlated well with the amplicon sequencing data. This study demonstrates the use of DNA sequencing to detect human faecal contamination in stormwater systems and the potential of tracing faecal contamination directly in the field.

Effect of Antibiotic-Mediated Microbiome Modulation on Rotavirus Vaccine Immunogenicity: A Human, Randomized-Control Proof-of-Concept Trial

Rotavirus vaccines (RVV) protect against childhood gastroenteritis caused by rotavirus (RV) but have decreased effectiveness in low- and middle-income settings. This proof-of-concept, randomized-controlled, open-label trial tested if microbiome modulation can improve RVV immunogenicity. Healthy adults were randomized and administered broad-spectrum (oral vancomycin, ciprofloxacin, metronidazole), narrow-spectrum (vancomycin), or no antibiotics and then vaccinated with RVV, 21 per group per protocol. Baseline anti-RV IgA was high in all subjects. Although antibiotics did not alter absolute anti-RV IgA titers, RVV immunogenicity was boosted at 7 days in the narrow-spectrum group. Further, antibiotics increased fecal shedding of RV while also rapidly altering gut bacterial beta diversity. Beta diversity associated with RVV immunogenicity boosting at day 7 and specific bacterial taxa that distinguish RVV boosters and RV shedders were identified. Despite the negative primary endpoint, this study demonstrates that microbiota modification alters the immune response to RVV and supports further exploration of microbiome manipulation to improve RVV immunogenicity.

Denoising the Denoisers: an independent evaluation of microbiome sequence error-correction approaches

High-depth sequencing of universal marker genes such as the 16S rRNA gene is a common strategy to profile microbial communities. Traditionally, sequence reads are clustered into operational taxonomic units (OTUs) at a defined identity threshold to avoid sequencing errors generating spurious taxonomic units. However, there have been numerous bioinformatic packages recently released that attempt to correct sequencing errors to determine real biological sequences at single nucleotide resolution by generating amplicon sequence variants (ASVs). As more researchers begin to use high resolution ASVs, there is a need for an in-depth and unbiased comparison of these novel “denoising” pipelines. In this study, we conduct a thorough comparison of three of the most widely-used denoising packages (DADA2, UNOISE3, and Deblur) as well as an open-reference 97% OTU clustering pipeline on mock, soil, and host-associated communities. We found from the mock community analyses that although they produced similar microbial compositions based on relative abundance, the approaches identified vastly different numbers of ASVs that significantly impact alpha diversity metrics. Our analysis on real datasets using recommended settings for each denoising pipeline also showed that the three packages were consistent in their per-sample compositions, resulting in only minor differences based on weighted UniFrac and Bray–Curtis dissimilarity. DADA2 tended to find more ASVs than the other two denoising pipelines when analyzing both the real soil data and two other host-associated datasets, suggesting that it could be better at finding rare organisms, but at the expense of possible false positives. The open-reference OTU clustering approach identified considerably more OTUs in comparison to the number of ASVs from the denoising pipelines in all datasets tested. The three denoising approaches were significantly different in their run times, with UNOISE3 running greater than 1,200 and 15 times faster than DADA2 and Deblur, respectively. Our findings indicate that, although all pipelines result in similar general community structure, the number of ASVs/OTUs and resulting alpha-diversity metrics varies considerably and should be considered when attempting to identify rare organisms from possible background noise.

Real-Time Imaging Revealed That Exoelectrogens from Wastewater Are Selected at the Center of a Gradient Electric Field

Exoelectrogens acclimated from the environment are the key to energy recovery from waste in bioelectrochemical systems. However, it is still unknown how these bacteria are selectively enriched on the electrode. Here we confirmed for the first time that the electric field (EF) intensity selects exoelectrogens from wastewater using an integrated electrovisual system with a gradient EF. Under the operating conditions ( I = 3 × 10^-3A), the EF intensity on the working electrode ranged from 6.00 V/cm at the center to 1.08 V/cm at the edge. A thick biofilm (88.9 μm) with spherical pink aggregates was observed at the center, while the color became gray at the edge (33.8 μm). The coverage of the biofilm also increased linearly with EF intensity from 0.42 at the edge (12 mm to the center) to 0.78 at the center. The biofilm at the center contained 76% Geobacter, which was 25% higher than that at the edge (60%). Geobacter anodireducens was the main species induced by the EF (50% at the center vs 24% at the edge). These results improve our fundamental knowledge of exoelectrogen acclimation and mixed electroactive biofilm formation, which has broader implications for energy recovery from waste and general understanding of microbial ecology.

Assessing and Interpreting the Within-Body Biogeography of Human Microbiome Diversity

A human body hosts a relatively independent microbiome including five major regional biomes (i.e., airway, oral, gut, skin, and urogenital). Each of them may possess different regional characteristics with important implications to our health and diseases (i.e., so-termed microbiome associated diseases). Nevertheless, these regional microbiomes are connected with each other through diffusions and migrations. Here, we investigate the within-body (intra-individual) distribution feature of microbiome diversity via diversity area relationship (DAR) modeling, which, to the best of our knowledge, has not been systematically studied previously. We utilized the Hill numbers for measuring alpha and beta-diversities and built 1,200 within-body DAR models with to date the most comprehensive human microbiome datasets of 18 sites from the human microbiome project (HMP) cohort. We established the intra-DAR profile (z-q pattern: the diversity scaling parameter z of the power law (PL) at diversity order q = 0-3), intra-PDO (pair-wise diversity overlap) profile (g-q), and intra-MAD (maximal accrual diversity) profile (D max-q) for the within-body biogeography of the human microbiome. These profiles constitute the "maps" of the within-body biogeography, and offer important insights on the within-body distribution of the human microbiome. Furthermore, we investigated the heterogeneity among individuals in their biogeography parameters and found that there is not an "average Joe" that can represent majority of individuals in a cohort or population. For example, we found that most individuals in the HMP cohort have relatively lower maximal accrual diversity (MAD) or in the "long tail" of the so-termed power law distribution. In the meantime, there are a small number of individuals in the cohort who possess disproportionally higher MAD values. These findings may have important implications for personalized medicine of the human microbiome associated diseases in practice, besides their theoretical significance in microbiome research such as establishing the baseline for the conservation of human microbiome.

Comparison of methods for phylogenetic B-cell lineage inference using time-resolved antibody repertoire simulations (AbSim)

Motivation

The evolution of antibody repertoires represents a hallmark feature of adaptive B-cell immunity. Recent advancements in high-throughput sequencing have dramatically increased the resolution to which we can measure the molecular diversity of antibody repertoires, thereby offering for the first time the possibility to capture the antigen-driven evolution of B cells. However, there does not exist a repertoire simulation framework yet that enables the comparison of commonly utilized phylogenetic methods with regard to their accuracy in inferring antibody evolution.

Results

Here, we developed AbSim, a time-resolved antibody repertoire simulation framework, which we exploited for testing the accuracy of methods for the phylogenetic reconstruction of B-cell lineages and antibody molecular evolution. AbSim enables the (i) simulation of intermediate stages of antibody sequence evolution and (ii) the modeling of immunologically relevant parameters such as duration of repertoire evolution, and the method and frequency of mutations. First, we validated that our repertoire simulation framework recreates replicates topological similarities observed in experimental sequencing data. Second, we leveraged Absim to show that current methods fail to a certain extent to predict the true phylogenetic tree correctly. Finally, we formulated simulation-validated guidelines for antibody evolution, which in the future will enable the development of accurate phylogenetic methods.

Availability and implementation

https://cran.r-project.org/web/packages/AbSim/index.html.

Contact

sai.reddy@ethz.ch.

Supplementary information

Supplementary data are available at Bioinformatics online.

Emergent simplicity in microbial community assembly

A major unresolved question in microbiome research is whether the complex taxonomic architectures observed in surveys of natural communities can be explained and predicted by fundamental, quantitative principles. Bridging theory and experiment is hampered by the multiplicity of ecological processes that simultaneously affect community assembly in natural ecosystems. We addressed this challenge by monitoring the assembly of hundreds of soil- and plant-derived microbiomes in well-controlled minimal synthetic media. Both the community-level function and the coarse-grained taxonomy of the resulting communities are highly predictable and governed by nutrient availability, despite substantial species variability. By generalizing classical ecological models to include widespread nonspecific cross-feeding, we show that these features are all emergent properties of the assembly of large microbial communities, explaining their ubiquity in natural microbiomes.

Comparison of Endophytic Microbial Community in Kiwifruit Plant Cultivars

The microbiome makes a significant contribution to plant health and endophytes may generate positive effects for the host. However, there is a limited knowledge available concerning the kiwifruit endophyte. Therefore, we discuss endophyte microbiome community structures among the kiwifruit cultivars. Total reads numbered 17620 in cv. Hayward, 11515 in cv. Haegeum and 13613 in cv. Jecygold. The number of OTUs was follows: 112 in cv. Hageum; and 87 in cvs. Hayward and Jecygold. Most of the identified OTUs were phylum Proteobacteria and it emerged that Actinobacteria, Firmicute and Bacteroidetes were mainly present. While the largest difference in Proteobacteria content is between cv. Haegeum and cv. Jecygold, they shared more OTUs than cv. Hayward. As well, this study revealed the presence of remarkably dominated OTU of Pseudomonas spp. in kiwifruit sap regardless of cultivars. To the best of our knowledge, this study is the first tone to investigate the kiwifruit endophyte-microbiome community.

Nitrogen fixation in a landrace of maize is supported by a mucilage-associated diazotrophic microbiota

Plants are associated with a complex microbiota that contributes to nutrient acquisition, plant growth, and plant defense. Nitrogen-fixing microbial associations are efficient and well characterized in legumes but are limited in cereals, including maize. We studied an indigenous landrace of maize grown in nitrogen-depleted soils in the Sierra Mixe region of Oaxaca, Mexico. This landrace is characterized by the extensive development of aerial roots that secrete a carbohydrate-rich mucilage. Analysis of the mucilage microbiota indicated that it was enriched in taxa for which many known species are diazotrophic, was enriched for homologs of genes encoding nitrogenase subunits, and harbored active nitrogenase activity as assessed by acetylene reduction and 15N2 incorporation assays. Field experiments in Sierra Mixe using 15N natural abundance or 15N-enrichment assessments over 5 years indicated that atmospheric nitrogen fixation contributed 29%-82% of the nitrogen nutrition of Sierra Mixe maize.

Stool Microbiota at Neutrophil Recovery Is Predictive for Severe Acute Graft vs Host Disease After Hematopoietic Cell Transplantation

Background

Graft-versus-host disease (GVHD) is common after allogeneic hematopoietic cell transplantation (HCT). Risk for death from GVHD has been associated with low bacterial diversity in the stool microbiota early after transplant; however, the specific species associated with GVHD risk remain poorly defined.

Methods

We prospectively collected serial weekly stool samples from 66 patients who underwent HCT, starting pre-transplantation and continuing weekly until 100 days post-transplant, a total of 694 observations in HCT recipients. We used 16S rRNA gene polymerase chain reaction with degenerate primers, followed by high-throughput sequencing to assess the relative abundance of sequence reads from bacterial taxa in stool samples over time.

Results

The gut microbiota was highly dynamic in HCT recipients, with loss and appearance of taxa common on short time scales. As in prior studies, GVHD was associated with lower alpha diversity of the stool microbiota. At neutrophil recovery post-HCT, the presence of oral Actinobacteria and oral Firmicutes in stool was positively correlated with subsequent GVHD; Lachnospiraceae were negatively correlated. A gradient of bacterial species (difference of the sum of the relative abundance of positive correlates minus the sum of the relative abundance of negative correlates) was most predictive (receiver operator characteristic area under the curve of 0.83) of subsequent severe acute GVHD.

Conclusions

The stool microbiota around the time of neutrophil recovery post-HCT is predictive of subsequent development of severe acute GVHD in this study.

Tucumã Oil Shifted Ruminal Fermentation, Reducing Methane Production and Altering the Microbiome but Decreased Substrate Digestibility Within a RUSITEC Fed a Mixed Hay - Concentrate Diet

Tucumã oil is sourced from the fruit pulp of the tucumã tree and contains high concentrations of unsaturated fatty acids and carotenoids. Due to these properties it may have the potential to decrease enteric methane (CH4) from ruminants when included in the diet. The objective of this study was to determine the effect of oil mechanically extracted from the fruit pulp of tucumã on fermentation characteristics, CH4 production and the microbial community using the rumen stimulation technique. Treatments consisted of a control diet (forage:concentrate; 70:30), and tucumã oil included at 0.5 or 1.0% (v/v). Addition of tucumã oil linearly decreased (P < 0.01) dry matter disappearance. Total gas (mL/d) and carbon dioxide (CO2) production (mL/d, mL/g DM) were unaffected (P ≥ 0.36) to increasing addition of tucumã oil where 0.5% (v/v) of Tucumã oil numerically increased both variables. Acetate and butyrate percentages of total VFA were linearly decreased (P ≤ 0.01) and propionate and valerate percentages of total VFA were linearly increased (P < 0.01) by increasing concentrations of tucumã oil added to the substrate. The ratio of acetate to propionate was linearly decreased (P < 0.01) with increasing concentration of tucumã oil. Methane production (mL/d) was linearly decreased (P = 0.04) with increasing addition of tucumã oil to the substrate. Tucumã oil reduced the bacterial richness and diversity when included at 1.0% (v/v) in both solid- and liquid- associated microbes. The abundance of the genera Fibrobacter and Rikenellaceae RC9 gut group were decreased and Pyramidobacter, Megasphaera, Anaerovibrio, and Selenomonas were enriched by the addition of 1.0% tucumã oil. In conclusion, tucumã oil resulted in the favorable shift in fermentation products away from acetate toward propionate, decreasing the production of CH4 when tucumã oil was included at 1.0% (v/v), however, substrate digestibility was also inhibited. The rumen microbiota was also altered by the addition of tucumã oil.

Microbial Community Changes in a Chlorinated Solvents Polluted Aquifer Over the Field Scale Treatment With Poly-3-Hydroxybutyrate as Amendment

This study investigated the organohalide-respiring bacteria (OHRB) and the supporting microbial populations operating in a pilot scale plant employing poly-3-hydroxybutyrate (PHB), a biodegradable polymer produced by bacteria from waste streams, for the in situ bioremediation of groundwater contaminated by chlorinated solvents. The bioremediation was performed in ground treatment units, including PHB reactors as slow release source of electron donors, where groundwater extracted from the wells flows through before the re-infiltration to the low permeability zones of the aquifer. The coupling of the biological treatment with groundwater recirculation allowed to drastically reducing the contamination level and the remediation time by efficiently stimulating the growth of autochthonous OHRB and enhancing the mobilization of the pollutants. Quantitative PCR performed along the external treatment unit showed that the PHB reactor may efficiently act as an external incubator to growing Dehalococcoides mccartyi, known to be capable of fully converting chlorinated ethenes to innocuous end-products. The slow release source of electron donors for the bioremediation process allowed the establishment of a stable population of D. mccartyi, mainly carrying bvcA and vcrA genes which are implicated in the metabolic conversion of vinyl chloride to harmless ethene. Next generation sequencing was performed to analyze the phylogenetic diversity of the groundwater microbiome before and after the bioremediation treatment and allowed the identification of the microorganisms working closely with organohalide-respiring bacteria.

Effects of Sub-Chronic MPTP Exposure on Behavioral and Cognitive Performance and the Microbiome of Wild-Type and mGlu8 Knockout Female and Male Mice

Motor dysfunction is a hallmark of Parkinson's disease (PD); however, non-motor symptoms such as gastrointestinal dysfunction often arise prior to motor symptoms. Alterations in the gut microbiome have been proposed as the earliest event in PD pathogenesis. PD symptoms often demonstrate sex differences. Glutamatergic neurotransmission has long been linked to PD pathology. Metabotropic glutamate receptors (mGlu), a family of G protein-coupled receptors, are divided into three groups, with group III mGlu receptors mainly localized presynaptically where they can inhibit glutamate release in the CNS as well as in the gut. Additionally, the gut microbiome can communicate with the CNS via the gut-brain axis. Here, we assessed whether deficiency of metabotropic glutamate receptor 8 (mGlu8), group III mGlu, modulates the effects of the neurotoxin, 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP), on behavioral and cognitive performance in female and male mice. We studied whether these effects are associated with changes in striatal tyrosine hydroxylase (TH) levels and the gut microbiome. Two-week sub-chronic MPTP increased activity of female and male wild-type (WT) and mGlu8 knockout (KO) mice in the open field. MPTP also showed genotype- and sex-dependent effects. MPTP increased the time WT, but not KO, females and males spent exploring objects. In WT mice, MPTP improved sensorimotor function in males but impaired it in females. Further, MPTP impaired cued fear memory in WT, but not KO, male mice. MPTP reduced striatal TH levels in WT and KO mice but these effects were only pronounced in males. MPTP treatment and genotype affected the diversity of the gut microbiome. In addition, there were significant associations between microbiome α-diversity and sensorimotor performance, as well as microbiome composition and fear learning. These results indicate that specific taxa may directly affect motor and fear learning or that the same physiological effects that enhance both forms of learning also alter diversity of the gut microbiome. MPTP's effect on motor and cognitive performance may then be, at least in part, be mediated by the gut microbiome. These data also support mGlu8 as a novel therapeutic target for PD and highlight the importance of including both sexes in preclinical studies.

Broadscale Ecological Patterns Are Robust to Use of Exact Sequence Variants versus Operational Taxonomic Units

Microbial ecologists have made exceptional improvements in our understanding of microbiomes in the last decade due to breakthroughs in sequencing technologies. These advances have wide-ranging implications for fields ranging from agriculture to human health. Due to limitations in databases, the majority of microbial ecology studies use a binning approach to approximate taxonomy based on DNA sequence similarity. There remains extensive debate on the best way to bin and approximate this taxonomy. Here we examine two popular approaches using a large field-based data set examining both bacteria and fungi and conclude that there are not major differences in the ecological outcomes. Thus, it appears that standard microbial community analyses are not overly sensitive to the particulars of binning approaches.

Recent discussion focuses on the best method for delineating microbial taxa, based on either exact sequence variants (ESVs) or traditional operational taxonomic units (OTUs) of marker gene sequences. We sought to test if the binning approach (ESVs versus 97% OTUs) affected the ecological conclusions of a large field study. The data set included sequences targeting all bacteria (16S rRNA) and fungi (internal transcribed spacer [ITS]), across multiple environments diverging markedly in abiotic conditions, over three collection times. Despite quantitative differences in microbial richness, we found that all α and β diversity metrics were highly positively correlated (r > 0.90) between samples analyzed with both approaches. Moreover, the community composition of the dominant taxa did not vary between approaches. Consequently, statistical inferences were nearly indistinguishable. Furthermore, ESVs only moderately increased the genetic resolution of fungal and bacterial diversity (1.3 and 2.1 times OTU richness, respectively). We conclude that for broadscale (e.g., all bacteria or all fungi) α and β diversity analyses, ESV or OTU methods will often reveal similar ecological results. Thus, while there are good reasons to employ ESVs, we need not question the validity of results based on OTUs.

IMPORTANCE Microbial ecologists have made exceptional improvements in our understanding of microbiomes in the last decade due to breakthroughs in sequencing technologies. These advances have wide-ranging implications for fields ranging from agriculture to human health. Due to limitations in databases, the majority of microbial ecology studies use a binning approach to approximate taxonomy based on DNA sequence similarity. There remains extensive debate on the best way to bin and approximate this taxonomy. Here we examine two popular approaches using a large field-based data set examining both bacteria and fungi and conclude that there are not major differences in the ecological outcomes. Thus, it appears that standard microbial community analyses are not overly sensitive to the particulars of binning approaches.

Impact of Edible Cricket Consumption on Gut Microbiota in Healthy Adults, a Double-blind, Randomized Crossover Trial

Edible insects are often considered a nutritious, protein-rich, environmentally sustainable alternative to traditional livestock with growing popularity among North American consumers. While the nutrient composition of several insects is characterized, all potential health impacts have not been evaluated. In addition to high protein levels, crickets contain chitin and other fibers that may influence gut health. In this study, we evaluated the effects of consuming 25 grams/day whole cricket powder on gut microbiota composition, while assessing safety and tolerability. Twenty healthy adults participated in this six-week, double-blind, crossover dietary intervention. Participants were randomized into two study arms and consumed either cricket-containing or control breakfast foods for 14 days, followed by a washout period and assignment to the opposite treatment. Blood and stool samples were collected at baseline and after each treatment period to assess liver function and microbiota changes. Results demonstrate cricket consumption is tolerable and non-toxic at the studied dose. Cricket powder supported growth of the probiotic bacterium, Bifidobacterium animalis, which increased 5.7-fold. Cricket consumption was also associated with reduced plasma TNF-α. These data suggest that eating crickets may improve gut health and reduce systemic inflammation; however, more research is needed to understand these effects and underlying mechanisms.

Selection of Appropriate Metagenome Taxonomic Classifiers for Ancient Microbiome Research

Ancient biomolecules from oral and gut microbiome samples have been shown to be preserved in the archaeological record. Studying ancient microbiome communities using metagenomic techniques offers a unique opportunity to reconstruct the evolutionary trajectories of microbial communities through time. DNA accumulates specific damage over time, which could potentially affect taxonomic classification and our ability to accurately reconstruct community assemblages. It is therefore necessary to assess whether ancient DNA (aDNA) damage patterns affect metagenomic taxonomic profiling. Here, we assessed biases in community structure, diversity, species detection, and relative abundance estimates by five popular metagenomic taxonomic classification programs using in silico-generated data sets with and without aDNA damage. Damage patterns had minimal impact on the taxonomic profiles produced by each program, while false-positive rates and biases were intrinsic to each program. Therefore, the most appropriate classification program is one that minimizes the biases related to the questions being addressed.

Metagenomics enables the study of complex microbial communities from myriad sources, including the remains of oral and gut microbiota preserved in archaeological dental calculus and paleofeces, respectively. While accurate taxonomic assignment is essential to this process, DNA damage characteristic of ancient samples (e.g., reduction in fragment size and cytosine deamination) may reduce the accuracy of read taxonomic assignment. Using a set of in silico-generated metagenomic data sets, we investigated how the addition of ancient DNA (aDNA) damage patterns influences microbial taxonomic assignment by five widely used profilers: QIIME/UCLUST, MetaPhlAn2, MIDAS, CLARK-S, and MALT. In silico-generated data sets were designed to mimic dental plaque, consisting of 40, 100, and 200 microbial species/strains, both with and without simulated aDNA damage patterns. Following taxonomic assignment, the profiles were evaluated for species presence/absence, relative abundance, alpha diversity, beta diversity, and specific taxonomic assignment biases. Unifrac metrics indicated that both MIDAS and MetaPhlAn2 reconstructed the most accurate community structure. QIIME/UCLUST, CLARK-S, and MALT had the highest number of inaccurate taxonomic assignments; false-positive rates were highest by CLARK-S and QIIME/UCLUST. Filtering out species present at <0.1% abundance greatly increased the accuracy of CLARK-S and MALT. All programs except CLARK-S failed to detect some species from the input file that were in their databases. The addition of ancient DNA damage resulted in minimal differences in species detection and relative abundance between simulated ancient and modern data sets for most programs. Overall, taxonomic profiling biases are program specific rather than damage dependent, and the choice of taxonomic classification program should be tailored to specific research questions.

IMPORTANCE Ancient biomolecules from oral and gut microbiome samples have been shown to be preserved in the archaeological record. Studying ancient microbiome communities using metagenomic techniques offers a unique opportunity to reconstruct the evolutionary trajectories of microbial communities through time. DNA accumulates specific damage over time, which could potentially affect taxonomic classification and our ability to accurately reconstruct community assemblages. It is therefore necessary to assess whether ancient DNA (aDNA) damage patterns affect metagenomic taxonomic profiling. Here, we assessed biases in community structure, diversity, species detection, and relative abundance estimates by five popular metagenomic taxonomic classification programs using in silico-generated data sets with and without aDNA damage. Damage patterns had minimal impact on the taxonomic profiles produced by each program, while false-positive rates and biases were intrinsic to each program. Therefore, the most appropriate classification program is one that minimizes the biases related to the questions being addressed.

Household triclosan and triclocarban effects on the infant and maternal microbiome

In 2016, the US Food and Drug Administration banned the use of specific microbicides in some household and personal wash products due to concerns that these chemicals might induce antibiotic resistance or disrupt human microbial communities. Triclosan and triclocarban (referred to as TCs) are the most common antimicrobials in household and personal care products, but the extent to which TC exposure perturbs microbial communities in humans, particularly during infant development, was unknown. We conducted a randomized intervention of TC‐containing household and personal care products during the first year following birth to characterize whether TC exposure from wash products perturbs microbial communities in mothers and their infants. Longitudinal survey of the gut microbiota using 16S ribosomal RNA amplicon sequencing showed that TC exposure from wash products did not induce global reconstruction or loss of microbial diversity of either infant or maternal gut microbiotas. Broadly antibiotic‐resistant species from the phylum Proteobacteria, however, were enriched in stool samples from mothers in TC households after the introduction of triclosan‐containing toothpaste. When compared by urinary triclosan level, agnostic to treatment arm, infants with higher triclosan levels also showed an enrichment of Proteobacteria species. Despite the minimal effects of TC exposure from wash products on the gut microbial community of infants and adults, detected taxonomic differences highlight the need for consumer safety testing of antimicrobial self‐care products on the human microbiome and on antibiotic resistance.

CasPER, a method for directed evolution in genomic contexts using mutagenesis and CRISPR/Cas9

Here we describe a method for robust directed evolution using mutagenesis of large sequence spaces in their genomic contexts. The method employs error-prone PCR and Cas9-mediated genome integration of mutant libraries of large-sized donor variants into single or multiple genomic sites with efficiencies reaching 98-99%. From sequencing of genome integrants, we determined that the mutation frequency along the donor fragments is maintained evenly and successfully integrated into the genomic target loci, indicating that there is no bias of mutational load towards the proximity of the double strand break. To validate the applicability of the method for directed evolution of metabolic gene products we engineered two essential enzymes in the mevalonate pathway of Saccharomyces cerevisiae with selected variants supporting up to 11-fold higher production of isoprenoids. Taken together, our method extends on existing CRISPR technologies by facilitating efficient mutagenesis of hundreds of nucleotides in cognate genomic contexts.

The cotton stainer's gut microbiota suppresses infection of a cotransmitted trypanosomatid parasite

The evolutionary and ecological success of many insects is attributed to mutualistic partnerships with bacteria that confer hosts with novel traits including food digestion, nutrient supplementation, detoxification of harmful compounds and defence against natural enemies. Dysdercus fasciatus firebugs (Hemiptera: Pyrrhocoridae), commonly known as cotton stainers, possess a simple but distinctive gut bacterial community including B vitamin-supplementing Coriobacteriaceae symbionts. In addition, their guts are often infested with the intestinal trypanosomatid parasite Leptomonas pyrrhocoris (Kinetoplastida: Trypanosomatidae). In this study, using experimental bioassays and fluorescence in situ hybridization (FISH), we report on the protective role of the D. fasciatus gut bacteria against L. pyrrhocoris. We artificially infected 2nd instars of dysbiotic and symbiotic insects with a parasite culture and measured parasite titres, developmental time and survival rates. Our results show that L. pyrrhocoris infection increases developmental time and slightly modifies the quantitative composition of the gut microbiota. More importantly, we found significantly higher parasite titres and a tendency towards lower survival rates in parasite-infected dysbiotic insects compared to symbiotic controls, indicating that the gut bacteria successfully interfere with the establishment or proliferation of L. pyrrhocoris. The colonization of symbiotic bacteria on the peritrophic matrix along the gut wall, as revealed by FISH, likely acts as a barrier blocking parasite attachment or entry into the hemolymph. Our findings show that in addition to being nutritionally important, D. fasciatus' gut bacteria complement the host's immune system in preventing parasite invasions and that a stable gut microbial community is integral for the host's health.

EukRef-Ciliophora: a manually curated, phylogeny-based database of small subunit rRNA gene sequences of ciliates

High-throughput sequencing (HTS) surveys, among the most common approaches currently used in environmental microbiology, require reliable reference databases to be correctly interpreted. The EukRef Initiative (eukref.org) is a community effort to manually screen available small subunit (SSU) rRNA gene sequences and produce a public, high-quality and informative framework of phylogeny-based taxonomic annotations. In the context of EukRef, we present a database for the monophyletic phylum Ciliophora, one of the most complex, diverse and ubiquitous protist groups. We retrieved more than 11 500 sequences of ciliates present in GenBank (28% from identified isolates and 72% from environmental surveys). Our approach included the inference of phylogenetic trees for every ciliate lineage and produced the largest SSU rRNA tree of the phylum Ciliophora to date. We flagged approximately 750 chimeric or low-quality sequences, improved the classification of 70% of GenBank entries and enriched environmental and literature metadata by 30%. The performance of EukRef-Ciliophora is superior to the current SILVA database in classifying HTS reads from a global marine survey. Comprehensive outputs are publicly available to make the new tool a useful guide for non-specialists and a quick reference for experts.

Variations in activities of sewer biofilms due to ferrous and ferric iron dosing

Addition of ferrous and ferric iron salts to wastewater is a commonly used practice for sulfide abatement in sewer force mains. When iron is added to wastewater where sulfate respiration takes place, it produces ferrous sulfide precipitates with the formed sulfide. The effect of iron addition has traditionally been focused on solely from the perspective of reaction stoichiometry. Possible influences on the microbial communities in biofilms growing in force mains have largely been neglected. In this study the activity and microbiome was examined in three pilot scale force mains conveying real wastewater, two subjected to iron treatment and one operated as an untreated control. Activity was measured on suspended biofilm samples extracted from the experimental setup. The microbiome of the biofilm was analyzed by V3 + V4 16S rDNA sequencing. Correlation analysis of chemical composition of the biofilms and activity measurements for operational taxonomic units of relevance to sulfide and methane production were performed. In conclusion, it was found that both ferrous and ferric treatment reduced sulfate reduction and methane production, and that both iron salts induced significant changes to force main biofilm microbiomes.

MicroPheno: predicting environments and host phenotypes from 16S rRNA gene sequencing using a k-mer based representation of shallow sub-samples

Motivation

Microbial communities play important roles in the function and maintenance of various biosystems, ranging from the human body to the environment. A major challenge in microbiome research is the classification of microbial communities of different environments or host phenotypes. The most common and cost-effective approach for such studies to date is 16S rRNA gene sequencing. Recent falls in sequencing costs have increased the demand for simple, efficient and accurate methods for rapid detection or diagnosis with proved applications in medicine, agriculture and forensic science. We describe a reference- and alignment-free approach for predicting environments and host phenotypes from 16S rRNA gene sequencing based on k-mer representations that benefits from a bootstrapping framework for investigating the sufficiency of shallow sub-samples. Deep learning methods as well as classical approaches were explored for predicting environments and host phenotypes.

Results

A k-mer distribution of shallow sub-samples outperformed Operational Taxonomic Unit (OTU) features in the tasks of body-site identification and Crohn's disease prediction. Aside from being more accurate, using k-mer features in shallow sub-samples allows (i) skipping computationally costly sequence alignments required in OTU-picking and (ii) provided a proof of concept for the sufficiency of shallow and short-length 16S rRNA sequencing for phenotype prediction. In addition, k-mer features predicted representative 16S rRNA gene sequences of 18 ecological environments, and 5 organismal environments with high macro-F1 scores of 0.88 and 0.87. For large datasets, deep learning outperformed classical methods such as Random Forest and Support Vector Machine.

Availability and implementation

The software and datasets are available at https://llp.berkeley.edu/micropheno.

Supplementary information

Supplementary data are available at Bioinformatics online.

Cold adaptation and replicable microbial community development during long-term low-temperature anaerobic digestion treatment of synthetic sewage

The development and activity of a cold-adapting microbial community was monitored during low-temperature anaerobic digestion (LtAD) treatment of wastewater. Two replicate hybrid anaerobic sludge bed-fixed-film reactors treated a synthetic sewage wastewater at 12°C, at organic loading rates of 0.25-1.0 kg chemical oxygen demand (COD) m-3 d-1, over 889 days. The inoculum was obtained from a full-scale anaerobic digestion reactor, which was operated at 37°C. Both LtAD reactors readily degraded the influent with COD removal efficiencies regularly exceeding 78% for both the total and soluble COD fractions. The biomass from both reactors was sampled temporally and tested for activity against hydrolytic and methanogenic substrates at 12°C and 37°C. Data indicated that significantly enhanced low-temperature hydrolytic and methanogenic activity developed in both systems. For example, the hydrolysis rate constant (k) at 12°C had increased 20-30-fold by comparison to the inoculum by day 500. Substrate affinity also increased for hydrolytic substrates at low temperature. Next generation sequencing demonstrated that a shift in a community structure occurred over the trial, involving a 1-log-fold change in 25 SEQS (OTU-free approach) from the inoculum. Microbial community structure changes and process performance were replicable in the LtAD reactors.

Predictive Modeling of Microbiome Data Using a Phylogeny-Regularized Generalized Linear Mixed Model

Recent human microbiome studies have revealed an essential role of the human microbiome in health and disease, opening up the possibility of building microbiome-based predictive models for individualized medicine. One unique characteristic of microbiome data is the existence of a phylogenetic tree that relates all the microbial species. It has frequently been observed that a cluster or clusters of bacteria at varying phylogenetic depths are associated with some clinical or biological outcome due to shared biological function (clustered signal). Moreover, in many cases, we observe a community-level change, where a large number of functionally interdependent species are associated with the outcome (dense signal). We thus develop "glmmTree," a prediction method based on a generalized linear mixed model framework, for capturing clustered and dense microbiome signals. glmmTree uses the similarity between microbiomes, which is defined based on the microbiome composition and the phylogenetic tree, to predict the outcome. The effects of other predictive variables (e.g., age, sex) can be incorporated readily in the regression framework. Additional tuning parameters enable a data-adaptive approach to capture signals at different phylogenetic depth and abundance level. Simulation studies and real data applications demonstrated that "glmmTree" outperformed existing methods in the dense and clustered signal scenarios.

WHAM!: a web-based visualization suite for user-defined analysis of metagenomic shotgun sequencing data

Background

Exploration of large data sets, such as shotgun metagenomic sequence or expression data, by biomedical experts and medical professionals remains as a major bottleneck in the scientific discovery process. Although tools for this purpose exist for 16S ribosomal RNA sequencing analysis, there is a growing but still insufficient number of user-friendly interactive visualization workflows for easy data exploration and figure generation. The development of such platforms for this purpose is necessary to accelerate and streamline microbiome laboratory research.

Results

We developed the Workflow Hub for Automated Metagenomic Exploration (WHAM!) as a web-based interactive tool capable of user-directed data visualization and statistical analysis of annotated shotgun metagenomic and metatranscriptomic data sets. WHAM! includes exploratory and hypothesis-based gene and taxa search modules for visualizing differences in microbial taxa and gene family expression across experimental groups, and for creating publication quality figures without the need for command line interface or in-house bioinformatics.

Conclusions

WHAM! is an interactive and customizable tool for downstream metagenomic and metatranscriptomic analysis providing a user-friendly interface allowing for easy data exploration by microbiome and ecological experts to facilitate discovery in multi-dimensional and large-scale data sets.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4870-z) contains supplementary material, which is available to authorized users.

Network hubs in root-associated fungal metacommunities

BACKGROUND

Although a number of recent studies have uncovered remarkable diversity of microbes associated with plants, understanding and managing dynamics of plant microbiomes remain major scientific challenges. In this respect, network analytical methods have provided a basis for exploring "hub" microbial species, which potentially organize community-scale processes of plant-microbe interactions.

METHODS

By compiling Illumina sequencing data of root-associated fungi in eight forest ecosystems across the Japanese Archipelago, we explored hubs within "metacommunity-scale" networks of plant-fungus associations. In total, the metadata included 8080 fungal operational taxonomic units (OTUs) detected from 227 local populations of 150 plant species/taxa.

RESULTS

Few fungal OTUs were common across all the eight forests. However, in each of the metacommunity-scale networks representing northern four localities or southern four localities, diverse mycorrhizal, endophytic, and pathogenic fungi were classified as "metacommunity hubs," which were detected from diverse host plant taxa throughout a climatic region. Specifically, Mortierella (Mortierellales), Cladophialophora (Chaetothyriales), Ilyonectria (Hypocreales), Pezicula (Helotiales), and Cadophora (incertae sedis) had broad geographic and host ranges across the northern (cool-temperate) region, while Saitozyma/Cryptococcus (Tremellales/Trichosporonales) and Mortierella as well as some arbuscular mycorrhizal fungi were placed at the central positions of the metacommunity-scale network representing warm-temperate and subtropical forests in southern Japan.

CONCLUSIONS

The network theoretical framework presented in this study will help us explore prospective fungi and bacteria, which have high potentials for agricultural application to diverse plant species within each climatic region. As some of those fungal taxa with broad geographic and host ranges have been known to promote the survival and growth of host plants, further studies elucidating their functional roles are awaited.

Exposure to Arsenic Alters the Microbiome of Larval Zebrafish

Exposure to environmental toxins such as heavy metals can perturb the development and stability of microbial communities associated with human or animal hosts. Widespread arsenic contamination in rivers and riparian habitats therefore presents environmental and health concerns for populations living near sources of contamination. To investigate how arsenic affects host microbiomes, we sequenced and characterized the microbiomes of twenty larval zebrafish exposed to three concentrations of arsenic that are found in contaminated water—low (10 ppb), medium (50 ppb), and high (100 ppb) for 20 days. We found that even a small concentration of arsenic changed the overall microbial composition, structure and diversity of microbial communities, causing dysbiosis in developing larval zebrafish microbiota. In addition, we found that a high concentration of arsenic also increased the abundance of a class 1 integron, an integrase-dependent system facilitating the horizontal transfer of genes conferring resistance to heavy metals and antibiotics.

Extraction and 16S rRNA Sequence Analysis of Microbiomes Associated with Rice Roots

Plant roots associate with a wide diversity of bacteria and archaea across the root-soil spectrum. The rhizosphere microbiota, the communities of microbes in the soil adjacent to the root, can contain up to 10 billion bacterial cells per gram of soil (Raynaud and Nunan, 2014) and can play important roles for the fitness of the host plant. Subsets of the rhizospheric microbiota can colonize the root surface (rhizoplane) and the root interior (endosphere), forming an intimate relationship with the host plant. Compositional analysis of these communities is important to develop tools in order to manipulate root-associated microbiota for increased crop productivity. Due to the reduced cost and increasing throughput of next-generation sequencing, major advances in deciphering these communities have recently been achieved, mainly through the use of amplicon sequencing of the 16S rRNA gene. Here we first present a protocol for dissecting the microbiota from various root compartments, developed using rice as a model. We next present a method for amplifying fragments of the 16S rRNA gene using a dual index approach. Finally, we present a simple workflow for analyzing the resulting sequencing data to make ecological inferences.

Alterations in the oral microbiome in HIV-infected participants after antiretroviral therapy administration are influenced by immune status

OBJECTIVE

To characterize the oral bacterial microbiome in HIV-infected participants at baseline and after 24 weeks of EFV/FTC/TDF.

DESIGN

Thirty-five participants co-enrolled in two AIDS Clinical Trials Group (ACTG) studies, A5272 and A5280, with paired saliva samples and complete data sets were assessed.

METHODS

Paired saliva samples were evaluated for bacterial microbiome using 16S rDNA PCR followed by Illumina sequencing. Diversity and differential abundance was compared between groups. A random forest classification scheme was used to determine the contribution of parameters in classifying participants' CD4+ T-cell count.

RESULTS

Bacterial communities demonstrated considerable variability both within participants and between timepoints, although they became more similar after 24 weeks of ART. At baseline, both the number of taxa detected and the average alpha diversity were variable between participants, but did not differ significantly based on CD4+ cell count, viral load or other factors. After 24 weeks of ART samples obtained from participants with persistently low CD4+ T-cell counts had significantly higher bacterial richness and diversity. Several differentially abundant taxa, including Porphyromonas species associated with periodontal disease, were identified, which discriminated between baseline and posttreatment samples. Analysis demonstrated that although inflammatory markers are important in untreated disease, the salivary microbiome may play an important role in CD4+ T-cell count recovery after ART.

CONCLUSION

Shifts in the oral microbiome after ART initiation are complex, and may play an important role in immune function and inflammatory disease.

Microbiome niche modification drives diurnal rumen community assembly, overpowering individual variability and diet effects

Niche modification is a process whereby the activity of organisms modifies their local environment creating new niches for other organisms. This process can have a substantial role in community assembly of gut microbial ecosystems due to their vast and complex metabolic activities. We studied the postprandial diurnal community oscillatory patterns of the rumen microbiome and showed that metabolites produced by the rumen microbiome condition its environment and lead to dramatic diurnal changes in community composition and function. After feeding, microbiome composition undergoes considerable change in its phylogenetic breadth manifested as a significant 3-5-fold change in the relative abundance of methanogenic archaea and main bacterial taxa such as Prevotella, in a manner that was independent of individual host variation and diet. These changes in community composition were accompanied by changes in pH and methane partial pressure, suggesting a strong functional connection. Notably, cross-incubation experiments combining metabolites and organisms from different diurnal time points showed that the metabolites released by microbes are sufficient to reproduce changes in community function comparable to those observed in vivo. These findings highlight microbiome niche modification as a deterministic process that drives diurnal community assembly via environmental filtering.

Endospores and other lysis-resistant bacteria comprise a widely shared core community within the human microbiota

Endospore-formers in the human microbiota are well adapted for host-to-host transmission, and an emerging consensus points to their role in determining health and disease states in the gut. The human gut, more than any other environment, encourages the maintenance of endospore formation, with recent culture-based work suggesting that over 50% of genera in the microbiome carry genes attributed to this trait. However, there has been limited work on the ecological role of endospores and other stress-resistant cellular states in the human gut. In fact, there is no data to indicate whether organisms with the genetic potential to form endospores actually form endospores in situ and how sporulation varies across individuals and over time. Here we applied a culture-independent protocol to enrich for endospores and other stress-resistant cells in human feces to identify variation in these states across people and within an individual over time. We see that cells with resistant states are more likely than those without to be shared among multiple individuals, which suggests that these resistant states are particularly adapted for cross-host dissemination. Furthermore, we use untargeted fecal metabolomics in 24 individuals and within a person over time to show that these organisms respond to shared environmental signals, and in particular, dietary fatty acids, that likely mediate colonization of recently disturbed human guts.

Sediment Microbial Communities Influenced by Cool Hydrothermal Fluid Migration

Cool hydrothermal systems (CHSs) are prevalent across the seafloor and discharge fluid volumes that rival oceanic input from rivers, yet the microbial ecology of these systems are poorly constrained. The Dorado Outcrop on the ridge flank of the Cocos Plate in the northeastern tropical Pacific Ocean is the first confirmed CHS, discharging minimally altered <15°C fluid from the shallow lithosphere through diffuse venting and seepage. In this paper, we characterize the resident sediment microbial communities influenced by cool hydrothermal advection, which is evident from nitrate and oxygen concentrations. 16S rRNA gene sequencing revealed that Thaumarchaea, Proteobacteria, and Planctomycetes were the most abundant phyla in all sediments across the system regardless of influence from seepage. Members of the Thaumarchaeota (Marine Group I), Alphaproteobacteria (Rhodospirillales), Nitrospirae, Nitrospina, Acidobacteria, and Gemmatimonadetes were enriched in the sediments influenced by CHS advection. Of the various geochemical parameters investigated, nitrate concentrations correlated best with microbial community structure, indicating structuring based on seepage of nitrate-rich fluids. A comparison of microbial communities from hydrothermal sediments, seafloor basalts, and local seawater at Dorado Outcrop showed differences that highlight the distinct niche space in CHS. Sediment microbial communities from Dorado Outcrop differ from those at previously characterized, warmer CHS sediment, but are similar to deep-sea sediment habitats with surficial ferromanganese nodules, such as the Clarion Clipperton Zone. We conclude that cool hydrothermal venting at seafloor outcrops can alter the local sedimentary oxidation–reduction pathways, which in turn influences the microbial communities within the fluid discharge affected sediment.

Assessing Cat Flea Microbiomes in Northern and Southern California by 16S rRNA Next-Generation Sequencing

Flea-borne diseases (FBDs) impact both human and animal health worldwide. Because adult fleas are obligately hematophagous and can harbor potential pathogens, fleas act as ectoparasites of vertebrates, as well as zoonotic disease vectors. Cat fleas (Ctenocephalides felis) are important vectors of two zoonotic bacterial genera listed as priority pathogens by the National Institute of Allergy and Infectious Diseases (NIAID-USA): Bartonella spp. and Rickettsia spp., causative agents of bartonelloses and rickettsioses, respectively. In this study, we introduce the first microbiome analysis of C. felis samples from California, determining the presence and abundance of relevant pathogenic genera by characterizing the cat flea microbiome through 16S rRNA next-generation sequencing (16S-NGS). Samples from both northern (NoCal) and southern (SoCal) California were assessed to expand current knowledge regarding FBDs in the state. We identified Rickettsia and Bartonella, as well as the endosymbiont Wolbachia, as the most abundant genera, followed by less abundant taxa. In comparison to our previous study screening Californian cat fleas for rickettsiae using PCR/digestion/sequencing of the ompB gene, the 16S-NGS approach applied herein showed a 95% level of agreement in detecting Rickettsia spp. There was no overall difference in microbiome diversity between NoCal and SoCal samples. Bacterial taxa identified by 16S-NGS in this study may help to improve epidemiological investigations, pathogen surveillance efforts, and clinical diagnostics of FBDs in California and elsewhere.

Distinct bacterial metacommunities inhabit the upper and lower respiratory tracts of healthy feedlot cattle and those diagnosed with bronchopneumonia

Specific nasopharyngeal bacterial communities can provide colonization resistance against respiratory pathogens in cattle. However, the role of bacterial communities of the lower airways in respiratory health remains largely unknown. Therefore, our objective was to compare nasopharyngeal and tracheal bacterial communities between healthy feedlot cattle and those with bronchopneumonia (BP). Deep nasal swabs and trans-tracheal aspiration samples were collected from steers with (n = 60) and without (n = 60) BP at 4 feedlots in Western Canada. After DNA extraction, 16S rRNA gene (V4) was amplified and sequenced. Alpha-diversity analysis revealed a lower bacterial diversity in the nasopharynx and trachea of steers with BP compared to healthy pen-mates. Bacterial communities present within the airways clustered into 4 distinct metacommunities that were associated with sampling locations and health status. Metacommunity 1, enriched with Mycoplasma bovis, Mannheimia haemolytica and Pasteurella multocida, was dominant in the nasopharynx and trachea of steers with BP. In contrast, metacommunity 3, enriched with Mycoplasma dispar, Lactococcus lactis and Lactobacillus casei, was mostly present in the trachea of healthy steers. Metacommunity 4, enriched with Corynebacterium, Jeotgalicoccus, Psychrobacter and Planomicrobium, was present in the nasopharynx only. Metacommunity 2, enriched with Histophilus somni, Moraxella and L. lactis, was present in both healthy and sick steers, but was primarily detected in one feedlot. We concluded that distinct bacterial metacommunities inhabited the nasopharynx and trachea of healthy feedlot cattle and those with BP. Because L. lactis and L. casei can inhibit M. haemolytica growth in vitro, their presence in healthy steers may have provided colonization resistance against bacterial respiratory pathogens.

Comparing Sanger sequencing and high-throughput metabarcoding for inferring photobiont diversity in lichens

The implementation of HTS (high-throughput sequencing) approaches is rapidly changing our understanding of the lichen symbiosis, by uncovering high bacterial and fungal diversity, which is often host-specific. Recently, HTS methods revealed the presence of multiple photobionts inside a single thallus in several lichen species. This differs from Sanger technology, which typically yields a single, unambiguous algal sequence per individual. Here we compared HTS and Sanger methods for estimating the diversity of green algal symbionts within lichen thalli using 240 lichen individuals belonging to two species of lichen-forming fungi. According to HTS data, Sanger technology consistently yielded the most abundant photobiont sequence in the sample. However, if the second most abundant photobiont exceeded 30% of the total HTS reads in a sample, Sanger sequencing generally failed. Our results suggest that most lichen individuals in the two analyzed species, Lasallia hispanica and L. pustulata, indeed contain a single, predominant green algal photobiont. We conclude that Sanger sequencing is a valid approach to detect the dominant photobionts in lichen individuals and populations. We discuss which research areas in lichen ecology and evolution will continue to benefit from Sanger sequencing, and which areas will profit from HTS approaches to assessing symbiont diversity.

The gut microbiome of nonhuman primates: Lessons in ecology and evolution

The mammalian gastrointestinal (GI) tract is home to trillions of bacteria that play a substantial role in host metabolism and immunity. While progress has been made in understanding the role that microbial communities play in human health and disease, much less attention has been given to host-associated microbiomes in nonhuman primates (NHPs). Here we review past and current research exploring the gut microbiome of NHPs. First, we summarize methods for characterization of the NHP gut microbiome. Then we discuss variation in gut microbiome composition and function across different NHP taxa. Finally, we highlight how studying the gut microbiome offers new insights into primate nutrition, physiology, and immune system function, as well as enhances our understanding of primate ecology and evolution. Microbiome approaches are useful tools for studying relevant issues in primate ecology. Further study of the gut microbiome of NHPs will offer new insight into primate ecology and evolution as well as human health.

Gut microbiome and magnetic resonance spectroscopy study of subjects at ultra-high risk for psychosis may support the membrane hypothesis

Background:

The microbiota-gut-brain axis and membrane dysfunction in the brain has attracted increasing attention in the field of psychiatric research. However, the possible interactive role of gut microbiota and brain function in the prodromal stage of schizophrenia has not been studied yet.

Methods:

To explore this, we collected fecal samples and performed Magnetic Resonance Spectroscopy (MRS) scans in 81 high risk (HR) subjects, 19 ultra-high risk (UHR) subjects and 69 health controls (HC). Then we analyzed the differences in gut microbiota and choline concentrations in the anterior cingulate cortex (ACC).

Results:

Presences of the orders Clostridiales, Lactobacillales and Bacteroidales were observed at increase levels in fecal samples of UHR subjects compared to the other two groups. The composition changes of gut microbiota indicate the increased production of Short Chain Fatty Acids (SCFAs), which could activate microglia and then disrupt membrane metabolism. Furthermore, this was confirmed by an increase of choline levels, a brain imaging marker of membrane dysfunction, which is also significantly elevated in UHR subjects compared to the HR and HC groups.

Conclusion:

Both gut microbiome and imaging studies of UHR subjects suggest the membrane dysfunction in the brain and hence might support the membrane hypothesis of schizophrenia.

Fecal Microbiome Among Nursing Home Residents with Advanced Dementia and Clostridium difficile

BACKGROUND/OBJECTIVES

Patients colonized with toxinogenic strains of Clostridium difficile have an increased risk of subsequent infection. Given the potential role of the gut microbiome in increasing the risk of C. difficile colonization, we assessed the diversity and composition of the gut microbiota among long-term care facility (LTCF) residents with advanced dementia colonized with C. difficile.

DESIGN

Retrospective analysis of rectal samples collected during a prospective observational study.

SETTING

Thirty-five nursing homes in Boston, Massachusetts.

PARTICIPANTS

Eighty-seven LTCF residents with advanced dementia.

MEASUREMENTS

Operational taxonomic units were identified using 16S rRNA sequencing. Samples positive for C. difficile were matched to negative controls in a 1:3 ratio and assessed for differences in alpha diversity, beta diversity, and differentially abundant features.

RESULTS

Clostridium difficile sequence variants were identified among 7/87 (8.04%) residents. No patient had evidence of C. difficile infection. Demographic characteristics and antimicrobial exposure were similar between the seven cases and 21 controls. The overall biodiversity among cases and controls was reduced with a median Shannon index of 3.2 (interquartile range 2.7-3.9), with no statistically significant differences between groups. The bacterial community structure was significantly different among residents with C. difficile colonization versus those without and included a predominance of Akkermansia spp., Dermabacter spp., Romboutsia spp., Meiothermus spp., Peptoclostridium spp., and Ruminococcaceae UGC 009.

CONCLUSION

LTCF residents with advanced dementia have substantial dysbiosis of their gut microbiome. Specific taxa characterized C. difficile colonization status.