Abstract:Comparing metagenomic samples is crucial for understanding microbial communities. For different groups of microbial communities, such as human gut metagenomic samples from patients with a certain disease and healthy controls, identifying group-specific sequences offers essential information for potential biomarker discovery. A sequence that is present, or rich, in one group, but absent, or scarce, in another group is considered "group-specific" in our study. Our main purpose is to discover group-specific sequence regions between control and case groups as disease-associated markers. We developed a long k-mer (k ≥ 30 bps)-based computational pipeline to detect group-specific sequences at strain resolution free from reference sequences, sequence alignments, and metagenome-wide de novo assembly. We called our method MetaGO: Group-specific oligonucleotide analysis for metagenomic samples. An open-source pipeline on Apache Spark was developed with parallel computing. We applied MetaGO to one simulated and three real metagenomic datasets to evaluate the discriminative capability of identified group-specific markers. In the simulated dataset, 99.11% of group-specific logical 40-mers covered 98.89% disease-specific regions from the disease-associated strain. In addition, 97.90% of group-specific numerical 40-mers covered 99.61 and 96.39% of differentially abundant genome and regions between two groups, respectively. For a large-scale metagenomic liver cirrhosis (LC)-associated dataset, we identified 37,647 group-specific 40-mer features. Any one of the features can predict disease status of the training samples with the average of sensitivity and specificity higher than 0.8. The random forests classification using the top 10 group-specific features yielded a higher AUC (from ∼0.8 to ∼0.9) than that of previous studies. All group-specific 40-mers were present in LC patients, but not healthy controls. All the assembled 11 LC-specific sequences can be mapped to two strains of Veillonella parvula: UTDB1-3 and DSM2008. The experiments on the other two real datasets related to Inflammatory Bowel Disease and Type 2 Diabetes in Women consistently demonstrated that MetaGO achieved better prediction accuracy with fewer features compared to previous studies. The experiments showed that MetaGO is a powerful tool for identifying group-specific k-mers, which would be clinically applicable for disease prediction. MetaGO is available at https://github.com/VVsmileyx/MetaGO.

MetaScope - Fast and accurate identification of microbes in metagenomic sequencing data

Fast and Simple Analysis of MiSeq Amplicon Sequencing Data with MetaAmp.

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

PathoScope 2.0: a complete computational framework for strain identification in environmental or clinical sequencing samples

MSC: a metagenomic sequence classification algorithm

Metaxa: a software tool for automated detection and discrimination among ribosomal small subunit (12S/16S/18S) sequences of archaea, bacteria, eukaryotes, mitochondria, and chloroplasts in metagenomes and environmental sequencing datasets

MTSv: rapid alignment-based taxonomic classification and high-confidence metagenomic analysis

SMeta, a Binning Tool Using Single-Cell Sequences to Aid Reconstructing Metageome Species Accurately

Clinical PathoScope: rapid alignment and filtration for accurate pathogen identification in clinical samples using unassembled sequencing data

Parallel-META 2.0: Enhanced Metagenomic Data Analysis with Functional Annotation, High Performance Computing and Advanced Visualization

Parallel-META 3: Comprehensive Taxonomical and Functional Analysis Platform for Efficient Comparison of Microbial Communities

A big data approach to metagenomics for all-food-sequencing

SMeta, a binning tool using single-cell sequences to aid in reconstructing species from metagenome accurately

Enhancing antimicrobial resistance detection with MetaGeneMiner: Targeted gene extraction from metagenomes

Accurate Reconstruction of Microbial Strains from Metagenomic Sequencing Using Representative Reference Genomes

A multi-amplicon 16S rRNA sequencing and analysis method for improved taxonomic profiling of bacterial communities

microbeMASST: a taxonomically informed mass spectrometry search tool for microbial metabolomics data

TAMA: improved metagenomic sequence classification through meta-analysis

MGS-Fast: Metagenomic shotgun data fast annotation using microbial gene catalogs

MetaComp: Comprehensive Analysis Software for Comparative Meta-Omics Including Comparative Metagenomics

Identifying Group-Specific Sequences for Microbial Communities Using Long K-Mer Sequence Signatures.