Partial convergence of the human vaginal and rectal maternal microbiota in late gestation and early post-partum.

The human vaginal and fecal microbiota change during pregnancy. Because of the proximity of these perineal sites and the evolutionarily conserved maternal-to-neonatal transmission of the microbiota, we hypothesized that the microbiota of these two sites (rectal and vaginal) converge during the last gestational trimester as part of the preparation for parturition. To test this hypothesis, we analyzed 16S rRNA sequences from vaginal introitus and rectal samples in 41 women at gestational ages 6 and 8 months, and at 2 months post-partum. The results show that the human vaginal and rectal bacterial microbiota converged during the last gestational trimester and into the 2nd month after birth, with a significant decrease in Lactobacillus species in both sites, as alpha diversity progressively increased in the vagina and decreased in the rectum. The microbiota convergence of the maternal vaginal-anal sites perinatally might hold significance for the inter-generational transmission of the maternal microbiota.

Microbiomes in Neverland.

Differentiating the effects of infant microbiota, developmental, and nutritional changes on immunological maturation during weaning is an ongoing challenge. In this issue of Cell Host & Microbe, Lubin and colleagues report a gnotobiotic mouse model that maintains neonatal-like microbiome composition into adulthood to help answer burning questions in this field.

Infant behavioral state and stool microbiome in infants receiving Lactocaseibacillus rhamnosus GG in formula: randomized controlled trial.

BACKGROUND: Our aim was to evaluate infant behavioral state, stool microbiome profile and calprotectin in infants with infantile colic receiving a partially hydrolyzed protein formula with or without added Lacticaseibacillus (formerly Lactobacillus) rhamnosus GG (LGG). METHODS: In this single-center, double-blind, controlled, parallel, prospective study, term infants (14-28 days of age) identified with colic (using modified Wessel's criteria: cried and/or fussed ≥ 3 h/day for ≥ 3 days/week, in a one-week period) were randomized to receive one of two formulas over a three-week feeding period: marketed partially hydrolyzed cow's milk-based infant formula (PHF, n = 35) or a similar formula with added LGG (PHF-LGG, n = 36). Parent-reported infant behavior was recorded at three time points (Study Days 2-4, 10-12, and 18-20). Duration (hours/day) of crying/fussing (averaged over each three-day period) was the primary outcome. Stool samples were collected at Baseline and Study End (Days 19-21) to determine stool LGG colonization (by qPCR) and microbial abundance (using 16S rRNA gene sequencing) and calprotectin (µg/g). RESULTS: Duration of crying/fussing (mean ± SE) decreased and awake/content behavior increased over time with no significant group differences over the course of the study. There were no group differences in the percentage of infants who experienced colic by study end. Colic decreased by Study End vs Baseline in both groups. Change in fecal calprotectin also was similar between groups. Comparing Study End vs Baseline, LGG abundance was greater in the PHF-LGG group (P < 0.001) whereas alpha diversity was greater in the PHF group (P = 0.022). Beta diversity was significantly different between PHF and PHF-LGG at Study End (P = 0.05). By study end, relative abundance of L. rhamnosus was higher in the PHF-LGG vs PHF group and vs Baseline. CONCLUSIONS: In this pilot study of infants with colic, both study formulas were well tolerated. Crying/fussing decreased and awake/content behavior increased in both study groups over the course of the study. Study results demonstrate a successful introduction of the probiotic to the microbiome. The partially hydrolyzed protein formula with added LGG was associated with significant changes in the gut microbiome. TRIAL REGISTRATION: ClinicalTrials.gov, ClinicalTrials.gov Identifier: NCT02340143 . Registered 16/01/2015.

Reproducible acquisition, management and meta-analysis of nucleotide sequence (meta)data using q2-fondue.

MOTIVATION: The volume of public nucleotide sequence data has blossomed over the past two decades and is ripe for re- and meta-analyses to enable novel discoveries. However, reproducible re-use and management of sequence datasets and associated metadata remain critical challenges. We created the open source Python package q2-fondue to enable user-friendly acquisition, re-use and management of public sequence (meta)data while adhering to open data principles. RESULTS: q2-fondue allows fully provenance-tracked programmatic access to and management of data from the NCBI Sequence Read Archive (SRA). Unlike other packages allowing download of sequence data from the SRA, q2-fondue enables full data provenance tracking from data download to final visualization, integrates with the QIIME 2 ecosystem, prevents data loss upon space exhaustion and allows download of (meta)data given a publication library. To highlight its manifold capabilities, we present executable demonstrations using publicly available amplicon, whole genome and metagenome datasets. AVAILABILITY AND IMPLEMENTATION: q2-fondue is available as an open-source BSD-3-licensed Python package at https://github.com/bokulich-lab/q2-fondue. Usage tutorials are available in the same repository. All Jupyter notebooks used in this article are available under https://github.com/bokulich-lab/q2-fondue-examples. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Multi-omics data integration reveals metabolome as the top predictor of the cervicovaginal microenvironment.

Emerging evidence suggests that host-microbe interaction in the cervicovaginal microenvironment contributes to cervical carcinogenesis, yet dissecting these complex interactions is challenging. Herein, we performed an integrated analysis of multiple "omics" datasets to develop predictive models of the cervicovaginal microenvironment and identify characteristic features of vaginal microbiome, genital inflammation and disease status. Microbiomes, vaginal pH, immunoproteomes and metabolomes were measured in cervicovaginal specimens collected from a cohort (n = 72) of Arizonan women with or without cervical neoplasm. Multi-omics integration methods, including neural networks (mmvec) and Random Forest supervised learning, were utilized to explore potential interactions and develop predictive models. Our integrated analyses revealed that immune and cancer biomarker concentrations were reliably predicted by Random Forest regressors trained on microbial and metabolic features, suggesting close correspondence between the vaginal microbiome, metabolome, and genital inflammation involved in cervical carcinogenesis. Furthermore, we show that features of the microbiome and host microenvironment, including metabolites, microbial taxa, and immune biomarkers are predictive of genital inflammation status, but only weakly to moderately predictive of cervical neoplastic disease status. Different feature classes were important for prediction of different phenotypes. Lipids (e.g. sphingolipids and long-chain unsaturated fatty acids) were strong predictors of genital inflammation, whereas predictions of vaginal microbiota and vaginal pH relied mostly on alterations in amino acid metabolism. Finally, we identified key immune biomarkers associated with the vaginal microbiota composition and vaginal pH (MIF), as well as genital inflammation (IL-6, IL-10, MIP-1α).

RESCRIPt: Reproducible sequence taxonomy reference database management.

Nucleotide sequence and taxonomy reference databases are critical resources for widespread applications including marker-gene and metagenome sequencing for microbiome analysis, diet metabarcoding, and environmental DNA (eDNA) surveys. Reproducibly generating, managing, using, and evaluating nucleotide sequence and taxonomy reference databases creates a significant bottleneck for researchers aiming to generate custom sequence databases. Furthermore, database composition drastically influences results, and lack of standardization limits cross-study comparisons. To address these challenges, we developed RESCRIPt, a Python 3 software package and QIIME 2 plugin for reproducible generation and management of reference sequence taxonomy databases, including dedicated functions that streamline creating databases from popular sources, and functions for evaluating, comparing, and interactively exploring qualitative and quantitative characteristics across reference databases. To highlight the breadth and capabilities of RESCRIPt, we provide several examples for working with popular databases for microbiome profiling (SILVA, Greengenes, NCBI-RefSeq, GTDB), eDNA and diet metabarcoding surveys (BOLD, GenBank), as well as for genome comparison. We show that bigger is not always better, and reference databases with standardized taxonomies and those that focus on type strains have quantitative advantages, though may not be appropriate for all use cases. Most databases appear to benefit from some curation (quality filtering), though sequence clustering appears detrimental to database quality. Finally, we demonstrate the breadth and extensibility of RESCRIPt for reproducible workflows with a comparison of global hepatitis genomes. RESCRIPt provides tools to democratize the process of reference database acquisition and management, enabling researchers to reproducibly and transparently create reference materials for diverse research applications. RESCRIPt is released under a permissive BSD-3 license at https://github.com/bokulich-lab/RESCRIPt.

Beating Naive Bayes at Taxonomic Classification of 16S rRNA Gene Sequences.

Naive Bayes classifiers (NBC) have dominated the field of taxonomic classification of amplicon sequences for over a decade. Apart from having runtime requirements that allow them to be trained and used on modest laptops, they have persistently provided class-topping classification accuracy. In this work we compare NBC with random forest classifiers, neural network classifiers, and a perfect classifier that can only fail when different species have identical sequences, and find that in some practical scenarios there is little scope for improving on NBC for taxonomic classification of 16S rRNA gene sequences. Further improvements in taxonomy classification are unlikely to come from novel algorithms alone, and will need to leverage other technological innovations, such as ecological frequency information.

Experiences and lessons learned from two virtual, hands-on microbiome bioinformatics workshops.

In October of 2020, in response to the Coronavirus Disease 2019 (COVID-19) pandemic, our team hosted our first fully online workshop teaching the QIIME 2 microbiome bioinformatics platform. We had 75 enrolled participants who joined from at least 25 different countries on 6 continents, and we had 22 instructors on 4 continents. In the 5-day workshop, participants worked hands-on with a cloud-based shared compute cluster that we deployed for this course. The event was well received, and participants provided feedback and suggestions in a postworkshop questionnaire. In January of 2021, we followed this workshop with a second fully online workshop, incorporating lessons from the first. Here, we present details on the technology and protocols that we used to run these workshops, focusing on the first workshop and then introducing changes made for the second workshop. We discuss what worked well, what didn't work well, and what we plan to do differently in future workshops.

Sources and Assembly of Microbial Communities in Vineyards as a Functional Component of Winegrowing.

Microbiomes are integral to viticulture and winemaking - collectively termed winegrowing - where diverse fungi and bacteria can exert positive and negative effects on grape health and wine quality. Wine is a fermented natural product, and the vineyard serves as a key point of entry for quality-modulating microbiota, particularly in wine fermentations that are conducted without the addition of exogenous yeasts. Thus, the sources and persistence of wine-relevant microbiota in vineyards critically impact its quality. Site-specific variations in microbiota within and between vineyards may contribute to regional wine characteristics. This includes distinctions in microbiomes and microbiota at the strain level, which can contribute to wine flavor and aroma, supporting the role of microbes in the accepted notion of terroir as a biological phenomenon. Little is known about the factors driving microbial biodiversity within and between vineyards, or those that influence annual assembly of the fruit microbiome. Fruit is a seasonally ephemeral, yet annually recurrent product of vineyards, and as such, understanding the sources of microbiota in vineyards is critical to the assessment of whether or not microbial terroir persists with inter-annual stability, and is a key factor in regional wine character, as stable as the geographic distances between vineyards. This review examines the potential sources and vectors of microbiota within vineyards, general rules governing plant microbiome assembly, and how these factors combine to influence plant-microbe interactions relevant to winemaking.

Measuring the microbiome: Best practices for developing and benchmarking microbiomics methods.

Microbiomes are integral components of diverse ecosystems, and increasingly recognized for their roles in the health of humans, animals, plants, and other hosts. Given their complexity (both in composition and function), the effective study of microbiomes (microbiomics) relies on the development, optimization, and validation of computational methods for analyzing microbial datasets, such as from marker-gene (e.g., 16S rRNA gene) and metagenome data. This review describes best practices for benchmarking and implementing computational methods (and software) for studying microbiomes, with particular focus on unique characteristics of microbiomes and microbiomics data that should be taken into account when designing and testing microbiomics methods.

A total crapshoot? Evaluating bioinformatic decisions in animal diet metabarcoding analyses.

Metabarcoding studies provide a powerful approach to estimate the diversity and abundance of organisms in mixed communities in nature. While strategies exist for optimizing sample and sequence library preparation, best practices for bioinformatic processing of amplicon sequence data are lacking in animal diet studies. Here we evaluate how decisions made in core bioinformatic processes, including sequence filtering, database design, and classification, can influence animal metabarcoding results. We show that denoising methods have lower error rates compared to traditional clustering methods, although these differences are largely mitigated by removing low-abundance sequence variants. We also found that available reference datasets from GenBank and BOLD for the animal marker gene cytochrome oxidase I (COI) can be complementary, and we discuss methods to improve existing databases to include versioned releases. Taxonomic classification methods can dramatically affect results. For example, the commonly used Barcode of Life Database (BOLD) Classification API assigned fewer names to samples from order through species levels using both a mock community and bat guano samples compared to all other classifiers (vsearch-SINTAX and q2-feature-classifier's BLAST + LCA, VSEARCH + LCA, and Naive Bayes classifiers). The lack of consensus on bioinformatics best practices limits comparisons among studies and may introduce biases. Our work suggests that biological mock communities offer a useful standard to evaluate the myriad computational decisions impacting animal metabarcoding accuracy. Further, these comparisons highlight the need for continual evaluations as new tools are adopted to ensure that the inferences drawn reflect meaningful biology instead of digital artifacts.

Phylogenetic farming: Can evolutionary history predict crop rotation via the soil microbiome?

Agriculture has long employed phylogenetic rules whereby farmers are encouraged to rotate taxonomically unrelated plants in shared soil. Although this forms a central tenet of sustainable agriculture, strangely, this on-farm "rule of thumb" has never been rigorously tested in a scientific framework. To experimentally evaluate the relationship between phylogenetic distance and crop performance, we used a plant-soil feedback approach whereby 35 crops and weeds varying in their relatedness to tomato (Solanum lycopersicum) were tested in a two-year field experiment. We used community profiling of the bacteria and fungi to determine the extent to which soil microbes contribute to phenotypic differences in crop growth. Overall, tomato yield was ca. 15% lower in soil previously cultivated with tomato; yet, past the species level there was no effect of phylogenetic distance on crop performance. Soil microbial communities, on the other hand, were compositionally more similar between close plant relatives. Random forest regression predicted log10 phylogenetic distance to tomato with moderate accuracy (R 2 = .52), primarily driven by bacteria in the genus Sphingobium. These data indicate that, beyond avoiding conspecifics, evolutionary history contributes little to understanding plant-soil feedbacks in agricultural fields; however, microbial legacies can be predicted by species identity and relatedness.

An Early Pandemic Analysis of SARS-CoV-2 Population Structure and Dynamics in Arizona.

In December of 2019, a novel coronavirus, SARS-CoV-2, emerged in the city of Wuhan, China, causing severe morbidity and mortality. Since then, the virus has swept across the globe, causing millions of confirmed infections and hundreds of thousands of deaths. To better understand the nature of the pandemic and the introduction and spread of the virus in Arizona, we sequenced viral genomes from clinical samples tested at the TGen North Clinical Laboratory, the Arizona Department of Health Services, and those collected as part of community surveillance projects at Arizona State University and the University of Arizona. Phylogenetic analysis of 84 genomes from across Arizona revealed a minimum of 11 distinct introductions inferred to have occurred during February and March. We show that >80% of our sequences descend from strains that were initially circulating widely in Europe but have since dominated the outbreak in the United States. In addition, we show that the first reported case of community transmission in Arizona descended from the Washington state outbreak that was discovered in late February. Notably, none of the observed transmission clusters are epidemiologically linked to the original travel-related case in the state, suggesting successful early isolation and quarantine. Finally, we use molecular clock analyses to demonstrate a lack of identifiable, widespread cryptic transmission in Arizona prior to the middle of February 2020.IMPORTANCE As the COVID-19 pandemic swept across the United States, there was great differential impact on local and regional communities. One of the earliest and hardest hit regions was in New York, while at the same time Arizona (for example) had low incidence. That situation has changed dramatically, with Arizona now having the highest rate of disease increase in the country. Understanding the roots of the pandemic during the initial months is essential as the pandemic continues and reaches new heights. Genomic analysis and phylogenetic modeling of SARS-COV-2 in Arizona can help to reconstruct population composition and predict the earliest undetected introductions. This foundational work represents the basis for future analysis and understanding as the pandemic continues.

QIIME 2 Enables Comprehensive End-to-End Analysis of Diverse Microbiome Data and Comparative Studies with Publicly Available Data.

QIIME 2 is a completely re-engineered microbiome bioinformatics platform based on the popular QIIME platform, which it has replaced. QIIME 2 facilitates comprehensive and fully reproducible microbiome data science, improving accessibility to diverse users by adding multiple user interfaces. QIIME 2 can be combined with Qiita, an open-source web-based platform, to re-use available data for meta-analysis. The following basic protocol describes how to install QIIME 2 on a single computer and analyze microbiome sequence data, from processing of raw DNA sequence reads through generating publishable interactive figures. These interactive figures allow readers of a study to interact with data with the same ease as its authors, advancing microbiome science transparency and reproducibility. We also show how plug-ins developed by the community to add analysis capabilities can be installed and used with QIIME 2, enhancing various aspects of microbiome analyses-e.g., improving taxonomic classification accuracy. Finally, we illustrate how users can perform meta-analyses combining different datasets using readily available public data through Qiita. In this tutorial, we analyze a subset of the Early Childhood Antibiotics and the Microbiome (ECAM) study, which tracked the microbiome composition and development of 43 infants in the United States from birth to 2 years of age, identifying microbiome associations with antibiotic exposure, delivery mode, and diet. For more information about QIIME 2, see https://qiime2.org. To troubleshoot or ask questions about QIIME 2 and microbiome analysis, join the active community at https://forum.qiime2.org. © 2020 The Authors. Basic Protocol: Using QIIME 2 with microbiome data Support Protocol: Further microbiome analyses.

Reproducibly sampling SARS-CoV-2 genomes across time, geography, and viral diversity.

The COVID-19 pandemic has led to a rapid accumulation of SARS-CoV-2 genomes, enabling genomic epidemiology on local and global scales. Collections of genomes from resources such as GISAID must be subsampled to enable computationally feasible phylogenetic and other analyses. We present genome-sampler, a software package that supports sampling collections of viral genomes across multiple axes including time of genome isolation, location of genome isolation, and viral diversity. The software is modular in design so that these or future sampling approaches can be applied independently and combined (or replaced with a random sampling approach) to facilitate custom workflows and benchmarking. genome-sampler is written as a QIIME 2 plugin, ensuring that its application is fully reproducible through QIIME 2's unique retrospective data provenance tracking system. genome-sampler can be installed in a conda environment on macOS or Linux systems. A complete default pipeline is available through a Snakemake workflow, so subsampling can be achieved using a single command. genome-sampler is open source, free for all to use, and available at https://caporasolab.us/genome-sampler. We hope that this will facilitate SARS-CoV-2 research and support evaluation of viral genome sampling approaches for genomic epidemiology.

Learning representations of microbe-metabolite interactions.

Integrating multiomics datasets is critical for microbiome research; however, inferring interactions across omics datasets has multiple statistical challenges. We solve this problem by using neural networks (https://github.com/biocore/mmvec) to estimate the conditional probability that each molecule is present given the presence of a specific microorganism. We show with known environmental (desert soil biocrust wetting) and clinical (cystic fibrosis lung) examples, our ability to recover microbe-metabolite relationships, and demonstrate how the method can discover relationships between microbially produced metabolites and inflammatory bowel disease.

Species abundance information improves sequence taxonomy classification accuracy.

Popular naive Bayes taxonomic classifiers for amplicon sequences assume that all species in the reference database are equally likely to be observed. We demonstrate that classification accuracy degrades linearly with the degree to which that assumption is violated, and in practice it is always violated. By incorporating environment-specific taxonomic abundance information, we demonstrate a significant increase in the species-level classification accuracy across common sample types. At the species level, overall average error rates decline from 25% to 14%, which is favourably comparable to the error rates that existing classifiers achieve at the genus level (16%). Our findings indicate that for most practical purposes, the assumption that reference species are equally likely to be observed is untenable. q2-clawback provides a straightforward alternative for samples from common environments.

Author Correction: Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.