Strain identification and quantitative analysis in microbial communities

Genome-resolved metagenomics: a game changer for microbiome medicine

Recent substantial evidence implicating commensal bacteria in human diseases has given rise to a new domain in biomedical research: microbiome medicine. This emerging field aims to understand and leverage the human microbiota and derivative molecules for disease prevention and treatment. Despite the complex and hierarchical organization of this ecosystem, most research over the years has relied on 16S amplicon sequencing, a legacy of bacterial phylogeny and taxonomy. Although advanced sequencing technologies have enabled cost-effective analysis of entire microbiota, translating the relatively short nucleotide information into the functional and taxonomic organization of the microbiome has posed challenges until recently. In the last decade, genome-resolved metagenomics, which aims to reconstruct microbial genomes directly from whole-metagenome sequencing data, has made significant strides and continues to unveil the mysteries of various human-associated microbial communities. There has been a rapid increase in the volume of whole metagenome sequencing data and in the compilation of novel metagenome-assembled genomes and protein sequences in public depositories. This review provides an overview of the capabilities and methods of genome-resolved metagenomics for studying the human microbiome, with a focus on investigating the prokaryotic microbiota of the human gut. Just as decoding the human genome and its variations marked the beginning of the genomic medicine era, unraveling the genomes of commensal microbes and their sequence variations is ushering us into the era of microbiome medicine. Genome-resolved metagenomics stands as a pivotal tool in this transition and can accelerate our journey toward achieving these scientific and medical milestones.

Sequencing-based analysis of microbiomes

Microbiomes occupy a range of niches and, in addition to having diverse compositions, they have varied functional roles that have an impact on agriculture, environmental sciences, and human health and disease. The study of microbiomes has been facilitated by recent technological and analytical advances, such as cheaper and higher-throughput DNA and RNA sequencing, improved long-read sequencing and innovative computational analysis methods. These advances are providing a deeper understanding of microbiomes at the genomic, transcriptional and translational level, generating insights into their function and composition at resolutions beyond the species level.

Impact of a Purified Microbiome Therapeutic on Abundance of Antimicrobial Resistance Genes in Patients With Recurrent Clostridioides difficile Infection

BACKGROUND

The gastrointestinal microbiota is an important line of defense against colonization with antimicrobial resistant (AR) bacteria. In this post hoc analysis of the phase 3 ECOSPOR III trial, we assessed impact of a microbiota-based oral therapeutic (fecal microbiota spores, live; VOWST Oral Spores [VOS], formerly SER-109]; Seres Therapeutics) compared with placebo, on AR gene (ARG) abundance in patients with recurrent Clostridioides difficile infection (rCDI).

METHODS

Adults with rCDI were randomized to receive VOS or placebo orally for 3 days following standard-of-care antibiotics. ARG and taxonomic profiles were generated using whole metagenomic sequencing of stool at baseline and weeks 1, 2, 8, and 24 posttreatment.

RESULTS

Baseline (n = 151) and serial posttreatment stool samples collected through 24 weeks (total N = 472) from 182 patients (59.9% female; mean age: 65.5 years) in ECOSPOR III as well as 68 stool samples obtained at a single time point from a healthy cohort were analyzed. Baseline ARG abundance was similar between arms and significantly elevated versus the healthy cohort. By week 1, there was a greater decline in ARG abundance in VOS versus placebo (P = .003) in association with marked decline of Proteobacteria and repletion of spore-forming Firmicutes, as compared with baseline. We observed abundance of Proteobacteria and non-spore-forming Firmicutes were associated with ARG abundance, while spore-forming Firmicutes abundance was negatively associated.

CONCLUSIONS

This proof-of-concept analysis suggests that microbiome remodeling with Firmicutes spores may be a potential novel approach to reduce ARG colonization in the gastrointestinal tract.

CyanoStrainChip: A Novel DNA Microarray Tool for High-Throughput Detection of Environmental Cyanobacteria at the Strain Level

Detecting cyanobacteria in environments is an important concern due to their crucial roles in ecosystems, and they can form blooms with the potential to harm humans and nonhuman entities. However, the most widely used methods for high-throughput detection of environmental cyanobacteria, such as 16S rRNA sequencing, typically provide above-species-level resolution, thereby disregarding intraspecific variation. To address this, we developed a novel DNA microarray tool, termed the CyanoStrainChip, that enables strain-level comprehensive profiling of environmental cyanobacteria. The CyanoStrainChip was designed to target 1277 strains; nearly all major groups of cyanobacteria are included by implementing 43,666 genome-wide, strain-specific probes. It demonstrated strong specificity by in vitro mock community experiments. The high correlation (Pearson's R > 0.97) between probe fluorescence intensities and the corresponding DNA amounts (ranging from 1-100 ng) indicated excellent quantitative capability. Consistent cyanobacterial profiles of field samples were observed by both the CyanoStrainChip and next-generation sequencing methods. Furthermore, CyanoStrainChip analysis of surface water samples in Lake Chaohu uncovered a high intraspecific variation of abundance change within the genus Microcystis between different severity levels of cyanobacterial blooms, highlighting two toxic Microcystis strains that are of critical concern for Lake Chaohu harmful blooms suppression. Overall, these results suggest a potential for CyanoStrainChip as a valuable tool for cyanobacterial ecological research and harmful bloom monitoring to supplement existing techniques.

Probiotic and postbiotic analytical methods: a perspective of available enumeration techniques

Probiotics are the largest non-herbal/traditional dietary supplements category worldwide. To be effective, a probiotic strain must be delivered viable at an adequate dose proven to deliver a health benefit. The objective of this article is to provide an overview of the various technologies available for probiotic enumeration, including a general description of each technology, their advantages and limitations, and their potential for the future of the probiotics industry. The current "gold standard" for analytical quantification of probiotics in the probiotic industry is the Plate Count method (PC). PC measures the bacterial cell's ability to proliferate into detectable colonies, thus PC relies on cultivability as a measure of viability. Although viability has widely been measured by cultivability, there has been agreement that the definition of viability is not limited to cultivability. For example, bacterial cells may exist in a state known as viable but not culturable (VBNC) where the cells lose cultivability but can maintain some of the characteristics of viable cells as well as probiotic properties. This led to questioning the association between viability and cultivability and the accuracy of PC in enumerating all the viable cells in probiotic products. PC has always been an estimate of the number of viable cells and not a true cell count. Additionally, newer probiotic categories such as Next Generation Probiotics (NGPs) are difficult to culture in routine laboratories as NGPs are often strict anaerobes with extreme sensitivity to atmospheric oxygen. Thus, accurate quantification using culture-based techniques will be complicated. Another emerging category of biotics is postbiotics, which are inanimate microorganisms, also often referred to as tyndallized or heat-killed bacteria. Obviously, culture dependent methods are not suitable for these products, and alternative methods are needed for their quantification. Different methodologies provide a more complete picture of a heterogeneous bacterial population versus PC focusing exclusively on the eventual multiplication of the cells. Alternative culture-independent techniques including real-time PCR, digital PCR and flow cytometry are discussed. These methods can measure viability beyond cultivability (i.e., by measuring cellular enzymatic activity, membrane integrity or membrane potential), and depending on how they are designed they can achieve strain-specific enumeration.

Maast: genotyping thousands of microbial strains efficiently

Existing single nucleotide polymorphism (SNP) genotyping algorithms do not scale for species with thousands of sequenced strains, nor do they account for conspecific redundancy. Here we present a bioinformatics tool, Maast, which empowers population genetic meta-analysis of microbes at an unrivaled scale. Maast implements a novel algorithm to heuristically identify a minimal set of diverse conspecific genomes, then constructs a reliable SNP panel for each species, and enables rapid and accurate genotyping using a hybrid of whole-genome alignment and k-mer exact matching. We demonstrate Maast's utility by genotyping thousands of Helicobacter pylori strains and tracking SARS-CoV-2 diversification.

Identifying species-specific k-mers for fast and accurate metagenotyping with Maast and GT-Pro

Genotyping single-nucleotide polymorphisms (SNPs) in microbiomes enables strain-level quantification. In this protocol, we describe a computational pipeline that performs fast and accurate SNP genotyping using metagenomic data. We first demonstrate how to use Maast to catalog SNPs from microbial genomes. Then we use GT-Pro to extract unique SNP-covering k-mers, optimize a data structure for storing these k-mers, and finally perform metagenotyping. For proof of concept, the protocol leverages public whole-genome sequences to metagenotype a synthetic community. For complete details on the use and execution of this protocol, please refer to Shi et al. (2022a)1 and Shi et al. (2022b).2.

Eco-friendly bioremediation of pollutants from contaminated sewage wastewater using special reference bacterial strain of Bacillus cereus SDN1 and their genotoxicological assessment in Allium cepa

The present study aimed to assess the Bacillus cereus SDN1 native bacterium's ability to clean up contaminated or polluted water. The isolated bacterium was identified by its morphological and biochemical characteristics, which were then confirmed at the genus level. Furthermore, the isolated B. cereus (NCBI accession No: MW828583) was identified genomically by PCR amplifying 16 s rDNA using a universal primer. The phylogenetic analysis of the rDNA sequence was analyzed to determine the taxonomic and evolutionary profile of the isolate of the previously identified Bacillus sp. Besides, B. cereus and the bacterial consortium were treated using sewage wastewater. After 15 days of treatment, the following pollutants or chemicals were reduced: total hardness particles removal varied from 63.33 % to 67.55 %, calcium removal varied from 90 % to 93.33 %, and total nitrate decreased range from 37.77 % to 22.22 %, respectively. Electrical conductivity ranged from 1809 mS/cm to 2500 mS/cm, and pH values ranged from 6.5 to 8.95. The outcome of in-situ remediation results suggested that B. cereus has a noticeable remediation efficiency to the suspended particles. A root tip test was also used to investigate the genotoxicity of treated and untreated sewage-contaminated waters on onion (Allium cepa) root cells. The highest chromosomal aberrations and mitotic inhibition were found in roots exposed to contaminated sewage water, and their results displayed chromosome abnormalities, including disorganized, sticky chain, disturbed metaphase, chromosomal displacement in anaphase, abnormal telophase, spindle disturbances, and binucleate cells observed in A. cepa exposed to untreated contaminated water. The study can thus be applied as a biomarker to detect the genotoxic impacts of sewage water pollution on biota. Furthermore, based on an identified bacterial consortium, this work offers a low-cost and eco-favorable method for treating household effluents.

Pitfalls of genotyping microbial communities with rapidly growing genome collections

Detecting genetic variants in metagenomic data is a priority for understanding the evolution, ecology, and functional characteristics of microbial communities. Many tools that perform this metagenotyping rely on aligning reads of unknown origin to a database of sequences from many species before calling variants. In this synthesis, we investigate how databases of increasingly diverse and closely related species have pushed the limits of current alignment algorithms, thereby degrading the performance of metagenotyping tools. We identify multi-mapping reads as a prevalent source of errors and illustrate a trade-off between retaining correct alignments versus limiting incorrect alignments, many of which map reads to the wrong species. Then we evaluate several actionable mitigation strategies and review emerging methods showing promise to further improve metagenotyping in response to the rapid growth in genome collections. Our results have implications beyond metagenotyping to the many tools in microbial genomics that depend upon accurate read mapping.

Statistical and Computational Methods for Microbial Strain Analysis

Microbial strains are interpreted as a lineage derived from a recent ancestor that have not experienced "too many" recombination events and can be successfully retrieved with culture-independent techniques using metagenomic sequencing. Such a strain variability has been increasingly shown to display additional phenotypic heterogeneities that affect host health, such as virulence, transmissibility, and antibiotics resistance. New statistical and computational methods have recently been developed to track the strains in samples based on shotgun metagenomics data either based on reference genome sequences or Metagenome-assembled genomes (MAGs). In this paper, we review some recent statistical methods for strain identifications based on frequency counts at a set of single nucleotide variants (SNVs) within a set of single-copy marker genes. These methods differ in terms of whether reference genome sequences are needed, how SNVs are called, what methods of deconvolution are used and whether the methods can be applied to multiple samples. We conclude our review with areas that require further research.

Genotyping Microbial Communities with MIDAS2: From Metagenomic Reads to Allele Tables

The Metagenomic Intra-Species Diversity Analysis System 2 (MIDAS2) is a scalable pipeline that identifies single nucleotide variants and gene copy number variants in metagenomes using comprehensive reference databases built from public microbial genome collections (metagenotyping). MIDAS2 is the first metagenotyping tool with functionality to control metagenomic read mapping filters and to customize the reference database to the microbial community, features that improve the precision and recall of detected variants. In this article we present four basic protocols for the most common use cases of MIDAS2, along with supporting protocols for installation and use. In addition, we provide in-depth guidance on adjusting command line parameters, editing the reference database, optimizing hardware utilization, and understanding the metagenotyping results. All the steps of metagenotyping, from raw sequencing reads to population genetic analysis, are demonstrated with example data in two downloadable sequencing libraries of single-end metagenomic reads representing a mixture of multiple bacterial species. This set of protocols empowers users to accurately genotype hundreds of species in thousands of samples, providing rich genetic data for studying the evolution and strain-level ecology of microbial communities. © 2022 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Species prescreening Basic Protocol 2: Download MIDAS reference database Basic Protocol 3: Population single nucleotide variant calling Basic Protocol 4: Pan-genome copy number variant calling Support Protocol 1: Installing MIDAS2 Support Protocol 2: Command line inputs Support Protocol 3: Metagenotyping with a custom collection of genomes Support Protocol 4: Metagenotyping with advanced parameters.