Supplementary Figure 1 Overview of the analysis pipeline. Supplementary Table 1 Details of conventionally raised and conventionalized mouse samples. Supplementary Discussion Expanded discussion of QIIME analyses presented in the main text; Sequencing of 16S rRNA gene amplicons; QIIME analysis notes; Expanded Figure 1 legend; Links to raw data and processed output from the runs with and without denoising.

/pdf/qiime-allows-analysis-of-high-throughput-community-2ob4mkjrat.pdf

QIIME allows analysis of high-throughput community sequencing data.

SILVA (from Latin silva, forest, http://www.arb-silva.de) is a comprehensive web resource for up to date, quality-controlled databases of aligned ribosomal RNA (rRNA) gene sequences from the Bacteria, Archaea and Eukaryota domains and supplementary online services. The referred database release 111 (July 2012) contains 3 194 778 small subunit and 288 717 large subunit rRNA gene sequences. Since the initial description of the project, substantial new features have been introduced, including advanced quality control procedures, an improved rRNA gene aligner, online tools for probe and primer evaluation and optimized browsing, searching and downloading on the website. Furthermore, the extensively curated SILVA taxonomy and the new non-redundant SILVA datasets provide an ideal reference for high-throughput classification of data from next-generation sequencing approaches.

The SILVA ribosomal RNA gene database project: improved data processing and web-based tools

We present the open-source software package DADA2 for modeling and correcting Illumina-sequenced amplicon errors (https://github.com/benjjneb/dada2). DADA2 infers sample sequences exactly and resolves differences of as little as 1 nucleotide. In several mock communities, DADA2 identified more real variants and output fewer spurious sequences than other methods. We applied DADA2 to vaginal samples from a cohort of pregnant women, revealing a diversity of previously undetected Lactobacillus crispatus variants.

/pdf/dada2-high-resolution-sample-inference-from-illumina-ygkdeyoxho.pdf

DADA2: High-resolution sample inference from Illumina amplicon data

Amplified marker-gene sequences can be used to understand microbial community structure, but they suffer from a high level of sequencing and amplification artifacts. The UPARSE pipeline reports operational taxonomic unit (OTU) sequences with ≤1% incorrect bases in artificial microbial community tests, compared with >3% incorrect bases commonly reported by other methods. The improved accuracy results in far fewer OTUs, consistently closer to the expected number of species in a community.

UPARSE: highly accurate OTU sequences from microbial amplicon reads

Background The analysis of microbial communities through DNA sequencing brings many challenges: the integration of different types of data with methods from ecology, genetics, phylogenetics, multivariate statistics, visualization and testing. With the increased breadth of experimental designs now being pursued, project-specific statistical analyses are often needed, and these analyses are often difficult (or impossible) for peer researchers to independently reproduce. The vast majority of the requisite tools for performing these analyses reproducibly are already implemented in R and its extensions (packages), but with limited support for high throughput microbiome census data. Results Here we describe a software project, phyloseq, dedicated to the object-oriented representation and analysis of microbiome census data in R. It supports importing data from a variety of common formats, as well as many analysis techniques. These include calibration, filtering, subsetting, agglomeration, multi-table comparisons, diversity analysis, parallelized Fast UniFrac, ordination methods, and production of publication-quality graphics; all in a manner that is easy to document, share, and modify. We show how to apply functions from other R packages to phyloseq-represented data, illustrating the availability of a large number of open source analysis techniques. We discuss the use of phyloseq with tools for reproducible research, a practice common in other fields but still rare in the analysis of highly parallel microbiome census data. We have made available all of the materials necessary to completely reproduce the analysis and figures included in this article, an example of best practices for reproducible research. Conclusions The phyloseq project for R is a new open-source software package, freely available on the web from both GitHub and Bioconductor.

/pdf/phyloseq-an-r-package-for-reproducible-interactive-analysis-1ez1xx3dlr.pdf

phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data.

The Dietary Guidelines for Americans advises on dietary intake to meet nutritional needs, promote health, and prevent diseases. Diet affects the intestinal microbiota and is increasingly linked to health. It is vital to investigate the relationships between diet quality and the microbiota to better understand the impact of nutrition on human health.This study aimed to investigate the differences in fecal microbiota composition in adults from the American Gut Project based on their adherence to the Dietary Guidelines for Americans.This study was a cross-sectional analysis of the 16S sequencing and food frequency data of a subset of adults (n = 432; age = 18-60 y; 65% female, 89% white) participating in the crowdsourced American Gut Project. The Healthy Eating Index-2015 assessed the compliance with Dietary Guideline recommendations. The cohort was divided into tertiles based on Healthy Eating Index-2015 scores, and differences in taxonomic abundances and diversity were compared between high and low scorers.The mean Total Score for low-scoring adults (58.1 ± 5.4) was comparable with the reported score of the average American adult (56.7). High scorers for the Total Score and components related to vegetables, grains, and dairy had greater alpha diversity than low scorers. High scorers in the fatty acid component had a lower alpha diversity than low scorers (95% CI: 0.35, 1.85). A positive log-fold difference in abundance of plant carbohydrate-metabolizing taxa in the families Lachnospiraceae and Ruminococcaceae was observed in high-scoring tertiles for Total Score, vegetable, fruit, and grain components (Benjamini-Hochberg; q < 0.05).Adults with greater compliance to the Dietary Guidelines demonstrated higher diversity in their fecal microbiota and greater abundance of bacteria capable of metabolizing complex carbohydrates, providing evidence on how Dietary Guidelines support the gut microbiota. 

Diet Quality and the Fecal Microbiota in Adults in the American Gut Project.

In a large cohort of 1,772 participants from the Hispanic Community Health Study/Study of Latinos with overlapping 16SV4 rRNA gene (bacterial amplicon), ITS1 (fungal amplicon), and shotgun sequencing data, we demonstrate that 16SV4 amplicon sequencing and shotgun metagenomics offer the same level of taxonomic accuracy for bacteria at the genus level even at shallow sequencing depths. In contrast, for fungal taxa, we did not observe meaningful agreements between shotgun and ITS1 amplicon results. Finally, we show that amplicon and shotgun data can be harmonized and pooled to yield larger microbiome datasets with excellent agreement (<1% effect size variance across three independent outcomes) using pooled amplicon/shotgun data compared to pure shotgun metagenomic analysis. Thus, there are multiple approaches to study the microbiome in epidemiological studies, and we provide a demonstration of a powerful pooling approach that will allow researchers to leverage the massive amount of amplicon sequencing data generated over the last two decades. 

https://doi.org/10.1016/j.crmeth.2022.100391

Comprehensive evaluation of shotgun metagenomics, amplicon sequencing, and harmonization of these platforms for epidemiological studies

Most experimental sciences now rely on computing, and biological sciences are no exception. As datasets get bigger, so do the computing costs, making proper optimization of the codes used by scientists increasingly important. Many of the codes developed in recent years are based on the Python-based NumPy, due to its ease of use and good performance characteristics. The composable nature of NumPy, however, does not generally play well with the multi-tier nature of modern CPUs, making any non-trivial multi-step algorithm limited by the external memory access speeds, which are hundreds of times slower than the CPU's compute capabilities. In order to fully utilize the CPU compute capabilities, one must keep the working memory footprint small enough to fit in the CPU caches, which requires splitting the problem into smaller portions and fusing together as many steps as possible. In this paper, we present changes based on these principles to two important functions in the scikit-bio library, principal coordinates analysis and the Mantel test, that resulted in over 100x speed improvement in these widely used, general-purpose tools.

Accelerating key bioinformatics tasks 100-fold by improving memory access

https://dl.acm.org/doi/pdf/10.1145/3437359.3465562

Daniel McDonald

Papers

Diet Quality and the Fecal Microbiota in Adults in the American Gut Project.

Comprehensive evaluation of shotgun metagenomics, amplicon sequencing, and harmonization of these platforms for epidemiological studies

Accelerating key bioinformatics tasks 100-fold by improving memory access

Accelerating Key Bioinformatics Tasks 100-fold by Improving Memory Access

Pre-Crystallisation for Improved Lpcvd Deposited Boron-Doped Poly-Si Hole-Selective Contacts