Reduced metagenome sequencing for strain-resolution taxonomic profiles

Abstract

Background Studies of shifts in microbial community composition has many applications. For studies at species or subspecies levels, the 16S amplicon sequencing lacks resolution and is often replaced by full shotgun sequencing. Due to higher costs, this restricts the number of samples sequenced. As an alternative to a full shotgun sequencing we have investigated the use of Reduced Metagenome Sequencing (RMS) to estimate the composition of a microbial community. This involves the use of double-digested restriction-associated DNA sequencing, which means only a smaller fraction of the genomes are sequenced. The read sets obtained by this approach have properties different from both amplicon and shotgun data, and analysis pipelines for both can either not be used at all or not explore the full potential of RMS data. Results We suggest a procedure for analyzing such data, based on fragment clustering and the use of a constrained ordinary least square de-convolution for estimating the relative abundance of all community members. Mock community datasets show the potential to clearly separate strains even when the 16S is 100% identical, and genome-wide differences is < 0.02, indicating RMS has a very high resolution. From a simulation study, we compare RMS to shotgun sequencing and show that we get improved abundance estimates when the community has many very closely related genomes. From a real dataset of infant guts, we show that RMS is capable of detecting a strain diversity gradient for Escherichia coli across time. Conclusion We find that RMS is a good alternative to either metabarcoding or shotgun sequencing when it comes to resolving microbial communities at the strain level. Like shotgun metagenomics, it requires a good database of reference genomes and is well suited for studies of the human gut or other communities where many reference genomes exist. A data analysis pipeline is offered, as an R package at https://github.com/larssnip/microRMS.

Exploration of the presence and abundance of multidrug resistance efflux genes in oil and gas environments

“…This stems from requiring knowledge of the protein sequence to accurately predict functionality, but our ability to predict folding is incomplete as we continue to find dissimilar domains with similar functions and vice versa as well as moonlighting functionalities [63]. Recently, methods have been developed to reduce the number of genomes sequenced in a metagenome for improved quality, but thse have not been widely used [64]. To accommodate the sequence diversity, functional annotation has been based on predicted ORFs and approximate start and stop codons, which leads to millions of hypothetical genes in rich metagenomes [65].…”

A Statistical Framework for Data Purification with Application to Microbiome Data Analysis

“…Our current study focuses on the disease-related microbial signatures by comparing different cohort populations. This may lead to spurious associations with diseases since different taxonomic compositions of microbes may contribute to similar functionality [25]. With the available metagenome and metatranscriptome data [26, 27], one can give a read-out of functional activities within microbial communities, not limited to taxa, thereby studying the disease-associated taxonomic and functional mechanisms simultaneously.…”

Readsynth: short-read simulation for consideration of composition-biases in reduced metagenome sequencing approaches

“…These communities used standardized concentrations of each community member, making them ideal for comparison with simulated reads. The first dataset from Snipen et al consists of a Human Microbiome Project mock community of 20 bacterial strains (BEI HM-782D) digested using the restriction enzymes EcoRI and MseI [ 5 ]. To simulate the abundance of each taxon, the input abundance table used the reciprocal of the ribosomal copy number determined from the ribosomal RNA operon copy number database [ 16 ].…”