Abstract:Shotgun metagenomics allows comprehensive sampling of the genomic information of microbes in a given environment and is a tool of choice for studying complex microbial systems. Mapping sequencing reads against a set of reference or metagenome-assembled genomes is in principle a simple and powerful approach to define the species-level composition of the microbial community under investigation. However, despite the widespread use of this approach, there is no established way to properly interpret the alignment results, with arbitrary relative abundance thresholds being routinely used to discriminate between present and absent species. Such an approach can be affected by significant biases, especially in the identification of rare species. Therefore, it is important to develop new metrics to overcome these biases. Here, we present Metapresence, a new tool to perform reliable identification of the species in metagenomic samples based on the distribution of mapped reads on the reference genomes. The analysis is based on two metrics describing the breadth of coverage and the genomic distance between consecutive reads. We demonstrate the high precision and wide applicability of the tool using data from various synthetic communities, a real mock community, and the gut microbiome of healthy individuals and antibiotic-associated-diarrhea patients. Overall, our results suggest that the proposed approach has a robust performance in hard-to-analyze microbial communities containing contaminated or closely related genomes in low abundance.IMPORTANCEDespite the prevalent use of genome-centric alignment-based methods to characterize microbial community composition, there lacks a standardized approach for accurately identifying the species within a sample. Currently, arbitrary relative abundance thresholds are commonly employed for this purpose. However, due to the inherent complexity of genome structure and biases associated with genome-centric approaches, this practice tends to be imprecise. Notably, it introduces significant biases, particularly in the identification of rare species. The method presented here addresses these limitations and contributes significantly to overcoming inaccuracies in precisely defining community composition, especially when dealing with rare members.

Estimating the total genome length of a metagenomic sample using k-mers

Sequence Repetitiveness Quantification and De Novo Repeat Detection by Weighted K-Mer Coverage.

Estimation of genomic characteristics by analyzing k-mer frequency in de novo genome projects

Accurate Genome Relative Abundance Estimation Based on Shotgun Metagenomic Reads

Estimation of genome size using k-mer frequencies from corrected long reads

Large-scale estimation of bacterial and archaeal DNA prevalence in metagenomes reveals biome-specific patterns

Towards complete representation of bacterial contents in metagenomic samples

Efficient De Novo Assembly and Recovery of Microbial Genomes from Complex Metagenomes Using a Reduced Set of k-mers

Diminishing Return for Increased Mappability with Longer Sequencing Reads: Implications of the k-mer Distributions in the Human Genome

k-mer-based approaches to bridging pangenomics and population genetics

Multiple Comparative Metagenomics using Multiset k-mer Counting

KMCP: accurate metagenomic profiling of both prokaryotic and viral populations by pseudo-mapping

Comprehensive Simulation of Metagenomic Sequencing Data with Non-Uniform Sampling Distribution

Large-scale k-mer-based analysis of the informational properties of genomes, comparative genomics and taxonomy

A high-precision genome size estimator based on the k-mer histogram correction

Multi-sample Estimation of Bacterial Composition Matrix in Metagenomics Data

Estimating the Length Distributions of Genomic Micro-satellites from Next Generation Sequencing Data.

Assessment of Kmer Degeneration Method for Complicated Genomes.

Metapresence: a tool for accurate species detection in metagenomics based on the genome-wide distribution of mapping reads

How Much Metagenomic Sequencing is Enough to Achieve a Given Goal?

Mapping-based genome size estimation