Reference-free structural variant detection in microbiomes via long-read co-assembly graphs.

MOTIVATION: The study of bacterial genome dynamics is vital for understanding the mechanisms underlying microbial adaptation, growth, and their impact on host phenotype. Structural variants (SVs), genomic alterations of 50 base pairs or more, play a pivotal role in driving evolutionary processes and maintaining genomic heterogeneity within bacterial populations. While SV detection in isolate genomes is relatively straightforward, metagenomes present broader challenges due to the absence of clear reference genomes and the presence of mixed strains. In response, our proposed method rhea, forgoes reference genomes and metagenome-assembled genomes (MAGs) by encompassing all metagenomic samples in a series (time or other metric) into a single co-assembly graph. The log fold change in graph coverage between successive samples is then calculated to call SVs that are thriving or declining. RESULTS: We show rhea to outperform existing methods for SV and horizontal gene transfer (HGT) detection in two simulated mock metagenomes, particularly as the simulated reads diverge from reference genomes and an increase in strain diversity is incorporated. We additionally demonstrate use cases for rhea on series metagenomic data of environmental and fermented food microbiomes to detect specific sequence alterations between successive time and temperature samples, suggesting host advantage. Our approach leverages previous work in assembly graph structural and coverage patterns to provide versatility in studying SVs across diverse and poorly characterized microbial communities for more comprehensive insights into microbial gene flux. AVAILABILITY AND IMPLEMENTATION: rhea is open source and available at: https://github.com/treangenlab/rhea.

Construction and characterization of a genome-scale ordered mutant collection of Bacteroides thetaiotaomicron.

BACKGROUND: Ordered transposon-insertion collections, in which specific transposon-insertion mutants are stored as monocultures in a genome-scale collection, represent a promising tool for genetic dissection of human gut microbiota members. However, publicly available collections are scarce and the construction methodology remains in early stages of development. RESULTS: Here, we describe the assembly of a genome-scale ordered collection of transposon-insertion mutants in the model gut anaerobe Bacteroides thetaiotaomicron VPI-5482 that we created as a resource for the research community. We used flow cytometry to sort single cells from a pooled library, located mutants within this initial progenitor collection by applying a pooling strategy with barcode sequencing, and re-arrayed specific mutants to create a condensed collection with single-insertion strains covering >2500 genes. To demonstrate the potential of the condensed collection for phenotypic screening, we analyzed growth dynamics and cell morphology. We identified both growth defects and altered cell shape in mutants disrupting sphingolipid synthesis and thiamine scavenging. Finally, we analyzed the process of assembling the B. theta condensed collection to identify inefficiencies that limited coverage. We demonstrate as part of this analysis that the process of assembling an ordered collection can be accurately modeled using barcode sequencing data. CONCLUSION: We expect that utilization of this ordered collection will accelerate research into B. theta physiology and that lessons learned while assembling the collection will inform future efforts to assemble ordered mutant collections for an increasing number of gut microbiota members.

Hydrogenotrophic methanogenesis in archaeal phylum Verstraetearchaeota reveals the shared ancestry of all methanogens.

Methanogenic archaea are major contributors to the global carbon cycle and were long thought to belong exclusively to the euryarchaeal phylum. Discovery of the methanogenesis gene cluster methyl-coenzyme M reductase (Mcr) in the Bathyarchaeota, and thereafter the Verstraetearchaeota, led to a paradigm shift, pushing back the evolutionary origin of methanogenesis to predate that of the Euryarchaeota. The methylotrophic methanogenesis found in the non-Euryarchaota distinguished itself from the predominantly hydrogenotrophic methanogens found in euryarchaeal orders as the former do not couple methanogenesis to carbon fixation through the reductive acetyl-CoA [Wood-Ljungdahl pathway (WLP)], which was interpreted as evidence for independent evolution of the two methanogenesis pathways. Here, we report the discovery of a complete and divergent hydrogenotrophic methanogenesis pathway in a thermophilic order of the Verstraetearchaeota, which we have named Candidatus Methanohydrogenales, as well as the presence of the WLP in the crenarchaeal order Desulfurococcales. Our findings support the ancient origin of hydrogenotrophic methanogenesis, suggest that methylotrophic methanogenesis might be a later adaptation of specific orders, and provide insight into how the transition from hydrogenotrophic to methylotrophic methanogenesis might have occurred.

Nucleic acid cleavage with a hyperthermophilic Cas9 from an uncultured Ignavibacterium.

Clustered regularly interspaced short palindromic repeats (CRISPR)-associated 9 (Cas9) systems have been effectively harnessed to engineer the genomes of organisms from across the tree of life. Nearly all currently characterized Cas9 proteins are derived from mesophilic bacteria, and canonical Cas9 systems are challenged by applications requiring enhanced stability or elevated temperatures. We discovered IgnaviCas9, a Cas9 protein from a hyperthermophilic Ignavibacterium identified through mini-metagenomic sequencing of samples from a hot spring. IgnaviCas9 is active at temperatures up to 100 °C in vitro, which enables DNA cleavage beyond the 44 °C limit of Streptococcus pyogenes Cas9 (SpyCas9) and the 70 °C limit of both Geobacillus stearothermophilus Cas9 (GeoCas9) and Geobacillus thermodenitrificans T12 Cas9 (ThermoCas9). As a potential application of this enzyme, we demonstrate that IgnaviCas9 can be used in bacterial RNA-seq library preparation to remove unwanted cDNA from 16s ribosomal rRNA without increasing the number of steps, thus underscoring the benefits provided by its exceptional thermostability in improving molecular biology and genomic workflows. IgnaviCas9 is an exciting addition to the CRISPR-Cas9 toolbox and expands its temperature range.

Single-cell RNA sequencing reveals intrinsic and extrinsic regulatory heterogeneity in yeast responding to stress.

From bacteria to humans, individual cells within isogenic populations can show significant variation in stress tolerance, but the nature of this heterogeneity is not clear. To investigate this, we used single-cell RNA sequencing to quantify transcript heterogeneity in single Saccharomyces cerevisiae cells treated with and without salt stress to explore population variation and identify cellular covariates that influence the stress-responsive transcriptome. Leveraging the extensive knowledge of yeast transcriptional regulation, we uncovered significant regulatory variation in individual yeast cells, both before and after stress. We also discovered that a subset of cells appears to decouple expression of ribosomal protein genes from the environmental stress response in a manner partly correlated with the cell cycle but unrelated to the yeast ultradian metabolic cycle. Live-cell imaging of cells expressing pairs of fluorescent regulators, including the transcription factor Msn2 with Dot6, Sfp1, or MAP kinase Hog1, revealed both coordinated and decoupled nucleocytoplasmic shuttling. Together with transcriptomic analysis, our results suggest that cells maintain a cellular filter against decoupled bursts of transcription factor activation but mount a stress response upon coordinated regulation, even in a subset of unstressed cells.

Long-term microfluidic tracking of coccoid cyanobacterial cells reveals robust control of division timing.

BACKGROUND: Cyanobacteria are important agents in global carbon and nitrogen cycling and hold great promise for biotechnological applications. Model organisms such as Synechocystis sp. and Synechococcus sp. have advanced our understanding of photosynthetic capacity and circadian behavior, mostly using population-level measurements in which the behavior of individuals cannot be monitored. Synechocystis sp. cells are small and divide slowly, requiring long-term experiments to track single cells. Thus, the cumulative effects of drift over long periods can cause difficulties in monitoring and quantifying cell growth and division dynamics. RESULTS: To overcome this challenge, we enhanced a microfluidic cell-culture device and developed an image analysis pipeline for robust lineage reconstruction. This allowed simultaneous tracking of many cells over multiple generations, and revealed that cells expand exponentially throughout their cell cycle. Generation times were highly correlated for sister cells, but not between mother and daughter cells. Relationships between birth size, division size, and generation time indicated that cell-size control was inconsistent with the "sizer" rule, where division timing is based on cell size, or the "timer" rule, where division occurs after a fixed time interval. Instead, single cell growth statistics were most consistent with the "adder" rule, in which division occurs after a constant increment in cell volume. Cells exposed to light-dark cycles exhibited growth and division only during the light period; dark phases pause but do not disrupt cell-cycle control. CONCLUSIONS: Our analyses revealed that the "adder" model can explain both the growth-related statistics of single Synechocystis cells and the correlation between sister cell generation times. We also observed rapid phenotypic response to light-dark transitions at the single cell level, highlighting the critical role of light in cyanobacterial cell-cycle control. Our findings suggest that by monitoring the growth kinetics of individual cells we can build testable models of circadian control of the cell cycle in cyanobacteria.

Reference-free Structural Variant Detection in Microbiomes via Long-read Coassembly Graphs.

Bacterial genome dynamics are vital for understanding the mechanisms underlying microbial adaptation, growth, and their broader impact on host phenotype. Structural variants (SVs), genomic alterations of 10 base pairs or more, play a pivotal role in driving evolutionary processes and maintaining genomic heterogeneity within bacterial populations. While SV detection in isolate genomes is relatively straightforward, metagenomes present broader challenges due to absence of clear reference genomes and presence of mixed strains. In response, our proposed method rhea, forgoes reference genomes and metagenome-assembled genomes (MAGs) by encompassing a single metagenome coassembly graph constructed from all samples in a series. The log fold change in graph coverage between subsequent samples is then calculated to call SVs that are thriving or declining throughout the series. We show rhea to outperform existing methods for SV and horizontal gene transfer (HGT) detection in two simulated mock metagenomes, which is particularly noticeable as the simulated reads diverge from reference genomes and an increase in strain diversity is incorporated. We additionally demonstrate use cases for rhea on series metagenomic data of environmental and fermented food microbiomes to detect specific sequence alterations between subsequent time and temperature samples, suggesting host advantage. Our innovative approach leverages raw read patterns rather than references or MAGs to include all sequencing reads in analysis, and thus provide versatility in studying SVs across diverse and poorly characterized microbial communities for more comprehensive insights into microbial genome dynamics.

Assembly of gut-derived bacterial communities follows "early-bird" resource utilization dynamics.

Diet can impact host health through changes to the gut microbiota, yet we lack mechanistic understanding linking nutrient availability and microbiota composition. Here, we use thousands of microbial communities cultured in vitro from human feces to uncover simple assembly rules and develop a predictive model of community composition upon addition of single nutrients from central carbon metabolism to a complex medium. Community membership was largely determined by the donor feces, whereas relative abundances were determined by the supplemental carbon source. The absolute abundance of most taxa was independent of the supplementing nutrient, due to the ability of fast-growing organisms to quickly exhaust their niche in the complex medium and then exploit and monopolize the supplemental carbon source. Relative abundances of dominant taxa could be predicted from the nutritional preferences and growth dynamics of species in isolation, and exceptions were consistent with strain-level variation in growth capabilities. Our study reveals that community assembly follows simple rules of nutrient utilization dynamics and provides a predictive framework for manipulating gut commensal communities through nutritional perturbations.

Prolonged delays in human microbiota transmission after a controlled antibiotic perturbation.

Humans constantly encounter new microbes, but few become long-term residents of the adult gut microbiome. Classical theories predict that colonization is determined by the availability of open niches, but it remains unclear whether other ecological barriers limit commensal colonization in natural settings. To disentangle these effects, we used a controlled perturbation with the antibiotic ciprofloxacin to investigate the dynamics of gut microbiome transmission in 22 households of healthy, cohabiting adults. Colonization was rare in three-quarters of antibiotic-taking subjects, whose resident strains rapidly recovered in the week after antibiotics ended. In contrast, the remaining antibiotic-taking subjects exhibited lasting responses, with extensive species losses and transient expansions of potential opportunistic pathogens. These subjects experienced elevated rates of commensal colonization, but only after long delays: many new colonizers underwent sudden, correlated expansions months after the antibiotic perturbation. Furthermore, strains that had previously transmitted between cohabiting partners rarely recolonized after antibiotic disruptions, showing that colonization displays substantial historical contingency. This work demonstrates that there remain substantial ecological barriers to colonization even after major microbiome disruptions, suggesting that dispersal interactions and priority effects limit the pace of community change.

Impact of a 7-day homogeneous diet on interpersonal variation in human gut microbiomes and metabolomes.

Gut microbiota metabolism of dietary compounds generates a vast array of microbiome-dependent metabolites (MDMs), which are highly variable between individuals. The uremic MDMs (uMDMs) phenylacetylglutamine (PAG), p-cresol sulfate (PCS), and indoxyl sulfate (IS) accumulate during renal failure and are associated with poor outcomes. Targeted dietary interventions may reduce toxic MDM generation; however, it is unclear if inter-individual differences in diet or gut microbiome dominantly contribute to MDM variance. Here, we use a 7-day homogeneous average American diet to standardize dietary precursor availability in 21 healthy individuals. During dietary homogeneity, the coefficient of variation in PAG, PCS, and IS (primary outcome) did not decrease, nor did inter-individual variation in most identified metabolites; other microbiome metrics showed no or modest responses to the intervention. Host identity and age are dominant contributors to variability in MDMs. These results highlight the potential need to pair dietary modification with microbial therapies to control MDM profiles.

The CIAMIB: a Large and Metabolically Diverse Collection of Inflammation-Associated Bacteria from the Murine Gut.

Gut inflammation directly impacts the growth and stability of commensal gut microbes and can lead to long-lasting changes in microbiota composition that can prolong or exacerbate disease states. While mouse models are used extensively to investigate the interplay between microbes and the inflamed state, the paucity of cultured mouse gut microbes has hindered efforts to determine causal relationships. To address this issue, we are assembling the Collection of Inflammation-Associated Mouse Intestinal Bacteria (CIAMIB). The initial release of this collection comprises 41 isolates of 39 unique bacterial species, covering 4 phyla and containing 10 previously uncultivated isolates, including 1 novel family and 7 novel genera. The collection significantly expands the number of available Muribaculaceae, Lachnospiraceae, and Coriobacteriaceae isolates and includes microbes from genera associated with inflammation, such as Prevotella and Klebsiella. We characterized the growth of CIAMIB isolates across a diverse range of nutritional conditions and predicted their metabolic potential and anaerobic fermentation capacity based on the genomes of these isolates. We also provide the first metabolic analysis of species within the genus Adlercreutzia, revealing these representatives to be nitrate-reducing and severely restricted in their ability to grow on carbohydrates. CIAMIB isolates are fully sequenced and available to the scientific community as a powerful tool to study host-microbiota interactions. IMPORTANCE Attempts to explore the role of the microbiota in animal physiology have resulted in large-scale efforts to cultivate the thousands of microbes that are associated with humans. In contrast, relatively few lab mouse-associated bacteria have been isolated, despite the fact that the overwhelming number of studies on the microbiota use laboratory mice that are colonized with microbes that are quite distinct from those in humans. Here, we report the results of a large-scale isolation of bacteria from the intestines of laboratory mice either prone to or suffering from gut inflammation. This collection comprises dozens of novel isolates, many of which represent the only cultured representatives of their genus or species. We report their basic growth characteristics and genomes and are making them widely available to the greater research community.

Establishment and characterization of stable, diverse, fecal-derived in vitro microbial communities that model the intestinal microbiota.

Efforts to probe the role of the gut microbiota in disease would benefit from a system in which patient-derived bacterial communities can be studied at scale. We addressed this by validating a strategy to propagate phylogenetically complex, diverse, stable, and highly reproducible stool-derived communities in vitro. We generated hundreds of in vitro communities cultured from diverse stool samples in various media; certain media generally preserved inoculum composition, and inocula from different subjects yielded source-specific community compositions. Upon colonization of germ-free mice, community composition was maintained, and the host proteome resembled the host from which the community was derived. Treatment with ciprofloxacin in vivo increased susceptibility to Salmonella invasion in vitro, and the in vitro response to ciprofloxacin was predictive of compositional changes observed in vivo, including the resilience and sensitivity of each Bacteroides species. These findings demonstrate that stool-derived in vitro communities can serve as a powerful system for microbiota research.

Quantifying rapid bacterial evolution and transmission within the mouse intestine.

Due to limitations on high-resolution strain tracking, selection dynamics during gut microbiota colonization and transmission between hosts remain mostly mysterious. Here, we introduced hundreds of barcoded Escherichia coli strains into germ-free mice and quantified strain-level dynamics and metagenomic changes. Mutations in genes involved in motility and metabolite utilization are reproducibly selected within days. Even with rapid selection, coprophagy enforced similar barcode distributions across co-housed mice. Whole-genome sequencing of hundreds of isolates revealed linked alleles that demonstrate between-host transmission. A population-genetics model predicts substantial fitness advantages for certain mutants and that migration accounted for ∼10% of the resident microbiota each day. Treatment with ciprofloxacin suggests interplay between selection and transmission. While initial colonization was mostly uniform, in two mice a bottleneck reduced diversity and selected for ciprofloxacin resistance in the absence of drug. These findings highlight the interplay between environmental transmission and rapid, deterministic selection during evolution of the intestinal microbiota.

Mini-Metagenomics and Nucleotide Composition Aid the Identification and Host Association of Novel Bacteriophage Sequences.

A broad spectrum of metagenomic and single cell sequencing techniques have become popular for dissecting environmental microbial diversity, leading to the characterization of thousands of novel microbial lineages. In addition to recovering bacterial and archaeal genomes, metagenomic assembly can also produce genomes of viruses that infect microbial cells. Because of their diversity, lack of marker genes, and small genome size, identifying novel bacteriophage sequences from metagenomic data is often challenging, especially when the objective is to establish phage-host relationships. The present work describes a computational approach that uses supervised learning to classify metagenomic contigs as phage or non-phage as well as assigning phage taxonomy based on tetranucleotide frequencies. Furthermore, the method assigns phage-host relationships using co-occurrence statistics derived from a recently developed mini-metagenomic experimental technique. This work evaluates method performance at identifying viral contigs and predicting taxonomic classification using publicly available references. Then, using two mini-metagenomic datasets, over 100 novel phage contigs from hot spring samples of Yellowstone National Park are identified and assigned to putative microbial hosts. Results of this work demonstrate the value of combining viral sequence identification with mini-metagenomic experimental methods to understand the microbial ecosystem.

Recovery of the Gut Microbiota after Antibiotics Depends on Host Diet, Community Context, and Environmental Reservoirs.

Antibiotics alter microbiota composition and increase infection susceptibility. However, the generalizable effects of antibiotics on and the contribution of environmental variables to gut commensals remain unclear. To address this, we tracked microbiota dynamics with high temporal and taxonomic resolution during antibiotic treatment in a controlled murine system by isolating variables such as diet, treatment history, and housing co-inhabitants. Human microbiotas were remarkably resilient and recovered during antibiotic treatment, with transient dominance of resistant Bacteroides and taxa-asymmetric diversity reduction. In certain cases, in vitro sensitivities were not predictive of in vivo responses, underscoring the significance of host and community context. A fiber-deficient diet exacerbated microbiota collapse and delayed recovery. Species replacement through cross housing after ciprofloxacin treatment established resilience to a second treatment. Single housing drastically disrupted recovery, highlighting the importance of environmental reservoirs. Our findings highlight deterministic microbiota adaptations to perturbations and the translational potential for modulating diet, sanitation, and microbiota composition during antibiotics.

Microfluidic-based mini-metagenomics enables discovery of novel microbial lineages from complex environmental samples.

Metagenomics and single-cell genomics have enabled genome discovery from unknown branches of life. However, extracting novel genomes from complex mixtures of metagenomic data can still be challenging and represents an ill-posed problem which is generally approached with ad hoc methods. Here we present a microfluidic-based mini-metagenomic method which offers a statistically rigorous approach to extract novel microbial genomes while preserving single-cell resolution. We used this approach to analyze two hot spring samples from Yellowstone National Park and extracted 29 new genomes, including three deeply branching lineages. The single-cell resolution enabled accurate quantification of genome function and abundance, down to 1% in relative abundance. Our analyses of genome level SNP distributions also revealed low to moderate environmental selection. The scale, resolution, and statistical power of microfluidic-based mini-metagenomics make it a powerful tool to dissect the genomic structure of microbial communities while effectively preserving the fundamental unit of biology, the single cell.