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.

InStrain enables population genomic analysis from metagenomic data and rigorous detection of identical microbial strains

Microbial strain-level population structure and genetic diversity from metagenomes

An integrated strain-level analytic pipeline utilizing longitudinal metagenomic data

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

Scalable Microbial Strain Inference in Metagenomic Data Using StrainFacts

Strain tracking in complex microbiomes using synteny analysis reveals per-species modes of evolution

Metagenomic strain detection with SameStr: identification of a persisting core gut microbiota transferable by fecal transplantation

StrainPanDA: Linked reconstruction of strain composition and gene content profiles via pangenome-based decomposition of metagenomic data

Accurate estimation of intraspecific microbial gene content variation in metagenomic data with MIDAS v3 and StrainPGC

Strainy: phasing and assembly of strain haplotypes from long-read metagenome sequencing

StrainSelect: A novel microbiome reference database that disambiguates all bacterial strains, genome assemblies and extant cultures worldwide

Towards Strain-Level Complexity: Sequencing Depth Required for Comprehensive Single-Nucleotide Polymorphism Analysis of the Human Gut Microbiome.

Assessing the performance of current strain resolution tools on long-read metagenomes

Average Nucleotide Identity based strain grouping allows identification of strain-specific genes in the pangenome

Strains, functions and dynamics in the expanded Human Microbiome Project

One species, many faces: The underappreciated importance of strain diversity

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

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

Intracellular spatial transcriptomic analysis toolkit (InSTAnT)

Towards estimating the number of strains that make up a natural bacterial population

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