Wei Shen 沈伟 Profile picture
Mar 9 13 tweets 5 min read
Can't wait to share the KMCP preprint with you, hope it's not too late. Thanks for reading and any feedback is welcome and appreciated. biorxiv.org/content/10.110…
We present a novel metagenomic profiling tool, KMCP, which not only allows for accurate taxonomic profiling of archaea, bacteria, and viral populations from metagenomic shotgun sequence data, but also provides confident pathogen detection for clinical samples of low depth.
We reimplement a modified COBS data structure to index reference genomes. The index size of KMCP is smaller than COBS, and the batch searching speed is increased by near 10 times. And the database construction is fast, taking only 25 min to index 47,894 genomes from GTDB.
To reduce the high false positives in traditional k-mer based methods, we introduce genomic positions to k-mers by splitting the reference genomes into chunks. We track both the origin and the approximate location in the reference for each query, denoted as pseudo-mapping.
KMCP utilizes genome coverage (chunks faction) to filter out spurious references, i.e., the reference genome chunks need to be uniformly covered by reads. Meanwhile, we require a matched genome having some high-confidence uniquely matched reads with high k-mer coverages.
Benchmarks on CAMI2 mouse gut data and Sun's simulated datasets showed KMCP has high accuracies close to MetaPhlAn3 and mOTUs3. Though mOTUs3 has good accuracy on prokaryotic metagenome datasets, it lacks the ability of viral detection.
Benchmarks on a mock virome dataset showed that KMCP outperformed Centrifuge, Bracken, and MetaPhlAn3, with high accuracies of both taxon identification and abundance estimation.
In pathogen detection, KMCP had an accuracy close to Kraken2, while generating shorter prediction lists and generating pathogens in priority positions by considering both sequence similarity and the genome chunks fraction, which help researchers rapidly interpret the reports.
We provide 3 pre-built databases for metagenomic profiling, including the prokaryotic database created with 47,894 archaea and bacteria from GTDB, the viral database with 27,936 virus genomes from Genbank, and the fungal database with 403 fungi genomes from RefSeq.
With the same reference genomes, KMCP generates much smaller databases (66 GB) than Bracken (304 GB) and Centrifuge (97 GB), and requires less (<60 GB) memory for searching than Kraken2/Bracken (250 GB) and Centrifuge (97 GB).
The mergeability of search results makes it flexible to build databases for different reference datasets and choose various databases to search. To update databases, users can either re-built the database with newly added genomes or only build an additional DB for the new genomes
KMCP searching speed is much slower than Bracken, Centrifuge, MetaPhlAn3, and mOTUs3. Fortunately, KMCP could utilize computer clusters to linearly accelerate the reads searching, with each node searching against a small database built with a partition of the reference genomes.
KMCP is implemented in #golang, with a single executable binary file available, and the source code is freely available under the MIT License. The latest version of KMCP can be obtained from the Bioconda channel of Conda package manager and github.com/shenwei356/kmcp

• • •

Missing some Tweet in this thread? You can try to force a refresh
 

Keep Current with Wei Shen 沈伟

Wei Shen 沈伟 Profile picture

Stay in touch and get notified when new unrolls are available from this author!

Read all threads

This Thread may be Removed Anytime!

PDF

Twitter may remove this content at anytime! Save it as PDF for later use!

Try unrolling a thread yourself!

how to unroll video
  1. Follow @ThreadReaderApp to mention us!

  2. From a Twitter thread mention us with a keyword "unroll"
@threadreaderapp unroll

Practice here first or read more on our help page!

Did Thread Reader help you today?

Support us! We are indie developers!


This site is made by just two indie developers on a laptop doing marketing, support and development! Read more about the story.

Become a Premium Member ($3/month or $30/year) and get exclusive features!

Become Premium

Don't want to be a Premium member but still want to support us?

Make a small donation by buying us coffee ($5) or help with server cost ($10)

Donate via Paypal

Or Donate anonymously using crypto!

Ethereum

0xfe58350B80634f60Fa6Dc149a72b4DFbc17D341E copy

Bitcoin

3ATGMxNzCUFzxpMCHL5sWSt4DVtS8UqXpi copy

Thank you for your support!

Follow Us on Twitter!

:(