Measurement of bacterial replication rates in microbial communities

Genome-centric investigation of anaerobic digestion using sustainable second and third generation substrates

Biogas production through co-digestion of second and third generation substrates is an environmentally sustainable approach. Green willow biomass, chicken manure waste and microalgae biomass substrates were combined in the anaerobic digestion experiments. Biochemical methane potential test showed that biogas yields of co-digestions were significantly higher compared to the yield when energy willow was the sole substrate. To scale up the experiment continuous stirred-tank reactors (CSRTs) are employed, digestion parameters are monitored. Furthermore, genome-centric metagenomics approach was employed to gain functional insight into the complex anaerobic decomposing process. This revealed the importance of Firmicutes, Actinobacteria, Proteobacteria and Bacteroidetes phyla as major bacterial participants, while Methanomicrobia and Methanobacteria represented the archaeal constituents of the communities. The bacterial phyla were shown to perform the carbohydrate hydrolysis. Among the representatives of long-chain carbohydrate hydrolysing microbes Bin_61: Clostridia is newly identified metagenome assembled genome (MAG) and Bin_13: DTU010 sp900018335 is common and abundant in all CSTRs. Methanogenesis was linked to the slow-growing members of the community, where hydrogenotrophic methanogen species Methanoculleus (Bin_10) and Methanobacterium (Bin_4) predominate. A sensitive balance between H₂ producers and consumers was shown to be critical for stable biomethane production and efficient waste biodegradation.

Virus-associated organosulfur metabolism in human and environmental systems

Viruses influence the fate of nutrients and human health by killing microorganisms and altering metabolic processes. Organosulfur metabolism and biologically derived hydrogen sulfide play dynamic roles in manifestation of diseases, infrastructure degradation, and essential biological processes. Although microbial organosulfur metabolism is well studied, the role of viruses in organosulfur metabolism is unknown. Here, we report the discovery of 39 gene families involved in organosulfur metabolism encoded by 3,749 viruses from diverse ecosystems, including human microbiomes. The viruses infect organisms from all three domains of life. Six gene families encode for enzymes that degrade organosulfur compounds into sulfide, whereas others manipulate organosulfur compounds and may influence sulfide production. We show that viral metabolic genes encode key enzymatic domains, are translated into protein, and are maintained after recombination, and sulfide provides a fitness advantage to viruses. Our results reveal viruses as drivers of organosulfur metabolism with important implications for human and environmental health.

Functional diversity of microbial communities in inactive seafloor sulfide deposits

The seafloor sulfide structures of inactive vents are known to host abundant and diverse microorganisms potentially supported by mineralogy of sulfides. However, little is known about the diversity and distribution of microbial functions. Here, we used genome-resolved metagenomics to predict microbial metabolic functions and the contribution of horizontal gene transfer to the functionality of microorganisms inhabiting several hydrothermally inactive seafloor deposits among globally distributed deep-sea vent fields. Despite of geographically distant vent fields, similar microbial community patterns were observed with the dominance of Gammaproteobacteria, Bacteroidota and previously overlooked Candidatus Patescibacteria. Metabolically flexible Gammaproteobacteria are major potential primary producers utilizing mainly sulfur, iron and hydrogen as electron donors coupled with oxygen and nitrate respiration for chemolithoautotrophic growth. In addition to heterotrophic microorganisms like free-living Bacteroidota, Ca. Patescibacteria potentially perform fermentative recycling of organic carbon. Finally, we provided evidence that many functional genes that are central to energy metabolism have been laterally transferred among members within the community and largely within the same class. Taken together, these findings shed light on microbial ecology and evolution in inactive seafloor sulfide deposits after the cessation of hydrothermal activities.

The sRNA Regulated Protein DdbA Is Involved in Development and Maintenance of the Chlamydia trachomatis EB Cell Form

The chlamydial small non coding RNA, IhtA, regulates the expression of both HctA and DdbA, the uncharacterized product of the C. trachomatis L2 CTL0322 gene. HctA is a small, highly basic, DNA binding protein that is expressed late in development and mediates the condensation of the genome during RB to EB differentiation. DdbA is conserved throughout the chlamydial lineage, and is predicted to express a small, basic, cytoplasmic protein. As it is common for sRNAs to regulate multiple mRNAs within the same physiological pathway, we hypothesize that DdbA, like HctA, is involved in RB to EB differentiation. Here, we show that DdbA is a DNA binding protein, however unlike HctA, DdbA does not contribute to genome condensation but instead likely has nuclease activity. Using a DdbA temperature sensitive mutant, we show that DdbAts creates inclusions indistinguishable from WT L2 in size and that early RB replication is likewise similar at the nonpermissive temperature. However, the number of DdbAts infectious progeny is dramatically lower than WT L2 overall, although production of EBs is initiated at a similar time. The expression of a late gene reporter construct followed live at 40°C indicates that late gene expression is severely compromised in the DdbAts strain. Viability assays, both in host cells and in axenic media indicate that the DdbAts strain is defective in the maintenance of EB infectivity. Additionally, using Whole Genome Sequencing we demonstrate that chromosome condensation is temporally separated from DNA replication during the RB to EB transition. Although DdbA does not appear to be directly involved in this process, our data suggest it is a DNA binding protein that is important in the production and maintenance of infectivity of the EB, perhaps by contributing to the remodeling of the EB chromosome.

Recurrent Urinary Tract Infections: Unraveling the Complicated Environment of Uncomplicated rUTIs

Urinary tract infections (UTIs) are frequent in humans, affecting the upper and lower urinary tract. Present diagnosis relies on the positive culture of uropathogenic bacteria from urine and clinical markers of inflammation of the urinary tract. The bladder is constantly challenged by adverse environmental stimuli which influence urinary tract physiology, contributing to a dysbiotic environment. Simultaneously, pathogens are primed by environmental stressors such as antibiotics, favoring recurrent UTIs (rUTIs), resulting in chronic illness. Due to different confounders for UTI onset, a greater understanding of the fundamental environmental mechanisms and microbial ecology of the human urinary tract is required. Such advancements could promote the tandem translation of bench and computational studies for precision treatments and clinical management of UTIs. Therefore, there is an urgent need to understand the ecological interactions of the human urogenital microbial communities which precede rUTIs. This review aims to outline the mechanistic aspects of rUTI ecology underlying dysbiosis between both the human microbiome and host physiology which predisposes humans to rUTIs. By assessing the applications of next generation and systems level methods, we also recommend novel approaches to elucidate the systemic consequences of rUTIs which requires an integrated approach for successful treatment. To this end, we will provide an outlook towards the so-called 'uncomplicated environment of UTIs', a holistic and systems view that applies ecological principles to define patient-specific UTIs. This perspective illustrates the need to withdraw from traditional reductionist perspectives in infection biology and instead, a move towards a systems-view revolving around patient-specific pathophysiology during UTIs.

Umibato: estimation of time-varying microbial interaction using continuous-time regression hidden Markov model

MOTIVATION

Accumulating evidence has highlighted the importance of microbial interaction networks. Methods have been developed for estimating microbial interaction networks, of which the generalized Lotka-Volterra equation (gLVE)-based method can estimate a directed interaction network. The previous gLVE-based method for estimating microbial interaction networks did not consider time-varying interactions.

RESULTS

In this study, we developed unsupervised learning-based microbial interaction inference method using Bayesian estimation (Umibato), a method for estimating time-varying microbial interactions. The Umibato algorithm comprises Gaussian process regression (GPR) and a new Bayesian probabilistic model, the continuous-time regression hidden Markov model (CTRHMM). Growth rates are estimated by GPR, and interaction networks are estimated by CTRHMM. CTRHMM can estimate time-varying interaction networks using interaction states, which are defined as hidden variables. Umibato outperformed the existing methods on synthetic datasets. In addition, it yielded reasonable estimations in experiments on a mouse gut microbiota dataset, thus providing novel insights into the relationship between consumed diets and the gut microbiota.

AVAILABILITY AND IMPLEMENTATION

The C++ and python source codes of the Umibato software are available at https://github.com/shion-h/Umibato.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Integrating Systems and Synthetic Biology to Understand and Engineer Microbiomes

Microbiomes are complex and ubiquitous networks of microorganisms whose seemingly limitless chemical transformations could be harnessed to benefit agriculture, medicine, and biotechnology. The spatial and temporal changes in microbiome composition and function are influenced by a multitude of molecular and ecological factors. This complexity yields both versatility and challenges in designing synthetic microbiomes and perturbing natural microbiomes in controlled, predictable ways. In this review, we describe factors that give rise to emergent spatial and temporal microbiome properties and the meta-omics and computational modeling tools that can be used to understand microbiomes at the cellular and system levels. We also describe strategies for designing and engineering microbiomes to enhance or build novel functions. Throughout the review, we discuss key knowledge and technology gaps for elucidating the networks and deciphering key control points for microbiome engineering, and highlight examples where multiple omics and modeling approaches can be integrated to address these gaps.

Replication Dynamics for Six Gram-Negative Bacterial Species during Bloodstream Infection

Bloodstream infections (BSI) are a major public health burden due to high mortality rates and the cost of treatment. The impact of BSI is further compounded by a rise in antibiotic resistance among Gram-negative species associated with these infections. Escherichia coli, Serratia marcescens, Klebsiella pneumoniae, Enterobacter hormaechei, Citrobacter freundii, and Acinetobacter baumannii are all common causes of BSI, which can be recapitulated in a murine model. The objective of this study was to characterize infection kinetics and bacterial replication rates during bacteremia for these six pathogens to gain a better understanding of bacterial physiology during infection. Temporal observations of bacterial burdens of the tested species demonstrated varied abilities to establish colonization in the spleen, liver, or kidney. K. pneumoniae and S. marcescens expanded rapidly in the liver and kidney, respectively. Other organisms, such as C. freundii and E. hormaechei, were steadily cleared from all three target organs throughout the infection. In situ replication rates measured by whole-genome sequencing of bacterial DNA recovered from murine spleens demonstrated that each species was capable of sustained replication at 24 h postinfection, and several species demonstrated <60-min generation times. The relatively short generation times observed in the spleen were in contrast to an overall decrease in bacterial burden for some species, suggesting that the rate of immune-mediated clearance exceeded replication. Furthermore, bacterial generation times measured in the murine spleen approximated those measured during growth in human serum cultures. Together, these findings provide insight into the infection kinetics of six medically important species during bacteremia.

Preterm Birth Is Correlated With Increased Oral Originated Microbiome in the Gut

Background

Preterm birth is one of the leading causes of perinatal morbidity and mortality. Gut microbiome dysbiosis is closely related to adverse pregnancy outcomes. However, the role of the gut microbiome in the pathogenesis of preterm birth remains poorly studied.

Method

We collected fecal samples from 41 women (cases presenting with threatened preterm labor =19, 11 of which delivered preterm; gestational age-matched no-labor controls, all of which delivered at term = 22) were recruited for the study. We performed 16S rRNA amplicon sequencing to compare the composition of the gut microbiome in threatened preterm labor cases and controls and among women who delivered preterm and at term. By annotating taxonomic biomarkers with the Human Oral Microbiome Database, we observed an increased abundance of potential oral-to-gut bacteria in preterm patients.

Results

Patients with preterm birth showed a distinct gut microbiome dysbiosis compared with those who delivered at term. Opportunistic pathogens, particularly Porphyromonas, Streptococcus, Fusobacterium, and Veillonella, were enriched, whereas Coprococcus and Gemmiger were markedly depleted in the preterm group. Most of the enriched bacteria were annotated oral bacteria using the Human Oral Microbiome Database. These potential oral-to-gut bacteria were correlated with clinical parameters that reflected maternal and fetal status.

Conclusions

This study suggests that patients who deliver preterm demonstrate altered gut microbiome that may contain higher common oral bacteria.

Helarchaeota and co-occurring sulfate-reducing bacteria in subseafloor sediments from the Costa Rica Margin

Deep sediments host many archaeal lineages, including the Asgard superphylum which contains lineages predicted to require syntrophic partnerships. Our knowledge about sedimentary archaeal diversity and their metabolic pathways and syntrophic partners is still very limited. We present here new genomes of Helarchaeota and the co-occurring sulfate-reducing bacteria (SRB) recovered from organic-rich sediments off Costa Rica Margin. Phylogenetic analyses revealed three new metagenome-assembled genomes (MAGs) affiliating with Helarchaeota, each of which has three variants of the methyl-CoM reductase-like (MCR-like) complex that may enable them to oxidize short-chain alkanes anaerobically. These Helarchaeota have no multi-heme cytochromes but have Group 3b and Group 3c [NiFe] hydrogenases, and formate dehydrogenase, and therefore have the capacity to transfer the reducing equivalents (in the forms of hydrogen and formate) generated from alkane oxidation to external partners. We also recovered five MAGs of SRB affiliated with the class of Desulfobacteria, two of which showed relative abundances (represented by genome coverages) positively correlated with those of the three Helarchaeota. Genome analysis suggested that these SRB bacteria have the capacity of H2 and formate utilization and could facilitate electron transfers from other organisms by means of these reduced substances. Their co-occurrence and metabolic features suggest that Helarchaeota may metabolize synergistically with some SRB, and together exert an important influence on the carbon cycle by mitigating the hydrocarbon emission from sediments to the overlying ocean.

Genomic adaptations enabling Acidithiobacillus distribution across wide-ranging hot spring temperatures and pHs

BACKGROUND

Terrestrial hot spring settings span a broad spectrum of physicochemistries. Physicochemical parameters, such as pH and temperature, are key factors influencing differences in microbial composition across diverse geothermal areas. Nonetheless, analysis of hot spring pools from the Taupo Volcanic Zone (TVZ), New Zealand, revealed that some members of the bacterial genus, Acidithiobacillus, are prevalent across wide ranges of hot spring pHs and temperatures. To determine the genomic attributes of Acidithiobacillus that inhabit such diverse conditions, we assembled the genomes of 19 uncultivated hot spring Acidithiobacillus strains from six geothermal areas and compared these to 37 publicly available Acidithiobacillus genomes from various habitats.

RESULTS

Analysis of 16S rRNA gene amplicons from 138 samples revealed that Acidithiobacillus comprised on average 11.4 ± 16.8% of hot spring prokaryotic communities, with three Acidithiobacillus amplicon sequence variants (ASVs) (TVZ_G1, TVZ_G2, TVZ_G3) accounting for > 90% of Acidithiobacillus in terms of relative abundance, and occurring in 126 out of 138 samples across wide ranges of temperature (17.5-92.9 °C) and pH (1.0-7.5). We recovered 19 environmental genomes belonging to each of these three ASVs, as well as a fourth related group (TVZ_G4). Based on genome average nucleotide identities, the four groups (TVZ_G1-TVZ_G4) constitute distinct species (ANI < 96.5%) of which three are novel Acidithiobacillus species (TVZ_G2-TVZ_G4) and one belongs to Acidithiobacillus caldus (TVZ_G1). All four TVZ Acidithiobacillus groups were found in hot springs with temperatures above the previously known limit for the genus (up to 40 °C higher), likely due to significantly higher proline and GC contents than other Acidithiobacillus species, which are known to increase thermostability. Results also indicate hot spring-associated Acidithiobacillus have undergone genome streamlining, likely due to thermal adaptation. Moreover, our data suggest that Acidithiobacillus prevalence across varied hot spring pHs is supported by distinct strategies, whereby TVZ_G2-TVZ_G4 regulate pH homeostasis mostly through Na⁺/H⁺ antiporters and proton-efflux ATPases, whereas TVZ_G1 mainly relies on amino acid decarboxylases.

CONCLUSIONS

This study provides insights into the distribution of Acidithiobacillus species across diverse hot spring physichochemistries and determines genomic features and adaptations that potentially enable Acidithiobacillus species to colonize a broad range of temperatures and pHs in geothermal environments. Video Abstract.

Genomic characterization of three novel Desulfobacterota classes expand the metabolic and phylogenetic diversity of the phylum

We report on the genomic characterization of three novel classes in the phylum Desulfobacterota. One class (proposed name Candidatus 'Anaeroferrophillalia') was characterized by heterotrophic growth capacity, either fermentatively or utilizing polysulfide, tetrathionate or thiosulfate as electron acceptors. In the absence of organic carbon sources, autotrophic growth via the Wood-Ljungdahl (WL) pathway and using hydrogen or Fe(II) as an electron donor is also inferred for members of the 'Anaeroferrophillalia'. The second class (proposed name Candidatus 'Anaeropigmentia') was characterized by its capacity for growth at low oxygen concentration, and the capacity to synthesize the methyl/alkyl carrier CoM, an ability that is prevalent in the archaeal but rare in the bacterial domain. Pigmentation is inferred from the capacity for carotenoid (lycopene) production. The third class (proposed name Candidatus 'Zymogenia') was characterized by fermentative heterotrophic growth capacity, broad substrate range and the adaptation of some of its members to hypersaline habitats. Analysis of the distribution pattern of all three classes showed their occurrence as rare community members in multiple habitats, with preferences for anaerobic terrestrial, freshwater and marine environments over oxygenated (e.g. pelagic ocean and agricultural land) settings. Special preference for some members of the class Candidatus 'Zymogenia' for hypersaline environments such as hypersaline microbial mats and lagoons was observed.

Microbial Hotspots in Lithic Microhabitats Inferred from DNA Fractionation and Metagenomics in the Atacama Desert

The existence of microbial activity hotspots in temperate regions of Earth is driven by soil heterogeneities, especially the temporal and spatial availability of nutrients. Here we investigate whether microbial activity hotspots also exist in lithic microhabitats in one of the most arid regions of the world, the Atacama Desert in Chile. While previous studies evaluated the total DNA fraction to elucidate the microbial communities, we here for the first time use a DNA separation approach on lithic microhabitats, together with metagenomics and other analysis methods (i.e., ATP, PLFA, and metabolite analysis) to specifically gain insights on the living and potentially active microbial community. Our results show that hypolith colonized rocks are microbial hotspots in the desert environment. In contrast, our data do not support such a conclusion for gypsum crust and salt rock environments, because only limited microbial activity could be observed. The hypolith community is dominated by phototrophs, mostly Cyanobacteria and Chloroflexi, at both study sites. The gypsum crusts are dominated by methylotrophs and heterotrophic phototrophs, mostly Chloroflexi, and the salt rocks (halite nodules) by phototrophic and halotolerant endoliths, mostly Cyanobacteria and Archaea. The major environmental constraints in the organic-poor arid and hyperarid Atacama Desert are water availability and UV irradiation, allowing phototrophs and other extremophiles to play a key role in desert ecology.

Adaptive differentiation and rapid evolution of a soil bacterium along a climate gradient

Increasing evidence suggests that evolutionary processes frequently shape ecological patterns; however, most microbiome studies thus far have focused on only the ecological responses of these communities. By using parallel field experiments and focusing in on a model soil bacterium, we showed that bacterial “species” are differentially adapted to local climates, leading to changes in their composition. Furthermore, we detected strain-level evolution, providing direct evidence that both ecological and evolutionary processes operate on annual timescales. The consideration of eco-evolutionary dynamics may therefore be important to understand the response of soil microbiomes to future environmental change.

Microbial community responses to environmental change are largely associated with ecological processes; however, the potential for microbes to rapidly evolve and adapt remains relatively unexplored in natural environments. To assess how ecological and evolutionary processes simultaneously alter the genetic diversity of a microbiome, we conducted two concurrent experiments in the leaf litter layer of soil over 18 mo across a climate gradient in Southern California. In the first experiment, we reciprocally transplanted microbial communities from five sites to test whether ecological shifts in ecotypes of the abundant bacterium, Curtobacterium, corresponded to past adaptive differentiation. In the transplanted communities, ecotypes converged toward that of the native communities growing on a common litter substrate. Moreover, these shifts were correlated with community-weighted mean trait values of the Curtobacterium ecotypes, indicating that some of the trait variation among ecotypes could be explained by local adaptation to climate conditions. In the second experiment, we transplanted an isogenic Curtobacterium strain and tracked genomic mutations associated with the sites across the same climate gradient. Using a combination of genomic and metagenomic approaches, we identified a variety of nonrandom, parallel mutations associated with transplantation, including mutations in genes related to nutrient acquisition, stress response, and exopolysaccharide production. Together, the field experiments demonstrate how both demographic shifts of previously adapted ecotypes and contemporary evolution can alter the diversity of a soil microbiome on the same timescale.

Spatiotemporal Variations in Growth Rate and Virulence Plasmid Copy Number during Yersinia pseudotuberculosis Infection

Pathogenic Yersinia spp. depend on the activity of a potent virulence plasmid-encoded ysc/yop type 3 secretion system (T3SS) to colonize hosts and cause disease. It was recently shown that Yersinia pseudotuberculosis upregulates the virulence plasmid copy number (PCN) during infection and that the resulting elevated gene dose of plasmid-encoded T3SS genes is essential for virulence.

Pathogenic Yersinia spp. depend on the activity of a potent virulence plasmid-encoded ysc/yop type 3 secretion system (T3SS) to colonize hosts and cause disease. It was recently shown that Yersinia pseudotuberculosis upregulates the virulence plasmid copy number (PCN) during infection and that the resulting elevated gene dose of plasmid-encoded T3SS genes is essential for virulence. When and how this novel regulatory mechanism is deployed and regulates the replication of the virulence plasmid during infection is unknown. In the present study, we applied droplet digital PCR (ddPCR) to investigate the dynamics of Y. pseudotuberculosis virulence PCN variations and growth rates in infected mouse organs. We demonstrated that both PCN and growth varied in different tissues and over time throughout the course of infection, indicating that the bacteria adapted to discrete microenvironments during infection. The PCN was highest in Peyer’s patches and cecum during the clonal invasive phase of the infection, while the highest growth rates were found in the draining mesenteric lymph nodes. In deeper, systemic organs, the PCN was lower and more modest growth rates were recorded. Our study indicates that increased gene dosage of the plasmid-encoded T3SS genes is most important early in the infection during invasion of the host. The described ddPCR approach will greatly simplify analyses of PCN, growth dynamics, and bacterial loads in infected tissues and will be readily applicable to other infection models.

Long-term evolution and short-term adaptation of microbiota strains and sub-strains in mice

Isobiotic mice, with an identical stable microbiota composition, potentially allow models of host-microbial mutualism to be studied over time and between different laboratories. To understand microbiota evolution in these models, we carried out a 6-year experiment in mice colonized with 12 representative taxa. Increased non-synonymous to synonymous mutation rates indicate positive selection in multiple taxa, particularly for genes annotated for nutrient acquisition or replication. Microbial sub-strains that evolved within a single taxon can stably coexist, consistent with niche partitioning of ecotypes in the complex intestinal environment. Dietary shifts trigger rapid transcriptional adaptation to macronutrient and micronutrient changes in individual taxa and alterations in taxa biomass. The proportions of different sub-strains are also rapidly altered after dietary shift. This indicates that microbial taxa within a mouse colony adapt to changes in the intestinal environment by long-term genomic positive selection and short-term effects of transcriptional reprogramming and adjustments in sub-strain proportions.

Taxonomic and functional analyses of intact microbial communities thriving in extreme, astrobiology-relevant, anoxic sites

BACKGROUND

Extreme terrestrial, analogue environments are widely used models to study the limits of life and to infer habitability of extraterrestrial settings. In contrast to Earth's ecosystems, potential extraterrestrial biotopes are usually characterized by a lack of oxygen.

METHODS

In the MASE project (Mars Analogues for Space Exploration), we selected representative anoxic analogue environments (permafrost, salt-mine, acidic lake and river, sulfur springs) for the comprehensive analysis of their microbial communities. We assessed the microbiome profile of intact cells by propidium monoazide-based amplicon and shotgun metagenome sequencing, supplemented with an extensive cultivation effort.

RESULTS

The information retrieved from microbiome analyses on the intact microbial community thriving in the MASE sites, together with the isolation of 31 model microorganisms and successful binning of 15 high-quality genomes allowed us to observe principle pathways, which pinpoint specific microbial functions in the MASE sites compared to moderate environments. The microorganisms were characterized by an impressive machinery to withstand physical and chemical pressures. All levels of our analyses revealed the strong and omnipresent dependency of the microbial communities on complex organic matter. Moreover, we identified an extremotolerant cosmopolitan group of 34 poly-extremophiles thriving in all sites.

CONCLUSIONS

Our results reveal the presence of a core microbiome and microbial taxonomic similarities between saline and acidic anoxic environments. Our work further emphasizes the importance of the environmental, terrestrial parameters for the functionality of a microbial community, but also reveals a high proportion of living microorganisms in extreme environments with a high adaptation potential within habitability borders. Video abstract.

Genome-resolved metagenomics reveals site-specific diversity of episymbiotic CPR bacteria and DPANN archaea in groundwater ecosystems

Candidate phyla radiation (CPR) bacteria and DPANN archaea are unisolated, small-celled symbionts that are often detected in groundwater. The effects of groundwater geochemistry on the abundance, distribution, taxonomic diversity and host association of CPR bacteria and DPANN archaea has not been studied. Here, we performed genome-resolved metagenomic analysis of one agricultural and seven pristine groundwater microbial communities and recovered 746 CPR and DPANN genomes in total. The pristine sites, which serve as local sources of drinking water, contained up to 31% CPR bacteria and 4% DPANN archaea. We observed little species-level overlap of metagenome-assembled genomes (MAGs) across the groundwater sites, indicating that CPR and DPANN communities may be differentiated according to physicochemical conditions and host populations. Cryogenic transmission electron microscopy imaging and genomic analyses enabled us to identify CPR and DPANN lineages that reproducibly attach to host cells and showed that the growth of CPR bacteria seems to be stimulated by attachment to host-cell surfaces. Our analysis reveals site-specific diversity of CPR bacteria and DPANN archaea that coexist with diverse hosts in groundwater aquifers. Given that CPR and DPANN organisms have been identified in human microbiomes and their presence is correlated with diseases such as periodontitis, our findings are relevant to considerations of drinking water quality and human health.

Metagenomics and electron microscopy are combined to analyse the diversity of episymbiotic CPR bacteria and DPANN archaea in eight groundwater communities.

Distribution of Mixotrophy and Desiccation Survival Mechanisms across Microbial Genomes in an Arid Biological Soil Crust Community

Desert surface soils devoid of plant cover are populated by a variety of microorganisms, many with yet unresolved physiologies and lifestyles. Nevertheless, a common feature vital for these microorganisms inhabiting arid soils is their ability to survive long drought periods and reactivate rapidly in rare incidents of rain. Chemolithotrophic processes such as oxidation of atmospheric hydrogen and carbon monoxide are suggested to be a widespread energy source to support dormancy and resuscitation in desert soil microorganisms. Here, we assessed the distribution of chemolithotrophic, phototrophic, and desiccation-related metabolic potential among microbial populations in arid biological soil crusts (BSCs) from the Negev Desert, Israel, via population-resolved metagenomic analysis. While the potential to utilize light and atmospheric hydrogen as additional energy sources was widespread, carbon monoxide oxidation was less common than expected. The ability to utilize continuously available energy sources might decrease the dependency of mixotrophic populations on organic storage compounds and carbon provided by the BSC-founding cyanobacteria. Several populations from five different phyla besides the cyanobacteria encoded CO2 fixation potential, indicating further potential independence from photoautotrophs. However, we also found population genomes with a strictly heterotrophic genetic repertoire. The highly abundant Rubrobacteraceae (Actinobacteriota) genomes showed particular specialization for this extreme habitat, different from their closest cultured relatives. Besides the ability to use light and hydrogen as energy sources, they encoded extensive O2 stress protection and unique DNA repair potential. The uncovered differences in metabolic potential between individual, co-occurring microbial populations enable predictions of their ecological niches and generation of hypotheses on the dynamics and interactions among them.IMPORTANCE This study represents a comprehensive community-wide genome-centered metagenome analysis of biological soil crust (BSC) communities in arid environments, providing insights into the distribution of genes encoding different energy generation mechanisms, as well as survival strategies, among populations in an arid soil ecosystem. It reveals the metabolic potential of several uncultured and previously unsequenced microbial genera, families, and orders, as well as differences in the metabolic potential between the most abundant BSC populations and their cultured relatives, highlighting once more the danger of inferring function on the basis of taxonomy. Assigning functional potential to individual populations allows for the generation of hypotheses on trophic interactions and activity patterns in arid soil microbial communities and represents the basis for future resuscitation and activity studies of the system, e.g., involving metatranscriptomics.

The distribution of cellular turnover in the human body

We integrated ubiquity, mass and lifespan of all major cell types to achieve a comprehensive quantitative description of cellular turnover. We found a total cellular mass turnover of 80 ± 20 grams per day, dominated by blood cells and gut epithelial cells. In terms of cell numbers, close to 90% of the (0.33 ± 0.02) × 10¹² cells per day turnover was blood cells.

Size Matters: Ultra-small and Filterable Microorganisms in the Environment

Ultra-small microorganisms are ubiquitous in Earth’s environments. Ultramicrobacteria, which are defined as having a cell volume of <0.1 μm³, are often numerically dominant in aqueous environments. Cultivated representatives among these bacteria, such as members of the marine SAR11 clade (e.g., “Candidatus Pelagibacter ubique”) and freshwater Actinobacteria and Betaproteobacteria, possess highly streamlined, small genomes and unique ecophysiological traits. Many ultramicrobacteria may pass through a 0.2-μm-pore-sized filter, which is commonly used for filter sterilization in various fields and processes. Cultivation efforts focusing on filterable small microorganisms revealed that filtered fractions contained not only ultramicrocells (i.e., miniaturized cells because of external factors) and ultramicrobacteria, but also slender filamentous bacteria sometimes with pleomorphic cells, including a special reference to members of Oligoflexia, the eighth class of the phylum Proteobacteria. Furthermore, the advent of culture-independent “omics” approaches to filterable microorganisms yielded the existence of candidate phyla radiation (CPR) bacteria (also referred to as “Ca. Patescibacteria”) and ultra-small members of DPANN (an acronym of the names of the first phyla included in this superphyla) archaea. Notably, certain groups in CPR and DPANN are predicted to have minimal or few biosynthetic capacities, as reflected by their extremely small genome sizes, or possess no known function. Therefore, filtered fractions contain a greater variety and complexity of microorganisms than previously expected. This review summarizes the broad diversity of overlooked filterable agents remaining in “sterile” (<0.2-μm filtered) environmental samples.

Geochemical transition zone powering microbial growth in subsurface sediments

No other environment hosts as many microbial cells as the marine sedimentary biosphere. While the majority of these cells are expected to be alive, they are speculated to be persisting in a state of maintenance without net growth due to extreme starvation. Here, we report evidence for in situ growth of anaerobic ammonium-oxidizing (anammox) bacteria in ∼80,000-y-old subsurface sediments from the Arctic Mid-Ocean Ridge. The growth is confined to the nitrate-ammonium transition zone (NATZ), a widespread geochemical transition zone where most of the upward ammonium flux from deep anoxic sediments is being consumed. In this zone the anammox bacteria abundances, assessed by quantification of marker genes, consistently displayed a four order of magnitude increase relative to adjacent layers in four cores. This subsurface cell increase coincides with a markedly higher power supply driven mainly by intensified anammox reaction rates, thereby providing a quantitative link between microbial proliferation and energy availability. The reconstructed draft genome of the dominant anammox bacterium showed an index of replication (iRep) of 1.32, suggesting that 32% of this population was actively replicating. The genome belongs to a Scalindua species which we name Candidatus Scalindua sediminis, so far exclusively found in marine sediments. It has the capacity to utilize urea and cyanate and a mixotrophic lifestyle. Our results demonstrate that specific microbial groups are not only able to survive unfavorable conditions over geological timescales, but can proliferate in situ when encountering ideal conditions with significant consequences for biogeochemical nitrogen cycling.

Proximal colon-derived O-glycosylated mucus encapsulates and modulates the microbiota

Colon mucus segregates the intestinal microbiota from host tissues, but how it organizes to function throughout the colon is unclear. In mice, we found that colon mucus consists of two distinct O-glycosylated entities of Muc2: a major form produced by the proximal colon, which encapsulates the fecal material including the microbiota, and a minor form derived from the distal colon, which adheres to the major form. The microbiota directs its own encapsulation by inducing Muc2 production from proximal colon goblet cells. In turn, O-glycans on proximal colon-derived Muc2 modulate the structure and function of the microbiota as well as transcription in the colon mucosa. Our work shows how proximal colon control of mucin production is an important element in the regulation of host-microbiota symbiosis.

Genome-Resolved Metagenomics and Antibiotic Resistance Genes Analysis in Reclaimed Water Distribution Systems

Water reuse is increasingly pursued to alleviate global water scarcity. However, the wastewater treatment process does not achieve full removal of biological contaminants from wastewater, hence microorganisms and their genetic elements can be disseminated into the reclaimed water distribution systems (RWDS). In this study, reclaimed water samples are investigated via metagenomics to assess their bacterial diversity, metagenome-assembled genomes (MAGs) and antibiotic resistance genes (ARGs) at both point of entry (POE) and point of use (POU) in 3 RWDS. The number of shared bacterial orders identified by metagenome was higher at the POE than POU among the three sites, indicating that specific conditions in RWDS can cause further differentiation in the microbial communities at the end of the distribution system. Two bacterial orders, namely Rhizobiales and Sphingomonadales, had high replication rates in two of the examined RWDS (i.e., site A and B), and were present in higher relative abundance in POU than at POE. In addition, MAG and ARG relative abundance exhibited a strong correlation (R2 = 0.58) in POU, indicating that bacteria present in POU may have a high incidence of ARG. Specifically, resistance genes associated with efflux pump mechanisms (e.g., adeF and qacH) increased in its relative abundance from POU to POE at two of the RWDS (i.e., site A and B). When correlated with the water quality data that suggests a significantly lower dissolved organic carbon (DOC) concentration at site D than the other two RWDS, the metagenomic data suggest that low DOC is needed to maintain the biological stability of reclaimed water along the distribution network.

Deep metagenomics examines the oral microbiome during dental caries, revealing novel taxa and co-occurrences with host molecules

Dental caries, the most common chronic infectious disease worldwide, has a complex etiology involving the interplay of microbial and host factors that are not completely understood. In this study, the oral microbiome and 38 host cytokines and chemokines were analyzed across 23 children with caries and 24 children with healthy dentition. De novo assembly of metagenomic sequencing obtained 527 metagenome-assembled genomes (MAGs), representing 150 bacterial species. Forty-two of these species had no genomes in public repositories, thereby representing novel taxa. These new genomes greatly expanded the known pangenomes of many oral clades, including the enigmatic Saccharibacteria clades G3 and G6, which had distinct functional repertoires compared to other oral Saccharibacteria. Saccharibacteria are understood to be obligate epibionts, which are dependent on host bacteria. These data suggest that the various Saccharibacteria clades may rely on their hosts for highly distinct metabolic requirements, which would have significant evolutionary and ecological implications. Across the study group, Rothia, Neisseria, and Haemophilus spp. were associated with good dental health, whereas Prevotella spp., Streptococcus mutans, and Human herpesvirus 4 (Epstein-Barr virus [EBV]) were more prevalent in children with caries. Finally, 10 of the host immunological markers were significantly elevated in the caries group, and co-occurrence analysis provided an atlas of potential relationships between microbes and host immunological molecules. Overall, this study illustrated the oral microbiome at an unprecedented resolution and contributed several leads for further study that will increase the understanding of caries pathogenesis and guide therapeutic development.

Thermogenic hydrocarbon biodegradation by diverse depth-stratified microbial populations at a Scotian Basin cold seep

At marine cold seeps, gaseous and liquid hydrocarbons migrate from deep subsurface origins to the sediment-water interface. Cold seep sediments are known to host taxonomically diverse microorganisms, but little is known about their metabolic potential and depth distribution in relation to hydrocarbon and electron acceptor availability. Here we combined geophysical, geochemical, metagenomic and metabolomic measurements to profile microbial activities at a newly discovered cold seep in the deep sea. Metagenomic profiling revealed compositional and functional differentiation between near-surface sediments and deeper subsurface layers. In both sulfate-rich and sulfate-depleted depths, various archaeal and bacterial community members are actively oxidizing thermogenic hydrocarbons anaerobically. Depth distributions of hydrocarbon-oxidizing archaea revealed that they are not necessarily associated with sulfate reduction, which is especially surprising for anaerobic ethane and butane oxidizers. Overall, these findings link subseafloor microbiomes to various biochemical mechanisms for the anaerobic degradation of deeply-sourced thermogenic hydrocarbons.

Microbial Evolution: Chlamydial Creatures from the Deep

A metagenomic study of marine sediments from a hydrothermal vent field in the Arctic Mid-Ocean Ridge revealed wider diversity amongst members of the phylum Chlamydiae than was previously known. Unlike known chlamydiae, some of the newly described marine-sediment species may be potentially free-living.

Optimal estimation of bacterial growth rates based on a permuted monotone matrix

Motivated by the problem of estimating bacterial growth rates for genome assemblies from shotgun metagenomic data, we consider the permuted monotone matrix model $Y=\Theta\Pi+Z$ where $Y\in \mathbb{R}^{n\times p}$ is observed, $\Theta\in \mathbb{R}^{n\times p}$ is an unknown approximately rank-one signal matrix with monotone rows, $\Pi \in \mathbb{R}^{p\times p}$ is an unknown permutation matrix, and $Z\in \mathbb{R}^{n\times p}$ is the noise matrix. In this article we study estimation of the extreme values associated with the signal matrix $\Theta$, including its first and last columns and their difference. Treating these estimation problems as compound decision problems, minimax rate-optimal estimators are constructed using the spectral column-sorting method. Numerical experiments on simulated and synthetic microbiome metagenomic data are conducted, demonstrating the superiority of the proposed methods over existing alternatives. The methods are illustrated by comparing the growth rates of gut bacteria in inflammatory bowel disease patients and control subjects.

Benchmarking microbial growth rate predictions from metagenomes

Growth rates are central to understanding microbial interactions and community dynamics. Metagenomic growth estimators have been developed, specifically codon usage bias (CUB) for maximum growth rates and “peak-to-trough ratio” (PTR) for in situ rates. Both were originally tested with pure cultures, but natural populations are more heterogeneous, especially in individual cell histories pertinent to PTR. To test these methods, we compared predictors with observed growth rates of freshly collected marine prokaryotes in unamended seawater. We prefiltered and diluted samples to remove grazers and greatly reduce virus infection, so net growth approximated gross growth. We sampled over 44 h for abundances and metagenomes, generating 101 metagenome-assembled genomes (MAGs), including Actinobacteria, Verrucomicrobia, SAR406, MGII archaea, etc. We tracked each MAG population by cell-abundance-normalized read recruitment, finding growth rates of 0 to 5.99 per day, the first reported rates for several groups, and used these rates as benchmarks. PTR, calculated by three methods, rarely correlated to growth (r ~−0.26–0.08), except for rapidly growing γ-Proteobacteria (r ~0.63–0.92), while CUB correlated moderately well to observed maximum growth rates (r = 0.57). This suggests that current PTR approaches poorly predict actual growth of most marine bacterial populations, but maximum growth rates can be approximated from genomic characteristics.

Spatial Metagenomics of Three Geothermal Sites in Pisciarelli Hot Spring Focusing on the Biochemical Resources of the Microbial Consortia

Terrestrial hot springs are of great interest to the general public and to scientists alike due to their unique and extreme conditions. These have been sought out by geochemists, astrobiologists, and microbiologists around the globe who are interested in their chemical properties, which provide a strong selective pressure on local microorganisms. Drivers of microbial community composition in these springs include temperature, pH, in-situ chemistry, and biogeography. Microbes in these communities have evolved strategies to thrive in these conditions by converting hot spring chemicals and organic matter into cellular energy. Following our previous metagenomic analysis of Pisciarelli hot springs (Naples, Italy), we report here the comparative metagenomic study of three novel sites, formed in Pisciarelli as result of recent geothermal activity. This study adds comprehensive information about phylogenetic diversity within Pisciarelli hot springs by peeking into possible mechanisms of adaptation to biogeochemical cycles, and high applicative potential of the entire set of genes involved in the carbohydrate metabolism in this environment (CAZome). This site is an excellent model for the study of biodiversity on Earth and biosignature identification, and for the study of the origin and limits of life.

Large freshwater phages with the potential to augment aerobic methane oxidation

There is growing evidence that phages with unusually large genomes are common across various microbiomes, but little is known about their genetic inventories or potential ecosystem impacts. In the present study, we reconstructed large phage genomes from freshwater lakes known to contain bacteria that oxidize methane. Of manually curated genomes, 22 (18 are complete), ranging from 159 kilobase (kb) to 527 kb in length, were found to encode the pmoC gene, an enzymatically critical subunit of the particulate methane monooxygenase, the predominant methane oxidation catalyst in nature. The phage-associated PmoC sequences show high similarity to (>90%), and affiliate phylogenetically with, those of coexisting bacterial methanotrophs, including members of Methyloparacoccus, Methylocystis and Methylobacter spp. In addition, pmoC-phage abundance patterns correlate with those of the coexisting bacterial methanotrophs, supporting host-phage relationships. Future work is needed to determine whether phage-associated PmoC has similar functions to additional copies of PmoC encoded in bacterial genomes, thus contributing to growth on methane. Transcriptomics data from Lake Rotsee (Switzerland) showed that some phage-associated pmoC genes were highly expressed in situ and, of interest, that the most rapidly growing methanotroph was infected by three pmoC-phages. Thus, augmentation of bacterial methane oxidation by pmoC-phages during infection could modulate the efflux of this potent greenhouse gas into the environment.

Ancestral Absence of Electron Transport Chains in Patescibacteria and DPANN

Recent discoveries suggest that the candidate superphyla Patescibacteria and DPANN constitute a large fraction of the phylogenetic diversity of Bacteria and Archaea. Their small genomes and limited coding potential have been hypothesized to be ancestral adaptations to obligate symbiotic lifestyles. To test this hypothesis, we performed cell-cell association, genomic, and phylogenetic analyses on 4,829 individual cells of Bacteria and Archaea from 46 globally distributed surface and subsurface field samples. This confirmed the ubiquity and abundance of Patescibacteria and DPANN in subsurface environments, the small size of their genomes and cells, and the divergence of their gene content from other Bacteria and Archaea. Our analyses suggest that most Patescibacteria and DPANN in the studied subsurface environments do not form specific physical associations with other microorganisms. These data also suggest that their unusual genomic features and prevalent auxotrophies may be a result of ancestral, minimal cellular energy transduction mechanisms that lack respiration, thus relying solely on fermentation for energy conservation.

Tools for successful proliferation: diverse strategies of nutrient acquisition by a benthic cyanobacterium

Freshwater cyanobacterial blooms have increased worldwide, channeling organic carbon into these systems, and threatening animal health through the production of cyanotoxins. Both toxic and nontoxic Microcoleus proliferations usually occur when there are moderate concentrations of dissolved inorganic nitrogen, but when phosphorus is scarce. In order to understand how Microcoleus establishes thick biofilms (or mats) on riverbeds under phosphorus-limiting conditions, we collected Microcoleus-dominated biofilms over a 19-day proliferation event for proteogenomics. A single pair of nitrogen-dependent Microcoleus species were consistently present in relatively high abundance, although each followed a unique metabolic trajectory. Neither possessed anatoxin gene clusters, and only very low concentrations of anatoxins (~2 µg kg-1) were detected, likely originating from rarer Microcoleus species also present. Proteome allocations were dominated by photosynthesizing cyanobacteria and diatoms, and data indicate biomass was actively recycled by Bacteroidetes and Myxococcales. Microcoleus likely acquired nutrients throughout the proliferation event by uptake of nitrate, urea, and inorganic and organic phosphorus. Both species also harbored genes that could be used for inorganic phosphate solubilization with pyrroloquinoline quinone cofactors produced by cohabiting Proteobacteria. Results indicate that Microcoleus are equipped with diverse mechanisms for nitrogen and phosphorus acquisition, enabling them to proliferate and out-compete others in low-phosphorus waters.

Comparative Genome-Centric Analysis of Freshwater and Marine ANAMMOX Cultures Suggests Functional Redundancy in Nitrogen Removal Processes

There is a lack of understanding of the interaction between anammox bacteria and the flanking microbial communities in both freshwater (non-saline) and marine (saline) ecosystems. Here, we present a comparative genome-based exploration of two different anammox bioreactors, through the analysis of 23 metagenome-assembled genomes (MAGs), 12 from freshwater anammox reactor (FWR), and 11 from marine anammox reactor (MWR). To understand the contribution of individual members to community functions, we applied the index of replication (iRep) to determine bacteria that are actively replicating. Using genomic content and iRep information, we provided a potential ecological role for the dominant members of the community based on the reactor operating conditions. In the non-saline system, anammox (Candidatus Brocadia sinica) and auxotrophic neighboring bacteria belonging to the phyla Ignavibacteriae and Chloroflexi might interact to reduce nitrate to nitrite for direct use by anammox bacteria. Whereas, in the saline reactor, anammox bacterium (Ca. Scalindua erythraensis) and flanking community belonging to phyla Planctomycetes (different than anammox bacteria)—which persistently growing in the system—may catabolize detritus and extracellular material and recycle nitrate to nitrite for direct use by anammox bacteria. Despite different microbial communities, there was functional redundancy in both ecosystems. These results signify the potential application of marine anammox bacteria for treating saline N-rich wastewaters.

What Is Metagenomics Teaching Us, and What Is Missed?

Shotgun metagenomic sequencing has revolutionized our ability to detect and characterize the diversity and function of complex microbial communities. In this review, we highlight the benefits of using metagenomics as well as the breadth of conclusions that can be made using currently available analytical tools, such as greater resolution of species and strains across phyla and functional content, while highlighting challenges of metagenomic data analysis. Major challenges remain in annotating function, given the dearth of functional databases for environmental bacteria compared to model organisms, and the technical difficulties of metagenome assembly and phasing in heterogeneous environmental samples. In the future, improvements and innovation in technology and methodology will lead to lowered costs. Data integration using multiple technological platforms will lead to a better understanding of how to harness metagenomes. Subsequently, we will be able not only to characterize complex microbiomes but also to manipulate communities to achieve prosperous outcomes for health, agriculture, and environmental sustainability.

Microbial succession during the transition from active to inactive stages of deep-sea hydrothermal vent sulfide chimneys

BACKGROUND

Deep-sea hydrothermal vents are highly productive biodiversity hotspots in the deep ocean supported by chemosynthetic microorganisms. Prominent features of these systems are sulfide chimneys emanating high-temperature hydrothermal fluids. While several studies have investigated the microbial diversity in both active and inactive sulfide chimneys that have been extinct for up to thousands of years, little is known about chimneys that have ceased activity more recently, as well as the microbial succession occurring during the transition from active to inactive chimneys.

RESULTS

Genome-resolved metagenomics was applied to an active and a recently extinct (~ 7 years) sulfide chimney from the 9-10° N hydrothermal vent field on the East Pacific Rise. Full-length 16S rRNA gene and a total of 173 high-quality metagenome assembled genomes (MAGs) were retrieved for comparative analysis. In the active chimney (L-vent), sulfide- and/or hydrogen-oxidizing Campylobacteria and Aquificae with the potential for denitrification were identified as the dominant community members and primary producers, fixing carbon through the reductive tricarboxylic acid (rTCA) cycle. In contrast, the microbiome of the recently extinct chimney (M-vent) was largely composed of heterotrophs from various bacterial phyla, including Delta-/Beta-/Alphaproteobacteria and Bacteroidetes. Gammaproteobacteria were identified as the main primary producers, using the oxidation of metal sulfides and/or iron oxidation coupled to nitrate reduction to fix carbon through the Calvin-Benson-Bassham (CBB) cycle. Further analysis revealed a phylogenetically distinct Nitrospirae cluster that has the potential to oxidize sulfide minerals coupled to oxygen and/or nitrite reduction, as well as for sulfate reduction, and that might serve as an indicator for the early stages of chimneys after venting has ceased.

CONCLUSIONS

This study sheds light on the composition, metabolic functions, and succession of microbial communities inhabiting deep-sea hydrothermal vent sulfide chimneys. Collectively, microbial succession during the life span of a chimney could be described to proceed from a "fluid-shaped" microbial community in newly formed and actively venting chimneys supported by the oxidation of reductants in the hydrothermal fluid to a "mineral-shaped" community supported by the oxidation of minerals after hydrothermal activity has ceased. Remarkably, the transition appears to occur within the first few years, after which the communities stay stable for thousands of years. Video Abstract.

VIBRANT: automated recovery, annotation and curation of microbial viruses, and evaluation of viral community function from genomic sequences

BACKGROUND

Viruses are central to microbial community structure in all environments. The ability to generate large metagenomic assemblies of mixed microbial and viral sequences provides the opportunity to tease apart complex microbiome dynamics, but these analyses are currently limited by the tools available for analyses of viral genomes and assessing their metabolic impacts on microbiomes.

DESIGN

Here we present VIBRANT, the first method to utilize a hybrid machine learning and protein similarity approach that is not reliant on sequence features for automated recovery and annotation of viruses, determination of genome quality and completeness, and characterization of viral community function from metagenomic assemblies. VIBRANT uses neural networks of protein signatures and a newly developed v-score metric that circumvents traditional boundaries to maximize identification of lytic viral genomes and integrated proviruses, including highly diverse viruses. VIBRANT highlights viral auxiliary metabolic genes and metabolic pathways, thereby serving as a user-friendly platform for evaluating viral community function. VIBRANT was trained and validated on reference virus datasets as well as microbiome and virome data.

RESULTS

VIBRANT showed superior performance in recovering higher quality viruses and concurrently reduced the false identification of non-viral genome fragments in comparison to other virus identification programs, specifically VirSorter, VirFinder, and MARVEL. When applied to 120,834 metagenome-derived viral sequences representing several human and natural environments, VIBRANT recovered an average of 94% of the viruses, whereas VirFinder, VirSorter, and MARVEL achieved less powerful performance, averaging 48%, 87%, and 71%, respectively. Similarly, VIBRANT identified more total viral sequence and proteins when applied to real metagenomes. When compared to PHASTER, Prophage Hunter, and VirSorter for the ability to extract integrated provirus regions from host scaffolds, VIBRANT performed comparably and even identified proviruses that the other programs did not. To demonstrate applications of VIBRANT, we studied viromes associated with Crohn's disease to show that specific viral groups, namely Enterobacteriales-like viruses, as well as putative dysbiosis associated viral proteins are more abundant compared to healthy individuals, providing a possible viral link to maintenance of diseased states.

CONCLUSIONS

The ability to accurately recover viruses and explore viral impacts on microbial community metabolism will greatly advance our understanding of microbiomes, host-microbe interactions, and ecosystem dynamics. Video Abstract.

Taxon-specific microbial growth and mortality patterns reveal distinct temporal population responses to rewetting in a California grassland soil

Microbial activity increases after rewetting dry soil, resulting in a pulse of carbon mineralization and nutrient availability. The biogeochemical responses to wet-up are reasonably well understood and known to be microbially mediated. Yet, the population level dynamics, and the resulting changes in microbial community patterns, are not well understood as ecological phenomena. Here, we used sequencing of 16S rRNA genes coupled with heavy water (H218O) DNA quantitative stable isotope probing to estimate population-specific rates of growth and mortality in response to a simulated wet-up event in a California annual grassland soil. Bacterial growth and mortality responded rapidly to wet-up, within 3 h, and continued throughout the 168 h incubation, with patterns of sequential growth observed at the phylum level. Of the 37 phyla detected in the prewet community, growth was found in 18 phyla while mortality was measured in 26 phyla. Rapid growth and mortality rates were measurable within 3 h of wet-up but had contrasting characteristics; growth at 3 h was dominated by select taxa within the Proteobacteria and Firmicutes, whereas mortality was taxonomically widespread. Furthermore, across the community, mortality exhibited density-independence, consistent with the indiscriminate shock resulting from dry-down and wet-up, whereas growth was density-dependent, consistent with control by competition or predation. Total aggregated growth across the community was highly correlated with total soil CO2 production. Together, these results illustrate how previously "invisible" population responses can translate quantitatively to emergent observations of ecosystem-scale biogeochemistry.

Novel clostridial lineages recovered from metagenomes of a hot oil reservoir

Oil reservoirs have been shown to house numerous microbial lineages that differ based on the in-situ pH, salinity and temperature of the subsurface environment. Lineages of Firmicutes, including Clostridiales, have been frequently detected in oil reservoirs, but are typically not considered impactful or relevant due to their spore-forming nature. Here we show, using metagenomics, a high temperature oil reservoir of marine salinity contains a microbial population that is predominantly from within the Order Clostridiales. These organisms form an oil-reservoir specific clade based on the phylogenies of both 16S rRNA genes and ribosomal proteins, which we propose to name ^UPetromonas tenebris, meaning they are single-celled organisms from dark rocks. Metagenome-assembled genomes (MAGs) of these Petromonas sp. were obtained and used to determine that these populations, while capable of spore-formation, were also likely replicating in situ in the reservoir. We compared these MAGs to closely related genomes and show that these subsurface Clostridiales differ, from the surface derived genomes, showing signatures of the ability to degrade plant-related compounds, whereas subsurface genomes only show the ability to process simple sugars. The estimation of in-situ replication from genomic data suggest that ^UPetromonas tenebris lineages are functional in-situ and may be specifically adapted to inhabit oil reservoirs.

Counting Replication Origins to Measure Growth of Pathogens

For the past several decades, the success of bacterial strains in infecting their host has been essentially ascribed to the presence of canonical virulence genes. While it is unclear how much growth rate impacts the outcome of an infection, it is long known that the efficacy of the most commonly used antibiotics is correlated to growth. This applies especially to β-lactams, whose efficacy is nearly abolished when cells grow very slowly. It is therefore reasonable to assume that a niche or genetic dependent change in growth rate could contribute to the variability in the outcome of antibiotic therapy. However, little is known about the growth rate of pathogens or their pathotypes in their host.

Temperature and Nutrient Levels Correspond with Lineage-Specific Microdiversification in the Ubiquitous and Abundant Freshwater Genus Limnohabitans

Most freshwater bacterial communities are characterized by a few dominant taxa that are often ubiquitous across freshwater biomes worldwide. Our understanding of the genomic diversity within these taxonomic groups is limited to a subset of taxa. Here, we investigated the genomic diversity that enables Limnohabitans, a freshwater genus key in funneling carbon from primary producers to higher trophic levels, to achieve abundance and ubiquity. We reconstructed eight putative Limnohabitans metagenome-assembled genomes (MAGs) from stations located along broad environmental gradients existing in Lake Michigan, part of Earth's largest surface freshwater system. De novo strain inference analysis resolved a total of 23 strains from these MAGs, which strongly partitioned into two habitat-specific clusters with cooccurring strains from different lineages. The largest number of strains belonged to the abundant LimB lineage, for which robust in situ strain delineation had not previously been achieved. Our data show that temperature and nutrient levels may be important environmental parameters associated with microdiversification within the Limnohabitans genus. In addition, strains predominant in low- and high-phosphorus conditions had larger genomic divergence than strains abundant under different temperatures. Comparative genomics and gene expression analysis yielded evidence for the ability of LimB populations to exhibit cellular motility and chemotaxis, a phenotype not yet associated with available Limnohabitans isolates. Our findings broaden historical marker gene-based surveys of Limnohabitans microdiversification and provide in situ evidence of genome diversity and its functional implications across freshwater gradients.IMPORTANCE Limnohabitans is an important bacterial taxonomic group for cycling carbon in freshwater ecosystems worldwide. Here, we examined the genomic diversity of different Limnohabitans lineages. We focused on the LimB lineage of this genus, which is globally distributed and often abundant, and its abundance has shown to be largely invariant to environmental change. Our data show that the LimB lineage is actually comprised of multiple cooccurring populations for which the composition and genomic characteristics are associated with variations in temperature and nutrient levels. The gene expression profiles of this lineage suggest the importance of chemotaxis and motility, traits that had not yet been associated with the Limnohabitans genus, in adapting to environmental conditions.

Vermamoeba vermiformis CDC-19 draft genome sequence reveals considerable gene trafficking including with candidate phyla radiation and giant viruses

Vermamoeba vermiformis is a predominant free-living amoeba in human environments and amongst the most common amoebae that can cause severe infections in humans. It is a niche for numerous amoeba-resisting microorganisms such as bacteria and giant viruses. Differences in the susceptibility to these giant viruses have been observed. V. vermiformis and amoeba-resisting microorganisms share a sympatric lifestyle that can promote exchanges of genetic material. This work analyzed the first draft genome sequence of a V. vermiformis strain (CDC-19) through comparative genomic, transcriptomic and phylogenetic analyses. The genome of V. vermiformis is 59.5 megabase pairs in size, and 22,483 genes were predicted. A high proportion (10% (n = 2,295)) of putative genes encoded proteins showed the highest sequence homology with a bacterial sequence. The expression of these genes was demonstrated for some bacterial homologous genes. In addition, for 30 genes, we detected best BLAST hits with members of the Candidate Phyla Radiation. Moreover, 185 genes (0.8%) best matched with giant viruses, mostly those related to the subfamily Klosneuvirinae (101 genes), in particular Bodo saltans virus (69 genes). Lateral sequence transfers between V. vermiformis and amoeba-resisting microorganisms were strengthened by Sanger sequencing, transcriptomic and phylogenetic analyses. This work provides important insights and genetic data for further studies about this amoeba and its interactions with microorganisms.

Metagenomic growth rate inferences of strains in situ

We developed a method for strain-level metagenomic estimation of growth rate (SMEG) for inferring growth rates of bacterial subspecies, or strains, from complex metagenomic samples. We applied our method, which is based on both reference strains and de novo approaches, to different gut metagenomic datasets, accurately identifying an outbreak-associated Escherichia coli strain and a previously unidentified association of an Akkermansia muciniphila strain in cancer immunotherapy responders. SMEG resolves strain-specific growth rates from mixtures of commensal or pathogenic strains to provide new insights into microbial interactions and disease associations at the strain level. SMEG is available for download at https://github.com/ohlab/SMEG.

Probabilistic model based on circular statistics for quantifying coverage depth dynamics originating from DNA replication

BACKGROUND

With the development of DNA sequencing technology, static omics profiling in microbial communities, such as taxonomic and functional gene composition determination, has become possible. Additionally, the recently proposed in situ growth rate estimation method allows the applicable range of current comparative metagenomics to be extended to dynamic profiling. However, with this method, the applicable target range is presently limited. Furthermore, the characteristics of coverage depth during replication have not been sufficiently investigated.

RESULTS

We developed a probabilistic model that mimics coverage depth dynamics. This statistical model explains the bias that occurs in the coverage depth due to DNA replication and errors that arise from coverage depth observation. Although our method requires a complete genome sequence, it involves a stable to low coverage depth (>0.01×). We also evaluated the estimation using real whole-genome sequence datasets and reproduced the growth dynamics observed in previous studies. By utilizing a circular distribution in the model, our method facilitates the quantification of unmeasured coverage depth features, including peakedness, skewness, and degree of density, around the replication origin. When we applied the model to time-series culture samples, the skewness parameter, which indicates the asymmetry, was stable over time; however, the peakedness and degree of density parameters, which indicate the concentration level at the replication origin, changed dynamically. Furthermore, we demonstrated the activity measurement of multiple replication origins in a single chromosome.

CONCLUSIONS

We devised a novel framework for quantifying coverage depth dynamics. Our study is expected to serve as a basis for replication activity estimation from a broader perspective using the statistical model.

Accurate and complete genomes from metagenomes

Genomes are an integral component of the biological information about an organism; thus, the more complete the genome, the more informative it is. Historically, bacterial and archaeal genomes were reconstructed from pure (monoclonal) cultures, and the first reported sequences were manually curated to completion. However, the bottleneck imposed by the requirement for isolates precluded genomic insights for the vast majority of microbial life. Shotgun sequencing of microbial communities, referred to initially as community genomics and subsequently as genome-resolved metagenomics, can circumvent this limitation by obtaining metagenome-assembled genomes (MAGs); but gaps, local assembly errors, chimeras, and contamination by fragments from other genomes limit the value of these genomes. Here, we discuss genome curation to improve and, in some cases, achieve complete (circularized, no gaps) MAGs (CMAGs). To date, few CMAGs have been generated, although notably some are from very complex systems such as soil and sediment. Through analysis of about 7000 published complete bacterial isolate genomes, we verify the value of cumulative GC skew in combination with other metrics to establish bacterial genome sequence accuracy. The analysis of cumulative GC skew identified potential misassemblies in some reference genomes of isolated bacteria and the repeat sequences that likely gave rise to them. We discuss methods that could be implemented in bioinformatic approaches for curation to ensure that metabolic and evolutionary analyses can be based on very high-quality genomes.

Marine Sediments Illuminate Chlamydiae Diversity and Evolution

The bacterial phylum Chlamydiae is so far composed of obligate symbionts of eukaryotic hosts. Well known for Chlamydiaceae, pathogens of humans and other animals, Chlamydiae also include so-called environmental lineages that primarily infect microbial eukaryotes. Environmental surveys indicate that Chlamydiae are found in a wider range of environments than anticipated previously. However, the vast majority of this chlamydial diversity has been underexplored, biasing our current understanding of their biology, ecological importance, and evolution. Here, we report that previously undetected and active chlamydial lineages dominate microbial communities in deep anoxic marine sediments taken from the Arctic Mid-Ocean Ridge. Reaching relative abundances of up to 43% of the bacterial community, and a maximum diversity of 163 different species-level taxonomic units, these Chlamydiae represent important community members. Using genome-resolved metagenomics, we reconstructed 24 draft chlamydial genomes, expanding by over a third the known genomic diversity in this phylum. Phylogenomic analyses revealed several novel clades across the phylum, including a previously unknown sister lineage of the Chlamydiaceae, providing new insights into the origin of pathogenicity in this family. We were unable to identify putative eukaryotic hosts for these marine sediment chlamydiae, despite identifying genomic features that may be indicative of host-association. The high abundance and genomic diversity of Chlamydiae in these anoxic marine sediments indicate that some members could play an important, and thus far overlooked, ecological role in such environments and may indicate alternate lifestyle strategies.

The Tempo and Mode of Angiosperm Mitochondrial Genome Divergence Inferred from Intraspecific Variation in Arabidopsis thaliana

The mechanisms of sequence divergence in angiosperm mitochondrial genomes have long been enigmatic. In particular, it is difficult to reconcile the rapid divergence of intergenic regions that can make non-coding sequences almost unrecognizable even among close relatives with the unusually high levels of sequence conservation found in genic regions. It has been hypothesized that different mutation and repair mechanisms act on genic and intergenic sequences or alternatively that mutational input is relatively constant but that selection has strikingly different effects on these respective regions. To test these alternative possibilities, we analyzed mtDNA divergence within Arabidopsis thaliana, including variants from the 1001 Genomes Project and changes accrued in published mutation accumulation (MA) lines. We found that base-substitution frequencies are relatively similar for intergenic regions and synonymous sites in coding regions, whereas indel and nonsynonymous substitutions rates are greatly depressed in coding regions, supporting a conventional model in which mutation/repair mechanisms are consistent throughout the genome but differentially filtered by selection. Most types of sequence and structural changes were undetectable in 10-generation MA lines, but we found significant shifts in relative copy number across mtDNA regions for lines grown under stressed vs. benign conditions. We confirmed quantitative variation in copy number across the A. thaliana mitogenome using both whole-genome sequencing and droplet digital PCR, further undermining the classic but oversimplified model of a circular angiosperm mtDNA structure. Our results suggest that copy number variation is one of the most fluid features of angiosperm mitochondrial genomes.