Abstract:With the development of sequencing technology and analytic tools, studying within-species variations enhances the understanding of microbial biological processes. Nevertheless, most existing methods designed for strain-level analysis lack the capability to concurrently assess both strain proportions and genome-wide single nucleotide variants (SNVs) across longitudinal metagenomic samples. In this study, we introduce LongStrain, an integrated pipeline for the analysis of large-scale metagenomic data from individuals with longitudinal or repeated samples. In LongStrain, we first utilize two efficient tools, Kraken2 and Bowtie2, for the taxonomic classification and alignment of sequencing reads, respectively. Subsequently, we propose to jointly model strain proportions and shared haplotypes across samples within individuals. This approach specifically targets tracking a primary strain and a secondary strain for each subject, providing their respective proportions and SNVs as output. With extensive simulation studies of a microbial community and single species, our results demonstrate that LongStrain is superior to two genotyping methods and two deconvolution methods across a majority of scenarios. Furthermore, we illustrate the potential applications of LongStrain in the real data analysis of The Environmental Determinants of Diabetes in the Young study and a gastric intestinal metaplasia microbiome study. In summary, the proposed analytic pipeline demonstrates marked statistical efficiency over the same type of methods and has great potential in understanding the genomic variants and dynamic changes at strain level. LongStrain and its tutorial are freely available online at https://github.com/BoyanZhou/LongStrain. Importance: The advancement in DNA-sequencing technology has enabled the high-resolution identification of microorganisms in microbial communities. Since different microbial strains within species may contain extreme phenotypic variability (e.g., nutrition metabolism, antibiotic resistance, and pathogen virulence), investigating within-species variations holds great scientific promise in understanding the underlying mechanism of microbial biological processes. To fully utilize the shared genomic variants across longitudinal metagenomics samples collected in microbiome studies, we develop an integrated analytic pipeline (LongStrain) for longitudinal metagenomics data. It concurrently leverages the information on proportions of mapped reads for individual strains and genome-wide SNVs to enhance the efficiency and accuracy of strain identification. Our method helps to understand strains' dynamic changes and their association with genome-wide variants. Given the fast-growing longitudinal studies of microbial communities, LongStrain which streamlines analyses of large-scale raw sequencing data should be of great value in microbiome research communities.

Strain/species identification in metagenomes using genome-specific markers

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

Strain-GeMS: Optimized Subspecies Identification from Microbiome Data Based on Accurate Variant Modeling.

microBiomeGSM: the identification of taxonomic biomarkers from metagenomic data using grouping, scoring and modeling (G-S-M) approach

Strain/Species-Specific Probe Design for Microbial Identification Microarrays

PStrain: an iterative microbial strains profiling algorithm for shotgun metagenomic sequencing data

Improved Metagenomic Taxonomic Profiling Using a Curated Core Gene-Based Bacterial Database Reveals Unrecognized Species in the Genus Streptococcus

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

Mining the unknown: a systems approach to metabolite identification combining genetic and metabolic information.

Toward Accurate Diagnosis and Surveillance of Bacterial Infections Using Enhanced Strain-Level Metagenomic Next-Generation Sequencing of Infected Body Fluids.

Identifying individual-specific microbial DNA fingerprints from skin microbiomes

Scalable Microbial Strain Inference in Metagenomic Data Using StrainFacts

Species-Specific Identification of Streptococcus Based on DNA Marker in 16S–23S Rdna Internal Transcribed Spacer

An integrated strain-level analytic pipeline utilizing longitudinal metagenomic data

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

Quantitative Detection of Bifidobacterium longum Strains in Feces Using Strain-Specific Primers

Strain-level Dissection of the Contribution of the Gut Microbiome to Human Metabolic Disease

Markers for Discriminating Campylobacter Concisus Genomospecies Using MALDI-TOF Analysis

Quantification Of Microbial Species In Solid State Fermentation Samples Using Signature Genomic Sequences

Cultivation of Common Bacterial Species and Strains from Human Skin, Oral, and Gut Microbiota

Microbial strain-level population structure and genetic diversity from metagenomes