Current challenges and best-practice protocols for microbiome analysis

Abstract Analyzing the microbiome of diverse species and environments using next-generation sequencing techniques has significantly enhanced our understanding on metabolic, physiological and ecological roles of environmental microorganisms. However, the analysis of the microbiome is affected by experimental conditions (e.g. sequencing errors and genomic repeats) and computationally intensive and cumbersome downstream analysis (e.g. quality control, assembly, binning and statistical analyses). Moreover, the introduction of new sequencing technologies and protocols led to a flood of new methodologies, which also have an immediate effect on the results of the analyses. The aim of this work is to review the most important workflows for 16S rRNA sequencing and shotgun and long-read metagenomics, as well as to provide best-practice protocols on experimental design, sample processing, sequencing, assembly, binning, annotation and visualization. To simplify and standardize the computational analysis, we provide a set of best-practice workflows for 16S rRNA and metagenomic sequencing data (available at https://github.com/grimmlab/MicrobiomeBestPracticeReview).

Download Full-text

Comparison between 16S rRNA and shotgun sequencing data for the taxonomic characterization of the gut microbiota

Scientific Reports ◽

10.1038/s41598-021-82726-y ◽

2021 ◽

Vol 11 (1) ◽

Author(s):

Francesco Durazzi ◽

Claudia Sala ◽

Gastone Castellani ◽

Gerardo Manfreda ◽

Daniel Remondini ◽

...

Keyword(s):

16S Rrna ◽

Gut Microbiota ◽

16S Rrna Gene ◽

Gene Sequencing ◽

16S Rrna Gene Sequencing ◽

Shotgun Sequencing ◽

Rrna Gene ◽

Metagenomic Sequencing ◽

Experimental Conditions ◽

Rrna Gene Sequencing

AbstractIn this paper we compared taxonomic results obtained by metataxonomics (16S rRNA gene sequencing) and metagenomics (whole shotgun metagenomic sequencing) to investigate their reliability for bacteria profiling, studying the chicken gut as a model system. The experimental conditions included two compartments of gastrointestinal tracts and two sampling times. We compared the relative abundance distributions obtained with the two sequencing strategies and then tested their capability to distinguish the experimental conditions. The results showed that 16S rRNA gene sequencing detects only part of the gut microbiota community revealed by shotgun sequencing. Specifically, when a sufficient number of reads is available, Shotgun sequencing has more power to identify less abundant taxa than 16S sequencing. Finally, we showed that the less abundant genera detected only by shotgun sequencing are biologically meaningful, being able to discriminate between the experimental conditions as much as the more abundant genera detected by both sequencing strategies.

Download Full-text

Yoghurt Consumption is Associated With Transient Changes in the Composition of the Human Gut Microbiome

10.21203/rs.3.rs-99718/v1 ◽

2020 ◽

Author(s):

Caroline Ivanne Le Roy ◽

Alexander Kurilshikov ◽

Emily Leeming ◽

Alessia Visconti ◽

Ruth Bowyer ◽

...

Keyword(s):

16S Rrna ◽

Gut Microbiota ◽

Visceral Fat ◽

Gut Microbiome ◽

Body Weight Gain ◽

Healthy Eating Index ◽

Rrna Gene ◽

Metagenomic Sequencing ◽

Sequencing Data ◽

Shotgun Metagenomic Sequencing

Abstract Background: Yoghurt contains live bacteria that could contribute via modulation of the gut microbiota to its reported beneficial effects such as reduced body weight gain and lower incidence of type 2 diabetes. To date, the association between yoghurt consumption and the composition of the gut microbiota is underexplored. Here we used clinical variables, metabolomics, 16S rRNA and shotgun metagenomic sequencing data collected on over 1000 predominantly female UK twins to define the link between the gut microbiota and yoghurt-associated health benefits. Results: According to food frequency questionnaires (FFQ), 73% of subjects consumed yoghurt. Consumers presented a healthier diet pattern (healthy eating index: beta = 2.17±0.34; P = 2.72x10-10) and improved metabolic health characterised by reduced visceral fat (beta = -28.18±11.71 g; P = 0.01). According to 16S rRNA gene analyses and whole shotgun metagenomic sequencing approach consistent taxonomic variations were observed with yoghurt consumption. More specifically, we identified higher abundance of species used as yoghurt starters Streptococcus thermophilus (beta = 0.41±0.051; P = 6.14x10-12) and sometimes added Bifidobacterium animalis subsp. lactis (beta = 0.30±0.052; P = 1.49x10-8) in the gut of yoghurt consumers. Replication in 1103 volunteers from the LifeLines-DEEP cohort confirmed the increase of S. thermophilus among yoghurt consumers. Using food records collected the day prior to faecal sampling we showed that increase in these two yoghurt bacteria could be transient. Metabolomics analysis revealed that B. animalis subsp. lactis was associated with 13 faecal metabolites including a 3-hydroxyoctanoic acid, known to be involved in the regulation of gut inflammation.Conclusions: Yoghurt consumption is associated with reduced visceral fat mass and changes in gut microbiome including transient increase of yoghurt-contained species (i.e. S. thermophilus and B. lactis).

Download Full-text

Best practices for variant calling in clinical sequencing

Genome Medicine ◽

10.1186/s13073-020-00791-w ◽

2020 ◽

Vol 12 (1) ◽

Author(s):

Daniel C. Koboldt

Keyword(s):

Best Practices ◽

Single Molecule ◽

Best Practice ◽

Variant Calling ◽

Clinical Samples ◽

Clinical Genetic ◽

Inherited Disorders ◽

Clinical Sequencing ◽

Sequencing Technologies ◽

Downstream Analysis

Abstract Next-generation sequencing technologies have enabled a dramatic expansion of clinical genetic testing both for inherited conditions and diseases such as cancer. Accurate variant calling in NGS data is a critical step upon which virtually all downstream analysis and interpretation processes rely. Just as NGS technologies have evolved considerably over the past 10 years, so too have the software tools and approaches for detecting sequence variants in clinical samples. In this review, I discuss the current best practices for variant calling in clinical sequencing studies, with a particular emphasis on trio sequencing for inherited disorders and somatic mutation detection in cancer patients. I describe the relative strengths and weaknesses of panel, exome, and whole-genome sequencing for variant detection. Recommended tools and strategies for calling variants of different classes are also provided, along with guidance on variant review, validation, and benchmarking to ensure optimal performance. Although NGS technologies are continually evolving, and new capabilities (such as long-read single-molecule sequencing) are emerging, the “best practice” principles in this review should be relevant to clinical variant calling in the long term.

Download Full-text

LRSDAY: Long-read Sequencing Data Analysis for Yeasts

10.1101/184572 ◽

2017 ◽

Author(s):

Jia-Xing Yue ◽

Gianni Liti

Keyword(s):

Genome Assembly ◽

Model Organism ◽

Sequencing Data ◽

Protein Coding ◽

Sequencing Technologies ◽

Long Reads ◽

Long Read ◽

Downstream Analysis ◽

Eukaryotic Organisms ◽

Genomic Regions

AbstractLong-read sequencing technologies have become increasingly popular in genome projects due to their strengths in resolving complex genomic regions. As a leading model organism with small genome size and great biotechnological importance, the budding yeast, Saccharomyces cerevisiae, has many isolates currently being sequenced with long reads. However, analyzing long-read sequencing data to produce high-quality genome assembly and annotation remains challenging. Here we present LRSDAY, the first one-stop solution to streamline this process. LRSDAY can produce chromosome-level end-to-end genome assembly and comprehensive annotations for various genomic features (including centromeres, protein-coding genes, tRNAs, transposable elements and telomere-associated elements) that are ready for downstream analysis. Although tailored for S. cerevisiae, we designed LRSDAY to be highly modular and customizable, making it adaptable for virtually any eukaryotic organisms. Applying LRSDAY to a S. cerevisiae strain takes ∼43 hrs to generate a complete and well-annotated genome from ∼100X Pacific Biosciences (PacBio) reads using four threads.

Download Full-text

VARIFI—Web-Based Automatic Variant Identification, Filtering and Annotation of Amplicon Sequencing Data

Journal of Personalized Medicine ◽

10.3390/jpm9010010 ◽

2019 ◽

Vol 9 (1) ◽

pp. 10

Author(s):

Milica Krunic ◽

Peter Venhuizen ◽

Leonhard Müllauer ◽

Bettina Kaserer ◽

Arndt von Haeseler

Keyword(s):

Genetic Variants ◽

Confidence Score ◽

Amplicon Sequencing ◽

Personal Genome ◽

Nucleotide Polymorphisms ◽

Sequencing Data ◽

Web Based ◽

Sequencing Errors ◽

Sequencing Technologies ◽

Variant Information

Fast and affordable benchtop sequencers are becoming more important in improving personalized medical treatment. Still, distinguishing genetic variants between healthy and diseased individuals from sequencing errors remains a challenge. Here we present VARIFI, a pipeline for finding reliable genetic variants (single nucleotide polymorphisms (SNPs) and insertions and deletions (indels)). We optimized parameters in VARIFI by analyzing more than 170 amplicon-sequenced cancer samples produced on the Personal Genome Machine (PGM). In contrast to existing pipelines, VARIFI combines different analysis methods and, based on their concordance, assigns a confidence score to each identified variant. Furthermore, VARIFI applies variant filters for biases associated with the sequencing technologies (e.g., incorrectly identified homopolymer-associated indels with Ion Torrent). VARIFI automatically extracts variant information from publicly available databases and incorporates methods for variant effect prediction. VARIFI requires little computational experience and no in-house compute power since the analyses are conducted on our server. VARIFI is a web-based tool available at varifi.cibiv.univie.ac.at.

Download Full-text

Deep Investigating the Changes of Gut Microbiome and Its Correlation With the Shifts of Host Serum Metabolome Around Parturition in Sows

Frontiers in Microbiology ◽

10.3389/fmicb.2021.729039 ◽

2021 ◽

Vol 12 ◽

Author(s):

Hao Fu ◽

Maozhang He ◽

Jinyuan Wu ◽

Yunyan Zhou ◽

Shanlin Ke ◽

...

Keyword(s):

16S Rrna ◽

Gut Microbiota ◽

16S Rrna Gene ◽

Gut Microbiome ◽

Bacterial Species ◽

Late Pregnancy ◽

Rrna Gene ◽

Metagenomic Sequencing ◽

Sequencing Data ◽

Functional Capacities

Parturition is a crucial event in the sow reproduction cycle, which accompanies by a series of physiological changes, including sex hormones, metabolism, and immunity. More and more studies have indicated the changes of the gut microbiota from pregnancy to parturition. However, what bacterial species and functional capacities of the gut microbiome are changed around parturition has been largely unknown, and the correlations between the changes of gut bacterial species and host metabolome were also uncovered. In this study, by combining 16S rRNA gene and shotgun metagenomic sequencing data, and the profiles of serum metabolome and fecal short-chain fatty acids (SCFAs), we investigated the changes of gut microbiome, serum metabolite features and fecal SCFAs from late pregnancy (LP) to postpartum (PO) stage. We found the significant changes of gut microbiota from LP to PO stage in both 16S rRNA gene sequencing and metagenomic sequencing analyses. The bacterial species from Lactobacillus, Streptococcus, and Clostridium were enriched at the LP stage, while the species from Bacteroides, Escherichia, and Campylobacter had higher abundances at the PO stage. Functional capacities of the gut microbiome were also significantly changed and associated with the shifts of gut bacteria. Untargeted metabolomic analyses revealed that the metabolite features related to taurine and hypotaurine metabolism, and arginine biosynthesis and metabolism were enriched at the LP stage, and positively associated with those bacterial species enriched at the LP stage, while the metabolite features associated with vitamin B6 and glycerophospholipid metabolism had higher abundances at the PO stage and were positively correlated with the bacteria enriched at the PO stage. Six kinds of SCFAs were measured in feces samples and showed higher concentrations at the LP stage. These results suggested that the changes of gut microbiome from LP to PO stage lead to the shifts of host lipid, amino acids and vitamin metabolism and SCFA production. The results from this study provided new insights for the changes of sow gut microbiome and host metabolism around parturition, and gave new knowledge for guiding the feeding and maternal care of sows from late pregnancy to lactation in the pig industry.

Download Full-text

SoupX removes ambient RNA contamination from droplet based single-cell RNA sequencing data

10.1101/303727 ◽

2018 ◽

Cited By ~ 58

Author(s):

Matthew D Young ◽

Sam Behjati

Keyword(s):

Single Cell ◽

Rna Sequencing ◽

Expression Profiles ◽

List Type ◽

Sequencing Data ◽

Sequencing Technologies ◽

Single Cell Rna Sequencing ◽

Biological Interpretation ◽

Downstream Analysis ◽

Rna Contamination

AbstractBackgroundDroplet based single-cell RNA sequence analyses assume all acquired RNAs are endogenous to cells. However, any cell free RNAs contained within the input solution are also captured by these assays. This sequencing of cell free RNA constitutes a background contamination that confounds the biological interpretation of single-cell transcriptomic data.ResultsWe demonstrate that contamination from this ‘soup’ of cell free RNAs is ubiquitous, with experiment-specific variations in composition and magnitude. We present a method, SoupX, for quantifying the extent of the contamination and estimating ‘background corrected’ cell expression profiles that seamlessly integrate with existing downstream analysis tools. Applying this method to several datasets using multiple droplet sequencing technologies, we demonstrate that its application improves biological interpretation of otherwise misleading data, as well as improving quality control metrics.ConclusionsWe present ‘SoupX’, a tool for removing ambient RNA contamination from droplet based single cell RNA sequencing experiments. This tool has broad applicability and its application can improve the biological utility of existing and future data sets.Key PointsThe signal from droplet based single cell RNA sequencing is ubiquitously contaminated by capture of ambient mRNA.SoupX is a method to quantify the abundance of these ambient mRNAs and remove them.Correcting for ambient mRNA contamination improves biological interpretation.

Download Full-text

Improved microbial community characterization of 16S rRNA via metagenome hybridization capture enrichment

10.1101/2020.12.18.423101 ◽

2020 ◽

Author(s):

Megan Sarah Beaudry ◽

Jincheng Wang ◽

Troy Kieran ◽

Jesse Thomas ◽

Natalia Juliana Bayona-Vasquez ◽

...

Keyword(s):

16S Rrna ◽

Sequence Similarity ◽

Rrna Gene ◽

Metagenomic Sequencing ◽

Sequencing Data ◽

16S Rrna Sequence ◽

Shotgun Metagenomics ◽

Reference Databases ◽

Rrna Sequences ◽

16S Rrna Sequences

Environmental microbial diversity is often investigated from a molecular perspective using 16S ribosomal RNA (rRNA) gene amplicons and shotgun metagenomics. While amplicon methods are fast, low-cost, and have curated reference databases, they can suffer from amplification bias and are limited in genomic scope. In contrast, shotgun metagenomic methods sample more genomic regions with fewer sequence acquisition biases. However, shotgun metagenomic sequencing is much more expensive (even with moderate sequencing depth) and computationally challenging. Here, we develop a set of 16S rRNA sequence capture baits that offer a potential middle ground with the advantages from both approaches for investigating microbial communities. These baits cover the diversity of all 16S rRNA sequences available in the Greengenes (v. 13.5) database, with no sequence having < 80% sequence similarity to at least one bait for all segments of 16S. The use of our baits provide comparable results to 16S amplicon libraries and shotgun metagenomic libraries when assigning taxonomic units from 16S sequences within the metagenomic reads. We demonstrate that 16S rRNA capture baits can be used on a range of microbial samples (i.e., mock communities and rodent fecal samples) to increase the proportion of 16S rRNA sequences (average >400-fold) and decrease analysis time to obtain consistent community assessments. Furthermore, our study reveals that bioinformatic methods used to analyze sequencing data may have a greater influence on estimates of community composition than library preparation method used, likely in part to the extent and curation of the reference databases considered.

Download Full-text

Diet induces reproducible alterations in the mouse and human gut microbiome

10.1101/541797 ◽

2019 ◽

Cited By ~ 1

Author(s):

Jordan E. Bisanz ◽

Vaibhav Upadhyay ◽

Jessie A. Turnbaugh ◽

Kimberly Ly ◽

Peter J. Turnbaugh

Keyword(s):

16S Rrna ◽

16S Rrna Gene ◽

Meta Analysis ◽

Alpha Diversity ◽

Operational Taxonomic Unit ◽

Molecular Testing ◽

Rrna Gene ◽

Metagenomic Sequencing ◽

Sequencing Data ◽

Dietary Interventions

SummaryThe degree to which diet reproducibly alters the human and mouse gut microbiota remains unclear. Here, we focus on the consumption of a high-fat diet (HFD), one of the most frequently studied dietary interventions in mice. We employed a subject-level meta-analysis framework for unbiased collection and analysis of publicly available 16S rRNA gene and metagenomic sequencing data from studies examining HFD in rodent models. In total, we re-analyzed 27 studies, 1101 samples, and 106 million reads mapping to 16S rRNA gene sequences. We report reproducible changes in gut microbial community structure both within and between studies, including a significant increase in the Firmicutes phylum and decrease in the Bacteroidetes phylum; however, reduced alpha diversity is not a consistent feature of HFD. Finer taxonomic analysis revealed that the strongest signal of HFD on microbiota species composition is Lactococcus spp., which we demonstrate is a common dietary contaminant through the molecular testing of dietary ingredients, culturing, microscopy, and germ-free mouse experiments. After in silico removal of Lactococcus spp., we employed machine learning to define a unique operational taxonomic unit (OTU)-based signature capable of predicting the dietary intake of mice and demonstrate that phylogenetic and gene-family transformations of this model are capable of accurately predicting human samples in controlled feeding settings (area under the receiver operator curve = 0.75 and 0.88 respectively). Together, these results demonstrate the utility of microbiome meta-analyses in identifying robust bacterial signals for mechanistic studies and creates a framework for the routine meta-analysis of microbiome studies in preclinical models.

Download Full-text

Yoghurt consumption is associated with transient changes in the composition of the human gut microbiome

10.21203/rs.3.rs-38248/v1 ◽

2020 ◽

Author(s):

Caroline Ivanne Le Roy ◽

Alexander Kurilshikov ◽

Emily Leeming ◽

Alessia Visconti ◽

Ruth Bowyer ◽

...

Keyword(s):

16S Rrna ◽

Gut Microbiota ◽

Visceral Fat ◽

Gut Microbiome ◽

Body Weight Gain ◽

Healthy Eating Index ◽

Rrna Gene ◽

Metagenomic Sequencing ◽

Sequencing Data ◽

Shotgun Metagenomic Sequencing

Abstract Background: Yoghurt contains live bacteria that could contribute via modulation of the gut microbiota to its reported beneficial effects such as reduced body weight gain and lower incidence of type 2 diabetes. To date, the association between yoghurt consumption and the composition of the gut microbiota is underexplored. Here we used clinical variables, metabolomics, 16S rRNA and shotgun metagenomic sequencing data collected on over 1000 predominantly female UK twins to define the link between the gut microbiota and yoghurt-associated health benefits. Results: According to food frequency questionnaires (FFQ), 73% of subjects consumed yoghurt. Consumers presented a healthier diet pattern (healthy eating index: beta = 2.17±0.34; P = 2.72x10 -10 ) and improved metabolic health characterised by reduced visceral fat (beta = -28.18±11.71 g; P = 0.01). According to 16S rRNA gene analyses and whole shotgun metagenomic sequencing approach consistent taxonomic variations were observed with yoghurt consumption. More specifically, we identified higher abundance of species used as yoghurt starters Streptococcus thermophilus (beta = 0.41±0.051; P = 6.14x10 -12 ) and sometimes added Bifidobacterium animalis subsp. lactis (beta = 0.30±0.052; P = 1.49x10 -8 ) in the gut of yoghurt consumers. Replication in 1103 volunteers from the LL-DEEP cohort confirmed the increase of S. thermophilus among yoghurt consumers. Using food records collected the day prior to faecal sampling we showed that increase in these two yoghurt bacteria could be transient. Metabolomics analysis revealed that B. animalis subsp. lactis was associated with 13 faecal metabolites including a 3-hydroxyoctanoic acid, known to be involved in the regulation of gut inflammation. Conclusions: Yoghurt consumption is associated with reduced visceral fat mass and changes in gut microbiome including transient increase of yoghurt-contained species ( i.e. S. thermophilus and B. lactis ).

Download Full-text