Integration of time-series meta-omics data reveals how microbial ecosystems respond to disturbance.

The development of reliable, mixed-culture biotechnological processes hinges on understanding how microbial ecosystems respond to disturbances. Here we reveal extensive phenotypic plasticity and niche complementarity in oleaginous microbial populations from a biological wastewater treatment plant. We perform meta-omics analyses (metagenomics, metatranscriptomics, metaproteomics and metabolomics) on in situ samples over 14 months at weekly intervals. Based on 1,364 de novo metagenome-assembled genomes, we uncover four distinct fundamental niche types. Throughout the time-series, we observe a major, transient shift in community structure, coinciding with substrate availability changes. Functional omics data reveals extensive variation in gene expression and substrate usage amongst community members. Ex situ bioreactor experiments confirm that responses occur within five hours of a pulse disturbance, demonstrating rapid adaptation by specific populations. Our results show that community resistance and resilience are a function of phenotypic plasticity and niche complementarity, and set the foundation for future ecological engineering efforts

MEGAN-LR: new algorithms allow accurate binning and easy interactive exploration of metagenomic long reads and contigs

Abstract Background There are numerous computational tools for taxonomic or functional analysis of microbiome samples, optimized to run on hundreds of millions of short, high quality sequencing reads. Programs such as MEGAN allow the user to interactively navigate these large datasets. Long read sequencing technologies continue to improve and produce increasing numbers of longer reads (of varying lengths in the range of 10k-1M bps, say), but of low quality. There is an increasing interest in using long reads in microbiome sequencing, and there is a need to adapt short read tools to long read datasets. Methods We describe a new LCA-based algorithm for taxonomic binning, and an interval-tree based algorithm for functional binning, that are explicitly designed for long reads and assembled contigs. We provide a new interactive tool for investigating the alignment of long reads against reference sequences. For taxonomic and functional binning, we propose to use LAST to compare long reads against the NCBI-nr protein reference database so as to obtain frame-shift aware alignments, and then to process the results using our new methods. Results All presented methods are implemented in the open source edition of MEGAN, and we refer to this new extension as MEGAN-LR (MEGAN long read). We evaluate the LAST+MEGAN-LR approach in a simulation study, and on a number of mock community datasets consisting of Nanopore reads, PacBio reads and assembled PacBio reads. We also illustrate the practical application on a Nanopore dataset that we sequenced from an anammox bio-rector community. Reviewers This article was reviewed by Nicola Segata together with Moreno Zolfo, Pete James Lockhart and Serghei Mangul. Conclusion This work extends the applicability of the widely-used metagenomic analysis software MEGAN to long reads. Our study suggests that the presented LAST+MEGAN-LR pipeline is sufficiently fast and accurate

Annotated bacterial chromosomes from frame-shift-corrected long-read metagenomic data

Abstract Background Short-read sequencing technologies have long been the work-horse of microbiome analysis. Continuing technological advances are making the application of long-read sequencing to metagenomic samples increasingly feasible. Results We demonstrate that whole bacterial chromosomes can be obtained from an enriched community, by application of MinION sequencing to a sample from an EBPR bioreactor, producing 6 Gb of sequence that assembles into multiple closed bacterial chromosomes. We provide a simple pipeline for processing such data, which includes a new approach to correcting erroneous frame-shifts. Conclusions Advances in long-read sequencing technology and corresponding algorithms will allow the routine extraction of whole chromosomes from environmental samples, providing a more detailed picture of individual members of a microbiome

Genome analysis of Pseudomonas sp. OF001 and Rubrivivax sp. A210 suggests multicopper oxidases catalyze manganese oxidation required for cylindrospermopsin transformation

Abstract Background Cylindrospermopsin is a highly persistent cyanobacterial secondary metabolite toxic to humans and other living organisms. Strain OF001 and A210 are manganese-oxidizing bacteria (MOB) able to transform cylindrospermopsin during the oxidation of Mn2+. So far, the enzymes involved in manganese oxidation in strain OF001 and A210 are unknown. Therefore, we analyze the genomes of two cylindrospermopsin-transforming MOB, Pseudomonas sp. OF001 and Rubrivivax sp. A210, to identify enzymes that could catalyze the oxidation of Mn2+. We also investigated specific metabolic features related to pollutant degradation and explored the metabolic potential of these two MOB with respect to the role they may play in biotechnological applications and/or in the environment. Results Strain OF001 encodes two multicopper oxidases and one haem peroxidase potentially involved in Mn2+ oxidation, with a high similarity to manganese-oxidizing enzymes described for Pseudomonas putida GB-1 (80, 83 and 42% respectively). Strain A210 encodes one multicopper oxidase potentially involved in Mn2+ oxidation, with a high similarity (59%) to the manganese-oxidizing multicopper oxidase in Leptothrix discophora SS-1. Strain OF001 and A210 have genes that might confer them the ability to remove aromatic compounds via the catechol meta- and ortho-cleavage pathway, respectively. Based on the genomic content, both strains may grow over a wide range of O2 concentrations, including microaerophilic conditions, fix nitrogen, and reduce nitrate and sulfate in an assimilatory fashion. Moreover, the strain A210 encodes genes which may convey the ability to reduce nitrate in a dissimilatory manner, and fix carbon via the Calvin cycle. Both MOB encode CRISPR-Cas systems, several predicted genomic islands, and phage proteins, which likely contribute to their genome plasticity. Conclusions The genomes of Pseudomonas sp. OF001 and Rubrivivax sp. A210 encode sequences with high similarity to already described MCOs which may catalyze manganese oxidation required for cylindrospermopsin transformation. Furthermore, the analysis of the general metabolism of two MOB strains may contribute to a better understanding of the niches of cylindrospermopsin-removing MOB in natural habitats and their implementation in biotechnological applications to treat water

The Isolate Caproiciproducens sp. 7D4C2 Produces n-Caproate at Mildly Acidic Conditions From Hexoses: Genome and rBOX Comparison With Related Strains and Chain-Elongating Bacteria

10.3389/fmicb.2020.594524Frontiers in Microbiology1159452

Comparative genomics of members of the genus Defluviicoccus with insights into their ecophysiological importance

Members of the genus Defluviicoccus occur often at high abundances in activated sludge wastewater treatment plants designed to remove phosphorus, where biomass is subjected to alternating anaerobic feed/aerobic famine conditions, believed to favor the proliferation of organisms like Ca. Accumulibacter and other phosphate-accumulating organisms (PAO), and Defluviicoccus. All have a capacity to assimilate readily metabolizable substrates and store them intracellularly during the anaerobic feed stage so that under the subsequent famine aerobic stage, these can be used to synthesize polyphosphate reserves by the PAO and glycogen by Defluviicoccus. Consequently, Defluviicoccus is described as a glycogen-accumulating organism or GAO. Because they share a similar anaerobic phenotype, it has been proposed that at high Defluviicoccus abundance, the PAO are out-competed for assimilable metabolites anaerobically, and hence aerobic P removal capacity is reduced. Several Defluviicoccus whole genome sequences have been published (Ca. Defluviicoccus tetraformis, Defluviicoccus GAO-HK, and Ca. Defluviicoccus seviourii). The available genomic data of these suggest marked metabolic differences between them, some of which have ecophysiological implications. Here, we describe the whole genome sequence of the type strain Defluviicoccus vanusT , the only cultured member of this genus, and a detailed comparative re-examination of all extant Defluviicoccus genomes. Each, with one exception, which appears not to be a member of this genus, contains the genes expected of GAO members, in possessing multiple copies of those for glycogen biosynthesis and catabolism, and anaerobic polyhydroxyalkanoate (PHA) synthesis. Both 16S rRNA and genome sequence data suggest that the current recognition of four clades is insufficient to embrace their phylogenetic biodiversity, but do not support the view that they should be re-classified into families other than their existing location in the Rhodospirillaceae. As expected, considerable variations were seen in the presence and numbers of genes encoding properties associated with key substrate assimilation and metabolic pathways. Two genomes also carried the pit gene for synthesis of the low-affinity phosphate transport protein, pit, considered by many to distinguish all PAO from GAO. The data re-emphasize the risks associated with extrapolating the data generated from a single Defluviicoccus population to embrace all members of that genus.Ministry of Education (MOE)National Research Foundation (NRF)Published versionThis work was supported by the Singapore National Research Foundation and Ministry of Education under the Research Centre of Excellence Programme

Dynamics of the fecal microbiome and antimicrobial resistome in commercial piglets during the weaning period

10.1038/s41598-021-97586-9Scientific Reports1111809

Identification of the mechanisms responsible for anomalies in the tropical lower thermosphere/ionosphere caused by the January 2009 sudden stratospheric warming

We apply the Entire Atmosphere GLobal (EAGLE) model to investigate the upper atmosphere response to the January 2009 sudden stratospheric warming (SSW) event. The model successfully reproduces neutral temperature and total electron content (TEC) observations. Using both model and observational data, we identify a cooling in the tropical lower thermosphere caused by the SSW. This cooling affects the zonal electric field close to the equator, leading to an enhanced vertical plasma drift. We demonstrate that along with a SSW-related wind disturbance, which is the main source to form a dynamo electric field in the ionosphere, perturbations of the ionospheric conductivity also make a significant contribution to the formation of the electric field response to SSW. The post-sunset TEC enhancement and pre-sunrise electron content reduction are revealed as a response to the 2009 SSW. We show that at post-sunset hours the SSW affects low-latitude TEC via a disturbance of the meridional electric field. We also show that the phase change of the semidiurnal migrating solar tide (SW2) in the neutral wind caused by the 2009 SSW at the altitude of the dynamo electric field generation has a crucial importance for the SW2 phase change in the zonal electric field. Such changes lead to the appearance of anomalous diurnal variability of the equatorial electromagnetic plasma drift and subsequent low-latitudinal TEC disturbances in agreement with available observations. Plain Language Summary – Entire Atmosphere GLobal model (EAGLE) interactively calculates the troposphere, stratosphere, mesosphere, thermosphere, and plasmasphere–ionosphere system states and their response to various natural and anthropogenic forcing. In this paper, we study the upper atmosphere response to the major sudden stratospheric warming that occurred in January 2009. Our results agree well with the observed evolution of the neutral temperature in the upper atmosphere and with low-latitude ionospheric disturbances over America. For the first time, we identify an SSW-related cooling in the tropical lower thermosphere that, in turn, could provide additional information for understanding the mechanisms for the generation of electric field disturbances observed at low latitudes. We show that the SSW-related vertical electromagnetic drift due to electric field disturbances is a key mechanism for interpretation of an observed anomalous diurnal development of the equatorial ionization anomaly during the 2009 SSW event. We demonstrate that the link between thermospheric winds and the ionospheric dynamo electric field during the SSW is attained through the modulation of the semidiurnal migrating solar tide

Recovery of complete genomes and non-chromosomal replicons from activated sludge enrichment microbial communities with long read metagenome sequencing

10.1038/s41522-021-00196-6npj Biofilms and Microbiomes712

Integrated time-resolved multi-omics for understanding microbial niche ecology

Microbial communities are strongly shaped by the niche breadths of their constituent populations. However, a detailed understanding of microbial niche ecology is typically lacking. Integrated multi-omic analyses of host- or environment-derived samples offer the prospect of resolving fundamental and realised niches in situ. In turn, this is considered a prerequisite for niche engineering in order to drive an individual population or a community towards a specific phenotype, e.g., improvement of a biotechnological process. Here, we sampled floating islets on the surface of an activated sludge tank in a time-series spanning 51 time-points over 14 months. Multi-omics datasets (metagenomics, metatranscriptomics, metaproteomics, and (meta-)metabolomics) were generated for all time-points. Leveraging nucleotide sequencing data, we analyzed the community structure and reconstructed genomes for specific populations of interest. Moreover, based on their metabolic potential, three major groups emerged, serving as proxies for their respective fundamental niches . Time-resolved linkage of the proteomic and transcriptomic data to the reconstructed genomes revealed a fine-grained picture of niche realization. In particular, environmental factors (temperature, metabolites, oxygen) were significantly associated with gene expression of individual populations. Furthermore, we subjected the community to controlled oxygen conditions (stable or dynamic) in a bioreactor experiment and measured the transcriptomic response. Our results suggest short-term adaptations of populations of interest with respect to lipid metabolism, among other pathways. In conclusion, our work demonstrates how longitudinal multi-omic datasets can be integrated in order to further our understanding of microbial niche ecology within a biotechnological process with potential applications beyond waste water treatment